第10章流域汇流.pptVIP

第五节 河道流量（洪水）演算 洪水波的形成：设某时刻河道水流处于某种状态，例如稳定流状态。如果由于暴雨径流、水电站运行或闸坝放水等原因，突然有一定水量在此时注入河道，则可以观察到，原来的水面因受到干扰而形成波动，这就是洪水波。 这种由于突然注人一定水量而在河道里增加的流量称为波流量。波流量从河道上游向下游运动，就产生了河道里的洪水波运动。 河道洪水演算方法 河道洪水演算可采用水力学和水文学两类方法。 水力学方法以圣维南方程组的求解为基础，适用于有准确的河道地形和河床观测数据的河段，当这些资料条件缺乏时，水文学方法就成为洪水演算的另一种重要方法。 在水文学中，马斯京根（Muskingum）法是河道洪水演算的一种重要方法。马斯京根法于1934年提出，并在美国马斯京根河上首先应用。 圣维南方程组 Saint-Venant system of equations 描述水道和其它具有自由表面的浅水体中渐变不恒定水流运动规律的偏微分方程组。由反映质量守恒定律的连续方程和反映动量守恒定律的运动方程组成。1871年由法国科学家圣维南提出而得名。 圣维南方程组的形式 　一维单宽水流情况下，圣维南方程组的典型形式为：  式中t为时间；s 为距水道某固定断面沿流程的距离；h、v、Z0分别为相应于s处过水断面的水深、断面平均流速和水底高程；Hf为由于摩阻损失而引起的能量比降；g为重力加速度；t和s为自变量；h和v为因变量；Z0、Hf可由s、h和v确定。 (1)式为连续方程，反映了水道中的水量平衡，即蓄量的变化率(第一项)应等于沿程流量的变化率(第二项)。 (2)式为运动方程。其中第一项反映某固定点的局地加速度，第二项反映由于流速的空间分布不均匀所引起的对流加速度。以上两项称为惯性项。第三项反映由于底坡引起的重力作用，称为重力项。第四项反映了水深的影响，称为压力项。第三、四项可合并为一项，即水面比降。第五项为水流内部及边界的摩阻损失。 方程组求解原理 将连续方程式简化为河段的水量平衡方程式，将动力方程式简化为槽蓄方程式，则圣维南方程组简化为具有下列形式： Idt—Qdt=dW W＝f(Q，I) 马斯京根法依据的基本原理为水量平衡方程和槽蓄方程，其形式为 Δt(I-Q)=w2-w1　　　 （1）  w=k［x·I+（1-x）Q］　　（2） 　　式中：I—上断面入流（m3/s）； 　　　　Q—下断面出流（m3/s）； 　　　　　Δt—计算时段长（h）； 　　　　W—河段的槽蓄量（m3/s.h）； 　　　　　K—槽蓄系数（h）； 　　　　　X—流量比重因子。 应用中的一个关键问题是模型参数k、x的估计，传统的方法有试错法、最小二乘法等。 马司京根演算公式 在洪水演算时，联解水量平衡方程式和马司京干槽蓄曲线方程式得； Q2＝C0I2＋C1I1＋C2Q1 公式的推导 I1?t/2＋I2?t/2－Q1?t/2－Q2?t/2＝KxI2＋KQ2－KxQ2－KxI1－KQ1＋KxQ1 Q2?t/2＋KQ2－KxQ2＝I2?t/2－KxI2＋I1?t/2＋KxI1－Q1?t/2＋KQ1－KxQ1 (?t/2＋K－Kx)Q2＝(?t/2－Kx)I2＋（?t/2＋Kx）I1＋（－?t/2＋K－Kx）Q1 得到： Q2＝C0I2＋C1I1＋C2Q1 由上2式得 水量平衡式 槽蓄曲线 对于一个河段，只要确定参数K、x值及选定演算时段?t后，可以求出C0、C1、C2，就能根据上站流量过程I(t)及下站起始流量计算出下站的流量过程Q(t)。 马司京根流量比重x一般情况 流量比重因素一般为0－0.5。 河段蓄水作用甚大时，如湖泊、水库，x值接近于零。 随楔蓄影响的增加，x值渐增，而以x＝0.5为极限。此时入流与出流之权重相等。 一般天然河道的x值，大都在0.1－0.4之间，但也有例外。如出现负值时，就难以用楔蓄的观点来解释，需要进一步分析和论证x的物理实质。 * 等流时面积分配线 * 3）出口流量过程的计算 假定把分成5块等流时面积等 （图3－21），现有 3个时段的均匀净雨量，根据等流时线的概念，第一块等流时面积 上的净雨量 ，在第一时段内流到出口断面，则第一时段内平均流量为 ： 第二时段内流出的水体为 上第二时段净雨 和 上第一时段上的净雨量 ，即 ，则第二时段内的平均流量为 为： * 同理： 利用求得的