武泗海 趙 虎 尹 成 賈 鵬

(①西南石油大學地球科學與技術(shù)學院，四川成都 610500； ②西南石油大學天然氣地質(zhì)四川省重點實驗室，四川成都 610500； ③川慶鉆探地球物理勘探公司，四川成都 610213)

·綜合研究·

廣角地震反射特征及反演研究

武泗海*①②趙 虎①②尹 成①②賈 鵬③

(①西南石油大學地球科學與技術(shù)學院，四川成都 610500； ②西南石油大學天然氣地質(zhì)四川省重點實驗室，四川成都 610500； ③川慶鉆探地球物理勘探公司，四川成都 610213)

武泗海，趙虎，尹成，賈鵬﹒廣角地震反射特征及反演研究.石油地球物理勘探，2017,52(5):1005-1015.

在小炮檢距情況下，地震反射波振幅隨著入射角的變化而變化，但相位基本保持不變。當入射角接近或大于臨界角時，地震反射波相位也隨之發(fā)生變化。振幅與相位都是波形信息的重要特征，因此，本文對廣角反射波特征進行了研究。首先以平層介質(zhì)為例，結(jié)合廣角地震記錄合成方法，研究了t-x域、τ-p域的反射波振幅/相位隨入射角變化的特征，并希望據(jù)此進行廣角的振幅/相位反演。但由于τ-p變換受采樣間隔等因素的影響，本文又利用逆時偏移成像提取共成像點道集的方法，獲取廣角反射信息，并將其應(yīng)用于彈性介質(zhì)參數(shù)的反演試驗。最后以平層介質(zhì)和局部Marmousi-2模型為例做了相關(guān)實驗，反演結(jié)果也驗證了廣角振幅與相位及聯(lián)合反演的有效性。

廣角反演τ-p變換 共成像點道集 逆時偏移

1 引言

在實際勘探工區(qū)中，如逆掩斷層區(qū)、海洋深水區(qū)、火山巖、碳酸鹽巖上覆層發(fā)育區(qū)，采集的資料經(jīng)常遇到超臨界角反射的問題。當入射角大于臨界角時，由于反射系數(shù)變?yōu)閺蛿?shù)產(chǎn)生附加相角，反射波相位產(chǎn)生變化，波形產(chǎn)生畸變。王志等[1]通過高斯射線正演模擬廣角反射的方法，討論了廣角反射特征及影響波形識別的主要因素； 孫成禹等[2]通過分析廣角反射數(shù)據(jù)的特征，研究廣角數(shù)據(jù)的校正方法及存在的問題； 孟憲軍等[3]通過研究地震反射系數(shù)相角變化特征，指出了臨界角前后的相角變化特征。

AVO/AVA技術(shù)根據(jù)振幅隨炮檢距變化的特征進行地下介質(zhì)的巖性信息識別[4,5]，但這僅僅考慮振幅大小及正負極性對反演的影響，實際上，反射波相位的變化也同樣重要。對于子波而言，振幅和相位對波形都具有非常重要的影響。此外，由于振幅在傳播過程中存在透射、波前擴散及反射臨界角等問題，振幅信息的提取精度也影響了反演效果，而相位的影響較小，其精度相對較高。因此，進一步研究相位變化特征對廣角反射資料的應(yīng)用具有一定實際意義。

在廣角反射存在的情況下，由于Zoeppritz方程是基于平面波入射的假設(shè)，這就導致平面波反射系數(shù)(PRCs)與實際地震反射特征并不完全一致。針對這一問題，Haase[6]和Ursenbach等[7]研究了在大炮檢距情況下的球面波AVO響應(yīng)，并提出了一種簡潔的平緩界面球面波反射系數(shù)(SRCs)的計算方法；Ayzenberg等[8]提出了有效反射系數(shù)(ERCs)的概念，將其應(yīng)用于彎曲界面的AVO分析。另外，前人就球面波場的平面波分解問題進行了研究，使Zoeppritz方程在廣角范圍依然適用，van der Baan等[9]研究了球面波的τ-p域分析方法，使平面波反射系數(shù)繼續(xù)適用于振幅隨慢度變化的反演。Zhu等[10]研究了廣角反射在τ-p域的響應(yīng)，通過反射率法模擬結(jié)果證實τ-p域反射特征與平面波反射系數(shù)基本一致。但τ-p變換也存在不足之處，為了避免假頻，對空間采樣提出了嚴格的要求。

