ABSTRACT
By solving a linear inverse problem, least-squares migration (LSM) can provide higher-precision images of complex subsurface structures than traditional adjoint-operator-based migration. LSM can be conducted in the data domain and the image domain (IDLSM). IDLSM in acoustic media has been extensively studied, commonly using single-parameter point-spread function (PSF) deconvolution. For high-accuracy imaging of complex media, we develop an image-domain elastic least-squares imaging method using multiparameter PSF deconvolution for vertically transversely isotropic (VTI) media. The multiparameter PSFs, including P-P, P-S, S-P, and S-S components, are calculated on a sparse grid using elastic VTI wave equation Born modeling and migration. The traditional migration result is then decomposed into individual local results, each corresponding to the size of a single PSF, through unit partitioning. Finally, a set of filters is derived from four elastic VTI PSF components to approximate the Hessian inverse, which are subsequently assembled to reconstruct an entire deconvolved image. The numerical experiments confirm that our method can significantly decrease multiparameter crosstalk artifacts and improve spatial resolution and amplitude fidelity.