大变形黏土防渗层中的污染物迁移和转化规律研究

国内湖泊疏浚污染底泥堆场一般以较厚的黏土层作为主要防渗层,由于在上覆底泥作用下黏土层会发生较大的固结变形,因此,在研究黏土防渗层中的污染物运移和转化规律时,应该考虑土体变形的影响。基于Gibson一维大变形固结理论和饱和多孔介质中的污染物对流扩散方程,建立了二者耦合的可变形多孔介质中污染物的运移和转化模型,其中首次考虑了土体自重和生物降解作用的影响。利用所建立模型的数值解,研究了在可变形黏土防渗层中的污染物运移和转化规律,同时分析了模型中不同项和主要参数的作用和影响。研究结果表明,土体大变形对黏土防渗层中污染物的运移有着较复杂的影响,一方面土体变形会加速污染物的运移;另一方面土体固结带来的渗透性减小会增加污染物的穿透时间,二者的不同作用取决于众多的影响因素,如土层厚度和吸附作用等。研究结果对于评估天然黏土防渗层对污染物的阻隔作用有重要的指导意义。     

污染泥堆场处置中自重固结与污染物迁移的耦合分析

基于一维固结方程和对流-弥散方程,提出了一种用于污染泥堆场处置中自重固结与污染物迁移耦合问题的求解方法,其中,自重固结方程由解析解直接求解,污染物迁移方程由显示差分法求解。在自重固结方面,新方法可考虑线性的压缩关系和渗透关系;在污染物迁移方面,新方法能够考虑对流、扩散、机械弥散、线性和非线性吸附以及孔隙率相依的有效扩散系数。新方法与CST1数值模型进行了对比,验证了新方法的正确性。基于所提出的解,分析了固结效应、压缩系数、渗透系数和污染物初始分布形式对污染物在淤泥中迁移过程的影响。计算结果表明:固结效应、压缩性、渗透性和污染物的初始分布形式对污染物的迁移过程有着显著的影响;固结效应能显著地加速污染物的流出;土体的压缩系数越大其污染物累计流出质量越大;土体的渗透系数越大其污染物累计流出质量越大;污染物的初始分布形式对污染物的流出速率影响较大。     

非等温分布条件下考虑半透水边界时饱和黏土的一维固结解析解

温度的变化会导致土体的物理力学性质改变,且在一些实际工程中,饱和黏土会处于非等温分布状态。为此,针对非等温分布条件下饱和黏土的一维固结问题,考虑了更具普遍性的半透水边界,通过某些假定推导了单级线性加荷形式下饱和黏土一维固结控制方程,并利用分离变量法求解得到了控制方程的解析解。通过将所提解析解分别与已有解析解和有限差分解展开对比分析,验证了所提解答的正确性。基于所提解析解,利用某一算例分析了温度梯度、半透水边界参数及加荷时间对固结性状的影响。结果表明：温度梯度 M 越大,土体的渗透性越大,土体的固结速率越快;半透水边界参数 R1和 R2越大,相同时间内土体的超孔隙水压力越小,土体的平均固结度越大;土体的平均固结度随加荷时间 tc 的增大而减小,这主要是由于加荷阶段所施加的外荷载小于最终荷载,但加荷时间tc的延长可一定程度减小土体中产生的最大超孔隙水压力。     

固结作用下粉土的持水性能与细观机制研究

针对传统仪器无法考虑岩土工程中温度和荷载等因素综合作用的缺陷,研制了一套多功能土-水特征曲线试验仪。以粉土为研究对象,开展了不同固结应力作用下粉土的土-水特征曲线试验,并测量其增减湿过程中土体的体积胀缩量,以便修正体积变化对其体积含水率的计算误差。结果表明,粉土在减湿过程中,土体的体积发生明显的收缩。其收缩量与固结应力水平有关,固结应力越小,则减湿引起的体缩量越大。但土体在增湿过程中的体积则基本保持不变。固结应力对粉土 土-水特征曲线的进气值、增减湿速率影响较大,固结压力越大,其进气值越大、减湿速率也越大。最后,为揭示固结应力对土的持水性能的影响机制,开展了不同固结应力作用下土体的细观试验。结果表明,固结应力主要改变了土体团粒间的大孔隙,而对黏土颗粒之间的孔隙影响较小。土体的持水性能与土体的孔隙大小和分布模式密切相关,大孔隙主要影响土体的进气值,而其孔隙分布模式则控制其增减湿湿速率。     

