Identification of aquifer and well parameters from step-drawdown tests

A new analysis technique has been proposed for interpreting transient step-drawdown test data. The proposed method is based on taking the derivative of the drawdown with respect to time for the entire pumping test period to eliminate the time-independent well-loss terms. The derivative function is subsequently integrated to obtain the time-dependent aquifer drawdown as a continuous function. The well-loss parameters are then obtained with higher accuracy once the aquifer behavior is identified. The proposed method is applicable for analyzing data obtained not only from ideal confined aquifers but also from other aquifer types (i.e. unconfined) and non-ideal aquifers (i.e. heterogeneous). The technique was tested for synthetically generated and field data; the proposed approach was noted to provide accurate aquifer and well-loss parameter estimates. The results of the proposed method were compared with those of some of the existing methods for analyzing step-drawdown test data and were found to be more reliable and robust.     

BP神经网络在BOULTON法确定潜水含水层参数中的应用

本文在分析具有迟后排水的N.S.Boulton第二潜水井流模型解析解的基础上,利用复合高斯求积法和学习速率、动量因子自适应的BP神经网络相结合对该模型数值求解,提出了确定潜水含水层系统参数的Boulton BP法。实例计算结果表明,用Boulton BP法求得参数正演出的降深—时间过程与抽水试验所得降深—时间过程拟合很好。该方法简单,快速,不需要将抽水试验所得降深—时间过程分为前、后两段分别求参,不仅简化潜水含水层的求参过程,而且有良好的求参效果。     

稳定井流直线图解法确定含水层参数

对于潜水井流,利用Dupuit公式计算的参数往往相差很大,难以直接选用。对Dupuit公式进行线性化,建立降深(s)与距离(r)的s-lnr直线关系,利用直线斜率、直线在x轴上的截距,可以直接求得唯一的渗透系数K和影响半径R。利用松花江河谷的承压水井抽水试验资料和洮儿河扇形地27个潜水抽水井及其观测孔的抽水试验资料,应用直线图解法分析计算,得到唯一的含水层参数T、K和R。将计算结果与直接利用Dupuit公式所计算的结果相比较,前者的求参效果较好。     

吉林省白城市地下水数值模拟及其变化趋势

为了给吉林省白城市地下水资源合理开发利用提供科学依据,运用地下水模拟技术对白城市地下水位变化趋势进行预测。首先在分析研究区水文地质条件的基础上建立水文地质概念模型,将含水系统概化为潜水含水层、越流层和承压含水层,并确定了系统的边界条件。运用地下水模型软件(GMS)对研究区地下水流进行了模拟和预测,预测结果表明在扩大开采的情况下,潜水含水层水位总体下降约1~2 m;随着开采量的增加和时间的增长,含水层在洮儿河扇形地、洮南、研究区东北部、镇赉、大安南部等地相继出现疏干现象,最大降幅达5~6 m。承压含水层在大安南部、镇赉北部等地水位下降较大,但是流场整体降幅不大,总体降幅在0.5 m以内。最后提出了地下水合理开发利用建议。     

Analytical solutions for constant-rate test in bounded confined aquifers with non-Darcian effect

This paper proposes a simplified analytical solution considering non-Darcian and wellbore storage effect to investigate the pumping flow in a confined aquifer with barrier and recharge boundaries. The mathematical modelling for the pumping-induced flow in aquifers with different boundaries is developed by employing image-well theory with the superposition principle, of which the non-Darcian effect is characterized by Izbash’s equation. The solutions are derived by Boltzmann and dimensionless transformations. Then, the non-Darcian effect and wellbore storage are especially investigated according to the proposed solution. The results show that the aquifer boundaries have non-negligible effects on pumping, and ignoring the wellbore storage can lead to an over-estimation of the drawdown in the first 10 minutes of pumping. The higher the degree of non-Darcian, the smaller the drawdown.     

侧向有限越流承压含水层中非完整井非稳定流模型及解析解

