Hydraulically Controlled Discrete Sampling from Open Boreholes

Groundwater sampling from open boreholes in fractured‐rock aquifers is particularly challenging because of mixing and dilution of fluid within the borehole from multiple fractures. This note presents an alternative to traditional sampling in open boreholes with packer assemblies. The alternative system called ZONFLO (zonal flow) is based on hydraulic control of borehole flow conditions. Fluid from discrete fractures zones are hydraulically isolated allowing for the collection of representative samples. In rough‐faced open boreholes and formations with less competent rock, hydraulic containment may offer an attractive alternative to physical containment with packers. Preliminary test results indicate a discrete zone can be effectively hydraulically isolated from other zones within a borehole for the purpose of groundwater sampling using this new method.     

Delineation of Ground Water Flow in Fractured Rock in the Southwestern Part of Manhattan, New York, Through Use of Advanced Borehole Geophysical Methods

Advanced borehole-geophysical methods were used to assess the geohydrology of fractured crystalline bedrock at five test boreholes in southwestern Manhattan Island, New York, in preparation for construction of a third water tunnel for New York City. The boreholes penetrated gneiss and other crystalline bedrock that has an overall southwest to northwest dipping foliation with a 60° dip. Most of the fractures encountered are either nearly horizontal or have moderate northwest dip azimuths. Fracture indexes range from 0.25 to 0.44 fracture per foot (0.3 m) of borehole.
Electromagnetic (EM) and heat-pulse flowmeter logs obtained under ambient and pumping conditions, together with other geophysical logs, indicate transmissive fracture zones in each borehole. Pumping tests of each borehole indicated transmissivity ranges from <2 to 360 ft²/day (0.2 to 33 m²/day). Ground water appears to flow within an interconnected fracture network toward the south and west within the study area. No correlation was indicated between the fracture index and the total borehole transmissivity.     

Cross-borehole flowmeter tests for transient heads in heterogeneous aquifers

Cross-borehole flowmeter tests have been proposed as an efficient method to investigate preferential flowpaths in heterogeneous aquifers, which is a major task in the characterization of fractured aquifers. Cross-borehole flowmeter tests are based on the idea that changing the pumping conditions in a given aquifer will modify the hydraulic head distribution in large-scale flowpaths, producing measurable changes in the vertical flow profiles in observation boreholes. However, inversion of flow measurements to derive flowpath geometry and connectivity and to characterize their hydraulic properties is still a subject of research. In this study, we propose a framework for cross-borehole flowmeter test interpretation that is based on a two-scale conceptual model: discrete fractures at the borehole scale and zones of interconnected fractures at the aquifer scale. We propose that the two problems may be solved independently. The first inverse problem consists of estimating the hydraulic head variations that drive the transient borehole flow observed in the cross-borehole flowmeter experiments. The second inverse problem is related to estimating the geometry and hydraulic properties of large-scale flowpaths in the region between pumping and observation wells that are compatible with the head variations deduced from the first problem. To solve the borehole-scale problem, we treat the transient flow data as a series of quasi-steady flow conditions and solve for the hydraulic head changes in individual fractures required to produce these data. The consistency of the method is verified using field experiments performed in a fractured-rock aquifer.     

Hydraulic head applications of flow logs in the study of heterogeneous aquifers

Permeability profiles derived from high-resolution flow logs in heterogeneous aquifers provide a limited sample of the most permeable beds or fractures determining the hydraulic properties of those aquifers. This paper demonstrates that flow logs can also be used to infer the large-scale properties of aquifers surrounding boreholes. The analysis is based on the interpretation of the hydraulic head values estimated from the flow log analysis. Pairs of quasi-steady flow profiles obtained under ambient conditions and while either pumping or injecting are used to estimate the hydraulic head in each water-producing zone. Although the analysis yields localized estimates of transmissivity for a few water-producing zones, the hydraulic head estimates apply to the far-field aquifers to which these zones are connected. The hydraulic head data are combined with information from other sources to identify the large-scale structure of heterogeneous aquifers. More complicated cross-borehole flow experiments are used to characterize the pattern of connection between large-scale aquifer units inferred from the hydraulic head estimates. The interpretation of hydraulic heads in situ under steady and transient conditions is illustrated by several case studies, including an example with heterogeneous permeable beds in an unconsolidated aquifer, and four examples with heterogeneous distributions of bedding planes and/or fractures in bedrock aquifers.     

New Method for Continuous Transmissivity Profiling in Fractured Rock

