Effect of nonionic surfactant partitioning on the dissolution kinetics of residual perchloroethylene in a model porous medium

At concentrations above the critical micelle concentration, surfactants can significantly enhance the solubilization of residual nonaqueous phase liquids (NAPL) and, for this reason, are the focus of research on surfactant-enhanced aquifer remediation (SEAR). As a consequence of their amphiphilic nature, surfactants may also partition to various extents between the organic and aqueous phases, thereby affecting SEAR performance. We report here on the observation and analysis of the effect of surfactant partitioning on the dissolution kinetics of residual perchloroethylene (PCE) by aqueous solutions (1000 mg/L) of the non-ionic surfactant Triton X-100 in a model porous medium. For this fluid system, batch equilibration experiments showed that the surfactant partitions strongly into the NAPL (NAPL-water partition coefficient equal to 12.5). Dynamic interfacial tension (IFT) measurements were employed to study surfactant diffusion and interfacial adsorption. The dynamic IFT measurements were consistent with partitioning of the surfactant between the two liquid phases. PCE dissolution experiments, conducted in a transparent glass micromodel using an aqueous surfactant solution, were contrasted to experiments using clean water. Surfactant partitioning was observed to delay significantly the onset of micellar solubilization of PCE, an observation reproduced by a numerical model. This effect is attributed to the reduction of surfactant concentration in the immediate vicinity of the NAPL-water interface, which accompanies transport of the surfactant into the NAPL. Accordingly, it is suggested that both the rate and the extent of diffusion of the surfactant into the NAPL affect the onset of and the driving force for micellar solubilization. While many surfactants do not readily partition in NAPL, this possibility must be considered when selecting non-ionic surfactants for the enhanced solubilization of residual chlorinated solvents in porous media.     

A functional relation for field-scale nonaqueous phase liquid dissolution developed using a pore network model

