一种地震噪声去除方法及装置.pdf

《一种地震噪声去除方法及装置.pdf》由会员分享，可在线阅读，更多相关《一种地震噪声去除方法及装置.pdf（18页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103984023 A (43)申请公布日 2014.08.13 CN 103984023 A (21)申请号 201410205484.3 (22)申请日 2014.05.15 G01V 1/36(2006.01) G06F 19/00(2011.01) (71)申请人 中国科学院地质与地球物理研究所 地址 100029 北京市朝阳区北土城西路 19 号 申请人 中国地质大学 （北京） (72)发明人 史小东 刘洪 孟小红 (74)专利代理机构 北京三友知识产权代理有限 公司 11127 代理人 贾磊 (54) 发明名称 一种地震噪声去除方法及装置 (57) 摘要 。

2、本发明属于地震资料处理领域， 具体涉及一 种地震噪声去除方法及装置。其中通过针对层状 介质层间多次波预测算子的快速计算， 得到层间 多次波预测的炮集数据。基于预测算子的层间多 次波预测技术可以一次预测所有层相关的层间多 次波， 利用层状介质的层间多次波和下覆层介质 在层状介质中产生的层间多次波计算得到地震噪 声， 并在炮集数据中减去所计算的层间多次波炮 集数据从而得到纯净的炮集数据。本发明的有益 效果在于能够大幅度的减小计算代价。 (51)Int.Cl. 权利要求书 2 页 说明书 8 页 附图 7 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书8页 附。

3、图7页 (10)申请公布号 CN 103984023 A CN 103984023 A 1/2 页 2 1. 一种地震噪声去除方法， 其特征在于包括， 步骤 101， 对炮集数据进行分离， 分为层状介质数据 Slayer和下覆层数据 Syu； 步骤 102， 对炮集数据中的层状介质数据 Slayer进行 -p 变换得到 STlayer， 利用 计算 STlayer每个 P 值对应的地震道的第 n 个地震数据的 预测算子， 其中， An1、 An2分别为第 n 个地震数据的预测算子的实部和虚部， fr、 fi 分别为 P 值对应的地震道中第 1 到第 n 个地震数据的傅里叶变换的实部和虚部， n。

4、 取 1 至 N 的正整 数， N 为 STlayer在时间方向地震数据的个数 ； 步骤 103， 利用 Anmultiple An(n)-An-1(n-1) 求取第 n 个地震数据的多次波预测算子 Anmultiple， 对 Anmultiple做反傅里叶变换， 得到时间域多次波预测算子 anmultiple， 令 anmultiple(n) 0, 得到 anint； 其中， anmultiple为第 n 个地震数据的多次波预测算子 ,anint为第 n 个地震数据 的层间多次波预测算子 ； 步骤 104， 按照公式预测所有 Slayer相关的 -p 域层间多次波 ； 其 中， di为第 i。

5、 个地震数据 , 为子波长度， * 表示数学上的褶积运算 ； 步骤 105， 计算层状介质数据中所有的层间多次波预测算子序列并求和记 为alayer， 对下覆层数据Syu进行-p变换得到STyu， 按照公式MyuSTyu*alayer得到下覆层介 质在层状介质中产生的层间多次波 Myu； 步骤 106， 将 Mlayer和 Myu进行相加， 得到 -p 域预测的层间多次波炮集数据， 然后进行 -p 反变换， 得到 X-T 域层间多次波的炮集数据 M ； 步骤 107， 在炮集数据中减去计算得到的层间多次波 M， 得到不含层间多次波的炮集数 据。 2. 根据权利要求 1 所述的一种地震噪声去除方。

6、法， 其特征在于， 在上述对炮集数据进 行分离之前还包括， 对所述炮集数据进行偏移， 确定层状介质的范围和状态， 若所述层状介 质状态为倾斜， 则将炮点和检波点延拓至和地层平行的倾斜面上。 3. 根据权利要求 1 所述的一种地震噪声去除方法， 其特征在于， 在上述在炮集数据中 减去层间多次波之前还包括， 对每一炮的炮集数据进行上述步骤 101- 步骤 106， 得到全部 炮集数据对应的层间多次波炮集数据。 4. 一种地震噪声去除装置， 其特征在于包括， 分离单元， 用于对炮集数据进行分离， 分为层状介质炮集数据 Slayer和下覆层炮集数据 Syu； 预测算子计算单元， 用于对炮集数据中的层状。

