Electrical Resistivity Imaging of a Permanganate Injection During In Situ Treatment of RDX‐Contaminated Groundwater

Groundwater beneath the former Nebraska Ordnance Plant (NOP) is contaminated with the explosive hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) and trichloroethene (TCE). Previous treatability experiments confirmed that permanganate could mineralize RDX in NOP aquifer material. The objective of this study was to determine the efficacy of permanganate to transform RDX in the field by monitoring a pilot‐scale in situ chemical oxidation (ISCO) demonstration. In this demonstration, electrical resistivity imaging (ERI) was used to create two‐dimensional (2‐D) images of the test site prior to, during, and after injecting sodium permanganate. The ISCO was performed by using an extraction‐injection well configuration to create a curtain of permanganate. Monitoring wells were positioned downgradient of the injection zone with the intent of capturing the permanganate‐RDX plume. Differencing between ERI taken preinjection and postinjection determined the initial distribution of the injected permanganate. ERI also quantitatively corroborated the hydraulic conductivity distribution across the site. Groundwater samples from 12 downgradient wells and 8 direct‐push profiles did not provide enough data to quantify the distribution and flow of the injected permanganate. ERI, however, showed that the permanganate injection flowed against the regional groundwater gradient and migrated below monitoring well screens. ERI combined with monitoring well samples helped explain the permanganate dynamics in downgradient wells and support the use of ERI as a means of monitoring ISCO injections.     

Pump-and-treat optimization using analytic element method flow models

Plume containment using pump-and-treat (PAT) technology continues to be a popular remediation technique for sites with extensive groundwater contamination. As such, optimization of PAT systems, where cost is minimized subject to various remediation constraints, is the focus of an important and growing body of research. While previous pump-and-treat optimization (PATO) studies have used discretized (finite element or finite difference) flow models, the present study examines the use of analytic element method (AEM) flow models. In a series of numerical experiments, two PATO problems adapted from the literature are optimized using a multi-algorithmic optimization software package coupled with an AEM flow model. The experiments apply several different optimization algorithms and explore the use of various pump-and-treat cost and constraint formulations. The results demonstrate that AEM models can be used to optimize the number, locations and pumping rates of wells in a pump-and-treat containment system. Furthermore, the results illustrate that a total outflux constraint placed along the plume boundary can be used to enforce plume containment. Such constraints are shown to be efficient and reliable alternatives to conventional particle tracking and gradient control techniques. Finally, the particle swarm optimization (PSO) technique is identified as an effective algorithm for solving pump-and-treat optimization problems. A parallel version of the PSO algorithm is shown to have linear speedup, suggesting that the algorithm is suitable for application to problems that are computationally demanding and involve large numbers of wells.     

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.     

Microtopographical and hydrophysical controls on subsurface flow and solute transport: A continuous solute release experiment in a subarctic bog

Resource extraction and transportation activities in subarctic Canada can result in the unintentional release of contaminants into the surrounding peatlands. In the event of a release, a thorough understanding of solute transport within the saturated zone is necessary to predict plume fate and the potential impacts on peatland ecosystems. To better characterize contaminant transport in these systems, approximately 13,000 L/day of sodium chloride tracer (200 mg/L) was released into a bog in the James Bay Lowland. The tracer was pumped into a fully penetrating well (1.5 m) between July 5 and August 18, 2015. Horizontal and vertical plume development was measured via in situ specific conductance and water table depth from an adaptive monitoring network. Over the spill period, the bulk of the plume travelled a lateral distance of 100 m in the direction of the slight regional groundwater and topographical slope. The plume shape was irregular and followed the hollows, indicating preferential flow paths due to the site microtopography. Saturated transport of the tracer occurred primarily at ~25 cm below ground surface (bgs), and at a discontinuous high hydraulic conductivity layer ~125 cm bgs due to a complex and heterogeneous vertical hydraulic conductivity profile. Plume measurement was confounded by a large amount of precipitation (233 mm over the study period) that temporarily diluted the tracer in the highly conductive upper peat layer. Longitudinal solute advection can be approximated using local water table information (i.e., depth and gradient); microtopography; and meteorological conditions. Vertical distribution of solute within the peat profile is far more complex due to the heterogeneous subsurface; characterization would be aided by a detailed understanding of the site‐specific peat profile; the degree of decomposition; and the type of contaminant (e.g., reactive/nonreactive). The results of this research highlight the difficulty of tracking a contaminant spill in bogs and provide a benchmark for the characterization of the short‐term fate of a plume in these complex systems.     

