A parametric transfer function methodology for analyzing reactive transport in nonuniform flow

We analyze reactive transport during in-situ bioremediation in a nonuniform flow field, involving multiple extraction and injection wells, by the method of transfer functions. Gamma distributions are used as parametric models of the transfer functions. Apparent parameters of classical transport models may be estimated from those of the gamma distributions by matching temporal moments. We demonstrate the method by application to measured data taken at a field experiment on bioremediation conducted in a multiple-well system in Oak Ridge, TN. Breakthrough curves (BTCs) of a conservative tracer (bromide) and a reactive compound (ethanol) are measured at multi-level sampling (MLS) wells and in extraction wells. The BTCs of both compounds are jointly analyzed to estimate the first-order degradation rate of ethanol. To quantify the tracer loss, we compare the approaches of using a scaling factor and a first-order decay term. Results show that by including a scaling factor both gamma distributions and inverse-Gaussian distributions (transfer functions according to the advection-dispersion equation) are suitable to approximate the transfer functions and estimate the reactive rate coefficients for both MLS and extraction wells. However, using a first-order decay term for tracer loss fails to describe the BTCs at the extraction well, which is affected by the nonuniform distribution of travel paths.     

3-D numerical evaluation of density effects on tracer tests

In this paper we present numerical simulations carried out to assess the importance of density-dependent flow on tracer plume development. The scenario considered in the study is characterized by a short-term tracer injection phase into a fully penetrating well and a natural hydraulic gradient. The scenario is thought to be typical for tracer tests conducted in the field. Using a reference case as a starting point, different model parameters were changed in order to determine their importance to density effects. The study is based on a three-dimensional model domain. Results were interpreted using concentration contours and a first moment analysis. Tracer injections of 0.036 kg per meter of saturated aquifer thickness do not cause significant density effects assuming hydraulic gradients of at least 0.1%. Higher tracer input masses, as used for geoelectrical investigations, may lead to buoyancy-induced flow in the early phase of a tracer test which in turn impacts further plume development. This also holds true for shallow aquifers. Results of simulations with different tracer injection rates and durations imply that the tracer input scenario has a negligible effect on density flow. Employing model cases with different realizations of a log conductivity random field, it could be shown that small variations of hydraulic conductivity in the vicinity of the tracer injection well have a major control on the local tracer distribution but do not mask effects of buoyancy-induced flow.     

3-D numerical evaluation of density effects on tracer tests

In this paper we present numerical simulations carried out to assess the importance of density-dependent flow on tracer plume development. The scenario considered in the study is characterized by a short-term tracer injection phase into a fully penetrating well and a natural hydraulic gradient. The scenario is thought to be typical for tracer tests conducted in the field. Using a reference case as a starting point, different model parameters were changed in order to determine their importance to density effects. The study is based on a three-dimensional model domain. Results were interpreted using concentration contours and a first moment analysis. Tracer injections of 0.036 kg per meter of saturated aquifer thickness do not cause significant density effects assuming hydraulic gradients of at least 0.1%. Higher tracer input masses, as used for geoelectrical investigations, may lead to buoyancy-induced flow in the early phase of a tracer test which in turn impacts further plume development. This also holds true for shallow aquifers. Results of simulations with different tracer injection rates and durations imply that the tracer input scenario has a negligible effect on density flow. Employing model cases with different realizations of a log conductivity random field, it could be shown that small variations of hydraulic conductivity in the vicinity of the tracer injection well have a major control on the local tracer distribution but do not mask effects of buoyancy-induced flow.     

Single-well reactive tracer test and stable isotope analysis for determination of microbial activity in a fast hydrocarbon-contaminated aquifer

Single-well reactive tracer tests, such as the push-pull test are useful tools for characterising in-situ bioattenuation processes in contaminated aquifers. However, the analytical models that are used to interpret push-pull data may be over-simplified, and potentially overlook important processes responsible for the frequent discrepancy between predicted and observed results obtained from push-pull tests. In this study, the limitations underlying the push-pull test methodology were investigated and were supported with results from a push-pull test conducted in a sulphate-reducing aquifer contaminated by crude oil. Poor (<7%) mass recovery was achieved when the push-pull test was performed in a fast-flowing aquifer, preventing a quantifiable reaction rate to be determined. Breakthrough curve data were unexplainable using simplified analytical models, but exhibited trends analogous with tests conducted by others, when >20% mass recoveries were achieved. Push-pull test data collected from sulphate-reducing aquifers indicate that the assumption of a well-mixed batch reactor system is incorrect and that reaction rates obtained from push-pull tests in such systems may be affected by the extraction regime implemented. Evidence of microbial respiration of the reactive tracer was provided by stable sulphur isotope analysis, from which an isotope fractionation factor of +9.9 +/- 8.1 per thousand was estimated. The stable isotope data support the argument that reaction rates calculated using push-pull tests are not uniformly distributed in space and time and are likely to be influenced by heterogeneities in the flow field.     