7、介质数据进行 -p 变换， 利用 权 利 要 求 书 CN 103984023 A 2 2/2 页 3 计算每个 P 值对应的地震道的第 n 个地震数据的预测算 子， 其中， An1、 An2分别为第 n 个地震数据对应的预测算子的实部和虚部， fr、 fi 分别为每道 地震数据中第 1 到第 n 个地震数据的傅里叶变换的实部和虚部， n 取 1 至 N 的正整数， N 为 Slayer在时间方向地震数据的个数 ； 层间多次波预测算子计算单元， 利用 Anmultiple An(n)-An-1(n-1) 求取第 n 个地震 数据的多次波预测算子 Anmultiple， 对 Anmultiple。

8、做反傅里叶变换， 得到时间域多次波预测算子 anmultiple， 令 anmultiple(n) 0, 得到 anint； 其中， anmultiple为第 n 个地震数据的多次波预测算 子 ,anint为第 n 个地震数据的层间多次波预测算子 ； 层状介质层间多次波预测单元， 按照公式预测所有层状介质相关 的层间多次波 Mlayer， 其中， di为第 i 个地震数据， 为子波长度， * 表示数学上的褶积运算 ； 下覆层层间多次波预测单元， 用于计算所述层状介质数据中所有的层间多次波预测算 子的序列并求和记为 alayer， 对该炮集数据的下覆层数据 Syu进行 -p 变换得到 STyu，。

9、 按照公式 Myu STyu*alayer得到下覆层介质在层状介质中产生的层间多次波 Myu； 层间多次波计算单元， 用于将该炮集数据的 Mlayer和 Myu进行相加， 得到 -p 域预测的 层间多次波炮集数据， 并进行 -p 反变换， 得到 X-T 域的层间多次波炮集数据 ； 解扰单元， 用于在炮集数据中减去所述层间多次波炮集数据。 5. 根据权利要求 4 所述的一种地震噪声去除装置， 其特征在于， 还包括偏移单元， 与所 述分离单元相连接， 对所述炮集数据进行偏移， 确定层状介质的范围和状态， 若所述层状介 质状态为倾斜， 则将炮点和检波点延拓至和地层平行的倾斜面上。 6. 根据权利要求。

10、 4 所述的一种地震噪声去除装置， 其特征在于， 还包括循环单元， 与所 述解扰单元相连接， 调用分离单元、 预测算子计算单元、 层间多次波预测算子计算单元、 层 状介质层间多次波预测单元、 下覆层层间多次波预测单元和层间多次波计算单元计算每一 炮中层间多次波炮集数据。 权 利 要 求 书 CN 103984023 A 3 1/8 页 4 一种地震噪声去除方法及装置 技术领域 0001 本发明属于地震资料处理领域， 具体属于一种地震噪声去除方法及装置。 背景技术 0002 地震是油气勘探最主要的手段， 在进行地震勘探时， 地震波在不同地层之间、 薄层 之间、 裂缝内部与裂缝之间都会发生多次散射。

11、或折射， 然后被检波器接受， 这类的波统称为 多次波或多次散射。 多次波处理不当容易混淆、 干扰一次波的能量， 甚至出现多次波同相轴 掩盖了一次波同相轴的情况， 在很大程度上降低了地震资料的偏移成像效果， 甚至导致对 地震资料解释和地下构造的错误认识， 直接影响着探井井位部署和勘探的成功率。 因此， 如 何有效地解决叠前多次波的压制问题是地震资料处理中的重点、 难点和热点。 0003 一般的， 地震波在地下介质传播过程中， 如果只有一次上行反射然后被接受点接 受， 称为一次波 ； 如果发生多次下行反射然后被接受点接受， 称为多次波。根据多次波下行 反射发生的位置差异， 多次波可分为表面多次波和。

