One-dimensional coupled model for large-deformation electroosmotic consolidation and heavy metal ion migration of silt
-
摘要: 电动处理技术在淤泥高效脱水和重金属污染物同步去除方面潜力巨大。基于分段线性差分法,将电渗固结和溶质迁移理论相结合,建立了一维大变形电渗固结-重金属离子迁移耦合模型ECT1。该模型可以模拟电场、渗流场和化学场的耦合作用,可以考虑电动处理过程中土性参数的非线性变化,以及吸附/解吸、中和/电离、沉淀/溶解等多种化学反应。分别采用Alshawabkeh污染土电动修复试验及数值模拟、高含水率重金属污染土电动处理试验,验证分析了ECT1模型的可靠性。通过算例分析,进一步阐述了高含水率污染淤泥电动处理过程中重金属迁移与电渗固结的耦合机理。Abstract: Dewatering and metal removal of silt can be obtained by applying electrokinetics. A coupled model for one-dimensional electroosmosis consolidation and ion migration, called ECT1, is proposed. It employs the piecewise linear finite difference method to simulate the ion migration and large-strain consolidation of soils under the coupling action of electric field, seepage field and chemical field, and to account for the nonlinear variation of physical and electro-chemical properties of soils, such as adsorption/desorption, neutralization/ionization, and precipitation/dissolution. The model is verified by Alshawabkeh's electrokinetic remediation experimental results for contaminated soil, numerical solutions and treatment tests on high-moisture content heavy metal-contaminated soil electrokinetic. The proposed model is applied to example studies to elucidate the coupling mechanism of heavy metal migration and consolidation in the electrokinetic treatment of high-moisture content contaminated silt.
-
Keywords:
- electroosmotic consolidation /
- ion migration /
- heavy metal /
- silt /
- electrokinetic treatment
-
0. 引言
随着中国工业化、城镇化发展进程的加快,港口航道泥沙淤积、河湖库底泥淤积现象严重[1]。淤积底泥大大降低了港口航道通航和调蓄洪水的能力,造成河湖库水质污染、水体恶化,对河道及周边环境造成污染[2]。为保证良好的水质及便利的通航条件,中国每年都清淤出数10亿m3以上的疏浚淤泥[3]。疏浚出的淤泥主要为细颗粒黏土,具有含水率高、不排水强度低、渗透性小等性质,难以在工程上直接利用,一般需要设置淤泥堆场[4]。由于高含水率和高黏粒含量,泥浆在堆场内天然自重落淤固结缓慢,导致长期占用大量土地资源,还易造成周边环境受到污染[5]。
为了加快土颗粒与水的分离,通过加入絮凝剂使土颗粒凝聚快速沉降[6]。章荣军等[7]开展室内絮凝-堆载预压试验,结果表明加入絮凝剂能够增大预压后淤泥的结构屈服应力。无机絮凝剂由于成本低,絮凝效果较好,已广泛应用于淤泥脱水处理中。袁国辉等[8]采用无机絮凝剂CaCl2优化传统电渗法,提高了淤泥后期的排水效率。刘飞禹等[9]利用无机絮凝剂FeCl3和Al2(SO4)3对电渗法处理河道疏浚淤泥进行预处理,明显改善电渗排水效果。上述研究大多集中在絮凝剂对疏浚淤泥沉降性能及对后续脱水工艺影响的研究,针对絮凝脱水后的水质与淤泥沉降性能之间的关联少有涉及,上覆水水质也是衡量絮凝药剂脱水性能的一项重要指标。在实际工程中,上覆水水质达到一定的排放标准才可以进行排放。
鉴于此,本文选取5种典型的无机絮凝剂对疏浚淤泥进行调理,开展一系列量筒试验,探讨絮凝剂种类及掺量对淤泥沉积特性及出水水质的影响,并初步确定最优絮凝剂及其最优掺量。通对过自重沉积固结稳定后沉降量、沉降速率、沉积稳定时间及上覆水浊度值进行分析,筛选出高效脱水且保证尾水排放达标的絮凝剂,以期为吹填淤泥快速落淤过程中的药剂选择和余水处理提供借鉴参考。
1. 试验材料与方案
1.1 试验材料
本文所用的疏浚淤泥取自西太湖宜兴周铁堆场,经网筛去除较大砾石和垃圾。通过一系列室内试验,测得淤泥的主要物理性质见表 1,粒径分布曲线见图 1。根据细粒土的分类标准,该淤泥属于高液限黏土(CH)。试验初步选用无机絮凝剂CaCl2、Ca(OH)2,PAC(聚合氯化铝),PAS(聚合硫酸铝),PAFSI(聚合硅酸铝铁)。
表 1 淤泥的基本物理性质Table 1. Basic physical properties of dredged materials含水率/% 液限/% 塑限/% 塑性指数 黏粒含量①/% 相对质量密度 有机质含量②/% 100.4 55.9 20.6 35.3 28.00 2.72 0.75 注:①粒径小于5 μm;②重铬酸钾容量法。 1.2 试验方案
试验前先将现场取回的淤泥过筛去除淤泥中较大的石子和其它杂质,配制成含水率为220%(4倍液限)的试验淤泥。采用下式计算出原状淤泥中需要添加蒸馏水的质量mw:
mw=w1−w01+w0m1。 (1) 式中:w0为原状淤泥含水率;w1为试验淤泥初始含水率;m1为原状淤泥的质量。加水混合后淤泥放在六连搅拌器上持续搅拌直至均匀。
为初步明确絮凝剂对高含水率疏浚淤泥自重沉积特性的影响,5种絮凝剂统一配制成100 mL目标浓度为1%的絮凝剂溶液。试验选用规格为500 mL的玻璃量筒。取400 mL的混合液倒入量筒,倒入过程中用玻璃棒进行引流,结束后用保鲜膜密闭封顶,以防水分散失。待泥水分离界面高度不再有明显变化时停止观察,记录不同时刻泥水分界面高度变化的相关试验数据,试验结束后检测余水浊度值。
根据前期沉降试验结果,得出CaCl2,PAS在沉积稳定时沉降量相对较大,分别降至187,198 mL。根据试验结果,初步筛选出CaCl2、PAS为最优絮凝剂,改变其掺量,具体试验方案见表 2。
表 2 量筒沉积试验方案Table 2. Schemes of measuring cylinder sedimentation tests加药种类 絮凝剂掺量/% CaCl2 0,0.5,0.8,0.9,1.0,1.2 PAS 0,0.1,0.2,0.5,0.7,1.0 2. 试验结果与分析
2.1 泥面沉降曲线
图 2为5种絮凝剂在目标浓度下泥水分离界面随时间变化的沉降曲线。在絮凝前期,PAFSI和PAC展现出良好的脱水性能,促使淤泥沉降效果较明显。直至72 h后,添加CaCl2和PAS的泥浆混合液沉降效果较好,最终沉降量较低。添加Ca(OH)2的泥浆混合液在前期絮凝沉降效果较好,在72 h后沉降基本不发生变化。根据初步沉降试验结果,以下试验选定最优絮凝剂为CaCl2和PAS,改变其掺量。
图 3为不同CaCl2掺量下泥水分界面随时间变化的沉降曲线。当CaCl2掺量为0.8%时,泥面沉降体积值达到最低,显著改善了疏浚淤泥的脱水性能。最终淤泥沉降体积值随CaCl2掺量的增加,呈先降低再上升的变化趋势。这是由于Ca2+通过离子交换和电中和作用,压缩黏土颗粒表面双电层厚度,使得颗粒间斥力降低,颗粒更易聚集沉降。随着CaCl2掺量的增加,黏土颗粒表面的负电荷全部中和,此时絮凝效果达到最佳。随着掺量进一步增加,黏土颗粒表面带正电荷,颗粒间的斥力增大,絮凝效果开始减弱。
图 4为不同PAS掺量下泥水分离界面沉降曲线。可以看出,当PAS掺量在0.2%时,泥面沉降量最低。随着PAS掺量的增加,土颗粒相互碰撞次数增多,发生吸附架桥,凝聚沉降作用。当PAS掺量达到一定时,颗粒间的排斥能为零,絮凝效果最好,最终沉降量较大。当PAS掺量的进一步增大时,水中存在较多聚合铝离子带正电,完全附着在黏土颗粒表面使得其带上正电荷,颗粒之间的斥力增加,絮凝效果开始减弱。
2.2 泥水分界面沉降速率
泥面沉降速率按照下式进行计算:
v=Si+1−Siti+1−ti。 (2) 式中:Si+1,Si分别为ti+1,ti时刻泥水分离界面沉降量。
对不同CaCl2,PAS掺量下淤泥沉降速率进行计算,结果如图 5所示。在不同絮凝剂掺量条件下,沉积速率的变化大致分为3个阶段。①第一阶段在泥水分界面刚出现时沉降速率较大,曲线呈快速下降趋势,淤泥沉降速率较大。②第二阶段细小的土颗粒絮凝成团胶结成大的聚体,絮凝团之间相互阻碍开始沉降,沉降速率快速降低直至一个稳定区间范围内。③第三阶段沉降速率变化极为平缓,随着时间的推移,淤泥自重沉降固结趋于稳定时,沉降速率趋近于零。由于在水下受到浮力作用,固结变形仍在缓慢进行。
2.3 沉降稳定时间Tc
疏浚淤泥在沉降过程中,泥水分离界面沉降速率不是固定的。参考Imai对泥浆自重沉降稳定时间Tc的定义[10],根据沉降过程中泥水分离界面随时间在半对数坐标系下的沉降变化曲线,计算不同絮凝剂掺量下泥浆自重沉降稳定时间Tc。泥浆沉降稳定时间确定方法示意图如图 6所示,具体结果见表 3。
表 3 不同絮凝剂掺量下淤泥自重沉积稳定时间Table 3. Stabilization time of self-weight sedimentation of sludge under different flocculant dosages絮凝剂 CaCl2/PAS 掺量/% 0 0.5/0.1 0.8/0.2 0.9/0.5 1.0/0.7 1.2/1.0 Tc/d 3.39 2.63/2.65 2.54/2.58 2.58/2.60 2.77/2.61 2.80/2.63 在添加0.8%CaCl2和0.2%PAS掺量下,混合液泥浆能较快地达到自重沉积稳定状态,沉降稳定时间较短,分别为2.54,2.6 d。
2.4 出水水质
对加入5种典型无机絮凝剂的混合泥浆液7 d沉降稳定后余水浊度值进行检测,结果如图 7所示。可以发现加入絮凝剂可大大降低余水的浊度值。CaCl2,PAS的去浊效果相对较优,且能够降低淤泥pH值。
在初步选定最优絮凝剂CaCl2,PAS的基础上,对不同掺量下混合泥浆自重沉降稳定后上覆水的浊度值进行测定。
结合图 8可发现,絮凝剂PAS在掺量为0.2%时,去浊率较高;絮凝剂CaCl2掺量为0.8%时可达到相同的去浊效果。达到相同的去浊效果时,PAS掺量较少。当CaCl2的掺量为0.8%时,去除率达到83.18%。当PAS的掺量为0.2%时,浊度去除率达到84%。这可能是由于PAS具有较大的比表面积,水解后产生较多的带正电粒子,能够与淤泥颗粒表面带负电粒子充分接触发生电性中和,使得颗粒失去稳定状态进而沉降下沉。
3. 结论
(1)通过5种典型无机絮凝剂的絮凝沉降试验,添加无机絮凝剂可大大加快高含水率疏浚淤泥的沉降速率,最终沉降量较低且余水相对清澈。对5种絮凝剂在目标浓度下泥浆自重沉积稳定后的沉降量进行观察,PAFSI,PAC和Ca(OH)2前期絮凝较快,CaCl2,PAS后期沉降较快,最终泥水分离界面沉降值较低。
(2)CaCl2,PAS表现出良好的脱水性能,CaCl2主要通过离子交换作用和电性中和发挥絮凝作用;PAS主要通过中和土颗粒表面负电荷发挥絮凝作用。CaCl2,PAS的最优掺量为0.8%,0.2%。在最优掺量下,浊度去除率分别为83.2%,84.5%,沉降稳定时间缩短至2.54,2.57 d。
-
表 1 ECT1排水边界条件
Table 1 Drainage boundary conditions of ECT1
条件 阳极 阴极 不透水 ite,b=ith,b=0 ite,1=ith,1=0 透水 ite,b=2(Vtm−VtRj)l0ith,b=2(hw2−htt,Rj)l0 ite,1=2Vt1l0ith,1=2(htt,1−hw1)l0 注:ite,b和ith,b为阳极边界处的电势梯度和水力梯度,ite,1和ith,1为阴极边界的电势梯度和水力梯度。 表 2 ECT1离子通量边界条件
Table 2 Species flux boundary conditions of ECT1
离子 阳极 阴极 Mem+ JtMem + ,b=ctMem + ,bJtw,b JtMem + ,1=ctMem + ,1Jtw,1 H+ JtH+,b=ctH+,bJtw,b+RtH+HtRj JtH + ,1=ctH + ,1Jtw,1 OH- JtOH−,b=ctOH−,bJtw,b JtOH−,1=ctOH−,1Jtw,1−RtOH−Ht1 注:下标b代表阳极边界,下标1代表阴极边界。以阳极边界为例,JtMem + ,b和ctMem + ,b分别为重金属离子Mem+在阳极边界(阳极电解室)处的通量和浓度,当阳极边界封闭或无该离子补充时,JtMem + ,b=0。 表 3 数值模拟输入参数[15]
Table 3 Input parameters of numerical simulation
参数 取值 水力渗透系数kh/(m·s-1) 6.79×10-10 电渗系数ke, 0/(m2·s-1·V-1) 1×10-9 压缩系数av/(MPa-1) 2.27×10-2 初始孔隙比e0 1.27 干重度γd/(kg·m-3) 1800 恒定电流密度I/(A·m-2) 1.33 模型几何边长L/m 0.7 初始铅浓度cPb,0/(mg·kg-1) 5320 初始酸碱性pH0 4 土样阳离子交换能力CEC/(mEq·100g-1) 1.06 土样曲折系数τ 0.44 表 4 扩散系数与离子迁移率
Table 4 Diffusion coefficient and ion mobility of species
物质种类 扩散系数
Di/(m2·s-1)离子电迁移率
ui/(m2·s-1·V-1)Pb2 + 9.45×10-10 7.36×10-8 H + 93.1×10-10 36.25×10-8 OH− 52.7×10-10 20.58×10-8 NO−3 19.0×10-10 7.44×10-8 表 5 基本物化参数
Table 5 Basic physical and chemical parameters of kaolin sample
物化参数 取值 初始含水率 79.2% 相对密度Gs 2.65 初始孔隙比e0 2.1 初始水力渗透系数kh,0/(m·s-1) 2.1×10-9 水力渗透指数Ck 0.9 初始电渗系数ke,0/(m2·sV-1) 5.8×10-9 初始酸碱性pH 8.32 压缩指数Cc 0.4 初始铅浓度cPb,0/(g·kg-1) 1 阳离子交换能力CEC/(mEq·100g-1) 2.07 曲折系数τ 0.45 表 6 算例输入参数
Table 6 Input parameters of example
参数 取值 初始水力渗透系数kh,0/(m·s-1) 2.08×10-9 水力渗透指数Ck 0.99 初始电渗系数ke,0/(m2·sV-1) 2.35×10-9 压缩指数Cc 0.5 初始孔隙比e0 1.9 初始孔隙率n0 0.66 初始相对密度Gs 2.62 阳极和阴极间电势差Vm/V 30 模型几何边长L×W×H(m×m×m) 1×1×1 初始铅浓度cPb,0/(g·kg-1) 1.9 初始酸碱性pH 6.6 阳离子交换能力CEC/(mEq·100g-1) 1.7 曲折系数τ 0.4 集总系数ε 0.19 -
[1] 周建, 魏利闯, 詹芳蕾, 等. 生物表面活性剂及其与柠檬酸联合用于污泥重金属电动修复[J]. 湖南大学学报(自然科学版), 2019, 46(6): 109–119. https://www.cnki.com.cn/Article/CJFDTOTAL-HNDX201906016.htm ZHOU Jian, WEI Li-chuang, ZHAN Fang-lei, et al. Electrokinetic repair of heavy metals in sludge by biosurfactant and its combination with citric acid[J]. Journal of Hunan University (Natural Sciences), 2019, 46(6): 109–119. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HNDX201906016.htm
[2] TUAN P A, SILLANPÄÄ M. Migration of ions and organic matter during electro-dewatering of anaerobic sludge[J]. Journal of Hazardous Materials, 2010, 173(1/2/3): 54–61.
[3] TANG X Q, LI Q Y, WANG Z H, et al. Improved isolation of cadmium from paddy soil by novel technology based on pore water drainage with graphite-contained electro-kinetic geosynthetics[J]. Environmental Science and Pollution Research International, 2018, 25(14): 14244–14253. doi: 10.1007/s11356-018-1664-4
[4] ESRIG M I. Pore pressure, consolidation and electrokinetics[J]. Journal of the SMFD, ASCE, 1968, 94(SM4): 899–921.
[5] HU L, WU W, WU H. Numerical model of electro-osmotic consolidation in clay[J]. Géotechnique, 2012, 62(6): 537–541. doi: 10.1680/geot.11.T.008
[6] 吴辉, 胡黎明. 考虑电导率变化的电渗固结模型[J]. 岩土工程学报, 2013, 35(4): 734–738. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201304020.htm WU Hui, HU Li-ming. Numerical simulation of electro-osmosis consolidation considering variation of electrical conductivity[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(4): 734–738. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201304020.htm
[7] 王军, 符洪涛, 蔡袁强, 等. 线性堆载下软黏土一维电渗固结理论与试验分析[J]. 岩石力学与工程学报, 2014, 33(1): 179–188. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201401021.htm WANG Jun, FU Hong-tao, CAI Yuan-qiang, et al. Analyses of one-dimensional electroosmotic consolidation theory and test of soft clay under linear load[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(1): 179–188. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201401021.htm
[8] FELDKAMP J R, BELHOMME G M. Large-strain electrokinetic consolidation: theory and experiment in one dimension[J]. Géotechnique, 1990, 40(4): 557–568. doi: 10.1680/geot.1990.40.4.557
[9] 王柳江, 刘斯宏, 王子健, 等. 堆载–电渗联合作用下的一维非线性大变形固结理论[J]. 工程力学, 2013, 30(12): 91–98. doi: 10.6052/j.issn.1000-4750.2012.04.0303 WANG Liu-jiang, LIU Si-hong, WANG Zi-jian, et al. A consolidation theory for one-dimensional large deformation problems under combined action of load and electroosmosis[J]. Engineering Mechanics, 2013, 30(12): 91–98. (in Chinese) doi: 10.6052/j.issn.1000-4750.2012.04.0303
[10] YUAN J, HICKS M A. Large deformation elastic electro-osmosis consolidation of clays[J]. Computers and Geotechnics, 2013, 54: 60–68. doi: 10.1016/j.compgeo.2013.05.012
[11] YUAN J, HICKS M A. Numerical simulation of elasto-plastic electro-osmosis consolidation at large strain[J]. Acta Geotechnica, 2016, 11(1): 127–143. doi: 10.1007/s11440-015-0366-z
[12] ZHOU Y D, DENG A, WANG C. Finite-difference model for one-dimensional electro-osmotic consolidation[J]. Computers and Geotechnics, 2013, 54: 152–165. doi: 10.1016/j.compgeo.2013.06.003
[13] 冯源. 城市污水污泥电动脱水机理试验研究及多场耦合作用理论分析[D]. 杭州: 浙江大学, 2012. FENG Yuan. Experimental Study on Electrokinetic Dewatering Mechanism of Sewage Sludge and Theoretical Analyses of Multi-Field Coupled Phenomenon[D]. Hangzhou: Zhejiang University, 2012. (in Chinese)
[14] CORAPCIOGLU M Y. Formulation of electro- chemicoosmotic processes in soils[J]. Transport in Porous Media, 1991, 6(4): 435–444.
[15] ALSHAWABKEH A N, ACAR Y B. Removal of contaminants from soils by electrokinetics: a theoretical treatise[J]. Journal of Environmental Science and Health Part A: Environmental Science and Engineering and Toxicology, 1992, 27(7): 1835–1861. doi: 10.1080/10934529209375828
[16] ALSHAWABKEH A N, ACAR Y B. Electrokinetic remediation: II theoretical model[J]. Journal of Geotechnical Engineering, 1996, 122(3): 186–196. doi: 10.1061/(ASCE)0733-9410(1996)122:3(186)
[17] YEUNG A T, DATLA S. Fundamental formulation of electrokinetic extraction of contaminants from soil[J]. Canadian Geotechnical Journal, 1995, 32(4): 569–583. doi: 10.1139/t95-060
[18] AL-HAMDAN A Z, REDDY K R. Electrokinetic remediation modeling incorporating geochemical effects[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2008, 134(1): 91–105. doi: 10.1061/(ASCE)1090-0241(2008)134:1(91)
[19] MASI M, CECCARINI A, IANNELLI R. Multi species reactive transport modelling of electrokinetic remediation of harbour sediments[J]. Journal of Hazardous Materials, 2017, 326: 187–196. doi: 10.1016/j.jhazmat.2016.12.032
[20] GOODISMAN J. Electrochemistry: Theoretical Foundations, Quantum and Statistical Mechanics, Thermodynamics, the Solid State[M]. New York: Wiley, 1987.
[21] YEUNG A T, HSU C N, MENON R M. Electrokinetic extraction of lead from kaolinites: I numerical modeling[J]. The Environmentalist, 2011, 31(1): 26–32. doi: 10.1007/s10669-010-9295-4
[22] KIM S O, KIM J J, KIM K W, et al. Models and experiments on electrokinetic removal of Pb(II) from kaolinite clay[J]. Separation Science and Technology, 2005, 39(8): 1927–1951. doi: 10.1081/SS-120030775
[23] ACAR Y B, ALSHAWABKEH A N. Electrokinetic remediation I: pilot-scale tests with lead-spiked kaolinite[J]. Journal of Geotechnical Engineering, 1996, 122(3): 173–185. doi: 10.1061/(ASCE)0733-9410(1996)122:3(173)
[24] YONG R N, WARKENTIN B P, PHADUNGCHEWIT Y, et al. Buffer capacity and lead retention in some clay materials[J]. Water, Air, and Soil Pollution, 1990, 53(1/2): 53–67.