《三维数字建模技术在某铜矿山中的应用.doc》由会员分享，可在线阅读，更多相关《三维数字建模技术在某铜矿山中的应用.doc（11页珍藏版）》请在三一文库上搜索。
1、三维数字建模技术在某铜矿山中的应用第43卷第3期2007年5月地质与勘探GEOLOGYANDPROSPECnNGVo1.43No.3May,2007三维数字建模技术在某铜矿山中的应用戴碧波,王李管,贾明涛,刘红兵,荆永滨,曾庆田(中南大学资源-9安全工程学院,长沙410083)摘要矿山三维数字建模与可视化技术研究,是数字矿山的核心组成部分,是现代矿山信息化研究的热点和重点.文章在对地质统计学和现代矿山数字建模技术研究基础上,结合国内某大型矿山,论述了矿床三维数字建模的过程和方法.结果表明,利用三维数字建模的方法来进行储量计算,品位估值,安排生产计划,将改变矿山传统的开采模式,是我国矿山现代化征
2、程中的必由之路.关键词地质统计学变异函数三维实体模型三维块体模型中图分类号P628文献标识码A文章编号04955331(2007)03009705矿床的三维数字模型是数字矿山的基础,是矿床的数字表征.采用三维数字技术建立矿床地质模型,进行矿产储量评价及矿山规划设计,对于我国矿业按国际通用方法评价矿床,适应市场经济的需要等无疑有十分重要的意义.文章在对地质统计学和英国DATAMINE软件技术的研究基础上,结合国内某大型铜矿,介绍了三维数字建模技术在矿山中的具体实践及应用.1矿山地质概况铜矿是我国一大型的地下开采矿山.共有铜,硫,铁矿体140多个,其中主矿体1个,单硫小矿体3个,未编号的零星小矿体
3、百多个.主矿体赋存于泥盆系上统五通组顶界和石炭系中上统层位中,属层控夕卡岩型铜矿床.主矿体长1810m,宽平均500m.矿体埋藏深,赋存标高为一690一1007m.矿体最大厚度100.67m,最小1.13m,平均厚34.16m,厚度变化系数80.65%,属较稳定型.矿体总体走向北东35.40.,倾向与背斜两翼产状一致,分别倾向北西和南东,矿体中部倾角较缓,一般均小于10.;而西北及东南边部较陡,一般为30.40.矿石中的金属矿物主要有磁黄铁矿,黄铁矿,黄铜矿等,主要的有用组份为铜.全矿床铜平均含量为1.01%,品位变化系数98.76%,属较均匀变化型.2DATAMINE简介Datamine是英
4、国MICL公司开发的一套广泛应用于地质,勘探,资源估计,储量计算,露天采矿和地下采矿设计的大型数字化矿山软件.它形成了一整套三维立体的和块体的建模工具,可将土建工程设计,三维模型建立,工程数据库构建等完全图形化,并解决复杂工程中境界优化的施工管理.它提供了与其他数据库和相关软件接口的功能,使该系统的数据可与其他数据库管理系统和相关软件相互查询和编辑,能够实现各种工程和矿体的三维立体显示和成图,并根据地质统计学的方法和原理提供进行矿体品位和储量估值的各种方法,同时还可以进行矿山的开采设计以及数字地形模型的建立.作为一套全面的集成软件系统,它具有先进,全面,高效,易于掌握等特点,极大地改进了测量工
5、程师,采矿工程师,地质工程师和高级管理人员之间的技术信息交流,使企业生产的各个环节在高效的管理控制之下,从而获取最大的经济效益.3建模基本参数结合现有的地质图件,钻孔数据等情况,按照矿体的形态,走向,倾向和空间分布特征确定建模的范围和基本参数.模型基本参数见表1.收稿日期20060717;修订日期20061010I责任编辑韩进国.基金项目十五国家重点科技攻关课题(编号:2004BA615Ao4)资助.第一作者简介戴碧波(1982年一),男,2004年毕业于中南大学,获学士学位,在读硕士生,现主要从事三维矿床建模研究工作.地质与勘探2007正表1模型参数表参数取值参数取值绝对模型原点坐枥m东向(
6、X):83880模型行数l10(2200m)北向(Y):21310绝对模型原点标高/ml100模型列数77(1540m)相对模型原点坐标/m东向(X):0模型层数68(1360m)北向(Y):0相对模型原点标高/m0模型转角(.)一37单元块尺寸/m202020模型旋转轴z轴4矿床三维数字建模根据datamine软件的功能和特点,实现矿山数字化技术步骤见图1.图1矿山数字化流程图4.1地质数据库及样品统计分析地质数据库由孔口坐标文件,孔口测斜文件,钻孔样品文件,钻孔岩性文件4个文件组成.它是进行地质解译,品位推估,储量计算与管理以及后续采矿设计等的基础;对样品统计分析的目的,一方面是为了掌握矿
7、床铜元素的分布情况,另一方面是指导后面品位推估时采用何种方法进行变异函数计算与分析.本次建模共收集和使用钻孔216个,有效样品数为7945个.根据地质统计学原理,为确保得到参数的无偏估计量,所有的样品数据应该落在相同的承载上,即同一类参数的地质样品段的承载应该一致.因此,在建立品位模型之前,需对样品进行组合.样品组合有多种方法,如按钻孔长度组合,按台阶高度组98合,混合组合等.本次研究选用按钻孔长度组合方法.组合样长度的确定要考虑多种因素,如原始样本长度,原始样本容量的大小,块体建模时单元块的尺寸等等.由于该矿床平均原始样品长度为1.5m,为了降低样品组合过程中可能导致的品位平均化程度,取组合
8、样长度仍为1.5m.采用按钻孔长度组合法组合样计算完成后,为了合理选择品位推估方法,确保估值精度,需进行组合样品统计分析了解其分布形式.铜组合样统计分析结果见表2.表2铜组合样品分析结果/%在品位估值时,为了消除特高品位样品的屏蔽效应,需限制样品品位区间.由于该铜矿床,品位高于6%的铜元素样品数目较少,这里把6%作为特异值,将品位大于6%的样品数据剔除.4.2.变异函数计算及拟合按球状模型用地质统计学方法计算了cu元素在走向,倾向,垂直(厚度)3个方向的变参数异函数2J.对该矿床Cu元素样品品位统计分析表明,Cu品位服从对数正态分布,因此,变异分析时,对铜样品品位进行了对数转换.变异分析参数及
9、结果见表3,实验变异函数及拟合的理论变异函数曲线见图24.表3铜品位变异函数分析参数及结果4.3交叉验证理论变异函数参数将用于后续的矿床品位推估或储量计算中,因此,理论变异函数参数取值的正误对品位估值结果的准确性具有非常大的影响.进行交叉验证的目的就是对理论变异函数参数的取值进行检验,判断应用这些参数进行品位估值时的估值效果.交叉验证的基本思想是:假定某已知样品点处的品位是不知道的,根据其周围的已知样品点,应用计算得到的理论变异函数参数对该点处的品位进行推估,然后计算该点处的实际品位与估计品位之第3期戴碧波等:三维数字建模技术在某铜矿山中的应用02O4O6O图2Cu元素沿走向方向实验及理论变异
10、函数曲线0204060图3Cu元素沿倾向方向实验及理论变异函数曲线变异函数曲线02O4O6O图4Cu元素沿垂直方向变异函数曲线间的误差.如果交叉验证的平均误差应趋近于0,误差的方差趋近于平均推估克立格方差,误差分布应属于正态分布且95%置信限应位于正负两倍的克立格方差范围内.则理论变异函数参数是准确的.交叉验证结果如表4,图5所示.由表4可知:变异函数参数比较准确,合理,可以用于进行克立格估值.Frequency.8-40图5交叉验证误差分布直方图表4交叉验证误差统计参数表误差均值0.0071误差方差平均克立格方差2个标准差范围内误差所占比例O.2O590.229895.94%4.4矿体圈定及
11、交互式修改矿体的圈定对矿床三维实体建模以及后续的矿产储量计算都起了关键性的作用.如果圈定的矿体接近实际情况,那么储量计算准确,对矿山生产具有指导作用,否则,将影响矿山的开采规模及生产.矿体圈定是在指在截取矿段J,确定了各见矿工程中的矿体边界(包括表内外储量及夹石的边界)基点之后,在勘探线剖面图上,结合地质条件及圈定原则,将这些基点合理地连接而圈出整个矿体的边界(图6).但是,由于地质条件的复杂性,地质工作人员根据实际情况可能需要对圈定的矿体边界进行调整,以反映实际生产需要.因此,通过采用datamine交互式的矿体边界点修改工具,对边界控制点的位置进行调整,以反映实际情况.图6圈定矿体边界地质
12、与勘探2007正4.5矿体三维实体建模矿体是由一系列相邻三角面,包裹成内外不透气的三维实体,其表面为不规则曲面.根据圈定的矿体剖面图,Datamine采用线框构模法.来建立矿体三维模型,线框构模法是一种表面构模技术,即把面上的点用直线连接起来,形成一系列多边形,然后把这些多边形拼接起来,形成一个多边形网格,以此来模拟矿体边界和空间形态(图7).图7矿体实体模型4.6矿体三维块体品位建模矿床的实体模型虽然给出了矿体的几何空间形态,但无法描绘出矿体内部的品位分布情况,因此必须在实体模型的基础上,建立矿床块体模型并应用地质统计学原理对模型进行品位赋值.Datamine块体品位建模的基本思想是将矿床在
13、三维空间内按照一定的尺寸划分为众多的单元块,然后对填满整个矿床范围内的单元块的品位根据已知的样品进行推估,并在此基础上进行储量的计算.块体模型属性可以是金属的品位,也可以是煤的质量或者比重等.块体的属性和图例都可以自动地用不同的颜色显示从而揭示矿体内部品位的具体分布情况(图8).本次建模参数见表5.图8Cu品位块体模型需要指出的是,搜索半径本应采用变异函数分析时计算出的变程,但由于模型规模较大,且勘探网度在首采区段以外较大,已知样品点较少,为估计出整个矿体范围内的品位,人为将其增大;但为了减少由于人为增大搜索半径导致的品位估值精度的降低1oo(在勘探网度较密区域),特意将参与估值的样品点数进行
14、限制;其他的变异函数参数采用计算得到的参数.表5铜元素品位建模参数表参数名称参数取值参数名称参数取值4.7矿石储量矿床品位模型建立后,便可以快捷地计算矿床储量,并提交储量报表(表6).同时可以根据需要计算任意边界的储量,大大节约了人力,物力和时间,提高了效率和精度.表6矿床储量统计表(按品位)1:3a.矿茹矿总计矿量:116252278.6平均品位:0.7154.8工程建模工程建模之前应将工程图纸进行数字化处理.工程图纸主要包括矿山原始的扫描图纸,还唷部分生产过程中的CAD图纸.工程图纸数字化工作主要包括两个部分:一是将原始扫描图纸在CAD环境下的矢量化;二是将已有的CAD图转换成后期所需要的
15、存储格式并进行3D转换及坐标修正,从而形成工程平剖面分布图.然后同样采用线框构模法建立工程实体模型(图9).5生产计划编制开采计划编制的主要流程包括:基础数据准备,第3期戴碧波等:三维数字建模技术在某铜矿山中的应用图9工程实体模型生成任务,流程优化与确定任务作业顺序,生产计划报表与可视化表达及动态更新与调整.基础数据的准备主要包括三维块体地质模型,掘进设计线,巷道固定断面形状,采场轮廓,不规则断面工程,设备台效,生产工序等.对设计进行分段,形成生产任务,同时根据基础数据,计算各分段任务的工程量,并为各分段任务添加各种属性与性质.流程优化指的是根据任务的工程量与设备台效确定各任务的作业时间,然后
16、根据计划技术经济指标目标,由工程的衔接关系及设备资源等约束条件,优化并确定作业顺序.生产计划报表主要包括设计统计表,虚拟评估数据报表,地质数据评估报表,品位与开采吨位统计表.同时,可由每一个任务的属性,通过一定的操作来生成各种所需的报表.输出各种数据报表,生成生产调度甘特图(图10),开采计划的生产过程动画(图l1).图10生产调度甘特图_图11开采计划动画演示中不同时期生产的状态图6结论在矿业中应用研究数字建模技术不仅十分必要,而且完全可行,这项技术的发展应用,解决了矿山地质工程建模的关键技术问题,矿床三维数字建模的研究和应用是一个复杂的系统工程,是地质学,数学,计算机技术,信息技术,采矿技
17、术相互融合的产物.利用三维数字建模的方法建立矿体模型可为进行储量计算,品位估值,安排生产计划,实时调度监控提供很好的平台,这样不仅对传统的管理手段提出了更高的要求,也将对我国矿山企业的数字化进程起到巨大的推动作用.参考文献1侯景儒,黄竞先.地质统计学及其在矿产储量计算中的应用M.北京:地质出版社,1982.2贾明涛,潘长良集成可视化矿床建模软件DMS在某矿山的应用J.中南工业大学.2000,31(5):396399.3贾明涛,王李管.三维变异函数的稳健统计学计算方法及其应用J.中南工业大学,1998,29(4)570574.4李艳,王恩德,鲍玉斌,等.基于钻孔数据的矿体三维可视化研究与实现J.
18、沈阳工业大学.2005,27(4):418421.5李赋屏,蔡劲宏,任建国.矿业软件在矿产储量评价中的应用J.桂林工学院.2005,25(1):2630.6南格立矿体线框模型及其建立方法J.有色矿山,2001,(5):14.7王李管,何昌盛,贾明,涛.三维地质体实体建模技术及其在工程中的应用.金属矿山,2006,356(2):5862.8陈爱兵,秦德先,张学书,等.基于MICROMINE矿床三维立体模型应用fJ.地质与探勘,2004,40(5):7780.APPLICATIoNOlF3DDIGITALMoDELGINACoPPERMINEDAIBibo,WANGLiguan,JIAMingta
19、o,LIUHongbing,JINGYongbin,ZENGQingtian(CollegeofResourcesandSecurityEngineering,CentralSouthUniversity,Changsha410083)Abstract:3DdigitalmodelingandvisualizationofminesisanimportantpartofDigitalMine,andalsoresealhhotspotofmodemmineinforma-tionization.Basedongeostatistiesanddigitalmodelingtechnologyin
20、modemmines,combinedwithalargescaledomesticmine,processandmethodof3Ddigitalmodelingarediscussed.Itisprovedthatusing3Ddigitalmodelingmethodtocalculatererve,estimategrade,scheduleproductionwillchangetraditionalexploitationmodeofmine,andisindispensabletOachieveminemodernization.Keywords:geologicalstatistics,va/lcefunction,3Dsolidmodel,3Dblockmodel,productionscheduling101