扎哈淖尔露天矿北帮边坡形态优化设计

2018-10-31白润才王东旭

中国地质灾害与防治学报 2018年5期

白润才，王东旭，王珍

(辽宁工程技术大学矿业学院，辽宁阜新 123000)

人类遇到的边坡种类众多，且每一种边坡都有其自身的特点，露天矿采场边坡的显著特点之一就是对地质条件无可选择的余地，不能因地质条件的不良而改变其位置[1-3]，因而其稳定性往往受一系列不良地质条件的影响，如地下水、软弱夹层及软岩等。近年来，国内外众多学者对复杂地质条件下的边坡稳定性分析及边坡空间形态优化设计做了大量研究，如王东等[4]以元宝山露天矿东帮含断层顺倾软岩边坡为工程背景，研究分析了断层作用下的顺倾软岩边坡稳定性及破坏机理；曹兰柱等[5]分析了边坡开挖过程中元宝山露天矿南帮边坡坡体内应力及位移随时间的边坡规律及变形破坏机理，最终确定南帮陡帮开采边坡空间形态；曾胜等[6]通过相似模拟分析了红砂岩顺层边坡在降雨条件下结构面强度及边坡稳定性变化规律。尽管众多学者对各种复杂地质条件下的边坡稳定性分析及边坡形态优化设计进行了一系列研究，但对于逆倾层状软岩边坡陡帮开采到界边坡形态优化设计却少之又少[7]。

扎哈淖尔露天煤矿境界呈窄长型，可采煤层多、厚度相对较小、倾角相对大，且矿田内断层多、构造复杂，岩体力学强度低，地下水较为丰富，存在泥岩弱层，遇水软化，易引起滑坡；该矿北帮(工作帮)设计到界边坡角27°，生产中由于上部岩层松软，在开采靠帮过程中边坡变形较大，目前北帮860水平以上边坡无法达到设计边坡角度，上部边坡角约为15°，导致采场下部台阶无法按照设计角度进行靠帮，随着采场重心逐步西移，东侧北帮没有到界不能进行内排，边坡暴露长度逐步增大，北帮边坡稳定性问题对采场北帮推进产生的不利影响将会更加明显。基于此，本文以扎哈淖尔露天煤矿北帮边坡为工程背景，采用刚体极限平衡与数值模拟相结合的方法，优化了北帮陡帮开采到界边坡形态。

1 边坡工程地质特征及计算方法的选取

1.1 边坡工程地质特征

扎哈淖尔露天北帮地层自上而下依次为：第四系、新生界第三系、晚侏罗-早白垩世霍林河群含煤岩系和晚侏罗世兴安岭群火山岩系共同构成。影响本区水文特征的含水层有两个，分别为第四系上更新、全新统(Q3+4)含细砂砾含水层；下部含煤段(J3～k1h3)煤、裂隙含水～岩风化孔隙含水带。研究区地下补给水主要来源有三部分，一是大气降水的直接渗入和南部火山岩地下水的侧向渗入，二是第四系(Q3+4)砂砾石层地下水直接补给，三是与第三系砂砾石含水层互为补给。工程地质补勘发现43煤底板的夹矸中高岭石含量高达43.5%～65.2%，考虑到43煤底部含水，可断定43煤层底板含有软弱夹层。图1为计算剖面平面位置图，地质剖面采场现状见图2红色线段，岩土体物理力学指标见表1。

图1 典型地质剖面平面位置图Fig.1 Plane graph of typical geologic sections

表1 岩土体物理力学指标Table 1 Soil and rock physical and mechanical indicators

1.2 边坡稳定性分析方法的选取

边坡稳定性分析经历百余年发展，主要形成了工程地质类比法、极限平衡分析法、数值模拟计算法以及现代数学分析法[8-11]。其中，工程地质类比法对岩土工作者经验性要求较高。极限平衡法是现阶段较为常用的计算方法，具有计算方法简单，计算结果准确等优势，但该方法未能充分考虑边坡岩土体内部应力-应变形态。针对于极限平衡分析法存在的不足，近年来随着计算机技术的高速发展，基于有限元、有限差分等技术手段开发形成的数值模拟软件大量应用于边坡稳定性分析中。该方法能够有效地解决岩土体应力-应变关系。现代数学分析法是通过数学手段求解岩土体动力学问题，但现阶段该方法尚未形成完整的理论体系[12-13]。

综合考虑上述各类方法的优缺点及其适用条件，针对于扎哈淖尔露天矿北帮工程现状，本文拟采用极限平衡分析法与数值模拟计算法，对扎哈淖尔露天矿北帮压煤回采过程中边坡稳定性进行深入研究，在保证边坡稳定性的前提下最大限度地回收底板煤。

2 北帮到界边坡形态优化设计

