利用MATLAB拟合vanGenuchten方程参数研究.pdfVIP

土 壤 (Soils), 2010, 42 (2): 268~274 利用 MATLAB 拟合 van Genuchten 方程参数的研究① * 杨改强， 霍丽娟， 杨国义， 李一菲， 钱天伟 （太原科技大学环境科学研究所，太原 030024 ） 摘 要： 选用 MATLAB 的 lsqcurvefit、nlinfit 、lsqnonlin、fminsearch、fminunc、fgoalattain、curvefit、nlintool 等 8 个函 数，分别拟合了 van Genuchten 方程参数，结果表明 lsqcurvefit、nlinfit 、curvefit 和 nlintool 等 4 个函数拟合效果较好，具有源代 码短、收敛性好等特点，特别是在实测数据较少的情况下，MATLAB 的拟合结果甚至优于RETC 和 REDT 。 关键词： van Genuchten ；土壤持水曲线；MATLAB ；拟合 中图分类号： S152.7 在土壤物理学、地下水水文学及土壤环境学中， 单形调优法作为优化手段, 编制了相关的计算程序 土壤的水力传导率是最基本的水动力学参数，而土壤 REDT ，对van Genuchten 方程的 4 个独立的未知参数 持水曲线（SWRC ）又是获得土壤的水力传导率常用 进行了优化拟合, 并与van Genuchten 等开发的 RETC 的手段[1] ，因此土壤持水曲线参数的确定有着非常重 程序的计算结果进行了对比，该软件对初值不需要太 要的意义。 高要求，但是需要一定数量的实测数据点，在实测数 迄今人们已提出了许多经验公式来描述土壤持水 据点较少时，程序运行也会出现不合理结果。 曲线，比较常用的有：Brooks-Corey 模型[2] 、Gardner 总之，目前国内外对于土壤持水曲线参数拟合已 - [3 4] [5] 模型 、van Genuchten 模型 和 Gardner-Russo 模 有较多的研究，但鉴于程序源代码或软件的不公开性 [6] [7] 型 等。徐绍辉等 对此 4 个模型的适应性进行了分 和程序本身的一些弱点，迄今为止还没有一种普遍适 析认为 van Genuchten 模型无论是对粗质地土壤，还 用的应用软件，使得很多土壤专业工作者不得不自己 是较黏质地的土壤，其拟合效果均较好；夏卫生等[8] 重新编制程序，花费了大量精力。而 MATLAB 软件是 通过对国内外土壤水动力学参数研究结果进行分析也 科技工作者广泛使用的一种软件，已经包括了大量的 得出，该模型不仅拟合效果较好，并能和土壤的机械 拟合函数，这就可能使我们的参数拟合的工作大大简 组成和土壤体积质量等联系起来，从土壤本身特性上 化，不需要再去考虑参数拟合的基本方法，花费大量 找到其含义。 时间编制拟合程序。因此，有必要对MATLAB 中的多 目前，对于土壤持水曲线的拟合最常见的为 van 种函数拟合土壤持水曲线参数的效果进行比较，从而 Genuchten 等[9]开发的 RETC 程序，该软件采用了最