HYDRUS-1D 软件应用

1. HYDRUS-1D软件应用 马文翠

2. 软件简介 • HYDRUS-1D是美国农业部盐渍土实验室开发的模拟非饱和土壤中的水、热、溶质运移的软件，它在模拟土壤中水分运动、盐分、污染物(如农药)和养分(如土壤氮素)运移方面得到广泛应用。软件分为1维、2维、3维三种，分别命名维Hydrus-1D、Hydrus-2D、Hydrus-3D。由于非饱和土壤水主要是1维垂向运动的形式，因此Hydrus-1D的应用非常广泛。

3. 1 m 1.5 m 根系层 1.细砂壤土 1.5 m 2.中砂 问题描述 • 以天津一块农田为背景，对2012年冬小麦收割期农田尺度下土壤水分运动进行数值分析 • 土壤剖面从地面算起的深度： 300cm • 土壤分层：分两层 • 准备模拟哪个时间段水分变化：0~30d • 根系分布：1m • 初始地下水位：地面以下6m

4. 软件界面 模拟计算 前处理 后处理

5. 主要过程

6. 水分运移方程 • 将坐标原点选在地面，取z轴向下为正，则一维饱和—非饱和带水分运移基本方程为Richards方程： • 式中：θ为含水率；k（θ）为水力传导度：h为负压水头；S 为植物根系吸水量，对裸露区为0。

7. 模型有2个土层 长度单位是cm 土壤剖面的深度是300 cm 几何信息

8. 质量平衡信息 土壤分区长度 [L] 分区水量 [L] 流入或流出量 [LT-1] 平均压力水头[L] 实际表面通量 [LT-1] (入渗/蒸发: -/+) 剖面底实际通量[LT-1] (流入/出: +/-) 整个区域的水质量平衡绝对误差[L] 整个区域的水质量平衡相对误差 [%]

9. 输出信息 选择输出时间点 输出5组结果，每6天一组

10. 最多迭代次数 含水量迭代精度 压力水头迭代精度 增大步长迭代次数信号 缩小步长的迭代次数信号 步长增大比例 步长缩小比例 用于生成水分特征曲线的散点值 最小吸力间距 最大吸力间距 迭代参数

11. 上、下边界条件 定压力水头 定通量 大气边界，可积水 大气边界，产流 变水头 变水头/通量 定压力水头 定通量 变水头 变通量 自由排水 深度排水 渗透面 水平排水 初始条件 使用压力水头 使用含水量 输入Etp和LAI 消光系数 最大积水厚度

12. 消光系数 利用比尔定律计算消光系数: 定义：植物群体的日同化量大都表示为冠层截获光合有效辐射(PAR)量的函数，而群体对PAR的截获量可通过到达冠层顶部的PAR在通过作物群体时的衰减进行计算，衰减程度用消光系数(K)表示。它取决于太阳角度、植被类型及叶片空间分布特征。 式中：TPAR是冠层底测得的瞬时光合有效辐射(μmol·m-2·s-1)，PAR是冠层顶测得的瞬时光合有效辐射，LAI是观测日的叶面积指数。 利用冠层上、下测得的光合有效辐射资料，根据其变化规律拟合方程求积分，得到入射PAR的日总量(μmol·m-2·d-1)和透射的PAR日总量，再根据当日的叶面积指数来计算日平均消光系数。

13. 潜在蒸散量处理 • 植被蒸腾量和土壤蒸发量分开处理，HYDRUS-1D推荐使用一个经验公式来把潜在蒸散量分割为蒸腾潜力和土壤蒸发潜力： • 其中ETp为潜在蒸散量（cm/d），Tp为潜在蒸腾量(cm/d)，Ep为土壤潜在蒸发量(cm/d)，LAI是叶面积指数，k为消光系数。

14. 上边界通量

15. 随时间变化边界条件

