張 林， 楊飛龍， 史朝陽(yáng)， 江 桂， 陳繼川

(1.陜西省地質(zhì)調(diào)查中心，西安 710068；2. 西安石油大學(xué) 地球科學(xué)與工程，西安 710065；3. 陜西省地質(zhì)調(diào)查院，西安 710054)

VTI介質(zhì)井間地震高斯束正演數(shù)值模擬

張 林1， 楊飛龍2， 史朝陽(yáng)1， 江 桂1， 陳繼川3

(1.陜西省地質(zhì)調(diào)查中心，西安 710068；2. 西安石油大學(xué) 地球科學(xué)與工程，西安 710065；3. 陜西省地質(zhì)調(diào)查院，西安 710054)

為了研究垂直橫向各向同性(VTI)介質(zhì)井間地震中,地震波的傳播方向及波場(chǎng)特征，采用高斯束方法對(duì)復(fù)雜構(gòu)造進(jìn)行波場(chǎng)正演模擬。分析VTI介質(zhì)中地震波的群相關(guān)系，并提出群相數(shù)據(jù)庫(kù)思想來(lái)解決地震波在VTI介質(zhì)中傳播的方向和速度問(wèn)題，改善了波場(chǎng)正演模擬的算法。同時(shí)高斯束正演方法不僅能夠解決復(fù)雜構(gòu)造的盲區(qū)問(wèn)題，且其動(dòng)力學(xué)特征更能反映各向異性介質(zhì)對(duì)地震波傳播的影響。較VTI介質(zhì)頻率-空間域彈性波正演模擬方法，正演結(jié)果一致且計(jì)算效率高，能有效分辨地震波場(chǎng)特征，為地震數(shù)據(jù)的處理與解釋提供了借鑒。

井間地震； 群相數(shù)據(jù)庫(kù)； 垂直橫向各向同性介質(zhì)； 高斯束； 正演模擬

0 引言

井間地震[1]因其高精度和高分辨率以及貼近目標(biāo)探測(cè)的特點(diǎn)和優(yōu)勢(shì)，將油氣勘探開(kāi)發(fā)帶入了精細(xì)勘探時(shí)代。井間地震勘探是將激發(fā)系統(tǒng)與接收系統(tǒng)分別放置于不同井中進(jìn)行數(shù)據(jù)采集，可以獲得多種類型的地震波場(chǎng)信息。豐富的波場(chǎng)信息為成像提供了大量的信息基礎(chǔ)，但同時(shí)多種波型也會(huì)產(chǎn)生相互干擾，使得地震波場(chǎng)難以分辨。地震波在各向異性介質(zhì)中傳播時(shí)，地震波場(chǎng)特征還會(huì)受到介質(zhì)彈性參數(shù)的變化而發(fā)生變化，對(duì)有效波場(chǎng)地識(shí)別與分離造成更大困難。地震波場(chǎng)的識(shí)別與分離方法有很多，其中井間地震正演數(shù)值模擬方法是最有效、最直觀的方法之一，不僅可以幫助井間地震任務(wù)的設(shè)計(jì)、制定適合的野外采集方案、指導(dǎo)野外采集工作的進(jìn)行，還可以識(shí)別井間復(fù)雜波場(chǎng)特征，指導(dǎo)地震數(shù)據(jù)處理工作的順利開(kāi)展。

