Sparseness‐constrained least‐squares inversion: Application to seismic wave reconstruction

Abstract

The spectrum of a discrete Fourier transform (DFT) is estimated by linear inversion, and used to produce desirable seismic traces with regular spatial sampling from an irregularly sampled data set. The essence of such a wavefield reconstruction method is to solve the DFT inverse problem with a particular constraint which imposes a sparseness criterion on the least-squares solution. A working definition for the sparseness constraint is presented to improve the stability and efficiency. Then a sparseness measurement is used to measure the relative sparseness of the two DFT spectra obtained from inversion with or without sparseness constraint. It is a pragmatic indicator about the magnitude of sparseness needed for wavefield reconstruction. For seismic trace regularization, an antialiasing condition must be fulfilled for the regularizing trace interval, whereas optimal trace coordinates in the output can be obtained by minimizing the distances between the newly generated traces and the original traces in the input. Application to real seismic data reveals the effectiveness of the technique and the significance of the sparseness constraint in the least-squares solution.

Efficient waveform tomography for lithospheric imaging: implications for realistic, two-dimensional acquisition geometries and low-frequency data

“…Trace interpolation has been developed as a solution to this aliasing problem through a number of different methods, including differential offset continuation (Fomel 2003), and minimum weighted norm interpolation (Liu & Sacchi 2004). Wang (2002) has developed interpolation in the frequency‐space domain for regular sampling, and in the frequency‐wavenumber domain for irregular sampling (Wang 2003).…”

An open-source Matlab code package for improved rank-reduction 3D seismic data denoising and reconstruction

“…However, the predictive filtering method can only be applied to regularly sampled seismic data. The second type is a transformed domain method (Candès et al, 2006a;Chen et al, 2014b;Gan et al, 2015b;Liu et al, 2015), which is based on compressive sensing theory (Candès et al, 2006b;Donoho, 2006) to achieve a successful recovery using highly incomplete available data (Sacchi et al, 1998;Wang, 2003;Chen et al, 2014a). Compressive sensing (CS) is a relatively new paradigm (Candès et al, 2006b;Donoho, 2006;Donoho et al, 2006;Gan et al, 2016;Liu et al, 2016) in signal processing that has recently received a lot of attention.…”

Curvelet-based POCS interpolation of nonuniformly sampled seismic records

“…By enforcing a Cauchy-Gaussian a priori sparseness within the Bayesian framework, Sacchi et al (1998) proposed a high-resolution Fourier transform to perform interpolation and extrapolation. Sparseness-constrained Fourier reconstruction was also successfully applied to the irregularly sampled seismic data, using least-squares criterion and minimum (weighted) norm constraint (Liu and Sacchi, 2004;Wang, 2003;Zwartjes and Gisolf, 2007;Zwartjes and Sacchi, 2007). Sparsity-promoting seismic trace restoration was even investigated in Radon domain (Kabir and Verschuur, 1995;Sacchi and Ulrych, 1995;Trad and Ulrych, 2002;Trad et al, 2003).…”