12、层间多次波。自由表面多次波是指在自 由表面至少发生一次下行反射的地震波， 层间多次波是指下行反射发生在自由表面以下的 反射界面的地震波。 当前， 表面多次波预测与衰减已有成熟的理论和方法， 且得到了大规模 的应用。而层间多次波的消除技术还不完善， 仍然面临着巨大的挑战。 0004 目前层间多次波衰减技术主要包括三类 ： 0005 (1) 基于时差差异去除层间多次波 0006 基于时差和倾角的差异进行滤波是一类重要的多次波去除方法。其原理是对 CMP( 共中心点道集 ) 道集进行 NMO( 正常时差校正 )， 校正后的道集中， 一次波被校平， 而多 次波未被校平。 实际资料处理中运用较多的Rad。

13、on(拉冬)变换+切除即属于此类方法。 一 般情况下， 基于时差的滤波方法消除由浅层地层产生的层间多次波效果不错， 而对于深层 层间多次波， 区别一次波和层间多次波更为复杂， 因此， 该方法的去除效果很大程度上依赖 处理员的解释。 0007 (2) 扩展 SRME( 表面相关的多次波衰减 ) 方法 0008 实际地震剖面上， 只能看出有限的地层产生的层间多次波， 因此， 一个合理的去除 层间多次波的扩展 SRME 算法就是利用全波场重建基准面方法把所有的炮点和检波点都延 拓至产生层间多次波的界面上， 于是地下的反射面变成了新的 “表面” ， 这样就可以重复利 用 SRME 算法去除层间多次波。。

14、扩展 SRME 的方法有模型驱动和数据驱动两种办法。 0009 (3) 逆散射级数法层间多次波去除技术 0010 Coates 和 Weglein(1996),Weglein 等 (1997) 都对逆散射级数层间多次波去除技 术进行过阐述， 它是一种可预测出所有层间多次波的方法。该预测方法来源于逆散射子级 数， 可解释为对所有可能成像的三个数据点进行求和， 该三个数据点具有 “低高低” 的 约束关系。该算法预测一阶所有层间多次波公式如下 ： 说 明 书 CN 103984023 A 4 2/8 页 5 0011 0012 对于多个层状介质而言， CFP( 共聚焦点预测 ) 算法很难确定产生层间。

15、多次波的界 面， 因而不适合处理多个层状介质产生的层间多次波。而逆散射技术在二维的情况下层间 多次波预测公式需要 5 重积分， 计算量非常大， 计算成本非常高， 因而很难应用于实际。 发明内容 0013 本发明的目的在于解决现有技术中去除地震噪声中的计算复杂， 计算成本高的问 题， 提出了一种地震噪声去除方法及装置， 高效的计算层状介质产生的层间多次波， 用以去 除地震噪声。 0014 一种地震噪声去除方法， 包括， 0015 步骤 101， 对炮集数据进行分离， 分为层状介质数据 Slayer和下覆层数据 Syu； 0016 步骤 102， 对炮集数据中的层状介质数据 Slayer进行 -p。

16、( 倾斜叠加 ) 变换得到 STlayer， 利用计算 STlayer每个 P 值对应的地震道的第 n 个 地震数据的预测算子， 其中， An1、 An2分别为第 n 个地震数据的预测算子的实部和虚部， fr、 fi 分别为 P 值对应的地震道中第 1 到第 n 个地震数据的傅里叶变换的实部和虚部， n 取 1 至 N 的正整数， N 为 STlayer在时间方向地震数据的个数 ； 0017 步骤103， 利用AnmultipleAn(n)-An-1(n-1)求取第n个地震数据的多次波预测算 子 Anmultiple， 对 Anmultiple做反傅里叶变换， 得到时间域多次波预测算子 anm。