A new method is presented to search for hydraulically transmissive features in open boreholes in bedrock. A flexible borehole liner made of a watertight, nylon fabric is filled with water to create a constant driving head to evert (reverse of invert) the liner down the hole so that the liner pushes the borehole water out into transmissive fractures or other permeable features. The descent rate is governed by the bulk transmissivity of the remaining permeable features below the liner. Initially, the liner descent rate or velocity is a measure of transmissivity (T) of the entire hole. As the everting liner passes and seals each permeable feature, changes in the liner velocity indicate the position of each feature and an estimate of T using the Thiem equation for steady radial flow. This method has been performed in boreholes with diameters ranging from 96 to 330 mm. Profiling commonly takes a few hours in holes 200‐ to 300‐m long. After arrival of the liner at the bottom of the hole, the liner acts as a seal preventing borehole cross connection between transmissive features at different depths. Liner removal allows the hole to be used for other purposes. The T values determined using this method in a dolostone aquifer were found to be similar to the values from injection tests using conventional straddle packers. This method is not a replacement for straddle‐packer hydraulic testing of specific zones where greater accuracy is desired; however, it is effective and efficient for scanning entire holes for transmissive features.     

Geophysical Assessment of Groundwater Potential: A Case Study from Mian Channu Area,Pakistan

An integrated study using geophysical method in combination with pumping tests and geochemical method was carried out to delineate groundwater potential zones in Mian Channu area of Pakistan. Vertical electrical soundings (VES) using Schlumberger configuration with maximum current electrode spacing (AB/2 = 200 m) were conducted at 50 stations and 10 pumping tests at borehole sites were performed in close proximity to 10 of the VES stations. The aim of this study is to establish a correlation between the hydraulic parameters obtained from geophysical method and pumping tests so that the aquifer potential can be estimated from the geoelectrical surface measurements where no pumping tests exist. The aquifer parameters, namely, transmissivity and hydraulic conductivity were estimated from Dar Zarrouyk parameters by interpreting the layer parameters such as true resistivities and thicknesses. Geoelectrical succession of five‐layer strata (i.e., topsoil, clay, clay sand, sand, and sand gravel) with sand as a dominant lithology was found in the study area. Physicochemical parameters interpreted by World Health Organization and Food and Agriculture Organization were well correlated with the aquifer parameters obtained by geoelectrical method and pumping tests. The aquifer potential zones identified by modeled resistivity, Dar Zarrouk parameters, pumped aquifer parameters, and physicochemical parameters reveal that sand and gravel sand with high values of transmissivity and hydraulic conductivity are highly promising water bearing layers in northwest of the study area. Strong correlation between estimated and pumped aquifer parameters suggest that, in case of sparse well data, geophysical technique is useful to estimate the hydraulic potential of the aquifer with varying lithology.     

Comparison of Methods for Sampling Dissolved Nitrogen in a Fractured Carbonate-Rock Aquifer

As part of an agricultural non-point-source study in the Conestoga River head waters area in Pennsylvania, different methods for collecting ground water samples from a fractured carbonate-rock aquifer were compared. Samples were collected from seven wells that had been cased to bedrock and drilled as open holes to the first significant water-bearing zone. All samples were analyzed for specific conductance, dissolved oxygen, and dissolved-nitrogen species. Water samples collected by a point sampler without pumping the well were compared to samples collected by a submersible pump and by a point sampler after pumping the well. Samples collected by using a point sampler, adjacent to major water-bearing zones in an open borehole without pumping the well, were not statistically different from samples collected from the pump discharge or from point samples collected adjacent to major water-bearing zones after pumping the well. Samples collected by using a point sampler without pumping the well at depths other than those adjacent to the water-bearing zones did not give the same results as the other methods, especially when the water samples were collected from within the well casings. It was concluded that, for the wells at this site, sampling adjacent to major water-bearing zones by using a point sampler without pumping the well provides samples that are as representative of aquifer conditions as samples collected from the pump discharge after reaching constant temperature and specific conductance, and by using a point sampler after pumping the well.     

Neogene aquifer properties specified through the interpretation of electrical sounding data,Salamiyeh region,central Syria

Twenty‐nine Schlumberger electrical soundings were carried out in the Salamiyeh region in Syria using a maximum current electrode separation of 1 km. Three soundings were made at existing boreholes for comparison. Aquifer parameters of hydraulic conductivity and transmissivity were obtained by analysing pumping test data from the existing boreholes. An empirical relationship between hydraulic conductivity determined from the pumping test and both resistivity and thickness of the Neogene aquifer has been established for these boreholes in order to calculate the geophysical hydraulic conductivity. A close agreement has been obtained between the computed hydraulic conductivity and that determined from the pumping test. The relationship established has, therefore, been generalized in the study area in order to evaluate hydraulic conductivity and transmissivity at all the points where geoelectrical measurements have been carried out. This generalization allows one to derive maps of the hydraulic conductivity and transmissivity in the study area based on geoelectrical measurements. These maps are important in future modelling processes oriented towards better exploitation of the aquifers. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

