Dictionary Learning and Sparse Representations for Denoising and Reconstruction of Marine Seismic Data

We have addressed the seismic data denoising problem, in which the noise is random and has an unknown spatiotemporally varying variance. In seismic data processing, random noise is often attenuated using transform-based methods. The success of these methods in denoising depends on the ability of the transform to efficiently describe the signal features in the data. Fixed transforms (e.g., wavelets, curvelets) do not adapt to the data and might fail to efficiently describe complex morphologies in the seismic data. Alternatively, dictionary learning methods adapt to the local morphology of the data and provide state-of-the-art denoising results. However, conventional denoising by dictionary learning requires a priori information on the noise variance, and it encounters difficulties when applied for denoising seismic data in which the noise variance is varying in space or time. Here, we propose a coherence-constrained dictionary learning (CDL) method for denoising that does not require any a priori information related to the signal or noise. To denoise a given window of a seismic section using CDL, overlapping small 2D patches are extracted and a dictionary of patch-size signals is trained to learn the elementary features embedded in the seismic signal. For each patch, using the learned dictionary, a sparse optimization problem is solved, and a sparse approximation of the patch is computed to attenuate the random noise. Unlike conventional dictionary learning, the sparsity of the approximation is constrained based on coherence such that it does not need a priori noise variance or signal sparsity information and is still optimal to filter out Gaussian random noise. The denoising performance of the CDL method is validated using synthetic and field data examples, and is compared with the KSVD and FX-Decon denoising. We found that CDL gives better denoising results than K-SVD and FX-Decon for removing noise when the variance varies in space or time.

[1]  Rajiv Kumar,et al.  Highly repeatable 3D compressive full-azimuth towed-streamer time-lapse acquisition — A numerical feasibility study at scale , 2017 .

[2]  Joel A. Tropp,et al.  Greed is good: algorithmic results for sparse approximation , 2004, IEEE Transactions on Information Theory.

[3]  S. Sahoo,et al.  Dictionary Training for Sparse Representation as Generalization of K-Means Clustering , 2013, IEEE Signal Processing Letters.

[4]  Michael Elad,et al.  Double Sparsity: Learning Sparse Dictionaries for Sparse Signal Approximation , 2010, IEEE Transactions on Signal Processing.

[5]  Kjersti Engan,et al.  Method of optimal directions for frame design , 1999, 1999 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings. ICASSP99 (Cat. No.99CH36258).

[6]  J. MacQueen Some methods for classification and analysis of multivariate observations , 1967 .

[7]  S. Mallat,et al.  Adaptive greedy approximations , 1997 .

[8]  Walter Söllner,et al.  A method of combining coherence-constrained sparse coding and dictionary learning for denoising , 2017 .

[9]  T. Elboth,et al.  Attenuation of Noise In Marine Seismic Data , 2009 .

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

[11]  Michael Elad,et al.  Optimally sparse representation in general (nonorthogonal) dictionaries via ℓ1 minimization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

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

[13]  Michael Elad,et al.  Multi-Scale Dictionary Learning Using Wavelets , 2011, IEEE Journal of Selected Topics in Signal Processing.

[14]  Y. C. Pati,et al.  Orthogonal matching pursuit: recursive function approximation with applications to wavelet decomposition , 1993, Proceedings of 27th Asilomar Conference on Signals, Systems and Computers.

[15]  Jianwei Ma,et al.  Seismic data restoration via data-driven tight frame , 2014 .

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

[17]  Dengliang Gao,et al.  Volume texture extraction for 3D seismic visualization and interpretation , 2003 .

[18]  Thomas Elboth,et al.  Noise in Marine Seismic Data , 2010 .

[19]  Charles C. Mosher,et al.  Wavelet Transform Methods For Geophysical Applications , 1994 .

[20]  Robert M. Haralick,et al.  Textural Features for Image Classification , 1973, IEEE Trans. Syst. Man Cybern..

[21]  Wei Wang,et al.  Ground-roll Separation By Sparsity And Morphological Diversity Promotion , 2010 .

[22]  Charles C. Mosher,et al.  Wavelet transform methods for geophysical applications , 1993, Optics & Photonics.

[23]  Bennett Eisenberg,et al.  Why Is the Sum of Independent Normal Random Variables Normal? , 2008 .

[24]  Walter Söllner,et al.  Dictionary learning for signal-to-noise ratio enhancement , 2015 .

[25]  Thomas Elboth,et al.  Flow and swell noise in marine seismic data , 2009 .

