土壤水热传输理论

,单击此处编辑母版标题样式,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,土壤水动力学 第,6,章 土壤水热传输,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,土壤水动力学,(Soil Water Dynamics),第,6,章 土壤水热传输,(Simultaneous Soil Moisture and Heat Transfer),毛 晓 敏,中国农业大学水利与土木工程学院,Tel: 13621343912,Email:,第,6,章 土壤水热传输,土壤热量传输理论,蒸发条件下土壤水热传输,土壤冻融过程中的水热传输,6.1,土壤热量传输理论,土壤热特性,土壤热量传输方程,土壤水热耦合传输方程,6.1.1,土壤热特性,土壤温度变化的机理和影响因素：,热能在土壤中的传输方式：,热传导,、辐射、热对流（水的流动、,潜热,和显热输送）,土壤热传导（,Thermal conduction,）,:,刚性物体的热传导方程,-Fouriers Law:,多孔介质：,6.1.1,土壤热特性,土壤温度变化的机理和影响因素：,土壤热特性参数（,Thermal properties,）,:,比热容：热量储存能力,热导率：热量传输能力,热扩散率：热量扩散能力,土壤比热容（,Specific heat capacity,）：,单位体积（质量）的土壤温度升降,1,所吸收或释放的热量，称为体积（质量）比热容,C,v,(J cm,-3,-1,),，,C,m,(J g,-1,-1,),:,C,v,=,t,C,m,比热容是土壤的一种广度性质，具有可加性，是土壤各组分,(,矿物质、有机质、水、冰、空气,),比热容之和,C,v,=,C,vm,m,C,vom,om,C,vw,w,C,vi,i,C,va,a,常温常压：,1.9 2.5 4.18 1.78 0.001,C,m,=,C,mm,x,m,C,mom,x,om,C,mw,x,w,C,mi,x,i,C,ma,x,a,常温常压：,0.73 1.9 4.18 1.94 1.00,一般情况下，,C,v,是土壤含水率的函数：,C,v,=1.925,m,2.510,om,4.18,w,忽略有机质时：,C,v,=1.925(1-,s,),4.18,w,一定容积,(,V,),的土壤，温度从,T,1,变化至,T,2,时所需要的热量：,Q,=,C,v,V,(,T,2,-,T,1,),热导率,(thermal conductivity),根据上述推导（,Fourier,Law,）：,q,h,=,-,K,h,T,q,h,：,热通量，,W cm,-2,K,h,：热导率，,W cm,-1,-1,热导率的影响因素,土壤组分及其形状、结构、比例,对于一定的土壤，主要取决于含水率,热导率的测定：,热导率的计算方法,1,：,de Vries (1963),半经验半理论的方法,、,X,、,K,h,分别为土壤各相的加权系数、体积含量、热导率，,下标,i,、,w,、,a,分别表示第,i,种固相组成、水、空气,计算时以水为连续介质，其加权系数为,1,，其它各相的加权系数取决于该相的平均温度梯度与水的平均温度梯度之比，并且与其形状有关,假设土壤颗粒为椭球，并且其半轴长,a,、,b,、,c,满足,a,=,b,=,nc,时，其加权系数,i,与其形状系数,g,i,有关：,n,1:,n,1:,n,=1:,g,i,=1/3,土壤空气计算导热率：,包括干空气、水汽的热导率,K,ha,、,K,hv,：,K,hav,=,K,ha,+,K,hv,土壤空气形状系数,：,g,a,=1/3-(1/3-0.035),a,/,v,a,、,v,分别为单位体积土壤中的空气体积、孔隙体积,热导率的计算方法,2,：经验方法,(Chung & Horton, 1987),K,h,=,b,1,+,b,2,+,b,3,0.5,热扩散率：,D,h,=,K,h,/,C,v,Fourier,热传导定律：,q,h,=,-,K,h,T,=,-,C,v,D,h,T,D,h,：,cm,2,/s,含水率,较小时，随,的增大而增大,含水率,较大时，随,的增大而减小,6.1.2,土壤热量传输方程,热流连续方程：,各向同性土壤热量传输基本方程：,一维热传导方程：,土壤含水率不随深度变化或其变化对热特性参数影响较小时，,C,v,、,K,h,、,D,h,可视为常数,6.1.3,土壤水热耦合传输方程,土壤水分运动对热量传输的影响：,土壤水分状况影响土壤热特性参数,土壤水分相变（蒸发与凝结、冻结与融化等）产生热源、汇,地表基质势影响水汽压，从而影响地表热量平衡,土壤温度对水分运动的影响,温度影响水的物理化学性质，从而引起土壤水分运动参数、基质势的变化,温度势所产生的水分运动,土壤水热耦合传输方程,一般可忽略温度变化对水分运动参数的影响、温度势的影响，先求解水流方程，然后求解温度方程,土壤水、汽、热耦合传输方程：,(,略,),6.2,蒸发条件下土壤水热迁移模型,土壤水热迁移方程,地表与大气之间的水热交换,定水位下潜水蒸发的数值模拟,6.2.1,土壤水热迁移方程,水分：,参数：,VG,模型,温度：,参数：经验模式,6.2.2,地表与大气之间的水热交换,地表能量平衡：,R,n,=,C,s,+,E,s,+,G,R,n,：净辐射,C,s,：显热消耗,E,s,：蒸发潜热,G,：地表的热通量,微气象学水热紊流扩散理论：,显热通量： 潜热通量,C,p,：空气定压体积比热容,：湿度计常数,T,1,、,T,a,分别为地表和大气温度,e,1,、,e,a,分别为地表和大气的水汽压,r,a,、,r,s,分别为空气动力学阻力和土壤阻力,地表水汽压：,M,为水汽摩尔质量,1,为土壤表面基质势,空气动力学阻力,r,a,：,z,：参考高度,d,：零平面位移,z,0,：粗糙高度,：,Karman,常数,u,：为参考高度处风速,S,t,：大气稳定度修正,土壤蒸发阻力,r,s,(s/m),：随地表含水率的减小而增大，经验模式：,r,s,=3.5(,s,/,),2.3,+33.5,土壤水热迁移的上边界条件：,根据地表能量平衡方程、地表与大气间水热扩散方程，得到土壤水热迁移的上边界条件,水分：通量边界,E,s,温度：三类边界,应用,-,定水位下潜水蒸发的数值模拟,中壤日潜水蒸发量,E,模拟与实测值比较,(H=0.75m),粉砂土日潜水蒸发量,E,模拟与实测值比较,(H=0.75m),6.3,土壤冻融过程中的水热传输,土壤冻融过程中的水热传输方程,地下水浅埋条件下土壤冻融过程中水热迁移的模拟,考虑地气间水热交换的越冬期土壤水热迁移,6.3.1,土壤冻融过程中的水热传输方程,考虑冻融过程中相变作用的土壤水分运动方程：,u,、,i,分别为土壤未冻水、冰的体积含量,t,、,Z,分别为时间、空间坐标（垂直向下为正）,D,(,u,),、,K,(,u,),分别为非饱和土壤水分扩散率、导水率,i,、,w,分别为冰、水的密度,考虑土壤冻融过程中相变作用的热量传输方程：,T,为土壤温度,C,v,、,分别为土壤体积热容量、热导率,L,为融化潜热,冻土中未冻水含量与负温的动态平衡：,土壤冻结特性曲线：,u,m,(T),由于相变作用使得水热运动强烈耦合，带来数值求解的困难。利用水热方程消除相变项：,利用冻结特性曲线，得到：,根据以上二式可得到：,简化为水热耦合方程：,相变热容量：,C,l,=,L,w,d,m,/d,T,等效体积热容量：,C,e,=,C,v,+,C,l,等效热导率：,e,=,+,D,(,u,),C,l,等效对流速度：,U,e,=,C,l,d,K,(,u,)/d,u,6.3.2,地下水浅埋条件下 土壤冻融过程中水热迁移的模拟,试验观测：内蒙古自治区巴彦淖尔盟临河市小召乡的内蒙水科所巴音实验场,潜水埋深,2m,时的模拟结果：,土壤冻融过程：,冻融过程中土壤含水率的变化特性,特征含水量：,土壤封冻时的初始含水量,0,该层土壤冻结前的谷值含水量,v,冻结后的稳定含水量,f,消融时的峰值含水量,p,消融后的含水量,e,6.3.3,考虑地气间水热交换的 越冬期土壤水热迁移,地表能量平衡：,R,n,=,H,+,LE,+,G,R,n,：净辐射,H,：显热消耗,L,E,：蒸发潜热,G,：地表的热通量,微气象学水热紊流扩散理论：,显热通量： 潜热通量,C,p,：空气定压体积比热容,：湿度计常数,T,1,、,T,a,分别为地表和大气温度,e,1,、,e,a,分别为地表和大气的水汽压,r,a,、,r,s,分别为空气动力学阻力和土壤阻力,地表水汽压：,M,为水汽摩尔质量,1,为土壤表面基质势,空气动力学阻力,r,a,：,Z,r,：参考高度,d,：零平面位移,Z,0,：粗糙高度,：,Karman,常数,Ua,：为参考高度处风速,P,1,、,P,2,：为大气稳定度修正因子,土壤蒸发阻力,r,s,(s/m),：随地表含水率的减小而增大，经验模式：,rs=-805+4140(s-o),土壤水热迁移的上边界条件：,根据地表能量平衡方程、地表与大气间水热扩散方程，得到土壤水热迁移的上边界条件,水分：通量边界,E,s,=LE/L,w,温度：三类边界,模拟结果水分,模拟结果温度,