地震波逆時(shí)偏移方法研究綜述_圖文

1、第33卷第3期2010年6月文章編號(hào) :1671-8585(2010 03-0153-07勘探地球物理進(jìn)展P rogr ess in Explor ation Geophysics Vol. 33, No. 3Jun. , 2010地震波逆時(shí)偏移方法研究綜述陳可洋(中國石油大慶油田有限責(zé)任公司勘探開發(fā)研究院, 黑龍江大慶163712摘要:地震波逆時(shí)偏移是目前精度較高的深度域成像方法, 不同學(xué)者從不同的角度對(duì)逆時(shí)偏移方法進(jìn)行了研究。因此, 在回顧地震波逆時(shí)偏移方法發(fā)展歷程的基礎(chǔ)上, 概述了當(dāng)前流行的幾種疊前、疊后逆時(shí)成像方法及其適用性, 通過模型實(shí)例分析了幾種疊后逆時(shí)成像方法存在的一些問題, 同

2、時(shí)對(duì)比了幾種疊前逆時(shí)成像方法對(duì)Mar m 2ousi 模型的逆時(shí)成像效果, 最后指出了逆時(shí)偏移方法的發(fā)展方向。關(guān)鍵詞:逆時(shí)偏移; 波動(dòng)方程; 爆炸反射面; 半速替換; 初至旅行時(shí); 相關(guān)法; 逆時(shí)成像中圖分類號(hào):P63114文獻(xiàn)標(biāo)識(shí)碼:A地震波逆時(shí)偏移法是現(xiàn)行偏移方法中最精確的一種成像方法, 是一種基于波動(dòng)理論的深度域偏移方法, 它與地震波正演數(shù)值模擬問題剛好相反, 將檢波器接收到的波場(chǎng)進(jìn)行逆時(shí)延拓, 利用逆時(shí)成像條件實(shí)現(xiàn)對(duì)地下各點(diǎn)的成像。由于逆時(shí)偏移方法采用全波場(chǎng)波動(dòng)方程, 并對(duì)波動(dòng)方程中的微分項(xiàng)進(jìn)行差分離散(如有限差分法、有限元法、偽譜法等 , 對(duì)波動(dòng)方程的近似較少(這取決于差分方程與微分

3、方程的逼近程度 , 因此不受傾角和偏移孔徑的限制, 可以有效地處理縱橫向存在劇烈變化的地球介質(zhì)物性特征(如速度、密度等 。Whitmore 1(1983 最早在第53屆SEG 年會(huì)上提出逆時(shí)偏移方法, Baysal 等(1983 提出了不同的逆時(shí)偏移概念, Loewenthal 和Mufti 3(1983 將逆時(shí)偏移方法應(yīng)用在空間-頻率域。Levin 等(1984 概括了逆時(shí)偏移的基本原理和實(shí)現(xiàn)方法。H ildebrand 等5(1987 將逆時(shí)偏移方法應(yīng)用于波阻抗成像, 取得了很好的效果。Levy 和Esmersoy642三維逆時(shí)成像。Yoon 等15(2004 指出逆時(shí)偏移方法需要解決的問

4、題, 如逆時(shí)偏移方法可以消除多次波效應(yīng), 但是無法壓制逆時(shí)偏移過程中的次生干擾。Zhang 等16(2007 提出真振幅逆時(shí)偏移的計(jì)算策略, Sava 等17(2005 提出了時(shí)移成像方法, 薛東川等(2008 采用有限元法實(shí)現(xiàn)了彈性波動(dòng)方程疊前逆時(shí)偏移, Guan 等19(2008 采用多步法高效實(shí)現(xiàn)了逆時(shí)偏移, Liu 等20(2008 提出了反假頻波動(dòng)方程正演和逆時(shí)偏移方法, Wards 等2118(2008 提出了相移時(shí)間步逆時(shí)偏移方法, Soubaras 等22(2008 提出了兩步法顯式匹配逆時(shí)偏移方法, Jones 23(2008 提出了逆時(shí)偏移的預(yù)處理方法, Guo 等(200

5、8 采用三角網(wǎng)格有限差分法實(shí)現(xiàn)了疊后逆時(shí)偏移。Chattopadhyay 等25(2008 和Costa 等26(2009 提出了時(shí)空移相關(guān)法疊前逆時(shí)成像方法, 實(shí)現(xiàn)了復(fù)雜介質(zhì)疊前逆時(shí)成像。陳可洋(2009 提出基于旅行時(shí)成像條件的高階交錯(cuò)網(wǎng)格有限差分疊前逆時(shí)偏移方法, Zhang 等31(2009 采用一步內(nèi)插法實(shí)現(xiàn)逆時(shí)深度偏移, Dus 2saud 等(2008 研究了逆時(shí)偏移的計(jì)算策略問題等。本文對(duì)逆時(shí)偏移方法進(jìn)行了數(shù)值計(jì)算和歸納總結(jié)。32273024(1988 研究了逆時(shí)波場(chǎng)外推、成像和反演問題。Zhu 和Lines 7(1994 比較了逆時(shí)偏移與克希霍夫積分偏移方法, 得出前者對(duì)Ma

6、r mousi 模型成像精度更高, 但耗時(shí)多的結(jié)論。Wu 等(1996 研究了三維高階有限差分法逆時(shí)偏移技術(shù), Causse 和U rsin 9(2000 進(jìn)行了粘彈性波動(dòng)方程試算, 證明逆時(shí)偏移方法對(duì)粘彈性波一樣適用。Sun 和McMechan (2001 進(jìn)行了標(biāo)量波動(dòng)方程的逆時(shí)深度偏移, 對(duì)縱波和橫波的成像結(jié)果表明, 彈性波逆時(shí)偏移比單純聲波效果更好。張美根等(2001 實(shí)現(xiàn)了各向異性彈性波疊前、疊后逆時(shí)深度偏移, 張會(huì)星等12(2002 實(shí)現(xiàn)了彈性波動(dòng)方程疊前逆時(shí)偏移, Mul 2der 等(2003 比較了單程和雙程波動(dòng)方程偏移, 1413111081 逆時(shí)成像方法文中采用的地震波正

7、演遞推算子和逆時(shí)延拓算子的具體計(jì)算公式可參見文獻(xiàn)27, 在此不再贅述。下面分述4種逆時(shí)成像條件及其適用性。收稿日期:2009-12-28; 改回日期:2010-02-04。作者簡(jiǎn)介:陳可洋(1983 , 男, 助理工程師, 碩士, 主要從事地震資154勘探地球物理進(jìn)展第33卷1. 1 基于爆炸反射界面原理和零時(shí)間的成像方法爆炸反射界面成像方法是最常用、最簡(jiǎn)單的一種成像方法, 適合于單程旅行時(shí)零偏移距剖面的疊后偏移成像, 不適合于疊前逆時(shí)成像。它把地下反射界面假設(shè)為爆炸源, 爆炸源的形狀、位置與反射界面的形狀和位置完全一致, 它所產(chǎn)生的脈沖波的強(qiáng)度和極性與界面反射系數(shù)的大小和正負(fù)一致, 所以稱地

8、下反射界面為爆炸反射界面。在零時(shí)刻所有爆炸反射界面同時(shí)起爆, 在地表接收地震波場(chǎng), 接著將地表接收到的波場(chǎng)以給定的模型速度沿逆時(shí)間方向延拓至零時(shí)刻, 此時(shí)的波場(chǎng)值就能夠正確描述地下反射界面的位置, 自動(dòng)實(shí)現(xiàn)疊后逆時(shí)深度成像。1. 2 半速替換原理和零時(shí)間的成像方法該方法同樣基于爆炸反射界面的基本假設(shè), 不同的是, 它考慮的是雙程旅行時(shí)零偏移距剖面偏移成像, 不考慮界面反射系數(shù)的大小與極性, 同樣不適用于疊前逆時(shí)成像。該方法將震源置于整個(gè)地表, 在零時(shí)刻所有震源同時(shí)起爆, 沿著速度模型向地下傳播, 遇地質(zhì)界面后反射至地表, 并被所有地表檢波器所接收。接著, 檢波器接收到的波場(chǎng)以給定介質(zhì)速度的一半

9、向地下傳播(這里采用半速替換的目的是使雙程旅行時(shí)等效為單程旅行時(shí), 因?yàn)檫@里只考慮反射界面至地表檢波器的距離, 它是地表 地下界面 地表總距離的一半(基于水平介質(zhì)假設(shè) 。速度減半后, 雙程旅行時(shí)就自動(dòng)等效于單程旅行時(shí) , 并沿著逆時(shí)間方向延拓至零時(shí)刻, 此時(shí)的波場(chǎng)已經(jīng)能夠正確地反映出地下反射界面的位置, 從而實(shí)現(xiàn)了疊后逆時(shí)成像。1. 3 初至旅行時(shí)成像該方法(可參見陳可洋提出的基于旅行時(shí)成像條件的疊前逆時(shí)成像方法 不同于以上兩種方法, 它通常適用于疊前聲波和彈性波逆時(shí)偏移成像, 本文將其推廣應(yīng)用于聲波疊后雙程旅行時(shí)零偏移距逆時(shí)偏移成像。它采用波場(chǎng)傳播的時(shí)間一致性成像原理, 即地震波到達(dá)地下某一

10、反射界面的時(shí)間與地表接收波場(chǎng)逆時(shí)延拓至該反射界面位置的旅行時(shí)一致, 否則該位置不在反射界面上, 由此篩選出符合條件的成像值, 組合后得到逆時(shí)偏移剖面。其具體過程可描述為:將震源置于地表(對(duì)疊后逆時(shí)偏移情況, 將震源置于整個(gè)地表, 而對(duì)于疊前逆時(shí)偏移情況, 則只置于點(diǎn)震源位置 , 在零時(shí)刻所有震源同時(shí)起爆, 沿著速度模型向地下傳播, 遇地質(zhì)界面后反射至地表被所有地表檢波器接收, 接著27件, 并將檢波器接收到的波場(chǎng)以給定介質(zhì)速度沿逆時(shí)間方向向地下延拓, 在延拓的每一個(gè)時(shí)間步中, 根據(jù)已知的初至旅行時(shí)成像條件, 篩選出具有相同旅行時(shí)的波場(chǎng)值作為其成像值, 從而實(shí)現(xiàn)逆時(shí)成像。1. 4 相干成像該方法

11、不同于以上3種方法, 它適合于處理具有雙程旅行時(shí)特征的疊前、疊后聲波和彈性波逆時(shí)偏移成像。它采用波場(chǎng)傳播的最大相干性成像原理, 即地震波到達(dá)地下某一反射界面的波場(chǎng)與地表接收到的波場(chǎng)逆時(shí)延拓至該反射界面的波場(chǎng)具有最大相干性。如果兩者的相干性不是最大, 那么該位置不在反射界面上, 由此記錄最大相干值并累加作為該位置處的成像值, 從而實(shí)現(xiàn)逆時(shí)偏移。其具體過程與初至旅行時(shí)成像方法相似, 差別在于不必求取初至旅行時(shí)成像條件, 而是對(duì)地表接收波場(chǎng)與震源波場(chǎng)在同一時(shí)刻作相干成像處理, 篩選出某一空間位置處的最大相干累加值作為其成像值, 從而實(shí)現(xiàn)逆時(shí)成像。目前典型的幾種相干成像方法如下(以二維為例 。1. 4

12、. 1 非時(shí)空移相關(guān)法這是最為常用的成像方法, 它是將震源波場(chǎng)(S 與接收波場(chǎng)(R 在同一空間位置(x, z 作互相關(guān)運(yùn)算再累加來估算該位置處的成像值I (x, z :I (x, z =E E Sst sXs(x, z , t *R s (x, z, t (1其頻率域表達(dá)式為:I (x, z =S (x, z, X R (x, z, X E Ess(2式中:x 和z 分別表示水平和深度軸, t 為時(shí)間, 下標(biāo)s 代表炮集個(gè)數(shù)。該成像方法對(duì)小阻抗差結(jié)構(gòu)的成像是有效的, 如圖1所示。式(2 與常規(guī)單程波偏移成像方法相比, 兩者的計(jì)算公式相同, 這很好地說明了單程波偏移方法和雙程波逆時(shí)偏移方法具有相

13、似的理論基礎(chǔ)。1. 4. 2 非時(shí)空移反射系數(shù)相關(guān)法該方法借鑒了照明分析的基本思想, 并可進(jìn)一步壓制次生干擾, 其計(jì)算公式如下:震源照明I 1(x, z =E S (x, z, t *R (x, z, t EE S (x, z, tssts2s t(3接收照明I 2(x, z =E S (x, z, t *R (x, z, t EE R (x, z, tssts2s t(4第3期陳可洋1地震波逆時(shí)偏移方法研究綜述155式(3 和式(4 在頻域也與常規(guī)單程波偏移算法計(jì)算公式相同。為了避免分母自相關(guān)值接近于零, 導(dǎo)致數(shù)值溢出(發(fā)散 , 對(duì)式(3 和式(4 的分母進(jìn)行修正處理, 增加一修正因子E (

14、也可稱為穩(wěn)定性因子 , 于是得到修正后的計(jì)算公式如下:震源照明I 1(x, z =接收照明I 2(x, z =E S (x, z, t *R (x, z, t EE R (x, z, t +Essts2s t(6式(5 和式(6 即為平滑逆時(shí)成像方程。Guit 2ton 等33(2007 指出, 在理論上, 修正因子E 取決于地震資料的信噪比, 而在實(shí)際計(jì)算中, 可通過選(5取不同的修正因子E 值作交互試驗(yàn), 分析其剩余譜 特征來實(shí)現(xiàn)參數(shù)的合理選取。EsE S s (x, z , t *R s (x, z, ttE St2s(x , z, t +E圖1 震源波場(chǎng)(a 、接收波場(chǎng)(b 和震源波場(chǎng)

15、與接收波場(chǎng)的相關(guān)處理(c1. 4. 3 時(shí)空移相關(guān)法時(shí)空移成像方法是通過震源波場(chǎng)與接收波場(chǎng)在空間方向和時(shí)間方向均作互相關(guān)運(yùn)算來估算最佳成像值的一種方法。它由Sava 和Fomel 17(2006 提出, 是1. 4. 1方法的一般情況, 其計(jì)算公式如下:I (m =子波作為震源, 聲波速度從上到下依次為2000m/s 、2500m/s 、3000m/s 、3500m/s, 時(shí)間步長為015ms, 縱橫向空間步長均為5m 。采用一階雙曲型聲波方程, 應(yīng)用高階交錯(cuò)網(wǎng)格有限差分法進(jìn)行離散, 差分精度為時(shí)間2階和空間10階, 在邊界處采用完全匹配層吸收邊界條件, 制作出不同零偏移距剖面, 并實(shí)現(xiàn)相應(yīng)的

16、疊后逆時(shí)深度偏移。(7EsE S s (m +h, t +S *tR s (m-h, t -S 在頻域可表示為 I (m =Es2i X SS s (m+h, X R s (m-h, X e E X(8式中:i=-1, m 代表空間坐標(biāo)位置(x, z , h 代表空間方向相關(guān)計(jì)算位移量(h x , h z , S 代表時(shí)間方向相關(guān)計(jì)算位移量。由該方程同樣可以推導(dǎo)出時(shí)空移反射系數(shù)相關(guān)法逆時(shí)成像方程。Sava 等17(2006 還指出時(shí)移比空移更為有效, 可以大大減小計(jì)算量和存儲(chǔ)量。2 實(shí)例分析2. 1 疊后逆時(shí)深度偏移以二維帶凹陷的層狀介質(zhì)模型為例, 如圖2所 10001000圖2 層狀介質(zhì)速度

17、模型圖3為基于爆炸反射界面原理和零時(shí)間成像156勘探地球物理進(jìn)展第33卷深度偏移剖面, 采用的主頻為30H z, 總旅行時(shí)間為1s; 圖4和圖5均為基于半速替換和零時(shí)間成像方法得到的零偏移距剖面和疊后逆時(shí)深度偏移剖面, 總旅行時(shí)間為2s, 采用的主頻分別為40H z 和10H z; 圖6是采用主頻40H z 、基于自激自收方式合成雙程旅行時(shí)零偏移距剖面, 再根據(jù)旅行時(shí)成像方法得到的疊后逆時(shí)深度偏移剖面。對(duì)比圖3a 、圖4a 、圖5a 可知, 盡管3種方法采用的主頻不同, 但剖面中均無明顯的數(shù)值頻散和邊界反射現(xiàn)象, 說明合成零偏移距剖面的計(jì)算精度較高。對(duì)比圖3b 、圖4b 、圖5b 和圖6可知,

18、 采用爆炸反射界面原理和零時(shí)間成像方法的疊后逆時(shí)偏移剖面能夠準(zhǔn)確地恢復(fù)地質(zhì)界面, 且信噪比較高。采用半速替換和零時(shí)間成像方法的疊后逆時(shí)偏移剖面也可以準(zhǔn)確地恢復(fù)水平層位, 但剖面中存在較強(qiáng)的頻散現(xiàn)象(即偽波動(dòng) , 其隨著主頻的降低而逐步減弱, 并且繞射波沒有收斂。究其原因, 首先是由于檢波器接收到的反射波是雙程旅行時(shí), 而繞射波傳播的路徑是從繞射點(diǎn)至地表檢波器, 因此是單程旅行時(shí), 采用半速替換后, 繞射波變成了/半單程旅行時(shí)0, 無法實(shí)現(xiàn)收斂。其次, 根據(jù)K =和N =(這里f $hK 為波長, v 為波速, f 為頻率, $h 為空間步長, N 為每波長所占有的空間節(jié)點(diǎn)數(shù) , 在給定頻率f

19、的情況下進(jìn)行半速替換, 波長K 必然減小, 而此時(shí)空間步長$h 大小不變, 那么每個(gè)波長所占空間樣點(diǎn)數(shù)必然減少一半。通常每個(gè)波長占8個(gè)采樣點(diǎn)時(shí)才能保證整個(gè)計(jì)算過程具有較弱的頻散現(xiàn)象(只有精細(xì)網(wǎng)格剖分才能最大程度壓制頻散, 而這必然帶來巨大的計(jì)算量 , 而這里半速替換后, 每個(gè)波長所占采樣點(diǎn)數(shù)減小為4個(gè), 從而導(dǎo)致了頻散現(xiàn)象顯著。由此可見, 該方法只適合于低頻和水平介質(zhì)的疊后逆時(shí)偏移。這也從側(cè)面說明, 該方法適合于CDP 域的疊前逆時(shí)成像, 因?yàn)?CDP 第3期陳可洋1地震波逆時(shí)偏移方法研究綜述157圖5 主頻10H z 時(shí), 基于自激自收方式合成的零偏移距剖面(a 及其半速替換后的疊后逆時(shí)深度

20、偏移剖面(b域波場(chǎng)的偏移處理是基于水平層狀介質(zhì)假設(shè)的。為了解決反射波和繞射波所具有的單、雙程旅行時(shí)問題, 本文采用垂直旅行時(shí)成像方法, 這樣可以保證繞射波準(zhǔn)確地收斂于繞射點(diǎn), 反射波準(zhǔn)確歸位(圖6 。由于不需要進(jìn)行半速替換, 因此偏移剖面無頻散現(xiàn)象。對(duì)比圖4b 、圖5b 和圖6也不難發(fā)現(xiàn), 剖面中均存在一些次生干擾, 這是基于雙程波 動(dòng)方程方法進(jìn)行逆時(shí)偏移所難免的。另外還需指出, 基于旅行時(shí)成像條件的疊后逆時(shí)偏移剖面中, 斷層會(huì)削弱部分深層同相軸的能量。2. 2 疊前逆時(shí)深度偏移采用3種逆時(shí)成像方法分別對(duì)Marmousi 速度模型(圖7 進(jìn)行疊前逆時(shí)成像試驗(yàn)。模型總大小為1915m 605m,

21、 空間縱橫向網(wǎng)格間距均為5m, 最小速度為1500m/s, 最大速度為5500m/s, 密度為1. 0g/cm 3, 時(shí)間步長為014ms, 滿足穩(wěn)定性條件。采用最大頻率為80H z 的零相位雷克子波作為震源, 采用右側(cè)放炮多次覆蓋的觀測(cè)系統(tǒng), 震源和檢波器均置于地表, 第1炮位于模型的420m 處, 炮間距為10m, 道間距為5m, 排列長度為420m, 共激發(fā)了151炮。圖8a圖8c 分別為應(yīng)用非時(shí)空移相關(guān)成像方法、非時(shí)空移反射系數(shù)相關(guān)成像方法和基于旅行時(shí)成像方法的疊前逆時(shí)深度偏移剖面。對(duì)比圖8與圖7可以看出, 各種地質(zhì)構(gòu)造形態(tài)成像清晰準(zhǔn)確, 3條右傾主干大斷層、向斜和背斜構(gòu)造及不整合界面

22、清晰可見, 局部的小斷塊也清晰可辨, 說明這幾種成像方法的逆時(shí)偏移噪聲能量較小, 計(jì)算精 度較高。圖6 根據(jù)零偏移距剖面(圖4a 并基于旅行時(shí)成像條件生成的疊后逆時(shí)深度剖面158 勘探地球物理進(jìn)展 第 33 卷 圖8 疊前逆時(shí)深度偏移剖面 a 非時(shí)空移相關(guān)成像方法; b 非時(shí)空移反射系數(shù)相關(guān)成像方法; c 基于旅行時(shí)成像條件 和存儲(chǔ)量均是非常龐大的, 特別是三維逆時(shí)偏移情 3 結(jié)束語 地震波逆時(shí)偏移方法是當(dāng)前勘探或計(jì)算地球 況, 這主要是由相關(guān)成像運(yùn)算造成的; 它的主頻 有限, 在保證逆時(shí)偏移無數(shù)值頻散的前提下, 要提 高逆時(shí)偏移剖面的主頻必然要減小空間網(wǎng)格步長, 這勢(shì)必增加更多額外的計(jì)算量和

23、存儲(chǔ)量。 綜上分析可知, 逆時(shí)偏移的工業(yè)化應(yīng)用必需建 立在快速、 高效的計(jì)算機(jī)水平的基礎(chǔ)上。 參 1 物理學(xué)界的研究熱點(diǎn), 該方法采用精確差分離散計(jì) 算方法來求解雙程波動(dòng)方程, 不需要進(jìn)行上下行波 分解, 能有效處理多值走時(shí)問題, 并可以通過簡(jiǎn)單 的直達(dá)波切除和窗口相關(guān)成像處理, 或通過應(yīng)用旅 行時(shí)成像方法拾取成像值來實(shí)現(xiàn)復(fù)雜介質(zhì)高精度 逆時(shí)偏移成像。 目前, 逆時(shí)偏移仍存在一些有待繼續(xù)研究的問 題, 如: 它是基于雙程地震波動(dòng)方程, 不可避免地 在逆時(shí)偏移過程中會(huì)產(chǎn)生次生干擾; 無論采用哪 種成像方法和數(shù)值離散計(jì)算方法, 它所需的計(jì)算量 考 文 獻(xiàn) Whitmore N D. Iterati

24、ve depth imaging by backwar d time pr opagation J . Expanded Abstracts of 53 rd An2 nual Inter national SEG Meeting, 1983, 382 384 2 3 Baysa l E, Kosloff D D, Sherwood J W C. Rever se time migration J . Geophysics, 1983, 48( 11 : 1 514 1 524 Loewenthal D, Mufti I R. Rever se 2time m igr ation in 第3期

25、 陳可洋1 地震波逆時(shí)偏移方法研究綜述 159 spatial frequency domain J . Geophysics, 1983, 48( 5 : 627 635 4 5 Levin S A. Pr inciples of r everse t ime migr at ion J . Ge2 ophysics, 1984, 49( 5 : 581 583 Hildebr and S T. R everse 2time dept h migration: Imped2 ance imaging condition J . Geophysics, 1987, 52 ( 8 : 21 20

26、 of 78t h Annual Internationa l SEG Meet ing, 2008, 2 341 2 345 Liu F Q, Zhang G Q, Mor ton S A, et a l. An ant i dis2 2 per sion wave equation for modeling and r everse 2time migration J . Expanded Abst racts of 78 t h Annual In2 ternational SEG Meeting, 2008, 2 277 2 281 War ds B D, Mar gr ave G F

27、 . Phase 2shift time 2stepping for rever se ime migration J . Expanded Abstracts of 2t 78 th Annual International SEG Meeting, 2008, 2 262 2 266 22 Soubar as R, Zhang Y. Two 2step explicit marching method for r everse time migration J . Expanded Ab2 st racts of 78 t h Annual Inter nat ional SEG Meet

28、ing, 2008, 2 272 2 276 Jones I F . P re 2processing consider at ions for r ever se time migrat ion J . Expanded Abstracts of 78 t h Annual International SEG Meeting, 2008, 2 297 2 301 24 Guo S J, Li Z C, Sun X D, et a l. Post2stack r everse 2time migration using a finite difference method based on t

29、ri2 angular gr ids J . Applied Geophysics, 2008, 5 ( 2 : 115 120 25 Chattopadhyay S, Mcmechan G A. Imaging conditions for prestack r everse 2time migrat ion J , Geophysics, 2008, 73( 3 : S81 S89 Costa J C, Silva Neto F A, Alcantara M R M, et a l. Ob2 liquity 2correction imaging condition for rever s

30、e time migration J . Geophysics, 2009, 74( 3 : S57 S66 27 28 陳可洋. 基于高階 有限差 分的波 動(dòng)方程 疊前逆 時(shí)偏移 方法 J . 石油物探, 2009, 48( 5 : 475 478 陳可洋, 楊微, 吳清嶺, 等. 幾種地震波疊后深度偏移方 法的比較 J . 勘 探地球 物理 進(jìn)展, 2009, 32( 4 : 257 260 陳可洋. 標(biāo)量聲波 波動(dòng)方 程高階 交錯(cuò)網(wǎng) 格有限 差分法 J . 中國海上油氣, 2009, 21( 4 : 232 236 陳可洋. 高階彈性 波動(dòng)方 程正演 模擬及 其逆時(shí) 偏移成 像研究 D

31、. 黑龍江大慶 : 大慶石油學(xué)院, 2009 Zhang Y, Zang G. One 2step extr apolation method for reverse time migr ation J . Geophysics, 2009, 74 ( 4 : A29 A33 32 Dussaud E, Symes W W, Williamson P, et al. Computa 2 tional strategies for r everse 2time migration J . Ex2 panded Abstracts of 78 t h Annual International SE

32、G Meeting, 2008, 2 267 2 271 33 Guitton A, Valenciano A, Bevc D, et al. Smoothing im 2 aging condition for shot2pr ole migr ation J . Geo2 physics, 2007, 72( 3 : S149 S154 1 060 1 064 6 Levy B C, Esmer soy C. Var iable background Bor n in2 version by wavefield back pr opagation J . Jour nal on Appli

33、ed Mathematics, 1988, 48( 4 : 30 32 7 8 Zhu J, Lines L. Imaging of complex subsur face str uc 2 tur es by VSP migration J . CJEG, 1994, 30( 1 : 73 83 Wu W, Lines L R, Lu H . Analysis of higher 2order, fi2 nite 2differ ence schemes in 32D reverse 2time migrat ion J . Geophysics, 1996, 61( 3 : 845 856

34、 Causse E, U rsin B. Viscoacoustic rever se ime migra2 2t tion J . Journal of Seismic Exploration, 2000, 9 ( 1 : 23 9 165 184 10 Sun R, M cMechan G A. Scalar rever se 2time depth mi2 gration of elast ic seismic data J . Geophysics, 2001, 66 ( 5 : 1 515 1 518 11 張美根, 王 妙月. 各 向異 性彈 性波 有限 元疊 前 逆時(shí) 偏 移 J

35、 . 地球物理學(xué)報(bào), 2001, 44( 5 : 711 719 12 張會(huì)星 , 寧書年. 彈性波動(dòng)方程疊前逆時(shí)偏移 J . 中國 礦業(yè)大學(xué)學(xué)報(bào), 2002, 31( 5 : 372 375 13 Mulder W, Plessix R. One 2way and two 2way wave 2e2 rd quation migration J . Expanded Abstr act s of 73 An2 nual International SEG Meeting, 2003, 1 292 1 295 Yoon K, Shin C, Suh S, et al. 3D reverse 2

36、time migra2 tion using acoustic wave equation: An experience with the SEG/ EAGE data set J . The Leading Edge, 2003, 15 22( 1 : 38 41 Yoon K, Marfurt K J, Houston U, et al. Challenges in rever se ime migration J . Expanded Abstr act s of 74t h 2t Annual Internationa l SEG Meeting, 2004, 1 057 1 060

37、16 Zhang Y, Sun J, Gray S. Reverse 2time migration: Am2 plitude and implementation issues J . Expanded Ab2 str act s of 77 th Annual Inter national SEG Meeting, 17 2007, 2 145 2 149 Sava P, Fomel S. T ime 2shift imaging condition J . Ex2 panded Abstracts of 75 t h Annual Inter nat ional SEG Meet ing

38、, 2005, 1 850 1 853 18 薛東川 , 王尚旭. 波動(dòng)方程有限元疊前逆時(shí)偏移 J . 石 油地球物理勘探, 2008, 43( 1 : 17 21 19 Guan H, Li Z, Wang B, et al. A multi step approach for 2 efficient reverse 2time migration J . Expanded Abstracts 26 14 29 30 31 ( 編輯: 戴春秋 第 33 卷第 3 期 2010 年 6 月 勘探地球物理進(jìn)展 P rogr ess in Explor ation Geophysics Vol.

39、33, No. 3 Jun. , 2010 ABSTRACT Review of seismic reverse time migr ation methods. Chen Keyang. PEG, 2010, 33( 3 : 153 159 Seismic rever se time migration is a method of high imaging precision in depth domain. Different realizations of reverse mi2 gration have been put forward. This paper review ed t

40、he history of seismic rever se time migration, and summarized several prev2 alent prestack and poststack reverse time migration schemes and their adaptability. We demonstrated different poststack re2 ver se time imaging conditions with model data, and compar ed the results fr om different reverse ti

41、me migration methods on Marmousi model. The outlook of r everse time migr ation was al2 so forecast. Keywords: reverse time migration; wave equation; exploding re2 flector inter face; half 2speed substitution; first arrival; correla2 tion method; reverse time imaging condition Chen K eyang, Explorat

42、ion and Development Research Institute of Daqing O ilfield Company Limited, Daqing 163712, China Review of quality control of seismic data processing in the course of integrated study of seismology & geology. Li Dawei, Yang Kai, Yang Shunjun. PEG, 2010, 33( 3 : 160 167 W ith more and more seismology

43、 & geology integrated work being implemented in geophysical exploration industry, the quality control of seismic data processing has been received great attention. In this paper, three typical QC methods used in industry: well2control Q C method based on VSP data, inver sion2 oriented QC which is ba

44、sed on rock physics modeling, and con2 sistency iented Q C based on a delicate processing to keep true 2or amplitude, fr equency, phase and waveform relatively, w er e re2 view ed in detail. Discussions on the future development of QC methods of seismic data processing were suggested. Keywords: inte

45、gr ated study of seismology & geology; quality control; w ell control seismic data processing Li Dawei, School of Ocean and Ear th Science, Tongji U niver si2 ty, Shanghai 200092, China Realization of time 2domain residual curvature migr ation velocity analysis in iCluster software platform. Zhang B

46、ing, Wang Hua2 zhong, Sun Chenglong, Liu Shaoyong. PEG, 2010, 33 ( 3 : 168 173 Residual cur vature analysis ( RCA based on common ima2 ging point gathers is one of the effective methods in impr oving the accuracy of migration velocity. By analyzing the relationship between the residual curvature of

47、Kirchhoff PSTM common imaging point gathers and migration velocity, we developed a PSTM velocity analysis module on iCluster seismic data pro2 cessing platform. We tested the module with model data and re2 al data. A ll the tests show that this module can adjust the RMS velocity field locally and qu

48、antificationally, and can wor k well with Kirchhoff prestack time migration. Keywords: r esidual curvatur e analysis; migr ation velocity analy2 sis; PSTM; iCluster Zhang Bing, School of Ocean and Ear th Science, Tongji Univer2 sity, Shanghai 200092, China Forward modeling based on anisotropic gas 2

49、bear ing sandstone model and characteristics analysis of AVO response. Zhou Hua2 ilai, Li Luming, Luo Shengxian, Wang Mingchun. PEG, 2010, 33 ( 3 : 174 178 In order to investigate the multi wave char acteristics of e2 2 lastic waves in anisotropic gas 2bearing sandstone r eser voirs to facilitate th

50、e inter pretation of multi wave seismic data, this pa 2 2 per adopts high2order staggered id finite difference method to 2gr simulate propagation of seismic waves in anisotropic gas 2bear 2 ing sandstone media based on anisotropic elastic wave equation. The effect of anisotropic coefficients on mult

51、i wavefield was an 2 2 alyzed, and the responses of the wavefield in isotropic and ani 2 sotropic gas2bear ing sandstone reservoir s were compared. Re 2 sults show that anisotropy significantly affects wavefield char 2 acteristics and AVO response, and the impact manifests different behaviors on P 2

52、waves and converted waves. The amplitude charac 2 teristics in seismic record is consistent with the AVO reponse, indi 2 cating that the combination of forward modeling and AV O anal 2 ysis can provide a solid base for identifying multi wavefield. 2 Keywords: multi wave and multi component; wavefiel

53、d char ac 2 2 2 teristics; AVO analysis; per fectly matched layer; forward modeling Zhou Huailai, State Key Laboratory of Oil and Gas Reser voir Geology and Exploitation and Key Lab of Ear th Explor ation & Information Techniques of Ministry of Education, Chengdu U2 niversity of Technology, Chengdu

54、610059, China Exper iment on ground resistance in high 2density electrical meas 2 ur ement. Zhang Lingyun, Liu Hongfu, Li Chengyou. PEG, 2010, 33( 3 : 179 183 A n experimental study of electrical prospecting with uni 2 form and non2uniform gr ound resistance on the same testing field w as introduced. The variability of appar ent resistance and the effect of ground resistance on sensiti