Impact of multi-purpose aquifer utilisation on a variable-density groundwater flow system in the Gippsland Basin, Australia

The Latrobe aquifer in the Gippsland Basin in southeastern Australia is a prime example for emerging resource conflicts in Australian sedimentary basins. The Latrobe Group forms a major freshwater aquifer in the onshore Gippsland Basin, and is an important reservoir for oil and gas in both onshore and offshore parts of the basin. The Latrobe Group and overlying formations contain substantial coal resources that are being mined in the onshore part of the basin. These may have coal-seam-gas potential and, in addition, the basin is considered prospective for its geothermal energy and CO₂ storage potential. The impacts of groundwater extraction related to coal-mine dewatering, public water supply, and petroleum production on the flow of variable-density formation water has been assessed using freshwater hydraulic heads and impelling force vectors. Groundwater flows from the northern and western edges towards the central part of the basin. Groundwater discharge occurs mainly offshore along the southern margin. Post-stress hydraulic heads show significant declines near the petroleum fields and in the coal mining areas. A hydrodynamic model of the Latrobe aquifer was used to simulate groundwater recovery in the Latrobe aquifer from different scenarios of cessation of groundwater and other fluid extractions.     

水力联系系数法定量评价含水层之间水力联系

为了定量研究与评价含水层之间地下水水力联系程度,首次提出了水力联系系数C（hydraulic connection coefficient）的概念。将水力联系系数C定义为观测孔目的含水层水位降深与该观测孔位置抽水含水层水位降深的比值。通过水力联系系数,可定量评价某含水层水平上同层之间和垂向上不同含水层之间的水力联系程度。依据鄂尔多斯盆地白垩系洛河组各含水层段的水力联系系数C值,将含水层之间水力联系分为5个等级。其中,0.000 0 ≤ C<0.062 5,0.062 5 ≤ C<0.125 0,0.125 0 ≤ C<0.250 0,0.250 0 ≤ C<0.500 0,C ≥ 0.500 0,分别代表水力联系等级为极弱、弱、中等、强、极强。以高家堡井田钻孔抽（放）水试验数据为例,采用水力联系系数和观测孔水位响应时间两个指标定量评价了区内巨厚层状非均质洛河组含水层内部水力联系。结果表明：洛河组中上段水平同层之间的水力联系系数分别为0.373 0、0.413 8,观测孔水位响应时间较短（约为5 min）,水力联系强;洛河组下段水平同层之间水力联系系数分别为0.440 1、0.491 1,观测孔水位响应时间较短（为9~20 min）,水力联系强;洛河组中上段与下段垂向水力联系系数分别为0.000 2、0.007 2、0.089 7,观测孔水位响应时间较长（大于60 min）,水力联系极弱至弱。     

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.     

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.     

Influence of Hydraulic Conductivity and Wellbore Design in the Fate and Transport of Nitrate in Multi-aquifer Systems

Nitrate concentrations in multi-aquifer systems are heavily affected by the presence of wellbores (active or abandoned) that are screened in several aquifers. The spatial variability of hydraulic conductivity in the confining layers has also an important impact on the concentrations. A synthetic three-dimensional flow and transport exercise was carried in a multi-aquifer system consisting of two aquifers separated by an aquitard in which 100 vertical wellbores had been drilled. To model the wellbores and the flow and transport connection between aquifers that they may induce, we assign a high vertical hydraulic conductivity and a low effective porosity to the cell blocks including the wells. With these parameters, a solute will travel quickly from one aquifer to the other without being stored in the well itself. The wellbores will act as preferential pathways, and the solute will move quickly between aquifers according to the hydrodynamic conditions. Not considering these preferential pathways could induce erroneous interpretations of the solute distribution in an aquifer. We also noted that when there are vertical wellbores that connect aquifers in a multi-aquifer system, low conductivity in the aquitard enhances the flow of solute through the wellbores. Time-varying pumping rates induce important fluctuations in nitrate concentrations; therefore, any estimate of the water quality of the aquifer will depend on the moment when the data has been recorded. Consequently, concentration maps obtained by interpolation of point samples are seldom a good indicator of the chemical status of groundwater bodies; alternatively, we recommend complementing the usual interpolated maps with numerical models to gain a true understanding of the spatial distribution of the solute concentration.     

