Hydrogeophysical Methods for Analyzing Aquifer Storage and Recovery Systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ～500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.     

Vertical Contaminant Profiling of Volatile Organic* in a Deep Fractured Basalt Aquifer

Volatile organic compounds delected in ground water from wells at Test Area North (TAN) at the Idaho National Engineering Laboratory (INEL) prompted RCRA facility investigations in 1989 and 1990 and a CERCLA-driven RI/FS in 1992. In order to address ground water treatment feasibility, one of the main objectives, of the 1992 remedial investigation was to determine the vertical extent of ground water contamination, where the principle contaminant, of concern is trichloroethylene (TCE). It was hypothesized that a sedimentary interbed at depth in the fractured basalt aquifer could be inhibiting vertical migration of contaminants to lower aquifers. Due to the high cost of drilling and installation of ground water monitoring wells at this facility (greater than $100,000 per well), a real time method was proposed for obtaining and analyzing ground water samples during drilling to allow accurate placement of well screens in zones of predicted VOC contamination. This method utilized an inflatable pump packer pressure transducer system interfaced with a datalogger and PC at land surface. This arrangement allowed for real lime monitoring of hydraulic head above and below the packer to detect leakage around the packer during pumping and enabled collection of head data during pumping for estimating hydrologic properties. Analytical results were obtained in about an hour from an on-site mobile laboratory equipped with a gas chromalograplvmass spectrometer (GC/MS). With the hydrologic and analytical results in hand, a decision was made to either complete the well or continue drilling to the next test zone. In almost every case, analytical results of ground water samples taken from the newly installed wells closely replicated the water quality of ground water samples obtained through the pump packer system.     

Optimizing Multiwell Aquifer Storage and Recovery Systems for Energy Use and Recovery Efficiency

As more aquifer storage and recovery (ASR) systems are employed for management of water resources, the skillful operation of multiwell ASR systems has become very important to improve their performance. In this study, we developed MODFLOW and MT3DMS models to simulate a multiwell ASR system in a synthetic aquifer to assess effects of hydrogeological and operational factors on the performance of the multiwell ASR system. We evaluated a simplified (dual well) ASR system in comparison with complex system (three-, four-, five-, and seven-well systems). Recovery and energy efficiencies were calculated using the model simulations. Factors such as higher hydraulic conductivity and longitudinal dispersivity significantly reduced the recovery and energy efficiencies of the system. In contrast, increasing the volume of recharged water increased the recovery efficiency; however, the energy efficiency was reduced. Recovery and energy efficiencies also plummet when there is an increase in the underlying regional gradient and the designed storage duration. Operating the system multiple times can yield higher volume of potable water, but the energy efficiency may not vary significantly after the second operating cycle. Single-well systems and multiwell systems exhibit similar responses to changes in physical factors, although operational factors have a more pronounced effect on the multiwell systems. One of the major findings was that fewer wells in a multiwell ASR system can yield higher volume of potable water and better output with respect to the electrical power being consumed. The results provide design engineers with guidelines for optimizing performance of the multiwell ASR systems.     

Arsenic Control During Aquifer Storage Recovery Cycle Tests in the Floridan Aquifer

Implementation of aquifer storage recovery (ASR) for water resource management in Florida is impeded by arsenic mobilization. Arsenic, released by pyrite oxidation during the recharge phase, sometimes results in groundwater concentrations that exceed the 10 µg/L criterion defined in the Safe Drinking Water Act. ASR was proposed as a major storage component for the Comprehensive Everglades Restoration Plan (CERP), in which excess surface water is stored during the wet season, and then distributed during the dry season for ecosystem restoration. To evaluate ASR system performance for CERP goals, three cycle tests were conducted, with extensive water‐quality monitoring in the Upper Floridan Aquifer (UFA) at the Kissimmee River ASR (KRASR) pilot system. During each cycle test, redox evolution from sub‐oxic to sulfate‐reducing conditions occurs in the UFA storage zone, as indicated by decreasing Fe²⁺/H₂S mass ratios. Arsenic, released by pyrite oxidation during recharge, is sequestered during storage and recovery by co‐precipitation with iron sulfide. Mineral saturation indices indicate that amorphous iron oxide (a sorption surface for arsenic) is stable only during oxic and sub‐oxic conditions of the recharge phase, but iron sulfide (which co‐precipitates arsenic) is stable during the sulfate‐reducing conditions of the storage and recovery phases. Resultant arsenic concentrations in recovered water are below the 10 µg/L regulatory criterion during cycle tests 2 and 3. The arsenic sequestration process is appropriate for other ASR systems that recharge treated surface water into a sulfate‐reducing aquifer.     

