6 地震处理之多次波衰减

合集下载

地震勘探中的常见地震干扰波及压制方法

地震勘探中的常见地震干扰波及压制方法论文提要在地震勘探中激发地震波时，由于激发、接收条件，自然环境和地表条件的影响，我们所采集到的地震数据中，既有有效波也有干扰波。

根据干扰波的物理特征、形成机理和形态，常把地震数据上的噪声分为规则噪声和随机噪声两大类。

规则噪声具有明显的运动学特征 ,如：面波、线性干扰、平行折射、声波、多次波干扰等，可以根据其运动学特征选择针对性的衰减方法；随机噪声是一种无规律的噪音，如：自然界风吹草动所产生的猝发脉冲、野值等。

为了提高地震勘探的精度，完成在各种复杂地区的勘探任务，使地震资料能更真实地反映地下的地质情况，如何突出有效波，压制干扰波就成为一个极其重要的问题。

通过暑假的实践，本论文中针对地震勘探中的常见地震干扰波进行总结、分类、衰减，并在国产软件GRISYS平台上，针对不同的干扰波进行分析，总结针对不同噪音的衰减方法。

正文一、规则干扰波规则干扰波是指有一定的主频和一定视速度的干扰波。

例如面波、声波、线性干扰波、多次波等。

下面就规则干扰波中的面波、声波、多次波和50Hz交流电干扰进行介绍。

（一）面波图1 面波的形成机理及实际地震记录上的面波从震源发出的波动分为两种: 一种是质点振动方向与传播方向一致的波，称为纵波。

另一种是质点振动方向与传播方向垂直的波，称为横波。

纵波的传播速度较快，在远离震源的地方这两种波动就分开，纵波先到，横波次之。

因此纵波又称P波，横波又称S波。

在没有边界的均匀无限介质中，只能有P波和S波存在，它们可以在三维空间中向任何方向传播，所以叫做体波。

但地球是有限的，有边界的，在界面附近，体波衍生出另一种形式的波，它们只能沿着界面传播，只要离开界面即很快衰减，这种波称为面波。

面波实际上是体波在地表衍生而成的次生波, 面波是一种很强并广泛存在的规则干扰波 ,在炮集上呈线性分布 ,其特征为低频、低速且振动延续时间长 ,严重影响中深层有效反射 ,大大降低地震资料的信噪比，如图1所示。

海洋资料多次波组合衰减技术及应用

海洋资料多次波组合衰减技术及应用张兴岩;朱江梅;杨薇;于宏;刘洋波【摘要】多次波的压制是海洋地震资料处理中的重要步骤.海洋地震勘探中的多次波一般都异常发育且较难压制,压制多次波的好坏直接决定了处理结果的质量.目前压制多次波的方法有很多种,但不同的方法基于不同的理论,具有不同的针对性.单一的方法压制多次波往往不能针对所有多次波,达不到理想的效果.笔者使用组合衰减多次波技术,综合了SRME、拉东变换、LIFT去噪对海洋实际地震数据的多次波进行压制,取得了较好的效果.%Antimultiple processing is an important step in marine seismic data processing. The multiple waves in marine seismic data are generally rather strong and very difficult to suppress, so the suppression of multiple waves decides the quality of processing results. There are many ways to suppress multiple waves; nevertheless, different approaches are based on different theories and have different targets. A single method can not suppress all multiple waves and hence fails to achieve the desired results. The authors employed the combined multiple wave attenuation technology, which combines SRME, radon transform and lift denoising to suppress multiple waves in actual marine seismic data, and has achieved good results.【期刊名称】《物探与化探》【年(卷),期】2011(035)004【总页数】5页(P511-515)【关键词】多次波;SRME;拉东变换;LIFT【作者】张兴岩;朱江梅;杨薇;于宏;刘洋波【作者单位】中海油能源发展钻采工程研究院地球物理研究所,广东湛江 524057;中海油能源发展钻采工程研究院地球物理研究所,广东湛江 524057;中海油能源发展钻采工程研究院地球物理研究所,广东湛江 524057;中海油能源发展钻采工程研究院地球物理研究所,广东湛江 524057;中海油能源发展钻采工程研究院地球物理研究所,广东湛江 524057【正文语种】中文【中图分类】P631.4在海洋地震勘探中，多次波的存在会影响地震成像的真实性和可靠性，并严重影响地震解释工作，多次波压制是提高地震资料处理成果质量的重要步骤之一。

海上地震资料多次波剔除方法研究

海上地震资料多次波剔除方法研究摘要：一直以来，消除多次波等相关噪音是海上地震资料处理要解决的问题。

对海上地震资料多次波剔除方法进行研究，以期为海上地震资料消除多次波提供参考。

关键词：海上地震资料多次波剔除中图分类号：p631.4 文献标识码：a 文章编号：1007-3973（2013）003-119-02在处理海上地震资料时，人们常被把没有消除的多次波错误地认为是一次波或是一次波中的一部分，这种处理方法极大程度地降低了海上地震资料信噪比。

目前，压制或消除多次波的方法很多，常见的有radon滤波、预测反褶积、 -p域反褶积或f-k域滤波等，这些滤波剔除法有其独有的优点，但不容忽视的是，上述方法若使用不当会加重对有效信号的损伤。

本文简要论述radon滤波、预测反褶积和f-k域滤波压制多次波的方法，并详细探讨聚束滤波压制多次波。

1 多次波概述在海上地震勘探时常常把弯曲界面反射波、断面反射波和绕射波定义为异常波，或把多次反射波称为异常波，但多次反射波和弯曲界面反射波、断面反射波和绕射波存在差异，多次反射波是一种单纯的干扰波，但其他几类波除了有干扰之外，还有可利用的一面。

本文讨论的多次波仅限于多次反射波。

地面和空气的分界面阻抗差别较明显，因此，地面和空气的分界面是一个良好的反射界面。

一般情况下，反射波传播至地面时，反射波可能从地面与空气的分界面发射向下传播，再次遇到反射界面的时候，可又一次发射回到地面，反射路径重复多次，就形成了多次反射波。

良好的反射界面是产生多次反射波的前提条件，倘若反射界面反射系数过小，经多次反射后，多次波就变得很微弱，所以，只有那些在发射系数较大的反射界面发生的多次反射波才可以被记录下来。

较为常见的多次波有层间多次波和表层多次波，层间多次波多发生于陆地，这里我们主要讨论表层多次波。

表层多次波多见于海上地震资料，分为水面和海底之间震荡产生的交混回响多次波和海底横向高差变化剧烈产生的绕射多次波。

地震资料处理流程与方法介绍(2)

实际处理中解决拉伸畸变的直接办法就是切除。
动校正前
动校正后
3、水平叠加
九、动校正、切除与叠加
叠加
同一反射点地震记录
叠加剖面
十、（短波长）剩余静校正
1、为什么要做剩余静校正
由于低速带的速度和厚度在横向上的变化，使野外表层参数测量不准确或无法测量，故使野外静校正后，爆炸点和接收点的静校正量还残存着或正或负的误差，这个误差称为“剩余静校正量”。
幅能量分布均匀合理。基本假设：近地表不均匀因素对地震记录影响十分复杂，把各种因素同时加以考
虑会使问题变得十分棘手，甚至无法解决。为了使问题简化并满足地表一致性要求，一般作如下假设：
（1）地表振幅影响因子对整道是一个常数，它是震源强度、表层衰减、检波器耦合等影响的总和系数。
（2）各振幅因子保持地表一致性原则。即不管波的传播路径如何，同一道集内所有道将具有同一补偿因子。如：同一炮的所有道将具有同一炮点的补偿因子，同一检波点所有道将具有同一检波点的补偿因子。
将野外磁带数据转换成处理系统格式，加载到磁盘上； 2、输入数据质量检查：炮号、道号、波形、道长、采样间隔等等。
二、置道头
道头：每个地震道的开始部分都有一个固定字节长度的空余段，这个空余段用来记录
描述本道各种属性的信息，称之为道头。如第8炮第2道，第126CMP等。
1、观测系统定义
模拟野外，定义一个相对坐标系，将野外的激发点、接收点的实际位置放到这个相对的坐标系中。 2、置道头
（3）输入数据为经准确的静校正、球面扩散、地层衰减补偿后的记录。 ——可以根据数据的具体情况，在处理的不同阶段多次使用。目前的流程大都使用一次。
2、地表一致性振幅补偿
五、振幅补偿