松藻矿区地下水动力场特征及其对煤层气富集的影响

松藻矿区是我国煤层气资源富集区之一。为揭示松藻矿区地下水对煤层气富集影响,利用煤田勘查获取的水文地质资料,基于含水层结构、水化学特征、抽水资料的系统分析,阐明了含煤地层上覆和下伏含水层地下水的空间展布特征,建立了地下水的分布模式,探讨了地下水动力场对煤层气富集的影响。研究表明:松藻矿区含水层和隔水层相互叠置,随着埋深的增加含水性减弱;含煤岩系与上覆、下伏含水层缺乏水力联系,为相互独立的含水系统;上覆长兴组和下伏茅口组地下水动力场呈南高北低展布;地下水动力场有近河流岩溶裂隙渗入型、断层沟通岩溶裂隙补给型、稳定承压弱含水型三种类型,以稳定承压弱含水型为主;含煤岩系的极弱含水性及其与上覆、下伏含水地层缺乏水力联系构成了区内水文地质控气基础,是研究区煤层含气性高的重要控制因素。      

The distribution of saline groundwater and its relation to the hydraulic conditions of aquifers and aquitards: examples from Israel

The effect of separation by aquitard layers on the distribution of saline groundwater in coastal aquifers has been demonstrated in two Israeli coastal aquifers: the Mediterranean and the Dead Sea aquifers. There is vertical separation in the Dead Sea area, even where the clayey aquitard layers are <1?m thick, exhibited by large differences in hydraulic head (2?C5?m), salinity (TDS of 50?C340?g/L) and chemical composition (e.g. Na/Cl range 0.28?C0.55). Similar features are found in the Mediterranean coastal aquifer, where the separating aquitard layers are thicker (??5?C10?m). Here, the different subaquifers host fresh and saline groundwater of different ages (tritium and ¹⁴C ages range from tens to thousands of years), as well as different chemical compositions. This high resolution of results can be obtained only by drilling without fluids; otherwise, the spatial information may lead to incorrect representation of the studied aquifer. This is especially important in saline systems where only partial flushing occurs and, thus, large variations in salinity and chemical composition are expected. The main factors controlling the salinity of groundwater in subaquifers in coastal aquifers are their connection to the sea or saline lakes, existence of brines, salinization and flushing rates, and separation by aquitard layers.     

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

In groundwater hydrology, aquitard heterogeneity is often less considered compared to aquifers, despite its significant impact on groundwater hydraulics and groundwater resources evaluation. A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity (K) and specific storage (Ss). The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or Ss. Two scenarios of analytical investigations are explored: one is the presence of aquitard interlayers with distinct K or Ss values, a common field-scale occurrence; another is an exponentially depth-decaying aquitard Ss, a regional-scale phenomenon supported by statistical analysis. Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion; a high-Ss interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion. Locations of low-K or high-Ss interlayers also significantly impact well hydraulics and groundwater budget. In the context of an exponentially depth-decaying aquitard Ss, a larger decay exponent can enhance aquifer drawdown. When using current models with a vertically homogeneous aquitard, half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard Ss should be used to calculate aquitard depletion and unconfined aquifer leakage.     

Geochemical indicators of interbasin groundwater flow within the southern Rio Grande Valley, southwestern USA

