Abstract:We have developed a new sparse-spike deconvolution (SSD) method based on Toeplitz-sparse matrix factorization (TSMF), a bilinear decomposition of a matrix into the product of a Toeplitz matrix and a sparse matrix, to address the problems of lateral continuity, effects of noise, and wavelet estimation error in SSD. Assuming the convolution model, a constant source wavelet, and the sparse reflectivity, a seismic profile can be considered as a matrix that is the product of a Toeplitz wavelet matrix and a sparse reflectivity matrix. Thus, we have developed an algorithm of TSMF to simultaneously deconvolve the seismic matrix into a wavelet matrix and a reflectivity matrix by alternatively solving two inversion subproblems related to the Toeplitz wavelet matrix and sparse reflectivity matrix, respectively. Because the seismic wavelet is usually compact and smooth, the fused Lasso was used to constrain the elements in the Toeplitz wavelet matrix. Moreover, due to the limitations of computer memory, large seismic data sets were divided into blocks, and the average of the source wavelets deconvolved from these blocks via TSMF-based SSD was used as the final estimation of the source wavelet for all blocks to deconvolve the reflectivity; thus, the lateral continuity of the seismic data can be maintained. The advantages of the proposed deconvolution method include using multiple traces to reduce the effect of random noise, tolerance to errors in the initial wavelet estimation, and the ability to preserve the complex structure of the seismic data without using any lateral constraints. Our tests on the synthetic seismic data from the Marmousi2 model and a section of field seismic data demonstrate that the proposed method can effectively derive the wavelet and reflectivity simultaneously from band-limited data with appropriate lateral coherence, even when the seismic data are contaminated by noise and the initial wavelet estimation is inaccurate.

A Sparse Spike Deconvolution Algorithm Based on a Recurrent Neural Network and the Iterative Shrinkage-Thresholding Algorithm

Sparse Deconvolution Method for Ultrasound Images Based on Automatic Estimation of Reference Signals.

High-precision Seismic Data Reconstruction with Multi-Domain Sparsity Constraints Based on Curvelet and High-Resolution Radon Transforms

A Nonstationary Sparse Spike Deconvolution with Anelastic Attenuation

Seismic Data Two-Step Recovery Approach Combining Sparsity-Promoting and Hyperbolic Radon Transform Methods

Adaptive step-size fast iterative shrinkage-thresholding algorithm and sparse-spike deconvolution

Seismic sparse-spike deconvolution via Toeplitz-sparse matrix factorization

Deep learning for multitrace sparse-spike deconvolution

Sparse Spike Deconvolution of Seismic Data Using Trust-Region Based Sqp Algorithm

Seismic Sparse-Spike Deconvolution Using Toeplitz-sparse Matrix Factorization

Blind Sparse-Spike Deconvolution with Thin Layers and Structure

Sparse-spike seismic inversion with semismooth newton algorithm solver

Borehole-driven sparse-spike deconvolution and seismic bandwidth enhancement

Sparse Seismic Data Reconstruction Based on a Convolutional Neural Network Algorithm

A Space-Time Adaptive Processing Method Based on Sparse Reconstruction of Reverberation Interference

Sparse Time–Frequency Analysis of Seismic Data: Sparse Representation to Unrolled Optimization

The Effectiveness and Application of the Gradient Sparse Regularization-Based Deconvolution Method

Frequency recovery of band-limited seismic data based on sparse spike train deconvolution

Sparse Seismic Deconvolution with Erratic Data

Application of sparse S transform network with knowledge distillation in seismic attenuation delineation

Deep Sparse Coding Using Optimized Linear Expansion of Thresholds