基于二階自適應同步擠壓S變換的時變子波提取方法

江雨濛,曹思遠,陳思遠,馬敏瑤,鄭 鐸,黃 芳,曹國明

(1.中國石油大學(北京)油氣資源與探測國家重點實驗室,北京102249;2.中國石油大學(北京)地球物理學院,北京102249;3.中國電子科技集團第二十九研究所,四川成都610036;4.中國石油大港油田公司,天津300280)

地震波在地下介質中傳播時,由于地層吸收衰減和噪聲干擾等因素的影響,子波高頻成分能量衰減更快,因此實際采集的地震數據往往是非穩態的,且地震子波的頻譜形態會隨傳播時間變化。由于地震子波是反射地震勘探的重要參數[1-2],也是地震資料處理和解釋的基礎,在很多環節中都起著關鍵作用,因此地震子波提取方法的精度直接影響后續地震資料處理及解釋的精度和準確性。傳統的地震子波提取方法大多基于時不變子波的假設條件,忽略了地震數據的非穩態特征,使得子波估計結果存在誤差,進而限制了后續反演結果的分辨率。為了滿足高精度地震勘探的需求,如何從非穩態地震數據中準確提取時變子波成為近年來研究的熱點之一[3-4]。

針對上述問題,前人提出了基于分段處理的時變子波估計方法[5-10],即首先將地震記錄劃分為若干段,并假設每一段內地震子波是穩態的,進而結合不同方法提取每一段的時不變子波。此類方法依賴于分段方法和分段長度的選擇,且每段提取的平均意義上的子波與實際子波存在一定的誤差,使得后續反演結果不精確[11]。近年來,為了消除時窗長度對提取子波的影響,提高子波估計精度,王蓉蓉等[12-13]和張漫漫等[14]相繼提出基于時頻譜模擬估計時變混合相位子波,實現逐點提取子波?！?br>