软弱基底内排土场支挡煤柱的力学效应

摘要:留设支档煤柱是提高软弱基底内排土场稳定性的重要手段之一。为了解决软弱基底内排土场支挡煤柱形态参数设计难题，实现矿山安全生产和经济效益兼顾，对软弱基底内排土场支挡煤柱的力学效应及形态设计方法进行研究。在构建软弱基底内排土场支挡煤柱力学模型的基础上，进行三维受力状态分析，将支挡煤柱视为在内排土场静止土压力及底界面剪切力双重作用下的固定梁，以摩尔-库伦准则作为失稳判据，推导获得了煤柱三维支挡效应的表达式；基于对煤柱三维支挡效应的力学成因类型及影响因素的认识，提出了煤柱支挡效应的二维等效方法，推导了底界面等效黏聚力cd与等效内摩擦角φd的表达式，从理论上揭示了支挡煤柱形态参数对支挡效应的影响。将该方法在胜利东二露天煤矿内排土场支档煤柱设计中进行了应用，并采用数值模拟方法进行验证。研究结果表明: 煤柱底界面的等效黏聚力cd及等效内摩擦系数tan φd与其顶宽a、下底角β,ω均呈正相关反比例函数关系，与高度H呈正相关倒差函数关系，与长度d呈负相关反比例函数关系；胜利东二露天煤矿内排土场边坡稳定系数Fs随支档煤柱顶宽a增加而增加，与煤柱高度H呈线性递增关系；煤柱的最佳尺寸为高度H=28 m，顶宽a=43 m，与传统的二维分析结果相比，可多才出原煤75 600 m3，经济效益显著。成果可用于优化设计软弱基底内排土场支挡煤柱形态参数，弥补采用传统方法设计造成的煤炭损失过大的不足，实现经济效益最大化。

关键词:软弱基底;内排土场;煤柱;支挡效应;等效黏聚力;等效内摩擦角

王东，李广贺，曹兰柱，等.软弱基底内排土场支挡煤柱的力学效应[J].煤炭学报,2019,44(3):933-940.doi:10.13225/j.cnki.jccs.2018.0422

WANG Dong,LI Guanghe,CAO Lanzhu,et al.Mechanical effect of supporting and retaining coal pillar on internal waste dump with weak basement[J].Journal of China Coal Society,2019,44(3):933-940.doi:10.13225/j.cnki.jccs.2018.0422

基金项目:国家自然科学基金资助项目(51874160,U1361211);辽宁省高等学校创新人才支持计划资助项目(LR2016039)

作者简介:王东(1978—)，男，辽宁法库人，教授，博士生导师。E-mail:lntu_wd@163.com

通讯作者:李广贺(1992—)，男，内蒙古通辽人，博士研究生。E-mail:li_7118@126.com

Mechanical effect of supporting and retaining coal pillar on internal waste dump with weak basement

(College of Mining Engineering,Liaoning Technical University,Fuxin 123000,China)

Abstract:Supporting and retaining coal pillar is one of the important means to improve the stability of internal waste dump with weak basement.To solve the problem of designing morphological parameters of supporting and retaining coal pillar on internal waste dump with weak basement,and to achieve the safe production and economic benefits of the mine,the mechanical effects and shape design methods were studied.Based on the mechanical model,the three-dimensional force state analysis was carried out.The supporting and retaining coal pillar was regarded as a fixed beam under the dual action of the static earth pressure of internal waste dump and the shear force at the bottom interface.The Mohr-Coulomb criterion was used as the instability criterion,authors derived the expression of three-dimensional supporting and retaining effect of coal pillar.Based on the type of mechanical genesis and influencing factors,a two-dimensional equivalent method for the supporting and retaining effect of coal pillars was proposed and the expression of cd and φd were derived.In theory,the influence law of morphological parameters on the supporting and retaining effect was revealed.The method was verified by numerical simulation the internal dump Shengli East No.2 open-pit mine.The results showed that cd and φd are positively inversely proportional to their top width a and lower bottom angles β and ω, positive correlat with the height H,and negative correlat with the length d The stability coefficient Fs of the internal dump of Shengli East No.2 open-pit mine increases with the increase of a,and increases linearly with H The optimal size of the coal pillar is 28 m and 43 m.Compared with the traditional two-dimensional analysis results,more than 75 600 m3 of raw coal can be produced,and the economic benefit is remarkable.The achievement optimizes the morphological parameters of supporting and retaining coal pillar to make up for the inadequacy of excessive coal losses caused by the traditional and to achieves a maximum economic benefit in mining operation.

