含水层和土壤的随机特征对水分运动的影响

将土壤或含水层的介质参数视为随机分布,建立了饱和非饱和水流运动的随机微分方程,利用随机数值求解方法,得到各种随机参数对水头分布和含水量空间变异性的影响。分析结果表明,随着渗透系数空间变异性和相关尺度的增大,含水层水头分布的变异性增大;对于区域内有抽水井分布的问题,在井点附近区域,z向水流速度空间变异性大,而在井点处变异性较小。当采用Gardner-Russo模型描述非饱和水分特征时,其非饱和土壤介质参数α的变异性对土壤负压和含水量变异性的影响大于饱和渗透系数,土壤含水量越小,含水量的空间变异性越大。初步讨论了根据随机数值方法研究地下水运动问题的可靠性分析方法。     

库水位上升条件下边坡渗流场模拟

库水位的上升将导致边坡体内渗流场的变化,进而引起边坡的失稳破坏。边坡体内的渗流场是饱和渗流与非饱和渗流共同作用的结果,通过对库水位上升过程中坡体内地下水位线的位置变化及相应的渗流场参数的计算表明:由于上部松散堆积体与基岩渗透系数的差异,基岩内地下水位线的抬升速度明显滞后于上部松散堆积体;体积含水量、压力水头、流速和水力梯度等渗流参量在基岩中的变化幅度较小,而在上部松散堆积体中的变化幅度则相对较大。这说明库水位上升对上部松散堆积体的影响相对来说还是比较明显的。     

泄洪雾雨作用下边坡饱和非饱和非稳定渗流研究

泄洪雾雨作用下的边坡岩体饱和非饱和渗流问题是目前我国西南地区在建和已建高坝所面临的一个共同问题。从理论上分析了雾雨作用下饱和非饱和介质中水分运动特征;根据饱和渗流和非饱和渗流的数学模型,建立起统一的岩体饱和非饱和非稳定渗流数学模型;编制了有地表入渗作用的饱和非饱和非稳定渗流三维有限元计算程序;对有限元求解中的容水度问题、非饱和水力参数问题以及初始水头场问题作了相应的优化处理;以实际工程为例进行了饱和非饱和渗流分析,计算结果表明,比天然降雨雨强大很多的雾雨入渗会形成对边坡稳定不利的暂态饱和区,并引起地下水位的抬高;暂态饱和区的大小和地下水位的增幅取决于边坡表面的护坡效应和排水、防渗帷幕等措施,因此应该加强边坡表面的保护以及合理的设置排水、帷幕等措施,以减少入渗和地下水位的抬高。     

基于EnKF综合水头和浓度观测数据推估地下水流模型参数

地下水反应运移模型具有参数个数众多,观测数据类型多样的特点。为了探究不同类型观测数据在反应运移模拟数据同化中的数据价值,构建了三氯乙烯降解反应运移模型的理想算例,基于水头和浓度两种类型观测数据,采用集合卡尔曼滤波方法推估渗透系数和贮水系数的非均质空间分布,讨论了影响同化结果的因素。结果表明:与仅同化水头数据的结果相比,联合同化水头和浓度观测数据推估渗透系数场和贮水系数场时具有更高的精度,在观测数据拟合和模型预测方面也有更好的表现。与目前溶质运移模型、非饱和流模型等地下水模型中的研究结果相似,数据同化结果受样本数量,观测井的数量和位置的影响,合理优化布置监测井和选择样本数量可有效改善数据同化效果并提高计算效率。     

多孔介质随机参数的相关性对非饱和水流的影响

基于van Genuchten-Mualem非饱和水分特征模型,建立了非饱和流运动的随机数值模型。将饱和水力传导度和孔隙大小分布参数视为服从对数正态分布的随机场,用Karhunen-Loeve展开分解,水头表示为混沌多项式展开。通过摄动方法得到一系列关于水头展开的偏微分方程,并用有限差分法进行求解。应用本文的模型分析了两随机场在统计不相关和完全相关模式下对水流随机分析的影响,结果表明两种模式下的水头均值相同,完全相关模式下的水头标准差较不相关模式下的明显偏小。     

基于流动电位正演模型的非饱和带水流过程监测

