Evaluation of methods for determination of hydraulic properties in an aquifer–aquitard system hydrologically connected to a river

Field and laboratory methods have been used to determine the hydraulic properties in a multiple-layer aquifer–aquitard system that is hydrologically connected to a river. First, hypothetical pumping tests in aquifer–aquitard systems were performed to evaluate the feasibility of MODFLOW-PEST in determining these parameters. Sensitivity analyses showed that: the horizontal hydraulic conductivity in the aquifer has the highest composite sensitivity; the vertical hydraulic conductivity has higher composite sensitivity than the horizontal hydraulic conductivity in the aquitard; and a partial penetration pumping well in an aquifer layer can improve the quality of the estimated parameters. This inverse approach was then used to analyze a pumping-recovery test conducted near the Platte River in southeastern Nebraska, USA. The hydraulic conductivities and specific yield were calculated for the aquitard and aquifer. The direct-push technique was used to generate sediment columns; permeameter tests on these columns produced the vertical hydraulic conductivities that are compatible with those obtained from the pumping-recovery test. Thus, the combination of the direct-push technique with permeameter tests provides a new method for estimation of vertical hydraulic conductivity. The hydraulic conductivity, determined from grain-size analysis, is smaller than the horizontal one but larger than the vertical one determined by the pumping-recovery test.     

Evaluation of unconfined aquifer parameters from flow to partially penetrating wells in Tailan River basin, China

Effective evaluation, management and abstraction of groundwater resources of any aquifer require accurate and reliable estimates of its hydraulic parameters. This study, therefore, looks at the determination of hydraulic parameters of an unconfined aquifer using both analytical and numerical approaches. A long-duration pumping test data obtained from an unconfined aquifer system within the Tailan River basin in Xinjiang Autonomous Region in the northwest of China is used, in this study, to investigate the best method for estimating the parameters of the aquifer. The pumping test was conducted by pumping from a radial collector well and measuring the response in nine observation wells; all the wells used in the test were partially penetrating. Using two well-known tools, namely AquiferTest and MODFLOW, as an aid for the analytical and numerical approaches, respectively, the parameters of the aquifer were determined and their outputs compared. The estimated horizontal hydraulic conductivity, vertical hydraulic conductivity, and specific yield for the analytical approach are 38.1–50.30 m/day, 3.02–9.05 m/day and 0.204–0.339, respectively, while the corresponding numerical estimates are 20.50–35.24 m/day, 0.10–3.40 m/day, and 0.27–0.31, respectively. Comparing the two, the numerical estimates were found to be more representative of the aquifer in the study area since it simulated the groundwater flow conditions of the pumping test in the aquifer system better than the analytical solution.     

Field experiments and numerical simulations of whirlpool foundation pit dewatering

A whirlpool foundation pit is a small-diameter, deep circular pit. Because of its depth and small diameter, a large drawdown is required, and a limited number of wells can be installed inside the pit. During excavation, partially penetrating wells inside and outside the foundation pit have to be installed to lower the water level when the aquifer is too thick. However, partially penetrating wells near partially penetrating curtains cannot be treated by analytical methods. Therefore, it is necessary to use numerical methods to predict dewatering during excavation. Field experiments were performed on whirlpool foundation pit 1880 of Baosteel Group, Shanghai, China, to obtain pumping rates and drawdown, pumping with a single well and two wells in the confined aquifer. The results indicate that the drawdown inside the pit induced by pumping wells outside the foundation pit was small, whereas it was large for pumping wells inside the pit. The pumping wells inside and outside the pit had to be combined to lower the water level. A three-dimensional numerical model was developed to simulate the dewatering process. The hydraulic conductivities of the confined aquifers were inversed by using the pumping tests. Operation schedules were simulated with the corrected model for different combinations of wells inside and outside the pit. The results suggest that different schedules and operation conditions affect drawdown. The monitored results during dewatering indicate that the simulation and field measurements were in agreement. The results can be applied to similar situations.     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Analysis of flow behavior in a landfill with cover soil of low hydraulic conductivity

 This paper presents the results of field tests of hydrologic parameters in a landfill and the results of numerical simulation to find the efficiency of the pumping method to reduce leachate levels in the landfill. The field hydraulic conductivity and storativity of waste and buried cover soils in the landfill are measured by pumping and slug tests. The hydrologic condition inside the landfill is first calibrated using the drawdown-time curve obtained from the pumping test, and the flow behavior of leachate during pumping in the landfill, when various layers of waste and buried cover soil exist, is analyzed through three-dimensional numerical simulation of flow. The results of the field investigation show that the buried cover soil of low hydraulic conductivity forms an impermeable layer preventing the downward flow of leachate and upward flow of landfill gas. The hydraulic conductivities of the pumping test and slug tests were quite close on the same order of magnitude. It was also possible to match the drawdown-time data of the field tests with those of the model using input data close to the hydrologic property obtained from the field tests. The numerical flow analysis showed that pumping was possible up to 120 tons/day for a single well without a drain, while the pumping rate could be increased to 300 tons/day for the same well with the drain. From the vertical section of the flow vector with a horizontal drain, the barrier role of buried cover soil is identified, which was proposed by examining the water contents of the disposed cover soil and waste in the field. Received: 15 May 1998 · Accepted: 4 January 1999     

