Modeling unsaturated flow in a layered formation under quasi-steady state conditions using geophysical data constraints

The identification of vadose zone flow parameters and solute travel time from the surface to the water table are key issues for the assessment of groundwater vulnerability. In this paper we use the results of time-lapse monitoring of the vadose zone in a UK consolidated sandstone aquifer using cross-hole zero-offset radar to assess and calibrate models of water flow in the vadose zone. The site under investigation is characterized by a layered structure, with permeable medium sandstone intercalated by finer, less permeable, laminated sandstone. Information on this structure is available from borehole geophysical (gamma-ray) logs. Monthly cross-hole radar monitoring was performed from August 1999 to February 2001, and shows small changes of moisture content over time and fairly large spatial variability with depth. One-dimensional Richards’ equation modeling of the infiltration process was performed under spatially heterogeneous, steady state conditions. Both layer structure and Richards’ equation parameters were simulated using a nested Monte Carlo approach, constrained via geostatistical analysis on the gamma-ray logs and on a priori information regarding the possible range of hydraulic parameters. The results of the Monte Carlo analysis show that, in order to match the radar-derived moisture content profiles, it is necessary to take into account the vertical scale of measurements, with an averaging window size of the order of the antenna length and the Fresnel zone width. Flow parameters cannot be uniquely identified, showing that the system is over parameterized with respect to the information content of the (nearly stationary) radar profiles. Estimates of travel time of water across the vadose zone are derived from the simulation results.     

Aquifer Classification for the UIC Program: Prototype Studies in New Mexico

Three case studies from New Mexico illustrate methods by which aquifers can be classified for purposes of the Federal Underground Injection Control program. The principal technique involves preparation of hydrogeologic maps or cross sections which display information on the permeability of rock units and the dissolved solids content of formation fluids. Because deep water wells are lacking in most areas, the analysis normally requires considerable interpretation of geological and geophysical logs collected by energy and mineral companies, plus use of a general model or concept about regional hydrogeology. Injection of waste fluids into aquifers containing water with less than 10,000 mg/1 dissolved solids is not allowed unless an exemption is justified by economic, engineering and other factors. Based on the case studies, regulatory exemptions will be possible for aquifers which are hydrocarbon or mineral-producing, or which are important for brine disposal purposes.     

A study of the geoelectrical properties of peatlands and their influence on ground-penetrating radar surveying1

Geophysical surveys and chemical analyses on cores were carried out in three Ontario peatlands, from which we have gained a better understanding of the peat properties that control the geophysical responses. The electrical conductivity depends linearly on the concentration of total dissolved solids in the peat pore waters and the pore waters in turn bear the ionic signatures of the underlying mineral sediments. The ionic concentration, and thus the electrical conductivity, increase linearly from the surface to basement. The average bulk electrical conductivity of peatlands at Ellice Marsh, near Stratford, and at Wally Creek Area Forest Drainage Project, near Cochrane, are of the order of 25 mS/m. The Mer Bleue peatland, near Ottawa, has extremely high electrical conductivity, reaching levels of up to 380 mS/m near the base of the peat. The Mer Bleue peatland water has correspondingly high values of total dissolved solids, which originate from the underlying Champlain Sea glaciomarine clays. The dielectric permittivity in peats is largely controlled by the bulk water content. Ground penetrating radar can detect changes in water content greater than 3%, occurring within a depth interval less than 15 cm. The principal peatland interfaces detected are the near-surface aerobic to anaerobic transition and the peat to mineral basement contact. The potential for the successful detection of the basement contact using the radar can be predicted using the radar instrument specifications, estimates of the peatland depth, and either the bulk peat or the peat pore water electrical conductivities. Predicted depths of penetration of up to 10 m for Ellice Marsh and Wally Creek exceed the observed depths of 1 to 2 m. At Mer Bleue, on the other hand, we observe that, as predicted, a 100 MHz signal will penetrate to the base of a 2 m thick peat but a 200 MHz signal will not.     

