CN114626225A - 一种基于设置拟合优度阈值的热史模拟方法、系统及设备 - Google Patents

Abstract

本发明涉及一种基于设置拟合优度阈值的热史模拟方法、系统及设备，该方法通过待模拟地质在垂直剖面上M组第一热史曲线集合，对M组第一热史曲线集合进行分析、处理，得到符合需求待模拟地质在垂直剖面上的第一热史模拟结果，通过设置不同阈值得到多个第一热史模拟结果，并从多个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果，该热史模拟结果准确率高，更符合实际，解决了现有采用低温年代学数据实现垂直剖面热史模拟方法存在模拟数据不准确的技术问题。

Description

一种基于设置拟合优度阈值的热史模拟方法、系统及设备

技术领域

本发明涉及地质技术领域，尤其涉及一种基于设置拟合优度阈值的热史模拟方法、系统及设备。

背景技术

近年来，基于低温年代学数据的热史模拟已成为近年来地质勘探领域的一项重要工具，对恢复岩石垂向移动的动态过程起到了关键作用。但是在热史模拟过程中，也存在一些问题，比如热史模拟的精度难以控制，尤其对于垂直剖面热史模拟就是一个亟待解决的问题。

发明内容

本发明实施例提供了一种基于设置拟合优度阈值的热史模拟方法、系统及设备，用于解决现有采用低温年代学数据实现垂直剖面热史模拟方法存在模拟数据不准确的技术问题。

为了实现上述目的，本发明实施例提供如下技术方案：

一种基于设置拟合优度阈值的热史模拟方法，包括以下步骤：

S1.获取待模拟地质在垂直剖面上M组第一热史曲线集合；

S2.采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个所述第一拟合优度；

S3.将M个所述第一拟合优度与优度阈值对比，筛选得到K个所述第一拟合优度；并对K个所述第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果；

S4.设置不同的优度阈值再按照步骤S3得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果；

其中，K小于M，K、M、A均为大于5的自然数，所述拟合优度公式为：

式中，GOF为第一拟合优度，age为热史曲线的模拟年龄，

为垂直剖面上所有样品实验测试年龄的平均误差，

为垂直剖面上样品模拟年龄与实验年龄之间差值的绝对值的平均。

优选地，获取待模拟地质在垂直剖面上M组第一热史曲线集合的步骤包括：

获取模拟地质在垂直剖面上顶端的样品作为首样品，所述首样品包括M条热史曲线；

以所述首样品的高度为基准，通过古地温梯度函数从模拟地质的垂直剖面上获取N-1个样品以及与每个所述样品对应的温差；

将所述首样品的M条热史曲线通过所述温差平移，得到对应所述样品的M条热史曲线；

从N个所述样品提取同一纵列的热史曲线组成一组第一热史曲线集合，每个所述样品包括M条热史曲线，则组成M组第一热史曲线集合。

优选地，得到每组第一热史曲线集合的第一拟合优度包括：对每组第一热史曲线集合通过U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法进行计算，得到对应的每组第一热史曲线集合的拟合优度，并从中选取拟合优度数值最小的作为该组第一热史曲线集合的第一拟合优度。

优选地，该基于设置拟合优度阈值的热史模拟方法包括：

获取与所述第二热史模拟结果对应的优度阈值作为最终优度阈值；

通过所述最终优度阈值对M组第一热史曲线集合的热史曲线进行筛选，得到M组第二热史曲线集合；

从M组所述第二热史曲线集合中提取每个样品的数条热史曲线；

在同一时间点上，对提取后每个样品中所有热史曲线上的每个数据求取平均，得到每个数据的平均值；

按时间顺序将所有所述数据的平均值连接成曲线，该曲线作为该样品的最终热史模拟结果。

优选地，该基于设置拟合优度阈值的热史模拟方法包括：将N个所述样品的最终热史模拟结果作为该待模拟地质的热史模拟结果。

优选地，该基于设置拟合优度阈值的热史模拟方法包括：获取待模拟地质在垂直剖面上至少10000组第一热史曲线集合。

本发明还提供一种基于设置拟合优度阈值的热史模拟系统，包括数据获取模块、拟合优度计算模块、筛选计算模块和第一输出结果模块；

