宋家文，陳小宏

(中國石油大學(xué)(北京)油氣資源與探測國家重點實驗室，北京102249)

地下介質(zhì)中的強反射界面會產(chǎn)生明顯的層間多次波，嚴重干擾目標層有效信號，并在成像剖面中產(chǎn)生虛假界面同相軸，甚至與真實的界面同相軸重合，導(dǎo)致錯誤的地震反演和解釋結(jié)果。因此，多次波通常被當作噪聲，在地震數(shù)據(jù)疊前處理階段進行壓制。目前，多次波壓制方法主要分為兩類[1]：一類是基于多次波與一次波特征和性質(zhì)差異的濾波法;另一類是基于波動方程的預(yù)測相減法。其中，工業(yè)界廣泛應(yīng)用的自由表面多次波壓制方法(Surface Related Multiple Elimination,SRME)[2-5]屬于后者。

與自由表面多次波不同，層間多次波的下行反射全部發(fā)生在地表或海面以下。當?shù)叵陆橘|(zhì)中存在鹽丘頂、底或海底等強反射界面時，地震數(shù)據(jù)中會產(chǎn)生明顯的層間多次波。相對于平坦的海水表面，海底或地下界面崎嶇不平，因此，層間多次波的形態(tài)往往比自由表面多次波更加復(fù)雜[6]。此外，由于在地層界面之間經(jīng)歷了多次反射，層間多次波與一次波的速度差異較小，基于多次波周期性與可分離性的濾波方法很難取得壓制效果。因此，波動方程預(yù)測相減法成為衰減層間多次波的首選。

反饋迭代法是一種基于波動方程的多次波預(yù)測相減方法，它以描述地震數(shù)據(jù)正演過程的反饋模型為理論基礎(chǔ)，利用地震波在自由表面、層界面間的傳播規(guī)律來預(yù)測多次波[7]，能夠適應(yīng)復(fù)雜的地下結(jié)構(gòu)。針對層間多次波，反饋迭代法一般采用層剝離的方式，按照反射界面由淺至深的順序，遞歸消除各界面的下行反射效應(yīng)。實際地下介質(zhì)中，能夠產(chǎn)生明顯層間多次波的強反射界面數(shù)量有限，因此可以在疊加剖面上人工選取基準面，把強反射界面劃分到不同的反射層，從而以反射層為基本物理單元預(yù)測層間多次波，然后通過自適應(yīng)相減實現(xiàn)層間多次波的壓制。

目前，反饋迭代層間多次波壓制方法主要分為兩種：基于延拓數(shù)據(jù)的共聚焦點(CFP)方法[8-9]與基于地表數(shù)據(jù)分離的Jakubowicz方法[10-11]。CFP方法需要將地表波場反向延拓至地下界面，構(gòu)建炮點在地表，檢波點在地下的CFP道集，以及炮點和檢波點分別位于地下的網(wǎng)格點道集，計算量較大，對地下速度模型有一定程度的依賴。相比之下，基于地表數(shù)據(jù)分離的方法直接利用地表觀測數(shù)據(jù)中分離出的三部分波場數(shù)據(jù)，進行地表一致性相關(guān)與褶積運算來預(yù)測層間多次波，從而避免了地表數(shù)據(jù)的反向延拓，具有計算效率高，不依賴地下速度的優(yōu)勢。

我們首先闡述基于地表數(shù)據(jù)分離的層間多次波壓制方法原理，并通過模擬算例展示該方法壓制層間多次波的具體步驟;然后針對地下介質(zhì)中存在多個強反射界面的地震數(shù)據(jù)，重點論述逐層壓制層間多次波的遞歸實現(xiàn)過程;最后用實際地震資料驗證了方法的有效性。

1 方法原理

二維地震數(shù)據(jù)可以排列成以炮點位置、檢波點位置和旅行時為坐標的空間-時間域數(shù)據(jù)體，沿時間方向作傅里葉變換，可以將數(shù)據(jù)體變換到頻率域，其中單頻切片對應(yīng)的矩陣稱為數(shù)據(jù)矩陣，一般用P表示[12]。Kinneging等[13]發(fā)展了三維地震數(shù)據(jù)的數(shù)據(jù)矩陣表示方法，極大地方便了地震資料處理技術(shù)的三維推廣。我們闡述的層間多次波預(yù)測方法是基于數(shù)據(jù)矩陣的，因此可適用于二維和三維地震數(shù)據(jù)。