双曲线radon变换衰减多次波技术在和15工区的应用

双曲线radon变换衰减多次波技术在和15工区的应用双曲线radon变换衰减多次波技术在和15工区的应用1. 介绍•双曲线radon变换衰减多次波技术是一种常用于勘探地震学领域的信号处理方法，通过消除或减弱地震数据中的多次反射波来提高地震图像的质量和解释能力。

•在和15工区中，双曲线radon变换衰减多次波技术已经被广泛应用于勘探活动中，下面将列举几个具体的应用示例。

2. 地层切割•双曲线radon变换衰减多次波技术可以通过对地震数据进行处理，将多次反射波从数据中去除，从而清晰地展示目标地层的信息。

•在和15工区的勘探活动中，地层切割是一项重要的任务，通过双曲线radon变换衰减多次波技术，可以将地震图像中的多次反射波减弱或消除，从而更好地识别目标地层。

3. 目标检测•双曲线radon变换衰减多次波技术可以帮助识别和定位潜在的油气储层。

通过去除地震图像中的多次反射波，可以提高图像的清晰度和对比度，从而更好地识别出潜在的目标层位。

•在和15工区的勘探活动中，目标检测是一项重要的任务，通过双曲线radon变换衰减多次波技术，可以辅助勘探人员更准确地定位和评估潜在的油气储层。

4. 缺陷识别•双曲线radon变换衰减多次波技术可以帮助勘探人员识别地质构造中的缺陷或异常。

通过分析地震数据中的多次反射波，可以发现一些异常反射波形态，并进一步推断可能存在的地质构造或油气藏。

•在和15工区的勘探活动中，缺陷识别对于油气勘探的成功至关重要。

双曲线radon变换衰减多次波技术可以帮助勘探人员在数据处理过程中更好地识别出地质中的缺陷。

5. 资源评估•双曲线radon变换衰减多次波技术可以提高地震图像的质量和清晰度，从而准确评估潜在的油气资源量。

通过去除多次反射波，可以获得更准确的地震数据，进而提高对油气资源量的估计准确度。

•在和15工区的勘探活动中，资源评估是一项关键任务，通过双曲线radon变换衰减多次波技术，可以帮助勘探人员更准确地评估油气资源的含量和分布情况。

几种简单实用的多次波去除方法及其应用

几种简单实用的多次波去除方法及其应用张恒超【摘要】多次波衰减是地震资料处理中的老问题,多次波的存在严重影响地震剖面成像的真实性和可靠性,干扰地震资料的解释,所以,在叠前应对多次波加以有效的消除,提高资料的品质和精度.而目前,多次波理论方面的研究及成果相当多,但对于实际生产有时却无能为力.因此,本文意在从原理和可操作性上,介绍几种在实际生产中既简单实用又省时有效的多次波去除方法及其应用效果,主要介绍了数据一致性预测反褶积、Radon变换、利用速度扫描工具及内切除等四种多次波去除方法及其应用效果.以提高资料的处理质量和效率.【期刊名称】《内蒙古石油化工》【年(卷),期】2011(037)021【总页数】4页(P21-24)【关键词】地震资料处理;多次波;去除;数据一致性预测反褶积;Radon变换;速度扫描;内切除【作者】张恒超【作者单位】中海油研究总院,北京,100027【正文语种】中文【中图分类】TP273由于地下强反射界面的存在,产生大量的多次波。

众所周知,多次波衰减是地震资料处理中的老问题,多次波也一直是处理中常见的干扰波。

由于地表及地下结构的变化,多次波的周期、频率、分布规律等具有多变性,并且常常和一次有效反射波相干涉,使得多次波问题复杂化。

多次波的存在,严重影响地震剖面成像的真实性和可靠性,干扰地震资料的解释,所以,在叠前应对多次波加以有效的消除,提高资料的品质和精度[1]。

多次波问题在我国一些地区还相当严重,因此,如何有效地压制多次波是地震勘探中的一个关键问题。

目前,多次波理论方面的研究及成果相当多,但对于实际生产有时却无能为力。

因此,本文意在介绍几种在实际生产中既简单实用又省时有效的多次波去除方法。

压制多次波的方法有很多,但其技术均是基于多次波下列特性之一:①一次波与多次波时差的区别(速度区别);②一次波与多次波在CM P叠加剖面上倾角的区别;③一次波与多次波频率成分的区别;④多次波的周期性。

利用地震干涉法衰减海底相关层间多次波

利用地震干涉法衰减海底相关层间多次波叶月明;姚根顺;赵昌垒;庄锡进;张华【摘要】在海洋地震勘探中,由于海底面具有较强的波阻抗差,当海底面以下存在高速反射层或散射体时,在海底与高速层(体)间会发育较强的海底相关层间多次波.由于这种相关层间多次波无特定的识别特征而难以去除,增加了地震资料处理的难度.为此,本文提出了基于地震干涉理论的海底相关层间多次波衰减方法和实现流程.首先介绍了基于互易理论的地震干涉法原理,阐述了其在非均匀介质中的近似条件和表达式;在此基础上提出了通过海底一次反射记录和海底以下构造的反射记录来构建海底相关层间多次波.该方法克服了基于共聚焦点(CFP)延拓方法对速度场的依赖性与逆散射级数(ISS)法计算效率低的缺点,是一种数据驱动型层间多次波衰减方法.Sigsbee2B模型测试及实际资料处理结果表明,海底相关层间多次波得到了较好的预测,叠前深度偏移结果中的偏移假象得到了明显压制,证明了文中方法的有效性.【期刊名称】《石油地球物理勘探》【年(卷),期】2015(050)002【总页数】7页(P225-231)【关键词】海底相关层间多次波;衰减;地震干涉法;互易理论;数据驱动;格林函数【作者】叶月明;姚根顺;赵昌垒;庄锡进;张华【作者单位】中国石油杭州地质研究院,浙江杭州310023;中国石油杭州地质研究院,浙江杭州310023;中国石油杭州地质研究院,浙江杭州310023;中国石油杭州地质研究院,浙江杭州310023;川庆钻探工程有限公司地球物理勘探公司技术发展中心,四川成都610213【正文语种】中文【中图分类】P6311 引言多次波分为地表相关多次波和层间多次波。

传统地表相关多次波衰减（SRME）方法可以较好地衰减自由表面相关多次波［1，2］，并且已经在业界得到了广泛应用。

由于地下界面多变性的特征，很难在地震剖面上有效识别层间多次波。

现阶段业界主要应用基于时差差异的方法进行层间多次波压制［3，4］，这种方法对于浅层适用性较好，但是难以区分深层层间多次波，而且对处理员的主观解释依赖性较强。

地震资料处理中的迭代法衰减多次波技术

