区域地下水溶质运移随机理论的研究与进展

在总结近年来国内外区域地下水溶质运动研究的基本理论、方法和部分成果的基础上,论述了溶质在大区域运动的主要影响因素为区域介质的空间变异性。首先总结了野外条件下饱和介质和非饱和介质土壤渗透性能的空间变异性结果,由于野外渗透介质严重的空间变异性,研究溶质在野外条件下的运动采用了随机理论方法。基于Lagrange方法和Euler方法,研究结果表明,在渗透系数为对数正态二阶平稳及一阶扰动近似条件下,平均浓度满足对流－弥散方程,方程中宏观弥散度决定于介质渗透性能的统计特征,总结了一系列宏观弥散系数的表达形式,在此基础上,指出了需要进一步研究的问题。     

地下水观测井井筒效应的多场耦合数值模拟

利用3种不同水流运移方程分别模拟井管附近不同区域的水流运动, 基于流量守恒原理实现不同流态区域边界的耦合, 建立了有代表性的观测井-含水层系统场景; 利用建立的耦合模型模拟了观测井-含水层系统中水头的分布, 基于模型模拟数据分析了观测井井筒存在对含水层局部水头分布及地下水水质采样和环境监测结果的影响; 还分析了地下水三维水流强度、观测井井径以及含水层介质参数等对井筒效应的影响规律: 井筒效应在粘土等渗透系数和比单位贮水系数相对较小的含水层介质中更为明显, 其影响随着三维水流强度及观测井井径的增加而增大; 进行了上述参数的敏感性分析, 指出对于同一参数其在不同区间的敏感性比例不同, 对于不同参数观测井井径的敏感性比例最大, 因此在地下水环境监测的工程实践中减小观测井井径是相对快速且有效提高监测和采样精度的方法.      

考虑参数空间变异性的非饱和土坡可靠度分析

在考虑多个土体参数空间变异性的基础上,提出了基于拉丁超立方抽样的非饱和土坡稳定可靠度分析的非侵入式随机有限元法。利用Hermite随机多项式展开拟合边坡安全系数与输入参数间的隐式函数关系,采用拉丁超立方抽样技术产生输入参数样本点,通过Karhunen-Loève展开方法离散土体渗透系数、有效黏聚力和内摩擦角随机场,并编写了计算程序NISFEM-KL-LHS。研究了该方法在稳定渗流条件下非饱和土坡可靠度分析中的应用。结果表明：非侵入式随机有限元法为考虑多个土体参数空间变异性的非饱和土坡可靠度问题提供了一种有效的分析工具。土体渗透系数空间变异性和坡面降雨强度对边坡地下水位和最危险滑动面位置均有明显的影响。当降雨强度与饱和渗透系数的比值大于0.01时,边坡失效概率急剧增加。当土体参数变异性或者参数间负相关性较大时,忽略土体参数空间变异性会明显高估边坡失效概率。     

渗透系数空间变异性研究

水文地质参数的空间变异性是随机理论研究的基础,而渗透系数是最为重要的水文地质参数。国外有关渗透参数空间变异性的研究工作已开展很多,但渗透系数究竟服从什么分布目前尚无确切答案。利用Borden含水层试验数据,对渗透系数的空间变异性进行探讨,结果表明若处理方法得当,渗透系数应服从对数正态分布。同时,还对今后野外开展含水层渗透系数空间变异性试验研究提出几点建议。     

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

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

地下水位在非饱和水流数据同化中的应用

为理解地下水位观测信息在非饱和水流数据同化中的数据价值,建立了基于地下水位动态观测信息的一维饱和-非饱和水流集合卡尔曼滤波,通过虚拟数值实验检验了地下水位观测信息在非饱和水力参数估计和水分校正中的潜在价值。研究结果表明:在以地下水位为唯一观测数据时,同时更新参数和水头比仅更新水头能更好地校正土壤剖面的水头分布;当多层单个水力参数未知时,地下水位观测可以为参数估计提供有效信息;当多层多个参数未知时,地下水位与多层多个参数之间的复杂关系导致观测信息难以估计出最优的(唯一的)参数值;地下水位可作为辅助信息,与含水量观测等信息联合运用改善参数估计和含水量预测精度。     

渗透系数的空间变异性对污染物运移的影响研究

随机水文地质学方法,较传统的确定性方法而言,是解决非均质含水层中水流和溶质运移问题的一种更为合理的手段。据以往研究,假设渗透系数场遵循对数正态分布,利用直接傅立叶变换方法来生成渗透系数随机场。应用基于随机理论的蒙特卡罗方法,来研究渗透系数的空间变异性对污染物运移结果的影响。实例研究表明,污染物在含水层中运移过程中污染羽的展布范围(二阶矩)随着渗透系数空间变异方差的增大而扩大,而污染羽在空间上的质心位置(一阶矩)基本不受方差的影响,仅取决于渗透系数随机场的均值大小。另外还分析了污染羽在各点的浓度变化方差和变异系数分别随渗透系数变异方差的变化状况。     

