非饱和土壤轻油污染多相流研究进展

轻油泄漏到地下后，在非饱和土壤中水、气、油三相相互作用、相互影响，共同构成多相流体系统。在简要分析轻油污染问题的产生、危害以及轻油泄漏在地表以下运移方式的基础上，介绍并总结了近些年来国内外土壤轻油污染有关的多相流研究方面的进展，包括国外的模型研究成果、国内外试验研究方法和研究成果，并指出目前研究中存在的问题和需要进一步探讨的方向。     

表面活性剂冲洗法治理非水相流体污染多相流研究进展

表面活性剂冲洗法是一种治理土壤与地下水非水相流体(NAPLs)污染的有效技术.在简要介绍表面活性剂冲洗法去除作用机理的基础上,总结了近些年来与表面活性剂冲洗法多相流相关的国内外试验研究和模型研究成果,主要包括表面活性剂作用所引起的混溶驱替和不混溶驱替多相流问题,并指出目前研究中存在的问题和进一步探讨的方向.     

表面活性剂冲洗法治理非水卡流体污染多相流研究进展

表面活性剂冲洗法是一种治理土壤与地下水非水相流体（NAPLs）污染的有效技术。在简要介绍表面活性剂冲洗法去除作用机理的基础上，总结了近些年来与表面活性剂冲洗法多相流相关的国内外试验研究和模型研究成果，主要包括表面活性剂作用所引起的混溶驱替和不混溶驱替多相流问题，并指出目前研究中存在的问题和进一步探讨的方向。     

表面活性剂冲洗修复多氯联苯污染土壤多相流研究

多氯联苯(PCBs)是一种具有持久性、抗生物降解性、脂肪溶性和明显的生物毒性等特性的持久性有机污染物,PCBs在土壤中难于准确定位、难被分解和强烈吸附,去除土壤中PCBs比较困难.表面活性剂冲洗法可以通过提高PCBs溶解度和降低水-PCBs界面张力来实现PCBs从土壤中去除;表面活性剂冲洗PCBs污染土壤涉及气相、水相、NAPLs相和固相等物质,是多相共存并相互发生作用的过程,发生相对渗透率、饱和度和毛细压力的变化;另外,为研究表面活性剂冲洗土壤中PCBs的去除机理,并降低PCBs对研究人员的危害,采用微观孔隙结构网络模型是一种较新颖的和效果显著的研究方法.开展表面活性剂冲洗PCBs污染土壤多相流研究,可以为PCBs污染场地修复提供理论基础和实验支持,并促进我国POPs履约工作的顺利进行.     

污染土壤毒性研究方法进展

综述了国内外有关污染土壤毒性研究方法 ,包括传统的方法如植物法、动物法和微生物法等和生态毒理学方法 ,并从生物、非生物和环境等方面论述了影响土壤毒性的因素 ,提出了当前土壤污染毒性研究方法中存在的问题 ,认为随着环境科学技术的发展 ,污染土壤毒性研究方法在环境保护中必将发挥越来越重要的作用     

重金属污染土壤的强化电动修复技术研究进展

在土壤重金属污染日益严重的背景下,寻找一种简单高效的土壤修复技术刻不容缓。电动修复技术由于操作简单、处理污染物多样、修复快速、成本低廉等优点,成为重金属污染土壤修复领域的研究热点。简述了中国土壤重金属污染的现状,电动修复的原理及其影响因素,综述了近年来国内外有关强化电动修复重金属污染土壤的研究进展及国内外应用实例,分析了电动修复技术中存在的问题,最后指明了强化电动修复技术未来的研究方向。     

土壤中邻苯二甲酸脂提取与检测方法的研究进展

邻苯二甲酸酯(PAEs)是一类环境激素类化合物,具有致畸、致癌和致突变的毒性。土壤中的PAEs可以通过食物链在人体内富集,威胁人类身体健康。当前土壤受PAEs污染现象日益普遍,选择合适的提取和检测方法评估土壤的PAEs污染程度显得至关重要。在前人研究成果的基础上,比较了索氏提取法、振荡提取法等传统方法与加速溶剂提取法等新型方法共7种土壤样品提取方法的优缺点,并综述了国内外检测土壤样品中PAEs含量的主要方法。展望了土壤中PAEs提取方法的发展趋势,为今后高效、快速检测土壤中的PAEs污染状况提供科学依据。     

