Using polymer solutions to enhance recovery of mobile coal tar and creosote DNAPLs

Direct pumping and enhanced recovery of coal tar and creosote dense, non-aqueous phase liquids (DNAPLs) from the subsurface have had mixed results because these DNAPLs are viscous fluids that can potentially alter aquifer wettability. To improve the inefficiencies associated with waterflooding, the research presented here considered the use of a polymer solution that can be added to the injected flood solution to increase the viscosity and decrease the velocity of the flooding solution. Results from one-dimensional, vertically oriented laboratory column experiments that evaluate the recovery of coal-derived DNAPL with both water and polymer flooding solutions are presented. The final DNAPL saturation remaining in the column was assessed in water and oil-wet systems for three viscous DNAPLs. Adding polymer to increase the aqueous solution viscosity did not have a significant impact in water-wet systems. A final DNAPL saturation of approximately 19% was achieved for both water and polymer floods. In contrast, the addition of polymer significantly improved recovery in oil-wet systems. The final saturation was over 40% in oil-wet systems after waterflooding, but approximately 19% with a polymer flushing solution. Although the final saturation produced with polymer flooding was similar between the oil- and water-wet systems, differences in the relative permeability and distribution of DNAPL in the porous matrix caused the DNAPL recovery to be much slower in the oil-wet system.     

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.     

Removal of PCB-DNAPL from a rough-walled fracture using alcohol/polymer flooding

Phase behaviour experiments employing PCB (Aroclor 1242)/alcohol/water systems were conducted with ethanol (EtOH) and n-propanol (nPA). Both exhibited an affinity for the aqueous phase within the entire two-phase region. As much as 88% by volume (88% vol.) EtOH and 80% vol. nPA were necessary to achieve full miscibility of the PCB in the aqueous phase. DNAPL-water interfacial tension (IFT) was reduced from 38.9 dyn/cm to 4.7 dyn/cm and 2.4 dyn/cm with 80% vol. EtOH and 76% vol. nPA. The addition of alcohol brought about 41% and 54% reductions in DNAPL viscosity at maximal concentrations of EtOH and nPA. Density of the PCB-DNAPL was relatively unaffected by the presence of alcohol. A series of seven experiments were conducted where successive slugs of nPA and xanthan gum polymer solutions were injected into a fractured shale sample. A 30% vol. nPA solution injected under a hydraulic gradient of 0.36 allowed enhanced PCB removal primarily through reduction of IFT and resulted in 72% DNAPL recovery. Several pore volumes of alcohol solution were necessary to displace all the potentially mobile non-wetting phase since the high-viscosity DNAPL was mobilized at a lower flow rate than the overall fluid velocity, illustrating non-piston displacement. The injection of a 95% vol. nPA alcohol solution, theoretically at a sufficient concentration to produce fully miscible displacement of the residual DNAPL at equilibrium, resulted in non-equilibrium partitioning of the PCB into the flushing solution, likely due to the high fluid velocities in the fracture. The injection of 200 pore volumes of 95% vol. nPA solution resulted in 94% DNAPL recovery. Alcohol floods operated below the miscibility envelope appear to be a valuable source zone remedial alternative where the objective is to reduce DNAPL mobility to zero, but it should be noted that DNAPL mobility is increased during the application of the technology and steps may need to be taken to prevent unwanted vertical mobilization.     

Dip-angle influence on areal DNAPL recovery by co-solvent flooding with and without pre-flooding

