[發明專利]Q值計算與反射系數反演方法、裝置、電子設備及介質

在審

• 申請號: 202010561698.X
• 申請日: 2020-06-18
• 公開（公告）號: CN113820747A
• 公開（公告）日: 2021-12-21
• 發明（設計）人: 鄭天發;鄧智勇;李金磊;李文成;熊治富;季春輝;劉勇;陳會霞
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勘探分公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創大成知識產權代理有限公司 11614
• 代理人: 高爽

1.一種Q值計算與反射系數反演方法，其特征在于，包括：

根據地層品質因子、頻散算子與褶積模型，計算衰減地震數據；

建立多道非穩態衰減地震數據的合成公式，并引入噪聲，獲得反演方程；

在所述反演方程中輸入地震數據，進行稀疏貝葉斯學習，計算多道反射系數；

給定Q值范圍，根據所述范圍內不同Q值對應的反射系數，計算對應的Lp范數；

通過Lp范數的極值位置確定最佳Q值，通過所述最佳Q值確定對應的多道反射系數。

2.根據權利要求1所述的Q值計算與反射系數反演方法，其中，所述衰減地震數據為：

其中，F^-1是傅里葉逆變換矩陣，M為頻率采樣點數，A是衰減矩陣，W為初始地震子波的離散采樣點構成的對角矩陣。

3.根據權利要求1所述的Q值計算與反射系數反演方法，其中，所述多道非穩態衰減地震數據的合成公式為：

S＝W_aR (2)

其中，W_a為計算參數，S為多道非穩態衰減地震記錄，R為多道反射系數矩陣。

4.根據權利要求1所述的Q值計算與反射系數反演方法，其中，所述反演方程為：

d＝Gm+n (3)

其中，d、G為計算參數，m為待求解的反演參數，d＝vec(S^T)，m＝vec(R^T)，vec表示矩陣的向量化計算，表示矩陣的Kronecker運算，I_L為單位矩陣，T表示矩陣轉置計算，n為噪聲。

5.根據權利要求1所述的Q值計算與反射系數反演方法，其中，通過公式(4)計算Lp范數：

其中，R為多道反射系數矩陣，p為表征L范數類型的值。

6.一種Q值計算與反射系數反演裝置，其特征在于，包括：

衰減地震數據計算模塊，根據地層品質因子、頻散算子與褶積模型，計算衰減地震數據；

反演方程建立模塊，建立多道非穩態衰減地震數據的合成公式，并引入噪聲，獲得反演方程；

反射系數計算模塊，在所述反演方程中輸入地震數據，進行稀疏貝葉斯學習，計算多道反射系數；

范數計算模塊，給定Q值范圍，根據所述范圍內不同Q值對應的反射系數，計算對應的Lp范數；

Q值確定模塊，通過Lp范數的極值位置確定最佳Q值，通過所述最佳Q值確定對應的多道反射系數。

7.根據權利要求6所述的Q值計算與反射系數反演裝置，其中，所述衰減地震數據為：

其中，F^-1是傅里葉逆變換矩陣，M為頻率采樣點數，A是衰減矩陣，W為初始地震子波的離散采樣點構成的對角矩陣。

8.根據權利要求6所述的Q值計算與反射系數反演裝置，其中，所述多道非穩態衰減地震數據的合成公式為：

S＝W_aR (2)

其中，W_a為計算參數，S為多道非穩態衰減地震記錄，R為多道反射系數矩陣；

所述反演方程為：

d＝Gm+n (3)

其中，d、G為計算參數，m為待求解的反演參數，d＝vec(S^T)，m＝vec(R^T)，vec表示矩陣的向量化計算，表示矩陣的Kronecker運算，I_L為單位矩陣，T表示矩陣轉置計算，n為噪聲。

9.一種電子設備，其特征在于，所述電子設備包括：

存儲器，存儲有可執行指令；

處理器，所述處理器運行所述存儲器中的所述可執行指令，以實現權利要求1-5中任一項所述的Q值計算與反射系數反演方法。

10.一種計算機可讀存儲介質，其特征在于，該計算機可讀存儲介質存儲有計算機程序，該計算機程序被處理器執行時實現權利要求1-5中任一項所述的Q值計算與反射系數反演方法。