技术摘要：
本发明提供了一种地震零值切片自动拾取方法及系统，包括：通过最大期望算法求解多个地震切片中每个地震切片具有不同聚类数目的混合高斯模型；根据贝叶斯信息准则选取每个地震切片的最优聚类数目；根据所有地震切片的最优聚类数目的统计分布确定零值切片的估计值和估计  全部
背景技术：
地震数据“高信噪比、高分辨率、高保真度”三高要求中，高分辨率是资料处理的基 本要求，地震数据的分辨率对后续处理以及综合解释有着重要意义。随着油气勘探不断向 复杂地区推进，分辨率的高低极大限制勘探工作的进程，尤其普通的处理技术不容易识别 薄层及复杂构造，提高分辨率的相关工作势在必行。 对于复杂构造而言，在大部分反射时间段，多个地质体的反射是重叠的，在地震切 片上会同时出现多个地质体的反射，但存在这样的时刻，在该时刻某一个地质体地震反射 的样点值为零，对应的地震切片上只出现其余几个地质体的反射，该时刻称为地质体的零 值时刻，相应切片称为零值切片。零值时刻的提取有助于识别地质体，如薄层砂体、河道，有 着提高地震勘探分辨率的潜能。但目前基本通过人工拾取零值切片，拾取过程耗时耗力，当 地震切片数量较少时，人为拾取勉强可以实现，但当地震切片数量为成千上万时，人工拾取 需要耗费大量的人力物力，时间长成本高。并且，对于传统的人工拾取方法，大多经过对原 始数据处理、加工和计算，但是并没有对原始数据本身的特征进行提炼分析。同时，当需要 计算的数据不满足先验假设及要求时，传统的人工拾取方法也会受到干扰，导致拾取结果 不准确。
技术实现要素：
本发明的一个目的在于提供一种地震零值切片自动拾取方法，快速、自动拾取零 值切片，在效率和精度上都高于人工，节约人力物力。本发明的另一个目的在于提供一种地 震零值切片自动拾取系统。本发明的再一个目的在于提供一种计算机设备。本发明的还一 个目的在于提供一种可读介质。 为了达到以上目的，本发明一方面公开了一种地震零值切片自动拾取方法，包括： 通过最大期望算法求解多个地震切片中每个地震切片具有不同聚类数目的混合 高斯模型； 根据贝叶斯信息准则选取每个地震切片的最优聚类数目； 根据所有地震切片的最优聚类数目的统计分布确定零值切片的估计值和/或估计 区间。 优选地，所述根据贝叶斯信息准则选取每个地震切片的最优聚类数目具体包括： 得到不同聚类数目的混合高斯模型的贝叶斯信息量； 确定所述贝叶斯信息量最小的聚类数目为每个地震切片的最优聚类数目。 优选地，所述根据所有地震切片的最优聚类数目的统计分布确定零值切片的估计 值和/或估计区间具体包括： 4 CN 111580156 A 说　明　书 2/11 页 根据最优聚类数目确定所述多个地震切片中的地质体数量； 根据所述地质体数量确定零值切片对应的聚类数目阈值； 最优聚类数目低于所述聚类数目阈值的连续地震切片为目标零值切片以确定零 值切片的估计值和/或估计区间。 优选地，所述方法进一步包括： 选取最优聚类数目低于阈值的目标地震切片的中间时刻作为零值切片的估计值； 选取最优聚类数目低于阈值的目标地震切片对应的时间范围作为零值切片的估 计区间。 本发明另一方面公开了一种地震零值切片自动拾取系统，包括： 模型建立模块，用于通过最大期望算法求解多个地震切片中每个地震切片具有不 同聚类数目的混合高斯模型； 切片聚类模块，用于根据贝叶斯信息准则选取每个地震切片的最优聚类数目； 零值时刻拾取模块，用于根据所有地震切片的最优聚类数目的统计分布确定零值 切片的估计值和/或估计区间。 优选地，所述切片聚类模块进一步用于得到不同聚类数目的混合高斯模型的贝叶 斯信息量，确定所述贝叶斯信息量最小的聚类数目为每个地震切片的最优聚类数目。 优选地，所述零值时刻拾取模块进一步用于根据最优聚类数目确定所述多个地震 切片中的地质体数量，根据所述地质体数量确定零值切片对应的聚类数目阈值，最优聚类 数目低于所述聚类数目阈值的连续地震切片为目标零值切片以确定零值切片的估计值和/ 或估计区间。 优选地，所述系统进一步包括估计模块； 所述估计模块用于选取最优聚类数目低于阈值的目标地震切片的中间时刻作为 零值切片的估计值，选取最优聚类数目低于阈值的目标地震切片对应的时间范围作为零值 切片的估计区间。 零值时刻拾取模块零值时刻拾取模块本发明还公开了一种计算机设备，包括存储 器、处理器以及存储在存储器上并可在处理器上运行的计算机程序， 所述处理器执行所述程序时实现如上所述方法。 本发明还公开了一种计算机可读介质，其上存储有计算机程序， 该程序被处理器执行时实现如上所述方法。 本发明通过采用不同聚类数目的混合高斯模型对地震切片的原始数据进行属性 聚类分析，得到每个地震切片的最优聚类数目，并通过多个地震切片的最优聚类数目可确 定多个地震切片中的零值切片。本发明通过采用混合高斯模型对原始数据进行属性聚类分 析，极大地挖掘原始数据本身的属性特征，对数据没有严格的要求与假设。实现零值切片的 自动拾取，可以有效取代人工，提高拾取效率和精度，进而提高薄层砂体、河道等地质体识 别的精度。 附图说明 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 5 CN 111580156 A 说　明　书 3/11 页 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。 图1示出本发明一种地震零值切片自动拾取方法具体实施例的流程图之一； 图2示出本发明一种地震零值切片自动拾取方法具体实施例的流程图之二； 图3示出本发明一种地震零值切片自动拾取方法具体实施例的流程图之三； 图4示出本发明一种地震零值切片自动拾取方法具体实施例的流程图之四； 图5示出本发明一种地震零值切片自动拾取方法具体实施例的流程图之五； 图6示出本发明一种地震零值切片自动拾取方法具体实施例的流程图之六； 图7示出本发明一种地震零值切片自动拾取方法具体实施例中一种零值切片的样 例图； 图8示出本发明一种地震零值切片自动拾取方法具体实施例中另一种零值切片的 样例图； 图9示出本发明一种地震零值切片自动拾取方法具体实施例中另一种非零值切片 的样例图； 图10示出本发明一种地震零值切片自动拾取方法具体实施例样本数据加入10％ 随机噪声示意图； 图11示出本发明一种地震零值切片自动拾取方法具体实施例样本数据加入10％ 随机噪声的最优聚类数目统计图； 图12示出本发明一种地震零值切片自动拾取方法具体实施例样本数据加入10％ 随机噪声的最优聚类数目随时间变化的统计图； 图13示出本发明一种地震零值切片自动拾取方法具体实施例实际三维数据体中 的一个切片示意图； 图14示出本发明一种地震零值切片自动拾取方法具体实施例实际三维数据体中 最优聚类数目随时间变化的统计图； 图15示出本发明一种地震零值切片自动拾取方法具体实施例实际三维数据体零 值时刻的估计值和估计区间与人工精细化处理结果的的对比图； 图16示出本发明一种地震零值切片自动拾取系统具体实施例的结构图之一； 图17示出本发明一种地震零值切片自动拾取系统具体实施例的结构图之二； 图18示出适于用来实现本发明实施例的计算机设备的结构示意图。