韓建光， 張曉波， 邢占濤， 陳　鵬， 王　赟， 于常青

(1.中國(guó)地質(zhì)科學(xué)院 地質(zhì)研究所，北京　100037；2.中國(guó)地質(zhì)科學(xué)院，北京　100037；3.中國(guó)華電集團(tuán)科學(xué)技術(shù)研究總院有限公司，北京　100160；4.成都理工大學(xué)　油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，成都　610059；5.中國(guó)地質(zhì)大學(xué)　地球物理與信息技術(shù)學(xué)院，北京　100083)

?

高斯束疊前深度偏移影響因素分析

韓建光1， 張曉波2， 邢占濤3， 陳鵬4， 王赟5， 于常青1

(1.中國(guó)地質(zhì)科學(xué)院 地質(zhì)研究所，北京100037；2.中國(guó)地質(zhì)科學(xué)院，北京100037；3.中國(guó)華電集團(tuán)科學(xué)技術(shù)研究總院有限公司，北京100160；4.成都理工大學(xué)油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，成都610059；5.中國(guó)地質(zhì)大學(xué)地球物理與信息技術(shù)學(xué)院，北京100083)

摘要：高斯束偏移方法是一種優(yōu)秀的偏移算法，不僅具有接近于波動(dòng)方程偏移方法的成像精度，而且保留了Kirchhoff偏移方法高效、靈活的優(yōu)點(diǎn)。高斯束偏移成像效果以及計(jì)算效率受許多因素影響，這里以二維共炮域高斯束疊前深度偏移方法為基礎(chǔ)，分析了初始束寬、成像角度及速度光滑程度對(duì)偏移效果和效率的影響機(jī)理，并通過(guò)洼陷模型、Marmousi模型以及Sigsbee 2B模型的偏移試算，對(duì)初始束寬、成像角度以及速度光滑程度對(duì)高斯束疊前深度偏移的影響進(jìn)行了分析。

關(guān)鍵詞：高斯束； 疊前深度偏移； 影響因素； 數(shù)值模型分析

0引言

高斯束方法將波動(dòng)方程與射線理論相結(jié)合，同時(shí)考慮了波的運(yùn)動(dòng)學(xué)特征和動(dòng)力學(xué)特征，通過(guò)運(yùn)動(dòng)學(xué)射線追蹤獲取射線路徑，動(dòng)力學(xué)追蹤獲取中心射線附近的能量分布。高斯束偏移方法是近年來(lái)發(fā)展的一種偏移算法，它不但克服了Kirchhoff 偏移方法無(wú)法解決的多值走時(shí)問(wèn)題，同時(shí)保留了Kirchhoff偏移方法高效、靈活的優(yōu)點(diǎn)，具有接近于波動(dòng)方程偏移的成像精度[1]。

在地震學(xué)領(lǐng)域，利用高斯束方法進(jìn)行波場(chǎng)計(jì)算[2-4]最早應(yīng)用于地震波場(chǎng)的正演模擬[5-6]，之后一些國(guó)外學(xué)者對(duì)高斯束偏移進(jìn)行了研究，Hill[7-8]先后提出了高斯束疊后偏移方法以及基于共偏移距道集的疊前深度偏移方法；Hale[9-10]詳細(xì)介紹了高斯束偏移與Kirchhoff偏移以及傾斜疊加相比的優(yōu)點(diǎn)，并給出了具體的計(jì)算公式以及優(yōu)化算法，分析了高斯束偏移的計(jì)算效率；Gray[11-12]針對(duì)Hill方法對(duì)觀測(cè)系統(tǒng)適應(yīng)性不足，提出了共炮域的疊前偏移方法，并在以后的研究中基于單程波真振幅偏移理論，提出了真振幅高斯束偏移方法；Popov[13-14]提出了一種新的真振幅高斯束疊前深度偏移理論，對(duì)地下復(fù)雜地質(zhì)構(gòu)造獲得了高質(zhì)量的成像結(jié)果。

高斯束偏移是一種優(yōu)秀的偏移算法，但實(shí)現(xiàn)過(guò)程比較復(fù)雜，需要一系列的數(shù)值計(jì)算方法以及優(yōu)化算法[10]。此外高斯束偏移成像受許多因素影響，如初始束寬的大小、成像角度的控制以及速度模型的光滑程度等。這里從成像效果以及計(jì)算效率兩個(gè)方面，通過(guò)不同的數(shù)值模型試算對(duì)其主要影響因素進(jìn)行研究分析。

