编码与解码框架下的局部平面波域浅水多次波压制方法

多次波压制是提高地震数据品质与成像结果分辨率的重要手段之一,浅海环境给多次波压制引入了额外的难点。浅海环境中主要的多次波是水体相关多次波,由于海底-海面路径短,且海面、海底都是强反射界面,因此该类多次波周期短、能量强,对有效波的干扰严重,采用常规多次波压制方法不能将其很好地压制。提出了针对浅水多次波的有效压制方法:将多次波的产生与压制用线性编码-解码系统进行描述,即将上行波在水层震荡产生多次波的过程视为编码过程,将剥离多次波的过程视为解码过程,并从线性系统求解角度给出了平面波域多次波压制算法。首先将水体相关多次波的预测转移到局部平面波域,利用稀疏τ-p变换减弱了采样不规则、近偏移距缺失造成的影响;然后利用射线理论描述水体介质中的地震波传播,高效生成多次波模型;最后在τ-p域进行匹配相减,并利用一次波与多次波的视速度差异,提高匹配精度。利用浅海OBC数据和拖缆数据验证了方法的有效性和实用性,通过对比压制浅水多次波前后的地震剖面和自相关结果,展示了方法压制浅水多次波的效果和应用前景。     

利用格林函数的多次迭代鬼波压制方法

在海洋宽频地震处理技术中,鬼波压制是拓宽地震资料频带、提高地震资料分辨率的关键,常用方法是在预测鬼波的基础上,将其从全波场中消除。显然精确地预测鬼波是鬼波压制技术的难点。依据格林函数理论,提出了基于鬼波传播的射线路径预测鬼波波场的多次迭代技术及流程,结合研发的曲波域匹配相减技术,有效压制了鬼波。由于预测鬼波波场较为复杂,导致在匹配相减时会出现损伤一次波或鬼波未彻底压制等情形。为此,进一步引入多次迭代思想,将匹配相减后得到的一次波的波场作为输入,再次通过格林函数预测鬼波并相减,得到更干净的一次波场;历经多次迭代可逐步提高鬼波预测的精度。模型数据和实际资料的处理结果表明,该方法在压制鬼波的同时能更好地保护一次波。     

与界面有关的层间多次波预测及压制

现今有关层间多次波的压制方法大多针对强反射界面向下反射的能量形成的层间多次波。本文通过理论分析和模型验证,扩展了与界面有关层间多次波的概念,提出了强反射界面向上反射的能量形成的层间多次波,即LLRIM多次波(Lower layer related internal multiple)。除上述概念以外,本文根据层间多次波压制理论,解剖了与界面有关的层间多次波预测算法,给出了理论推导过程以及时域和频域表达式。通过分析和实验验证,界面向上反射形成的层间多次波的预测可以由界面向下反射形成的层间多次波的预测方法类推得到,本文给出了预测表达式和模拟实验结果,证明了方法的有效性。     

海上单道地震资料中多次波的衰减

海洋单道地震资料中存在严重多次波,但针对单道地震资料识别和压制多次波的方法不多。通过剖面观察、自相关分析等手段将海洋地质调查单道地震资料中存在的多次波分为两大类。针对各自的特征,利用预测反褶积和模型拟合方法对它们进行了衰减试验。试验表明,预测反褶积用来衰减周期性强的短周期多次波比较有效;对于不规则的崎岖构造产生的短周期多次波,预测反褶积会有二次残留,可以使用自适应预测反褶积进行压制,但会造成数据空白。模型拟合方法适用于压制与一次波有明显区别的长周期多次波,其中参数L_(NC)值对多次波衰减的效果有影响,值太大则难以识别多次波;值太小,会使有效波被衰减。对于多次波与一次反射差别明显(反向斜交)的情况,L_(NC)值可以取小一些使得多次波衰减更干净;对于多次波与一次反射差别不明显(同向斜交)的情况,L_(NC)值要取大一些,以免造成有效反射被误认为多次波而衰减掉。     

