深埋小净距多线平行盾构掘进相互作用分析

付 钊，柯宁静，卢康明，郭萧阳，张孟喜

（1.上海大学土木工程系，上海 200072；2.上海隧道工程有限公司盾构工程分公司，上海 200023）

硬X 射线自由电子激光装置位于上海张江综合性国家科学中心核心区域，装置总长约3.1 km，包含超导直线加速器隧道、波荡器隧道、光束线隧道等结构。该装置建成后，将成为世界上最高效和最先进的自由电子激光用户装置之一。硬X 射线自由电子激光具有更高的亮度、更短的脉冲结构和更好的相干性，能提供的X 射线峰值亮度比第三代同步辐射光源高109倍。作为世界顶级的科研基础设施，硬X 射线自由电子激光装置将刺激多类学科的发展，为物理、化学、生命科学、材料科学、能源科学等多学科提供高分辨成像、超快过程探索、先进结构解析等尖端研究手段。该装置的土建部分由多条平行盾构隧道组成。在深埋小净距的特殊情况下，多线盾构隧道的施工必定会给各条隧道，尤其是首挖隧道带来重大影响，也会对随后的科学实验装置的使用产生不利影响。因此，对于深埋小净距多线隧道盾构推进时，各条隧道，尤其是首挖隧道的断面变形历时曲线和收敛规律有必要进行深入研究。

目前多线平行隧道的研究主要集中于双线。针对双线隧道开挖引起的地表沉降、隧道变形等问题，国内外学者通过理论研究[1-6]、数值模拟[7-14]、模型试验[15-16]和现场监测[17-19]等不同方面进行了深入的探讨。对于三线隧道，王明年等[20]通过模型试验对软弱围岩下三孔小间距浅埋暗挖隧道的地表沉降控制技术进行了研究；陈越峰等[21]针对上海某地铁近距离三线并行盾构区间隧道的施工进行现场监测分析研究，总结归纳了盾构施工对周围深层土体水平位移、深层土体沉降、孔隙水压力的基本影响规律；赵建华[22]采用数值模拟分析方法研究了三孔隧洞不同施工开挖顺序工况下的地表沉降、初期支护位移、围岩塑性区和初期支护受力情况，并进行对比分析。

由于目前对多线平行盾构的研究主要集中在双线，而对深埋三线及以上平行盾构的研究较少。本文通过结合上海硬X 射线激光装置土建工程，探究了深埋小净距多线平行盾构隧道在不同开挖顺序条件下，各盾构隧道，尤其是首挖隧道的变形情况，对其断面的变形收敛规律进行了研究，提出了合适的隧道开挖次序；探究了不同施工参数对首挖隧道管片变形的影响，为后续施工的顺利推进提供了保障。

1 工程概况

上海硬X 射线自由电子激光装置土建工程位于上海市浦东新区张江园区祖冲之路与华夏中路之间，西侧为罗山路高架和磁悬浮。项目主要由长约3.2 km的地下隧道、5 个竖井及竖井附近的地面设施组成。地下隧道段整体呈南北走向。拟建工程地下隧道分为加速器段、波荡器段、光束线段和超长光束线段。加速器段由1 条隧道构成，长度约1 285 m；波荡器段由3 条隧道构成，长度各约400 m；光束线段由3 条隧道构成，长度各约300 m；超长光束线由3 条隧道构成，长度各约785 m。其对应的平面示意图如图1所示，盾构之间的关系如图2所示。隧道均采用拼装式混凝土衬砌管片，外径为6.6 m，内径为5.9 m，管片厚度为350 mm，环宽1.5 m。管片混凝土强度等级为C55，抗渗等级P12，错缝拼装。隧道间距为4 m，盾构机覆土深度30 m，C/D约为4.55。工程采用3 台土压式平衡盾构机施工，相隔一段时间渐次出发。本工程自2019年开始施工，现阶段已完成超深工作井工程，盾构段还未开始施工。施工方下阶段将会按照本文所提供的开挖次序进行工程测试。盾构推进的准备工作正在有条不紊的推进当中。

图1 硬X 射线平行盾构工程平面图Fig.1 Layout of the hard X-ray parallel shield engineering

本区间地层由黏性土、粉土、砂性土交互沉积而成，按土层结构、成因及性状特性将其划分为7 层，土层分布情况如图2所示，由地勘报告可得各土层物理力学参数见表1。地下水为潜水，埋深2 m。