在海洋地震勘探中，海水表面的下行反射作用會產(chǎn)生強烈的自由表面多次波，因此，地震數(shù)據(jù)可表示為一次波P0(含有層間多次波)與自由表面多次波M0之和，即

(1)

Verschuur等[2]提出了衰減自由表面多次波的SRME方法，該方法通過引入表面算子A，得到了完全數(shù)據(jù)驅(qū)動的自由表面多次波預(yù)測公式：

(2)

其中,A=S-1R-，S表示震源子波，R-為海水表面的下行反射系數(shù)矩陣，一般R-=-I。

圖1描述了公式(2)所示的自由表面多次波預(yù)測過程，由炮點k激發(fā)，檢波點j接收的自由表面多次波，可通過紅線所示的共炮點道集P與藍線所示的共檢波點道集P0，沿地表進行一致性褶積后疊加得到，表面算子A的作用相當于對預(yù)測得到的自由表面多次波進行震源子波反褶積。實際應(yīng)用中子波未知，預(yù)測得到的多次波與真實地震記錄存在振幅或相位上的差異，因此，必須從地震記錄中自適應(yīng)減去，以達到衰減自由表面多次波的目的。

圖1 自由表面多次波預(yù)測示意圖解

按照界面由淺至深的順序逐步消除多次波的過程可描述為

(3)

其中，Pn-1表示消除界面z0≤z≤zn-1產(chǎn)生的多次波以后的地震數(shù)據(jù);類似地，Pn表示消除由界面z0≤z≤zn產(chǎn)生的多次波;Mn表示與界面zn有關(guān)的層間多次波。注意，當n=0時，公式(3)等價于公式(1)，即自由表面多次波可以看作層間多次波的一種特例，因此可以將自由表面多次波的預(yù)測公式(2)加以推廣，用來預(yù)測層間多次波。

Jakubowicz[10]提出層間多次波可由地表觀測數(shù)據(jù)中直接分離出的3部分波場數(shù)據(jù)，經(jīng)過地表一致性相關(guān)與褶積運算來預(yù)測，其預(yù)測公式為

(4)

圖2 層間多次波預(yù)測示意圖解

(5)

2 模擬算例

我們通過二維模擬數(shù)據(jù)來展示壓制層間多次波的具體步驟。圖3a給出了正演模擬所用的鹽丘速度模型，海底和鹽丘頂界反射系數(shù)較大，是兩個主要的層間多次波下行反射界面。模擬數(shù)據(jù)通過固定排列采集，炮點和檢波點間隔均為15m，總共模擬了121炮，每炮121道。正演模擬時，通過加吸收邊界，地表產(chǎn)生的自由表面多次波被衰減，因此模擬數(shù)據(jù)中只包含一次波與層間多次波。圖3b展示了模擬數(shù)據(jù)的零偏移距剖面，其中，箭頭所示的同相軸是由海底下行反射產(chǎn)生的層間多次波，其振幅甚至高于深層的一次反射;尤其是1.0s附近的層間多次波與有效信號發(fā)生了嚴重干涉，會降低后續(xù)反演與解釋結(jié)果的可靠性。

圖3 鹽丘速度模型(a)和模擬數(shù)據(jù)的零偏移距剖面(b)

由海底下行反射產(chǎn)生的層間多次波可按照以下步驟消除。首先，在零偏移距剖面(圖3b)上拾取如黃色虛線所示的基準面z1。然后，利用基準面t0時間與上覆地層NMO速度，在疊前炮集上分離基準面上、下地層的反射波場。其中，波場分離方法主要有兩種：一種是直接利用基準面在疊前炮集上對應(yīng)的時距雙曲線進行簡單切除，適用于基準面上、下地層反射同相軸不交叉的情況;另一種是首先通過f-k域反向波場延拓，將地表波場延拓至基準面，然后以零時間為界分離延拓波場，并分別正向延拓至地表，得到基準面上、下地層的反射波場，采用該方法可以分離相互交叉的同相軸。圖4a 展示了模擬數(shù)據(jù)的一個炮集記錄，利用波場延拓的方法將原始炮集分離成上、下兩部分，得到如圖4b與圖4c所示的分別來自基準面z1上、下地層界面的反射波場，可以看出，海底繞射和遠偏移距的交叉同相軸也得到了很好的分離。