ｔｉｏｎｓ，ｔｈｉｓｍｅｔｈｏｄｃａｎａｔｔｅｎｕａｔｅｍｕｌｔｉｐｌｅｗａｖｅｓ，ａｔｔｈｅｓａｍｅｔｉｍｅｔｒｉｅｓｔｏｐｒｏｔｅｃｔｔｈｅｅｆｆｅｃ — ｔｉｖｅｓｉｇｎａｌ，ｗｈｉｃｈｃａｎａｃｈｉｅｖｅｔｈｅｐｕｒｐｏｓｅｔｈａｔｃａｎｓｕｐｐｒｅｓｓｍｕｌｔｉｐｌｅｗａｖｅｓａｎｄｋｅｅｐｔｈｅ
同时，把噪声衰减到可以接受的程度。这种思路无论是用于去噪还是去除多次波，都有良好的应用前景。
关键词：多次波；剩余多次波；迭代法；地震资料处理
中图分类号：Ｐ６３１．４
文献标识码：Ａ
ｅｆｆｅｃｔｉｖｅｓｉｇｎａｌｓｉｎｍａｘｉｍｕｍ．Ｉｔｉｓｎｏｔｊｕｓｔ３ｓｉｎｇｌｅｍｏｄｕｌｅ，ｂｕｔａ“ ｓｉｇｎａｌｎｏｉｓｅｓｅｐａｒａｔｉｏｎ ”
ｆｒｏｍｔｈｅｒｅｍｏｖｅｄｎｏｉｓｅａｎｄｒｅｃｏｎｓｔｒｕｃｔｄａｔａｗｉｔｈｔｈｅｅｆｆｅｃｔｉｖｅｓｉｇｎａ１．Ｂａｓｅｄｏｎｍｕｌｔｉ — ｉｔｅｒａ
Ａｂｓｔｒａｃｔ：Ｍｕｌｔｉｐｌｅａｔｔｅｎｕａｔｉｏｎｔｅｃｈｎｉｑｕｅｉｓｏｎｅｏｆｔｈｅｋｅｙｔｅｃｈｎｏｌｏｇｉｅｓｉｎｔｈｅｃｏｎｖｅｎｔｉｏｎａｌｐｒｏｃｅｓｓｉｎｇｏｆｓｅｉｓｍｉｃｄａｔａ．Ａｔｐｒｅｓｅｎｔ，ｍｏｓｔｏｆｔｈｅｍｅｔｈｏｄｓｆｏｒｔｈｅａｔｔｅｎｕａｔｉｏｎｏｆｍｕｌｔｉｐｌｅｗａｖｅｓｗｉｌｌｄａｍａｇｅｔｈｅｅｆｆｅｃｔｉｖｅｓｉｇｎａ１．Ｍｕｌｔｉｐｌｅａｔｔｅｎｕａｔｉｏｎｔｅｃｈｎｉｑｕｅｂｙｉｔｅｒａｔｉｖｅｍｅｔｈｏｄｉｓｕｓｅｄａｓａｍｅａｎｓｔｏｓｕｐｐｒｅｓｓｍｕｌｔｉｐｌｅｗａｖｅｓ，ｗｈｉｃｈｃａｎｅｘｔｒａｃｔｔｈｅｗｅａｋｅｆｆｅｃｔｉｖｅｓｉｇｎａｌ

地震如何利用地震波幅度衰减地震灾害程度

地震如何利用地震波幅度衰减地震灾害程度地震是一种自然灾害，经常给人们的生活带来极大的危害和损失。

然而，地震波的传播和幅度衰减特性使得我们有机会通过科学手段来减轻地震灾害的程度。

本文将探讨地震波的幅度衰减原理以及利用幅度衰减来减轻地震灾害的方法。

一、地震波的幅度衰减原理地震波是在地震破裂源处产生，并在地壳中传播的振动波动。

地震波幅度衰减是指地震波振幅随距离的增加而逐渐减小的现象。

这种衰减现象受到多种因素的影响，包括地震波类型、地壳介质性质和路径长度等。

地震波的幅度衰减主要受到能量耗散和波传播衰减两个方面的影响。

能量耗散是指地震波在传播过程中由于多次反射、散射和吸收而减小能量。

波传播衰减是指地震波在介质中传播过程中由于散射、衍射和衰减等因素而逐渐减小振幅。

二、利用地震波幅度衰减减轻地震灾害的方法1. 地震预测与监测利用地震波幅度衰减的特性，科学家们可以通过监测和记录地震波在不同距离下的振幅变化来预测地震发生的可能性以及地震烈度的分布情况。

这可以帮助人们及时采取措施来减轻地震灾害的程度。

2. 建筑物抗震设计在建筑物的设计过程中，可以利用地震波的幅度衰减特性来选择适当的抗震设计参数和结构形式。

例如，对于地震烈度较大的地区，设计师可以采用更加耐震的结构形式和材料，以减轻地震对建筑物的破坏程度。

3. 地质灾害防治地震波的幅度衰减特性也可以应用于地质灾害的防治工作中。

例如，山区地震波传播中的幅度衰减可以帮助人们预测滑坡、泥石流等地质灾害的危险程度，并及时采取相应的防范措施。

4. 救援和应急响应当地震发生后，救援人员可以利用地震波幅度衰减的特性来判断受灾区域和程度。

这有助于救援人员更加准确地评估灾情，并有针对性地展开救援工作，最大限度地减少人员伤亡和财产损失。

5. 科学研究与技术创新地震波幅度衰减的研究还可以推动科技领域的发展与创新。

在使用新材料、新技术进行地震波传播模拟与仿真、构建地震波传播衰减模型等方面，科学家们可以不断改进地震灾害预测与减灾的能力。

多次波去除

地震数据中多次波去除当反射波向上传播遇到某个具有较大波阻抗差的界面时，反射波有可能从这个界面反射，向下传播，当遇到反射界面时，再次发生反射返回地（海）面，于是就形成了多次波。

包括一次多次波、多次多次波，层间多次波，而层间多次波混杂在有效波之间，最不容易识别和消除。

除的基本方法是利用多次波的可预测性以及多次波和有效波的速度差异进行消除。

在CDP 叠加过程中，针对多次波衰减，我们常用三种方法（老）：1.基于相同的到达时间，一次波与多次波正常时差有明显差异，我们分别有一次波波与多次波速度进行动校正和叠加，然后相见，把多次波从一次波叠加剖面中消除，即减去法消除多次波。

2.基于多次波旅行时间在零炮间距道上，具有周期性的特点，选择合适的步长，采用预测反褶积的方法。

3.基于增益平衡补偿原理，提高一次波强度，削弱多次波。

近几年提出的方法1.FK域压制多次波对数据进行NMO校正，校正速度介于一次波和多次波之间，这样一次波校正过头，多次波校正不足，校正后记录在时空域一次波上饶，多次波仍然下弯。

把校正后的数据进行二维傅里叶变换，由于一次波与多次波视速度的符号不同，可以在FK域分别位于零波数轴的两侧。

设计FK滤波器，对数据进行滤波，然后对滤波后的数据进行二维反傅里叶变换，返回时空域。

用相同的速度进行反动校正。

2.波动方程外推法压制海底多次波利用波动方程对波场进行外推，压制海底产生的多次波，首先对对水面上采集的数据波场外推到水底面，然后去掉多次反射，再把结果外推到水面上，因此，他的中心思想是通过波动方程来预测多次波预测。