扎哈淖尔露天矿北帮采用纵采方式进行降深，边坡的暴露时间小于10年，按照GB 50197—2005《煤炭工业露天矿设计规范》要求，安全储备系数宜为1.0～1.2，考虑到北帮北侧为排土场，同时北帮地表布设连续系统，最终确定北帮安全储备系数为1.2。

该矿北帮设计到界边坡采用组合台阶形式，工作平盘宽度为40 m，保安平盘宽度为10 m，台阶高度为12 m，台阶坡面角为70°，原设计到界边坡角为27°。生产中由于上部岩层松软，在开采靠帮过程中边坡变形较大，目前北帮11煤以上边坡无法达到设计边坡角度，上部边坡角约为15°，导致采场下部台阶无法按照设计角度进行靠帮，大量原煤无法采出。为了能够安全经济地回采北帮下部压煤，以提高露天矿经济效益，需重新设计能够满足安全要求的最优边坡形态。考虑会发生的潜在滑坡模式为圆弧滑动和以43煤层底板弱层为底界面的切层-顺层滑动，笔者拟采用逐阶段的优化方法，分别对各剖面11煤底板以上、11煤底板至43煤底板两个阶段稳定边坡形态参数进行了优化设计。结果如图2及表2所示。

图2 各剖面边坡形态优化结果图Fig.2 The optimized slope of each profile

表2 各剖面边坡形态优化结果图Table 2 The optimized slope of each profile

由计算结果可知，边坡稳定性随边坡角度的增加逐渐降低，1剖面、2剖面、3剖面、4剖面11底板以上边坡潜在滑坡模式为剪切圆弧滑动，稳定帮坡角分别为22°、21°、21°、21°；43煤底板以上边坡潜在滑坡模式为剪切圆弧-弱层相结合的组合滑动，11煤底板至43煤底板间稳定帮坡角分别为37°、36°、36°、36°，整体稳定帮坡角分别为29°、30°、28°、28°，边坡高度分别为276 m、348 m、343.5 m、353 m。考虑到还有爆破震动等多种因素会对边坡稳定性产生影响，为了保证北帮稳定性，推荐陡帮角度为36°，11底板以上边坡角为21°，边坡形态为双折线型凸边坡，这样既可以保证北帮边坡具有足够的稳定性又能够尽可能的回收资源。按照上述所确定的最终边坡形态，采用平面算量方法，较原设计整体边坡角27°共计可多回采压煤9.827 6×106t，剥离2.764 24×106m3。

3 数值模拟

在边坡稳定性分析领域应用最早的数值模拟方法就是有限单元法，并且也是截止到今使用较为广泛的一种的数值方法。这种方法解决了极限平衡法中刚体假设的弊端，通过应力、应变规律分析其边坡的破坏变形过程，可以确定其最危险滑裂的位置范围等等。因此，笔者采用大型岩土工程分析软件FLAC3D对扎哈淖尔露天矿北帮陡帮开采边坡进行数值模拟研究，分析它的滑坡机理和变形破坏特征。模型边界条件为：模型顶部和边坡面为自由边界、两侧及底部边界固定不动、采用重力加载。应用ANSYS软件建立边坡计算模型并提取节点数据和单元数据，通过ANSYS-TO-FLAC3D软件将提取的节点数据和单元数据导入到FLAC3D中进行迭代计算求解。

基于强度折减理论，对岩土体物理力学参数(黏聚力c和内摩擦角φ)进行循环折减，计算终了时各剖面折减系数分别为1.256、1.267、1.251、1.243，边坡处于失稳或破坏状态，因此确定1、2、3、4剖面边坡稳定系数分别为1.256、1.267、1.251、1.243。图3、图4分别为边坡临界失稳状态下各剖面水平与竖直应力云图。从图中可知边坡岩土体在自重应力与泊松效应的作用下，水平应力与竖直应力均随埋深的增加逐渐增大且顺着坡面方向呈层状一直延伸到坡脚处。在边坡表层附近出现有拉应力，由边坡岩土体自身性质可知其受压但不受拉，因此在边坡表面附近岩土体易发生拉张破裂，说明边坡岩土体变形破坏方式为张拉-剪切复合型破坏。

图3 水平应力云图Fig.3 The horizontal stress nephogram

图4 竖直应力云图Fig.4 The vertical stress nephogram

图5和图6分别为各剖面边坡临界失稳时的总位移图、剪切应变增量云图。由图5可以看出，边坡滑坡模式主要受43煤底板弱层的影响，易发生以煤层底板弱层为底界面的切层-顺层滑动，边坡失稳时滑面的形状及位置，其滑坡模式及滑面位置与刚体极限平衡计算结果基本吻合。由图6可以看出剪切应变增量显著的区域集中在煤层底板弱层，而边坡岩土体塑性屈服是边坡失稳的直接原因，因此剪切破坏对岩土体的破坏类型起到决定性的作用。