Increased groundwater withdrawals for the growing population in the Rio Grande Valley and likely alteration of recharge to local aquifers with climate change necessitates an understanding of the groundwater connection between the Jornada del Muerto Basin and the adjoining and more heavily used aquifer in the Mesilla Basin. Separating the Jornada and Mesilla aquifers is a buried bedrock high from Tertiary intrusions. This bedrock high or divide restricts and/or retards interbasin flow from the Jornada aquifer into the Mesilla aquifer. The potentiometric surface of the southern Jornada aquifer near part of the bedrock high indicates a flow direction away from the divide because of a previously identified damming effect, but a groundwater outlet from the southern Jornada aquifer is necessary to balance inputs from the overall Jornada aquifer. Differences in geochemical constituents (major ions, δD, δ¹⁸O, δ³⁴S, and ⁸⁷Sr/⁸⁶Sr) indicate a deeper connection between the two aquifers through the Tertiary intrusions where Jornada water is geochemically altered because of a geothermal influence. Jornada groundwater likely is migrating through the bedrock high in deeper pathways formed by faults of the Jornada Fault Zone, in addition to Jornada water that overtops the bedrock high as previously identified as the only connection between the two aquifers. Increased groundwater withdrawals and lowering of the potentiometric surface of the Jornada aquifer may alter this contribution ratio with less overtopping of the bedrock high and a continued deeper flowpath contribution that could potentially increase salinity values in the Mesilla Basin near the divide.     

Mapping an aquitard breach using shear-wave seismic reflection

In multi-layered hydrostratigraphic systems, aquitard breaches caused by faulting or paleo-erosion can allow substantial quantities of water of differing quality to be exchanged between aquifers. Seismic reflection technology was used to map the extent and orientation of an aquitard breach connecting a shallow alluvial aquifer to the deeper semi-confined Memphis aquifer in southwestern Tennessee, USA. Geophysical well logs indicate the presence of the aquitard at borehole locations that define the beginning and end points on two seismic survey lines, which intersect at a borehole where the aquitard is absent. A SE–NW-oriented paleochannel, 350 m wide and approximately 35–40 m deep, is interpreted from the seismic reflection surveys. The paleochannel cuts through the aquitard and into the upper part of the Memphis aquifer, thus creating a hydraulic connection between the shallow unconfined and deeper, semi-confined aquifers. The results indicate the potential of the shear-wave seismic reflection methods to resolve shallow breaches through fine-grained aquitards given availability of sufficient well control.     

Regional hydrostratigraphy and groundwater flow modeling in the arsenic-affected areas of the western Bengal basin, West Bengal, India

The first documented interpretation of the regional-scale hydrostratigraphy and groundwater flow is presented for a ～21,000-km² area of the arsenic-affected districts of West Bengal [Murshidabad, Nadia, North 24 Parganas and South 24 Parganas (including Calcutta)], India. A hydrostratigraphic model demonstrates the presence of a continuous, semi-confined sand aquifer underlain by a thick clay aquitard. The aquifer thickens toward the east and south. In the south, discontinuous clay layers locally divide the near-surface aquifer into several deeper, laterally connected, confined aquifers. Eight 22-layer model scenarios of regional groundwater flow were developed based on the observed topography, seasonal conditions, and inferred hydrostratigraphy. The models suggest the existence of seasonally variable, regional, north–south flow across the basin prior to the onset of extensive pumping in the 1970s. Pumping has severely distorted the flow pattern, inducing high vertical hydraulic gradients across wide cones of depression. Pumping has also increased total recharge (including irrigational return flow), inflow from rivers, and sea water intrusion. Consequently, downward flow of arsenic contaminated shallow groundwater appears to have resulted in contamination of previously safe aquifers by a combination of mechanical mixing and changes in chemical equilibrium.     

云应盆地东北部含水层结构特征及地下水转化模式

云应盆地东北部属鄂北贫水地区,赋存于古近系—第四系含水层中的地下水是当地生产、生活用水的主要来源,亟需查明含水层的结构、含水层间地下水的转化关系等基本条件,为研究区内合理开发利用地下水提供依据。本研究通过野外水文地质调查、水文地质钻探工作,将研究区划分为单层含水层与双层含水层结构两个亚区（6个小区）。并通过地下水水位动态长期监测,获取了区内不同含水层的水位动态变化特征,分析各含水层之间的水力联系,建立了区域地下水转化的概念模式,即:研究区地下水以接受山前降雨入渗及风化裂隙水侧向径流补给为主,主要以水平径流的形式经古近系孔隙-裂隙含水层及第四系孔隙承压含水层往澴水方向运移,而后进入第四系孔隙潜水含水层。地下水和地表水在不同季节补排模式不同,雨季地表水（澴水）补给地下水,旱季地下水向地表水（澴水）排泄。古近系孔隙-裂隙水与上覆第四系孔隙水联系密切互为补给,共同构成具有统一水力联系的垂向多层结构的含水系统。独特的含水层结构决定了区内地下水接受降水补给的条件较差,地下水可开采资源量总体较贫乏,建议重点利用区域地表水资源,适度开发地下水资源,推进农业节水灌溉工程,实现水资源可持续利用。     