逆時偏移(RTM)作為目前成像精度最高的偏移技術(shù)，在提高復雜介質(zhì)的成像質(zhì)量方面表現(xiàn)出優(yōu)越的性能，而且RTM輸出的共成像點道集還能提供振幅和相位信息，為后續(xù)屬性解釋提供依據(jù)[11-13]。目前，RTM技術(shù)已從二維到三維，從聲波方程逐步發(fā)展到彈性波方程。許多學者針對RTM中的成像條件、邊界條件和成像噪聲壓制方面的問題做了大量的研究工作。薛東川[14]對比了兩大類成像條件：入射波到達時間類與入射波/反射波的相關(guān)類成像條件的特點、計算量及存儲量； 張智等[15]提出了彈性波逆時偏移的穩(wěn)定激發(fā)振幅成像條件，與歸一化成像條件相比，該成像條件具有更好的成像能力和空間分辨率； 陳康等[16]提出了頻率域修正的四階拉普拉斯濾波方法，該方法在改善成像噪聲壓制效果的同時，還保持了子波的振幅與相位。此外，為解決RTM波場存儲和計算量大的問題，王保利等[17]分析了目前常用的幾類存儲方案，并結(jié)合存儲邊界策略和檢查點技術(shù)，提出了有效存儲邊界策略，在保證波場延拓精度的同時，大幅降低了RTM的存儲需求； 王維紅等[18]提出了針對VSP RTM的PML邊界存儲方法，有效地節(jié)約了存儲空間； 李博等[19]利用CUDA編程架構(gòu)實現(xiàn)CPU/GPU異構(gòu)并行技術(shù)，解決了逆時偏移的計算量大的瓶頸，并給出了偏移算法的實施策略。

本文首先在研究廣角反射波波形變化特征的基礎(chǔ)上，利用τ-p變換進行平面波分解后提取振幅與相位信息，對比分析t-x域、τ-p域的反射波振幅/相位隨入射角變化的特征。由于τ-p變換的限制，在此基礎(chǔ)上，本文利用RTM提取角度域共成像點道集，并以平層和局部Marmousi-2模型為例說明和實現(xiàn)了廣角AVA/PVA反演，反演結(jié)果也進一步驗證了廣角反演的有效性。

2 廣角反射波變化特征研究

Zoeppritz方程精確地描述了平面波在彈性分界面的反射及透射變化特征，由于方程繁雜，其中不僅包含反射、透射變化，還存在橫波及轉(zhuǎn)換波等不同波型的變化。因此，本文以PP波反射情況為研究對象，為方便起見，采用以下PP波反射系數(shù)公式[20]

Rpp=

(1)

其中

A=(r4-Q)2T1T3

B=(r4-Q-1)2T0T1T2T3

由于廣角反射波僅出現(xiàn)在臨界角之后，根據(jù)Castagna等[21]提出的四類AVO響應(yīng)，本文從中選取Ⅰ類和Ⅲ類為例分析廣角反射波的變化特征,其對應(yīng)模型的彈性參數(shù)如表1所示。

表1 Ⅰ類和Ⅲ類AVO模型的彈性參數(shù)

根據(jù)式(1)計算上述兩類模型的反射系數(shù)—振幅及相位隨入射角變化曲線，如圖1所示，通過對比可以得出以下結(jié)論。

(1)Ⅰ類模型的AVO響應(yīng)曲線：在小入射角情況下，振幅隨著入射角增大而變化，相位基本保持不變。當入射角大于第一臨界角時，振幅存在突變；相位則隨著入射角增加而變?yōu)樨撝? 在第二臨界角位置，由于相位角度范圍限制在-180°～180°,相位產(chǎn)生突變。

(2)Ⅲ類模型的AVO相應(yīng)曲線：由于該模型中不存在臨界角，相位基本保持不變，只有振幅隨著入射角的增加而變化。

為了直觀地分析廣角地震反射特征，下面給出根據(jù)Ⅰ類模型生成的廣角反射記錄，如圖2所示。

