综贯经验-分析-数值模拟的 Flysch 地层隧道稳定性评估：以泽尼察隧道为例

《Geotechnics》：Integrated Empirical–Analytical–Numerical Assessment of Tunnel Stability in Flysch: A Case Study of the Zenica Tunnel Ekrem Bekta?evi?, Luka Crnogorac, Kemal Guti?, Vancho Adjiski, Rade Tokali? and Ahmed Mu?ija

为解决复杂地质条件下隧道稳定性评估难题，研究人员开展了一项针对 Flysch 地层中公路隧道稳定性的综合研究。该研究以泽尼察隧道为案例，整合了岩体分级、收敛-约束法（CCM）和非线性有限元（FEM）分析，建立了稳定性系数 Ns与隧道收敛之间的强指数关系，并识别出 RMR < 25 和 Ns> 5 为显著挤压行为的临界阈值，为类似地质条件下的隧道稳定性评估与支护设计提供了更可靠的方法论。

　　

在群山峻岭中开凿隧道，是现代基建征服天堑的壮举，却也伴随着巨大的技术挑战。尤其当隧道穿越复杂多变的地质体时，工程师们仿佛在与“地球的脾气”进行一场无声的博弈。其中，Flysch地层（一种由砂岩、泥灰岩、粘土岩等频繁互层构成的复杂沉积序列）因其显著的异质性和各向异性，常常成为隧道工程的“梦魇”。在这样的地层中，岩体质量可能在短距离内剧烈变化，传统的单一评估方法往往顾此失彼，难以精确预测隧道开挖后的变形行为和支持需求。那么，如何才能更可靠地评估这类“难缠”地层中隧道的稳定性，从而指导更安全、更经济的支护设计呢？近期发表在《Geotechnics》上的一项研究，为我们提供了一套集成化的解决方案。

为了回答上述问题，以Ekrem Bekta?evi?、Luka Crnogorac、Kemal Guti?、Vanjo Adjiski、Rade Tokali?和Ahmed Mu?ija组成的研究团队，选择波黑Corridor Vc走廊上的泽尼察公路隧道（Zenica road Tunnel）作为典型案例，开展了一项综合性研究。他们创造性地构建了一个融合了经验分级、分析模型和数值模拟的“三维一体”研究框架。具体而言，该研究将国际上广泛应用的岩体评分（RMR）系统、用于评估隧道与围岩相互作用的收敛-约束法（Convergence–Confinement Method, CCM），以及基于Hoek-Brown破坏准则的非线性二维有限元（FEM）分析结合在一起。通过这套组合方法，研究人员不仅量化了不同地质条件下的隧道稳定性，更重要的是，揭示了控制隧道变形的关键阈值，为工程设计决策提供了清晰的科学依据。

主要技术方法：本研究采用了一种集成的分析框架。首先，基于现场地质勘察和实验室测试数据，应用岩体评分（RMR）系统对隧道沿线的岩体质量进行分类。其次，运用收敛-约束法（CCM）进行解析计算，核心参数包括稳定性系数（N_s= σ₀/σ_cm，即初始应力与岩体强度之比）和临界应力释放度（λ_cr）。最后，使用PLAXIS 2D软件进行非线性有限元（FEM）数值模拟，以验证解析结果并更精细地分析应力-应变场。研究对象为泽尼察隧道右线隧道（RTT）超过三分之二长度的代表性断面，覆盖了从Ⅲ类到Ⅴ类的不同岩体质量。

研究结果

3.1. 稳定性系数N_s的解析评估

通过计算隧道沿线各里程的稳定性系数N_s，研究发现其值在0.13到7.66之间大幅波动，清晰地反映了飞复理石地层的极度不均一性。当N_s< 1时，岩体主要表现为弹性行为；N_s> 1则进入弹塑性行为区，开始形成塑性区。最为关键的发现是，在岩体质量最差的V类岩体（RMR ≤ 20）区段，特别是里程0+689.21处，N_s值高达7.66，标志着岩体已进入高度不稳定状态。这表明N_s系数能够可靠地综合反映上覆岩层厚度和岩体质量的双重影响。

3.2. 相对变形分析与稳定性阈值

