Three-dimensional irregular seismic data reconstruction via low-rank matrix completion

ABSTRACTWe have developed a new algorithm for the reconstruction of seismic traces randomly missing from a uniform grid of a 3D seismic volume. Several algorithms have been developed for such reconstructions, based on properties of the seismic wavefields and on signal processing concepts, such as sparse signal representation in a transform domain. We have investigated a novel approach, originally introduced for noise removal, which is based on the premise that for suitable representation of the seismic data as matrices or tensors, the rank of the seismic data (computed by singular value decomposition) increases with noise or missing traces. Thus, we apply low-rank matrix completion (MC) with a designed texture-patch transformation to 3D seismic data reconstruction. Low-rank components capture geometrically meaningful structures in seismic data that encompass conventional local features such as events and dips. The low-rank MC is based on nuclear-norm minimization. An efficient L1-norm minimizing algorithm...

[1]  Andrea Montanari,et al.  Message passing algorithms for compressed sensing: II. analysis and validation , 2009, 2010 IEEE Information Theory Workshop on Information Theory (ITW 2010, Cairo).

[2]  Felix J. Herrmann,et al.  Application of randomized sampling schemes to curvelet-based sparsity-promoting seismic data recovery , 2013 .

[3]  Mauricio D. Sacchi,et al.  f-x adaptive seismic-trace interpolation , 2009 .

[4]  Kristopher A. Innanen,et al.  Seismic data interpolation using a fast generalized Fourier transform , 2011 .

[5]  Mauricio D. Sacchi,et al.  FX Singular Spectrum Analysis , 2009 .

[6]  Necati Gulunay,et al.  Seismic trace interpolation in the Fourier transform domain , 2003 .

[7]  Mauricio D. Sacchi,et al.  Convergence improvement and noise attenuation considerations for beyond alias projection onto convex sets reconstruction , 2013 .

[8]  Milton J. Porsani,et al.  Seismic trace interpolation using half-step prediction filters , 1999 .

[9]  Ruslan Salakhutdinov,et al.  Practical Large-Scale Optimization for Max-norm Regularization , 2010, NIPS.

[10]  Andrea Montanari,et al.  Message-passing algorithms for compressed sensing , 2009, Proceedings of the National Academy of Sciences.

[11]  Mauricio D. Sacchi,et al.  Beyond alias hierarchical scale curvelet interpolation of regularly and irregularly sampled seismic data , 2010 .

[12]  Lynn Burroughs,et al.  Rank-Reduction-Based Trace Interpolation , 2010 .

[13]  Felix J. Herrmann,et al.  Higher Dimensional Blue-noise Sampling Schemes For Curvelet-based Seismic Data Recovery , 2009 .

[14]  Rabab Kreidieh Ward,et al.  Some empirical advances in matrix completion , 2011, Signal Process..

[15]  Yin Zhang,et al.  Solving a low-rank factorization model for matrix completion by a nonlinear successive over-relaxation algorithm , 2012, Mathematical Programming Computation.

[16]  Side Jin,et al.  5D seismic data regularization by a damped least-norm Fourier inversion , 2010 .

[17]  Stanley Osher,et al.  A Low Patch-Rank Interpretation of Texture , 2013, SIAM J. Imaging Sci..

[18]  Tamara G. Kolda,et al.  MATLAB tensor classes for fast algorithm prototyping. , 2004 .

[19]  Mauricio D. Sacchi,et al.  POCS Method for Seismic Data Reconstruction of Irregularly Sampled Data , 2013 .

[20]  M. Sacchi,et al.  Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis , 2011 .

[21]  Jon F. Claerbout,et al.  Interpolation with smoothly nonstationary prediction-error filters , 1999 .

[22]  Yu Zhang,et al.  Antileakage Fourier transform for seismic data regularization , 2005 .

[23]  Ryota Tomioka,et al.  Estimation of low-rank tensors via convex optimization , 2010, 1010.0789.

[24]  N. Kreimer,et al.  Evaluation of a New 5D Seismic Volume Reconstruction Method: Tensor Completion versus Fourier Reconstruction , 2011 .

[25]  Five Dimensional Interpolation , 2008 .

[26]  S. Osher,et al.  Seismic data reconstruction via matrix completion , 2013 .

[27]  Mauricio D. Sacchi,et al.  Minimum weighted norm interpolation of seismic records , 2004 .

[28]  S. Spitz.  Seismic trace interpolation in the F-X domain , 1991 .

[29]  Pablo A. Parrilo,et al.  Guaranteed Minimum-Rank Solutions of Linear Matrix Equations via Nuclear Norm Minimization , 2007, SIAM Rev..

[30]  Tad J. Ulrych,et al.  Application of singular value decomposition to vertical seismic profiling , 1988 .

[31]  Mauricio D. Sacchi,et al.  A fast reduced-rank interpolation method for prestack seismic volumes that depend on four spatial dimensions , 2013 .

[32]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

[33]  D. J. Verschuur,et al.  Restoration of missing offsets by parabolic Radon transform1 , 1995 .

[34]  Jieping Ye,et al.  Tensor Completion for Estimating Missing Values in Visual Data , 2013, IEEE Trans. Pattern Anal. Mach. Intell..

[35]  J. Suykens,et al.  Nuclear Norms for Tensors and Their Use for Convex Multilinear Estimation , 2011 .

[36]  C. Brezinski.  Interpolation and Extrapolation , 2001 .

[37]  Tamara G. Kolda,et al.  Categories and Subject Descriptors: G.4 [Mathematics of Computing]: Mathematical Software— , 2022 .

[38]  Mauricio D. Sacchi,et al.  Fourier Reconstruction of Nonuniformly Sampled, Aliased Seismic Data , 2022 .

[39]  Emmanuel J. Candès,et al.  Exact Matrix Completion via Convex Optimization , 2009, Found. Comput. Math..

[40]  S. Trickett.  F-xy eigenimage noise suppression , 2003 .

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

[42]  Daniel Trad,et al.  Five-dimensional interpolation: Recovering from acquisition constraints , 2009 .

[43]  Felix J. Herrmann,et al.  Non-parametric seismic data recovery with curvelet frames , 2008 .

[44]  Stewart Trickett,et al.  F-xy Cadzow Noise Suppression , 2008 .

[45]  B. Recht,et al.  Tensor completion and low-n-rank tensor recovery via convex optimization , 2011 .

[46]  Gerlind Plonka-Hoch,et al.  Compressive Video Sampling With Approximate Message Passing Decoding , 2012, IEEE Transactions on Circuits and Systems for Video Technology.

[47]  Yu Zhang,et al.  Antileakage Fourier transform for seismic data regularization in higher dimensions , 2010 .

[48]  Nadia Kreimer,et al.  A tensor higher-order singular value decomposition for prestack seismic data noise reduction and interpolation , 2012 .

[49]  R. Abma,et al.  3D interpolation of irregular data with a POCS algorithm , 2006 .

[50]  Joshua Ronen,et al.  Wave‐equation trace interpolation , 1987 .