Calculation method for ray velocity in the general TI medium using Thomsen's parameter model

ZiKang XU, GuangNan HUANG, DaYu ZHANG, HongXing LI, ZhenWei ZHONG

Prog Geophy ›› 2026, Vol. 41 ›› Issue (2) : 823-833.

PDF(1670 KB)
Home Journals Progress in Geophysics
Progress in Geophysics

Abbreviation (ISO4): Prog Geophy      Editor in chief:

About  /  Aim & scope  /  Editorial board  /  Indexed  /  Contact  / 
PDF(1670 KB)
Prog Geophy ›› 2026, Vol. 41 ›› Issue (2) : 823-833. DOI: 10.6038/pg2026JJ0174

Calculation method for ray velocity in the general TI medium using Thomsen's parameter model

Author information +
History +

Abstract

The group velocity in a general Transversely Isotropic (TI) medium is a function of the phase slowness direction angle. However, in practical calculations, it is often necessary to compute the group velocity along a specific ray direction. Thus, the corresponding phase slowness direction angle must first be derived from the ray direction, and then the group velocity in that direction can be calculated. Traditional ray velocity calculation methods for general TI media typically discretize the phase slowness angle with small angular intervals, compute the ray direction angle and ray velocity for each discrete angle, and then use interpolation based on the ray direction angle-ray velocity database to obtain the ray velocity for the actual ray direction. Nevertheless, this approach suffers from time-consuming computations and insufficient accuracy. This paper proposes a calculation method for ray velocity in the context of the general TI medium Thomsen parameter model. Firstly, the phase slowness vector is represented using the phase slowness direction vector and phase velocity. Subsequently, a system of equations is directly established between the ray direction and the phase slowness direction angle. Finally, the phase slowness direction angle obtained from solving the system of equations is substituted into the group velocity formula to derive the ray velocity for that ray direction. In the numerical simulations section, this ray velocity calculation method was used to obtain correct traveltime isosurfaces for homogeneous Vertical Transverse Isotropy (VTI) and Tilted Transverse Isotropy (TTI) medium models. Additionally, when applying the ray velocity formula to the traveltime calculation of an overthrust TTI model, the traveltime results were identical to those obtained from an elastic parameter model. Furthermore, by introducing an interpolation method to solve the phase slowness direction angle along the ray direction and applying it to the shortest path ray tracing method, the corresponding seismic wave traveltime field was obtained, providing a reliable comparison for the traveltime results of qP-waves and qSV-waves in the overthrust TTI model. Comparative tests show that the ray velocity calculation method for the Thomsen parameter model of general TI media proposed in this paper is correct, and compared with traditional interpolation methods, this method significantly improves the computational efficiency.

Key words

General TI Media / Thomsen Model / Ray Velocity / Phase Slowness

Cite this article

Download Citations
ZiKang XU , GuangNan HUANG , DaYu ZHANG , et al . Calculation method for ray velocity in the general TI medium using Thomsen's parameter model[J]. Progress in Geophysics. 2026, 41(2): 823-833 https://doi.org/10.6038/pg2026JJ0174

References

