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

间歇性变围压循环荷载作用下饱和软黏土变形特性

黄珏皓, 王宏超, 陈健, 付晓东, 阎晓玲, 马超

黄珏皓, 王宏超, 陈健, 付晓东, 阎晓玲, 马超. 间歇性变围压循环荷载作用下饱和软黏土变形特性[J]. 岩土工程学报, 2023, 45(S1): 67-70. DOI: 10.11779/CJGE2023S10038
引用本文: 黄珏皓, 王宏超, 陈健, 付晓东, 阎晓玲, 马超. 间歇性变围压循环荷载作用下饱和软黏土变形特性[J]. 岩土工程学报, 2023, 45(S1): 67-70. DOI: 10.11779/CJGE2023S10038
HUANG Juehao, WANG Hongchao, CHEN Jian, FU Xiaodong, YAN Xiaoling, MA Chao. Effects of intermittent cyclic loading with cyclic confining pressure on deformation behaviors of saturated clay[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(S1): 67-70. DOI: 10.11779/CJGE2023S10038
Citation: HUANG Juehao, WANG Hongchao, CHEN Jian, FU Xiaodong, YAN Xiaoling, MA Chao. Effects of intermittent cyclic loading with cyclic confining pressure on deformation behaviors of saturated clay[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(S1): 67-70. DOI: 10.11779/CJGE2023S10038

间歇性变围压循环荷载作用下饱和软黏土变形特性  English Version

基金项目: 

国家自然科学基金面上基金项目 52079135

中国科学院青年创新促进会项目 2021325

详细信息
    作者简介:

    黄珏皓(1991—),男,副研究员,主要从事软土动力特性及本构模型方面的研究工作。E-mail:jhhuang@whrsm.ac.cn

  • 中图分类号: TU442

Effects of intermittent cyclic loading with cyclic confining pressure on deformation behaviors of saturated clay

  • 摘要: 交通荷载由循环加载和间歇停振两部分组成,间歇停振对土体动力特性的影响不容忽视。交通荷载引起的应力场包含循环变化的轴向偏应力和围压。因此,通过对软黏土开展一系列间歇性变围压循环加载试验,对饱和软黏土在间歇性加载和循环围压作用下的变形特性进行了研究,分析了循环围压和间歇停振阶段排水条件对软黏土累积轴向应变的影响。研究结果表明,不同试验条件下软黏土累积轴向应变增量随加载阶段的变化趋势保持一致,即累积轴向应变增量随加载阶段的增加逐渐减小。累积轴向应变增量随循环围压的增大而减小。间歇停振阶段不同排水条件下,试样累积轴向变形变化规律也有所区别:部分排水条件下,循环加载阶段产生的超孔隙水压力在间歇停振阶段消散,试样累积轴向应变在间歇停振阶段有一定程度的增长;不排水条件下,循环加载阶段产生的一部分变形在间歇停振阶段得以慢慢恢复,试样累积轴向应变有所减小。上述研究成果可加深对间歇性循环荷载作用下软黏土变形规律的认识。
    Abstract: Traffic loading is composed of cyclic loading and intermittent periods, in which the influences of intermittent period on the dynamic characteristics of soil cannot be ignored. The stress field in subgrade soil caused by the traffic loading includes cyclic deviator stress and cyclic confining pressure. Therefore, a series of cyclic triaxial tests in intermittent cyclic loading mode with cyclic confining pressure are carried out to study the deformation behaviors of saturated clay, and the effects of the cyclic confining pressure and drained conditions during the intermittent period are investigated. The results show that the relationship between the accumulated axial strain increment and the loading stages are similar under different test conditions: the accumulated axial strain increment gradually decreases as the loading stage increases. Additionally, the variation of the accumulated axial strain increment is also affected by the cyclic confining pressure, in which the accumulated axial strain increment decreases with the increasing cyclic confining pressure. Besides that, the variations of accumulated axial strain are also different under different drained conditions during the intermittent period: the excess pore water pressure induced in cyclic loading period dissipates under partially drained conditions during the intermittent period, leading to that the accumulated axial strain of soft clay increases in a certain extent during the intermittent period; part of the deformation produced during the cyclic loading period can recover under undrained conditions during the intermittent period, and the accumulated axial strain of soft clay decreases. The above test results may deepen the understanding of the deformation behaviors of soft clay under intermittent cyclic loading.
  • 对于垃圾填埋场、危险废物处置场等污染场地的防渗阻隔问题,底部衬垫系统是污染场地防渗系统的重要组成部分,主要用于阻止和延缓渗滤液中污染物的向下扩散与迁移[1-4]。压实黏土是底部衬垫系统的重要组成材料,由于其具有较低的渗透性和较好的化学相容性,因而被广泛用于污染场地的防渗系统中。尽管如此,危害性较大的污染物仍可通过对流、扩散等方式击穿压实黏土衬垫,从而对周围土体及地下水造成污染[1-9]。为此,许多学者对污染物在压实黏土衬垫中的运移问题展开了研究[5-9]

    在污染物一维运移的解析理论方面,一些学者探究了压实黏土衬垫中污染物的运移问题[6-7]。如陈云敏等[6]求解得到了污染物在成层黏土衬垫中一维扩散运移的解析解;张文杰等[7]建立了污染物在压实黏土衬垫中的一维对流–扩散–吸附解析模型。与此同时,针对不同情况,一些学者提出了污染物在包含压实黏土衬垫的复合衬垫系统中一维运移的解析解[5, 8-9]。如Foose[5]在考虑扩散–吸附情况下,获得了有机污染物在土工膜与黏土层组成的复合衬垫中一维瞬态运移的解析解;Pu等[9]针对底部衬垫系统为复合衬垫系统的情况,采用解析方法研究了有机污染物在压实黏土衬垫中的运移问题。然而,上述研究均忽略了温度对污染物运移过程的影响。

    现场测定数据显示,污染场地(如垃圾填埋场)内污染土体中发生的化学反应会使得场地内部的温度高达90℃,而场地外部土体和地下水的温度一般较低,因此压实黏土衬垫中存在温度梯度[10-12]。温度梯度的存在不仅会使得污染物在压实黏土衬垫中发生热扩散,还会影响土体的物理力学性质以及与污染物运移过程相关的参数,进而影响污染物运移过程[12-18]。对于非等温分布条件下压实黏土衬垫中污染物的运移问题,目前相关的理论和试验研究较少[15-17]

    在非等温分布条件下污染物运移问题的试验研究中,Rosanne等[19-20]通过开展试验研究指出,压实黏土中存在因温度梯度引起的污染物热扩散现象。在理论方面,已有学者开展了考虑不同情况的相关解析理论研究[12-17]。如吴珣等[12]通过某些简化假定推导得到了非等温分布条件下有机污染物一维运移的解析解;张春华[14]通过假定污染物的扩散系数等参数不随温度发生变化,发展了考虑热扩散、降解等作用下有机污染物在压实黏土衬垫中一维运移的解析解;Yan等[15]建立了考虑分子扩散–吸附–热扩散情况下污染物在多孔介质中一维运移的解析模型;Peng等[16-17]研究了考虑热扩散情况下复合衬垫系统中污染物的运移问题,并提出了双层复合衬垫和三层复合衬垫中污染物一维运移的解析模型。但上述研究均未考虑土体固结对污染物运移过程的影响。

    在底部衬垫系统服役过程中,由于其上部通常会堆填大量的固废垃圾体,这使得压实黏土衬垫会受到较大竖向应力,从而使得压实黏土衬垫发生固结变形[21-23],进而影响污染物运移过程。因此,需考虑土体固结对污染物运移过程的影响。目前,有关压实黏土衬垫中固结与污染物运移的耦合问题已有较多的研究[24-26]。然而,这些研究大多忽略了温度的影响。田改垒等[27]基于小变形固结理论建立了可考虑热效应的渗透–扩散–固结耦合模型,但该模型忽略了温度变化对污染物运移过程中相关参数的影响,且所对应的工况与污染场地中的工况不一致。此外,为获得污染物在压实黏土衬垫中一维运移时的解析解,上述研究中多采用线性等温吸附方程来描述土颗粒的吸附特征,这与黏土颗粒实际的非线性特性存在区别[28-30]

    针对以上研究情况,本文在已有研究的基础上,建立了非等温分布条件下压实黏土衬垫中固结与污染物运移的耦合模型,并采用有限差分法对该耦合模型进行了求解。通过将耦合模型的计算结果分别与Rosanne等[19]所开展的热扩散试验结果、张文杰等[7]所求解析解的计算结果、Alshawabkeh等[31]所提数值模型的计算结果进行对比分析,验证了所建耦合模型的正确性。基于所建耦合模型,分析了相关参数变化对污染物运移过程的影响。

    图 1为非等温分布条件下压实黏土衬垫(CCL)中固结与污染物运移耦合模型的计算简图。图 1中,从上至下依次为堆填体、渗滤液收集和移除系统、渗滤液、CCL、水平导流层及不透水基底层。这里假定CCL始终处于饱和状态,其厚度为L,坐标系z由CCL上边界垂直向下,hw为CCL上部的渗滤液水头;C0为污染物浓度,这里假定污染物浓度保持恒定[16-17, 30]。为研究非等温分布条件对固结及污染物运移过程的影响,同时考虑到温度的传导速度较快,这里假定CCL中的温度T随深度呈线性分布,TaTb分别为上部和底部的温度[12-18]。此外,考虑到污染场地在运行过程会不断堆填固废垃圾体,从而使得CCL受到的竖向应力不断增大,这里设定堆填体产生的加荷速率为Q,并记堆填总时间(即加荷时间)为tc

    图  1  耦合模型的计算简图
    Figure  1.  Calculation diagram for coupling model

    由于底部衬垫系统的水平面积较大,CCL的固结问题可近似为一维问题。已有试验研究表明,CCL的渗透系数会随温度的升高而增大[32-33],但温度变化对土体压缩曲线和压缩模量的影响较小[34-35]。基于试验研究结果,作如下基本假定:①CCL是均质、各向同性的,且处于饱和状态;②CCL中固相(土颗粒)和液相(孔隙水)均不可压缩,土体变形固结是由于液相的排出所引起;③液相的排出会导致土体的孔隙率减小,但土体的变形量远小于土体的厚度,即小变形假定成立;④液相中污染物的浓度很低,忽略污染物浓度变化对土体物理力学性质的影响;⑤土中液相的渗流符合Darcy定律;⑥温度变化不会导致CCL的压缩性改变,仅影响CCL的渗透系数。

    为研究污染物在CCL中的一维运移问题,在上述基本假定的基础上,还需补充以下3点假定:①污染物在CCL中的扩散遵循Fick第二定律,土颗粒对污染物的吸附状态已达到平衡吸附状态;②污染物为单一离子或有机物;③污染物的运移方式考虑分子扩散、对流、机械弥散、吸附和热扩散。

    在CCL一维固结过程中,液相和固相的质量守恒方程有[36]

    (nvtcρf)z=(nρf)t,
    (1a)
    [(1n)vsρs]z=[(1n)ρs]t,
    (1b)

    式中,n为CCL的孔隙率,vtc为液相中因土体固结产生的流动速度,ρf为液相的密度,vsρs分别为固相的移动速度和密度,t为时间。

    由于固相和液相的密度保持不变,根据式(1a),(1b)有:

    (nvtc)z=nt,
    (2a)
    [(1n)vs]z=(1n)t
    (2b)

    根据式(2)可得

    vsz+vdcz=0,
    (3)

    式中,vdc为由于土体中因固结引起的达西流速,vdc=n(vtcvs)

    根据达西定律,土体中液相的流动速度为

    vd=vdc+vdh,
    (4a)
    vdc=kvγwuexz,
    (4b)
    vdh=kvγwuhz,
    (4c)

    式中,kv为液相的渗透系数,γw为液相的重度,uex为因外荷载作用产生的超孔隙水压力,其会导致土体发生固结作用,uh为渗滤液水头产生的超孔隙水压力,其会导致液相发生流动,但不会使得土体发生固结变形[36]vdh为因渗滤液水头产生的达西流速。

    根据CCL中所作的非等温分布假定可得

    T(z)=Ta+Mz,
    (5)

    式中,M为CCL中的温度梯度(℃/m),M= (TbTa)/(TbTa)LL

    对于液相的渗透系数kv,参考文献[37],渗透系数kv与温度T的关系可近似表达为

    kv=(0.029T+0.420)kv(T0=20C),
    (6)

    式中,kv(T0)为温度为T0时的渗透系数。

    结合式(5),(6),渗透系数kv与深度z关系为

    kv(z)=kv(T0)(Nz+Qa),
    (7)

    式中,N=0.029MN为渗透系数梯度(1/m),Qa=0.029Ta+0.420kv(0)=kv(T0)Qakv(0)为CCL上部的渗透系数。

    对于因渗滤液水头产生的达西流速vdh,其不产生固结效应,这里假定渗滤液水头hw保持不变,并在该水头下达到了稳定渗流状态,则有

    vdh=kv(z)γwduhdz
    (8)

    对式(8)进行一定的变换,并利用式(7)的表达式进行积分可得

    vdh=Nkv(T0)hwln[(N/Qa)L+1]
    (9)

    在CCL中,根据液相中的渗流连续性条件可得

    vdcz=mvσt,
    (10)

    式中,mv为体积压缩系数,σ为有效应力。

    根据式(4b),(10)可进一步写为

    z(kv(z)γwuexz)=mvσt
    (11)

    根据土体的有效原理可知:

    σ=σuexuhus,
    (12)

    式中,σ为总应力,us为静水压力。

    由于渗滤液水头下的渗流已达稳定渗流状态,且静水压力不随时间变化,则根据式(12)和上部荷载的加荷条件,式(11)可进一步改写为

    z(kv(z)γwuexz)={mvuextmvQ(t<tc)mvuext(ttc) 。 
    (13)

    式(13)即为非等温分布条件下CCL的一维固结控制方程。

    当CCL发生固结变形时,液相会发生流动,固相(土颗粒)也会发生相对移动。根据质量守恒定律,考虑分子扩散、对流、机械弥散、吸附和热扩散时污染物的一维运移方程可写为[12-17, 36]

    (nC)t+[(1n)ρsS]t=z(nDCz)(nvtC)z[1n)vsρsS]z+z(nDTTz),
    (14)

    式中,C(z,t)为液相中污染物的浓度,S为单位质量土颗粒所吸附的污染物质量,D为水动力弥散系数,D=De+DmdDe为有效扩散系数,Dmd=αLvtDmd为机械弥散系数,αL为纵向弥散度,vt=vtc+vthvt为液相的流动速度,vth为液相中因渗滤液水头产生的流动速度,vth=vdh/vdhnnDT为污染物的热扩散系数。

    根据污染物的一维运移机理,以垂直向下方向为正方向的污染物通量J(z,t)的表达式可写为

    J=nDCz+[nvtC+(1n)vsρsS]nDTTz
    (15)

    对于SC的关系,参考已有研究结果,这里采用Freundlich非线性平衡吸附模型来描述固相对污染物的非线性吸附特性[28-30],即

    S=KfCF,
    (16)

    式中,Kf为与吸附能力相关的Freundlich吸附系数,F为拟合常数。

    当拟合常数F=1时,则式(16)可退化为经典的线性等温吸附模型,即

    S=KfC=KdC,
    (17)

    式中,Kd为线性吸附系数。

    对于有效扩散系数De。一方面,有研究指出有效扩散系数De会随孔隙率n的减小而减小[28, 30];另一方面,有效扩散系数De会随温度T的增大而增大[15, 38-39]。综合考虑两者的影响,De可表示为

    De=D0(T0)nβTT0μT0μT,
    (18)

    式中,D0(T0)为温度为T0时的自由扩散系数,β为经验指数,μT0μT分别为温度T0T时污染物的动力黏度。

    基于文献[15, 3839]的研究,De可近似表示为

    De=D0(T0)nβ[1+A(TT0)],
    (19)

    式中,A为与有效扩散系数相关的温度系数。

    将式(5)代入式(19)可得

    De=D0(T0)nβ(Gz+H),
    (20)

    式中,G=AMH=1+A(TaT0)

    对于污染物的热扩散系数DT,参考文献[13~17]和式(20),其表达式为

    DT=CSTDe=CSTD0(T0)nβ(Gz+H),
    (21)

    式中,ST为索雷特系数。

    根据上述关系式,式(14)可进一步改写为

    [nC+(1n)ρsKfCF]t=[D0(T0)nβ+1(Gz+H)+nαLvt].2Cz2+{[D0(T0)nβ+1(Gz+H)+nαLvt]z[nvt+(1n)vsρsKfCF1]+MSTD0(T0)nβ+1(Gz+H)}CzC{[nvt+(1n)vsρsKfCF1]MSTD0(T0)nβ+1(Gz+H)}z
    (22)

    同理,根据关系式,污染物通量J(z,t)的表达式可进一步改写为

    J=[D0(T0)nβ+1(Gz+H)+nαLvt]Cz+[nvtC+(1n)vsρsKfCFCMSTD0(T0)nβ+1(Gz+H)]
    (23)

    为便于分析,这里以渗滤液水头下的渗流达到稳定时刻为初始时刻,并假定在该时刻开始堆填固废垃圾体。因此,CCL一维固结的初始条件可写为

    uex(z,0)=0
    (24)

    对于CCL一维固结的边界条件,这里可将上下边界条件均视为完全排水边界情况,则有:

    uex(0,t)=0
    (25)
    uex(L,t)=0
    (26)

    对于污染物在CCL中的一维运移问题,CCL在初始时刻未受到污染,因而可假定其初始条件为

    C(z,0)=0
    (27)

    由于CCL的上边界直接与渗滤液接触,而下边界为水平导流层,因而上下边界条件可写为[6, 16-17]

    C(0,t)=C0
    (28)
    C(L,t)=0
    (29)

    上述即为非等温分布条件下压实黏土衬垫(CCL)中考虑固结与污染物运移的耦合模型,该耦合模型包括CCL的一维固结控制方程式(13)和CCL中污染物的一维运移方程式(22),以及相应的初始条件和边界条件式(24)~(29)。

    对于非等温分布条件下CCL中固结和污染物运移的耦合模型,由于耦合模型中控制方程的系数均是变化的,属于变系数偏微分方程,因而以下将采用有限差分法对耦合模型进行求解。

    根据式(7),对控制方程式(13)进行整理可得

    Cv0(Nz+Qa)2uexz2+Cv0Nuexz={uextQ(t<tc)uext(ttc),
    (30)

    式中,Cv0=kv(T0)mvγwCv0为温度为T0时CCL的固结系数。

    ΔzΔt分别为计算空间步长和时间步长,并将空间坐标和时间坐标分别进行I等分和J2等分,则有zi=iΔzΔz=L/Iti=iΔtj=0,1,2,,J1, J2Δt=tc/J1=t0/J2t0为给定时间。因此,控制方程式(30)的Crank- Nicholson型隐式差分格式可写为

    Cv0(NΔzi+Qa)(ujexi+12ujexi+ujexi12Δz2+uj+1exi+12uj+1exi+uj+1exi12Δz2)+Cv0Nujexi+1ujexi12Δz={uj+1exiujexiΔtQ(jJ1)uj+1exiujexiΔt(j>J1)
    (31)

    相应的初始条件和边界条件可写为

    u0exi=0
    (32)
    ujex0=0
    (33)
    ujexI=0
    (34)

    利用式(29)中的隐式差分格式以及相应的求解条件式(32)~(34),即可对非等温分布条件下CCL中的一维固结过程进行求解。

    对于CCL中固相的移动速度vs,对式(3)进行整理,两边积分可得

    vs(L,t)vs(z,t)=vdc(z,t)vdc(L,t)
    (35)

    由于CCL底部固相不发生移动(即vs(L,t)=0),因此任意位置固相的移动速度vs可表达为

    vs(z,t)=vdc(L,t)vdc(z,t)
    (36)

    基于有限差分解的计算结果和上述关系式,即可确定达西速度vdc、移动速度vs和孔隙率n等参数。与此同时,利用达西速度vdcvdh可以分别确定对应液相的流动速度vtcvth

    参考式(13)中的有限差分方法,对式(22)展开有限差分,其Crank-Nicholson型隐式差分格式可写为

    TjiCj+1iCjiΔt=Xji(Cj+1i+12Cj+1i+Cj+1i12Δz2)+Xji(Cji+12Cji+Cji12Δz2)+YjiCji+1Cji12ΔzZjiCji,
    (37)

    其中,

    Tji=nji+(1nji)ρsKfCjF1i,
    (38)
    Xji=D0(T0)njβ+1i(GΔzi+H)+njiαLvjti,
    (39)
    Yji=MSTD0(T0)njβ+1i(GΔzi+H)[njivjti+(1nji)vjsiρsKfCjF1i]+12Δz{D0(T0)nji+1β+1[GΔz(i+1)+H]+nji+1αLvjti+1D0(T0)nji1β+1[GΔz(i1)+H]nji1αLvjti1},
    (40)
    Zji=12Δz{[nji+1vjti+1+(1nji+1)vjsi+1ρsKfCjF1i+1]MSTD0(T0)nji+1β+1[GΔz(i+1)+H][nji1vjti1+(1nji1)vjsi1ρsKfCji1F1]+MSTD0(T0)njβ+1i1.[GΔz(i1)+H]} 。 
    (41)

    相应的初始条件和边界条件可写为

    C0i=0
    (42)
    Cj0=C0
    (43)
    CjI=0 
    (44)

    对于非等温分布条件下CCL中污染物的一维运移问题,利用差分方程式(37)及相应的求解条件式(42)~(44)即可进行求解。

    对于污染物通量J,其差分格式可写为

    Jji=XjiCji+1Cji12Δz+ΔjiCji,
    (45)

    式中,

    Δji=[njivjti+(1nji)vjsiρsKfCjF1iMSTD0(T0)njβ+1i(GΔzi+H)]
    (46)

    需要说明的是,对于边界处的差分格式,可利用相邻节点展开有限差分。

    基于差分方程式(31),(37),可利用Matlab软件编制相应的程序,从而对CCL中考虑固结与污染物运移的耦合问题展开运算。

    为验证CCL中固结与污染物运移耦合模型的正确性,以下将耦合模型的有限差分解分别与Rosanne等[19]热扩散试验结果、张文杰等[7]所求解析解、Alshawabkeh等[31]所提数值模型进行对比分析。

    Rosanne等[19]对非等温分布条件下压实黏土中溶质的一维运移问题开展了室内试验研究。在该试验中,压实黏土厚度为2 mm,孔隙率为0.6,溶质为NaCl溶液,其初始浓度为1 moL/m3。压实黏土上下边界的温度分别为12,38℃,且上下边界的浓度分别为1.21,8.83 moL/m3。参考其他文献的模拟方法[20],这里采用在24.9℃时测定得到的有效扩散系数De=6.0×1011 m2/s及A=0.025/C来近似描述溶质的扩散特性。在该试验中,测定得到的索雷特系数ST=0.13/C。需说明的是,对于该热扩散试验,尽管所建耦合模型的初始条件和边界条件等与该试验的相关条件存在一些不同之处,但这些不同之处均可通过调整相关参数以实现对热扩散试验过程的模拟。

    图 2所示为所建耦合模型有限差分解的计算结果和Rosanne等[19]热扩散试验结果的对比情况。从图 2可以看出,不同时间下所建耦合模型计算所得的溶质浓度分布情况与Rosanne等[19]的试验结果基本吻合,这一定程度说明了所建耦合模型的正确性。同时,对比结果也表明本文所建耦合模型可较为准确预测溶质在压实黏土中的热扩散过程。

    图  2  耦合模型计算结果和Rosanne等[19]热扩散试验结果对比
    Figure  2.  Comparison between calculated results by proposed coupling model and test results by Rosanne et al.

    张文杰等[7]研究了污染物在CCL中的一维运移问题,并通过假定土颗粒的吸附模型为线性吸附模型推导得到了相应的解析解,但该解答中未考虑土体固结和温度的影响,因而相关系数均假定为常数。为进一步验证所建耦合模型的正确性,将所建耦合模型的有限差分解与张文杰等[7]所求的解析解进行对比。除指定参数外,取如下基本物理力学参数:厚度L=1.0 m,体积压缩系数mv=0.05 1/MPa,固相密度ρs=2.76 kg/m3,液相密度ρf=1.0 kg/m3,初始孔隙率n0=0.42,加荷速率Q=0.3125 kPa/d,加荷时间tc=3200 d,参考温度T0= 20℃,渗透系数kv(T0)= 2.96×10-10 m/s,渗滤液水头hw=1.0 m,底部温度Tb=20℃,温度梯度M=-30℃/m,自由扩散系数D0(T0)= 8.60×10-10 m2/s,经验指数β=1.82,Freundlich吸附系数Kf=0.63 cm3/g,拟合常数F=0.8,索雷特系数ST=0.05/℃,温度系数A=0.025/℃,污染物浓度C0=10 mg/L,纵向弥散度αL= 0.02 m。需说明的是,计算参数的选定参考了相关文献中所采用的数据及相关试验的测定结果[7, 13-17, 28, 30]。污染物在CCL中运移时的初始条件和边界条件取本文所建耦合模型中采用的相应条件,这里不再赘述。

    图 3所示为本文所建耦合模型有限差分解的计算结果和张文杰等[7]所求解析解的计算结果的对比情况。从图 3可以看出,所建耦合模型计算所得的浓度分布情况与张文杰等[7]所求解析解的计算结果十分一致,这说明所建耦合模型可用于预测考虑对流–扩散–吸附时CCL中污染物的一维运移过程。

    图  3  本文耦合模型与张文杰等[7]解析解计算结果的对比
    Figure  3.  Comparison between proposed coupling model and analytical solution proposed by Zhang et al.

    Alshawabkeh等[31]建立了考虑小变形固结与污染物运移的耦合模型(简称Alshawabkeh耦合模型)。相比于本文所建耦合模型,Alshawabkeh耦合模型忽略了固结变形引起的相关参数变化以及机械弥散等因素的影响,且该研究并未考虑温度引起的热扩散对污染物运移过程的影响。为验证本文所提方法用于求解耦合问题的正确性,可将本文所建耦合模型的计算结果与Alshawabkeh耦合模型的计算结果进行对比。除指定参数外,这里设定相关计算参数同3.2节中的参数。需说明的是,由于Alshawabkeh耦合模型中所考虑的吸附模型为等温线性吸附模型[31],因而这里取Kd=KfF=1,其他条件同本文所建耦合模型时所采用的相应条件,这里不再赘述。

    图 4所示为本文所建耦合模型的计算结果和Alshawabkeh耦合模型的计算结果的对比情况。从图 4可以看出,所建耦合模型计算所得污染物通量随时间的变化情况与Alshawabkeh耦合模型的计算结果十分一致,这说明所建耦合模型可用于计算土体固结与污染物运移的耦合问题。

    图  4  本文耦合模型与Alshawabkeh耦合模型计算结果的对比
    Figure  4.  Comparison between proposed coupling model and coupling model proposed by Alshawabkeh et al.

    为认识非等温分布条件下CCL中考虑固结时污染物的运移特性,除指定参数外,这里以3.2节中的参数为例,探究温度梯度、加荷速率及Freundlich吸附系数等因素对污染物在CCL中运移规律的影响。

    图 5描述了不同温度梯度M下污染物浓度随深度的分布规律和污染物底部通量随时间的变化规律。从图 5(a)可知,M的绝对值越大,相同时间内同一深度处污染物的浓度越高,这说明CCL中的温度梯度会加快运移速率。从图 5(b)可以看出,污染物底部通量随M绝对值的增大出现明显增大,如t=60000 d时,相比于M=0 ℃/m时,M为–10,–20,–30,–40 ℃/m时的底部通量分别增大了26.3%,54.2%,84.3%和117%。此外,底部通量变化规律表明,M绝对值的增大会使得污染物运移过程更早趋于稳定。

    图  5  温度梯度M对污染物运移规律的影响
    Figure  5.  Influence of temperature gradient on law of contaminant transport

    出现上述规律的原因在于:温度梯度M绝对值的增大不仅会加快热扩散过程,还会增大污染物在CCL中运移时的有效扩散系数,进而加快了污染物的运移速率,并增大底部通量。分析表明,温度梯度M会加快污染物在CCL中的运移速率,且在一定温度梯度下,污染物的底部通量可达不考虑温度梯度时底部通量的2倍以上。因此,实际工程中应考虑CCL中的温度梯度对污染物运移过程的影响[12-17]

    图 6给出了不同加荷速率Q下污染物浓度随深度的分布规律和污染物底部通量随时间的变化规律。从图 6(a)可看出,在运移初期(如t≤2.0×104 d),相同时间内同一深度处污染物浓度会随Q的增大而减小,这反映了CCL的固结效应一定程度会减慢运移速率;但运移过程达到稳态时(如t=6.0×104 d时),加荷速率Q的增大反而使得污染物的浓度增大。从图 6(b)可知,加荷速率Q越大,相同时间内污染物底部通量越小,如t=6.0×104 d时,相比于不考虑固结效应(Q=0.0 kPa/d)时的情况,Q=0.3125 kPa/d时(对应最终荷载为1 MPa)染物的底部通量减小了约10.0%。与此同时,从图 6(b)中可以看出,不考虑固结效应会低估污染物的运移过程达到稳态的时间。

    图  6  加荷速率Q对污染物运移规律的影响
    Figure  6.  Influences of loading rate on law of contaminant transport

    上述规律表明,CCL的固结效应会减慢污染物的运移速率,并减小运移时的底部通量,这与Pu等分析得到的规律一致[30],但同时,固结效应也会使得运移过程达到稳态时的浓度增大。对于运移速率减慢和底部通量减小的原因,这主要是由于固结效应会使得土体孔隙率减小,进而使得与污染物运移相关的系数减小,如有效扩散系数和热扩散系数减小。对于运移过程达到稳态时浓度增大的原因,这主要是由于当考虑固结效应时,相关参数的比值会发生变化,如对流速度与水动力弥散系数的比值会增大,因而改变了运移过程达到稳态时污染物浓度的分布规律[7, 14]

    图 7给出了不同Freundlich吸附系数Kf下污染物浓度随深度的分布规律和污染物底部通量随时间的变化规律。从图 7中可知,Freundlich吸附系数Kf对污染物运移过程有着显著影响,在运移初期(如t=1.5×104 d),Kf的增大使得相同时间内同一深度处污染物的浓度显著降低,同时也使得相同时间内污染物的底部通量减小。但随着运移过程的进行,当运移过程达到稳态时(如t=6.0×104 d),不同Kf下污染物的浓度分布曲线几乎完全重合,且不同Kf下底部通量趋于同一稳定值。若定义污染物底部通量达0.445 mg/m2/d时所需的时间为运移稳态时间(图 7(b)),则相比不考虑吸附作用时,Freundlich吸附系数Kf=0.63 cm3/g时的运移稳态时间增大了约200%。

    图  7  Freundlich吸附系数Kf对污染物运移规律的影响
    Figure  7.  Influences of Freundlich adsorption coefficient on law of contaminant transport

    出现上述运移规律的主要原因:CCL的吸附作用主要对污染物运移过程起着阻滞作用,因而主要影响运移速率快慢[28-30]。当运移过程达到稳态时,吸附作用对浓度分布规律和底部通量的影响可忽略。研究表明,CCL对污染物的吸附作用有利于延缓运移过程,且与不考虑吸附作用时相比,考虑吸附作用可使得运移过程达到稳态时所需的时间延长3倍及以上。因此,在选择CCL时,应考虑其吸附特性[28, 30]

    (1)CCL中的温度梯度会加快污染物运移速率,相同时间内污染物的浓度和底部通量会随温度梯度M绝对值的增大而增大,且在一定温度梯度下,底部通量可达不考虑温度梯度时底部通量的2倍以上。因此,实际工程中应考虑温度梯度对运移过程的影响。

    (2)加荷速率Q(CL的固结效应)对污染物运移过程的影响是多方面的。一方面,Q的增大会使得土体的孔隙率减小,从而减慢了运移速率,也减小了底部通量;另一方面,Q的增大也会使得污染物运移过程达到稳态时的浓度增大。因此,实际工程中应综合考虑固结效应对运移过程的影响。

    (3)Freundlich吸附系数Kf的增大会减慢污染物在CCL中的运移过程,且与不考虑吸附作用时相比,考虑吸附作用可使得污染物运移过程达到稳态时所需的时间延长3倍及以上。当污染物的运移过程达到稳态时,土颗粒对污染物的吸附作用可忽略。在选择CCL时,应考虑土体的吸附特性。

  • 图  1   多阶段间歇性变围压循环加载试验示意图

    Figure  1.   Schematic illustration of intermittent cyclic loading

    图  2   不同应力路径斜率累积应变变化曲线

    Figure  2.   Variations of accumulated axial strain under intermittent cyclic loading with different values of η

    图  3   不同应力路径斜率累积应变增量随加载阶段变化曲线

    Figure  3.   Variations of accumulated axial strain increment versus number of loading stages with various values of η

    图  4   循环加载和间歇停振阶段累积轴向应变变化曲线

    Figure  4.   Variations of accumulated axial strain under cyclic loading and intermittent periods

    表  1   试验黏土基本物理性质

    Table  1   Physical properties of test clay

    重度γ/(kN·m-3) 含水率w/% 液限wL/% 塑限wp/% 塑性指数Ip 渗透系数K/(10-7cm·s-1) 孔隙比e
    17.60 48.60 51.90 19.80 32.10 2.26 1.30
    下载: 导出CSV

    表  2   循环三轴试验方案

    Table  2   Programs of cyclic triaxial tests

    停振阶段排水条件 编号 η 加载次数 停振时长/s
    不排水 U01~U04 0.33,1.00,1.50,2.00 1000×4 3600
    部分排水 P01~P04
    下载: 导出CSV

    表  3   停振阶段不同排水条件下各加载阶段对应累积应变

    Table  3   Accumulated axial strains of each loading stage under different intermittent drained conditions

    停振阶段排水条件 加载阶段 循环加载阶段完成后对应累积轴向应变/% 停振阶段完成后对应累积轴向应变/%
    不排水 1 0.668 0.624
    2 0.188 0.147
    3 0.139 0.100
    4 0.119 0.083
    部分排水 1 0.753 0.789
    2 0.290 0.300
    3 0.194 0.187
    4 0.148 0.126
    下载: 导出CSV
  • [1]

    NIE R S, LI Y F, LENG W M, et al. Deformation characteristics of fine-grained soil under cyclic loading with intermittence[J]. Acta Geotech, 2020, 15: 3041-3054 doi: 10.1007/s11440-020-00955-3

    [2]

    LEI H Y, LIU M, FENG S X, et al. Cyclic behavior of Tianjin soft clay under intermittent combined-frequency cyclic loading[J]. International Journal of Geomechanics, 2020, 20(10): 04020186. doi: 10.1061/(ASCE)GM.1943-5622.0001805

    [3]

    SAKAI A, SAMANG L, MIURA N. Partially-drained cyclic behavior and its application to the settlement of a low embankment road on silty-clay[J]. Soils and Foundations, 2003, 43(1): 33-46. doi: 10.3208/sandf.43.33

    [4]

    CAI Y Q, GU C A, WANG J, et al. One-way cyclic triaxial behavior of saturated clay: comparison between constant and variable confining pressure[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2013, 139(5): 797-809. doi: 10.1061/(ASCE)GT.1943-5606.0000760

    [5]

    HUANG J H, CHEN J, KE W H, et al. Damping ratio evolution of saturated Ningbo clays under cyclic confining pressure[J]. Soil Dynamics and Earthquake Engineering, 2021, 143: 106581. doi: 10.1016/j.soildyn.2021.106581

    [6]

    GU C, WANG J, CAI Y Q, et al. Undrained cyclic triaxial behavior of saturated clays under variable confining pressure[J]. Soil Dynamics and Earthquake Engineering, 2012, 40: 118-128. doi: 10.1016/j.soildyn.2012.03.011

    [7]

    SUN L, GU C, WANG P. Effects of cyclic confining pressure on the deformation characteristics of natural soft clay[J]. Soil Dynamics and Earthquake Engineering, 2015, 78: 99-109. doi: 10.1016/j.soildyn.2015.07.010

  • 期刊类型引用(3)

    1. 薛强,杜延军,胡黎明,詹良通,李江山. 环境土力学与工程研究进展. 土木工程学报. 2025(03): 83-112 . 百度学术
    2. 张春华,黄江东,邓正定,谢海建,邓通发. 热渗透作用下有机污染物在双人工复合衬层中的一维运移模型. 岩土工程学报. 2024(06): 1254-1262 . 本站查看
    3. 江文豪,李江山,冯晨. 考虑力学-化学荷载下压实黏土垫层中一维非线性固结与污染物运移耦合模型. 岩土工程学报. 2023(11): 2289-2298 . 本站查看

    其他类型引用(2)

图(4)  /  表(3)
计量
  • 文章访问数:  140
  • HTML全文浏览量:  26
  • PDF下载量:  52
  • 被引次数: 5
出版历程
  • 收稿日期:  2023-07-04
  • 网络出版日期:  2023-11-23
  • 刊出日期:  2023-10-31

目录

/

返回文章
返回