High-Accuracy Simulation of Rayleigh Waves Using Fractional Viscoelastic Wave Equation

Abstract

The propagation of Rayleigh waves is usually accompanied by dispersion, which becomes more complex with inherent attenuation. The accurate simulation of Rayleigh waves in attenuation media is crucial for understanding wave mechanisms, layer thickness identification, and parameter inversion. Although the vacuum formalism or stress image method (SIM) combined with the generalized standard linear solid (GSLS) is widely used to implement the numerical simulation of Rayleigh waves in attenuation media, this type of method still has its limitations. First, the GSLS model cannot split the velocity dispersion and amplitude attenuation term, thus limiting its application in the Q-compensated reverse time migration/full waveform inversion. In addition, GSLS-model-based wave equation is usually numerically solved using staggered-grid finite-difference (SGFD) method, which may result in the numerical dispersion due to the harsh stability condition and poses complexity and computational burden. To overcome these issues, we propose a high-accuracy Rayleigh-waves simulation scheme that involves the integration of the fractional viscoelastic wave equation and vacuum formalism. The proposed scheme not only decouples the amplitude attenuation and velocity dispersion but also significantly suppresses the numerical dispersion of Rayleigh waves under the same grid sizes. We first use a homogeneous elastic model to demonstrate the accuracy in comparison with the analytical solutions, and the correctness for a viscoelastic half-space model is verified by comparing the phase velocities with the dispersive images generated by the phase shift transformation. We then simulate several two-dimensional synthetic models to analyze the effectiveness and applicability of the proposed method. The results show that the proposed method uses twice as many spatial step sizes and takes 0.6 times that of the GSLS method (solved by the SGFD method) when achieved at 95% accuracy.

Quasi-P-Wave Reverse Time Migration in TTI Media with a Generalized Fractional Convolution Stencil

“…[19] proposed a pure quasi-P-wave equation that could fundamentally avoid S-wave artifacts. However, this pure quasi-P-wave equation involved a space-fractional pseudo-differential operator that depended on the anisotropy parameters, making it unsuitable for direct resolution using conventional solvers such as the finitedifference method and the pseudo-spectral method [20][21][22]. In order to solve the pure quasi-P-wave equation efficiently, Refs.…”