楊兆斌,張銘記,賈正虹

(1.東方地球物理公司研究院,河北涿州 072750;2.東方地球物理公司 東部勘探事業(yè)部長慶經(jīng)理部,陜西西安 710021;3.東方地球物理公司 物探技術(shù)研究中心,河北涿州 072750)

時頻～炮檢距域面波衰減方法及應(yīng)用

楊兆斌1,張銘記2,賈正虹3

(1.東方地球物理公司研究院,河北涿州 072750;2.東方地球物理公司 東部勘探事業(yè)部長慶經(jīng)理部,陜西西安 710021;3.東方地球物理公司 物探技術(shù)研究中心,河北涿州 072750)

面波以低頻強能量為特征,數(shù)據(jù)處理中常用的面波衰減技術(shù)有:頻率濾波、F-K濾波和速度濾波等。由于面波特殊的頻散現(xiàn)象,使得上述方法的應(yīng)用效果受到一定的影響。基于面波的特征隨炮檢距變化這一特點,提出一種時頻～炮檢距域面波衰減方法:①對單炮記錄中的面波區(qū)域做Fourier變換或時頻分析,了解面波頻率隨時間和炮檢距的變化規(guī)律;②設(shè)計隨炮檢距變化的濾波器,并對地震數(shù)據(jù)進行處理;③采用減去法分析所去除的面波部份,迭代調(diào)整濾波器的參數(shù);④確定最佳的隨時間～炮檢距變化壓制面波參數(shù),以保證在壓制面波的同時,不損失反射波的低頻信息。這些方法在點接收的高密度數(shù)據(jù)處理中經(jīng)過應(yīng)用,已達到了預(yù)期的效果。

面波衰減;時頻～炮檢距域;高密度;點接收;濾波

0 前言

面波是由震源產(chǎn)生的沿地面?zhèn)鞑サ牟?其垂直分量由瑞雷波組成。面波的強度及傳播規(guī)律與激發(fā)參數(shù)有關(guān),也與近地表的環(huán)境和條件相關(guān)。對于陸上地震勘探中最常見的干擾波,在采集設(shè)計中通過優(yōu)選激發(fā)條件和多檢波器組合來衰減面波。在數(shù)據(jù)分析方面,學(xué)者們進行了相關(guān)的研究[1、2],如采用分頻方式研究面波的頻散現(xiàn)象,即相速度與頻率的關(guān)系,為利用面波進行近地表結(jié)構(gòu)研究提供了依據(jù)。在數(shù)據(jù)處理中,面波壓制技術(shù)取得了明顯的進展:Saatcilar等[3]認為面波是振幅與頻率的調(diào)制后的波場,采用線性頻率調(diào)制濾波方法壓制面波,對信號的影響較小;李衛(wèi)忠等[4]提出一種波場變換方法,利用中值濾波分離面波和有效波,在壓制面波的同時,保護了低頻反射波的能量;蔡希玲[5]提出的分頻自適應(yīng)檢測異常噪聲的技術(shù),可以用于衰減強能量的面波,其分頻技術(shù)的運用提高了壓制噪聲中信號的保真度;王在民等[6]利用小波變換進行分頻,在低頻帶內(nèi)用平均能量法識別和壓制面波,取得了良好的效果。

目前,隨著高密度地震勘探的興起,野外采用了數(shù)字檢波器和點接收的方式。由于野外不組合原始數(shù)據(jù)中的面波及其發(fā)育,因此,我們應(yīng)重新認識數(shù)據(jù)中面波的特點和分布規(guī)律,并采取有效方法衰減面波。

作者在本文中簡要分析了在自由表面附近產(chǎn)生的瑞雷面波及其特點,說明了面波的頻散特性及波形隨傳播距離變化的規(guī)律,為壓制面波的方法研究提供了理論基礎(chǔ)。同時通過實際數(shù)據(jù),展示了目前高密度采集數(shù)據(jù)中的面波特征及其影響。在此基礎(chǔ)上,作者提出時頻～炮檢距域的面波壓制方法,盡量保證在衰減面波的同時,減少對低頻反射信號的損失,這在實際數(shù)據(jù)的應(yīng)用中,取得了較好的效果。