A pore network model with cubic chambers and rectangular tubes was used to estimate the nonaqueous phase liquid (NAPL) dissolution rate coefficient, Kdissai, and NAPL/water total specific interfacial area, ai. Kdissai was computed as a function of modified Peclet number (Pe') for various NAPL saturations (SN) and ai during drainage and imbibition and during dissolution without displacement. The largest contributor to ai was the interfacial area in the water-filled corners of chambers and tubes containing NAPL. When Kdissai was divided by ai, the resulting curves of dissolution coefficient, Kdiss versus Pe' suggested that an approximate value of Kdiss could be obtained as a weak function of hysteresis or SN. Spatially and temporally variable maps of Kdissai calculated using the network model were used in field-scale simulations of NAPL dissolution. These simulations were compared to simulations using a constant value of Kdissai and the empirical correlation of Powers et al. [Water Resour. Res. 30(2) (1994b) 321]. Overall, a methodology was developed for incorporating pore-scale processes into field-scale prediction of NAPL dissolution.     

Characterization of pore scale NAPL morphology in homogeneous sands as a function of grain size and NAPL dissolution

In this study, we investigate pore scale morphology of nonaqueous phase liquids (NAPLs) trapped in different pore sizes using tracer techniques. Specific interfacial area and saturation of NAPL trapped in homogeneous sands were measured using the interfacial and partitioning tracer techniques. The observed NAPL-water interfacial areas increased in a log-linear fashion with decreasing sand grain size, but showed no clear trend with residual NAPL saturation formed in the various grain sizes. The measured values were used to calculate the NAPL morphology index, which characterizes the spatial NAPL distribution within the pore space. The NAPL morphology indices, increased exponentially with decreasing grain size, indicating that the NAPL becomes smaller, but more blobs. For a fixed grain size, the specific interfacial area and saturation of the NAPL were measured following changes caused by dissolution using alcohol. The observed interfacial areas showed a decrease linearly as a function of the NAPL saturation.     

Evaluation of simplified mass transfer models to simulate the impacts of source zone architecture on nonaqueous phase liquid dissolution in heterogeneous porous media

Nonaqueous phase liquid (NAPL) dissolution was studied in three-dimensional (3D) heterogeneous experimental aquifers (25.5 cm x 9 cm x 8.5 cm) with two different longitudinal correlation lengths (2.1 cm and 1.1 cm) and initial spill volumes (22.5 ml and 10.5 ml). Spatial and temporal distributions of NAPL during dissolution were measured using magnetic resonance imaging (MRI). At high NAPL spill volume, average effluent concentrations initially increased during dissolution, as NAPL pools transitioned to NAPL ganglia, and then decreased as the total NAPL-water interfacial area decreased over time. Experimental results were used to test six dissolution models: (i and ii) a one-dimensional (1D) model using either specific NAPL-water interfacial area values estimated from MR images at each time step (i.e., 1D quasi-steady state model), or an empirical mass transfer (Sh') correlation (i.e., 1D transient model), (iii and iv) a multiple analytical source superposition technique (MASST) using either the NAPL distribution determined from MR images at each time step (i.e., MASST steady state model), or the NAPL distribution determined from mass balance calculations (i.e., MASST transient model), (v) an equilibrium streamtube model, and (vi) a 3D grid-scale pool dissolution model (PDM) with a dispersive mass flux term. The 1D quasi-steady state model and 3D PDM captured effluent concentration values most closely, including some concentration fluctuations due to changes in the extent of flow reduction. The 1D transient, MASST steady state and transient, and streamtube models all showed a monotonic decrease in effluent concentration values over time, and the streamtube model was the most computationally efficient. Changes during dissolution of the effective NAPL-water interfacial area estimated from imaging data are similar to changes in effluent concentration values. The 1D steady state model incorporates estimates of the effective NAPL-water interfacial area directly at each time point; the 3D PDM does so indirectly through mass balance and a relative permeability function, which causes reduced water flow through high saturation NAPL regions. Hence, when model accuracy is required, the results indicate that a surrogate of this effective interfacial area is required. Approaches to include this surrogate in the MASST and streamtube models are recommended.     

Locally-calibrated light transmission visualization methods to quantify nonaqueous phase liquid mass in porous media

Five locally-calibrated light transmission visualization (LTV) methods were tested to quantify nonaqueous phase liquid (NAPL) mass and mass reduction in porous media. Tetrachloroethylene (PCE) was released into a two-dimensional laboratory flow chamber packed with water-saturated sand which was then flushed with a surfactant solution (2% Tween 80) until all of the PCE had been dissolved. In all the LTV methods employed here, the water phase was dyed, rather than the more common approach of dyeing the NAPL phase, such that the light adsorption characteristics of NAPL did not change as dissolution progressed. Also, none of the methods used here required the use of external calibration chambers. The five visualization approaches evaluated included three methods developed from previously published models, a binary method, and a novel multiple wavelength method that has the advantage of not requiring any assumptions about the intra-pore interface structure between the various phases (sand/water/NAPL). The new multiple wavelength method is also expected to be applicable to any translucent porous media containing two immiscible fluids (e.g., water-air, water-NAPL). Results from the sand-water-PCE system evaluated here showed that the model that assumes wetting media of uniform pore size (Model C of Niemet and Selker, 2001) and the multiple wavelength model with no interface structure assumptions were able to accurately quantify PCE mass reduction during surfactant flushing. The average mass recoveries from these two imaging methods were greater than 95% for domain-average NAPL saturations of approximately 2.6x10(-2), and were approximately 90% during seven cycles of surfactant flushing that sequentially reduced the average NAPL saturation to 7.5x10(-4).     

A micromodel analysis of factors influencing NAPL removal by surfactant foam flooding

A methodology to study the trichloroethylene (TCE) and dodecane removal in porous media by surfactant foams (SF) was presented by using etched-glass micromodels. The purpose of this work was to systematically evaluate the impact of various physicochemical factors such as gas fraction (GF), surfactant concentration, pore structure and nonaqueous phase liquid (NAPL) types on NAPL removal during SF flooding. The TCE displacement by SF was dependent on the gas fraction of SF. Low GFs (50% and 66%) were more efficient for TCE removal and sweep efficiencies than a high GF (85%). An increase in TCE removal was observed with increasing surfactant concentration at a fixed GF. TCE removal by SF flooding appeared to be dependent more to the value of Capillary number rather than to the concentration of surfactant solution. The effect of the pore heterogeneity was evaluated by employing two different types of micromodels. The Capillary number is an important parameter in the determination of sweep efficiency or gas saturation of SF in a nonhomogeneous porous medium. However, the TCE removal from a nonhomogeneous porous medium may not be associated with sweep efficiency. The initial configuration of residual TCE blobs in a nonhomogeneous porous medium would also be influential in displacing TCE. Sweep efficiencies and pressure responses of two NAPL systems (TCE and dodecane) were monitored to evaluate foam stability when the foam contacts the NAPLs. Stable foam contacting with TCE is implied, while it appears that dodecane cause the SF to collapse. All results indicate that the Capillary number (a ratio of viscous forces to capillary forces) is the most important parameter for TCE removal by SF flooding. Micromodel visualizations of water, surfactant and SF floods were showed and also discussed.     

Refinement of the density-modified displacement method for efficient treatment of tetrachloroethene source zones

A novel method to remediate dense nonaqueous phase liquid (DNAPL) source zones that incorporates in situ density conversion of DNAPL via alcohol partitioning followed by displacement with a low interfacial tension (IFT) surfactant flood has been developed. Previous studies demonstrated the ability of the density-modified displacement (DMD) method to recover chlorobenzene (CB) and trichloroethene (TCE) from heterogeneous porous media without downward migration of the dissolved plume or free product. However, the extent of alcohol (n-butanol) partitioning required for in situ density conversion of high-density NAPLs, such as tetrachloroethene (PCE), could limit the utility of the DMD method. Hence, the objective of this study was to compare the efficacy of two n-butanol delivery approaches: an aqueous solution of 6% (wt) n-butanol and a surfactant-stabilized macroemulsion containing 15% (vol) n-butanol in water, to achieve density reduction of PCE-NAPL in two-dimensional (2-D) aquifer cells. Results of liquid-liquid equilibrium studies indicated that density conversion of PCE relative to water occurred at an n-butanol mole fraction of 0.56, equivalent to approximately 5 ml n-butanol per 1 ml of PCE when in equilibrium with an aqueous solution. In 2-D aquifer cell studies, density conversion of PCE was realized using both n-butanol preflood solutions, with effluent NAPL samples exhibiting density reductions ranging from 0.51 to 0.70 g/ml. Although the overall PCE mass recoveries were similar (91% and 93%) regardless of the n-butanol delivery method, the surfactant-stabilized macroemulsion preflood removed approximately 50% of the PCE mass. In addition, only 1.2 pore volumes of the macroemulsion solution were required to achieve in situ density conversion of PCE, compared to 6.4 pore volumes of the 6% (wt) n-butanol solution. These findings demonstrate that use of the DMD method with a surfactant-stabilized macroemulsion containing n-butanol holds promise as an effective source zone remediation technology, allowing for efficient recovery of PCE-DNAPL while mitigating downward migration of the dissolved plume and free product.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: dissolved-plume retardation

A natural gradient emplaced-source (ES) controlled field experiment was conducted at the Borden aquifer research site, Ontario, to study the transport of dissolved plumes emanating from residual dense nonaqueous-phase liquid (DNAPL) source zones. The specific objective of the work presented here is to determine the effects of solute and co-solute concentrations on sorption and retardation of dissolved chlorinated solvent-contaminant plumes. The ES field experiment comprised a controlled emplacement of a residual multicomponent DNAPL below the groundwater table and intensive monitoring of dissolved-phase plumes of trichloromethane (TCM), trichloroethylene (TCE), and perchloroethylene (PCE) plumes continuously generated in the aquifer down gradient from gradual source dissolution. Estimates of plume retardation (and dispersion) were obtained from 3-D numerical simulations that incorporated transient source input and flow regimes monitored during the test. PCE, the most retarded solute, surprisingly exhibited a retardation factor approximately 3 times lower than observed in a previous Borden tracer test by Mackay et al. [Water Resour. Res. 22 (1986) 2017] conducted approximately 150 m away. Also, an absence of temporal trend in PCE retardation contrasted with the previous Borden test. Supporting laboratory studies on ES site core indicated that sorption was nonlinear and competitive, i.e. reduced sorption of PCE was observed in the presence of TCE. Consideration of the effects of relatively high co-solute (TCE) concentration (competitive sorption) in addition to PCE concentration effects (nonlinear sorption) was necessary to yield laboratory-based PCE retardation estimates consistent with the field plume values. Concentration- and co-solute-based sorption and retardation analysis was also applied to the previous low-concentration pulse injection test of Mackay et al. [Water Resour. Res. 22 (1986) 2017] and was able to successfully predict the temporal field retardation trends observed in that test. While it is acknowledged that other "nonideal transport" effects may contribute, our analysis predicts differences in the PCE retardation magnitude and trend between the two experiments that are consistent with field observations based on the marked solute concentration differences that resulted from contrasting source conditions. Solute and co-solute concentration effects have heretofore received little attention, but may have wide significance in aquifers contaminated by point-source pollutants because many plumes contain mixed solutes over wide concentration ranges in strata that are likely subject to nonlinear sorption.     

Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 2. Numerical simulation

A numerical model of surfactant enhanced solubilization was developed and applied to the simulation of nonaqueous phase liquid recovery in two-dimensional heterogeneous laboratory sand tank systems. Model parameters were derived from independent, small-scale, batch and column experiments. These parameters included viscosity, density, solubilization capacity, surfactant sorption, interfacial tension, permeability, capillary retention functions, and interphase mass transfer correlations. Model predictive capability was assessed for the evaluation of the micellar solubilization of tetrachloroethylene (PCE) in the two-dimensional systems. Predicted effluent concentrations and mass recovery agreed reasonably well with measured values. Accurate prediction of enhanced solubilization behavior in the sand tanks was found to require the incorporation of pore-scale, system-dependent, interphase mass transfer limitations, including an explicit representation of specific interfacial contact area. Predicted effluent concentrations and mass recovery were also found to depend strongly upon the initial NAPL entrapment configuration. Numerical results collectively indicate that enhanced solubilization processes in heterogeneous, laboratory sand tank systems can be successfully simulated using independently measured soil parameters and column-measured mass transfer coefficients, provided that permeability and NAPL distributions are accurately known. This implies that the accuracy of model predictions at the field scale will be constrained by our ability to quantify soil heterogeneity and NAPL distribution.     

Effect of source variability and transport processes on carbon isotope ratios of TCE and PCE in two sandy aquifers

Chlorinated ethenes often migrate over extended distances in aquifers and may originate from different sources. The aim of this study was to determine whether stable carbon isotope ratios remain constant during dissolution and transport of chlorinated ethenes and whether the ratios can be used to link plumes to their sources. Detailed depth-discrete delineation of the carbon isotope ratio in a tetrachloroethene (PCE) plume and in a trichloroethene (TCE) plume was done along cross-sections orthogonal to groundwater flow in two sandy aquifers in the Province of Ontario, Canada. At the TCE site, TCE concentrations up to solubility were measured in one high concentration zone close to the bottom of the aquifer from where dense non-aqueous phase liquid (DNAPL) was collected. A laboratory experiment using the DNAPL indicated that only very small carbon isotope fractionation occurs during dissolution of TCE (0.26 per thousand), which is consistent with field observations. At most sampling points, the delta(13)C of dissolved TCE was similar to that of the DNAPL except for a few sampling points at the bottom of the aquifer close to the underlying aquitard. At these points, a (13)C enrichment of up to 2.4 per thousand was observed, which was likely due to biodegradation and possibly preferential diffusion of TCE with (12)C into the aquitard. In contrast to the TCE site, several distinct zones of high concentration were observed at the PCE site and from zones to zone, the delta(13)C values varied substantially from -24.3 per thousand to -33.6 per thousand. Comparison of the delta(13)C values in the high concentration zones made it possible to divide the plume in the three different domains, each probably representing a different episode and location of DNAPL release. The three different zones could still be distinguished 220 m from the DNAPL sources. This demonstrates that carbon isotope ratios can be used to differentiate between different zones in chlorinated ethene plumes and to link plume zones to their sources. In addition, subtle variations in delta(13)C at plume fringes provided insight into mechanisms of plume spreading in transverse vertical direction. These variations were identified because of the high-resolution provided by the monitoring network.     

The dissolution of benzene, toluene, m-xylene and naphthalene from a residually trapped non-aqueous phase liquid under mass transfer limited conditions

The results of dissolution experiments for benzene, toluene, m-xylene and naphthalene (BTXN) from a relatively insoluble oil phase (tridecane), residually trapped in a non-sorbing porous medium, are described. This mixture was chosen to simulate dissolution of soluble aromatic compounds from a petroleum hydrocarbon mixture, e.g., crude oil, for which a large fraction of the mixture is relatively insoluble. The experiments were carried out at a small source length to interstitial velocity ratio, L/v, so that dissolution would be mass transfer limited (MTL). When fitted to data for toluene, a multiregion mass transfer model was found to predict the experimental data satisfactorily for the other components without adjustment of the mass transfer rate parameters. These results indicate that the dissolution process can be generalized for various hydrophobic organic compounds present in a multicomponent non-aqueous phase liquid (NAPL) when mass transfer limitations are present. This also suggests that dissolution data obtained for one compound can be useful for predicting the dissolution histories for other compounds from petroleum hydrocarbon mixtures.     

Comparison of PCE and TCE disappearance in heated volatile organic analysis vials and flame-sealed ampules

The rates of hydrolysis reported for tetrachloroethylene (PCE) and trichloroethylene (TCE) at elevated temperatures range over two orders-of-magnitude, where some of the variability may be due to the presence of a gas phase. Recent studies suggest that volatile organic analysis (VOA) vials provide a low-cost and readily available zero headspace system for measuring aqueous-phase hydrolysis rates. This work involved measuring rates of PCE and TCE disappearance and the corresponding appearance of dechlorination products in water-filled VOA vials and flame-sealed ampules incubated at 21 and 55 °C for up to 95.5 days. While PCE and TCE concentrations readily decreased in the VOA vials to yield first-order half lives of 11.2 days for PCE and 21.1 days for TCE at 55 °C, concentrations of anticipated dechlorination products, including chloride, remained constant or were not detected. The rate of PCE disappearance was 34 times faster in VOA vials at 55 °C compared to values obtained with flame-sealed ampules containing PCE-contaminated water. In addition, the concentration of TCE increased slightly in flame-sealed ampules incubated at 55 °C, while a decrease in TCE levels was observed in the VOA vials. The observed losses of PCE and TCE in the VOA vials were attributed to diffusion and sorption in the septa, rather than to dechlorination. These findings demonstrate that VOA vials are not suitable for measuring rates of volatile organic compound hydrolysis at elevated temperatures.     

Mass flux from a non-aqueous phase liquid pool considering spontaneous expansion of a discontinuous gas phase

The partitioning of non-aqueous phase liquid (NAPL) compounds to a discontinuous gas phase results in the repeated spontaneous expansion, snap-off, and vertical mobilization of the gas phase. This mechanism has the potential to significantly affect the mass transfer processes that control the dissolution of NAPL pools by increasing the vertical transport of NAPL mass and increasing the total mass transfer rate from the surface of the pool. The extent to which this mechanism affects mass transfer from a NAPL pool depends on the rate of expansion and the mass of NAPL compound in the gas phase. This study used well-controlled bench-scale experiments under no-flow conditions to quantify for the first time the expansion of a discontinuous gas phase in the presence of NAPL. Air bubbles placed in glass vials containing NAPL increased significantly in volume, from a radius of 1.0 mm to 2.0 mm over 215 days in the presence of tetrachloroethene (PCE), and from a radius of 1.2 mm to 2.3 mm over 22 days in the presence of trans-1,2-dichloroethene (tDCE). A one-dimensional mass transfer model, fit to the experimental data, showed that this expansion could result in a mass flux from the NAPL pool that was similar in magnitude to the mass flux expected for the dissolution of a NAPL pool in a two-fluid (NAPL and water) system. Conditions favouring the significant effect of a discontinuous gas phase on mass transfer were identified as groundwater velocities less than ~0.01 m/day, and a gas phase that covers greater than ~10% of the pool surface area and is located within ~0.01 m of the pool surface. Under these conditions the mass transfer via a discontinuous gas phase is expected to affect, for example, efforts to locate NAPL source zones using aqueous concentration data, and predict the lifetime and risk associated with NAPL source zones in a way that is not currently included in the common conceptual models used to assess NAPL-contaminated sites.     

Bioenhancement of NAPL pool dissolution: experimental evaluation

Experiments were conducted to quantify nonaqueous phase liquid (NAPL) pool dissolution and its enhancement by in situ biodegradation. The experiments were performed using square cross-section, glass-bead packed column reactors with a small pool of a toluene-in-dodecane mixture (toluene mole fraction, X(tol) approximately 0.02 or 0.09). Experimental quasi-steady-state toluene dissolution fluxes were determined using a 14C-carbon mass-balance approach during water flushing with and without biodegradation. The experiments demonstrated a statistically significant bioenhancement of the toluene dissolution flux of up to roughly twofold at average pore water velocities of approximately 0.1 and 1 m/day when the toluene mole fraction was low ( approximately 0.02); however, little or no bioenhancement was observed with the higher mole fraction ( approximately 0.09). Although it cannot be determined conclusively, the weight of evidence based on biomass measurements and model analyses suggests that the reduced bioenhancement for the high mole fraction was due to higher dissolved toluene concentrations, which may have caused toxicity effects. Importantly, even though NAPL dissolution was not bioenhanced in every case, the biodegradation reduced toluene concentrations to low levels in the reactor effluents.     

A three-layer diffusion-cell to examine bio-enhanced dissolution of chloroethene dense non-aqueous phase liquid

Microbial reductive dechlorination of trichloroethene (TCE) and perchloroethene (PCE) in the vicinity of their dense non-aqueous phase liquid (DNAPL) has been shown to accelerate DNAPL dissolution. A three-layer diffusion-cell was developed to quantify this bio-enhanced dissolution and to measure the conditions near the DNAPL interface. The 12 cm long diffusion-cell setup consists of a 5.5 cm central porous layer (sand), a lower 3.5 cm DNAPL layer and a top 3 cm water layer. The water layer is frequently refreshed to remove chloroethenes at the upper boundary of the porous layer, while the DNAPL layer maintains the saturated chloroethene concentration at the lower boundary. Two abiotic and two biotic diffusion-cells with TCE DNAPL were tested. In the abiotic diffusion-cells, a linear steady state TCE concentration profile between the DNAPL and the water layer developed beyond 21 d. In the biotic diffusion-cells, TCE was completely converted into cis-dichloroethene (cis-DCE) at 2.5 cm distance of the DNAPL. Dechlorination was likely inhibited up to a distance of 1.5 cm from the DNAPL, as in this part the TCE concentration exceeded the culture’s maximum tolerable concentration (2.5 mM). The DNAPL dissolution fluxes were calculated from the TCE concentration gradient, measured at the interface of the DNAPL layer and the porous layer. Biotic fluxes were a factor 2.4 (standard deviation 0.2) larger than abiotic dissolution fluxes. This diffusion-cell setup can be used to study the factors affecting the bio-enhanced dissolution of DNAPL and to assess bioaugmentation, pH buffer addition and donor delivery strategies for source zones.     

Characterization of DNAPL from the U.S. DOE Savannah River Site

The composition of chlorinated hydrocarbon DNAPLs (dense non-aqueous phase liquids) from field sites can be substantially different than the material originally purchased for use as a solvent. Waste management practices at the U.S. Department of Energy's (DOE) Savannah River Site (SRS) included co-disposal of a wide range of organic and inorganic wastes. In 1991, a clear, orange-colored DNAPL was found in two wells near the SRS M-area settling basin. Waste effluent from the fuel and target fabrication facilities that were discharged to this settling basin included acids, caustics, metals and chlorinated solvents. The characterization of the SRS DNAPL suggests that numerous constituents partitioned into the DNAPL during its use as a solvent, co-disposal and ultimate migration through the subsurface. Trace constituents in the DNAPL include metals, from processing operations or co-disposal practices and subsurface minerals, high molecular weight hydrocarbons and alkyl esters, and acids. This complex mixture results in DNAPL-water interfacial properties that are substantially different than would be expected from a simple mixture of PCE and TCE. Under conditions when there is a high DNAPL to water volume ratio, a semi-rigid film accumulates on water droplets suspended in the DNAPL. It is concluded that the array of precipitated metal species comprising this film contributes to the interfacial tension that is over an order of magnitude lower than expected for a "clean" PCE/TCE mixture.     

零价铁修复受三氯乙烯污染地下水的实验研究

通过批实验和柱实验研究了三氯乙烯（TCE）初始浓度、四氯乙烯(PCE)等对零价铁去除三氯乙烯的影响，并建立了三氯乙烯降解的反应动力学方程。结果表明：（1）零价铁对TCE具有较好的降解效果，反应符合准一级反应动力学方程，表观反应速率常数随TCE浓度的增加而减小；（2）在铁粉充足的条件下，TCE初始浓度对降解效果影响不显著，且TCE去除率皆可达到90%以上；（3）PCE的存在抑制了TCE的脱氯反应。PCE和TCE共存时，TCE的最大去除率仅为64.2%；TCE脱氯反应的表观反应速率明显降低，反应半衰期由TCE单独存在时的6.8～9.7 h增大到66 h～346.5 h。     

Secondary imbibition in NAPL-invaded mixed-wet sediments

A previously developed pore network model is used here to study the spontaneous and forced secondary imbibition of a NAPL-invaded sediment, as in the displacement of NAPL by waterflooding a mixed-wet soil. We use a 3D disordered pore network with a realistic representation of pore geometry and connectivity, and a quasi-static displacement model that fully describes the pore-scale physics. After primary drainage (NAPL displacing water) up to a maximum capillary pressure, we simulate secondary imbibition (water displacing NAPL). We conduct a parametric study of imbibition by varying systematically the controlling parameters: the advancing contact angles, the fraction of NAPL-wet pores, the interfacial tension, and the initial water saturation. Once the secondary imbibition is completed, the controlling displacement mechanisms, capillary pressures, relative permeabilities, and trapped NAPL saturations are reported. It is assumed that NAPL migrates into an initially strongly water-wet sediment, i.e., the receding contact angles are very small. However, depending on the surface mineralogy and chemical compositions of the immiscible fluid phases, the wettability of pore interiors is altered while the neighborhoods of pore corners remain strongly water-wet-resulting in a mixed-wet sediment. Here, we compare three different levels of wettability alteration: water-wet (advancing contact angles (20 degrees to 55 degrees), intermediate-wet (55 degrees to 120 degrees), and NAPL-wet (120 degrees to 155 degrees). The range of advancing contact angles and the fraction of NAPL-wet pores have dramatic effects on the NAPL-water capillary pressures and relative permeabilities. The spatially inhomogeneous interfacial tension has a minor impact on the trapped NAPL saturation and relative permeability to NAPL, and a slight effect on the relative permeability to water. The initial water saturation has a slight effect on the two-phase flow characteristics of water-wet sediments; however, with more NAPL-wet pores in the sediment, it starts to have a profound effect on the water and NAPL relative permeabilities.     

Push-pull partitioning tracer tests using radon-222 to quantify non-aqueous phase liquid contamination

Naturally occurring radon in groundwater can be used as an in situ partitioning tracer for locating and quantifying non-aqueous phase liquid (NAPL) contamination in the subsurface. When combined with the single-well, push-pull test, this methodology has the potential to provide a low-cost alternative to inter-well partitioning tracer tests. During a push-pull test, a known volume of test solution (radon-free water containing a conservative tracer) is first injected ("pushed") into a well; flow is then reversed and the test solution/groundwater mixture is extracted ("pulled") from the same well. In the presence of NAPL radon transport is retarded relative to the conservative tracer. Assuming linear equilibrium partitioning, retardation factors for radon can be used to estimate NAPL saturations. The utility of this methodology was evaluated in laboratory and field settings. Laboratory push-pull tests were conducted in both non-contaminated and trichloroethene NAPL (TCE)-contaminated sediment. The methodology was then applied in wells located in non-contaminated and light non-aqueous phase liquid (LNAPL)-contaminated portions of an aquifer at a former petroleum refinery. The method of temporal moments and an approximate analytical solution to the governing transport equations were used to interpret breakthrough curves and estimate radon retardation factors; estimated retardation factors were then used to calculate TCE saturations. Numerical simulations were used to further investigate the behavior of the breakthrough curves. The laboratory and field push-pull tests demonstrated that radon retardation does occur in the presence of TCE and LNAPL and that radon retardation can be used to calculate TCE saturations. Laboratory injection-phase test results in TCE-contaminated sediment yielded radon retardation factors ranging from 1.1 to 1.5, resulting in calculated TCE saturations ranging from 0.2 to 0.9%. Laboratory extraction-phase test results in the same sediment yielded a radon retardation factor of 5.0, with a calculated TCE saturation of 6.5%. Numerical simulation breakthrough curves provided reasonably good matches to the approximate analytical solution breakthrough curves. However, non-equilibrium radon partitioning and heterogeneous TCE distributions may affect the retardation factors and TCE saturation estimates.     

A set of constitutive relationships accounting for residual NAPL in the unsaturated zone

Although laboratory experiments show that non-aqueous phase liquid (NAPL) is retained in the unsaturated zone, no existing multiphase flow model has been developed to account for residual NAPL after NAPL drainage in the unsaturated zone. We developed a static constitutive set of saturation-capillary pressure relationships for water, NAPL and air that accounts for both this residual NAPL and entrapped NAPL. The set of constitutive relationships is formulated similarly to the set of scaled relationships that is frequently applied in continuum models. The new set consists of three fluid-phase systems: a three-phase system and a two-phase system, that both comply with the original constitutive model, and a newly introduced residual NAPL system. The new system can be added relatively easily to the original two- and three-phase systems. Entrapment is included in the model. The constitutive relationships of the non-drainable residual NAPL system are based on qualitative fluid behavior derived from a pore scale model. The pore scale model reveals that the amount of residual NAPL depends on the spreading coefficient and the water saturation. Furthermore, residual NAPL is history-dependent. At the continuum scale, a critical NAPL pressure head defines the transition from free, mobile NAPL to residual NAPL. Although the Pc-S relationships for water and total liquid are not independent in case of residual NAPL, two two-phase Pc-S relations can represent a three-phase residual system of Pc-S relations. A newly introduced parameter, referred to as the residual oil pressure head, reflects the mutual dependency of water and oil. Example calculations show consistent behavior of the constitutive model. Entrapment and retention in the unsaturated zone cooperate to retain NAPL. Moreover, the results of our constitutive model are in agreement with experimental observations.