17、ultiple， 令 anmultiple(n) 0, 得到 anint； 其中， anmultiple为第 n 个地震数据的多次波预测算子 ,anint为第 n 个地震数 据的层间多次波预测算子 ； 0018 步骤 104， 按照公式预测所有 STlayer相关的 -p 域层间多次 波 ； 其中， di为第 i 个地震数据 , 为子波长度， 一般等于 25 ；* 表示数学上的褶积运算。 0019 步骤 105， 计算层状介质数据中所有的层间多次波预测算子序列并求和 记为 alayer， 对下覆层数据 Syu进行 -p( 倾斜叠加 ) 变换得到 STyu， 按照公式 Myu STyu*alay。

18、er得到下覆层介质在层状介质中产生的层间多次波 Myu； 0020 步骤 106， 将 Mlayer和 Myu进行相加， 得到 -p 域预测的层间多次波， 然后进行 -p 反变换， 得到 X-T 域 ( 时间偏移距域 ) 层间多次波的炮集数据 M ； 说 明 书 CN 103984023 A 5 3/8 页 6 0021 步骤 107， 在炮集数据中减去计算得到的层间多次波 M， 得到不含层间多次波的炮 集数据。 0022 根据本发明实施例所述一种地震噪声去除方法的一个进一步的方面， 在上述对炮 集数据进行分离之前还包括， 对所述炮集数据进行偏移， 确定层状介质的范围和状态 ( 水 平或倾斜 。

19、)， 若所述层状介质状态为倾斜， 则将炮点和检波点延拓至和地层平行的倾斜面 上。 0023 根据本发明实施例所述一种地震噪声去除方法的再一个进一步的方面， 在上述在 炮集数据中减去所述层间多次波炮集数据之前还包括， 对每一炮的炮集数据进行上述步骤 101- 步骤 106， 得到全部炮集数据的时间域的层间多次波炮集数据。 0024 本发明实施例还提供了一种地震噪声去除装置， 包括， 0025 分离单元， 用于对炮集数据进行分离， 分为层状介质数据 Slayer和下覆层数据 Syu； 0026 预测算子计算单元， 用于对炮集数据中的层状介质数据进行 -p( 倾斜叠加 ) 变 换， 利用计算每个P值。

20、对应的地震道的第n个地震数据的 预测算子， 其中， An1、 An2分别为第 n 个地震数据对应的预测算子的实部和虚部， fr、 fi 分别 为每道地震数据中第 1 到第 n 个地震数据的傅里叶变换的实部和虚部， n 取 1 至 N 的正整 数， N 为 Slayer在时间方向地震数据的个数 ； 0027 层间多次波预测算子计算单元， 利用 Anmultiple An(n)-An-1(n-1) 求取第 n 个地 震数据的多次波预测算子 Anmultiple， 对 Anmultiple做反傅里叶变换， 得到时间域多次波预测算 子 anmultiple， 令 anmultiple(n) 0, 得到。

21、 anint； 其中， anmultiple为第 n 个地震数据的多次波预测算 子 ,anint为第 n 个地震数据的层间多次波预测算子 ； 0028 层状介质层间多次波预测单元， 按照公式预测所有层状介质 相关的层间多次波 Mlayer， 其中， di为第 i 个地震数据， 为子波长度， 一般等于 25 ；* 表 示数学上的褶积运算。 0029 下覆层层间多次波预测单元， 用于计算所述层状介质数据中所有的层间多次波预 测算子的序列并求和记为 alayer， 对该炮集数据的下覆层数据 Syu进行 -p( 倾斜 叠加 ) 变换得到 STyu， 按照公式 Myu STyu*alayer得到下覆层介。

22、质在层状介质中产生的层间 多次波 Myu； 0030 层间多次波计算单元， 用于将该炮集数据的 Mlayer和 Myu进行相加， 得到 -p 域预 测的层间多次波炮集数据， 并进行 -p 反变换， 得到 X-T 域 ( 时间偏移距域 ) 层间多次波 炮集数据 ； 0031 解扰单元， 用于在炮集数据中减去所述层间多次波炮集数据。 0032 根据本发明实施例所述一种地震噪声去除装置的一个进一步的方面， 还包括偏移 单元， 与所述分离单元相连接， 对所述炮集数据进行偏移， 确定层状介质的范围和状态 ( 水 说 明 书 CN 103984023 A 6 4/8 页 7 平或倾斜 )， 若所述层状介质。

