土壤水动力学基本方程.ppt

关于土壤水动力学基本方程现在学习的是第1页，共78页2.1 2.1 土壤水流概述土壤水流概述2.1.1 2.1.1 毛细管中的毛细管中的层流层流运动运动yuFddyudd单位面积切应力Poiseuill(普氏）定律：取流股流股圆柱面压力为212pppyp作用流股侧面积的切应力为yL22yLpyudd负号出现是由于 随 的增长而减少。uy液体质点作有条不紊的线性运动，彼此互不混杂。图2-1 圆管中的层流运动现在学习的是第2页，共78页2.1 2.1 土壤水流概述土壤水流概述2.1.1 2.1.1 毛细管中的层流运动毛细管中的层流运动由此2yLpdydu积分得 42yCLpyu积分常数 ： 求得0uRy4RC2C得 224yRLpyu在 处，0yLpRu42max单位时段通过细管的流量LpRyyuuQR8240dd通过细管的平均流速paRLpRu228细管形状参数压力梯度现在学习的是第3页，共78页2.1 2.1 土壤水流概述土壤水流概述2.1.2 2.1.2 土壤水流简化模型土壤水流简化模型为什么需要简化模型Q简化QJw平均“流速”实际流量横截面面积通量（L/T) ：水流并不是在整个横截面 上进行；真实水流通道大于表观长度。wJ压力势或基质势与实际土壤孔隙中的势能相等。现在学习的是第4页，共78页2.2 2.2 饱和土壤水运动的达西定律饱和土壤水运动的达西定律 饱和流：水分充满土壤孔隙的水流。饱和流：水分充满土壤孔隙的水流。2.2.1 Darcy 定律wJ通量 ：QJw21HH 1-1和2-2之间产生一个力F2121zzHHzHFDarcy研究：wJF 引进一个比例常数 ，称作土壤导水率（soil water conductivity)sKQ水力传导度现在学习的是第5页，共78页2.2.2 Darcy 2.2.2 Darcy 定律的适用范围定律的适用范围Darcy定律只适应土壤水流为层流的情况。水流的两种流动形态对颗粒极细的土壤，如粘土，水流表现出非Newton流(Bingham流)性质。屈服点临界梯度实际上，Darcy定律在绝大多数情况下可应用于土壤水流计算，只是在粗砂或粘质土壤情况下要注意Darcy定律的限定。2.2 2.2 饱和土壤水运动的达西定律饱和土壤水运动的达西定律 现在学习的是第6页，共78页2.2.3 2.2.3 导水率导水率导水率 综合反映了多孔介质和流体的某些物理性质。sK（1） 计算导水率 的公式sK（2） 实验室测定实验室测定仪器的基本原理与Darcy原始实验装置大致类似，只是多采用自下而上的水流。定水头法：Why?HQKs2.2 2.2 饱和土壤水动的达西定律饱和土壤水动的达西定律 现在学习的是第7页，共78页2.2 2.2 饱和土壤水运动饱和土壤水运动 的达西定律的达西定律 2.2.3 2.2.3 导水率导水率（2） 实验室测定变水头法：示意图br2R2L1t 时刻立管水柱深 ， 时刻为 ，时刻为 。通量1b2b t tb2t tRtbrtRtbJdddrd2w222土柱底部000111Hpz土柱顶部 tbLHtbpLz222 LtbLKzzHHKtRtbrJss212122ddw现在学习的是第8页，共78页2.22.2饱和土壤水运动的达西定律饱和土壤水运动的达西定律 2.2.3 2.2.3 导水率导水率（2） 实验室测定变水头法：整理得： tKtbLtbRrLsdd22积分得： 212122ttsbbdtKtbLtbRrLd 122212lnttRLbLbLrKs变更立管水深，可求得一系列 ，在误差范围内求其平均值。sK无论定水头法还是变水头法必须考虑土样的代表性。现在学习的是第9页，共78页2.2.3 2.2.3 导水率导水率（3） 田间现场测定双环法：tWKs入渗量(cm3)测定时段内环横截面积双环法一般只能测定地表土壤导水率，用其他仪器，如Guelph仪可测其他深度土壤的导水率。导水率大致范围6cm/d 很小 616cm/d 低 1640cm/d 中 40100cm/d 高 100cm/d 很高 外环的作用？2.22.2饱和土壤水运动的达西定律饱和土壤水运动的达西定律 现在学习的是第10页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 绝大多数田间和植物根区的土壤水流过程都处绝大多数田间和植物根区的土壤水流过程都处在非饱和状态。非饱和流研究为土壤物理学最在非饱和状态。非饱和流研究为土壤物理学最活跃的研究领域之一。活跃的研究领域之一。2.3.1 非饱和流与饱和流的比较（1）土壤水流驱动力不同（2） 导水率的差别（3） 土壤孔隙对饱和水流和非饱和水流影响的差别现在学习的是第11页，共78页非饱和流与饱和流的比较：共同之处：都服从热力学第二定律，都是从水势高的地方向水势低的地方运动。不同之处：土壤水流的驱动力不同。 饱和流的驱动力是重力势和压力势； 非饱和流的是重力势和基质势。导水率差异 非饱和导水率远低于饱和导水率；当基质势从0降低到-100kpa时，导水率可降低几个数量级，只相当于饱和导水率的十万分之一。土壤空隙的影响土壤。在高吸力下，粘土的非饱和导水率比砂土高。现在学习的是第12页，共78页饱和流与非饱和流（饱和流与非饱和流（1） 饱和流：土壤空隙全部充满水时的流动。饱和流：土壤空隙全部充满水时的流动。n发生情形发生情形： 1. 大量持续降水和稻田淹灌时，垂直向下；大量持续降水和稻田淹灌时，垂直向下； 2. 地下泉水涌出，垂直向上；地下泉水涌出，垂直向上； 3. 平原水库库底周围，水平方向。平原水库库底周围，水平方向。n推动力：推动力：重力势梯度和压力势梯度重力势梯度和压力势梯度n影响因素：影响因素：soil texture and structure现在学习的是第13页，共78页饱和流与非饱和流（饱和流与非饱和流（2） 非饱和流：土壤空隙未全部充满水时的流动。非饱和流：土壤空隙未全部充满水时的流动。n发生情形：发生情形：大多数情况大多数情况n推动力：推动力：基模势梯度和重力势梯度基模势梯度和重力势梯度现在学习的是第14页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.2 Buckingham-Darcy2.3.2 Buckingham-Darcy通量定律通量定律Buckingham对Darcy定律描述土壤非饱和流提出修正的两个基本假设： 土壤非饱和流驱动力是基质势与重力势之和的梯度； 非饱和土壤水流的导水率是土壤含水量或基质势的函数。以水势头为单位，Buckingham-Darcy通量定律可写成： 1zhhKzzhhKzHhKJw符号相同,向上为正难点现在学习的是第15页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.2 Buckingham-Darcy2.3.2 Buckingham-Darcy通量定律通量定律Buckingham-Darcy通量定律也可写成： 1zhhKzzhhKzHhKJw问题：两种写法是否会影响计算结果？基质势 是土壤深度 和时间 的函数，所以用偏微分 表示：hztzhztzhtzzhzhzhzt,lim0偏微分方程用以对非稳态流的数学描述，如是稳态流，上式变为常微分方程。符号相反,向下为正现在学习的是第16页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.32.3.3非饱和导水率的数学表达非饱和导水率的数学表达在饱和含水量附近，砂质土壤的导水率高于粘质土壤。随着吸力的发展，砂质土壤大孔隙排空，流径增加，由此其导水率低于粘质土壤。几个非饱和导水率经验公式： mSbaSK 1msCSKSK mssKK rsrMsMKK112121现在学习的是第17页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.32.3.3非饱和导水率的数学表达非饱和导水率的数学表达毛细管模型将土壤概化为一束不同数量和尺寸的毛细管，并假定毛管束的特征曲线与所代表的实际土壤的特征曲线相同。以下几方面与实际土壤不同： 毛细管有相同长度； 水流边界与实际土壤不同； 毛细管半径完全控制着水膜厚度； 水流是稳态的。 每根毛细管都是连通的；现在学习的是第18页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.32.3.3非饱和导水率的数学表达非饱和导水率的数学表达由普氏定律：LpRQ84HgpwcjjLHgRQ84w w半径为 毛细管的流量jR模型的毛细管长度，实际土壤表观长度为 ， 。LcLL通过模型的总流量MjjjcjMjjTRNLHQNQ1418 8w wg g毛管束中半径为 的毛管数量。jR毛管束中不同尺寸毛细管的数量。通过毛管束的通量MjjjcTRnLHgQJ148ww ，毛管中半径为 的毛细管单位面积的数量。jjNn jR现在学习的是第19页，共78页用于计算非饱和导水率的毛细管模型2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.32.3.3非饱和导水率的数学表达非饱和导水率的数学表达如图，将实际土壤水特征曲线 分成等宽为 的若干份。 h于是有,2,21sshhhh假定：当 时， 的毛管都排空。jhh jRR jjghRw2毛管半径由确定。如假定毛管横截面面积 ，则单位长度毛管排水量为：2R211211RnRn如此，2jjRn代入毛管通量计算式MjjjcRnLHgJ148ww现在学习的是第20页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.32.3.3非饱和导水率的数学表达非饱和导水率的数学表达zHLLhgLHLLhgLHLLhgLHRgRRLHgJcMjjcMjjcMjjcMjjMjjjc12212212212142120121288w ww ww ww ww ww w设MjjshgK12212w为弯曲度。cLL现在学习的是第21页，共78页2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 2.3.32.3.3非饱和导水率的数学表达非饱和导水率的数学表达含水量为 ，s最大半径为 的毛管排空。1RMjjshgK22212w对一般情况1, 2 , 112122MihgiKMijjsw又MjjMijjMjjMijjsshhhghgKiK1212122122111212ww212111,2,11Mj ijssMjjhKiKiMh 例题2.1现在学习的是第22页，共78页 一个重要概念：通量不随时间变化的土壤水流称做稳态流，或恒定流；通量随时间变化的土壤水流称作非稳态流，或非恒定流，或瞬态流。讨论：如Darcy实验示意图。当土壤水流达到稳态水流后，哪些土壤水运动要素不随时间发生变化？水位在这种情况下，土壤水流是否能达到稳态流？2.3.4 稳定状态下的非饱和流问题稳定状态下的非饱和流问题2.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第23页，共78页2.3.4 稳定状态下的非饱和流问题稳定状态下的非饱和流问题非饱和流在大多数情况下都处在非稳定状态，只有在一些理想条件下才可近似将非饱和流看作稳态流。当基质势差 在土柱两端（ ）保持不变，则土柱水流最终达到稳态。Buckingham-Darcy 定律可写成：12hhh12zzL 1zhhKJddwBuckingham-Darcy通量定律 1zhhKJw w 为常微分，因为只取决于 ，而与 无关。zhddzt2.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第24页，共78页（1） Darcy定律的积分形式定律的积分形式 1zhhKJddw 1zhhKJddw zhhKJddw 1 zhKJhd dd dw w1当 2211,zhhzhh 2121211zzzhKJhzzhhddw2.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第25页，共78页（2） 地下水稳定时的土壤蒸发地下水稳定时的土壤蒸发虽然田间水分蒸发不是一个稳态过程，但稳定的地下水位向裸地土壤表面蒸发，在一段时间，大气蒸发条件相对稳定，可近似看作稳态蒸发。例题2.2（3） 稳态向下的土壤水流向下的稳态水流在田间几乎不会出现，但在某些情况下，如频繁灌水或降雨，可近似地将向下的田间水流看作向下运动的稳态流问题。例题2.32.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第26页，共78页（3） 稳态向下的土壤水流稳态向下的土壤水流对任何通量 ，当地下水位不太浅时，接近地表的基质势趋于常数。i重要结论：当水流以常量下渗时， 趋于0，水流只在重力梯度下运动。可得：zh dd hKJw w（4） 非饱和导水率的稳态实验室测定例题2.42.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第27页，共78页2.3.52.3.5 测定非饱和导水率的瞬时剖面法测定非饱和导水率的瞬时剖面法实验室测定：如图，土柱进水底面为参照面，向上为正。 1zSKzzSKzzhKJw 1zSwJKzwJtwddJzt zJzJztzwwd002.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第28页，共78页2.3.5 测定非饱和导水率的瞬时剖面法 zztJzJ0wwd0is在每次测定间隔 确定情况下：tztztztz1d0ztzd0实际是图中两条曲线的面积。由供水的Mariotte瓶刻度读出2.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第29页，共78页2.3.5 测定非饱和导水率的瞬时剖面法测定非饱和导水率的瞬时剖面法田间测定：田间测定一般不会有实验室测定那样的通过Mariotte瓶测得的 ，因此计算水流通量的公式写成更一般形式。0wJNiztJztJJNiizziiiii, 2 , 11111wwwdNi, 2 , 1式中， 为测定点的编号，即张力计（或TDR）埋设根数。 瞬时剖面法测定关键是要事先确定一个截面的通量，可通过地表辅膜或零通量面法确定。 1zSwJK由 确定 ，在得到足够数据后，可选线型拟合。 K与毛管模型法比较2.3非饱和土壤水运动的达西定律非饱和土壤水运动的达西定律 现在学习的是第30页，共78页2.4 2.4 非饱和土壤水运动的基本方程非饱和土壤水运动的基本方程zyx,如图，取一块长，宽，高为 微小土体。根据物质守恒原理。对一维垂直流：在 时段进入土体的水量=t在 时段离开土体的水量+t在 时段储存在土体中水容量的增量+t在 时段由植物根系吸收而失去的水量t2.4.1质量守恒与基本方程的推导现在学习的是第31页，共78页2.4.12.4.1质量守恒与基本方程的推导质量守恒与基本方程的推导进入土体的水容量=tyxttzyxJ21,w离开土体的水容量=tyxttzzyxJJ21,ww储存在土体中水容量的增量=zyxtzzyxttzzyx,21,21,土体中由植物吸收的水容量=tzyxrw源汇项现在学习的是第32页，共78页2.4.12.4.1质量守恒与基本方程的推导质量守恒与基本方程的推导由物质守恒原理可得：tyxttzzyxJtyxttzyxJ21,21,wwzyxtzzyxzyxttzzyx,21,21,tzyxrw两边同除以tzyxzttzyxJttzzyxJ21,21,ww0,21,21,wrttzzyxttzzyx2.4 非饱和土壤水运动的基本方程现在学习的是第33页，共78页2.4.12.4.1质量守恒与基本方程的推导质量守恒与基本方程的推导当,0,0tz有：0wwzrtJ对更一般的三维情况，有：0wwwwrtzJyJxJzyxkJjJiJJzyxwwwwkji,zyx,分别是方向三个单位矢量。当土体中无植物根系存在， 。0wr以上推导两个基本假设：水是不可压缩的； 土壤基质在水流过程中保持不变。2.4 非饱和土壤水运动的基本方程现在学习的是第34页，共78页2.4.22.4.2基本方程的各种形式基本方程的各种形式（1 1）以含水率）以含水率为因变量的基本方程为因变量的基本方程（2 2）以基质势）以基质势h h为因变量的基本方程为因变量的基本方程（3 3）RichardsRichards方程的柱坐标形式方程的柱坐标形式2.4 非饱和土壤水运动的基本方程现在学习的是第35页，共78页2.4.2基本方程的各种形式基本方程的各种形式(1) Richards方程的含水量形式（ 方程）Darcy定律的 变量形式： 可以直接写成 的函数：K KhK 由微分变换zhzhdd比水容量 hCCdd KzhKzhhKJddw1 KzDKzCK代入连续性方程： zKzDzt现在学习的是第36页，共78页 zKzDzt 方程，又称扩散型方程。上式是一个二阶非线性偏微分方程。非线性是指所求函数 又是 和 的变量。DK注意：土壤水运动不是扩散运动， 的引入只是一种数学处理方式。 D扩散运动的驱动力是浓度梯度，而土壤水运动的驱动力是土水势梯度。扩散型方程的优点是，与 相比， 的变化范围要小得多；其缺点是扩散型方程只能用在均质土壤剖面上。 K DWhy?(1) Richards方程的含水量形式（ 方程）2.4.2基本方程的各种形式基本方程的各种形式现在学习的是第37页，共78页(2) Richards方程的基质势形式（ 方程）h由导数性质，有： ddhhC hthtt也是比水容量， hhCdd代入连续性方程： zhKzhhKzthhC上式称作Richards方程的基质势形式，又称 方程。也是一个二阶非线性偏微分方程。从理论上讲，可以用在非均质土壤剖面的水流问题。h无论 方程还是 方程，一般都忽略土壤水的滞后作用。h2.4.2基本方程的各种形式基本方程的各种形式现在学习的是第38页，共78页(3) Richards方程的柱坐标形式方程的柱坐标形式如图，以 为垂直轴的Darcy定律可表示为：z zHKJHKrJrHKJzwwwr1基本假设：水是不可压缩的；土壤基质在水流过程中不变形；土壤是各向同性的（isotropic）。2.4.2基本方程的各种形式现在学习的是第39页，共78页 方向：r在 时段，流入流出的水量差：ttzrrJrJrrww在 时段，流入单元体的水量：ttzrJr rw w入流面积流出水量：tzrrrrJJrrww出流面积(3) Richards方程的柱坐标形式方程的柱坐标形式2.4.2基本方程的各种形式现在学习的是第40页，共78页 方向： 时段流入水量：ttzrJw入流面积流出水量：tzrJJww出流面积在 时段，流入流出的水量差：ttzrJw 方向：z在 时段，流入流出的水量差：ttzrrJzzw面积(3) Richards方程的柱坐标形式方程的柱坐标形式2.4.2基本方程的各种形式现在学习的是第41页，共78页在 时段流入和流出单元体总的水量差：ttzrzJrJrJrJzrrwwww单元体体积为：zrrr2略去高阶无穷小量后为：zrr单元体内水分增量为：tzrrt由物质守恒原理，得连续性方程：zJJrrJJrtzrrwwww11(3) Richards方程的柱坐标形式方程的柱坐标形式2.4.2基本方程的各种形式现在学习的是第42页，共78页将柱坐标系Darcy定律代入连续性方程，得： zHKJHKrJrHKJzwwwr1 zHKzHKrrHKrrHrKt21柱坐标Richards方程常写成： zKzhKzhKrrhrKrrt211(3) Richards方程的柱坐标形式方程的柱坐标形式2.4.2基本方程的各种形式现在学习的是第43页，共78页同样，柱坐标下的Richards方程也可写出 方程和 方程：h zKzDzDrrrDrrt211 方程： 方程：h zhKzhhKzhhKrrhhrKrrthhC211对 在 方向无变化，同时具有轴对称特点，方程可以简化为：z rhrKrrt1或 rrDrrt1 例题2.5(3) Richards方程的柱坐标形式方程的柱坐标形式2.4.2基本方程的各种形式现在学习的是第44页，共78页2.4.3 土壤水运动方程的定解条件 基本方程定解问题 初始条件 定解条件 边界条件(1) 初始条件（initial condition） 表示所研究问题的初始状态：梁昆淼，数学物理方法，人民教育出版社，1979对于方程：对于 方程：用角标“0”表示初始已知量现在学习的是第45页，共78页(2) 边界条件（boundary conditions） 边界条件一般分为三种：一、二、三类边界。 一类边界条件(Dirichlet)：变量已知边界对于方程：对于 方程：用角标“1”表示第一类边界上的值。 在一维垂向土壤水分运动中，一类边界的情况发生：压力入渗（地表有薄层积水）；(1) 强烈蒸发（表土达到风干含水率）。现在学习的是第46页，共78页 二类边界条件(Neumann)：水流通量已知边界对于方程：对于 方程：用角标“ ”表示第二类边界。 在一维垂向土壤水分运动中，二类边界的情况发生：降雨、灌水入渗、蒸发强度已知的边界；(1) 不透水边界和无蒸发入渗的边界，此时 =0 。现在学习的是第47页，共78页 三类边界条件(Cauchy)：水流通量随边界 上的变 量或 而变化的情况。对于方程：用角标“ ”表示第三类边界。321fzf通式：时，则为第三类边界。时，则为第二类边界。时，则为第一类边界。000321其中，fVariable(, )对于 方程：现在学习的是第48页，共78页1. 试推导垂直一维土壤水流的基本方程，并写出如右图所示情况AB剖面土壤水分运移的定解问题。（地表入渗强度为q(t)，B点地下水位保持不变，AB剖面土壤均质）。B Azq(t)第2章 作业现在学习的是第49页，共78页2.52.5土壤水运动其它求解方法土壤水运动其它求解方法-通量法通量法tyxzqqqxyz =简化为zqtz 积分得 zzq zq zdzt21,zzzzz t dzz t dzQ z Q z=对上式积分得 确定某一断面处的通量，主要应用达西定律，其方法有零通量法，表面通量法和定位通量法，统称为土壤水分运动通量法。定义现在学习的是第50页，共78页零通量面（零通量面（ZFP-Zero Flux Plane）在土壤较为潮湿情况下，剖面上部由于蒸发水分向上运动；而剖面下部水分在重力作用下向下运动。在水分向上和向下运动交接处，必然会出现零通量面。由Darcy定律： zHKzHhKJwzHZFP0如图，ZFP处切线斜率垂直于 轴，即 H0zHZFP是一个复杂问题，有许多这方面的研究。例题2.62.5土壤水运动其它求解方法土壤水运动其它求解方法-通量法通量法现在学习的是第51页，共78页2.5.12.5.1零通量面与零通量法零通量面与零通量法2.5土壤水运动其它求解方法-通量法零通量面的类型单一聚合型零通量面(该情况下土壤水分由上下两侧向零通量面迁移 ) 单一发散型零通量面 (该情况下土壤水分由零通量面向上下两侧迁移 )具有多个零通量面 (发生在间隔降雨，入渗和蒸发交替出现的情况下 )现在学习的是第52页，共78页2.52.5土壤水运动其它求解方法土壤水运动其它求解方法-通量法通量法2.5.1零通量面与零通量法(a)单一聚合型(b)单一发散型(c)多零通量面现在学习的是第53页，共78页sQ零通量面位置不变时水势与含水率的分布 0012,HHzzz t dzz t dz=sQ 即图中abcd的面积。当含水率减少时， 0，表示通量向上，土壤水分蒸发；反之， 0， 表示通量向下，水分向下层土壤入渗。sQsQ2.52.5土壤水运动其它求解方法土壤水运动其它求解方法-通量法通量法2.5.1零通量面与零通量法现在学习的是第54页，共78页00g2100Q,zzz t dzz t dz 即图中ade的面积。当土壤含水量减少时， ，表示通量向下，即潜水接受补给；反之， 0，表示通量向上，意味着蒸发时潜水有消耗。 gQgQgQ2.52.5土壤水运动其它求解方法土壤水运动其它求解方法-通量法通量法2.5.1零通量面与零通量法现在学习的是第55页，共78页2.5.2.1表面通量法 表面通量法是以地表处的入渗量或蒸发量作为已知条件的。 21,HHszzQ zQz t dzz t dz2.5.2.2定位通量法211mmmq zKz 21,zzzzQ zQ zz t dzz t dz2.52.5土壤水运动其它求解方法土壤水运动其它求解方法-通量法通量法2.5.2表面通量法现在学习的是第56页，共78页2.6 2.6 土壤中溶质运移的基本方程土壤中溶质运移的基本方程 土壤中的溶质处在一个物理，化学和生物土壤中的溶质处在一个物理，化学和生物的的相互联系和连续变化的系统中。本节相互联系和连续变化的系统中。本节侧重分析溶质运移的侧重分析溶质运移的物理作用物理作用，并认为，并认为溶质运移主要是通过溶质运移主要是通过对流和水动力弥散对流和水动力弥散两种作用实现的。两种作用实现的。现在学习的是第57页，共78页2.6 土壤中溶质运移的基本方程土壤中溶质运移的基本方程2.6.1 土壤溶质的形成形式土壤溶质的形成形式 严格而言，自然界中的土壤水分都是以溶液形式存在的。以溶液中的溶质为考察对象，考察某一深度土体土壤溶质的形成形式：进入溶液的溶质：1.以降水，灌溉和径流形式从土壤表层进入土壤；2.原先在土壤颗粒上的物质会被土壤溶液所溶解或解吸附而进入土壤溶液； 3.土壤空气中的一些物质也会进入土壤溶液。 离开溶液的溶质：1.溶质会从研究土体的下界面离开研究土体； 2.溶液中一部分溶质会沉降到土壤颗粒上，也可能被土壤颗粒吸附而脱离土壤溶液； 3.一些溶质挥发脱离土壤溶液而进入土壤空气中。还有一个重要形成过程是溶液中的溶质与周围环境中的某些物质或自身发生各种化学或生物反应，使土壤溶液中的溶质增加或减少。 现在学习的是第58页，共78页2.6 土壤中溶质运移的基本方程土壤中溶质运移的基本方程2.6.2 溶质质量守恒公式溶质质量守恒公式取一个六面单元土体，单元体的体积 ，从时刻 到时刻 ，有：zyxVttt在 时段进入单元体的溶质质量=在 时段离开单元体溶质的质量+在 时段单元体储存的溶质质量的增加+在时段 单元体由于化学、生物反应或植物吸收从单元 体消失的溶质tttt如果我们假设溶质只在 方向流动，则上式四项可写成： z现在学习的是第59页，共78页2.6 土壤中溶质运移的基本方程土壤中溶质运移的基本方程tyxttzyxJ21,s styxttzzyxJ21,s szyxtzzyxCttzzyxC,2121T TT Ttzyxttzzyxr2121,s s总的溶质通量（单位时间单位面积流过的溶质质量） 总的溶质浓度（单位容积土壤的溶质质量） 单位容积的溶质的反应速率（单位时间单位土壤容积溶质的损失） 上式除以 ，整理可得： tzyx0212121212121ttzzyxrttzzyxCttzzyxCzttzyxJttzzyxJ,s sT TT Ts ss s取 ， ，得0z0t0s ss sT TrzJtC上式称作一维溶质守恒方程。 现在学习的是第60页，共78页2.6.3溶质运移的对流和水动力弥散2.6 2.6 土壤中溶质运移的基本方程土壤中溶质运移的基本方程2.6.3.1溶质的对流运移 对流是指在土壤水分运动的过程中，同时携带着溶质运移。单位时间内通过土壤单位横截面积的溶质质量称为溶质通量，溶质的对流通量记为 。单位体积土壤水溶液中所含有的溶质质量，称为溶质的浓度，记为 。 cJccJqccJv c现在学习的是第61页，共78页2.6.1溶质运移的对流和水动力弥散2.6.3.2溶质的分子扩散溶质的分子扩散是由于分子的不规则热运动引起的，其趋势是溶质由浓度高处向浓度低处运移，以求最后浓度的均匀。自由水中溶质的分子扩散通量符合Fick第一定律，即00dcJDz 在土壤中，溶质的分子扩散规律同样符合Fick第一定律,dscJDz 现在学习的是第62页，共78页若将土壤孔隙设想为均匀的圆形毛管，半径为若将土壤孔隙设想为均匀的圆形毛管，半径为 ，管轴，管轴线与土水势梯度方向一致，此时，管内半径为的任一线与土水势梯度方向一致，此时，管内半径为的任一点的流速点的流速 可表示为可表示为2.6.1溶质运移的对流和水动力弥散2.6.3.3溶质的机械弥散及水动力弥散aR v r 2214av rRrL28aRvL以上两式即为Poiseuille方程 现在学习的是第63页，共78页n由于土壤颗粒和孔隙在微观尺度上的不均匀性，由于土壤颗粒和孔隙在微观尺度上的不均匀性，溶液在流动过程中，溶质不断被分细后进入更为溶液在流动过程中，溶质不断被分细后进入更为纤细的通道，纤细的通道，每个细孔中流速的方向和大小都不一每个细孔中流速的方向和大小都不一样样，正是这种原因使溶质在流动过程中逐渐分，正是这种原因使溶质在流动过程中逐渐分散并占有越来越大的渗流区域范围。溶质的这散并占有越来越大的渗流区域范围。溶质的这种运移现象称为种运移现象称为机械弥散机械弥散。 2.6.1溶质运移的对流和水动力弥散2.6.3.3溶质的机械弥散及水动力弥散现在学习的是第64页，共78页由机械弥散引起的溶质通量可写成类似的表达式由机械弥散引起的溶质通量可写成类似的表达式: hhcJDvz 分子扩散和机械弥散的机理是不同的，但上式与 的表达相似，而且一般都同时存在，实际上难以区分，因此，将分子扩散与机械弥散结合，称为水动力弥散。水动力弥散所引起的溶质通量可表示为: ,DhsshccJDvDDvzz 2.6.1溶质运移的对流和水动力弥散2.6.3.3溶质的机械弥散及水动力弥散dscJDz 现在学习的是第65页，共78页2.6. 4溶质运移的基本方程 2.6 2.6 土壤中溶质运移的基本方程土壤中溶质运移的基本方程溶质运移的对流和水动力弥散作用，决定了溶质的总通量 为对流通量 和水动力弥散通量 之和。 JcJDdhJJJ即,shcJDvqcz cJtz 由和,shcqccDvtzzz,shcccDvqtzzz得一维的溶质运移基本方程 现在学习的是第66页，共78页 当太阳辐射能源源源不断地到达地表，除一部分加热当太阳辐射能源源源不断地到达地表，除一部分加热近地面空气外，大部分均被土壤所吸收。随着表土温度的近地面空气外，大部分均被土壤所吸收。随着表土温度的提高，热量逐渐流入土壤深层，称为提高，热量逐渐流入土壤深层，称为正值交换正值交换。在冬季或。在冬季或夜间，很少辐射能到达地面时，突土壤深层储存的热量流夜间，很少辐射能到达地面时，突土壤深层储存的热量流向土表，称为向土表，称为负值交换负值交换。土壤的正值和负值交换统称为。土壤的正值和负值交换统称为土壤土壤的热量交换的热量交换，它是决定土壤温度的基本因素。，它是决定土壤温度的基本因素。2.7 2.7 土壤中热流基本方程土壤中热流基本方程现在学习的是第67页，共78页在 方向进入单元体内的热量： zyxJz T T流出的热量为：在 方向流进和流入的热量差为： zzyxzJzT TzyxzTz22接下页0zTJzzT TzyxzJyxJzzT TT T将Fourie公式代入：2.7 土壤中热流基本方程土壤中热流基本方程2.7.1 土壤热流基本方程现在学习的是第68页，共78页同时考虑 和 方向上的热量变化，单元体总的热量变化可表示为： xyzyxyTxTzTyxz 222222根据热量平衡式，单位时间内单元体热量的变化为： TVCQv vtTzyxCQv v如假设 ，两式相等，即可得到下式：yxz222222yTxTzTtTCv v接下页2.7 土壤中热流基本方程土壤中热流基本方程2.7.1 土壤热流基本方程现在学习的是第69页，共78页重写上式：222222yTxTzTtTCv v如只考虑一维情况，则：22zTtTCv v用土壤热扩散率表示，则有：22zTatT以上三式为非稳定条件下土壤热运动偏微分方程，也就是土壤热流的基本方程。2.7 土壤中热流基本方程土壤中热流基本方程2.7.1 土壤热流基本方程现在学习的是第70页，共78页2.7.2 土壤热流基本方程的解 （1）解析解tATT2sina a年平均温度 波的振幅 右图示为土表10深度多年月平均温度的观测值与正弦模型曲线。土壤表层温度和时间的正弦模型曲线如下式： 时间正弦波周期 2.7 土壤中热流基本方程土壤中热流基本方程现在学习的是第71页，共78页2.7.2 土壤热流基本方程的解土壤热流基本方程的解 （1）解析解定解问题： 22zTatT200tztATtzTsin,a azTtzTa a,lim对定解问题作Laplace变换，得： 220zpzTazTpzTp, 10ftfpttfd dd d由 22pzftzf,2 22 22 2z zz z由Laplace变换2.7 土壤中热流基本方程土壤中热流基本方程现在学习的是第72页，共78页（1）解析解解析解220papTpTa a,pTpzTza a,lim由Laplace变换后的定解问题得： aTpzTapzpzTa a,d,d22上式为一个二阶线性非齐次常微分方程，其通解是：pTapzCapzCpzTa aexpexp,21 5pcc 由zTtzTa a,lim 32121pgcpfctgctfc 422ptsin由00tztATtzTsin,a a2.7 土壤中热流基本方程土壤中热流基本方程现在学习的是第73页，共78页（1）解析解解析解将下边界条件用于通解，有01CpTapzCapzCpzTa aexpexp,21将上边界条件用于通解，有222paC故求得常微分方程的解为： pTpzTza a,lim220papTpTa a,pTapzpapzTa aexp,22逆变换后可求得原方程之解：aadzdztdzATtzT20sinexp,a aLaplace逆变换是一个较为复杂的数学过程，但可以通过查Laplace变换表求得。2.7 土壤中热流基本方程土壤中热流基本方程现在学习的是第74页，共78页（1）解析解解析解dzAexp 土壤中任何深度的温度变化都呈现周期的正弦波形，正弦波的振幅为 。相对于表层的正弦波，深 度为 的正弦波向后推移 个步长， 。 zttzt所有深度的正弦波都有相同的周期 , 。 2所有深度的年平均温度都是 。 a aT求解过程有几点基本限制条件： 只有热传导 ；热扩散率为常数； 土壤表面的温度近似于正弦波型变化。 问题：为什么正弦模型用于年温度变化要比用于日温度变化更接近实际观测值？2.7 土壤中热流基本方程土壤中热流基本方程现在学习的是第75页，共78页2.7.2.22.7.2.2数值解数值解 如上所述，土壤热流基本方程的解析解只能在特定的基本如上所述，土壤热流基本方程的解析解只能在特定的基本条件下求得。对大多数实际情况，像水流基本方程一样，方程条件下求得。对大多数实际情况，像水流基本方程一样，方程可用数值法求解。可用数值法求解。 2.7. 2土壤热流基本方程的解2.7.2.3水、热耦合问题 一个思路是将热和水的问题在较小步长情况下分别单独求解，再回过头来作一定修正，这样反复交替进行，直到求得较满意的结果。 第2章结束现在学习的是第76页，共78页现在学习的是第77页，共78页感谢大家观看现在学习的是第78页，共78页