徐 陽(yáng),羅明璋,王 智,曾 磊

(1.長(zhǎng)江大學(xué)電子信息學(xué)院,湖北荊州434023;2.長(zhǎng)江大學(xué)城市建設(shè)學(xué)院,湖北荊州434023)

面波是陸上地震勘探中的主要干擾波之一,尤其在我國(guó)西部山地地區(qū),復(fù)雜的地質(zhì)構(gòu)造使得面波干擾現(xiàn)象尤為突出。面波是一種具有頻散特性的規(guī)則干擾波,在炮集上呈線性分布,其特點(diǎn)是低頻率、強(qiáng)振幅,而且衰減較慢,嚴(yán)重影響中、深層的有效反射,降低了地震資料的信噪比。在疊前地震數(shù)據(jù)處理過程中,必須有效壓制面波干擾,壓制效果直接影響到地震數(shù)據(jù)的后續(xù)處理,乃至偏移成像效果,最終影響地震資料的地質(zhì)解釋結(jié)果。

為了有效壓制面波干擾,研究人員提出了各種解決方案。李文杰等[1]提出了基于頻率衰減特性的面波壓制方法。先在地震記錄上確定面波區(qū)域,然后針對(duì)不同的炮檢距選用不同的低截頻率,在選定的面波區(qū)域里進(jìn)行高通濾波,有效地衰減了面波區(qū)域里的面波能量。但是低截濾波或高通濾波會(huì)嚴(yán)重?fù)p失深層的低頻有效反射波信息。岳龍等[2]研究了連續(xù)小波變換在面波壓制中的應(yīng)用,根據(jù)面波和有效波在小波域能量和分布區(qū)域的不同,提出了利用連續(xù)小波變換自動(dòng)壓制面波的方法;陳文超等[3]提出了基于連續(xù)小波變換的自適應(yīng)面波壓制方法,根據(jù)面波干擾的影響隨炮檢距變化等特點(diǎn),采用連續(xù)小波變換實(shí)現(xiàn)自適應(yīng)面波衰減?；谛〔ㄗ儞Q實(shí)現(xiàn)面波壓制的有效性,很大程度上取決于小波函數(shù),選取不同的小波函數(shù)會(huì)取得不同的面波壓制效果。李楊等[4]提出了基于廣義S變換(GST)的面波壓制技術(shù)。采用廣義S變換將時(shí)間域地震信號(hào)變換到時(shí)頻域,在時(shí)頻域內(nèi)確定面波干擾區(qū)域,利用時(shí)頻濾波方法壓制面波,達(dá)到分頻分時(shí)衰減面波的目的。由于面波的實(shí)際頻率難以精確標(biāo)定,為了徹底去除面波,在時(shí)頻域內(nèi)確定面波干擾區(qū)域時(shí)往往存在擴(kuò)大化的現(xiàn)象,導(dǎo)致在衰減面波的同時(shí)有效反射波也會(huì)受到損害。傳統(tǒng)壓制面波的方法還有F-K濾波切除、τ-p域變換[5-6]等,這些方法壓制面波都有一定的效果,但也都有其局限性[7]。譬如F-K濾波要求有規(guī)則的空間采樣間隔,因而復(fù)雜條件下面波去除效果不佳,而且F-K濾波是將若干道記錄作為一個(gè)數(shù)據(jù)整體進(jìn)行處理,在切除面波的同時(shí),也會(huì)將包含在面波記錄中的有效低頻信息一并去除,造成對(duì)低頻有效反射波的損害。τ-p域變換是依據(jù)有效波和干擾波視速度符號(hào)和大小的不同來壓制干擾,由于面波的速度和頻率從淺到深都會(huì)發(fā)生變化,經(jīng)過τ-p域變換后,面波數(shù)據(jù)不會(huì)是一個(gè)點(diǎn),因此很難將其完全去除。近年來出現(xiàn)了一些新的采用多種變換聯(lián)合壓制面波的方法。馬見青等[8]提出了采用S變換和TT變換聯(lián)合壓制面波的方法。包乾宗等[9]提出了采用連續(xù)小波變換和脊波變換聯(lián)合壓制面波的方法。孔慶豐[10]提出了基于Hilbert-Huang變換壓制面波的方法。畢云云等[11]提出了基于離散曲波變換字典和二維局部離散余弦變換字典組合的面波壓制方法。這些利用聯(lián)合變換壓制面波的方法取得了較好效果。

