基于数值仿真的超大断面扁平结构隧道洞口偏压优化分析

Bias Optimization Analysis of Super-Large Section Flat-Structure Tunnel Entrance Based on Numerical Simulation

摘要

针对超大断面扁平结构矿山法公路隧道洞口段偏压问题展开研究，利用数值仿真技术对所提出的反压回填加固方案进行分析，且在数值模型中考虑了双线隧道不同开挖顺序对偏压问题的影响。数值仿真结果表明：反压回填能够有效解决浅埋偏压隧道的拱顶塌落、不均匀沉降、水平变形大和不对称应力分布等问题，验证了反压回填方案的可行性。综合评价结果表明：反压回填后先开挖深埋侧隧道对围岩变形的影响更小。在此基础上，结合现场施工方案深入分析了超大断面扁平结构隧道的围岩应力分布特征，现场施工监测数据与数值模拟结果具有较高的匹配度，验证了数值模型的准确性，以及所提出优化施工方案的可行性。该研究结果可为超大断面扁平结构隧道洞口偏压段的设计和施工提供理论支撑和参考。

Abstract

This paper examines the issue of asymmetrical pressure at the entrance of a super-large section flat structure mine tunnel.The study employs numerical simulations，utilizing the finite difference software FLAC^3D，to analyze a proposed backfilling scheme.The simulations consider the impact of varying excavation sequences in double-line tunnels on the asymmetrical pressure problem.Consequently，numerical models were developed to examine four distinct excavation strategies：Unbackfilled with the left line excavated first，unbackfilled with the right line excavated first，backfilled with the left line excavated first，and backfilled with the right line excavated first.By analyzing the distribution of plastic zones，as well as the vertical and horizontal displa-cements of the tunnels under these varying excavation strategies，the feasibility of the backfilling approach was assessed and validated.The backfilling scheme effectively addresses the issue of uniform settlement of the surrounding rock under biased conditions，thereby mitigating the horizontal displacement of the tunnel vault and the horizontal convergence of the sidewall，ultimately enhancing the stability of the surrounding rock.The stress distribution within the tunnel’s surrounding rock follows the pattern of “arch foot>sidewall>arch floor>vault” with stress concentrations occurring at the arch foot and the vault being prone to tensile failure，which reflects the stress characteristics of the surrounding rock in super-large section flat structure tunnels.The implementation of a backfilling scheme mitigates issues related to loosening at the vault and asymmetrical stress distribution at the arch footings.Based on an analysis of surrounding rock deformation，the concept of secondary deformation is introduced，allowing for a comparative assessment of the secondary impacts of various excavation sequences on the tunnel.The findings indicate that，for the tunnels under study，the optimal excavation sequence post-backfilling is to first excavate the left line followed by the right line.Conversely，the optimal sequence prior to backfilling involves first excavating the right line and then the left line.The results of the final comprehensive evaluation indicate that initiating excavation on the deeply buried side of the tunnel post-backfilling significantly reduces the deformation of the surrounding rock.Based on numerical simulation outcomes，appropriate reinforcement measures are recommended for practical construction applications.Furthermore，an in-depth analysis of the stress distribution characteristics of the surrounding rocks and supporting structures in super-large section flat structure tunnels is conducted in conjunction with the on-site construction plan.The monitoring data closely align with the numerical simulation results regarding the deformation observed around the periphery of the tunnel’s cross-section.This congruence confirms the accuracy of the numerical model developed in this study and underscores the viability of the proposed optimized construction scheme.The findings presented herein offer theoretical support and serve as a data reference for the design and construction of the entrance bias section in super-large section flat structure tunnels.

Graphical abstract

关键词

超大断面 / 扁平结构 / 双线隧道 / 洞口偏压 / 软弱围岩 / 数值模拟

Key words

super-large section / flat structure / double-line tunnel / entrance bias / weak surrounding rock / numerical simulation

