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

Abstract

Abstract

Incorporating anisotropy and complex topography is necessary to perform traveltime tomography in complex land environments while being a computational challenge when traveltimes are computed with finite‐difference eikonal solvers. Previous studies have taken this challenge by computing traveltimes in transverse isotropic media involving complex topography with a finite‐difference eikonal equation solver on a curvilinear grid. In this approach, the source singularity, which is a major issue in eikonal solvers, is managed with the elliptical multiplicative factorization method, where the total traveltime field is decomposed into an elliptical base traveltime map, which has a known analytical expression and an unknown perturbation field. However, the group velocity curve can deviate significantly from an ellipse in anellipitically anisotropic media. In this case, the elliptical base traveltime field differs significantly from the anelliptical counterpart, leading to potentially suboptimal traveltime solutions, even though it helps to mitigate the detrimental effects of the source singularity. To overcome this issue, we develop a more accurate topography‐dependent eikonal solver in transverse isotropic media that relies on anelliptical factorization. To achieve this, we first define the coordinate transform from the Cartesian to the curvilinear coordinate system, which provides the necessary framework to implement the topography‐dependent transverse isotropic finite‐difference eikonal solver with arbitrary source and receiver positioning. Then, we develop a semi‐analytical method for the computation of the topography‐dependent anelliptical base traveltime field. Finally, we efficiently solve the resulting quadratic elliptical equation using the fast sweeping method and a quartic anelliptical source term through fixed‐point iteration. We assess the computational efficiency, stability and accuracy of the new eikonal solver against the solver based on elliptical factorization using several transverse isotropic numerical examples. We conclude that this new solver provides a versatile and accurate forward engine for traveltime tomography in complex geological environments such as foothills and thrust belts. It can also be used in marine environments involving complex bathymetry when tomography is applied to redatumed data on the sea bottom.

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

Article metrics loading...

/content/journals/10.1111/1365-2478.13656
2025-01-26
2026-02-14

  • From This Site

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

Full text loading...

