, 顾晓薇, 胥孝川东北大学 资源与土木工程学院，辽宁 沈阳 110819
收稿日期：2017-05-12
基金项目：国家自然科学基金资助项目(51474049, 51674062)；国家自然科学基金青年基金资助项目(51604061)；辽宁省自然科学基金资助项目(2014020040)；中央高校基本科研业务费专项资金资助项目(N160104009)。
作者简介：王青(1962-)，男，内蒙古兴和人，东北大学教授，博士生导师;
顾晓薇(1971-)，女，辽宁凤城人，东北大学教授，博士生导师。
摘要：针对“剥岩高峰”对矿山生产带来的不利影响，研究了设备配置更新策略与露天开采外委剥离的优化问题.基于铲运设备的运营成本与生产能力随设备役龄变化的关系，建立了决策单台设备退役时间的动态经济寿命计算模型，并根据此模型给出了具体的设备配置算法步骤.设置一系列不同的矿山自剥离量与外委剥离量的分界水平，计算各分界水平条件下的外委单价可行区间，然后在区间内确定外委商谈价格并分析得出与之对应的总成本净现值最低的外委方案.实例计算证明，采用该方法减少设备投资的成本要高于外委岩石剥离成本.
关键词：剥岩高峰外委经济寿命设备配置更新单价可行区间
Optimization on Stripping Plan of Outsourced Open-Pit Mines
WANG Qing, ZHANG Zhi-tong
, GU Xiao-wei, XU Xiao-chuanSchool of Resources & Civil Engineering, Northeastern University, Shenyang 110819, China
Corresponding author: ZHANG Zhi-tong, E-mail: zhangzhitong8776 @163.com
Abstract: Since the "stripping peak" has a negative effect on mine production, the tactics of equipment update and the optimization on stripping plan of outsourced open-pit mines were studied. Based on the relation between the cost of stripping equipment, the production capacity and the changing of the service age, a dynamic economic life model was established to decide the retire time of a single equipment. Meanwhile, the configuration algorithm steps of a specific equipment was also given out. First, a series of different levels of mine stripping and outsourced stripping boundaries were set to calculate the feasible range of unit cost for the outsourced program, then the outsourced contract price was determined in the feasible range; Finally, based on the contract price, the net present value of the total cost was analyzed to get the optimized outsourced plan. It is proved that the cost of reducing equipment investment by this method is higher than the cost of outsourced rock stripping.
Key Words: stripping peakoutsourcingeconomic lifeequipment updatefeasible range of unit cost
露天开采是一个复杂的系统工程，影响开采经济效益的要素众多，其中，设计方案要素和主体设备的配置对露天开采的总体经济效益起决定性作用，也是露天开采优化领域的研究重点[1].露天矿生产过程中，生产剥采比随深度一般都经历一个由小到大, 达到高峰后再逐渐减小的变化过程.其中，岩石剥离高峰期会导致矿山短期内需要集中大量的设备并配备足够多的工作人员.然而，这些设备和人员在高峰过后又要削减.这样短期的增减会导致设备利用率降低、基建费用增大、生产成本增高, 并且使露天矿的生产组织工作复杂化.因此，为了得到较稳定的生产剥采比，在生产计划中常进行剥采比的均衡.然而，对于一定的矿体形态、最终境界和开采方式，剥采比均衡的结果往往是将剥离高峰期的岩石提前剥离.由于资金的时间价值，大量提前剥离会降低总体经济效益[2].
在20世纪80年代初期, 云庆夏教授首次提出了动态存储规划均衡生产剥采比的理论方法, 该方法考虑了剥岩量和剥岩费用随时间的动态变化及剥岩费用的时间价值, 以露天开采总剥离费用最小为原则均衡生产剥采比[3].但该方法需要获取的未知量，要通过传统的手工设计方法是很难实现的, 因此, 至今尚未在生产设计中得到实际应用[4].Alarie等就露天矿调度优化方法做了系统的综述[5]，但露天矿装运系统的运营效率仍有改进空间，优化方法的研究仍在继续[6].?elebi对于露天开采计划优化的问题，建立了设备的选择和评估模型[7].但是，该动态规划模型没有考虑剥离设备的运营成本与生产能力随设备役龄的变化问题，也没有考虑设备更新策略.Hartman等是当前设备更新问题研究的最前沿与集大成者，基本涉及到设备更新问题的各个研究角度[8-9].Sevim等应用动态规划研究了露天矿开采中外委(合同)剥离的优化问题，确定每年最佳合同剥离量以及矿山自己应拥有的剥离生产能力[10].本文结合某露天矿山的采剥计划，在剥离的高峰期雇佣公司外部的作业队承担部分岩石的剥离工作，并在剥离计划优化中纳入了设备的更新问题.
1 按动态经济寿命进行岩石采剥与运输设备配置计划的算法从设备的投资看，使用寿命越长，同样的投资可以完成的总作业量越多，那么分摊到单位作业量上的投资费用就越低.从设备运营成本看，使用寿命越长，年运营成本越高，而年作业量(生产效率)越低，导致分摊到单位作业量上的运营成本越高.单独考虑一台设备更新时，最好的标准是在设备的整个服役期间内单位产量的总费用现值最低，而现值的计算要将设备运行期内每一年的年生产能力、年运行成本和设备残值计算在内，同时要考虑成本的变化对计算带来的影响.考虑到实际情况，如果设备更新太早，服役年数太少，总产量低，单位产量分摊的投资费用很高，使单位产量的总费用也高；如果更新太晚，服役年数太长，总产量高且分摊到单位产量的投资低，但运行成本很高，单位产量的总费用还是高.因此，有一个最佳的更新时间，即最佳服役年数.
使用寿命期内分摊到单位作业量上的费用现值，定义为单位折现成本，它是使用寿命n的函数，记为c(n)；使单位折现成本c(n)达到最小时的使用寿命称为动态经济寿命ne.
 | (1) |
