img

QQ群聊

img

官方微信

  • CN 62-1112/TF 
  • ISSN 1005-2518 
  • 创刊于1988年
高级检索

黄金科学技术, 2019, 27(3): 458-465 doi: 10.11872/j.issn.1005-2518.2019.03.458

采选技术与矿山管理

基于混合整数规划法的露天矿配矿品位控制

涂鸿渐,1,2, 王李管,1,2, 陈鑫1,2, 任助理1,2, 张炬1,2

1. 中南大学资源与安全工程学院,湖南 长沙 410083

2. 中南大学数字矿山研究中心,湖南 长沙 410083

Ore Grade Control Blending in Open-pit Mine Based on Mixed Integer Programming

TU Hongjian,1,2, WANG Liguan,1,2, CHEN Xin1,2, REN Zhuli1,2, ZHANG Ju1,2

1. School of Resources and Safety Engineering,Central South University,Changsha 410083,Hunan,China

2. Center of Digital Mine Research,Central South University,Changsha 410083,Hunan,China

通讯作者: 王李管(1964-),男,山西乡宁人,教授,从事数字矿山研究工作。liguan_wang@163.com

收稿日期: 2018-10-08   修回日期: 2019-02-22   网络出版日期: 2019-07-08

基金资助: 国家重点研发计划项目“深部金属矿集约化连续采矿理论与技术”.  2017YFC0602905

Received: 2018-10-08   Revised: 2019-02-22   Online: 2019-07-08

作者简介 About authors

涂鸿渐(1994-),男,江西上饶人,硕士研究生,从事数字矿山研究工作hongjian_tu@163.com , E-mail:hongjian_tu@163.com

摘要

为了解决露天矿配矿耗时长、成果精度低及与卡调系统不关联的问题,对露天矿配矿进行系统科学分析,提出了基于混合整数规划的露天矿自动化配矿系统,实现对露天矿品位的快速准确控制。以炮孔岩样化验结果和矿山块段模型数据为基础,运用泰森多边形方法对配矿单元进行估值划分,综合考虑矿量、品位、设备和空间位置等约束条件,建立了以出矿品位波动最小为目标函数的露天矿配矿混合整数规划模型。并在Visual Studio环境下用C++语言调用CPLEX求解器进行求解,进而得到用于卡调系统的露天矿配矿开采和派车优化方案。该方法应用于安徽某大型石灰岩露天矿,经实例验证,此方法可在1 s内输出配矿结果,同时配矿指标均在目标值上下1.5%范围内波动,对指导矿山的实际生产有重要意义。

关键词: 露天矿 ; 配矿 ; 混合整数规划 ; 品位控制 ; 石灰石 ; 块段模型 ; 泰森多边形法

Abstract

Rational planning and management of ore quality is an important issue to improve the resource and economic benefits of the mine.Resources and quality control of ore should be run through all aspects of the entire production management process.In order to solve problems existed in ore blending for open-pit mine,such as long-time ore blending,low accuracy of planned results and not associated with truck dispatching system,the open-pit mine automation ore blending system was proposed.It is based on Mixed Integer Programming and would make the ore grade control in open-pit mine more accurately and rapidly,so rational mining and fine control of production targets for the open-pit mine automation ore blending would be achieved.It was put forward based on the current status quo of stope,exploration engineering data and resource model,on the premise of production quality control,and mixed integer programming algorithm is used as the theoretical method.This paper taking a large limestone open-pit in Anhui Province as research object.In the first place,thiessen polygon method was applied to estimation value and divide the ore blending model based on rock sample of blasthole and mine block model.Then a mixed integer programming model was presented according to the special ore blending requirements of muck piles in open pit mineto solve the optimization problem of ore blending.The model takes the minimum fluctuation of the ore grade as the objective function,and takes into account a variety of constraints,such as the muck sub-piles mining constraint,vehicle dispatching manner constraint,shovel car mining constraint,number of shovel car in each pile constraint,maximum grade of ore dumping constraint,minimum grade of ore dumping constraint .At last,the automatic ore distribution system of open pit mine is established based on iMine spatial database technology with C++ language programming in Visual Studio 2012,by solving the mixed integer programming model of open-pit ore blending by IBM ILOG CPLEX Optimization Studio, the open-pit ore blending optimization scheme of mine and truck dispatching can be achieved.The method was applied to Quanjiao Conch limestone mine,several months of practical application proved that this model is scientific and feasible.The ore blending method can output result in 1 s,and the maximum deviation of each index of ore blending fluctuates from the target value within the range of 1.5%,the average deviation of each index from the target value is less than 0.2%,and the coefficient of variation is less than 1%.It meets the requirement of mine grade deviation control,so it is of important significance for guiding the actual production of open-pit mine.

