非饱和土壤水盐运移二维数值模拟.pdfVIP

非饱和土壤水盐运移二维数值模拟* 1, 2 2 2 1 刘炳成 ，刘伟 ，孔力 ，李庆领 1. 青岛科技大学机电工程学院，青岛（266061 ） 2.华中科技大学能源与动力工程学院，武汉 (430074) E-mail ：bingchengliu@ 摘 要：本文应用多孔介质的流动和热质传输理论，对土壤水盐运移进行了机理分析，建立 了描述分层土壤中热、湿、盐迁移的数学模型并进行了二维数值模拟，得到了兴隆粘壤土中 温度场、盐含量场和湿度场的动态分布规律，探讨了环境条件对土壤盐分迁移的影响。 关键词：土壤；盐渍化；干饱和层；水盐运移 土壤盐渍化是一种有害溶性盐含量过多导致的土壤中毒现象，已经成为当今世界灌溉农 业中的一个严重问题。全球有各种盐渍土约9.5 亿公顷，占全球陆地面积的10%，广泛分布 于100 多个国家和地区。我国各类盐渍土约有8180 万公顷，其中次生盐渍化土壤占30% 以 上，还有1733 万公顷土壤存在潜在盐渍化。 土壤水盐运移的理论研究始于Darcy 定律，Richards 首先将其用于非饱和土壤水的流动, 并推导出相应的土壤水运动的定量研究[11,12] 。Cardner 和Bresler 对土壤与溶质间的相互作用 进行研究，并总结了考虑作物吸附等因素时非饱和对水盐运动的规律[13] 。可动—不可动两 区模型[22]考虑了溶质运移过程中的化学和物理非平衡性，改善了对溶质迁移特征的描述。 Skopp[23]在两区模型的基础上提出了两流区模型。此外，Hutson[24]和Gurner[25]先后提出了多 流区模型，趋于更加细致的描述土壤孔隙流速分布。80 年代后期张蔚榛[27]等将水盐平衡理 论和数值模拟方法应用于区域水盐预测预报，使我国的农田水盐运动研究进入了一个新的研 究阶段。本文应用多孔介质的流动和热质传输理论，考虑干饱和土壤层的影响，对土壤中热、 湿、盐迁移进行了二维数值模拟，探讨了土壤中温度场、盐含量场和湿度场的动态分布规律。 1 圆柱坐标下土壤水盐运移二维非稳态数学模型 1.1 数学方程 （1）非饱和区 连续性方程 ∂(ρ ε ) ∂ 1 ∂ l l 液相： + (ρ ε u ) + (ρ ε v ) −m (1) l l l l l l ∂t ∂x r ∂r ∂(ρ ε ) ∂ 1 ∂ v g + [ ( + )] + [ ( )] 水蒸汽： ρ ε u u ρ ε v +v m (2) ∂t ∂x v g g v r ∂r v g g v ∂(ρ ε ) ∂ 1 ∂ g g 气相混合物： + (ρ ε u ) + (ρ ε v ) m (3) ∂t ∂x g g g r ∂r g g g 盐分对流弥散方程（考虑不动水体和吸附或解吸作用）： ∂(εlm sm ) ∂(εlim sim ) ∂ ∂sm 1 ∂