饱和-非饱和文献题录

ANN-SFLA based parameter estimation method for an unsaturated-saturated simulation model

【作者】：Das, M.; Bhattacharjya, R. K.; Kartha, S. A.
【来源】：MODELING EARTH SYSTEMS AND ENVIRONMENT, 2024, 10(1), 751-765
【摘要】：对饱和和非饱和带的地下水含水层进行数值模拟需要了解控制水流的水力参数。

然而，这些参数可能不容易获得，并且需要进行估计。

可以通过使用逆优化模型来估计参数，其中该模型最小化了观测到的和模拟的水头之间的误差函数。

由于参数估计是一个非凸问题，多个解满足所施加的约束，从而导致解的非唯一性。

另一方面，由于数值流动模型中的非线性，当与优化模型耦合时，模拟需要很高的计算时间。

本文提出了一种利用元启发式Shuffled Frog Leaping算法（SFLA）估计非饱和和饱和流量参数的新技术。此外，在模拟中，人工神经网络（ANN）被独特地结合在一起，以减少预测水头的计算时间。ANN-SFLA模型成功地估计了假设的三维地下水含水层模拟模型的非饱和和饱和参数。所提出的模型在计算时间和性能预测方面的高效性反映了其有效性。此外，使用方差分解技术进行全局灵敏度分析，以确定每个流量参数的相对重要性。
【关键词】：Non-convex; Hydraulic parameters; Constraints; Groundwater

【标题】：Application of Newton iteration to the numerical solution of variable-density groundwater flow in saturated-unsaturated coastal zones

【作者】：Yao, M.; Yu, S.; Li, H.
【来源】：FRONTIERS IN MARINE SCIENCE, 2023, 10(1127036),
【摘要】：在解决海岸带变饱和和变密度地下水流动问题时，高度非线性的水盐耦合方程组可能值得更多的关注。经典的Picard格式计算速度慢、精度低，难以满足用户的实际需求。在这里，我们开发了一种新的基于牛顿格式的沿海地下水流动问题的数值解，并通过分析海水入侵的情况，比较了不同数值方法的优缺点。仿真结果表明，变密度效应显著延长了模型的计算时间，但与Picard格式相比，牛顿格式仍具有计算速度快、收敛性好的优点，尤其是在潮汐波动高频大振幅、含水层坡度陡、网格粗的情况下。这项研究揭示了牛顿格式在解决海岸变饱和和密度地下水流动问题方面的计算特征，为海岸地下水流动问题的数值求解提供了新的思路和见解。
【关键词】：Newton scheme; Picard scheme; numerical solution; variable density flow; saturated-unsaturated coastal zones

【标题】：A Richards? equation-based model for wave-resolving simulation of variably-saturated beach groundwater flow dynamics
【作者】：Clement, J.; Sous, D.; Bouchette, F. (...)
【来源】：JOURNAL OF HYDROLOGY, 2023, 619(129344),
【摘要】：

本研究引入了一个基于理查兹方程的模型来描述可变饱和海滩地下水流量。

表面波的传播是通过相位分辨非静力SWASH程序计算的。SWASH数据用于在海滩表面建立合适的动态边界条件，以迫使Richards方程。后者采用加权不连续伽辽金法和自适应网格细化相结合的方法求解。

通过与一个瞬态地下水位补给问题的实验室实验的比较，验证了该模型的有效性。

然后，考虑BARDEX II原型规模实验来评估海滩地下水动力学的模型能力。

对三个不同泻湖水位的障壁滩进行了研究。无波动条件下的稳态结果显示出极好的一致性。

瞬态波模拟根据压头、饱和度、地下水位位置和地下水速度进行评估，用于时间平均、斜波分解和频谱分析。

通过与实验数据以及基于达西方程的模型的比较，结果为海滩地下水建模带来了有趣的见解。

首次通过修改沉积物相对重量来评估地下水对床沙动力学的影响。

【关键词】：Swash zone; Barrier beach; Numerical modelling; Unsaturated porous media; Seepage; Adaptive discontinuous Galerkin method

【标题】：Three-dimensional tunnel face stability considering the steady-state seepage in saturated and unsaturated regions with changing water levels
【作者】：Hou, C.; Zhang, Z.; Yang, X.
【来源】：COMPUTERS AND GEOTECHNICS, 2022, 146(104741),
【摘要】：

对于在水位变化的非饱和土中开挖的隧道，由于开挖，地下水位以上可能存在稳态非饱和流，地下水位以下可能存在朝向隧道表面的饱和流。

水位的变化会影响水下饱和流的不利影响以及非饱和流对剪切强度的贡献。

为了找出非饱和和饱和流效应的区域和大小的规律性，在极限分析的运动学方法与数值模拟相结合的情况下，建立了一个分析框架。

采用基于吸力应力的方程来表示非饱和流对剪切强度的贡献。

对地下水位以下饱和流的发展进行了数值模拟，并在三维离散旋转破坏机制上提取和插值得到的孔隙压力分布。

建立了包含非饱和流表观内聚力和饱和流影响的孔隙压力的功率平衡方程，以优化临界面压力。

通过与以往文献的比较和流体力学分析的数值模拟，验证了所提出的分析框架。通过参数分析，研究了地下水位、水流类型、非饱和特性和水力滞后对掌子面稳定性的影响。结果表明，当非饱和带足够大时，只需考虑非饱和流的影响以及各种非饱和特性引起的差异。该工作可为地下水位变化的非饱和土隧道开挖初期支护设计提供良好的参考。
【关键词】：Face stability; Kinematic approach; Unsaturated flow; Saturated flow; Suction stress

Generalized mass-conservative finite volume framework for unified saturated-unsaturated subsurface flow

【作者】：Dey, S.; Dhar, A.
【来源】：JOURNAL OF HYDROLOGY, 2022, 605(127309),
【摘要】：

由于理查兹方程在非饱和流动区的强非线性，多孔介质中饱和-非饱和流动的数值模拟非常具有挑战性。

稳健有效的数值模型对于准确模拟区域尺度地下系统的流动动力学至关重要。

在本文中，我们为基于压头和混合形式的Richards方程提出了广义质量守恒有限体积求解器，称为subsurfaceFlowFOAM，该求解器是使用OpenFOAM（R）框架开发的。

由于土壤-水特征关系的固有非线性和水文条件随时间的急剧变化，Richards方程的求解经常遇到数值不稳定性和收敛问题。

在这项工作中，我们提出了一种稳定和自适应的时间步长算法，该算法根据解决方案的要求调整时间步长，以确保平滑收敛。

我们通过域分解方法进一步实现了并行计算技术，以加快区域尺度建模中涉及的大规模模拟的计算过程。次表面FlowFOAM的数值结果与文献中提出的各种累积质量平衡误差不足的次表面流动问题的数据非常一致。求解器的求解精度和计算效率确保了其在解决多维区域尺度地下流问题中的潜在适用性。
【关键词】：Subsurface flow; Richards Equation; OpenFOAM (R); Stabilized and adaptive time-stepping algorithm; Parallel computation

【标题】：A saturated-unsaturated coupling model for groundwater flowing into seepage wells: a modeling study for groundwater development in river basins
【作者】：Ke, X.; Wang, W.; Xu, X. (...)
【来源】：ENVIRONMENTAL EARTH SCIENCES, 2021, 80(71121),

【摘要】：

渗水井作为一种主要且受欢迎的地下水开采结构，通常用于将河水输送到含水层中以获取水资源。

它可以帮助改善供需之间的不平衡，特别是在缺水地区。

由于抽水造成的大量抽水可能导致河流与地下水断开连接，形成不饱和带，严重影响渗水井的效率。

然而，目前大多数非管井提取结构模型只考虑饱和流动，没有考虑非饱和条件。

为了解决这一局限性，利用井管和含水层之间的交换流速作为耦合点，建立了饱和-非饱和耦合模型。

此外，还利用物理模拟试验数据对模型进行了评价。统计结果表明，对于具有较强泄流能力的河流，该模型可以很好地估计水位下降和抽水率，均方根偏差分别为0.0114m和0.0079Ls（-1）；对于泄漏能力较弱的河流，分别为0.0129 m和0.0099 L s（-1）。还讨论了河流与含水层断开的临界水位下降以及非饱和带的变化。本研究为合理降深的渗水井设计提供了重要信息，同时能够预测潜在的水量，同时有助于保护地下水环境。
【关键词】：Saturated-unsaturated coupling model; Groundwater; Unsaturated zone; Drawdown; Physical simulation test

【标题】：Efficient mass conservative numerical model for solving variably saturated groundwater flow
【作者】：Aharmouch, A.; Amaziane, B.
【来源】：JOURNAL OF HYDROLOGY, 2021, 603(126976B),

【摘要】：

提出了一种多维质量守恒数值方法，特别适用于有限的计算资源，用于求解瞬态变饱和地下水流动问题。

Richards方程在空间上采用有限元法离散，在时间上采用隐式Euler格式离散，其中使用质量守恒和质量集中技术来保持数值模拟的稳定。

此外，所涉及的刚度矩阵和质量矩阵以一种确保较少计算工作量的方式进行近似。

为了确认此代码的准确性和效率，我们使用一个、两个和三维问题的基准测试对其进行了验证。

该模型也应用于一个实际的现场案例问题，清楚地表明了它的优越性。该代码对每个问题都取得了可靠的结果。

【关键词】：Groundwater flow; Finite element; Modeling; Richards equation; Variably saturated flow

【标题】：Space-time polyharmonic radial polynomial basis functions for modeling saturated and unsaturated flows
【作者】：Ku, C.; Hong, L.; Liu, C. (...)
【来源】：ENGINEERING WITH COMPUTERS, 2021,
【摘要】：

在本文中，我们提出了一种新的无网格方法，该方法涉及使用时空多谐径向多项式基函数来模拟多孔介质中的饱和和非饱和流动。

在本研究中，使用无网格配置方法在时空域中发展了时空多谐径向多项式基函数。

该域包含三组配置点，即内部点、源点和边界点，用于控制方程的空间和时间离散化。

由于初始和边界数据是可访问的时空边界，因此地下水流动问题的解是通过求解时空域中的逆边值问题来近似的，而不使用传统的时间推进格式。

研究了饱和和非饱和流动问题，以证明所提出方法的稳健性。将使用所提出的方法获得的结果与使用传统的多谐样条径向基函数获得的结果进行了比较。所提出的时空多谐径向多项式基函数获得了高精度的解。此外，在求解饱和和非饱和流问题时，所提出的函数的精度和稳定性都高于传统的时间推进格式。
【关键词】：Space-time; Meshless; Polyharmonic; Radial basis function; Saturated; Unsaturated

【标题】：Simulating contaminant transport in unsaturated and saturated groundwater zones
【作者】：Sarma, R.; Singh, S. K.
【来源】：WATER ENVIRONMENT RESEARCH, 2021, 93(9), 1496-1509
【摘要】：随着人口的增加，对清洁和充足水的需求正在迅速增长。这种对水日益增长的需求需要对水量和水质进行测量。仿真建模在过去二十年中变得越来越流行，主要是因为它们的预测能力。本文综述了地下水资源管理中的模拟建模方法，重点介绍了用于模拟污染物通过含水层系统传输的模型。最近的研究论文整合了非饱和带和饱和带的模型，也进行了研究和描述。综合模型需要评估地下水区域之间的复杂相互更准确的建模方法。鼓励在区域一级应用此类模型，以做出最佳的地下水管理决策。从业者观点在作用以及水和溶质在其中的运动。因此，与独立开发的饱和和非饱和带模型相比，集成模型提供了过去的几十年里，地下水系统中污染物传输的建模已经得到了巨大的应用。存在许多独立模拟非饱和和饱和区域中的流动和溶质传输的模型。最近，重点放在开发需要较少输入和运行时间的高级耦合建模方法上。
【关键词】：contaminant transport modeling; coupled models; groundwater modeling; saturated zone; unsaturated zone

【标题】：High-resolution shock-capturing numerical simulations of three-phase immiscible fluids from the unsaturated to the saturated zone
【作者】：Feo, A.; Celico, F.
【来源】：SCIENTIFIC REPORTS, 2021, 11(52121),
【摘要】：非饱和和饱和多孔含水层中不混溶污染物流体流动的数值模拟对于正确管理地下水资源在许多科学领域具有重要意义。我们提出了一个高分辨率的数值模型，模拟多孔含水层中非饱和和饱和带中的三相不混溶流体流动。我们对每个相使用耦合守恒质量方程，并研究多相流体流动的动力学，作为不同相的饱和度、毛细管压力、渗透率和孔隙度、初始条件和边界条件的函数。为了处理由偏微分方程的非线性引起的尖锋，并准确地传播流体组分的尖锋。我们使用高分辨率的冲击捕获方法来处理由于毛细管压力和渗透率引起的不连续性，这些不连续性取决于三个不同相的饱和度。该问题数值求解的主要方法是基于离散（在空间中）变量的（完全）显式演化。由于显式方法要求时间步长足够小，因此这种条件非常有限，特别是对于长时间积分。随着当今多核计算机计算速度和容量的提高，可以使用高性能计算来详细模拟污染物的命运流。
【关键词】：

【标题】：An advanced discrete fracture model for variably saturated flow in fractured porous media
【作者】：Koohbor, B.; Fahs, M.; Hoteit, H. (...)
【来源】：ADVANCES IN WATER RESOURCES, 2020, 140(103602),
【摘要】：

用离散裂缝矩阵（DFM）模型对裂隙多孔介质中的变饱和流（VSF）进行精确建模是一个具有计算挑战性的问题。

DFM模型在大空间和时间尺度的实际现场研究中对VSF的适用性通常是有限的，这不仅是因为它需要详细的裂缝特征，而且还因为它涉及过多的计算工作。

我们发展了一种有效的数值格式来求解离散裂隙多孔介质中的Richards方程。该方案将空间离散化的混合混合混合有限元法与时间积分的线法相结合。裂缝被建模为二维多孔基质中的低维界面（1D）。我们为裂缝开发了一种新的质量块集（ML）技术，以消除解决方案中的非物理振荡和收敛问题，这显著提高了效率，实现了更大的现场应用。针对实验室规模的地下水位补给问题，通过商业模拟器验证了所提出的新方案。针对裂隙干土中水渗透的一个具有挑战性的问题，检验了所开发方案的计算效率，并与标准数值技术进行了比较。我们证明了ML技术对于提高鲁棒性和效率至关重要，这优于我们测试的常用方法。然后，在一项关于气候变化对黎巴嫩阿萨尔岩溶含水层/泉水系统地下水资源影响的研究中证明了我们方法的适用性。模拟，包括气候变化情景下的补给预测，进行了大约80年，直到2099年。这项研究证明了我们提出的方案在处理涉及大时间和空间尺度、高可变补给的实际现场案例时的适用性。我们的结果表明，地下水位对裂缝的存在很敏感，忽视裂缝会导致高估可用地下水量。所提出的数值方法适用于DFM模型，可以扩展到不同的二维和三维有限元框架。
【关键词】：Variably saturated flow; Fractured Porous media; Discrete fracture-matrix model; Mass lumping; Mixed hybrid finite element method; Method of Lines; climate change

【标题】：Simulation of variably-saturated flow in fractured porous media using composite element method
【作者】：Hou, X.; Chen, S.
【来源】：ROCK AND SOIL MECHANICS, 2020, 41(4), 1437-1446
【摘要】：

采用复合单元法建立了裂隙多孔介质中变饱和流动的数值模拟模型。

该模型具有以下功能。

裂缝不需要离散成特定的单元，而是根据其几何位置插入多孔基质单元中，形成复合单元。

在复合单元内，分别建立了裂缝流动和多孔基质流动的计算方程，这些方程可以通过裂缝-基质界面连接起来，并集成到复合单元方程中。