为了获取非饱和带水流过程的信息,借助流动电位正演模型,通过数值实验探讨非降雨和降雨两种条件下非饱和带流动电位和水流过程的关系,然后用南京中山植物园试验场地野外观测的流动电位和张力数据加以对比和验证。野外试验表明:流动电位可以有效地反映非饱和带水流过程。在夏季无降雨入渗的条件下,日周期变化的地表地下温度差导致水分的运动,流动电位准确地指示了非饱和带含水量和毛细压力的变化情况,从而指示出了水分运移的方向;在夏季有降雨入渗的条件下,降雨锋面推进之处,含水量和流动电位同时有明显的响应,进而根据不同位置的流动电位对降雨入渗响应的时刻差,直接求出入渗锋面的推进速度。     

岩体裂隙网络渗流变水温影响分析

考虑岩体裂隙渗流变水温影响,推导单裂隙变水温水流近似解析解和有限元解,在此基础上分别建立裂隙二维网络变水温渗流数值求解方程,分别对应裂隙网络变水温渗流分析的近似解析法和子结构法。分析变温水流运动规律发现：（1）单裂隙内水流水头与水力坡降成非线性关系,当水流由高温区向低温区流动时,水头分布曲线为凸曲线,此时按线性渗流简化水头整体偏小;当由低温区向高温区流动时,水头分布曲线为凹曲线,此时按线性渗流简化水头整体偏大。（2）单裂隙内,高水温处水力坡降小,低水温处水力坡降大;裂隙平均水温越高,流速越快;裂隙网络内存在与裂隙宽度相似的温度偏流效应,即交叉节点水流有偏向水流温度高的裂隙流动的趋势。在温度较高和温度梯度较大的区域,应该考虑水流温度变化对渗流场的影响。     

单环入渗试验与数值反演法结合推求土壤水力参数

单环入渗试验广泛应用于野外测定土壤水力参数。文章在三维非饱和带水分模拟程序VSF基础上,加入了单环入渗水流计算模块,并将改进的VSF与参数优化程序LEVMAR相耦合,构建了适用于单环入渗试验的土壤水力特征参数反演程序VSF-LEVMAR。利用大型室内土槽单环定水头入渗及变水头入渗试验资料,反求水平及垂向饱和入渗系数等土壤水力参数。应用结果表明:该方法既可以利用单环定水头入渗试验单独推求吸湿过程中的土壤水力参数,也可以结合单环定水头入渗与变水头入渗试验同时推求吸湿及脱湿过程中的土壤水力参数。     

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

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

均质非饱和土一维稳态流水头垂直分布的解析解

渗流场水头分布计算是进行渗流量和渗流水力坡降计算的基础,准确、有效地求取渗流场水头分布是渗流计算的关键环节。对均质非饱和土体一维稳态流的流动方程进行分析,考虑到渗透系数是与基质吸力相关的函数,通过数学变换,给出了稳定渗流场的解析通式,并基于渗透性函数中的Gardner模型,给出了非饱和土一维稳态流水头垂直分布的解析解。该解析通式表明,均质非饱和土一维稳态流水头垂直分布主要受地表水头、深度和流动率3个因素控制。分别计算了一维稳态蒸发条件下粉土和黏土两种典型土类水头沿垂直方向的分布。计算结果表明：稳态蒸发条件下粉土层和黏土层内的水头分布表现出相似的变化规律,即自地表至地下水位处随着土层深度的增加,水头分布呈现出加速递减的趋势;在相同的蒸发条件下,对于相同深度处的黏土和粉土而言,黏土层内水头更高些;对同一种土类而言,在较大的蒸发状态下同一深度处土层内水头更高。反之,则较低。     

降水入渗补给过程中优先流的确定

优先流是降水、灌溉水等入渗补给地下水的主要形式之一, 流速快, 流动路径复杂, 难以定量描述.针对优先流难以定量描述的问题, 以郑州地中渗透仪观测资料为基础, 探讨了新乡亚砂土等试筒降水入渗过程及其中的优先流补给量比例.根据土壤的水力性质、气候等资料建立不存在优先流的数值模拟模型来刻画降水入渗补给过程, 通过模拟获得的地下水入渗补给量与实测地下水入渗补给量的历时曲线, 将大于模拟值的实测值视为优先流的量及确定其在总补给量中所占的比例.结果表明, 优先流占总补给量的比例在10%~80%之间; 随着土壤粘性增加, 优先流所占比例呈增加趋势; 随地下水位埋深的增大, 优先流所占比例呈逐渐下降趋势.      

Numerical simulation of groundwater flow for a coastal plain in Japan: data collection and model calibration

