CN117251739B - 用于砂岩型铀资源定量预测的样本集的构建方法 - Google Patents

Abstract

本发明实施例涉及铀资源预测技术领域，具体涉及一种用于砂岩型铀资源定量预测的样本集的构建方法。样本集构建方法包括如下步骤：获取砂岩型铀矿研究区中的多个地质类图件信息，其中，每个地质类图件信息包含与空间相关的具有地质属性的数据信息；获取砂岩型铀矿研究区中的钻孔信息，其中，钻孔信息包含钻孔所在位置以及是否有矿；基于地质类图件信息，确定钻孔所在位置对应的地质属性的数据信息，将钻孔对应的地质属性的数据信息以及是否有矿，作为一组样本。本发明实施例提出的样本集构建方法，有利于在基于机器学习的砂岩型铀矿资源预测工作中提供科学、有效的样本数据。

Description

用于砂岩型铀资源定量预测的样本集的构建方法

技术领域

本发明涉及铀资源预测技术领域，具体涉及一种用于砂岩型铀资源定量预测的样本集的构建方法。

背景技术

机器学习是大数据分析和知识发现切实有效的工具，并将有助于砂岩型铀矿资源预测评价。借助地理信息系统制图技术（GIS）技术和机器学习算法，提取地质、物探、化探、遥感等多源信息，并对其进行集成研究，能够挖掘、识别隐性的和深层次的找矿信息。进行砂岩型铀矿资源预测评价的首要工作就是建立科学有效的样本集。

发明内容

为此，本发明实施例提出了一种用于砂岩型铀资源定量预测的样本集的构建方法。构建方法包括如下步骤：获取砂岩型铀矿研究区中的多个地质类图件信息，其中，每个地质类图件信息包含与空间相关的具有地质属性的数据信息；获取砂岩型铀矿研究区中的钻孔信息，其中，钻孔信息包含钻孔所在位置以及是否有矿；基于地质类图件信息，确定钻孔所在位置对应的地质属性的数据信息，将钻孔对应的地质属性的数据信息以及是否有矿，作为一组样本。

本发明实施例提出的样本集构建方法，有利于在基于机器学习的砂岩型铀矿资源预测工作中提供科学、有效的样本数据。

附图说明

图1示出了本申请实施例中构建样本集方法的流程图；

图2示出了本申请实施例中确定钻孔所在位置对应的地质属性的数据信息的流程图；

图3示出了本申请一个实施例中未修正的等值线数据图；

图4示出了本申请实施例中等值线图数据转化的离散点数据的示意图；

图5示出了将图4所示离散点数据进行插值处理得到的覆盖砂岩型铀矿研究区的连续性数据；

图6示出了本申请实施例中的样本集构建的示意性原理图。

具体实施方式

为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请实施例的附图，对本申请的技术方案进行清楚、完整地描述。显然，所描述的实施例是本申请的一个实施例，而不是全部的实施例。基于所描述的本申请的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。

需要说明的是，除非另外定义，本申请使用的技术术语或者科学术语应当为本申请所属领域内具有一般技能的人士所理解的通常意义。若全文中涉及“第一”、“第二”等描述，则该“第一”、“第二”等描述仅用于区别类似的对象，而不能理解为指示或暗示其相对重要性、先后次序或者隐含指明所指示的技术特征的数量，应该理解为“第一”、“第二”等描述的数据在适当情况下可以互换。若全文中出现“和/或”，其含义为包括三个并列方案，以“A和/或B”为例，包括A方案，或B方案，或A和B同时满足的方案。

本申请实施例提供一种用于砂岩型铀资源定量预测的样本集的构建方法。参见图1，本申请实施例的构建方法可以包括：步骤S100至步骤S300。

步骤S100：获取砂岩型铀矿研究区中的多个地质类图件信息，其中，每个地质类图件信息包含与空间相关的具有地质属性的数据信息。

步骤S200：获取砂岩型铀矿研究区中的钻孔信息，其中，钻孔信息包含钻孔所在位置以及是否有矿。

步骤S300：基于地质类图件信息，确定钻孔所在位置对应的地质属性的数据信息，将钻孔对应的地质属性的数据信息以及是否有矿，作为一组样本。

