基于.RETC与MATLAB的土壤水力学参数计算.ppt

基于RETC与MATLAB的土壤水力学参数计算 汇报人：Iris Liu 2013.4.19 主要内容 土壤水力学参数简介 模型精度评价 土壤水分特征曲线的确定方法 RETC软件和MATLAB软件介绍 实例演示 土壤水力学参数简介 土壤水分的动态变化是土壤-植物-大气连续体（SPAC）系统的中心环节和联系纽带 。利用土壤水动力学模型模拟土壤中水分的动态变化有助于进一步认识SPAC系统的水分运移规律和水分迁移的定量关系。而计算土壤水分和水文通量最基本的方法就是非饱和土壤水分运动基本方程Richards方程 ： 取渗透系数K(θ)和比水容量C(θ)的比值，定义为水的扩散系数D(θ)，即: 对于垂向一维流，可变为： 考虑一维垂向的土壤水分运动和植物根系吸水，采用以 含水率 θ 为因变量的 Richards 方程表达式为: 求解非饱和土壤水分运动方程进而预报非饱和土壤水分运动必须首先获得土壤水动力学参数。下面简单介绍这几个参数： 1）土壤水分特征曲线：是土壤含水量和基质势的关系曲线，反映了土壤水的能量和数量关系。 2）非饱和水力传导度：是反映土壤水分在压力水头差作用下流动的性能，一般在饱和土壤中导水率称为渗透系数。 3）土壤扩散率D(θ)是单位含水量梯度下非饱和流的通量，其与土壤含水率θ或土壤基质势的函数关系通过试验测定。 4）比水容量是指单位基质势的变化引起含水量的变化，数值上等于土壤水分特征曲线斜率的倒数 。 本次主要介绍土壤水分特征曲线的求解。 土壤水分特征曲线的确定方法 土壤水分特征曲线的确定方法 直接方法 间接方法 负压计法 砂性漏斗法 压力板法 稳定土壤含水率法 经验公式法 土壤转换函数方法 分形方法 人工神经网络方法 在计算土壤水分特征曲线中经验公式法中应用较多的有以下几种:Brooks-Corey模型、Gardner模型、Van Genuchten模型（简称VG模型），Gardner-Russo模型。各模型表达形式如下： Brooks-Corey模型 式中：Se—饱和度； θ—体积含水量，cm3/cm3； h—压力水头，cm； θs—饱和含水量，cm3/cm3； α—进气压力值的倒数，即α=1/hb； λ—大于零的正常数，反应土壤空隙大小分布。 Van Genuchten模型： 式中：α，n，m—经验拟合参数，m=1-1/n。 其它同上。 Gardner模型： 式中：a、b为大于0的正常数。其它同上。 在所有这些模型中，Van Genuchten模型以其与实测数据曲线拟合程度较好而得到广泛应用。 VG模型参数的求解可以采用RETC软件、MATLAB软件和SPSS软件进行。本次主要介绍利用RETC软件、MATLAB软件进行求解。 RETC软件 RETC软件是由美国国家盐改中心实验室研制开发的模型，它根据1913个不同岩性的颗粒组成、干密度、土壤水分特征曲线、饱和渗透系数等实测数据，应用土壤传递函数模型（PTFs法）和神经网络技术建立了土壤水分特征曲线参数和饱和渗透系数与土壤颗粒组成、干密度之间的函数关系，从而能够从获取的某一方面资料中预测出其它难以求取的数据。 RETC和MATLAB软件介绍 土壤各粒径百分比 实测h—θ值 土壤水分特征曲线 土壤容重（Bulk Density） MATLAB软件 MATLAB软件是由Math Works公司于1984年推出的，它可以实现数值分析、优化、统计、偏微分方程数值解等领域的计算和图形显示功能的软件包。其函数的形式分类成库，使用时可以直接调用。本次分别调用了matlab中的非线性拟合函数lsqcurvefit和非线性回归函数nlinfit函数和nlintool函数求解VG模型中的四个参数θs、θr、α和n，并对其进行了比较分析。 模型精度评价 求参效果的好坏，需要对得到的拟合方程的精度进行检验和评价。通常检验的指标有拟合优度（R2）、均方根误差（RMSE）