软黏土层一维有限应变固结的超静孔压消散研究

根据土力学固结理论计算分析软黏土层固结过程的超静孔隙水压力值,确定软黏土体固结过程的强度增长,对排水固结法处理软土地基至关重要。软黏土层固结过程中土体变形较大时,有限应变固结理论和小应变固结理论计算分析软黏土固结所得结果差异较大。利用非线性有限元法及程序,通过对软黏土层固结工程算例的计算结果分析,研究了有限应变固结理论和小应变固结理论计算分析软黏土层一维固结超静孔压值消散的差异;探讨了软黏土体一维固结过程中,几何非线性、土体渗透性变化和压缩性变化对超静孔隙水压力消散的影响。研究结果表明,当土体的变形较大时,有限应变固结理论计算出的超静孔压要比小应变固结理论得到的值消散的更快。考虑土体固结过程中渗透性的变化时,超静孔压消散变慢;可用软黏土渗透性变化指数ck 反映渗透性变化对超静孔压消散的影响,渗透性变化指数ck值越小、超静孔压消散越慢。固结过程中软黏土压缩性的大小及变化也影响超静孔压的消散,可用软黏土的压缩指数cc反映固结过程中压缩性的大小及变化对超静孔压消散的影响,软黏土的压缩指数cc越小,固结过程软黏土层中的超静孔压消散越快。     

非饱和土化学-水力-力学耦合行为数值模拟

在描述非饱和土力学本构行为的Gens-Alonso模型[1]的基础上,结合化学浓度对前固结压力的软化公式[2],并考虑了化学浓度对黏聚力的软化作用,建立了描述非饱和土化学-水力-力学耦合行为的本构模型。对包括化学-水力-力学本构关系、流体流动和污染物传输的控制方程的平衡方程组进行有限元离散,并采用此非饱和土化学-水力-力学耦合行为的数值模型对非饱和土中隧道周围化学-水力-力学耦合过程进行了模拟,对纯粹的力学响应、弹性化学力学响应、弹塑性化学力学响应 3 种情况进行了分析,数值结果显示了化学浓度作用下隧道周围的体积变化量和应力状态比单纯的力学响应有较大的提高,隧道表面的收缩量也有明显的增大。数值模拟的结果显示了该方法对模拟非饱和土复杂的耦合行为的可行性。     

地下水开采引发的土体变形对污染物迁移影响研究

地下水开采引发的土体变形对污染物迁移具有重要影响。结合地下水渗流理论、太沙基一维固结理论和溶质运移对流-弥散理论,建立了变形土体中污染物迁移三维耦合数值模型,考虑了孔隙度、渗透系数和水动力弥散系数随有效应力的动态变化。通过对实际案例进行模拟,分析了地下水开采引发的土体变形对污染物迁移规律的影响。结果表明,地下水位下降,土体有效应力增加,土体骨架压缩变形,使得孔隙度、渗透系数减小,导致水动力弥散系数增幅减小,从而延缓了污染物迁移的对流扩散过程。     

高饱和条件下水-力-化耦合模型及数值模拟

填埋场底部黏土垫层特性对其长期防渗隔污服役性能的有效发挥具有极其重要的作用。针对实际堆场中黏土屏障呈现出的高饱和状态,将液相(孔隙水)与气相(闭塞气泡)视为混合流体,通过分别建立土体应力平衡方程、混合流体质量守恒方程及溶质质量守恒方程,综合考虑土颗粒-孔隙流体-溶质间的相互作用机制,推导得到了高饱和度条件下,溶质在黏土防渗层中运移的水-力-化全耦合模型,能够实现多物理场耦合作用时土层变形量、混合流体压力及溶质溶度随时空分布的直接精确求解。采用多场耦合有限元分析软件COMSOL对所建模型开展数值模拟,结果表明,模型结果与Peters所得结果吻合较好;黏土垫层中可压缩性气体的存在延缓了超孔隙流体压力的消散,加大了土层的沉降量,对溶质在防渗黏土层中的运移进程起到了显著的阻滞效果。     

