High resolution Radon transform based on the reweighted-iterative soft threshold algorithm
XUE Yaru1, GUO Mengjun1, FENG Luyu1, MA Jitao2, CHEN Xiaohong2
1. College of Information Science and Engineering, China University of Petroleum(Beijing), Beijing 102249, China; 2. College of Geophysics, China University of Petroleum(Beijing), Beijing 102249, China
Abstract:The resolution of Radon transform is the key to seismic data processing. The iterative weighting method based on Bayes inversion improves the re-solution of Radon transform, but its convergence rate is low. In light of the strong correlation between Radon transform spaces, the convergence rate of the iterative soft threshold algorithm applied to Radon transform inversion is also low, and the resolution is poor. In this paper, the iterative reweighted least squares algorithm is embedded into the iterative soft threshold algorithm to form a reweighted-iterative soft threshold algorithm. The idea of weighted matrix in high-resolution Radon transform is introduced, and the prior information of Radon parameters is employed to constrain the inversion error function, overcoming the disadvantages of slow convergence and low resolution of the iterative soft threshold algorithm. Synthetic records and real seismic data processing show that this method improves the resolution of Radon transform and achieves good performance in multiple separation and noise suppression.
Claerbout J F,Johnson A G.Extrapolation of time-dependent waveforms along their path of propagation[J].Geophysical Journal of the Royal Astronomical Society, 1971,26(1-4):285-293.
[2]
Hampson D.Inverse velocity stacking for multiple elimination[J].Journal of the Canadian Society of Exploration Geophysicists,1986,22(1):44-45.
[3]
Sacchi M D,Ulrych T.High resolution velocity ga-thers and offset space reconstruction[J].Geophysics,1995,60(4):1169-1177.
[4]
Sacchi M D,Porsani M J.Fast high resolution parabolic Radon transform[C].SEG Technical Program Expanded Abstracts,1999,18:1477-1480.
[5]
Herrmann P.De-aliased,high-resolution radon transforms[C].SEG Technical Program Expanded Abstracts,2000,19:1953-1956.
[6]
刘仕友,马继涛,孙万元,等.基于低频约束的高分辨率Radon变换多次波压制方法研究[J].物探化探计算技术,2019,41(3):293-298.LIU Shiyou,MA Jitao,SUN Wanyuan,et al.Multiple attenuation using low frequency constrained high-resolution parabolic Radon transform[J].Techniques for Geophysical and Geochemical Exploration,2019,41(3):293-298.
[7]
刘喜武,刘洪,李幼铭.高分辨率Radon变换方法及其在地震信号处理中的应用[J].地球物理学进展,2004,19(1):8-15.LIU Xiwu,LIU Hong,LI Youming.High resolution Radon transform and its application in seismic signal processing[J].Progress in Geophysics,2004,19(1):8-15.
[8]
王维红,裴江云,张剑锋.加权抛物Radon变换叠前地震数据重建[J].地球物理学报,2007,50(3):851-859.WANG Weihong,PEI Jiangyun,ZHANG Jianfeng,et al.Prestack seismic data reconstruction using weighted parabolic Radon transform[J].Chinese Journal of Geophysics,2007,50(3):851-859.
[9]
薛亚茹,唐欢欢,陈小宏.高阶高分辨率Radon变换地震数据重建方法[J].石油地球物理勘探,2014,49(1):95-100.XUE Yaru,TANG Huanhuan,CHEN Xiaohong.Reconstruction method of seismic data using high-order high resolution Radon transform[J].Oil Geophysical Prospecting,2014,49(1):95-100.
[10]
王亮亮,毛伟建,唐欢欢,等.快速3D抛物Radon变换地震数据保幅重建[J].地球物理学报,2017,60(7):2801-2812.WANG Liangliang,MAO Weijian,TANG Huanhuan,et al.Amplitude preserved seismic data reconstruction by fast 3D parabolic Radon transform[J].Chinese Journal of Geophysics,2017,60(7):2801-2812.
[11]
薛亚茹,王敏,陈小宏.基于SL0的高分辨率Radon变换及数据重建[J].石油地球物理勘探,2018,53(1):1-7.XUE Yaru,WANG Min,CHEN Xiaohong.High re-solution Radon transform based on SL0 and its application in data reconstruction[J].Oil Geophysical Prospecting,2018,53(1):1-7.
[12]
Cary P W.The simplest discrete Radon transform[C].SEG Technical Program Expanded Abstracts,1998,17:1999-2002.
[13]
Schonewille M,Aaron P,Allen T.Applications of time-domain high-resolution Radon demultiple[C].SEG Technical Program Expanded Abstracts,2007,26:1-5.
[14]
巩向博,韩立国,李洪建.各向异性Radon变换及其在多次波压制中的应用[J].地球物理学报,2014, 57(9):2928-2936.GONG Xiangbo,HAN Liguo,LI Hongjian.Anisotropic Radon transform and its application to demultiple[J].Chinese Journal of Geophysics,2014,57(9):2928-2936.
[15]
Trad D,Ulrych T and Sacchi M.Latest views of the sparse Radon transform[J].Geophysics,2003,68(10):386-399.
[16]
Lu W K.An accelerated sparse time-invariant Radon transform in the mixed frequency-time domain based on iterative 2D model shrinkage[J].Geophysics,2013,78(4):V147-V155.
[17]
Wang B,Zhang Y,Lu W,et al.A robust and efficient sparse time-invariant Radon transform in the mixed time-frequency domain[J].IEEE Transactions on Geoscience and Remote Sensing,2019,57(10):7558-7566.
[18]
巩向博,韩立国,王升超.混合域高分辨率双曲Radon变换及其在多次波压制中的应用[J].石油物探,2016,55(5):711-718.GONG Xiangbo,HAN Liguo,WANG Shengchao.High-resolution hyperbolic Radon transform and its application in multiple suppression[J].Geophysical Prospecting for Petroleum,2016,55(5):711-718.
[19]
Zhu L C,Liu E T,McClellan J H.Seismic data denoising through multiscale and sparsity-promoting dictionary learning[J].Geophysics,2015,80(6):WD45-WD57.
[20]
Zhu L C,Liu E T,McClellan J H.Joint seismic data denoising and interpolation with double-sparsity dictionary learning[J].Journal of Geophysics and Engineering,2017,14(4):802-810.
[21]
李慧,韩立国,张良,等.Radon域下的K-SVD字典的混采分离[J].世界地质,2019,38(1):256-267.LI Hui,HAN Liguo,ZHANG Liang,et al.Separation of blended seismic acquisition with K-SVD dictionary in Radon domain[J].Global Geology,2019,38(1):256-267.
[22]
兰南英,张繁昌,张益明,等.快速结构字典学习三维地震数据重建方法[J].石油地球物理勘探,2020,55(1):1-9.LAN Nanying,ZHANG Fanchang,ZHANG Yiming,et al.3D seismic data reconstruction based on a fast structure dictionary learning method[J].Oil Geophy-sical Prospecting,2020,55(1):1-9.
[23]
王雄文,王华忠.基于压缩感知的高分辨率平面波分解方法研究[J].地球物理学报,2014,57(9):2946-2960.WANG Xiongwen,WANG Huazhong.A research of high-resolution plane-wave decomposition based on compressed sensing[J].Chinese Journal of Geophy-sics,2014,57(9):2946-2960.
[24]
Latif A,Mousa W A.An efficient undersampled high-resolution Radon transform for exploration seismic data processing[J].IEEE Transactions on Geoscience and Remote Sensing,2017,55(2):1010-1024.
[25]
魏亚杰,张盼,许卓.基于稀疏约束反演的三维混采数据分离[J].地球物理学报,2019,62(10):4000-4009.WEI Yajie,ZHANG Pan,XU Zhuo.Separation of 3D blending seismic data based on sparse constrained inversion[J].Chinese Journal of Geophysics,2019,62(10):4000-4009.
[26]
赵子越,李振春,张敏.利用压缩感知技术的离散正交S变换地震数据重建[J].石油地球物理勘探,2020,55(1):29-35.ZHAO Ziyue,LI Zhenchun,ZHANG Min.Seismic data reconstruction using discrete orthonormal S-transform based on compressive sensing[J].Oil Geo-physical Prospecting,2020,55(1):29-35.
[27]
Daubechies I,Defrise M,Mol C D.An iterative thresholding algorithm for linear inverse problems with a sparsity constraint[J].Communications on Pure and Applied Mathematics,2010,57(11):1413-1457.
[28]
Donoho D L.De-noising by soft thresholding[J].IEEE Transactions on Information Theory,2006,41(3):613-627.
[29]
Beck A,Teboulle M.A fast iterative shrinkage-thresholding algorithm for linear inverse problems[J].SIAM Journal on Imaging Sciences,2009,2(1):183-202.
[30]
潘树林,闫柯,李凌云,等.自适应步长FISTA算法稀疏脉冲反褶积[J].石油地球物理勘探,2019,54(4):737-743.PAN Shulin,YAN Ke,LI Lingyun,et al.Sparse-spike deconvolution based on adaptive step FISTA algorithm[J].Oil Geophysical Prospecting,2019,54(4):737-743.
[31]
Scales J A,Gersztenkorn A,Treitel S.Fast Lp solution of large,sparse,linear systems:Application to seismic travel time tomography[J].Journal of Computational Physics,1988,75(2):314-333.
[32]
Chen Z H,Lu W K.Non-iterative high resolution Radon transform[C].Extended Abstracts of 73rd EAGE Conference and Exhibition,2011,4803-4807.
[33]
Daubechies I,Devore R,Fornasier M,et al.Iteratively reweighted least squares minimization for sparse recovery[J].Communications on Pure and Applied Mathematics,2010,63(1):1-38.