所述数据获取模块，用于获取待模拟地质在垂直剖面上M组第一热史曲线集合；

所述拟合优度计算模块，用于采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个所述第一拟合优度；

所述筛选计算模块，用于将M个所述第一拟合优度与优度阈值对比，筛选得到K个所述第一拟合优度；并对K个所述第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果；

所述第一输出结果模块，用于设置不同的优度阈值再通过所述筛选计算模块得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果；

其中，K小于M，K、M、A均为大于5的自然数，所述拟合优度公式为：

式中，GOF为第一拟合优度，age为热史曲线的模拟年龄，

为垂直剖面上所有样品实验测试年龄的平均误差，

为垂直剖面上样品模拟年龄与实验年龄之间差值的绝对值的平均。

优选地，所述数据获取模块包括第一样品获取子模块、第二样品获取子模块、第一数据获取子模块和第二数据获取子模块；

所述第一样品获取子模块，用于获取模拟地质在垂直剖面上顶端的样品作为首样品，所述首样品包括M条热史曲线；

所述第二样品获取子模块，用于以所述首样品的高度为基准，通过古地温梯度函数从模拟地质的垂直剖面上获取N-1个样品以及与每个所述样品对应的温差；

所述第一数据获取子模块，用于将所述首样品的M条热史曲线通过所述温差平移，得到对应所述样品的M条热史曲线；

所述第二数据获取子模块，用于从N个所述样品提取同一纵列的热史曲线组成一组第一热史曲线集合，每个所述样品包括M条热史曲线，则组成M组第一热史曲线集合。

优选地，该基于设置拟合优度阈值的热史模拟系统包括第二输出结果模块，所述第二输出结果模块用于获取与所述第二热史模拟结果对应的优度阈值作为最终优度阈值，通过所述最终优度阈值对M组第一热史曲线集合的热史曲线进行筛选，得到M组第二热史曲线集合，从M组所述第二热史曲线集合中提取每个样品的数条热史曲线，在同一时间点上，对提取后每个样品中所有热史曲线上的每个数据求取平均，得到每个数据的平均值，按时间顺序将所有所述数据的平均值连接成曲线，该曲线作为该样品的最终热史模拟结果。

本发明还提供一种基于设置拟合优度阈值的热史模拟设备，包括处理器以及存储器；

所述存储器，用于存储程序代码，并将所述程序代码传输给所述处理器；

所述处理器，用于根据所述程序代码中的指令执行上述所述的基于设置拟合优度阈值的热史模拟方法。

从以上技术方案可以看出，本发明实施例具有以下优点：该基于设置拟合优度阈值的热史模拟方法、系统及设备，该方法包括S1.获取待模拟地质在垂直剖面上M组第一热史曲线集合；S2.采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个第一拟合优度；S3.将M个第一拟合优度与优度阈值对比，筛选得到K个第一拟合优度；并对K个第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果；S4.设置不同的优度阈值再按照步骤S3得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果。该基于设置拟合优度阈值的热史模拟方法通过待模拟地质在垂直剖面上M组第一热史曲线集合，对M组第一热史曲线集合进行分析、处理，得到符合需求待模拟地质在垂直剖面上的第一热史模拟结果，通过设置不同阈值得到多个第一热史模拟结果，并从多个第一热史模拟结果筛选出符合需求待模拟地质在垂直剖面上的第二热史模拟结果，该热史模拟结果准确率高，解决了现有采用低温年代学数据实现垂直剖面热史模拟方法存在模拟数据不准确的技术问题。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。

图1为本发明实施例所述的基于设置拟合优度阈值的热史模拟方法的步骤流程图；

图2为本发明实施例所述的基于设置拟合优度阈值的热史模拟方法的步骤流程图；

图3为本发明实施例的基于设置拟合优度阈值的热史模拟方法一组热史曲线集合的示意图；

图4为本发明实施例的基于设置拟合优度阈值的热史模拟方法拟合优度的分布直方图；

图5为本发明实施例的基于设置拟合优度阈值的热史模拟方法拟合优度的散点图；

图6为本发明实施例所述的基于设置拟合优度阈值的热史模拟系统的框架图。

具体实施方式

为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

