第五章 河道洪水演算及实时洪水预报.docVIP

第五章 河道洪水演算及实时洪水预报 河道洪水演算，是以河槽洪水波运动理论为基础，由河段上游断面的水位、流量过程预报下游断面的水位、流量过程。本文着重介绍马斯京根洪水演算方法以及简化的水力学方法。 5.1 马斯京根演算法 马斯京根演算法是美国麦卡锡(G. T. McCarthy)于1938年在美国马斯京根河上使用的流量演算方法。经过几十年的应用和发展，已形成了许多不同的应用形式。下面介绍主要的演算形式。 该法将河段水流圣维南方程组中的连续方程简化为水量平衡方程，把动力方程简化为马斯京根法的河槽蓄泄方程，对简化的方程组联解，得到演算方程。 5.1.1 基本原理 该法的基本原理，就是根据入流和起始条件，通过逐时段求解河段的水量平衡方程和槽泄方程，计算出流过程。 在无区间入流情况下，河段某一时段的水量平衡方程为 (5-1) 式中：、分别为时段初、末的河段入流量；、分别为时段初、末的河段出流量；、分别为时段初、末的河段蓄量。 河段蓄水量与泄流量关系的蓄泄方程，一般可概括为 (5-2) 式中：为河段任一流量对应的槽蓄量。 根据建立蓄泄方程的方法不同，流量演算法可分为马斯京根法、特征河长发等。马斯京根法就是按照马斯京根蓄泄方程建立的流量演算方法。 5.1.2 马斯京根流量演算方程 马斯京根蓄泄方程可写为 (5-3) 式中：为蓄量参数，也是稳定流情况下的河段传播时间；称为流量比重因子；为示储流量。 联立求解式(5-2)和(5-3)，得到马斯京根流量演算公式为 (5-4) 其中： (5-5) (5-6) 式中：、和为马斯京根洪水演算方法的演算系数，，都是、和的函数。对于某一河段而言，只要确定了、和，便可求得、和。于是，由入流过程和初始条件，通过式(5-4)，逐时段演算，就可得到出流过程。 马斯京根演算法的参数、和，可以根据上、下游断面的实测流量过程，用最小二乘法计算出。 从式(5-5)可知，当时，，对是负效应，容易在出流过程线的起涨段出现负流量；但时，，对是负效应，容易在出流过程线的退水段出现负流量。所以要求。 5.1.3 马斯京根连续演算法 (1) 基本公式 为了避免出现负出流等不合理现象，保证上、下断面的流量在计算时段内呈线性变化和在任何时刻流量在时段内沿程呈线性变化，一般要求。1962年赵人俊教授提出了马斯京根分段连续演算法。将演算河段分成个子河段后，每个子河段参数、与未分河段时的参数、的关系为 (5-7) (5-8) 分段连续演算的每段推流公式仍是式(5-4)，但其中的系数采用(5-7)和(5-8)来代替。也可以利用马斯京根汇流系数来进行流量演算。 (2) 计算实例 黄河花园口至夹河滩河段长度为105.4km，采用，分3段进行演算，每段的马斯京根演算参数，。 通过马斯京根连续演算，将花园口1982年7月30日20时～8月4日16日时的流量过程演算到夹河滩。 由已知条件得到：，，，，。 输入数据文件为 3 4.2 0.1 4 1982 7 30 20 8 4 16 30 1140 1650 4780 6240 6320 6190 5850 5580 5400 5580 5520 7020 8710 10400 11500 11400 13000 14500 15300 15200 15100 14300 13100 12800 11700 11100 10600 8670 7930 7350 输出文件为 1440 1444 1532 1965 2945 4174 5143 5666 5823 5770 5662 5619 5789 6359 7370 8616 9801 10818 11836 12929 13892 14519 14735 14543 14041 13391 12681 11942 11131 10175 5.1.4 非线性马斯京根演算法 (1) 基本演算公式 随着马斯京根法的应用，人们发现参数和不是常数，而是随流量变化的，从而人们开始了对非线性马斯京根方法的研究。 水量平衡方程为 (5-9) 槽蓄方程为