不排水黏土中隧道开挖引发地层移动的概率分析：经验方法再探讨

《Soils and Foundations》：Probabilistic analysis of tunneling-induced ground movements in undrained clay: empirical methods revisited

　　

隧道开挖引起的地层移动可能对城市既有基础设施造成损坏。地层对隧道开挖的响应已有广泛研究，但土体刚度空间变异性对地层移动的影响尚不清楚。本文针对不排水黏土（undrained clay）中隧道上方绿地（greenfield）地层移动进行概率分析。将有限元法（FEM）与随机场理论（random field theory）相结合，在蒙特卡洛模拟（Monte?Carlo simulation, MCS）框架下，获取不同埋深隧道在黏土不排水刚度（undrained stiffness, Eu）具空间变异性时的随机地层移动剖面。结果表明，黏土刚度的空间变异会增大地层移动，尤其对浅埋隧道影响显著。利用计算结果及汇编的案例历史数据，结合现有经验预测方法重新评估了与地层移动分布和量值相关的参数。给出了系统失效概率（system failure probability），即最大地层移动超过允许值的概率。

《不排水黏土中隧道开挖引发地层移动的概率分析：经验方法再探讨》论文解读

本文发表于《Soils and Foundations》。既有隧道开挖引发的地层移动（ground movement）研究中，经验方法（如Peck公式、O'Reilly & New公式等）和解析解多基于确定性假设，即将黏土的不排水弹性模量E_u取为均一常数，忽略了天然土体固有的空间变异性（spatial variability）。而土体刚度变异会直接影响隧道周围变形场的发育，传统确定性方法可能低估地表最大沉降与水平位移。此外，已有考虑土体随机性的研究多集中于砂土排水条件，针对不排水（constant?volume）黏土中隧道诱发变形的概率特性研究较为缺乏。因此，研究人员开展了将有限元分析（finite element analysis, FEA）与随机场理论（random field theory）耦合的蒙特卡洛模拟（Monte?Carlo simulation, MCS）研究，定量评价不排水黏土E_u的空间相关性对地表竖向及水平向绿地（greenfield）地层移动统计特征的影响，并基于数值结果与30组现场案例数据库重新校核经典经验方法中的沟槽宽度参数K（trough width parameter）及体积损失V_L（volume loss），给出以最大地表沉降或水平位移超越允许值为判据的系统失效概率（system failure probability, p_f），为城市浅埋盾构隧道初步设计提供概率依据。

主要关键技术方法：

研究人员建立二维平面应变有限元模型（OPTUM软件），隧道直径D＝5 m，埋深比H/D＝1、2、3、5，土域尺寸12D×8D，侧边界限竖向自由、底边界固定。不排水泊松比ν_u＝0.5，E_u服从对数正态分布（lognormal distribution），均值μ_Eu＝20 MPa，变异系数COV_Eu＝30％，垂直与水平相关距离θ_v＝4 m、θ_h＝50 m（指数型相关函数），采用Karhunen?Loève展开生成随机场并映射至自适应加密网格节点。开挖模拟采用收敛—约束法（convergence–confinement method），应力释放因子λ＝0.1（对应体积损失约1％）。对每个H/D进行400次MCS以获得地表沉降与水平位移的统计特征，并以Kolmogorov?Smirnov检验验证最大值服从正态分布；同时汇编全球30例黏土中隧道地表沉降案例历史数据库用于经验参数校核与对比。

研究结果：

2. Empirical estimation of ground movements——竖向与水平向经验方法回顾

综述了Peck(1969)、O'Reilly & New(1982)、Mair(1979)、Loganathan & Poulos(1998)等给出的拐点距x_i与覆土深度H的关系式（x_i＝KH或幂律形式2x_i/D＝a(H/D)^b），以及Schmidt(1987)建议的水平最大位移与竖向最大位移比h_max/v_max＝0.36(H/D)^0.2，明确本研究待校核的经验参数。

4. Variability of ground movements——地层移动的变异性分析

MCS得到的地表沉降槽仍近似符合高斯分布（Gaussian distribution），但最大沉降不一定严格位于隧道轴线且剖面呈非对称，源于随机场中较软区域产生更大变形。随H/D增大（深埋），地表沉降变异性减小——深层隧道因更大体积的空间平均效应削弱E_u局部波动影响。v_max与h_max的均值μ_vmax、μ_hmax均大于均匀土体的确定性解v_max,det、h_max,det，μ_vmax/v_max,det增幅1.49％～13.82％，表明忽略空间变异性会低估地层移动。v_max与h_max的COV（13.48％～16.01％）小于E_u的COV（30％），符合局部变形区空间平均效应。KS检验确认v_max、h_max近似服从正态分布。

5. Evaluation of empirical method——经验方法参数再评估

基于数值结果与30例现场数据库拟合得：x_i＝0.51H（线性），与Mair(1979)及O'Reilly & New(1982)吻合；归一化关系2x_i/D＝1.68(H/D)^0.62，在H/D＞3时接近Peck(1969)与Loganathan & Poulos(1998)建议值。由v_max–x_i关系反算体积损失V_L≈1.02％，与FE模型中λ＝0.1所产生体积损失（0.8％～1.6％）一致。h_max/v_max在H/D＝2附近出现峰值后略降，变化范围0.19～0.29，整体低于Cording & Hansmire(1975)建议上限0.4。

6. Probability of failure——失效概率分析

假定v_max、h_max服从正态分布N(μ,σ2)，定义失效为隧道诱发最大移动超过允许值α·v_max,det（或α·h_max,det），α为调整系数。给出不同H/D及α下p_f值表：当α≥1.6时p_f＜1％。以FHWA(2009)建议的允许沉降v_cri＝10 mm为限值，H/D≥3.57时p_f≤1％，浅埋（H/D＝1）时p_f显著升高。表明考虑土体空间变异性后，浅埋隧道需更保守的允许位移控制或更大的安全余量。

讨论与结论（翻译浓缩）：

研究表明，黏土不排水刚度E_u的空间变异性使地表最大沉降与水平位移均值较均匀土体增大1.5％～13.8％，且地表移动剖面呈非对称性；空间变异性对浅埋隧道（小H/D）影响更为显著。基于随机有限元结果与30组现场案例重新标定得拐点距经验式x_i＝0.51H（K≈0.51），体积损失V_L≈1.0％～1.02％与模型输入一致，h_max/v_max典型值0.19～0.29。考虑刚度空间变异给出以允许位移为阈值的系统失效概率，当调整系数α＞1.6可使p_f＜1％。研究成果可为黏土中城市浅埋盾构隧道初步设计阶段的概率化地层移动预测与风险评估提供参考；研究限于线弹性不排水黏土模型，后续可拓展至更复杂本构关系。