Wavefield reconstruction is a crucial step in the seismic processing flow. For instance, unsuccessful interpolation leads to erroneous multiple predictions that adversely affect the performance of multiple elimination, and to imaging artifacts. We present a new non-parametric transform-based reconstruction method that exploits the compression of seismic data b the recently developed curvelet transform. The elements of this transform, called curvelets, are multi-dimensional, multi-scale, and multi-directional. They locally resemble wavefronts present in the data, which leads to a compressible representation for seismic data. This compression enables us to formulate a new curvelet-based seismic data recovery algorithm through sparsity-promoting inversion (CRSI). The concept of sparsity-promoting inversion is in itself not new to geophysics. However, the recent insights from the field of "compressed sensing" are new since they clearly identify the three main ingredients that go into a successful formulation of a reconstruction problem, namely a sparsifying transform, a sub-Nyquist sampling strategy that subdues coherent aliases in the sparsifying domain, and a data-consistent sparsity-promoting program. After a brief overview of the curvelet transform and our seismic-oriented extension to the fast discrete curvelet transform, we detail the CRSI formulation and illustrate its performance on synthetic and read datasets. Then, we introduce a sub-Nyquist sampling scheme, termed jittered undersampling, and show that, for the same amount of data acquired, jittered data are best interpolated using CRSI compared to regular or random undersampled data. We also discuss the large-scale one-norm solver involved in CRSI. Finally, we extend CRSI formulation to other geophysical applications and present results on multiple removal and migration-amplitude recovery.

Wavefield reconstruction is a crucial step in the seismic processing flow. For instance, unsuccessful interpolation leads to erroneous multiple predictions that adversely affect the performance of multiple elimination, and to imaging artifacts. We present a new non-parametric transform-based reconstruction method that exploits the compression of seismic data b the recently developed curvelet transform. The elements of this transform, called curvelets, are multi-dimensional, multi-scale, and multi-directional. They locally resemble wavefronts present in the data, which leads to a compressible representation for seismic data. This compression enables us to formulate a new curvelet-based seismic data recovery algorithm through sparsity-promoting inversion (CRSI). The concept of sparsity-promoting inversion is in itself not new to geophysics. However, the recent insights from the field of "compressed sensing" are new since they clearly identify the three main ingredients that go into a successful formulation of a reconstruction problem, namely a sparsifying transform, a sub-Nyquist sampling strategy that subdues coherent aliases in the sparsifying domain, and a data-consistent sparsity-promoting program. After a brief overview of the curvelet transform and our seismic-oriented extension to the fast discrete curvelet transform, we detail the CRSI formulation and illustrate its performance on synthetic and read datasets. Then, we introduce a sub-Nyquist sampling scheme, termed jittered undersampling, and show that, for the same amount of data acquired, jittered data are best interpolated using CRSI compared to regular or random undersampled data. We also discuss the large-scale one-norm solver involved in CRSI. Finally, we extend CRSI formulation to other geophysical applications and present results on multiple removal and migration-amplitude recovery.

Wavefield reconstruction is a crucial step in the seismic processing flow. For instance, unsuccessful interpolation leads to erroneous multiple predictions that adversely affect the performance of multiple elimination, and to imaging artifacts. We present a new non-parametric transform-based reconstruction method that exploits the compression of seismic data b the recently developed curvelet transform. The elements of this transform, called curvelets, are multi-dimensional, multi-scale, and multi-directional. They locally resemble wavefronts present in the data, which leads to a compressible representation for seismic data. This compression enables us to formulate a new curvelet-based seismic data recovery algorithm through sparsity-promoting inversion (CRSI). The concept of sparsity-promoting inversion is in itself not new to geophysics. However, the recent insights from the field of "compressed sensing" are new since they clearly identify the three main ingredients that go into a successful formulation of a reconstruction problem, namely a sparsifying transform, a sub-Nyquist sampling strategy that subdues coherent aliases in the sparsifying domain, and a data-consistent sparsity-promoting program. After a brief overview of the curvelet transform and our seismic-oriented extension to the fast discrete curvelet transform, we detail the CRSI formulation and illustrate its performance on synthetic and read datasets. Then, we introduce a sub-Nyquist sampling scheme, termed jittered undersampling, and show that, for the same amount of data acquired, jittered data are best interpolated using CRSI compared to regular or random undersampled data. We also discuss the large-scale one-norm solver involved in CRSI. Finally, we extend CRSI formulation to other geophysical applications and present results on multiple removal and migration-amplitude recovery. / Science, Faculty of / Earth, Ocean and Atmospheric Sciences, Department of / Graduate