Keywords: open-pit ; ore blending ; mixed integer programming ; grade control ; limestone ; block model ; thiessen polygon method

PDF (3180KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

涂鸿渐, 王李管, 陈鑫, 任助理, 张炬. 基于混合整数规划法的露天矿配矿品位控制[J]. 黄金科学技术, 2019, 27(3): 458-465 doi:10.11872/j.issn.1005-2518.2019.03.458

TU Hongjian, WANG Liguan, CHEN Xin, REN Zhuli, ZHANG Ju. Ore Grade Control Blending in Open-pit Mine Based on Mixed Integer Programming[J]. Gold Science and Technology, 2019, 27(3): 458-465 doi:10.11872/j.issn.1005-2518.2019.03.458

矿山生产计划是矿山工作的重要组成部分,其目标是制定符合某些预定目标的开采计划[1,2]。根据计划时间跨度可将其分为长期、中期、短期生产计划。矿山生产计划应考虑生产能力、采区顺序、回采约束和品位约束等矿区可能存在的限制条件[3]。由于矿山工程数据在精度和完整性上存在不足,采矿工程领域充满不确定性。例如,由于空间品位存在不确定性导致开采计划的制定变得极具挑战[4,5,6]。这些不确定性使合理的开采计划变得更加重要[7,8]

配矿是矿山在制定生产计划时为实现矿石质量均衡并达到选厂品位要求而采用的重要措施与手段。为了实现科学高效的生产配矿,许多学者对露天矿的配矿计划进行了研究。例如,将模拟退火法和线性规划模型应用于矿山短期生产计划编制和生产调度中[9,10],提出了启发式露天矿短期计划调度方法[11]和用于多采场卡车分配的启发式算法[12]。除应用数学方法之外,一些研究人员已经将模拟仿真应用于短期计划配矿优化中,包括模拟仿真技术在露天矿山生产计划表制定[13]和短期生产计划评估中的应用[14,15]。国内的配矿研究也取得了诸多成果,例如:黄启富等[16]使用粒子群算法在矿山配矿中实现了降本增效;杨珊等[17]运用运筹学优化软件Xpress-MP对配矿模型进行计算,确定出最优的配矿方案;王李管等[18]和吴丽春等[19]运用线性规划和目标规划制定露天矿配矿计划。

然而上述研究多数仍处于试验研究阶段,难以满足实际生产的需求,特别是石灰石矿配矿对开采回收率、采矿贫化率等要求严格,需要更精确的配矿计划,并且需要将结果与卡调系统实现实时对接,快速有效地调整配矿计划以满足矿山生产的需要。针对此种情况,运用混合整数规划法,选取目标函数及有关约束条件,结合矿山现场真实数据,构建适用于露天矿品位控制的配矿计划模型。

1 配矿原理与流程

1.1 石灰石配矿原理

配矿是指在矿山生产过程中,根据出矿质量标准和矿床资源综合利用经济指标,按计划搭配各品位的矿石,使混合品位达到要求,从而提高矿山经济效益的措施和手段。

在石灰石配矿中,参与决策的因素主要有爆破作业后可以进行回采工作的爆堆及其分组、在作业地点出矿的电铲及数量、运输矿石的卡车及数量和卸矿后进行矿石混合的卸点情况。其中,爆堆具有不同的矿石品位,铲车和卡车有固定的数量和一定的生产能力,卸点需要混合矿石品位满足矿山的要求。因此,矿山在生产时需要通过配矿来决定具体的派车方案,使得在计划时间内的出矿混合品位达到目标要求[20,21,22]

1.2 配矿单元划分与估值

在进行配矿计划编制之前,建立合理有效的配矿单元品位分布模型是确保配矿计划可行性的必要条件。在露天采场生产过程中,通常会对每个爆堆分组的炮孔爆破岩样进行化验分析,结合生产勘探过程中得到的块段模型品位信息,应用泰森多边形法对配矿单元进行划分和估值[23,24,25]

在配矿单元中,依据炮孔现场采样分析和荧光分析仪化验结果的编号划分爆堆分组,同时按各分组中炮孔点的位置将配矿单元所在的数据平面分割成若干个由泰森多边形表示的子区域[26,27,28]。其中,每个泰森多边形是由相应的炮孔点与周围的所有邻域炮孔点间作垂直平分线后将各垂直平分线依次连接组合而成,由此通过离散的炮孔化验品位可估计各爆堆分组中的元素品位,划分估值后的配矿单元模型如图1所示。

图1

图1   配矿单元模型

Fig.1   Element model of ore blending


1.3 露天矿配矿混合整数规划模型构建

(1)目标函数。根据露天矿生产配矿的原则,选取配矿周期内出矿品位波动最小为目标函数,以确保出矿品位均匀,进而降低后续工作成本,提升企业经济效益,其目标函数的数学表达式为

MinmizeeEλe×(ge-+ge+)

式中:e为矿石参与配矿元素的索引;E为金属元素的集合,eEge-0为元素e小于需求品位的偏差值;ge+0为元素e大于需求品位的偏差值;λe>0为元素e的权重系数,当进行多元素配矿时,权重系数决定了元素满足选矿品位的优先级。

(2)决策变量。通过对露天矿爆堆出矿过程进行分析,以露天矿爆堆和分组是否出矿为配矿计划模型的二元决策变量,选取第i爆堆第j分组第s种铲车第n种出矿方式(xi,j,s,n),元素e品位小于最小品位要求的偏差(ge-)及元素e品位大于最大品位要求的偏差(ge+)为配矿模型的决策变量,定义当第i爆堆第j分组第s种型号铲车第n种出矿时变量xi,j,s,n取值为1,否则值为0;第i爆堆第j分组出矿时变量xi,j取值为1,否则值为0。

xi,j,s,n=1,ijsn0
xi,j=1,ij0

(3)约束条件。①变量约束,二元决策变量,松弛变量非负。

xi,j,s,n=0,1iD,jG,sS,nNs
xi,j=0,1iD,jG
ge-0eE
ge+0eE

式中:D为爆堆(泛指爆堆或者堆场)的集合;G为分组(索引到每个爆堆)的集合;S为铲车编号的集合;Ns为铲车s配矿方案的集合。

②分组出矿和派车方式约束,每个分组最多只有一种铲车出矿方式。

snxi,j,s,n=xi,jiD,jG
xi,jnxi,j,s,niD,jG,sS

③铲车出矿约束,每个铲车只能在一个分组出矿。

ijnxi,j,s,n=1iD,jG

④每个爆堆最多铲车约束。

jxi,jKiiD

式中:Ki为爆堆i最多分配的铲车数。

⑤卸矿点的最大、最小品位约束。

ijsnxi,j,s,n×λj×DGs,n×GGi,j,eijsnxi,j,s,n×λj×DGs,n+ge-ge̲eE
ijsnxi,j,s,n×λj×DGs,n×GGi,j,eijsnxi,j,s,n×λj×DGs,n-ge+ge¯eE

式中:λj为每个分组j的装矿条件系数;DGs,n为第s种铲车第n种的出矿量;GGi,j,e为爆堆i分组j元素e品位集合;ge¯ge̲为卸矿点对元素e的最大、最小品位。

⑥卸矿点的出矿量约束。

dm̲ijsnxi,j,s,n×λj×DGs,n×Tdm¯

式中:T为配矿周期;dm¯dm̲分别为卸矿点的最大、最小出矿量。

本配矿数学模型是针对多元素的通用型配矿模型,考虑了元素配矿的优先级,以各元素品位偏差之和最小为目标函数,当短期计划出现偏差时,即使配矿条件达不到品位指标,也能够根据现有约束求得最接近品位要求的配矿方案。

1.4 配矿系统流程

配矿系统流程如图2所示,具体步骤为:(1)生成爆堆分组模型,根据炮孔化验分组数据和爆破参数生成爆堆分组,并计算每个爆堆分组的矿量及各元素的品位;(2)设置回采参数和设备信息,包括卸矿点目标出矿量、卸矿点目标品位、铲车数目、铲车生产能力和爆堆允许派车范围等,确定爆堆和出矿相关约束;(3)建立混合整数规划模型,结合爆堆分组及相关设置参数,在配矿系统生成配矿数学模型;(4)求解配矿混合整数规划模型,计算目标函数的品位偏差,若无最优解根据需求返回修改相关参数;(5)输出最优解作为最佳配矿方案并实时监测,当计划执行出现偏差时自动调整并重新修改规划模型。

图2

图2   配矿功能系统流程

Fig.2   Flow diagram of ore blending system


2 工程应用

2.1 工程概况

某大型石灰石矿位于全椒县十字乡310°方向约10 km处,距全椒县城335°方向约14 km,距滁州市城区约17 km。如图3所示,该矿为露天开采,年采矿能力为600×104 t,配料能力为30×104 t。矿区采场狭长,矿体中间品位高,两边品位低,品位上下变异小,平面略大,采场面积约为1 km2。矿山主要设备有钻机1台,挖掘机4台,运输车辆13台。目前有生产工作面5个,每天正常采矿生产面3个,矿山每天早上按生产指标进行生产车辆配比,调度车辆分别到挖掘机区域运矿,车辆完成装载后运送到下料口进行破碎。

图3

图3   某大型石灰石矿山三维地质图

Fig.3   3D geological map of one large limestone mine

1.地形等高线;2.采场坡顶线;3.采场坡底线;4.三维地表;5.采场圈定范围


2.2 数据准备

根据该矿的情况,利用现场采集的炮孔数据,生成当日参与配矿的BD2017-15006N,BD201709-11227S,BD201709-11229S,BD201709-11230S,BD201711-10006S,BD201711-10022M,BD201711-11201M,BD201711-11210S等共25个爆堆,如图4所示。

图4

图4   配矿爆堆分布图

Fig.4   Muck piles distribution map of ore blending

1.采场坡顶线;2.采场坡底线;3.高镁爆堆分组;4.高钙爆堆分组;5.高硅爆堆分组;6.泥夹石爆堆分组


根据炮孔化验分组对每个爆堆内部按顺序进行分组,矿石体重为2.7 t/m3,利用炮孔数据进行估值,表1为估值完成后各爆堆分组矿量品位信息。

表1   爆堆分组矿量品位信息

Table 1  Ore quantity and grade information of muck sub-piles

爆堆配矿单元CaO/%MgO/%R2OSM矿石量/t
BD2017-1500060145.1130.7960.2625.3915 837.47
BD201709-11227S0153.1220.9210.1453.48916 443.42
0252.4130.8920.1373.18214 902.94
BD201709-11229S0142.1042.1180.7533.49712 839.54
0251.5271.3040.1744.0045 954.49
0353.8580.8850.1193.60616 981.37
BD201712-15012M0142.3450.9760.4605.74215 983.21
0247.4740.7070.4115.14313 921.77
0348.4171.2790.4333.28512 642.13

新窗口打开| 下载CSV


2.3 模型求解

设置出矿目标矿量为(21 000±1 000)t/班,目标品位CaO为46.20%、MgO为1.48%、R2O[0.658×K2O+Na2O]为0.45、SM[SiO2/(Al2O3+Fe2O3)]为3.40,铲运设备PC400反铲能力为380~420 t/h、CAT988装载机能力为450~500 t/h、PC1250正铲为500~600 t/h。针对配矿数据,利用C++语言编写露天矿配矿优化的混合整数规划模型,在Microsoft Visual Studio 2012中调用CPLEX采用分支定界算法进行快速求解,操作环境主要为Microsoft Visual Studio 2012,IBM ILOG CPLEX Optimization Studio V12.7.0,Windows 7 Build 7700.0.winmain.100122-1900,Intel Core i5-4590,Corsair 4GB。为了使运行结果更加清晰,对模型进行求解一周的配矿结果,模型平均求解时长为0.825 s。表2为某日实际配矿方案,包含爆堆分组位置、卸矿位置、铲运设备、出矿品位和出矿量等。

表2   配矿方案

Table 2  Ore blending scheme

铲车型号爆堆配矿单元矿石类型CaO/%MgO/%R2OSM矿石量/t铲车/辆
总计石灰石46.2211.4800.4543.37821 110.00703
01#PC400BD2017-15006N01高硅45.1130.7960.2625.3913 040.00101
02#PC400BD201712-10009S01高镁44.0012.5070.2643.0003 300.00110
03#PC400BD201712-10014M01高镁37.8372.6231.1922.6093 360.00112
05#PC400BD201712-15012M03高钙48.4171.2790.4333.2853 360.00112
CAT988BD201709-11230S02高钙53.5320.6810.1162.5454 000.00133
PC1250BD201712-10022S02泥夹石46.8171.1850.4993.7084 050.00135

新窗口打开| 下载CSV


2.4 结果分析

统计一周内主要指标CaO,MgO,R2O[0.658×K2O+Na2O],SM[SiO2/(Al2O3+Fe2O3)]的数值波动情况,如图5所示。根据结果计算各成分的标准离差率,如图6所示。

图5

图5   配矿指标数值波动图

Fig.5   Numerical wave diagram of ore blending target


图6

图6   配矿标准离差率图

Fig.6   Coefficient of variance diagram of ore blending


图5可知,单次配矿结果主要指标CaO、MgO、R2O和SM的最大偏差分别为0.1083%、1.0815%、1.2892%和0.7064%,可见指标平均数值越大,模型所得计划值与目标值偏差越小,各指标偏差均小于2%,符合出矿要求。图6所示为配矿指标平均值与目标值之比以及各指标配矿标准差所得变异系数。由图6可知,CaO、MgO、R2O和SM的平均偏差值分别为0.0075%、0.1636%、0.0612%和0.0191%,变异系数分别为0.0702%、0.6413%、0.4421%和0.6095%,表明配矿方案中各指标值均在目标值上下小范围波动,具有较高的准确性和较好的稳定性。

3 结论

(1)提出了一种用于露天矿山配矿计划编制的混合整数规划模型。在总结前人研究成果的基础上,综合考虑生产运输的全过程,将配矿中存在的相互约束融入到模型中,使模型相较于目前配矿方法更加精细、普适性更好。

(2)在构建模型时,考虑了露天矿生产的基本原则,保证爆堆和设备的合理分配,在保证矿山出矿量满足要求的情况下,选择配矿时期内矿石品位波动范围最小作为目标函数,使放矿计划得以优化。

(3)将模型运用于全椒海螺石灰石矿配矿计划编制中,以班次为单位统计量,在保证出矿需求的情况下,1 s之内就能得到合理的配矿方案,同时将结果上传至卡调系统,极大地提高了矿山配矿效率。

(4)设备的行进路线对露天矿配矿结果具有重要影响,在本模型中将设备与配矿单元在空间中直接关联,后续将进一步研究基于采场线路的配矿方案,使配矿结果更加精确,同时为保证模型在大规模规划中的配矿效率,后续将把GPU加速引入露天矿配矿模型求解中。

参考文献

Espinoza D Goycoolea M Moreno E et al .

MineLib:A library of open pit mining problems

[J].Annals of Operations Research,20132061):93-114.

