Refinement of Bouwer‐Rice Baildown Test Analysis

The recoverability of light nonaqueous phase liquids (LNAPL) in the subsurface can be estimated using LNAPL transmissivity. LNAPL transmissivity is analogous to aquifer transmissivity in that it represents the volume of LNAPL that flows through a unit width of a porous medium for a unit gradient in a unit time. Methods for estimating LNAPL transmissivity from baildown test data have been modified from the Bouwer and Rice (1976) slug test method by Lundy and Zimmerman (1996) and Huntley (2000). The primary assumptions when estimating LNAPL transmissivity with the Bouwer‐Rice method include, a quasi‐steady‐state model for recharge to the well (the model assumes steady‐state radial flow to the well with rate dependent well drawdown and no storage effects) and that the ratio of change in LNAPL drawdown to change in LNAPL thickness at the well (ds/db) is constant. This ratio will be referred to as the j‐ratio. Rather than having to meet boundary conditions for a predetermined j‐ratio value, each baildown test provides the data to estimate the j‐ratio value that is unique to that data set. This calculation methodology in turn results in the Bouwer‐Rice method being applicable to a wider range of baildown tests where it is not required for the potentiometric surface or LNAPL/water interface to remain constant.     

Identification and Assessment of Confined and Perched LNAPL Conditions

Although confined and perched light nonaqueous phase liquids (LNAPLs) have previously been recognized, the majority of technical LNAPL literature focuses on unconfined LNAPL. Little information exists regarding the appropriate use of LNAPL distribution and transmissivity data to distinguish between confined, perched, and unconfined LNAPL hydrogeological scenarios. This paper describes three case histories that illustrate how the observed behavior of LNAPL can be used to identify the hydrogeologic condition of LNAPL at a given site and improved methods for calculating LNAPL drawdown based on these hydrogeologic conditions. The assessment methodology uses routinely available data such as fluid gauging, boring lo, laser‐induced fluorescence, visual observations of soil cores, and LNAPL baildown testing. Identification of the correct LNAPL hydrogeologic condition results in more accurate LNAPL conceptual site models, improved estimates of LNAPL recovery rates and volumes, more appropriate technology applications, and improved accuracy of LNAPL remediation metrics such as LNAPL transmissivity.     

Evaluating an Analytical Model to Predict Subsurface LNAPL Distributions and Transmissivity from Current and Historic Fluid Levels in Groundwater Wells: Comparing Results to Numerical Simulations

A recent analytical model predicts free, entrapped, and residual LNAPL saturations and the LNAPL transmissivity in the subsurface from current and historic fluid levels in groundwater wells. As such, the model accounts for effects of fluid level fluctuations in a well. The model was developed to predict LNAPL specific volumes and transmissivities from current fluid level measurements in wells and either recorded historic fluid level fluctuations in wells or estimates. An assumption is made in the model that the predictions are not dependent on whether the historic highest or lowest fluid level elevations in a well occur first. To test the assumption, we conduct two simulations with a modified multiphase flow numerical code TMVOC that incorporates relative permeability‐saturation‐capillary head relations employed in the model. In one simulation, the initial condition is for fluid levels in a well at the historic highest elevations. In the other simulation, the initial condition is for fluid levels in a well at the historic lowest elevations. We change the boundary conditions so both historical conditions occur followed by generating the current condition. Results from the numerical simulations are compared to model predictions and show the assumption in the analytical model is reasonable. The analytical model can be used to develop/refine conceptual site models and for assessing potential LNAPL recovery endpoints, especially on sites with fluctuating fluid levels in wells.     

LNAPL Recovery Endpoints: Lessons Learnt Through Modeling,Experiments, and Field Trials

