The ℓ1−2-norm-regularized basis pursuit seismic inversion based on the exact Zoeppritz equation

Abstract

Sparsity constraints have been widely adopted in the regularization of ill-posed problems to obtain subsurface properties with sparseness features. However, the target parameters are generally not sparsely distributed, and sparsity constraints lead to results that are missing information. Smooth constraints (e.g., [Formula: see text] norm) lead to the insufficient resolution of the inversion results. To overcome this issue, an effective solution is to convert the target parameters to a sparse representation, which can then be solved with sparsity constraints. For the estimation of elastic parameters, a high-resolution and reliable seismic basis pursuit inversion is developed based on the exact Zoeppritz equation. Furthermore, the [Formula: see text] norm is developed as a constraint, wherein a regularized function is minimized with the alternating direction method of the multipliers algorithm. Numerical examples and real data applications demonstrate that our method can not only improve the accuracy of the inversion results, especially the S-wave velocity and density information but also increase the resolution of the inversion results. Furthermore, the [Formula: see text]-norm constraint has better noise suppression demonstrating great potential in practical applications.