包气带水文地质学-第五章-土壤—植物—大气系统中的水流运动课件.ppt

《包气带水文地质学-第五章-土壤—植物—大气系统中的水流运动课件.ppt》由会员分享，可在线阅读，更多相关《包气带水文地质学-第五章-土壤—植物—大气系统中的水流运动课件.ppt（39页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、SPAC1第五章第五章 土壤土壤植物植物大气系统中的大气系统中的水流运动水流运动Chapter 5 Water Moves in SPAC5.1 概述概述5.2 植物体中水流植物体中水流5.3 腾发量的估算腾发量的估算5.4 植物根系吸水规律植物根系吸水规律25.1 概述概述SPAC系统的概念系统的概念l长期以来，土壤物理学、植物生理学及气象学等学科，在各自领域里，对水长期以来，土壤物理学、植物生理学及气象学等学科，在各自领域里，对水分运动规律分别进行研究，把一个本来统一的整体，人为地隔裂为一个个孤分运动规律分别进行研究，把一个本来统一的整体，人为地隔裂为一个个孤立的部分进行研究。立的部分进行

2、研究。l田间水分循环中的土壤、植物与大气是一个物理上统一的动态系统田间水分循环中的土壤、植物与大气是一个物理上统一的动态系统，各种过各种过程相互关连，程相互关连，Philip，J.R.(1966)把这个系统称为把这个系统称为土壤土壤植物植物大气连续体大气连续体(Soil-Plant-Atmosphere Continuum)，简称简称SPAC系统系统。l在这个系统中，水总是从能量高处向能量低处运动，在这个系统中，水总是从能量高处向能量低处运动，“水势水势”理论在土壤、理论在土壤、植物和大气中普遍通用。植物和大气中普遍通用。lSPAC系统的提出是土壤系统的提出是土壤植物植物水分关系研究方面的一次

3、重要突破，它水分关系研究方面的一次重要突破，它不仅统一了不仅统一了SPAC中的能量关系，为分析和研究系统中的水分运输、能量转化中的能量关系，为分析和研究系统中的水分运输、能量转化的动态过程提供了方便，而且将进一步促进农田水利与土壤物理、自然地理、的动态过程提供了方便，而且将进一步促进农田水利与土壤物理、自然地理、植物生理、农业气象、农田生态、农业水文等边缘学科的相互渗透。植物生理、农业气象、农田生态、农业水文等边缘学科的相互渗透。水势决定水势决定SPAC中水的运动方向中水的运动方向lSPAC中的水分由水势高处向水势低中的水分由水势高处向水势低处运动处运动。lSPAC中的水流中的水流土壤中的水分

4、运动土壤中的水分运动（土面蒸发与入渗）。根系吸水过程根系吸水过程。决定于土壤的特性和根系吸水性能。水在植物体内运移、输送过程水在植物体内运移、输送过程。决定于植物体的输水能力。植物体内水分向大气散发过程植物体内水分向大气散发过程。决定于叶面水汽压与大气水汽压差。4 SPAC系统水流运动复杂系统水流运动复杂lSPAC系统中水流运动是一个十分复杂的过程系统中水流运动是一个十分复杂的过程。液态水的运动土壤和植物体中的水流；气态水的扩散干土层和大气中的水汽扩散；水的相变土面蒸发和叶面蒸腾过程中液态水到气态水的相变。l目前，对目前，对SPAC系统的水流运动研究尚处于初始阶段，许系统的水流运动研究尚处于初

5、始阶段，许多理论问题没有解决，有待深入研究。多理论问题没有解决，有待深入研究。5植物体内水势植物体内水势l溶质势（溶质势（s s）植物根系吸水的机理是因液泡膜(也称原生质膜)为半透膜，液泡膜内的细胞液是含有多种成分溶质的复杂溶液。由于溶质的存在，与纯水相比，降低了自由能，其能量相对差值为溶质势溶质势，也称，也称渗透压势渗透压势。在势能梯度作用下，就产生吸水或失水的过程，s s 00。细胞溶质势的高低，主要取决于液泡内溶液的浓度，温带生长的大多数作物细胞的溶质势为-10X105-20X105Pa之间。5.2 植物体中水流植物体中水流6l压力势（压力势（p p）当水进入细胞，细胞因吸水而膨胀，此时

