VTI巖石界面產(chǎn)生的非均勻折射P波極化狀態(tài)分析

法　林　張曉琳　席冰潔　周　鑫　牟錦鵬　楊　宏

(①西安郵電大學(xué)電子工程學(xué)院，陜西西安 710121； ②西安郵電大學(xué)通信與信息工程學(xué)院，陜西西安 710121)

1　引言

彈性各向異性在地球內(nèi)部是一種常見的現(xiàn)象，許多地球物理學(xué)家已觀察到各向異性介質(zhì)中的極化現(xiàn)象[1-4]并對反射/折射系數(shù)、極化等方面進行了廣泛的研究[5-11]，其研究成果已應(yīng)用到勘探地震數(shù)據(jù)的反演解釋中[12-15]。人們通常把沉積巖(例如頁巖)視為橫向各向同性(TI)介質(zhì)，理論上VTI介質(zhì)力學(xué)性質(zhì)可用六角晶體的彈性張量進行描述[3，4]。

在地球物理勘探中，巖石的各向異性對聲波的極化、傳播和反射/折射有很大影響，正確地評估巖石各向異性的這些影響對地球物理勘探的正演建模和地震勘探數(shù)據(jù)反演至關(guān)重要。在許多已發(fā)表的文獻中，可見各向異性巖層中傳播的均勻平面波的極化特性方面的研究[1，4，16-19]，以及在各向異性巖層界面上產(chǎn)生的非均勻折射P波的極化特性方面的研究，例如，對于存在和不存在異常入射角的兩種情況，法林等[10，20]推導(dǎo)了非均勻折射P波極化系數(shù)的嚴(yán)格數(shù)學(xué)表達式。

對于存在異常入射角的VTI介質(zhì)界面的情況，在異常入射角對非均勻折射P波橢圓極化狀態(tài)的影響下能夠正確地計算出反射/折射系數(shù)的幅度及其相位與入射角之間的變化關(guān)系。如果忽略異常入射角引起的非均勻折射P波對極化狀態(tài)的影響，則不能正確地在過異常入射角區(qū)域計算出VTI巖石界面的PP波的反射系數(shù)以及反射P波相對入射P波在界面上產(chǎn)生的渡越時間。因而對于大角度入射而獲得的地震勘探數(shù)據(jù)就不能進行精確的AVO反演分析，對上述地震勘探數(shù)據(jù)進行時—深轉(zhuǎn)換時也會引起較大的誤差。

本文對各向異性巖石界面上產(chǎn)生的非均勻折射P波的極化狀態(tài)進行了數(shù)值計算，用類似“Poincaré球”的可視化方式描述了巖石各向異性和異常入射角對其橢圓極化狀態(tài)的影響。在VTI介質(zhì)界面的反射系數(shù)和極化系數(shù)及前異常入射角區(qū)域非均勻模式轉(zhuǎn)換波的極化系數(shù)表達式的基礎(chǔ)上，推導(dǎo)了過異常入射角區(qū)域?qū)?yīng)的非均勻模式轉(zhuǎn)換波的極化系數(shù)的表達式，并首次建立了相應(yīng)的橢圓極化方程，同時計算了入射角與橢圓極化軌跡之間的關(guān)系，并論證了這個異常入射角[10]不僅可以使橢圓極化軌跡的形狀發(fā)生突變，而且還會使橢圓極化旋轉(zhuǎn)方向發(fā)生改變。

2　理論背景

對于VTI介質(zhì)，可以用六角晶系固體的彈性剛度張量

(1)

描述其主要力學(xué)性質(zhì)[4，21]。式中Cij與各向異性巖石參數(shù)之間的關(guān)系[22]為

(2)

式中：ε、δ和γ是巖石各向異性參數(shù)；α和β分別是P波和SV波的垂直相速度；ρ是介質(zhì)密度。

對于角頻率為ω的簡諧P波，其在x-z平面上傳播并入射到VTI介質(zhì)界面上(圖1)，入射P波、反射P波、折射P波、反射SV波和折射SV波的質(zhì)點位移矢量的表達式為

(3)

式中：θ為入射角、反射角或折射角；R0定義為1，R表示反射或折射系數(shù)；ux和uz為極化系數(shù)；k為波數(shù)。

圖1　　　　各向異性巖石界面模型的入射P波和模

對于不存在外力的情況，入射P波和在界面上產(chǎn)生的模式轉(zhuǎn)換波的極化系數(shù)ux和uz可以通過求解簡化的Christoffel方程