Characterization of an unstable rock mass based on borehole logs and diverse borehole radar data

Unstable rocky slopes are major hazards to the growing number of people that live and travel though mountainous regions. To construct effective barriers to falling rock, it is necessary to know the positions, dimensions and shapes of structures along which failure may occur. To investigate an unstable mountain slope distinguished by numerous open fracture zones, we have taken advantage of three moderately deep (51.0–120.8 m) boreholes to acquire geophysical logs and record single-hole radar, vertical radar profiling (VRP) and crosshole radar data. We observed spallation zones, displacements and borehole radar velocity and amplitude anomalies at 16 of the 46 discontinuities identified in the borehole optical televiewer images. The results of the VRP and crosshole experiments were disappointing at our study site; the source of only one VRP reflection was determined and the crosshole velocity and amplitude tomograms were remarkably featureless. In contrast, much useful structural information was provided by the single-hole radar experiments. Radar reflections were recorded from many surface and borehole fracture zones, demonstrating that the strong electrical property contrasts of these features extended some distance into the adjacent rock mass. The single-hole radar data suggested possible connections between 6 surface and 4 borehole fractures and led to the discovery of 5 additional near-surface fracture zones. Of particular importance, they supplied key details on the subsurface geometries and minimum subsurface lengths of 8 of the 10 previously known surface fracture zones and all of the newly discovered ones. The vast majority of surface fracture zones extended at least 40–60 m into the subsurface, demonstrating that their depth and surface dimensions are comparable.     

Time‐Lapse Electrical Geophysical Monitoring of Amendment‐Based Biostimulation

Biostimulation is increasingly used to accelerate microbial remediation of recalcitrant groundwater contaminants. Effective application of biostimulation requires successful emplacement of amendment in the contaminant target zone. Verification of remediation performance requires postemplacement assessment and contaminant monitoring. Sampling‐based approaches are expensive and provide low‐density spatial and temporal information. Time‐lapse electrical resistivity tomography (ERT) is an effective geophysical method for determining temporal changes in subsurface electrical conductivity. Because remedial amendments and biostimulation‐related biogeochemical processes often change subsurface electrical conductivity, ERT can complement and enhance sampling‐based approaches for assessing emplacement and monitoring biostimulation‐based remediation. Field studies demonstrating the ability of time‐lapse ERT to monitor amendment emplacement and behavior were performed during a biostimulation remediation effort conducted at the Department of Defense Reutilization and Marketing Office (DRMO) Yard, in Brandywine, Maryland, United States. Geochemical fluid sampling was used to calibrate a petrophysical relation in order to predict groundwater indicators of amendment distribution. The petrophysical relations were field validated by comparing predictions to sequestered fluid sample results, thus demonstrating the potential of electrical geophysics for quantitative assessment of amendment‐related geochemical properties. Crosshole radar zero‐offset profile and borehole geophysical logging were also performed to augment the data set and validate interpretation. In addition to delineating amendment transport in the first 10 months after emplacement, the time‐lapse ERT results show later changes in bulk electrical properties interpreted as mineral precipitation. Results support the use of more cost‐effective surface‐based ERT in conjunction with limited field sampling to improve spatial and temporal monitoring of amendment emplacement and remediation performance.     

Application of Electromagnetic Logging to Contamination Investigations in Glacial Sand-and-Gravel Aquifers

