基于数值模拟的含水层参数识别

利用地下水流数值模拟软件Modflow耦合大型抽水试验和长观孔水位两个模型的方法,对含水层参数进行了识别与校验,该方法求得的含水层参数避免了只利用抽水试验模型校验所带来的不确定性。同时利用抽水实验数据进行抽水与恢复的全过程曲线拟合,水位降深模拟相对误差小于6％,并且对均衡水量、泉水的溢出量以及末流场拟合程度进行了分析,均衡水量模拟相对误差最大为2．6％,泉水溢出量模拟相对误差最大为8．14％,末流场拟合程度也较好,从而提高了含水层参数的准确性。实例研究表明,这种求参方法精度高、切实可行,可以推广到其他的含水层参数确定中。     

Review: “Jacob’s Zoo”—how using Jacob’s method for aquifer testing leads to more intuitive understanding of aquifer characteristics

The interpretation of aquifer responses to pumping tests is an important tool for assessing aquifer geometry and properties, which are critical in the assessment of water resources or in environmental remediation. However, the responses of aquifers, measured by time-drawdown relationships in monitoring wells, are nonunique solutions that are affected by many factors. Jacob’s Zoo is a collection of graphical interpretations that allows students and practitioners to develop an intuitive feel for how natural hydrogeological systems work, and develop a set of skills that provide a better understanding of aquifer properties far beyond interpretation of pumping tests. Jacob’s Zoo, based on the work of Jacob (1950), fosters a deeper understanding, although few practitioners realize the full utility of the method. Jacob CE (1950) Flow of groundwater, In: Rouse H (ed) Engineering Hydraulics, Wiley, New York. P 321–386.

     

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

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

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.     

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

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

无观测孔双井干扰抽水求取承压含水层导水系数方法探讨

煤矿疏降水工程需要以含水层组参数为基础进行合理布设。结合某矿区内的非稳定流双孔干扰抽水试验,在无观测孔情况下,根据抽水井资料分别利用通用直线法、水位恢复法、解析法及优化拟合四种方法确定承压含水层的导水系数,并利用后期疏降水孔资料对上述方法进行验证,结果表明通用直线法和水位恢复法计算结果可靠。     

A case history of field pumping tests in a deep gravel formation in the Taipei Basin,Taiwan

33 large-diameter wells embedded in 2-m thick, 63-m deep diaphragm walls were constructed to reduce both the uplift pressures and the groundwater inflow during the excavations. As the actual thickness of the pumped aquifer is unknown, the installed wells are regarded as partial penetration wells. Single-well and multi-well pumping tests were conducted in the deep gravel formation of Taipei Basin to derive the hydraulic parameters and to investigate the drawdown characteristics at both the construction and remote sites. However, the tidal effect on the drawdown of both the pumping well and nearby observation wells was found significant. Additionally, wellbore storage, skin, and leakage need to be taken into account for deriving the hydraulic parameters. Hence, a method to remove these five factors influencing the drawdown curve is developed, which takes advantage from the late-time characteristics of drawdown data and the early-time behavior of drawdown. Some currently available semi-log graphic techniques are therefore proven applicable for parameter determination. Validity of the proposed method is verified by the good agreement between the calculated and the measured drawdown of both the pumping well and observation well.     

Automatic parameter estimation applied on a groundwater model: The problem of structure identification

The properties of a Pleistocene aquifer located in NW Germany were evaluated by means of an inverse two-dimensional numerical model in the vertical plane. The steady-state drawdown recorded at the end of a pumping test was used to calibrate the model with a Gauss-Newton parameter estimation method. The results were analyzed within a formal statistical procedure by means of eigenvalues and eigenvectors. Although a mathematically meaningful value for the standard error can be obtained, the resulting set of permeabilities depends strongly upon the choice of the zonation pattern.     

Multi-observation well aquifer test case study: is recovery coincident with the cessation of pumping?

Due to increases in water demand, the City of Kenedy, TX, USA must expand their small drinking water supply in the Gulf Coast aquifer system. Groundwater wells owned by the City of Kenedy, Karnes County, TX were examined to estimate properties for the Jasper aquifer. Conditions of four wells were assessed, after which two wells were rehabilitated and used as pumping wells in aquifer tests. Aquifer tests show that recovery in observation wells was not coincident with the cessation of pumping. Post-pumping data were selectively excluded so that only recovery data were used for analyses. Transmissivity for the Jasper aquifer ranges from 102 to 242 m² d^?1, and storativity ranges from 6.9E?05 to 3.3E?04. Transmissivity computed from recovery data was approximately 25 % higher than transmissivity computed from time-drawdown data. Field measured specific capacities and drawdowns were compared to theoretical specific capacities and drawdowns to calculate pumping well efficiencies in the range of 52.2–99.4 %. This study indicates that water demand for the City of Kenedy could be met by incorporating the tested wells into the water supply system. Future studies should be designed to estimate groundwater recharge rates and a complete water balance for computing a sustainable maximum annual yield.     

Re-evaluating the US Geological Survey’s pumping tests (1967) in the Punjab region of Pakistan for use in groundwater studies