6、，原生质体向细当水进入细胞，细胞因吸水而膨胀，此时，原生质体向细胞壁产生压力，称为胞壁产生压力，称为膨压膨压；相应地，细胞壁也将对内产生；相应地，细胞壁也将对内产生压力，称为压力，称为壁压壁压，膨压和壁压大小相等，方向相反。，膨压和壁压大小相等，方向相反。原生原生质体所受壁压为正压力，质体所受壁压为正压力，该压力引起的水势称为该压力引起的水势称为压力势压力势，为为正值正值。l基质势（基质势（m m）细胞原生质体中的细胞质，是一种亲水胶体，对水分有较细胞原生质体中的细胞质，是一种亲水胶体，对水分有较强的吸附作用，从而降低了其中水分的自由能，由此引起强的吸附作用，从而降低了其中水分的自由能，由此引

7、起的水势叫基质势的水势叫基质势,m m 0 0。l 植物体内细胞的总水势可用下式表示：植物体内细胞的总水势可用下式表示：7l当细胞胶体被水当细胞胶体被水饱和饱和时，时，基质势趋近于零基质势趋近于零，仅，仅-0.1X 105 Pa左右，可忽略不左右，可忽略不计。此时，植物细胞水势可表示为：计。此时，植物细胞水势可表示为：l干燥的种子，细胞未形成液泡，则溶质势和压力势均为零，仅有基质势。干燥的种子，细胞未形成液泡，则溶质势和压力势均为零，仅有基质势。l若细胞继续吸水，若细胞继续吸水，压力势等于溶质势时，细胞水势等于零压力势等于溶质势时，细胞水势等于零。则细胞失去吸水。则细胞失去吸水能力。能力。l上

8、述细胞中水势的变化是由各分势变化引起的，而各分势的变化并非独立的，上述细胞中水势的变化是由各分势变化引起的，而各分势的变化并非独立的，而是相互关联的。而是相互关联的。l植物根系内的生理活动使根系从土壤中吸水，并从根部上升的压力称为植物根系内的生理活动使根系从土壤中吸水，并从根部上升的压力称为根压根压。根压根压使植物吸水的现象称为使植物吸水的现象称为主动吸水主动吸水。l植物植物蒸腾作用蒸腾作用引起植物根部从土壤中吸水，补充植物体消耗水分的现象称为引起植物根部从土壤中吸水，补充植物体消耗水分的现象称为被动吸水被动吸水。8 植物体中水流植物体中水流l植物由根系从土壤中吸收水分，极少植物由根系从土壤中

9、吸收水分，极少量用于各种代谢作用，主要消耗于蒸量用于各种代谢作用，主要消耗于蒸腾作用。腾作用。l在充分供水条件下，外界蒸发条件基在充分供水条件下，外界蒸发条件基本不变，可假定流经植物体内的水流本不变，可假定流经植物体内的水流为稳定流。即叶面蒸腾强度等于植物为稳定流。即叶面蒸腾强度等于植物体内输水速度等于植物根部对土壤水体内输水速度等于植物根部对土壤水分的吸收速度，可用下式表示：分的吸收速度，可用下式表示：临界叶水势临界叶水势-30-30barbar9为了便于分析和考虑问题，通常将这个连续体统用一个等为了便于分析和考虑问题，通常将这个连续体统用一个等效电路来代替效电路来代替 105.3 腾发量的

10、估算腾发量的估算腾发量概述腾发量概述l 腾发量腾发量是单位时间、单位面积上土面蒸发和作物蒸腾量之和。是单位时间、单位面积上土面蒸发和作物蒸腾量之和。它是分析计算作物生长条件下土壤水分运动的基本资料。它是分析计算作物生长条件下土壤水分运动的基本资料。l影响腾发量的因素影响腾发量的因素供给能量（液供给能量（液气）气）外界条件外界条件水汽输送水汽输送土壤特性、植物本身特性土壤特性、植物本身特性土壤水分和植物吸水性能则是控制腾发的第三个因素。土壤水分和植物吸水性能则是控制腾发的第三个因素。土壤含水率越高，输水能力越高，可向植物根系充分供水。土壤含水率越高，输水能力越高，可向植物根系充分供水。植物不同生

