基于非饱和模拟的海绵体土壤蓄渗能力评价

在武汉市长江新城起步区开展了海绵体土壤蓄渗能力评价,基于10个双环入渗试验和土壤颗分资料,利用HYDRUS-1D构建了非饱和水流数值模拟模型,对2018年降雨蒸发条件下海绵体土壤中的水分运移过程进行模拟,定量计算了海绵体土壤的蓄水容量和入渗水量。结果显示：非饱和数值模拟很好的仿真了降雨入渗过程,可准确计算海绵体土壤的蓄水容量和入渗能力。     

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

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

鄂西恩施地区非饱和土降雨入渗规律

在分析鄂西恩施地区降雨特征的基础上,选取6组典型降雨条件,设计并进行了室内人工降雨土柱非饱和入渗试验。通过实测土体体积含水率,揭示了鄂西恩施地区非饱和土降雨入渗规律,获得了降雨入渗影响区体积含水率变化规律、入渗前锋运移规律以及降雨强度和历时对入渗的影响规律。通过GeoStudio软件对室内降雨试验进行了数值模拟,数值模拟结果与试验结果基本相同,进一步验证了降雨入渗规律。同时依据实测结果,分析了数值模拟的不足之处,为通过数值模拟方法研究非饱和土降雨入渗奠定了基础。     

基于FLAC3D平台的边坡非饱和降雨入渗分析

为研究降雨入渗条件下三维边坡渗流场的变化过程,在全面研究FLAC3D软件渗流分析模块功能及算法的基础上,通过编写FISH函数完善了FLAC3D软件的非饱和渗流计算功能。基于降雨入渗机制的分析,编写了降雨入渗及停雨出渗边界的FISH函数,实现边坡三维降雨入渗过程的模拟。通过降雨条件下算例边坡的饱和非饱和渗流场变化过程的数值计算,研究了降雨入渗引起的边坡暂态饱和区变化过程,并与已有研究成果进行对比分析,从而验证了自编FISH函数实现边坡三维非饱和渗流计算结果的正确性。研究成果表明,文中提出的边坡三维非饱和降雨入渗分析是可行的,同时也为三维边坡非饱和渗流场的研究提供了一条有益的途径。     

淮北平原降雨入渗补给系数随地下水埋深变化特征

水文地质参数对地下水资源评价起着至关重要的作用。其中,降雨入渗补给系数是影响浅层地下水水量、水质的重要参数。它对研究区域水量转化和水量平衡也十分重要。但是由于受降雨量、土壤类型、植被、地下水埋深等诸多因素的影响,准确判断降雨入渗补给系数存在很大困难。如果没有考虑这些因素的影响,尤其是降雨量和地下水埋深的影响,所推求的降雨入渗补给系数就会存在较大误差。结合安徽省淮北平原区五道沟水文实验站观测的降雨量、地下水补给量、地下水水位资料,利用两种不同的方法推求了不同降雨量等级的次降雨入渗补给系数。根据统计学理论研究了不同降雨量条件下,次降雨入渗补给系数随地下水埋深变化的分布规律,建立了次降雨入渗补给系数与地下水埋深的回归模型,并进行了相应的检验。研究表明,在控制地下水埋深的条件下,次降雨入渗补给系数随地下水埋深的变化符合指数分布;在地下水位自由变动的条件下符合伽玛分布。     

降雨入渗条件下非饱和土边坡稳定分析

针对降雨入渗土坡的稳定问题,建立一个考虑水力渗透系数特征曲线、土-水特征曲线以及修正的Mohr-Coulomb破坏准则的非饱和土流固耦合有限元计算模型,进行雨水入渗下非饱和土边坡渗流场和应力场耦合的数值模拟,得到非饱和土边坡变形与应力的若干重要规律。研究成果为降雨入渗条件下非饱和土边坡的稳定分析提供了基础。     

降雨入渗条件下粗粒土路堤暂态饱和区发展规律及稳定性研究

