Experimental investigation on evolution process of suffusion in gap-graded cohesionless soil
-
摘要: 渗蚀是指土体内部细颗粒在渗流作用下克服粒间作用力逐步脱离,在粗颗粒骨架孔隙中发生迁移,并可能引起土骨架应力重分布和变形的侵蚀过程。利用自制局部孔压可测的三轴渗蚀装置,研究了细颗粒含量和相对密度对不连续级配无黏性土渗蚀行为的影响,并根据渗流沿程局部水力梯度的时空演变揭示了土体渗蚀发展的演变特征。结果表明:不连续级配无黏性土启动和破坏水力梯度随细颗粒含量/相对密度的增大而增大;相对密度越大侵蚀程度越小,且等向应力下相对密度增大到一定值后,土体会由渗流不稳定状态转变为稳定状态;土体渗蚀启动的内部表现是局部水力梯度的突变及渗流沿程上的不均匀分布。土体渗蚀导致细颗粒流失,土体孔隙比增大,在等向应力条件下,会引发体积收缩现象。Abstract: The suffusion involves selective erosion and gradual migration of fine particles through the voids of soil skeleton formed by coarse particles under seepage flow. As a result, redistribution of soil skeleton stress and deformation of soil may be induced. In this study, a series of suffusion tests are carried out using the triaxial erosion apparatus with measurable local pore pressure. The effects of the initial fine particle content and initial relative density on the suffusion of a gap-graded cohesionless soil are investigated. According to the spatial-temporal evolution of local hydraulic gradients along seepage path, the evolution process of the suffusion is revealed. Test results show that both the initiation and the failure hydraulic gradients of the gap-graded cohesionless soil increase with the increase of the fine particle content and relative density. The cumulative loss of fine particles decreases significantly with the increase of the relative density. When the relative density increases to a certain value under isotropic stress condition, the gap-graded cohesionless soil will change from an unstable state of seepage to a stable one. Additionally, the internal manifestation of suffusion initiation of soil is the mutation and uneven distribution of local hydraulic gradient along seepage path. The suffusion will cause the loss of fine particles, as well as the increase of the void ratio. Under the isotropic stress condition, the volume shrinkage is induced.
-
0. 引言
对于垃圾填埋场、危险废物处置场等污染场地的防渗阻隔问题,底部衬垫系统是污染场地防渗系统的重要组成部分,主要用于阻止和延缓渗滤液中污染物的向下扩散与迁移[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. 压实黏土衬垫中的耦合模型
1.1 计算简图
如图 1为非等温分布条件下压实黏土衬垫(CCL)中固结与污染物运移耦合模型的计算简图。图 1中,从上至下依次为堆填体、渗滤液收集和移除系统、渗滤液、CCL、水平导流层及不透水基底层。这里假定CCL始终处于饱和状态,其厚度为L,坐标系z由CCL上边界垂直向下,hw为CCL上部的渗滤液水头;C0为污染物浓度,这里假定污染物浓度保持恒定[16-17, 30]。为研究非等温分布条件对固结及污染物运移过程的影响,同时考虑到温度的传导速度较快,这里假定CCL中的温度T随深度呈线性分布,Ta和Tb分别为上部和底部的温度[12-18]。此外,考虑到污染场地在运行过程会不断堆填固废垃圾体,从而使得CCL受到的竖向应力不断增大,这里设定堆填体产生的加荷速率为Q,并记堆填总时间(即加荷时间)为tc。
1.2 基本假定
由于底部衬垫系统的水平面积较大,CCL的固结问题可近似为一维问题。已有试验研究表明,CCL的渗透系数会随温度的升高而增大[32-33],但温度变化对土体压缩曲线和压缩模量的影响较小[34-35]。基于试验研究结果,作如下基本假定:①CCL是均质、各向同性的,且处于饱和状态;②CCL中固相(土颗粒)和液相(孔隙水)均不可压缩,土体变形固结是由于液相的排出所引起;③液相的排出会导致土体的孔隙率减小,但土体的变形量远小于土体的厚度,即小变形假定成立;④液相中污染物的浓度很低,忽略污染物浓度变化对土体物理力学性质的影响;⑤土中液相的渗流符合Darcy定律;⑥温度变化不会导致CCL的压缩性改变,仅影响CCL的渗透系数。
为研究污染物在CCL中的一维运移问题,在上述基本假定的基础上,还需补充以下3点假定:①污染物在CCL中的扩散遵循Fick第二定律,土颗粒对污染物的吸附状态已达到平衡吸附状态;②污染物为单一离子或有机物;③污染物的运移方式考虑分子扩散、对流、机械弥散、吸附和热扩散。
1.3 CCL的一维固结控制方程
在CCL一维固结过程中,液相和固相的质量守恒方程有[36]:
−∂(nvtcρf)∂z=∂(nρf)∂t, (1a) −∂[(1−n)vsρs]∂z=∂[(1−n)ρs]∂t, (1b) 式中,n为CCL的孔隙率,vtc为液相中因土体固结产生的流动速度,ρf为液相的密度,vs,ρs分别为固相的移动速度和密度,t为时间。
由于固相和液相的密度保持不变,根据式(1a),(1b)有:
−∂(nvtc)∂z=∂n∂t, (2a) −∂[(1−n)vs]∂z=∂(1−n)∂t。 (2b) 根据式(2)可得
∂vs∂z+∂vdc∂z=0, (3) 式中,vdc为由于土体中因固结引起的达西流速,vdc=n(vtc−vs)。
根据达西定律,土体中液相的流动速度为
vd=vdc+vdh, (4a) vdc=−kvγw∂uex∂z, (4b) vdh=−kvγw∂uh∂z, (4c) 式中,kv为液相的渗透系数,γw为液相的重度,uex为因外荷载作用产生的超孔隙水压力,其会导致土体发生固结作用,uh为渗滤液水头产生的超孔隙水压力,其会导致液相发生流动,但不会使得土体发生固结变形[36],vdh为因渗滤液水头产生的达西流速。
根据CCL中所作的非等温分布假定可得
T(z)=Ta+Mz, (5) 式中,M为CCL中的温度梯度(℃/m),M= (Tb−Ta)/(Tb−Ta)LL。
对于液相的渗透系数kv,参考文献[37],渗透系数kv与温度T的关系可近似表达为
kv=(0.029T+0.420)⋅kv(T0=20∘C), (6) 式中,kv(T0)为温度为T0时的渗透系数。
结合式(5),(6),渗透系数kv与深度z关系为
kv(z)=kv(T0)⋅(Nz+Qa), (7) 式中,N=0.029M,N为渗透系数梯度(1/m),Qa=0.029Ta+0.420,kv(0)=kv(T0)Qa,kv(0)为CCL上部的渗透系数。
对于因渗滤液水头产生的达西流速vdh,其不产生固结效应,这里假定渗滤液水头hw保持不变,并在该水头下达到了稳定渗流状态,则有
vdh=−kv(z)γwduhdz。 (8) 对式(8)进行一定的变换,并利用式(7)的表达式进行积分可得
vdh=Nkv(T0)hwln[(N/Qa)L+1]。 (9) 在CCL中,根据液相中的渗流连续性条件可得
∂vdc∂z=mv∂σ′∂t, (10) 式中,mv为体积压缩系数,σ′为有效应力。
根据式(4b),(10)可进一步写为
∂∂z(kv(z)γw∂uex∂z)=−mv∂σ′∂t。 (11) 根据土体的有效原理可知:
σ′=σ−uex−uh−us, (12) 式中,σ为总应力,us为静水压力。
由于渗滤液水头下的渗流已达稳定渗流状态,且静水压力不随时间变化,则根据式(12)和上部荷载的加荷条件,式(11)可进一步改写为
∂∂z(kv(z)γw∂uex∂z)={mv∂uex∂t−mvQ(t<tc)mv∂uex∂t(t⩾tc) 。 (13) 式(13)即为非等温分布条件下CCL的一维固结控制方程。
1.4 CCL中污染物的一维运移方程
当CCL发生固结变形时,液相会发生流动,固相(土颗粒)也会发生相对移动。根据质量守恒定律,考虑分子扩散、对流、机械弥散、吸附和热扩散时污染物的一维运移方程可写为[12-17, 36]
∂(nC)∂t+∂[(1−n)ρsS]∂t=∂∂z(nD∂C∂z)−∂(nvtC)∂z−∂[1−n)vsρsS]∂z+∂∂z(nDT∂T∂z), (14) 式中,C(z,t)为液相中污染物的浓度,S为单位质量土颗粒所吸附的污染物质量,D为水动力弥散系数,D=De+Dmd,De为有效扩散系数,Dmd=αLvt,Dmd为机械弥散系数,αL为纵向弥散度,vt=vtc+vth,vt为液相的流动速度,vth为液相中因渗滤液水头产生的流动速度,vth=vdh/vdhnn,DT为污染物的热扩散系数。
根据污染物的一维运移机理,以垂直向下方向为正方向的污染物通量J(z,t)的表达式可写为
J=−nD∂C∂z+[nvtC+(1−n)vsρsS]−nDT∂T∂z。 (15) 对于S和C的关系,参考已有研究结果,这里采用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分别为温度T0和T时污染物的动力黏度。
De=D0(T0)nβ[1+A(T−T0)], (19) 式中,A为与有效扩散系数相关的温度系数。
将式(5)代入式(19)可得
De=D0(T0)nβ(Gz+H), (20) 式中,G=AM,H=1+A(Ta−T0)。
对于污染物的热扩散系数DT,参考文献[13~17]和式(20),其表达式为
DT=CSTDe=CSTD0(T0)nβ(Gz+H), (21) 式中,ST为索雷特系数。
根据上述关系式,式(14)可进一步改写为
∂[nC+(1−n)ρsKfCF]∂t=[D0(T0)nβ+1(Gz+H)+nαLvt].∂2C∂z2+{∂[D0(T0)nβ+1(Gz+H)+nαLvt]∂z−[nvt+(1−n)vsρsKfCF−1]+MSTD0(T0)nβ+1(Gz+H)}∂C∂z−C∂{[nvt+(1−n)vsρsKfCF−1]−MSTD0(T0)nβ+1(Gz+H)}∂z。 (22) 同理,根据关系式,污染物通量J(z,t)的表达式可进一步改写为
J=−[D0(T0)nβ+1(Gz+H)+nαLvt]∂C∂z+[nvtC+(1−n)vsρsKfCF−CMSTD0(T0)nβ+1(Gz+H)]。 (23) 1.5 初始条件和边界条件
为便于分析,这里以渗滤液水头下的渗流达到稳定时刻为初始时刻,并假定在该时刻开始堆填固废垃圾体。因此,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)。
2. 耦合模型的有限差分解
对于非等温分布条件下CCL中固结和污染物运移的耦合模型,由于耦合模型中控制方程的系数均是变化的,属于变系数偏微分方程,因而以下将采用有限差分法对耦合模型进行求解。
2.1 一维固结控制方程的有限差分解
根据式(7),对控制方程式(13)进行整理可得
Cv0(Nz+Qa)∂2uex∂z2+Cv0N∂uex∂z={∂uex∂t−Q(t<tc)∂uex∂t(t⩾tc), (30) 式中,Cv0=kv(T0)mvγw,Cv0为温度为T0时CCL的固结系数。
设Δz,Δt分别为计算空间步长和时间步长,并将空间坐标和时间坐标分别进行I等分和J2等分,则有zi=iΔz,Δz=L/I,ti=iΔt,j=0,1,2,⋯,J1,⋯ J2,Δt=tc/J1=t0/J2,t0为给定时间。因此,控制方程式(30)的Crank- Nicholson型隐式差分格式可写为
Cv0(NΔzi+Qa)(ujexi+1−2ujexi+ujexi−12Δz2+uj+1exi+1−2uj+1exi+uj+1exi−12Δz2)+Cv0Nujexi+1−ujexi−12Δz={uj+1exi−ujexiΔt−Q(j⩽J1)uj+1exi−ujexiΔ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等参数。与此同时,利用达西速度vdc,vdh可以分别确定对应液相的流动速度vtc,vth。
2.2 一维运移方程的有限差分解
参考式(13)中的有限差分方法,对式(22)展开有限差分,其Crank-Nicholson型隐式差分格式可写为
TjiCj+1i−CjiΔt=Xji(Cj+1i+1−2Cj+1i+Cj+1i−12Δz2)+Xji(Cji+1−2Cji+Cji−12Δz2)+YjiCji+1−Cji−12Δz−ZjiCji, (37) 其中,
Tji=nji+(1−nji)ρsKfCjF−1i, (38) Xji=D0(T0)njβ+1i(GΔzi+H)+njiαLvjti, (39) Yji=MSTD0(T0)njβ+1i(GΔzi+H)−[njivjti+(1−nji)vjsiρsKfCjF−1i]+12Δz{D0(T0)nji+1β+1[GΔz(i+1)+H]+nji+1αLvjti+1−D0(T0)nji−1β+1[GΔz(i−1)+H]−nji−1αLvjti−1}, (40) Zji=12Δz{[nji+1vjti+1+(1−nji+1)vjsi+1ρsKfCjF−1i+1]−MSTD0(T0)nji+1β+1[GΔz(i+1)+H]−[nji−1vjti−1+(1−nji−1)vjsi−1ρsKfCji−1F−1]+MSTD0(T0)njβ+1i−1.[GΔz(i−1)+H]} 。 (41) 相应的初始条件和边界条件可写为
C0i=0, (42) Cj0=C0, (43) CjI=0 。 (44) 对于非等温分布条件下CCL中污染物的一维运移问题,利用差分方程式(37)及相应的求解条件式(42)~(44)即可进行求解。
对于污染物通量J,其差分格式可写为
Jji=−XjiCji+1−Cji−12Δz+ΔjiCji, (45) 式中,
Δji=[njivjti+(1−nji)vjsiρsKfCjF−1i−MSTD0(T0)njβ+1i(GΔzi+H)]。 (46) 需要说明的是,对于边界处的差分格式,可利用相邻节点展开有限差分。
基于差分方程式(31),(37),可利用Matlab软件编制相应的程序,从而对CCL中考虑固结与污染物运移的耦合问题展开运算。
3. 耦合模型的验证
为验证CCL中固结与污染物运移耦合模型的正确性,以下将耦合模型的有限差分解分别与Rosanne等[19]热扩散试验结果、张文杰等[7]所求解析解、Alshawabkeh等[31]所提数值模型进行对比分析。
3.1 与热扩散试验结果的对比
Rosanne等[19]对非等温分布条件下压实黏土中溶质的一维运移问题开展了室内试验研究。在该试验中,压实黏土厚度为2 mm,孔隙率为0.6,溶质为NaCl溶液,其初始浓度为1 moL/m3。压实黏土上下边界的温度分别为12,38℃,且上下边界的浓度分别为1.21,8.83 moL/m3。参考其他文献的模拟方法[20],这里采用在24.9℃时测定得到的有效扩散系数De=6.0×10−11 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.3.2 与解析解的对比
张文杰等[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.3.3 与已有耦合模型的对比
Alshawabkeh等[31]建立了考虑小变形固结与污染物运移的耦合模型(简称Alshawabkeh耦合模型)。相比于本文所建耦合模型,Alshawabkeh耦合模型忽略了固结变形引起的相关参数变化以及机械弥散等因素的影响,且该研究并未考虑温度引起的热扩散对污染物运移过程的影响。为验证本文所提方法用于求解耦合问题的正确性,可将本文所建耦合模型的计算结果与Alshawabkeh耦合模型的计算结果进行对比。除指定参数外,这里设定相关计算参数同3.2节中的参数。需说明的是,由于Alshawabkeh耦合模型中所考虑的吸附模型为等温线性吸附模型[31],因而这里取Kd=Kf和F=1,其他条件同本文所建耦合模型时所采用的相应条件,这里不再赘述。
图 4所示为本文所建耦合模型的计算结果和Alshawabkeh耦合模型的计算结果的对比情况。从图 4可以看出,所建耦合模型计算所得污染物通量随时间的变化情况与Alshawabkeh耦合模型的计算结果十分一致,这说明所建耦合模型可用于计算土体固结与污染物运移的耦合问题。
4. 算例分析
为认识非等温分布条件下CCL中考虑固结时污染物的运移特性,除指定参数外,这里以3.2节中的参数为例,探究温度梯度、加荷速率及Freundlich吸附系数等因素对污染物在CCL中运移规律的影响。
4.1 温度梯度的影响
图 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绝对值的增大会使得污染物运移过程更早趋于稳定。
出现上述规律的原因在于:温度梯度M绝对值的增大不仅会加快热扩散过程,还会增大污染物在CCL中运移时的有效扩散系数,进而加快了污染物的运移速率,并增大底部通量。分析表明,温度梯度M会加快污染物在CCL中的运移速率,且在一定温度梯度下,污染物的底部通量可达不考虑温度梯度时底部通量的2倍以上。因此,实际工程中应考虑CCL中的温度梯度对污染物运移过程的影响[12-17]。
4.2 加荷速率的影响
图 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)中可以看出,不考虑固结效应会低估污染物的运移过程达到稳态的时间。
上述规律表明,CCL的固结效应会减慢污染物的运移速率,并减小运移时的底部通量,这与Pu等分析得到的规律一致[30],但同时,固结效应也会使得运移过程达到稳态时的浓度增大。对于运移速率减慢和底部通量减小的原因,这主要是由于固结效应会使得土体孔隙率减小,进而使得与污染物运移相关的系数减小,如有效扩散系数和热扩散系数减小。对于运移过程达到稳态时浓度增大的原因,这主要是由于当考虑固结效应时,相关参数的比值会发生变化,如对流速度与水动力弥散系数的比值会增大,因而改变了运移过程达到稳态时污染物浓度的分布规律[7, 14]。
4.3 Freundlich吸附系数的影响
图 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%。
出现上述运移规律的主要原因:CCL的吸附作用主要对污染物运移过程起着阻滞作用,因而主要影响运移速率快慢[28-30]。当运移过程达到稳态时,吸附作用对浓度分布规律和底部通量的影响可忽略。研究表明,CCL对污染物的吸附作用有利于延缓运移过程,且与不考虑吸附作用时相比,考虑吸附作用可使得运移过程达到稳态时所需的时间延长3倍及以上。因此,在选择CCL时,应考虑其吸附特性[28, 30]。
5. 结论
(1)CCL中的温度梯度会加快污染物运移速率,相同时间内污染物的浓度和底部通量会随温度梯度M绝对值的增大而增大,且在一定温度梯度下,底部通量可达不考虑温度梯度时底部通量的2倍以上。因此,实际工程中应考虑温度梯度对运移过程的影响。
(2)加荷速率Q(CL的固结效应)对污染物运移过程的影响是多方面的。一方面,Q的增大会使得土体的孔隙率减小,从而减慢了运移速率,也减小了底部通量;另一方面,Q的增大也会使得污染物运移过程达到稳态时的浓度增大。因此,实际工程中应综合考虑固结效应对运移过程的影响。
(3)Freundlich吸附系数Kf的增大会减慢污染物在CCL中的运移过程,且与不考虑吸附作用时相比,考虑吸附作用可使得污染物运移过程达到稳态时所需的时间延长3倍及以上。当污染物的运移过程达到稳态时,土颗粒对污染物的吸附作用可忽略。在选择CCL时,应考虑土体的吸附特性。
-
表 1 试验材料基本物性参数
Table 1 Physical properties of test soils
基本物
性参数细粒
含量FC/%颗粒相对质量密度Gs 不均匀系数Cu 曲率系数Cc 最大孔隙比emax 最小孔隙比emin 细粒组 100 2.65 1.45 0.96 1.306 0.706 粗粒组 0 2.65 1.16 0.98 1.053 0.697 FC20 20 2.65 19.20 15.41 0.764 0.428 FC30 30 2.65 21.69 0.12 0.748 0.378 FC40 40 2.65 22.83 0.10 0.729 0.336 表 2 试验工况
Table 2 Test program
试验系列 方案编号 FC /% 干密度ρd/
(g·cm-3)相对密度Dr 初始孔隙比e 有效围压/kPa 固结后孔隙比e0 内部稳定性判定 Kézdi K-L FC20-D8-C100 20 0.87 0.461 U U FC30-D8-C100 30 0.78 0.462 U U 系列Ⅰ FC40-D8-C100 40 1.8 0.65 0.472 100 0.466 U U FC30-D8-C100-R1 30 0.75 0.462 U U FC30-D8-C100-R2 30 0.75 0.462 U U FC30-D6-C100 1.6 0.25 0.656 0.600 U U FC30-D7-C100 1.7 0.50 0.559 0.518 U U 系列Ⅱ FC30-D8-C100 30 1.8 0.75 0.472 100 0.462 U U FC30-D9-C100 1.9 0.95 0.395 0.388 U U FC30-D9-C100-R1 1.9 0.95 0.395 0.388 U U 注:方案编号中,FC为细颗粒含量;D为干密度;C(confining pressure)为等向应力即围压数值;R为重复性试验。内部稳定性判定中,U表示内部不稳定。Kézdi法中dc15/df85 > 4为内部不稳定,其中dc15为土体粗粒组15%质量百分数所对应粒径;df85为土体细粒组85%质量百分数所对应粒径。K-L法(Kenney和Lau方法)中(H/F)min<1为内部不稳定,其中H为土体任意粒径d对应质量百分数;F为土体任意粒径d~4d之间颗粒对应质量百分数。 表 3 可重复性试验结果
Table 3 Results of repeatability tests
试验工况 iin Cv/% if Cv/% k0/(cm·s-1) Cv/% kin/(cm·s-1) Cv/% kav/(cm·s-1) Cv/% FC30-D8-C100 1.15 0.8 2.27 4.2 0.0381 4.5 0.0365 0.5 0.0378 1.8 FC30-D8-C100-R1 1.17 2.46 0.0392 0.0366 0.0375 FC30-D8-C100-R2 1.15 2.24 0.0423 0.0369 0.0389 -
[1] FOSTER M, FELL R, SPANNAGLE M. The statistics of embankment dam failures and accidents[J]. Canadian Geotechnical Journal, 2000, 37(5): 1000-1024. doi: 10.1139/t00-030
[2] Iternational Commiss on Large Dams. Internal Erosion of Existing Dams, Levees and Dikes, and Their Toundations (ICOLD Bulletin 164)[R]. Paris: Iternational Commiss on Large Dams, 2017.
[3] KÉZDI, A. Soil Physics[M]. Amsterdam: Elsevier Scientific Publishing Company, 1979.
[4] 李伟一, 钱建固, 尹振宇, 等. 间断级配砂土渗流侵蚀现象的CFD-DEM模拟[J]. 岩土力学, 2021, 42(11): 3191-3201. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202111028.htm LI Weiyi, QIAN Jiangu, YIN Zhenyu, et al. Simulation of seepage erosion in gap graded sand soil using CFD-DEM[J]. Rock and Soil Mechanics, 2021, 42(11): 3191-3201. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202111028.htm
[5] OUEIDAT M, BENAMAR A, BENNABI A. Effect of fine particles and soil heterogeneity on the initiation of suffusion[J]. Geotechnical and Geological Engineering, 2021, 39(3): 2359-2371. doi: 10.1007/s10706-020-01632-8
[6] 宋林辉, 黄强, 闫迪, 等. 水力梯度对黏土渗透性影响的试验研究[J]. 岩土工程学报, 2018, 40(9): 1635-1641. doi: 10.11779/CJGE201809009 SONG Linhui, HUANG Qiang, YAN Di, et al. Experimental study on effect of hydraulic gradient on permeability of clay[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(9): 1635-1641. (in Chinese) doi: 10.11779/CJGE201809009
[7] PACHIDEH V, MIR MOHAMMAD HOSSEINI S M. A new physical model for studying flow direction and other influencing parameters on the internal erosion of soils[J]. Geotechnical Testing Journal, 2019, 42(6): 20170301. doi: 10.1520/GTJ20170301
[8] CHANG D S, ZHANG L M. Critical hydraulic gradients of internal erosion under complex stress states[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2013, 139(9): 1454-1467. doi: 10.1061/(ASCE)GT.1943-5606.0000871
[9] 陈勇, 闵泽鑫, 夏振尧, 等. 渗流作用下粉土质砂潜蚀演化特征与预测模型[J/OL]. 土木与环境工程学报(中英文), 2023: 1-9. CHEN Yong, MING Zexin, XIA Zhenyao, et al. Evolution characteristics and prediction model on suffusion of silty-sand subjected to seepage[J/OL]. Journal of Civil and Environmental Engineering, 2023: 1-9. https://kns.cnki.net/kcms/detail/50.1218.TU.20210816.1145.003.html. (in Chinese)
[10] 刘杰, 谢定松. 砾石土渗透稳定特性试验研究[J]. 岩土力学, 2012, 33(9): 2632-2638. doi: 10.16285/j.rsm.2012.09.009 LIU Jie, XIE Dingsong. Research on seepage stability experiment of gravelly soil[J]. Rock and Soil Mechanics, 2012, 33(9): 2632-2638. (in Chinese) doi: 10.16285/j.rsm.2012.09.009
[11] 陈生水, 凌华, 米占宽, 等. 大石峡砂砾石坝料渗透特性及其影响因素研究[J]. 岩土工程学报, 2019, 41(1): 26-31. doi: 10.11779/CJGE201901002 CHEN Shengshui, LING Hua, MI Zhankuan, et al. Experimental study on permeability and its influencing factors for sandy gravel of Dashixia Dam[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(1): 26-31. (in Chinese) doi: 10.11779/CJGE201901002
[12] 陈群, 谷宏海, 何昌荣. 砾石土防渗料-反滤料联合抗渗试验[J]. 四川大学学报(工程科学版), 2012, 44(1): 13-18. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201201004.htm CHEN Qun, GU Honghai, HE Changrong. Combination seepage failure test of gravelly soil and the filter[J]. Journal of Sichuan University (Engineering Science Edition), 2012, 44(1): 13-18. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201201004.htm
[13] LIANG Y, YEH T C J, ZHA Y Y, et al. Onset of suffusion in gap-graded soils under upward seepage[J]. Soils and Foundations, 2017, 57(5): 849-860. doi: 10.1016/j.sandf.2017.08.017
[14] LIU K W, QIU R Z, SU Q, et al. Suffusion response of well graded gravels in roadbed of non-ballasted high speed railway[J]. Construction and Building Materials, 2021, 284: 122848. doi: 10.1016/j.conbuildmat.2021.122848
[15] 袁涛, 蒋中明, 刘德谦, 等. 粗粒土渗透损伤特性试验研究[J]. 岩土力学, 2018, 39(4): 1311-1316, 1336. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201804022.htm YUAN Tao, JIANG Zhongming, LIU Deqian, et al. Experiment on the seepage damage coarse grain soil[J]. Rock and Soil Mechanics, 2018, 39(4): 1311-1316, 1336. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201804022.htm
[16] KE L, TAKAHASHI A. Strength reduction of cohesionless soil due to internal erosion induced by one-dimensional upward seepage flow[J]. Soils and Foundations, 2012, 52(4): 698-711. doi: 10.1016/j.sandf.2012.07.010
[17] WANG J J, QIU Z F. Anisotropic hydraulic conductivity and critical hydraulic gradient of a crushed sandstone–mudstone particle mixture[J]. Marine Georesources & Geotechnology, 2017, 35(1): 89-97.
[18] MOFFAT R, FANNIN R J, GARNER S J. Spatial and temporal progression of internal erosion in cohesionless soil[J]. Canadian Geotechnical Journal, 2011, 48(3): 399-412. doi: 10.1139/T10-071
[19] 谷敬云, 罗玉龙, 张兴杰, 等. 基于平面激光诱导荧光的潜蚀可视化试验装置及其初步应用[J]. 岩石力学与工程学报, 2021, 40(6): 1287-1296. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202106019.htm GU Jingyun, LUO Yulong, ZHANG Xingjie, et al. A suffusion visualization apparatus based on planar laser induced fluorescence and the preliminary application[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(6): 1287-1296. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202106019.htm
[20] 梁越, 代磊, 魏琦. 基于透明土和粒子示踪技术的渗流侵蚀试验研究[J]. 岩土工程学报, 2022, 44(6): 1133-1140. doi: 10.11779/CJGE202206018 LIANG Yue, DAI Lei, WEI Qi. Experimental study on seepage erosion based on transparent soil and particle tracing technology[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(6): 1133-1140. (in Chinese) doi: 10.11779/CJGE202206018
[21] 常东升, 张利民. 土体渗透稳定性判定准则[J]. 岩土力学, 2011, 32(增刊1): 253-259. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S1045.htm CHANG Dongsheng, ZHANG Limin. Internal stability criteria for soils[J]. Rock and Soil Mechanics, 2011, 32(S1): 253-259. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2011S1045.htm
[22] 朱秦, 苏立君, 刘振宇, 等. 颗粒迁移作用下宽级配土渗透性研究[J]. 岩土力学, 2021, 42(1): 125-134. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202101014.htm ZHU Qin, SU Lijun, LIU Zhenyu, et al. Study of seepage in wide-grading soils with particles migration[J]. Rock and Soil Mechanics, 2021, 42(1): 125-134. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202101014.htm
[23] 周健, 姚志雄, 张刚. 基于散体介质理论的砂土管涌机制研究[J]. 岩石力学与工程学报, 2008, 27(4): 749-756. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200804016.htm ZHOU Jian, YAO Zhixiong, ZHANG Gang. Research on piping mechanism in sandy soils based on discrete element theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(4): 749-756. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200804016.htm
[24] KE L, TAKAHASHI A. Experimental investigations on suffusion characteristics and its mechanical consequences on saturated cohesionless soil[J]. Soils and Foundations, 2014, 54(4): 713-730.
[25] 陈锐, 谭润锵, 赵燕茹, 等. 一种用于乳胶膜开孔后的密封装置及其使用方法: CN111042097A[P]. 2020-04-21. CHEN Rui, TAN Runqiang, ZHAO Yanru, et al. Sealing Device Used after Hole Forming of Latex Film and Application Method of Sealing Device: CN111042097A[P]. 2020-04-21. (in Chinese)
[26] KENNEY T C, LAU D. Internal stability of granular filters[J]. Canadian Geotechnical Journal, 1985, 22(2): 215-225.
[27] CHANG D S, ZHANG L M, et al. A stress-controlled erosion apparatus for studying internal erosion in soils[J]. Geotechnical Testing Journal, 2011, 34(6): 103889.
[28] CHEN C, ZHANG L M, ZHU H. A photographic method for measuring soil deformations during internal erosion under triaxial stress conditions[J]. Geotechnical Testing Journal, 2018, 41(1): 20170031.
[29] RICHARDS K S, REDDY K R, et al. True triaxial piping test apparatus for evaluation of piping potential in earth structures[J]. Geotechnical Testing Journal, 2010, 33(1): 102246.
[30] LIANG Y, ZENG C, WANG J J, et al. Constant gradient erosion apparatus for appraisal of piping behavior in upward seepage flow[J]. Geotechnical Testing Journal, 2017, 40(4): 630-642.
[31] DENG G, ZHANG L L, CHEN R, et al. Experimental investigation on suffusion characteristics of cohesionless soils along horizontal seepage flow under controlled vertical stress[J]. Frontiers in Earth Science, 2020, 8: 195.
-
期刊类型引用(3)
1. 薛强,杜延军,胡黎明,詹良通,李江山. 环境土力学与工程研究进展. 土木工程学报. 2025(03): 83-112 . 百度学术
2. 张春华,黄江东,邓正定,谢海建,邓通发. 热渗透作用下有机污染物在双人工复合衬层中的一维运移模型. 岩土工程学报. 2024(06): 1254-1262 . 本站查看
3. 江文豪,李江山,冯晨. 考虑力学-化学荷载下压实黏土垫层中一维非线性固结与污染物运移耦合模型. 岩土工程学报. 2023(11): 2289-2298 . 本站查看
其他类型引用(2)
-
其他相关附件