井間地震正演數(shù)值模擬方法主要包括波動(dòng)方程類和射線類正演方法。波動(dòng)方程類正演方法的優(yōu)點(diǎn)在于能夠獲得較完整的波場(chǎng)信息，但是計(jì)算效率低，尤其是面對(duì)大規(guī)模地震勘探及三維地質(zhì)模型。張文波等[2]使用交錯(cuò)網(wǎng)格有限差分方法對(duì)井中地震彈性波進(jìn)行正演模擬，以射線理論為基礎(chǔ)的正演方法在計(jì)算效率上都優(yōu)于波動(dòng)方程類正演方法，但是面對(duì)物性變化較大的復(fù)雜構(gòu)造，波場(chǎng)信息不完整，出現(xiàn)盲區(qū)。cerveny等[3]提出了使用高斯射線束方法進(jìn)行地震波場(chǎng)正演數(shù)值模擬，高斯束是對(duì)波動(dòng)方程的高頻近似，它是將地震波場(chǎng)分解到一定頻率范圍內(nèi)的射線束上來(lái)實(shí)現(xiàn)波場(chǎng)的數(shù)值模擬，不僅包含波場(chǎng)傳播的運(yùn)動(dòng)學(xué)特征，還具有地震波傳播的動(dòng)力學(xué)變化。并且具有高效、精確等特點(diǎn)，適用于復(fù)雜地質(zhì)構(gòu)造，對(duì)復(fù)雜構(gòu)造的盲區(qū)、臨界區(qū)等具有較好效果。楊飛龍等[4]將高斯束正演方法應(yīng)用在各向同性介質(zhì)的井間地震波場(chǎng)正演數(shù)值模擬當(dāng)中。針對(duì)射線類正演方法計(jì)算效率快、波動(dòng)方程類正演方法精度高的優(yōu)點(diǎn)，筆者在研究VTI介質(zhì)井間地震正演時(shí)，結(jié)合二者優(yōu)點(diǎn)采用高斯射線束方法進(jìn)行正演數(shù)值模擬。

1 方法技術(shù)

1.1 群相數(shù)據(jù)庫(kù)

在地震勘探中，我們所研究的地球介質(zhì)因地質(zhì)作用影響出現(xiàn)各向異性，介質(zhì)的各向異性表現(xiàn)在地震波的傳播速度與衰減方向上。各向異性介質(zhì)中，地震波沿著同一點(diǎn)出發(fā)的不同方向上介質(zhì)物理性質(zhì)不同，主要體現(xiàn)在群速度和相速度上，由于各向異性的作用使得群速度與相速度分離，導(dǎo)致群角與相角也產(chǎn)生分離。

各向異性介質(zhì)中，群速度代表了地震波能量的傳播速度，體現(xiàn)在波場(chǎng)動(dòng)力學(xué)特征，也與地震波的射線路徑有關(guān)[5]；相速度表示的是波矢量的傳播速度，方向與波前面垂直，它是求解簡(jiǎn)諧平面波解時(shí)假定的速度，通過(guò)求解Christoffel方程獲取[6]。群角是指介質(zhì)的對(duì)稱軸與射線傳播方向之間的夾角；而相角指的是介質(zhì)的對(duì)稱軸方向和波矢量方向的夾角，群速度的方向就是群角的方向，相速度的方向與相角的方向一致(圖1)。從圖1中可以看到，波矢量的方向垂直于波前面，相速度的方向即為波前面的傳播方向，群速度的方向與波射線的方向一致。波前面在各向異性介質(zhì)下為一橢球面，當(dāng)相角不等于群角時(shí)，相速度與群速度也不相等。

圖1 各向異性介質(zhì)的群相關(guān)系Fig.1 Group-phase relations of anisotropic medium

Daley等[7]根據(jù)VTI介質(zhì)的應(yīng)力與應(yīng)變的關(guān)系，推導(dǎo)出地震波在橫向各向同性介質(zhì)中傳播的相速度，使用Thomsen參數(shù)來(lái)表示彈性參數(shù)并將介質(zhì)近似為弱各向異性介質(zhì)，其地震波傳播的相速度可表示為：

(1)

式中：θ為相角；ε為P波的各向異性參數(shù)，ε越大，介質(zhì)的各向異性越強(qiáng)；δ為影響VTI介質(zhì)中垂直對(duì)稱軸附近縱波的速度的參數(shù)；γ表示的是橫波的各向異性參數(shù)，γ與各向異性的強(qiáng)度成正比，當(dāng)γ=0時(shí)，表示介質(zhì)中不存在橫波的各向異性。

Berryman[8]提出了使用相速度計(jì)算群速度的計(jì)算公式：

(2)

當(dāng)入射角為0°或者90°時(shí)，群速度與相速度相等。

地震波在各向異性介質(zhì)中傳播時(shí)，波矢量方向與地震波傳播方向不一致，因此群角和相角也不相等。Thomsen[8]給出了弱各向異性介質(zhì)中地震波傳播的群角和相角之間的關(guān)系為：