循环荷载作用下超固结软黏土软化-孔压模型研究

随着循环次数的增加,循环荷载作用下正常固结软黏土孔隙水压力逐渐增加,而孔压的增加将导致土体发生软化现象。对于超固结软黏土,循环荷载作用下在循环初期将产生负孔压,而此时也发生循环软化现象,这显然与有效应力原理是相矛盾的。以往的研究往往通过建立软化指数与循环次数的关系来描述土体的循环软化特性,从而不能反映土体的循环过程中残余孔压变化对土体循环软化的影响。通过对杭州饱和软黏土进行应力控制的循环三轴试验,对循环荷载作用下超固结软黏土的软化特性及孔压发展规律及两者关系进行了研究。在试验的基础上建立了超固结软黏土循环软化-孔压模型,该模型反映了残余孔压增长对超固结软黏土循环软化特性的影响规律。     

湿度影响下莫高窟壁画地仗层吸附水及吸力变化特征研究

莫高窟气候干燥,窟内壁画地仗层常处于低含水率、高吸力的状态,受到洞窟环境湿度波动的影响较大。为了进一步了解湿度影响下地仗层土体中吸附水以及吸力变化特征,本文通过蒸汽吸附法测试莫高窟不同地仗层高吸力段内的土水特征曲线,并且根据测试结果分析土体中由于范德华力与毛细凝聚作用产生的不同吸力的变化特征;同时进一步利用蒸汽吸附法测试了含NaCl地仗层土水特征曲线并探讨水汽吸附过程中渗透吸力特征;此外利用热重分析方法测试地仗层土体中在范德华力与毛细凝聚作用下所吸附水分的特征。结果表明:地仗层水汽吸附过程中土体内大部分基质吸力来自于毛细凝聚作用,土颗粒分子间范德华力仅在含水率极低的情况下为土体提供比较明显的吸力,湿度影响下地仗层水汽吸附过程中土体内渗透吸力仅在外界环境湿度大于地仗层中盐分的潮解临界湿度时才比较明显;地仗层土体所吸附的水分大部分以弱结合水形式存在,并且地仗层中澄板土含量越高,弱结合水临界含水率越大。结论可以为湿度影响下的壁画病害机理以及地仗层中水汽运移等问题的研究提供理论基础。     

Parametric sensitivity analysis of coupled mechanical consolidation and contaminant transport through clay barriers

In this paper, an extensive parametric sensitivity analysis of coupled consolidation and solute transport in composite landfill liner systems has been undertaken. The analysis incorporates results of more than 3000 simulations for various combinations of barrier thickness, waste loading rate, initial void ratio, compression index, hydraulic conductivity and dispersion coefficient. However, it is noted that to limit the extent of the study a constant coefficient of consolidation is assumed in the analysis presented here, though this assumption is easily relaxed. Results of the parametric sensitivity analysis are succinctly presented using dimensionless plots, which allow the comparison of results for a large number of parameter values, and so the clear identification of the most important determinants on contaminant transport through the liner system. The dimensionless plots demonstrate a pessimum (for which the ‘breakthrough time’ is minimised). Numerical results reveal that in cases of extreme liner compressibility an order of magnitude reduction in contaminant transit time may arise due to coupling between solute transport and consolidation, while for barriers of low compressibility and porosity (such as well-engineered composite compacted clay landfill liners), it is found that the contaminant transit time may still be reduced by more than 30%. The numerical results suggest that the use of coupled consolidation–contaminant transport models are sometimes required for informed and conservative landfill liner design.     

Theoretical investigation of the effects of consolidation on contaminant transport through clay barriers