23、状态为倾斜， 则将炮点和检波点延拓至和地层平行的倾斜面 上。 0033 根据本发明实施例所述一种地震噪声去除装置的再一个进一步的方面， 还包括循 环单元， 与所述解扰单元相连接， 调用分离单元、 预测算子计算单元、 层间多次波预测算子 计算单元、 层状介质层间多次波预测单元、 下覆层层间多次波预测单元和层间多次波计算 单元计算每一炮的炮集数据中时间域的层间多次波炮集数据。 0034 通过上述实施例的方法和装置， 和 CFP 相比， 本发明方法完全由数据驱动 ； 采用 李代数积分求解预测算子， 进而得到层间多次波预测算子， 运算量由原来的 n*n 降低为 n*log(n)， 适合大规模数据处理，。

24、 可以一次预测所有上覆层层状介质相关的层间多次波。 附图说明 0035 结合以下附图阅读对实施例的详细描述， 本发明的上述特征和优点， 以及额外的 特征和优点， 将会更加清楚。 0036 图 1 所示为本发明实施例一种地震噪声去除方法的流程图 ； 0037 图 2 所示为本发明实施例一种地震噪声去除装置的结构示意图 ； 0038 图 3 所示为本发明实施例一种地震噪声去除的具体流程图 ； 0039 图 4 所示为本发明实施例层状介质部分炮集数据的层间多次波预测流程图 ； 0040 图 5 所示为本发明实施例计算下覆介质部分炮集数据在层状介质中产生的层间 多次波的流程图 ； 0041 图 6a 。

25、所示为本发明实施例炮集数据正演模型示意图 ； 0042 图 6b 所示为炮集数据示意图 ； 0043 图 6c 为本发明实施例去除层间多次波后的炮集数据示意图 ； 0044 图 6d 所示为本发明实施例 -p 域的 alayer的示意图 ； 0045 图 6e 所示为本发明实施例 -p 层间多次波预测结果示意图 ； 0046 图 6f 所示为本发明实施例 X-T 域层间多次波预测结果示意图。 具体实施方式 0047 下面的描述可以使任何本领域技术人员利用本发明。 具体实施例和应用中所提供 的描述信息仅为示例。 这里所描述的实施例的各种延伸和组合对于本领域的技术人员是显 而易见的， 在不脱离本发。

26、明的实质和范围的情况下， 本发明定义的一般原则可以应用到其 他实施例和应用中。 因此， 本发明不只限于所示的实施例， 本发明涵盖与本文所示原理和特 征相一致的最大范围。 0048 如图 1 所示为本发明实施例一种地震噪声去除方法的流程图。 0049 包括步骤 101， 对炮集数据进行分离， 分为层状介质炮集数据 Slayer和下覆层炮集 数据 Syu。其中层状介质炮集数据是指在地下呈层状排列的介质数据， 所述下覆层炮集数据 是指在地下呈不规则状的介质数据。 0050 对炮集数据中的层状介质数据 Slayer进行 -p( 倾斜叠加 ) 变换得到 STlayer， 利用 说 明 书 CN 1039。

27、84023 A 7 5/8 页 8 计算 STlayer每个 P 值对应的地震道的第 n 个地震数据的 预测算子， 其中， An1、 An2分别为第 n 个地震数据的预测算子的实部和虚部， fr、 fi 分别为 P 值对应的地震道中第 1 到第 n 个地震数据的傅里叶变换的实部和虚部， n 取 1 至 N 的正整 数， N 为 STlayer在时间方向地震数据的个数 ； 0051 步骤103， 利用AnmultipleAn(n)-An-1(n-1)求取第n个地震数据的多次波预测算 子 Anmultiple， 对 Anmultiple做反傅里叶变换， 得到时间域多次波预测算子 anmultipl。

