常规越流承压含水层影响半径计算方法之改进

影响半径是一个有争议的概念,针对不同含水层和计算模型各有不同的计算方式,而在工程应用中常使用理论依据不甚明确的经验公式。针对可形成稳定流的越流承压含水层,展开对确定其影响半径的常用方法即3个及以上观测孔稳定流阶段降深-距离直线图解法的分析,认为这一方法容易造成工程师对影响半径的误判,指出其易于误读之处并对使用该法提出了修正建议。     

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

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

越流承压含水层不稳定流问题有限差分法的稳定条件

证明了一维承压越流含水层系统不稳定流问题差分格式的稳定性条件,显式差分格式的稳定性条件是△ｔ（４Ｔ／（μ（△ｘ）＾２）＋ｂ／μ）≤２隐式格式是无条件稳定的,以河间地块一维承压越流含水层的为例,根据给出的稳定性条件进行了数值计算,结果表明,在满足稳定性条件时,对计算结果符合实际。     

海潮引起有越流的滨海承压含水层地下水头变化的数值模拟

海潮波动可以引起海岸带有越流的承压含水层地下水头发生波动。建立了基于有限差分法的滨海地区一维承压含水层地下水运动数值模型。通过将潮汐波动概化为正弦波,分别对初始水头水平及线性倾斜的承压含水层模拟了滨海地区有越流的承压含水层地下水头随潮汐波动的变化。通过对两种情形下的变化比较,结果表明,受海潮影响的滨海承压含水层地下水头与海潮有相似的波动特征,但变幅减小,受海潮的影响程度与离海岸的距离有关,随着离海岸距离的增加,地下水头的变幅及潮汐效率呈负指数函数衰减,水头倾斜情形下变幅更小,潮汐效率更小,滞后时间更短,地下水头对海潮的滞后时间随距离呈线性增加。     

承压含水层井流-盖层弯曲效应的解析理论

指出承压含水层盖层的弯曲变形与开采井周围的径向地下水运动存在相互作用, 而这一效应在传统的井流理论中没有被认识到.通过引入弹性薄板理论, 建立了无越流的承压含水层井流-顶板弯曲效应的解析模型, 同时考虑了含水层和水的压缩性, 结果表明Theis井流方程给出的抽水降深偏小.在此基础上推导了有越流承压含水层井流-盖层弯曲效应的偏微分方程, 求出了解析解, 并与传统理论的结果进行了对比, 表明Hantush-Jacob公式计算的降深也是偏小的.在抽水井附近和抽水初期, 传统理论可能导致显著的相对误差.      

华北平原深层淡水在开采条件下接受上覆咸水越流补给——以天津平原为例

分析物探测井曲线和地下水位动态发现，天津平原深层孔隙地下淡水在开采条件下接受上覆咸水的越流补给，而且咸水底界未大量向下位移。本文讨论了咸质含水组底部钙化粘土微观粘粒与水溶液的相互作用，提出了“滤盐层”的概念。     

越流含水层中抽水井附近非达西流两区模型近似解析解

构建了越流含水层中抽水井附近非达西流动的两区模型,即距离抽水井较近的区域为非达西流,而相对较远区域为达西流,两区之间的临界半径可根据临界雷诺数确定.采用线性化近似方法和Laplace变换相结合分别得到了非达西流区域和达西流区域的水位降深在拉氏空间下的解析解,应用数值Laplace逆变换-Stehfest方法得到其在实空...     

非线性多含水层地下水越流系统水资源优化管理模型研究

提出了非线性多含水层地下水越流系统水资源管理的方法。内容包括非线性系统水资源管理的基本理论，算法原理，程序软件设计和理论检验，并把这一方法应用到唐山市以防治岩溶地面塌陷为目的水资源管理中。     

用水位恢复数据反演越流承压含水层参数

水文地质勘测中,由于水位恢复阶段没有人力和机械因素干扰,由测量数据可以画出平滑的曲线,更适宜于分析水文地质参数,但经典非稳定流Theis公式和Theis水位恢复法所依据的若干假设条件在实际中难以满足。针对越流承压含水层具有补给条件而很快稳定以及非稳定抽水阶段流量常“抖动”的特点,论证了利用稳定流量和Hantush-Jacob越流模型计算的可行性,并用扩展的卡尔曼滤波器反演水文地质参数,经文献中历史数据和实地抽水试验检验结果表明,由该方法反演的预测值和实测值吻合良好。     

承压含水层非完整反滤回灌井的稳定流计算

反滤回灌井是地下水库中一种有效的回灌设施,它由回灌井和井口反滤层组成。非完整反滤回灌井井流运动由砂反滤层的竖向流和非完整回灌井井流组合而成。文中提出了承压含水层非完整反滤回灌井稳定流的计算方法,分析了反滤层、井损和回灌堵塞等因素对回灌量的影响,并结合实例证明了反滤回灌井回灌量大大低于普通回灌井的回灌量。     

不同地连墙插入深度下承压含水层减压降水对既有隧道的影响

当基坑需进行未截断条件下的承压含水层降水时,承压水抽降产生的影响范围远大于基坑开挖的影响范围,其对基坑外既有隧道的影响值得重视。进行了承压层减压降水对既有盾构隧道影响的有限元仿真模拟,并分析了采用地下连续墙插入承压含水层时,不同插入深度对隧道变形的影响,结果表明：随着地下连续墙插入承压含水层深度的增加,直至完全截断承压含水层,既有隧道沉降逐渐减小,但受地下连续墙变形的影响,隧道在水平方向仍存在较大位移。因此,当既有隧道与减压井净距较小时,即使承压含水层被完全截断,也应重视减压降水对隧道水平方向变形造成的影响。此外,长期减压降水将会引起隧道产生可观的沉降。因此,应尽量避免承压含水层未被完全截断条件下进行长期减压降水。     

