盾构隧道局部不同部位漏水对隧道和地表沉降的影响.docx
- 文档编号:29356856
- 上传时间:2023-07-22
- 格式:DOCX
- 页数:9
- 大小:24.23KB
盾构隧道局部不同部位漏水对隧道和地表沉降的影响.docx
《盾构隧道局部不同部位漏水对隧道和地表沉降的影响.docx》由会员分享,可在线阅读,更多相关《盾构隧道局部不同部位漏水对隧道和地表沉降的影响.docx(9页珍藏版)》请在冰豆网上搜索。
盾构隧道局部不同部位漏水对隧道和地表沉降的影响
盾构隧道局部不同部位漏水对隧道和地表沉降的影响
(1.天津大学滨海土木工程结构与安全教育部重点实验室,天津300072;2.天津大学水利工程仿真与安全国家重点实验室,天津300072;3.天津市市政工程设计研究院,天津300051)摘要:
针对富水地区中隧道管片局部在隧道不同区域发生漏水灾害情况,建立三维流固耦合模型,对比分析了不同位置发生漏水以及在不同漏水量下地表沉降、隧道变形规律。
研究表明隧道漏水位置越靠近隧道底部时,地表沉降最大值越小,沉降槽呈现宽而浅,漏水量越大,沉降越明显;隧道局部漏水会使隧道纵向沉降不均匀,使得隧道环缝张开,漏水位置越靠近隧道底部,环缝张开量越大,张开量与流量之间成正比,从而验证了隧道漏水与隧道变形耦合发展的机理。
关键词:
隧道工程;漏水位置;灾害;沉降;隧道变形0引言在饱和软土地区,盾构隧道不均匀沉降引起隧道发生纵向弯曲,导致隧道环缝张开,引发的隧道渗水的情况十分常见。
隧道渗水引起隧道进一步变形,影响隧道运营安全,甚至引发更为严重的漏水漏砂事故,国内的金山引水隧道和前苏联的圣彼得堡1号线隧道都曾因隧道漏水漏砂而被迫终止运营[1-2]。
大多数学者[3-5]采用二维模型进行数值模拟计算,但隧道发生局部渗水主要集中在环缝、纵缝、螺栓孔等位置,隧道漏水问题实质是一个复杂的三维问题,二维模型难以反映隧道的纵向变形情况,在隧道纵向受力上,忽略了隧道自身对该环隧道管片的约束作用;在水的补给上,二维模型中水的补给全部来自于漏点横向两侧,忽略了来自隧道纵向上的补给。
针对二维模型的不足,有的学者[6]采用三维模型进行了渗水模拟,从二维到三维,学者们针对隧道的长期渗水行为研究成果斐然,但对于隧道局部渗水流量的增大导致隧道纵向变形,从而引起灾害的情况关注较少。
此外,隧道发生漏水灾害后在短期即可对隧道和土体造成较大影响,目前缺乏相关研究。
笔者建立了三维模型对隧道局部漏水问题进行流固耦合分析,分析了隧道不同漏水位置、流量大小对地表沉降及隧道变形的影响。
本文仅研究盾构隧道局部漏水(土颗粒没有漏出)这种情况下的影响,更为极端的土体流失的情况将另撰文研究。
1盾构隧道局部渗漏的数值模拟1.1三维数值计算模型天津市地铁2号线建国道—天津站区间隧道采用土压平衡盾构施工,隧道顶部的平均埋深为10m,盾构管片厚度350mm,管片宽度1.2m,管片外径6.2m,混凝土管片的弹性模量为34.5GPa,泊松比为0.16,本文以该工程为背景建立数值模型。
采用FLAC3D有限差分软件,建立可考虑流固耦合的三维模型,如图1所示。
模型计算域横向(X轴向)长度为120m,纵向(Y轴向)长度为120m,深度(Z轴向)为53m,模型坐标原点取在模型地表中心处,隧道中心轴线沿Y向,位于X=0m,Z=-13m地表之下13m处,经过多次试算,该计算域尺寸基本满足边界效应的影响要求。
同时,模型网格在隧道附近进行了加密。
图1有限差分计算模型
Fig.1Modeloffinitedifferencemethodmodel本文选择M-C模型进行隧道漏水分析和计算,J.H.Shin等[3-5]文章的分析结果也验证了M-C模型对于隧道渗水引发土体变形计算问题的适应性,故本文采用M-C模型进行隧道的漏水分析和计算。
按照所选定的本构模型,数值计算所需要的土层物理力学参数如表1所示。
表1土层参数表
Table1Parametersofsoillayers土层编号土层名称土层埋深/m密度/(103kg·m-3)弹性模量E/MPa孔隙率渗透系数/(10-4cm·s-1)黏聚力/kPa内摩擦角/(°)①1杂填土21.394.960.485.5560010①2粉质黏土41.419.200.480.05781510②2粉土221.6717.400.400.48001529③1粉质黏土361.5716.800.420.02403020④1褐黄色粉质黏土531.3817.000.330.002531191.2隧道刚度选取盾构隧道管片接头较难模拟,很多学者采用刚度等效折减的方法模拟接头的作用[7-8]。
本文同样引入横向刚度有效率和纵向刚度有效率来考虑隧道接头对其刚度的影响,即隧道刚度各向异性[7]。
根据各向同性材料的弹性矩阵对称性可知:
在隧道横截面上,弹性模量通过横截面上的横向刚度有效率ηxy进行折减,而在隧道纵向上,弹性模量通过纵向刚度有效率ηz进行折减,即:
横断面上的剪切模量与横断面弹性模量有关,剪切模量Gxy与Ex和Ey同比例折减;而纵向上的剪切模量与弹性模量相互独立,根据何川等[8]的研究,可采用纵向上的抗剪切刚度折减系数进行折减:
式中:
下标x,y和z为各向同性材料的3个主方向;E为原混凝土的弹性模量;G为原混凝土的剪切模量;kq为纵向接头螺栓的剪切刚度;A为管片环横截面面积;ls为单环管片衬砌的幅宽。
根据钟小春等人[9]的研究,横向刚度有效率的取值在0.4~0.8左右,黄宏伟等[10]也通过模型试验证实隧道横向刚度有效率在0.7左右,本文中模型横向刚度有效率取值为0.7。
隧道衬砌层纵向变形主要由结构弯曲变形引起,钟小春等[11]研究表明纵向刚度有效率取值和隧道横向刚度有效率、环缝处的螺栓数量、管片环宽度,纵向螺栓预紧力等因素有关。
本文横向刚度有效率取值为0.7,天津地铁二号线隧道结构参数如表2所示,根据钟小春的研究成果,得到纵向刚度有效率在0.1左右,故本文中模型纵向刚度有效率取值为0.1。
模型中隧道参数见表3。
表2盾构隧道的主要结构参数
Table2Structuralparametersofshieldsegment外径D/m6.25.51.234.51.48304001720670内径d/m环宽ls/m弹性模量E/GPa剪切模量G/GPa螺栓直径d/mm长度l/mm数量/个螺栓弹性模量El/GPa螺栓预应力Nl/MPa表3衬砌管片参数
Table3Parametersoftunnelsegment密度/(g·m-3)vEx/GPaEy/GPaEz/GPaGxy/GPaGyz/GPaGzx/GPa2.50.224.224.23.51.040.170.171.3隧道漏水的模拟及工况潘海泽等[12]根据衬砌渗漏水量的大小及渗漏状态将隧道渗漏水灾害分为渗润、滴水、流水、喷水、淋瀑5个等级,在其调查的95座隧道之中有7%的隧道存在缝漏喷水现象,根据经验,一处喷水的漏水量可达1m3/d,隧道在某一局部区域产生5处喷水即可致灾,即产生5m3/d漏水量,以此漏水量为基准建立隧道局部漏水的灾害环境。
本文采用Shin等[3]提出的等效衬砌渗透系数的方法来实现漏水,通过改变衬砌的渗透系数来模拟漏水。
衬砌渗透系数是其劣化、裂缝等的综合反映。
通过试算的方法,不断调整衬砌与土体渗透系数比,最终确定不同工况下漏水位置处衬砌的渗透系数取值。
隧道发生漏水的位置主要集中在隧道环缝、纵缝以及螺栓接头等地方[4],为了反映隧道局部不同位置的漏水情况,将隧道漏水位置选在隧道模型中部一环管片(Y=-0.6m到Y=0.6m)的4个区域,如图1和图2所示。
图2渗漏隧道分区
Fig.2Differentzonesofleakagetunnel针对在这4个区域发生相同的漏水量以及在同一区域发生不同的漏水量建立了7种工况,如表4所示。
表4计算工况
Table4Calculatedworkingconditions工况漏水区域衬砌与土体渗透系数比漏水量/(m3·d-1)1Ⅰ区0.215.02Ⅱ区0.175.03Ⅲ区0.145.04Ⅳ区0.135.05Ⅰ区0.407.56Ⅰ区0.5910.07Ⅰ区0.8112.51.4计算分析步骤及边界条件设置本文主要研究隧道建成之后运营期间发生漏水所产生的影响,在前期首先模拟盾构隧道的开挖,主要的分析步骤如下:
1)初始地应力平衡;2)模拟隧道开挖;3)设置漏水区域,通过调整衬砌渗透系数确定漏水量,进行漏水模拟。
模型位移边界条件为:
地表为自由边界;模型土体的侧向边界和底面边界约束各自法向位移,模型隧道始发和接收端也约束其转动位移。
模型中地下水位位于地表以下1.5m,并在模型侧面施加了孔压边界条件,模型侧面水位以下施加梯度为9.8kPa/m的静水压力以模拟水源补给;模型底面和隧道不透水部分设置为不透水边界;隧道发生漏水区域在漏水期间设置为透水边界,并将该部分隧道内侧孔压设置为0。
2隧道局部漏水的影响分析2.1隧道局部漏水对地表沉降的影响图3是隧道局部不同位置发生漏水时地表最大沉降随时间的变化情况,图4为200d时隧道局部不同位置漏水时横断面处地表沉降曲线。
图3不同漏水位置地表沉降最大值随时间变化曲线
Fig.3Variationsofmaximumsurfacesettlementwithtimeatdifferentleakagepositions图4不同漏水位置地表沉降槽
Fig.4Surfacesettlementtroughsatdifferentleakagepositions综合这两张图,可以得出以下规律:
1)隧道局部发生漏水之后,前72h地表迅速发生沉降,之后地表沉降缓慢发展,漏水灾害的治理应该在72h内展开[5]。
2)漏水量相同时,漏水位置越靠近隧道顶部,地表最大沉降量越大,隧道局部Ⅰ区漏水相比于Ⅳ区漏水时地表最大沉降提高16%,分析其原因,在三维模型之中考虑了隧道纵向的约束作用,隧道不能轻易随着土动,当漏水位置降低,使得更多的土体固结发生在隧道下侧,这部分土体的变形需绕过隧道传递到地表,隧道阻碍土体变形的向上传递。
而二维模型无法考虑隧道的纵向约束,往往隧道随着土体运动,从而得出相反结论(漏点越低,沉降值越大[10-11]),与隧道实际工作状态不符。
这种隧道对土体变形的扩散效应最终使得隧道局部漏水位置越低,地表沉降最大值反而越小。
3)地表沉降槽形状接近peck曲线,隧道位置靠近隧道底部时,地表沉降槽呈现宽而浅,而当漏水位置靠近隧道顶部时,地表沉降槽呈现窄而深,如果上方存在建筑物或者刚性路面,更可能因为不均匀沉降而发生开裂。
4)由于漏水位置的偏离,会使得沉降槽最低点偏离隧道中心竖轴,由于隧道的纵向刚度的影响,隧道下方的土体固结变形需绕过隧道传递到地表,所以工况4和工况3相比于工况1和工况2沉降槽最低点偏离隧道中心竖轴更远,这也进一步佐证了第2)条中土体变形扩散的推论。
图5是不同漏水量引起的地表最大沉降随时间的变化情况。
图5不同漏水流量下地表沉降最大值随时间变化曲线
Fig.5Variationsofmaximumsurfacesettlementwithtimeatdifferentleakageflowrates图6、图7为隧道局部在不同漏水量下第200天时横断面处地表沉降曲线和最大沉降值。
图6不同漏水量地表沉降槽
Fig.6Surfacesettlementtroughsatdifferentleakageflowrates图7不同漏水量下地表最大沉降值
Fig.7Maximumsurfacesettlementatdifferentleakageflowrates综合图5~图7,可以得出以下规律:
1)漏水量增大,会使得地表沉降发展更迅速,同时沉降达到稳定所需的时间越长。
2)漏水量增大,地表沉降越大,漏水影响的范围也越广,工况6相比于工况1漏水量提高1倍,地表沉降最大值提高2.4,接近等比例,吴怀娜等[13]也得到类似的规律。
相比隧道局部不同漏水位置对地表沉降的影响,漏水量差异造成的影响更加明显。
2.2隧道局部漏水对隧道的影响图8和图9分别为隧道局部不同位置漏水和不同漏水量下隧道最大沉降值随时间变化曲线,结合图3和图5,可以看出隧道沉降的发展相比于地表沉降更平缓一些,漏水位置越靠下,隧道下部固结的土体越多,使得隧道沉降越大,工况4相比于工况1,隧道沉降最大值增大112%;当流量增大时,隧道下部土体固结值增大,使得隧道最大沉降量也增大。
图8不同漏水位置隧道沉降最大值随时间变化曲线
Fig.8Variationsofmaximumtunnelsettlementwithtimeatdifferentleakagepositions图9不同漏水流量下漏水隧道沉降值随时间变化曲线
Fig.9Variationsofmaximumtunnelsettlementwithtimeatdifferentleakageflowrates不同工况下隧道与地表之间沉降最大值之间的比值(沉降比)如表5所示,各工况的沉降比位于0.33~80之间,漏水位置越低,沉降比越大,这与刘印[4]二维渗漏水模型的计算结果一致。
随着漏水量增大,沉降比增加并不明显。
表5隧道与地表之间沉降最大值的比值
Table5Ratioofmaximumtunnelsettlementandmaximumsurfacesettlement工况隧道沉降/mm地表沉降/mm比例1-2.09-6.440.332-2.78-6.260.443-3.63-5.860.624-4.43-5.560.805-3.99-11.800.346-5.37-15.430.357-6.46-18.100.38隧道发生漏水后,引起隧道周围土体的固结,最终引发隧道的下沉,图10和图11分别为隧道局部不同位置漏水和不同漏水量下隧道中心轴线的沉降曲线,隧道局部漏水会使得隧道纵向沉降不均匀,漏水位置越靠下,隧道沉降的范围越广;漏水量增大,隧道沉降的范围越广。
图10不同位置漏水隧道纵向沉降曲线
Fig.10Tunnellongitudinalsettlementcurveatdifferentleakagepositions图11不同流量下隧道纵向沉降曲线
Fig.11Tunnellongitudinalsettlementcurveatdifferentleakageflowrates隧道的不均匀沉降不仅影响列车的运营,这种隧道曲率上的变化更会使得环缝进一步张开[13],刘印等[14]也通过隧道漏水机理分析验证了隧道曲率对于隧道环缝张开量的显著影响。
隧道漏水引起的隧道纵向曲率半径的变化见表6,在笔者计算的几种工况下,最小曲率半径为44.5km,继续发展可能超过《地铁隧道保护条例》规定由外界因素引起的隧道曲率半径不小于15km的限制限值。
表6不同工况下隧道曲率半径和环缝张开量
Table6Curvatureradiusandjointopenningoftunnelatdifferentworkingconditions工况漏水横断面处隧道曲率半径/km环缝张开量/mm工况1116.30.05工况268.10.09工况356.80.11工况447.70.13工况568.90.09工况650.10.13工况744.50.14假设此时隧道环缝螺栓还处于弹性状态,将表2中的隧道结构参数以及表6中漏水位置处隧道的曲率半径代入式(4)[21]中:
式中:
δj为与相对应的环缝接头螺栓的最大张开量;R为纵向曲率半径;ls为环宽;x和φ分别为中性轴的位置和角度。
得到不同工况下环缝接头螺栓的最大张开量如表6所示。
当隧道局部发生漏水后,隧道沉降使得隧道环缝最大张开量增大,当漏水位置越靠近隧道底部时,环缝张开量越大,当环缝张开量越大,漏水量会增加,对比工况1和工况5~7可知,漏水量增大会使得接头螺栓张开量越大,进一步恶化漏水环境,从而验证了隧道漏水与隧道变形耦合发展的机理。
3结语本文通过建立三维模型,分析了在灾害环境下隧道局部一环管片的不同位置发生漏水以及不同漏水量下地表沉降、隧道沉降的规律,主要结论如下:
1)由于隧道纵向刚度产生的阻碍作用,使得隧道局部漏水位置越低,地表沉降最大值越小。
地表最大沉降值与漏水量成线性正比的关系,漏水量增大,地表受影响的范围也扩大,同时沉降值达到稳定所需要的时间也越久;灾害发生后72h内地表沉降迅速发展,这段时间进行治理最为有效。
2)隧道沉降达到稳定相比地表沉降达到稳定更缓慢,且漏水位置越靠近隧道底部,漏水量越大,沉降达到稳定时间越长;漏水位置越靠近隧道底部,隧道沉降值越大,沉降值与漏水量之间成正比。
3)隧道局部漏水使得隧道纵向发生不均匀沉降,漏水位置越靠近隧道底部,环缝张开量越大,造成漏水量增大,漏水量增大,使得环缝张开量更大,隧道漏水与隧道变形之间恶化耦合发展。
参考文献:
[1]LIAOSM,PENGFL,SHENSL.Analysisofshearingeffectontunnelinducedbyloadtransferalonglongitudinaldirection[J].TunnelingandUndergroundSpaceTechnology,2008,23(4):
241-430.[2]白云,肖晓春,胡向东.国内外重大地下工程事故与修复技术[M].北京:
中国建筑工业出版社,2012.BAIYun,XIAOXiao-chun,HUXiang-dong.Domesticandinternationalimportantengineeringaccidentsandrepairtechnology[M].Beijing:
ChinaArchitecture&BuildingPress,2012.[3]SHINJH,ADDENBROOKETI,POTTSDM.Anumericalstudyoftheeffectofgroundwatermovementonlong-termtunnelbehaviour[J].Géotechnique,2002,52(6):
391-403.[4]刘印,张冬梅,黄宏伟.盾构隧道局部长期渗水对隧道变形及地表沉降的影响分析[J].岩土力学,2013,34
(1):
290-298.LIUYin,ZHANGDong-mei,HUANGHong-wei.Influenceoflong-termpartialdrainageofshieldtunnelontunneldeformationandsurfacesettlement[J].RockandSoilMechanics,2013,34
(1):
290-298.[5]郑刚,戴轩.灾害环境下隧道不同部位漏水对于周围土体及平行隧道的影响研究[J].岩石力学与工程学报,2015,34(S1):
3196-3207.ZHENGGang,DAIXuan.Influenceofdifferentleakagepositionsoftunnelonsurroundingsoilsandparalleltunnelunderdisasterenvironment[J].ChineseJournalofRockMechanicsandEngineering,2015,34(S1):
3196-3207.[6]WONGSAROJJ,SOGAK,MAIRRJ.Modellingoflong-termgroundresponsetotunnellingunderStJames′sPark,London[J].Géotechnique,2007,57
(1):
75-90.[7]曾东洋.盾构隧道衬砌结构力学行为及施工对环境的影响研究[D].成都:
西南交通大学,2005.ZENGDong-yangResearchonthemechanicalbehaviorofshieldtunnelliningstructureandinfluenceontheenvironmentinconstruction[D].Chengdu:
SouthwestJiaotongUniversity,2005.[8]何川,苏宗贤,曾东洋.盾构隧道施工对已建平行隧道变形和附加内力的影响研究[J].岩石力学与工程学报,2007,26(10):
2063-2069.HEChuan,SUZong-xian,ZENGDong-yang.Researchoninfluenceofshieldtunnelconstructionondeformationandsecondaryinnerforceofconstructedparalleltunnel[J].ChineseJournalofRockMechanicsandEngineering,2007,26(10):
2063-2069.[9]钟小春,朱伟,季亚平,等.盾构衬砌管片环弯曲等效刚度的一种确定方法[J].地质与勘探,2003,39(S):
185-189.ZHONGXiao-chun,ZHUWei,JIYa-ping,etal.Themethodtoconfirmtheeffectivebendingstiffnessofshield-driventunnellining[J].GeologyandProspecting,2003,39(S):
185-189.[10]黄宏伟,徐凌,严佳梁,等.盾构隧道横向刚度有效率研究[J].岩土工程学报,2006,28
(1):
11-18.HUANGHong-wei,XULing,YANJia-liang,etal.Studyontransverseeffectiverigidityratioofshieldtunnels[J].ChineseJournalofGeotechnicalEngineering,2006,28
(1):
11-18.[11]钟小春,张金荣,秦建设,等.盾构隧道纵向等效弯曲刚度的简化计算模型及影响因素分析[J].岩土力学,2011,32
(1):
132-136.ZHONGXiao-chun,ZHANGJin-rong,QINJian-she,etal.Simplifiedcalculationmodelforlongitudinalequivalentbendingstiffnessofshieldtunnelanditsinfluencefactors′analysis[J].RockandSoilMechanics,2011,32
(1):
132-136.[12]潘海泽,黄涛,杨海静,等.运营隧道渗漏水灾害分类和等级评定方法[J].干旱区地理,2009,32
(1):
145-161.PANHai-ze,HUANGTao,YANGHai-jing,eta1.Classificationandgradingassessmentmethodofleakagedisastersinrunningtunnel[J].AridLandGeography,2009,32
(1):
145-161.[13]郑永来,韩文星,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 盾构 隧道 局部 不同 部位 漏水 地表 沉降 影响