廣義S變換具有可變形態(tài)的窗函數(shù),可以獲得比S變換更好的局部頻譜以及更好的時(shí)頻聚集性。利用這一優(yōu)勢(shì),對(duì)時(shí)域內(nèi)含有面波干擾的地震記錄逐道進(jìn)行廣義S變換,將單道地震記錄的一維時(shí)間域信號(hào)映射成二維時(shí)頻圖。具有低頻特征的面波信號(hào),其能量在時(shí)頻圖上沿時(shí)間軸聚集,呈長(zhǎng)條帶狀分布,反射波能量則沿頻率軸聚集,并朝頻率軸正方向收斂。根據(jù)面波與反射波各自能量聚集特征,對(duì)時(shí)頻圖進(jìn)行分析,確定出面波信號(hào)的時(shí)頻區(qū)域,針對(duì)該區(qū)域設(shè)計(jì)時(shí)頻濾波器,將位于該區(qū)域內(nèi)的信號(hào)清零,得到去除面波之后的地震記錄,再與原始地震記錄相減,可將面波從整個(gè)地震記錄中分離出來[12-13]。由于面波與有效反射波在低頻范圍內(nèi)可能存在重疊,在用廣義S變換分離面波的過程中,會(huì)損害重疊區(qū)內(nèi)低頻有效反射波。為減少這種損害,選擇有利于面波壓制的小波函數(shù),將分離得到的面波記錄再次進(jìn)行二維離散小波分解,得到4個(gè)小波系數(shù)分量,對(duì)面波主要集中的低頻高波數(shù)分量進(jìn)行高通濾波,濾除面波成分,然后再進(jìn)行小波系數(shù)分量重構(gòu),重構(gòu)結(jié)果中將不再包含面波成分,但有效反射波信號(hào)得以保留。最后將重構(gòu)結(jié)果疊加到已實(shí)現(xiàn)面波分離的地震記錄中,實(shí)現(xiàn)在壓制面波干擾的同時(shí),最大限度地減少對(duì)低頻有效反射波的損害。本文在前人研究的基礎(chǔ)上提出了一種新的廣義S變換與二維離散小波變換聯(lián)合壓制面波的方法,并通過理論模型測(cè)試和實(shí)際地震資料處理驗(yàn)證了本文方法的有效性。

1 基本原理

1.1 廣義S變換

S變換是由STOCKWELL等[14]于1996年提出的一種加時(shí)窗傅里葉變換,是介于短時(shí)傅里葉變換和連續(xù)小波變換之間的一種非平穩(wěn)信號(hào)分析處理方法。傳統(tǒng)S變換采用形態(tài)固定的高斯窗函數(shù),給實(shí)際應(yīng)用帶來一些不便。許多學(xué)者對(duì)時(shí)窗函數(shù)進(jìn)行改造,提出了廣義S變換[15-17]。MANSINHA等[18]提出了如下可變形態(tài)窗函數(shù):

(1)

式中:t為時(shí)間;f為頻率;k為時(shí)窗尺度。當(dāng)k取值較大時(shí),時(shí)窗變寬,可提高信號(hào)的頻率分辨率,當(dāng)k取值較小時(shí),時(shí)窗變窄,可提高信號(hào)的時(shí)間分辨率,從而得到可調(diào)節(jié)時(shí)間、頻率分辨率的廣義S變換:

(2)

式中:SGST(τ,f)為時(shí)域信號(hào)h(t)的廣義S變換結(jié)果;f為頻率參數(shù);k為時(shí)窗尺度參數(shù);τ為時(shí)窗的時(shí)間中點(diǎn)。

時(shí)域信號(hào)h(t)可通過GST反變換進(jìn)行重構(gòu),反變換公式為:

