[发明专利]一种地震零值切片自动拾取方法及系统

说明书

本发明提供了一种地震零值切片自动拾取方法及系统，包括：通过最大期望算法求解多个地震切片中每个地震切片具有不同聚类数目的混合高斯模型；根据贝叶斯信息准则选取每个地震切片的最优聚类数目；根据所有地震切片的最优聚类数目的统计分布确定零值切片的估计值和估计区间，本发明可快速、自动拾取零值切片，在效率和精度上都高于人工，节约人力物力。

技术领域

本发明涉及勘探地球物理技术领域，尤其涉及一种地震零值切片自动拾取方法及系统。

背景技术

地震数据“高信噪比、高分辨率、高保真度”三高要求中，高分辨率是资料处理的基本要求，地震数据的分辨率对后续处理以及综合解释有着重要意义。随着油气勘探不断向复杂地区推进，分辨率的高低极大限制勘探工作的进程，尤其普通的处理技术不容易识别薄层及复杂构造，提高分辨率的相关工作势在必行。

对于复杂构造而言，在大部分反射时间段，多个地质体的反射是重叠的，在地震切片上会同时出现多个地质体的反射，但存在这样的时刻，在该时刻某一个地质体地震反射的样点值为零，对应的地震切片上只出现其余几个地质体的反射，该时刻称为地质体的零值时刻，相应切片称为零值切片。零值时刻的提取有助于识别地质体，如薄层砂体、河道，有着提高地震勘探分辨率的潜能。但目前基本通过人工拾取零值切片，拾取过程耗时耗力，当地震切片数量较少时，人为拾取勉强可以实现，但当地震切片数量为成千上万时，人工拾取需要耗费大量的人力物力，时间长成本高。并且，对于传统的人工拾取方法，大多经过对原始数据处理、加工和计算，但是并没有对原始数据本身的特征进行提炼分析。同时，当需要计算的数据不满足先验假设及要求时，传统的人工拾取方法也会受到干扰，导致拾取结果不准确。

发明内容

本发明的一个目的在于提供一种地震零值切片自动拾取方法，快速、自动拾取零值切片，在效率和精度上都高于人工，节约人力物力。本发明的另一个目的在于提供一种地震零值切片自动拾取系统。本发明的再一个目的在于提供一种计算机设备。本发明的还一个目的在于提供一种可读介质。

为了达到以上目的，本发明一方面公开了一种地震零值切片自动拾取方法，包括：

通过最大期望算法求解多个地震切片中每个地震切片具有不同聚类数目的混合高斯模型；

根据贝叶斯信息准则选取每个地震切片的最优聚类数目；

根据所有地震切片的最优聚类数目的统计分布确定零值切片的估计值和/或估计区间。

优选地，所述根据贝叶斯信息准则选取每个地震切片的最优聚类数目具体包括：

得到不同聚类数目的混合高斯模型的贝叶斯信息量；

确定所述贝叶斯信息量最小的聚类数目为每个地震切片的最优聚类数目。

优选地，所述根据所有地震切片的最优聚类数目的统计分布确定零值切片的估计值和/或估计区间具体包括：

根据最优聚类数目确定所述多个地震切片中的地质体数量；

根据所述地质体数量确定零值切片对应的聚类数目阈值；

最优聚类数目低于所述聚类数目阈值的连续地震切片为目标零值切片以确定零值切片的估计值和/或估计区间。

优选地，所述方法进一步包括：

选取最优聚类数目低于阈值的目标地震切片的中间时刻作为零值切片的估计值；

选取最优聚类数目低于阈值的目标地震切片对应的时间范围作为零值切片的估计区间。

本发明另一方面公开了一种地震零值切片自动拾取系统，包括：

模型建立模块，用于通过最大期望算法求解多个地震切片中每个地震切片具有不同聚类数目的混合高斯模型；