解建建,王緒本

(成都理工大學(xué) 地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,四川成都 610059)

基于單程波真振幅分步傅里葉疊前深度偏移方法

解建建,王緒本

(成都理工大學(xué) 地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,四川成都 610059)

這里將單程波真振幅方程與分步傅里葉算子(SSF)相結(jié)合,同時(shí)還結(jié)合了保幅算法和分步傅里葉算法的優(yōu)點(diǎn),因此該方法具有計(jì)算量小,占內(nèi)存少,能處理橫向變化的速度等優(yōu)點(diǎn)。并且克服了傅里葉有限差分方法偏移后的振幅都有很大的偏差的不足。與目前廣泛應(yīng)用的常規(guī)的分步有限差分疊前深度偏移相比,具有成像精度高,保持地震波動(dòng)力學(xué)特征等優(yōu)點(diǎn)。在M armousi模型上成功地進(jìn)行了真振幅分步傅里疊前深度偏移處理,取得了理想的成像效果。

單程波動(dòng)方程;真振幅;分步傅里葉算子;疊前深度偏移

0 前言

為了得到較好的成像效果,人們致力于研究基于波場(chǎng)延拓理論的波動(dòng)方程疊前深度偏移算法。根據(jù)采用的波場(chǎng)延拓算子的不同,近年來(lái)研制出的波動(dòng)方程疊前深度偏移算法有多種:運(yùn)用波數(shù)～頻率域的相移法;相移加插值PSPI[1]算法;分步傅氏SSF算法[2];頻率～空間域的有限差分FD算法等[3]。在成像過(guò)程中,相移法和有限差分算法各有優(yōu)缺點(diǎn)。分步傅立葉偏移方法在處理速度橫向變化的能力與計(jì)算效率等方面,較之“相移加插值”方法明顯存在優(yōu)勢(shì)。該方法的地震成像沒(méi)有傾角限制,但對(duì)復(fù)雜地質(zhì)體(如超覆、逆掩斷層、鹽丘等)成像,其內(nèi)部速度場(chǎng)橫向變化非常劇烈,成像的誤差非常大。嚴(yán)格地說(shuō),相移法只適用于垂向速度變化,不允許速度的橫向變化。由于分步傅立葉方法沒(méi)有考慮速度場(chǎng)二階以上的擾動(dòng),因此,它很難使復(fù)雜地質(zhì)體正確成像。有限差分算法能夠保持地震波的動(dòng)力學(xué)特征,允許速度的任意變化,但存在偏移角度的限制。那么,如何才能有一種高計(jì)算效率,且成像精度高的偏移算法呢?

最近幾年來(lái),人們都是在原有方法的基礎(chǔ)上,在提高計(jì)算效率和保幅方面做進(jìn)一步研究。作者在本文依賴傳統(tǒng)的SSF偏移算子高計(jì)算效率的優(yōu)勢(shì),提出了基于單程波方程的真振幅分步傅里葉偏移算法,這樣既提高了計(jì)算效率,又解決了SSF算子在強(qiáng)橫向介質(zhì)中不適應(yīng)的問(wèn)題,并通過(guò)M armousi模型試算來(lái)驗(yàn)證了此算法。

1 真振幅分步傅里葉偏移原理

1.1 單程波方程

張關(guān)泉[4]證明了光滑不均勻介質(zhì)中,全聲波方程近似為

并由此提出了下面的單程波分解方程:

為了實(shí)現(xiàn)真振幅共炮偏移,張宇[5]提出了單程波方程和邊界條件

和

1.2 真振幅偏移算子

以上行波方程(3)為例:

(1)以參考速度做相移法。在空間頻率域中方程(5)分二步實(shí)現(xiàn):

其中 kx、ky分別為x、y方向的空間波數(shù);v0(z)為隨深度變化的參考速度。

方程(6)的求解可以寫成:

式(8)為變量z的常微分方程,經(jīng)求解可得:

(2)薄鏡校正項(xiàng)。求解該常微分方程(9)得到:

為了得到正確的偏移振幅,在偏移之前必須進(jìn)行球面擴(kuò)散校正。式(10)和式(11)就是在SSF的基礎(chǔ)上,針對(duì)散射項(xiàng)所做的偏移算子,能夠正確地顯示出振幅的信息。

2 模型試算

為了驗(yàn)證本文中方法的保幅性,首先我們?cè)O(shè)計(jì)了一個(gè)簡(jiǎn)單的2D凹槽模型。圖1是用傳統(tǒng)SSF偏移成像剖面,可見(jiàn)能夠給出正確的成像位置,但是該方法對(duì)構(gòu)造成像而言精度比較低。

圖1 傳統(tǒng)SSF偏移成像Fig.1 Trad itional SSFm igration im aging

圖2為真振幅SSF偏移成像剖面,它不但給出了正確的成像位置,而且整張剖面的成像振幅比較均衡。該方法對(duì)于傾斜界面和最下面的平層的成像,比傳統(tǒng)SSF方法的成像效果好,精度更高,而且對(duì)構(gòu)造成像更加清晰,這說(shuō)明該方法能有效處理強(qiáng)間斷速度變化。