(3)

式中:H(f)為h(t)的傅里葉變換。

由于GST與其傅里葉譜保持直接聯(lián)系,在進(jìn)行正、反變換的過程中,可以實(shí)現(xiàn)信號(hào)的無損轉(zhuǎn)換,這也為實(shí)際應(yīng)用帶來極大方便。GST能夠標(biāo)定信號(hào)中不同頻率成分在各個(gè)時(shí)刻的頻譜,精確表達(dá)有效信號(hào)和干擾信號(hào)在時(shí)頻域內(nèi)各自的特征,因而可以利用GST來進(jìn)行時(shí)頻濾波。一般時(shí)域信號(hào)可表示為有效信號(hào)s(t)與干擾信號(hào)n(t)之和:

(4)

上式兩邊同時(shí)取GST:

(5)

對(duì)由(5)式得到的時(shí)頻圖進(jìn)行分析,根據(jù)干擾信號(hào)的時(shí)間和頻率范圍,確定其時(shí)頻區(qū)域D,針對(duì)該區(qū)域設(shè)計(jì)時(shí)頻濾波器:

(6)

用時(shí)頻濾波器對(duì)(5)式進(jìn)行濾波處理,使得SGST(n)項(xiàng)為0,保留SGST(s)項(xiàng),經(jīng)GST反變換就可以得到去除干擾之后的有效時(shí)域信號(hào)。

1.2 二維離散小波變換

對(duì)L2(R)空間滿足允許條件的小波母函數(shù)Ψ(x)進(jìn)行伸縮或平移,可得一維連續(xù)小波函數(shù):

(7)

式中:a為與頻率相關(guān)的尺度因子;τ為與時(shí)間相關(guān)的平移因子。如果取a=2j,τ=2jk,其中j,k∈Z,則可得到一維離散小波函數(shù):

(8)

同理,可定義一維離散尺度函數(shù):

(9)

將一維小波函數(shù)推廣,可定義二維離散小波函數(shù)及二維離散尺度函數(shù)。若h(x)和g(x)分別為與Φj,k(x)和Ψj,k(x)對(duì)應(yīng)的鏡像濾波器函數(shù),可得到二維離散小波分解的Mallat算法公式[19]:

(10)

實(shí)際地震記錄是二維數(shù)據(jù),水平方向表示偏移距,垂直方向表示時(shí)間,利用Mallat公式進(jìn)行二維離散小波分解,可得到如圖1所示的4個(gè)小波系數(shù)分量:低頻低波數(shù)分量LL、低頻高波數(shù)分量LH、高頻低波數(shù)分量HL和高頻高波數(shù)分量HH,LL分量還可以繼續(xù)進(jìn)行下一層分解。二維離散小波變換對(duì)信號(hào)采取二倍抽取方式,隨著分解層數(shù)的增加,同一類型分量的頻率和波數(shù)會(huì)成倍數(shù)遞減,例如LL3的頻率和波數(shù)是LL2的一半,而LL2的頻率和波數(shù)又是LL1的一半。

圖1 二維離散小波分解示意

1.3 聯(lián)合壓制面波的步驟

利用GST和二維離散小波變換聯(lián)合壓制面波的具體步驟如下。

1) 從單炮地震記錄中按順序依次抽取每一道信號(hào)進(jìn)行廣義S變換,得到該道信號(hào)對(duì)應(yīng)的二維時(shí)頻圖。

2) 面波和有效反射波在二維時(shí)頻圖上的能量聚集和圖像顯示明顯不同,面波能量在低頻范圍內(nèi)沿時(shí)間軸聚集,其圖像呈長(zhǎng)條帶狀分布,反射波能量則沿頻率軸聚集,其圖像向頻率軸正方向收斂,據(jù)此可在時(shí)頻圖上找出面波區(qū)域。

3) 針對(duì)面波區(qū)域設(shè)計(jì)時(shí)頻濾波器,將面波對(duì)應(yīng)區(qū)域內(nèi)的數(shù)據(jù)充0,并進(jìn)行GST反變換,得到去除面波后的單道時(shí)域信號(hào)。對(duì)整個(gè)單炮記錄重復(fù)步驟1)至步驟3),即可得到整個(gè)記錄去除面波后的結(jié)果。

