宋战平,王凯蒙,王 涛,王军保,唐坤尧
(1. 西安建筑科技大学 土木工程学院,陕西 西安 710055;2. 西安建筑科技大学 陕西省岩土与地下空间工程重点实验室,陕西 西安 710055;3. 中国铁建大桥工程局集团有限公司,天津 300300)
近年来,伴随着国家城镇化进程的不断加速,城市地铁隧道下穿既有建(构)筑物在工程建设中不可避免.国内外的专家学者结合具体的工程实例,在理论分析、数值计算、试验研究以及现场监测等方面对隧道的开挖方法和辅助施工措施进行了大量的研究,取得了卓越的研究成果.
袁驰宇[1]使用数值模拟软件分析研究了山区浅埋地铁隧道下穿既有运营高速公路的四种常用施工工法,并通过对比数值模拟和隧道监控量测的数据,推荐浅埋隧道下穿施工采用CRD法.张著[2]结合背阴坡隧道现场监控量测结果,模拟并对比了不同施工方法的下穿隧道有限元计算模型,结果表明山区公路路基的稳定性随着开挖进尺的减小而提高,侧壁导坑法效果最好,其次是台阶法,全断面法对稳定性影响最大.研究了下穿隧道位置对于公路路基地表沉降的影响;刘松涛[3]等以湘桂铁路石头岗隧道工程实例,使用数值分析软件对多种施工方法进行了分析研究,并结合实际施工结果,验证了上部中隔墙法较好地控制了高速公路路面沉降.在理论分析方面,为解决浅埋隧道周边应力场及位移场的精确性,提出了复变函数法、随机介质理论等计算方法,结合收集到的监控量测数据,得到了符合工程实际情况的解析解[4-6].但从相关研究文献可以看出,由于现场情况复杂,理论分析有着局限性,对于山区浅埋隧道下穿既有构筑物的绝大多数研究是通过测量人员定期监测现场布置的监控测量点或使用自动化采集装置收集数据,通过监控测量隧道关键断面及部位在开挖过程中的沉降值以及收敛位移值.施工方法主要是选用双侧壁导坑法、CRD法、CD法等在控制隧道沉降和收敛方面比较好的工法,保证施工和隧道的安全,同时结合地质超前预报、大管棚超前支护等施工方法控制隧道以及既有建(构)筑物的位移变形情况[7-10].但是,采用监测方式进行研究存在着不足,监测点受地形及施工影响布置不够精准,人员或自动采集装置收集的监测数据通常存在着一定地滞后性及误差性等,隧道在开挖施工过程中不能对围岩变形情况以及施工附近区域的地表沉降、原有建筑物结构的安全进行精准控制.且当前浅埋隧道下穿既有构筑物的数值模型多为二维平面模拟,由于实际情况较为复杂,在建立模型时不得不进行一定程度的简化,而且计算时没有考虑在隧道开挖和支护过程中的三维空间效应[10-13].
本文以贵阳地铁1号线下麦西车站YD1K1+395~+505段浅埋地铁隧道下穿既有高速公路项目为工程背景,使用Midas GTS/NX有限元软件建立三维有限元计算模型,通过模拟隧道在施工时采用四种不同的施工工法,同时考虑软件里隧道开挖过程中的三维空间效应,比较不同工法下的位移场变化和支护结构受力情况,得到模拟施工的最优工法.运用数值模拟软件计算不同施工工法下隧道变形及高速路地表沉降量,并结合现场监控量测数据分析发现CRD法为最优施工工法.本文为西南山区浅埋地铁隧道的施工安全和正常运营的高速公路的行车安全予以技术指导,同时进一步为西南山区浅埋地铁隧道下穿既有高速公路及交通涵洞提供施工上的理论支持.
下穿环城高速公路段浅埋地铁隧道全长528 m,线路设计为单洞双线,其中共有110 m隧道下穿高速公路路堤及涵洞出口,隧道与高速公路成50°夹角相交,如图1所示.其中:YD1K1+395~+456段共61 m下穿公路路堤,YD1K1+456~+505段共49 m下穿涵洞出口地段及附近地面,其平面布置见图2.该段浅埋隧道下穿高速公路所处的地质条件较差,岩土勘察报告显示下穿段地质主要为页岩和第四系土层,岩石破碎程度大属于V级围岩.高速公路路面距离下穿地铁隧道拱顶高度为23 m,涵洞形式为拱涵,其洞口孔径和高度为别为4 m、5 m,翼墙受地质条件所限露出地表处的高度为2~5 m,翼墙墙背后为高速公路路堤.该段下穿地铁隧道埋深较浅,与涵洞
图1 隧道位置图Fig.1 The location of the tunnel
图2 隧道—高速路—涵洞位置关系示意图Fig.2 The relationship between tunnel, highway, and culvert
处的最小水平距离只有2.2 m,下穿隧道在涵洞出口处最小埋深只有0.8 m,在涵洞出口右侧20 m范围内埋深在2.5 m以内.因为涵洞为通向拌合厂的出口,为保证施工进度,不能封闭路段.
下穿地铁隧道线路方向与既有高速公路成50°角斜交,通过建立隧道的三维几何模型,反映出隧道和高速公路的实际位置关系,同时也能更好和高效地分析隧道施工过程中对高速公路地面的影响.同时考虑到现场地质条件的复杂性,模型对土层和建筑材料进行了一定简化,将隧道围岩和土体都视为各向同性体,屈服准则采用Mohr-Column准则,使用3D实体单元进行模拟,土体参数来源于地质勘察报告.在实际施工中,锚杆使用钢筋制作,初次支护、二次衬砌主要为钢筋拱架和混凝土,故屈服准则均为线弹性,采用2D单元划分,初期支护的等效弹性模量使用模拟法计算得出[14-15].计算参数详见表1.
同时为了尽可能减少使用模拟软件建立模型时产生的尺寸效应,减轻模型边界约束的影响,提高软件分析计算时的精度,取模型的边界为地铁隧道跨度的3~5倍.计算模型在X、Y、Z方向上分别取171 m×100 m×100 m,即为实际工程的YD1K1+327~YD1K1+498段,具体的隧道模拟模型如图3所示.计算模型边界约束使用软件提供的自动约束功能添加,模型四周和地面为法向约束,地表不添加约束为自由面.高速公路车辆荷载的取值按规范取公路I级荷载10.5 kN/m.
表1 模型的力学参数
图3 计算模型Fig.3 The model of computation
山区浅埋下穿隧道施工时,通常选用的施工工法为:CD法、CRD法、双侧壁导坑法这类强调小断面强支护的工法.本文在以上三种工法的基础上增加了台阶法进行了模拟及计算.为了尽量减少数值软件模拟的边界条件在计算时对结果的影响,选取隧道中间断面即X=85.5 m处进行分析比较,从而得出最有效且适合的施工方法.
运用Midas软件模拟施工时,台阶法施工步骤为:先施工上台阶部分,由于隧道围岩是V级围岩,上台阶长度考虑取1~2倍隧洞宽,模拟时上台阶每轮进尺取20 m,再施工下台阶,上下台阶间距取20 m,最后施作仰拱部分;CD法施工步骤为:先开挖1-3隧道左侧断面,然后开挖4-6隧道右侧断面;CRD法施工步骤为:隧道左侧断面和隧道右侧断面交替开挖;双侧壁导坑法施工与CD法类似,其步骤为:先施作隧道左侧超前导洞,再施作隧道右侧导洞,在此期间保留中间核心土区域,最后再进行隧道中间断面的开挖施工.除台阶法外,模拟其余三种施工方法时,每一施工步的进尺取3 m.具体步骤见图4.
图4 观测点位置图Fig.4 Location of observation points
3.1.1 垂直方向位移
由表2和图5可以发现,在四种不同的施工工法下,隧道施工至相同地点时,隧道四周各观测点竖向位移的监测数据:双侧壁导坑法的拱顶沉降值在完成施工后为-3.217 cm,左、右拱肩的沉降值分别为-0.821 cm和-2.303 cm,通过比较可知,其各监测点的沉降值是四种施工工法中最小的.其次是CRD法和CD法,比较8个监测点的沉降值,这两种工法的沉降值大小相当,台阶法的沉降值最大,拱顶沉降达到了-3.995 cm.从左右边墙处的竖向位移值来看,四种工法近乎没有差别,左边墙位移值在1 cm左右,右边墙位移在-0.2 cm左右.隧道底部的情况与拱顶相似,采用台阶法施工时隧道拱底隆起值最大,达到了6.640 cm,CRD法和CD法较小,而且两种施工方法在监测点所得到的位移值基本没有差异,隆起值都在6.5 cm左右,双侧壁导坑法最小,隆起值只有6.074 cm.
通过分析观测点所得的拱顶沉降值,其余三种施工方案同双侧壁导坑法相比,台阶法的沉降值增大了24.2 %,CD法和CRD法的沉降值分别大了19.8 %、17.8 %;而在拱底隆起值方面,三种施工方法分别增大了9.3 %、7.3 %、7.2 %,工法排序与上文相同.可以看出,在控制隧道拱顶沉降和拱底隆起方面,台阶法都不能取得理想的施工效果,CD法和CRD法较为接近双侧壁导坑法,且两种工法在各观测点所得到的数据相差不大.通过比较可知,采用双侧壁导坑法施工效果最好,其次是CRD法和CD法,考虑到施工复杂程度和经济成本,推荐采用CRD法施工.
图5 不同工法下各观测点竖向位移图Fig.5 The vertical displacement of observation points under different construction methods
台阶法CD法CRD法双侧壁导坑法拱顶1-3.995-3.854-3.790-3.217左拱肩2-1.459-1.070-1.069-0.821右拱肩3-2.942-2.968-2.935-2.303左边墙41.0381.1041.1140.952右边墙5-0.218-0.171-0.158-0.218左拱脚63.5633.3463.3462.829右拱脚72.6662.7502.7422.108拱底86.6406.5226.5146.074
3.1.2 水平方向位移
表3和图6为浅埋隧道在四种施工工法下,模拟隧道施工至指定位置并施作支护措施后,得到的8个监测点水平位移的具体数值情况:通过比较可知,台阶法、CD法、CRD法这三种施工工法引起的水平位移值都比双侧壁导坑法要小,综合比较这三种施工工法,台阶法引起的水平位移值在8个监测点中,只有左边墙4号点的数值小于其余两种工法,其余各观测点的水平位移均大于CD法、CRD法.CRD法所引起的水平位移值是四种工法中最小的,CD法次之.
表3 各施工方法各观测点水平位移情况表(单位:cm)
图6 不同工法下各观测点水平位移图Fig.6 The horizontal displacement observation points under different construction methods
分析图6可知,8个监测点中,只在监测点左拱肩2、拱底8处,双侧壁导坑法施工所导致的水平位移值比台阶法引起的水平位移值小,而在其余监测点处,双侧壁导坑法控制水平位移的效果都是最差的.为了更直观的比较四种工法在控制隧道水平位移时的优劣,选择左边墙4号监测点作为研究对象.双侧壁导坑法在此观测点的水平位移值是8个监测点上最大的,位移值达到了0.826 cm.在这一监测点上,其他三种施工工法同双侧壁导坑法相比较而言:台阶法减小了21.4 %,CD法减小了17.1 %,CRD法减小了14.8 %.
通过上文的分析可知,四种施工工法在控制隧道水平收敛方面,双侧壁导坑法和台阶法表现出的控制效果较差,分析其原因,由于台阶法施工时掌子面开挖面积较大,施工对隧道扰动比较大,造成水平位移较大.双侧壁导坑法分部小断面保留核心土开挖,保证了施工时隧道的稳定,但其施工步骤较多,且工艺复杂,造成水平位移不能得到较好的控制.CD法和CRD法控制效果无明显差异,使用CD法开挖时造成的水平位移平均值比双侧壁导坑法小了20 %左右,说明这两种施工方法在控制隧道水平位移方面较好.
3.1.3 地表沉降
图7为四种开挖工法的沉降曲线,四条沉降曲线都是两边沉降值小,中间沉降值大,表明地面有沉降形成.对比分析各施工方法,地表沉降最大值出现在5号观测点,台阶法的沉降最大,沉降值达到了-1.371 cm.在5号观测点,其他三种工法同台阶法进行比较可知,监测点处的地表沉降值分别减小了15.2 %、17.6 %、30.1 %.其中,双侧壁导坑法在观测点CD1处沉降最小,只有-0.517 cm.控制高速公路地表沉降方面,四种施工工法都可以达到规范控制要求,双侧壁导坑法在控制地表沉降方面效果最好,其次是CRD法和CD法,CRD法的沉降值相较而言更小,这两种工法与双侧壁导坑法在相同监测点上的沉降值相差很小,台阶法的效果最差.
控制地表沉降是下穿施工的一个重点,模拟结果表明,四种施工方法引起的地表沉降在0.517~1.371 cm之间,对高速公路的影响轻微.综合隧洞垂直、水平位移的模拟结果,施工方法的最优选择是双侧壁导坑法,但考虑到现场施工进度和工程造价,选择CRD法施工符合现场实际情况,且施工效果较好.
表4 不同工法下地表各测点沉降值(单位:cm)
图7 不同工法下高速公路各测点沉降曲线图Fig.7 Settlement curve of measuring points under different construction methods
图8为四种施工工法在隧道掌子面开挖完成后的剖面图,其塑性区的分布无明显差异,只是塑性区的占比和数值大小略有差异.各工法施工后的塑性区大小占模型截面总面积的情况:按图8中各施工工法的排列顺序来看,塑性区占比分别为39.2 %、39.5 %、48.5 %、48.5 %.双侧壁导坑法相较于其他三种工法,每次开挖面积较小,施工步骤复杂,造成开挖过程中对围岩的扰动较大,掌子面处的围岩应力需要较长的时间才能完成重分布,分析图8后可看出,双侧壁导坑法施工完成后,塑性区从隧道左上方地面经隧道右下方贯通至围岩底部.
图8 各工法下围岩塑性区云图Fig.8 The plastic zone of surrounding rock under different construction methods
各施工工法在完成初期喷混后,隧道支护结构的最大主应力分布云图见图9,初期支护出现最大主应力的位置都位于隧道中部,出现区域位于隧道下穿高速公路施工处.比较四种施工方法的拱顶初喷最大主应力云图的区域大小,台阶法云图面积最大且数值最大,CD法和CRD法次之,双侧壁导坑法最小.表5为初期喷混的第一主应力值,结合图9分析,台阶法的拉应力、压应力的极值相较于其他三种施工工法都是最大的,三种工法按受力效果的排序同上.比较出现最大值的台阶法和最小值的双侧壁导坑法,在最大拉应力和最大压应力上,分别增加了41.1 %和58.2 %.
图9 不同工法下初期喷混最大主应力云图Fig.9 The maximum principal stress of initial shotcrete under different construction methods
表5 初期喷射混凝土第一主应力受力值(单位:MPa)
通过分析比较四种施工工法的模拟结果后,可以发现:隧道下穿浅埋地层施工时,台阶法的结果是最不理想的,由于施工时开挖面积较大对围岩扰动最为严重,其在施工过程中在拱顶和地表沉降的数值是最大的,表明其控制效果远不如其他几种工法.双侧壁导坑法同其他三种工法相比其在控制地表和拱顶竖直方向上表现最好,在控制隧道的水平位移方面较差.从模拟结果来看,CD法与CRD法在各项数据上都逊于双侧壁导坑法,但这两种工法没有明显的差异,CRD法要略微优于CD法.
从围岩塑性区的云图面积上来看,双侧壁导坑法的塑性区占比最高.而在完成初期喷混后拱顶的受力情况方面,双侧壁导坑法最优,CRD法和CD法略为次之,台阶法施工后拱顶的最大拉应力值、最大压应力值是这些工法中最大的.
在实际施工的复杂程度方面:台阶法的施工工序较少,施工技术成熟,相较于其他工法而言比较简单,成本和效率上最好控制.CD法需要在分步施工时添加临时支撑,但较于CRD法和双侧壁导坑法还是较为方便.但在浅埋地层中开挖隧道,CD法在控制围岩变形方面稍差于这两种工法,但在施工速度和成本控制方面优于CRD法.双侧壁导坑的施工效果最好,但施工进度不快,其施工成本造价也高.
在实际施工时要考虑工期和成本,必须在满足工程本身安全可靠和质量的基础上,尽可能加快施工速度.通过上述比较,考虑到此段地铁隧道下穿既有高速公路及公路涵洞,为了保证高速公路的行车安全及正常运营,推荐此下穿段使用CRD法施工.使用CRD法施工能够很好的控制高速公路地表沉降、隧道的收敛变形速率、围岩应力每步施工时都会得到一定地释放,使隧道保持相对稳定的状态.同时,分步小断面开挖并快速初喷和添加支护,使开挖部分封闭成环,能改善隧道初支的受力情况,保障隧道施工人员的安全.
拱顶下沉观测断面所在的里程为YD1K1+450,拱顶沉降观测使用的仪器为GFS-1测微器.图10为模拟施工后得到的数据和现场监测数据,模拟施工和现场施工都选用的是CRD法施工,观测点位置在施工至断面前选择隧道拱顶上方地表,通过后选择隧道拱顶.由图中现场监测沉降曲线可知,在距离监测断面30 m的地方进行施工时,隧道拱顶地表的沉降开始逐渐增大,施工至YD1K1+450断面这段时间内,沉降值增大到10 mm左右后,沉降在隧道施工通过监测断面后基本没有增加,随着隧道接着施工,沉降值达到了13 mm左右,继续推进至15 m后,现场监测的沉降值急剧增加,推进到16 m后,拱顶沉降又再次保持稳定,在下穿施工期间,拱顶的总计沉降达到了39.8 mm,沉降变化的平均速率达到了2.487 5 mm/d.
图10 YD1K1+450处沉降变化量曲线图Fig.10 The settlement variation curve of YD1K1+450 mileage
模拟沉降曲线同现场监测曲线相比,其模拟施工至监测断面时,沉降值只有4.12 mm,再施工2 m后,沉降曲线同现场监测曲线相吻合,在15 m处,沉降值没有突然急剧增加,而是继续稳定在13 mm左右,同实际沉降曲线相背离.出现差异的原因是由于掌子面推进至此时,因天气原因暂时停止施工,复工后,隧道因受降水影响,拱顶沉降突然变大,模拟施工时未考虑降雨因素,导致监测数据与模拟数据出现巨大差异,故其沉降曲线波动很小,同实际情况相比有很大差异.
模拟沉降曲线的走势在未下雨前,与实际沉降曲线走势基本一致,说明实际施工中使用CRD法取得了预期效果.
图11为地表监测点布置图.从图12可看出,浅埋地铁隧道下穿高速公路施工,下穿隧道右线隧道上方的高速公路路面沉降比左线上方的路面大,监测点中405-3处的沉降值最大,为86.07 mm.390、395处的监测点路面沉降在27~36 mm之间,属于施工扰动的正常范围之内.395至410监测点数据显示,有8个监测点处的路面沉降值在50 mm及以上,高速公路正常运营已受到影响.出现过大沉降的原因是,下穿隧道右线隧道靠近已有涵洞,受隧道开挖扰动的影响较大.
图11 监测点布置图Fig.11 The layout of measuring points
图12 12月监测点累计沉降量Fig.12 The settlement of measuring points in December
从图13可看出,沉降最大处的405-3监测点在下穿高速公路施工完成后,其沉降值继续增加,从86.07 mm增大到87.62 mm,累计沉降1.55 mm.沉降基本保持稳定的原因有两个,一个是因为下穿施工前高速路基已受到扰动,路基沉降已基本完成,二是CRD施工后沉降得到了有效控制,在施工过程中路面没有发生过大沉降,保证了高速公路的行车安全.
图13 405-3监测点累计沉降量Fig.13 The cumulative settlement at 405-3 measuring point
经过以上的数值模拟计算和分析后,能得到下列结论:
(1)对于浅埋地铁隧道下穿既有高速公路的施工时,应分步小断面开挖,减少掌子面开挖面积,从而减轻施工对隧道周围土体的扰动.在控制地表沉降方面,CRD法仅次于双侧壁导坑法,但二者相差较小,台阶法引起的沉降最大.
(2)分析四种工法施工时隧道支护结构的受力状态,按照初喷施工后隧道拱顶的拉应力区域大小来排序为:台阶法、CD法、CRD法、双侧壁导坑法.CD法和CRD法完成初期喷混施工后,隧道拱顶处的受力状态几乎相等,台阶法中隧道支护的拉、压应力值为上述施工方案中的最大值.
(3)从施工工艺的复杂程度上看,台阶法的施工技术最为成熟且造价较低,双侧壁导坑法施工最为复杂,且施工造价高.为了尽量减少对正常运营的高速公路的通行产生影响,同时保证隧道和高速公路的安全,减少地表沉降.推荐下麦西浅埋地铁隧道下穿高速公路施工时采用CRD工法.
(4)通过对现场监测数据的分析,CRD工法在控制拱顶沉降、地表沉降方面起到了预期作用,所量测得到的沉降值在合理范围之内.