[本文引用: 1]

Rahmanpour M Osanloo M .

Resilient decision making in open pit short-term production planning in presence of geologic uncertainty

[J].International Journal of Engineering,Transactions B:Applications,2016297):1022-1028.

[本文引用: 1]

Caccetta L Hill S P .

An application of branch and cut to open pit mine scheduling

[J].Journal of Global Optimization,2003272/3):349-365.

[本文引用: 1]

Azimi Y Osanloo M Esfahanipour A .

An uncertainty based multi-criteria ranking system for open pit mining cut-off grade strategy selection

[J].Resources Policy,2013382):212-223.

[本文引用: 1]

Gholamnejad J Osanloo M .

Incorporation of ore grade uncertainty into the push back design process

[J].The Journal of the Southern African Institute of Mining and Metallurgy,20071073):177-185.

[本文引用: 1]

Rahmanpour M Osanloo M .

Determination of value at risk for long-term production planning in open pit mines in the presence of price uncertainty

[J].The Journal of the Southern African Institute of Mining and Metallurgy,20161163):229-236.

[本文引用: 1]

Newman A M Rubio E Caro R et al .

A review of operations research in mine planning

[J].Interfaces,2010403):222-245.

[本文引用: 1]

Osanloo M Gholamnejad J Karimi B .

