1887
Volume 73, Issue 8
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Raypath tracing is a commonly used technique in geophysics, employed to simulate and analyse seismic wave propagation paths from source to receiver in complex media. In isotropic media, raypaths can be obtained by tracing from the receiver point along directions perpendicular to the wavefront towards the source point, based on the Fermat principle, because in isotropic media, the ray direction aligns with the ray gradient direction. In an anisotropic medium, the ray direction generally differs from the ray gradient direction, rendering the conventional tracing method inaccurate. Solving raypaths using Hamilton's canonical equations is a powerful method. However, in anisotropic media, the complex dependence of wave velocity on the propagation direction complicates the Hamiltonian function, significantly increasing computational complexity. To address this problem, we have derived a scheme based on the relationship between the group velocity vector and the slowness vector in anisotropic media. Firstly, the slowness vector is derived from the traveltime obtained through the eikonal equation, followed by the computation of the group velocity vector. Then, the raypath is determined by tracing back from the receiver point using the group velocity components to the source point. The efficiency and accuracy of our approach are validated through three numerical experiments.

© 2025 European Association of Geoscientists & Engineers. © 2025 European Association of Geoscientists & Engineers.
Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.70087
2025-10-17
2025-11-09

  • From This Site

    /content/journals/10.1111/1365-2478.70087
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. A Sethian, J.1996. “A Fast Marching Level Set Method for Monotonically Advancing Fronts.” Proceedings of the National Academy of Sciences of the United States of America93: 1591–1595.
     [Google Scholar]
  2. B Waheed, U.2020. “A Fast‐Marching Eikonal Solver for Tilted Transversely Isotropic Media.” Geophysics85: S385–S393.
     [Google Scholar]
  3. Bóna, A., M.Slawinski, and P.Smith. 2009. “Ray Tracing by Simulated Annealing: Bending Method.” Geophysics74: T25–T32.
     [Google Scholar]
  4. Bouteiller, P. L., M.Benjemaa, L.Métivier, and J.Virieux. 2018. “An Accurate Discontinuous Galerkin Method for Solving Point‐Source Eikonal Equation in 2D Heterogeneous Anisotropic Media.” Geophysical Journal International212: 1498–1522.
     [Google Scholar]
  5. Bregman, N. D., R. C.Bailey, and C. H.Chapman. 1989. “Crosshole Seismic Tomography.” Geophysics54: 200–215.
     [Google Scholar]
  6. Červený, V.2001. Seismic Ray Method. Cambridge University Press.
     [Google Scholar]
  7. E Vidale, J.1988. “Finite‐Difference Calculation of Travel Times.” Bulletin of the Seismological Society of America78: 2062–2076.
     [Google Scholar]
  8. E Vidale, J.1990. “Finite‐Difference Calculation of Traveltimes in Three Dimensions.” Geophysics55: 521–526.
     [Google Scholar]
  9. Fomel, S., S.Luo, and H.Zhao. 2009. “Fast Sweeping Method for the Factored Eikonal Equation.” Journal of Computational Physics228: 6440–6455.
     [Google Scholar]
  10. Golikov, P., and A.Stovas. 2012. “Traveltime Parameters in Tilted Transversely Isotropic Media.” Geophysics77: C43–C55.
     [Google Scholar]
  11. Grubas, S., A.Duchkov, and G.Loginov. 2022. “Neural Eikonal Solver: Improving Accuracy of Physics‐Informed Neural Networks for Solving Eikonal Equation in Case of Caustics.” Journal of Computational Physics474: 111789.
     [Google Scholar]
  12. Han, S., W.Zhang, and J.Zhang. 2017. “Calculating Qp‐Wave Traveltimes in 2‐D TTI Media by High‐Order Fast Sweeping Methods With a Numerical Quartic Equation Solver.” Geophysical Journal International210: 1560–1569.
     [Google Scholar]
  13. Huang, G., and S.Luo. 2020. “Hybrid Fast Sweeping Methods for Anisotropic Eikonal Equation in Two‐Dimensional Tilted Transversely Isotropic Media.” Journal of Scientific Computing00: 84.
     [Google Scholar]
  14. Huang, G., S.Luo, J.Deng, and V.Vavryčuk. 2020. “Traveltime Calculations for qP, qSV, and qSH Waves in Two‐Dimensional Tilted Transversely Isotropic Media.” Journal of Geophysical Research: Solid Earth125: e2019JB018868.
     [Google Scholar]
  15. Koketsu, K., and S.Sekine. 1998. “Pseudo‐Bending Method for Three‐Dimensional Seismic Ray Tracing in a Spherical Earth With Discontinuities.” Geophysical Journal International132: 339–346.
     [Google Scholar]
  16. Lan, H., and L.Chen. 2018. “An Upwind Fast Sweeping Scheme for Calculating Seismic Wave First‐Arrival Travel Times for Models With an Irregular Free Surface.” Geophysical Prospecting66: 327–341.
     [Google Scholar]
  17. Lu, Y., and W.Zhang. 2021. “A Fast Sweeping Method for Calculating qP‐wave Traveltimes in 3D Vertical Transversely Isotropic Media Using a Quadratic Equation.” Geophysical Journal International227: 2121–2136.
     [Google Scholar]
  18. Ludlam, J., K.Tant, V.Dolean, and A.Curtis. 2023. “Travel Times and Ray Paths for Acoustic and Elastic Waves in Generally Anisotropic Media.” Journal of Computational Physics494: 112500.
     [Google Scholar]
  19. Luo, S., and J.Qian. 2012. “Fast Sweeping Methods for Factored Anisotropic Eikonal Equations: Multiplicative and Additive Factors.” Journal of Scientific Computing52: 360–382.
     [Google Scholar]
  20. Ma, T., and Z.Zhang. 2014. “Calculating Ray Paths for First‐Arrival Travel Times Using a Topography‐Dependent Eikonal Equation Solver.” Bulletin of the Seismological Society of America104: 1501–1517.
     [Google Scholar]
  21. Qian, J., Y.Zhang, and H.Zhao. 2007a. “A Fast Sweeping Method for Static Convex Hamilton‐Jacobi Equations.” Journal of Scientific Computing31: 237–271.
     [Google Scholar]
  22. Qian, J., Y.Zhang, and H.Zhao. 2007b. “Fast Sweeping Methods for Eikonal Equations on Triangular Meshes.” SIAM Journal on Numerical Analysis45: 83–107.
     [Google Scholar]
  23. R Cassell, B.1982. “A Method for Calculating Synthetic Seismograms in Laterally Varying Media.” Geophysical Journal International69: 339–354.
     [Google Scholar]
  24. Rawlinson, N., J.Hauser, and M.Sambridge. 2008. “Seismic Ray Tracing and Wavefront Tracking in Laterally Heterogeneous Media.” Advances in Geophysics49: 203–273.
     [Google Scholar]
  25. Sadri, M., and M. A.Riahi. 2010. “Ray Tracing and Amplitude Calculation in Anisotropic Layered Media.” Geophysical Journal International180: 1170–1180.
     [Google Scholar]
  26. Stovas, A., and T.Alkhalifah. 2013. “Mapping of moveout in tilted transversely isotropic media.” Geophysical Prospecting61, no. 6: 1171–1177.
     [Google Scholar]
  27. Stovas, A., Y.Roganov, and V.Roganov. 2021. “Geometrical Characteristics of Phase and Group Velocity Surfaces in Anisotropic Media.” Geophysical Prospecting69: 53–69.
     [Google Scholar]
  28. Thomsen, L.1986. “Weak Elastic Anisotropy.” Geophysics51: 1954–1957.
     [Google Scholar]
  29. Um, J., and C.Thurber. 1987. “A Fast Algorithm for Two‐Point Seismic Ray Tracing.” The Bulletin of the Seismological Society of America77: 972–986.
     [Google Scholar]
  30. Vavryčuk, V.2006. “Calculation of the Slowness Vector From the Ray Vector in Anisotropic Media.” Proceedings of the Royal Society of London Series A462: 883–896.
     [Google Scholar]
  31. Waheed, U. B., and T.Alkhalifah. 2017. “A Fast Sweeping Algorithm for Accurate Solution of the Tilted Transversely Isotropic Eikonal Equation Using Factorization.” Geophysics82: WB1–WB8.
     [Google Scholar]
  32. Waheed, U. B., C. E.Yarman, and G.Flagg. 2015. “An Iterative, Fast‐Sweeping‐Based Eikonal Solver for 3D Tilted Anisotropic Media.” Geophysics80: C49–C58.
     [Google Scholar]
  33. Wang, Y.2013. “Simultaneous Computation of Seismic Slowness Paths and the Traveltime Field in Anisotropic Media.” Geophysical Journal International195: 1141–1148.
     [Google Scholar]
  34. Wang, Y.2014. “Seismic Ray Tracing in Anisotropic Media: A Modified Newton Algorithm for Solving Highly Nonlinear Systems.” Geophysics79: T1–T7.
     [Google Scholar]
  35. Yaskevich, S., A. A.Duchkov, and Y.Ivanov. 2020. “Approximate Traveltime Inversion in Downhole Microseismic Monitoring.” Geophysical Prospecting68: 918–925.
     [Google Scholar]
  36. Zhang, K., H.Yang, Z.Li, and M.Hu. 2024. “Double‐Difference Constrained Reflection Tomography in Two‐Dimensional Elastic Media.” Geophysical Prospecting72, no. 6: 2245–2260.
     [Google Scholar]
  37. Zhang, L., and B.Zhou. 2018. “Calculation of Slowness Vectors from Ray Directions for qP‐, qSV‐, and qSH‐waves in Tilted Transversely Isotropic Media.” Geophysics83: C153–C160.
     [Google Scholar]
  38. Zhang, Q., X.Ma, and Y.Nie. 2021. “An Iterative Fast Sweeping Method for the Eikonal Equation in 2D Anisotropic Media on Unstructured Triangular Meshes.” Geophysics86: U49–U61.
     [Google Scholar]
  39. Zhao, H.2005. “A Fast Sweeping Method for Eikonal Equations.” Mathematics of Computation74: 603–628.
     [Google Scholar]
  40. Zhao, H.2007. “Parallel Implementations of the Fast Sweeping Method.” Journal of Computational Mathematics25: 421–429.
     [Google Scholar]
  41. Zhou, X., H.Lan, L.Chen, et al. 2021. “An Iterative Factored Topography‐Dependent Eikonal Solver for Anisotropic Media.” Geophysics86: U121–U134.
     [Google Scholar]
/content/journals/10.1111/1365-2478.70087
Loading
An Efficient Method for Calculating Raypaths of First‐Arrival Traveltimes in Transversely Isotropic Media
Geophysical Prospecting 73, e70087 (2025); https://doi.org/10.1111/1365-2478.70087
/content/journals/10.1111/1365-2478.70087
/content/journals/10.1111/1365-2478.70087
Loading

Data & Media loading...

  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error