沟灌土壤湿润体空间矩特征参数估算模型.pdf

第２８卷 第６期 水 科 学 进 展 Ｖｏｌ．２８，Ｎｏ．６ ２０１７年１１月　 　 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ Ｎｏｖ．，２０１７ ＤＯＩ： １０１４０４２／ ｊｃｎｋｉ３２１３０９２０１７０６００４ 沟灌土壤湿润体空间矩特征参数估算模型 １ ２ １ １ 聂卫波 ，马孝义 ， 聂坤堃 ， 费良军 （１ 西安理工大学水资源研究所， 陕西 西安　 ７１００４８； ２ 西北农林科技大学旱区农业水土工程教育部重点实验室， 陕西 杨凌　 ７１２１００） 摘要： 为进一步增加矩分析理论在研究沟灌土壤湿润体运移特性方面的实用性， 采用数值模拟和理论分析相结合 的方法， 定量分析了土壤初始含水率、 入渗水深、 沟底宽和边坡系数等因素对沟灌自由入渗条件下土壤湿润体空 间矩特征参数的影响， 并建立了其与主要影响因素间的函数关系式。 结果表明， 土壤初始含水率、 入渗水深和沟 底宽对垂直向土壤水分分布的质量重心影响较大； 水平向湿润锋的平均离散度主要受土壤初始含水率、 入渗水深、 沟底宽和边坡系数的影响， 而垂直向对土壤初始含水率和入渗水深变化的敏感性较高； 选取不同土壤质地对所建 模型的可靠性进行了验证， 结果表明所建模型估算不同入渗时刻的空间矩特征参数值与根据ＨＹＤＲＵＳ软件模拟结 果所得计算值具有高度一致性， 所有验证组合条件下两者相对误差绝对值均值均低于８５％， 且无显著性差异， 说 明所建模型估算沟灌土壤湿润体空间矩特征参数具有高度可靠性。 研究结果可为沟灌灌水方案设计和管理提供理 论基础和技术支撑。 关键词： 沟灌； 土壤湿润体； 矩分析； 特征参数 中图分类号： Ｓ２７５３　 　 　 文献标志码： Ａ　 　 　 文章编号： １００１⁃６７９１（２０１７）０６⁃０８２９⁃１０ 沟灌土壤湿润体运移特性是灌水方案设计和管理的重要基础， 但田间入渗过程中湿润体难以直接观测， 使其较难量化表征， 导致了问题研究的复杂性。 为此， 众多学者提出了多种土壤入渗过程中湿润体运移特性 ［１］ 预测方法， 大致可分为３类， 即解析法、 数值模拟法和经验公式法 。 解析法主要以Ｇｒｅｅｎ⁃Ａｍｐｔ模型为基 ［２⁃３］ 础， 将一维土壤水分运动假定为活塞流， 对湿润体运移距离进行预测 ； 而对于二维和三维土壤水分运动， 较难直接采用该方法进行预测。 数值模拟法主要以土壤水分运动方程为基础， 通过求解该方程来模拟土壤湿 润体运移特性， 目前已有较为成熟的商用软件（如ＨＹＤＲＵＳ软件）， 可适用于多种土壤入渗形式， 并已取得 ［４⁃７］ 了较多的研究成果 ； 但该方法求解过程复杂， 且有时还存在计算结果不收敛和水量平衡误差较大等问题， 这在一定程度上限制了土壤入渗特性的数值模拟研究。 经验公式法主要以土壤入渗试验资料为基础， 通过分 析不同因素对土壤湿润体运移特性的影响， 建立估算模型预测湿润体运移距离［１，８⁃１１］ ； 但沟灌入渗影响因素 众多， 采用试验数据拟合所得湿润锋运移距离预测模型， 一般只在特定的试验条件下成立， 较难推广应用。 近年来， 部分学者通过引入矩分析理论对土壤入渗特性进行研究［１２⁃１４］ ， 其优点在于可较为简单和直观地描 述土壤入渗过程中湿润体范围和入渗水量分布， 避免了复杂繁琐的数值计算过程， 以便于进行灌水方案设计 时控制土壤湿润体范围与作物根系分布相匹配。 但以往研究多为定性分析， 对于如何确定不同入渗时刻沟灌 土壤湿润体空间矩特征参数， 即湿润体水分分布的质量重心和空间分布的离散程度值， 未提出相应的估算方 法， 使得该方法的应用存在一定的局限性。 基于此， 本文在已有研究成果的基础上作进一步拓展， 采用数值 模拟和理论分析相结合的方法， 通过定量分析不同土壤质地条件下初始含水率、 入渗水深、 沟底宽和边坡系