Multidirectional-vector-based elastic reverse time migration and angle-domain common-image gathers with approximate wavefield decomposition of P- and S-waves

Elastic reverse time migration (E-RTM) has limitations when the migration velocities contain strong contrasts. First, the traditional scheme of P/S-wave mode separation is based on Helmholtz's equations, which ignore the conversion between P-and S-waves at the current separation time. Thus, it contains an implicit assumption of the constant shear modulus and requires smoothing the heterogeneous model to approximately satisfy a locally constant condition. Second, the vector-based imaging condition needs to use the reflection-image normal, and it also cannot give the correct polarity of the PP image in all possible conditions. Third, the angle-domain common-image gathers (ADCIGs) calculated using the Poynting vectors (PVs) do not consider the wave interferences that happen at each reflector. Therefore, smooth models are often used for E-RTM. We relax this condition by proposing an improved data flow that involves three new contributions. The first contribution is an improved system of P/S-wave mode separation that considers the converted wave generated at the current time, and thus it does not require the constant-shear-modulus assumption. The second contribution is the new elastic imaging conditions based on multidirectional vectors; they can give the correct image polarity in all possible conditions without knowledge of the reflection-image normal. The third contribution is two methods to calculate multidirectional propagation vectors (PRVs) for RTM images and ADCIGs: One is the elastic multidirectional PV, and the other uses the sign of wavenumber-over-frequency (k/omega) ratio obtained from an amplitude-preserved approximate-propagation-angle-based wavefield decomposition to convert the particle velocities into multidirectional PRVs. The robustness of the improved data flow is determined by several 2D numerical examples. Extension of the schemes into 3D and amplitude-preserved imaging conditions is also possible.