It is important to estimate what light nonaqueous phase liquid (LNAPL) recovery can be practicably achieved from subsurface environments. Over the last decade, research to address this included a broad field program, laboratory measurements and experimentation, and modeling approaches. Here, we consolidate key findings from the research in the context of current literature and understanding, with a focus on a well-validated, multiphase multicomponent modeling approach to achieve estimates of reasonable endpoints for LNAPL recovery. Simple analytical models can provide approximate saturation distributions and estimates of LNAPL recoverability via transmissivity approximation, but are insufficient to predict LNAPL saturation- and composition-based recovery endpoints for various recovery technologies. This is because they cannot account for multiphase, multicomponent fate and transport and key processes such as hysteresis. Recent advances to improve estimates of the fraction of recoverable LNAPL and its transmissivity are summarized. These advances include further development and application of a well-validated model to characterize active LNAPL recovery endpoints. We present key factors that affect the determination of LNAPL recovery endpoints, and outline how recovery endpoints are affected by natural source zone depletion (NSZD—currently gaining acceptance as a LNAPL remediation option). Major factors include geo-physical characteristics of the formation, magnitude of an LNAPL release and partitioning properties of the key LNAPL constituents of concern. Based on the capabilities of the validated model, the paper also provides a basis to optimize LNAPL recovery efforts.     

Magnetic Resonance Imaging of Nonaqueous Phase Liquid Distribution in Unconsolidated Soils Using Discriminatory Freezing

Organic contaminants present as nonaqueous phase liquids (NAPLs) in the subsurface often pose a long-term risk to human health and the environment. Investigating the distribution of NAPLs in porous media remains a major challenge in risk assessment and management of contaminated sites. Conventional soil coring and monitoring wells have been widely used over past decades as the primary means of subsurface investigation to determine NAPL extent. Known limitations of conventional approaches have led us to explore an alternative or a complementary technique to provide high-quality information of NAPL source zone architecture. This work advances an imaging tool for a variety of organic NAPL contaminants in unconsolidated soils through magnetic resonance imaging (MRI) of frozen cores. Using trichloroethylene (TCE) and o-xylene as model species, we illustrate that discriminatory freezing of water, while keeping the NAPL in a liquid state, enables high-resolution qualitative delineation of NAPL distribution within porous media. This novel approach may help improve site conceptual models and consequentially lead to highly tailored, more efficient remedial measures.     

A practical model for mobile, residual, and entrapped NAPL in water-wet porous media

Flow of nonvolatile nonaqueous phase liquid (NAPL) and aqueous phases that account for mobile, entrapped, and residual NAPL in variably saturated water-wet porous media is modeled and compared against results from detailed laboratory experiments. Residual saturation formation in the vadose zone is a process that is often ignored in multifluid flow simulators, which might cause an overestimation of the volume of NAPL that reaches the ground water. Mobile NAPL is defined as being continuous in the pore space and flows under a pressure gradient or gravitational body force. Entrapped NAPL is defined as being occluded by the aqueous phase, occurring as immobile ganglia surrounded by aqueous phase in the pore space and formed when NAPL is replaced by the aqueous phase. Residual NAPL is defined as immobile, nonwater entrapped NAPL that does not drain from the pore spaces and is conceptualized as being either continuous or discontinuous. Free NAPL comprises mobile and residual NAPL. The numerical model is formulated on mass conservation equations for oil and water, transported via NAPL and aqueous phases through variably saturated porous media. To account for phase transitions, a primary variable switching scheme is implemented for the oil-mass conservation equation over three phase conditions: (1) aqueous or aqueous-gas with dissolved oil, (2) aqueous or aqueous-gas with entrapped NAPL, and (3) aqueous or aqueous gas with free NAPL. Two laboratory-scale column experiments are modeled to verify the numerical model. Comparisons between the numerical simulations and experiments demonstrate the necessity to include the residual NAPL formation process in multifluid flow simulators.     

Transport and fate of nonaqueous phase liquid (NAPL) in variably saturated porous media with evolving scales of heterogeneity

 A stochastic simulation is performed to study multiphase flow and contaminant transport in fractal porous media with evolving scales of heterogeneity. Numerical simulations of residual NAPL mass transfer and subsequent transport of dissolved and/or volatilized NAPL mass in variably saturated media are carried out in conjunction with Monte Carlo techniques. The impact of fractal dimension, plume scale and anisotropy (stratification) of fractal media on relative dispersivities is investigated and discussed. The results indicate the significance of evolving scale of porous media heterogeneity to the NAPL transport in the subsurface. In general, the fractal porous media enhance the dispersivities of NAPL mass plume transport in both the water phase and the gas phase while the influence on the water phase is more significant. The porous media with larger fractal dimension have larger relative dispersivities. The aqueous horizontal dispersivity exhibits a most significant increase against the plume scale.     

