孟凡可， 安黎明， 李永康

(1.華北水利水電大學(xué) 地球科學(xué)與工程學(xué)院, 鄭州 450046; 2.河南省航空物探遙感中心, 鄭州 450002)

多次波是海洋地震勘探中最突出的問題之一，多次波的存在嚴重影響地震波成像的真實性和可靠性，如何有效壓制多次波并突出有效波是海洋地震資料處理的關(guān)鍵之一[1]。多次波的衰減方法大致可以分為兩大類[2]。在國外，最開始由Robinson[3]引入維納濾波來壓制多次波,這種方法對于1-D,特別是在線性Radon域效果明顯，Backus[4]利用預(yù)測反褶積來消除海上鳴震，Alam等[5]提出了對常規(guī)中點(common middle point，CMP)道集進行傾斜疊加變換的方法。目前在國內(nèi)多采用預(yù)測反褶積和表面相關(guān)多重消除(surface related multiple elimination，SRME)對多次波進行壓制，并能取得不錯的效果。但二者均具有一定局限性，反褶積是建立在子波是最小相位的時不變的、高信噪比、白噪假設(shè)等3個假設(shè)的基礎(chǔ)上的，并且在時間域隨著偏移據(jù)的增大，多次波的周期性減弱[6-8]。SRME則嚴重依賴于近偏移距信息并且同樣要求數(shù)據(jù)要有較高的信噪比[9-10]。為了解決常規(guī)“串聯(lián)”式衰減多次波帶來的弊端，本文開發(fā)了聯(lián)合多次波建模壓制技術(shù)。

1 鶯歌海盆地多次波概況

鶯歌海盆地早期勘探開發(fā)主要集中在1 200 ms左右的鶯歌海盆地東方某氣田，該區(qū)域壓力主要隨埋深變化，為常規(guī)壓力地質(zhì)環(huán)境[7]。隨著鶯歌海盆地勘探開發(fā)往中深層復(fù)雜區(qū)域黃流組、梅山組開展，開始面臨高溫高壓地質(zhì)環(huán)境的挑戰(zhàn)[8]。經(jīng)研究發(fā)現(xiàn)，鶯歌海盆地沉降速率很快,沉積物的快速堆積使得巖石極易形成欠壓實，沉積物中的大量流體無法排出，在上覆地層的壓力下，孔隙流體形成超壓，使速度降低，造成速度倒轉(zhuǎn)(圖1)。有效地層與多次波在速度上基本難以區(qū)分，對多次波的壓制帶來很大的挑戰(zhàn)(圖2)。如東方某氣田自總體開發(fā)方案(overall development plan,ODP)實施前后儲量變化大，主要原因在于后期水道改造，儲層橫縱向變化快，橫向看似連通，但縱向多期疊置關(guān)系復(fù)雜，之前處理結(jié)果與鉆井結(jié)果不吻合。而東方區(qū)、樂東區(qū)等區(qū)塊儲層層系相似、沉積模式可類比，對于該區(qū)的地震解釋來說微小的變化對該區(qū)的解釋開發(fā)方案都有著巨大的影響。因此多次波衰減在鶯歌海盆地高溫高壓區(qū)地震資料處理中起到了關(guān)鍵的作用[11-12]。

圖1 鶯歌海盆地樂東區(qū)典型多次波特征剖面

圖2 鶯歌海盆地樂東區(qū)典型速度倒轉(zhuǎn)特征

2 聯(lián)合多次波建模匹配衰減技術(shù)

短程多次波是地震波從某一深部界面反射回來再在海面(或地面)向下反射，然后又在某一個較淺的界面發(fā)生反射所形成的多次波[13]。短程多次波會使一次反射波的振幅、頻率和相位發(fā)生畸變，而且較難發(fā)現(xiàn)和識別。所以，在確定地震記錄上的同相軸波形時，準確地識別短程多次波顯得非常重要。通過研究分析得知，鶯歌海盆地高溫高壓區(qū)自由表面多次波中又以水層多次波干擾為主，剩余自由表面多次波干擾和層間多次波等類型[14-16]。本文主要開展針對鶯歌海盆地高溫高壓區(qū)水層多次波、自由表面多次波及層間多次波模型預(yù)測方法的研究[17-19]。

