河道流量演算与洪水预报课件.ppt

1、第4章 河道流量演算与洪水预报圣维南方程组0LQtA连续性方程动量方程fSLVgVtVgtz1水面坡度局地惯性项迁移惯性项摩阻坡度惯性项01LQBtz或说明：底宽 b水面宽 B水深 h1m（边坡系数）A水面宽：过水断面面积：湿润周：水力半径：流量模数：3/222112)(12)(2ARnKmhbhbmhARmhbhbmhAmhbBn 河床粗糙度对固定河床均是水深h的函数 水力学模型的核心是圣维南方程的求解。圣维南方程是双曲拟线性偏微分方程，目前还无法求得其精确解析解，在实际应用中常采用数值近似解。数值近似方法主要有：-特征线法-直接差分-瞬时流态法-微幅波理论法-有限单元法特征线法：特征线法：

2、 这一方法是根据偏微分方程理论，先将基本方程变换为特征线的常微分方程组，然后对该微分方程进行离散，再结合初始条件和边界条件求数值解或图解。这种方法物理概念明确，数学分析严谨，计算结果精度较高。差分法：差分法： 将基本方程组直接离散化，进而联解由此得到的一组代数方程组。依据离散化时采用的数值格式不同，可将直接差分法分为显式差分法和隐式差分法显式差分法和隐式差分法两种。显式差显式差分法是根据前一时刻的已知值逐点分别求解下一时刻的未知值，计算过程简单分法是根据前一时刻的已知值逐点分别求解下一时刻的未知值，计算过程简单，但稳定性差，计算时间步长限制较多，步长较大时，计算可能不稳定，精度也难以保证；隐式

3、差分法不能直接由前一时刻求解下一时刻的值，必须同时对所有节点列出隐式差分法不能直接由前一时刻求解下一时刻的值，必须同时对所有节点列出差分方程而求解大型代数方程组，计算较为复杂差分方程而求解大型代数方程组，计算较为复杂，但稳定性好，计算时间步长可以取得较大，计算速度快。Preissmann计算方法-四点隐格式四点隐格式)(21)(2),(1111njnjnjnjfffftxfxffxffxfnjnjnjnj1111)1 (tfffftfnjnjnjnj21111下标代表空间步长上标表示时间步长为权重系数, (01)(x,t)令令 ,并记并记 ,则上式为：则上式为：)(21)(2),(11njnj

4、jjfffftxfxffxffxfnjnjjj11tfftfjj21Preissmann计算方法-四点隐格式四点隐格式0)()(0122KQQxzgAAQxtQxQBtz利用Preissmann格式，上式变为：0)()()(2211111xQQxQQBBBBtzzjjjjjjjjjj0)(|)(|2)1 ()(|)(|2)(1)()(2)(2)()()1 ()()(222111121111211111111111112121211121111njnjnjnjnjnjnjnjnjnjnjnjnjnjnjnjnjnjjjnjnjjjnjnjnjnnjnjnjnjjKQQAKQQAgKQQAKQQA

5、gzzxzzxAAgAAgAQAQxAQAQxtQQjj0)(|)(|2)1 ()(| )()(| )(2)(1)()(2)(2)()()1 ()()(222111121211111111111121221121111njnjnjnjnjnjnjnjjnjjnjjnjjnjjnjjnjjnjjnjnjnjjjnjnjjjnjnjnjnjnjjnjjnjjnjjKQQAKQQAgKKQQQQAAKKQQQQAAgzzxzzxAAgAAgAQAQxAAQQAAQQxtQQjjPreissmann计算方法-四点隐格式四点隐格式利用下面关系式上式线性化：)1 (1)1 (11jjnjjjnjjnjAA

6、AAAAAA利用利用在在的泰勒展开的泰勒展开)21 ()(1)1 ()(1)(12222njjnjnjjnjjnjKKKKKKKK2)1 (1)(xxf在在的展开的展开jnjnjjnjQQQQQ2)()(22jnjnjnjjnjjnjQQQQQQQQ|2| )(jnjjnjnjjzBzdzdAAjnjnjjzdzdKKjnjnjjzdzdBB0)()(1)(4)(111111xQQxQQBBBBBBtzzjjjjjjjjjjjjPreissmann计算方法-四点隐格式四点隐格式线性化后，连续性方程变为：0)()(1 4)(111111xQQxQQBBBBBBtzzjjjjjjjjjjjj0)(

7、)()(44)(112111111xQQxQQBBzdzdBzdzdBtBBtzzjjjjjjjnjnjjnjnjjjjjjjjjjjjjjEzDQCzBQA1111111)(411njnjjBBxtAnjnjnjnjnnjjdzdBBBxQQtBj2111)()(41)(411jjjBBxtCnjnjnjnjnjnjjdzdBBBxQQtD112111)()(41)()(4111njnjnjnjjQQBBxtEMake Presentation much more funWPS官方微博kingsoftwpsjjjjjjjjjEzDQCzBQA212122222)(|241njnjnjnjnj

8、jKQAtgAQxtAPreissmann计算方法-四点隐格式四点隐格式线性化后，动量方程变为：2)(|)()()()(2211222njnjnjnjnjnjnjnjnjnjnjnjnjnjnjjdzdKKABKQQtgzzgAAgAQBxtB2111112)(|241njnjnjnjnjjKQAtgAQxtC2)(|)()()()(2211111212112njnjnjnjnjnjnjnjnjnjnjnjnjnjnjjdzdKKABKQQtgzzgAAgAQBxtD)(|)(|)()( 2)( 2(2211111121212njnjnjnjnjnjnjnjnjnjnjnjnjnjnjnjjK

9、QQAKQQAtgzzAAgAQAQxtEPreissmann计算方法-四点隐格式四点隐格式通常边界条件分为以下通常边界条件分为以下3 3种：种：1. 给出水位变化过程：给出水位变化过程：2. 给出流量变化过程：3. 给出水位流量关系：以给出水位变化过程为例，方程组的矩阵形式如下：nmnmmmnnmmmmmmmmmmmzzEEEEEEzzzQzQzQzQDCBADCBADCBADCBADCBADCBA11211221221111111322111, 21, 21, 21, 21, 11, 11, 11, 122222222121212122121212111111111100000000000

10、0000000010至此，将复杂的偏微分方程组偏微分方程组化为了简单的线性代数方程组线性代数方程组这是偏微分方程运算的有限差分运算的基本流程和方法！这是偏微分方程运算的有限差分运算的基本流程和方法！随着遥感、计算机技术的进步，对圣维南方程的越来越高。根据对动力方程的不同简化，河道里的洪水波可分为：根据对动力方程的不同简化，河道里的洪水波可分为：1）运动波运动波2）扩散波扩散波3）惯性波）惯性波4）动力波）动力波1 1）运动波）运动波在动力方程中，对于山区性的河道，河底比降较大，在动力方程中，对于山区性的河道，河底比降较大，惯性项与附加比降项惯性项与附加比降项都可忽略。都可忽略。 0LQutQ特

11、点特点：水位：水位- -流量、流量流量、流量- -过水断面面积、波速过水断面面积、波速- -流量流量关系均为单一线；波速不变的条件下，流量在传播过关系均为单一线；波速不变的条件下，流量在传播过程中只位移而不衰减。程中只位移而不衰减。发生条件发生条件：只有在陡坡的情况下，才有可能：只有在陡坡的情况下，才有可能, ,而满足运动波的条件。而满足运动波的条件。 0SS 2 2）扩散波）扩散波在动力方程中，对于一般的天然河道水流，在动力方程中，对于一般的天然河道水流，惯性项较惯性项较其它项要小两个数量级，通常忽略其它项要小两个数量级，通常忽略。常用的。常用的流量演算流量演算水文学方法都忽略惯性项，且常将

12、动力方程简化为槽水文学方法都忽略惯性项，且常将动力方程简化为槽蓄方程，属于扩散波蓄方程，属于扩散波。 SSKLhSKQLL或或 式中式中-恒定流流量恒定流流量; ; -附加比降附加比降 S0-恒定流比降恒定流比降, ,一般可近似等于河底比降一般可近似等于河底比降 001SSQQ0QS扩散波的特点：扩散波的特点： 1）水位流量关系为多值函数关系。）水位流量关系为多值函数关系。 2）洪水在传播过程中，既要位移，又要坦化。）洪水在传播过程中，既要位移，又要坦化。 3）波速）波速 。流量流量Q和过水断面面积和过水断面面积A关关系有绳套，故对应某一传播流量的波速并非单值系有绳套，故对应某一传播流量的波速

13、并非单值AQu/3）惯性波）惯性波当当i=if=0时，即时，即水面比降为水面比降为0，没有摩阻损失，没有摩阻损失，水水深沿程变化完全是由惯性项引起的。深沿程变化完全是由惯性项引起的。4 4）动力波）动力波动力方程中各项均不忽略所描述的洪水波为动力动力方程中各项均不忽略所描述的洪水波为动力波波。对于受潮汐、闸、坝等严重影响的河段要用。对于受潮汐、闸、坝等严重影响的河段要用动力波进行演算。动力波进行演算。水量平衡方程和槽蓄方程水量平衡方程和槽蓄方程 LtAQtALQ0LLLtAQ00对连续性方程沿河长积分，可导出河段的水量平衡方程的微分形式：dttdWttWtOtIttWLAtLtAtItOtQt

14、LQQLLL)()()()()()()(), 0(),(000对河长L积分：dttdWttWtOtI)()()()(I,OtI(t)O(t)tdWI,OtI(t)O(t)tWt1t2河段水量平衡方程的差分形式：河段水量平衡方程的差分形式：I1I2Q1Q2122121)(21)(21WWtQQtII槽蓄方程槽蓄方程),( QIfW 河段的槽蓄量取决于和段中的水位沿程分布情况，即水面曲线的形状。但是，河段每一断面的水位与流量又存在一定的关系。当把河段的槽蓄量表示为入流量和出流量的函数时：称为河段的称为河段的槽蓄方程槽蓄方程I： 河段的入流流量；Q： 河段的出流流量；W：河段的蓄量当把河段的槽蓄量表

15、示为出流量的函数时：称为河段的称为河段的蓄泄方程蓄泄方程)(QfW 122121)(21)(21WWtQQtII水量平衡方程：槽蓄方程：)()(2211QfWQfW当已知河段入流量过程当已知河段入流量过程，根据水量平衡方程和槽蓄方程，即可求得，根据水量平衡方程和槽蓄方程，即可求得Q2值和值和W2值，值，对河段预报而言，对河段预报而言，Q2即为预报值；即为预报值；若逐时段连续计算，即可得到下断面的若逐时段连续计算，即可得到下断面的出流量过程出流量过程Q(t)矩形水槽稳定流： W=L*A=K*V*A=K*Q天然河道稳定流 H下Q0单一； H下W单一； QW单一天然河道不稳定流 出现绳套关系第二节第

16、二节 特特征河长法征河长法 特征河长特征河长（抵偿河长）的概念有前苏联著名水文学家加里宁和米留柯夫于1958年首次提出。 苏联水文学家。1916年11月10日生于巴库，1975年1月2日卒于莫斯科。1937年毕业于哈尔科夫水文气象学院，1951年获地理科学博士学位，自1954年起任教授，1970年当选为科学院通讯院士。19371942年，先后在国立水文研究所和哈尔科夫水文气象学院从事科研和教学工作。19421961年，在苏联水文气象总局中央预报研究所任高级研究员、水文预报研究处处长，兼任敖德萨水文气象学院教授。1961年在莫斯科大学地理系任教，1963年起担任陆地水文教研室主任。他是苏联科学院

17、水问题研究所创始人之一。他多年担任气象与水文杂志编委及水资源杂志副主编。 加里宁曾长期从事春汛和雨洪形成过程春汛和雨洪形成过程的基本研究，提出总入流概念，开辟了不依靠降水量资料计算产流量的途径，并创立用河网蓄水量和三角级数汇流曲线进行洪水预报的方法。1958年他与.米留柯夫共同发表特征河长概念，得出河槽非恒定流的近似计算方法，并应用分段连续演算的方法推求汇流曲线，于1963年进一步推导出河槽瞬时单位线。特征河长概念还在流域汇流计算和水位流量关系单值化等实际工作中得到应用。他对全球河川径流变化及水量交换的总规律进行过研究，提出了全球水文问题和水资源宏观管理的新课题。加里宁倡议并参加了应用空间信息

18、进行水文研究的工作。在苏联和国际水文界中，他较早倡议并参加应用电子计算机进行水文过程数学模拟的研究。 他的主要著作有短期水情预报方法原理短期水情预报方法原理(1952)、水体非恒定流近似计、水体非恒定流近似计算（合著，算（合著，1958）、水文预报（合著，）、水文预报（合著，1960）及全球水文学问题）及全球水文学问题(1968)。有些著作已在其他国家翻译出版。),(SzfQ由水力学可知，河段中任一断面的由水力学可知，河段中任一断面的流量是水位和水面流量是水位和水面比降的函数比降的函数：LZZS0S0+SS0-S上中下Z下Q 假设中断面水文不变；假设中断面水文不变； 漲洪时（蓝线），上断面先涨

19、，下断面后涨，漲洪时（蓝线），上断面先涨，下断面后涨，下断面水位比稳定水流降低下断面水位比稳定水流降低z，使，使得下断面的流量比稳定流时减少得下断面的流量比稳定流时减少；但由于这时水面比降比稳定流时；但由于这时水面比降比稳定流时增加了增加了S,这又会这又会使得通过下断面的流量比稳定流时增加使得通过下断面的流量比稳定流时增加；落洪时，由于上断面先落，下断面后落，落洪时，由于上断面先落，下断面后落，情况与涨洪时相反。情况与涨洪时相反。)(QfW 寻找这样一个河段长，在其下端面处，由于水位变化引起的流量变化正好与由于水面比降变化引起的流量变化相互抵偿，以致河段的槽蓄量与其下端面流量呈呈单值关系，即：

20、则该河长称为特征河长（抵偿河长）dsSQdzzQdQ对Q=Q(z,S）求全微分：根据特征河长定义：000)(0222,2,0QzSQLSSQSLzQSQSQSKQSLdzdQl 特征河长与河道的水力要素，即流量、比降和水位-流量关系坡度有关，是河道水力特征的综合参数；l河道的水力特征又决定了河道洪水波运动的特点基于特征河长的流量演算：基于特征河长的流量演算：在演算河段长等于特征河长时，假定蓄量W和出流Q存在线性关系。槽蓄方程：QKQfWl)(Kl为常数，特征河长的传播时间122121)(21)(21WWtQQtII2211QKWQKWll水量平衡方程：槽蓄方程：lllKtKtCKttCQCII

21、CQ5 . 05 . 0,5 . 05 . 0)(20121202第三节 马斯京根法马斯京根法是由G.T麦卡锡于1938年提出的，因首先应用于美国的马斯京根河而得名；我国从20世纪50年代开始对该方法进行深入的研究，并逐步进行改进：-1962年，华东水利学院（河海大学）提出马斯京根法有限差分解的河网单位线；-1985年，华东水利学院（河海大学）提出马斯京根法非线性解以及矩阵解。槽蓄方程5 . 0 , 0)1 ()1 (xQxxIQKQQxxIKW),(QIfW 对于一般的槽蓄是：而马斯京根方法假设槽蓄方程是：K表示蓄量流量关系曲线的坡度；x表示上、下断面流量在槽蓄量中的相对权重。如果河槽调蓄作

22、用大，则x小；反之，x大。例如，对水库而言，入流量不起作用，x=0;若入流与出流的影响相同，则x=0.5.马斯京根法的实际是马斯京根法的实际是：通过调整比重因子x，使流量与槽蓄量呈单一关系联解水量平衡方程与槽蓄方程)1 (),()(21)(21122121OxxIKOIfWWWtOOtII)1 ()1 ()(21)(2111222121OxxIKOxxIKtOOtII15 . 05 . 05 . 05 . 05 . 05 . 02102101211202CCCKxKtKxKtCKxKtKxtCKxKtKxtCOCICICO只要确定参数K,x,t，便可求出C0,C1,C2。根据上断面流量过程I(

23、t)及下断面起始流量。便可迭代计算下端面的流量过程Q(t)时间（月时间（月.日日 时：分）时：分）万县实测入流量万县实测入流量I0.26I20.48I10.26Q1宜昌演算出流宜昌演算出流Q27.1 14:0019900000228007.2 08:0024300631895525928217987.3 2:003880010088116645667274197.3 20:005000013000186247129387537.4 14:0053800139882400010076480647.5 8:0050800132082582412497515297.6 2:0043400112842

24、438413398490667.6 20:003510091262083212757427157.7 14:002690069941684811106349487.8 8:00224005824129129086278227.9 2:00196005096107527234230827.9 20:001790046549408600120063【例】 已知长江万县-宜昌河段的x=0.15，K=t=18h，求C0，C1，C2值。26. 05 . 05 . 048. 05 . 05 . 026. 05 . 05 . 0210KxKtKxKtCKxKtKxtCKxKtKxtCK，x参数的物理意义)1

25、 (KQOxxIKW为波速度) () (0QCQCLdQdWK槽蓄方程：从方程来看，就是调整x，使得河段蓄水量W与示蓄流量Q成单一的曲线，K就是这个曲线的斜率： Q： Q - W等价于等价于Q0 W，相当于河槽蓄相当于河槽蓄量量W下的恒定流流量下的恒定流流量 K：等价于恒定流状态下的河段传播时间，等价于恒定流状态下的河段传播时间，即即 dQdWK经过分析推导，可以得到：Llxx21x1-水面线形状参数，反映了楔蓄的大小。当水面为直线 x1=1/2;l-特征河长000)(QzSQlLlLlxx22121000)(QzSQl由上面公式可以看出：由上面公式可以看出：1. x由两部分组成，x1代表水面

26、曲线形状，反映了楔蓄的影响；L/l为按特征河长划分的河段数，反映了河槽的调蓄主要2.由于l0，故xL,x0.3. 在上游河道，S0较大，l较小，河道的调蓄能力小，x较大；4. 在下游河道，S0较小，l较大，河道的调蓄能力大，x较小；5. 对于一般的河道，上游的x较大；下游的x较大。K，x的推求方法 试算法：试算法：假定假定x值值 相应相应Q值值 Q W关系关系，取最密集于一直线时的取最密集于一直线时的x值值，量取关系线斜率为量取关系线斜率为K值值 分析法：分析法： 数学优化方法数学优化方法uLKLlx221【例】 已知长江万县-宜昌河段的洪水实测资料，求K,x值。时间时间（月（月.日日 时：分

27、）时：分）万县实测入万县实测入流量流量I宜昌演宜昌演算出流算出流Q区间区间径流径流量量q修正实测修正实测出流量出流量Qr=Q-qQ=I-QrQSQ=Qr+x(I-Qr)x=0.10 x=0.25x=0.157.1 14:00199007.2 08:002430023700600231001200730002322023400232807.3 2:003880027000160025400134001340073002674028750274107.3 20:005000037800120036600134009850207003794039950386107.4 14:005380048400

28、9004750063002850305504813049075484457.5 8:00508005190050051400-600-3200334005134051250513107.6 2:00434004960040049200-5800-6650302004862047750483307.6 20:00351004300040042600-7500-7900235504185040725414757.7 14:00269003560040035200-8300-7450156503437033125339557.8 8:00224002930030029000-6600-5450820

29、0283402735028010试算法：试算法：【例】 已知长江万县-宜昌河段的洪水实测资料，求K, x值。时间时间（月（月.日日 时：分）时：分）万县实测入万县实测入流量流量I宜昌演宜昌演算出流算出流Q区间区间径流径流量量q修正实测修正实测出流量出流量Qr=Q-qQ=I-QrQSQ=Qr+x(I-Qr)x=0.10 x=0.25x=0.157.1 14:00199007.2 08:002430023700600231001200730002322023400232807.3 2:003880027000160025400134001340073002674028750274107.3 20:

30、005000037800120036600134009850207003794039950386107.4 14:0053800484009004750063002850305504813049075484457.5 8:00508005190050051400-600-3200334005134051250513107.6 2:00434004960040049200-5800-6650302004862047750483307.6 20:00351004300040042600-7500-7900235504185040725414757.7 14:00269003560040035200

31、-8300-7450156503437033125339557.8 8:00224002930030029000-6600-545082002834027350280107.9 2:00196002420030023900-4300-2150275023470228252325500.511.522.533.522.533.544.555.5S(104m3/s)Q(!04m3/s) x=0.10 x=0.25x=0.15dQdWK=dW/dQ=18h则由x=0.15,k=18h，t=18h。可由：26. 05 . 05 . 048. 05 . 05 . 026. 05 . 05 . 0210K

32、xKtKxKtCKxKtKxtCKxKtKxtC利用迭代公式算出出流量过程！分析法：马斯京根法的参数（K,x）及其槽蓄方程Q-W关系的物理意义已经明确，就为推求参数提供了水力学方法。Llx221000)(QzSQlQ【例】 沅水沅陵站至王家河站河段（L=112Km），用水力学方法推求K，x。v(m/s)Q(m3/s)王家河沅陵C0=1.66VK=l/C0第4节 河道相应水位（流量）预报相应水位相应水位（流量）：（流量）：在河段同次洪水过程线上，处于同一位在河段同次洪水过程线上，处于同一位相点上、下站的水位相点上、下站的水位（流量）（流量） 相应水位（流量）预报相应水位（流量）预报：利用某时刻上

33、游站利用某时刻上游站的水位（流量）预报一定时间（如传播时间的水位（流量）预报一定时间（如传播时间 ）后下游站的水位（流量）：后下游站的水位（流量）：)QQ(fQttt，区间因素，上，下，下，)QQ(ftt上，下，1 1、基本原理、基本原理1、上、上、下游流量基本关系下游流量基本关系QL：上下站相应流量的差值，它随上、下站流量的大小和附加比降不同而异，其实质是反映洪水波变形中的坦化作用。 在制定相应水位法的预报方案时,一般采取水位过程线上的特征点,如洪峰等,作出该特征点的相应水位关系曲线与传播时间曲线。 简单的相应水位法简单的相应水位法 在无支流汇入的河段上,若影响洪水波传播的因素比较单纯,可用

34、简单的相应水位法。- 方案：根据上游站和下游站的实测水位过程线,摘录相应的特征点即洪峰水位值及其出现时间(见表),并绘制相应洪峰水位相关曲线及其传播时间曲线作为预报方案。如果已知上游洪峰水位可以根据该曲线预测如果已知上游洪峰水位可以根据该曲线预测出洪出洪峰到达下游时间峰到达下游时间、流量（水位）流量（水位）Q Q图 23上游-下游洪峰传播关系曲线(b)(a涪江涪江桥（上）左家岩（下）洪峰相应水位关系曲线)(b(a)(h1AO21OI1LQO(m3/s)m,t+m,t/s)3(mI（2）以下游站同时水位为参数的相应水位法）以下游站同时水位为参数的相应水位法下游站同时水位下游站同时水位Z Z下下,

35、t,t：上游站水位Z上,t出现时刻的下游水位,它与Z上,t一起能反映t时刻的水面比降变化；同时，也间接地反映区间入流和断面冲淤以及回水顶托等因素的影响。-预报方案预报方案：制作预报方案时，以下游站同时水位Z下,t为参数作等值线,分别绘制Z上,t Z下,t Z下,t+和Z上,t Z下,t 相关曲线,如图。-作业预报作业预报：按t时刻的水位Z上,t及Z下,t，按图箭头方向查得Z下,t+和，从而预报出t+时下游将要出现的洪峰水位Z上,t+ 由图中可看出，上游水位上游水位z上上，t 、下游同时水位、下游同时水位z下，下，t这两个因素反映了河道内水面坡度水面坡度 和河槽内涨河槽内涨水前的底水水前的底水

36、以下游同时水位下游同时水位Z下，下，t为参数预报方案：图中：以下游洪峰水位Zm下，t+为横轴，以上游洪峰水位Zm 上，t为纵轴，以下游同时水位Z下，t为参数ttmtmZZfZ,下上下以上游站涨差为参数的水位相关法以上游站涨差为参数的水位相关法上述各种洪峰水位预报方案,可近似地用来预报下游站的洪水过程。但由于它们没有反映洪水过程中附加比降的变化等因素,使预报的洪水过程常常有比较大的系统误差。为克服这种缺点,可用以上游站水位涨差为参数的水位相关法。洪水波通过某一断面时,波前的附加比降为正,使涨水过程的涨率Z上/t为正;波后的附加比降为负,使落水过程的涨率为负。因此, 水位(流量)过程线的涨(落)率

37、在很大程度上反映了附加比降和水面比降的大小。以涨率涨率 为参数预报方案：图中：以下游洪峰水位Zm下，t+为横轴横轴，以下游同时水位Z下，t 为纵轴纵轴，以上游总涨差Z上为参数参数tZttmZZfZ,下上下以支流的相应水位支流的相应水位为参数预报方案：it , i ,t , 1,上上下ZZfZmtm洪水过程预报 摘取相应水位（除峰、谷等特征点以外）的方法 1 洪波展开法 2 为常数的方法 3 现时校正法 1.洪波展开法2. 为常数的方法可以证明： 是水面坡度S的函数tZ,下实际中，常以河段平均传播时间 摘取 上、下游站的“相应水位” 预报方案为：),(t ,t ,下上下ZZfZt 两者精度是相同

38、的t ,下上下ZZfZttt ,下上下ZZfZtt也可用时段水位涨差为参数建立预报方案 预报方案为：),(t ,t ,下上下ZZfZt3. 现时校正法 具体做法: 作业预报时,先按平均线得预报值,参照前一时 刻下游站实测值偏离预报值得差值,修正按平均线所得的预报值现时校正法，目前在作业预报中广泛使用 前一时刻,预报值A，实测值B，误差AB。那么,在本时刻查得预报曲线为C点时,就要参照AB的误差支流汇入-预报方案预报方案：上图是长江万县水文站宜昌水文站河段以Q上为参数的水位预报方案。-作业预报作业预报：已知t时刻的Z上,t(或Z下,t)、Z上(或Q上)，在图上查出预报的下游水位Z下,t+。tin

39、ititQQfQ,1,)(上上下（4 4）合成流量法）合成流量法 在有支流汇入的河段，按照上游干、支流各站的传播时间，把各站同时刻到达下游站的流量叠加起来得合成流量，然后建立合成流量与下游站相应流量的关系曲线，进行预报的方法称为合成流量法。该法预报下游站流量的关系式为: ：为上游干、支流各站相应流量；i:上游干、支流各站到下游站的洪水传播时间；n:上游干、支流的测站数目。根据上式的关系该法的预见期取决于上游各站中传播时间最短的一个。一般情况下,上游各站中以干流站的流量为最大,从预报精度的要求出发,常常用它的传播时间作为预报方案的预见期。第五节 回水和感潮河段的水位（流量）预报受变动回水顶托影响

40、河段的水位（流量）预报在干、支流会河、湖汇合附近的河段，上游来水与支流或湖泊来水之间在干、支流会河、湖汇合附近的河段，上游来水与支流或湖泊来水之间相互干扰，常发生回水顶托，影响洪水波运行变化特征。建立预报方案相互干扰，常发生回水顶托，影响洪水波运行变化特征。建立预报方案是，要分析上游来水和回水顶托这两项因素及其作用程度。是，要分析上游来水和回水顶托这两项因素及其作用程度。对于上游干流来水影响为主的河段，可先建立上、下游站相应水位（流对于上游干流来水影响为主的河段，可先建立上、下游站相应水位（流量）关系，用反映回水顶托的要素作为参数；量）关系，用反映回水顶托的要素作为参数；如支流发生大流量洪水时对干流有回水顶托的影响，除建立干流河段上、如支流发生大流量洪水时对干流有回水顶托的影响，除建立干流河段上、下游站相应水位关系外，还应建立支流来水量与回水影响量之间的关系，下游站相应水位关系外，还应建立支流来水量与回水影响量之间的关系，用以修正预报值。用以修正预报值。感潮河段的水位预报 感潮河段：受海洋潮汐影响的入海河流的下游河段。 感潮河段的水情变化除与上游来水有关外，还受潮汐波动的影响，若遇台风和暴雨，常使感潮河段水位抬高