• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

非饱和土空间分数阶渗流模型的有限差分方法研究

王睿, 周宏伟, 卓壮, 薛东杰, 杨帅

王睿, 周宏伟, 卓壮, 薛东杰, 杨帅. 非饱和土空间分数阶渗流模型的有限差分方法研究[J]. 岩土工程学报, 2020, 42(9): 1759-1764. DOI: 10.11779/CJGE202009021
引用本文: 王睿, 周宏伟, 卓壮, 薛东杰, 杨帅. 非饱和土空间分数阶渗流模型的有限差分方法研究[J]. 岩土工程学报, 2020, 42(9): 1759-1764. DOI: 10.11779/CJGE202009021
WANG Rui, ZHOU Hong-wei, ZHUO Zhuang, XUE Dong-jie, YANG Shuai. Finite difference method for space-fractional seepage process in unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(9): 1759-1764. DOI: 10.11779/CJGE202009021
Citation: WANG Rui, ZHOU Hong-wei, ZHUO Zhuang, XUE Dong-jie, YANG Shuai. Finite difference method for space-fractional seepage process in unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(9): 1759-1764. DOI: 10.11779/CJGE202009021

非饱和土空间分数阶渗流模型的有限差分方法研究  English Version

基金项目: 

国家自然科学基金项目 51674266

“十三五”国家重点研发计划项目 2016YFC0600704

详细信息
    作者简介:

    王睿(1983—),男,博士,主要从事岩土力学、渗流力学等方面的科研工作。E-mail:wangruicumtb@163.com

    通讯作者:

    周宏伟, E-mail:zhw@cumtb.edu.cn

  • 中图分类号: TU441

