Radon transform based on greedy fast iterative shrinkage threshold and its application in multiple suppression
ZHANG Quan1,2,3, LEI Qin1, LIN Baiyue1, PENG Bo1,2, LIU Shuyan1,2
1. School of Computer Science, Southwest Petroleum University, Chengdu, Sichuan 610500, China; 2. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University), Chengdu, Sichuan 610500, China; 3. School of Information and Communication Engineering, University of Electronic Science and Technology of China, Chengdu, Sichuan 611731, China
Abstract:In seismic exploration, multiples seriously affect the interpretation accuracy of seismic data, and effective suppression of multiples is important in seismic data processing. The parabolic Radon transform is a common method to suppress multiples. The iterative shrinkage thresholding algorithm (ISTA) is the most widely used method in the industry to obtain the solution to the inverse problem of the parabolic Radon transform. It has excellent computational accuracy and efficiency, but for massive seismic data, the processing efficiency still needs to be improved. To improve the convergence rate of the parabolic Radon transform, this study proposes greedy fast ISTA (Greedy FISTA) to processing of inversion problem for Radom transform suppressing multiple, and constract an accelerated sparse time-invariant Randon transform in the mixed frequency-time domain based on fast interative shrinkage-thesholding algorithn(SRTGFIS). Unlike ISTA, Greedy FISTA takes the weighted sum of the results of the previous two iterations as the iteration starting point, and it introduces restart conditions and convergence conditions to reduce the oscillation period in the iteration process and accelerate the calculation. The multiple suppression experiments with synthetic and real data show that compared with ISTA and FISTA, the proposed algorithm has a great improvement in convergence efficiency and a slight improvement in convergence accuracy.
张全, 雷芩, 林柏栎, 彭博, 刘书妍. 基于贪婪—快速迭代收缩阈值的Radon变换及其在多次波压制中的应用[J]. 石油地球物理勘探, 2022, 57(6): 1332-1341.
ZHANG Quan, LEI Qin, LIN Baiyue, PENG Bo, LIU Shuyan. Radon transform based on greedy fast iterative shrinkage threshold and its application in multiple suppression. Oil Geophysical Prospecting, 2022, 57(6): 1332-1341.
RADON J. On the determination of functions from their integrals along certain manifolds[J]. Mathematische Physikalische, 1917,69: 262-267.
[2]
CLAERBOUT J F,JOHNSON A G. Extrapolation of time-dependent waveforms along their path of propagation[J]. Geophysical Journal International,1971,26(1-4): 285-293.
[3]
THORSON J R,CLAERBOUT J F. Velocity-stack and slant-stack stochastic inversion[J]. Geophysics,1985,50(12): 2727-2741.
[4]
李远钦,杜友福. 双曲Radon变换及其在地震资料处理中的应用[J]. 江汉石油学院学报,1990,12(2): 16-20.LI Yuanqin,DU Youfu. Hyperbola Radon transform and its application in seismic data processing[J]. Journal of Jianghan Petroleum Institute,1990,12(2): 16-20.
[5]
SACCHI M D,ULRYCH T J. High-resolution velo-city gathers and offset space reconstruction[J]. Geophysics,1995,60(4): 1169-1177.
[6]
SACCHI M D,PORSANI M. Fast high resolution parabolic Radon transform[C]. SEG Technical Program Expanded Abstracts,1999,18: 1477-1480.
[7]
牛滨华,孙春岩,张中杰,等. 多项式Radon变换[J]. 地球物理学报,2001,44(2): 263-271.NIU Binhua,SUN Chunyan,ZHANG Zhongjie,et al. Polynomial Radon transform[J]. Chinese Journal of Geophysics,2001,44(2): 263-271.
[8]
刘喜武,刘洪,李幼铭. 高分辨率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.
[9]
ETHAN J N,MATTHIAS G I. Amplitude preservation of Radon-based multiple-removal filter[J]. Geophysics,2006,71(5): 123-12.
[10]
刘保童,朱光明. 一种频率域提高Radon变换分辨率的方法[J]. 西安科技大学学报,2006,26(1): 112-116,129.LIU Baotong,ZHU Guangming. An approach to improve the resolution of Radon transform in frequency domain[J]. Journal of Xi’an University of Science and Technology,2006,26(1): 112-116,129.
[11]
SCHONEWILLE M,AARON P A. Applications of time-domain high-resolution Radon demultiple[J]. SEG Technical Program Expanded Abstracts,2007,26: 2565-2569.
[12]
薛亚茹,陈小宏,马继涛. 多方向正交多项式变换压制多次波[J]. 地球物理学报,2012,55(10): 3450-3458.XUE Yaru,CHEN Xiaohong,MA Jitao. Multiples attenuation based on directional orthogonal polyno-mial transform[J]. Chinese Journal of Geophysics,2012,55(10): 3450-3458.
[13]
LU Wenkai. An accelerated sparse time-invariant Radon transform in the mixed frequency-time domain based on iterative 2D model shrinkage[J]. Geophy-sics,2013,78(4): V147-V155.
[14]
薛亚茹,杨静,钱步仁. 基于高阶稀疏Radon变换的预测多次波自适应相减方法[J]. 中国石油大学学报(自然科学版),2015,39(1): 43-49.XUE Yaru,YANG Jing,QIAN Buren. Multiples prediction and adaptive subtraction with high-order sparse Radon transform[J]. Journal of China University of Petroleum (Edition of Natural Sciences),2015,39(1): 43-49.
[15]
SUN Wenzhi,LI Zhenchun,QU Yingming,et al. Multiple attenuation using λ-f domain high-order and high-resolution Radon transform based on SL0 norm[J]. Applied Geophysics,2019,16(4): 473-482.
[16]
薛亚茹,郭蒙军,冯璐瑜,等. 应用加权迭代软阈值算法的高分辨率Radon变换[J]. 石油地球物理勘探,2021,56(4): 736-744,757.XUE Yaru,GUO Mengjun,FENG Luyu,et al. High resolution Radon transform based on the reweighted-iterative soft threshold algorithm[J]. Oil Geophysical Prospecting,2021,56(4): 736-744,757.
[17]
孙宁娜,曾同生,戴晓峰,等. 基于多窗口联合优化的多次波自适应相减方法[J]. 石油物探,2022,61(3): 463-472.SUN Ningna,ZENG Tongsheng,DAI Xiaofeng,et al. Adaptive multiple subtraction based on multi-windows joint optimization[J]. Geophysical Prospecting for Petroleum,2022,61(3): 463-472.
[18]
马继涛,刘仕友,廖震. 三维高精度保幅Radon变换多次波压制方法[J]. 石油地球物理勘探,2022,57(3): 582-592.MA Jitao,LIU Shiyou,LIAO Zhen. Research on multiple attenuation using 3D high-precision amplitude-preserving Radon transform[J]. Oil Geophysical Prospecting,2022,57(3): 582-592.
[19]
HAMPSON D. Inverse velocity stacking for multiple elimination[C]. SEG Technical Program Expanded Abstracts,1986,5: 422-424.
[20]
熊登,赵伟,张剑锋. 混合域高分辨率抛物Radon变换及在衰减多次波中的应用[J]. 地球物理学报,2009,52(4): 1068-1077.XIONG Deng,ZHAO Wei,ZHANG Jianfeng. Hybrid-domain high-resolution parabolic Radon transform and its application to demultiple[J]. Chinese Journal of Geophysics,2009,52(4): 1068-1077.
[21]
ABBAD B,URSIN B,PORSANI M J. A fast,modified parabolic Radon transform[J]. Geophysics,2011,76(1): V11-V24.
刘小舟,胡天跃,刘韬,等. 数据增广的编解码卷积网络地震层间多次波压制方法[J]. 石油地球物理勘探,2022,57(4): 757-767.LIU Xiaozhou,HU Tianyue,LIU Tao,et al. Seismic internal multiple suppression method with encoder-decoder convolutional network based on data augmentation[J]. Oil Geophysical Prospecting,2022,57(4): 757-767.
[24]
LIANG J,SCHNLIEB C B. Improving FISTA: faster,smarter and greedier[DB/OL]. (2021-01-20). https://arxiv.org/abs/1811.01430v1.
[25]
曾华会,王孝,杨维,等. 分级组合多次波压制技术——以玛湖地区为例[J]. 石油地球物理勘探,2018,53(增刊2): 13-19.ZENG Huahui,WANG Xiao,YANG Wei,et al. Stepped-combination multiples suppression: an exam-ple in Mahu Area[J]. Oil Geophysical Prospecting,2018,53(S2): 13-19.
[26]
LI Z X,LI Z C. Accelerated parabolic Radon domain 2D adaptive multiple subtraction with fast iterative shrinkage thresholding algorithm and its application in parabolic Radon domain hybrid demultiple method[J]. Journal of Applied Geophysics,2017,143(1): 86-102.