多粒度矿体三维建模不确定性表达与传递方法.pdf

《多粒度矿体三维建模不确定性表达与传递方法.pdf》由会员分享，可在线阅读，更多相关《多粒度矿体三维建模不确定性表达与传递方法.pdf（12页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201911253301.4 (22)申请日 2019.12.09 (71)申请人 湖北省国土资源厅信息中心 地址 430000 湖北省武汉市武昌区公正路 27号 申请人 武汉蚁图时空科技有限公司 (72)发明人 李江徐江嬿徐凯 (74)专利代理机构 武汉知产时代知识产权代理 有限公司 42238 代理人 万文广 (51)Int.Cl. G06T 17/00(2006.01) (54)发明名称 一种多粒度矿体三维建模不确定性表达与 传递方法 (57)摘要 本发明提供了一种多粒度矿。

2、体三维建模不 确定性表达与传递方法， 该方法基于不确定性推 理网络和地质统计学克里格算法及序惯高斯模 拟算法， 通过对建模过程中因误差产生的推理型 不确定性和计算型不确定性的归纳、 总结和算法 解析， 结合矿体三维模型建立过程中的误差分布 与传递状况， 建立了多粒度矿体三维属性模型的 不确定性表达与传递模型框架， 实现了对矿体三 维模型不确定性的定量分析。 本发明的有益效果 是： 本发明所提出的技术方案实现了对矿体三维 模型不确定性的表达和传递过程的描述， 对矿山 企业利用三维模型进行设计、 规划和生产提供了 理论依据。 权利要求书3页 说明书7页 附图1页 CN 110930504 A 20。

3、20.03.27 CN 110930504 A 1.一种多粒度矿体三维建模不确定性表达与传递方法， 其特征在于： 包括以下步骤： S101： 获取多源矿山建模数据； S102： 对所述多源矿山建模数据产生的误差进行不确定性赋值， 形成矿山多源数据不 确定性集合， 进而对所述矿山多源数据不确定性集合内的不确定性赋值进行合成计算， 生 成建模数据可信度D_CF； S103： 对矿山建模过程误差进行可信度赋值， 形成结构建模过程不确定性集合， 进而对 所述结构建模过程不确定性集合内不确定性赋值进行合成计算， 生成结构建模过程可信度 SMP_CF； S104： 对所述建模数据可信度D_CF和所述建模过。

4、程数据可信度SMP_CF进行叠加计算， 得到结构模型可信度SM_CF； S105： 根据矿体数据基础和矿体三维模型需求的不同， 对建立三维属性模型按模拟建 模和估值建模两种方式进行多粒度特征建模并根据建模过程分别计算不确定性， 即克里格 估值可信度KC_CF和模拟计算可信度SC_CF； S106： 将所述结构模型可信度SM_CF分别与所述克里格估值可信度KC_CF和所述模拟计 算可信度SC_CF进行叠加计算， 得到多粒度矿山三维属性模型可信度AM_CF。 2.如权利要求1所述的一种多粒度矿体三维建模不确定性表达与传递方法， 其特征在 于： 步骤S101中， 所述多源矿山建模数据包括： 钻孔位置。

5、坐标、 钻孔深度、 钻孔样品段起点坐 标、 终点坐标、 钻孔样品长度和钻孔组合样长度； 所述多源矿山建模数据根据钻孔、 地质剖 面地质资料整理获得。 3.如权利要求1所述的一种多粒度矿体三维建模不确定性表达与传递方法， 其特征在 于： 步骤S102中， 所述多源矿山建模数据产生的误差包括矿山勘测过程中产生的误差； 对所 述多源矿山建模数据产生的误差进行不确定性赋值， 形成矿山多源数据不确定性集合， 进 而对所述矿山多源数据不确定性集合内的不确定性赋值进行合成计算， 生成建模数据可信 度D_CF； 具体包括： S201： 对不同类型数据的误差分别进行不确定性赋值； 其中， 对某一个类型的数据产生。

