Simplified implicit finite-difference method of spatial derivative using explicit schemes with optimized constant coefficients based on L1 norm

Abstract

The explicit finite-difference (EFD) method is widely used in numerical simulation of seismic wave propagation to approximate spatial derivatives. However, the traditional and optimized high-order EFD methods suffer from the saturation effect, which seriously restricts the improvement of numerical accuracy. In contrast, the implicit FD (IFD) method approximates the spatial derivatives in the form of rational functions and thus can obtain much higher numerical accuracy with relatively low orders; however, its computational cost is expensive due to the need to invert a multidiagonal matrix. We derive an explicit strategy for the IFD method to reduce the computational cost by constructing the IFD method with the discrete Fourier matrix; then, we transform the inversion of the multidiagonal matrix into an explicit matrix multiplication; next, we construct an objective function based on the [Formula: see text] norm to reduce approximation error of the IFD method. This explicit strategy of the IFD method can avoid inverting the multidiagonal matrix, thus improving the computational efficiency. This constant coefficient optimization method reduces the approximation error in the medium-wavenumber range at the cost of tolerable deviation (smaller than 0.0001) in the low-wavenumber range. For the 2D Marmousi model, the root-mean-square error of the numerical results obtained by this method is one-fifth that of the traditional IFD method with the same order (i.e., 5/3) and one-third that of the traditional EFD method with much higher orders (i.e., 72). The significant reduction of numerical error makes the developed method promising for numerical simulation in large-scale models, especially for long-time simulations.