Explicit and implicit coupling during solute transport through clay membrane barriers

Simulations of salt (KCl) flux through a 1-m-thick clay membrane barrier (CMB) based on coupled solute transport theory are compared to simulated fluxes based on traditional advective-dispersive transport theory. The simulations are based on measured values for the effective salt-diffusion coefficient (Ds*) and chemico-osmotic efficiency coefficient (omega) for a bentonite-based barrier material subjected to KCl solutions. The results indicate that the exit salt flux is reduced due to both explicit coupling (hyperfiltration and chemico-osmotic counter-advection) and an implicit coupling effect resulting from the decrease in Ds* due to a decrease in the apparent tortuosity factor, tau a, with an increase in omega. Implicit coupling is shown to be more significant than explicit coupling for reducing and retarding salt flux through a CMB under diffusion-dominated conditions. Failure to account for the implicit coupling effect may result in unrealistic results, such as the existence of salt flux through a perfect (ideal) clay membrane (i.e., omega=1).     

Theory for reactive solute transport through clay membrane barriers

The theoretical development for one-dimensional, coupled migration of solutes with different ionic mobilities through clay soils that behave as ion-restrictive membranes, referred to as clay membrane barriers (CMBs), is presented. The transport formulation is based on principles of irreversible thermodynamics and accounts explicitly for coupling effects of hyperfiltration (ultrafiltration) and chemico-osmotic counter-advection associated with clay membrane behavior in the absence of electrical current. Since, by definition, no solute can enter a "perfect" or "ideal" membrane, the concept of an implicit coupling effect, such that the effective salt-diffusion coefficient, Ds* approaches zero as the chemico-osmotic efficiency coefficient, omega approaches unity is introduced. The theoretical development also illustrates that, even in the absence of membrane behavior, traditional advective-dispersive transport theory based on a constant value of Ds* for the solutes may not be appropriate for simulating transient transport in reactive (ion exchanging) systems. This potential limitation is illustrated through simulations for solute mass flux involving the migration of a binary salt solution (KCl) through a clay barrier with exchange sites saturated with a single exchangeable cation (e.g., Na+) that enters the pore solution upon ion exchange with the salt cation (K+).     

Sorption of tylosin on clay minerals

The equilibrium sorption of tylosin (TYL) on kaolinite and montmorillonite was measured at different solution pH using batch reactor systems. The results showed that all the sorption isotherms were nonlinear and that the nonlinearity decreased as the solution pH increased for a given clay. At a specific aqueous concentration, the single-point sorption distribution coefficient (K_D) of TYL decreased rapidly as the solution pH increased. A speciation-dependent sorption model that accounted for the contributions of the cationic and neutral forms of TYL fit the data well, suggesting that the sorption may be dominated by both ion exchange and hydrophobic interactions. The isotherm data also fit well to a dual mode model that quantifies the contributions of a site-limiting Langmuir component (ion exchange) and a non-specific linear partitioning component (hydrophobic interactions). X-ray diffraction analyses revealed that the interlayers of montmorillonite were expanded due to the uptake of TYL. TYL molecules likely form a monolayer surface coverage.     

Availability and bio-accessibility of metals in the clay fraction of urban soils of Sevilla

The availability of Cd, Cr, Cu, Ni, Mn, Pb and Zn present in the finest size particles of urban soils is studied by comparing the concentrations in the clay fraction with those extracted from the whole soil by either single-extraction or sequential extraction method. Many metals are preferentially present in the finest particles as compared to coarser fractions. This is true for most metals studied, except Mn and, perhaps, Cd. Those metals present in the clay fraction are often in easily bio-accessible forms, especially Cu, Pb and Zn. The results suggest that bio-accessible forms of these three metals are distributed among the three sequential fractions, and even the fraction considered as ‘residual’ is also bio-accessible to a significant extent. The statistical analysis shows some distinctions among metals that are compared to the ‘urban’, ‘natural’, or intermediate behaviour of the various metals as proposed earlier in the literature.     

Magnetic resonance imaging of nonaqueous phase liquid during soil vapor extraction in heterogeneous porous media

