Efficient reverse time migration method in tilted transversely isotropic media based on a pure pseudoacoustic wave equation

Abstract

In general, velocity anisotropy in shale media has been widely observed in laboratory and field work, which means that disregarding this characteristic can lead to inaccurate imaging locations when data are imaged with reverse time migration (RTM). Wavefields simulated with the conventional coupled pseudoacoustic wave equation may introduce S-wave noise and this equation is only valid in transversely isotropic media ([Formula: see text]). Certain decoupled qP-wave equations require the use of the pseudospectral method, which makes them computationally inefficient. To address these issues, we develop a new pure qP acoustic wave equation based on the acoustic assumption, which can be solved more efficiently using the finite-difference (FD) method. This equation can also be used in the forward-modeling process of RTM in tilted transversely isotropic (TTI) media. First, we perform a Taylor expansion of the root term in the pure qP-wave dispersion relation. This leads to an anisotropic dispersion relation that is decomposed into an elliptical anisotropic background factor and a circular correction factor. Second, we obtain the pure qP-wave equation in TTI media without a pseudodifferential operator. The new equation can be efficiently solved using FD methods and can be applied to RTM in TTI media with strong anisotropy. Our method indicates greater tolerance to numerical errors and is better suited for strong anisotropy, as compared with previously published methods. Numerical examples indicate the high kinematic and phase accuracy of our pure qP-wave equation along with its stability in TTI media characterized by ([Formula: see text]). By using a sag model and an overthrust TTI model, we determine the efficiency and accuracy of our TTI RTM.