6、 的误差进行不确定性赋值的方法为： 将该类型的数据多次测量所得结果的平均值与被测量 的真值之差进行比较并转换为该类型数据的可信度CF(Ei)， 计算公式如公式(1)所示： 上式中， Ei代表第i个类型的矿山建模数据； CF(Ei)代表该类型数据的可信度值， 取值范 围为0,1； xj代表该类型数据的第j个测量值的真值； 为第j个测量值的多次测量平均值； i1,2,n， 且不同的取值代表不同类型的矿山建模数据； S202： 将各类型数据的可信度值CF(Ei)组成的集合作为矿山多源数据不确定性集合； 并 对所述矿山多源数据不确定性集合内的所有可信度值进行合成计算， 得到建模数据可信度 D_CF； 。

7、其中， 以单一可信度值的最小值为合成计算结果， 计算公式如公式(2)所示： D_CF(E)minCF(E1),CF(E2),.CF(En) (2)。 上式中， n为矿山建模数据的类型总数， D_CF(E)即为D_CF。 4.如权利要求3所述的一种多粒度矿体三维建模不确定性表达与传递方法， 其特征在 权利要求书 1/3 页 2 CN 110930504 A 2 于： 步骤S103中， 所述矿山建模过程包括： 矿体圈定、 剖面轮廓线提取和重构、 剖面多边形对 应以及矿体剖分建模过程； 对矿山建模过程误差进行可信度赋值包括对矿体圈定、 剖面轮廓线提取和重构、 剖面 多边形对应以及矿体剖分建模过程中产。

8、生的误差分别进行可信度赋值， 进而将各建模过程 对应的可信度值组成的集合作为结构建模过程不确定性集合； 其中， 根据专家经验及具体 情况分别对矿体圈定、 剖面轮廓线提取和重构、 剖面多边形对应以及矿体剖分建模过程中 产生的误差分别进行可信度赋值； 对所述结构建模过程不确定性集合中的各可信度值进行合成计算， 生成结构建模过程 可信度SMP_CF； 其中， 以单一可信度值的最小值为合成计算结果， 计算公式如公式(3)所示： SMP_CFminCF(B1),CF(B2),CF(Bn) (3) 上式中， Bi代表第i类建模过程数据； CF(Bi)代表第i类建模过程数据对应的可信度值； i 1,2,n，。

9、 n为建模数据总的类型个数。 5.如权利要求4所述的一种多粒度矿体三维建模不确定性表达与传递方法， 其特征在 于： 步骤S104中， 结构模型可信度SM_CF由建模数据可信度D_CF与结构建模过程可信度SMP_ CF相乘得到， 计算公式如公式(4)所示： SM_CFD_CF*SMP_CF (4)。 6.如权利要求1所述的一种多粒度矿体三维建模不确定性表达与传递方法， 其特征在 于： 步骤S105中， 克里格估值可信度KC_CF的计算方法， 包括以下步骤： S301： 利用步骤S101获取的多源建模数据作为建模样本， 进行实验变差函数计算和理 论模型函数拟合， 以获得的变异函数理论模型为依据设置。

10、搜索邻域的方向和变程， 并依据 所述搜索邻域的方向和变程， 确定包含在内的参与插值计算的矿体样本数目n和样本品位 值； S302： 根据所述矿体样本数目n和样本品位值， 对样本进行克里格插值计算， 构建矿山 三维属性模型， 并得到克里格估计方差； S303： 根据矿山企业对矿体三维模型预定的精度可接受范围， 得到矿体三维模型的无 偏误差区间为(1- ) 2,(1+ )2， 其中， 2为无偏差母体方差， 为矿山规划与生产部门预 先设定的误差调整系数； S304： 根据克里格估计方差和S301步骤中确定的参与插值计算的矿体样本数目n计算 可靠性概率p， p值即为克里格估值可信度KC_CF； 克里格。