Remediation of Heterogeneous Aquifers Subject to Uncertainty

Optimal cost pump-and-treat ground water remediation designs for containment of a contaminated aquifer are often developed using deterministic ground water models to predict ground water flow. Uncertainty in hydraulic conductivity fields used in these models results in remediation designs that are unreliable. The degree to which uncertainty contributes to the reliability of remediation designs as measured by the characterization of the uncertainty is shown to differ depending upon the geologic environments of the models. This conclusion is drawn from the optimal design costs for multiple deterministic models generated to represent the uncertainty of four distinct models with different geologic environments. A multi scenario approach that includes uncertainty into the remediation design called the deterministic method for optimization subject to uncertainty (DMOU) is applied to these distinct models. It is found that the DMOU is a method for determining a remediation design subject to uncertainty that requires minimal postprocessing efforts. Preprocessing, however, is required for the application of the DMOU to unique problems. In the ground water remediation design problems, the orientation of geologic facies with respect to the orientation of flow patterns, pumping well locations, and constraint locations are shown to affect the preprocessing, the solutions to the DMOU problems, and the computational efficiency of the DMOU approach. The results of the DMOU are compared to the results of a statistical analysis of the effects of the uncertainty on remediation designs. This comparison validates the efficacy of the DMOU and illustrates the computational advantages to using the DMOU over statistical measures.     

Imaging hyporheic zone solute transport using electrical resistivity

Traditional characterization of hyporheic processes relies upon modelling observed in‐stream and subsurface breakthrough curves to estimate hyporheic zone size and infer exchange rates. Solute data integrate upstream behaviour and lack spatial coverage, limiting our ability to accurately quantify spatially heterogeneous exchange dynamics. Here, we demonstrate the application of near‐surface electrical resistivity imaging (ERI) methods, coupled with experiments using an electrically conductive stream tracer (dissolved NaCl), to provide in situ imaging of spatial and temporal dynamics of hyporheic exchange. Tracer‐labelled water in the stream enters the hyporheic zone, reducing electrical resistivity in the subsurface (to which subsurface ERI is sensitive). Comparison of background measurements with those recording tracer presence provides distributed characterization of hyporheic area (in this application, ∼0·5 m²). Results demonstrate the first application of ERI for two‐dimensional imaging of stream‐aquifer exchange and hyporheic extent. Future application of this technique will greatly enhance our ability to quantify processes controlling solute transport and fate in hyporheic zones, and provide data necessary to inform more complete numerical models. Copyright © 2010 John Wiley & Sons, Ltd.     

In Situ Chemical Oxidation of RDX-Contaminated Groundwater with Permanganate at the Nebraska Ordnance Plant

Groundwater beneath the former Nebraska Ordnance Plant (NOP) is contaminated with the explosive hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). The current pump and treat facility is preventing offsite migration but does not offer a short-term solution. Our objective was to quantify the effectiveness of permanganate to degrade RDX in situ. This was accomplished by performing laboratory treatability experiments, aquifer characterization, and a pilot-scale in situ chemical oxidation (ISCO) demonstration. Treatability experiments confirmed that permanganate could mineralize RDX in the presence of NOP aquifer solids. The pilot-scale ISCO demonstration was performed using an extraction-injection well configuration to create a curtain of permanganate between two injection wells. RDX destruction was then quantified as the RDX-permanganate plume migrated downgradient through a monitoring well field. Electrical resistivity imaging (ERI) was used to identify the subsurface distribution of permanganate after injection. Results showed that RDX concentrations temporally decreased in wells closest to the injection wells by 70% to 80%. Observed degradation rates (0.12 and 0.087/d) were lower than those observed under laboratory batch conditions at 11.5 °C (0.20/d) and resulted from lower than projected permanganate concentrations. Both ERI and spatial electrical conductivity measurements verified that permanganate distribution was not uniform throughout the 6.1-m (20 feet) well screens and that groundwater sampling captured both treated and nontreated groundwater during pumping. Although heterogeneous flow paths precluded a uniform permanganate distribution, pilot-scale results provided proof-of-concept that permanganate can degrade RDX in situ and support permanganate as a possible remedial treatment for RDX-contaminated groundwater.     

Geophysical mapping of hyporheic processes controlled by sedimentary architecture within compound bar deposits