(4)

獲得。式中

(5)

其中：A1=A11-A44；A2=A44-A33；A3=A13+A44；vm是五種波的相速度解[20]。在不存在異常入射角的情況下，對于入射、反射和折射P波式(5)中Qm取“+”，對于反射、折射SV波Qm取“-”；在存在異常入射角的情況下，在過異常入射角區(qū)域，對于折射P波Qm取“-”。

3　模型建立

根據(jù)測量出的巖石各向異性參數(shù)[22-24]，經(jīng)大量計算后發(fā)現(xiàn)對于大多數(shù)各向異性巖石構(gòu)成的界面不存在異常入射角，而對某些特殊各向異性巖石界面存在著異常入射角的現(xiàn)象[10]。選擇各向異性頁巖(A頁巖)、泰勒砂巖(T砂巖)以及油頁巖(O頁巖)構(gòu)成兩種各向異性巖石界面模型：以A頁巖作為入射介質(zhì)、T砂巖作為折射介質(zhì)，構(gòu)成第一種各向異性巖石界面模型；以A頁巖作為入射介質(zhì)、O頁巖為折射介質(zhì)，構(gòu)成第二種各向異性巖石界面模型。三種巖石的彈性參數(shù)如表1所示[22]。

表1　三種VTI介質(zhì)巖石的彈性參數(shù)

3.1　模型1：不存在異常入射角情況

考慮不存在異常入射角的各向異性巖石界面系統(tǒng)，重點是對應(yīng)折射P波的過臨界角入射。由于巖石介質(zhì)橫向各向同性，有Wx,2=Wy,2，則模型可簡化為二維。由斯奈爾定理，當(dāng)sinθ2變?yōu)榇笥?的實數(shù)時，式(3)可以轉(zhuǎn)換為質(zhì)點位移分量表達式形式，即折射波質(zhì)點位移的兩個獨立分量可表示為

(6)

(7)

瞬時質(zhì)點位移的表達式為

(8)

由上式可得非均勻折射P波的橢圓極化方程為

(9)

式中

(10)

其中上標(biāo)“*”表示復(fù)共軛。將橢圓極化軌跡末端與x軸之間的夾角稱為非均勻折射P波的橢圓極化旋轉(zhuǎn)角，可表示為

tan(ωt-k0xsinθ0+φ2)]

(11)

上式表明橢圓極化旋轉(zhuǎn)角隨著時間的增大而增大，因此非均勻折射P波是右旋橢圓極化波。

對于給定的各向異性巖石界面，可以得到非均勻折射P波的橢圓極化軌跡與入射角之間的關(guān)系。對于A頁巖和T砂巖組成的界面模型，不存在對應(yīng)折射P波的異常入射角，僅存在一個對應(yīng)折射P波的入射臨界角(θ0C=48.34°)。在超過臨界角區(qū)域，非均勻折射P波是右旋橢圓極化波，計算的非均勻折射P波的橢圓極化狀態(tài)如圖2a所示，可以看出，橢圓極化軌跡的大小相對于入射角的增大而減??；圖2b為不同入射角的橢圓極化軌跡在Wx,2-Wz,2平面的投影；圖2c是橢圓極化初始旋轉(zhuǎn)角與入射角的關(guān)系，其絕對值隨著入射角的增大而單調(diào)遞增，當(dāng)入射角接近90°時初始旋轉(zhuǎn)角趨近于-90°。

3.2　模型2：存在異常入射角情況

對于A頁巖作為入射介質(zhì)、O頁巖作為折射介質(zhì)構(gòu)成界面模型2，存在一個入射臨界角(θ0C=32.99°)和一個異常入射角(θ0A=62.16°)。在前入射角區(qū)域，即θ0∈(θ0C，θ0A)，非均勻折射P波的橢圓極化狀態(tài)的理論分析同模型1，非均勻折射P波的極化狀態(tài)數(shù)學(xué)表達式相同，是一個右旋橢圓極化波。

在過異常入射角區(qū)域θ0∈(θ0A，90°)，由斯奈爾定理可知，當(dāng)cosθ2變?yōu)榇笥?的實數(shù)時，式(3)可以轉(zhuǎn)換為質(zhì)點位移分量表達式形式，即非均勻折射P波質(zhì)點位移的兩個獨立分量可以表示為

(12)

即

(13)

