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

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

基于VSF的非饱和-饱和带耦合过程的三维水流数值模拟

包气带作为"四水转化"的关键带,在地下水流与溶质运移模拟中起着重要作用。VSF是一款基于MODFLOW-2000开发的用于模拟非饱和-饱和带水分变化和运移过程的数值模拟程序。本文首先介绍了VSF程序的设计原理与结构,并重点介绍了VSF区别于MODFLOW-2000的用以概化非饱和带水分运移和边界条件设置的子程序包模块,包括Richards方程渗流(REF)、浸润面(SPF)、地表滞水(PND)、地表蒸发(SEV)和植物根区蒸散发(RZE)等子程序包。本文以REF、SPF、PND和SEV等子程序包为例,将VSF与MODFLOW-2000耦合应用于一个非饱和-饱和带水流运移过程的模拟算例。结果表明,VSF可与MODFLOW耦合实现饱和带水位与非饱和带饱和度的同步计算,并能刻画模型中单元格饱和状态的变化过程,其多种实用的边界条件可提供较为全面的模拟条件,是模拟实际包气带-饱和带耦合问题的有效工具。     

岩体离散裂隙网络的非饱和渗流数值分析

针对裂隙岩体的非饱和渗流问题,基于离散裂隙网络模型并结合非饱和Darcy定律、Richards方程、非饱和本构模型以及Signorini型饱和-非饱和互补溢出边界,提出了离散裂隙网络非饱和渗流问题的数学模型。采用有限单元法建立了裂隙网络非饱和渗流模型的数值求解格式和对应的迭代算法。通过与矩形坝稳定渗流、一维竖直裂隙非饱和入渗以及室内二维瞬态排水渗流的试验、数值及理论结果对比分析,验证了文中算法的有效性;根据流量等效原则,指出了裂隙网络模型应用于求解连续介质非饱和渗流问题的有效性。验证了该算法对于求解裂隙边坡降雨入渗问题的可靠性,揭示了降雨入渗过程裂隙网络流量分布的非均匀性及裂隙产状对降雨入渗流动具有重要的控制作用。     

利用算子分裂迎风均衡格式解对流为主溶质运移问题

水污染模拟问题是水流问题与溶质运移问题的耦合问题.各种常见的数值解法在以对流为主溶质运移问题的求解中都会遇到困难,如用有限单元法或有限差分法时,会产生数值弥散与过量这两类误差.引入算子分裂迎风均衡格式法求解对流为主的水污染模拟问题,较好地克服了数值弥散和数值解出现振荡问题,该格式具有良好的稳定性、单调性及守恒性特点.     

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

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

非饱和土热湿耦合传输问题的数值解——理论及一维问题的解

本文介绍了求解非饱和土壤中热量和水分耦合传输问题的一种数值方法——积分有限差分方法(IFDM)。基于菲利普-迪弗瑞斯(Philip-De Vries)多孔介质热湿耦合流动模型,采用积分有限差分方法,编制了求解热湿传输问题的计算程序 HM1,该程序可用于求解各类工程中遇到的多孔介质一维传热传湿问题。文中还给出了计算结果与实测结果相比较的实例。     

基于流网单元的污染物优势运移数值模型

垃圾填埋场滤出液、入侵海水、核废物及生产生活废水等污染物随地下水的迁移威胁人类生存。地下水渗流的随机性导致溶质运移问题更加复杂。根据流网特点,利用流线与水头等势线对求解域进行离散,基于质量守恒原理建立流网单元内溶质浓度求解的隐式有限体积差分格式。基于多孔介质孔隙流速分布规律,利用蒙特卡洛法建立流管单元随机流速场进行溶质运移过程的数值模拟,最后根据数值模拟和模型试验的结果对比,验证了数值模拟方法的准确性。基于流网单元的数值模型中沿流线方向的物质交换由对流和扩散共同作用,而流管间的物质交换只有扩散作用,因此,可在不使用弥散系数下进行污染物运移的模拟。引入随机方法确定流管内流速为研究非均匀流场中染污物的优势迁移提供了新的思路。     