容易理解，在步骤S200中的钻孔信息中包含了钻孔所在位置以及从该钻孔中获取的岩心是否成矿的信息。当该钻孔中获取的岩心成矿时，对应的钻孔信息为有矿；当该钻孔中获取的岩心未成矿时，对应的钻孔信息为无矿。

在一些实施例中，本申请实施例的构建方法还包括：基于砂岩型铀矿研究区的地质信息、矿床的主成矿期、砂岩型铀矿成因以及成矿模式，确定参与样本集构建的多个地质属性。

可以根据研究区的地质背景和成矿背景分析，确定参与样本集构建的多个地质属性。例如，分析认为矿体是在某个成矿时期的古河道中产生的，可以将沉积相、成矿期地层埋深等作为参与样本集构建的地质属性，可以获得对应沉积相图层、成矿期地层埋深图层等。

本申请的发明人发现，砂岩型铀矿最主要的几个地质属性是砂体厚度、含砂率、灰砂率以及沉积相。

在一些实施例中，地质属性可以包括砂体厚度、含砂率、灰砂率中的至少一种。本申请的发明人发现，当地质属性包含砂体厚度、含砂率、灰砂率中的至少一种时，构建的样本集的有效性较高。

在本申请的一个实施例中，地质属性包括砂体厚度、含砂率、灰砂率以及沉积相。本申请的发明人发现，这样构建的样本集的有效性更高。

在一些实施例中，地质类图件可以为包含地质属性的等值线数据的砂岩铀矿地质类图件和/或包含离散型数据的砂岩铀矿地质类图件。

图6示出了本申请实施例中的样本集构建的示意性原理图，图中示出了4个地质类图件信息，4个地质类图件信息分别对应含砂率图层20、灰砂率图层30、砂体厚度图层40以及沉积相图层50，4个地质类图件信息对应的地质属性分别为含砂率、灰砂率、砂体厚度以及沉积相。

其中，含砂率图层20、灰砂率图层30、砂体厚度图层40分别为包含等值线数据的砂岩铀矿地质类图件，即图层中分别包括与空间相关的等值线信息。在这些图层中，位于不同空间位置处（如对应不同坐标）的地质属性的数值相同的点共同组成等值线。

沉积相图层50属于包含离散型数据的砂岩铀矿地质类图件。对于这一类地质类图件，图层中不同空间位置对应的相同地质属性对应的名称种类不同（不同的名称种类对应沉积相图层50中的不同颜色）。例如，对于沉积相地质属性，其可能的名称种类包括河道、心滩、泛滥平原等。可以为离散型的地质属性的不同名称种类设定对应的数值属性，例如将河道沉积相设定为1，心滩沉积相设定为2等。

在本申请实施例中，对于地质类图件信息，可以利用格网对其图层进行划分。如此，基于构建的格网，可以获得每个格网对应的地质属性的数据信息。这样，对于具有钻孔的单个格网，则可以获得该格网对应的对应地质属性的数据信息。

在一些实施例中，当地质类图件为包含地质属性的等值线数据的砂岩铀矿地质类图件时，可以根据等值线数据得到覆盖砂岩型铀矿研究区的连续性数据（即数字表面模型数据图层），以便于在将图层进行格网化时，能够保证每个格网中具有对应地质属性的数值。

本申请实施例通过将等值线图数据转换为离散点数据，再将离散点数据进行插值处理，能够得到覆盖砂岩型铀矿研究区的连续性数据。

具体地，将等值线图数据转换为离散点数据的示意图可以参见图4；对离散点数据进行插值处理，得到覆盖砂岩型铀矿研究区的连续性数据可以参见图5。

容易理解，每条等值线上的点除了具有等值线对应的数值外，还具有与数值对应的线属性。以砂体厚度图层为例，其每条砂体厚度等值线上的各点具有砂体厚度的数值和砂体厚度线属性。例如不同砂体厚度等值线对应的数值可以为30、50、86，将等值线转为离散点后，这些离散点的属性与等值线的属性相同，均为砂体厚度。