为了分析降雨入渗条件下粗粒土路堤暂态饱和区发展规律及稳定性,结合工程实际及气象资料,通过拟定不同的降雨计算方案,采用饱和-非饱和渗流理论对粗粒土路堤边坡在不同降雨条件下雨水入渗过程进行数值模拟.研究结果表明:降雨过程中路堤边坡表面会形成一定范围的暂态饱和区,暂态饱和区的形成及发展过程与降雨强度及降雨时间密切相关;降雨入渗仅会引起路堤边坡坡面以下水位线的大幅上升,对路堤中部地下水位线影响较小;降雨停止后,路堤上部暂态饱和区迅速消散,而路堤中下部暂态饱和区消散过程则表现出明显的滞后特征,但总体而言暂态饱和区的面积呈下降趋势;降雨入渗将导致入渗区域内基质吸力的丧失与含水率的增加,是造成路堤边坡稳定性降低的主要因素,路堤的失稳形式表现为整体滑移与表面局部坍塌相结合.     

降雨条件下岩土饱和-非饱和渗流分析

基于Buckley-Leverett两相渗流方程,提出了新的岩土饱和-非饱和渗流数学模型,利用有限差分方法给出了隐式压力显示饱和度数值求解方法,编制了饱和-非饱和渗流计算程序.结合工程实例,模拟了降雨入渗条件下边坡岩体渗流场孔隙压力变化和含水饱和度变化,模拟结果验证了所提出的模型对饱和-非饱和渗流的有效性.     

数值模拟在反演矿区水文地质条件中的应用

地下水数值模拟技术已经普遍应用于矿区水文地质问题的研究,以涌水量预测居多,但对于矿区水文地质条件的反演目前仍处于探索阶段.针对马坑矿区开采条件下关键水文地质条件认识不清的事实,以矿区岩溶裂隙地下水为研究对象,建立了三维地下水水流数值模型,并利用地下水位动态及涌水量资料对模型进行了识别,反演了矿区目前存在的若干重大水文地质问题,对矿区的边界条件、充水来源及途径、地表入渗等问题有了进一步的认识.该方法也是对数值模拟技术应用的一种新探索.     

降雨入渗条件下非饱和路堤变形与边坡的稳定数值模拟

对于常常处于非饱和状态的路堤在降雨入渗条件下变形问题的数值模拟研究,由于涉及到非饱和土力学,一直是个热点和难点。现有的强度折减有限元法一般仅对强度指标进行折减,但折减后得到的路堤变形场与实际变形场存在一定差异。论文分析认为,在降雨入渗影响下非饱和路堤不仅强度参数会降低,而且变形参数也会随着应力水平和强度参数的变化而变化。采用基于变模量的弹塑性强度折减法,对不同降雨入渗深度下路堤的变形规律进行数值模拟,并与室内模型试验结果进行比较;同时采用强度折减法和极限平衡法中的Bishop法,分析了不同降雨入渗深度对非饱和路堤稳定性和滑弧变化特征的影响规律。     

Influence of model conceptualisation on one-dimensional recharge quantification: Uley South,South Australia

Model conceptualisation is a key source of uncertainty in one-dimensional recharge modelling. The effects of different conceptualisations on transient recharge predictions for the semi-arid Uley South Basin, South Australia, were investigated. One-dimensional unsaturated zone modelling was used to quantify the effect of variations of (1) lithological complexity of the unsaturated zone, and (2) representation of preferential flow pathways. The simulations considered ranges of water-table depths, vegetation characteristics, and top soil thicknesses representative for the study area. Complex lithological profiles were more sensitive to the selected vegetation characteristics and water-table depth. Scenarios considering runoff infiltration into, and preferential flow through sinkholes resulted in higher and faster recharge rates. A comparison of modelled and field-based recharge estimates indicated that: (1) the model simulated plausible recharge rates, (2) only the models with preferential flow correctly reproduced the timing of recharge, and (3) preferential flow is probably redistributed in the unsaturated zone rather than passing to the water table directly. Because different but equally plausible conceptual models produce widely varying recharge rates, field-based recharge estimates are essential to constrain the modelling results.     

美国Sand Hills地区地下水数值模拟及水量平衡分析

利用地下水数值模型MODFLOW和非饱和带水平衡模型对处于半干旱半湿润沙丘地区(Sand Hills)地下水位进行了模拟,并分析了含水层补排水量,河流与地下水补排关系,以及区域水平衡过程。揭示了独特沙丘地形和土壤特性对地下水补排量的影响。模拟结果表明,入渗率大、非饱和带厚的沙丘有利于降水入渗补给,减少了地下水蒸散发损失。加上下覆含水层具有良好的地下水储水空间,是该地区储存丰富的地下水量,以维持河流稳定流量,供给众多湖泊和湿地的原因。该研究对我国地下水资源评价和生态环境脆弱地区水资源保护具有指导意义。     