圖1 Ⅰ類和Ⅲ類AVO模型的反射振幅(a)與相位(b)隨入射角變化曲線

圖2 Ⅰ類AVO模型的合成廣角反射記錄

從圖中可以看出：廣角記錄所表現(xiàn)的特征與平面波反射系數(shù)特征一致，廣角的反射振幅比臨界前強，且存在明顯的相位變化特征。

3 廣角反射的τ -p域分析

在小入射角情況下，平面波反射系數(shù)與球面波反射系數(shù)差異較小； 但當入射角大于臨界角時，平面波反射波系數(shù)(PRC)不再適用于遠炮檢距t-x域球面波記錄。正如前文所述，通過將平面波分解，可以為廣角信息在τ-p域利用平面波反射系數(shù)進行反演提供可能。本文通過波動方程正演模擬得到Ⅰ類模型的平層介質(zhì)共炮點t-x域記錄(圖3)，并利用τ-p變換實現(xiàn)了平面波分解，得到對應(yīng)的τ-p域記錄(圖4)。

通過提取τ-p域不同入射角的反射波記錄，將提取波形按角度排列得到如圖5所示的角道集，并結(jié)合前文合成的理論角道集，可以看出，兩者無論是臨界角所處的位置，還是振幅和相位變化特征都基本吻合。

與t-x域相比，τ-p域廣角反射振幅與相位更接近于平面波反射系數(shù)。圖5中道集的振幅與相位隨入射角變化的特征曲線分別如圖6、圖7所示。由于在計算反正切值時的角度范圍限制，導致了圖1中相位在第二臨界角附近產(chǎn)生突變。為了降低該突變對反演結(jié)果的影響，圖7中進行了相位展開處理，使得相位變化更加連續(xù)，并給出了與t-x域提取振幅與相位及理論反射波振幅和相位曲線對比。

圖3 Ⅰ類AVO模型共炮點記錄的垂向速度分量(局部)

圖4 Ⅰ類AVO模型共炮點記錄的τ -p域波場

圖5 τ -p域選排后的PP波廣角反射記錄(0°～70°)

圖6 Ⅰ類AVO模型的廣角振幅隨入射角變化曲線

圖7 Ⅰ類AVO模型的廣角相位隨入射角變化曲線

從圖中不難發(fā)現(xiàn)：在入射角較小時，t-x域反射波振幅和相位與理論曲線基本一致，但隨著入射角的增大，提取的振幅和相位曲線與平面波理論反射曲線存在明顯誤差。與之相反，τ-p域提取的振幅和相位信息在廣角則能保持較高的精度。由于大角度的反射的振幅更強，相位變化特征也更加明顯，這在一定程度上也有利于提高介質(zhì)參數(shù)反演的穩(wěn)定性和精度。

4 RTM角道集AVA/PVA反演

通過廣角反射的τ-p域分析，在大入射角情況下，τ-p域的反射特征與理論結(jié)果基本一致。但τ-p變換要求規(guī)則的空間采樣間隔和橫向為均勻介質(zhì)[10]，這與實際情況并不完全符合。因此，本文利用RTM方法獲取角度域共成像點道集(ADCIGs)，并將其應(yīng)用于廣角AVA/PVA反演。

由于RTM是基于雙程波波動方程，與射線類和單程波類成像方法相比，RTM成像結(jié)果更為精確。因此，RTM是目前主流的提取ADCIGs的方法，基于RTM提取ADCIGs主要分為基于方向矢量類方法、局部平面波分解類方法和局部移動成像條件法的三大類[22]。

為了獲取可靠的角度域共成像點道集，并將其應(yīng)用于廣角AVA/PVA反演，本文采用Zhang等[23]提出的ADCIGs提取方法。ADCIGs提取方法的主要難點在于共成像點位置的局部入射角的計算，即波場傳播方向角度與反射界面角度的法向之差。

4.1 局部入射角的計算

由于平面P波的位移、速度及加速度的極化方向平行于傳播方向，因此波場傳播方向可以通過P波波場的梯度得到，對于彈性波波動方程，P波波場可以通過對波場求散度進行分離，即

(2)

式中：U表示地震波場的位移；λ和μ表示彈性參數(shù)；ρ表示密度。