海潮引起初始水位倾斜的滨海承压含水层地下水位 变化的数值模拟

海潮波动可以引起海岸带地下水位发生波动。建立了基于有限差分法的滨海地区一维承压含水层地下水运动数值模型。将潮汐波动概化为正弦波,模拟了滨海地区地下水位随潮汐波动的变化。通过与初始水位水平的承压含水层水位变化的比较表明,受海潮影响的滨海承压含水层地下水位与海潮有相似的波动特征,但变幅减小,受海潮的影响程度与离海岸的距离有关,随着离海岸距离的增加,地下水位的变幅及潮汐效率呈负指数函数衰减,但比前者变化程度稍缓,地下水位对海潮的滞后时间随距离呈线性增加。     

振荡试验确定倾斜承压含水层水文地质参数

为评价振荡试验在获取倾斜承压含水层水文地质参数过程中含水层倾角对计算结果的影响,基于Kipp模型推导了考虑承压含水层倾角的振荡试验解析解模型。倾斜承压含水层量纲一贮水系数和量纲一导水系数不变条件下,测试井-含水层系统的量纲一阻尼系数随着倾角的增加而增加。修正模型计算的量纲一阻尼系数与Kipp模型计算的量纲一阻尼系数偏差10%以上时,必须考虑倾角的影响并且承压含水层倾角对水头响应曲线造成影响临界角度约为34°。修正模型和Kipp模型对现场振荡试验水头响应曲线分析结果的差异表明,忽略含水层倾角会对倾斜承压含水层导水系数和贮水系数的计算结果产生较大误差。     

Numerical modeling of groundwater inflow from a confined aquifer into Sangan open pit mine,northeast Iran

Numerical models are useful in the evaluation of the interaction between groundwater systems and mining activities. They can be successfully used to predict the quantity of inflow into open pits and to design an appropriate dewatering system. In this paper, a two-dimensional axi-symmetric finite element model called SEEP/W has been used to predict the groundwater inflow into Sangan open pit mine (anomaly north C). The Sangan iron mine is located at 280 km south-east of Mashhad, Iran, in arid and warm climate conditions wherein precipitation is generally limited. The water inflow to the pit is mainly from a confined aquifer, mainly by horizontal flow in the upper layers and vertical flow at the pit bottom. The results of the numerical model of the ground water inflow are presented and compared with those obtained from Theis, Cooper-Jacob and Jacob-Lohman analytical solutions. Ground water inflow monitoring was also carried out in a trial excavation at the Sangan mine in order to calibrate the model. The model was then used to predict groundwater inflow into Sangan open pit mine during its advancement. This model provides valuable information for designing an appropriate dewatering system.     

Conceptual model of a coastal aquifer system in northern Greece and assessment of saline vulnerability due to seawater intrusion conditions

This paper refers to the development of a conceptual model for the management of a coastal aquifer in northern Greece. The research presents the interpretation and analysis of the quantitative (groundwater level recordings and design of piezometric maps) regime and the formation of the upcone within the area of investigation. Additionally it provides the elaboration of the results of chemical analyses of groundwater samples (physicochemical parameters, major chemical constituents and heavy metals and trace elements) of the area which were taken in three successive irrigation periods (July–August 2003, July–August 2004 and July 2005), in order to identify areas of aquifer vulnerability. The study identifies the areas where ion exchange phenomena occur, as well as the parts of the aquifer where the qualitative degradation of the aquifer system is enhanced. The paper, finally, assesses the lack of any scientific groundwater resources management of the area by the local water authorities, as well as the current practices of the existing pumping conditions scheme as applied by groundwater users.An erratum to this article can be found at     

Sensitivity analysis for the EPIK method of vulnerability assessment in a small karstic aquifer, southern Belgium

Applying the EPIK parametric method, a vulnerability assessment has been made for a small karstic groundwater system in southern Belgium. The aquifer is a karstified limestone of Devonian age. A map of intrinsic vulnerability of the aquifer and of the local water-supply system shows three vulnerability areas. A parameter-balance study and a sensitivity analysis were performed to evaluate the influence of single parameters on aquifer-vulnerability assessment using the EPIK method. This approach provides a methodology for the evaluation of vulnerability mapping and for more reliable interpretation of vulnerability indices for karst groundwater resources. Received, March 1999/Revised, December 1999, February 2000/Accepted, February 2000     

A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.     

A comparative study of four schemes for groundwater vulnerability mapping in a diffuse flow carbonate aquifer under Mediterranean climatic conditions

This paper shows the results of a comparative study involving application of the vulnerability mapping methods known as AVI, GOD, DRASTIC and EPIK to a pilot carbonate massif in southern Spain, namely the Torremolinos aquifer. The main objectives of the study were to determine which methods are most suitable for diffuse flow carbonate aquifers such as in southern Spain, and to evaluate variations in the degree of vulnerability associated to the rainfall variations that normally occur in a Mediterranean climate. According to three of the above methods, the aquifer is moderately vulnerable, but the AVI method evaluated it as highly vulnerable—this, however, is improbable. The vulnerability maps reflect the great importance of geology-related parameters (mainly those concerned with lithology) and, to a lesser degree, that of the depth of the groundwater table which is related to the rainfall. After this latter parameter, it is possible to distinguish between humid and dry climatic situations; thus, vulnerability increases in a humid year, especially according to the GOD and AVI methods. In conclusion, the GOD method seems the most adequate of the methods applied in this work for vulnerability mapping of diffuse flow carbonate aquifers in the Mediterranean domains.