Long-term open pit mine production planning:A review of models and algorithms

[J].International Journal of Mining,Reclamation and Environment,2008221):3-35.

[本文引用: 1]

Kumral M Dowd P A .

Adelaide research and scholarship:Short-term mine production scheduling for industrial minerals using multi-objective simulated annealing

[C]//Bandopadhyay S.International Symposium on the Application of Computers and Operations Research in the Minerals Industry.Colorado:Society for Mining,Metallurgy and Exploration,2002731-741.

[本文引用: 1]

Desaulniers G Langevin A Riopel D et al .

Dispatching and conflict-free routing of automated guided vehicles:An exact approach

[J].International Journal of Flexible Manufacturing Systems,2003154):309-331.

[本文引用: 1]

Asad M W A .

Implementing a blending optimization model for short-range production planning of cement quarry operation

[J].Journal of Mining Science,2010465):525-535.

[本文引用: 1]

Souza M J F Coelho I M Ribas S et al .

A hybrid heuristic algorithm for the open-pit-mining operational planning problem

[J].European Journal of Operational Research,20102072):1041-1051.

[本文引用: 1]

Fioroni M M Franzese L A G Bianchi T J et al .

Concurrent simulation and optimization models for mining planning

[C]//Proceedings of the 2010 Winter Simulation Conference.New York:IEEE,2008759-767.