1 瑞雷面波及其特點

在彈性分界面上形成的反射波和折射波,隨著時間的增加向整個彈性空間的介質(zhì)內(nèi)傳播,統(tǒng)稱為體波,是地震勘探研究的重點。相對于體波而言,在彈性分界面還存在著一種波動,從能量上來說它們只分布在彈性分界面附近,被稱為面波。分布在自由表面附近的面波被稱為瑞雷(Rayleihg)面波,瑞雷面波包含P波與SV波二種成份,其形成可以用波動方程做嚴格的推導(dǎo)[7]。

經(jīng)研究結(jié)果表明[8],在自由表面附近,瑞雷波傳播速度與頻率無關(guān),以面波形式傳播的振動,在傳播過程中不受頻散畸變。但是,實際介質(zhì)表面并非完全“自由”,而是在其上還存在一個非彈性的覆蓋層(例如地面上的表土層)。當考慮到此覆蓋層的質(zhì)量時,瑞雷波的性質(zhì)會有所變化。這時,瑞雷波變?yōu)轭l散波,傳播速度成為頻率的函數(shù)。疏松覆蓋層的存在是導(dǎo)致瑞雷面波頻散的原因,即常規(guī)地震勘探中觀測到的瑞雷波,總是具有很明顯的頻散特征的原因。

瑞雷面波的傳播速度,取決于地面介質(zhì)的彈性常數(shù),但通常小于橫波的傳播速度。由于介質(zhì)的彈性參數(shù)隨深度而增大,所以瑞雷波的傳播速度也隨深度的增加而增大。

由理論分析可知,瑞雷波的質(zhì)點運動軌跡在xoz平面內(nèi)為一個逆進的橢圓,是一種橢圓極化波。它的波剖面為一擺線,故又稱為地滾波。瑞雷波是一種低速、低頻波,其頻譜不包含尖銳值,占據(jù)一定的頻帶范圍。對實際介質(zhì)來說,在地表面附近,特別是地震勘探中的降速帶底面附近,介質(zhì)的彈性常數(shù)變化急劇,因而瑞雷面波的速度隨波長而改變,具有頻散性質(zhì),波形也隨傳播距離而變化。

2 高密度采集數(shù)據(jù)中的面波

在工程勘探中,面波作為有效波來研究近地表特征的變化,解決工程地質(zhì)中的問題。對于油氣勘探中觀測的地震數(shù)據(jù),面波是被作為一種主要的干擾波。通過對面波的提取、分析與研究,概括起來,面波具有以下特征:

(1)能量強、頻率低。

(2)沿地表水平方向傳播,其振幅在深度方向上按指數(shù)規(guī)律衰減,其能量分布在僅限于距地面為二倍瑞雷面波的波長范圍的薄層以內(nèi)。

(3)瑞雷面波質(zhì)點位移的軌跡是一個橢圓,其垂直位移分量比水平位移分量超前π/2相位。

(4)對于半無限空間彈性介質(zhì),瑞雷面波沒有頻散現(xiàn)象,不同頻率的瑞雷面波傳播速度相同。

(5)對于層狀介質(zhì),在自由表面的瑞雷面波具有頻散現(xiàn)象,即傳播速度與波的頻率有關(guān)。

(6)視速度有規(guī)律地隨傳播速度的變化而變化,反映在地震記錄中,同相軸并非是一簇直線,而是一簇曲線。

因此,在中深層地震勘探中,面波的影響范圍較大,必須有效地消除其影響。

在高密度勘探中,面波有可能被無假頻采樣,故面波在原始記錄中有較好的連續(xù)性。圖1(a)(見下頁)為華北地區(qū)高密度試驗線中單炮記錄,道距為5m,數(shù)字檢波器單點接收,面波被充分采樣,具有明顯的特征和很強的能量,是該區(qū)資料中最典型的干擾波。如果不作面波壓制,按照常規(guī)的方法疊加,其結(jié)果如圖1(b)(見下頁)所示,即使在疊加剖面上,面波仍有極強的能量與規(guī)律性。