工程建设中当距离抽水井r=rb处水位基本没有变化或不受抽水影响时,或当此处存在止水帷幕时,含水层系统视为侧向有限延伸,rb为有限半径。为此,构建更加符合工程实际的侧向有限延伸的典型弱透水层-承压水层系统中非完整井非稳定流计算模型,同时考虑井径和井储效应的影响,应用Laplace变换和分离变量法得到了水位降深在拉氏空间下的解析解,并应用拉氏数值逆变换Stehfest法得到真实空间下的水位降深。新建立的解析解可以进一步退化为诸多已有解,并进一步将其与已知解和有限元数值解进行对比,验证了所得解的正确性和可靠性。基于新建解重点分析了侧向边界和井的完整性对承压水层水位降深的影响。结果表明:含水层系统的侧向有限边界仅对抽水后期的水位降深影响明显,含水层系统侧向无限延伸情况下的水位降深要大于情形1(在r=rb处为定水头边界)且明显小于情形2(在r=rb处为不透水边界)下的水位降深,rb越小,两者之间的误差越大;抽水井的完整性对整个抽水期间不同情形下的水位降深均有明显的影响,承压含水层顶板处的水位降深随着抽水井滤管的长度和埋深的增加而减小。     

Determination of the parameter pattern and values for a one-dimensional multi-zone unconfined aquifer

In an aquifer, heterogeneity plays an important role in governing groundwater flow. Hence, aquifer characterization should involve both the pattern and values of the hydrogeological parameters. A new analytical solution describing the one-dimensional groundwater flow in a multi-zone unconfined aquifer is presented, and a methodology developed from the analytical solution and a heuristic approach for determining the pattern and values of the aquifer parameters are proposed. The analytical solution demonstrates that the hydraulic head varies spatially and is influenced by aquifer heterogeneity. Simulated annealing, a heuristic approach, is incorporated with the solution to simultaneously identify the pattern and values of the hydraulic conductivity for a horizontal multi-zone unconfined aquifer. This approach may be used to give an approximate result for a two-dimensional problem by dividing the model area into a number of transects along the transverse direction, identifying the parameter values along the longitudinal direction for each transect, and then smoothing the identified results.     

确定含水层参数的全程曲线拟合法

含水层参数的确定是地下水资源评价、数值模拟、开发利用与保护以及科学管理的重要基础,抽水试验法是确定含水层参数的主要方法。根据含水层抽水试验的泰斯（Theis）公式,运用优化理论,结合计算机技术,提出利用抽水试验数据确定含水层参数的全程曲线拟合法。该方法充分考虑抽水过程中水位快速下降段、水位缓慢下降段和水位恢复段的观测资料,满足计算水位与实测水位之间的误差平方和最小、计算的水位过程线与抽水过程的实测曲线最佳拟合。实例研究表明,该方法速度快、精度高、求参结果最优,既适用于承压含水层求参,又可推广到潜水含水层求参。     

变异条件下内蒙古呼包平原地下水演化趋势

重点探讨气候变异对内蒙古呼包平原地下水资源、地下水流场和动态的影响。在介绍呼包平原水文地质条件的基础上,运用地下水模型软件PMWIN 5.1建立了研究区双层地下水流模型,并根据1998年地下水动态观测和统测资料对模型进行识别。通过对近51 a降水量的分析,确定将连续出现的2个特枯年作为地下水的变异条件,运用地下水模型对变异条件下地下水状态进行模拟,结果表明：2年间对地下水的严重超采,使地下水亏空量达75 296.34 万m³,潜水和承压含水层的地下水位大幅度下降,潜水位最大下降4 m左右,承压水位最大下降17 m左右;漏斗中心主要集中在开采量大的呼和浩特和包头市区。     

The Basin of Mexico aquifer system: regional groundwater level dynamics and database development

The aquifer system of the Basin of Mexico is the main source of water supply to the Mexico City Metropolitan Zone. Management of the Basin’s water resources requires improved understanding of regional groundwater flow patterns, for which large amounts of data are required. The current study analyses the regional dynamics of the potentiometric groundwater level using a new database called the Basin of Mexico Hydrogeological Database (BMHDB). To foster the development of a regional view of the aquifer system, data on climatological, borehole and runoff variables are part of the BMHDB. The structure and development of the BMHDB are briefly explained and then the database is used to analyze the consequences of groundwater extraction on the aquifer’s confinement conditions using lithology data. The regional analysis shows that the largest drawdown rates are located north of Mexico City, in Ecatepec (a region that has not yet received attention in hydrogeological studies), due to wells that were drilled as a temporary solution to Mexico City’s water-supply problem. It is evident that the aquifer has changed from a confined to an unconfined condition in some areas, a factor that is responsible for the large subsidence rates (40 cm/year) in some regions.     

