Analysis of a vertical dipole tracer test in highly fractured rock

The results of a vertical dipole tracer experiment performed in highly fractured rocks of the Clare Valley, South Australia, are presented. The injection and withdrawal piezometers were both screened over 3 m and were separated by 6 m (midpoint to midpoint). Due to the long screen length, several fracture sets were intersected, some of which do not connect the two piezometers. Dissolved helium and bromide were injected into the dipole flow field for 75 minutes, followed by an additional 510 minutes of flushing. The breakthrough of helium was retarded relative to bromide, as was expected due to the greater aqueous diffusion coefficient of helium. Also, only -25% of the total mass injected of both tracers was recovered. Modeling of the tracer transport was accomplished using an analytical one-dimensional flow and transport model for flow through a fracture with diffusion into the matrix. The assumptions made include: streamlines connecting the injection and withdrawal point can be modeled as a dipole of equal strength, flow along each streamline is one dimensional, and there is a constant Peclet number for each streamline. In contrast to many other field tracer studies performed in fractured rock, the actual travel length between piezometers was not known. Modeling was accomplished by fitting the characteristics of the tracer breakthrough curves (BTCs), such as arrival times of the peak concentration and the center of mass. The important steps were to determine the fracture aperture (240 microm) based on the parameters that influence the rate of matrix diffusion (this controls the arrival time of the peak concentration); estimating the travel distance (11 m) by fitting the time of arrival of the centers of mass of the tracers; and estimating fracture dispersivity (0.5 m) by fitting the times that the inflection points occurred on the front and back limbs of the BTCs. This method works even though there was dilution in the withdrawal well, the amount of which can be estimated by determining the value that the modeled concentrations need to be reduced to fit the data (approximately 50%). The use of two tracers with different diffusion coefficients was not necessary, but it provides important checks in the modeling process because the apparent retardation between the two tracers is evidence of matrix diffusion and the BTCs of both tracers need to be accurately modeled by the best fit parameters.     

Application of ERT,Saline Tracer and Numerical Studies to Delineate Preferential Paths in Fractured Granites

Accurate quantification of in situ heterogeneity and flow processes through fractured geologic media remains elusive for hydrogeologists due to the complexity in fracture characterization and its multiscale behavior. In this research, we demonstrated the efficacy of tracer-electrical resistivity tomography (ERT) experiments combined with numerical simulations to characterize heterogeneity and delineate preferential flow paths in a fractured granite aquifer. A series of natural gradient saline tracer experiments were conducted from a depth window of 18 to 22 m in an injection well (IW) located inside the Indian Institute of Technology Hyderabad campus. Tracer migration was monitored in a time-lapse mode using two cross-sectional surface ERT profiles placed in the direction of flow gradient. ERT data quality was improved by considering stacking, reciprocal measurements, resolution indicators, and geophysical logs. Dynamic changes in subsurface electrical properties inferred via resistivity anomalies were used to highlight preferential flow paths of the study area. Temporal changes in electrical resistivity and tracer concentration were monitored along the vertical in an observation well located at 48 m to the east of the IW. ERT-derived tracer breakthrough curves were in agreement with geochemical sample measurements. Fracture geometry and hydraulic properties derived from ERT and pumping tests were further used to evaluate two mathematical conceptualizations that are relevant to fractured aquifers. Results of numerical analysis conclude that dual continuum model that combines matrix and fracture systems through a flow exchange term has outperformed equivalent continuum model in reproducing tracer concentrations at the monitoring wells (evident by a decrease in RMSE from 199 to 65 mg/L). A sensitivity analysis on model simulations conclude that spatial variability in hydraulic conductivity, local-scale dispersion, and flow exchange at fracture-matrix interface have a profound effect on model simulations.     

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.     

Evaluation of Combined Direct-Push Methods Used for Aquifer Model Generation

Most established methods to characterize aquifer structure and hydraulic conductivities of hydrostratigraphical units are not capable of delivering sufficient information in the spatial resolution that is desired for sophisticated numerical contaminant transport modeling and adapted remediation design. With hydraulic investigation methods based on the direct-push (DP) technology such as DP slug tests, DP injection logging, and the hydraulic profiling tool, it is possible to rapidly delineate hydrogeological structures and estimate their hydraulic conductivity in shallow unconsolidated aquifers without the need for wells. A combined application of these tools was used for the investigation of a contaminated German refinery site and for the setup of hydraulic aquifer models. The quality of DP investigation and the models was evaluated by comparisons of tracer transport simulations using these models and measured breakthroughs of two natural gradient tracer tests. Model scenarios considering the information of all tools together showed good reproduction of the measured breakthroughs, indicating the suitability of the approach and a minor impact of potential technical limitations. Using the DP slug tests alone yielded significantly higher deviations for the determined hydraulic conductivities compared to considering two or three of the tools. Realistic aquifer models developed on basis of such combined DP investigation approaches can help optimize remediation concepts or identify flow regimes for aquifers with a complex structure.     

