張書杰，薛霆虓

(1.桂林理工大學(xué)地球科學(xué)學(xué)院，廣西 桂林 541004；2.廣西隱伏金屬礦產(chǎn)勘查重點實驗室，廣西 桂林 541004)

0 引 言

地震偏移成像技術(shù)一般分為有基于波場延拓的波動方程偏移方法和基于射線理論的Kirchhoff偏移方法兩類。波動方程偏移方法是對波動方程的進行分解，能較好的適應(yīng)速度模型的橫向變化，但由于采用單程波的近似方程，在波場傳播方向受到傾角限制，對陡傾角構(gòu)造成像能力較差，尤其對于鹽丘下回折波成像存在很大的困難。Kirchhoff偏移方法有很高的計算效率和實用性，在當(dāng)前實際生產(chǎn)中廣泛應(yīng)用，但Kirchhoff偏移方法的分辨率會隨著深度增加變差且缺乏相應(yīng)的振幅信息。隨著油氣勘探難度的增大和目標(biāo)復(fù)雜性的增加，這時迫切需要發(fā)展一種新的偏移方法來彌補這些不足。

逆時偏移方法是基于雙程波波動方程，通過對波場正向和反向延拓來實現(xiàn)偏移。逆時偏移是較Kirchhoff偏移方法而言能有效利用全波場信息，較單程波偏移方法而言，無傳播方向限制，不需要對波場進行上、下行波的分離，而且不受傾角限制，并能較好的成像回轉(zhuǎn)波、多次波，理論上在給定的速度場足夠準(zhǔn)確的情況下能對任意復(fù)雜介質(zhì)精準(zhǔn)成像，是一種具有極高成像精度的偏移方法。

1 逆時偏移噪音的產(chǎn)生機制

逆時偏移有著與生俱來的三大缺陷：計算量大、存儲量大和存在虛假成像。如何有效壓制虛假的噪聲是獲得正確逆時偏移成像結(jié)果的關(guān)鍵所在。

根據(jù)互相關(guān)成像條件公式

(1)

t=t1+t2

只要滿足：震源正向延拓波場的旅行時間t1與檢波點反向延拓波場的旅行時間t2之和，等于在從震源出發(fā)至檢波點的總旅行時間t，就滿足互相關(guān)成像條件，因此都可以成像。如圖1所示的兩條射線的路徑上，任何一點的走時之和都是等于t的，因此在這兩條射線的每個點上，都滿足互相關(guān)成像條件，在偏移時都會成像，而真實的成像點只在反射界面處，這就是低頻噪音產(chǎn)生的原因。

方框處A是震源正向延拓波場的與檢波點反向延拓波場夾角等于180°的成像噪音產(chǎn)生的位置，圓圈處B是震源正向延拓波場的與檢波點反向延拓波場夾角小于180°的成像噪音產(chǎn)生的位置。

圖1 逆時偏移低頻噪音產(chǎn)生的原因

2 基于波場分離的逆時偏移去噪方法

Liu[1]提出通過一種波場分離的方法，提取所需要的波場分量，只用在有效路徑上進行互相關(guān)成像，來達到去除噪音的目的。其主要思路就是直接將波場拆分，就可以對波場分量進行人工選擇，來壓制反射波能量。利用波場的傳播路徑將真實的成像點和虛假的成像點區(qū)別開來，震源波場和檢波點波場在相反的傳播方向所成的像才是真實的界面，而在相同的波場方向所成的像是虛假成像，如果只在真實的成像點處使用成像條件進行成像，這樣就從成因上消除了低頻噪音。

在逆時偏移中使用的是雙程波動方程，它包含了上行波和下行波兩個組成部分，如果在z方向進行上、下行波的分解，則震源波場和檢波點波場可分解為

S(x，t)=Sz+(x，t)+Sz-(x，t)

(2)

R(x，t)=Rz+(x，t)+Rz-(x，t)

(3)

其中，“+”號代表下行波場，“-”號代表上行波場。將方程(2)和(3)代入互相關(guān)成像條件中(1)中，得到新的成像條件(4)