A two-dimensional (2D) laboratory model was used to study effects of gravity on areal recovery of a representative dense non-aqueous phase liquid (DNAPL) contaminant by an alcohol pre-flood and co-solvent flood in dipping aquifers. Recent studies have demonstrated that injection of alcohol and co-solvent solutions can be used to reduce in-situ the density of DNAPL globules and displace the contaminant from the source zone. However, contact with aqueous alcohol reduces interfacial tension and causes DNAPL swelling, thus facilitating risk of uncontrolled downward DNAPL migration. The 2D laboratory model was operated with constant background gradient flow and a DNAPL spill was simulated using tetrachloroethene (PCE). The spill was dispersed to a trapped, immobile PCE saturation by a water flood. Areal PCE recovery was studied using a double-triangle well pattern to simulate a remediation scheme consisting of an alcohol pre-flood using aqueous isobutanol ( approximately 10% vol.) followed by a co-solvent flood using a solution of ethylene glycol (65%) and 1-propanol (35%). Experiments were conducted with the 2D model oriented in the horizontal plane and compared to experiments at the 15 degrees and 30 degrees dip-angle orientations. Injection was applied either in the downward or upward direction of flow. Experimental results were compared to theoretical predictions for flood front stability and used to evaluate effects of gravity on areal PCE recovery. Sensitivity experiments were performed to evaluate effects of the alcohol pre-flood on PCE areal recovery. For experiments conducted with the alcohol pre-flood and the 2D model oriented in the horizontal plane, results indicate that 89-93% of source zone PCE was recovered. With injection oriented downward, results indicate that areal PCE recovery was 70-77% for a 15 degrees dip angle and 57-59% for a 30 degrees dip angle. With injection oriented upward, results indicate that areal PCE recovery was 57-60% at the 30 degrees dip angle, which was similar to PCE recovery for injection in the downward flow direction. Lower areal PCE recovery at greater dip angles in either direction of flow was attributed to DNAPL swelling and migration, flood front instabilities and bypassing of the displaced fluid past the extraction wells during the alcohol pre-flood. Additional results demonstrate that the use of an alcohol pre-flood can be beneficial in improving DNAPL recovery in the horizontal orientation, but pre-flooding may reduce areal recovery efficiency in dip-angle orientations. This study also demonstrates the use of theoretical perturbation (fingering) analysis in predicting NAPL recovery efficiency for flooding processes in remediating aquifers with dip angles.     

Effects of initial saturation on properties modification and displacement of tetrachloroethene with aqueous isobutanol

Packed column experiments were conducted to study effects of initial saturation of tetrachloroethene (PCE) in the range of 1.0-14% pore volume (PV) on mobilization and downward migration of the non-aqueous phase liquid (NAPL) product upon contact with aqueous isobutanol ( approximately 10 vol.%). This study focused on the consequences of swelling beyond residual saturation. Columns were packed with mixtures of neat PCE, water and glass beads and waterflooded to establish a desired homogeneous residual saturation, and then flooded with aqueous isobutanol under controlled hydraulic conditions. Results showed a critical saturation of approximately 8% PV for these packed column experimental conditions. At low initial PCE saturations (<8% PV), experimental results showed reduced risk of NAPL-product migration upon contact with aqueous isobutanol. At higher initial PCE saturations (>8% PV), results showed NAPL-product mobilization and downward migration which was attributed to interfacial tension (IFT) reduction, swelling of the NAPL-product, and reduced density modification. Packed column results were compared with good agreement to theoretical predictions of NAPL-product mobilization using the total trapping number, N(T). In addition to the packed column study, preliminary batch experiments were conducted to study the effects of PCE volumetric fraction in the range of 0.5-20% on density, viscosity, and IFT modification as a function of time following contact with aqueous isobutanol ( approximately 10 vol.%). Modified NAPL-product fluid properties approached equilibrium within approximately 2 h of contact for density and viscosity. IFT reduction occurred immediately as expected. Measured fluid properties were compared with good agreement to theoretical equilibrium predictions based on UNIQUAC. Overall, this study demonstrates the importance of initial DNAPL saturation, and the associated risk of downward NAPL-product migration, in applying alcohol flooding for remediation of DNAPL contaminated ground water sites.     

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.     

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.     

Removal of DNAPL contamination from the saturated zone by the combined effect of vertical upward flushing and density reduction

A common aspect of innovative remediation techniques is that they tend to reduce the interfacial tension between the aqueous and non-aqueous phase liquids, resulting in mobilization of the organic contaminant. This complicates the remediation of aquifers, contaminated with Dense Non-Aqueous Phase Liquids (DNAPLs), as they are likely to migrate downwards, deeper into the aquifer and into finer layers. A possible solution is the use of swelling alcohols, which tend to reduce the density difference between the aqueous phase and the DNAPL. To avoid premature mobilization upon the initial contact between the DNAPL and the alcohol, several researchers have proposed the use of vertical upward flow of the alcohol. In this paper, we present an equation, which describes the upward mobilization of both continuous and discontinuous DNAPLs and so the important parameters governing the upward controlled mobilization of the DNAPL. The need and required magnitude of this specific discharge was investigated by conducting four column experiments in which the initial density of the DNAPL and the permeability was varied. It was shown that the required flow velocities increase with the permeability of the porous medium and the initial density difference between the aqueous phase and the DNAPL. Whenever the specific discharge falls below the critical value, the DNAPL moves downward. A second set of column experiments looked at the impact of permeability of porous medium on the solubilization and mobilization of DNAPL during alcohol flooding. Columns, packed with coarse or fine sand, containing a residual trichloroethylene (TCE) or perchloroethylene (PCE) saturation were flushed with the alcohol mixture at a fixed specific discharge rate. The induced pressure gradients in the aqueous phase, which were higher in the fine sand, resulted for this porous medium in extensive mobilization of the DNAPL against the direction of the buoyancy force. The density of the first NAPL coming out of the top of the fine sand was close to that of the pure DNAPL. In the coarser sand, the pressure gradients were sufficient to prevent downward migration of the DNAPL, but upward mobilization was minimal. The predominant removal mechanism in this case was the much slower solubilization.     

