土壤水运动分析

,单击此处编辑母版标题样式,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,土壤水动力学 第,4,章 土壤水运动分析,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,土壤水动力学,(Soil Water Dynamics),第,4,章 土壤水运动分析,(Analysis of Soil Water Movement),毛 晓 敏,中国农业大学水利与土木工程学院,Tel: 13621343912,Email:,第,4,章 土壤水运动分析,入渗条件下的土壤水分运动,土壤水分再分布,蒸发条件下的土壤水分运动,4.1,入渗条件下的土壤水分运动,入渗过程概述,Green-Ampt,入渗模型,水平入渗条件下的,Philip,解法,垂直入渗条件下的,Philip,解法,入渗条件下的,Parlange,解法,入渗公式小结,4.1.1,入渗过程概述,入渗：,水分进入土壤的过程,研究意义：,水文、农田水利、水资源评价、农业及环境学,入渗类型：,一维垂直入渗,：降水、灌溉入渗,二维垂直、水平入渗：河道、,渠道入渗,三维垂直、水平入渗：水库、湖泊入渗,入渗条件下土壤水分运动,非饱和土壤水运动,饱和非饱和土壤水运动,积水条件下的干土入渗：,分区：,饱和区,过渡区,传导区,湿润区,饱和区、过渡区,一般不存在,积水条件下的干土入渗：,积水后，表土含水率很快增加到,0,(,s,),地表处含水率梯度由大变小，,t,足够大时地表含水率不变,地表入渗率逐渐减小,湿润锋不断下移，含水率变化平缓,描述土壤入渗过程的物理量：,入渗率,i,：单位时间内通过单位面积的入渗水量（地表水通量），,mm/min, mm/h, mm/d,累积入渗量,I,：从入渗开始到某一时刻通过单位面积的总水量，,mm,i,与,I,的关系：,i,(,t,)=d,I,(,t,)/d,t,入渗过程的影响因素：,供水速率,P,（降水强度,）、土壤入渗能力,f,P,f,，入渗率,I,取决于入渗能力,f,（剖面控制），,超过,f,的部分产生,地表积水、径流,（超渗产流）,地表入渗的边界条件：,通量边界：积水前,积水边界：积水后,饱和非饱和流动,h,(0,t,)=,h,0,一类边界（灌溉模型）：地表含水率接近饱和、不积水：,(0,t,)=,0,4.1.2 Green-Ampt,入渗模型,初始干燥土壤在薄层积水时的入渗：,活塞模型：存在明显湿润锋面，将土壤分为湿土（饱和含水率）、干土（初始含水率）,Green-Ampt,模型（水平）：,地 表：水势,H,湿润锋：水势,-,s,f,入渗率：,i,=,K,s,(,s,f,+,H,)/,z,f,累积入渗量：,I,=(,s,-,i,),z,f,i,=d,I,/d,t,：,积分,(,t,=0,z,f,=0),：,得到：,Green-Ampt,模型（垂直）,:,地 表：水势,H,湿润锋：水势,-(,s,f,+,z,f,),入渗率：,i,=,K,s,(,s,f,+,z,f,+,H,)/,z,f,累积入渗量：,I,=(,s,-,i,),z,f,i,=d,I,/d,t,：,积分,(,t,=0,z,f,=0),：,已知,s,、,i,、,K,s,、,s,f,、,H,：,z,f,t,，,i,，,I,课后作业（,3,）：,土壤初始含水率为,0.08m,3,/m,3,饱和含水率,0.35m,3,/m,3,，饱和导水率,1m/d,，湿润锋处的基质吸力为,0.05m,。试根据,Green-Ampt,模型，计算绘出当垂直入渗的入渗水头为,20cm,时，湿润锋、入渗率和累积入渗量随时间的变化。,Green-Ampt,模型的简化情况：,地表积水很浅或入渗时间较长（,z,f,较大）时，,H,可以忽略,：,i,=,K,s,1+,(,s,-,i,),s,f,/,I,t,较小时：,湿润区等效水分扩散率：,Green-Ampt,模型是,D,(,),近似为,函数时的解,应用：湿润锋处吸力,s,f,难以测定,推广：,初始含水率不均匀,(x) =,1 x=0,0 x0,D(,)=a,(, - ,s,),4.1.3,水平入渗条件下的,Philip,解法,定解问题：,求解思路：,利用,Boltzmann,变换,将,PDE,转换为,ODE,利用迭代方法求解,ODE,Boltzmann,变换：,x,(,t,),方程：,分离变量法：,x,(,t,)=,(,),s,(,t,),式中左端为,t,的函数,右端为,的函数,const,第,2,式积分：,s,(,t,)=2,a,(,t,+,c,1,),1/2,x,(,t,)=,(,) 2,a,(,t,+,c,1,),1/2,Boltzmann,变换,(continue),：,令,(,),=,(2,a,),1/2,(,),：,x,(,t,)=,(,) (,t,+,c,1,),1/2,入渗边界条件：,(,0,) =0,初始条件：,(,i,) =,x,/,c,1,1/2,上式左端为一数，右端为一函数,c,：,(,i,) =0,，,i,=,0,，稳态,c,=0,：,x,(,t,)=,(,),t,1/2,Boltzmann,变换,土壤水运动,ODE,：,积分：,应用：,已知土壤水分扩散率,D,(,),，推求含水率分布,或,x,Philip,迭代法,根据水平土柱入渗试验结果（,x,），反推土壤水分扩散率,D,(,),作业（,3,）：根据水平土柱试验结果求,D(,),X(cm),2,4,8,12,16,20,24,26,28,(cm3/cm3),0.496,0.460,0.448,0.442,0.438,0.429,0.418,0.395,0.382,X(cm),29,30,31,32,33,34,35,36,(cm3/cm3),0.366,0.352,0.344,0.312,0.278,0.085,0.068,0.065,计算结果采用表格、绘图方式表示。,采样时间,t,470min,4.1.4,垂直入渗条件下的,Philip,解法,定解问题：,以,z,为变量的基本方程：,级数解：,边界条件：,i,(,0,) =0,，,i,=1, 2, 3, ,初始条件：,1,(,i,) =,将级数解代入方程，按,t,i,/2,合并，可得到关于,i,(,),的常微分方程（,ODE,）,其中,1,(,),为水平入渗解中的,(,),迭代计算：,i,(,) =0,，,i,=1, 2, 3, ,适用于入渗时间,t,较小的情况,实际应用：取级数前,2,项，得到,Philip,入渗公式：,I,(,t,)=,St,1/2,+,At,i,(,t,)=0.5,St,-1/2,+,A,吸渗率：,稳定入渗率：,4.1.5,入渗条件下的,Parlange,解法,基本方程：,已知第,p,次迭代结果,z,p,(,t,),，求导得到：,对,由,i,到,积分，假定,i,较小，从而,D,(,i,)0,，,K,(,i,)0,：,再对,积分一次，积分限由,z,至,0,（地表），相应的含水率由,至,0,，得到第,p,+1,次迭代结果,z,p,+1,(,t,),：,第一类边界条件,第二类边界条件,4.1.6,入渗公式小结,Green-Ampt,垂直入渗公式：,i,=,K,s,1+ (,s,-,i,),s,f,/,I,i,=,i,c,+,b,/,I,具有一定物理基础，确定,K,s,、,s,f,比较困难,对非均质土壤、初始含水率不均匀也可应用,可用于降水入渗分析（通量控制积水）,Philip,入渗公式：,I,(,t,)=,St,1/2,+,At,i,(,t,)=0.5,St,-1/2,+,A,适用范围：均质土壤垂直入渗,Smith,入渗公式：,考斯加可夫,(Kostiakov),入渗公式：,i,(,t,)=,Bt,-,a,经验性公式，水平入渗,Horton,入渗公式：,有关入渗的其它一些问题：,层状土的入渗：入渗率降低（相对于单一表土）；手指流（,finger flow,）,二维,/,三维入渗：一般用数值方法分析,滴灌、渗灌条件下的入渗,河、渠入渗,湖、库入渗,4.2,土壤水的再分布,再分布：降水、灌溉之后，水分在土壤剖面中的运动,表层水分由于蒸发、继续向下运动而减小,深层水分继续向下运动，湿润锋下移,田间持水量：,研究方法：以数值计算方法为主,4.3,蒸发条件下的土壤水分运动,土面蒸发概述,地下水位恒定条件下均质土壤的稳定蒸发,蒸发条件下的土壤水分非稳定运动,4.3.1,土面蒸发概述,陆面蒸发：,土面蒸发：,soil evaporation,植物蒸腾（散发）：,plant transpiration,腾发（蒸散）：,evapotranspiration,蒸发强度的决定性因素：,外界条件：辐射、气温、湿度、风速等气象因素决定了大气蒸发能力,水分蒸发过程中的能量供给,蒸发面水汽向大气的扩散过程,潜在蒸发强度：土壤供水充分时，由大气蒸发能力决定的最大可能蒸发强度，一般用水面蒸发强度表示,土壤条件：土壤表层的含水率状况决定了土壤供水能力,土面蒸发的阶段性,(,大气蒸发能力不变,),稳定蒸发阶段：,AB,，,k,；蒸发强度,E,0,；,k,取决于土质和大气蒸发能力，毛管破裂点（,50%,70%,田持）,土面蒸发强度随含水率变化阶段：,BC,。随着含水率的降低，向表层土壤的供水减小，地表水汽压降低，蒸发强度减小,水汽扩散阶段：表土含水率很低时，形成干土层，干土层底部蒸发水汽扩散,蒸发强度的日变化对土壤水分的影响,变化模式：夜晚为,0,，白天按余弦规律变化,表层土壤水分也会出现周期性波动，但波动深度不大,DL,Daytime length,SN,Solar noon,蒸发边界条件：,稳定蒸发阶段：通量边界,E=E,0,，,k,（可进一步考虑日内变化）,蒸发降低阶段：通量边界,E=E,0,f,(,),=E,0,(,a,+,b,),，,c,k,水汽扩散阶段：不考虑蒸发面下移，按一类边界处理：,=,c,蒸发机理：,道尔顿蒸发定律,：,道尔顿通过实验，,1802,年提出了反映蒸发面的蒸发速率与影响蒸发诸因素的关系式,W=C,（,E-e,）,/p,a,式中的,W,为水面蒸发速率；（,E - e,）为空气的饱和差，其中,E,为水面温度下的饱和水汽压，,e,为水面上空气的实际水汽压；,p,a,为气压；,C,为与风速有关的比例系数。上式表明：水面的蒸发速率，与水面上空气的饱和水汽压同实际水汽压的差值成正比，而与水面上的气压成反比，并随水面上的风速加大而增大。,蒸发机理：,空气动力学原理（紊流扩散理论）,：,近地表的空气处于紊流状态，空气质点的流径相互不平行，且以极不规则的方式相互掺混。因此，能将其中的微粒及运动属性迅速弥散，这种作用称为紊流扩散。,蒸发时水汽的运移可看做在紊流扩散作用下，垂直方向上水汽输送的结果：,比湿 用水汽压,表示，,并令交换系数,A=,K,v,其中,K,v,为水汽输送的运动粘性系数，则得到水汽输送方程：,交换系数,A,不随输送量属性而变，因此,Kv,等于紊流运动粘性系数,K,m,即,根据,2,个高度处的有关气象资料（风速、水汽压,、气温,），确定水汽,、显热,的垂直输送通量,4.3.2,地下水位恒定条件下土壤 的稳定蒸发,稳定蒸发条件下土壤含水率及吸力分布,含水率：,吸力：,取不同,K,(,s,),、,E,，可得到,z,s,分布,=,E,/,a,1,，,=,a,2,+1,m,=1, 1.5, 2, 3, 4(,粘土砂土,),时的积分结果,:,m=1,m,=2:,恒定水位下潜水极限蒸发强度：,随着,E,的增加，地表含水率减小、吸力增大,吸力 时，土面稳定蒸发强度达到最大可能值，作为潜水极限蒸发强度,E,max,近似：,m,=1, 1.5, 2, 3, 4(,粘土砂土,),时，,c,=3.77, 2.46, 1.76, 1.52,E,max,与土质、,H,有关,E,H,0,H,1,H,2,H,3,潜水稳定蒸发及其经验公式：,潜水蒸发是由液态水在非饱和土壤中向上运移和水汽从土面向大气散发两个过程组成。,潜水蒸发量的大小及其变化规律对于浅层地下水资源评价、土壤盐碱化成因分析和控制等十分重要。,潜水稳定蒸发强度取决于大气蒸发能力（通常以水面蒸发强度表示）和土壤供水能力（取决于土质和潜水埋深）。,潜水蒸发经验公式：,C,=,E,g,/,E,20,幂函数型公式：,C,=,aH,-b,指数型公式：,C,=,a,e,-,bH,阿维里扬诺夫公式：,C,=,a,(1-,H,/,H,max,),-b,清华公式：,E,=,E,max,1,exp(-,E,20,/,E,max,),反,Logistic,公式：,实例：新疆叶尔羌河均衡场,砂壤土,粉砂土,恒定地下水位下层状土壤的稳定蒸发,层间含水率可能不连续，基质势连续,潜水蒸发强度、各层的土壤水分运动通量和地表蒸发强度相等,可以运用以上分析，自下而上推求各层土壤基质势（含水率）分布。,(s),s,4.3.3,蒸发条件下的非稳定土壤水运动,解析解、半解析解（自学）,数值解（土壤大气间水热交换模型，以后再讲）,