明君，王建立，劉建輝，王志勇，李潔麗，劉國昌

1.中海石油（中國）有限公司天津分公司，天津 300459 2.中國石油大學(xué)（北京）地球物理學(xué)院，北京 102249

地震繞射波是地下構(gòu)造細(xì)節(jié)和異常的反映，地下構(gòu)造的不連續(xù)點(diǎn)會產(chǎn)生繞射波，如裂縫、斷層、溶洞和地層尖滅等。在原始地震勘探數(shù)據(jù)中，既存在繞射波也存在反射波，反射波是地質(zhì)背景的綜合反映，繞射波則是地質(zhì)細(xì)節(jié)的小構(gòu)造的反映。相對反射波而言，繞射波的能量往往很弱，甚至難以識別，導(dǎo)致繞射信息常常被高能量的反射信息所覆蓋，而繞射波信息是提高地震勘探分辨率的重要信息載體。在地震資料處理解釋中，正確識別斷層、尖滅和小尺寸散射體等地質(zhì)不連續(xù)性是一個(gè)重要問題，繞射波包含有關(guān)產(chǎn)生繞射波的介質(zhì)結(jié)構(gòu)和組成的有價(jià)值信息。然而，在標(biāo)準(zhǔn)地震資料處理中，繞射波常被視為噪聲而被抑制[1-3]。把繞射波從原始地震勘探記錄中分離出來進(jìn)行單獨(dú)成像，可以加強(qiáng)小型地震異常體的成像，從而提高復(fù)雜斷塊的地震成像精度。

在過去的幾十年，許多地球物理學(xué)者對地震繞射波分離方法做了相關(guān)研究。1988年Kanasewich和Phadke[4]提出了一種適用于共偏移距道集的繞射波走時(shí)曲線直接拾取繞射波的方法來分離繞射波。1992年Claerbout[5]首次提出了平面波預(yù)測濾波器的概念，并在信噪分離方面得到了很好的應(yīng)用。Fomel[6]進(jìn)一步完善了平面波預(yù)測濾波器，并提出了平面預(yù)測濾波器的殘差可以作為繞射波的近似。Nowak等[7]根據(jù)反射波和繞射波同相軸的幾何差異，采用加權(quán)雙曲Radon濾波器分離地震繞射波。Berkovitch等[8]利用繞射波多次聚焦疊加剖面將繞射波聚焦到繞射點(diǎn)位置，同時(shí)將反射波反聚焦到整個(gè)道集記錄，從而分離繞射波。Koren和Ravve[9]提出了方位相關(guān)反射波屬性，通過構(gòu)建加權(quán)疊加濾波器將繞射波和反射波分離。Klokov等[10]根據(jù)角度域共成像點(diǎn)道集中反射波和繞射波的差異，首先壓制去除反射頂點(diǎn)附近反射能量，然后利用混合Radon變換去除剩余反射波能量，進(jìn)而分離出繞射波。Moser和Howard[11]提出基于深度域的反射聚焦，從疊前數(shù)據(jù)中過濾反射進(jìn)而對繞射波分離成像的方法。Decker等[12]使用人工模擬數(shù)據(jù)和從皮森盆地獲得的野外實(shí)際數(shù)據(jù)進(jìn)行繞射波分離，評估解決繞射特征的多種方法，證明平面波預(yù)測濾波器方法能使圖像保存更多繞射能量和較少的遷移偽影，但不足是無法消除斜率不連續(xù)的反射變量。上述繞射波分離算法主要是利用了繞射波和反射波在不同地震數(shù)據(jù)道集中的差異，通過數(shù)學(xué)變換或者信號預(yù)測的方法將繞射波識別和分離出來。在國內(nèi)研究方面，趙娟娟等[13]利用F-K濾波方法分離成像繞射波，通過F-K域斜率濾波去除反射波從而分離出繞射波信息。馬永生[14]從油藏中具有代表性的地質(zhì)模型出發(fā)，分析VSP中繞射波的基本特點(diǎn)，理論上總結(jié)了繞射波運(yùn)動學(xué)的部分認(rèn)識規(guī)律。黃建平等[15]著重介紹了平面波預(yù)測濾波器及其在疊前域、疊后域分離并成像繞射波的方法原理，探討了繞射波分離成像方法的應(yīng)用前景和改進(jìn)方向。朱生旺[16]提出一種局部傾角濾波和預(yù)測反演聯(lián)合起來分離繞射波的方法。蔣波等[17]提出一種基于反射波層拉平的繞射波分離與成像方法用于碳酸鹽巖縫洞儲集體地質(zhì)目標(biāo)，提高對溶洞、裂縫等特殊異常體的刻畫精度。李正偉等[18]根據(jù)傾角域共成像點(diǎn)道集中菲涅耳帶變化特點(diǎn)，通過傾角-偏移距道集中精確切除菲涅耳帶的方法來分離繞射波。本文在平面波分解繞射波分離[6]基礎(chǔ)上，分析了繞射波與反射波信號特征的差異，深入研究了平面波預(yù)測的繞射波分離方法的影響因素，并探討了噪聲和平面波預(yù)測濾波器平滑參數(shù)對地震繞射波分離結(jié)果的影響，為在實(shí)際地震數(shù)據(jù)處理中應(yīng)用平面波分解方法進(jìn)行繞射波分離提供了理論指導(dǎo)。

