1887
Volume 53, Issue 4
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

For the high-precision seismic imaging of complex structures on irregular topographies, this study proposes a Chebyshev spectral-element reverse time migration (CSE-RTM) method coupled with the implicit Newmark time integral method. The approach is first applied to a numerical example to consider factors affecting numerical accuracy, such as the time step and the pattern of subdivision of the computational domain. The results show that the smaller the time step is, the higher the numerical accuracy is and that numerical accuracy improves when the grid nodes tend to be distributed uniformly in the computational domain. The spectral element method combines the advantages of the stronger boundary-related adaptability of the finite element method with the high accuracy and fast convergence of the spectral method. This method is applied to reverse time migration imaging, and tests of the model show that its accuracy of imaging is higher than that of the conventional finite-difference reverse time migration method in case of complex structures in areas with irregular topography. To improve the efficiency of parallel computing of the CSE-RTM, an iterative algorithm is proposed to solve the hierarchical equilibrium sparse linear equations (HELPs) for it based on the idea of the degrees of cohesion of freedom and local relaxation. This algorithm has the advantage that the efficiency of parallel computation does not decrease with an increasing number of processors while the number of iterations required for convergence is controlled. Owing to the fast convergence of the spectral method on a sparse grid coupled with the parallel computational scheme of HELPs, the efficiency of calculation of the CSE-RTM is significantly improved.

© 2021 Australian Society of Exploration Geophysicists
Loading

Article metrics loading...

/content/journals/10.1080/08123985.2021.1986387
2022-07-04
2026-01-23

  • From This Site

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

Full text loading...

References

  1. Fang, Q.1992. Studies on the accuracy of finite element analysis of implicit Newmark method for wave propagation problems. Explosion and Shock Waves12 no. 1: 45–53.
     [Google Scholar]
  2. Fischer, P.F., and A.T.Patera. 1989. Spectral element methods for the incompressible Navier-stokes equations. York: ASME. 71–143. doi:10.1007/978‑1‑4613‑0535‑4_3
     https://doi.org/10.1007/978-1-4613-0535-4_3
  3. Komatitsch, D., and J.Tromp. 2010. Introduction to the spectral element method for three-dimensional seismic wave propagation:. Geophysical Journal of the Royal Astronomical139: 806–22. doi:10.1046/j.1365‑246x.1999.00967.x.
     https://doi.org/10.1046/j.1365-246x.1999.00967.x [Google Scholar]
  4. Komatitsch, D., J.Vilotte, R.Vai, J. M.Castillo‐Covarrubias, and F.J.Sánchez‐Sesma. 1999. The spectral element method for elastic wave equations: application to 2D and 3D seismic problems. International Journal for Numerical Methods in Engineering 45: 1136–64. doi:10.1002/(SICI)1097‑0207(19990730)45:9<1139::AIDNME617>3.0.CO;2.
     https://doi.org/10.1002/(SICI)1097-0207(19990730)45:9<1139::AIDNME617>3.0.CO;2 [Google Scholar]
  5. Lin, W.J., C.Su, and G.Seriani. 2018. The poly-grid spectral element method and its implementation in high performance computing. Journal of Applied Acoustics37 no. 1: 42–51.
     [Google Scholar]
  6. Patera, A.T.1984. A spectral element method for fluid dynamics: laminar flow in a channel expansion. Journal of Computational Physics54 no. 3: 468–88. doi:10.1016/0021‑999(84)90128‑1.
     https://doi.org/10.1016/0021-999(84)90128-1 [Google Scholar]
  7. Priolo, E., and G.Seriani. 1991. A numerical investigation of Chebyshev spectral element method for acoustic wave propagation. In Proceedings of the 13th IMACS Conference on Computation and Applied Mathematics, 551–6 Dublin.
  8. Shan, R., and W.Zhu. 2013. Parallel SEM seismic wave simulation based on GPU. Coal Geology of China25 no. 3: 52–4.
     [Google Scholar]
  9. Shao, X.M., and Z.H.L.Lan. 1995. Absorbing boundary conditions for anisotropic elastic wave equations. Chinese Journal of Geophysics38 no. supp.: 57–73.
     [Google Scholar]
  10. Wang, M.C.2003. Finite element method. Beijing: Tsinghua University Press, 468–522.
  11. Wang, X.M., G.Seriani, and W.J.Lin. 2007. Using spectral element method to calculate some theory problems about elastic wave field. Science in China Press 37 no. 1: 41–59.
     [Google Scholar]
  12. Wang, T.K., R.H.Li, X.F.Li, M.G.Zhang, and G.H.Long. 2007. Numerical spectralelement modeling for seismic wave propagation in transversely isotropic medium. Progress in Geophysics 22 no. 3: 778–84.
     [Google Scholar]
  13. Wang, X.M., G.Seriani, and W.J.Lin. 2007. Using spectral element method to calculate some theory problems about elastic wave field. Science in China Press37 no. 1: 41–59.
     [Google Scholar]
/content/journals/10.1080/08123985.2021.1986387
Loading
Simulating and imaging wavefield of irregular topography based on spectral element method
Exploration Geophysics 53, 398 (2022); https://doi.org/10.1080/08123985.2021.1986387
/content/journals/10.1080/08123985.2021.1986387
/content/journals/10.1080/08123985.2021.1986387
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): algorithm; depth migration; imaging; reverse time; Wave equation

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