在二維情況下，其傳播方向余弦(γ1,γ2)可表示為

(3)

(4)

一般情況下，波場傳播方向角度α為

(5)

對于反射界面法向，可以將其RTM成像結(jié)果進行Hilbert變換，再利用變換結(jié)果的瞬時相位的空間梯度(瞬時波數(shù))得到[24]。假定s(x,z)為偏移圖像，與其Hilbert變換q(x,z)構(gòu)成的復信號c(x,z)可表示為

c(x,z)=s(x,z)+iq(x,z)

(6)

對于復數(shù)圖像信息c(x,z)，對其瞬時相位φ求取空間梯度可以得到瞬時波數(shù)，空間梯度計算公式為

(7)

式中w1和w2分別表示反射體瞬時波數(shù)的水平和垂直分量。由此可估算反射體的法向角度β。在二維情況下，反射體法向角度β為

(8)

通過式(5)和式(8)可求得波場傳播角度α及反射體法向角度β。如圖8所示，在二維情況下，MM′表示反射界面，則成像點處的入射角θ為

θ=α-β

(9)

圖8 局部入射角計算示意圖

4.2 RTM角道集的提取

為了獲取較精確的反射波振幅與相位信息,成像條件的選擇至關(guān)重要。根據(jù)Sandip等[25]對成像條件的保幅性研究，振幅比成像條件具有更高的成像分辨率，能夠真實地反映反射系數(shù)的大小。因此，采用振幅比成像條件，其成像時刻為成像點的激發(fā)時刻。

利用彈性波一階應(yīng)力—速度方程式的高階差分計算進行震源波場的正推和檢波點波場的逆時延拓。在震源/檢波點的波場延拓過程中，利用散度算子提取彈性波場中P波分量用于成像[26]，與常規(guī)的水平/垂直分量成像相比，其成像振幅更加精確。在每一個共炮點記錄的波場反傳的同時，一方面根據(jù)成像時刻及振幅比成像條件計算成像點的成像值；另一方面，需要保存當前炮在該成像點成像時刻的傳播角度α,為下一步提取角道集做準備。

需要注意的是，在成像過程中，幾何效應(yīng)導致波場在垂向上存在拉伸畸變[27,28]，嚴重降低了廣角信息的分辨率和精度。因此在進行成像之前首先對子波進行了拉伸校正，經(jīng)過校正后的廣角信息有了明顯改善。最后，通過對所有炮的偏移結(jié)果進行疊加，得到最終偏移成像結(jié)果。

圖9 RTM角道集計算的流程

根據(jù)式(8)，結(jié)合最終成像剖面可以求得每個成像點位置的法向角度β，進而求得局部入射角θ。在每一炮內(nèi)，根據(jù)入射角進行角度分選，提取成像值，按照入射角度進行選排，即可得到ADCIGs。ADCIGs提取的總體流程如圖9所示。一般情況下，由于角道集的離散角度與真實反射角存在差異， 因此需要對抽取的角道集進行平滑處理，其平滑公式為

(10)

式中：E表示抽取的角度域共成像點道集；kθ表示角度道中的離散入射角； Image表示成像值；θ表示當前時刻成像點的入射角；σ表示高斯函數(shù)的方差；Num_shot表示偏移總炮數(shù)。

4.3 AVA/PVA反演策略

本文利用模擬退火法的反演方法進行廣角AVA/PVA反演試驗工作，模擬退火法反演有不依賴于初始模型、不需要求解雅克比矩陣、且不易陷入局部最優(yōu)解等優(yōu)點，在地球物理資料反演中有許多成功實例。

根據(jù)式(1),PP波反射系數(shù)RPP可以表示為入射角θ及r1、r2、r3、r4的函數(shù)

RPP=RPP(θ,r1,r2,r3,r4)

(11)

由于r1、r2、r3分別表示利用vP1規(guī)范化后的vP2、vS1、vS2；r4是利用ρ1規(guī)范化后的ρ2，因此式(11)中只有入射角θ為獨立變量。通常情況下，可通過層析成像獲得上層介質(zhì)縱波速度vP1，而ρ1則可根據(jù)加德納(Gardner)公式計算，或從實驗室測定得到。