Key words:weak basement;internal waste dump;coal pillar;supporting and retaining effect;equivalent cohesion;equivalent internal friction angle

内排土场是近水平、缓倾斜煤层露天开采的必然产物，其稳定性关系到矿山安全高效生产。国内如霍林河、宝日希勒、安家岭等大型露天煤矿的内排土场均发生过滑坡或大变形[1-4]，其主要诱发因素是基底岩体遇水软化形成的演化弱层[5-6]。由于这类排土场基底岩体亲水性强，采用常规的疏干排水措施的效果很难量化，仍要承担一定的滑坡风险[7];若进行基底处理，不仅需较大的剥离成本，而且当基底弱层较厚时，很难获得理想的治理效果[8]。留设支挡煤柱是改善软弱基底内排土场稳定性的重要手段之一，但必然造成一定的煤炭损失。因此，确定的合理支挡煤柱形态参数是解决该问题的关键。

传统的煤柱支挡效应研究仅考虑支挡煤柱与基底间抗剪力，应用二维刚体极限平衡法对支挡煤柱形态进行优化，使软弱基底内排土场的稳定性达到安全储备系数要求，没有考虑支挡煤柱在端帮两侧的抗剪力，忽略煤柱的三维支挡效应，易造成支挡煤柱形态参数过大，资源损失严重;现有的三维刚体极限平衡法通常将滑体近似为椭球体[9-11]，其滑坡模式与软弱基底内排土场边坡不符;数值模拟方法本身的复杂性难以在工程中推广[12]。因此，迫切需要提出一种新的科学手段，用于解决软弱基底内排土场支挡煤柱形态参数设计难题，对实现露天矿安全、经济开采有着十分重要的意义。

笔者在三维受力状态分析的基础上，推导获得了煤柱三维支挡效应的表达式;基于对煤柱三维支挡效应的力学成因类型及影响因素的认识，提出了煤柱支挡效应的二维等效方法，从理论上揭示支挡煤柱形态参数对支挡效应的影响规律，并以胜利东二露天煤矿软弱基底内排土场为工程背景，优化了软弱基底内排土场支挡煤柱的形态参数。

1 支挡煤柱的力学效应分析

1.1 力学模型的建立

受矿山采运设备规格与剥采工程需求限制，内排土场支档煤柱的空间形态可近似为具有梯形截面的四棱柱，从力学角度则可视为两端嵌入端帮的固定梁(图1(a))。显然，支挡煤柱受内排土场后方静止土压力F,侧界面正应力σ1，σ2,抗剪力T1，T2以及底界面支撑力N、抗剪力S与重力W的作用(图1(b))，煤柱能够保持稳定且发挥最大支挡效应的条件为F<S+T1+T2，即支挡效应的大小由T1，T2和 S的合力大小决定。其中，由于煤柱宽度与两侧端帮长度比值较小，可忽略两侧端帮三角形垂直载荷作用产生的附加力[13]，T1，T2的大小仅由煤柱的性质与侧界面形态参数决定，且T1=T2;底界面抗剪力S由煤柱基底岩层性质与底界面形态参数决定。对于特定的矿山，煤柱及其基底岩层性质可视为已知量，因此，煤柱的支挡效应仅与其空间形态参数β，ω，H，a，d有关。其中，β，ω分别为煤柱下底角，(°);H，a，d分别为煤柱高度、顶宽与长度，m。

图1 支档煤柱空间形态与受力分析模型
Fig.1 Space shape and force analysis model of supporting and retaining coal pillar

1.2 煤柱三维支挡效应解析

由于煤柱的支挡效应仅与其空间形态参数相关，结合煤柱自重应力的分布规律，以煤柱底界面一侧角点作为原点，将梯形断面划分为3个区域，命名为Ⅰ区、Ⅱ区、Ⅲ区，建立随 x1点与原点距离变化的煤柱高度函数HⅠ(x),HⅡ(x),HⅢ(x)，如图1(c)所示。

根据煤柱高度函数H(x)求出侧界面任意一点(x,y)的竖向自重应力σh(x):

然后，利用静止土压力计算方法[13]求出侧界面任意一点(x,y)的水平向的自重应力σ1(x):

最后，利用微积分原理，便可获得各区域侧界面抗剪力TⅠ,TⅡ,TⅢ，求和即获得支挡煤柱侧界面抗剪力T1。

同样，可基于图1(c)将煤柱分区求算得出底界面任意一点(x,y)的正应力σd(x):

式中，φj为煤柱底板岩层的内摩擦角，(°);cj为煤柱底板岩层的黏聚力，kPa。

