基于波动方程预测和双曲Radon变换联合压制表面多次波

基于波动方程预测的表面多次波压制方法可处理复杂地下介质的地震资料,但计算成本较高.基于滤波的多次波压制方法计算效率较高,但其成功应用仅局限于一次波和多次波有明显时差差别的地震数据,对来自速度逆转等复杂介质数据则较难获得满意的压制效果.本文将波动方程预测的反馈迭代法和滤波法有效结合,采用GPU(图形处理器)和CPU协同并行加速计算粗略预测表面多次波,随后在双曲Radon域比较分析原始数据和预测的多次波,设计合理有效的Butterworth型自适应滤波器,滤出原始数据Radon域中的多次波能量,进行Radon反变换后,在时空域将多次波从原始数据中减去,得多次波压制结果.文中对理论模拟的单炮数据、复杂的SMAART模型以及实际地震数据进行了计算,结果表明,结合基于波动方程预测和双曲Radon变换的方法有效突破了两种方法各自的局限性,可高效高精度地压制复杂地下介质的表面多次波.     

海上多次波压制与成像方法研究进展

海洋地震勘探中多次波非常发育,如何合理处理多次波是一个很重要的课题.本文首先对海上多次波的产生机理、分类及识别方式进行了回顾,然后对目前国内外发展的多次波处理技术进行了总结.文中从将多次波作为噪音和将多次波视为有效信息,即多次波压制与成像这两个方面分别进行了总结.多次波在一次波成像中的严重干扰推动了多次波压制技术的不断发展,主要发展了基于滤波理论和基于波动理论的两大类多次波压制方法,但是综合考虑计算效率和多次波压制精度,近年来,将两种方法相结合成为一种新的发展趋势;而多次波作为来自地下界面的真实反射也包含了丰富的地下信息,目前,诸多地球物理学家分别从对多次波进行单独成像和对一次波和多次波联合成像的角度对多次波的有效利用开展了大量研究,结果显示利用多次波可以获取比一次波更多的地下信息,但是该技术发展还不够成熟,这将需要广大地球物理工作者继续开展深入的研究.     

多次波压制的研究现状与进展

目前地震资料处理技术主要利用一次波对地下地质体进行成像或反演,多次波的存在会严重影响地震成像、反演与解释结果.因此,多次波通常作为噪声在地震数据叠前处理阶段进行压制.为了合理选择不同的多次波压制方法,提高地震资料信噪比和保真度,本文系统综述了多次波压制的国内外研究现状,简要介绍了滤波法、预测相减法、稀疏反演法的基本原理,并围绕理论假设条件与实际数据要求,分析比较了各种方法的适用性与优缺点.基于波动方程的预测相减法克服了滤波法关于地震信号统计特征的假设,增强了对复杂地震数据的适应性,但最小二乘自适应相减会损害有效信号.稀疏反演法采用全波形反演的方法估计一次波与多次波,避免了自适应相减步骤,能够更加保真地恢复有效信号.本文通过实际算例重点阐述了稀疏反演法相对于反馈迭代法在多次波压制中的突破性进展,同时指明了稀疏反演法存在的问题并展望了未来的发展方向.     

基于稀疏反演三维表面多次波压制方法

三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.     

分层层间多次波预测与自适应压制方法研究

在地震勘探中,层间多次波在地震剖面上形成的虚假界面会严重影响对地质构造的解释.由于层间多次波在走时、频率和形态上与一次波的差异较小,因此较难对层间多次波进行识别和压制.本文使用分层层间多次波预测方法与基于2D卷积信号盲分离的多次波自适应相减方法相结合的策略对层间多次波进行压制处理.首先推导并阐述了分层层间多次波预测方法的理论基础,然后对地下介质中存在的多个多次波产生界面实现了层间多次波的分层预测,最后采用基于2D卷积信号盲分离的多次波自适应相减方法对层间多次波进行压制,以更好的保护有效信号.分层预测的方法通过分离层间多次波产生层位,对目标层位产生的所有层间多次波直接进行预测,该方法避免了共聚焦点方法构建共聚焦点道集的聚焦运算,降低了计算成本,同时该方法是基于数据驱动的,算法容易实现.模型数据和实际地震资料试算结果表明,本文方法能直接利用地表数据对层间多次波进行压制,并有效的保护一次波．     

