A hydrogeophysical model of the relationship between geoelectric and hydraulic parameters of anisotropic aquifers

An electrical resistivity method has been used to determine aquifer parameters in the Ganga-Yamuna interfluve in northern India. An existing relationship between the geoelectrical and hydraulic parameters has been modified for the case of an anisotropic aquifer. The hydrogeological framework in the upper part of the Ganga-Yamuna interfluve is evaluated by using existing relationships between hydraulic parameters and geoelectrical parameters for alluvial aquifers. On the basis of aquifer geometry, the area has been divided into two hydraulic units: the western Yamuna flood plain and the Ganga flood plain towards the east. The resistivity data collected in parts of the study area are first interpreted in terms of true resistivity and thicknesses of subsurface layers. The electrical parameters (resistivity and thicknesses) are subsequently correlated with the available pumping test data. Distinct correlations between transmissivity and modified transverse resistance are obtained for the two hydraulic units. A four-parameter model consisting of hydraulic conductivity, modified longitudinal resistivity, modified transverse resistance and hydraulic anisotropy is presented for the anisotropic aquifer underlain by conductive fine grained sediments. The model has been validated at a number of locations, where aquifer parameters are known from pumping test data.     

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.     

江汉平原典型含水层水文地质参数反演

含水层参数对于定性分析区域地下水资源评价、数值模拟及预报、开发利用与保护及科学管理具有重要意义,而江汉平原水文地质参数求解的研究却相对较少.2015年江汉平原1:50 000水文地质调查项目完成杨林尾图幅与陆溪口图幅,选择杨林尾图幅和陆溪口图幅中的4个代表性钻孔,分别做了3组抽水试验和一组压水试验.对抽水数据使用Dupuit公式法结合抽水稳定阶段数据求解含水层参数,同时利用Aquifer Test软件中Theis标准曲线法、直线图解法对抽水试验中非稳定条件下抽水数据对含水层参数求解,得到孔隙水含水层渗透系数及弹性给水度;利用压水试验工程规范求解基岩裂隙含水层的渗透系数.探讨了利用非稳定抽水试验条件下求解含水层参数方法的可行性及该方法的优势.计算结果表明:杨林尾镇浅层(20.2~64.55 m)含水层渗透系数为0.075 m/d,弹性给水度为5.8×10-2;深层(138~160 m)含水层渗透系数为9.89 m/d,弹性给水度为2.3×10-5;陆溪口镇浅层(19.4~36.4 m)含水层渗透系数1.26 m/d,弹性给水度为1.1×10-1;基岩渗透系数为0.012 m/d.通过结果对比分析发现对于单孔非稳定抽水试验,对前期水位降深数据筛选分析,同样可以利用非稳定井流理论反演含水层参数,结果比较可靠.      

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

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

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

含水层的水文地质参数是进行地下水资源计算、地下水污染防控等所必需的基础数据,结合数值模拟技术进行含水层参数反演很有必要.按照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 aquifer hydraulic parameters from surface geophysical measurements: a case study of the Upper Cretaceous aquifer, central Sinai, Egypt

The integration of geophysical data with direct hydrogeological measurements can provide a minimally invasive approach to characterize the subsurface at a variety of resolutions and over many spatial scales. The field of hydrogeophysics has attracted much attention during the last two decades. In this domain, the geophysical data inverted to geophysical models are interpreted in terms of the hydrogeology to serve as a basis for the definition of hydraulic models in the areas of interest. The hydraulic conductivity (K) value measured in a reference borehole has been combined with the electrical conductivity obtained from nearby geo-electromagnetic sounding data in the Cenomanian (Upper Cretaceous) aquifer, central Sinai, Egypt. The resulting relation was interpreted with Dar Zarrouk parameters to infer the transmissivity variations at other vertical electrical sounding locations, where K values are unknown. Coincident transient electromagnetic data have been adopted to increase accuracy while interpreting the aquifer geoelectrical properties. The results indicate that the transmissivity values in the aquifer of interest vary from 2,446 to 9,694 m²/day, and K varies from 12.9 to 57.0 m/day throughout the studied area.     

A comparative study of hydraulic conductivity estimations using geostatistics

