《第二章土壤水分运动基本方程2》由会员分享，可在线阅读，更多相关《第二章土壤水分运动基本方程2（24页珍藏版）》请在金锄头文库上搜索。
1、.第二章 土壤水分运动基本方程如前所述，达西定律是由达西（Darcy，Henry 1856）通过饱和砂柱渗透试验得出，后由Richards（1931）将其扩伸至非饱和水流中，并规定导水率为土壤负压h的函数，即 （2-2-1）式中：为水势梯度； k（h）为导水率，是土壤负压h的函数； q 为水流通量或流速。Richards方程垂向一维方程为注意：H=hz，垂直坐标向上为“+”；向下时为“”。由于k（h）受滞后影响较大，上式仅适用于单纯的吸湿或脱湿过程。若将导水率作为容积含水率函数，即以k（）代替人k（h），则可避免滞后作用的影响。一般说来达西定律对饱和与非饱和水流均可适用，即水流通量与势能梯度成
2、正比。但在饱和土壤中，压力为正值，其总水头包括了由该点在地下水面以下深度来确定的静水压力（正值）和相对于基准面高度来确定的位置水头，总水头为压力水头和位置水头之和，水由总水头高处向低处流动。在非饱和土壤中，基质势为负值，土水势在不考虑溶质势、温度势及气压势时，只包括重力势和基质势。因此，总水头常以负压水头和位置水头之和来表示。一维Richards方程的几种形式：根据（K=CD）得： 精品.第一节 直角坐标系中土壤水分运动基本方程一、基本方程的推导土壤水分运动一般遵循达西定律，且符合质量守恒的连续性原理。土壤水分运动基本方程可通过达西定律和连续方程进行推导。如图2-2-1所示，从土壤中取出微分单
3、元体abcdefgh，其体积为，由于该立方体很小，在各个面上的每一点流速可以看成是相等的，设其流速为，在tt+t时段内，流入立方体的质量为（3个面流入）： （2-2-2）流出立方体的质量为（3个面流出）： （2-2-3）式中：水的密度；分别表示微分体x、y、z方向长度；，分别表示水流经微分体后，其流速在x、y、z方向的变化值。由式(2一22）、式（223）之差可求得流入和流出立方体的质量差： （224）设为立方体内土壤含水率，则在t时间内立方体内质量变化又可写为 （225）根据质量平衡原理（流入量流出量储存量变化量），式（324）、式（325）应相等，即精品. （2-2-6）根据达西定律得：，
4、 （2-2-7）式中k（）土壤水力传导度，为含水率的函数；H总土水势，为基质势与重力势之和（Hhz）。因此，式（226）可以写作以下形式： （2-2-8）上式可以简写为 （2-2-9）式（228）或式（229）为土壤水分运动基本方程。在饱和土壤中，含水量和基质势均为常量。水力传导度也为常量，常称渗透系数，则方程（228）可写为 （2210）或写作 （2210） （2211）式中：2拉普拉斯算子。式（2210）或式（2210）为饱和土壤水流的拉普拉斯方程。二、基本方程的不同形式为运用基本方程分析各种实际问题的方便，可将基本方程改写为多种表达形式。为简便起见，以下均以一维垂向土壤水分运动为例，给出
5、基本方程的不同表达形式。（一）以含水率为变量的基本方程由式（228）可得一维垂向土壤水分运动的基本方程为 （2212）式中：H总土水势； z为水流方向坐标，取z向上为正。因为H=h十z，所以上式可写作精品. （2213）式(2213）为以为变量的基本方程，将代入式（2213）得：令，则式（2213）可以写成（一维垂向土壤水分运动方程）： （2214）在水平运动的情况下，重力项等于0，所以，其形式与Fick扩散定律相同。式（2214）具有扩散方程的形式，故将D（）称为扩散度。 （2214）Fick定律：自由水中溶质的分子扩散通量符合Fick定律：式中：J为溶质的扩散通量； D为溶质的扩散系数；
6、为溶质的浓度梯度。（二）以基质势h为变量的基本方程由于 ，则式（2214）可以写成： （2215）式中：c（h）比水容量（也称容水度），c（h）=，表示单位基质势变化时含水率变化。（三）以参数v为因变量的基本方程采用Kirchhoff变换，令则 精品.由式（2-2-15）得： （2216）式中hc土壤的进气值，即土壤含水率开始小于饱和含水率时的负压值。另外，；在非饱和区： 在饱和区： 且因为 ，所以 ；则方程式（2216）为：（四）以位置坐标z为变量的土壤水运动方程以z为变量，则z为、t的函数，z（，t）为未知函数。已知=（z，t），当处，可以解出z= z（，t），即14对z，t分别求导数：，
7、精品.于是 及将以上式子代入方程（2214）得： （2217）（五）以参数u为因变量的土壤水运动方程定义式中：初始含水率； ； 饱和含水率。由式（2-2-14）得： 将代入上式得： 所以 （2-2-18）以上各式中式（2214）、式（2215）是二种经常采用的形式，形式的选定取决于要解决问题的边界条件和初始条件。精品.以含水率为因变量的基本方程常用于求解均质土层或全剖面为非饱和流动问题，这种方程形式对于层状土壤或求解饱和非饱和流问题不适用；以负压水头h为因变量的基本方程是应用较多的一种形式，可适用于饱和非饱和水流求解及层状土壤的水分运动分析计算，但由于非饱和土壤水的导水率k（h）及容水度c(h
8、)，受滞后影响较大，计算中参数选取不当会造成较大误差；以v，u为因变量基本方程实际上分别相当于以负压水头h和含水率为因变量的基本方程，在某些情况下由于经代换后方程较为简单，易于求解；以坐标为因变量的基本方程根据定解条件需要求解较简单的土壤水分运动问题。以上为直角坐标系中土壤水分运动的基本方程，求解某些土壤水分运动问题时，采用柱坐标系可能更方便。第二节 柱坐标系中土壤水分运动基本方程在推导柱坐标系中的基本方程时，方法同直角坐标系，同样可用达西定律与连续方程相结合的方法导出。若以z轴为轴的柱坐标系，根据达西定律，在此坐标系中可表示为：式中：r、z分别表示柱半径，角坐标和垂直坐标（图222） q r
9、、q、qz相应于r、z三个方向的通量； H总水势。下面利用质量守恒来推导连续方程。t时段内，在r方向的流入量为zt，流出量为精品.，则流入与流出量之差（忽略高阶无穷小量）为 （2-2-19）同理，在方向流入流出量之差为 （2-2-20）在z方向土壤水分流入流出量之差为 （2-2-21）上述三个方向流入和流出单元体的水量差总计为 （2-2-22）单元体体积应为 ，略去高阶无穷小量后为rzt，在t时间内单元体内水分增量为 （2-2-23）根据质量守恒原理武（3222）应与式（3223）相等，即 （2-2-24）式（2224）为柱坐标系中土壤水分运动的连续方程。将式（2218）代入上式，即得柱坐标系
10、中土壤水分运动基本方程： （2-2-25）以总水势H=h+z，水容度c=，以及导水率k（），扩散度D（）等代入，基本方程可表示为 （2-2-26）对于平面轴时称问题，上式可改写为 （2-2-26）精品.同理可推得以x（或y）轴为轴的柱坐标系的基本方程： （2-2-27）关于球坐标系中的基本方程应用较少，推导方法同上，这里不再论述。第三节 土壤水分运动基本方程的定解条件土壤水运动基本方程的定解条件包括初始条件和边界条件，为了简单起见，将以直角坐标系中基本方程常用形式为例进行论述。（一）初始条件相应于前式（1）、式（15）的初始条件分别以下式表达： （2-2-28） （2-2-29）脚标“i”表示初始已知量。初始条件：t=0时剖面上、h的分布已知。（二）边界条件边界条件一般有一类边界、二类边界、三类边界三种。1一类边界条件（变量已知的边界1）对干式（2214）、式（2215）的一类边界的表达式为 （2-2-30） （2-2-31）脚标“0”均表示一类边界上的值；z0为一类边界的坐标。在一维垂向土壤水分运动中，一类边界的情况发生在压力入渗（地表形成水层）时，地表含水率达到饱和含水率，或当强烈蒸发时，表土达到风干土含水率的情况。2二类边界条件（边界2上水流通量已知的情况）相应于式（221