基于过程模拟的李官堡水源地地下水污染预警

基于过程模拟的地下水污染预警是基于包气带和饱和带过程模拟的耦合实现的,且融合了涵盖整个地下水系统的四个预警指标。选取浑河冲洪积扇的李官堡水源地为例,基于Hydrus-1D和Visual Modflow分别进行包气带过程模拟和饱和带过程模拟,针对潜层和承压层的复杂程度分别制定了各自的预警临界值;经分析可知,随着时间的推移,地表污染物持续进入到地下水中,污染晕范围不断扩大,地下水中最大浓度值和水源井浓度值也持续增加,且承压含水层的预警级别要高于潜水含水层;同时基于潜层和承压层中污染物最大浓度值与时间分别近似呈线性关系和指数关系可进行长时间尺度的预警;并提出了零级预警区和一级预警区管理措施以监测和预防为主,二级预警区、三级预警区和四级预警区管理措施以控制和监测为主。     

双重介质含水层求参方法探讨

双重介质模型能较好地刻画天然裂隙含水层的非均质性及液体流动和溶质迁移特征。在石油开采、核废物处置安全评价和环境水文地质等方面受到广泛重视。但是,双重介质模型的实际应用远远落后于其理论发展。其原因除了双重介质模型复杂、计算工作量大以外,模型参数确定的理论和测试方法很不完善,难以取得两套较客观的模型参数是另一重要原因。本文就双重介质含水层参数的确定问题,从其介质特征和水流运动特征入手,探讨了双重介质含水层中抽水和注水条件下的井流动力学方程及其求解方法,并讨论了其适用条件及存在的问题,最后,用实际资料说明了求参过程,检验了公式。     

地下水位变化对透—阻型岩溶塌陷影响的分析

文章从一维地下水运动和渗透力学的角度,分析比较潜水位上升与承压水位下降对岩溶地区透—阻型盖层中阻水层渗透稳定性的影响,重点讨论了承压水位下降速度(降速)与下降深度（降深）对阻水层中渗透坡降的影响,结果表明:地下水位变化（潜水位上升或承压水位下降）产生的非稳定渗流不利于岩溶洞穴开口上方阻水层的稳定,承压水位的下降对岩溶开口附近处阻水层中渗透力的影响远大于潜水位的变化;在承压层水位最大降深确定的条件下,承压水位下降速度愈快,岩溶开口附近处阻水层中向下渗透力愈大。因此,在覆盖型岩溶地区抽取地下水时,为了减缓或避免覆盖型塌陷的发生,应同时控制好最大降深和最大开采速度。      

抽水试验和长观水位联合模拟确定含水层参数

为了更加准确地确定含水层参数,提出了利用地下水流数值模拟软件Modflow耦合大型抽水试验和长观孔水位两个模型的方法,对含水层参数进行了识别与校验。该方法求得的含水层参数避免了只利用抽水试验模型校验所带来的不确定性,同时利用抽水实验数据进行抽水与恢复的全过程曲线拟合,水位降深模拟相对误差小于6%,提高了含水层参数的准确性。     

Anisotropy and bedding effects on the hydro geological regime in a confined aquifer to design an appropriate dewatering system

Prediction of groundwater inflow into mining excavations is very important in order to design an effective dewatering system to keep the mine workings dry and create prolonged cone of depression. The effects of anisotropy ratio and bedding on the hydraulic head and drawdown curves of a dewatering test carried out in a fully penetrating well in a confined aquifer have been investigated. An existing numerical finite element model has been used to perform the simulations. The results of the numerical model are compared to those from analytical Jacob and Lohman solution for estimating hydraulic heads and drawdown curves. It was found that the anisotropy ratio and bedding should not have a significant effect on drawdown and the quantity of inflow into a confined aquifer. It was further found that taking the simultaneous effects of anisotropy and bedding into account reduces the differences in the results of analytical and numerical methods. Comparison of the field data and model predictions showed that, the modelling results for a three layer anisotropic aquifer fit well to the field data than those results obtained for a single layer aquifer and the relative error decreased from 4.81 % to 2.98 %.     

基坑工程井点降水分析计算

针对目前应用广泛的解析计算方法,本文分析了基坑工程井点降水的特点,按照有限元理论,在考虑不同边界条件、坑内外是否有水力联系等的基础上分析了降水中潜水、承压水地下水水位随时间的下降过程。因这种分析方法可考虑不同的基坑形状、计算边佃条件、井类型、含水层类型等,故具有较强的适用性。     

