全 文 ：收稿日期 : 2004205230
基金项目 : 教育部优秀青年教师资助项目、人事部留学回国人员科技择优项目重点类课题和国家“十五”西部攻关
专项 (2001BA901A34)及北京林业大学研究生院 2003 年度研究生培养基金项目联合资助
作者简介 : 张强 (1979 —) ,男 ,辽宁人 ,硕士 ,主要研究方向 :土壤水.3 责任作者 : 孙向阳 (1965 —) ,男 ,山东人 ,博士 ,教授 ,主要研究方向 :土壤生态.
林业科学研究　2004 ,17 (增刊) :105～108
Forest Research
　　文章编号 :100121498 (2004)增刊20105204
土壤水分研究进展
张　强1 , 孙向阳 3 , 张广才2
(11 北京林业大学水土保持学院 ,教育部水土保持与荒漠化防治重点实验室 ,北京　100083 ;
21 宁夏林业研究所 ,宁夏 银川　750004)
摘要 :论述了国内外在土壤水分入渗方程模拟、土壤水分运动参数和 SPAC 理论等方面的研究进展。
文中指出沙地水分研究所采用的主要观点由形态学观点转到了能态学观点 ;土壤水分运动参数的测
定方面做了大量的工作 ,取得了较好的成绩 ;随着计算机的迅速发展 ,数值模拟得以广泛的发展。
关键词 :土壤水分运动参数 ;SPAC;数值模拟
中图分类号 :S15217 　　　文献标识码 :A
我国沙漠面积广大 ,约有 100 万 km2 ,而干旱缺水的沙漠生态环境对多数植物的生长造成
威胁 ,也是形成沙区植被稀疏、种类单纯的根本原因。沙地作为一种裸露的生态系统 ,特别是
干旱地区沙地 ,其制约因素中最重要的因素是沙地的水分分布状况及其量值 ,沙地水分状况及
其平衡也是沙区开发利用最重要的依据。随着测试手段和计算机应用的发展 ,以及学科间的
相互渗透 ,土壤水分的研究由经验到理论、由定性到定量发生了质的转变。
1 　土壤水分研究观点的转变
长期以来形态学观点是沙地水分研究所采用的主要观点。这种观点指导人们研究分析了
沙土的水分物理常数 ;确定了沙地水分运动的简单过程 ,加深了对于沙地水分的认识。然而 ,
形态学观点过分的偏重于静态地看待沙地水分运动 ,研究多停留在水分总量的平衡上面 ,将连
续动态变化进行人为地分割 ,缺乏与动力环境的有机结合 ,故难于从本质机理上揭示沙地水分
动力学规律 ,致使沙地水分研究难有较大的突破。
20 世纪出现了土壤水分研究的能态观点 ,30 年代以来 ,特别是 50 年代以后 ,土壤水能态
研究以及在此基础上用数学物理方法定量研究水分运动取得了一系列进展 ,并形成相对独立
的学科 ———土壤水动力学[1 ] 。80 年代末期 ,中国科学院沙坡头开放试验研究站分别与清华大
学水利水电科学系、瑞典Lund 大学水资源工程系进行合作 ,开展了沙地水分动力学的基础性
研究[2 ,3 ] ,从而将我国的沙地水分研究推向新的阶段。
2 　国外土壤水分研究现状
Green2Ampt 模型的修正和 Philip 和 Parlange 入渗方程的求解是国外当前土壤水分入渗研
究的两大主要方面[4 ] 。经过修正的 Green2Ampt 模型能较好地说明非均质土壤的降水入渗过
程。但在应用该模型时 ,如何测定湿润锋处土壤水势的问题未能很好地解决。土壤水分运动
方程的求解问题 ,目前最有效的方法是通过数值计算 ,即有限差分法和有限元法。但是利用数
值计算法所需要处理的数据和运算量相当大 ,因此 ,应用计算机来进行模拟势在必行。
土壤水分运动参数 ,包括水分特征曲线、导水率和扩散率。它表征了土壤本身的性质和特
征 ,是应用数学物理方法定量分析研究土壤水运动的必不可少的资料 ,也是进行土壤水分运动
数学模拟必不可少的依据。推求和测定非饱和导水参数的方法一般分为稳态法和瞬态法两
类。由于非饱和导水率田间的广泛变异、测定技术的复杂、耗时和庞大的费用 ,限制了田间水
分运动的数值模拟研究。因此 ,很多研究者采用比较容易获得的土壤水分特征曲线建立模型 ,
求解土壤非饱和导水率。Marshall 根据 Childs Collis2Georg 建立的方程式发展了仅用水分特征
曲线即可求出非饱和导水率的模型。此后 ,Millington 和 Quirk 对该模型进行了修正 ,并得到了
广泛的应用[5 ] 。Brooks 和 Corey 及 Jeppson 以 Burdine 的模型为基础 ,建立了预测非饱和导水率
的简单解析式。后来 ,Mualem 从对土壤水分特征曲线和土壤饱和导水率的认识出发 ,对上述模
型进行扩展 ,得出一预测导水率的新的积分模型 ,只要为土壤水分特征曲线找到合适的表达
式 ,代入 Mualem 模型就可获得简单的解析式。
3 　国内土壤水分研究现状
20 世纪 50 年代中期至 60 年代中期 ,以原苏联 A1A1Pone 为代表的形态水分研究观点和方
法全面、系统地介绍到中国 ,为我国土壤水分研究的发展起到了积极的推动作用。1977 年在
杭州举行的第一次全国土壤物理学术讨论会上 ,土壤水分的能量概念首次被介绍到国内 ,标志
着我国的土壤水分研究已步入一个崭新的阶段[6 ] 。
随后 ,SPAC理论也在 80 年代初被介绍到国内 ,它把土壤2植物2大气看作是一个连续体 ,水
在该系统中运转传输 ,就像链环一样互相衔接 ,并可用统一的能量指标 ———水势来表示各个部
位甚至不同介质之间能量水平的变化和相互关系[7 ] 。陈建耀等[8 ]利用中国科学院禹城综合试
验站内的大型蒸渗仪、波文比、水力蒸发器等仪器 ,获得了大量水平衡因子的试验数据和 SPAC
模型中的有关参数 ,以大型蒸渗仪实测值为基准 ,验证了农田土壤2植物2大气连续体模型的模
拟值 ,并主要就蒸散和潜水蒸发量 ,对实测与模拟值作了比较分析 ,探讨了两者差异的原因。
邵明安等[9 ]以植物根系吸水的人工模拟试验所测得的资料为依据 ,运用水流的电模拟原理 ,定
量分析了 SPAC中水流阻力各分量的大小、变化规律及其相对重要性。康绍忠等[10 ]在对 SPAC
水分传输机理研究的基础上 ,提出了包括根区土壤水分动态模拟、作物根系吸水模拟和蒸发蒸
腾模拟 3 个子系统的 SPAC 水分传输动态模拟模型 ,设计了 SPAC 水分传输模拟的计算机软
件。
土壤水分运动参数 ,包括水分特征曲线、导水率和扩散率。现在国内已普遍利用水分特征
曲线来研究 115 MPa 以下不同吸力的土壤持水性能 ,而特征曲线的斜率即比水容量 ,可表示吸
力变化时吸入或释出的水量。目前我国主要土类大都已具有小于 115 MPa 的持水曲线和土壤
601 林 　业 　科 　学 　研 　究 第 17 卷
基本性质对持水性影响的资料[11～13 ] 。最近人们对小于 011 MPa 的低吸力段研究较多 ,因为它
所吸持的水分活动性最强 ,且大部分是易效水。如庄季屏、王伟[14 ]对东北 5 种耕作土壤小于
011 MPa 的低吸力段持水性能进行了研究 ,还对影响低吸力段持水量的因素、所吸持水分对植
物的有效程度、比水容量及脱水速度进行了研究。数理建模是研究土壤水分运动的重要手段。
杨诗秀等[15 ]用 FORTRAN 语言建立了一个被饱和水一维流动的数值计算模型 ,并作了实验室
验证。该模型用以模拟均质土壤、地下水埋藏很深 ,且在不同的初始条件与边界条件下的水分
运动 ,可用来计算土壤剖面入渗、蒸发、蒸腾和再分配以及当这些现象交替出现时的水分运动
过程。周择福、洪玲霞[16 ]由达西定理和能量守恒原理推导了土壤水分入渗的数学模型 ,水平
土柱法实测了模型中的基本运动参数 :土壤水分扩散率 D (θ) ,推求了土壤水分非饱和导水率
K(θ) ,经过计算机用有限差分法模拟了 6 块不同林地的土壤水分入渗过程 ,实地试验检验了
模拟结果 ,得到了很好的效果。吴长文、王礼先[17 ]采用 4 种常见的入渗模型对北京密云水库
库区的 8 种主要林分的坡面土壤进行了入渗模拟 ,结果 Horton 入渗模型可方便地计算出林地
土壤的有效贮水量 ,Smith2P 模型能较好地描述在降雨条件下的入渗过程。周维博[18 ]通过野
外试验 ,求得了计算裸露土壤蒸发条件下土壤水分运动的上边界条件 ,并用数值模拟的方法 ,
通过实例 ,分析了降雨蒸发条件下层状土壤剖面土壤水分运移和相互转化关系。邵明安等[19 ]
使用积分方法求解了一维水平非饱和土壤水分运动问题 ,根据其解建立了推求非饱和土壤水
分运动参数的简单入渗法 ,用以推求 Van Genuchten 特征曲线模型中的参数α和 n。α和 n 是
根据湿润区的特征长度、吸渗率和土壤的饱和导水率 ( ks) 来确定的 ,而非饱和导水率可由α、n
和 ks 确定。这一新的简单入渗法是基于 Richards 方程和土壤导水特征的闭合型方程。简单入
渗法提供了利用瞬态水流方法来确定土壤水分特征曲线而替代通常的平衡法。目前测定非饱
和导水率应用较广的实验室方法是“稳定流”法 ,在田间原位测定时则用“内排水”法 ;至于扩散
率的测定 ,通常是在实验室内用水平土柱法测得。
由于缺乏对沙地水文过程 ,尤其对于土2气界面、土2根界面的深入理解 ,还难以对沙地水
分动态给以全面的认识。沙地凝结水是土2气相互作用的产物 ,在干旱或极端干旱区有重要的
研究价值 ,它对植物生存和局部地区短历时生态环境有重要意义[20 ] 。我国也开展了沙地凝结
水的研究 ,但仍有待于进一步研究。
参考文献 :
[1 ] 雷志栋 ,杨诗秀 ,谢森传. 土壤水动力学[M] . 北京 :清华大学出版社 ,1988
[2 ] 杨诗秀 ,雷志栋 ,段新杰 ,等. 应用土壤水动力学原理研究沙区土壤水分动态特征[A] . 见 :刘家琼 :中国科学院沙坡头沙
漠试验研究站年报 1991 —1992[ R] . 兰州 :甘肃科学技术出版社 ,1993
[3 ] 刘元波 ,杨诗秀. 积水条件下沙地入渗湿润锋速度变化规律[A] . 见 :崔书红. 沙漠生态系统[M] . 兰州 :甘肃科学技术出版
社 ,1994
[4 ] 马履一. 国内外土壤水分研究现状与进展[J ] . 世界林业研究 ,1997 (5) :26～31
[ 5 ] Bruce R R. Hydraulic conductivity evaluation of the soil profile from soil water retention relations[J ] . Soil Sci Am Proc ,1972 ,36 :555～
561
[6 ] 朱祖祥. 土壤水分的能量概念及其意义[J ] . 土壤学进展 ,1979 (1) :1～2
[7 ] 庄季屏. 四十年来的中国土壤水分研究[J ] . 土壤学报 ,1989 ,26 (3) :241～247
[8 ] 陈建耀 ,刘昌明 ,吴凯. 利用大型蒸渗仪模拟土壤2植物2大气连续体水分蒸散[J ] . 应用生态学报 ,1999 ,10 (1) :45～48
[9 ] 邵明安 ,陈志雄. SPAC中的水分运动[J ] . 西北水土保持研究所集刊 ,1991 (13) :3～12
701增刊 张强等 :土壤水分研究进展
[10 ] 康绍忠 ,刘晓明 ,高新科 ,等. 土壤2植物2大气连续体水分传输的计算机模拟[J ] . 水利学报 ,1992 (3) :1～12
[11 ] 陈志雄 ,汪仁真. 中国几种主要土壤的持水性质[J ] . 土壤学报 ,1979 ,14 (3) :277～281
[12 ] 杨金楼 ,朱连龙 ,朱济成. 上海地区土壤持水性的研究[J ] . 土壤学报 ,1982 ,19 (4) :331～341
[13 ] 张景略 ,苗付山. 黄泛平原不同质地土壤的持水特性[J ] . 土壤学报 ,1985 ,22 (4) 350～355
[14 ] 庄季屏 ,王伟. 土壤低吸力段持水性能及其与早期土壤干旱的关系研究[J ] . 土壤学报 ,1986 ,23 (4) : 306～313
[15 ] 杨诗秀 ,雷志栋 ,谢森传. 均质土壤一维非饱和流动通用程序[J ] . 土壤学报 ,1985 (1) :24～34
[16 ] 周择福 ,洪玲霞. 不同林地土壤水分入渗和入渗模拟的研究[J ] . 林业科学 ,1997 ,33 (1) :9～17
[17 ] 吴长文 ,王礼先. 林地土壤的入渗及其模拟分析[J ] . 水土保持研究 ,1995 ,2 (1) :71～75
[18 ] 周维博. 降雨入渗和蒸发条件下野外层状土壤水分运动的数值模拟[J ] . 水利学报 ,1991 (9) :32～36
[19 ] 邵明安 ,王全九. 推求土壤水分运动参数的简单入渗法 (理论分析) [J ] . 土壤学报 ,2000 ,37 (1) :1～7
[20 ] 刘元波 ,陈荷生 ,高前兆. 沙地水分动力学研究新视角[J ] . 中国沙漠 ,1997 ,17 (1) :95～97
Advance in Soil Water Research
ZHANG Qiang1 , SUN Xiang2yang1 , ZHANG Guang2cai2
(11Beijing Forestry University ,Beijing 　100083 ,China ;
21Ningxia Forestry Institute ,Yinchuan 　750004 ,Ningxia ,China)
Abstract :The research at home and abroad on the simulation of soil water infiltration equation soil water move2
ment parameters and the theory of SPAC were reviewed. The authors indicated that the main standpoint had been
moved to energy from morphology ,and much effect had been made on the measurement of soil water movement
parameters and great achievements had been gained. With the rapid development of computer technique ,exten2
sive progress in quantitative simulation had been made.
Key words :soil water movement parameter ;SPAC;quantitative simulation
801 林 　业 　科 　学 　研 　究 第 17 卷