A non-balanced staggered-grid finite-difference scheme for the first-order elastic wave-equation modeling

Wenquan Liang,Yanfei Wang,Ursula Iturrarán-Viveros

DOI: https://doi.org/10.48550/arXiv.2012.07143

2020-12-14

Abstract:We introduce an efficient and accurate staggered-grid finite-difference (SGFD) method to solve the two-dimensional elastic wave equation. We use a coupled first-order stress-velocity formulation. In the standard implementation of SGFD method the same SGFD operator is used to approximate the spatial derivatives. However, we propose a numerical method based on mixed SGFD operators which happen to be more efficient with similar accuracy in comparison to uniform SGFD operator. We refer the proposed method as the non-balanced SGFD numerical scheme which means combining high-order SGFD operators with second-order SGFD operators. A very care attention is directed at the derivation of the SGFD operator coefficients. The correctness of proposed scheme is proven by dispersion analysis. Through SGFD modeling examples, we verify/demonstrate that the proposed non-balanced operator offers a similar level of accuracy with a cheaper computation cost compared to the more expensive balanced SGFD method.

Numerical Analysis,Geophysics

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to improve the efficiency and accuracy of the numerical simulation of the two - dimensional elastic wave equation while reducing the computational cost. Specifically:

1. **Problems of existing methods**:
- The traditional Staggered - Grid Finite - Difference (SGFD) method has the problem of numerical dispersion when simulating elastic wave propagation.
- Using high - order SGFD operators can improve accuracy, but it will increase computational cost and memory requirements.

2. **The proposed new method**:
- The paper proposes a Non - balanced SGFD scheme, which calculates spatial derivatives by combining SGFD operators of different lengths (long operators and short operators).
- This method aims to significantly reduce the computational cost while maintaining a similar accuracy to the traditional uniform SGFD method.

3. **Specific implementation**:
- **Non - balanced SGFD scheme**: Use a mixed SGFD operator, that is, use high - order and second - order SGFD operators simultaneously in the same time step to calculate spatial derivatives.
- **Coefficient calculation**: Through carefully designing SGFD coefficients, ensure that the error introduced by the short operator is offset by the long operator, thereby reducing numerical dispersion and improving simulation accuracy.
- **Stability analysis**: Derive the stability conditions of the non - balanced SGFD scheme and verify its effectiveness through numerical experiments.

4. **Verification and comparison**:
- The paper verifies the effectiveness of the non - balanced SGFD scheme through numerical experiments, including simulating uniform models and complex salt body models.
- The experimental results show that the non - balanced SGFD scheme can reduce the computation time by about 40% while maintaining high accuracy.

In summary, the main goal of this paper is to provide a more efficient and accurate numerical solution method for the elastic wave equation by introducing the non - balanced SGFD scheme, in order to address the problems of high computational cost and numerical dispersion in traditional methods.