Use of tracer tests to evaluate the impact of enhanced-solubilization flushing on in-situ biodegradation

Tracer tests were conducted to evaluate the effect of a complexing sugar flush (CSF) on in-situ biodegradation potential at a site contaminated by jet fuel, solvents, and other organic compounds. Technical-grade hydroxypropyl-beta-cyclodextrin was used during the CSF study, which was conducted in a hydraulically isolated cell emplaced in a surficial aquifer. In-situ biodegradation potential was assessed with the use of tracer tests, which were conducted prior to and immediately following the CSF study. Ethanol, hexanol, and benzoate were used as the biodegradable tracers, while bromide was used as a nonreactive tracer. The results indicate that the biodegradation of benzoate was similar for both tracer tests. Conversely, the biodegradation of ethanol (23% increase) and hexanol (41% increase) was greater for the post-CSF tracer test. In addition, analysis of core samples collected from within the test cell indicates that the population density of aerobic jet-fuel degraders increased in the vicinity of the injection wells during the CSF. These results indicate that the cyclodextrin flush did not deleteriously affect the indigenous microbial community. This study illustrates that tracer tests can be used to evaluate the impact of remediation activities on in-situ biodegradation potential.     

Tracer test at El Berrocal site

Tracer tests provide highly valuable information about the transport properties of saturated rocks which is essential to the characterization of a potential radioactive waste repository site. In the frame of El Berrocal project, a set of tracer tests was performed in a complex geometry of inclined boreholes, combined with highly fractured transmissive zones. The aims of the tracer test programme were to gain experience, knowledge and insight into field transport experiments. To achieve this a detailed programme of new instrumentation design, site characterization and laboratory tasks was developed. For field monitoring a new electronic system was developed. The system is able to measure up to 256 parameters per borehole, with surface equipment to control pumping rates and physical and chemical parameters at both injection and extraction boreholes. The experiments progressed from single borehole dilution tests under both natural flow and forced gradient conditions to convergent flow tracer tests. Dilution tests helped to discriminate the most suitable borehole sections at which to inject tracers. The tracers were selected by the results of the laboratory programme. Uranine (fluorescein), DTPA-gadolinium (diethylenetriaminopentacetic acid, gadolinium (III)), and deuterium were injected simultaneously in one borehole section and recovered at another borehole 20 m away, pumping at a flowrate of 0.1 1 min⁻¹. First results showed a thickness porosity of 1.2 × 10⁻³ m and a longitudinal dispersivity of 17.0 m using uranine data acquired over a period of 4 d, at which point the recovery concentration had reached a maximum. However, gadolinium and deuterium appeared to travel faster, arriving at peak values after only 2 d of injection.     

Fluorescent dye imaging of the volume sampled by single well forced-gradient tracer tests evaluated in a laboratory-scale aquifer physical model

This study presents a new method to visualise forced-gradient tracer tests in 2-D using a laboratory-scale aquifer physical model. Experiments were designed to investigate the volume of aquifer sampled in vertical dipole flow tracer tests (DFTT) and push-pull tests (PPT), using a miniature monitoring well and straddle packer arrangement equipped with solute injection and recovery chambers. These tests have previously been used to estimate bulk aquifer hydraulic and transport properties for the evaluation of natural attenuation and other remediation approaches. Experiments were performed in a silica glass bead-filled box, using a fluorescent tracer (fluorescein) to deduce conservative solute transport paths. Digital images of fluorescein transport were captured under ultraviolet light and processed to analyse tracer plume geometry and obtain point-concentration breakthrough histories. Inorganic anion mixtures were also used to obtain conventional tracer breakthrough histories. Concentration data from the conservative tracer breakthrough curves was compared with the digital images and a well characterised numerical model. The results show that the peak tracer breakthrough response in dipole flow tracer tests samples a zone of aquifer close to the well screen, while the sampling volume of push-pull tests is limited by the length of the straddle packers used. The effective sampling volume of these single well forced-gradient tests in isotropic conditions can be estimated with simple equations. The experimental approach offers the opportunity to evaluate under controlled conditions the theoretical basis, design and performance of DFTTs and PPTs in porous media in relation to measured flow and transport properties.     

