謝 飛，李 佩，黃中玉，3，魏修成，3，朱成宏

（1.中國石油化工股份有限公司石油勘探開發(fā)研究院，北京 100083；2.中國地質(zhì)大學(xué)（北京），北京 100083；3.中國石油化工集團(tuán)公司多波地震重點實驗室，北京 100083）

在復(fù)雜介質(zhì)條件下，疊前深度偏移技術(shù)能提供很好的構(gòu)造成像結(jié)果。Kirchhoff積分疊前深度偏移方法以其靈活的目標(biāo)處理功能，高效的數(shù)值計算和對空間非規(guī)則采樣數(shù)據(jù)的適應(yīng)能力等特點，在石油工業(yè)界得到了廣泛應(yīng)用。然而，對于由復(fù)雜構(gòu)造引起的多值走時、陰影、焦散等問題，采用Kirchhoff偏移方法成像效果往往不如單程波方法。高斯束偏移方法是Kirchhoff偏移方法的改進(jìn)，既具備Kirchhoff偏移方法的靈活性和適應(yīng)性，又能解決焦散、多值走時等問題，其成像效果堪比單程波算子，同時還能克服單程波算子的傾角限制，成像回轉(zhuǎn)波［1－6］。

Hill［7］提出了疊后高斯束偏移方法，通過局部的傾斜疊加將波場分解為高斯束。對每個高斯波束偏移成像、疊加得到最終的偏移結(jié)果。Hale［8－9］詳細(xì)分析了Kirchhoff偏移、傾斜疊加偏移和高斯束偏移的關(guān)系，指出高斯束偏移與Kirchhoff偏移、傾斜疊加偏移相比的優(yōu)點，并給出了具體的計算公式，分析了其計算效率。

Hill［10］給出了疊前高斯束偏移的方法。在共偏移距偏移中，通過對所有偏移距射線參數(shù)進(jìn)行掃描找到使總走時虛部最小的炮點射線參數(shù)和檢波點射線參數(shù)的組合。這種掃描可以在稀疏的偏移成像點網(wǎng)格上執(zhí)行，然后將值再內(nèi)插到較密集的成像網(wǎng)格點上，因此具有較高的計算效率和精度。

Gray［11］指出由Hill提出的共偏移距偏移方法對于較規(guī)則的觀測系統(tǒng)是合適的，但是對于其它觀測系統(tǒng)的記錄，或者是地表高程變化大、近地表速度變化劇烈時不利于成像。Gray提出炮集高斯束偏移更適合陸地起伏地表的地震數(shù)據(jù)。

1 高斯束偏移方法原理

假設(shè)地震波場滿足標(biāo)量波動方程：

則可將由r′點源激發(fā)傳播到r點的地震波場表示為瑞雷積分形式：

Hill［10］的研究表明，點源波場的格林函數(shù)可用高斯射線束的積分形式表示：

其中，uGB（r，r′，p′；ω）為在r′點激發(fā)、方向為p′的射線束。因此，當(dāng)將（2）式中的格林函數(shù)用（3）式的高斯束形式表示時，就可以得到利用高斯束表示的地震波場。

當(dāng)射線束的中心不在r′點而在鄰近的r0點時，由于在高斯束發(fā)射的初始位置波束是特定傳播方向的平面波波場，Hill［10］引入相移因子exp［－iωp′（r′－r0）］，則r′點到r點的格林函數(shù)表示為

在克?；舴蚍e分公式（2）中使用格林函數(shù)（4）時，應(yīng)該將公式（2）中的積分范圍劃分為小的空間積分窗，并且保證在計算格林函數(shù)（4）時，空間積分窗的中心r0點與積分窗內(nèi)的任意點r′相距不遠(yuǎn)，能保持初始高斯束的平面波性質(zhì)，為此，Hill［10］引入了如下高斯函數(shù)：

其中，ωl是參考頻率是射線束寬度；a是網(wǎng)格點最小間距。并得到了地震數(shù)據(jù)局部傾斜疊加的頻率域表達(dá)形式：

