1、 专业知识分享版使命：加速中国职业化进程 联系电话：0755-86153458摘要：针对地铁联络通道冻结施工时常面临的冻胀过大且冻胀率较难测定的施工难题，依托杭州某复杂地层内越江隧道联络通道冻结施工实例，采用数值分析不同假定冻胀率对冻土体内力和主隧道变形的影响，通过冻土体内力变化制定出冻土体的临界冻胀率，同时对主隧道竖向变形进行监测，并结合主隧道计算竖向变形反馈出实际冻胀率，最后根据实际冻胀率和临界冻胀率的比较，采取相应施工措施。研究结果表明：此类地层联络通道冻结施工时，冻土体临界冻胀率为2%，可采用可观测的主隧道竖向变形监测来反馈冻土体的实际冻胀状态，当冻结和暗挖阶段主隧道竖向变形量分别小于
2、7.5 mm 和4.88 mm 时，实际冻胀率为1%2%，冻胀对冻结效果的影响尚在可控范围内，故除常规泄压孔减胀措施，无需再采用加强减胀措施。关键词：冻结法；联络通道；冻胀率；三维数值模拟；监测自人工冻结法于1883年成功应用于阿尔巴里德煤矿深 103 m 的9号井凿井工程后，在此后1个多世纪里，其被国内外不同国家广泛应用于煤矿、隧道、建筑基础和地铁等领域，而在地铁联络通道施工中，对于强度低、自稳性差、含水量高的软土地层，采用冻结加固是行之有效的方法，据不完全统计，上海地铁98%的联络通道和全部越江隧道联络通道及泵站均采用冻结法施工 1 。人工冻结土体过程是包括冰水相变、潜热释放和水分迁移的时
3、变和相变的物理力学问题，而联络通道冻结法施工时，常面临成孔偏斜超限、冻结冻胀过大 2 、开挖蠕变、融沉过大等问题，其中以由水分迁移引发的土体冻胀对冻结效果的影响最为突出，如上海复兴东路隧道某联络通道因冻胀应力过大，通道开挖后冻土体变形过大，冻结管断裂，影响冻结效果 2 。众多学者对联络通道的人工冻结冻胀问题进行了研究，岳丰田等 3 依托大连路隧道联络通道工程进行室内原状土冻结的模型试验，土层为强度较低、稳定性较差的淤泥质粘土和粉质粘土，试验得到有荷条件下冻土体达到设计厚度时冻胀率为0.76%；仇培云等 4 对大连路隧道联络通道附近主隧道竖向变形进行跟踪监测并得出隧道变形在冻结初期缓慢发展且有波
4、动，交圈过程迅速上升并达最大，且后续下降的结论；陈沙，杨文武等 5 依托崇明长江隧道联络通道工程采用 FLAC3D 对冻结过程进行三维力学分析，联络通道所在土层为灰色砂质粉土，得出随着冻土体的冻胀率增大，冻土体应力相应增加，假定冻胀率不超过6%时，冻土体计算强度值均小于设计值。国内对联络通道冻结体冻胀问题的现有研究多集中于室内试验、数值模拟和现场监测的手段，此3种研究手段各有优势和不足，室内模型实验可定性反映冻胀变化，但由于和实际状态的差异性，较难得出冻土体实际冻胀状态；现场监测虽能得到冻胀压力和主隧道的变形规律，但无法直接对冻土体的实际冻胀状态进行定量判断；而数值模拟的准确性受制于参数的取值
5、和模拟方法的选取。同时现有研究尚未给出联络通道实际冻结施工时能定量判断冻土体冻胀状态的有效方法，而由于联络通道冻结加固施工自身存在较大的风险，尤其当联络通道周边地层为复杂地层时，施工风险更大，若能定量判断冻土体冻胀状态，则能更好地控制冻结效果。以国内某越江隧道复杂地层条件下联络通道冻结加固的工程实例为背景，发挥室内试验、数值模拟和现场监测3 种研究手段各自优势，首先对联络通道所在土层分别进行室内冻土试验测定其特性和强度，为数值模拟所用参数提供依据；并通过 ansys 建立三维有限元模型，根据规范规定的联络通道周边冻土体冻胀率范围 6 ，分析不同假定冻胀率对冻土体内力和主隧道变形的影响，通过冻土
6、体内力设计值与实际值的差值（即实际安全系数）与规范规定的安全系数 7 对比结果制定出冻土体临界冻胀率；同时对冻结法施工各阶段主隧道变形进行监测，结合主隧道计算变形量反馈出冻土体的实际冻胀率，最后将临界冻胀率与实际冻胀率进行比较，实际冻胀率超出临界冻胀率时，加强减胀措施来减小冻胀，控制冻结效果。1工程背景介绍杭州某越江隧道主隧道外径6.2 m、内径为5.5 m，左右线通道中心标高均为-23.522 m，左右两侧专业知识分享版使命：加速中国职业化进程 联系电话：0755-86153458主隧道中心距离为12 m，联络通道由喇叭口、水平通道和通道下方的泵房3部分组成，联络通道结构从外到内主要由支护层
7、和结构层组成，联络通道的开挖轮廓高（外放约100 mm）4.455 m，宽4.3 m，局部喇叭口处高5.895 m，宽 4.9 m， ；泵房开挖轮廓约长 4.6 m（沿联络通道方向） ，宽 4.1 m，深 4.6 m，如图 1 所示。联络通道所在土层主要为(1) 2淤泥质粉质粘土、(9) 1a 粉质粘土、 (12) 2细砂和(12) 4圆砾层，粉砂层和圆砾层均为承压水层，水头为0.23 MPa，而圆砾层所含圆砾直径为0.22 cm，含量占30%40% ，所含卵石粒径2 6 cm 为主，含量 2030 ，圆砾层渗透系数大、强度较高、压缩性低；联络通道上部的淤泥质粉质粘土渗透系数小、强度较低、压缩
8、性高，如表1所示。由图1 和表 1可知，联络通道区域上部土层为 (8)2淤泥质粉质粘土、下部土层为(12) 2细砂和(12) 4圆砾，属于典型的上软下硬复杂地层，下部的(12) 2细砂和(12) 4圆砾层渗透系数较大，承压水头较高。2联络通道所处地层室内冻土试验由于联络通道周边土层冻结加固后所形成的冻土体为复杂的四相体系，与未冻土的性质差异较大，为较准确地给定数值模拟所需的地层参数，对联络通道所处地层进行室内冻土试验，测定冻土特性和强度。冻土体试验温度对结果影响较大，因为随温度降低，冰晶含量相对增多，冰土混合物强度提高，使冻土强度呈非线性规律提高 8 ，而本工程冻土墙设计平均温度 -10 ，由
9、于试验温度低于-10 条件下冻土的强度要高于试验温度等于-10 条件下相同冻土土样的强度，基于安全考虑，不同地层制成的冻土专业知识分享版使命：加速中国职业化进程 联系电话：0755-86153458试样均制备为-10 的冻土土样，如表2 所示。由表2 可知，在相同外部条件下， 2淤泥质粉质粘土的抗压强度与抗拉强度均小于其他土层，分别为3.56 MPa 和 2.3 MPa，基于安全考虑，将 3.56 MPa 和 2.3 MPa 分别作为冻土体抗压强度和抗拉强度设计值，而冻土体抗剪强度设计值按照经验取1.5 MPa，并参考 旁通道冻结法技术规程 6 III 类冻结壁强度检验安全系数规定：抗压不小于
10、2.0 MPa，抗拉不小于3.0 MPa，抗剪不小于2.0 MPa，当冻土体设计强度值与实际强度值比值大于检验安全系数时，冻土体处于安全状态，而其比值小于检验安全系数时，冻土体存在风险。3冻胀对冻结效果影响的数值分析3.1三维计算模型及参数选取建立三维有限元模型（如图2图3所示） ，在对称界面约束其法向位移，模型总尺寸为25 m20 m50 m，采用空间8节点六面体单元进行离散，共计划分 104 151个单元。冻土体温度设定为-10 ，不同土层的弹性模量和泊松比等计算参数采用冻土特性和强度测试结果，对于主隧道管片结构，选取 C50钢筋混凝土参数，而对于联络通道结构，按联络通道设计取值，如表3所
11、示。河床以上静水压力以附加面荷载施加，承压水头也以附加面荷载作用在承压含水层顶板和底板界面处。专业知识分享版使命：加速中国职业化进程 联系电话：0755-86153458为分析冻胀对冻结效果的影响，采用反应冻胀状态的冻胀率指标，借鉴建筑地基基础设计规范4 关于圆砾、砂土及粘土冻土冻胀率范围的规定，数值计算冻胀率上限考虑为6%，且由于冻胀率小于1%为不冻胀范围，该范围对冻结效果无影响，故冻胀率下限考虑为1%。3.2计算结果及分析由于数值计算冻胀率范围为1% 6%，故分别对1%，2% ，3%，4%，6%五种不同的冻胀率对主隧道竖向变形和冻土体内力的影响分别进行计算，冻结阶段4%冻胀率对主隧道和冻土
12、体的影响如图4所示。在联络通道冻结施工阶段，由图4（a）4 （c）知，冻土体待开挖的水平通道侧墙与主隧道交界位置出现拉应力最大值，达到0.264 MPa，呈明显应力集中，分布范围不大，冻土体大部分区域拉应力小于0.42 MPa；而冻土体待开挖的水平通道拱底和主隧道交界位置出现压应力最大值，达到2.27 MPa，呈明显应力集中，分布范围较小，冻土体大部分区域压应力小于2.05 MPa；冻土体待开挖的水平通道拱底和主隧道交界位置出现剪应力最大值，达到1.12 MPa，呈明显应力集中，分布范围较小，冻土体大部分区域剪应力小于1.08 MPa。而由图4（d）可知，4%冻胀率下主隧道竖向变形最大位置为靠
13、近联络通道位置且与冻土帷幕交界的主隧道侧壁，将竖直向上假定为正方向，隧道竖向变形最大正值为+22.1 mm，主隧道拱底竖向变形量为+12.8 mm，隧道整体有上浮趋势。分别对1% ，2% ，3%，4%，6%5 种不同的冻胀率对主隧道竖向变形和冻土体内力的影响进行计算，暗挖阶段4% 冻胀率对主隧道和冻土体的影响如图5 所示。专业知识分享版使命：加速中国职业化进程 联系电话：0755-86153458在联络通道暗挖施工阶段，由图5（a）（c）知，泵房外底部角点区域的冻土体出现拉应力最大值，达到1.53 MPa，呈明显应力集中，分布范围不大，冻土体大部分区域拉应力小于 1.16 MPa；而泵房外底部
14、底边区域的冻土体出现压应力最大值，达到2.4 MPa，呈明显应力集中，分布范围较小，冻土体大部分区域压应力小于2.18 MPa；泵房外底部角点区域的冻土体也出现剪应力最大值，达到 1.9 MPa，呈明显应力集中，分布范围不大，冻土体大部分区域剪应力小于1.53 MPa，剪应力最大值已经超出冻土体设计抗剪强度。而由图5（d）可知，暗挖施工后 4%冻胀率下主隧道竖向变形最大位置发生转移，主隧道竖向变形有正有负，竖向变形最大正值在靠近联络通道的主隧道拱底处，最大正值为+7.7 mm。冻结和暗挖不同阶段不同假定冻胀率对主隧道变形和冻土体应力的影响如表4所示。冻结及暗挖阶段不同冻胀率下强度安全系数变化如
15、图6所示。由表4 、图6可知，随着假定冻胀率增加，冻土体的应力值均相应增加，且联络通道附近主隧道拱底竖向变形量也随之增大，在3% 冻胀率下，冻结和暗挖不同阶段冻土体抗拉和抗压安全系数均小于规范规定的安全检验安全系数，但冻土体暗挖阶段最大剪应力为0.97 MPa，暗挖阶段抗剪安全系数为1.55，小于相应的强度检验安全系数，而在2%冻胀率下，不同阶段的冻土体强度安全系数均大于相应的强度检验安全系数，此时暗挖阶段抗剪安全系数为2.11，略大于抗剪安全系数2 ，故2%冻胀率为冻土体临界冻胀率。专业知识分享版使命：加速中国职业化进程 联系电话：0755-861534584主隧道竖向变形监测对联络通道实际
16、施工过程中冻结管成孔、冻结及暗挖不同阶段主隧道竖向变形进行监测，主隧道竖向变形监测点沿着联络通道中心线近似对称分布，相邻监测点相隔2.44.8 m，且随着离中心线的距离增加，相邻监测点间距增大，并以联络通道中心线为原点，以此原点沿着掘进方向为正，背离掘进方向为负值，监测点具体布置如图7所示。冻结和暗挖阶段隧道拱底不同测点竖向变形增量如图8所示。为分析冻土体冻胀对主隧道的影响，需剔除成孔施工阶段主隧道的已有主隧道竖向变形量，故两阶段主隧道拱底测点的变形增量分别为冻结阶段、暗挖阶段主隧道竖向变形量与成孔阶段主隧道竖向变形量的差值。由图8 可知，冻结阶段和暗挖阶段主隧道拱底测点的竖向变形增量均呈现“
17、距联络通道越近，量值越大” 的规律。冻结阶段联络通道附近主隧道测点拱底竖向变形最大增量为+7.5 mm，暗挖阶段联络通道附近主隧道测点拱底竖向变形最大增量为+4.88 mm，且两阶段竖向变形最大值均出现在同一测点，受暗挖阶段联络通道及泵房的开挖卸载作用，联络通道附近测点竖向位移量要小于冻结阶段的变形量，最大差专业知识分享版使命：加速中国职业化进程 联系电话：0755-86153458值为2.62 mm。通过冻结和暗挖阶段联络通道附近主隧道拱底竖向变形量+7.5 mm、+4.88 mm 与表3计算得出的不同假定冻胀率下主隧道的变形量的对比，可得竖向变形实测值小于2%假定冻胀率下主隧道的变形量，且
18、大于1%假定冻胀率下主隧道的变形量值，故冻土体的实际冻胀率范围为1%2%，小于临界冻胀率2%。5结语本研究以冻土强度及特性室内试验、数值模拟和现场实测为研究手段，将室内试验所得冻土参数应用于数值计算分析，得出冻土体的临界冻胀率，并根据现场实测主隧道拱底竖向变形量与数值计算结论进行对比，得出冻土体的实际冻胀率，最后通过临界冻胀率与实测冻胀率的对比分析，确定是否在常规减胀措的基础上，进一步采用加强减胀措施，研究表明：1）在此类地层进行联络通道冻结施工时，冻土体的临界冻胀率为 2%。2）冻结阶段和暗挖阶段主隧道拱底测点的实测竖向变形增量均呈现“距联络通道越近，主隧道拱底竖向变形越大”的规律，冻土体的实际冻胀率在1%2% 之间，小于临界冻胀率，故除常规泄压孔减胀措施，无需再采用加强减胀措施。3）在类似地层进行联络通道冻结施工时，采用可观测的主隧道竖向变形监测，反馈冻土体的冻胀状态，若冻结阶段主隧道竖向变形量小于+7.5 mm，且暗挖阶段主隧道竖向变形量小于+4.88 mm，则冻胀对冻结效果的影响尚在可控范围内。