Controlled release,blind tests of DNAPL characterization using partitioning tracers

The partitioning tracer technique for dense nonaqueous phase liquid (DNAPL) characterization was evaluated in an isolated test cell, in which controlled releases of perchloroethylene (PCE) had occurred. Four partitioning tracer tests were conducted, two using an inverted, double five-spot pumping pattern, and two using vertical circulation wells. Two of the four tests were conducted prior to remedial activities, and two were conducted after. Each test was conducted as a "blind test" where researchers conducting the partitioning tracer tests had no knowledge of the volume, method of release, nor resulting spatial distribution of DNAPL. Multiple partitioning tracers were used in each test, and the DNAPL volume estimates varied significantly within each test based on the different partitioning tracers. The tracers with large partitioning coefficients generally predicted a smaller volume of PCE than that expected based on the actual release volume. However, these predictions were made for low DNAPL saturations (average saturation was approximately 0.003), under conditions near the limits of the method's application. Furthermore, there were several factors that may have hindered prediction accuracy, including tracer degradation and remedial fluid interference.     

A multi-directional tracer test in the fractured Chalk aquifer of E. Yorkshire, UK

A multi-borehole radial tracer test has been conducted in the confined Chalk aquifer of E. Yorkshire, UK. Three different tracer dyes were injected into three injection boreholes and a central borehole, 25 m from the injection boreholes, was pumped at 330 m(3)/d for 8 days. The breakthrough curves show that initial breakthrough and peak times were fairly similar for all dyes but that recoveries varied markedly from 9 to 57%. The breakthrough curves show a steep rise to a peak and long tail, typical of dual porosity aquifers. The breakthrough curves were simulated using a 1D dual porosity model. Model input parameters were constrained to acceptable ranges determined from estimations of matrix porosity and diffusion coefficient, fracture spacing, initial breakthrough times and bulk transmissivity of the aquifer. The model gave equivalent hydraulic apertures for fractures in the range 363-384 microm, dispersivities of 1 to 5 m and matrix block sizes of 6 to 9 cm. Modelling suggests that matrix block size is the primary controlling parameter for solute transport in the aquifer, particularly for recovery. The observed breakthrough curves suggest results from single injection-borehole tracer tests in the Chalk may give initial breakthrough and peak times reasonably representative of the aquifer but that recovery is highly variable and sensitive to injection and abstraction borehole location. Consideration of aquifer heterogeneity suggests that high recoveries may be indicative of a high flow pathway adjacent, but not necessarily connected, to the injection and abstraction boreholes whereas low recoveries may indicate more distributed flow through many fractures of similar aperture.     

Laboratory longitudinal diffusion tests: 2. Parameter estimation by inverse analysis

This study focuses on the verification of test interpretations for different state analyses of diffusion experiments. Part 1 of this study identified that steady, quasi-steady and equilibrium state analyses for the through- and in-diffusion tests with solution reservoirs are generally feasible where the tracer is not highly sorptive. In Part 2 we investigate parameter identifiability in transient-state analysis of reservoir concentration variation using a numerical approach. For increased generality, the analytical models, objective functions and Jacobian matrix necessary for inverse analysis of transient-state data are reformulated using unified dimensionless parameters. In these dimensionless forms, the number of unknown parameters is reduced and a single dimensionless parameter represents the sorption property. The dimensionless objective functions are evaluated for individual test methods and parameter identifiability is discussed in relation to the sorption property. The effects of multiple minima and measurement error on parameter identifiability are also investigated. The main findings are that inverse problems for inlet and outlet reservoir concentration analyses are generally unstable and well-posed, respectively. Where the tracer is sorptive, the inverse problem for the inlet reservoir concentration analysis may have multiple minima. When insufficient measurement data is collected, multiple solutions may result and this should be taken into consideration when inversely analyzing data including that of inlet reservoir concentration. Verification of test interpretation by cross-checking different state analyses is feasible where the tracer is not highly sorptive. In an actual experiment, test interpretation validity is demonstrated through consistency between theory and practice for different state analyses.     

Diffusive partitioning tracer test for the quantification of nonaqueous phase liquid (NAPL) in the vadose zone: Performance evaluation for heterogeneous NAPL distribution

