冻融期地下水位变化情况下土壤水分运动的初步研究农业工程学报.PDFVIP

第 15 卷　第 2 期 农 业 工 程 学 报 V o l. 15　N o. 2 　　　 　　　　　　　 　　　　　　　 64 1999 年 6 月 T ransactions of the CSA E Jun. 1999 冻融期地下水位变化情况下土壤水分运动的初步研究 ① 尚松浩 　雷志栋　杨诗秀 王　义　赵冬梅 ( ) 　　 ( ) 清华大学 内蒙古自治区水利科学研究院 摘　要　建立了地下水浅埋条件下土壤冻融期水热耦合迁移模型, 并采用混合型R ichards 方程对 冻融过程中的土壤饱和—非饱和水分流动问题进行了初步研究。模拟结果与试验资料的对比表明 所采用的模型是可靠的。根据模拟结果, 冻结过程中潜水位最大降深与初始潜水埋深间的关系可用 反L ogistic 公式表示。对于内蒙古河套灌区, 考虑融化期春小麦需水要求的土壤封冻前的适宜地下 水位约为 1. 5 m～ 2. 0 m 。 关键词　土壤冻融　水热迁移　土壤饱和—非饱和流动　适宜地下水位　数值模拟 地下水位较浅时, 非饱和土壤水分与地下水是相互联系的。一般情况下, 地下水的补给与 ( ) [ 1 ] 消耗及相应潜水位的升降主要由降水 灌溉 、蒸散发所引起 。而在冻融期, 土壤冻结引起地 下水向上补给土壤水, 地下水位下降; 土壤融化引起土壤水下渗补给地下水, 地下水位升高。在 地下水位变化情况下需考虑冻融过程中的土壤饱和—非饱和流动。 对冻融期土壤水热迁移已进行了较多的研究[ 2～ 5 ] 。本文采用混合型R ichards 方程, 并结合 冻土热迁移方程对冻融过程中土壤饱和—非饱和流动进行了数值模拟, 分析了潜水位的变化 规律。在此基础上, 结合内蒙河套灌区的农业生产实际, 分析了冻前适宜潜水埋深问题。 1　地下水浅埋条件下冻融过程中土壤水热耦合迁移模型及其求解 考虑均质各向同性的垂直土柱的一维冻结问题。在忽略水汽迁移、热量对流作用的情况 ( 下, 非饱和土壤水分方程可以用 方程来描述。 方程包括 型 以含水率为自 R ichards R ichards ) ( ) ( ) 变量 、 型 以基质势为自变量 和混合型 方程同时包含含水率、基质势 等三种类型。在研究 h 土壤饱和—非饱和流动问题时, 第一种类型不适用, 第二种类型在数值计算中不容易满足水量 [ 6 ] 平衡关系 。本文采用混合型R ichards 方程来描述冻融过程中土壤水分运动 u h i i = [K ( ) - K ( ) ] - ( 1)