[發明專利]一種基于自適應增益限反Q濾波的地震資料處理方法

有效

• 申請號: 201910988175.0
• 申請日: 2019-10-17
• 公開（公告）號: CN110579805B
• 公開（公告）日: 2021-03-12
• 發明（設計）人: 李勇;張固瀾;段景;何承杰;李彪;杜皓;羅帆;詹熠宗;魏憲濱;王佳;徐峰
• 申請（專利權）人: 西南石油大學
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 成都方圓聿聯專利代理事務所(普通合伙) 51241
• 代理人: 李鵬

1.一種基于自適應增益限反Q濾波的地震資料處理方法，其特征在于：包括以下步驟：

S1.輸入地震資料包含的所有地震數據；

S2.選擇一道地震數據作為輸入信號，通過時頻分析求取輸入信號不同時間點有效頻帶的截止頻率；

S3.通過不同時間有效頻帶的截止頻率算出對應的自適應增益限，得到的振幅補償函數；

S4.若輸入信號信噪比較高，選用頻率域的實現方式實現吸收衰減補償，若輸入信號的信噪比較低，選用時頻域的實現方式實現吸收衰減補償；

S5.對吸收衰減補償后的結果沿著頻率方向求和，得到補償后的時間域信號；

S6.重復步驟S2～S5，依次對地震資料中包含的每一道地震數據進行補償。

2.根據權利要求1所述的一種基于自適應增益限反Q濾波的地震資料處理方法，其特征在于：所述步驟S2包括以下子步驟：

S201.選擇一道地震數據作為輸入信號，進行時頻分析，得到輸入信號的時頻譜S(t,f)；

S202.計算t時刻時頻譜S(t,f)的模G_t(f)，并求出t時刻的時頻譜模G_t(f)的最大值G_t(f_p)：

其中，Max[]表示取最大值；|·|表示取模運算；

S203.計算地震資料有效頻帶的截止頻率f_d(t)：

其中，d是G_t(f)的最大動態范圍，單位分貝；L_t(f)表示所有滿足的頻率成分。

3.根據權利要求1所述的一種基于自適應增益限反Q濾波的地震資料處理方法，其特征在于：所述步驟S3包括以下子步驟：

S301.計算自適應增益限C(t)和穩定因子δ₃：

式中，p，a，b是用戶預設的參數，p為實數且p＞1；a、b為實數且b≥a≥0；Q為地層品質因子，將步驟S2得到的f_d(t)作為頻率f帶入中進行計算，得到γ[t,f_d(t)]，由于f_d(t)為有效頻帶的截止頻率，即希望進行吸收衰減補償的有效頻帶的最大頻率，故令γ[t,f_d(t)]表征任意時間t的振幅補償函數的最大值，即自適應增益限C(t)；

S302.獲取得到的振幅補償函數：

其中，B₃(t,f)為自適應增益限反Q濾波振幅補償函數。

4.根據權利要求1所述的一種基于自適應增益限反Q濾波的地震資料處理方法，其特征在于：

所述步驟S4中，判斷輸入信號信噪比高低的方式如下：

S201.測量輸入信號中包含的信號能量P_signal和噪音能量P_noise；

S202.計算輸入信號的信噪比SNR：

S203.判斷信噪比SNR是否大于預設閾值；

若是，則輸入信號的信噪比較高；

若否，則輸入信號的信噪比較低。

5.根據權利要求1所述的一種基于自適應增益限反Q濾波的地震資料處理方法，其特征在于：所述步驟S4中，選用頻率域的實現方法進行吸收衰減補償的公式如下：

A(0,f)為輸入信號的頻譜，h₁(t)為頻率域處理后的時域信號，Re[]表示取復信號的實部；

選用時頻域的實現方法行吸收衰減補償的公式如下：

式中，N_t是Re[H(t,f)]≠0的總數，h₂(t)為時頻域處理后的時域信號。