1887
Volume 51, Issue 3
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

ABSTRACT

To improve the efficiency of seismic waveform inversion, we construct a novel high-order (6th-order) frequency-domain nearly analytic discrete method for forward modelling in inversion processes. Compared with some existing numerical schemes, this method is more powerful in suppressing numerical dispersion and enhancing the accuracy of computational results based on coarse discrete grids. We describe the discretisation of the frequency-domain wave equation, investigate the efficiency of wave-field simulation and perform numerical dispersion analysis. The inversion efficiency can also be improved when the inexact rotated block triangular preconditioned Krylov iteration solvers are incorporated into the corresponding linear system. Moreover, the inversion behaviour of frequency-domain NAD methods is also systematically discussed as one of its important applications. We take three representative media models as examples and the reliable inversion results are obtained reflecting the effectiveness of the proposed methods.

Loading

Article metrics loading...

/content/journals/10.1080/08123985.2019.1699786
2020-05-03
2026-01-18
Loading full text...

Full text loading...

References

  1. Cerjan, C., D. Kosloff, R. Kosloff, and M. Reshef 1985 A nonreflecting boundary condition for discrete acoustic and elastic wave equations. Geophysics50, no. 4: 705–8.
    [Google Scholar]
  2. Engquist, B., and A. Majda 1977 Absorbing boundary conditions for numerical simulation of waves. Proceedings of the National Academy of Sciences74, no. 5: 1765–6.
    [Google Scholar]
  3. Erlangga, Y.A., C. Vuik, and C.W. Oosterlee 2004 On a class of preconditioners for solving the Helmholtz equation. Applied Numerical Mathematics50, no. 3–4: 409–25.
    [Google Scholar]
  4. Ernst, O.G., and M.J. Gander 2012 Why it is difficult to solve Helmholtz problems with classical iterative methods. In Numerical analysis of multiscale problems, 325–63. Berlin, Heidelberg: Springer.
    [Google Scholar]
  5. Gander, M.J., L. Halpern, and F. Magoules 2007 An optimized Schwarz method with two-sided Robin transmission conditions for the Helmholtz equation. International Journal for Numerical Methods in Fluids55, no. 2: 163–75.
    [Google Scholar]
  6. Gander, M.J., and F. Nataf 2005 An incomplete LU preconditioner for problems in acoustics. Journal of Computational Acoustics13, no. 3: 455–76.
    [Google Scholar]
  7. Hager, W.W., and H. Zhang 2006 A survey of nonlinear conjugate gradient methods. Pacific Journal of Optimization2, no. 1: 35–58.
    [Google Scholar]
  8. Jo, C.H., C. Shin, and J.H. Suh 1996 An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator. Geophysics61, no. 2: 529–37.
    [Google Scholar]
  9. Kolb, P., F. Collino, and P. Lailly 1986 Pre-stack inversion of a 1-D medium. Proceedings of the IEEE74, no. 3: 498–508.
    [Google Scholar]
  10. Komatitsch, D., and J. Tromp 2003 A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation. Geophysical Journal International154, no. 1: 146–53.
    [Google Scholar]
  11. Lailly, P. 1983 The seismic inverse problem as a sequence of before stack migrations. In Conference on inverse scattering: Theory and application, ed. J Bee Bednar, 206–20. Philadelphia: Society for Industrial and Applied Mathematics.
    [Google Scholar]
  12. Laird, A.L., and M.B. Giles 2002 Preconditioned iterative solution of the 2D Helmholtz equation.
  13. Lang, C., and Z.R. Ren 2015 Inexact rotated block triangular preconditioners for a class of block two-by-two matrices. Journal of Engineering Mathematics93, no. 1: 87–98.
    [Google Scholar]
  14. Lang, C., and D.H. Yang 2016 A nearly analytic discrete method for solving the acoustic-wave equations in the frequency domain. Geophysics82, no. 1: T43–57.
    [Google Scholar]
  15. Liu, Q.H. 1997 The PSTD algorithm: A time-domain method requiring only two cells per wavelength. Microwave and Optical Technology Letters15, no. 3: 158–65.
    [Google Scholar]
  16. Liu, S.L., C. Lang, H. Yang, and W.S. Wang 2018 A developed nearly analytic discrete method for forward modeling in the frequency domain. Journal of Applied Geophysics149: 25–34.
    [Google Scholar]
  17. Liu, S.L., X.F. Li, W.S. Wang, and Y.S. Liu 2014 A mixed-grid finite element method with PML absorbing boundary conditions for seismic wave modeling. Journal of Geophysics and Engineering11, no. 5: 055009.
    [Google Scholar]
  18. Liu, S.L., X.F. Li, W.S. Wang, Y.S. Liu, M.G. Zhang, and H. Zhang 2014 A new kind of optimal second-order symplectic scheme for seismic wave simulations. Science China Earth Sciences57, no. 4: 751–8.
    [Google Scholar]
  19. Liu, S.L., X.F. Li, W.S. Wang, and T. Zhu 2015 Source wavefield reconstruction using a linear combination of the boundary wavefield in reverse time migration. Geophysics80, no. 6: S203–12.
    [Google Scholar]
  20. Liu, S.L., I. Suardi, D.H. Yang, S. Wei, and P. Tong 2018 Teleseismic traveltime tomography of Northern Sumatra. Geophysical Research Letters45, no. 24: 13231–9.
    [Google Scholar]
  21. Liu, S.L., D.H. Yang, X.P. Dong, X. Liu, and Y. Zhang 2017 Element-by-element parallel spectral-element methods for 3-D teleseismic wave modeling. Solid Earth8, no. 5: 969–86.
    [Google Scholar]
  22. Liu, S.L., D.H. Yang, and C. Lang 2017 Modified symplectic schemes with nearly-analytic discrete operators for acoustic wave simulations. Computer Physics Communications213: 52–63.
    [Google Scholar]
  23. Liu, S.L., D.H. Yang, and J. Ma 2017 A modified symplectic PRK scheme for seismic wave modeling. Computers & Geosciences99: 28–36.
    [Google Scholar]
  24. Marfurt, K.J. 1984 Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations. Geophysics49, no. 5: 533–49.
    [Google Scholar]
  25. Moczo, P., J.O.A. Robertsson, and L. Eisner 2007 The finite-difference time-domain method for modeling of seismic wave propagation. Advances in Geophysics48: 421–516.
    [Google Scholar]
  26. Mora, P. 1987 Nonlinear two-dimensional elastic inversion of multioffset seismic data. Geophysics52, no. 9: 1211–28.
    [Google Scholar]
  27. Nocedal, J., and S.J. Wright 2006Numerical optimization. New York: Springer Science & Business Media.
  28. Pratt, R.G. 1990 Frequency-domain elastic wave modeling by finite differences: A tool for crosshole seismic imaging. Geophysics55, no. 5: 626–32.
    [Google Scholar]
  29. Pratt, R.G. 1999 Seismic waveform inversion in the frequency domain, part 1: theory and verification in a physical scale model. Geophysics64, no. 3: 888–901.
    [Google Scholar]
  30. Pratt, R.G., C. Shin, and G.J. Hick 1998 Gauss–Newton and full Newton methods in frequency–space seismic waveform inversion. Geophysical Journal International133, no. 2: 341–62.
    [Google Scholar]
  31. Saad, Y. 2003Iterative methods for sparse linear systems. Society for Industrial and Applied Mathematics.
  32. Saad, Y., and M.H. Schultz 1986 GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM Journal on Scientific Computing7, no. 3: 856–69.
    [Google Scholar]
  33. Shin, C., S. Jang, and D.J. Min 2001 Improved amplitude preservation for prestack depth migration by inverse scattering theory. Geophysical Prospecting49, no. 5: 592–606.
    [Google Scholar]
  34. Sleijpen, G.L., and D.R. Fokkema 1993 BiCGSTAB(l) for linear equations involving unsymmetric matrices with complex spectrum. Electronic Transactions on Numerical Analysis1, no. 11: 2000.
    [Google Scholar]
  35. Song, Z.M., and P.R. Williamson 1995 Frequency-domain acoustic-wave modeling and inversion of crosshole data: part I 2.5-D modeling method. Geophysics60, no. 3: 784–95.
    [Google Scholar]
  36. Tape, C., Q. Liu, and J. Tromp 2007 Finite-frequency tomography using adjoint methods – methodology and examples using membrane surface waves. Geophysical Journal International168, no. 3: 1105–29.
    [Google Scholar]
  37. Tarantola, A. 1984 Inversion of seismic reflection data in the acoustic approximation. Geophysics49, no. 8: 1259–66.
    [Google Scholar]
  38. Tarantola, A. 1986 A strategy for nonlinear elastic inversion of seismic reflection data. Geophysics51, no. 10: 1893–903.
    [Google Scholar]
  39. Tong, P., D.H. Yang, B.L. Hua, and M.X. Wang 2013 A high-order stereo-modeling method for solving wave equations. Bulletin of the Seismological Society of America103, no. 2A: 811–33.
    [Google Scholar]
  40. Versteeg, R. 1994 The Marmousi experience: Velocity model determination on a synthetic complex data set. The Leading Edge13, no. 9: 927–36.
    [Google Scholar]
  41. Wang, M.X., D.H. Yang, and G.J. Song 2012 Semi-analytical solutions and numerical simulations of 2 D SH wave equation. Diqiu Wuli Xuebao55, no. 3: 914–24.
    [Google Scholar]
  42. Yang, D.H., M. Lu, R.S. Wu, and J.M. Peng 2004 An optimal nearly analytic discrete method for 2D acoustic and elastic wave equations. Bulletin of the Seismological Society of America94, no. 5: 1982–92.
    [Google Scholar]
  43. Yang, D.H., J.M. Peng, M. Lu, and T. Terlaky 2006 Optimal nearly analytic discrete approximation to the scalar wave equation. Bulletin of the Seismological Society of America96, no. 3: 1114–30.
    [Google Scholar]
  44. Yang, D.H., J.W. Teng, Z.J. Zhang, and E.R. Liu 2003 A nearly analytic discrete method for acoustic and elastic wave equations in anisotropic media. Bulletin of the Seismological Society of America93, no. 2: 882–890.
    [Google Scholar]
  45. Zhang, J.H., and Z.X. Yao 2012 Optimized finite-difference operator for broadband seismic wave modeling. Geophysics78, no. 1: A13–A18.
    [Google Scholar]
/content/journals/10.1080/08123985.2019.1699786
Loading
/content/journals/10.1080/08123985.2019.1699786
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): 2D modelling; finite difference; frequency-domain; full waveform; 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