基于解析法和数值法的非稳定流抽水试验参数反演

含水层的水文地质参数是进行地下水资源计算、地下水污染防控等所必需的基础数据,结合数值模拟技术进行含水层参数反演很有必要.按照1:5万水文地质调查规范在江汉平原仙桃市杨林尾镇复兴水厂不同含水层位开展抽水试验,包括深层含水层单孔抽水试验以及浅层含水层中群孔(2孔)抽水试验.对于单孔抽水试验,应用第1类越流系统井流理论进行参数反演;对于群孔抽水试验,推导了特定综合井函数,并利用特定标准曲线匹配法和直线图解法求解了含水层参数.随后利用FEFLOW软件建立了相应数值模型,拟合了含水层参数.结果表明:浅层含水层的渗透系数变化范围为21.66~54.00 m/d,贮水率变化范围为1.28×10-5~8.00×10-4 m-1;深层含水层渗透系数变化范围为1.27~7.00 m/d,贮水率变化范围为3.90×10-6~5.00×10-6 m-1.对于深层承压含水层而言,越流补给量较大.采用数值模拟方法结合抽水试验数据求参,综合考虑了含水层结构,拟合效果好,所得结果更加可靠.      

Estimation of regional-scale groundwater flow properties in the Bengal Basin of India and Bangladesh

Quantitative evaluation of management strategies for long-term supply of safe groundwater for drinking from the Bengal Basin aquifer (India and Bangladesh) requires estimation of the large-scale hydrogeologic properties that control flow. The Basin consists of a stratified, heterogeneous sequence of sediments with aquitards that may separate aquifers locally, but evidence does not support existence of regional confining units. Considered at a large scale, the Basin may be aptly described as a single aquifer with higher horizontal than vertical hydraulic conductivity. Though data are sparse, estimation of regional-scale aquifer properties is possible from three existing data types: hydraulic heads, ¹⁴C concentrations, and driller logs. Estimation is carried out with inverse groundwater modeling using measured heads, by model calibration using estimated water ages based on ¹⁴C, and by statistical analysis of driller logs. Similar estimates of hydraulic conductivities result from all three data types; a resulting typical value of vertical anisotropy (ratio of horizontal to vertical conductivity) is 10⁴. The vertical anisotropy estimate is supported by simulation of flow through geostatistical fields consistent with driller log data. The high estimated value of vertical anisotropy in hydraulic conductivity indicates that even disconnected aquitards, if numerous, can strongly control the equivalent hydraulic parameters of an aquifer system.     

Regional fracture network permeability using outcrop scale measurements

Estimating the hydraulic properties of fractured aquifers is challenging due to the complexity of structural discontinuities that can generally be measured at a small scale, either in core or in outcrop, but influence groundwater flow over a range of scales. This modeling study uses fracture scanline data obtained from surface bedrock exposures to derive estimates of permeability that can be used to represent the fractured rock matrix within regional scale flow models. The model is developed using PETREL, which traditionally benefits from high resolution data sets obtained during oil and gas exploration, including for example seismic data, and borehole logging data (both lithological and geophysical). The technique consists of interpreting scanline fracture data, and using these data to generate representative Discrete Fracture Network (DFN) models for each field set. The DFN models are then upscaled to provide an effective hydraulic conductivity tensor that represents the fractured rock matrix. For each field site, the upscaled hydraulic conductivities are compared with estimates derived from pumping tests to validate the model. A hydraulic conductivity field is generated for the study region that captures the spatial variability of fracture networks in pseudo-three dimensions from scanline data. Hydraulic conductivities estimated using this approach compare well with those estimated from pumping test data. The study results suggest that such an approach may be feasible for taking small scale fracture data and upscaling these to represent the aquifer matrix hydraulic properties needed for regional groundwater modeling.     

