[發明專利]地震時頻反射系數反演方法及裝置

在審

• 申請號: 201910998469.1
• 申請日: 2019-10-21
• 公開（公告）號: CN112764094A
• 公開（公告）日: 2021-05-07
• 發明（設計）人: 李紅兵;潘豪杰;蔡生娟;李勇根
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 周曉飛;許曼

1.一種地震時頻反射系數反演方法，其特征在于，包括：

獲取地震資料；所述地震資料包括疊后地震數據；

基于所述地震資料提取地震子波數據；

基于廣義S變換對所述疊后地震數據進行時頻分解，得到地震記錄的時頻譜；

基于廣義S變換對所述地震子波數據進行時頻分解，得到地震子波的時頻譜；

根據所述地震記錄的時頻譜和所述地震子波的時頻譜生成地震時頻反射系數。

2.根據權利要求1所述的地震時頻反射系數反演方法，其特征在于，根據所述地震記錄的時頻譜和所述地震子波的時頻譜生成地震時頻反射系數，包括：

根據所述地震記錄的時頻譜和所述地震子波的時頻譜確定時頻反射系數反演公式；

根據所述時頻反射系數反演公式生成地震時頻反射系數。

3.根據權利要求2所述的地震時頻反射系數反演方法，其特征在于，根據所述地震記錄的時頻譜和所述地震子波的時頻譜確定時頻反射系數反演公式，包括：

確定時頻反射系數反演公式為：S(f_j)＝W(f_j)ER(f_j)；

其中，S(f_j)為地震時頻譜中對應于頻率f_j(j＝1，…，N)的單頻簡諧波分量復信號，W(f_j)為由子波時頻譜中對應于頻率f_j的單頻簡諧波分量復數信號構建的子波簡諧波分量復矩陣，E為由元素exp(-i2πf₀jΔ)構成的對角復矩陣，R(f_j)為頻率f_j的簡諧波分量反射系數序列，N為樣點個數，Δ為時間采樣率，i為虛數，即

4.根據權利要求2所述的地震時頻反射系數反演方法，其特征在于，根據所述時頻反射系數反演公式生成地震時頻反射系數，包括：

根據L1-L2混合范數反演方法對所述時頻反射系數反演公式求解，以生成地震時頻反射系數。

5.根據權利要求1所述的地震時頻反射系數反演方法，其特征在于，基于所述地震資料提取地震子波數據，包括：

獲取測井資料和井旁地震記錄數據；

基于所述地震資料提取第一地震子波數據；

根據所述測井資料和所述第一地震子波數據生成第一地震記錄數據；

若所述第一地震記錄數據與所述井旁地震記錄數據相同或差值在預設范圍內，將所述第一地震子波數據作為地震子波數據。

6.根據權利要求1-5任一項所述的地震時頻反射系數反演方法，其特征在于，還包括：

根據所述地震時頻反射系數生成第一交會圖，以根據所述第一交會圖得到地震預測結果。

7.根據權利要求1-5任一項所述的地震時頻反射系數反演方法，其特征在于，還包括：

根據所述地震時頻反射系數生成頻率道集數據；

根據所述頻率道集數據確定零頻率的反射系數和梯度值；

根據所述零頻率的反射系數和所述梯度值生成第二交會圖，以根據所述第二交會圖得到地震預測結果。

8.根據權利要求7所述的地震時頻反射系數反演方法，其特征在于，根據所述頻率道集數據確定零頻率的反射系數和梯度值，包括：

根據以下公式確定零頻率的反射系數和梯度值：R(f)＝A+Gf^1/2，

其中，R(f)為頻率道集數據中的分頻反射率，f為頻率，A為截距，G為梯度。