资源描述:
收稿日期 2016 -10 -17 基金项目 国家自然科学基金青年基金资助项目51604061ꎻ 国家自然科学基金资助项目51474049ꎬ 51674062ꎻ 辽宁省自然科 学基金资助项目2014020040ꎻ 教育部高等学校博士学科点专项科研基金资助项目20130042110012ꎻ 中央高校基本 科研业务费专项资金资助项目N160104009. 作者简介 胥孝川1986 - ꎬ男ꎬ四川绵阳人ꎬ东北大学博士后研究人员ꎻ 顾晓薇1971 - ꎬ女ꎬ辽宁凤城人ꎬ东北大学教授ꎬ博士生 导师. 第39卷第4期 2018 年 4 月 东北 大 学 学 报 自 然 科 学 版 Journal of Northeastern UniversityNatural Science Vol. 39ꎬNo. 4 Apr.2 0 1 8 doi 10.12068/ j. issn.1005 -3026.2018.04.023 多因素对露天矿最终境界设计的影响 胥孝川1ꎬ 顾晓薇1ꎬ 王 青1ꎬ 刘剑平1ꎬ2 1 东北大学 资源与土木工程学院ꎬ 辽宁 沈阳 110819ꎻ 2 东煤沈阳基础工程公司ꎬ 辽宁 沈阳 110016 摘 要 以技术经济参数及矿体赋存条件为依据ꎬ基于锥体排除法生成利润值最大的无约束境界ꎻ根据 地表构筑物在矿区的地理位置ꎬ建立境界优化范围约束线ꎬ在约束范围内优化最终境界并计算境界外矿石采 用地下开采的利润值ꎻ建立矿山开采生态成本计算模型ꎬ分别计算两种优化方案的生态成本. 研究结果表明 当只考虑经济因素时ꎬ如果地表构筑物的移动成本低于 0 99 亿元净现值ꎬ矿山适合采用无约束境界开采ꎻ 当综合考虑生态成本时ꎬ应该选择在受约束境界内开采. 关 键 词 露天矿ꎻ最终境界ꎻ生态成本ꎻ利润ꎻ地表约束 中图分类号 TD 8 文献标志码 A 文章编号 1005 -3026201804 -0570 -05 Influence of Multi ̄factors on Ultimate Pit Design of Open Pit Mine XU Xiao ̄chuan1ꎬ GU Xiao ̄wei1ꎬ WANG Qing 1ꎬ LIU Jian ̄ping1ꎬ2 1 School of Resources & Civil Engineeringꎬ Northeastern Universityꎬ Shenyang 110819ꎬ Chinaꎻ 2 Dongmei Foundation Company of Shenyangꎬ Shenyang 110016ꎬ China. Corresponding author GU Xiao ̄weiꎬ E ̄mail 493212755@ qq. com Abstract Based on the ecological ̄economic parameters and occurrence of orebodyꎬ the cone exclusive method was used to design the profit ̄maximum pit without constraint. According to the location of surface structures at the ore districtꎬ constraint lines were built to restrain optimization regionꎬ and the profit was obtained from ores outside constraint lines by underground mining. An ecological costs model of mine was built to calculate the ecological cost of two optimization schemes. The results showed that under only considering economic conditionꎬ it’s suitable for a unrestrained pit to mine if the moving cost of surface structures is less than 0 99 108 NPV. The restrained pit is the best choice while ecological costs are considered. Key words open pit mineꎻ ultimate pitꎻ ecological costꎻ profitꎻ surface constrains 矿产资源开采作为一个高风险行业ꎬ受到市 场行情及地质条件等不确定性因素影响ꎬ其中对 露天矿境界的确定影响最大. 露天矿最终境界设 计不仅涉及到矿区开采范围ꎬ同时还关系到矿岩 总量、生产能力大小、需要准备的尾矿库面积及排 土场面积等ꎬ这些都直接关系到矿山最终盈利大 小. 对露天矿境界优化的影响因素ꎬ国内外专家从 不同方面进行了研究. 文献[1 -2]考虑品位分布 的不确定性ꎬ对露天矿最终境界进行设计ꎻ文献 [3]分析了金属价格的变化特性ꎬ基于风险成本 最小化策略建立境界优化算法ꎻ文献[4]从资金 的时间价值考虑境界优化问题ꎻ文献[5]综合考 虑不同地质条件对境界的影响ꎬ并用神经网络算 法求解最优境界ꎻ文献[6]考虑不同方位帮坡角 的变化及复杂边坡角情况ꎬ采用浮动圆锥法优化 境界ꎻ文献[7]将采区、岩性及矿石品位等因素变 化而引起的生产成本及选矿回收率的变化反映在 境界优化中. 但在境界优化过程中ꎬ将环境成本作 为影响境界确定的因素的研究还相对较少. 虽然 目前国内外对矿山环境问题高度重视ꎬ如文献 [8]中提到加拿大颁布了采矿与矿物可持续开 发指南与探矿环保指导ꎻ文献[9]将生命周期评 价法用于矿山开采的整个生产流程ꎬ但很少有将 环境成本用于指导矿山设计. 鉴于此ꎬ本文以司家 营露天矿为例ꎬ分别从经济角度和生态环境角度 评价受地质条件约束下的最终境界选择. 1 境界优化算法 露天矿境界优化方法很多ꎬ包括浮锥移动法、 图论法、线性规划法及动态规划法等ꎬ本文使用的 锥体排除法是正锥排除法ꎬ即锥体顶点向上ꎬ排除 的是不盈利部分ꎬ改善了锥体排除法中锥体重叠 的影响ꎬ设置了多级优化. 其基本思想就是从模型 底部开始逐步排除那些满足一定几何约束条件且 不盈利的部分ꎬ剩下部分构成了最终境界. 根据矿 山地质资料ꎬ矿区范围及受高压线等地表构筑物 约束所圈定的开采范围如图 1 所示ꎬ为 -30 m 水 平的品位模型ꎬ图中空白部分表示目前该水平已 开采ꎬ空白范围内的深色部分代表矿石已采 出ꎬ根据剖面线 I - I′得到的剖面如图 2 所示简 化后的规则块体矿床模型. 图 1 -30 m 品位模型 Fig 1 Grade model on -30 m level 图 2 锥体排除法示意图 Fig 2 Sketch map of cone exclusive method 在图 2 中ꎬ首先从最低层 L =1 开始ꎬ构造 一个足够大的锥体ꎬ锥体的锥面倾角在不同方位 等于矿山设计的最终帮坡角ꎬ将锥顶点置于模块 的中心ꎬ判断是否有地表约束线. 如果有ꎬ则要判 断该模块是否位于图 2 中虚线圈定的范围内ꎬ如 果在虚线内ꎬ计算锥体内矿石量和岩石量. 根据采 矿、剥岩、选矿单位成本ꎬ精矿售价ꎬ采、选回收率 及精矿品位等判断采出锥体内的矿石和岩石是否 盈利ꎬ如果不盈利ꎬ则将锥体内的矿石和岩石都排 除掉即将图 2 中受该锥体影响的所有模块的底 部标高提升到锥面位置ꎻ如果盈利小于零ꎬ不作 任何处理ꎻ如果模块不在虚线内ꎬ则不用进行计算 和判断ꎬ直接删除ꎬ然后将锥顶点沿同一水平移动 到下一个模块中心ꎬ直到遍历完该水平所有的模 块ꎬ再将锥体向上移动一个水平根据需要ꎬ可以 是一个台阶高度ꎬ也可以是半个台阶高度或者1/4 台阶高度ꎬ判断这一水平的所有模块ꎻ当锥体从 最低水平移动到最高水平地表后ꎬ就完成了一 轮扫描ꎬ如果在此轮排除过程中ꎬ有锥体被排除ꎬ 则再进行下一轮扫描ꎬ否则算法结束ꎬ剩下模块组 成的境界即为最终境界. 2 生态成本计算模型 矿山生态成本就是对矿山开采造成的环境问 题的具体量化. 不同的采矿设计必然引起不同的 环境问题ꎬ比如最终境界圈定的范围大小将会直 接影响到矿区地表植被破坏的面积ꎬ剥离的岩石 越多ꎬ那么需要的排土场面积越大ꎬ产生的环境影 响也就越深远. 考虑到矿山的实际情况及数据的 可获得性ꎬ这里将矿山生态成本分为直接经济损 失、复垦成本、能源消耗的生态成本及外部生态价 值损失. 2 1 直接经济损失 对于以生物质生产为目的的土地ꎬ其直接生 态价值是能够提供的生物产品粮食、肉类、奶 类、林产品等在特定市场条件下的价值. 除市场 价格外ꎬ这一经济价值取决于土地的生物质生产 力生态系统的主要功能之一. 矿山对土地的直 接占压和挖损造成了土地直接经济价值的损失ꎬ 因此把它纳入生态成本是完全合理与必要的. 地 价是土地能够在未来为拥有者或使用者带来的直 接经济价值在当时当地市场条件下的一种体现 Cz= AV - C .1 式中Cz为直接生态价值ꎬ元hm -2ꎻA 为植被破 坏面 积ꎬ hm2ꎻ V 为 单 位 面 积 的 年 总 收 益ꎬ 元 hm -2 a -1ꎻ C 为 单 位 面 积 的 年 总 费 175第 4 期 胥孝川等 多因素对露天矿最终境界设计的影响 用ꎬ元 hm -2 a -1. 2 2 复垦成本 复垦成本是把矿山生产所破坏的各类土地ꎬ 按复垦方案中确定的复垦方向恢复到具有破坏前 或期望的生态生产力所需要的费用ꎬ包括复垦工 程结束后一个时期的养护费用 Cr= A cr.2 式中Cr为复垦和养护总成本ꎬ104元ꎻcr为单位 面积复垦和养护的费用ꎬ104元 hm -2. 2 3 能源消耗的碳排放成本 矿山生产中化石能源的消耗及爆破活动都会 向大气排放 SO2ꎬCO2及 NOx等ꎬ其中 CO2是导 致温室效应的主要气体成分. 处理 CO2排放有两 个途径①根据林地具有吸收 CO2的特性ꎬ通过 种植当量的林地以吸收矿山每年排放的 CO2ꎻ② 碳捕 捉 和 贮 存 carbon capture and storageꎬ CCS. 一些国家正在进行火电厂 CCS 的实验ꎬ国 际能源署International Energy Agencyꎬ IEA 也 在对 CCS 的成本和效率进行评估. 所以ꎬ能耗的 生态成本可依据能耗的碳排放量和碳捕捉与贮存 成本来计算. 碳捕捉法只需要知道由于能源消耗产生的 CO2ꎬ然后使用相关技术进行处理 CE= QN+ Qf+ Qe ce.3 式中CE为露天煤矿能耗生态成本总现值ꎬ104 元ꎻce为碳捕捉成本ꎬ元 t -1ꎻQ N 为爆破作业产生 的 CO2排放量ꎬ 104t ꎻQf为化石能源消耗的 CO2 排放量ꎬ 104tꎻQe为电消耗的 CO2排放量ꎬ104t. 2 4 外部生态价值损失 外部生态价值是土地可提供的除了带来直接 生态价值的生态服务之外的其他生态服务的价 值ꎬ这些生态服务没有直接进入经济系统或矿山 企业的财务核算体系ꎬ故称之为“外部生态价 值”. 不同类型土地的生态服务功能不同ꎬ其外部 生态价值也不同. 土地及其上的植被被破坏后ꎬ可 以认为生态系统失去了提供生态服务的功能ꎬ也 就失去了这部分价值ꎬ所以它是生态成本的重要 组成部分. 不同生态系统的生态功能不一样ꎬ以林 地为例ꎬ具有吸收二氧化碳、释放氧气的功能ꎬ同 时还具有分泌杀菌素、保持水土促进养分循环、调 节气候、净化空气等功能. 所以ꎬ外部生态价值损 失可以根据这些生态功能进行估算 Ceco= n + l A VfEE.4 式中Ceco为外部生态价值损失ꎬ104元ꎻn 为矿山 开采寿命ꎬaꎻl 为矿山开采结束后复垦和生态恢 复时间ꎬaꎻVfEE为单位面积林地的生态服务价值ꎬ 104元 hm -2. 3 最终境界设计 司家营露天矿位于河北省滦县ꎬ初期在东邦 设计的主运输干道由于边坡不稳定需要改设ꎬ考 虑将主运输干道布置在西帮ꎬ但西帮由于高压线 及公路约束ꎬ无法扩帮ꎬ势必造成原设计境界底部 抬高ꎬ损失矿量. 根据矿山目前的地质资料ꎬ如果 只移动高压线ꎬ境界范围变化不大ꎻ但在高压线和 公路同时移动时ꎬ最终境界相当于在 1 个无约束 状态进行优化ꎬ境界范围和矿岩量可能会发生很 大变化ꎬ产生的环境影响也会相应增加. 哪种方案 最好ꎬ需要从经济角度和环境保护角度进行评价. 本文从两方面进行考虑①考虑高压线等地表构 筑物约束条件ꎬ进行最终境界优化ꎬ即有约束优 化ꎻ②不考虑任何地表构筑物的约束ꎬ进行最终境 界优化ꎬ即无约束优化. 并分别计算两方面的生态 成本ꎬ最终选出生态 - 经济最优的境界设计方案. 目前ꎬ矿区分两个采场进行生产ꎬ如图 1 所 示. 其中 I 采场属于新建采场ꎬ各项技术经济数据 还不够稳定ꎬ所以对两个采场均使用 II 采场的数 据进行优化. 境界优化技术经济参数如表 1 所示ꎬ 其中废石混入率为 10% ꎬ废石品位为 0ꎬ精矿品位 为 65% ꎬ精矿售价为 560 元 t -1. 表 1 境界优化技术经济参数 Table 1 Technical and economic parameters used for optimizing ultimate pit 采矿成本/ 元 t -1 选矿成本/ 元 t -1 剥离成本/ 元 t -1 选矿回收率/ %采矿回收率/ % 氧化矿原生矿氧化矿原生矿剥岩剥土氧化矿原生矿氧化矿原生矿 656564627572819494 矿山不同方位的最终帮坡角如表 2 所示ꎬ考 虑如果在西帮布置坑线ꎬ帮坡角必须变缓ꎬ故在 II 采场西帮方位 180根据实际情况变缓 7相 对初始帮坡角设计ꎬ为了保证相邻方位过渡不 至于突变ꎬ所以在方位 162 5取 102 5方位与 180方位帮坡角的平均值ꎬ同样ꎬ在方位 202 5 275东北大学学报自然科学版 第 39 卷 取 180与 275的平均值. 表 2 帮坡角 Table 2 Wall slope angle 方位角030 0102 5162 5180202 5275342 5 初始帮坡角39 040 542 044 547 046 045 042 0 帮坡角39 040 542 041 040 042 545 042 0 基于上述参数设置ꎬ以地表范围约束线为限 制条件ꎬ分别考虑有约束和无约束的境界优化ꎬ优 化结果如表 3 所示. 表 3 优化结果 Table 3 Optimization results 约束 矿石量/ 亿 t岩石量/ 亿 t精矿量/ 亿 t境界范围/ hm2I 采场最终深度/ m II 采场最终深度/ m利润/ 亿元 无6 157 521 048 41 974 0421-562-412175 8 有5 261 414 721 51 678 9345-517-382168 5 由表 3 可知ꎬ受地表构筑物影响ꎬ境界矿岩 量、开采范围和开采深度均发生了较大变化. 相对 于无约束移动地表构筑物条件下的境界优化 结果ꎬ有约束不移动地表构筑物条件下优化得 到的境界中矿石量考虑废石混入率及矿石回 采率减少 0 896 1 亿 tꎬ下降 14 55% ꎻ岩石量减 少 6 326 9 亿 tꎬ下降 30 06% ꎻ矿岩总量减少 27% ꎻ精矿量矿山最终产品ꎬ考虑选矿回收率 减少 0 295 1 亿 tꎬ下降 15% ꎻ境界范围ꎬ即采区地 表最终破坏面积缩小 76 hm2ꎬ减小 18% ꎻI 采场ꎬ 即采场北区ꎬ开采深度减少 45 m上升 3 个台 阶ꎻII 采场ꎬ即采场南区ꎬ开采深度减少 30 m上 升2 个台阶. 由表3 可知ꎬ只从利润值角度考虑ꎬ 当铁精矿市场价格为 560 元/ t 时ꎬ如果高压线及 公路移动的成本小于 7 3 亿元净现值为 0 99 亿ꎬ其中生产能力 2 000 万 t/ aꎬ折现率7% ꎬ就值 得移动ꎬ否则就不移动. 优化结果的三维实体图如 图 3 所示. 图 3 三维境界模型 Fig 3 3D model of ultimate pit 根据上述优化结果ꎬ可以初步估计在有约束 和无约束方案下各自的生态成本. 矿山目前的占 地主要包括采场、尾矿库和排岩场等ꎬ其中矿山剥 离的岩石一部分40% 是外委给运输公司进行 填海ꎬ另外的部分用于修筑尾矿库和铁路ꎬ所以司 家营没有排土场. 根据选矿工艺ꎬ平均 1t 精矿的 尾矿量为 2 45 t. 参考王青等[10]关于尾矿库占地 面积计算方法尾矿品位为 6 5% ꎬ尾矿库深度为 70 mꎬ尾矿库形态系数为 1 5ꎬ尾矿沉实后的体重 为 1 75 t/ m3ꎬ 可 得 尾 矿 库 的 面 积 分 别 为 592 2 hm2无约束和 503 67 hm2有约束ꎬ司 家营 采 场 面 积 分 别 为 421 hm2 无 约 束 和 345 hm2无约束. 对于露天矿ꎬCO2的来源包括爆破、采场耗 电及选矿厂耗电. 根据炸药反应化学式ꎬ一般1 t 铵油炸药的 CO2放量为 3 tꎬ1 度电需要标准煤 0 404 kgꎬ1t 标准煤的 CO2放量为 2 49 t. 根据统 计ꎬ司家营矿露天采场每万吨采剥量的炸药消耗 为 2 2 tꎬ耗电为 0 35 kw h/ tꎻ选矿厂入选矿石耗 电为 28 07 kW h/ t. 区域复垦成本为 40 万元/ hm2假设土地复 垦和生态恢复总共时间为 5aꎬ复垦时间发生在矿 山开采结束后ꎬ每年复垦面积相同ꎬ征地成本 120 万元/ hm2ꎬ林地外部生态价值损失为 2 31 万元/ hm2[11]年矿山生产能力为 2 000 万 tꎬ所以 两种方案的外部生态价值损失会持续到第 36a 和 第 31aꎬ根据美国最新标准知碳排放社会成本大 概为 260 元/ t. 将上述参数设置代入式1 式4ꎬ可知司 家营露天矿在两种境界设计方案下的生态成本分 别为 83 32 亿元无约束ꎬ其中直接经济损失 24 32 亿元ꎬ复垦成本 4 05 亿元ꎬ能源消耗的碳 排放成本 47 70 亿元ꎬ外生生态价值损失 7 26 亿 375第 4 期 胥孝川等 多因素对露天矿最终境界设计的影响 元和 69 32 亿元有约束ꎬ其中直接经济损失 20 37 亿元ꎬ复垦成本 3 39 亿元ꎬ能源消耗的碳 排放成本 40 46 亿元ꎬ外生生态价值损失 5 10 亿 元. 综合考虑两种方案的境界利润值及对应的 矿山生态成本ꎬ得到考虑生态成本的矿山生态 - 经济整体利润分别为 92 48 亿元无约束 和 99 18 亿元有约束. 可以看出ꎬ当考虑生态成本 时ꎬ无约束境界的利润值反而比有约束的少了 6 71 亿元. 即ꎬ如果考虑环境问题的话ꎬ不值得移 动高压线等地标构造物. 上述两种方案的比较只是以当前既定高压 线、公路等为约束条件进行讨论ꎬ即不是从优化角 度考虑布置的ꎬ假设构筑物可以任意移动ꎬ那么很 有可能存在一个更合适的构筑物布置位置ꎬ使得 综合生态成本的总盈利更高ꎬ即比当前构筑物所 在位置得到的综合盈利更大. 4 结 论 1 受地表构筑物影响ꎬ最终境界在开采范 围、采剥量及开采深度等都发生了很大变化. 矿岩 量减少了 27% 矿石量减少 0 9 亿 tꎬ精矿量减 少 15% ꎬ境界开采范围缩少 18% ꎬ北区开采深度 减少 45 m3 个台阶ꎬ南区减少 30 m2 个台 阶. 其中对境界起主要影响的因素是西帮高 压线. 2 地表构筑物是否值得移动受生态环境影 响较大. 不考虑生态成本时ꎬ在当前技术经济参数 下ꎬ无约束境界盈利比受约束境界盈利高出 7 3 亿元净现值 0 99 亿元ꎬ即如果地表构筑物移 动成本低于该净现值ꎬ就值得移动ꎻ考虑生态成本 时ꎬ有约束境界利润值反而比无约束境界高出 6 71 亿元ꎬ即地表构筑物不值得移动. 参考文献 [ 1 ] Kent Mꎬ PeattieRꎬ ChamberlainV. Incorporatinggrade uncertainty in the decision to expand the main pit at the Navachab Gold MineꎬNamibiaꎬthrough the use of stochastic simulation [ C ] / / Orebody Modeling and Strategic Mine Planning.Perth AustralasianInstituteofMiningand Metallurgyꎬ2007207 -216. [ 2 ] Dimitrakopoulos RꎬMartinez LꎬRamazan S. Optimizing open pit design with simulated orebodies and Whittle four ̄x a maximum upside/ minimumdownsideapproach [ C ] / / Orebody Modeling and Strategic Mine Planning. Perth Australasian InstituteofMiningandMetallurgyꎬ 2007 201 -206. [ 3 ] Akbari A DꎬOsanloo Mꎬ Shirazi M A. Ultimate pit limit UPL determinationthroughminimizingriskcosts associated with price uncertainty[J]. Gospodarka Surowcami Mineralnymi ̄Mineral Resources Managementꎬ2008ꎬ244 157 -170. [ 4 ] Latorre Eꎬ Golosinski T S. Definition of economic limits taking into consideration time value of money [ J]. CIM Journalꎬ2011ꎬ23162 -170. [ 5 ] Frimpong Sꎬ Asa Eꎬ Szymanski J. Intelligent modeling advances in open pit mine design and optimization research [J]. International Journal of Surface MiningꎬReclamation and Environmentꎬ2002ꎬ162134 -143. [ 6 ] 张彤炜ꎬ王李管ꎬ龚元翔. 露天开采的境界优化算法研究及 应用[J]. 金属矿山ꎬ2008230 -34. Zhang Tong ̄weiꎬWang Li ̄guanꎬGong Yuan ̄xiang. Study onalgorithmsforopen ̄pitlimitoptimizationandits application[J]. Metal Mineꎬ2008230 -34. [ 7 ] 黄勇ꎬ胡庆雄ꎬ赵兴宽. 多种因素变化的露天矿境界优化 [J]. 金属矿山ꎬ20111235 -38. Huang Yongꎬ Hu Qing ̄xiongꎬ Zhao Xing ̄kuan. Open pit optimization with various changeable factors [ J]. Metal Mineꎬ20111235 -38. [ 8 ] Odell C J. Integration of sustainability into the mine design process[ D ].VancouverUniversityofBritish Columbiaꎬ2004. [ 9 ] Azapagic AꎬClift R. Life cycle assessment as a tool for improving process performance a case study on boron products [J]. International Journal of Life Cycle Assessmentꎬ 1999ꎬ43133 -142. [10]王青ꎬ胥孝川ꎬ顾晓薇ꎬ等. 金属露天矿开采的生态足迹和 生态成本[J]. 资源科学ꎬ2012ꎬ34112133 -2138. Wang QingꎬXu Xiao ̄chuanꎬGu Xiao ̄weiꎬet al. Ecological footprint and ecological costs in open ̄pit metal mining[J]. Resources Scienceꎬ2012ꎬ34112133 -2138. [11]余新晓ꎬ鲁绍伟ꎬ靳芳ꎬ等. 中国森林生态系统服务功能价 值评估[J]. 生态学报ꎬ2005ꎬ2582096 -2102. Yu Xin ̄xiaoꎬLu Shao ̄weiꎬJin Fangꎬet al. The assessment of the forest ecosystem services evaluation in China[J]. Acta Ecological Sinicaꎬ2005ꎬ2582096 -2102. 475东北大学学报自然科学版 第 39 卷
展开阅读全文