Using a three-dimensional finite element model, this study characterizes groundwater flow in a costal plain of the Seto Inland Sea, Japan. The model characterization involved taking field data describing the aquifer system and translating this information into input variables that the model code uses to solve governing equations of flow. Geological geometry and the number of aquifers have been analyzed based on a large amount of geological, hydrogeological and topographical data. Results of study demonstrate a high correlation between the ground surface elevation and the groundwater level in the shallow coastal aquifer. For calibrating the numerical groundwater model, the groundwater flow was simulated in steady state. In addition, the groundwater level and trend in the transient state has also been elucidated. The numerical result provides excellent visual representations of groundwater flow, presenting resource managers and decision makers with a clear understanding of the nature of the types of groundwater flow pathways. Results build a base for further analysis under different future scenarios.     

GIS-based spatial and temporal investigation of groundwater level fluctuations under rice-wheat ecosystem over Haryana

Groundwater irrigation is the most predominant method used across India and about 50% of the total irrigated area is dependent on it. The state of Haryana has witnessed a spectacular increase in agricultural production in the last few decades, and is largely dependent on groundwater for irrigation. Groundwater mining for irrigation has become a reality in the state over the years as the number of tube wells has increased from 0.02 million in 1966 to 0.73 million in 2012, showing alarming signs of over-exploitation. The impact of increased groundwater irrigation on groundwater levels has not been studied both spatially and temporally. Therefore, this study has been undertaken to investigate the groundwater level fluctuations in the state using geographical information system (GIS) from the groundwater level data of 893 observation wells obtained from Groundwater Cell, Department of Agriculture, Government of Haryana, Panchkula, for the period 2004-12. Many researchers have applied GIS to reveal the spatial and temporal structure of groundwater level fluctuation and as a management and decision tool. The analysis of results indicated a mix of negative and positive trends in the groundwater levels. However, the negative trends were much more pronounced than positive ones. Groundwater level in the state ranged between 0.16 to 65.97 m from the ground surface and per cent area with groundwater level depth more than 10 meter (critical category) was about 56% in 2004 and has increased to 64% by 2012. The average annual decline in groundwater level was observed to be above 32 cm/year, with the strongest decline (108.9 cm/year) in Kurukshetra district.     

Groundwater under stress: the importance of management

The most significant water problem facing the United States and the world is the scarcity of water. Scarcity is expected to intensify during the twenty-first century due to global population growth, economic growth and the need to protect environmental assets. In many regions, climate change may cause scarcity to become more acute. Groundwater has the capacity to buffer extreme hydrologic events because recharge is not directly tied to precipitation and run-off events. Groundwater storage offers many opportunities to increase storage and enhance the general water supply. The value of groundwater will increase as scarcity intensifies and effective groundwater management will be required if groundwater values are to be reaped. Effective groundwater management will require attention to certain well-established economic principles. Persistent overdraft is always self terminating. Individualistically competitive modes of extraction are inefficient. It is virtually always cheaper to prevent pollution in the first place than to remediate it. Groundwater management schemes appear to be most effective when they are locally developed and managed. Such schemes can regulate extractions and/or recharge either directly or indirectly. Evidence from the United States suggests that direct regulations can be made to work if they focus on recharge and developing supplementary sources of water. Evidence from South Asia suggests that regulation of extractions and provision of recharge water is most effective when indirect measures are used in which groundwater and ancillary variables such as electricity are co-managed. This evidence suggests that as supplementary sources of supply become scarce, innovative and indirect groundwater management schemes are likely to become more common.     

Estimation of unrecorded groundwater abstractions in Jordan through regional groundwater modelling

Jordan suffers from water scarcity and groundwater covers the majority of Jordan’s water supply. Therefore, there is an urgent need to manage this resource conscientiously. A regional numerical groundwater flow model, developed as part of a decision support system for the country of Jordan, allows for quantification of the overexploitation of groundwater resources and enables determination of the extent of unrecorded agricultural groundwater abstraction. Groundwater in Jordan is abstracted from three main aquifers partly separated by aquitards. With updated geological, structural, and hydrogeological data available in the country, a regional numerical groundwater flow model for the whole of Jordan and the southernmost part of Syria was developed using MODFLOW. It was first calibrated for a steady-state condition using data from the 1960s, when groundwater abstraction was negligible. After transient calibration using groundwater level measurements from all aquifers, model results reproduce the large groundwater-level declines experienced in the last decades, which have led to the drying out of numerous springs. They show a reversal of groundwater flow directions in some regions, due to over-abstraction, and demonstrate that documented abstractions are not sufficient to cause the observed groundwater-level decline. Only after considering irrigation water demand derived from remote sensing data, the model is able to simulate these declines. Illegal abstractions can be quantified and predictive scenarios show the potential impact of different management strategies on future groundwater resources.

     