New combined log and tracer test interpretation method for identifying transfers in fissured aquifers

Tracer tests represent the most appropriate approach for assessing hydrodispersive parameters such as transversal and longitudinal dispersivities or kinematic porosity on an aquifer scale. They are generally carried out by injecting a tracer in a borehole and measuring its concentration over time in neighboring boreholes by extracted volume sampling or downhole measurements. Logging is one of the most suitable methods for evaluating fissured reservoirs. But short circuits between fractures with different hydraulic potential through boreholes induce mixing phenomena that cannot be avoided without packers. This mixing can shift the breakthrough curves deduced from the logs for each producing fracture and distort determination of their location.
The method proposed in this paper aims at measuring the flow rate and the solute breakthrough for hydraulically active fractures, in open boreholes. It involves estimating a velocity profile along the borehole column by the analysis of two successive logs: a shift function according to depth is thus determined by comparison between log portions on each successive one. The velocity gradients reflect the inward or outward flow rates produced by each fracture. On the basis of these flow rates, it is possible to determine the mixing effects inside the borehole and then to plot unbiased breakthrough curves for each producing fracture.
This method was applied at a granitic site in the eastern Pyrenees. In spite of some questionable limitations, the results showed that the method seems adapted to situations with many fractures. The precise hydraulic pattern which is obtained at the borehole scale is discussed in terms of a dual porosity model. Furthermore, interpretation of the breakthrough curves for fractures corrected for mixing effects revealed that Peclet numbers are strongly underestimated if this phenomenon is not considered.     

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.     

Biased monitoring of fresh water-salt water mixing zone in coastal aquifers

In coastal aquifers, significant vertical hydraulic gradients are formed where fresh water and underlying salt water discharge together upward to the seafloor. Monitoring boreholes may act as "short circuits" along these vertical gradients, connecting between the higher and the lower hydraulic head zones. When a sea tide is introduced, the fluctuations of both the water table and the depth of the mixing zone are also biased due to this effect. This problem is intensified in places of long-screen monitoring boreholes, which are common in many places in the world. For example, all approximately 500 boreholes of the fresh water-salt water mixing zone in the coastal aquifer of Israel are installed with 10 to 50 m long screens. We present field measurements of these fluctuations, along with a three-dimensional numerical model. We find that the in-well fluctuation magnitude of the mixing zone is an order of magnitude larger than that in the porous media of the actual aquifer. The primary parameters that affect the magnitude of this bias are the anisotropy of the aquifer conductivity and the borehole hydraulic parameters. With no sea tide, borehole interference is higher for the anisotropic case because the vertical hydraulic gradients are high. When tides are introduced, the amplitude of the mixing zone fluctuation is higher for the isotropic case because the overall effective hydraulic conductivity is greater than the conductivity in the anisotropic case. In the aquifer, the fresh water-salt water mixing zone fluctuations are dampened, and tens of meters inland from the shoreline, the fluctuations are on the order of few centimeters.     

A Model to Simulate Pumping Tests in Wells with High Loss and Backflow

Traditional methods of analyzing pumping tests in single wells fail when the well loss is very high due to a low transmissivity skin. Because of the restricted rate at which water can enter a high loss well from the aquifer, well casing storage becomes a significant factor. Additionally, if a slug of water enters the well from the pump column immediately after the pump is switched off, it has a long‐lasting significant effect on the recovering water level in the well because it cannot be absorbed rapidly by the aquifer. A theoretical model is derived here that simulates the water level in a well in these circumstances. In the model, the continuously changing rate of water inflow from the aquifer to the well is approximated by a step function with a finite difference time step. It is demonstrated by a real example that the model can be applied easily to analyze pumping tests, including tests with a varying pumping rate. The analysis confirms suspected high well loss, calculates the unknown rate of backflow, and determines the aquifer's transmissivity.     

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.     

Use of Groundwater Level Fluctuations near an Operating Water Supply Well to Estimate Aquifer Transmissivity