污染土壤毒性研究方法进展

综述了国内外有关污染土壤毒性研究方法，包括传统的方法如植物法、动物法和微生物法等和生态毒理学方法，并从生物、非生物和环境等方面论述了影响土壤毒性的因素，提出了当前土壤污染毒性研究方法中存在的问题，认为随着环境科学技术的发展，污染土壤毒性研究方法在环境保护中必将发挥越来越重要的作用。     

土壤柴油污染修复的抽气提取去除实验研究

为得到土壤气相抽提(SVE)去除柴油的优化条件,进行了一维土柱抽气提取去除柴油污染物的实验研究,研究不同初始含水率、不同抽气量对污染土壤中柴油去除率的影响及不同深度残留柴油的变化规律。结果表明：在本实验模拟的范围内,抽气量越大,SVE处理效果越好;初始含水率越低,处理效果越好;此外,不同深度去除率变化的规律基本上是随深度的增大而减小。实验结果可为土壤轻油污染实际治理提供实验数据基础。     

中国土壤石油污染的危害及治理对策

近年来中国陆上石油泄漏事故频发,严重污染土壤,对生态环境和经济建设等造成很大的负面影响.在分析土壤石油污染危害的基础上,着重探讨了土壤石油污染治理的技术措施,提出了综合防治的政策性建议.     

Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zone

A quantitative two-dimensional laboratory experiment was conducted to investigate the immiscible flow of a light non-aqueous phase liquid (LNAPL) in the vadose zone. An image analysis technique was used to determine the two-dimensional saturation distribution of LNAPL, water and air during LNAPL infiltration and redistribution. Vertical water saturation variations were also continuously monitored with miniature resistivity probes. LNAPL and water pressures were measured using hydrophobic and hydrophilic tensiometers. This study is limited to homogeneous geological conditions, but the unique experimental methods developed will be used to examine more complex systems. The pressure measurements and the quantification of the saturation distribution of all the fluids in the entire flow domain under transient conditions provide quantitative data essential for testing the predictive capability of numerical models. The data are used to examine the adequacy of the constitutive pressure-saturation relations that are used in multiphase flow models. The results indicate that refinement of these commonly used hydraulic relations is needed for accurate model prediction. It is noted in particular that, in three-fluid phase systems, models should account for the existence of a residual NAPL saturation occurring after NAPL drainage. This is of notable importance because residual NAPL can act as a non negligible persistent source of contamination.     

Semi-analytical solution of two-dimensional sharp interface LNAPL transport models

In this paper, we present semi-analytical solutions for two-dimensional equations governing transport of Light Non-Aqueous Phase Liquids (LNAPL) in unconfined aquifers. The proposed model is based on sharp interface displacement and steady groundwater flow assumptions, where both the water–LNAPL interface and the LNAPL–air interface are represented as sharp interfaces. In the case of steady groundwater flow, these equations can be reduced to a two-dimensional nonlinear solute transport equation, with the LNAPL thickness in the free product lens being the primary unknown variable. The linearized form of this solute transport equation falls into the category of two-dimensional transport equation with time-dependent dispersion coefficients. This equation can be solved analytically for an infinite domain region. In this paper, the general form of the analytical solution for the transport equation, as well as the solutions for some specific cases are presented. To demonstrate the utility of the proposed solution, numerical results obtained for two example problems are discussed and presented comparatively with a finite-element solution and other more restrictive solutions available in the literature. Although the solutions discussed in this paper have some simplifying assumptions, such as sharp-interfaces between fluid phases, steady groundwater flow and homogeneous aquifer properties, the semi-analytical solutions presented in this study may be used effectively as bench mark solutions in evaluating LNAPL migration in the subsurface. These solutions are simple and cost effective to implement and may be used in the calibration of other more complex numerical solutions that can be found in the literature.     

Modeling dissolution and volatilization of LNAPL sources migrating on the groundwater table

