改变入渗强度的地下水流模式实验

自Tóth提出地下水流系统理论以来,大多学者基于定水头上边界的解析解或数值模拟研究地下水流系统发育特征,尤其是多级水流系统模式。认为复杂的盆地"地形势"就会发育多级地下水流系统,简单地形发育单一地下水流系统。利用自主研发的地下水流系统砂槽模型进行系列实验发现,实验条件下(模型尺寸、介质和多个河谷不变)水流系统的发育受控于降水入渗强度。得出随着降水入渗强度的增大,水流系统模式由单一的区域水流系统,到复杂的多级次水流系统,再到局部的水流系统的结论;盆地多个势汇中,只有实际成为地下水排泄点的汇势才能影响地下水流系统模式,而其他的势汇只是可能的潜在势汇。实验表明,用定流量上边界(降水入渗强度)能够更好地揭示盆地地下水流系统模式发育与变化机理。     

通量上边界与水头上边界方法的地下水流系统模拟对比

水头上边界方法（简称水头法,GHB）给定了潜水面形态和固化了排泄点,限制了不同地下水流系统模式的形成与转化。分别用数值法进行了通量上边界（简称通量法,FUB）与水头上边界的地下水流系统对比模拟。结果表明：嵌套式多级地下水流系统（Tóth典型模式）在运用水头法和通量法进行系统转化模拟时,得出的水流模式可能相似或完全不同;通量法在条件（盆地形态、入渗强度等）改变时潜水面能够自动形成,从而得出不同变化条件下的水流系统特征;水头法由于给定了潜水面和固化了排泄点,在改变盆地其他因素时,盆地补给（排泄）也发生同步改变,此时地下水流模式不是单因素变化的结果,因此在给定条件下不能得出完整的地下水流系统变化模式。基于通量法与水头法在地下水流系统模拟中的优势与不同,在进行盆地地下水流系统理论和实际研究时,应该综合2种方法的特点,结合实际资料条件进行方法的选取与应用。     

不同地下水流系统模式渗流场和温度场的互相影响

在地热资源丰富的地区,需要研究不同地下水流系统发育模式下渗流场和温度场的互相影响.基于二维潜水盆地多源汇的数值模拟和室内砂箱实验,改变降雨入渗强度,通过砂箱底部加温研究上下边界不同温度差条件下的渗流场和温度场的变化.研究结果表明:①随着降雨入渗强度加大,地下水流速增大,地下水流系统由单一区域系统向复杂的局部+区域、局部...     

通量上边界渗透系数随埋深增加指数衰减的地下水流系统

地下水流系统理论是当代水文地质学的核心概念框架,研究不同控制因素对盆地地下水流系统发育模式的影响具有重要意义。近年来的研究表明利用通量上边界能够更好地揭示盆地不同要素对地下水流系统模式转化的影响。基于通量上边界,采用数值模拟方法,研究稳定流条件下,渗透系数随埋深呈指数衰减的非均质含水层对盆地地下水流系统模式转化的影响。结果表明:随着渗透系数随埋深指数衰减程度的加大,盆地潜水面整体抬升,盆地上部地下水流速增大,水力梯度增大;同时,盆地地下水流系统由复杂的多级次水流系统到单一的局部水流系统,顺向局部水流系统占据的空间消减,而逆向局部水流系统占据的空间增大,流速近似为零的局部滞留区域向左下方移动。     

盆地地下水流系统形成与影响因素分析

目前区域（盆地）地下水流系统模拟研究中,常用的定水头与通量两种上边界条件刻画方法与实际条件存在差距。通过对比两种方法的差异和各自适用条件,采用解析法讨论地下水位的形成控制机制,提出了改进后的变通量上边界数值模型,并以鄂尔多斯盆地北部白垩系地下水流系统为例分析了盆地地下水流系统的形成与影响因素。研究表明,鄂尔多斯盆地北部白垩系水流系统地下水位受地形、补给条件和渗透系数三者共同控制,同时特有的气候、地形和岩性组合通过控制地下水位影响地下水流系统的发育演化。采用变通量上边界法探讨上边界条件改变对盆地水流系统的影响,对深刻认识区域地下水流系统形成演化机制,揭示地下水系统与上边界气候变化、植被生态变化之间的相互作用关系具有一定优势。     

氢氧同位素在地下水流系统的重分布：从高程效应到深度效应

