Sparseness-constrained data continuation with frames: Applications to missing traces and aliased signals in 2/3-D

We present a robust iterative sparseness-constrained interpolation algorithm using 2-/3-D curvelet frames and Fourier-like transforms that exploits continuity along reflectors in seismic data. By choosing generic transforms, we circumvent the necessity to make parametric assumptions (e.g. through linear/parabolic Radon or demigration) regarding the shape of events in seismic data. Simulation and real data examples for data with moderately sized gaps demonstrate that our algorithm provides interpolated traces that accurately reproduce the wavelet shape as well as the AVO behavior. Our method also shows good results for de-aliasing judged by the behavior of the (f − k)-spectrum before and after regularization.

[1]  L. Karlovitz.  Construction of nearest points in the Lp, p even, and L∞ norms. I , 1970 .

[2]  Douglas W. Oldenburg,et al.  Wavelet estimation and deconvolution , 1981 .

[3]  Mauricio D. Sacchi,et al.  Estimation of the discrete Fourier transform, a linear inversion approach , 1996 .

[4]  A. W. F. Volker,et al.  Reconstruction As Efficient Alternative For Least Squares Migration , 2000 .

[5]  A. Duijndam,et al.  Parabolic Radon transform, sampling and efficiency , 2001 .

[6]  Emmanuel J. Candès,et al.  New multiscale transforms, minimum total variation synthesis: applications to edge-preserving image reconstruction , 2002, Signal Process..

[7]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

[8]  Mauricio D. Sacchi,et al.  Latest views of the sparse Radon transform , 2003 .

[9]  Daniel Trad,et al.  Interpolation and multiple attenuation with migration operators , 2003 .

[10]  D. Donoho,et al.  Redundant Multiscale Transforms and Their Application for Morphological Component Separation , 2004 .

[11]  F. Herrmann,et al.  Robust Curvelet-Domain Primary-Multiple Separation with Sparseness Constraints , 2005 .

[12]  Felix J. Herrmann,et al.  Robust Curvelet-Domain Data Continuation with Sparseness Constraints , 2005 .

[13]  D. Donoho,et al.  Simultaneous cartoon and texture image inpainting using morphological component analysis (MCA) , 2005 .

[14]  Felix J. Herrmann,et al.  Application of Stable Signal Recovery to Seismic Data Interpolation , 2006 .