Electromagnetic (EM) logging provides an efficient method for high-resolution, vertical delineation of electrically conductive contamination in glacial sand-and-gravel aquifers. LM. gamma, and lithologic logs and specific conductance data from sand-and-gravel aquifers at five sites in the northeastern United States were analyzed to define the relation of KM conductivity to aquifer lithology and water quality. Municipal waste disposal, septic waste discharge, or highway deicing salt application at these sites has caused contaminant plumes in which the dissolved solids concentration and specific conductance of ground water exceed background levels by as much as 10 to 20 limes.
The major hydrogeologic factors that affected KM log response at the five sites were the dissolved solids concentration of the ground water and the silt and clay content in the aquifer. KM conductivity of sand and gravel with uncontaminated water ranged from less than 5 to about 10 millisiemens per meter (mS/m); that of silt and clay zones ranged from about 15 to 45 mS/m: and that of the more highly contaminated zones in sand and gravel ranged from about 10 to more than 80 mS/m. Specific conductance of water samples from screened intervals in sand and gravel at selected monitoring well installations was significantly correlated with KM conductivity.
CM logging can be used in glacial sand-and-gravel aquifer investigations to (1) determine optimum depths for the placement of monitoring well screens: (2) provide a nearly continuous vertical profile of specific conductance to complement depth-specific water quality samples; and (3) identify temporal changes in water quality through sequential logging. Detailed lithologic or gamma logs, preferably both, need to be collected along with the F.M logs to define zones in which elevated EM conductivity is caused by the presence of sill and clay beds rather than contamination.     

On Earth’s Mantle Constitution and Structure from Joint Analysis of Geophysical and Laboratory-Based Data: An Example

Determining Earth’s structure is a fundamental goal of Earth science, and geophysical methods play a prominent role in investigating Earth’s interior. Geochemical, cosmochemical, and petrological analyses of terrestrial samples and meteoritic material provide equally important insights. Complementary information comes from high-pressure mineral physics and chemistry, i.e., use of sophisticated experimental techniques and numerical methods that are capable of attaining or simulating physical properties at very high pressures and temperatures, thereby allowing recovered samples from Earth’s crust and mantle to be analyzed in the laboratory or simulated computationally at the conditions that prevail in Earth’s mantle and core. This is particularly important given that the vast bulk of Earth’s interior is geochemically unsampled. This paper describes a quantitative approach that combines data and results from mineral physics, petrological analyses of mantle minerals, and geophysical inverse calculations, in order to map geophysical data directly for mantle composition (major element chemistry and water content) and thermal state. We illustrate the methodology by inverting a set of long-period electromagnetic response functions beneath six geomagnetic stations that cover a range of geological settings for major element chemistry, water content, and thermal state of the mantle. The results indicate that interior structure and constitution of the mantle can be well-retrieved given a specific set of measurements describing (1) the conductivity of mantle minerals, (2) the partitioning behavior of water between major upper mantle and transition-zone minerals, and (3) the ability of nominally anhydrous minerals to store water in their crystal structures. Specifically, upper mantle water contents determined here bracket the ranges obtained from analyses of natural samples, whereas transition-zone water concentration is an order-of-magnitude greater than that of the upper mantle and appears to vary laterally underneath the investigated locations.     

Mapping a Brine Plume Using Surface Geophysical Methods in Conjunction with Ground Water Quality Data

Oil field brine was applied to a gravel roadbed at an instrumented study site in Newark, Ohio, to simulate the use of brine as a deicer on roads in certain areas of Ohio. Brine was applied on a weekly basis eight times during the winter of 1988 as part of the deicing simulation. Eleven wells were installed at the site prior to brine application to permit collection of ground water quality samples. Surface geophysical methods — electrical resistivity and electromagnetic conductivity — were used to map the resulting brine plume. The accuracy of the methods was evaluated by comparing geophysical and ground water quality data. The presence of brine in ground water resulted in a decrease in resistivity and an increase in conductivity. Specific conductance measured in the field was used as a general indicator of the presence of the brine plume in ground water. Chloride concentration was an indicator of brine in the ground water. Results of the surface geophysical surveys correlated best with chloride and dissolved solids concentrations, and with specific conductance in ground water. The surface geophysical methods were found to be useful for qualitative interpretations of ground water quality changes resulting from the application of brine on roads.     

Geoelectrical and hydrochemical investigations for characterizing the salt water intrusion in the Khanasser valley, northern Syria

