1887
Volume 72, Issue 5
  • E-ISSN: 1365-2478

Abstract

Abstract

We have developed a numerical scheme for the second‐order acoustic wave equation based on the Lie product formula and Taylor‐series expansion. The scheme has been derived from the analytical solution of the wave equation and in the approximation of the time derivative for a wavefield. Through these two equations, we obtained the first‐order differential equation in time, where the time evolution operator of the analytic solution of this differential equation is written as a product of exponential matrices. The new numerical solution using a Lie product formula may be combined with Taylor‐series, Chebyshev, Hermite and Legendre polynomial expansion or any other expansion for the cosine function. We use the proposed scheme combined with the second‐ or fourth‐order Taylor approximations to propagate the wavefields in a recursive procedure, in a stable manner, accurately and efficiently with even larger time steps than the conventional finite‐difference method. Moreover, our numerical scheme has provided results with the same quality as the rapid expansion method but requiring fewer computations of the Laplacian operator per time step. The numerical results have shown that the proposed scheme is efficient and accurate in seismic modelling, reverse time migration and least‐squares reverse time migration.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.13481
2024-05-21
2026-01-25
Loading full text...

Full text loading...

References

  1. Araujo, E. & Pestana, R. (2022) Viscoelastic and viscoacoustic modelling using the Lie product formula. Geophysical Prospecting, 70, 1295–1312.
    [Google Scholar]
  2. Araujo, E.S. & Pestana, R.C. (2019) Time‐stepping wave‐equation solution for seismic modeling using a multiple‐angle formula and Taylor expansion. Geophysics, 84, T299–T311.
    [Google Scholar]
  3. Araujo, E.S. & Pestana, R.C. (2021) Time evolution of the first‐order linear acoustic/elastic wave equation using Lie product formula and Taylor expansion. Geophysical Prospecting, 69, 70–84. https://doi.org/10.1111/1365‐2478.13033
    [Google Scholar]
  4. Araujo, E.S., Pestana, R.C. & Santos, A. W.G. (2014) Symplectic scheme and the Poynting vector in reverse‐time migration. Geophysics, 79, S163–S172.
    [Google Scholar]
  5. Carcione, J.M., Kosloff, D. & Kosloff, R. (1988) Wave propagation simulation in a linear viscoacoustic medium. Geophysics, 93, 393–407.
    [Google Scholar]
  6. Cerjan, C., Kosloff, D., Kosloff, R. & Reshef, M. (1985) A nonreflecting boundary condition for discrete acoustic and elastic wave equations. Geophysics, 50, 705–708.
    [Google Scholar]
  7. Chen, J. (2009) Lax‐Wendroff and Nyström methods for seismic modelling. Geophysical Prospecting, 57, 931–941.
    [Google Scholar]
  8. Chen, J.M. (2007) High‐order time discretizations in seismic modeling. Geophysics, 72, SM115–SM122.
    [Google Scholar]
  9. Cherifi, C., Cherifi, H., Karsai, M. & Musolesi, M. (2017) Complex networks & their applications VI, Proceedings of Complex Networks 2017 (The Sixth International Conference on Complex Networks and Their Applications). Springer.
    [Google Scholar]
  10. Dablain, M.A. (1986) The application of high‐order differencing to the scalar wave equation. Geophysics, 51, 54–66.
    [Google Scholar]
  11. Defez, E., Sastre, J., Ibáñez, J. & Ruiz, P.A. (2009) Computing matrix functions solving coupled differential models. Mathematical and Computer Modelling, 50, 831–839.
    [Google Scholar]
  12. Etgen, J. (1986) High‐order finite‐difference reverse time migration with the 2‐way non‐reflecting. SEP48, 133–146.
    [Google Scholar]
  13. Fomel, S.L., Ying, L. & Song, X. (2013) Seismic wave extrapolation using low rank symbol approximation. Geophysical Prospecting, 61, 526–536.
    [Google Scholar]
  14. Geodynamics . (2023). seismic_cpml sharewave. https://github.com/geodynamics/seismic_cpml (accessed 4 March 2023).
  15. Hall, B.C. (2015) Lie groups, Lie algebras, and representations: An elementary introduction. Berlin: Springer.
    [Google Scholar]
  16. Kosloff, D., Filho, A., Tessmer, E. & Behle, A. (1989) Numerical solution of the acoustic and elastic wave equation by new rapid expansion method. Geophysical Prospecting, 37, 383–394.
    [Google Scholar]
  17. Pestana, R.C. & Stoffa, P.L. (2010) Time evolution of the wave equation using rapid expansion method. Geophysics, 75, 121–131.
    [Google Scholar]
  18. Plessix, R.E. (2006) A review of the adjoint‐state method for computing the gradient of a functional with geophysical applications. Geophysical Journal International, 167, 495–503.
    [Google Scholar]
  19. Soubaras, R. & Zhang, Y. (2008). Two‐steps explicit marching method for reverse time migration, F041, 70th EAGE Conference and Exhibition, Rome/Italy: Houten, the Netherlands: EAGE.
    [Google Scholar]
  20. Stickel, E.U. (1994) A splitting method for the calculation of the matrix exponential. Analysis, 14, 103–112.
    [Google Scholar]
  21. Tal‐Ezer, H. (1986) Spectral methods in time for hyperbolic problems. SIPM Journal on Numerical Analysis, 23, 11–20.
    [Google Scholar]
  22. Zhang, Y. & Zhang, G. (2009) One‐step extrapolation method for reverse time migration. Geophysics, 74, A29–A33.
    [Google Scholar]
  23. Zhu, F., Huang, J. & Yu, H. (2018) Least‐squares Fourier finite‐difference pre‐stack depth migration for VTI media. Journal of Geophysics and Engineering, 15, 421–437.
    [Google Scholar]
/content/journals/10.1111/1365-2478.13481
Loading
/content/journals/10.1111/1365-2478.13481
Loading

Data & Media loading...

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