Viscoacoustic reverse time migration with a robust space-wavenumber domain attenuation compensation operator

Intrinsic attenuation gives rise to phase dispersion and amplitude loss during seismic wave propagation. Not correcting these effects in seismic imaging can result in inaccurate reflector locations. dimmed amplitudes, and degraded spatial resolution. In reverse time migration (RTM), attenuation compensation can be implemented by reversing the sign of the dissipation term and keeping the dispersion term unchanged for backward wave-field extrapolation. Although this Q-compensated RTM scheme can effectively correct attenuation effects, amplitude amplification during backpropagation might lead to numerical instabilities, especially for field data with strong high-frequency noise. To mitigate this problem, we have developed a robust space-wavenumber compensation operator and applied it to viscoacoustic RTM. By analyzing the dispersion-only and viscoacoustic Green's functions, we obtain an analytical solution for the attenuation compensation operator in a homogeneous medium. Because it is a time-frequency operator, to apply it directly in viscoacoustic RTM requires access to the extrapolated wavefields within a certain time window. To avoid storing the wavefields and improve the computational efficiency, we use an approximated dispersion relation and convert the time-frequency operator to an equivalent space-wavenumber operator, which allows us to implement attenuation compensation on the fly during wavefield extrapolation. The hybrid-domain property of the operator enables us to account for the wavenumber-dependent compensation. A similar strategy also can be applied to the migrated images as a poststack processing approach, which is more efficient than the prestack compensation. Two synthetic and one land field data set examples demonstrate the feasibility and adaptability of our method.