Laboratory investigation of flux reduction from dense non-aqueous phase liquid (DNAPL) partial source zone remediation by enhanced dissolution

This study investigated the benefits of partial removal of dense nonaqueous phase liquid (DNAPL) source zones using enhanced dissolution in eight laboratory scale experiments. The benefits were assessed by characterizing the relationship between reductions in DNAPL mass and the corresponding reduction in contaminant mass flux. Four flushing agents were evaluated in eight controlled laboratory experiments to examine the effects of displacement fluid property contrasts and associated override and underride on contaminant flux reduction (R(j)) vs. mass reduction (R(m)) relationships (R(j)(R(m))): 1) 50% ethanol/50% water (less dense than water), 2) 40% ethyl-lactate/60% water (more dense than water), 3) 18% ethanol/26% ethyl-lactate/56% water (neutrally buoyant), and 4) 2% Tween-80 surfactant (also neutrally buoyant). For each DNAPL architecture evaluated, replicate experiments were conducted where source zone dissolution was conducted with a single flushing event to remove most of the DNAPL from the system, and with multiple shorter-duration floods to determine the path of the R(j)(R(m)) relationship. All of the single-flushing experiments exhibited similar R(j)(R(m)) relationships indicating that override and underride effects associated with cosolvents did not significantly affect the remediation performance of the agents. The R(j)(R(m)) relationship of the multiple injection experiments for the cosolvents with a density contrast with water tended to be less desirable in the sense that there was less R(j) for a given R(m). UTCHEM simulations supported the observations from the laboratory experiments and demonstrated the capability of this model to predict R(j)(R(m)) relationships for non-uniformly distributed NAPL sources.     

Modeling field-scale cosolvent flooding for DNAPL source zone remediation

A three-dimensional, compositional, multiphase flow simulator was used to model a field-scale test of DNAPL removal by cosolvent flooding. The DNAPL at this site was tetrachloroethylene (PCE), and the flooding solution was an ethanol/water mixture, with up to 95% ethanol. The numerical model, UTCHEM accounts for the equilibrium phase behavior and multiphase flow of a ternary ethanol-PCE-water system. Simulations of enhanced cosolvent flooding using a kinetic interphase mass transfer approach show that when a very high concentration of alcohol is injected, the DNAPL/water/alcohol mixture forms a single phase and local mass transfer limitations become irrelevant. The field simulations were carried out in three steps. At the first level, a simple uncalibrated layered model is developed. This model is capable of roughly reproducing the production well concentrations of alcohol, but not of PCE. A more refined (but uncalibrated) permeability model is able to accurately simulate the breakthrough concentrations of injected alcohol from the production wells, but is unable to accurately predict the PCE removal. The final model uses a calibration of the initial PCE distribution to get good matches with the PCE effluent curves from the extraction wells. It is evident that the effectiveness of DNAPL source zone remediation is mainly affected by characteristics of the spatial heterogeneity of porous media and the variable (and unknown) DNAPL distribution. The inherent uncertainty in the DNAPL distribution at real field sites means that some form of calibration of the initial contaminant distribution will almost always be required to match contaminant effluent breakthrough curves.     

Displacement and sweep efficiencies in a DNAPL recovery test using micellar and polymer solutions injected in a five-spot pattern