Hyporheic exchange influences water quality and controls numerous physical, chemical, and biological processes. Despite its importance, hyporheic exchange and the associated dynamics of solute mixing are often difficult to characterize due to spatial (e.g., sedimentary heterogeneity) and temporal (e.g., river stage fluctuation) variabilities. This study coupled geophysical techniques with physical and chemical sediment analyses to map sedimentary architecture and quantify its influence on hyporheic exchange dynamics within a compound bar deposit in a gravel-dominated river system in southwestern Ohio. Electromagnetic induction (EMI) was used to quantify variability in electrical conductivity within the compound bar. EMI informed locations of electrode placement for time-lapse electrical resistivity imaging (ERI) surveys, which were used to examine changes in electrical resistivity driven by hyporheic exchange. Both geophysical methods revealed a zone of high electrical conductivity in the center of the bar, identified as a fine-grained cross-bar channel fill. The zone acts as a baffle to flow, evidenced by stable electrical conditions measured by time-lapse ERI over the study period. Large changes in electrical resistivity throughout the survey period indicate preferential flowpaths through higher permeability sands and gravels. Grain size analyses confirmed sedimentological interpretations of geophysical data. Loss on ignition and x-ray fluorescence identified zones with higher organic matter content that are locations for potentially enhanced geochemical activity within the cross-bar channel fill. Differences in the physical and geochemical characteristics of cross-bar channel fills play an important role in hyporheic flow dynamics and nutrient processing within riverbed sediments. These findings enhance our understanding of the applications of geophysical methods in mapping riverbed heterogeneity and highlight the importance of accurately representing geomorphologic features and heterogeneity when studying hyporheic exchange processes.     

Tracking CO2 plume migration during geologic sequestration using a probabilistic history matching approach

Tracking the migration of the CO₂ plume is essential in order to better manage the operation of geologic sequestration of CO₂. However, the large cost of most monitoring technologies, such as time-lapse seismic, limits their application. We investigated the application of a probabilistic history matching methodology using routine measurements of injection data, which are affected by the presence of large-scale heterogeneities, as an inexpensive alternative to track the migration of CO₂ plume in an aquifer. The approach is demonstrated first through a synthetic example in which the ability to detect the presence of flow barriers was tested. In a second example, we applied our method to the In Salah field, one of the largest geological sequestration projects in the world, where the main direction of high permeability features was inferred. The accuracy and reproducibility of the results show that our approach for assisted history matching is an economic and viable option for plume monitoring during geologic CO₂ sequestration.     

Natural Gradient Tracer Test to Evaluate Natural Attenuation of MTBE Under Anaerobic Conditions

A natural gradient tracer test using perdeuterated MTBE was conducted in an anaerobic aquifer to determine the relative importance of dispersion and degradation in reducing MTBE concentrations in ground water. Preliminary ground water chemistry and hydraulic conductivity data were used to place the tracer within an existing dissolved MTBE plume at Port Hueneme, California. Following one year of transport, the tracer plume was characterized in detail.
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m²/day and 0.08 m²/day, respectively). A mass-balance analysis, carried out after one year of transport, accounted for 110% of the injected mass and indicated that no significant mass loss occurred. The plume structure created by zones of higher and lower hydraulic conductivity at the site was complex, consisting of several localized areas of high tracer concentration in a lower concentration plume. This is important because the aquifer has generally been characterized as exhibiting fairly minor heterogeneity. In addition, the tracer plume followed a curved flowpath that deviated from the more macroscopic direction of ground water flow inferred from local ground water elevation measurements and the behavior of the existing plume. Understanding the mass balance, plume structure, curvature of the tracer plume, and consequently natural attenuation behavior required the detailed sampling approach employed in this study. These data imply that a detailed understanding of site hydrogeology and an extensive sampling network may be critical for the correct interpretation of monitored natural attenuation of MTBE.     

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.     

Electrokinetic Migration of Nitrate Through Heterogeneous Granular Porous Media

This study investigates and quantifies the influence of physical heterogeneity in granular porous media, represented by materials with different hydraulic conductivity, on the migration of nitrate, used as an amendment to enhance bioremediation, under an electric field. Laboratory experiments were conducted in a bench‐scale test cell under a low applied direct current using glass bead and clay mixes and synthetic groundwater to represent ideal conditions. The experiments included bromide tracer tests in homogeneous settings to deduce controls on electrokinetic transport of inorganic solutes in the different materials, and comparison of nitrate migration under homogeneous and heterogeneous scenarios. The results indicate that physical heterogeneity of subsurface materials, represented by a contrast between a higher‐hydraulic conductivity and lower‐hydraulic conductivity material normal to the direction of the applied electric field exerts the following controls on nitrate migration: (1) a spatial change in nitrate migration rate due to changes in effective ionic mobility and subsequent accumulation of nitrate at the interface between these materials; and (2) a spatial change in the voltage gradient distribution across the hydraulic conductivity contrast, due to the inverse relationship with effective ionic mobility. These factors will contribute to higher mass transport of nitrate through low hydraulic conductivity zones in heterogeneous porous media, relative to homogeneous host materials. Overall electrokinetic migration of amendments such as nitrate can be increased in heterogeneous granular porous media to enhance the in situ bioremediation of organic contaminants present in low hydraulic conductivity zones.     