4) 用原始記錄與去除面波后的結(jié)果相減,得到分離出來的面波記錄。

5) 對(duì)分離出來的面波記錄,進(jìn)行二維離散小波分解,得到四個(gè)小波系數(shù)分量,對(duì)面波成分主要集中的低頻高波數(shù)分量LH進(jìn)行高通濾波,其它分量不變。如果低頻低波數(shù)分量LL中還含有面波成分,則對(duì)該分量進(jìn)行下一層次的二維離散小波分解,同時(shí)對(duì)其中的低頻高波數(shù)分量LH進(jìn)行高通濾波,重復(fù)進(jìn)行,直到低頻低波數(shù)分量LL中不再含有面波成分。

6) 從最后一層開始逐層進(jìn)行二維離散小波分量重構(gòu),重構(gòu)結(jié)果中將不再包含面波成分。將重構(gòu)結(jié)果疊加到原始記錄去除了面波之后的結(jié)果中,即完成了對(duì)面波干擾的聯(lián)合壓制。

上述過程中,在GST環(huán)節(jié)要注意面波區(qū)域的劃分,它關(guān)系到面波分離是否徹底。在二維離散小波變換環(huán)節(jié),要注意選用合適的小波函數(shù),本文通過分析和反復(fù)試驗(yàn),結(jié)合面波信號(hào)特征和小波函數(shù)的性質(zhì),從理論分析和實(shí)際資料處理兩方面考慮,優(yōu)選出了Sym5小波函數(shù),實(shí)際應(yīng)用表明,Sym5小波函數(shù)的支撐集、濾波器長(zhǎng)度、消失矩等都適合于對(duì)面波干擾的壓制[20-21]。

2 理論模型測(cè)試

圖2為采用波場(chǎng)模擬方式生成的理論模型合成記錄,單震源中間激發(fā),二維均勻各向同性彈性波介質(zhì),采用變系數(shù)吸收邊界條件。理論模型設(shè)計(jì)了2層反射層,上層為水平層,下層為傾斜層,采樣間隔為1ms,共750個(gè)采樣點(diǎn),70道。面波頻率25Hz,反射波頻率30Hz。整個(gè)記錄不含隨機(jī)噪聲。

圖3為理論模型第6道數(shù)據(jù)進(jìn)行GST后得到的時(shí)頻圖,根據(jù)面波與反射波各自在時(shí)頻圖上的能量聚集特點(diǎn),很容易區(qū)別出來。由于面波頻率與反射波頻率很接近,從圖3a可見它們?cè)跁r(shí)頻圖的低頻范圍內(nèi)內(nèi)存在重疊區(qū)域,因此在進(jìn)行時(shí)頻濾波的時(shí)候,對(duì)面波區(qū)域充0必然會(huì)影響到低頻反射波,如圖3b所示。

圖2 理論模型合成記錄

圖3 理論模型第6道的GST時(shí)頻圖a 面波區(qū)充0前; b 面波區(qū)充0后

圖4為對(duì)理論模型全部70道數(shù)據(jù)分別進(jìn)行GST,并進(jìn)行面波區(qū)域充0,將面波信號(hào)分離之后得到的結(jié)果,與圖2對(duì)比可以發(fā)現(xiàn),低頻反射波的同相軸形狀發(fā)生了改變,這正是由于面波區(qū)域充0影響了低頻反射波的結(jié)果,換句話說,低頻反射波信號(hào)在面波分離過程中受到了損害。圖5為理論模型利用GST分離出來的面波記錄,可以看到,其中除了面波信號(hào)被分離了出來之外,還有一部分低頻反射波信號(hào)也被分離了出來。

圖4 理論模型面波分離之后的結(jié)果