淄博市大武水源地地下水水位总体预测

本文采用水均衡法对淄博市大武水源地地下水位变动进行研究，并对降水量与人渗补给量的关系以及平均给水度随地下水位的变化关系进行了分析。通过历史资料的验证，表明本文所建立的水均衡模型对于地下水位变动的预测是可靠的，而且简单可行，随后根据不同的地下水开采量等对未来地下水位的变动趋势进行了预测。     

Modeling root water uptake when calculating unsaturated flow in the vadose zone and groundwater recharge

This paper considers models of the transpiration root water uptake and principles for substantiating their parameters when calculating unsaturated flow in the vadose zone and estimating infiltration groundwater recharge. The possible use of generalized parameters for different landscape and soil conditions is substantiated based on analyzing the natural patterns of water balance formation in the vadose zone and the sensitivity of transpiration models. Recommendations are given for selecting the models of transpiration water uptake and substantiating their parameters, which enables one to minimize the errors in estimating infiltration groundwater recharge.     

Review: Current and emerging methods for catchment-scale modelling of recharge and evapotranspiration from shallow groundwater

A review is provided of the current and emerging methods for modelling catchment-scale recharge and evapotranspiration (ET) in shallow groundwater systems. With increasing availability of data, such as remotely sensed reflectance and land-surface temperature data, it is now possible to model groundwater recharge and ET with more physically realistic complexity and greater levels of confidence. The conceptual representation of recharge and ET in groundwater models is critical in areas with shallow groundwater. The depth dependence of recharge and vegetation water-use feedback requires additional calibration to fluxes as well as heads. Explicit definition of gross recharge vs. net recharge, and groundwater ET vs. unsaturated zone ET, in preparing model inputs and reporting model results is necessary to avoid double accounting in the water balance. Methods for modelling recharge and ET include (1) use of simple surface boundary conditions for groundwater flow models, (2) coupling saturated groundwater models with one-dimensional unsaturated-zone models, and (3) more complex fully-coupled surface-unsaturated-saturated conceptualisations. Model emulation provides a means for including complex model behaviours with lower computational effort. A precise ET surface input is essential for accurate model outputs, and the model conceptualisation depends on the spatial and temporal scales under investigation. Using remote sensing information for recharge and ET inputs in model calibration or in model–data fusion is an area for future research development. Improved use of uncertainty analysis to provide probability bounds for groundwater model outputs, understanding model sensitivity and parameter dependence, and guidance for further field-data acquisition are also areas for future research.     

A coupled groundwater-flow-modelling and vulnerability-mapping methodology for karstic terrain management

A coupled groundwater-flow-modelling and vulnerability-mapping methodology for the management of karst aquifers with spatial variability is developed. The methodology takes into consideration the duality of flow and recharge in karst and introduces a simple method to integrate the effect of temporal storage in the unsaturated zone. In order to investigate the applicability of the developed methodology, simulation results are validated against available field measurement data. The criteria maps from the PaPRIKa vulnerability-mapping method are used to document the groundwater flow model. The FEFLOW model is employed for the simulation of the saturated zone of Palaikastro-Chochlakies karst aquifer, in the island of Crete, Greece, for the hydrological years 2010–2012. The simulated water table reproduces typical karst characteristics, such as steep slopes and preferred drain axes, and is in good agreement with field observations. Selected calculated error indicators—Nash-Sutcliffe efficiency (NSE), root mean squared error (RMSE) and model efficiency (E′)—are within acceptable value ranges. Results indicate that different storage processes take place in different parts of the aquifer. The north-central part seems to be more sensitive to diffuse recharge, while the southern part is affected primarily by precipitation events. Sensitivity analysis is performed on the parameters of hydraulic conductivity and specific yield. The methodology is used to estimate the feasibility of artificial aquifer recharge (AAR) at the study area. Based on the developed methodology, guidelines were provided for the selection of the appropriate AAR scenario that has positive impact on the water table.     

Estimation of groundwater recharge of shallow aquifer on humid environment in Yaounde,Cameroon using hybrid water-fluctuation and hydrochemistry methods

