土壤水动力学SWD4土壤水运动分析ppt课件.ppt
土壤水动力学 (Soil Water Dynamics)第4章 土壤水运动分析 (Analysis of Soil Water Movement)毛 晓 敏 中国农业大学水利与土木工程学院Tel: 13621343912Email: 土壤水动力学 第4章 土壤水运动分析第4章 土壤水运动分析n入渗条件下的土壤水分运动n土壤水分再分布n蒸发条件下的土壤水分运动 入渗补给 Rg 潜水蒸发 Eg 灌溉 I 降水 P 蒸发 E 地表 地下水面 土壤水动力学 第4章 土壤水运动分析4.1 入渗条件下的土壤水分运动n入渗过程概述nGreen-Ampt入渗模型n水平入渗条件下的Philip解法n垂直入渗条件下的Philip解法n入渗条件下的Parlange解法n入渗公式小结土壤水动力学 第4章 土壤水运动分析4.1.1 入渗过程概述n入渗:水分进入土壤的过程n研究意义:水文、农田水利、水资源评价、农业及环境学n入渗类型:u一维垂直入渗:降水、灌溉入渗u二维垂直、水平入渗:河道、渠道入渗u三维垂直、水平入渗:水库、湖泊入渗n入渗条件下土壤水分运动u非饱和土壤水运动u饱和非饱和土壤水运动土壤水动力学 第4章 土壤水运动分析n积水条件下的干土入渗:u分区:t饱和区t过渡区t传导区t湿润区u饱和区、过渡区 一般不存在0 i s饱和区过渡区传导区湿润区湿润锋z土壤水动力学 第4章 土壤水运动分析n积水条件下的干土入渗:u积水后,表土含水率很快增加到0 (s )u地表处含水率梯度由大变小,t足够大时地表含水率不变u地表入渗率逐渐减小u湿润锋不断下移,含水率变化平缓0 i 0 st1t2t3z土壤水动力学 第4章 土壤水运动分析n描述土壤入渗过程的物理量:u入渗率i:单位时间内通过单位面积的入渗水量(地表水通量),mm/min, mm/h, mm/du累积入渗量I:从入渗开始到某一时刻通过单位面积的总水量,mmui与I的关系:i(t)=dI(t)/dt0)(z)(), 0()(zKDtqtiLtztztI0d)0 ,(),()(土壤水动力学 第4章 土壤水运动分析n入渗过程的影响因素:u供水速率P(降水强度)、土壤入渗能力fuPf,入渗率I 取决于入渗能力f (剖面控制), 超过f的部分产生 地表积水、径流(超渗产流)0c c积水点K0积水后a b径流积水前btpR0tpti(t)土壤水动力学 第4章 土壤水运动分析n地表入渗的边界条件:u通量边界:积水前u积水边界:积水后饱和非饱和流动 h(0, t)=h0u一类边界(灌溉模型):地表含水率接近饱和、不积水:(0, t)=0 00)()(RzKtD土壤水动力学 第4章 土壤水运动分析4.1.2 Green-Ampt入渗模型n初始干燥土壤在薄层积水时的入渗:u活塞模型:存在明显湿润锋面,将土壤分为湿土(饱和含水率)、干土(初始含水率)0 i stt1zm=-sf 湿土 =s 干土 =i H zf 土壤水动力学 第4章 土壤水运动分析nGreen-Ampt模型(水平):u地 表:水势 Hu湿润锋:水势 - sfu入渗率: i=Ks(sf+H)/zfu累积入渗量: I=(s-i) zfui=dI/dt:u积分(t=0, zf=0):u得到:ffissfzHsKtzddtDtHsKzisfsf022tDzIisfis02)()(tDdtdIiis2)(/0土壤水动力学 第4章 土壤水运动分析nGreen-Ampt模型(垂直):u地 表:水势Hu湿润锋:水势-(sf+zf)u入渗率: i=Ks(sf+zf+H)/zfu累积入渗量: I=(s-i) zfui=dI/dt:u积分(t=0, zf=0):u已知s、i、Ks、sf、H: zft, i , I 湿土 =s 干土 =i H zf fffissfzHszKtzddHsHszHszKtfffffsisln)(土壤水动力学 第4章 土壤水运动分析课后作业(3):n土壤初始含水率为0.08m3/m3,饱和含水率0.35m3/m3,饱和导水率1m/d,湿润锋处的基质吸力为0.05m。试根据Green-Ampt模型,计算绘出当垂直入渗的入渗水头为20cm时,湿润锋、入渗率和累积入渗量随时间的变化。土壤水动力学 第4章 土壤水运动分析nGreen-Ampt模型的简化情况:u地表积水很浅或入渗时间较长( zf较大)时,H可以忽略: i=Ks1+ (s-i) sf /Iut较小时: 湿润区等效水分扩散率: )/1ln(ffffsisszszKtfffissfzDzHsKtz/ddisfsHsKD)2/()(2)(2tDitDItDzisisf土壤水动力学 第4章 土壤水运动分析nGreen-Ampt模型是D()近似为函数时的解n应用:湿润锋处吸力sf难以测定n推广:u初始含水率不均匀D含水率土壤水扩散率D函数近似(x) =1 x=00 x0D()=a( - s) 土壤水动力学 第4章 土壤水运动分析4.1.3 水平入渗条件下的Philip解法n定解问题:n求解思路:u利用Boltzmann变换 将PDE转换为ODEu利用迭代方法求解 ODExtxtxDxtii0, 00)(0土壤水动力学 第4章 土壤水运动分析nBoltzmann变换:ux(, t) 方程:u分离变量法:x(, t)=() s(t)u u式中左端为t的函数,右端为的函数 constu第2式积分:s(t)=2a(t+c1)1/2 x(, t)=() 2a(t+c1)1/2xDtx)(d)(d)(dd)(1d)(d)(DttstsaDd)(d)(dd)(1attstsd)(d)(土壤水动力学 第4章 土壤水运动分析nBoltzmann变换(continue):u令()= (2a)1/2(): x(, t)=() (t+c1)1/2u入渗边界条件: (0) =0u 初始条件: (i) =x/ c11/2 u上式左端为一数,右端为一函数uc : (i) =0, i =0,稳态uc=0: x(, t)=() t1/2Boltzmann变换u土壤水运动ODE:u积分:2)()(dd)(ddDiDd)()(dd)(2土壤水动力学 第4章 土壤水运动分析n应用:u已知土壤水分扩散率D(),推求含水率分布或x Philip迭代法u根据水平土柱入渗试验结果( x ),反推土壤水分扩散率D()iDd)()(dd)(2土壤水动力学 第4章 土壤水运动分析作业(3):根据水平土柱试验结果求D()X(cm)248121620242628(cm3/cm3)0.4960.4600.4480.4420.4380.4290.4180.3950.382X(cm)2930313233343536(cm3/cm3)0.3660.3520.3440.3120.2780.0850.0680.065计算结果采用表格、绘图方式表示。计算结果采用表格、绘图方式表示。采样时间采样时间t470min 土壤水动力学 第4章 土壤水运动分析4.1.4 垂直入渗条件下的Philip解法n定解问题:n以z为变量的基本方程:LzztzKzDztii00)()(0d)(d),()(),(KtzDttz0dddd2222zDztzzKzD土壤水动力学 第4章 土壤水运动分析n级数解:u边界条件:i(0) =0,i=1, 2, 3, u初始条件: 1(i) =n将级数解代入方程,按t i/2合并,可得到关于i() 的常微分方程(ODE)u其中1() 为水平入渗解中的()u迭代计算: i() =0,i=1, 2, 3, u适用于入渗时间t较小的情况12/2/242/332/222/11)(.)()()()(),(iiittttttz土壤水动力学 第4章 土壤水运动分析n实际应用:取级数前2项,得到Philip入渗公式:uI(t)=St1/2+Atui(t)=0.5St-1/2+A u吸渗率:u稳定入渗率: 0d)(1iS)(d)(02KAi土壤水动力学 第4章 土壤水运动分析4.1.5 入渗条件下的Parlange解法n基本方程:n已知第p次迭代结果zp(, t),求导得到:n对由i到积分,假定i较小,从而D(i)0,K(i)0:n再对积分一次,积分限由z至0(地表),相应的含水率由至0,得到第p+1次迭代结果zp+1(, t): d)(d),()(),(KtzDttzttzKtzDpp),(d)(d),()(1ittzKDtzppd),()()(),(10dd),()()(),(1ittzKDtzpp土壤水动力学 第4章 土壤水运动分析u第一类边界条件u第二类边界条件土壤水动力学 第4章 土壤水运动分析4.1.6 入渗公式小结nGreen-Ampt垂直入渗公式:ui=Ks1+ (s-i) sf /I, i=ic+ b /Iu具有一定物理基础,确定Ks、 sf比较困难u对非均质土壤、初始含水率不均匀也可应用u可用于降水入渗分析(通量控制积水))/1ln(ffffsisszszKt土壤水动力学 第4章 土壤水运动分析nPhilip入渗公式:uI(t)=St1/2+At, i(t)=0.5St-1/2+Au适用范围:均质土壤垂直入渗nSmith入渗公式:t0 tp ti R ippttttAiittRi)(0土壤水动力学 第4章 土壤水运动分析n考斯加可夫(Kostiakov)入渗公式:ui(t)=Bt-au经验性公式,水平入渗nHorton入渗公式:)exp()(0tiiiicc土壤水动力学 第4章 土壤水运动分析n有关入渗的其它一些问题:u层状土的入渗:入渗率降低(相对于单一表土);手指流(finger flow)u二维/三维入渗:一般用数值方法分析t滴灌、渗灌条件下的入渗t河、渠入渗t湖、库入渗土壤水动力学 第4章 土壤水运动分析4.2 土壤水的再分布n再分布:降水、灌溉之后,水分在土壤剖面中的运动u表层水分由于蒸发、继续向下运动而减小u深层水分继续向下运动,湿润锋下移n田间持水量:n研究方法:以数值计算方法为主土壤水动力学 第4章 土壤水运动分析4.3 蒸发条件下的土壤水分运动n土面蒸发概述n地下水位恒定条件下均质土壤的稳定蒸发n蒸发条件下的土壤水分非稳定运动土壤水动力学 第4章 土壤水运动分析4.3.1 土面蒸发概述n陆面蒸发:u土面蒸发:soil evaporationu植物蒸腾(散发):plant transpirationu腾发(蒸散):evapotranspiration土壤水动力学 第4章 土壤水运动分析n蒸发强度的决定性因素:u外界条件:辐射、气温、湿度、风速等气象因素决定了大气蒸发能力t水分蒸发过程中的能量供给t蒸发面水汽向大气的扩散过程t潜在蒸发强度:土壤供水充分时,由大气蒸发能力决定的最大可能蒸发强度,一般用水面蒸发强度表示u土壤条件:土壤表层的含水率状况决定了土壤供水能力土壤水动力学 第4章 土壤水运动分析n土面蒸发的阶段性(大气蒸发能力不变)u稳定蒸发阶段:AB,k;蒸发强度E0;k取决于土质和大气蒸发能力,毛管破裂点(50%70%田持)CAB tE0 ECABk E0 E土壤水动力学 第4章 土壤水运动分析u土面蒸发强度随含水率变化阶段:BC。随着含水率的降低,向表层土壤的供水减小,地表水汽压降低,蒸发强度减小u水汽扩散阶段:表土含水率很低时,形成干土层,干土层底部蒸发水汽扩散CAB tE0 ECABk E0 E土壤水动力学 第4章 土壤水运动分析n蒸发强度的日变化对土壤水分的影响u变化模式:夜晚为0,白天按余弦规律变化u表层土壤水分也会出现周期性波动,但波动深度不大DLDaytime lengthSNSolar noonDL/ )SN(cosDL2)(tEtEdEd=6mm/d00.20.40.60.8104812162024tE土壤水动力学 第4章 土壤水运动分析n蒸发边界条件:u稳定蒸发阶段:通量边界tE=E0,k(可进一步考虑日内变化)u蒸发降低阶段:通量边界tE=E0f()=E0(a+b) , cku水汽扩散阶段:不考虑蒸发面下移,按一类边界处理:t=c土壤水动力学 第4章 土壤水运动分析蒸发机理:n道尔顿蒸发定律道尔顿蒸发定律:n道尔顿通过实验,1802年提出了反映蒸发面的蒸发速率与影响蒸发诸因素的关系式 W=C(E-e)/pan 式中的W为水面蒸发速率;(E - e)为空气的饱和差,其中E为水面温度下的饱和水汽压,e为水面上空气的实际水汽压; pa为气压;C为与风速有关的比例系数。上式表明:水面的蒸发速率,与水面上空气的饱和水汽压同实际水汽压的差值成正比,而与水面上的气压成反比,并随水面上的风速加大而增大。 土壤水动力学 第4章 土壤水运动分析蒸发机理:n空气动力学原理(紊流扩散理论)空气动力学原理(紊流扩散理论):n 近地表的空气处于紊流状态,空气质点的流径相互不平行,且以极不规则的方式相互掺混。因此,能将其中的微粒及运动属性迅速弥散,这种作用称为紊流扩散。n蒸发时水汽的运移可看做在紊流扩散作用下,垂直方向上水汽输送的结果:n比湿 用水汽压 表示, , 并令交换系数A=Kv其中Kv为水汽输送的运动粘性系数,则得到水汽输送方程:zqAEqeapq/ezeKpEva土壤水动力学 第4章 土壤水运动分析n交换系数A不随输送量属性而变,因此Kv等于紊流运动粘性系数Km 即zeKpEva21212122mvlnKKzzzzuuka212a12122aln- Ezzpeeuuk根据2个高度处的有关气象资料(风速、水汽压、气温),确定水汽、显热的垂直输送通量 土壤水动力学 第4章 土壤水运动分析4.3.2 地下水位恒定条件下土壤 的稳定蒸发n稳定蒸发条件下土壤含水率及吸力分布u含水率:u吸力:0,)(dd)(zEKzDsd)()(sEKDz0, 01dd)(zsEzssKssKEzsd)(/110土壤水动力学 第4章 土壤水运动分析u取不同K(s)、E,可得到zs分布t =E/a1,=a2+1m=1, 1.5, 2, 3, 4(粘土砂土)时的积分结果:m=1 m=2: 21)(asasKmsszsmd10)tan(/zs 1esz土壤水动力学 第4章 土壤水运动分析n恒定水位下潜水极限蒸发强度:u随着E的增加,地表含水率减小、吸力增大u吸力 时,土面稳定蒸发强度达到最大可能值,作为潜水极限蒸发强度Emaxu近似:um=1, 1.5, 2, 3, 4(粘土砂土)时, c=3.77, 2.46, 1.76, 1.52uEmax与土质、H有关msmasKasasK)/(1)(21mmmmmscHHmmaHmmaKE)/sin(/)/sin(/1maxEH0H1H2H3土壤水动力学 第4章 土壤水运动分析潜水稳定蒸发及其经验公式:u潜水蒸发是由液态水在非饱和土壤中向上运移和水汽从土面向大气散发两个过程组成。u潜水蒸发量的大小及其变化规律对于浅层地下水资源评价、土壤盐碱化成因分析和控制等十分重要。 u潜水稳定蒸发强度取决于大气蒸发能力(通常以水面蒸发强度表示)和土壤供水能力(取决于土质和潜水埋深)。土壤水动力学 第4章 土壤水运动分析 潜水蒸发经验公式:C=Eg/E20u幂函数型公式: C=aH-bu指数型公式: C=ae-bH u阿维里扬诺夫公式: C=a(1-H/Hmax)-bu清华公式:E=Emax1exp(-E20/Emax)u反Logistic公式:rHbeKC1土壤水动力学 第4章 土壤水运动分析u实例:新疆叶尔羌河均衡场砂壤土粉砂土土壤水动力学 第4章 土壤水运动分析n恒定地下水位下层状土壤的稳定蒸发u层间含水率可能不连续,基质势连续u潜水蒸发强度、各层的土壤水分运动通量和地表蒸发强度相等u可以运用以上分析,自下而上推求各层土壤基质势(含水率)分布。ssKEzsd)(/110 szsds0(s)s土壤水动力学 第4章 土壤水运动分析4.3.3 蒸发条件下的非稳定土壤水运动n解析解、半解析解(自学)n数值解(土壤大气间水热交换模型,以后再讲)