3-D Q-Compensated Image-Domain Least-Squares Reverse Time Migration Through Point Spread Functions

Least-squares reverse time migration (LSRTM) has the potential to retrieve a high-resolution subsurface image. However, the standard acoustic LSRTM approach may produce a blurred image, if directly applying it to attenuated seismic recordings. In this letter, we developed a novel 3D Q-compensated image-domain LSRTM approach, denoted as Q-IDLSRTM. The Hessian matrix in the proposed approach is efficiently estimated from the point spread functions (PSFs) which are calculated by a combination of viscoacoustic Born modeling and reverse time migration (RTM) based on the generalized standard linear solid (GSLS) wave equation. The major advantage of the proposed image-domain inversion is that it is much faster than data-domain inversion. The L1 norm constraint and total variation (TV) regularization are used to produce a sparse solution and maintain the structural continuity of the inverted image. We determine the effectiveness of the proposed approach with a part of the 3D Overthrust model and the resulting images demonstrate the ability of our approach to image subsurface structures with enhanced resolution and balanced amplitude relative to the RTM image and inverted image from the acoustic image-domain LSRTM approach.