11、长阶段，蒸腾量不同。植物不同生长阶段，蒸腾量不同。11l腾发量的类型腾发量的类型参考作物潜在腾发量：参考作物潜在腾发量：指地面完全覆盖，高度均匀一致，充分供水条件指地面完全覆盖，高度均匀一致，充分供水条件下，短的绿草地下，短的绿草地(如牧草如牧草)上腾发的水量，主要上腾发的水量，主要受气象条件影响受气象条件影响。FAOFAO修改的参照腾发量的定义修改的参照腾发量的定义为：参照腾发量为一种假想的参照作物冠层为：参照腾发量为一种假想的参照作物冠层的腾发速率，假设作物高度为的腾发速率，假设作物高度为0.12m0.12m，固定的叶面阻力为，固定的叶面阻力为70s/m70s/m，反射率，反射率为为0.2

12、30.23，非常类似于表面开阔，高度一致，生长旺盛，完全遮盖地面而，非常类似于表面开阔，高度一致，生长旺盛，完全遮盖地面而不缺水的绿色草地的蒸腾和蒸发量。不缺水的绿色草地的蒸腾和蒸发量。实际作物潜在腾发量：实际作物潜在腾发量：作物在充分供水作物在充分供水，生长良好条件下的最大腾发量。生长良好条件下的最大腾发量。它除它除受气象条件影响受气象条件影响外，还外，还与作物及生育阶段有关与作物及生育阶段有关。作物实际腾发量：作物实际腾发量：在田间实际条件下作物的蒸发量。它除在田间实际条件下作物的蒸发量。它除受气象条件、受气象条件、作物及其生育阶段影响作物及其生育阶段影响外，还外，还受土壤水分状况的影响受

13、土壤水分状况的影响。12l作物蒸腾和土壤蒸发的区别作物蒸腾和土壤蒸发的区别不同点不同点：作物的蒸腾不仅决定于外界条件、土壤性质和土壤：作物的蒸腾不仅决定于外界条件、土壤性质和土壤含水量，同时也与作物有关。含水量，同时也与作物有关。共同点共同点，在于它，在于它也具有两个也具有两个阶段阶段。在土壤含水率高于临界含水在土壤含水率高于临界含水率时，作物蒸腾也仅与外界率时，作物蒸腾也仅与外界条件有关，而在临界含水率条件有关，而在临界含水率以下时，蒸腾强度同样随水以下时，蒸腾强度同样随水分消耗而不断下降，其数值分消耗而不断下降，其数值与有效含水率与有效含水率(可被作物利用可被作物利用的含水率的含水率)成直

14、线关系。成直线关系。13 潜在腾发量的估算方法潜在腾发量的估算方法l空气动力学法空气动力学法（根据紊流扩散理论计算腾发强度）；（根据紊流扩散理论计算腾发强度）；l能量平衡法能量平衡法（利用能量平衡原理估算腾发量）；（利用能量平衡原理估算腾发量）；l综合法综合法Penman-Monteith公式公式（以上两种方法结合）；（以上两种方法结合）；lFAO-PPP-17修正式修正式l经验公式法经验公式法。14l空气动力学法空气动力学法空气动力学法空气动力学法根据紊流扩散理论计算腾发强度，又称紊流扩根据紊流扩散理论计算腾发强度，又称紊流扩散法。散法。由于近地表的空气处于紊流流态，根据紊流扩散理论：由于近

15、地表的空气处于紊流流态，根据紊流扩散理论：ET单位时间，单位面积上的腾发量单位时间，单位面积上的腾发量g/(cm2s)空气的密度空气的密度(g/cm3)在在20oC、1013 mbar时为时为 1.2X10-3g/cm3；水的分子量与空气的分子量之比水的分子量与空气的分子量之比(为为0.622)；Pa大气压大气压(mbar)；el，e2相应于离地面高度为相应于离地面高度为z1和和z2处的水汽压处的水汽压(mbar)；u1,u2相应于离地面高度为相应于离地面高度为z1和和z2处的平均风速，具速度量纲处的平均风速，具速度量纲(m/s)；K是冯是冯.卡曼卡曼(Von Karman)常数常数(0.41

16、)。15l能量平衡法能量平衡法能量平衡方法能量平衡方法利用能量平衡原理估算腾发量。将腾发看作能利用能量平衡原理估算腾发量。将腾发看作能量消耗过程，通过平衡计算，求出腾发所消耗的能量，再将量消耗过程，通过平衡计算，求出腾发所消耗的能量，再将其折算成水量，用以计算作物田间耗水量。其折算成水量，用以计算作物田间耗水量。计算腾发强度计算腾发强度ET的计算式为的计算式为L汽化潜热汽化潜热(J/g)；鲍文比；鲍文比；Rn 土壤表面的净辐射土壤表面的净辐射J/(cm2d)；G土壤表面的热通量J/(cm2d)；湿度计常数，与气压和温度有关，可近似取值为湿度计常数，与气压和温度有关，可近似取值为0.66mbar