超大断面扁平结构隧道是指开挖断面的面积超过140 m²、扁平比小于0.65（含仰拱时）或小于0.55（不含仰拱时）的椭圆扁平状隧道（曲海锋，2009；彭波，2021）。近年来，我国公路隧道建设数量和规模显著增加，为缓解高密度的交通问题，（超）大断面扁平结构隧道形式得到了广泛应用，如：郑州市韩门隧道，断面面积约为200 m²，扁平比为0.55（皇民等，2022）；广州市北二环高速公路的木强隧道，断面面积为280 m²，扁平比仅为0.33（廖清炯等，2001）；重庆地铁5号线隧道断面面积达388 m²，扁平比为0.61（金星亮等，2022）。超大断面低扁平比的隧道断面形式能够显著提高断面面积的利用率，但同时也加剧了隧道围岩应力分布的复杂性，从而给隧道施工带来了更多的挑战和不利影响，如：隧道拱顶易松弛，降低了围岩承载拱的稳定性；拱脚处应力集中，对两侧地基造成较大的压力作用，支护结构也受到较大的轴力和弯矩作用（段慧玲等，2009；蒋坤等，2019）。尤其对于山岭隧道而言，沟谷或傍山隧道的洞口段往往存在严重的偏压问题，在不平衡的侧向压力作用下隧洞易产生较大的水平变形和不均匀沉降，加之洞口段埋深较浅，围岩软弱，开挖过程中容易引起隧洞坍塌和山体滑坡等安全事故，大幅增加了洞口段的施工难度（高爱民，2018；刘文杰，2020）。因此，分析超大断面扁平结构隧道偏压工况下的受力特征及优化施工技术，对于增强隧道在施工和使用中的安全性具有重要意义。

山岭隧道的偏压问题得到了学者的广泛关注，已有研究通过模型试验（周晓军等，2006；代树林等，2022）、理论计算（邱业建等，2015；杜建明等，2021）和数值模拟（宋战平等，2017；林锦腾等，2020；任纪磊，2023）等方法，对隧道偏压进行了深入分析。偏压荷载主要是由斜坡地形或顺层构造引起的，且随着隧道埋深的增加，偏压影响逐渐减弱（周晓军等，2006；刘小军等，2011）。对于浅埋的偏压隧道，围岩受到的侧向挤压作用明显，因此水平侧向压力系数和水平压力占比系数等考虑侧向压力影响的概念被引入浅埋偏压隧道围岩压力的计算方法中，形成了较为统一的计算模型（刘翔等，2017；徐冲等，2022）。为了提高浅埋偏压隧道的稳定性，研究人员探索了不同的施工方案，指出先开挖浅埋侧和双侧壁导坑法等方案能够有效减小偏压隧道的变形（陈贵红，2012；韩立志，2020）。针对特定地形条件下的浅埋偏压隧道，还可以采取超前注浆支护、反压回填和长锚杆支护等措施来控制隧道变形（李志星等，2011；胡炜等，2020）。此外，也有学者提出通过挡土墙、超前管棚和非对称支护结构等措施来控制偏压隧道中围岩的不规则变形和受力（马杲宇等，2006；潘文韬等，2021；宋战平等，2022）。

综上所述，现有文献从理论计算、形成原因、工程危害和施工工艺等方面对偏压隧道进行了诸多研究，对于常见的隧道偏压问题具有较好的借鉴意义。然而，有关偏压施工优化措施的研究多停留在工程实践层面，对于其作用机理的深入分析，还需结合数值模拟等手段进行针对性研究。此外，当前有关隧道偏压问题的研究，较少涉及超大断面且具有高扁平比形式的隧道，对这种越来越广泛采用的隧道结构与浅埋偏压问题耦合所导致的围岩破坏模式和工程危害的了解不够深入，难以提出适用的施工优化措施，因此对于此类隧道偏压问题的研究亟待完善。本文以深汕合作区通港大道2号隧道作为背景工程，将现场监测数据与数值仿真结果相结合，分析超大断面扁平结构隧道在偏压荷载作用下的围岩压力和变形分布规律，探究不同施工方案的作用机理和可行性，最终提出适用于浅埋偏压超大断面扁平结构隧道施工方案的优化措施，为相关隧道的设计和施工提供理论支撑和数据参考。

1 工程概况

1.1 基本概况

深汕合作区通港大道2号隧道位于广东省深汕特别合作区西南侧、环境科技产业园的西侧，是一条穿越山岭的双线矿山法隧道，南端洞口段两线隧道净距为20.24 m。隧道右线起讫里程桩号为YK6+739.708~YK7+098.708，长度为359 m，左线起讫里程桩号为ZK6+688~ZK7+089，长度为401 m。隧道（单洞）最大净宽为18.762 m，最大净高为11.542 m，断面面积约为232 m²，扁平比为0.615，属于超大断面扁平结构隧道。

1.2 地质概况

隧址地质情况较为复杂，隧洞穿越岩层种类较多，其中洞口段多为强风化围岩，岩石完整性和自稳能力较差，属于Ⅴ级围岩。由图1可知，隧道南端入口处左线埋深远高于右线埋深，属于浅埋偏压段。此外，隧道开挖面积较大，扁平率较低，在偏压作用下结构受力情况复杂，且右线上方覆土较薄，存在大量孤石群，开挖施工易导致覆土松动，发生坍塌且难以施作支护措施。综上所述，本隧道洞口段的结构稳定性受多个不利因素的叠加作用，在偏压的影响下隧道进洞施工存在较大的难度及安全风险。因此，本隧道在洞口段施工的复杂程度具有较好的代表性，针对不同的施工方案，对其开展缓解偏压作用机理和效果分析，对于后续类似工程施工优化措施的选取具有较好的借鉴意义。

