井周流场,水跃值与井内水位降深的关系

本文以单井稳定流抽水为例，从井固流场与地下水向集水井幅射流的涌水量方程入手，阐明井结构，井周泥浆堵塞和井内水位降深的关系。证明了水跃值是井周泥浆堵塞降深与井周三维流域或紊流引起的二次项降深（简称“井损“）之代数和，并得出了水跃值不同组成部分随井内降深变化而变化的规律。从而提出了消除水跃值，利用裘布依涌水量方程进行水文地质参数计算。     

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.     

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.     

Forward Coupled Modeling and Assessment of Gravity Anomalies Caused by Pumping Tests in Unconfined Aquifers Under Unsteady-State Conditions

Analysis of drawdown data from pumping tests is the most common method for identification of aquifer properties and estimation of hydraulic parameters in hydrogeology. Limited hard hydraulic data from wells, expensive and intrusive, can be supplemented with geophysical spatially distributed information. Gravity, sensible to mass variations, can provide information about hydraulic parameters and storage coefficients of aquifers. A satisfactory use of gravity requires the computation of temporal evolution of gravity anomalies related with the dynamic hydraulic process. For simplified and ideal problems, the calculation can be accomplished with analytical solutions. However, the limitations of analytical solutions in capturing some relevant characteristics, like heterogeneities, that require the use of numerical models, are well known. This article shows how the time-dependent gravity anomaly caused by a pumping test can be calculated in a unique code, in a coupled manner and in a single run. The methodology avoids the use and transference of data between various codes and can be utilized in a coupled hydrogeophysical inversion procedure. Additionally, it is flexible enough to include a wide range of characteristics that can be encountered in any pumping test, and is an adequate framework for field survey planning via analysis and study of time-dependent gravity curves.     

垃圾填埋场抽水试验及降水方案设计

垃圾填埋场中的渗滤液水位过高会引发一系列环境和稳定问题,工程上可用竖井抽水降低渗滤液水位。通过在填埋场现场进行抽水试验,确定垃圾土的渗透系数和抽水影响半径,在此基础上对填埋场降水的瞬态流问题进行有限元模拟,分析了抽水井口径和间距对填埋场降水的影响,提出了降水方案的设计步骤和方法。抽水试验表明,现场垃圾的渗透系数约为3.6×10-4cm/s,抽水影响半径约为20m。数值分析表明,井径的变化对于降水效果影响不大,而合理选择抽水井间距对降水十分关键。进行抽水方案设计时,应首先根据工期和降水幅度要求计算井间距,按井的出水速度选择水泵,再根据水泵确定井径,最后根据井径和过滤层形式确定钻孔尺寸并选择钻机。     

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.     

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

含水层参数对于定性分析区域地下水资源评价、数值模拟及预报、开发利用与保护及科学管理具有重要意义,而江汉平原水文地质参数求解的研究却相对较少.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.通过结果对比分析发现对于单孔非稳定抽水试验,对前期水位降深数据筛选分析,同样可以利用非稳定井流理论反演含水层参数,结果比较可靠.      

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

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

Parameters affecting the occurrence of quicksand and the drying up of large diameter wells that gain water from the bottom: a case study from Iran

Large diameter fully cased wells that gain water from the bottom are often dug in sandy and collapsible aquifers. They have cylindrical vertical walls lined with brick or concrete. The well bottom is partially filled with aquifer material through which the flow is vertically upward. When the vertical hydraulic gradient reaches a critical value, quicksand occurs and the well structure can be destroyed. Another difficulty encountered is drawdown in the wellbore and the drying up of the well. To overcome these problems, the flow around and beneath these wells is numerically simulated. The simulation results are used to investigate the effect of well and aquifer parameters on quicksand and drawdown. For practical purposes, the dimensionless drawdown-time and dimensionless vertical gradient-time curves are developed. It was found that the ratio of filling material thickness to well radius affects the shape of these type curves. The type curves may be used to predict the time after pumping commences when quicksand occurs and the well dries up. They are also useful to design the safe pumping rate and duration as well as the optimum well radius. These are demonstrated by analyzing the pumping test data from a case study in the arid Chah Kutah region, southern Iran.     

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.     