加速型改进迭代法在非饱和土渗流中的应用研究

Richards方程常用于非饱和土渗流问题,并且应用广泛。在数值求解中,对Richards方程线性化,进而采用有限差分法进行数值离散以及迭代计算。其中传统的迭代法比如Jacobi迭代、Gauss-Seidel迭代法（GS）和连续超松驰迭代法（successive over-relaxation method,简称SOR）迭代收敛率较慢,尤其在离散空间步长较小以及离散时间步长较大时。因此,采用整体校正法以及多步预处理法对传统迭代法进行改进,提出一种基于整体校正法的多步预处理Gauss-Seidel迭代法（improved Gauss-Seidel iterative method with multistep preconditioner based on the integral correction method,简称ICMP(m)-GS）求解Richards方程导出的线性方程组。通过非饱和渗流算例,并与传统迭代法和解析解对比,对改进算法的收敛率和加速效果进行了验证。结果表明,提出的ICMP(m)-GS可以很大程度地改善线性方程组的病态性,相较于常规方法GS,SOR以及单一改进方法,ICMP(m)-GS具有更快的收敛率,更高的计算效率和计算精度。该方法可以为非饱和土渗流的数值模拟提供一定参考。     

含单裂隙非饱和带中轻非水相流体修复的数值模拟

针对非饱和带中油类污染物时空分布的研究,室内实验很难定量分析运移机理,野外检测成本高且破坏地层。数值模拟法作为一种应用广且成熟的方法,可以用来分析油类污染物在非饱和带中的运移规律。为了研究单井抽提及原位冲洗修复时,含单裂隙非饱和带中轻非水相流体（Light non-aqueous phase liquids,LNAPL）的时空变化规律,建立了数值模型,分析不同条件下LNAPL的修复效果及时空变化规律。模拟结果发现,LNAPL注入时优先流入裂隙,停止注入时优先流出裂隙。单井抽提修复模拟表明,抽提流量越大,修复效率越高。原位冲洗技术能有效补充地下水,防止产生新的环境问题;注水井起到“冲洗”及稀释污染物的作用,模拟最优方案修复面积达到96%,修复率达到75%,LNAPL饱和度控制在约0.05;对比分析发现,注水井布设在污染物范围的上边界时修复效果最好,能有效“冲洗”污染物并携带至抽提井中抽出地表。该研究为受轻油污染的土壤及地下水修复提供了科学的理论依据及有效的评估方法。     

非饱和土中水平入渗方程显式求解

在求解非饱和态土中水分入渗问题时,水力函数是体积含水率或者吸力的函数,致使其控制方程呈现出强非线性的特征,进而使得其求解变得十分困难。基于水分在土体介质中流动耗时取极值路径的选择这一假定,引入时间泛函,基于变分法原理将水平入渗问题转化为泛函极值问题。通过求解Euler–Lagrange方程,结合边界条件,得到非线性瞬态水平入渗问题的显式解析解。结合Brooks-Corey型水力函数,显式地求解出该类型非饱和态土的体积含水率发展分布规律。通过计算4种不同类型土体的水平入渗规律,将求解结果与已有结果以及数值结果进行对比,验证了该方法的有效性。结果表明:体积含水率分布与位置距离和湿润峰距离比值呈幂函数关系,指数取决于土-水特征曲线的形状参数;初始条件与边界条件会对体积含水率分布造成不同程度的影响。     

Error Estimates for a Time Discretization Method for the Richards' Equation

We present a numerical analysis of a time discretization method applied to Richards' equation. Written in its saturation-based form, this nonlinear parabolic equation models water flow into unsaturated porous media. Depending on the soil parameters, the diffusion coefficient may vanish or explode, leading to degeneracy in the original parabolic equation. The numerical approach is based on an implicit Euler time discretization scheme and includes a regularization step, combined with the Kirchhoff transform. Convergence is shown by obtaining error estimates in terms of the time step and of the regularization parameter.     

