NMR Logging to Estimate Hydraulic Conductivity in Unconsolidated Aquifers

Nuclear magnetic resonance (NMR) logging provides a new means of estimating the hydraulic conductivity (K) of unconsolidated aquifers. The estimation of K from the measured NMR parameters can be performed using the Schlumberger‐Doll Research (SDR) equation, which is based on the Kozeny–Carman equation and initially developed for obtaining permeability from NMR logging in petroleum reservoirs. The SDR equation includes empirically determined constants. Decades of research for petroleum applications have resulted in standard values for these constants that can provide accurate estimates of permeability in consolidated formations. The question we asked: Can standard values for the constants be defined for hydrogeologic applications that would yield accurate estimates of K in unconsolidated aquifers? Working at 10 locations at three field sites in Kansas and Washington, USA, we acquired NMR and K data using direct‐push methods over a 10‐ to 20‐m depth interval in the shallow subsurface. Analysis of pairs of NMR and K data revealed that we could dramatically improve K estimates by replacing the standard petroleum constants with new constants, optimal for estimating K in the unconsolidated materials at the field sites. Most significant was the finding that there was little change in the SDR constants between sites. This suggests that we can define a new set of constants that can be used to obtain high resolution, cost‐effective estimates of K from NMR logging in unconsolidated aquifers. This significant result has the potential to change dramatically the approach to determining K for hydrogeologic applications.     

Assessment of NMR Logging for Estimating Hydraulic Conductivity in Glacial Aquifers

Glacial aquifers are an important source of groundwater in the United States and require accurate characterization to make informed management decisions. One parameter that is crucial for understanding the movement of groundwater is hydraulic conductivity, K. Nuclear magnetic resonance (NMR) logging measures the NMR response associated with the water in geological materials. By utilizing an external magnetic field to manipulate the nuclear spins associated with ¹H, the time-varying decay of the nuclear magnetization is measured. This logging method could provide an effective way to estimate K at submeter vertical resolution, but the models that relate NMR measurements to K require calibration. At two field sites in a glacial aquifer in central Wisconsin, we collected a total of four NMR logs and obtained measurements of K in their immediate vicinity with a direct-push permeameter (DPP). Using a bootstrap algorithm to calibrate the Schlumberger-Doll Research (SDR) NMR-K model, we estimated K to within a factor of 5 of the DPP measurements. The lowest levels of accuracy occurred in the lower-K (K < 10⁻⁴ m/s) intervals. We also evaluated the applicability of prior SDR model calibrations. We found the NMR calibration parameters varied with K, suggesting the SDR model does not incorporate all the properties of the pore space that control K. Thus, the expected range of K in an aquifer may need to be considered during calibration of NMR-K models. This study is the first step toward establishing NMR logging as an effective method for estimating K in glacial aquifers.     

Hydraulic Conductivity Calibration of Logging NMR in a Granite Aquifer,Laramie Range,Wyoming

In granite aquifers, fractures can provide both storage volume and conduits for groundwater. Characterization of fracture hydraulic conductivity (K) in such aquifers is important for predicting flow rate and calibrating models. Nuclear magnetic resonance (NMR) well logging is a method to quickly obtain near-borehole hydraulic conductivity (i.e., K_NMR) at high-vertical resolution. On the other hand, FLUTe flexible liner technology can produce a K profile at comparable resolution but requires a fluid driving force between borehole and formation. For three boreholes completed in a fractured granite, we jointly interpreted logging NMR data and FLUTe K estimates to calibrate an empirical equation for translating borehole NMR data to K estimates. For over 90% of the depth intervals investigated from these boreholes, the estimated K_NMR are within one order of magnitude of K_FLUTe. The empirical parameters obtained from calibrating the NMR data suggest that “intermediate diffusion” and/or “slow diffusion” during the NMR relaxation time may occur in the flowing fractures when hydraulic aperture are sufficiently large. For each borehole, “intermediate diffusion” dominates the relaxation time, therefore assuming “fast diffusion” in the interpretation of NMR data from fractured rock may lead to inaccurate K_NMR estimates. We also compare calibrations using inexpensive slug tests that suggest reliable K_NMR estimates for fractured rock may be achieved using limited calibration against borehole hydraulic measurements.     

Hydraulic Conductivity from Nuclear Magnetic Resonance Logs in Sediments with Elevated Magnetic Susceptibilities