式中因子χ2=iωcosθ2/v2是正數(shù)[12，13]。由上式可得折射P波的橢圓極化方程

(14)

其中

(15)

圖2　界面模型1的非均勻折射P波的極化狀態(tài)

(a)折射P波的橢圓極化軌跡與入射角的關(guān)系； (b)橢圓極化軌跡投影到Wx,2-Wz,2平面上； (c)在坐標(biāo)原點處的橢圓極化的初始旋轉(zhuǎn)角

同理，非均勻折射P波的橢圓極化旋轉(zhuǎn)角可表示為

×

cot(ωt-k0xsinθ0+φ2)}

(16)

由上式可知，橢圓極化旋轉(zhuǎn)角隨著時間增加而減小，橢圓極化軌跡隨著時間沿順時針方向旋轉(zhuǎn)，因此當(dāng)θ0∈(θ0A，90°)，非均勻折射P波是左旋橢圓極化波。

圖3為界面模型2的非均勻折射P波的橢圓極化狀態(tài)。從圖3a可以看出，橢圓極化軌跡的大小隨入射角增加而減小，在異常入射角處，非均勻折射P波的橢圓極化軌跡突然增大，在過異常入射角區(qū)域θ0∈(θ0A，90°)，極化軌跡也是隨著入射角的增大而

圖3　界面模型2的非均勻折射P波的極化狀態(tài)

(a)橢圓極化軌跡的左旋偏振和右旋偏振； (b)橢圓極化軌跡投影到Wx,2-Wz,2平面上； (c)在坐標(biāo)原點處的橢圓極化旋轉(zhuǎn)的初始狀態(tài)減小。圖3b為不同入射角得到的橢圓極化軌跡在Wx,2-Wz,2平面上的投影，可以看出，在小于異常入射角區(qū)域非均勻折射P波為左旋橢圓極化波，在大于異常入射角區(qū)域非均勻折射P波為右旋橢圓極化波。從圖3c可以看出，非均勻折射P波的橢圓極化旋轉(zhuǎn)初始角(絕對值)隨入射角的增大而增大，當(dāng)入射角向異常入射角的方向增大時，橢圓極化旋轉(zhuǎn)初始角接近-90°； 在異常入射角處，橢圓極化旋轉(zhuǎn)初始角(絕對值)突然減小，其旋轉(zhuǎn)方向發(fā)生改變；在過異常入射角區(qū)域θ0∈(θ0A，90°)，橢圓極化旋轉(zhuǎn)初始角相對較小，且變化較為緩慢(圖3c中線段CD)。

4　結(jié)束語

基于兩種真實的VTI介質(zhì)界面模型，即A頁巖與T砂巖之間界面以及A頁巖與O頁巖之間界面，本文推導(dǎo)了兩種界面模型的非均勻折射P波的橢圓極化方程，分析了在上述兩種界面上產(chǎn)生的非均勻折射P波的極化狀態(tài)。在此基礎(chǔ)上，給出了上述兩種界面產(chǎn)生的非均勻折射P波的三維橢圓極化狀態(tài)圖(Poincaré球)。對于上述兩種界面模型，其橢圓極化軌跡隨著入射角的增大而減小。

對于存在異常入射角的VTI介質(zhì)界面(界面模型2)，在異常入射角處，非均勻折射P波橢圓極化不但經(jīng)歷了形狀突變、軌跡銳增和橢圓極化旋轉(zhuǎn)初相角銳減的變化過程，而且由右旋橢圓極化波突變?yōu)樽笮龣E圓極化波。

[1]Crampin S，Stephen R A，McGonigle R.The polarization of P waves in anisotropic media.Geophysical Journal International，1982，68(2)：477-485.

[2]Tsvankin I D，Chesnokov E M.Synthetic of body waves seismograms from point sources in anisotropic media.Journal of Geophysical Research:Solid Earth，1990，95(B7)：11317-11331.

[4]Carcione J M.Wave Fields in Real Media：Wave Pro-pagation in Anisotropic，Anelastic and Porous Media.Elsevier,Amsterdam，2007.

[5]Fa L，Castagna J P，Dong H.An accurately fast algorithm of calculating reflection/transmission coefficients.Science China Physics Mechanics & Astronomy，2008，51(7)：823-846.

[6]Fa L，Castagna J P，Zeng Z et al.Effects of anisotropy on time-depth relation in transversely isotropic me-dium with a vertical axis of symmetry.China Science Bulletin，2010，55(21)：2243-2251.

