Three key factors in seismic data reconstruction based on compressive sensing
Wen Rui1, Liu Guochang1, Ran Yang2
1. State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum(Beijing), Beijing 102249, China; 2. Daqing Logging and Testing Service Company, Daqing Oilfield Company, PetroChina, Daqing, Heilongjiang 163000, China
Abstract:Seismic data reconstruction is significant for seismic data processing and imaging.Reconstruction methods based on compressive sensing are widely used.In these methods,the sparsity transform,the iterative algorithm,and the threshold model affect the final reconstruction performance and computation efficiency.For the sparsity transform,we analyze the influence of Fourier transform,Curvelet transform and Seislet transform in seismic data reconstruction.For the iterative algorithm,we discuss the projection onto convex sets (POCS),the iterative hard thresholding (IHT),Bregman,and the joint reconstruction by sparsity-promoting inversion (JRSI) reconstruction methods and analyze the advantages and disadvantages of these four methods.For the threshold model,we work on the linear threshold model,the exponential threshold model,and the data-driven threshold model.We analyze how these three key factors affect the reconstruction results through synthetic and real data examples.At the end we draw some conclusions and propose suggestions for practical seismic data reconstruction.
温睿, 刘国昌, 冉扬. 压缩感知地震数据重建中的三个关键因素分析[J]. 石油地球物理勘探, 2018, 53(4): 682-693.
Wen Rui, Liu Guochang, Ran Yang. Three key factors in seismic data reconstruction based on compressive sensing. Oil Geophysical Prospecting, 2018, 53(4): 682-693.
Spitz S.Seismic trace interpolation in the F-X domain.Geophysics,1991,56(6):785-794.
[4]
Crawley S.Seismic Trace Interpolation with Nonstationary Prediction-error Filters[D].Stanford University,2000.
[5]
Naghizadeh M,Sacchi M D.F-x adaptive seismic-trace interpolation.Geophysics,2009,74(1):V9-V16.
[6]
Liu G C,Chen X H.Seismic data interpolation using frequency domain complex nonstationary autoregression.Geophysical Prospecting,2018,66(3):478-497.
[7]
Trickett S,Burroughs L,Milton A et al.Rank reduction-based trace interpolation.SEG Technical Program Expanded Abstracts,2010,29:1989-1992.
[8]
Gao J J,Cheng J,Sacchi M D.A new 5D seismic reconstruction method based on a parallel square matrix factorization algorithm.SEG Technical Program Expanded Abstracts,2015,34:3784-3788.
[9]
Donoho D L.Compressed sensing.IEEE Transactions on Information Theory,2006,52(4):1289-1306.
[10]
Candès E J,Tao T.Near optimal signal recovery from random projections:universal encoding strategies.IEEE Transactions on Information Theory,2006,52(12):5406-5425.
[11]
Duijndam A W,Schonewille M A,Hindriks K.Reconstruction of band-limited signals irregularly sampled along one spatial direction.Geophysics,1999,64(2):524-538.
[12]
Trad D,Ulrych T,Sacchi M D.Accurate interpolation with high resolution time-variant Radon transforms.Geophysics, 2002,67(2):644-656.
[13]
Hargreaves N,ver West B,Wombell R et al.Multiple attenuation using an apex-shifted radon transform.SEG Technical Program Expanded Abstracts,2003,22:1929-1932.
[14]
Thomson D,Herrmann F J,Hennenfent G et al.A parallel windowed fast-discrete curvelet transform applied to seismic processing.SEG Technical Program Expanded Abstracts,2006,25:2767-2771.
[15]
刘国昌,陈小宏,郭志峰等.基于Curvelet变换的缺失地震数据插值方法.石油地球物理勘探,2011,46(2):237-246.Liu Guochang,Chen Xiaohong,Guo Zhifeng et al.Missing seismic data rebuilding by interpolation based on Curvelet transform.OGP,2011,46(2):237-246.
[16]
刘财,李鹏,刘洋等.基于Seislet变换的反假频迭代数据插值方法.地球物理学报,2013,56(5):1619-1627.Liu Cai,Li Peng,Liu Yang et al.Iterative interpolation beyond aliasing using Seislet transform.Chinese Journal of Geophysics,2013,56(5):1619-1627.
[17]
Gan S W,Wang S D,Chen Y K et al.De-aliased seismic data interpolation using seislet transform with low-frequency constraint.SEG Technical Program Expanded Abstracts,2015,34:3900-3904.
[18]
张良,韩立国,许德鑫等.基于压缩感知技术的Shearlet变换重建地震数据.石油地球物理勘探,2017,52(2):220-225.Zhang Liang,Han Liguo,Xu Dexin et al.Seismic data reconstruction with Shearlet transform based on compressed sensing technology.OGP,2017,52(2):220-225.
[19]
Bregman L.The method of successive projection for finding a common point of convex sets (Theorems for determining common point of convex sets by method of successive projection).Soviet Math,1965,6(3):487-490.
[20]
Stark H,Oskoui P.High-resolution image recovery from image plane arrays,using convex projections.Journal of the Optical Society of America A (Optics,Image Science and Vision),1989,6(11):1715-1726.
[21]
Abma R,Kabir N.3D interpolation of irregular data with a POCS algorithm.Geophysics,2006,71(6):E91-E97.
[22]
Herrmann F J,Hennenfent G.Non-parametric seismic data recovery with Curvelet frames.Geophysical Journal of International,2008,173(1):233-248.
[23]
Osher S,Mao Y,Dong B et al.Fast Linearized Bregman Iteration for Compressed Sensing and Sparse Denoising.UCLA CAM Report,2008,8-37.
[24]
Ma J W.Improved iterative curvelet thresholding for compressed sensing.IEEE Transactions on Instrumentation and Measurement,2011,60(1):126-136.
[25]
勾福岩,刘财,刘洋等.基于Bregman迭代的复杂地震波场稀疏域插值方法.应用地球物理,2014,11(3):277-288.Gou Fuyan,Liu Cai,Liu Yang et al.Complex seismic wavefield interpolation based on the Bregman iteration method in the sparse transform domain.Applied Geophysics,2014,11(3):277-288.
[26]
Bai L S,Lu H Y,Liu Y et al.A fast joint seismic data reconstruction by sparsity-promoting inversion.Geophysical Prospecting,2017,65(4):926-940.
[27]
Galloway E,Sacchi M D.POCS method for seismic data reconstruction of irregularly sampled data.CSPG/CSPE Geo Convention,2007,555.
[28]
张华,陈小宏.基于Jitter采样和曲波变换的三维地震数据重建.地球物理学报,2013,56(5):1637-1649.Zhang Hua,Chen Xiaohong.Seismic data reconstruction based on Jittered sampling and curvelet transform.Chinese Journal of Geophysics,56(5):1637-1649.
[29]
王本锋,陈小宏,李景叶等.POCS联合改进的Jitter采样理论曲波域地震数据重建.石油地球物理勘探,2015,50(1):20-28.Wang Benfeng,Chen Xiaohong,Li Jingye et al.Seismic data reconstruction based on POCS and improved Jittered sampling in the curvelet domain.OGP,2015,50(1):20-28.
[30]
Gao J J,Stanton A,Naghizadeh M et al.Convergence improvement and noise attenuation consideration for beyond alias projection onto convex sets reconstruction.Geophysical Prospecting,2013,61(S1):138-151.
[31]
葛子建,李景叶,潘树林等.一种快速收敛的抗噪POCS地震数据重构方法.应用地球物理,2015,12(2):169-178.Ge Zijian,Li Jingye,Pan Shulin et al.A fast-convergence POCS seismic denoising and reconstruction method.Applied Geophysics,2015,12(2):169-178.
[32]
Clearbout J.Image Estimation by Example:Geophysical Soundings Image Construction Multidimensional Auto Regression.Stanford University,2008.
[33]
Blumensath T,Davies M E.Iterative thresholding for sparse approximations.Journal of Fourier Analysis and Applications,2008,14(5-6):629-654.
[34]
Osher S,Burger M,Goldfarb D et a1.An iterated re-gularization method for total variatian-based image restoration.Siam Journal on Multiscale Modeling & Simulation,2005,4(2):460-489.
[35]
Xu J,Osher S.Iterative regularization and nonlinear inverse scale space applied to wavelet-based denoi-sing.IEEE Transactions on Image Processing,2006,12(2):534-544.
[36]
Burger M,Gilboa G,Osher S et al.Nonlinear inverse scale space methods.Communications in Mathematical Sciences, 2006,4(1):175-208.
[37]
He L,Chang T C,Osher S et al.MR Image Reconstruction by Using the Iterative Refinement Method and Nonlinear Inverse Scale Space Methods.UCLA CAM Report,2006,6-35.
Candès E,Demanet L,Donoho D et al.Fast discrete curvelet transforms.Multiscale Modeling & Simulation,2005, 5(3):861-899.
[40]
Fomel S.Towards the seislet transform.SEG Technical Program Expanded Abstracts,2006,25:2847-2851.
[41]
王伟,陈双廷,王宝彬等.五维规则化技术研究与应用.石油地球物理勘探,2017,52(增刊1):28-33.Wang Wei,Chen Shuangting,Wang Baobin et al.Application of 5D regularization in seismic data proces-sing.OGP,2017,52(S1):28-33.
[42]
Hennenfent G,Herrmann F J.Seismic denoising with nonuniformly sampled curvelets.Computing in Science and Engineering,2006,28(3):16-25.