高精度Radon变换方法及去噪分析.pdfVIP

高精度Radon变换方法及去噪分析 郭全仕张卫华黄华昌 李佩 中国石化石油勘探开发研究院南京石油物探研究所，江苏南京，210014 摘要简要介绍了线性、抛物线和双曲线Radon变换方法，以及求解反问题的广义逆算法和两种 高精度算法，即具有较高速度分辨率的稀疏约束反演算法和具有较高波场恢复性能的自适应波场 分离算法；通过台成数据模型比较并证实了上述方法在扫描速度分辨能力和波场恢复方面的优 势。同时将这些方法应用于地面地震的多次渡压制和井问地震的管波消除．取得了令人满意的 效果。 关键词Radon变换稀疏约束反演自适应算法波场分离多次波管波 在—P变换…出现后的几年间，逐渐应用到速度分析方面，如压制多次波、纵横波分离、 偏移和反演等。由于它是采用直线变换，变换后所反映的速度特征不明显。常见的叠加速度 谱分析方法采用了双曲线扫描获取相干能量，但对速度的分辨率较低，而且这个过程是不可逆 利用叠加速度场到偏移距道集反演方法，实现了t一*域的双曲线Radon变换，但由于该算法 涉及到超大规模稀疏矩阵的求逆，计算量过大，即使是现在的计算机也难以满足其运算的需 求，难以应用于实际资料的处理。其后Hampson”o对算法进行改进，引入了抛物线型近似算 法，将通常的双曲线型时距曲线，用没有物理意义的抛物线型参数曲线代替，再将问题转化到 频率域，转变成对单个频率分别求解。这就使求逆矩阵的维数大幅度减小，计算量呈数量级降 低，也使双曲线Radon变换方法能够用于实际地震资料地处理”’“…。 虽然HⅢpson算法极大提高了计算效率，但由于地震信号在时间和空问上都是有限的离 度点不是严格收敛到固定位置，有效波和干扰渡在Radon域相互交叉，如果采用切除曲线直接 分开有效波和干扰波，势必造成这两种波能量间的相互渗透。为了解决这种能量相互渗透造 Radon域自适应分离的高精度算法。作为比较，广义逆算法本文称之为常规算法。 在地震去噪效果分析方面，前人的研究工作oo“主要集中于速度扫描问隔和扫描范围引 起的空间假频问题。本文将讨论前述算法的速度分辨率和波场恢复性能，并通过理论模型和 实际资料进行验证和试验。 1 Radon变换 在频率域，Radon变换可分解为单个频率的变换形式，即 rm=L8d 【d’=Lm 一72～ 度分别为偏移距道数和速度点的采样道数，变换矩阵算子L决定了变换的类型。三种Radon 变换算子￡的元素定义为 ，e。“8 线性变换 毛={e一…何 抛物线变换 (2) Le～’嘶 双曲线变换 这里；是偏移距，∞是变换矩阵对应的角频率，双曲线变换中的∞’是地震道经过l一‘2转换后 的Fourier频率，与速度口有关的扫描参数a定义为 P=1／v．q=1／v2 式(1)实际上定义的是共轭和正交映射问题。算子p为共轭转置算子。从式(2)中可以 看出，该算子只有相位移动，而不改变谱的能量，也没有加强噪音。也就是通过(1)式的正变 换和反变换，最后得到 d，=LL8d 件下变换工才是酉矩阵呢?推算表明，只有采用密集的速度网格(却—加或向—加)，G才接近 单位矩阵f。 遗憾的是，减小速度离散间隔同样会增加计算量，我们不可能无限增加速度采样数量。在 实际计算中，偏移距和速度的采样数总是有限的，并且在某些频率段，G还是病态矩阵。在这 种情形下，通常采用阻尼最小二乘方法求取G的广义逆，取G7=一+G，有 m=LH(心4-L矿)。1d (3) 式中¨是约束因子，在某种意义上也可以将其看成是道集记录中的随机噪音估计。 由于Radon反演问题经常是病态的，它们通常不满足存在性、唯一性及稳定性条件。这时 可以用正则化方法求解，也就是通过添加某些类型的先验信息，使之满足波场重建所需要的唯 一及稳定条件。先验信息经常以规定