A compact high order Alternating Direction Implicit method for three-dimensional acoustic wave equation with variable coefficient

Abstract: Efficient and accurate numerical simulation of seismic wave propagation is important in various Geophysical applications such as seismic full waveform inversion (FWI) problem. However, due to the large size of the physical domain and requirement on low numerical dispersion, many existing numerical methods are inefficient for numerical modelling of seismic wave propagation in a heterogeneous media. Despite the great efforts that have been devoted during the past decades, it still remains a challenging task in the development of efficient and accurate finite difference method for the multi-dimensional acoustic wave equation with variable velocity. In this paper, we proposed a Pad\'{e} approximation based finite difference scheme for solving the acoustic wave equation in three-dimensional heterogeneous media. The new method is obtained by combining the Pad\'{e} approximation and a novel algebraic manipulation. The efficiency of the new algorithm is further improved through the Alternative Directional Implicit (ADI) method. The stability of the new algorithm has been theoretically proved by the energy method. The new method is conditionally stable with a better Courant - Friedrichs - Lewy condition (CFL) condition, which has been verified numerically. Extensive numerical examples have been solved, which demonstrated that the new method is accurate, efficient and stable.