抛物线Radon变换方法在压制地震多次波中的应用

多次波的存在会降低地震资料的信噪比,影响地震资料处理效果和后续的地质解释精度,压制多次波干扰是地震资料处理中的一个重要环节。利用有限差分方法正演模拟几个典型地质模型的地震波场,之后进行动校正,再通过抛物线Radon变换方法把t-x域一次反射波和多次波变换到τ-q域,由于存在速度的差异,Radon域记录中的一次反射波和多次波的能量互相分开,在Radon域对多次波的能量切除,保留反射波的有效信息,达到压制多次波的效果,提高资料的信噪比。通过对实际地震数据进行抛物线Radon变换结果也进一步验证了该方法在压制多次波中的应用效果。      

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

波路径偏移压制层间多次波的理论与应用

消除层间多次波是地震勘探资料处理研究领域的难题,尤其对于实际资料的处理,到目前为止还很难找到一种完全有效的方法. 本文给出了仅对一次波成像既波路径偏移方法压制层间多次波方法,在共炮道集和共检波点道集分别计算炮点射线的入射角和检波点射线的出射角,由此计算的角度作为射线追踪的初始角度,计算地震波射线的传播路径. 结合由程函方程计算的走时表,判断偏移范围是反射波还是多次波. 在前期偏移过程压制多次波的理论研究基础上,本文主要研究波路径偏移消除多次波的应用部分. 为了进一步说明效果的有效性,计算了在单炮和共成像点道集压制层间多次波,给出了实际资料的压制多次波的偏移结果.     

基于波动方程三维表面多次波预测方法研究

与传统的二维表面多次波预测算法相比,基于波动方程的全三维表面多次波预测方法无需对地下介质做简单近似,其更符合地震波在地下介质中传播的真实状况,是地震资料处理中解决多次波预测问题的强有力工具.本文从三维多次波预测的基本理论出发,给出了全三维多次波预测算法的预测矩阵表示、计算方法以及实现条件,采用GPU(图形处理器)加速全三维表面多次波预测,较传统的CPU串行计算,GPU并行预测表面多次波的计算效率约提高165倍.文中分别利用二维和三维表面多次波预测算法对理论模拟的含表面多次波的三维地震数据进行多次波预测计算,对比分析结果表明,相比于二维算法,文中所述的基于波动方程的全三维表面多次波预测效果明显改善,其计算精度更高,辅以合理有效的自适应相减算法,可获得高精度的地震勘探资料表面多次波压制数据.     

基于深层神经网络压制多次波

有效压制多次波一直是地震勘探中的难点问题.尽管已发展了多种多次波压制方法,但仍存在多次波压制不全、计算耗时长等缺陷,使得应对复杂地质地震数据多次波压制具有挑战性.为了突破现有多次波压制方法的局限性,本文提出了一种基于深层神经网络的多次波压制方法,采用的深层神经网络是一种改进的具有卷积编码器和卷积解码器的U-net网络.不同于常规方法依赖于滤波或波动理论,该方法仅依赖于大量训练数据.训练数据以含多次波的原始地震数据作为输入,不含多次波的地震数据作为输出,通过最小化损失函数来优化神经网络参数.训练成功的网络模型具备较好地分离多次波和一次波的能力,可直接用来快速压制地震数据中的多次波,避免了常规方法涉及的大规模计算.工业界模型数据测试结果表明,本文提出的深层神经网络方法能有效压制复杂地质地震数据中的多次波,同时还具有较高的泛化能力和多次波压制效率.