本申请实施例提供了一种基于设置拟合优度阈值的热史模拟方法、系统及设备，用于解决了现有采用低温年代学数据实现垂直剖面热史模拟方法存在模拟数据不准确的技术问题。

实施例一：

图1为本发明实施例所述的基于设置拟合优度阈值的热史模拟方法的步骤流程图。

如图1所示，本发明实施例提供了一种基于设置拟合优度阈值的热史模拟方法，包括以下步骤：

S1.获取待模拟地质在垂直剖面上M组第一热史曲线集合。

需要说明的是，主要是获取待模拟地质在垂直剖面上M组第一热史曲线集合，为后续步骤S2至步骤S4提供数据基础。

S2.采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个第一拟合优度。

在本发明实施例中，拟合优度公式为：

式中，GOF为第一拟合优度，age为热史曲线的模拟年龄，

为垂直剖面上所有样品实验测试年龄的平均误差，

为垂直剖面上样品模拟年龄与实验年龄之间差值的绝对值的平均。

需要说明的是，主要是通过拟合优度公式计算每组第一热史曲线集合的拟合优度，得到M个第一拟合优度，为步骤S3提供数据。

S3.将M个第一拟合优度与优度阈值对比，筛选得到K个第一拟合优度；并对K个第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果。

需要说明的是，主要是将步骤S2中得到的M个第一拟合优度与其中一个优度阈值比较，将大于优度阈值的第一拟合优度筛选出，得到K个大于优度阈值的第一拟合优度，之后对筛选出的K个第一拟合优度求取其平均值，该平均值作为待模拟地质在垂直剖面上的第一热史模拟结果。其中，K小于M，K、M、A均为大于5的自然数。优度阈值可以根据需求设置，此处不做限定。

S4.设置不同的优度阈值再按照步骤S3得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果。其中，A为大于1的自然数。

需要说明的是，主要是设置不同的优度阈值再通过步骤S3得到A个第一热史模拟结果，从A个第一热史模拟结果选取其拟合优度最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果。在本发明实施例中，不断重复设置不同的优度阀值，可以获取不同优度阀值对应的第一热史模拟结果，然后将得到A个第一热史模拟结果进行对比，从A个第一热史模拟结果选取最为合适的拟合优度阀值，确保拟合优度的精度。由于拟合优度越高代表热史模拟精度越高，拟合优度数值越小而说明平均热史越平滑。在保证获取足够数量的热史曲线情况下(一般至少获得大约10组热史曲线)，取最大的可接受的优度阀值作为最为合适的拟合优度阀值，可根据实际需求选取最大拟合优度数值的第一热史模拟结果作为最为适合的拟合优度阀值，此处不做具体限定。

在本发明实施例中，该基于设置拟合优度阈值的热史模拟方法包括：将得到的M个第一拟合优度绘制频数分布直方图和散点图，便于查看。

需要说明的是，可以将得到的M个第一拟合优度在Excel中绘制频数分布直方图和散点图，通过分布直方图和散点图，可以更为直观地筛选出不同优度阀值所能筛选得到的热史曲线。

本发明提供的一种基于设置拟合优度阈值的热史模拟方法，包括S1.获取待模拟地质在垂直剖面上M组第一热史曲线集合；S2.采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个第一拟合优度；S3.将M个第一拟合优度与优度阈值对比，筛选得到K个第一拟合优度；并对K个第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果；S4.设置不同的优度阈值再按照步骤S3得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果。该基于设置拟合优度阈值的热史模拟方法通过待模拟地质在垂直剖面上M组第一热史曲线集合，对M组第一热史曲线集合进行分析、处理，得到符合需求待模拟地质在垂直剖面上的第一热史模拟结果，通过设置不同阈值得到多个第一热史模拟结果，并从多个第一热史模拟结果筛选出符合需求待模拟地质在垂直剖面上的第二热史模拟结果，该热史模拟结果准确率高，解决了现有采用低温年代学数据实现垂直剖面热史模拟方法存在模拟数据不准确的技术问题。

在本发明的一个实施例中，在步骤S1中，获取待模拟地质在垂直剖面上M组第一热史曲线集合的步骤包括：

