1、富水砂性地层盾构开挖流摘要：运用 FLAC3D 三维非线性有限差分法计算程序，考虑了隧道开挖、管片拼装、盾尾注浆、浆液固结等施工步骤，对某地区盾构通过富水砂层地段进行了流-固耦合与非耦合分析，并与实测结果进行了对比，得出了此地层下盾构推进过程中地表沉降变化规律，揭示了盾构掘进施工对周围环境的影响特征。计算方法和研究结果可供类似盾构工程参考。关键词：盾构；富水砂层；流固耦合；沉降监测 中图分类号：P539 文献标识码： A Research on Shield Excavation Considering Fluid Solid Coupling Analysis at Watered Sand
2、y Stratum Zhang Peng, Yao Shengfeng （Anhui urban construction detection technology Co,LTD, He Fei,230088） Abstract: Using three-dimensional nonlinear finite difference program FLAC3D, consider the tunnel excavation, the segment assembling, tail shield grouting, construction steps, such as slurry con
3、solidation of shield in one region with rich water sand section coupled fluid-solid interaction was carried out with the comparison and analysis, and compared with the measured results, obtained the change rule of surface settlement of shield of the process, to reveal the influence of shield tunneli
4、ng construction on the surrounding environment. Calculation method and the results can be used for reference in the similar shield engineering. Key words: shield ; Watered sandy stratum ; fluid solid coupling subsidence monitoring 引言 随着我国轨道交通工程建设的规模和难度越来越大，工程事故时有发生，大量的城市地下工程事故一般都与地下水有关，其主要原因是隧道开挖后，地
5、下渗流场的改变导致地应力场的调整，引起隧道围岩变形或者失稳1。研究表明2，在地下工程施工过程中，存在地下水的渗流场和地应力场耦合作用问题，主要表现为岩土的变形引起岩土渗透性能的改变，导致流体孔隙压力发生改变；另外，流体孔隙压力的改变使得岩土的应力状态发生变化，同时改变了岩土的物理力学性质。本文将流固耦合思想引入到隧道开挖地应力场及渗流场分析中，利用有限差分软件 FLAC3D 对盾构在富水砂层中掘进引起的周围土体位移场进行对比研究，并与监测结果进行对比，揭示盾构掘进对周围环境的影响特征。为隧道工程的方案设计提供了一定帮助，从而为后续施工安全提供技术指导。 工程及地质概况 本文研究对象为某城市轨道
6、交通一盾构区间，起于艾溪湖东站，止于太子殿站，为单圆地下铁道，选用德国海瑞克土压平衡盾构机进行隧道掘进，开挖荒径 6.28m，采用 C50 钢筋混凝土管片，宽 1.2m，厚0.30m。本文选取该区间某一段运用 FLAC3D 有限差分程序进行本段沉降的数值模拟计算。 场地地层由人工填土（Q4ml） 、第四系上更新统冲积层（Q3al） ，下部为第三系新余群（Exn）基岩。隧道主要穿越地层为4 粗砂层和5砾砂层；隧道基地为圆砾和强风化粉砂质泥岩；隧道覆土厚度：8.973m21.543m。场地地下水类型可分为上层滞水、松散岩类孔隙水、红色碎屑岩类裂隙溶隙水三种类型。上层滞水主要赋存于浅部素填土层之中，
7、无上覆隔水层，下部粉质粘土层为其隔水底板。松散岩类孔隙水主要赋存于上更新统(Q3al)冲积砂砾卵石层中，为潜水，地下水位埋深16.1518.90m。红色碎屑岩类裂隙溶隙水主要赋存于第三系新余群含钙粉砂岩与钙质泥岩层段。 计算模型 模型建立 本文计算模拟太子殿站至艾溪湖东站上行线地铁盾构掘进过程，方向自东向西。研究段隧道在上行线 2 号联络通道处（中心里程 SK25+600） ，隧道直线通过，纵坡 27，研究段长 60m，考虑到方便建模，将研究段隧道按直线段处理。该段隧道平均埋深 17 米。 2）模型建立与单元剖分 在隧道轴线是平面应变问题，根据圣维南原理，将模型尺寸取为地铁隧道开挖洞径的 35
8、 倍最合理。考虑到盾构隧道开挖土工环境的对称性，取隧道及周围地层的一半进行建模。模型的尺寸为：30m（x 方向)60m(y 方向)36m(z 方向)。建立三维模型如图 1 所示，共 13923 个节点，12300 个单元。 图 1 三维计算模型 本构模型、结构单元的选取及材料性能的模拟 假设条件 地表面和各土层均呈均质水平层状分布； 用弹性实体单元模拟作为隧道支护结构的衬砌管片； 计算时在盾构作业开挖面处施加一定的表面力，以模拟土仓压力；在等代层外侧土体施加表面力来模拟注浆压力，土仓压力、注浆压力大小按照实际施工参数确定； 盾构推进一步的长度为 1.2m； 将土层视为各向同性、连续的弹塑性材料
9、，材料塑性屈服准则采用Mohr-Coulomb 弹塑性屈服准则。 计算流-固耦合效应时，采用各项同性渗流模型，将土层视为多孔介质，流体在孔隙介质中流动服从达西定律，同时满足 Biot 方程。 计算参数 管片模拟：C50 钢筋混凝土管片，宽 1.2m，厚 0.30m， ，弹性模量为34.5GPa，泊松比为 0.2，密度为 2500kg/m3。 注浆层模拟：依张云等3提出的“等代层”概念。等代层厚度约为14cm(厚度为盾构机开挖半径和隧道衬砌外层半径之差)，注浆材料的强度会随时间而增强，为了模拟注浆材料不同阶段的力学性质，根据实际浆液凝固过程，分为以下 2 种：液体注浆材料视为低刚度材料，弹模取
10、5MPa，重度为 16.5kN/m3，泊松比为 0.2；短期固化注浆材料取 2h固化后浆体的弹模值，即 48MPa，重度为 19.5kN/m3，泊松比为 0.25。 土层参数依据本工程的勘察报告，在富水砂层中模拟盾构开挖计算过程中所采用的有关计算参数如表 1。 边界条件 力学边界采用齐次边界条件，对于渗流边界，除地下水位处设定为零孔压边界外，其余边界均为不透水边界。 隧道开挖后渗流边界条件：未拼装管片时，认为开挖边界为自由透水边界；管片拼装后，管片和注浆层视为渗透率极低的各项同性渗流模型。 表 1 各土岩层物理力学参数 土层名称 含水量 % 密度 /(kg/m3) 粘聚力 c/kpa 摩擦角
11、/() 弹性 模量 E/Mpa 剪切模量 E/Mpa 孔隙率 渗透率 k/(m2Pa-1s-1) 模拟土 层深度/m 素填土 25 1850 5 10 2.80 0.91 0.44 5.10E-09 2 粉质粘土 23.5 1930 55 22 25.62 7.81 0.43 3.62E-12 4 细砂 24.2 1950 0 30 43.08 8.59 0.43 3.67E-09 3 中砂 23.2 1960 0 32 43.38 9.42 0.40 2.95E-08 3 粗砂 20 2020 0 35 50.31 10.09 0.36 1.24E-07 5 砾砂 21.4 2000 0 3
12、6 42.42 10.94 0.39 1.24E-07 6 圆砾 23 2050 0 38 53.07 11.70 0.39 1.24E-07 5 强风化粉 砂质泥岩 23 2200 55 23 60.00 18.99 0.18 3.42E-10 8 计算结果对比分析 盾构沿着 y 轴方向掘进，文中分析工况为模拟盾构从 y=0m 掘进至y=54m，分析中，实现了盾构推进、盾尾间隙、土体应力释放、管片拼装、盾尾注浆、浆液固结等过程模拟。盾构推进 48m 地表沉降云图如图 2 所示。 图 2 地表沉降云图（盾构推进 48m) 纵向沉降位移分析 图 3 为盾构机掘进至 36m 位置时，隧道地表轴线沉
13、降曲线。从图中可以得出：由于盾构机掘进，扰动了周围的土层，隧道周围的土体向隧道内部的位移。隧道上方的土体产生向下的位移，即出现沉降。 图 3 地表沿隧道轴向沉降曲线 由图 3 可见，在开挖面前方 20m 与后方 10m 范围内，隧道中轴线地表纵向沉降的监测值与计算值基本吻合；在盾构开挖面到达前 20m 处，地面开始产生沉降，地面沉降主要发生在盾构通过以及盾尾脱出后1015m 范围内。 开挖面前方土体的沉降值最大不超过 6mm，且离开挖面越远位移量越小，距开挖面 20m 以外的土体几乎不产生位移，开挖面后方，地表沉降表现为距开挖面越远沉降量越大。 横向沉降位移分析 模型地表横向监测断面，取在纵向
14、 Y=12m 的位置。盾构与监测断面的几何位置 Y 表示盾构机刀盘到监测断面的水平距离，距离 Y 为正值表示盾构机没有到达监测断面，Y 为负值表示盾构机刀盘已穿过监测断面。 图 4 地表监测断面沉降槽变化曲线(计算值) 图 3 和图 4 表明：对于盾构隧道施工，隧道开挖诱发地层变形依然存在，且呈现明显的空间效应；沿隧道轴线方向的地表纵向沉降曲线，最大沉降值稳定为-10.25mm(图 3)。图 4 可见地表监测断面距开挖面不同距离时的横向沉降槽的空间分布特性：随着盾构的逐渐推进地表沉降值逐渐增大，沉降槽的宽度也在变大。离开开挖面 36m 和 42m 时监测断面的最大沉降分别为-9.58mm 和-
15、9.65mm。图 3 和图 4 对比可知，沿隧道轴向，离开开挖面达到 20m 左右时，地表沉降已基本趋于稳定。 图 5 地表横向沉降的计算值与监测值对比 图 5 为地表横向沉降的计算值与监测值对比，图中计算值与实测值反映的规律基本一致，在隧道中轴线右 510m 范围内计算值与监测值非常接近，两者的差值不超过 2mm；在与隧道中轴线距离超过 10m 的范围，计算结果与监测结果相比，量值和影响范围均偏大，沉降曲线的最大绝对误差达 4mm，经分析认为这是由于计算模型简化以及计算中未能有效地评价盾尾同步注浆对地层损失的补偿作用造成的。 结论 本文在富水砂层盾构隧道掘进地面沉降问题研究中，数值模拟计算与
16、现场沉降监测数据基本一致，证明了数值模拟研究盾构隧道掘进引发地面沉降问题的合理性。 盾构在富水砂层推进过程中，在隧道开挖面前方，离开挖面越远位移量越小，距开挖面 20m 以外的土体几乎不产生位移，开挖面后方，地表沉降表现为距开挖面越远量值越大，但趋势渐缓。 模拟盾构隧道埋深 17m，计算出盾构开挖横向沉降槽影响范围大致在2025m，与监测数据较吻合，且离隧道纵断面中心线越远，盾构施工引起的地表沉降越小。 参考文献 纪佑军，刘建军，程林松考虑流固耦合的隧道开挖数值模拟J岩土力学，2011，32(4)：12291234 苑莲菊，李振全，武胜忠等工程渗流力学及应用M北京中国建材工业出版社，2001 张云，殷宗泽，徐水福盾构法隧适引起的地表变形分析J岩石力学与工程学报，2002，21(3)：388392 陈育民，徐鼎平FLAC/FLAC3D 基础与下程实例M北京：中国水利水电出版社，2009. 孙书伟，林杭，任连伟. FLAC3D 在岩土工程中的应用M北京：中国水利水电出版社，2011 作者简介：张鹏(1986.2)，男，工程师，从事城市轨道交通风险咨询工作及建筑工程的检测、科研工作