為了減小對(duì)低頻反射波信號(hào)的損害,將如圖5所示從理論模型中分離出來的面波記錄再次進(jìn)行二維離散小波分解,第一層分解得到4個(gè)小波系數(shù)分量如圖6所示,可見面波成分主要集中在低頻低波數(shù)分量LL1和低頻高波數(shù)分量LH1中。對(duì)于第一層LL1分量中還含有的面波,再進(jìn)行第二層分解,同理還可進(jìn)行第三、第四層分解,每一層分解都會(huì)得到LL、LH、HL、HH 4個(gè)小波系數(shù)分量。對(duì)各層的LH分量進(jìn)行高通濾波,濾除其中的面波成分,然后將最后一層的4個(gè)小波系數(shù)分量重構(gòu),作為前一層的LL分量,依次進(jìn)行,直到第一層為止。重構(gòu)結(jié)果如圖7所示,可見圖上已經(jīng)沒有了面波成分,有效反射波得以保留。將重構(gòu)結(jié)果疊加到理論模型分離面波之后的記錄(圖4)中,得到的最終結(jié)果如圖8所示,可見面波成分得到很好的壓制,同時(shí)有效反射波信號(hào)的同相軸得到較大改善,更加接近于原始理論模型合成記錄(圖2),說明本文方法在壓制面波干擾的同時(shí)起到了保護(hù)低頻反射波的作用。

圖5 理論模型被分離出來的面波記錄

圖6 理論模型分離出來的面波記錄進(jìn)行第一層二維離散小波分解的結(jié)果a LL1分量; b HL1分量; c LH1分量; d HH1分量

圖7 理論模型分離出來的面波記錄經(jīng)二維離散小波變換重構(gòu)結(jié)果

圖8 理論模型經(jīng)過GST與二維離散小波變換聯(lián)合壓制面波的結(jié)果

3 實(shí)際資料處理

圖9為西部某山地原始單炮地震記錄,共184道,每道采樣點(diǎn)數(shù)3500,采樣間隔2ms,圖上可見典型的“掃帚狀”面波發(fā)育,其下面的有效反射波信號(hào)被完全遮蓋,降低了整個(gè)記錄的信噪比。圖10a為從該地震記錄中抽取第61道信號(hào)進(jìn)行GST得到的二維時(shí)頻圖,圖10b為在時(shí)頻圖上找出面波區(qū)域,并進(jìn)行數(shù)據(jù)充0后的結(jié)果。對(duì)記錄中每一道都進(jìn)行如上相同的處理,就完成了整個(gè)記錄中的面波分離。圖11為采用GST從原始記錄中分離面波后的結(jié)果,圖12為GST分離出來的面波記錄。在利用GST進(jìn)行面波分離過程中,為了使面波徹底分離,往往會(huì)擴(kuò)大數(shù)據(jù)充0區(qū)域,從而導(dǎo)致低頻有效反射波受到損害。為了減少這種損害,對(duì)GST分離出來的面波記錄進(jìn)行二維離散小波分解,將二維地震數(shù)據(jù)轉(zhuǎn)換到時(shí)間、空間、頻率、波數(shù)四維域中,得到4個(gè)小波系數(shù)分量,對(duì)包含面波成分的低頻高波數(shù)分量進(jìn)行高通濾波,濾除面波,保留有效反射波,然后再進(jìn)行小波系數(shù)分量重構(gòu),最后將重構(gòu)結(jié)果疊加到由GST分離面波之后的記錄中,實(shí)現(xiàn)在壓制面波的同時(shí)減少對(duì)低頻有效反射波的損害。

圖9 原始單炮地震記錄

圖10 原始地震記錄第61道GST變換時(shí)頻圖a 面波區(qū)充0前; b 面波區(qū)充0后