=Iz1+Iz2+Iz3+Iz4

(4)

經(jīng)過真實的反射界面時，入射波場和反射波場的傳播方向在界面的法線方向會發(fā)生變化。因此震源波場和檢波點波場在反射界面上都會產(chǎn)生相對于界面向上傳播的波和相對于界面向下傳播的波。按照震源波場和檢波點波場的傳播方向，將其波場分離成其各自的傳播方向不同的波場，只有兩者傳播方向相反的波場才會在成像時產(chǎn)生真實的成像點。也就是說，只需要Ι21和Ι22兩項就可以實現(xiàn)正確的逆時偏移成像。在相同傳播方向所成的像是虛假的成像點，因此Ι23和Ι24這兩項其成像結(jié)果是低頻噪音，所以應(yīng)該將其刪去。這樣便可以將(4)式改寫成新的成像條件(5)式：

(5)

為了成功的應(yīng)用這個成像條件，需要將震源波場和檢波點波場拆分成上行波場分量和下行波場分量。本文采用的是二維傅里葉變換的方法，先將震源和檢波點波場轉(zhuǎn)換到頻率波數(shù)域再對其進行上、下行波的分離，然后通過反變換到時間域，這樣就可以成功的利用波場分離的方法來的壓制逆時偏移噪音。

3 模型計算

圖2是Marmousi 2模型的速度模型示意圖。該模型上部有低速覆蓋層，地層中有含氣透鏡體模型和油氣儲層構(gòu)造，模型中部有三條大的斷層，下部有高速的橫向侵入巖，模型中部陡傾有斷層下的背斜儲層構(gòu)造。Marmousi2模型的地層數(shù)有199個，眾多的反射界面，陡傾斷層，強烈的橫向、縱向速度變化和內(nèi)部的背斜構(gòu)造都是成像的難點。

圖2 Marmousi 2速度模型示意圖

圖3為去噪前Marmousi 2模型逆時偏移成像結(jié)果，低頻噪音基本覆蓋整個模型，尤其模型上部高振幅的假象最為嚴重，遮掩了下部構(gòu)造，影響了成像質(zhì)量。

圖3 去噪前Marmousi 2模型成像結(jié)果

圖4是波場分離法去噪后成像效果圖。使用波場分離方法后成像噪音壓制情況良好，圖像的低頻噪音的影響基本消除。去噪之后各地層反射界面顯示準(zhǔn)確，含氣透鏡體模型、三條大的斷層和下部的高速的橫向侵入巖更清晰，地層特征細節(jié)更為明顯。

圖4 波場分離法去噪后成像效果圖

4 小 結(jié)

通過波場分離方法在Marmousi 2標(biāo)準(zhǔn)模型上的測試，證明了該方法的有效性，由于在成因上消除了噪音，因此比其他去噪方法更具優(yōu)勢，所以該方法在油氣資料的地震勘探中有著較好的實用價值。但是要運用兩次二維傅里葉變換一定程度上相應(yīng)的也增加了計算量和存儲量。

單一的去噪方法并不能滿足各方面需求，根據(jù)不同的實現(xiàn)環(huán)境和成像條件，綜合考慮不同方法的優(yōu)缺點，對噪音壓制方法進行合理的搭配，構(gòu)建在不同實現(xiàn)條件和成像需求情況下逆時偏移去噪方法才是最佳抉擇策略[2]。

參 考 文 獻

[1]Liu F Q，Zhang G Q，Morton S A，et al.Reverse time migration using one-way wave field imaging condition[A].In：77rd Annual International Meeting[C].SEG Expanded Abstracts，2007，2170～2173.

[2]杜啟振，朱釔同，張明強，等.疊前逆時深度偏移低頻噪聲壓制策略研究[J].地球物理學(xué)報，2013，56(7)：2391～2401.

[3]丁亮，劉洋.逆時偏移成像技術(shù)研究進展[J].地球物理學(xué)進展，2011，26(3)：1085～1100.

[4]許璐，孟小紅，劉國峰.逆時偏移去噪方法研究進展[J].地球物理學(xué)進展，2012，27(4)：1548～1556.