[7]Fa L，Zhao M，Liu Y et al.Polarization of plane wave propagating inside elastic hexagonal crystal solids.Science China Physics Mechanics & Astronomy，2014，57(2)：251-262.

[8]徐文才，楊國權(quán)，李振春等.橫向各向同性介質(zhì)擬聲波一階速度—應(yīng)力方程.石油地球物理勘探，2016，51(1)：87-96.

Xu Wencai,Yang Guoquan,Li Zhenchun et al.First order velocity-stress equation in TI media.OGP,2016,51(1):87-96.

[9]Fa L，Zhao J，Han Y et al.The influence of rock anisotropy on elliptical-polarization state of inhomogenously refracted P-wave.Science China Physics Mechanics & Astronomy，2016，59(4)：263-272.

[10]法林，聶琳，成洪申.用經(jīng)典定理證明各向異性巖石界面異常入射角的存在.石油地球物理勘探，2017，52(2)：248-252.

Fa Lin,Nie Lin,Cheng Hongshen.Existence of ano-malous incident-angle on anisotropy rock interface verified by classic laws.OGP,2017,52(2):248-252.

[11]張恒磊，胡哲，胡祥云等.基于反射波各向異性特征的保真去噪方法.石油地球物理勘探，2017，52(2)：233-241.

Zhang Henglei,Hu Zhe,Hu Xiangyun et al.Seismic fidelity de-noising with reflection anisotropy.OGP,2017,52(2):233-241.

[12]羊屋三維處理、解釋一體化方法研究組.縱波VTI介質(zhì)理論與應(yīng)用研究.石油地球物理勘探，2002，37(3)：267-275.

Research Group for 3-D Processing and Interpretation Integrated Method.Research on P-wave theory in VTI media and its application.OGP,2002,37(3):267-275.

[13]羊屋三維處理、解釋一體化方法研究組.VTI介質(zhì)的AVO理論與應(yīng)用研究.石油地球物理勘探，2002，37(4)：363-371.

Research Group of Method for Yangwu 3-D Integrating Processing into Interpretation.Study of AVO theory and its application in VTI medium.OGP,2002,37(4):363-371.

[14]劉瑞合，趙金玉，印興耀等.VTI介質(zhì)各向異性參數(shù)層析反演策略與應(yīng)用.石油地球物理勘探，2017，52(3)：484-490.

Liu Ruihe,Zhao Jinyu,Yin Xingyao et al.Strategy of anisotropic parameter tomography inversion in VTI medium.OGP,2017,52(3):484-490.

[15]李勤，李慶春，張林.VTI介質(zhì)多波各向異性參數(shù)分析.石油地球物理勘探，2014，49(3)：503-507.

Li Qin,Li Qingchun,Zhang Lin.Anisotropic parameter analysis of multi-component data in VTI media.OGP,2014,49(3):503-507.

[16]Hosten B.Reflection and transmission of acoustic plane waves on an immersed orthotropic and visco-elastic solid layer.Journal of the Acoustical Society of America，1991，89(6)：2745-2753.

[17]Lanceleur P，Ribeiro H，Belleval J D.The use of inhomogenous waves in the reflection-transmission problem at a plane interface between two anisotropic media.Journal of the Acoustical Society of America，1993，93(4)：1882-1892.

[18]Helbig K，Schoenberg M.Anomalous polarization of elastic waves in transversely isotropic media.Journal of the Acoustical Society of America，1987，81(5)：1235-1245.

[19]Daley P D，Hron F.Reflection and transmission coefficients for transversely isotropic media.Bulletin of the Seismological Society of America，1977，67(3)：661-675.

[20]Fa L，Brown R L，Castagna J P.Anomalous post-critical refraction behavior for certain transversely isotropic media.Journal of the Acoustical Society of America，2006，120(6)：3479-3492.

[21]Auld B A.Acoustic Fields and Waves in Solids.John Wiley and Sons，New York，1973，107-124.

[22]Thomsen L.Weak elastic anisotropy.Geophysics，1986，51(10)：1954-1966.

[23]Wang Z.Seismic anisotropy in sedimentary rocks，part 2：Laboratory data.Geophysics,2002,67(5)：1423-1440.

[24]Vernik L，Nur A.Ultrasonic velocity and anisotropy of hydrocarbon source rocks.Geophysics，1992，57(5)：727-735.