Fast online orthonormal dictionary learning for efficient full waveform inversion

Full waveform inversion (FWI) delivers high-resolution images of a subsurface medium property by minimizing itera-tively the misfit between observed and simulated seismic data, and is commonly used by the oil and gas industry for geophysical exploration. FWI is a challenging problem because seismic surveys cover ever larger areas of interest and collect massive volumes of data. The dimensionality of the problem and the heterogeneity of the medium both stress the need for faster algorithms, so sparse regularization techniques can be used to accelerate and improve imaging results. In this paper, we propose a compressive sensing method for the FWI problem by exploiting the sparsity of geological model perturbations over learned dictionaries. Based on stochastic approximations, the dictionaries are updated itera-tively to adapt changing models during FWI iterations. Meanwhile, the dictionaries are kept orthonormal in order to maintain the corresponding transform in a fast and compact manner so that these transforms do not introduce extra computational overhead to FWI. Establishing such a sparsity regu-larization on the model enables us to significantly reduce the workload by only collecting 0.625% of the field data without introducing subsampling artifacts. Hence, the computational burden of large-scale FWI problems can be greatly reduced.

[1]  Sheng Chen,et al.  Orthogonal least squares methods and their application to non-linear system identification , 1989 .

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

[3]  William Grimson,et al.  Object recognition by computer - the role of geometric constraints , 1991 .

[4]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[5]  Léon Bottou,et al.  The Tradeoffs of Large Scale Learning , 2007, NIPS.

[6]  Michael Elad,et al.  Analysis K-SVD: A Dictionary-Learning Algorithm for the Analysis Sparse Model , 2013, IEEE Transactions on Signal Processing.

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

[8]  Mark W. Schmidt,et al.  Optimizing Costly Functions with Simple Constraints: A Limited-Memory Projected Quasi-Newton Algorithm , 2009, AISTATS.

[9]  Mike E. Davies,et al.  Iterative Hard Thresholding for Compressed Sensing , 2008, ArXiv.

[10]  Yücel Altunbasak,et al.  Approximation and Compression With Sparse Orthonormal Transforms , 2015, IEEE Transactions on Image Processing.

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

[12]  Andrzej Cichocki,et al.  Inpainting color images in learned dictionary , 2012, 2012 Proceedings of the 20th European Signal Processing Conference (EUSIPCO).

[13]  P. Schönemann,et al.  A generalized solution of the orthogonal procrustes problem , 1966 .

[14]  Manuel Nieto Vesperinas,et al.  Scattering And Diffraction In Physical Optics , 2006 .

[15]  Joel A. Tropp,et al.  Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit , 2007, IEEE Transactions on Information Theory.

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

[17]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

[18]  Feng Liu,et al.  Image Denoising Via Sparse and Redundant Representations Over Learned Dictionaries in Wavelet Domain , 2009, 2009 Fifth International Conference on Image and Graphics.

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