1高斯束方法原理

1.1基本原理

高斯束方法將波動(dòng)方程與射線理論相結(jié)合，在射線中心坐標(biāo)系下求解波動(dòng)方程。二維各向同性介質(zhì)中，中心射線坐標(biāo)系下高斯束頻率域的表達(dá)式為[2- 3]

1.2高斯束特點(diǎn)

高斯束是波動(dòng)方程在射線坐標(biāo)系下的高頻近似解，具有波動(dòng)方程的相關(guān)特性。與傳統(tǒng)的射線方法不同，高斯束方法不僅在中心射線上有波場(chǎng)能量，而且可以計(jì)算射線周圍的波場(chǎng)值，其形態(tài)受初始束寬和頻率的影響。

圖1為不同初始束寬高斯束在均勻介質(zhì)中傳播的瞬時(shí)波場(chǎng)圖，介質(zhì)中的波速為3 000 m/s，子波頻率為25 Hz，初始束寬分別為150 m、260 m、340 m和500 m。從圖1中可以看出，初始束寬的選取影響著波前面的形態(tài)。初始寬度越小，高斯束波前曲率變化劇烈，波場(chǎng)擴(kuò)散越快；初始寬度越大，高斯束波前曲率變化較為平緩，波場(chǎng)擴(kuò)散越慢。

圖1　初始束寬不同的高斯束在均勻介質(zhì)中傳播的瞬時(shí)波場(chǎng)圖Fig．1　Gaussian beam wavefield in a homogeneous medium with different initial beam width(a)150 m ；(b)260 m ；(c)340 m ；(d)500 m

圖2為不同頻率高斯束在均勻介質(zhì)中傳播的瞬時(shí)波場(chǎng)圖，高斯束的初始寬度相同，均為260 m，頻率分別為20 Hz、25 Hz、30 Hz和40 Hz。從圖2中可以看到，高斯束的波前曲率與頻率無(wú)關(guān)，頻率的大小影響著高斯束在中心射線周圍能量的分布。

結(jié)果表明，頻率越高高斯束能量越集中于中心射線附近，即頻率越高沿垂直于中心射線方向衰減越快，高斯束的寬度越窄。

圖2　不同頻率的高斯束在均勻介質(zhì)中傳播的瞬時(shí)波場(chǎng)圖Fig．2　Gaussian beam wavefield in a homogeneous medium with different frequencies(a)20 Hz ；(b)25 Hz ；(c)30 Hz ；(d)40 Hz

2高斯束疊前深度偏移原理

高斯束疊前偏移成像是由震源處正向延拓的波場(chǎng)與束中心點(diǎn)處反向延拓波場(chǎng)的互相關(guān)獲得。在二維各向同性介質(zhì)中，假設(shè)xs=(xs,0)和xr=(xr,0)分別為震源和接收點(diǎn)。根據(jù)Hill[8]提出的共偏移距域高斯束疊前偏移公式，給出了二維共炮域高斯束疊前偏移公式為式(2)。

(3)

式中：ωr為參考頻率；L0表示高斯束的初始束寬；u(xr,xs,ω)為地震記錄頻譜。

用A、T分別表示高斯束的復(fù)振幅值和旅行時(shí)，則高斯束表達(dá)式可以表示為式(4)。

uGB(x,x0,p,ω)=Aexp(iωT)

(4)

將高斯束表達(dá)式(4)代入式(2)，則高斯束疊前深度偏移公式可以表示為式(5)。

(5)

(6)

3高斯束偏移影響因素分析

3.1初始束寬的影響

