非平衡-非线性吸附情况下填埋场污染物运移分析

采用Langmuir等温吸附线方程描述非线性吸附性能, 基于改进的混合元方法，通过数值计算与分析，探讨了非平衡吸附条件下污染物运移过程及其机理。计算结果表明: 当考虑非平衡、非线性吸附性能时, 污染物穿透曲线即浓度的时程变化曲线尖锐而狭窄、峰值点前移, “拖长尾”现象不明显, 由此说明，土颗粒对污染物的非平衡、非线性吸附使得污染物的穿透能力增强, 滞留能力下降。进一步的变动参数比较分析表明：Langmuir等温线方程中的参数B、压实粘土衬里的渗透性及地下水渗流速度对污染物运移过程具有显著的影响。     

考虑非等温环境下污染物在黏土中的运移解析模型

利用广义积分变换法推导了非等温条件下污染物在压实黏土衬垫中的运移解析解。所提出的解析解考虑了分子扩散、对流、吸附和热扩散的耦合效应,并考虑了压实黏土的渗透系数、分布系数和有效扩散系数3个参数随温度的变化。通过与热扩散试验结果、已有的解析解和基于COMSOL软件的数值模型的对比,验证了所提出的解析解。利用经过验证的解析解,研究了非等温环境以及压实黏土的渗透系数、有效扩散系数和分布系数3个参数随温度的变化对甲苯在压实黏土衬垫系统中迁移的影响。结果表明,非等温环境以及压实黏土的渗透系数、有效扩散系数和分布系数随温度的变化均对甲苯在压实黏土衬垫中迁移有着显著的影响。不考虑非等温环境的影响将极大低估污染物的溢出量和污染物的击穿时间。既有解析解忽略压实黏土的渗透系数和有效扩散系数随温度的变化会极大低估甲苯的流出速率,而忽略压实黏土分布系数随温度的变化会极大高估甲苯击穿衬垫系统的时间和达到稳态的时间。所提出的解析解能够考虑热扩散作用,同时考虑了压实黏土的渗透系数、分布系数和有效扩散系数3个参数随温度的变化,较既有解析解更贴近工程实际,能够为压实黏土衬垫系统的设计和服役性能评价提供指导和借鉴。     

非等温条件下有机污染物在粘土衬垫中的扩散分析

填埋场中垃圾降解作用产生的热量在黏土衬垫中的运移会形成温度场,温度场会对污染物的运移产生影响。文章在多孔介质热力学理论和污染物运移理论相结合的基础上,主要考虑温度场对污染物扩散的影响,建立了非等温条件下有机污染物在黏土衬垫中的一维扩散模型,并采用分离变量法得到了模型的解析解。基于文中模型,对比了非等温条件和等温条件下有机污染物在黏土衬垫中的运移。非等温条件的黏土衬垫底部浓度比等温条件下的浓度大,非等温模型比等温模型更偏于安全。解析解可为实验数据的拟合及分析、衬垫的设计提供参考。进一步的参数分析表明:增大阻滞因子Rd和黏土厚度L可以减小黏土衬垫底部由温度梯度引起的有机污染物浓度。     

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

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

非等温条件下有机污染物在黏土衬垫中的扩散分析

填埋场中垃圾降解作用产生的热量在黏土衬垫中的运移会形成温度场,温度场会对污染物的运移产生影响。文章在多孔介质热力学理论和污染物运移理论相结合的基础上,主要考虑温度场对污染物扩散的影响,建立了非等温条件下有机污染物在黏土衬垫中的一维扩散模型,并采用分离变量法得到了模型的解析解。基于文中模型,对比了非等温条件和等温条件下有机污染物在黏土衬垫中的运移。非等温条件的黏土衬垫底部浓度比等温条件下的浓度大,非等温模型比等温模型更偏于安全。解析解可为实验数据的拟合及分析、衬垫的设计提供参考。进一步的参数分析表明:增大阻滞因子Rd和黏土厚度L可以减小黏土衬垫底部由温度梯度引起的有机污染物浓度。     

非等温分布条件下三层复合衬垫中有机污染物一维瞬态运移规律研究

针对非等温分布条件下由土工膜(geomembrane,简称GMB)+土工复合膨润土衬垫(geosynthetic clay liner,简称GCL)+压实黏土衬垫(compacted clay liner,简称CCL)组成的3层复合衬垫中有机污染物一维瞬态运移问题,考虑了对流、扩散、机械弥散、吸附、降解和热扩散等因素,建立了相应的数学模型。该模型可考虑扩散系数和渗透系数等参数随温度变化。利用有限差分法,获得了该模型的数值解。通过将所建模型计算结果分别与试验结果、已有解析模型计算结果和COMSOL软件计算结果进行对比,验证了所建模型的正确性。基于所定义的击穿时间t_b,以甲苯作为代表性有机污染物分析和讨论了不同因素对运移行为的影响。结果表明：(1)非等温分布条件会使渗透系数、扩散系数和线性吸附系数等参数发生变化,其中渗透系数和扩散系数的变化会加快运移过程,尤其是扩散系数,而线性吸附系数的变化会减慢运移速率;(2)复合衬垫上部和底部温度差ΔT的增大会使击穿时间t_b减小,底部运移通量J_b增大;当ΔT为0、10、20、30、40 ...     