1 平面波預(yù)測繞射波分離

地震記錄既含有反射波同相軸也含有繞射波同相軸，在點(diǎn)震源激發(fā)時(shí)，反射波和繞射波波形在炮記錄中都是雙曲線，不易于分離反射波和繞射波。在平面波震源激發(fā)時(shí)，來自反射界面的反射波是平面波，而來自繞射體的繞射波是雙曲線。因此，在平面波震源條件下反射波和繞射波具有明顯的幾何差異，更有利于繞射波分離。在實(shí)際地震勘探中大都是點(diǎn)源激發(fā)，反射波和繞射波的同相軸不易區(qū)分，但我們根據(jù)爆炸反射面理論將共偏移距剖面近似為平面波震源獲得的地震記錄，這樣繞射波和反射波有較好的幾何差異，根據(jù)其幾何差異進(jìn)行繞射波分離。

平面波預(yù)測濾波器利用未知的局部傾角構(gòu)造最優(yōu)非穩(wěn)態(tài)預(yù)測誤差濾波器，根據(jù)相鄰地震道的數(shù)據(jù)，在局部傾角平滑變化約束下，使預(yù)測誤差最小化以此預(yù)測目標(biāo)道數(shù)據(jù)。通過最優(yōu)化方法不斷迭代可獲得局部傾角場，傾角場能夠反映反射波場幾何信息。

局部平面波的微分方程表示為[6]：

其 中，P(t,x) 是 波 場，x、t分 別 是 炮 檢 距 和 旅 行 時(shí)，σ是局部傾角，是x和t的函數(shù)，局部傾角是旅行時(shí)對炮檢距的導(dǎo)數(shù)。

Fomel給出了Z變換域平面波預(yù)測的形式[6]：

利用平面波預(yù)測濾波器預(yù)測分離繞射波的過程中，關(guān)鍵一步是反射波同相軸局部傾角的估計(jì)。如果局部傾角已知的情況下，可以用C(σ)表示地震數(shù)據(jù)與二維濾波器的褶積算子。因此，在最小平方意義下，可通過求解如下最小二乘目標(biāo)函數(shù)估計(jì)局部傾角：

其中，d為已知的地震數(shù)據(jù)。優(yōu)化問題（3）可以利用高斯—牛頓線性化迭代的方法求解：

