有压引水隧洞多层柔性叠合衬砌对断层蠕滑错动的适应性研究

徐文韬; 徐红超; 石长征; 伍鹤皋; 张超; 李端正

doi:10.11779/CJGE20230103

有压引水隧洞多层柔性叠合衬砌对断层蠕滑错动的适应性研究 English Version

1.
武汉大学水资源工程与调度全国重点实验室, 湖北武汉 430072
2.
云南省滇中引水工程有限公司, 云南昆明 650051

基金项目:

云南省重大科研专项计划项目 202102AF080001

详细信息

作者简介:
徐文韬(1996—)，男，博士研究生，主要从事地下工程与压力管道方面的研究工作。Email：xwt_whu@whu.edu.cn

通讯作者:
伍鹤皋, E-mail：wbf1988@vip.sina.com

中图分类号: TU43
计量
- 文章访问数: 237
- HTML全文浏览量: 23
- PDF下载量: 70
出版历程
- 收稿日期: 2023-02-09
- 网络出版日期: 2024-05-14
- 刊出日期: 2024-04-30

Adaptation of multi-layer flexible stacked linings of pressure water transmission tunnels subjected to fault creep

1.
State Key Laboratory of Water Resources Engineering and Management, Wuhan University, Wuhan 430072, China
2.
Central Yunnan Provincial Water Diversion Project Co., Ltd., Kunming 650051, China

摘要

摘要: 长距离输水隧洞工程建设很难避免穿越活动断层的问题，有压隧洞采用何种衬砌结构型式穿越活动断层迄今为止研究还不够深入。针对目前隧洞过活断层常用抗断措施存在的问题，结合有压引水隧洞自身的特点，提出了一种适合有压引水隧洞穿越活动断层的多层柔性叠合衬砌结构型式。依托某工程的实践情况，采用有限元数值分析方法论证了多层柔性叠合衬砌结构的可行性和合理性；同时分析了垫层厚度、垫层弹模对隧洞结构内力响应和抗断性能的影响。计算结果表明：垫层厚度越大，伸缩节总变形量越大，对于隧洞衬砌适应断层错动越有利；但当垫层厚度过大时，将对衬砌的受力产生不利影响，因此建议根据工程设防的断层错动量选择合适厚度的垫层，不必采用过厚的垫层。垫层弹模的变化对钢管应力的影响较为显著，综合考虑钢衬应力、混凝土损伤以及对断层错动的适应性能，建议选取3~5 MPa的垫层弹模最为合适。研究结果可为有压输水隧洞应对断层蠕滑变形时的结构设计提供参考。
- 有压隧洞 /
- 衬砌结构 /
- 断层错动 /
- 柔性接头 /
- 垫层 /
- 数值模拟
Abstract: Crossing active faults is a difficult problem to avoid when building long-distance water transmission tunnels, and the structural form of pressure tunnels crossing active faults has not been thoroughly studied. In this study, a multi-layer flexible stacked lining suitable for pressure water transmission tunnels crossing active faults is proposed aiming at addressing the problems of the common anti-fracture measures for the tunnels crossing active faults as well as their characteristics. Based on the actual situation of a project, the finite unit method is used to analyze the stress and deformation of the structures under the condition of using different tunnel lining structures in the project. According to the calculated results, the reasonableness of the multi-layer flexible stacked lining structures is confirmed. The numerical simulation is used to assess the impact of cushion thickness and elastic modulus on the internal force response and fracture resistance of the tunnel structures. The findings demonstrate that the thicker the cushion thickness the larger the total deformation of the expansion joints. It is more favorable for the tunnel lining to adapt to fault dislocation. However, when the cushion thickness increases further, there is a certain adverse effect on the force of the lining, so it is recommended to choose the appropriate thickness of the cushion layer according to the amount of fault dislocation. It is not necessary to use a too-thick cushion layer. The change of elastic modulus of cushion has a significant effect on the stress of steel pipes. Considering the stress of steel lining, the damage of concrete and the adaptability to fault dislocation, it is recommended that the elastic modulus of the cushion be selected as 3 to 5 MPa. The research results can provide some references for the design of pressure water transmission tunnels crossing fault zones.
- pressure tunnel /
- lining structure /
- fault creep /
- flexible joint /
- cushion /
- numerical simulation

HTML全文

0. 引言

近年来,盾构法应用越来越广泛,尤其在跨江隧道和城市地铁建设中。但是随着大直径盾构隧道的广泛应用,施工期管片上浮问题普遍而突出,过大的上浮量会导致管片错台、环缝张开、环间螺栓剪断和渗漏水等安全事故的发生,严重影响隧道成型质量。因此为了保证隧道施工与运营的安全,有必要对盾构隧道施工期管片的上浮机理进行研究,以及时采取科学合理的应对性措施。