[本文引用: 1]

Stevanovic D Kolonja B Stankovic R et al .

Application of stochastic models for mine planning and coal quality control

[J].Thermal Science,2014184):1361-1372.

[本文引用: 1]

Shishvan M S Benndorf J .

The effect of geological uncertainty on achieving short-term targets:A quantitative approach using stochastic process simulation

[J].Journal of the Southern African Institute of Mining and Metallurgy,20151163):259-264.

[本文引用: 1]

黄启富陈建宏 .

基于PSO的矿山企业动态配矿优化研究

[J].计算机工程,2011378):175-177180.

[本文引用: 1]

Huang Qifu Chen Jianhong .

Research on dynamic mine ore blending optimization based on particle swarm optimization in mining enterprises

[J].Computer Engineering2011378):175-177180.

[本文引用: 1]

杨珊陈建宏杨海洋 .

基于Xpress-MP堆积型铝土矿堆场配矿优化研究

[J].金属矿山,2010393):9-11.

[本文引用: 1]

Yang Shan Chen Jianhong Yang Haiyang et al .

Optimization research of accumulated bauxite ore blending in yard based on Xpress-MP

[J].Metal Mine2010393):9-11.

[本文引用: 1]

王李管宋华强毕林 .

基于目标规划的露天矿多元素配矿优化

