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. \copyright2005 Society of Exploration Geophysicists