盾构隧道的上浮问题是一个复杂的三维力学问题^[1],对于其上浮的原因,已有研究^[2]表明一般需要具备两大必要条件：①管片上浮变形发生的空间,②引起管片上浮的力。国内外学者从横向和纵向上分别对施工期盾构隧道的上浮问题进行研究,且已取得一定的成果。盾构隧道横向（局部）上浮机理的理论研究日益趋向成熟,主要包括均质圆环法^[3]、修正惯用法^[4-5]、多铰圆环法^[6]、梁-弹簧模型法^[7]、梁–接头模型法^[8]等。近年来盾构隧道纵向性能的研究受到越来越多的关注,但其研究目前尚处于起步阶段,主要有小泉淳^[7]提出的纵向梁–弹簧模型和志波尤纪夫^[9]提出的纵向等效连续化梁模型,之后的学者基于这两种理论做了大量的工作,叶飞等^[10]基于等效连续梁理论,将隧道和土体的相互作用等效为土弹簧,研究了同步注浆产生的不同形式上浮力作用下隧道结构的纵向变形特性。但是该模型没有考虑同步注浆过程中浆液硬化的时变性和上覆土体基床系数的各异性。朱令等^[11]、杨方勤^[12]、舒瑶等^[13]等考虑了同步注浆过程中浆液黏度时变性及上覆土体基床系数各异性的影响,但需要通过有限元数值模拟实现,难以被施工和设计人员接受及应用。

在前人研究的基础上,考虑浆液黏度时变性、上覆土体基床系数各异性以及施工步荷载叠加效应的影响,基于弹性地基梁的矩阵传递法理论,建立施工期盾构隧道纵向上浮分析模型,并进行了试验验证。

1. 隧道纵向上浮模型的构建

在盾尾进行同步注浆时,盾构隧道管片的受力为一个复杂的三维问题,涉及到地层-浆液-管片的相互作用,但是目前大多数的纵向分析都不考虑横向变形的影响,将隧道的三维受力简化为纵向的一维问题进行求解,最为经典的是日本学者志波由纪夫提出的纵向等效连续梁模型,由于该模型概念清楚、计算简单,近年来被广泛应用,本文的隧道纵向上浮模型的也是在等效连续梁模型的基础上构建,如图1所示。该模型可以分成两部分,第一部分为弹性地基上长度为L的有限长梁,在该段长度范围内,鉴于浆液压力的耗散和浆液黏度的时变特性^[14],假定上浮力呈现出三角形线性变化。随着浆液的逐渐硬化,上覆土体基床系数也从零逐渐增大,最后在距盾尾约为L处趋于一个稳定的值K,一般认为在该段范围内,基床系数近似于三角形线性变化。浆液未凝固区长度L根据隧道受到上浮力作用范围来确定。第二部分为离盾尾更远的弹性地基上的半无限长梁^[15-16],该段起点为距盾尾为L处,从该处起认为浆液固结硬化完成,因此基床系数不再发生变化,此外浆液压力已耗散至零,半无限长梁不受上浮力的作用。弹性地基有限长梁的左端认为是铰接以模拟盾尾处对管片的约束^[12],半无限长梁右端无穷远处变形、内力为零。有限长梁和半无限长梁在x=L处内力、变形连续,基于上述边界条件和矩阵传递法理论,求解上述耦合梁模型的内力和变形。

图 1 盾构隧道纵向上浮分析模型

Figure 1. Longitudinal uplift analysis model for shield tunnel

下载: 全尺寸图片幻灯片

盾构隧道的施工为一个动态的过程,随着盾构机的向前掘进一环,则同步注浆产生的上浮力向前移动一环,地基基床系数变化区间也向前移动一环。如此循环往复可模拟盾构机不断向前开挖的过程。施工步示意图如图2所示。

图 2 施工动态示意图

Figure 2. Diagram of dynamic construction

下载: 全尺寸图片幻灯片

2. 模型参数确定

2.1 等效地层抗力系数

假定盾构隧道同步注浆材料及地层为弹性体,半无限长梁区域中的等效地层抗力系数示意图如图3所示。

图 3 等效地层抗力系数示意图

Figure 3. Stiffness coefficient of equivalent stratigraphic spring

下载: 全尺寸图片幻灯片

利用Muir Wood理论解^[17],计算注浆层与地层共同作用下的等效地层抗力系数：

$k = \frac{3 E_{0}}{(1 + ν) (5 - 6 v) R c}$ ,

(1)

$\begin{matrix} E_{0} = | - \frac{1}{D_{c} + 2 H \tan θ} + \frac{1}{D_{c}} | / [\frac{1}{E_{0b}} (- \frac{1}{D_{c} + 2 H_{b} \tan θ} + \frac{1}{D_{c}}) + \\ \frac{1}{E_{0g}} (- \frac{1}{D_{c} + 2 H_{g} \tan θ} + \frac{1}{D_{c} + 2 H_{b} \tan θ})] \end{matrix}$ ,

(2)

式中, $ν$ 为泊松比, $R_{c}$ 为管片形心线的半径, $D_{c}$ 为管片形心线的直径, $D_{0}$ 为管片外径, $E_{0}$ 为考虑同步注浆影响的等效地层变形模量, $E_{ob}$ 为壁后注浆层的变形模量, $E_{og}$ 为靠近注浆层的地层变形模量, $H$ 为注浆影响范围, $H_{b}$ 为同步注浆层的厚度, $H_{g}$ 为影响范围内地层的厚度, $θ$ 为荷载的分布宽度（一般取值30°）。

2.2 等效抗弯刚度

等效刚度法将每环衬砌环视为均质圆环,在纵向上基于刚度等效的思想把由管片和接头构成的盾构等效为具有相同结构刚度和特性的均质连续梁,其等效刚度大小^[18-19]]为

$E I_{eq} = \frac{\cos^{3} φ}{\cos φ + (\frac{π}{2} + φ) \sin φ} E_{c} I_{c}$ ,

(3)

$φ + \cot φ = π (\frac{1}{2} + \frac{n K_{b} l}{E_{c} A_{c}})$ 。