圖2(a)(見下頁)為我國中部地區(qū)某煤田高密度采集的原始單炮記錄,圖2上部的曲線為炮檢距,面波的出現(xiàn)規(guī)律與炮檢距分布相關(guān)。圖2(b)為圖2(a)對應(yīng)的F-X譜,從頻譜中也能夠看到面波的特征與炮檢距的關(guān)系。面波嚴重干擾了近炮檢距的有效信息,影響AVO分析、疊前信息提取及疊前偏移等方法的應(yīng)用效果。

高密度勘探的目的,是為了獲得更高的高分辨率。為了提高分辨率,必須拓寬有效信號的頻帶寬度,向低頻拓展的范圍比向高頻端拓展同樣的范圍所帶來的益處更大。這是因為向低頻拓展更易于提高地震子波頻帶的相對寬度,從而減少子波的周期數(shù),在波阻抗反演中,更有利于薄儲層分辨和解釋。此外,強能量面波的存在,會給后續(xù)的振幅補償,地表一致性處理等帶來負面的影響〔9]。由此可見,在高密度數(shù)據(jù)處理中,采取有效方法壓制面波是一項重要的工作。

3 時頻～炮檢距域面波壓制方法及實現(xiàn)步驟

常用的消除面波的方法包括一維高通濾波、二維F-K濾波、小波變換壓制面波和波場分離等。這些方法利用了面波與有效波在某一特征上的差異來壓制面波,這在實際應(yīng)用中均有效果,但也存在著一些局限性。我們根據(jù)高密度采集數(shù)據(jù)中面波的特點,即面波的振幅、頻率、速度是隨時間和炮檢距方向變化的,提出頻率隨時間和炮檢距變化的面波壓制方法,較好地保護了低頻有效信號,也提高了壓制面波的準確度。

圖1 高密度地震記錄Fig.1 H igh density seism ic data

圖2 煤田高密度三維數(shù)據(jù)及其頻譜Fig.2 3D seism ic dataw ith high density from coalm ine

為了更細致地研究高密度采集數(shù)據(jù)中面波的分布特點,這里選出三維數(shù)據(jù)中含有面波的排列分析,并進行Fourier變換,僅顯示50 Hz以內(nèi)的中低頻成份。圖3為七個排列數(shù)據(jù)的頻譜,在近排列的小炮檢距上,面波能量強,頻帶寬。隨著炮檢距的增大,面波能量減弱,頻率降低。所以,有必要以時變、空變、頻變的方式壓制面波。