Finite difference method for space-fractional seepage process in unsaturated soil

  • 摘要: 地下水在非饱和土体中的输运过程及渗流特性具有重要的研究意义。首先,从反常扩散方程的角度分析地下水在非饱和土中的输运过程,引入Conformable导数,得到了一维情况下非饱和土空间分数阶渗流方程;其次,用全隐式形式的有限差分法求解渗流方程离散格式,得到了求解的迭代矩阵;最后,根据已有论文的试验数据对渗流模型的有效性进行了验证,并对公式中的参数进行了敏感性分析。结果表明,新模型能够较好地描述地下水在非饱和土体中的输运过程,Conformable导数具有较好的适用性。
    Abstract: The transport process and seepage properties of underground water in unsaturated soil are of great research significance. Firstly, the transport process of underground water in unsaturated soil in an anomalous diffusion perspective is analyzed, and by involving the Conformable derivative, the space-fractional anomalous diffusion equation in one-dimensional situation is derived. Secondly, the discretization schemes of this seepage equation are derived by using the full-implicit finite difference method. The iterative matrix of numerical modelling is derived as well. In the end, the effectiveness of the numerical solution in the space-fractional anomalous diffusion equation is verified by comparing the numerical solution with the experimental data in the existing paper. In addition, the sensitivity of the parameters in the seepage equation is analyzed. The results show that the new model can well characterize the transport process of underground water in unsaturated soil, and the Conformable derivative is suitable for this method.
  • 自20世纪七八十年代以来,竖向隔离屏障被广泛地应用到坝体、垃圾填埋场和城市工业污染场地的工程防渗措施之中[1-2]。目前中国大多采用水泥深层搅拌桩为主的竖向防渗墙阻隔填埋场渗滤液或污染场地内的污染地下水迁移。在填埋场和污染场地使用的防渗墙中,水泥-膨润土(CB)作为竖向阻隔屏障材料在英法等欧洲国家在工程中使用较为普遍[3],美国则更广泛使用土-膨润土系(SB)竖向阻隔屏障[4]。对比上述3种常见的竖向阻隔墙,CB阻隔墙属于刚性材料,具有强度。在实际工程中能够有效地阻隔场地污染物的向外迁移,但在阻隔污染物的过程中容易受到重金属或有机污染物对水泥的水化抑制以及硫酸盐的腐蚀作用。SB隔离屏障虽然造价较低,具有良好的防渗性能但由于较低强度的特性导致其不能适用于工程屏障强度有要求的场景[5]。基于上述隔离墙的特点,国内外学者提出了一种使用粒化高炉矿粉(GGBS)替代约80%水泥作为竖向阻隔屏障材料。与水泥相比,GGBS价格适中,且是一种低碳、绿色可持续材料[6-7]。Haha等[8]和Yi等[9]研究发现GGBS水化速度较慢,在活性MgO的激发作用下可以明显提高其水化速度,且采用MgO激发矿粉加固土,可显著提高土体的力学性能。作者所在课题组[10-13]将活性氧化镁碱激发矿粉用于竖向阻隔屏障制备了氧化镁碱激发矿粉-膨润土-土(MSBS)竖向阻隔屏障,并在其防渗性能、强度特征、干湿循环耐久性、微观特征等方面做了研究。

    近些年来,工程屏障(膨润土防水毯(GCL)、压实黏土衬垫、竖向阻隔屏障)的化学渗透膜效应在竖向隔离屏障阻滞污染物溶质运移方面受到了广泛的关注,笔者所在课题组研究了聚阴离子纤维素(PAC)和六偏磷酸钠(SHMP)改性膨润土/砂竖向阻隔屏障材料的化学渗透膜效应和重金属Pb2+的扩散特征,结果表明,改性的屏障材料针对Pb2+的有效扩散系数低于未改性的屏障材料[14-17]。化学渗透指的是由于屏障材料呈现半透膜效应,溶剂沿溶质浓度梯度相反方向迁移的现象,也即从溶质低浓度位置迁移至溶质高浓度位置的现象,进而导致溶质在屏障的浓度扩散运移受到约束[17]。常用化学渗透膜效率系数ω来评价屏障材料对溶质的阻滞能力。ω位于0~1,0,1分别表示材料不具备化学渗透膜效应、具有理想化学渗透膜效应[18]。研究表明[19]指出化学渗透膜效应对于评价屏障阻隔污染物溶质运移具有重要的意义。

    目前关于MSBS阻隔材料的化学渗透膜效应和溶质扩散运移特征报道很少。基于此,研究新型MgO激发高炉矿粉-膨润土竖向隔离墙的膜效应和溶质运移扩散特征,并对比分析与其他工程屏障材料的差异特征则很有必要。本文主要采用MgO作为碱激发剂用以激发粒化高炉矿渣制备成氧化镁碱激发矿粉-膨润土竖向阻隔防渗墙(MSBS),采用土柱化学渗透试验,研究在重金属Pb2+离子作用下MSBS阻隔材料的膜效应特性和溶质扩散运移规律。

    MSBS阻隔材料是由碱激发粒化高炉矿粉、砂土和膨润土按照一定比例混合而成。试验用土采用南京河漫滩砂土,膨润土为商用的钠化改良膨润土,产地江苏镇江。其中砂土过1.0 mm筛,膨润土过0.075 mm筛,两者基本物理指标见表 1所示。GGBS产自南京某水泥厂,活性氧化镁为上海某化工厂生产,两者主要的物理化学指标见表 2所示。

    表  1  试验用土的基本土工参数
    Table  1.  Properties of soil used in this study
    基本指标 天然含水率/% 重度Gs 塑限wP/% 液限wL/% 黏粒含量/% 粉粒含量/% 砂粒含量/% 比表面积/(m2·g-1) pH 蒙脱石含量/%
    南京砂土 4.81 2.62 5.62 14.18 80.2 7.32
    镇江膨润土 11.20 2.66 54 103 99.00 1.00 378.5 8.60 66.9
    下载: 导出CSV 
    | 显示表格
    表  2  粒化高炉矿粉基本物理-化学特征
    Table  2.  Properties of GGBS
    基本指标 数值
    GGBS MgO
    颜色 灰白色粉末 白色粉末状
    含量/% 77.6
    活性/s 90~100
    比表面积/(m2·g-1) 0.2564 28.0230
    密度/(g·mL-1) 3.58
    碱度 1.871
    pH 12.21 10.71
    下载: 导出CSV 
    | 显示表格

    试样制备方法:首先将混合干料按照设计比例使用台式搅拌机均匀拌和5 min,按照Wu等[12]研究建议土、膨润土、GGBS和MgO的质量比设置为100︰10︰9︰1;随后向阻隔材料中加入适量的水并控制浆液的坍落度在150 mm左右,试验过程中保证屏障材料各组分不变,最终屏障材料的含水率为35.4%;之后将屏障材料分层依次填入环刀内,使其高度达到10 mm, 并用密封袋密封置于温度22℃、湿度95%的标准环境中养护90 d,确保屏障材料水化反应基本完成,各力学指标基本稳定[12]

    本次试验采用Pb(NO3)2溶液作为典型重金属污染液。溶液制备方法:根据现有报道的工业污染场地地下水Pb2+污染程度以及《地下水环境质量标准》规定[20-21],本文将0.1 mg/L的Pb2+饮用水标准设置为污染液的初始浓度,之后按照初始浓度的5,10,20,100,500倍配制6种不同浓度的Pb2+离子污染液。根据ASTM D1125[22]和ASTM E70[23]相关规范测得不同污染液的物理化学指标pH、电导率(EC)和浓度,见表 3

    表  3  Pb(NO3)2污染液基本化学性质
    Table  3.  Basic chemical properties of Pb(NO3)2 liquids
    溶液 Pb2+浓度c0/(mg·L-1) EC/(μS·cm-1) pH
    Pb(NO3)2 0.1 44.5 6.84
    0.5 51.6 659
    1 84.8 6.45
    2 120.1 6.31
    10 249.6 5.96
    50 487.4 5.63
    下载: 导出CSV 
    | 显示表格

    化学渗透试验测试装置主要由测试腔室、压差测试系统、数据采集系统和溶液供给和收集系统组成,能够实现试样的有效扩散和膜效率系数的测定,具体示意图见图 1[16]。测试腔室由上下聚四氟乙烯(PTFE)顶部和底座、1个直径61.8 mm,高度20 mm的不锈钢环刀组成。压力传感器(HM50-1-A-F1-W2)为某公司生产的新型防腐蚀压力传感器,量程范围0~100 kPa,精度0.001 kPa。数据采集器为某公司设计生产型号为DataTaker DT80G Series 3。

    图  1  测试装置示意图[16]
    Figure  1.  Schematic diagram of testing apparatus[16]

    将制备的试样放置到80 kPa下的真空腔室中养护24 h直至试样完全饱和。在腔室两端放上透水石和浸湿的滤纸,将制备好的试样装入到四周涂抹凡士林的环刀中放到滤纸上方,在其顶部再放置一层滤纸和透水石,完成整个装样流程。在试样准备阶段完成后开始化学渗透试验,主要步骤包括冲刷阶段、基准压差测试阶段和化学渗透压差测试阶段。①冲刷阶段:在试样底部施加20 kPa压力的条件下注入煮沸的去离子水(DIW),测试收集到的渗滤液的体积和EC值,当EC值小于0.5倍的初始电导率时,可认为已完成本阶段试验。试验过程中注入去离子水的目的是一方面为了去除自身可能残留的金属阳离子防止其对后续试验参数的测定造成扰动,另一方面可以保证试样充分饱和[24]。②基准压差测试阶段:完成第一阶段试验后,以0.0252 mL/min的速度[25]从试样两端注入去离子水同时每分钟测定一次两端的基准压差,测试顶部和底部收集到液体的pH、EC和离子浓度,当两端测得的EC、离子浓度、pH和压差读数变化范围在±1%以内,即认为稳定可开始下一阶段试验。③化学渗透压差测试阶段:将上一阶段试样顶部注入的去离子水依次更换0.1 mg/L(4.8×10-4 mM)、0.5 mg/L(0.00241 mM)、1 mg/L(0.00483 mM)、2 mg/L(0.00965 mM)、10 mg/L(0.04826 mM)、50 mg/L(0.241 mM)等6种不同测试浓度的Pb(NO3)2污染液,底部去离子水的注入速率保持不变。试验过程中定期测试两端收集到渗滤液,待渗滤液离子浓度、pH、EC以及试样两端的压差稳定后,将顶部注入的溶液替换为更高浓度的污染液并重复上一次试验,直到全部浓度的污染液均测试完成后结束化学渗透压差测试。需要注意的是试验中所用掺量配比的MSBS阻隔材料在自来水中的渗透系数为2.86×10-10 m/s,满足设计要求[12]

    化学渗透膜效应通常来用化学渗透膜效率系数ω来评价,该系数主要采用封闭边界条件和无边界条件两种方法测试[26-27]。其中在封闭条件中,试样两端产生的化学渗透压差,可以抵消化学渗透压流动的趋势,膜效应根据压力积累大小和边界上的溶液浓度来确定的。通常来说封闭系统边界条件适用于:①试样恒定浓度差;②试样恒定体积;③零流量液体渗流通过试样[28]。每次试验中保持样品上下两端浓度差恒定和测试试样的体积恒定,通过测量化学渗透压差ΔP来计算试样化学渗透膜效率系数,同时兼测溶质的有效扩散系数D*和阻滞因子Rd

    试样的化学渗透膜效率系数ω主要根据试样两端实际测得化学渗透压差与理论计算所得值的比值[29]

    ω=ΔPΔπ,
    (1)

    式中,∆P为试样两端实际测得的化学渗透压力差(kPa),∆π为试样两端理论化学渗透压力值(kPa),其数值一般使用van’t Hoff方程[24]计算在理想状态下化学渗透压力差。

    Δπ=vRTNi=1Δci,
    (2)

    式中,v为电解质的离子数,在本试验Pb(NO3)2污染液中取3,R为理想气体常数,数值8.3145 J/(mol/K),T为绝对温度,室温25℃下等于298.5 K,∆Ci为两侧溶质浓度差(mol/L),N为溶质离子种类。

    图 2为污染液溶质的扩散过程,可以看出当单位面积累积溶质通量随着时间线性变化时,溶质扩散达到稳定状态。假设为C0t=C0C0b=0(C0tC0b分别为试样顶部和底部的进液浓度),的理想边界条件可使用下式求算有效扩散系数D*[30]

    D=Lnc0top×dQtdt,
    (3)
    图  2  扩散达到稳定状态单位面积累积溶质通量与时间关系[31]
    Figure  2.  Schematic diagram of steady state diffusion results[31]

    式中,D*为溶质的有效扩散系数(m2/s),L为试样高度(m),n为试样孔隙率,c0-top为试样顶部的初始溶液的溶质浓度,Qt为单位面积累积溶质通量(mg/m2),

    Qt=1ANtj1Δmj=1ANtj=1cbottom,jΔVj,
    (4)

    式中,∆mj为∆t时间间隔内收集到渗滤液中溶质质量的增量(g),∆V为∆t间隔内收集的渗滤液的体积增量(L),cbottom为一定时间间隔内底部收集液中的溶质浓度(mol/L)。

    阻滞因子Rd可由以下公式求算[30]

    Rd=6DL2×TL,
    (5)

    式中,TLQtt曲线横坐标截距,kQt随着时间变化的直线段斜率。

    图 3为MSBS试样两端收集到的渗滤液电导率EC和pH值与时间之间的关系。由图可知,试样的初始pH介于9.5~10,电导率介于220~240 μS/cm,随着时间的增加,渗滤液中pH和EC逐渐减少,pH在10天后逐渐接近稳定值6.85,EC在13 d之后逐渐接近稳定值40 μS/cm。可以发现,当试验时长为6 d时,此时的EC值低于初始EC值的50%,已达到冲刷试验的终止条件,可进行下一阶段试验。

    图  3  去离子水冲刷阶段渗滤液pH/EC随时间的变化关系
    Figure  3.  Variation of pH/electrical conductivity of effluent solution with time at flushing stage using DIW

    图 4为MSBS试样两端收集到Pb(NO3)2污染液pH随着时间的变化。可以看出第一阶段的基准压差阶段(图中的DIW冲刷阶段)试样顶部和底部的pH相差不大,介于6.6~6.9。在之后的化学渗透压差阶段,pH随着试验时间的增加不断减小,两端的pH值差值不断增大,且顶部的下降速率明显低于底部,最终两端的pH均趋于稳定。分析上述现象原因是由于顶部注入液由低浓度更换到高浓度后,底部仍保持原速率的蒸馏水注入,这时试样两端的Pb2+污染液浓度差发生了显著的变化,在化学势差的作用下顶部Pb2+向底部扩散。在整个过程中底部Pb2+浓度的增速明显滞后于顶部收集液中的离子浓度,两端收集液中pH均呈现下降趋势,且顶部下降速率明显高于底部。

    图  4  MSBS屏障试样上下两端收集液pH随时间的变化关系
    Figure  4.  Variation of pH in liquids collected from upper and lower ends of MSBS specimen with time

    图 5为MSBS试样两端收集到Pb(NO3)2污染液EC随着时间的变化。可以看出第一阶段的基准压差阶段(图中的DIW冲刷阶段)试样顶部和底部的EC介于11~40 μS/cm,在Pb2+污染液浓度为0.5~2 mg/L的初始化学渗透压差阶段,试样两端的收集液EC值基本没有明显变化。在化学渗透压差中后阶段,当第5、第11、第18、第25、第32天更换更高浓度的Pb(NO3)2污染液时(对应冲洗Pb2+浓度分别为0.5,1,2,10,50 mg/L),试样两端的收集液EC值迅速增大,且顶部收集液的EC增速明显高于底部,两端EC差值不断加剧。在两端收集液的电导率增大到一定峰值后都逐渐趋于稳定不变。分析原因是由于,跟pH变化原理类似,在顶部注入液由低浓度更换到更高浓度时,顶部收集液的EC值迅速增大,之后在化学势差的作用下,顶部高浓度Pb2+经试样扩散到底部,但扩散速率较慢,因此底部收集液中EC缓慢增加。最后随着冲刷时间的增加,两端的收集液浓度达到平衡,EC值也趋于稳定。

    图  5  MSBS屏障试样上下两端收集液EC随时间的变化关系
    Figure  5.  Variation of electrical conductivity in liquids collected from upper and lower ends of MSBS specimen with time

    图 6为MSBS试样两端收集到Pb(NO3)2污染液浓度随着时间的变化。可以看出第一阶段的基准压差阶段(图中的DIW冲刷阶段),试样本身不含有重金属离子,因此收集液中的Pb2+浓度为0。在化学渗透差压阶段,随着试样顶部不同浓度Pb(NO32污染液的注入,试样顶部Pb2+浓度迅速增大,底部Pb2+浓度变化明显滞后于顶部收集液,且增长缓慢,试样上下两端的Pb2+浓度差也在注入更高浓度污染液的情况下不断加剧。在试验冲刷前期第5天时,当Pb2+浓度由0.1 mg/L更换到0.5 mg/L后,顶部的收集液Pb浓度由0.098 mg/L迅速增大到0.344 mg/L,并在5 d后逐渐稳定至0.5 mg/L。在25 d和32 d注入液浓度更换到10 mg/L和50 mg/L后,顶部收集液迅速增大至7.51 mg/L和48.98 mg/L,增幅远大于前3次更换注入液时的浓度变化。整个冲刷阶段过程中,底部收集液的浓度在前3次更换更高浓度的注入液时基本没有变化,仅仅在最后2次更换注入液时才有了显著的变化,但增幅明显小于顶部收集液浓度的变化。

    图  6  MSBS屏障试样上下两端收集液中Pb2+浓度与时间的变化关系
    Figure  6.  Variation of Pb2+ concentration in liquids collected from upper and lower ends of MSBS specimen with time

    随着试样顶底高浓度中Pb2+的注入,两端的污染液浓度差不断增大,实测压差也发生着变化。图 7为Pb(NO3)2污染液作用下试样两端的收集液中化学渗透压差随时间变化规律。可以看出基准压差阶段(图中的DIW冲刷阶段),试样两端的压差保持10Pa不变。Shackelford等[32]研究认为基准差压测试阶段腔室两端的化学渗透压力差主要由于以下3种原因造成:①试样两端透水石的孔隙结构不同导致两端实际的注液压力存在差异;②试样两端的透水石上检测到了残留的金属阳离子;③试样并非完全均匀, 上下两端的结构差异导致。在试验阶段,当腔室顶部注入一定浓度的Pb(NO3)2污染液后,试样实测压差先是急剧增大到一定数值后快速下降,最后趋于稳定。该现象可能是由于顶部注入的污染液由低浓度更换到高浓度后,试样两端的压力差显著增大,在化学势压差的作用下Pb2+逐渐扩散到底部,底部离子浓度缓慢增加,两端压力差逐渐降低,直至污染物离子扩散趋于稳定,化学渗透压差不再变化。

    图  7  MSBS屏障试样上下两端化学渗透压差随时间变化关系
    Figure  7.  Variation of chemico-osmotic pressure of MSBS specimen with time

    图 8表示试样的化学渗透膜效率系数随着时间的变化。可以看出,膜效率系数随着冲刷时间(注入溶液的浓度)的增加而逐级递减。随着Pb(NO3)2污染液在试样中扩散逐渐稳定,试样的膜效率系数也不再变化。在前3次注入液Pb2+浓度提高到0.5, 1, 2 mg/L更高浓度时,膜效率系数每次变化较小,分别降低了20.02%,18.75%和11.54%。而注入液更换到更高浓度10,50 mg/L时,膜效率系数显著下降相应降低了39.13%,44.57%。分析原因在于:当Pb2+浓度显著增加后,试样中膨润土中的蒙脱石黏土矿物的双电层厚度显著压缩,膨润土颗粒出现团聚现象,同时MSBS水化产物对Pb2+的吸附能力也减弱,上述影响程度随Pb2+浓度升高而愈加显著。对于某一给定Pb2+浓度,随着试验时间增加,膨润土中的蒙脱石黏土矿物的双电层厚度不再变化,MSBS水化产物对Pb2+吸附不再继续受影响。因此,膜效率系数也随试验时间而稳定。

    图  8  试样化学渗透膜效率系数随着时间的变化
    Figure  8.  Change of chemico-osmotic efficiency coefficient of MSBS specimen with time

    图 9统计了国内外学者[14, 25, 33-36]使用不同岩土工程材料在Pb2+,Zn2+,Ca+,Na+,K+作用下化学渗透膜效率随污染物浓度变化的规律。其中梅丹兵使用的是12.5%膨润土的砂-膨润土竖向屏障(SB),张润使用的是2%SHMP改性的10%膨润土的砂–膨润土竖向屏障(SHMP-SB),傅贤雷使用的是0.6%PAC改性的10%膨润土的砂–膨润土竖向屏障(PAC-SB),Tang使用的是95%Fukakusa Clay与5%钠基膨润土组成的改性击实黏土,Shackelford和Malusis使用的是膨润土防水毯(GCL)。以上掺量均为干重比。结果发现,屏障材料的化学渗透膜效率系数基本都随着溶质(阳离子)浓度的增加而减小。在溶质浓度较低时,SB屏障材料保持着较高化学渗透膜效率;然而随着浓度的开始升高,膜效率系数急剧下降,直到污染液的浓度超过5 mmol/L后,膜效率下降速率明显变缓,最后呈现趋于一致的趋势。与其他屏障材料一样,MSBS屏障材料在重金属污染液作用下具有明显的化学渗透膜效应。GCL的化学渗透膜效率系数在研究的浓度范围内高于SB以及改性黏土,这是由于GCL相对于SB和改性黏土具有更高的膨润土含量(SB为5%~10%,改性黏土为5%,GCL为100%)。

    图  9  不同岩土工程材料化学渗透膜效率系数随Pb2+浓度的变化
    Figure  9.  Relationship between chemico-osmotic efficiency coefficient and concentration of Pb2+ for different barrier materials

    图 10表示6种Pb2+离子浓度作用下MSBS竖向阻隔屏障材料单位面积累计Pb2+质量通量Q随化学渗透时间t的关系。不同浓度的Pb2+作用下,屏障材料单位面积累计通量都与冲刷时间成正比,此时可认为Pb2+扩散已稳定。对比6种浓度下Qt拟合直线的斜率k可以发现:随着Pb2+浓度增加,斜率k随着Pb2+浓度增加而增大,斜率由0.037增加到20.01,扩大了540倍;扩散稳定时间TL也由2.35 d减小到2.07 d。

    图  10  MSBS屏障材料单位面积累计Pb2+质量通量Q随化学渗透时间t的关系
    Figure  10.  Variation of unit area accumulated quality flux of MSBS barrier specimen with time

    图 1112分别表示MSBS竖向屏障中Pb2+污染物离子时有效扩散系数D*和屏障材料的阻滞因子Rd随污染物浓度的变化规律,并比对了Malusis等[28]使用的SB竖向屏障和Henning等[29]使用的SHMP-SB竖向屏障。由图 11可知SB竖向阻隔防污屏障中Pb2+的有效扩散系数随离子浓度增加显著增大,而在MSBS中Pb2+的有效扩散系数略有增加,且Pb2+通过土-膨润土竖向屏障(SB)的有D*明显高于MSBS竖向阻隔屏障。分析原因可能是由于这是由于固化90d后MSBS已经充分水化,水化生成的水合硅酸钙C-S-H和类水滑石相(Mg6Al2CO3(OH)16•4H2O,简称Ht)凝胶物质能够有效地填充MSBS的孔隙空间,降低了孔隙之间的连通性,形成的密实的微结构能够阻碍重金属离子的通过,降低屏障材料的扩散系数[11, 37-38]

    图  11  不同竖向阻隔材料有效扩散系数随Pb2+浓度的变化
    Figure  11.  Relationship between effective diffusion coefficient and concentration of Pb2+for different vertical cutoff wall backfills
    图  12  不同竖向阻隔材料阻滞因子随Pb2+浓度的变化
    Figure  12.  Relationship between retardation factor and concentration of Pb2+ for different cutoff wall backfills

    图 12比较了不同Pb2+浓度下MSBS和传统SB阻隔屏障材料的阻滞因子Rd。可以看出MSBS竖向阻隔屏障对于重金属Pb2+Rd明显低于SB竖向屏障。随着Pb2+浓度的增加,MSBS和传统SB的Rd呈现不同的变化规律。SB屏障材料的Rd随着污染物离子浓度的增大显著减少,而MSBS屏障材料的Rd基本不随污染物离子浓度的变化而变化。分析原因可能是由于MSBS中膨润土相对于SB掺量较少,且膨润土与MgO-GGBS的碱性体系下黏土矿物(蒙脱石为主)发生溶蚀,导致MSBS中实际膨润土含量更低。因此对于重金属的吸附、固定和络合弱于SB。但MSBS自身掺有一定量的GGBS,其水化产物生成的凝胶材料能够对膨润土颗粒分子起胶接、包裹作用,这一定程度上能够保护膨润土颗粒避免被金属阳离子“攻击”,使得其对于Pb2+的阻滞效果基本不随污染物浓度的变化影响,在一定浓度的Pb2+溶液中的稳定性要明显好于其他阻隔材料。

    本文采用MgO作为碱激发剂用以激发粒化高炉矿粉制备成氧化镁碱激发矿渣–膨润土新型竖向阻隔屏障(MSBS)。通过土柱化学渗透试验,研究新型隔离屏障阻隔重金属Pb2+污染物的化学渗透膜效应和溶质扩散运移特性,主要得出3点结论。

    (1)随着试样顶部注入浓度更高的Pb(NO3)2溶液,MSBS屏障材料两端的化学渗透压差压差显著增大,达到一定峰值后缓慢下降并逐渐稳定,最终试样两端渗透压差逐渐平衡。

    (2)MSBS具有良好的半透膜效应,化学渗透膜效率系数随着Pb2+浓度的增加而减小。同一种金属阳离子作用下,各类SB竖向屏障材料的化学渗透膜效应系数在较高浓度(约5 mM)下趋于一致。

    (3)随着Pb2+浓度的增加,MSBS针对Pb2+的有效扩散系数略有增加,但变化幅度明显低于SB竖向阻隔屏障材料。此外,MSBS的阻滞因子随着Pb2+浓度的增加保持相对稳定,而SB竖向阻隔屏障对于Pb2+的阻滞因子则随Pb2+浓度增加而显著降低。

  • 图  1   隐式有限差分法示意图

    Figure  1.   Full-implicit finite difference method

    图  2   相同时间不同分数阶阶次α对于反常扩散的影响

    Figure  2.   Effects of fractional order α on anomalous diffusion under same time

    图  3   含水率θ随时间t、位置x变化图

    Figure  3.   Fitting results of changes of moisture content with time and displacement

    图  4   不同时间拟合结果与试验数据对比图

    Figure  4.   Comparison between fitting results and experimental data of changes of moisture content under different time

    图  5   空间分数阶渗流模型与Richards方程拟合结果对比

    Figure  5.   Comparison of fitting results between space-fractional order and Richard's equation

  • [1] 陈正汉. 非饱和土与特殊土力学的基本理论研究[J]. 岩土工程学报, 2014, 36(2): 201-272. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201402002.htm

    CHEN Zheng-han. On basic theories of unsaturated soils and special soils[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(2): 201-272. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201402002.htm

    [2]

    BUCKINGHAM E. Studies on the movement of soilmoisture[J]. US Department of Agriculture, Bureau of Soils B, 1907, 38: 1-61.

    [3]

    RICHARDS L A. Capillary conduction of liquids through porous medium[J]. Physics, 1931, 1: 318-333. doi: 10.1063/1.1745010

    [4]

    RICHARDS L A. The usefulness of capillary potential to soil-moisture and plant investigations[J]. J Agr Research, 1928, 37: 719-742.

    [5]

    NEUMAN S P. Universal scaling of hydraulic conductivities and dispersivities in geologic media[J]. Water Resources Research, 1990, 26(8): 1749-1758. doi: 10.1029/WR026i008p01749

    [6]

    HATANO Y, HATANO N. Dispersive transport of ions in column experiments: an explanation of long‐tailed profiles[J]. Water Resources Research, 1998, 34(5): 1027-1033. doi: 10.1029/98WR00214

    [7]

    HAGGERTY R, MCKENNA S A, MEIGS L C. On the late‐time behavior of tracer test breakthrough curves[J]. Water Resources Research, 2000, 36(12): 3467-3479.

    [8] 孙洪广, 陈文, 蔡行. 空间分数阶导数“反常”扩散方程数值算法的比较[J]. 计算物理, 2009, 26(5): 719-724. doi: 10.3969/j.issn.1001-246X.2009.05.011

    SUN Hong-guang, CHEN Wen, CAI Xing. Comparative study of numerical algorithms for ‘anomalous’ diffusion equation with spatial fractional derivatives[J]. Chinese Journal of Computational Physics, 2009, 26(5): 719-724. (in Chinese) doi: 10.3969/j.issn.1001-246X.2009.05.011

    [9] 常福宣, 陈进, 黄薇. 反常扩散与分数阶对流–扩散方程[J]. 物理学报, 2005, 54(3): 1113-1117. doi: 10.3321/j.issn:1000-3290.2005.03.020

    CHANG Fu-xuan, CHEN Jin, HUANG Wei, et al. Anomalous diffusion and fractional advection-diffusion equation[J]. Acta Physica Sinica, 2005, 54(3): 1113-1117. (in Chinese) doi: 10.3321/j.issn:1000-3290.2005.03.020

    [10] 孙洪广, 常爱莲, 陈文, 等. 反常扩散:分数阶导数建模及其在环境流动中的应用[J]. 中国科学:物理学力学天文学, 2015, 45(10): 104702. https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK201510002.htm

    SUN Hong-guang, CHANG Ai-lian, CHEN Wen, et al. Anomalous diffusion: fractional derivative equation models and applications in environmental flows[J]. Scientia Sinica, 2015, 45(10): 104702. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK201510002.htm

    [11]

    ZHOU H W, YANG S, ZHANG S Q. Conformable derivative approach to anomalous diffusion[J]. Physica A: Statistical Mechanics and its Applications, 2018, 491: 1001-1013. doi: 10.1016/j.physa.2017.09.101

    [12]

    KHALIL R, AL HORANI M, YOUSEF A, et al. A new definition of fractional derivative[J]. Journal of Computational and Applied Mathematics, 2014, 264: 65-70. doi: 10.1016/j.cam.2014.01.002

    [13]

    PACHEPSKY Y, TIMLIN D, RAWLS W. Generalized Richards' equation to simulate water transport in unsaturated soils[J]. Journal of Hydrology, 2003, 272(1/2/3/4): 3-13.

    [14]

    PACHEPSKY Y A, TIMLIN D. Water transport in soils as in fractal media[J]. Journal of Hydrology, 1998, 204: 98-107.

    [15]

    MILLER E E, MILLER R D. Physical theory for capillary flow phenomena[J]. Journal of Applied Physics, 1956, 27(4): 324-332. doi: 10.1063/1.1722370

    [16]

    OZISIK N. Finite Difference Methods in Heat Transfer[M]. Boca Raton: CRC Press, 1994.

    [17]

    EL ABD A, CZACHOR A, MILCZAREK J. Neutron radiography determination of water diffusivity in fired clay brick[J]. Applied Radiation and Isotopes, 2009, 67(4): 556-559.

    [18] 陈正汉, 谢定义, 王永胜. 非饱和土的水气运动规律及其工程性质研究[J]. 岩土工程学报, 1993, 15(3): 9-20. doi: 10.3321/j.issn:1000-4548.1993.03.002

    CHEN Zheng-han, XIE Ding-yi, WANG Yong-sheng, et al. Experimental study of laws of fulid motion, suction and pore pressures in unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 1993, 15(3): 9-20. (in Chinese) doi: 10.3321/j.issn:1000-4548.1993.03.002

    [19]

    PARLANCE M B, PRASAD S N, PARLANGE J Y, et al. Extension of the Heaslet‐Alksne technique to arbitrary soil water diffusivities[J]. Water Resources Research, 1992, 28(10): 2793-2797. doi: 10.1029/92WR01683

    [20] 梅岭, 姜朋明, 李鹏, 等. 非饱和土的土水特征曲线试验研究[J]. 岩土工程学报, 2013, 35(增刊1): 124-128. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2013S1019.htm

    MEI Ling, JIANG Peng-ming, LI Peng, et al. Soil-water characteristic curve tests on unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(S1): 124-128. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2013S1019.htm

  • 期刊类型引用(5)

    1. 邓圣义,周永伟,蒲诃夫,李育超,闵明. 阻隔墙中考虑半透膜效应的污染物迁移数值解. 浙江大学学报(工学版). 2025(02): 342-350 . 百度学术
    2. 薛强,杜延军,胡黎明,詹良通,李江山. 环境土力学与工程研究进展. 土木工程学报. 2025(03): 83-112 . 百度学术
    3. 陈杰,唐占元,安之焕,高健,马渊博. 矿粉固化黄土的力学性能及机理研究. 公路. 2025(03): 33-41 . 百度学术
    4. 刘松玉,蔡国军,张炜,周宏磊,邓永锋. 岩土工程勘察、测试与评价进展. 土木工程学报. 2024(10): 108-124 . 百度学术
    5. 张志红,杨灏闻,郑九州. 高岭土-膨润土化学渗透膜效应试验及微观机理分析. 岩土工程学报. 2023(09): 1963-1970 . 本站查看

    其他类型引用(3)

图(5)
计量
  • 文章访问数:  324
  • HTML全文浏览量:  27
  • PDF下载量:  186
  • 被引次数: 8
出版历程
  • 收稿日期:  2018-01-28
  • 网络出版日期:  2022-12-07
  • 刊出日期:  2020-08-31

目录

/

返回文章
返回