(4)

式中 $φ$ 为弯曲后中性轴的位置角； $K_{b}$ 为连接螺栓的线刚度, $K_{b} = E_{b} A_{b} / l_{b}$ , $E_{b}$ 为螺栓弹性模量, $A_{b}$ 为螺栓横截面积, $l_{b}$ 为螺栓计算长度, $n$ 为环间连接螺栓个数； $I_{c}$ 为隧道纵向惯性矩； $E_{c}$ 为管片弹性模量； $A_{c}$ 为隧道横截面面积； $l_{c}$ 为相邻两管片环中心线间的距离,其大小近似于管片环宽 $B$ 。

2.3 浆液未凝固区长度

同步注浆中浆液的凝固硬化为浆液压力逐渐耗散的复杂物理—化学过程,等效连续梁模型中浆液未凝固区长度主要和同步注浆工艺、浆液性质、地层特性以及注浆压力等因素有关。以往的研究^[11-12,15]中,浆液未凝固区的长度取值多基于现场的经验,并未给出合理的取值依据。梁禹等^[14]考虑浆液压力耗散、黏度时变性和地层特性等的影响,得到盾构隧道纵向浆液压力分布规律。浆液纵向扩散的最远距离可表征为浆液未凝固区长度,可表示为

$L = \frac{P_{G 0}}{B}$ ,

(5)

$P_{G 0} (t) = P_{0} - \frac{2 G}{r} \frac{1}{1 - ν} \frac{n_{i} - n_{e}}{1 - n_{i}} x (t)$ 。

(6)

式中 $G$ 为围岩土体剪切模量； $r$ 为围岩扰动层厚度； $ν$ 为扰动层泊松比； $n_{i}$ 为浆液孔隙比； $n_{e}$ 为浆液凝固硬化后浆饼孔隙率；通过室内试验获得, $P_{0}$ 为初始注浆压力, $P_{G0}$ 为考虑浆液压力耗散效应的初始浆液压力。 $x (t)$ 为浆饼的厚度,浆饼的厚度通过以下浆液固结方程求得

$(\frac{x}{k} + \frac{F}{k_{s} (t)}) \frac{d x}{d t} + C x = D h_{0}$ ,

(7)

式中, $x$ 为浆饼厚度, $F$ , $C$ , $D$ 为与隧道和浆液特性相关的系数, $k$ 为浆饼的渗透系数, $k_{s} (t)$ 为浆液向地层的渗透系数。考虑浆液黏度时变性影响时, $k_{s} (t)$ 可表示为

$k_{s} (t) = \frac{k_{w}}{β_{0}} e^{- ξ t}$ ,

(8)

式中, $β_{0}$ 为浆液初始黏度与水的黏度比, $β_{0} = u_{0} / u_{w}$ 。 $k_{w}$ 为水在地层中的渗透系数, $ξ$ 为跟浆液孔隙率相关的参数,通过室内试验获得。

2.4 浆液最大上浮力

最大上浮力取值为浆液压力及管片自重等荷载的合力,叶飞等^[10]提出了上浮力的分布形式近似为扇形分布。衬砌环受压密注浆作用,注浆浆液由于重力作用在管片下方聚集,形成较大的集中力或分布力,压密阶段注浆压力形成的向上的合力：

$F = \int_{- θ}^{θ} B P R_{0} \cos α d α = 2 B P R_{0} \sin θ (0 \leq θ \leq \frac{π}{2})$ 。

(9)

式中 $B$ 为浆液压力作用宽度,一般取1 m； $P$ 为该处管片所受浆液压力平均值； $R_{0}$ 为管片外径； $θ$ 为注浆浆液泡分布区域边界与竖向的夹角,一般取 $π / 2$ 。

3. 基于矩阵传递法的模型求解

将施工期的盾构隧道纵向上浮模型在x=L处分成两部分,左边为基床系数线性变化的有限长梁,梁的左端视为铰接,受线性荷载作用。右边为不受上浮力作用、基床系数固定的半无限长梁。对于这类复杂变基床系数的耦合弹性地基梁问题,目前的研究多基于有限元数值模拟手段对其进行求解,还没有形成更具体的理论对其进行分析。本章基于矩阵传递法理论^[20]对变基床系数条件下的弹性地基梁问题求解,克服了这两类地基梁计算方法不相同的困难,力学概念清晰、明确。

刘庆潭等^[20]通过拉普拉斯变换得到初参数解的形式：

(10)

式中, $v$ 为挠度, $θ$ 为转角, $M$ 为弯矩, $Q$ 为剪力, $K$ 为基床系数, $q$ 为线性荷载初值, $Δ q / Δ L$ 为线性荷载斜率。 $φ_{1}$ , $φ_{2}$ , $φ_{3}$ , $φ_{4}$ 为双曲函数和三角函数的组合形式,表示以下形式：

$\begin{array}{l} φ_{1} = c o h (β x) \cos (β x), \\ φ_{2} = c o h (β x) \sin (β x) + \sinh (β x) \cos (β x), \\ φ_{3} = c o h (β x) \cos (β x), \\ φ_{4} = c o h (β x) \sin (β x) - \sinh (β x) \cos (β x), \end{array}}$

(11)

式中, $β$ 为地基梁特征系数,与地基梁抗弯刚度和基床系数有关。这种初参数解的形式也称为矩阵传递法,该方法对处理弹性地基上有限长梁的问题具有十分明显的优势。为了接下来方便计算,令式（11）的传递矩阵为 $T$ ,荷载影响项矩阵为 $F q$ 。