17、/oC；T1，T2 分别为两个不同高度分别为两个不同高度z1、z2处的绝对温度处的绝对温度；e1，e2为两个不同高度为两个不同高度zl、z2处的水汽压处的水汽压(mbar)。16l 综合法综合法彭曼公式彭曼公式彭曼公式是在紊流扩散法与能量平衡法基础上建立彭曼公式是在紊流扩散法与能量平衡法基础上建立的，它的，它是上述两种方法综合结果。由英国是上述两种方法综合结果。由英国Penman(1948)提出，是提出，是目前计算腾发量中应用较为广泛的方法目前计算腾发量中应用较为广泛的方法。假定假定：蒸发面水汽压力为饱和水汽压力，该压力仅是温度的函数，蒸发面水汽压力为饱和水汽压力，该压力仅是温度的函数，求得的

18、腾发量为潜在腾发量求得的腾发量为潜在腾发量ET；蒸发面处的气温蒸发面处的气温t与观测高程与观测高程z2(=2m)处气温相等。处气温相等。17es2观测高度观测高度z2面上的饱和水汽压；面上的饱和水汽压；es0蒸发面上的水汽压；蒸发面上的水汽压；T2,T0分别为观测面上和蒸发面上的绝对温度；分别为观测面上和蒸发面上的绝对温度；湿度计常数湿度计常数(mbar/)；Ea干燥力，蒸发面上温度等于气温时的蒸发量干燥力，蒸发面上温度等于气温时的蒸发量(mm/d).u1,u2z2高度上风速高度上风速(ms)；es1,e2分别表示观测面上饱和水汽压和量测时实际水汽压分别表示观测面上饱和水汽压和量测时实际水汽压

19、(一般以一般以mbar计计)。上式中，湿度计常数与气压有关，上式中，湿度计常数与气压有关，es,是温度的函数，在一个标准大气压是温度的函数，在一个标准大气压(1013mbar)时有关数值如表所示。时有关数值如表所示。18lFAO-PPP-17彭曼修正式彭曼修正式19l 经验公式法经验公式法以太阳辐射以太阳辐射Rs为指标的估算腾发量的公式为指标的估算腾发量的公式 LETpa1 Rs+a2 a1,a2经验系数，随作物和地点而异，由田间试验资料确定。经验系数，随作物和地点而异，由田间试验资料确定。式中式中ETpi第第i个月份的潜在蒸发量个月份的潜在蒸发量(mm/d)；T相应月份的日气温平均值相应月份

20、的日气温平均值()；Pdi该月份日平均白昼小时占全年昼长小时的百分数；该月份日平均白昼小时占全年昼长小时的百分数；Ci经验系数。经验系数。以气温为指标的腾发量公式以气温为指标的腾发量公式 可用时段内平均气温或积温表示，可用时段内平均气温或积温表示，1977年年FAO推荐公式：推荐公式：20以空气湿度为指标估算腾发量的公式以空气湿度为指标估算腾发量的公式ETp(es-e)式中式中 经验系数；经验系数；es蒸发面上饱和水汽压；蒸发面上饱和水汽压；e空气中实际水汽压。空气中实际水汽压。水面蒸发量法水面蒸发量法以水面蒸发量来反映综合气象因素，其一般形式为以水面蒸发量来反映综合气象因素，其一般形式为ET

21、p E0式中式中:E0水面蒸发量；水面蒸发量；经验常数，决定于参考作物和蒸发皿的类型。经验常数，决定于参考作物和蒸发皿的类型。l以上公式经验系数，均需由试验资料确定。应根据当地条件以上公式经验系数，均需由试验资料确定。应根据当地条件选用合适公式，以实测试验资料确定经验系数。选用合适公式，以实测试验资料确定经验系数。21 作物实际腾发量的计算作物实际腾发量的计算作物系数作物系数Kc反映了作物和土壤因素。反映了作物和土壤因素。作作物物因因素素:作作物物种种类类、作作物物生生理理特特性性、发发育育阶阶段段、播播种种日日期期、生生育期长短等。育期长短等。地理位置地理位置:不同地区作物系数变化规律不同，