Two‐scale modeling of solute dispersion in unsaturated double‐porosity media: Homogenization and experimental validation

A two‐scale modeling of solute transport in double‐porosity (DP) media under unsaturated water flow conditions is presented. The macroscopic model was developed by applying the asymptotic homogenization method. It is based on theoretical and empirical considerations dealing with the orders of magnitude of characteristic quantities involved in the process. For this purpose a physical model that mimics the behavior of DP medium was built. The resulting two‐equation model relies on a coupling exchange term between micro‐ and macro‐porosity subdomains associated with local non‐equilibrium solute concentrations. The model was numerically implemented (Comsol Multiphysics^®) to simulate the macroscopic one‐dimensional physical process taking place into the porous medium of 3D periodic microstructure. A series of dispersion experiments of NaCl solution under unsaturated steady‐state flow conditions were performed. The experimental results were used first to calibrate the dispersion coefficient of the model, and second to validate it through two other independent experiments. The excellent agreement between the numerical simulations and the measurements of the time evolution of the non‐symmetrical breakthrough curves provides a proof of predictive capacity of the developed model. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling Two-Phase Incompressible Flow in Porous Media Using Mixed Hybrid and Discontinuous Finite Elements

In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors.     

Inexact Newton methods and the method of lines for solving Richards' equation in two space dimensions

Richards' equation (RE) is often used to model flow in unsaturated porous media. This model captures physical effects, such as sharp fronts in fluid pressures and saturations, which are present in more complex models of multiphase flow. The numerical solution of RE is difficult not only because of these physical effects but also because of the mathematical problems that arise in dealing with the nonlinearities. The method of lines has been shown to be very effective for solving RE in one space dimension. When solving RE in two space dimensions, direct methods for solving the linearized problem for the Newton step are impractical. In this work, we show how the method of lines and Newton-iterative methods, which solve linear equations with iterative methods, can be applied to RE in two space dimensions. We present theoretical results on convergence and use that theory to design an adaptive method for computation of the linear tolerance. Numerical results show the method to be effective and robust compared with an existing approach.     

New class of solutions for water infiltration problems in unsaturated soils

This paper presents the results of approximate analytical solutions to Richards' equation, which governs the problem of unsaturated flow in porous media. The existing methods generally fall within the category of numerical and analytical methods, often having many restrictions for practical situations. In the present study, two approximate analytical methods known as the differential transform method (DTM) and homotopy perturbation method (HPM) were employed to find analytical solutions to Richards' equation. The methods were found to be robust in finding solutions practically identical to those from the existing analytical and numerical methods. Two representative examples were considered in order to evaluate the accuracy of the solutions obtained by DTM and HPM, revealing high level of accuracy in both cases.     

Steady infiltration from buried point source into heterogeneous cross‐anisotropic unsaturated soil

The paper presents the analytical solution for the steady‐state infiltration from a buried point source into two types of heterogeneous cross‐anisotropic unsaturated half‐spaces. In the first case, the heterogeneity of the soil is modelled by an exponential relationship between the hydraulic conductivity and the soil depth. In the second case, the heterogeneous soil is represented by a multilayered half‐space where each layer is homogeneous. The hydraulic conductivity varies exponentially with moisture potential and this leads to the linearization of the Richards equation governing unsaturated flow. The analytical solution is obtained by using the Hankel integral transform. For the multilayered case, the combination of a special forward and backward transfer matrix techniques makes the numerical evaluation of the solution very accurate and efficient. The correctness of both formulations is validated by comparison with alternative solutions for two different cases. The results from typical cases are presented to illustrate the influence on the flow field of the cross‐anisotropic hydraulic conductivity, the soil heterogeneity and the depth of the source. Copyright © 2004 John Wiley & Sons, Ltd.     

基于田间试验的土壤水分运动模型选择

