第38卷第6期2017年6月东北大学学报（自然科学版）Journal of Northeastern University(Natural Science)Vol.38 ,No.6Jun.20 17doi：10. 3969/j.issn. 1005 -3026. 2017.06.024基于混合整数规划法的地下矿采掘计划编制刘晓明，徐志强，陈鑫，钟德云(中南大学资源与安全工程学院，湖南长沙410083)摘要：矿山采掘计划编制对企业总体经济效益具有深远的影响，手动编制技术已难以满足矿山的需 求.为有效提髙计划编制的髙效性和合理性，通过对无底柱分段崩落法的研究分析，建立了以出矿品位波动最 小为目标函数的混合整数规划（MIP)模型.模型包含了生产开采的全过程，以保证井下各工序的合理衔接.对 北洺河铁矿采掘计划编制的试验结果表明，该方法具有编制时间短、工序衔接更合理、计划结果较手动方法更 优越的特点.研究表明采用该方法能有效提髙计划的合理性，使企业获得更大经济效益.关键词：地下矿山；无底柱分段崩落法;采掘计划编制;混合整数规划；品位波动最小中图分类号：TD 214 文献标志码：A文章编号：1005 -3026(2017)06 -0880 -06Production Scheduling in Underground Mining Using Mixed Integer ProgrammingLIU Xiao-ming，XU Zhi-qiang，CHEN Xin，ZHONG De-jun(1.School of Resources and Safety Engineering，Central South University，Changsha 410083，China. Corresponding author:XU Zhi-qiang，E-mail: 1184609635@)A b stra ct：Production scheduling has a profound influence on the overall economic benefits o f the m ining enterprise,and the manual technology has d ifficu lty to meet the needs o f mine.To improve the efficiency and rationality o f scheduling,we studied the non-pillar sublevel caving m ethod，and presented a mixed integer programming (M IP)m odel.This model makes m inimum fluctuations o f ore grade as the objective fu nctio n，and it contains the whole process o f m ining to ensure all procedure could be more reasonable connected.In the test o f the production scheduling o f the Beiminghe iron m in e，the results show that it spent shorter time and the procedure was more reasonable，what's m ore，the results are better than the manual method.The study shows that this method can effectively improve the rationality o f schedule and make the enterprise get more economic benefits.K ey w o rd s：underground m in e；non-pillar sublevel caving m ethod;production scheduling; mixed integer programming；m inim um fluctuations o f ore grade地下矿采掘计划编制是矿山生产经营管理中 一项极其重要的内容，计划编制是否合理，将直接 影响到矿山企业的总体经济效益[|-2].传统的手 动编制方法不仅耗时长、强度大,而且准确性差，修改难度大[2].究其原因，主要是在编制计划过 程中,需要综合考虑各生产工序和采场之间的时 空间制约性及其连续性，同时还应尽可能使企业 经济效益最大化.随着计算机三维可视化技术的 迅速发展，地下矿三维可视化技术给矿山采掘计划编制提供了一个良好的环境.关于自动编制矿山采掘计划，国内外学者运 用多种方法对其进行了研究[2-9]，大体可分为三 类:计算机模拟技术、人工智能和数学规划法.文 献[1]应用计算机模拟技术对冬瓜山铜矿未来两 年的回采过程进行了模拟分析，不足是需要有预 演的计划方案.文献[2]运用人工智能中的遗传 算法对矿山采掘计划进行了研究分析，但方法存 在初始群体选取困难和可能因迭代次数不够而错收稿日期：2015-12-28基金项目：金属矿山安全与健康国家重点实验室开放基金资助项目（2015- JSKSSYS -01).作者简介：刘晓明（1982 -),男，江西新余人，中南大学讲师，博士.第6期刘晓明等：基于混合整数规划法的地下矿采掘计划编制881失最优方案的风险.数学规划法通过建立数学模 型来求解目标函数进而得到矿山采掘计划，如 D6n al等[9]提出了以贴现金属量最大为目标函数 的回采顺序优化模型;Kumral[4,6]基于启发式和 混合整数规划法提出了以净现值最大为目标函数 的数学模型;Nehring[3]提出了新的优化矿山短期 生产调度和设备布置模型.这些模型的求解结果 均是在井下其他工序能很好相互衔接的基础上得 到的，然而矿山的生产开采是一个系统性工程，井 下采掘计划有其固有的特殊性和复杂性，单一考虑 矿块回采顺序最优，难以保证井下其他工序的进度 能跟得上矿块的回采顺序，导致工序之间无法衔接 而影响进度计划的实施.因此，一个完备的采掘计 划应当包含矿山生产过程的各工序.本文综合分析了无底柱分段崩落法的生产开 采，提出了一个包含采掘生产全过程的混合整数 规划模型，以每月出矿品位波动最小为目标函数. 将该模型运用于北洺河铁矿2016〜2018年采掘 计划的编制.1计划编制混合整数规划模型1.1计划编制原理矿山在回采井下矿石之前，需要先完成开拓、探矿、采准和切割等准备工作.计划的目的就是将 这些工作在时间上按一定的先后顺序进行排列，并尽可能使矿山企业获得最大效益.而以无底柱 分段崩落法开采的矿山井下生产开采可划分为掘 进、中深孔和回采三个阶段,掘进阶段又可分为系 统工程掘进和采准掘进.系统工程掘进主要指开 拓巷道、探矿巷道以及溜井的掘进;采准掘进是指 采场进路、切割巷和切割井的掘进;中深孔阶段是 在采场内钻凿中深孔；回采阶段包括落矿和矿石 运搬等回采作业.因此可根据各阶段开采特点和 工序的回采顺序构建出混合整数规划计划模型. 1.2混合整数规划模型1.2.1目标函数在计划模型中，选取计划期内每月出矿品位 波动最小为目标函数以确保出矿品位均匀，进而 降低配矿和选矿成本，提升企业经济效益.m in Z +<?,2) .(1)teT式中：r为时间f的集合;定义心和分别为f时 期出矿品位高于和低于入选品位的偏差量，且均 取值为正.1.2.2决策变量通过对无底柱分段崩落法开采过程的划分，选 取回采、中深孔、采准掘进和系统工程掘进为计划 模型的二元决策变量，并定义当一个采场或一段系 统工程在t时期开始推进时,取值为1，否则为0.H a t{1，a采场在t时期开始回采；= 0，其他.(2) Z. = L，a采场在t时期开始中深孔；，其他.(3)J a t = {0，a采场在t时期开始采准掘进；，其他.(4)Z(t = {0，系统工程/在t时期开始推进；，其他.(5) 1.2.3约束条件回采约束：1)逻辑约束:计划期内a采场只开采一次.X H a t ^ 1，V a e A，t e T.(6)teT式中:A为采场a的集合.2) 变量非负，二元变量.Hat=1 或 0，V a e A，t e T .(7)3) 每月出矿品位约束.X X8a H a,'J/p + 8a- 8n=8-a e A t'eTt—h a+1^t^ tV a e A，t e T.(8)式中：h a为a采场回采持续时间;8为矿石的人选 品位;8a为a采场矿石平均品位;Y为a采场每月 采下的矿石量;P为矿山每月采下矿石总量.4) 分段回采设备数量约束.X X Hat'矣U v e V，t e T.(9)a eAv t'eTt—h a+1备t备t式中：y为分段v的集合;Av为分段v内采场a的集合;乂,为分段v在t时期最多能容纳的回采设备数量.5)回采设备总量约束.X X H at'=N h，V e T，(10)ae A t e Tt—h a+1备t备t式中:N h为回采设备总数.6) 水平约束.Ha t 矣X Ha't'，Va e A，a'e A al，t e T.t eTt'忘t—h a'(11)式中:A a,为水平约束矩阵，规定当采场a'必须优 先于采场a回采时取值为1，否则为0.882东北大学学报（自然科学版)第38卷7)垂直约束.K X H d,t,,Vfl e A,a"e A a v,t e T. (12)t'eTt'^t-0.5hd"式中:AaV S垂直约束矩阵,并规定当采场a"必须优先于采场a回采时取值为1,否则为0.8)回采设备调度约束.22X H a t^X X X I，v'^va eA v't'eT v'aeAv't"eTt—ha+1备t备t t—h a备t"备t—1V t e T,v e V,(13)中深孔约束：1)逻辑约束:计划周期内a采场中深孔只钻 凿一•次.X L 矣 1, V a e A,f e T.(14)te T2)变量非负，二元变量.Zat =1 或 0, V a e A,f e T.(15)3) 采场内中深孔超前回采约束.H t(X Z t,V a e A，t e T.(16)t'eTt'^t—Z a式中:Z a为a采场钻进中深孔的时长.4)中深孔设备数量约束.X X Z a t,=N z,V t e T.(17)ae A t e Tt—Za+1矣t'矣t式中:N z为中深孔设备总数.5) 中深孔设备调度约束.X X X Zat'<v'忘v a e A v't'eTt—Za+1矣t'矣tX X X Zat,,V t e T,v e V.(18)v'忘v aeAv't"e Tt—Z a矣t"矣t—1采准掘进约束：1)逻辑约束:计划期内a采场采准掘进只掘 进一次.X人t 矣 1, V a e A,f e T.(19)te T2)变量非负，二元变量.人t =1 或 0, V a e A,f e T.(20)3) 采准掘进超前中深孔约束.Zat 矣X 人t',V a e A,t e 厂（21)t'eTt'^t—./a式中:九为a采场进行采准掘进的时长.4) 采准掘进设备数量约束.X X人t'=N.,V t e 厂（22)ae A t e Tt—Z a+1^t'^t式中:N.为采准掘进设备总数.5) 采准掘进设备调度约束.X X X人,'<v'忘v a e A v't'e Tt—.Z a+1^t'^t X X X人t",V t e T,v e V.(23) v'忘v aeA v't"e Tt—./a^t，^t—1系统工程约束：1) 逻辑约束:计划期内系统工程/只掘进一次.X^；t ^1,V f e r,/e L.(24)te T式中:L为系统工程/的集合.2) 变量非负，二元变量.X/t=1 或0, V r e T J e L .(25)3) 系统工程超前采准掘进约束.人t 矣X X/t',V a e A,t e T,/ e 尤/.(26)t eTt备t—文/式中:x,为掘进系统工程/的时长;X a,为a采场的 采准与系统工程/之间的约束矩阵.当系统工程/优先于a采场的采准时取值为1,否则为0.4)系统工程间相互约束.&矣X Z r t',V t e T,/e L,/'e X u . (27)t eTt备t—文/式中:X…为系统工程之间的约束矩阵.规定当系 统工程/'优先于/掘进时取值为1,否则为0.5) 系统工程设备数量约束.X X X W=N,,V t e T.(28)/eL t'eTt—文/+1备t备t式中:K为系统工程的设备总数.2模型分析2. 1约束条件分析在创建的计划模型中，选取井下各生产工序 在时空间上存在的相互关系为约束条件.式（8)是用于衡量t时期内出矿品位与入选品位之间的 偏差量，所得值心，^2将为目标函数的求解提供 基础数据，其余约束条件可分为以下五类：1)逻辑约束:式(6) ,(14) ,(19) ,(24)限的是每个采场和每段系统工程在计划期内不能重 复开采，而式（7) ,(15) ,(20) ,(25)则是保证当 其参与计划编制时，取值为1,否则为0.式（27 )是系统工程之间的约束关系，其定义如图1所示，即按照推进方向只有当1号系统工程完成后才可 以进行2号和3号的掘进，也只有当3号系统工 程完成后才可以进行5号的掘进，该约束目的是 防止编制的计划出现逻辑上的错误.2)设备数量约束:式(9)和（10)都是对回设备数量的约束，不同之处在于式（9)限制的是 同一分段内回采设备的数量，以防止回采设备在 同一分段内过度集中，而式（10)限制的则是整个第6期刘晓明等：基于混合整数规划法的地下矿采掘计划编制883矿山可用于回采作业的设备数量.同理，式（17), (22),(28)分别限制的是矿山用于推进中深孔、 采准掘进以及系统工程掘进的设备数量.图1系统工程之间约束关系Fig . 1 Constraints relationship between systemsengineering 3)空间约束：由无底柱分段崩落法开采特点 可知，回采中为降低矿石的损失贫化以及保证生 产安全，开采采场在水平和垂直方向上均会受到 一定的约束.因此模型中分别用式（11)和式（12) 对水平和垂直方向进行限制（图2),并定义水平 方向上限制采场a 开采的采场c 和^未完成时, 采场a 不能开采;垂直方向上限制采场^开采的 采场a 未完成其50%时，采场6不能开采./采场a \/■ \采场6图2采场空间约束关系Fig . 2 Space constraints of stopes 4)时间约束:井下各工序的开采在时间上都 具有一定的逻辑性，即一个采场在开始回采时，中 深孔已完成;一个采场在开始钻凿中深孔时，采准 掘进已完成;一个采场在开始采准掘进时，相关的系统工程已完成.因此模型中利用式（16)，(21)， (26)来保证各工序在时间上的连续性.5)设备调度约束：由于无底柱分段崩落法开 采的矿山其分段的开采顺序为下行式，因此为避 免设备出现向上调度的问题，模型中分别用式 (13)，(18)，(23)来限制回采、中深孔、采准掘进设备 的调度方向,其原理是:对于第一个分段:f 时期生产 采场数之和不大于f -1时期生产采场数之和;对于 第一、二个分段:f 时期生产采场数总和不大于f -1 时期生产采场数总和;对于第一、二、三个分段:f 时 期生产采场数总和不大于f -1时期生产采场数总和， 并依此类推直至所有计划分段.通过上述约束条件，能够很好地描述出生产 开采过程中井下各工序在时空间上的相互关系，也能够有效保证工序间的合理衔接，因此模型的 约束条件是可行的、合理的和完备的.2.2回采顺序优化路线通过对无底柱分段崩落法采掘生产过程及计 划编制理论的研究分析，采场回采顺序优化的技 术路线如图3所示.( 开始 ）_______________I______________1 .~~ t .I 采矿方法分析~1 丨生n r a 编制原厕~■ . t . 1~「 提取%学关系~~II 建立混合整k 规划模型I I 基础数据准备I~~ T '|在M A T L A B +Y A L M IP 職下编写観语言|? '调用C PLEX 求解模型♦I 回采顺序优化方案确定I( 结’束 )图3回采顺序优化路线Fig . 3 Optimization route of extraction sequence3工程实例3. 1 工程概况北洺河铁矿是五矿邯郸矿业有限公司的一座 大中型地下矿山.矿床产于燕山期闪长岩与奥陶 纪石灰岩接触带，为接触交代型磁铁矿床,矿物以 磁铁矿、黄铁矿为主，伴有少量赤铁矿、褐铁矿、 硫、钴等.矿床埋藏深度为265〜679 m ，铁矿石平均品位为49.79%，入选品位为40%，设计年产量为180万t ，采用无底柱分段崩落法开采.矿山在2015年末回采、中深孔、掘进的生产 现状为：回采和中深孔主要集中于-140和 -155 m 分段;采准掘进主要集中于-140，- 155， -170 m 分段;系统工程除-170 m 分段剩余少量外，主要集中于-185 m 分段系统工程的推进.矿山2016年采掘计划目标为217万t 铁矿石，其中主采场为181. 2万t ，外委矿为14. 4万t ，掘进带 矿为21. 4万t .井下各工序设备数量及其生产能 力如表1所示.表1设备及其生产能力Table 1 Equipments and production capacity工序名称设备数量/台设备型号月产能备注回采6T 0R 0400E 3. 1 万 ti 台备用中深孔5SIMBA 8 000 m 1台备用采准掘进5YT 23 型113 m1 台备用系统工程27655 型135 m-采场c采场^884东北大学学报（自然科学版)第38卷3. 2数据准备现利用创建的计划模型编制矿山2016〜2018年的采掘计划，预计参与回采的采场有-140, -155, - 170, - 185, -200 m分段，共 56个采场，计划周期为36个月，S P4 = [1,2, —，56]式5,6,…,n1.2, …,111.2, …,131.2, …,14_1，2, (11)T[1,2,…，36].a采场矿石的平均品位心=[0. 44,0. 47, …,0.35]，矿石入选品位g =0.40.系统工程以一 个周期内系统工程的推进长度进行量化，可得 乙=[1,2,…，46]，x, = [1,1，…,1]，而中深孔和采 准掘进均以采场为单元进行量化,所以采场各工 序剩余工作量的持续时间心=[1,8,…，4]，Za = [0,0,…,3]，九=[0,0,…,4].回采过程中一个采 场只配备1台回采设备，因此采场每月采下的矿 石量为1台回采设备的生产能力，即Y=3. 1万t，矿山每月采下矿石总量等于所有回采设备生产能力的总和,尸=15.5万t.分段v在f时期最多能 容纳的回采设备数乂, =3,由表1可知，每月用于 回采、中深孔、采准掘进的设备数量分别为5台、4 台、4台，即= 5，% = 4，% =4,用于掘进系统工 程的设备在矿山实际运用中有时为1台,有时为2 台,因此有1矣M矣2.根据水平约束、垂直约束、采场 与系统工程间约束以及系统工程之间相互约束的定 义并查分段平面图可知：a1 Aala56,1X al = i_^56,1a1，56，Aa v ='^1,1•"^1,56 " a56,56__^56,1"* 056,56-S1,46，x u=_屮1,1…屮1,46 _ ^56,46_-屮46,1"* ^46,46_其中:a^，，<p取值为1或0.3. 3模型求解在MATLAB + YALMIP环境中编写完模型 语言后，调用CPLEX求解器求解.计算机为: Windows7操作系统，Intel(R)i5处理器,4 G B内存.软件版本为：皿八丁[八32014，匸戸[丑乂12.5.1，YALMIP2013.模型求解时长8 095 s，变量个数为 7 704,约束个数为21 306,月最大品位偏差量为 0.035,目标函数值为0.267.采场回采顺序H a，如 表2所示.表2采场回采顺序求解结果Table2 Results of stopes extraction sequence H at采场a时间/月123 4 •■33343536 -140 m5#10000000-140 m6#10000000-140 m7#00000000-200 m9#00000000-200m10#00000000- 200m11#00000000根据Ha,的求解结果以及各采场的回采持续 时长,可编制出北洺河铁矿2016〜2018年采场回 采网络衔接甘特图，如图4所示.图4中，红、绿、黄、蓝、紫色分别代表5台回 采设备的调度情况,可知首采地段分别为-140 m 分段的5号、6号、8号采场以及-155 m分段的1号、4号采场，设备调度方向均为自上而下，且在 任何时间f内都有5台设备在回采作业,计划周 期内没有出现设备闲置的问题.3.4结果对比分析设备参数和生产能力将混合整数规划模型求 得的每月出矿品位偏差量心和g,2值与手动编制 的品位偏差量进行对比,并用折线图表示出来，其 结果如图5所示.再根据每月出矿品位，便可求得 北洺河铁矿计划期内的年回收金属量（图6).1 246 8 10 1213 TT ¥¥19--w—1T 1T1T201碎W月0^03月|〇4月|〇5月|〇6月|〇7月|〇8月|明月|l0月|l l月|l2月2017年01月丨02月丨03月|〇4月丨05月丨06月丨07月丨08月|〇9月丨10月丨11月丨12月2018年01月|一03月|〇4j j| 05_— 07月丨08月丨09月丨10月丨11小2月■ T O R0400E-1■ T O R0400E-2口 T O R0400E-3□ T O R0400E-4■ T O R0400E-5图4回采网络衔接甘特图Fig.4 Gantt chart of stopes extraction sequence注：图中只展示了计划期内参与开采的采场，未参与开采的采场未加以显示.第6期刘晓明等：基于混合整数规划法的地下矿采掘计划编制885由图5可知，当采用手动编制进度计划时，出矿品位整体波动很大，在前两年内，出矿品位均高于人选品位，且上向偏差在2%左右，到第三年时,全年采出矿石的品位又整体低于入选品位，下向偏差也基本维持在1% .而利用混合整数规划模型编制的采掘计划，前期7个月品位偏差较大,后期基本稳定在入选品位40%上下小范围波动，出现品位偏差较大的原因主要与2015年底备采矿量有关，因备采矿量品位固定，且需优先回采后才能回采新生成的备采采场，另外从前期回采持续时间来看，7个月时长能很好地与铁矿山备采储量保有期相呼应.图6为两种不同方法在回采计划期内金属量的回收情况，可知两种方法均能满足矿山年产金属量的要求.2016年2017年2018年三年回采总和图6回收金属量Fig.6 Quantity of metal recovered综合图5和图6的计算结果,可分别计算出两种方法在品位和金属量上偏差量的期望与标准差，结果如表3所示.通过数据对比可知，利用混合整数规划模型编制的采掘计划出矿品位更均衡，计划编制时间更短，因此利用混合整数规划模型编制的采掘计划较手动编制计划更优.表3偏差量的期望与标准差Table3 Expectation and standard deviation ofdeviations%编制方法品位金属量期望 标准差期望标准差手动编制 1.09 1.6116.9624.88 M IP编制0.42 1.08 6.5416.764结 论1)提出了一种编制无底柱分段崩落法采计划的混合整数规划模型.在总结前人研究成果的基础上,综合考虑生产开采的全过程，巧妙地将工序间存在的相互关系融入到模型中，最后通过求解模型得到矿山的采掘计划.2)模型充分考虑了生产开采的全过程，有避免了计划在实施过程中工序之间难以合理衔接的问题，使编制出来的生产计划可执行性更强.3)将模型运用于北洺河铁矿采掘计划编制的试验结果表明，利用混合整数规划模型编制矿山的采掘计划比手动编制计划月出矿品位和金属量波动小，且计划编制时间短、强度小、准确度高，因此编制结果更能满足矿山的生产需求.参考文献：[1]贾明涛，潘长良,谢学斌.深井矿山回采过程模拟及调控研究[J].湘潭矿业学院学报,2003,18(3) :9 - 12.(Jia Ming-tao，Pan Chang-liang，Xie Xue-bin. The simulationand control of the exploitation processes in a deeper deposit[J ]. Journal of Xiangtan Mining Institute，2003 ,18 ( 3 )：9-12.)[2]云庆夏，高文炜，陈永锋，等.优化地下采掘进度计划的进化算法[J].金属矿山，2002,31 (4) :11 -13.(Yun Qing-xia，Gao Wen-wei，Chen Yong-feng，et al.Evolutionary algorithm for optimizing the underground miningscheduling[ J]. Metal M ine，2002，31(4) ：11 -13.)[3 ] Nehring E T P K. Dynamic short term production schedulingand machine allocation in underground mining usingmathematical programming [ J ]. Mining Technology ,2010,119(4) ：212 -220.[ 4] Kumral M. Production planning of mines：optimisation of block sequencing and destination [ J ]. International Journal ofMining，Reclamation and Environment，2012，26 ( 2 ):93 - 103.[5] Yashar P，Hooman A，Dwayne D T. A multi-step approach forblock-cave production scheduling optimization [ J ] .International Journal o f Mining Science and Technology，2013,23(5) ：739 -750.[ 6] Kumral M. Multi-period mine planning with multi-process routes [ J ]. International Journal o f Mining Science andTechnology，2013,23(3) ：317 -321.[7] Topal M. Implementing a production schedule at LKAB' sKiruna mine [ J ]. Interfaces，2004，34(2) ：124 -134.[8] Martinez M A，Newman A M. A solution approach foroptimizing long- and short-term production scheduling atLKAB ’ s Kiruna mine [ J ]. European Journal of OperationalResearch ,2011，211(1) ：184 - 197.[9] Donal 0，Sullivan A N. Optimization-based heuristics forunderground mine scheduling [ J ] . European Journal ofOperational Research ,2015，241(3) ：248-259.。
