应用加权迭代软阈值算法的高分辨率Radon变换

Radon变换的分辨率是其进行地震数据处理的关键因素。基于Bayes反演理论的迭代加权方法改善了Radon变换的分辨率,但其收敛速度较慢;由于Radon变换空间的强相关性,常用的迭代软阈值方法(ISTA)应用于Radon变换反演时的收敛速度也较慢,且分辨率较低。为此,将迭代加权最小二乘法嵌入ISTA中,形成了加权ISTA算法。该方法引入高分辨Radon变换中加权矩阵的思想,利用Radon参数的先验信息约束反演误差函数,克服了ISTA收敛速度慢、分辨率低的缺点。合成记录和实际地震资料处理结果表明,该方法提高了Radon变换分辨率,在多次波分离、噪声压制中效果显著。

线性同相轴波场分离的高分辨率τ-p变换法

基于最小二乘τ-p变换和τ-p域模型稀疏分布的假设,本文给出高分辨率τ-p变换的推导及其模型空间域的离散采样公式,同时给出了保振幅线性同相轴波场分离的算法流程.在求解本文给出的高分辨率τ-p正变换时,由于待求解的矩阵不具备最小平方法所具有的Toeplitz结构,故采用Cholesky分解法进行计算.本文模拟了井间地震和阵列声波测井中的Stoneley上下行波的分离算法过程,高分辨率正反τ-p变换且滤波所得结果显示本文算法误差小和保振幅的特点.对于在τ-p域距离很近或时间域同相轴近于水平的线性波场,高分辨率算法的聚焦作用使得所分离波场畸变小,体现本文算法精度高的优点.理论模型试算表明本文给出的高分辨率τ-p变换线性波场分离算法具有稳定性、精度高和保振幅的特点.

对比抛物Radon正变换几种矩阵求解方法

抛物Radon变换法(Parabolic Radon Transform)在地震资料处理中有广泛的应用。PRT可对不同频率的地震数据解耦处理,这一特点使得抛物Radon变换的计算效率比双曲Radon变换有数量级上的提高。在频率域求解时,需要对每一个频率成份求解同样大小的线性方程组。求解抛物Radon正变换的计算方法主要有Levinson递推法、共轭梯度法、Cholesky分解法和直接矩阵求逆法。最小平方抛物Radon正变换所形成的矩阵具有Toeplitz结构,可采用Levinson递推法进行计算。高分辨率抛物Radon正变换所形成矩阵的Toeplitz结构被破坏,一般采用共轭梯度法或Cholesky分解法进行求解。这里详细推导了复Toeplitz矩阵的Levinson递推算法,并分别对求解方程的四种方法进行了讨论,最后给出抛物Radon正变换求解的数值算例,并对所给出的四种方程求解方法的计算效率及计算精度进行了对比。

面波基模式频散曲线提取解精度对比研究

针对不同方法提取面波基模式频散曲线精度问题,笔者分别采用τ-p变换、频率分解法、F-K变换、高分辨率线性拉东变换(High-Resolution Linear Radon Transform,简称HRLRT)对六层递增型地质模型合成瑞雷波记录进行频散能量成像,并按频散能量最大值拾取基模式频散曲线。为定量评价理论模型基模式频散曲线提取解与解析解的接近程度,引入了均方差与相关系数两种评价参数。评价结果表明,高分辨率线性拉东变换提取基模式频散曲线精度最高,均方差为11.167 8,相关系数为0.994 9;F-K变换提取基模式频散曲线精度最低,均方差为195.274,相关系数为0.515 2。

海上地震资料多次波压制技术研究与应用

本文围绕深海地震资料处理中的最大难点-多次波干扰展开攻关研究。根据海上地震资料多次波传播路径的不同,将其划分为虚反射、表层多次波、微屈多次波、简单多次波、层间多次波5种,又依据它们在地震剖面上是否会产生独立的同向轴将其归纳为短程多次波和长程多次波两大类别,合理透彻地分析了海上多次波的产生机制和特点,总结归纳出当前各种海上地震资料多次波压制方法的优缺点和适用条件:一是"τ-p域反褶积"-虽然远近道多次波都可去除,但会产生空间假频现象;二是"高分辨率Radon分解法"-去长程多次波效果较好,近道无能为力。三是"模型法"-适用于有效构造和多次波干扰形态差异大的地区。四是"表面相关法"-远近道多次波去除效果都不错,但算法过于耗时,现有机器能力不适合大规模使用。