Introduction to hydromechanical well tests in fractured rock aquifers

This article introduces hydromechanical well tests as a viable field method for characterizing fractured rock aquifers. These tests involve measuring and analyzing small displacements along with pressure transients. Recent developments in equipment and analyses have simplified hydromechanical well tests, and this article describes initial field results and interpretations during slug and constant-rate pumping tests conducted at a site underlain by fractured biotite gneiss in South Carolina. The field data are characterized by displacements of 0.3 μm to more than 10 μm during head changes up to 10 m. Displacements are a hysteretic function of hydraulic head in the wellbore, with displacements late in a well test always exceeding those at similar wellbore pressures early in the test. Displacement measurements show that hydraulic aperture changes during well tests, and both scaling analyses and field data suggest that T changed by a few percent per meter of drawdown during slug and pumping tests at our field site. Preliminary analyses suggest that displacement data can be used to improve estimates of storativity and to reduce nonuniqueness during hydraulic well tests involving single wells.     

Estimating Flow Using Tracers and Hydraulics in Synthetic Heterogeneous Aquifers

Regional ground water flow is most usually estimated using Darcy's law, with hydraulic conductivities estimated from pumping tests, but can also be estimated using ground water residence times derived from radioactive tracers. The two methods agree reasonably well in relatively homogeneous aquifers but it is not clear which is likely to produce more reliable estimates of ground water flow rates in heterogeneous systems. The aim of this paper is to compare bias and uncertainty of tracer and hydraulic approaches to assess ground water flow in heterogeneous aquifers. Synthetic two-dimensional aquifers with different levels of heterogeneity (correlation lengths, variances) are used to simulate ground water flow, pumping tests, and transport of radioactive tracers. Results show that bias and uncertainty of flow rates increase with the variance of the hydraulic conductivity for both methods. The bias resulting from the nonlinearity of the concentration–time relationship can be reduced by choosing a tracer with a decay rate similar to the mean ground water residence time. The bias on flow rates estimated from pumping tests is reduced when performing long duration tests. The uncertainty on ground water flow is minimized when the sampling volume is large compared to the correlation length. For tracers, the uncertainty is related to the ratio of correlation length to the distance between sampling wells. For pumping tests, it is related to the ratio of correlation length to the pumping test's radius of influence. In regional systems, it may be easier to minimize this ratio for tracers than for pumping tests.     

Characterization of a multilayer aquifer using open well dilution tests

An approach to characterization of multilayer aquifer systems using open well borehole dilution is described. The approach involves measuring observation well flow velocities while a nearby extraction well is pumped by introducing a saline tracer into observation wells and collecting dilution vs. depth profiles. Inspection of tracer profile evolution allows discrete permeable layers within the aquifer to be identified. Dilution profiles for well sections between permeable layers are then converted into vertical borehole flow velocities and their evolution, using an analytic solution to the advection-dispersion equation applied to borehole flow. The dilution approach is potentially able to measure much smaller flow velocities that would be detectable using flowmeters. Vertical flow velocity data from the observation wells are then matched to those generated using a hydraulic model of the aquifer system, "shorted" by the observation wells, to yield the hydraulic properties of the constituent layers. Observation well flow monitoring of pumping tests represents a cost-effective alternative or preliminary approach to pump testing each layer of a multilayer aquifer system separately using straddle packers or screened wells and requires no prior knowledge of permeable layer depths and thicknesses. The modification described here, of using tracer dilution rather than flowmeter logging to obtain well flow velocities, allows the approach to be extended to greater well separations, thus characterizing a larger volume of the aquifer. An example of the application of this approach to a multilayer Chalk Aquifer in Yorkshire, Northeast England, is presented.     

Measurement of LNAPL Flux Using Single‐Well Intermittent Mixing Tracer Dilution Tests