分析左半段弹性地基梁有限长梁变基床系数问题时,基于面积等效的思想将线性变化的三角形地基等效成若干段如图4所示,有限长梁左端的约束视为铰接,因此4个初参数的 $w_{0}$ , $M_{0}$ 被指定,剩下两个初参数 $θ_{0}$ , $Q_{0}$ 待求,参数的传递从左边铰接端开始,矩阵传递的过程可表示为

图 4 三角形基床转化为阶梯形基床

Figure 4. Conversion of triangular subgrade bed into stepped one

下载: 全尺寸图片幻灯片

$\begin{array}{l} S_{1} = T_{0} S_{0} + F q_{0}, \\ S_{2} = T_{1} S_{1} + F q_{1}, \\ S_{3} = T_{2} S_{2} + F q_{2}, \\ ... \\ S_{L} = T_{L - 1} S_{L - 1} + F q_{L - 1}, \end{array}}$

(12)

式中,S_n为第n段梁右端的内力矩阵,T_n-1为第n段梁的传递矩阵,Fq_n-1为第n段梁的荷载附加项矩阵。

从式（12）可以看出, $x = L$ 处的变形及内力可以表示为 $θ (0)$ , $Q (0)$ 的线性表达,再结合半无限长梁的无穷远端的边界条件就可以求解出 $θ (0)$ , $Q (0)$ 。结合式（11）,半无限长梁任意处的竖向位移可表示为

$\begin{array}{l} v = φ_{1} w_{L} - \frac{1}{2 β} φ_{2} θ_{L} - \frac{2 β^{2}}{K} φ_{3} M_{L} - \frac{β}{K} φ_{4} Q_{L} + \\ \frac{q}{K} (1 - φ_{1}) + \frac{Δ q}{K Δ L} (x - \frac{φ_{2}}{2 β}) \end{array}$ 。

(13)

引入4个衰减函数：

$\begin{array}{l} W_{1} = e^{- β x} \cos (β x), \\ W_{2} = e^{- β x} (\sin (β x) - \cos (β x)), \\ W_{3} = e^{- β x} \sin (β x), \\ W_{4} = e^{- β x} (\sin (β x) + \cos (β x)), \end{array}}$

(14)

4个增长函数：

$\begin{array}{l} U_{1} = e^{β x} \cos (β x), \\ U_{2} = e^{β x} (\sin (β x) - \cos (β x)), \\ U_{3} = e^{β x} \sin (β x), \\ U_{4} = e^{β x} (\sin (β x) + \cos (β x)), \end{array}}$

(15)

衰减函数、增长函数具有以下特性：

$\begin{array}{l} W (x = + \infty) = 0, \\ U (x = + \infty) = \infty 。 \end{array}}$

(16)

同时 $φ_{1}$ , $φ_{2}$ , $φ_{3}$ , $φ_{4}$ 和衰减函数、增长函数之间有着如下的关系：

$\begin{array}{l} φ_{1} = \frac{U_{1}}{2} + \frac{W_{1}}{2}, \\ φ_{2} = \frac{U_{4}}{2} + \frac{W_{2}}{2}, \\ φ_{3} = \frac{U_{3}}{2} - \frac{W_{3}}{2}, \\ φ_{4} = \frac{W_{4}}{2} + \frac{U_{2}}{2} 。 \end{array}}$

(17)

将式（17）代入式（13）并且结合无穷远处衰减函数的性质可以得到

$\begin{array}{l} v = \frac{U_{1}}{2} v_{L} - \frac{U_{4}}{4 β} θ_{L} - \frac{β^{2} U_{3}}{K} M_{L} - \frac{β U_{2}}{2 K} Q_{L} + \\ \frac{q}{K} (1 - \frac{U_{1}}{2}) + \frac{Δ q}{Δ L K} (x - \frac{U_{4}}{4 β}) \end{array}$ 。

(18)

将 $U_{1}$ , $U_{2}$ , $U_{3}$ , $U_{4}$ 按式（11）的形式代入式（18）可以得

$\begin{array}{l} v = e^{β x} (S_{1} \cos β x + S_{2} \sin β x) + \frac{q}{K} + \frac{Δ q}{Δ L K} x, \\ S_{1} = \frac{v_{L}}{2} - \frac{θ_{L}}{4 β} + \frac{β}{2 K} Q_{L} - \frac{q}{2 K} - \frac{Δ q}{4 K β Δ L}, \\ S_{2} = - \frac{θ_{L}}{4 β} - \frac{β^{2} M_{L}}{K} - \frac{β Q_{L}}{2 K} - \frac{Δ q}{4 K β Δ L} 。 \end{array}}$

(19)

由于半无限长弹性地基梁段不受上浮力作用,因此式（19）的 $q$ , $Δ q / Δ L$ 的值为0。根据半无限长梁的性质可知无穷远处的挠度 $v$ 等于0,由于 $x \to + \infty$ , $\sin α x$ 和 $\cos α x$ 不恒为0,所以要使式（19）中 $v$ 的收敛为0,必须使 $\sin α x$ 和 $\cos α x$ 的系数为零。因此可以得到 $S_{1}$ 和 $S_{2}$ ,可表示为