获取模拟地质在垂直剖面上顶端的样品作为首样品，首样品包括M条热史曲线；

以首样品的高度为基准，通过古地温梯度函数从模拟地质的垂直剖面上获取N-1个样品以及与每个样品对应的温差；

将首样品的M条热史曲线通过温差平移，得到对应样品的M条热史曲线；

从N个样品提取同一纵列的热史曲线组成一组第一热史曲线集合，每个样品包括M条热史曲线，则组成M组第一热史曲线集合。其中，N为大于1的自然数。

需要说明的是，主要是对获取待模拟地质在垂直剖面上M组第一热史曲线集合方法进行描述，在本实施例中，以岩石作为待模拟地质进行举例说明，选择垂直剖面上的顶部样品作为首样品，基于时间与温度接受范围内选取其中首样品的至少10000条热史曲线，并采用蒙特卡洛方法从大于10000条热史曲线中搜索得到M条热史曲线，其它垂直剖面上的样品通过与首样品的高差乘以古地温梯度函数，获得该样品与首样品的温差，基于温差在M条热史曲线的基础上上移、下移或变形得到与其它N-1个样品对应的M条热史曲线。其中，温差＝古地温梯度*高度，现今地温梯度与古地温梯度相等。

在本发明的一个实施例中，该基于设置拟合优度阈值的热史模拟方法得到每组第一热史曲线集合的第一拟合优度除了通过拟合优度公式计算得到，还可以包括：对每组第一热史曲线集合通过U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法进行计算，得到对应的每组第一热史曲线集合的拟合优度，并从中选取拟合优度数值最小的作为该组第一热史曲线集合的第一拟合优度。

需要说明的是，该基于设置拟合优度阈值的热史模拟方法通过U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法进行计算，并从中选取上述U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法等几种方式所得拟合优度中数值最小的作为该组第一热史曲线集合的第一拟合优度，例如通过U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法这4种方法得来的拟合优度分别为GOF1～GOF4，那么以它们的最小作为这一组曲线的第一拟合优度，即是GOF＝Min(GOF1、…、GOF4)。通过U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法进行计算得到对应的每组第一热史曲线集合的拟合优度，是本领域比较成熟的技术，此实施例中不再对U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法如何计算拟合优度的内容进行详细描述。

图2为本发明实施例所述的基于设置拟合优度阈值的热史模拟方法的步骤流程图。

在本发明的一个实施例中，该基于设置拟合优度阈值的热史模拟方法包括：步骤S5，步骤S5的内容为：

获取与第二热史模拟结果对应的优度阈值作为最终优度阈值；

通过最终优度阈值对M组第一热史曲线集合的热史曲线进行筛选，得到M组第二热史曲线集合；

从M组第二热史曲线集合中提取每个样品的数条热史曲线；

在同一时间点上，对提取后每个样品中所有热史曲线上的每个数据求取平均，得到每个数据的平均值；

按时间顺序将所有数据的平均值连接成曲线，该曲线作为该样品的最终热史模拟结果。其中，将N个样品的最终热史模拟结果作为该待模拟地质的热史模拟结果。

需要说明的是，主要是通过步骤S1至步骤S4，得到最终优度阈值，通过最终优度阈值从M组热史曲线集合中进行筛选，得到筛选后热史曲线组成的M组第二热史曲线集合，之后从M组第二热史曲线集合提取N个样品对应数条热史曲线，之后在同一时间点上，对每个样品的数条热史曲线上对应的数据求取平均，得到每个数据的平均值；按时间顺序将所有数据的平均值连接成曲线，该曲线作为该样品的最终热史模拟结果。其中，待模拟地质在垂直剖面上所有样品的最终热史模拟结果作为待模拟地质输出的热史模拟结果。

在本发明实施例中，该基于设置拟合优度阈值的热史模拟方法通过步骤S1和步骤S2计算得到第一拟合优度，将M个第一拟合优度分别与优度阈值对比，当修改不同数值的优度阈值时，不需要重新计算每组热史曲线的拟合优度，直接通过优度阈值与得到的M个第一拟合优度进行对比，筛选出符合优度阈值的热史曲线，以使该基于设置拟合优度阈值的热史模拟方法简单化。

