A method for complex traveltime calculation based on the complex eikonal equation and the modified fast marching
Huang Xingguo1,2, Sun Jianguo1,2, Sun Zhangqing1,2, Wang Qianlong1,2
1. College of Geo-exploration Science and Technology, Jilin University, Changchun, Jilin 130026, China;
2. Laboratory for Integrated Geophysical Interpretation Theory, the Ministry of Land and Resources, Changchun, Jilin 130026, China
Abstract:It is difficult to calculate the complex traveltime by directly solving the complex eikonal equation,because this equation constitutes of the real part and the imaginary part.We propose an optimization method by combining the perturbation theory with the intersection of product.This optimization solves first complex eikonal equation,then calculates complex traveltime.To calculate the complex traveltime of bending ray,we introduce non-uniform grid upwind difference method and a modified fast marching method to deal with the area near the central ray.Numerical test results show that the proposed method and strategies in this paper are effective.
黄兴国, 孙建国, 孙章庆, 王乾龙. 基于复程函方程和改进的快速推进法的复旅行时计算方法[J]. 石油地球物理勘探, 2016, 51(6): 1109-1118,1147.
Huang Xingguo, Sun Jianguo, Sun Zhangqing, Wang Qianlong. A method for complex traveltime calculation based on the complex eikonal equation and the modified fast marching. OGP, 2016, 51(6): 1109-1118,1147.
Felsen L B. Evanescent waves. Journal of the Optical Society of America,1976,66(8):751-760.
[2]
Felsen L B. Geometrical theory of diffraction,evanescent waves,complex rays,and Gaussian beams. Geophysical Journal of the Royal Astronomical Society,1984,79(1):77-88.
[3]
?ervený V,Popov M M. Computation of wave fields in inhomogeneous media-Gaussian beam approach. Geophysical Journal of the Royal Astronomical Society,1982,70(1):109-128.
[4]
?ervený V,Klimeš L,PšeníkI. Paraxial ray approximations in the computation of seismic wavefields in inhomogeneous media. Geophysical Journal International,1984,79(1):89-104.
[5]
?ervený V. Gaussian beam synthetic seismograms.Geophysics,1985,58(1-3):44-72.
[6]
?ervený V. Seismic Ray Theory. Combridge University Press,London,2001.
[7]
Magnanini R,Talenti G. On complex-valued solutions to a 2D Eikonal equation, Part One:qualitative properties. Contemporary Mathematics,1999,283(1):203-229.
[8]
Li S,Fomel S,Vladimirsky A. Improving wave-equation fidelity of Gaussian beams by solving the complex eikonal equation. SEG Technical Program Expanded Abstracts,2011,30:3829-3834.
[9]
Klimeš M,Klimeš L. Perturbation expansions of complex-valued traveltime along real-valued reference rays. Geophysical Journal International,2011,186(2):751-759.
[10]
Vavry uk V. On numerically solving the complex eikonal equation using real ray-tracing methods:A comparison with the exact analytical solution. Geophysics,2012,77(4):T109-T116.
[11]
Sun J G,Sun Z Q,Han F X. A finite difference scheme for solving the eikonal equation including surface topography. Geophysics,2011,76(4):T53-T63.
[12]
Magnanini R,Talenti G. On complex-valued solutions to a 2D eikonal equation, Part three:analysis of a Bäcklund transformation. Applicable Analysis,2006,85(1-3):249-276.
[13]
Huang X,Sun J,Sun Z. Local algorithm for computing complex travel time based on the complex eikonal equation. Physical Review E,2016,93(4):043307.
[14]
Sethian J A. Fast marching methods. Society for Industrial and Applied Mathematics,1999,41(2):199-235.
[15]
Sethian J A.3-D traveltime computation using the fast marching method. Geophysics,1999,64(2):516-523.
Sun Zhangqing,Sun Jianguo,Han Fuxing. The method of computation of traveltime based on linear interpolation and narrow band technique. OGP,2009,44(4):436-441.
[17]
Wu R S. Gaussian beams,complex rays,and the analytic extension of the Green's function in smoothly inhomogeneous media. Geophysical Journal International,1985,83(1):93-110.