子域估值法在甲玛铜多金属矿资源量估算中的应用

在传统品位插值过程中，一个矿化域一般只使用单个搜索椭球体。然而对于某些矿体形态复杂或者局部矿化变化强烈的矿化域时，单一的搜索椭球体并不能满足对该矿化域进行合理估值，从而影响了资源储量估算可信度和资源分类，以及后续的采矿设计和排产等技术工作。以甲玛铜多金属矿为研究对象，利用Surpac三维软件子域估值法对甲玛铜多金属矿进行估值和评价，寻求解决这种局部形态或矿化变化强烈的矿体资源储量估算方法，子域估值法是一种很好的方法，更能确保资源储量可信度。

1 地质数据库

地质数据库是数据库技术在矿山地质工作中的实际应用，是矿山矿产资源储量估算的基础，是采矿设计的依据

，是矿山生产管理的有效手段。地质数据的完整性和可靠性直接影响矿山生产经营和决策

。因此，建立地质数据库是实现地质资源信息数字化管理的最基础工作

。

通过收集甲玛多金属矿历年以来地质钻孔、探槽、坑内钻、巷道刻线等数据，获得此次地质资源数据库地质属性信息，并建立相应的数据表，以存储相应的地质资源信息

。

2 矿化域三维模型

2.1 矿化域范围圈定

在充分研究甲玛铜多金属矿床地质特征、成矿理论和控矿因素的基础上，根据矿床地质特征圈定矿化域。在同一矿化域内，以Cu≥0.3%、Mo≥0.03%、Pb≥0.3%、Zn≥0.5%作为边界品位进行混圈。若一样品中任意一种金属元素品位大于或等于上述边界品位时，该样品则圈入矿化域内。

在图书馆创客空间运动如火如荼开展的热潮中，学者们对创客空间建设的必要性进行了探讨，主要分为两种观点，一种观点认为建设创客空间非常必要，但也有学者从不同角度提出冷思考。

2.2 建立矿化域三维实体模型

以勘探线剖面进行矿化域地质解译，按地质规范进行外推。完成整个矿化域解译后，再以三角网化工具手动连结三角网，最终形成验证通过可用于块体模型估算的矿化域实体模型。

3 数据处理

3.1 样长统计分析

提取矿内样品经统计分析得到样本数107777个，样长均值1.67m、中值为2.00m，结合甲玛矿山地质勘探及生产取样实际，得到工程组合样长为2m。

式中：H为测压管水头，m；c为水锤波速度，m/s；x，t分别为水锤波沿直线传播的距离和时间；g为当地重力加速度，m/s2；f(v)是速度的非线性函数。

3.2 特异值处理

根据对勘探工程组合数据分析得到Cu品位最小值0.0001%、最大值32.694%、Cu均值0.80%，标准差0.99，Sichel-t 0.94，Cu元素的品位累积频率曲线总体上呈一条直线,说明Cu元素分布基本上符合正态分布,同时，Cu元素品位累积频率曲线(图1)上部又都存在少量离散点，且偏离直线方向,这说明样品中存在特高品位

。通过数据转换得到Cu(log10)累计频率直方图（图2），偏度-0.75，峰度4.40，基本符合正态分布条件，使用数理统计的方法进行特高品位下限值处理方法，95%的置信区间

(CI)的上限值=Cu均值+1.96标准差=0.8+1.96×0.99，得到Cu特高品位下限值2.74%。

由计算结果可知，渗透系数的变化与排水量呈正比关系，其中有限元计算结果的比例系数为3.66，方法一计算结果的比例系数为3.13。而渗透系数的变化对剩余水头高度则没有影响。

通过线串运算的方法截取Cu品位的特异值，统计分析Cu品位最小值0.0001%、最大值2.74%、均值0.74%、Sichel-t值为0.89%，特高品位处理合理，变化系数95%，Cu元素在矿化域分布均匀。

4 子域划分

（1）陡Cu1子域变异函数及椭球体

子域变异函数研究：

应用实体内外区分点功能把组合样线文件按矿化域产状分割成三个独立的矿化域，以每个子矿化域为独立矿化域进行变异函数分析研究，通过变异函数球状模型拟合，得到各子域变异函数的块金值、基台值、变程及椭球体参数。

2017年6月21日的国务院会议中，李克强总理提出“分享经济”对于新常态下的中国经济发展具有重要意义。在当今的中国“分享经济”充分利用互联网技术，不断的创造出新的产品和服务模式，为供给侧改革、化解产能、提供工作岗位做出了重要贡献。与“分享经济”相匹配，人才市场也出现了许多“分享经济”的特征，出现了全新的雇佣模式。李克强总理提到，他在南方省市考察时发现当地出现了“零工市场”，工人们没有固定的雇主，工作按小时收费，一天可以服务多个雇主。

在以往甲玛铜多金属矿资源储量估算时，一个矿化域使用一个椭球体参数，而矿体具有明显的陡→缓→陡的变化特征(图3)。因此，使用单一椭球体(图4)不能代表整个矿化域产状，不能合理对整个矿化域插值评价，从而影响资源储量估算精度。为保证资源储量估算可信度，在研究以往地质资料的基础上