复合单元方程与传统的有限元方程格式相同，因此传统的有限单元方程的求解技术可以应用于复合单元方程。利用欠松弛迭代、质量矩阵集块和自适应时间步进技术，开发了裂隙多孔介质中变饱和流动的计算程序。通过对一维干土入渗和复杂裂隙含水层渗流的模拟，验证了该模型的合理性和适用性。模拟结果为进一步了解非饱和裂隙含水层的地下水流动特征提供了理论依据。
【关键词】：fractured porous media; variably-saturated flow; composite element method; fracture-matrix interaction

【标题】：Analytical Solutions for Water Infiltration into Unsaturated-Semi-Saturated Soils Under Different Water Content Distributions on the Top Boundary
【作者】：Sanayei, H. R. Z.; Talebbeydokhti, N.; Rakhshandehroo, G. R.
【来源】：IRANIAN JOURNAL OF SCIENCE AND TECHNOLOGY-TRANSACTIONS OF CIVIL ENGINEERING, 2019, 43(4), 747-760
【摘要】：地表水向非饱和饱和土壤渗透的建模和模拟在水文循环过程研究中具有重要意义。非饱和饱和土中的瞬态流体流动一般用Richards方程来描述。该方程是一个非线性抛物型偏微分方程，存在许多数值解和有限解析解。本文选择了含水量和导水率的指数模型，使非线性方程变为线性方程，然后使用变量分离和傅立叶级数展开技术推导了二维水渗透到部分和/或部分顶部边界上含水量不同分布的非饱和半饱和土中的解析解。在矩形土域中研究了二维水渗透，其中高含水量区被视为半饱和区，残余含水量区在初始时被视为非饱和区。本文共调查了三个案例。在情况1中，考虑了顶部边界部分的含水量均匀分布，而侧边缘没有流动边界。同时，假设非饱和-半饱和域的水平地下水位为初始条件。在情况2中，研究了倾斜的地下水位初始条件和顶部边界部分的含水量均匀分布。最后，在情况3中，获得了含水量在顶部边界和倾斜地下水位初始条件下的两个部分上的斜坡分布的解析解。为了说明推导方程的使用，从数值解中获得含水量值，并将其与两种情况下的分析解中的含水量值进行比较，显示出小于2%的误差。这些分析解可作为数值方法验证和效率评估的基准。
【关键词】：Analytical solution; Water infiltration; Unsaturated-semi-saturated soil; Two dimension; Richards' equation; Inclined water table

【标题】：An improved method for the calculation of unsaturated-saturated water flow by coupling the FEM and FDM
【作者】：Gao, Y.; Pu, S.; Zheng, C. (...)
【来源】：SCIENTIFIC REPORTS, 2019, 9(14995),
【摘要】：非饱和土壤和饱和地下水含水层中水运动的数值模拟对于水资源管理模拟非常重要。开发有效的非饱和流和饱和流耦合数值算法一直是水文建模中的一个长期挑战，尤其是在区域尺度建模中。本研究将有限元法（FEM）和有限差分法（FDM）相结合，提出了一种新的求解理查兹方程的方法，即FE-FDM，用于模拟非饱和饱和水流。控制方程在水平方向上采用有限元离散，在垂直方向上采用FDM离散。该方法结合了有限元法在域离散化（尤其是在复杂计算域）方面的优势，以及FDM在建模简单高效方面的优势。通过三个测试案例验证了新方法的有效性，表明FE-FDM解是准确的，适用于区域尺度的问题。在测试案例中，将FE-FDM解与实验数据和用常用软件包（包括FEFLOW和COMSOL）计算的数值结果进行了比较。本研究验证了FE-FDM适用于非饱和饱和带的水流模拟。
【关键词】：

【标题】：A comprehensive quasi-3-D model for regional-scale unsaturated-saturated water flow
【作者】：Mao, W.; Zhu, Y.; Dai, H. (...)
【来源】：HYDROLOGY AND EARTH SYSTEM SCIENCES, 2019, 23(8), 3481-3502
【摘要】：对于区域尺度上非饱和饱和流的计算高效建模，考虑一维（1-D）土壤水流和三维地下水流的准三维（3-D）格式是一种替代方法。然而，由于土壤水分建模所需的高度非线性和密集的输入数据以及耦合方案的可靠性，区域尺度问题在实践中仍然具有挑战性。本研究开发了一种新的准三维模型，该模型与具有MOD-FLOW三维水动力学模型的UBMOD一维土壤水分平衡模型相耦合。提出了一种以垂直净补给和非饱和带深度为交换信息的迭代方案的新实现方法。开发了一个建模框架来组织土壤-水模型和地下水模型的耦合方案，并处理前后信息。通过使用两项已发表的研究评估了耦合模型的优势和劣势。对比结果表明，该耦合模型在计算精度和质量平衡误差方面都是令人满意的。讨论了空间和时间离散化以及应力周期对模型精度的影响。此外，在一项真实世界的研究中，耦合模型被用于评估地下水补给。测量的地下水位和土壤含水量用于校准模型参数，2年示踪实验的地下水补给数据用于评估补给估算。现场应用进一步证明了该模型的实用性。所开发的模型和建模框架为在区域尺度上评估非饱和饱和流系统提供了一个方便而灵活的工具。
【关键词】：

【标题】：On the uncertainty of initial condition and initialization approaches in variably saturated flow modeling
【作者】：Yu, D.; Yang, J.; Shi, L. (...)
【来源】：HYDROLOGY AND EARTH SYSTEM SCIENCES, 2019, 23(7), 2897-2914
【摘要】：土壤水分运动直接影响环境、农业和水文。土壤水分运动的模拟需要准确确定模型参数以及初始和边界条件。然而，在模拟开始时很难获得准确的初始土壤湿度或基质电位剖面，因此有必要从任意初始条件运行模拟模型，直到初始条件（UIC）的不确定性减小，这通常被称为“升温”。在本文中，我们比较了两种常用的UIC量化方法（一种是基于在多个水文年中递归运行单个模拟，另一种是以实现各种初始条件的蒙特卡洛模拟为基础），并确定了不同土壤质地、气象条件和土壤剖面长度所需的预热时间t（wu）（通过预热模型消除UIC所需的最小时间）。然后，我们在两个数据同化框架（即集成卡尔曼滤波器和迭代集成平滑器）内分析了不同初始条件对参数估计的影响，并评估了几种现有的模型初始化方法，这些方法使用可用数据来检索初始土壤水分剖面。我们的结果表明，蒙特卡罗模拟和多年的递归模拟都可以证明UIC的时间行为，并建议使用一个通用阈值来确定t（wu）。此外，还定量地确定了变饱和流建模的t（wu）与模型设置（土壤质地、气象条件和土壤剖面长度）之间的关系。此外，为了获得更好的参数和状态估计性能，我们在同化数据之前提出了一个预热期。
【关键词】：

【标题】：Development and evaluation of a variably saturated flow model in the global E3SM Land Model (ELM) version 1.0
【作者】：Bisht, G.; Riley, W. J.; Hammond, G. E. (...)
【来源】：GEOSCIENTIFIC MODEL DEVELOPMENT, 2018, 11(10), 4085-4102
【摘要】：改进近地表土壤水分和地下水水文的全球尺度模型表示对于准确模拟陆地过程和预测气候变化对水资源的影响非常重要。大多数现有的地表模型，包括我们在这里修改的默认E3SM陆地模型（ELMv0），通常使用不同的公式来计算渗流带和潜水带中的水传输。Clark等人（2015）将变饱和Richards方程流动模型确定为改进土壤水分和浅层地下水动力学耦合模拟的重要能力。在这项工作中，我们在ELMv1中开发了变饱和流模型（VSFM），以统一非饱和区和饱和区土壤水文过程的处理。VSFM在三个基准问题上进行了测试，并根据观测结果和现有的基准模型（PFLOTRAN）对结果进行了评估。ELMv1 VSFM的地下排水参数f（d）经过校准，以匹配观测约束和空间显式的全球地下水位深度（WTD）乘积。79%的全局1.9度x 2.5度网格单元获得了最佳的空间显式f（d）值，而其余21%的全局网格单元预测的WTD比观测约束的估计更深。与使用默认f（d）值的预测进行比较表明，校准显著改善了预测，主要是通过允许更深的WTD。使用国际土地模型基准包（ILAMB）进行的模型评估表明，WTD预测的改进不会降低任何其他指标的模型技能。我们评估了VSFM模型的计算性能，发现该模型比具有最佳处理器布局的默认ELMv0贵约30%。VSFM的模块化软件设计不仅提供了为一系列问题设置配置模型的灵活性，还允许独立于ELM代码构建模型，从而能够针对其他模型直接测试模型的物理特性。
【关键词】：

【标题】：Analysis of three-dimensional unsaturated-saturated flow induced by localized recharge in unconfined aquifers
【作者】：Chang, C.; Huang, C.; Yeh, H.
【来源】：HYDROLOGY AND EARTH SYSTEM SCIENCES, 2018, 22(7), 3951-3963
【摘要】：在地下水补给过程中，地表水通常通过上覆的非饱和带进入含水层。到目前为止，由于局部补给，对耦合的非饱和-饱和流动模型的解析解的发展很少给予关注。本文建立了一个数学模型来描述地表局部补给的无侧限含水层中的三维瞬态非饱和饱和流。该模型包含非饱和流的Richards方程、饱和地层的流动方程和描述非饱和土壤特性的Gardner本构模型。这两个流动方程通过水头和流量在地下水位的连续性条件进行耦合。采用拉普拉斯变换和傅立叶余弦变换的方法，导出了耦合流动模型的半解析解。进行敏感性分析，以探索水头对每个含水层参数变化的响应。提出了一种定量工具来评估补给效率，表示从补给到含水层的水的百分比。我们发现，如果满足与非饱和土壤性质和初始含水层厚度相关的两个标准，则非饱和流对饱和水头的影响可以忽略不计。根据本解决方案预测的水头分布与有限差分模拟的水头分布非常匹配。本解决方案的预测也与加利福尼亚州弗雷斯诺县人工补给池现场实验的观测数据非常一致。
【关键词】：

【标题】：Transient Water Flow and Nitrate Movement Simulation in Partially Saturated Zone
【作者】：Jahangeer; Gupta, P. K.; Yadav, B. K.
【来源】：JOURNAL OF IRRIGATION AND DRAINAGE ENGINEERING, 2017, 143(0401704812),
【摘要】：渗流带水流对土壤生态和地下水开发具有重要意义。对于可持续灌溉管理、水土保持措施、污染风险评估、地下水补给和应用修复技术，也需要全面了解通过非饱和带的水流和溶质输送。因此，在本研究中，进行了实际和模拟实验，以评估排水通量、补给流入对水分流动和硝酸盐在地下运动的作用。采用非线性Richards方程模拟了非饱和带的水流，该方程的参数与土壤含水量和压头有关。在本研究中，评估了以下情况下土壤水分和硝酸盐通过可变饱和区的传输：（1）研究区域的排水通量，（2）研究区域中的补给通量，（3）非饱和土柱中的瞬态渗透，以及（4）二维（2D）砂箱设置中的硝酸盐运动。在研究领域中，使用单一孔隙模型来模拟土壤水力参数。结果表明，在底部边界施加排水通量后，地下水位位置随时间显著下降。在恒定的补给流量下，研究区域内的地下水位位置缓慢增加。当施加恒定的渗透通量时，土壤水分锋向地下水位移动。此外，硝酸盐的运动以平流通量为主，并在垂直方向上受到补给通量的显著影响。（C） 2017年美国土木工程师学会。
【关键词】：Unsaturated zone; Transient water flow; HYDRUS 2/3D; Water recharge; Tank experiments

【标题】：Saturated-unsaturated groundwater modeling using 3D Richards equation with a coordinate transform of nonorthogonal grids
【作者】：Deng, B.; Wang, J.
【来源】：APPLIED MATHEMATICAL MODELLING, 2017, 50, 39-52
【摘要】：复杂地形下的饱和-非饱和流通常使用Richards方程求解。由于编码的简单性，Richards方程的离散化通常采用有限差分法或有限体积法。复杂的次表面边界几何形状导致曲线网格系统中的非正交网格，这给离散化和网格生成带来了困难。本文开发了一种垂直坐标变换，使计算域在垂直方向上规则。结果，曲线表面的网格可以成功地转换为计算网格，该计算网格允许以高效的计算和更简单的编码来求解理查兹方程。将笛卡尔坐标系中的各向异性理查兹方程转换为任意坐标系中，并进一步表示为适用于正交坐标系的形式。将广义第三边界条件转换为适用于正交坐标系的形式。采用有限体积法求解正交坐标系下的理查兹方程。通过四个实例对现有正交坐标系进行了验证。正交坐标系的计算结果与Ansys Fluent在笛卡尔坐标系中求解的地下流情况的结果非常一致。对于山坡的耦合情况，计算结果与实验数据吻合良好。目前V型集水区和平板情况的结果与HydroGeoSphere和PAWS的结果非常一致。该算法可以改进复杂域水文模型中理查兹方程的网格生成。（C） 2017爱思唯尔股份有限公司保留所有权利。
【关键词】：Richards equation; Orthogonal coordinate transform; Finite volume method; Groundwater modeling; Hydrological model

【标题】：Regional Quasi-Three-Dimensional Unsaturated-Saturated Water Flow Model Based on a Vertical-Horizontal Splitting Concept
【作者】：Zhu, Y.; Shi, L.; Wu, J. (...)
【来源】：WATER, 2016, 8(1955),
【摘要】：由于三维（3-D）非饱和饱和水流方程的高度非线性，对于大规模应用来说，使用完全3-D的数值模型在计算上是昂贵的。本文基于垂直/水平分裂（VHS）概念，建立了一种新的非饱和饱和水流模型，将三维非饱和饱和Richards方程分解为二维水平方程和一维垂直方程。导出了三角棱镜单元中平均水头梯度的水平面，以将三维方程分解为二维方程。然后用水平衡法计算由二维方程表示的平均水头梯度水平面内的横向流量。采用有限差分法对一维垂向方程进行离散。通过将这两个方程耦合到一个具有单个矩阵的统一非线性系统中，同时求解。通过将建模结果与Hydrus1D、SWMS2D和FEFLOW的建模结果进行比较，使用三个综合案例对所开发的模型代码进行评估。我们进一步将该模型应用于区域尺度建模，模拟地下水位波动，以评估模型在复杂条件下的适用性。所提出的建模方法被发现相对于测量是准确的。
【关键词】：unsaturated-saturated model; vertical-horizontal splitting concept; water balance method; finite difference method; regional scale modeling

【标题】：A new quasi-3D unsaturated-saturated hydrogeologic model of the Plateau de Saclay (France)
【作者】：Renard, F.; Tognelli, A.
【来源】：JOURNAL OF HYDROLOGY, 2016, 535, 495-508
【摘要】：为占地74公里的萨克雷高原（法国巴黎西南20公里）开发了一个新的水文地质模型。它基于枫丹白露金沙含水层的2D模型，结合非饱和带的1D模型，考虑到地下水渗透场的空间变异性和非饱和带厚度。本文基于瞬态流动模拟，重点估算了三重参数（渗透率、渗透性和有效孔隙度）。首先，通过反演不同渗透值的平均水头数据来获得渗透率。然后，通过使用瞬态流模拟和拟合1970-2014年期间几个测压计的水头测量值，确定渗透和有效孔隙度。渗透量估计在125至150 mm/yr之间，平均渗透率在2.10（-5）至6.10（-5）m/s之间，有效孔隙度在20%至30%之间。此外，对不饱和带的作用进行了量化，在示踪剂运输的情况下，诱导延迟估计约为20年。（C） 2016爱思唯尔B.V.保留所有权利。
【关键词】：Hydrogeology; Saclay; Permeability; Porosity; Vadose zone; Infiltration