$\begin{matrix} \frac{v_{L}}{2} - \frac{θ_{L}}{4 β} + \frac{β}{2 K} Q_{L} = \frac{q}{2 K} + \frac{Δ q}{4 K β Δ L} - \\ \frac{θ_{L}}{4 β} - \frac{β^{2} M_{L}}{K} - \frac{β Q_{L}}{2 K} = \frac{Δ q}{4 K β Δ L} \end{matrix}$ 。

(20)

式（20）为半无限长梁的边界条件,从式（10）知 $v_{L}$ , $θ_{L}$ , $M_{L}$ , $Q_{L}$ 为 $θ_{0}$ , $Q_{0}$ 的线性表达,所以求解上述二元一次方程组即可求解出 $θ_{0}$ , $Q_{0}$ ,从而可以得到各段梁的变形和内力理论解。

4. 工程验证

长沙市南湖路隧道采用盾构法施工,盾构隧道埋深为10.79 m,采用双管单层形式,分为南北两线穿越湘江,工程位置示意图如图5所示。选取的典型横断面如图6所示,地层从上到下依次为淤泥质土、强风化砾岩。

图 5 越江隧道工程位置图

Figure 5. Schematic diagram of tunnel across Xiangjiang River

下载: 全尺寸图片幻灯片

图 6 典型断面示意图

Figure 6. Schematic diagram of typical section

下载: 全尺寸图片幻灯片

4.1 参数选取

隧道管片外径5.65 m,管片内径5.15 m,管片厚度0.5 m,管片幅宽2 m,环向分块9块（即由6个标准块、2个邻接块和1个封顶块组成）,衬砌环采用错缝的方式拼接,每环选用12根M24螺栓,环缝选用10根M24螺栓。盾构机外径5.815 m,盾构注浆层厚度0.165 m,初始注浆压力为0.18 MPa,采用典型4孔注浆。模型的主要参数见表1,2所示。

表 1 浆液主要参数

Table 1. Key parameters of grout

初始注浆压力/MPa	注浆层厚度/m	浆饼渗透系数/(m·s^-1)	浆液孔隙率	浆饼孔隙率	浆体变形模量/MPa
0.18	0.165	4.7×10^-8	0.425	0.417	15

下载: 导出CSV

| 显示表格

表 2 地基土及混凝土部分参数

Table 2. Partial parameters of foundation soil and concrete

地层	混凝土密度/(kg·m^-3)	土体变形模量/MPa	土体泊松比	土体渗透系数/(m·d^-1)
强风化砾岩	2500	18.5	0.3	1.85×10^-6

下载: 导出CSV

| 显示表格

4.2 模型计算结果

（1）浆液纵向压力分布规律

基于式（6）可以得到同步注浆作用下浆液压力耗散规律,如图7所示,注浆一般持续2～3 h,选择 $t$ =7500 s与 $t$ =0 s时浆液压力进行对比。从图7可知,考虑浆液黏度变化时 $t$ =7500 s的浆液压力相对于 $t$ =0 s的浆液压力衰减10%,不考虑浆液黏度变化时 $t$ =7500 s的浆液压力相对于 $t$ =0 s时的浆液压力衰减27.2%,因此考虑浆液黏度时变性的影响时,最终的衰减幅度更小。浆液黏度增大,其流动性和地层的渗透性变小,进而会影响浆液压力耗散。

图 7 浆液压力耗散曲线

Figure 7. Dissipation curves of grout pressure

下载: 全尺寸图片幻灯片

现场对注浆口附近浆液压力变化进行了监测,结果如图8所示,通过对比浆液设定压力和实测浆液压力的数据可知,注浆口处浆液压力存在一定的衰减,衰减幅度大致在16.6%～24.5%。现场的实测结果表明了纵向浆液扩散衰减理论的合理性,因此在分析同步注浆下盾构隧道纵向的力学行为时,应考虑浆液耗散效应带来的影响。

图 8 浆液压力现场实测

Figure 8. Field measurement of groug pressure

下载: 全尺寸图片幻灯片

基于浆液压力耗散曲线计算结果,可以得到同步注浆过程中浆液压力沿纵向的扩散规律及浆液未凝固区长度,如图9,10所示,浆液压力沿纵向的分布呈线性衰减趋势,若不考虑浆液黏度时变性的影响,同步注浆过程中浆液压力的耗散幅度为0.049 MPa,占初始预设浆液压力的27.7%,浆液未凝固区长度为7.53 m。若考虑浆液黏度时变性的影响,浆液压力耗散幅度为0.018 MPa,占初始预设浆液压力的10%,浆液未凝固区长度为9.32 m,相对于不考虑浆液黏度时变性影响时,其长度增加1.79 m。

图 9 不考虑浆液黏度时变性的纵向浆液压力分布

Figure 9. Distribution of longitudinal pressure without consideration of time-varying property of grout viscosity

下载: 全尺寸图片幻灯片

图 10 考虑浆液黏度时变性的纵向浆液压力分布

Figure 10. Distribution of longitudinal pressure with consideration of time-varying property of grout viscosity

下载: 全尺寸图片幻灯片

（2）典型断面预测值

基于式（1）可以计算得到壁后注浆与地层共同作用下的等效地层抗力系数 $k$ =3.64×10² kN/m³,地层抗力系数乘以隧道宽度为基床系数K=4.11×10³ kN/m²,基于式（2）可以计算得到弹性地基上等效连续梁的等效抗弯刚度(EI)_eq=9.08×10⁸kN·m²。