The stability of subsurface Light Nonaqueous Phase Liquids (LNAPLs) is a key factor driving expectations for remedial measures at LNAPL sites. The conventional approach to resolving LNAPL stability has been to apply Darcy's Equation. This paper explores an alternative approach wherein single‐well tracer dilution tests with intermittent mixing are used to resolve LNAPL stability. As a first step, an implicit solution for single‐well intermittent mixing tracer dilution tests is derived. This includes key assumptions and limits on the allowable time between intermittent mixing events. Second, single‐well tracer dilution tests with intermittent mixing are conducted under conditions of known LNAPL flux. This includes a laboratory sand tank study and two field tests at active LNAPL recovery wells. Results from the sand tank studies indicate that LNAPL fluxes in wells can be transformed into formation fluxes using corrections for (1) LNAPL thicknesses in the well and formation and (2) convergence of flow to the well. Using the apparent convergence factor from the sand tank experiment, the average error between the known and measured LNAPL fluxes is 4%. Results from the field studies show nearly identical known and measured LNAPL fluxes at one well. At the second well the measured fluxes appear to exceed the known value by a factor of two. Agreement between the known and measured LNAPL fluxes, within a factor of two, indicates that single‐well tracer dilution tests with intermittent mixing can be a viable means of resolving LNAPL stability.     

An Unsteady State Tracer Method for Characterizing Fractures in Bedrock Wells

Evaluating contaminants impacting wells in fractured crystalline rock requires knowledge of the individual fractures contributing water. This typically involves using a sequence of tools including downhole geophysics, flow meters, and straddle packers. In conjunction with each other these methods are expensive, time consuming, and can be logistically difficult to implement. This study demonstrates an unsteady state tracer method as a cost‐effective alternative for gathering fracture information in wells. The method entails introducing tracer dye throughout the well, inducing fracture flow into the well by conducting a slug test and then profiling the tracer concentration in the well to locate water contributing fractures where the dye has been diluted. By monitoring the development of the dilution zones within the wellbore with time, the transmissivity and the hydraulic head of the water contributing fractures can be determined. Ambient flow conditions and the contaminant concentration within the fractures can also be determined from the tracer dilution. This method was tested on a large physical model well and a bedrock well. The model well was used to test the theory underlying the method and to refine method logistics. The approach located the fracture and generated transmissivity values that were in excellent agreement with those calculated by slug testing. For the bedrock well tested, two major active fractures were located. Fracture location and ambient well conditions matched results from conventional methods. Estimates of transmissivity values by the tracer method were within an order of magnitude of those calculated using heat‐pulse flow meter data.     

Heat Tracing in a Fractured Aquifer with Injection of Hot and Cold Water

Heat as a tracer in fractured porous aquifers is more sensitive to fracture-matrix processes than a solute tracer. Temperature evolution as a function of time can be used to differentiate fracture and matrix characteristics. Experimental hot (50 °C) and cold (10 °C) water injections were performed in a weathered and fractured granite aquifer where the natural background temperature is 30 °C. The tailing of the hot and cold breakthrough curves, observed under different hydraulic conditions, was characterized in a log–log plot of time vs. normalized temperature difference, also converted to a residence time distribution (normalized). Dimensionless tail slopes close to 1.5 were observed for hot and cold breakthrough curves, compared to solute tracer tests showing slopes between 2 and 3. This stronger thermal diffusive behavior is explained by heat conduction. Using a process-based numerical model, the impact of heat conduction toward and from the porous rock matrix on groundwater heat transport was explored. Fracture aperture was adjusted depending on the actual hydraulic conditions. Water density and viscosity were considered temperature dependent. The model simulated the increase or reduction of the energy level in the fracture-matrix system and satisfactorily reproduced breakthrough curves tail slopes. This study shows the feasibility and utility of cold water tracer tests in hot fractured aquifers to boost and characterize the thermal matrix diffusion from the matrix toward the flowing groundwater in the fractures. This can be used as complementary information to solute tracer tests that are largely influenced by strong advection in the fractures.     

The Effect of Hydraulic Gradient and Pattern of Conduit Systems on Tracing Tests: Bench-Scale Modeling

Tracer breakthrough curves provide valuable information about the traced media, especially in inherently heterogeneous karst aquifers. In order to study the effect of variations in hydraulic gradient and conduit systems on breakthrough curves, a bench scale karst model was constructed. The bench scale karst model contains both matrix and a conduit. Eight tracing tests were conducted under a wide range of hydraulic gradients from 1 to greater than 5 for branchwork and network-conduit systems. Sampling points at varying distances from the injection point were utilized. Results demonstrate that mean tracer velocities, tracer mass recovery and linear rising slope of the breakthrough curves were directly controlled by hydraulic gradient. As hydraulic gradient increased, both one half the time for peak concentration and one fifth the time for peak concentration decreased. The results demonstrate the variations in one half the time for peak concentration and one fifth the time for peak concentration of the descending limb for different sampling points under differing hydraulic gradients are mainly controlled by the interactions of advection with dispersion. The results are discussed from three perspectives: different conduit systems, different hydraulic-gradient conditions, and different sampling points. The research confirmed the undeniable role of hydrogeological setting (i.e., hydraulic gradient and conduit system) on the shape of the breakthrough curve. The extracted parameters (mobile-fluid velocity, tracer-mass recovery, linear rising limb, one half the time for peak concentration, and one fifth the time for peak concentration) allow for differentiating hydrogeological settings and enhance interpretations the tracing tests in karst aquifers.     

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.     

