MICROMINE培训-资源／储量估算.doc

8 资源/储量估算8.1 资源/储量估算的工业指标本矿区的工业指标是经过出春黄金设计院、吉林省第二地质调查所、天池工贸有限公司共同研究确定的，并参考了本地区敦化大石河钼矿等矿山，所采用的工业指标而确定，其指标如下：（1）矿石质量条件边界品位：0.03单工程最低工业品位：0.05（2）矿床开采技术条件最小可采厚度： 4.0米夹石剔除厚度： 8.0米经济合理剥采比： 5米3/米3当品位较高而矿体厚度达不到可采厚度时，采用米百分值（0.2米%）圈矿。8.2 资源/储量估算方法的选择及其依据本次报告没有采用传统手工资源量估算方法（地质断面法、地质块断法）进行资源量估算，而是采用了澳大利亚Micromine公司地质软件，对本矿床的钼矿体进行了圈定和储量估算，本次储量估算采用的是距离反比加权法（IDW）进行资源量估算的。 本次资源储量估算采用的软件为澳大利亚Micromine公司的三维矿产资源评价软件MICROMINE 11.03版，该软件已经通过国土资源部认证，认证书见附件。同时用该软件的封闭多面体估算法（Polygonal Section Estimate）对估算结果进行验证8.2.1 资源资源量估算的方法和原理（一）距离反比法 距离反比加权插值法（Inverse Distance Weighting）首先是由气象学家和地质工作者提出的，后来由于D.Shepard 的工作被称为谢别德法（Shepard）方法。它的基本原理是设平面上分布一系列离散点，己知其位置坐标（xi，yi）和属性值zi（i= 1，2，n）， p（x，y）为任一格网点，根据周围离散点的属性值，通过距离反比加权插值求P 点属性值。距离反比加权插值法综合了泰森多边形的邻近点法和多元回归法的渐变方法的长处，它假设P点的属性值是在局部邻域内中所有数据点的距离反比加权平均值，可以进行确切的或者圆滑的方式插值。周围点与P 点因分布位置的差异，对P （z）影响不同，我们把这种影响称为权函数W i（x， y），方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降。对于一个较大的方次，较近的数据点被给定一个较高的权重份额；对于一个较小的方次，权重比较均匀地分配给各数据点。计算一个格网结点时，给予一个特定数据点的权值，与指定方次的结点到观测点的距离倒数成比例。当计算一个格网结点时，配给的权重是一个分数，所有权重的总和等于1.0。当一个观测点与一个格网结点重合时，该观测点被给予一个实际为1.0的权重，所有其它观测点被给予一个几乎为0.0 的权重。换言之，该结点被赋给与观测点一致的值，这就是一个准确插值。权函数主要与距离有关，有时也与方向有关，若在P点周围四个方向上均匀取点，那么可不考虑方向因素，这时：式中：，表示由离散点（xi，yi）至P（x，y）点的距离。P（z）为要求的待插点的值。权函数， 本次储量估算u值取2（距离平方成反比）。（二）封闭多面体估算法封闭多面体估算法计算的步骤是，首先根据圈定的矿体模型（三角形网）的体积，按以下过程进行储量估算，估算的结果较精确。1. 确定三角网的最小Z值（最低海拔标高），将该值作为所有参与体积计算的立体三角形的基准平面； 2. 对于每个三角形，计算其与基准平面之间的体积； 3. 确定三角形和基准平面之间的体积是位于模型之内还是模型之外，通常根据每个三角形的方向来进行判断； 4. 如果在模型以内，就将其加到总体积中；如果在模型以外，就将其从总体积中减掉。 然后对模型内的所有样品使用简单平均或系数加权的方法得到总的品位和比重。如果样品在模型内间隔均匀，并且使用样长加权计算，而且选择了忽略缺失区间的话，那么三角网格模型的品位应该与块模型非常相似。如果样品间隔不 是非常均匀，并且有很多探槽和坑道的话，那么由于线框内的样品聚集，线框品位和块模型品位之间可能会存在差异。 最后，用模型的体积乘以比重得到矿石量，再用矿石量乘以品位得到金属量。 8.2.2 地质软件工作原理及流程本次采用的资源资源量估算流程见图81。图81 资源储量估算工作流程图（一）数据准备及数据格式1、准备数据内容（1）地形数据地形数据采集的是本矿区实测的1：2000地质地形图，其坐标范围是：东坐标：4258660042588600北坐标：49052004906800（2）测量数据 矿区钻孔工程井口坐标、测斜数据，编录基线拐点坐标测量数据。 （3）样品数据 矿区钻孔的样品分析数据。 （4）岩性数据 矿区钻孔的岩性数据。2、数据采集（1）地形数据 矿区的坐标系统为北京坐标系，三度带坐标，分带号为XX。本次资源量估算没有收集到原始的地形测量成果，故采用矿区已有MAPCAD/GIS格式的地形线文件数据，对地形线进行高程赋值，经校正、转换导入到MICROMINE软件中，地形数据文件数据项为：北坐标，东坐标，相对高程，线ID号。 （2）测量数据 本次采集了96个钻孔的井口坐标及测斜数据，编录基线拐点坐标测量数据。钻孔工程井口坐标数据文件，数据项为：工程编号，北坐标，东坐标，相对高程，深度。钻孔测斜文件数据项：工程编号，深度，方位角，倾角。 （3）样品数据 本次资源储量估算使用了来自矿区钻孔工程的全部样品数据，样品总数为22188件，分析元素以Mo为主。（4）岩性数据 本次岩性数据收集了矿区全部钻孔的岩性数据，其数据文件数据项：工程编号，自，至，层厚，岩石名称。 在MICROMINE中建立了钻孔数据库。应用Micromine进行矿体地质域的圈定和资源储量估算，需要三类基本数据：工程坐标文件、工程测斜文件和样品分析结果列表。此次工作所使用的数据类型齐全、相互匹配，能够满足资源储量估算的需要。（二）数据检查由原始数据输入Micromine系统后，经过了“针对工程的数据校验”、“最大值、最小值检查”和在三维视图中对勘查工程的数据校验，在进行地质体及矿体圈定过程中发现错误并及时修正，保证了在数据转抄、转换和导入过程中对错误及时更正。（三）地质解译 根据剖面端点坐标，剖面视域范围，见表81，由人工按工业指标圈定剖面（见图82），在此基础上，圈定、连接、标注各地质要素（如地层界线、岩体界线、地层或岩体代号、断层（带）界线、矿体界线等）。矿区共布置15条剖面，剖面间距100米。剖面端点坐标、方位、剖面的前后视野宽度表表81 勘探 线号左端点右端点方位限制视野宽度北坐标东坐标北坐标东坐标前后总宽度1949059674258688149064074258706523202040154905930425869744906497425872112320204011490577242587017490652342587329232020407490557342587041490640142587387232020403490535542587058490660242587580232020400490541042587190490643142587617232020404490522842587222490653542587769232020408490537942587394490639642587820232020401249053414258748649064144258793523202040164905348425875974906248425879752320204020490539342587725490608642588015232020402449053514258781549062354258818623202040284905445425879644906072425882272320204032490540942588058490603142588318232020403649053724258815149059924258841023202040图82 剖面上矿体地质域的圈连（四）建立线框模型在MICROMINE的三维视图环境中调入各剖面圈定的多边形地质体界线、矿体界线，建立三维地质实体或矿体实体。图83 在三维视图中圈连矿体地质域（五）标识采样数据 根据建立的三维矿体实体，在样品分析表中对各实体(矿体)包围的样品分别标记。图84 样品分析表中矿体的标记（六）样品组合样品组合均将取样样品的样长加权分割成1米的单元，准备进行储量估算，样长加权分割组合见下图：图85 样品分割组合（七）建立空块段模型 根据矿体实体在空间的分布范围，矿体实体在走向、倾向的变化及开采段高等因素确定矿块划分规格，将矿体实体划分为若干个立方体小块。 图86 建立空块段模型及参数（八）界定空块模型 利用建立的矿体模型，其中以剔除夹石或采空区通过的矿块，来界定空块模型。图87 界定空块模型范围（九）插值 根据各矿体实体确定的数据搜索范围用距离反比法或其它方法对划分的若干个立方体矿体块进行估值。图88 矿体块品位估值参数（十）资源量级别分类 由于按搜索椭球体的不同半径、搜索次数、工程数量等因素确定资源量级别，会出现高级别储量块体分布在工程网度低的范围内，不符合我国储量类别的划分。因此本次储量估算，是按勘探规范要求的勘探网度疏密程度，人工圈定划分储量级别的，符合我国目前勘探规范要求。（十一）资源储量可靠性验证用另外一种估算方法封闭多面体法，估算全部或部分矿体实体资源量，与所采用的估算方法所得出资源量结果进行对比，评述其估算可靠程度。图89 封闭多面体法资源量估算结果8.3 资源/储量估算参数的确定8.3.1 单工程矿体厚度 根据样品长度、钻孔方位、矿体产状，用下列厚度计算公式先计算单样品的真厚度，各样品真厚度之和为单工程矿体真厚度。厚度计算公式：m = L（cossincossincos）式中m为单样真厚度，L为样长，为矿体倾角，钻孔倾角，为钻孔方位角与矿体倾向之差；当钻孔方向与矿体倾向相反时式中取“+”，反之取“-”。各矿体厚度变化系数的厚度值，即采用此公式计算得到。8.3.2 矿体平均品位 本次采用适于机算的距离反比加权法，该方法将矿体划分为若干个立方体块，根据各立方体块的估值结果，对各小立方体块进行统计分析以各小立方体的体积加权求得矿体的平均品位。8.3.3 钼元素特高品位的确定与处理根据铜、铅、锌、银、镍、钼矿地质勘查规范（DZ/T02142002），将单样品位值高于矿床（体）平均品位六至八倍的样品确定为特高品位样。分别根据各矿体品位变化系数的大小来确定具体的特高品位样：变化系数100，品位变化程度为均匀，特高品位下限值取平均品位的6倍；变化系数100160，品位变化程度为较均匀，特高品位下限值取平均品位的7倍；变化系数160，品位变化程度为不均匀，特高品位下限值取平均品位的8倍。用各矿体中钼的特高品位下限值代替各自的特高品位。特高品位处理结果见表82特高品位处理结果一览表表82 序号勘探线号工程号样品数平均品位样 号处理情况原始值替换值1234568.3.4 样品组合根据样品的长度绘制直方图，见图810。图810 样长直方图 由直方图的结果可知样长的平均值为1.0米，该值应为样品组合的长度。由于原始样品中样长为1.0米的样品居多数，因此组合样长选择为1.0米以便尽量使用样品的原始信息。8.3.5 矿石体重本区矿石体重采用小体重法，共测试小体重样品件，计算平均体重为 2.7t/m3，详见表83矿 石 体 重 表表8-3 序号矿石类型矿石体重（吨/米3）123456789平 均8.4 矿体的圈定根据工业指标、化学分析结果、矿体产出特征、勘查类型、工程控制程度等因素综合考虑，矿体圈定原则如下：8.4.1 单工程矿体圈定本次采用单指标圈定矿体，主要考虑Mo元素，其它伴生元素品位低达不到综合利用价值。即凡样品中Mo任一品位等于或大于边界品位者均视为矿体，应用MICROMINE软件在屏幕上进行交互式圈定。1、当单工程矿体内部有大于与小于边界品位的样品相间分布，但其厚度小于夹石剔除厚度时也作为矿体一并圈入。2、夹石剔除厚度按工业指标的规定执行即矿体中连续厚度8.0米，钼品位0.03%部份作夹石，并在软件中形成夹石实体模型。此为单工程人工夹石处理部分，尚有另外一部分为软件估算的夹石各矿块钼的估值品位小于矿体边界品位的矿块的集合。8.4.2 矿体的圈定及资源量估算边界的确定 根据矿体赋存部位、产出空间位置及控矿因素、工程控制情况等全面综合分析的基础上圈连矿体及确定资源量估算边界。（一）矿体的圈定在单工程圈定的基础上，主要以矿体空间分布规律、控矿条件为主导因素，重点考虑矿体的产出位置及三度空间上的对应关系进行圈定。在软件中矿体圈定形成的边界为轮廓线（软件中一种线的数据类型）。（二）资源/资源量估算边界的确定本次在机内圈定的矿体边界与资源量估算边界为同一边界。处理原则如下所述：1、见矿工程边缘有未见矿工程控制，尖推两工程间距的二分之一为矿体边界。2、见矿工程边缘无工程控制，走向外推相邻勘探线距的一半，最多50米为矿体边界，平行勘探线方向外推100米为矿体边界，以上外推均采用尖推。8.4.3 建立实体模型矿区共建立了两类模型，分述如下：（一）矿体实体模型矿区全部矿体的实体模型，见图811图8-11 矿体空间实体及钻孔空间位置分布图（二）地面模型 在软件中建立了矿区DTM地面模型，主要用于剖面地形绘制，确定剥采地面边界。图8-12 地形实体空间位置分布图8.5 矿块模型把矿体分布的空间范围划分为小的矩形块以进行品位插值，根据勘探线距、开采段高及矿体形态的复杂程度确定块尺寸的大小，以勘探线距、开采段高为块大小的整数倍。本区勘探线距100米，段高为15米，矿体的形态复杂程度为中等，产状较缓，选择块尺寸为202015米（北东高），划分子块尺寸为101015米以保证实体模型边界的精确度。 把品位插值的搜索空间分为8个扇区，并限定每个扇区最多5个点，总点数至少2个点，防止某一方向过多的样品对一个点估值。块模型的参数见图813。图8-13 块模型参数8.6 品位插值 8.6.1 插值的元素 在矿体的块模型中对中心位用距离立方反比法进行插值，钼元素用处理特高品位后的品位。8.6.2 搜索椭球的定义 矿床的勘查类型是类型，勘查类型为类型的钼矿床的控制的勘查间距规范要求为100200100200米。根据矿床的地质特征和矿体的空间分布规律，矿床的控制的工程勘查间距确定为80200100200米。 根据矿体的总体走向、倾向、倾角确定搜索椭球的半径、方位角、倾角、倾伏角和长度因子，将基本的搜索椭球的轴定义为矿床的控制最小工程勘间距75米，将搜索椭球轴的因子定义为间距75米。图8-14 块模型和矿体品位插值的基本搜索椭球位置图对不满足基本搜索椭球条件的块，通过改变搜索椭球为基本搜索半径的2倍的方法使所有的块都能估算出钼的品位值。8.6.3 距离平方反比插值对矿体块模型用距离平方反比法按基本搜索椭球参数对Mo元素品位进行估值，如有矿块的钼品位值为空时，依次改变搜索半径为75米、150米进行搜索，直至所有的块的钼品位都估算出结果。在估值时对每一个块都记录估值次数、参与估值的工程数、样品数、和样品品位的标准离差。8.7 资源/储量分类 本次储量估算没有依据块模型品位估值时的估算次数和参与估值的工程数对块模型的地质可靠程度进行分类的，其资源类别分类与我国目前的分类有差异。因此，资源类别分类是按照我国铜、铅、锌、银、镍、钼矿地质勘查规范（DZ/T02142002）要求对资源储量进行分类的。矿区的勘探类型为类型，探明的内蕴经济资源量（331）按80100米100米工程网度估算；控制的内蕴经济资源量（332）按160200米100米工程网度估算；推断的内蕴经济资源量（333）按160200米200米工程网度估算以及外推部分。8.8 资源/资源量估算结果根据勘探类型，以及矿体控制工程的疏密,将矿体储量划分不同的资源储量类别（如331、332、333、低品位等），分别进行资源储量估算，并对块模型中钼的吨位、品位和金属量的进行了统计。其储量估算汇总表（见表84）。通过前述的资源储量估算方法，估算矿床总资源量为：矿石量31374.14万吨、钼金属量244782吨、平均品位0.078%。其中估算探明的内蕴经济资源量（331）矿石量10237.18万吨、钼金属资源量90087吨、平均品位0.088%；控制的内蕴经济资源量（332）矿石量9389.52万吨、钼金属资源量88036吨、平均品位0.094%；推断的内蕴经济资源量（333）矿石量3441.28万吨、钼金属资源量32695吨、平均品位0.095%；低品位矿石量8306.14万吨、钼金属资源量33964吨、平均品位0.041%。331+332矿石量占总矿石储量的62.55%、钼金属量占总资源储量的72.77%。季德钼矿储量估算结果表表84 储量类别体积矿石量体重品位金属量（m3）（t）（t/m3）（%）（t）331379155001023718502.70.0889008733234776000938952002.70.0948803633312745500344128502.70.09532695合计854370002306799002.70.091210818低品位30763500830614502.70.04133964总计1162005003137413502.70.0782447828.9 资源储量可靠性验证一、全局检验 对全矿段用MICROMINE软件中的距离平方反比法得到的钼品位吨位报告和用三角网多边形估算得到的钼品位吨位报告进行比较，结果见表85。（应采用其