Modeling a Ground Water Circulation Well Alternative

Ground water circulation wells (GCWs) provide an appealing alternative to typical pump-and-treat ground water remediation systems because of the inherent resource-conservative nature of the GCW systems. GCW performance prediction is challenging because the consideration of extraction and recharge in a single well is unusual for most practitioners, the technology is relatively new, and a meaningful body of literature has not been published. A three-part evaluation process using state-of-the-practice numerical ground water flow and mass transport models was developed for application during GCW pilot studies at the former Nebraska Ordnance Plant site. A small-scale ground water flow model was developed during the pilot study planning process to predict the system performance and to locate performance-measuring monitoring wells. Key predictions included the capture zone predicted to develop upgradient of the GCW, the downgradierit GCW recharge zone, and the circulation zone centered on the GCW. The flow model was subsequently verified using ground water elevation data and contaminant concentration data collected during pilot study operation. Aquifer parameters were reestimated as a result of the verification process. Those parameter values were used as input to a larger scale model, which was used to develop a remedial alternative consisting of multiple GCW systems.     

Encapsulated Cell Bioremediation: Evaluation on the Basis of Particle Tracer Tests

Microencapsulation of degradative organisms enhances microorganism survivability (Stormo and Crawford 1994). The use of encapsulated cell microbeads for in situ biodegradation depends not only on microorganism survival but also on microbead transport characteristics. Two forced-gradient, recirculating-loop tracer experiments were conducted to evaluate the feasibility of encapsulated cell transport and bioremediation on the basis of polystyrene microsphere transport results. The tracer tests were conducted in a shallow, confined, unconsolidated, heterogeneous, sedimentary aquifer using bromide ion and 2 μm, 5 μn, and 15 μm microsphere tracers. Significant differences were observed in the transport of bromide solute and polystyrene microspheres. Microspheres reached peak concentrations in monitoring wells before bromide, which was thought to reflect the influence of aquifer heterogeneity. Greater decreases in microsphere C/C_o ratios were observed with distance from the injection wells than in bromide C/C_o ratios, which was attributed to particle filtration and/or settling. Several methods might be considered for introducing encapsulated cell microbeads into a subsurface environment, including direct injection into a contaminated aquifer zone, injection through a recirculating ground water flow system, or emplacement in a subsurface microbial curtain in advance of a plume. However, the in situ use of encapsulated cells in an aquifer is probably limited to aquifers containing sufficiently large pore spaces, allowing passage of at least some encapsulated cells. The use of encapsulated cells may also be limited by differences in solute and microbead transport patterns and flowpath clogging by larger encapsulated cell microbeads.     

Imaging Pathways in Fractured Rock Using Three‐Dimensional Electrical Resistivity Tomography

Major challenges exist in delineating bedrock fracture zones because these cause abrupt changes in geological and hydrogeological properties over small distances. Borehole observations cannot sufficiently capture heterogeneity in these systems. Geophysical techniques offer the potential to image properties and processes in between boreholes. We used three‐dimensional cross borehole electrical resistivity tomography (ERT) in a 9 m (diameter) × 15 m well field to capture high‐resolution flow and transport processes in a fractured mudstone contaminated by chlorinated solvents, primarily trichloroethylene. Conductive (sodium bromide) and resistive (deionized water) injections were monitored in seven boreholes. Electrode arrays with isolation packers and fluid sampling ports were designed to enable acquisition of ERT measurements during pulsed tracer injections. Fracture zone locations and hydraulic pathways inferred from hydraulic head drawdown data were compared with electrical conductivity distributions from ERT measurements. Static ERT imaging has limited resolution to decipher individual fractures; however, these images showed alternating conductive and resistive zones, consistent with alternating laminated and massive mudstone units at the site. Tracer evolution and migration was clearly revealed in time‐lapse ERT images and supported by in situ borehole vertical apparent conductivity profiles collected during the pulsed tracer test. While water samples provided important local information at the extraction borehole, ERT delineated tracer migration over spatial scales capturing the primary hydrogeological heterogeneity controlling flow and transport. The fate of these tracer injections at this scale could not have been quantified using borehole logging and/or borehole sampling methods alone.     

