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

Abstract

In seismic modeling and reverse time migration (RTM), incorporating anisotropy is crucial for accurate wavefield modeling and high-quality images. Due to the trade-off between computational cost and simulation accuracy, the pure quasi-P-wave equation has good accuracy to describe wave propagation in tilted transverse isotropic (TTI) media. However, it involves a fractional pseudo-differential operator that depends on the anisotropy parameters, making it unsuitable for resolution using conventional solvers for fractional operators. To address this issue, we propose a novel pure quasi-P-wave equation with a generalized fractional convolution operator in TTI media. First, we decompose the conventional pure quasi-P-wave equation into an elliptical anisotropy equation and a fractional pseudo-differential correction term. Then, we use a generalized fractional convolution stencil to approximate the spatial-domain pseudo-differential term through the solution of an inverse problem. The proposed approximation method is accurate, and the wavefield modeling method based on it also accurately describes quasi-P-wave propagation in TTI media. Moreover, it only increases the computational cost for calculating mixed partial derivatives compared to those in vertical transverse isotropic (VTI) media. Finally, the proposed wavefield modeling method is utilized in RTM to correct the anisotropic effects in seismic imaging. Numerical RTM experiments demonstrate the flexibility and viability of the proposed method.