[發明專利]深度域地震子波的確定方法和裝置

有效

• 申請號: 201710301335.0
• 申請日: 2017-05-02
• 公開（公告）號: CN107229075B
• 公開（公告）日: 2019-06-11
• 發明（設計）人: 葉月明;張金陵;徐美茹;莊錫進
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/48
• 分類號: G01V1/48
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 李輝;徐煥

1.一種深度域地震子波的確定方法，其特征在于，包括：

獲取待測區域的聲波測井數據；

根據所述聲波測井數據，確定深度域反射系數序列和井位置的深度域成像數據；

根據所述深度域反射系數序列和所述井位置的深度域成像數據，確定深度域地震子波；

其中，根據所述深度域反射系數序列和所述井位置的深度域成像數據，確定深度域地震子波，包括：

根據所述深度域反射系數序列，確定深度域反射系數序列的自相關參數；

根據所述深度域反射系數序列和所述井位置的深度域成像數據，確定深度域反射系數序列和深度域成像數據的互相關參數；

根據所述自相關參數和所述互相關參數，確定匹配算子，并將所述匹配算子作為所述深度域地震子波，其中，深度域成像數據為地震成像數據。

2.根據權利要求1所述的方法，其特征在于，根據聲波測井數據，確定深度域反射系數序列，包括：

根據所述聲波測井數據，獲取速度場數據和密度場數據；

根據所述速度場數據和所述密度場數據，確定所述深度域反射系數序列。

3.根據權利要求2所述的方法，其特征在于，根據所述速度場數據和所述密度場數據，確定所述深度域反射系數序列，包括：

按照以下公式確定所述深度域反射系數序列：

其中，R^d(z_i)為深度域反射系數序列中的第i個采樣點的深度域反射系數，υ(z_i)為速度場數據中的第i個采樣點的速度場數據，ρ(z_i)為密度場數據中的第i個采樣點的密度場數據，z_i為測井數據中的第i個采樣點的深度值，i為測井數據中深度方向上的采樣點的編號。

4.根據權利要求1所述的方法，其特征在于，根據所述深度域反射系數序列，確定深度域反射系數序列的自相關參數，包括：按照以下公式確定所述自相關參數：

其中，為第i個采樣點的自相關參數，R^d(n)為深度為n時的深度域反射系數，R^d(n+z_i)為深度為第i個采樣點的深度和n之和時的深度域反射系數，n為深度變量，取值為1至l范圍內的自然數，z_i為測井數據中的第i個采樣點的深度值，l為算子的長度。

5.根據權利要求1所述的方法，其特征在于，根據所述深度域反射系數序列和所述井位置的深度域成像數據，確定深度域反射系數序列和深度域成像數據的互相關參數，包括：按照以下公式確定所述互相關參數：

其中，μ_Rx(z_i)為第i個采樣點的互相關參數，x^d(n)為深度為n時的成像數據，R^d(n+z)為深度為第i個采樣點的深度和n之和時的深度域反射系數，n為深度變量，取值為1至l范圍內的自然數，z_i為測井數據中的第i個采樣點的深度值，l為算子的長度。

6.根據權利要求1所述的方法，其特征在于，根據所述自相關參數和所述互相關參數，確定匹配算子，包括：

根據所述自相關參數和所述互相關參數，構建求解矩陣方程；

通過求解所述求解矩陣方程，確定所述匹配算子。

7.根據權利要求6所述的方法，其特征在于，根據所述自相關參數和所述互相關參數，構建求解矩陣方程，包括：

根據所述自相關參數和所述互相關參數，構建以下方程作為所述求解矩陣方程：

其中，為深度為z時的自相關參數，μ_Rx(z)為深度為z時的互相關參數，w^d(z)為深度為z時的匹配算子，l為匹配算子的長度，z為深度，z的取值為0至l范圍內的自然數。