综合利用多种方法确定潜水含水层渗透系数

为确定潜水含水层的渗透系数,施工了一个主抽水井和两个观测井,采用了潜水含水层稳定流完整井多孔、单孔抽水试验的公式法、非稳定流的完整井的半对数直线图解法和基于抽水试验资料处理软件的纽曼模型求参法三种方法进行计算,结果发现各计算方法的结果相差不大,能够相互验证。最后,选取半对数直线图解法的计算结果作为最终计算结果。该研究对潜水含水层渗透系数的计算有一定的参考价值。     

水力层析法与克立金法估算非均质含水层渗透系数场比较

含水层非均质性空间分布特征的识别,是准确模拟地下水流和污染物运移的前提.基于室内非均质含水层砂箱实验,分别利用水力层析法和克立金插值法刻画了非均质含水层渗透系数场.研究结果表明:(1) 水力层析法与克立金法相比,不仅可以更好地刻画非均质含水层渗透系数场,还可以更高精度地预测地下水流过程;(2) 水力层析抽水实验中,通过增加抽水实验组数可以有效地提高水力层析参数反演的精度,但是抽水实验组数增加到一定程度以后,再增加抽水组数不会显著提升参数反演的效果.后续需要进一步研究水力层析抽水实验合适的组数,进一步对抽水井进行优化布设.      

Modeling approaches and strategies for data-scarce aquifers: example of the Dar es Salaam aquifer in Tanzania

Management of groundwater resources can be improved by using groundwater models to perform risk analyses and to improve development strategies, but a lack of extensive basic data often limits the implementation of sophisticated models. Dar es Salaam in Tanzania is an example of a city where increasing groundwater use in a Pleistocene aquifer is causing groundwater-related problems such as saline intrusion along the coastline, lowering of water-table levels, and contamination of pumping wells. The lack of a water-level monitoring network introduces a problem for basic data collection and model calibration and validation. As a replacement, local water-supply wells were used for measuring groundwater depth, and well-top heights were estimated from a regional digital elevation model to recalculate water depths to hydraulic heads. These were used to draw a regional piezometric map. Hydraulic parameters were estimated from short-time pumping tests in the local wells, but variation in hydraulic conductivity was attributed to uncertainty in well characteristics (information often unavailable) and not to aquifer heterogeneity. A MODFLOW model was calibrated with a homogeneous hydraulic conductivity field and a sensitivity analysis between the conductivity and aquifer recharge showed that average annual recharge will likely be in the range 80–100 mm/year.     

Analysis of Steady-State Flow To Multiscreened Wells Under Natural Vertical Flow Conditions

Water levels measured at multiscreened wells in unconfined aquifers may not coincide, in general, with the elevation of the water table. The presence of vertical gradients (as often is the case in recharge areas) or the existence of confining layers may cause the water levels to differ from local hydraulic heads in the aquifer. In these cases, a misinterpretation of water levels may lead to the erroneous conclusion that observed drawdowns are provoked by overpumping. In this paper, we analyze the effect that a natural vertical gradient has on water levels in wells screened over their entire saturated thickness. As one would expect, it is observed that, even without pumping, the water level in the wells lies below the water table. Type curves relating the steady-state drawdown to the vertical gradient and to the hydraulic conductivity anisotropy are presented. These curves were obtained using a groundwater flow numerical model (FREESURF: Neuman and Witherspoon, 1970). The theoretical results are checked with field data from deep wells in the detrital Madrid aquifer. In this particular aquifer, it is observed that the effect of vertical gradients is important both in terms of drawdowns and flow rates.

---

     

Re-evaluating the US Geological Survey’s pumping tests (1967) in the Punjab region of Pakistan for use in groundwater studies