Three approaches for estimating the hydraulic conductivity (K) of the Trifa aquifer, Morocco were investigated: (1) kriging of the K values obtained from pumping tests, (2) cokriging of the pumping test data with electrical resistivity data as a secondary variable, and (3) cokriging of the pumping test data with the slope of the water table. Gauss-transformed values of the variables are used because they provide more robust variograms and transformed values of the primary and secondary variables show correlations higher than the raw values, which is beneficial in cokriging. In cokriging with electrical resistivity, two zones are considered since the geological deposits are different from the north to the south of the aquifer, which is reflected in different correlations between the variables. Comparison of the three approaches is based mainly on the estimation errors, and to a lesser degree on the cross-validations of the corresponding variogram models and general considerations, like the measurements’ reliability and aquifer make-up. The best-estimated K is given by cokriging with the slope of the water table and is therefore preferred for further use in groundwater flow modeling. Thus, electrical resistivity or the slope of the water table can both be used as secondary variables to estimate K, especially in heterogeneous aquifers with lateral variations in lithology, as is the case of the Trifa aquifer.     

Vertical movement of water in a High Plains Aquifer induced by a pumping well

Field observation and numerical simulations were carried out to evaluate the hydraulic relationship between the shallow and deep aquifer of a High Plains Aquifer system, in which shallow and deep aquifers are separated by an aquitard. Pumping from the lower aquifer resulted in a small drawdown in the upper aquifer and a larger drawdown in the aquitard; pumping from the shallow aquifer caused a small drawdown in the aquitard and the deep aquifer. Analysis of pumping test data gives the values of the hydraulic conductivity of the aquitard and the deep aquifer. Long-term observation of groundwater levels in the shallow and deep aquifers showed that a strong downward hydraulic gradient was maintained during an irrigation season. Numerical simulations were used to calculate the induced leakage of water from the shallow to the deep aquifer. Water budget analyses suggested that after pumping continues for a couple of days, the leakage from the overlying layers begins to supply the majority of the withdrawal from the deep aquifer. However, the induced leakage from the upper shallow aquifer can travel only a few meters into the aquitard, and it can not reach the lower aquifer during a 90 day pumping period. The major portion of the induced leakage occurred during the pumping period, but a small leakage can continue as a residual effect after the pumping period. The vertical hydraulic conductivity of the aquitard plays a major role in partitioning the ratio of the induced leakage for the pumping and after-pumping periods.     

利用滨海地下水系统对双频潮汐的响应识别含水层压力传导系数

含水层的压力传导系数是地下水溶质运移及扩散的重要参数,通常压力传导系数由抽水试验确定,即通过抽水给地下水系统一个扰动,然后监测地下水水位的响应,进而计算压力传导系数.对滨海地区含水层来说,潮汐波动就是天然的抽(注)水试验.潮汐由不同频率的分潮构成,前两项主频潮汐基本可代表潮汐主要特征.利用地下水水位对潮汐的响应,识别含...     

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.     

Hydrogeological characterization of aquifer in palla flood plain of Delhi using integrated approach

The Yamuna flood plains spread across the northern part of Indian subcontinent is home to millions of people. The ever-growing population in these plains make it difficult to sustain the demand of freshwater resources. However, the productive aquifers of flood plains could provide solution for these issues. In this context, it is necessary to understand the aquifer characteristics. Thus, the paper attempts to characterize the aquifer in Palla area of the flood plain using integrated approach. Besides, grain size analysis and site-litholog study, the nature of aquifer material was also ascertained from bulk mineralogy of the sediments using X-Ray Diffraction. The aquifer parameters were estimated with help of long duration pumping test data. Moreover, the effect of pumping on salinity variation and hydrochemical facies evolution was also examined. The sand dominant, unconfined aquifer was estimated to have horizontal hydraulic conductivity in the range of 25 m/day and vertical hydraulic conductivity of 6–7 m /day. While the specific yield of the aquifer was estimated in the range of 0.07–0.08. It is observed that under conducive active flood plain environment, the given sand mineralogy at the site does not allow salinity increase in groundwater even after more than a decade of groundwater pumping. In fact, over years, hydrochemical facies have evolved towards Bicarbonate type. These things put together make the active flood plain aquifer a sustainable groundwater resource.     

