深部软岩巷道开挖变形影响因素的数值分析.docVIP

深部软岩巷道开挖变形影响因素的 数值分析 许广 唐又驰 （辽宁工程技术大学资源与环境工程学院 辽宁 阜新123000） 摘要: 根据深部软岩巷道变形力学机理，采用不连续变形分析（DDA）软件对影响巷道开挖变形的岩体力学因素，包括地应力、围岩破裂面倾角、围岩块体弹性模量、围岩块体泊松比、围岩破裂面粘结力、围岩破裂面内摩擦角、围岩破裂面抗拉强度变化引起围岩变形位移的规律进行数值模拟研究分析。结果表明地应力是影响围岩变形的最重要因素，影响围岩变形的诸多因素指标存在临界值,得出了影响围岩变形位移的围岩特性因素重要性顺序，提出深部软岩巷道围岩控制的关键在于及时加固围岩以提高其残余强度，优化支护形式，在围岩中形成关键承载结构以阻止变形区的扩展。 中图分类号：TD353 文献标识码：A 关键词： 深部软岩 ；变形； DDA ；位移 前言 煤炭在我国国民经济中占有举足轻重的作用，约占能源一次生产和消费结构的70%左右。我国煤炭资源丰富，己探明的煤炭资源量占世界总量的11.1％，但我国的煤炭资源埋深在1000m以下的为2.95×1012t，占煤炭资源总量的53％。随着大规模矿山的开采，我国煤矿开采深度以每年8-12米的速度增加，预计未来20年我国很多煤矿将进入1000-1500的开采深度。随着开采深度的增加，各种工程灾害问题也日益增多，其中深部软岩巷道围岩变形破坏现象尤为突出。矿井进入深部开采后，在浅部呈现中硬岩特征的岩体转化为高应力软岩，这些巷道共同的特点是在各种因素（如地应力、动压影响、地质构造及岩体成份）等影响下围岩节理裂隙发育，松散破碎，矿压显现强烈，变形量大、破坏范围大，呈持续流变形态，巷道冒顶、片帮和底膨层出不穷，严重制约矿井安全生产。为有效防治深部软岩巷道变形带来的灾害问题，众多学者对深部巷道变形问题进行了研究[1-8]。但是由于岩体本身的复杂性，许多问题还有待进一步研究。本文则在前人研究成果的基础上依据巷道围岩变形的非连续特性，应用非连续变形分析（DDA）软件[9，10]，对影响深部高应力软岩巷道开挖变形的岩体力学因素进行数值模拟分析，以期为工程实践提供理论支持。 1 深部软岩巷道开挖变形力学机理 巷道开挖前，地下岩体处于三轴压力平衡状态；一旦开挖，这个平衡系统就会遭到破坏，围岩应力会重新调整[9]。当巷道围岩应力小于围岩强度，则围岩处于自行稳定状态；当巷道围岩应力大于围岩强度，围岩则发生变形破坏，甚至处于不稳定状态或失稳状态[11]。由于深部巷道周边调整后的围岩应力完全达到了围岩强度，所以围岩发生破坏，导致围岩应力降低，而最大主应力向围岩内部转移，直至达到新的三向应力平衡状态为止，此时围岩中将出现一个破裂区，这个由于应力作用产生的环状破裂区称为巷道围岩松动圈[12]。它有一个发生、发展和稳定的过程，稳定后的松动圈厚度综合反映了围岩应力、围岩强度等因素共同作用的结果。松动圈以外是塑性区和弹性区，该“三区”的力学行为与岩石全应力应变曲线中的三个阶段相对应，其中巷道围岩弹性区和塑性区对应于全应力应变曲线峰前弹性段塑性段，破裂区对应于峰值后残余强度段。大量研究结果揭示，巷道支护的主要对象是松动圈形成发展过程中的碎胀变形压力。松动圈较小，围岩碎胀变形也email:1075348477@ 图1 巷道围岩三区与岩石全应力应变曲线 的对应关系 Fig.1 Corresponding relation between the three zones of surrounding rock and curve of stress vs strain Ⅰ－松动圈(残余强度区)；Ⅱ－松动圈(软化区)；Ⅲ－塑性区；Ⅳ－弹性区 2 DDA软件的计算优势 DDA（Discontinuous Deformation Analysis具有更加严密的数学依据和更加完备的理论基础。m，高3m；模拟巷道为可缩量大的整体式支护。地应力加载方式采用改变初始应力和集中荷载的方式，模拟巷道开采深度的变化。开挖巷道为先加载后开挖，并假定块体破裂面相互（节理面与层理面）呈90°夹角。采用平面应变模型计算，模拟范围为长×宽=15 m×15m，边界条件如图2所示。在巷道顶、底板及两帮布置4个围岩表面位移观测点（见图2，其中1、3点为上下帮观测点，2、4点为顶底板观测点）。选用数值模拟参数见表1。模拟方案为通过分别改变地应力、围岩破裂面倾角、围岩块体弹性模量、围岩块体泊松比、围岩破裂面粘结力、围岩破裂面内摩擦角、围岩破裂面抗拉强度的参数值来模拟巷道围岩的变形规律。 图2 模型边界条件 Fig.2 Boundary conditions of model 表1 数值模拟围岩参数表 Table 1 Simulation parameters of