图3为本发明实施例的基于设置拟合优度阈值的热史模拟方法一组热史曲线集合的示意图，图4为本发明实施例的基于设置拟合优度阈值的热史模拟方法拟合优度的分布直方图，图5为本发明实施例的基于设置拟合优度阈值的热史模拟方法拟合优度的散点图。

在本发明实施例中，以钻井取垂直剖面且在垂直剖面上获取6个样品为案例，对每个样品进行磷灰石U-Th/He定年及其中5个样品进行磷灰石裂变径迹年龄测定、其中又测试获得3个样品的磷灰石裂变径迹围限径迹长度分布数据，然后利用这些数据进行综合评价，进行热史模拟，如图3所示，该基于设置拟合优度阈值的热史模拟方法通过获取一组10000条热史曲线组成的热史曲线集合。按照通常设置的拟合优度阀值(0.05)而言，没有一组热史曲线其拟合优度会大于这个阈值，也就是说没有热史曲线为可以被保留，那么此实验测试很可能就得不到热史模拟结果了，这种现象经常会使得地质学家不理解，实测数据越多，满足拟合优度的难度就越大，看上去反而变得更加不准确。显然出现这种矛盾的原因就在于，对拟合优度的分布不了解，过于迷信拟合优度的通用阀值(例如0.05)，没有能够根据实际情况灵活调整阀值，因此才得不到合适的模拟结果。对此我们通过绘制拟合优度的直方图(图4)和散点图(图5)，可以看到所有热史曲线组的拟合优度(GOF)的分布，在0.01～0.04之间仍然分布有10个左右的拟合优度值，虽然没有0.05这个拟合优度值高，但仍然具有参考价值，因此根据拟合优度分布图，重新设置了拟合优度的阀值，即0.01。最终取>优度阈值为0.01的热史曲线，将其平均获得第一热史模拟结果。与该样品实际所在地区地质演化相对比，其热史模拟结果与实际是比较符合的。图3中虚线分别代表顶、底样品的温度的上限及下限，图中框架表示根据地质背景获得的大致的热史曲线拐点范围，有助于约束随机搜索的范围并提高模拟可靠性。如图4和图5所示，通过拟合优度的直方图及散点图，可以更改不同的优度阀值来获取更多或者较少的热史曲线数量，进而获得更符合用户实际需要的热史模拟结果。

实施例二：

图6为本发明实施例的基于设置拟合优度阈值的热史模拟系统的框架图。

如图6所示，本发明实施例还提供一种基于设置拟合优度阈值的热史模拟系统，包括数据获取模块10、拟合优度计算模块20、筛选计算模块30和第一输出结果模块40；

数据获取模块10，用于获取待模拟地质在垂直剖面上M组第一热史曲线集合；

拟合优度计算模块20，用于采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个第一拟合优度；

筛选计算模块30，用于将M个第一拟合优度与优度阈值对比，筛选得到K个第一拟合优度；并对K个第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果；

第一输出结果模块40，用于设置不同的优度阈值再通过筛选计算模块得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果；

其中，K小于M，K、M、A、A均为大于5的自然数，拟合优度公式为：

式中，GOF为第一拟合优度，age为热史曲线的模拟年龄，

为垂直剖面上所有样品实验测试年龄的平均误差，

为垂直剖面上样品模拟年龄与实验年龄之间差值的绝对值的平均。

在本发明实施例中，数据获取模块10包括第一样品获取子模块、第二样品获取子模块、第一数据获取子模块和第二数据获取子模块；

第一样品获取子模块，用于获取模拟地质在垂直剖面上顶端的样品作为首样品，首样品包括M条热史曲线；

第二样品获取子模块，用于以首样品的高度为基准，通过古地温梯度函数从模拟地质的垂直剖面上获取N-1个样品以及与每个样品对应的温差；

第一数据获取子模块，用于将首样品的M条热史曲线通过温差平移，得到对应样品的M条热史曲线；

第二数据获取子模块，用于从N个样品提取同一纵列的热史曲线组成一组第一热史曲线集合，每个样品包括M条热史曲线，则组成M组第一热史曲线集合。

