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

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

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]利用弹性力学Mindlin解，考虑顶进的正面附加推力、管节与土体间的摩擦力、土体损失等施工影响因素，提出适用于施工阶段的地面变形理论解答公式。李辉等^[2]针对矩形顶管施工，提出考虑摩擦力、正面顶推力、土体损失三者共同作用的地层位移计算方法，变形量叠加得到的地层变形预测模型可用于类似工程。焦义等^[3]为研究圆形顶管施工穿越特殊地层时对地表变形的影响，根据Mindlin解考虑顶进正面附加压力和侧向摩擦力，得到地表变形计算公式。张志伟等^[4]运用Mindlin解和随机介质理论，推导得到矩形顶管顶进期间地表隆沉变形的解析解，工程实例的运用证明所提方法的变形预测可行性。当前对于顶进施工引起地表变形的理论计算公式所考虑的施工因素较少，需要提出一种适用于大断面管幕箱涵顶进施工，可综合考虑多种施工因素和土体三维变形的计算公式，便于实现对地表变形的及时控制，为变形预测提供理论依据。

依托合肥市某大断面管幕箱涵顶进工程，通过理论计算、数值模拟和实测分析等方法，研究顶进施工过程中地表变形规律，提出管幕箱涵顶进施工引起的地表变形的理论计算公式，通过工程实例验证理论公式的合理性和准确性。

1. 土体受力模型与计算假定

Mindlin理论解是在布辛内斯克（Boussinnesq·J）解答的基础上推导出竖向或水平向荷载作用在半无限弹性体内部时体内任意一点的应力应变的数学解答，被广泛用于顶进施工对土体变形的影响研究^[5-6]。其将土体看作是各向同性半无限弹性体，利用弹性力学的Mindlin解^[7]，综合考虑各施工影响因素和具体受力特性，推导出管幕箱涵顶进施工过程引起的土体变形的三维理论解。

管幕箱涵顶进施工过程中开挖面上所受正面推进力作用如图 1所示，将土体视为半无限空间弹性体，以工作面竖向轴线与地表交点为坐标原点，管幕箱涵顶进方向为x轴，将x轴顺时针旋转90°得到y轴，从地表处开始向下扩展得到z轴。管幕箱涵顶进开挖面中心位置至地面原点的垂直距离为c，据此可利用Mindlin解答计算地层中任意一点M（x，y，z）的竖向位移。

图 1 水平集中力作用示意

Figure 1. Schematic diagram of action of horizontally concentrated force

下载: 全尺寸图片幻灯片

根据Mindlin解答，在弹性半无限空间内地表以下c深度处作用一个水平集中力时，土体中任一点M（x，y，z）在x，y，z轴上的位移分量如下：

$u = \frac{p}{{16{\text{π }}G(1 - \mu )}}\left[ {\frac{{3 - 4\mu }}{M} + \frac{1}{N} + \frac{{{x^2}}}{{{M^3}}} + \frac{{{\text{(}}3 - 4\mu {\text{)}}}}{{{N^3}}}} \right. + \\ ~~~~\frac{{2cz}}{{{N^3}}}\left( {1 - \frac{{3{x^2}}}{{{N^2}}}} \right) + \left. {\frac{{4(1 - \mu )(1 - 2\mu )}}{{N + z + c}}\left( {1 - \frac{{{x^2}}}{{N(N + z + c)}}} \right)} \right] \text{，}$

(1)

$v = \frac{{pxy}}{{16{\text{π }}G(1 - \mu )}}\left[ {\frac{1}{{{M^3}}} + \frac{{(3 - 4\mu )}}{{{N^3}}} - \frac{{6cz}}{{{N^5}}} \cdot } \right. \\ ~~~~~~~~\left. {\left( {1 - \frac{{3{x^2}}}{{{N^2}}}} \right) - \frac{{4(1 - \mu )(1 - 2\mu )}}{{N{{(N + z + c)}^2}}}} \right] \text{，}$

(2)

$w = \frac{{px}}{{16{\text{π }}G(1 - \mu )}}\left[ {\frac{{z - c}}{{{M^3}}} + \frac{{(3 - 4\mu )(z - c)}}{{{N^3}}} - } \right. \\~~~~~~~~~~~~~~~ \left. {\frac{{6cz(z + c)}}{{{N^5}}} + \frac{{4(1 - \mu )(1 - 2\mu )}}{{N(N + z + c)}}} \right] 。$

