Richards.docVIP

Richards方程是由谁提出的？是用来干嘛的？ Richards方程最先是由Richards这个人在1931年研究流体通过多孔介质中毛细管传导作用时推导出来的[1]。 对于各向同性的土壤、不可压缩的液体、三维情形的非饱和水流运动的控制方程即Richards方程： 式中：为含水量， 为时间，为渗透系数，为非饱和土壤的总土水势，，， 表示坐标轴方向。 非饱和土壤水分运动方程几种形式 描述非饱和土中水的流动控制方程一般采用 Richards微分方程。对于一维非饱和土壤水分运动问题，Richards方程常常表示为3种标准的形式:基于压力水头，基于含水率和基于位置坐标形式。 以负压水头为因变量的基本方程[2]： 式中：为水的密度，为重力加速度，为土壤含水量与基质吸力关系曲线的斜率，为负压水头（基质势），为时间，为渗透系数，为土壤的深度。 以含水率为因变量的基本方程[2]： 式中：为含水量，为时间，为扩散率，为土壤的深度，为渗透系数。其中， 为基质势。 以位置坐标为因变量的基本方程[2]： 当=0时，可以得到如下的水分运动方程： 式中，符号意义同上。该方程主要用于解析或半解析方法对非饱和土流动的求解。 我所研究非饱和分层土壤水分运动准备采用的方程 式中： ，为第层土壤的压力水头（基质势）（cm），土壤的水力传导率（cmh-1）,为时间（h），为土壤的深度（cm）。 求解方程（5） 第一步先推求土水特征曲线 推求土水特征曲线的方法有很多，归纳起来可分为直接测定方法和间接推求法。 直接测定方法。即为在实验室中通过测定土壤一系列含水量与其对应的基质势的测定方法统称。主要可以分为张力计法[3]、压力膜法、离心机法、砂芯漏斗法、平衡水汽压法。 间接推求法是指利用粒径分布、容重、有机质含量等土壤基本性质，通过某种经验的或半经验的关系建立起来的预测水分特征曲线的一大类方法。按照其基本原理和构建方法，间接方法大致可以分为3类，即：土壤转换函数方法[4](PTF，Pedotransfer functions)、物理—经验方法(Physico empirical method)、分形几何方法(Fractal method)。 土壤水分特征曲线是高度非线性的函数，受到多种环境和土壤因素的影响，准确测定这一关系非常困难。而且绝大多数方法成本较高而且十分繁杂和费时，同时实验结果常存在明显的不确定性，土壤在水平向和垂向的空间变异性也限制了直接测定方法在实际工作中的应用。因此本文准备采用间接方法来推求土水特征曲线。本文准备采用常用的Brooks-Corey 土水特征曲线。 Brooks-Corey 模型[5-6]中的含水量与压力水头的函数关系式如下： 式中：为体积含水率；h为压力水头[L](为基质吸力)；为饱和含水率；为残留含水率；为进气压力值[L]；为孔隙分布的指数参数。 将上述关系式（6）代入Burdine[6]提出的关系式 根据上述两个公式可以推导出渗透系数的表达式 式中：为渗透系数[LT-1]； 为饱和渗透系数[LT-1]；其余符号意义同前。 第二步率定Brooks-Corey 模型中的参数（、、、） 先通过实验测量得到土壤的含水量与压力水头（基质势）一些实验数据，然后再通过粒子群算法优化Brooks-Corey模型中的参数，最终得到Brooks-Corey模型土水特征曲线。 第三步将Brooks-Corey模型中渗透系数表达式（8）代入公式（4）得到只以基质势为因变量的水分运动方程。 第四步用有限差分法解方程（9） 用有限差分法解此方程之前，先要将研究区域z划分成有限个小区，还得将时间t分段。 然后再把方程化成以为因变量的差分形式： 式中：为的函数，为的函数，，为的函数，，表示空间网格顺序号，表示时间网格顺序号。 第五步确定初始条件和边界条件 初始条件： 边界条件： 上边界条件 下边界条件 饱和-非饱和土壤中的溶质迁移转化方程几种形式 描述非饱和土壤一维溶质运移的基本方程[7]： 式中， 为土壤洁液浓度；其溶质可为盐分、化肥、农药中的有害物质、重金属等；为土壤容积含水量；为水动力弥散系数； 为水分通量；为时间坐标( min) 。 描述污染物质在饱和及非饱和土壤中迁移转化的基本方程[8]： 式