11、估值可信度KC_CF的详细计算步骤， 包括： 1)设Xi为来自全部建模数据的钻孔样本子集， 为样本均值， 统计量 2X12+X22+Xn2 是自由度为n的 2分布， 记为 2 2(n)， 样本方差即克里格估计方差， 对于克里格估计方差 S2和母体方差 2， 根据 预先定义可接受精度范围的母体方差 2， 误差区间 (1- ) 2,(1+ )2； 2)根据如公式(5)所示的 2(n)分布的概率密度函数： 权利要求书 2/3 页 3 CN 110930504 A 3 计算获得克里格估计方差S2对母体误差的 2的精度表达， 如公式(6)所示： (1- ) 2S2(1+ )2 (6) 3)根据步骤2)的。

12、公式进行变换得到公式(1- )(n-1) 2(1+ )(n-1)， 进而得到如公 式(7)所示的概率公式： 由上式可得到不同体品位样品数目n下克里格估计方差S2落在其母体误差区间 2的概 率p， 将概率p表达为0,1， 即得到克里格估值可信度KC_CF。 7.如权利要求6所述的一种多粒度矿体三维建模不确定性表达与传递方法， 其特征在 于： 步骤S105中， 模拟计算可信度SC_CF的计算方法， 包括以下步骤： S401： 根据矿体当前的勘探阶段进行分类， 依据类别的不同生成与真实矿体具备不同 置信度h的多粒度矿体三维模型； 其中， 所述勘探阶段的类别包括勘查阶段、 详查阶段和生 产阶段； h取。

13、值按勘查阶段、 详查阶段和生产阶段分别赋值为0.05、 0.1和0.25； S402： 对矿体三维模型的多源矿体建模数据进行克里格估值计算， 并依据上述S301 S304的克里格估值可信度的计算过程， 计算得到克里格估值可信度KC_CF； S403： 将多粒度矿山三维模型预先划定的置信度h转换为多粒度约束可信度MGC_CF； 转 换公式如公式(8)所示： MGC_CF1-h (8) S404： 对KC_CF和MGC_CF进行叠加计算， 生成模拟计算可信度SC_CF； 叠加计算公式如公 式(9)所示： SC_CFKC_CF*MGC_CF (9)。 8.如权利要求7所述的一种多粒度矿体三维建模不确。

14、定性表达与传递方法， 其特征在 于： 步骤S106中， 多粒度矿山三维属性模型可信度AM_CF的叠加计算公式包括： 估值建模方式下的属性模型可信度计算公式如公式(10)： AM_CFSM_CF*KC_CF (10) 模拟建模方式下的属性模型可信度计算公式如公式(11)： AM_CFSM_CF*SC_CF (11)。 权利要求书 3/3 页 4 CN 110930504 A 4 一种多粒度矿体三维建模不确定性表达与传递方法 技术领域 0001 本发明涉及矿山三维地质建模技术领域， 尤其涉及一种多粒度矿体三维建模不确 定性表达与传递方法。 背景技术 0002 矿体三维模型的建模过程受矿山结构的复杂。

15、性、 地质工作者对矿体认知的不完备 性及勘探数据的稀疏性等多方面影响， 不确定性成为其固有属性并贯穿于矿体建模过程的 始终。 目前针对不确定性在三维矿体建模整个过程中的传播和定量描述仍然是尚未解决的 难题， 尤其是误差计算方法、 修正机制及其实现技术等仍处于探索阶段。 发明内容 0003 为了解决矿山三维模型不确定性表达及其传递过程研究的不足， 本发明提出一种 多粒度矿体三维建模不确定性表达与传递方法， 对建模过程中分布的多源误差进行了归纳 及分类， 对建模过程中计算型误差的主要来源克里格估值和序惯高斯模拟计算等进行了算 法解析， 利用不确定性推理网络和C-F可信度模型实现了矿体三维模型不确定。