根据前文得出的纵向浆液压力分布规律和式（9）,再减去管片自重可以得到纵向上浮力分布形式：

$\begin{array}{l} F_{1} (x) = 3643.94 - 393.24 x, \\ F_{2} (x) = 2536.51 - 393.24 x, \\ F_{3} (x) = 3227.18 - 393.24 x 。 \end{array}}$

(21)

式（21）分别为不考虑浆液压力耗散（下文简称为工况F₁）、考虑浆液压力耗散但不考虑浆液黏度时变性（工况F₂）、考虑浆液压力耗散和浆液黏度时变性（工况F₃）的3种盾构隧道纵向上浮力分布形式。将上述3种形式的上浮力作为盾构隧道纵向上浮分析模型的输入荷载,代入盾构隧道纵向上浮分析模型,可以得到一个施工步下的上浮量如图11所示,从盾构铰接端到无穷远处,上浮量先增大后减少,在6～7环（约12 m）达到最大,后迅速下降,约在20环（40 m）处降为零,随后出现一段数值较小的负位移,最后约40环（80 m）左右管片的位移重新为0,之后一直保持不变。工况F₁的上浮量最大,为3.3 mm；工况F₃的上浮量次之,为2.4 mm；工况F₂的上浮量最小,为1.7 mm。这是因为不考虑浆液压力耗散效应时,初始注浆压力最大,浆液未凝固区最长,同步注浆上浮力作用的范围最长；而工况F₃和F₂相比,浆液未凝固区长度更长,上浮力作用范围更长。

图 11 隧道纵向上浮量

Figure 11. Longitudinal uplift amount of tunnel

下载: 全尺寸图片幻灯片

盾构机实际施工时不断向前掘进,还应考虑施工步的影响时,将典型断面对应的不同施工步下的上浮曲线进行叠加,得到目标环的上浮量变化规律,如图12所示,典型断面对应的目标环脱出盾构后,上浮量持续增加,约向前掘进40环（80 m）后目标环的上浮量趋于稳定,不再变化。工况F₁的最终上浮量最大,为36.1 mm；工况F₃次之,最终上浮量为25.3 mm；工况F₂最终上浮量最小,为17.6 mm。

图 12 典型断面累计上浮量

Figure 12. Cumulative uplift amount of a typical section

下载: 全尺寸图片幻灯片

现场对典型断面的目标环拱顶和拱底的上浮量进行了监测,现场实测数据如图13所示,目标环脱出盾尾以后,第一天上浮量达15 mm。随后上浮量增长速率下降,最终约10 d后保持不变,拱顶的最终上浮量为28.7 mm,拱底的最终上浮量为26.8 mm,取两者的平均值（27.75 mm）作为目标环的最终上浮量。通过和前面3种情况下的对比可以发现,同时考浆液压力耗散效应和浆液黏度时变性影响的最终上浮量和现场实测数据最为接近,误差为8.8%,其他两种工况下的上浮量预测值误差较大,分别为30%和36.6%。

图 13 典型断面现场实测上浮量

Figure 13. Measured uplift amount of a typical section

下载: 全尺寸图片幻灯片

5. 结论与展望

基于弹性地基梁的矩阵传递法理论,综合考虑了浆液压力耗散效应、浆液黏度时变性、上覆土体基床系数各异性和施工步叠加效应的影响,建立施工期盾构隧道纵向上浮分析模型,并基于工程案例实测结果进行了验证,可得出以下结论：

（1）模型将施工期的盾构隧道视为变基床系数的有限长弹性地基梁和基床系数不变的半无限弹性地基长梁组成,通过初参数矩阵传递的形式,实现了这两类地基梁计算方法的耦合,同时解决了考虑浆液压力耗散效应、浆液黏度时变特性、基床系数各异性和施工步叠加效应等多因素影响下对上浮量精细化预测的问题。

（2）以长沙市南湖路隧道工程为案例,进行了实测数据与计算值的对比,发现计算结果和现场实测较为吻合,同时在盾构隧道施工过程中,同步注浆作用下引起的管片上浮量从盾尾处到无穷远先增大后减小,距盾尾约6～7环处上浮量最大,之后上浮量下降,远离盾尾超40环的管片不再上浮。研究成果可为后续类似盾构隧道上浮量的控制提供参考。

（3）建立的施工期盾构隧道纵向上浮分析模型没有考虑管片接头的弯曲非线性特性,始终认为纵向等效抗弯刚度为固定值,而纵向等效抗弯刚随弯矩、轴力变化呈现非线性动态变化的特性,因此为了能反映隧道纵向更为精确的内力和变形分布,还需建立更为精细化的上浮分析模型。

图 1 近场区主要断裂带延伸示意图

Figure 1. Schematic diagram of extension of main fracture zone

下载: 全尺寸图片幻灯片

图 2 两种隧洞衬砌结构示意图

Figure 2. Schematic diagram of two kinds of tunnel lining structures

下载: 全尺寸图片幻灯片

图 3 有限元整体模型图

Figure 3. Overall finite element model

下载: 全尺寸图片幻灯片

图 4 混凝土材料本构参数

Figure 4. Parameters of constitutive model

下载: 全尺寸图片幻灯片

图 5 边界条件示意图

Figure 5. Boundary conditions

下载: 全尺寸图片幻灯片

图 6 断层错动情况下钢管合位移图

Figure 6. Resultant displacements of pipes under fault condition

下载: 全尺寸图片幻灯片

图 7 不同断层错动量下钢管的应力变化曲线