(3)

式中：x为到集中力作用点的水平距离，以水平力作用方向为正；y为到集中力作用点的横向水平距离；z为到地面的竖直距离，以竖直向下为正；c为管幕箱涵开挖面中心位置到地面原点的竖直距离；M为离水平力作用点的空间距离；N为离水平力作用点关于地表坐标的对称点的空间距离；p为水平集中力；G为土体剪切模量， $G = (1 - 2\mu {K_0}){E_{{\text{s0}}}}/\left[ {2(1 + \mu )} \right]$ ，其中，E_s0为土的压缩模量，K₀为静止土压力系数，μ为土体泊松比，三者根据工程详细勘察报告取值，下同。

用弹性理论分析土体时，可以用等效均质地基代替多层地基，将多层地基各土层的弹性模量加权平均值作为等效均质地基的弹性模量，该做法产生的误差并不大^[8]。基于该结论，可以将Mindlin解答用于成层土。

针对管幕箱涵顶进过程中由正面顶推力、管节与土体之间摩阻力、土体损失3种因素引起的土体沉降变形，根据Mindlin解答，推导得到考虑三维空间效应的大断面管幕箱涵顶进施工引起的土体沉降计算公式。将土体视为半无限空间弹性体，管幕箱涵顶进作用下土体受力模型简图如图 2所示。

图 2 顶进作用下土体受力模型

Figure 2. Schematic diagram of model for forces on soil under jacking action

下载: 全尺寸图片幻灯片

为简化后续计算，做如下假定：①土体为各向同性线弹性半无限空间体；②土体为不排水固结，土体变形仅在施工期间发生；③管幕箱涵顶进正前方的开挖面为荷载作用面，将正面推进力近似等效为矩形均布荷载；④顶进施工过程中管幕箱涵始终保持直线顶进，不存在纠偏等操作；⑤箱涵管节外表面与周围土体之间的摩擦力为均匀分布荷载；⑥管幕箱涵顶进施工仅考虑空间位置上的改变，不考虑时间效应。

2. 地表变形解析解答

2.1 正面推进力引起的土体变形

管幕箱涵顶进施工由正面推进力引起的土体变形为管幕顶进施工与箱涵顶进施工两部分土体变形量之和。研究圆形断面钢管幕正面推进力引起的土体变形时假设顶进正前方的开挖面上的顶推力为均布荷载，在荷载作用的圆截面内取微分面积 $r{\text{d}}r{\text{d}}\theta$ ，根据Mindlin解答，由圆钢管幕均布正面推进力引起的土体中任意一点的竖向位移为

$\begin{aligned} & w_1=\frac{p_1 x}{16 \pi G(1-\mu)} \\ & \int_0^{2 \pi} \int_0^{R / 2}\left[\frac{z-h+r \sin \theta}{M_1^3}+\frac{(3-4 \mu)(z-h+r \sin \theta)}{N_1^3}-\right. \\ & \left.\frac{6 z(h-r \sin \theta)(z+h-r \sin \theta)}{N_1^5}+\frac{4(1-\mu)(1-2 \mu)}{N_1\left(N_1+z+h-r \sin \theta\right)}\right] r \mathrm{d} r \mathrm{d} \theta_0 \end{aligned}$

(4)

式中：x为到开挖面的水平距离，以顶进方向为正；y为与开挖面竖向中轴线之间的横向水平距离；z为到地面的竖直距离，以向下为正；h为开挖面中心点与地面之间的垂直距离；R为钢管半径；p₁为正面推进力；G为土体剪切模量； ${M_1} =$ $[{x^2} + {(y + r\cos \theta )^2} + (z - h +$ $r\mathrm{sin}\theta {)}^{2}{]}^{\frac{1}{2}}；{N}_{1}={[{x}^{2}+{(y+r\mathrm{cos}\theta )}^{2}+{(z+h-r\mathrm{sin}\theta )}^{2}]}^{\frac{1}{2}}。$

在管幕箱涵法顶进过程中，箱涵在管幕群的保护下顶进，降低了对周围土体的扰动，减小了地表变形，所以研究箱涵顶进过程对地表变形的影响应考虑到管幕对其的保护作用。基于当层法原理，将既有钢管幕结构等效为当层厚度的与周围土体相同的地层^[9]。

