A fast algorithm for computing complex traveltime based on L-BFGS method
HUANG Xing-Guo1,2, SUN Jian-Guo1,2, SUN Zhang-Qing1,2, WANG Qian-Long1,2
1. College for Geoexploration Science and Technology, Jilin University, Changchun 130026;
2. Laboratory for Integrated Geophysical Interpretation Theory of the Ministry for Land and Resources of China—Laboratory for Wave Theory and Imaging Technology, Changchun 130026
Complex traveltime of seismic wave plays an important role in dealing with caustics with which geometric ray theory faces.To obtain a higher accuracy and a more efficient method for calculating the complex traveltime,the authors introduce the L-BFGS optimization method to calculate the imaginary slowness.In addition,orthogonal condition of the real part and imaginary part of complex traveltime are taken as the objective function,and then a new L-BFGS-FMM method is obtained for computing complex traveltime.The numerical examples demonstrate the effectiveness of L-BFGS-FMM used in this paper.An analysis of accuracy and efficiency of the results obtained by using the analytical method,dynamic ray tracing method,Gauss-Newton-conjugate gradient method and L-BFGS-FMM method shows that the L-BFGS-FMM method introduced in this paper is superior to all the other methods,and hence the method introduced in this paper can meet the large-scale computing requirement in practical applications.
黄兴国, 孙建国, 孙章庆, 王乾龙. 基于L-BFGS理论求解复程函方程的地震波复走时计算方法[J]. 物探与化探, 2016, 40(5): 961-967.
HUANG Xing-Guo, SUN Jian-Guo, SUN Zhang-Qing, WANG Qian-Long. A fast algorithm for computing complex traveltime based on L-BFGS method. Geophysical and Geochemical Exploration, 2016, 40(5): 961-967.
[1] Felsen L B.Evanescent waves[J].Journal of the Optical Society of America,1976,66(8):751-760.
[2] ?ervený V,Popov M MI.Pšen?ík I.Computation of wave fields in inhomogeneous media[J].Geophysics,1982,70(1):109-128.
[3] Popov M M.A new method of computation of wave fields using Gaussian beams[J].Wave Motion,1982,4(1):85-97.
[4] Hill N R.Gaussian beam migration[J].Geophysics,1990,55(11):1416-1428.
[5] Hill N R.Prestack Gaussian beam depth migration[J].Geophysics,2001,66(4):1240-1250.
[6] Magnanini R,Talenti G.On complex-valued solutions to a 2D Eikonal equation. Part One:qualitative properties[J]. Contemporary Mathematics,1999,283:203-229.
[7] Li S,Fomel S,Vladimirsky A.Improving wave-equation fidelity of Gaussian beams by solving the complex eikonal equation[C]//San Antonio:Expanded Abstracts of 81stSEG Annual Internet Meeting,2011:3829-3834.
[8] Sethian J A,Popovici A M.3-D traveltime computation using the fast marching method[J].Geophysics,1999,64(2):516-523.
[9] Huang X,Sun J,Sun Z.Local algorithm for computing complex travel time based on the complex eikonal equation[J]. Physical Review E,2016,93(4):043307(1)-043307(10).
[10] Perry A.A class of conjugate gradient algorithms with a two-step variable-metric memory[J].Discussion Papers from Northwestern University,Center for Mathematical Studies in Economics and Management Science,1977,269.
[11] Nocedal J.Updating quasi-Newton matrices with limited storage[J].Mathematics of computation,1980,35(151):773-782.
[12] Liu D C,Nocedal J.On the limited memory BFGS method for large scale optimization[J]. Mathematical programming,1989,45:503-528.
[13] Bleistein N.Mathematical methods for wave phenomena[M].Orlando:Academic Press,1984.
[14] Magnanini R,Talenti G.On complex-valued solutions to a two-dimensional eikonal equation. Ⅱ. Existence theorems[J].SIAM journal on mathematical analysis,2003,34(4):805-835.
[15] Schubert L K.Modification of a quasi-Newton method for nonlinear equations with a sparse Jacobian[J].Mathematics of Computation,1970,24(109):27-30.
[16] Deuflhard P.A modified Newton method for the solution of ill-conditioned systems of nonlinear equations with application to multiple shooting[J].Numerische Mathematik,1974,22(4):289-315.
[17] Dennis,Jr J E,Moré J J.Quasi-Newton methods,motivation and theory[J].SIAM review,1977,19(1):46-89.
[18] Nocedal J,Wright S.Numerical optimization[M].New York:Springer Science & Business Media,2006.
[19] Wu R S.Gaussian beams,complex rays,and the analytic extension of the Green's function in smoothly inhomogeneous media[J]. Geophysical Journal International,1985,83(1):93-110.