[J].东北大学学报(自然科学版),2017387):1031-1036.

[本文引用: 1]

Wang Liguan Song Huaqiang Bi Lin et al .

Optimization of open pit multielement ore blending based on goal programming

[J].Journal of Northeastern University(Natural Science)2017387):1031-1036.

[本文引用: 1]

吴丽春王李管彭平安 .

露天矿配矿优化方法研究

[J].矿冶工程,2012324):8-12.

[本文引用: 1]

Wu Lichun Wang Liguan Peng Ping’an et al .

Optimization methods for ore blending in open-pit mine

[J].Mining and Metallurgical Engineering2012324):8-12.

[本文引用: 1]

任助理毕林王李管 .

基于混合整数规划法的自然崩落法放矿计划优化

[J].工程科学学报,2017391):23-30.

[本文引用: 1]

Ren Zhuli Bi Lin Wang Liguan et al .

Optimization of drawing scheduling based on mixed integer programming in block cave mining

[J].Chinese Journal of Engineering2017391):23-30.

[本文引用: 1]

Rodrigo M Enrico Z Fredy K et al .

Availability-based simulation and optimization modeling framework for open-pit mine truck allocation under dynamic constraints

[J].International Journal of Mining Science and Technology,2013231):113-119.

[本文引用: 1]

Coelho V N Souza M J F Coelho I M et al .

Multi-objective approaches for the open-pit mining operational planning problem

[J].Electronic Notes in Discrete Mathematics,2012392):233-240.

[本文引用: 1]

Du D Z Frank H .

Computing In Euclidean Geometry

[M].EvanstonWorld Scientific1992.

[本文引用: 1]

Boissonnat J D Nielsen F , Nock R .

On bregman voronoi diagrams

[J].Discrete and Computational Geometry,2010442):281-307.

[本文引用: 1]

陈丽胡乃联 .

基于Voronoi图的爆区品位和矿量计算

[J].爆破,2011282):39-41.

[本文引用: 1]

Chen Li Hu Nailian .

Method of calculating ore grade and quantity of blast zone based on Voronoi diagram

[J].Blasting2011282):39-41.

[本文引用: 1]

杨桦胡乃联孙晓 .

基于炮孔化验数据的资源模型动态更新技术研究

[J].矿业研究与开发,2013332):105-109.

[本文引用: 1]

Yang Hua Hu Nailian Sun Xiao et al .

Dynamic update technology of geological resource model based on blasthole data

[J].Mining Research and Development2013332):105-109.

[本文引用: 1]

Ghasemi E Amnieh H B Bagherpour R .

Assessment of backbreak due to blasting operation in open pit mines:A case study

[J].Environmental Earth Sciences,2016757):552-558.

[本文引用: 1]

Bye A .

The strategic and tactical value of a 3D geotechnical model for mining optimization,Anglo Platinum,Sandsloot open pit

[J].The Journal of the South African Institute of Mining and Metallurgy,20061064):97-104.

[本文引用: 1]

/