利用Claerbout［12］的相關(guān)成像原理，將炮點波場和檢波點波場均按照克?；舴蚍e分公式向下延拓，并將兩個波場進(jìn)行互相關(guān)可得到偏移成像公式：

將偏移成像公式（7）中的格林函數(shù)用（4）式的高斯束表示，并結(jié)合高斯函數(shù)（5）式，就可以得到以高斯束形式表達(dá)的偏移成像公式：

其中，C0是與高斯函數(shù)相關(guān)的系數(shù)。

Hale［8］對克?；舴?、傾斜疊加和高斯束3種偏移方法進(jìn)行了詳細(xì)比較，揭示了三者間的聯(lián)系與區(qū)別，認(rèn)為在時間域做高斯束偏移要比在頻率域效率高，并給出了單個高斯束對地下成像貢獻(xiàn)的時間域表達(dá)式：

其中，AR和AI分別對應(yīng)復(fù)值振幅的實部和虛部；τR和τI分別為復(fù)值走時的實部和虛部。復(fù)值振幅和復(fù)值走時均可通過高斯束射線追蹤得到。px為射線參數(shù)，代表傾斜疊加形成的平面波傳播方向；（τR，τI）為高斯窗內(nèi)地震數(shù)據(jù)的傾斜疊加結(jié)果；（τR，τI）為（τR，τI）的Hilbert變換。（τR，τI）可由下式得到：

其中，Bj（ω，px）為第j個高斯窗內(nèi)通過傾斜疊加得到的頻率域結(jié)果，其具體表達(dá)式為

式中：C為與高斯窗函數(shù)相關(guān)的振幅系數(shù)；Fj（ω，kx）為地震記錄fj（x，t）的二維Fourier變換。

關(guān)于高斯窗函數(shù)的空間間隔與波束射線參數(shù)采樣間隔的關(guān)系，Hill［7，10］和Hale［8］均進(jìn)行了詳細(xì)的討論，我們采用了Hill［7，10］給出的參數(shù)選取準(zhǔn)則。

2 基于高斯射線束的角道集

角度域共成像點道集在AVA 分析和偏移速度分析方面有著非常廣闊的應(yīng)用前景。郭朝斌等［13］和李振春等［14］給出了利用粗網(wǎng)格走時信息計算角度域共成像點道集的方法。我們利用高斯束中心射線所攜帶的豐富的速度、方向信息，便捷直觀地給出了角度域共成像點道集的計算方法。

在中心射線坐標(biāo)系下，某一網(wǎng)格點處的走時可由對應(yīng)的中心射線通過二階泰勒展開得到［15］（圖1）：

其中，τ0為中心射線走時；M可由動力學(xué)射線追蹤得到。

利用坐標(biāo)變換，將中心射線坐標(biāo)系轉(zhuǎn)換為直角坐標(biāo)系，（13）式變?yōu)?/p>

其中，

式中：θ表示射線方向與Z軸的夾角。由（14）式對走時分別計算X方向和Z方向的偏導(dǎo)數(shù)，得到慢度px和pz：

其中，

得到成像點處的水平和垂直慢度后，我們可以分別計算出來自炮點的入射波場的射線角度φs 和來自射線束中心點的反射波場的射線角度φr：

在分別計算出φs 和φr 后，可通過計算二者之間的夾角來得到孔徑角。對于常規(guī)縱波勘探而言，反射角γ為孔徑角的一半，即

圖1 由中心射線計算成像點走時和射線方向原理示意圖解

3 模型測試

我們首先利用Marmousi國際標(biāo)準(zhǔn)模型對基于高斯束偏移的角度域共成像點道集算法進(jìn)行了測試。