An integrated approach of geoelectrical and hydrochemical investigation surveys was proposed for indicating contact regions between saline and fresh groundwater in the Khanasser valley region, northern Syria. The qualitative and quantitative interpretations of 34 vertical electrical soundings (VES) enable to characterize the salt water intrusion laterally and vertically. The established iso-apparent resistivity maps for different AB/2 spacings obviously indicate the presence of a lowresistivity (less than 4 Ohm·m) zone related to the salt water intrusion in the Quaternary and Paleogene deposits. The different hydrochemical and geophysical parameters, such as electrical resistivity, total dissolved solids (TDS) and major ions concentrations used to characterize the salt water intrusion gave almost similar results in locating and mapping the different boundaries of the groundwater salinity. The proposed approach is useful for mapping the interface between different groundwater qualities, and can be therefore used to successfully characterize the salt water intrusion phenomenon in other semi-arid regions. The application of such an approach is a powerful tool and can be used for water resource management in the water scarce areas.     

Direct-Push Electrical Conductivity Logging for High-Resolution Hydrostratigraphic Characterization

Fine-scale hydrostratigraphic features often play a critical role in controlling ground water flow and contaminant transport. Unfortunately, many conventional drilling- and geophysics-based approaches are rarely capable of describing these features at the level of detail needed for contaminant predictions and remediation designs. Previous work has shown that direct-push electrical conductivity (EC) logging can provide information about site hydrostratigraphy at a scale of relevance for contaminant transport investigations in many unconsolidated settings. In this study, we evaluate the resolution and quality of that information at a well-studied research site that is underlain by highly stratified alluvial sediments. Geologic and hydrologic data, conventional geophysical logs, and particle-size analyses are used to demonstrate the capability of direct-push EC logging for the delineation of fine-scale hydrostratigraphic features in saturated unconsolidated formations. When variations in pore-fluid chemistry are small, the electrical conductivity of saturated media is primarily a function of clay content, and hydrostratigraphic features can be described at a level of detail (<2.5 cm in thickness) that has not previously been possible in the absence of continuous cores. Series of direct-push EC logs can be used to map the lateral continuity of layers with non-negligible clay content and to develop important new insights into flow and transport at a site. However, in sand and gravel intervals with negligible clay, EC logging provides little information about hydrostratigraphic features. As with all electrical logging methods, some site-specific information about the relative importance of fluid and sediment contributions to electrical conductivity is needed. Ongoing research is directed at developing direct-push methods that allow EC logging, water sampling, and hydraulic testing to be done concurrently.     

淮北临涣矿采煤沉陷区不同水体水化学特征及其影响因素

为研究淮北临涣矿采煤沉陷区不同水体的补给水源及溶质来源,在现场调查的基础上,系统采集丰水期、平水期、枯水期沉陷区积水、地表河水和浅层地下水样进行测试分析,采用Piper三线图、Gibbs图和因子分析方法,对不同水体水化学特征及其影响因素进行讨论.结果表明:地表水水体总溶解性固体(TDS)质量浓度表现为枯水期丰水期平水期,浅层地下水表现为枯水期平水期丰水期,地表水TDS质量浓度明显高于浅层地下水.地表水中主要阴阳离子为Na~+、Cl~-和SO_4~(2-),水化学类型主要为SO_4~(2-)-Cl~--Na~+型;浅层地下水离子以HCO_3~-、Ca~(2+)和Mg~(2+)为主,表现为HCO_3~--Ca~(2+)-Mg~(2+)型.结合Gibbs图和因子分析可知,地表水受蒸发作用、地表径流以及采煤活动等因素影响,浅层地下水在一定程度上体现出大气降水和地表水补给的特点,受岩石风化作用影响较为明显.     

Deriving TDS values in coarse sediments from long normal and electromagnetic logs