倾斜界面的多次波在二维平面波域中的选择性衰减方法

在许多地质条件下,气-水界面处的强反射为地震记录中的多次波提供了大部分的能量,在此,我们阐述一种自由界面的多次波衰减方法,就是运用众所周知的恒定内插技术推导出一个地震记录的反射算子模型,再利用它进行反射波衰减。当然这种方法是在二维平面波域中运行的;而且针对不同的射线参数,考虑到地质构造的横向变化。在这种横向变化是平滑的条件下,地震数据在二维平面波域中压缩得很好,而且最终的有限带宽矩阵也大大减少了计算时间。这种处理方法的一个重要特点就是它的灵活性,也就是允许从任意选择的反射中剔除多次波。 为了得到没有多次波的资料,波动理论基础上的多次波衰减方法就是试图估算出源函数或者地下界面反射系数。而我们的方法就是综合利用了上述两种方法。首先利用反射系数预测多次波旅行时间,然后寻求一个源函数预测振幅。模拟和实测数据都表明,这种方法用来衰减界面多次波时是稳定而成功的。     

倾斜界面的多次波在二维平面波域中的选择性衰减方法

在许多地质条件下,气－水界面处的强反射为地震记录中的多次波提供了大部分的能量,在此,我们阐述一种自由界面的多次波衰减方法,就是运用众所周知的恒定内插技术推导出一个地震记录的反射算子模型,再利用它进行反射波衰减。当然这种方法是在二维平面波域中运行的,而且针对不同的射线参数,考虑到地质构造的横向变化。在这种横向变化是平滑的条件下,地震数据在二维平面波域中压缩得很好,而且最终的有限带宽矩阵也大大减少了计算时间,这种处理方法的一个重要特点就是它的灵活性,也就是允许从任意选择的反射中剔除多次波。为了得到没有多次波的资料,波动理论基础上的多次波衰减方法就是试图估算出源函数或者地下界面反射系数。而我们的方法就是综合利用了上述两种方法。首先利用反射系数预测多次波旅行时间,然后寻求一个源函数预测振幅。模拟和实测数据都表明,这种方法用来衰减界面多次波时是稳定而成功的。     

海上地震资料处理中的组合压制多次波技术

本文针对澳大利亚海上地震资料中全程多次波和层间多次波发育的特点,分析了各类多次波压制方法的优缺点,提出了针对该区海上地震资料多次波压制策略:①在第一次速度分析后应用Radon变换压制大部分全程多次波;②应用预测反褶积压制大部分层间多次波;③在第二次速度分析后,应用高精度Radon变换压制剩余多次波;④在叠前偏移CRP道集上应用Radon变换和内切除压制小时差层间多次波。应用上述组合方式可以有效地压制多次波,突出一次有效反射波。     

多种多次波压制技术提高可燃冰地震识别

厄瓜多尔西海岸大陆斜坡地区发育大面积的可燃冰,其地震反射特征与海底反射极其相似。该区属 于典型的深水区域且海底崎岖,其海洋地震资料存在自由表面多次波、层间多次波和绕射多次波等多种 类型波。强能量多次波的存在,降低了资料的信噪比,影响了一次反射波的准确成像,也严重影响了可燃 冰的地震识别。采用自由表面多次波压制（SRME）技术、高分辨率拉冬变换多次波压制技术和绕射多次波 压制技术等串联组合,解决了该区多种类型多次波的压制问题,提高了利用地震资料识别可燃冰的能力, 并预测了该区可燃冰的分布范围。     

浅水多次波衰减

浅水海域地震资料通常会有较为严重的鸣震等短周期多次波干扰,一般采用预测反褶积进行压制,但是在地质构造较复杂地区其效果并不理想。而表层相关多次波压制方法(SRME)(Surface related multiple elimination)方法,由于海底反射信息较差,很难压制浅水海底相关等短周期多次波。本文使用数据驱动的浅水多次波压制方法,该方法同时基于SRME褶积方式与多道预测算子,在浅水复杂构造区压制短周期多次波方面可以取得比预测反褶积更好的效果。在实际资料应用中,使用本文方法与预测反褶积组合的方式,其多次波压制效果更为理想,优于单独使用其中任何一种方法。     

