Vertical Hydraulic Gradient Estimation in Clay Till,Using MiHPT Advanced Direct-Push Technology

The vertical transport of contaminants from source areas is employed in many risk assessment models and screening tools in order to estimate the contaminant mass discharge (CMD) into underlying aquifers. The key parameters for estimating CMD are the contaminant source area and concentration, and the vertical water flux, the latter of which depends on the vertical hydraulic conductivity and the vertical hydraulic gradient in the subsurface. This study focuses on advancing the use of the combined membrane interface probe hydraulic profiling tool (MiHPT) to investigate the vertical hydraulic gradient across a clay till overlying a sandy aquifer at a contaminated site in Denmark. Only the HPT is necessary for the estimate of vertical hydraulic gradient. The hydraulic head, clay till thickness, and resulting vertical hydraulic gradients found using the MiHPT compared well with observations from nearby nested wells. The parameter with the largest discrepancy was the thickness of the clay till. The advantage of the MiHPT is its relatively quick depth discrete access to information regarding subsurface permeability, vertical hydraulic gradients and contaminant distribution across a site. In this case study, performance of additional dissipationtests during the HPT log to acquire determination of the vertical hydraulic gradient increased the cost by 3% compared to standard HPT logs.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

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.     

Vertically Integrated Hydraulic Conductivity: A New Parameter for Groundwater-Surface Water Analysis

Numerical modeling of groundwater-surface water interactions provides vital information necessary for determining the extent of nutrient transport, quantifying water budgets, and delineating zones of ecological support. The hydrologic data that drive these models are often collected at disparate scales and subsequently incorporated into numerical models through upscaling techniques such as piecewise constancy or geostatistical methods. However, these techniques either use basic interpolation methods, which often simplifies the system of interest, or utilize complex statistical methods that are computationally expensive, time consuming, and generate complex subsurface configurations. We propose a bulk parameter termed “vertically integrated hydraulic conductivity” (K_V), and defined as the depth-integrated resistance to fluid flow sensed at the groundwater-surface water interface, as an alternative to hydraulic conductivity when investigating vertical fluxes across the groundwater-surface water interface. This bulk parameter replaces complex subsurface configurations in situations dominated by vertical fluxes and where heterogeneity is not of primary importance. To demonstrate the utility of K_V, we extracted synthetic temperature time series data from a forward numerical model under a variety of scenarios and used those data to quantify vertical fluxes using the amplitude ratio method. These quantified vertical fluxes and the applied hydraulic head gradient were subsequently input into Darcy's Law and used to quantify K_V. This K_V was then directly compared to the equivalent hydraulic conductivity (K_T) assuming an infinitely extending layer. Vertically integrated hydraulic conductivity allows for more accurate and robust flow modeling across the groundwater-surface water interface in instances where complex heterogeneities are not of primary concern.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Hydraulic conductivity in upland blanket peat: measurement and variability

A key parameter used in wetland hydrological and landform development models is hydraulic conductivity. Head recovery tests are often used to measure hydraulic conductivity, but the calculation techniques are usually confined to rigid soil theory. This is despite reports demonstrating the misapplication of rigid soil theory to non‐rigid soils such as peats. Although values of hydraulic conductivity calculated using compressible techniques have been presented for fenland peats, these data have never, to the authors' knowledge, been compared with such calculations in other peat types. Head recovery tests (slug withdrawal) were performed on piezometers at depths ranging from 10 to 80 cm from the surface on north Pennines blanket peats. Results were obtained using both rigid and compressible soil theories, thus allowing comparison of the two techniques. Compressible soil theory gives values for hydraulic conductivity that are typically a factor of five times less than rigid soil calculations. Hydraulic conductivity is often assumed to decrease with depth in upland peats, but at the study site in the northern Pennines it was not found to vary significantly with depth within the range of peat depths sampled. The variance within depth categories was not significantly different to the variance between depth categories showing that individual peat layers did not have characteristic hydraulic conductivity values. Thus, large lateral and vertical differences in hydraulic conductivity over short distances create problems for modelling but may help account for the high frequency of preferential flow pathways within what is otherwise a low matrix hydraulic conductivity peat. Hydraulic conductivity was found to vary significantly between sampling sites, demonstrating that hillslope‐ or catchment‐scale variability may be more important than plot‐scale variability. Values for compressibility of the peats are also reported. These generally decline with depth, and they also vary significantly between sampling sites. There are implications for the way in which measurements of hydraulic conductivity and other properties of blanket peat are interpreted, as the effects of environmental change in one part of a peat catchment may be very different to those in another. Copyright © 2002 John Wiley & Sons, Ltd.     

The Influence of Topology on Hydraulic Conductivity in a Sand-and-Gravel Aquifer