按照设备动态经济寿命进行更新配置，WK-10B型电铲负责铲矿作业，TR-100型卡车负责运输作业，两种设备各自进行更新配置.以电铲设备为例，计算步骤如下：
步骤1??输入已知的露天矿采剥计划，包括矿山需自己剥离的岩石量(Wi)和矿山开采寿命(N)，计算电铲、卡车设备的动态经济寿命ne.
步骤2??令i=1，首先根据矿山第1年采剥计划确定所需电铲设备的数目，满足：
 | (2) |
步骤3??令i=2，矿山第1年采剥计划确定所需电铲设备数目n1，设备的役龄相应为1年.满足：
 | (3) |
步骤4??令i=i+1，如果i < ne+1，执行步骤4，否则转入步骤5.
当矿山开采到第i年时，满足：
 | (4) |
步骤5??如果i=ne+1，即设备达到动态经济寿命，则该设备应该选择退役，令电铲设备的生产能力Pi(xi, k)=0.
步骤6??当i>ne+1时，这时候达到动态经济寿命的设备已经退役了，其他设备仍按照上面步骤4进行计算，直至i=N，即矿山开采结束年份.
步骤7??算法结束，输出该生产计划设备配置方案，计算出在矿山整个服务期内总的矿山剥岩和运输的成本现值.
2 矿山自剥离与外委剥离最优方案的确定方法2.1 优化原理不进行生产剥采比均衡，基建后每年的剥岩量变化如图 1所示.按照一定的变化步长，分别设置不同的自剥离与外委剥离界线值Q1，Q2，Q3.如果矿山自己剥离与外委剥离的分界线设置在Q1水平，那么Q1线以上的岩石量就是需要外委剥离的部分，外委开始于第t1年，结束于第t2年，每年的外委剥离量和外委剥离总量都可计算出来，留给矿山自己每年需要承担的剥离量也随之而定；如果这一分界线从Q1变为Q2，那么，每年的外委剥离量和外委剥离总量都相应地增加，外委的时间段增长；反之，每年的外委剥离量和外委剥离总量都相应减少，外委的时间段缩短.总的优化思路是，以一定的增量设置一系列不同水平的矿山自剥与外委剥离的分界水平，计算各分界水平条件下外委剥离部分的成本现值与矿山自剥离部分的运营成本现值之和，成本现值总额最低者即为最优方案.
图 1(Fig. 1)
 | 图 1 剥岩量随时间变化的曲线Fig.1 Changing of the stripping rock with time |