16、性表达及传 递框架， 为矿山企业在评估、 规划及生产等方面提供了理论依据和决策支持能力。 0004 一种多粒度矿体三维建模不确定性表达与传递方法， 具体包括： 0005 S101： 获取多源矿山建模数据； 0006 S102： 对所述多源矿山建模数据产生的误差进行不确定性赋值， 形成矿山多源数 据不确定性集合， 进而对所述矿山多源数据不确定性集合内的不确定性赋值进行合成计 算， 生成建模数据可信度D_CF； 0007 S103： 对矿山建模过程误差进行可信度赋值， 形成结构建模过程不确定性集合， 进 而对所述结构建模过程不确定性集合内不确定性赋值进行合成计算， 生成结构建模过程可 信度SMP_。

17、CF； 0008 S104： 对所述建模数据可信度D_CF和所述建模过程数据可信度SMP_CF进行叠加计 算， 得到结构模型可信度SM_CF； 0009 S105： 根据矿体数据基础和矿体三维模型需求的不同， 对建立三维属性模型按模 拟建模和估值建模两种方式进行多粒度特征建模并根据建模过程分别计算不确定性， 即克 里格估值可信度KC_CF和模拟计算可信度SC_CF； 0010 S106： 将所述结构模型可信度SM_CF分别与所述克里格估值可信度KC_CF和所述模 拟计算可信度SC_CF进行叠加计算， 得到多粒度矿山三维属性模型可信度AM_CF。 0011 进一步地， 步骤S101中， 所述多源。

18、矿山建模数据包括： 钻孔位置坐标、 钻孔深度、 钻 孔样品段起点坐标、 终点坐标、 钻孔样品长度和钻孔组合样长度； 所述多源矿山建模数据根 据钻孔、 地质剖面地质资料整理获得。 0012 进一步地， 步骤S102中， 所述多源矿山建模数据产生的误差包括矿山勘测过程中 产生的误差； 对所述多源矿山建模数据产生的误差进行不确定性赋值， 形成矿山多源数据 说明书 1/7 页 5 CN 110930504 A 5 不确定性集合， 进而对所述矿山多源数据不确定性集合内的不确定性赋值进行合成计算， 生成建模数据可信度D_CF； 具体包括： 0013 S201： 对不同类型数据的误差分别进行不确定性赋值； 。

19、其中， 对某一个类型的数据 产生的误差进行不确定性赋值的方法为： 将该类型的数据多次测量所得结果的平均值与被 测量的真值之差进行比较并转换为该类型数据的可信度CF(Ei)， 计算公式如公式(1)所示： 0014 0015 上式中， Ei代表第i个类型的矿山建模数据； CF(Ei)代表该类型数据的可信度值， 取 值范围为0,1； xj代表该类型数据的第j个测量值的真值；为第j个测量值的多次测量平 均值； i1,2,n， 且不同的取值代表不同类型的矿山建模数据； 0016 S202： 将各类型数据的可信度值CF(Ei)组成的集合作为矿山多源数据不确定性集 合； 并对所述矿山多源数据不确定性集合内的。

20、所有可信度值进行合成计算， 得到建模数据 可信度D_CF； 其中， 以单一可信度值的最小值为合成计算结果， 计算公式如公式(2)所示： 0017 D_CF(E)minCF(E1),CF(E2),.CF(En) (2)。 0018 上式中， n为矿山建模数据的类型总数， D_CF(E)即为D_CF。 0019 进一步地， 步骤S103中， 所述矿山建模过程包括： 矿体圈定、 剖面轮廓线提取和重 构、 剖面多边形对应以及矿体剖分建模过程； 0020 对矿山建模过程误差进行可信度赋值包括对矿体圈定、 剖面轮廓线提取和重构、 剖面多边形对应以及矿体剖分建模过程中产生的误差分别进行可信度赋值， 进而将各。