Figure 7. Curves of stress of steel pipes under different fault displacement dislocations

下载: 全尺寸图片幻灯片

图 8 不同断层错动量下钢管Mises应力分布云图

Figure 8. Distribution of stress of steel pipes under different fault displacement dislocations

下载: 全尺寸图片幻灯片

图 9 混凝土压缩损伤图

Figure 9. Compressive damages of concrete

下载: 全尺寸图片幻灯片

图 10 断层错动情况下钢管合位移图

Figure 10. Resultant displacements of pipes under fault conditions

下载: 全尺寸图片幻灯片

图 11 管道断面中心位移沿管轴线的分布曲线

Figure 11. Distribution curves of central displacement of pipeline sections along pipe axis

下载: 全尺寸图片幻灯片

图 12 垫层径向变形

Figure 12. Radial deformations of cushion

下载: 全尺寸图片幻灯片

图 13 不同衬砌型式钢管的应力变化曲线

Figure 13. Curves of stress of steel pipes with different lining types

下载: 全尺寸图片幻灯片

图 14 不同衬砌型式混凝土的相对损伤比

Figure 14. Relative damage ratios of concrete for different lining types

下载: 全尺寸图片幻灯片

图 15 不同垫层厚度钢管的Mises应力变化曲线

Figure 15. Curves of Mises stress for steel pipes with different cushion thicknesses

下载: 全尺寸图片幻灯片

图 16 伸缩节变形量随垫层厚度的变化关系

Figure 16. Relationship between deformation of expansion joint and cushion thickness

下载: 全尺寸图片幻灯片

图 17 管节首尾断面顶底轴向位移差曲线

Figure 17. Curves of top-bottom axial displacement difference for head and tail sections of pipe sections

下载: 全尺寸图片幻灯片

图 18 不同垫层厚度混凝土的损伤云图

Figure 18. Damages of concrete with different cushion thicknesses

下载: 全尺寸图片幻灯片

图 19 不同垫层弹模钢管的Mises应力变化曲线

Figure 19. Curves of Mises stress of steel pipes with different cushion elastic moduli

下载: 全尺寸图片幻灯片

图 20 伸缩节变形量随垫层弹模的变化关系

Figure 20. Variation of deformation of expansion joint with change in elastic modulus of cushion

下载: 全尺寸图片幻灯片

表 1 围岩材料参数表

Table 1 Material parameters of surrounding rock

名称	围岩类别
名称	Q6	Q5
天然重度/(g·cm^-3)	2.15~2.19	2.43~2.45
内摩擦角/(°)	22~30	34~41
泊松比	0.34~0.38	0.31~0.34
岩体变形模量/GPa	0.3~1.1	2~4
黏聚力/MPa	0.14~0.26	0.41~0.82

下载: 导出CSV

表 2 Q355钢材材料参数表

Table 2 Material parameters of Q355 steel

钢号	壁厚/mm	允许应力/MPa	屈服强度/MPa
Q355	24	210	345

下载: 导出CSV

表 3 垫层和混凝土材料参数表

Table 3 Material parameters of cushion and concrete

材料	重度/(10^-5 N·mm^-3)	泊松比	弹性模量/MPa	轴心抗压/抗拉强度设计值/MPa
C20混凝土	2.4	0.20	2.55×10⁴	9.6/1.1
柔性连接段	2.4	0.20	128	—
垫层	0.5	0.16	3	—

下载: 导出CSV

表 4 不同衬砌型式内容表

Table 4 Content of different lining types

方案	结构型式
1	钢管、混凝土
2	钢管节段、混凝土、伸缩节
3	钢管、混凝土、软垫层
4	钢管、混凝土节段、柔性连接
5	钢管节段、伸缩节、软垫层、混凝土节段、柔性连接

下载: 导出CSV

参考文献(15)

