对土壤水分入渗特性有重要影响.DOC

应用HYDRUS-1D模拟砂质夹层土壤入渗特性 范严伟1,2, 毕贵权1, 赵文举1 （1. 兰州理工大学能源与动力工程学院, 兰州 730050 ; 2.西部灾害与环境力学教育部重点实验室/兰州大学土木工程与力学学院, 兰州 730000） 摘 要 ：依据非饱和土壤水分运动理论，采用HYDRUS-1D软件，对砂质夹层土壤入渗特性进行数值模拟，分析各因素对砂质夹层土壤入渗规律的影响。结果表明：砂质夹层结构对土壤入渗特性有较大影响，具有暂时的阻水和减渗作用；湿润锋穿过砂层上界面后，入渗过程变为稳渗阶段，稳渗率主要受砂层质地、砂层埋深和压力水头影响，与土壤初始含水率和砂层厚度无关；砂质夹层土壤剖面水分分布不连续，上层土壤基本饱和，砂层土壤未饱和，土壤剖面含水率主要受砂层质地、砂层埋深和砂层厚度的影响。研究结果可为农业水资源利用及工程防渗技术提供理论依据。 关键词 ：砂质夹层; 影响因素; 入渗特性; 数值模拟; HYDRUS-1D 中图分类号 ：S152.7 由于风沙、水文和地质等外力作用，田间土壤剖面普遍存在砂质夹层结构[1-2]。夹层结构会改变土壤水分的入渗性能和分布状况，对土壤入渗特性有重要影响[3-6]。与均质土壤相比，砂质夹层土壤的入渗特性主要受砂层质地、埋深、厚度以及土壤初始含水率和压力水头等因素影响。许多学者采用室内试验方法，研究了砂层质地、埋深和厚度对土壤入渗特性的影响[7-9]。随着计算机技术的发展，数值模拟方法越来越多地应用于土壤水分运动的研究中，成为定量分析土壤水分入渗特性的有效工具[10-12]。其中，美国农业部盐渍土实验室开发的HYDRUS-1D软件，可对一维饱和-非饱和土壤中水、热及溶质运移进行数值模拟[13]。HYDRUS-1D应用广泛，具有良好的适用性，可用来模拟层状土壤水分入渗特性[14-16]。本文依据非饱和土壤水分运动理论，采用HYDRUS-1D软件，对不同影响因素组合下的砂质夹层土壤入渗特性进行数值模拟，分析各影响因素条件下砂质夹层土壤的入渗规律，为农田水资源高效利用和工程防渗技术提供理论依据。 1 层状土一维垂直入渗数学模型 1.1 基本方程 HYDRUS-1D中，一维饱和-非饱和土壤水运动采用Richards方程进行描述。其表达式为： （1） 式中：为垂向坐标，规定向下为正，cm；为土壤比水容量，1/cm；为在饱和区和非饱和区分别为压力水头和基质势，cm；为入渗时间，min；为土壤非饱和导水率，cm/min。 式（1）中涉及的土壤水分特征曲线和土壤非饱和导水率采用van Genuchten-Mualem模型拟合[17]。即 （2） （3） 式中：为土壤残余含水率，cm3/cm3；为土壤饱和含水率，cm3/cm3；、和为土壤物理特性有关的拟合参数；；，=1-1/；为土壤饱和导水率，cm/min；=0.5。 图1 砂质夹层土壤计算简图 Fig.1 Calculation diagram of soil with sand interlayer 1.2 定解条件 图1为砂质夹层土壤水分运动模拟计算简图。试验开始时，土壤水分剖面为稳定剖面，计算域内各点土水势相等。即 （4） 式中：为土壤初始总水势，cm； 试验中，土柱上边界保持恒定水头入渗，下边界入渗水量未到达。即 （5） （6） 式中：为压力水头，cm；为土柱高度（比湿润锋所湿润范围大）；cm。 1.3 数值求解 利用HYDRUS-1D软件，对层状土一维垂直入渗的数学模型进行求解。模拟不同砂层质地、初始含水率、砂层厚度、砂层埋深和压力水头条件下的砂质夹层土壤入渗特性。求解过程中，采用Galerkin有限元法对土壤剖面进行空间离散，采用隐式差分格式进行时间离散。模拟土层深度为100 cm，时间步长为0.01 min，空间步长为1 cm，模拟历时200 min。土壤质地通过van Genuchten-Mualem模型中的参数来体现，土壤初始含水率、砂层埋深和砂层厚度通过初始条件来设定，而压力水头通过边界条件来实现。 为保证研究土壤的广泛性和成果的普适性。模拟中土壤的van Genuchten-Mualem模型水力特性参数分别取自参考文献[6]、[13]和[15]。见表1。 表1 不同土质van Genuchten-Mualem模型水力特性参数 Table.1 Hydraulic parameters in van Genuchten-Mualem model of differ