Derivative‐Assisted Classification of Fractured Zones Crossing a Deep Borehole

In this study, the derivative analysis using the derivative of drawdown with respect to log‐time was utilized to determine candidates for hydraulic conductor domains (HCDs). At a 500‐m deep borehole in the study site, the fractured rocks crossing the borehole were first classified in fractured and nonfractured zones by core logging and geophysical loggings, such as acoustic televiewing, density, and flow loggings. After conducting the hydraulic tests such as constant head withdrawal and recovery tests at the fractured zones and the nonfractured zones, the derivative analyses were carried out, of which the results were evaluated to determine the candidates for HCDs. For the nonfractured zones, the diagnostic plot has only a big hump indicating poor connection of the background fractures to the permeable geologic media, while those of the candidates for HCDs show various flow regimes. On the basis of these results, the candidates for HCDs among the fractured zones were determined. From discussion on the results, the combination of the spacing analysis and derivative analysis following a hydraulic test is recommended for determining the candidates for HCDs rather than other geophysical loggings.     

Numerical Modeling of Methane Leakage from a Faulty Natural Gas Well into Fractured Tight Formations

Horizontal drilling and hydraulic fracturing have enabled hydrocarbon recovery from unconventional reservoirs, but led to natural gas contamination of shallow groundwaters. We describe and apply numerical models of gas‐phase migration associated with leaking natural gas wells. Three leakage scenarios are simulated: (1) high‐pressure natural gas pulse released into a fractured aquifer; (2) continuous slow leakage into a tilted fractured formation; and (3) continuous slow leakage into an unfractured aquifer with fluvial channels, to facilitate a generalized evaluation of natural gas transport from faulty natural gas wells. High‐pressure pulses of gas leakage into sparsely fractured media are needed to produce the extensive and rapid lateral spreading of free gas previously observed in field studies. Transport in fractures explains how methane can travel vastly different distances and directions laterally away from a leaking well, which leads to variable levels of methane contamination in nearby groundwater wells. Lower rates of methane leakage (≤1 Mcf/day) produce shorter length scales of gas transport than determined by the high‐pressure scenario or field studies, unless aquifers have low vertical permeabilities (≤1 millidarcy) and fractures and bedding planes have sufficient tilt (～10°) to allow a lateral buoyancy component. Similarly, in fractured rock aquifers or where permeability is controlled by channelized fluvial deposits, lateral flow is not sufficiently developed to explain fast‐developing gas contamination (0‐3 months) or large length scales (～1 km) documented in field studies. Thus, current efforts to evaluate the frequency, mechanism, and impacts of natural gas leakage from faulty natural gas wells likely underestimate contributions from small‐volume, low‐pressure leakage events.     

Identifying Aquifer Type in Fractured Rock Aquifers using Harmonic Analysis

Determining aquifer type, unconfined, semi‐confined, or confined, by drilling or performing pumping tests has inherent problems (i.e., cost and complex field issues) while sometimes yielding inconclusive results. An improved method to cost‐effectively determine aquifer type would be beneficial for hydraulic mapping of complex aquifer systems like fractured rock aquifers. Earth tides are known to influence water levels in wells penetrating confined aquifers or unconfined thick, low‐porosity aquifers. Water‐level fluctuations in wells tapping confined and unconfined aquifers are also influenced by changes in barometric pressure. Harmonic analyses of water‐level fluctuations of a thick (～1000 m) carbonate aquifer located in south‐central Oklahoma (Arbuckle‐Simpson aquifer) were utilized in nine wells to identify aquifer type by evaluating the influence of earth tides and barometric‐pressure variations using signal identification. On the basis of the results, portions of the aquifer responded hydraulically as each type of aquifer even though there was no significant variation in lithostratigraphy. The aquifer type was depth dependent with confined conditions becoming more prevalent with depth. The results demonstrate that harmonic analysis is an accurate and low‐cost method to determine aquifer type.     

The Horizontal Reactive Media Treatment Well (HRX Well®) for Passive In‐Situ Remediation