A field experiment consisting of geophysical logging and tracer testing was conducted in a single well that penetrated a sand-and-gravel aquifer at the U.S. Geological Survey Toxic Substances Hydrology research site on Cape Cod, Massachusetts. Geophysical logs and flowmeter/pumping measurements were obtained to estimate vertical profiles of porosity ϕ, hydraulic conductivity K, temperature, and bulk electrical conductivity under background, freshwater conditions. Saline-tracer fluid was then injected into the well for 2 h and its radial migration into the surrounding deposits was monitored by recording an electromagnetic-induction log every 10 min. The field data are analyzed and interpreted primarily through the use of Archie's (1942) law to investigate the role of topological factors such as pore geometry and connectivity, and grain size and packing configuration in regulating fluid flow through these coarse-grained materials. The logs reveal no significant correlation between K and ϕ, and imply that groundwater models that link these two properties may not be useful at this site. Rather, it is the distribution and connectivity of the fluid phase as defined by formation factor F, cementation index m, and tortuosity α that primarily control the hydraulic conductivity. Results show that F correlates well with K, thereby indicating that induction logs provide qualitative information on the distribution of hydraulic conductivity. A comparison of α, which incorporates porosity data, with K produces only a slightly better correlation and further emphasizes the weak influence of the bulk value of ϕ on K.     

DOWN HOLE FLOWMETER ANALYSIS USING AN ASSOCIATED CALIPER LOG

Abstract From flowmeter interpretation studies using wells with no available caliper data has emerged information which is also pertinent to flowmeter analysis in conjunction with a caliper log. A simple method for interpreting flowmeter and caliper log pairs incorporating this insight has been developed. Location of constant-flow regions, inflow and outflow zones is done by manually comparing the shapes of the caliper and flowmeter logs. A flowrate log is then produced using quantitative volumetric flowrates calculated for several depths in the well. The log contains all the significant flow information without spurious noise and would be suitable for converting to an apparent hydraulic conductivity log.     

Stochastic ground water flow simulation with a fracture zone continuum model

A method is presented for incorporating the hydraulic effects of vertical fracture zones into two-dimensional cell-based continuum models of ground water flow and particle tracking. High hydraulic conductivity features are used in the model to represent fracture zones. For fracture zones that are not coincident with model rows or columns, an adjustment is required for the hydraulic conductivity value entered into the model cells to compensate for the longer flowpath through the model grid. A similar adjustment is also required for simulated travel times through model cells. A travel time error of less than 8% can occur for particles moving through fractures with certain orientations. The fracture zone continuum model uses stochastically generated fracture zone networks and Monte Carlo analysis to quantify uncertainties with simulated advective travel times. An approach is also presented for converting an equivalent continuum model into a fracture zone continuum model by establishing the contribution of matrix block transmissivity to the bulk transmissivity of the aquifer. The methods are used for a case study in west-central Florida to quantify advective travel times from a potential wetland rehydration site to a municipal supply wellfield. Uncertainties in advective travel times are assumed to result from the presence of vertical fracture zones, commonly observed on aerial photographs as photolineaments.     

Signatures in flowing fluid electric conductivity logs

Flowing fluid electric conductivity logging provides a means to determine hydrologic properties of fractures, fracture zones, or other permeable layers intersecting a borehole in saturated rock. The method involves analyzing the time-evolution of fluid electric conductivity (FEC) logs obtained while the well is being pumped and yields information on the location, hydraulic transmissivity, and salinity of permeable layers. The original analysis method was restricted to the case in which flows from the permeable layers or fractures were directed into the borehole (inflow). Recently, the method was adapted to permit treatment of both inflow and outflow, including analysis of natural regional flow in the permeable layer. A numerical model simulates flow and transport in the wellbore during flowing FEC logging, and fracture properties are determined by optimizing the match between simulation results and observed FEC logs. This can be a laborious trial-and-error procedure, especially when both inflow and outflow points are present. Improved analyses methods are needed. One possible tactic would be to develop an automated inverse method, but this paper takes a more elementary approach and focuses on identifying the signatures that various inflow and outflow features create in flowing FEC logs. The physical insight obtained provides a basis for more efficient analysis of these logs, both for the present trial and error approach and for a potential future automated inverse approach. Inflow points produce distinctive signatures in the FEC logs themselves, enabling the determination of location, inflow rate, and ion concentration. Identifying outflow locations and flow rates typically requires a more complicated integral method, which is also presented in this paper.     

Multilayer analytic element modeling of radial collector wells

