马斯京根法洪水演算公式参数确定方法的商榷

1、 文章编号:100724929(2009 0820057202马斯京根法洪水演算公式参数确定方法的商榷常蒲婷1, 杨侃1, 赵建刚2, 沈雪娇1(1. 河海大学水文水资源与水利工程科学国家重点实验室, 江苏南京210098;2. 华北水利水电学院水利学院, 河南郑州450011摘要:针对近年来最小二乘法用于马斯京根法出现的一种趋势, , 剖析了马斯京根法进行洪水演算的原理, 指出了马斯京根法中参数的物理意义, 从理论上论证了康吉演算法推算具有衰减特性的洪水波运动必须满足的条件, 时造成的参数物理意义的模糊, , 所得结论愿与读者商榷。关键词:流量演算法; 马斯京根法; 最小二乘法; 中图分类号

2、:TV11文献标识码Method for Parameters of Method for Flood C alculation FormulaPu 2ting 1,YANG K an 1, ZHAO Jian 2gang 2,SHEN Xue 2jiao 1(1. State Key Laboratory of Hydrology 2Water Resources and Hydraulic Engineering , Hohai University ,Nanjing 210098,China ; 2. Hydraulic Engineering Department ,North Chi

3、na University of Water Conservancy and Electric Power ,Zhengzhou 450011,China Abstract :Aiming at the trend of the least square method used in Muskingum method in recent years , starting f rom the hydraulic theo 2ry of the flow algorithms method , this paper analyzes the flooding calculations princi

4、ple of the Muskingum method and points out the physical meaning of parameters of Muskingum method. Then through theoretical explanation of the Muskingum method , this paper theoretically demonstrates the conditions which should be met when the Cunge algorithms method is used for flood with attenuati

5、on characteristics. It also points out the f uzzy of the parameters physical meaning and possible arising problems when the least square method is used in the Muskingum method , and the method is not able to guarantee the establishment of the conditions. K ey w ords :flow algorithms method ; Musking

6、um method ; least square method ; determination of parameters收稿日期:2009201229作者简介:常蒲婷(19842 , 女, 硕士研究生, 主要从事水资源规划与管理方面的研究工作。1马斯京根法简介马斯京根(Muskingum 法是McCarthy 于1934年提出, 并在美国马斯京根河上首先应用的一种用于河道洪水演算的方法。由于该法在数学上比较简单, 在一般河道的洪水演算中效果较好, 因而到现在仍在世界范围内被广泛应用。马斯京根法属于流量演算法的一种, 流量演算法是在圣唯南(Saint -Venant 方程组简化的基础上, 根据

7、河段的水量平衡原理与槽蓄关系把河段上游断面的入流量过程演算为下游断面的出流量过程的方法, 也就是通过河槽调蓄作用的计算反映河道洪水波运动的变化规律, 它是河道非恒定流计算中一种水文学的近似简解方法1。一般情况下, 槽蓄关系(W =f (O 因不同洪水波受附加比降的影响而呈现顺时针绳套、逆时针绳套和单值关系, 如果槽蓄关系呈单值线性函数形式, 则流量演算可大为简化。这就是流量演算法的一般思路。具体到马斯京根法来说, 它是以连续方程表示的水量平衡方程和以槽蓄方程代替的动量方程来描述的, 即:d t=I -O (1 W =k (xI +(1-x O (2式中:I 和O 分别为研究河段的上段面入流和下

8、段面出流; W 为该河段内的蓄水量; k 和x 则为参数, k 值是槽蓄曲线的坡度, 从物理意义上讲, k 值等于相应蓄水量W 下稳定流状态河段传播时间; x 为流量比重因子; x I +(1-x O 为示储流量, 用Q 表示。75节水灌溉2009年第8期 (1 、(2 两式的差分解为:O 2=C 0I 2+C 1I 1+C 2O 1(3其中:C 0=k -kx +0. 5t C 1=k -kx +0. 5t C 2=k -kx +0. 5t C 0+C 1+C 2=1. 0(4上述就是大家所熟悉的马斯京根法河道洪水演算公式。对于给定的河段, 演算时段长t 确定后, 待定参数有k 和x 。马斯

9、京根法假定x 、k 是常数, 通过反复调整x , 使得河段蓄量W 与示储流量Q 成单一的直线, k 就是这个曲线的斜率:k =d W /d Q 。通常步骤为:假设不同的x 计算示储流量Q , 直至W 和Q 呈单值关系或近似为单值关系为止。但实际上x 并不是常数, W -Q 不是直线, 而是曲线, 所以也是变化的。这需要从k 、x 的物理意义来解释。现已证明:k 值等于在相应蓄量下恒定流状态下的传播时间, 这是k 的物理概念。显然k 随恒定流流量而变化, 取k 为常数是有误差的; 而x 表示本演算河段的河槽调蓄能力, 显然对于不同数量级的洪水x 也是不同的。文献3、4, , C 01+C 2=1

10、. 0为约束条件, C 0、C 1、C 2, , C 0、C 1、C 2反算得到, 这就是文献的一般思路。由以上论述可见:文献3、4中是将马斯京根法作为黑箱模型来处理, 这种先优化C 0、C 1、C 2而后反算k 、x 的做法不但不能保证示储流量Q 和蓄水量W 呈单值关系, 而且使得k 和x 本已揭示清楚的物理意义丧失, 使得本已具有明确物理意义的马斯京根法重新沦为“黑箱模型”, 这种计算方法是值得商榷的。2马斯京根法的理论解释早在1969年, 法国水力学家康吉(Cunge 在研究运动波方程时发现了数值扩散现象, 他采用四点线性隐式差分格式, 对运动波方程进行差分, 取二阶近似就是扩散波方程2

11、。四点线性隐式差分为(见图1 :9tn+1n n+1nt (5 9xn+1n+1n n x (6式中:Q 的上标为时间, 下标为断面位置, 例如Q n +1j 表示断面j在n +1时刻的流量; X 和Y 分别为时间和空间的差分格式的权重系数; t 为时间步长; x 为空间步长。将式(4 、(5 代入运动波方程, 并取Y =0. 5, 令K =x/C k (C k 为波速 , 并解出Qn +1j +1, 得:Qn+1j+1=C 0Q n j+C 1Q n+1j+C 2Qn j+1(7其中:C 0=K -KX +0. 5t C 1=K -KX -0. 5t C 2=K -KX +0. 5tC 0+

12、C 1+C 2= 1图1偏心差分格式公式(7 其实就是马斯京根法的推流公式。而欲使该式能推算具有衰减特性的洪水波运动, :X =2k (8; k ; x 为河段长; X 为时间。与洪水波扩散系数, 波速, 河段长等因素有为正值, 所以X 总是满足X 0. 5。通过上述分析可以看到:马斯京根法已经由经验试错上升到了具有一定水力学理论支撑的理论水平。至于k 、x 两个参数更是已经具有了明确的物理意义, 而且只有在X 0. 5的条件下, 采用马斯京根法进行洪水演算才能得到稳定解。文献3、4中, 采用公式(4 反向推导得到了k 、x 的计算公式, 若要使得X 0. 5成立的必要条件是:由文献3、4中采

13、用的最小二乘优化算法优化得到的C 0、C 1、C 2和采用传统的试错法得到的k 、x 进而计算得到的C 0、C 1、C 2分别相等或者误差在可接受范围内。而这一点由于河道汇流的非线性现象的存在并不是总能满足的, 因此, 文献3、4采用的方法并不能从理论上保证X 0. 5这个条件的成立, 而只能通过数值计算偶然满足, 关于此点本文作者愿与文献3、4的作者商榷。3结语本文作者通过流量演算法的水力学理论出发, 阐述了采用流量演算法进行河段洪水演算的基本原理, 进而通过马斯京根法的理论解释-康吉演算法, 从两个方面论述文献3、4解决问题方法中存在的不足之处, 具体是:不但使得k 、x 两个参数的物理意义丧失, 而且由C 0、C 1、C 23个参数反算得到的的表达式并不能从理论上保证满足X 0. 5这个条件, 因此应用最小二乘法于马斯京根法公式参数的确定, 这样的计算方法还称为马斯京根法是不确切的, 称之为黑箱模型更为贴切; 加上水文事件的随机性, 甚至可能导致错误的计算结果,