在進(jìn)行高斯束偏移時(shí)，需要選取一系列的參數(shù)，如初始束寬L0、束中心間隔α、射線參數(shù)采樣間隔Δp等。作為波動(dòng)方程的高頻漸進(jìn)解，高斯束在傳播路徑上寬度越窄，其精度越高。但是如果初始寬度過(guò)小，高斯束在傳播過(guò)程波場(chǎng)會(huì)迅速擴(kuò)散，如圖1所示。將每個(gè)束中心點(diǎn)附近的道集通過(guò)局部?jī)A斜疊加分解為不同出射方向的局部平面波進(jìn)行波場(chǎng)延拓，是高斯束偏移方法的核心算法之一，能夠有效地減少計(jì)算量。在保證波場(chǎng)充分覆蓋的前提下，適當(dāng)?shù)卦龃笫行拈g隔，可以提高偏移的計(jì)算效率。高斯束疊前偏移需要在震源和束中心點(diǎn)處分別以不同的射線參數(shù)出射高斯束進(jìn)行波場(chǎng)計(jì)算，為了保證偏移成像的準(zhǔn)確性，初始射線參數(shù)的間隔必須足夠小以保證具有充足的射線覆蓋。相比于束中心間隔，初始射線參數(shù)間隔的大小對(duì)高斯束偏移的計(jì)算時(shí)間有著顯著的影響。因此，如何選取上述參數(shù)對(duì)高斯束偏移具有重要的影響。一般情況下可以參照Hill[7-8]所給定的選取準(zhǔn)則

L0=2πVa/ωr

(7)

(8)

(9)

式中：Va為偏移速度的幾何平均值；ωh為最高有效頻率。初始束寬的大小決定束中心間隔和射線參數(shù)采樣間隔的取值。

采用洼陷模型模擬數(shù)據(jù)對(duì)初始束寬，對(duì)高斯束

偏移的影響進(jìn)行了測(cè)試，洼陷速度模型如圖3所示。正演模擬采用有限差分方法，震源為30 Hz的Ricker子波，共81炮，炮間距為50 m，每炮301檢波點(diǎn)雙邊接收，道間距為10 m，采樣時(shí)間為2.6 s，采樣間隔為2 ms。圖4為初始束寬分別為100 m、200 m、300 m、500 m的洼陷模型偏移結(jié)果，從圖4中可以看到，不同的初始束寬高斯束偏移成像結(jié)果基本相同。雖然初始束寬的選取對(duì)單一高斯束周圍波場(chǎng)計(jì)算產(chǎn)生較大的影響，但高斯束偏移過(guò)程中通過(guò)高斯束疊加成像均可以得到準(zhǔn)確的偏移結(jié)果，說(shuō)明高斯束偏移方法對(duì)初始束寬大小不敏感。

圖3　洼陷速度模型Fig．3　Sag velocity model

圖4　不同初始寬度的洼陷模型高斯束疊前深度偏移結(jié)果Fig．4　Gaussian beam prestack depth migration results with different initial beam width(a)初始束寬為 100 m的成像結(jié)果；(b)初始束寬為 200 m的成像結(jié)果；(c)初始束寬為300 m的成像結(jié)果；(d)初始束寬為500 m的成像結(jié)果

