Time-domain numerical solution of the 3D acoustic wave equation is essential part of seismic imaging procedures such as FWI and RTM. Due to numerical dispersion error, we need to design accurate numerical schemes. Although high-order schemes are effective in treating the problem, but are not computationally efficient. In this study we suggest an accurate method for solution of 3D acoustic wave equation based on second-order implicit finite-difference scheme. To achieve this goal, we combine three different version of Laplacian operators along with mass matrix un-lumping. By minimizing dispersion error, we obtain optimized coefficients corresponding to each CFL number. The accurate second-order scheme leads to a compact stencil with 27 nodal points. According to the numerical tests this scheme can simulate waves propagation in different direction accurately for different CFL numbers. Due to compact characteristics of the used 27-point scheme it requires minimum computational cost comparing to the high-order ones.


Article metrics loading...

Loading full text...

Full text loading...


  1. Chen, J.B
    [2014] A 27-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method. Geophysical Prospecting, 62, 258–277.
    [Google Scholar]
  2. Hustedt, B., Operto, S. and Virieux, J.
    [2004] Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modelling. Geophysical Journal International, 157, 1269–1296.
    [Google Scholar]
  3. Jo, C.H, Shin, C. and Suh, J.H
    [1996] An optimal 9-point, finite-difference, frequency-space, 2-Dscalar wave extrapolator. Geophysics, 61(2), 529 –537.
    [Google Scholar]
  4. Marfurt, K.
    , [1984] Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations. Geophysics, 49, 533–549.
    [Google Scholar]
  5. Operto, S., Virieux, J., Amestoy, P., L'Excellent, J.-Y., Giraud, L. and Ali, H.B.H.
    [2007] 3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study. Geophysics, 72, no. 5, SM195–SM211.
    [Google Scholar]
  6. Park, H., Park, Y., Lee, J. and Shin, C.
    [2017] A 27-point finite difference method using a symmetric perfectly matched layer for the 3D Helmholtz equation. SEG International Exposition and 87th Annual Meeting.
    [Google Scholar]
  7. Shin, C. and Sohn, H.
    [1998] A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator. Geophysics, 63(1), 289–296.
    [Google Scholar]

Data & Media loading...

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