Consolidation of clayey contaminant barriers such as landfill liners has been postulated as a cause of early breakthrough of contaminants. In this paper we theoretically investigate this proposition. For this purpose a sophisticated one‐dimensional, large‐deformation model of coupled mechanical consolidation and solute transport is employed. This new model is a generalization of existing coupled consolidation and solute transport models described in the literature. It takes into account both non‐linearities in geometry as well as constitutive relations. The latter relate the compressibility, hydraulic conductivity and coefficient of effective diffusivity to the deformation of the soil. The model is applied to a case study of a clay liner and geomembrane system. Results obtained from numerical solution of the model equations are compared with those from various simplified models, including a ‘diffusion only’ (i.e. a rigid soil) model traditionally used in contaminant barrier design. For barriers incorporating low compressibility soils (as for many well compacted clays), there is little difference between contaminant transit (i.e. breakthrough) times predicted by the two models. However, for contaminant barriers incorporating more compressible soils, consolidation is shown to significantly accelerate transport. These results indicate the potential importance of accounting for the effects of soil consolidation and highlight the limitations of existing models when modelling solute transport through composite barriers utilizing soft soils. Based on these limited results, we suggest a possible way of taking into account soil consolidation using simplified models. Copyright © 2008 John Wiley & Sons, Ltd.     

One‐dimensional contaminant transport through a deforming porous medium: theory and a solution for a quasi‐steady‐state problem

Contaminant migration through soil is usually modelled mathematically using the dispersion–advection equation. This type of model finds application when planning the remediation of contaminated land, predicting the movement of polluted groundwater and designing engineered landfills. Usually the analysis assumes that the porous media through which the contaminant migrates is stationary. However, the construction of landfills on clay soils means that the soil beneath the landfill will undergo time‐dependent deformation as the soil consolidates. To date, there are no published data on the effect a deforming porous media may have on contaminant transport beneath a landfill; indeed, there appears to be no theory of contaminant migration through a deforming soil. In this paper, a one‐dimensional theory of contaminant migration through a saturated deforming porous media is developed based on a small and large strain analysis of a consolidating soil and conservation of contaminant mass. By selection of suitable parameters, the new transport equation reduces to the familiar one‐dimensional dispersion–advection equation for a saturated soil with linear, reversible, equilibrium controlled sorption of the contaminant onto the soil skeleton. Analytic solutions to a quasi‐steady‐state contaminant transport problem for a deforming media are presented, and a preliminary assessment made of the potential importance of soil deformation on the results of a contaminant migration analysis. Copyright © 2000 John Wiley & Sons, Ltd.     

复合衬层中变系数有机污染物迁移规律分析

土工膜和下伏压实黏土组成的复合衬层已被广泛用作填埋场的防渗屏障系统。有机污染物在复合衬层中迁移时,其主要运移机制是扩散作用。假设有机污染物在土工膜中稳态扩散,并假设在下伏衬层中扩散系数为迁移距离的线性函数,且考虑降解作用的影响,建立了有机污染物在复合衬层中的一维扩散模型,针对零浓度下边界条件,获得了模型解析解。基于该解析解,分析讨论了相关参数的敏感性。结果表明,相关参数对计算结果影响很大,降解半衰期对污染物运移同样存在较大影响。该解析模型仅适用于有机污染物在两层复合衬垫中的迁移问题,可为填埋场的初步设计提供参考     

Efficacy of Fine-Grained Soil as Landfill Liner Material for Containment of Chrome Tannery Sludge

This paper evaluates the potential use of a fine-grained soil obtained from a site in West Bengal, India, as a suitable landfill liner material for the containment of hexavalent chromium from tanning waste sludge. The physico-chemical properties of the soil were determined. The soil adsorption affinity for hexavalent chromium was also assessed through adsorption batch and breakthrough column tests. The zero point charge (pH_zpc) of the soil was found to be 7.3. The batch kinetics and column tests results indicated that the soil liner possesses a relatively good hexavalent chromium adsorption capacity at natural or slightly alkaline condition. The adsorption results showed that the hexavalent chromium uptake by the soil follows both Langmuir and Freundlich adsorption isotherms. This study also illustrated that the hexavalent chromium breakthrough curve in the column experiment reached equilibrium concentration after 3.5 pore volumes (900 h). Overall, this study showed that the fine-grained soil has the potential for usage as a landfill liner or as a component of a landfill barrier system to prevent chromium contamination from the tannery waste disposal.     

Solute transfer during consolidation based on a solid‐fluid‐solute coupling model