28、e， 令 anmultiple(n) 0, 得到 anint； 其中， anmultiple为第 n 个地震数据的多次波预测算子 ,anint为第 n 个地震数 据的层间多次波预测算子 ； 0052 步骤 104， 按照公式预测所有 STlayer相关的 -p 域层间多次 波 ； 其中， di为第 i 个地震数据 , 为子波长度， 一般等于 25 ；* 表示数学上的褶积运算。 0053 步骤 105， 计算层状介质数据中所有的层间多次波预测算子序列并求和 记为 alayer， 对下覆层数据 Syu进行 -p( 倾斜叠加 ) 变换得到 STyu， 按照公式 Myu STyu*alayer得到下覆。

29、层介质在层状介质中产生的层间多次波 Myu； 0054 步骤 106， 将 Mlayer和 Myu进行相加， 得到 -p 域预测的层间多次波， 然后进行 -p 反变换， 得到 X-T 域 ( 时间偏移距域 ) 层间多次波的炮集数据 M ； 0055 步骤 107， 在炮集数据中减去计算得到的层间多次波 M， 得到不含层间多次波的炮 集数据。 0056 作为本发明的一个实施例， 在上述步骤 101 之前还包括， 对所述炮集数据进行偏 移， 确定层状介质的范围和状态， 若所述层状介质状态为倾斜， 则将炮点和检波点延拓至和 地层平行的倾斜面上。 0057 作为本发明的一个实施例， 在上述步骤 107。

30、 之前还包括， 对每一炮的炮集数据进 行上述步骤 101- 步骤 106， 得到全部炮集数据对应的层间多次波炮集数据。 0058 通过上述实施例的方法， 和 CFP 相比， 本发明方法完全由数据驱动 ； 采用李代数积 分求解预测算子， 进而得到层间多次波预测算子， 运算量由原来的 n*n 降低为 n*log(n)， 适 合大规模数据处理， 可以一次预测所有上覆层层状介质相关的层间多次波。 0059 如图 2 所示为本发明实施例一种地震噪声去除装置的结构示意图。 0060 包括分离单元 201， 用于对炮集数据进行分离， 分为层状介质数据 Slayer和下覆层 数据 Syu。其中层状介质数据是指。

31、在地下呈层状排列的介质数据， 所述下覆层数据是指在地 下呈不规则状的介质数据。 0061 预测算子计算单元， 用于对炮集数据中的层状介质数据进行 -p( 倾斜叠加 ) 变 说 明 书 CN 103984023 A 8 6/8 页 9 换， 利用计算每个P值对应的地震道的第n个地震数据的 预测算子， 其中， An1、 An2分别为第 n 个地震数据对应的预测算子的实部和虚部， fr、 fi 分别 为每 P 值地震道数据中第 1 到第 n 个地震数据的傅里叶变换的实部和虚部， n 取 1 至 N 的 正整数， N 为 Slayer在时间方向上地震数据的个数 ； 0062 层间多次波预测算子计算单元。

32、， 利用 Anmultiple An(n)-An-1(n-1) 求取第 n 个地 震数据的多次波预测算子 Anmultiple， 对 Anmultiple做反傅里叶变换， 得到时间域多次波预测算 子 anmultiple， 令 anmultiple(n) 0, 得到 anint； 其中， anmultiple为第 n 个地震数据的多次波预测算 子 ,anint为第 n 个地震数据的层间多次波预测算子 ； 0063 层状介质层间多次波预测单元， 按照公式预测所有层状介质 相关的层间多次波 Mlayer， 其中， di为第 i 个地震数据， 为子波长度， 一般等于 25 ；* 表 示数学上的褶积运。

33、算。 0064 下覆层层间多次波预测单元， 用于计算所述层状介质数据中所有的层间多次波预 测算子的序列并求和记为 alayer， 对该炮集数据的下覆层数据 Syu进行 -p( 倾斜 叠加 ) 变换得到 STyu， 按照公式 Myu STyu*alayer得到下覆层介质在层状介质中产生的层间 多次波 Myu； 0065 层间多次波计算单元， 用于将该炮集数据的 Mlayer和 Myu进行相加， 得到 -p 域预 测的层间多次波炮集数据， 并进行 -p 反变换， 得到 X-T 域 ( 时间偏移距域 ) 层间多次波 炮集数据 ； 0066 解扰单元， 用于在炮集数据中减去所述层间多次波炮集数据。 0。

