Viscoacoustic reverse time migration with stable and effective two-way attenuation compensation

Abstract

Intrinsic attenuation in seismic wave propagation leads to amplitude dissipation and phase dispersion in seismic data. The attenuation (∼1/ Q) effect is a common cause for inaccurate image locations and dimmed image energy in reverse time migration (RTM). Conventional viscoacoustic or viscoelastic RTM methods implement attenuation compensation by reversing the sign of the amplitude loss term and keeping the dispersion term unchanged for forward and backward wavefield extrapolation, which requires the decoupled viscoacoustic wave equation and gives rise to the numerical instability issue. To address these problems, we develop a stable and effective two-way attenuation-compensated viscoacoustic RTM method. Our method uses a new two-way normalized crosscorrelation imaging condition to compensate the attenuation effect. In the new imaging condition, only the attenuated forward wavefield of sources and the viscoacoustic Green’s functions of receivers are required, which eliminates the instability source arising from the operation of reversing the sign of the dissipation term in the viscoacoustic wave equation. Because no modifications are made to the viscoacoustic wave equation, the new imaging condition can be calculated without any additional computational cost and can be applied to decoupled and coupled viscoacoustic wave equation-based RTMs, which is more flexible than other attenuation compensation strategies. Two synthetic tests and one field data example are presented to demonstrate the stability and effectiveness of our method.