当层法原理示意图如图 3所示，上下地层厚度分别为h₁，h₂，弹性模量分别为E₁，E₂，依据当层法将两地层等效为具有同样力学参数的当层土，等效以后的相关力学参数为E₁，将下层土等效厚度用当层厚度 ${h'_2}$ 表示，与 ${h_2}$ 的关系满足

${h'_2} = {h_2}{\left( {\frac{{{E_2}}}{{{E_1}}}} \right)^a} 。$

(5)

图 3 当层法原理示意

Figure 3. Schematic diagram of principle of layered method

下载: 全尺寸图片幻灯片

式中：a为当层指数，按照经验取值为0.33^[10]。

得到既有上排钢管幕当层法等效后的当层土厚度为13.63 m，如图 4所示。在计算箱涵顶进对地表变形的影响时，计算所用土层分布比实际土层分布多一层h₄=13.63 m的土层，即钢管幕当层转换后的土层，计算中箱涵中心到地表距离h’比实际距离h多13.63 m，这样可以考虑管幕对箱涵顶进的影响。

图 4 管幕当层法置换示意

Figure 4. Schematic diagram of replacement of curtain by layered method

下载: 全尺寸图片幻灯片

研究矩形断面箱涵正面推进力引起的土体变形时假设顶进正前方的开挖面上的顶推力为均布荷载，利用微积分思想，在开挖面上选取一个微分单元（ ${\text{d}}y'{\text{d}}z'$ ），微分单元的坐标表示为（ $y'$ ， $z'$ ），其上作用的微集中力为 $p{\text{d}}y'{\text{d}}z'$ ，最后利用Mindlin解答对整个开挖面进行积分，即得到正面推进力引起的土体中任意点的竖向位移。根据Mindlin解答，由箱涵均布正面推进力引起的土体中任意一点的竖向位移为

$\begin{gathered} u_{z 1}=\int_{-A / 2}^{A / 2} \int_{-B / 2}^{B / 2} \frac{P x}{16 \pi G(1-\mu)}\left[\frac{z-h^{\prime}-z^{\prime}}{R_1^3}+\frac{(3-4 \mu)\left(z-h^{\prime}-z^{\prime}\right)}{R_2^3}-\right. \\ \left.\frac{6 z\left(h^{\prime}+z^{\prime}\right)\left(z+h^{\prime}+z^{\prime}\right)}{R_2^5}+\frac{4(1-\mu)(1-2 \mu)}{R_2\left(R_2+z+z^{\prime}+h^{\prime}\right)}\right] \mathrm{d} y^{\prime} \mathrm{d} z^{\prime} 。 \end{gathered}$

(6)

式中：A为箱涵矩形断面的长度；B为箱涵矩形断面的宽度； $x'$ 为当到开挖面的水平距离，以顶进方向为正； $y'$ 为与开挖面竖向中轴线之间的横向水平距离； $z'$ 为到地面的竖直距离，以向下为正； $h'$ 为当层法转换后开挖面中心到地表的垂直距离；R₁为离水平力作用点的空间距离；R₂为离水平力作用点关于地表坐标的对称点的空间距离；P为水平集中力；G为土体剪切模量。

2.2 摩擦力引起的土体变形

管幕箱涵顶进过程中由摩擦力引起的土体变形为管幕顶进与箱涵顶进两个过程中摩擦力分别引起的变形量之和，研究圆形断面管幕管节外表面摩擦力引起的土体变形时假设顶进过程中管幕箱涵管节外表面与周围土体之间的摩擦力为均布荷载。在荷载作用的管节表面取微分面积 $R{\text{d}}l{\text{d}}\theta$ ，根据Mindlin解答，由钢管幕摩擦力引起的土体中任意一点（x，y，z）的竖向位移为

$\begin{aligned} w_2= & \frac{p_2 R}{16 \pi G(1-\mu)} \int_0^{2 \pi} \int_0^L(x+l) . \\ & {\left[\frac{z-h+R \sin \theta}{M_2^3}+\frac{(3-4 \mu)(z-h+R \sin \theta)}{N_2^3}-\right.} \\ & \frac{6 z(h-R \sin \theta)(z+h-R \sin \theta)}{N_2^5}+ \\ & \left.\frac{4(1-\mu)(1-2 \mu)}{N_2\left(N_2+z+h-R \sin \theta\right)}\right] \mathrm{d} / \mathrm{d} \theta 。 \end{aligned}$

(7)

式中：x为到开挖面的水平距离，以顶进方向为正；y为与开挖面竖向中轴线之间的横向水平距离；z为到地面的竖直距离，以向下为正；L为管节长度；R为钢管半径；p₂为管节与土体之间的单位面积摩擦力；h为开挖面中心到地表的垂直距离；G为土体剪切模量。 ${M_2} = \sqrt {{{(x + l)}^2} + {{(y + R\cos \theta )}^2} + {{(z - h + R\sin \theta )}^2}}$ ， ${N_2} =$ $\sqrt {{{(x + l)}^2} + {{(y + R\cos \theta )}^2} + {{(z + h - R\sin \theta )}^2}}$ 。

研究矩形断面箱涵管节摩擦力引起的土体变形时假设顶进过程中箱涵管节外表面与周围土体之间的摩擦力为均布荷载，均布荷载的大小与其所在位置有关。对于矩形断面箱涵，在其顶进过程中4个外表面均与周围土体接触，其中左右两侧接触面上的摩擦力对称相等，上下两个接触面上的摩擦力因应力状态不同需要分别计算。

左右接触面上摩擦力所引起的土体竖向位移为

$\begin{aligned} & u_{z 2}=\int_{-B / 2}^{B / 2} \int_0^L \frac{P_{\mathrm{f}}(x+L)}{16 \pi G(1-\mu)}\left[\frac{z-h^{\prime}-z^{\prime}}{R_1^3}+\frac{(3-4 \mu)\left(z-h^{\prime}-z^{\prime}\right)}{R_2^3}-\right. \\ & \left.~~\frac{6 z\left(h^{\prime}+z^{\prime}\right)\left(z+h+z^{\prime}\right)}{R_2^5}+\frac{4(1-\mu)(1-2 \mu)}{R_2\left(R_2+z+z^{\prime}+h^{\prime}\right)}\right] \mathrm{d} z^{\prime} \mathrm{d} l 。 \end{aligned}$

(8)

上接触面的摩擦力所引起的土体竖向位移为

$\begin{aligned} & u_{z 3}=\int_{-A / 2}^{A / 2} \int_0^L \frac{P_{\mathrm{t1}}(x+L)}{16 \pi G(1-\mu)}\left[\frac{z-h^{\prime}+0.5 B}{R_1^3}+\right. \\ & \frac{(3-4 \mu)\left(z-h^{\prime}+0.5 B\right)}{R_2^3}-\frac{6 z\left(h^{\prime}-0.5 B\right)\left(z+h^{\prime}-0.5 B\right)}{R_2^5}+ \\ & \left.~~~~~~~~~~~~~~~~~~~~\frac{4(1-\mu)(1-2 \mu)}{R_2\left(R_2+z+h^{\prime}-0.5 B\right)}\right] \mathrm{d} y^{\prime} \mathrm{d} l \quad 。 \end{aligned}$

(9)

下接触面的摩擦力所引起的土体竖向位移为

${u_{z4}} = \int_{{{ - A} \mathord{\left/ {\vphantom {{ - A} 2}} \right. } 2}}^{{A \mathord{\left/ {\vphantom {A 2}} \right. } 2}} {\int_0^L {\frac{{{P_{t2}}(x + L)}}{{16\pi G(1 - \mu )}}\left[ {\frac{{z - h' - 0.5B}}{{R_1^3}} + } \right.} } \\ \frac{{(3 - 4\mu )(z - h' - 0.5B)}}{{R_2^3}} - \frac{{6z(h' + 0.5B)(z + h' + 0.5B)}}{{R_2^5}} + \\ \left. ~~~~~~~~~~~~{\frac{{4(1 - \mu )(1 - 2\mu )}}{{{R_2}({R_2} + z + h' + 0.5B)}}} \right]dy'dl 。$

(10)

式中：A为箱涵矩形断面的长度；B为箱涵矩形断面的宽度；L为箱涵管节长度； $H'$ 为当层法转换后箱涵开挖面中心到地表的垂直距离；R₁为离水平力作用点的空间距离；R₂为离水平力作用点关于地表坐标的对称点的空间距离； ${P_f}$ 为箱涵左右接触面的摩擦力； ${P_{t1}}$ 为箱涵上接触面的摩擦力； ${P_{t2}}$ 为箱涵下接触面的摩擦力；G为土体剪切模量。

2.3 土体损失引起的土体变形

现有研究中提出的经验公式多依托单线盾构工程，而本文依托于管幕工程，需要考虑群管顶进对地表变形的影响。在平行顶进钢管的过程中，如果两个钢管的轴线距离较小，先施工钢管对周围土体的扰动会使得后施工顶管对土体的扰动变大，导致沉降槽宽度增大，管幕顶进施工过程中多根钢管顶进加剧了对地表土体的影响，沉降槽宽度系数i增大，应当对理论公式进行修正。

Peck沉降槽宽度系数计算公式为 $i = R{(h/2R)^n}$ ， n=0.8~1.0，土质越软则n值越大。Loganathan^[12]给出的计算式为 $i = 0.6(h + R)$ ，即使是软土（n=1.0），Peck计算式所得结果仍小于Loganathan计算式。魏纲^[11]认为黏性土地区i值与 $R + h\tan ({45^ \circ } - {\varphi \mathord{\left/ {\vphantom {\varphi 2}} \right. } 2})$ 呈线性关系，故将管幕顶进沉降槽宽度系数i计算公式修正为

$i=0.6\left[R+h \tan \left(45^{\circ}-\varphi / 2\right)\right] 。$

(11)

Loganathan^[12]于1998年拟合得到考虑钢管周围土体的非等量径向移动的土体损失为

${V_l} = {V_{loss}}\exp \left( { - \frac{{1.38{y^2}}}{{{{(h + R)}^2}}}} \right) 。$

(12)

钢管单位长度土体损失量为

${V_{loss}} = \pi {R^2} - \pi {(R - g/2)^2} \text{，}$

(13)

$g = 2R(1 - \sqrt {1 - \eta } ) 。$

(14)

考虑土体泊松比时的钢管轴线上方最大沉降量计算公式：

${S_{1\max }} = 4(1 - \mu )\frac{{{V_{loss}}}}{{\pi h}} 。$

(15)

本文提出管幕顶进施工由土体损失引起的横向地表沉降修正Sagaseta公式为

${S_1} = \frac{{2(1 - \mu ){V_l}}}{\pi }\frac{h}{{{y^2} + 4 \times 0.36{{\left[ {R + h\tan ({{45}^ \circ } - {\varphi \mathord{\left/ {\vphantom {\varphi 2}} \right. } 2})} \right]}^2}}} \cdot \\\left\{1-\frac{x}{\sqrt{{x}^{2}+{y}^{2}+4\times 0.36{\left[R+h\mathrm{tan}(45°-\varphi /2)\right]}^{2}}}\right\}\mathrm{exp}\left(-\frac{1.38{y}^{2}}{{(h+R)}^{2}}\right)。$

(16)

式中：R为掘进机外径；g为等效土体损失参数；η为土体损失百分率；x为到开挖面的水平距离，以顶进方向为正；y为与开挖面竖向中轴线之间的横向水平距离； ${V_l}$ 为钢管单位长度非等量径向移动的土体损失量；S₁为距离钢管中轴线x处的地表沉降量。

式（16）既考虑钢管顶进过程中发生的非等量径向移动，即钢管并不会悬空处于开挖空间的中心位置，而是会落到开挖空间的底部边界，所引起的地表沉降大于钢管发生等量径向移动引起的地表沉降，又考虑到土体泊松比对地表沉降的影响，同时考虑到群管顶进会增大地表沉降槽宽度，对群管顶进的沉降槽宽度系数进行修正，保证计算公式更符合实际施工情况。

在管幕箱涵法顶进过程中，箱涵在钢管幕群的保护下进行顶进，所引起的土体损失小于单线矩形顶管施工所引起的土体损失，即地层损失率更小，地表变形沉降更小。由于地层损失率η主要与工程地质情况、水文地质情况、隧道施工方法、施工技术水平以及工程管理经验等因素有关，因此参数的取值离散性较大且依赖于地区经验。计算中应选取较小的地层损失率η值，合肥地区黏土地层损失率通常为0.5%~2.5%，本工程中箱涵顶进时选为1.0%。

计算箱涵顶进土体损失对地表变形的影响时，应使用当层法转换后的土层分布，将钢管幕等效为与箱涵穿越土层具有同样力学参数的当层土，计算中箱涵中心到地表距离比实际距离多13.63 m，以考虑钢管幕的保护作用。本文提出矩形箱涵顶进施工由土体损失引起的修正Sagaseta横向地表沉降公式为

${S}_{2}=\frac{2(1-\mu )\eta AB{h}^{\prime }}{\pi ({y}^{2}+h{\text{'}}^{2})}\left[1-\frac{x}{\sqrt{{x}^{2}+{y}^{2}+h{\text{'}}^{2}}}\right]\mathrm{exp}\left(-\frac{1.38{y}^{2}}{{\left({h}^{\prime }+\frac{B}{2}\right)}^{2}}\right)。$

(17)

式中：S₂为距离箱涵中轴线y处的地表沉降量；x为到开挖面的水平距离，以顶进方向为正；y为到箱涵竖向中轴线的水平距离；i为地表沉降槽宽度系数；V_s为箱涵单位长度土体损失量；η为地层损失率；A为箱涵矩形断面的长度；B为箱涵矩形断面的宽度； $h'$ 为当层法转换后箱涵中心到地表的垂直距离。

2.4 地表变形理论公式

管幕箱涵顶进施工引起地表变形的因素主要是正面推进力、管节四周摩擦力和土体损失三个因素，现有研究认为各因素对变形的影响是独立的，变形可以进行叠加^[13]，地表的竖向总变形量是管幕与箱涵所引起的变形量之和，总计算公式为

${\omega _{pipe}} = {\omega _1} + {\omega _2} + {S_1} \text{，}$

(18)

${u_{box}} = {u_{z1}} + 2{u_{z2}} + {u_{z3}} + {u_{z4}} + {S_2} \text{，}$

(19)

$u = {\omega _{pipe}} + {u_{box}} 。$

(20)

式中： ${\omega _{pipe}}$ 为钢管顶进引起的地表变形；ω₁为钢管顶推力引起的地表变形；ω₂为钢管摩擦力引起的地表变形；S₁为钢管土体损失引起的地表变形；u_box为箱涵顶进引起的地表变形；u_z1为箱涵顶推力引起的地表变形；u_z2为箱涵左右两侧接触面上摩擦力引起的地表变形；u_z3为箱涵上接触面摩擦力引起的地表变形；u_z4为箱涵下接触面摩擦力引起的地表变形；S₂为箱涵土体损失引起的地表变形；u为管幕箱涵顶进结束引起的地表变形。

该公式结合了由于土体损失引起的土体变形，叠加正面推进力、管节上下侧和左右侧摩擦力和土体损失引起的土体变形，得到的三维变形量与实际情况产生的变形量较接近，可以合理反映管幕箱涵顶进动态施工过程引起的土体变形，同时公式中涉及的参数易确定，公式运用较为方便。

3. 算例分析

为验证第二节所提出的管幕箱涵顶进引起地表变形理论公式的合理性，依托合肥大断面管幕箱涵顶进施工工程实例，通过编写计算程序求解地表各点位移值。根据地表变形理论公式分别计算三个监测断面处管幕和箱涵的正面顶推力、摩擦力和土体损失对地表沉降的影响。

隧道断面为矩形，主体为单箱两孔矩形框架结构，结构全宽32.8 m，全高12.05 m，其中箱涵断面宽29.4m，高8.75 m。根据设计要求，箱涵周围共计施工44根Φ1.6 m×16 mm钢管，钢管长135 m，其中上、下部各16根，钢管间距184 mm，左、右边墙各6根，钢管间距209 mm。管幕与箱涵上部及左右两侧间隙10 cm，下部管幕与箱涵紧贴。

模型宽度135 m，高度60 m，长度135 m，共816123个节点，1036147个单元。为了保证计算效率，将箱涵管节设置为5 m一节进行顶进施工模拟，共计27段管节。土层参数见表 1，管幕箱涵结构参数见表 2，管幕箱涵监测断面及测点布置示意如图 5所示。

表 1 土层基本参数

Table 1. Basic parameters of soil layers

序号	名称	厚度/m	重度/(kN·m^-³)	K₀	泊松比	孔隙比	黏聚力/kPa	摩擦角/(°)	压缩模量/MPa
1	①杂素填土	2	19.0	0.65	0.35	0.907	10.0	8.0	2.50
2	⑤₁黏土	3	19.8	0.50	0.33	0.729	73.3	14.4	13.86
3	⑤₂黏土	16	20.0	0.44	0.31	0.696	79.2	15.0	15.93
4	⑥全风化土	13	21.3	0.37	0.28	0.692	29.7	15.9	11.15
5	⑦强风化泥质砂岩	26	23.0	0.25	0.23	0.628	12.0	48.0	45.00

下载: 导出CSV

| 显示表格

表 2 管幕箱涵结构参数

Table 2. Structural parameters of curtain box culvert

结构	材料	重度/ (kN·m^-³)	弹性模量/ GPa	泊松比
箱涵	C50混凝土	25.0	34.5	0.16
钢管	Q235钢	78.5	210.0	0.30

下载: 导出CSV

| 显示表格

图 5 管幕箱涵监测点及监测断面布置

Figure 5. Schematic diagram of layout of measuring points of curtain box culvert

下载: 全尺寸图片幻灯片

在监测断面1处，管幕顶进引起的地表变形理论值、数值模拟计算值和实测值对比如图 6所示。在顶推力、摩擦力、土体损失共同作用下，管幕顶进引起的最大沉降为27.92 mm，箱涵顶进引起的最大沉降为40.75 mm，地表最大沉降量发生在箱涵中轴线处，其理论值为68.67 mm，数值模拟计算值为72.12 mm，监测值为66.92 mm，3个最大沉降值较为接近，整体变形趋势较为一致，数值差距在合理范围内，因为现场施工较为复杂，解析解答和数值模拟很难全面准确反映实际情况。

图 6 断面一管幕箱涵顶进结束地表变形

Figure 6. Comparison of surface deformations at end of jacking of curtain box culvert in Section 1

下载: 全尺寸图片幻灯片

在监测断面2处，管幕顶进引起的地表变形理论值、数值模拟计算值和实测值对比如图 7所示。在顶推力、摩擦力、土体损失共同作用下，管幕顶进引起的最大沉降为25.12 mm，箱涵顶进引起的最大沉降为35.15 mm，地表最大沉降量发生在箱涵中轴线处，其理论值为60.27 mm，数值模拟计算值为66.97 mm，监测值为61.95 mm。理论值与数值模拟计算值整体变形趋势较为一致，理论值相比数值模拟计算值较小，因为数值模拟中考虑到了实际顶进施工中有8 m超挖。由于现场实际工程有跑道存在，一定程度上抑制了变形的发展，所以变形实测值较小。

图 7 断面二管幕箱涵顶进结束地表变形

Figure 7. Comparison of surface deformations at end of jacking of curtain box culvert in Section 2

下载: 全尺寸图片幻灯片

在监测断面3处，管幕顶进引起的地表变形理论值、数值模拟计算值和实测值对比如图 8所示。在顶推力、摩擦力、土体损失共同作用下，管幕顶进引起的最大沉降为19.86 mm，箱涵顶进引起的最大沉降为28.12 mm，地表最大沉降量发生在箱涵中轴线处，其理论值为47.98 mm，数值模拟计算值为57.59 mm，监测值为52.08 mm。理论值与数值模拟计算值变化趋势一致，理论值较小于数值模拟计算值，原因是有限元模拟时考虑了实际施工中的超挖，二者差距在合理范围内。

图 8 断面三管幕箱涵顶进结束地表变形

Figure 8. Comparison of surface deformations at end of jacking of curtain box culvert in Section 3

下载: 全尺寸图片幻灯片

总体上，断面一处的地表变形大于断面二和断面三，始发位置处地表变形最大，距开挖面越近地表沉降越小，与实际顶进引起的变形规律相符。3个断面处的地表变形理论值、数值模拟计算值和实测值数值较为接近，整体变化趋势一致，曲线相关性较高，验证了所提出的综合考虑顶进施工因素影响的管幕箱涵顶进地表变形理论解答公式的合理性，理论解答与数值模拟可较好地综合预测大断面管幕箱涵顶进对地表变形的影响并指导实际施工。

4. 结论

针对现有顶进施工引起地表变形的理论计算公式所考虑的施工因素较少的问题，提出了一种适用于大断面管幕箱涵顶进施工，可综合考虑多种施工因素和土体三维变形的计算公式，并通过工程实例验证了理论公式的准确性。

（1）根据理论计算、数值模拟和现场监测数据，在管幕箱涵施工过程中，地表最大沉降量发生在箱涵中线处，向两侧逐渐减小。

（2）运用Mindlin解答和修正Sagaseta公式计算得到管幕箱涵顶进施工的地表变形计算公式。将复合地层简化为简单的均匀地层，通过积分方法得到地层沉降计算结果，计算快速、应用简单，而由于对复杂地层的简化，计算所得的地层沉降与实际相比偏大。通过工程实例计算分析地表沉降变化规律，与数值模拟计算结果和监测值对比，验证了理论公式的合理性。

（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