图2 小净距平行盾构隧道横剖面图（单位：m）Fig.2 Cross section of a parallel shield tunnel（unit：m）

表1 土体物理力学参数Table 1 Mechanical parameters of soils

2 多线盾构施工模拟

2.1 多线盾构施工模拟情况

为了更准确地模拟多线平行盾构的开挖情况，根据对称原则以及实际的盾构开挖情况（以3#-4#、4#-5#区间为例，图3）采用2 种不同的开挖方式，开挖示意图如图4所示。其中5#、6#、7#隧道分别对应左、中、右线隧道。

图3 实际盾构掘进情况Fig.3 Actual conditions of shield tunneling

图4 开挖方式示意图Fig.4 Schematic diagram of the excavation method

（1）左线-右线-中线

考虑到实际多条隧道开挖间隔以及模型的尺寸大小，以

Δ=120 m 为距离进行盾构隧道开挖。整个多线隧道开挖分为4 步：①首先开挖左线隧道；②左线隧道开挖 Δ距离后，右线开挖；开挖60 m 后，左线隧道开挖完成；③右线隧道开挖 Δ距离之后，中线开挖；开挖60 m 后，右线隧道开挖完成；④中线隧道开挖完成。

（2）中线-右线-左线

与工况1 类似，以 Δ=120 m 为距离进行盾构隧道开挖。整个多线隧道开挖分为4 步：①首先开挖中线隧道；②中线隧道开挖 Δ距离后，右线开挖；开挖60 m 后，中线隧道开挖完成；③右线隧道开挖 Δ距离之后，左线开挖；开挖60 m 后，右线隧道开挖完成；④左线隧道开挖完成。

2.2 三维有限元模型建立

根据隧道间的穿越位置关系，考虑盾构掘进对土体的扰动及其影响范围，利用有限元软件ABAQUS 构建三维有限元模型，模拟实际工况。模型的尺寸以及盾构相对位置示意图见图5。模型边界条件：所有侧面限制其法向位移；底面设置为3 个方向全约束；上表面取至地表，为自由面，保证模型的竖向变形不受边界条件约束。盾构隧道的开挖采用生死单元法模拟。

将衬砌管片壁后受扰动土体和注浆情况概化为一均质、等厚的等代层进行计算，根据文献[13]的方法，将等代层的厚度取为0.14 m。

2.3 模型参数选择

土体的本构模型采用Mohr-Coulomb 模型，各土层物理力学参数见表1。衬砌管片采用弹性模型。C55混凝土根据规范取弹性模量为35.5 GPa，泊松比0.17。为考虑环与环之间连接对隧道整体刚度的降低，将模型中衬砌弹性模量取为C55 的85%，即折减后衬砌弹性模量E= 30.2 GPa。模拟盾构机施工过程时，开挖面支护压力采用0.73 倍的土层静止土压力，并考虑应力梯度，将开挖面的支护压力设置为梯形，上小下大；为了简化施工步骤，将6 环设置为1 个开挖步，共计46 个开挖步。开挖后立即施加支护力，并添加等代层和衬砌管片。

不施加地下水时，土体与混凝土管片单元均采用实体单元C3D8，模型共计75 840 个实体单元。施加地下水时，采用设置孔隙水压力边界的方法进行施做。在土层2 m 深处设置孔隙水压力为0，然后以梯度的方式向下进行设置，梯度值为水的重度。分析步采用Soils 稳态。土体采用C3D8P 单元，管片和等代层依然采用C3D8 单元。

3 盾构施工影响数值模拟结果分析

3.1 开挖过程中首挖隧道管片变形分析

为简化隧道断面管片变形的描述，取首挖隧道第10 节管片，即图5中LC、CC、RC 位置处管片为研究对象。选取内径顶部处点A、右端点C、底部处点D 和左端点B 的变形进行研究（图6）。

图5 盾构开挖有限元模型Fig.5 Finite element model of shield excavation

图6 管片截面计算测点位置Fig.6 Situation of the segment section

根据各点竖向和水平变形增量数据与隧道开挖历程的关系，可以得到不同施工阶段首挖隧道管片的变形增量曲线。

为了更合理地描述管片各特征点的位移增量随时间历程的变化，选取模型的历程分析步为横坐标，用以描述整个模型开挖的过程。由于本节探究的是后挖隧道对首挖隧道管片的变形影响，因此时间历程由第二条隧道开挖起始，直至所有隧道开挖完成。