最后，利用微积分原理，便可获得各区域底界面抗剪力SⅠ，SⅡ，SⅢ，求和为煤柱底界面抗剪力S。

2 煤柱支挡效应二维等效方法

目前，二维刚体极限平衡法是排土场边坡稳定性分析中最有效、应用最广泛的方法，为此，如何将获得的煤柱三维支挡效应进行二维等效，是能够将成果实现工程应用的关键。显然，煤柱的三维支挡效应为侧界面与底界面的抗剪力合力，而采用二维刚体极限平衡法时，煤柱的支挡效应仅体现为底界面的抗剪力，因此，只要找到合适的方法将侧界面上的抗剪力等效至底界面，即可实现三维支挡效应的二维等效。

分析式(17)可知，煤柱三维支挡效应实质是一剪切反力，表现为不同方位的黏聚力和内摩擦阻力的合力，分别由抗剪强度参数黏聚力与内摩擦角提供。前者仅与煤柱的形态参数有关，后者则可能由于排弃工程的发展而受到煤柱上方垂直外载荷的影响，这与二维刚体极限平衡法中条块底界面剪切反力的组成类型与影响因素是一致的。因此，可以通过合并同类项的方式，求解底界面的等效黏聚力cd及等效内摩擦角φd或等效内摩擦系数tan φd来实现煤柱支挡效应的二维等效。

分析式(19)可知，等效黏聚力cd及等效内摩擦系数tan φd与煤柱顶宽a、下底角β、ω及长度d均呈反比例函数关系，但与煤柱顶宽a、下底角β，ω呈正相关，与煤柱长度d呈负相关;当d=∞时，煤柱侧界面抗剪力可忽略，表现为二维支挡效应。需要注意的是，式(19)不能直接体现等效黏聚力cd及等效内摩擦系数tan φd随煤柱高度H的变化规律，可将式(19)转化为

分析式(20)可知，等效黏聚力cd及等效内摩擦系数tan φd与煤柱高度H呈倒差函数关系，单调递增。

综上分析可知，除煤柱长度d以外的各个煤柱形态参数值的增大均有利于提高其支挡效应。

3 工程实例

3.1 工程背景

胜利东二号露天煤矿内排土场基底倾角2°～3°，正常作业参数为平盘宽度50 m、坡面角33°、台阶高度24 m。岩性以泥岩为主，高岭石、蒙脱石等黏土类矿物成分较高，遇水软化严重，属典型的软弱基底内排土场。由于南帮滑体不断向内排土场呈近似流体状运动，828水平以下的内排空间几乎完全被滑体占据;另一方面，地下水的长期作用使靠近西帮的局部内排基底严重泥化，必将进一步降低内排土场稳定性。为节省滑体清理费用，继续稳定内排，支挡煤柱合理的形态参数是解决问题的关键。内排土场典型工程地质剖面和剥采工程现状如图2，3所示，各岩土体物理力学指标见表1。

图2 内排土场典型工程地质剖面
Fig.2 Typical engineering geologic profile of the inner dump

图3 胜利东二露天煤矿剥采工程平面
Fig.3 Mining and stripping engineering plan of Shengli East No.2 open-pit mine

表1 岩土体物理力学指标
Table 1 Geotechnical physical mechanics index

3.2 煤柱支挡效应影响因素分析

式(19)表明，影响煤柱支挡效应影响因素众多，笔者基于MATLAB与origin软件，采用控制变量法对单一因素进行研究[14]，揭示支挡煤柱形态参数对支挡效应的影响规律，分析结果如图4所示。

从图4可以看出，工程实例揭示的支挡煤柱形态参数对支挡效应的影响规律与理论分析具有一致性;支挡煤柱形态参数对等效黏聚力cd影响幅度较大，对等效内摩擦系数tan φd次之;支挡煤柱形态参数对等效黏聚力cd与等效内摩擦系数tan φd的敏感性大小顺序为H>d>a>β，ω，说明H对内排土场稳定性影响最大。

图4 煤柱支挡效应与形态参数的关系曲线
Fig.4 Curve of relationship between the supporting and retaining effect and the shape parameters

3.3 支挡煤柱形态参数确定

参照《煤炭工业露天矿设计规范》(GB50197―2015)[15]，并综合考虑煤炭资源损失量、边坡服务年限、重要程度、岩土体指标掌握程度以及潜在危害，确定胜利东二号露天煤矿内排土场边坡的安全储备系数为1.20。