【标题】：Simulating water flow in variably saturated soils: a comparison of a 3D model with approximation-based formulations
【作者】：Hopp, L.; Fatichi, S.; Ivanov, V. Y.
【来源】：HYDROLOGY RESEARCH, 2016, 47(2), 274-290
【摘要】：在水文模型中，经常使用理查兹方程来描述可变饱和流量，无论是在全三维（3D）实现中，还是使用基于垂直流量的1D理查兹方程和其他二维流量近似的准3D框架。然而，尚不清楚这些近似在哪种配置或在哪种边界条件下可以产生足够的估计。在这项研究中，采用准3D方法的两种配方以全3D模型（HYDRUS-3D）为基准。公式为：使用Dupuit Forchheimer假设的实时综合盆地模拟器+VEGetation Generator for Interactive Evolution（tRIBS+VEGGIE）模型和实现运动学方法的Tethys&Chloris（T&C）模型。在基于事件的山坡模拟中，研究了域坡度、山坡大小、事件大小和初始水分条件对模拟径流和土壤水分动力学的影响。Dupuit Forchheimer假设（tRIBS VEGGIE）与HYDRUS-3D解决方案的偏差仅适用于初始干燥土壤的模拟。使用运动学方法（T&C）会导致与3D解决方案的偏差，主要针对小山坡区域和缓坡角。这尤其适用于地下径流生产和地表径流生产之间的划分，T&C偏向后者。对于所研究的所有其他情况，更简单的公式提供了3D模型的合理近似值。
【关键词】：hillslope hydrology; hydrologic modeling; parsimonious approximation; surface-subsurface hydrology

【标题】：PHT3D-UZF: A Reactive Transport Model for Variably-Saturated Porous Media
【作者】：Wu, M. Z.; Post, V. E. A.; Salmon, S. U. (...)
【来源】：GROUNDWATER, 2016, 54(1), 23-34
【摘要】：开发了基于MODFLOW/MT3DMS的反应传输模型PHT3D的修改版本，以将当前的反应传输能力扩展到地下系统的可变饱和组分，并结合气体物质的扩散反应传输。该代码被称为PHT3D-UZF，包含了由MODFLOW的非饱和区流量（UZF1）包计算的通量项。采用了与UZF-MT3DMS中使用的方法类似的体积平均方法。PHT3D-UZF内基于PHREEQC的化学过程计算与UZF1的分析求解方法相结合，可以进行全面的反应迁移研究（即生物地球化学转化），这些研究共同涉及饱和和不饱和带过程。UZF1适用于区域规模的应用，仅模拟不饱和带内的向下通量。通过与现有数值模型的模拟结果进行比较，对该模型进行了验证。对几个基准问题进行了比较，这些问题涵盖了一系列重要的水文和反应输送过程。定义了一个二维模拟场景，以说明在沙质酸性硫酸盐土壤环境中脱水后的地球化学演化。其他潜在应用包括模拟可变饱和系统中的生物地球化学过程，跟踪农业污染物、营养物质、天然和异生有机化合物以及药物等微污染物的运输和去向，以及同位素模式的演变。
【关键词】：

【标题】：Modeling Variably Saturated Multispecies Reactive Groundwater Solute Transport with MODFLOW-UZF and RT3D
【作者】：Bailey, R. T.; Morway, E. D.; Niswonger, R. G. (...)
【来源】：GROUND WATER, 2013, 51(5), 752-761
【摘要】：建立了一个能够模拟变饱和多孔介质中多相反应性溶质迁移的数值模型。该模型由反应传输模型RT3D（三维反应传输）的修改版本组成，该模型与不饱和区流（UZF1）包和MODFLOW相关联。该模型被称为UZF-RT3D，根据已公布的分析基准以及其他已公布的污染物传输模型进行测试，包括HYDRUS-1D、VS2DT和SUTRA，以及CATHY和TRAN3D的耦合流动和传输建模系统。探讨了一维、二维和三维变饱和系统的比较。虽然包括了几个测试案例来验证UZF-RT3D中可变饱和传输的正确实现，但也包括了其他案例来证明该代码在模型运行时间和处理可变饱和地下系统中多种相互作用物种的反应动力学方面的有用性。由于UZF1依赖于忽略Richards方程中扩散项的非饱和流的运动波近似，因此UZF-RT3D可用于UZF1公式合理的大型含水层系统，即毛细管压力梯度可忽略，土壤参数可被视为均质。减少的模型运行时间和包括特定地点化学物质和化学反应的能力使UZF-RT3D成为一个有吸引力的模型，用于有效模拟可变饱和大规模地下系统中的多谱反应传输。
【关键词】：

【标题】：PARALLEL MODELING OF THREE-DIMENSIONAL VARIABLY SATURATED GROUND WATER FLOWS WITH UNSTRUCTURED MESH USING OPEN SOURCE FINITE VOLUME PLATFORM OPENFOAM
【作者】：Liu, X.
【来源】：ENGINEERING APPLICATIONS OF COMPUTATIONAL FLUID MECHANICS, 2013, 7(2), 223-238
【摘要】：我们报道了基于多孔介质中可变饱和流的多物理平台OpenFOAM的三维开源模型suGWFoam的开发。利用有限体积法和改进的Picard迭代法在非结构网格上对非线性Richards方程进行空间离散，该方法本质上保证了质量守恒。实现了Richards方程不同形式之间的切换算法，并提供了适用于不同模拟场景的三种收敛准则选择。线性化方程得到了一个对称矩阵系统，该系统可以用以下Krylov子空间线性求解器之一求解：预条件共轭梯度（PCG）、不完全Cholesky预条件共轭斜率（ICCG）和几何凝聚代数多重网格（GAMG）。运行时可选择各种预处理器。通过域分解的并行计算在大规模模拟中显示出良好的加速效果。与其他可用的地下水代码相比，该模型的优点在于OpenFOAM的体系结构。该模型建立在平台的模块化和层次化设计之上，其顶层代码不到几百行，但具有三维非结构化网格、并行计算以及多种空间和时间离散化方案的完整功能。这将水资源界的研究人员从繁琐而漫长的编码过程中解放出来。使用OpenFOAM中现有的流求解器（如通道流和波浪），将地下水流与地表水流体动力学耦合起来也非常容易。该模型功能齐全，验证表明与文献中的数据在广泛的应用中非常一致。它既可用于基础研究，也可用于现实世界中的应用。
【关键词】：ground water; aquifer; unsaturated flow; three-dimensional models; finite volume method; parallel computing

【标题】：Modeling coupled unsaturated and saturated nitrate distribution of the aquifer Westliches Leibnitzer Feld, Austria
【作者】：Klammler, G.; Kupfersberger, H.; Rock, G. (...)
【来源】：ENVIRONMENTAL EARTH SCIENCES, 2013, 69(2), 663-678
【摘要】：奥地利Westliches Leibnitzer Feld含水层是为超过100000名居民提供区域和超区域饮用水的重要资源，但该地区也提供了良好的农业条件。这种双重用途涉及冲突（例如，氮浸出造成的非点源地下水污染），必须协调这些冲突才能实现可持续共存。在含水层规模上，数值模型是模拟地下水数量和质量行为的最先进工具，也是实施地下水保护措施的决策支持系统。虽然完全和迭代耦合的模拟模型考虑了饱和和非饱和区之间的反馈，但沙质土壤条件和根区下的地下水深度允许在调查区域使用非饱和水流和硝酸盐传输模型SIMWASER/STOTRASIM与FEFLOW的单向顺序耦合。考虑到地表水体、农业、住宅和森林地区向地下水的水和氮输入是分开的，第一次模拟结果与观测到的地下水位相匹配，但总体上低估了硝酸盐浓度。因此，模拟了假设地表氮输入较高的多个场景，以与测量的硝酸盐浓度收敛。初步结果表明，农田对地下水的氮输入主要受农田贡献的影响。
【关键词】：Modeling/modelling; Coupled model; Nitrogen/nitrate; Unsaturated/vadose zone; Saturated zone; Aquifer scale

【标题】：Modeling Variably Saturated Subsurface Solute Transport with MODFLOW-UZF and MT3DMS
【作者】：Morway, E. D.; Niswonger, R. G.; Langevin, C. D. (...)
【来源】：GROUND WATER, 2013, 51(2), 237-251
【摘要】：通过结合为MODFLOW开发的非饱和带流（UZF1）软件包，对MT3DMS地下水溶质运移模型进行了修改，以模拟非饱和带中的溶质运移。修改后的MT3DMS代码使用体积平均法，将基于拉格朗日量的UZF1流体通量和储存变化映射到固定网格上。该链接模型被称为UZF-MT3DMS，根据分析求解的已发布基准以及其他已发布代码进行测试，最常见的是美国地质调查局的变饱和二维流动和传输模型。一组测试用例的结果表明，修改后的代码准确地模拟了非饱和带中的溶质平流、扩散和反应。还研究了二维和三维模拟，以确保正确模拟非饱和-饱和区的相互作用。因为UZF1解是解析的，所以可以进行大规模的流动和输运研究，而不需要理查兹方程数值解所需的计算和数据负担。结果表明，使用UZF-MT3DMS可以显著节省仿真运行时间，这是在参数估计和不确定性分析期间需要数百或数千次模型运行时的一个重要发展。三维可变饱和流动和输运模拟显示，UZF-MT3DMS的运行时间不到依赖Richards方程的模型所需时间的十分之一。考虑到其准确性和效率，以及MODFLOW和MT3DMS的广泛使用，在这个熟悉的建模框架中增加的非饱和区运输能力将有利于广泛的用户船舶。
【关键词】：

Assessment of a 3D Variably Saturated Subsurface Model

【作者】：Caviedes-Voullieme, D.; Murillo, J.; Garcia-Navarro, P.
【会议】：PROCEEDINGS OF THE 35TH IAHR WORLD CONGRESS, VOLS III AND IV 2013 Int Assoc Hydro Environm Engn & Res, Chengdu, PEOPLES R CHINA
【主办单位】：
【摘要】：一些环境、灌溉和工程应用对三维地下水流的精确模拟特别感兴趣。理查兹方程是最常用于单相、等温、可变饱和流动的物理模型。这是一个需要数值求解的非线性偏微分方程，其计算成本可能很高，特别是在需要大量单元的3D域中。在这项工作中，提出了一个三维理查兹方程求解的数值模型。该方案建立在有限体积框架内，并在时间上隐式集成。隐式格式在理论上是无条件稳定的，尽管精度可能会受到大时间步长的影响。该模型处理地表的自然干燥或洪水条件，以及土壤中的饱和或非饱和条件。数值格式是质量守恒的和精确的。它也很容易退化为2D和1D模型。针对1D和3D地下水分析解决方案以及实验案例和合成测试案例进行基准测试。通过这些测试，对模型的准确性进行了测试，并对效率进行了研究。使用六面体和四面体单元研究了对网格分辨率和结构的敏感性。还研究了对时间步长的敏感性。最后，评估了模型的计算成本，并评估了不同建模选择的效率。
【关键词】：3D Richards equation; Finite Volume; Mesh topology; Unsteady flow in porous media

1介绍

理查兹方程(理查兹，1931)描述了可变饱和多孔介质中的流动。它已被广泛解决使用有限差分(克莱门特，1994;Dogan, 2005)或有限元素(Diersch, 1999;Paniconi, 1994;Panday, 1993;Simunek, 2012;Wu, 2010)，因为方程的椭圆/抛物线性质。有限体积法(FVM)也是求解理查兹方程的一个有趣的选择，因为它可以很容易地控制质量守恒。已经研究了形式数学问题(Eymard et al . 1999,2001)，已经分析了一维中的守恒、稳定性和精度(caviedes - voulli me et al . 2012)。Manzini et al(2004)制定了二维有限体积方案，并在异质土壤和植物根系上进行了测试。

结构网格通常用于求解多维域的理查兹方程。而非结构化网格的使用要少得多。在某些情况下，为了准确地表示几何形状的表面和体积以及流动方向的可变性，非结构化网格很有意义，并且包括在本工作中，以评估具有此类网格类型的模型的性能。

本文提出了一种三维有限体积格式来求解结构化和非结构化网格中的理查兹方程。详细描述了该方案，包括一般细胞拓扑的边界条件的制定。此外，它还探讨了三维有限体积模型对解析解、一维实验测试用例和二维异质测试用例的适用性和敏感性。特别注意了仿真中的误差行为和计算时间。

【标题】：A fully coupled numerical modeling for regional unsaturated-saturated water flow
【作者】：Zhu, Y.; Shi, L.; Lin, L. (...)
【来源】：JOURNAL OF HYDROLOGY, 2012, 475, 188-203
【摘要】：开发数值高效的非饱和流和饱和流耦合算法是地下水文建模中的一个长期挑战，尤其是在区域尺度建模中。在本研究中，提出了一种新的区域尺度非饱和-饱和水流耦合数值模拟方案。根据空间分布的输入、流动特征和地形条件，将建模领域划分为地平线上的子区域。每个子区域的非饱和区由单个一维土柱表示。采用水平衡分析方法建立了地下水三维模型。通过非饱和土柱和饱和含水层之间的垂直流动，非饱和带和饱和带在空间和时间上隐含地耦合，因为非饱和区和饱和带中的水头被整合在一个矩阵方程中。通过将模拟结果与广泛使用的软件（包括Hydrorus1D、SWMS2D、FEFLOW和HydroGeoSphere）的模拟结果进行比较，验证了耦合方案，并在几个假设示例中评估了计算效率。在实际应用中，数值结果表明，该耦合模型可以在相当小的计算成本下获得令人满意的模拟结果。与现有模型相比，新的数值格式更适合于复杂域几何和交替充放电通量的区域尺度模型。然而，由于方法开发中涉及的假设，耦合方法有其固有的局限性，在非饱和带中侧向流占主导地位的情况下应谨慎使用。（C） 2012爱思唯尔B.V.保留所有权利。
【关键词】：Unsaturated-saturated flow; Large-scale water flow simulation; Model verification; Richards' equation; Water balance analysis method

【标题】：Improved predictions of saturated and unsaturated zone drawdowns in a heterogeneous unconfined aquifer via transient hydraulic tomography: Laboratory sandbox experiments
【作者】：Berg, S. J.; Illman, W. A.
【来源】：JOURNAL OF HYDROLOGY, 2012, 470, 172-183
【摘要】：对无侧限含水层抽水试验的解释在很大程度上是基于忽略含水层异质性的分析解决方案。在这项研究中，我们研究了是否可以通过包括瞬态水力层析成像（THT）中的导水率（K）和比蓄水量（S-S）信息来改进非均质无侧限含水层及其上方非饱和带的水位下降响应预测。我们还研究了这些预测是否受到通过实验室悬柱实验或现场排水曲线校准估计的非饱和流量参数的影响。为了研究这些问题，我们设计并进行了实验室沙箱实验，以表征非均质无侧限含水层的饱和和非饱和特性。具体而言，我们在完全饱和条件下进行了抽水试验，并通过将介质处理为均质或非均质来解释压降响应。然后，我们进行了另一次抽水试验，使地下水位下降，类似于无侧限含水层的抽水试验。使用可变饱和流模型进行的模拟表明：（1）饱和和非饱和带的均质参数很难预测非均质无侧限含水层的响应；（2） 通过THT获得的非均匀饱和水力参数分布显著改善了无侧限含水层饱和区的水位下降预测；和（3）考虑非饱和带参数的不均匀性，在非饱和带的预测中产生了微小的改进，但在饱和带没有。这些结果似乎支持了Mao等人的发现。（2011）非饱和带的空间变异性在无侧限含水层抽水过程中观察到的S形水位下降时间曲线的形成中起着次要作用。Crown版权所有（C）2012爱思唯尔B.V.出版。保留所有权利。
【关键词】：Unconfined aquifer; Unsaturated zone; Subsurface heterogeneity; Hydraulic tomography; Model calibration; Model validation