[26]  Klaus Mosegaard,et al.  Seismic Texture Classification: A Computer-aided Approach to Stratigraphic Analysis , 1995 .

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

[28]  Yangkang Chen,et al.  Double Sparsity Dictionary for Seismic Noise Attenuation , 2016 .

[29]  Deanna Needell,et al.  CoSaMP: Iterative signal recovery from incomplete and inaccurate samples , 2008, ArXiv.

[30]  R. Neelamani,et al.  Coherent and random noise attenuation using the curvelet transform , 2008 .

[31]  S. M. Doherty,et al.  Seismic Data Analysis: Processing, Inversion, and Interpretation of Seismic Data , 2000 .

[32]  Don R. Hush,et al.  Network constraints and multi-objective optimization for one-class classification , 1996, Neural Networks.

[33]  J. Kovacevic,et al.  Life Beyond Bases: The Advent of Frames (Part I) , 2007, IEEE Signal Processing Magazine.

[35]  A. Bruckstein,et al.  K-SVD : An Algorithm for Designing of Overcomplete Dictionaries for Sparse Representation , 2005 .

[36]  Bhaskar D. Rao,et al.  Sparse signal reconstruction from limited data using FOCUSS: a re-weighted minimum norm algorithm , 1997, IEEE Trans. Signal Process..

[37]  Cássio Fraga Dantas,et al.  Learning Dictionaries as a Sum of Kronecker Products , 2017, IEEE Signal Processing Letters.

[38]  Margaret Yu Seismic interference noise elimination - a multidomain 3D filtering approach , 2011 .

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

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

[41]  Sergey Fomel,et al.  OC-seislet: Seislet transform construction with differential offset continuation , 2010 .

[42]  Jean-Luc Starck,et al.  Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit , 2012, IEEE Transactions on Information Theory.

[43]  James H. McClellan,et al.  Seismic data denoising through multiscale and sparsity-promoting dictionary learning , 2015 .

[44]  Vicente Oropeza,et al.  The Singular Spectrum Analysis method and its application to seismic data denoising and reconstruction , 2010 .

[45]  Necati Gulunay,et al.  FXDECON and complex wiener prediction filter , 1986 .

[46]  Michael Elad,et al.  Analysis versus synthesis in signal priors , 2006, 2006 14th European Signal Processing Conference.

[47]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[48]  Stanley Osher,et al.  Monte Carlo data-driven tight frame for seismic data recovery , 2016 .

[49]  Michael Elad,et al.  Efficient Implementation of the K-SVD Algorithm using Batch Orthogonal Matching Pursuit , 2008 .

[50]  Mauricio D. Sacchi,et al.  On sampling functions and Fourier reconstruction methods , 2010 .

[51]  Stéphane Mallat,et al.  A Wavelet Tour of Signal Processing - The Sparse Way, 3rd Edition , 2008 .

[52]  B. P. West,et al.  Interactive seismic facies classification using textural attributes and neural networks , 2002 .

[53]  Tadeusz J. Ulrych,et al.  Information-Based Inversion and Processing with Applications , 2011 .

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

[55]  Michael Elad,et al.  Analysis and synthesis sparse modeling methods image processing , 2011 .

[56]  A. Day,et al.  Wavefield-separation methods for dual-sensor towed-streamer data , 2013 .

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

[58]  S. Berman Limit Theorems for the Maximum Term in Stationary Sequences , 1964 .

[59]  Mauricio D. Sacchi,et al.  Interpolation and denoising of high-dimensional seismic data by learning a tight frame , 2015 .

[60]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[61]  Nasser Kazemi,et al.  Attenuation of swell noise in marine streamer data via nonnegative matrix factorization , 2016 .

[62]  P. Hubral Computing true amplitude reflections in a laterally inhomogeneous earth , 1983 .

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

[64]  Walter Söllner,et al.  Sparsity promoting morphological decomposition for coherent noise suppression: Application to streamer vibration related noise , 2016 .

[65]  Michael Elad,et al.  Why Simple Shrinkage Is Still Relevant for Redundant Representations? , 2006, IEEE Transactions on Information Theory.

[66]  Felix J. Herrmann,et al.  Randomized sampling and sparsity: Getting more information from fewer samples , 2010 .

[67]  A. Gisolf,et al.  Fourier reconstruction with sparse inversion , 2007 .

[68]  Philippe Herrmann,et al.  De-aliased, High-Resolution Radon Transforms , 2000 .

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

[70]  Michael Elad,et al.  Sparse and Redundant Representations - From Theory to Applications in Signal and Image Processing , 2010 .