式中：λ—支护应力比，本次λ取为0.5，0.6，0.7，0.8；

σs—开挖面中心顶推力；

σ0—隧道中心原始地层静止土压力。

注浆的弹性模量受材料种类、硬化时间以及周围土体的影响较大，对管片的变形收敛具有一定的影响，因此也对其进行了参数的敏感性分析，这里取弹性模量4，6，8，10 MPa 进行分析。由表2可以看出有无地下水对管片收敛的影响较小，因此选取无地下水时的情况进行参数分析。

3.1.1 工况1“左-右-中”管片变形分析

选取LC 管片顶部处与底部处的竖直U3 变形增量和左端点与右端点的水平U1 变形增量绘制变形增量图，进行管片变形分析。

初到广州的五建迅速瞄准和服务华南市场。但与此同时，五建面临着依赖传统工程建设项目、效益增长方式单一、市场投资领域变化、生产方式转变、资源环境和生产要素制约以及劳动力成本上升等问题。眼见传统工程建设的路越来越窄，未来如何生存已经成为不得不思考的问题。“我们清醒地认识到，只有加快企业转型升级步伐，开发新的效益增长领域，企业才能焕发出新的生机和活力。”五建党委书记蒋德军说。

（1）地下水的影响

由图7可知，随着后2 条隧道的开挖，首挖隧道管片顶部和底部的竖向位移增量呈现一小一大“同向双峰”的变化形式，且双峰出现的位置皆位于另外2 条隧道开挖面临近所选取管片的时间。峰值不同的原因是由于后2 条隧道与首挖隧道的不同间距所导致的。由图7（a）（b）可知，顶部和底部点的位移增量皆为负值，即首挖隧道在另外2 条隧道开挖的影响下，隧道管片整体不断下移，这是由于后2 条隧道开挖后应力释放的结果。在地下水的影响下，由于水压的作用，导致隧道的竖向变形增量变大，A 点竖向变形最大值由无水时的-0.066 mm 变为-0.121 mm；D 点的竖向变形最大值由无水时的-0.026 mm 变为-0.072 mm。由图7（c）（d）可知，在地下水的影响下，管片水平变形增量略有增大，尤其是右侧隧道开挖导致的LC 管片变形增量显著增大。B 点水平变形增量由0.111 mm变为0.219 mm，C 点的水平变形增量由0.116 mm 变为0.23 mm。综上可知，地下水对工程有着一定的影响。

图7 地下水对管片变形的影响（工况1）（E=10 MPa，λ=0.7）Fig.7 Influence of groundwater on the deformation of pipe segment（case 1）（E=10 MPa,λ=0.7）

（2）注浆弹性模量的影响

由图8（a）（b）可知，在4 种工况下，随着注浆弹性模量的增大，隧道管片顶部和底部的竖向变形增量随之减小。由图8（c）（d）可知，在4 种工况下，随着注浆弹性模量的增大，右侧隧道开挖对首挖隧道的水平变形影响随之减小，而中部隧道开挖对首挖隧道水平变形的影响先增大后略微减小。

图8 注浆弹性模量对管片变形的影响（工况1）（λ=0.7）Fig.8 Influence of elastic modulus of grouting on the deformation of pipe segment（case 1）（λ=0.7）

（3）顶推力的影响

如图9所示，在4 种工况下，随着开挖面支护力的增加，A 点和D 点的竖向位移增量不断减小，B 点和C 点的水平位移增量不断增大。但支护力在0.6P（P为隧道中心原始地层静止土压力）和0.7P时隧道管片变形增量基本一致，因此将支护力维持在0.6P～0.7P之间时，隧道管片变形比较稳定，更容易控制。

图9 顶推力对管片变形的影响（工况1）（E=10 MPa）Fig.9 Influence of supporting force of excavation face on the deformation of pipe segment（case 1）（E=10 MPa）

3.1.2 工况2“中-右-左”管片变形分析

选取CC 管片顶部处与底部处的竖直U3 变形增量和左端点与右端点的水平U1 变形增量绘制变形增量图如图10所示。

（1）地下水的影响

由图10（a）（b）可知，在没有地下水的影响下，由于土体开挖导致的应力释放使得首挖的中部隧道整体上浮，双峰的增量都为正值。而在有地下水的影响下，使得首挖的中部隧道整体下沉，双峰的增量均为负值。由图10（c）（d）可知，在有地下水的影响下，管片左右两侧的水平变形增量相较于无地下水时显著增大，左端点最大值由0.295 mm 变为0.427 mm，右端点最大值由0.325 mm 变为0.473 mm。地下水对衬砌的变形有着一定的影响。