In 1967, the US Geological Survey (USGS) published the results of 141 pumping tests carried out throughout the Pakistani Punjab to establish representative hydraulic parameters of its large aquifer. Many authors have since concluded that the USGS had over-estimated the horizontal hydraulic conductivity (k _r) by 25–100 %, leaving vertical anisotropy and aquifer depth unresolved. No test wells have ever been drilled below 450 m to reach the base of the aquifer, although petroleum explorations mention depths between 1,500 and 4,500 m. After comparison and re-evaluation of all related papers, this study concludes that the USGS interpretation was correct, that its hydraulic values still stand without change, and that the USGS’s applied distance drawdown interpretation is valid to prevent influence of partial penetration on the results. This study also uniquely resolved vertical anisotropy and aquifer thickness by using early- and late-time drawdowns separately and proper scaling of the coordinates, which has often been omitted. With appropriate scaling, all interpretations match the data. The representative hydraulic aquifer values are: k _r?=?65 m/d, vertical anisotropy k _r/k _z?=?25 and aquifer depth 500–1,500 m. The conclusion is that these values can be used, at least as first estimates, for groundwater studies in the Pakistani Punjab.     

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.     

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.     

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.     

Determination of the rate and volume of leakage using the slopes of time-drawdown data

Quantification of leakage is very important in the selection and design of the remediation systems of leaky aquifers that receive contaminated leakage. This is an approach for the calculation of leakage using only two slopes of time-drawdown data. These slopes represent before and after the start of leakage, and are applied to four examples. Results generally agree with those determined by the Hantush approach. Comparison of the two approaches, however, shows that the Hantush approach quantifies leakage using three aquifer parameters (transmissivity, storativity, and leakage factor), the value of which depend on the pumping test method used; it assumes constant hydraulic head in the aquifer supplying leakage, which may not be valid under field conditions; and it ignores differences between the viscosities of the leakage water and the aquifer water, which influence the leakage rate. The proposed approach is free from all three limitations.     

可地浸砂岩型铀矿床抽水试验方法对比及有关问题探讨

本文对新疆伊犁盆地库捷尔太铀矿床抽水试验及参数计算方法做了较为详细的对比分析，总结出适合于可地浸砂岩型铀矿床不同工作阶段求取 含水层渗透系数的抽水试验方法；阐述了在自喷井中用扬水试验求取含水层渗透系数方法的可行性及其在实际应用时应注意的问题；提出利用数量固定不可调的深井泵进行变流量抽水以控制不同的降深，解决了多年来因空气压缩机性能所限难以进行长时间抽水的问题，并对抽水试验的延续时间、降深次数以及降深等参数提出了具体建议。     

地下水加坡理论在井筒涌水量预算中的应用

地下水加坡理论可有效的消除抽水试验引起的无效降深,结合云驾岭煤矿井筒检查孔抽水试验资料,对该井简检查孔涌水量进行了预算,预算副井筒啮。砂岩含水层涌水量为523．03m^3／h,实际涌水量为455m^3／h．误差率14．95％,可见运用加坡理论计算的矿井涌水量与井筒实际涌水量基本相符。     

Comparative pump test analysis of a multiple bore array in a thin unconfined aquifer

Standard aquifer tests have been applied to the results of a 13‐day pump test on a circular array of 12 bores in a thin unconfined aquifer in northern Tasmania. During pumping the array behaves in a manner similar to a large diameter bore. By assuming that it is also hydrologically equivalent to a central imaginary bore, all the techniques produce ‘acceptable’ values for transmissivity and specific yield. Values of T in the range 56–99 m²/day and S in the range 0.16–0.29 are indicated. Independent assessment of S using a method developed from first principles yields a value as high as 0.33. No further refinement is possible, although drawdowns are accurate to the nearest millimetre and were measured in 13 observation bores. Each individual method of analysis produces apparently precise results, but the inherent imprecision of all standard techniques is obvious when a comparative approach is adopted. Although T and S are probably reasonably constant throughout the aquifer, ‘correct’ values were only approached after long pumping times. The duration of the pump test remains a more critical factor in calculating T and S than the method of analysis. A simple empirical prediction theory has been developed (at this stage for the Tas‐manian aquifer only) which describes the drawdown configuration around any number of bores for any spatial arrangement, and which is independent of the imprecisely known T and S values obtained from formal analysis.     

Global curve-fitting for determining the hydrogeological parameters of leaky confined aquifers by transient flow pumping test

Proper management of groundwater resources requires an accurate evaluation of the parameters (hydraulic properties) that control the movement and storage of groundwater. Hydrogeological parameters are the basis of groundwater evaluation, modeling, and management and so on. A global curve-fitting method incorporating pumping test data and water table recovery data was introduced in the present study. The principal and procedures of the method were elucidated in detail. The drawdown and recovery data from two sets of transient flow pumping test conducted in no. 2 water source site of Shizuishan city were used to verify the calculation accuracy of the proposed method. The hydrogeological parameters were also estimated with traditional type curve-fitting method on the basis of formula derived by Hantush and Jacob. The hydrogeological parameters calculated by the two methods were compared and the results show that the parameters obtained by the global curve-fitting method are a little bigger than but very close to those obtained by the traditional type curve-fitting method. The proposed method which possesses three major advantages is feasible and reliable in aquifer parameter identification. A comparative study on various methods for parameter identification is required and expected in future study.     

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.