根据剥采工程现状可知，下底角ω=28°;由于南北端帮位置固定，煤柱长度d=108 m;为满足剥采工程作业要求需在816、804水平留有15 m宽的运输平盘，台阶坡面角为65°，因此下底角β=36°。

考虑到内排空间的充分利用，内排土场按正常作业参数跟进。此时，对于任一给定的煤柱高度H，通过式(19)可获得不同顶宽a所对应的等效黏聚力cd及等效内摩擦系数tan φd，进而可通过内排土场稳定性二维分析确定稳定系数恰好满足安全储备系数的顶宽，即最佳顶宽a。笔者分别对是否考虑支挡煤柱侧界面抗剪力T1，T2两种条件下，确定了煤柱高度为28，30，32，34，36 m时的最佳顶宽，并计算了相应的煤柱体积V，见表2。列举H=36 m的内排土场稳定性二维分析结果如图6，7所示。考虑T1，T2时不同高度下稳定系数与顶宽关系曲线如图5所示。

表2 煤柱体积计算
Table 2 Calculation of coal pillar volume

图5 不同高度下稳定系数与顶宽关系曲线
Fig.5 Relationship curves between Fs and aunder different heights

图6 H=36 m时不同煤柱顶宽的稳定系数计算结果
(不考虑T 1,T 2)
Fig.6 Factor of stability results of different a when H=
36 m (not consider T 1,T 2)

图7 H=36 m时不同煤柱顶宽的稳定系数计算结果
(考虑T1,T2)
Fig.7 Factor of stability results of different a when H=36 m (consider T1, T2)

分析图5可知，考虑支挡煤柱侧界面抗剪力T1,T2时内排土场边坡稳定系数Fs随顶宽a增加而增加，与煤柱高度H呈线性递增关系;分析表2可知，在满足同样的安全储备系数前提下，是否考虑支挡煤柱侧界面抗剪力T1，T2两种条件下，支挡煤柱高度为28 m，顶宽分别为43，68 m时支挡煤柱体积最小，前者较后者能够使露天煤矿多回采煤炭资源约为75 600 m3的，按胜利东二露天煤矿当前吨煤纯利润70元计算，为露天煤矿创造近700万元的经济效益，由此可见，煤柱的三维支挡效应尤为重要，不能被忽略。

3.4 支挡煤柱空间破坏位置与稳定性数值模拟

为了验证支挡煤柱空间破坏位置与内排土场边坡稳定性，笔者基于支挡煤柱高度为28 m，顶宽为43 m的最优煤柱形态构建数值模拟模型(图8)。模型加载方式为自重加载，坡面为自由面，四周与底部分别采用水平和垂直位移约束。数值模拟计算结果显示:内排土场边坡整体稳定性系数Fs=1.208，与二维等效方法计算结果具有一致性。内排土场边坡三维位移云图如图9所示。

图8 数值模拟模型
Fig.8 Numerical simulation model

图9 内排土场边坡三维位移云图
Fig.9 3D displacement contours of internal waste dump

图10 支挡煤柱的空间破坏位置示意
Fig.10 Space damage position of supporting and retaining coal pillar

分析图9可知，支挡煤柱的空间破坏位置如图10中垂直面AB，DE，由此证明将支挡煤柱视为上下长度相等具有合理性;受排弃物侧向压力作用，临空面上部位移明显，向煤柱侧逐渐减小，说明煤柱具有较好的支挡作用。

4 结论

(1)将支挡煤柱视为在内排土场静止土压力及底界面剪切力双重作用下的固定梁，以摩尔-库伦准则作为失稳判据，推导了煤柱三维支挡效应的表达式，表明该效应体现为不同方位的黏聚力和内摩擦阻力的剪切合力，侧界面抗剪力是煤柱支挡效应的重要组成部分，不能被忽略。

(2)通过分析煤柱三维支挡效应的力学成因类型及影响因素，提出了煤柱支挡效应的二维等效方法，推导了底界面等效黏聚力cd与等效内摩擦角φd的表达式，可与经典的剩余推力法相结合，优化设计软弱基底内排土场支挡煤柱的形态参数;该方法计算的结果与数值模拟结果进行对比，表明该方法具有合理性。

(3)煤柱的等效黏聚力cd及等效内摩擦系数tan φd与其顶宽a、下底角β、ω均呈正相关反比例函数关系，与高度H呈正相关倒差函数关系，与长度d呈负相关反比例函数关系。因此，除煤柱长度以外，各个形态参数的增大均有利于提高其支档效应。

