_11_6_特别重大坍塌事故矿区采场稳定性三维数值模拟的分析.docVIP

河北邢台石膏矿11.6特别重大坍塌事故 矿区采场稳定性三维数值模拟分析 马海涛 （中国安全生产科学研究院） 摘 要 本文采用岩土力学数值计算有限差分方法，建立了河北邢台11.6矿难塌陷区采场的三维数学模型，反演了该采区的开采历史和围岩破坏、围岩应力分布、地表沉陷、监测点位移变化等情况，从力学角度阐述本次事故发生发展的过程，对地下采场的安全稳定性给出了客观的评价。为事故处理责任的认定和避免同类事故的发生提出了参考性意见。 关键词 石膏矿 采场稳定性 数值分析 坍陷 3D- Numerical Simulation of Stope Stability in 11.6 Accident of Gypsum Mine Collapse in Xingtai Hebei Province Ma Haitao (China Academy of Safety Science And Technology) Abstract：finite difference method in rock-geotechnical mechanics, 3D numerical simulation model of gypsum mine in collapse field has been set up to judge the stability of underground mining area. Mining history, damage of wall rock, stress distribution, settlement of ground and displacement of monitor points have been studied within the numerical model. The results would give a reference for the responsibility confirmation, and make it possible to avoid some similar accidents. Keyword：cave-in 前言 2005年11月6日河北省邢台市境内的邢台县会宁镇尚汪庄康立石膏矿发生坍塌事故波及太行、林旺两个石膏矿塌陷区直径约米，米×米范围内地面不同程度出现裂缝，地面一幢二层楼坍塌造成33人死亡特别重大安全生产事故模型侧面限制水平移动，模型底面限制垂直移动。计算中采用莫尔-库仑(Mohr-Coulomb)屈服准则判断岩体的破坏 (1) 式中，(1、(3 分别是最大和最小主应力，c，( 分别是粘结力和摩擦角。当fs0时，材料将发生剪切破坏。在通常应力状态下，体的抗拉强度很低，因此可根据抗拉强度准则（(3((T）判断岩体是否产生拉破坏。 图1塌陷矿区采掘工程平面图 图2 三维数学模型网格图 图3 地层模型图 表1 计算参数 岩层 弹性模量 泊松比 体积模量 剪切模量 密度 内聚力 内摩擦角 抗拉强度 单位 MPa — MPa MPa Kg/m3 MPa 度 MPa 表土 50.00 0.30 41.67 19.23 1800 0.08 20.00 0.05 第四系沉积岩 800.00 0.23 493.83 325.20 2100 0.32 21.50 0.26 灰岩 50000.00 0.25 33333.33 20000.00 2400 8.40 40.00 4.70 石膏 70000.00 0.30 58333.33 26923.08 2360 5.40 38.60 2.80 矿井隔离矿柱稳定分析 设置矿井隔离矿柱的目的，是将矿井划分成独立的生产单元，使矿井在生产过程中尽量避免受相邻矿井开采活动的影响（包括开采扰动、岩体塌陷、通风、排水、运输和提升、水火灾害控制等）。因此，在矿山行业规范中又将其称之为“保安矿柱”[4]。 矿体开采后，上覆岩体重量向矿柱上转移，造成矿柱中产生应力集中，矿柱的边缘部分由于处于侧向无约束的单向受力（(3=0）状态产生破坏而丧失承载能力，此范围称为屈服区（塑性区），屈服区向内为弹性区（核区）。大量研究表明，欲使矿柱保持稳定，矿柱须有弹性核区的存在（图4）。 图4 矿柱受力状态分区 从图1的井田规划中可看出，除了塌陷核心区的康立与林旺矿之间没有设立保安矿柱外，其它各矿间均规划了宽度在20米以上的矿柱。计算结果证明，20米矿柱可以保证矿区的安全生产。根据最终计算的破坏区域来看（图5），由于事故矿区对井田保安矿柱进行超采，保安矿柱尺寸远远小于20米，大部分进入塑性屈服状态，已不具备稳定承载能力。由于第二石膏矿与事故矿区间的保安矿柱保存较为