式中， ?σ 為傾角增量， σ0為 初始傾角，C(σ0)為對應(yīng)初 始 傾 角 σ0估 計(jì) 的 褶 積 算 子。C′(σ0)為C(σ0)的導(dǎo)數(shù)值。為了保證傾角場穩(wěn)定可靠，在反演求解過程中引入正則化約束項(xiàng)，控制傾角場沿著空間和時(shí)間兩個(gè)方向平滑。由于平面預(yù)測是針對平面波形的預(yù)測，而繞射波難于滿足其預(yù)測方程，因此一旦估計(jì)出了傾角場，將傾角參數(shù)代入預(yù)測方程（4），預(yù)測的殘差就是繞射波的近似。

圖1給出了平面波預(yù)測繞射波分離算法的流程圖。在平面波預(yù)測繞射波分離方法中，傾角場估計(jì)的正則化平滑參數(shù)是一個(gè)重要的參數(shù)，它控制著傾角的空間平滑性，從而約束著預(yù)測的范圍。同時(shí)，實(shí)際地震勘探中噪聲的存在也會影響傾角的估計(jì)和繞射波的分離，因此需要對噪聲和正則化參數(shù)進(jìn)行深入的分析討論。

圖1 繞射波分離算法流程圖Fig.1 A flowchart for separation of diffraction wave

2 噪音水平和光滑半徑影響因素

以地震勘探中被廣泛采用的Sigsbee 2A模型為例對繞射波分離算法的影響因素進(jìn)行分析。在繞射波分離的過程中，傾角場預(yù)測結(jié)果對繞射波分解效果影響較大，而噪聲和平面波預(yù)測平滑參數(shù)均對傾角場的估計(jì)產(chǎn)生影響。圖2?4中，在目標(biāo)數(shù)據(jù)中增添了不同方差的均值為0的正態(tài)分布隨機(jī)噪音，給出了不同噪聲水平下的傾角場估計(jì)結(jié)果和繞射波分離結(jié)果。隨著噪音水平的增大，原始數(shù)據(jù)共偏移距地震數(shù)據(jù)變得復(fù)雜，估計(jì)的傾角場準(zhǔn)確性降低，尤其在低信噪比情況下傾角場估計(jì)的誤差變大，但是高信噪情況下傾角場估計(jì)的結(jié)果較好，盡管在原始數(shù)據(jù)中可以觀察到噪聲，但估計(jì)的傾角幾乎和沒有噪聲數(shù)據(jù)估計(jì)的傾角一致，其主要原因是正則化平滑技術(shù)約束了傾角場的平滑性，提高估計(jì)的穩(wěn)定性和可靠性。在高信噪比情況下，盡管傾角場估計(jì)比較準(zhǔn)確，但是繞射波的分離仍然受到了噪聲影響，其原因?yàn)樵谄矫娌A(yù)測中我們只把可預(yù)測的平面波提取出來，而噪聲也被留在了殘差中，而殘差正是繞射波的估計(jì)，因此在分離的繞射波中含有噪聲。噪聲水平越強(qiáng)，分離的繞射波效果越差，在低信噪比條件下，傾角場的估計(jì)也變差，繞射波分離結(jié)果幾乎被淹沒在噪聲和假象下，無法得到較好的繞射波信息。由于噪聲會降低傾角場估計(jì)的精度，同時(shí)噪聲也會影響分離算法的效果，因此，在利用平面波預(yù)測法進(jìn)行繞射波分離時(shí)，必須要先壓制隨機(jī)噪音，從而提高分離結(jié)果的精度。

圖2 不含噪聲地震數(shù)據(jù)的傾角場估計(jì)和繞射波分離a.原始含有繞射波的數(shù)據(jù)，b.估計(jì)的傾角場，c.分離的繞射波。Fig.2 The dip field estimation and diffraction separation of seismic data without noisea.The original data containing diffraction wave,b.estimated dip field,c.separated diffraction wave.

圖3 高信噪比地震數(shù)據(jù)的傾角場估計(jì)和繞射波分離（噪聲方差為0.001）a.原始含有繞射波的數(shù)據(jù)，b.估計(jì)的傾角場，c.分離的繞射波。Fig.3 The dip field estimation and diffraction separation of high SNR seismic data（noise variance is 0.001）a.The original data containing diffraction wave,b.estimated dip field,c.separated diffraction wave.

