抗噪性加权抛物线Radon变换的地震道重建

叠前CMP道集中,走时为双曲线的地震反射同相轴经过部分动校正( NMO)之后变为近似抛物线,当道集中缺失某些地震道时,可利用抛物线Radon变换通过迭代计算进行地震道的重建和恢复. 笔者基于抛物线Radon变换的基本原理,利用加权抛物线Radon变换的计算方法,进行叠前地震道的恢复和重建,并且在数据域对其加权系数进行改进和优化,模型数据和实际资料的试算结果表明,改进后的计算方法即使在较强的随机噪声下依然可以稳建同相轴,且避免了常规方法在重建过程中产生虚假同相轴的缺点,与原始的加权抛物线Radon变换的计算方法相比,新方法具有压制随机噪声的优点,即具有较好的抗噪性,计算精度也有所提高,且算法稳健、高效,具有良好的应用前景.     

压制噪声的高分辨率Radon变换法

Radon变换在地震处理中有诸多应用,在地震同相轴识别和估计方面有良好效果。对最小平方意义下反演和基于贝叶斯原理的稀疏约束反演的Radon变换域结果进行了对比,后者的变换效果更好,收敛能量的分辨率高,即压制了离散运算中的截断效应。又通过加噪模型数据,对比了高分辨率Radon变换方法和Radon域内的二维蒙版滤波方法,两者压制规则干扰和随机噪声的效果较好并且结果相似,但后者计算效率较高。在实际数据处理中,比较了FK滤波、高分辨率线性和抛物线Radon变换,以及二维蒙版滤波方法的去噪效果和计算效率,选择二维蒙版滤波方法最优。     

地震资料处理中的Radon变换及反Q滤波

Radon变换利用信号的横向相关性和信号根幅随炮检距变化的特征,提取有效反射信号同相轴,同时提取的零炮检距地震剖面,保留了地震记录的高频信息,达到了提高分辨率的目的。叠前反Q滤波是以FUTTERMAN衰减模型为理论基础,校正地震子波相位拉伸,补偿地震波频率和振幅的损失,达到提高地震资料分辨率和信噪比及改善同相轴连续性的目的。     

基于双曲Radon变换的混采地震数据重建

采用多级中值滤波分离混采数据,可以很好地保持单炮的有效信息;对分离后的单炮数据,采用双曲Radon变换进行地震数据道重建,可以在Radon域获得极高的分辨率从而准确地重建地震数据。单炮道重建时,通过稀疏约束控制Radon变换来提高地震数据的重建精度。双曲Radon变换计算耗时较长,在算子求解时,本文采用快速迭代收缩阈值算法(FISTA),明显加快收敛速度,提高计算效率。模拟数据和实际数据表明,采用多级中值滤波和双曲Radon变换的方法重建混采地震数据可以获得较高的重建精度。     

基于低频约束的高分辨率Radon变换多次波压制方法研究

抛物线Radon变换是地震资料处理中常用的一种多次波压制方法,常规Radon变换由于采用最小二乘算法的原因,对不同的曲率值采用相同的加权值进行求解,分辨率较低。为提高Radon域数据的分辨率,应采用对不同曲率值应用不同的加权值的方法。一般的思路是利用前一次迭代的结果得到加权矩阵,该算法需要迭代进行,运算量较大。这里拟采用低频约束的方式提高Radon变换的分辨率,即利用前一个频率的运算结果对下一个频率的计算进行约束,该算法分辨率高,计算速度快。模拟数据和实际数据的测试表明,本文方法在多次波压制处理中,可以快速有效地完成多次波去除。     

基于脊波域的低信噪比地震资料处理技术

在小波分析的基础上发展起来的脊波变换方法,能更好地处理含线状变化特征的信号。针对地震资料中的同相轴信息,尝试利用脊波变换方法对其进行处理,提高剖面信噪比,突出同相轴信息。在对某工区实际地震资料的处理中发现,处理后的地震剖面同相轴品质及连续性有了明显改善,信噪比增强,分辨率相应提高,体现了该方法相对常规小波分析方法的优越性。     

抛物Radon域地震干涉技术

当地震干涉技术应用在不完整数据体上时,缺失的地震记录会在合成的虚炮集记录上造成大量的假象。为减少假象的产生,笔者将抛物Radon变换与地震干涉技术相结合,在抛物Radon域中应用地震干涉技术。本文通过数值模拟验证了抛物Radon域地震干涉技术的有效性,并将该方法应用在不完整数据体上。结果表明,抛物Radon域地震干涉技术可以有效减少假象的产生,并能够压制缺失的地震数据的影响,提高虚炮集记录的质量。     