This study examined the application of slim-hole nuclear magnetic resonance (NMR) tools to estimate hydraulic conductivity (K_NMR) in an unconsolidated aquifer that contains a range of grain sizes (silt to gravel) and high and variable magnetic susceptibilities (MS) (10⁻⁴ to 10⁻² SI). A K calibration dataset was acquired at 1-m intervals in three fully screened wells, and compared to K_NMR estimates using the Schlumberger-Doll research (SDR) equation with published empirical constants developed from previous studies in unconsolidated sediments. While K_NMR using published constants was within an order of magnitude of K, the agreement, overprediction, or underprediction of K_NMR varied with the MS distribution in each well. An examination of the effects of MS on NMR data and site-specific empirical constants indicated that the exponent on T_2ML (n-value in the SDR equation, representing the diffusion regime) was found to have the greatest influence on K_NMR estimation accuracy, while NMR porosity did not improve the prediction of K. K_NMR was further improved by integrating an MS log into the NMR analyses. A first approach detrended T_2ML for the effects of MS prior to calculating K_NMR, and a second approach introduced an MS term into the SDR equation. Both were found to produce similar refinements of K_NMR in intervals of elevated MS. This study found that low frequency NMR logging with short echo times shows promise for sites with moderate to elevated MS levels, and recommends a workflow that examines parameter relationships and integrates MS logs into the estimation of K_NMR.     

The usefulness of outcrop analogue air permeameter measurements for analyzing aquifer heterogeneity: quantifying outcrop hydraulic conductivity and its spatial variability

Saturated hydraulic conductivity (K) is one of the most important parameters determining groundwater flow and contaminant transport in both unsaturated and saturated porous media. Although several well‐established laboratory methods exist for determining K, in situ measurements of this parameter remain very complex and scale dependent. Often, the limited accessibility of subsurface sediments for sampling means an additional impediment to our ability to quantify subsurface K heterogeneity. One potential solution is the use of outcrops as analogues for subsurface sediments. This paper investigates the use of air permeameter measurements on outcrops of unconsolidated sediments to quantify K and its spatial heterogeneity on a broad range of sediment types. The Neogene aquifer in northern Belgium is used as a case study for this purpose. To characterize the variability in K, 511 small‐scale air permeability measurements were performed on outcrop sediments representative over five of the aquifer's lithostratigraphic units. From these measurements, outcrop‐scale equivalent K tensors were calculated using numerical upscaling techniques. Validation of the air permeameter‐based K values by comparison with laboratory constant head K measurements reveals a correlation of 0.93. Overall, the results indicate that hand‐held air permeameters are very efficient and accurate tools to characterize saturated K, as well as its small‐scale variability and anisotropy on a broad range of unconsolidated sediments. The studied outcrops further provided a qualitative understanding of aquifer hydrostratigraphy and quantitative estimates about K variability at the centimetre‐scale to metre‐scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Geostatistical analysis and conditional simulation for estimating the spatial variability of hydraulic conductivity in the Choushui River alluvial fan,Taiwan

This work evaluated the spatial variability and distribution of heterogeneous hydraulic conductivity (K) in the Choushui River alluvial fan in Taiwan, using ordinary kriging (OK) and mean and individual sequential Gaussian simulations (SGS). A baseline flow model constructed by upscaling parameters was inversely calibrated to determine the pumping and recharge rates. Simulated heads using different K realizations were then compared with historically measured heads. A global/local simulated error between simulated and measured heads was analysed to assess the different spatial variabilities of various estimated K distributions. The results of a MODFLOW simulation indicate that the OK realization had the smallest sum of absolute mean simulation errors (SAMSE) and the SGS realizations preserved the spatial variability of the measured K fields. Moreover, the SAMSE increases as the spatial variability of the K field increases. The OK realization yields small local simulation errors in the measured K field of moderate magnitude, whereas the SGS realizations have small local simulation errors in the measured K fields, with high and low values. The OK realization of K can be applied to perform a deterministic inverse calibration. The mean SGS method is suggested for constructing a K field when the application focuses on extreme values of estimated parameters and small calibration errors, such as in a simulation of contaminant transport in heterogeneous aquifers. The individual SGS realization is useful in stochastically assessing the spatial uncertainty of highly heterogeneous aquifers. Copyright © 2004 John Wiley & Sons, Ltd.     

The estimation of shear-wave statics using in situ measurements in marine near-surface sediments