Aquifer Storage and Recovery Using Saline Aquifers: Hydrogeological Controls and Opportunities

Aquifer storage and recovery (ASR) is a valuable tool for managing variations in the supply and demand of freshwater, but system performance is highly dependent upon system-specific hydrogeological conditions including the salinity of the storage-zone native groundwater. ASR systems using storage zones containing saline (>10,000 mg/L of total dissolved solids) groundwater tend to have relatively low recovery efficiencies (REs). However, the drawbacks of low REs may be offset by lesser treatment requirements and may be of secondary importance where the stored water (e.g., excess reclaimed, surface, and storm waters) would otherwise go to waste and pose disposal costs. Density-dependent, solute-transport modeling results demonstrate that the RE of ASR systems using a saline storage zone is most strongly controlled by parameters controlling free convection (e.g., horizontal hydraulic conductivity) and mixing of recharged and native groundwater (e.g., dispersivity and aquifer heterogeneity). Preferred storage zone conditions are moderate hydraulic conductivities (5 to 20 m/d), low degrees of aquifer heterogeneity, and primary porosity-dominated siliclastic and limestones lithologies with effective porosities greater than 5%. Where hydrogeological conditions are less favorable, operational options are available to improve RE, such as preferential recovery from the top of the storage zone. Injection of large volumes of excess water currently not needed into saline aquifers could create valuable water resources that could be tapped in the future during times of greater need.     

Relative Recovery of Thermal Energy and Fresh Water in Aquifer Storage and Recovery Systems

This paper explores the relationship between thermal energy and fresh water recoveries from an aquifer storage recovery (ASR) well in a brackish confined aquifer. It reveals the spatial and temporal distributions of temperature and conservative solutes between injected and recovered water. The evaluation is based on a review of processes affecting heat and solute transport in a homogeneous aquifer. In this simplified analysis, it is assumed that the aquifer is sufficiently anisotropic to inhibit density‐affected flow, flow is axisymmetric, and the analysis is limited to a single ASR cycle. Results show that the radial extent of fresh water at the end of injection is greater than that of the temperature change due to the heating or cooling of the geological matrix as well as the interstitial water. While solutes progress only marginally into low permeability aquitards by diffusion, conduction of heat into aquitards above and below is more substantial. Consequently, the heat recovery is less than the solute recovery when the volume of the recovered water is lower than the injection volume. When the full volume of injected water is recovered the temperature mixing ratio divided by the solute mixing ratio for recovered water ranges from 0.95 to 0.6 for ratios of maximum plume radius to aquifer thickness of 0.6 to 4.6. This work is intended to assist conceptual design for dual use of ASR for conjunctive storage of water and thermal energy to maximize the potential benefits.     

Storage and Drainage Characteristics of a Highly Heterogeneous Karst Aquifer in Houzhai Basin

The characteristics of karst aquifers are difficult to be determined due to their heterogeneous physical properties and lack of hydrogeological information. In this case study, we applied two methods for a comparative analysis of storage and drainage characteristics in upstream, midstream, and downstream of Houzhai cave stream basin. In the first method, Minimum Smoothed Method (MSM) is used to determine the proportion of baseflow to the total flow (Baseflow Index, BFI). In the second method, a bicarbonate‐base two‐end member mixing model is used to quantify the slow flow component and fast flow component. For both methods, slow flow and quick flow are quantified at three sampling sites, which provide useful information for the analysis of storage and drainage characteristics. The results from flow separation method and hydrogeochemical analysis show a consistently increasing trend of the proportion of slow flow to total flow from the upstream to downstream which indicates that the voids of highly conductive conduits and well‐connected fissures decrease along the flow paths in the Houzhai cave stream basin in southwest China. The upstream areas have a low proportion of baseflow which indicates a high drainage capacity due to high permeable conduits and well‐connected fissures. The downstream areas, on the contrary, have a high proportion of baseflow which indicates a high storage capacity and slow infiltration due to the predominant presence of matrix and poorly‐connected fissures. These numerical methods provide alternative ways to investigate the storage and drainage characteristics of karst aquifers where direct measurement are not available.     