Soil vapor extraction (SVE) is commonly used to remediate nonaqueous phase liquids (NAPLs) from the vadose zone. This paper aims to determine the effect of grain size heterogeneity on the removal of NAPL in porous media during SVE. Magnetic resonance imaging (MRI) was used to observe and quantify the amount and location of NAPL in flow-through columns filled with silica gel grains. MRI is unique because it is nondestructive, allowing three-dimensional images to be taken of the phases as a function of space and time. Columns were packed with silica gel in three ways: coarse grains (250-550 microm) only, fine grains (32-63 microm) only, and a core of fine grains surrounded by a shell of coarse grains. Columns saturated with water were drained under a constant suction head, contaminated with decane, and then drained to different decane saturations. Each column was then continuously purged with water-saturated nitrogen gas and images were taken intermittently. Results showed that at residual saturation, a sharp volatilization front moved through the columns filled with either coarse-grain or fine-grain silica gel. In the heterogeneous columns, the volatilization front in the core lagged just behind the shell because gas flow was greater through the shell and decane in the core diffused outward to the shell. When decane saturation in the core was above residual saturation, decane volatilization occurred near the inlet, the relative decane saturation throughout the core dropped uniformly, and decane in the core flowed in the liquid phase to the shell to replenish volatilized decane. These results indicate that NAPL trapped in low-permeability zones can flow to replenish areas where NAPL is lost due to SVE. However, when residual NAPL saturation is reached, NAPL flow no longer occurs and diffusion limits removal from low-permeability zones.     

An analytical model for the interpretation of pulse injection experiments performed for testing the spatial variability of clay formations

This paper presents an analytical model to describe pulse injection experiments. This model solves the advection-diffusion equation while taking into account back diffusion from the clay core to the inlet and from the outlet to the clay core. In most analytical models, back diffusion is neglected. For sufficiently high Péclet numbers, this is a good approximation. However, in experiments where the Péclet number is low, back diffusion is important and must be taken into account. An additional advantage of the present model is that both concentration and flux are conserved at the inlet and at the outlet of the clay core. This model is used to fit pulse injection experiments with iodide and tritiated water (HTO) in clay cores. The (new) model is required for fitting the experimental results since in clay layers advection is very slow leading to a low Péclet number. The experiments are performed on clay cores taken from different depths from the Boom Clay and the Ypres Clay layer under the site of the nuclear power plant of Doel (Belgium). The quality of all fits is excellent and the obtained parameter values are coherent. For HTO, the fitted value for the diffusion accessible porosity is consistent with measurements of the water content in Ypres Clay cores. In both types of clays, the apparent diffusion coefficient at zero flow is between 10(-10) and 2 x 10(-10) m(2)/s for iodide and between 2 x 10(-10) and 3 x 10(-10) m(2)/s for HTO. The dispersion length is in the order of 10(-3) m. The average value for the diffusion accessible porosity is between 0.35 and 0.4 for HTO and between 0.2 and 0.25 for iodide.     

Obtaining the porewater composition of a clay rock by modeling the in- and out-diffusion of anions and cations from an in-situ experiment

A borehole in the Callovo–Oxfordian clay rock in ANDRA's underground research facility was sampled during 1 year and chemically analyzed. Diffusion between porewater and the borehole solution resulted in concentration changes which were modeled with PHREEQC's multicomponent diffusion module. In the model, the clay rock's pore space is divided in free porewater (electrically neutral) and diffuse double layer water (devoid of anions). Diffusion is calculated separately for the two domains, and individually for all the solute species while a zero-charge flux is maintained. We explain how the finite difference formulas for radial diffusion can be translated into mixing factors for solutions. Operator splitting is used to calculate advective flow and chemical reactions such as ion exchange and calcite dissolution and precipitation. The ion exchange reaction is formulated in the form of surface complexation, which allows distributing charge over the fixed sites and the diffuse double layer. The charge distribution affects pH when calcite dissolves, and modeling of the experimental data shows that about 7% of the cation exchange capacity resides in the diffuse double layer. The model calculates the observed concentration changes very well and provides an estimate of the pristine porewater composition in the clay rock.     

Characterisation of HTO diffusion properties by an in situ tracer experiment in Opalinus clay at Mont Terri