(3)

將式(3)分別寫成P波、SV波及SH波的形式：

(4)

式中：vp、vsv、δ、ε、γ為Thomsen參數(shù)；θ是相角；φ是群角。

地震波在各向異性介質(zhì)中傳播時(shí)，群角為射線角，當(dāng)射線遇到地層界面發(fā)生反射或者透射時(shí)地震波的入射角用相角來(lái)表示。式(4)為相角計(jì)算群角的計(jì)算公式，可是射線追蹤過(guò)程中我們常遇到已知群角計(jì)算相角的問(wèn)題。Byun等[9]均對(duì)該問(wèn)題進(jìn)行了研究，提出了采用相角計(jì)算群角的近似公式。趙愛(ài)華等[10]將弱各向異性介質(zhì)中相角與群角的關(guān)系表示為θ=φ-g(φ)，由群角計(jì)算相角需要進(jìn)行近似取舍，并且運(yùn)算過(guò)程比較復(fù)雜，計(jì)算效率低?；诖?，本次研究中仍以式(4)為基礎(chǔ)，首先計(jì)算出任意相角對(duì)應(yīng)的群角大小，然后計(jì)算出所對(duì)應(yīng)的相速度和群速度，將它們保存在一個(gè)數(shù)據(jù)庫(kù)中，在此稱為群相數(shù)據(jù)庫(kù)。如圖2所示，為群相數(shù)據(jù)庫(kù)中群角、相角、群速度和相速度之間的關(guān)系。當(dāng)進(jìn)行地震波場(chǎng)正演模擬時(shí)，射線以某一群角入射，遇到地層界面發(fā)生反射或透射時(shí)讀取群數(shù)據(jù)庫(kù)中與之對(duì)應(yīng)的相角，進(jìn)行SNELL定理計(jì)算。高斯束正演模擬時(shí)每個(gè)檢波器的能量是由多條射線能量加權(quán)而來(lái)，在VTI介質(zhì)高斯束正演時(shí),需要進(jìn)行多次SNELL運(yùn)算，若提前計(jì)算相角對(duì)應(yīng)的群角大小，便可在運(yùn)算時(shí)直接調(diào)用數(shù)據(jù)庫(kù)中的信息，避免運(yùn)算過(guò)程中產(chǎn)生的誤差，同時(shí)會(huì)提高運(yùn)算的效率。

圖2 群相關(guān)系示意圖Fig.2 Group-phase relations diagram

圖3 VTI介質(zhì)地震波傳播示意圖Fig.3 Seismic wave propagation diagram in VTI medium

1.2 反射與透射

如圖3所示，地震波在VTI介質(zhì)中傳播時(shí)在界面處發(fā)生透射、反射，仍然遵循斯奈爾定律[11]。但是對(duì)于VTI介質(zhì)，地震波傳播的速度和入射角度都存在群相分離現(xiàn)象，因此較各向同性介質(zhì)來(lái)說(shuō)斯奈爾定律更加復(fù)雜。

在VTI介質(zhì)中，介質(zhì)的對(duì)稱軸方向與垂直方向一致，斯奈爾定律可以表示為：

(5)

式中：P為射線參數(shù)；α為入射角；θ為相角；v(θ)為相速度。

在各向異性介質(zhì)中，可以將反射系數(shù)與透射系數(shù)分為兩部分：①各向同性介質(zhì)性質(zhì)的反射、透射系數(shù)項(xiàng);②各向異性介質(zhì)性質(zhì)的反射、透射系數(shù)項(xiàng)?？梢员硎緸槭?6)。

(6)

Ripp、Tipp、Rips和Tips分別為P波的反射系數(shù)、透射系數(shù)，SV波的反射系數(shù)和透射系數(shù)。

以PP波(P波入射P波反射)為例，反射系數(shù)的各向同性項(xiàng)可以表示為式(7)。

(7)

反射系數(shù)的各向異性項(xiàng)可以表示為：

(8)

其中：vp01和vs01為界面以上的P波速度和SV波速度；vp02和vs02為界面以下的P波速度和SV波速度；ρ1和ρ2分別為界面上下兩種介質(zhì)的密度。