圖2給出了Mamousi速度模型。圖3給出了采用本文方法對Mamousi模型數(shù)據(jù)進(jìn)行偏移得到的結(jié)果，可見成像效果良好。在輸出時我們選擇的角度范圍為［－60°，60°］，但由于觀測系統(tǒng)為右邊放炮左邊接收，所以在［0，60°］基本沒有能量，因此，我們只顯示［－60°，0］的角度域道集（圖4）。圖4a對應(yīng)圖3中最左邊的三角形位置，該處構(gòu)造相對簡單，呈現(xiàn)層狀結(jié)構(gòu)，對應(yīng)的角道集也被拉平；圖4b對應(yīng)圖3中間的三角形位置，圖4c對應(yīng)圖3中最右邊的三角形位置，這兩處對應(yīng)的構(gòu)造復(fù)雜，橫向變速劇烈，角道集也同樣得到了拉平。由此可見，即使在如此復(fù)雜的構(gòu)造條件下，算子的成像精度依然很高，保證了最終高質(zhì)量的成像結(jié)果。

圖2 Marmousi模型

我們設(shè)計了一個包含直立斷層的復(fù)雜模型，該模型尺寸水平方向為27 000m，垂直方向為4 000m，其速度場如圖5所示。利用有限差分技術(shù)模擬炮集記錄，觀測系統(tǒng)為中間放炮兩邊接收，檢波點間隔為50m，炮點間隔為100m，偏移距范圍為［－6 000m，6 000m］，共110炮，第1炮炮點位置在8 000m處。

圖5 包含直立斷層的復(fù)雜構(gòu)造速度模型

圖6a和圖6b分別給出了采用二階廣義屏偏移方法和高斯束偏移方法對該復(fù)雜模型的炮集記錄進(jìn)行偏移得到的結(jié)果。對比圖6a和圖6b可以看出，兩種方法整體上均能對構(gòu)造有比較好的成像效果，但是在陡傾角區(qū)域，高斯束偏移方法成像效果更好（圖6中紅色箭頭、紅色方框、紅色橢圓所示位置），二階廣義屏偏移方法無法對直立斷層進(jìn)行成像，而高斯束偏移方法能夠清晰、準(zhǔn)確地對直立斷層進(jìn)行成像。

圖7給出了圖6中紅色方框內(nèi)偏移結(jié)果的局部放大顯示。對比圖7a和圖7b可見，廣義屏算法在陡傾角構(gòu)造部位引入了較強(qiáng)的偏移噪聲，而且沒有將能量歸位到正確的位置上；而高斯束偏移結(jié)果信噪比高，同時將能量歸位到了正確的傾角構(gòu)造位置上。但由于觀測系統(tǒng)的原因，高斯束方法也未成像出中間的一段陡傾角構(gòu)造。圖8給出了圖6中紅色橢圓內(nèi)偏移結(jié)果的局部放大顯示。對比圖8a和圖8b可見，對于陡傾角構(gòu)造，高斯束的偏移結(jié)果要比廣義屏算法能量聚焦更好。

同時，我們對上述復(fù)雜構(gòu)造模型采用二階廣義屏偏移和高斯束偏移的計算效率進(jìn)行了對比，其結(jié)果如表1所示?？梢?，高斯束偏移的計算效率要明顯高于二階廣義屏偏移的計算效率。

表1 高斯束偏移與二階廣義屏偏移有關(guān)參數(shù)

4 應(yīng)用效果分析

我們將高斯束疊前深度偏移算法應(yīng)用于西亞某地區(qū)的實際資料成像處理。研究區(qū)地質(zhì)條件復(fù)雜，速度橫向變化劇烈，存在高陡構(gòu)造，而且原始資料品質(zhì)差，經(jīng)過前期處理的偏移輸入道集如圖9a所示。由圖9a可見，該區(qū)地震資料信噪比較低，尤其是在中深層，有效信號被噪聲淹沒。圖9b為建立的深度域速度模型。圖10a和圖10b分別給出了對研究區(qū)地震資料進(jìn)行克?；舴蚱坪透咚故频慕Y(jié)果；圖11a和圖11b分別給出了CMP50和CMP250處高斯束偏移生成的角度域共成像點道集。對比圖10a和圖10b可見，高斯束偏移結(jié)果信噪比較高，同相軸連續(xù)性更好，特別是在陡傾角區(qū)域（圖中橢圓區(qū)域），高斯束偏移方法實現(xiàn)了對高陡構(gòu)造的優(yōu)質(zhì)成像；在深層由于輸入道集信噪比低，克?；舴蚍椒乖肽芰ο鄬^弱，成像效果較差，而高斯束偏移剖面中深層信噪比較高，同相軸連續(xù)性較好。

