谈高水压山岭隧道衬砌水压力的计算

谈高水压山岭隧道衬砌水压力的计算仇文革高新强（西南交通大学土木工程学院成都 610031）[摘要] 在高水压岩溶区修建隧道在国内外都是个技术难题，采用理论分析方法和等效连续介质模型用数值方法分析了隧道渗流场的分布，对地下水在围岩、注浆圈、衬砌上水压力分布规律及其作用系数进行了分析。

研究结果表明：衬砌背后水压力的大小与围岩、注浆圈、衬砌的厚度、渗透系数和隧道的控制排水量关系密切。

[关键词] 山岭隧道；衬砌水压力；水压力作用系数；水压力计算1 引言山岭隧道地下水流是处于地下水流系统之中，隧道所处的水力单元不同隧道位置的渗流场也有所不同。

在山体中开挖隧道后地下水的所有补给、排泄状态见图1。

尽管隧道围岩裂隙发育，注浆堵水后仍可将围岩、注浆圈、衬砌看成是符合达西定律的串联水流网络；若围岩中存在大的岩溶管道，可以视为并联水流网络；若两种都存在可以视为混合水流网络。

图1 隧道区地下水补给与排泄的示意图隧道开挖后根据地下水补给量（Q g）和隧道最大排水量（Q p）的大小关系可以产生3种演变模式（见图2）。

(b) Qg <<Qp (b) Qg <Qp(c) Qg >>Qp图2 隧道衬砌水压力3种模式（a ）种模式，地下水补给量（Q g ）远远小于隧道最大排水量（Q p ），地下水水位降低到隧道以下，衬砌上的水压力等于0。

（b ）种模式，地下水补给量（Q g ）小于隧道最大排水量（Q p ），地下水水位降低到一定程度后处于稳定状态，隧道处的水头小于原始水头0H γ，衬砌上有可能存在一定的水压力。

（c ）种模式，地下水补给量（Q g ）远远大于隧道最大排水量（Q p ），地下水水位处于稳定状态，若全封堵，隧道处的水头基本等于原始水头0H γ，在排水状态下，衬砌上也有可能存在一定的水压力。

本文主要讨论这种条件下衬砌上的水压力。

2 均质围岩中水压力的计算 2.1 浅埋隧道隧道处于浅埋时，计算衬砌水压力和最大涌水量可以分别采用（1）式和（2）式计算[3]（二维非轴对称流）。

）（w r HKQ H P 000ln 2πγ−= （1）wr KH Q 00max ln 2π=（2） 式中，γ为水的容重，N/m 3；P 为作用于衬砌外表面的水压力，Pa ； K 为岩层等效渗透系数，m/s ； H 0为隧道原始静水头，m ； γw 为开挖隧道断面等价圆半径，m ；h w 为开挖隧道断面上的水头，m ； Q max 为隧道最大涌水量，m 3/s 。

图3 围岩、注浆圈和衬砌各层渗流模型（c ）衬砌（b ）注浆圈2.2 深埋隧道在一些山岭隧道的修建中，隧道的埋置深度大，隧道穿越富水区，水头高（如圆梁山隧道和锦屏交通隧道），如何在不影响生态环境的条件下，安全、可靠地进行隧道设计、施工和运营一个关键的问题就是如何确定衬砌上水压力大小[1~2]。

目前，在隧道衬砌水荷载的计算中，铁路、交通部门还没有制定统一的规范，大多还是参照水工隧洞设计规范和经验方法，有关水荷载的论述也散见于各部门和学科的专著及专业杂志上。

水工隧洞中衬砌水荷载一般包括内水压力和外水压力两部分[4]，二者作用对象均为衬砌（对于围岩可直接称为水压力）。

外水压力是和有压隧洞中内水压力相对而言的，而铁路、公路隧道一般不存在内水压力，通常简称为“水压力”[5]，故其衬砌水荷载与外水压力指的是同一概念。