A partitioning tracer test based on gas-phase diffusion in the vadose zone yields estimates of the residual nonaqueous phase liquid (NAPL) saturation. The present paper investigates this technique further by studying diffusive tracer breakthrough curves in the vadose zone for a heterogeneous NAPL distribution. Tracer experiments were performed in a lysimeter with a horizontal layer of artificial kerosene embedded in unsaturated sand. Tracer disappearance curves at the injection point and tracer breakthrough curves at some distance from the injection point were measured inside and outside of the NAPL layer. A numerical code was used to generate independent model predictions based on the physicochemical sand, NAPL, and tracer properties. The measured and modeled tracer breakthrough curves were in good agreement confirming the validity of important modeling assumptions such as negligible sorption of chlorofluorocarbon (CFC) tracers to the uncontaminated sand and their fast reversible partitioning between the soil air and the NAPL phase. Subsequently, the model was used to investigate different configurations of NAPL contamination. The experimental and model results show that the tracer disappearance curves of a single-well diffusive partitioning tracer test (DPTT) are dominated by the near-field presence of NAPL around the tip of the soil gas probe. In contrast, breakthrough curves of inter-well tracer tests reflect the NAPL saturation in between the probes, although there is no unique interpretation of the tracer signals if the NAPL distribution is heterogeneous. Numerical modeling is useful for the planning of a DPTT application. Simulations suggest that several cubic meters of soil can be investigated with a single inter-well partitioning tracer test of 24-hour duration by placing the injection point in the center of the investigated soil volume and probes at up to 1 m distance for the monitoring of gaseous tracers.     

Convergent-flow tracer tests in heterogeneous media: combined experimental-numerical analysis for determination of equivalent transport parameters

In modeling transport within naturally heterogeneous aquifers, it is usually assumed that the transport equations valid at local scales can also be applied at larger scales. At larger scales, the heterogeneous domain is represented by an equivalent homogeneous medium. Convergent-flow tracer tests constitute one of the most frequently used field tests to estimate effective input parameters of equivalent homogeneous aquifers. Traditionally, statistical approaches applied to groundwater flow and solute transport have provided tools to estimate these equivalent parameters. These approaches are based on a number of simplifications including the assumption that the point transmissivity values follow a multilog-normal random function. Several investigators have found that this assumption may not be valid in many field cases. In order to study the applicability of the equivalent homogeneous formulation in a nontraditional stochastic field, a number of experimental and numerical studies were conducted. The results are used to determine the apparent values of porosity and dispersivity that would be obtained if convergent-flow tracer tests were conducted in a deterministically generated heterogeneous transmissivity field displaying anisotropy in the correlation structure. It is shown that in this particular heterogeneous media, apparent porosity strongly depends on connectivity rather than on transmissivity. This dependence on connectivity questions the theoretical results obtained in continuum equivalent fields to estimate effective porosity.     

An application of the theory of kinematics of mixing to the study of tropospheric dispersion

The most common technique used for numerical simulations of tracer mixing is that of the numerical solution of the advection–diffusion equation with the unresolved fluxes parameterized using the similarity theory. Despite correct predictions of the overall directions of transport, models based on a numerical solution of the advection–diffusion equation lack sufficient accuracy to correctly reproduce the coupling of mixing with small scale processes which are sensitive to the microstructure of the tracer distribution. The objective of this paper is to revisit the basic formalism employed in numerical models used to investigate atmospheric tracers. The main mathematical method proposed here is the theory of kinematics of mixing which could be applied effectively for simulations of atmospheric transport processes. At the beginning of the paper, we introduce simple mathematical transformations in order to demonstrate how complex topological structures are created by mixing processes. These idealistic flow systems are essential to explain transport properties of much more complex three-dimensional geophysical flows. An example of the application of the kinematics of mixing to the analysis of tracer transport on a planetary scale is presented in the following sections. The complex filamentary structures simulated in the numerical experiment are evaluated using some commonly applied statistical measures in order to compare the results with the data published in the literature. The results of the experiment are also analysed with the help of simple conceptual models of fluid filaments. The microstructure of the tracer distribution introduced in the paper is essential to increase our understanding of atmospheric transport and to develop more realistic parameterizations of small-scale mixing. The presented results could also be used to improve calculations of the coupling between microphysical processes and tracer mixing.     

Modeling field-scale multiple tracer injection at a low-level waste disposal site in fractured rocks: effect of multiscale heterogeneity and source term uncertainty on conceptual understanding of mass transfer processes