污染场地六价铬迁移转化机制与数值模拟研究

随着全球工业化迅猛发展,土壤和地下水六价铬污染日益严重。基于某铁合金厂铬渣场地现场调查与采样分析,开展铬渣场地土样吸附、渗透和弥散试验,研究六价铬在粉质黏土土样中的吸附特性和迁移规律,建立考虑对流-弥散-吸附的六价铬迁移三维动力学模型,结合数值软件获取污染源位于场地上、下游时地下水中六价铬迁移分布特征,并揭示弥散度?和分配系数 对六价铬时空分布的影响。试验结果表明,粉质黏土对六价铬吸附符合Langmuir等温吸附模型,最大吸附量为466.6 mg/kg;蒸馏水和160 mg/L 六价铬溶液入渗下粉质黏土渗透系数约为6.5×10–7～6.7×10–7 cm/s,1 000 mg/L六价铬溶液的渗透系数增大至4.4×10–6 cm/s;粉质黏土水动力弥散系数D为1.4×10–4 m2/d,计算得到阻滞因子 为4.2～10。数值模拟结果表明,场地下游受到六价铬污染时,即使不考虑分子扩散作用,上游仍有被污染的风险,污染程度取决于含水层的弥散度;考虑含水层对六价铬吸附时,土体分配系数越大,六价铬污染羽分布范围越小,在预测地下水中六价铬浓度分布时应重点考虑六价铬吸附等转化过程。     

有机污染物在复合衬垫中的运移分析

有机污染物在垃圾填埋场复合衬垫中有两种不同的运移方式:土工膜上缺陷处污染物的运移和完好复合衬垫中污染物的扩散。文章分别用有限差分法和解析法对两种方式下污染物的累积量进行计算,分析了两种方式对有机污染物运移所做的贡献。研究表明:GCL(人工合成材料粘土衬垫)复合衬垫以有机污染物在完好复合衬垫中的扩散为主,可以忽略缺陷处污染物的运移;而CCL(压实黏土衬垫)复合衬垫中两种运移方式所做的贡献在同一数量级,需要考虑两种方式的共同作用。在GCL复合衬垫下铺设粘土层能有效改善GCL复合衬垫对有机污染物的阻止效果。     

基于有机污染物扩散的复合衬垫厚度的概率设计

有机污染物在复合衬垫中的运移以扩散方式为主。目前关于运移参数不确定性及其对衬垫性能影响的研究还不多,复合衬垫厚度的设计也没有考虑不确定性的影响。总结了运移参数已有的成果,采用蒙特卡罗法研究了半无限条件下运移参数的变异性对衬垫底部污染物质量浓度和通量的影响,并对衬垫厚度进行了设计。研究表明:土工膜中有机污染物的分配系数和扩散系数的变异性对质量浓度和通量的影响可以忽略,而土中有机污染物的扩散系数和阻滞系数的变异性对其影响则非常显著。确定性方法计算的结果偏于危险,有必要加强对土中有机污染物的扩散系数和阻滞系数概率特性的研究。以甲苯为例进行黏土衬垫厚度设计,结果表明:概率计算方法得到的黏土衬垫厚度约为确定性方法的1.71~1.81倍,采用概率方法进行衬垫设计更合理。     

基于混合物理论的水力-化学-力学污染物输运模型研究

由于废弃物的堆积以及填埋场上覆盖层的重量,衬垫因此承受一定的压力,发生力学固结;由于黏土本身的结构特性,当污染物通过黏土垫层时可能会发生化学固结,而固结作用可直接对污染物的对流输运产生影响,也会引起土层体积和结构的变化,从而改变其固有的输运性质。从混合物理论出发,建立了水力-化学-力学作用下的溶质输运理论框架,统一地描述了水力-化学-力学作用下变形、水的吸附或解吸附、对流和扩散现象。将输运系数考虑为有效孔隙率的函数,由此反映固结对输运参数的影响。对所建立的数学模型进行无量纲化处理,并使用有限元软件COMSOL Multiphysics进行求解,最后进行了参数分析。研究结果表明,垫层土体中吸附水和自由水之间的化学势差异越大,层间吸附水的解吸附量越大;当土体的软硬性质不同时,当层间吸附水解吸附后,土体固结程度不同,从而对污染物输运过程的影响不同。而当垫层上覆荷载增大时,土体宏观孔隙的扩张的限制作用增强,从而对土体产生一种修复作用,抑制污染物在土体中的输运。     

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     

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.     