圖11 角度域共成像點道集

5 結(jié)束語

高斯束偏移方法保留了射線方法靈活高效、無傾角限制的優(yōu)勢，同時克服了焦散、陰影區(qū)、多值走時等問題；基于射線中心坐標(biāo)系的動力學(xué)射線追蹤得到的慢度信息可直接用于計算共成像點角度道集。復(fù)雜模型和實際數(shù)據(jù)的測試結(jié)果驗證了高斯束偏移方法在陡傾角構(gòu)造、低信噪比資料成像方面的優(yōu)勢。

本研究僅限于各向同性介質(zhì)的構(gòu)造成像，在振幅保持以及各向異性高斯束偏移方面有待更深入的研究。

［1］Virgilio M，Simone R，Daniele C.Gaussian beam migration：prestack，common－shot，TTI，true amplitude in angular domain［J］.Expanded Abstracts of 78thAnnual Internat SEG Mtg，2008，2236－2239

［2］Cockshott I.Specular beam migration—a low cost 3D prestack depth migration［J］.Expanded Abstracts of 76thAnnual Internat SEG Mtg，2006，2634－2637

［3］Gao F，Zhang P，Dirks V.Fast beam migration—a step toward interactive imageing［J］.Expanded Abstracts of 75thAnnual Internat SEG Mtg，2005，2470－2474

［4］ˇCervenyV，Popov M M，Pˇsenˇcík I.Computation of wave fields in inhomogeneous media—gaussian beam approach［J］.Geophysical Journal Royal Astronomical Society，1982，70：109－128

［5］ˇCervenyV，Pˇsenˇcík I.Gaussian beams in elastic 2－D laterally varying layered structures［J］.Geophysical Journal Royal Astronomical Society，1984，78：65－91

［6］Costa C A，Raz S，Kosloff D.Gaussian beam migration［J］.Expanded Abstracts of 59thAnnual Internat SEG Mtg，1989，1169－1171

［7］Hill N R.Gaussian beam migration［J］.Geophysics，1990，55（11）：1416－1428

［8］Hale D.Migration by the Kirchhoff，slant stack，and Gaussian beam methods［R］.Colorado：CWP Annual Project Review Meeting，1992

［9］Hale D.Computational aspects of Gaussian beam migration［R］.Colorado：CWP Annual Project Review Meeting，1992

［10］Hill N R.Prestack Gaussian beam depth migration［J］.Geophysics，2001，66（4）：1240－1250

［11］Gray S H.Gaussian beam migration of common－shot records［J］.Geophysics，2005，70（4）：S71－S77

［12］Claerbout J F.Imaging the earth's interior［M］.Oxford：Blackwell Scientific Publication，1985：376－378

［13］郭朝斌，李振春，岳玉波.高斯束成像技術(shù)及其應(yīng)用［J］.石油物探，2011，50（1）：38－44，58 Guo C B，Li Z C，Yue Y B.Gaussian beam migration and its application［J］.Geophysical Prospecting for Petroleum，2011，50（1）：38－44，58

［14］李振春，岳玉波，郭朝斌，等.高斯波束共角度保幅深度偏移［J］.石油地球物理勘探，2010，45（3）：360－365 Li Z C，Yue Y B，Guo C B，et al.Gaussian beam common angle preserved－amplitude migration［J］.Oil Geophysical Prospecting，2010，45（3）：360－365

［15］ˇCerveny V.Seismic ray theory［M］.New York：Cambridge University Press，2001：234－400