如图db是下行波，ub是上行波，ur是海底下面反射界面上产生的上行波，我们可以看到在上行波ub中，即包含了下行波db在水底界面产生的多次反射，也包含了上行波ur在水底界面上的投射能量，前者是我们要消除的多次波，后者是我们要保留的一次反射则ur’=ub-r*ur上式中r代表水底界面的反射系数，当r只是一个标量的时候，*可以理解为一个乘号，当r 为一个反射系数序列的时候，*表示一个褶积符号。

浅水多次波衰减技术在多道地震数据处理中的应用

DOI: 10.16562/ki.0256-1492.2018101202浅水多次波衰减技术在多道地震数据处理中的应用杨佳佳1,2，潘军1,2，栾锡武1,2，颜中辉1,2，刘鸿1,21. 自然资源部油气资源和环境地质重点实验室，中国地质调查局青岛海洋地质研究所，青岛 2660712. 海洋国家实验室海洋矿产资源评价与探测技术功能实验室，青岛 266071摘要：拖缆采集的海上多道地震数据受海水间传播的多次波影响，往往在共炮点道集上发育周期性的强振幅干扰波，这些噪音会掩盖海底以下及强反射层以下地层的面貌，严重影响多道地震数据的成像效果。

为了消除海上多道地震数据中存在的海水间传播的多次波对地震数据的干扰，首先分析浅水多次波在海水间多次震荡的产生机制，采用了τ-P 域静校正延迟技术对其进行压制，处理的核心是在τ-P 域对周期性多次波模型进行预测，再通过自适应相减以去除地震数据中的多次波。

实际资料处理结果表明，方法对海水间震荡的浅水多次波具有很好的压制效果，经过浅水多次波去除后的叠加剖面信噪比得到有效提高，剖面的品质得以提升。

关键词：浅水多次波衰减；多道地震；τ-P 变换中图分类号：P738 文献标识码：AApplication of attenuation technology to shallow water multiples in multi-channel seismic data processingYANG Jiajia 1,2, PAN Jun 1,2, LUAN Xiwu 1,2, YAN Zhonghui 1,2, LIU Hong 1,21. Key Laboratory of Marine Environmental Geology, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China2. Laboratory for Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, ChinaAbstract: The multi-channel seismic data acquired by cables are often affected by multiple waves propagating between seawater and periodic strong amplitude interferences on common shot gathers, which include the seabed and the strong reflections below it. Some solutions to eliminate the interference of multiple waves are discussed in this paper. The generation mechanism of the shallow water multiple waves i.e. the multiple oscillations between sea-level and seabed are firstly discussed and analyzed. Then statics technique is used to suppress the multiples in the τ-p domain according to the vibrating characteristics of the multiples in the sea. The core of the method is to predict the cycling multiples in the τ-p domain and attenuate them by adaptive subtract. The results demonstrate that the method is effective in suppression of shallow water multiple waves caused by sea-water oscillation, the S/N ratio is obviously improved and the profile is better than that before the processing.Key words: shallow water de-multiple; multi-channel seismic; τ-p transform海底孕育了丰富的油气资源，近年来海洋油气勘探得到迅速发展，海上拖缆采集是高效的获取海洋地震资料的手段，由于海水面是一个强波阻抗界面，地震波传播到海水面产生的反极性强反射的多次波[1]是海洋地震资料处理中面临的一大难题。

Radon变换在浅层地震数据多次波衰减中的应用

