水力层析法与克立金法估算非均质含水层渗透系数场比较

含水层非均质性空间分布特征的识别,是准确模拟地下水流和污染物运移的前提.基于室内非均质含水层砂箱实验,分别利用水力层析法和克立金插值法刻画了非均质含水层渗透系数场.研究结果表明:(1) 水力层析法与克立金法相比,不仅可以更好地刻画非均质含水层渗透系数场,还可以更高精度地预测地下水流过程;(2) 水力层析抽水实验中,通过增加抽水实验组数可以有效地提高水力层析参数反演的精度,但是抽水实验组数增加到一定程度以后,再增加抽水组数不会显著提升参数反演的效果.后续需要进一步研究水力层析抽水实验合适的组数,进一步对抽水井进行优化布设.      

含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响

为探讨含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响,基于非均质含水层砂箱实验,分别用传统等效均质模型、克立金插值和水力层析刻画含水层渗透系数场,并探讨了先验信息对水力层析结果的影响.将不同方法估算的渗透系数场用以预测地下水流和溶质运移过程,以此判断不同方法估算结果的优劣,分析含水层非均质性对地下水流和溶质运移的影响.结果表明:与克立金插值法相比,水力层析法可以更好地刻画含水层非均质性,较准确地预测地下水流和溶质运移过程;钻孔岩心渗透系数样本值作为先验信息可以提高水力层析法估算结果的精度;传统等效均质模型无法准确预测地下水流和溶质运移过程.含水层非均质性的增强将导致溶质污染羽分布形态和运移路径的空间变异性增强,并且优势通道直接决定溶质的分布及运移路径.      

基于水力走时反演刻画裂隙含水层非均质性

获取含水层水力参数空间非均质分布信息是研究地下水渗流、地下水污染物运移等诸多地下水问题的重要基础。然而,受常规勘察技术所限,含水层的非均质性难以直接刻画,尤其在裂隙介质含水层中,水力参数分布的非均质性更加突出,进一步增加了刻画的难度。针对该问题,本研究首先通过德国哥廷根大学北校区内的Neutra试验场地进行的64次抽水试验,获取了该场地不同深度观测点的水头响应曲线,然后使用解析方法对64组数据进行分析和参数评估,同时采用水力走时的方法对井间水力参数分布进行反演计算,最后将得到的结果分别与经典解析解参数估计和热示踪试验结果进行对比验证。结果表明,解析解参数估计的结果虽然能够在一定程度上展现含水层的垂向非均质性,但是无法刻画井间含水层水力参数的非均质分布;与热示踪试验结果的对比验证了基于水力走时反演的水力层析法在刻画裂隙介质水力参数的空间非均质分布的可靠性。     

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

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

运用水力层析法对含水层裂隙带成像

含水层裂隙带的识别,是准确模拟计算地下水流和污染扩散的前提。本研究通过数值实验,运用近年来发展起来的水力地质层析方法,对含水层裂隙带成像。文章设定一个17m×17m的二维垂直饱和裂隙含水层,含水层有2口斜井。在斜井与裂隙带的交叉点分别进行抽水试验,由感应器收集水位响应数据,运用VSAFT2对裂隙带进行识别。数值实验证明,水力层析法结合地下水流的逆算方法——连续线性算法(SSLE),可以比较准确地识别非均质含水层的参数分布;非稳定流条件下的识别结果好于稳定流的结果;信息齐全条件下的结果好于资料缺乏时的结果。     

水力层析法刻画非均质含水层K与S采样时间优化设计

水力层析法(hydraulic tomography,简称HT)是一种高精度刻画含水层非均质性的抽水试验新技术,通过在不同地点及深度进行多组抽水试验并收集含水层水位变化资料,进行地下水流数值反演,对水文地质参数的空间分布做详细的识别。由于观测水位数据随着时间变化与渗透系数K、储水系数S的相关性不断变化,因此运用不同时刻水位数据反演参数K与S精度不一样。基于水力层析理论,采用剖面二维承压含水层数值模型研究了时间采样对估计参数K与S精度的影响。研究结果表明,随着采样时间的增大,渗透系数K的估计逐渐变好,最后趋于一个稳定值;在K已知条件下,采样时间t0(t0定义为JACOB直线图解法中s-lgt曲线直线段的延长部分与lgt轴的交点所对应的时间)时刻之前,储水系数S的估计随着采样时间的增加而变好,t0之后随着采样时间的增加反演效果变差。     

含水层层状非均质对地下水流系统的影响

