2.3 岩石力学与稳定性分析

岩石力学对许多矿业工程师来说比较晦涩：公式多、理论体系复杂，同一个岩体在不同模型中还可能使用不同参数。学习本章不要求把重点放在本构模型推导或复杂手算上，也不必为了记住每一个经验公式而纠结。矿业工程师更需要建立一种工程判断能力：看得懂岩石和结构面，知道什么条件可能失稳，知道该收集哪些数据，并能把分析结论转化为可实施的采矿方案。

阅读时可以始终围绕四个问题：

1. 岩石和岩体有什么性质？ 关注强度、变形、节理、断层、风化、含水、方向性和应力，而不只看一个 UCS 或 RMR 数字。
2. 工程可能怎样破坏？ 是台阶楔体、整体边坡滑动、巷道掉块、高应力片帮、采场失稳，还是矿柱超载？破坏模式决定分析方法。
3. 工程方案怎样控制风险？ 台阶角、跨度、采矿顺序、排水、爆破、支护、充填和人员暴露，往往比追求更复杂的计算更直接。
4. 现场表现是否支持设计？ 用揭露、位移、孔压、支护载荷和微震验证预测；若实际行为不同，就更新模型和方案。

因此，本章中的公式主要用于理解参数关系、检查数量级和阅读设计报告。真正需要掌握的是“岩石性质 → 破坏模式 → 分析方法 → 工程措施 → 监测反馈”这条主线。专业岩石力学工程师负责建立和签发模型，矿业工程师则需要提供可靠的地质与生产条件、识别不合理假设，并把稳定性要求落实到采场布置和日常作业中。

岩石力学设计是一条完整的证据链：工程地质模型 → 试验与结构面数据 → 岩体参数 → 破坏模式 → 分析模型 → 设计措施 → 监测验证。其中任一环节失真，复杂软件也只会给出精确但不可靠的结果。

使用边界：本章用于学习参数推导、模型选择和结果复核方法，不能替代注册或具备相应资质的岩土/采矿工程师签发的设计。最终安全系数、失效概率、支护参数和监测触发值，应服从项目所在地现行法规、批准的设计准则及现场风险评估。

2.3.1 岩石强度试验

学习目标

• 按工程问题选择 UCS、三轴、巴西劈裂、点荷载、直剪与变形试验
• 区分完整岩、结构面和岩体参数，避免把试件强度直接用于工程尺度
• 由三轴数据拟合 Mohr–Coulomb 或 Hoek–Brown 包络，并报告适用应力范围
• 建立取样、制样、试验、统计和设计值选取的 QA/QC 链条

核心概念

试验对象与设计用途

试验计划与样本代表性

试验数量不能只按“每种岩性几件”安排。至少同时覆盖：主要岩性、风化等级、蚀变带、矿岩接触带、软弱夹层、含水状态、不同取芯方向和预期应力区间。样本位置应保留钻孔号、深度、方位、回收率、RQD、照片和破坏模式，使结果能回到三维地质模型。

公式

