Least-squares reverse time migration (LSRTM) is a powerful tool in seeking broadband-wavenumber reflectivity images. It produces better images over reverse-time migration (RTM) at the expense of computational cost. The Hessian effect can be measured in the image domain with the point-spread function (PSF). Here, we try to measure the Hessian effect in the data domain with the so-called trace-spread function (TSF). The difference between PSF and TSF is that the former originates from LTL in the image domain while the latter from LLT in the data domain. By comparing the TSFs with their original corresponding traces (or beams), we can design adaptive matching filters for preconditioning to alleviate the Hessian effect. However, the full TSF matrix is expensive. In this article, we propose a multiscale solution, which first has a diagonal approximation to LLT in beams, and then handle the full submatrix composed of the one-beam traces using the Sherman-Morrison formula. The preconditioned beams are superimposed into a ``deblurred'' data for remigration. Through synthetic and real data examples, we see that: 1) single-step data-domain LSRTM can yield deblurred RTM images via adaptive matching filters and 2) the beam-by-beam consideration outperforms the trace-by-trace one.
|Original language||English (US)|
|Number of pages||7|
|Journal||IEEE Transactions on Geoscience and Remote Sensing|
|State||Published - Nov 13 2019|
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The authors would like to thank S. Operto for the inspiration of the beam-by-beam consideration. They would also like to thank G. Schuster and D. Peter for their helpful discussion.