Shear‐wave statics in marine seismic exploration data are routinely too large to be estimated using conventional techniques. Near‐surface unconsolidated sediments are often characterized by low values of V_s and steep velocity gradients. Minor variations in sediment properties at these depths correspond to variations in the shear‐wave velocity and will produce significant static shifts. It is suggested that a significant proportion of the shear‐wave statics solution can be estimated by performing a separate high‐resolution survey to target near‐surface unconsolidated sediments. Love‐wave, shear‐wave refraction and geotechnical measurements were individually used to form high‐resolution near‐surface shear‐wave velocity models to estimate the shear‐wave statics for a designated survey line. Comparisons with predicted statics revealed that shear‐wave statics could not be estimated using a velocity model predicted by substituting geotechnical measurements into empirical relationships. Empirical relationships represent a vast simplification of the factors that control V_s and are therefore not sufficiently sensitive to estimate shear‐wave statics. Refraction measurements are potentially sensitive to short‐wavelength variations in sediment properties when combined with accurate navigational data. Statics estimated from Love‐wave data are less sensitive, and sometimes smoothed in appearance, since interpreted velocity values represent an average both laterally and vertically over the receiver array and the frequency–depth sensitivity range, respectively. For the survey site, statics estimated from near‐surface irregularities using shear‐wave refraction measurements represent almost half the total statics solution. More often, this proportion will be greater when bedrock relief is less.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Evaluation of Hydraulic Properties of Aquitards Using Earthquake‐Triggered Groundwater Variation

The hydraulic properties of aquitards are not easily obtained because monitoring wells are usually installed in aquifers for groundwater resources management. Earthquake‐induced crust stress (strain) triggers groundwater level variations over a short period of time in a large area. These groundwater anomalies can be used to investigate aquifer systems. This study uses a poroelastic model to fit the postseismic variations of groundwater level triggered by the Chi‐Chi earthquake to evaluate the hydraulic properties of aquitards in the Jhoushuei River alluvial fan (JRAF), Taiwan. Six of the adopted eight wells with depths of 70 to 130 m showed good agreement with the recovery theory. The mean hydraulic conductivities (K) of the aquifers for the eight wells are 1.62 × 10^?4 to 9.06 × 10^?4 m/s, and the thicknesses are 18.8 to 46.1 m. The thicknesses of the aquitards are 11.3 to 42.0 m. Under the isotropic assumption for K, the estimated values of K for the aquitards are 3.0 × 10^?8 to 2.1 × 10^?6 m/s, corresponding to a silty medium. The results match the values obtained for the geological material of the drilling core and those reported in previous studies. The estimated values were combined with those given in previous studies to determine the distribution of K in the first two aquitards in the JRAF. The distribution patterns of the aquitards reflect the sedimentary environments and fit the geological material. The proposed technique can be used to evaluate the K value of aquitards using inverse methods. The inversion results can be used in hydrogeological analyses, contaminant modeling, and subsidence evaluation.     

Depth‐dependent hydraulic conductivity distribution patterns of a streambed

Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (K_v) of the channel sediments and the aquifer materials. This method was applied to produce K_v profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the K_v profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of K_v in the channel sediments. This depth‐dependent pattern of K_v distribution for the channel sediments is considered to be produced by hyporheic processes. This K_v‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 John Wiley & Sons, Ltd.     

A Tube Seepage Meter for In Situ Measurement of Seepage Rate and Groundwater Sampling

We designed and evaluated a “tube seepage meter” for point measurements of vertical seepage rates (q), collecting groundwater samples, and estimating vertical hydraulic conductivity (K) in streambeds. Laboratory testing in artificial streambeds show that seepage rates from the tube seepage meter agreed well with expected values. Results of field testing of the tube seepage meter in a sandy‐bottom stream with a mean seepage rate of about 0.5 m/day agreed well with Darcian estimates (vertical hydraulic conductivity times head gradient) when averaged over multiple measurements. The uncertainties in q and K were evaluated with a Monte Carlo method and are typically 20% and 60%, respectively, for field data, and depend on the magnitude of the hydraulic gradient and the uncertainty in head measurements. The primary advantages of the tube seepage meter are its small footprint, concurrent and colocated assessments of q and K, and that it can also be configured as a self‐purging groundwater‐sampling device.     

Prediction of Remediation of a Heterogeneous Aquifer: A Case Study