地震記錄中面波信號(hào)在時(shí)間-空間域的分布范圍與有效波不同,主要集中在低頻低波數(shù)分量LL和低頻高波數(shù)分量LH中。對(duì)由GST分離出來的面波記錄(圖12)進(jìn)行二維離散小波分解,得到第一層分量如圖13所示,可見LL1和LH1分量中存在面波成分。對(duì)包含面波成分的LL1分量還可以進(jìn)行第二層分解。將各層分解得到的LH分量進(jìn)行高通濾波,濾除面波成分,然后從最后一層開始進(jìn)行小波系數(shù)分量重構(gòu),并作為前一層的LL分量,依次進(jìn)行,直到第一層為止。重構(gòu)結(jié)果如圖14所示,可見明顯同相軸,它們正是在GST分離面波過程中丟失的低頻反射波。其中一些同相軸的“蚯蚓化”是由于多層小波變換后的邊界效應(yīng)和累積誤差所導(dǎo)致,選擇適當(dāng)?shù)男〔ê瘮?shù)及分解層數(shù)可以降低累積誤差,重構(gòu)結(jié)果與由GST分離面波之后的結(jié)果記錄疊加后“蚯蚓化”現(xiàn)象可以得到改善,但依然存在,這是本文方法聯(lián)合壓制面波的局限性。

圖11 原始記錄GST分離面波后的結(jié)果

圖12 GST分離出來的面波記錄

圖13 對(duì)分離出來的面波記錄進(jìn)行第一層二維離散小波分解的結(jié)果a LL1分量; b HL1分量; c LH1分量; d HH1分量

將重構(gòu)結(jié)果疊加到分離面波之后的記錄(圖11)中,實(shí)現(xiàn)聯(lián)合壓制面波的最終結(jié)果如圖15所示,與圖11 對(duì)比可以看到,圖15中有效反射波的能量得到加強(qiáng),說明本文方法在壓制面波干擾的同時(shí)能夠保護(hù)低頻有效反射波。

作為對(duì)比,采用傳統(tǒng)F-K濾波法對(duì)圖9所示原始地震記錄進(jìn)行面波壓制。圖16為原始單炮地震記錄的F-K譜,可見面波與有效反射波在低頻范圍內(nèi)存在重疊區(qū)。圖中淡藍(lán)色折線以上部分是面波干擾的壓制范圍,只要將該范圍內(nèi)的數(shù)據(jù)充0切除,就可以實(shí)現(xiàn)面波壓制,但是重疊區(qū)內(nèi)面波干擾與有效信號(hào)無法完全分開,在進(jìn)行F-K濾波切除時(shí)必然會(huì)損失掉部分有效反射波。圖17為原始地震記錄進(jìn)行F-K濾波后的結(jié)果,由于重疊區(qū)的存在,切除面波的同時(shí)也切除了部分有效反射波,導(dǎo)致原來被面波覆蓋的區(qū)域內(nèi)同相軸不夠清晰。

圖14 分離面波記錄的二維離散小波重構(gòu)結(jié)果

圖15 GST與二維離散小波變換聯(lián)合壓制面波結(jié)果

為了進(jìn)一步說明本文方法的優(yōu)越性,分別對(duì)傳統(tǒng)F-K濾波、單純GST面波分離以及本文提出的GST與二維離散小波變換聯(lián)合壓制面波的結(jié)果與原始記錄一起作頻譜分析,如圖18所示。

圖18a為采用傳統(tǒng)F-K濾波壓制面波之后的頻譜,由于面波干擾和有效信號(hào)在低頻區(qū)域內(nèi)不可能完全分開,在進(jìn)行面波壓制時(shí),不可避免會(huì)損害重疊部分的有效反射波。

圖18b為采用GST進(jìn)行面波分離后的頻譜,觀察圖中藍(lán)色譜線可見,頻率在10Hz以下的面波被壓制掉了,但是,該范圍內(nèi)有效反射波也一起被壓制掉了,因此單純采用GST分離面波的方法來壓制面波干擾將造成低頻區(qū)內(nèi)有效反射波的損失。

圖16 原始單炮地震記錄的F-K譜

圖17 F-K濾波壓制面波的結(jié)果

圖18c為采用GST+二維離散小波變換聯(lián)合壓制面波后的頻譜,觀察圖中藍(lán)色譜線可見,頻率在10Hz以下的強(qiáng)振幅面波信號(hào)被壓制掉了,但該頻率范圍內(nèi)低振幅信號(hào)的能量得以保留,說明本方法在壓制面波的同時(shí),可以減少對(duì)低頻區(qū)內(nèi)有效反射波的損害。