21、建模 过程对应的可信度值组成的集合作为结构建模过程不确定性集合； 其中， 根据专家经验及 具体情况分别对矿体圈定、 剖面轮廓线提取和重构、 剖面多边形对应以及矿体剖分建模过 程中产生的误差分别进行可信度赋值； 0021 对所述结构建模过程不确定性集合中的各可信度值进行合成计算， 生成结构建模 过程可信度SMP_CF； 0022 其中， 以单一可信度值的最小值为合成计算结果， 计算公式如公式(3)所示： 0023 SMP_CFminCF(B1),CF(B2),CF(Bn) (3) 0024 上式中， Bi代表第i类建模过程数据； CF(Bi)代表第i类建模过程数据对应的可信度 值； i1,2,n。

22、， n为建模数据总的类型个数。 0025 进一步地， 步骤S104中， 结构模型可信度SM_CF由建模数据可信度D_CF与结构建模 过程可信度SMP_CF相乘得到， 计算公式如公式(4)所示： 0026 SM_CFD_CF*SMP_CF (4)。 0027 进一步地， 步骤S105中， 克里格估值可信度KC_CF的计算方法， 包括以下步骤： 0028 S301： 利用步骤S101获取的多源建模数据作为建模样本， 进行实验变差函数计算 和理论模型函数拟合， 以获得的变异函数理论模型为依据设置搜索邻域的方向和变程， 并 依据所述搜索邻域的方向和变程， 确定包含在内的参与插值计算的矿体样本数目n和样。

23、本 品位值； 0029 S302： 根据所述矿体样本数目n和样本品位值， 对样本进行克里格插值计算， 构建 矿山三维属性模型， 并得到克里格估计方差； 0030 S303： 根据矿山企业对矿体三维模型预定的精度可接受范围， 得矿体三维模型的 说明书 2/7 页 6 CN 110930504 A 6 无偏误差区间为(1- ) 2,(1+ )2， 其中， 2为无偏差的母体方差， 为矿山规划与生产部 门预先设定的误差调整系数； 0031 S304： 根据克里格估计方差和S301步骤中确定的参与插值计算的矿体样本数目n 计算可靠性概率p， p值即为克里格估值可信度KC_CF； 0032 克里格估值可信。

24、度KC_CF的详细计算步骤， 包括： 00331)设Xi为来自全部建模数据的钻孔样本子集， 为样本均值， 统计量 2X12+X22 +Xn2是自由度为n的 2分布， 记为x2 2(n)， 样本方差即克里格估计方差， 对于克里格估 计方差x2和母体方差 2， 根据 预先定义可接受精度范围的母体方差 2， 误 差区间(1- ) 2,(1+ )2； 0034 2)根据如公式(5)所示的 2(n)分布的概率密度函数： 0035 0036 计算获得克里格估计方差S2对母体误差的 2的精度表达， 如公式(6)所示： 0037 (1- ) 2S2(1+ )2 (6) 0038 3)根据步骤2)的公式进行变换。

25、得到公式(1- )(n-1) 2(1+ )(n-1)， 进而得到 如公式(7)所示的概率公式： 0039 0040 由上式可得到不同体品位样品数目n下克里格估计方差S2落在其母体误差区间 2 的概率p， 将概率p表达为0,1， 即得到克里格估值可信度KC_CF。 0041 进一步地， 步骤S105中， 模拟计算可信度SC_CF的计算方法， 包括以下步骤： 0042 S401： 根据矿体当前的勘探阶段进行分类， 依据类别的不同生成与真实矿体具备 不同置信度h的多粒度矿体三维模型； 其中， 所述勘探阶段的类别包括勘查阶段、 详查阶段 和生产阶段； h取值按勘查阶段、 详查阶段和生产阶段分别赋值为0。