References

  1. Alkhalifah, T. (1995) Gaussian beam depth migration for anisotropic media. Geophysics, 60, 1474–1484.
     [Google Scholar]
  2. Alkhalifah, T. (2000) An acoustic wave equation for anisotropic media. Geophysics, 65, 1239–1250.
     [Google Scholar]
  3. Carcione, J., Kosloff, D., Behle, A. & Seriani, G. (1992) A spectral scheme for wave propagation simulation in 3D elastic‐anisotropic media. Geophysics, 57, 1593–1607.
     [Google Scholar]
  4. David Kerlick, G. & Klopfer, G.H. (1982) Assessing the quality of curvilinear coordinate meshes by decomposing the Jacobian matrix. Applied Mathematics and Computation, 10–11, 787–807.
     [Google Scholar]
  5. Detrixhe, M., Gibou, F. & Min, C. (2013) A parallel fast sweeping method for the eikonal equation. Journal of Computational Physics, 237, 46–55.
     [Google Scholar]
  6. Fomel, S., Luo, S. & Zhao, H.‐K. (2009) Fast sweeping method for the factored eikonal equation. Journal of Computational Physics, 228, 6440–6455.
     [Google Scholar]
  7. Gao, K. & Huang, L. (2008) A hybrid eikonal solver for accurate first‐arrival traveltime computation in anisotropic media with strong contrasts. arXiv preprint arXiv:2008.07684v1.
  8. Gray, S.H., Etgen, J., Dellinger, J. & Whitmore, D. (2001) Seismic migration problems and solutions. Geophysics, 66, 1622–1640.
     [Google Scholar]
  9. Gray, S.H. & Marfurt, K.J. (1995) Migration from topography: improving the near‐surface image. Canadian Journal of Exploration Geophysics, 31, 18–24.
     [Google Scholar]
  10. Guo, G., Lan, H., Zhou, X., Liu, Y., Waheed, U.B. & Chen, J. (2022) Topography‐dependent eikonal tomography based on the fast‐sweeping scheme and the adjoint‐state technique. Geophysics, 87, U29–U41.
     [Google Scholar]
  11. Han, S., Zhang, W. & Zhang, J. (2017) Calculating qP‐wave travel times in 2D TTI media by high‐order fast sweeping methods with a numerical quartic equation solver. Geophysical Journal International, 210, 1560–1569.
     [Google Scholar]
  12. Hao, Q. & Stovas, A. (2016) Analytic calculation of phase and group velocities of P‐waves in orthorhombic media. Geophysics, 81, C79–C97.
     [Google Scholar]
  13. Hicks, G.J. (2002) Arbitrary source and receiver positioning in finite‐difference schemes using Kaiser windowed sinc functions. Geophysics, 67, 156–166.
     [Google Scholar]
  14. Huang, G., Luo, S., Deng, J. & Vavryčuk, V. (2020) Traveltime calculations for qP, qSV, and qSH waves in two‐dimensional tilted transversely isotropic media. Journal of Geophysical Research: Solid Earth, 125, e2019JB018868.
     [Google Scholar]
  15. Jin, C. & Zhang, J. (2018) Stereotomography of seismic data acquired on undulant topography. Geophysics, 83, U35–U41.
     [Google Scholar]
  16. Lan, H., Chen, J. & Zhang, Z. (2014) A fast sweeping scheme for calculating P wave first‐arrival travel times in transversely isotropic media with an irregular surface. Pure and Applied Geophysics, 171, 2199–2208.
     [Google Scholar]
  17. Lan, H. & Zhang, Z. (2013) Topography‐dependent eikonal equation and its solver for calculating first‐arrival traveltimes with an irregular surface. Geophysical Journal International, 193, 1010–1026.
     [Google Scholar]
  18. Le Bouteiller, P., Benjemaa, M., Métivier, L. & Virieux, J. (2018) An accurate discontinuous Galerkin method for solving point–source Eikonal equation in 2‐D heterogeneous anisotropic media. Geophysical Journal International, 212, 1498–1522.
     [Google Scholar]
  19. Le Bouteiller, P., Benjemaa, M., Métivier, L. & Virieux, J. (2019) A discontinuous Galerkin fast‐sweeping Eikonal solver for fast and accurate traveltime computation in 3D tilted anisotropic media. Geophysics, 84, C107–C118.
     [Google Scholar]
  20. Lecomte, I. & Kaschwich, T. (2008) Closer to real earth in reservoir characterization: a 3D isotropic/anisotropic PSDM simulator. In 78th Annual SEG meeting and exposition, expanded abstracts, Houston, TX: Society of Exploration Geophysics, pp. 1570–1574.
  21. Lu, Y., Zhang, J., Yang, K., Yang, J. & Li, Z. (2024) A fast solution for the eikonal equation based on quadratic function in weakly tilted transversely isotropic media. IEEE Transactions on Geoscience and Remote Sensing, 62, 1–10.
     [Google Scholar]
  22. Lu, Y. & Zhang, W. (2021) A fast sweeping method for calculating qP‐wave traveltimes in 3‐D vertical transversely isotropic media using a quadratic equation. Geophysical Journal International, 227, 2121–2136.
     [Google Scholar]
  23. Luo, S. & Qian, J. (2012) Fast sweeping method for factored anisotropic eikonal equations: multiplicative and additive factors. Journal of Scientific Computing, 52, 360–382.
     [Google Scholar]
  24. Noble, M., Gesret, A. & Belayouni, N. (2014) Accurate 3‐D finite difference computation of travel time in strongly heterogeneous media. Geophysical Journal International, 199, 1572–1585.
     [Google Scholar]
  25. Pica, A. (1997) Fast and accurate finite difference solution of the 3D eikonal equation parameterized in celerity. In 67th Annual international meeting, expanded abstracts. Houston, TX: Society of Exploration Geophysics, pp. 1774–1777.
     [Google Scholar]
  26. Popovici, A.M. & Sethian, J. (2002) 3D imaging using higher order fast marching traveltimes. Geophysics, 67, 604–609.
     [Google Scholar]
  27. Qian, J. & Symes, W. (2002) An adaptive finite‐difference method for traveltimes and amplitudes. Geophysics, 67, 167–176.
     [Google Scholar]
  28. Qian, J. & Symes, W.W. (2001) Paraxial eikonal solvers for anisotropic quasi‐P travel times. Journal of Computational Physics, 173, 256–278.
     [Google Scholar]
  29. Rawlinson, N. & Sambridge, M. (2004) Multiple reflection and transmission phases in complex layered media using a multistage fast marching method. Geophysics, 69(5), 1328–1350.
     [Google Scholar]
  30. Sadarjoen, I.A., De Leeuw, W.C. & Post, F.H. (2001) Visualization techniques for curvilinear grids. Technical report. Report 95‐138. Delft, the Netherlands: TU Delft ‐ Delft University of Technology.
  31. Sambolian, S., Górszczyk, A., Operto, S., Ribodetti, A. & Tavakoli, B. (2021) Mitigating the ill‐posedness of first‐arrival traveltime tomography with slopes: application to the eastern Nankai Trough OBS dataset (Japan). Geophysical Journal International, 227, 898–921.
     [Google Scholar]
  32. Sambolian, S., Operto, S., Ribodetti, A., Tavakoli, B. & Virieux, J. (2019) Parsimonious slope tomography based on eikonal solvers and the adjoint‐state method. Geophysical Journal International, 218, 456–478.
     [Google Scholar]
  33. Sambolian, S., Operto, S., Ribodetti, A. & Virieux, J. (2021) Consistent seismic event location and subsurface parameters inversion through slope tomography: a variable‐projection approach. Geophysical Journal International, 224, 1956–1979.
     [Google Scholar]
  34. Shragge, J. (2016) Acoustic wave propagation in tilted transversely isotropic media: incorporating topography. Geophysics, 81, C265–C278.
     [Google Scholar]
  35. Tavakoli F, B., Operto, S., Ribodetti, A. & Virieux, J. (2017) Slope tomography based on eikonal solvers and the adjoint‐state method. Geophysical Journal International, 209(3), 1629–1647.
     [Google Scholar]
  36. Tavakoli F, B., Operto, S., Ribodetti, A. & Virieux, J. (2019) Matrix‐free anisotropic slope tomography: theory and application. Geophysics, 84(1), R35–R57.
     [Google Scholar]
  37. Tavakoli F, B., Ribodetti, A., Virieux, J. & Operto, S. (2015) An iterative factored eikonal solver for TTI media. In SEG technical program expanded abstracts 2015, 85th Annual SEG Meeting, New Orleans, LA. Houston, TX: Society of Exploration Geophysicists, pp. 3576–3581.
     [Google Scholar]
  38. Thomsen, L.A. (1986) Weak elastic anisotropy. Geophysics, 51, 1954–1966.
     [Google Scholar]
  39. Tsvankin, I. (2005) Seismic signature and analysis of reflection data in anisotropic media. In Seismic Exploration, 2nd edition. Handbook of Geophysical Exploration, volume 29. Amsterdam, the Netherlands: Elsevier, pp. 21–22.
     [Google Scholar]
  40. Tsvankin, I., Gaiser, J., Grechka, V., van der Baan, M. & Thomsen, L. (2010) Seismic anisotropy in exploration and reservoir characterization: An overview. Geophysics, 75, 75A15–75A29.
     [Google Scholar]
  41. Vidale, D. (1988) Finite‐difference calculation of travel time. Bulletin of the Seismological Society of America, 78, 2062–2076.
     [Google Scholar]
  42. Waheed, U. & Alkhalifah, T. (2017) A fast sweeping algorithm for accurate solution of the tilted transversely isotropic eikonal equation using factorization. Geophysics, 82.
     [Google Scholar]
  43. Waheed, U.B., Yarman, C.E. & Flagg, G. (2015) An iterative, fast‐sweeping‐based eikonal solver for 3D tilted anisotropic media. Geophysics, 80(3), C49–C58.
     [Google Scholar]
  44. Zhang, Q., Ma, X. & Nie, Y. (2021) An iterative fast sweeping method for the eikonal equation in 2D anisotropic media on unstructured triangular meshes. Geophysics, 86, U49–U61.
     [Google Scholar]
  45. Zhao, H. (2005) A fast sweeping method for eikonal equations. Mathematics of Computation, 74, 603–627.
     [Google Scholar]
  46. Zhou, X., Lan, H., Chen, L., Guo, G., Lei, Y., Waheed, U.B. & Pan, S. (2021) An iterative factored topography‐dependent eikonal solver for anisotropic media. Geophysics, 86, U121–U134.
     [Google Scholar]
  47. Zhou, X., Lan, H., Chen, L., Guo, G., Waheed, U.B. & Badal, J. (2023) A topography‐dependent eikonal solver for accurate and efficient computation of traveltimes and their derivatives in 3D heterogeneous media. Geophysics, 88, U17–U29.
     [Google Scholar]
/content/journals/10.1111/1365-2478.13656
Loading
Topography‐dependent eikonal solver in tilted transverse isotropic media with anelliptical factorization
Geophysical Prospecting 73, 595 (2025); https://doi.org/10.1111/1365-2478.13656
/content/journals/10.1111/1365-2478.13656
/content/journals/10.1111/1365-2478.13656
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): anisotropy; computing aspects; finite difference methods; model building; modelling; seismic anisotropy; traveltime; two‐dimensional finite difference modelling

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