Reverse time migration (RTM) for attenuating media should take amplitude compensation and phase correction into consideration. However, attenuation compensation during seismic propagation suffers from numerical instability because of the boosted high-frequency ambient noise. We have developed a novel adaptive stabilization method for [Formula: see text]-compensated RTM ([Formula: see text]-RTM), which exhibits superior properties of time variance and [Formula: see text] dependence over conventional low-pass filtering-based method. We derive the stabilization operator by first analytically deriving [Formula: see text]-space Green’s functions for a constant-[Formula: see text] wave equation with decoupled fractional Laplacians and its compensated equation. The time propagator of Green’s function for the viscoacoustic wave equation decreases exponentially, whereas that of the compensated equation is exponentially divergent at a high wavenumber, and it is not stable after the wave is extrapolated for a long time. Therefore, the Green’s functions theoretically explain how the numerical instability existing in [Formula: see text]-RTM arises and shed light on how to overcome this problem pertinently. The stabilization factor required in the proposed method can be explicitly identified by the specified gain limit according to an empirical formula. The [Formula: see text]-RTM results for noise-free data using low-pass filtering and adaptive stabilization are compared over a simple five-layer model and the BP gas chimney model to verify the superiority of the proposed approach in terms of fidelity and stability. The [Formula: see text]-RTM result for noisy data from the BP gas chimney model further demonstrates that our method enjoys a better antinoise performance and helps significantly to enhance the resolution of seismic images.
Green’s Functions for the Wave Equation