26、.05、 0.1和0.25； 0043 S402： 对矿体三维模型的多源矿体建模数据进行克里格估值计算， 并依据上述 S301S304的克里格估值可信度的计算过程， 计算得到克里格估值可信度KC_CF； 0044 S403： 将多粒度矿山三维模型预先划定的置信度h转换为多粒度约束可信度MGC_ CF； 转换公式如公式(8)所示： 0045 MGC_CF1-h (8) 0046 S404： 对KC_CF和MGC_CF进行叠加计算， 生成模拟计算可信度SC_CF； 叠加计算公式 如公式(9)所示： 0047 SC_CFKC_CF*MGC_CF (9)。 0048 进一步地， 步骤S106中， 多粒。

27、度矿山三维属性模型可信度AM_CF的叠加计算公式包 括： 0049 估值建模方式下的属性模型可信度计算公式如公式(10)： 说明书 3/7 页 7 CN 110930504 A 7 0050 AM_CFSM_CF*KC_CF (10) 0051 模拟建模方式下的属性模型可信度计算公式如公式(11)： 0052 AM_CFSM_CF*SC_CF (11)。 0053 本发明提供的技术方案带来的有益效果是： 本发明所提出的技术方案根据多粒度 矿体三维模型构建过程中的误差分布情况， 利用不确定性推理网络和可信度C-F模型， 通过 对克里格估值算法和序惯高斯模拟算法的深入解析， 实现了对矿体三维模型不。

28、确定性的表 达和传递过程的描述， 对矿山企业利用三维模型进行设计、 规划和生产提供了理论依据。 附图说明 0054 下面将结合附图及实施例对本发明作进一步说明， 附图中： 0055 图1是本发明实施例中一种多粒度矿体三维建模不确定性表达与传递方法的流程 图。 具体实施方式 0056 为了对本发明的技术特征、 目的和效果有更加清楚的理解， 现对照附图详细说明 本发明的具体实施方式。 0057 本发明的实施例提供了一种多粒度矿体三维建模不确定性表达与传递方法。 0058 请参考图1， 图1是本发明实施例中一种多粒度矿体三维建模不确定性表达与传递 方法的流程图， 具体包括如下步骤： 0059 S10。

29、1： 获取多源矿山建模数据； 0060 S102： 对所述多源矿山建模数据产生的误差进行不确定性赋值， 形成矿山多源数 据不确定性集合， 进而对所述矿山多源数据不确定性集合内的不确定性赋值进行合成计 算， 生成建模数据可信度D_CF； 0061 S103： 对矿山建模过程误差进行可信度赋值， 形成结构建模过程不确定性集合， 进 而对所述结构建模过程不确定性集合内不确定性赋值进行合成计算， 生成结构建模过程可 信度SMP_CF； 0062 S104： 对所述建模数据可信度D_CF和所述建模过程数据可信度SMP_CF进行叠加计 算， 得到结构模型可信度SM_CF； 0063 S105： 根据矿体数。

30、据基础和矿体三维模型需求的不同， 对建立三维属性模型按模 拟建模和估值建模两种方式进行多粒度特征建模并根据建模过程分别计算不确定性， 即克 里格估值可信度KC_CF和模拟计算可信度SC_CF； 0064 S106： 将所述结构模型可信度SM_CF分别与所述克里格估值可信度KC_CF和所述模 拟计算可信度SC_CF进行叠加计算， 得到多粒度矿山三维属性模型可信度AM_CF。 0065 步骤S101中， 所述多源矿山建模数据包括： 钻孔位置坐标、 钻孔深度、 钻孔样品段 起点坐标、 终点坐标、 钻孔样品长度和钻孔组合样长度等类型的数据； 所述多源矿山建模数 据根据钻孔、 地质剖面等地质资料整理获得。

