Abstract:Conventional finite-difference methods will encounter strong numerical dispersion on coarse grids or when higher frequencies are used. The nearly analytic discrete method (NADM) and its improved version, developed recently, can suppress numerical dispersion efficiently. However, the NADM is imperfect, especially in long-time seismic wave propagation simulation. To overcome this limitation, by minimizing the energy-error function, we have developed a modified optimal nearly analytical discretized (MNAD) method. The stencil of MNAD for the 2D wave equation is a novel diamond stencil. MNAD has two major advantages: The first one is that MNAD can suppress numerical dispersion effectively on coarse grids, which significantly improved computational efficiency and reduced memory demands. The second advantage is that the energy error of MNAD is much smaller after long-time simulation. The numerical dispersion analysis shows that the maximum phase velocity error was 5.92% even if only two sampling points were adopted in each minimum wavelength. To simulate a wavefield without visible numerical dispersion, the computational speed of MNAD, measured by CPU time, was approximately 4.32 times and 1.43 times comparing with the fourth-order Lax-Wendroff correction (LWC) method and the optimal nearly analytical discretized method (ONADM), respectively. MNAD also shows good numerical stability. Its CFL condition increased 24.3% comparing with that of ONADM, from 0.523 to 0.650. The total energy error was less than 1.5% after 300-s simulation whereas the error of other numerical schemes, such as the fourth-order and eighth-order LWC, etc., was up to more than 10% under the same computational parameters. Numerical results showed that our MNAD method performed well in computational efficiency, simulation accuracy, and numerical stability, as well as provided a useful tool for large-scale, long-time seismic wave propagation simulation.

Optimally Accurate Nearly Analytic Discrete Scheme for Wave-Field Simulation in 3D Anisotropic Media

An Optimal Nearly Analytic Discrete Method for 2d Acoustic and Elastic Wave Equations

A Nearly Analytical Discrete Method For Wave-Field Simulations In 2d Porous Media

A Nearly Analytic Discrete Method for Acoustic and Elastic Wave Equations in Anisotropic Media

An Improved Nearly Analytical Discrete Method: an Efficient Tool to Simulate the Seismic Response of 2-D Porous Structures

An Improved Optimal Nearly Analytical Discretized Method for 2D Scalar Wave Equation in Heterogeneous Media Based on the Modified Nearly Analytical Discrete Operator

Optimal Nearly Analytic Discrete Approximation to the Scalar Wave Equation

A Modified Nearly Analytic Discrete Method and Wavefield Simulations in Transversely Isotropic Media

Elastic Wave Modelling Method Based on the Displacement–velocity Fields: an Improving Nearly Analytic Discrete Approximation

A modified NAD algorithm with minimum numerical dispersion for simulation of anisotropic wave propagation

An Optimal Nearly Analytic Discrete-Weighted Runge-Kutta Discontinuous Galerkin Hybrid Method for Acoustic Wavefield Modeling

An Efficient and high accuracy finite-difference scheme for the acoustic wave equation in 3D heterogeneous media

A WNAD method for seismic stress-field modeling in heterogeneous media

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

FDwave3D: a MATLAB solver for the 3D anisotropic wave equation using the finite-difference method

Parallel WNAD algorithm for solving 3D elastic equation and its wavefield simulations in TI media

Simulation of Acoustic Wavefields in Heterogeneous Media: A Robust Method for Automatic Suppression of Numerical Dispersion

A 2.5d Finite-Difference Algorithm for Elastic Wave Modeling Using Near-Optimal Quadratures

A modified time-space optimized symplectic method for solving 3D elastic wave equations

Discontinuous Galerkin Modeling of 3D Arbitrary Anisotropic Q

An NAD Scheme with Wavenumber Error Optimized for 2D Scalar Wave Equation