孔隙网络模型在土壤水文学中的应用研究进展-水科学进展.PDF

第 18 卷 第 6 期 水 科 学 进 展 Vol18 ,No6 2007 年 11 月 　　 ADVANCES IN WATER SCIENCE 　 Nov. , 2007 　 孔隙网络模型在土壤水文学中的应用研究进展 吕　菲 , 刘建立 ( 中国科学院南京土壤研究所 , 江苏 南京 　210008) 摘要 : 非饱和土壤的水力性质是采用模型定量模拟水分和溶质在非饱和带中运动的最重要的物理参数 , 可以用网 络模型来预测 。孔隙网络模型的主要优点在于可以对发生在土壤孔隙尺度上的物理 、化学过程进行直观的表达和 模拟 。目前 , 国内外研究者已在样品图像获取 、图像分析 、模型建立以及求解等方面取得了一定的进展 。在对其进 行综合评述的基础上 , 指出了已有网络模型在建立过程中存在的不足之处 , 并进一步明确了今后的研究方向。 关 　键 　词 : 孔隙尺度 ; 网络模型 ; 土壤水力性质 ; 土壤水文学 ( ) 中图分类号 : S1527 ; G35311 　　　文献标识码 : A 　　　文章编号 : 2007 非饱和土壤的水力性质包括土壤水分特征曲线和水力传导率函数 , 是采用模型定量模拟水分和溶质在非饱 和带中运动的最重要的物理参数 , 它与孔隙的几何形态和拓扑性质有着直接的联系[ 1] 。孔隙结构特征与有效水 力性质之间的关系可以用网络模型来模拟 。孔隙网络模型的主要优点在于可以对发生在土壤孔隙尺度上的物 理 、化学过程进行直观的表达和模拟 , 并能直接根据孔隙空间的拓扑性质预测土壤水力性质 。 早期的网络模型将土壤孔隙理想化成圆柱状的毛细管束 , 并通过简单的水分运动模拟来预测连续体尺度上 的介质性质 , 如非饱和水力传导率函数 。但是 , 这些模型均未考虑孔隙之间的相互连通状况 , 因此并不能代表 实际的土壤介质孔隙空间。Fatt 最早提出了以二维或三维网络模型来表征孔隙的结构特征[2～4 ] , 之后这一方法 在石油工程等领域得到了进一步的发展[5 ,6 ] 。近些年来 , 该方法 已逐渐应用到土壤物理学 、水文学等领域 中[7～9 ] , 尤其是多孔介质中的单相流与多相流研究 , 如相对渗透性[8～10 ] ; 两相系统中孔隙结构对相对渗透性 及毛细管压力滞后性的影响[ 11] ; 介质渗透性的预测[ 12 ] ; 毛细管压力 、饱和度及界面面积之间的关系[ 13 ] ; 释水 与吸湿过程[ 14 ,15 ] ; 相的分布 、界面面积及质量传递[ 16 ] 等 。 ( ) 建立孔隙网络模型最关键的步骤是如何识别土壤孔隙结构 孔隙形态与连通状况 并将之定量反映到网络模 型中。由于实际土壤在孔隙形态学上的复杂性 , 为了便于网络模型的建立与求解 , 孔隙体和孔喉通常被简化成 简单的几何形状 : 孔隙体多以球形或立方体表示 , 孔喉则以圆柱形或其他具有代表性的非环形管道表示 。网络 模型可以是二维或者三维的 , 但是二维模型在反映实际介质三维