基于不同形式Richards方程可建立不同适用范围和计算精度的数值模型,针对具体情况下如何选择合适模型的问题,以武汉大学农田水利试验场田间入渗试验为例,选用6种模型(Picard-h模型、Picard-θ模型、Picard-mix模型、Ross模型、动力波模型和水均衡模型),运用贝叶斯模型平均(BMA)方法进行了模型选择的计算;针对BMA方法无法考虑模型计算效率的缺点,进一步提出了可同时考虑模型计算精度与计算效率的改进BMA方法。计算结果表明,在本田间尺度问题中,Ross模型排序最高,说明其兼具高精度与高效率,改进BMA方法可增加高计算效率模型被选中的概率,使模型选择更加全面合理。     

非饱和渗流模拟中非均匀空间网格的改进方法

Richards方程在非饱和渗流模拟及其他相关领域应用广泛。在数值求解过程中,可以采用有限差分方法进行数值离散并迭代求解,为了获得较可靠的数值解,常规的均匀网格空间步长往往是较小的。在一些不利数值条件下,如入渗于干燥土壤,迭代计算费时甚至精度也不能得到很好改善。因此,文章提出Chebyshev空间网格改进方法,结合有限差分方法对Richards方程进行数值离散以获得线性方程组,并通过经典的Picard迭代方法进行迭代求解线性方程组以得到Richards方程的数值解。通过均质土和分层土2个不利情况下的非饱和渗流算例,又结合模型解析解和软件Hydrus-1D,对比研究了改进网格方法与均匀网格方法获得数值解的精度。结果表明,提出的Chebyshev网格方法相较于传统的均匀网格,可以在较少的节点数下获得较高的数值精度,又具有较小的计算开销,有较好的应用前景。     

非饱和渗流问题的自适应欠松弛变量变换方法

在使用有限元方法求解非饱和土渗流问题时,土-水特征曲线和渗透率函数的强烈非线性经常会造成计算中出现迭代不收敛、计算误差大等问题。基于变量变换的思想,结合时间步长自适应技术提出了一种求解非饱和渗流问题的新方法--欠松弛RFT变换方法（ATUR1）。ATUR1方法通过变量变换,大大降低了Richards方程中未知数在空间和时间上的非线性程度,从而改善这种非线性所带来的计算收敛困难和精度差等问题。欠松弛技术的引入减少了迭代过程中的振荡现象,进一步提高了非线性迭代计算的效率。时间步长自适应技术则有效地控制整个计算过程的误差。数值算例结果说明,ATUR1可以有效地提高计算效率和精度,是一种准确有效的计算方法。     

海河流域不同下垫面土壤水分动态模拟研究

针对海河流域不同的下垫面类型,选取密云(果园林地)、大兴(城郊农田)、馆陶(平原农田)3个观测站,建立垂直方向上以含水率θ为因变量、含根系吸水项的非饱和土壤水分运动数值计算模型。该模型以一维Richards方程为基础(以下简称RE模型),采用实测的降水和蒸散数据作为模型的上边界条件,运用全隐式有限差分法,分别对不同生长期内的土壤水分进行数值模拟,得到时间序列的土壤水分廓线,并分别采用成熟软件HYDRUS-1D的模拟结果和各观测站实测土壤水分对RE模型进行交叉验证和直接验证。结果表明RE模型能够很好地模拟海河流域不同下垫面土壤水分动态变化过程,3个站模拟结果与实测土壤水分数据的均方根误差(RMSE)分别为0.03127,0.0359和0.0409 cm3/cm3。与HYDRUS-1D软件模拟结果(其与观测值的RMSE分别为0.03759,0.0647和0.0467 cm3/cm3)相比,RE模型模拟的土壤水分具有更高的精度,也显示出RE模型的可靠性。探讨3个站土壤水分的时空变异规律及其影响因子并以大兴站为例,通过优化RE模型参数,探讨犁底层对土壤水分模拟结果的影响,进一步改善RE模型的模拟精度。