1.3 施工方案

综合考虑洞口偏压、软弱围岩和孤石等不利因素，采用回填反压平衡法对洞口进行加固处理，即：隧道开挖工法采用双侧壁导坑法，移除洞口上方孤石后，辅以超前小导管及上覆水泥土回填反压平衡进行施工，同时为初期支护结构预留足够的埋深空间。根据上述方案，洞口段隧道横截面的加固方案设计如图2所示。

2 数值模型

2.1 几何尺寸和材料参数

本研究中采用有限差分程序FLAC^3D建立数值模型，选取左线里程ZK6+740~ZK6+760（对应右线里程YK6+760~YK6+780）作为模拟段。为减少计算时间和降低边界效应的影响，模型左侧边界距左洞洞壁最短距离为60 m，约为6倍洞径，右侧边界距右洞洞壁最短距离为60 m，下侧边界距离拱底30 m，约为3倍洞径，模型上部边界根据实际地形进行重构，双洞埋深分别为五级围岩段的平均埋深（左洞为17.63 m，右洞为2.63 m）。模型的左侧和右侧约束x方向的位移，前方和后方约束y方向的位移，下侧边界约束z方向的位移，上侧地表为自由端，坐标系如图3所示。为验证上述加固施工方案的效果，建立了2套数值模型，即不进行任何加固处理的原始洞口段模型［图3（a）］和反压回填处理后的洞口段模型［图3（b）］，用以对比分析偏压处理施工措施的加固效果。

根据现场勘查结果，隧道南端临近洞口端场地内的揭露地层主要为粉质黏土、强风化和全风化花岗岩，围岩级别为Ⅴ级。回填材料采用5%水泥土，在数值模型中结合现场试验数据和相关文献（曾胜华等，2010；韩立志，2020）选取Ⅴ级围岩和水泥土参数，围岩和回填土均采用实体单元进行划分，服从Mohr-Coulomb破坏准则，具体材料参数见表1。

2.2 模拟方案

首先基于原始洞口段几何模型（即不进行反压回填处理），按照先开挖左线后开挖右线的顺序，先模拟未施加反压回填的原始洞口段开挖（模型“未反压—先左线”），以评估偏压对隧道稳定性的影响，再模拟施加反压回填后的开挖过程（模型“反压—先左线”），用以验证反压回填措施的加固效果。此外，考虑到双线隧道开挖顺序可能影响围岩压力和变形，在上述2种模型的基础上，分别按先挖右线再挖左线的顺序进行模拟（模型“未反压—先右线”和“反压—先右线”），结合数值模拟结果判定缓解偏压问题最有利的开挖顺序。为了直观地体现偏压问题对隧道围岩产生的影响，并兼顾数值模型的计算效率，上述模型仅考虑隧洞成型后的计算结果，忽略了施工过程和支护措施。图4所示为上述4种模拟方案，在图中圆点位置处设置了数据监测点。

3 数值模拟结果分析

3.1 塑性区分布

由模拟结果的塑性区分布（图5）可知，如不进行反压回填处理，开挖后左线隧洞塑性区主要分布在洞壁和拱脚处，而右线隧洞开挖后塑性区主要集中在左下方拱脚和拱顶的右上方，且由于右线隧洞上方覆土厚度较浅，拱顶的塑性区由隧道断面边界连续扩展至地面（图5，模型“未反压—先左线”），因此在偏压条件下右线隧道拱顶的薄弱破碎围岩容易引发隧洞坍塌。

经过反压回填后，右线隧洞的塑性区范围明显减少，拱顶处的塑性区也并没有连续扩展至地表（图5，模型“反压—先左线”）。通过对比先开挖左洞和先开挖右洞2种情况，发现不同开挖顺序对塑性区分布影响不大。

3.2 竖直位移分析

通过对比不同模拟方案的竖直位移云图（图6），发现左洞的拱顶沉降变形和拱底隆起变形均远高于右洞，不同施工方案对拱顶沉降的影响较小。选取整条隧洞中轴线上5个监测点记录的拱顶沉降量的平均值作为对比，可知经过反压处理后的左洞拱顶沉降基本保持不变，降幅仅为0.9%（图7）。而右线由于覆土厚度较薄，拱顶沉降一直处于较低水平，在反压回填后降幅约为17.6%。