图10 地下水对管片变形的影响（工况2）（E=10 MPa,λ=0.7）Fig.10 Influence of groundwater on the deformation of pipe segment（case 2）（E=10 MPa,λ=0.7）

（2）注浆弹性模量的影响

如图11（a）（b）所示，在4 种工况下，随着注浆弹性模量的增大，隧道管片顶点和底部点的竖向变形增量绝对值不断减小，直至10 MPa 时变为正值。这是由于注浆弹性模量过低，会使得衬砌管片直接承受来自周围的围岩压力作用，导致衬砌发生过量变形。而在图11（c）（d）中，管片左右两侧的水平变形会由于等代层的刚度变化而产生先增大后减小的变化趋势。

图11 注浆层弹性模量对管片变形的影响（工况2）（λ=0.7）Fig.11 Influence of elastic modulus of grouting layer on the deformation of pipe segment（case 2）（λ=0.7）

（3）顶推力的影响

由图12（a）（b）可知，在4 种工况下，随着开挖面支护力的增大，管片顶部和底部的竖向位移增量基本不变。由图12（c）（d）可知，随着支护力的增大，管片左右两侧点的水平位移增量不断增大，但增大幅度很小。当支护力在0.6P～0.7P之间时，隧道管片变形增量基本一致，因此将支护力维持在0.6P～0.7P之间时，隧道管片变形比较稳定，更容易控制。

图12 顶推力对管片变形的影响（工况2）（E=10 MPa）Fig.12 Influence of supporting force of excavation face on the deformation of pipe segment（case 2）（E=10 MPa）

3.2 隧道横断面管片收敛情况分析

由以上分析可知，工况1 与工况2 隧道收敛变形情况有所不同。为了便于管片变形收敛情况的比较，选取首挖隧道第10 环处管片断面进行分析。分析时分别考虑了有无地下水的影响，2 种工况下管片的收敛变形图如图13所示，收敛变形值见表2。

结合图13与表2中的数据可知，对于工况1，在无地下水情况下，隧道管片整体向右下方偏移，隧道趋于扁平形状，管片顶、底之间的相对位移量为0.07 mm，两侧拱腰之间管片的相对位移值0.08 mm。对于工况2，隧道管片整体向左上方偏移，表现为隧道趋于扁平，管片顶、底之间的相对位移量为0.1 mm，两侧拱腰之间管片的相对位移值为0.08 mm。当存在地下水时，工况1 的管片收敛情况变化很小，而工况2 的收敛变化较大，竖向收敛值达到0.18 mm。因此，对比这几种情况，第一种开挖方式，即“左-右-中”的开挖方式对首挖隧道影响最小。

图13 隧道无地下水时管片收敛变形Fig.13 Convergence deformation of tunnel segment without water

表2 不同工况下管片收敛变形Table 2 Convergence deformation of segments under different conditions /mm

4 结论

（1）首挖隧道为左线隧道时，随着右线和中线隧道的开挖，首挖隧道管片变形增量呈“一小一大”双峰分布，这是由于后续两线隧道距离首挖隧道距离远近的不同导致；首挖隧道为中部隧道时，随着右线和左线隧道的开挖，首挖隧道管片竖直方向变形增量呈现“同侧相同双峰”分布，水平方向变形增量呈现“异侧相同双峰”分布，表明该工况下首挖隧道的管片变形不受两侧隧道开挖次序的影响。

（2）随着注浆弹性模量的增大，隧道管片竖向变形增量不断减小，水平变形增量有一个先增大后减小的变化；随着开挖面支护力的增大，隧道管片的竖向变形增量绝对值不断减少，水平变形增量不断增大。当顶推力在0.6P～0.7P之间时，隧道管片变形增量基本一致，因此将顶推力维持在0.6P～0.7P之间时，隧道管片变形比较稳定，更容易控制。

（3）有地下水存在时，管片变形增量的变化规律基本不变，但管片整体位移较无地下水时增大。

（4）在对比2 种工况隧道开挖完成后的隧道管片收敛情况后，工况1 的收敛情况更小，因此开挖时建议采用左线-右线-中线的开挖方式，降低隧道开挖对已建隧道管片的影响。