We developed a method to estimate aquifer transmissivity from the hydraulic-head data associated with the normal cyclic operation of a water supply well thus avoiding the need for interrupting the water supply associated with a traditional aquifer test. The method is based on an analytical solution that relates the aquifer's transmissivity to the standard deviation of the hydraulic-head fluctuations in one or more observation wells that are due to the periodic pumping of the production well. We analyzed the resulting analytical solution and demonstrated that when the observation wells are located near the pumping well, the solution has a simple, Dupuit like form. Numerical analysis demonstrates that the analytical solution can also be used for a quasi-periodic pumping of the supply well. Simulation of cyclic pumping in a statistically heterogeneous medium confirms that the method is suitable for analyzing the transmissivity of weakly or moderately heterogeneous aquifers. If only one observation well is available, and the shift in the phase of hydraulic-head oscillations between the pumping well and the observation well is not identifiable. Prior knowledge of aquifer's hydraulic diffusivity is required to obtain the value of the aquifer transmissivity.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

Using Hydrochemical Fades to Delineate Ground Water Flowpaths in Fractured Shale

Large differences in chemistry between sampling points separated In short vertical intervals are often observed in contaminant plumes in both granular and fractured aquifers. However, most regional models assume that such differences will be reduced by dispersive mixing during transport. At a field site located in a discharge area on the Oak Ridge Reservation, Tennessee, ground water flows along discrete flowpaths, as evidenced by the presence of four distinct water types—Ca-HCO₃, Ca-Na-HCO₃, and Na-Ca-HCO₃, and Na-Ca-HCO₃-S0₄—in samples collected from shallow (< 3D in) multilevel wells. The preservation of distinct chemical signatures suggests that ground water must he contained in discrete flow zones during much of its transport time. The chemical composition of the water types can be explained primarily by strata-bound flow over varying flowpath lengths and secondarily by mixing of waters during cross-formational flow in a discharge zone. The hydrochemical facies identified by correlation of water types between the boreholes indicate the general orientation of ground water How paths. These inferred flowpaths are oblique to the orientation of the measured hydraulic gradient and are more closely aligned with bedding and the calculated flow direction. Results of this study indicate that discrete multilevel sampling for analysis of major ions, in addition to information gathered from tracer tests, borehole flow tests. and visual core observations, can provide valuable information on flow directions and preferential flowpaths for contaminant transport.     

Induced Temperature Gradients to Examine Groundwater Flowpaths in Open Boreholes

Techniques for characterizing the hydraulic properties and groundwater flow processes of aquifers are essential to design hydrogeologic conceptual models. In this study, rapid time series temperature profiles within open‐groundwater wells in fractured rock were measured using fiber optic distributed temperature sensing (FO‐DTS). To identify zones of active groundwater flow, two continuous electrical heating cables were installed alongside a FO‐DTS cable to heat the column of water within the well and to create a temperature difference between the ambient temperature of the groundwater in the aquifer and that within the well. Additional tests were performed to examine the effects of pumping on hydraulic fracture interconnectivity around the well and to identify zones of increased groundwater flow. High‐ and low‐resolution FO‐DTS cable configurations were examined to test the sensitivities of the technique and compared with downhole video footage and geophysical logging to confirm the zones of active groundwater flow. Two examples are presented to demonstrate the usefulness of this new technique for rapid characterization of fracture zones in open boreholes. The combination of the FO‐DTS and heating cable has excellent scope as a rapid appraisal tool for borehole construction design and improving hydrogeologic conceptual models.     

Fracture Flow Characterization with Low-Noise Spontaneous Potential Logging

Geophysical well logging has been applied for fracture characterization in crystalline terrains by physical properties measurements and borehole wall imaging. Some of these methods can be applied to monitor pumping tests to identify fractures contributing to groundwater flow and, with this, determine hydraulic conductivity and transmissivity along the well. We present a procedure to identify fractures contributing to groundwater flow using spontaneous potential measurements generated by electrokinetic processes when the borehole water head is lowered and then monitored while recovering. The electrokinetic model for flow through a tabular gap is used to interpret the measured data and determine the water head difference that drives the flow through the fracture. We present preliminary results at a test site in crystalline rocks on the campus of the University of São Paulo.     

Flow in an aquifer charged with hot water from a fault zone

Many geothermal anomalies are intersected by vertical fault zones (narrow zones of fractured material with large effective permeability). These conduits are probably responsible for much of the upwelling of hot water from depth. This paper considers a shallow aquifer intersected by a vertical fault. The fluid flow in the aquifer is numerically modeled as a two-dimensional problem. It is observed that the temperature distribution in the aquifer is governed primarily by lateral flow of hot water supplied from the intersecting vertical fault and only secondarily by conduction. The numerical results also provide a possible explanation for the local temperature maxima and inversions occasionally observed in borehole measurements. The present model is an alternative to that based on mushroom-shaped isotherm distributions found in high Rayleigh number large-scale circulation cell calculations.