In 1967, the US Geological Survey (USGS) published the results of 141 pumping tests carried out throughout the Pakistani Punjab to establish representative hydraulic parameters of its large aquifer. Many authors have since concluded that the USGS had over-estimated the horizontal hydraulic conductivity (k _r) by 25–100 %, leaving vertical anisotropy and aquifer depth unresolved. No test wells have ever been drilled below 450 m to reach the base of the aquifer, although petroleum explorations mention depths between 1,500 and 4,500 m. After comparison and re-evaluation of all related papers, this study concludes that the USGS interpretation was correct, that its hydraulic values still stand without change, and that the USGS’s applied distance drawdown interpretation is valid to prevent influence of partial penetration on the results. This study also uniquely resolved vertical anisotropy and aquifer thickness by using early- and late-time drawdowns separately and proper scaling of the coordinates, which has often been omitted. With appropriate scaling, all interpretations match the data. The representative hydraulic aquifer values are: k _r?=?65 m/d, vertical anisotropy k _r/k _z?=?25 and aquifer depth 500–1,500 m. The conclusion is that these values can be used, at least as first estimates, for groundwater studies in the Pakistani Punjab.     

基于多源数据的弱透水层水文地质参数反演研究

弱透水层是含水层系统的重要组成部分,其水文地质参数的准确获取一直以来都是研究难点,传统室内试验难以克服应力状态及环境变化对土层参数的影响,因此如何在原位状态下评价其水文地质参数有着重要意义。选取呼和浩特盆地某淤泥层为研究对象,基于抽水试验和溶质运移试验获得的多源数据,采用基于控制体积有限差分法的非结构化网格（Unstructured Grid）方法,建立三维地下水流-溶质运移耦合模型,对弱透水层水文地质参数进行反演。结果显示:模型计算的抽水井中水位和溶质浓度变化规律与实际观测数据能够较好拟合,模型识别出的抽水含水层渗透参数（4.8 m/d）与通过解析法计算得出的结果（4.17 m/d）相近,建立的数值模型符合实际水文地质条件;反演得到淤泥质弱透水层垂向渗透系数为1.2×10?4 m/d、储水率为1.0×10?5 m?1。基于多源数据的参数反演方法可为弱透水层参数研究提供借鉴。     

Interactions of diffuse and focused allogenic recharge in an eogenetic karst aquifer (Florida, USA)

The karstic upper Floridan aquifer in north-central Florida (USA) is recharged by both diffuse and allogenic recharge. To understand how recharged water moves within the aquifer, water levels and specific conductivities were monitored and slug tests were conducted in wells installed in the aquifer surrounding the Santa Fe River Sink and Rise. Results indicate that diffuse recharge does not mix rapidly within the aquifer but instead flows horizontally. Stratification may be aided by the high matrix porosity of the eogenetic karst aquifer. Purging wells for sample collection perturbed conductivity for several days, reflecting mixing of the stratified water and rendering collection of representative samples difficult. Interpretive numerical simulations suggest that diffuse recharge impacts the intrusion of allogenic water from the conduit by increasing hydraulic head in the surrounding aquifer and thereby reducing influx to the aquifer from the conduit. In turn, the increase of head within the conduits affects flow paths of diffuse recharge by moving newly recharged water vertically as the water table rises and falls. This movement may result in a broad vertical zone of dissolution at the water table above the conduit system, with thinner and more focused water-table dissolution at greater distance from the conduit.     

Characteristics of Valadares Island aquifer,Paranaguá coastal plain,Brazil

. Aquifer characterization in terms of geologic–geomorphological considerations, pumping test and hydrochemistry has been carried out to evaluate the groundwater potential of one of the delicate coastal island systems in Parana State, Brazil. The Valadares Island, lying in the Atlantic coastal zone, represents a stabilized beach ridge having a thickness of more than 20 m unconsolidated quartz-rich sand. Field and laboratory investigations suggest that this homogeneous coastal aquifer can partially meet the local domestic demand of fresh water even though it is in hydraulic continuity with the estuarine and tidal systems of Paranaguá bay. Efficiency of the aquifer to hold and release water was determined by conducting a pumping test in one of the partially penetrating wells in the island for 24 h. The physical, chemical and biological analyses of samples from 12 5-cm diameter wells of the area have indicated that groundwater needs purification for drinking purpose. It is advocated that the groundwater potential of the island aquifer is worth considering for sustained but controlled exploitation.     

