技术摘要：
本发明提供一种多次波的压制方法、设备以及系统，涉及能源勘探技术领域。所述方法包括：获取地震数据；将所述地震数据分选到共中心点域，形成共中心点道集；根据所述共中心点道集确定共中心点道集的多次波最终模型；将所述共中心点道集与所述共中心点道集的多次波最终  全部
背景技术：
多次波是地震资料中普遍存在的一种干扰波，其类型多、波场复杂，并且与有效波 有非常相似的性质。因此，多次波的压制一直是地震资料处理中的难题。 现有技术中，去除多次波时首先需对多次波进行识别。在地震资料处理工作中，可 以从道集、速度谱、叠加剖面以及偏移剖面中进行识别。多次波识别的越清晰，其压制效果 往往也越好。对于多次波的压制方法经过近几十年的研究得到了长足的发展，基本可分为 如下三类： ①基于多次波的周期性来压制多次波，如共中心点叠加技术、拉东变换技术。 ②基于多次波与有效波的速度差异压制多次波，如τ-p域预测反褶积技术。 ③多次波模型减去法，如FK倾角滤波技术。 由于多次波的周期、速度都有时变、空变的特征，且大部分压制方法都是基于数据 或模型驱动。上述单一方法的应用无法完整的压制多次波。 因此，如何研究和开发出一种试验方案，其能够完整的压制多次波是本领域亟待 解决的技术难题。
技术实现要素：
本发明提供了一种多次波的压制方法、设备以及系统，采用地震地质一体化理念， 对多次波进行多数据域多方法的组合去除，从而克服单一方法去噪不完善的局限性，实现 保真去除多次波提高成像分辨率的目的。 本发明的目的是，提供一种多次波的压制方法，所述方法包括： 获取地震数据； 将所述地震数据分选到共中心点域，形成共中心点道集； 根据所述共中心点道集确定共中心点道集的多次波最终模型； 将所述共中心点道集与所述共中心点道集的多次波最终模型组合运算，得到去除 多次波后的共中心点道集。 在本发明的优选实施方式中，根据所述共中心点道集确定共中心点道集的多次波 最终模型包括： 根据所述共中心点道集确定共中心点道集的多次波初始模型； 根据所述共中心点道集的多次波初始模型确定多次波速度； 利用所述多次波速度对所述共中心点道集进行多次波动校正； 利用倾角滤波在所述共中心点道集的多次波初始模型上对所述多次波进行预测， 得到共中心点道集的多次波最终模型。 4 CN 111596358 A 说　明　书 2/9 页 在本发明的优选实施方式中，根据所述共中心点道集确定共中心点道集的多次波 初始模型包括： 利用一次波速度对所述共中心点道集进行动校正； 将动校正后的共中心点道集进行分选，形成共偏移距道集； 在所述共偏移距道集中，分析产生多次波的目标层； 根据所述共偏移距道集结合所述目标层确定目标区域； 利用倾角滤波对所述目标区域里的多次波进行预测，得到所述共偏移距道集的多 次波初始模型； 将所述共偏移距道集的多次波初始模型分选，得到共中心点道集的多次波初始模 型。 在本发明的优选实施方式中，根据所述共中心点道集确定共中心点道集的多次波 初始模型还包括当所述目标层存在地层倾角时，对所述目标层进行层拉平处理。 在本发明的优选实施方式中，根据所述共中心点道集的多次波初始模型确定多次 波速度包括： 利用一次波速度对所述共中心点道集的多次波初始模型进行反动校正； 根据反动校正后的共中心点道集的多次波初始模型制作速度谱； 根据所述速度谱进行多次波速度解释，得到多次波速度。 在本发明的优选实施方式中，根据所述共中心点道集确定共中心点道集的多次波 最终模型还包括：对所述共中心点道集的多次波最终模型进行多次波反动校正。 本发明的目的是，提供一种多次波的压制系统，所述系统包括： 地震数据获取模块，用于获取地震数据； 地震数据分选模块，用于将所述地震数据分选到共中心点域，形成共中心点道集； 最终模型确定模块，用于根据所述共中心点道集确定共中心点道集的多次波最终 模型； 多次波压制模块，用于将所述共中心点道集与所述共中心点道集的多次波最终模 型组合运算，得到去除多次波后的共中心点道集。 在本发明的优选实施方式中，所述最终模型确定模块包括： 初始模型确定模块，用于根据所述共中心点道集确定共中心点道集的多次波初始 模型； 多次波速度确定模块，用于根据所述共中心点道集的多次波初始模型确定多次波 速度； 多次波动校正模块，用于利用所述多次波速度对所述共中心点道集进行多次波动 校正； 多次波预测模块，用于利用倾角滤波在所述共中心点道集的多次波初始模型上对 所述多次波进行预测，得到共中心点道集的多次波最终模型。 在本发明的优选实施方式中，所述初始模型确定模块包括： 动校正模块，用于利用一次波速度对所述共中心点道集进行动校正； 道集分选模块，用于将动校正后的共中心点道集进行分选，形成共偏移距道集； 目标层分析模块，用于在所述共偏移距道集中，分析产生多次波的目标层； 5 CN 111596358 A 说　明　书 3/9 页 目标区域确定模块，用于根据所述共偏移距道集结合所述目标层确定目标区域； 模型确定模块，用于利用倾角滤波对所述目标区域里的多次波进行预测，得到所 述共偏移距道集的多次波初始模型； 初始模型分选模块，用于将所述共偏移距道集的多次波初始模型分选，得到共中 心点道集的多次波初始模型。 在本发明的优选实施方式中，所述初始模型确定模块还包括： 拉平处理模块，用于当所述目标层存在地层倾角时，对所述目标层进行层拉平处 理。 在本发明的优选实施方式中，所述多次波速度确定模块包括： 初始模型反动校正模块，用于利用一次波速度对所述共中心点道集的多次波初始 模型进行反动校正； 速度谱制作模块，用于根据反动校正后的共中心点道集的多次波初始模型制作速 度谱； 速度解释模块，用于根据所述速度谱进行多次波速度解释，得到多次波速度。 在本发明的优选实施方式中，所述最终模型确定模块还包括： 最终模型反动校正模块，用于对所述共中心点道集的多次波最终模型进行多次波 反动校正。 本发明的目的之一是，提供一种计算机设备，包括：适于实现各指令的处理器以及 存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行一种多次波 的压制方法。 本发明的目的之一是，提供一种计算机可读存储介质，存储有计算机程序，所述计 算机程序用于执行一种多次波的压制方法。 本发明的有益效果在于，提供了一种多次波的压制方法、设备以及系统，采用地震 地质一体化理念，对多次波进行多数据域多方法的组合去除，从而克服单一方法去噪不完 善的局限性，实现保真去除多次波提高成像分辨率的目的，具有保真压制多次波的目的能 力，在有效去除多次波干扰的前提下不损伤有效信号，大幅通过成像分辨率，对地质构造的 刻画更加准确。 为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例， 并配合所附图式，作详细说明如下。 附图说明 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。 图1为本发明实施例提供的一种多次波的压制方法的流程图； 图2为图1中的步骤S103的具体流程图； 图3为图2中的步骤S201的具体流程图； 图4为图2中的步骤S202的具体流程图； 6 CN 111596358 A 说　明　书 4/9 页 图5为本发明实施例提供的一种多次波的压制系统的结构框图； 图6为本发明实施例提供的一种多次波的压制系统中的最终模型确定模块的结构 框图； 图7为本发明实施例提供的一种多次波的压制系统中的初始模型确定模块的结构 框图； 图8为本发明实施例提供的一种多次波的压制系统中的多次波速度确定模块的结 构框图； 图9为本发明提供的具体实施例中地震数据分选形成CMP道集的示意图； 图10为本发明提供的具体实施例中一次波动校正后的CMP道集的示意图； 图11为本发明提供的具体实施例中地震数据分选形成共偏移距道集的示意图； 图12、图13为本发明提供的具体实施例中共偏移距道集多次波预测示意图； 图14为本发明提供的具体实施例中多次波速度解释的示意图； 图15为本发明提供的具体实施例中CMP域最终多次波模型的示意图； 图16为本发明提供的具体实施例中多次波压制前的CMP道集示意图； 图17为本发明提供的具体实施例中多次波压制后的CMP道集示意图； 图18为本发明提供的具体实施例中多次波压制前偏移剖面的示意图； 图19为本发明提供的具体实施例中多次波压制后偏移剖面的示意图。