大气降水的氢氧同位素含量具有高程效应,降水入渗后参与地下水循环,其高程效应如何受地下水流系统的影响转化为地下水氢氧同位素的深度效应？现有研究对于这个问题缺少定量认识。文章构建单向倾斜盆地和双峰波状盆地的稳态地下水循环理论模型,采用MODFLOW模拟剖面二维地下水流场、采用MT3DMS模拟重同位素分子的对流-弥散过程,得到地下水D和¹⁸O含量的空间分布,探讨了氢氧同位素高程效应在地下水流系统转化为深度效应的机理。结果表明：在单斜盆地,补给区大气降水D和¹⁸O含量的高程效应转化为排泄区地下水δD和δ¹⁸O值随埋深增大而指数型衰减的深度效应;在双峰波状盆地,当含水层渗透性相对入渗强度较大时(K₀/w=1 000),仅发育一个区域地下水流系统,在区域地下水的排泄区δD和δ¹⁸O随埋深增大呈现S形曲线分布;当含水层渗透性相对入渗强度较小时(K0/w=250),双峰波状盆地发育多个局部地下水流系统,区域地下水的排泄区δD和δ¹⁸O随埋深增大呈现S形曲线,而局部地下水排...     

基于GMS的岩溶山区三维地下水流模式识别

以高岚河东部多条支沟切割形成的连排型河间地块为研究区,基于1∶5万岩溶水文地质调查成果,利用GIS技术和径流分割,量化岩溶含水系统顶底板高程、获取入渗补给系数及地下水径流量等数据,运用三维地下水流数值模拟及粒子追踪地下水流线的方法,得出了研究区的地下水流模式。其结果表明:研究区地下水流模式主要受区域尺度的势源和势汇控制,局部的地形起伏对地下水位的影响不明显;支沟切割较深的凉伞沟、滩淤河流域,水流驱动力大,有利于局部水流系统的发育;随着东部补给高程的增加,流程逐渐增长,多发育局部水流系统,至山脊附近多发育排泄至高岚河的中间水流系统。当研究区年降雨量由中等的1 021.1 mm降为多年最低的725.5 mm时,由东部补给区排泄至高岚河的中间水流系统增加;在此降雨强度区间,不具备发育越过河间地块的中间或区域地下水流系统。      

红碱淖—沟岔地区地下水流系统研究

本文根据对红碱淖－沟岔地区地下水形成的背景条件、水化学特征、环境同位素组成的分析和剖面地下水势的计算，把研究区划分为２个区域地下水流系统和３个局部地下水流系统，对不同等级地下水流系统的边界、系统的循环作了详细探讨，深刻揭示了这种复合地形区、均质含水介质中地下水流的循环模式。     

盆地多级次地下水流系统盐分运移实验模拟

地下水年龄和滞留时间包含了地下水循环和演化的重要信息,被广泛用于盆地地下水循环模式的研究.使用多级次地下水流系统演示仪,实验模拟了三级水流系统模式中地下水年龄及滞留时间分布.研究发现,盆地底部、区域流线的下游、盆地滞留区最晚响应;浅部的局部水流系统稳定后的浓度值相对较低;中间水流系统相对深部的区域水流系统也较低,滞留区...     

基于砂槽模型研究不同水流密度下盆地地下水流系统

在干旱半干旱地区,沉积盆地是人类生产、生活用水的主要载体,为查明在水流自身密度改变条件下盆地地下水流系统的变化特征,需在变密度条件下对盆地地下水流系统进行模拟研究。本文采用砂槽物理模型,在改变区域性水流密度条件下模拟研究地下水流系统的变化规律。试验结果显示:随着区域性水流密度的增加,区域水流系统的渗流速度和循环量均减少,且流线径流距离和径流深度也随之减小;而局部水流系统的渗流速度和循环量均增加。结果表明区域性水流密度的增加会抑制区域水流系统的发育,对局部水流系统有一定增强作用。本研究着重强调了区域性水流自身密度变化对盆地地下水流系统所产生的影响,并通过物理试验模拟得出在改变区域性水流密度条件下盆地地下水流系统的变化规律。     

鄂尔多斯白垩系地下水盆地水文地质概念模型

水文地质概念模型是地下水模拟与水资源评价的基础。为了构建地下水数值模拟模型,以鄂尔多斯白垩系地下水盆地为研究对象,以丰富翔实的勘查资料为依据,从白垩系地下水系统边界条件、水文地质结构、水流系统、水文地质参数、源汇项等方面进行概化,联合采用ArcGIS、RS、GMS同位素、基流分割、地质统计学、随机模拟等技术方法,建立了一组水文地质概念模型及其相应的空间离散表达,并讨论了地下水模拟面临的问题与发展趋势。     

Groundwater modeling applied to production well management in a contaminated aquifer

