技术摘要：
本发明提供一种识别多次波的方法和装置，该方法包括：获取地震工区的时间层位数据和测井数据；根据时间层位数据和测井数据，确定地震工区的波阻抗体；将波阻抗体转换为一次波反射率；根据一次波反射率和测井数据，正演合成一次波地震记录和全波地震记录；根据一次波地  全部
背景技术：
在石油勘探中，地震波通过地下岩层传播，当遇到强反射界面时，则产生强能量的 反射波。反射波返回向上传播遇到良好的强反射界面时，再次反射并向下传播，如此往返就 形成了多次波。因此通常情况下，我国各个探区都不同程度地存在多次波，特别是在诸如四 川盆地、塔里木盆地和鄂尔多斯盆地，构造平缓、强反射界面较多，多次波干扰更为严重。这 类多次波与深层一次波振幅相当，动校正时差较小，与深层一次波相干涉，难以识别，影响 地震数据偏移质量、地震属性提取以及反演等地震解释研究的可靠性。 由于多次波成因复杂，在地震资料处理过程中识别和压制难度大，会造成有效波 反射能量畸变、使储层预测结果错误，给油气勘探和开发带来很大的风险，因此，有必要对 其分布进行识别。 目前地震处理解释中，通常基于多次波的速度和周期性特征在速度谱、叠加剖面、 地震叠前道集上识别多次波：在叠前阶段，速度谱存在低速能量团、CMP道集上出现动校不 足，随偏移距增大具有下拉现象、深层和浅层具有相似旅行周期；在叠后阶段，通过叠加地 震剖面上强反射界面和强反射界面之下地震反射的形态和时间识别多次波；当地下反射界 面倾角较小时，如果在强反射界面下方存在另一组与其形态一致，且出现时间小于其2倍的 反射同相轴时，就可以判断该组同相轴为多次反射波；对于倾斜地层，叠加剖面上二次反射 波的倾角约为其一次反射波的2倍。 但在现有技术中，都只能实现单个观测点(一个叠前地震道集或一个地震速度谱) 和一条观测线(叠加地震剖面)上的多次波识别，无法得到多次波在地震工区平面上的整体 发育程度和分布情况。由于地下构造和地层厚度变化很快，导致多次波在平面上也具有强 烈的非均质特点，单个点和单条线的结果显然无法代表整个地震工区的情况，满足不了多 次波识别和压制效果评价的需要。
技术实现要素：
本发明实施例提供一种识别多次波的方法，利用反射率正演模拟方法得到多次波 在某个地质层在平面上的能量分布，其结果可用于指导处理人员合理压制地震数据中存在 的多次波，也可以用于解释人员评价地震资料品质，给储层预测结果进行可靠性评价，降低 开发勘探中井位部署的风险，该方法包括： 获取地震工区的时间层位数据和测井数据； 根据时间层位数据和测井数据，确定地震工区的波阻抗体； 将波阻抗体转换为一次波反射率； 根据一次波反射率和测井数据，正演合成一次波地震记录和全波地震记录； 根据一次波地震记录和全波地震记录，从地震工区中确定不同干扰等级的多次波 5 CN 111610563 A 说　明　书 2/9 页 干扰区域。 本发明实施例还提供一种识别多次波的装置，包括： 数据获取模块，用于获取地震工区的时间层位数据和测井数据； 波阻抗体确定模块，用于根据时间层位数据和测井数据，确定地震工区的波阻抗 体； 一次波反射率转换模块，用于将波阻抗体转换为一次波反射率； 正演合成模块，用于根据一次波反射率和测井数据，正演合成一次波地震记录和 全波地震记录； 多次波干扰区域确定模块，用于根据一次波地震记录和全波地震记录，从地震工 区中确定不同干扰等级的多次波干扰区域。 本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并 可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述识别多次波 的方法。 本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有 执行实现上述识别多次波的方法的计算机程序。 本发明实施例提供的一种识别多次波的方法和装置，通过建立波阻抗体，利用反 射率正演模拟的方法正演合成一次波地震记录和全波地震记录；通过对全波地震记录和一 次波地震记录分析，得到多次波在某个地质层在平面上的能量分布，相对现有技术中多次 波单点、单线的识别方法，本发明实施例能够对整个地震工区多次波的整体发育程度和分 布情况进行预测，用于分析潜在多次波特征，应用在多次波压制处理中，从而实现多次波压 制处理参数优选，一定程度上减少了压制多次波的不确定性；本发明实施例还能够根据划 分的不同干扰等级多次波干扰区域，评价地震资料品质，划分出地震资料不可靠的区域，从 而减小储层预测的多解性，降低井位部署的风险。 附图说明 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。在附图中： 图1为本发明实施例一种识别多次波的方法的示意图。 图2为本发明实施例一种识别多次波的方法的应用实例建立的波阻抗体在某条地 震测线上显示的结果示意图。 图3为本发明实施例一种识别多次波的方法的应用实例的一次波地震记录地震剖 面图。 图4为本发明实施例一种识别多次波的方法的确定全波反射率的流程图。 图5为本发明实施例一种识别多次波的方法的应用实例的全波地震记录地震剖面 图。 图6为本发明实施例一种识别多次波的方法的应用实例中相同位置的多次波地震 记录图。 6 CN 111610563 A 说　明　书 3/9 页 图7为本发明实施例一种识别多次波的方法的应用实例获取的信噪比属性切片 图。 图8为本发明实施例一种识别多次波的装置示意图。