Hydrogeologic characteristics of the alluvial tuff aquifer of northern Sahand Mountain slopes, Tabriz, Iran

The Tabriz area is a densely populated area of northwestern Iran (more than 1.5 million in population) with a large proportion of its drinking, domestic, industrial and agricultural water supplied from groundwater resources. The average rate of drinking and industrial water use in the city of Tabriz is about 3.45 m³ s^–1. The Plio-Pleistocene unconfined alluvial tuff aquifer (about 1,275 km²), the most important aquifer in the area, has been known for many years as a reliable resource. The greatest estimated thickness of the alluvial tuff lies in the Saidabad area, with 350 m thickness. There are 994 deep and 284 shallow active pumping wells and 83 qanats operate in the alluvial tuff aquifer. The total water withdrawal from all these artificial discharge points has been measured at 72, 3.8 and 17 million m3/year, respectively. Analytical and numerical methods have been applied to the constant rate pumping test data from the Saidabad wellfield (eight pumping and three observation wells). The values of electrical conductivity in the groundwater of alluvial tuff aquifer range from 203 to 960 μS cm^–1 and bicarbonate type water dominates.     

基于广义径向流模型的非均质孔隙含水层井流试验分析

Theis模型、Dupuit模型等经典井流模型分析非均质含水层井流试验数据有一定的局限性,获取的参数不适合表征非均质含水层特性。而GRF模型可以获取含水层流动特性的数据。相比于Theis模型、Dupuit模型,GRF模型更能表征非均质含水层特性。以黄石东湖新村棋盘洲长江大桥的抽水试验数据为例,采用Theis模型和GRF模型计算含水层渗透系数,结合实际水文地质条件,对比分析不同方法计算的水文地质参数,并计算含水层水流维数和表观压力传导系数（K_f /S_sf）。结果表明:研究区含水层为细砂夹条带状黏土透镜体的非均质含水层,采用GRF模型计算结果更符合实际情况,渗透系数为 4.09×10?3cm/s;含水层水流维数为1.61,地下水为双线性流动状态,含水层对抽水试验的响应主要受黏土条带控制;观测井和抽水主井的K_f /S_sf呈非线性相关,进一步验证了含水层的非均质性。在非均质孔隙含水层中,应用多孔联合非稳定GRF井流试验方法不仅能确定水文地质参数,并且能丰富对含水层特性的认知。     

Development of a simple method for determining the influence radius of a pumping well in steady-state condition

Influence radius of a pumping well is a crucial parameter for hydrogeologists and engineers. Knowing the radius of influence for a designed drawdown enables one to calculate the pumping rate required to layout a project foundation that may need lowering of groundwater level to a certain depth due to dewatering operation. In addition, this is important for hydrogeologists to determine ground water contamination flow paths and contributing recharge area for domestic water supply and aquifer management purposes. Empirical formulas that usually neglect vital parameters to determine the influence radius accurately have been traditionally utilized due to lack of adequate methods. In this study, a physically based method, which incorporates aquifer hydraulic gradient for determining the influence radius of a pumping well in steady-state flow condition, was developed. It utilizes Darcy and Dupuit laws to calculate the influence radius, where Darcy's law and Dupuit equation, in steady-state condition, represent the inflow and the outflow of the pumping well, respectively. In an untraditional manner, this method can be also used to determine aquifer hydraulic conductivity as an alternative to other pumping test methods with high degree of accuracy. The developed method is easy to use; where a simple mathematical calculator may be used to calculate the influence radius and the pumping rate or hydraulic conductivity. By comparing the results from this method with the MODFLOW numerical model outputs with different simulated scenarios, it is realized that this method is much superior and more advantageous than other commonly used empirical methods.     

A numerical investigation to illustrate the consequences of hydraulic connections between granular and fractured-rock aquifers