A study of environmental chloride and groundwater balance has been carried out in order to estimate their relative value for measuring average groundwater recharge under a humid climatic environment with a relatively shallow water table. The hybrid water fluctuation method allowed the split of the hydrologic year into two seasons of recharge (wet season) and no recharge (dry season) to appraise specific yield during the dry season and, second, to estimate recharge from the water table rise during the wet season. This well elaborated and suitable method has then been used as a standard to assess the effectiveness of the chloride method under forest humid climatic environment. Effective specific yield of 0.08 was obtained for the study area. It reflects an effective basin-wide process and is insensitive to local heterogeneities in the aquifer system. The hybrid water fluctuation method gives an average recharge value of 87.14 mm/year at the basin scale, which represents 5.7% of the annual rainfall. Recharge value estimated based on the chloride method varies between 16.24 and 236.95 mm/year with an average value of 108.45 mm/year. It represents 7% of the mean annual precipitation. The discrepancy observed between recharge value estimated by the hybrid water fluctuation and the chloride mass balance methods appears to be very important, which could imply the ineffectiveness of the chloride mass balance method for this present humid environment.     

Estimating recharge at Yucca Mountain, Nevada, USA: comparison of methods

Obtaining values of net infiltration, groundwater travel time, and recharge is necessary at the Yucca Mountain site, Nevada, USA, in order to evaluate the expected performance of a potential repository as a containment system for high-level radioactive waste. However, the geologic complexities of this site, its low precipitation and net infiltration, with numerous mechanisms operating simultaneously to move water through the system, provide many challenges for the estimation of the spatial distribution of recharge. A variety of methods appropriate for arid environments has been applied, including water-balance techniques, calculations using Darcy's law in the unsaturated zone, a soil-physics method applied to neutron-hole water-content data, inverse modeling of thermal profiles in boreholes extending through the thick unsaturated zone, chloride mass balance, atmospheric radionuclides, and empirical approaches. These methods indicate that near-surface infiltration rates at Yucca Mountain are highly variable in time and space, with local (point) values ranging from zero to several hundred millimeters per year. Spatially distributed net-infiltration values average 5 mm/year, with the highest values approaching 20 mm/year near Yucca Crest. Site-scale recharge estimates range from less than 1 to about 12 mm/year. These results have been incorporated into a site-scale model that has been calibrated using these data sets that reflect infiltration processes acting on highly variable temporal and spatial scales. The modeling study predicts highly non-uniform recharge at the water table, distributed significantly differently from the non-uniform infiltration pattern at the surface. Electronic Publication     

The Guarani Aquifer System: estimation of recharge along the Uruguay–Brazil border

The cities of Rivera and Santana do Livramento are located on the outcropping area of the sandstone Guarani Aquifer on the Brazil–Uruguay border, where the aquifer is being increasingly exploited. Therefore, recharge estimates are needed to address sustainability. First, a conceptual model of the area was developed. A multilayer, heterogeneous and anisotropic groundwater-flow model was built to validate the conceptual model and to estimate recharge. A field campaign was conducted to collect water samples and monitor water levels used for model calibration. Field data revealed that there exists vertical gradients between confining basalts and underlying sandstones, suggesting basalts could indirectly recharge sandstone in fractured areas. Simulated downward flow between them was a small amount within the global water budget. Calibrated recharge rates over basalts and over outcropping sandstones were 1.3 and 8.1% of mean annual precipitation, respectively. A big portion of sandstone recharge would be drained by streams. The application of a water balance yielded a recharge of 8.5% of average annual precipitation. The numerical model and the water balance yielded similar recharge values consistent with determinations from previous authors in the area and other regions of the aquifer, providing an upper bound for recharge in this transboundary aquifer.     

疏干开采条件下晋祠岩溶水系统的水箱模型

根据晋祠岩溶水系统的水均衡要素,建立晋祠岩溶水系统水箱模型及水均衡方程。使用晋祠岩溶水系统的补给量、排泄量和岩溶水位等数据建立晋祠岩溶水贮量的亏缺量与岩溶水位间的关系,采用非线性回归分析建立泉流量与岩溶水位的关系,并利用这两个关系替换水均衡方程中的未知量,从而求解所建水箱模型。通过1981至2006年的岩溶水位、亏缺量、泉流量数据检验所建模型,并以绝对误差的绝对值及线性相关系数评价模型的模拟效果。结果表明岩溶水位、亏缺量、泉流量模拟结果的误差均小于0.5,多数误差值小于0.01,线性相关系数均大于0.7,说明所建模型能够较准确地在以年为单位的时间尺度上模拟泉流量、岩溶水位和亏缺量的变化。该模型考虑了人工补给量、管井开采量及煤矿排水量等人为因素,可用于宏观的地下水资源保护及管理。