Seismic trace reconstruction approach based on noise immunity weighted parabolic radon transform
JING Shuang-Wei1, CAO Si-Yuan1, ZHANG Hao-Ran1, HE Hai-Lei2, YUAN Dian1, HE Yuan1, CUI Zhen1, ZHANG Jian-Hong1, GUO Qi1
1. School of Geophysics and Information Engineering, China University of Petroleum(Beijing), Beijing 102200, China;
2. Changqing Geophysical Department of Eastern Geophysical Company, Xi'an 710000, China
In the pre-stack CMP gathers,the events of seismic reflection approximate a hyperbola.After a partial NMO correction,the seismic events approximate parabola.When the gathers miss some seismic traces,parabolic Radon transform algorithm can be used by the iterative to reconstruct and recover the missing seismic traces.Based on the basic principle of parabolic Radon transform,the authors put forward the calculation method of weighted parabolic Radon transform for the reconstruction and recovery of missing pre-stack seismic traces,and improved and optimized the weighted coefficients in the data domain.Test results of model data and real data show that the improved calculation method still has a robust reconstruction result even under the strong random noise,and the method does not produce false events,compared with the calculation method of the original weighted parabolic Radon transform,the method has the advantage of suppressing random noise,this is,the advantage of noise immunity.In addition,the calculation accuracy is also improved,and the algorithm has the characteristics of steadiness and high efficiency.
[1] 张保卫.Radon变换及其在地震数据处理中的应用[D].长安大学,2007.[2] Larner K,Gibson B,Rothman D.Trace interpolation and the design of seismic surveys[J].Geophysics,1981,46(4):407.[3] Spitz S.Seismic trace interpolation in F-X domain[J].Geophysics,1991,56(6):785-794.[4] Bardan V.Trace interpolation in seismic data processing[J].Geophysical Processing,1987,35(3):343-358.[5] Schonewill M A,Duijndam A J W.Parabolic Radon transform,sampling and efficiency[J].Geophysics,2001,66(2):667-678.[6] 黄新武,吴津,牛滨华.基于抛物线拉东变换的地震道重构[J].中国矿业大学学报,2003,32(5):534-539.[7] 黄新武,吴津,牛滨华.抛物线Radon变换中的参数采样与假频[J].石油大学学报:自然科学版,2003,27(2):27-32.[8] Sacchi M D,Porsanio M.Fasthigh resolution parabolic Radon Transform[C]//Expanded Abstracts of 69th SEG Annual Internet Meeting,1999.[9] 刘喜武,刘洪高.高分辨率Radon变换方法及其在地震信号处理中的应用[J].地球物理学进展,2004,19(1):8-15.[10] 王维红,刘洪.道均衡抛物线Radon变换法地震道重建[J].石油地球物理勘探,2005,40(5):518-522.[11] 王维红,裴江云,张剑锋.加权抛物Radon变换叠前地震数据重建[J].地球物理学报,2007,50(3):851-859.[12] 石颖,柯璇,王维红.一种加权抛物线Radon变换地震波场重建方法[J].物探与化探,2012,36(5):846-850.[13] 姜岩,王维红.对比抛物Radon正变换几种矩阵求解方法[J].物探化探计算技术,2007,29(6):555-559.[14] 王维红,刘洪.抛物线Radon变换法近偏移距波场外推[J].地球物理学进展,2005,20(2):289-293.[15] 曾有良,乐友善,单启铜.基于高分辨率Radon变换的VSP波场分离方法[J].石油物探,2007,46(2):115-119.