2.1 原理及方法

當前常規(guī)的多次波衰減主要采用淺水多次波衰減SDM，波動方程SRME、三維褶積SRME，以及各種技術(shù)組合衰減多次波，即“串聯(lián)式”多次波壓制，已經(jīng)可以相對較好地建立多次波模型衰減多次波，但是從理論上來說，仍然存在著一定的弊端(圖3)。新開發(fā)的多次波聯(lián)合建模匹配衰減技術(shù)(JMPAS)，既是一項針對技術(shù)，也是一種處理思想：基于原始波場盡可能地預(yù)測出準確的各種類型多次波，通過一次性匹配減去對多次波進行最大程度衰減，同時避免了多次匹配而產(chǎn)生的噪聲[20-21]。

圖3 串聯(lián)模式下衰減多次波存在的弊端分析

假設(shè)地質(zhì)模型為水平層狀介質(zhì)，自由表面為吸收邊界，假設(shè)x0(t)為不含自由表面多次波的信號，當x0(t)在自由界面向下發(fā)射時，該信號可以看作一個新的震源，這個過程也可以描述為自由表面反射產(chǎn)生的脈沖和整個脈沖序列褶積產(chǎn)生一階多次波，負號代表自由表面向下的反射。

設(shè)x(t)為含有地表相關(guān)多次波的總體響應(yīng)，通過反饋模型可以推出：

x(t)=x0(t)1-x0(t)2+x0(t)3-…

(1)

通過反饋模型還可得到脈沖響應(yīng)和所有多次波的隱式關(guān)系。設(shè)δ(t)代表震源函數(shù)及下行響應(yīng)，x(t)代表總體反射響應(yīng)，x0(t)表示來自介質(zhì)的響應(yīng)，則x(t)可以表示為

x(t)=x0(t)[δ(t)-x(t)]=x0(t)-x0(t)x(t)

(2)

式(2)可以理解為自由表面發(fā)生的一次反射響應(yīng)和總體響應(yīng)進行褶積后就可以預(yù)測出地表相關(guān)多次波。在極淺水域，實際接收到的水底響應(yīng)模糊，不利于水體多次波模型的建立。

假設(shè)x(t)不含水體響應(yīng)，設(shè)h(t)代表水體響應(yīng)，則x(t)可以表示為

x(t)=x0(t){δ(t)-[x(t)+h(t)]}=

x0(t)-x0(t)[x(t)+h(t)]

(3)

式中，

(4)

將式(1)變換到頻率域，得

(5)

式(5)可以表示為

x(f)=x0(f)-x0(f)x(f)

(6)

式(6)可以表示為

x(f)=x0(f){1-[x(f)+Δ(f)]}

(7)

上式給出了一次反射波和整個數(shù)據(jù)之間的關(guān)系，寫成x0的表達式，則為

x0(f)=x(f){1-[x(f)+Δ(f)]}-1

(8)

這就是利用數(shù)據(jù)的總體響應(yīng)與海底響應(yīng)來計算去除多次波后的響應(yīng)方程。

圖4 聯(lián)合多次波建模思路示意圖

聯(lián)合多次波建模匹配衰減技術(shù)(JMPAS)建模思路如圖4所示，其具有以下幾點技術(shù)優(yōu)勢：①多種傳播路徑的多次波模型預(yù)測更加準確；②該方法主要是基于海底模型與數(shù)據(jù)本身聯(lián)合預(yù)測；③聯(lián)合衰減在最大限度壓制多次波的前提下提高保真度。

根據(jù)原理通過正演對聯(lián)合匹配衰減技術(shù)進行論證?；陴椥苑匠淘O(shè)計并正演一個淺水含多次波地震模型數(shù)據(jù)，模型淺層存在低速異常?；诖四Ｐ蛯Χ啻尾ǚ椒ㄟM行測試，通過測試對比可以看到常規(guī)的串聯(lián)方法多次波衰減殘留較多，通過聯(lián)合多次波建模衰減多次波殘留較少，如圖5所示。

