A source fast sweeping method for solving elliptically anisotropic eikonal equation
CUI Ningcheng1,2, HUANG Guangnan1,2, LI Hongxing1,2, ZHANG Hua1,2, ZHANG Xiaofeng1,2, XIAO Kun1,2
1. State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang, Jiangxi 330013, China;
2. Department of Geophysics, School of Geophysics and Measurement-control Technology, East China University of Technology, Nanchang, Jiangxi 330013, China
Abstract:The seismic traveltime calculation is an important part of ray tracing,tomography,and seismic migration.Conventional traveltime fast sweeping methods scan directly the global computational domain without considering the scanning validity.In a real calculation,the scanning without the source point is invalid due to the lack of source information and cannot calculate right traveltime.In this paper,a source fast sweeping method is proposed,which moves the starting point of scanning to the source without invalid scanning,and improves the calculation efficiency.The source fast sweeping method in elliptical anisotropy medium is realized based on the factored eikonal equations and upwind difference scheme.Numerical simulation results show that the proposed source fast sweeping method obtains the same results as conventional methods,but it achieves higher calculation efficiency.
Stoep D M.Velocity anisotropy measurements in wells[J].Geophysics,1966,31(5):900-916.
[2]
Crampin S.Suggestions for a consistent terminology for seismic anisotropy[J].Geophysical Prospecting,1989,37(7):753-770.
[3]
Fomel S.On anelliptic approximations for qP velocities in VTI media[J].Geophysical Prospecting,2004,52(3):247-259.
[4]
Alkhalifah T,Fomel S.The basic components of residual migration in VTI media using anisotropy continuation[J].Journal of Petroleum Exploration and Production Technology,2011,1(1):17-22.
[5]
黄国娇,孙江兵,白超英,等.三维TI介质中多波旅行时层析成像[J].石油地球物理勘探,2018,53(1):63-72.HUANG Guojiao,SUN Jiangbing,BAI Chaoying,et al.Seismic multi-wave traveltime tomography in 3D TI media[J].Oil Geophysical Prospecting,2018,53(1):63-72.
[6]
刘财,邓馨卉,郭智奇,等.基于岩石物理的页岩储层各向异性表征[J].石油地球物理勘探,2018,53(2):339-346.LIU Cai,DENG Xinhui,GUO Zhiqi,et al.Shale reservoir anisotropic characterization based on rock physics[J].Oil Geophysical Prospecting,2018,53(2):339-346.
程玖兵,康玮,王腾飞.各向异性介质qP波传播描述Ⅰ:伪纯模式波动方程[J].地球物理学报,2013,56(10):3474-3486.CHENG Jiubing,KANG Wei,WANG Tengfei.Description of qP-wave propagation in anisotropic media,Part Ⅰ:Pseudo-pure-mode wave equations[J].Chinese Journal of Geophysics,2013,56(10):3474-3486.
[9]
程玖兵,陈茂根,王腾飞,等.各向异性介质qP波传播描述Ⅱ:分离纯模式标量波[J].地球物理学报,2014,57(10):3389-3401.CHENG Jiubing,CHEN Maogen,WANG Tengfei,et al.Description of qP-wave propagation in anisotropic media,Part Ⅱ:Separation of pure-mode scalar waves[J].Chinese Journal of Geophysics,2014,57(10):3389-3401.
[10]
Huang G,Zhou B,Li H,et al.2D seismic reflection tomography in strongly anisotropic media[J].Journal of Geophysics and Engineering,2014,11(6):1-8.
[11]
Huang G,Zhou B,Li H,et al.Seismic traveltime inversion based on tomographic equation without integral terms[J].Computers & Geosciences,2017,104:29-34.
[12]
Taillandier C,Noble M,Chauris H,et al.First-arrival traveltime tomography based on the adjoint-state method[J].Geophysics,2009,74(6):WCB1-WCB10.
[13]
Vidale J E.Finite-difference calculation of travel times[J].Bulletin of the Seismological Society of America,1988,78(6):2062-2076.
[14]
Van Trier J,Symes W W.Upwind finite-difference calculation of traveltimes[J].Geophysics,1991,56(6):812-821.
[15]
Qin F,Luo Y,Olsen K B,et al.Finite-difference solution of the eikonal equation along expanding wave-fronts[J].Geophysics,1992,57(3):478-487.
[16]
Sethian J A,Popovici A M.3-D traveltime computation using the fast marching method[J].Geophysics,1999,64(2):516-523.
[17]
谢春,刘玉柱,董良国,等.伴随状态法初至波旅行时层析[J].石油地球物理勘探,2014,49(5):877-883.XIE Chun,LIU Yuzhu,DONG Liangguo,et al.First arrival traveltime tomography based on the adjoint state method[J].Oil Geophysical Prospecting,2014,49(5):877-883.
[18]
黄兴国,孙建国,孙章庆,等.基于复程函方程和改进的快速推进法的复旅行时计算方法[J].石油地球物理勘探,2016,51(6):1109-1118.HUANG Xingguo,SUN Jianguo,SUN Zhangqing,et al.Complex traveltime calculation method based on complex function equation and improved fast mar-ching method[J].Oil Geophysical Prospecting,2016,51(6):1109-1118.
[19]
孙章庆,孙建国,韩复兴.三维起伏地表条件下地震波旅行时计算的不等距迎风差分法[J].地球物理学报,2012,55(7):2441-2449.SUN Zhangqing,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.
[20]
Zhao H.A fast sweeping method for Eikonal equations[J].Mathematics of Computation,2004,74(250):603-627.
[21]
Zhang Y,Zhao H,Qian J,et al.High order fast swee-ping methods for static Hamilton-Jacobi equations[J].Journal of Scientific Computing,2006,29(1):25-56.
[22]
Qian J,Zhang Y,Zhao H,et al.Fast sweeping methods for eikonal equations on triangular meshes[J].SIAM Journal on Numerical Analysis,2007,45(1):83-107.
[23]
Fomel S,Luo S,Zhao H,et al.Fast sweeping method for the factored eikonal equation[J].Journal of Computational Physics,2009,228(17):6440-6455.
[24]
Qian J,Zhang Y,Zhao H,et al.A fast sweeping method for static convex Hamilton-Jacobi equations[J].Journal of Scientific Computing,2007,31(1):237-271.
[25]
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.
[26]
Waheed U B,Yarman C E,Flagg G,et al.An iterative,fast-sweeping-based eikonal solver for 3D tilted anisotropic media[J].Geophysics,2015,80(3):C49-C58.
[27]
Waheed U B,Alkhalifah T.A fast sweeping algorithm for accurate solution of the tilted transversely isotro-pic eikonal equation using factorization[J].Geophy-sics,2017,82(6):WB1-WB8.
Alkhalifah T.An acoustic wave equation for orthorhombic anisotropy[J].Geophysics,2003,68(4):1169-1172.
[30]
Bakulin A,Grechka V,Tsvankin I,et al.Estimation of fracture parameters from reflection seismic data-Part Ⅰ:HTI model due to a single fracture set[J].Geophysics,2000,65(6):1788-1802.
[31]
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.
[32]
Tsai Y R,Cheng L,Osher S,et al.Fast sweeping algorithms for a class of Hamilton-Jacobi equations[J].SIAM Journal on Numerical Analysis,2003,41(2):673-694.