Compact Smoothness and Relative Sparsity Algorithm for High-Resolution Wavelet and Reflectivity Inversion of Seismic Data

Wavelet and reflectivity inversion (WRI) is an important issue in seismic data processing. To overcome the ill-posedness of WRI inversion with more efficient parameter selection and better lateral continuity of reflectivities, we propose a new WRI algorithm named compact smoothness and relative sparsity (CSRS) algorithm, where a normalized compact constraint and a normalized smooth regularization are proposed for the wavelet inversion, and a relative sparsity constraint is proposed for the reflectivity inversion. The proposed constraints and regularization make the parameters of WRI easy to be selected. The proposed relative sparsity constraint can lead to a reflectivity profile with good lateral continuity, as it can be suitable for various seismic data with a fixed sparsity parameter. We also propose an efficient algorithm for solving corresponding WRI optimization problem. The whole WRI problem is divided into reflectivity inversion subproblem and wavelet inversion subproblem by using alternating iterative method, where the initial wavelet is estimated by smoothing the absolute amplitude spectrum of averaged seismic data. The proximal algorithm is applied to solve both reflectivity inversion subproblem and wavelet inversion subproblem. By replacing Toeplitz matrix multiplication with the fast Fourier transform (FFT) and using compact wavelet, our algorithm can be efficient for 3-D seismic data. The numerical examples on 2-D synthetic data, 2-D offshore field data, and 3-D onshore field data demonstrate that, compared to the Toeplitz-sparse matrix factorization (TSMF) algorithm, the CSRS algorithm with fixed default parameters can get high-resolution reflectivities with better lateral continuity, and requires much less computational time.

Automated summary

This paper proposes a new algorithm called compact smoothness and relative sparsity (CSRS) algorithm for wavelet and reflectivity inversion (WRI) in seismic data processing. The algorithm introduces normalized constraints and regularization to overcome the ill-posedness of WRI inversion and improve the lateral continuity of reflectivities. It also presents an efficient algorithm for solving the optimization problem of WRI by dividing it into reflectivity inversion and wavelet inversion subproblems.