各向异性介质弹性波频散关系在波场分解、速度分析和偏移成像等处理过程中起着至关重要的作用.相比于各向同性介质,拥有更多独立参数的各向异性介质弹性波频散方程显得尤为复杂(杜泽源等,2019;李小斌和阎建国,2021),对应的时空域波动方程也很复杂.倾斜横向各向同性(TTI)介质以贴近实际地球介质情况而被广泛地研究.在观测坐标系下的TTI介质刚度矩阵存在13个不同的弹性参数,而在本构坐标系下实际存在5个独立的弹性参数,这使得描述相、群速度以及频散关系的方程较为复杂(金世勋等,2018).诸如弱各向异性近似(Thomsen, 1986)、小角度近似(Dellinger and Muir, 1988)和椭圆各向异性近似(Cohen, 1997)等近似方法相继被提出,但这些理论都有各自的假设条件.许多学者基于这些近似理论进行了进一步的研究和改进(Dellinger et al., 1993; Alkhalifah and Tsvankin, 1995; Alkhalifah, 1998, 2000; Schoenberg and De Hoop, 2000; Fomel, 2004; Ursin and Stovas, 2006; Fomel and Stovas, 2010; Farra and Pšenčík, 2013; Sripanich and Fomel, 2015).
TTI介质弹性波频散关系方程可以通过求解Christoffel方程得到,利用不同的近似方法对其进行近似处理,并利用傅里叶逆变换将频率-波数域算子变换到时空域,可以得到对应的解耦波动方程. Alkhalifah和Tsvankin(1995)、Alkhalifah(1998, 2000) 开创性地提出了声学假设近似方法,即将沿对称轴的横波速度近似视做零值,推导了一种简化的VTI介质弹性波频散关系方程,进而得到qP波方程.多种纯qP方程相继被提出(Klíe and Toro, 2001; Zhou et al., 2006; Hestholm, 2007; Du et al., 2008, 2010; Duveneck et al., 2008; Chu et al., 2011, 2013; Bloot et al., 2013; Schleicher and Costa, 2016; Li and Zhu, 2018; Xu et al., 2020),但是声学假设条件下的波动方程存在两组共轭解,分别对应qP波和退化qSV波,而后者是不被期望的“干扰”波,在地震波数值模拟的过程中,随着迭代次数的增加,退化qSV波会产生严重的数值频散问题,特别在各向异性参数ε<δ的情况下,描述退化qSV波的解发散,误差随时间累计,导致数值解不稳定.Jin和Stovas(2018)定义一组新的参数用于描述VTI介质中退化qSV波运动学特征,并分析了其传播特征.为了能完全消除声学近似中的qSV波,Klíe和Toro(2001)在此基础上进行改进,消除了描述退化qSV波的一组解析解,得到纯qP波波动方程.此外也可以通过引入辅助变量,实现对VTI介质qP波方程的降阶处理,使其更容易实现(Zhou et al., 2006; Du et al., 2008; Chu et al., 2011).Liu等(2009)基于VTI介质声学近似方程,解耦得到了一种稳定的纯qP波和退化qSV波的波动方程,并应用于逆时偏移中;Fletcher等(2009)基于坐标旋转法推导了一般TTI介质qP波、qSV波的波动方程,并指出若人为设置TTI介质对称轴方向的横波速度足够大,则可以移除qSV的波面三角区,进而使波场能稳定传播;Chu等(2011)对Fletcher推导的TTI介质拟声波方程进行因式分解,利用一阶和高阶Taylor展开近似qP波方程,并指出采用一阶Taylor展开近似的qP波方程足以准确运用于大部分实际应用中;梁锴等(2009)从Thomsen弱各向异性近似和声学假设近似出发,对TTI介质弹性波波动方程进行了分解,导出TTI介质弱各向异性条件近似qP波和qSV波波动方程以及声学假设近似的qP波波动方程;黄翼坚等(2011)对VTI介质弹性波精确相速度表达式进行多项式近似,得到了VTI介质纯qP波方程;杨鹏等(2017)利用近似展开法,通过分解偏微分算子得到了TI介质纯qP波动方程;张庆朝等(2019)根据弱各向异性近似的假设,推导了任意空间取向TI介质弹性波相速度近似表达式,并进一步导出了qP波和qSV波的解耦波动方程,波场模拟结果显示该波动方程有着较好的稳定性;Chu等(2013)利用伪谱法对TTI介质拟声波方程(Chu et al., 2011)进行了数值模拟,误差分析表明一阶Taylor近似适合较弱的各向异性介质,强各向异性介质中需要采用更高阶的近似,以控制误差在允许范围内;杜启振等(2015)采用伪谱法和有限差分法相结合的混合法求解VTI介质纯qP方程,极大地提升了运算效率.
本文在TTI介质弹性波精确频散关系的基础上,利用近似的配方法(梁锴等,2018;孙上饶等,2021),推导TTI介质qP波和qSV波近似频散关系,并将频率-波数域算子反变换到时-空域,导出了qP波和qSV波近似解耦的波动方程.使用两组模型参数对近似频散关系方程进行数值计算,得到了三维频散关系曲面和二维曲线,验证了近似频散关系的有效性.随后考察了近似式与精确式的差异项数值在不同各向异性强度下的变化情况,并分析了XOZ面内近似频散关系曲线的相对误差分布.最后,使用有限差分方法求解TTI介质纯qP波和纯qSV波近似波动方程,模拟了qP波和qSV波在均匀、层状及复杂TTI介质中的传播.数值模拟结果表明,在η<0和各向异性倾角变化较大的介质中,纯qP波和纯qSV波近似波动方程依然可以保持稳定.
1 方法原理 1.1 TTI介质弹性波近似频散关系方程TTI介质弹性波相速度和频散关系可通过将平面波位移方程代入到Christoffel方程求解得到.TTI介质qP波和qSV波精确频散关系方程可表示为:
(1) |
其中,D={[(1+2η)(1+2δ)vP02-vS02](K12+K22)-fvP02K32}2+4fvP02[(1+2δ)vP02-vS02](K12+K22)K32,f=1-vS02/vP02,K1=kxcosθ0+kzsinθ0,K2=ky,K3=-kxsinθ0+kzcosθ0,k为波数,而kx=ksinθcosφ,ky=ksinθsinφ,kz=kcosθ分别表示沿x、y、z方向的波数,θ为传播极角,φ为传播方位角;vP0和vS0分别为沿对称轴方向的qP波和qSV波相速度,ε和δ为Thomsen各向异性参数,
依照近似配方法的思想(梁锴等, 2018)消除D项的一次根号,简化频散关系方程(1),其中式(1)中D项可改写为:
(2) |
令式中a=[(1+2η)(1+2δ)vP02-vS02](K12+K22)+fvP02K32,为了统一量纲,将b构造为:
(3) |
其中
(4) |
令式(4)中右边第二项近似为零,即:
(5) |
那么,D项可近似表示为D≈(a+b)2=Da.将D的近似值Da代入到精确qP波和qSV波频散关系(1)中,即可得到相对应的基于近似配方法的频散关系方程:
(6) |
其中,
经运算知K12+K22+K32=kx2+ky2+kz2.当θ0=0时,TTI介质退化为VTI介质,进而qP波和qSV波近似频散关系方程退化为:
(7) |
其中
TTI介质解耦的波动方程可以通过对TTI介质近似频散关系方程进行傅里叶逆变换得到.将频散关系方程(6)左右两边同时乘以频率-波数域内的qP波波场
(8) |
(9) |
其中PP(t, x, y, z)和PSV(t, x, y, z)分别为qP波和qSV波时空域的波场,AiP和AiSV(i=1~9)为方程系数,具体表达形式见附录A.
解耦波动方程(8)和(9)既包含对时间和空间的混合偏导,也存在混合空间导数,这使得其计算量大于单一变量的空间导数(Fletcher et al., 2009),不利于波动方程的数值求解.考虑波动方程(8)的二维形式,并引入辅助变量QP,变换得到等效的纯qP波波动方程为:
(10) |
(11) |
同理引入辅助变量QSV,对方程(9)进行变换,得到等效qSV波波动方程为:
(12) |
(13) |
可以利用方程(10)—(13)分别进行纯qP波和纯qSV波的数值模拟.
2 数值示例 2.1 数值试算和误差分析Vernik和Liu(1997)对多组含油气地层进行了地震波相速度和各向异性参数的测量,本文从测量数据中选取两组中强各向异性页岩作为TTI介质模型,参数如表 1所示.
通过对比qP波和qSV波的精确和近似频散关系中的D项,定义差异项ΔD为:
(14) |
利用表 1的model 1的速度参数和倾角θ0对ΔD进行数值计算,在|σ|值固定不变的情况下,仅改变参数η和δ的值,结果如图 2所示,图中空心箭头处为沿对称轴方向传播,实心箭头处为垂直对称轴方向传播.参数|η|越小,差异ΔD越小,并且在对称轴传播方向附近的误差相对垂直对称轴传播方向小;在沿对称轴和垂直对称轴的传播方向,近似值与精确值相等.特别地,当η=0时,介质呈现椭圆各向异性,近似频散关系方程(5)等于精确式(1).
对TTI介质qP波和qSV波频散关系进行三维数值计算,结果如图 3所示.图 3a、b为model 1的精确qP波和qSV波频散关系曲面,图 3e、f为model 1的近似qP波和qSV波频散关系曲面,图 3c、d为model 2的精确qP波和qSV波频散关系曲面,图 3g、h为model 2的近似qP波和qSV波频散关系曲面.两组模型中的近似qP波和qSV波频散关系曲面(图 3e—h)与精确频散关系曲面(图 3a—d)方位特征较为相似,其中qP波的精度明显高于qSV波.
为考察近似值与精确值的误差,这里讨论频散曲面在XOZ平面内的频散关系曲线的相对误差分布情况.在观测坐标系XOZ平面内,基于近似配方法的频散关系方程(6)可表示为:
(15) |
定义相对误差(RE)为:
(16) |
根据表 1中两组模型参数可得到二维TTI介质弹性波频散关系曲线,如图 4所示,图中短实线表示TTI介质对称轴方向,qP波和qSV波近似频散关系曲线与精确值在对称轴和垂直对称轴方向误差为零,且误差关于上述两个传播方向呈现对称分布的形式.对比图 4a、b可知,qP波近似频散关系精度整体上高于qSV波,原因在于qSV波近似频散关系的精度在较大程度上受组合参数σ的控制,该参数中纵横波速度比起到放大ε和δ之间差值的作用,从而导致了误差的增大.图 5较为直观地呈现了qP波和qSV波近似频散关系曲线相对误差在θ∈[0, 180°] 范围内的分布情况,本例中沿对称轴(空心箭头处)和垂直对称轴(实心箭头处)的两个传播方向对应的相对误差为零,并且相对误差在这两个传播方向附近一定角度范围内很小,这同时也契合差异项ΔD的变化趋势(图 2).
为验证文中导出的TTI介质qP波和qSV波近似波动方程(10)—(13)的有效性,本文利用有限差分法对其进行均匀、层状和复杂介质波场数值模拟.首先,设置均匀介质模型大小为401×401,空间间隔为Δx=Δz=10 m,时间采样间隔为Δt=1 ms. 模型参数为:vP0=3.00 km·s-1,vS0=1.73 km·s-1,η=0.143,δ=0.2,σ=0.601,θ0=π/4;震源子波是主频为25 Hz的雷克子波,位于模型中心点.高阶空间有限差分近似有利于压制数值频散,但会增加运算成本,合理选择差分近似的阶数可以在适当增加运算时间的情况下得到更好的数值模拟结果.考虑到泊松方程的求解效率问题,本文对方程(10)—(13)取空间十阶、时间二阶精度的差分近似进行波场模拟.图 6展示了TTI介质弹性波以及解耦的qP波和qSV波波场快照.对比图 6a、c可知纯qP波的运动学特征基本与弹性波中qP波的特征相似,即使在大偏移距处,两者的波场特征也很相似;对比图 6b、c可知,相比于弹性波qSV波场,纯qSV波的波前面三角区不明显,但其运动学特征整体保持较好的相似度.
其次,对层状各向异性介质模型进行数值模拟.模型大小为401×401,空间步长为Δx=Δz=10 m,时间采样间隔为Δt=1 ms,模型参数如表 2所示.注意到表 2中layer 1和layer 2为η>0的TTI介质模型,layer 3是η<0的TTI介质模型,这样设置是为了测试纯qP波和纯qSV波波动方程的稳定性.震源子波与均匀介质模型中相同,位于(2 km, 0.1 km)处,检波点设置于深度为0.1 km的水平面上,如图 7a所示;图 7b—d为0.8 s时刻的TTI介质地震波波场快照.通过对比可知,纯qP波(图 7b)和纯qSV波(图 7c)入射到界面时只产生反射、透射的同类波,不产生转换波.相对常规弹性波场(图 7d),前两者有着更加清晰的波场,利于研究同类型波的传播特征.此外,数值模拟结果表明,纯qP波方程在η<0和各向异性倾角变化较大情况下,仍然可以保持较好的稳定性.
最后,本文对BP盐丘模型进行数值模拟;模型大小为601×601,空间间隔和时间采样间隔与上述均匀介质模型一致,模型参数如图 8所示.震源子波与均匀介质模型中相同,位于(3 km, 3 km)处.数值模拟结果表明,纯qP波波场中只可见透射和反射qP波的信息,其波前面与弹性波qP波的波前面位置几乎一致;相应地,纯qSV波中只可见透射和反射qSV波的信息;从相速度误差分析中可见,qSV波的误差大于qP波的误差,因此纯qSV波的波前面与弹性波SV波波前面存在一定的差异,但总体上在可接受的范围内(图 9).
本文在TTI介质弹性波精确频散关系的基础上,利用近似配方法推导了qP波和qSV波近似频散关系,通过傅里叶逆变换将其从频率-波数域变换到时-空域,进而导出了纯qP波和纯qSV波近似解耦波动方程.
(1) 理论分析和数值计算表明,近似和精确频散关系的差异项ΔD与|η|呈正相关,在椭圆各向异性介质中差异项为零,即近似式与精确式等效.
(2) 近似频散关系曲线相对误差整体较小,且在沿对称轴和垂直对称轴的传播方向上相对误差为零.
(3) 参数σ中纵横波速度比使得qSV波近似频散关系曲线相对误差普遍较qP波大.
(4) 数值结果表明,基于近似配方法解耦的纯qP波和纯qSV波波动方程在η<0和对称轴倾角变化较大的介质模型中依然可以保持稳定;相对于弹性波波场,两者均不产生转换波,波场更加清晰.
本文推导的TTI介质解耦波动方程可以被进一步推广到任意空间取向的TI介质中,此时需要从三维空间的角度讨论解耦纵横波的传播特征.本文的解耦波动方程可以用于各向异性介质正演模拟和逆时偏移算法中,但因用常规有限差分方法求解泊松方程的效率较低,可以考虑采用伪谱法、低秩分解法、有限差分法和伪谱法的混合法等效率较高的方法求解纯qP波和纯qSV波方程,以提高正演模拟和偏移成像的效率.
附录A在qP波和qSV波近似频散关系(6)两边同时乘以波场PP(ω, kx, ky, kz)和PSV(ω, kx, ky, kz),再进行傅里叶逆变换(ω→i∂/∂t, kx→i∂/∂x, ky→i∂/∂y, kz→i∂/∂z),化简合并后得到三维TTI介质纯qP波和纯qSV波的近似波动方程:
(A1) |
(A2) |
其中,式(A1)的系数AiP(i=1, 2, …, 9)为:
(A3) |
式(A2)的系数AiSV(i=1, 2, …, 9)为:
(A4) |
Alkhalifah T. 1998. Acoustic approximations for processing in transversely isotropic media. Geophysics, 63(2): 623-631. DOI:10.1190/1.1444361 |
Alkhalifah T. 2000. An acoustic wave equation for anisotropic media. Geophysics, 65(4): 1239-1250. DOI:10.1190/1.1444815 |
Alkhalifah T, Tsvankin I. 1995. Velocity analysis for transversely isotropic media. Geophysics, 60(5): 1550-1566. DOI:10.1190/1.1443888 |
Bloot R, Schleicher J, Santos L T. 2013. On the elastic wave equation in weakly anisotropic VTI media. Geophysical Journal International, 192(3): 1144-1155. DOI:10.1093/gji/ggs066 |
Chu C L, Macy B K, Anno P D. 2011. Approximation of pure acoustic seismic wave propagation in TTI media. Geophysics, 76(5): WB97-WB107. DOI:10.1190/geo2011-0092.1 |
Chu C L, Macy B K, Anno P D. 2013. Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method. Geophysical Prospecting, 61(3): 556-567. DOI:10.1111/j.1365-2478.2012.01077.x |
Cohen J K. 1997. Analytic study of the effective parameters for determination of the NMO velocity function in transversely isotropic media. Geophysics, 62(6): 1855-1866. DOI:10.1190/1.1444286 |
Dellinger J, Muir F. 1988. Imaging reflections in elliptically anisotropic media. Geophysics, 53(12): 1616-1618. DOI:10.1190/1.1442446 |
Dellinger J, Muir F, Karrenbach M. 1993. Anelliptic approximations for TI media. Journal of Seismic Exploration, 2: 23-40. |
Du Q Z, Guo C F, Gong X F. 2015. Hybrid PS/FD numerical simulation and stability analysis of pure P-wave propagation in VTI media. Chinese Journal of Geophysics (in Chinese), 58(4): 1290-1304. DOI:10.6038/cjg20150417 |
Du X, Fletcher R P, Fowler P J. 2008. A new pseudo-acoustic wave equation for VTI media. //78th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts.
|
Du X, Fletcher R P, Fowler P J. 2010. Pure P-wave propagators versus pseudo-acoustic propagators for RTM in VTI media. 80th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts.
|
Du Z Y, Wu G C, LI Y S. 2019. Three-dimensional elastic wave forward modeling in orthorhombic media based on least squares. Progress in Geophysics (in Chinese), 34(1): 69-79. DOI:10.6038/pg2019BB0566 |
Duveneck E, Milcik P, Bakker P M, et al. 2008. Acoustic VTI wave equations and their application for anisotropic reverse-time migration. //78th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts, 2186-2190.
|
Farra V, Pšenčík I. 2013. Moveout approximations for P-and SV waves in VTI media. Geophysics, 78(5): WC81-WC92c. DOI:10.1190/geo2012-0408.1 |
Fletcher R P, Du X, Fowler P J. 2009. Reverse time migration in tilted transversely isotropic (TTI) media. Geophysics, 74(6): WCA179-WCA187. DOI:10.1190/1.3269902 |
Fomel S. 2004. On anelliptic approximations for qP velocities in VTI media. Geophysical Prospecting, 52(3): 247-259. DOI:10.1111/j.1365-2478.2004.00413.x |
Fomel S, Stovas A. 2010. Generalized nonhyperbolic moveout approximation. Geophysics, 75(2): U9-U18. DOI:10.1190/1.3334323 |
Hestholm S. 2007. Acoustic VTI modeling using high-order finite-differences. //77th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts, 139-143.
|
Huang Y J, Zhu G M, Liu C Y. 2011. An approximate acoustic wave equation for VTI media. Chinese Journal of Geophysics (in Chinese), 54(8): 2117-2123. DOI:10.3969/j.issn.0001-5733.2011.08.019 |
Jin S, Stovas A. 2018. S-wave kinematics in acoustic transversely isotropic media with a vertical symmetry axis. Geophysical Prospecting, 66(6): 1123-1137. DOI:10.1111/1365-2478.12635 |
Kim S H, Yin X Y, Liang K. 2018. Decomposition of stiffness matrix for TTI media. Progress in Geophysics (in Chinese), 33(2): 644-652. DOI:10.6038/pg2018BB0240 |
Klíe H, Toro W. 2001. A new acoustic wave equation for modeling in anisotropic media. //71st Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts, 1171-1174.
|
Li X B, Yan J G. 2021. Seismicwavefield simulation and azimuth feature analysis for fractured media based on Bond transformation. Progress in Geophysics (in Chinese), 36(5): 2034-2044. DOI:10.6038/pg2021EE0401 |
Li X Y, Zhu H J. 2018. A finite-difference approach for solving pure quasi-P-wave equations in transversely isotropic and orthorhombic media. Geophysics, 83(4): C161-C172. DOI:10.1190/geo2017-0405.1 |
Liang K, Jin S X, Yin X Y. 2018. General approximation of phase velocity for elastic waves in TTI media. Geophysical Prospecting for Petroleum (in Chinese), 57(6): 813-822. |
Liang K, Wu G C, Yin X Y, et al. 2009. Wave equation decomposition in 3-D TTI medium. Oil Geophysical Prospecting (in Chinese), 44(1): 19-27. |
Liu F, Scott A M, Jiang S, et al. 2009. Decoupled wave equations for P and SV waves in an acoustic VTI media. //71nd Annual Conference and Exhibition, EAGE, Extended Abstracts, 2844-2848.
|
Schleicher J, Costa J C. 2016. A separable strong-anisotropy approximation for pure qP-wave propagation in transversely isotropic media. Geophysics, 81(6): C337-C354. DOI:10.1190/geo2016-0138.1 |
Schoenberg M A, De Hoop M V. 2000. Approximate dispersion relations for qP-qSV-waves in transversely isotropic media. Geophysics, 65(3): 919-933. DOI:10.1190/1.1444788 |
Sripanich Y, Fomel S. 2015. On anelliptic approximations for qP velocities in transversely isotropic and orthorhombic media. Geophysics, 80(5): C89-C105. DOI:10.1190/geo2014-0534.1 |
Sun S R, Liang K, Yin X Y, et al. 2021. Approximate 3D phase and group velocity for elastic wave in TTI media based on an approximate match method. Oil Geophysical Prospecting (in Chinese), 56(3): 496-504. |
Thomsen L. 1986. Weak elastic anisotropy. Geophysics, 51(10): 1954-1966. DOI:10.1190/1.1442051 |
Tsvankin I, Thomsen L. 1994. Nonhyperbolic reflection moveout in anisotropic media. Geophysics, 59(8): 1290-1304. DOI:10.1190/1.1443686 |
Ursin B, Stovas A. 2006. Traveltime approximations for a layered transversely isotropic medium. Geophysics, 71(2): D23-D33. DOI:10.1190/1.2187716 |
Vernik L, Liu X Z. 1997. Velocity anisotropy in shales: A petrophysical study. Geophysics, 62(2): 521-532. DOI:10.1190/1.1444162 |
Xu S B, Stovas A, Alkhalifah T, et al. 2020. New acoustic approximation for transversely isotropic media with a vertical symmetry axis. Geophysics, 85(1): C1-C12. DOI:10.1190/geo2019-0100.1 |
Yang P, Li Z C, Gu B L. 2017. Pure quasi-P wave forward modeling method in TI media and its application to RTM. Chinese Journal of Geophysics (in Chinese), 60(11): 4447-4467. DOI:10.6038/cjg20171130 |
Zhang Q C, Zhu G W, He D K, et al. 2019. A decoupled qP-wave and qSV-wave equation in TI media with arbitrary spatial oriented symmetry axis. Chinese Journal of Geophysics (in Chinese), 62(11): 4353-4366. DOI:10.6038/cjg2019M0528 |
Zhou H, Zhang G, Bloor R. 2006. An anisotropic acoustic wave equation for VTI media. //76th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts.
|
杜启振, 郭成锋, 公绪飞. 2015. VTI介质纯P波混合法正演模拟及稳定性分析. 地球物理学报, 58(4): 1290-1304. DOI:10.6038/cjg20150417 |
杜泽源, 吴国忱, 李雨生. 2019. 基于最小二乘的三维正交介质弹性波高精度正演模拟. 地球物理学进展, 34(1): 69-79. DOI:10.6038/pg2019BB0566 |
黄翼坚, 朱光明, 刘池洋. 2011. 一个近似的VTI介质声波方程. 地球物理学报, 54(8): 2117-2123. DOI:10.3969/j.issn.0001-5733.2011.08.019 |
金世勋, 印兴耀, 梁锴. 2018. TTI介质刚度矩阵的分解. 地球物理学进展, 33(2): 644-652. DOI:10.6038/pg2018BB0240 |
李小斌, 阎建国. 2021. 基于Bond变换的裂缝介质波场方位特征模拟及分析. 地球物理学进展, 36(5): 2034-2044. DOI:10.6038/pg2021EE0401 |
梁锴, 金世勋, 印兴耀. 2018. TTI介质弹性波相速度的一般性近似. 石油物探, 57(6): 813-822. DOI:10.3969/j.issn.1000-1441.2018.06.003 |
梁锴, 吴国忱, 印兴耀, 等. 2009. 三维TTI介质波动方程分解. 石油地球物理勘探, 44(1): 19-27. DOI:10.3321/j.issn:1000-7210.2009.01.005 |
孙上饶, 梁锴, 印兴耀, 等. 2021. 三维TTI介质弹性波相、群速度的近似配方表征. 石油地球物理勘探, 56(3): 496-504. |
杨鹏, 李振春, 谷丙洛. 2017. 一种TI介质纯qP波正演方法及其在逆时偏移中的应用. 地球物理学报, 60(11): 4447-4467. DOI:10.6038/cjg20171130 |
张庆朝, 朱国维, 何登科, 等. 2019. 任意空间取向TI介质qP和qSV波解耦的波动方程. 地球物理学报, 62(11): 4353-4366. DOI:10.6038/cjg2019M0528 |