在本发明实施例中，该基于设置拟合优度阈值的热史模拟系统包括第二输出结果模块50，第二输出结果模块50用于获取与第二热史模拟结果对应的优度阈值作为最终优度阈值，通过最终优度阈值对M组第一热史曲线集合的热史曲线进行筛选，得到M组第二热史曲线集合，从M组第二热史曲线集合中提取每个样品的数条热史曲线，在同一时间点上，对提取后每个样品中所有热史曲线上的每个数据求取平均，得到每个数据的平均值，按时间顺序将所有数据的平均值连接成曲线，该曲线作为该样品的最终热史模拟结果。

需要说明的是，实施例二系统中的模块对应于实施例一方法中的步骤，实施例一方法的步骤内容已在实施例一中详细阐述了，在此实施例二中不再对模块的内容进行详细阐述。

实施例三：

本发明实施例提供了一种基于设置拟合优度阈值的热史模拟设备，包括处理器以及存储器；

存储器，用于存储程序代码，并将程序代码传输给处理器；

处理器，用于根据程序代码中的指令执行上述的基于设置拟合优度阈值的热史模拟方法。

需要说明的是，处理器用于根据所程序代码中的指令执行上述的一种基于设置拟合优度阈值的热史模拟方法实施例中的步骤。或者，处理器执行计算机程序时实现上述各系统/装置实施例中各模块/单元的功能。

示例性的，计算机程序可以被分割成一个或多个模块/单元，一个或者多个模块/单元被存储在存储器中，并由处理器执行，以完成本申请。一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述计算机程序在终端设备中的执行过程。

终端设备可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。终端设备可包括，但不仅限于，处理器、存储器。本领域技术人员可以理解，并不构成对终端设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如终端设备还可以包括输入输出设备、网络接入设备、总线等。

所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

存储器可以是终端设备的内部存储单元，例如终端设备的硬盘或内存。存储器也可以是终端设备的外部存储设备，例如终端设备上配备的插接式硬盘，智能存储卡(SmartMedia Card，SMC)，安全数字(Secure Digital，SD)卡，闪存卡(Flash Card)等。进一步地，存储器还可以既包括终端设备的内部存储单元也包括外部存储设备。存储器用于存储计算机程序以及终端设备所需的其他程序和数据。存储器还可以用于暂时地存储已经输出或者将要输出的数据。

所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种基于设置拟合优度阈值的热史模拟方法，其特征在于，包括以下步骤：

S1.获取待模拟地质在垂直剖面上M组第一热史曲线集合；

S2.采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个所述第一拟合优度；

S3.将M个所述第一拟合优度与优度阈值对比，筛选得到K个所述第一拟合优度；并对K个所述第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果；

S4.设置不同的优度阈值再按照步骤S3得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果；

其中，K小于M，K、M、A均为大于5的自然数，所述拟合优度公式为：

式中，GOF为第一拟合优度，age为热史曲线的模拟年龄，

为垂直剖面上所有样品实验测试年龄的平均误差，

为垂直剖面上样品模拟年龄与实验年龄之间差值的绝对值的平均。

2.根据权利要求1所述的基于设置拟合优度阈值的热史模拟方法，其特征在于，获取待模拟地质在垂直剖面上M组第一热史曲线集合的步骤包括：

获取模拟地质在垂直剖面上顶端的样品作为首样品，所述首样品包括M条热史曲线；

以所述首样品的高度为基准，通过古地温梯度函数从模拟地质的垂直剖面上获取N-1个样品以及与每个所述样品对应的温差；

将所述首样品的M条热史曲线通过所述温差平移，得到对应所述样品的M条热史曲线；

从N个所述样品提取同一纵列的热史曲线组成一组第一热史曲线集合，每个所述样品包括M条热史曲线，则组成M组第一热史曲线集合。

3.根据权利要求1所述的基于设置拟合优度阈值的热史模拟方法，其特征在于，得到每组第一热史曲线集合的第一拟合优度包括：对每组第一热史曲线集合通过U-Th/He、裂变径迹年龄、裂变径迹长度和镜质体反射率Ro方法进行计算，得到对应的每组第一热史曲线集合的拟合优度，并从中选取拟合优度数值最小的作为该组第一热史曲线集合的第一拟合优度。

4.根据权利要求2所述的基于设置拟合优度阈值的热史模拟方法，其特征在于，包括：