This paper presents an empirical relationship of quantitatively linked electromagnetic (EM) borehole recordings of the total dissolved solids (TDS) in pore water in the Quaternary deposits of the Belgian coastal plain. First, the long normal (LN) logs are linked to EM logs, then the already developed relationships between LN resistivity measurements and the TDS values are rewritten for EM recordings. The main parameter in these equations is the formation factor, which is derived from ground water analyses and LN logs through Archie's law. The EM recording has several advantages compared to the LN logs. The EM analysis allows measuring in PVC-cased wells and is not hindered by the invasion zone around the well. Furthermore, it has a high vertical resolution. LN logs can be measured only once, after drilling a well; EM recordings can be repeated several times in monitoring wells, which allows the gathering of time-dependent data over a complete vertical cross section. Such data could be obtained with LN logs only in wells with screens over the full-depth interval, which causes a hydraulic short circuit. This short circuit can result in a large artificial flow through the well between different levels, resulting in a salinity profile, which is no longer representative for the studied site. Remediation against short circuiting is a reduction of the screened interval, which strongly reduces the gathered information. The application of the derived equations is one of setting up a monitoring network along the Belgian coast to monitor the trend in salinity levels and comparing present salinity levels with older LN recordings to investigate the salinity changes in the last 30 years. Deep wells already present in the Belgiancoastal plain can then be used to monitor both the fresh water head changes and the salt water evolution. The technique has also been used for parameter identification for which real concentration measurements were needed.     

Hydrochemical aspects of the negro river,patagonia, argentina

Water samples collected in the Negro River drainage basin (Patagonia, Argentina) during two different field trips—the first in May 1972 and the second in January 1973—were analysed for the major dissolved inorganic constituents. Q-mode factor analyses demonstrate that two factors account for over 98 per cent of the variance in both sample suites. They represent: a) the direct influence of solution of evaporites and associated carbonate rocks in an arid environment-type; b) the influence of the weathering products of metamorphic, igneous rocks, and varied non-calcareous sediments. R-mode factor analyses of the two sets of data suggest that although the main mechanisms controlling water chemistry were essentially the same, the principal sources of dissolved solids differed significantly in each separate occasion. Further, the major mechanisms controlling its water chemistry are within a wide spectrum of the so-called ‘rock dominance’.     

Push-pull tests for assessing in situ aerobic cometabolism

Three types of single-well push-pull tests were developed for use in assessing the feasibility of in situ aerobic cometabolism of chlorinated aliphatic hydrocarbons (CAHs). These included transport tests, biostimulation tests, and activity tests. Transport tests are conducted to evaluate the mobility of solutes used in subsequent tests. These included bromide or chloride (conservative tracers), propane (growth substrate), ethylene, propylene (CAH surrogates), dissolved oxygen (electron acceptor), and nitrate (a minor nutrient). Tests were conducted at an experimental wellfield of Oregon State University. At this site, extraction phase breakthrough curves for all solutes were similar, indicating apparent conservative transport of the dissolved gases and nitrate prior to biostimulation. Biostimulation tests were conducted to stimulate propane-utilizing activity of indigenous microorganisms and consisted of sequential injections of site ground water containing dissolved propane and oxygen. Biostimulation was detected by the increase in rates of propane and oxygen utilization after each injection. Activity tests were conducted to quantify rates of substrate utilization and to confirm that CAH-transforming activity had likely been stimulated. In particular, the transformation of injected CAH surrogates ethylene and propylene to the cometabolic byproducts ethylene oxide and propylene oxide provided evidence that activity of the monooxygenase enzyme system, responsible for aerobic cometabolic transformations of CAHs, had likely also been stimulated. Estimated zero-order transformation rates decreased in the order propane > ethylene > propylene. The series of push-pull tests developed and field tested in this study should prove useful for conducting rapid, low-cost feasibility assessments for in situ aerobic cometabolism of CAHs.     

Comparison of formation and fluid-column logs in a heterogeneous basalt aquifer

