Denoising method based on sparse hyperbolic Radon transform in the frequency domain
JIA Chun-Mei1, JIANG Guo-Qing2, LIU Zhi-Cheng1, XU Lu1
1. Research Institute of Geophysical Prospecting, SINOPEC, Nanjing 211103, China;
2. Basic Geological Research Institute, Geological Survey of Jiangsu Province, Nanjing 210018, China
Radon transform is a kind of technical method widely used in seismic data processing,but there are problems of low resolution and low computational efficiency in practical application.On the basis of phase shift hyperbola,the authors deduced time-invariant hyperbolic integral path,proposed sparse hyperbolic Radon transform based on Bayes principle,provided parameter selection principle of discrete transform,and introduced dynamic time windows to ensure the condition that reflection waves of different layers have the best convergence.Sparse hyperbolic Radon transform in the frequency domain has higher resolution and it can be realized in the frequency domain,so the computational efficiency can be drastically improved.Model data test and actual application result show that sparse hyperbolic Radon transform in the frequency domain can effectively remove regular interference,suppress random noise and improve the signal-to-noise ratio;in addition,it also has better fidelity and amplitude preservation.
[1] Claerbout J F, Johnson A G. Extrapolation of time-dependent waveforms along their path of propagation[J]. Geophysical Journal International,1971,26(1):285-293.
[2] Sacchi Mauricio D, Ulrych Tadeusz J. High-resolution velocity gathers and offset space reconstruction[J]. Geophysics,1995,60 (4):1169-1177.
[3] Cary P W. The simplest discrete Radon transform[C]//Expanded Abstracts of 68th SEG Annual Internet Meeting,1998:1999-2002.
[4] Sacchi M D, Porsani M. Fast high-resolution parabolic Radon transforms[C]//Expanded Abstracts of 69th SEG Annual Internet Meeting,1999:1477-1480.
[5] Nowak E, Imhof M. Amplitude preservation of Radon-based multiple-removal filters[J]. Geophysics, 2006, 71(5):123-126.
[6] Yu Z, Ferguson J, McMechan G, et al. Wavelet-Radon domain dealiasing and interpolation of seismic data[J]. Geophysics, 2007, 72(2):41-49.
[7] Wang J F.Greedy least-squares and its application in Radon transforms[C]//2009 CSPG CSEG CWLS Convention, 2009:5-8.
[8] 王维红, 裴江云, 张剑锋. 加权抛物Radon变换叠前地震数据重建[J]. 地球物理学报,2007, 50(3):851-859.
[9] 熊登, 赵伟, 张剑锋. 混合域高分辨率抛物Radon变换及在衰减多次波中的应用[J]. 地球物理学报,2009,52(4):1068-1077.
[10] 石颖, 王维红. 基于波动方程预测和双曲Radon变换联合压制表面多次波[J]. 地球物理学报, 2012,55(9):3115-3125.
[11] 鲁娥, 李庆春. 混合Radon变换地震噪声压制的应用[J]. 物探与化探,2013,37(4):706-710.
[12] 张振波, 轩义华. 高分辨率抛物线拉冬变换多次波压制技术[J]. 物探与化探,2014,38(5):981-988.
[13] Richard J Castle. A theory of normal moveout[J]. Geophysics,1994,59(6):983-999.
[14] 石颖, 柯璇, 王维红. 一种加权抛物Radon变换地震波场重建方法[J]. 物探与化探,2012,36(5):846-850.
[15] 渥·伊尔马滋.地震资料分析——地震资料处理、反演和解释[M].北京:石油工业出版社,2006.
[16] Schonewille M A, Duijndam A J W. Parabolic Radon transform, sampling and efficiency[J]. Geophysics,2001,66(2):667-678.
[17] Trad D, Ulrych T, Sacchi M D. Latest views of the sparse Radon transform[J]. Geophysics,2003, 68:386-399.
[18] Zhang J F, Wang W H. Multiple attenuation: an approach by incorporating multiple prediction with Radon transform[J]. Journal of Seismic Exploration,2006,15:81-99.