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

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

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

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

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

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

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

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

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

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

Radon变换多次波压制方法及应用研究

多次波在地震资料中是普遍存在的。Radon变换是现行商业地震资料多次波压制处理软件应用最多的模块之一。这里在详细阐述Radon变换基本原理和离散采样计算的基础上,给出了理论模拟地震数据的多次波压制算法和计算结果。进而将Radon变换多次波压制方法应用于松辽盆地地震资料处理中,在对实际的CRP道集合理设计滤波切除函数后,应用本文的Radon变换算法,得到多次波压制的地震剖面。对比多次波压制前、后的叠前时间偏移地震剖面,压制多次波后地震剖面的信噪比很高,地质构造特征清晰,可应用于后续的地质构造解释,地震属性提取和油气储层预测。该方法的计算效率较高,算法易于实现。理论和实际地震数据的多次波压制结果表明,该算法具有高精度和实用性强的特点。     

应用小波变换和C3相干算法检测地震剖面同相轴

为了提高地震剖面的信噪比和分辨率，使反射同相轴容易被识别和追踪，这里提出了一种应用小波算法及C3相干算法检测地震反射同相轴的新方法。利用小波算法提高地震剖面分辨率，减少大部分随机噪声，再运用C3相关算法和倾角扫描法检测地震剖面同相轴。应用该方法对模型数据进行了试算，证明了可行性，然后检测实际地震剖面同相轴，处理后的地震剖面同相轴品质及连续性都有了明显改善，信噪比增强，分辨率相应提高。     

频域稀疏双曲Radon变换去噪方法

Radon变换是地震数据处理中广泛应用的一种技术方法,但实际应用中存在分辨率及运算效率低的问题。笔者在相移双曲线的基础上,推导了时不变双曲积分路径,提出了基于贝叶斯原理的稀疏双曲Radon变换,给出离散变换的参数选取准则,并且为保证不同反射层的反射波具有最佳收敛,在变换中引入了动态时窗。频域稀疏双曲Radon变换不仅具有较高的分辨率,而且可以在频域实现,计算效率大幅提高。模型数据测试和实际资料应用结果表明,频域稀疏双曲Radon变换可以有效去除规则干扰波、压制随机噪声,提高数据信噪比,且具有较好的保真度和保幅性。     

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

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

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

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

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

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

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

数学变换是地震勘探资料处理方法中的重要手段之一,从数学模型的求解到地震数据压缩的很多方面都离不开数学变换。着重概述了现今地震勘探方法中几种常用数学变换的基本原理及方法,包括傅立叶变换、小波变换、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系统中的应用表明:不同方法对多次波的衰减效果各有优劣,组合滤波方法能达到更好的效果。     

拉东域基于最小二乘反演的斜缆数据鬼波压制方法

斜缆采集资料中鬼波的陷频特征限制了地震记录的频带宽度和分辨率,给地震资料反演、解释带来困难。通过压制鬼波,可获得高分辨率的宽频数据。基于拉东域的水平拖缆鬼波压制方法,结合斜缆鬼波随偏移距变化的特点,推导出一次波和鬼波在频率-拉东域的逆变换算子,建立检波器处总波场和海表面处上行波场之间新的关系式,利用最小二乘反演精确求解获得海表面处上行波场,并延拓得到拖缆处鬼波压制后的记录。通过考虑鬼波延迟时间受偏移距和出射角的影响,弥补鬼波延迟时间估计中存在的误差,无需进行反演迭代求取最优鬼波延迟时间,提高了计算效率。合成数据及海上实际斜缆数据测试结果表明,研究方法能较好地压制鬼波,达到拓宽地震记录频带的目的。