【标题】：Comparison between iteration schemes for three-dimensional coordinate-transformed saturated-unsaturated flow model
【作者】：An, H.; Ichikawa, Y.; Tachikawa, Y. (...)
【来源】：JOURNAL OF HYDROLOGY, 2012, 470, 212-226
【摘要】：本文比较了三维坐标变换饱和-非饱和流模型的三种不同迭代方法。皮卡德迭代法和牛顿迭代法是求解理查兹方程的常用方法。Picard方法易于实现且具有成本效益（在单独迭代的基础上）。然而，它的收敛速度比牛顿方法慢。另一方面，尽管牛顿方法收敛得更快，但与Picard方法相比，它的实现更复杂，每次迭代消耗更多的CPU资源。这两种方法在饱和-非饱和流有限元模型中的比较在以前的研究中已经得到了很好的评价。然而，在坐标变换的有限差分模型（FDM）中，两种迭代方法可能表现出不同的行为。此外，Newton-Krylov方法可能是坐标变换FDM的合适替代方法，因为它需要评估19点模板矩阵。19点模板的形成是一个相当复杂和费力的过程。相反，Newton-Krylov方法计算矩阵向量积，可以通过计算原始非线性函数的差来容易地近似矩阵向量积。在这方面，Newton-Krylov方法可能是坐标变换FDM最合适的迭代方法。然而，该方法涉及在Newton-Krylov方法中的每次Krylov迭代中进行近似的额外成本。在本文中，我们评估了三种迭代方法Picard的效率和鲁棒性。使用三维坐标变换的FDM模拟通过多孔介质的饱和-非饱和流动的牛顿和牛顿-Krylov方法。（C） 2012爱思唯尔B.V.保留所有权利。
【关键词】：Richards equation; Saturated-unsaturated flow; Finite difference model (FDM); Coordinate transformation method; Newton-Krylov method

【标题】：Source Screening Module for Contaminant Transport Analysis Through Vadose and Saturated Zones
【作者】：Bedekar, V.; Neville, C.; Tonkin, M.
【来源】：GROUND WATER, 2012, 50(6), 954-958
【摘要】：在复杂的地点，渗流带内可能存在许多潜在的污染物来源。筛选水平分析有助于确定哪些潜在来源区域应成为详细调查和分析的重点。已经开发了一个源筛选模块（SSM），以支持对渗流带废物场地对地下水质量构成的威胁的初步评估。该工具实现了分析解决方案，以模拟污染物通过非饱和区和饱和区的传输，从而预测潜在地下水受体的时变浓度。SSM在一个用户友好的、基于Microsoft Excel（TM）的界面内实现的单个模拟中集成了多个传输过程。
【关键词】：

【标题】：A Quasi-Three-Dimensional Variably Saturated Groundwater Flow Model for Climate Modeling
【作者】：Xie, Z.; Di, Z.; Luo, Z. (...)
【来源】：JOURNAL OF HYDROMETEOROLOGY, 2012, 13(1), 27-46
【摘要】：在本研究中，通过将三维土壤水和地下水流近似划分为非饱和垂直土壤水流和水平地下水流，建立了准三维变饱和地下水流模型，以模拟土壤水、地下水和植被之间的相互作用。所开发的模型包括以地下水位为移动边界的一维非饱和土壤水流模型，该模型使用基于水平域中离散网格单元形成的垂直土柱的自适应网格结构，水平域的二维地下水流动模型，以及连接域中水平网格单元的两个组件的界面模型。利用该模型进行了综合实验，测试了模型参数对河流高程、地表导水率和地表流量的敏感性，实验结果表明了该模型在不同条件下的稳健性。将该模型与全三维方案的模拟结果进行了比较，验证了其可行性和有效性。以塔里木河下游河道输水为例，对所建立的应苏段地下水位高程模拟模型进行了验证。最后，利用该模型对塔里木河流域进行了数值试验，探讨了大尺度高起伏地形下的地下水潜流问题。结果表明，该模型能较好地模拟地下水位。
【关键词】：

【标题】：Slurry wall containment performance: monitoring and modeling of unsaturated and saturated flow
【作者】：Pedretti, D.; Masetti, M.; Marangoni, T. (...)
【来源】：ENVIRONMENTAL MONITORING AND ASSESSMENT, 2012, 184(2), 607-624
【摘要】：执行了一项为期2年的具体计划，对渗流区和饱和区进行监测和测试，并结合数值分析，以评估用于控制污染区域的泥浆墙系统的整体性能。尽管存在局部物理限制（泥浆墙嵌入平均2米厚的弱透水层），但至少20年来，在米兰（意大利）市附近供应饮用水的含水层中观察到高浓度的氯化溶剂（高达110毫克/升（千分之一欧元））。监测和现场试验的结果已用于执行非饱和饱和数值模型。这些结果产生了必要的定量信息，既可用于确定该地区不同介质的水力特性，也可用于校准和验证数值模型。调查含水层较浅部分的回填材料极大地影响了封装区域的自然补给。从潮湿期到干旱期的瞬态模拟突出了横向屏障泄漏之间的比率变化，这支持了通过安全壳系统的水损失的特定场景。监测和建模的结合允许在不使用任何侵入性技术的情况下对物理约束的总体性能进行可靠的估计。
【关键词】：Slurry wall; Groundwater monitoring; Numerical modeling; Leakage; Milan (Italy)

【标题】：COUPLING OF UNSATURATED AND SATURATED SOIL ZONE MODELS TO SIMULATE THE WATER BALANCE IN MINING AREAS
【作者】：Meyer, M.; Sallwey, J.; Blankenburg, R. (...)
【会议】：12TH INTERNATIONAL MULTIDISCIPLINARY SCIENTIFIC GEOCONFERENCE, SGEM 2012, VOL. III 2012 Albena, BULGARIA
【主办单位】：
【摘要】：该项目的目的是将非饱和和饱和土壤带的模型结合起来，以一种面向过程的方式模拟水平衡，特别是矿区的地下水补给，特别是考虑到不断变化的气候条件。两个单独的模拟软件是耦合模型的基础：（1）莱比锡IBGW的PCGeofim（R），它使用有限体积法模拟三维地下水流动，特别是采矿区的地下水流动；（2）德累斯顿大学废物管理和污染场地处理研究所的PCSiWaPro（R）是一个二维模型，通过求解非饱和土壤区瞬态流动的Richard方程来模拟水平衡。必须为耦合定义交换过程和参数，以及两个模拟软件的不同尺寸和类型的离散化。使用蒸渗计数据进行评估，以验证建模的地下水补给通量与测量数据。
【关键词】：water balance; groundwater; unsaturated zone; coupling

【标题】：Coupling of unsaturated and saturated soil zone models to estimate groundwater recharge in mining areas
【作者】：Sallwey, J.; Meyer, M.; Graeber, P. W. (...)
【会议】：GROUNDWATER MODELING AND MANAGEMENT UNDER UNCERTAINTY 2012 Kuwait Inst Sci Res, KUWAIT
【主办单位】：
【摘要】：对于露天矿区的模拟和管理来说，地下水补给是影响水平衡的一个因素。以过程为导向的方式看待它将有助于考虑不断变化的气候条件。为了实现面向过程的地下水补给计算，将非饱和和饱和土壤带的模型耦合起来。两个单独的模拟工具是耦合模型的基础：PCGeofim（R）和PCSiWaPro（R），前者使用有限体积法模拟三维地下水流动，特别是矿区地下水流动；后者通过求解非饱和土壤区瞬态流动的RICHARDS方程来模拟水平衡。对于耦合交换，定义了过程和参数，实现了地下水补给计算的新程序，并适应了两个模拟软件的不同离散化。使用蒸渗计数据进行评估，以验证模拟的地下水补给通量与测量数据。
【关键词】：

【标题】：Vegetation controls on variably saturated processes between surface water and groundwater and their impact on the state of connection
【作者】：Banks, E. W.; Brunner, P.; Simmons, C. T.
【来源】：WATER RESOURCES RESEARCH, 2011, 47(W11517),
【摘要】：包气带在地表水-地下水相互作用中发挥着重要作用，并对生物地球化学、生态和次水文过程产生了强烈影响。它也是控制地表水和地下水之间连接状态的不饱和带的存在。尽管最近在水文地质变量如何影响地表水-地下水相互作用方面取得了进展，但对降水和蒸散的水文气候影响的了解有限。更具体地说，需要对植被变化可能发生的变化进行基于物理的理解。虽然植被的存在可能会导致河床下的不饱和带发育并改变连接状态，这在质量上似乎是显而易见的，但迄今为止尚未得到定量证明。此外，到目前为止，还没有探索根系灭绝深度、地形和土地清理的影响等变量的影响。在这项研究中，使用完全耦合的、基于物理的二维瞬态均质模型来模拟土地清理和植被重建对常年河流系统连接状态的影响。模拟表明，植被的存在会在河流和含水层之间形成一个不饱和区，并影响连接状态，从集水区清除根深蒂固的植被可能会对连接状态和水资源状况产生重大影响。
【关键词】：

【标题】：Variably saturated flow and multicomponent biogeochemical reactive transport modeling of a uranium bioremediation field experiment
【作者】：Yabusaki, S. B.; Fang, Y.; Williams, K. H. (...)
【来源】：JOURNAL OF CONTAMINANT HYDROLOGY, 2011, 126(3-4), 271-290
【摘要】：2008年原位铀生物修复现场实验的三维耦合可变饱和流和生物地球化学反应迁移模型用于更好地了解在脉冲乙酸盐修正、季节地下水位变化、空间可变物理（导水率、孔隙率）和地球化学（反应表面积）材料特性下控制铀行为的迁移和生物地球化工反应的相互作用。虽然2008年在科罗拉多州里夫进行的大锈乙酸盐生物刺激现场实验的模拟与之前在里夫IFRC现场进行的现场实验中确定的行为基本一致，但额外的过程和特性细节提供了一些新的见解。这项工作的一个主要结论是，当金属还原菌可获得超过硫酸盐还原菌需求的乙酸盐时，铀的生物还原最为有效。蛋白质组学分析中发现的最初少量缓慢生长的硫酸盐还原菌的加入，使生物硫化物促进的Fe（III）矿物溶解，从而增加了Fe（II）的来源。实验期间地下水位的下降显著降低了含水层的饱和厚度，并在新的不饱和包气带中产生了反应物和产物以及未缓解的铀。高渗透性砂砾石结构导致注入井附近的局部高流速，增加了醋酸盐稀释度。在降级位置，这些结构为乙酸盐输送创造了优先的流动路径，增强了TEAP反应性和辅助反应的局部区域。相反，与渗透率较低的岩相（如细粒岩相）和包气带相关的较小输送速率限制了乙酸盐的进入和反应。然而，一旦乙酸盐进入，这些相同的区域限制了随后的乙酸盐稀释，并提供了更长的停留时间，当末端电子供体和受体不受限制时，这导致了更高浓度的TEAP反应产物。最后，基于岩相的孔隙度和反应表面积变化表明会影响含水铀浓度分布，其中精细岩相的局部影响对U（VI）表面络合作用的影响最大。在自然复杂的Rifle IFRC地下系统中，对综合生物地球化学反应网络以及控制铀行为的空间和时间可变过程、特性和条件进行建模的能力需要一个地下模拟器，该模拟器可以使用大规模并行计算机的大内存和计算性能。在这种情况下，eSTOMP模拟器在128个处理器核上运行12小时，用于模拟110天的现场实验和50天的生物刺激后行为。（C） 2011爱思唯尔B.V.保留所有权利。
【关键词】：Uranium; Bioremediation; Reactive transport modeling

【标题】：A new Iterative Alternating Direction Implicit (IADI) algorithm for multi-dimensional saturated-unsaturated flow
【作者】：An, H.; Ichikawa, Y.; Tachikawa, Y. (...)
【来源】：JOURNAL OF HYDROLOGY, 2011, 408(1-2), 127-139
【摘要】：饱和-非饱和地下流的数值模拟广泛应用于科学和工程的许多分支，计算机技术的快速发展不仅使一维模拟成为可能，而且使多维模拟成为可能。然而，多维地下流模拟仍然会给计算资源带来沉重的负担，特别是对于长周期宽区域的模拟。迭代交替方向隐式（IADI）格式在计算成本和算法简单性方面具有一定的优势。然而，目前很少使用它，因为它偶尔会引起数值不稳定性和收敛困难。另一个原因是原始IADI方案无法进行三维模拟。本研究提出了一种求解饱和-非饱和流动方程的高级IADI算法：该高级格式比原始IADI格式在数值上更稳定，可用于三维模拟。通过测试模拟对所提出方案的性能进行了评估。在所有的模拟中，新方法被证明比通过改进的Picard迭代方法线性化的全隐式格式更快，同时仍然产生非常相似的结果。（C） 2011爱思唯尔版权所有。保留所有权利。
【关键词】：Richards equation; Saturated-unsaturated flow; ADI scheme; Finite difference model

【标题】：Two-way coupling of unsaturated-saturated flow by integrating the SWAT and MODFLOW models with application in an irrigation district in arid region of West China
【作者】：Luo, Y.; Sophocleous, M.
【来源】：JOURNAL OF ARID LAND, 2011, 3(3), 164-173
【摘要】：本文在SWATMOD99地表水综合模型的基础上，实现了SWAT2000和MODFLOW96模型非饱和-饱和流相互作用的双向耦合，并在内蒙古河套灌区的应用。对SWAT2000和MODFLOW模型进行了重大修订和增强，用于模拟HID的详细水文预算和耦合的非饱和和饱和相互作用以及灌溉渠水文。给出并讨论了季节性地下水向浅层地下水补给和蒸发的模拟结果，以及该地区的年度水量预算。结果表明，当地下水较浅时，非饱和-饱和相互作用双向耦合的必要性，以及综合利用地表水和地下水模型的信息对耦合相互作用进行更具物理基础的评估的可行性。
【关键词】：SWAT; MODFLOW; SWATMOD; coupled unsaturated-saturated flow; irrigation canals; groundwater recharge; groundwater evaporation

【标题】：Stochastic finite element modelling of water flow in variably saturated heterogeneous soils
【作者】：Nezhad, M. M.; Javadi, A. A.; Abbasi, F.
【来源】：INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, 2011, 35(12), 1389-1408
【摘要】：水流在很大程度上受到过程所经过的域的特性的影响。人们普遍认为，地球物质在空间中的点与点之间有着极端的变化。因此，这种不均匀性导致土壤水力特性的高度变化。为了开发土壤中迁移过程的准确预测模型，应考虑这种可变性的影响。在本研究中，建立了二维随机有限元流动模型，用于模拟非饱和土中的水流。在该模型中，采用空间域有限元法和时域有限差分格式求解了由经典Richard方程的扰动谱随机方法得到的水流随机偏微分控制方程。数学模型中嵌入的有效水文参数取决于毛细管张力水头的时间导数；这使得考虑由于土壤水文特性的大规模变化而引起的滞后成为可能。该模型还能够模拟渗透和蒸发事件以及陆地表面边界条件从一种类型的事件到另一种类型事件的快速变化，基于该模型中用于实现陆地表面边界状态的方案。该模型通过分层蒸渗仪试验获得的数据进行了验证。该模型还用于模拟长灌溉沟下的水流。与确定性模型相比，用该模型获得的结果与实验测量结果显示出更好的一致性。这种一致性的可能原因是，在所开发的模型中，考虑了土壤性质变化的影响以及参数滞后对水流和含水量再分配的影响。版权所有（C）2010 John Wiley&Sons，有限公司。
【关键词】：stochastic finite element method; heterogeneous soil; water flow; unsaturated soil

【标题】：Modeling Unsaturated-Saturated Flow and Nickel Transport in Fractured Rocks
【作者】：Masciopinto, C.; Caputo, M. C.
【来源】：VADOSE ZONE JOURNAL, 2011, 10(3), 1045-1057
【摘要】：这项研究调查了2001年被未经授权的污泥废物处理污染的地点（意大利南部的Altamura）断裂岩石中可变饱和流和Ni迁移的模型。时间推移电阻率层析成像（ERT）和渗透仪实验结果用于约束非饱和流模型中的近表面边界条件。考虑到土壤和岩石表面的不规则性，塑料环被用作渗透计，因为它的实验装置非常通用，适用于许多不同的地质条件。通过使延时ERT成为减少计算不确定性的有效工具，所提出的方法允许从层析成像切换到由模型获得的水压等值线图。模拟结果预测了镍在Altamura遗址包气带中的浓度和停留时间。这些结果被用于连续研究Ni向深层裂隙含水层的水平迁移。模拟提供了地下水中明显的Ni路径，以及距离污染源10.9km的下游井中的预期浓度。结果与所收集的采样数据一致，证实地下水被Ni污染了16个月。污染羽流始于2004年冬天，在地下水中以10 m/d的平均速度向大海移动。
【关键词】：        