通過(guò)模型試驗(yàn)可知，高斯束偏移的成像效果對(duì)初始寬度的大小并不是特別的敏感，只要給定的初始束寬值與式(7)計(jì)算得到的值相差不大，基本都可以得到較為理想的成像結(jié)果。但是由于初始束寬的變化導(dǎo)致束中心間隔和射線參數(shù)采樣間隔不同，對(duì)計(jì)算效率產(chǎn)生一定的影響。因此，可以在保證偏移精度不受影響的前提下，適當(dāng)?shù)恼{(diào)整上面的參數(shù)，以提高偏移計(jì)算效率。

3.2成像角度的影響

除了上述參數(shù)之外，成像角度的控制也是影響高斯束偏移的重要因素。成像角度過(guò)大可能會(huì)導(dǎo)致不必要的偏移噪音，而且需要花費(fèi)更多的計(jì)算時(shí)間,成像角度過(guò)小會(huì)限制某些陡傾構(gòu)造的成像。

作者通過(guò)對(duì)Marmousi模型的試算，分析了不同成像角度對(duì)高斯束偏移成像的影響，Marmousi速度模型如圖5所示。正演記錄共240 炮，炮間距為25 m，每炮96 道，道間距為25 m，最小偏移距為200 m，時(shí)間采樣間隔為4 ms，采樣點(diǎn)數(shù)為750。圖6(a)和圖6(b)分別為最大成像角度為100°和150°時(shí)的高斯束疊前深度偏移成像結(jié)果。從圖6可以看出，對(duì)于最大成像角度為100°時(shí)所得到的偏移結(jié)果，由于受成像角度限制，斷層附近的陡傾構(gòu)造未能成像(矩形所示部分)。而當(dāng)最大成像角度為150°時(shí)，斷層附近陡傾構(gòu)造獲得較好的成像，高斯束偏移準(zhǔn)確地恢復(fù)了模型復(fù)雜的構(gòu)造形態(tài)。通過(guò)Marmousi模型測(cè)試可知，成像角度的控制對(duì)于高斯束疊前深度偏移具有較大的影響，可以根據(jù)偏移的需要，適當(dāng)?shù)卣{(diào)整成像角度的大小。對(duì)于存在陡傾構(gòu)造的模型，可以適當(dāng)增大成像角度，以保證陡傾構(gòu)造能夠獲得成像。

圖5　Marmousi速度模型Fig．5　Marmousi velocity model

圖6　不同成像角度的Marmousi模型高斯束疊前深度偏移結(jié)果Fig．6　Gaussian beam prestack depth migration results with different imaging angle(a)最大成像角度為100°的成像結(jié)果；(b)最大成像角度為150°的成像結(jié)果

3.3速度光滑的影響

高斯束疊前深度偏移方法的關(guān)鍵，在于運(yùn)動(dòng)學(xué)射線追蹤和動(dòng)力學(xué)射線追蹤。射線理論是地震波場(chǎng)的一種高頻近似，其應(yīng)用前提條件需滿足地震波速度在一個(gè)波長(zhǎng)范圍內(nèi)的相對(duì)變化很小，而實(shí)際的速度模型可能具有不連續(xù)的速度突變點(diǎn)、間斷面或很強(qiáng)的速度梯度。為了滿足射線理論的應(yīng)用前提，需要對(duì)速度模型進(jìn)行光滑處理以減小速度的相對(duì)變化，以保證射線追蹤的可行性和有效性。因此，速度模型的光滑處理，對(duì)基于射線理論的高斯束疊前深度偏移具有較大影響。

通過(guò)對(duì)Sigsbee 2B模型試算分析速度光滑處理對(duì)高斯束偏移成像效果的影響。Sigsbee 2B層速度模型如圖7(a)所示，模型中包含鹽丘、斷層以及繞射體等復(fù)雜構(gòu)造，圖7(b)為Sigsbee 2B偏移速度模型。圖8(a)和圖8(b)分別為采用原始速度以及速度光滑模型偏移得到的成像結(jié)果。從圖8中可以看到，采用原始速度進(jìn)行高斯束偏移，由于鹽丘構(gòu)造內(nèi)部速度與周圍速度差異較大，射線很難穿過(guò)高速鹽丘體，導(dǎo)致鹽丘體下部邊界未能成像。而速度模型進(jìn)行光滑處理之后再進(jìn)行偏移成像，鹽丘構(gòu)造下部界面取得了較好的聚焦成像，模型整體的成像效果有了明顯的提高。通過(guò)Sigsbee 2B模型測(cè)試可知，對(duì)速度模型進(jìn)行適當(dāng)?shù)墓饣幚?，可以明顯的改善高斯束偏移成像效果。

圖7　Sigsbee 2B速度模型Fig．7　Sigsbee 2B velocity model(a)層速度模型；(b)偏移速度模型

圖8　Sigsbee 2B模型高斯束疊前深度偏移結(jié)果Fig．8　Gaussian beam prestack depth migration results of Sigsbee 2B model(a)原始速度模型高斯束偏移成像結(jié)果；(b)速度光滑模型高斯束偏移成像結(jié)果

4結(jié)論

高斯束偏移成像效果以及計(jì)算效率受許多因素影響，作者分別通過(guò)洼陷模型、Marmousi模型以及Sigsbee 2B模型的偏移試算，系統(tǒng)地分析了初始束寬、成像角度以及速度光滑對(duì)高斯束疊前深度偏移的影響。雖然初始束寬的大小對(duì)高斯束偏移成像質(zhì)量影響不大，但初始束寬的大小決定束中心間隔和射線參數(shù)采樣間隔的取值，因此對(duì)計(jì)算效率產(chǎn)生較大的影響。成像角度對(duì)高斯束疊前深度偏移具有較大的影響，成像角度過(guò)大可能會(huì)導(dǎo)致不必要的偏移噪音，而且需要花費(fèi)更多的計(jì)算時(shí)間；成像角度過(guò)小會(huì)導(dǎo)致某些陡傾構(gòu)造無(wú)法成像。因此，需要綜合考慮高斯束偏移的成像效果和計(jì)算效率，適當(dāng)?shù)卣{(diào)整上面的參數(shù)。此外，速度光滑處理可以減小速度模型的相對(duì)變化，保證射線追蹤的可行性和有效性，對(duì)高斯束疊前深度偏移質(zhì)量有較大的影響，因此在高斯束偏移過(guò)程中應(yīng)對(duì)速度模型進(jìn)行適當(dāng)?shù)墓饣幚怼?/p>

