Cost-effective elastic-wave modeling is the key to practical elastic reverse time migration and full-waveform inversion implementations. We have developed an efficient elastic pseudo-analytical finite-difference (PAFD) scheme for elastic-wave extrapolation. The elastic PAFD scheme is based on a modified pseudo-spectral method, k-space method, in which a pseudo-analytical operator is used to ensure the high accuracy of elastic-wave extrapolation. However, the k-space method is motivated for a pure wave mode, and thus its application in coupled first-order elastic-wave equations may cause the elastic pseudo-analytical operators to suffer from crosstalk between the P- and S-wavefields. The approaches presented attempt to overcome these shortcomings by introducing two improvements to achieve the goal. This is done, first, by performing a predictor-corrector strategy in first-order elastic-wave equations to eliminate those errors during wave extrapolation. Considering the massive computational cost in the spectral domain, we have developed an efficient elastic PAFD implementation, in which an innovative model-adaptive finite-difference coefficient-predicted scheme is provided to reduce the computational cost of elastic pseudo-analytical operator differencing. Dispersion analysis demonstrates the flexibility with varying velocity and superior performance of our PAFD scheme for spatial and temporal dispersion suppression than the existing Taylor-expansion-based scheme. Under the same simulation parameters, several numerical examples prove that the elastic PAFD scheme can provide more accurate simulation results, whereas the conventional scheme suffers from spatial or temporal dispersion errors, even in complex heterogeneous media.
C-N Difference Schemes for Dissipative Symmetric Regularized Long Wave Equations with Damping Term