然而，偏压作用会导致右线隧洞右侧具有较大范围的隆起区域（图6中的放大区域），说明右线隧洞左右两侧的沉降不对称，极有可能导致隧洞往沉降大的一侧（左侧）倾斜，使隧洞轮廓发生倾斜倒塌变形。经过反压回填后，隆起区域明显减少，右线隧洞左右两侧围岩的竖向位移趋于对称分布，缓解了偏压的作用效果。

3.3 水平位移分析

选取整条隧洞中轴线上拱顶处数据监测点记录的水平位移的平均值作为对比（图8），可知若不采取反压回填措施，在偏压条件下，双线隧洞的拱顶围岩均会产生较大的沿x轴正方向的水平位移。而回填水泥土的重力产生反向挤压效果，在一定程度上平衡了挤压效果，使得左线隧洞的拱顶水平位移分别从1.00 mm和1.01 mm减少至0.49 mm和0.53 mm，降幅近一半；而右线隧洞由于受到水泥土重力的反向挤压作用，产生了沿x轴负方向的位移。由此可见，反压回填措施能够有效抵消偏压力的作用，限制拱顶围岩的水平位移。此外，开挖顺序对水平位移的影响并不明显。

根据隧道开挖后洞壁两侧数据监测点记录的向洞内收敛的水平位移值（图9），可计算出隧洞拱腰处水平收敛位移量。由图9可知，反压回填对拱腰水平收敛量产生的影响较为显著，表现为：先开挖左洞和先开挖右洞2种情况下右线隧洞的水平收敛分别从3.83 mm和3.78 mm降低至2.95 mm和2.91 mm，降幅约为23%，而左线隧洞的水平收敛分别从2.46 mm和2.45 mm增加至3.32 mm和3.31 mm，增幅约为35%。上述变化特征说明进行反压回填后，偏压力的作用效果被削弱，双线隧洞的收敛变形量趋于一致，且避免了右线隧洞过大的水平收敛。

3.4 围岩应力分布特征

从布置的监测点（图4）中选取拱顶、拱底、拱腰及拱脚处的最大主应力值进行分析。由表2可知，4种模型中的最大主应力大小均呈现出“ 拱脚>拱腰>拱底>拱顶”的变化规律，且拱脚处的应力相较其他部位高出1~2个数量级，因此需加强拱脚处的支护。而拱顶处压应力相对较小，有时甚至处于受拉状态，导致围岩易发生松弛，因此对拱顶处的围岩进行加固（如超前小导管注浆）是十分必要的。若不进行反压回填，右线隧洞的左右侧拱脚应力值相差较大，左侧拱脚处最大主应力比右侧拱脚处大67.9%（先开挖左洞）和67.7%（先开挖右洞），即左侧围岩承受更多的荷载，对围岩的稳定性构成不利影响。经过反压回填后，两侧拱脚应力相差百分比缩小至11.3%和11.7%，上述变化均有利于维持围岩的稳定性。此外，通过监测发现，开挖顺序对最大主应力的影响较小。

3.5 不同开挖顺序影响分析

通过分析不同开挖顺序模型中各分析步的数据，发现不同开挖顺序对隧道围岩的影响主要发生在交错开挖过程中，后挖隧洞会影响先挖隧洞的围岩变形和支护结构的受力情况。为了分析后挖隧洞的开挖过程对先挖隧洞周边围岩变形的影响程度，将先开挖隧洞受到后开挖隧洞的影响所产生的变形量记为“二次变形量”{\gamma }_i，将不同位置处二次变形量的总和\gamma作为反映隧洞周边变形受影响程度的指标，即：

{\gamma }_i=|a_{\mathrm{2}-}{a}_{\mathrm{1}}|(i=\mathrm{1,2},\mathrm{3}\dots )（1）

\gamma ={\sum }_{i=\mathrm{1}}^n{\gamma }_i（2）

式中：a_{\mathrm{1}}为第一步开挖后，先挖隧洞的某一项指标数值（如拱顶沉降位移）；a_{\mathrm{2}}为第二步开挖后，先挖隧洞相应的指标数值。

本文选取拱顶沉降、拱顶水平位移和拱腰水平收敛3个指标进行统计（即n=3）。4种模型的各项指标数值，以及相应的二次变形量{\gamma }_i和二次变形量总和γ均列于表3中。

由表3可知，如不进行反压回填，先挖左线隧洞的二次变形量总和比先挖右线隧洞的二次变形量总和大，即先开挖右线隧洞相对更安全，这与以往研究得出的先开挖浅埋侧隧洞的结论相符（刘小军等，2011）。若采取反压回填措施，先开挖左线隧洞的二次变形量总和小于先开挖右线隧洞的情况。

