多道自适应匹配滤波方法压制表面多次波

表面多次波压制是海洋地震数据处理的关键环节,基于波动方程的多次波压制方法可由多次波模型预测和自适应相减两步实现.因此,除了采用有效的算法实现高精度的多次波预测外,多次波的自适应相减也是改善多次波压制效果的重要手段.文中采用多道自适应滤波器完成表面多次波的自适应相减,利用数据驱动的基于波动方程预测多次波模型的高频重建道、Hilbert变换道、以及它们相应的上下延拓道改善预测多次波模型的振幅、相位、旅行时以及频带信息,使得修正后多次波模型与实际地震数据中的多次波更好地匹配.实际数据测试表明,文中所述方法可有效地完成表面多次波的自适应相减.通过讨论分析交叠时窗和滤波器的长度对自适应相减效果的影响,表明合理地利用时窗以及选择滤波器长度将有效地改善自适应相减的效果.     

反馈迭代法压制表面多次波效果分析

地震数据处理通常假定反射数据仅由一次波组成,地震勘探资料中的表面多次波通常被视为相干噪音予以压制.基于波动方程预测的反馈迭代法为数据驱动的方法,可有效压制复杂地下介质的表面多次波,通常可由级数展开法和迭代法两种方法实现.文中分别从理论阐述以及应用效果分析方面对基于波动方程预测的反馈迭代法的两种实现途径进行研究,并采用GPU/CPU协同并行加速计算预测表面多次波,改善以往CPU计算效率低的状况.文中对含表面多次波的炮数据和复杂的SMAART模型进行了计算,并将利用级数展开法与迭代法压制表面多次波的效果进行对比分析,结果表明,迭代形式反馈迭代法压制表面多次波的方法效果更佳.由于实现迭代法和级数展开法需要完成的自适应匹配相减的次数不同,迭代法的计算成本略高于级数展开法.     

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

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

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

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

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

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

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

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

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

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

基于多道卷积信号盲分离的多次波自适应相减方法

本文将多次波自适应相减问题表示为一个多道卷积信号的盲分离问题.利用2D卷积核来表示预测多次波和实际多次波之间的差异,并采用分离出的一次波信号的非高斯性最大化作为优化目标,我们提出一种基于多道卷积信号盲分离的多次波自适应相减算法.为了求解上述非线性优化问题,所提方法将其转化为一个迭代线性优化问题,采用迭代最小二乘方法加以实现.由于采用了多道卷积信号盲分离模型,所提方法能够适应预测和真实多次波之间在时间及空间上的变化.通过对简单模型数据、Pluto数据和实际数据进行处理,验证了所提算法的有效性.     

表面多次波压制的研究进展:回顾与展望

在绝大多数的地震资料处理中,将多次波视为相干噪音进行压制.表面多次波的压制方法主要分为两大类:基于波动方程预测的预测减去法和基于信号处理的滤波法.本文系统地回顾了两类多次波压制方法,分析了各种方法处理各类介质的优势及存在问题.本文也阐述了基于波动方程预测的自适应相减方法的发展.分析了多次波压制方法的发展趋势,随着GPU...     

泰勒展开近似 带连续性约束的L1范数方法用于多次波自适应相减(英文)

一次波L1范数最小化的多次波自适应相减方法,简称L1方法,是基于匹配滤波器设计的多次波自适应相减算法中的一种常用方法.当一次波和多次波混杂在一起时,L1方法有时会伤害一次波,导致一次波同相轴的连续性变差.本文利用预测误差滤波器度量一次波同相轴的连续性,在L1方法的皋础上,提出一种能够在压制多次波的同时,尽量保持一次波同相轴连续性的多次波自适戍相减算法,简称连续性约束L1方法.利用Pluto模犁数据进行多次波相减的结果表明,连续性约束L1方法能够在有效压制多次波的同时,更好地保护一次波.     

基于表层多次波数据的近道地震数据插值方法研究(英文)

本文提出基于原始含表层多次波数据实现叠前共炮集地震数据插值。相对于利用相邻道的信息变换或外推插值用于缺失的地震数据重建,本文方法利用表层多次波数据互相关构建准一次波,将蕴含在表层多次波数据中的,而在采集记录中表现为缺失的近炮检距信息提取出来,并在滑动时间空间窗内采用最小二乘匹配滤波和均方根振幅校正方法进行准一次波校正而后用于数据插值重建。本文方法适用于表层多次波比较发育,同时又存在数据缺失尤其是近炮检距数据缺失情况。方法易于实现,不需多次波和一次波的提取,利用多次波中蕴含的信息实现缺失的地震数据弥补,为含有表层多次波的数据进行近炮检距地震信息的插值重建提供了一个很好的思路。     

Estimating primaries by sparse inversion of the 3D curvelet transform and the L1-norm constraint

In this paper, we built upon the estimating primaries by sparse inversion (EPSI) method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and L1-norm regularization, and use alternating optimization to directly estimate the primary reflection coefficients and source wavelet. The 3D curvelet transform is used as a sparseness constraint when inverting the primary reflection coefficients, which results in avoiding the prediction subtraction process in the surface-related multiples elimination (SRME) method. The proposed method not only reduces the damage to the effective waves but also improves the elimination of multiples. It is also a wave equationbased method for elimination of surface multiple reflections, which effectively removes surface multiples under complex submarine conditions.     

