摘要
地质勘探对于资源开发具有不可或缺的重要性，准确的资源储量报告是体现勘查成果的核心环节，而基于三维地质建模和地质统计学的资源量估算方法受到了地勘行业一线技术人员的日益重视。依据2020 年《固体矿产资源量估算规程》标准，本文详细阐述了使用 3DMine Plus矿业工程软件在现行规范要求下进行资源量估算的工作原理和流程，包括矿体实体模型构建、数据处理、变异函数分析、块体估值等关键环节。在新疆东戈壁钼矿床补充勘探报告中，利用矿区最新勘查数据，采用了与变异函数紧密结合的距离幂次反比法对资源量进行估算，并与传统的平行断面法估算结果进行了比较分析，二者矿石量、金属量、平均品位相对误差分别为5. 81%、2. 67%、3. 67%。验证了3DMine Plus软件在资源量估算模块中的简便性和准确性，也证明了其在地质勘查和矿山生产管理等实际地质工作中的可靠性和实用性。
关键词
Abstract
In the geological prospecting industry， it pays a more and more important attention for front-line technicians to use the resource estimation method based on 3D geological modeling and geostatistics method for recently years， which attributing to geological exploration plays an indispensable role in resource development， and it is the fundamental aspect to reflect the exploration results for an accurate resource reserve report. This paper provides a comprehensive explanation of the operational principles and processes involved in utilizing 3D Mine Plus mining engineering software for resource estimation， adhering to the current regulatory requirements， which includes crucial steps such as constructing a solid model of the ore body， processing data， conducting variance function analysis， and performing block valuation， according to the standards set forth Solid Mineral Resources Estimation Regulations in the 2020. This study uses the latest exploration data in the supplementary exploration report of Donggebi Mo deposite in Xinjiang， in comparison with traditional parallel section method， employing the Inversed Distance Weighting method which is closely integrated with variance function analysis， the results reveal relative errors of 5. 81% for ore content，2. 67% for metal content， and 3. 67% for average grade respectively. It can be concluded that the simplicity and accuracy of resource estimation module by using 3D Mine Plus software have been verified while its reliability and practicality in geological work such as geological exploration and mine production management have also been demonstrated.
0 引言
新疆东戈壁钼矿床是笔者所在单位于 2006 年投入勘查、2010 年探明的一处特大型钼矿床，提交钼金属量 50.8 万 t，平均品位 0.115%（河南省地矿局第二地质勘查院，2010①）。该矿床品位高、埋藏浅，适宜大规模露采。随着近年来国家日益重视矿产资源安全，叠加钼矿产品价格上行因素，该矿山开发计划被提上日程。2023年，笔者所在单位对东戈壁钼矿床开展了补充勘探。
对地勘人员而言，资源储量估算是矿产勘查工作的最重要成果，同时精准快速地估算矿体品位和储量对矿山生产管理和规划具有重要意义（张炎等，2011）。中国固体矿产资源资源量估算常用的方法包括几何法、地质统计学法、距离幂次反比法、 SD法等。随着计算机技术的发展，矿业工程软件的操作性和实用性得到长足进步，近年来新颁布的地质勘查规范中倡导使用软件估算资源储量的频次明显增加。本文以新疆东戈壁特大型钼矿床为例，详细阐述在现行勘查规范要求下（中华人民共和国自然资源部，2020a），以满足资源储量报告提交条件为前提，3DMine Plus软件在该矿床三维建模和资源量估算中的工作原理及流程，以期为行业同仁在生产实践中提供部分借鉴。
1 区域地质背景
矿区构造位置位于塔里木板块北缘，康古尔断裂与雅满苏断裂之间（图1），以雅满苏大断裂为界，本区横跨黄山—秋格明塔什复理石岩带和阿奇山 —雅满苏岛弧带，出露地层主要以晚古生代火山-沉积地层为主，构造主要表现为一系列走向北东或近东西的断裂带、复式褶皱等，以褶皱为主，断裂次之。华力西构造回旋期，中酸性—酸性岩浆活动频繁，以花岗岩类为主，闪长岩类次之，多呈长条带状和轴状岩株产出，破坏了褶皱、断裂构造的完整性，侵入的最新地层为二叠系。矿区位于康古尔塔格铜、镍、金、钼、稀有金属成矿带，矿化分布广泛，种类多样，大小矿床（点）星罗棋布（张永等，2017）。
2 矿床地质特征
东戈壁钼矿床产在隐伏斑状花岗岩东西两侧外接触带上，赋矿岩石为石炭系干墩组浅变质碎屑岩，少部分为变质辉绿岩、变质安山岩，矿化对围岩没有选择性，矿体与围岩界线一般不清楚，依靠分析结果而定。依据产出位置可分东、西2个矿段，以东矿段为主矿体（图2）。
勘查工作控制主矿体沿勘查线走向最大长度 1448 m，垂直勘查线走向最大长度 1534 m，平面投影面积 1.58 km2。主矿体最大厚度 417.84 m，平均厚度 149.31 m，平均品位 0.113%。矿体总体产状平缓，由中心部位向周边缓慢倾斜，南部及西部倾角 0°~5°，东北部变陡倾至 30°。根据矿化品位圈出的矿体呈近似层状—透镜状，矿体中心部位不含或少含低品位矿及夹石，向周边低品位矿及夹石增多，矿体分枝变多、变薄。低品位矿呈透镜状、似层状分布于工业矿中，或与工业矿、夹石互层、相间产出，向矿体边部厚度加大，矿体整体形态偏复杂（图3）。
3 资源量估算的方法选择及流程
东戈壁钼矿床主矿体中心区域以70 m×70 m工程间距控制，工程数量多（126个钻孔）、样品数量多（17971 件化学样品），矿体规模大，形态复杂，采用传统资源量估算方法的弊端是费时费力、估算过程易受人为干扰。为提高效率、准确刻画矿体三维真实形态，减小资源量估算偏差，本文资源储量估算方法选择距离幂次反比法，该方法操作流程较为简单，受原始样品数据空间分布和空间构型特征影响小（谢徽等，2023a），相比克里格方法，更适合地勘技术人员使用。同时与平行断面法进行估算结果对比，确保本次估算的准确性和可靠性。
图1东天山大地构造略图
Ⅰ—准噶尔板块；Ⅰ1—小热泉子—大南湖岛弧带；Ⅱ—塔里木板块；Ⅱ1—秋格明塔什—黄山复理石岩带；Ⅱ2—阿齐山—雅满苏岛弧带；Ⅱ3— 卡瓦布拉克地块；Ⅱ4—帕尔岗地块；Ⅱ5—兴地地块；Ⅱ6—北山华里西期裂谷；Ⅱ7—巴伦台地块；Ⅱ8—南天山碰撞带；Ⅱ9—库鲁克塔格地块； F1—康古尔塔格断裂；F2—雅满苏断裂；F3—沙泉子断裂；F4—拱拜子断裂；F5—帕尔岗断裂；F6—卡瓦布拉克断裂；F7—辛格尔断裂；F8—兴地断裂；F9—星星峡断裂；F10—阿拉塔格断裂；F11—车尔臣断裂
图2东戈壁钼矿区地质简图
采用距离幂次反比法进行资源量估算流程见图4，其中最为关键的3个步骤为：（1）根据工程指标体系或矿块指标体系圈定估算域；（2）依据矿体产状、变异函数模型和勘查工程间距来确定搜索椭球体的产状、轴比及半径；（3）采用距离幂次反比法对块体进行估值。
4 数据准备
本文所建立钻孔数据库利用的勘查工程数据包括126个钻孔的孔口坐标、地质编录和测斜数据、 17971件化学样品的测试数据、岩性数据和1∶10000地形测量数据。数据在 Excel表格分类整理完毕后导入3DMine软件形成钻孔数据库，数据库按照定位表、测斜表、化验表和岩性表分开导入（赵少攀等， 2018），各表数据主要结构见表1，钻孔分布见图5。
图3戈壁钼矿区H07勘探线剖面图
图4距离幂次反比法资源量估算流程图
5 估算流程
5.1 地质解译建立矿化域或者矿体模型
根据最新勘查规范（中华人民共和国自然资源部，2020a），使用距离幂次反比法估算资源量时，可根据矿体和矿化域范围进行估值确定。前者与几何法一致，采用工程指标体系圈定矿体边界，建立矿体模型；后者采用矿块指标体系，即首先根据地质矿化规律和经济评价合理确定一个品位（一般介于地质矿化品位和工程指标体系中的边界品位之间），以此品位边界圈定矿化域，建立地质模型，继而对矿化域内单元块进行品位估值，划分出矿石和废石，统计资源量。
规范中对如何选用 2 种指标体系的判别条件是：当品位分布均匀、矿体与围岩是渐变关系时，采用矿块指标体系偏优；相反条件下，采用工程指标体系偏优。而对于矿化域的品位边界，可从样品分布直方图、累计分布曲线、概率图中寻找矿化拐点确定，或通过类比法确定，或采用边界品位的 1/2或根据尾矿品位值确定。事实上，矿块指标体系在20 世纪 80 年代由国外引入，写入新规范的时间不长，针对上述矿块指标体系使用中的 2 个问题，缺乏具体的明确规定，因此对国内勘查生产一线人员造成一定困扰，目前该指标体系使用普及度并不高。
东戈壁钼矿矿床类型及矿化特征并不同于典型的斑岩型矿床，表现为如下方面：矿体仅赋存于斑岩体外围围岩中；矿石组构主要呈脉状、细脉状，钼矿化与石英脉、钾长石-石英脉高度相关；矿体与围岩存在一定突变性。根据以上特征，分析东戈壁钼矿矿体边界属性应介于硬边界和软边界之间（高帮飞，2018），结合新规范要求，本文估算工作仍采用工程指标体系圈定矿体。矿体圈定从单工程开始，按照单工程→剖面→三维矿体顺序，依次圈连，步骤如下：
表1钻孔数据表结构
图5东戈壁钼矿钻孔分布图
（1）根据数据库结合勘查线生成矿区剖面图，按照资源量估算工业指标和矿体圈定原则，根据矿体赋存部位、产状及变化规律、矿石类型、品级等因素，将相邻工程对应的见矿部位进行连接（图6a）；
（2）根据勘查线剖面图上矿体圈定线框，把线框在空间上展开（图6b）并将相邻剖面对应的线框进行实体连接（图6c）；
（3）根据上述确定的实体连接及外推原则，组建矿床三维矿体模型（图6d）。
通过以上方法的连接，分别建立矿体模型、夹石模型、岩性模型。
5.2 样品数据统计分析及特异值处理
对样品统计分析及特异值处理时，应以矿化域或矿体为单位进行，因为将矿体品位值作为区域化变量对待，对其研究应限制在一个品位值具有相关性的空间里。实际工作中通常先在钻孔数据库化验表中建立一个属性，根据矿体圈定结果，人工标记出矿体或者矿化域内的样品（卢大超和付友山， 2010）。例如 ZK0020 钻孔（图7），H27-32 样品段和 H42-47 样品段为矿体，其余为围岩，建立一个 zone 属性，将矿体标记为1，围岩或夹石不标记。然后在 3DMine 软件平台中以 zone=1 为约束条件提取样品点，并保存所需要的样品点文件。
图6矿床三维实体模型构建流程
a—H04勘查线剖面上圈连矿体；b—横断面联合剖面示意图；c—相邻剖面矿体连接；d—矿床三维模型
Mo作为一种低丰度元素，其原始品位值服从偏正态分布（图8a），因此统计时需对 Mo 原始品位值进行对数处理：采用先放大一定倍数，然后进行对数变换，再加入一个第三参数进行调偏的方法，使数据更好地服从标准正态分布。本研究采用 Mo原始品位放大1000倍，加入+3值调偏后的对数值来进行统计分析。处理公式如下：
(1)
式（1）中：A1为处理后的样品 Mo品位值（%），A 为Mo原始品位值（%）。
对数处理后样品统计特征值见表2、图8b。
特高品位处理是资源量估算中的重要一环，其目的是消除特高品位对资源量估算的夸大影响，尤其是空间上分布离散的特高品位值影响更为突出 （谢徽等，2022）。对于如何确定和处理特高品位，目前尚无通用标准，根据勘查规范，几何法、地质统计学法、SD法都有各自的处理方法。行业内实践中常用的确定特高品位方法有：平均品位的 6~8 倍法 （品位变化系数法）、累计概率分布曲线法、分位数法、标准差法（3σ 准则法）等；处理方法一般有代替法和截取法。
本研究针对以上常用 4 种方法进行了对比，各方法确定的特高品位原始 Mo品位下限值和需处理样品个数见表3。
特高品位是贵金属和部分有色金属矿产的重要特征之一，应处理但影响范围又不宜过大，事实上，不管用何种数学方法确定特高品位下限后，地质技术人员都应从地质的角度综合判定是否需要处理。根据东戈壁钼矿矿体及矿化特征，综合 4 种方法数据对比，本文选择品位变化系数法：将单样品位值高于（工业）矿体加权平均品位（0.150%）7 倍（即品位≥1. 050%）的样品视为特高品位，处理的方法是用单工程样品段矿体平均品位代替该样品品位参与计算。
表2不同类型样品基本统计分析
注：Mo为对数处理后的值。
图7提取矿体（矿化域）内样品数据流程图
对特高品位处理后的合理性检验也是必须进行的一步工作（杨丰铭等，2022），经验证，处理后样品算数平均值略小于西舍尔估值，相比特高品位处理前（表2），二者数值更接近，同时偏度和变异系数进一步减小；处理后东矿段主矿体金属量损失量 18796 t，占处理前总金属量的 5.31%，说明特高品位处理结果合理。
特高品位处理完毕后，须进行等样长组合，作为下步估值的依据。本研究选择2 m作为组合样长度的原因是：统计所有钻孔样品（17971 件）样长数据后，样长 2 m 的样品数量占比最大，达 60%，组合前后样段长度变量统计特征能最大程度保持一致。
图8主矿体样品品位（%）统计直方图
a—原始样品Mo统计分布直方图；b—对数处理后样品Mo统计分布直方图
表3特高品位处理对比
5.3 估值参数的确定
5.3.1 块模型参数的确定
块体模型是矿床品位及资源储量估算的基础。建立块体模型的基本思想是将矿床在三维空间内按照一定的尺寸划分为一定的单元块，然后对矿体 （矿化域）范围内的单元块的品位根据已知的样品进行估值，并在此基础上进行资源储量估算。单元块尺寸的选取应当依据以下参数：矿体规模、勘探工程间距、品位变化程度、采矿设计的最小开采段高、克里格效率等。块大小必须能反映品位分布，块不宜太大，否则不能准确显示品位变化；块也不能太小，否则数据文件太大，并且会增大估值方差，品位估值变得不可靠（李赋屏等，2005；刘宇英和李勇，2019）。对于单元块XY平面尺寸，根据经验参考值，一般选择本矿区最密勘查工程间距 1/4~1/8（杜菊民等，2022），Z 尺寸宜是中段或者台阶高度的约数，且能与组合样长度匹配，包括 1~N 个完整的组合样。
东戈壁钼矿主矿体沿走向、倾向长度基本相等，探求探明资源量的工程间距为70 m×70 m，设计台阶高度 12 m，组合样样长 2 m。本次块模型单元块标准 XY尺寸选择探明工程间距的 1/4、1/5、1/7，Z 尺寸统一选择 2 m，即 17.5 m×17.5 m×2 m，15 m× 15 m×2 m，10 m×10 m×2 m 共 3 种尺寸进行对比，采用统一估值参数进行普通克里格法估值，对比 Mo 品位、克里格条件偏斜、克里格效率等指标（表4），最终确定块体尺寸取工程间距的 1/4（17.50 m），次级块尺寸取标准块尺寸的1/2。块体模型完成后，添加矿岩类型、比重、元素品位、资源量级别、估值利用钻孔数等属性（滕正双和于红萍，2015）。
表4矿体块模尺寸估值参数对比
利用矿体范围的约束，求出矿体内块体的体积，与矿体的实体体积进行对比分析（表5），体积相对偏差小于 1%，表明块体模型尺寸参数设置的合理性。
5.3.2 变异函数模型参数的确定
变异函数模型是地质统计学的重要工具，其反映样品品位的空间相关性，求变异函数主要是为了找到在某个方向上矿体以及品位具有最好的连续性（赵向东等，2021）。变异函数的模型有很多，例如球状模型、指数模型、高斯模型、孔穴效应模型等，在实际工作中最常使用的是前两种模型。变异函数模型的构建主要分为 3 步：（1）计算变异函数； （2）结构分析；（3）交叉验证。
表5矿体体积与块体体积对比
计算变异函数，就是为了获取矿体三维变异函数的8个参数，包括块金值、基台值、变程；椭球体走向、倾向、倾角；主轴与次轴之比、主轴与短轴之比。变异函数计算过程是由软件系统自行完成的，计算结果的好坏，则受地质技术人员实操经验以及参数调整影响，仍具有一定主观臆断性。具体操作过程中，块金值通过在小滞后距（步长）尺度下，计算沿孔变异函数来确定；基台值通过选择全向变异函数进行拟合来确定，并确定最优变程（谢徽等， 2023b）；各轴向变异函数则通过一系列计算得到的沿水平方向样品值的实验变异函数，结合矿体产状特征，经多次试算，选择出合适的那一组进行理论变异函数的拟合得出。
东戈壁钼矿区的勘查类型属Ⅰ类型偏复杂型，东矿段主矿体探明的勘查间距为 70 m×70 m，矿体走向北西部、南部近水平，东北部倾角约 30°，整体平均为 25°，在模型上矿体又具有侧倾角，整体为 20°。首先在小滞后距（步长）尺度下，计算沿孔变异函数来确定块金常数，经过变异函数的拟合，块金常数取 0.90；全向变异函数显示基台值 0.40，最优变程为140 m（图9a）；然后确定各轴向变异函数，设置12个扇区，通过调整面参数、角度误差限、搜索距离以及步长等参数进行试算，多次对比和筛选后，确定主轴在 145°方向，倾角-35°（图9b）；次轴在 235°方向，倾角 0°（图9c）；短轴在 325°方向，倾角-55°（图9d）时，函数曲线最符合正态分布，也代表着矿体的最大连续性方向，并采用球状模型对获得的各个方向的实验变异函数进行拟合（李晓晖等， 2023）。
图9矿体变异函数分析
a—全向变异函数分析；b—主轴变异函数分析；c—次轴变异函数分析；d—短轴变异函数分析
本文建立的变异函数模型具有几何异向性，即具有相同的基台值而变程不同，对不同方向的变异函数进行套合，不同方向上的变异性之差可以用变程之比来表示：K=a1/a2。K 代表各项异性比，表示在 a1 方向上距离为 h 的两点间的平均变异程度与在a2方向上距离为Kh的两点间的平均变异程度相同。求得主轴与次轴之比为1.20，主轴与短轴之比为1.95。
经交叉验证，品位估计值与实际值的误差均值为-0. 0021，接近于“0”；误差方差与克里格方差比值趋近于“1”（图10），符合交叉验证的2个理论判断依据，可视为合格。
5.3.3 搜索椭球体参数的确定
通过上述对东矿段主矿体 Mo元素开展的变异函数研究，确定出椭球体三轴的方位、轴比，一般情况下这些参数与变异函数确定的各轴产状与变程比相同，确定完成后，应将新绘制的搜索椭球体投到矿体中从各方向进行观察，判断是否合适。主轴搜索半径应根据工程间距、样品数和估算域规模确定（中华人民共和国自然资源部，2020b），东戈壁钼矿勘查类型为Ⅰ类偏复杂，主矿体中心区域以 70 m×70 m 工程间距控制，样品数量足够，所以第一次估值主轴搜索半径为70 m，第二次和第三次估值主轴搜索半径依次按倍数递增，分别为 140 m、280 m。由于样品沿钻孔方向密集，为防止丛聚效应，保证估值结果准确，一般采用八分圆搜索方式，样品个数也应多次实验对比，不宜过多或过少（谢徽等， 2023b）。本文估值搜索椭球体产状、分区数、样品数设置情况见图11所示。
图10交叉验证残差图
图11距离幂次反比法估值搜索参数设置图
5.3.4 距离幂次的确定
距离反比加权插值法对于品位变化系数高，矿体零星而分散的矿种，需要选择较大的幂。距离幂次反比法的幂次不同，则有不同的适用范围和估值效果，当幂次等于 2 时，称为“距离平方反比法”，当幂次等于3时，称为“距离立方反比法”。
距离幂次反比法计算公式如下：
(2)
式（2）中：xb为待估块的内插品位（%），xi为已知点的品位（%），di为已知点到待估块质心点的距离 （m）；m为幂次。
本文通过交叉验证确定各矿体估值幂次，在幂次 2. 0~3. 0 之间，选择估计均值和原始均值的差异均值趋近于“0”、差异方差和差异平方均值比值趋近于“1”的幂次，最终各矿体估算幂次确定为 2（表6）。
表6估算幂次确定
5.4 资源量类型确定
根据勘查规范（中华人民共和国自然资源部， 2020b），资源量分类应根据搜索半径、工程数、样品数等因素确定。而在实际操作中，往往第三个参数的作用被淡化，用前两个参数来划分资源量类型。
本研究依据块模型品位估值时的运算次数（搜索半径）和参与估值的的工程数进行资源量分类。勘探程度的工程网度为 70 m×70 m，所以当搜索半径为 70 m 时，能够搜索到≥3个工程的块，表明块在矿体走向和倾向上均有工程控制且达到了探明资源量的要求，划分为探明资源量；当搜索半径为 70 m，只能搜索到 2个工程，或搜索半径为 140 m，能搜索到≥3 个工程的块，能达到控制资源量的要求，划分为控制资源量；当搜索半径为70 m，只能搜索到1 个工程，或搜索半径为 140 m，只能搜索到＜3 个工程，或搜索半径为280 m能够搜索到工程的块，划分为推断资源量。具体参数见表7。
根据勘查规范要求，采用地质统计学方法估算资源量时，应对矿体（矿化域）内的所有子块体品位及其他属性进行估值。但是在实际操作中，按照上述 3 次估值后，往往会出现估算域边部的少量块体无法满足估值条件而未估上值的情况，这时需要进行第 4 次估值，放宽搜索条件，例如扩大搜索半径、减少最少样品数，完成对该部分未估值块体的估值。
表7资源量划分椭球体搜索参数
块体约束后分别进行推断资源量、控制资源量和探明资源量的赋值（图12）。对于“孤岛”现象，即不同类型的资源量在空间上相互穿插或包含，应根据工程控制程度和地质可靠程度人为调整，以利于矿山建设设计和生产。
图12矿体资源量分类块体
5.5 资源量估算结果报告
根据已赋值完成的块体模型，设置资源量报告参数，分矿体、品级和资源量级别类型，以实体为约束生成资源量报表（图13）。
5.6 估算结果可靠性验证
（1）全局验证法
全局验证的方法是沿某些方向按照一定的间距划分一系列区域，分别计算落在每一个区域中参与估值的组合样品平均值和块体模型估值的平均值，绘成 2 条曲线，通过 2 条曲线的吻合度，判断估值的合理性和可靠性。
对 H23~H24 勘查线进行全局验证，共选择 13 条勘查线，7795 个组合样品点，扩展距离 35 m。通过图14可知，块体属性（估值品位）曲线与组合样品 （已知品位）曲线吻合度较好。
图13块体报告设置
图14矿体H23~H24勘查线全局验证曲线图
（2）局部验证法
对H04勘查线KZK0408钻孔的部分样品品位与其附近的块模型品位进行对比（图15），可以看出距离幂次反比法对矿块进行了较均匀的估值，估值合理。
（3）与几何法对比验算
采用距离幂次反比法估算东矿段主矿体资源量钼矿石量29625. 0万t，金属量335062 t，平均品位 0.113%。
用平行断面法与距离幂次反比法得到的资源量结果进行验证，平行断面法求得东矿段主矿体资源量钼矿石量31453.1万t，金属量344259 t，平均品位0.109%，二者对比表见表8。
两者钼矿石量相对误差为 5.81%，金属量相对误差为 2.67%，误差较小。由此可见，本文选用距离幂次反比法估算方法正确，估算结果可靠。
6 结论
（1）本文详细阐述了运用3DMine Plus软件在新疆东戈壁钼矿三维地质建模及资源量估算过程中的生产实践应用，侧重对数据库准备、圈定矿体（矿化域）、数据处理、估值参数确定过程中遇到的一些常见问题展开了论述。结果表明，利用3DMine Plus 矿业工程软件编制的资源储量报告能够满足 2020年《固体矿床资源量估算规程》要求。
图15KZK0408钻孔样品品位与矿块模型局部结果校验图
表8距离幂次反比法与垂直断面法资源量结果对比
（2）采用全局验证和局部验证，显示估值方法合理，估值结果可靠；并与传统的平行断面法估算结果进行了比较分析，二者矿石量、金属量、平均品位相对误差分别为 5.81%、2.67%、3.67%。验证了 3DMine Plus 软件在资源量估算模块中的简便性和准确性，也证明了其在地质勘查和矿山生产管理等实际地质工作中的可靠性和实用性。
注释
① 河南省地矿局第二地质勘查院.2011. 新疆东戈壁钼矿勘探报告[R].