A vertically averaged two-dimensional model was developed to describe areal spreading and migration of light nonaqueous-phase liquids (LNAPLs) introduced into the subsurface by spills or leaks from underground storage tanks. The NAPL transport model was coupled with two-dimensional contaminant transport models to predict contamination of soil gas and groundwater resulting from a LNAPL migrating on the water table. Numerical solutions were obtained by using the finite-difference method. Simulations and sensitivity analyses were conducted with a LNAPL of pure benzene to study LNAPL migration and groundwater contamination. The model was applied to subsurface contamination by jet fuel. Results indicated that LNAPL migration were affected mostly by volatilization. The generation and movement of the dissolved plume was affected by the geology of the site and the free-product plume. Most of the spilled mass remained as a free LNAPL phase 20 years after the spill. The migration of LNAPL for such a long period resulted in the contamination of both groundwater and a large volume of soil.     

A pragmatic approach for estimation of source-zone emissions at LNAPL contaminated sites

When considering natural attenuation as a remediation strategy at a site contaminated by a light non-aqueous phase liquid (LNAPL), it is important to consider the emission of contaminants from the source zone. A quantification of source-zone emissions is essential both for comparison with down-gradient mass fluxes to provide an estimate of fractional mass flux reduction, as well as for estimating the source lifetime. Because the spatial distribution of LNAPL at a field site is strongly dependent on both the spill circumstances and the heterogeneity of the geologic materials, which can be problematic for in-situ determination, alternative methods for estimating source-zone emissions are needed. In this work, a three-dimensional multiphase flow and transport modelling approach is used to investigate the relationship between the lateral extent of an LNAPL body and the emission of contaminants to groundwater at a contaminated site. For simulations involving an LNAPL release in an aquifer comprised of heterogeneous porosity and permeability distributions that were generated geostatistically, it is shown that a simple linear relationship exists between the lateral extent of the LNAPL body in the capillary fringe and the emission to the aqueous phase. The parameters describing the relationship are found to be linear functions of the groundwater flow velocity and the vertical infiltration rate. This site-specific relationship provides a simple method to estimate contaminant emissions to groundwater at LNAPL contaminated sites.     

Multispectral image analysis method to determine dynamic fluid saturation distribution in two-dimensional three-fluid phase flow laboratory experiments

The need for measuring dynamic fluid saturation distribution in multi-dimensional three-fluid phase flow experiments is hampered by lack of appropriate techniques to monitor full field transient flow phenomena. There is no conventional technique able to measure dynamic three-fluid phase saturation at several array points of the flow field at the same time. A multispectral image analysis technique was developed to determine dynamic NAPL, water and air saturation distribution in two-dimensional three-fluid phase laboratory experiments. Using a digital near-infrared camera, images of sand samples with various degrees of NAPL, water and air saturation were taken, under constant lighting conditions and within three narrow spectral bands of the visible and near-infrared spectrum. It was shown that the optical density defined for the reflected luminous intensity was a linear function of the NAPL and the water saturation for each spectral band and for any two and three-fluid phase systems. This allowed the definition of dimensionless lump reflection coefficients for the NAPL and the water phase within each spectral band. Consequently, at any given time, two images taken within two different spectral bands provided two linear equations which could be solved for the water and the NAPL saturation. The method was applied to two-dimensional three-phase flow experiments, which were conducted to investigate the migration and the distribution of LNAPL in the vadose zone. The method was used to obtain continuous, quantitative and dynamic full field mapping of the NAPL saturation as well as the variation of the water and the air saturation during NAPL flow. The method provides a non-destructive and non-intrusive tool for studying multiphase flow for which rapid changes in fluid saturation in the entire flow domain is difficult to measure using conventional techniques.     

On the selection of primary variables in numerical formulation for modeling multiphase flow in porous media