【标题】：Study on unsteady saturated-unsaturated flow subjected to variable soil water characteristic curve
【作者】：Cao, P.; Tao, Y.; Shu, L.
【会议】：HYDROLOGICAL CYCLE AND WATER RESOURCES SUSTAINABILITY IN CHANGING ENVIRONMENTS 2011 Hohai Univ, Nanjing, PEOPLES R CHINA
【主办单位】：
【摘要】：基于Richards方程，提出了一个有限差分程序，用于模拟河流-含水层系统饱和-非饱和多孔介质中的瞬态非稳态和非线性渗流。为了响应地下水位波动引起的非饱和含水层结构的变化，采用可变土壤-水特征曲线来模拟地下水渗流。垂直渗流段自下而上分为饱和区、过渡区和垂直渗流区。水平流和垂直流在饱和区和过渡区同步发生，垂直渗流区只有渗流的垂直运动。空间离散化采用有限差分和线性形状函数，时间离散化采用Picard格式。自动生成响应于一般网格结构的基本信息，从而获得适应无渗表面的自动均匀网格。基于砂箱试验，讨论了自由面渗流的地下水位和径流量对土壤水分特征曲线的响应。结果表明，与固定的土壤水分特征曲线相比，可变曲线提供了更准确的模拟结果。
【关键词】：saturated-unsaturated flow; soil water characteristic curve; salt-box experiment; Richards equation; finite difference

【标题】：Three-dimensional finite difference saturated-unsaturated flow modeling with nonorthogonal grids using a coordinate transformation method
【作者】：An, H.; Ichikawa, Y.; Tachikawa, Y. (...)
【来源】：WATER RESOURCES RESEARCH, 2010, 46(W11521),
【摘要】：多孔介质中饱和-非饱和流动的研究在科学和工程的许多分支中都引起了人们的兴趣。在现有的各种数值模拟方法中，有限差分法具有离散化简单的优点。这种方法已被广泛用于模拟饱和-非饱和流动。然而，几何复杂流场的模拟需要在传统的有限差分模型中使用高分辨率网格，因为传统的有限分离散化假设正交坐标系。这使得有限差分模型在计算上不如可以处理非正交网格的其他数值模型（例如有限元模型和有限体积模型）有效。为了克服这一缺点，我们使用坐标变换方法，建立了一个模拟饱和-非饱和流的多维有限差分模型；该模型可以处理非正交网格。为了便于编码，通过坐标变换方法导出的交叉导数项被明确地评估。因此，与具有正交网格的传统有限差分模型一样，在迭代计算中使用7点模板作为隐式项。我们通过进行测试模拟来评估所提出模型的性能。然后，我们将模拟结果与密集网格解进行比较，以评估所提出模型的数值精度。为了检查所提出的模型的性能，我们将使用所提出的模式获得的模拟结果与使用（1）完全隐式考虑所有项的模式、（2）有限元模式和（3）具有高分辨率正交网格的传统有限差分模式获得的结果进行了比较。
【关键词】：

【标题】：Two-dimensional finite elements model for selenium transport in saturated and unsaturated zones
【作者】：Tayfur, G.; Tanji, K. K.; Baba, A.
【来源】：ENVIRONMENTAL MONITORING AND ASSESSMENT, 2010, 169(1-4), 509-518
【摘要】：建立了二维有限元模型，模拟了饱和和非饱和土壤中硒的二维迁移。该模型考虑了植物对水、硒酸盐、亚硒酸盐和硒蛋氨酸的吸收。它还考虑了硒酸盐、亚硒酸盐和硒代蛋氨酸的吸附和解吸、氧化和还原、挥发以及化学和生物转化。除了模拟水流、硒酸盐、亚硒酸盐和硒蛋氨酸的传输，该模型还模拟了有机硒和气态硒的传输。将所开发的模型应用于模拟两个不同的现场观测数据。对观测数据的模拟令人满意，平均绝对误差为48.5μg/l，平均相对误差为8.9%。
【关键词】：Finite elements method; Unsaturated zone; Selenium (Se); Selenate; Selenite; Selenomethionine; Organic Se; Gaseous Se; Modeling

【标题】：Assignment of groundwater level at the wastewater collecting pipe for a landfill groundwater simulation: utilization of the two-dimensional saturated-unsaturated groundwater numerical model
【作者】：Dang, T. H.; Jinno, K.; Tsutsumi, A.
【来源】：ENVIRONMENTAL EARTH SCIENCES, 2010, 60(8), 1575-1582
【摘要】：为了了解填埋作业对地下水流动行为的影响，建立了二维水平地下水模拟模型。该模型节省了计算机内存和时间，并与三维地下水流模型进行了比较。然而，最大的困难是研究现场收集管边界的分配。因此，采用二维垂直模型来计算集水管上方地下水位的变化。本文重点研究了通过应用二维垂直模型的结果来验证集水管边界分配的有效性。将二维水平模型与补给模型相结合，求解地下水流动的偏微分方程。采用有限差分法和迭代逐次过松弛法。总结了整个垃圾填埋场渗滤液收集的排水量，并与处理废水的平均排水量进行了比较。研究表明，二维垂直模型的结果是有效的，可以应用于二维水平地下水流动模拟。
【关键词】：Collecting pipe; Landfill; Groundwater flow model

【标题】：Depth-specific water sampling system for detailed sampling in the unsaturated zone, the capillary fringe and the saturated zone
【作者】：Walther, M.; Fank, J.; Reimann, T.
【来源】：GRUNDWASSER, 2010, 15(1), 19-32
【摘要】：非饱和带和饱和带之间过渡区的流动和输送过程会影响地下水的质量和性质。整个过渡区是一个物理和生物化学高度活跃的反应室，必须考虑特殊的水力特性。使用现有的测量技术，无法完全理解当前的过程，或者另一方面，需要付出巨大的努力。因此，过渡带的影响往往被忽视（Ronen和Sorek，《不饱和带——被忽视的自然组成部分》，以色列水利委员会水质司，特拉维夫，2005年）。利用实验室实验和数值模拟，对特定深度测量技术的创新概念进行了评估。结果表明，经过对实际设计的一些修改，该装置能够以高的时间和空间分辨率对非饱和区和饱和区之间的整个过渡区进行采样，包括毛细管条纹。
【关键词】：Measuring devices; Depth-specific sampling; Unsaturated zone; SUIR; Capillary fringe; Upper groundwater body; Multi-layer suction cup sampler

【标题】：MOVEMENT OF DENSE PLUMES IN VARIABLY SATURATED POROUS MEDIA: NUMERICAL MODEL AND RESULTS
【作者】：Graf, T.; Simmons, C. T.; Boufadel, M. C. (...)
【会议】：PROCEEDINGS OF THE XVIII INTERNATIONAL CONFERENCE ON COMPUTATIONAL METHODS IN WATER RESOURCES (CMWR 2010) 2010 Barcelona, SPAIN
【主办单位】：
【摘要】：卫生填埋场意外泄漏的渗滤液的总溶解固体浓度可达40000 mg/L。因此，渗滤液的密度明显高于不饱和和饱和地下的水。因此，溢出土壤或在处置场底部释放的渗滤液将通过变密度流输送通过非饱和土壤区，并最终到达饱和地下水区。为了更好地了解地下羽流的传输，Simmons等人（2002）在充满沙子的玻璃罐中，在完全饱和和可变饱和的流动条件下进行了实验室实验。在本研究中，建立了一个新的数值模型来模拟变饱和条件下稠密流体的运动。新模型基于现有的HydroGeoSphere模型，并使用Boufadel等人（1999）定义的Elder（1967）问题的非饱和流版本对多孔介质中的二维非饱和饱和变密度流动和输运进行了验证。新模型可以替代地将第一级或第二级的Oberbeck-Bossinesq（OB）假设应用于可变密度流。在模拟Simmons等人（2002）进行的实验室实验时，进行了数值模拟，重点关注OB假设的有用性。模拟表明，当应用OB2时，盐指倾向于更早地聚结，并且OB1更适合模拟变密度流。
【关键词】：Density; unsaturated; Oberbeck; Boussinesq; HydroGeoSphere

【标题】：Patterns and dynamics of river-aquifer exchange with variably-saturated flow using a fully-coupled model
【作者】：Frei, S.; Fleckenstein, J. H.; Kollet, S. J. (...)
【来源】：JOURNAL OF HYDROLOGY, 2009, 375(3-4), 383-393
【摘要】：利用基于并行物理的地表-地下模型PARFLOW研究了具有深地下水位的非均质冲积河流含水层系统中河流含水层交换的空间模式和时间动力学。将两个尺度的含水层非均质性纳入模型。根据条件地质统计指标模拟来表示冲积层的水文相结构。利用平行高斯模拟创建了水力相体内水力传导率（K）的亚尺度变异性。在蒙特卡罗框架下研究了亚尺度非均质性的影响。模拟了30天流量事件中河流含水层交换的动力学和模式。模拟结果表明，在以高导流体相为特征的优先流动区，河道和深层地下水位之间迅速形成饱和连接。当河流与低电导率的水文相相交时，河流形态正下方的浅层栖息饱和区，但深层地下水位的渗流仍然不饱和，渗流率较低。优先流区虽然只占河道的50%左右，但占总渗漏量的98%以上。地下水补给通过这些区域最有效地实现。K-sat的亚尺度变化略微增加了渗流量，但没有改变一般的渗流模式（优先流动区与上层滞流区）。总的来说，得出的结论是，典型的冲积非均质性（水文相结构）是覆盖深层地下水位的河流中河流含水层交换的重要控制因素。模拟的模式和动力学与现场观测结果以及之前使用更简单模型进行的建模研究的结果一致。冲积层的异质性导致河流含水层交换的不同模式和动态，对地下水补给和河岸带的管理（例如，通过饱和带的河道-洪泛平原连通性）有影响。（C） 2009爱思唯尔B.V.保留所有权利。
【关键词】：River-aquifer interactions; Heterogeneity; Physically-based surface-subsurface model; Variably-saturated flow; PARFLOW; Groundwater recharge

【标题】：An Efficient Lumped Mixed Hybrid Finite Element Formulation for Variably Saturated Groundwater Flow
【作者】：Belfort, B.; Ramasomanana, F.; Younes, A. (...)
【来源】：VADOSE ZONE JOURNAL, 2009, 8(2), 352-362
【摘要】：渗流带渗透的精确数值模拟仍然是一个挑战，尤其是当对非常尖锐的锋面进行建模时。在这项工作中，我们使用了混合混合有限元（MHFE）方法，该方法允许同时逼近压头和速度，并可以处理具有高度非均匀渗透率的一般不规则网格。然而，对于处理非饱和水流的许多问题，MHFE解表现出显著的非物理振荡。为了避免这种现象，我们用MHFE方法开发了一种有效的质量集块格式来求解混合形式的Richards方程。在这项工作中，比较了标准和集中MHFE公式，并确定了集中公式减少一维和二维渗透问题的非物理振荡的能力。基于M-矩阵性质的理论分析，保证了离散最大值原理，以及本研究中进行的实际测试案例，强调了使用锐角三角测量来完全消除非物理振荡的优势。事实上，与标准方法相反，集总公式满足M-矩阵性质，而对要使用的时间步长没有任何约束。
【关键词】：

【标题】：A state-dependent parameterization of saturated-unsaturated zone interaction
【作者】：Bogaart, P. W.; Teuling, A. J.; Troch, P. A.
【来源】：WATER RESOURCES RESEARCH, 2008, 44(W1142311),
【摘要】：地下水通过毛细上升作为根区水分的重要来源，其相关性越来越得到人们的认可。这在一定程度上反映在许多当前的地表方案中，这些方案越来越多地通过根区和地下水系统之间的双向相互作用通量来取代单向（即向下）排水。这种双向饱和-不饱和相互作用通量的完全物理正确实现需要使用高度非线性的Richards方程进行瞬态模拟，这是一种计算要求很高的方法。我们测试了一个经典的简单近似，该近似将根区地下水相互作用通量计算为准稳态条件下的达西方程的向下排水通量和向上毛细管上升通量的净效应。我们发现，对于湿根区和/或浅层地下水，该近似值内的误差很大，与相互作用通量本身的大小相同。我们提出了一种新的基于达西方程的通量闭合形式参数化，该参数化同时考虑了根区土壤湿度和地下水位深度。该参数化的参数值列出了11种不同的、广泛应用的土壤质地描述。所提出的方法的高数值效率使其适用于要求苛刻的应用，例如蒙特卡罗框架或高空间分辨率。
【关键词】：

【标题】：Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions
【作者】：Park, C.; Beyer, C.; Bauer, S. (...)
【来源】：ENVIRONMENTAL GEOLOGY, 2008, 55(8), 1755-1766
【摘要】：对地下系统中的传输过程进行精确有效的数值模拟是许多现场调查或修复研究的先决条件。随机游动粒子跟踪（RWPT）方法过去曾被引入，以克服模拟多孔介质中的传播过程（如平流主导的质量传输）时的数值困难。RWPT方法精度的关键是地下水速度场的数值计算精度。在本文中，使用了Yeh（Water Resur Res 7:1216-12251981）最初提出的基于全局节点的速度计算方法。该方法有三个方面的改进：（1）扩展到非结构化网格，（2）显著提高计算效率，（3）扩展到饱和（地下水）和非饱和系统（土壤-水）。利用文献中的数值基准实例对新的RWPT方法进行了测试，并将其用于饱和和非饱和地下水中污染物迁移的两个现场应用。为了评估模型的平流输运，通过比较几种已发表的粒子路径线或流线的结果来证明速度场的准确性。给定所选的测试问题，发现基于全局节点的速度估计与CK方法一样准确（Cordes和Kinzelbach in Water Res 28（11）:2903-2911992），但不如高度非均质介质中流动的混合或混合混合有限元方法准确。为了评估平流扩散输运，Hassan和Mohamed（J Hydrol 275（3-4）:242-2602003）研究了一个输运问题，并使用不同数量的粒子进行了评估。结果表明，与公布的数字相比，使用所提出的方法将给定问题所需的粒子数量减少了约两个数量级，而不会损失浓度等值线的准确性。
【关键词】：velocity approximation; particle tracking; random walk; transport in porous media; contaminant leaching from road constructions

【标题】：Generation of two-dimensionally variable saturated hydraulic conductivity fields:: Model theory, verification and computer program
【作者】：Taskinen, A.; Sirvio, H.; Bruen, M.
【来源】：COMPUTERS & GEOSCIENCES, 2008, 34(8), 876-890
【摘要】：利用具有二维相关结构的对数正态分布随机场的一组比例因子，从不同的来源汇编了空间可变K（s）场的生成方法理论。通过推导一般情况下模型参数的计算公式，进一步发展了它。该方法的用户需要定义K（s）字段的平均值、变异系数和相关长度。各向同性场和各向异性场都可以产生。通过分析由一组不同参数值执行的示例生成，验证了提供给读者的方法和计算机程序。由于生成场中设定值的具体化是令人满意的准确的，并且它们的分布与模型结构相一致，因此该方法和计算机程序可以被认为是模拟空间可变导水率场的有用和可靠的工具。然而，该方法处理随机向量的协方差结构的方式限制了要生成K（s）值的字段大小或节点数量。（C） 2007爱思唯尔有限公司保留所有权利。
【关键词】：spatial variability; random field; generation; Monte Carlo simulation; spatial correlation; log-normal distribution

