Muskingum法及其分段连续演算的若干理论探讨.pdfVIP

第 卷 第 期 水 科 学 进 展 ， *1 3 EF9*1 =E93 年 月 ， )++) ** 2;2=?@ A= B2C?D @A?=? =EGH )++) !#$%’( 法及其分段连续演算的若干理论探讨 芮孝芳 （河海大学水资源环境学院，江苏 南京 )*++,- ） 摘要：根据水力学原理和洪水波运动理论，通过对 法的关键参数 与特征河长、扩散 !#$%’( ! 波动力方程和运动波数值扩散之间关系的分析，给出了 更为全面的物理解释；证明了 ! !#$ %’( 法槽蓄方程是扩散波动力方程近似的表达；指出了!#$%’( 法演算公式在一定条件下是 扩散波方程的二阶精度解。讨论了 法的使用条件和分段连续演算的必要性；应用 变 !#$%’( / 换方法导出了!#$%’( 法的分段连续演算的汇流系数公式。 关 键 词：!#$%’( 法；特征河长；扩散波方程；运动波数值扩散；分段连续演算 中图分类号： 文献标识码： 文章编号： （ ） 0 11- 2 *++* 34,* )++) +3 3-) +4 法自 世纪 年代中期创立以来，已在世界上众多河流的洪水演算中得到了 !#$%’( )+ 1+ 广泛的应用，其特点是计算简单、所需资料少、当满足其使用条件时精度令人满意。但该法创 立伊始，其槽蓄方程的表达式纯属一种经验假设，其中关键参数 ，由于形式上似乎起着流量 ! “权重”的作用，故传统上称为“流量比重因子”。直到 世纪 年代之前，国外许多学者一 )+ 6+ ［］ ［］ 直认为 的值应在 之间 * ) ! + 8 * 或者在 之间 。中国学者在 世纪 年代通过大 + 8 +96 )+ 6+ 8 3+ 量实践发现，欲使 法的计算精度令人满意， 的取值没有发现大于 的，并可取 !#$%’( ! +96 负值，即 。而且发现， 值的大小反映了洪水波在传播过程中坦化变形的程度，与河 ! +96 ! ! 段长度、河底比降、干支流顶托等因素有关（表 ）。这不是偶然的，文献［ ］曾对此作过一定的 * 1 研究。本文将在前人研究的基础上，根据洪水波运动理论，对!#$%’( 法的理论基础再作进 一步探讨，给出 的更全