Soil washing with micellar solutions is a promising alternative for the remediation of DNAPL source zones. As with any flushing technology, the success of soil washing with micellar solutions depends in a very large part on the ability of the solution to contact the contaminant (sweep efficiency) and then on the efficiency of contaminant removal once this contact is made (displacement efficiency). We report here on a field test where a micellar solution was used to recover a DNAPL in an open five-spot pattern in which polymer solutions were also injected before and after the washing solution to improve sweep efficiency. The washing solution formulation was optimised in the laboratory prior to the test to obtain good dissolution capacity. For a high-concentration and low-volume soil flushing remediation test such as the one performed (0.8 pore volumes of actual washing solution injected), slug sizing of the washing solution is critical. It was evaluated by an analytical solution. In a five-spot pattern, the displacement efficiency of the washing solution was observed to vary in the porous medium as a function of the radial distance from the injection well because: (1) the volume of the washing solution flowing through a section of the test cell changes (maximum close to the injection well and minimal at the pumping wells); (2) the in situ velocity changes (maximum at the wells and minimum between the wells) and; (3) the contact time of the washing solution with the NAPL changes as a function of the distance from the injection well. The relative importance of the recovery mechanisms, mobilisation and dissolution, was also observed to vary in the test cell. The reduced velocity increased the contact time of the washing solution with the DNAPL enhancing its dissolution, but the decrease of the capillary number caused less mobilisation. The washing process is much more extensive around the injection well. The use of an injection-pumping pattern allowing a complete sweep of the remediated area is essential. Following a comprehensive characterisation, modeling is an efficient tool to design the injection-pumping scheme and to optimise injection and pumping rates providing the best areal sweep. The vertical sweep can be controlled by using a polymer solution (Xanthan gum). The polymer solution also has a positive effect on front stability between the solutions injected. The injection rate of the polymer solution that follows the washing solution must be kept minimal initially to prevent dilution of the washing solution by fingering.     

DNAPL source depletion: linking architecture and flux response

The relationship between dense non-aqueous phase liquid (DNAPL) mass reduction and contaminant mass flux was investigated experimentally in four model source zones. The flow cell design for the experiments featured a segmented extraction well that allowed for analysis of spatially resolved flux information. This flux information was coupled with image analysis of the NAPL spatial distribution to investigate the relationship between flux and the up-gradient NAPL architecture. Results indicate that in the systems studied, the relationship between DNAPL mass reduction and contaminant mass flux was primarily controlled by the NAPL architecture. A specific definition of NAPL architecture was employed where the source zone is resolved into a collection of streamtubes with spatial variability in NAPL saturation along each streamtube integrated and transformed into an effective NAPL content for each streamtube. The distribution of NAPL contents among the streamtubes (NAPL architecture) controlled dissolution dynamics. Two simplified models, a streamtube model and an effective Damkohler number model, were investigated for their ability to simulate dissolution dynamics.     

Equilibrium nonaqueous phase liquid pool geometry in coarse soils with discrete textural interfaces

This paper presents a model for the geometry of nonaqueous phase liquid (NAPL) pools and mounds in homogeneous soils and soils with discrete textural interfaces. It is shown that the concepts of capillary pressure-saturation curve hysteresis and entry pressures are integral to the complete conceptualization of pool and mound geometry. Unless hysteresis is included in the analysis, light NAPL (LNAPL) in homogeneous soils cannot exist in pools at all, and dense NAPL (DNAPL) will not mound on horizontal textural interfaces unless lateral confining boundaries are present. The proposed model also implies that remobilization of DNAPL pools will occur at lower hydraulic gradients than those predicted with previous models. Comparing predicted and experimental DNAPL and LNAPL pool thicknesses and the location of an LNAPL lens with respect to the top of the capillary fringe validate the model.     

Mobilization and entry of DNAPL pools into finer sand media by cosolvents: two-dimensional chamber studies

Two-dimensional chamber studies were conducted to determine qualitative and quantitative performance of cosolvents targeted at pooled dense non-aqueous phase liquid (DNAPL) (perchlorethylene, PCE) residing above a fine-grain capillary barrier. Downward mobilization of DNAPL, up gradient along an overriding cosolvent front, was observed. This produced significant pooling above a fine-grain layer that in some cases lead to entry into the capillary barrier beneath. Entry pressure calculations using physical and hydrogeologic parameters provided an excellent prediction of breakthrough of DNAPL into the capillary barrier. Calculations predict approximately 0.5 m of DNAPL would be necessary to enter a Beit Netofa clay, under extreme cosolvent flooding conditions (100% ethanol). Gradient injection of cosolvent did not appear to provide any benefit suggesting a rapid decrease in interfacial tension (IFT) compared to the rate of DNAPL solubilization. Use of a partitioning alcohol (tertiary butyl alcohol, TBA) resulted in DNAPL swelling and reduced entry into the capillary barrier. However, the trapping of flushing solution, containing PCE, could potentially lead to longer remediation times.     