Radon域下的K-SVD字典的混采分离

基于地震反演来分离混采数据的方法,笔者提出用K次迭代奇异值分解(K-SVD)来更新Radon域下混叠信号中的字典原子的方法:在同步源采集和地震稀疏反演的背景下,将混合地震信号的分离视为稀疏反演问题,将共检波点域的数据表示在Radon域内,此时有效信号同相轴收敛;对数据阈值滤波后进行分块字典学习,进一步稀疏地表示地震数据;最后,固定字典,更新恢复信号和稀疏系数完成分离。模拟和实际资料处理结果表明:该方法对于混采数据的分离相对中值滤波、小波变换等更有效、分离质量明显提升,可应用于实际混叠数据中。     

线性Radon域地震干涉层间多次波预测方法

为避免时空域地震干涉技术预测层间多次波过程中出现假象干扰,影响地震资料处理。笔者结合线性Radon变换与地震干涉技术,在线性Radon域中利用地震干涉技术预测层间多次波。通过数值模拟验证了线性Radon域地震干涉层间多次波预测方法的有效性,并将该方法应用到复杂速度模型中。结果表明,线性Radon域地震干涉技术可以有效减少时空域地震干涉技术预测层间多次波存在的假象,压缩数据体积大小,提高计算效率,并能够适用于复杂地质环境。     

串联去噪技术在金属矿地震勘探中的应用

金属矿地震勘探中原始数据的信噪比很低,多种干扰混杂在记录中,常规去噪手段难以达到理想效果。针对金属矿地震数据中强线性干扰,设计了f-k 滤波、Radon 变换方法二级衰减; 针对强随机噪声,设计Radon 变换以及Curvelet 域组合变换法二级衰减,组成了f-k 滤波法、Radon 变换及Curvelet 域组合变换法串联去噪方法。应用于金昌金属矿多类混杂强噪声的地震实际数据中,大幅提高了数据的信噪比,经验证本文方法实用、有效。     

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

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

数学变换方法在地震勘探中的应用

数学变换是地震勘探资料处理方法中的重要手段之一,从数学模型的求解到地震数据压缩的很多方面都离不开数学变换。着重概述了现今地震勘探方法中几种常用数学变换的基本原理及方法,包括傅立叶变换、小波变换、K-L变换及R adon变换等,并介绍了各种方法在地震勘探领域的应用及发展趋势。     

面波频散能量谱计算方法

面波频散能量谱成像精度是近地表横波速度计算的重要影响因素,有效提取面波频散曲线是后续反演地层结构及介质物性参数的重要环节。文中论述了τ-p变换法（τ为垂直波慢度,p为水平波慢度）、高分辨率线性Radon变换法、频率分解法、相移法、f-k变换法（f为频率,k为波数）、SFK变换法和矢量波数变换法等7种不同面波频散能量谱计算方法的基本原理;通过上述方法,对三层速度递增水平层状介质模型正演模拟合成Rayleigh波记录进行频散能量谱计算,并与理论频散曲线进行比较;根据所得频散能量谱,分析不同方法在数据处理过程中的特点与应用差异。结果表明：在提取频散曲线过程中,高分辨率线性Radon变换法和矢量波数变换法提取精度较高,f-k变换法提取精度最低,低频信息与高频信息无法有效识别;SFK变换法、矢量波数变换法、高分辨率线性Radon变换法抗噪性能强;τ-p变换法、f-k变换法、相移法计算速度优于其他方法。     

基于F-K和Radon变换的多次波衰减方法

多次波干扰是地震勘探中最重要的问题之一,其衰减技术也是地震数据处理的难题之一。本文主要介绍了多次波的形成、分类及其特征,并在 ProMax 系统上对合成地震记录的多次波进行了压制:首先采用合适的速度对时-空域的地震记录进行 NMO 校正;然后,分别用 F -K 变换和 Radon 变换方法将其变换到 F -K 域和 Radon域,对多次波所占据的主要区域进行充零衰减;最后,通过对应的反变换将其变换到时-空域,从而达到多次波衰减的目的。通过在 ProMax 系统中的应用表明:不同方法对多次波的衰减效果各有优劣,组合滤波方法能达到更好的效果。