Roof stability of underhand slice stopping with cemented filling
-
摘要: 胶结充填体顶板的稳定性是决定下向分层胶结采矿法回采工作面安全的重要因素.以某矿山为工程背景，采用现场监测和数值计算相结合的手段，通过对其应力、位移、塑性区以及钢筋的受力特征进行分析，研究固化剂胶结尾砂充填体充当下分层回采顶板时的稳定性.结果表明：钢筋对充填体整体的稳定性起到加固作用，采用质量比为1:6和1:12的新型固化剂胶结分级尾砂材料明显优于质量比为1:4和1:8的42.5级水泥胶结分级尾砂材料；结合现场监测和数值计算结果，表明采用新型固化剂胶结分级尾砂充填，可以满足下分层回采顶板稳定性的要求.Abstract: The stability of roof is the main factor to determine the safety of mining face in the underhand slice stopping with cemented filling method. Taking a mine as engineering background, the stress, displacement, plastic zone and the force characteristics of steel were analyzed to study the roof stability of cemented classifying tailing by new cemented agent using the means of combining field monitoring with numerical methods. The result shows that the steel reinforcement has a great effect to backfill of the whole body; and when the mass ratio are 1 to 6 and 1 to 12, the new cemented agent materials are batter than the 42.5 grade cement agent material with the mass ratio of 1 to 4 and 1 to 8; according to the numerical calculation results and field monitoring, the new cemented classifying tailings can meet the requirements of filling stope roof stability.
-
Keywords:
- cemented classifying tailings /
- roof stability /
- stress /
- displacement /
- numerical methods
-
0 矿井通风系统可靠性分析的特点
矿井通风系统可靠性的高低对矿井安全生产有着直接的影响.国内外自20世纪80年代以来，已取得了初步的研究成果[1-4].但由于问题本身的复杂性，还有许多问题尚未解决.矿井通风系统可靠性研究应在充分考虑以下几个特点的基础上进行.
(1)矿井通风系统是由风机、分支和通风构筑物三者组成的有机整体，在讨论矿井通风系统可靠性时，应充分考虑各单元的特点以及各单元之间的相互联系和相互制约；
(2)各通风构筑物单元在使用过程中的可靠性是相对的、动态变化的，其在某一时刻的可靠度不仅依赖于当时构筑物所处的状态，也与预计的剩余服役期有关；
(3)各分支的风量在系统结构不变时，不仅与风阻和风压有关，而且与风机和通风构筑物当时所处的状态有关.
1 分支的可靠性
分支可靠性定义为：在矿井通风网络结构不变的条件下，在规定的时间内某一分支i的风量值能够保持在一个合理的范围之内的概率，称为该分支的可靠度，记为Ri.并相应地称该分支在任意时刻能保持其风量值在一个合理的范围之内的概率为该分支的有效度，记为Ai.
以Qi(t)表示实际情况下第i分支在t时刻的风量值，Qi(0)表示在风机和通风构筑物都稳定可靠的前提下第i分支的风量值.根据安全规程规定或分支用风需求确定两个风量限定值Wi(1)和Wi(2)，当第i分支满足
(1) 则第i分支正常，否则第i分支为不正常.因此第i分支的有效度Ai可表示为
(2) 通常，矿井通风网络的分支可靠度不仅是指风流的稳定程度或抗干扰能力，它与下列因素密切相关：①风流稳定性，角联网络分支造成的风流反向的可能性，风机的不稳定运行，井下提升运输设备因素等对风流运动的瞬时扰动现象；②调节设施工作的可靠性以及失效状态可能产生的危害程度和范围；③分支中的有毒性、窒息性和爆炸性气体与粉尘含量，一旦超限将使局部网络暂时失效；④漏风情况，风流的质量状况以及气候条件等.
由文献[5-8]可推导出，Qi(t)保持在闭区间[Qi(0)-Wi(1)，Qi(0)+Wi(2)]内的概率(即第i分支的可靠度)为
(3) 
(4) 
(5) Ri值即为在综合考虑了分支风阻的随机扰动及通风构筑物的可靠性等因素后，第i条分支的风量保持在[Qi(0)-Wi(1)，Qi(0)+Wi(2)]范围内的概率，即第i条分支的可靠度.
2 风机的可靠性
可用一般可修系统的可靠性计算方法来评价风机的可靠性[9-12].
对于单台风机若其寿命遵从指数分布1-e-λt，t≥0，故障后的修理时间也遵从指数分布1-e-μt，t≥0.当两个随机变量相互独立，且故障部件修复后的寿命分布与新部件相同时，单台风机的有效度为
(6) 对于主扇，根据矿山安全规程的要求，都需要设备用风机(一般只有电机备用).可将主扇单元看作用两个相同型号部件组成的冷储备系统，若这两台风机寿命和故障后修理时间均服从负指数分布，则主扇单元的有效度为
(7) 一般情况下，若某矿备有n台相同型号的主扇，正常情况下需要k台主扇工作，当工作的某台主扇发生故障后，则将剩余的n-k台中的某台启动，此时可将主扇单元视为k/n(G)表决系统，若每台主扇的寿命分布和故障后的修理时间分布均服从负指数分布，则主扇单元的有效度为
(8) 3 通风构筑物的可靠性
通风构筑物的类别不同，其可靠度的定义和计算方法亦不同.
3.1 截断风流的通风构筑物可靠性
对于风门、挡风墙、风帘截断风流的通风构筑物，其主要功能就是将风流截断，阻止风流通过.其可靠度定义为：对于截断风流的通风构筑物，在某一时刻t完成截断风流的能力，称为该构筑物在t时刻的可靠度.若Rs1(t)表示这类构筑物在t时刻的可靠度，则有
(9) 式(9)中，Qi(t)为截断风流的通风构筑物在t时刻的漏风量，m3/s；QSj(t)为截断风流的通风构筑物入风侧的设计风量，m3/s.其值为网络解算所得的风量值，此值不随构筑物的服役期t的变化而变化，而只与t时通风系统的网络结构有关.
新建截断风流的通风构筑物总有一定的漏风量，其值可通过实测而得，以Qi(0)表示.则新建截断风流的通风构筑物的可靠度为
(10) 该通风构筑物在使用过程中，随着其服役期t的增长，由于不断受到外部环境及构筑物本身的逐渐老化等因素的影响，在不维修的情况下，其漏风量会逐渐增加.在t时构筑物的漏风量Qi(t)可以按式(11)评估
(11) 式(11)中，Ci(t)是与构筑物在服役期内的经历有关的系数.
3.2 通过风流的通风构筑物可靠性
风桥等通过风流的通风构筑物的功能就是要使风流顺利地通过.可将其可靠度定义为：对于通过风流的通风构筑物在某一时刻t，其能使风流顺利通过的能力称为该构筑物在t时刻的可靠度.若以RS2(t)表示这类构筑物在t时刻的可靠度，则有
(12) 式(12)中，Qi(t)为通过风流的通风构筑物入风侧的风量，m3/s，即为网络解算所得之风量值；Q0(t)为该通风构筑物出风侧的风量值，m3/s.
新建构筑物在风流通过的过程中总有一定的损失，即有
(13) 即其风量损失为
(14) 在不维修的情况下，随着构筑物服役期t的增长，QS(t)逐渐增大.QS(t)可以实测，也可以按式(15)评估
(15) 式(15)中，CS(t)是与构筑物在使用期内的经历等因素有关的系数.
3.3 调节风流的通风构筑物可靠性
对于风窗等调节风流的通风构筑物，其可靠度定义为：对于调节风流的通风构筑物在某时刻t，其能通过的风量符合设计要求的程度，称为该构筑物的可靠度.并以RS3(t)表示.据此定义有
(16) 式(16)中，Qp(t)为t时刻该类构筑物的实际过风量，m3/s；Qd为设计构筑物的过风量，m3/s.Qp(t)值可以实测获得，也可利用与式(15)相似的公式评估而得.
4 矿井通风网络的可靠性
对于某一矿井通风网络，以m表示其节点个数，P表示其通路矩阵，A表示其关联矩阵，I为单位矩阵.若令C=A+I，取矩阵C中元素Cm1的余子式，并采用逻辑乘和逻辑加的运算方法对这一余子式展开，则得到从起始节点(进风井口)到终止节点(排风井口)的全部有向通路.
若以pij描述分支ej是否在第i通路上，并令
(17) 则可求得该矿井通风网络的通路矩阵为
(18) 式(18)中，k为该网络的通路数；n为该网络的分支数.
若以Rnet表示矿井通风网络的可靠度，则有
(19) 式(19)中，pi表示矿井通风网络的第i条通路.
式(19)是网络可靠度的计算式，它与网络中进风点和出风点间的最小通路密切相关.但由于最小通路间的不相斥性，具有大量公共边，使得网络可靠度的计算相当困难.为了改进算法，采用不交型布尔代数算法.通过不交型布尔代数的化简、展开吸收、归并等运算，使得最小通路间不含有公共边.这样对于通路，可按互不相容事件求概率，得到网络可靠度.
5 矿井通风系统可靠度优化分配
采用一个可行的状态S=(R，Q，H)，使得S在满足按需分风和经济性的条件下，可靠度达到最大，这个问题可用下列非线性规划进行描述[13].
(20) 式(20)中，S1和S2分别为固定风量分支集和风机分支集；bij为该网络回路矩阵B中的元素；Ei为第i台风机的目标流线分支集；Hfi为第i台风机的压力上限值.
式(20)的目标函数比较复杂，直接求解困难较大，应设法把目标函数进行简化.由式(4)和式(5)不难看出，要使得每个Ri都尽可能的大，只有让均值μ*i尽可能地接近Qi(0)，让方差σ*i尽可能小，而由μj(1)、μj(2)、σj(1)和σj(2)的不确定性，也只有让每个αij和βij的绝对值尽可能地小才行.根据已有研究结果[7-8]，矿井通风网络的阻力定律、风量平衡定律和风压平衡定律可以综合写成下面的方程组.
(21) 由隐函数存在定理知，如矿井通风网络处于平衡状态S0的Jacobi行列式不为零，即
(22) 则式(21)在S0的某一邻域内存在唯一的连续偏导数解[14].
(23) 由此可知，当J的绝对值充分大时，可使每个αij和βij的绝对值尽可能地小.式(20)的目标函数可以改为对J进行极大化.各分支风量一定时，调节各分支的风阻极大化J的优化模型如下：设已知各分支的风量为Q1(0)，Q2(0)，…，Qn(0)，至少存在一棵生成树，使其所有余树枝风量不为零.从经济方面考虑，约定风机的压力上限分别为Hf1，Hf2，…，HfF，则这时的优化模型为
(24) 式(24)中，ΔRj为分支j的风阻增量；J(ΔR1，ΔR2，…，ΔRn)为Jacobi行列式；S3为不可调分支[15]集；Ei和S2同前，M=106，且
(25) 
(26) 
(27) 
(28) 
(29) 6 应用
图 1为某矿井通风网络，10条分支，6个节点.风机工作方式为抽出式，型号为4-72-11№20，转速为315r/min，风流上行流动.各分支对应的风阻值(N·s2/m8)为式(30)，该风机的H-Q特性方程为式(31)，分支4、8、9为固定风量分支，该网络无火灾时的风量分配为式(32)
(30) 
(31) 
(32) 现分析火灾时各分支的风流稳定性.选一棵生成树，如图 2所示(实线部分).取Wi(1)=0.2Qi(0)，Wi(2)=0.3Qi(0)，μi(1)=μi(2)=0，σi(1)=0.04，σi(2)=1，i=1，2，…，10.上机计算求得各分支可靠度
(33) 由可靠度知分支5风流最不稳定，分支9风流最稳定；固定风量分支中，分支4最不稳定，分支9最稳定.
当分支7发生火灾时，主干分支中的风量增加，旁侧分支中的风量减小甚至反向.与主干分支相交的分支4、6、8、9中，由可靠度可知分支4风流首先反向，分支6次之，分支8再次之，这与分支可靠度是对应的.解算结果证明由分支可靠度来判断风流稳定性是正确的.尽管分支5没有反向，但它的风流稳定程度受干扰最大，例如当Δt=1300℃时，分支5中风量变化率为401.40%，其他分支中风量变化率值均小于它.实际上，风流稳定性是风向和风量稳定性的总称.
-
![]()
图 2 -310 m中段第3分层5’分条和5分条测点布置
Fig 2. -310 m Middle third hierarchical 5'points and 5 points of the measuring points laid
![]()
图 7 采用质量比为1:6和1:12固化剂胶结充填Y方向应力云图
Fig 7. Y direction stress cloud of cemented filling with mass ratio of 1:6 and 1:12 cemented agent
![]()
图 8 采用质量比为1:6和1:12固化剂胶结充填Y方向位移云图
Fig 8. Y direction displacement cloud of cemented filling with mass ratio are 1:6 and 1:12 cemented agent
![]()
图 9 采用质量比为1:6和1:12固化剂胶结充填塑性区云图
Fig 9. Plastic zone cloud of cemented filling with mass ratio of 1:6 and 1:12 cemented agent
表 1 第3分层5’分条应变计编号与布点编号对照
Table 1 Third hierarchical 5'points of strain gage ID number and distribution table
布点位置 测点方向 应变计编号 上部 沿进路方向 35’B01 垂直进路方向 35’B02 垂直作业面方向 35’B03 中部 沿进路方向 35’B04 垂直进路方向 35’B05 垂直作业面方向 35’B06 下部 沿进路方向 35’B07 垂直进路方向 2FB08 垂直作业面方向 2FB09 
下载: 导出CSV
表 2 第3分层5’分条钢筋计编号与布点编号对照
Table 2 Third hierarchical 5'points of reinforcement meter ID number and distribution table
布点位置 钢筋计编号 1号 2号 3号 4号 5号 6号 上部 35’G01 35'G04 35'G07 35'G10 35'G13 35'G16 中部 35’G02 35'G05 35'G08 35'G11 35'G14 35'G17 下部 35’G03 35'G06 35'G09 35'G12 35'G15 35'G18
下载: 导出CSV
表 3 钢筋力学参数
Table 3 Steel reinforced mechanical parameters
参数名称 直径
/mm长度
/m屈服点
/MPa抗拉强度
/MPa弹性模量
/MPa参数指标 16.0 2.4 335 490 200
下载: 导出CSV
表 4 岩体及充填介质力学参数
Table 4 Mechanical parameters of rock mass and filling medium
类型 度/(kN·m-3) 性模量/GPa 泊松比 度/MPa 结力/MPa /(°) 围岩 25.87 20.83 0.20 2.50 3.00 44 矿体 33.89 10.63 0.26 1.20 2.20 36 1:12固化剂充填体 19.50 0.45 0.29 0.40 0.50 17 1:6固化剂充填体 19.50 0.73 0.25 0.60 0.70 19 1:8水泥充填体 20.00 0.42 0.27 0.36 0.50 16 1:4水泥充填体 20.00 0.69 0.24 0.50 0.60 18
下载: 导出CSV
表 5 采用质量比为1:6和1:12固化剂胶结充填（加竖筋）
Table 5 Using mass ratio of 1:6 and 1:12 curing agent cemented agent filling and vertical reinforcement
开挖步骤 垂直方向应力/kPa 垂直方向位移/m 塑性区分布 1 200 0.004 5 无 2 500 0.008 0 无 3 500 0.008 0 侧帮充填矿柱有小块塑性区 4 500 0.010 0 侧帮充填矿柱有小块塑性区 5 600 0.010 0 侧帮充填矿柱有小块塑性区 6 600 0.010 0 侧帮充填矿柱有小块塑性区
下载: 导出CSV
-
[1] 李建波.金属矿床地下开采[M].北京:冶金工业出版社, 2011. [2] 周科平, 夏明, 肖雄, 等.开采环境再造深孔诱导崩落模型试验研究[J].岩石力学与工程学报, 2009, 28(11):2328-2335. doi: 10.3321/j.issn:1000-6915.2009.11.023 [3] 李一帆, 张建明, 邓飞, 等.深部采空区尾砂胶结充填体强度特性试验研究[J].岩土力学, 2005, 26(6):865-868. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200506007.htm [4] 谢开维, 何哲祥, 张常青.下向分层进路充填采矿法中采用尾砂胶结充填的研究[J].黄金, 2000(2):15-18. http://www.cnki.com.cn/Article/CJFDTOTAL-HJZZ200002004.htm [5] 刘贤俊, 何文, 谢勇, 等.新型尾砂胶结剂胶结分级尾砂充填体强度研究[J].有色金属科学与工程, 2015, 6(2):72-79. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201502014 [6] 何文, 高忠, 赵奎, 等.胶结充填体顶板裂缝雷达响应特征的正演研究[J].矿业研究与开发, 2015(1):51-55. http://www.cnki.com.cn/Article/CJFDTOTAL-KYYK201501012.htm [7] 刘志祥, 李夕兵.尾砂胶结充填体力学试验及损伤研究[J].金属矿山, 2004(11):22-24. doi: 10.3321/j.issn:1001-1250.2004.11.007 [8] 周瑞林.分级尾砂胶结充填在武山铜矿的应用[J].有色金属 (矿山部分), 2010, 62(3):1-2. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKU201003000.htm [9] 孙德民, 任建平, 焦华喆, 等.某矿全尾砂胶结充填物料性能研究[J].金属矿山, 2012, 41(3):6-9. http://www.cnki.com.cn/Article/CJFDTOTAL-ZXLJ201312055.htm [10] 李夕兵, 刘志祥.高阶段尾砂胶结充填体力学研究与博弈树配比优化[J].安全与环境学报, 2004, 4(4):87-90. http://www.cnki.com.cn/Article/CJFDTOTAL-AQHJ200404025.htm [11] 曹帅, 宋卫东, 薛改利, 等.考虑分层特性的尾砂胶结充填体强度折减试验研究[J].岩土力学, 2015, 36(10):2869-2876. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201510018.htm [12] 徐文彬, 杜建华, 宋卫东, 等.超细全尾砂材料胶凝成岩机理试验[J].岩土力学, 2013(8):2295-2302. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201308033.htm [13] 祝丽萍, 倪文, 黄迪, 等.粉煤灰全尾砂胶结充填料[J].北京科技大学学报, 2011, 33(10):1190-1196. http://www.cnki.com.cn/Article/CJFDTOTAL-BJKD201110002.htm [14] 卢平.制约胶结充填采矿法发展的若干充填体力学问题[J].黄金, 1994(7):18-22. http://www.cnki.com.cn/Article/CJFDTOTAL-HJZZ407.005.htm [15] 刘志祥, 兰明, 肖思友, 等.尾砂胶结充填体的损伤破坏及其与岩体的合理匹配[J].中国有色金属学报:英文版, 2015(3):954-959. http://www.cnki.com.cn/Article/CJFDTOTAL-ZYSY201503033.htm [16] 赵奎, 王晓军, 刘洪兴, 等.布筋尾砂胶结充填体顶板力学性状试验研究[J].岩土力学, 2011, 32(1):9-14. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201101004.htm [17] 王健, 高卫红, 龚囱, 等.下向进路胶结充填顶板稳定性数值模拟[J].江西有色金属, 2008, 22(4):4-7. http://www.cnki.com.cn/Article/CJFDTOTAL-JXYS200804003.htm [18] МАЛЕТИН, 陈鼎懿.胶结充填体的稳定性[J].有色矿山, 1987(10):24-25. http://www.cnki.com.cn/Article/CJFDTOTAL-KYYK404.008.htm [19] 董璐, 高谦, 南世卿, 等.超细全尾砂新型胶结充填料水化机理与性能[J].中南大学学报 (自然科学版), 2013, 44(4):1571-1577. http://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201304041.htm [20] 惠林.冬瓜山铜矿新型胶结材料的研究[J].金属矿山, 2008(9):159-160. http://www.cnki.com.cn/Article/CJFDTOTAL-JSKS200809052.htm
计量
- 文章访问数: 84
- HTML全文浏览量: 14
- PDF下载量: 1
