适于储量计算矿体建模方法.docVIP

适于储量计算矿体建模方法 　　[摘要]储量计算工作量大，计算繁琐。传统方法是在二维图形的基础上，采用几何平均的方法计算完成。根据钻孔穿过矿体具有对称性，将同一矿体相邻剖面见矿钻孔进行连接，构建不规则三角网（TIN），将TIN中的三角形逐个沿钻孔向下扩展生成GTP，在此基础上以矿体表面为约束将矿体模型进行块段剖分，建立三维栅格数据模型，为矿体内部品位插值提供基础。该建模流程简单实用，对矿体储量计算具有很强的实用价值。 　　[关键词]矿体 储量计算 三维建模 栅格 TIN 　　中图分类号：TD-0文献标识码：A文章编号：1671－7597（2009）0420109－02 　　 　　一、引言 　　 　　储量计算工作量大，计算繁琐。传统方法是在二维图形的基础上，采用几何平均的方法计算完成，其实质是以规则形体的体积近似计算矿体的体积。该方法计算精度低，不利于计算自动化。利用三维模拟技术，在矿体三维模型基础上计算储量，可以将抽象的数字信息转化为直观的、易于理解的三维图形信息传达给地质人员。 　　矿产勘查过程中，勘探工程（钻孔、探槽等）按照一定距离沿勘探线布置，矿体的空间位置和形态则根据地质成矿规律以及相邻勘探线剖面上钻孔中见矿的样品圈定，也就是说矿体边界和形态的确定以勘探线剖面图为依据。根据这一特点，利用计算机模拟三维矿体形态过程中，可以用勘探线剖面为约束条件，生成三维矿体的线框模型。首先在相邻两个剖面之间矿体的上表面构建Delaunay三角网，然后根据钻孔穿过矿体具有对称性的特点，将TIN中的三角形逐个沿钻孔向下扩展生成GTP，相邻两剖面连接示意图如图2所示。 　　GTP[3-6]是广义三棱柱（Generalized Tri-Prism,GTP）的简称，GTP是由上、下两个三角形和三个侧面四边形所组成的空间单元，如图1（a）。利用GTP进行矿体建模时，其上下底面的三角形集合可以表达矿体的表面，侧面四边形可以表达相邻矿块之间的关系，而GTP体元则可表达矿体的内部实体。针对矿体的尖灭、分叉等情况，GTP体元可以发生退化，图1（b）为一条棱边P3P6退化为一点的情况，图1（c）为两条棱边P2P5和P3P6退化为一点的情况。 　　GTP模型建立以后，通过对剖分建立块段模型，这样就建立了矿体三维栅格模型，本文根据地质钻孔数据的特点提出利用钻孔数据以相邻剖面为约束条件，生成矿体GTP模型，在此基础上建立三维栅格模型，在矿体栅格模型的基础上计算矿体储量。为地质体内部属性变量（成矿元素品位）插值估算提供了依据。 　　 　　二、三维矿体建模过程 　　 　　（一）特异值处理 　　在地质勘探领域中特异值通常指成矿元素的特高品位值，特异值对储量计算会产生很大的影响。地质勘探中的特异值是实际存在的，为了不损失信息量，不能把特异值简单的剔除，而通常采用固定的上限值代替特异值。常用的成熟处理方法主要有估计邻域法和影响系数法。本文采用估计邻域法处理特异值。估计邻域法是1980年Hawkins提出的一种识别和处理特异值的方法，计算公式如下[7]： 　　上式中I为识别特异值的统计量，服从自由度为1和∞的F分布，当I大于3.84时，可疑值G被确定为特异值，即表示在95％的置信区间上确定G为特异值。式中m为不包含G的邻域内其它样品信息值的均值，n为不包含G的邻域内的样品数， ２为邻域内观测值的平均方差。当识别出特异值后，将特异值G以下限值GL代替。GL的值采用下式计算： 　　（二）样品组合 　　为了使矿体空间品位估值更具科学性和合理性，需要对样品进行重新组合，组合后的样品品位按照长度加权平均计算。 　　（三）单工程圈定 　　单工程圈定即确定单个钻孔中矿体的边界，传统方法多采用双指标进行圈定，如果用地质统计学估算储量也可采用单指标圈定。双指标圈定方法的控制参数有边界品位、工业品位、最小可采厚度和夹石剔除厚度。单指标圈定的控制参数为经济品位、最小可采厚度和夹石剔除厚度。双指标圈定把介于边界品位和工业品位之间的矿圈为表外矿。当品位介于工业品位和边界品位之间的样品参与工程加权平均品位计算后，若平均值不低于工业品位且厚度大于可采厚度时，则这些样品圈入工业矿体，反之若平均值低于工业品位则圈入表外矿，在圈定过程中把圈入矿体的样品作为一个组合样。单工程圈定后，对每个组合样进行属性赋值，目的是确定剖面连接规则，地质人员通过地质解译能够确定在同一剖面上哪些组合样属于同一矿体，对同一矿体的组合样赋给相同的属性值。 　　（四）剖面生成 　　一个矿体在同一钻孔中可能包含多个组合样，不同组合样间为夹石，剖面的连接过程包括矿体的连接和夹石的连接。首先判断相邻的两个钻孔之间是否存在相同的矿体，如果存在则同一矿体的最顶端组合样的上端点相互连接，最底端组合样的下端点