3.6 综合评价

基于上述分析结果，不同施工方案对隧道拱顶沉降的影响较小，而对拱顶水平位移、拱腰水平收敛和左右侧拱脚主应力差值的影响较大，因此，最终的方案决策需根据上述4种受影响较大的指标进行综合评价。为避免主观因素的影响，本文采用基于客观权重的熵权法对各方案进行综合评价。各指标的具体数值（取各项指标的最大绝对值）如表4所示。

熵权法计算流程如下：

（1）根据表4数据建立原始数据矩阵 X （m个指标，n个方案），按照式（3）对其进行标准化处理，得到标准化后的数据矩阵 Y。

y_{ij}=\left\{\begin{array}{c}\frac{x_{ij}-\mathrm{m}\mathrm{i}{\mathrm{n}}_j\left\{x_{ij}\right\}}{\mathrm{m}\mathrm{a}{\mathrm{x}}_j\{x_{ij}\}-\mathrm{m}\mathrm{i}{\mathrm{n}}_j\left\{x_{ij}\right\}},\mathrm{正向}\mathrm{型指}\mathrm{标}(\mathrm{数值}\mathrm{越大}\mathrm{越优})\\ \frac{\mathrm{m}\mathrm{a}{\mathrm{x}}_j\left\{x_{ij}\right\}-x_{ij}}{\mathrm{m}\mathrm{a}{\mathrm{x}}_j\{x_{\mathrm{i}\mathrm{j}}\}-\mathrm{m}\mathrm{i}{\mathrm{n}}_j\left\{x_{ij}\right\}},\mathrm{负向}\mathrm{型指}\mathrm{标}(\mathrm{数值}\mathrm{越小}\mathrm{越优})\end{array}\right.（3）

（2）计算第j项指标中第i个方案的值所占的比重，表示为

p_{ij}=\frac{y_{ij}}{\overline{{\sum }_{i=\mathrm{1}}^n{y}_{ij}}}（4）

（3）计算第j项指标的熵值，表示为

e_j=-\frac{\mathrm{1}}{\mathrm{l}\mathrm{n}(n)}{\sum }_{i=\mathrm{1}}^n{p}_{ij}ln(p_{ij})（5）

式中：满足\frac{\mathrm{1}}{\mathrm{l}\mathrm{n}(n)}>\mathrm{0}，e_j>\mathrm{0}。当p_{ij}=\mathrm{0}时，令p_{ij}\mathrm{l}\mathrm{n}(p_{ij})=\mathrm{0}。

（4）计算第j项指标的信息熵冗余度，表示为

d_j=\mathrm{1}-e_j（6）

（5）计算各指标的权重，表示为

w_j=\frac{d_j}{\overline{{\sum }_{j=\mathrm{1}}^m{d}_j}}（7）

（6）计算综合评分，表示为

F_i={\sum }_{j=\mathrm{1}}^m{y}_{ij}{w}_j（8）

综合评分越高，说明方案效果越优。计算得到的各指标客观权重如表5所示，各方案综合评分如表6所示。

由表6可知，反压回填的方案明显优于未反压回填的方案，且反压回填方案中，先挖左线的效果更佳，更能保证围岩的稳定性。

4 现场监测和实况模型

4.1 模拟方案

为了进一步指导实际施工，并验证数值模拟结果，以模型反压—先左线为基础，依据实际现场施工情况，进一步细化模型，具体包括加入超前支护和初期支护措施，以及对洞口段开挖工法（双侧壁导坑法）进行全过程的数值仿真，同时在施工现场开展了隧道拱顶沉降和水平收敛的监测，用于验证数值模型的合理性。细化后的模型（实际开挖模型）如图10所示，开挖和支护工序模拟如图11所示，模型中双侧壁导坑法的施工和支护步序详见表7。在模型中，超前支护的加固作用是通过加强注浆区域的材料属性来实现的（赵立财，2014），加强后的材料参数见表8。

4.2 位移分析

选取左线ZK6+750（对应右线YK6+770）断面上拱顶沉降和拱腰水平收敛数据的现场监测值与数值仿真结果进行对比，结果如图12所示。由图12可知，数值模拟结果与实际监测数据接近，左线的拱顶沉降误差为0.34 mm，水平收敛误差为0.89 mm，右线模拟结果的拱顶沉降和水平收敛与实际监测数据几乎一致，说明细化后的模型能够较为准确地反映实际施工情况，进而对施工过程进行实时指导和反馈。此外，在数值模拟结果中，左右线的水平收敛基本一致，说明反压回填措施有效削减了偏压作用（如图9“未反压”模型结果），而实际监测中由于围岩复杂的非均匀分布、节理发育程度、地下水分布和施工质量等因素，导致左线与右线的水平收敛存在一定的差距（左线为2.4 mm，右线为1.6 mm）。

4.3 围岩应力分析

选取左线ZK6+750（对应右线YK6+770）断面上隧道各部位的最大主应力并进行绘图，如图13所示。由图13可知，左线隧洞和右线隧洞最大主应力的最大值仍出现在拱脚处，拱顶和拱底的最大主应力较小。结合表2分析可知，相较一次挖通隧洞且无支护时的“反压—先左线”模型结果，采用双侧壁导坑法和及时的初期支护措施后，拱顶处的最大主应力大小分别增加了83.62 kPa（左线）和49.61 kPa（右线），拱底处的最大主应力大小也增加了41.12 kPa（左线）和32.11 kPa（右线），降低了拱顶和拱底处松弛和受拉的风险；而拱脚处的最大主应力均减小至250 kPa以下，缓解了拱脚处的应力集中问题。

此外，隧道初期支护结构受力情况如图14所示，初喷混凝土受拉区域集中分布在拱脚和左洞拱肩处，最大拉应力为1.45 MPa，小于C25混凝土的抗拉强度（约为2.7 MPa），受拉破坏风险较低。拱腰处钢拱架压应力集中，最大压应力为36.67 MPa，小于钢拱架抗压强度（约为215 MPa）。拱脚处锚杆处于受压状态，拱顶处锚杆处于受拉状态，且由于埋深较浅，锚杆的最大轴力仅为30.89 N（受拉），远小于D25中空注浆锚杆的极限拉力（约为180 kN）。上述情况说明，实际施工采用的双侧壁导坑开挖法和初期支护措施在本工程超大断面扁平结构隧道开挖作业中，能够减小对围岩的干扰，有利于维持围岩的稳定性。

5 结论

针对超大断面扁平结构隧道洞口段的偏压问题，提出了反压回填加固方案，通过数值模拟探究了反压回填措施的作用机理，验证了反压回填措施的可行性以及先挖左线的合理性。数值模拟结果与实际监测数据高度匹配，为实际施工提供了数据参考和合理建议。本研究得出如下结论：

（1）在浅埋偏压影响下，隧洞拱顶产生较大的水平位移，浅埋侧洞口由于覆土较薄，易受损，导致隧洞两侧沉降不均匀，进而造成两侧拱脚应力相差较大，易引起隧洞向沉降大的一侧倾斜。

（2）相比于未反压回填的结果，反压回填后隧道的拱顶水平位移减少了近一半，平衡了左右隧洞的水平收敛，以及浅埋侧隧洞两侧的不均匀沉降和不对称拱脚应力，增强了围岩稳定性，同时回填区域为初期支护结构提供了足够的埋深空间。

（3）本文所研究隧道的围岩应力分布呈现“拱脚>拱腰>拱底>拱顶”的规律，在拱脚处压应力集中，拱顶易松弛塌落，需要对隧道上方围岩进行加固（如超前小导管注浆）。

（4）根据本文提出的“二次变形量”的对比结果和综合评价结果，反压回填后先开挖深埋侧的方案效果最优。

参考文献

[1] Chen Guihong，2012.A study of the construction method for a multi-arch tunnel with unsymmetrical pressure［J］.Modern Tunnelling Technology，49（1）：84-88，95.

[2] Dai Shulin， Yang Youyuan， Xin Yulong，et al，2022.Model experimental investigation on the effect of terrain bias on highway tunnel bias degree［J］.Science Technology and Engineering，22（7）：2901-2906.

[3] Du Jianming， Fang Qian， Lu Hai，et al，2021.Calculation method for surrounding rock pressure of shallow tunnel with asymmetrical pressure of variable slopes［J］.Journal of Central South University（Science and Technology），52（11）：4088-4098.

[4] Duan Huilin， Zhang Lin，2009.Comparative study of rational excavation methods for large-span highway tunnels［J］.China Civil Engineering Journal，42（9）：114-119.

[5] Gao Aimin，2018.Statistics of construction safety accidents and analysis of risk factors in tunnel portal section［J］.Railway Energy Saving and Environmental Protection and Occupational Safety and Health，8（4）：220-224.

[6] Han Lizhi，2020.Numerical simulation and research on excavation of shallow buried biased tunnel with both sides heading method［J］.Highway，65（5）：324-328.

[7] Hu Wei， Yu Yu， Tan Xinrong，et al，2020.Study on failure characteristics and supporting measures of tunnel embedded in geologically bedding strata under unsymmetrical pressure［J］.Chinese Journal of Underground Space and Engineering，16（Supp.1）：450-456，510.

[8] Huang Min， Jin Muyan，2022.College on the deformation predi-ction and control of weak rock using core soil method in Hanmen large section tunnel［J］.Journal of Henan University of Engineering（Natural Science Edition），34（2）：19-23.

[9] Jiang Kun， Zhang Bingbing，2019.Optimization analysis on flatness ratio of oversized-section tunnel［J］.Urban Roads Bridges and Flood Control，（8）：291-294.

[10] Jin Xingliang， Hu Xibo， Liang Bin，2022.Study on force characteristics of composite lining structure of shallow buried large cross section flat tunnel［J］.Guangdong Architecture Civil Engineering，29（2）：63-65.

[11] Li Zhixing， Zhang Yuhua， Shi Xiaomin，2011.Excavation technology of shallow tunnel under unsymmetrical pressure［J］.Construction Technology，40（3）：60-62.

[12] Liao Qingjiong， He Feng，2001.Construction method of biased shallow buried section of Muqiang Tunnel［J］.Highways and Automotive Applications，（5）：39-40.

[13] Lin Jinteng， Qiu Hongzhi， Wang Yongzhen，2020.Deformation analysis of surrounding rock in super shallow buried section of mountain tunnel［J］.Science Technology and Engineering，20（24）：9982-9987.

[14] Liu Wenjie，2020.Analysis and identification of risk causes in shallow buried bias tunnel construction［J］.Sichuan Architecture，40（3）：320-322.

[15] Liu Xiang， Fang Qian，2017.Calculation method for loose surrounding rock pressure of shallow tunnel under asymmetrical pressure on variable slope［J］.China Railway Science，38（2）：64-70.

[16] Liu Xiaojun， Zhang Yongxing，2011.Analysis of reasonable excavation sequence and stress characteristics of portal section of shallow tunnel with unsymmetrical loadings［J］.Chinese Journal of Rock Mechanics and Engineering，30（Supp.1）：3066-3073.

[17] Ma Gaoyu， Chen Ziquan， Yang Wenbo，et al，2019.Analysis on mechanical properties of supporting structure of shallow buried tunnel passing through loose accumulation body under asymmetric load［J］.Subgrade Engineering，（4）：20-26.

[18] Pan Wentao， Wu Fangyin， He Chuan，et al，2021.Construction methods and asymmetric design optimization for shallow-buried tunnels subjected to unsymmetrical loads［J］.Tunnel Construction，41（Supp.1）：352-361.

[19] Peng Bo，2021.Research on numerical simulation and optimi-zation of excavation method for flat structure tunnel with super-large section［J］.Northern Communications，（10）：78-82.

[20] Qiu Yejian， Peng Limin， Lei Mingfeng，2015.Upper bound solutions for surrounding rock pressure of shallow bias tunnel［J］.China Civil Engineering Journal，48（6）：106-113.

[21] Qu Haifeng，2009.An overview on construction and research of super-large cross-section tunnels with low height-span ratios［J］.Tunnel Construction，29（2）：166-171.

[22] Ren Jilei，2023.Construction optimisation and force analysis of entrance sections of large-section shallow bias tunnels［J］.Transpo World，（21）：161-163.

[23] Song Zhanping， Pei Jiafeng， Pan Hongwei，et al，2022.Analysis on the support effect of advance large pipe umbrella at portal section of shallow tunnel under unsymmetrical load［J］.Modern Tunnelling Technology，59（6）：86-96.

[24] Song Zhanping， Wang Tong， Zhou Jianjun，et al，2017.Analysis of construction optimization and mechanical characteristic of shallow large section tunnel［J］.Chinese Journal of Underground Space and Engineering，13（2）：459-468.

[25] Xu Chong， Wang Qian，2022.A unified model and calculation method of surrounding rock pressure for shallow buried tunnel［J］.Journal of Railway Engineering Society，39（8）：46-52.

[26] Zeng Shenghua， Zeng Juan，2010.Experimental study on stren-gth property of low proportional cement soil［J］.Subgrade Engineering，（4）：17-19.

[27] Zhao Licai，2014.Research on excavation effect of grouting pipe-roof in tunnel passing through water-rich weak zone［J］.Chinese Journal of Underground Space and Engineering，10（Supp.2）：1853-1859.

[28] Zhou Xiaojun， Gao Yang， Li Zelong，et al，2006.Experimental study on the uneven rock pressure and its distribution applied on a tunnel embedded in geologically bedding strata［J］.Modern Tunnelling Technology，（1）：12-21.

[29] 陈贵红，2012.偏压连拱隧道施工方法研究［J］.现代隧道技术，49（1）：84-88，95.

[30] 代树林，杨有源，辛毓龙，等，2022.地形偏压对公路隧道偏压程度影响模型试验研究［J］.科学技术与工程，22（7）：2901-2906.

[31] 杜建明，房倩，海路，等，2021.地表变坡下浅埋偏压隧道围岩压力计算方法［J］.中南大学学报（自然科学版），52（11）：4088-4098.

[32] 段慧玲，张林，2009.大跨度公路隧道合理开挖方法对比研究［J］.土木工程学报，42（9）：114-119.

[33] 高爱民，2018.隧道洞口段施工安全事故统计与风险因素分析［J］.铁路节能环保与安全卫生，8（4）：220-224.

[34] 韩立志，2020.浅埋偏压隧道双侧壁导坑法开挖数值模拟研究［J］.公路，65（5）：324-328.

[35] 胡炜，喻渝，谭信荣，等，2020.顺层偏压隧道破坏特征及支护措施研究［J］.地下空间与工程学报，16（增1）：450-456，510.

[36] 皇民，金沐妍，2022.韩门大断面隧道预留核心土法施工中软弱围岩变形预测与控制［J］.河南工程学院学报（自然科学版），34（2）：19-23.

[37] 蒋坤，张兵兵，2019.超大断面隧道扁平率优化分析［J］.城市道桥与防洪，（8）：291-294.

[38] 金星亮，胡曦波，梁斌，2022.浅埋超大断面扁平隧道复合式衬砌结构力学特性分析［J］.广东土木与建筑，29（2）：63-65.

[39] 李志星，张宇华，石效民，2011.浅埋偏压隧道开挖技术［J］.施工技术，40（3）：60-62.

[40] 廖清炯，何峰，2001.木强隧道偏压浅埋段施工方法［J］.公路与汽运，（5）：39-40.

[41] 林锦腾，邱洪志，王永祯，2020.山岭隧道超浅埋段开挖围岩变形分析［J］.科学技术与工程，20（24）：9982-9987.

[42] 刘文杰，2020.浅埋偏压隧道施工风险源分析及其识别方法［J］.四川建筑，40（3）：320-322.

[43] 刘翔，房倩，2017.变坡面浅埋偏压隧道松动围岩压力计算方法［J］.中国铁道科学，38（2）：64-70.

[44] 刘小军，张永兴，2011.浅埋偏压隧道洞口段合理开挖工序及受力特征分析［J］.岩石力学与工程学报，30（增1）：3066-3073.

[45] 马杲宇，陈子全，杨文波，等，2019.浅埋偏压隧道穿越松散堆积体支护结构受力特性分析［J］.路基工程，（4）：20-26.

[46] 潘文韬，吴枋胤，何川，等，2021.浅埋偏压隧道施工工法研究与非对称设计优化［J］.隧道建设（中英文），41（增1）：352-361.

[47] 彭波，2021.超大断面扁平结构隧道开挖方法数值模拟及优化研究［J］.北方交通，（10）：78-82.

[48] 邱业建，彭立敏，雷明锋，2015.浅埋偏压隧道围岩压力上限法解析解［J］.土木工程学报，48（6）：106-113.

[49] 曲海锋，2009.扁平特大断面隧道修筑及研究概述［J］.隧道建设，29（2）：166-171.

[50] 任纪磊，2023.大断面浅埋偏压隧道洞口段施工优化及受力分析［J］.交通世界，（21）：161-163.

[51] 宋战平，裴佳锋，潘红伟，等，2022.浅埋偏压隧道洞口段超前大管棚支护效果分析［J］.现代隧道技术，59（6）：86-96.

[52] 宋战平，王童，周建军，等，2017.浅埋偏压大断面隧道施工优化及受力特征分析［J］.地下空间与工程学报，13（2）：459-468.

[53] 徐冲，王乾，2022.浅埋隧道围岩压力的统一表达模型及计算方法［J］.铁道工程学报，39（8）：46-52.

[54] 曾胜华，曾娟，2010.低掺量水泥土强度特性试验研究［J］.路基工程，（4）：17-19.

[55] 赵立财，2014.注浆管棚在隧道穿越富水软弱带开挖效应研究［J］.地下空间与工程学报，10（增2）：1853-1859.

[56] 周晓军，高杨，李泽龙，等，2006.地质顺层偏压隧道围岩压力及其分布特点的试验研究［J］.现代隧道技术，（1）：12-21.

基金资助

上海隧道工程有限公司专项研究科研项目“软弱地层超大断面扁平结构矿山法隧道施工关键技术研究”(2022-sk-19)