露天矿境界是判定矿山资源储量及服务年限的依据,是矿山生产设计规划的基础。随着国内外露天矿山的发展,露天境界的圈定原则与方法不断丰富,主要方法包括:传统手工法、动锥法、LG图论法、网络流法、境界动态优化法等,目前在我国的露天矿开采设计中,广泛采用n0≤nj的原则确定境界[1],即以境界剥采比不大于经济合理剥采比为原则。传统的手工方法是按照统一的经济合理剥采比,采用边帮线段比法圈定境界[2-3],因其方便、快捷、易于掌握的特点被广泛采用[4];动锥法是一种用系统模拟技术来解决露天矿开采境界的方法[5],基本出发点是采用相互交错或者相互重叠的锥体来模拟最终开采境界[6],圆锥体密集程度决定模拟露天矿的真实度,具有主观操作性;LG图论法是指将露天矿的各个方块视作节点,方块之间的开采关系用有向弧表示,利用图论搜索最大开采闭包,在有向网络图中求出盈利最大的开采境界,是具有严格数学逻辑的最终境界优化方法[7],但存在着方法复杂、计算耗时长、处理复杂边坡角困难等缺点[8];网络流法是在建立矿床经济模型的基础上,根据露天开采的几何约束条件,形成矿岩网络,应用最大流原理和标号算法,从网络中找出经济效益最佳的矿岩块段组合,作为最优露天开采境界[9];境界动态优化法是指通过编排进度计划等方式将时间属性植入矿床块体模型[10],构建集经济时间序列预测、矿岩时间属性赋值和动态经济指标计算为一体的境界全动态优化方法[11]。通过分析可知,上述方法的核心目的都是相同的,即:最大程度地实现露天开采矿山经济效益最大化。
对于倾斜煤层烧变岩侵蚀露天矿,为了充分发挥露天开采优势,尽可能实现矿山初期经济效益最大化,首采区通常选在煤层埋藏较浅、煤层赋存稳定的地段[12],即对露天矿进行分区开采,以下境界指分区境界。新疆属典型的大陆性温带荒漠气候,气候干燥,煤层赋存较浅[13],干旱、半干旱区岩层和裸露的煤层经风化大量吸热,当内存温度达到一定高温时(80℃以上)就可以引起煤本身增温燃烧;高山隆升使岩体煤层裸露,煤系地层受到剥蚀煤层出露受到干热氧化作用形成火烧区[14]。由于新疆特殊的地质及气候条件,导致该地区煤田的火烧区点多面广,燃烧速度快,煤炭储量损失多,经济损失严重。火烧岩主要受煤层烘烤以及大地构造的影响[15],一般垂向上是下部裂隙发育较强,上部岩体裂隙较弱[16],且倾角大,烧变岩的存在对边坡稳定性产生影响,在矿山开采过程中,为了保证边帮稳定,需将其全部清掉,这样便大幅增加了开采剥离量,使得境界内平均剥采比变大,甚至导致平均剥采比大于经济合理剥采比,致使按境界剥采比小于或等于经济合理剥采比原则圈定的境界不合理。因此,对于烧变岩侵蚀露天矿,不能简单地按照境界剥采比小于或等于经济合理剥采比的原则圈定境界,应综合考虑平均剥采比的变化来圈定其最优境界。
笔者分析了新疆煤田烧变岩赋存特点,并以神新能源黑山露天煤矿为背景,分析烧变岩的侵蚀对露天矿境界圈定的影响;通过构建模拟境界位置,分析境界剥采比与平均剥采比随分区境界外扩的变化规律,确定倾斜煤层烧变岩侵蚀露天矿境界圈定的原则及适用条件;以边帮线段比法为依据,采用拥有自主知识产权SMCAD软件构建了黑山露天煤矿三维仿真境界;在真三维的情况下,在圈定的仿真境界内填充煤、岩层状属性块体,采用块量反验的方式校验境界位置的合理性,确定黑山露天煤矿的最优境界。
黑山露天煤矿位于天山中段以北的山间谷地,地势北高南低,西高东低。勘探区内发育的可采煤层包括:9煤、11煤、13-2煤;其中,9煤平均厚度为9.70 m,密度1.31 t/m3、11煤平均厚度为9.55 m,密度1.32 t/m3、13-2煤平均厚度为29.95 m,密度1.34 t/m3,煤层稳定赋存倾角为12°。烧变岩是地表附近的厚煤层在地质历史时期自燃,使围岩受到烘烤或烧烤而形成的特殊类型的岩体[17],黑山露天煤矿露头处围岩历经13-2煤自燃区的烧烤蚀变作用,形成了一条近东西向的烧变岩带,地表出露的烧变岩主要由13-2号煤层自燃后形成,宽200~1 200 m,长约18 km,垂深为60~240 m,呈带状横贯全区,是典型的倾斜煤层烧变岩侵蚀露天矿。
通过分析可知:黑山露天矿采场北帮烧变岩边帮构造形态为一宽缓的NE向不对称背斜,倾角(13°~17°),是典型顺倾边坡;且烧变岩底板存在软弱夹层,由于弱层的抗剪强度较低,极大地降低了岩土体的抗滑稳定性,使得黑山露天矿北帮边坡滑坡风险加大,对于这种底板含有软弱夹层的松散均质顺倾边坡,其滑坡模式应由弱层和圆弧面组合控制,由此可以确定,黑山露天矿采场北帮边坡的潜在滑坡模式主要为切层-顺层滑动,即边坡可能发生沿着以圆弧或以弱层为底界面的切层-顺层滑动[18]。为了确保北帮烧变岩边坡稳定性,黑山露天矿北帮烧变岩侵蚀边坡安全储备系数取1.2。依据Mohr-Coulomb准则,同时考虑剪切破坏和拉伸破坏,选取典型滑面借助FLAC3D3.0对北帮烧变岩进行边坡稳定性分析,基于强度折减原理,循环折减岩体强度直到边坡刚好达到临界失稳状态[19-20],以此获得边坡的稳定系数Fs,如图1所示、表1所示。其中,烧变岩容重24.7 kN/m3,黏聚力0.01 MPa,内摩擦角11°。
图1 北帮烧变岩边坡稳定性分析
Fig.1 Stability analysis of north slope burnt rock slope
表1 北帮烧变岩典型滑面稳定性分析
Table 1 Stability analysis of typical sliding surface of burnt rock in north slope
典型滑面滑面1滑面2滑面3滑面4滑面5稳定系数1.2051.2031.2111.2091.213
分析图1、表1并结合边坡稳定性分析可知,黑山露天矿在开采13-2煤的同时为了保证北帮边坡稳定,必须将火烧岩全部清除,并对烧变岩底板的软弱岩层进行处理,然后再以烧变岩底板作为北帮开采境界。
露天矿圈定境界的原则包括:①境界剥采比nk小于或等于经济合理剥采比nj;②平均剥采比等于经济合理剥采比;③生产剥采比等于经济合理剥采比。由于原则①是从全矿床或全露天矿成本及盈利等指标最优为基础得出的,而且此原则在露天矿设计过程中具有较强的可操作性[21];因此,各露天矿在圈定开采境界时,以原则①为依据设计露天开采境界。
图2 顶帮境界模拟位置
Fig.2 Simulation drawing of top slope boundary
由于黑山露天矿北帮已确定以烧变岩底板为界,端帮境界一帮依据矿权界决定;因此,顶帮境界位置是影响其境界大小的关键。当顶帮境界位置满足原则①,即nk≤nj时为确定的合理顶帮境界;由图2剖面信息可知,13-2煤层露头处存在大量烧变岩,在开采13-2煤的同时须将其全部清除,使得此境界的剥离负荷量大大增加,平均剥采比np满足:np>nk且np>nj;所以该顶帮境界位置不能作为合理的境界位置。因此,对于受烧变岩侵蚀严重的露天矿,不可简单按照原则①圈定境界,应当综合考虑平均剥采比的变化来确定最优境界。
为了确定倾斜煤层烧变岩侵蚀露天矿顶帮分区境界外扩过程中,根据境界剥采比与平均剥采比的变化规律,得出倾斜煤层烧变岩侵蚀露天矿境界圈定原则,依照黑山露天矿到界边帮参数构建模拟最终帮,不断模拟降深单位深度Δh(即一个台阶高度),外扩水平距离为x构建模拟境界位置,如图3所示,计算不同境界位置i处的境界剥采比nk,i与平均剥采比np,i满足以下条件:
(1)
(2)
(3)
式中:θ为煤层倾角,(°);e为总境界数;f为总水平数;Vi为境界位置i剥离总量,m3;Mi为境界位置i可采煤量,m3;Lij为第i个境界位置第j水平岩两端帮之间的距离,m;Mij为第i个境界位置第j水平煤距端帮的距离,m;xij为第i个境界位置第j水平顶底帮之间的距离,m;;lij为第i个境界位置第j水平岩边帮厚度,m;mij为第i个境界位置第j水平煤边帮厚度,m;γ为煤容重,kN/m3。
图3 烧变岩侵蚀露天矿模拟境界位置
Fig.3 Simulation boundary location of open-pit mine by eroded burnt rock
黑山露天矿设计台阶高度Δh为10 m,将其代入式(1),可知:外扩水平距离Δx为47 m。依照上述方法,构建模拟境界位置,由式(2)、式(3)计算境界剥采比nk,i与平均剥采比np,i,见表2。
鉴于表2中模拟境界位置与境界剥采比nk,i、平均剥采比np,i为连续函数关系,通过回归分析,确保境界位置与境界剥采比nk,i、平均剥采比np,i近似符合多项式分布,即nk,i满足式(4),np,i满足式(5),得
nk(x)=-1.381×10-4+1.039x-1.543x2+1.097x3-0.334x4+0.055x5-0.055x6
(4)
np(x)=10.778-2.019x-0.457/x+0.162x2+0.652x2-0.005x3-0.101x3
(5)
其中,式(4)、式(5)的相关性系数R2都达到0.99,趋于1,可知本回归分析可靠度满足要求。为了确定烧变岩侵蚀露天矿境界圈定的原则,以nk(x)、np(x)为基础,通过绘制境界剥采比与平均剥采比随境界位置变化关系曲线,建立经济合理剥采比评价模型,如图4所示。
表2 境界剥采比与平均剥采比变化量
Table 2 Change scale of boundary stripping ratio and the average stripping ratio
参数模拟境界位置01234567891011境界剥采比/(m3·t-1)0.000.320.821.613.545.166.237.869.5310.3611.2012.43平均剥采比/(m3·t-1)无穷大9.547.046.094.853.953.413.062.932.031.451.29
图4 境界剥采比与平均剥采比变化曲线
Fig.4 Change curves of boundary stripping ratio and average stripping ratio
分析图4可知:随着顶帮境界不断外扩,使得境界剥采比nk不断加大,平均剥采比np逐渐减小。在境界位置x0处,有n0为4.4 m3/t,使得nk=np=n0成立。给定未知经济合理剥采比nj,该评价模型必然存在以下3种假设:
1)假设一:nj>n0,曲线有2个交点,存在境界位置X1、X2二个评价位置。①在境界位置X1处,有nk≤nj且np=nj,但是在境界位置X1之后,平均剥采比还在不断减小,且仍满足nk≤nj且np≤nj;因此,境界位置X1不能作为合理境界位置;②在境界位置X2处,有nk=nj且np≤nj,在境界位置X2之前,满足 nk<nj且np<nj;在境界位置X2之后,满足:nk>nj且np<nj;因此,境界位置X2是合理境界位置。
2)假设二:nj=n0,曲线有且仅有一个交点,存在境界位置X0一个评价位置。在境界位置X0处,有nk=np=nj,并且在境界位置X0之前,有np>nj,在境界位置X0之后,有nk>nj;因此,境界位置X0是合理境界位置。
3)假设三:nj<n0,曲线有2个交点,存在境界位置X3、X42个评价位置。①在境界位置X3处,虽有nk=nj;但是使得np>nj成立;因此,境界位置B1不能作为合理境界位置。②在境界位置X4处,虽有np=nj;但是使得nk>nj成立;因此,境界位置X4不能作为合理境界位置。
综上所述,只有当经济合理剥采比nj满足:nj≥n0时,烧变岩侵蚀露天矿存在以境界剥采比小于或等于经济合理剥采比的原则圈定的合理境界;当经济合理剥采比满足:nj<n0时,烧变岩侵蚀露天矿不存在经济合理境界。
经济合理剥采比的确定方法基本上可分2类:①是成本相等法;②是价格法。前者适用于露天和地下开采的贫化损失相同或近似,矿石资源丰富且不很贵重时,但从实际生产看,露天开采和地下开采在贫化损失方面相差很大,露天开采采出率可达90%以上,而地下开采回采率仅为50%左右[18]。因此,各露天矿在圈定开采境界时,以市场条件为依托,采用后者计算经济合理剥采比,即
(6)
式中:dl为原煤售价,元/t;a为纯采煤成本,元/t;b为纯剥离成本,元/m3。
通过市场调研及成本预算可知:黑山露天煤矿2019年一季度煤炭售价dl为:132.61元;采煤成本a为12.52元/t;剥离成本b为:14.90元/m3。将其代入式(6),得黑山露天煤矿经济合理剥采比nj为8.06m3/t。
综上可知,黑山露天矿经济合理剥采比nj满足:nj>n0=4.4m3/t,可以按照原则①圈定境界。依据黑山露天煤矿提供的勘探线地质剖面图,选取典型剖面采用边帮线段比法确定剖面最优边帮位置,按照黑山露天矿到界边帮参数绘制虚拟最终帮;根据煤、岩层层位信息通过三维坐标转换将剖面对应到真实的三维空间地理位置;在真三维空间背景下,将所绘制的各剖面虚拟最终帮按照标高信息通过空间位置关系赋值依次拼接起来,构建依附实际的黑山露天煤矿三维仿真境界,如图5所示。但是,由于煤层赋存情况有时会出现发育不稳甚至出现尖灭,三维仿真顶帮境界只能保证剖面位置满足nk≤nj,对于整个顶帮境界而言是不确定的;因此,需要对三维仿真顶帮境界的最优性进行校验。
图5 黑山露天矿三维模拟境界
Fig.5 3D simulation boundary of Heishan Open-pit Coal Mine
将露天矿三维模拟境界以地表界为边界约束条件生成DTM面模型,这样所圈定的仿真境界相当于是一个内空的盒子,以盒子边界为约束采用块体填充满整个开采境界,并建立矿岩属性数据库(包括:岩性、容重等);将建立好的各煤层顶底板面模型约束块体,以块体质心坐标(xgh、ygh、zgh)为约束标准,假如块体质心划入块体约束条件之内,块体计入算量,反之剔除块体,将矿岩属性数据赋予层状块体模型,构建层状属性块体地质模型,如式(6)、式(7)、图6。境界内的块体总量相当于就是整个采区要开采的矿岩总量,煤层底板与煤层顶板之间块量相当于是煤量,剩余量为岩量,如式(8)、式(9)。假定储量计算过程母块尺寸大小为:L×W×H(长×宽×高),子块大小为:uL×vW×wH,u、v、w为块体缩放比例;块体的大小直接影响储量计算快慢,块体越大计算速度越快,块体越小计算速度越慢,但是块体越小精度越高,通常母块大小为5×5×2,缩放比例取0.5[22-24]。
bij(xgh,ygh,zgh)=Lgh×Wgh×Hgh
(6)
(7)
(8)
(9)
式中:xg,min第g水平x方向约束条件最小值;xg,max第g水平x方向约束条件最大值;xgh第g水平第h个块体x方向坐标;yg,min第g水平y方向约束条件最小值;yg,max第g水平y方向约束条件最大值;ygh第g水平第h个块体y方向坐标;zg,min第g水平z方向约束条件最小值;zg,max第g水平z方向约束条件最大值;zgh第g水平第h个块体z方向坐标。Vr为境界内总岩量,m3;m为岩块体数量,个;bg为第g水平第h个块体体积,m3;Pc为境界内总煤量,t;n为煤块体个数,个;Lgh为第g水平第h个块体长度,m;Wgh为第g水平第h个块体宽度,m;Hgh为第g水平第h个块体高度,m。
图6 开采境界块量约束示意
Fig.6 Constraints schematic of mining boundary block quantity
针对三维仿真顶帮境界位置的模糊不确定性,笔者提出以块量反验的方法校验顶帮境界的合理性。依据境界剥采比的定义:下降单位深度(或一个台阶高度增加)的岩矿量之比。将得到的顶帮境界里缩水平距离Δx,依照内缩顶帮与原顶帮之间的矿岩量之比得到块量里收境界剥采比为同样,将得到的顶帮外扩水平距Δx,依照外扩顶帮与原顶帮之间的矿岩量得到块量外扩境界剥采比的如果有则说明所圈定的顶帮位置合理的;如果有则说明所圈定的顶帮位置过大;如果有则说明所圈定的顶帮位置过小。其中,为第i个里收境界的块量境界剥采比;为第i个外扩境界的块量境界剥采比。
图7 块量境界剥采比约束示意
Fig.7 Constraints schematic of block quantity boundary stripping ratio
依据上述方法,得到黑山露天煤矿煤岩层层状块体模型,如图8所示。
图8 黑山露天煤矿煤岩层层状块体模型
Fig.8 Coal rock layer block model of Heishan Open-pit Coal Mine
将黑山露天矿圈定的顶帮境界不断里收水平距离Δx为47 m,得到其块量里收境界剥采比同样,将圈定的顶帮分区境界外扩水平距离47m,得到其外块量扩境界剥采比可知:黑山露天矿顶帮块量境界剥采比变化规律,见表3和如图9所示。
表3 块量境界剥采比变化
Table 3 Blockquantity boundary stripping ratio change
境界位置煤量/万t剥离量/万m3境界剥采比/(m3·t-1)里收31 115.136 208.095.57里收21 047.036 994.376.68里收11 040.957 255.836.97设计1 058.037 787.887.36外扩11 117.748 994.838.05外扩21 093.4810 077.589.22外扩31 051.1311 363.3810.81
图9 块量境界剥采比变化曲线
Fig.9 Curves of block quantity boundary stripping ratio
分析表3、图9可知:外扩境界1块量外扩境界剥采比满足:说明通过剖面边帮线段比法圈定的黑山露天矿仿真境界变小,外扩境界1才是最为合理的境界位置。
1)倾斜煤层烧变岩侵蚀露天矿,随着境界的扩大,满足:境界剥采比逐渐增加,平均剥采比逐渐变小;只有当经济合理剥采比满足:nj≥n0时,存在以境界剥采比小于或等于经济合理剥采比的原则圈定合理境界;当nj<n0时,烧变岩侵蚀露天矿不存在合理境界,该原则同样适用新疆地区普遍存在的火烧岩侵蚀露天矿。
2)针对新疆煤田火烧区特点,以及煤、岩层赋存情况的随机性,使得采用剖面边帮线段比法通过三维坐标转换确定的仿真境界具有模糊不确定性;通过构建煤岩层层状属性块体模型,采用块量反验方法得到其块量里收境界剥采比与块量外扩境界剥采比,能够从三维角度实现反验设计境界合理性的目的。
3)通过分析可知,黑山露天煤矿经济合理剥采比为8.06 m3/t,满足:n0=4.4 m3/t≤nj,即存在以境界剥采比小于或等于经济合理剥采比原则圈定的合理境界;由块量反验可知,因此,外扩境界1是黑山露天煤矿的最优境界。
[1] 赵文斌,李小武.谈谈用传统方法确定露天境界[J].有色金属设计,2009,36(4):37-41.
ZHAO Wenbin,LI Xiaowu.Talking about determining the open pitlimit by traditional method[J].Nonferrous Metals Design,2009,36(4):37-41.
[2] ALFORD C G,WHITTLEJ.Application of Lerchs-Grossman pit optimization to the design of open pit mines[J].Symposia Series-Australasian Institute of Mining and Metallurgy,1986,12(8):201-207.
[3] HERNER B.Neutral sizing system allows fasteasy conversion for acid machine[J].Pulp and Paper,2001,64(1):154-163.
[4] 刘 闯,白润才,刘光伟,等.复合煤层露天矿下部煤层开采境界优化[J].煤炭学报,2015,40(9):2085-2090.
LIU Chuang,BAI Runcai,LIU Guangwei,et al.Optimization of mining limits of deep-seated coal seams in multi-seams surface coal mine[J].Journal of China Coal Society,2015,40(9):2085-2090.
[5] 贾海波.基于实体模型的露天矿境界圈定动锥法研究[D].沈阳:东北大学,2008:32-34.
[6] 马 帅.基于浮动圆锥法的三道庄露天矿境界优化[D].西安:西安建筑科技大学,2013:54-67.
[7] 张彤炜.基于图论学的露天开采境界优化算法研究及程序设计[D].长沙:中南大学,2008:44-52.
[8] 张延凯,李克庆,胡乃联,等.露天矿境界优化 LG 算法初始有向图生成研究[J].煤炭学报,2015,40(S2):367-373.
ZHANG Yankai,LI Keqing,HU Nailian.et al.A study of generating the initial directed graph for LG Algorithm of open-pit optimization[J].Journal of China Coal Society,2015,40(S2):367-373.
[9] 陈孝华.用网络最大流标号算法进行露天开采境界的优化设计[J].云南冶金,1982(3):3-6.
CHEN Xiaohua.Network maximal flow for optimal design of open-pit mining limits[J].Yunnan Metallurgy,1982(3):3-6.
[10] 杨 彪,罗周全,廖江南,等.露天矿山开采境界动态优化圈定[J].煤炭学报,2010,35(3):368-372.
YANG Biao,LUO Zhouquan,LIAO Jiangnan,et al.Dynamic optimization for open pit mining limit[J].Journal of China Coal Society,2010,35(3):368-372.
[11] 杨 彪,罗周全,陆 广,等.基于经济时间序列预测的露天矿开采境界动态优化[J].煤炭学报,2011,36(1):29-33.
YANG Biao,LUO Zhouquan,LU Guang,et al.Open pit mining limit dynamic optimization based on economic time series forecasting[J].Journal of China Coal Society,2011,36(1):29-33.
[12] 曹 博,陶亚彬,白润才,等.倾斜煤层窄长露天矿分区开采分期境界优化[J].重庆大学学报,2019,42(4):101-110.
CAO Bo,TAO Yabin,BAI Runcai,et al.Block system optimization by stages for inclined seam long narrow open-pit mine[J].Journal of Chongqing University,2019,42(4):101-110.
[13] 田 军,李海洲,夏 冬.新疆白砾滩露天煤矿烧变岩特征及边坡稳定性分析[J].煤矿安全,2013,44(7):223-225.
TIAN Jun,LI Haizhou,XIA Dong.Burnt rock characteristics and slope stability analysis in Xinjiang Bailitan open pit mine[J].Safety in Coal Mines,2013,44(7):223-225.
[14] 张秀山.新疆煤田火烧区特征及防治对策[J].新疆地质,2001(2):150-152.
ZHANG Xiushan.Characteristics of coal burning area in Xinjiang coalfield and its prevention and control measures[J].Xinjiang Geology,2001(2):150-152.
[15] 刘 侃.新疆准南煤田某煤矿火烧区的水文地质特征[J].内蒙古煤炭经济,2017(19):141-144.
LIU Kan.Hydrogeological characteristics of a coal mine burning area in Zhunnan Coalfield Xinjiang[J].Inner Mongolia Coal Economy,2017(19):141-144.
[16] 范立民,仵 拨,向茂西,等.我国西部保水采煤区受保护烧变岩含水层研究[J].煤炭科学技,2016,44(8):1-6.
FAN Limin,WU Bo,XIANG Moxi,et al.Study on protective burnt rock aquifer in water preserved coal mining area of western China[J].Coal Science and Technology,2016,44(8):1-6.
[17] 贾文凯,梁小山,张 萌.李家沟煤矿烧变岩的分布规律研究及其水文地质意义[J].科学技术创新,2018(10):13-14.
JIA Wenkai,LIANG Xiaoshan,ZHANG Meng.Study on distribution law of burnt rocks in Lijiagou coal mine and its hydrogeological significance[J].Scientific and Technological Innovation,2018(10):13-14.
[18] 曹兰柱,白继元,王 东,等.含顺倾弱层软岩边坡三维稳定性影响因素分析[J].金属矿山,2017(1):174-178.
CAO Lanzhu,BAI Jiyuan,WANG Dong,et al.Analysis on three dimensional stability influence factors of soft rock slope with dip weak layer[J].Metal Mine,2017(1):174-178.
[19] 陈建平.边坡强度折减法判据及点安全系数研究及工程应用[D].赣州:江西理工大学,2012:56-63.
[20] 袁 灿.基于强度折减法的浅变质岩风化层边坡稳定性研究[D].武汉:武汉理工大学,2008:30-56.
[21] 王 岩,王为举.露天煤矿开采境界与经济剥采比确定[J].现代矿业,2011,31(12):83-85.
WANG Yan,WANG Weiju.Open-pit mining boundary and economic stripping ratio is certain[J].Modern Mining,2011,31(12):83-85.
[22] 黄俊歆,王李管,毕 林,等.改进的露天境界优化几何约束模型及其应用[J].重庆大学学报,2010,33(12):78-83.
HUANG Junxin,WANG Liguan,BI Lin,et al.The improved geometrical constraint model of open-pit mine boundary optimum and its application[J].Journal of Chongqing University,2010,33(12):78-83.
[23] 陈三明,何玉州,罗文敏,等.基于地质三维模拟的攀枝花铁矿块体模型储量估算[J].桂林理工大学学报,2013,33(4):610-615.
CHEN Sanming,HE Yuzhou,LUO Wenmin,et al.Based on geological 3D simulation of Panzhihua iron ore reserve estimation block model[J].Journal of Guilin University of Technology,2013,33(4):610-615.
[24] 王洪江,杨柳华,吴爱祥,等.基于块体模型的储量估算方法[J].有色金属,2014,66(6):87-91.
WANG Hongjiang,YANG Liuhua,WU Aixiang,et al.Application of reserves estimation method based on block model[J].Non-ferrous Metal,2014,66(6):87-91.