一维畦灌地表水流-土壤水动力学耦合模型Ⅰ建模 Coupled surface-subsurface hydrodynamic model for border irrigationⅠ. Model development.pdfVIP

水 利 学 报 2013年5月 SHUILI XUEBAO 第44卷 第5期 文章编号：0559-9350（2013）05-0570-08 一维畦灌地表水流-土壤水动力学耦合模型 Ⅰ：建模 1，2 1，2 1，2 1，2 董勤各 ，许 迪 ，章少辉 ，白美健 （1. 中国水利水电科学研究院 水利研究所，北京 100048；2. 国家节水灌溉北京工程技术研究中心，北京 100048） 摘要：畦灌地表水流-土壤水动力学耦合模型可从物理机制上综合模拟地面灌溉过程中地表水与土壤水系统之间 的互动关系，为开展区域尺度地面灌溉系统性能评价提供了基础工具。利用二阶时空离散精度的混合数值解法求 解一维全水动力学模型，建立了一维畦灌地表水流运动模型，并利用二阶时间离散精度有限差分法和构建的四阶 空间离散精度有限体积法解算一维土壤水动力学模型。同时，以地表水和土壤水质量守恒方程交替迭代计算的地 表水深收敛值为基础，利用动量守恒方程获得地表水流速场，实现了地表水-土壤水系统间的有效迭代耦合，构 建起基于数值解法和迭代耦合模式同步改进的一维畦灌地表水流-土壤水动力学耦合模型。 关键词：畦灌；迭代耦合；数值解法；地表水深；质量守恒；动量守恒 中图分类号：S275.3；S147.3 文献标识码：A 1 研究背景 利用畦灌地表水流-土壤水动力学耦合模型可从物理机制上建立起地表水与土壤水系统之间的互 动关系，有效模拟灌溉前后饱和-非饱和土壤水动态变化过程，为开展区域尺度地面灌溉系统性能评 价提供了基础工具。畦灌地表水流-土壤水动力学耦合模型通常包括外耦合、迭代耦合、完全耦合等 3种模式［1-3］。外耦合模式是以地表水流模型模拟获得的地面水深作为上边界条件，直接赋给土壤水 动力学模型或简化的入渗模型，两者间无迭代过程［4-6］。外耦合模式虽然计算效率较高，但因未考虑 土壤水对地表水的反馈作用，致使模拟精度较低。迭代耦合模式是在任一时间步内，以入渗水深或 地表水深作为计算收敛的判别指标，交替迭代地表水流和土壤水运动控制方程，直到计算结果符合 预设的模拟精度时为止。迭代耦合模式的计算效率相对较低，但模拟精度得到改善，目前为人们普 遍采用［7-9］。但在地表畦灌下，迭代耦合模式的计算区域在水平向与垂向的空间尺度间约存在两个量 级的差异，这降低了该模式的模拟精度［8-9］。若采用混合精度的空间离散格式，可在一定程度上改善 ［10］ ［9，11］ 因量级差异引起的较差模拟精度问题 。另外，从地表水与土壤水运动控制方程可知 ，前者存在 着质量和动量守恒，后者仅有质量守恒，但现有迭代耦合模式却未详细区分该差异，这或许也是导 致现有迭代耦合模式具有较低计算效率的原因之一。完全耦合模式是同步解算地表-土壤水动力学模 型［1-3，12-13］，其虽具有震荡敏感性低和质量守恒误差小等优点，但因系数矩阵规模较大，具体求解较为 困难［1-3］，故至今鲜见其用于实际灌溉模拟。 本文针对现有畦灌地表水流-土壤水动力学迭代耦合模式存在的水平向与垂向空间尺度量级差异 降低了模拟精度、未详细区分地表水-土壤水运动质量和动量守恒过程可能降低计算效率的主要问 题，在采用二阶时空离散精度混合数值解法求解畦灌地表水流运动模型基础上，利用二阶时间离散 精度有限差分法和构建的四阶空间离散精度有限体积法求解土壤水动力学模型，通过地表水-土壤水 收稿日期：2012-04-06 基金项目：国家863计划重点项目课题（2011AA100505）；国家自然科学基金资助项目51209277） 作者简介：董勤各（1982-），男，河南南阳人，博士生，主要从事灌水管理方法与技术、水土资源高效利用研究。