22、除赤道附近不同地区作物系数变化规律不同，除赤道附近KcKc基本基本为常数外，其他地区作物系数都随作物的生长期而变化。为常数外，其他地区作物系数都随作物的生长期而变化。KcKc值通常参考联合国粮农组织灌溉及排水从书第值通常参考联合国粮农组织灌溉及排水从书第2424卷提出的作卷提出的作物生育四阶段标准划分：初期、发展期、中期、后期物生育四阶段标准划分：初期、发展期、中期、后期。l根据潜在腾发量估算作物实际腾发量根据潜在腾发量估算作物实际腾发量计算出的潜在腾发量，必须根据计算出的潜在腾发量，必须根据作物系数作物系数Kc和和土壤供水系数土壤供水系数K 加以修正，才能求得实际的腾发量，计算公式为：加以修

23、正，才能求得实际的腾发量，计算公式为：ET(t)=ETp Kc K 22部分作物生长阶段天数（部分作物生长阶段天数（d）23部部分分作作物物在在生生长长中中期期和和完完熟熟时时的的Kc值值24华北地区冬小麦作物系数华北地区冬小麦作物系数Kc 地区地区播种播种-分蘖分蘖分蘖分蘖-返青返青返青返青-拨节拨节拨节拨节-抽穗抽穗抽穗抽穗-成熟成熟全期全期产量水平产量水平（kg/亩）亩）河北临西河北临西1.450.880.700.900.650.86300河北望都河北望都0.580.580.981.520.960.88350山西山西0.860.490.851.031.070.86380山东山东1.050

24、.320.761.151.220.92-0.99河南新乡河南新乡1.39-1.210.51-0.920.92-0.690.871.220.97300天津天津0.88300-350北京彰化北京彰化0.670.340.471.471.530.9240025华北地区夏玉米作物系数华北地区夏玉米作物系数Kc地区地区播种播种-拨节拨节 拨节拨节-抽穗抽穗 抽穗抽穗-收获收获全期全期产量水平产量水平(kg/亩亩)河北临西河北临西0.960.310.740.71472河北藁城河北藁城0.650.370.950.81367-465河南新乡河南新乡0.560.661.32-0.940.7740526华北地区棉花

25、作物系数华北地区棉花作物系数Kc 地区地区播种播种-出苗出苗出苗出苗-现蕾现蕾现蕾现蕾-开花开花开花开花-吐絮吐絮吐絮吐絮-收获收获全期全期产量水平产量水平(kg/亩亩)河北临西河北临西0.320.581.011.160.760.77河北藁城河北藁城0.500.500.910.950.440.72河南新乡河南新乡(夏棉夏棉)0.66-0.76河南新乡河南新乡0.130.360.891.431.500.96144(皮棉皮棉)农田腾发量随土壤含水率变化，一般可分为两个阶段，农田腾发量随土壤含水率变化，一般可分为两个阶段，当土壤含水率高于临界含水率时当土壤含水率高于临界含水率时，农田腾发量大小主要取

26、决于，农田腾发量大小主要取决于气象条件和作物生长状况，不随土壤含水率而变化，称为稳定气象条件和作物生长状况，不随土壤含水率而变化，称为稳定腾发阶段。此时，作物实际腾发量与潜在腾发量的关系为：腾发阶段。此时，作物实际腾发量与潜在腾发量的关系为：ET土壤充分供水时，土壤充分供水时，t时段内平均腾发量；时段内平均腾发量；P t时段内降雨时段内降雨P；I t时段内时段内灌水总量；灌水总量；Q t时段内根层以下深层渗漏量；时段内根层以下深层渗漏量；W t时段内，计算深层渗漏量时段内，计算深层渗漏量所取断面以上土体所消耗的水量所取断面以上土体所消耗的水量(可用时段始末土壤含水率剖面求得可用时段始末土壤含水

27、率剖面求得)。此时此时ET为作物根层平均土壤含水率高于临界含水率时的农田为作物根层平均土壤含水率高于临界含水率时的农田实际腾发量，可根据田间试验资料，用水量平衡方法求得：实际腾发量，可根据田间试验资料，用水量平衡方法求得：lKc 的确定方法的确定方法28 根据临西试验站实测资料通过平衡计算求得平均作物根据临西试验站实测资料通过平衡计算求得平均作物系数与相对时间系数与相对时间(作物生长历时与全生育期历时的比值作物生长历时与全生育期历时的比值)关系曲线。关系曲线。29土壤供水修正系数土壤供水修正系数K ，与土壤质地及土壤含水量有关。，与土壤质地及土壤含水量有关。当土壤含水率低于某一临界值时，由于土

