水利学报.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 蓄水坑灌条件下土壤水分运动的数值模拟 周青云 孙西欢 康绍忠 中国农业大学中国农业水问题研究中心 北京 太原理工大学水利科学与工程学院 山西太原 摘要 本文以非饱和土壤水分运动理论为基础 根据蓄水坑灌法土壤水分入渗的特点 建立了蓄水坑灌条件土壤水 分运动的数学模型 并用有限差分法进行求解 同时用试验资料对模型进行验证 并用模型计算结果分析了蓄水 坑灌条件下土壤初始含水率对土壤水分分布的影响 能为蓄水坑灌法的合理设计提供了一定的参考 关键词 蓄水坑灌法 土壤水分运动 数值模拟 有限差分法 中图分类号 文献标识码 蓄水坑灌法是一种将解决干旱问题和水土流失问题结合在一起适合于山丘区果林灌溉的新方 法 蓄水坑灌法有别于传统的地面灌溉 它主要是根据果树栽培的特点和我国北方地区干旱缺水的 实际情况 依据果树的灌水定额 环果树树冠投影的外围挖若干个深为 的坑 对果树的中深根 层进行灌水 不但供水快而且可以引导果树根系向下深扎从而提高其抗旱能力 目前对蓄水坑灌法的 研究尚处于初级阶段 研究蓄水坑灌法的土壤水分运动规律是进一步认识该灌水技术的基础 随着计 算机技术的发展 数值计算越来越多地应用于农田灌水技术的研究中 数值模拟能够更快 更准确地提 供土壤水分运动的相关信息 是合理确定灌水技术要素的有效方法 本文以非饱和土壤水分运动理论为基础 针对蓄水坑灌法土壤水分入渗的特点 建立土壤水分运动 的数学模型 并对模拟结果进行分析验证 数学模型 单株树蓄水坑布置如图 所示 假定土壤均质且各向同性 因此 在图 中 可视为零通量 面 零通量面两侧没有水量交换 在入渗过程中 水分由水势高处向水势低处运动 即水分由蓄水坑向外 扩展 在入渗初期 相临蓄水坑的湿润土体尚未连在一起 此时的入渗属于自由入渗阶段 当相临蓄水坑 的湿润土体相交时 湿润锋的推进受到相邻湿润体含水率分布状况的影响 属于干扰入渗阶段 在干扰 入渗阶段 相临蓄水坑的湿润土体相交处 即 处两侧水势相等 没有水分运动 可视为零通量面 由 于蓄水坑排列及入渗的对称性 可选择具有代表性的单元进行模拟计算 即通过模拟图 中阴影部分的 土体来研究蓄水坑灌法的入渗特性 图 为蓄水多坑入渗时的计算区域示意图 其中 半圆柱 为蓄水坑 为坑 内水面 为树干位置 为地面 为计算土体的下边界 为零通量面 距离 为两 个蓄水坑中心距离的一半 水分沿蓄水坑壁向周围土体入渗 基本方程 由于蓄水坑边壁区域在入渗初期可能会出现正的压力 因此本模型选用以负压水头 为因变量的非饱和土壤水分运动基本方程 假定土壤均质且各向同性 入渗水流为连续介质且不可压 收稿日期 基金项目 国家自然科学基金项目 作者简介 周青云 女 山西河津人 博士生 主要从事节水灌溉研究 缩 在土壤水分运动过程中土壤骨架不变形 其基本方程在柱坐标系下可描述为 式中 规定 向下为正 为半径 为角度 为负压水头 为比水容量 ? 为 非饱和土壤的导水率 ? 果树 蓄水坑 环状沟 田间输水沟 图 土壤入渗计算区域 图 单株树蓄水坑布置 定解条件 初始条件 式中 为土壤的初始负压水头 为土壤的初始含水率 ? 为土壤水分特征曲线的 经验常数 边界条件 由于计算区域较大 故 边界和 边界可认为水分无法到达 则 边界