A comparative study of data- and image- domain LSRTM under velocity-impedance parametrization

TL;DR

This study compares data- and image-domain LSRTM for seismic imaging, revealing that data-domain approaches better reconstruct reflectivity, while image-domain methods face limitations due to Hessian sampling issues, with implications for multiparameter inversion.

Contribution

The paper provides a comprehensive comparison of multiparameter LSRTM in data and image domains, highlighting the advantages of data-domain inversion and the challenges in image-domain multiparameter reconstruction.

Findings

Data-domain LSRTM yields superior reflectivity reconstruction.

Image-domain multiparameter inversion suffers from Hessian sampling limitations.

Monoparameter impedance inversion performs well in the image domain.

Abstract

Least-squares reverse time migration (LSRTM) is one of the classic seismic imaging methods to reconstruct model perturbations within a known reference medium. It can be computed in either data or image domain using different methods by solving a linear inverse problem, whereas a careful comparison analysis of them is lacking in the literature. In this article, we present a comparative study for multiparameter LSRTM in data- and image- domain in the framework of SMIwiz open software. Different from conventional LSRTM for recovering only velocity perturbation with variable density, we focus on simultaneous reconstruction of velocity and impedance perturbations after logorithmic scaling, using the first-order velocity-pressure formulation of acoustic wave equation. The first 3D data-domain LSRTM example has been performed to validate our implementation, involving expensive repetition of Born modelling and migration over a number of iterations. As a more cost-effective alternative, the image-domain LSRTM is implemented using point spread function (PSF) and nonstationary deblurring filter. Dramatic disctinctions between data and image domain methods are discovered with 2D Marmousi test: (1) The data-domain multiparameter inversion provides much better reconstruction of reflectivity images than image-domain approaches, thanks to the complete use of Hessian in Krylov space; (2) The poor multiparameter image-domain inversion highlights the limitation of incomplete Hessian sampling and strong parameter crosstalks, making it difficult to work in practice; (3) In contrast, monoparameter image-domain inversion for seismic impedance is found to work well. These observations have been further validated on Viking Graben Line 12 dataset.