第４４卷㊀第１－２期２０２０年６月地质学刊ＪｏｕｒｎａｌｏｆＧｅｏｌｏｇｙＶｏｌ．４４Ｎｏ．１－２Ｊｕｎ．，２０２０ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７４－３６３６．２０２０．ｈ１．０１５Ｒａｄｏｎ变换在浅层地震数据多次波衰减中的应用罗传根，贺剑波，刘屹立，王振宇（江苏省地质勘查技术院，江苏南京２１００４９）摘要：在浅层地震勘探中，多次波常会造成构造假象，给地震资料解释工作带来极大的困难㊂基于信号理论的叠前和叠后预测反褶积㊁Ｒａｄｏｎ变换㊁ｆ－ｋ滤波㊁聚束滤波等多次波衰减方法，以及基于波动方程理论的波场外推法㊁逆散射级数法㊁反馈迭代法等多次波预测相减法，是目前去除多次波的主要手段，各种方法都有一定的适用性㊂主要探讨了Ｒａｄｏｎ变换方法在浅层地震资料多次波衰减中的应用效果，模型试算和实际资料应用证明，当１次波与多次波的剩余时差＞３０ｍｓ时，Ｒａｄｏｎ滤波方法效果最好㊂关键词：浅层地震；多次波；Ｒａｄｏｎ变换；数据处理；模型试算；应用实例中图分类号：Ｐ６３１４＋４３㊀㊀㊀文献标识码：Ａ㊀㊀㊀文章编号：１６７４－３６３６（２０２０）０１－０１３２－０９收稿日期：２０１９－１０－１０；修回日期：２０２０－０２－１１；编辑：蒋艳基金项目：江苏省地质勘查基金项目南通市城市地质调查（苏财建［２０１７］１４０）作者简介：罗传根（１９６５㊀），男，高级工程师，勘查技术与工程专业，主要从事浅层地震勘探工作，Ｅ⁃ｍａｉｌ：７０６４０９６０７＠ｑｑ．ｃｏｍ通信作者：贺剑波（１９８８㊀），男，工程师，硕士，地球探测与信息技术专业，Ｅ⁃ｍａｉｌ：３６１７９１８８４＠ｑｑ．ｃｏｍ０㊀引㊀言正确识别并区分干扰波和有效波是开展地震资料处理工作的前提，１次反射波地震勘探中通常将面波㊁声波㊁浅层折射波㊁多次波等作为干扰波来处理，其中面波㊁声波㊁浅层折射波等规则干扰波的频率及速度等与１次反射波差异明显，因此比较容易去除㊂但多次反射波的产生与地下地层结构㊁地质构造形态有关，产生机制复杂，种类繁多，因此增加了资料处理时的识别与去除难度，且多次波的存在会在速度谱上形成假能量团，对速度分析造成误导，偏移叠加剖面上会出现构造假象，进而影响后续资料解释的正确性，所以多次波的衰减与去除是目前地球物理工作研究的重点之一㊂目前常用的多次波压制方法主要基于信号理论或波动方程理论：基于信号理论的多次波压制方法有叠前和叠后预测反褶积㊁Ｒａｄｏｎ变换㊁ｆ－ｋ滤波㊁聚束滤波等，主要根据１次波和多次波之间的时差㊁视速度㊁周期性不同等特征进行分离去除；基于波动方程理论的多次波压制方法有波场外推法㊁逆散射级数法㊁反馈迭代法等，主要利用地震波的动力学和运动学特征等传播信息进行多次波预测压制（刘建辉，２０１０）㊂对浅层地震资料中的多次波进行衰减处理主要利用Ｒａｄｏｎ变换，国外多应用于多次波衰减（Ｈａｍｐｓｏｎ，１９８６；Ｄｏｕｇｌａｓｅｔａｌ．，１９９２），即根据多次波与１次波同相轴在Ｒａｄｏｎ空间的差异性分布，实现多次波与１次波的分离，通过切除多次波能量团并反变换回数据空间，达到压制多次波的目的，我国也进行了这方面的深入研究（朱生旺等，２００２；刘喜武等，２００４；张军华等，２００４；黄新武等，２００５；熊登等，２００９）㊂目前，Ｒａｄｏｎ变换去除多次波的应用主要集中在石油勘探㊁煤田勘探等深层反射地震资料处理方面（褚玉环等，２００９），应用效果有待进一步验证㊂１㊀Ｒａｄｏｎ变换的基本理论Ｒａｄｏｎ变换将地震数据从时空域映射到Ｒａｄｏｎ域空间，１次波和多次波因其速度差异形成不同的能量团，进而可以在Ｒａｄｏｎ空间中去除多次波能量团㊂地震学中通常将Ｒａｄｏｎ变换定义为下式：ｕ（τ，ｐ）＝ʏ＋ɕ－ɕｄ（ｔ，ｘ）ｄｘ（１）第４４卷㊀第１－２期罗传根㊀等：Ｒａｄｏｎ变换在浅层地震数据多次波衰减中的应用１３３㊀式（１）中，ｄ（ｔ，ｘ）为时间域地震数据；ｕ（τ，ｐ）表示Ｒａｄｏｎ域的数据值；ｘ为空间变量值，２Ｄ空间中表示Ｘ方向的偏移距信息；ｐ为曲线的曲率，线性情况下ｐ＝ｓｉｎθｖ－１，即射线参数，抛物线情况下ｐ＝１／２ｔ０ｖ２；τ为时间轴的截距值㊂ｔ为双程旅行时间值，若ｔ＝τ＋ｐｘ，则为线性Ｒａｄｏｎ变换；若ｔ＝τ＋ｐｘ２，则为抛物线型Ｒａｄｏｎ变换；若ｔ＝τ２＋ｐ２ｘ２，则为双曲线型Ｒａｄｏｎ变换㊂将双曲积分算子ｔ＝τ２＋ｐ２ｘ２代入式（１）并离散化，得到２Ｄ双曲型Ｒａｄｏｎ正变换的离散表达式：ｕ（τ，ｐ）＝ðｘ＝ｘｍａｘｘ＝ｘｍｉｎｄ（ｔ２＝τ２＋ｘ２ｐ２，ｘ）（２）相应得到式（２）反变换的离散公式：ｄ（ｔ，ｘ）＝ðｐ＝ｐｍａｘｐ＝ｐｍｉｎｕ（τ２＝ｔ２－ｘ２ｐ２，ｐ）（３）将式（２）㊁（３）写成矩阵表达式：ｕ＝ＬＴｄ（４）ｄᶄ＝Ｌｕ（５）式中，ｄ为ｔ－ｘ域原始数据，ｄᶄ为Ｒａｄｏｎ反变换后的ｔ－ｘ域数据，ｕ为Ｒａｄｏｎ域数据，Ｌ为正变换算子，Ｔ为矩阵的共轭转置，Ｌ和ＬＴ通常不可逆，常用最小平方（最小二乘）法计算其广义逆㊂定义ｅ＝ｄ－ｄᶄ＝ｄ－Ｌｕ为剩余量，则ｅ的最小平方误差为：ｓ＝ｅＴｅ（６）对于式（６）中ｕ的求解，即最小平方误差ｓ的求解分２种情况：①超定情况下，ｕ＝（ＬＴＬ）－１ＬＴｄ；②欠定情况下，ｕ＝ＬＴ（ＬＬＴ）－１ｄ㊂其中，（ＬＴＬ）－１和（ＬＬＴ）－１２个算子在变换中可以保持较高的分辨率，起反褶积作用㊂实际情况下，为了正变换的稳定性，会在变换算子中加入白噪，即（ＬＴＬ＋αＩ）－１和（ＬＬＴ＋αＩ）－１㊂其中，Ｉ是单位向量，α一般取ＬＨＬ最大特征值平方根的１％（Ｈ为矩阵的共轭）㊂２㊀Ｒａｄｏｎ滤波Ｒａｄｏｎ滤波是利用多次波和１次波的速度差异将输入数据变换到Ｒａｄｏｎ域，实现１次波与多次波的分离，再将多次波能量团切除，最终得到时空域去除多次波的数据，实现过程为动校正＋Ｒａｄｏｎ正变换＋Ｒａｄｏｎ域多次波去除＋Ｒａｄｏｎ反变换＋反动校正的数据映射（图１）㊂图１㊀Ｒａｄｏｎ滤波流程图Ｆｉｇ．１㊀ＦｌｏｗｃｈａｒｔｏｆＲａｄｏｎｆｉｌｔｅｒｉｎｇ３㊀模型试算３．１㊀多次波ＣＭＰ道集试演首先正演一个发育２组多次波的共中心点（ＣＭＰ）道集来验证方法的有效性㊂２５０ｍｓ处发育１组反射波，速度为３０００ｍ／ｓ；３００ｍｓ处发育１组多次波，速度为２８００ｍ／ｓ；５００ｍｓ处发育１组多次波，速度为２５００ｍ／ｓ㊂３组反射波子波频带范围均为５１５４０５０Ｈｚ，２００道接收，道间距为１０ｍ㊂正演记录经反射波速度（３０００ｍ／ｓ）动校正后（图２）显示，该ＣＭＰ记录中反射波同相轴被校平，而多次波同相轴欠校正㊂对动校正后的ＣＭＰ道集进行Ｒａｄｏｎ变换（图３）㊂１次反射波（２５０ｍｓ）出现在剩余时差（Ｍｏｖｅｏｕｔ）为０附近，多次波（３００㊁５００ｍｓ）出现在其相应的剩余时差附近㊂１次反射波和多次波在Ｒａｄｏｎ域实现了分离，沿着两者能量团的分界线拾取１条界面（图３中红线所示），就能将两者分离开来㊂Ｒａｄｏｎ变换模拟出的ＣＭＰ道集和Ｒａｄｏｎ滤波去除多次波以后的ＣＭＰ道集（图４）显示，Ｒａｄｏｎ变换模拟出的ＣＭＰ道集上的１次反射波和多次波与原始输入道集一致，说明Ｒａｄｏｎ变换的参数是合理的，Ｒａｄｏｎ滤波后的道集显示多次波得到了明显的衰减㊂１３４㊀地㊀质㊀学㊀刊２０２０年图２㊀动校正后含多次波ＣＭＰ道集Ｆｉｇ．２㊀ＣＭＰｇａｔｈｅｒｓｗｉｔｈｍｕｌｔｉｐｌｅｗａｖｅｓａｆｔｅｒｎｏｒｍａｌｍｏｖｅｏｕｔ图３㊀Ｒａｄｏｎ域时间－剩余时差关系Ｆｉｇ．３㊀ＲｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎｔｉｍｅａｎｄｒｅｍａｉｎｉｎｇｔｉｍｅｍｏｖｅｏｕｔｄｉａｇｒａｍｉｎＲａｄｏｎｄｏｍａｉｎ图４㊀Ｒａｄｏｎ变换模拟道集及Ｒａｄｏｎ滤波后道集Ｆｉｇ．４㊀ＭｏｄｅｌｅｄｇａｔｈｅｒｕｓｉｎｇＲａｄｏｎｔｒａｎｓｆｏｒｍａｎｄｇａｔｈｅｒｓａｆｔｅｒＲａｄｏｎｆｉｌｔｅｒｉｎｇ第４４卷㊀第１－２期罗传根㊀等：Ｒａｄｏｎ变换在浅层地震数据多次波衰减中的应用１３５㊀３．２㊀有效反射波和多次波ＣＭＰ道集试演再正演１个含２组有效反射波和２组多次波的共中心点（ＣＭＰ）道集来验证方法的有效性㊂２５０㊁８００ｍｓ处分别发育１组反射波，速度为３０００ｍ／ｓ；３００ｍｓ处发育１组多次波，速度为２８００ｍ／ｓ；８００ｍｓ处发育１组多次波，速度为２５００ｍ／ｓ㊂４组反射波子波频带范围均为５１５４０５０Ｈｚ，２００道接收，道间距为１０ｍ㊂正演记录经反射波速度（３０００ｍ／ｓ）动校正后（图５）显示，该ＣＭＰ记录中反射波同相轴被校平，而多次波同相轴欠校正㊂对动校正后的ＣＭＰ道集进行Ｒａｄｏｎ变换（图６）㊂１次反射波（２５０㊁８００ｍｓ）出现在时差（Ｍｏｖｅｏｕｔ）为０附近，多次波（３００㊁８００ｍｓ）出现在其相应的剩余时差附近㊂１次反射波和多次波在Ｒａｄｏｎ域实现了分离，沿着两者能量团的分界线拾取１条界面（图６中红线所示），就能将两者分离开来㊂Ｒａｄｏｎ变换模拟出的ＣＭＰ道集和Ｒａｄｏｎ滤波去除多次波后的ＣＭＰ道集（图７）显示，Ｒａｄｏｎ变换模拟出的ＣＭＰ道集上的１次反射波和多次波与原始输入道集一致，说明Ｒａｄｏｎ变换的参数是合理的㊂从Ｒａｄｏｎ滤波后的道集上可以看到多次波得到了明显的衰减㊂４㊀实际测试利用上述Ｒａｄｏｎ变换流程对陆上某工作区的二维资料进行多次波衰减㊂该条测线采用可控震源车作为震源，扫描频带为１０１２０Ｈｚ，炮间距为１２ｍ，道间距为４ｍ，２８８道接收，ＣＭＰ范围为２０００４０００，采样数为３００１点，采样率为０．５ｍｓ㊂采用初步拾取的速度对原始ＣＭＰ道集进行动校正后的道集（图８）显示，２００ｍｓ附近的１次反射波同相轴被校平，而３００ｍｓ附近的多次波同相轴欠校正㊂采用上述流程进行Ｒａｄｏｎ滤波后的道集（图９）显示，３００ｍｓ附近的多次波得到了较好的衰减㊂对比基于原始ＣＭＰ道集的速度分析得到的速度谱（图１０㊁图１１）和基于Ｒａｄｏｎ滤波后的ＣＭＰ道集得到的速度谱（图１２㊁图１３），发现Ｒａｄｏｎ滤波后得到的速度谱上多次波的能量团消失，１次反射波的能量团更加明显，速度谱更加干净，减少了速度谱拾取中的陷阱㊂利用初步拾取的速度ｖ０得到的叠加时间剖面（图１４）显示，３００ｍｓ附近发育１组基岩多次波，振幅强㊁连续性好，给地震资料解释带来较大的干扰㊂对原始道集利用速度ｖ０进行Ｒａｄｏｎ滤波后再利用速度ｖ０得到的叠加时间剖面（图１５）显示，３００ｍｓ附近发育的基岩多次波得到了一定程度的衰减，由于Ｒａｄｏｎ滤波后速度谱上多次波能量团消失，速度分析更加精确，因此对Ｒａｄｏｎ滤波后的道集重新进行速度分析得到速度ｖ１，基于ｖ１得到的叠加时间剖面（图１６）显示，叠加剖面上基岩多次波得到明显衰减，基岩反射波振幅强㊁连续性好，有利于层位的追踪与划分㊂５㊀结㊀论模型试算及实际测试均验证了Ｒａｄｏｎ滤波方法对浅层地震多次波衰减处理的有效性㊂（１）Ｒａｄｏｎ滤波法衰减多次波是基于１次波与多次波同相轴的曲率差异，因此输入数据应该是与构造信息无关的道集，如ＣＭＰ道集㊂当１次波与多次波的剩余时差＞３０ｍｓ时，Ｒａｄｏｎ滤波方法效果最好，但随着目的层埋深变浅，其效果变差，因此用于解译浅层地震资料难度较大㊂（２）Ｒａｄｏｎ滤波法衰减多次波的效果与在Ｒａｄｏｎ域拾取的切除线密切相关：切除线靠近１次波能量团时，多次波衰减效果好，但１次波能量受损；切除线靠近多次波能量团时，多次波衰减不充分㊂通常，沿２０５０ｍｓ时差线拾取是合适的㊂（３）由于Ｒａｄｏｎ滤波的理论仍然是基于地下水平层假设，因此该方法在实际应用时主要适用于地下构造比较简单的情况；对于地下构造复杂的情况，可尝试多种多次波衰减的方法相互组合，以便取得最佳效果㊂参考文献褚玉环，刘清华，王宝，等，２００９．拉冬变换压制多次波技术在大庆探区的应用［Ｊ］．大庆石油地质与开发，２８（６）：３２５－３２７．黄新武，孙春岩，牛滨华，等，２００５．基于数据一致性预测域压制自由表面多次波［Ｊ］．地球物理学报，４８（１）：１７３－１８０．刘喜武，刘洪，李幼铭，２００４．高分辨率Ｒａｄｏｎ变换方法及其在地震信号处理中的应用［Ｊ］．地球物理学进展，１９（１）：８－１５．１３６㊀地㊀质㊀学㊀刊２０２０年图５㊀动校正后含多次波ＣＭＰ道集Ｆｉｇ．５㊀ＣＭＰｇａｔｈｅｒｓｗｉｔｈｍｕｌｔｉｐｌｅｗａｖｅａｆｔｅｒｎｏｒｍａｌｍｏｖｅｏｕｔ图６㊀Ｒａｄｏｎ域时间－剩余时差关系Ｆｉｇ．６㊀ＲｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎｔｉｍｅａｎｄｒｅｍａｉｎｉｎｇｔｉｍｅｍｏｖｅｏｕｔｄｉａｇｒａｍｉｎＲａｄｏｎｄｏｍａｉｎ图７㊀Ｒａｄｏｎ变换模拟道集及Ｒａｄｏｎ滤波后道集Ｆｉｇ．７㊀ＭｏｄｅｌｅｄｇａｔｈｅｒｕｓｉｎｇＲａｄｏｎｔｒａｎｓｆｏｒｍａｎｄＣＭＰｇａｔｈｅｒｓａｆｔｅｒＲａｄｏｎｆｉｌｔｅｒｉｎｇ㊀第４４卷㊀第１－２期罗传根㊀等：Ｒａｄｏｎ变换在浅层地震数据多次波衰减中的应用１３７图８㊀原始ＣＭＰ道集动校正后示意图Ｆｉｇ．８㊀ＯｒｉｇｉｎａｌＣＭＰｇａｔｈｅｒｓａｆｔｅｒｎｏｒｍａｌｍｏｖｅｏｕｔ图９㊀Ｒａｄｏｎ滤波后道集示意图Ｆｉｇ．９㊀ＧａｔｈｅｒｓａｆｔｅｒＲａｄｏｎｆｉｌｔｅｒｉｎｇ图１０㊀原始ＣＭＰ２１００处速度谱Ｆｉｇ．１０㊀ＶｅｌｏｃｉｔｙｓｐｅｃｔｒｕｍａｔｔｈｅｏｒｉｇｉｎａｌＣＭＰ２１００１３８㊀地㊀质㊀学㊀刊２０２０年图１１㊀Ｒａｄｏｎ滤波后ＣＭＰ２１００处速度谱Ｆｉｇ．１１㊀ＶｅｌｏｃｉｔｙｓｐｅｃｔｒｕｍａｔＣＭＰ２１００ａｆｔｅｒＲａｄｏｎｆｉｌｔｅｒｉｎｇ图１２㊀原始ＣＭＰ２８００处速度谱Ｆｉｇ．１２㊀ＶｅｌｏｃｉｔｙｓｐｅｃｔｒｕｍａｔｔｈｅｏｒｉｇｉｎａｌＣＭＰ２８００图１３㊀Ｒａｄｏｎ滤波后ＣＭＰ２８００处速度谱Ｆｉｇ．１３㊀ＶｅｌｏｃｉｔｙｓｐｅｃｔｒｕｍａｔＣＭＰ２８００ａｆｔｅｒＲａｄｏｎｆｉｌｔｅｒｉｎｇ㊀第４４卷㊀第１－２期罗传根㊀等：Ｒａｄｏｎ变换在浅层地震数据多次波衰减中的应用１３９图１４㊀初始速度叠加剖面Ｆｉｇ．１４㊀Ｓｔａｃｋｅｄｐｒｏｆｉｌｅｗｉｔｈｉｎｉｔｉａｌｖｅｌｏｃｉｔｙ图１５㊀Ｒａｄｏｎ滤波后初始速度叠加剖面Ｆｉｇ．１５㊀ＳｔａｃｋｅｄｐｒｏｆｉｌｅｗｉｔｈｉｎｉｔｉａｌｖｅｌｏｃｉｔｙａｆｔｅｒＲａｄｏｎｆｉｌｔｅｒｉｎｇ图１６㊀Ｒａｄｏｎ滤波后速度更新后叠加剖面Ｆｉｇ．１６㊀ＳｔａｃｋｅｄｐｒｏｆｉｌｅｗｉｔｈｕｐｄａｔｅｄｖｅｌｏｃｉｔｙａｆｔｅｒＲａｄｏｎｆｉｌｔｅｒｉｎｇ１４０㊀地㊀质㊀学㊀刊２０２０年刘建辉，２０１０．基于波动理论压制多次波方法研究［Ｄ］．山东青岛：中国石油大学（华东）．刘学通，李振春，张凯，等，２０１３．基于３Ｄ高分辨率Ｒａｄｏｎ变换的多次波压制方法研究［Ｃ］／／中国地球物理２０１３：第十九专题论文集．北京：中国地球物理学会：６１３．熊登，赵伟，张剑锋，２００９．混合域高分辨率抛物Ｒａｄｏｎ变换及其在衰减多次波中的应用［Ｊ］．石油地球物理学报，５２（４）：１０６８－１０７７．朱生旺，魏修成，李峰，等，２００２．用抛物线Ｒａｄｏｎ变换稀疏解分离和压制多次波［Ｊ］．石油地球物理勘探，３７（２）：１１０－１１５．张军华，吕宁，雷凌，等，２００４．抛物线拉冬变换消除多次波的应用要素分析［Ｊ］．石油地球物理勘探，３９（４）：３９８－４０５．ＤｏｕｇｌａｓＪＦ，ＣｈａｒｌｅｓＣＭ，１９９２．Ｓｕｐｐｒｅｓｓｉｏｎｏｆｍｕｌｔｉｐｌｅｒｅｆｌｅｃ⁃ｔｉｏｎｓｕｓｉｎｇｔｈｅＲａｄｏｎｔｒａｎｓｆｏｒｍ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，５７（３）：３８６－３９５．ＨａｍｐｓｏｎＤ，１９８６．Ｉｎｖｅｒｓｅｖｅｌｏｃｉｔｙｓｔａｃｋｉｎｇｆｏｒｍｕｌｔｉｐｌｅｅｌｉｍｉｎａ⁃ｔｉｏｎ［Ｊ］．ＪｏｕｒｎａｌｏｆｔｈｅＣａｎａｄｉａｎＳｏｃｉｅｔｙｏｆＥｘｐｌｏｒａｔｉｏｎＧｅ⁃ｏｐｈｙｓｉｃｉｓｔｓ，２２（１）：４２２－４２４．ＡｐｐｌｉｃａｔｉｏｎｏｆＲａｄｏｎｔｒａｎｓｆｏｒｍｉｎｓｈａｌｌｏｗｓｅｉｓｍｉｃｄａｔａｐｒｏｃｅｓｓｉｎｇｗｉｔｈｍｕｌｔｉｐｌｅｓａｔｔｅｎｕａｔｉｏｎＬｕｏＣｈｕａｎｇｅｎ，ＨｅＪｉａｎｂｏ，ＬｉｕＹｉｌｉ，ＷａｎｇＺｈｅｎｙｕ（ＧｅｏｌｏｇｉｃａｌＥｘｐｌｏｒａｔｉｏｎＴｅｃｈｎｏｌｏｇｙＩｎｓｔｉｔｕｔｅｏｆＪｉａｎｇｓｕＰｒｏｖｉｎｃｅ，Ｎａｎｊｉｎｇ２１００４９，Ｊｉａｎｇｓｕ，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｉｎｓｈａｌｌｏｗｓｅｉｓｍｉｃｅｘｐｌｏｒａｔｉｏｎ，ｍｕｌｔｉｐｌｅｓａｌｗａｙｓｃａｕｓｅｓｔｒｕｃｔｕｒａｌａｒｔｉｆａｃｔｓ，ｗｈｉｃｈｂｒｉｎｇｓｇｒｅａｔｄｉｆｆｉｃｕｌｔｉｅｓｔｏｓｅｉｓｍｉｃｄａｔａｉｎｔｅｒｐｒｅｔａｔｉｏｎ．Ｔｈｅｓｉｇｎａｌ⁃ｔｈｅｏｒｙ－ｂａｓｅｄｍｕｌｔｉｐｌｅａｔｔｅｎｕａｔｉｏｎｍｅｔｈｏｄｓｉｎｃｌｕｄｉｎｇｐｒｅ⁃ｓｔａｃｋａｎｄｐｏｓｔ⁃ｓｔａｃｋｐｒｅｄｉｃｔｉｏｎｄｅｃｏｎｖｏｌｕｔｉｏｎ，Ｒａｄｏｎｔｒａｎｓｆｏｒｍ，ｆ⁃ｋｄｏｍａｉｎｆｉｌｔｅｒｉｎｇ，ｂｕｎｃｈｉｎｇｆｉｌｔｅｒ，ａｎｄｗａｖｅ⁃ｅｑｕａｔｉｏｎ⁃ｔｈｅｏｒｙ－ｂａｓｅｄｗａｖｅｆｉｅｌｄｅｘｔｒａｐｏｌａｔｉｏｎ，ｉｎｖｅｒｓｅｓｃａｔｔｅｒｉｎｇｓｅ⁃ｒｉｅｓｍｅｔｈｏｄ，ｆｅｅｄｂａｃｋｉｔｅｒａｔｉｏｎａｒｅｔｈｅｍａｉｎｍｅｔｈｏｄｓｏｆｍｕｌｔｉｐｌｅｓｅｌｉｍｉｎａｔｉｏｎ．Ｂｕｔｅａｃｈｍｅｔｈｏｄｈａｓｃｅｒｔａｉｎｌｉｍｉｔｓｉｎｉｔｓａｐｐｌｉｃａｂｉｌｉｔｙ．ＴｈｉｓｐａｐｅｒｍａｉｎｌｙｅｘｐｌｏｒｅｓｔｈｅａｐｐｌｉｃａｔｉｏｎｏｆＲａｄｏｎｔｒａｎｓｆｏｒｍｉｎｔｈｅｍｕｌｔｉｐｌｅ⁃ｗａｖｅａｔｔｅｎｕａｔｉｏｎｏｆｓｈａｌｌｏｗｓｅｉｓｍｉｃｄａｔａ．ＴｈｅｍｏｄｅｌｔｒｉａｌａｎｄａｃｔｕａｌｄａｔａａｐｐｌｉｃａｔｉｏｎｐｒｏｖｅｔｈａｔｔｈｅＲａｄｏｎｆｉｌｔｅｒｉｎｇｍｅｔｈｏｄｗｏｒｋｓｂｅｓｔｗｈｅｎｔｈｅｒｅｓｉｄｕａｌｔｉｍｅｄｉｆｆｅｒｅｎｃｅｂｅｔｗｅｅｎｔｈｅｐｒｉｍａｒｙａｎｄｍｕｌｔｉｐｌｅｗａｖｅｓｉｓｇｒｅａｔｅｒｔｈａｎ３０ｍｓ．Ｋｅｙｗｏｒｄｓ：ｓｈａｌｌｏｗｓｅｉｓｍｉｃ；ｍｕｌｔｉｐｌｅｓ；Ｒａｄｏｎｔｒａｎｓｆｏｒｍ；ｄａｔａｐｒｏｃｅｓｓｉｎｇ；ｍｏｄｅｌｔｅｓｔ；ａｐｐｌｉｃａｔｉｏｎｅｘａｍｐｌｅ。

