One-dimensional analytical model for contaminant transport through CCL under thermal diffusion and its application
-
摘要: 建立了污染物在对流、分子扩散和热扩散作用下压实黏土衬垫（CCL）的一维运移解析模型，并获得解析解。无量纲设计曲线表明对流作用对污染物在CCL中运移的影响至关重要；然而随着热扩散作用的增强，对流作用的影响会有所减弱。参数敏感性分析结果表明，当渗滤液水头达到3 m时，10 a的底部浓度和通量分别是无对流情况下的3.5~4.9倍和5.9~15.1倍；当热扩散作用较强（M=-5）时，10 a底部浓度和通量是无热扩散作用下的2.6倍和3.5倍；温度升高会增大土的渗透系数，从而影响污染物在CCL中的运移，在衬垫系统设计过程中应予以考虑。提出了衬垫厚度设计简化计算方法，并以西安江村沟填埋场为例，对CCL进行了简化设计；结果表明：若要保证渗滤液中的Cl-和As击穿CCL的时间大于50 a，则CCL的厚度分别需要不小于11.16 m和1.75 m。Abstract: A one-dimensional analytical model as well as the analytical solution for the contaminant transport in compacted clay liner (CCL) is proposed. The dimensionless design curves show that the advection has significant effects on the contaminant transport in CCL. However, with the increases of thermal diffusion, the effects of the advection will decrease. The results of parameter sensitive analysis show when the leachate head reaches 3 m, the 10-year bottom concentration and flux increase by the factors of 3.5~4.9 and 5.9~15.1, respectively, compared with those without the advection. When the thermal diffusion is great enough (M=-5), the bottom concentration and flux increase by the factors of 2.6 and 3.5, respectively, compared with those without the thermal diffusion. The effects of the temperature on the permeability coefficient of CCL should be considered in the design of landfill clay liners. A simplified method for determination of the thickness of the landfill liner is proposed. The simplified method is used to design the liner system at the Jiangcungou Landfill site in Xi 'an. When considering the chloride (Cl-) and arsenic (As) as the index contaminants, the breakthrough time will be longer than 50 years in the case with the thickness of CCL larger than 11.16 m and 1.75 m, respectively.
-
Keywords:
- analytical model /
- temperature /
- thermal diffusion /
- leachate /
- simplified design
-
0. 引言
近年来，中国大力发展南海地区基础设施建设，钙质砂作为沿海地区的主要岩土体受到广泛使用。有别于陆源石英砂，钙质砂作为难溶类碳酸盐（多为碳酸钙）物质的集合体，具有形状不规则、内孔隙高、颗粒易破碎、棱角度高等性质[1-3]。国内外学者对于钙质砂的力学特性开展了大量研究：张晨阳等[4]分析了细粒含量对钙质砂力学特性的影响，认为钙质砂峰值偏应力与剪胀性随细粒含量的增加而逐渐减小；黄宏翔等[5]得出在相同条件下，钙质砂的抗剪强度和残余强度均大于石英砂；Rui等[6]发现不规则颗粒形状的钙质砂表现出较慢的累积孔隙压力速率，显示出更高的抗液化性。对于钙质砂的颗粒破碎特性，研究表明在低围压情况下就会发生颗粒破碎，颗粒形状对其力学性质的影响较大，在高围压情况下钙质砂颗粒破碎情况加剧，颗粒破碎对其力学性质影响愈加显著[7]。Coop等[8]指出钙质砂的颗粒破碎程度受轴向应力大小及初始颗粒形状影响；汪轶群等[9]在三轴试验中发现随着初始围压的增大，钙质砂颗粒破碎程度加大，土样整体剪胀趋势减小，且破碎引起的能量耗散增加。
上述研究为以钙质砂的颗粒特性为切入点，对其力学性质进行了研究。但主要以钙质砂作为纯净砂进行分析，而仅考虑某一粒径、某一种类砂土多与实际工程情况不符，存在一定的局限性。中国南海地区广泛分布源于岛礁中的碳酸盐土体和沿岸地区的非碳酸盐类土所结合的碳酸盐混合土[10]，其力学性质势必由两种土体共同决定。李小刚等[11]通过三轴试验发现，混合砂峰值强度和相对破碎率随石英砂含量的增加而减小，黏聚力及峰值内摩擦角随碳酸盐含量的增加而增加。Lade等[12]指出土体中粉土、粉砂等细颗粒的含量是海底斜坡、土坝和路堤等产生液化破坏的主要原因。含细颗粒土体的试验同样表明细颗粒含量、粒径分布和颗粒形状特征是影响砂土抗剪强度和液化反应的主要因素[13-15]，且细颗粒含量与粉土特性对砂土静态液化行为的影响存在耦合关系[16]。此外，砂土中不同颗粒特征的粗颗粒和细颗粒相互作用，对混合土体的整体力学性质也产生显著影响[17-18]。
综上所述，细颗粒对钙质混合土体力学性质影响显著，而现有研究中多数仅考虑钙质细粒含量对钙质砂力学特性的影响，在对钙质混合砂的研究中，细颗粒特征不同对混合砂的影响程度亦不同。基于此，本文采用角形石英粉和圆形玻璃珠作为细颗粒，分别与钙质砂混合，对钙质混合砂进行不同轴向压力条件下的环剪试验，探究不同细颗粒含量及形状对钙质混合砂峰值强度、残余强度和轴向变形的影响，对钙质混合砂的颗粒破碎特征进行研究分析，进一步探讨颗粒破碎特征对钙质混合砂性质的影响规律，从而为南海地区工程建设提供参考。
1. 试验材料及方案
1.1 试验材料
试验所用钙质砂取自南海区域某岛礁附近。试验前先将砂样烘干，采用标准筛进行筛分，得到粒径范围为0.1~1 mm的钙质砂颗粒。玻璃珠因具有与土体相似的矿物学特征[19]，多用于模拟砂土，其形状接近理想的球形；石英粉则更具棱角和不规则性。因此，采用粒径小于0.075 mm的石英粉及玻璃珠模拟混合砂中不同形状的石英细颗粒。3类基本砂样的物理指标参数见表 1，颗粒级配曲线见图 1。按试样总质量的0%，10%，20%，30%，40%，将石英粉、玻璃珠细颗粒分别与钙质砂均匀混合得到9种试样：含石英粉试样，由J10，J20，J30，J40表示；含玻璃珠试样，由Y10，Y20，Y30，Y40表示；0表示纯钙质砂试样。
表 1 基本砂样的物性指标Table 1. Physical properties of basic sand samples试样 圆度 宽长比 D10/
μmD50/
μmD60/
μm不均匀系数Cu 曲率系数Cc 颗粒相对质量密度 钙质砂 0.517 0.541 374.5 637.7 643.5 1.72 1.28 2.75 玻璃珠 0.924 0.893 35.1 50.8 58.6 1.67 0.92 2.64 石英粉 0.635 0.746 27.2 45.7 51.3 1.89 1.05 2.66 注：①圆度为单颗粒投影面积与等周长圆面积之比；宽长比为单颗粒最小Feret直径与最大Feret直径之比。②圆度、宽长比取采样（1000粒）均值。 ![]()
图 1 3种砂样颗粒级配曲线Figure 1. Grain-size distribution curves of three kinds of sand particles为定量化表示土体形态特征，多通过颗粒形状扫描仪器，根据颗粒轮廓图像建立二维形状参数（圆度、宽长比、粗糙度等），能够较好地表述颗粒形状与其力学性质之间的关系[20-22]。对于二元颗粒混合物，参考Yang等[23]提出的组合形状参数概念，认为二元混合物组合形状参数由各颗粒材料的形状特征值及其对应的质量占比决定，组合圆度、组合宽长比越接近于1，表明颗粒整体越接近于圆形、越规整。如钙质混合砂J30的组合圆度、组合宽长比计算分别为
Sc−J30=Sc−J×Fc+Sc−G×(1−Fc)=0.635×30%+0.517×70%=0.552, (1) Ar−J30=Ar−J×Fc+Ar−G×(1−Fc)=0.746×30%+0.541×70%=0.603 (2) 图像颗粒分析仪如图 2所示，扫描基本砂样轮廓图如图 3所示，混合砂样特性指标如表 2所示。
表 2 混合砂样特性指标Table 2. Characteristics of mixed sand samples砂样 相对密实度Dr/% 最大孔隙比emax 最小孔隙比emin 组合圆度Sc 组合宽长比Ar 0 60 1.109 0.812 0.517 0.541 J10 60 1.010 0.707 0.529 0.562 J20 60 0.893 0.553 0.541 0.582 J30 60 0.794 0.467 0.552 0.603 J40 60 0.743 0.430 0.564 0.623 Y10 60 0.977 0.724 0.558 0.541 Y20 60 0.856 0.605 0.598 0.576 Y30 60 0.744 0.512 0.639 0.611 Y40 60 0.685 0.491 0.680 0.647 1.2 试验方案
试验采用HJ-1型环剪仪，其试样尺寸为外径100 mm，内径60 mm，高20 mm。天然状态下土体的相对密实度在20%~80%，且天然密实度多为中密状态[24]，故考虑实际工程情况，控制钙质混合砂的相对密实度Dr =60%（中密），进行不同轴向压力条件下的环剪试验。试验步骤如下：①采用漏斗法将试样装入容器。首先，将环剪仪的上、下剪切盒固定在一起，然后将试样缓慢倒入漏斗之中，通过漏斗将砂样以零高度的落差缓慢倒入盒中，最后将试样抚平，使其高度一致。②采用浸水法饱和试样。试样装好后倒入蒸馏水，使其完全浸没，浸泡4 h确保试样饱和度达到要求。浸泡完成后将仪器组装完成。③逐级施加100，200，300 kPa的轴向压力，待固结完成后，以1°/min的剪切速率进行环向剪切试验。
2. 试验结果分析
2.1 剪切应力-剪切位移曲线分析
试验结束后，得到不同细颗粒含量及形状条件下的剪切应力-剪切位移曲线，如图 4所示。由图 4（a）可得，纯钙质砂的剪切应力随着剪切位移的增加而逐渐增加，呈近线性增长，随后，剪切应力的增长趋势变缓。当轴向压力为100，200，300 kPa时，剪切应力分别在剪切位移为10.5，11.9，13.6 mm时达到峰值强度，其值分别为100.4，177.9，251 kPa；剪切位移分别在28，30，31 mm处得到土样的残余强度，其值分别为88.23，159.4，213.51 kPa。由各应力-位移图线得：各组试样的抗剪强度均随轴向压力的增加而增大；当剪切位移相同时，轴向压力越大，试样的剪切应力越大；同时，随轴向压力的增加，各组试样峰值强度出现延后性，且各组试样的应力-位移曲线均呈现软化特性。
对不同细颗粒含量及形状条件下钙质混合砂的峰值强度进行统计并绘制图 5。如图 5（a）所示，随着细颗粒含量的增加，混合砂峰值强度逐渐降低。与纯钙质砂相比，细颗粒含量为10%时，混合砂的峰值强度略小于纯钙质砂，其中，含石英粉混合砂峰值强度最大降幅为12%，含玻璃珠混合砂峰值强度最大降幅为16%，此时细颗粒填充土体骨架的孔隙部分，土体骨架仍由钙质砂颗粒组成；细颗粒含量为20%时，含石英粉混合砂峰值强度最大降幅为17%，含玻璃珠混合砂峰值强度最大降幅为29%；细颗粒含量为30%时，峰值强度持续降低，含石英粉混合砂峰值强度最大降幅为22%，含玻璃珠混合砂峰值强度最大降幅扩大至36%；当细颗粒含量增加至40%，峰值强度发生大幅度降低，含石英粉混合砂峰值强度最大降幅为30%，含玻璃珠混合砂峰值强度最大降幅则高达47%，表明此时土体骨架由钙质砂颗粒承担转变为钙质砂、细颗粒共同承担。此外，如图 5（b），（c）所示，峰值强度受细颗粒形状影响，随组合圆度和组合宽长比增加而逐渐减小，含玻璃珠混合砂出现更大降幅。结合剪切应力-剪切位移曲线（图 4）中含玻璃珠混合砂的曲线波动更为明显，表明细颗粒形状的改变影响了骨架颗粒间的接触状态以及力链传递，进而降低土体骨架的稳定性，圆形细颗粒的加入使得土体骨架的稳定性更差。各组试样峰值强度及下降幅度如表 3所示。
![]()
图 5 不同细颗粒含量及形状下的峰值强度Figure 5. Peak intensities at different contents and shapes of fine particles表 3 不同轴向压力下试样的峰值强度及下降幅度Table 3. Peak strengths and decrease amplitudes of samples under different axial pressures轴压 类别 0 J10 J20 J30 J40 Y10 Y20 Y30 Y40 100 kPa 峰值强度/kPa 100.4 93.6 90.7 78.0 70.2 84.8 71.7 64.3 53.6 降幅/% — 7 10 22 30 16 29 36 47 200 kPa 峰值强度/kPa 177.9 156.5 148.7 142.3 139.4 148.7 133.6 120.9 105.3 降幅/% — 12 16 20 22 16 25 32 41 300 kPa 峰值强度/kPa 251.0 241.3 208.1 198.4 188.6 232.0 201.3 185.2 158.9 降幅/% — 4 17 21 25 8 20 26 37 2.2 不同细粒含量及形状下的强度指标特征
图 6为混合砂内摩擦角随组合形状参数的变化规律。纯钙质砂的内摩擦角正切值为0.74（φ=36.4°），混合砂内摩擦角正切值随细颗粒含量、组合圆度及组合宽长比增加均表现线性减小趋势。在细颗粒含量为40%时，组合宽长比和组合圆度达到最大值，混合砂的整体颗粒形状趋于规整，含石英粉混合砂内摩擦角正切值达到最小值0.57（φ=29.7°），其降幅为23%；含玻璃珠混合砂的内摩擦角正切值为0.53（φ= 27.8°），其降幅为28%，此时表现出更小的抗剪强度。原因在于，细颗粒的填充降低了土体骨架的稳定性，在粗颗粒间起到润滑的作用，降低了颗粒间的摩擦，因而内摩擦角逐渐减小；玻璃珠因颗粒形状圆滑，随着细颗粒含量的增加，混合砂的组合形状参数逐渐接近圆形颗粒的形状参数，颗粒之间运动产生的逐渐摩擦力较小，而石英粉的颗粒形状较不规则，在颗粒间的相对运动中，不利于滚动，需要较多的力进行滑动，故含玻璃珠混合砂的内摩擦角呈现更大的降幅。
![]()
图 6 不同细颗粒含量及形状下的内摩擦角Figure 6. Angles of internal friction at different contents and shapesof fine particles钙质砂为无黏性土，但由于具有不规则的颗粒形状，颗粒之间的嵌固咬合作用产生了不可忽视的表观黏聚力[25]。图 7为混合砂表观黏聚力随颗粒含量及组合形状参数的变化规律。试验发现，随着细颗粒含量的增加，混合砂表观黏聚力由纯钙质砂时的28 kPa开始逐渐降低，含石英粉混合砂的表观黏聚力降幅较小，最终降至17.2 kPa，而含玻璃珠的混合砂表观黏聚力呈现出较大的下降趋势，最终仅为1.2 kPa。其原因是，随着细颗粒含量的增加，钙质砂颗粒所形成的骨架逐渐被细颗粒填充，减少了钙质砂颗粒之间的接触，进而降低了钙质砂颗粒之间的嵌固咬合作用；当细颗粒含量为40%时，土体骨架由钙质砂及细颗粒共同承担，由于外形更接近于球形的玻璃珠细颗粒并不具有咬合能力，导致含玻璃珠混合砂的表观黏聚力大幅下降。
![]()
图 7 不同细颗粒含量及形状下的表观黏聚力Figure 7. Apparent cohesions at different contents and shapes of fine particles2.3 不同细粒含量及形状下的残余强度特征
为定量化表述钙质混合砂的残余强度特性，采用软化系数IB分析混合砂在剪切试验后的应变软化程度。IB越大，说明砂样的强度在峰值后衰减程度越大，软化越明显，试样的变形稳定性越低。其公式为[26]
IB=σdp−σdrσdp (3) 式中：σdp为剪切破坏时的峰值强度（kPa）；σdr为剪切破坏后残余强度（kPa）。
图 8为不同细颗粒含量和形状条件下混合砂的软化系数。随着轴向压力的增大，各组混合砂的软化系数均呈现先减小后增加的趋势，轴向压力为200 kPa时，各组试样软化系数最小，试样相对稳定。细颗粒含量为10%时软化系数达到最大值，混合砂的强度在峰值后变化最大，软化程度明显。细颗粒含量超过10%后，混合砂的软化系数随着细颗粒含量的增加而减小，且在细颗粒含量为40%时软化系数达到最小值。含石英粉混合砂的软化系数明显小于含玻璃珠混合砂，表明含石英粉混合砂的变形稳定性较高，剪切应力在峰值强度后的衰减较少。混合砂中玻璃珠含量为10%时，软化最为明显，软化系数达到最大值，其数值均在0.14以上，其原因主要是由于细颗粒含量的增加，充填在钙质砂颗粒所组成骨架的孔隙中，而玻璃珠因其颗粒形状较为圆滑，在颗粒运动中起到了润滑作用，使得含玻璃珠的钙质砂混合砂受到较小的力时就能发生颗粒运动，因此试样的变形稳定性较低，在强度达到峰值后有着较大的衰减程度。与张晨阳等[4]研究结果类似，在细粒含量为10%时稳定性最差，后续随着细颗粒含量的不断增加，钙质砂的软化系数逐渐降低，试样稳定性逐渐提高。
![]()
图 8 不同轴向压力条件下试样的脆性指标Figure 8. Brittleness indexes of samples at different axial pressures2.4 不同细粒含量及形状下的变形特征
钙质混合砂在环剪过程中的轴向位移变化可在一定程度上反映出其剪胀特征。图 9为轴向压力100 kPa时环剪试验的轴向位移变形结果。从图 9中可以看出，纯钙质砂在剪切开始时出现小幅度剪缩，随着剪切位移不断增加，出现明显的剪胀现象，最大剪胀量为-0.438 mm。由图 9（a）可知，当含石英粉的混合砂中细颗粒含量分别为10%，20%，30%，40%时，最大剪缩量为0.181，0.203，0.26，0.536 mm。表明随着细颗粒含量的增加，混合砂的变形特征逐渐由剪胀状态转变为剪缩状态。由图 9（b）可知，当含玻璃珠混合砂中细颗粒含量为10%，20%，30%，40%时，最大剪缩量为分别为0.261，0.331，0.421，0.604 mm。表明同一围压下随着细颗粒含量的增加，试样的剪胀性逐渐减小，剪缩状态逐渐明显，当细颗粒含量为10%~30%时，两种不同细颗粒的混合砂的轴向变形都较为接近，相差幅度较小，当细颗粒含量增加至40%时，其轴向变形的增长趋势明显增加，且含玻璃珠的混合砂的剪缩量比含石英粉的混合砂具有明显的增长趋势，且剪缩特征更加明显。
为研究不同轴向压力条件下混合砂的轴向变形情况，对细颗粒含量为20%和40%的混合砂的轴向变形进行汇总分析，如图 10所示。由图 10（a）可知，细颗粒含量为20%的混合砂在轴向压力为100 kPa时，轴向位移先逐渐增大，随后发生小幅度的减小，最后趋于稳定，两种混合砂的最大剪缩量分别为0.203，0.262 mm；当轴向压力为200 kPa时，混合砂均呈现出明显的剪缩特征，剪缩量随剪切位移增加而增加，两种混合砂的最大剪缩量分别为0.41，0.662 mm；轴向压力为300 kPa时，混合砂均呈现出较大的剪缩特征，随着剪切位移的增加，两种混合砂的最大剪缩量分别为0.723，0.834 mm。如图 10（b）所示，细颗粒含量为40%时混合砂均表现出明显的剪缩特征。随着剪切位移的增加，剪缩量呈现出不断上升的趋势，且随着轴向压力的增加，两种混合砂的最大剪缩量分别达到0.536，0.604，0.718，0.843，0.974，1.181 mm。表明随着轴向压力不断增加，混合砂的最大剪缩量逐渐增加，呈现出明显的剪缩特征，是因为随着轴向压力不断增加，试样在剪切过程中发生较大程度的颗粒破碎现象，大颗粒逐渐减少，致使混合砂的剪缩现象显著。
![]()
图 10 细颗粒含量为20%和40%的钙质混合砂的轴向变形−剪切位移关系曲线Figure 10. Axial deformation-shear displacement curves for calcareous mixed sand with fine particles contents of 20% and 40%3. 颗粒破碎特征
钙质混合砂特殊的剪切特征主要是由于其颗粒特征和颗粒破碎等多重因素共同作用的结果，为了对不同细颗粒含量及形状条件下混合砂的剪切过程进行细观分析，需建立细颗粒含量及形状与颗粒破碎之间的联系。对剪切后试样进行颗粒级配试验，得到剪切前后各组试样的颗粒级配曲线，采用Hardin[27]所提出的相对破碎指标Br进行颗粒破碎的定量化评价，相对颗粒破碎率的计算原理如图 11所示，其计算公式如下：
Br=BtBp0。 (4) 式中：Bp0为初始破碎势，即试样试验前级配曲线的积分；Bt为总破碎势，即试验前后的整体破碎势之差。
图 12为不同细颗粒含量及形状条件下钙质混合砂的破碎关系图。由图 12可知，各组混合砂的相对破碎率均小于3%。在轴向压力为100，200，300 kPa条件下，试样的最大破碎率分别为0.9%，1.3%，2.4%。表明随轴向压力的不断增加，更多的钙质砂颗粒发生破碎，混合砂的颗粒破碎度呈现出不断上升的趋势。在细颗粒含量0~40%的增长过程中，当轴向应力为100 kPa时，两种混合砂的相对破碎率由0.9%分别减小至0.41%和0.29%；当轴向应力为200 kPa时，两种混合砂的相对破碎率由1.3%分别减小至0.87%和0.78%；当轴向应力为300 kPa时，两种混合砂的相对破碎率由2.4%分别减小至1.88%和1.78%，该现象是因为细颗粒含量的增加降低了钙质砂之间的咬合程度，颗粒破碎率也随之减小。与含玻璃珠混合砂的相对破碎率相比，含石英粉混合砂的相对破碎率更高，二者差值在0.07%~0.16%，其主要原因是石英粉的颗粒形状与玻璃珠相比更为不规则，棱角度较高，对于钙质砂颗粒之间的咬合降低程度较小，致使含石英粉的钙质混合砂有着较高的颗粒破碎率。
![]()
图 12 不同细颗粒含量及形状条件下钙质混合砂的破碎关系图Figure 12. Crushing relationship of calcareous mixed sand under different contents and shapes of fine particles图 13为钙质混合砂的相对破碎率随组合形状参数变化趋势图。由图 13可知，混合砂的相对破碎率随着细颗粒含量的增加不断减小，与组合圆度及组合宽长比均表现出线性负相关关系，表明当混合砂整体的组合形状参数越偏离球形颗粒，颗粒轮廓呈现枝状时，其破碎程度越高。细颗粒含量的增加使得混合砂整体的组合形状参数向球形颗粒靠近，包围在钙质砂周围细颗粒增加，减少了钙质砂颗粒之间的接触，起到润滑作用，使得颗粒破碎程度不断降低。玻璃珠含量的不断增加，表现出更为明显的润滑作用，且混合砂整体的组合形状参数更加趋近球形颗粒的形状参数，含玻璃珠混合砂的颗粒破碎程度明显较小。
![]()
图 13 相对破碎率与组合形状参数关系图Figure 13. Relationship between relative crushing rate and composite shape parameters4. 结论
通过对不同细颗粒形状及含量条件下的钙质混合砂进行环剪试验，分析细颗粒形状及含量间的差异对钙质混合砂力学特性及破碎特性的影响，得到以下5点结论。
（1）钙质混合砂的应力-位移曲线均呈现软化特性，含圆形玻璃珠混合砂的峰值强度明显小于含角形石英粉混合砂的峰值强度，随着细颗粒含量不断增加，混合砂的峰值强度逐渐降低，且含玻璃珠混合砂峰值强度有着更大的降低趋势。
（2）钙质混合砂的强度指标随细颗粒含量、组合圆度、组合宽长比增加均呈现减小趋势。由于玻璃珠表面圆滑，致使含玻璃珠混合砂表现出更低的内摩擦角和表观黏聚力。
（3）钙质混合砂的软化系数随轴向压力的增加先减小后增大，在200 kPa时达到最小值，并随细颗粒含量增加表现出先增大后减小的变化趋势，当细颗粒含量为10%时，软化系数最大，当细颗粒含量超过10%后，混合砂的软化系数逐渐减小，且含石英粉试样的软化系数明显低于含玻璃珠试样。
（4）纯钙质砂表现出明显的剪胀，且随细颗粒含量和轴向压力增加，混合砂的剪缩量逐渐增大，当细颗粒含量达到40%时呈现出明显的剪缩，含玻璃珠混合砂的剪缩增量更大。
（5）在不同轴压下，钙质混合砂的相对破碎率与组合圆度及组合宽长比均呈线性负相关，且轴压的增加会造成相对破碎率的增加。随细颗粒含量的增加，混合砂组合形状参数趋于球形形状参数，颗粒破碎程度减小，含玻璃珠混合砂具有更小的相对破碎率。
-
![]()
图 2 本文解和数值法的结果对比
Figure 2. Comparison between proposed analytical solutions and numerical results
![]()
图 4 渗滤液水头对污染物在CCL中运移的影响(ST=1×10-2K-1，ΔT=-40 K)
Figure 4. Effects of leachate head on contaminant transport in CCL (ST=1×10-2 K-1, ΔT=-40 K)
![]()
图 5 渗滤液水头对污染物在CCL中运移的影响(ST=5×10-2K-1，ΔT=-100 K)
Figure 5. Effects of leachate head on contaminant transport in CCL (ST=5×10-2 K-1, ΔT=-100 K)
![]()
图 6 热扩散作用对污染物在CCL中运移的影响(hw=0.3 m)
Figure 6. Effects of thermal diffusion on contaminant transport in CCL (hw=0.3 m)
![]()
图 7 热扩散作用对污染物在CCL中运移的影响(hw=3 m)
Figure 7. Effects of thermal diffusion on contaminant transport in CCL(hw=3 m)
![]()
图 8 不同温度下渗透系数对污染物在膨润土中运移的影响
Figure 8. Effects of permeability coefficient on contaminant transport in bentonite liner under different temperatures
![]()
图 9 污染物击穿CCL衬垫时厚度随渗滤液水头变化的线性拟合结果
Figure 9. Linear fitting results of relationship between thickness of CCL and leachate head based on breakthrough of contaminant
表 1 污染物在CCL中运移的参数取值表
Table 1 Parameters of contaminant transport in CCL
参数 符号 取值 渗滤液中污染物浓度 C0/(mg·L-1) 100 渗滤液水头 hw/m 0.3/3.0 有效扩散系数 D*/(m2·s-1) 4×10-10 黏土厚度 L /m 2.0 索雷特系数 ST /K-1 3×10-2 黏土层孔隙度 n 0.4 黏土层渗透系数 k /(m·s-1) 1×10-9 阻滞因子 Rd 1.0 黏土层两端温差 ΔT /K -50 
下载: 导出CSV
表 2 不同温度污染物在膨润土中运移的参数取值
Table 2 Transport parameters for contaminant in bentonite liner under different temperatures
膨润土平均温度
T0/℃膨润土渗透系数
k/(m·s-1)渗滤液水头
hw/m有效扩散系数
D*/(m2·s-1)膨润土厚度
L/m30 2.12×10-13 0.3 3.5×10-12 0.1 60 4.29×10-13 90 6.48×10-13
下载: 导出CSV
-
[1] 陈云敏. 环境土工基本理论及工程应用[J]. 岩土工程学报, 2014, 36(1): 1-46. doi: 10.11779/CJGE201401001 CHEN Yunmin. A fundamental theory of environmental geotechnics and its application[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(1): 1-46. (in Chinese) doi: 10.11779/CJGE201401001
[2] LEI M F, LIN D Y, LIU J W, et al. Modified chloride diffusion model for concrete under the coupling effect of mechanical load and chloride salt environment[J]. AIP Advances, 2018, 8(3): 035029. doi: 10.1063/1.5027540
[3] YU C, LIU J F, MA J J, et al. Study on transport and transformation of contaminant through layered soil with large deformation[J]. Environmental Science and Pollution Research, 2018, 25(13): 12764-12779. doi: 10.1007/s11356-018-1325-7
[4] PU H F, FOX P J, SHACKELFORD C D. Assessment of consolidation-induced contaminant transport for compacted clay liner systems[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 142(3): 04015091. doi: 10.1061/(ASCE)GT.1943-5606.0001426
[5] 田改垒, 张志红. 考虑热效应的污染物在土中扩散、渗透和固结耦合模型[J]. 岩土工程学报, 2022, 44(2): 278-287. doi: 10.11779/CJGE202202009 TIAN Gailei, ZHANG Zhihong. The coupled model of contaminant diffusion, osmosis and consolidation in soil considering thermal effect[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(2): 278-287. (in Chinese) doi: 10.11779/CJGE202202009
[6] YAN H X, XIE H J, WU J W, et al. Analytical model for transient coupled consolidation and contaminant transport in landfill liner system[J]. Computers and Geotechnics, 2021, 138: 104345. doi: 10.1016/j.compgeo.2021.104345
[7] YU C, WANG H, FANG D F, et al. Semi-analytical solution to one-dimensional advective-dispersive-reactive transport equation using homotopy analysis method[J]. Journal of Hydrology, 2018, 565: 422-428. doi: 10.1016/j.jhydrol.2018.08.041
[8] STOPPIELLO M G, LOFRANO G, CAROTENUTO M, et al. A comparative assessment of analytical fate and transport models of organic contaminants in unsaturated soils[J]. Sustainability, 2020, 12(7): 2949. doi: 10.3390/su12072949
[9] DING X H, FENG S J, ZHENG Q T, et al. A two-dimensional analytical model for organic contaminants transport in a transition layer-cutoff wall-aquifer system[J]. Computers and Geotechnics, 2020, 128: 103816. doi: 10.1016/j.compgeo.2020.103816
[10] 谢海建. 成层介质污染物的运移机理及衬垫系统防污性能研究[D]. 杭州: 浙江大学, 2008. XIE Haijian. A Study on Contaminant Transport in Layered Media and the Performance of Landfill Liner Systems[D]. Hangzhou: Zhejiang University, 2008. (in Chinese)
[11] 陈云敏, 谢海建, 柯瀚, 等. 层状土中污染物的一维扩散解析解[J]. 岩土工程学报, 2006, 28(4): 521-524. doi: 10.3321/j.issn:1000-4548.2006.04.018 CHEN Yunmin, XIE Haijian, KE Han, et al. Analytical solution of contaminant diffusion through multi-layered soils[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(4): 521-524. (in Chinese) doi: 10.3321/j.issn:1000-4548.2006.04.018
[12] 张文杰, 赵培, 贾文强. 一维对流-扩散试验各种边界条件及其统一形式解析解[J]. 岩土力学, 2015, 36(10): 2759-2764. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201510003.htm ZHANG Wenjie, ZHAO Pei, JIA Wenqiang. Boundary conditions of one-dimensional convection- diffusion column tests and unified analytical solution[J]. Rock and Soil Mechanics, 2015, 36(10): 2759-2764. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201510003.htm
[13] HARRIS A P F, MCDERMOTT C I, KOLDITZ O, et al. Modelling groundwater flow changes due to thermal effects of radioactive waste disposal at a hypothetical repository site near Sellafield, UK[J]. Environmental Earth Sciences, 2015, 74(2): 1589-1602. doi: 10.1007/s12665-015-4156-6
[14] PENG M Q, FENG S J, CHEN H X, et al. Analytical model for organic contaminant transport through GMB/CCL composite liner with finite thickness considering adsorption, diffusion and thermodiffusion[J]. Waste Management, 2021, 120: 448-458. doi: 10.1016/j.wasman.2020.10.004
[15] XIE H J, ZHANG C H, SEDIGHI M, et al. An analytical model for diffusion of chemicals under thermal effects in semi-infinite porous media[J]. Computers and Geotechnics, 2015, 69: 329-337. doi: 10.1016/j.compgeo.2015.06.012
[16] THOMAS H R, SEDIGHI M, VARDON P J. Diffusive reactive transport of multicomponent chemicals under coupled thermal, hydraulic, chemical and mechanical conditions[J]. Geotechnical and Geological Engineering, 2012, 30(4): 841-857. doi: 10.1007/s10706-012-9502-9
[17] SEDIGHI M, THOMAS H, VARDON P. Modelling thermal impacts on reactive transport processes related to multicomponent chemicals in compacted clays[C]//2nd International Symposium on Computational Geomechanics(ComGeo Ⅱ). Dubrovnik, 2011: 538-546.
[18] THOMAS H R, SEDIGHI M. Thermal Effects on Chemical Diffusion in Multicomponent Ionic Systems[M]//Engineering Geology for Society and Territory-Volume 6. Cham: Springer International Publishing, 2014: 525-528.
[19] FOOSE G J. Transit-time design for diffusion through composite liners[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2002, 128(7): 590-601. doi: 10.1061/(ASCE)1090-0241(2002)128:7(590)
[20] CARSLAW H S, JAEGER J C. Conduction of Heat in Solids[M]. 2d ed. Oxford: Clarendon Press, 1959.
[21] 张春华, 吴家葳, 陈赟, 等. 基于污染物击穿时间的填埋场复合衬垫厚度简化设计方法[J]. 岩土工程学报, 2020, 42(10): 1841-1848. doi: 10.11779/CJGE202010009 ZHANG Chunhua, WU Jiawei, CHEN Yun, et al. Simplified method for determination of thickness of composite liners based on contaminant breakthrough time[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(10): 1841-1848. (in Chinese) doi: 10.11779/CJGE202010009
[22] XIE H J, JIANG Y S, ZHANG C H, et al. An analytical model for volatile organic compound transport through a composite liner consisting of a geomembrane, a GCL, and a soil liner[J]. Environmental Science and Pollution Research, 2015, 22(4): 2824-2836. doi: 10.1007/s11356-014-3565-5
[23] HAO Z S, BARLAZ M A, DUCOSTE J J. Finite-element modeling of landfills to estimate heat generation, transport, and accumulation[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2020, 146(12): 04020134.
[24] PARK S, YOON S, KWON S, et al. Temperature effect on the thermal and hydraulic conductivity of Korean bentonite buffer material[J]. Progress in Nuclear Energy, 2021, 137: 103759.
[25] FOX P M, TINNACHER R M, CHESHIRE M C, et al. Effects of bentonite heating on U(VI)adsorption[J]. Applied Geochemistry, 2019, 109: 104392.
[26] ZHAN T L T, GUAN C, XIE H J, et al. Vertical migration of leachate pollutants in clayey soils beneath an uncontrolled landfill at Huainan, China: a field and theoretical investigation[J]. Science of the Total Environment, 2014, 470/471: 290-298.
[27] 何晓晓. 江村沟垃圾渗滤液尾水的深度处理研究[D]. 西安: 长安大学, 2013. HE Xiaoxiao. Deep Research on the Tail Water of Leachate from Jiang Cungou MSW Landfill[D]. Xi'an: Changan University, 2013. (in Chinese)
[28] 陈云敏, 刘晓成, 徐文杰, 等. 填埋生活垃圾稳定化特征与可开采性分析: 以中国第一代卫生填埋场为例[J]. 中国科学: 技术科学, 2019, 49(2): 199-211. https://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201902010.htm CHEN Yunmin, LIU Xiaocheng, XU Wenjie, et al. Analysis on stabilization characteristics and exploitability of landfilled municipal solid waste: case of a typical landfill in China[J]. Scientia Sinica(Technologica), 2019, 49(2): 199-211. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201902010.htm
[29] 张春华. 填埋场复合衬垫污染物热扩散运移规律及其优化设计方法[D]. 杭州: 浙江大学, 2018. ZHANG Chunhua. Mechanisms for Contaminant Transport in Landfill Composite Liners under Thermal Effect and Its Optimization Design Method[D]. Hangzhou: Zhejiang University, 2018. (in Chinese)
[30] ROWE R. Short-and long-term leakage through composite liners: The 7th Arthur Casagrande Lecture 1[J]. Canadian Geotechnical Journal, 2012, 49(2): 141-169.
-
期刊类型引用(7)
1. 刘志，杨阳. 黄原胶联合微生物加固改善混合砂蓄水特性试验研究. 土木与环境工程学报(中英文). 2025(03): 49-57 . 
百度学术
2. 梁越，冉裕星，许彬，张鑫强，何慧汝. 细颗粒含量影响渗流侵蚀规律的细观机理研究. 岩土工程学报. 2025(05): 1099-1106 . 本站查看
3. 高霞，梅雯艳，吴强，崔祥龙. 球度和扁平度对含瓦斯水合物煤体宏细观力学性质影响研究. 中国安全生产科学技术. 2025(05): 86-94 . 百度学术
4. 蔡国庆，刁显锋，杨芮，王北辰，高帅，刘韬. 基于CFD-DEM的流-固耦合数值建模方法研究进展. 哈尔滨工业大学学报. 2024(01): 17-32 . 百度学术
5. 程建毅，郭晓军，李泳. 泥石流物源土体标度分布参数与粘聚力的关系. 山地学报. 2024(03): 401-410 . 百度学术
6. 孙增春，刘汉龙，肖杨. 砂-粉混合料的分数阶塑性本构模型. 岩土工程学报. 2024(08): 1596-1604 . 本站查看
7. 孙建强. 颗粒形状对黄土干密度与抗剪强度的影响. 黑龙江工程学院学报. 2024(06): 8-15+23 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 287
- HTML全文浏览量: 61
- PDF下载量: 55
- 被引次数: 12