16. 最小压力水头值(hCritA) • 即地面土壤达到最干燥状态时的压力水头。理论上，土壤十分干燥时，吸力很大，而液态孔隙水的压强很小，与空气湿度保持平衡关系，有 • hA为最小压力水头，Hr为空气绝对湿度，RT/Mg为空气的摩尔气体常数。公式需要近地面空气湿度，一般可取饱和水汽湿度，因2 cm以下土壤空气湿度往往饱和，随温度发生变化，可据近地面气温变化推算地表土壤空气湿度，再换算成压力水头。HYDRUS中需输入hA的绝对值，缺省值为 • hCritA=|hA|=106 cm • 此数值只对土壤蒸发起作用。建议：hCritA 对应的土壤含水量应至少比残余含水量大0.005，在模拟根系吸水时，hA还应低于P3。否则，当根系吸水的临界压力水头和地面蒸发的最小压力水头满足时,会导致回流现象，不合理。

17. 模拟30 d内的变化 时间单位 时间步长信息 上边界是随时间变化的，每天一组数据，共30组数据 自动处理蒸腾量在每天24小时内的变化 时间信息

18. 蒸腾量处理 其中Tp(t)是瞬时潜在蒸腾量，t为时间。模型假设早上6点之前以及晚上18点-24点的蒸腾量总和只占全天蒸腾量的1%。 HYDRUS-1D 可以使用一个经验公式来处理每天24小时潜在蒸腾量的变化，设某天的潜在蒸腾量为 （cm/d），则

19. 单孔介质模型 双重介质模型 双重孔隙度 双重渗透性 主要在溶质运移问题中使用 吸湿和疏干滞后过程模型 土壤水力模型

20. van Genuchten-Mualem公式处理水分特征曲线 θs – 饱和含水量 [-] θr – 残余含水量 [-] a, m, n ，l- 经验参数 [1/L], [-], [-],[-] Se – 有效含水量 [-] Ks – 渗透系数 [L/T] K (h) – 压力水头为 h时的非饱和导水率 [L/T]

21. Sandy loam sand 编号 r s (cm1) n Ks (cm/d) l 土壤数据库 水力参数 θr： 土壤残余含水量 θs： 土壤饱和含水量 Alpha、 N：土壤持水函数中的参数[L-1] Ks：饱和导水率，即渗透系数 [LT-1] l：　导水函数中的曲度参数 [-]

22. 神经网络预测模块 99-101% 0.5-2.0 0.0-1.0

23. 水分胁迫模型 盐分胁迫模型 补偿吸水域值 （临界湿润度） 根系吸盐模型 根吸水

24. 根吸水模型 其中， S(x)是吸水强度函数(cm/(cm·d))， (h)是水分胁迫函数， b(x)是根系吸水分配（密度）函数，描述根系吸水的空间变异， Tp是潜在蒸腾量(cm/d)，Ta是实际蒸腾量(cm/d)，h为土壤压力水头(cm)，LR为根系层的深度。

25. 吸水强度 (cm/d) =1 r2H 胁迫函数 r2L =0 P3 P2L P2H Popt P0 压力水头(pressure head, cm) 两个临界吸水强度 数据库 根吸水参数

26. 根系吸水分配函数b(x) • 确定根系随深度的分布，b(x)满足： • 其中，LR为根系层厚度, x是深度。本算例中根系层厚度1 m，假设根系分布函数为线性，并有 • , x1 • 容易证明上式满足积分为1的条件。选择Conditions/ Root Distribution，使用Edit condition，划定根系层范围(100 cm)，顶部数值设为0，底部设为2，让程序自动进行线性差值。

27. Conditions Root Distribution 根系分布函数 使用工具条 根系分布

28. 初始条件 • 选择Conditions/ Initial Condition。本算例中初始地下水位位于地下600 cm，假设土壤剖面初始状态是静力平衡态，则模型底部压力水头为-300 cm，地面压力水头为-600 cm。使用Edit Condition，设置顶部和底部压力水头，并让程序自动插值形成初始条件。

29. 加入观察点 • 选择Conditions/ Observation Points，使用Insert工具条添加若干观察点。

30. 土壤剖面数据汇总

31. 运行模型，查看结果

36. Thank You !