区域尺度上含水层非均质具有复杂的结构性和随机性,难以准确刻画,造成非均质对区域地下水流系统的影响机制研究不够深入。本文以鄂尔多斯盆地白垩系地下水流系统为研究实例,选择典型剖面,采用剖面二维随机数值模拟方法,通过对比不同非均质刻画方法下地下水流场的变化,探讨含水层层状非均质对地下水流系统的影响机制。结果显示,均质条件下模型各向异性（含水层水平和垂向渗透系数比值Kh/Kv）取值为1000时,地下水流场与实际条件较为接近;非均质条件下,渗透系数方差取值0.91,水平相关长度取值5000 m,Kh/Kv取值150时,接近实际条件。研究表明,在大尺度地下水流模拟研究中,采用水平相关长度、渗透系数方差和各向异性值三个变量生成的随机场能很好地刻画含水层的层状非均质特征及其对水流系统的影响控制作用。由于含水层不同尺度层状非均质的叠加效应,采用均质各向异性介质等效概化含水层层状非均质性会造成等效各向异性值偏大失真的效应。     

基于EnKF融合地球物理数据刻画含水层非均质性

含水层非均质性的刻画是模拟地下水中污染物运移的关键。以渗透系数为研究对象,构建了综合集合卡尔曼滤波方法、有效电阻率模型与地下水运移模型的同化框架,通过融合地球物理观测数据与污染物浓度观测数据来推估渗透系数的空间分布。基于理想算例,验证了该同化框架刻画含水层非均质渗透系数场的有效性,并针对不同初始参数信息与观测类型对比了耦合与非耦合水文地球物理方法的适用性。研究结果表明:基于集合卡尔曼滤波方法同化多种类型的观测数据,可有效地推估非均质参数空间分布。当初始信息较准确时,耦合方法的参数推估精度更高;初始信息存在偏差时,非耦合方法有更好的同化效果。由于非耦合方法计算成本较低且对初始信息缺失时适用性更强,在实际应用中可先基于非耦合方法初步估计参数,再利用耦合方法进一步提高参数推估精度。融合多种类型观测数据可有效提高参数推估效果。     

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

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

抽水含水层对非抽水含水层水位影响的效应分析

利用自动水位计记录的抽水试验过程中不同层位含水组之间的动态水位,分析了在抽水试验最初阶段时抽水含水组与非抽水含水组之间的水力联系。研究结果表明:(1)抽水含水层对非抽水含水层存在水位影响效应,在抽水初期,水位影响效应较大;当抽水含水层水位达到稳定或缓慢下降时,水位影响效应随之消失。(2)单位降深水位影响效应显示,上部含水层距离抽水含水层越远,水位影响效应越大;下部含水层距离抽水含水层越远,水位影响效应越小。(3)利用三维地下水流数值模拟方法和应力应变对含水层水位动态影响的原理对水位影响效应机理进行分析,发现水位影响效应与含水层的固体骨架压缩系数(α)、水体积压缩系数(β)、孔隙度(n)和边界控制系数(Cm)等参数有关。     

Automatic estimation of aquifer parameters using long-term water supply pumping and injection records

Analyses are presented of long-term hydrographs perturbed by variable pumping/injection events in a confined aquifer at a municipal water-supply well field in the Region of Waterloo, Ontario (Canada). Such records are typically not considered for aquifer test analysis. Here, the water-level variations are fingerprinted to pumping/injection rate changes using the Theis model implemented in the WELLS code coupled with PEST. Analyses of these records yield a set of transmissivity (T) and storativity (S) estimates between each monitoring and production borehole. These individual estimates are found to poorly predict water-level variations at nearby monitoring boreholes not used in the calibration effort. On the other hand, the geometric means of the individual T and S estimates are similar to those obtained from previous pumping tests conducted at the same site and adequately predict water-level variations in other boreholes. The analyses reveal that long-term municipal water-level records are amenable to analyses using a simple analytical solution to estimate aquifer parameters. However, uniform parameters estimated with analytical solutions should be considered as first rough estimates. More accurate hydraulic parameters should be obtained by calibrating a three-dimensional numerical model that rigorously captures the complexities of the site with these data.     

A comparison of estimated and calculated effective porosity

 Effective porosity in solute-transport analyses is usually estimated rather than calculated from tracer tests in the field or laboratory. Calculated values of effective porosity in the laboratory on three different textured samples were compared to estimates derived from particle-size distributions and soil–water characteristic curves. The agreement was poor and it seems that no clear relationships exist between effective porosity calculated from laboratory tracer tests and effective porosity estimated from particle-size distributions and soil–water characteristic curves. A field tracer test in a sand-and-gravel aquifer produced a calculated effective porosity of approximately 0.17. By comparison, estimates of effective porosity from textural data, moisture retention, and published values were approximately 50–90% greater than the field calibrated value. Thus, estimation of effective porosity for chemical transport is highly dependent on the chosen transport model and is best obtained by laboratory or field tracer tests. Received, March 1997 · Revised, August 1997 · Accepted, August 1997     

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.     