Measurement of LNAPL flow using single-well tracer dilution techniques

This paper describes the use of single-well tracer dilution techniques to resolve the rate of light nonaqueous phase liquid (LNAPL) flow through wells and the adjacent geologic formation. Laboratory studies are presented in which a fluorescing tracer is added to LNAPL in wells. An in-well mixer keeps the tracer well mixed in the LNAPL. Tracer concentrations in LNAPL are measured through time using a fiber optic cable and a spectrometer. Results indicate that the rate of tracer depletion is proportional to the rate of LNAPL flow through the well and the adjacent formation. Tracer dilution methods are demonstrated for vertically averaged LNAPL Darcy velocities of 0.00048 to 0.11 m/d and LNAPL thicknesses of 9 to 24 cm. Over the range of conditions studied, results agree closely with steady-state LNAPL flow rates imposed by pumping. A key parameter for estimating LNAPL flow rates in the formation is the flow convergence factor alpha. Measured convergence factors for 0.030-inch wire wrap, 0.030-inch-slotted polyvinyl chloride (PVC), and 0.010-inch-slotted PVC are 1.7, 0.91, and 0.79, respectively. In addition, methods for using tracer dilution data to determine formation transmissivity to LNAPL are presented. Results suggest that single-well tracer dilution techniques are a viable approach for measuring in situ LNAPL flow and formation transmissivity to LNAPL.     

The influence of wettability on NAPL dissolution fingering

Laboratory experiments and numerical simulations in homogeneous porous media were used to investigate the influence of porous medium wettability on the formation and growth of preferential dissolution pathways, dissolution fingers, during nonaqueous phase liquid (NAPL) dissolution. As the porous medium became increasingly NAPL-wet, dissolution fingers grew wider and slower. This result was observed in physical experiments with 0% and 100% NAPL-wet conditions and confirmed with numerical simulations at these and intermediate wettabilities. A previously derived expression for an upscaled mass transfer rate coefficient that accounts for the growth of dissolution fingers was used to quantify the effect of fingering on overall NAPL removal rates. For the test cases evaluated, NAPL dissolution fingering controlled the overall rate of NAPL dissolution after the dissolution front moved 4 cm in 0% NAPL-wet conditions and 18 cm in 100% NAPL-wet conditions. Thus, even in completely NAPL-wet media dissolution fingering may control the overall rate of NAPL dissolution after relatively short travel distances. The importance of NAPL dissolution fingering in heterogeneous systems with spatially varying NAPL saturations, though, remains an important question for future work.     

Acoustic and electromagnetic properties of soils saturated with salt water and NAPL

Electrical, seismic, and electromagnetic methods can be used for noninvasive determination of subsurface physical and chemical properties. In particular, we consider the evaluation of water salinity and the detection of surface contaminants. Most of the relevant properties are represented by electric conductivity, P-wave velocity, and dielectric permittivity. Hence, it is important to obtain relationships between these measurable physical quantities and soil composition, saturation, and frequency. Conductivity in the geoelectric frequency range is obtained with Pride's model for a porous rock. (The model considers salinity and permeability.) White's model of patchy saturation is used to calculate the P-wave velocity and attenuation. Four cases are considered: light nonaqueous phase liquid (LNAPL) pockets in water, dense nonaqueous phase liquid (DNAPL) pockets in water, LNAPL pockets in air, and DNAPL pockets in air. The size of the pockets (or pools), with respect to the signal wavelength, is modeled by the theory. The electromagnetic properties in the GPR frequency range are obtained by using the Hanai–Bruggeman equation for two solids (sand and clay grains) and two fluids (LNAPL or DNAPL in water or air). The Hanai–Bruggeman exponent (1/3 for spherical particles) is used as a fitting parameter and evaluated for a sand/clay mixture saturated with water.Pride's model predicts increasing conductivity for increasing salinity and decreasing permeability. The best-fit exponent of the Hanai–Bruggeman equation for a sand/clay mixture saturated with water is 0.61, indicating that the shape of the grains has a significant influence on the electromagnetic properties. At radar frequencies, it is possible to distinguish between a water-saturated medium and a NAPL-saturated medium, but LNAPL- and DNAPL-saturated media have very similar electromagnetic properties. The type of contaminant can be better distinguished from the acoustic properties. P-wave velocity increases with frequency, and has dissimilar behaviour for wet and dry soils.     

