Anisotropic VSP reverse-time migration with optimal pure aeoustic wave and complex wavefield separation

Ignoring anisotropy characteristic of subsurface media may lead to misplaced images and low resolution of the target for the reverse-time migration (RTM). The mature anisotropic RTM methods are mainly based on the pseudoacoustic wave approximation. Although these schemes have high computational efficiency, most of pseudo-acoustic wave equations (PWEs) inevitably encounter SV-wave artifacts or instability for anisotropic modeling and imaging. To improve the anisotropic RTM quality, we develop a combination of optimal pure acoustic wave and complex wavefield separation to conduct anisotropic RTM for both surface and vertical seismic profiling (VSP) acquisition geometries. Among the proposed scheme, we derive an optimal pure acoustic wave dispersion relation, and solve the corresponding wave equation by incorporating finite-difference and Poisson solver. The modified equation can remove SV-wave artifacts and instability of PWEs. Wavefield separation approach can choose desired wavefield components along different directions to carry out the final imaging, which can effectively suppress low-frequency imaging noise. Moreover, the hybrid absorbing boundary condition is adopted to suppress artificial boundary reflections during wavefield extrapolation. Basic theory and modeling examples demonstrate that the developed schemes can generate RTM results with high accuracy.