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

解建建,王緒本

(成都理工大學 地球探測與信息技術教育部重點實驗室,四川成都 610059)

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

解建建,王緒本

(成都理工大學 地球探測與信息技術教育部重點實驗室,四川成都 610059)

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

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

0 前言

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

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

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

1.1 單程波方程

張關泉[4]證明了光滑不均勻介質中,全聲波方程近似為

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

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

和

1.2 真振幅偏移算子

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

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

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

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

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

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

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

2 模型試算

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

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

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

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

3 結論

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

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

圖3 M armousi模型速度場Fig.3 Velocity field ofM armousimodel

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

圖5 M armousi模型真振幅SSF偏移成像結果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]張關泉.波動方程的上行波和下行波的耦合方程組[J].應用數學學報,1993,18(2):251.

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

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

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

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

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

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

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

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

P 631.4

A

1001—1749(2011)01—0036—04

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

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