Natural or artificial hydraulic connections between a granular aquifer in contact with a fractured-rock aquifer can have significant physical and chemical impacts at both a local and a regional scale. In this study, numerical simulations are conducted in order to illustrate the hydrogeological consequences of such hydraulic relationships. The numerical investigation, based on a conceptual model, focuses on the effects of the hydraulic connections when conducting a pumping test in a well that is set into a granular confined aquifer overlying a fractured-rock aquifer which presents a few fractures directly connected to the granular aquifer. It is illustrated that when interpreting the pumping test with the conventional methods consisting of plotting the drawdown versus time, a bias is introduced on the estimation of the transmissivity of the granular aquifer due to groundwater flowing from the fractured-rock aquifer via connecting fractures. However, it is underlined that plotting drawdown log-derivative versus time helps to diagnose the existence of these hydraulic relationships and therefore avoids committing a bias on the transmissivity estimation of the granular aquifer. Numerical simulation results also illustrate that hydraulic connections between the two aquifers can have significant impacts on the hydrogeochemical signature of the granular aquifer under investigation.     

抽水含水层对非抽水含水层水位影响的效应分析

利用自动水位计记录的抽水试验过程中不同层位含水组之间的动态水位,分析了在抽水试验最初阶段时抽水含水组与非抽水含水组之间的水力联系。研究结果表明:(1)抽水含水层对非抽水含水层存在水位影响效应,在抽水初期,水位影响效应较大;当抽水含水层水位达到稳定或缓慢下降时,水位影响效应随之消失。(2)单位降深水位影响效应显示,上部含水层距离抽水含水层越远,水位影响效应越大;下部含水层距离抽水含水层越远,水位影响效应越小。(3)利用三维地下水流数值模拟方法和应力应变对含水层水位动态影响的原理对水位影响效应机理进行分析,发现水位影响效应与含水层的固体骨架压缩系数(α)、水体积压缩系数(β)、孔隙度(n)和边界控制系数(Cm)等参数有关。     

Quantitative analysis of springs

Growing demand for groundwater resources and stringent environmental concerns has led to large groundwater investigations, including characterization of aquifer systems that are hydraulically connected to springs. A pumping test is one of the most reliable means of quantifying hydraulic characteristics and the response of natural springs to pumping because it yields results that, in general, are representative of a larger area than are results from a single point observation. Recharge to the aquifer sustaining discharge from springs must be evaluated prior to the utilization of springs. The spring hydrograph is analyzed, as the shape of a hydrograph is a reflection of the response of the aquifer to recharge. The form and rate of recession provide significant information on the storage, lithological composition, and structural characteristics of the aquifer system sustaining the spring. Water tracing techniques have been developed and used over a period of centuries to delineate catchment boundaries, estimate groundwater flow velocities, determine areas of recharge, and identify sources of pollution of spring water.     

Inverse modelling of aquifer parameters in basaltic rock with the help of pumping test method using MODFLOW software

The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems. As groundwater flows through an aquifer, its composition and temperature may variation dependent on the aquifer condition through which it flows. Thus, hydrologic investigations can also provide useful information about the subsurface geology of a region. But because such studies investigate processes that follow under the Earth's shallow, obtaining the information necessary to answer these questions is not continuously easy. Springs, which discharge groundwater table directly, afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity (T) and storativity (S) are vital for the evaluation of groundwater resources. There are several methods to estimate the accurate aquifer parameters (i.e. hydrograph analysis, pumping test, etc.). In initial days, these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory. The simultaneous information on the hydraulic behavior of the well (borehole) that provides on this method, the reservoir and the reservoir boundaries, are important for efficient aquifer and well data management and analysis. The most common in-situ test is pumping test performed on wells, which involves the measurement of the fall and increase of groundwater level with respect to time. The alteration in groundwater level (drawdown/recovery) is caused due to pumping of water from the well. Theis (1935) was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer. It is essential to know the transmissivity (T = Kb, where b is the aquifer thickness; pumping flow rate, Q = TW (dh/dl) flow through an aquifer) and storativity (confined aquifer: S = bS_s, unconfined: S = S_y), for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer. The determination of aquifer's parameters is an important basis for groundwater resources evaluation, numerical simulation, development and protection as well as scientific management. For determining aquifer's parameters, pumping test is a main method. A case study shows that these techniques have been fast speed and high correctness. The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.     

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.