1.3 高斯束正演

高斯束正演是將波場(chǎng)分解到具有一定頻率范圍的射線束上實(shí)現(xiàn)地震波場(chǎng)的數(shù)值模擬，它是將動(dòng)力學(xué)方程集中在射線附近的高頻漸近時(shí)間調(diào)和解。高斯束正演包括運(yùn)動(dòng)學(xué)射線追蹤、動(dòng)力學(xué)射線追蹤和波場(chǎng)疊加三個(gè)階段。

高斯射線束運(yùn)動(dòng)學(xué)追蹤即求解從震源發(fā)出的所有射線經(jīng)過(guò)地層反射、透射后最終到達(dá)接收井的射線路徑、旅行時(shí)和中心射線的振幅。運(yùn)動(dòng)學(xué)射線追蹤方法很多，常用的是求解程函方程的射線追蹤方法。

(9)

式中：τ為位移；v為速度；x、z為水平和垂直坐標(biāo)，其中：x、y、z為直角坐標(biāo)系下空間坐標(biāo)，τ為旅行時(shí)，v為地震波傳播速度。已知各反射界面的分布函數(shù)，炮點(diǎn)的坐標(biāo)和出射角，即射線的函數(shù)表達(dá)式是已知的，當(dāng)按一定角度步長(zhǎng)打出一條條射線后，根據(jù)界面函數(shù)和射線函數(shù)關(guān)系可以求出兩者交點(diǎn)的坐標(biāo)，然后判斷交點(diǎn)是否在我們所設(shè)計(jì)模型的有效范圍內(nèi)。如果在就記錄此交點(diǎn)坐標(biāo)，然后利用Snell定律求反射或透射角，進(jìn)一步求出生成射線方程，再利用上述方法求與其他界面的交點(diǎn)，直至到達(dá)接收井上，求出整條射線路徑的所有控制點(diǎn)坐標(biāo)，最后把各控制點(diǎn)坐標(biāo)輸出到存儲(chǔ)器中；如果不在就直接進(jìn)行下一個(gè)角度的射線追蹤。

經(jīng)過(guò)運(yùn)動(dòng)學(xué)追蹤，就有了中心射線，在此基礎(chǔ)上可以進(jìn)行這高斯射線束的動(dòng)力學(xué)射線追蹤，即求得p、q的函數(shù)值。函數(shù)p(s)、q(s)在高斯射線束中起著非常重要的作用，它們決定了高斯射線束能量的分布狀態(tài)，也表征沿射線傳播方向的高頻地震波場(chǎng)動(dòng)力學(xué)特征。

用U(R,t)表示地震波傳播到接收點(diǎn)R的波場(chǎng)，根據(jù)傅里葉變化得到時(shí)間域波場(chǎng)：

(10)

寫成離散形式為：

(11)

式中：φ為從震源發(fā)出射線的入射角；φ0和φN分別為起始入射角和終止入射角；g為高斯波包；Δφ為入射角的間隔。

根據(jù)文獻(xiàn)[3]，可知波包g的近似解析表達(dá)式為式(12)。

g(R,φ)= (2πfm)1/2|Aφ|exp{-[2πfm(t-θ)/

γ]2+(2πfmG/γ)-2πfmG}·

(12)

式中:f*=fm·(1-4πfmG/γ2)，為高斯波包的主頻。

2 數(shù)值模擬

圖4 VTI介質(zhì)井間地震高斯束正演地質(zhì)模型Fig.4 Geological model of cross-well seismic Gaussian beam forward in VTI medium

表1 VTI介質(zhì)井間地震地質(zhì)模型參數(shù)

表2 井間地震觀測(cè)系統(tǒng)參數(shù)

圖5 VTI介質(zhì)井間地震高斯束正演射線路徑Fig.5 The ray path of cross-well seismic Gauss beam forward in VTI medium

如圖4所示，為VTI介質(zhì)井間地震高斯束正演地質(zhì)模型，模型參數(shù)如表1所示。井間地震高斯束正演的觀測(cè)系統(tǒng)如表2所示。