Use of mathematical modeling to investigate inter-aquifer contamination by organic-rich water through an unplugged well (central Wielkopolska, Poland)

Hydrogeological research was performed to establish the origin of organic-rich water in a Quaternary aquifer (Wielkopolska Buried Valley aquifer, Poland). The spatial distribution of such water was determined by means of multilevel sampling with a packer. The nature, spatial distribution and chemistry of the organic-rich water suggest hydraulic connection between the Quaternary aquifer and the underlying Neogene aquifer. This connection may be enabled by an old abandoned, improperly plugged well, and must therefore be considered as an artificial hydraulic connection between hydrogeologically separated aquifers. For the verification of this hypothesis, a groundwater contaminant transport model was constructed. The results of contaminant transport modeling allow the contaminant transport parameters (transverse and longitudinal dispersivity) to be identified in the field. The numerical groundwater model was used for the prognosis of organic-rich water remediation, and for specification of water-resource management principles in the region where organic-rich water occurs. The results of this research enable the verification of assumptions regarding complex groundwater flow conditions in the boundary zone of a buried valley.     

How basin model results enable the study of multi-layer aquifer response to pumping: the Paris Basin, France

The Albian aquifer of the Paris Basin (France) has been exploited since 1841 and shows drastic drawdown. A three-dimensional (3D) groundwater flow model is used to study the hydrodynamic response of the multi-layered aquifers to pumping activity in the Albian, at basin scale over 167 years. This 3D flow model uses geometry and hydrodynamic parameter distributions that are inherited from a genetic approach through basin modelling, the basin model creating a geometric pattern of hydrodynamic properties constrained by geological history. The paper aims to promote the use of the basin model approach (long time scale, 248 Ma) for the study of deep-aquifer response to anthropogenic perturbation (short time scale, 167 years) in situations for which hydrodynamic data are scarce but geological data are numerous. The results show that parameter distribution is insufficient to reproduce the Albian aquifer behaviour, notably highlighting a different meaning of the specific storage coefficient between basin modelling and groundwater-flow modelling. Dividing the storage coefficient by 100 and including available transmissivity data significantly improved the model/data comparison. The potential impact on a deep aquitard is then discussed. This study sheds light on the advantages and limitations of the basin model approach for groundwater-flow modelling in 3D.     

地面沉降对含水层参数及给水能力的影响研究

地面沉降主要由深层地下水开采造成,含水层的压缩释水是深层地下水开采量主要构成来源,同时随其压缩变形孔隙比减小造成储水系数、渗透系数的减小,对弱透水层非弹性释水量、越流补给量造成影响。地面沉降过程中的水文地质参数非线性变化及对承压含水层系统的反馈作用也成为水文地质领域的前沿问题。为对地下水开采量、沉降量、地质参数变化以及给水能力变化之间的关系做一个较为定量定性的探究,以含水层压缩过程中的物理机制为依据,并基于沧州地区深层承压含水层的地下水位变化过程和水文地质参数,采用数学手段构建储水系数和渗透系数变化的一维非线性沉降模型。模拟结果显示随承压水头下降,储水系数最高可减小77%,含水层系统的给水能力和储水能力会随压缩变形减小50%甚至更多。研究成果为深入认识水文地质参数与应力变化相关关系、科学评价承压含水层地下水储水调节能力有重要借鉴意义。     