Analysis of pumping test data for determining unconfined-aquifer parameters: Composite analysis or not?

Recently, composite analysis (CA), which simultaneously analyzes all drawdown data from multiple observation wells, has been applied to determine the hydraulic parameters of an unconfined aquifer. Moench (1994) claimed that the value of specific yield (S _y) determined from non-composite analysis (nonCA) is sometimes unrealistically low as compared with that obtained by water-balance calculation, and results from CA are better representative of aquifer properties than those from nonCA. To examine the validity of this assertion, the drawdown data from a pumping test conducted at Cape Cod, Massachusetts, USA, were analyzed using both nonCA and CA methods. The results show that the mean estimates of hydraulic conductivity and S _y determined from CA are close to those determined from nonCA. In some cases the analysis based on CA also results in low estimates of S _y as compared with those determined based on nonCA. A hypothetical case study is presented, which examines the effect of measurement errors on the estimated parameters. The results indicate that the CA method also gives poorer estimates of S _y than the nonCA method if the pumping test data contain measurement errors. Moench AF (1994) Specific yield as determined by type-curve analysis of aquifer-test data. Ground Water, 32(6):949–957.     

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

含水层的水文地质参数是进行地下水资源计算、地下水污染防控等所必需的基础数据,结合数值模拟技术进行含水层参数反演很有必要.按照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.对于深层承压含水层而言,越流补给量较大.采用数值模拟方法结合抽水试验数据求参,综合考虑了含水层结构,拟合效果好,所得结果更加可靠.      

小范围群井与单井抽水试验推求水文地质参数的比较分析

小范围岩溶区水文地质参数的推求是一项困难的工作.过去常常定性选取抽水孔与观测孔,利用单井抽水试验资料进行推求,这样往往造成推求误差过大.该文通过建立多层次多目标系统模糊决策模型,分别定量的优选了干扰井群抽水试验的最佳观测孔和单井抽水试验的最佳抽水试验单井;并运用了井群概化的方法,采用干扰井群抽水试验资料推求小范围水文地质参数,并且与单井抽水 试验资料推求的参数进行了比较分析,经分析应用方法取得了较好的效果.     

Hydraulic tomography analysis of municipal-well operation data with geology-based groundwater models

The sustainable management of groundwater resources is essential to municipalities worldwide due to increasing water demand. Planning for the optimized use of groundwater resources requires reliable estimation of hydraulic parameters such as hydraulic conductivity (K) and specific storage (S_s). However, estimation of hydraulic parameters can be difficult with dedicated pumping tests while municipal wells are in operation. In this study, the K and S_s of a highly heterogeneous, multi-aquifer/aquitard system are estimated through the inverse modeling of water-level data from observation wells collected during municipal well operations. In particular, four different geological models are calibrated by coupling HydroGeoSphere (HGS) with the parameter estimation code PEST. The joint analysis of water-level records resulting from fluctuating pumping and injection operations amounts to a hydraulic tomography (HT) analysis. The four geological models are well calibrated and yield reliable estimates that are consistent with previously studies. Overall, this research reveals that: (1) the HT analysis of municipal well records is feasible and yields reliable K and S_s estimates for individual geological units where drawdown records are available; (2) these estimates are obtained at the scale of intended use, unlike small-scale estimates typically obtained through other characterization methods; (3) the HT analysis can be conducted using existing data, which leads to substantial cost savings; and (4) data collected during municipal well operations can be used in the development of new groundwater models or in the calibration of existing groundwater models, thus they are valuable and should be archived.

     

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

为了定量研究与评价含水层之间地下水水力联系程度,首次提出了水力联系系数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）,水力联系极弱至弱。     

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.