基于贝叶斯理论的Shearlet域一次波与多次波分离

多次波的去除效果很大程度上决定了地震成像的质量,尽管表面相关多次波的压制(SRME)研究方法取得了很大进步,但预测出来的多次波在相位和振幅上的混乱仍然是一个难题,因此减去的效果好坏尤为关键。常规的最小平方匹配方法并不能很好地处理这种预测的误差,利用Shearlet变换的多尺度、多方向特性,可以使预测的信号在Shearlet域变得稀疏、圆滑。以此为前提,结合贝叶斯概率最大化理论,建立了一个有效的一次波与多次波分离算法。以SRME预测的多次波为例,从贝叶斯预测角度出发,通过精确求解最优化问题就可以达到一次波与多次波分离的目的,该算法能够更好地控制预测信号和实际信号的误差。理论和实际数据的实验表明,相比最小平方匹配减去方法,该方法可以有效地提高减去效果,更好地压制多次波以及高频混叠,并且可以更好地估计一次波。     

τ-p域水体模型驱动压制浅水区水层多次波

在近海油气地震勘探中,压制在海面及海底之间的短周期水层多次波一直是海洋地震资料处理的难点,传统上采用预测反褶积或SRME压制此类多次波存在很大局限性。因此,本文提出了一种τ-p域水体模型驱动压制浅水区水层多次波方法,该方法在τ-p域用已知水层参数(水速和水深)构建初始模型,并利用波场延拓方法来预测水层多次波模型,然后采用自适应相减从原始数据消除多次波。通过理论模型试算和实际数据应用,表明此法能有效压制浅水区水层多次波。     

塔里木盆地孔雀河地区多次波压制方法与效果

针对孔雀河多次波的特点,应用稀疏PRT压制多次波方法,利用多次反射波的特征及其与一次反射波的差异,加强对数据的抛物线分离,更精确、更有效地实现了多次波的分离与压制,达到压制孔雀河地震资料目的层强多次波干扰,恢复弱有效反射的目的。     

基于数据一致性原理预测与压制自由表面多次波的效果分析

 本文基于数据一致性预测与压制自由表面多次波方法原理,研制了迭代自适应压制多次波方法,用于模型数据和海洋实际地震数据的多次波压制处理。模型试算和实际海洋地震数据处理结果表明,与现有的常规多次波消除方法相比,此方法收敛速度快,并能很好地保持有效波能量。在应用此法时,要特别注意以下几点：①基于数据一致性原理,通过时间—空间域褶积预测自由表面多次波,多次波预测算子可以是含多次波的原始数据,也可利用已部分消除多次波的数据（如抛物线拉东变换滤波方法压制多次波后的数据）,这样能更快地获得较好的多次波压制效果;②在一定简化条件下,预测与压制自由表面多次波可表达为一个线性最小平方反演的过程,实现过程是一个迭代自适应过程,只需要很少的几次迭代就能获得较好的多次波压制效果;③对于地下地层结构非常复杂时的多次波压制,文中方法在更有效地压制各种自由表面相关多次波的同时,还可以压制部分层间多次波,同时也能更好地保持一次波有效能量。     

消除多次反射的一种 新方法——“消去”法

一般地说长排列和多次复盖能有效地压制多次波,然而,当多次反射的能量比一次反射的能量大得多时,在最后的剖面上就出现大量多次波剩余干扰能量,以致难于分辨出真反射波,如图1的左部所示。如果野外资料用相应于多次反射的动校正函数作了正常时差校正,就能获得一张非常连续的剖面。这     

Radon变换压制层间多次波技术在高石梯-磨溪地区的应用