一种组合式衰减地震资料中多次波的方法

一种组合式衰减地震资料中多次波的方法摘要：在海上地震勘探中我们最常见的干扰波就是多次波。

多次波的存在，严重地干扰了地震记录，妨碍我们对有效波的辨认。

在剖面上存在较强的多次波时，如果在解释中不能正确地把多次波识别出来，就会造成错误的地质解释。

所以，在海洋地震资料的处理中，多次波的衰减是极为重要的一个环节，对于中、深层海洋地震资料，本文采用如下方法压制多次波，先进行SRME 方法去除长周期多次波能量，再应用tau-p 域内预测反褶积去除占大部分的短周期多次波，最后用高精度Rado 变换去除剩余的多次波，实际资料结果表明，该方法有不错的效果。

Abstract: In the offshore seismic exploration, interference wave is the most common multiples wave. The existence of multiple wavesseriously interfere with the seismic record, prevent us identify effective waves. when there is a strong multiple wave in the section, if themultiple wave can not be correctly identified in interpretation, the wrong geological interpretation will caused. Therefore, in dealing withmarine seismic data, the attenuation of multiple wave is an extremely important aspect. For the middle and deep marine seismic data, thispaper uses the following method for multiple wave attenuation, firstly performing SRME to remove the long-period multiple wave energy,then applying the tau-p domain predictive deconvolution to remove most of the short-period multiple wave，finally using high-precisionRado transform to remove the remaining multiples wave. The actual data results show that the method has a good effect.关键词：地震资料处理；多次波压制；tau-p 域预测反褶积；SRME；Rado 变换Key words: seismic data processing；multiple wave suppression；tau-p domain predictive deconvolution；SRME；Rado transformation中图分类号院P631.4 文献标识码院A 文章编号院1006-4311（2014）07-0323-021 多次波类型概述地震记录接受的反射波中，除了接受到有效反射波以外，还包含各种类型的多次波，地下条件越复杂，多次波的类型就越复杂。

