Compressed-sensing application - Pre-stack kirchhoff migration

Least-squares migration is a linearized form of waveform inversion that aims to enhance the spatial resolution of the subsurface reflectivity distribution and reduce the migration artifacts due to limited recording aperture, coarse sampling of sources and receivers, and low subsurface illumination. Least-squares migration, however, due to the nature of its minimization process, tends to produce smoothed and dispersed versions of the reflectivity of the subsurface. Assuming that the subsurface reflectivity distribution is sparse, we propose the addition of a non-quadratic L1-norm penalty term on the model space in the objective function. This aims to preserve the sparse nature of the subsurface reflectivity series and enhance resolution. We further use a compressed-sensing algorithm to solve the linear system, which utilizes the sparsity assumption to produce highly resolved migrated images. Thus, the Kirchhoff migration implementation is formulated as a Basis Pursuit denoise (BPDN) problem to obtain the sparse reflectivity model. Applications on synthetic data show that reflectivity models obtained using this compressed-sensing algorithm are highly accurate with optimal resolution.