A long-term single borehole diffusion experiment using tritiated water as tracer was carried out in Opalinus clay, an argillaceous rock formation that is accessible at the Mont Terri Underground Research Laboratory, situated in the Swiss Jura. The tracer was diluted in reconstituted formation water and introduced into a packed-off section of a borehole located in saturated rock. Pressure in this interval was maintained equal to the pore pressure of the surrounding rock in order to prevent any hydraulic gradient around the borehole and to avoid advective transport processes. The evolution of the tracer concentration in the injection system was monitored over time. After 1 year of diffusion, the claystone surrounding the interval was retrieved by overcoring the whole borehole and packer system, and by an adjacent oblique borehole. Compressed air was used as drilling fluid to reduce rock disturbances. The recovered overcore was sampled along profiles perpendicular to the borehole wall with a view to determining the tracer-concentration profiles in the rock. To avoid further evaporation of tritiated water, subsamples were immediately transferred into polyethylene bottles and disaggregated by adding a known amount of tracer-free water. Fifteen profiles were determined and showed a decreasing tracer concentration with distance into the rock. The pore-water contents were constant along those profiles, confirming that only very little water was lost during overcoring operations. The evolution of tritium-tracer concentration in the injection system over time and in situ profiles were interpreted with a 3-D numerical simulation of the experiment. That allowed for the identification of the transport parameters (orthotropic diffusion tensor and porosity) by minimising the relative quadratic error between the experimental and simulated data. The fitting is good and the results are consistent with data obtained on drill-core samples. The result of tritiated water is discussed regarding (1) the potential effect of mechanical and/or chemical disturbances around the injection borehole and (2) the specific behaviour of tritiated water.     

Design of aquifer remediation systems: (2) estimating site-specific performance and benefits of partial source removal

A Lagrangian stochastic model is proposed as a tool that can be utilized in forecasting remedial performance and estimating the benefits (in terms of flux and mass reduction) derived from a source zone remedial effort. The stochastic functional relationships that describe the hydraulic "structure" and non-aqueous phase liquid (NAPL) "architecture" have been described in a companion paper (Enfield, C.G., Wood, A.L., Espinoza, F.P., Brooks, M.C., Annable, M., Rao, P.S.C., this issue. Design of aquifer remediation systems: (1) describing hydraulic structure and NAPL architecture using tracers. J. Contam. Hydrol.). The previously defined functions were used along with the properties of the remedial fluids to describe remedial performance. There are two objectives for this paper. First, is to show that a simple analytic element model can be used to give a reasonable estimate of system performance. This is accomplished by comparing forecast performance to observed performance. The second objective is to display the model output in terms of change in mass flux and mass removal as a function of pore volumes of remedial fluid injected. The modelling results suggest that short term benefits are obtained and related to mass reduction at the sites where the model was tested.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

A partially coupled, fraction-by-fraction modelling approach to the subsurface migration of gasoline spills

The subsurface spreading behaviour of gasoline, as well as several other common soil- and groundwater pollutants (e.g. diesel, creosote), is complicated by the fact that it is a mixture of hundreds of different constituents, behaving differently with respect to e.g. dissolution, volatilisation, adsorption and biodegradation. Especially for scenarios where the non-aqueous phase liquid (NAPL) phase is highly mobile, such as for sudden spills in connection with accidents, it is necessary to simultaneously analyse the migration of the NAPL and its individual components in order to assess risks and environmental impacts. Although a few fully coupled, multi-phase, multi-constituent models exist, such models are highly complex and may be time consuming to use. A new, somewhat simplified methodology for modelling the subsurface migration of gasoline while taking its multi-constituent nature into account is therefore introduced here. Constituents with similar properties are grouped together into eight fractions. The migration of each fraction in the aqueous and gaseous phases as well as adsorption is modelled separately using a single-constituent multi-phase flow model, while the movement of the free-phase gasoline is essentially the same for all fractions. The modelling is done stepwise to allow updating of the free-phase gasoline composition at certain time intervals. The output is the concentration of the eight different fractions in the aqueous, gaseous, free gasoline and solid phases with time. The approach is evaluated by comparing it to a fully coupled multi-phase, multi-constituent numerical simulator in the modelling of a typical accident-type spill scenario, based on a tanker accident in northern Sweden. Here the PCFF method produces results similar to those of the more sophisticated, fully coupled model. The benefit of the method is that it is easy to use and can be applied to any single-constituent multi-phase numerical simulator, which in turn may have different strengths in incorporating various processes. The results demonstrate that the different fractions have significantly different migration behaviours and although the methodology involves some simplifications, it is a considerable improvement compared to modelling the gasoline constituents completely individually or as one single mixture.     

