土壤水动力学复习笔记

学习好资料学习好资料欢迎下载土壤水动力学是许多学科的基础，它的研究涉及农田水利学、水文学、地下水文学、水文地质学、土壤物理学、环境科学等学科。1）合理开发和科学管理水资源;2）调控农田墒情，促进农业节水；3）土壤改良和水土环境的改善。土壤各个指标，计算意义，相互关系。土壤一是由矿物质和生物紧密结合的固相、液相和气相三相共存的一个复杂的、多相的、非均匀多孔介质体系。定性指标一质地、结构。定量指标-孔隙度、密度、含水率、饱和度等。⑶含水率。体积含水率：0v二Vw/VO重量（质量）含水率：0g二mw/ms饱和度：w=Vw/Vv贮水深度：h=H0（量刚为L）主要测定方法：称重法（烘干法）、核技术测量：中子仪，Y射线仪、电磁测量：时域反射仪TDR）、核磁共振测量、热脉冲测量、遥感测量：大面积地表含水率；⑷水分常数。吸湿水，束缚在土粒表面的水汽，最大吸湿量（吸湿常数）薄膜水，吸湿水外层连续水膜，最大分子持水量（薄膜水不能被植物吸收时）凋萎系数；毛管水，土壤孔隙（毛管）水气界面为一弯月面，分毛管上升水、毛管悬着水，田间持水量（毛管悬着水达到最大）田持；重力水，大孔隙中的水，饱和含水率。农业生产中常用的水分常数：田间持水量field（moisture）capacity）：农田土壤某一深度内保持吸湿水、膜状水和毛管悬着水的最大水量。凋萎系数Wiltingcoefficient、：土壤中的水分不能被根系吸收、植物开始发生永久凋萎时的土壤含水率，也称凋萎含水率或萎蔫点。土壤有效含水量（availablewatercontentofsoil）:土壤中能被作物吸收利用的水量，即田间持水量与凋萎系数之间的土壤含水量。土壤含水率与水分常数的应用：估计水分对植物生长的影响；计算灌溉水量；根据土壤水分的动态变化估算腾发量（地面蒸发+植物蒸腾）⑸土水势（Soilwaterpotential）：可逆、等温地从特定高度和大气压下的纯水池转移极少量水到土壤中某一点时单位数量纯水所做的功。定义土壤中任一点的单位数量土壤水分的吉氏自由能与标准参考状态下自由能的差值为该点的总土水势妙二〃p+/T+〃s+〃m+/g土水势二压力势+温度势+溶质势+基质势+重力势，标准状态下p=0，将单位数量的水分从标准状态移动到另一状态时，如果环境对土壤水做功屮>0;如果土壤水对环境做功，屮<0。重力势：地球重力场对土壤水作用的结果；大小：取决于土壤水相对于参考面的高度；势能：Eg=±Mgz。压力势：压力场中压力差的存在而引起的；大小：取决于水压与大气压之差；势能：Ep=VAp;饱和地下水：地下水面以下深度h，屮p$0;非饱和土壤水：气孔连通：p=0，封闭未充水孔隙：气压势（目前考虑较少、基质势：由土壤基质对土壤水的吸持作用（毛管作用、吸附作用）所产生，自由水被土壤吸持后,自由能降低，土水势减小（0—负值），屮m<0，饱和一非饱和水分运动研究中：负压势h土壤基质对土壤水分吸持作用：与土壤含水量有关，屮m〜0（土壤水分特征曲线）；基质势的测定：张力计（负压计）法，砂性漏斗法，压力仪法，离心机法，稳定土壤水分剖面法。溶质势：土壤溶液中溶质对土壤水分综合作用的结果，渗透试验（半透膜试验）：溶质的存在降低了水的势能ips<0单位体积土壤水的溶质势：屮s二-（c/M）RTc:溶液浓度（g/cm3）;M:溶质摩尔质量（g/mol），c/M:溶液摩尔浓度（mol/cm3），R：摩尔气体常数，8.314MPacm3/molK，T：热力学温度（K）。土壤：不存在半透膜，一般不考虑溶质势；需要考虑溶质势的情况：植物根系吸水，植物细胞渗流，水汽扩散温度势：温差；大小：屮T=-SeAT，Se：单位数量土壤水分的熵，不易定量描述；温差对土壤水分运动影响不大，因此一般不考虑温度势；温度对土壤水分运动的影响：温度影响水的物理化学性质（粘滞性、表面张力、渗透压等从而影响基质势、溶质势和土壤水分运动参数，温度决定水的相变和热特性参数。其他分势：荷载势：土壤承受荷载，湿润势：膨胀土在饱和状态下产生的土水势。小结，一般情况下：土壤水：二屮m土z，地下水：屮=h±z，存在半透膜时：考虑溶质势。