Multiple factors may affect the scale-up of laboratory multi-tracer injection into structured porous media to the field. Under transient flow conditions and with multiscale heterogeneities in the field, previous attempts to scale-up laboratory experiments have not answered definitely the questions about the governing mechanisms and the spatial extent of the influence of small-scale mass transfer processes such as matrix diffusion. The objective of this research is to investigate the effects of multiscale heterogeneity, mechanistic and site model conceptualization, and source term density effect on elucidating and interpreting tracer movement in the field. Tracer release and monitoring information previously obtained in a field campaign of multiple, conservative tracer injection under natural hydraulic gradients at a low-level waste disposal site in eastern Tennessee, United States, is used for the research. A suite of two-pore-domain, or fracture-matrix, groundwater flow and transport models are calibrated and used to conduct model parameter and prediction uncertainty analyses. These efforts are facilitated by a novel nested Latin-hypercube sampling technique. Our results verify, at field scale, a multiple-pore-domain, multiscale mechanistic conceptual model that was used previously to interpret only laboratory observations. The results also suggest that, integrated over the entire field site, mass flux rates attributable to small-scale mass transfer are comparable to that of field-scale solute transport. The uncertainty analyses show that fracture spacing is the most important model parameter and model prediction uncertainty is relatively higher at the interface between the preferred flow path and its parent bedrock. The comparisons of site conceptual models indicate that the effect of matrix diffusion may be confined to the immediate neighborhood of the preferential flow path. Finally, because the relatively large amount of tracer needed for field studies, it is likely that source term density effect may exaggerate or obscure the effect of matrix diffusion on the movement of tracers from the preferred flow path into the bedrock.     

Field-scale colloid migration experiments in a granite fracture

An understanding of particle migration in fractured rock, required to assess the potential for colloid-facilitated transport of radionuclides, can best be evaluated when the results of laboratory experiments are demonstrated in the field. Field-scale migration experiments with silica colloids were carried out at AECL's Underground Research Laboratory (URL), located in southern Manitoba, to develop the methodology for large-scale migration experiments and to determine whether colloid transport is possible over distances up to 17 m. In addition, these experiments were designed to evaluate the effects of flow rate and flow path geometry, and to determine whether colloid tracers could be used to provide additional information on subsurface transport to that provided by conservative tracers alone. The colloid migration studies were carried out as part of AECL's Transport Properties in Highly Fractured Rock Experiment, the objective of which was to develop and demonstrate methods for evaluating the solute transport characteristics of zones of highly fractured rock. The experiments were carried out within fracture zone 2 as two-well recirculating, two-well non-recirculating, and convergent flow tests, using injection rates of 5 and 10 1 min⁻¹. Silica colloids with a 20 nm size were used because they are potentially mobile due to their stability, small size and negative surface charge. The shapes of elution profiles for colloids and conservative tracers were similar, demonstrating that colloids can migrate over distances of 17 m. The local region of drawdown towards the URL shaft affected colloid migration and, to a lesser extent, conservative tracer migration within the flow field established by the two-well tracer tests. These results indicate that stable colloids, with sizes as small as 20 nm, have different migration properties from dissolved conservative tracers. - 1997 Atomic Energy of Canada.     

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.     

Kinetic limitations on tracer partitioning in ganglia dominated source zones

Quantification of the relationship between dense nonaqueous phase liquid (DNAPL) source strength, source longevity and spatial distribution is increasingly recognized as important for effective remedial design. Partitioning tracers are one tool that may permit interrogation of DNAPL architecture. Tracer data are commonly analyzed under the assumption of linear, equilibrium partitioning, although the appropriateness of these assumptions has not been fully explored. Here we focus on elucidating the nonlinear and nonequilibrium partitioning behavior of three selected alcohol tracers - 1-pentanol, 1-hexanol and 2-octanol in a series of batch and column experiments. Liquid-liquid equilibria for systems comprising water, TCE and the selected alcohol illustrate the nonlinear distribution of alcohol between the aqueous and organic phases. Complete quantification of these equilibria facilitates delineation of the limits of applicability of the linear partitioning assumption, and assessment of potential inaccuracies associated with measurement of partition coefficients at a single concentration. Column experiments were conducted under conditions of non-equilibrium to evaluate the kinetics of the reversible absorption of the selected tracers in a sandy medium containing a uniform entrapped saturation of TCE-DNAPL. Experimental tracer breakthrough data were used, in conjunction with mathematical models and batch measurements, to evaluate alternative hypotheses for observed deviations from linear equilibrium partitioning behavior. Analyses suggest that, although all tracers accumulate at the TCE-DNAPL/aqueous interface, surface accumulation does not influence transport at concentrations typically employed for tracer tests. Moreover, results reveal that the kinetics of the reversible absorption process are well described using existing mass transfer correlations originally developed to model aqueous boundary layer resistance for pure-component NAPL dissolution.     

