一种基于数字地质模型的矿床预测方法技术

本发明专利技术属于成矿预测技术领域，具体涉及一种基于数字地质模型的矿床预测方法，步骤一：根据需要建立研究区数字地质模型；根据研究区实际情况建立数字地质模型，包括：构造模型、矿体模型和岩体模型；步骤二：根据数字地质模型，建立块体模型和各预测要素块体模型；步骤三：计算出证据权法的后验概率值和信息量值，并把证据权法的后验概率值和信息量值结合起来，进行分析，求出有利区间；步骤四：将赋予块体模型的后验概率值和信息量值，筛选出有利区间，在次基础上，结合地质认识，圈定靶区。圈定靶区。圈定靶区。

【技术实现步骤摘要】
一种基于数字地质模型的矿床预测方法


[0001]本专利技术属于成矿预测
，具体涉及一种基于数字地质模型的矿床预测方法。

技术介绍

[0002]进入21世纪，越来越多的矿床进行了三维地质模型的构建，得以直观查看各种预测要素与矿体的空间分布关系，以便准确地预测已知矿体外围和深部延伸趋势。同时矿产资源预测与评价也经历了从二维GIS技术支撑的预测评价向三维成矿预测转换的过程。三维成矿预测可采用的数学模型很多，如证据权模型和信息量模型，但是，现有技术中无论单一采用哪种数学模型，都会存在局限性，因此，为了避免产生这种局限，需要设计一种采用多种数学模型综合的矿床预测方法，将证据权模型和信息量模型综合起来，求出有利区间，利用GOCAD软件，将证据权值和信息量值赋予块体模型，以此筛选出有利区间，圈定靶区，来解决现有技术中采用单一数学模型进行矿床预测局限性的技术问题。

技术实现思路