参见图2，在步骤S300中，基于地质类图件信息，确定钻孔所在位置对应的地质属性的数据信息，包括步骤S310至步骤S340。

步骤S310：基于地质类图件信息，获取地质属性的等值线图数据。

步骤S320：将等值线图数据转换为离散点数据。

步骤S330：对离散点数据进行插值处理，得到覆盖砂岩型铀矿研究区的连续性数据图层。

步骤S340：基于连续性数据图层，确定钻孔所在位置对应的地质属性的数据信息。

在一些实施例中，在步骤S320将等值线图数据转换为离散点数据之前，还包括：对等值线图数据进行修正。

参见图3，某些情况下，地质类图件信息中的地质属性等值线图数据可能存在重叠、断点、曲线不圆滑等问题。本申请的发明人发现，这些问题会导致构建的样本集的有效性降低，因此，需要将等值线图进行修正。

可以通过地理信息系统制图技术（GIS）技术，根据空间拓扑学检查对断线、重复区域进行修正，对不圆滑曲线进行点加密，从而得到圆滑曲线的等值线要素图。

具体地，修正过程包括将等值线图数据进行可视化处理，人工查看数据质量，并通过GIS工具检查是否存在重叠、断点区域，如果存在则通过软件进行修复。此外，还可以对不圆滑的曲线进行加密点处理，使其变得圆滑，从而在后续步骤中将等值线转化为离散点时，圆滑的曲线转化成的离散点的总体形态更能符合等值线的形态。

在一些实施例中，在步骤S330中，可以采用克里金法对离散点数据进行插值处理对离散点数据进行插值处理，以得到覆盖砂岩型铀矿研究区的连续性数据。

获得离散点数据后，图中等值线之间的部分仍然存在空白区域，需要通过插值处理在整个研究区范围内进行插值计算，为空白区域添加数据，从而获得覆盖砂岩型铀矿研究区的连续性数据图层。该连续性数据图层可以为GIS中的数字表面模型（DSM）。在该模型中，整个表面按照一定的格网尺寸做成覆盖整个表面的数据，任何一个位置都有一个数值。

克里金法是一种地质学中常用的插值方法，该方法以区域化变量为基础，以变异函数为基本工具，对未知样点进行满足地质认识的最优化估计。图5示出了本申请一个实施例中进行插值处理后的预测区数字表面模型。克里金法插值后得到的结果整体趋势与地质认识一致，但是不会严格地遵循数据的绝对值。用克里金法进行插值，相比于其他的插值方法，准确性更高。

在一些实施例中，在步骤S340中，基于连续性数据图层，确定钻孔所在位置对应的地质属性的数据信息，包括：基于钻孔信息构建格网图层；基于连续性数据图层和钻孔信息，根据空间位置的匹配，获得每个钻孔所在位置对应的格网以及格网对应的地质属性的数据信息，从而确定钻孔所在位置对应的地质属性的数据信息。

在一些实施例中，构建的格网图层仅包含序列号ID与基本图形要素，该基本图形要素为格网，基本图形要素本身包含了空间位置信息。其中，序列号ID用于每一个样本进行编号，以便对应。

可以基于空间位置将数据图层的数据连接到该格网图层中。具体地，根据空间位置的匹配将数据图层对应的地质属性添加到格网图层的图层属性中，并将数据图层的属性值添加到格网图层对应的属性列中，以形成新的格网图层。此时，该新格网图层，不仅包含原有的序列号ID和基本图形要素，还包括其他数据图层的属性值，作为其属性1、属性2等，以此类推。

可以基于空间位置的匹配，将每个钻孔信息添加到新的格网图层的目标列中。根据该新的格网图层，即可确定样本集。

可以将钻孔信息进行格网化，创建的格网图层中包含无钻孔信息格网和有钻孔信息的格网，其中有钻孔信息格网对应钻孔有矿或钻孔无矿，将有钻孔信息的格网的数据作为样本集进行数据提取。由此可见，在本申请实施例中，研究区内钻孔的数量决定了样本集的数量。