The relationship between the water table and the surface flow of a losing stream, lower Medano Creek, Great Sand Dunes National Monument, Colorado

 Proposed groundwater withdrawals in the San Luis Valley of Colorado may lower the water table in Great Sand Dunes National Monument. In response, the National Park Service initiated a study that has produced a generalized conceptual model of the hydrologic system in order to assess whether a lowering of the water table might decrease the surface flow of lower Medano Creek. Based upon information obtained during the drilling of several boreholes, there appear to be five important hydrostratigraphic units underlying lower Medano Creek within the upper 30 m of the ground surface: 1. a perched aquifer overlying an aquitard located between about 5 and 6 m below the ground surface; 2. the aquitard itself; 3. an unconfined aquifer located between the upper and lower aquitards; 4. an aquitard located between about 27 and 29 m below the ground surface; and 5. a confined underlying the lower aquitard. Because the areal extent of the aquitards cannot be determined from the borehole data, a detailed conceptual model of the hydrogeologic system underlying lower Medano Creek cannot be developed. However, a generalized conceptual model can be envisioned that consists of a complex system of interlayered aquifers and leaky aquitards, with each aquifer having a unique hydraulic head. Water levels in the perched aquifer rise rapidly to their annual maximum levels in response to the arrival of the flow terminus of Medano Creek during the spring runoff event, and the location of the flow terminus is directly dependent upon the discharge of the creek. Water levels in the deeper, non-perched aquifers do not appear to fluctuate significantly in response to the arrival of the flow terminus, demonstrating that it is unlikely that the proposed groundwater withdrawals will decrease the surface flow of lower Medano Creek. Received: 27 December 1995 · Accepted: 20 February 1996     

孟加拉国巴拉普库利亚煤矿含水层水力联系研究

为进一步查明孟加拉国巴拉普库利亚煤矿水文地质特征及含水层间水力联系,为矿井水害防治提供理论依据,以19组水质化验数据为基础,结合含水层及隔水层空间展布特征、井田构造特征、水位历时曲线、水化学类型、氢氧同位素特征等,综合分析新近系UDT含水层、Ⅵ煤顶板含水层、Ⅵ煤含水层之间的水力联系。结果表明,井田北部LDT隔水层局部缺失,为UDT含水层水向含煤地层补给提供了条件;井田北翼煤层顶板含水层与UDT含水层水位变化规律密切相关,且水位相近,初步证明两者存在水力联系;各含水层水均为HCO₃-Ca·Na·Mg型,均为低矿化度水,进一步证明各含水层间存在水力循环;聚类分析结果表明各含水层水质存在一定关联度,推演含水层间水力联系程度;氢(δD)氧(δ18O)同位素特征点分布于全球大气降水线附近,表明大气降水是各含水层共同的补给水源。研究成果可以指导孟加拉国巴拉普库利亚煤矿Ⅵ煤层开采时水害防治方向。      

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.     

Impact of groundwater withdrawals on the interaction of multi-layered aquifers in the Viterbo geothermal area (central Italy)

The impact of groundwater withdrawals on the interaction between multi-layered aquifers with different water qualities in the Viterbo geothermal area (central Italy) was studied. In this area, deep thermal waters are used to supply thermal spas and public pools. A shallow overlying aquifer carries cold and fresh water, used for irrigation and the local drinking-water supply. Starting with a conceptual hydrogeological model, two simplified numerical models were implemented: a steady-state flow model of the entire groundwater system, and a steady-state flow and heat transport model of a representative area, which included complex interactions between the aquifers. The impact of increased withdrawals associated with potential future development of the thermal aquifer must be considered in terms of the water temperature of the existing thermal sources. However, withdrawals from the shallow aquifer might also influence the discharge of thermal sources and quality of the water withdrawn from the shallow wells. The exploitation of the two aquifers is dependent on the hydraulic conductivity and thickness of the intervening aquitard, which maintains the delicate hydrogeological equilibrium. Effective methods to control this equilibrium include monitoring the vertical gradient between the two aquifers and the residual discharge of natural thermal springs.     

含水层结构和地下水面形态的虚拟模拟——以苏锡常地区为例

利用Visual C++作为虚拟苏锡常地区第四系含水层结构和地下水面形态的开发工具,采用OpenGL对它们进行三维可视化表达,研究者可采用不同的显示方式对含水层和地下水面进行观察,了解各地层间的空间分布特征、含水层间的补给关系、地下水面的形态和水平方向上地下水补给区和排泄区的分布情况,以及地下水面的动态演示过程。