《武汉工程大学学报》 2013年09期
34-38
出版日期：2013-10-10
ISSN:1674-2869
CN:42-1779/TQ
大型露天磷矿数字矿床系统
0引言 进入21世纪以来，随着科技进步、矿产资源消费的急剧增长和开采加工难度的日益增大，采矿逐渐走向数字化和智能化.“数字矿山”作为信息时代的必然产物，在矿业中发挥出越来越大的作用，是矿业发展的目标和方向［1\|3］.在20世纪70年代，人们开始将AUTOCAD技术二次开发进行三维地质建模及可视化技术研究.20世纪90年代后国际上出现了一大批三维矿山专业软件，主要有MicroMine、Surpac、DataMine、Vulcan、MineScape、MineSoft、MineSight等，而国内的一些矿山企业和科研单位从国外引进了MicroMine、Surpac、DataMine等并进行了二次研究开发［4\|7］，如崔灿［8］、石慧［9］等利用Surpac软件对矿山数字化进行了研究，钟德云［10］、徐恒雷［11］、郭强［12］等利用3DMine软件建立了矿山的三维数字模型，肖英才［13］、蒋京名［14］对DIMINE软件工程实体建模进行了研究. 笔者依据云南磷化集团昆阳磷矿地质报告、地形地质CAD图、剖面图等相关资料，研究昆阳磷矿矿床分布空间关系，建立“孔口表”、“测斜表”、“样品表”等数据信息表，依据DIMINIE软件地质数据库的要求检验其是否符合格式要求或是否在录入过程中出现进行错误校正，然后依次将上述信息导入DIMINE软件，生成完整的地质数据库.最后利用地形地质CAD图和钻孔数据库建立昆阳磷矿四采区的地形模型、三维实体模型，通过块段模型，并对不同品级矿石赋予各块段不同的体重值，运用距离幂估值地质统计法进行储量估算，实现资源储量估值.1昆阳磷矿概况 昆阳磷矿位处滇池西南部的香条村背斜南翼东段，矿床类型为寒武纪早期形成的沉积型磷块岩矿床.矿层总体呈向南缓倾斜的单斜构造.矿层（体)分为上、下两层矿，中间为夹层.上矿层工业矿体厚5.00～7.00 m，最厚15.30 m，平均厚度6.18 m，厚度变化系数为34.60%，属稳定型，厚度变化小；下矿层工业矿体厚度3.00～4.00 m，最厚7.95 m，平均3.25 m，厚度变化系数57.70%，属较稳定型，变化中等；夹层厚度一般0.5～2.0 m，东部最厚可达7.30 m，西部最薄0.05 m，全区层位稳定. 矿石自然类型主要有白云质磷块岩、条带（纹）状磷块岩、粒屑磷块岩、含硅质条带磷块岩、砾状磷块岩和生物碎屑磷块岩等六个自然类型. 昆阳磷矿Ⅰ品级矿石工业类型以混合型为主，少部分为硅质及硅酸盐型；Ⅱ品级矿石工业类型以碳酸盐型及混合型为主，少部分为硅质及硅酸岩型;Ⅲ品级矿石工业类型以碳酸盐型为主，部分为混合型.2矿山三维地质建模2.1地质数据库 地质数据库是三维地质建模的重要组成部分，地质勘探数据是用来进行建模的基础资料，完备的地质数据库是进行地质解译、品位推估、储量计算与管理以及后续采矿设计等的基础，其准确性以及多数据信息的变异性分析是非常重要的.昆阳磷矿地质数据库数据主要来源有：昆阳磷矿钻孔资料、钻孔编录资料及钻孔综合成果表、昆阳磷矿地形地质图等. 地质数据库的建库思路是首先收集整理建立地质数据库所需要的大量数据资料，并根据DIMINE软件系统格式将地质报告中钻孔数据按照“孔口表”、“测斜表”、“样品表”等分别录入不同的EXCEL表格文件中.对于已经录入的数据信息，依据DIMINIE软件地质数据库的要求检验其是否符合格式要求或是否在录入过程中出现进行错误校正.然后依次将上述信息导入DIMINE软件，生成完整的地质数据库.2.2地表地形模型 地表地形模型是用来虚拟地形和表面，一般由若干地形线和地形特征散点组成，考虑到每个点的坐标值，将所有点联成若干相邻的三角面，然后形成一个三角网随着地面起伏变化的单层模型. 昆阳磷矿三维地表模型的建立过程如下：第9期柴修伟，等：大型露天磷矿数字矿床系统武汉工程大学学报第35卷①先在CAD中将地形等高线、坐标网和等高线的值描绘和标记出来，将处理好的CAD图形再导入到DIMINE软件系统中，然后在DIMINE视图管理器环境下把平面坐标转换为它的相对地理坐标. ②对等高线进行赋值. ③在DIMINE视图管理器环境下，进行创建数字地面模型（Digital Terrain Model，以下简称：DTM）的工作，生成的DTM也可进一步在等高线之间进行加密，生成比较圆滑的三维地形模型，如图1所示.图1昆阳磷矿四采区数字地表模型Fig.1Digital terrain model of forth mining area of Kunyang mine2.3三维地质模型 三维地质模型的建立一般有两种方法，一是通过CAD或MapGis对地质剖面图进行矢量化，导入DIMINE进行建模；二是通过调入钻孔数据库，并对矿体边界作地质解释，然后利用剖面圈定矿体，对昆阳磷矿四采区的三维地质模型采用的是第二种建模方式. 四采区的三维地质实体模型包括地形模型（含公路模型）、断层模型、矿体实体模型等，通过调用全部模型或部分模型，可以显示各个模型之间的总体关系或局部关系，实现矿床三维可视化表达的目标.图2是昆阳磷矿四采区三维地质模型.图2昆阳磷矿四采区三维地质模型Fig.23D geologic model of forth mining area of Kunyang mine3资源储量估算 资源储量估算是通过已知采样点数据，利用一定的数学方法如距离幂次反比法、普通克里格法等，对预测矿体进行空间插值，以便求得未知矿体的平均品位、体积、吨位和矿石量及金属量等一些特性.DIMINE软件利用块段方法通过对块段模型的子块进行估值来实现资源储量估算，估算方法有最近距离法、距离幂次反比法和克里格法等. 块段模型是矿床品位推估及储量计算的基础，块段模型的基本思想是将矿床在三维空间内按照一定的尺寸划分为众多的单元块，然后对填满整个矿床范围内的单元块的品位根据已知的样品进行推估，并在此基础上进行储量的估算.采用距离幂估值，对昆阳磷矿四采区的上层矿（1z1\|3)和下层矿（1z1\|1)块段模型应用“模型配色”功能，以P2O5品位高低配色，效果如图3所示.图3距离幂估值后块段模型P2O5品位配色效果图（左：1z1\|3,右：1z1\|1)Fig.3Block model with P2O5 grade color matching after power distance valuation（L：1z1\|3,R：1z1\|1)在各块段模型中填加“体重”字段，并采用函数计算方式，根据不同品级矿石赋予各块段不同的体重值.块段模型中体重字段的赋值函数如下：风化带上层矿： if(above equal(P2O5,25),if(above equal(P2O5,30),2.17,2.168),2.664) 风化带下层矿： if(above equal(P2O5,25),if(above equal(P2O5,30),2.212,2.249),2.696) 原生带上矿层： if(above equal(P2O5,25),if(above equal(P2O5,30),2.456,2.597),2.664) 原生带下矿层： if(above equal(P2O5,25),if(above equal(P2O5,30),2.665,2.599),2.696) 赋予体重属性后，应用地质统计法进行储量估算，分高程段统计，对上层矿（1z1\|3)和下层矿（1z1\|1)矿体按中段估算资源储量，估算结果以Excel表格显示.由于数据表格较大，仅截取上层矿（1z1\|3）1980中段数据，如表1所示.表1上层矿（1z1\|3) 1980中段资源储量统计表Table 1Statistical table of resource reserve of 1z1\|3, 1980 midpiece中段 1980矿石品级 0.00-8.00 8.00-15.00 15.00-25.00 25.00-30.00 30.00-100.00 汇总估值距离 无 无无 无 无 无体积/m3 164 096 1 568 160 191 824 0 0 1 924 080体重/（t/m3） 2.664 000 034 2.664 000 034 2.664 000 034 0 0 1.598 400 021矿石量/t 437 151.75 4 177 578.25 511 019.156 3 0 0 5 125 749.156品位 P2O5/% 5.553 991 795 11.387 274 74 16.003 721 24 0 0 6.588 997 555 Fe2O3/% 0.862 469 852 0.845 281 124 0.983 616 948 0 0 0.538 273 585 AI2O3/% 0.954 506 218 0.855 262 995 0.960 224 032 0 0 0.553 998 649 MgO/% 5.566 897 869 5.165 808 201 2.971 955 776 0 0 2.740 932 369 CaO/% 36.238 269 81 36.868 789 67 36.651 008 61 0 0 21.951 613 62 SiO2/% 19.063 381 2 18.987 707 14 22.448 1620 8 0 0 12.099 850 08 F/%0.0179 891 01 0.236 335 412 0.000 241 865 0 0 0.050 913 276 CO2/% 12.032 241 82 11.050 330 16 6.955 164 433 0 0 6.007 547 283续表1中段 1980金属量 P2O5/t 24 279.371 09 475 712.312 5 81 782.078 13 0 0 581 773.761 7 Fe2O3/t 3 770.302 002 35 312.281 25 5 026.470 703 0 0 44 109.053 96 AI2O3/t 4 172.640 625 35 729.281 25 4 906.928 711 0 0 44 808.850 59 MgO/t 24 335.791 02 215 805.671 9 15 187.262 7 0 0 255 328.725 6 CaO/t 158 416.234 4 1 540 222.5 187 293.671 9 0 0 1 885 932.406 SiO2/t 83 335.906 25 793 226.312 5 114 714.398 4 0 0 991 276.617 2 F/t 78.639 671 33 9 873.096 68 1.235 976 696 0 0 9 952.972 328 CO2/t 52 599.152 34 461 636.218 8 35 542.222 66 0 0 549 777.593 8通过对不同矿层矿石分高程段进行储量估算，可清楚了解不同水平中段的矿石体积﹑体重﹑矿石量等基本参数，以及矿石中各组分的品位及其金属量等特征，为矿山资源的可持续管理提供重要的参考依据.4结语 a.以DIMINE软件系统为平台，依据昆阳磷矿钻孔资料、钻孔编录资料及钻孔综合成果表、昆阳磷矿地形地质图等，建立了昆阳磷矿地质数据库，为三维地质建模和储量估算提供了基础，并可对钻孔数据库进行新数据的录入、编辑，使数据库不断更新. b.通过软件系统提供的一些工具和软件实体的布尔运算功能，建立了矿山地表地形模型和三维地质模型，实现昆阳磷矿四采区的三维可视化功能，进而可以实现对三维模型不同方向的剖面切割. c.通过建立块段模型，并对不同品级矿石赋予各块段不同的体重值，对赋予组分、体重等信息的块段模型，运用距离幂估值地质统计法进行储量估算，分高程段统计，对上层矿（1z1\|3)和下层矿（1z1\|1)矿体按中段估算，实现资源储量估值.致谢 本研究得到了云南磷化集团有限公司磷资源开发利用工程技术研究分公司和昆阳磷矿的支持与帮助，在此表示衷心感谢！