[1] The influence of aquifer property correlation on multiphase fluid migration and entrapment was explored through the use of correlated and uncorrelated porosity, permeability, and capillary pressure-saturation (P c -Sat) parameter fields in a crosssectional numerical multiphase flow model. Data collected from core samples in a nonuniform sandy aquifer were used to generate three-dimensional aquifer parameter fields. Porosity was assumed to be uniform or simulated using sequential Gaussian simulation (SGS). Permeability (k) was modeled independently of porosity using SGS as well as simulated geostatistical indicator classes derived from measured grain size distribution curves. Retention characteristics were assigned employing Leverett scaling of a representative P c -Sat curve to the geostatistical k fields or, alternatively, on the basis of simulated indicator classes and porosity values. Ensemble dense nonaqueous phase liquid (DNAPL) infiltration and entrapment behavior for a hypothetical tetrachloroethylene (PCE) spill was simulated in four sets of two-dimensional profiles extracted from these realizations. Comparisons of saturation profiles and spatial moments from point source DNAPL infiltration simulations suggest that choices involving the geostatistical algorithm used to model k and the incorporation of variable versus uniform porosity have a smaller influence than choices involving the scaling of capillary retention properties to k. From these simulations it is apparent that the degree of spatial correlation in P c -Sat parameters exerts a controlling influence on predicted DNAPL spreading and redistribution in saturated aquifers. The resultant distribution of mass within a DNAPL source zone will have implications for DNAPL recovery and subsequent mass fluxes in remediation operations.

[1] Organic liquid saturation distributions resulting from a simulated tetrachloroethene (PCE) spill were generated with alternative models of spatially varying aquifer properties for a statistically homogeneous, nonuniform sand aquifer. The distributions were analyzed to quantify DNAPL source zone characteristics and then incorporated as initial conditions for simulated PCE recovery using surfactant-enhanced aquifer remediation (SEAR). The predicted evolution of the spatial distribution of DNAPL saturations or source zone ''architectures'' and associated remediation efficiencies are strongly influenced by the spatial correlation of aquifer parameters and multiphase flow constitutive relationships. Model predictions suggest that removal of 60 to 99% of entrapped PCE can reduce dissolved contaminant concentration and mass flux under natural gradient conditions by approximately two orders of magnitude. Aqueous phase contaminant flux, however, does not vary consistently as a function of the percentage of DNAPL removed, and notable differences in flux evolution were observed for models incorporating correlated versus uncorrelated capillary entry pressure and permeability fields. Simulation results demonstrate that the application of alternative models of aquifer parameter spatial variability can influence predicted DNAPL infiltration, entrapment, and recovery, even for relatively homogeneous aquifers of the type investigated here. Results also demonstrate potential benefits, in the form of reduced mass flux, accruing from partial mass removal that may not be readily predicted from analyses relying on simplified conceptual models for DNAPL source zone architecture or aquifer flow fields.

A pilot-scale demonstration of surfactant-enhanced aquifer remediation (SEAR) was conducted to recover dense nonaqueous phase liquid (DNAPL) tetrachloroethene (PCE) from a sandy glacial outwash aquifer underlying a former dry cleaning facility at the Bachman Road site in Oscoda, MI. Part one of this two-part paper describes site characterization efforts and a comprehensive approach to SEAR test design, effectively integrating laboratory and modeling studies. Aquifer coring and drive point sampling suggested the presence of PCE-DNAPL in a zone beneath an occupied building. A narrow PCE plume emanating from the vicinity of this building discharges into Lake Huron. The shallow unconfined aquifer, characterized by relatively homogeneous fine-medium sand deposits, an underlying clay layer, and the absence of significant PCE transformation products, was judged suitable for the demonstration of SEAR. Tween 80 was selected for application based upon its favorable solubilization performance in batch and two-dimensional sand tank treatability studies, biodegradation potential, and regulatory acceptance. Three-dimensional flow and transport models were employed to develop a robust design for surfactant delivery and recovery. Physical and fiscal constraints led to an unusual hydraulic design, in which surfactant was flushed across the regional groundwater gradient, facilitating the delivery of concentrations of Tween 80 exceeding 1% (wt) throughout the treatment zone. The potential influence of small-scale heterogeneity on PCE-DNAPL distribution and SEAR performance was assessed through numerical simulations incorporating geostatistical permeability fields based upon available core data. For the examined conditions simulated PCE recoveries ranged from 94to 99%. The effluent treatment system design consisted of low-profile air strippers coupled with carbon adsorption to trap off-gas PCE and discharge of treated aqueous effluent to a local wastewater treatment plant. The systematic and comprehensive design methodology described herein may serve as a template for application at other DNAPL sites.

Dense nonaqueous phase liquid (DNAPL) source zones comprise persistent sources of groundwater contamination that are recalcitrant to complete remediation using conventional (e.g., pump and treat) or emerging (e.g., surfactant flushing) technologies. Increased attention to the assessment of the benefits of partial mass removal from such contaminant source zones has intensified efforts to model multiphase flow and transport behavior. This paper describes the simulated recovery of a tetrachloroethene (PCE) spill in a statistically homogeneous but nonuniform aquifer, incorporating nonuniformity in both nonaqueous phase liquid saturation and pore velocities. We developed a ganglia-to-pool metric to quantify DNAPL source-zone architecture, and explored the correlation of this metric with dissolved mass flux behavior in response to partial DNAPL mass removal. Dissolution of 20%-70% of PCE mass from models exhibiting low ganglia-to-pool ratios resulted in a larger predicted reduction of dissolved contaminant mass flux than models with high ganglia-to-pool ratios. Results of this study suggest that DNAPL source-zone characterization at field sites with homogeneous, nonuniform aquifers would benefit from inclusion of an estimate of the overall ganglia-to-pool ratio. Simulations demonstrate that flux reduction behavior depends on the source-zone architecture, which is not readily predictable using a priori assumptions about the spatial correlation of physical aquifer parameters. Model results further suggest that stochastic investigations of DNAPL source remediation at field sites should avoid reliance upon Leverett scaling of capillary entry pressures to permeability fields, which can artificially narrow the range of simulated behaviors.

No abstract