Compressive sensing method with Huber norm minimization constraint on reconstruction error
LI Zhong-xiao1, LI Yongqiang2,3, GU Bingluo2,3, LI Zhenchun2,3
1. School of Electronic Information, Qingdao University, Qingdao, Shandong 266071, China; 2. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 3. Laboratory for Marine Resources, Qingdao National Laboratory for Marine Science and Techno-logy, Qingdao, Shandong 266071, China
Abstract:Seismic data contain strong noise outliers.The compressive sensing(CS) method based on L2 norm minimization constraint on reconstruction error assumes that the reconstruction error satisfies Gaussian distribution.Therefore,the CS method above cannot remove super-Gaussian noise outliers.To better remove outliers and improve interpolation accuracy,Huber norm was used instead of L2 norm to implement the minimization constraint on reconstruction error.The minimization constraint of Huber norm is equivalent to the L1 norm minimization constraint on large reconstruction error (noise outlier) and the L2 norm minimization constraint on small reconstruction error (Gaussian random noise).Therefore,the proposed method is robust when dealing with noise outlier.Furthermore,theoretical pseudo seismic data were introduced to convert the Huber norm minimization problem to the L2 norm minimization problem,in order to solve the Huber-L0 minimization problem of the proposed CS method based on Huber norm minimization constraint on construction error.Additionally,the affection of Gaussian noise intensity,noise outlier intensity and parameter selection on interpolation accuracy is tested.The processing results of synthetic and field data demonstrated that the proposed CS method based on the Huber norm minimization constraint of the construction error can better remove the noise outliers and preserve effective signals compared with the CS method based on the L2 norm minimization constraint of the construction error.
Abma R,Kabir N.3D interpolation of irregular data with a POCS algorithm[J].Geophysics,2006,71(6):E91-E97.
[3]
Berkhout A J,Verschuur D J.Estimation of multiple scattering by iterative inversion,Part I:Theoretical consideration[J].Geophysics,1997,62(5):1586-1595.
[4]
Li Z X,Lu W K.Adaptive multiple subtraction based on 3D blind separation of convolved mixtures[J].Geophysics,2013,78(6):V251-V266.
[5]
Li Z X,Li Z C,Lu W K.Multichannel predictive deconvolution based on the fast iterative shrinkage thresholding algorithm[J].Geophysics,2016,81(1):V17-V30.
[6]
Li Z X,Li Z C.Accelerated 3D blind separation of convolved mixtures based on the fast iterative shrin-kage thresholding algorithm for adaptive multiple subtraction[J].Geophysics,2018,83(2):V99-V113.
[7]
Plessix R E,Muler W A.Frequency-domain finite-difference amplitude-preserving migration[J].Geophysical Journal International,2004,157(3):975-987.
[8]
Hamid H,Pidlisecky A.Multitrace impedance inversion with lateral constraints[J].Geophysics,2015,80(6):M101-M111.
[9]
Candès E J,Romberg J K,Tao T.Stable signal recovery from incomplete and inaccurate measurements[J].Communications on Pure and Applied Mathematics,2006,59(8):1207-1223.
[10]
Donoho D L.Compressed sensing[J].IEEE Transactions on Information Theory,2006,52(4):1289-1306.
[11]
Ma J W.Compressed sensing by iterative thresholding of geometric wavelets:a comparing study[J].International Journal of Wavelets,Multiresolution and Information Processing,2011,9(1):63-77.
[12]
白兰淑,刘伊克,卢回忆,等.基于压缩感知的Curvelet域联合迭代地震数据重建[J].地球物理学报,2014,57(9):2937-2945.BAI Lanshu,LIU Yike,LU Huiyi,et al.Curvelet-domain joint iterative seismic data reconstruction based on compressed sensing[J].Chinese Journal of Geophysics,2014,57(9):2937-2945.
[13]
Gan S W,Wang S D,Chen Y K,et al.Compressive sensing for seismic data reconstruction via fast projection onto convex sets based on seislet transform[J].Journal of Applied Geophysics,2016,130(5):194-208.
[14]
Wang B F,Wu R S,Chen X H,et al.Simultaneous seismic data interpolation and denoising with a new adaptive method based on dreamlet transform[J].Geophysical Journal International,2015,201(2):1180-1192.
[15]
陈祖斌,王丽芝,宋杨,等.基于压缩感知的小波域地震数据实时压缩与高精度重构[J].石油地球物理勘探,2018,53(4):674-681.CHEN Zubin,WANG Lizhi,SONG Yang,et al.Seismic data real-time compression and high-precision reconstruction in the wavelet domain based on the compressed sensing[J].Oil Geophysical Prospecting,2018,53(4):674-681.
[16]
Wang B F.An efficient POCS interpolation method in the frequency-space domain[J].IEEE Geoscience and Remote Sensing Letters,2016,13(9):1384-1387.
[17]
王新全,耿瑜,WU Rushan,等.基于压缩感知的Dream-let域数据重构方法及应用[J].石油地球物理勘探,2015,50(3):399-404.WANG Xinquan,GENG Yu,WU Rushan,et al.Seismic data reconstruction in Dreamlet domain based on compressive sensing[J].Oil Geophysical Prospecting,2015,50(3):399-404.
[18]
张良,韩立国,许德鑫,等.基于压缩感知技术的Shearlet变换重建地震数据[J].石油地球物理勘探,2017,52(2):220-225.ZHANG Liang,HAN Liguo,XU Dexin,et al.Seismic data reconstruction with Shearlet transform based on compressed sensing technology[J].Oil Geophysical Prospecting,2017,52(2):220-225.
[19]
刘财,李鹏,刘洋,等.基于Seislet变换的反假频迭代数据插值方法[J].地球物理学报,2013,56(5):1619-1627.LIU Cai,LI Peng,LIU Yang,et al.Iterative data interpolation beyond aliasing using Seislet transform[J].Chinese Journal of Geophysics,2013,56(5):1619-1627.
[20]
Chen K,Sacchi M D.Robust reduced-rand filtering for erratic seismic noise attenuation[J].Geophysics,2015,80(1):V1-V11.
[21]
Raphael S,Viguier G,Gondoin R,et al.Multidimensional simultaneous random plus erratic noise attenuation and interpolation for seismic data by joint low-rank and sparse inversion[J].Geophysics,2015,80(6):WD129-WD141.
[22]
Zhao Q,Du Q Z,Gong X F,et al.Signal-preserving erratic noise attenuation via iterative robust sparsity-promoting filter[J].IEEE Transactions on Geoscience and Remote Sensing,2018,56(6):3547-3560.
[23]
Guitton A,Symes W W.Robust inversion of seismic data using the Huber norm[J].Geophysics,2003,68(4):1310-1319.
[24]
Huber P J,Ronchetti E M.Robust Statistics(Second Edition)[M].John Wiley & Sons,New Jersey,2011,693.
[25]
Wong R K W,Lee T C M.Matrix completion with noisy entries and outliers[J].Journal of Machine Learning Research,2017,18(1):1-25.
[26]
Daubechies I,Defrise M,De M C.An iterative thre-sholding algorithm for linear inverse problems with a sparsity constraint[J].Communications on Pure and Applied Mathematics,2004,57(11):1413-1457.
[27]
温睿,刘国昌,冉扬.压缩感知地震数据重建中的三个关键因素分析[J].石油地球物理勘探,2018,53(4):682-693.WEN Rui,LIU Guochang,RAN Yang.Three key factors in seismic data reconstruction based on compressive sensing[J].Oil Geophysical Prospecting,2018,53(4):682-693.
[28]
Wang F J,Hua Q S,Liu C S.Boundary function me-thod for inverse geometry problem in two-dimensional anisotropic heat conduction equation[J].Applied Mathematics Letters,2018,doi:10.1016/j.aml.2018.05.004.
[29]
Holland P W, Welsh R E.Robust regression using iteratively reweighted least-squares[J].Communications in Statistics-Theory and Methods,1977,6(9):813-827.
[30]
Huber P J.Robust estimation of a location parameter[J].The Annals of Mathematical Statistics,1964,35(1):73-101.
[31]
Gao J J,Stanton A,Naghizadeh M,et al.Convergence improvement and noise attenuation considerations for beyond alias projection non convex sets reconstruction[J].Geophysical Prospecting,2013,61(S1):138-151.
[32]
Hennenfent G,Herrmann F J.Simply denoise:Wavefield reconstruction via jittered undersampling[J].Geo-physics,2008,73(3):V19-V28.
[33]
Gao J J,Chen X H,Li J Y,et al.Irregular seismic data reconstruction based on exponential threshold model of POCS method[J].Applied Geophysics,2010,7(3):229-238.
[34]
Yang P L,Gao J H,Chen W C.On analysis-based two-step interpolation methods for randomly sampled seismic data[J].Computers & Geosciences,2013,51(2):449-461.