34、067 作为本发明的一个实施例， 还包括偏移单元 208， 与所述分离单元 201 相连接， 对 所述炮集数据进行偏移， 确定层状介质的状态， 若所述层状介质状态为倾斜， 则将炮点和检 波点延拓至和地层平行的倾斜面上。 0068 作为本发明的一个实施例， 还包括循环单元 209， 与所述解扰单元 207 相连接， 调 用分离单元 201、 预测算子计算单元 202、 层间多次波预测算子计算单元 203、 层状介质层间 多次波预测单元204、 下覆层介质层间多次波预测单元205和层间多次波计算单元206计算 每一炮的炮集数据中时间偏移距域的层间多次波， 得到全部炮集的层间多次波的干扰。 0069。

35、 通过上述实施例的装置， 和 CFP 相比， 本发明方法完全由数据驱动 ； 采用李代数积 分求解预测算子， 进而得到层间多次波预测算子， 运算量由原来的 n*n 降低为 n*log(n)， 适 合大规模数据处理， 可以一次预测所有上覆层层状介质相关的层间多次波。 0070 如图 3 所示为本发明实施例一种地震噪声去除的具体流程图。 0071 包括步骤 301， 流程开始， 如图 6a 所示为本发明实施例炮集数据正演模型示意图， 包括 5 个水平层介质， 下覆层为起伏介质 ； 如图 6b 所示为炮集数据示意图。 0072 步骤 302， 对炮集数据进行偏移或输出偏移距剖面， 确定层状介质范围， 。

36、确定层状 说 明 书 CN 103984023 A 9 7/8 页 10 介质是否倾斜， 如果倾斜则进入步骤 303， 否则进入步骤 304。 0073 步骤 303， 将炮点和检波点延拓至和地层平行的倾斜面。 0074 步骤 304， 对叠前的炮集数据进行分离 ： 层状介质炮集数据部分 (Slayer) 和下覆层 介质炮集数据 (Syu)。 0075 步骤 305， 针对于层状介质炮集数据进行 sub1 如图 4 所示的步骤处理。 0076 步骤 306， 针对于下覆层介质炮集数据进行 sub2 如图 5 所示的步骤处理。 0077 步骤307， 将sub1与sub2的计算结果相加， 得到-。

37、p域层间多次波预测的炮集数 据(如图6e所示为本发明实施例-p层间多次波预测结果示意图)， 然后进行Tau-p反变 换， 得到该炮集数据对应的 X-T 域 ( 时间偏移距域 ) 层间多次波炮集数据 ( 如图 6f 所示为 本发明实施例时间域层间多次波预测结果示意图 )。 0078 步骤 308， 利用炮集数据减去相应的层间多次波， 得到该炮纯净的炮集数据 ( 如图 6c 为本发明实施例去除层间多次波后的炮集数据示意图 )。 0079 步骤 309， 是否处理完所有炮的炮集数据， 如果处理完毕则进入步骤 310， 否则返 回步骤 304。 0080 步骤 310， 结束。 0081 如图 4 所。

38、示为本发明实施例层状介质部分炮集数据的层间多次波预测流程图。 0082 包括步骤 401， 对层状介质部分的炮集数据进行 -p 变换。 0083 步骤 402， 对每一个 P 值所对应的地震道， 依次求取每个地震数据的预测算子， 多 次波预测算子， 然后得到每个地震数据的层间多次波预测算子。 0084 在 本 步 骤 中，对 每 一 个 P 值 对 应 的 地 震 道 数 据 按 照 公 式 求取第n(n1至N,N为Slayer在时间方向地震数据的个 数)个地震数据的预测算子， 按照公式AnmultipleAn(n)-An-1(n-1)求取第n个地震数据的 多次波预测算子 Anmultiple。