A new multilayer approach is presented for the modeling of ground water flow to radial collector wells. The approach allows for the inclusion of all aspects of the unique boundary condition along the lateral arms of a collector well, including skin effect and internal friction losses due to flow in the arms. The hydraulic conductivity may differ between horizontal layers within the aquifer, and vertical anisotropy can be taken into account. The approach is based on the multilayer analytic element method, such that regional flow and local three-dimensional detail may be simulated simultaneously and accurately within one regional model. Horizontal flow inside a layer is computed analytically, while vertical flow is approximated with a standard finite-difference scheme. Results obtained with the proposed approach compare well to results obtained with three-dimensional analytic element solutions for flow in unconfined aquifers. The presented approach may be applied to predict the yield of a collector well in a regional setting and to compute the origin and residence time, and thus the quality, of water pumped by the collector well. As an example, the addition of three lateral arms to a collector well that already has three laterals is investigated. The new arms are added at an elevation of 2 m above the existing laterals. The yield increase of the collector well is computed as a function of the lengths of the three new arms.     

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.     

Parameter equivalence for the Gardner and van Genuchten soil hydraulic conductivity functions for steady vertical flow with inclusions

The analytic element method is well suited for the Gardner hydraulic conductivity function, but is limited in describing real soils. Therefore, parameter equivalence between the van Genuchten and Gardner hydraulic conductivity functions is explored for the case of steady vertical flow through a homogeneous medium with a single inclusion, i.e., a binary soil. The inclusion has different hydraulic parameters than the background medium. Equivalence is established using three methods: (1) effective capillary drive; (2) capillary length; (3) and a least-squares optimization method that aims to fit a Gardner function to a corresponding van Genuchten function by minimizing the difference in log conductivity over a specified pressure range. Comparisons between hydraulic models are made based on scatterplots of pressure head and the vertical Darcian flux obtained using a finite-element numerical solution with both constitutive relations. For applicability of an equivalent Gardner function over a broad range of pressure heads, the crossover pressure must be maintained between the two parametric functions. The crossover pressure is defined as the pressure in which the hydraulic conductivity of the inclusion is equal to the background. It can be shown that a hybrid methodology of preserving the crossover pressure exactly and using the effective capillary drive will result in hydraulic parameters that are easily obtained and provide good agreement between the conductivity functions of the GR model to the VG model.     

Steady‐State Homogeneous Approximations of Vertical Velocity from EC Profiles

Electrical conductivity (EC) logs were obtained by both open‐borehole logging and passive multilevel sampling (MLS) in an observation borehole penetrating the Coastal Aquifer in Tel Aviv, Israel. Homogeneous vertical velocities for a 70‐m thick subaquifer were approximated from each profile using a steady‐state advection‐diffusion model. The open‐borehole log led to an overestimation of the steady‐state upward advective flux of deep brines (vertical velocity of 0.95 cm/yr as compared to 0.07 cm/yr for the MLS profile). The combination of depth‐dependent data and the suggested simple modeling approach comprises a method for assessing the vertical location of salinity sources and the nature of salt transport from them (i.e., advective vs. diffusive). However, in this case, the easily obtained open‐borehole logs should not be used for collecting depth‐dependent data.     

Application of harmonic analysis of water levels to determine vertical hydraulic conductivities in clay-rich aquitards

A harmonic analysis method was used to determine vertical hydraulic conductivities (Kv) in geologic media between vertically separated piezometers using water level measurements. In this method, each water level time series was filtered and then decomposed using harmonic analysis into a sum of trigonometric components. The phase and amplitude of each harmonic function were calculated. These data were used to estimate Kv values between vertically separated data sets assuming one-dimensional transient flow. The method was applied to water level data collected from nested piezometers at two thick clay-rich till aquitards in Saskatchewan, Canada. At one site, routine water levels were measured in 12 piezometers (installed between 1 and 29 m below ground surface) since installation (1995). At the other site, water levels were measured in seven piezometers (installed between 4 and 53 m below ground surface) since installation (1998-1999). The Kv calculated using harmonic analysis decreased with depth below the water table at both sites, approaching matrix estimates of hydraulic conductivity between 10 and 11 m and between 21 and 43 m below ground surface. These depths reflected the depth of extensive vertical fracturing at the sites and showed that the depth of fracturing may be site specific.     

Modeling Steady Sea Water Intrusion with Single‐Density Groundwater Codes

Steady interface flow in heterogeneous aquifer systems is simulated with single‐density groundwater codes by using transformed values for the hydraulic conductivity and thickness of the aquifers and aquitards. For example, unconfined interface flow may be simulated with a transformed model by setting the base of the aquifer to sea level and by multiplying the hydraulic conductivity with 41 (for sea water density of 1025 kg/m³). Similar transformations are derived for unconfined interface flow with a finite aquifer base and for confined multi‐aquifer interface flow. The head and flow distribution are identical in the transformed and original model domains. The location of the interface is obtained through application of the Ghyben‐Herzberg formula. The transformed problem may be solved with a single‐density code that is able to simulate unconfined flow where the saturated thickness is a linear function of the head and, depending on the boundary conditions, the code needs to be able to simulate dry cells where the saturated thickness is zero. For multi‐aquifer interface flow, an additional requirement is that the code must be able to handle vertical leakage in situations where flow in an aquifer is unconfined while there is also flow in the aquifer directly above it. Specific examples and limitations are discussed for the application of the approach with MODFLOW. Comparisons between exact interface flow solutions and MODFLOW solutions of the transformed model domain show good agreement. The presented approach is an efficient alternative to running transient sea water intrusion models until steady state is reached.     