傾角估計(jì)問題屬于反問題范疇，存在多解性和不適定性，為了解決這一問題，需要采用正則化技術(shù)對估計(jì)的局部傾角場進(jìn)行約束。在反問題正則化中光滑半徑參數(shù)控制著正則化程度，下面討論在傾角估計(jì)過程中，沿著水平方向平滑和沿垂直方向平滑半徑大小對傾角預(yù)測和繞射波分離結(jié)果的影響。圖5顯示了不同平滑半徑下傾角估計(jì)和繞射波分離的結(jié)果。隨著水平或垂直方向光滑半徑的不斷增大，傾角估計(jì)值的范圍越來越小，傾角估計(jì)的結(jié)果在水平或垂直方向上越來越平滑，隨著光滑半徑的進(jìn)一步增大，理論上當(dāng)光滑半徑趨于無窮大時(shí)，傾角估計(jì)結(jié)果也會趨于與坐標(biāo)軸平行。從對應(yīng)的繞射波分離結(jié)果來看，隨著光滑半徑的增大，傾角估計(jì)值變小，繞射波分離結(jié)果顯得越來越干凈，信噪比明顯提高。具體原因有兩方面：一是大的光滑半徑將部分噪音平滑掉了，減少了噪音對分離結(jié)果的干擾；二是大的光滑半徑將小的繞射波忽略了，這樣大尺度的繞射波就顯得很干凈。光滑半徑太小會使傾角估計(jì)的范圍大于理論值，使繞射波分離結(jié)果中出現(xiàn)許多全波場中沒有的小尺度繞射波，使得繞射波分離結(jié)果信噪比降低，利用這些小尺度繞射波偏移成像后可能出現(xiàn)許多小尺度地質(zhì)構(gòu)造體假象，即造成偏移噪音和偏移假象[19]。

圖4 低信噪比地震數(shù)據(jù)的傾角估計(jì)和繞射波分離（噪聲方差為0.025）a.原始含有繞射波的數(shù)據(jù)，b.估計(jì)的傾角場，c.分離的繞射波。Fig.4 The dip estimation and diffraction separation of low SNR seismic data（noise variance is 0.025）a.The original data containing diffraction wave,b.estimated dip field,c.separated diffraction wave.