Suppressing multiples using an adaptive multichannel filter based on L1-norm

Adaptive subtraction is an important link for removing surface-related multiples in the wave equation-based method. In this paper, we propose an adaptive multichannel subtraction method based on the L1-norm. We achieve enhanced compensation for the mismatch between the input seismogram and the predicted multiples in terms of the amplitude, phase, frequency band, and travel time. Unlike the conventional L2-norm, the proposed method does not rely on the assumption that the primary and the multiples are orthogonal, and also takes advantage of the fact that the L1-norm is more robust when dealing with outliers. In addition, we propose a frequency band extension via modulation to reconstruct the high frequencies to compensate for the frequency misalignment. We present a parallel computing scheme to accelerate the subtraction algorithm on graphic processing units (GPUs), which significantly reduces the computational cost. The synthetic and field seismic data tests show that the proposed method effectively suppresses the multiples.     

不同阶次自由表面相关多次波预测与成像方法

海水与空气间的强波阻抗界面使得海洋地震数据普遍发育自由表面相关多次波,多次波信息的利用是提高海洋地震资料成像品质的新突破点.近年来发展了一系列多次波成像方法,干涉假象是制约其应用推广的关键问题之一.为了避免假象影响,本文提出了不同阶次自由表面相关多次波预测与成像方法,首先,修改了传统SRME(表面相关多次波衰减)方法中的边界条件,通过多次波升阶次与匹配相减的方法预测出不同阶次自由表面相关多次波;其次,基于单程波偏移算子和"面炮"偏移策略,以一次反射波或第(N-1)阶自由表面相关多次波为下行波场正向延拓,以第1阶多次或第N阶多次波为上行波场逆向延拓,并在每一层互相关成像得到第1阶或N阶多次波单独成像.本方法避免了低阶多次波和高阶多次波产生的相关假象,且相对于全波算子的偏移方法具有较高的计算效率,增强了多次波成像方法的实用性.单层模型和三层模型测试验证了本方法的正确性,并在我国某深海探区实际资料处理中得到了成功应用.相对于传统一次波成像,分阶次多次波成像具有更高的照明均衡度、垂向分辨率和信噪比.本研究表明,海洋多次波成像是一次波成像的有力补充,对于稳定海底沉积的深海地区,具有一定的应用前景.     

基于Marchenko理论压制自由表面多次波方法研究

常规虚源点Marchenko自聚焦多次波预测方法只适用于预测不含自由表面的多次波模型,局限于压制层间多次波,该方法在构建上下行格林函数场前,必须从反射响应中去除所有与表面相关的多次波.本文对构建上下行Marchenko格林函数方程进行改进,得到了包含一次波、层间多次波和自由表面多次波的格林函数,利用改进的Marchenko自聚焦预测方法预测自由表面多次波.本文利用水平层状模型数据及SMARRT模型数据证明,改进后的Marchenko法预测海底相关的自由表面多次波效果较为理想,该方法避免了常规SRME自由表面多次波预测方法需要近道重构的缺陷,能够有效提高地震资料的信噪比和分辨率.     

基于相关迭代的非因果匹配滤波器多次波压制方法

针对常规多次波压制(SRME)方法存在的缺陷进行了改进.首先在多次波模型道预测阶段,采用数据相关和迭代更新的策略预测多次波模型道,该方法降低了SRME方法采用数据空间褶积对输入数据要求严格的限制,提高了对近偏移距缺失和空间假频数据的适应性.在匹配相减阶段,本文设计了一种非因果非平稳的匹配滤波器,该滤波器可以对整道多次波模型进行处理,而且当预测多次波模型道滞后于实际多次波时,也能够对多次波模型道进行很好的匹配.模型数据和实际数据试算证明该方法在多次波模型道预测和匹配相减阶段的优越性.     

逆子波域消除多次波方法研究

SRMA（与表面相关多次波的衰减）算法包含预测和相减两步.相减算法中,当多次波与反射波同相轴相交时,如何有效减去多次波、保留反射波,是面临的主要问题.通过分析非正交性对滤波器（逆子波）的影响,可以证明：逆子波因非正交性产生的误差呈近似的高斯分布.在此基础上,本文提出了在逆子波域（单道自适应相减滤波的滤波算子的集合）,利用其误差的概率分布特征,对逆子波进行估计,用逆子波的估计对逆子波进行校正来消除多次波的方法.其步骤为：首先用SRMA方法预测出表面多次波,并对每一单炮进行单道自适应相减,得到逆子波,形成逆子波域;其次,在逆子波域采用中值滤波,提取接近真实逆子波的逆子波估计;第三,在逆子波域用逆子波估计对畸变的逆子波进行校正;最后采用校正后的逆子波来衰减多次波.通过简单模型和SMARRT模型的测试,该方法不仅能够有效减去多次波,而且在相交的区域,能够保持反射波同相轴的连续性并恢复其正确的振幅.