The migration of contaminant through soil is usually modeled using the advection‐dispersion equation and assumes that the porous media is stationary without introducing a constitutive equation to represent soil structure. Consequently, time‐dependent deformation induced by soil consolidation or physical remediation is not considered, despite the need to consider these variables during planning for the remediation of contaminated ground, the prediction of contaminated groundwater movement, and the design of engineered landfills. This study focuses on the numerical modeling of solute transfer during consolidation as a first step to resolve some of these issues. We combine a coupling theory‐based mass conservation law for soil‐fluid‐solute phases with finite element modeling to simulate solute transfer during deformation and groundwater convection. We also assessed the sensitivity of solute transfer to the initial boundary conditions. The modeling shows the migration of solute toward the ground surface as a result of ground settlement and the dissipation of excess pore water pressure. The form of solute transport is dependent on the ground conditions, including factors such as the loading schedule, contamination depth, and water content. The results indicate that an understanding of the interaction between coupling phases is essential in predicting solute transfer in ground deformation and could provide an appropriate approach to ground management for soil remediation.     

渗透系数的变异性对压实黏土衬垫性能的影响

压实黏土衬垫是填埋场等环境岩土工程中常用的屏障材料,其厚度和渗透系数是主要的设计参数。考虑渗透系数的变异性,用衬垫底部污染物相对浓度和通量等指标分析了设计参数对衬垫有效性的影响。研究表明,渗透系数的变异性对衬垫性能有很大的影响：当衬垫厚度较小或渗透系数均值较大时,衬垫底部出现高浓度（接近于1）的概率很大,这个浓度区间可能比相对浓度的均值高很多,对周围环境污染有较大的威胁。从非反应性溶质的运移来看,我国规范中对衬垫厚度和渗透系数的规定是合理的。相对于污染物的通量,厚度和渗透系数对衬垫底部相对浓度的影响更为显著,相对浓度更适合作为评价衬垫性能的指标。     

The 14-year Performance of a Compacted Clay Liner used as Part of a Composite Liner System for a Leachate Lagoon

The migration of contaminants through a 2.9 m thick compacted clay liner (CCL) for a landfill leachate lagoon is examined 14 years after construction. The clay liner formed the lower portion of the composite liner system but the geomembrane (GM) was found to have defects that had allowed leachate to migrate between the GM and CCL. Chloride, sodium, potassium, calcium and magnesium pore water profiles through the CCL are examined. It is shown that chloride migrated approximately 1.7 m into the CCL during the 14 years of the lagoon operation, sodium approximately 1.2 m, and potassium 0.7 m. Diffusion and sorption data from laboratory diffusion testing are utilized in combination with a finite layer contaminant transport model to predict field contaminant migration profiles through the composite liner system and to establish the time of ‘failure’ of the geomembrane at sometime between 0 and 6 years after installation. Relatively high sorptive uptake of potassium by the CCL soil is observed from batch testing and diffusion testing with field data suggesting an even larger amount of sorption. It is hypothesized that organic sludge matter at the base of the lagoon is responsible for potassium uptake from the leachate. This field case highlights the importance of the compacted clay liner as part of the composite liner system in acting as a diffusion barrier during the lifetime of the lagoon as well as using relatively non-conservative contaminants such as chloride and sodium to estimate geomembrane ‘failure’ times     

CCL吸附特性及孔隙率降低对污染物运移的影响

假定孔隙均匀地分布于土体的物质空间内和土骨架对污染物的吸附特性服从平衡线性,对基本体积质量关系进行分析,提出了由于土体对污染物的吸附而引起的孔隙率降低的估算公式。在考虑土体孔隙率变化的条件下,建立了污染物一维运移的控制方程,并考虑垃圾生物降解效应、压实黏土衬里（CCL）防渗层、下覆有限厚度含水层等实际情况,确定了初始条件和边界条件。对所建立的初边值问题进行了数值求解,且对某假想填埋场情况进行了变动参数与对比计算,结果表明,由于土颗粒对污染物的吸附所引起的孔隙率降低,显著地降低了污染物对压实黏土衬里的穿透能力。与常孔隙率情况相比,CCL中污染物的峰值浓度降低近10 %,含水层中污染物浓度降低更显著。当考虑土体孔隙率变化时,弥散对污染物运移具有控制作用,分布系数对污染物的运移具有重要影响。