作者在文中用M armousi模型驗(yàn)證該方法,能夠有效地解決復(fù)雜地質(zhì)模型。M armousi模型構(gòu)造復(fù)雜,存在劇烈的橫向變速,因此M armousi模型是國(guó)際公認(rèn)的檢驗(yàn)疊前深度偏移方法的理想模型。圖3是M armousi模型速度場(chǎng),這里所采用的數(shù)據(jù)集的參數(shù)為:共有240個(gè)炮集,每炮96道,炮間距25m,道間距25m,采樣率4m s。下頁(yè)圖4是在準(zhǔn)確的速度模型下,用傳統(tǒng)的SSF疊前深度偏移成像的結(jié)果。下頁(yè)圖5是在準(zhǔn)確的速度模型下,得到的真振幅SSF疊前深度偏移成像的結(jié)果。從二者的對(duì)比中可以更加明顯地看出,真振幅SSF偏移成像的結(jié)果中模型淺層的三個(gè)大的斷層,以及中深層的背斜和低速體等構(gòu)造形態(tài),都能夠得到更好的體現(xiàn),刻畫出了波場(chǎng)隨速度場(chǎng)變化的細(xì)節(jié),整個(gè)偏移剖面具有質(zhì)量很高的成像效果。通過(guò)以上針對(duì)M armousi模型的偏移成像結(jié)果的對(duì)比,可以看出與傳統(tǒng)SSF偏移算子相比,真振幅SSF偏移算子的優(yōu)勢(shì)所在。

3 結(jié)論

真振幅SSF偏移是基于傳統(tǒng)的SSF的一種偏移方法,通過(guò)對(duì)凹槽模型和M armousi模型的多個(gè)炮點(diǎn)的偏移成像結(jié)果對(duì)比,真振幅SSF疊前深度偏移方法是實(shí)現(xiàn)較復(fù)雜地質(zhì)構(gòu)造成像的一種快速有效的方法,它解決了傳統(tǒng)SSF方法不能保證實(shí)現(xiàn)保幅偏移成像的缺點(diǎn),在提高計(jì)算效率的同時(shí),也加強(qiáng)了速度場(chǎng)橫向變化的適應(yīng)性。在降低計(jì)算量和提高成像質(zhì)量方面將會(huì)有很大的優(yōu)越性。

圖2 真振幅SSF偏移成像Fig.2 True amp litude SSFm igration im aging

圖3 M armousi模型速度場(chǎng)Fig.3 Velocity field ofM armousimodel

圖4 M armousi模型傳統(tǒng)的SSF偏移成像結(jié)果Fig.4 Traditional SSFm igration im aging ofM armousimodel

圖5 M armousi模型真振幅SSF偏移成像結(jié)果Fig.5 True amp litude SSFm igration im aging ofM armousimodel

[1]GAZDAG J.W ave equationm igration w ith the phaseshiftm ethod[J].Geophysics,1978,(43):1342.

[2]STOFFA P L,FORKEMA J T,DE LUNA FRE～RE.ETC.Sp ilt-step Fourierm igration[J].Geophsics,1990,55(2):410.

[3]R ISTOW D.RUH 1 T.Fourier finite-differencem igration[J].Geophysic,1994,59(12):1822.

[4]張關(guān)泉.波動(dòng)方程的上行波和下行波的耦合方程組[J].應(yīng)用數(shù)學(xué)學(xué)報(bào),1993,18(2):251.

[5]賀振華,王才經(jīng),李建朝,等.反射地震資料偏穆處理與反演方法[M].重慶:重慶大學(xué)出版社,1989.

[6]張宇.振幅保真的單程波方程偏移理論[J].地球物理學(xué)報(bào),2006,49(5):1410.

[7]李振春,朱緒峰,韓文功,等.真振幅偏移方法綜述[J].勘探地球物理進(jìn)展,2008,2(31):12.

[8]崔興福,張關(guān)泉,吳雅麗.三維非均勻介質(zhì)中保幅地震偏移算子研究[J].地球物理學(xué)報(bào),2004,47(3):509.

[9]賀振華,王才經(jīng),李建朝,等.反射地震資料偏移處理與反演方法[M].重慶:重慶大學(xué)出版社,1989.

[10]張宇,徐升,張關(guān)泉.等.真振幅全傾角單程波方程偏移方法[J].石油物探,2007,46(6):582.

[11]呂彬,王宇超,李斐,等.保幅型裂步傅里葉疊前深度偏移方法探討[J].巖性油氣藏,2007,19(3):101.

[12]董偉,徐兆濤,分步傅立葉波動(dòng)方程三維疊前深度偏移算法及應(yīng)用[J].油氣地質(zhì)與采收率,2008,(增刊):16.

P 631.4

A

1001—1749(2011)01—0036—04

2010-07-26 改回日期:2010-11-01

解建建(1985-),女,碩士,現(xiàn)主要從事波動(dòng)方程偏移成像方面的研究。