非均质土壤中二维非饱和土壤水分运动的随机分析

以一水平向具有空间变异而垂直向为均质的二维非饱和流动区域中的均匀入渗问题为例,应用Monte Carlo随机模拟方法分析了土壤水分变量的随机统计特性及其一、二阶矩的时空分布规律。在随机分析过程中,将所研究的流动区域土壤水力特性的空间变异性以随空间位置变化的标定参数δ(x)表示,并将标定参数δ(x)视为一维随机空间函数的实现,应用随机生成模型来生成参数δ(x)的随机样本。通过随机模拟分析得到:垂直方向上负压水头方差与平均负压水头近似呈一线性关系;随入渗时间的延长,不同深度处表征负压水头空间变异结构的自相关函数趋于一稳定结构。     

含水层渗透系数空间变异性对地下水数值模拟的影响

地下水数值模拟是目前定量研究地下水水量和水质的重要手段。使用基于随机理论的MonteCarlo方法来进行地下水数值模拟。这种方法能较好地考虑水文地质参数的空间变异性。主要将MonteCarlo方法和确定性模型模拟方法的模拟结果在渗透系数场、水头场、速度场和浓度场等方面进行了比较。结果表明:在模拟三维非均质含水层中的溶质运移问题时,充分考虑了含水层渗透系数空间变异性的MonteCarlo法比确定性方法更为有效,模拟精度提高了很多,且对模拟误差及误差来源有合理的数学解释。     

相间相互作用原理与土壤水动力学基本方程

进一步阐明了多相介质力学分析的相间相互作用原理.应用相间相互作用原理重新推导了饱和土壤和非饱和土壤水分运动的控制方程,使两者具有统一的表达形式和理论基础.说明了达西定律的物理意义,在此基础上给出了非饱和土壤导水系数与饱和土壤渗透系数之间关系的表达式,该表达式在饱和条件下退化为饱和土的渗透系数.引用大连地区亚粘土和硅微粉在非饱和稳态渗流条件下的渗透试验结果验证了导水系数和饱和土的渗透系数之间的关系.     

Stochastic Modeling of Variably Saturated Transient Flow in Fractal Porous Media

This work presents the application of a Monte Carlo simulation method to perform an statistical analysis of transient variably saturated flow in an hypothetical random porous media. For each realization of the stochastic soil parameters entering as coefficients in Richards' flow equation, the pressure head and the flow field are computed using a mixed finite element procedure for the spatial discretization combined with a backward Euler and a modified Picard iteration in time. The hybridization of the mixed method provides a novel way for evaluating hydraulic conductivity on interelement boundaries. The proposed methodology can handle both large variability and fractal structure in the hydraulic parameters. The saturated conductivity K _s and the shape parameter _vg in the van Genuchten model are treated as stochastic fractal functions known as fractional Brownian motion (fBm) or fractional Gaussian noise (fGn). The statistical moments of the pressure head, water content, and flow components are obtained by averaging realizations of the fractal parameters in Monte Carlo fashion. A numerical example showing the application of the proposed methodology to characterize groundwater flow in highly heterogeneous soils is presented.     

Consolidation in spatially random unsaturated soils based on coupled flow‐deformation simulation

This paper integrates random field simulation of soil spatial variability with numerical modeling of coupled flow and deformation to investigate consolidation in spatially random unsaturated soil. The spatial variability of soil properties is simulated using the covariance matrix decomposition method. The random soil properties are imported into an interactive multiphysics software COMSOL to solve the governing partial differential equations. The effects of the spatial variability of Young's modulus and saturated permeability together with unsaturated hydraulic parameters on the dissipation of excess pore water pressure and settlement are investigated using an example of consolidation in a saturated‐unsaturated soil column because of loading. It is found that the surface settlement and the pore water pressure profile during the process of consolidation are significantly affected by the spatially varying Young's modulus. The mean value of the settlement of the spatially random soil is more than 100% greater than that of the deterministic case, and the surface settlement is subject to large uncertainty, which implies that consolidation settlement is difficult to predict accurately based on the conventional deterministic approach. The uncertainty of the settlement increases with the scale of fluctuation because of the averaging effect of spatial variability. The effects of spatial variability of saturated permeability k_sat and air entry parameters are much less significant than that of elastic modulus. The spatial variability of air entry value parameters affects the uncertainties of settlement and excess pore pressure mostly in the unsaturated zone. Copyright © 2016 John Wiley & Sons, Ltd.     

Stochastic finite element modelling of water flow in variably saturated heterogeneous soils