剥离高峰一般都不是出现在矿山开采的前几年中，如果设定的某个矿山自有剥离能力是一个等于分界水平的常数，仍然意味着大量的提前剥岩，考虑到资金的时间价值，这样做会降低矿山的整体经济效益[11]；如果只考虑矿山单位剥岩成本，就不能反映短期的剥离高峰带来的不经济性.这种不经济性主要体现在剥离设备数量的突然大量增加，高峰过后又有多台剥离设备在未到“退役年龄”时提前退役.因此，矿山的自有剥离能力必须转化为矿山剥离设备的配置，并考虑设备退役和更新问题，使所配置设备的生产能力尽量沿着对应于分界水平的自剥量变化，如图 1所示，对应于分界水平取为Q1值的自剥离量即为曲线abcd.
2.2 优化算法的步骤首先，按照矿山采剥计划，计算每年需要的岩石剥离量，设置矿山自有剥离与外委剥离分界水平初始值Q0及变化步长ΔQ，把矿山的自有剥离能力转化为矿山剥离设备的配置，并考虑设备退役和更新问题.
第一年：根据当年需要的岩石剥离量，为矿山配置设备，使其超前剥离量等于或大于0.
第二年：第一年的剥岩和运输设备每台的役龄增加1年，成为第2年初的生产设备，如果有达到动态经济寿命的设备，根据具体情况，更新或直接退役，目的是使前2年的实际剥离量减去前2年的需要剥离的岩石量大于或等于0.如果把所有到了动态经济寿命的设备都更新还不够，就得投入一台或多台新设备.
第三年：第二年的设备每台的役龄增加1年，与第二年相同的算法，得到第三年的设备队伍以及其累积超前剥离量.
如此每年进行下去，直到最后一年.
在上述计算过程中，为矿山自己的剥离量配置设备后，每年外委剥离岩石量及其发生的时间段都可以相应地确定出来.其中，必须保证任何一年的矿山自有剥离量与外委剥离量之和不低于当年的累计需求剥离量，否则，矿石就无法达到生产目标.计算这一方案的运营成本现值总额, 然后，根据之前确定的外委界线变化步长ΔQ，改变分界水平，重复上述步骤，得出最佳方案.本文应用的矿山自有剥离量与外委剥离量分界线的优化算法中，降低了矿山的基建投资成本，充分考虑到了资金的时间价值对矿山整体经济效益的影响；提高了矿山设备的使用效率，使矿山的总投资额大大减少.
3 实例应用以某露天矿设备运营数据为基准.电铲型号为WK-10B，购置价格I为1 885万元，技术寿命N为20年；卡车型号为TR100，购置价格I为800万元，技术寿命N为12年.折现率取6.8%，年均真实成本上升率取1.5%.结合式(1)计算单位折现成本随使用寿命n的变化函数c(n)，求得电铲和卡车的动态单位成本现值.其中，WK-10B型电铲的动态经济寿命为15年，TR100型卡车的动态经济寿命为8年.
根据当前技术经济参数，对该露天矿进行开采计划优化，优化结果如图 2所示[12].岩石剥离量的高峰出现在第4年到第12年.
图 2(Fig. 2)
 | 图 2 某露天矿开采计划优化结果Fig.2 Mining optimization result for the open-pit mine |
本算例按照动态经济寿命进行更新，卡车设备配置方案如表 1所示.表中“自剥离与外委剥离界线”表示矿山自行剥离岩石量的上限，例如将该界线设置为1 000万t时，矿山自行剥离岩石量的上限为1 000万t；该界线设为3 800万t时，表示矿山完全采取自行剥离的方案，即外委岩石量为零.
表 1(Table 1)
 表 1 卡车设备配置方案Table 1 Dispatching plan of transportation equipment
| 表 1 卡车设备配置方案 Table 1 Dispatching plan of transportation equipment | ||||||||||||||||||||||||||||||||||||
如表 1所示，“+”表示投入新设备；“-”表示退役设备；“±”表示退役设备的同时又投入了新设备；符号后面的数字为更新的台数.更新台数与投入设备的工作台数以及更新时间是一一对应的.例如表 1中，当“自剥离与外委剥离界线”取1 000万t时，更新台数“(-2, +1)”表示在开采的第11年年初，退役2台设备的同时，投入了1台新设备，此时的工作台数为3台.
在矿山自剥离与外委剥离最优方案的计算模型中，设置完成不同的自剥离与外委剥离界线的值后，可计算出各分界水平条件下，除去外委剥离岩石的部分，矿山完全靠自身剥岩的成本现值为Cz, i；矿山不采取外委方式，完全靠自身进行剥离工作时的剥岩成本现值为Cw, i.如式(5)所示，为矿山外委方案的最高外委单价，因为若承包商提出的外委价格高于该值，那么外委出对应方案中的岩石量所需的费用要高于矿山完全靠自身剥离时的运营成本，这样外委的方式就不能起到降低矿山建设投资的作用了.
 | (5) |
