稳定井流抽水含水层参数计算及富水性评价

针对目前单孔稳定流求参存在的问题,本文在分析新生界松散含水层条件及三次降深抽水过程基础上,利用其抽水试验恢复阶段的数据,分别求得各含水层多个参数,其值真实反映了含水层的实际情况。利用多元回归方法,求得降深与流量关系,通过其系数值大小分析,间接得出各含水层的富水性程度,为地下水的勘探与评价提供一定借鉴。     

Effect of wellbore storage and finite thickness skin on flow to a partially penetrating well in a phreatic aquifer

An analytical model is presented for the analysis of constant flux tests conducted in a phreatic aquifer having a partially penetrating well with a finite thickness skin. The solution is derived in the Laplace transform domain for the drawdown in the pumping well, skin and formation regions. The time-domain solution in terms of the aquifer drawdown is then obtained from the numerical inversion of the Laplace transform and presented as dimensionless drawdown–time curves. The derived solution is used to investigate the effects of the hydraulic conductivity contrast between the skin and formation, in addition to wellbore storage, skin thickness, delayed yield, partial penetration and distance to the observation well. The results of the developed solution were compared with those from an existing solution for the case of an infinitesimally thin skin. The latter solution can never approximate that for the developed finite skin. Dimensionless drawdown–time curves were compared with the other published results for a confined aquifer. Positive skin effects are reflected in the early time and disappear in the intermediate and late time aquifer responses. But in the case of negative skin this is reversed and the negative skin also tends to disguise the wellbore storage effect. A thick negative skin lowers the overall drawdown in the aquifer and leads to more persistent delayed drainage. Partial penetration increases the drawdown in the case of a positive skin; however its effect is masked by the negative skin. The influence of a negative skin is pronounced over a broad range of radial distances. At distant observation points the influence of a positive skin is too small to be reflected in early and intermediate time pumping test data and consequently the type curve takes its asymptotic form.     

Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases,southern Tunisia: implications of return flow

Abstract

Major ions and stable isotopes in groundwaters of the Plio-Quaternary shallow aquifer of the Djerid oases, southern Tunisia, were investigated to elucidate the origin of groundwater recharge and the mineralization processes. It has been demonstrated that the groundwater composition is mainly controlled by the water–rock interaction, the encroachment of brines from the Chotts and the return flow of irrigation waters. The isotopically depleted groundwater samples suggest that the recharge waters derive from an old palaeoclimatic origin. However, the enriched groundwater samples reflect the presence of evaporated recharge water. Furthermore, the large negative deuterium-excess values indicate the effect of secondary evaporation processes, probably related to the return flow of irrigation waters pumped from the underlying aquifer.

Editor D. Koutsoyiannis; Associate editor E. Custodio

Citation Tarki, M., Dassi, L. and Jedoui, Y., 2012. Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases, southern Tunisia: implications of return flow. Hydrological Sciences Journal, 57 (4), 790–804.     

Integrated assessment of lateral flow, density effects and dispersion in aquifer storage and recovery

Aquifer storage and recovery (ASR) involves the injection of freshwater into an aquifer for later recovery and use. This paper investigates three major factors leading to reduction in performance of ASR systems in brackish or saline aquifers: lateral flow, density-driven flow and dispersive mixing. Previous analyses of aquifer storage and recovery (ASR) have considered at most two of the above processes, but never all three together, and none have considered lateral flow and density effects together. In this analysis, four dimensionless parameters are defined to give an approximate characterisation of lateral flow, dispersive mixing, mixed convection (density effects during pumping) and free convection (density effects during storage). An extensive set of numerical models spanning a wide parameter range is then used to develop a predictive framework using the dimensionless numbers. If the sum of the four dimensionless numbers (denoted R_ASR) exceeds 10, the ASR operation is likely to fail with no recoverable freshwater, while if R_ASR < 0.1, the ASR operation is likely to provide at least some recovery of freshwater. The predictive framework is tested using limited data available from ASR field sites, broadly lending support to the framework. This study has several important implications. Firstly, the lack of completeness of field data sets in the literature must be rectified if we are to properly characterise mixed-convective flow processes in ASR operations. Once data are available, the dimensionless numbers can be used to identify suitable ASR sites and the desirable operational conditions that maximise recovery efficiencies.     

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Estimation of irrigation return flow on monthly time resolution using SWAT model under limited data availability

ABSTRACT

This study is an assessment of the return flow ratio of an irrigation abstraction using the flow records of a downstream stream-gauging station, with the example of the Kozluk scheme irrigated by diverted water from Garzan Creek flowing through the southeastern region of Turkey. In the planning reports of the major dam projects of the region, an unverified return ratio was assumed in eliminating the influence of this irrigation on the flow measurements of Garzan Creek. The correctness of this assumed return ratio is evaluated by analysing the monthly streamflow measurements of the Besiri station through a Soil and Water Assessment Tool (SWAT) model constructed with the coarse-scale topography, land use and soil data from open source databases. The results show the necessity of irrigation project-based return flow analyses using regional fine-scale datasets, instead of rule-of-thumb assumptions, to determine the effects of irrigation activities on flow regimes more accurately.     

Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. II: modelling the effects on groundwater flow

Anisotropy and heterogeneity of hydraulic conductivity (K) are seldom considered in models of mire hydrology. We investigated the effect of anisotropy and heterogeneity on groundwater flow in bog peat using a steady‐state groundwater model. In five model simulations, four sets of K data were used. The first set comprised measured K values from an anisotropic and heterogeneous bog peat. These data were aggregated to produce the following simplified data sets: an isotropic and heterogeneous distribution of K; an isotropic and homogeneous distribution; and an anisotropic and homogeneous distribution. We demonstrate that, where anisotropy and heterogeneity exist, groundwater flow in bog peat is complex. Fine‐scale variations in K have the potential to influence patterns and rates of groundwater flow. However, for our data at least, it is heterogeneity and not anisotropy that has the greater influence on producing complex patterns of groundwater flow. We also demonstrate that patterns and rates of groundwater flow are simplified and reduced when measured K values are aggregated to create a more uniform distribution of K. For example, when measured K values are aggregated to produce isotropy and homogeneity, the rate of modelled seepage is reduced by 28%. We also show that when measured K values are used, the presence of a drainage ditch can increase seepage through a modelled cross‐section. Our work has implications for the accurate interpretation of hydraulic head data obtained from peat soils, and also the understanding of the effect of drainage ditches on patterns and rates of groundwater flow. Copyright © 2002 John Wiley & Sons, Ltd.     

由GPS资料讨论中国大陆块体运动与应变模型的建立及检验

目前很多文献基于弹性力学假设,以小形变、均匀应变场为前提,利用GPS资料建立了中国大陆块体的整体旋转与均匀应变模型,并对模型进行了相关检验.但上述对模型的检验侧重从模型的拟合效果上进行分析和比较,忽略了对模型的适用性和假设条件合理性的检验.实际上,不同块体由于自身构造和外部条件的不同,其运动状况与应变状态并不完全一致,模型的适用性需进一步的验证.本文提出了两步检验方法对所建立的模型在中国大陆6个主要块体上进行了检验：第一步对应变参数进行显著性检验,第二步对应变均匀性进行检验,结果表明：根据所用数据精度,华南块体的刚性程度在95% 的置信水平上可以认为显著的,建立模型时不必考虑内部形变;其他块体内部形变显著,整体旋转与均匀应变模型并不能反映真实的地壳运动状态,必须考虑其形变的不均匀性.