Estimation of primary drainage three-phase relative permeability for organic liquid transport in the vadose zone

The modeling of transport of organic liquid contaminants through the vadose zone often requires three-phase relative permeabilities. Since these are difficult to measure, predictive models are usually used. The objective of this study is to assess the ability of eight common models to predict the drainage relative permeability to oil in a three-phase system (water-oil-air). A comparison of the models' estimates using data set from Oak [Oak, M.J., 1990. Three-phase relative permeability of water-wet Berea. In: Seventh Symposium on Enhanced Oil Recovery, Paper SPE/Doe 20183. Tulsa, OK, April 22-25] showed that they provide very different predictions for the same system. The goodness of the models does not increase with the amount of data or computation that the models require. Also, the calculations showed how different interpretations of the models and of the terminology associated with them can significantly impact the predictions. Thus, considerable error may be introduced into the simulations of organic liquid transport in the vadose zone depending on the selection and interpretation of the three-phase relative permeability model.     

Partitioning of aromatic and oxygenated constituents into water from regular and ethanol-blended gasolines

Variability in gasoline-water partitioning of major aromatic constituents (benzene, toluene, ethylbenzene, and xylenes (BTEX)) and methyl tert-butyl ether (MTBE) were examined for regular and ethanol-blended gasolines. By use of a two-phase liquid-liquid equilibrium model, the distribution of nonpolar solutes between fuel phase and water was related to principles of equilibrium. The models derived using Raoult's law convention for activity coefficients and liquid solubility is presented. The observed inverse log-log linear dependence of K_fw values on aqueous solubility, could be well predicted by assuming gasoline to be an ideal solvent mixture. Oxygenated additives (i.e., ethanol and MTBE), in the low percent range (below 5%), were shown to have minimal or negligible cosolvent effects on hydrocarbon partitioning. In the case of high fuel-to-water ratio (e.g., 1:1) or near contaminant source zone, the cosolvent effect of oxygenated gasoline with high content of ethanol (e.g., E85) will be environmentally significant.     

SEAR技术修复土壤和地下水中NAPL污染的研究进展

就SEAR技术修复土壤及地下水中NAPL污染的原理及发展现状进行了综述.SEAR技术可以快速有效地去除土壤和地下水中的NAPL污染源,适于多种污染物.该技术通过增溶和增流2种途径提高NAPL污染物的去除率.表面活性剂的选择和微乳液体系的调配是SEAR技术实施的关键环节.将SEAR技术用于高浓度NAPL污染源的治理,并与生物修复和自然降解相结合,是经济高效的治理方案.     

SEAR技术修复土壤和地下水中NAPL污染的研究进展

就SEAR技术修复土壤及地下水中NAPL污染的原理及发展现状进行了综述.SEAR技术可以快速有效地去除土壤和地下水中的NAPL污染源,适于多种污染物.该技术通过增溶和增流2种途径提高NAPL污染物的去除率.表面活性剂的选择和微乳液体系的调配是SEAR技术实施的关键环节.将SEAR技术用于高浓度NAPL污染源的治理,并与生物修复和自然降解相结合,是经济高效的治理方案.     

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.     

Investigation of surfactant-enhanced mass removal and flux reduction in 3D correlated permeability fields using magnetic resonance imaging

Magnetic resonance imaging (MRI) was used to visualize the NAPL source zone architecture before and after surfactant-enhanced NAPL dissolution in three-dimensional (3D) heterogeneously packed flowcells characterized by different longitudinal correlation lengths: 2.1 cm (aquifer 1) and 1.1 cm (aquifer 2). Surfactant flowpaths were determined by imaging the breakthrough of a paramagnetic tracer (MnCl(2)) analyzed by the method of moments. In both experimental aquifers, preferential flow occurred in high permeability materials with low NAPL saturations, and NAPL was preferentially removed from the top of the aquifers with low saturation. Alternate flushing with water and two surfactant pulses (5-6 pore volumes each) resulted in approximately 63% of NAPL mass removal from both aquifers. However, overall reduction in mass flux (Mass Flux 1) exiting the flowcell was lower in aquifer 2 (68%) than in aquifer 1 (81%), and local effluent concentrations were found to increase by as high as 120 times at local sampling ports from aquifer 2 after surfactant flushing. 3D MRI images of NAPL revealed that NAPL migrated downward and created additional NAPL source zones in previously uncontaminated areas at the bottom of the aquifers. The additional NAPL source zones were created in the direction transverse to flow in aquifer 2, which explains the higher mass flux relative to aquifer 1. Analysis using a total trapping number indicates that mobilization of NAPL trapped in the two coarsest sand fractions is possible when saturation is below 0.5 and 0.4, respectively. Results from this study highlight the potential impacts of porous media heterogeneity and NAPL source zone architecture on advanced in-situ flushing technologies.     