外委的最低单价应该为各分界水平条件下，Cz, i与对应的矿山靠自剥离的岩石量的比值.因为若承包商的定价小于该值，那么最经济的方式应该是把所有的剥岩工作都给承包商完成，这样研究就变得没有意义.给出不同的矿山自剥离与外委剥离界线条件下的外委单价可行区间，以及在可行区间内取不同外委商谈价格时，各方案的运营成本现值总额如表 2所示.
表 2(Table 2)
表 2 不同外委方案运营成本现值对比Table 2 Comparison of the cost present value with different outsourced projects
| 表 2 不同外委方案运营成本现值对比 Table 2 Comparison of the cost present value with different outsourced projects | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
在计算出的外委单价区间内，针对确定出的外委商谈价格，有与之对应的唯一的最佳外委剥离界线值.例如，当外委商谈价格定为3.5元/t时，自剥离与外委剥离界线设置在1 500万t水平时，是最佳的方案，该方案的运营成本现值比矿山完全靠自身剥离的方式的成本要节省约12.82%.由此可知，最佳的外委剥离界线值需要根据矿山与合同商共同商订的外委单价进行确定.
4 结论1) 根据矿山设备的动态经济寿命模型，分别计算出某矿山WK-10B型电铲设备的动态经济寿命是15年，TR-100型卡车的动态经济寿命是8年.
2) 结合某矿山的采剥计划, 给定矿山自行剥离量与外委剥离量分界线的变化步长，确定出一系列岩石自行剥离量与外委剥离量的分界水平，计算得出与之对应的各方案的矿山铲运设备的更新策略以及自有设备运营成本现值总额.
3) 实例计算结果表明，一方面，不同的岩石自行剥离量与外委剥离量的分界水平对应的外委价格可行区间不同，要根据各方案对应的可行价格区间来确定外委商谈价格；另一方面，外委商谈价格确定后(分别取外委商谈价格为2.5, 3.0, 3.5元/ t)，得出的总成本净现值最低的岩石外委与自行剥离分界线，即最优方案也随之确定(分别为1 000, 1 000, 1 500万t).
参考文献
| [1] | Newman A M, Rubio E, Caro R, et al. A review of operations research in mine planning[J].Interfaces, 2010, 40(3): 222–245.DOI:10.1287/inte.1090.0492 |
| [2] | 陈俊智. 露天矿生产剥采比的优化与控制[D]. 昆明: 昆明理工大学, 2001. ( Chen Jun-zhi. Optimization and regulation for stripping-ration of open-pits mines[D]. Kunming: Kunming University of Science and Technology, 2001.http://cdmd.cnki.com.cn/Article/CDMD-10674-2001002758.htm) |
| [3] | 云庆夏. 均衡露天矿剥离工作的新方法[J].有色金属(矿山部分), 1984(4): 15–18. ( Yun Qing-xia. A new method to balance open-pit mine stripping work[J].Nonferrous Metals(Mining Section), 1984(4): 15–18.) |
| [4] | 左瑞, 露天矿生产剥采比的优化调控[D]. 武汉: 武汉理工大学, 2011. ( Zuo Rui. Optimization and regulation of open-pit stripping ratio[D]. Wuhan: Wuhan University of Technology, 2011.http://cdmd.cnki.com.cn/Article/CDMD-10497-1011104975.htm) |
| [5] | Alarie S, Gamache M. Overview of solution strategies used in truck dispatching systems for open pit mines[J].International Journal of Surface Mining, Reclamation and Environment, 2002, 16(1): 59–76.DOI:10.1076/ijsm.16.1.59.3408 |
| [6] | Ercelebi S G, Bascetin A. Optimization of shovel-truck system for surface mining[J].Journal of the Southern African Institute of Mining and Metallurgy, 2009, 109(7): 433–439. |
| [7] | ?elebi N. An equipment selection and cost analysis system for open-pit coal mines[J].International Journal of Surface Mining Reclamation and Environment, 1998, 12(4): 181–187.DOI:10.1080/09208118908944042 |
| [8] | Hartman J C, Rogers J. Dynamic programming approaches for equipment replacement problems with continuous and discontinuous technological change[J].IMA Journal of Management Mathematics, 2006, 17(2): 26–30. |
| [9] | Rogers J L, Hartman J C. Equipment replacement under continuous and discontinuous technological change[J].IMA Journal of Management Mathematics, 2005, 16(1): 23–36.DOI:10.1093/imaman/dph027 |
| [10] | Sevim H, Wang Q, de Tomi G. Economics of contracting overburden removal[C]//Proceedings of the International Symposium APCOM. Berlin, 1990: 573-584. |
| [11] | Latorre E, Golosinski T S. Definition of economic limits taking into consideration of time value of money[J]. Canadian Institute of Mining, Metallurgy and Petroleum (CIM) Journal, 2011, 2(3): 162–170. |
| [12] | Gu X W, Wang Q. Dynamic phase-mining optimization in open-pit metal mines[J].Transactions of Nonferrous Metals Society of China, 2010, 20(10): 1974–1980.DOI:10.1016/S1003-6326(09)60404-0 |