39、， 对 Anmultiple做反傅里叶变换， 得到时间域多次波预测算子 anmultiple， 令 anmultiple(n) 0, 得到 anint； 其中， anmultiple为第 n 个地震数据的多次波预测算子 ,anint为 第 n 个地震数据的层间多次波预测算子 ； 0085 步骤 403， 将每个地震数据与对应的层间多次波预测算子进行褶积， 得到该炮集数 据的层状介质部分的层间多次波 Mlayer。 0086 在本步骤中， 按照公式预测所有 -p 域层状介质相关的层间 多次波 Mlayer。其中， di为第 i 个地震数据。 * 表示数学上的褶积运算。 0087 步骤 404， 。

40、结束。 0088 如图 5 所示为本发明实施例计算下覆介质部分炮集数据录在层状介质中产生的 层间多次波的流程图。 0089 包括步骤 501， 对下覆层介质部分的炮集数据进行 -p 变换。 0090 步骤 502， 对于每一个 P 值所对应的地震数据， 计算所述层状介质数据中所有的层 说 明 书 CN 103984023 A 10 8/8 页 11 间多次波预测算子的序列并求和记为alayer。 (如图6d所示为本发明实施例-p 域的 alayer的示意图 ) 0091 步骤 503， 将进行了 -p 变换的下覆层介质部分的地震数据与上述总的层间多次 波预测算子进行褶积， 得到下覆层介质部分的。

41、炮集数据在层状介质里产生的层间 多次波 Myu。 0092 步骤 504， 结束。 0093 通过上述的方法及装置， 针对多层层状介质 ( 包括水平和倾斜的情况 ) 的特点， 本发明提出了层间多次波预测序列， 层间多次波可以简单的表示为有效波与该层和上覆层 对应的层间多次波预测算子之和的褶积， 从而实现快速的层间多次波预测。同时层间多次 波预测序列采用李代数积分表示， 使得层间多次波预测序列的计算量由原来的 n*n 降低为 n*log(n)， 基于预测算子的层间多次波预测技术可以一次预测所有层相关的层间多次波， 因此解决了 CFP 算法需要指出产生层间多次波的 层 的问题 ； 同时， 该技术预。

42、测层间多次 波的表达式仅为两重积分， 和 ISS 技术的 5 重积分相比， 运算量大幅度减少， 更适合大规模 的地震数据处理。 0094 本发明可以以任何适当的形式实现， 包括硬件、 软件、 固件或它们的任意组合。本 发明可以根据情况有选择的部分实现， 比如计算机软件执行于一个或多个数据处理器以及 数字信号处理器。 本文的每个实施例的元素和组件可以在物理上、 功能上、 逻辑上以任何适 当的方式实现。事实上， 一个功能可以在独立单元中、 在一组单元中、 或作为其他功能单元 的一部分来实现。 因此， 该系统和方法既可以在独立单元中实现， 也可以在物理上和功能上 分布于不同的单元和处理器之间。 00。

43、95 在相关领域中的技术人员将会认识到， 本发明的实施例有许多可能的修改和组 合， 虽然形式略有不同， 仍采用相同的基本机制和方法。为了解释的目的， 前述描述参考了 几个特定的实施例。 然而， 上述的说明性讨论不旨在穷举或限制本文所发明的精确形式。 前 文所示， 许多修改和变化是可能的。 所选和所描述的实施例， 用以解释本发明的原理及其实 际应用， 用以使本领域技术人员能够最好地利用本发明和各个实施例的针对特定应用的修 改、 变形。 说 明 书 CN 103984023 A 11 1/7 页 12 图 1 说 明 书 附 图 CN 103984023 A 12 2/7 页 13 图 2 说 明 书 附 图 CN 103984023 A 13 3/7 页 14 图 3 说 明 书 附 图 CN 103984023 A 14 4/7 页 15 图 4 说 明 书 附 图 CN 103984023 A 15 5/7 页 16 图 5 说 明 书 附 图 CN 103984023 A 16 6/7 页 17 图 6a 图 6b 图 6c 图 6d 说 明 书 附 图 CN 103984023 A 17 7/7 页 18 图 6e 图 6f 说 明 书 附 图 CN 103984023 A 18 。