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.

What problem does this paper attempt to address?

This paper attempts to solve the problems encountered in reverse - time migration (RTM) of quasi - P - waves in tilted transversely isotropic (TTI) media. Specifically, the paper focuses on how to reduce the computational cost while maintaining high accuracy and effectively avoid the generation of S - wave artifacts, thereby improving the quality of seismic imaging. The following are the main research objectives of the paper: 1. **Improve the accuracy of the quasi - P - wave equation**: In traditional methods, when dealing with the quasi - P - wave equation in TTI media, due to the involvement of fractional - order pseudo - differential operators, the calculation is complex and S - wave artifacts are easily generated. The paper proposes a new pure quasi - P - wave equation. By introducing a generalized fractional convolution template to approximate the fractional - order pseudo - differential operator, the accuracy of the equation is improved. 2. **Reduce the computational cost**: The new method not only improves the computational accuracy but also reduces the computational cost through local convolution and the finite - difference method. This makes it more efficient when dealing with large - scale seismic data. 3. **Eliminate S - wave artifacts**: The traditional coupled quasi - P - wave equation will generate S - wave artifacts when simulating the wave field in TTI media, which affects the imaging quality. The method proposed in the paper can fundamentally avoid these artifacts and ensure the accuracy of the imaging results. 4. **Verify the effectiveness of the method**: Through numerical experiments, the paper shows the performance of the proposed method in different models (such as the homogeneous model, the nappe structure TTI model, and the 2007 BP TTI model), proving its feasibility and effectiveness in practical applications. ### Main content of the paper 1. **Introduction**: - Introduces the importance of reverse - time migration (RTM) in seismic imaging and its advantages in dealing with complex geological structures. - Emphasizes that ignoring the influence of anisotropy in strongly anisotropic regions (such as subsalt sediments and shale reservoirs) will lead to blurry imaging and positional deviation. 2. **Method**: - Proposes a new pure quasi - P - wave equation, which is decomposed into an elliptic anisotropy equation and a fractional - order pseudo - differential correction term. - Uses a generalized fractional convolution template to approximate the pseudo - differential term in the spatial domain and determines the convolution template by solving a nonlinear inverse problem. - Derives in detail the new pure quasi - P - wave equation and its forms in the time and spatial domains. 3. **Numerical implementation**: - Prepares a library of 932,331 fractional convolution templates of size 11×11. - Describes the specific steps of pure quasi - P - wave field modeling in the TTI model, including loading the convolution template, convolution operation, finite - difference calculation, and wave field update. 4. **Numerical experiments**: - Verifies the effectiveness of the proposed method through four numerical experiments (homogeneous model, nappe structure TTI model, 2007 BP TTI model). - Shows the superior performance of the proposed method compared to traditional methods in wave field simulation and RTM imaging in different models. ### Conclusion The pure quasi - P - wave equation and its approximation method proposed in the paper can effectively avoid S - wave artifacts, improve imaging quality and computational efficiency when dealing with wave field modeling and reverse - time migration in TTI media. The numerical experiment results further verify the feasibility and effectiveness of this method.

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