基于Radon变换的地震数据多次波衰减方法调研

基于Radon变换的地震数据多次波衰减方法调研姓名：董雅竹学号：SA170070062017.12.06目录背景理论基础方法原理分辨率优化结论背景Radon变换可以有效的实现多次波的衰减。

减轻多次波对一次有效波能量的干扰，增强地震数据的可靠性及地震成像的真实性。

为后续的地震成像、反演解释和钻井定位等提供良好的数据基础。

1234多次波是一种在强反射界面之间进行多次反射的干扰波。

多次波衰减方法的研究经历了从一维到二维、三维空间，从研究简单的层状介质到研究复杂的地下构造，从需要大量的先验信息到不需要先验信息的发展过程。

迄今为止，多次波的衰减方法主要可以分为两大类：一类是滤波法，另一类是波动方程预测相减法（如SRME 法）。

Radon 变换属于滤波法的一种，运算效率高，容易拓展到三维空间，对于远炮检距的地震数据压制效果好。

123理论基础线性Radon变换线性Radon正变换线性Radon反变换（1）（2）（3）（4）线性Radon变换线性Radon正变换示意图抛物线型Radon变换抛物线型Radon正变换抛物线型Radon反变换（1）（3）（2）（4）抛物线型Radon变换抛物线型Radon正变换示意图方法原理抛物线型Radon变换衰减多次波1234再将频率域数据进行Radon变换转换到Radon域，预测出对应的地震波中多次波的部分并去除。