Analysis of a gas-phase partitioning tracer test conducted in an unsaturated fractured-clay formation

The gas-phase partitioning tracer method was used to estimate non-aqueous phase liquid (NAPL), water, and air saturations in the vadose zone at a chlorinated-solvent contaminated field site in Tucson, AZ. The tracer test was conducted in a fractured-clay system that is the confining layer for the underlying regional aquifer. Three suites of three tracers were injected into wells located 14, 24, and 24 m from a single, central extraction well. The tracers comprised noble gases (traditionally thought to be nonsorbing), alkanes (primarily water partitioning), perfluorides (primarily NAPL partitioning), and halons (both NAPL and water partitioning). Observations of vacuum response were consistent with flow in a fractured system. The halon tracers exhibited the greatest amount of retardation, and helium and the perfluoride tracers the least. The alkane tracers were unexpectedly more retarded than the perfluoride tracers, indicating low NAPL saturations and high water saturations. An NAPL saturation of 0.01, water saturation of 0.215, and gas saturation of 0.775 was estimated based on analysis of the suite of tracers comprising helium, perfluoromethylcyclohexane and dibromodifluoromethane, which was considered to be the most robust set. The estimated saturations compare reasonably well to independently determined values.     

Infiltration of PCE in a system containing spatial wettability variations

A two-dimensional infiltration experiment was conducted to investigate and quantify the effect of spatial wettability variations on DNAPL migration and entrapment in saturated sands. Experimental observations of tetrachloroethylene (PCE) infiltration showed that organic-wet sand lenses acted as very effective capillary barriers, retaining PCE and inhibiting its downward migration. A multiphase numerical simulator was used to model this sand box experiment. The simulator incorporates wettability-modified van Genuchten and Brooks-Corey capillary pressure/saturation relationships as well as Burdine and Mualem relative permeability relationships. PCE mass distributions, estimated by image analysis of digital photographs taken during the infiltration event, were compared to simulation results. Although both relative permeability models were qualitatively able to predict the PCE retention in the organic-wet layers, simulations with the Mualem model failed to capture the observed rate of PCE migration. A traditional multiphase simulator, incorporating water-wet capillary retention relations, failed to predict both PCE pathways and retention behavior. This study illustrates the potential influence of subsurface wettability variations on DNAPL migration and entrapment and supports the use of modified capillary relations in conjunction with the Burdine model in multiphase flow simulators.     

Hot water flushing for immiscible displacement of a viscous NAPL

Thermal remediation techniques, such as hot water flooding, are emerging technologies that have been proposed for the removal of nonaqueous phase liquids (NAPLs) from the subsurface. In this study a combined laboratory and modeling investigation was conducted to determine if hot water flooding techniques would improve NAPL mass removal compared to ambient temperature water flushing. Two experiments were conducted in a bench scale two-dimensional sandbox (55 cmx45 cmx1.3 cm) and NAPL saturations were quantified using a light transmission apparatus. In these immiscible displacement experiments the aqueous phase, at 22 degrees C and 50 degrees C, displaced a zone with initial NAPL saturations on the order of 85%. The interfacial tension and viscosity of the selected light NAPL, Voltesso 35, are strongly temperature-dependent. Experimental results suggest that hot water flooding reduced the size of the high NAPL saturation zone, in comparison to the cold water flood, and yielded greater NAPL mass recovery (75% NAPL removal vs. 64%). Hot water flooding did not, however, result in lower residual NAPL saturations. A numerical simulator was modified to include simultaneous flow of water and organic phases, energy transport, temperature and pressure. Model predictions of mass removal and NAPL saturation profiles compared well with observed behavior. A sensitivity analysis indicates that the utility of hot water flooding improves with the increasing temperature dependence of NAPL hydraulic properties.