Abstract:The time-space-domain finite-difference method (FDM) is widely used in forward modeling of wave equations. Conventional explicit FDMs (EFDMs) (with high order in space and the second order in time) would result in apparent temporal and spatial dispersion for high frequencies and large time steps. Moreover, the saturation effect of Taylor expansion seriously restricts the improvement of bandwidth coverage and efficiency of EFDMs. We have developed a variable-order optimization scheme of the implicit FDM (IFDM) to improve the computational efficiency and numerical accuracy of forward modeling. Then, we applied time-dispersion transforms (TDTs) to filter out the temporal dispersion generated by the second-order temporal approximations. Our method greatly alleviates the saturation effect of the high-order spatial finite-difference (FD) operators. Dispersion analysis indicates that the optimized coefficients of the IFDM based on the Remez algorithm can achieve the widest bandwidth (close to the Nyquist wavenumber), which corresponds to the shortest length of the spatial FD operator under a given error threshold. Our method has great potential to approach the highest spectral accuracy but with minimal increase in computational cost. Numerical experiments indicate that the combination of the variable-order optimization of IFDM and TDTs can significantly reduce numerical dispersion. Compared with traditional methods, our scheme is more advantageous for numerical simulation of large-scale geologic models because it has the least amount of calculation burden under the same accuracy requirements.

A Robust Method of Computing Finite Difference Coefficients Based on Vandermonde Matrix

A Novel Matrix-Free Finite Element Method for Time-Harmonic Maxwell’s Equations

A Compact High-Order Finite-Difference Method with Optimized Coefficients for 2D Acoustic Wave Equation

Splitting Algorithms for the High-Order Compact Finite-Difference Schemes in Wave-Equation Modeling

Wavefield simulation of acoustic VTI wave equation based on adaptive-coefficient finite-difference frequency-domain method

Finite-difference Modelling in Two-Dimensional Anisotropic Media Using a Flux-Corrected Transport Technique

A finite‐difference‐based technique for numerically efficient computation of capacitance matrices for 2‐dimensional multi‐conductor problems involving thin dielectric coating

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

Approximating and preconditioning the stiffness matrix in the GoFD approximation of the fractional Laplacian

A stabilized radial basis-finite difference (RBF-FD) method with hybrid kernels

Efficient and Stable Finite Difference Modelling of Acoustic Wave Propagation in Variable-density Media

Variable-order optimal implicit finite-difference schemes for explicit time-marching solutions to wave equations

Seismic modeling of scalar seismic wave propagation with finitedifference scheme in frequency-space domain.

Adaptive Finite Differencing in High Accuracy Electronic Structure Calculations

An Optimized Time-Space-domain Finite Difference Method with Piecewise Constant Interpolation Coefficients for Scalar Wave Propagation

Optimal finite-difference schemes for elastic wave based on improved cosine-combined window function

An introduction to finite element methods for inverse coefficient problems in elliptic PDEs

An accurate and efficient multiscale finite-difference frequency-domain method for the scalar Helmholtz equation

The Convolutional Differentiator Method for Numerical Modelling of Acoustic and Elastic Wavefields

Numerical simulation of acoustic waves propagation by finite element method based on optimized matrices