貴州煤田地質(zhì)局地質(zhì)勘察研究院 黃 鑫
貴州交通科學(xué)研究院有限責(zé)任公司 梁 棟 夏曉勇
復(fù)雜地區(qū)地震資料疊前去噪技術(shù)
貴州煤田地質(zhì)局地質(zhì)勘察研究院 黃 鑫
貴州交通科學(xué)研究院有限責(zé)任公司 梁 棟 夏曉勇
在對復(fù)雜地區(qū)地震資料進(jìn)行處理時(shí),噪聲的出現(xiàn)會降低速度分析的精度,影響反褶積效果和靜校正的精度,破壞有效波的連續(xù)性,嚴(yán)重影響地震剖面的質(zhì)量。這些地區(qū)的地震資料信噪比較低、干擾嚴(yán)重且干擾波種類繁多,但高分辨率處理是以高信噪比為前提的,因此,為了獲得高信噪比的疊前地震數(shù)據(jù),需要在保真的前提下最大限度地壓制干擾波,突出有效波。
地震噪聲可能是相干噪聲,也可能是非相干噪聲。相干噪聲至少可以在一些地震道上追蹤出來,而非相干噪聲在所有的地震道上是不會有相似之處的,因此根據(jù)相鄰地震道的特征不能預(yù)測出某一道究竟有什么。相干噪聲與非相干噪聲之間的區(qū)別通常與地震記錄有關(guān),如果檢波器之間的距離變小,非相干噪聲可能看起來像相干噪聲,不過在地震記錄中討論非相干噪聲時(shí),并不考慮檢波器之間的距離。
相干噪聲可以分為2類:一類是能量基本上是沿水平方向傳播的相干噪聲,另一類是能量基本上以垂直方向到達(dá)測線的相干噪聲。這兩種噪聲之間的區(qū)別在于前者具有可重復(fù)性,而后者卻沒有。也就是說,其區(qū)別在于當(dāng)震源重復(fù)放炮時(shí),在同一地震道的同一時(shí)刻觀測到的是否為同一種噪聲。非相干噪聲通常是指隨機(jī)噪聲(空間的隨機(jī)性),隨機(jī)噪聲并不只是具有不可預(yù)測性,還具有確定的統(tǒng)計(jì)性,而且在大多數(shù)情況下,地震噪聲并不是完全隨機(jī)的。除了多次波外,根據(jù)每一道前面的地震記錄一般不會知道后面的反射波是什么樣的,因此通常認(rèn)為在時(shí)間軸上地震噪聲是趨向于隨機(jī)的。這些噪聲的相干性、傳播的方向性與可重復(fù)性,構(gòu)成了改善地震記錄質(zhì)量的大多數(shù)方法的基礎(chǔ)。現(xiàn)借用Tadeusz等所用的一個(gè)簡單的表達(dá)式來描述地震記錄。
式中,D表示地震記錄,S表示記錄中的有效成分,Nc與Nr分別表示規(guī)則干擾和不規(guī)則干擾。那么,地震資料的信噪比S/N可以表示為:
也可以寫作:
式中,N = Nc+ Nr。式(2)表明,提高信噪比應(yīng)減小分母項(xiàng),即通常說的壓制干擾,式(3)表明,提高信噪比還可以從加強(qiáng)有效波的能量來實(shí)現(xiàn),即有效波識別加強(qiáng),疊前去噪處理技術(shù)就從這兩方面著手來設(shè)計(jì)的。
1. 相干噪聲壓制。相干噪聲的能量通常相對較強(qiáng),它的存在直接影響反褶積效果,降低速度分析精度,并最終降低疊加剖面的信噪比。相干噪聲的共同特點(diǎn)是它的空間分布具有規(guī)律性,多數(shù)情況下比研究人員需要的有效信號還更有規(guī)律。相干噪聲的存在嚴(yán)重破壞了原始記錄的面貌。它的視速度較低,但頻率與有效反射波的差異不是很明顯,因而不能用帶通濾波的方法去壓制它。
目前對于相干噪聲的消除,主要方法有F-X域相干噪聲壓制、F-K域?yàn)V波、多道傾角濾波和τ-P域去噪等。這些方法主要利用相干噪聲與有效信號的視速度差異,有的也兼顧頻率的差異。
(1)F-X(頻率-空間)域相干噪聲壓制(FXCNS)。F-X域相干噪聲壓制用于壓制炮點(diǎn)產(chǎn)生的規(guī)則干擾,常見的聲波、面波都是它壓制的對象。F-X域相干噪聲壓制的基本思想是在F-X域用最小平方法估算出指定速度和頻率范圍內(nèi)的相干噪聲,然后從原來的記錄中將其減去。由于該方法是在F-X域?qū)崿F(xiàn),對空間采樣的規(guī)則性沒有特別的要求,它是先估算相干噪聲然后在將其減去,所以對有效信號的損害較小,利于保持記錄的波形特征。
FXCNS的主要參數(shù)是速度和頻率范圍,它們要確定壓制哪些相干噪聲。視速度范圍可從記錄上估量,頻率范圍可通過頻譜分析確定。
(2)F-K域?yàn)V波。F-K域?yàn)V波也叫二維濾波。它的出現(xiàn)旨在克服由于僅在頻域切除而對有效波的損失過多的缺點(diǎn),而另外增加一個(gè)波數(shù)域(K域)加以限制,以減小切除過程中有效波的損失。對于一炮地震記錄而言,當(dāng)按道沿T方向做一維傅氏變換時(shí)得到各道的頻率譜,當(dāng)按時(shí)間沿X方向做一維傅氏變換時(shí)得到各時(shí)刻的波數(shù)譜。頻率是波在單位時(shí)間內(nèi)振動的次數(shù),與此相同,波數(shù)就是波在單位距離內(nèi)振動的次數(shù)。由此可見,在波長不變的情況下,波的速度越小,波數(shù)也越小。所以,當(dāng)將一炮T-X域的地震數(shù)據(jù)通過二維傅氏變換轉(zhuǎn)換到F-K域后,有效波和相干噪聲由于視速度的差別,它們的能量將分布在不同的扇形區(qū)域,因而很容易分離開來。F-K域?yàn)V波的實(shí)現(xiàn)就是在F-K域?qū)⑾喔稍肼暷芰糠植紖^(qū)域的數(shù)據(jù)切除后再反變換回T-X域(時(shí)-空域),從而達(dá)到壓制相干噪聲的目的。
(3)多道傾角濾波。多道傾角濾波的濾波器在F-K域設(shè)計(jì)。與F-K域?yàn)V波的原理相同,相干噪聲與有效波由于視速度的差異在F-K域?qū)⒎植荚诓煌膮^(qū)域,設(shè)計(jì)濾波器也是在F-K平面上指定要去除或保留的區(qū)域,要去除的區(qū)域置為0,要保留的區(qū)域置為1,交界處附近有一斜坡漸變帶。這個(gè)區(qū)域可以是扇形也可以是多邊形,它是多道傾角濾波器的F-K域響應(yīng)。將其反變換到T-X域即得T-X域的濾波函數(shù)。
多道傾角濾波的實(shí)現(xiàn)則既可以在F-K域也可在F-X域或T-X域進(jìn)行。如果在F-K域進(jìn)行,只需將輸入記錄變換到F-K域后再乘以濾波器的F-K域響應(yīng)即可變換到T-X域即可。 如果在F-X域進(jìn)行,則要將T-X域的濾波函數(shù)和輸入記錄都變換到F-X域,二者相乘后再變換到T-X域。如果在T-X域進(jìn)行,則簡單地將T-X域的濾波函數(shù)與輸入記錄褶積即可。這三種方式?jīng)]有本質(zhì)的區(qū)別,由于計(jì)算方法不同,處理效果略有差異,花費(fèi)時(shí)間也不一樣。F-K域費(fèi)時(shí)較少,但它的假頻效應(yīng)也較明顯。
(4)τ-P域去噪。τ-P域去噪技術(shù)是又一個(gè)去除相干噪聲的手段,它所用的主要數(shù)學(xué)工具是拉冬變換。拉冬變換將T-X域的數(shù)據(jù)轉(zhuǎn)換到τ-P域。在τ-P域中,一個(gè)視速度的同相軸能量集中在一個(gè)小的局部范圍,不同視速度的同相軸能量集中在不同的位置。τ-P域去噪的做法是將那些不需要的同相軸能量從τ-P域中分離出來,并用反拉冬變換將其反變換到T-X域,然后從原來的輸入記錄將其減去。
2. 隨機(jī)噪聲衰減。在地震記錄中,隨機(jī)噪聲的存在是不可避免的。人們甚至不接受那些完全沒有隨機(jī)噪聲的記錄,認(rèn)為它不真實(shí)。在資料處理中,完全沒有隨機(jī)噪聲的記錄相當(dāng)于理論合成記錄,由于某些處理手段和計(jì)算方法的近似性,理論合成記錄的處理結(jié)果將那些近似誤差顯露無遺,所以適量隨機(jī)噪聲的存在并沒有壞處。但過于嚴(yán)重的隨機(jī)噪聲也會影響處理效果,甚至使有的處理無法下手,因此必須衰減。
現(xiàn)在處理軟件中所用的隨機(jī)噪聲衰減模塊均是基于預(yù)測原理:在F-X域(對于三維是F-X-Y域)中,地震數(shù)據(jù)將分為沿X方向可預(yù)測和不可預(yù)測兩部分,相干同相軸(包括信號和規(guī)則噪聲)是可預(yù)測的,而隨機(jī)噪聲是不可預(yù)測的。于是隨機(jī)噪聲衰減的做法是對輸入的T-X域數(shù)據(jù)的各道做傅氏變換,得到F-X域數(shù)據(jù);對每個(gè)頻率成分沿X方向做預(yù)測濾波;變換回T-X域即得隨機(jī)噪聲衰減后的結(jié)果。
疊前去噪的重點(diǎn)是相干噪聲的壓制,本文所述的幾種相干噪聲的壓制方法各有優(yōu)點(diǎn):F-X域相干噪聲壓制和F-K域?yàn)V波使用較為方便,為實(shí)際資料處理中的首選方法;τ-P域去噪在衰減多次波方面有其獨(dú)到的手段;而多道傾角濾波除了可用于疊前相干噪聲壓制以外,還可用于疊后的修飾性處理。需要說明的是,無論哪種相干噪聲壓制方法都不可避免地或多或少要傷害到有效信號,因此不宜過量,否則,雖然疊前記錄上的相干噪聲看起來去得很干凈,而疊加結(jié)果反而不好。疊前隨機(jī)噪聲的壓制不是必須的。目前野外采集都是多次覆蓋,疊加是壓制隨機(jī)噪聲的最好辦法。對于復(fù)雜地區(qū)的低信噪比資料,由于有效信號的空間分布缺乏同相性,可能被當(dāng)做隨機(jī)噪聲而被衰減。如果非做不可,也要適度。