28、壤供水能力有限，腾当土壤含水率低于某一临界值时，由于土壤供水能力有限，腾发量受到抑制，农田腾发量将随土壤含水率而变化，此阶段为发量受到抑制，农田腾发量将随土壤含水率而变化，此阶段为非稳定腾发阶段非稳定腾发阶段。M.Jenson公式：公式：Thornthwaite公式：公式：经验公式经验公式30l 根据水面蒸发量估算作物实际蒸腾量根据水面蒸发量估算作物实际蒸腾量水面蒸发量是气象因素中最容易获得的第一手资料，而潜水面蒸发量是气象因素中最容易获得的第一手资料，而潜在腾发量则需要通过试验或计算求得。根据水面蒸发量估在腾发量则需要通过试验或计算求得。根据水面蒸发量估算腾发量可采用下式：算腾发量可采用下式

29、：ET/E0a +b 式中：式中：a,b经验常数，视作物种类、生长阶段和土壤条件而定。经验常数，视作物种类、生长阶段和土壤条件而定。不同阶段采用的常数不同阶段采用的常数a,b应根据试验资料分别测定。若在各应根据试验资料分别测定。若在各生长阶段采用相同常数，则需根据作物系数加以修正。生长阶段采用相同常数，则需根据作物系数加以修正。3132l根系吸水研究有根系吸水研究有微观微观和和宏观宏观两种途径。两种途径。l微观模型微观模型采取单根径向流模型描述根系吸水过程。采取单根径向流模型描述根系吸水过程。l宏观模型宏观模型是忽略水分向单根的流动，将整个土是忽略水分向单根的流动，将整个土根系统根系统看作一个

30、连续体，在土壤水分运动中引入一个以一定模式看作一个连续体，在土壤水分运动中引入一个以一定模式分布的根系吸水速率项（源汇项）加以体现。分布的根系吸水速率项（源汇项）加以体现。在目前开展土壤植物大气连续体的模拟研究中，根系吸水在水流连续在目前开展土壤植物大气连续体的模拟研究中，根系吸水在水流连续方程中常作为一个吸水项处理，即广泛采用根系吸水的宏观模型。方程中常作为一个吸水项处理，即广泛采用根系吸水的宏观模型。l根据建立模型时经验和理论侧重的不同根据建立模型时经验和理论侧重的不同,可以将根系吸水可以将根系吸水模型分为模型分为3 3类类:经验模型经验模型、半经验半理论模型半经验半理论模型和和理论模型理

31、论模型。根系吸水模型根系吸水模型5.4 5.4 植物根系吸水规律植物根系吸水规律33经验模型（微观模型）经验模型（微观模型）经验模型是以假定植物根系分布均匀经验模型是以假定植物根系分布均匀,土壤水势均一为基础的土壤水势均一为基础的,物理学与物理学与生物学意义正确生物学意义正确,但是参数难于测定但是参数难于测定,不便于实际应用。不便于实际应用。GardnerGardner提出的单根吸水模型提出的单根吸水模型MolzMolz等人将根系吸水过程分为土壤中水分向根表面流动和水等人将根系吸水过程分为土壤中水分向根表面流动和水分在根组织内的流动两个过程分在根组织内的流动两个过程,从而考虑根的水力特性从而考

