Abstract: The real geological medium is close to viscoelastic media in which seismic wave propagate with dispersion and attenuation. It is important to compensate the unwanted viscous effects in seismic processing. It is more accurate and physically more consistent to mitigate these effects in a wave-equation-based prestack depth migration. Historically, reverse time migration (RTM) based on directly solving the two-way wave equation has provided a superior way to image complex geologic regions. However, instability usually arises when considering compensation for absorption. Most researchers conduct high frequency filtering in wavenumber domain before or during wave field extrapolate in RTM to ensure stability. In this paper, we use viscoacoustic wave equation derived by Bai et al. (2013) to do Q-RTM and stabilize extrapolate by adding a regularization term. Compared with direct filtering, the regularization parameter can be space-varying. So this is suitable for severely variational regions. And we also find that source normalized cross-correlation imaging condition is more suitable in Q-RTM.