[發明專利]一種時頻域可控震源近地表吸收補償方法

有效

• 申請號: 201210323261.8
• 申請日: 2012-09-04
• 公開（公告）號: CN103675901A
• 公開（公告）日: 2014-03-26
• 發明（設計）人: 凌云;高軍;張汝杰;郭向宇;李少英
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京市中實友知識產權代理有限責任公司 11013
• 代理人: 劉天語;任潔

1.一種時頻域可控震源近地表吸收補償的方法，特點是通過以下步驟實現：

1）使用可控震源激發并記錄地震波數據；

2）將地震波數據在時頻域進行分解；所述的分解是在頻率域采用余弦函數構造的分頻濾波函數對地震波數據進行分頻處理；

3）對步驟2）分解后的各個頻段的分頻數據進行共炮集頻率域中值濾波，產生每一個頻段的炮集模型道；

4）采用高階e指數函數作為近地表吸收衰減的擬合函數，利用炮集模型道數據擬合求取每一炮各個頻段的統計吸收衰減函數；

5）將掃描信號在時頻域進行分解為不同的頻段，利用吸收衰減函數在相應頻段形成全頻帶的掃描信號，處理得到含近地表吸收的掃描信號，再與原始掃描信號相關，得到含近地表吸收的可控震源模擬子波；

6）將步驟5）得到的每一炮的含近地表吸收的可控震源模擬子波，應用最小平方算法求出反褶積因子，應用至對應的炮集數據；

7）對所有炮數據進行步驟1）-6）的處理，實現時頻域可控震源近地表吸收補償。

2.根據權利要求1的方法，特點是步驟2）所述的分頻濾波函數是在頻率域將0-f_N折疊頻率平均分為17～25個頻段，在每個頻段用余弦函數產生一個濾波器，相鄰每個濾波器重疊一半，構造覆蓋全部頻率的多個濾波器。

3.根據權利要求1的方法，特點是步驟3）所述的中值濾波是將某一炮點或檢波點各個地震道按采樣點排列后取中值，獲得炮集模型道。

4.根據權利要求1的方法，特點是步驟4）所述的擬合計算公式如下：

ϵ=Σt=1tN{lnMIDj=1,2,···MA[Xij(t,fk)]-αi(t,fk)}2---(2)]]>

式中:A[·]為振幅；—中值統計濾波；X_ij(t,f_k)為分頻后的炮集數據；f_k為第k個頻帶；t_n為采樣時間；i—炮號；j—道號；

其中:α_i(t,f_k)＝α₀(f_k)+α₁(f_k)t+…+α_p(f_k)t^p????（3）

式中：α_i(t,f_k)為第i炮的擬合函數；p為擬合階數，為3～5；α_p為擬合系數；f_k為第k個頻段；k=1，2，…，N，k是指步驟2）所述分頻段的序號，N是頻段的個數。

5.根據權利要求1的方法，特點是步驟5）所述的分解是將掃描信號按步驟2）的方法分解為不同的頻段。