Landfill Base Liners: Assessment of Material Equivalency and Impact to Groundwater

This paper gathered available flow and transport solutions and used them for two composite liners, consisting of geomembrane (GM) overlying either a compacted clay liner (CCL) or a geosynthetic clay liner (GCL). Its aim is to provide a guiding framework for the possible choices of (a) approaches to bottom liner design, (b) respective analytical solutions to flow and transport equations, as well as (c) parameters required for each type of solution. On the basis of the obtained results, the following recommendations are made. When the goal of analysis is to determine material equivalency, leachate flow rate is an adequate key parameter for GM-CCL composite liners. For GM-GCL composite liners, it is necessary to compute contaminant concentration or mass flux, considering (a) transport through defects for inorganic contaminants and (b) diffusion and the contribution of any available attenuation layer for organic contaminants. When the goal of analysis is to assess impact to groundwater, it is advised to calculate both discharge rate and contaminant mass flux regardless of liner type. The critical parameter for the transport calculations is the retardation factor of the contaminant, for the case of CCLs, while the results for GCLs are much less sensitive to this parameter.     

考虑热扩散的黏土衬垫中污染物运移简化计算

以单层压实黏土衬垫为研究对象,在传统对流弥散方程的基础上考虑温度梯度下污染物的热扩散运移,得到简化条件下半无限边界的解析解。研究表明,考虑热扩散效应所得的衬垫底部相对浓度和通量均大于传统对流弥散方程所得结果,S_T越大时相对浓度和通量也越大,当S_T=0.5 K-1时稳定通量达到S_T=0时稳定通量2.5倍。热扩散通量是对流通量的-AS_TD/v倍,两者均随时间的增加而增大;分子扩散通量则先增加达到峰值后减小趋于零。同时考虑温度对运移参数的影响和热扩散效应时,温度升高使衬垫底部相对浓度增大,导致击穿时间提前,污染物通量显著增大;当温度从30℃升高至70℃时,击穿时间从8.5 a提前至5.1 a,稳定通量增大约1.7倍。在进行衬垫设计时不考虑温度的影响将使污染风险大幅增加。     

Numerical simulation on the cracking and failure law of compacted clay lining in landfill closure cover system

According to the mass conservation equation of solid phase, water transport, air migration and enthalpy balance equations in a deformable and unsaturated porous medium, and the mechanics model, a thermal‐hydro‐mechanical (THM) coupling mathematical model is developed in order to accurately predict the crack propagation of the compacted clay lining (CCL) in a landfill closure cover system. Numerical simulation analysis of the cracking and failure of the CCL in a landfill closure cover system under complex conditions (landfill gas pressure, variation of soil surface temperature, moisture loss of coverage soil, and uneven settlement of landfilled MSW (municipal solid waste)) is performed. The results show that the cracks are discovered firstly on the soil surface and then extend to the deeper cover system and became bigger. The cracks had broken through the closure cover system at 30 days. Over the depth range of 0.06 m, the change rate of settlement is found to be the highest because of the strong influence of temperature. The rise of temperature around the cracks is fast when the CCL is broken through by the cracks. There is a hollow area near the cracks over the depth range from 0.15 to 0.7 m, and the temperature stabilize at around 11°C without any obvious reduction. Under the condition of evaporation, the main crack in the CCL is similar to a pumping well, which causes more rapid moisture content reduction around the cracks. Vertical deformation increases linearly with depth in the process of deformation and cracking of CCL induced by uneven settlement of landfilled MSW. Copyright © 2014 John Wiley & Sons, Ltd.     

A general approach to the modelling of contaminant transport through composite landfill liners with intact or leaking geomembranes

Finite‐element models of contaminant transport through composite landfill liners require highly refined meshes around the interface between the geomembrane and the clay layer, especially if leakage through holes in the geomembrane is considered. In addition, no general formulation for transport through leaking geomembranes can be found in the literature. The paper develops a general approach to time‐dependent contaminant migration through composite liners with intact or leaking geomembranes. Equations are derived for various combinations of system conditions including Dirichlet and Neumann boundary conditions in the waste, constant mass of contaminants in the waste, steady state or transient transport in the geomembrane, and steady state or transient seepage velocities in the mineral liner. The effect of the geomembrane on transport in the soil is converted into an equivalent boundary condition applicable at the top of the clay layer. Hence, only the media underlying the top geomembrane are explicitly represented in the numerical model, yielding a computationally efficient algorithm. The new formulation is validated in conjunction with finite‐layer, finite‐element and boundary‐element methods, by comparing its predictions to those of more conventional approaches which represent the geomembrane explicitly. The scope of the method is illustrated by modelling a landfill liner with a geomembrane leaking in five locations. Copyright © 2007 John Wiley & Sons, Ltd.