Co-injection of air and steam for the prevention of the downward migration of DNAPLs during steam enhanced extraction: an experimental evaluation of optimum injection ratio predictions

When steam is injected into soil containing a dense volatile non-aqueous phase liquid contaminant, the DNAPL vaporized within the heated soil region condenses and accumulates ahead of the steam condensation front. If enough DNAPL accumulates, gravitational forces can overcome trapping forces allowing the liquid contaminant to flow downward. By injecting air with steam, a portion of the DNAPL vapor remains suspended in equilibrium with the air, decreasing liquid contaminant accumulation ahead of the steam condensation front, and thus reducing the possibility of downward migration. In a previous work, a theoretical model was developed to predict the optimum injection ratio of air to steam that would eliminate accumulation of DNAPL ahead of the temperature front and thus minimize the potential for downward migration. In this work, the theoretical model is summarized, and an experiment is presented in order to evaluate the optimum injection ratio prediction. In the experiment, a two-dimensional water saturated sand pack is contaminated with a known mass of TCE (DNAPL). The system is then remediated by co-injecting air and steam at the predicted optimum injection ratio, calculated based on the average contaminant soil concentration in the sand pack. Results for the co-injection of air and steam are compared to results for the injection of pure steam or pure air. Injection at the predicted optimum injection ratio for a volumetric average NAPL saturation, reduced accumulation of the contaminant ahead of the condensation front by over 90%, as compared to steam injection alone. This indicates that the optimum injection ratio prediction is a valuable tool for limiting the spreading of DNAPL during steam-enhanced extraction. Injection at the optimum injection ratio resulted in earlier recovery of contaminant than for steam injection alone. Co-injection of steam and air is also shown to result in much higher recovery rates than air injection alone.     

DNAPL remediation with in situ chemical oxidation using potassium permanganate. II. Increasing removal efficiency by dissolving Mn oxide precipitates

In situ chemical oxidation (ISCO) schemes using MnO4- have been effective in destroying chlorinated organic solvents dissolved in ground water. Laboratory experiments and field pilot tests reveal that the precipitation of Mn oxide, one of the reaction products, causes a reduction of permeability, which can lead to flow bypassing and inefficiency of the scheme. Without a solution to this problem of plugging, it is difficult to remove DNAPL from the subsurface completely. In a companion paper, we showed with batch experiments that Mn oxide can be dissolved rapidly with certain organic acids. This study utilizes 2-D flow-tank experiments to examine the possibility of nearly complete DNAPL removal by ISCO with MnO4-, when organic acids are used to remove Mn oxide. The experiments were conducted in a small 2-D glass flow tank containing a lenticular silica-sand medium. Blue-dyed trichloroethylene (TCE) provided residual, the perched and pooled DNAPL. KMnO4 at 200 mg/l was flushed through the DNAPL horizontally. Once plugging reduced permeability and prevented further delivery of the oxidant, citric or oxalic acids were pumped into the flow tank to dissolve the Mn oxide precipitates. Organic ligands removed the Mn oxide precipitates relatively quickly, and permitted another cycle of MnO4- flooding. Cycles of MnO4-/acid flooding continued until all of the visible DNAPL was removed. The experiments were monitored with chemical analysis and visualization. A mass-balance calculation indicated that by the end of the experiments, all the DNAPL was removed. The results show also how heterogeneity adds complexity to initial redistribution of DNAPL, and to the efficiency of the chemical flooding.     

Design of aquifer remediation systems: (1) describing hydraulic structure and NAPL architecture using tracers

Aquifer heterogeneity (structure) and NAPL distribution (architecture) are described based on tracer data. An inverse modelling approach that estimates the hydraulic structure and NAPL architecture based on a Lagrangian stochastic model where the hydraulic structure is described by one or more populations of lognormally distributed travel times and the NAPL architecture is selected from eight possible assumed distributions. Optimization of the model parameters for each tested realization is based on the minimization of the sum of the square residuals between the log of measured tracer data and model predictions for the same temporal observation. For a given NAPL architecture the error is reduced with each added population. Model selection was based on a fitness which penalized models for increasing complexity. The technique is demonstrated under a range of hydrologic and contaminant settings using data from three small field-scale tracer tests: the first implementation at an LNAPL site using a line-drive flow pattern, the second at a DNAPL site with an inverted five-spot flow pattern, and the third at the same DNAPL site using a vertical circulation flow pattern. The Lagrangian model was capable of accurately duplicating experimentally derived tracer breakthrough curves, with a correlation coefficient of 0.97 or better. Furthermore, the model estimate of the NAPL volume is similar to the estimates based on moment analysis of field data.     