Contaminant plumes whose characteristic length is smaller than the horizontal integral scale of the hydraulic conductivity, K, are abundant in shallow, phreatic aquifers. In such cases, the aquifer can be regarded as layered, with K being only a function of the vertical coordinate. The heterogeneity of K has a critical role upon the efficiency of remediation of such sites, for example, by Pump and Treat schemes. The expected efficiency is a random variable, with uncertainty. Quantifying this uncertainty can be of great importance to decision making. In this study, we focus on a case study in the coastal aquifer of Israel and compare two different approaches for constructing realizations of K: continuous and indicator. We observe a significant difference between the constructed realizations, which results in a considerable difference in the predicted remediation efficiency and its uncertainty. Furthermore, we study the effect of conditioning the realizations by a rather limited number of K data points. We find that the conditioning results in a major reduction of the uncertainty. In addition, we compare the results of the transport model to a simplified semi‐analytical solution that is based on assuming radial flow. We find a good agreement with the three‐dimensional numerical model. This result illustrates that the simplified solution can be used for prediction of the remediation efficiency when the flow at the plume vicinity can be regarded as radial.     

Negative correlation between porosity and hydraulic conductivity in sand-and-gravel aquifers at Cape Cod,Massachusetts, USA

Although it may be intuitive to think of the hydraulic conductivity K of unconsolidated, coarse-grained sediments as increasing monotonically with increasing porosity Φ, studies have documented a negative correlation between these two parameters under certain grain-size distributions and packing arrangements. This is confirmed at two sites on Cape Cod, Massachusetts, USA, where groundwater investigations were conducted in sand-and-gravel aquifers specifically to examine the interdependency of several aquifer properties using measurements from four geophysical well logs. Along with K and Φ, the electrical resistivity R₀ and the natural gamma activity γ of saturated deposits were determined as functions of depth. Qualitative examination of results from the first site implies a negative correlation between K and Φ that is substantiated by a rigorous multivariate analysis of log data collected from the second site. A principal components analysis describes an over-determined system of inversion equations, with approximately 92% of the cumulative proportion of the total variance being accounted for by only three of the four eigenvectors. A subsequent R-mode factor analysis projects directional trends among the four variables (K, Φ, R₀ and γ), and a negative correlation between K and Φ emerges as the primary result.     

Combining 3D Hydraulic Tomography with Tracer Tests for Improved Transport Characterization

Hydraulic tomography (HT) is a method for resolving the spatial distribution of hydraulic parameters to some extent, but many details important for solute transport usually remain unresolved. We present a methodology to improve solute transport predictions by combining data from HT with the breakthrough curve (BTC) of a single forced‐gradient tracer test. We estimated the three dimensional (3D) hydraulic‐conductivity field in an alluvial aquifer by inverting tomographic pumping tests performed at the Hydrogeological Research Site Lauswiesen close to Tübingen, Germany, using a regularized pilot‐point method. We compared the estimated parameter field to available profiles of hydraulic‐conductivity variations from direct‐push injection logging (DPIL), and validated the hydraulic‐conductivity field with hydraulic‐head measurements of tests not used in the inversion. After validation, spatially uniform parameters for dual‐domain transport were estimated by fitting tracer data collected during a forced‐gradient tracer test. The dual‐domain assumption was used to parameterize effects of the unresolved heterogeneity of the aquifer and deemed necessary to fit the shape of the BTC using reasonable parameter values. The estimated hydraulic‐conductivity field and transport parameters were subsequently used to successfully predict a second independent tracer test. Our work provides an efficient and practical approach to predict solute transport in heterogeneous aquifers without performing elaborate field tracer tests with a tomographic layout.     

Sensitivity analysis of hydraulic conductivity and Manning's n parameters lead to new method to scale effective hydraulic conductivity across model resolutions

Hydrological modelling is an important tool for research, policy, and management, but uncertainty remains about parameters transferability from field observations made at small scale to models at the catchment scale and larger. This uncertainty compels the need to develop parameter relationships that are translatable across scale. In this study, we compare the changes to modelled processes as resolution is coarsened from 100‐m to 1‐km in a topographically complex, 255‐km² Colorado River headwater catchment. We conducted a sensitivity analysis for hydraulic conductivity (K) and Manning's n parameters across four orders of magnitude. Results showed that K acts as a moderator between surface and subsurface contributions to streamflow, whereas n moderates the duration of high intensity, infiltration‐excess flow. The parametric sensitivity analysis informed development of a new method to scale effective hydraulic conductivity across modelling resolutions in order to compensate for the loss of topographic gradients as resolution is coarsened. A similar mathematical relationship between n and lateral resolution changes was not found, possibly because n is also sensitive to time discretization. This research provides an approach to translate hydraulic conductivity parameters from a calibrated coarse model to higher resolutions where the number of simulations are limited by computational demand.     