No description available.

In this paper, we present a nonlinear curvelet-based sparsity-promoting formulation for three problems related to seismic noise, namely the ’good’, corresponding to noise generated by random sampling; the ’bad’, corresponding to coherent noise for which (inaccurate) predictions exist and the ’ugly’ for which no predictions exist. We will show that the compressive capabilities of curvelets on seismic data and images can be used to tackle these three categories of noise-related problems.

In this overview of the recent Curvelet Reconstruction with Sparsity-promoting Inversion (CRSI) method, we present our latest 2-D and 3-D interpolation results on both synthetic and real datasets. We compare these results to interpolated data using other existing methods. Finally, we discuss the challenges related to sparsity-promoting solvers for the large-scale problems the industry faces.

Incomplete data represents a major challenge for a successful prediction and subsequent removal of multiples. In this paper, a new method will be represented that tackles this challenge in a two-step approach. During the first step, the recenly developed curvelet-based recovery by sparsity-promoting inversion (CRSI) is applied to the data, followed by a prediction of the primaries. During the second high-resolution step, the estimated primaries are used to improve the frequency content of the recovered data by combining the focal transform, defined in terms of the estimated primaries, with the curvelet transform. This focused curvelet transform leads to an improved recovery, which can subsequently be used as input for a second stage of multiple prediction and primary-multiple separation.

n this paper, two different applications of phase transitions to exploration seismology will be discussed. The first application concerns a phase diagram ruling the recovery conditions for seismic data volumes from incomplete and noisy data while the second phase transition describes the behavior of bi-compositional mixtures as a function of the volume fraction. In both cases, the phase transitions are the result of randomness in large system of equations in combination with nonlinearity. The seismic recovery problem from incomplete data involves the inversion of a rectangular matrix. Recent results from the field of "compressive sensing" provide the conditions for a successful recovery of functions that are sparse in some basis (wavelet) or frame (curvelet) representation, by means of a sparsity ($\ell_1$-norm) promoting nonlinear program. The conditions for a successful recovery depend on a certain randomness of the matrix and on two parameters that express the matrix' aspect ratio and the ratio of the number of nonzero entries in the coefficient vector for the sparse signal representation over the number of measurements. It appears that the ensemble average for the success rate for the recovery of the sparse transformed data vector by a nonlinear sparsity promoting program, can be described by a phase transition, demarcating the regions for the two ratios for which recovery of the sparse entries is likely to be successful or likely to fail. Consistent with other phase transition phenomena, the larger the system the sharper the transition. The randomness in this example is related to the construction of the matrix, which for the recovery of spike trains corresponds to the randomly restricted Fourier matrix. It is shown, that these ideas can be extended to the curvelet recovery by sparsity-promoting inversion (CRSI) . The second application of phase transitions in exploration seismology concerns the upscaling problem. To counter the intrinsic smoothing of singularities by conventional equivalent medium upscaling theory, a percolation-based nonlinear switch model is proposed. In this model, the transport properties of bi-compositional mixture models for rocks undergo a sudden change in the macroscopic transport properties as soon as the volume fraction of the stronger material reaches a critical point. At this critical point, the stronger material forms a connected cluster, which leads to the creation of a cusp-like singularity in the elastic moduli, which in turn give rise to specular reflections. In this model, the reflectivity is no longer explicitly due to singularities in the rocks composition. Instead, singularities are created whenever the volume fraction exceeds the critical point. We will show that this concept can be used for a singularity-preserved lithological upscaling.

In this paper, we present a nonlinear curvelet-based sparsity-promoting formulation for three problems in seismic processing and imaging namely, seismic data regularization from data with large percentages of traces missing; seismic amplitude recovery for subsalt images obtained by reverse-time migration and primary-multiple separation, given an inaccurate multiple prediction. We argue why these nonlinear formulations are beneficial.

Incomplete data represents a major challenge for a successful prediction and subsequent removal of multiples. In this paper, a new method will be represented that tackles this challenge in a two-step approach. During the first step, the recenly developed curvelet-based recovery by sparsity-promoting inversion (CRSI) is applied to the data, followed by a prediction of the primaries. During the second high-resolution step, the estimated primaries are used to improve the frequency content of the recovered data by combining the focal transform, defined in terms of the estimated primaries, with the curvelet transform. This focused curvelet transform leads to an improved recovery, which can subsequently be used as input for a second stage of multiple prediction and primary-multiple separation.