最后，按照Jakubowicz方法的預(yù)測公式(5)，傳播路徑穿過基準面z1至少4次的層間多次波，可由基準面以下的反射波場與基準面以上的反射波場二維互相關(guān)得到的結(jié)果，再與基準面以下的反射波場進行一致性褶積來預(yù)測。圖5a展示了預(yù)測得到的與基準面z1相關(guān)的層間多次波，通過對比可以發(fā)現(xiàn)，其走時和相位與原始炮集(圖4a)非常吻合，但是存在振幅上的差異。通過常規(guī)最小二乘自適應(yīng)相減得到如圖5b所示的衰減掉與基準面z1相關(guān)的層間多次波后的炮集記錄，可以看出，深層的有效信號恢復(fù)效果理想，層間多次波得到了較好的消除(圖5c)。

圖4 模擬數(shù)據(jù)原始炮集(a)和分離出的來自基準面z1以上(b)及以下(c)的反射波場

圖5 預(yù)測得到的與z1相關(guān)的層間多次波(a)和與z1相關(guān)層間多次波衰減后的炮集(b)以及消除掉的層間多次波(c)

圖6a是衰減了與基準面z1相關(guān)的層間多次波后的零偏移距剖面，與模擬原始數(shù)據(jù)(圖3b)對比，可以清楚地看到深層有效信號得到了明顯的恢復(fù)，而圖3b中兩個黃色箭頭所示的層間多次波得到了有效地衰減。但仔細觀察衰減層間多次波后的結(jié)果(圖6a)可以發(fā)現(xiàn)，1.3s附近(圖6a中紅色箭頭所示)的層間多次波仍然存在于地震數(shù)據(jù)中，因為它是由鹽丘頂界面下行反射產(chǎn)生的，其傳播路徑經(jīng)過基準面z1的次數(shù)少于4次。因此，對于地下介質(zhì)中存在多個強反射界面的地震數(shù)據(jù)，需要按照界面由淺至深的順序拾取多個基準面，依次遞歸消除層間多次波。

鹽丘頂界面產(chǎn)生的層間多次波可通過3步進行預(yù)測：①在零偏移距剖面上鹽丘頂、底界面之間選取新的基準面z2(圖6a中紅色虛線所示);②以衰減了與基準面z1相關(guān)層間多次波后的炮集(圖5b)為輸入，分離介于基準面z1與z2之間的反射同相軸，以及z2以下的反射波場;③對分離出的數(shù)據(jù)沿地表進行一致性相關(guān)與褶積，預(yù)測與z2有關(guān)的層間多次波，并從地震記錄中自適應(yīng)減去。圖6b 展示了遞歸消除與z2有關(guān)層間多次波后的零偏移距剖面。對比圖6a與圖6b 可以發(fā)現(xiàn)，鹽丘頂界產(chǎn)生的1.3s附近的層間多次波已經(jīng)被消除，在CMP道集(圖6c)上清晰顯示了這部分能量。

圖6 消除了與z1有關(guān)層間多次波后的零偏移距剖面(a)和遞歸消除了與z2有關(guān)層間多次波后的零偏移距剖面(b)以及遞歸消除掉的層間多次波(c)

圖7展示了模擬數(shù)據(jù)層間多次波衰減前、后的CMP道集。原始CMP道集(圖7a)中除了0.5s

附近振幅較強的層間多次波外，1.0s以下的有效信號亦與層間多次波嚴重干涉在一起，難以分辨;

圖7 衰減層間多次波前(a)、后(b)的模擬數(shù)據(jù)CMP道集

而經(jīng)過2次層間多次波衰減后的CMP道集(圖7b)中，深層的弱一次反射得到了很好的恢復(fù)。對比多次波衰減前、后的速度譜(圖8a與圖8b)可以發(fā)現(xiàn)，層間多次波對應(yīng)的能量團(黑色箭頭所示)被消除，原本淹沒在層間多次波中的有效信號(紅色箭頭所示)得以突顯。

圖8 層間多次波衰減前(a)、后(b)速度譜

3 實際資料處理