Selecting the proper primary variables is a critical step in efficiently modeling the highly nonlinear problem of multiphase subsurface flow in a heterogeneous porous-fractured media. Current simulation and ground modeling techniques consist of (1) spatial discretization of mass and/or heat conservation equations using finite difference or finite element methods; (2) fully implicit time discretization; (3) solving the nonlinear, discrete algebraic equations using a Newton iterative scheme. Previous modeling efforts indicate that the choice of primary variables for a Newton iteration not only impacts computational performance of a numerical code, but may also determine the feasibility of a numerical modeling study in many field applications. This paper presents an analysis and general recommendations for selecting primary variables in simulating multiphase, subsurface flow for one-active phase (Richards' equation), two-phase (gas and liquid) and three-phase (gas, water and nonaqueous phase liquid or NAPL) conditions. In many cases, a dynamic variable switching or variable substitution scheme may have to be used in order to achieve optimal numerical performance and robustness. The selection of primary variables depends in general on the sensitivity of the system of equations to the variables selected at given phase and flow conditions. We will present a series of numerical tests and large-scale field simulation examples, including modeling one (active)-phase, two-phase and three-phase flow problems in multi-dimensional, porous-fractured subsurface systems.     

A field trial to assess the performance of CO2-supersaturated water injection for residual volatile LNAPL recovery

A pilot scale field trial was conducted to evaluate the recovery of volatile, light non-aqueous phase liquids (LNAPLs) using a novel remediation method termed supersaturated water injection (SWI). SWI uses a patented technology to efficiently dissolve high concentrations of CO₂ into water at elevated pressures. This water is injected into the subsurface resulting in the nucleation of CO₂ bubbles at and away from the injection point. The nucleating bubbles coalesce, rise and volatilize residual LNAPL ganglia. In this study, an LNAPL composed of 103 kg of volatile pentane and hexane, and 30 kg of non-volatile Soltrol was emplaced below the water table at residual saturation. The SWI technology removed 78% of the pentane and 50% of the less volatile hexane. Contaminant mass was still being removed when the system was shut down for practical reasons. The mass removed is comparable to that expected for air sparging but a much smaller volume of gas was injected using the SWI system.     

地下水中轻质有机污染物(LNAPL)透镜体研究

在二维砂槽模型中模拟了轻质油在均匀多孔介质地下水非饱和区中的运移过程。模拟结果表明,地下水毛细区是轻质油污染的重点区,除了ＬＮＡＰＬ的残留以外,进入地下水饱和的ＬＮＡＰＬ终将被地下水顶托回到毛细区中,毛细区以上的约大多民将在重力作用下进入毛细区中,试验中观察到达稳定状态时ＬＮＡＰＬ透镜体的上边缘略微高出毛细区。利用多孔介质毛细管模型,建立了利用界面张力、接触角、介质特征孔隙直径等物理量估算不同位置     

Persistence of LNAPL sources: relationship between risk reduction and LNAPL recovery

Light nonaqueous phase liquids (LNAPLs), such as fuels, are the source of much soil and groundwater contamination. Though the mobility of LNAPLs is limited in many environments, dissolved-phase components, such as benzene, can produce groundwater plumes that are more mobile than the LNAPL source. In such a setting, it is commonly assumed that recovery of the LNAPL will result in a reduction in risk associated with the dissolved phase. This paper synthesizes several existing multiphase and chemical transport solutions into a single linked methodology that predicts concentrations of soluble constituents within and downgradient of LNAPL source zones from dissolution of those constituents into groundwater flowing through and below LNAPL sources. This approach has been applied to a variety of LNAPL spill conditions. For biodegradable compounds, these analyses show that the period of time where the dissolved-phase plume is expanding is very small compared to the duration of most LNAPL sources, and that the downgradient extent is generally less than about 100 m for BTEX compounds. Therefore, the risk to receptors, as measured by the maximum downgradient extent of dissolved-phase plume or the maximum concentration of these compounds at a downgradient receptor, is generally unrelated to the longevity of the LNAPL sources. The maximum downgradient extent of the dissolved-phase plume is determined almost entirely by the groundwater velocity and the biodegradation rate. These analyses further demonstrate that recovery of LNAPL by hydraulic methods is often ineffective at reducing risk. Except in coarse-grained soils or intermediate soils with significant LNAPL saturations, free-product recovery approaches do not result in significant reductions in the longevity of downgradient dissolved-phase contamination. Further, for biodegradable constituents, remediation does not result in a near-term decrease in the downgradient extent of contamination. Cleanup methods that act to change the composition of the LNAPL source are more effective at reducing the downgradient concentrations, particularly for fine-grained soils or when LNAPL saturations are low.