因此，在應(yīng)用平面波預(yù)測進(jìn)行傾角估計(jì)和繞射波分離時(shí)應(yīng)選擇合適的光滑半徑。光滑半徑較小時(shí)，算法對小型地質(zhì)體的繞射波識別更敏感。光滑半徑較大時(shí)，小尺度繞射波沒有很好地分離出來，不利于小尺度地質(zhì)體的識別與成像。但較大的光滑半徑對大型的繞射波聚焦頂點(diǎn)的識別非常清晰。在處理實(shí)際數(shù)據(jù)時(shí)，應(yīng)根據(jù)具體的研究目的對平滑參數(shù)做適當(dāng)?shù)恼{(diào)整。例如，識別的目標(biāo)是地下小尺度構(gòu)造體、小型縫洞時(shí)，可以選擇較小的光滑半徑；若研究的目標(biāo)是大型裂縫、斷層、不整合面以及地層缺失時(shí)，應(yīng)當(dāng)選擇較大的光滑半徑。

3 地震數(shù)據(jù)繞射波分離算例

為了說明平面波預(yù)測地震繞射波分離的效果，我們將其應(yīng)用于一個(gè)模擬數(shù)據(jù)和兩個(gè)實(shí)際數(shù)據(jù)繞射波分離中。圖6a是一個(gè)模擬的共偏移距道集，可以清楚地看到繞射波和反射波互相干擾。圖6b是預(yù)測出的傾角場，這里的光滑半徑為5。圖6c是經(jīng)過分離后的結(jié)果，可以看出分離的繞射波中反射波同相軸幾乎不存在，繞射波信息明顯，該算例表明平面波預(yù)測濾波器可以很好地預(yù)測出反射波信息并對其進(jìn)行壓制。

然后，我們將其應(yīng)用于實(shí)際海洋地震勘探的共偏移距疊前道集（圖7a），可以發(fā)現(xiàn)該數(shù)據(jù)繞射波非常發(fā)育，從疊前共偏移距道集上可以清楚看到一些繞射波，尤其是對于3.5 s處復(fù)雜斷塊引起的繞射波非常清楚，我們首先估計(jì)其局部傾角場（圖7b），光滑半徑確定為5。在反射波發(fā)育的地方，傾角場比較平滑且數(shù)值較小，而在繞射波發(fā)育的地方傾角較大，從分離的繞射波結(jié)果（圖7c）可以看出平面波預(yù)測方法可以很好地分離繞射波，這些被分離出的繞射波可以很好地進(jìn)行斷塊構(gòu)造的成像[20]及儲層表征[21]。

圖5 不同平滑半徑對傾角估計(jì)及繞射分離結(jié)果的影響a-b是光滑半徑為3的結(jié)果，c-d是光滑半徑為10的結(jié)果，e-f是光滑半徑為50的結(jié)果。Fig.5 Influence of different smoothing radius on dip angle estimation and diffraction separation results.a-b is the result of smooth radius of 3,c-d is the result of smooth radius of 10,e-f is the result of smooth radius of 50.

最后，我們對疊后的實(shí)際地震數(shù)據(jù)進(jìn)行繞射波分離（圖8a），由于疊加會破壞一些繞射波的信息，因此在疊后的剖面上繞射波能量較弱，但是疊后的剖面信噪比較高，對估計(jì)傾角場比較有利。圖8b為估計(jì)的局部傾角場信息，由于疊加對繞射波的壓制，同時(shí)因?yàn)樵摂?shù)據(jù)剖面上的繞射波分布散亂，為了得到正確的傾角，此處的橫向光滑半徑為50，縱向光滑半徑為20，圖8c顯示了疊后地震數(shù)據(jù)繞射波分離后的結(jié)果，可以發(fā)現(xiàn)繞射波很好地被分離出來。模擬和實(shí)際數(shù)據(jù)表明平面波預(yù)測方法可以有效分離繞射波。

圖6 模擬的具有繞射波的共偏移距道集a.原始數(shù)據(jù)，b.預(yù)測傾角場，c.分離的繞射波。Fig.6 The synthetic common offset gather with diffracted wavesa.original data,b.prediction of dip field,c.separated diffracted waves.

圖7 實(shí)際海洋地震數(shù)據(jù)的共偏移距道集a.原始數(shù)據(jù)，b.估計(jì)的傾角場，c.分離的繞射波。Fig.7 The common offset gather of real marine seismic dataa.original data,b.estimated dip field,c.separated diffracted waves.

4 結(jié)論

平面波預(yù)測濾波器在反射波同相軸光滑連續(xù)的條件下，可以有效地估計(jì)地震數(shù)據(jù)的局部傾角場，進(jìn)而分離出繞射波，該方法對隱藏在強(qiáng)反射能量下的小型繞射波信息起到很好的分離效果，為后期繞射波單獨(dú)成像提供繞射波信息。通過對噪聲和平滑參數(shù)對地震繞射波分離結(jié)果的影響分析，發(fā)現(xiàn)噪聲水平嚴(yán)重影響繞射波分離效果，噪音過大分離結(jié)果會不準(zhǔn)確；平面波預(yù)測濾波器光滑半徑偏小會引入噪音和假象，光滑半徑偏大會導(dǎo)致繞射波分離不徹底，小斷塊地震繞射波響應(yīng)難以分離。在繞射波分離算法中，振幅保真性也是一個(gè)重要問題，因此提高算法的抗噪性和保真性是今后研究的重要方向。

圖8 實(shí)際地震數(shù)據(jù)疊加剖面a.原始數(shù)據(jù)，b.估計(jì)的局部傾角，c.分離的繞射波。Fig.8 The field poststack seismic sectiona.original data,b.estimated local dip field,c.separated diffraction wave.