一维畦灌地表水流-土壤水动力学耦合模型Ⅱ：验证.pdfVIP

水 利 学 报 2013年6月 SHUILI XUEBAO 第44卷 第6期 文章编号：0559-9350（2013）06-0726-08 一维畦灌地表水流-土壤水动力学耦合模型Ⅱ：验证 1，2 1，2 1，2 1，2 董勤各 ，许 迪 ，章少辉 ，白美健 （1. 中国水利水电科学研究院 水利研究所，北京 100048；2. 国家节水灌溉北京工程技术研究中心，北京 100048） 摘要：根据室内外试验观测数据和模拟结果，对比分析了基于混合数值解法求解的一维畦灌地表水流运动模型以 及基于二阶时间离散精度有限差分法-四阶空间离散精度有限体积法解算的一维土壤水动力学模型的模拟效果， 检验了现有迭代耦合模式改善后的模拟性能。结果表明，利用混合数值解法求解一维畦灌地表水流运动模型具有 比Roe有限体积法更好的数值计算稳定性和收敛性；与一阶时间离散精度有限差分法-二阶空间离散精度有限差 分法及一阶时间离散精度有限差分法-二阶空间离散精度有限体积法相比，采用二阶时间离散精度有限差分法- 四阶空间离散精度有限体积法解算一维土壤水动力学模型具备较好的模拟精度、质量平衡性和收敛速率。与此同 时，改善后的迭代耦合模式呈现出较好的模拟精度和质量守恒性，且计算效率平均提高0.7倍以上。基于数值解 法和迭代耦合模式同步改进后构建的一维畦灌地表水流-土壤水动力学耦合模型，具有提高模拟精度和计算效率 的良好能力。 关键词：畦灌；迭代耦合模型；数值解法；模型验证；模拟效果 中图分类号：S275.3；S147.3 文献标识码：A 1 研究背景 利用畦灌地表水流-土壤水动力学耦合模型可从物理机制上综合模拟地表水与土壤水系统间的互 动关系，有效模拟地面灌溉前后饱和-非饱和土壤水动态变化过程［1-2］，为开展区域尺度地面灌溉系 统性能评价提供了基础工具。畦灌地表水流-土壤水动力学耦合模型通常包括外耦合、迭代耦合、完 全耦合等3种模式［3-4］。其中迭代耦合模式是在相同时间步内，以入渗水深或地面水深作为数值计算 收敛的判别指标，交替迭代求解地表水流运动模型和土壤水动力学模型，直到计算结果符合预设的 模拟精度为止。虽然该模式目前为人们普遍采用，但也存在着有效提高数值模拟精度和计算效率的 问题［5-6］。为此，作者构建起基于数值模拟解法和迭代耦合模式同步改进的一维畦灌地表水流-土壤 水动力学耦合模型［7］，以达到提高数值模拟精度和计算效率的目的。 本文基于室内外试验观测数据，从改进数值解法和迭代耦合模式入手，比较分析基于混合数值 解法［8-9］求解的一维畦灌地表水流运动模型以及基于二阶时间离散精度有限差分法-四阶空间离散精 度有限体积法（FD-FV-2/4）解算的一维土壤水动力学模型的模拟效果，检验经改善后的迭代耦合模 式（ICE-W）的模拟性能，综合评价构建的一维畦灌地表水流-土壤水动力学耦合模型。 2 室内外试验与对比的数值解法及迭代耦合模式 2.1 室内外试验 2.1.1 室内土柱试验 供试土壤分别来自中国水利水电科学研究院北京大兴试验基地和山东东营市 收稿日期：2012-04-06 基金项目：国家863计划重点项目课题（2011AA100505）；国家自然科学基金资助项目51209277） 作者简介：董勤各（1982-），男，河南南阳人，博士生，主要从事灌水管理方法与技术、水土资源高效利用研究。 E-mail：qgdong2011@163.com — 726 — 3 3 麻湾灌区油坊村，经风干后，磨碎过2mm筛备用，体积含水率0.02cm/cm 。试验土柱高80cm、直径 25cm，将土壤分8层（每层高10cm）装入有机玻璃管内，打毛层间土壤表面［10］，干容重ρb 见表1。试 验时土柱上边界条件为定压水头，由马氏瓶供水维持在2cm水深，下边界为自由排水条件，试验时 间120min。待试验结束后，沿土柱每隔10cm采集土壤样本，利用烘干法测定土壤重量含水量，并计 算土壤体积含水率。根据土壤干容重和土壤质地等数据，利用基于人工神经网络方法建立的Rose