2D Laplace-Fourier domain acoustic wave equation modeling with an optimal finite-difference method

Laplace-Fourier (L-F) domain finite-difference (FD) forward modeling is an important foundation for L-F domain full-waveform inversion (FWI). An optimal modeling method can improve the efficiency and accuracy of FWI. A flexible FD stencil, which requires pairing and centrosymmetricity of the involved gridpoints, is used on the basis of the 2D L-F domain acoustic wave equation. The L-F domain numerical dispersion analysis is then performed by minimizing the phase error of the normalized numerical phase and attenuation propagation velocities to obtain the optimization coefficients. An optimal FD forward modeling method is finally developed for the L-F domain acoustic wave equation and applied to the traditional standard 9-point scheme and 7- and 9-point schemes, where the latter two schemes are used in discontinuous-grid FD modeling. Numerical experiments show that the optimal L-F domain FD modeling method not only has high accuracy but can also be applied to equal and unequal directional sampling intervals and discontinuous-grid FD modeling to reduce computational cost.