【标题】：Simple models for subsurface solute transport that combine unsaturated and saturated zone pathways
【作者】：Connell, L. D.
【来源】：JOURNAL OF HYDROLOGY, 2007, 332(3-4), 361-373
【摘要】：提出了一系列半解析和解析解，用于将表面释放的污染物垂直通过非饱和区输送到地下水位，然后在地下水中。通过非饱和带的垂直迁移由平流-扩散方程的一维解表示，该方程用于计算地下水位处的溶质通量。这与拉普拉斯空间中三维平流-扩散方程的解相耦合，以表示随后的地下水横向输送。细长羽流近似用于推导饱和区三维输运的简明方程，并推导了点、尖和矩形区域源的方程形式。提出了三种不同的方法来描述不饱和区的迁移，它们在植物根区的表示方式上有所不同。复射耦合输运关系在拉普拉斯空间中，并使用de Hoog算法进行瞬态计算。建立了稳态溶质输运的解析表达式。在一系列假设模拟中，将瞬态结果与使用SUTRA的等效数值分析进行了比较。所开发的方法被证明工作良好，并为数值分析提供了一种有效的替代方法，特别是在对稳态感兴趣的情况下。（c） 2006爱思唯尔B.V.保留所有权利。
【关键词】：

【标题】：Unsaturated-zone leaching and saturated-zone mixing model in heterogeneous layers
【作者】：Lee, S. S.
【会议】：Experimental Unsaturated Soil Mechanics 2007 Bauhaus Univ Weimar, Weimar, GERMANY
【主办单位】：
【摘要】：开发了一个筛选水平模型，用于模拟污染物通过非饱和区的迁移以及随后在饱和区内的混合。该一维有限差分模型模拟了污染物的液相平流、液相和气相扩散、吸附和衰变的传输过程。该模型通过使用简单的质量平衡技术将非饱和区渗滤液与地下水混合来估计饱和区的污染物浓度。该模型可以成为对地下污染物的潜在影响进行初步评估的有用工具。该模型可以处理土柱的垂直非均质性和初始污染物浓度的不均匀性。通过将其模拟结果与分析溶液和实验室土柱实验进行比较，对其进行了验证。此外，该模型还通过三种不同土壤样本量的渥太华石英砂和480ppm盐水作为地下水污染源的实验室实验进行了验证。
【关键词】：unsaturated-zone; saturated-zone; leaching model; groundwater; modelling; GUI

【标题】：Quantifying colloid retention in partially saturated porous media
【作者】：Zevi, Y.; Dathe, A.; Gao, B. (...)
【来源】：WATER RESOURCES RESEARCH, 2006, 42(W12S0312),
【摘要】：胶体污染物复合物和胶体大小的病原体通过土壤向地下水的迁移令人担忧。与仅提供积分信号的穿透曲线相比，孔尺度胶体行为的可视化和量化将能够更好地描述和模拟单个表面的保留机制。我们测试了两种方法来量化在孔尺度图像序列中观察到的胶体运动和保留。在用静态图像进行初步测试后，检查了合成微珠悬浮液通过不饱和砂的三组图像。区域程序（在ImageJ中实现）和布尔程序（在KS400中实现）对于初始测试图像和三个图像系列中的总胶体覆盖区域产生了几乎相同的结果。由于电子噪声导致像素级亮度波动，布尔过程倾向于低估附着的胶体计数，反之则高估移动胶体计数。区域程序对附着胶体的高估误差较小。孔尺度上胶体保留的可靠量化可用于提高当前对胶体在不饱和多孔介质中传输机制的理解。例如，Langmuir等温线很好地描述了单个空气/水弯月面/固体界面处的附着计数。
【关键词】：

【标题】：Operator-splitting errors in coupled reactive transport codes for transient variably saturated flow and contaminant transport in layered soil profiles
【作者】：Jacques, D.; Simunek, J.; Mallants, D. (...)
【来源】：JOURNAL OF CONTAMINANT HYDROLOGY, 2006, 88(3-4), 197-218
【摘要】：将地下流动和传输过程与（生物）地球化学反应相结合的一种可能方法是通过算子分裂方法将两个完全独立的代码耦合，一个用于可变饱和流动和溶质传输，另一个用于平衡和动力学生物地球化学反应。本文评估了涉及瞬态大气边界条件（降水、蒸散）和层状土壤剖面的典型土壤环境问题的多组分系统算子分裂方法的准确性。最近开发的HP1程序用于求解耦合输运和化学方程。对于稳态流动条件，发现精度主要是所采用的空间离散化的函数，而在较小程度上是时间离散化的。对于瞬态流动情况，精度以复杂的方式取决于网格离散化、时间步长和主要流动条件（渗透与蒸发）。虽然在稳态流量或主要渗透期，较细的网格尺寸减少了数值误差，但当在高蒸散需求期（导致土壤表面附近的高压头梯度）使用较细的格栅尺寸时，误差有时会略有增加（通常小于50%）。这表明，在具有高蒸发边界条件的时期，操作员分裂误差是最显著的。通过使用渗透过程中的性能指数（网格Peclet和Courant数的乘积）或蒸散过程中的最大时间步长来约束时间步长，可以减少算子分裂误差。几个测试问题被用来为最优的空间和时间离散化提供指导。（c） 2006爱思唯尔B.V.保留所有权利。
【关键词】：coupled reactive transport; unsaturated flow modelling; geochemical reactions; operator-splitting errors

【标题】：Changes in air saturation and air-water interfacial area during surfactant-enhanced air sparging in saturated sand
【作者】：Kim, H.; Choi, K.; Moon, J. (...)
【来源】：JOURNAL OF CONTAMINANT HYDROLOGY, 2006, 88(1-2), 23-35
【摘要】：在进行空气喷射以去除含水层中的挥发性有机污染物之前，降低地下水的表面张力可以大大提高空气含量和空气喷射的影响程度。然而，关于水饱和度、由空气喷射引起的空气-水接触面积和水的表面张力之间的函数关系的详细信息尚不可用。在本研究中，使用实验室规模的砂填充柱研究了在空气喷射前向地下水中加入水溶性阴离子表面活性剂（十二烷基苯磺酸钠）对空气-水界面面积和水饱和度的影响。已经发现，水饱和度随着水的表面张力的降低而降低，直到达到这种趋势被逆转的点，使得水饱和度随着表面张力的进一步降低而增加。在45.4 dyn/cm的表面张力下实现了0.58的最低水饱和度，这被认为是本研究中使用的初始水饱和砂的最大去饱和的最佳表面张力。使用气体界面示踪剂正癸烷测量了由于空气喷射而在砂柱中产生的空气-水接触面积，发现该面积随着水饱和度的降低而单调增加。这项研究的结果为表面活性剂增强空气喷射去除含水层中的挥发性污染物提供了有用的设计信息。（c） 2006爱思唯尔B.V.保留所有权利。
【关键词】：air sparging; surfactant; aquifer; remediation; contamination

【标题】：Simulation of groundwater recharge dynamics in partially saturated fractured soils incorporating spatially variable fracture apertures
【作者】：Cey, E.; Rudolph, D.; Therrien, R.
【来源】：WATER RESOURCES RESEARCH, 2006, 42(W094139),
【摘要】：[1]本文采用数值模拟方法分析了部分饱和裂隙低渗透土壤在地下水补给过程中的复杂流动行为。目的是了解裂隙-基质相互作用和不同物理参数对流动过程的影响。三维离散裂缝模型用于模拟向包含单个垂直粗糙壁裂缝的变饱和多孔块体的渗透。结果表明，流动系统动力学对基质性质，即导水率和前期含水量，比裂缝性质更敏感。高度的裂缝-基质相互作用限制了优先沿裂缝向下流动，并降低了裂缝粗糙度的影响。研究结果表明，在模拟水流通过非饱和土壤时，可以忽略单个裂缝内的孔隙变异性。然而，裂缝在水的运动中起着重要作用，并控制着表面应用示踪剂的垂直通量。
【关键词】：

【标题】：Simulating environmental tracer transport in unsaturated-saturated porous media
【作者】：Onnis, G. A.; Althaus, R.; Klump, S. (...)
【会议】：CALIBRATION AND RELIABILITY IN GROUNDWATER MODELLING: FROM UNCERTAINTY TO DECISION MAKING 2006 The Hague, NETHERLANDS
【主办单位】：
【摘要】：自20世纪50年代初以来，人类活动已经将环境示踪剂释放到大气中。它们为研究地下动力学，特别是确定地下水停留时间提供了强大的工具。当考虑深层地下水位（地表下>10英寸）时，这些示踪剂的非饱和带输送成为正确评估地下水年龄的关键问题。穿越渗流带可能会延迟10年以上。在本工作中，我们研究了四种不同环境示踪剂H-3、He-3、Kr-85和SF6的地下动力学，特别强调了示踪剂在非饱和带中的传输。所开发的方法应用于Baltenswil含水层（瑞士）。
【关键词】：environmental tracers; groundwater; groundwater age dating; transport; unsaturated zone

【标题】：Estimation of groundwater consumption by phreatophytes using diurnal water table fluctuations: A saturated-unsaturated flow assessment
【作者】：Loheide, S. P.; Butler, J. J.; Gorelick, S. M.
【来源】：WATER RESOURCES RESEARCH, 2005, 41(W070307),
【摘要】：[1] 在许多干旱和半干旱环境中，潜水植物的地下水消耗量很难测量，但却是水预算的重要组成部分。在过去的70年里，潜水植物对地下水的消耗利用是使用一种方法来估计的，该方法分析了经过地下水位筛选的水井的水位线的日趋势（White，1932）。用这种方法获得的估计值的可靠性从未使用饱和-非饱和流模拟进行过严格评估。我们对常见的流动几何形状和一系列水力特性进行了这样的评估。结果表明，White方法的主要误差来源是比产量估计的不确定性。如果忽略排水时间和地下水位深度对特定产量的影响，White方法通常会大大高估地下水的蒸发蒸腾耗量。我们利用现成的比产量的概念作为估计适合与White方法一起使用的比产量值的基础。定义了基于沉积物质地估算现成特定产量的指南。将这些指南与White方法结合使用，应能更准确地量化地下水的蒸发蒸腾耗量。
【关键词】：

【标题】：Stochastic modeling of variably saturated transient flow in fractal porous media
【作者】：Guarracino, L.; Santos, J. E.
【来源】：MATHEMATICAL GEOLOGY, 2004, 36(2), 217-238
【摘要】：这项工作介绍了应用蒙特卡罗模拟方法对假设的随机多孔介质中的瞬态变饱和流进行统计分析。对于作为系数输入Richards流动方程的随机土壤参数的每次实现，使用混合有限元程序进行空间离散化，并结合向后的Euler和时间上的修正Picard迭代来计算压头和流场。混合法的杂交为评价单元间边界的导水性提供了一种新的方法。所提出的方法既能处理水力参数的大变异性，又能处理水力指标的分形结构。van Genuchten模型中的饱和电导率K-s和形状参数alpha（vg）被视为随机分形函数，称为分数布朗运动（fBm）或分数高斯噪声（fGn）。压头、含水量和流量分量的统计矩是通过以蒙特卡罗方式对分形参数进行平均来获得的。给出了一个数值例子，说明了所提出的方法在描述高度非均质土壤中的地下水流动方面的应用。
【关键词】：unsteady unsaturated flow; finite elements; heterogeneous soils

【标题】：Efficient steady-state solution techniques for variably saturated groundwater flow
【作者】：Farthing, M. W.; Kees, C. E.; Coffey, T. S. (...)
【来源】：ADVANCES IN WATER RESOURCES, 2003, 26(8), 833-849
【摘要】：我们考虑使用Richards方程（RE）模拟稳态变饱和地下水流动。数值求解RE的困难是众所周知的。大多数RE离散化方法都会导致非线性系统庞大且难以求解。稳态问题的非线性系统的解可能特别具有挑战性，因为稳态解的良好初始猜测通常很难获得，并且由此产生的线性系统可能难以缩放。像Picard迭代或牛顿方法的变体这样的常见方法有其优点，但在某些条件下与标准全球化技术相比表现不佳。一段时间以来，伪瞬态延拓在计算流体动力学中被用于获得牛顿方法和标准线搜索策略失败的问题的稳态解。在这里，我们研究了伪瞬态延拓以及牛顿方法与标准全球化技术相结合用于异构域中的稳态问题。我们研究了直接和预条件Krylov迭代线性求解器的方法性能。然后，我们提出了在一系列条件下获得RE稳态解的稳健有效方法的建议。爱思唯尔有限公司出版。
【关键词】：Richards' equation; steady-state; pseudo-transient continuation

        

【标题】：Solute spreading in nonstationary flows in bounded, heterogeneous unsaturated-saturated media
【作者】：Lu, Z. M.; Zhang, D. X.
【来源】：WATER RESOURCES RESEARCH, 2003, 39(10493),
【摘要】：[1] 在随机溶质（平流）传输模型中，通常假设速度场是平稳的（统计上均匀的）或平均流是单向的。在这项研究中，我们使用拉格朗日方法，建立了一个通用的随机模型，用于随机非均质多孔介质中变饱和流的输运。模型中的平均流量是多向的，速度场可以是非平稳的（具有位置相关统计）。速度场的非平稳性可能是由介质性质或复杂边界配置的统计非均匀性引起的。粒子的平均位置是使用平均拉格朗日速度来确定的。粒子扩散（位移协变量）用满足时变动力学方程的状态转移矩阵表示，该方程的系数矩阵是平均拉格朗日速度场的导数。在静止速度场的特殊情况下，转换矩阵变为相同矩阵，我们的模型简化为Dagan[1984]的著名模型。对于非平稳但单向的流场，我们的模型简化为Butera和Tanda[1999]以及Sun和Zhang[2000]的模型。通过与以下三种情况的蒙特卡罗模拟进行比较，检验了输运模型的有效性：平均重力主导流、具有地下水位边界的非饱和流和饱和-非饱和流中的输运。我们的模型结果与蒙特卡罗模拟结果非常一致。
【关键词】：solute transport; nonstationary flow; heterogeneity; stochastic processes; saturated/unsaturated flow

【标题】：Simulating interactions between saturated and unsaturated storage in a conceptual runoff model
【作者】：Seibert, J.; Rodhe, A.; Bishop, K.
【来源】：HYDROLOGICAL PROCESSES, 2003, 17(2), 379-390
【摘要】：对集水区范围内化学成分通量建模的需求越来越大。概念径流模型为此类建模任务提供了基础，前提是它们能够捕捉基本的水文过程。然而，大多数概念模型并没有完全解决不饱和和饱和储存之间的相互作用。这可能会导致对浅层地下水流域进行不切实际的模拟，其中很大一部分土壤体积可能会导致非饱和和饱和蓄水，具体取决于地下水位。向饱和蓄水中添加少量水将导致储存在非饱和区的大量水变为“饱和”。非饱和蓄水的最大体积也随着地下水位的上升而减少，即饱和蓄水量的增加。在这项研究中，提出了一个新的模型概念，其中特别强调饱和和非饱和储存之间的相互作用。总蓄水量分为两个隔间，代表饱和蓄水量和非饱和蓄水量，边界根据隔间的水量预算上下移动。地下水动力学显示出沿所研究的山坡的明显模式。靠近溪流的地区的地下水位具有类似于径流的动态，而进一步上坡的水位对降雨的反应是延迟的。为了在模型中表示这些差异，将山坡细分为河岸水库和上坡水库。将新模型的性能与更简单的模型变体的性能进行了比较，没有空间区分，也有或没有允许不饱和和饱和储存之间相互作用的新公式。新的模型方法为模拟径流和地下水动力学提供了最好的结果。山坡的细分是性能提高的主要原因。为了进一步测试模型结构，模拟流中的O-18浓度，并将其与测量值进行比较。版权所有（C）2003 John Wiley Sons，有限公司。
【关键词】：conceptual modelling; hillslope hydrology; saturated-unsaturated zones interactions; water storage; O-18