圖18 采用不同方法壓制面波的頻譜對(duì)比a 傳統(tǒng)F-K濾波; b GST面波分離; c GST+二維離散小波變換聯(lián)合壓制面波

4 結(jié)論

利用廣義S變換將地震記錄信號(hào)從時(shí)間域轉(zhuǎn)換到時(shí)頻域,采用時(shí)頻濾波可以方便地實(shí)現(xiàn)面波分離,但分離過程可能損傷低頻有效反射波。為了避免或減少這種損傷,將GST分離出來的面波記錄進(jìn)行二維離散小波分解,對(duì)面波成分主要集中的低頻高波數(shù)分量進(jìn)行高通濾波,濾除面波,保留有效反射波,然后再進(jìn)行小波系數(shù)分量重構(gòu),最后將重構(gòu)結(jié)果疊加到已分離面波的地震記錄中,實(shí)現(xiàn)在壓制面波干擾的同時(shí),最大限度地減少對(duì)低頻有效反射波的損害,整個(gè)操作僅對(duì)地震記錄中的面波區(qū)域及面波成分進(jìn)行分離與壓制,其它域不受影響。通過理論模型測(cè)試及實(shí)際地震資料處理驗(yàn)證,本文提出的這種聯(lián)合壓制面波方法切實(shí)可行,面波壓制效果明顯,在壓制過程中能較好地保護(hù)低頻有效反射波,提高了地震記錄的信噪比。

參 考 文 獻(xiàn)

[1] 李文杰,魏修成,寧俊瑞,等.基于頻率衰減特性的面波壓制方法[J].石油物探,2008,47(3):225-227

LI W J,WEI X C,NING J R,et al.A method for suppressing surface wave by considering the characteristic of frequency attenuation[J].Geophysical Prospecting for Petroleum,2008,47(3):225-227

[2] 岳龍,劉懷山,尹燕欣,等.基于連續(xù)小波變換的面波衰減方法研究[J].石油物探,2016,55(2):214-222

YUE L,LIU H S,YIN Y X,et al.Attenuation of ground roll based on continuous wavelet transform[J].Geophysical Prospecting for Petroleum,2016,55(2):214-222

[3] 陳文超,高靜懷,包乾宗.基于連續(xù)小波變換的自適應(yīng)面波壓制方法[J].地球物理學(xué)報(bào),2009,52(11):2854-2861

CHEN W C,GAO J H,BAO Q Z.Adaptive attenuation of ground roll via continuous wavelet transform[J].Chinese Journal of Geophysics,2009,52(11):2854-2861

[4] 李楊,董守華,杜慶順.基于廣義S變換的面波壓制技術(shù)研究[J].能源技術(shù)與管理,2009(3):116-118

LI Y,DONG S H,DU Q S.Surface wave suppression technique based on generalized S-transform analysis[J].Energy Technology and Management,2009(3):116-118

[5] TURNET G.Aliasing in the tau-p transform and removal of spatially aliased coherent noise[J].Geophysics,1990,55(11):1496-1503

[6] SPITZER R,NITSCHE F O,GREEN A G.Reducing source-generated noise in shallow seismic data using linear and hyperbolic transformations[J].Geophysics,2001,66(5):1612-1621

[7] 胡國(guó)斌,王智杰.幾種面波壓制方法的討論[J].石化技術(shù),2016,23(3):65-66

HU G B,WANG Z J.Discussion on several surface wave suppression methods[J].Petrochemical Industry Technology,2016,23(3):65-66

[8] 馬見青,李慶春.S變換和TT變換聯(lián)合壓制地震面波[J].物探化探計(jì)算技術(shù),2011,33(3):252-257

MA J Q,LI Q C.Joint S transformation and TT-transformation method of surface wave suppression[J].Computing techniques for geophysical and geochemical exploration,2011,33(3):252-257