⑹基质势（土壤水吸力）随土壤含水率而变化，其关系称为土壤水分特征曲SOilWatercharacteristiccurve）或土壤持水曲线（Soilwaterretentioncurve）土壤水分特征曲线反映了土壤水和土壤基质间的相互作用：高吸力：以吸附作用为主，低吸力：以毛管作用为主。基本形状，土壤饱和：s二-屮m=0，对土壤施加一定吸力，吸力较小时无水分排出，含水率维持饱和值，当吸力增加至某一临界值a后，土壤水分开始排出，含水率减小。sa称为进气值；一般地，粗质地的砂性土壤或结构良好的土壤，其a较小；相反，细质地的粘性土壤，其Sa较大；吸力不断增大，含水率不断减小。土壤质地：相同吸力下，不同土质的含水率差别较大；土壤的粘粒含量越高，同一吸力条件下土壤的含水率越大；相同含水率对植物的有效性不同。土壤结构：土壤压实后，孔隙度减小，大孔隙减少，饱和含水率减小，中等孔隙增加，小孔隙变化不大；一般的，土壤的粘粒含量越高，同一吸力条件下土壤的含水率越大。温度：温度升高：水的粘滞性、表面张力下降,基质势增大（基质吸力减小;）土壤水分变化历史：脱湿，吸湿。土壤水分特征曲线并非单值曲线，主脱湿曲线（土壤从饱和到干燥过程曲线称为主脱湿线湿曲线（土壤从干燥到饱和过程曲线称为主吸湿线）、扫描曲线（土壤从部分湿润开始排水或廓分干燥到吸湿过程线）。相同吸力下脱湿过程的含水率较大；砂土滞后现象明显。滞后现象在粗质的土壤于低吸力范围内最为明显；孔隙排水时的吸力较它们充水时吸力大得多。土壤水分特征曲线的应用，进行土壤吸力和含水率e之间的换算。间接地反映出土壤中孔隙大小的分布。可用来分析不同质地土壤的持水性和土壤水分的有效性。应用数理方法定量分析土壤水分运动时，水分特征曲线和比水容量都是必不可少的重要参数。[7]Darcy定律的微分形心=_kdH三维：q=-KgradH=-K▽HsssdL梯度：...QH.QH.aH.qradH^14kaxayazDarcy定律的适用范围：小Re数层流：与粘滞力相比，惯性力作用可以忽略不计。在紊流状态下，通量与水势梯度呈非线性关系；对于颗粒极细的土壤：克服一定的初始水头差才能发生流动；一般情况下，Darcy定律有效。饱和导水率Ks:综合反映了多孔介质对流体流动的阻碍作用，多孔介质的基质特征：质地、结构•…流体物理性质：粘滞性、密度•…实验室测定：现场测定：双环入渗试验，Guelph渗透仪，抽水试验。1907年：EdgarBuchkingham将Darcy定律推广到非饱和土壤水：饱和:q二-K▽H*非饱和：q二-K（切▽屮,q二-K（0）▽屮q二-K（切▽（^土z），驱动力：土水势（重力势+基质势）梯度不能笼统地说水由高处流向低处，或湿处流向干处。导水率：小于饱和导水率，是基质势（含水率）的函数。非饱和导水率随基质势（含水率）的减小而减小的原因：部分孔隙充气，随着含水率的降低，实际过水面积减小；随着含水率的降低，较大孔隙排水，土壤水在较小的孔隙流动，水流阻力增大，实际流速减小；小孔隙弯曲程度增加非饱和导水率的影响因素：与土壤质地有关，（Ex:Miller&Gardner,1962湿润情况下：砂性土K＞粘性土K，干燥情况下：砂性土＜粘性土K;与土壤结构有关，（Ex:渠底夯实；农田表面结皮），土壤干容重增大，K减小；K与含水率关系受滞后作用影响较小，但与基质势（or吸力）关系则受滞后影响。⑻Richards方程的不同形式，混合形式:aK(⑻Richards方程的不同形式，混合形式:aK(0)ayyrK(0)aK(0)az不同形式基本方程的特点：混合方程是一般形式；方程：数学处理，适用于均质非饱和土壤，扩散率D的变化比K小；屮m方程：可用于饱和一非饱和流动、土壤分层等情况；K的变化范围大，数值计算时需要特别处理以保证质量守恒；以或z为因变量的基本方程：简单情况下的解析解和半解析解。初始条件：所研究问题初始状态，即初始时刻自变量在研究区域上的分布对于8型方程，需已知对于屮型方程，需已知(比沪Q-0恢川“或写作：h(x,y,z,O)=hO(x,y,z)m边界条件，一般分为三类：第一类边界条件(变量已知边界Dirichlet条件)对于8型方程，0(x,y,z,t)=9(x,y,z,t)(x,y,z)ier对于屮型方程，h(x,y,z,t)=h(x,y,z,t)(x,y,z)eV厂为一类边界区域举例：地面薄层积水入渗时，地表可视为一类边界条件；土壤下边界若选在潜水面处，潜水位不变时常视为一类边界条件第二类边界条件(水流通量已知边界Neuman条件)q=-K⑼▽讪±z)-D(8)旳+K(8)=q(t)m~zr2♦垂向一维：-K(h)ah+K(h)=q(t)r2举例：通量为零的情况：如不透水边界、无蒸发入渗的边界通量已知的情况：降雨、灌溉、蒸发强度已知第三类边界条件(水流通量随边界上的变量变化而变化的情况)dfca+af=a1^z23♦举例：当土壤蒸发强度为表土含水率(基质势)函数的情况-D(-D(8)旳+K(8)=a8+bah-K(h)二+K(h)二af(h)+baz3零通量面法：零通量面-当水势梯度为0,该处通量为0,则该处为零通量面分类：单一聚合型零通量面；单一发散型零通量面；多个零通量面。由于零通量面为已知通量(零)断面，若1和t2时段内零通量位置不变，则根据两时刻的土壤含水率观测值，可计算出时段内任一断面处流过的土壤水通量。⑼入渗：水分进入土壤的过程。积水条件下的干土入渗：分区饱和区，过渡区，传导区，湿润区；饱和区、过渡区一般不存在；积水条件下的干土入渗：积水后，表土含水率很快增加到e0(<es)；地表处含水率梯度由大变小,，足够大时地表含水率不变；地表入渗率逐渐减小；湿润锋不断下移，含水率变化平缓。入渗率单位时间内通过单位面积的入渗水量(地表水通量mm/min,mm/h,mm/d，累积入渗量I：从入渗开始到某一时刻通过单位面积的总水量，mm。入渗过程的影响因素：供水速率P(降水强度…)、土壤入渗能力f，P<f，入渗率i取决于供水速率P(通量控制)；P>f，入渗率I取决于入渗能力f(剖面控制)，超过f的部分产生地表积水、径流(超渗产流)地表

入渗的边界条件：通量边界：积水前-D(0)「爲+K(e=RL贡J|z=00积水边界：积水后——饱和一非饱和流动h(0,t)=h一类边界(灌溉模型：地表含水率接近饱和、不积水：(0,t)=^GreenAmpt入渗模型初始干燥土壤在薄层积水时的入渗：'活塞模型：存在明显湿润锋面，将土壤分为湿土(饱和含水率)干土(初始含水率).水平，地表：水势H,湿润锋：水势一sf，入渗率:i二Ks(sf+H)/zf，累积入渗量：I=(0s-0i)zf，i=dl/dt：dzks+Hz2s+H-cf=—积分(t=o，zf=o)：f=Ks0f-0t=D/i=dl/dt=i=dl/dt=(0-0)耳得到：I=(0-0)z=(0-0).2DTsifsi'0垂直•地表水势H，垂直•地表水势H，湿润锋：水势+z)，入渗率i=K(s+z+H)/z，累积入渗量:Ksit=0-0si=dl/dt：积分(t=0,z=0)：z+s+H_ff，z.fz―(s+H)lnLff=呼尸入渗公式小结Green-Ampt垂直入渗公式：i=K[1+僧｛)s/I],t0-0[S…]i=i+b/Isi7—sln(1+z/s)czffffs具有一定物理基础，确定C)s比较困难；对非均质土壤、初始含水率不均匀也可应用；可用于降水入渗分析(通量控制一积水)■Philip入渗公式：I(t)=St1/2+At,i(t)=0.5St-1/2+A适用范围：均质土壤垂直入渗fi=Rt<tSmith入渗公式｛ij+A(t-1)-at>tp»0p■考斯加可夫(Kostiakov)入渗公式:i(t)二Bt-a，经验性公式，水平入渗Horton入渗公式'='+(i-i)exP-pt)c0c[10]蒸发强度的决定性因素：外界条件：辐射、气温、湿度、风速等气象因素决定了大气蒸发能力；水分蒸发过程中的能量供给；蒸发面水汽向大气的扩散过程；潜在蒸发强度：土壤供水充分时，由大气蒸发能力决定的最大可能蒸发强度，一般用水面蒸发强度表示土壤条件：土壤表层的含水率状况决定了土壤供水能力。土面蒸发的阶段性(大气蒸发能力不变)稳定蒸发阶段：ab，e>ek；蒸发强度e；ek取决于土质和大气蒸发能力，毛管破裂0点(50%〜70%田持)土面蒸发强度随含水率变化阶段：BC。随着含水率的降低，表层土壤的供水减小，地表水汽压降低，蒸发强度减小

水汽扩散阶段：表土含水率很低时，形成干土层，干土层底部蒸发一水汽扩散潜水蒸发是由液态水在非饱和土壤中向上运移和水汽从土面向大气散发两个过程组成。潜水蒸发量的大小及其变化规律对于浅层地下水资源评价、土壤盐碱化成因分析和控制等十分重要。潜水稳定蒸发强度取决于大气蒸发能力(通常以水面蒸发强度表示)和土壤供水能力(取决于土质和潜水埋深[11]对流(Convection)：溶质随土壤水分的运动而运移，溶质对流通量:J=qc=vOc，q：土c壤水通量，c：土壤水溶液的质量浓度,v=q/6土壤水溶液的平均孔隙流速。对流作用下溶质的穿透时间：t=L/v=LO/qL分子扩散(Moleculardiffusion)：不同浓度的溶质间在浓度梯度的作用下产生的质量交换原因：分子不规则热运动(布朗运动)趋势：高浓度f低浓度，溶质的分子扩散符合Fick第一定律:土壤水溶液：JFick第一定律:土壤水溶液：JJd分子扩散通量；Ds扩散系数溶质扩散系数Ds:随土壤含水率的减小而减小机械弥散(mechanicaldispersion)：又称对流弥散(convectiondispersion)，由于土壤水流和土壤孔隙结构相互作用使溶质在土壤孔隙中运移和分散的过程。产生原因：土壤孔隙不均匀，流速方向和大小不同，使得溶质在水流运动中分散Qc溶质的机械弥散通量：J=-D(V)，机械弥散系数Dh(v)：Dh(v)二入|v|hhQZ经验常数(dispersivity)入与土壤质地、结构有关；机械弥散系数还与空间尺度有关水动力弥散(hydrodynamicdispersion)：综合考虑分子扩散与机械弥散溶质的水动力弥散通量：J=—[D(°)+D(V)1=-D(^0)]shshshQZOZ水动力弥散系数(综合扩散一弥散系数)：d(v,e)，当土壤水流速度较大时，机械弥散作用超过扩散作用；土壤溶液静止时，只有分子扩散作用土壤溶质通量：J=J+Jc=qc土壤溶质通量：J=J+Jc=qc-D(如)shshOcQZ溶质运移连续方程(一维，不考虑源汇项)0(0c)_QJQ^'QZ一维溶质运移方程：根据水流连续方程:Qc1Q(qc)QZQz对流一弥散方程：CDE土壤水分运动对溶质运移的影响Q(土壤水分运动对溶质运移的影响Q(0cQ+Q丿水流通量q影响溶质对流；土壤含水率/、水流通量q(土壤平均孔隙流速v=q/&)影响溶质的水动力弥散；土壤溶质的源汇项、动态储存一般与有关。溶质对土壤水分运动的影响溶质势：屮=-cRT/(pwgM)水流通量：_-K(0)『Q屮m+fQ屮s+*sCQzo丿f选择系数：=1:完全选择，溶质不能通过半透膜)；=0：溶质能通过，溶质势对水流无影响；粘土层存在一定的半透膜作用，非粘土中一般可以忽略溶

质势的影响。溶质的存在导致土壤导水率的降低。溶质可能影响密度，产生密度流（如海水入侵）土壤比热容（Specificheatcapacity）：单位体积（质量）的土壤，温度升降1°C所吸收或释放的热量，称为体积（质量）比热容，比热容是土壤的一种广度性质，具有可加性,是土壤各组分（矿物质、有机质、水、冰、空气比热容之和。热导率：热量传输能力，数值上等于温度梯度为时的热通量。：q=-K▽T，q:热通hhh量，Wcm-2，K:热导率，Wcm-iC-ih热扩散率：D=K/c，Fourier热传导定律：q=-K^T=-CD▽T，D:cm2/s随着含水率的增加；热导率的提高较比热容的提高要快，'但当含水率较高时情况恰好相反。因此含水率0较小时，随e的增大而增大，含水率e较大时，随e的增大而减小。土壤水分运动对热量传输的影响：土壤水分状况影响土壤热特性参数；土壤水分相变（蒸发与凝结、冻结与融化等）产生热源、汇；地表基质势影响水汽压，从而影响地表热量平衡土壤温度对水分运动的影响温度影响水的物理化学性质，从而引起土壤水分运动参数、基质势的变化;温度势所产生的水分运动。地表能量平衡：R=C+Xe+GR:净辐射C:显热消耗，九E：蒸发潜热；G：地表的nssnss热通量早期对各种水分传输过程单独研究，概念不统一后来，用水势的概念综合研究水分传输,将土壤一植物一大气视为连续体，Philip于1966年提出SPAC的概念，土壤一植物一大气连续体：Soil-Plant-AtmosphereContinuum。SPAC系统的三个重要环节：根系吸水过程，水在植物体内的输运过程，植物体内的水分向大气扩散过程。蒸腾（Transpiration：主要通过叶片气孔散失。水分在气孔腔和细胞间隙的叶肉表面蒸发（水面蒸发）水汽达到饱和；水汽由气孔腔经气孔向叶面气层扩散；水汽从叶面气层向大气扩散。腾发的影响因素：能量：辐射，潜热消耗L=2.5MJ/kg;水汽输送：风速、温度梯度、湿度梯度；其它影响因素：植物生长状况、土壤供水状况直接测定法：根据实测水分剖面，利用零定位）通量法计算；称重式蒸渗仪；涡度相关R=LE+H+G=LE+nLETR=LE+H+G=LE+nLETLE涡度相关是指某种物质的垂直通量，即这种物质的浓度与其垂直速度的协方差。涡度相关法已经成为直接测定大气与植物群落气体交换通量的通用标准方法。能量平衡与紊流扩散理论结合，只需一个高度的气象资料oPenman公式（计算潜在腾发量）水面f湿润表面LEtp发量）水面f湿润表面LEtp丄EPeku(；丫naE2A-e片a22s2a2学习好资料学习好资料欢迎下载22+1PI其中Ea常称为干燥力，反映了空气湿