陆上地震资料中,层间多次波与一次波的速度差异较小,一次波的速度拾取困难,利用Radon变换压制多次波往往难以取得好的应用效果。在四川盆地高石梯-磨溪地区的地震资料中,震旦系-寒武系存在多次波的干扰,致使有效反射的振幅畸变,弱的有效反射被掩盖,造成测井合成地震记录与实际地震资料匹配差。该多次波具有干扰能量强、与有效波速度差异小的特点。基于处理和解释一体化的研究思路,在速度谱上利用时间域的构造层位来约束一次波速度的拾取过程,以减小多次波干扰对速度拾取的影响,提高平面上和纵向上速度拾取的精度;利用不同权系数的叠加速度进行扫描处理和Radon变换压制多次波,确定最佳匹配权系数,实现更高精度的多次波压制。在高石梯-磨溪地区,三维地震资料经压制多次波处理后,地震剖面上多次波的能量得到有效衰减,波组特征更为清楚,测井合成地震记录和井旁地震道相关性大幅提高。应用效果表明,这一技术在压制多次波的同时,能够有效地保护一次波能量,减少压制多次波的不确定性,提高地震资料的信噪比和成像效果,是压制陆上层间多次波的有效技术。     

电磁场延拓成像中对多次波和透过波的压制方法

与地震勘探一样，大地电磁勘探也存在着多次波与透过波干扰，研究这一问题对于大地电磁勘探资料的地质解释非常重要。本文从理论上分析了大地电磁拟地震解释中多次波的存在，指出了大地电磁测控电磁波反射函数中不存在多次波干扰而仅有上透过波干扰，导出了电磁波反射函数与上透过波干扰的数学关系，指出了压制上透过波干扰的选频办法，最后，用大地电磁拟地震解释法和电磁波反射函数延拓成像法对理论模型进行了计算，证实了前面的理论。     

地震反射数据的全波动方程向下延拓

本文提出一种压制地震记录中多次反射的新方法。它是建立在应用全声波方程(双曲线形式)作为向下延拓算子的向下延拓方法基础上的。我们证明记录波场的向下延拓可绘成反射率函数曲线,而没有多次反射波。该方法成功地应用于单道平面波以分解(倾斜叠加)合成记录和野外数据。     

反射地震资料全波方程的向下延拓

研制出一种压制地震记录中的多次反射的新方法。这种新方法建立在向下延拓步骤上，这种向下延拓步骤使用全波方程式（双曲形）作为向下延拓因子。我们将演示记录波场的向下延拓测绘出多次反射同相轴的射率函数的情景。本方法被成功地应用于平面波分解的（倾斜迭加）合成资料及野外资料的合个道中。     

应用逆散射级数波场预测和2D卷积盲分离压制层间多次波

多次波压制是地震数据处理的重要环节。基于波动方程的预测相减法是压制层间多次波的常用方法,它包括层间多次波预测以及自适应相减两个步骤。文中提出了一种改进的层间多次波压制方法,通过引入阶跃函数、有限积分区间等假设条件简化逆散射级数公式,实现层间多次波的预测;然后引入2D卷积盲分离方法实现层间多次波的自适应匹配相减。理论模型试算及实际资料应用效果表明,该方法能在不依赖速度模型前提下高效实现层间多次波预测;采用的多道卷积盲分离自适应相减方法,比传统L₂范数最小化自适应相减方法和单道卷积盲分离自适应相减方法,能在有效压制层间多次波的同时,更好地保护一次波。     

用最佳加权改进迭加质量

理想的地震剖面只存在一次反射波,即震波在界而上只反射一次。然而,实际上的记录剖面或多或少会被多次波所扰乱,即被反射的波往往多于一次。最糟糕的是这些多次波可以在界面上与一次波同时到达并具有足够的振幅与一次波相混淆。简单的共反射点迭加通常是足以压制多次波的,但在有些地区对多次波作些补充压制也是必要的。在迭加前对资料作适当的加权可以使结果进一步改善。本文叙述了一种确定最佳加权的方法,并展示了处理的实际结果。最初的权决定于被迭加的一次波和多次波模型。然后逐步变更各道的权值,直至找