Recovery of seismic wavefields based on compressive sensing by an l1-norm constrained trust region method and the piecewise random subsampling

SUMMARY Due to the influence of variations in landform, geophysical data acquisition is usually subsampled. Reconstruction of the seismic wavefield from subsampled data is an ill-posed inverse problem. Compressive sensing (CS) can be used to recover the original geophysical data from the subsampled data. In this paper, we consider the wavefield reconstruction problem as a CS and propose a piecewise random subsampling scheme based on the wavelet transform. The proposed sampling scheme overcomes the disadvantages of uncontrolled random sampling. In computation, an l1-norm constrained trust region method is developed to solve the CS problem. Numerical results demonstrate that the proposed sampling technique and the trust region approach are robust in solving the ill-posed CS problem and can greatly improve the quality of wavefield recovery.

[1]  David L Donoho,et al.  Compressed sensing , 2006, IEEE Transactions on Information Theory.

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

[3]  Xue Feng,et al.  Retrieval of the aerosol particle size distribution function by incorporating a priori information , 2007 .

[4]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .

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

[6]  Naihua Xiu,et al.  Regularizing active set method for nonnegatively constrained ill-posed multichannel image restoration problem. , 2009, Applied optics.

[7]  F. Herrmann,et al.  Simply denoise: Wavefield reconstruction via jittered undersampling , 2008 .

[8]  Stephen P. Boyd,et al.  An Interior-Point Method for Large-Scale $\ell_1$-Regularized Least Squares , 2007, IEEE Journal of Selected Topics in Signal Processing.

[9]  Shiqian Ma,et al.  Projected Barzilai–Borwein method for large-scale nonnegative image restoration , 2007 .

[10]  Changchun Yang,et al.  A review on restoration of seismic wavefields based on regularization and compressive sensing , 2011 .

[11]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[12]  Felix J. Herrmann,et al.  Curvelet-based seismic data processing : A multiscale and nonlinear approach , 2008 .

[13]  Mário A. T. Figueiredo,et al.  Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems , 2007, IEEE Journal of Selected Topics in Signal Processing.

[14]  Changchun Yang,et al.  ON TIKHONOV REGULARIZATION AND COMPRESSIVE SENSING FOR SEISMIC SIGNAL PROCESSING , 2012 .

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

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

[17]  Roger Fletcher,et al.  Projected Barzilai-Borwein methods for large-scale box-constrained quadratic programming , 2005, Numerische Mathematik.

[18]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[19]  Ya-Xiang Yuan,et al.  Convergence and regularity of trust region methods for nonlinear ill-posed inverse problems , 2005 .

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