對(duì)該VTI介質(zhì)的井間地質(zhì)模型使用上述觀測(cè)系統(tǒng)進(jìn)行高斯束正演，得到如圖5所示的射線路徑。井間地震的波場(chǎng)較地面地震和VSP更加復(fù)雜，在VTI介質(zhì)中，地震波傳播的射線路徑也不像均勻介質(zhì)中那么有序，地震波傳播時(shí)遵循VTI介質(zhì)下SNELL定律。圖5中可以看到井間地震波場(chǎng)正演包含直達(dá)P波、直達(dá)S波，以及上行反射P波、上行反射S波與下行反射P波、下行反射S波。在實(shí)際的井間地震數(shù)據(jù)處理中，識(shí)別上下行反射波場(chǎng)是在波場(chǎng)分離中重要的工作之一。透過(guò)正演研究，能夠給井間地震數(shù)據(jù)處理提供一個(gè)依據(jù)，指導(dǎo)數(shù)據(jù)處理進(jìn)一步開(kāi)展。

圖6 VTI介質(zhì)井間地震高斯束正演射線路徑Fig.6 The ray path of cross-well seismic Gauss beam forward in VTI medium

圖7 VTI介質(zhì)井間地震高斯束正演波場(chǎng)記錄Fig.7 The wave record of cross-well seismic Gauss beam forward in VTI medium

圖6為第二炮井間地震高斯束正演的射線路徑，圖7為其相應(yīng)的波場(chǎng)記錄。從圖6中可以看到，在VTI介質(zhì)中地震波傳播時(shí)產(chǎn)生的射線分布不均勻，尤其面對(duì)復(fù)雜地質(zhì)構(gòu)造的時(shí)候會(huì)產(chǎn)生射線的扭轉(zhuǎn)，在地震記錄上也表現(xiàn)為反轉(zhuǎn)的現(xiàn)象。圖7的波場(chǎng)記錄中直達(dá)P波記錄中有空缺的地方，也出現(xiàn)一小段斜率變化較大的地方，這是因?yàn)樵赩TI介質(zhì)中射線傳播時(shí)遇到界面后反射透視所遵循的VTI介質(zhì)SNELL定律。從圖7中的波場(chǎng)記錄上可以看到，VTI介質(zhì)的井間地震波場(chǎng)更加復(fù)雜，在復(fù)雜地質(zhì)構(gòu)造中地震波場(chǎng)表現(xiàn)的沒(méi)有那么整齊，出現(xiàn)了一些小的錯(cuò)段。在實(shí)際井間地震數(shù)據(jù)中，由于井間地震數(shù)據(jù)的復(fù)雜性，再加上復(fù)雜構(gòu)造，對(duì)于識(shí)別分辨波場(chǎng)十分困難，透過(guò)VTI介質(zhì)井間高斯束正演數(shù)值模擬，可以有效模擬實(shí)際地下構(gòu)造井間復(fù)雜波場(chǎng)特征，幫助我們分辨所需波場(chǎng)，進(jìn)行下一步地震數(shù)據(jù)處理。

在實(shí)際井間地震數(shù)據(jù)處理時(shí)，波場(chǎng)分離工作是處理的重要任務(wù)。通過(guò)不同波場(chǎng)類型的正演，可以幫助我們有效地分辨井間地震波場(chǎng)傳播的特點(diǎn)，為進(jìn)一步井間地震數(shù)據(jù)處理提供依據(jù)。圖8(d)中的下行反射P波的時(shí)距曲線上出現(xiàn)一段斜率不同的下行反射P波波場(chǎng)特征，在接收井深度410 m～510 m范圍里。造成此種現(xiàn)象的原因是，地震波在各向異性介質(zhì)中傳播時(shí)受到不同方向上速度差異的影響，在實(shí)際的井間地震資料上有很多類似這種由于地層的各向異性和復(fù)雜構(gòu)造造成的波場(chǎng)特點(diǎn)。

圖8 VTI介質(zhì)井間地震上下行反射波場(chǎng)特征Fig.8 The up and down reflection wave field characteristics of cross-well seismic Gauss beam forward in VTI medium(a)上行S波的射線路徑；(b)對(duì)應(yīng)的正演記錄；(c)下行P波的射線路徑；(d)對(duì)應(yīng)的正演記錄