非混溶饱和两相渗流与孔隙介质耦合作用的理论研究——数学模型

基于连续介质力学的一般理论,在黑油多相渗流理论模型的基础上,推导了非混溶饱和两相渗流与变形孔隙介质耦合作用的数学模型方程,并对流固耦合效应进行了初步的讨论     

Influence of mass transfer resistance on detection of nonaqueous phase liquids with partitioning tracer tests

Partitioning interwell tracer tests (PITTs) are a relatively new technique for measuring the amount of nonaqueous phase liquid (NAPL) within saturated porous media. In this work we examined the influence of mass transfer limitations on the accuracy of measured NAPL from PITTs. Two mathematical models were used along with laboratory column experiments to explore the influence of tracer partition coefficient, tracer detection limit, and injected tracer mass on NAPL measurements. When dimensionless mass transfer coefficients were small, NAPL measurement errors decreased with decreasing tracer partition coefficient, decreasing tracer detection limit, and increasing injected tracer mass. Extrapolating breakthrough curves exponentially reduced but did not eliminate systematic errors in NAPL measurement. Although transport in a single stream tube was used in the mathematical models and laboratory experiments, the results from this simplified domain were supported by data taken from a three-dimensional computational experiment, where the NAPL resided as large pool. Based on these results, we suggest guidelines for interpreting tracer breakthrough data to ascertain the importance of mass transfer limitations on NAPL measurements.     

An equivalent medium model for wave simulation in fractured porous rocks

Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.     

A tracer dilution method for fracture characterization in bedrock wells

This investigation was undertaken to develop an integrated method of downhole fracture characterization using a tracer. The method presented can be used to locate water-bearing fractures that intersect the well, to determine the ambient fracture flow rate and hydraulic head, and to calculate fracture transmissivity. The method was tested in two fractured crystalline bedrock wells located at the University of Connecticut in Storrs. The method entails injecting a tracer (uranine dye) into the well, while at the same time water is pumped out of the well. After steady-state conditions are reached, a borehole tracer concentration profile is developed. The dilution of the tracer is used to locate the inflowing fractures and to determine their flow rate. The fracture flow rate, plus the drawdown in the well, is then used to determine the fracture hydraulic head, transmissivity, and ambient flow rate.     

Analytical and Numerical Modeling of Solute Intrusion,Recovery, and Rebound in Fractured Bedrock

Contaminated groundwater in fractured bedrock can expose ecosystems to undesired levels of risk for extended periods due to prolonged back-diffusion from rock matrix to permeable fractures. Therefore, it is key to characterize the diffusive mass loading (intrusion) of contaminants into the rock matrix for successful management of contaminated bedrock sites. Even the most detailed site characterization techniques often fail to delineate contamination in rock matrix. This study presents a set of analytical solutions to estimate diffusive mass intrusion into matrix blocks, it is recovered by pumping and concentration rebound when pumping ceases. The analytical models were validated by comparing the results with (1) numerical model results using the same model parameters and (2) observed chloride mass recovery, rebound concentration, and concentration in pumped groundwater at a highly fractured bedrock site in Alberta, Canada. It is also demonstrated that the analytical solutions can be used to estimate the total mass stored in the fractured bedrock prior to any remediation thereby providing insights into site contamination history. The predictive results of the analytical models clearly show that successful remediation by pumping depends largely on diffusive intrusion period. The results of initial mass from the analytical model was used to successfully calibrate a three-dimensional discrete fracture network numerical model further highlighting the utility of the simple analytical solutions in supplementing the more detailed site numerical modeling. Overall, the study shows the utility of simple analytical methods to support long-term management of a contaminated fractured bedrock site including site investigations and complex numerical modeling.     

Bioremediation in Fractured Rock: 1. Modeling to Inform Design,Monitoring, and Expectations

Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes.     