若深埋隧道围岩完整，裂隙不发育时，可以假定围岩为均质各向同性，根据地下水动力学理论，以无限体含水层中的竖井进行理论分析（可视为轴对称流）。

根据图3的隧道渗流模型推导得到计算衬砌背后水压力的计算公式⎟⎟⎟⎟⎟⎠⎞⎜⎜⎜⎜⎜⎝⎛++−+==g r r l l g g l l lr w l w l r r K K r r K K r r r r H H H H P ln ln ln ln )(1111γγ （3） 式中： ——每延米隧道排水水量； 0Q ——衬砌内表面水头； 1H ——衬砌背后水头； l H ——注浆圈圈外表面水头； g H ——围岩表面水头； r H ——围岩渗透系数； r K ——注浆圈渗透系数； g K ——衬砌渗透系数； l K ——衬砌内半径； 1r ——衬砌外半径； l r ——注浆圈半径； g r ——围岩半径；r rw γ——水的容重；——衬砌背后水压力； l P ——注浆圈外表面水压力。

g P 铁路隧道衬砌内表面的水压力认为等于0，围岩的半径等于地下水的静水头，则（3）式可以表示为：1H r r 0H gr l l g g l l lw l w l r H K K r r K K r r r r H H P 0110ln ln ln ln++==γγ （4）若衬砌为复合式衬砌，衬砌背后设防水布、防水布后设置无纺布渗透层和盲管排水系统，此时设衬砌不透水，从注浆圈渗出的地下水只从排水系统排出。

⎟⎟⎠⎞⎜⎜⎝⎛+−=g l g g r rw w l r H r r K K K QH P 000ln ln 2πγγ （5） （4）、（5）式中符号的含义同（3）式。

3 裂隙围岩中水压力的计算 3.1 计算模型裂隙围岩中水压力的计算以圆梁山隧道毛坝向斜区[6]为例进行阐述，该岩溶管道在接近地表的浅部极其发育。

岩溶发育高程的下限在海拔900左右，900米以下岩溶现象不发育，而隧道的高程在550米，深孔钻探稳定水位为999.69米，说明在高程900米到999.69米以上地下水运动以岩溶管道和溶隙流为主，具有非线性特性，而在550米高程到900米高程，溶蚀现象不明显，和隧道发生直接水力联系的地下水系统应该是裂隙系统，900米高程以上岩溶管道系统对裂隙系统起补给源的作用，和隧道不直接发生水力联系。

针对研究的区域可以近似忽略900米高程以上地下水运动的非线性特性并将其处理为线性流，通过使用等效岩体水力学参数将其概化为等效连续介质模型（式6），对隧道周围渗流场进行数值模拟。

⎪⎪⎪⎩⎪⎪⎪⎨⎧≥Γ∈′=∂∂+∂∂≥Γ∈=Ω′∈′=Ω′∈≥∂∂=−−+⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛∂∂∂∂+⎟⎠⎞⎜⎝⎛∂∂∂∂′02011100***,),(),,,(),cos(),cos(,),(),,,(),,(),(),,(),,(),(,,),(t t z x t z x q z n z H T x n xH T t t z x t z x H t z x H z x z x H t z x H z x t t t H S z z x x Q z H K z x H K x M zz xx b ii i zz xx δ （6）式中，H ——岩溶地下水水头； **,yy xx K K 和——分别为等效裂隙介质的等效渗透系数和等效储水系数；*S M ′——岩溶介质含水层厚度（对于非承压水，使用自由水头）；i Q ——第i 口开采井抽水量；δ——狄拉克函数；*0H ——初始水头； *1H ——已知定水头；1Γ、——第一、二类边界；2Γ*1*1,yy yy xx xx K M T K M T ⋅=⋅=——等效导水系数；Ω′——地下水渗流区域。

3.2 边界条件分析隧道毛坝向斜段的地层依次为S 、D 的泥岩和砂质泥岩，与的石灰岩，和的石灰岩、的薄层泥岩、石灰岩等。

根据深孔抽水试验，与、的导水性较好，试验中，水位恢复快，为该区的主要富水岩层；和则富水、导水性较差，试验抽水后，水位不能恢复；q l P +1m P 1c P 2d T 1w P 2q l P +1m P 1d T 1c P 2w P 2S 、D 的地层基本可以看作隔水层。

隧道穿越地层主要是与的石灰岩，和的薄层泥岩、石灰岩等。

由于地层的封隔，除断层导入外，和与没有产生直接的水力联系。

隧道附近有北西向平移断层经过，其中两条和隧道走向近于平行的平移断层和规模较大，从1：10000工程地质图上分析，基本切穿了以上地层，考虑到上层的富水岩层由于断层的导入作用，在分析隧道穿越的与地层渗流场时，将这两条断层选作定水头边界。

当与隧道围岩向隧道排水时，假定通过断层补给这两个岩层，并保持自己的水头不变。

据隧道地勘报告报告，900m 标高以上岩体岩溶比较发育，岩溶水的快速径流与补给是可能形成定水头边界的。

q l P +1m P 1c P 2w P 2w P 2d T 1q l P +1m P 16F 10F d T 1q l P +1m P 1q l P +1m P 1d T 1S 、D 的地层则处理为与隧道围岩渗流区东西和底板的隔水边界。

上边界取为水头边界条件，相对隧道中心位置水头为490m ，下边界、左右边界取为不透水边界；隧道衬砌内边界水压力为0。

q l P +1m P 13.3 地下水渗流模型岩体水力学参数分析主要根据深孔压水试验和抽水试验成果与在野外根据不同标高节理裂隙统计计算的渗透张量给出。

深孔压水试验和抽水试验成果各层渗透系数见表1和2，在野外根据不同标高节理裂隙统计计算的渗透张量见表3。

计算分析采用表4中渗透系数。

1_1−Y Z 1_1−Y Z深孔压水试验成果表 表11_1−Y Z 深孔抽水试验成果表 表2试验段层试验段厚度M（m ）静止水位（m ）降深S （m ）流量Q （L/s ）单位流量q （L/s•m ）渗透系数 （m/d ）d T 1 96.66 139.34 27.18 0.00669 0.00652 0.0160 c P 2 108.96 126.85 35.98 0.0113 0.0114 0.0000561 w P 2 131.96 121.27 50.19 0.0110 0.0112 0.0000202 m P 1120.33 127.70 34.72 0.0129 0.0148 0.0147根据不同标高节理裂隙统计计算的渗透张量 表3有限元分析模型围岩渗透参数分区表 表4围岩位置标高（m ）渗透系数（m/d ）第一层 800-900 0.9443第二层 700-800 0.4372第三层 600-700 0.2096第四层 480-600 0.10243.4 有限元模型和计算工况划分毛坝向斜段隧道衬砌为复合式抗水压衬砌，采用圆形断面型钢混凝土衬砌。

有限元模型和网格划分见图4。

分析工况按隧道是否设排水系统、注浆圈的渗透系数、厚度、衬砌渗透系数等分表5中46种工况。

表5中，围岩采用等效连续介质模拟，工况1-1～1-31注浆圈、衬砌采用各项同性均匀介质模拟；工况2-1～2-8中衬砌看成是不透水材料，衬砌背后通过设置渗透系数较大的透水材料模拟排水系统；工况3-1～3-7中也将衬砌看成不透水材料，在排水孔位置设置渗透系数很大的材料模拟。

（a ）整体 （b ）隧道周围图4 有限元模型和网格划分等效连续介质有限元分析工况表 表5工况编号 排水型式注浆圈 厚度L g /m注浆圈渗透 系数K g /cm/S衬砌渗透 系数K l /cm/S 1-1 水通过衬砌渗流进隧道 0 —— 4.6×10-81-2 水通过衬砌渗流进隧道 3 4.60×10-7 4.6×10-81-3 水通过衬砌渗流进隧道 5 4.60×10-4 4.6×10-41-4 水通过衬砌渗流进隧道 5 4.60×10-4 4.6×10-51-5 水通过衬砌渗流进隧道 5 4.60×10-4 4.6×10-61-6 水通过衬砌渗流进隧道 5 4.60×10-4 4.6×10-71-7 水通过衬砌渗流进隧道 5 4.60×10-4 4.6×10-81-8 水通过衬砌渗流进隧道 5 4.60×10-4 4.6×10-91-9 水通过衬砌渗流进隧道 5 4.60×10-4 4.6×10-101-10 水通过衬砌渗流进隧道 5 4.60×10-5 4.6×10-41-11 水通过衬砌渗流进隧道 5 4.60×10-5 4.6×10-51-12 水通过衬砌渗流进隧道 5 4.60×10-5 4.6×10-61-13 水通过衬砌渗流进隧道 5 4.60×10-5 4.6×10-71-14 水通过衬砌渗流进隧道 5 4.60×10-5 4.6×10-81-15 水通过衬砌渗流进隧道 5 4.60×10-5 4.6×10-91-16水通过衬砌渗流进隧道5 4.60×10-54.6×10-10续上工况编号排水型式注浆圈厚度L g/m注浆圈渗透系数K g/cm/S衬砌渗透系数K l /cm/S1-17 水通过衬砌渗流进隧道 5 4.60×10-6 4.6×10-41-18 水通过衬砌渗流进隧道 5 4.60×10-6 4.6×10-51-19 水通过衬砌渗流进隧道 5 4.60×10-6 4.6×10-61-20 水通过衬砌渗流进隧道 5 4.60×10-6 4.6×10-71-21 水通过衬砌渗流进隧道 5 4.60×10-6 4.6×10-81-22 水通过衬砌渗流进隧道 5 4.60×10-6 4.6×10-91-23 水通过衬砌渗流进隧道 5 4.60×10-6 4.6×10-101-24 水通过衬砌渗流进隧道 5 4.60×10-7 4.6×10-41-25 水通过衬砌渗流进隧道 5 4.60×10-7 4.6×10-51-26 水通过衬砌渗流进隧道 5 4.60×10-7 4.6×10-61-27 水通过衬砌渗流进隧道 5 4.60×10-7 4.6×10-71-28 水通过衬砌渗流进隧道 5 4.60×10-7 4.6×10-81-29 水通过衬砌渗流进隧道 5 4.60×10-7 4.6×10-91-30 水通过衬砌渗流进隧道 5 4.60×10-7 4.6×10-101-31 水通过衬砌渗流进隧道8 4.60×10-7 4.6×10-82-1 水通过排水系统进隧道0 —— —— 2-2 水通过排水系统进隧道 3 4.60×10-7——2-3 水通过排水系统进隧道 5 4.60×10-4——2-4 水通过排水系统进隧道 5 4.60×10-5——2-5 水通过排水系统进隧道 5 4.60×10-6——2-6 水通过排水系统进隧道 5 4.60×10-7——2-7 水通过排水系统进隧道 5 1.30×10-6——2-8 水通过排水系统进隧道8 4.60×10-7——3-1 仅在衬砌墙脚位置设排水孔0 —— —— 3-2 仅在衬砌墙脚位置设排水孔 3 4.60×10-7——3-3 仅在衬砌墙脚位置设排水孔 5 4.60×10-4——3-4 仅在衬砌墙脚位置设排水孔 5 4.60×10-5——3-5 仅在衬砌墙脚位置设排水孔 5 4.60×10-6——3-6 仅在衬砌墙脚位置设排水孔 5 4.60×10-7——3-7 仅在衬砌墙脚位置设排水孔8 4.60×10-7——3.5 计算结果及分析各工况计算得到的隧道控制排水量、注浆圈外表面、衬砌背后水压力作用系数结果见表6。