无限承压含水层中主孔涌水多孔观测定降深井流试验水文地质参数计算

利用自溢井(孔) 做放水试验确定含水层水文地质参数, 是简便而经济的确定含水层水文地质参数的试验方法.涌水试验基本上属于定降深井流试验.在无限承压含水层中进行涌水试验, 使涌水孔(井) 降深s_w保持不变, 涌水量随时间的延续而减少, 除涌水孔(井) 本身以外, 含水层中任意点的水头H随时间的延续而降低.在单孔(井) 定降深井流试验利用涌水孔1/Q-lg t直线图解法计算水文地质参数的基础上, 建立一孔(主孔) 涌水多孔观测水位降低的直线图解法, 即s/Q-lg t直线图解法.      

露头处具有淤泥层的滨海含水层系统中海潮引起的水头波动

滨海地区社会与经济的发展引发了各种各样的水文地质问题, 因此对滨海地区的水文地质条件, 尤其是地下水与海水之间的水力联系的研究尤为重要.本文考察了海底露头处具有淤泥层的滨海含水层系统中由海潮引起的水头波动.该系统由潜水含水层、承压含水层和介于其间的弱透水层构成.建立了该系统的数学模型, 得到了该模型的解析解.该解析解包含2个无量纲参数: 弱透水层的相对越流系数和淤泥层的相对透水系数.解析解表明, 淤泥层使各处的水头波幅被缩小了一个常数倍(波幅缩减因子), 并使各处的波动相位产生了一个不超过45°的正位移常数(时滞).该时滞对半日潮不超过1.5h, 对全日潮不超过3h.波幅缩减因子和正位移常数只与弱透水层的相对越流系数和淤泥层的相对透水系数有关.当这两个无量纲参数取某些特殊值时, 本次研究的解和前人考虑的几种较简单情形所对应的解一致.承压含水层中水头波幅随着淤泥层的相对透水系数的递增而严格递增, 随弱透水层的相对越流系数的递增而严格递减; 波动相位随着这2个无量纲参数的递增均严格递减.分析表明弱透水层的越流和淤泥层的存在均对承压含水层水头波动有显著影响.      

The Investigation of Mass Transfer in the Karasu Karstic Aquifer, Konya, Turkey

In this study, the changes in the chemical composition of the groundwater along a flow path were examined by using the water samples collected from unconfined, semi-confined and confined parts of the Karasu karstic aquifer. It was determined that transport of bicarbonate, calcium, and magnesium was dominant in unconfined and semi-confined parts of the aquifer, whereas calcite and dolomite precipitate in the confined parts. On the other hand, gypsum dissolution is present in all parts of the aquifer. In addition, the computed saturation indices explain the occurrences and precipitation of travertines in the Goksu Valley, which is the discharge area for the aquifer. Electronic Publication     

Assessing hydrochemical evolution of groundwater in limestone terrain via principal component analysis

This paper describes the use of multivariate statistical analysis to trace hydrochemical evolution in a limestone terrain at Zagros region, Iran. The study area includes a deep confined aquifer, overlaid by an unconfined aquifer. The method involves the use of principal component analysis (PCA) to assess and evaluate the hydrochemical evolution based on chemical and isotope variables of 12 piezometers drilled in both the unconfined and confined aquifers. First PCA on all variables shows that water–rock interaction under different conditions with respect to the atmospheric CO₂ is the main process responsible for chemical constituents. As a result, combinations of several ratios such as Ca/TDS, SO₄/TDS and Mg/TDS with physico-chemical and isotope variables reveal different hydrochemical evolution trend in the aquifers. Second PCA on the selective samples and variables reveals that displacement of the unconfined samples from dry to wet season follows a refreshing trend towards river samples that is characterized by reducing electrical conductivity and increasing sulphate and tritium contents. However, the refreshing trend cannot be traced in the confined aquifer samples suggesting no recharge from river to the confined aquifer. Third PCA reveals that, chemical composition of water samples in the unconfined aquifer tends to have considerable difference from each other in the end of recharge period. In contrast, the confined aquifer samples have a tendency to show similar chemical composition during recharge period in comparison to end of dry period. This difference is caused by different mechanism of recharge in the unconfined aquifer (through the whole aquifer surface) and the confined aquifer (through the limited recharge area).