，结合矿化域产状将其分为三个子域。以Cu元素为例，分为陡Cu1、缓Cu、陡Cu2三个子域（图5）。

当前，解决人民日益增长的美好生活需求和不平衡不充分的发展之间的矛盾，改善企业社区居民生活条件已是势在必行。改造企业社区，政府发挥主导作用的同时，积极鼓励多元力量参与企业创新型社区的转型和建设。利用、改造、发展、引导，多管齐下，实现共享共治，打造全新的泰煤家园！

（2）缓Cu子域变异函数及椭球体

分别对三个独立子矿化域插值,组合样与估值结果验证，得到整体插值Cu(绿色线)及各子域插值Cu(黄色线)Swath曲线(图9)。

1.1 临床资料 选取2016年1月至2018年1月在我院呼吸内科门诊就诊及住院的长期吸烟哮喘患者。入选条件：符合中华医学会制定的支气管哮喘防治指南哮喘诊断标准；年龄18-60岁之间；吸烟史≥10包年。两组患者治疗前一般资料比较差异无统计学意义（P＞0.05）。

（3）陡Cu2子域变异函数及椭球体

（4）模型验证

两组干预前腕背伸、肘屈曲、肩外展、膝屈伸活动度组间无显著差异(P>0.05);干预后观察组各项指标均明显优于对照组(P<0.05)。见表1。

以上三个子域变异函数通过交叉验证得到满足理想结果的三个条件：一是，真值与估计值误差的均值趋近于“0”；二是，误差方差与克里格估计方差比值趋近于“1”

；三是，两个标准差之间的置信区间大于等于95%。

陡Cu2子域样本数2448个，Cu品位最小值0.0001%、最大值2.74%、均值0.46%、Sichel-t0.65%，变化系数124%，偏度2.21，峰度8.30，通过log10对数转换后统计分析偏度0.70、峰度3.04,基本符合正太分布。通过对缓Cu2变异函数球状模型拟合（图8）和交叉验证(表3)，得到陡Cu2椭球体参数，方位75°，倾伏角0°，倾角45°，主/次轴比1.11，主/最小轴比5.77。

缓Cu子域样本数56972个，Cu品位最小值0.0001%、最大值2.74%、均值0.78%、Sichel-t0.95%，变化系数88.8%，偏度1.22，峰度3.85，通过log10对数转换后统计分析偏度1、峰度5，基本符合正太分布。通过对缓Cu变异函数球状模型拟合（图7）和交叉验证(表2)，得到缓Cu子域椭球体参数，方位120°，倾伏角0°，倾角45°，主/次轴比1.58，主/最小轴比5.07。

陡Cu1子域样本数33782个，Cu品位最小值0.0001%、最大值2.74%、均值0.70%、Sichel-t0.73%，变化系数104%，偏度1.50，峰度4.4，通过log10数据转换后统计分析偏度0.54、峰度3.79,符合正太分布。通过对陡Cu1变异函数球状模型拟合（图6）和交叉验证(表1)，得到了Cu椭球体参数，方位105°，倾伏角0°，倾角45°，主/次轴比1.28，主/最小轴比1.74。

从Swath曲线得出，使用子域估值法估值，估值结果更接近原始组合样，有利于资源储量估算可信度和资源分类，以及后续的采矿设计和排产等技术工作。

（5）资源储量对比。

由表5得出，用子域估值法估值，以Cu 0.3%为边界品位，从矿石体积、矿石量、金属量都相应减少，Cu资源储量体积减少0.023%，吨位减少0.024%，金属量减少0.005%。

5 结论

（1）在甲玛铜多金属矿根据矿体产状分割成若干个子域，分别进行子域变异函数研究及插值得到较好的应用。

（2）对于解决这种局部形态和矿化变化剧烈的矿体，子域估值法是一种很好的方法，更能确保资源储量可信度。

（3）子域估值法有效剔除了废石量，提高了资源储量估算的精度，能够更好的指导后续的采矿生产活动，增加经济效益。

（4）子域界限的划分难度和工作量比较大，对于特别复杂矿体的估算耗时且容易出错不适用。

[1]李雪峰,宋志伟,吕海栋,等.一种矿山地质资源储量管理方法:,CN106529755A[P].2017.

[2]张国亮,张海峰,周勇.基于Surpac的某珍珠岩矿三维建模与露天境界圈定[J].采矿技术,2012,012(001):11-13.

[3]车海龙.Micromine软件在某金矿床资源储量估算中的应用[J].吉林地质,2020,39(3):7.

[4]郑文宝,黎枫佶,唐菊兴,等.基于Micromine软件下地质统计学在甲玛矽卡岩型铜多金属矿储量计算中的应用[J].地质与勘探,2011,47(4):726-736.

[5]尚琼.安徽省沙溪铜金矿床凤台山矿段三维建模及成矿元素空间分析[D].合肥工业大学,2014.

[6]周娣.滑绳溜车保护在多绳摩擦提机中的应用[J].中国金属通报,2020(4):2.

[7]唐攀,唐菊兴,唐晓倩,等.传统方法和地质统计学在矿产资源/储量分类中的对比分析[J].金属矿山,2013.