具体地，获得监督学习所需要的正、负样本，需要对本区域内的钻孔信息按照是否含矿，设定目标属性值，与其他属性值共同作为样本标识。为了获得用于砂岩型铀资源定量预测所需要的正、负样本，需要对研究区内的钻孔信息按照是否含矿，划分为0（无矿孔）和1（有矿孔），同样按照空间位置匹配方式，添加到新格网图层的目标属性一列，作为监督学习所需要的正、负样本标识，参见表1。

表1 新的格网图层的格网数据

如此，基于构建的新的格网图层，可以获得每个格网对应的不同地质属性的数据信息。而对于具有钻孔的单个格网，则可以获得该格网对应的不同地质属性的数据信息，从而获得一组样本。从而，根据该新格网图层，能够获得用于砂岩型铀资源定量预测的样本集。

在一些实施例中，在构建格网图层时，需要确定单个格网的尺寸。单个格网的尺寸可以根据经验确定。为了保证样本数据的准确有效，需要保证每个格网中只包含一个钻孔。因此，在本申请的一些实施例中，可以基于钻孔信息确定格网图层中单个格网的尺寸。

在一些实施例中，基于钻孔信息确定单个格网的尺寸，具体可以包括：基于钻孔信息确定两个钻孔之间的最小间距；根据最小间距确定单个格网的尺寸。

单个格网的尺寸可以与两个钻孔之间的最小间距基本相同。“基本相同”可以理解为两者之差不超过两个钻孔之间的最小间距的10%。

容易理解，确定钻孔最小间距，是为了确定“单个格网的尺寸是多少”才能确保“每个网格只含有1个钻孔”。两个钻孔之间的最小间距可以通过GIS的空间分析工具确定，之后要看格网划分后，是否实现了“每个格网只含有1个钻孔”。可以人工手动调整格网的大小，在确保“每个网格只含有1个钻孔”的前提下，出于减少计算量的考虑，网格可以尽量大一些。

例如，当基于钻孔信息确定相邻两个钻孔之间的最小间距为198米时，单个格网的尺寸可以为198米*198米，或者，也可以使单个格网的尺寸为200米*200米。根据最小间距确定单个格网的尺寸（或称为格网单元的尺寸），能够保证每个格网中只包含一个钻孔。

在一些实施例中，在步骤S320中，将等值线图数据转换为离散点数据，包括：根据单个格网的尺寸，确定从等值线图数据中提取的两个相邻点数据的间隔。

选取离散点时，两个相邻点数据的间隔可以根据构建的格网的单个格网的尺寸确定。例如，在一些实施例中，两个相邻点数据的间隔可以为单个格网尺寸的1/2~1/3。本申请的发明人发现，按照这样的方式，能够满足格网化后数据精度要求，有利于提高样本集的有效性。例如当单个格网的尺寸为200米时，确定从等值线图数据中提取的两个相邻点数据的间隔为70米~100米之间。

参见图6，在一些实施例中，例如，可以创建单个格网尺寸为200米*200米的格网图层。格网图层可以包括ID和基本图形要素。

基于空间位置的匹配，将含砂率图层20、灰砂率图层30、砂体厚度图层40以及沉积相图层50的地质属性添加到格网图层对应的图层属性表中，并将数据图层的属性值添加到格网图层对应的属性列中，以形成新的格网图层10。新的格网图层10对应的格网数据参见表1。新的格网图层10不仅包括原有的ID和基本图形要素，还包含从含砂率图层20、灰砂率图层30、砂体厚度图层40以及沉积相图层50添加的属性值。

参见图6，基于空间位置的匹配，含砂率图层20、灰砂率图层30、砂体厚度图层40以及沉积相图层50中，每个图层对应相同ID的格网的地质属性均能够添加到对应格网图层10相同位置的格网中。

容易理解，对于沉积相图层50的图层数据，其为离散型数据，可以为离散型的地质属性的不同名称种类设定对应的数值属性，例如将河道沉积相设定为1，心滩沉积相设定为2等。在基于空间位置的匹配时，可以将对应的对应数值属性赋值到相应的空间格网。

在一些实施例中，对于每组样本，基于钻孔是否有矿为每组样本添加标签，例如，对于钻孔有矿的样本，添加标签作为正样本，对于钻孔无矿的样本，添加标签作为负样本。

