《冬小麦根层土壤水量平衡的系统动力学模型 Soil water balance model in root zone of winter wheat based on system dynamics》.pdfVIP

第26卷第10期 农业工程学报 、，01．26No．10 2010焦 oftheCSAE Oct．20lO21 10月 Transactions 冬小麦根层土壤水量平衡的系统动力学模型 陈皓锐，黄介生※，伍靖伟，杨金忠 武汉大学水资源与水电工程科学国家重点实验室，武汉430072 摘要：为获取冬小麦根系层水量转化情况，该文采用系统动力学的建模思想和Vensim软件构建了冬小麦一维逐同土壤 水量平衡模型。模型将2m土层概化为十个串联的水箱，计算了灌溉降雨后的土壤水分下渗、土壤蒸发、作物蒸腾、毛 管上升补给和水分重分配等物理过程。利用河北省石津灌区军齐干渠北二支一斗渠2007--2009年两季冬小麦的田间试验 资料对模型进行了率定和验证，结果显示率定期和验证期的平均残差比例和分散均方根比例均在15％以内。三种极端条 件测试和六种参数的敏感性测试以及与Hydrus．1模型的比较表明模型假定合理，没有发牛结构性错误。对灌区两季冬小 麦生育期的土壤水分转化进行模拟，结果表明降雨和灌溉是主要供水水源，毛管水上升量很小，底部渗漏较大，而土壤 储水量变化很小。 关键词：土壤，水，系统动力学，根层土壤水量平衡，冬小麦 doi：10．3969巧．issn．1002—6819．2010．10．004 中图分类号：$271 文献标志码：A 文章编号：1002—6819 2010 一10-0021-08 陈皓锐，黄介生，伍靖伟，等．冬小麦根层土壤水量平衡的系统动力学模型[J】．农业工程学报，2010，26 10 ：21--28． Chen a1．Soilwaterbalancemodelinrootzoneofwinterwheatbasedon Jingwei，et system Haomi，HuangJiesheng，Wu with ofthe Chinese dynamics[J]．TransactionsCSAE，2010，26 10 ：21--28． inEnglishabstract 视化模型作用结构的编程语言相比，系统动力学方法又 0引言 可以提供诸多的政策分析工具使其功能上更为强大。而 且随着研究的深入，当模拟系统与真实水循环系统的复 系统动力学 systemdynamics，SD 是麻省理1 学 院经济学家J．W．Forrester教授于1950年提出的一项模拟 杂性越接近时，系统动力学方法基于“反馈+流量图”的 决策技术【11。该方法综合了系统论、信息论、控制论、决 可视化建模思想和自带的政策分析工具带来的高效优势 策论和仿真等多类研究成果，其分析对象主要针对复杂 会越来越得以凸显。 非线性的社会经济宏观大系统，伴随着WORLDII模型 本文按照系统动力学建模思想构建了一维的逐日冬 及以此为基础的《世界动力学》【2J，WORLDIII模型及以 小麦根系层土壤水量平衡模型。采用河北省石津灌区军 此为基础的《增长的极限》【3】和《趋向全球的平衡》14]的 齐干渠北二支一斗渠2a的田间实测资料对模型进行了率 发表以及美国国家模型的应用【5j，系统动力学方法于1970 定和验证，利用Vensim软件自带工具验证了模型在极端 —1980年中期逐渐成熟起来，并开始应用于能源、交通、 条件和参数变动下的行为有效性及田间水分转化。 生态、环境等多个领域。水资源领域也在20世纪90年 1 模型描述 代中后期引入该方法并在水资源开发决策、水资源承载 能力分析、水资源可持续管理等宏观研究上涌现了一大 把冬小麦根系发育所在的0～2m土层视作独立系 批研究成果睁引，但在水循环尤其是土壤水循环的模拟仿 统，按照20 真方面使用该方法尚不多见，本文即采用系统动力学思 设置为天。第l层土壤入流有灌溉水量、降雨量、下层 想对冬小麦根系层土壤水量平衡进行模拟以丰富本领域 土壤重分配入流的水量，出流有下渗垒第2层的水量、 的研究手段。 土壤蒸发和根系吸水；第2～9层土壤入流有上层土壤的 与常用的编程语言如Fortran．、Matlab等相比，使用 下渗水量、下层土壤重分配入流的水量，出流有根系吸 系统动力学方法建模的优点在于其结合了图论理论，使 水、渗入下层的水量和重分配至上层土壤的水量；10层 得模型结构透明，而且该方法还可以提供大量的控件方 土壤入流有上层土壤的下渗水量、地下水毛管上升补给 便流量图的绘制，加快了建模的速度；而与VB等能够可 水量，出流有流出2m下边界的深层渗漏、重分配至第9 层的水量。基于上述分析并结合水量平衡原理，各层土 收稿日期：2009．12-02修订R期：2010-09-05 壤逐日土壤水量平衡方程如下。 基金项目：国家重点基础研究发展计划 973 项目 2006CIM03406 资助 第l层：