2. 围岩力学性质与塑性区计算参数获取
塑性区计算的可靠性,首先取决于输入参数是否真实反映巷道围岩的力学状态。围岩塑性区半径、形态与分布范围,本质上由三组参数共同控制:一是岩石强度参数(黏聚力 c、内摩擦角 φ、单轴抗压强度 σc、抗拉强度 σt),二是变形参数(弹性模量 E、泊松比 ν),三是地应力参数(垂直应力 σv、水平主应力 σH、σh 及其方向)。本节围绕这三组参数的测试方法、取值原则与工程处理展开,重点说明如何把“试验室或现场得到的原始数据”转化为“可直接代入塑性区计算模型的参数”。
一、岩石强度参数:从岩块试验到岩体强度
强度参数通常通过室内岩石力学试验获得。单轴压缩试验可得到 σc 与 E、ν;三轴压缩试验在围压 σ3 分级加载下得到一组(σ1, σ3)破坏点,按莫尔–库仑准则拟合即可求黏聚力与内摩擦角。以常见的砂岩为例(示例数据,仅作教学演示):围压取 0、5、10、20 MPa 四级,对应峰值轴压约 60、95、125、180 MPa,在 σ–τ 平面上作破坏包络线,可得 c≈12 MPa、φ≈35°。若试样离散性大,应至少取 3~5 个试样,剔除明显异常点后再拟合,并给出参数区间而非单一值。
关键取值原则是:试验室测的是岩块强度,塑性区计算需要的是岩体强度。岩体被结构面切割,其强度显著低于岩块,必须进行折减。工程上常用地质强度指标 GSI 或 RMR 折减:由 σc 与 GSI 估算岩体强度,或按 Hoek–Brown 准则换算等效 c、φ。对于节理发育的软岩,折减系数可能低至 0.2~0.4;完整硬岩则可取 0.6~0.8。折减后还应结合现场变形、片帮、底鼓等显现反演校核,避免“参数偏强导致塑性区算小、支护偏弱”的常见失误。
二、变形参数:加载路径与尺度效应
弹性模量 E 应取应力–应变曲线中近似直线段的割线或切线模量,并注明取值区间(如 50% σc 处割线模量)。泊松比 ν 由横向与轴向应变之比确定,软岩中常达 0.3~0.4,硬岩约 0.2~0.25。需注意,塑性区计算中若采用弹塑性模型,E 影响的是弹性区应力重分布与位移量级,而 c、φ 决定塑性区边界,二者不可混用。对于层状或各向异性明显的围岩,应分别给出平行与垂直层理方向的模量,并在模型中按实际方位赋值。
三、地应力参数:原位测试与方向判定
地应力无法由岩块试验获得,必须依靠原位测试或区域资料反演。常用方法包括水压致裂法(可同时给出主应力大小与方向,适合深孔)、套芯应力解除法(精度高但工序复杂)、以及声发射凯泽效应估算法。垂直应力一般按上覆岩层自重估算 σv=γH,水平应力则通过侧压系数 k=σh/σv 表达,构造发育区 k 常大于 1,甚至达 2~3。方向判定至关重要:最大水平主应力方向与巷道轴线的夹角直接决定塑性区是否偏帮、顶板是否对称破坏,因此测试报告中必须给出主应力方位角,并在计算模型中建立一致的坐标系(规定 x 轴指向最大水平主应力方向,或明确标注方位)。
四、取值原则与工程处理要点
第一,分层分带取值。巷道穿过多种岩层时,应按顶板、两帮、底板分别给出参数,不能用一个“综合值”代替全部围岩。第二,给出区间并做敏感性分析。塑性区半径对 c、φ 和侧压系数 k 十分敏感,建议对关键参数取低、中、高三档分别计算,观察塑性区半径变化幅度,若变化超过 30%,说明参数不确定性大,应补充测试或加强支护冗余。第三,统一坐标系与符号约定。计算前明确压应力为正、x 轴指向最大水平主应力方向,并在报告中注明,避免方向混淆导致塑性区形态判断错误。第四,用现场显现反演校核。若计算塑性区远小于实际片帮、底鼓范围,应检查强度折减是否不足、地应力是否偏低,必要时按实测变形反演修正参数。
综上,参数获取不是简单的“做试验、填数字”,而是一条从岩块试验、岩体折减、地应力测试到模型赋值与反演校核的完整链条。只有把强度、变形、地应力三组参数放在同一坐标系和同一尺度下处理,塑性区计算结果才具备工程可信度,才能为支护设计提供可靠依据。