CN103308943A - 一种海洋地震资料处理中层间多次波衰减的方法及装置 - Google Patents
一种海洋地震资料处理中层间多次波衰减的方法及装置 Download PDFInfo
- Publication number
- CN103308943A CN103308943A CN2013101724932A CN201310172493A CN103308943A CN 103308943 A CN103308943 A CN 103308943A CN 2013101724932 A CN2013101724932 A CN 2013101724932A CN 201310172493 A CN201310172493 A CN 201310172493A CN 103308943 A CN103308943 A CN 103308943A
- Authority
- CN
- China
- Prior art keywords
- wave field
- interface
- field item
- item
- seismologic record
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012545 processing Methods 0.000 title abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 76
- 239000011229 interlayer Substances 0.000 claims abstract description 18
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000012217 deletion Methods 0.000 claims description 7
- 230000037430 deletion Effects 0.000 claims description 7
- 230000004907 flux Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 abstract 2
- 239000010410 layer Substances 0.000 description 8
- 230000002596 correlated effect Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000011295 pitch Substances 0.000 description 1
- 230000001698 pyrogenic effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种海洋地震资料处理中层间多次波衰减的方法及装置,包括:对原始地震记录进行预处理,得到预处理后的地震记录PPR(x,z0;t);将动校正后的地震记录进行叠加,在叠加后的剖面上获取产生层间多次波的强反射界面L的位置;在所述动校正后的地震记录中拾取界面L产生的一次反射波,并将该一次反射波进行反动校,得到第二波场项P2(x,z0;t);根据所述界面L的位置,删除所述动校正后的地震记录中该界面L以上的反射波场,并进行反动校,得到第一波场项P1(x,z0;t);令第三波场项P3(x,z0;t)等于第一波场项;根据第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w);在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t);由于该方案通过纯数据驱动进行层间多次波预测和衰减,因而不依赖与速度场的精度,提高了海洋地震资料的处理质量。
Description
技术领域
本发明涉及地震资料处理技术领域,特别涉及一种海洋地震资料处理中层间多次波衰减的方法及装置。
背景技术
随着油气勘探的深入,海洋油气勘探受到了越来越多的重视,海洋地震资料处理成果是海洋油气勘探的重要依据,而多次波是影响海洋地震资料处理的最为关键因素之一。由于海平面和海底都存在强速度差的反射界面,反射系数很大,所以表面相关多次波非常严重,严重的影响了地震资料处理的可靠性,降低了数据的信噪比,干扰了地震资料解释人员对构造的判断。在海上地震资料中,地表相关多次波是其主要形式,当前阶段已经有多种地表相关多次波衰减方法,其中应用最为广泛的就是基于反馈环的地表相关多次波衰减(SRME)方法,该方法已经在工业界得到了广泛的应用。
但是,当地下存在较强反射界面时,就会产生能量不容忽视的层间多次波,如海洋资料中海底与海底下盐丘顶部、陆上资料中表层或浅层火成岩之间,均会有层间多次波与有效波相干涉。与地表面相关多次波相比,层间多次波反射路径更为复杂,规律性更弱,因而在剖面上非常不易识别。
为了消除层间多次波对地震资料的影响,当前工业界主要采用基于周期性的预测反褶积和基于多次波与一次波速度差异拉东域多次波衰减方法,但是当多次波周期性不明显或多次波与一次波速度差异小的情况,这类方法处理效果非常有限。目前现有其它较为先进的层间多次波去除方法有王德利等(2011)提出的CFP技术层间多次波预测及小曲线Curvelet域相减方法,该方法的根源要追溯到Berkhout教授最早提出的共聚焦点(CFP)成像技术。该技术具有一定的抗噪能力,但是由于是基于波场延拓的方法,需要具有较高的速度精度。另一类方法是Weglein等人提出完全数据驱动型的逆散射去层间多次波方法,该方法可从逆散射子序列中一次预测出所有的层间多次波,但运行成本非常高,在工业界至今还没有应用。
层间多次波的衰减方法是当前世界性难题,随着油气勘探开发的不断深入,该项技术的研究有着重要的应用价值,对提高地震资料处理品质有着重要意义。
发明内容
为克服现有技术中存在的问题,本发明提供一种海洋地震资料处理中层间多次波衰减的方法及装置,以实现纯数据驱动进行层间多次波预测和衰减,从而不依赖与速度场的精度,提高海洋地震资料的处理质量。
为达到上述目的,本发明一方面提供了一种海洋地震资料处理中层间多次波衰减的方法,包括:
对原始地震记录进行预处理,得到预处理后的地震记录PPR(x,z0;t);
对所述预处理后的地震记录PPR(x,z0;t)进行动校正,得到动校正后的地震记录;
将所述动校正后的地震记录进行叠加,在叠加后的剖面上获取产生层间多次波的强反射界面L的位置;
根据所述界面L的位置,在所述动校正后的地震记录中拾取界面L产生的一次反射波,并将该一次反射波进行反动校,得到第二波场项P2(x,z0;t);
根据所述界面L的位置,删除所述动校正后的地震记录中该界面L以上的反射波场,并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校,得到第一波场项P1(x,z0;t);
将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校,得到第三波场项P3(x,z0;t);
将所述第一波场项P1(x,z0;t)、第二波场项P2(x,z0;t)和第三波场项P3(x,z0;t)进行频域变换,得到频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w);
根据所述频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w);
将所述层间多次波PLM(x,z0;w)进行时域变换,得到时间域的层间多次波PLM(x,z0;t);
在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t)。
本发明另一方面提供了一种海洋地震资料处理中层间多次波衰减的装置,包括:
预处理单元,用于对原始地震记录进行预处理,得到预处理后的地震记录PPR(x,z0;t);
动校正单元,用于对所述预处理后的地震记录PPR(x,z0;t)进行动校正,得到动校正后的地震记录;
界面L拾取单元,用于将所述动校正后的地震记录进行叠加,在叠加后的剖面上获取产生层间多次波的强反射界面L的位置;
第二波场项计算单元,用于根据所述界面L的位置,在所述动校正后的地震记录中拾取界面L产生的一次反射波,并将该一次反射波进行反动校,得到第二波场项P2(x,z0;t);
第一波场项计算单元,用于根据所述界面L的位置,删除所述动校正后的地震记录中该界面L以上的反射波场,并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校,得到第一波场项P1(x,z0;t);
第三波场项计算单元,用于将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校,得到第三波场项P3(x,z0;t);
频域变换单元,用于将所述第一波场项P1(x,z0;t)、第二波场项P2(x,z0;t)和第三波场项P3(x,z0;t)进行频域变换,得到频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w);
层间多次波预测单元,用于根据所述频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w);
时域变换单元,用于将所述层间多次波PLM(x,z0;w)进行时域变换,得到时间域的层间多次波PLM(x,z0;t);
多次波衰减单元,用于在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t)。
本发明能够达到的有益效果:通过纯数据驱动型的层间多次波预测和衰减,能够不依赖与速度场的精度,提高了海洋地震资料的处理质量。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明的原理示意图;
图2是具体实施方式的简单模型速度场;
图3中,a是包含层间多次波的预处理后的地震记录;b是产生层间多次波高速层的一次反射波,也就是第二波场项;c是产生层间多次波高速层以下层界面的反射波场,是第一波场项,也近似为第三波场项;
图4中,a是包含层间多次波的原始单炮记录;b是匹配衰减完层间多次波后的记录;c是预测出的层间多次波;
图5是本发明实施例一种海洋地震资料处理中层间多次波衰减的方法的流程图;
图6是本发明实施例一种海洋地震资料处理中层间多次波衰减的装置的结构图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
图1是本发明的简单原理图,地表有实际激发震源S,接收记录R,地下存在强反射界面L。层间多次波传播路径是SO1OO2R,S1与R1分别是虚震源和虚接收点。将波场分为三个波场项来构建该层间多次波传播路径,第一波场项传播路径是SO1R1,第二波场项S1OR1和第三波场项S1O2R,这三个波场项都是可以在地表接收到的,层间多次波SO1OO2R可以通过第一波场项SO1R1与第三波场项S1O2R的褶积再与第二波场项S1OR1的相关来得到。
基于强反射界面间会产生可识别的层间多次波这一认识,建立如图2所示的理论模型,该模型第二层为一高速层,速度为4000m/s,层间多次波主要由该层顶界面与其下反射层产生,采用二阶声波有限差分模拟合成记录,子波主频20Hz,每炮101道接收,道距20米,中间放炮,炮距20米,共模拟单炮101炮,如图3中,(a)为2000米处激发的单炮合成记录,其中,M1,M2和M3均为与该强反射层相关的层间多次波。为突出层间多次波,模拟时在表面加吸收边界,该记录中不存在自由表面多次波。按照本发明的方法原理,利用(1)式预测高速层顶界面产生的层间多次波,需要从去掉干扰波后的数据中提取产生层间多次波界面的一次反射波也就是第二波场项,以及该界面以下的所有反射波即第一波场项。为防止直接切除损伤部分有效波,先将炮集记录抽成CMP道集,动校正后选择矩形窗切除,然后反动校抽成炮集,即可获得所需要的两类波场,图3中,(b)是层间多次波相关界面的一次反射波P2(x,z0;t),将该界面以上所有反射切除后的波场为第二波场项P1(x,z0;t),即图3中的(c),同时也作为第三波场项P3(x,z0;t)。
图4中的(a)是根据式(2)相减后得到的衰减层间多次波后的记录PME(x,z0;t),图4中的(b)是根据式(1)预测出的与4000m/s速度层相关的层间多次波,从图中可以看出,层间多次波M1,M2和M3均得到了有效压制。
具体通过以下实施例进行说明:
实施例一
如图5所示,为本实施例一种海洋地震资料处理中层间多次波衰减的方法的流程图,包括以下步骤:
步骤501,对原始地震记录进行预处理,得到预处理后的地震记录PPR(x,z0;t)。
对原始地震记录进行预处理,包括面波去除、随机噪音衰减、野值去除和地表相关多次波衰减,得到预处理后的地震记录PPR(x,z0;t)。
步骤502,对所述预处理后的地震记录PPR(x,z0;t)进行动校正,得到动校正后的地震记录;
将整个炮集动校正的目的是为了更容易的拾取处界面L的一次反射,拾取出来后再反动校回去就相当于从炮集数据中提取出了L层的一次反射。
步骤503,将所述动校正后的地震记录进行叠加,在叠加后的剖面上获取产生层间多次波的强反射界面L的位置。
步骤504,根据所述界面L的位置,在所述动校正后的地震记录中拾取界面L产生的一次反射波,并将该一次反射波进行反动校,得到第二波场项P2(x,z0;t)。
步骤505,根据所述界面L的位置,删除所述动校正后的地震记录中该界面L以上的反射波场,并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校,得到第一波场项P1(x,z0;t)。
步骤506,将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校,得到第三波场项P3(x,z0;t)。
步骤507,将所述第一波场项P1(x,z0;t)、第二波场项P2(x,z0;t)和第三波场项P3(x,z0;t)进行频域变换,得到频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)。
步骤508,根据所述频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w)。
基于地震干涉测量理论,在波场运动学中,褶积过程可以认为是两个波场传播路径的求和,而相关是两个传播路径的相减。在时间域的褶积等价于在频率域的乘积,在时间域的相关等价于在频率域一个信号与另外一个信号的复共轭的乘积。在频率域,定义界面L以下层的所有反射波场为第一波场项P1(x,z0;w),定义界面L的一次反射波场为第二波场项P2(x,z0;w),实际生产中无法确定界面L以下所有构造的一次反射P3(x,z0;w),所以近似第三波场项P3(x,z0;w)等于第一波场项P1(x,z0;w)。在频率域,第一波场项与第二波场项复共轭的乘积再与第三波场项的乘积就构建了界面L相关的层间多次波场PLM(x,z0;w)。公式表示为:
其中,x是空间横向位置坐标,z0代表波场在地表接收,w代表频率域,PLM(x,z0;w)是预测出的与界面L相关的层间多次波。
步骤509,将所述层间多次波PLM(x,z0;w)进行时域变换,得到时间域的层间多次波PLM(x,z0;t)。
步骤510,在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t)。
预测的层间多次波PLM(x,z0;t)与炮集数据PPR(x,z0;t)需要进行能量匹配,在相同的数量级后才能进行匹配相减,基于最小二乘匹配滤波方法,得到匹配因子,然后从炮集数据Ppr(x,z0;t)中减去匹配滤波后的层间多次波,即得到衰减掉的与L层相关层间多次波记录,数学描述如下:
其中,t是记录时间,A(x,τ)是匹配因子,NT是整个记录时间的长度。
实施例二
如图6所示,为本实施例一种海洋地震资料处理中层间多次波衰减的装置的结构图,包括以下单元:
预处理单元601,用于对原始地震记录进行预处理,得到预处理后的地震记录PPR(x,z0;t);
对原始地震记录进行预处理,包括面波去除、随机噪音衰减、野值去除和地表相关多次波衰减,得到预处理后的地震记录PPR(x,z0;t)。
动校正单元602,用于对所述预处理后的地震记录PPR(x,z0;t)进行动校正,得到动校正后的地震记录;
将整个炮集动校正的目的是为了更容易的拾取处界面L的一次反射,拾取出来后再反动校回去就相当于从炮集数据中提取出了L层的一次反射。
界面L拾取单元603,用于将所述动校正后的地震记录进行叠加,在叠加后的剖面上获取产生层间多次波的强反射界面L的位置;
第二波场项计算单元604,用于根据所述界面L的位置,在所述动校正后的地震记录中拾取界面L产生的一次反射波,并将该一次反射波进行反动校,得到第二波场项P2(x,z0;t);
第一波场项计算单元605,用于根据所述界面L的位置,删除所述动校正后的地震记录中该界面L以上的反射波场,并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校,得到第一波场项P1(x,z0;t);
第三波场项计算单元606,用于将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校,得到第三波场项P3(x,z0;t);
频域变换单元607,用于将所述第一波场项P1(x,z0;t)、第二波场项P2(x,z0;t)和第三波场项P3(x,z0;t)进行频域变换,得到频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w);
层间多次波预测单元608,用于根据所述频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w);
基于地震干涉测量理论,在波场运动学中,褶积过程可以认为是两个波场传播路径的求和,而相关是两个传播路径的相减。在时间域的褶积等价于在频率域的乘积,在时间域的相关等价于在频率域一个信号与另外一个信号的复共轭的乘积。在频率域,定义界面L以下层的所有反射波场为第一波场项P1(x,z0;w),定义界面L的一次反射波场为第二波场项P2(x,z0;w),实际生产中无法确定界面L以下所有构造的一次反射P3(x,z0;w),所以近似第三波场项P3(x,z0;w)等于第一波场项P1(x,z0;w)。在频率域,第一波场项与第二波场项复共轭的乘积再与第三波场项的乘积就构建了界面L相关的层间多次波场PLM(x,z0;w)。公式表示为:
其中,x是空间横向位置坐标,z0代表波场在地表接收,w代表频率域,PLM(x,z0;w)是预测出的与界面L相关的层间多次波。
时域变换单元609,用于将所述层间多次波PLM(x,z0;w)进行时域变换,得到时间域的层间多次波PLM(x,z0;t);
多次波衰减单元610,用于在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t)。
预测的层间多次波PLM(x,z0;t)与炮集数据PPR(x,z0;t)需要进行能量匹配,在相同的数量级后才能进行匹配相减,基于最小二乘匹配滤波方法,得到匹配因子,然后从炮集数据Ppr(x,z0;t)中减去匹配滤波后的层间多次波,即得到衰减掉的与L层相关层间多次波记录,数学描述如下:
其中,t是记录时间,A(x,τ)是匹配因子,NT是整个记录时间的长度。
本发明的有益效果是:该方法是纯数据驱动型,不依赖与速度场的精度,是一种具有较高应用价值的层间多次波衰减方法。
上述实施例能够达到的有益效果:通过纯数据驱动型的层间多次波预测和衰减,能够不依赖与速度场的精度,提高了海洋地震资料的处理质量。
本领域一般技术人员在此设计思想之下所做任何不具有创造性的改造,均应视为在本发明的保护范围之内。
Claims (8)
1.一种海洋地震资料处理中层间多次波衰减的方法,其特征在于,包括:
对原始地震记录进行预处理,得到预处理后的地震记录PPR(x,z0;t);
对所述预处理后的地震记录PPR(x,z0;t)进行动校正,得到动校正后的地震记录;
将所述动校正后的地震记录进行叠加,在叠加后的剖面上获取产生层间多次波的强反射界面L的位置;
根据所述界面L的位置,在所述动校正后的地震记录中拾取界面L产生的一次反射波,并将该一次反射波进行反动校,得到第二波场项P2(x,z0;t);
根据所述界面L的位置,删除所述动校正后的地震记录中该界面L以上的反射波场,并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校,得到第一波场项P1(x,z0;t);
将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校,得到第三波场项P3(x,z0;t);
将所述第一波场项P1(x,z0;t)、第二波场项P2(x,z0;t)和第三波场项P3(x,z0;t)进行频域变换,得到频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w);
根据所述频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w);
将所述层间多次波PLM(x,z0;w)进行时域变换,得到时间域的层间多次波PLM(x,z0;t);
在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t)。
2.如权利要求1所述海洋地震资料处理中层间多次波衰减的方法,其特征在于,所述根据所述频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w)具体为:
其中,x是空间横向位置坐标,z0代表波场在地表接收,w代表频率域。
3.如权利要求1所述海洋地震资料处理中层间多次波衰减的方法,其特征在于,所述将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校,得到第三波场项P3(x,z0;t)的算法为:
P3(x,z0;t)=P1(x,z0;t)。
4.如权利要求1所述海洋地震资料处理中层间多次波衰减的方法,其特征在于,所述在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t)具体为:
根据最小二乘法在所述预处理后的地震记录PPR(x,z0;t)中获取能量匹配因子A(x,τ);
根据所述匹配因子对所述层间多次波PLM(x,z0;t)进行能量匹配,并在所述预处理后的地震记录PPR(x,z0;t)中衰减掉进行能量匹配后的层间多次波,得到衰减掉层间多次波的记录PME(x,z0;t),具体为:
其中,t是记录时间,NT是整个记录时间的长度。
5.一种海洋地震资料处理中层间多次波衰减的装置,其特征在于,包括:
预处理单元,用于对原始地震记录进行预处理,得到预处理后的地震记录PPR(x,z0;t);
动校正单元,用于对所述预处理后的地震记录PPR(x,z0;t)进行动校正,得到动校正后的地震记录;
界面L拾取单元,用于将所述动校正后的地震记录进行叠加,在叠加后的剖面上获取产生层间多次波的强反射界面L的位置;
第二波场项计算单元,用于根据所述界面L的位置,在所述动校正后的地震记录中拾取界面L产生的一次反射波,并将该一次反射波进行反动校,得到第二波场项P2(x,z0;t);
第一波场项计算单元,用于根据所述界面L的位置,删除所述动校正后的地震记录中该界面L以上的反射波场,并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校,得到第一波场项P1(x,z0;t);
第三波场项计算单元,用于将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校,得到第三波场项P3(x,z0;t);
频域变换单元,用于将所述第一波场项P1(x,z0;t)、第二波场项P2(x,z0;t)和第三波场项P3(x,z0;t)进行频域变换,得到频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w);
层间多次波预测单元,用于根据所述频率域的第一波场项P1(x,z0;w)、第二波场项P2(x,z0;w)和第三波场项P3(x,z0;w)计算得到与所述界面L相关的层间多次波PLM(x,z0;w);
时域变换单元,用于将所述层间多次波PLM(x,z0;w)进行时域变换,得到时间域的层间多次波PLM(x,z0;t);
多次波衰减单元,用于在所述预处理后的地震记录PPR(x,z0;t)中衰减掉所述时间域的层间多次波PLM(x,z0;t),得到衰减掉层间多次波的记录PME(x,z0;t)。
6.如权利要求5所述海洋地震资料处理中层间多次波衰减的装置,其特征在于,所述层间多次波预测单元根据以下公式计算层间多次波PLM(x,z0;w):
其中,x是空间横向位置坐标,z0代表波场在地表接收,w代表频率域。
7.如权利要求5所述海洋地震资料处理中层间多次波衰减的装置,其特征在于,所述第三波场项计算单元计算第三波场项P3(x,z0;t)的算法为:
P3(x,z0;t)=P1(x,z0;t)。
8.如权利要求5所述海洋地震资料处理中层间多次波衰减的装置,其特征在于,所述多次波衰减单元根据以下公式进行衰减:
其中,t是记录时间,NT是整个记录时间的长度,A(x,τ)是根据最小二乘法在所述预处理后的地震记录PPR(x,z0;t)中获取的能量匹配因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310172493.2A CN103308943B (zh) | 2013-05-10 | 2013-05-10 | 一种海洋地震资料处理中层间多次波衰减的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310172493.2A CN103308943B (zh) | 2013-05-10 | 2013-05-10 | 一种海洋地震资料处理中层间多次波衰减的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103308943A true CN103308943A (zh) | 2013-09-18 |
CN103308943B CN103308943B (zh) | 2015-11-18 |
Family
ID=49134363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310172493.2A Active CN103308943B (zh) | 2013-05-10 | 2013-05-10 | 一种海洋地震资料处理中层间多次波衰减的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103308943B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984023A (zh) * | 2014-05-15 | 2014-08-13 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN104977615A (zh) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | 一种基于模型统计拾取的深水obc资料多次波压制方法 |
CN105259574A (zh) * | 2015-10-12 | 2016-01-20 | 中国石油大学(华东) | 基于一次波稀疏约束的多道预测反褶积方法 |
CN105572735A (zh) * | 2014-10-17 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种提高叠前深度偏移成像精度的方法及装置 |
CN106199699A (zh) * | 2016-06-30 | 2016-12-07 | 恒泰艾普(北京)能源科技研究院有限公司 | 利用传输矩阵法清除鬼波的方法 |
CN106324669A (zh) * | 2016-08-22 | 2017-01-11 | 中国石油大学(华东) | 一种地震勘探数据中各阶表层多次波分离方法 |
CN104360394B (zh) * | 2014-12-01 | 2017-02-22 | 中国海洋石油总公司 | 一种海上变深度缆地震数据鬼波压制的方法 |
CN108196304A (zh) * | 2018-01-09 | 2018-06-22 | 中国石油天然气股份有限公司 | 一种多次波压制方法及装置 |
CN109425897A (zh) * | 2017-08-30 | 2019-03-05 | 中国石油化工股份有限公司 | 消除地震资料野值干扰的方法及*** |
CN110873893A (zh) * | 2018-09-03 | 2020-03-10 | 中国石油化工股份有限公司 | 深水自由表面多次波预测和压制方法及其*** |
CN111308554A (zh) * | 2020-03-10 | 2020-06-19 | 清华大学 | 一种基于多次波生成层自适应提取的层间多次波预测方法 |
CN111551993A (zh) * | 2019-02-12 | 2020-08-18 | 中国石油天然气股份有限公司 | 压制鸣震的方法及装置 |
CN117214950A (zh) * | 2023-09-25 | 2023-12-12 | 中海石油(中国)有限公司上海分公司 | 一种多次波压制方法、装置、设备及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6101448A (en) * | 1998-01-15 | 2000-08-08 | Schlumberger Technology Corporation | Multiple attenuation of multi-component sea-bottom data |
US6219621B1 (en) * | 1998-06-30 | 2001-04-17 | Exxonmobil Upstream Research Co. | Sparse hyperbolic inversion of seismic data |
US6832161B1 (en) * | 2003-08-15 | 2004-12-14 | Westerngeco, L.L.C. | Method for attenuating water layer multiples |
CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
CN102169189A (zh) * | 2011-01-19 | 2011-08-31 | 中国海洋石油总公司 | 深水层间多次波消除方法 |
CN102636811A (zh) * | 2012-04-10 | 2012-08-15 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种海上二维地震资料中多次波的消除方法 |
CN102741710A (zh) * | 2009-11-20 | 2012-10-17 | 雪佛龙美国公司 | 考虑波束方位的用于衰减地震数据中的多次波的***和方法 |
-
2013
- 2013-05-10 CN CN201310172493.2A patent/CN103308943B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6101448A (en) * | 1998-01-15 | 2000-08-08 | Schlumberger Technology Corporation | Multiple attenuation of multi-component sea-bottom data |
US6219621B1 (en) * | 1998-06-30 | 2001-04-17 | Exxonmobil Upstream Research Co. | Sparse hyperbolic inversion of seismic data |
US6832161B1 (en) * | 2003-08-15 | 2004-12-14 | Westerngeco, L.L.C. | Method for attenuating water layer multiples |
RU2006107997A (ru) * | 2003-08-15 | 2006-07-10 | Вестернджеко, Л.Л.С. (Us) | Способ ослабления многократных волн от границ водного слоя |
CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
CN102741710A (zh) * | 2009-11-20 | 2012-10-17 | 雪佛龙美国公司 | 考虑波束方位的用于衰减地震数据中的多次波的***和方法 |
CN102169189A (zh) * | 2011-01-19 | 2011-08-31 | 中国海洋石油总公司 | 深水层间多次波消除方法 |
CN102636811A (zh) * | 2012-04-10 | 2012-08-15 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种海上二维地震资料中多次波的消除方法 |
Non-Patent Citations (3)
Title |
---|
李鹏等: "多次波问题的研究进展", 《地球物理学进展》 * |
牛滨华等: "波动方程压制多次波的技术方法", 《地学前缘》 * |
黄新武等: "基于数据一致性原理预测与压制自由表面多次波的效果分析", 《石油地球物理勘探》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984023A (zh) * | 2014-05-15 | 2014-08-13 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN103984023B (zh) * | 2014-05-15 | 2017-01-25 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN105572735A (zh) * | 2014-10-17 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种提高叠前深度偏移成像精度的方法及装置 |
CN105572735B (zh) * | 2014-10-17 | 2018-02-02 | 中国石油天然气股份有限公司 | 一种提高叠前深度偏移成像精度的方法及装置 |
CN104360394B (zh) * | 2014-12-01 | 2017-02-22 | 中国海洋石油总公司 | 一种海上变深度缆地震数据鬼波压制的方法 |
CN104977615A (zh) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | 一种基于模型统计拾取的深水obc资料多次波压制方法 |
CN105259574A (zh) * | 2015-10-12 | 2016-01-20 | 中国石油大学(华东) | 基于一次波稀疏约束的多道预测反褶积方法 |
CN106199699A (zh) * | 2016-06-30 | 2016-12-07 | 恒泰艾普(北京)能源科技研究院有限公司 | 利用传输矩阵法清除鬼波的方法 |
CN106324669A (zh) * | 2016-08-22 | 2017-01-11 | 中国石油大学(华东) | 一种地震勘探数据中各阶表层多次波分离方法 |
CN109425897A (zh) * | 2017-08-30 | 2019-03-05 | 中国石油化工股份有限公司 | 消除地震资料野值干扰的方法及*** |
CN109425897B (zh) * | 2017-08-30 | 2020-04-07 | 中国石油化工股份有限公司 | 消除地震资料野值干扰的方法及*** |
CN108196304A (zh) * | 2018-01-09 | 2018-06-22 | 中国石油天然气股份有限公司 | 一种多次波压制方法及装置 |
CN108196304B (zh) * | 2018-01-09 | 2019-07-09 | 中国石油天然气股份有限公司 | 一种多次波压制方法及装置 |
CN110873893A (zh) * | 2018-09-03 | 2020-03-10 | 中国石油化工股份有限公司 | 深水自由表面多次波预测和压制方法及其*** |
CN111551993A (zh) * | 2019-02-12 | 2020-08-18 | 中国石油天然气股份有限公司 | 压制鸣震的方法及装置 |
CN111308554A (zh) * | 2020-03-10 | 2020-06-19 | 清华大学 | 一种基于多次波生成层自适应提取的层间多次波预测方法 |
CN117214950A (zh) * | 2023-09-25 | 2023-12-12 | 中海石油(中国)有限公司上海分公司 | 一种多次波压制方法、装置、设备及存储介质 |
CN117214950B (zh) * | 2023-09-25 | 2024-05-28 | 中海石油(中国)有限公司上海分公司 | 一种多次波压制方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103308943B (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103308943B (zh) | 一种海洋地震资料处理中层间多次波衰减的方法及装置 | |
US9075163B2 (en) | Interferometric seismic data processing | |
CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
CN102590859B (zh) | 垂向各向异性介质准p波方程逆时偏移方法 | |
CN106932824B (zh) | 陆地地震勘探资料的降维自适应层间多次波压制方法 | |
CN101598809A (zh) | 一种自适应消除线性规则噪声以及多次波干扰的方法 | |
CN103424777B (zh) | 一种提高地震成像分辨率的方法 | |
Williams et al. | Scholte wave inversion and passive source imaging with ocean-bottom DAS | |
CN101620276A (zh) | 用于在地震数据中衰减多次反射的方法 | |
EP2622376A2 (en) | Interferometric seismic data processing for a towed marine survey | |
CN104755963A (zh) | 用于衰减地震数据中的噪声的***和方法 | |
WO2012021218A2 (en) | Attenuating internal multiples from seismic data | |
CN104345343B (zh) | 一种复杂海底相关的层间多次波预测方法 | |
US20110199858A1 (en) | Estimating internal multiples in seismic data | |
Mutter et al. | The Mohorovicic discontinuity in ocean basins: Some observations from seismic data | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN103558633A (zh) | 一种基于稀疏反演的层间多次波压制方法 | |
CN104977615B (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN104570095A (zh) | 一种基于Radon变换消除斜缆虚反射的方法 | |
CN102269824B (zh) | 一种地震数据子波相位转换处理的方法 | |
CN108957522B (zh) | 一种基于克希霍夫积分的海底地震仪数据多次波预测方法 | |
Yang et al. | Crosswell frequency-domain reverse time migration imaging with wavefield decomposition | |
EP3387466B1 (en) | Efficient internal multiple prediction methods | |
US10564305B2 (en) | Efficient internal multiple prediction methods | |
CN115236730A (zh) | 一种层间多次波傅里叶有限差分的地震波场偏移成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |