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

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

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.     

Effects of source zone heterogeneity on surfactant-enhanced NAPL dissolution and resulting remediation end-points

The effectiveness of removal of nonaqueous phase liquids (NAPLs) from the entrapment source zone of the subsurface has been limited by soil heterogeneity and the inability to locate all entrapped sources. The goal of this study was to demonstrate the uncertainty of degree of source removal associated with aquifer heterogeneity. In this demonstration, source zone NAPL removal using surfactant-enhanced dissolution was considered. Model components that simulate the processes of natural dissolution in aqueous phase and surfactant-enhanced dissolution were incorporated into an existing code of contaminant transport. The dissolution modules of the simulator used previously developed Gilland-Sherwood type phenomenological models of NAPL dissolution to estimate mass transfer coefficients that are upscaleable to multidimensional flow conditions found at field sites. The model was used to simulate the mass removal from 10 NAPL entrapment zone configurations based on previously conducted two-dimensional tank experiments. These entrapment zones represent the NAPL distribution in spatially correlated random fields of aquifer hydraulic conductivity. The numerical simulations representing two-dimensional conditions show that effectiveness of mass removal depends on the aquifer heterogeneity that controls the NAPL entrapment and delivery of the surfactant to the locations of entrapped NAPLs. Flow bypassing resulting from heterogeneity and the reduction of relative permeability due to NAPL entrapment reduces the delivery efficiency of the surfactant, thus prolonging the remediation time to achieve desired end-point NAPL saturations and downstream dissolved concentrations. In some extreme cases, the injected surfactant completely bypassed the NAPL source zones. It was also found that mass depletion rates for different NAPL source configurations vary significantly. The study shows that heterogeneity result in uncertainties in the mass removal and achievable end-points that are directly related to dissolved contaminant plume development downstream of the NAPL entrapment zone.     

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.     

Flow behavior and residual saturation formation of liquid carbon tetrachloride in unsaturated heterogeneous porous media

The formation of residual, discontinuous nonaqueous phase liquids (NAPLs) in the vadose zone is a process that is not well understood. To obtain data that can be used to study the development of a residual NAPL saturation in the vadose zone and to test current corresponding models, detailed transient experiments were conducted in intermediate-scale columns and flow cell. The column experiments were conducted to determine residual carbon tetrachloride (CCl(4)) saturations of two sands and to evaluate the effect of CCl(4) vapors on the water distribution. In the intermediate-scale flow cell experiment, a rectangular zone of the fine-grained sand was packed in an otherwise medium-grained matrix. A limited amount of CCl(4) was injected from a small source and allowed to redistribute until a pseudo steady state situation had developed. A dual-energy gamma radiation system was used to determine fluid saturations at numerous locations. The experiments clearly demonstrated the formation of residual CCl(4) saturations in both sands. Simulations with an established multifluid flow simulator show the shortcomings of current relative permeability-saturation-capillary pressure (k-S-P) models. The results indicate that nonspreading behavior of NAPLs should be implemented in simulators to account for the formation of residual saturations.     

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).     

Impact of nonaqueous phase liquid (NAPL) source zone architecture on mass removal mechanisms in strongly layered heterogeneous porous media during soil vapor extraction

An existing multiphase flow simulator was modified in order to determine the effects of four mechanisms on NAPL mass removal in a strongly layered heterogeneous vadose zone during soil vapor extraction (SVE): a) NAPL flow, b) diffusion and dispersion from low permeability zones, c) slow desorption from sediment grains, and d) rate-limited dissolution of trapped NAPL. The impacts of water and NAPL saturation distribution, NAPL-type (i.e., free, residual, or trapped) distribution, and spatial heterogeneity of the permeability field on these mechanisms were evaluated. Two different initial source zone architectures (one with and one without trapped NAPL) were considered and these architectures were used to evaluate seven different SVE scenarios. For all runs, slow diffusion from low permeability zones that gas flow bypassed was a dominant factor for diminished SVE effectiveness at later times. This effect was more significant at high water saturation due to the decrease of gas-phase relative permeability. Transverse dispersion contributed to fast NAPL mass removal from the low permeability layer in both source zone architectures, but longitudinal dispersion did not affect overall mass removal time. Both slow desorption from sediment grains and rate-limited mass transfer from trapped NAPL only marginally affected removal times. However, mass transfer from trapped NAPL did affect mass removal at later time, as well as the NAPL distribution. NAPL flow from low to high permeability zones contributed to faster mass removal from the low permeability layer, and this effect increased when water infiltration was eliminated. These simulations indicate that if trapped NAPL exists in heterogeneous porous media, mass transfer can be improved by delivering gas directly to zones with trapped NAPL and by lowering the water content, which increases the gas relative permeability and changes trapped NAPL to free NAPL.     

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 constitutive model for air-NAPL-water flow in the vadose zone accounting for immobile, non-occluded (residual) NAPL in strongly water-wet porous media

A hysteretic constitutive model describing relations among relative permeabilities, saturations, and pressures in fluid systems consisting of air, nonaqueous-phase liquid (NAPL), and water is modified to account for NAPL that is postulated to be immobile in small pores and pore wedges and as films or lenses on water surfaces. A direct outcome of the model is prediction of the NAPL saturation that remains in the vadose zone after long drainage periods (residual NAPL). Using the modified model, water and NAPL (free, entrapped by water, and residual) saturations can be predicted from the capillary pressures and the water and total-liquid saturation-path histories. Relations between relative permeabilities and saturations are modified to account for the residual NAPL by adjusting the limits of integration in the integral expression used for predicting the NAPL relative permeability. When all of the NAPL is either residual or entrapped (i.e., no free NAPL), then the NAPL relative permeability will be zero. We model residual NAPL using concepts similar to those used to model residual water. As an initial test of the constitutive model, we compare predictions to published measurements of residual NAPL. Furthermore, we present results using the modified constitutive theory for a scenario involving NAPL imbibition and drainage.     

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.     

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.     

Effect of water content on transient nonequilibrium NAPL-gas mass transfer during soil vapor extraction

The effect of water content on the volatilization of nonaqueous phase liquid (NAPL) in unsaturated soils was characterized by one-dimensional venting experiments conducted to evaluate the lumped mass transfer coefficient. An empirical correlation based upon the modified Sherwood number, Peclet number, and normalized mean grain size was used to estimate initial lumped mass transfer coefficients over a range of water content. The effects of water content on the soil vapor extraction SVE process have been investigated through experimentation and mathematical modeling. The experimental results indicated that a rate-limited NAPL-gas mass transfer occurred in water-wet soils. A severe mass transfer limitation was observed at 61.0% water saturation where the normalized effluent gas concentrations fell below 1.0 almost immediately, declined exponentially from the initiation of venting, and showed long tailing. This result was attributed to the reduction of interfacial area between the NAPL and mobile gas phases due to the increased water content. A transient mathematical model describing the change of the lumped mass transfer coefficient was used. Simulations showed that the nonequilibrium mass transfer process could be characterized by the exponent beta, a parameter which described the reduction of the specific area available for NAPL volatilization. The nonequilibrium mass transfer limitations were controlled by the soil mean grain size and pore gas velocity, were well described by beta values below 1.0 at low water saturation, and were well predicted with beta values greater than 1.0 at high water saturation.     

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.     

Performance assessment of NAPL remediation in heterogeneous alluvium

Over the last few years, more than 40 partitioning interwell tracer tests (PITTs) have been conducted at many different sites to measure nonaqueous phase liquid (NAPL) saturations in the subsurface. While the main goal of these PITTs was to estimate the NAPL volume in the subsurface, some were specifically conducted to assess the performance of remedial actions involving NAPL removal. In this paper, we present a quantitative approach to assess the performance of remedial actions to recover NAPL that can be used to assess any NAPL removal technology. It combines the use of PITTs (to estimate the NAPL volume in the swept pore volume between injection and extraction wells of a test area) with the use of several cores to determine the vertical NAPL distribution in the subsurface. We illustrate the effectiveness of such an approach by assessing the performance of a surfactant/foam flood conducted at Hill Air Force Base, UT, to remove a TCE-rich NAPL from alluvium with permeability contrasts as high as one order of magnitude. In addition, we compare the NAPL volumes determined by the PITTs with volumes estimated through geostatistical interpolation of aquifer sediment core data collected with a vertical frequency of 5-10 cm and a lateral borehole spacing of 0.15 m. We demonstrate the use of several innovations including the explicit estimation of not only the errors associated with NAPL volumes and saturations derived from PITTs but also the heterogeneity of the aquifer sediments based upon permeability estimates. Most importantly, we demonstrate the reliability of the     

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

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

A proposed model to include a residual NAPL saturation in a hysteretic capillary pressure-saturation relationship

A residual non-aqueous phase liquid (NAPL) present in the vadose zone can act as a contaminant source for many years as the compounds of concern partition to infiltrating groundwater and air contained in the soil voids. Current pressure-saturation-relative permeability relationships do not include a residual NAPL saturation term in their formulation. This paper presents the results of series of two- and three-phase pressure cell experiments conducted to evaluate the residual NAPL saturation and its impact on the pressure-saturation relationship. A model was proposed to incorporate a residual NAPL saturation term into an existing hysteretic three-phase parametric model developed by Parker and Lenhard [Water Resour. Res. 23(12) (1987) 2187], Lenhard and Parker [Water Resour. Res. 23(12) (1987) 2197] and Lenhard [J. Contam. Hydrol. 9 (1992) 243]. The experimental results indicated that the magnitude of the residual NAPL saturation was a function of the maximum total liquid saturation reached and the water saturation. The proposed model to incorporate a residual NAPL saturation term is similar in form to the entrapment model proposed by Parker and Lenhard, which was based on an expression presented by Land [Soc. Pet. Eng. J. (June 1968) 149].     

Three-dimensional visualization and quantification of non-aqueous phase liquid volumes in natural porous media using a medical X-ray Computed Tomography scanner

This study demonstrates the capabilities of a typical medical X-ray Computed Tomography (CT) scanner to non-destructively quantify non-aqueous phase liquid (NAPL) volumes, saturation levels, and three-dimensional spatial distributions in packed soil columns. Columns packed with homogeneous sand, heterogeneous sand, or natural soil, were saturated with water and injected with known quantities of gasoline or tetrachloroethene and scanned. A methodology based on image subtraction was implemented for computing soil porosity and NAPL volumes in each 0.35 mm x 0.35 mm x 1 mm voxel of the columns. Elimination of sample positioning errors and instrument drift artifacts was essential for obtaining reliable estimates of above parameters. The CT data-derived total NAPL volume was in agreement with the measured NAPL volumes injected into the columns. CT data-derived NAPL volume is subject to a 2.6% error for PCE and a 15.5% error for gasoline, at average NAPL saturations as low as 5%, and is mainly due to instrument noise. Non-uniform distributions of NAPL due to preferential flow, and accumulation of NAPL above finer-grained layers could be observed from the data on 3-D distributions of NAPL volume fractions.     

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.     

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.     

Predictive modelling of NAPL injection tests in variable aperture spatially correlated fractures.

In preparation for a field experiment where a NAPL will be injected into a fractured sandstone aquifer, a 2D invasion percolation model of DNAPL migration in a single horizontal fracture with varying aperture has been developed. This simulation investigated the effect of spatially correlated fracture aperture fields on pressure-saturation relationships as a function of variable aperture mean, standard deviation, and spatial correlation statistics under hydrostatic conditions. Results from spatially correlated variable aperture fields can be significantly different from random fields. Longer ranges decreased entry pressures and higher standard deviations decreased nonwetting phase saturations. Mean aperture is the major control on displacement pressure, entry pressure and the form of the pressure-saturation curve. Simulation results using statistical parameters for a variable aperture natural sandstone fracture as described by Yeo et al. [International Journal of Rock Mechanics and Mining Sciences 35 (1998) 1051] closely resemble a Brooks-Corey pressure-saturation function, and exhibit a distinct entry pressure followed by a rapid increase in nonwetting phase saturation. Graphical estimates of entry pressure provide a good approximation of the critical aperture controlling the formation of a continuous nonwetting phase pathway in a variable aperture fracture. Consequently, we show that multiphase flow concepts developed for porous media can successfully be applied to variable aperture fractures. Entry pressure correlates well to the mean aperture in these simulations when using a Gaussian distribution of fracture aperture. Interpreted aperture distributions from NAPL injection experiments do not fit the true distribution well at low nonwetting phase saturations, but do at higher saturations above the entry pressure. Consequently, true, mechanical aperture variation within a fracture plane cannot be determined from NAPL injection experiments either in the field or laboratory.