31、。 0066 步骤S102中， 所述多源矿山建模数据产生的误差包括矿山勘测过程中产生的误 差； 对所述多源矿山建模数据产生的误差进行不确定性赋值， 形成矿山多源数据不确定性 集合， 进而对所述矿山多源数据不确定性集合内的不确定性赋值进行合成计算， 生成建模 说明书 4/7 页 8 CN 110930504 A 8 数据可信度D_CF； 具体包括： 0067 S201： 对不同类型数据的误差分别进行不确定性赋值； 其中， 对某一个类型的数据 产生的误差进行不确定性赋值的方法为： 将该类型的数据多次测量所得结果的平均值与被 测量的真值之差进行比较并转换为该类型数据的可信度CF(Ei)， 计算公式如。

32、公式(1)所示： 0068 0069 上式中， Ei代表第i个类型的矿山建模数据； CF(Ei)代表该类型数据的可信度值， 取 值范围为0,1； xj代表该类型数据的第j个测量值的真值；为第j个测量值的多次测量平 均值； i1,2,n， 且不同的取值代表不同类型的矿山建模数据； 0070 S202： 将各类型数据的可信度值CF(Ei)组成的集合作为矿山多源数据不确定性集 合； 并对所述矿山多源数据不确定性集合内的所有可信度值进行合成计算， 得到建模数据 可信度D_CF； 其中， 以单一可信度值的最小值为合成计算结果， 计算公式如公式(2)所示： 0071 D_CF(E)minCF(E1),CF。

33、(E2),.CF(En) (2)。 0072 上式中， n为矿山建模数据的类型总数， D_CF(E)即为D_CF。 0073 步骤S103中， 所述矿山建模过程包括： 矿体圈定、 剖面轮廓线提取和重构、 剖面多 边形对应、 矿体剖分等建模过程； 0074 对矿山建模过程误差进行可信度赋值包括对矿体圈定、 剖面轮廓线提取和重构、 剖面多边形对应以及矿体剖分等建模过程中产生的误差分别进行可信度赋值， 进而将各类 建模过程对应的可信度值组成的集合作为结构建模过程不确定性集合； 其中， 根据专家经 验及具体情况分别对矿体圈定、 剖面轮廓线提取和重构、 剖面多边形对应以及矿体剖分等 建模过程中产生的误差。

34、分别进行可信度赋值， 取值范围为0,1； 0075 对所述结构建模过程不确定性集合中的各可信度值进行合成计算， 生成结构建模 过程可信度SMP_CF； 0076 其中， 以单一可信度值的最小值为合成计算结果， 计算公式如公式(3)所示： 0077 SMP_CFminBF(B1),BF(B2),BF(Bn) (3) 0078 上式中， Bi代表第i类建模过程数据； BF(Bi)代表第i类建模过程数据对应的可信度 值； i1,2,n， n为建模数据总的类型个数。 0079 在本发明实施例中， 以构建矿体表面模型过程中的剖面轮廓线建立为例， 当矿体 形态简单、 内部无夹石分支、 轮廓线点列摆列规律且。

35、点数差别小时， 解译结果具有较高的可 信度， 可赋值为0.80.9， 对于断层、 夹石和分枝等情况的存在的复杂的矿体， 该类型矿体 表面模型拟合度差， 可赋予较低可信度， 如0.20.3。 0080 步骤S104中， 结构模型可信度SM_CF由建模数据可信度D_CF与结构建模过程可信 度SMP_CF相乘得到， 计算公式如公式(4)所示： 0081 SM_CFD_CF*SMP_CF (4)。 0082 步骤S105中， 克里格估值可信度KC_CF的计算方法， 包括以下步骤： 0083 S301： 利用步骤S101获取的多源建模数据作为建模样本， 进行实验变差函数计算 和理论模型函数拟合， 以获得。

36、的变异函数理论模型为依据设置搜索邻域的方向和变程， 并 依据所述搜索邻域的方向和变程， 确定包含在内的参与插值计算的矿体样本数目n和样本 品位值； 说明书 5/7 页 9 CN 110930504 A 9 0084 S302： 根据所述矿体样本数目n和样本品位值， 对样本进行克里格插值计算， 构建 矿山三维属性模型， 并得到克里格估计方差； 0085 S303： 根据矿山企业对矿体三维模型预定的精度可接受范围， 得矿体三维模型的 无偏误差区间为(1- ) 2,(1+ )2， 其中， 2为无偏差的母体方差， 为矿山规划与生产部 门预先设定的误差调整系数； 0086 S304： 根据克里格估计方差。

