- Access by Xinjiang University
Local algorithm for computing complex travel time based on the complex eikonal equation
Phys. Rev. E 93, 043307 – Published 8 April, 2016
DOI: https://doi.org/10.1103/PhysRevE.93.043307
Abstract
The traditional algorithm for computing the complex travel time, e.g., dynamic ray tracing method, is based on the paraxial ray approximation, which exploits the second-order Taylor expansion. Consequently, the computed results are strongly dependent on the width of the ray tube and, in regions with dramatic velocity variations, it is difficult for the method to account for the velocity variations. When solving the complex eikonal equation, the paraxial ray approximation can be avoided and no second-order Taylor expansion is required. However, this process is time consuming. In this case, we may replace the global computation of the whole model with local computation by taking both sides of the ray as curved boundaries of the evanescent wave. For a given ray, the imaginary part of the complex travel time should be zero on the central ray. To satisfy this condition, the central ray should be taken as a curved boundary. We propose a nonuniform grid-based finite difference scheme to solve the curved boundary problem. In addition, we apply the limited-memory Broyden-Fletcher-Goldfarb-Shanno technology for obtaining the imaginary slowness used to compute the complex travel time. The numerical experiments show that the proposed method is accurate. We examine the effectiveness of the algorithm for the complex travel time by comparing the results with those from the dynamic ray tracing method and the Gauss-Newton Conjugate Gradient fast marching method.
Physics Subject Headings (PhySH)
Article Text
References (40)
- L. B. Felsen, J. Opt. Soc. Am. 66, 751 (1976).
- L. B. Felsen, Geophys. J. Roy. Astronom. Soc. 79, 77 (1984).
- V. Červený and M. M. Popov, Geophys. 70, 109 (1982).
- N. Hill, Geophys. 55, 1416 (1990).
- N. Hill, Geophys. 66, 1240 (2001).
- S. Gray, Geophys. 70, S71 (2005).
- V. Červený, L. Klimeš, and I. Pšenčík, Geophys. J. Int. 79, 89 (1984).
- V. Červený, Geophys. 58, 44 (1985).
- V. Červený, Seismic Ray Theory (Cambridge University Press, London, 2001).
- R. Magnanini and G. Talenti, Contemp. Math. 283, 203 (1999).
- R. Magnanini and G. Talenti, Appl. Anal. 85, 249 (2006).
- W. Y. D. Wang and G. A. Deschamps, Proc. IEEE. 62, 1541 (1974).
- D. J. Hearn and E. S. Krebes, Geophys. 55 (1990).
- D. J. Hearn and E. S. Krebes, Pure Appl. Geophys. 132 (1990).
- L. H. T. Le, E. S. Krebes, and G. E. Quiroga-Goode, Geophys. J. Int. 116, 598 (1994).
- T. F. Zhu and K. Y. Chun, Geophys. J. Int. 119, 269 (1994).
- C. J. Thomson, Studia Geophysica et Geodaetica 41, 345 (1997).
- S. J. Chapman, J. M. H. Lawry, J. R. Ockendon, and R. H. Tew, SIAM Rev. 41, 417 (1999).
- Y. A. Kravtsov, G. W. Forbes, and A. A. Asatryan, Prog. Opt. 39, 1 (1999).
- R. A. Egorchenkov and Y. A. Kravtsov, JOSA A 18, 650 (2001).
- Y. A. Kravtsov, Alpha Science (Harrow, UK, 2005).
- A. Hanyga and M. Seredyńska, Pure Appl. Geophys. 157, 679 (2000).
- D. Gajewski and I. Pšencík, Geophys. 57, 27 (1992).
- V. Vavrycuk, Geophys. J. Int. 175 (2008).
- V. Vavrycuk, Geophys. 5, 73 (2008).
- V. Vavrycuk, Geophys. J. Int. 181, 1665 (2010).
- M. Klimeš and L. Klimeš, Geophys. J. Int. 186, 751 (2011).
- V. Vavryčuk, Geophys. 77, T109 (2012).
- V. Červený, Geophys. J. Int. 73, 389 (1983).
- S. Li, S. Fomel, and A. Vladimirsky, SEG Annual Meeting, San Antonio, Texas (SEG, Tulsa, 2011).
- J. A. Sethian, Soc. Ind. Appl. Math. 41, 199 (1999).
- J. A. Sethian, Geophys. 64, 516 (1999).
- J. M. Perry, Discussion Papers. 269 (1977).
- J. Nocedal, Math. Computat. 35, 773 (1980).
- D. C. Liu and J. Nocedal, Math. Prog. 1, 45 (1989).
- R. Magnanini and G. Talenti, SIAM J. Math. Anal. 34, 805 (2003).
- J. Nocedal and S. Wright, Numerical Optimization (Springer-Verlag, New York, 1999).
- J. G. Sun, Z. Q. Sun, and F. X. Han, Geophys. 76, T53 (2011).
- C. Gao, J. G. Sun, and P. Qi, Chinese J. Geophys. 58, 1333 (2015).
- R. S. Wu, Geophys. J. Int. 83, 93 (1985).