A study was conducted to evaluate production strategies for a well field system near a source of groundwater contamination. Numerical modeling of groundwater flow was employed to generate hydraulic head configurations for different production scenarios. For a given scenario, an evaluation of contamination susceptibility was made by comparing head distributions in two aquifer units to the positions of the contaminant source and discharging water supply wells. The results of this study suggest that groundwater flow modeling can be a useful technique for planning the production of water supply wells in aquifers at risk of contamination from anthropogenic pollution sources.     

地下水流系统理论与研究方法的发展

地下水流系统理论的提出,推动了现代水文地质学的发展。以Tóth经典地下水流系统理论建立方法为基础,综述了基于Tóth方法的地下水流系统的模拟成果,分析了Tóth方法得出的水流模式与控制因素的关系,以及地下水流系统理论从概念到实际应用的发展。同时,对中国地质大学(武汉)提出的地下水流系统通量上边界模拟方法也进行了系统论述,在物理模拟实验与数值模拟基础上,认为通量上边界分析方法是对Tóth方法的改进与完善,有利于对地下水流系统发育的物理机制理解;该方法能够更全面认识地下水流系统模式及其转化,定量出各影响因素对地下水流模式的控制关系。最后指出地下水流系统理论是当代水文地质学的核心概念框架,应该重视地下水流系统理论物理机制和数学模拟方法的研究,加强新技术方法的引入,拓宽其应用领域的研究等。     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

An analytical study on groundwater flow in drainage basins with horizontal wells

Analytical studies on release/capture zones are often limited to a uniform background groundwater flow. In fact, for basin-scale problems, the undulating water table would lead to the development of hierarchically nested flow systems, which are more complex than a uniform flow. Under the premise that the water table is a replica of undulating topography and hardly influenced by wells, an analytical solution of hydraulic head is derived for a two-dimensional cross section of a drainage basin with horizontal injection/pumping wells. Based on the analytical solution, distributions of hydraulic head, stagnation points and flow systems (including release/capture zones) are explored. The superposition of injection/pumping wells onto the background flow field leads to the development of new internal stagnation points and new flow systems (including release/capture zones). Generally speaking, the existence of n injection/pumping wells would result in up to n new internal stagnation points and up to 2n new flow systems (including release/capture zones). The analytical study presented, which integrates traditional well hydraulics with the theory of regional groundwater flow, is useful in understanding basin-scale groundwater flow influenced by human activities.     

Hydrogeological model of the Baltic Artesian Basin

The Baltic Artesian Basin (BAB) is a complex multi-layered hydrogeological system in the south-eastern Baltic covering about 480,000 km². The aim of this study is to develop a closed hydrogeological mathematical model for the BAB. Heterogeneous geological data from different sources were used to build the geometry of the model, i.e. geological maps and stratigraphic information from around 20,000 boreholes. The finite element method was used for the calculation of the steady-state three-dimensional (3D) flow of unconfined groundwater. The 24-layer model was divided into about 1,000,000 finite elements. A simple recharge model was applied to describe the rate of infiltration, and the discharge was set at the water-supply wells. Variable hydraulic conductivities were used for the upper (Quaternary) deposits, while constant hydraulic conductivity values were assumed for the deeper layers. The model was calibrated on the statistically weighted borehole water-level measurements, applying L-BFGS-B (automatic parameter optimization method) for the hydraulic conductivities of each layer. The principal flows inside the BAB and the integral flow parameters were analyzed. The modeling results suggest that deeper aquifers are characterized by strong southeast–northwest groundwater flow, which is altered by the local topography in the upper, active water-exchange aquifers.     

Using lithologic modeling techniques for aquifer characterization and groundwater flow modeling of the Sohag area, Egypt

Three dimensional lithologic modeling techniques have been used for detailed characterization and groundwater flow modeling of the Quaternary aquifer system of the Sohag area, Egypt. Well log data were used for building the lithologic model using RockWorks. A groundwater flow model, facilitated by MODFLOW 2000, was built using results of the lithologic model. The obtained lithologic model honored the well log data and revealed a complex sedimentary system, which is mainly composed of six lithologic categories: clay, clay and sand, fine sand, coarse sand, sand and gravel, and gravel. Inter-fingering and presence of lenses are the main characteristics of the sedimentary basin in the study area. A wide range of hydraulic conductivities is present, which vary spatially and control the groundwater flow. Heterogeneity of the aquifer system is spatially represented where different hydraulic conductivities are found in the different directions. Sandy layers tend to be connected. Hydraulic continuity is represented by inter-fingering and connection of sandy materials within the aquifer system.