反变换回到时空域，得到衰减多次波后的地震数据。

对动校正后的数据进行傅里叶变换，将其由时空域转化到频率域。

抛物线型Radon变换衰减多次波CMP道集Radon域预测出的多次波压制多次波后的CMP道集在实际中，由于是在有限范围内进行离散运算，进行Radon变换时存在端点效应以及截断效应。

Radon 变换的分辨率问题端点效应是由同相轴的两个端点产生的。

截断效应是由于Radon 变换在沿着不同的积分路径进行积分时，每一个反射同相轴对于所有不同的路径都有贡献而产生的。

1 2 3 这两种效应会严重影响变换的精度，导致变换后在 Radon 域中的能量发散严重，形成“剪刀状尾巴”，分辨率降低。

新疆和田河多次波衰减研究及应用

新疆和田河多次波衰减研究及应用摘要：塔里木盆地和田河东地区受其地质条件特殊，地震剖面中深层多次波发育，严重影响了有效反射波叠加成像，必须最大限度地消除多次波。

关键词：新疆和田河地震资料处理多次波衰减一、前言和田河东区块地处塔克拉玛干沙漠腹地，工区行政区划隶属于新疆维吾尔自治区和田地区。

在构造上位于巴楚隆起东部斜坡的吐木休克断裂带和玛扎塔格断裂带之间。

根据区域及相关钻探资料，该区具有从古生界到中、新生界的巨厚沉积，区自上而下发育第四系、新近系（上新统、中新统）、二叠系、石炭系、泥盆系、志留系、奥陶系、寒武系、震旦系地层等，缺失侏罗系、白垩系，部分缺失古近系，基底为前震旦系变质岩。

从二维地震剖面看，本区地震反射信息较为丰富，主要有T90、T81、T74、T70、T60、T54、T50、T21、T20等反射层。

本区多次波的产生主要是在二叠系顶界反射，一个较强的波阻抗界面，这套强的反射与地表形成全程多次反射，同时，下面地层的反射波与这套地层间又形成层间多次反射。

多次波的存在，造成深层信噪比低，严重影响了有效反射波叠加成像，必须最大限度地消除多次波。

目前消除多次波，尤其是有效消除层间多次波仍然是地球物理界的一大难题。

在项目研究过程中，首先通过多种手段识别多次波，分析其产生机理及特性、分布区域，再根据多次波的特性选择有针对性的技术进行压制。

二、和田河东地区多次波特征和识别如果能够识别出多次波，压制多次波的过程就会有的放。

多次波的识别对地震剖面运用合理有效的多次波分析手段，认识并掌握地震数据中多次波的类型和性质。

根据多次波的特点，在识别和田河多次波时采用以下方法。

1.利用速度谱和扫描速度谱的方法，分析数据中多次波的分布、发育及视速度情况等。

2.利用不同的偏移距叠加剖面分析方法，分析多次波能量分布情况。

3.利用CMP道集识别多次波，在经过动校正后的CMP道集上，有效波被拉平，而多次波由于具有较低的速度校正不足，在剖面上表现为抛物线形态。