乌鲁木齐河流域柴窝堡盆地与河谷区地下水流模拟

乌鲁木齐河流域柴窝堡盆地与河谷区是乌鲁木齐市的重要供水水源地。为了更好地为地下水资源的可持续开发提供依据,本文尝试将几何形状及含水层厚度差异悬殊的柴窝堡盆地与河谷区联合起来,运用GMS模型系统建立了区域地下水稳定流和非稳定流数值模型。通过对观测孔地下水位过程线和地下水流场进行拟合,校正了研究区的渗透系数、给水度和储水率等水文地质参数。并运用模型分析了研究区的地下水流系统、地下水均衡量、地下水位变化趋势及地下水位监测网的设计。结果表明研究区地下水流主要存在两个径流排泄区:沿乌鲁木齐主河道的径流排泄区和以柴窝堡湖和大小盐湖为排泄中心的径流区;地下水储存量相对较大,能够调节季节性变化与地下水开采的影响;应加强主要河流洪积扇补给区和山前侧向补给带地下水位的监测与排泄区泉水流量的观测。     

Structure and parameter schematization of the vadose zone for groundwater recharge modeling

On the basis of published soil profile data typical retention curves and hydraulic conductivity functions were obtained that depend on soil depth, lithological characters, and different surface landscape conditions. The corresponding average soil hydraulic parameters for M.Th.Van Genuchten equations were estimated. These average soil parameters are suitable for groundwater recharge estimations based on unsaturated flow modeling.     

Numerical simulation of groundwater flow and contamination transport in Nahavand plain aquifer,west of Iran

Numerical simulation of groundwater flow used for the estimation of hydraulic and hydrologic parameters which is an important tool for the management of aquifers. This study presents the results of a mathematical model developed for the simulation of groundwater flow in Nahavand plain aquifer in the southwest Hamadan province. For this purpose Groundwater Modeling Software (GMS) was used which supports the MODFLOW-2000 code. After gathering required data such as the hydrological, hydrogeological and topography maps, a 3D hydrogeological model of plain was constructed with borehole and surface elevation data. Then MODFLOW was used for simulation of flow. After initial simulation of the flow, the model was calibrated in steady state with trial-and-error and parameter estimation methods the observed head of groundwater table monitoring data of 1997. Results of calibration show that error between observed head and computed head is in allowable range. Also results of computed head with model show that groundwater flow is in the direction of the dominate slope (southeast to northwest). Finally MODPATH code which simulates advective transport of particles was used for estimation of flow path and source of contaminants.     

Environmental impacts of groundwater abstraction in Neresnica brook catchment (Slovak Republic)

The main goal of the research was a detailed study of the consequences, caused by groundwater abstraction in the Neresnica brook catchment (Slovak Republic) since the water supply source Podzamcok was put in operation in 1973. This goal was supplemented by research of groundwater origin, groundwater chemical composition and its changes. Field measurements, analyses of water samples and statistical evaluation of hydrometeorological data were used to gain and process the data. It was showed that overexploitation of groundwater in the catchment resulted in distinct diminishing of stream flow discharges in low flow periods, groundwater level dropping and disturbance of hydraulic connection between surface and groundwater. Groundwater quality was quite stable during the whole period of observation. Groundwater modeling was proposed to be used for the re-estimation of utilizable groundwater amounts in the catchment.     

Uncertainty analysis of groundwater model based on data assimilation北大核心CSCD

黑河流域中下游地下水系统受上游冰冻圈融水和降雨的补给,由气候变暖导致的冰冻圈萎缩致使中下游地下水系统的稳定性面临更多的风险。地下水模型是地下水系统稳定性评估的有效手段,但是地下水模型参数往往存在较大的不确定性。为此,本文提出了基于数据同化算法的不确定性分析方法,通过包含观测资料信息减小模型不确定性。采用所提方法分析了(基于MODFLOW构建)黑河流域中游地下水模型中13个参数的不确定性,讨论了算法超参数的影响及其最优取值,分析了地下水模型参数的不确定性。实验结果证明数据同化算法可有效减小地下水模型参数的不确定性,观测资料的种类与数量对参数不确定性的减小起到重要作用;不同地下水模型参数的不确定性不同,地表水与地下水相互作用频繁的区域参数不确定性较大;含水层渗透系数、含水层给水度以及灌溉回流系数对模型输出的地下水位输出影响显著,河床水力传导系数对模型输出的河流流量影响较大。本研究将为地下水研究提供更加可靠的模型方法,为西北内流区地下水哺育的绿洲生态系统稳定可持续研究提供重要支撑。