Establishment and application of two-dimensional axisymmetric temperature field model for high-level radioactive waste disposal repository
-
摘要: 高放废物处置库温度场长期演化特征是处置库缓冲层安全性能评估和设计、处置容器空间优化布局的重要依据之一。建立了处置单元二维轴对称双层热分析模型，以揭示处置库多重屏障系统近场温度演化规律。对传热控制方程应用拉普拉斯和有限傅里叶正弦变换，得到缓冲层和围岩层温度的拉普拉斯域解，借助Crump方法对拉普拉斯域解进行数值反演获得对应的时域解；通过与线热源解和数值解的对比验证了模型的正确性，分析了处置单元近场多屏障系统温度场的时空分布特征，探讨了不同参数对缓冲层峰值温度的影响；借助模型半解析解确定最小处置容器间距以及对原位试验结果进行预测。结果表明：①与线热源解相比，模型半解析解可以更准确地模拟缓冲层温度变化；②当处置隧道间距超过40 m时，其对缓冲层峰值温度影响较小；③在处置隧道间距取40 m时，最小处置容器间距为7.7 m；④模型半解析解可以较好地预测原位试验结果。Abstract: The long-term evolution of temperature field in high-level radioactive waste disposal repository is one of the important bases for the safety evaluation and design of buffer layer and the spatial optimization of disposal container. A two-dimensional axisymmetric thermal analysis model with two layers of disposal unit is established to reveal the near-field temperature evolution of multi-barrier system. The Laplace and finite Fourier sine transforms are applied to the governing equations for heat transfer to obtain the temperature solutions to the buffer layer and the surrounding rock layer in the Laplace domain. The corresponding solutions in time domain are obtained by numerical inversion of the Laplace domain solutions with the help of the Crump method. The correctness of the model is verified by comparing with the linear heat source solution and the numerical solution. The temporal and spatial distribution characteristics of the temperature field in the near-field multi-barrier system of the disposal unit are analyzed, and the effects of different parameters on the peak temperature of the buffer layer are discussed. The semi-analytical solution of the model is used to determine the minimum disposal container spacing and predict the results of the in-situ tests. The results show that the semi-analytical solution of the model can simulate the temperature change of buffer layer more accurately by comparing with the linear heat source solution. When the tunnel spacing exceeds 40 m, it has small effects on the peak temperature of buffer layer. When the tunnel spacing is 40 m, the container spacing is 7.7 m. The semi-analytical solution of the model can well predict the results of the in-situ test.
-
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类略慢但泥水分离更为彻底的特点，其也适合在压滤机领域应用，未来可在有关方向做更深入研究。
-
![]()
图 2 处置单元多重屏障系统数值分析模型
Figure 2. Numerical model for multi-barrier system in disposal unit
![]()
图 3 本文半解析解与其他方法的对比
Figure 3. Comparison between semi-analytical method in this study and other methods
![]()
图 4 不同影响半径下缓冲层温度变化
Figure 4. Temperature variation of buffer layer under different influence radii
![]()
图 5 处置单元多重屏障系统近场温度分布图
Figure 5. Temperature distribution of multi-barrier system in disposal unit
![]()
图 6 处置隧道间距对缓冲层峰值温度的影响
Figure 6. Effects of tunnel spacing on peak temperature of buffer layer
![]()
图 7 处置容器间距对缓冲层峰值温度的影响
Figure 7. Effects of container spacing on peak temperature of buffer layer
![]()
图 8 缓冲层厚度对缓冲层峰值温度的影响
Figure 8. Effects of thickness of buffer layer on its peak temperature
![]()
图 9 缓冲层热传导系数对缓冲层峰值温度的影响
Figure 9. Effects of thermal conductivity of buffer layer on its peak temperature
![]()
图 10 围岩热传导系数对缓冲层峰值温度的影响
Figure 10. Effects of thermal conductivity of rock on peak temperature of buffer layer
![]()
图 11 缓冲层峰值温度随处置容器间距变化图
Figure 11. Variation of peak temperature of buffer layer with different container spacings
![]()
图 14 原位试验结果与不同解之间的对比
Figure 14. Comparison between in-situ test results and different solutions
表 1 参数取值
Table 1 Values of parameters
参数 量值 参数 量值 λB/(W·(K·m)-1) 0.68 z2/m 502.625 λR/(W·(K·m)-1) 2.50 t0/℃ 10 cB/(J·(K·kg)-1) 1014 t1/℃ 34 cR/(J·(K·kg)-1) 832 P(0)/W 1545 ρB/(kg·m-3) 1 770 ϕ 0.91 ρR/(kg·m-3) 2 630 a1 0.049 a/m 0.525 a2 0.696 b/m 0.825 a3 -0.059 h/m 5.25 a4 0.271 L/m 1000 a5 0.027 z0/m 0 a6 -0.010 z1/m 487.375 a7 0.026 
下载: 导出CSV
表 2 ATLAS Ⅲ试验加热和冷却信息
Table 2 Heating and cooling information of ATLAS Ⅲ heating tests
步骤 阶段 热功率/W 日期 天数序号 持续/d 1 加热 0→400 2007-04-02→2007-04-05 0→4 4 2 稳定 400 2007-04-06→2007-05-20 4→49 45 3 加热 400→900 2007-05-21→2007-05-25 49→54 4 4 稳定 900 2007-05-26→2007-07-29 54→120 66 5 加热 900→1400 2007-07-30→2007-08-03 120→125 5 6 稳定 1400 2007-08-04→2008-04-16 125→381 256 7 冷却 0 2008-04-17→2009-11-02 381→945 564
下载: 导出CSV
表 3 参数取值
Table 3 Values of parameters
参数 量值 参数 量值 λG/(W·(K·m)-1) 16.3 ρBoom/(kg·m-3) 2000 λBoom/(W·(K·m)-1) 1.53 a/m 0.16 cG/(J·(K·kg)-1) 500 b/m 0.19 cBoom/(J·(K·kg)-1) 1434 t0/℃ 16 ρG/(kg·m-3) 7850 注：表中下标G和Boom分别表示不锈钢钢管和Boom黏土。
下载: 导出CSV
-
[1] 王驹, 陈伟明, 苏锐, 等. 高放废物地质处置及其若干关键科学问题[J]. 岩石力学与工程学报, 2006, 25(4): 801-812. doi: 10.3321/j.issn:1000-6915.2006.04.015 WANG Ju, CHEN Weiming, SU Rui, et al. Geological disposal of high-level radioactive waste and its key scientific issues[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(4): 801-812. (in Chinese) doi: 10.3321/j.issn:1000-6915.2006.04.015
[2] 叶为民, 王琼, 潘虹, 等. 高压实高庙子膨润土的热传导性能[J]. 岩土工程学报, 2010, 32(6): 821-826. http://cge.nhri.cn/cn/article/id/13418 YE Weimin, WANG Qiong, PAN Hong, et al. Thermal conductivity of compacted GMZ01 bentonite[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(6): 821-826. (in Chinese) http://cge.nhri.cn/cn/article/id/13418
[3] 陈航, 张虎元, 郭永强, 等. 混合型缓冲回填材料导热性能测试与预测研究[J]. 岩石力学与工程学报, 2014, 33(增刊2): 4312-4320. CHEN Hang, ZHANG Huyuan, GUO Yongqiang, et al. Measurement and prediction of thermal properties of bentonite-sand mixtures as buffer backfilling materials for high level radioactive waste[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S2): 4312-4320. (in Chinese)
[4] 谈云志, 李辉, 王培荣, 等. 膨润土受热作用后的水-力性能研究[J]. 岩土力学, 2019, 40(2): 489-496. TAN Yunzhi, LI Hui, WANG Peirong, et al. Hydro-mechanical performances of bentonite respond to heat-treated history[J]. Rock and Soil Mechanics, 2019, 40(2): 489-496. (in Chinese)
[5] ZHANG J R, SUN D A, YU H H, et al. Swelling of unsaturated GMZ07 bentonite at different temperatures[J]. Bulletin of Engineering Geology and the Environment, 2020, 79(2): 959-969. doi: 10.1007/s10064-019-01595-y
[6] 陈正汉, 秦冰. 缓冲/回填材料的热-水-力耦合特性及其应用[M]. 北京: 科学出版社, 2017. CHEN Zhenghan, QIN Bing. Thermal-Hydraulic-Mechanical Coupling Characteristics of Buffer/Backfill Materials and its Application[M]. Beijing: Science Press, 2017. (in Chinese)
[7] 陈正汉, 郭楠. 非饱和土与特殊土力学及工程应用研究的新进展[J]. 岩土力学, 2019, 40(1): 1-54. CHEN Zhenghan, GUO Nan. New developments of mechanics and application for unsaturated soils and special soils[J]. Rock and Soil Mechanics, 2019, 40(1): 1-54. (in Chinese)
[8] 陈正汉. 非饱和土与特殊土力学[M]. 北京: 中国建筑工业出版社, 2022. CHEN Zhenghan. Mechanics for Unsaturated and Special Soils[M]. Beijing: China Architecture & Building Press, 2022. (in Chinese)
[9] CHO W J, KIM G Y. Reconsideration of thermal criteria for Korean spent fuel repository[J]. Annals of Nuclear Energy, 2016, 88: 73-82. doi: 10.1016/j.anucene.2015.09.012
[10] 刘月妙, 王驹, 蔡美峰, 等. 热-力耦合条件下高放废物处置室间距研究[J]. 铀矿地质, 2009, 25(6): 373-379. doi: 10.3969/j.issn.1000-0658.2009.06.009 LIU Yuemiao, WANG Ju, CAI Meifeng, et al. Study on disposal pit space for high-level radioactive waste in thermal-mechanical coupling conditiion[J]. Uranium Geology, 2009, 25(6): 373-379. (in Chinese) doi: 10.3969/j.issn.1000-0658.2009.06.009
[11] CARSLAW H S, JAEGER J C. Conduction of Heat in Solids[M]. 2d ed. Oxford: Clarendon Press, 1959.
[12] CLAESSON J, PROBERT T. Thermoelastic stress due to a rectangular heat source in a semi-infinite medium. Presentation of an analytical solution[J]. Engineering Geology, 1998, 49(3/4): 223-229.
[13] HÖKMARK H, FÄLTH B. Thermal dimensioning of the deep repository[R]. Stockholm: Svensk Kärnbränslehantering AB, 2003.
[14] HÖKMARK H, LÖNNQVIST M, KRISTENSSON O. Strategy for thermal dimensioning of the final repository for spent nuclear fuel[R]. Stockholm: Svensk Kärnbränslehantering AB, 2009.
[15] 刘东东, 项彦勇. 高放射核废处置库温度场的分布线热源解析模型[J]. 岩石力学与工程学报, 2019, 38(增刊1): 2816-2822. LIU Dongdong, XIANG Yanyong. A distributed line heat-source analytical model for the temperature field of a high level nuclear waste repository[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(S1): 2816-2822. (in Chinese)
[16] IKONEN K. Thermal analysis of repository for spent EPR-type fuel[R]. Olkiluoto: Posiva Oy, 2005.
[17] HARTLEY L, HOCH A, JACKSON P, et al. Groundwater flow and transport modeling during the temperate period for the SR-Can assessment. Forsmark area-version 1.2[R]. Stockholm: Svensk Kärnbränslehantering AB, 2006.
[18] 陈永贵, 贾灵艳, 叶为民, 等. 施工接缝对缓冲材料水-力特性影响研究进展[J]. 岩土工程学报, 2017, 39(1): 138-147. doi: 10.11779/CJGE201701012 CHEN Yonggui, JIA Lingyan, YE Weimin, et al. Advances in hydro-mechanical behaviors of buffer materials under effect of technological gaps[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(1): 138-147. (in Chinese) doi: 10.11779/CJGE201701012
[19] XU X, HE L Q, SUN D A, et al. Validation of the fully-analytical solution with temperature superposition for the nuclear waste repository[J]. Nuclear Engineering and Design, 2023, 409: 112367. doi: 10.1016/j.nucengdes.2023.112367
[20] 陶文铨. 传热学[M]. 5版. 北京: 高等教育出版社, 2019. TAO Wenquan. Heat Transfer[M]. 5th ed. Beijing: Higher Education Press, 2019. (in Chinese)
[21] HE L Q, ZHOU X Y, SUN D A. Fully analytical solution in time and space domains on temperature in multi-barrier nuclear waste repository[J]. Computers and Geotechnics, 2023, 154: 105164. doi: 10.1016/j.compgeo.2022.105164
[22] CRUMP K S. Numerical inversion of Laplace transforms using a Fourier series approximation[J]. Journal of the ACM, 1976, 23(1): 89-96. doi: 10.1145/321921.321931
[23] 徐国庆. 国际高放废物处置研发工作在花岗岩地区的进展[J]. 世界核地质科学, 2016, 33(2): 119-124. doi: 10.3969/j.issn.1672-0636.2016.02.010 XU Guoqing. Abroad progress in R & D work on high-level radioactive waste disposal in granite areas[J]. World Nuclear Geoscience, 2016, 33(2): 119-124. (in Chinese) doi: 10.3969/j.issn.1672-0636.2016.02.010
[24] CHEN G J, SILLEN X, VERSTRICHT J, et al. ATLAS Ⅲ in situ heating test in boom clay: field data, observation and interpretation[J]. Computers and Geotechnics, 2011, 38(5): 683-696. doi: 10.1016/j.compgeo.2011.04.001
[25] DE BRUYN D, LABAT S. The second phase of ATLAS: the continuation of a running THM test in the HADES underground research facility at Mol[J]. Engineering Geology, 2002, 64(2/3): 309-316.
-
其他相关附件
计量
- 文章访问数: 306
- HTML全文浏览量: 46
- PDF下载量: 81
