Failure analysis of tunnel face influenced by overload and longitudinal slope
-
摘要: 针对地表超载和纵向坡度影响下的隧道开挖面稳定性问题，提出源于构造破坏机制的初始三角形网格更新策略，应用刚体平动运动单元上限分析（UB-RTME）方法开展非线性规划运算，获得地表临界超载系数σs/c图表与滑移线网破坏模式。结果表明：①以网格更新得到的UB-RTME滑移线网与Optum G2自适应加密网格，均可反映隧道开挖面破坏影响范围及承载关键区域，且UB-RTME以有效间断线及单元速度确定破坏范围，可揭示相对破坏范围S/D2定量数据及滑移线网破坏形态与影响因素的关联规律。②临界超载系数σs/c随纵向坡度角θ增大而显著降低，随埋深比C/D增大而增大，随重度系数γD/c增大而不断减小，随内摩擦角ϕ增大而非线性增大。③对于C/D和ϕ数值较大的情况，网格密度将显著影响UB-RTME运算得到的σs/c数值解精度；网格密度增加后，σs/c计算结果精度显著提高，且体现破坏特征的网状滑移线变得更加光滑。Abstract: Aiming at the stability of excavation face of tunnels under the influence of surface overload and longitudinal slope, an initial triangular mesh updating strategy derived from the structural failure mechanism is proposed. The computation is conducted by using the upper bound method with the rigid translatory moving elements (UB-RTME) with nonlinear programming to generate stability charts for the load factor σs/c and obtain the slip line failure modes. The results show that the range affected by the failure modes and the major failure domain are given by both the UB-RTME slip line mesh and the Optum G2 adaptive encryption mesh. The range of failure area is determined by the effective discontinuities and element velocity characterized by the UB-RTME. The quantitative data of the failure area S/D2 and the relationship between the slip line failure modes and the influencing factors are provided. The load factor σs/c decreases significantly with the increase of the longitudinal slope angle θ, increases with the increasing buried depth ratio C/D. It decreases with the increase of the gravity coefficient γD/c while increases nonlinearly with the increase in the internal friction angle ϕ. For the case with large values of C/D and ϕ, the mesh density will significantly affect the accuracy of the σs/c numerical solution obtained by the UB-RTME operation. With the increase of the mesh density, the accuracy of the calculated results of σs/c is significantly improved, and the slip line reflecting the failure characteristics becomes smoother.
-
Keywords:
- surface overload /
- longitudinal slope /
- tunnel face /
- stability analysis /
- mesh update /
- failure mode /
- slip line network
-
0. 引言
泥水平衡盾构通过向密封土仓内注入泥浆来平衡盾构开挖面的土水压力，广泛适用于大直径隧道盾构施工工程。盾构泥浆在长时间掘进过程中比重与黏度会不断增大，为了确保盾构开挖面的稳定，需要对泥水进行分离处理，进而控制泥浆指标[1-2]。现有的盾构泥水分离技术存在分离效率低、分离后泥渣含水高、废液污染等问题。因此，如何高效、经济的通过泥水分离处理盾构泥浆已成为泥水盾构施工过程中的一项重大技术难题。
盾构泥水分离技术是多步骤协同完成的复杂工艺过程，包括振动筛分、多级旋流、压滤与离心等步骤[3-4]。这些泥水分离步骤原理不尽相同，但都需要絮凝剂的添加这一关键步骤。絮凝剂可以使泥浆中的胶体和颗粒物凝聚、沉降和去除，为后续的脱水等步骤创造条件[5]。絮凝剂种类多种多样，在泥水分离方面运用最多的是一些传统的无机或者有机絮凝剂[6-8]。李春林等[9]针对泥水盾构穿越黏粉土地层产生的大量高含水率废弃泥浆脱水难的问题，分别进行了阴离子型聚丙烯酰胺（APAM）、阳离子型聚丙烯酰胺（CPAM）和非离子型聚丙烯酰胺（NPAM）等PAM类有机絮凝剂的泥浆絮凝试验，总结了泥水盾构中PAM类有机絮凝剂的作用规律与使用方案。王东星等[10]利用自制真空预压抽滤装置对5种絮凝剂进行脱水试验，发现有机试剂APAM的絮凝效果最好，可以快速实现淤泥脱水减容，无机试剂HCA影响效果有限，而无机试剂PAFSI和PAC絮凝效果相近。除了以上的传统有机或者无机絮凝剂，近年来也有很多的新型复合絮凝剂运用在泥水分离方面[11-12]。郭利芳等[13]研发了一种新型复合絮凝剂（S003），通过脱水固化试验，发现与PFS和PAM相比，S003处理的疏浚底泥脱水速度更快，具有更好的除浊效果。梁高杰等[14]为改善进口铝土矿溶出赤泥沉降性能，自制一种新型氧肟酸型絮凝剂，泥水分离后上清液澄清度以及沉降速度优于传统工业絮凝剂。
目前，虽然针对泥水分离中的絮凝试剂研究较多，但鲜有进行泥浆沉淀后的脱水性能研究。经调研，现场泥浆脱水后的泥渣存在含水率过高、难以及时外运等问题。同时，分离后的废浆也存在含泥量高、外排不达标和污染环境等问题。针对上述问题，本文依托广东湛江湾海底隧道工程，研发了一种新型H型脱水试剂，使得分离后的泥渣含水率低、废浆含泥量低，并通过室内絮凝脱水试验、微观机理分析和现场试验等，对比分析了H型试剂和PAM类絮凝剂分离效果的差异性及H型试剂的适用性。
1. 工程概况
1.1 基本情况
广东湛江湾海底隧道是广湛高铁全线控制性工程。如图 1所示，全长8230 m，隧道最大埋深约50 m。隧道开挖直径14.3 m，是目前独头掘进最长的大直径跨海高铁隧道。盾构长距离穿越第四系全新统粉质黏土，地层的整体的均匀性较差。
1.2 泥水处理工艺
广湛铁路湛江湾海底隧道工程泥水处理系统设置在地面，由泥浆循环系统、泥水分离系统、制调浆系统以及泥水监测系统4个部分组成。
现场的废弃泥浆处理离心工艺流程与离心设备如图 2所示，本研究的现场试验将基于该离心设备。离心机启动后，利用泥浆泵将废弃泥浆从泥浆池中抽吸至4台离心机中，当泥浆泵向离心机中进泥时，打开试剂泵，将絮凝试剂抽吸至泥浆管道中，使试剂与泥浆在管道中混合，再进入离心机中脱水处理，达到泥水分离的效果。离心机设备的处理能力为60~100 m³/h。
1.3 问题分析与解决思路
现阶段（2022年7月至今）盾构穿越区域主要为砂质泥岩，现场离心分离效果一般，废浆即使处理后仍然存在比重黏度大、难于沉淀等问题。废浆性质稳定，几乎很难在自然状态下沉淀。为解决上述的问题，本文研发的一种新型H型脱水试剂。试剂由本团队通过大量现场试验研配而成，其性质以及原现场使用的脱水试剂的详细信息如表 1所示。
表 1 试剂信息Table 1. Information of reagent材料 原现场试剂 H型试剂 组成 一种离子型聚丙烯酰胺，化学式为(C3H5NO)n，分子量：1500万，产地：广州 由改性铝系高聚物，铁系高聚物以及有机高分子聚合物（分子量580万）3种试剂按照4∶8∶1的比例复配而成，需要现配现用 性质 无色无味，高黏性，呈水玻璃状，有少许气泡，遇光发生降解 呈淡黄色，无异味，偏碱性（ph8.1），低黏性（26 mPa.s），有红色絮状沉淀 成本 15000元/t 4200元/t 实图 
2. 室内试验
2.1 室内试验设计
首先采集现场分离前的废弃泥浆（位于旋流后的沉淀池内）。运输过程中，由于蒸发等原因导致水分流失，泥浆密度为1.091g/cm3，含固率为13.4%。为了使泥浆和现场接近，将泥浆的含固率稀释为1.47%。后续将先后开展絮凝剂最优添加量、种类筛选试验、泥浆压滤脱水试验以及机理试验，对比分析H型试剂和现场常用的PAM类絮凝剂在分离效果和作用机理方面的差异性。
（1）絮凝剂最优添加量、种类筛选试验
试验采用絮凝沉降试验的方式，以浊度仪（量程1000NTU）测量的上清液浊度以及量筒读出的泥水分离体积比作为评价指标，确定每种絮凝剂的最优添加量（添加量=试剂干重/泥浆中固体干重）。对4种絮凝剂按顺序完成以下试验：取泥浆置于5个1000 mL试管中，将不同添加量的絮凝剂溶液分别加入盛有50 mL蒸馏水的烧杯中，每种试剂的添加量来源于预试验与工程实际，具体添加量见表 2。将50 mL絮凝剂溶液和泥浆混合后快速（100 r/min）搅拌1 min，静置沉淀1 h后测试上清液的浊度与泥水分离体积比（体积比=上清液体积/泥浆总体积）。
表 2 最优添加量筛选结果汇总Table 2. Summary of optimal additive filtering results种类 组号 添加量/% 体积比/% 浊度/NTU 阳离子型聚丙烯酰胺CPAM 1 0.05 21.67 489.6 2 0.10 30.83 200.3 3 0.25 31.50 189.7 4 0.50 24.50 176.2 5 1.00 23.33 253.8 阴离子型聚丙烯酰胺APAM 6 0.05 12.83 279.1 7 0.10 23.33 187.2 8 0.25 29.33 135.7 9 0.50 35.50 299.3 10 1.00 35.00 366.3 非离子型聚丙烯酰胺NPAM 11 0.05 2.17 — 12 0.10 8.54 — 13 0.25 20.50 293.6 14 0.50 6.26 — 15 1.00 2.22 — H型试剂 16 0.05 12.58 26.2 17 0.10 45.69 10.5 18 0.25 51.67 8.6 19 0.50 46.58 270.3 20 1.00 47.16 560.3 注：无数据说明分离的上清液的体积过少，浊度仪无法测量 基于每种絮凝剂的最优添加量，确定最优絮凝剂种类。按照上述试验类似的方法，每隔20 min记录一次上清液的浊度与泥水分离体积比，共120 min。
（2）泥浆压滤脱水试验
将筛选出的絮凝剂作用后的泥浆采用抽真空负压的方式进行脱水试验[15-17]。利用图 3的简易脱水装置，取量筒下端的泥水分离后的300 mL泥浆，将泥浆与试剂充分搅拌混合后倒入脱水装置上方的玻璃器皿中，并在下方垫上滤纸，连接气泵，对加药处理后的泥浆进行抽滤脱水，记录60 min内脱除的水量，每5 min记录一次。最后测量泥渣的含水率。
（3）絮凝效果机理分析试验
采用扫描电镜SEM（包括EDS面扫）、Zeta电位分析仪、马尔文2000激光粒度仪、旋转黏度仪等仪器，研究试验中出现的絮凝效果差异机理。
2.2 试验结果
（1）最优添加量
为确定4种絮凝剂的最优添加量，分别进行多组絮凝沉淀试验。表 2为1 h后泥浆添加絮凝剂后试验结果。图 4~7为各组絮凝效果图及对应的浊度与体积比变化图。
图 4为CPAM作用后的泥浆。随着投加量的增加，泥水分离体积比呈现先增后减的趋势。添加量在0.1%~0.25%时达到最大体积比31.50%。上清液浊度曲线近似为一条抛物线，添加量在0.1%~0.5%时浊度较低。综合两者，CPAM与泥浆作用的最优添加量选择在0.1%~0.25%。
图 5为APAM作用后的泥浆。投加量超过0.5%后，体积比基本保持在32%上下，此时达到了最大的泥浆反应添加量。泥浆上清液浊度呈现先减后增的趋势，添加量在0.25%时上清液浊度最低，为135.7 NTU。综合两者，APAM与泥浆作用的最优添加量选择在0.25%~0.5%。
图 6为NPAM作用后的泥浆。除0.25%组，其余组泥水分离体积比变化不明显。上清液浊度由于没有分离出足量的清液，只测取了0.25%添加量的一组。NPAM的絮凝效果相比其他絮凝剂的作用微弱，因此不选择NPAM作为本工程的絮凝剂。
图 7为H型试剂作用后的泥浆。除0.05%组，其余组的泥水分离体积比都很明显。0.25%添加量的上清液浊度最低，约8.6NTU。在1%的添加量时，非常浑浊，原因可能是过量的絮凝剂没有被溶解。因此H型试剂的最优添加量选择在0.1%~0.25%。
（2）最优种类
基于以上的最优添加量继续开展最优种类筛选试验。CPAM、APAM、H型试剂与泥浆作用的最优添加量分别选择0.1%，0.25%，0.1%。图 8为20 min与2 h时3种絮凝剂分别作用后的效果图。图 9为浊度与泥水分离体积比随时间变化的趋势图。
从图 8中可以看出，20 min前CPAM与APAM的作用效果相近。泥浆中加入CPAM、APAM后，出现了明显的泥水分离现象，但上清液仍然较浑浊。如图 9（a）所示，即使60 min后，浊度也高达230 NTU。H型试剂前期絮凝沉降速度较慢，但随着时间的推移，泥浆絮凝的效果越来越明显，且在40 min后，泥水分离体积比超过了CPAM与APAM。H型试剂作用的上清液很清澈，保持在20~40 NTU。综上可知H型试剂的泥水分离效果最好，可降低约85%的上清液浊度，CPAM次之，APAM最差。
（3）泥浆压滤脱水试验结果
图 10为添加3种絮凝剂后滤液体积随时间变化的趋势图。滤液体积反映了泥浆脱水的难易程度。图 12左侧为3种试剂压滤脱水产出的泥渣表观图。
APAM与CPAM作用后的泥浆有着类似的压滤脱水效果，60 min内总共压滤出了约78 mL的液体。作用后的泥渣呈豆腐块状，触感黏滞且疏松多水，这种泥渣不利于直接运输处置。H型试剂的压滤脱水效果明显好于另外两者：其脱水效率在25 min左右达到最大。60 min内总共压滤出了约96.23 mL的液体。H型试剂作用后的泥渣呈流砂状，触感也较干，方便直接运输处理。表 3为60 min后测量的各组泥渣含水率。使用H型试剂可以有效提高泥浆压滤脱水的效率，降低18%左右的泥饼含水率。
表 3 压滤泥渣含水率Table 3. Water contents of filter press sludge絮凝剂种类 APAM CPAM H型试剂 泥渣含水率/% 35.5 33.2 17.4 EDS可以半定性地分析泥渣的元素含量。图 11为EDS测量的泥渣元素含量图。CPAM与APAM处理后的泥渣元素占比相似，但上清液的pH值呈碱性（约9.3），H型试剂处理后的泥渣，Al、Ca等元素含量都提高了3%~6%。这表明H型试剂的部分金属元素被泥浆中的黏土矿物吸附，且通过测试，压滤排出的尾水pH呈中性（约7.8），重金属元素也未超标，符合《污水综合排放标准》[18]中外排条件。
2.3 絮凝效果差异的机理分析
下面采用微观机理法，分析3种絮凝剂在分离速度、上清液浊度与脱水滤液体积差异的原因机理。
图 12为各组处理后的泥渣以及其SEM扫描电镜微观图（1000倍）。图 13为采用粒度仪测量出的絮凝剂絮凝后的粒度分布图。图 12（a）可以看出原泥浆的结构均匀，颗粒普遍较小，10 μm以下的颗粒占77.63%，此状态下的泥浆性质稳定，难于沉淀；图 12（b），（c）中，加入APAM与CPAM后，颗粒间通过吸附连接，出现了直径超过50 μm（分别占比28.35%与47.21%）的较大颗粒，颗粒间为点点接触，该状态下的泥渣呈现板块豆腐状；图 12（d）中，加入H型试剂后，小颗粒出现了聚集，其大小（直径超过50 μm占35.88%）介于原泥浆与PAM类泥浆，颗粒呈现片状结构，这是H型试剂中的无机矿物落入泥浆中，营造出了矿物骨架，使得颗粒间点面或者面面接触，该状态下的泥渣呈现类似砂状。小颗粒聚集成团是泥浆絮凝的基本原理之一，它是由于絮凝剂颗粒溶于水后，会形成网状结构，通过吸附桥架作用和网捕作用的方式聚集泥浆颗粒形成聚团（简称：桥网作用）[19]。在重力作用下，大颗粒的沉降速度快于小颗粒，这便是PAM类组泥浆前期絮凝速度快于H型试剂组泥浆的直接原因。相同的现象也可以在图 13粒径级配图中定量的解释，加入絮凝剂后，3类试剂的曲线相对原泥浆整体向右偏移，桥网作用使得粒径＜10 μm的细粒更容易聚集成大颗粒。偏移的趋势CPAM＞H型试剂＞APAM，CPAM与H型试剂的桥网作用强于APAM，原因是只有APAM携带阴离子，它会引起高分子絮凝剂长链周围阴离子之间静电斥力增强，造成长链之间相互胶连，絮凝剂吸附架桥能力降低[20]。由上可知，桥网作用的强弱：CPAM＞H型试剂＞APAM。
图 14为采用Zeta电位仪测量的3种絮凝剂不同添加量对应的泥浆Zeta电位变化趋势图。Zeta电位的数值反映了分散体系的稳定性，Zeta电位绝对值越大，体系越稳定。胶体稳定性与Zeta电位的关系如表 4[7]。原泥浆的ζ电位绝对值高达35mV，属于中等稳定性，泥浆难于自然沉淀。泥浆中加入各类絮凝剂后，ζ电位都发生了变化，这涉及到絮凝的另外一个机理：当絮凝剂溶于水后产生的多核羟基络合物等与土颗粒会发生氢键吸附和静电吸引[21]，这可以有效地中和土颗粒表面的负电，导致ζ电位下降，进而减少了土颗粒间的静电引力，使泥浆更容易沉淀（简称电中和作用）。添加APAM后，ζ电位绝对值先下降到32.1 mV后上升，原因是前期由于桥网作用，屏蔽了一部分的负电荷，但APAM带负电，会导致ζ电位慢慢上升，并且APAM的最优掺量0.2%也在电位曲线转折附近，这说明APAM在此时的电中和作用最强，但是相比另外两类絮凝剂，APAM产生絮凝的主要原因是桥网作用。添加CPAM与H型试剂后，可以看出Zeta电位很快发生了下降，但是H型试剂的整体降低趋势大于CPAM，在添加浓度为0.3%时，H型试剂的ζ电位绝对值降低幅度比CAPM平均高45%，因此H型试剂的电中和能力要强于CPAM。结合图 11可知，造成上述现象主要因为H型试剂处理后的泥渣中Al、Ca等浓度都高于CPAM，因此相同浓度下，H型试剂中的离子价态更高。由上可知，电中和作用的强弱：H型试剂＞CPAM＞APAM。
表 4 胶体中悬浮颗粒的稳定性与Zeta电位的关系Table 4. Relationship between stability of suspended particles and Zeta electric potential稳定性能 Zeta电位/mV 最高量聚集和沉淀 0~+3 强聚集和沉淀的范围 +5~-5 聚集开始 -5~-15 分散体的灵敏值 -16~-30 中等稳定性 -31~-40 相当好的稳定性 -41~-60 很好的稳定性 -61~-80 上文通过桥网作用与电中和作用两个机理角度，解释了3种絮凝剂在泥水分离效果方面差异的原因。下面采用微观分析与泥浆黏滞性理论结合的方式分析3种絮凝剂在压滤效果差异的原因。
图 15为压滤脱水后的泥渣（H型试剂处理）经烘干喷金后的某区域EDS能谱分析及其对应的元素含量。泥渣表面分布着较高含量的的Ca、Al、Si等元素。Ca来源于H型试剂，而Al、Si分别来源于泥浆黏土矿物中的Al2O3与SiO2。Ca在碱性环境下可与SiO2、Al2O3等发生如下的火山灰反应[20]：
该反应生成了不溶于水的硅酸钙与铝酸钙胶状物质，这些物质附着在泥浆颗粒表面形成纤维状的晶体，如图 15中的EDS面扫局部放大图。纤维晶体互相吸附连接形成了泥浆颗粒间的排水通道骨架。这些含有Ca的坚硬通道骨架不仅扩大了排水通道，还起到了一定的支撑作用，在压滤的作用下，延迟了排水通道的闭合，间接释放了更多的水分，提高了泥浆压滤脱水的效率。
图 16为泥浆在添加APAM、CPAM以及H型试剂后的黏度变化图。试验初始泥浆具有高达1536 mPa.s的黏性。黏性的出现是由于泥浆颗粒具有负电势，其硅酸盐矿物晶体表面具有负电势的氧原子会吸引水分子中具有正电势的氢原子，从而形成一层结合牢固的具有较高黏滞性水分子层[22]，黏滞性越强，泥浆颗粒之间就会有越强的连接吸附作用，泥浆在宏观上就会表现越高的黏性。加入APAM与CPAM后，溶液的黏度反而上升，可能是由于高分子链没有完全展开，发生了缠结，其相互运动受到了空间上的阻碍[21]。如图 17（b）所示，较高的泥浆黏滞性不仅会堵塞泥浆颗粒间排水的通道，降低排水效率，而且还不利于微小的泥浆絮团沉降，导致了APAM与CPAM作用后的泥浆上清液浊度低于H型试剂。但是加入H型试剂后，泥浆的黏度降低到了1422 mPa.s，如图 17（c）所示，由于H型絮凝剂分子具有正电势，一方面中和并屏蔽了泥浆颗粒表面的负电势，另一方面可以减少具有高黏滞性的水分子层的形成，进而增加了泥浆颗粒之间的结合水释放。同时，由于泥浆的黏性下降，泥浆分子之间因相对运动或流动所产生的内摩擦力减小，润滑度上升，能够使得泥浆中的水分在抽滤负压作用下，更快通过颗粒间的渗水通道排出。
综上，絮凝效果的优劣同时受桥网与电中和等作用的影响，是一个复杂的过程。H型试剂兼具了较好的桥网作用与电中和作用，同时营造了渗水性更强的排水通道，降低了泥浆的黏滞性，提高了泥浆压滤脱水的效率，因此H型试剂的使用效果优于APAM与CPAM，它是适用于本工程泥水分离的优秀絮凝试剂。
3. 现场生产运用
3.1 现场原试剂对比
湛江现场首先进行了原试剂（原现场脱水试剂）与H型试剂的泥浆处理效果对比预试验。图 18为两种试剂静置2 h后的处理效果对比图。H型试剂处理后的泥浆上清液更为清澈，且沉淀后的泥浆呈紧密的砂状，原现场脱水试剂处理后的泥浆上清液较浑，且沉淀后的泥浆呈疏松的絮状。
图 19为2种试剂的试验数据。H型试剂的前期沉降速度较慢，60~80内，H型试剂的沉降体积超过原现场脱水试剂，这与前文室内试验结果相似。H型试剂的上清液浊度保持着较低水平，在18NTU左右，远低于原现场脱水试剂。H型的泥水分离也更为彻底，2 h后上清液的体积为263.5 mL。
3.2 生产试验设计
现场分别采用原现场脱水试剂和H型试剂对20 m3泥浆进行处理试验。将粉状试剂在溶药箱（容积约900 L）中溶解，采用环流的方式对掺入外加剂的溶液进行搅拌混合，同时使用高压水枪冲解凝结成团的药团使其充分溶解。待试剂溶解后，溶液通过试剂泵输入进容积20 m3的小型泥浆池。试验过程中不断记录泥浆浓度等参数。之后将废弃泥浆输入至小型泥浆池中，在泥浆抽吸过程中，同时将溶解好的液体试剂抽吸至泥浆池进泥口，使试剂与泥浆充分混合。分别采集加入两种试剂后的泥浆脱水泥渣和尾水进行对比分析。图 20为泥水分离现场试验的全流程实拍图。
3.3 试验结果
（1）进泥参数
表 5，6记录了两组试验的泥浆指标。由于试验准备需要2 d，含固率出现了波动。为了保证室内试验与现场试验的结果可比性，同时考虑现场试验中药池底部未抽取的药剂，根据预试验的结果，添加量需要上浮20%左右。
表 5 原现场脱水试剂生产试验时的泥浆指标Table 5. Slurry indices during PAM agent production trials泥浆指标 试验阶段 平均值 进泥前期 试验过程中 密度/(g·cm-3) 1.048 1.008 1.0076 1.02 1.02 含固率/% 7.36 1.28 1.21 3.15 3.25 表 6 H型试剂生产试验时的泥浆指标Table 6. Slurry indices during production tests with H-agent泥浆指标 试验阶段 平均值 进泥前期 试验过程中 密度/(g·cm-3) 1.0272 1.012 1.006 1.0048 1.01 含固率/% 4.25 1.9 0.96 0.77 1.97 （2）出泥结果
图 21为离心机运转过程中，经过离心脱水后的尾水（左H型试剂，右原现场脱水试剂）。可以看出H型试剂相比与原现场脱水试剂作用后的泥浆，其颜色偏淡黄，而原现场脱水试剂偏灰色。灰色为原本泥浆的颜色，因此添加原现场脱水试剂的尾水仍有一定的砂浆沉淀。原现场脱水试剂絮凝离心后的效果差于H型试剂。
图 22为试验过程中采集的尾水（左H型试剂首，右原现场脱水试剂）。静置2 min后的上清液可以看出添加H型试剂的泥浆出现了明显的矾花，而加入原现场脱水试剂的泥浆仍未析出上清液，因此H型试剂的泥水分离效果更好。试验还进行了上清液的PH值测量，发现H型试剂处理后的泥浆PH值偏中性，而原现场脱水试剂偏碱性，因此H型试剂处理后的上清液更符合排放要求。
图 23为泥浆经过离心脱水后的泥渣（左H型试剂，右原现场脱水试剂）。采用原现场脱水试剂试剂处理后的泥渣，筛分后状态类似于豆腐脑状态，在空气中长时间晾晒也无法完成脱水，其含水率仍然较高。采用H型试剂处理后的泥渣，泥样颗粒感明显，且晾晒后容易失水变干，方便渣土的运输处置。
对两种试剂处理后的泥渣含固率进行定量的测量，如表 7所示。H型试剂生产的脱水泥样含固率达到了70.3%，原现场脱水试剂生产的脱水泥样含固率仅为39.8%。由此可知，虽然原现场脱水试剂可以使低浓度泥浆中出现明显的矾花，但沉积在底部的浓缩泥浆含水率仍较高，且其中的水分难以脱除。
表 7 加入两种试剂生产的脱水泥样含固率Table 7. Solids content of de-cemented samples produced by试剂名称 脱水泥样含固率/% 原现场脱水试剂 39.8 H型试剂 70.3 3.4 试剂用量和价格分析
表 8所示为H型试剂与原现场脱水试剂的现场试剂用量与成本。使用H型试剂处理泥浆，虽然需要更高的用量，但其单价4200元/t相比原现场脱水试剂的15000元/t更低。综合成本H型试剂也更低，约39.94元/m³。综上，与原现场脱水试剂相比，现场使用0.5%的H型试剂，按照9.5 kg/t的用量处理泥浆，可提高30.5%的脱水泥样含固率，降低20%的工程成本。
表 8 试剂用量和价格分析Table 8. Analysis of dosage and prices of reagentsH型试剂 原现场脱水试剂 用量/
(kg·m-³)单价/
(元·m-³)用量/
(kg·m-³)单价/
(元·m-³)9.51 39.94 3.19 47.92 注：①H型试剂按照当前市场价4200元/t计算、原现场脱水试剂单价按照15000元/t计算；②成本已考虑由于现场设备等原因未溶解的试剂。 4. 结论与展望
针对广湛铁路湛江湾海底隧道工程盾构废弃泥浆处理工艺中的不足，开展了泥浆处理现场试验研究，得到以下3点结论。
（1）室内试验中，CPAM、APAM、H型试剂与泥浆作用的最优添加量分别为0.1%，0.25%，0.1%。添加H型试剂后的泥浆前期絮凝速度要慢于CPAM与APAM，但脱水后的泥渣含水率降低了约18.1%，且上清液非常清澈，可降低约85%的上清液浊度，ph7.8也符合直接外排条件。
（2）絮凝机理试验中，CPAM与APAM处理后的泥浆颗粒呈现黏粒状的微观结构，而H型试剂呈现片状微观结构。3类试剂的桥网作用与电中和作用使得粒径＜10 μm的细粒更容易聚集成大颗粒，加快了泥浆的絮凝速度，超过50 μm的颗粒分别占比28.35%，35.88%，47.21%。泥浆桥网作用的强弱：CPAM＞H型试剂＞APAM，电中和作用的强弱：H型试剂＞CPAM＞APAM。H型试剂兼具了较好的桥网作用与电中和作用，也可以营造渗水性更强的Ca排水通道，同时降低泥浆的黏滞性（下降约7.42%），提高泥浆压滤脱水的效率。
（3）现场生产试验中，与原现场脱水试剂相比，使用0.1%的H型试剂，按照9.5kg每吨的用量处理泥浆，可提高30.5%的脱水泥样含固率，降低20%的工程成本（约39.94元/m³）。采用H型试剂处理后的泥渣，泥样颗粒感明显，晾晒后更容易失水变干，方便渣土的运输处置。综上，H型试剂是适用于本工程泥水分离的优秀絮凝试剂。此外，针对H型试剂初始反应相比PAM类略慢但泥水分离更为彻底的特点，其也适合在压滤机领域应用，未来可在有关方向做更深入研究。
-
![]()
图 1 考虑超载和纵坡效应的隧道开挖面稳定性分析模型
Figure 1. Stability analysis of model for tunnel face influenced by overload and longitudinal slope
![]()
图 2 运动单元上限分析网格更新过程样例图(θ =0°, C/D =3, γD/c =0, ϕ =15°)
Figure 2. Sample charts of mesh updates with UB-RTME(θ =0°, C/D =3, γD/c =0, ϕ =15°)
![]()
图 3 非结构化初始网格更新样例图(θ =0°, C/D =3, γD/c =0, ϕ =15°)
Figure 3. Sample charts of unstructured initial mesh update(θ =0°, C/D =3, γD/c =0, ϕ =15°)
![]()
图 4 超载和纵坡影响下的隧道开挖面临界超载系数σs/c上限解曲线图
Figure 4. Curves of variation of upper bound solution of σs/c for tunnel face influenced by overload and longitudinal slope
![]()
图 5 隧道开挖面滑移线网破坏模式及网格变形特征图(ϕ=15°, γD/c=2, C/D=3, θ=0°)
Figure 5. Slip line failure mechanisms and mesh deforming characteristics for tunnel face (ϕ=15°, γD/c=2, C/D=3, θ=0°)
![]()
图 6 埋深比C/D影响下的隧道开挖面滑移线网破坏模式
Figure 6. Slip line failure mechanisms for tunnel face influenced by C/D
![]()
图 7 重度系数γD/c和内摩擦角ϕ影响下的隧道开挖面滑移线网破坏模式
Figure 7. Slip line failure mechanisms of tunnel face influenced by γD/c and ϕ
![]()
图 8 网格疏密影响下的隧道开挖面滑移线网破坏模式(θ=0°, ϕ=25°, C/D=2, γD/c=2)
Figure 8. Slip line failure mechanisms for tunnel face influenced by mesh density (θ=0°, ϕ=25°, C/D=2, γD/c=2)
表 1 隧道开挖面稳定性分析计算参数
Table 1 Dimensionless parameters for stability of tunnel face
计算参数 参数取值 C/D 1, 2, 3, 4 γD/c 0, 1, 2, 3 ϕ/(°) 10, 15, 20, 22.5, 25 θ/(°) -10, 0, 10 
下载: 导出CSV
表 2 运动单元上限分析网格更新操作
Table 2 Mesh update operation of UB-RTME method
序号 名称 网格原状 网格更新 操作说明 1 移动节点 
修改坐标调整节点位置 2 新增节点 
在公共边按长度比例新增节点，添加间断线 3 单元缩减 
合并单元缩减单元 4 对角互换 
相邻单元公共边互换 5 单元加密 
单元一分为四加密 6 成层加密 
沿指定方向成层加密
下载: 导出CSV
表 3 隧道开挖面破坏临界超载系数σs/c上限解
Table 3 Upper bound solutions of critical load factor σs/c for failure of tunnel face
C/D ϕ/(°) γD/c (θ=10°) γD/c (θ=0°) γD/c (θ=-10°) 0 1 2 3 0 1 2 3 0 1 2 3 10.0 7.47 5.50 3.50 1.47 8.17 6.25 4.29 2.29 9.01 7.16 5.25 3.30 15.0 11.33 8.90 6.38 3.79 12.62 10.29 7.87 5.35 14.24 12.07 9.78 7.37 1 20.0 19.22 15.93 12.49 8.87 21.93 18.89 15.65 12.20 25.60 22.93 20.01 16.84 22.5 26.40 22.40 18.18 13.70 30.67 27.04 23.13 18.92 36.36 33.33 29.94 26.20 25.0 38.06 32.99 27.58 21.80 44.90 40.42 35.52 30.19 54.60 51.11 47.08 42.54 10.0 11.34 7.87 4.30 0.63 12.11 8.75 5.27 1.67 13.01 9.78 6.42 2.93 15.0 19.45 14.92 10.16 5.13 21.15 16.89 12.35 7.53 23.16 19.25 15.03 10.50 2 20.0 39.13 32.45 25.26 17.49 43.61 37.63 31.05 23.84 49.07 44.06 38.35 31.97 22.5 60.05 51.35 41.89 31.56 67.97 60.44 52.02 42.67 77.78 71.94 65.06 57.17 25.0 98.40 86.30 73.02 58.41 113.40 103.36 91.95 79.09 132.34 125.40 116.84 106.70 10.0 14.43 9.44 4.27 -1.14 15.24 10.41 5.37 0.10 16.19 11.57 6.69 1.58 15.0 26.89 20.20 13.06 5.36 28.91 22.66 15.81 8.56 31.22 25.53 19.27 12.44 3 20.0 60.84 50.41 38.98 26.40 66.93 57.75 47.45 35.97 74.10 66.61 57.85 47.87 22.5 100.59 86.36 70.61 53.12 112.27 100.30 86.61 71.11 126.21 117.33 106.52 93.82 25.0 179.87 158.75 135.21 108.88 204.11 187.32 167.74 145.25 233.42 222.73 208.87 191.93 10.0 17.06 10.55 3.72 -3.50 17.88 11.60 4.96 -2.06 18.90 12.86 6.43 -0.37 15.0 33.88 24.97 15.36 4.84 36.15 27.85 18.79 8.84 38.77 31.24 22.87 13.62 4 20.0 83.79 69.34 53.33 35.48 91.38 78.83 64.54 48.40 100.30 90.21 78.21 64.33 22.5 146.57 126.17 103.33 77.66 161.81 144.98 125.43 103.01 180.38 168.30 153.21 135.18 25.0 279.81 248.03 212.19 171.67 313.31 288.78 259.67 228.22 355.36 340.68 320.89 296.10
下载: 导出CSV
表 4 隧道开挖面破坏临界重度系数γcrD/c对比数据
Table 4 Comparison of critical gravity load factor γcrD/c for tunnel face
C/D ϕ/(°) 文献[12] 变形单元 文献[13] 运动单元 Optum G2结果 本文解 γcrD/c σs/c 1 10 4.00 4.04 3.97 4.00 0.255 1 15 4.89 4.88 4.77 4.90 0.264 1 20 5.90 5.89 5.73 5.99 0.247 2 10 3.48 3.41 3.37 3.40 0.199 2 15 4.39 4.36 4.28 4.40 0.180 2 20 5.62 5.55 5.43 5.67 0.137 3 10 3.05 2.99 2.96 3.00 0.102 3 15 4.08 4.01 3.95 4.05 0.112 3 20 5.37 5.34 5.24 5.45 0.193
下载: 导出CSV
表 5 隧道开挖面相对破坏范围S/D2计算结果
Table 5 Computated results for failure area S/D2 of tunnel face
C/D ϕ/(°) γD/c (θ=10°) γD/c (θ=0°) γD/c (θ=-10°) 0 1 2 3 0 1 2 3 0 1 2 3 1 10 2.17 2.05 1.93 1.81 2.36 2.21 2.06 1.93 2.54 2.35 2.19 2.04 15 2.16 2.03 1.90 1.76 2.35 2.19 2.04 1.89 2.59 2.39 2.21 2.05 20 2.15 2.04 1.92 1.79 2.35 2.21 2.08 1.94 2.57 2.41 2.25 2.11 25 2.17 2.08 1.99 1.88 2.37 2.26 2.16 2.05 2.60 2.47 2.36 2.24 2 10 5.70 5.23 4.78 4.29 6.17 5.60 5.08 4.56 6.73 6.05 5.47 4.90 15 5.74 5.30 4.85 4.37 6.21 5.69 5.20 4.67 6.76 6.15 5.62 5.08 20 5.78 5.44 5.08 4.69 6.24 5.84 5.45 5.05 6.81 6.35 5.92 5.52 25 5.86 5.63 5.38 5.11 6.33 6.06 5.79 5.52 6.89 6.59 6.30 6.03 3 10 10.97 9.85 8.75 7.47 11.77 10.46 9.26 7.97 12.68 11.15 9.87 8.57 15 11.08 10.11 9.12 7.97 11.85 10.74 9.70 8.60 12.83 11.57 10.41 9.32 20 11.18 10.50 9.79 8.99 11.96 11.20 10.44 9.66 12.94 12.05 11.24 10.46 25 11.37 10.94 10.50 10.01 12.15 11.69 11.23 10.75 13.15 12.61 12.11 11.64 4 10 17.99 15.89 13.79 10.86 19.22 16.78 14.51 11.77 20.53 17.75 15.40 12.81 15 18.20 16.48 14.71 12.46 19.29 17.41 15.55 13.52 20.83 18.64 16.64 14.71 20 18.40 17.27 16.08 14.71 19.62 18.35 17.11 15.81 21.01 19.59 18.27 17.02 25 18.72 18.07 17.40 16.63 19.93 19.19 18.48 18.48 21.38 20.56 19.81 19.09
下载: 导出CSV
-
[1] 严松宏, 李国良. 隧道与地下工程荷载计算的研究现状[J]. 兰州交通大学学报, 2020, 39(2): 1-7. doi: 10.3969/j.issn.1001-4373.2020.02.001 YAN Songhong, LI Guoliang. Research status of the load calculation of tunnel and underground engineering[J]. Journal of Lanzhou Jiaotong University, 2020, 39(2): 1-7. (in Chinese) doi: 10.3969/j.issn.1001-4373.2020.02.001
[2] 李修磊, 李金凤, 杨超. 黏土地层浅埋盾构隧道开挖面三维稳定性上限分析[J]. 铁道学报, 2021, 43(4): 166-174. doi: 10.3969/j.issn.1001-8360.2021.04.021 LI Xiulei, LI Jinfeng, YANG Chao. Upper-bound solutions for three-dimensional excavation face stability of shallow shield tunnel in cohesive-frictional soils[J]. Journal of the China Railway Society, 2021, 43(4): 166-174. (in Chinese) doi: 10.3969/j.issn.1001-8360.2021.04.021
[3] 张明告, 周顺华, 黄大维, 等. 地表超载对地铁盾构隧道的影响分析[J]. 岩土力学, 2016, 37(8): 2271-2278. ZHANG Minggao, ZHOU Shunhua, HUANG Dawei, et al. Analysis of influence of surface surcharge on subway shield tunnel under[J]. Rock and Soil Mechanics, 2016, 37(8): 2271-2278. (in Chinese)
[4] 黄大维, 周顺华, 冯青松, 等. 地表均布超载作用下盾构隧道上覆土层竖向土压力转移分析[J]. 岩土力学, 2019, 40(6): 2213-2220. HUANG Dawei, ZHOU Shunhua, FENG Qingsong, et al. Analysis for vertical earth pressure transference on overlaying soils of shield tunnel under uniform surface surcharge[J]. Rock and Soil Mechanics, 2019, 40(6): 2213-2220. (in Chinese)
[5] 曾明生, 王文法, 康海波, 等. 基于不同破坏模式的盾构隧道掌子面支护力上限分析[J]. 铁道科学与工程学报, 2022, 19(8): 2360-2368. ZENG Mingsheng, WANG Wenfa, KANG Haibo, et al. Upper limit analysis of support force on shield tunnel face based on different failure modes[J]. Journal of Railway Science and Engineering, 2022, 19(8): 2360-2368. (in Chinese)
[6] HUANG M S, SONG C X. Upper-bound stability analysis of a plane strain heading in non-homogeneous clay[J]. Tunnelling and Underground Space Technology, 2013, 38: 213-223. doi: 10.1016/j.tust.2013.07.012
[7] 黄茂松, 宋春霞, 王浩然. 基于上限定理的软土隧道开挖面稳定性分析[J]. 防灾减灾工程学报, 2014, 34(3): 330-335. HUANG Maosong, SONG Chunxia, WANG Haoran. Upper bound limit analysis for face stability of shield tunnel[J]. Journal of Disaster Prevention and Mitigation Engineering, 2014, 34(3): 330-335. (in Chinese)
[8] 黄茂松, 宋春霞, 吕玺琳. 非均质黏土地基隧道环向开挖面稳定上限分析[J]. 岩土工程学报, 2013, 35(8): 1504-1512. http://cge.nhri.cn/article/id/15259 HUANG Maosong, SONG Chunxia, LÜ Xilin. Upper bound analysis for stability of a circular tunnel in heterogeneous clay[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(8): 1504-1512. (in Chinese) http://cge.nhri.cn/article/id/15259
[9] SURAPARB K, BOONCHAI U. Design equation for stability of a circular tunnel in anisotropic and heterogeneous clay[J]. Underground Space, 2022, 7(1): 76-93. doi: 10.1016/j.undsp.2021.05.003
[10] 阳军生, 张箭, 杨峰. 浅埋隧道掌子面稳定性二维自适应上限有限元分析[J]. 岩土力学, 2015, 36(1): 257-264. YANG Junsheng, ZHANG Jian, YANG Feng. Stability analysis of shallow tunnel face using two-dimensional finite element upper bound solution with mesh adaptation[J]. Rock and Soil Mechanics, 2015, 36(1): 257-264. (in Chinese)
[11] YANG F, ZHANG J, YANG J S, et al. Stability analysis of unlined elliptical tunnel using finite element upper-bound method with rigid translatory moving elements[J]. Tunnelling and Underground Space Technology Incorporating Trenchless Technology Research, 2015, 50: 13-22.
[12] YANG F, ZHANG J, ZHAO L H, et al. Upper-bound finite element analysis of stability of tunnel face subjected to surcharge loading in cohesive-frictional soil[J]. KSCE Journal of Civil Engineering, 2016, 20(6): 2270-2279. doi: 10.1007/s12205-015-0067-z
[13] 孙雁军, 阳军生, 罗静静, 等. 隧道工作面稳定性与滑移线网破坏模式研究[J]. 岩土工程学报, 2019, 41(7): 1374-1380. doi: 10.11779/CJGE201907024 SUN Yanjun, YANG Junsheng, LUO Jingjing, et al. Stability and mesh-like collapse mechanism of tunnel face[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(7): 1374-1380. (in Chinese) doi: 10.11779/CJGE201907024
[14] 杨峰, 何诗华, 吴遥杰, 等. 非均质黏土地层隧道开挖面稳定运动单元上限有限元分析[J]. 岩土力学, 2020, 41(4): 1412-1419, 1436. YANG Feng, HE Shihua, WU Yaojie, et al. Tunnel face stability analysis by the upper-bound finite element method with rigid translatory moving element in heterogeneous clay[J]. Rock and Soil Mechanics, 2020, 41(4): 1412-1419, 1436. (in Chinese)
[15] WENG X L, SUN Y F, YAN B H, et al. Centrifuge testing and numerical modeling of tunnel face stability considering longitudinal slope angle and steady state seepage in soft clay[J]. Tunnelling and Underground Space Technology, 2020, 101: 103406. doi: 10.1016/j.tust.2020.103406
[16] CHENG C, JIA P J, ZHAO W, et al. Experimental and analytical study of shield tunnel face in dense sand strata considering different longitudinal inclination[J]. Tunnelling and Underground Space Technology, 2021, 113: 103950. doi: 10.1016/j.tust.2021.103950
[17] 马辉, 王飞, 刘泽挂, 等. 深埋特长大坡度斜井变形规律与影响因素分析[J]. 地下空间与工程学报, 2020, 16(增刊2): 950-956. MA Hui, WANG Fei, LIU Zegua, et al. Analysis of deformation law and influencing factors of deep-buried extra-long inclined shaft with large slope[J]. Chinese Journal of Underground Space and Engineering, 2020, 16(S2): 950-956. (in Chinese)
[18] SHI X M, LIU B G, TANNANT D, et al. Influence of consolidation settlement on the stability of inclined TBM tunnels in a coal mine[J]. Tunnelling and Underground Space Technology, 2017, 69: 64-71. doi: 10.1016/j.tust.2017.06.013
[19] 雷明锋, 彭立敏, 施成华, 等. 迎坡条件下盾构隧道开挖面极限支护力计算与分析[J]. 岩土工程学报, 2010, 32(3): 488-492. http://cge.nhri.cn/article/id/12413 LEI Mingfeng, PENG Limin, SHI Chenghua, et al. Calculation and analysis of limit support force of shield tunnel excavation face under facing-slope conditions[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(3): 488-492. (in Chinese) http://cge.nhri.cn/article/id/12413
[20] 周峻, 杨子松, 彭芳乐. 上坡条件下盾构开挖面极限支护压力研究[J]. 地下空间与工程学报, 2011, 7(5): 914-918. doi: 10.3969/j.issn.1673-0836.2011.05.016 ZHOU Jun, YANG Zisong, PENG Fangle. A study on the support pressure limit of the excavation face of shield tunnel under upslope conditions[J]. Chinese Journal of Underground Space and Engineering, 2011, 7(5): 914-918. (in Chinese) doi: 10.3969/j.issn.1673-0836.2011.05.016
[21] 赵智涛, 曹伍富, 王霆. 隧道爬坡开挖与水平开挖掌子面稳定性对比分析[J]. 现代隧道技术, 2018, 55(6): 94-100. ZHAO Zhitao, CAO Wufu, WANG Ting. Contrastive analysis of working face stability in cases of climbing excavation and horizontal excavation[J]. Modern Tunnelling Technology, 2018, 55(6): 94-100. (in Chinese)
[22] ZHAO L H, LI D J, LI L, et al. Three-dimensional stability analysis of a longitudinally inclined shallow tunnel face[J]. Computers and Geotechnics, 2017, 87: 32-48. doi: 10.1016/j.compgeo.2017.01.015
[23] 程诚, 赵文, 王迎超, 等. 密砂地层盾构隧道纵坡开挖面稳定性理论分析[J]. 中国公路学报, 2023, 36(4): 157-168. doi: 10.3969/j.issn.1001-7372.2023.04.014 CHENG Cheng, ZHAO Wen, WANG Yingchao, et al. Theoretical analysis of longitudinally-inclined shield tunnel face stability in dense sand stratum[J]. China Journal of Highway and Transport, 2023, 36(4): 157-168. (in Chinese) doi: 10.3969/j.issn.1001-7372.2023.04.014
[24] 张箭, 戚瑞宇, 宗晶瑶, 等. 盾构隧道环向开挖面破坏机制及剪胀效应研究[J]. 岩土力学, 2022, 43(7): 1833-1844. ZHANG Jian, QI Ruiyu, ZONG Jingyao, et al. Failure mechanism of shield tunnel circumferential excavation face and the influence of the dilatancy effect on the tunnel stability[J]. Rock and Soil Mechanics, 2022, 43(7): 1833-1844. (in Chinese)
-
其他相关附件
计量
- 文章访问数: 321
- HTML全文浏览量: 65
- PDF下载量: 64