[0003]本专利技术针对现有技术中,提出一种基于数字地质模型的矿床预测方法，用于解决采用单一数学模型进行矿床预测具有局限性的技术问题。
[0004]本专利技术的技术方案如下所述：
[0005]一种基于数字地质模型的矿床预测方法，包括以下步骤：
[0006]步骤一：根据需要建立研究区数字地质模型；根据研究区实际情况建立数字地质模型，包括：构造模型、矿体模型和岩体模型；
[0007]步骤二：根据数字地质模型，建立块体模型和各预测要素块体模型；
[0008]步骤三：计算出证据权法的后验概率值和信息量值，并把证据权法的后验概率值和信息量值结合起来，进行分析，求出有利区间；
[0009]步骤3.1：计算证据权法的后验概率值；
[0010]步骤3.2：计算三维信息量值；
[0011]步骤3.3：将证据权法的后验概率值和信息量值结合分析，求出有利区间；
[0012]步骤3.4：分析证据权法的后验概率值，求出后验概率值的有利区间；根据步骤3.1证据权法计算的后验概率值，统计各区间中含已知矿体的矿块数；
[0013]步骤3.5：分析信息量值，求出信息量值的有利区间；
[0014]步骤3.6：综合分析证据权法的后验概率值和信息量值，求出最终有利区间；
[0015]步骤四：将赋予块体模型的后验概率值和信息量值，筛选出有利区间，在次基础上，结合地质认识，圈定靶区。
[0016]所述步骤3.1：计算证据权法的后验概率值，还包括：
[0017]在GOCAD软件中对已约束的块体进行属性赋值；
[0018]计算先验概率：根据已知矿体的分布，统计与分析各预测要素与矿体在空间上的
关系，计算各块体内各证据因子的先验概率；
[0019]计算权重：结合证据权公式(1)、公式(2)，进一步计算该研究区的权重值；
[0020][0021][0022]其中，W+代表各证据因子在预测区域中存在时的权重值，W-代表各证据因子在预测区域中不存在时的权重值，B为有因子存在的单元数，为因子
[0023]缺失的单元数，D为含矿单元数，为不含矿点的单元数。
[0024]预测要素证据因子层与已知矿床的证据因子层的相关程度的计算公式如下公式(3)：
[0025]C＝W
+-W-............(3)
[0026]其中，C表示预测要素证据因子层与已知矿床的证据因子层的相关程度，则C；当C大于0时，C的值与该证据因子的成矿有利度成正相关，当C为0时，说明该证据因子与成矿没有关系；当C小于0时，C的值与该证据因子的成矿有利度成负相关；
[0027]计算后验概率：在先验概率的基础上，通过公式(4)计算得到后验概率；
[0028][0029]其中：L(D/B)＝L(D)+W
+
，后验概率值越大，矿体成矿越有利；
[0030]所述步骤3.2：计算三维信息量值，包括：通过收集各类地质要素单变量在研究区的空间分布情况；对研究区预测要素定量化分析、统计，将矿床各预测要素的属性值进行二值化处理，通过三维信息量方法计算得出矿床区域内各找矿要素的信息量值；
[0031]计算成矿要素信息量值的条件概率计算公式如下公式(5)：
[0032][0033]公式IA(B)表示找矿标志A中含有矿体B的信息量值的大小；P(A)表示研究范围内出现找矿标志A的概率值；P(A/B)表示在B矿存在时找矿标志A出现的概率值；
[0034]通过频率值进行估算，计算结果同概率值计算结果等同；故在此使用频率值代替概率值进行估算，计算公式如下公式(6)：
[0035][0036]其中，Nj表示同时包含找矿标志A与矿体的单元块数；N表示含矿单元块总数；Sj表示具有找矿标志A的单元块数；S表示所含单元块总数。
[0037]所述步骤3.5，分析信息量值，求出信息量值的有利区间包括：根据研究区预测要素信息量分析各个预测要素对成矿的影响；按照信息量值大小，结合研究区地质背景和成矿规律划分信息量值区，综合统计信息量区间内含矿单元数。
[0038]所述步骤3.6，综合分析证据权法的后验概率值和信息量值，求出最终有利区间，包括：在后验概率区间下，结合信息量值研究结果，计算不同信息量区间中的矿块比例、块体比例以及含矿比率；
[0039]所述矿块比例为：信息量区间内含已知矿块数量/总已知矿块数量；
[0040]所述块体比例为：信息量区间内块体总数量/研究区块体总数量；
[0041]所述含矿比率为：矿块比例/块体比例。
[0042]本专利技术的有益效果为：
[0043]本专利技术设计提出一种基于数字地质模型的矿床预测方法，将计算出的证据权法中的后验概率值和信息量法中的信息量值通过数据处理后，将两个值结合分析，求出有利区间，在次基础上进行预测，圈定靶区，提高了预测的准确性和可视性。本专利技术解决了在矿床预测过程中，使用单一方法的局限性，利用证据权法的后验概率和信息量，两种数值，结合分析，求出有利区间，缩小并三维展示预测范围，在此基础上再圈定靶区，提高预测准确性和可视性。
附图说明
[0044]图1为本专利技术设计的一种基于数字地质模型的矿床预测方法的流程框图；
[0045]图2为本专利技术实施例所述的研究区数字地质模型；
[0046]图3为本专利技术实施例所述的研究区后验概率含已知矿块统计图；
[0047]图4为本专利技术实施例所述研究区信息量值叠合已知矿体统计结果图；
[0048]图5为本专利技术实施例所述的研究区信息量统计分析图；
[0049]图6为本专利技术实施例所述的靶区分布位置图。
具体实施方式
[0050]下面结合附图和实施例对本专利技术的一种基于数字地质模型的矿床预测方法进行详细说明。
[0051]一种基于数字地质模型的矿床预测方法，包括以下步骤：
[0052]步骤一：根据需要建立研究区数字地质模型；根据研究区实际情况建立数字地质模型，包括：本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于数字地质模型的矿床预测方法，其特征在于：包括以下步骤：步骤一：根据需要建立研究区数字地质模型；根据研究区实际情况建立数字地质模型，包括：构造模型、矿体模型和岩体模型；步骤二：根据数字地质模型，建立块体模型和各预测要素块体模型；步骤三：计算出证据权法的后验概率值和信息量值，并把证据权法的后验概率值和信息量值结合起来，进行分析，求出有利区间；步骤3.1：计算证据权法的后验概率值；步骤3.2：计算三维信息量值；步骤3.3：将证据权法的后验概率值和信息量值结合分析，求出有利区间；步骤3.4：分析证据权法的后验概率值，求出后验概率值的有利区间；根据步骤3.1证据权法计算的后验概率值，统计各区间中含已知矿体的矿块数；步骤3.5：分析信息量值，求出信息量值的有利区间；步骤3.6：综合分析证据权法的后验概率值和信息量值，求出最终有利区间；步骤四：将赋予块体模型的后验概率值和信息量值，筛选出有利区间，在次基础上，结合地质认识，圈定靶区。2.根据权利要求1所述的一种基于数字地质模型的矿床预测方法，其特征在于：所述步骤3.1：计算证据权法的后验概率值，还包括：在GOCAD软件中对已约束的块体进行属性赋值；计算先验概率：根据已知矿体的分布，统计与分析各预测要素与矿体在空间上的关系，计算各块体内各证据因子的先验概率；计算权重：结合证据权公式(1)、公式(2)，进一步计算该研究区的权重值；计算权重：结合证据权公式(1)、公式(2)，进一步计算该研究区的权重值；其中，W+代表各证据因子在预测区域中存在时的权重值，W-代表各证据因子在预测区域中不存在时的权重值，B为有因子存在的单元数，为因子缺失的单元数，D为含矿单元数，为不含矿点的单元数。预测要素证据因子层与已知矿床的证据因子层的相关程度的计算公式如下公式(3)：C＝W
+-W-…………
(3)其中，C表示预测要素证据因子层与已知矿床的证据因子层的相关程度，则C；当C大于0时，C的值与该证据因子的成矿有利度成正相关，当C为0时，说明该证据...

【专利技术属性】
技术研发人员：白芸，朱鹏飞，李晓翠，孔维豪，曹珂，孙璐，刘琳莹，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：