内蒙古河套灌区土壤冻融过程水热耦合迁移模拟研究.docxVIP

内蒙古河套灌区土壤冻融过程水热耦合迁移模拟研究 1 土壤水热耦合迁移 内蒙古河套灌区是我国第三个灌区，位于河套平原中部。该灌溉的特点是高热资源和干旱。由于独特的黄河水灌溉，该灌溉的热量资源潜力得到了充分利用。这是内蒙古的一个重要粮食基地。河套地区是季节性寒冷地区，土壤在一年半以上后会发生变化。由于各种原因，该地区灌水量大，排水停滞，地下水位高。在土壤冻结过程中，地下水被转移到冻结区，土壤盐分被转移到冻土层。随着春季土壤蒸发，地表土壤盐分积累，甚至水分盐渍，这对春季小麦的栽培和发芽期的生长非常有害。因此，深入研究固河地区土壤冻结过程中的水分转移规律对农业生产至关重要。 土壤冻融过程中的水热迁移是冻土物理学研究的基本内容.土壤冻结过程中的热量传输、水分迁移相互作用、相互影响.在研究土壤冻结过程中, 只有综合考虑水、热的耦合作用, 才能准确反映其运动的客观物理机制, 从而掌握其运动规律.目前, 用于描述冻结过程中土壤水热耦合迁移的模型主要有两类:基于多孔介质中液态水分的粘性流动、热平衡原理的Harland模型 (Harland, 1973) 和应用不可逆热力学原理描述土壤水热通量的热力学模型 (Kunget al., 1986) . Wang Yi (1993) 利用地中渗透仪对此问题进行了现场试验研究.本文以内蒙古河套灌区越冬期田间土壤墒情变化为背景, 应用在Harland模型基础上建立的水热耦合迁移模型 (尚松浩等, 1997) , 对河套灌区地下水浅埋条件下整个冻融过程中的土壤水热迁移进行了模拟.在模拟的基础上, 分析了越冬期自然条件下土壤水热迁移规律. 2 土壤温度与热耦合 考虑均质各向同性的垂直土柱的一维冻结问题.在忽略水汽迁移、热量对流作用的情况下, 非饱和土壤水分方程为 (雷志栋等, 1988) : ?θu?t=??z(D?θu?z?K)?ρsρw?θs?t(1) 式中:θu, θs分别为土壤未冻水、冰的体积含量;t为时间;z为空间坐标 (以垂直向下为正) ;D(θu) ,K(θu) 分别为土壤非饱和水分扩散率及导水率;ρs, ρw分别为冰、水密度. 在冻土中, 未冻水含量与负温保持动态平衡, 这一关系表示了冻土中水、热运动间的相互联系: θu=θm(T)(2) 式中:θm(T) 为相应土壤负温条件下可能的最大未冻水含量. 土壤温度的变化用冻土水热耦合方程来描述 (尚松浩等, 1997) : Ce?T?t=??z(λe?T?z)?Ue?T?t(3) 式中:T为土壤温度;Ce, λe,Ue分别为冻结土壤的等效体积热容量、热导率、热对流速度, 其值分别为: Ce=Cv+Cl(4a)λe=λ+D(θu)?Cl(4b)Ue=Cl(dK(θ)/dθ)(4c)Cl=Lfρw(dθm/dT)(4d) 式中:Cv,Cl, λ分别为土壤体积热容量、相变热容量、热导率;Lf为融化潜热. 在非冻土中,Cl=0, 方程 (3) 即为非冻土的热传导方程.因此方程 (3) 既适用于土壤冻结区, 也适用于土壤未冻结区. 对上述模型的求解采用全隐式有限差线性化迭代方法 (尚松浩等, 1997) . 3 土壤水分变化过程 田间试验于1990年11月15日至1991年4月30日 (共166 d) 在内蒙古水利科学研究所巴音实验场的地中渗透仪中进行, 该实验场位于内蒙古自治区巴彦淖尔盟临河市小召乡境内.地中渗透仪为圆形, 直径113cm (面积1 m2) , 深度约为3.2 m, 下端封闭.试验过程中, 通过不断补水或吸水来保持地下水位的恒定, 潜水埋深分别为1.5 m, 2.0 m和2.5 m.试验用土为扰动土, 土样取用河套灌区具有代表性的土壤砂质壤土.试验过程中的观测项目主要包括土壤含水量、温度.土壤水分的动态变化用中子仪来测定.表层土壤的温度用地温计测量, 深层地温用热敏电阻温度计测量.另外, 试验期间同时还进行了气温观测. 在土壤冻结期, 土壤蒸发比较微弱, 在计算过程中忽略土壤蒸发, 将水分运动上边界视为零通量边界.水分运动下边界取为地下水位处, 取其含水率为略小于饱和值的0.435.土壤上边界温度取地表温度的回归曲线. 在上述边界条件下模拟了潜水埋深为2 m情况下土壤冻融过程, 图1, 2为模拟结果与试验观测值的比较.图1表示潜水位以上土壤贮水量的动态变化过程, 土壤贮水量的计算值与实测值吻合较好.在土壤融化期计算的贮水量偏高, 这主要是由于在计算过程中忽略土壤蒸发造成的. 土壤贮水量在冻结期增加而在融化期减小是由土壤冻融过程中水分迁移的特殊规律造成的.在冻结期温度梯度所引起的基质势梯度作用下, 地下水不断向上层补给并向冻结区迁移, 使得地下水位以上土壤贮水量不断增加.从3月上、中旬起, 土壤从冻土层下部开始向上融化, 融化后的水分又逐渐回