The dependence of stream depletion by seasonal pumping on various hydraulic characteristics and engineering factors

Compensation pumping is used to alleviate deficiencies in streamflow discharge during dry seasons. Short-term groundwater pumping can use aquifer storage instead of catchment-zone water until the drawdown reaches the edge of the stream. The capacitance is a complex, dimensionless parameter of an aquifer system that defines the delayed effect on streamflow when there is groundwater pumping. This parameter is a function of aquifer hydraulic characteristics, pumping time, and distance between the well and stream edge; the latter can involve stream leakance and vertical leakance of an associated aquitard. Three typical hydraulic cases of combined water systems (major catchment-zone wells close to the stream and compensation pumping wells) were classified depending on their capacitance structure (i.e. the relationship between surface water and groundwater): (1) perfect hydraulic connection between the stream and aquifer; (2) imperfect hydraulic connection between the stream and aquifer; and (3) essentially imperfect hydraulic connection between the stream and the underlying confined aquifer. The impact of various hydraulic characteristics and engineering factors on stream depletion was examined by conceptual and numerical modeling. To predict the suitability and efficiency of a combined water system application, regression tests were undertaken on unit stream depletion and capacitance, and power dependencies were defined.     

Evaluation of a pumping test of the Snake River Plain aquifer using axial-flow numerical modeling

The Snake River Plain aquifer in southeast Idaho is hosted in a thick sequence of layered basalts and interbedded sediments. The degree to which the layering impedes vertical flow has not been well understood, yet is a feature that may exert a substantial control on the movement of contaminants. An axial-flow numerical model, RADFLOW, was calibrated to pumping test data collected by a straddle-packer system deployed at 23 depth intervals in four observation wells to evaluate conceptual models and estimate properties of the Snake River Plain aquifer at the Idaho National Engineering and Environmental Laboratory. A delayed water-table response observed in intervals beneath a sediment interbed was best reproduced with a three-layer simulation. The results demonstrate the hydraulic significance of this interbed as a semi-confining layer. Vertical hydraulic conductivity of the sediment interbed was estimated to be about three orders of magnitude less than vertical hydraulic conductivity of the lower basalt and upper basalt units. The numerical model was capable of representing aquifer conceptual models that could not be represented with any single analytical technique. The model proved to be a useful tool for evaluating alternative conceptual models and estimating aquifer properties in this application. Electronic Publication     

Estimating hydraulic properties from tidal attenuation in the Northern Guam Lens Aquifer, territory of Guam, USA

Tidal-signal attenuations are analyzed to compute hydraulic diffusivities and estimate regional hydraulic conductivities of the Northern Guam Lens Aquifer, Territory of Guam (Pacific Ocean), USA. The results indicate a significant tidal-damping effect at the coastal boundary. Hydraulic diffusivities computed using a simple analytical solution for well responses to tidal forcings near the periphery of the island are two orders of magnitude lower than for wells in the island’s interior. Based on assigned specific yields of ～0.01–0.4, estimated hydraulic conductivities are ～20–800 m/day for peripheral wells, and ～2,000–90,000 m/day for interior wells. The lower conductivity of the peripheral rocks relative to the interior rocks may best be explained by the effects of karst evolution: (1) dissolutional enhancement of horizontal hydraulic conductivity in the interior; (2) case-hardening and concurrent reduction of local hydraulic conductivity in the cliffs and steeply inclined rocks of the periphery; and (3) the stronger influence of higher-conductivity regional-scale features in the interior relative to the periphery. A simple numerical model calibrated with measured water levels and tidal response estimates values for hydraulic conductivity and storage parameters consistent with the analytical solution. The study demonstrates how simple techniques can be useful for characterizing regional aquifer properties.     

Estimation of hydraulic characteristics from electrical resistivity data in coastal aquifers of southern India

The applicability of geophysical techniques has been examined for evaluating aquifer properties like transmissivity and hydraulic conductivity of coastal aquifers, Tuticorin, Tamil Nadu. The pumping test data of 10 wells are interpreted by using forward modelling to obtain the aquifer characteristics in the study area. The available vertical electrical soundings (VES) data in the vicinity of the sites of pumping test have been interpreted; and true resistivity and thickness are determined at each site in the study area. Empirical relationships are established for estimating the hydraulic parameters from the electrical data.