将得到的新格网数据转换为机器学习适用的表格类数据，以便于提高样本集数据适用范围，可以提供给任何一款软件或直接采用代码进行计算分析。将新格网数据中的ID和地质属性1、地质属性2等属性导出为包含n个地质属性列和m个数据记录的csv文件。最终完成用于构建机器学习的连续型数据样本集。

对于本发明的实施例，还需要说明的是，在不冲突的情况下，本发明的实施例及实施例中的特征可以相互组合以得到新的实施例。以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，本发明的保护范围应以权利要求的保护范围为准。

Claims (9)

1.一种用于砂岩型铀资源定量预测的样本集的构建方法，包括：

获取砂岩型铀矿研究区中的多个地质类图件信息，其中，每个所述地质类图件信息包含与空间相关的具有地质属性的数据信息；

获取所述砂岩型铀矿研究区中的钻孔信息，其中，所述钻孔信息包含钻孔所在位置以及是否有矿；

基于所述地质类图件信息，确定钻孔所在位置对应的地质属性的数据信息，将所述钻孔对应的地质属性的数据信息以及是否有矿，作为一组样本；

所述基于所述地质类图件信息，确定钻孔所在位置对应的地质属性的数据信息，包括：当所述地质类图件为包含地质属性的等值线数据的砂岩铀矿地质类图件时，基于所述地质类图件信息，获取所述地质属性的等值线图数据；将所述等值线图数据转换为离散点数据；对所述离散点数据进行插值处理，得到覆盖所述砂岩型铀矿研究区的连续性数据图层；基于所述连续性数据图层，确定钻孔所在位置对应的地质属性的数据信息；

所述基于所述连续性数据图层，确定钻孔所在位置对应的地质属性的数据信息，包括：基于所述钻孔信息构建格网图层；基于所述连续性数据图层和所述钻孔信息，根据空间位置的匹配，获得每个钻孔所在位置对应的格网以及所述格网对应的地质属性的数据信息，从而确定钻孔所在位置对应的地质属性的数据信息；

所述地质属性包括砂体厚度、含砂率、灰砂率以及沉积相；对应的所述地质类图件信息分别为含砂率图层、灰砂率图层、砂体厚度图层以及沉积相图层，其中，所述含砂率图层、所述灰砂率图层、所述砂体厚度图层分别为包含等值线数据的砂岩铀矿地质类图；所述沉积相图层中不同空间位置对应的相同地质属性对应的名称种类不同，所述方法包括：为不同名称种类设定对应的数值属性，在基于空间位置的匹配时，将对应的对应数值属性赋值到相应的空间格网。

2.根据权利要求1所述的构建方法，其中，所述基于所述钻孔信息构建格网图层，包括：

基于所述钻孔信息确定所述格网图层中单个格网的尺寸。

3.根据权利要求2所述的构建方法，其中，所述基于所述钻孔信息确定单个格网的尺寸，包括：

基于所述钻孔信息确定两个钻孔之间的最小间距；

根据所述最小间距确定单个格网的尺寸。

4.根据权利要求2所述的构建方法，其中，所述将所述等值线图数据转换为离散点数据，包括：

根据所述单个格网的尺寸，确定从所述等值线图数据中提取的两个相邻点数据的间隔。

5.根据权利要求4所述的构建方法，其中，以所述单个格网尺寸的1/2~1/3为间隔，从所述等值线图数据中提取点数据。

6.根据权利要求1所述的构建方法，其中，采用克里金法对所述离散点数据进行插值处理。

7.根据权利要求1所述的构建方法，还包括：

对于每组所述样本，基于钻孔是否有矿为每组所述样本添加标签。

8.根据权利要求1所述的构建方法，其中，在将所述等值线图数据转换为离散点数据之前，还包括：

对所述等值线图数据进行修正。

9.根据权利要求1所述的构建方法，还包括：

基于所述砂岩型铀矿研究区的地质信息、矿床的主成矿期、砂岩型铀矿成因以及成矿模式，确定参与所述样本集构建的多个地质属性。