土壤持水曲线van Genuchten模型求参的一种新方法.pdfVIP

维普资讯 第44卷 第 6期 土 壤 学 报 Vo1．44．No．6 2007年 11月 ACTA PEDOLOGICA SINICA NOV．，2007 土壤持水 曲线vanGenuchten模型求参的一种新方法* 刘贤赵 李嘉竹 张振华 (鲁东大学地理与规划学院，山东烟台 264025) A NEW METHOD TO ESTⅡ ATING PARAM ETERSOF VAN GENUCHTEN M ODEL FOR SOⅡ W ATER RETENTION LiuXianzhao LiJiazhu ZhangZhenhua (CollegeofGeographyandPlanning，LudongUniversity，Yantai，Shandong 264025，China) 关键词 土壤持水 曲线 ；vanGenuchten模型；DPS数据处理系统；Marquardt方法 中图分类号 S152．7 文献标识码 A 土壤持水曲线是研究土壤水动力学性质必不可 tem，简称 DPS)数据处理系统求解此问题 ，则使该 问 少的重要参数。在 已经建立 的众多数学模型中，van 题变得非常容易。 Genuchten模型以其与实测数据拟合程度好而得到 广泛的应用 ，而运用该模型的关键是如何获得其 中 1 vanGenuchten模型 的4个参数。仅就我 国而言，不少学者投入大量 的 精力发展 了确定 vanGenuchten模型参数 的方法。 vanGenuchten模型是表征非饱和土壤持水能力 邵明安等u J基于一维土壤水分运动的Richards方 的重要曲线，是获得其他土壤水动力参数及土壤水 程提 出了推求土壤 vanGenuchten模 型和 Brooks． 分参数的基础。它的具体表达形式为[m] Corey模型参数的简单人渗法 ；王金生等 将最小二 丹 一 丹 。 二 乘法和非线性单纯形法相结合拟合了vanGenuchten … 。[1+(d ·h)] 模型参数 ；徐绍辉等 J也借助最小二乘法并结合 Pi． 式 中，0为土壤体积含水量 (cm3emI3)，h为负压水 card迭代法拟合 了砂质黏壤土的vanGenuchten模型 头(一amH20)， 、0，分别为土壤饱和含水量和残余 参数；李春友等 J也利用单纯形调优法拟合 van 含水量 (cm3amI3)，a为与土壤平均孔隙半径有关 Genuchten模型的参数 ；魏义长等 J运用 Matlab编程 的参数 (am )，n为曲线形状参数 ，m=1—1／n。鉴 软件对辽西淋溶褐土 vanGenuchten模型的参数进 于土壤水分测定方法的习惯表示 ，本文以土壤水吸 行了推导估算。尽管上述方法均得到了较好的拟合 力值 (+)替代压力水头 (一)，以质量含水量 (gg ) 结果，但这些方法要么借助于土柱入渗试验 ，要么其 替代体积含水量 (cm3amI3)，应用 DPS数据处理系 算法需要编程，或者需要专用软件 ，存在耗时、费力 统进行模型参数求解 。 和利用效率低等问题。 Matlab软件l8J虽然功能强大，但其工作界面对 2 vanGenuchten模型参数求解 我国科学工作者来说较难适应 ，而且要求非常专业 的数学知识和较高的外语水平，算法需要编程 ，这在 DPS数据处理系统是 由唐启义于 1997年设计 一 定程度上 限制 了