5D and 4D pre-stack seismic data completion using tensor nuclear norm (TNN)

In this paper we present novel strategies for completion of 5D pre-stack seismic data, viewed as a 5D tensor or as a set of 4D tensors across temporal frequencies. In contrast to existing complexity penalized algorithms for seismic data completion, which employ matrix analogues of tensor decompositions such as HOSVD or use overlapped Schatten norms from different unfoldings (or matricization) of the tensors, our approach uses a recently proposed decomposition called tensor SVD or tSVD for short, proposed in [Kilmer and Martin (2011)]. We show that seismic data exhibits low complexity under tSVD, i.e. is compressible under tSVD representation, and we subsequently propose a new complexity penalized algorithm for pre-stack seismic data completion under missing traces. This complexity measure which we call the Tensor Nuclear Norm (TNN) is motivated by algebraic properties of the tSVD. We test the performance of the proposed algorithms on synthetic and real data and show that missing data can be reliably recovered under heavy down-sampling.

[1]  Hassan Mansour,et al.  Seismic Data Interpolation and Denoising Using SVD-free Low-rank Matrix Factorization , 2013, London 2013, 75th eage conference en exhibition incorporating SPE Europec.

[2]  Misha Elena Kilmer,et al.  Third-Order Tensors as Operators on Matrices: A Theoretical and Computational Framework with Applications in Imaging , 2013, SIAM J. Matrix Anal. Appl..

[3]  Carla D. Martin,et al.  An Order-p Tensor Factorization with Applications in Imaging , 2013, SIAM J. Sci. Comput..

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

[5]  Emmanuel J. Candès,et al.  Near-Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? , 2004, IEEE Transactions on Information Theory.

[6]  N. Kreimer,et al.  Tensor completion based on nuclear norm minimization for 5D seismic data reconstruction , 2013 .

[7]  Nadia Kreimer,et al.  Reconstruction of seismic data via tensor completion , 2012, 2012 IEEE Statistical Signal Processing Workshop (SSP).

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

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

[10]  Nadia Kreimer,et al.  Tensor Completion via Nuclear Norm Minimization for 5D Seismic Data Reconstruction , 2012 .

[11]  Venkatesh Saligrama,et al.  Information Theoretic Bounds for Compressed Sensing , 2008, IEEE Transactions on Information Theory.

[12]  Nadia Kreimer,et al.  A Tensor Higher-order Singular Value Decomposition (HOSVD) For Pre-stack Simultaneous Noise-reduction And Interpolation , 2011 .

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

[14]  Mauricio D. Sacchi,et al.  A Fast Rank Reduction Method For the Reconstruction of 5D Seismic Volumes , 2011 .

[15]  M. Kilmer,et al.  Factorization strategies for third-order tensors , 2011 .

[16]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[17]  D. Trad.  Five Dimensional Interpolation , 2008 .

[18]  G. Watson.  Characterization of the subdifferential of some matrix norms , 1992 .

[19]  Felix J. Herrmann,et al.  Hierarchical Tucker Tensor Optimization - Applications to 4D Seismic Data Interpolation , 2013 .

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

[21]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[22]  Felix J. Herrmann,et al.  Seismic denoising with nonuniformly sampled curvelets , 2006, Computing in Science & Engineering.

[23]  Venkatesh Saligrama,et al.  Thresholded Basis Pursuit: LP Algorithm for Order-Wise Optimal Support Recovery for Sparse and Approximately Sparse Signals From Noisy Random Measurements , 2011, IEEE Transactions on Information Theory.