【标题】：Formulation and application of the unsaturated/saturated catchment models SUSCAT and WEC-C
【作者】：Beverly, C. R.; Croton, J. T.
【来源】：HYDROLOGICAL PROCESSES, 2002, 16(12), 2369-2394
【摘要】：本文描述了一个非饱和/饱和耦合模型框架的制定和应用，该框架旨在研究采矿对集水区水量和地下水动力学的影响。在有限元（SUSCAT）和有限差分（WEC-C）求解方案中实现了模型概念化，并发现给出了类似的结果。该模型框架模拟了地表水和地下水的耦合系统，其中基于物理的求解方案用于模拟通过非饱和带的一维运动，分布式模型用于模拟二维饱和地下水流动。每个土柱包括一系列层，每一层都连接到相邻的细胞。当一层中的任何单元形成饱和厚度时，都会考虑地下横向流动。所提供的模拟结果基于西澳大利亚达令山脉的一个集水区，该集水区被逐步开采并随后恢复。结果预测，由于采矿，矿区下方的地下水系统将出现2至4米的峰值上升。采矿六年后，植被恢复后，地下水上升已降至模拟未开采水位以上1米。据估计，采矿导致的相应流量增加在21毫米/年时达到峰值，在矿区植被恢复八年后降至7.4毫米/年。两个模型得出的模拟地下水响应和径流结果与观测数据一致。版权所有（C）2002 John Wiley Sons，有限公司。
【关键词】：unsaturated flow; recharge; stream flow; groundwater models; distributed models; catchment water balance

【标题】：Solute transport through a heterogeneous coupled vadose-saturated zone system with temporally random rainfall
【作者】：Foussereau, X.; Graham, W. D.; Akpoji, G. A. (...)
【来源】：WATER RESOURCES RESEARCH, 2001, 37(6), 1577-1588
【摘要】：The transport of nonreactive solutes through a coupled, two-dimensional, randomly heterogeneous vadose-saturated zone system subject to temporally random rainfall is predicted by Monte Carlo simulation and compared with previously published analytic results for three different rainfall patterns. The relative contributions of the uncertain inputs (i.e., rainfall and saturated conductivity) to the prediction uncertainty of solute transport are quantified in terms of the statistical moments of the pore water velocity, the plume spatial moments, and solute flux breakthrough curves at downstream control planes. Results show that the mean and variance of the saturated zone pore water velocity were approximately equivalent for the cases of uniform and random rainfall and were well predicted by the analytical relationships developed by Rubin and Bellin [1994]. As a result, the mean plume displacement, estimated by the trajectory of the mean plume center of mass, was found to be nearly identical for these cases. In the temporally random rainfall case, the saturated zone mean plume experienced more spread in the direction of mean flow at early times. However, the asymptotic rates of spatial spreading of the mean solute plumes were found to be approximately equivalent for the uniform and random rainfall cases and well predicted by the approximate expressions for longitudinal macrodispersivity in nonuniform flow proposed by Destouni and Graham [1995]. Random rainfall and random soil properties increased prediction uncertainty of the solute plume behavior in the vadose zone by an order of magnitude when compared with the uniform rain and random soil case. This effect was reduced considerably when the solute entered the saturated zone, where random rainfall produced only slightly larger prediction uncertainty than the uniform rainfall case. The analytic model developed by Destouni and Graham [1995] accurately predicted the temporal breakthrough of the mean solute plume at saturated zone control planes for all cases, if transport through the unsaturated zone accounted for the effects of temporally random rainfall using the methodology developed by Foussereau et al. [2000a, 2000b]. Results of this work indicate that for the humid climates studied here, uncertain rainfall patterns dominate transport prediction uncertainty in the shallow unsaturated zone, while uncertain solute breakthrough to the saturated zone and uncertain hydraulic conductivity dominate prediction uncertainty in the saturated zone.
【关键词】：

【标题】：Biofiltration of residual fertilizer nitrate and atrazine by Rhizobium meliloti in saturated and unsaturated sterile soil columns
【作者】：Mehmannavaz, R.; Prasher, S. O.; Markarian, N. (...)
【来源】：ENVIRONMENTAL SCIENCE & TECHNOLOGY, 2001, 35(8), 1610-1615
【摘要】：This study was undertaken to investigate whether microbial bioaugmentation of subsurface soil with subsurface irrigation could be used as a biofiltration/biocontrol technology for agricultural pollutants. Nine Plexiglas columns, 458 mm long x 139 mm in diameter, were packed with a sterilized sandy loam soil. Subsurface irrigation, through a controlled water table management system, was used to deliver bacteria, Rhizobium meliloti A-025, to the soil and to maintain aerobic (unsaturated) or anaerobic (saturated) conditions in the columns. Nitrate and atrazine, a fertilizer and a corn herbicide, were applied to the soil surface, and leaching was affected by simulated rainfall events. The soil and drainage waters were analyzed for nitrate and atrazine residues after each rainfall simulation throughout the experimental period during which the soil was kept saturated for a total of 80 days and unsaturated for a total of 70 days. The monitoring of transport and survival of the implanted bacterial strain (A-025) showed that subsurface irrigation was successful in introducing and transporting the bacteria throughout the soil columns. During the saturated period, significantly more (95% probability) nitrate-N leached into the drainage waters from the control columns than from the bioaugmented columns; the increase being 450% or more for the abiotic control columns. The amount of atrazine that leached into the drainage waters during the unsaturated period was also significantly more from control columns as opposed to bioaugmented columns, with the increase being 262%.
【关键词】：

【标题】：A fully coupled finite element analysis of water-table fluctuation and land deformation in partially saturated soils due to surface loading
【作者】：Kim, J. M.
【来源】：INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2000, 49(9), 1101-1119
【摘要】：针对部分饱和土壤中地表荷载引起的地下水位波动和土地变形，建立了一个全耦合的数值模型。该数值模型是基于变饱和多孔介质中地下水渗流的孔隙弹性控制方程和Galerkin有限元方法建立的。针对具有六种不同初始地下水位高程的土柱表面荷载的一维固结问题，验证了数值模型。然后将该数值模型应用于建筑工地部分饱和土壤表面荷载的二维固结问题。这两个问题的数值模拟结果表明，部分饱和土壤的地下水位波动，地下水位以上的非饱和区对部分饱和土壤在表面荷载下的固结行为有显著影响。这种影响是由部分机械加载应力的永久吸收以及固体骨架变形场与非饱和带中的地下水流场之间由于其部分饱和而产生的弱流体力学耦合引起的。版权所有（C）2000 John Wiley&Sons，有限公司。
【关键词】：unsaturated soil; surface loading; water-table fluctuation; land deformation; poroelasticity theory; finite element simulation

【标题】：Saturated area formation on nonconvergent hillslope topography with shallow soils: A numerical investigation
【作者】：Ogden, F. L.; Watts, B. A.
【来源】：WATER RESOURCES RESEARCH, 2000, 36(7), 1795-1804
【摘要】：饱和区形成的预测对于浅层、高透水性土壤流域的水文建模非常重要。这项模拟研究考察了湿润过程中山坡特性和降雨量对饱和源区演变的相对重要性。这项研究的重点是具有恒定坡度的均匀、不收敛的山坡地形。二维可变饱和地下水模型VS2D[Healy，1990；Lappala等人1993]用于模拟稳定均匀降雨作用下饱和区域的形成。该研究方法系统地改变了四种山坡特性，即不透水层深度、坡长、坡角、平均饱和导水率和降雨量。结果表明，在平衡时饱和的山坡长度的分数是参数Phi的函数，参数Phi定义为降雨量乘以斜坡长度，再除以斜坡角度、土壤厚度和饱和导水率。从平衡时间的角度分析了饱和区的时间演化；非饱和带的非线性和山坡性质的差异导致了非唯一的平衡时间关系。当倾向于引起表面饱和的因素与倾向于消散表面饱和的因子近似平衡时，达到平衡的时间最大。润湿过程中表面饱和度的时间演变遵循多种轨迹。在不同性质的斜坡上，等效的山坡土壤蓄水率会导致广泛的地表饱和条件。土壤较浅、坡角较小的山坡最容易形成饱和区。然而，当1<Phi<6时，陡峭山坡上的表面饱和度可能会发生突然变化。
【关键词】：

【标题】：Modeling of pyrite oxidation in saturated and unsaturated subsurface flow systems
【作者】：Xu, T. F.; White, S. P.; Pruess, K. (...)
【来源】：TRANSPORT IN POROUS MEDIA, 2000, 39(1), 25-56
【摘要】：黄铁矿氧化（FeS2）导致金属的酸化和迁移。可变饱和反应流系统中黄铁矿氧化的数学建模具有挑战性，因为（1）它是通过多相流和传输过程的复杂相互作用发生的，以及（2）在不同的氧化还原条件下，关键物种的水浓度变化超过几十个数量级。在这里，我们提出了一个氧化还原过程的通用多相反应传输模型。在TOUGREACT和TOUGH2-CHEM仿真代码中分别使用顺序迭代和同时求解进行了两种替代实现。这两个代码都用于模拟具有简单一维流动和反应条件的完全饱和和可变饱和黄铁矿氧化问题。两种代码的结果表明，在环境条件下，由于反应导致的氧分压降低对流体流动的影响并不显著。然而，必须注意的是，当流体流动和化学反应强烈耦合时，例如当地热储层发生沸腾时，这可能是至关重要的。完全同步的方法有一个完整的过程描述。序列迭代方法被发现在计算上更有效。氧气扩散过程在非饱和条件下黄铁矿氧化的化学演化中起主导作用。以二维裂隙岩体为例，说明了在复杂的流动和地球化学条件下黄铁矿的氧化作用。这个例子表明，在可变饱和流动系统中，黄铁矿氧化对水文地球化学演化有很大影响。预测的原生岩石矿物的蚀变和次生矿物组合的发育与现场观测一致。这个例子是氧化风化过程的原型，对地球科学、工程和环境应用具有广泛意义。
【关键词】：pyrite oxidation; acid mine drainage; saturated-unsaturated media; oxygen transport; numerical approach; reactive chemical transport

【标题】：Coupled unsaturated/saturated groundwater models assessing heavy metal transport
【作者】：Nützmann, G.; Holzbecher, E.; Ginzel, G. (...)
【会议】：COMPUTATIONAL METHODS IN WATER RESOURCES, VOLS 1 AND 2: COMPUTATIONAL METHODS FOR SUBSURFACE FLOW AND TRANSPORT 2000 UNIV CALGARY, FAC SCI, CALGARY, CANADA
【主办单位】：
【摘要】：在实际应用中，部分饱和多孔介质的三维问题通常可以归结为非饱和带的垂直流和饱和含水层的水平流。本文提出了一种使用一维非饱和和二维饱和输运模型模拟重金属在非均质区域中输运的方法。两个模型之间的耦合是使用一个新的实现界面进行的，该界面必须以地下水模型可以处理的一般形式转换非饱和带重金属浸出的穿透曲线。证明了所提出的模型的适用性，并显示出对复杂生态水文系统模拟的可能推广。
【关键词】：

【标题】：Numerical modeling of volatile and slightly soluble contaminant flow in unsaturated and saturated zone
【作者】：Itoh, K.; Tosaka, H.; Matsubara, O. (...)
【会议】：GROUND WATER UPDATES 2000 OMIYA, JAPAN
【主办单位】：
【摘要】：为了预测几种化学污染物对土壤和地下水的污染，作者开发了考虑污染物向空气中溶解和挥发的多相型数值模拟技术，将其建模为不同相之间的瞬态扩散。进行了实验室实验，以检查挥发性物质在不饱和区的行为。根据实验室实验和数值模拟的比较结果，我们可以重现测量的瞬态浓度，并验证了瞬态挥发模型。
【关键词】：groundwater contamination; volatilization; dissolution; numerical simulation

【标题】：A numerical model for density-and-viscosity-dependent flows in two-dimensional variably saturated porous media
【作者】：Boufadel, M. C.; Suidan, M. T.; Venosa, A. D.
【来源】：JOURNAL OF CONTAMINANT HYDROLOGY, 1999, 37(1-2), 1-20
【摘要】：我们提出了二维变饱和介质中水流和溶质传输的公式，该公式考虑了溶质对水密度和粘度的影响。控制方程采用无量纲形式，取决于六组无量纲参数。这些方程在空间上使用Galerkin有限元公式进行离散，并在时间上使用具有质量块的反向Euler格式进行积分。采用改进的Picard方法对水流方程进行线性化。通过与已发表的数值结果进行比较，验证了由此产生的数值模型，即MARUN模型。然后，它被用于研究海水浓度（约30 g l（-1））下的海滩水力学，以提供营养物质，对海滩上的石油泄漏进行生物修复。我们在假设海滩的矩形截面上进行的数值模拟表明，在质地细、各向异性比低和/或物理分散性低的土壤中，非饱和区的浮力非常重要。在这种情况下，在淡水溶液中向受污染的海滩施用溶解的营养物质优于在海水溶液中施用。在我们考虑的情况下，浓度产生的粘度效应相对于浓度产生的密度效应可以忽略不计。（C） 1999 Elsevier Science B.V.保留所有权利。
【关键词】：beaches; interface; freshwater; saltwater; bioremediation

【标题】：A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case
【作者】：Yakirevich, A.; Borisov, V.; Sorek, S.
【来源】：ADVANCES IN WATER RESOURCES, 1998, 21(8), 679-689
【摘要】：提出了一种准三维（quasi 3-D）模型，用于模拟土壤非饱和区和饱和区的地下水流动和溶质运移。该模型基于非饱和带的垂直方式和基本水平的地下水流动的假设。非饱和带的一维Richards方程在浸润面与饱和带的二维流动方程耦合。后者是通过对含水层厚度范围内饱和区的三维流动方程进行平均而获得的。与渗漏潜水含水层的Boussinesq方程不同，所开发的模型不包含具有特定产量的储存项和自然补给的源项。相反，它包括潜水表面的水通量项，理查兹方程通过该项与地下水流动方程相联系。饱和区的垂直水通量是根据流体质量平衡方程计算的，而饱和区的水平水通量则是根据达西定律计算的。采用三维输运方程模拟溶质运移过程。提出了一种求解一般拟三维情形的数值算法。所开发的方法以准二维横截面情况（QUASI2D）为例。对三个合成问题的模拟表明，QUASI2D和两个全二维流和传输代码（SUTRA和2DSOIL）获得的结果之间具有良好的一致性。然而，使用QUASI2D代码的模拟比使用SUTRA和2DSOIL代码的模拟快几倍。（C） 1998爱思唯尔科学有限公司。保留所有权利。
【关键词】：subsurface flow; transport; unsaturated and saturated zones; phreatic aquifer; finite difference solution; quasi three-dimensional approach

【标题】：Sensitivity and moment analyses of head in variably saturated regimes
【作者】：Li, B.; Yeh, T.
【来源】：ADVANCES IN WATER RESOURCES, 1998, 21(6), 477-485
【摘要】：提出了一种近似多维多孔介质中变饱和流水头统计矩的数值方法。该近似依赖于有限元流模型的一阶泰勒级数展开和变饱和流的伴随态数值方法来评估灵敏度。该方法可用于分析与具有任何类型的边界和初始条件的变饱和多孔介质中稳态或瞬态流压头预测相关的不确定性。从而减轻了随机分析方法的局限性，如谱/扰动方法和耗时的蒙特卡罗模拟技术。举例说明了该方法的实用性，并研究了变饱和流态下水头变化的时间演变。结果表明，地下水位的波动会对水头变化的传播产生显著影响。（C） 1998爱思唯尔科学有限公司。保留所有权利。
【关键词】：