Matrix diffusion-derived plume attenuation in fractured bedrock

Matrix diffusion can attenuate the rate of plume migration in fractured bedrock relative to the rate of ground water flow for both conservative and nonconservative solutes of interest. In a system of parallel, equally spaced constant aperture fractures subject to steady-state ground water flow and an infinite source width, the degree of plume attenuation increases with time and travel distance, eventually reaching an asymptotic level. The asymptotic degree of plume attenuation in the absence of degradation can be predicted by a plume attenuation factor, beta, which is readily estimated as R' (phi(m)/phi(f)), where R' is the retardation factor in the matrix, phi(m) is the matrix porosity, and phi(f) is the fracture porosity. This dual-porosity relationship can also be thought of as the ratio of primary to secondary porosity. Beta represents the rate of ground water flow in fractures relative to the rate of plume advance. For the conditions examined in this study, beta increases with greater matrix porosity, greater matrix fraction organic carbon, larger fracture spacing, and smaller fracture aperture. These concepts are illustrated using a case study where dense nonaqueous phase liquid in fractured sandstone produced a dissolved-phase trichloroethylene (TCE) plume approximately 300 m in length. Transport parameters such as matrix porosity, fracture porosity, hydraulic gradient, and the matrix retardation factor were characterized at the site through field investigations. In the fractured sandstone bedrock examined in this study, the asymptotic plume attenuation factors (beta values) for conservative and nonconservative solutes (i.e., chloride and TCE) were predicted to be approximately 800 and 12,210, respectively. Quantitative analyses demonstrate that a porous media (single-porosity) solute transport model is not appropriate for simulating contaminant transport in fractured sandstone where matrix diffusion occurs. Rather, simulations need to be conducted with either a discrete fracture model that explicitly incorporates matrix diffusion, or a dual-continuum model that accounts for mass transfer between mobile and immobile zones. Simulations also demonstrate that back diffusion from the matrix to fractures will likely be the time-limiting factor in reaching ground water cleanup goals in some fractured bedrock environments.     

Characterization of Crystalline Bedrock Contaminated by Dense Nonaqueous Liquid

A systematic hydrogeologic site characterization has been completed in a fractured rock flow system, with the objective of identifying contaminant migration and fate pathways from a historical release of 1,1,1-trichloroethane (TCA). The study integrated hydrogeologic analysis techniques such as borehole geophysical logging, pumping test analysis, and hydrochemical facies analysis to study the impact of a dense nonaqueous phase liquid (DNAPL) in a sparsely fractured crystalline bedrock. The assessment methodology can be divided into two parts: (1) characterization of the source area, where DNAPL is acting as a residual source of TCA, and (2) characterization of the downgradient plume. Reduction in DNAPL mass in the source area has resulted in significant and sustained reductions in downgradient concentrations, suggesting that remediation of fractured crystalline bedrock contaminated with DNAPL is possible and not "technically infeasible."     

Simulation of nonaqueous phase groundwater contamination

Simulation of groundwater contamination by nonaqueous phase contaminants requires solution of the multiphase flow equations for a porous medium. A three phase, two dimensional, one component nonaqueous model is developed. The contaminant can partition between the air, water and nonaqueous phase. An adaptive implicit method is used to discretize the equations. A variable substitution method, which takes into account subregions where the nonaqueous phase is nonexistent, is employed to more efficiently solve the model equations. A detailed discussion of the boundary conditions is also presented.     

轻非水相液体污染过程的高密度电阻率成像法室内监测

为探讨高密度电阻率成像法监测多孔介质中轻非水相液体迁移过程的有效性,本文通过三维砂槽进行了非饱和带中轻非水相液体的污染试验,并利用高密度电阻率成像法进行了同步的动态监测.试验之后,将砂槽层层挖开,通过数码成像,获取了污染区域的实际范围与形状.结果表明,由高密度电阻率成像法圈定的污染区域在范围与形状上都与实际的结果比较接近,并可通过三维电阻率相对值的时间变化明显的看出轻非水相液体的污染过程.这说明利用高密度电阻率成像法对非饱和多孔介质中轻非水相液体的空间分布范围进行圈定并监测其迁移过程是完全可行的.