---

     

Determining the optimal pumping duration of transient pumping tests for estimating hydraulic properties of leaky aquifers using global curve-fitting method: a simulation approach

Global curve-fitting method (GCFM) is regarded as an effective approach in hydrogeological parameter estimation, as it integrates and uses pumping data and water recovery data of a transient pumping test for parameter estimation. The impacts of pumping duration on hydrogeological parameter estimation by GCFM were investigated in the present study using 2 in situ pumping tests and 24 simulated transient pumping tests. Empirical formulas for determining the optimal pumping duration were derived. The study results suggest that pumping duration will have impacts on the accuracy of hydrogeological parameter estimation. When pumping duration is longer than a certain period, relative errors of hydrogeological parameter estimation keep relatively stable within an acceptable limit. Therefore, it is unnecessary to continue the pumping for a very long time after the groundwater level has become stable. When the change rate of drawdown over time (γ) in an observation well located within a distance of 100 m to the pumping well reaches 0.134, the pumping can be stopped. If there are more than one observation wells in a pumping test, the smallest γ value should be selected to determine the optimal pumping duration. This research is meaningful in the instruction of pumping tests, and will reduce test costs greatly.     

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

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井流试验方法不仅能确定水文地质参数,并且能丰富对含水层特性的认知。     

Flowing fluid electrical conductivity logging of a deep borehole during and following drilling: estimation of transmissivity,water salinity and hydraulic head of conductive zones

Flowing fluid electrical conductivity (FFEC) logging is a hydrogeologic testing method that is usually conducted in an existing borehole. However, for the 2,500-m deep COSC-1 borehole, drilled at Åre, central Sweden, it was done within the drilling period during a scheduled 1-day break, thus having a negligible impact on the drilling schedule, yet providing important information on depths of hydraulically conductive zones and their transmissivities and salinities. This paper presents a reanalysis of this set of data together with a new FFEC logging data set obtained soon after drilling was completed, also over a period of 1 day, but with a different pumping rate and water-level drawdown. Their joint analysis not only results in better estimates of transmissivity and salinity in the conducting fractures intercepted by the borehole, but also yields the hydraulic head values of these fractures, an important piece of information for the understanding of hydraulic structure of the subsurface. Two additional FFEC logging tests were done about 1 year later, and are used to confirm and refine this analysis. Results show that from 250 to 2,000 m depths, there are seven distinct hydraulically conductive zones with different hydraulic heads and low transmissivity values. For the final test, conducted with a much smaller water-level drawdown, inflow ceased from some of the conductive zones, confirming that their hydraulic heads are below the hydraulic head measured in the wellbore under non-pumped conditions. The challenges accompanying 1-day FFEC logging are summarized, along with lessons learned in addressing them.     

基于遍历搜索算法的水文地质参数优化求解

为了解决传统的Theis配线法在求参过程中无法利用全部抽水试验数据、手动配线效率较低、求参过程受人为主观因素影响较大等问题,利用抽水试验数据,以Theis公式和叠加原理为理论基础,以理论计算降深与实测降深的Nach-Sutcliffe效率系数值达到最大为目标函数,基于Matlab软件编程,再利用遍历搜索算法对导水系数（T）、贮水系数（S）在其对应范围内进行遍历,实现了水文地质参数的自动精确求解。将遍历搜索算法应用于2个抽水试验实例的水文地质参数求解,并与传统求参方法的计算结果进行比较分析。结果表明：利用遍历搜索算法求参的计算结果与传统求参方法的计算结果相近,表明利用遍历搜索算法求解水文地质参数有效可靠;2个抽水实例利用遍历搜索算法得到的参数对应的理论计算降深与实测降深的Nach-Sutcliffe效率系数值分别为0.996 5、0.970 8,且相比传统的求参方法而言更接近1.000 0,说明其拟合程度更好。