Brown R J , Lawton D C , Cheadle S P . Scaled physical modelling of anisotropic wave propagation: multioffset profiles over an orthorhombic medium. Geophysical Journal International, 1991, 107 (3): 693- 702.
Every A G , Kim K Y . Time domain dynamic response functions of elastically anisotropic solids. The Journal of the Acoustical Society of America, 1994, 95 (5): 2505- 2516.
Farra V . Improved first-order approximation of group velocities in weakly anisotropic media. Studia Geophysica et Geodaetica, 2004, 48 (1): 199- 213.
Fomel S . On anelliptic approximations for qP velocities in VTI media. Geophysical Prospecting, 2004, 52 (3): 247- 259.
Grechka V . Ray-direction velocities in VTI media. Geophysics, 2013, 78 (1): F1- F5.
Guo L P , Yang Q Y , Li Z C , et al. First arrival ray tracing in complex anisotropic medium. Geophysical Prospecting for Petroleum, 2016, 55 (1): 18- 24.
Helbig K , Thomsen L . 75th Anniversary Paper: 75-plus years of anisotropy in exploration and reservoir seismics: A historical review of concepts and methods. Geophysics, 2005, 70 (6): 9ND- 23ND.
Huang G N , Liu Y , Li H X , et al. Research on sensitivity of phase velocity with respect to elastic moduli parameters in TTI media. Geophysical Prospecting for Petroleum, 2014, 53 (3): 262- 271.
Huang G N , Zhou B , Deng J Z , et al. Traveltime tomography of qP reflection waves in anisotropic TI media. Chinese Journal of Geophysics, 2015a, 58 (6): 2035- 2045.
Huang G N , Zhou B , Deng J Z , et al. Traveltime tomography of qP reflection waves in anisotropic TI media. Chinese Journal of Geophysics, 2015, 58 (3): 302- 312.
Huang G N , Deng J Z , Li H X , et al. Nonlinear traveltime tomography of P-SV and P-SH waves. Oil Geophysical Prospecting, 2015b, 50 (6): 1127- 1133.
Huang G N , Luo S T , Deng J Z , et al. Traveltime calculations for qP, qSV, and qSH waves in two-dimensional tilted transversely isotropic media. Journal of Geophysical Research: Solid Earth, 2020, 125 (8): e2019JB018868
Jiang X Y , Zhang Y , Gan L D , et al. Seismic techniques for predicting fractures in granite buried hills. Oil Geophysical Prospecting, 2020, 55 (3): 694- 704.
Li L. 2008. Analytical expressions of phase and group velocities for qP waves propagating in strong anisotropic media of arbitrary symmetry[Master's thesis](in Chinese). Beijing: Institute of Geology, China Earthquake Administration.
Li L , Zhang G Z , Zhang J J , et al. Fracture weakness inversion based on reflection amplitude Fourier coefficient. Journal of China University of Petroleum (Edition of Natural Science), 2020, 44 (6): 36- 45.
Liang K , Cao D P , Sun S R , et al. Extended anisotropic linear approximation for seismic wave group velocity in 3D VTI media. Journal of China University of Petroleum (Edition of Natural Science), 2021, 45 (6): 51- 59.
Lin Y , Guo X Y , Sun X E , et al. Studies on the influence of anisotropy in seismic exploration. Oil Geophysical Prospecting, 2010, 45 (4): 606- 623.
Liu Y M , Yin X Y , Zhang S Y , et al. Recent advances in wide-azimuth seismic exploration. Oil Geophysical Prospecting, 2014, 49 (3): 596- 610.
Ma D T , Zhu G M , Fan T E . A search algorithm for first break traveltime in 2D TTI medium. Oil Geophysical Prospecting, 2011, 46 (5): 710- 714. 710-714, 836
Meng Y J , Shao J B , Li S P , et al. Application of geophysical methods to the exploration of nickel ore in Ngara area, Tanzania. Progress in Geophysics, 2025, 40 (3): 1004- 1015.
Pšenčík I , Vavryčuk V . Approximate relation between the ray vector and the wave normal in weakly anisotropic media. Studia Geophysica et Geodaetica, 2002, 46 (4): 793- 807.
Qi H , Huang G N , Chen H , et al. Research on source singularity correction for finite difference travel time calculation with shortest path method. Progress in Geophysics, 2025, 40 (3): 1105- 1114.
Song L P , Every A G . Approximate formulae for acoustic wave group slownesses in weakly orthorhombic media. Journal of Physics D: Applied Physics, 2000, 33 (17): L81
Tai L X , Gao Y , Liu G , et al. Crustal seismic anisotropy in the southeastern margin of Tibetan Plateau by ChinArray data: shear-wave splitting from temporary observations of the first phase. Chinese Journal of Geophysics, 2015, 58 (11): 4079- 4091.
Tan S J , Yu C Q , Nie F J , et al. Exploring the relationship between shallow sandstone-type uranium deposits and deep oil and gas in the Pengyang area of the Ordos Basin based on seismic data. Progress in Geophysics, 2025, 40 (2): 524- 540.
Teng J W , Wang G J , Zhang B M , et al. The seismic anisotropy and geodynamics of Earth's interior media. Progress in Geophysics, 2000, 15 (1): 1- 35.
Thomsen L . Weak elastic anisotropy. Geophysics, 1986, 51 (10): 1954- 1966.
Vavryčuk V . Elastodynamic and elastostatic Green tensors for homogeneous weak transversely isotropic media. Geophysical Journal International, 1997, 130 (3): 786- 800.
Vavryčuk V . Calculation of the slowness vector from the ray vector in anisotropic media. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2006, 462 (2067): 883- 896.
Wang Y , Yang D H , Yin C C , et al. Anisotropic geophysics and vector field. Chinese Science Bulletin, 2017, 62 (23): 2595- 2605.
Yin X Y , Liu Q . Anisotropic rock physics modeling of tight sandstone and applications. Journal of China University of Petroleum (Edition of Natural Science), 2016, 40 (2): 52- 58.
Zhang L , Meng Y , Qi L S , et al. Research progress of fault identification technology based on seismic data. Progress in Geophysics, 2025, 40 (1): 208- 219.
Zhang L J , Zhou B . Calculation of slowness vectors from ray directions for qP-, qSV-, and qSH-waves in tilted transversely isotropic media. Geophysics, 2018, 83 (4): C153- C160.
Zhang Z J . A review of the seismic anisotropy and its applications. Progress in Geophysics, 2002, 17 (2): 281- 293.
Zhao A H , Ding Z F . New approximate expressions of Seismic group velocities for weakly anisotropic media. Progress in Geophysics, 2005, 20 (4): 916- 919.
立鹏 , 勤勇 , 振春 , 等. 复杂各向异性介质初至波射线追踪. 石油物探, 2016, 55 (1): 18- 24.
光南 , , 红星 , 等. TTI介质相速度对弹性模量参数的敏感性研究. 石油物探, 2014, 53 (3): 262- 271.
光南 , Zhou B , 居智 , 等. 各向异性TI介质qP反射波走时层析成像. 地球物理学报, 2015a, 58 (6): 2035- 2045.
光南 , 居智 , 红星 , 等. P-SV波和P-SH波非线性走时层析成像. 石油地球物理勘探, 2015b, 50 (6): 1127- 1133.
晓宇 , , 利灯 , 等. 花岗岩潜山裂缝地震预测技术. 石油地球物理勘探, 2020, 55 (3): 694- 704.
李磊. 2008. 任意对称系强各向异性介质准纵波相速度和群速度解析表达式[硕士论文]. 北京: 中国地震局地质研究所.
, 广智 , 佳佳 , 等. 基于反射振幅傅里叶系数的裂缝弱度反演方法. 中国石油大学学报(自然科学版), 2020, 44 (6): 36- 45.
, 丹平 , 上饶 , 等. VTI介质地震波群速度三维扩展各向异性线性近似表征. 中国石油大学学报(自然科学版), 2021, 45 (6): 51- 59.
, 向宇 , 祥娥 , 等. 地震勘探中的各向异性影响问题研究. 石油地球物理勘探, 2010, 45 (4): 606- 623.
依谋 , 兴耀 , 三元 , 等. 宽方位地震勘探技术新进展. 石油地球物理勘探, 2014, 49 (3): 596- 610.
德堂 , 光明 , 廷恩 . 二维TTI介质中初至波走时的搜索算法. 石油地球物理勘探, 2011, 46 (5): 710- 714. 710-714, 836
有杰 , 江波 , 水平 , 等. 地球物理方法在坦桑尼亚恩加拉地区镍矿勘查中的应用研究. 地球物理学进展, 2025, 40 (3): 1004- 1015.
, 光南 , , 等. 结合最短路径法的有限差分走时计算源点奇异性校正研究. 地球物理学进展, 2025, 40 (3): 1105- 1114.
龄雪 , , , 等. 利用中国地震科学台阵研究青藏高原东南缘地壳各向异性: 第一期观测资料的剪切波分裂特征. 地球物理学报, 2015, 58 (11): 4079- 4091.
顺佳 , 常青 , 逢君 , 等. 应用地震资料探讨鄂尔多斯盆地彭阳地区浅层砂岩型铀矿与深层油气关系. 地球物理学进展, 2025, 40 (2): 524- 540.
吉文 , 光杰 , 秉铭 , 等. 地球内部各圈层介质的地震各向异性与地球动力学. 地球物理学进展, 2000, 15 (1): 1- 35.
, 顶辉 , 长春 , 等. 各向异性地球物理与矢量场. 科学通报, 2017, 62 (23): 2595- 2605.
兴耀 , . 致密储层各向异性地震岩石物理建模及应用. 中国石油大学学报(自然科学版), 2016, 40 (2): 52- 58.
, , 丽莎 , 等. 利用地震资料进行断层识别技术的研究进展. 地球物理学进展, 2025, 40 (1): 208- 219.
中杰 . 地震各向异性研究进展. 地球物理学进展, 2002, 17 (2): 281- 293.
爱华 , 志峰 . 一种弱各向异性介质地震波群速度的近似表示新方法. 地球物理学进展, 2005, 20 (4): 916- 919.

感谢审稿专家提出的修改意见和编辑部的大力支持!

RIGHTS & PERMISSIONS

Copyright ©2026 Progress in Geophysics. All rights reserved.
PDF(1670 KB)

Accesses

Citation

Detail

Sections
Recommended

/