(4)工程实例分析表明，胜利东二露天煤矿内排土场边坡稳定系数Fs随顶宽a增加而增加，与煤柱高度H呈线性递增关系，当煤柱高度H=28 m，顶宽a=43 m时，既满足安全要求，又能实现经济效益最大化。

[1] 唐文亮,彭洪阁,马力,等.露天矿内排土场滑坡成因及治理措施研究[J].煤炭技术,2016,35(8):166-168.

TANG Wenliang,PENG Hongge,MA Li,et al.Study on cause and control measures of internal waste dump landslide in open-pit[J].Coal Technology,2016,35(8):166-168.

[2] 李江,张利忠,徐振东.宝日希勒露天矿内排土场边坡稳定分析[J].露天采矿技术,2009(3):29-31.

LI Jiang,ZHANG Lizhong,XU Zhendong.The slope stability analysis of internal dump in Baorixile open-pit mine[J].Opencast Mining Technology,2009(3):29-31.

[3] 鲁冰,陈吉望,孙凯,等.抚顺东露天矿内排土场长远发展规划[J].露天采矿技术,2015(12):15-18,25.

LU Bing,CHEN Jiwang,SUN Kai,et al.Long-term development plan of Fushun east open-pit mine[J].Opencast Mining Technology,2015(12):15-18,25.

[4] 杨秀.安家岭露天矿内排土场南帮边坡稳定性研究[D].阜新:辽宁工程技术大学,2012.

YANG Xiu.Study on the stability of south slope internal dump of Anjialing open-pit mine[D].Fuxin:Liaoning Technical University,2012.

[5] 王家臣,陈冲.软弱基底排土场边坡稳定性三维反演分析[J].中国矿业大学学报,2017,46(3):474-479.

WANG Jiachen,CHEN Chong.Three-dimensional back analysis for stability of slope dumped on weak basement[J].Journal of China University of Mining & Technology,2017,46(3):474-479.

[6] 杨巨文,曹兰柱,张利忠,等.露天矿软岩基底内排土场变形、破坏规律[J].辽宁工程技术大学学报,2007,26(S2):1-3.

YANG Juwen,CAO Lanzhu,ZHANG Lizhong,et al.Deforming and failing law of in pit soft rock basement dump in surface mine[J].Liaoning Technical University,2007,26(S2):1-3.

[7] 纪玉石.黄土基底排土场失稳机理与稳定控制技术研究[D].沈阳:东北大学,2013.

JI Yushi.The instability Mechanism and stability control technology Research on loess substrate Dump[D].Shenyang:Northeastern University,2013.

[8] 张建平,王俊.能量法在软弱基底排土场处理中的应用[J].煤炭技术,2015,34(6):204-205.

ZHANG Jianping,WANG Jun.Energy method application in mine weak basal processing[J].Coal Technology,2015,34(6):204-205.

[9] HUANG C C,TSAI C C.New method for 3D and asymmetric slope stability analysis[J].ASCE.Journal of Geotechnical and Environmental Engineering,2000,126(10):917-927.

[10] 陈祖煜,弥宏亮,汪小刚.边坡稳定三维分析的极限平衡方法[J].岩土工程学报,2001,23(5):525-529.

CHEN Zuyu,MI Hongliang,WANG Xiaogang.A three-dimensional limit equilibrium method for slope stability analysis[J].Chinese Journal of Geotechnical engineering,2001,23(5):524-529.

[11] 朱大勇,丁秀丽,刘华丽,等.对称边坡三维稳定性计算方法[J].岩石力学与工程学报,2007,26(1):22-27.

ZHU Dayong,DING Xiuli,LIU Huali,et al.Method of three-dimensional stability analysis of a symmetrical slope[J].Chinese Journal of Rock Mechanics and Engineering,2007,26(1):22-27.

[12] 王浩然,黄茂松,刘怡林.含软弱夹层边坡的三维稳定性极限分析[J].岩土力学,2013,34(S2):156-160.

WANG Haoran,HUANG Maosong,LIU Yilin.Three-dimensional stability analysis of slope with weak interlayer[J].Rock and Soil Mechanics,2013,34(S2):156-160.

[13] 邵全,韦敏才.土力学与基础工程[M].重庆:重庆大学出版社,1998.

[14] 周伟,韩流,舒继森,等.凹形边坡圆弧滑动模式稳定性控制机理研究[J].中国矿业大学学报,2016,45(1):70-76.

ZHOU Wei,HAN Liu,SHU Jisen,et al.Research on stability control mechanism of concave slope with circular sliding mode[J].Journal of China University of Mining & Technology,2016,45(1):70-76.