最小压力水头值.ppt

软件简介 HYDRUS-1D是美国农业部盐渍土实验室开发的模拟非饱和土壤中的水、热、溶质运移的软件，它在模拟土壤中水分运动、盐分、污染物(如农药)和养分(如土壤氮素)运移方面得到广泛应用。软件分为1维、2维、3维三种，分别命名维Hydrus-1D、Hydrus-2D、Hydrus-3D。由于非饱和土壤水主要是1维垂向运动的形式，因此Hydrus-1D的应用非常广泛。 问题描述 以天津一块农田为背景，对2012年冬小麦收割期农田尺度下土壤水分运动进行数值分析 土壤剖面从地面算起的深度： 300cm 土壤分层：分两层 准备模拟哪个时间段水分变化：0~30d 根系分布：1m 初始地下水位：地面以下6m 主要过程 水分运移方程 将坐标原点选在地面，取z轴向下为正，则一维饱和—非饱和带水分运移基本方程为Richards方程： 式中：θ为含水率；k（θ）为水力传导度：h为负压水头；S 为植物根系吸水量，对裸露区为0。 几何信息 土壤分区长度 [L] 分区水量 [L] 流入或流出量 [LT-1] 平均压力水头[L] 实际表面通量 [LT-1] (入渗/蒸发: -/+) 剖面底实际通量[LT-1] (流入/出: +/-) 整个区域的水质量平衡绝对误差[L] 整个区域的水质量平衡相对误差 [%] 输出信息 迭代参数 上、下边界条件 消光系数 利用比尔定律计算消光系数: 植被蒸腾量和土壤蒸发量分开处理，HYDRUS-1D推荐使用一个经验公式来把潜在蒸散量分割为蒸腾潜力和土壤蒸发潜力： 其中ETp为潜在蒸散量（cm/d），Tp为潜在蒸腾量(cm/d)，Ep为土壤潜在蒸发量(cm/d)，LAI是叶面积指数，k为消光系数。 上边界通量 随时间变化边界条件 即地面土壤达到最干燥状态时的压力水头。理论上，土壤十分干燥时，吸力很大，而液态孔隙水的压强很小，与空气湿度保持平衡关系，有 hA为最小压力水头，Hr为空气绝对湿度，RT/Mg为空气的摩尔气体常数。公式需要近地面空气湿度，一般可取饱和水汽湿度，因2 cm以下土壤空气湿度往往饱和，随温度发生变化，可据近地面气温变化推算地表土壤空气湿度，再换算成压力水头。HYDRUS中需输入hA的绝对值，缺省值为 hCritA=|hA|=106 cm 此数值只对土壤蒸发起作用。建议：hCritA 对应的土壤含水量应至少比残余含水量大0.005，在模拟根系吸水时，hA 还应低于P3。否则，当根系吸水的临界压力水头和地面蒸发的最小压力水头满足时,会导致回流现象，不合理。 时间信息 其中Tp(t)是瞬时潜在蒸腾量，t为时间。模型假设早上6点之前以及晚上18点-24点的蒸腾量总和只占全天蒸腾量的1%。 HYDRUS-1D 可以使用一个经验公式来处理每天24小时潜在蒸腾量的变化，设某天的潜在蒸腾量为 （cm/d），则 土壤水力模型 水力参数 神经网络预测模块 根吸水 根吸水模型 其中， S(x)是吸水强度函数(cm/(cm·d))， ?(h)是水分胁迫函数， b(x)是根系吸水分配（密度）函数，描述根系吸水的空间变异， Tp是潜在蒸腾量(cm/d)，Ta是实际蒸腾量(cm/d)，h为土壤压力水头(cm)，LR为根系层的深度。 根吸水参数 根系吸水分配函数b(x) 确定根系随深度的分布，b(x)满足： 其中，LR为根系层厚度, x是深度。本算例中根系层厚度1 m，假设根系分布函数为线性，并有 , x?1 容易证明上式满足积分为1的条件。选择Conditions/ Root Distribution，使用Edit condition，划定根系层范围(100 cm)，顶部数值设为0，底部设为2，让程序自动进行线性差值。 根系分布 初始条件 选择Conditions/ Initial Condition。本算例中初始地下水位位于地下600 cm，假设土壤剖面初始状态是静力平衡态，则模型底部压力水头为-300 cm，地面压力水头为-600 cm。使用Edit Condition，设置顶部和底部压力水头，并让程序自动插值形成初始条件。 加入观察点 选择Conditions/ Observation Points，使用Insert工具条添加若干观察点。 土壤剖面数据汇总 运行模型，查看结果 P0：Value of the pressure head below which roots start to extract water from the soil. Popt：Value of the pressur