Recovery of Injected Freshwater from a Brackish Aquifer with a Multiwell System

Herein we propose a multiple injection and recovery well system strategically operated for freshwater storage in a brackish aquifer. With the system we call aquifer storage transfer and recovery (ASTR) by using four injection and two production wells, we are capable of achieving both high recovery efficiency of injected freshwater and attenuation of contaminants through adequately long residence times and travel distances within the aquifer. The usual aquifer storage and recovery (ASR) scheme, in which a single well is used for injection and recovery, does not warrant consistent treatment of injected water due to the shorter minimum residence times and travel distances. We tested the design and operation of the system over 3 years in a layered heterogeneous limestone aquifer in Salisbury, South Australia. We demonstrate how a combination of detailed aquifer characterization and solute transport modeling can be used to maintain acceptable salinity of recovered water for its intended use along with natural treatment of recharge water. ASTR can be used to reduce treatment costs and take advantage of aquifers with impaired water quality that might locally not be otherwise beneficially used.     

Cyclodextrin-Enhanced Solubilization of Organic Contaminants with Implications for Aquifer Remediation

Reagents that enhance the aqueous solubility of non-aqueous phase organic liquid (NAPL) contaminants are under investigation for use in enhanced subsurface remediation technologies. Cyclodextrin, a glucose-based molecule, is such a reagent. In this paper, laboratory experiments and numerical model simulations are used to evaluate and understand the potential remediation performance of cyclodextrin. Physical properties of cyclodextrin solutions such as density, viscosity, and NAPL-aqueous inter-facial tension are measured. Our analysis indicates that no serious obstacles exist related to fluid properties that would prevent the use of cyclodextrin solutions for subsurface NAPL remediation. Cyclodextrin-enhanced solubilization for a large suite of typical ground water contaminants is measured in the laboratory, and the results are related to the physicochemical properties of the organic compounds. The most-hydrophobic contaminants experience a larger relative solubility enhancement than the less-hydrophobic contaminants but have lower aqueous-phase apparent solubilities. Numerical model simulations of enhanced-solubilization flushing of NAPL-contaminated soil demonstrate that the more-hydrophilic compounds exhibit the greatest mass-removal rates due to their greater apparent solubilities, and thus are initially more effectively removed from soil by enhanced-solubilization-flushing reagents. However, the relatively more hydrophobic contaminants exhibit a greater improvement in contaminant mass-removal (compared with water flushing) than that exhibited for the relatively hydrophilic contaminants.     

222Rn as Natural Tracer for LNAPL Recovery in a Crude Oil‐Contaminated Aquifer

The objective of this study was to investigate whether ²²²Rn in groundwater can be used as a tracer for light non‐aqueous phase liquid (LNAPL) quantification at a field site treated by dual‐phase LNAPL removal. After the break of a pipeline, 5 ha of soil in the nature reserve Coussouls de Crau in southern France was contaminated by 5100 m³ of crude oil. Part of this oil seeped into the underlying gravel aquifer and formed a floating oil body of about 3.9 ha. The remediation consists of plume management by hydraulic groundwater barriers and LNAPL extraction in the source zone. ²²²Rn measurements were performed in 21 wells in and outside the source zone during 15 months. In uncontaminated groundwater, the radon activity was relatively constant and remained always >11 Bq/L. The variability of radon activity measurements in wells affected by the pump‐and‐skim system was consistent with the measurements in wells that were not impacted by the system. The mean activities in wells in the source zone were, in general, significantly lower than in wells upgradient of the source zone, owing to partitioning of ²²²Rn into the oil phase. The lowest activities were found in zones with high non‐aqueous phase liquid (NAPL) recovery. LNAPL saturations around each recovery well were furthermore calculated during a period of high groundwater level, using a laboratory‐determined crude oil–water partitioning coefficient of 38.5 ± 2.9. This yielded an estimated volume of residual crude oil of 309 ± 93 m³ below the capillary fringe. We find that ²²²Rn is a useful and cheap groundwater tracer for finding zones of good LNAPL recovery in an aquifer treated by dual‐phase LNAPL removal, but that quantification of NAPL saturation using Rn is highly uncertain.