Search

Blind sparse-spike deconvolution is a widely used method to estimate seismic wavelets and sparse reflectivity in the shape of spikes based on the convolution model. To increase the vertical resolution and lateral continuity of the estimated reflectivity, we further improve the sparse-spike deconvolution by introducing the atomic norm minimization and structural regularization, respectively. Specifically, we use the atomic norm minimization to estimate the reflector locations, which are further used as position constraints in the sparse-spike deconvolution. By doing this, we can vertically separate highly thin layers through the sparse deconvolution. In addition, the seismic structural orientations are estimated from the seismic image to construct a structure-guided regularization in the deconvolution to preserve the lateral continuity of reflectivities. Our improvements are suitable for most types of sparse-spike deconvolution approaches. The sparse-spike deconvolution method with Toeplitz-sparse matrix factorization (TSMF) is used as an example to demonstrate the effectiveness of our improvements. Synthetic and real examples show that our methods perform better than TSMF in estimating the reflectivity of thin layers and preserving the lateral continuities.

Seismic wavelet estimation and deconvolution are essential for high-resolution seismic processing. Because of the influence of absorption and scattering, the frequency and phase of the seismic wavelet change with time during wave propagation, leading to a time-varying seismic wavelet. To obtain reflectivity coefficients with more accurate relative amplitudes, we should compute a nonstationary deconvolution of this seismogram, which might be difficult to solve. We have extended sparse spike deconvolution via Toeplitz-sparse matrix factorization to a nonstationary sparse spike deconvolution approach with anelastic attenuation. We do this by separating our model into subproblems in each of which the wavelet estimation problem is solved by the classic sparse optimization algorithms. We find numerical examples that illustrate the parameter setting, noisy seismogram, and the estimation error of the [Formula: see text] value to validate the effectiveness of our extended approach. More importantly, taking advantage of the high accuracy of the estimated [Formula: see text] value, we obtain better performance than with the stationary Toeplitz-sparse spike deconvolution approach in real seismic data.