37、和S301步骤中确定的参与插值计算的矿体样本数目n 计算可靠性概率p， p值即为克里格估值可信度KC_CF； 0087 克里格估值可信度KC_CF的详细计算步骤， 包括： 00881)设Xi为来自全部建模数据的钻孔样本子集， 为样本均值， 统计量 2X12+X22 +Xn2是自由度为n的 2分布， 记为 2 2(n)， 样本方差即克里格估计方差， 对于克里格估 计方差S2和母体方差 2， 根据 预先定义可接受精度范围的母体方差 2， 误 差区间(1- ) 2,(1+ )2； 0089 2)根据如公式(5)所示的 2(n)分布的概率密度函数： 0090 0091 计算获得克里格估计方差S2对母体。

38、误差的 2的精度表达， 如公式(6)所示： 0092 (1- ) 2S2(1+ )2 (6) 0093 3)根据步骤2)的公式进行变换得到公式(1- )(n-1) 2(1+ )(n-1)， 进而得到 如公式(7)所示的概率公式： 0094 0095 由上式可得到不同体品位样品数目n下克里格估计方差S2落在其母体误差区间 2 的概率p， 将概率p表达为0,1， 即得到克里格估值可信度KC_CF。 0096 模拟计算可信度SC_CF的计算方法， 包括以下步骤： 0097 S401： 根据矿体当前的勘探阶段进行分类， 依据类别的不同生成与真实矿体具备 不同置信度h的多粒度矿体三维模型； 其中， 所述。

39、勘探阶段的类别包括勘查阶段、 详查阶段 和生产阶段； h取值按勘查阶段、 详查阶段和生产阶段分别赋值为0.05、 0.1和0.25； 0098 S402： 对矿体三维模型的多源矿体建模数据进行克里格估值计算， 并依据上述 S301S304的克里格估值可信度的计算过程， 计算得到克里格估值可信度KC_CF； 0099 S403： 将多粒度矿山三维模型预先划定的置信度h转换为多粒度约束可信度MGC_ CF； 转换公式如公式(8)所示： 0100 MGC_CF1-h (8) 0101 S404： 对KC_CF和MGC_CF进行叠加计算， 生成模拟计算可信度SC_CF； 叠加计算公式 如公式(9)所示。

40、： 0102 SC_CFKC_CF*MGC_CF (9)。 说明书 6/7 页 10 CN 110930504 A 10 0103 步骤S106中， 多粒度矿山三维属性模型可信度AM_CF的叠加计算公式包括： 0104 估值建模方式下的属性模型可信度计算公式如公式(10)： 0105 AM_CFSM_CF*KC_CF (10) 0106 模拟建模方式下的属性模型可信度计算公式如公式(11)： 0107 AM_CFSM_CF*SC_CF (11)。 0108 本发明的有益效果是： 本发明所提出的技术方案根据多粒度矿体三维模型构建过 程中的误差分布情况， 利用不确定性推理网络和可信度C-F模型， 通过对克里格估值算法和 序惯高斯模拟算法的深入解析， 实现了对矿体三维模型不确定性的表达和传递过程的描 述， 对矿山企业利用三维模型进行设计、 规划和生产提供了理论依据。 0109 以上所述仅为本发明的较佳实施例， 并不用以限制本发明， 凡在本发明的精神和 原则之内， 所作的任何修改、 等同替换、 改进等， 均应包含在本发明的保护范围之内。 说明书 7/7 页 11 CN 110930504 A 11 图1 说明书附图 1/1 页 12 CN 110930504 A 12 。