A new tracer technique for monitoring groundwater fluxes: the Finite Volume Point Dilution Method

Quantification of pollutant mass fluxes is essential for assessing the impact of contaminated sites on their surrounding environment, particularly on adjacent surface water bodies. In this context, it is essential to quantify but also to be able to monitor the variations with time of Darcy fluxes in relation with changes in hydrogeological conditions and groundwater - surface water interactions. A new tracer technique is proposed that generalizes the single-well point dilution method to the case of finite volumes of tracer fluid and water flush. It is called the Finite Volume Point Dilution Method (FVPDM). It is based on an analytical solution derived from a mathematical model proposed recently to accurately model tracer injection into a well. Using a non-dimensional formulation of the analytical solution, a sensitivity analysis is performed on the concentration evolution in the injection well, according to tracer injection conditions and well-aquifer interactions. Based on this analysis, optimised field techniques and interpretation methods are proposed. The new tracer technique is easier to implement in the field than the classical point dilution method while it further allows monitoring temporal changes of the magnitude of estimated Darcy fluxes, which is not the case for the former technique. The new technique was applied to two experimental sites with contrasting objectives, geological and hydrogeological conditions, and field equipment facilities. In both cases, field tracer concentrations monitored in the injection wells were used to fit the calculated modelled concentrations by adjusting the apparent Darcy flux crossing the well screens. Modelling results are very satisfactory and indicate that the methodology is efficient and accurate, with a wide range of potential applications in different environments and experimental conditions, including the monitoring with time of changes in Darcy fluxes.     

Superposition of borehole-to-surface voltage residuals for Vadose Zone plume delineation

An injected tracer field experiment was conducted at the University of Idaho Ground Water Field Laboratory to evaluate the application of borehole-to-surface voltage measurements for delineation of the tracer distribution in partially saturated, fractured basalt. A tap water tracer was injected into a fracture-dominated, salt-water plume formed during a previous salt-water injection experiment. The tap water tracer was injected into a central injection well under constant hydraulic head for 34 days. The injection well was surrounded by seven test boreholes. Each borehole contained several copper wire electrodes for borehole-to-surface potential measurements between a surface grid of 224 copper sulfate, porous pot electrodes. Eight pole-pole, borehole-to-surface voltage data sets were acquired during each measurement period by energization of a selected electrode in each of the eight boreholes. Predicted voltages for a uniform earth (homogeneous and isotropic) potential model (finite difference) were subtracted from each data set (for its respective current source location), and the voltage residuals superposed to create new data sets with greater measurement sensitivity and coverage, to aid in interpretation. These data sets were collected over four measurement periods during tap water injection and four measurement periods during the subsequent 64-day drainage phase. The data were interpreted with the use of three-dimensional models and by comparisons with other electrical and hydrological observations. Results indicate that superposition of multiple data sets of voltage residuals significantly improved the lateral resolution of subsurface bulk resistivity changes that occurred over time.     

Transport of volatile compounds in porous media in the presence of a trapped gas phase

The presence of an immobile gaseous phase in an otherwise-saturated porous medium affects the transport of volatile compounds. The linear theory of partitioning tracers suggests that a volatile tracer introduced into such a system should be retarded with a constant retardation factor. Using high concentrations, however, the saturation of the gaseous phase will change as an effect of the tracer test itself. Competitive gas transfer among all volatile compounds and the change of saturation may lead to tracer concentrations that are temporarily higher than those injected. We analyze the system in the framework of the coherence theory by Helfferich [Soc. Pet. Eng. J. 21 (1) (1981) 51]. The governing equations are formulated as functions of total concentration, i.e., the mass of solutes in all phases per pore volume. Neglecting dispersion and mass-transfer kinetics, we derive the characteristic form of the resulting system of hyperbolic equations. In a system with N volatile compounds, a variation of the concentrations splits up into N waves, each traveling with its own characteristic velocity. If the presence of a gaseous phase is sustained, one wave will be a standing one. We perform numerical model calculations for tracers with various Henry's law coefficients and show that the results agree with the semi-analytical solution obtained by coherence theory.