同樣,也可以采用頻率掃描的方法,分析面波的能量、頻率隨炮檢距的變化規(guī)律。運用低通與高通相結(jié)合的方法,來分析面波的規(guī)律,濾波參數(shù)和濾波因子如表1所示。針對數(shù)據(jù)中面波的不同特征,濾波參數(shù)可做適當?shù)恼{(diào)整。

通過對面波區(qū)域的頻率分析與分頻掃描,確定出不同炮檢距范圍內(nèi)的面波特征,就可以有針對性地設(shè)計濾波器,對數(shù)據(jù)中的面波進行處理。本文中的方法考慮了面波能量、頻率隨炮檢距變化的特點,在面波出現(xiàn)的區(qū)域里,對于不同的時間T、不同的炮檢距,設(shè)計出不同的頻率濾波器,以不同的濾波參數(shù)壓制面波,較好地壓制了不同頻率、不同強度的面波干擾。

圖3 地震數(shù)據(jù)的F-X譜Fig.3 F-X spectrum of seism ic data

表1 面波分析的參數(shù)和濾波因子(LP為低通濾波,HP為高通濾波)Tab.1 The filterparam eters for ground ro ll analysis(LP indicates low pass filtering and HP for high pass filtering)

圖4(見下頁)為不同炮檢距對應(yīng)的濾波因子。在炮檢距20m～1 000m的范圍內(nèi),設(shè)計了八個不同參數(shù)的濾波器,以適應(yīng)面波頻率的變化。在沒有面波的炮檢距范圍,不做任何處理,這樣就可以研究面波隨時間的變化規(guī)律,設(shè)計時變的濾波器,在沒有面波的時間段,如炮檢距較大的區(qū)域,保持原數(shù)據(jù)的特征。在面波出現(xiàn)的區(qū)域內(nèi),根據(jù)面波的頻帶范圍,應(yīng)用合適的參數(shù)衰減面波,這樣就實現(xiàn)了T-F-X域的面波壓制,較好地保持了數(shù)據(jù)中的低頻成份,最大限度地壓制面波干擾。

實現(xiàn)步驟:

(1)對單炮記錄中的面波區(qū)域做Fourier變換或時頻分析,了解面波的頻帶范圍及頻率隨時間和炮檢距的變化規(guī)律。

(2)根據(jù)面波的特征,設(shè)計時間隨炮檢距變化的濾波器,并對地震數(shù)據(jù)進行處理。

(3)采用減去法觀察所去除的面波部份,通過對去除后的低頻數(shù)據(jù)分布規(guī)律的分析,調(diào)整濾波器的參數(shù)。

(4)對面波進行迭代處理,確定最佳的T-FX域的濾波器,保證在壓制面波的同時,不損失反射波的低頻信息 。

4 應(yīng)用實例

時頻～炮檢距域面波衰減方法經(jīng)軟件實現(xiàn)后,可應(yīng)用于實際高密度數(shù)據(jù)處理中。

圖5(a)為煤田高密度三維數(shù)據(jù)近排列的原始單炮記錄;圖5(b)為壓制面波后的結(jié)果;圖5(c)為去除的面波,不含低頻反射信號。

圖6(a)(見下頁)為該三維數(shù)據(jù)中遠排列的原始單炮記錄;圖6(b)為壓制面波后的結(jié)果;圖6(c)是去除的面波。

圖4 不同炮檢距對應(yīng)的濾波因子(隨序號的增加,炮檢距增大)Fig.4 The filters in differentoffsetswhich increase as the sequence num bers increase

圖5 單炮記錄及去噪效果Fig.5 Shot data and ground ro llattenuation

由于該方法根據(jù)不同炮檢距的面波特點做了針對性處理,使遠排列、近排列的面波都得到了較好地壓制,所以對三維數(shù)據(jù)有良好的適應(yīng)性。

5 結(jié)論與認識

在高密度點接收的地震數(shù)據(jù)中,面波非常發(fā)育,頻散現(xiàn)象明顯,表現(xiàn)為視速度變化的一簇曲線,嚴重影響了地震數(shù)據(jù)的品質(zhì),所以有效地壓制面波是高密度數(shù)據(jù)中的首要工作。作者在本文采用時頻～炮檢距域壓制面波方法,可以根據(jù)數(shù)據(jù)中面波的能量、頻率隨時間和炮檢距變化的特點設(shè)計濾波器,具有良好的針對性,能夠合理地保持有效信號的低頻成份,對于高密度地震數(shù)據(jù)處理,低頻勘探等方面具有現(xiàn)實意義。

[1]N IGEL A,ANSTEY.W hatever happened to ground ro ll[J].The Leading Edge,2003,(3):40.

[2]CHOON B PARK,R ICHARD D M ILLER,JIANGHA I X IA.M u ltichannel analysis of surface waves[J].Geophysics,1999,64(3):800.

[3]SAATC ILAR R,NEZIH I CAN ITEZ N.A m ethod of ground roll elim ination[J].Geophysics,1988,57(3):894.

[4]李衛(wèi)忠,張明振,王成禮,等.壓制面波的波場分離方法[J].石油地球物理勘探,1998,35(5):679.

[5]蔡希玲.聲波和強能量干擾的分配自適應(yīng)檢測與壓制方法[J].石油地球物理勘探,1999,34(4):373.

[6]王在民,李莉,朱鵬飛,等.用平均能量法分頻識別和壓制面波[J].新疆石油地質(zhì),2005,26(6):653.

[7]單娜琳,程志平,劉云楨.工程地震勘探[M].北京:冶金工業(yè)出版社,2006.

[8]孫成禹.地震波理論與方法[M].北京:中國石油大學(xué)出版社,2007.

[9]蔡希玲,賀振華,黃德濟.異常噪聲對疊前多道處理的影響[J].石油物探,2002,41(1):84.

P 315.3+1

A

1001—1749(2011)01—0030—06

2010-07-01 改回日期:2010-09-16

楊兆斌(1962-),男,工程師,從事地震數(shù)據(jù)處理工作。