參考文獻(xiàn)：

[1]岳玉波，李振春，錢忠平，等．復(fù)雜地表?xiàng)l件下保幅高斯束偏移[J]．地球物理學(xué)報(bào)，2012，55(4)：1376-1383.

YUE Y B,LI Z C,QIAN Z P,et al.Amplitude-preserved Gaussian beam migration under complex topographic conditions[J].Chinese Journal of Geophysics,2012,55(4):1376-1383.(In Chinese)

[4]Popov M M.A new method of computation of wave fields using Gaussian beams[J].Wave Motion,1982(4):85-97.

[6]NOWACK R L.Calculation of Synthetic Seismograms with Gaussian Beams[J].Pure and Applied Geophysics,2003,160:487-507.

[7]HILL N R.Gaussian beam migration[J].Geophysics,1990,55:1416-1428.

[8]HILL N R.Prestack Gaussian-beam depth migration[J].Geophysics,2001,66:1240-1250.

[9]HALE D.Migration by the Kirchhoff,slant stack and Gaussian beam methods[R].Colorado:CWP Annual Project Review Meeting,1992.

[10]HALE D. Computational aspects of Gaussian beam migration[R].Colorado:CWP Annual Project Review Meeting,1992.

[11]GRAY S H.Gaussian beam migration of common-shot records[J].Geophysics, 2005,70(4):S71-S77.

[12]GRAY S H,BLEISTEIN N.True-amplitude Gaussian-beam migration[J].Geophysics,2009,74(2):S11-S23.

[13]POPOV M M,SEMTCHENOK N M,POPOV P M,et al.Reverse time migration with Gaussian beams and velocity analysis applications[J].Extended Abstracts of 70thAnnual Internat EAGE Mtg,2008：F048.

[14]POPOV M M,SEMTCHENOK N M,POPOV P M,et al.Depth migration by the Gaussian beam summation method[J].Geophysics, 2010,75(2):S81-S93.

收稿日期：2015-04-09改回日期：2015-08-03

基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(41425017，41374131，41374128)；科技部863課題(2013AA064201)

作者簡(jiǎn)介：韓建光(1987-)，男，博士，主要從事地震偏移成像與多分量地震技術(shù)研究,E-mail:hanjianguang613@163.com。

文章編號(hào)：1001-1749(2016)03-0403-07

中圖分類號(hào)：P 631.4

文獻(xiàn)標(biāo)志碼：A

DOI：10.3969/j.issn.1001-1749.2016.03.17

The analysis on influence factors of imaging precision in Gaussian beam prestack depth migration

HAN Jian-guang1, ZHANG Xiao-bo2, XING Zhan-tao3,CHEN Peng4, WANG Yun5, YU Chang-qing1

(1.Institute of Geology,Chinese Academy of Geological Sciences,Beijing100037,China;2.Chinese Academy of Geological Sciences,Beijing100037,China;3.China Huadian Science and Technology Institute,Beijing100160,China;4.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation Chengdu University of Technology,Chengdu610059,China；5.School of Geophysics and Information Technology,China University of Geoscience,Beijing100083,China)

Abstract:Gaussian beam migration is an elegant and efficient depth migration method with accuracy comparable to wave-equation migration and efficiency and flexibility comparable to Kirchhoff migration. Gaussian beam migration imaging quality and computation efficiency is influenced by many factors. In this paper, based on Gaussian beam prestack depth migration algorithm of common-shot records in two-dimensional media, the influence mechanism of the initial beam width, imaging angle and velocity smooth on effectiveness and efficiency of migration are studied. Finally, the effects of the initial beam width, imaging angle and velocity smooth on Gaussian beam prestack depth migration are studied by testing of sub-sag model, Marmousi model and Sigsbee 2B model.

Key words:Gaussian beam; prestack depth migration; influence factor; numerical model analysis