當已知vP1和ρ1參數(shù)時，對其他四個參數(shù)進行反演，在L2范數(shù)意義下，目標函數(shù)φ(x)可以表示為

(12)

5 模型試算

5.1 層狀介質(zhì)模型

為了說明RTM角道集抽取和廣角反射振幅與相位信息提取方法的有效性，以層狀介質(zhì)模型為例(圖10)，模型尺寸為8km×2km，上/下介質(zhì)縱、橫波速度(km/s)及密度(g/cm3)分別為4.0、1.54、2.3和2.0、0.88、2.0。常規(guī)的逆時偏移包括震源波場的正傳和檢波點的逆時延拓兩個波場模擬過程，采用高階交錯網(wǎng)格有限差分法實現(xiàn)波場的模擬，其中炮點與檢波器均勻布設(shè)在水平地表。圖11為炮點位于4km處的單炮記錄的速度水平分量與垂直分量。在切除直達波后，該記錄作為波場逆時延拓的邊界條件。

圖10 單一層狀介質(zhì)模型

圖11 切除直達波后的單炮記錄水平分量(a)及垂直分量(b)

通過前文所述方法得到層狀介質(zhì)的單炮偏移結(jié)果(圖12a)及經(jīng)過子波拉伸校正[29]的結(jié)果(圖12b)。然后根據(jù)入射角信息，從多炮偏移結(jié)果中提取共成像點的不同角度反射信息獲取ADCIGs(圖13)。

圖12 PP波偏移結(jié)果拉伸校正前(a)、后(b)對比

圖13 水平位置在4km處的ADCIGs

由于廣角子波波形存在畸變，如果提取峰值振幅作為反射系數(shù)，會存在較大的誤差。此處采用求取子波包絡(luò)峰值的方法，降低波形變化帶來的誤差。另外，相位則利用互相關(guān)法計算，通常垂直入射/反射的波形常常也并不是零相位，因此，通常以零度入射的相位作為參照子波進行相位的相關(guān)計算。通過對ADCIGs中提取的振幅信息和相位隨入射角變化曲線(圖14)可以看出：實際振幅在小角度與理論值基本吻合，大角度由于球面波效應(yīng)，兩者的差異逐漸增加；相位則基本保持不變。這與前文所述Ⅰ類AVO曲線特征基本一致。

圖14 ADCIGs提取的振幅和相位信息與理論值對比

5.2 Marmousi-2模型

為了驗證本文提出的廣角AVA/PVA聯(lián)合反演的有效性，以局部Marmousi-2模型(圖15)為例,將模型的尺寸設(shè)為2.0km(寬)×0.88km(高)。由于相位變化僅存在于臨界角后，因此將目的層上、下層介質(zhì)的縱波速度(km/s)、橫波速度(km/s)和密度(g/cm3)分別填充為2.0、0.88、2.0和4.0、1.54、2.3。逆時偏移的模擬參數(shù)分別為：震源為40Hz主頻Ricker子波;采樣間隔為0.4ms;空間采樣間隔均為4m；記錄時間為1400ms；炮數(shù)為60個;檢波點個數(shù)為501個，均勻分布于地表。

在震源正向延拓過程中，根據(jù)每個成像點的激發(fā)時間，結(jié)合式(5)獲得震源縱波傳播方向角度。圖16為震源位于1.2km處的P波傳播方向角度。

在檢波點逆時延拓過程中，根據(jù)成像點的激發(fā)時刻，應(yīng)用振幅比成像條件進行成像，完成單炮偏移，通過疊加所有單炮結(jié)果獲得最終的偏移剖面，再由式(7)獲得其復圖像信號相位得到瞬時波數(shù)，進而計算反射體法向角度(圖17)。已知震源傳播方向角度和反射體法向角度，即可計算在當前炮中該模型的入射角分布情況(圖18)。

圖15 Marmouis-2模型(紅色框區(qū)為示例模型)

圖16 水平距離1km處震源的P波傳播方向角度

圖17 局部Marmousi-2模型的反射體法向角度

