Multi-phase seismic traveltime computation in 2D/3D undulated surface model using factored eikonal equation
ZHANG Yun1, LI Xihai1, BAI Chaoying2,3, NIU Chao1, WANG Yiting1, ZENG Xiaoniu1
1. Rocket Force University of Engineering, Xi'an, Shaanxi 710025, China; 2. Department of Geophysics, College of Geological Engineering and Geomatics, Chang'an University, Xi'an, Shaanxi 710054, China; 3. Institute of Computational Geophysics, Chang'an University, Xi'an, Shaanxi 710054, China
Abstract:The seismic traveltime computation scheme in undulating surface conditions is a basic tool to study the underground structures of such surface areas. The fast marching method (FMM) and the fast sweeping method (FSM) are both developed based on solving the eikonal equation with finite difference. They have the problem of source singularity due to the high curvature of the wavefront around the source. Previous studies show that the com-putational efficiency of FMM is higher than FSM for complex models. Thus, this paper employs the FMM to solve the factorization equation and avoid source singularity. Specifically, the original eikonal equation can be transformed into the factored eikonal equation, in which the seismic traveltime can be regarded as the product of a distance function T 0 and a correction factor T1 of traveltime. The correction factor T1 of traveltime can be solved by the FMM algorithm and then the distance function T0 is multiplied to obtain the traveltime(T). To address the problem that the upwind finite difference formula with even grid spacing is not applicable to surface/interface undulation, this paper constructs an upwind finite difference formula with uneven grid spacing. Finally, the multistage computational technique is adopted to propose a computation method for global multi-phase seismic traveltime. The simulation tests indicate that the new algorithm solves the source singularity of FMM, significantly improves the computational accuracy and efficiency of the original algorithm, and can accurately calculate the multi-phase seismic traveltime.
RAWLINSON N, KENNETT B L. Teleseismic tomography of the upper mantle beneath the southern Lachlan Orogen, Australia[J]. Physics of the Earth and Planetary Interiors, 2008, 167(1-2):84-97.
[2]
张风雪, 李永华, 吴庆举, 等. FMTT方法研究华北及邻区上地幔P波速度结构[J]. 地球物理学报, 2011, 54(5):1233-1242.ZHANG Fengxue, LI Yonghua, WU Qingju, et al. The P wave velocity structure of upper mantle beneath the North China and surrounding regions from FMTT[J]. Chinese Journal of Geophysics, 2011, 54(5):1233-1242.
[3]
左建军, 林松辉, 孔庆丰, 等. 井间地震直达波和反射波联合层析成像及应用[J]. 石油地球物理勘探, 2011, 46(2):226-231.ZUO Jianjun, LIN Songhui, KONG Qingfeng, et al. Imaging by a joint tomography of direct wave traveltime and reflection wave traveltime in crosswell seismic and its application[J]. Oil Geophysical Prospecting, 2011, 46(2):226-231.
BAI C, GREENHALGH S. 3D local earthquake hypocenter determination with an irregular shortest-path method[J]. Bulletin of the Seismological Society of America, 2006, 96(6):2257-2268.
[6]
赵爱华, 丁志峰, 白志明. 基于射线追踪技术计算地震定位中震源轨迹的改进方法[J]. 地球物理学报, 2015, 58(9):3272-3285.ZHAO Aihua, DING Zhifeng, BAI Zhiming. Improvement of the ray-tracing based method calculating hypocentralloci for earthquake location[J]. Chinese Journal of Geophysics, 2015, 58(9):3272-3285.
[7]
刘磊, 宋维琪, 杨小慧, 等. 剪张源约束下的单井微地震震源机制反演[J]. 石油地球物理勘探, 2021, 56(5):1093-1104.LIU Lei, SONG Weiqi, YANG Xiaohui, et al. Single-well microseismic focal mechanism inversion under shear-tensile source constraint[J]. Oil Geophysical Prospecting, 2021, 56(5):1093-1104.
[8]
JULIAN B R, GUBBINS D. Three dimensional seismic ray tracing[J]. Journal of Geophysics, 1977, 43(1):95-113.
[9]
UM J, THURBER C. A fast algorithm for two-point seismic ray tracing[J]. Bulletin of the Seismological Society of America, 1987, 77(3):972-986.
[10]
MOSER T J. Shortest path calculation of seismic rays[J]. Geophysics, 1991, 56(1):59-67.
[11]
BAI C, GREENHALGH S, ZHOU B. 3D ray tracing with a modified shortest path method[J]. Geophysics, 2007, 72(4):T27-T36.
[12]
LEUNG S, QIAN J. An adjoint state method for threedimensional transmission traveltime tomography using first-arrivals[J]. Communications in Mathematical Sciences, 2006, 4(1):249-266.
[13]
TONG P. Adjoint-state traveltime tomography:eikonal equation-based methods and application to the Anza Area in Southern California[J]. Journal of Geophysics Research:Solid Earth, 2021, 126(5):e2021JB021818.
[14]
谢春, 刘玉柱, 董良国, 等. 伴随状态法初至波走时层析[J]. 石油地球物理勘探, 2014, 49(5):877-883.XIE Chun, LIU Yuzhu, DONG Liangguo, et al. First arrival tomography based on the adjoint state method[J]. Oil Geophysical Prospecting, 2014, 49(5):877-883.
[15]
LUO S, QIAN J. Fast sweeping methods for factored anisotropic eikonal equations:multiplicative and additive factors[J]. Journal of Scientific Computing, 2012, 52(2):360-382.
[16]
HAN S, ZHANG W, ZHANG J. Calculating qP-wave traveltimes in 2-D TTI media by high-order fast sweeping methods with a numerical quartic equation solver[J]. Geophysical Journal International, 2017, 210(3):1560-1569.
[17]
崔宁城, 黄光南, 李红星, 等. 求解椭圆各向异性介质程函方程的源点快速扫描算法[J]. 石油地球物理勘探, 2019, 54(4):796-804.CUI Ningcheng, HUANG Guangnan, LI Hongxing, et al. A source fast sweeping method for solving elliptically anisotropic eikonal equation[J]. Oil Geophysical Prospecting, 2019, 54(4):796-804.
[18]
HUANG G, LUO S, DENG J, et al. Traveltime calculations for qP, qSV, and qSH waves in two-dimensional tilted transversely isotropic media[J]. Journal of Geophysical Research:Solid Earth, 2020, 125(8):e2019 JB018868.
[19]
HUANG G, LUO S. Hybrid fast sweeping methods for anisotropic eikonal equation in two-dimensional tilted transversely isotropic media[J]. Journal of Scientific Computing, 2020, 84(2):32.
[20]
芦永明, 张伟. 求解二次方程的二维VTI介质qSV波和qSH波走时快速扫描算法[J]. 地球物理学报, 2022, 65(1):244-255.LU Yongming, ZHANG Wei. Calculating the traveltimes of qSV and qSH waves using fast sweeping method for 2D VTI media by solving quadratic equation[J]. Chinese Journal of Geophysics, 2022, 65(1):244-255.
[21]
SETHIAN J A, POPOVICI A M. 3-D traveltime computation using the fast marching method[J]. Geophysics, 1999, 64(2):516-523.
[22]
ZHAO H. A fast sweeping method for Eikonal equations[J]. Mathematics of Computation, 2005, 74(250):603-627.
[23]
SETHIAN J A, VLADIMIRSKY A. Ordered upwind methods for static Hamilton-Jacobi equations:theory and algorithms[J]. SIAM Journal on Numerical Analysis, 2003, 41(1), 325-363
[24]
QIAO B, PAN Z, HUANG W, et al. An adaptive finitedifference method for accurate simulation of first-arrival traveltimes in heterogeneous media[J]. Applied Mathematics and Computation, 2021, 394:125792.
[25]
WAHEED U B. A fast marching eikonal solver for tilted transversely isotropic media[J]. Geophysics, 2020, 85(6):S385-S393.
[26]
兰海强, 张智, 徐涛, 等. 地震波走时场模拟的快速推进法和快速扫描法比较研究[J]. 地球物理学进展, 2012, 27(5):1863-1870.LAN Haiqiang, ZHANG Zhi, XU Tao, et al. A comparative study on the fast marching and fast sweeping methods in the calculation of first-arrival traveltime field[J]. Progress in Geophysics, 2012, 27(5):1863-1870.
[27]
GREMAUD P A, KUSTER C M. Computational study of fast methods for the eikonal equation[J]. SIAM Journal on Scientific Computing, 2006, 27(6):1803-1816.
[28]
CHACON A, VLADIMIRSKY A. Fast two-scale methods for eikonal equations[J]. SIAM Journal on Scientific Computing, 2012, 34(2):A547-A578.
[29]
LI S, VLADIMIRSKY A, FOMEL S. First-break traveltime tomography with the double-square-root eikonal equation[J]. Geophysics, 2013, 78(6):U89-U101.
[30]
BENAICHOUCHE A, NOBEL M. GESRET A. First arrival traveltime tomography using the fast marching method and the adjoint state technique[C]. Extended Abstracts of 77th EAGE Conference and Exhibition, 2015, 1-5.
[31]
RAWLINSON N, SAMBRIDGE M. Multiple reflection and transmission phases in complex layered media using multistage fast marching method[J]. Geophysics, 2004, 69(5):1338-1350.
[32]
DE KOOL M, RAWLINSON N, SAMBRIDGE M. A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media[J]. Geophysical Journal International, 2006, 167(1):253-270.
[33]
孙章庆, 孙建国, 韩复兴. 三维起伏地表条件下地震波走时计算的不等距迎风差分法[J]. 地球物理学报, 2012, 55(7):2441-2449.SUN Zhanqing, SUN Jianguo, HAN Fuxing. Traveltime computation using the upwind finite difference method with nonuniform grid spacing in a 3D undulating surface condition[J]. Chinese Journal of Geophysics, 2012, 55(7):2441-2449.
[34]
KIM S, FOLIE D. The group marching method:An O(N) level set eikonal solver[C]. SEG Technical Program Expanded Abstracts, 2000, 19:2297-2300.
[35]
孙章庆, 孙建国, 王雪秋, 等. 三维复杂山地多级次群推进迎风混合法多波型走时计算[J]. 地球物理学报, 2017, 60(5):1861-1873.SUN Zhangqing, SUN Jianguo, WANG Xueqiu, et al. Computation of multiple seismic traveltime in mountainous areas with complex 3D conditions using the multistage group marching upwind hybrid method[J]. Chinese Journal of Geophysics, 2017, 60(5):1861-1873.
[36]
ALKHALIFAH T, FOMEL S. Implementing the fast marching eikonal solver:Spherical versus Cartesian coordinates[J]. Geophysical Prospecting, 2001, 49(2):165-178.
[37]
FOMEL S, LUO S, ZHAO H. Fast sweeping method for the factored eikonal equation[J]. Journal of Computational Physics, 2009, 228(17):6440-6455.
[38]
周小乐, 兰海强, 陈凌, 等. 曲线坐标系因式分解程函方程及其走时计算[J]. 地球物理学报, 2020, 63(2):638-651.ZHOU Xiaole, LAN Haiqiang, CHEN Ling, et al. The factored eikonal equation in curvilinear coordinate system and its numerical solution[J]. Chinese Journal of Geophysics, 2020, 63(2):638-651.
[39]
BAI C, HUANG G, ZHAO R. 2-D/3-D irregular shortest-path ray tracing for multiple arrivals and its applications[J]. Geophysical Journal International, 2010, 183(3):1596-1612.