(三)一矿防水围岩厚度综合计算0811.doc

平煤一矿系列资料之（三） 平煤一矿北三回风井（573-605m） 防水围岩安全厚度综合计算报告 中 国 矿 业 大 学 2013．08．12 目 录 1 基本情况介绍 2 2 井筒松动圈厚度计算 2 2.1 理论计算 2 2.2 数值模拟方法 3 3 防水围岩厚度计算 5 3.1 自推导公式计算 5 3.2 维亚洛夫—扎列茨基公式计算 5 3.3 拉麦公式计算 6 4 结 论 6 1 基本情况介绍 一矿北三回风立井由平煤建工集团建井一处负责施工，井筒设计深度1075m，净径6.5m，采用冻结法施工，设计冻结深度667m，井筒已施工至726.5m位置。井筒591m处井壁在壁后注浆的过程中压裂形成直径0.6m的孔洞，造成大量涌水涌砂现象，后采用构筑止浆垫，地面割管注浆治理，目前井筒排水清淤已经完成，正在进行割管注浆效果检验及补充注浆工作。 关于检查、补注范围的确定，需要计算防水围岩的安全厚度。根据工程设计原则，只要在安全厚度范围内，围岩不存在注浆盲区，即已形成防水帷幕，则认为注浆达到预期效果，可以进行揭止浆垫的工作，当然设计时需适当考虑一定的安全系数。 通过施工钻孔进行注浆质量的监测，本质上应是对防水围岩厚度范围内的注浆效果进行验证和补注，这样一来计算防水围岩的安全厚度是进行探水钻孔布设的前提。 2 井筒松动圈厚度计算 2.1 理论计算 图1 塑性区计算模型 在岩体中开挖地下硐室时，将在围岩中出现一个塑性松动圈。松动圈计算参考下式进行。 式中：R1——松动圈半径，m；R0——硐室半径，m； σ0——原始应力，MPa；Pi——支护强度，MPa；φ——岩石内摩擦角，°；Cm——岩石粘聚力，MPa；L——松动圈厚度，m。 参数取值：硐室半径R0取荒径4.05m，原始应力σ0=λrh=0.3*25000*600=4.5MPa，支护强度Pi取0，岩石内摩擦角φ取30°，岩石粘聚力Cm取1.5MPa，将上述取值代入公式： 经过理论计算，井筒600m位置围岩松动圈厚度0.68m。 2.2 数值模拟方法 数值计算采用模拟软件ADIAN进行，ADIAN模拟软件是基于有限元方法开发的，基本原理是依据地质模型建立数值模型，然后对数值模型模块进行赋值，然后软件自行计算，通过计算结果的提取分析得出所需的结果。 本工程所建立的数值计算模型如图2所示：模型为二维模型，长×宽=50×50，井筒荒径预设为4.05m，围岩基本力学参数：弹性模量E=3×108，内摩擦角ψ=30°，粘聚力c=1.5×106MPa。 图2 井筒塑性区数值计算模型 (a)应变云图 (b)塑性区发育云图 图3 数值计算结果 计算结果如图3所示，截取了井筒开挖后周围应变大小及塑性区发育范围两幅图，从图中可以看出，井筒开挖后，井筒近距离围岩的应变明显呈现突变现象，塑性区发育特征明显，从图3(b)的统计来看，塑性区厚度为0.7m左右，与理论计算结果较为吻合。 根据理论计算井筒600m处围岩松动圈厚度为0.68m，数值计算松动圈厚度为0.7m，两者综合取围岩松动圈厚度为0.7m。 3 防水围岩厚度计算 3.1 自推导公式计算 图4 理论计算模型 理论计算模型如图4所示，假设安全围岩厚度为L，围岩水压5.5MPa，井筒荒径4.05m，松动圈厚度0.7m，则计算安全围岩层内径为4.75m。在计算围岩对称取两个微计算单元A、A’，则其所受的纵向应力σ=(4.75*2+L*2)*5.5/(2L)。根据岩石抗压强度理论，在其安全保护层范围内微单元所受的应力不得大于岩石抗压强度σp，则纵向应力σ=(4.75*2+L*2)*5.5/(2L)σp，取壁后围岩抗压强度为35MPa，则(4.75*2+L*2)*5.5/(2L)35，L2.9m。 3.2 维亚洛夫—扎列茨基公式计算 维亚洛夫—扎列茨基公式为： 式中：E—安全防水围岩厚度，m；K—安全系数，取1.5；ζ—固约系数，取0.2；h—模板段高，取4.5m；P—水压，P=5.5Mpa；σ(t)—围岩长时强度，σ(t)==17.5Mpa。 经计算，600m段安全防水围岩厚度取2.93m。 3.3 拉麦公式计算 拉麦公式为： =3.00 m 式中a—井筒掘进半径，此取值应加上松动圈厚0.7m，a=4.05+0.7=4.75m；P—水压，P=5.5Mpa；σ-围岩许用抗压强度，取35/2=17.5Mpa。 经计算，600m段安全防水围岩厚度取3.00m。 4 结 论 通过两种方法确定的围岩松动圈厚度为0.7m。通过三种方法分别算的安