在模擬算例中，基于地表數(shù)據(jù)分離的層間多次波壓制方法取得了理想的效果。我們將該方法應(yīng)用到墨西哥灣實際地震資料中，驗證其有效性。圖9a 是墨西哥灣深水資料的原始疊加剖面，自由表面多次波出現(xiàn)在3.5s以后，因此，這里我們只考慮壓制層間多次波。

在圖9a展示的原始疊加剖面上，1.7～2.3s存在一系列反射界面，其中海底與鹽丘頂界的反射系數(shù)較大，可以產(chǎn)生大量的層間多次波。選取介于強反射界面之間的基準面z1與z2，可將地下反射界面分成3個反射層，遞歸應(yīng)用基于地表數(shù)據(jù)分離的層間多次波壓制方法可依次消除由反射層z0

圖9 墨西哥灣實際地震資料基于地表數(shù)據(jù)分離的層間多次波壓制a 原始疊加剖面; b 衰減層間多次波后的疊加剖面; c 衰減掉的層間多次波

4 結(jié)束語

基于地表數(shù)據(jù)分離的層間多次波壓制方法以層界面為基本物理單元預(yù)測層間多次波，不依賴地下速度，能適應(yīng)復(fù)雜地層。該方法的本質(zhì)與SRME相同，只是增加了一次互相關(guān)運算，計算量約是SRME壓制自由表面多次波的兩倍，具有較高的工業(yè)應(yīng)用價值。當?shù)叵陆橘|(zhì)中存在多個強反射界面時，基于單一基準面的預(yù)測方法會產(chǎn)生多次波泄漏，因此需要合理選擇多個基準面對地下反射界面進行層劃分，按照由淺至深的順序逐層遞歸地消除其下行反射產(chǎn)生的層間多次波。模擬數(shù)據(jù)與墨西哥灣實際資料的層間多次波遞歸壓制處理結(jié)果表明，該方法可以有效地衰減層間多次波，恢復(fù)有效信號。

參 考 文 獻

[1] Weglein A.Multiple attenuation:an overview of recent advances and the road ahead (1999)[J].The Leading Edge,1999,18(1):40-44

[2] Verschuur D J,Berkhout A J,Wapenaar C P A.Adaptive surface-related multiple elimination[J].Geophysics,1992,57(9):1166-1177

[3] Berkhout A J,Verschuur D J.Estimation of multiple scattering by iterative inversion,part I:theoretical considerations[J].Geophysics,1997,62(5):1586-1595

[4] Verschuur D J,Berkhout A J.Estimation of multiple scattering by iterative inversion,part II:practical aspects and examples[J].Geophysics,1997,62(5):1596-1611

[5] Dragoset B,Verschuur E,Moore I,et al.A perspective on 3D surface-related multiple elimination[J].Geophysics,2010,75(5):75A245-75A261

[6] Verschuur D J.Seismic multiple removal techniques-past,present and future[M].Houten:EAGE Publications BV,2006:1-191

[7] Verschuur D J,Matson K H,Berkhout A J,et al.A comparison of the feedback and inverse scattering internal multiple attenuation methods[J].Expanded Abstracts of 61stAnnual Conference,1999,1-5

[8] Berkhout A J,Verschuur D J.Removal of internal multiples with the common-focus-point (CFP) approach:part 1—explanation of the theory[J].Geophysics,2005,70(3):V45-V60

[9] Verschuur D J,Berkhout A J.Removal of internal multiples with the common-focus-point (CFP) approach:part 2—application strategies and data examples[J].Geophysics,2005,70(3):V61-V72

[10] Jakubowicz H.Wave equation prediction and removal of interbed multiples[J].Expanded Abstracts of 68thAnnual Internat SEG Mtg,1998,1527-1530

[11] 馬繼濤,陳小宏,姚逢昌,等.基于一次波逆時理論的層間多次波衰減方法[J].石油地球物理勘探,2013,48(2):181-186

Ma J T,Chen X H,Yao F C,et al.Internal multiple attenuation using time-inverted primary[J].Oil Geophysical Prospecting,2013,48(2):181-186

[12] Berkhout A J.Seismic migration,imaging of acoustic energy by wave field extrapolation,A:theoretical aspects[M].3rded.Amsterdam:Elsevier Science Limited,1982:1-445

[13] Kinneging N K,Budejicky V,Wapenaar C P A,et al.Efficient 2D and 3D shot record redatuming[J].Geophysical Prospecting,1989,37(5):493-530