32、虑根的水力特性,提出提出的单根吸水土的单根吸水土-根系统水流运动模型。根系统水流运动模型。r为土壤中某点到根中心的径向距离为土壤中某点到根中心的径向距离;rr为根为根半径半径;为土壤含水量为土壤含水量,为土壤水势为土壤水势;0、0为初始值为初始值;()、()为土壤导水率和扩为土壤导水率和扩散率散率;为单位根长的吸水速率。为单位根长的吸水速率。m为土壤基质势为土壤基质势;c为比水容量为比水容量;rs为相邻为相邻2条根间距的中点到根轴线的距离条根间距的中点到根轴线的距离;t为作物根组织水势为作物根组织水势;Dt为根组织水分扩散系数为根组织水分扩散系数;re为根内皮层半径。为根内皮层半径。34经验模

33、型（微观模型）特点经验模型（微观模型）特点微观模型对定量描述根系区的微观区域土壤水分运微观模型对定量描述根系区的微观区域土壤水分运动规律、分析根吸水的机理、根水势和土水势的关动规律、分析根吸水的机理、根水势和土水势的关系以及蒸腾时土壤和根水势变化特征有一定的作用系以及蒸腾时土壤和根水势变化特征有一定的作用；没有考虑根层不同位置上根系密度不同所带来根系没有考虑根层不同位置上根系密度不同所带来根系吸水速率的差异吸水速率的差异,也没有考虑根系吸水分布随根系生也没有考虑根系吸水分布随根系生长的变化长的变化;将建立在单根吸水层次上模型应用于宏观整个根系将建立在单根吸水层次上模型应用于宏观整个根系系统的吸

34、水却是不可行的系统的吸水却是不可行的，所以其应用和研究受到所以其应用和研究受到一定的限制一定的限制。35Selim和和Iskandar(1978)提出蒸腾量在根系层内按各层单提出蒸腾量在根系层内按各层单位体积根长位体积根长L与导水率与导水率k的乘积分配的公式的乘积分配的公式Molz(1981)采用蒸腾量在根系层内按采用蒸腾量在根系层内按 L(hs-hr)的乘的乘积进行分配的公式积进行分配的公式半经验半理论模型半经验半理论模型Molz和和Remson(1970)提出蒸腾量在根系层内按各层单位提出蒸腾量在根系层内按各层单位体积根长体积根长L与扩散度与扩散度D的乘积分配的乘积分配36Feddes，K

35、owalik，Zaradny(1978)提出以土壤负压为基提出以土壤负压为基础的表达式础的表达式h1,h2,h3,h4分别为相应于r，d，an，s时的土壤负压值。h1h2h3h40(h)37Molz和和Remson(1976)提出作物蒸腾量提出作物蒸腾量(根系吸水总量根系吸水总量)在根层内按在根层内按4:3:2:l的比例分配。即根系吸水率在根层的比例分配。即根系吸水率在根层内按线性分布的公式内按线性分布的公式T作物的蒸腾量；作物的蒸腾量；Lr有效根层的深度。有效根层的深度。38半经验半理论模型的特点半经验半理论模型的特点半经验半理论模型能较好地考虑植物生长对水分吸收的影响半经验半理论模型能较好

36、地考虑植物生长对水分吸收的影响,以及植物根系阻力与水流推动力和水流速度的关系。其结果以及植物根系阻力与水流推动力和水流速度的关系。其结果可直接应用到田间。可直接应用到田间。但大多数半经验半理论模型没有考虑根系微区的水势梯度但大多数半经验半理论模型没有考虑根系微区的水势梯度,将将根系描述成固定的阻力网络根系描述成固定的阻力网络,不允许根系传导性随时间变化不允许根系传导性随时间变化,因此限制了该类模型的应用。因此限制了该类模型的应用。39理论模型理论模型理论模型是以电流的产生与电压差成正比的相似原理为基理论模型是以电流的产生与电压差成正比的相似原理为基础础,利用根系吸水速率与土壤和根系水势的差的关系来建模利用根系吸水速率与土壤和根系水势的差的关系来建模的。最早是的。最早是Vanden Honert(1948),具有代表性的是具有代表性的是Hillel模型模型:S(z,t)为土壤水势为土壤水势;r(t)为根水势为根水势;Rs为土壤阻力为土壤阻力,Rr为根阻力。为根阻力。理论模型为根系模型模拟研究提供了新的思路理论模型为根系模型模拟研究提供了新的思路,但是植物根但是植物根系阻力等有待进一步研究系阻力等有待进一步研究,特别是植物根系阻力与水流推动特别是植物根系阻力与水流推动力、水流速度的相关关系。力、水流速度的相关关系。