Compressed wavefield extrapolation with curvelets

An explicit algorithm for the extrapolation of one-way wavefields is proposed which combines recent developments in information theory and theoretical signal processing with the physics of wave propagation. Because of excessive memory requirements, explicit formulations for wave propagation have proven to be a challenge in 3-D. By using ideas from “compressed sensing”, we are able to formulate the (inverse) wavefield extrapolation problem on small subsets of the data volume, thereby reducing the size of the operators. According to compressed sensing theory, signals can successfully be recovered from an imcomplete set of measurements when the measurement basis is incoherent with the representation in which the wavefield is sparse. In this new approach, the eigenfunctions of the Helmholtz operator are recognized as a basis that is incoherent with curvelets that are known to compress seismic wavefields. By casting the wavefield extrapolation problem in this framework, wavefields can successfully be extrapolated in the modal domain via a computationally cheaper operatoion. A proof of principle for the “compressed sensing” method is given for wavefield extrapolation in 2-D. The results show that our method is stable and produces identical results compared to the direct application of the full extrapolation operator.

[1]  F. Herrmann,et al.  Non-linear primary-multiple separation with directional curvelet frames , 2007 .

[2]  Huub Douma,et al.  Leading-order seismic imaging using curvelets , 2007 .

[3]  A. Gisolf,et al.  Fourier reconstruction of marine-streamer data in four spatial coordinates , 2006 .

[4]  Laurent Demanet,et al.  Fast Discrete Curvelet Transforms , 2006, Multiscale Model. Simul..

[5]  Robert G. Fuller,et al.  Numerical Computations in US Undergraduate Physics Courses , 2006, Computing in Science & Engineering.

[6]  Cornelis Vuik,et al.  A new iterative solver for the time-harmonic wave equation , 2006 .

[7]  Yaakov Tsaig,et al.  Extensions of compressed sensing , 2006, Signal Process..

[8]  Hervé Chauris,et al.  Seismic imaging in the curvelet domain and its implications for the curvelet design , 2006 .

[9]  David L Donoho,et al.  Compressed sensing , 2006, IEEE Transactions on Information Theory.

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

[11]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[12]  René-Édouard Plessix,et al.  How to choose a subset of frequencies in frequency-domain finite-difference migration , 2004 .

[13]  E. Candès,et al.  The curvelet representation of wave propagators is optimally sparse , 2004, math/0407210.

[14]  R. Plessix,et al.  Frequency-domain finite-difference amplitude-preserving migration , 2004 .

[15]  Mi Lu,et al.  Accelerating seismic migration using FPGA-based coprocessor platform , 2004, 12th Annual IEEE Symposium on Field-Programmable Custom Computing Machines.

[16]  A. Guitton,et al.  Adaptive subtraction of multiples using the L1‐norm , 2004 .

[17]  Robert D. Nowak,et al.  An EM algorithm for wavelet-based image restoration , 2003, IEEE Trans. Image Process..

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

[19]  E. Candès,et al.  Recovering edges in ill-posed inverse problems: optimality of curvelet frames , 2002 .

[20]  Ru-Shan Wu,et al.  Generalization of the phase-screen approximation for the scattering of acoustic waves , 2000 .

[21]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

[22]  Kees Wapenaar,et al.  Modal expansion of one‐way operators in laterally varying media , 1998 .

[23]  Mauricio D. Sacchi,et al.  Interpolation and extrapolation using a high-resolution discrete Fourier transform , 1998, IEEE Trans. Signal Process..

[24]  F. J. Dessing,et al.  A wavelet transform approach to seismic processing , 1997 .

[25]  Patrick Joly,et al.  Splitting of Operators, Alternate Directions, and Paraxial Approximations for the Three-Dimensional Wave Equation , 1995, SIAM J. Sci. Comput..

[26]  Mauricio D. Sacchi,et al.  Minimum entropy deconvolution with frequency-domain constraints , 1994 .

[27]  D. Hale,et al.  Imaging salt with turning seismic waves , 1992 .

[28]  J. Wang,et al.  Subsurface imaging using magnetotelluric data , 1988 .

[29]  T. Ulrych,et al.  Analytic minimum entropy deconvolution , 1982 .

[30]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[31]  A. J. Berkhout,et al.  Seismic Migration: Imaging of Acoustic Energy by Wave Field Extrapolation , 1980 .

[32]  J. Claerbout.  Toward a unified theory of reflector mapping , 1971 .