基于N_s和σ_cm/σ₀比值，研究预测了隧道轮廓的相对径向变形ε(%)。结果表明，变形范围从最小的0.04%到最大的11.73%。根据既定的稳定性图表标准（变形小于1%为“微小支护问题”，1-2.5%为“轻微挤压问题”，2.5-5%为“严重问题”，5-10%为“非常严重的挤压问题”，大于10%为“极端挤压问题”），研究成功地将隧道各段按风险等级进行了分类。例如，在浅埋段（H < 100米），变形保持在1%以下，处于弹性稳定状态；而在前述N_s=7.66的极端区段，预测变形高达11.73%，属于“极端挤压问题”，支护系统有失效风险。这一分类为差异化支护设计提供了直接依据。

3.3. 临界应力释放度与隧道稳定性的关系

在收敛-约束法框架内，临界应力释放度λ_cr是岩体行为从稳定向不稳定转变的阈值。研究发现，在不稳定区段（N_s> 5），λ_cr值很高（如0.87），这意味着绝大部分潜在变形在开挖面附近就迅速释放，要求支护系统必须快速安装并生效。相反，在稳定的Ⅲ类岩体区段，λ_cr值较低（如0.26），应力释放过程平缓，对支护的即时性要求较低。这一定量关系为确定最佳支护时机提供了关键参数。

3.4. 数值验证与隧道稳定性的空间分析

研究通过PLAXIS 2D软件对解析模型进行了数值验证。有限元模拟结果显示，最大位移（29.0至35.2毫米）恰好出现在RMR ≤ 20的V类岩体区段，与解析模型识别的高N_s值关键区完全吻合。统计验证显示，解析预测与数值模拟结果之间具有极强的相关性，决定系数R²介于0.89至0.96之间，证明了所建立模型的内部一致性和高可靠性。比较发现，收敛-约束法（CCM）的估算相对保守，预测的变形值通常比有限元（FEM）结果更高，但两者在识别关键不稳定区域上结论完全一致。这表明CCM适用于全局稳定性评估和设计基准的建立，而FEM则能更真实地模拟局部位移和塑性区分布，用于优化支护细节。

3.5. 研究的局限性

研究也坦诚了其局限性。首先，采用的二维有限元模型和β-法（应力释放法）是对三维开挖面效应的近似，在高度各向异性的Flysch地层中，可能对塑性区分布的模拟存在细微偏差。其次，分析主要关注岩体的初始弹塑性响应，未考虑软弱岩体在高地应力下可能存在的长期蠕变和强度衰减的时间效应。此外，地下水的影响仅通过静力降低有效应力的方式处理。未来的研究需要引入完全耦合的水-力模型和三维模拟，以更精确地定义Flysch层理与非对称轮廓位移之间的相互作用。

研究结论与意义

本项关于泽尼察隧道的研究，雄辩地证明了在异质Flysch地层中进行地下结构稳定性预测时，采用集成化方法具有至关重要的意义。研究摒弃了依赖单一方法的传统思路，通过将经验分级、解析计算和数值验证三者有机结合，全面揭示了岩体力学行为的内在规律。研究表明，设计方案的可靠性恰恰源于这些方法之间的相互印证与补充，从而显著降低了在Flysch这类复杂地质环境中进行隧道工程建设的固有风险。

本研究的核心科学贡献在于量化了一个明确的行为边界：当岩体评分（RMR）指数低于25，且稳定性系数（N_s）大于5时，标志着隧道围岩行为从可控的、近似线性的响应，向非线性的、显著的挤压行为发生临界转变。尽管传统的收敛-约束法（CCM）倾向于给出保守的、偏大的理论位移预测，但在工程实践中，这恰恰构成了一个重要的安全缓冲。另一方面，研究所采用的数值模型证实，支护系统的适时激活，能够有效地将一个理论上的不稳定岩体，转变为一个受控的岩土工程体系，并将安全系数稳定在1.38至1.45的可靠范围内。

研究揭示的高统计相关性（R²> 0.89）表明，所建立的模型不仅适用于泽尼察隧道这一具体案例，更为迪纳拉山脉（Dinarides）乃至全球类似地质环境下的隧道设计提供了一个普适性更强的分析框架。这项研究倡导的是一种基于早期预测参数的主动性设计方法，而非仅仅依赖施工过程中的被动监测和调整。这种前瞻性的方法，不仅确保了在最挑战性的岩体类别中开挖隧道的稳定性，也为在Flysch地层中进行更可预测、更经济，且至关重要的——更安全的隧道工程建设，奠定了坚实的科学基础。