基于田块尺度含水率观测的土壤水力参数多模型反演.pdfVIP

第 31卷 第 6期 农 业 工 程 学 报 Vbl|31 NO．6 2015仨 3月 Transactionsofthe Chinese Society ofA culturalEngineering Mar．2015 135 基于田块尺度含水率观测的土壤水力参数多模型反演 刘 昭，杨文元，查元源，杨金忠※ (武汉大学水资源与水电工程科学国家重点实验室， 武汉 430072) 摘 要：该文利用反演工具UCODE与4种Richards方程数值模型 (Ross模型、Picrad一0模型、Picrad—mix模型和Picrad-h 模型)进行耦合，构建了4种不同的反演模型。基于田块尺度含水率观测数据，分别用4种反演模型优化了研究田块的 土壤参数。研究结果表明，4种模型的反演精度依次为：Ross模型、Picard-0模型、Picard．mix模型和 Picard．h模型，但 差异并不显著，反演效率以Ross模型最优。随着网格的加密，各种模型所反演参数的模拟精度改善不明显。本文还讨论 了土壤水运动中 “异参同效性”现象，并提出 “参数 曲线带”的概念——即由反求的同效参数土壤水分曲线和水力传导 度曲线形成的包络图。随着模拟精度要求越高，同效参数越少，“参数曲线带”越窄，并认为反求的同效参数曲线在含水 率观测信息较多的地方交汇。 关键词：土壤；水分；模型；Richards方程；数值反演；异参同效；参数 曲线带 doi：10．3969／j．issn．1002—6819．2015．06．019 中图分类号：S152．7 文献标志码：A 文章编号：1002—6819(2015)一06—135—10 刘 昭，杨文元，查元源，等．基于田块尺度含水率观测的土壤水力参数多模型反演 [J]．农业工程学报，2015，31(6)： 135— 144． Liu Zhao，YangWenyuan，ZhaYuanyuan，eta1．Multipleinversemodelsforestimating soilhydraulicparametersbasedon field—scalemoisturecontentobservation[J]．TransactionsoftheChineseSocietyofAgriculturalEngineering(Transactionsofthe CSAE)，2015，31(6)：135—144．(inChinesewithEnglishabstract) 压和含水率)的观测，在合理范围内调控参数，使得模型 0 引 言 模拟结果与观测结果达到最佳匹配效果6【】。这种方法又称 Richards方程是一种基于达西定律的土壤水动力方 为模型的校正，是获取参数的重要方法。 程，具有坚实的物理基础，因此被广泛用于模拟包气带 由于 Richards方程具有高度非线性，大多解析解仅 土壤水分运动。方程 中负压和水力传导度是土壤含水率 在特定的边界和简化后的参数条件下才能推导，因此其 的非线性函数，土壤水力参数就是描述这种非线性关系 数值解法更具有普适性。在数值方法中，根据 Richards 模型中的经验或物理常数。常用的土壤水负压和水力传 方程的主变量 (含水量、负压等)的选择、数值离散方 导度与土壤含水率的关系式有：Brooks．Corey模型lJJ， 法 (有限单元法、有限差分法)、离散格式质量守恒、 Gardner．Russo模型 【0】和 vanGenuchten—Mualem[】模型 时间步长的处理 (迭代方法、非迭代方法等)等因素， 等。例如 Brooks．Corey模型中包含了五个独立的参数 。 形成了不同的 Richards方程的数值求解模型。例如，应 但是，如何获得这些参数是一个既困难又很重要的问题。 用广泛的土壤水盐运动模拟 的软件——_HYDRU