A new in‐situ remediation concept termed a Horizontal Reactive Media Treatment Well (HRX Well®) is presented that utilizes horizontal wells filled with reactive media to passively treat contaminated groundwater in‐situ. The approach involves the use of large‐diameter directionally drilled horizontal wells filled with granular reactive media generally installed parallel to the direction of groundwater flow. The design leverages natural “flow‐focusing” behavior induced by the high in‐well hydraulic conductivity of the reactive media relative to the aquifer hydraulic conductivity to passively capture and treat proportionally large volumes of groundwater within the well. Clean groundwater then exits the horizontal well along its downgradient sections. Many different types of solid granular reactive media are already available (e.g., zero valent iron, activated carbon, ion exchange resins, zeolite, apatite, chitin); therefore, this concept could be used to address a wide range of contaminants. Three‐dimensional flow and transport simulations were completed to assess the general hydraulic performance, capture zones, residence times, effects of aquifer heterogeneity, and treatment effectiveness of the concept. The results demonstrate that capture and treatment widths of up to tens of feet can be achieved for many aquifer settings, and that reductions in downgradient concentrations and contaminant mass flux are nearly immediate. For a representative example, the predicted treatment zone width for the HRX Well is approximately 27 to 44 feet, and contaminant concentrations immediately downgradient of the HRX Well decreased an order of magnitude within 10 days. A series of laboratory‐scale physical tests (i.e., tank tests) were completed that further demonstrate the concept and confirm model prediction performance. For example, the breakthrough time, peak concentration and total mass recovery of methylene blue (reactive tracer) was about 2, 35, and 20 times (respectively) less than chloride (conservative tracer) at the outlet of the tank‐scale HRX Well.     

Leaky Aquifer Models to Simulate the Well Flow in Fractured Aquifers with Linear Interporosity Flow

Following Hemker and Maas (1987) the models of two or three leaky aquifers are applied to simulate the flow to vertical wells operating in the fractured or dual porosity aquifers. The software WellTest (WT) (Székely 2015) is used for calculating the drawdown and discharge rate variation. The comparative analysis with the independent analytical solutions by Boulton and Streltsova-Adams (1978), Warren and Root (1963), Kazemi et al. (1969) concluded with acceptable agreement between the WT simulation and the alternate calculation methods. The selected field tests have been conducted in fractured limestone aquifers. The pumping test west of Copenhagen shows an example of fractured aquifer with considerable negative skin effect at the well face. The flowing well Wafra W1 in Kuwait operates in the two-zone aquifer exhibiting sufficient vertical recharge via leakage beyond a circular domain of estimated radius of 2460 m.     

A Tracer Test for Detecting Cross Contamination Along a Monitoring Well Column

A tracer test was used to evaluate whether cross contamination exists along a monitoring well completed through a shallow ground water system in fractured clay and screened in a sand and gravel aquifer. The fractured clay is separated from the sand and gravel deposit by a layer of highly plastic unfractured clay. A natural vertical downward hydraulic gradient of approximately 0.5 exists between the shallow system and the sand and gravel aquifer. Ground water contamination was detected in an adjacent monitoring well screened in the fractured clay and in the monitoring well screened in the sand and gravel deposit. No ground water contamination was apparent in an intermediate well screened in the unfractured clay layer. A tracer of sodium bromide was injected into a shallow boring near the monitoring wells. The tracer was detected in the monitoring well in the sand and gravel aquifer after three to seven days. The bromide concentration continued to increase in this well with time while the concentration in the shallow boring declined. This trend of tracer concentration indicates the tracer has in fact migrated downward and possibly traveled along the well column.     

Detecting a Defective Casing Seal at the Top of a Bedrock Aquifer

An improperly sealed casing can produce a direct hydraulic connection between two or more originally isolated aquifers with important consequences regarding groundwater quantity and quality. A recent study by Richard et al. (2014) investigated a monitoring well installed in a fractured rock aquifer with a defective casing seal at the soil–bedrock interface. A hydraulic short circuit was detected that produced some leakage between the rock and the overlying deposits. A falling‐head permeability test performed in this well showed that the usual method of data interpretation is not valid in this particular case due to the presence of a piezometric error. This error is the direct result of the preferential flow originating from the hydraulic short circuit and the subsequent re‐equilibration of the piezometric levels of both aquifers in the vicinity of the inlet and the outlet of the defective seal. Numerical simulations of groundwater circulation around the well support the observed impact of the hydraulic short circuit on the results of the falling‐head permeability test. These observations demonstrate that a properly designed falling‐head permeability test may be useful in the detection of defective casing seals.