Enhancing 3D migration images resolution: A fast and robust implicit point spread function deconvolution method based on wavenumber domain representation

Image-domain least-squares migration has been proven effective in enhancing the spatial resolution of migration images. However, due to the intrinsic complexity of the 3D problem, the conventional local-stationary deblurring method cannot achieve a good balance between computational efficiency and deblurring accuracy. The explicit point spread function (PSF) deconvolution method has been proposed to improve the resolution of migration images on a point-wise bias with substantial computational and storage costs, especially for the 3D PSF deconvolution. To achieve high-resolution imaging with reduced costs, we introduce an implicit PSF deconvolution method based on wavenumber domain representation to enhance the resolution of 3D migration images. Using a deep-learning optimizer, coupled with a corresponding modified hyperparameter, we can achieve a fast convergence for the solution of the 3D PSF deconvolution operator. A loss interpolation technique is introduced to obtain non-sensitive PSF deconvolution operators for dense interpolation to reduce high-frequency artifacts. 3D synthetic and large-scale field dataset results demonstrate that our approach can produce high-resolution images with reduced migration artifacts and balanced amplitude. Additionally, our proposed method outperforms existing methods with at least three times computational cost reduction and two orders of magnitude reduction in storage cost for the 3D PSF deconvolution. It demonstrates that the proposed approach can be a cost-effective and practical high-resolution imaging tool for large-scale 3D datasets. • Implicit PSF deconvolution enhances 3D migration image resolution effectively. • Deep-learning optimizer ensures fast convergence for 3D PSF deconvolution. • Loss interpolation reduces artifacts and improves dense PSF operator interpolation. • 3D tests show high-resolution images with fewer artifacts, balanced amplitudes, and 3x computational and 100x storage cost reductions.