Investigation of the Contamination of Fractured Formations by Non-Newtonian Oil Pollutants

In many practical applications, non-aqueous phase liquid (NAPL) pollutants exhibiting a clearly non-Newtonian rheological behavior (e.g. crude oil, suspensions of engine oils, asphalt, creosote, etc.) may migrate through fractured formations and contaminate aquifers. The present work is the first step toward the development of non-Darcian models concerning the non-linear NAPL flow in single fractures, and detemiination of the coupled effects of non-Newtonian NAPL rheology and flow rate on the transient immiscible displacement of an aqueous phase by a NAPL. Initially, a protocol is developed for the preparation and rheological characterization of synthetic non-Newtonian NAPLs, which are based on waxy oils. Then, an artificial transparent glass-etched single fracture of controlled morphology is fabricated and used for the measurement of the non-linear pressure gradient--superficial velocity relationship for the flow of NAPL of varying rheology. Pore network simulations and effective medium approximation (EMA) are used for the interpretation of the experimental results and derivation of an analytic non-Darcian one-phase flow model. Visualization experiments of the immiscible displacement of an aqueous phase by Newtonian and non-Newtonian NAPLs are performed on the artificial fracture under controlled values of the viscosity ratio and capillary number (ratio of viscous to capillary forces). Comparative study of the Newtonian and non-Newtonian NAPL flow patterns allows us to evaluate the interactive effects of NAPL rheology, flow rates and fracture morphology on the spatial and temporal distribution of such liquid pollutants within single fractures     

Investigation of the Contamination of Fractured Formations by Non-Newtonian Oil Pollutants

In many practical applications, non-aqueous phase liquid (NAPL) pollutants exhibiting a clearly non-Newtonian rheological behavior (e.g. crude oil, suspensions of engine oils, asphalt, creosote, etc.) may migrate through fractured formations and contaminate aquifers. The present work is the first step toward the development of non-Darcian models concerning the non-linear NAPL flow in single fractures, and determination of the coupled effects of non-Newtonian NAPL rheology and flow rate on the transient immiscible displacement of an aqueous phase by a NAPL. Initially, a protocol is developed for the preparation and rheological characterization of synthetic non-Newtonian NAPLs, which are based on waxy oils. Then, an artificial transparent glass-etched single fracture of controlled morphology is fabricated and used for the measurement of the non-linear pressure gradient—superficial velocity relationship for the flow of NAPL of varying rheology. Pore network simulations and effective medium approximation (EMA) are used for the interpretation of the experimental results and derivation of an analytic non-Darcian one-phase flow model. Visualization experiments of the immiscible displacement of an aqueous phase by Newtonian and non-Newtonian NAPLs are performed on the artificial fracture under controlled values of the viscosity ratio and capillary number (ratio of viscous to capillary forces). Comparative study of the Newtonian and non-Newtonian NAPL flow patterns allows us to evaluate the interactive effects of NAPL rheology, flow rates and fracture morphology on the spatial and temporal distribution of such liquid pollutants within single fractures.     

Relative velocities of DNAPL and aqueous phase plume migration

Numerical simulation is used to examine the relative velocities of DNAPL and aqueous phase plumes in sandy aquifers where lateral spreading of DNAPL has occurred at the base of the aquifer. The scenario being modeled is one where a permeable aquifer is underlain by a sloping aquitard, which results in lateral migration of the DNAPL down the slope, in addition to lateral migration of an aqueous phase plume subject to a specified hydraulic gradient. A sensitivity analysis is presented to the impacts of both DNAPL properties and geologic properties. The most important chemical properties governing the relative velocities of the DNAPL and the shallow aqueous phase plume are the DNAPL viscosity and the aqueous component soil-water partition coefficient (Kd). The dip of the underlying aquitard was found to be relatively unimportant, at least for the range of values studied. The scenario under consideration can be important in conceptual model development and remedial design, as in certain cases DNAPL could be migrating in areas without the evidence of a well-developed aqueous phase plume. The implication of this work is that the absence of a shallow aqueous phase plume directly downgradient of a DNAPL source zone does not rule out the possibility of deep occurrences of DNAPL beyond the shallow monitoring well network. A further finding of this study is that the occurrence of a highly sorbing compound in groundwater at virtually any concentration may indicate the immediate upgradient presence of residual or pooled DNAPL.