A Method for Estimating Field-Scale Mass Transfer Rate Parameters and Assessing Aquifer Cleanup Times

A general method for estimating ground-water solute mass transfer rate parameters from field test data is presented. The method entails matching solute concentration and hydraulic head data collected during the recovery phase of a pumping test through application of a simulation-regression technique. Estimation of hydraulic conductivity and mass transfer rate parameter values is performed by fitting model simulations to the data. Parameter estimates are utilized to assess cleanup times for pump-and-treat aquifer remediation scenarios. Uncertainty in the cleanup time estimate is evaluated using statistical information obtained with the parameter estimation technique. Application of the method is demonstrated using a hypothetical ground-water flow and solute transport system. Simulations of field testing, parameter estimation, and remedial time frames are performed to evaluate the usefulness of the method. Sets of random noise that signify potential field and laboratory measurement errors are combined with the hypothetical data to provide rigorous testing of the method. Field tests are simulated using ranges of values for data noise, the mass transfer rate parameters, the test pumping rates, and the duration of recovery monitoring to evaluate their respective influence on parameter and cleanup time estimates. The demonstration indicates the method is capable of yielding accurate estimates of the solute mass transfer rate parameters. When the parameter values for the hypothetical system are well estimated, cleanup time predictions are shown to be more accurate than when calculated using the local equilibrium assumption.     

Ground Water Sampling Bias Observed in Shallow, Conventional Wells

A previous field demonstration project on nitrate-based bioremediation of a fuel-contaminated aquifer used short-screened clustered well points in addition to shallow (10 foot), conventional monitoring wells to monitor the progress of remediation during surface application of recharge. These well systems were placed in the center and at one edge of each of two treatment cells. One cell received recharge amended with nitrate (nitrate cell), and the other received unamended recharge (control cell). Data from the clustered well points were averaged to provide a mean estimate for comparison with the associated conventional monitoring well.
Conservative tracer profiles were similar for each of the four systems, with better fits obtained for well systems located at the edge of the treatment cells. However, aromatic hydrocarbon and electron acceptor profiles varied greatly for the two center well systems, with the conventional monitoring well data suggesting that remediation was proceeding at a much more rapid rate than indicated by the cluster well points. Later tests with an electromagnetic borehole flowmeter demonstrated a significant vertical flow through the well-bore of the conventional monitoring well under simulated operating conditions. This created an artifact during sampling, thought to arise from preferential flow of recharge water from the water table to deeper portions of the contaminated zone resulting in several effects, including an actual decreased residence time of water sampled by the conventional well. These data provide additional evidence that conventional monitoring wells may be inadequate for monitoring remediation in the presence of significant vertical hydraulic gradients, even for fairly shallow homogeneous aquifers.     

Polymer-Enhanced DNAPL Flushing from Low-Permeability Media: An Experimental Study

Remediation of dense nonaqueous phase liquids (DNAPLs) is recognized as one of the most difficult problems associated with ground water pollution. The pump-and-treat technique, usually consisting of a continuous operation of extraction-injection wells, is widely used for ground water remediation. In a stratified or otherwise heterogeneous aquifer, however, this technique suffers from tailing and rebound problems, which limit its cleanup efficiency and result in higher operation costs. The tailing and rebound is usually due to slow diffusion of contaminants out of lower- permeability heterogeneities into the flow regime of the higher-permeability zone. In this study, we conduct bench-scale experiments to investigate a novel polymer system and injection method to improve the pump-and-treat technique for DNAPL trapped in a layer of porous media that has a relatively low permeability compared to the surrounding media. This technique might be useful, for example, to remove DNAPL from these low-permeability zones after removal of DNAPL from the higher-permeability zones by a more traditional remediation method. The polymer system consists of a mixture of anionic and cationic polyacrylamides in solution and the injection method is based on flow-induced polymer adsorption, called bridging adsorption. The study includes single and parallel-column experiments. The measured polymer penetration depths were compared with values predicted from a numerical simulation, which was developed previously by the authors of this paper. The experiments and simulations show that the polymer injection leads to a modification of the permeability contrast that favors a more efficient pump-and-treat process. These results suggest that additional research to upscale the technology to pilot scales is warranted.     

A 3D ERT study of solute transport in a large experimental tank

A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m³. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.