Cost Comparisons of Aquifer Heterogeneity Characterization Methods

The characterization of heterogeneity in hydraulic conductivity (K) is a major challenge for subsurface remediation projects. There are a number of field studies that compare the K estimates obtained using various techniques, but to our knowledge, no field‐based studies exists that compare the performance of estimated K heterogeneity fields or the associated characterization costs. In this paper, we compare the costs of characterizing the three‐dimensional K heterogeneity and its uncertainty estimates of a glaciofluvial aquifer‐aquitard sequence at a 15 m × 15 m × 18 m field site situated on the University of Waterloo campus. We compare geostatistical analysis of high resolution permeameter K data obtained from repacked core samples in five boreholes and hydraulic tomography analysis of four pumping tests consisting of up to 41 monitoring points per test. Aside from the comparison of costs, we also assess the performance of each method by predicting several pumping tests. Our analysis reveals that hydraulic tomography is somewhat more costly than the geostatistical analysis of high resolution permeameter K data due to the higher capital costs associated with the method. However, the equipment may be reused at other sites; hence these costs may be recovered over the life of the equipment. More significantly, hydraulic tomography is able to capture the most important features of the aquifer‐aquitard sequence leading to more accurate predictions of independent pumping tests. This suggests that more robust remediation systems may be designed if site characterization is performed with hydraulic tomography.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

A Joint Analytic Method for Estimating Aquitard Hydraulic Parameters

The vertical hydraulic conductivity (K_v), elastic (S_ske), and inelastic (S_skv) skeletal specific storage of aquitards are three of the most critical parameters in land subsidence investigations. Two new analytic methods are proposed to estimate the three parameters. The first analytic method is based on a new concept of delay time ratio for estimating K_v and S_ske of an aquitard subject to long‐term stable, cyclic hydraulic head changes at boundaries. The second analytic method estimates the S_skv of the aquitard subject to linearly declining hydraulic heads at boundaries. Both methods are based on analytical solutions for flow within the aquitard, and they are jointly employed to obtain the three parameter estimates. This joint analytic method is applied to estimate the K_v, S_ske, and S_skv of a 34.54‐m thick aquitard for which the deformation progress has been recorded by an extensometer located in Shanghai, China. The estimated results are then calibrated by PEST (Doherty 2005), a parameter estimation code coupled with a one‐dimensional aquitard‐drainage model. The K_v and S_ske estimated by the joint analytic method are quite close to those estimated via inverse modeling and performed much better in simulating elastic deformation than the estimates obtained from the stress‐strain diagram method of Ye and Xue (2005). The newly proposed joint analytic method is an effective tool that provides reasonable initial values for calibrating land subsidence models.     

Ring permeametry: design,operation and error analysis

Assessment of slope stability, soil management or contaminant transport problems usually requires numerous, yet accurate point measurements of permeability. This technical note describes a new method for the rapid field assessment of permeability in near-surface soils or unconsolidated sediments. The procedure is known as ‘ring permeametry’ and is an ex situ core-based method giving measurements which can be guaranteed to be stratum-specific, unlike measurements from some in situ techniques. The potential sources of precision and bias error within the method are quantified and their effect on the uncertainty of permeability estimates is illustrated. © 1997 John Wiley & Sons, Ltd.     

An Integrated Laboratory Method to Measure and Verify Directional Hydraulic Conductivity in Fine‐to‐Medium Sandy Sediments

The constant‐head permeameter test (CHPT) is widely used in sandy samples as a standard method in the laboratory to investigate hydraulic conductivity (K). However, it neither can be used to consistently determine directional hydraulic conductivity (DHC) nor guarantee the comparability of measured K values of samples with different sizes. Therefore, this paper proposes an integrated laboratory method, called modified CHPT (MCHPT), for the efficient determination and verification of consistent DHC values in fine‐to‐medium sandy sediments, based on a new methodological framework. A precise and standardized procedure for preparing the experimental setup of MCHPT was conducted, based on the integrated experimental setup of CHPT and tracer tests. Moreover, a formula was yielded for the time‐optimized sample saturation control. In comparison with grain size‐based methods, the validity of consistent K_h and K_v values determined by MCHPT was convincing.