Water flow is greatly influenced by the characteristics of the domain through which the process occurs. It is generally accepted that earth materials have extreme variations from point to point in space. Consequently, this heterogeneity results in high variation in hydraulic properties of soil. In order to develop an accurate predictive model for transport processes in soil, the effects of this variability should be considered. In this study a two‐dimensional stochastic finite element flow model was developed for simulation of water flow through unsaturated soils. In this model, the stochastic partial differential governing equation of water flow, obtained from implementation of the perturbation‐spectral stochastic method on classical Richard's equation, was solved using a finite element method in the space domain and a finite difference scheme in the time domain. The effective hydrological parameters embedded in the mathematical model depend on time derivatives of capillary tension head; this makes possible to consider the hysteresis due to large‐scale variability of soil hydrological properties. The model is also capable of simulating infiltration and evaporation events and rapid change in the land surface boundary condition from one type event to another, based on a scheme used in the model for implementation of land surface boundary condition. The model was validated with the data obtained from a layered lysimeter test. The model was also used to simulate water flow under a long irrigation furrow. The results obtained with this model show better agreement with experimental measurements in comparison with a deterministic model. The possible reason for this agreement is that in the developed model, the influence of the variability of the properties of soil and effects of parameter hysteresis on water flow and water content redistribution are considered. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluation of probabilistic flow in two unsaturated soils

A variably saturated flow model is coupled to a first-order reliability algorithm to simulate unsaturated flow in two soils. The unsaturated soil properties are considered as uncertain variables with means, standard deviations, and marginal probability distributions. Thus, each simulation constitutes an unsaturated probability flow event. Sensitivities of the uncertain variables are estimated for each event. The unsaturated hydraulic properties of a fine-textured soil and a coarse-textured soil are used. The properties are based on the van Genuchten model. The flow domain has a recharge surface, a seepage boundary along the bottom, and a no-flow boundary along the sides. The uncertain variables are saturated water content, residual water content, van Genuchten model parameters alpha (α) and n, and saturated hydraulic conductivity. The objective is to evaluate the significance of each uncertain variable to the probabilistic flow. Under wet conditions, saturated water content and residual water content are the most significant uncertain variables in the sand. For dry conditions in the sand, however, the van Genuchten model parameters α and n are the most significant. Model parameter n and saturated hydraulic conductivity are the most significant for the wet clay loam. Saturated water content is most significant for the dry clay loam. Electronic Publication     

Three-dimensional finite element modelling for consolidation due to groundwater withdrawal in a desaturating anisotropic aquifer system

A poroelastic numerical model is presented to evaluate three-dimensional consolidation due to groundwater withdrawal from desaturating anisotropic porous media. This numerical model is developed based on the fully coupled governing equations for groundwater flow in deforming variably saturated porous media and the Galerkin finite element method. Two different cases of unsaturated aquifers are simulated for the purpose of comparison: a cross-anisotropic soil aquifer, and a corresponding isotropic soil aquifer composed of a geometrically averaged equivalent material. The numerical simulation results show that the anisotropy has a significant effect on the shapes of three-dimensional hydraulic head distribution and displacement vector fields. Such an effect of anisotropy is caused by the uneven partitioning of the hydraulic pumping stress between the vertical and horizontal directions in both groundwater flow field and solid skeleton deformation field. Copyright © 1999 John Wiley & Sons, Ltd.     

非饱和土壤水分运动参数空间变异性研究进展与展望

非饱和土壤水分运动参数的空间特征是科学认识大尺度土壤水分动态变化的基础和先决条件。在参考国内外大量文献的基础上,对非饱和土壤水分运动参数的研究进展从模型研究到空间尺度研究进行了分析与评述,指出机理模型研究是未来土壤水分运动参数空间变异性研究的方向,而尺度研究目前研究的精度不够;在对土壤水分运动参数的空间变异性研究各类方法分析和归纳后,指出土壤水分运动参数变异性研究的方法包括直接测定方法和间接估计方法两种,并对各种方法的优缺点及适用条件进行了总结;最后指出土壤水分运动参数空间变异性研究存在的问题主要集中在土壤水分运动参数的空间变异性结果的标准化问题和不同尺度土壤水分运动参数的相互转化等问题。     

随机单裂隙饱和/非饱和渗透系数研究

裂隙易在降雨作用下诱发滑坡等灾害,裂隙的饱和/非饱和渗透特性是研究此类问题的关键。利用精密数控机床制作随机粗糙裂隙面,并研制了一套仪器进行此随机粗糙裂隙的渗流试验,得到了裂隙的饱和渗透系数,然后通过间接方法预测此裂隙的非饱和渗透系数。研究发现,当裂隙平均开度为0.4 mm时,其饱和渗透系数为0.1 m/s。通过立方定律得到的水力等效隙宽为0.35 mm,小于其平均隙宽。同时裂隙的渗透系数与平均隙宽的平方成正比,这与立方定律的趋势相一致。研究得到了不同隙宽裂隙的非饱和渗透系数函数,当基质吸力小于进气值时,渗透系数为一常数,即为饱和渗透系数;当基质吸力大于进气值时,裂隙板的渗透系数急剧减小。当裂隙板的基质吸力达到其残余含水率对应的吸力值时,裂隙板的渗透系数基本稳定。在此情况下,基质吸力的继续增加对渗透系数的影响非常小,很难使渗透系数减小。