Deep observation boreholes in the vicinity of active production wells in Honolulu, Hawaii, exhibit the anomalous condition that fluid-column electrical conductivity logs and apparent profiles of pore-water electrical conductivity derived from induction conductivity logs are nearly identical if a formation factor of 12.5 is assumed. This condition is documented in three boreholes where fluid-column logs clearly indicate the presence of strong borehole flow induced by withdrawal from partially penetrating water-supply wells. This result appears to contradict the basic principles of conductivity-log interpretation. Flow conditions in one of these boreholes was investigated in detail by obtaining flow profiles under two water production conditions using the electromagnetic flowmeter. The flow-log interpretation demonstrates that the fluid-column log resembles the induction log because the amount of inflow to the borehole increases systematically upward through the transition zone between deeper salt water and shallower fresh water. This condition allows the properties of the fluid column to approximate the properties of water entering the borehole as soon as the upflow stream encounters that producing zone. Because this condition occurs in all three boreholes investigated, the similarity of induction and fluid-column logs is probably not a coincidence, and may relate to aquifer response under the influence of pumping from production wells.     

Delineation of sulphide ore-zones by borehole radar tomography at Hellyer Mine, Australia

Velocity and absorption tomograms are the two most common forms of presentation of radar tomographic data. However, mining personnel, geophysicists included, are often unfamiliar with radar velocity and absorption. In this paper, general formulae are introduced, relating velocity and attenuation coefficient to conductivity and dielectric constant. The formulae are valid for lossy media as well as high-resistivity materials. The transformation of velocity and absorption to conductivity and dielectric constant is illustrated via application of the formulae to radar tomograms from the Hellyer zinc–lead–silver mine, Tasmania, Australia. The resulting conductivity and dielectric constant tomograms constructed at Hellyer demonstrated the potential of radar tomography to delineate sulphide ore zones.     

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.     

Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods

Abstract

The resolution of the freshwater and saline water aquifers in a coastal terrain (Mahanadi Basin, India) is updated. We analysed electrical borehole log data at four sites and compared the water resistivity regime of the freshwater and saline water zones obtained from electrical borehole logging, with the resistivity regime obtained by interpreting vertical electrical sounding (VES) data. The multilayer VES data interpretation is modified to a simple model, containing only the freshwater zone and the saline water zone. The composite geophysical parameters of the freshwater and saline water zones, in particular the resistivity and longitudinal unit conductance regime, are identified. The resolution obtained from the composite geophysical data analyses is very clear and convincing. The composite longitudinal unit conductance regime of the saline water zones is very high compared to that of the freshwater zones. This makes the identification of the two aquifers easy and increases its reliability. A technique which enables analysis of composite geophysical data of freshwater and saline water zones at VES sites in the vicinity of the borehole log sites is proposed. The significance of longitudinal unit conductance in resolving the freshwater and saline water aquifers is illustrated graphically. The proposed technique is validated by correlating the longitudinal unit conductance and resistivity with the total dissolved solids. The efficiency of the technique is validated by carrying out discriminant function analysis.

Citation Hodlur, G. K., Dhakate, R., Sirisha, T. & Panaskar, D. B. (2010) Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods. Hydrol. Sci. J. 55(3), 414–434.     

Detection of Contaminant Plumes by Borehole Geophysical Logging

Two borehole geophysical methods—electromagnetic induction and natural gamma radiation logs—were used to vertically delineate landfill leachate plumes in a glacial aquifer. Geophysical logs of monitoring wells near two land-fills in a glacial aquifer in west-central Vermont show that borehole geophysical methods can aid in interpretation of geologic logs and placement of monitoring well screens to sample landfill leachate plumes.
Zones of high electrical conductance were delineated from the electromagnetic log in wells near two landfills. Some of these zones were found to correlate with silt and clay units on the basis of drilling and gamma logs. Monitoring wells were screened specifically in zones of high electrical conductivity that did not correlate to a silt or clay unit. Zones of high electrical conductivity that did not correlate to a silt or clay unit were caused by the presence of ground water with a high specific conductance, generally from 1000 to 2370 μS/cm (microsiemens per centimeter at 25 degrees Celsius). Ambient ground water in the study area has a specific conductance of approximately 200 to 400 μS/cm. Landfill leachate plumes were found to be approximately 5 to 20 feet thick and to be near the water table surface.