[1]	杨启贵, 张传健, 颜天佑, 等. 长距离调水工程建设与安全运行集成研究及应用[J]. 岩土工程学报, 2022, 44(7): 1188-1210. doi: 10.11779/CJGE202207002 YANG Qigui, ZHANG Chuanjian, YAN Tianyou, et al. Integrated research and application of construction and safe operation of long-distance water transfer projects[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(7): 1188-1210. (in Chinese) doi: 10.11779/CJGE202207002
[2]	何川, 李林, 张景, 等. 隧道穿越断层破碎带震害机理研究[J]. 岩土工程学报, 2014, 36(3): 427-434. doi: 10.11779/CJGE201403004 HE Chuan, LI Lin, ZHANG Jing, et al. Seismic damage mechanism of tunnels through fault zones[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(3): 427-434. (in Chinese) doi: 10.11779/CJGE201403004
[3]	RUSSO M, GERMANI G, AMBERG W. Design and construction of large tunnel through active faults: a recent application[C]// Proceedings of the International Conference of Tunnelling & Underground Space Use. Istanbul, Turkey, 2002.
[4]	周光新, 盛谦, 崔臻, 等. 走滑断层错动影响下跨活断层铰接隧洞破坏机制模型试验[J]. 岩土力学, 2022, 43(1): 37-50. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202201004.htm ZHOU Guangxin, SHENG Qian, CUI Zhen, et al. Model test of failure mechanism of tunnel with flexible joint crossing active fault under strike-slip fault dislocation[J]. Rock and Soil Mechanics, 2022, 43(1): 37-50. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202201004.htm
[5]	崔光耀, 王李斌, 王明年, 等. 隧道纤维混凝土衬砌抗错断性能模型试验研究[J]. 振动与冲击, 2019, 38(13): 50-56, 80. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ201913009.htm CUI Guangyao, WANG Libin, WANG Mingnian, et al. Model tests for anti-breaking performance of a fiber reinforced concrete tunnel lining[J]. Journal of Vibration and Shock, 2019, 38(13): 50-56, 80. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ201913009.htm
[6]	蒋树屏, 李鹏, 林志. 穿越活动断层区隧道的抗断设计对策[J]. 重庆交通大学学报(自然科学版), 2008, 27(6): 1034-1036, 1041. https://www.cnki.com.cn/Article/CJFDTOTAL-CQJT200806007.htm JIANG Shuping, LI Peng, LIN Zhi. Design strategies of breaking resistance of tunnels crossing active faults zone[J]. Journal of Chongqing Jiaotong University (Natural Science), 2008, 27(6): 1034-1036, 1041. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CQJT200806007.htm
[7]	任兴普, 李晓彬, 李卫功. 洗马河二级赛珠水电站引水隧洞跨越活断层设计[J]. 中国水运(下半月), 2015, 15(2): 167-168, 200. https://www.cnki.com.cn/Article/CJFDTOTAL-ZSUX201502078.htm REN Xingpu, LI Xiaobin, LI Weigong. Design of diversion tunnel crossing active fault in Ximahe second-stage Saizhu hydropower station[J]. China Water Transport, 2015, 15(2): 167-168, 200. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZSUX201502078.htm
[8]	赵坤, 陈卫忠, 赵武胜, 等. 逆断层错动作用下隧道衬砌铰接设计参数研究[J]. 岩石力学与工程学报, 2018, 37(增刊1): 3411-3421. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2018S1031.htm ZHAO Kun, CHEN Weizhong, ZHAO Wusheng, et al. Study on parameters of articulated design of tunnel lining under reverse fault dislocation[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(S1): 3411-3421. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2018S1031.htm
[9]	苏凯, 常智慧, 崔金鹏, 等. 深埋隧洞开挖数值模拟分析的纵向模型范围研究[J]. 岩土力学, 2016, 37(增刊2): 706-714. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2016S2090.htm SU Kai, CHANG Zhihui, CUI Jinpeng, et al. Study of model longitudinal range in numerical simulation of deep tunnel excavations[J]. Rock and Soil Mechanics, 2016, 37(S2): 706-714. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2016S2090.htm
[10]	石长征, 石雅竹, 伍鹤皋, 等. 长距离埋地钢管中波纹管伸缩节的作用研究[J]. 特种结构, 2021, 38(5): 64-68. https://www.cnki.com.cn/Article/CJFDTOTAL-TZJG202105014.htm SHI Changzheng, SHI Yazhu, WU Hegao, et al. Function of bellows expansion joint in long-distance buried steel pipe[J]. Special Structures, 2021, 38(5): 64-68. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TZJG202105014.htm
[11]	史鑫宇, 姚燕, 王玲, 等. 基于单轴拉压模拟的CDP模型参数影响[J]. 建筑结构, 2021, 51(增刊2): 999-1007. https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG2021S2171.htm SHI Xinyu, YAO Yan, WANG Ling, et al. The influence of CDP model parameters based on the numerical simulation of uniaxial loading test[J]. Building Structure, 2021, 51(S2): 999-1007. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG2021S2171.htm
[12]	程展林, 左永振, 姜景山, 等. 粗粒料试验中界面摩阻力的试验研究[J]. 岩土工程学报, 2009, 31(3): 331-334. doi: 10.3321/j.issn:1000-4548.2009.03.004 CHENG Zhanlin, ZUO Yongzhen, JIANG Jingshan, et al. Experimental research on interface friction in granular materials[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(3): 331-334. (in Chinese) doi: 10.3321/j.issn:1000-4548.2009.03.004
[13]	刘学增, 王煦霖, 林亮伦. 60°倾角正断层黏滑错动对山岭隧道影响的试验研究[J]. 土木工程学报, 2014, 47(2): 121-128. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201402016.htm LIU Xuezeng, WANG Xulin, LIN Lianglun. Model experimental study on influence of normal fault with 60° dip angle stick-slip dislocation on mountain tunnel[J]. China Civil Engineering Journal, 2014, 47(2): 121-128. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201402016.htm
[14]	马亚丽娜, 崔臻, 盛谦, 等. 正断层错动对围岩-衬砌体系响应影响的离散-连续耦合模拟研究[J]. 岩土工程学报, 2020, 42(11): 2088-2097. doi: 10.11779/CJGE202011014 MA Yalina, CUI Zhen, SHENG Qian, et al. Influences of normal fault dislocation on response of surrounding rock and lining system based on discrete-continuous coupling simulation[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(11): 2088-2097. (in Chinese) doi: 10.11779/CJGE202011014
[15]	周光新, 崔臻, 盛谦, 等. 活动断裂错动位移模式对隧洞变形与内力的影响研究[J]. 防灾减灾工程学报, 2021, 41(6): 1323-1330, 1349. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK202106022.htm ZHOU Guangxin, CUI Zhen, SHENG Qian, et al. Study on the deformation and internal force of the tunnel under the displacement pattern of the active fault zone[J]. Journal of Disaster Prevention and Mitigation Engineering, 2021, 41(6): 1323-1330, 1349. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK202106022.htm