Robust dictionary learning for erratic noise-corrupted seismic data reconstruction

Dictionary learning methods adaptively train their bases from the given data in an iterative manner; hence, they can capture more detailed features and achieve sparser representation than a method that uses a fixed basis. However, there also exists a good chance of erratic noise corrupting the dictionary because of the insufficiency of the L1-norm regularization in distinguishing the signal and erratic noise, causing the inaccuracy of both dictionary learning and its application. A robust dictionary learning method is proposed to reconstruct the missing data even in the presence of strong erratic noise. Specifically, a projected operator, which is constructed with the robust Huber misfit, is used to damp erratic noise to a low-amplitude level. In the dictionary learning step, erratic noise is gradually reduced to an acceptable level with the help of this projected operator, from which the signal prototypes (or atoms) can be safely extracted in each iteration. Then, the learned dictionary is used to iteratively estimate the signal from the noisy and under-sampled data, performing noise attenuation as well as data interpolation at non-recording locations. We test the proposed dictionary learning method using 2D and 3D noisy seismic examples and compare it with other state-of-the-art methods. Numerical results demonstrate its superior effectiveness at recovering missing data and increasing signal-to-noise ratio in the presence of erratic noise.

[1]  P. Holland,et al.  Robust regression using iteratively reweighted least-squares , 1977 .

[2]  W. Symes,et al.  Robust inversion of seismic data using the Huber norm , 2003 .

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

[4]  Michael Elad,et al.  Image Denoising Via Sparse and Redundant Representations Over Learned Dictionaries , 2006, IEEE Transactions on Image Processing.

[5]  Yangkang Chen,et al.  Fast dictionary learning for noise attenuation of multidimensional seismic data , 2017, Geophysical Journal International.

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

[7]  G. Plonka,et al.  Seismic data interpolation and denoising by learning a tensor tight frame , 2017 .

[8]  Mauricio D. Sacchi,et al.  Minimum weighted norm wavefield reconstruction for AVA imaging , 2005 .

[9]  Jian-Feng Cai,et al.  Data-driven tight frame construction and image denoising , 2014 .

[10]  Mauricio D. Sacchi,et al.  Robust f‐x projection filtering for simultaneous random and erratic seismic noise attenuation , 2017 .

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

[12]  J. Claerbout,et al.  Robust Modeling With Erratic Data , 1973 .

[13]  Felix J. Herrmann,et al.  Fighting the Curse of Dimensionality: Compressive Sensing in Exploration Seismology , 2012, IEEE Signal Processing Magazine.

[14]  Qiang Zhao,et al.  Constrained data-driven tight frame for robust seismic data reconstruction , 2017 .

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

[16]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

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

[18]  Yi Ma,et al.  Robust principal component analysis? , 2009, JACM.

[19]  Anamitra Makur,et al.  Enhancing Image Denoising by Controlling Noise Incursion in Learned Dictionaries , 2015, IEEE Signal Processing Letters.

[20]  Weilin Huang,et al.  Simultaneous denoising and reconstruction of 5-D seismic data via damped rank-reduction method , 2016 .

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

[22]  Kuldeep Kumar,et al.  Robust Statistics, 2nd edn , 2011 .

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

[24]  Jingjie Cao,et al.  Recovery of seismic wavefields based on compressive sensing by an l1-norm constrained trust region method and the piecewise random subsampling , 2011, Geophysical Journal International.

[25]  Guillermo Sapiro,et al.  Online dictionary learning for sparse coding , 2009, ICML '09.

[26]  Qiang Zhao,et al.  Iterative Double Laplacian-Scaled Low-Rank Optimization for Under-Sampled and Noisy Signal Recovery , 2019, IEEE Transactions on Geoscience and Remote Sensing.

[27]  M. Elad,et al.  $rm K$-SVD: An Algorithm for Designing Overcomplete Dictionaries for Sparse Representation , 2006, IEEE Transactions on Signal Processing.

[28]  B. Ripley,et al.  Robust Statistics , 2018, Encyclopedia of Mathematical Geosciences.

[29]  Yangkang Chen,et al.  Data-driven multitask sparse dictionary learning for noise attenuation of 3D seismic data , 2017 .

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

[31]  Thomas C. M. Lee,et al.  Matrix Completion with Noisy Entries and Outliers , 2015, J. Mach. Learn. Res..

[32]  David J. Field,et al.  Emergence of simple-cell receptive field properties by learning a sparse code for natural images , 1996, Nature.

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

[34]  Qiang Zhao,et al.  Signal-Preserving Erratic Noise Attenuation via Iterative Robust Sparsity-Promoting Filter , 2018, IEEE Transactions on Geoscience and Remote Sensing.

[35]  Lynn Burroughs,et al.  Robust Rank-Reduction Filtering for Erratic Noise , 2012 .

[36]  Lina Liu,et al.  Sparse graph-regularized dictionary learning for suppressing random seismic noiseGraph dictionary learning , 2018 .

[37]  J.-C. Pesquet,et al.  A Douglas–Rachford Splitting Approach to Nonsmooth Convex Variational Signal Recovery , 2007, IEEE Journal of Selected Topics in Signal Processing.

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

[39]  Mauricio D. Sacchi,et al.  Robust reduced-rank filtering for erratic seismic noise attenuation , 2015 .

[40]  David Le Meur,et al.  Multidimensional simultaneous random plus erratic noise attenuation and interpolation for seismic data by joint low-rank and sparse inversion , 2015 .

[41]  G. H. F. Gardner,et al.  Effects of irregular sampling on 3‐D prestack migration , 1994 .