在多炮偏移中，不同位置的炮點可以提供目標點的不同角度反射信息。在確定該角度后，結(jié)合式(10)可以從不同位置炮點的偏移結(jié)果中抽取角度域共成像點道集。圖19所示分別為原始模型與修改模型的5組角道集，在原始模型中(圖19a)，由于速度差異較小， 目標深度上層反射層次清晰。為了凸顯廣角反射特征，修改后的模型(圖19b)在目標深度速度差異較大，對上層角道集產(chǎn)生了一定的影響。

圖18 局部Marmousi-2模型的入射角

從目標層位的角道集中可以提取廣角反射波的振幅和相位信息，如圖19所示。本節(jié)將其應(yīng)用于中廣角反演中，通常情況下，廣角畸變的存在導致難以獲取精確的大角度反射信息。本文假定在上層縱波速度vP1及密度ρ1已知的條件下，利用廣角振幅和相位進行其他彈性參數(shù)的反演。

將從角道集中提取的廣角振幅與相位信息帶入式(12)進行迭代反演。本文將AVA、PVA及聯(lián)合反演的結(jié)果進行對比，如圖20所示。由于受到所提取ADCIGs質(zhì)量的影響，本文僅采用入射角為0°～50°的振幅與相位數(shù)據(jù)。

圖19 模型修改前(a)、后(b)的ADCISs，角度范圍為-90°～90°

圖20 模型的AVA/PVA反演結(jié)果

表2 三種反演結(jié)果相對誤差對比

表2給出了上述三類方法的反演結(jié)果與模型真實參數(shù)的相對誤差值。通過對比分析反演結(jié)果可知：反射波振幅主要取決于上、下界面的波阻抗比值，而AVA反演又依賴于振幅信息，因此vP2和ρ2的相對誤差呈符號相反的特征，但是波阻抗的反演結(jié)果較為準確。由于振幅信息受透射及波前擴散影響較大，AVA反演結(jié)果并不十分理想； 相比較而言，相位信息較為可靠，其對vS2反演的效果非常好。同時，研究還發(fā)現(xiàn)，由于相位對密度變化極不敏感，因而導致ρ2的反演誤差較大；在AVA/PVA聯(lián)合反演中，觀測數(shù)據(jù)與反演結(jié)果擬合較好，通過表2也可以看出，聯(lián)合反演的結(jié)果較兩者都更接近于模型參數(shù)。

6 結(jié)論和認識

通過對廣角地震反射波振幅/相位的分析和研究，得到以下結(jié)論和認識。

(1)在反射波臨界角存在的情況下，隨入射角的增大，反射波振幅隨之變化的同時，相位也并不總是零相位，在大入射角時，甚至會出現(xiàn)極性反轉(zhuǎn)的情況。

(2)地震波形信息包含了振幅與相位，特別是在廣角條件下，相位特征更明顯。在波形反演中隱含了反射波的振幅和相位信息，本文則以AVA/PVA反演的方式利用了廣角反射的波形信息。

(3)通過τ-p域的分析結(jié)果可以看出，由于Zoeppritz方程以平面波入射為假設(shè)條件，其并不適用于廣角地震反射波特征的分析。本文通過模型說明了廣角τ-p域分析的可行性，但由于τ-p變換尚存不足之處，τ-p域的廣角反演還有待于進一步研究。

(4)本文利用逆時偏移獲取角度域共成像點道集，并將提取的廣角振幅與相位信息用于AVA/PVA反演，模型實例驗證了該方法的可行性。

(5)反演結(jié)果對比發(fā)現(xiàn)：AVA/PVA的反演結(jié)果各有所長，AVA對縱波速度/密度反演效果較好，PVA則對橫波速度反演較好，聯(lián)合反演在一定程度上彌補了兩者的不足。由此可見，在多波多分量采集、大炮檢距和超臨界角采集中，廣角反射的資料具有一定的應(yīng)用價值。

基于逆時偏移提取的角道集質(zhì)量直接影響了廣角反演的結(jié)果，尤其是在復雜模型中，激發(fā)時刻的計算及基于振幅比成像條件也存在不穩(wěn)定性，而且廣角地震反射特征復雜，存在嚴重的拉伸效應(yīng)；另外，折射首波很可能與一次反射波相混淆。這就導致了基于方向矢量類的反射角計算方法并不能完全適用。因此，逆時偏移方法中入射角的計算和基于相關(guān)成像條件下的拉伸校正、角道集提取和保幅性的研究是下一步的工作方向。

