不同侧压系数下矩形巷道逐步开挖的稳定性分析.docxVIP

? ? 不同侧压系数下矩形巷道逐步开挖的稳定性分析 ? ? 栗东平，冯瑞广 （1.河北工程大学 土木工程学院，河北 邯郸 056038；2.河北工程大学 计算力学与工程应用研究中心，河北 邯郸 056038） 深部煤层处于高地应力、高地温、高孔隙水压力的地质环境。因为对深部巷道围岩的力学特性认识不够充分，不能制定合理的巷道围岩支护措施，导致层裂、岩爆、设备损坏和生产受阻等问题频频发生[1-4]。目前，许多学者对深部巷道围岩的稳定性进行了大量的研究。如：郭建卿等[5]通过FLAC3D软件建模分析侧压力系数对巷道变形及周边应力分布的影响规律；林志斌等[6]基于深埋圆形巷道工程的研究，得出了岩体总裂纹数与侧压力系数呈抛物线关系；范磊等[7]基于Mohr-Coulomb 和Drucker-Prager 强度准则，并采用数值模拟软件分析得出：随着侧压力系数的变化，围岩应力的分布状态基本不变；龙翼等[8]采用试验与数值模拟相结合的方法，对不同地应力场下爆破扰动诱发巷道破裂现象进行研究，得出了地应力场控制围岩破坏特征，爆破扰动加剧围岩破坏速度的结论；董春亮等[9]、王新丰等[10-11]基于数值模拟、层次分析及相似模拟试验等方法，得出了深部开挖巷道围岩变形破坏特征；王学滨等[12-14]基于三维非均质模型对巷道围岩进行研究，得出不同轴压比、内摩擦角对分区破裂化的影响规律；Chen 等[15]、Wei等[16]根据现场实测资料建立颗粒流模型，研究了侧压力对巷道围岩稳定性的影响规律及围岩的变形损伤特征和机理；Yu 等[17]、Meng 等[18]通过数值模拟，研究了深部巷道的应力、位移、和塑性区分布的演化规律。上述研究成果表明，深部岩体地应力场对巷道围岩的破坏特征的研究，有助于煤层深部开采围岩稳定性分析及安全控制措施的提出。巷道开挖扰动是围岩变形和破坏的主要因素，目前对巷道逐步开挖和不同侧压力对围岩稳定性的研究较少，巷道逐步开挖对围岩稳定性的影响机理有待于进一步研究。为此，根据某矿实际资料建立数值模型，采用快速应力边界法建立初始应力场[19]，之后进行巷道逐步开挖分析，并与相似模拟试验[20]结果进行对比分析，验证数值模型的可行性；最后通过修改模型的侧压力系数，对不同侧压系数下巷道逐步开挖塑性区和主应力差的演化规律进行研究。 1 数值模型的建立 1.1 参数选取 根据某矿轨道下山岩层数据确定深部巷道开挖模型参数[20]。建立数值模型尺寸为长（x 方向）×宽（y方向）×高（z 方向）：50 m×50 m×47 m。模型共被离散为82 000 个单元和87 822 个节点。初始边界条件为xoz 平面、yoz 平面和xoy 平面的下表面采用速度约束（固定速度为0 m/s），上边界为自由平面。在模型上边界处添加上覆岩层的自重应力，约为25.23 MPa。模型本构关系选为Mohr-Coulomb 强度准则，数值模型如图1。巷道形状为矩形，截面尺寸宽（x 方向）×高（z 方向）：5 m×4 m。岩层物理力学参数见表1。 图1 数值模型Fig.1 Numerical model 表1 岩层物理力学参数Table 1 Physical and mechanical parameters of rock strata 1.2 模拟方案 考虑到模型具有对称的特点，只对巷道的一侧进行监测。因此在模型y=25 m 平面上布置测点，主要布置点在巷道顶板、帮部和底板，各测线监测点间距1 m，测点布置如图2。计算过程分为2 步。 图2 测点布置图Fig.2 Layout of measuring points 1）根据实测数据建立模型并划分网格，对不同岩层进行赋值，采用S-B 法生成初始地应力场。 2）将模型初始速度、初始位移设置为0，对模型xoy 下边界、xoz 边界和yoz 边界添加速度约束，之后对巷道进行逐步开挖，每步开挖5 m，共计开挖10步，其中每2 次开挖间隔进行1 次静力平衡。 开挖10 次之后，形成贯穿模型的巷道。 模型计算中共采取5 个计算方案，侧压力系数分别取值为0.5、0.8、1.1、1.4、1.7，其他参数取值相同。 2 计算结果验证 文献［13］进行了相似模拟试验，试验原型为某矿轨道下山，巷道埋深1 000 m，掘进断面为高4 m宽5 m 的矩形截面，模型上边界处施加25.23 MPa的垂直应力，水平边界处施加27.75 MPa 的水平应力，其他条件详见上文参数。为与试验结果进行对比验证，在模型中依据试验测点位置，等比例布置监测点。数值模拟结果与文献中试验结果数据对比见表2 和表3。 表2 最大位移对比Table 2 Comparison of maximum displacement 表3 部分测点应力对比Table 3 Stress c