第6章-下渗

降水蒸发径流水汽输送入渗蒸发降水水汽输送水文循环：水圈中的各种水体在太阳辐射和地心引力作用下通过不断蒸发、水汽输送、凝结、降落、下渗、地面和地下径流的往复循环过程。（降雨、蒸发、径流、下渗）内因是根据，外因是条件，内因通过外因起作用水分由海洋输送到大陆，又回到海洋的循环称为大循环。水分由海洋输送到海洋，又回到海洋，叫海洋小循环。水分由陆地输送到陆地，又回到陆地，叫内陆小循环。水量平衡原理要正确理解水量平衡方程必须掌握两个关键点:（1）计算时段（针对一定长度的时间段）（2）研究区域(针对某一空间尺度)水量平衡是水文循环遵循物质不灭定律的具体表现，水量平衡方程式实际就是水量收支平衡关系式，具体表达式为：流域水量平衡非闭合流域：In:P+qiOut:q+E+Rs+Rg+qo(P+qi)－(q+E+Rs+Rg+qo)=△W闭合流域：qi=qo=0若q=0，则In:POut:E+Rs+RgP－(E+Rs+Rg)=△W令R=Rs+Rg,则P-(E+R)=△W多年平均：P=R+EqoqiRgRsqPE流域与周围区域的地下水交换降雨量（深）：指一定时段内降落在某一点或某一面积上的总雨量，用深度表示，以mm计。降雨历时：降雨从某时刻到另一时刻所经历的时间称为降雨历时；一次降雨从开始到结束所经历的时间称为次降雨历时，以min,h或d计。降雨强度：单位时间内的降雨量称为降雨强度，以mm/min或mm/h计。降雨面积：降雨笼罩范围的水平投影面积称为降雨面积，km2计。暴雨中心：暴雨集中的较小的局部地区，称为暴雨中心。降水降水过程线降水累计过程线等雨量线降水特征综合曲线降雨资料的代表性、一致性和可靠性利用双累积曲线检验降水资料的一致性。双累计曲线是指被检验雨量站的累积降雨量与其周围若干雨量站平均值的累计雨量的相关曲线。流域平均降水量的计算算数平均法泰森(Thiessen)多边形法等雨量线法距离平方倒数法土壤水自由水｛束缚水｛吸湿水膜状水毛管水重力水｛渗透重力水支持重力水｛毛管悬着水毛管上升水｛土壤水类型最大吸湿水量：在饱和空气中，干燥土粒能够吸附的最大水量。又称吸湿系数最大分子持水量：土粒分子力所结合水分的最大量，薄膜水厚度达最大值。凋萎含水量（凋萎系数）：植物根系的吸力约为15个大气压，对于土粒吸附的吸力大于该值的水分，植物则无法利用。当土壤水分低于这时的含水量时，植物将缺水而凋萎死亡，该土壤含水量称为凋萎含水量。土壤水分常数毛管断裂含水量：湿润的土壤逐渐干燥时，毛管悬着水的连续状态开始断裂，此时的土壤含水量称毛管断裂含水量。土壤含水量低于该值后，土壤中的水分只能以水汽和薄膜水的形式向蒸发面运移，约为田持的65%。田间持水量：土壤能够保持而不在重力作用下流走的最大含水量，称田间持水量。这时继续下渗的雨水，将补给潜水，形成地下径流。饱和含水量：土壤中的孔隙全部被水充满情况下的土壤含水量。土壤水分常数(一)土水势土壤水在各种力作用下，与同样温度、高度和大气压条件的纯自由水相比，其自由能降低，这个差值即为土水势（）。优点：可作为判断各种土壤水分能态的统一标准和尺度重力势ψg取决于水分在重力场中的位置压力势ψp包括气压势和静水压势基质势ψm是指由分子力和毛管力引起的势能的总称。基质势总是为负值溶质势ψs又称渗透势，负值。土壤溶质浓度越高，溶质势越低。温度势ψrBA10CM2CM4CM如图：已知某饱和土柱，各数据见图，又知传导度K=310-8，求A、B间的运动方向及速度大小。解：对Aψg=10cm,ψp=2cm,Φ=12cm对Bψg=0cm,ψp=6cm,Φ=6cm故水流方向A—BV=-310-8（6-12）/10=1.810-8m/s第六章下渗下渗的物理过程下渗理论及经验下渗曲线天然条件下的下渗土壤水的再分布第一节下渗的物理过程一、与下渗有关的基本概念土壤含水率沿深度方向的变化曲线称为土壤水分剖面（土壤含水率垂向分布），它描述了土壤含水量在深度方向上的分布情况。水分透过土壤层沿垂直和水平方向渗入土壤中的现象称为下渗。下渗率f：又称下渗强度。指单位面积上、单位时间内渗入土壤中的水量。下渗能力fp：又称下渗容量。指在充分供水条件下的下渗率。累积下渗量F：入渗开始后一段时间内，通过单位面积下渗到土壤中总水量。dF/dt=f一、与下渗有关的基本概念下渗曲线确切地说，应称下渗能力曲线，指地面充分供水条件下下渗率随时间的变化过程线。累计下渗曲线：累计下渗量随时间的变化曲线。实际下渗率：累计下渗曲线上任一点切线的斜率为某一时刻的下渗率稳定下渗率：随时间的推移下渗率逐渐减小，最后趋于较稳定的值并持续下降，称此时的下渗率为稳定下渗率。0.020.040.060.080.0100.0050100150时间(min)累积下渗量(mm)累积下渗量曲线0.001.002.003.004.005.00050100150时间(min)下渗率(mm/min)下渗率曲线二、下渗的物理过程根据水分所受作用力及运动特征，干燥土壤在充分供水条件下的下渗分三阶段：渗润阶段：主要受分子力作用，入渗水成为薄膜水，当土壤含水量达到最大分子持水量时结束。渗漏阶段：主要受毛管力、重力作用，入渗水主要成为毛管水，当土壤含水量达到饱和含水量时结束。渗透阶段：受重力作用，入渗水成为自由重力水向下渗出。下渗过程中土壤水分动态及分布规律讨论条件：积水下渗、无地下水、初始含水量均匀、供水时间长含水量(%)深度(m)饱和含水量田间持水量风干土饱和带过渡带水分传递带湿润带湿润锋饱和带过渡带水分传递带湿润带湿润锋三、下渗容量与土壤水分剖面的关系0(,)npsFztdKt0(,)nppsdFdfztdKdtdt在供水条件充分大的条件下已知某地土壤累积下渗能力曲线及某次降雨的累积雨量曲线如图，累积下渗量曲线F(t)=36t0.5+0.4t,累积降雨量曲线为P(t)=9.4t(1)求该地的下渗能力曲线f(t)及该次降雨强度过程(2)由该次降雨求出产流开始的时间F(mm)P(t)F(t)t(分钟)解(1)f(t)=dF(t)/dt=18t-0.5+0.4I(t)=dP(t)/dt=9.4(mm/分钟)(2)f(t)=I(t)18t-0.5+0.4=9.4t=4分第二节下渗理论与公式1非饱和下渗理论第二节下渗理论与公式1非饱和下渗理论条件1.忽略重力；2.供水充分、表面无积水；3.均质半无限土柱，初始土壤含水量分布均匀。定解问题（1）、忽略重力作用的下渗isittzzDzt),(),0()0,(])([边界条件初始条件方程第二节下渗理论与公式1非饱和下渗理论问题求解①当D()=D为常数时，问题变为：isittzzDt),(),0()0,(22令y(z,t)=(z,t)-i,则：0),(),0(0)0,(022tyytyzyzyDtyis)],([~~)(~)0,()],([)()0()0()0()]([)]([22221121tzyLydzydzyLyPzytzyPLtyLffPfPtfLPtfLnnnnn以z为参数，将y(z,t)关于t作拉氏变换:（1）、忽略重力作用的下渗第二节下渗理论与公式1非饱和下渗理论（1）、忽略重力作用的下渗问题求解（续）求逆变换得：)2(0Dtzerfcyy拉氏变换中，象函数PaeP1对应的原函数为)2(taerfc0),(~),0(~0)0,(~0~~022tyPytyzyyPdzydD经拉氏变换后问题为:zDPePyy0~解为:Dtz2deerfcyy220)(第二节下渗理论与公式1非饱和下渗理论（1）、忽略重力作用的下渗问题求解（续）得到原问题的解为：detzisi2)(2),(2104000)()1)(()()(|2tDDteDzDzDVfiszDtziszzzp下渗能力曲线形状为：fpt问题求解1.当D()=D为常数时，问题变为：isittzzDt),(),0()0,(22下渗能力曲线形状为：fpt21)(tDfisp(求解过程不展开)第二节下渗理论与公式1非饱和下渗理论（1）、忽略重力作用的下渗第二节下渗理论与公式1非饱和下渗理论（1）、忽略重力作用的下渗②当D为的函数时，求解过程不展开，结论：212121)()()(StdttdFtfSttFp第二节下渗理论与公式1非饱和下渗理论（2）、考虑重力作用的下渗定解问题：Z向下为正isittzzkzDzt),(),0()0,()(])([边界条件初始条件方程结论：BStdttdFtfBtSttFp212121)()()(fpt第二节下渗理论与公式2饱和下渗理论一、基本假定1.半无限土柱，初始土壤含水量分布均匀。2.地面积水深hp；3.下渗锋面以上是饱和的，=S，K=KS；4.下渗锋面以下为初始土壤含水量，吸力hS。SBPAhzhSiABZ下渗锋面位置第二节下渗理论与公式2饱和下渗理论二、公式推导fp～Z的关系充分供水条件下单位时间单位面积上入渗水量=fpZ——下渗锋面位移。dtdzfzFzhKhhzhzhKqfzhzhKdzKdzdKqispissssppssppssABss)()()1()(z(t)的解21)(2)21()(])(21)([)(3121)1ln()]1ln([)()1ln(11)111()(11)(]1[)(2232thKzthzKthzhhzhzKuuuuthzhzKuduudtdzhzKdtdzzhKzhKdtdzisssssissssssissssisssisssisssis取前两项：因两边积分：：而即：第二节下渗理论与公式2饱和下渗理论二、公式推导第二节下渗理论与公式2饱和下渗理论二、公式推导fp(t)的解212)(thKtzisss21)(5.0)(thKKtfsisssp饱和下渗理论和非饱和下渗理论推得的下渗曲线均为t-1/2的函数，为下渗经验公式的提出奠定了理论基础。]1[)(zhKtfssp第二节下渗理论与公式2饱和下渗理论例题内的入渗总量。，和入渗强度时的求入渗锋面位置公式：已知水平方向入渗的min30min30)2(2121iStthKSAmptGreenfitftfmmSIsmShKdxdKiShmmthKSitfffttffitftf82.2)(3/1083.70)2(05.7)2(172121入渗总量）（入渗强度：距离为，入渗点至入渗锋面的，入渗锋面入渗点）入渗锋面位置（解：Kostiakov公式Horton公式Green-Ampt公式Philip公式Holtan公式Smith公式Smith-Parlange公式第二节下渗理论与公式3经验下渗曲线batft0)(ti0t)(titKostiakov公式该公式是Kostiakovt提出的，表达式为：时，只有在水平吸渗的情况下才出现，垂直入渗条件下，显然不符合。但在实际情况下，只要能确定出t的期限，使用该公式还是比较简便而且较为准确。2/12taf霍顿公式tkcceffff)(0Green-Ampt公式Green-Ampt公式是研究初始干燥土壤在薄层积水条件下入渗问题。基本假定是，入渗时存在明确的水平湿润锋面同时具有固定