[1] 王志,賀振華,黃德濟等.廣角反射波場特征研究及正演模擬分析.地球物理學進展,2003,18(1):116-121. Wang Zhi,He Zhenhua,Huang Deji et al.The wave-field features researching of wide-angle reflection and the analyzing of forward modeling.Progress in Geophysics,2003,18(1):116-121.

[2] 孫成禹,倪長寬,李勝軍等.廣角地震反射數(shù)據(jù)特征及校正方法研究.石油地球物理勘探,2007,42(1)：24-29. Sun Chengyu,Ni Changkuan,Li Shengjun et al.Feature of wide-angle reflection data and correction me-thod.OGP,2007,42(1):24-29.

[3] 孟憲軍,劉福平,王玉梅.地震反射系數(shù)的相角變化.石油地球物理勘探,2010,45(增刊1):121-124. Meng Xianjun,Liu Fuping,Wang Yumei.The phase angle change for seismic reflection coefficient.OGP,2010,45 (S1):121-124.

[4] 代榮獲,張繁昌,劉漢卿.疊前地震資料AVA變參數(shù)反演方法.地球物理學進展,2015,30(1):261-266. Dai Ronghuo,Zhang Fanchang,Liu Hanqing.AVA inversion using replacing parameters method for pre-stack seismic data.Progress in Geophysics,2015,30(1):261-266.

[5] 代榮獲,張繁昌,劉漢卿等.基于貝葉斯理論的逐次迭代非線性AVA反演方法.吉林大學學報(地球科學版),2014,44(6):2026-2033. Dai Ronghuo,Zhang Fanchang,Liu Hanqing et al.Non-linear pre-stack seismic AVA inversion based on Bayesian theory using successive iteration method.Journal of Jilin University(Earth Science Edition),2014,44(6):2026-2033.

[6] Haase A B.Plane waves,spherical waves and angle-dependent P-wave reflectivity in elastic VTI-models.SEG Technical Program Expanded Abstracts,2002,21:304-307.

[7] Ursenbach C P,Haase A B.Efficient spherical-wave AVO modeling.The Leading Edge,2007,26(12):1584-1589.

[8] Ayzenberg M,Tsvankin I,Aizenberg A et al.Effective reflection coefficients for curved interfaces in transversely isotropic media.Geophysics,2009,74(5):33-53.

[9] van der Baan M and Smit D.Amplitude analysis of isotropic P-wave reflections.Geophysics,2006,71(6):93-103.

[10] Zhu Xinfa,McMechan G.Analysis of phase variation with ray parameter in theτ-pdomain.SEG Technical Program Expanded Abstracts,2012,31:1-5.

[11] 王保利,高靜懷,陳文超等.逆時偏移中用Poynting矢量高效地提取角道集.地球物理學報,2013,56(1):262-268. Wang Baoli,Gao Jinghuai,Chen Wenchao et al.Extracting efficiently angle gathers using Poynting vector during reverse time migration.Chinese Journal of Geophysics,2013,56(1):262-268.

[12] 楊午陽,張厚柱,撒利明等.逆時偏移關(guān)鍵問題探討.石油地球物理勘探，2016,51(6):1251-1262. Yang Wuyang,Zhang Houzhu,Sa Liming et al.Some issues about reverse time migration.OGP,2016,51(6):1251-1262.

[13] 吳國忱,秦海旭.各向同性彈性波隨機邊界逆時偏移.地球物理學進展,2014,29(4):1815-1821. Wu Guochen,Qin Haixu.Elastic reverse time migration in isotropic medium based on random boundary.Progress in Geophysics,2014,29(4):1815-1821.

[14] 薛東川.幾種疊前逆時偏移成像條件的比較.石油地球地球物理勘探,2013,48(2):222-227. Xue Dongchuan.Imaging condition comparison of pre-stack reverse time migration.OGP,2013,48(2):222-227.

[15] 張智,劉有山,徐濤等.彈性波逆時偏移中的穩(wěn)定激發(fā)振幅成像條件.地球物理學報,2013，56(10):3523-3533. Zhang Zhi,Liu Youshan,Xu Tao et al.A stable excitation amplitude imaging condition for reverse time migration in elastic wave equation.Chinese Journal of Geophysics,2013,56(10):3523-3533.