圖5 串聯(lián)模式與聯(lián)合建模在地震模型中的應(yīng)用對比

2.2 實際數(shù)據(jù)應(yīng)用

以下為多次波聯(lián)合建模匹配衰減技術(shù)(JMPAS)在鶯歌海盆地東方區(qū)和樂東區(qū)的實際應(yīng)用效果。圖6為東方區(qū)原始疊加，從剖面上來看多次波信息掩蓋了目的層的有效信號，無法看清具體形態(tài)和地質(zhì)特征。圖7可以看到通過串聯(lián)方法衰減多次波后，目標段地質(zhì)特征清晰，多次波基本衰減干凈，但仍存在一些殘留多次波。圖8通過聯(lián)合多次波建模衰減后整個剖面多次波衰減干凈，多次波殘留多次波少。

圖6 鶯歌海東方某場區(qū)實際地震資料原始波場疊加

從圖9可以看出，聯(lián)合建模技術(shù)道集上的多次波壓制效果明顯優(yōu)于串聯(lián)模式方法，多次波衰減徹底，有效信號的保真性更高。泥巖層段的多次波干擾也明顯沒有殘留。從圖10可以看出，聯(lián)合建模衰減后的自相關(guān)譜分析多次波衰減更好，頻譜上受多次波干擾造成的頻譜周期性的“抖動”更少。從圖11可以看出，在樂東強超壓區(qū)應(yīng)用新開發(fā)的聯(lián)合建模技術(shù)，疊加剖面上中深層多次波被明顯壓制后，真實構(gòu)造得以恢復(fù)，構(gòu)造樣式和地層的可解釋性增強，T51、T60界面更容易追蹤解釋。三亞組砂體振幅特征明顯。崖城組不整合接觸關(guān)系清晰。從圖12可以看出，道集上的多次波壓制效果明顯優(yōu)于串聯(lián)模式方法，近道殘留多次波衰減徹底，有效信號的保真性更高。

圖10 鶯歌海東方某場區(qū)實際地震資料多次波壓制方法自相關(guān)及平譜對比

圖11 鶯歌海樂東區(qū)地震資料多次波壓制方法疊加對比

圖12 鶯歌海樂東區(qū)地震資料多次波壓制方法道集對比

3 結(jié)論

鶯歌海盆地中深層高溫高壓區(qū)速度倒轉(zhuǎn)，主要目的層黃流組、梅山組地層受多次波干擾現(xiàn)象嚴重。之前的處理技術(shù)并沒有專門針對中深層高溫高壓區(qū)復(fù)雜多次波的壓制技術(shù)，中深層的處理效果一直不太理想。隨著勘探開發(fā)逐漸轉(zhuǎn)向中深層高溫高壓區(qū)，速度倒轉(zhuǎn)區(qū)的多次波干擾必然要被重視。針對鶯歌海淺水區(qū)高溫高壓地震資料的特征，首先，以自由表面多次波為重點，針對二維SRME、傳統(tǒng)三維SRME的不足，開發(fā)了基于反饋環(huán)理論的三維廣義SRME多次波預(yù)測壓制技術(shù)。然后針對常規(guī)“串聯(lián)式”多次波壓制技術(shù)的不足，新開發(fā)了基于原始波場的多路徑多次波聯(lián)合建模匹配衰減技術(shù)(JMPAS)。該技術(shù)既是一項針對技術(shù)，也是一種處理思想：基于原始波場盡可能地預(yù)測出準確的各種類型多次波，通過一次性匹配減去對多次波進行最大程度衰減，同時避免了多次匹配而產(chǎn)生的噪聲。通過該技術(shù)在鶯歌海盆地高溫高壓區(qū)東方區(qū)和樂東區(qū)的應(yīng)用，取得了相當顯著的效果，為中深層的勘探開發(fā)提供了高質(zhì)量的地震資料。