【标题】：Partially saturated hydrological flows: Numerical experiments and analyses
【作者】：Ababou, R.; Tregarot, G.; Larabi, A.
【会议】：COMPUTATIONAL METHODS IN CONTAMINATION AND REMEDIATION OF WATER RESOURCES: PROCEEDINGS OF 12TH INTERNATIONAL CONFERENCE ON COMPUTATIONAL METHODS IN WATER RESOURCES, VOL 1 1998 IRAKLION, GREECE
【主办单位】：
【摘要】：水文流以不同的形式出现在地面和地面以下：地表径流、径流、渗透、蒸发、无侧限地下水流、上层滞水水位、各种水文界面处的未决、渗透和渗出。本文的重点是涉及自由表面（地下水位）的饱和/非饱和耦合流。我们的方法基于变饱和流的混合公式和Darcy Buckingham通量梯度定律的广义形式。使用三维有限体积代码（Bigflow）对几个测试案例进行了建模：1）含有粘土透镜体的沙质土壤中的池塘渗透（上层滞水和反向地下水位）；2） 降雨期间/之后，地下水向溪流的补给/排放；3） 饱和锋通过河岸的横向推进（Polubarinova Kochina的湿润锋问题，1962）；4） 细多孔层上大孔油藏地下水位的垂直排水。试验1表示，由于河流水位上升（洪水），泛滥平原上的积水渗透。测试2可用于分析连接到地下水位的流的排放速率。试验3可以被视为一种模拟河流中洪水波通过期间河岸蓄水的方法。试验4是覆盖阻碍多孔层的柱中1D地下水位排水的部分饱和模型。对数值结果进行了分析，并与解析解进行了比较。一些模拟结果（测试4）也已经用三维有限元代码Geoprof进行了基准测试。这项工作是更广泛的研究工作的一部分，旨在制定、参数化和建模通过地表/地下水文盆地的空间分布耦合流。
【关键词】：

【标题】：A case study simulation of DBCP groundwater contamination in Fresno County, California - 2. Transport in the saturated subsurface
【作者】：Loague, K.; Abrams, R. H.; Davis, S. N. (...)
【来源】：JOURNAL OF CONTAMINANT HYDROLOGY, 1998, 29(2), 137-163
【摘要】：本文是模拟加利福尼亚州弗雷斯诺县因长期非点源应用1,2-二溴-3-氯丙烷（DBCP）而造成的区域性地下水污染的多页系列文章的第二部分。在本研究的第一部分[League，K.，Lloyd，D.，Nguyen，A.，Davis，S.N，Abrams，R.H.，1997。加利福尼亚州弗雷斯诺县DBCP地下水污染的案例研究模拟：1。通过非饱和地下进行浸出。J.康塔明。Hydrol。（本期）]，我们的重点是不饱和带和DBCP浸出到地下水位。本文的重点是评估DBCP在饱和地下的传输。这里，地下水流量代码MODFLOW与溶质运移代码MT3D相耦合，以模拟弗雷斯诺县研究区地下水位以下35年内DBCP的潜在命运和运移。本文的结果说明了研究区DBCP羽流的长期区域尺度演变。模拟的DBCP浓度仅限于相对较浅的较年轻沉积物，通常远低于为该化学品设定的最大污染物水平。这里提供的模拟表明，DBCP的非点源应用不是研究区域中观察到的热点的原因。在本文的最后，我们讨论了本文和上述配套论文（Loague等人，1997）中报道的弗雷斯诺案例研究的含义和扩展。“消除不可能的东西。剩下的东西，无论多么不可能，都必须是真相。”——夏洛克·福尔摩斯（C）1998年爱思唯尔科学B.V。
【关键词】：nonpoint sources (NPS) of groundwater contamination 1,2-dibromo-3-chloropropane (DBCP); Fresno County, CA, USA

【标题】：Monotonicity considerations for saturated-unsaturated subsurface flow
【作者】：Forsyth, P. A.; Kropinski, M. C.
【来源】：SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1997, 18(5), 1328-1354
【摘要】：证明了单调性是保证离散非线性不饱和Aom方程不产生新的非物理局部极大值和极小值的充分条件。推导了迁移率项的各种加权类型的单调性条件。基本离散化是有限元类型的，但结果可以扩展到有限体积离散化。中心权重仅是条件单调的，而上游权重是无条件单调的。给出了高度异构问题的示例计算。
【关键词】：saturated-unsaturated flow; monotonicity

【标题】：Conditional stochastic analysis for multiphase transport in randomly heterogeneous, variably saturated media
【作者】：Zhang, D. X.
【来源】：TRANSPORT IN POROUS MEDIA, 1997, 27(3), 265-287
【摘要】：地下水有机物污染最近已成为一个日益引起人们关注的健康和环境问题。通常，必须研究非水相液体（NAPL）、水和空气的多相系统，以便真实地描述这些材料在地下的运动。已经开发了许多模型来研究多孔介质中的多相流和/或多相输运。然而，使用模型来研究介质非均质性对这种流动和传输的影响只是最近才发生的事情。对于饱和和非饱和介质中的单相流动和输运，已经证明了介质非均质性的研究适用于随机分析。在本文中，我们将单相输运的欧拉-拉格朗日随机理论推广到在流动稳定且相处于局部化学平衡的条件下，随机非均匀介质中的多相多相输运问题。我们提出了理论表达式来描述任意相中任意物种随机浓度的前两个条件矩。尽管它们揭示了一些基本性质，并有助于深入了解问题的本质，但如果没有高分辨率蒙特卡罗模拟或近似（闭包），就无法评估这些表达式。因此，我们提出了两组可行的近似，一组是弱近似，另一组是线性化的伪菲基近似。前者产生第一条件矩的非线性积分微分方程，而后者产生线性微分方程。然后，可以从弱近似或伪菲基近似的显式表达式中计算二阶矩。
【关键词】：stochastic analysis; conditional probability; multiphase transport; heterogeneous porous media

【标题】：Linearized cosimulation of hydraulic conductivity, pressure head, and flux in saturated and unsaturated, heterogeneous porous media
【作者】：Harter, T.; Gutjahr, A. L.; Yeh, T.
【来源】：JOURNAL OF HYDROLOGY, 1996, 183(1-2), 169-190
【摘要】：开发了一种有效的共模拟器，用于在饱和或非饱和条件下生成随机水力特性场和相关流态。该联合模拟器结合了基于快速傅立叶变换技术的谱随机场发生器和饱和和非饱和多孔介质中流动的一阶扰动/谱解。由于模拟器的一阶近似，模拟了几种不同变率的地质介质中的流态，以研究模拟器的准确性。对于轻度和中度非均质地质介质，发现模拟器在压头场和通量分布方面非常准确。此外，该模拟器的执行时间大大小于任何经典数值模拟器。然而，随着地质介质变得高度不均匀，模拟器的精度会降低。
【关键词】：

【标题】：MODELING MULTICOMPONENT REACTIVE CHEMICAL-TRANSPORT IN NONISOTHERMAL UNSATURATED SATURATED SOILS .2. NUMERICAL SIMULATIONS
【作者】：WU, G.; CHIENG, S. T.
【来源】：TRANSACTIONS OF THE ASAE, 1995, 38(3), 827-838
【摘要】：用Galerkin有限元法求解了饱和/非饱和土壤条件下水分、热量和活性化学物质同时传递的偏微分方程组（Wu和Chieng，1995）。用迭代Newton-Raphson方法求解了化学反应的非线性代数方程。化学输运方程是根据每个组分物种的总分析浓度来制定的，并且可以顺序求解。实现了求解线性方程组的三种方法：一维和二维小问题的带状矩阵法、二维和三维大问题的迭代法和稀疏矩阵法。给出了两种高效可靠的迭代方法：对称矩阵方程的共轭梯度法和一般非对称矩阵方程双共轭梯度平方法。使用开发的计算机程序解决了一个典型的地下排水问题。运行了两个模拟来研究程序的行为。两种模拟都产生了合理和令人满意的结果。
【关键词】：MATHEMATICAL MODEL; NUMERICAL SIMULATION; MULTICOMPONENT; WATER FLOW; HEAT TRANSFER; SOLUTE TRANSPORT; CHEMICAL REACTIONS

【标题】：STOCHASTIC SIMULATION OF THE REGIONAL PESTICIDE TRANSPORT INCLUDING THE UNSATURATED AND THE SATURATED ZONE
【作者】：FRANKE, H. J.; TEUTSCH, G.
【来源】：ECOLOGICAL MODELLING, 1994, 75, 529-539
【摘要】：欧洲共同体新的饮用水标准允许农药总浓度最高为0.5微克/升。地方、州和联邦一级的水利工程和环境机构被迫严格评估各种自然条件下农药的污染潜力。因此，开始了大量的农药研究，主要是在实验室规模使用柱或蒸渗仪实验，或在集水区规模使用地下水质量数据。然而，也需要采用耦合方法，将实验室数据纳入区域集水区研究，同时考虑到主导交通参数的空间和时间变异性和/或不确定性。本文提出了一种随机方法，将一维农药运移数值模型与三维地下水流动和运移数值模型相耦合。通过在单个农田的尺度上对一维参数集进行多次实现，考虑了土壤参数的可变性。含水层（饱和带）的区域异质性可以通过钻孔和结构数据来确定地描述，也可以使用随机场生成器（如Turning Bands方法）来描述。根据文献中的一组土壤数据，分析了不同土壤参数分布在假设地下水抽取井的浓度突破曲线方面的重要性。这项分析旨在优化数据收集策略，用于未来的农药渗滤液风险评估。使用所描述的耦合随机建模方法，估计了生产井中农药浓度的敏感性，这是土壤和含水层输送特性变化的结果。
【关键词】：GROUNDWATER; PESTICIDES

【标题】：CHEMICAL FATE AND TRANSPORT IN A DOMESTIC SEPTIC SYSTEM - APPLICATION OF A VARIABLY SATURATED MODEL FOR CHEMICAL MOVEMENT
【作者】：SHUTTER, S. B.; SUDICKY, E. A.; ROBERTSON, W. D.
【来源】：ENVIRONMENTAL TOXICOLOGY AND CHEMISTRY, 1994, 13(2), 223-231
【摘要】：对单户化粪池系统对地下水质量的影响进行了详细的数学建模研究。这项工作有助于提高对与化粪池系统相关的地下过程的概念理解，并深入了解化粪池系统衍生化学品的命运。与化粪池系统相关的化学羽流的二维饱和-非饱和模拟在密切匹配的现场观测中取得了成功。结果表明，羽流受到有限的分散（或扩散），反映了含水层物质的均匀性。沿地下水位附近羽流上边缘的急剧浓度梯度表明垂直混合较弱，并意味着进入化粪池系统含水层降级区的补给水稀释有限。校准的分散度纵向值（α（L））为10cm，横向垂直值（a（r））为1mm。模拟消费洗涤剂产品中存在的活性化学成分进入化粪池系统的运动，以评估吸附和生物降解的影响。发现洗涤剂化合物线性烷基苯磺酸钠（LAS）和次氮基三乙酸（NTA）的迁移主要受这些过程的控制，该模型预测的有限移动与现场观察到的一致。对于不饱和区中半衰期等于15天和分布系数（K-d）值等于1.5 L/kg的情况，预计LAS从瓷砖排水沟移动<0.5 m。然而，半衰期明显较长的化学物质或在土壤基质中分配有限的化学物质可能会迁移得更远。
【关键词】：SEPTIC SYSTEM; GROUND WATER; CONTAMINANTS; MODEL; SAND

【标题】：A NUMERICAL-MODEL FOR WATER-FLOW AND CHEMICAL-TRANSPORT IN VARIABLY SATURATED POROUS-MEDIA
【作者】：YEH, T.; SRIVASTAVA, R.; GUZMAN, A. (...)
【来源】：GROUND WATER, 1993, 31(4), 634-644
【摘要】：建立了一个二维数值模型，用于模拟水在变饱和多孔介质中的流动和化学物质传输。采用Galerkin有限元技术，采用Picard迭代法或牛顿迭代法求解非线性流动方程。将伽辽金技术分别应用于达西方程，得到了一个连续的速度场。使用具有一阶损失项的双位点吸附-解吸模型来描述反应性溶质的化学行为。传输方程的平流部分采用一步后向粒子跟踪求解，色散部分采用正则伽辽金有限元技术求解。为了节省计算机内存和执行时间，采用预条件共轭梯度法迭代求解线性联立方程组。将该模型应用于几个流动和输运问题，并将数值结果与观测值和解析值进行了比较。该模型被发现与分析值和观测值非常相似，甚至在非常尖锐的锋面附近也是如此。
【关键词】：

【标题】：IMPROVED 3-DIMENSIONAL FINITE-ELEMENT TECHNIQUES FOR FIELD SIMULATION OF VARIABLY SATURATED FLOW AND TRANSPORT
【作者】：PANDAY, S.; HUYAKORN, P. S.; THERRIEN, R. (...)
【来源】：JOURNAL OF CONTAMINANT HYDROLOGY, 1993, 12(1-2), 3-33
【摘要】：对高度非均质地下介质中饱和-非饱和地下水流动和污染物迁移的精确三维模拟需要求解数值技术的极端灵活性。模拟大规模复杂现场问题的主要局限性是可用的计算能力，以及校准、验证和执行这些系统的灵敏度分析所需的大量数据管理工作。提出了一种有效的有限元模型，它克服了解决复杂问题时遇到的许多挑战。采用Galerkin和上游加权残差法求解流动和输运方程。对于高度非线性的流动问题，使用弦斜率Newton-Raphson线性化和相对渗透率的上游加权的选项是必要的。否则，使用鲁棒性较差的Picard方案可以实现计算节省。欠松弛公式进一步缓解了非线性方程组求解中的收敛困难。六面体正交曲线单元用于准确有效地离散具有不规则地层几何形状的区域，过渡单元用于对有限元网格进行分级，仅在高梯度区域进行更精细的离散。使用影响系数公式生成元素矩阵，从而避免了昂贵的数值积分。提供了7点、11点和27点晶格连接结构，作为提高求解精度或提高计算速度和存储要求的选项。Galerkin和上游加权方案在不过度涂抹浓度的情况下缓解了求解输运方程的稳定性问题。预处理共轭梯度（PCG）和ORTHOMIN技术分别用于求解由元素流和输运方程的组装产生的大型对称和非对称矩阵方程，如果情况允许，则从计算中消除重复的矩阵组装和分解。讨论了整个模型的组装，并参考先前发表的工作，了解已经测试和实施的方案的细节。本文仅详细介绍了这些方案的进一步发展。通过选定的模拟实例证明了该模型的应用，包括评估多层无侧限含水层系统上方浅层废物处理设计的水分运动和污染物迁移。
【关键词】：

【标题】：UNSATURATED AND SATURATED FLOW IN RESPONSE TO PUMPING OF AN UNCONFINED AQUIFER - NUMERICAL INVESTIGATION OF DELAYED DRAINAGE
【作者】：AKINDUNNI, F. F.; GILLHAM, R. W.
【来源】：GROUND WATER, 1992, 30(6), 873-884
【摘要】：在这项研究中，使用数值方法来检验Nwankwor等人提出的毛细管边缘假说。（1992）解释抽水过程中无侧限含水层储存特性的变化。根据现场观测，该假设表明，无侧限含水层对抽水施加的应力的响应在很大程度上受移动地下水位上方形成的垂直水力梯度的大小控制。在通过一维滞回实验室实验和二维轴对称变饱和半分析模型进行验证后，使用数值解模拟Nwankwor等人的抽水试验。（1992）。用于模拟的数值参数是从独立于泵送试验的实验中选择的。在不同深度和距抽油井不同径向距离的情况下，数值结果与现场数据之间有很好的一致性。从数字上获得的地下水位以上垂直梯度幅度的变化以及由此产生的表观比产量值的变化也显示出与现场数据相似的趋势。数值结果和现场结果之间的良好一致性提供了令人信服的证据，证明现场观测和由此产生的关于比产量时间变化的毛细管条纹假说确实符合可变饱和流的物理原理。通过将早期数据拟合到类型曲线中经常获得的高蓄水性值与水力梯度增加的时期相对应，可能与含水层的固有特性（如压缩性）无关。与Nwankwor等人的“体积平衡法”一致。（1984），结果还表明，无侧限含水层的最终比产量与可排水孔隙度的关系比通常从类型曲线分析中获得的低值更大。
【关键词】：