[16] 陳康,吳國忱.逆時偏移拉普拉斯算子濾波改進算法.石油地球物理勘探,2012,47(2):249-255. Chen Kang,Wu Guochen.An improved filter algorithm based on Laplace operator in reverse-time migration.OGP,2012,47(2):249-255.

[17] 王保利,高靜懷,陳文超等.地震疊前逆時偏移的有效邊界存儲策略.地球物理學報,2012,55(7):2412-2421. Wang Baoli,Gao Jinghuai,Chen Wenchao et al.Efficient boundary storage strategies for seismic reverse time migration.Chinese Journal of Geophysics,2012,55(7):2412-2421.

[18] 王維紅,郭雪豹,石穎等．VSP逆時偏移及其存儲策略研究．石油地球物理勘探,2015,50(1)：77-83. Wang Weihong,Guo Xuebao,Shi Ying et al．VSP reverse time migration and its data storage strategy．OGP,2015,50(1)：77-83．

[19] 李博,劉紅偉,劉國峰等.地震疊前逆時偏移算法的CPU/GPU實施對策.地球物理學報,2010,53(12):2938-2943. Li Bo,Liu Hongwei,Liu Guofeng et al.Computational strategy of seismic reverse time migration on CPU/GPU.Chinese Journal of Geophysics,2010,53(12):2938-2943.

[21] Castagna J P and Swan H W.Principles of AVO cross-plotting.The Leading Edge,1997,16(4):337-344.

[22] Jin H,McMechan G A and Guan H M.Comparison of methods for extracting ADCIGs from RTM.Geophysics,2014,79(3):89-103.

[23] Zhang Q S and McMechan G A.Direct vector-field method to obtain angle-domain common-image gathers from isotropic acoustic and elastic reverse time migration.Geophysics,2011,76(5):135-149.

[24] Havlicek J P,Havlicek J W,Mamuya N D et al.Skewed 2D Hilbert transforms and computed AM-FM models.IEEE International Conference on Image Proce-ssing,1998,602-606.

[25] Chattopadhyay S,McMechan G A.Imaging conditions for pre-stack reverse-time migration.Geophysics,2008,73(3):81-89.

[26] 李振春,張華,劉慶敏等.彈性波交錯網(wǎng)格高階有限差分法波場分離數(shù)值模擬.石油地球物理勘探,2007,42(5):249-255. Li Zhenchun,Zhang Hua,Liu Qingmin et al.Numeric simulation of elastic wavefield separation by staggering grid high-order finite difference algorithm.OGP,2007,42(5):249-255.

[27] Tygel M,Schleicher J and Hubral P.Pulse distortion in depth migration.Geophysics,1994,59(10):1561-1569.

[28] Zhu Xinfa and McMechan G A.Stretch-free migration imaging condition.SEG Technical Program Expanded Abstracts,2013,32:3882-3887.

[29] 楊佳佳，何兵壽，陳婷．逆時深度偏移中的子波拉伸校正．石油地球物理勘探，2016，51(1)：135-140. Yang Jiajia，He Bingshou，Chen Ting．Wavelet stretch correction for reverse-time migration．OGP，2016，51(1)：135-140．

(本文編輯：劉英)

武泗海 碩士研究生，1991年生；2015年畢業(yè)于西南石油大學勘查技術(shù)與工程專業(yè)，獲學士學位；現(xiàn)在西南石油大學地球科學與技術(shù)學院攻讀碩士學位，主要從事地震采集、偏移成像與高性能計算的研究。

1000-7210(2017)05-1005-11

P631

A

10.13810/j.cnki.issn.1000-7210.2017.05.014

*四川省成都市新都區(qū)新都大道8號西南石油大學地球科學與技術(shù)學院，610500。Email：wsh_18980796471@sina.com

本文于2016年9月20日收到，最終修改稿于2017年8月6日收到。

本項研究受國家重大專項(2016ZX05024-001-003)、中國石油天然氣集團公司科學研究與技術(shù)開發(fā)項目“深層與非常規(guī)物探新方法新技術(shù)”(2016A-33)聯(lián)合資助。