高分辨率抛物线拉冬变换多次波压制技术

对于孔径过小,采样不充分的地震数据,抛物线拉冬变换的多次波压制效果不理想,且存在空间假频问题。综合Sacchi、Mauricio、Todd Mojesky等对抛物线拉冬变换压制多次波的分析,阐明了高分辨率抛物线拉冬变换原理。通过理论模型试验和实际地震数据的分析证明,高分辨率抛物线拉冬变换可以突破离散采样和有限孔径的限制,实现对地震数据中的多次波及空间假频的快速、有效压制。     

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

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

抛物线拉冬变换算法在近偏移距道数据恢复中的应用

在地震数据处理中,恢复缺失的近偏移距道数据是一项非常重要的工作。本文简述了一种利用抛物线拉冬变换算法来实现空间插值的方法。该方法首先要对原始的共炮点道集（或ＣＭＰ道集）进行动校正,然后通过选择正确的曲率参数对包含零偏移距道的数据作抛物线拉冬变换,再经反变换之后就可将原来缺失的数据正确地恢复出来。文中用理论合成记录验证了该方法的有效性。     

利用小波变换进行地震数据压缩

根据Ｍ带小波变换理论,结合地震数据的特征,将Ｍ带小波变换应用于地震数据的压缩处理,提出了地震数据的压缩编码算法。实验结果表明,该方法简单,对地震数据压缩行之有效。     

基于复值曲波变换的地震数据重建方法

传统的地震数据采样必须遵循奈奎斯特采样定理,而野外数据采样可能由于地震道缺失或者勘探成本限制,不一定满足采样定理的要求,因此需要进行叠前数据重建,以满足后续处理的要求。为此,笔者提出一种基于复值曲波变换和凸集投影算法(POCS)的地震数据重建方法。首先引入凸集投影算法和能够刻画地震数据局部化特征的复值曲波变换,在重建过程中针对传统阈值参数收敛速度较慢的缺点,提出了指数平方根衰减规律的阈值参数,并采用硬阈值算子对二维地震数据进行重建,从而降低了迭代次数,提高了重建精度。通过理论数据的研究表明,该方法重建效果显著,最后将该技术应用于实际地震勘探资料,获得较好的应用效果。     

拉东变换压制多次波在新疆巴楚地区的应用

多次波是地震勘探中一种常见的干扰波,巴楚地区浅层强反射界面导致了多次波异常发育,严重影响了地震数据的成像,甚至产生了构造假象,压制多次波对于新疆巴楚地区地震数据的构造精确成像具有重要的理论意义与实用价值。结合巴楚地区地质结构特征,对地震采集的原始数据进行深入剖析,同时根据多次波自身的特点,在试验了多种技术方法之后,认为拉东变换压制多次波方法能够在该区取得地震数据构造成像满意的效果。给出了在不同数据域准确识别多次波的几种基本方法和判断准则,同时详细阐述了拉东变换压制多次波的基本原理。通过工区实际资料的处理,系统地演示了该技术方法运用的过程、结果及关键点。     

Shearlet域稀疏约束地震数据重建

在地震数据处理流程中,通常对不规则的、稀疏的或者缺失的地震数据进行插值处理,通过插值方法来避免多次波的预测错误和成像假频等现象,使地震数据处理更加精准。Shearlet变换是一种多尺度变换,具有最佳的稀疏性、方向性以及局部化特性。将Shearlet变换与基于Landweber加速下降迭代方法结合起来对地震数据进行插值,在保证求解精度的同时提高了计算效率。信号和噪声在Shearlet域具有不同的分布特点,通过阈值法压制随机噪声,可提高算法的抗噪性。此外,采用jitter采样的方式,更好地压制了假频信息。理论和实际地震数据验证了该方法的有效性。     

小波变换与地震奇异性属性

地震属性是由地震数据提取的有关地震波的几何形态、运动学特征、动力学特征和统计学特征等组成。这里主要讨论了一种新的地震属性(奇异性属性)计算方法。由于地震数据携带着大量的奇异性信息,因此,这一新的地震属性是在偏移地震数据的小波变换和奇异性分析的基础上形成的。小波变换使检测数据中的局部奇异程度成为可能。为能够估计出偏移地震数据每一采样点的局部奇异规律,从小波变换系数的角度来提取地震信号的奇异指数。这种新属性实行单道处理,不需要子波和速度信息,它给出了地下分层情况和断层位置。在地震数据奇异性属性的基础上,用小波变换来划分地层旋回,从而提高了小波变换在地震解释中的应用。     

应用时空变—倾角扫描叠加KL变换提高地震资料信噪比

使用KL变换可去除地震资料中随机噪声和相干噪声，提高地震资料信噪比。但常规KL变换仅能加强水平方向同相轴，对倾斜或弯曲同相轴处理效果较差，且在处理大量数据时计算成本很高，在实际生产中难以广泛应用。一些文献中针对常规KL变换的两个缺陷进行了改进，即使用倾角扫描叠加KL变换和采用数据分块技术。作者在本文中将这两个方面的改进有机地结合起来，提出应用时空变－倾角扫描叠加KL变换对地震资料进行处理，并指出使用该方法时需要注意的一些问题。理论模型和实际资料处理效果表明，使用本方法可以克服常规KL变换的缺陷，有效改进剖面处理效果。     

Radon变换的MATLAB实现

Ｒａｄｏｎ变换（包括线性、抛物Ｒａｄｏｎ变换）是地震数据处理中的一种强有力工具。作者在文中阐述了其原理及最小平方算法,并给出了ＭＡＴＬＡＢ语言缩写的源码,试算结果表明了该算法的有效性和程序的正确性,该程序可直接当作工具来使用。     

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

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

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

提高Radon变换分辨率的研究一直是地震数据处理方法研究的热点之一。常用的提高分辨率的算法是在频率域进行的,然而此方法在频域中计算的模型权重是偶然耦合的,会对所有同相轴施加相同的权重,导致高能量同相轴地震道生成假像。此文给出了3种时间域提高Radon变换分辨率的方法,IST(iterative shrinkage thresholding)、FIST(fast iterative shrinkage thresholding)和SRTIS(sparse radon transform iterative shrinkage),并对3种算法的收敛效率和计算效果进行了对比分析。理论和实际数据的测试表明,SRTIS方法在效率和效果的对比上均优于其他两种,并且具有较好的去除多次波能力。     

双曲Radon变换参数分析

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

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

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

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

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

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

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

地震波初至拾取的分形研究

在地震数据处理过程中 ,传统的方法都是以振幅、频率、相位等的变化为依据来拾取地震初至波。作者突破了这些传统依据 ,根据分形的原理 ,用时窗的变化来反映地震道曲线的局部形态特征 ,通过一元线性回归模型统计出曲线的形态特征参数———分数维 ,然后由分数维的变化来判断初至地震波 ,得出了一种高效的拾取初至地震波的新方法。针对该方法的计算过程开发了处理软件模块。理论的模拟证明该方法非常适用于目前越来越庞大而复杂的地震勘探数据资料的处理。通过实际的试算表明软件模块的适用性强。该方法精度高 ,计算速度极快 ,抗噪能力强 ,在低信噪比的情况下也能快速精确地拾取初至地震波。     

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

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

峨眉山柿子坪热田地表氡气异常特征

本文通过柿子坪地区地表氡气测量,对热田地表氡气异常特征作了初步分析。该区地表氡气高值异常呈环状产出,并严格受构造控制。同时指出地表氡气测量是热田勘查的有效途径。     

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

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