Using data assimilation method to calibrate a heterogeneous conductivity field conditioning on transient flow test data

A data assimilation method is developed to calibrate a heterogeneous hydraulic conductivity field conditioning on transient pumping test data. The ensemble Kalman filter (EnKF) approach is used to update model parameters such as hydraulic conductivity and model variables such as hydraulic head using available data. A synthetical two-dimensional flow case is used to assess the capability of the EnKF method to calibrate a heterogeneous conductivity field by assimilating transient flow data from observation wells under different hydraulic boundary conditions. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating continuous hydraulic head measurements and the hydraulic boundary condition will significantly affect the simulation results. For our cases, after a few data assimilation steps, the assimilated conductivity field with four Neumann boundaries matches the real field well while the assimilated conductivity field with mixed Dirichlet and Neumann boundaries does not. We found in our cases that the ensemble size should be 300 or larger for the numerical simulation. The number and the locations of the observation wells will significantly affect the hydraulic conductivity field calibration.     

An Evaluation of Interpolation Methodologies for Generating Three-Dimensional Hydraulic Property Distributions from Measured Data

The process of attempting to model ground-water systems requires a good understanding of the spatial variation of aquifer hydraulic properties. The capabilities of the more recent innovative flowmeters such as the electromagnetic and heat pulse flowmeters provide the sensitivity to measure ambient flows and pump-induced flows. These flowmeters provide the measurements of pump-induced vertical flows which are analyzed to obtain vertical variations in horizontal hydraulic conductivity, K(z). With discrete areal K-values, K(x, y), and vertical profiles of K, provided by multiwell testing, the essential elements are present to produce a three-dimensional hydraulic conductivity field. The advent of these new flow measuring devices has contributed much to the motivation behind this paper. This paper presents the results of applying deterministic and stochastic methodology to the three-dimensional interpolation of hydraulic properties, specifically, hydraulic conductivity, K. Three of the approaches applied in this paper are deterministic in nature, inverse-distance weighting, inverse-distance-squared weighting, and ordinary kriging, while the fourth is a stochastic approach based on self-affine fractals. All of the methods are applied to measured data collected from 14 wells at a site in the United States near Mobile, Alabama. The three-dimensional K-distributions generated by each of the methods are used as inputs to an advective based transport model with the resulting model output compared to a two-well tracer study run previously at the same site.     

Geostatistics applied to cross-well reflection seismic for imaging carbonate aquifers

Cross-well seismic reflection data, acquired from a carbonate aquifer at Port Mayaca test site near the eastern boundary of Lake Okeechobee in Martin County, Florida, are used to delineate flow units in the region intercepted by two wells. The interwell impedance determined by inversion from the seismic reflection data allows us to visualize the major boundaries between the hydraulic units. The hydraulic (flow) unit properties are based on the integration of well logs and the carbonate structure, which consists of isolated vuggy carbonate units and interconnected vug systems within the carbonate matrix. The vuggy and matrix porosity logs based on Formation Micro-Imager (FMI) data provide information about highly permeable conduits at well locations. The integration of the inverted impedance and well logs using geostatistics helps us to assess the resolution of the cross-well seismic method for detecting conduits and to determine whether these conduits are continuous or discontinuous between wells. A productive water zone of the aquifer outlined by the well logs was selected for analysis and interpretation. The ELAN (Elemental Log Analysis) porosity from two wells was selected as primary data and the reflection seismic-based impedance as secondary data. The direct and cross variograms along the vertical wells capture nested structures associated with periodic carbonate units, which correspond to connected flow units between the wells. Alternatively, the horizontal variogram of impedance (secondary data) provides scale lengths that correspond to irregular boundary shapes of flow units. The ELAN porosity image obtained by cokriging exhibits three similar flow units at different depths. These units are thin conduits developed in the first well and, at about the middle of the interwell separation region, these conduits connect to thicker flow units that are intercepted by the second well. In addition, a high impedance zone (low porosity) at a depth of about 275 m, after being converted to ELAN porosity, is characterized as a more confined low porosity structure. This continuous zone corresponds to a permeability barrier in the carbonate aquifer that separates the three connected conduits observed in the cokriging image. In the zones above and below this permeability barrier, the water production is very high, which agrees with water well observations at the Port Mayaca aquifer.