获取与所述第二热史模拟结果对应的优度阈值作为最终优度阈值；

通过所述最终优度阈值对M组第一热史曲线集合的热史曲线进行筛选，得到M组第二热史曲线集合；

从M组所述第二热史曲线集合中提取每个样品的数条热史曲线；

在同一时间点上，对提取后每个样品中所有热史曲线上的每个数据求取平均，得到每个数据的平均值；

按时间顺序将所有所述数据的平均值连接成曲线，该曲线作为该样品的最终热史模拟结果。

5.根据权利要求4所述的基于设置拟合优度阈值的热史模拟方法，其特征在于，包括：将N个所述样品的最终热史模拟结果作为该待模拟地质的热史模拟结果。

6.根据权利要求1所述的基于设置拟合优度阈值的热史模拟方法，其特征在于，包括：获取待模拟地质在垂直剖面上至少10000组第一热史曲线集合。

7.一种基于设置拟合优度阈值的热史模拟系统，其特征在于，包括数据获取模块、拟合优度计算模块、筛选计算模块和第一输出结果模块；

所述数据获取模块，用于获取待模拟地质在垂直剖面上M组第一热史曲线集合；

所述拟合优度计算模块，用于采用拟合优度公式对每组第一热史曲线集合计算，得到每组第一热史曲线集合的第一拟合优度，则M组第一热史曲线集合对应得到M个所述第一拟合优度；

所述筛选计算模块，用于将M个所述第一拟合优度与优度阈值对比，筛选得到K个所述第一拟合优度；并对K个所述第一拟合优度求取均值，得到待模拟地质在垂直剖面上的第一热史模拟结果；

所述第一输出结果模块，用于设置不同的优度阈值再通过所述筛选计算模块得到A个第一热史模拟结果，从A个第一热史模拟结果筛选出最为合适的第一热史模拟结果作为待模拟地质在垂直剖面上的第二热史模拟结果；

其中，K小于M，K、M、A均为大于5的自然数，所述拟合优度公式为：

式中，GOF为第一拟合优度，age为热史曲线的模拟年龄，

为垂直剖面上所有样品实验测试年龄的平均误差，

为垂直剖面上样品模拟年龄与实验年龄之间差值的绝对值的平均。

8.根据权利要求7所述的基于设置拟合优度阈值的热史模拟系统，其特征在于，所述数据获取模块包括第一样品获取子模块、第二样品获取子模块、第一数据获取子模块和第二数据获取子模块；

所述第一样品获取子模块，用于获取模拟地质在垂直剖面上顶端的样品作为首样品，所述首样品包括M条热史曲线；

所述第二样品获取子模块，用于以所述首样品的高度为基准，通过古地温梯度函数从模拟地质的垂直剖面上获取N-1个样品以及与每个所述样品对应的温差；

所述第一数据获取子模块，用于将所述首样品的M条热史曲线通过所述温差平移，得到对应所述样品的M条热史曲线；

所述第二数据获取子模块，用于从N个所述样品提取同一纵列的热史曲线组成一组第一热史曲线集合，每个所述样品包括M条热史曲线，则组成M组第一热史曲线集合。

9.根据权利要求8所述的基于设置拟合优度阈值的热史模拟系统，其特征在于，包括第二输出结果模块，所述第二输出结果模块用于获取与所述第二热史模拟结果对应的优度阈值作为最终优度阈值，通过所述最终优度阈值对M组第一热史曲线集合的热史曲线进行筛选，得到M组第二热史曲线集合，从M组所述第二热史曲线集合中提取每个样品的数条热史曲线，在同一时间点上，对提取后每个样品中所有热史曲线上的每个数据求取平均，得到每个数据的平均值，按时间顺序将所有所述数据的平均值连接成曲线，该曲线作为该样品的最终热史模拟结果。

10.一种基于设置拟合优度阈值的热史模拟设备，其特征在于，包括处理器以及存储器；

所述存储器，用于存储程序代码，并将所述程序代码传输给所述处理器；

所述处理器，用于根据所述程序代码中的指令执行如权利要求1-6任意一项所述的基于设置拟合优度阈值的热史模拟方法。

Images