Relevance Vector Machines with Uncertainty Measure for Seismic Bayesian Compressive Sensing and Survey Design

Seismic data acquisition in remote locations involves sampling using regular grids of receivers in a field. Extracting the maximum possible information from fewer measurements is cost-effective and often necessary due to malfunctions or terrain limitations. Compressive Sensing (CS) is an emerging framework that allows reconstruction of sparse signals from fewer measurements than conventional sampling rates. In seismic CS, the utilization of sparse solvers has proven to be successful, however, algorithms lack predictive uncertainties. We apply the Relevance Vector Machine (RVM) to seismic CS and propose a novel utilization of multi-scale dictionaries of basis functions that capture different variations in the data. Furthermore, we propose the use of a new predictive uncertainty measure using the information from the neighbours of each estimation to produce accurate uncertainty maps. We apply the RVM to different seismic signals and obtain state-of-the-art reconstruction accuracy. Using the RVM and its predictive uncertainty map, it is possible to quantify risk associated with seismic data acquisition and at the same time guide future survey design.

[1]  George Eastman House,et al.  Sparse Bayesian Learning and the Relevan e Ve tor Ma hine , 2001 .

[2]  Mauricio D. Sacchi,et al.  Accurate interpolation with high-resolution time-variant Radon transforms , 2002 .

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

[4]  Michael E. Tipping,et al.  Analysis of Sparse Bayesian Learning , 2001, NIPS.

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

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

[7]  王彦飞,et al.  Accelerating seismic interpolation with a gradient projection method based on tight frame property of curvelet , 2015 .

[8]  A. Stanton,et al.  Mitigating Artifacts in Projection Onto Convex Sets Interpolation , 2015 .

[9]  Michael E. Tipping,et al.  Fast Marginal Likelihood Maximisation for Sparse Bayesian Models , 2003 .

[10]  Nadia Kreimer,et al.  A Comparison of 5D Reconstruction Methods , 2012 .

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

[12]  Carl E. Rasmussen,et al.  Healing the relevance vector machine through augmentation , 2005, ICML.

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

[14]  David Zhang,et al.  FSIM: A Feature Similarity Index for Image Quality Assessment , 2011, IEEE Transactions on Image Processing.

[15]  Balas K. Natarajan,et al.  Sparse Approximate Solutions to Linear Systems , 1995, SIAM J. Comput..

[16]  M. Oristaglio SEAM Phase II—Land Seismic Challenges , 2012 .

[17]  Lawrence Carin,et al.  Bayesian Compressive Sensing , 2008, IEEE Transactions on Signal Processing.

[18]  Mauricio D. Sacchi,et al.  High-resolution prestack seismic inversion using a hybrid FISTA least-squares strategy , 2013 .

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

[20]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

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

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

[23]  Michael P. Friedlander,et al.  Probing the Pareto Frontier for Basis Pursuit Solutions , 2008, SIAM J. Sci. Comput..