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

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

混合Radon变换地震噪声压制的应用

如何有效地从地震资料中剔除相干噪声一直是地震勘探十分关心的问题,尤其是多次波的压制和消除。笔者采用相似函数压制端点效应和截断效应的混合Radon变换从地震剖面中一次性分离面波等线性噪声,再应用自适应滤波技术在Radon域识别并剔除多次波,并使用双曲Radon反变换计算获得有效反射地震信号。理论模型及实测资料试算结果表明,该方法对多次波压制效果良好,同时,应用混合Radon变换可简化数据中噪声压制处理流程。     

压制截断效应的Radon变换在波场分离中的应用

由于离散运算,Radon变换方法中存在截断效应,传统的最小二乘约束反演并不能有效抑制这种效应。介绍了一种频域空间稀疏约束算法,在反演迭代过程中,根据前一次迭代的结果,通过贝叶斯原理将加权矩阵与前一次迭代的结果联系起来,得到新的加权矩阵;然后求解这个加权矩阵方程,得到频率域的稀疏解。比较了用最小二乘反演和压制截断效应的Radon域变换的效果,后者效果明显优于前者。模拟了一个5层的水平地层的地质模型,用这种压制截断效应的Radon变换进行了VSP资料的上下行波场的分离,得到了较好的效果。     

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

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

二维小波变换与中值滤波联合去噪方法研究

由于随机噪声是一种频带较宽的干扰波,因此依靠单一的去噪处理方式往往难以获得清晰反映目标体的地震信息。小波变换能够较好的去除高斯噪声,保留有效波中、高频成分,提高记录的信噪比,但去除脉冲噪声的效果却并不理想;中值滤波具有良好的边缘保持特性,虽低频去噪声效果有限,但去除脉冲噪声效果明显。因此可利用二维小波变换与中值滤波优势互补的方法,进行叠前去噪处理,达到去除宽频随机噪声的目的。首先运用二维小波变换的理论,采用自适应门限阀值方法进行去噪,同时结合中值滤波方法联合去噪。模型与实际数据的应用效果表明,联合去噪方法可有效压制噪声能量、保留高频有效信号、提高地震记录信噪比与分辨率。     

基于稀疏反演算法的高分辨率Radon变换及其在多次波压制中的应用

针对目前双曲Radon变换Radon域能量团能量泄漏的问题,提出一种基于稀疏约束反演算法的高分辨率Radon变换进行多次波的压制。该方法在常规双曲Radon变换的基础上加入稀疏约束条件,应用共轭导向梯度(CGG)反演算法进行求解,能有效地提高Radon域速度谱的聚焦效果,有利于分离一次波和多次波。模型和实际资料试算证明了该方法的有效性和实用性。     

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

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

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

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

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

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

稳健主成分分析在地震资料异常值噪声压制中的应用

提高地震数据的信噪比是地震资料处理的重要目标之一。传统的地震去噪方法虽然可以有效压制随机噪声,但对非高斯分布的异常值噪声压制效果欠佳。本研究展示了一种基于稳健主成分分析的地震数据异常值噪声压制方法。该方法在频率-空间域通过对地震数据实施稳健低秩近似来求取理想无噪声数据。在目标函数构建方面,采用核范数最小化模型求取理想的低秩近似数据,并使用l₁范数最小化模型来估计异常值噪声。此外,运用增广拉格朗日乘子法求解该反演问题。最后,模型数据和实际资料的去噪结果验证了本研究方法的有效性,与传统F-XY域预测滤波法去噪结果进行对比,也显示本研究方法在有效压制异常值噪声的同时能更好地保护有效波能量。     

基于Radon域阈值迭代的同时震源波场分离

由于同时震源数据包含串扰很难直接成像,因此需要对混叠的炮记录进行分离,为了将混叠的波场进行分离,笔者建立了基于稀疏域反演的阈值迭代模型,将混合采集数据的共检波点道集转换到Radon域进行稀疏表示,利用阈值算子进行约束,通过迭代反演得到分离结果。实验算例表明通过较少的迭代即可获得信噪比较高的分离结果,验证了反演方法的有效性,并且该算法对随机噪声有很强的压制效果。     

利用Radon变换进行三度体重磁异常反演

从Ｒａｄｏｎ变换角度分析解释了重磁三度体异常和二度体异常之间的等价关系，提出了利用Ｒａｄｏｎ正变换将三度体异常转化为二度体异常，对二度体异常做一维反演，然后将反演结果通过反Ｒａｄｏｎ变换实现三维场源图象重建的思想，并付诸实施，在微机上通过模型检验，证明该方法实用、可行。     

深拖平缆采集资料中鬼波衰减方法探讨

鬼波衰减是深拖平缆采集资料处理中非常重要的一步,也是实现宽频处理的关键步骤。为此,笔者分析了深拖平缆资料特点,推导了鬼波的时距曲线方程,在此基础上,介绍了线性Radon域最小平方反演方法衰减鬼波的基本原理,该方法首先通过最小平方方法求得无鬼波的Radon域数据,再应用反变换算子中的一次波算子进行反变换,从而实现了对鬼波的压制。采用该理论对模型数据与实际地震资料进行了试处理,结果表明,深拖平缆资料通过本文的鬼波衰减技术处理后,陷频点得到了补偿,有效频带得到了拓宽,达到了宽频处理的效果。     

面波频散能量谱计算方法

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

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

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

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

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

基于迭代阈值收缩的高分辨率Radon变换方法效果对比

提高Radon变换分辨率的研究一直是地震数据处理方法研究的热点之一.常用的提高分辨率的算法是在频率域进行的,然而此方法在频域中计算的模型权重是偶然耦合的,会对所有同相轴施加相同的权重,导致高能量同相轴地震道生成假像.此文给出了3种时间域提高Radon变换分辨率的方法,IST(iterative shrinkage th...     

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

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

双曲Radon变换参数分析

为避免双曲Radon变换在消除多次波和提取速度差异的过程中出现假频现象,对地震资料处理造成一定的影响。笔者在时间域对最小二乘双曲Radon变换中的曲率参数和最小二乘迭代次数进行了研究,计算出平衡准确度和计算量的最佳参数选取范围为,曲率参数在30~40之间,迭代次数在6~10之间。理论模型测试和实际资料处理结果表明,适当的参数选取能够在保证准确度的同时极大减小计算量,有利于双曲Radon变换的实际应用。