Determination of porosity and storage capacity of a calcareous aquifer (France) by correlation and spectral analyses of time series

Time-series analyses were used to investigate the relationships between barometric pressure changes, earth tides, and water-level fluctuations in a confined aquifer. The method was applied to data from the fractured aquifer at the Hydrogeological Experimental Site in Poitiers (France) and used to yield estimates of the aquifer’s storage capacity, porosity and barometric efficiency. The aim is to address relevance of these analyses for an aquifer showing both fracture draining and confined karstic flow in thin strata. Cross-correlation and spectral analysis are used to compare water-level head and atmospheric-pressure fluctuations. Porosity and storage capacity are calculated using this method and compared to results from petrophysical measurements and hydraulic pumping tests, respectively. The storage coefficients calculated by the time-series analyses are in agreement with those obtained by interpretation of the interference pumping tests. Conversely, porosity values calculated by time-series analyses are underestimated compared to those obtained by other methods.     

Water table fluctuation analyses and associated empirical approach to predict spatial distribution of water table at Yeşilköy/AgiosAndronikos aquifer

The groundwater of the deep Ye?ilköy aquifer is the only water resource for agricultural and domestic consumption at the Karpaz Peninsula of Cyprus, which stretches approximately 100 km from the northeast of capital Nicosia to the northern tip of Cyprus. During the last decade, over-pumping and following dry periods have depleted the groundwater resources and the water surface elevation of the aquifer has dropped. The aim of this study is to understand the behavior of the Ye?ilköy aquifer in the last decade for the proper management of groundwater resources. This has been achieved based on well survey and field survey studies, monitoring programs followed by pumping tests, and safe yield analysis. Most of the research effort has been focused on field and well survey studies to quantify agricultural water consumption and abstraction rates from the aquifer. A long-term groundwater level monitoring program, short-term continuous groundwater monitoring and pumping tests provided information for the regression analyses while deriving a sixth order polynomial relationship between the period parameter and the head parameter. The equation was helpful to predict the short-term behavior of the water level when the present hydrogeological conditions prevail. The pumping test results satisfied the hydraulic properties of calcarenite formation yielding T = 1,782 m²/day and S = 0.0012. The results of safe yield analysis show that the annual deficit of the aquifer is 0.496 million cubic meters (MCM), which is equivalent to a 0.6 m drop in groundwater levels per year. Finally, the resultant annual safe yield of the aquifer is estimated as 0.84 MCM.     

基于局域化EnKF同时同化地下水水流和溶质运移数据的渗透系数场识别研究

局域化改进集合卡尔曼滤波（EnKF）可以克服EnKF方法在使用小集合时,对参数识别精度较低的缺陷,其能同化 地下水位观测数据有效识别渗透系数场。实际工作中,溶质运移数据也较容易获得。崔凯鹏（2013）尝试增加溶质运移 数据以改进只同化水流数据对渗透系数的估计结果,但是精度提高有限。本文在其基础上修改模型,进一步增加溶质注 入井,探究同时同化水头和溶质运移数据,对渗透系数场识别效果,之后对比了局域化EnKF与非局域化EnKF参数识别结 果,并分析了溶质影响范围与参数识别的关系。结果表明：同时同化溶质运移和水头资料,比同化单一种类观测数据识别 的渗透系数精度更高;相同实现数目下,局域化EnKF比EnKF对渗透系数场的估计结果与真实场更为接近;仅考虑溶质影 响范围内的渗透系数,同化水头数据在最后时刻参数识别结果好于同化溶质运移数据参数识别结果,但差别不大。     

Estimating aquifer thickness using multiple pumping tests

A method to estimate aquifer thickness and hydraulic conductivity has been developed, consisting of multiple pumping tests. The method requires short-duration pumping cycles on an unconfined aquifer with significant seasonal water-table fluctuations. The interpretation of several pumping tests at a site in India under various initial conditions provides information on the change in hydrodynamic parameters in relation to the initial water-table level. The transmissivity linearly decreases compared with the initial water level, suggesting a homogeneous distribution of hydraulic conductivity with depth. The hydraulic conductivity is estimated from the slope of this linear relationship. The extrapolation of the relationship between transmissivity and water level provides an estimate of the aquifer thickness that is in good agreement with geophysical investigations. The hydraulically active part of the aquifer is located in both the shallow weathered and the underlying densely fractured zones of the crystalline basement. However, no significant relationship is found between the aquifer storage coefficient and initial water level. This new method contributes to filling the methodological gap between single pumping tests and hydraulic tomography, in providing information on the variation of the global transmissivity according to depth. It can be applied to any unconfined aquifer experiencing large seasonal water-table fluctuations and short pumping cycles.