[9] 包乾宗,李慶春,陳文超,等.聯(lián)合連續(xù)小波變換和脊波變換的面波衰減方法及應(yīng)用[J].石油物探,2011,50(4):367-372

BAO Q Z,LI Q C,CHEN W C,et al.Continuous wavelet transform and ridgelet joint ground-roll attenuation method and its application[J].Geophysical Prospecting for Petroleum,2011,50(4):367-372

[10] 孔慶豐.基于Hilbert-Huang變換的面波壓制方法研究[J].石油物探,2012,51(5):446-450

KONG Q F.Surface wave suppression method based on Hilbert-Huang transform[J].Geophysical Prospecting for Petroleum,2012,51(5):446-450

[11] 畢云云,汪金菊,徐小紅,等.基于離散曲波變換字典和二維局部離散余弦變換字典組合的面波壓制[J].石油物探,2017,56(2):222-231

BI Y Y,WANG J J,XU X H,et al.Ground roll attenuation based on the combination of discrete curvelet transform dictionary and two-dimensional local discrete coine transform dictionary[J].Geophysical Prospecting for Petroleum,2017,56(2):222-231

[12] 劉霞,徐濤,段玉波,等.基于廣義S變換的時(shí)頻濾波技術(shù)研究[J].自動(dòng)化技術(shù)與應(yīng)用,2012,31(2):15-19

LIU X,XU T,DUAN Y B,et al.Study on time-frequency filtering technology based on generalized S-transform[J].Techniques of Automation and Applications,2012,31(2):15-19

[13] 隆文韜,張猛.基于廣義S變換的面波壓制方法[J].油氣地球物理,2015,13(1):23-26

LONG W T,ZHNG M.The surface wave suppression based on generalized S-transform[J].Petroleum Geophysics,2015,13(1):23-26

[14] STOCKWELL R G,MANSINHA L,LOWE R P.Localization of the complex spectrum:the S transform[J].IEEE Transtions on Signal Process,1996,44(4):998-1001

[15] 高靜懷,陳文超,李幼銘,等.廣義S變換與薄互層地震響應(yīng)分析[J].地球物理學(xué)報(bào),2003,46(4):526-532

GAO J H,CHEN W C,LI Y M,et al.Generalized S transform and seismic response analysis of the interbeds[J].Chinese Journal of Geophysics,2003,46(4):526-532

[16] 周竹生,陳友良.含可變因子的廣義S變換及其時(shí)頻濾波[J].煤田地質(zhì)與勘探,2011,39(6):63-66

ZHOU Z S,CHEN Y L.Generalized S-transform with variable-factor and its time-frequency filtering[J].Coal Geology & Exploration,2011,39(6):63-66

[17] 王德營(yíng),李振春,王姣.S域時(shí)頻濾波分析與改進(jìn)[J].石油物探,2015,54(4):396-403

WANG D Y,LI Z C,WANG J.The analysis and improvement on time-frequency filtering in S-transform domain[J].Geophysical Prospecting for Petroleum,2015,54(4):396-403

[18] MANSINHA L,STOCKWELL R G,LOWE R P.Local S-spectrum analysis of 1-D and 2-D data[J].Physics of the Earth and Planetary Interiors,1997,103(3):329-336

[19] MALLAT S.A wavelet tour of signal processing(Second Edition)[M].Beijing:China Machine Press,2003:303-312

[20] 張華,陳小宏,楊海燕.地震信號(hào)去噪的最優(yōu)小波基選取方法[J].石油地球物理勘探,2011,46(1):70-75

ZHANG H,CHEN X H,YANG H Y.Optimistic waveletbasis selection in seismic signal noise elimination[J].Oil Geophysical Prospecting,2011,46(1):70-75

[21] 楊學(xué)亭,劉財(cái),劉洋,等.基于連續(xù)小波變換的時(shí)頻域地震波能量衰減補(bǔ)償[J].石油物探,2014,53(5):523-529

YANG X T,LIU C,LIU Y,et al.The attenuation compensation of seismic wave energy in time-frequency domain based on continuous wavelet transform[J].Geophysical Prospecting for Petroleum,2014,53(5):523-529