3 不同井間地震正演方法對(duì)比研究

圖9 VTI介質(zhì)地質(zhì)模型Fig.9 Geological model in VTI medium

為了驗(yàn)證VTI介質(zhì)井間地震高斯束正演算法的準(zhǔn)確性，選擇簡(jiǎn)單的層狀地質(zhì)模型將其與VTI介質(zhì)頻率-空間域彈性波正演模擬方法進(jìn)行對(duì)比。建立如圖9所示的VSP勘探地質(zhì)模型(VSP勘探是一種特殊的井間地震勘探)，模型參數(shù)如表3所示。炮點(diǎn)放置于地面300 m處，偏移距為300 m。檢波點(diǎn)安置在接收井(直井)上，道間距為5 m，共120道。分別使用高斯束方法和波動(dòng)方程方法進(jìn)行正演模擬，得到如圖10所示的波場(chǎng)特征。既包含縱波震源產(chǎn)生的波場(chǎng)，也有橫波震源產(chǎn)生的波場(chǎng)。對(duì)比其縱波震源產(chǎn)生的波場(chǎng)特征，可以看出，VTI介質(zhì)井間地震高斯束正演的波場(chǎng)形態(tài)、旅行時(shí)與高階交錯(cuò)網(wǎng)格彈性波數(shù)值模擬結(jié)果一致，驗(yàn)證了該方法的準(zhǔn)確性。波動(dòng)方程正演模擬方法不僅具有地震波傳播的動(dòng)力學(xué)特征，也能清晰反映地震波在VTI介質(zhì)中傳播的運(yùn)動(dòng)學(xué)特點(diǎn)。將兩種方法對(duì)同一地質(zhì)模型正演效率進(jìn)行對(duì)比，如表4所示，可見(jiàn)在相同電腦配置下，高斯束正演方法的計(jì)算速度要比波動(dòng)方程法正演效率高近乎50倍。在此將兩種方法進(jìn)行對(duì)比，更加說(shuō)明本文所提出的群相數(shù)據(jù)庫(kù)思想應(yīng)用在VTI介質(zhì)井間地震高斯束正演模擬中，能夠改善VTI介質(zhì)正演模擬算法。

圖10 不同正演方法波場(chǎng)對(duì)比示意圖Fig.10 The contrast diagram of wave field between two different forward methods(a)VTI介質(zhì)高斯束正演的波場(chǎng)特征；(b)高階交錯(cuò)網(wǎng)格彈性波數(shù)值模擬方法正演的結(jié)果

表3 VTI介質(zhì)井間地震地質(zhì)模型參數(shù)

表4 單炮正演時(shí)間統(tǒng)計(jì)表

4 結(jié)論

筆者研究了VTI介質(zhì)井間地震高斯束正演方法，提出了使用群相數(shù)據(jù)庫(kù)思想來(lái)解決群角計(jì)算相角的問(wèn)題，有效地改善了VTI介質(zhì)的射線正演計(jì)算方法。高斯束正演方法不僅能解決復(fù)雜地質(zhì)構(gòu)造產(chǎn)生的盲區(qū)問(wèn)題，且作為射線類正演方法具有較高的運(yùn)算效率。將其與井間地震高階交錯(cuò)網(wǎng)格彈性波數(shù)值模擬方法進(jìn)行對(duì)比，運(yùn)算結(jié)果一致，證明筆者所研究的VTI介質(zhì)井間地震高斯束正演方法的正確性，為在井間地震復(fù)雜構(gòu)造采集設(shè)計(jì)與波場(chǎng)分離工作中提供了借鑒。

[1] 何惺華.井間地震[M].北京：石油工業(yè)出版社，2008.

HE X H. Cross-well seismic[M]. Beijing: Petroleum Industry Press，2008.(In Chinese)

[2] 張文波. 井間地震交錯(cuò)網(wǎng)格高階差分?jǐn)?shù)值模擬及逆時(shí)偏移成像研究[D].西安:長(zhǎng)安大學(xué)，2005.