• UCS 点荷载估算：$UCS \approx k \times Is(50)$（k 需本地标定）
• Mohr–Coulomb 剪切强度：$τ_{\mathrm{f}} = c + \sigma_{\mathrm{n}}' \times \tan\phi$
• Mohr–Coulomb 主应力式：$\sigma_{\mathrm{1}}' = A\sigma_{\mathrm{3}}' + B$，其中 $A = (1 + \sin\phi )/(1 - \sin\phi )$，$B = 2c \cos\phi /(1 - \sin\phi )$
• 完整岩 Hoek–Brown：$\sigma_{\mathrm{1}}' = \sigma_{\mathrm{3}}' + \sigma_{\mathrm{ci}} \times (m_{\mathrm{i}} \times \sigma_{\mathrm{3}}'/\sigma_{\mathrm{ci}} + 1)^0.5$
• 巴西劈裂抗拉强度：$\sigma_{\mathrm{t}} = 2P / (\pi Dt)$（P 为破坏载荷，D 为圆盘直径，t 为厚度）
• 割线模量：$E_{\mathrm{sec}} = \Delta \sigma / \Delta ε$；报告所取应力区间，避免不同实验室定义不一致

经验关系如 $BTS \approx UCS/8-UCS/12$ 或 $UCS \approx kIs(50)$只能用于筛查。k、BTS/UCS 和 E/UCS 比值均受岩性、含水、结构和试验方法影响，设计前应由同源配对样本建立本地回归及预测区间。

工具与标准

• 常用软件/仪器：伺服压力/三轴试验机、结构面直剪仪、点荷载仪、施密特锤、超声波测试；RocData 或统计软件拟合强度包络
• 相关标准：现行 ISRM 建议方法、GB/T 50266《工程岩体试验方法标准》及项目实验室 QA/QC 程序

操作步骤

1. 从决策反推试验：明确参数将用于台阶、总体边坡、巷道、采场、矿柱还是数值模型。
2. 建立地工域：按岩性、结构、风化、蚀变和含水条件划分样本组，不按钻孔平均。
3. 制定取样矩阵：覆盖最弱域、方向性与预期应力范围；记录无法制样的弱岩比例。
4. 试样验收：量测尺寸、端面、含水、密度、波速和可见缺陷；超差样本保留但单独标识。
5. 执行试验：记录完整应力—应变曲线、峰后响应和破坏照片，而非只抄最大载荷。
6. 拟合参数：比较 Mohr–Coulomb 与 Hoek–Brown，检查残差、异常值和适用围压区间。
7. 选取设计值：按后果、样本量和变异性选特征值/保守分位数，并做上下界敏感性。
8. 回写模型：结果关联至地工域、深度、方向和版本；新揭露或监测异常触发更新。

岩石强度试验方案流程

试验方法选择

可复算示例：由三轴主应力估算 c、φ

某完整岩在两个围压下的示例破坏点为 $(\sigma_{\mathrm{3}}, \sigma_{\mathrm{1}}) = (5, 38) MPa$ 与 (15, 70) MPa。仅用于说明计算，正式设计应采用更多围压点和重复样。

1. 主应力直线斜率：$A = (70 - 38)/(15 - 5) = 3.2$。
2. $\sin\phi = (A - 1)/(A + 1) = 2.2/4.2 = 0.524$，故 $\phi \approx 31.6^{\circ}$。
3. 截距：$B = 38 - 3.2 \times 5 = 22 MPa$。
4. $c = B(1 - \sin\phi )/(2cos\phi ) \approx 6.15 MPa$。

这个 c、φ 只代表该完整岩试件在所测围压范围内的线性等效值，不能直接赋给破碎岩体或天然结构面。若将拟合范围扩展到更高围压，等效 c、φ 也会改变，因此报告必须注明拟合区间。

数据质量与异常值处理

知识延伸

代表性通常比试验精度更重要：一台高精度试验机不能补救取样偏差。软弱破碎带最难取得标准试件，却往往控制边坡和巷道稳定。设计报告应同时披露“成功制样样本”和“无法制样区段”，必要时用原位试验、结构面模型或保守下界补充。

峰值强度不等于长期强度：已有剪切位移的断层、黏土充填面和长期受水软化的边坡，应关注残余强度。深部脆性岩则可能在峰值后突然掉落，弹性模型会低估剥落和岩爆风险。

强度与刚度承担不同角色：极限平衡主要使用抗剪强度；有限元/差分模型还需要 E、ν、胀角、抗拉强度、软化和渗流参数。只输入 c、φ 而随意取 E，会得到不可信的位移和塑性区范围。

常见误区

• 只试验完整的新鲜岩芯，弱夹层与破碎带没有参数
• 干燥岩芯代表长期饱和或雨季工况
• 用单次 UCS、均值或全球默认 k 代表整个地工域
• 用完整岩 c、φ 直接计算结构面滑动或岩体安全系数
• 三轴围压远高于浅部边坡应力，却将拟合结果外推至低围压
• 忽略层理方向、峰后软化、抗拉强度和变形参数
• 删除低值异常点但不调查其地质意义

关联章节

• 2.3.2 岩体分级 — UCS 输入 RMR
• 1.4.2 岩芯编录 — RQD 与取样
• 2.1.2 台阶与边坡参数 — 坡面角设计

自测要点

1. 完整岩、结构面和岩体强度分别服务哪些稳定性问题？
2. 为什么三轴试验的围压范围必须接近工程有效应力范围？
3. 若弱岩芯无法制成 UCS 试件，怎样避免数据库只保留“强者样本”？
4. 复算示例中的 c、φ，并解释为什么不能直接用于断层滑动分析。

2.3.2 岩体分级（RMR、Q、GSI）

学习目标

• 正确采集 RQD、节理组数、粗糙度、充填、地下水和应力信息
• 对同一岩体计算 RMR、Q、GSI，并解释三种体系为何可能给出不同判断
• 将 GSI、扰动系数 D、完整岩参数转换为 Hoek–Brown 岩体参数
• 明确经验分级、设计参数与支护决策之间的边界

核心概念

分级前先建立地工域

RMR、Q 和 GSI 都不应按整条钻孔或整面边坡平均。建议先建立地工域，再在域内统计分布与下界。常见分域维度包括：

• 岩性及蚀变：新鲜侵入岩、泥化带、矿岩接触带、构造角砾岩。
• 结构：块状、层状、片理化、断层破碎；节理组产状与工程方向的关系。
• 风化和水：新鲜、弱风化、强风化；干燥、渗水和承压水区。
• 应力与扰动：浅表松弛区、深部高应力区、爆破损伤带、采动影响区。
• 工程尺度：台阶、跨台阶边坡、巷道交叉口、采场顶板和矿柱可能由不同域控制。

参数与指标

RMR 分项（按所采用版本查表）

Q 体系参数

公式

• RQD：RQD = Σ(长度≥10 cm 的完整岩芯段) / 评价区间长度 × 100%
• RMR：$RMR = R_{1} + R_{2} + R_{3} + R_{4} + R_{5} \pm adjustment$
• Q 值：$Q = (RQD / Jn) \times (Jr / Ja) \times (Jw / SRF)$
• Q 体系等效尺寸：$D_{\mathrm{e}} = Span / ESR$（跨度或洞室高度除以开挖支护比 ESR；支护图需同时使用 Q 与 D_e）
• 经验相关式（仅供诊断）：$RMR \approx 9 \times \ln(Q) + 44$
• Hoek–Brown 岩体准则：$\sigma_{\mathrm{1}}' = \sigma_{\mathrm{3}}' + \sigma_{\mathrm{ci}} \times (m_{\mathrm{b}} \times \sigma_{\mathrm{3}}'/\sigma_{\mathrm{ci}} + s)^a$
• 岩体参数：$m_{\mathrm{b}} = m_{\mathrm{i}} \times \exp((GSI - 100)/(28 - 14D))$
• 岩体参数：$s = \exp((GSI - 100)/(9 - 3D))$
• 岩体参数：$a = 0.5 + (\exp(-GSI/15) - \exp(-20/3))/6$

$RMR \approx 9ln(Q)+44$ 是跨项目数据拟合，不是质量控制公式。RMR 与 Q 输入侧重不同，差异较大可能是真实反映地下水、应力或节理特征，也可能来自评分错误；应回到原始观测复核，而不是人为调整分数使两者一致。

工具与标准

• 常用软件/仪器：DIPS/同类赤平投影工具、RocData、三维地质建模软件、数字岩芯和结构面测绘系统
• 相关标准：GB/T 50218；ISRM 岩体表征建议

操作步骤

1. 定义评价尺度：明确是 1.5 m 岩芯回次、10 m 巷道段、单台阶还是跨台阶岩体。
2. 统一编录规则：对断裂、机械破碎、岩芯拼接、节理组、粗糙度、充填和水建立照片图谱。
3. 建立地工域：把地质、结构、强度、水和扰动整合为三维分区。
4. 独立计算体系：分别计算 RMR 与 Q，GSI 由结构和表面条件评价；保留原始分项。
5. 调查差异：用经验相关式和相邻区间作异常诊断，不为“对齐分数”改数据。
6. 导出强度参数：选 σ_ci、m_i、GSI、D，计算 m_b、s、a；记录每一项来源与范围。
7. 按应力范围等效：若软件需要 Mohr–Coulomb c、φ，应在目标围压区间拟合，不能全范围取一组常数。
8. 验证与更新：用揭露、变形、掉块、支护载荷和监测数据反演分域与参数。

岩体分级与参数导出流程

三大分级体系对照

可复算示例：GSI 与扰动对岩体参数的影响

设完整岩 $\sigma_{\mathrm{ci}} = 100 MPa$、$m_{\mathrm{i}} = 10$、$GSI = 55$。

两组参数来自同一个 GSI，但扰动使 m_b、s 显著降低。D 不是用来“调模型直到 FS 合格”的自由旋钮，应由爆破质量、开挖方法、损伤范围和现场观察支持，并通过敏感性分析展示不确定性。

从 Hoek–Brown 转换 c、φ 时的检查

Hoek–Brown 是非线性的，而 Mohr–Coulomb 是直线。软件输出的等效 c、φ 会随边坡高度、隧道埋深和所选围压上限改变。报告至少写明：

• 采用的 σ_ci、m_i、GSI、D 及其范围；
• 等效拟合的最小/最大有效围压；
• c、φ 是峰值、残余还是应变软化后的值；
• 参数用于岩体、断层、节理还是爆破损伤带；
• 地下水与有效应力如何处理，避免把孔压影响重复折减进强度。

知识延伸

分级是结构化描述，不是最终设计：RMR/Q 经验表建立于特定跨度、开挖和支护条件。大型采场、深部高应力巷道、永久高边坡和复杂断层必须用运动学、结构块体、极限平衡或数值方法专项验算。

RQD 有方向性和阈值效应：与钻孔轴线近平行的节理可能被漏掉；9.9 cm 与 10.1 cm 岩芯段地质意义接近，RQD 评分却不同。因此应同时报告裂隙频率、岩芯回收率和结构面产状，不能只保存一个 RQD 数字。

GSI 需要范围而非伪精确单值：例如 GSI 45–55 比 GSI=50 更诚实。把上下界传播到 m_b、s、a 和 FS/位移结果，才能看见参数不确定性。与 2.2.3 支护 和 2.1.2 边坡 形成“分级—参数—分析—监测”闭环。

常见误区

• 固定用 1.5 m 作所有 RQD 分母，忽略真实评价区间
• 钻探机械破碎也计为天然裂隙，或岩芯拼接后人为提高 RQD
• 用 RMR–Q 经验式强行修改评分，使两者“看起来一致”
• GSI 取单值但没有照片、结构描述和上下界
• D 被当作模型校准旋钮，而没有爆破/扰动证据
• 等效 c、φ 未注明拟合应力范围
• 直接用 RMR/Q 支护表签发大型采场或高边坡设计

关联章节

• 2.2.3 支护设计 — RMR 支护映射
• 2.3.3 边坡稳定性 — 岩体 c、φ
• 1.1.3 构造地质分析 — 结构面数据

自测要点

1. RMR、Q 和 GSI 的输入重点与主要用途分别是什么？
2. 为什么 RMR–Q 相关式不能作为评分验收标准？
3. RQD 应采用什么分母？为什么还必须同时记录节理方向和裂隙频率？
4. 复算 GSI 示例中的 m_b、s、a，并解释 D 从 0 增至 0.7 的工程含义。

2.3.3 边坡稳定性分析

学习目标

• 区分台阶、跨台阶和总体边坡的控制破坏模式与分析尺度
• 用赤平投影筛选平面、楔形和倾倒破坏，并理解运动学分析的边界
• 为圆弧、非圆弧、结构控制和复杂三维边坡选择合适的分析方法
• 建立地下水、地震、爆破损伤、采动顺序与长期软化工况
• 解释 FS、失效概率、位移和风险后果，形成可审计的设计决策

核心概念

三个设计尺度不能混用

台阶分析合格不能证明总体边坡合格；总体 FS 合格也不能防止局部楔体越过窄平台。设计必须同时满足局部可控和整体稳定。

分析输入与证据来源

设计接受准则

FS、P_f、允许位移和监测阈值应在分析前写入“设计接受准则”，并与边坡寿命、工况、后果和数据可信度匹配。永久高后果边坡通常比短期、低后果工作台阶要求更高裕度，但本章不提供可直接套用的统一数值。

接受准则至少回答：

• 分析对象是台阶、跨台阶还是总体边坡，服务到哪个开采阶段；
• 正常、雨季、极端地下水、地震和异常排水工况分别采用什么判据；
• FS 是均值参数、特征值还是下界参数下的结果；
• 是否还要求 P_f、位移、塑性区连通性、落石捕集率或 TARP 监测能力；
• 边际达标时需要哪些治理、观测法条件和复核频率。

公式

• 有效应力：$\sigma_{\mathrm{n}}' = \sigma_{\mathrm{n}} - u$
• 单一平面滑动简式：$FS = [cA + (Wcos\beta - U)\tan\phi ] / (Wsin\beta )$
• 锚索/外力的一般表达：在分子加入沿滑面抗滑分量，并在法向力中考虑其法向分量；由设计软件按几何分解
• 孔压比（筛查指标）：$r_{\mathrm{u}} = U / W$；正式分析优先使用实测孔压场或渗流模型
• 可靠度指标一阶近似：$\beta_{\mathrm{R}} = (\mu_{\mathrm{R}} - \mu_{\mathrm{S}}) / \sqrt{\sigma_{\mathrm{R}}^{2} + \sigma_{\mathrm{S}}^{2}}$
• 失效概率：$P_{\mathrm{f}} = P(g(X) \leq 0)$，其中极限状态函数可取 $g = FS - 1$
• 强度折减系数：$c_{\mathrm{reduced}} = c/SRF$，$\tan\phi_{\mathrm{reduced}} = \tan\phi /SRF$

平面滑动简式假设滑体刚性、单一平面、力学参数均一且未显式考虑张裂缝水压力、地震和三维约束，只适合手算复核。正式分析应按实际几何、孔压、外载和结构组合建模。

工具与标准

• 运动学/块体：DIPS、RocPlane、Swedge、UnWedge 或同类工具
• 极限平衡：Slide2/Slide3、SLOPE/W、同类 Bishop、Spencer、Morgenstern–Price 求解器
• 数值分析：RS2/RS3、FLAC2D/3D、3DEC/UDEC、PLAXIS；按连续体或不连续体机制选择
• 水文与监测：地下水有限元、测压管/孔压计、边坡雷达、棱镜、测斜与无人机点云
• 相关标准：项目所在地现行边坡、矿山安全、地震和监测规范，以及经批准的矿山边坡设计准则

第一步：运动学筛选

运动学分析只回答“几何上能否移动”，不计算岩桥、黏聚力、地下水和支护后的实际稳定裕度。

坡向容差、摩擦锥和横向限制应做敏感性。结构面产状分散时，用概率运动学统计“可动组合比例”，不要只画一条平均大圆。

第二步：选择力学模型

至少用一种独立方法复核高后果边坡。例如 LEM 给出 FS，数值模型检查塑性区是否沿预期滑面贯通；若两者破坏机制完全不同，应先回查模型，不应只比较数字是否相近。

第三步：建立分析工况矩阵

地下水模型应与水量平衡相容：边界补给、断层导水、排水孔流量和测压水位不能互相矛盾。仅在剖面上随手画一条“潜水线”通常是边坡分析最大的隐藏风险。

操作步骤

1. 明确设计问题：对象、寿命、后果、工况、接受准则和决策截止时间。
2. 审查地工模型：地形、地工域、断层、节理组、地下水和爆破损伤是否达到所需尺度。
3. 运动学筛选：按边坡分区和结构面分布识别平面、楔形与倾倒模式。
4. 选取代表剖面/三维域：覆盖最高、最陡、最湿、结构最不利及关键设施下方区域。
5. 建立基准模型：先复现当前几何、孔压和已观测行为，再预测未来阶段。
6. 运行工况矩阵：正常、雨季、排水失效、地震、长期软化和闭坑恢复水位。
7. 开展敏感性与概率分析：优先变化孔压、结构连续性、残余强度和地工域边界，而非只变化 c、φ。
8. 独立复核机制：比较 LEM 与数值结果、二维与三维、峰值与残余、排水前后破坏面。
9. 设计治理措施：优先消除暴露和控制水，再改几何、加固与设置监测/TARP。
10. 形成可追溯结论：记录控制滑面、关键假设、裕度、触发条件和需在施工中验证的预测。

边坡稳定性分析流程

破坏模式与分析方法

可复算示例：地下水对平面滑动的影响

假设刚性滑体重量 $W = 25 MN$，滑面倾角 $\beta = 38^{\circ}$，面积 $A = 800 m^{2}$，结构面 $c = 10 kPa = 0.01 MN/m^{2}$、$\phi = 32^{\circ}$。

干燥工况 U=0：

• 黏聚抗力：$cA = 0.01 \times 800 = 8.0 MN$
• 摩擦抗力：$Wcos\beta \tan\phi = 25 \times cos_{38}^{\circ} \times tan_{32}^{\circ} \approx 12.31 MN$
• 下滑力：$Wsin\beta = 25 \times sin_{38}^{\circ} \approx 15.39 MN$
• $FS_{\mathrm{dry}} = (8.0 + 12.31)/15.39 \approx 1.32$

滑面总水压力 U=4 MN：

• 有效法向力由约 19.70 MN 降至 15.70 MN
• $FS_{\mathrm{wet}} = [8.0 + (19.70 - 4.0)tan_{32}^{\circ}]/15.39 \approx 1.16$

水压力未改变几何与岩石重量，却使 FS 从约 1.32 降至 1.16。该例说明边坡分析中孔压往往比微调坡角或完整岩 UCS 更敏感。正式模型还需考虑张裂缝水、非均匀孔压、地震和三维约束。

敏感性与概率分析怎样做得有用

1. 先用 tornado 图或局部敏感性识别 FS/P_f 对哪些变量最敏感。
2. 对地质不确定性使用情景：断层连续/不连续、岩桥高/低、排水有效/失效。
3. 对可统计变量使用分布：φ、c、孔压、GSI、结构面产状；检查分布截断和相关性。
4. c 与 φ、GSI 与 D、降雨与孔压通常并非独立，不能默认全部独立随机抽样。
5. 报告 FS 分布、P_f、控制破坏面是否切换，以及尾部样本对应的物理情景。
6. P_f 低不代表后果可接受；应与人员暴露和设施后果结合为风险判断。

模型审查清单

稳定治理的优先顺序

1. 暴露控制：撤离人员、改道、设置落石和危险区隔离，先降低后果。
2. 水控制：截排地表水、水平排水孔、减压井、坡脚排水和维护能力验证。
3. 几何控制：卸载坡顶、压脚、降低跨台阶/总体角、加宽平台、调整坡道位置。
4. 采序控制：分期开挖、保留临时压脚、限制同时暴露长度、优化爆破方向和单段能量。
5. 加固与防护：锚索、抗滑桩、网喷、挡石设施；必须验证施工可达性、耐久性和失效模式。
6. 观测法：只有在传感器覆盖、数据可靠、TARP 权责和撤离时间均能满足时，监测才能作为设计控制的一部分。

知识延伸

FS 不是安全概率：FS=1.5 可能来自高置信度参数，也可能来自极少数据和乐观水位；两者风险完全不同。反之，FS 较低但有高可靠排水、实时雷达、无人暴露和快速撤离能力的短期工作台阶，其运营决策也不能只看一个数字。

二维模型的陷阱：二维剖面默认滑体沿垂直剖面无限延伸，可能忽略三维端部约束，也可能错过斜交断层形成的更不利滑动方向。沟谷、凸坡、有限长度断层、坡道切口和不规则坑壁应考虑三维验证。

爆破损伤应有空间范围：不应简单把整个边坡 c 降低固定百分比。更合理的做法是根据预裂半孔率、超挖、声波/成像、爆破记录和裂缝观察建立近坡面损伤带，并对厚度与强度做敏感性。

模型应预测可观测量：除 FS 外，模型还应预测位移方向、变形集中区、孔压响应和裂缝位置。监测若与预测不符，应更新模型假设，而不是只调整报警阈值。

常见误区

• 只分析总体边坡，忽略台阶楔体与落石；或反过来以台阶合格代替总体分析
• 赤平投影显示“可滑”便直接认定失稳，未检查迹长、岩桥和水
• 所有地工域使用同一 c、φ，断层和爆破损伤带没有单独参数
• 用一条主观潜水线代表复杂地下水，未用测压数据校准
• 只搜索圆弧滑面，漏掉弱层复合滑面和结构控制破坏
• 只报告最小 FS，不显示控制滑面、参数来源和工况
• 二维最安全剖面被选为代表剖面，未覆盖最弱三维组合
• 概率分析假设所有变量独立、正态且无物理上下限
• 数值模型未做网格、边界、本构和施工阶段敏感性
• 监测被当作补偿不合格设计的理由，但没有足够预警时间和撤离能力

关联章节

• 1.1.3 构造地质分析 — 结构面输入
• 2.1.2 台阶与边坡参数 — 几何迭代
• 2.3.4 地压监测 — 变形验证

自测要点

1. 台阶、跨台阶与总体边坡分别最容易漏掉哪些破坏模式？
2. 平面、楔形和倾倒破坏在赤平投影上的筛选逻辑是什么？
3. 为什么运动学可行不等于力学失稳？
4. Bishop、Spencer/Morgenstern–Price、FEM 强度折减和 DEM 各适合什么问题？
5. 复算地下水平面滑动示例，并说明 U 增加为何只降低摩擦抗力而不直接降低 cA。
6. 一个高后果边坡的设计工况矩阵至少应包括哪些水位、地震和长期情景？
7. 为什么 P_f 较低仍不能脱离失效后果直接判断风险可接受？

2.3.4 地压监测

学习目标

• 从预测破坏模式反推传感器类型、位置、量程、频率和冗余
• 区分位移、速度、加速度、孔压、支护载荷与微震趋势的工程意义
• 建立可执行的 TARP（触发—行动—响应计划），明确数据质量和授权责任
• 用多源监测反演边坡、巷道、采场和崩落区模型

核心概念

监测方法选择

没有一种仪器可以独立覆盖全部失效机制。高后果区域通常采用“面状表面监测 + 深部位移 + 孔压 + 人工巡检”的互补组合，并配置不同供电/通信路径。

公式

• 位移速度：$v_{\mathrm{i}} = (d_{\mathrm{i}} - d_{i-1})/(t_{\mathrm{i}} - t_{i-1})$
• 加速度：$a_{\mathrm{i}} = (v_{\mathrm{i}} - v_{i-1})/(t_{\mathrm{i}} - t_{i-1})$
• 收敛速率：$C_{\mathrm{rate}} = (C_{2} - C_{1}) / \Delta t$
• 反速度：1/v 对时间的趋势可用于特定渐进破坏的经验外推，但只在持续加速、机制稳定和数据质量可靠时使用
• 地震矩（概念）：$M_{0} = G \times A \times \Delta u$，G 为剪切模量、A 为破裂面积、Δu 为平均滑移
• 数据可用率：Availability = 有效数据时长 / 应采集时长 × 100%

短窗口会放大噪声，长窗口会延迟报警。速度与加速度应同时报告计算窗口、滤波方式、缺失值处理和单位，TARP 不能只写一个脱离算法的“mm/d”。

工具与标准

• 常用软件/仪器：边坡雷达、自动全站仪、GNSS、测斜/位移/孔压系统、微震阵列、LiDAR/无人机、InSAR 平台和集中告警系统
• 相关标准：经批准的矿山地压/边坡监测方案、仪器校准规程、数据治理标准、应急预案和项目 TARP

操作步骤

1. 写出失效假设：预计滑面、位移方向、先出现孔压还是位移、可用预警时间和受威胁对象。
2. 定义监测目标：设计验证、长期趋势、运营预警、支护验收或事故调查不可混为一套阈值。
3. 布置互补传感器：在驱动区、剪切带、稳定参考区和受体区布点；检查视线、量程、供电、通信和维护可达性。
4. 建立基线：覆盖足够的正常生产、爆破、降雨/融雪和温度变化；记录初始零点与安装状态。
5. 建立数据 QA：自动检查跳点、漂移、缺测、基准点移动、时间同步和传感器饱和。
6. 制定 TARP：每一级写清触发组合、数据确认、责任人、时限、停工/撤离范围和恢复条件。
7. 运行与交叉验证：将位移、孔压、降雨、爆破、开挖和支护施工画在统一时间轴。
8. 反演模型：比较预测与实测位移方向、深度、孔压和时序，更新参数、结构或地下水边界。
9. 事件后复盘：确认阈值是否过早/过晚、仪器是否覆盖控制区、响应动作是否按时完成。

监测工作流

TARP 模板：阈值必须现场标定

实际 TARP 可以包含数值阈值，但应来自基线统计、模型预测、设备噪声、历史事件和所需撤离提前量。不能把其他矿山的 1、2、5 mm/d 直接复制为本站阈值。

数据异常与真实变形的区分

可复算示例：窗口选择如何改变报警

某棱镜位移为第 0 天 10 mm、第 1 天 11 mm、第 2 天 12 mm、第 3 天 16 mm、第 4 天 22 mm。

• 日速度依次为 1、1、4、6 mm/d，显示最近两天明显加速。
• 若只计算 4 天平均速度：$(22-10)/4 = 3 mm/d$，会掩盖当前 6 mm/d 的增长。
• 若只看最后两个点，单次测量噪声又可能被放大。

因此 TARP 常同时使用短窗口速度、较长窗口趋势、加速度和多源确认。窗口长度应与破坏发展时间和测量噪声共同确定。

知识延伸

监测系统也是安全关键设备：需对覆盖率、延迟、误报、漏报、备用电源、通信冗余和人员值守进行完整性管理。传感器离线不是“没有位移”，而是失去一层风险控制。

领先指标与滞后指标：孔压上升、微震迁移和支护载荷变化可能早于可见位移；大裂缝和掉块通常已是较晚信号。监测设计应争取足够行动时间，而不是只记录破坏发生。

反分析不等于任意调参：只有当多个观测量能被同一参数集解释时，反演才有意义。若只把 E 调低以匹配一个棱镜位移，却无法解释滑面深度和孔压响应，模型仍未校准。

常见误区

• 先购买仪器再思考破坏模式，结果关键滑面和盲区没有覆盖
• 把其他矿山的位移速度阈值直接复制到本站 TARP
• 只设位移触发，不监测孔压、加速度、支护载荷和数据可用率
• 仅坡顶棱镜无深部测斜，无法判断深部滑面
• 植被、陡坡或快速变形区只依赖 InSAR，无地面验证
• 单个传感器报警即下结论，或反过来等待所有传感器都报警才行动
• 传感器离线时默认风险未变化，继续高暴露作业
• 治理、改坡或仪器更换后沿用旧基线和旧阈值
• TARP 写了颜色但没有责任人、时间限制、撤离范围和复工条件

关联章节

• 2.3.3 边坡稳定性分析 — 模型校准
• 2.1.2 台阶与边坡参数 — 治理改坡
• 7.2 安全管理 — TARP 制度

自测要点

1. 位移、速度、加速度和孔压在失稳演化中分别提供什么信息？
2. 如何用雷达、棱镜、测斜和孔压组成互补的露天边坡监测系统？
3. 为什么传感器离线本身应成为 TARP 的触发条件？
4. 复算窗口示例，并比较 1 天和 4 天速度对运营决策的影响。
5. TARP 的数值阈值应由哪些现场证据共同确定？

2.3.5 地下开挖、采场与矿柱稳定性

学习目标

• 识别地下开挖中的结构块体、脆性剥落、挤压、膨胀和岩爆机制
• 从原位应力、岩体分级、跨度和采矿顺序选择分析方法
• 用矿柱承载、应力路径、强度折减和收敛约束校核稳定性
• 将分析结果转化为开挖顺序、支护时机和监测要求

核心概念

地下稳定性问题与方法

参数与公式

• 自重垂直应力初估：$\sigma_{\mathrm{v}} \approx ∫\rho g dz \approx \rho gH$；构造应力显著时不能只用覆岩自重
• 水平/垂直应力比：$k = \sigma_{\mathrm{h}}/\sigma_{\mathrm{v}}$，应由水压致裂、套孔解除或其他原位测试与反演约束
• 矿柱安全系数：$FS_{\mathrm{pillar}} = Pillar strength / Pillar demand$
• 面积分担法柱载荷（规则采场初估）：$\sigma_{\mathrm{p}} \approx \sigma_{\mathrm{v}} \times A_{\mathrm{tributary}}/A_{\mathrm{pillar}}$
• 水力半径：HR = 暴露面积/暴露周长，用于采场稳定图的几何尺度
• 收敛—约束平衡：围岩反应曲线与支护特性曲线交点给出平衡压力和位移
• 强度折减：用于连续体模型估算整体失稳裕度，但不替代显式楔体和动力分析

面积分担法假设规则几何、均匀载荷和柱间相互作用简化，仅适合初筛。孤立柱、非规则采空区、倾斜矿体、充填体和回采顺序必须用三维模型与现场应力/变形验证。

工具与标准

• 结构与采场：DIPS、UnWedge、Mathews/Potvin 稳定图工具
• 应力与连续体：Map3D、RS2/RS3、FLAC2D/3D、Abaqus 等
• 显式结构/崩落：UDEC/3DEC、离散元或混合连续—不连续方法
• 监测：收敛、MPBX、多点应力计、锚杆/锚索测力、微震、扫描和损伤测绘
• 设计依据：项目所在地地下矿山、支护、动力灾害和职业安全规范，以及矿山批准的地压管理计划

操作步骤

1. 定义稳定对象：巷道、交叉口、硐室、采场、矿柱、充填体或区域采动系统。
2. 建立三维地工模型：岩性、断层、节理、地工域、采空区、充填和相邻工程。
3. 确定原位应力：结合测量、区域构造和开挖损伤反演，给出方向与不确定范围。
4. 先判机制：结构控制、应力控制、重力控制、挤压或动力破坏决定方法选择。
5. 做经验初筛：RMR/Q、稳定图、面积分担法或弹性应力模型用于快速识别风险区。
6. 建立阶段模型：按真实开挖、充填和支护安装顺序模拟，不能只计算最终空区。
7. 校核局部与整体：楔体、支护构件、矿柱、采场帮和区域应力转移分别检查。
8. 敏感性与失效路径：变化原位应力、弱层、峰后软化、支护时机和充填刚度。
9. 转化为施工控制：最大无支护暴露、支护类型与时机、采序、禁入区和 TARP。
10. 用监测反演：收敛、损伤深度、支护载荷和微震位置与模型预测比较并更新。

地下稳定性分析流程

原位应力与模型边界检查

采场与矿柱稳定的系统思维

矿柱“平均 FS 合格”并不表示采区稳定。需要同时检查：

• 柱边峰值应力与剥落后有效柱宽是否持续缩小；
• 相邻采场回采是否把载荷转移到剩余矿柱和主运输巷道；
• 充填体接顶、刚度和龄期是否能按计划分担载荷；
• 局部断层是否切断矿柱或形成贯通滑移面；
• 回采顺序是否造成高能量释放和微震簇迁移；
• 永久矿柱的长期蠕变、爆破损伤、浸水和闭坑水位恢复。

可复算示例：面积分担法初筛

规则房柱布置中，柱宽 6 m、房宽 9 m，按方形单元估算：

• 单元面积：$A_{\mathrm{tributary}} = (6 + 9)^{2} = 225 m^{2}$
• 柱面积：$A_{\mathrm{pillar}} = 6^{2} = 36 m^{2}$
• 若原位垂直应力 $\sigma_{\mathrm{v}} = 12 MPa$，则平均柱应力初估：$\sigma_{\mathrm{p}} = 12 \times 225/36 = 75 MPa$
• 若采用经尺寸效应修正的柱强度为 105 MPa，则 $FS_{\mathrm{pillar}} \approx 105/75 = 1.40$

这个结果只能作为筛查：它没有考虑柱边剥落、非均匀载荷、断层切割、邻区回采、水平应力和充填承载。若后果高或几何不规则，应进入三维阶段模型。

支护设计如何与稳定分析衔接

支护不是事后附加荷载。模型应体现：

• 锚杆/锚索安装时围岩已经发生多少位移；
• 锚固段是否穿过潜在楔体或塑性区进入稳定岩体；
• 喷层、网和钢架提供面约束还是承担大变形；
• 构件是否具有足够延性，能否适应峰后位移而不脆断；
• 交叉口、断层、采场口和爆破损伤区是否单独加强；
• 拉拔试验、喷层厚度、预应力和竣工位置是否进入 QA 数据库。

详细支护选型见 2.2.3 支护设计。

常见误区

• 只用 RMR/Q 支护表，不分析楔体、应力集中和采矿顺序
• 原位应力完全采用 $\rho gH$，忽略构造应力与方向
• 二维巷道剖面代表三维交叉口、采场端部和矿柱系统
• 模型一次挖出最终采空区，未模拟支护、充填与阶段开采
• 只看弹性最大应力，不检查峰后软化、损伤深度和能量释放
• 矿柱强度未做宽高比与尺度修正，柱载荷只用平均值
• 数值位移与现场坐标、基准和施工时序无法对应
• 微震事件率上升被简单等同于即将岩爆，未检查目录完整性和空间迁移

关联章节

• 2.2.2 采矿方法选择 — 采场几何、充填与回采顺序
• 2.2.3 支护设计 — 锚杆、喷层与钢架
• 2.3.1 岩石强度试验 — 强度、刚度与峰后参数
• 2.3.2 岩体分级 — RMR、Q、GSI 参数入口
• 2.3.4 地压监测 — 收敛、支护载荷与微震反演
• 7.2 安全管理 — 地压风险与作业控制

自测要点

1. 结构控制楔体、脆性剥落和挤压大变形的证据分别是什么？
2. 为什么地下稳定模型必须按开挖、支护和充填顺序分阶段计算？
3. 面积分担法适用于哪些条件，又忽略了哪些关键三维效应？
4. 复算矿柱示例；若房宽增至 12 m，平均柱应力和 FS 如何变化？
5. 收敛—约束法中，支护安装时机为什么会改变最终平衡点？