高分辨率Radon变换方法及其在地震信号处理中应用.pdfVIP

1 9 1 Vol. 19 No. 1 2004 1(008~ 015) PROGRESS IN GEOP YSICS January 2004 Radon 刘喜武, 刘 洪, 李幼铭 ( , 100029) Radon 变换方法在地震资料处理中广泛采用, 在地震同相轴识别和估计方面具有良好效果. 论是倾斜叠 加,还是广义Radon 变换方法, 一般采用最小二乘反演方法实现. 目前, 在提高反演算法的效率和分辨率方面仍值得 研究.本文从倾斜叠加的定义出发, 阐明Radon 变换分辨率问题的来源和解决办法. 采用最小二乘反演方法研究高分 辨率抛物线Radon 变换和双曲Radon 变换时, 给出稀疏约束预条件共轭梯度法求解的高分辨率Radon 变换的实现方 法, 同阻尼最小二乘方法相比, 分辨率和精度明显提高,文中给出了模型算例.根据有效波和多次波NMO 后剩余时差 不同, 采用高分辨率抛物线和双曲Radon变换可以压制多次波,分别给出了方法原理, 最后给出应用实例. 研究表明, 稀疏约束预条件共轭梯度法可以有效实现高分辨率Radon 变换; 数值算例表明,算法计算效率和精度较高, 可以更好 地实现多次波压制. Radon 变换, 最小二乘反演, 分辨率,稀疏约束, 共轭梯度,地震, 多次波,压制 P315 A 1004-2903(2004)01-0008-08 High resolution radon transform and its application in seismic signal processing LIU X-i wu, LIU ong, LI You-ming (I nstitute of Geology Geophysics, ChineseA cademy of Sciences, Beij ing 100029, China) Abstract Radon transformation is employed commonly in seismic dataprocessing, especially in the events identification and estimation. The least square inversion scheme is the most ordinary method to implement the Radon transform, both for slant stack and for general Radon transform. Although many and many study have been made, there are a lot of works needed doing in enhancing thetransform revolution and the computation efficiency. In this paper ,Parabolic Radon Trans- form ( PRT) and yperbolic Radon Transform( RT) havebeenstudied by least square inversion . In the implementation of the methods, sparse constraint preconditions conjugate gradient method has been applied for enhancing theresolut