ZHANG W B. The studies of staggered grid high-order finite difference numerical simulation and reverse-time migration imaging in cross-well seismic[D].Xi an: Chang’an University，2005.(In Chinese)

[4] 楊飛龍, 孫淵, 李緒宣,等. 基于高斯射線束的斜井VSP正演方法[J]. 地球物理學(xué)進(jìn)展, 2014, 29(6)：2791-2799.

YANG F L, SUN Y, LI X X, et al. Deviated hole VSP forward method based on the Gaussian beam[J]. Progress in Geophysics, 2014,29(6)：2791-2799.(In Chinese)

[5] 杜麗英, 劉國(guó)明, 杜麗娟,等. VTI介質(zhì)中地震波反射波合成記錄的方法研究[J]. 地球物理學(xué)進(jìn)展, 2001, 2 (16): 58-64.

DU L Y, LIU G M, DU L J, et al. The method of seismic and reflection wave synthetic seismo gram in VTI media[J]. Progress in Geophysics,2001, 2 (16): 58-64. (In Chinese)

[6] DALEY P F, HRON F. Reflection and transmission coefficients for transversely isotropic media [J]. Bulletin of the Seismological Society of America, 1977, 67 (3): 661-675.

[7] BERRYMAN J G. Long-wave elastic anisotropy in transversely isotropic media [J]. Geophysics, 1979, 44 (5): 896-917.

[8] THOMSEN,L.Weadly elostic anisotropy[J].Geophysics,1986,51(10):1954-1966.

[9] BYUN B,CORRIGAN D,GAISER J.Anisotropic velocity analysis for lithology discrimination [J].Geophysics,1989,54 (12):1564-1574.

[10] 趙愛(ài)華, 丁志峰. 一種弱各向異性介質(zhì)地震波群速度的近似表示新方法[J]. 地球物理學(xué)進(jìn)展, 2005, 20 (4): 916-919.

ZHAO A H, DING Z F. New approximate expressions of seismic group velocities for weakly anisotropic media[J]. Progress in Geophysics,2005, 20 (4): 916-919. (In Chinese)

[11] 李建國(guó), 李彥鵬, 郭曉玲. VTI 介質(zhì)試射射線追蹤 [J]. 石油地球物理勘探, 2010, (4): 491-496.

LI J G, LI Y P, GUO X L.VTI medium test-firing ray tracing[J].OGP,2010(4): 491-496. (In Chinese)

Gaussianbeamforwardnumericalsimulationofcross-wellseismicinVTImedium

ZHANG Lin1, YANG Feilong2, SHI Zhaoyang1, JIANG Gui1, CHEN Jichuan3

(1.Shaanxi Center of Geological Survey, Xi'an 710068， China；2.School of Earth Sciences and Engineering, Xi'an Shiyou University, Xi'an 710065， China；3.Shaanxi Geological Survey Institute, Xi'an 710054， China)

The Gaussian beam method is used in wave field forward modeling in complex structure, in order to study the direction of seismic wave and the wave field characteristics in vertical transversely isotropic(VTI) medium in cross-well seismic. The group-phase relations of seismic wave are analyzed in VTI media, the group-phase database thought is put forward to solve the problems of direction and velocity of seismic wave in VTI media, and the algorithm of wave field forward modeling is improved. At the same time, Gaussian beam forward modeling method is not only solving the problem of blind area of complex structure, and its dynamic characteristics can reflect the anisotropy medium effects on seismic wave propagation. Compared with elastic wave forward modeling method in frequency-spatial domain in VTI media, the forward modeling results are consistent and the computational efficiency is higher. The characteristics of seismic wave field can be effective separation. The method provides accurate basis for seismic data processing and interpretation.

cross-well seismic; group-phase database; vertical transversely isotropic medium; gaussian beam; forward modeling

2017-03-16 改回日期： 2017-06-09

張林(1978-)，男，高級(jí)工程師，主要從事綜合地球物理勘探，E-mail：136751912@qq.com。

陳繼川(1970-)，男，高級(jí)工程師，主要從事物探、遙感、測(cè)繪等方面研究，E-mail：527379546@qq.com。

1001-1749(2017)06-0791-08

P 631.4

A

10.3969/j.issn.1001-1749.2017.06.12