由流量资料推求设计洪水

1、第9章 水文预报内容简介研究对象本章研究水文现象的客观规律，利用现时已经掌握的水文、气象资料，预报水文要素未来变化过程。研究内容1.短期洪水预报；2.枯水预报；3.施工水文预报； 4.水文实时预报方法。研究目的在防汛工作中，及时准确的水文预报，是防汛抗洪指挥决策的重要科学依据；在水能、水资源的合理调度、开发利用和保护以及航运等工作中，都需要有水文预报作指导。第9.1节 概 述内容提要1. 水文预报的重要作用；2. 水文预报的分类； 3. 水文预报工作的基本程序学习要求掌握预见期的定义及水文预报工作的基本程序。水文预报的重要作用可靠的洪水预报对防止洪水灾害具有特别重要的作用。例如在河流防洪抢险中

2、，需要及时预报出防洪地点即将出现的洪峰水位、流量，以便在洪峰到来之前，迅速加高加固堤防、转移可能受淹的群众和物资，动用必要的防洪设施等，把洪水灾害减小到最低限度。图为1998年长江沙市水位预报与实测情况。图9.1.1 1998年长江沙市水位预报与实测情况在水库管理中，可以利用洪水预报，使上游来的洪水与区间洪水的高峰段彼此错开（称错峰），即下游洪水很大时，水库把上游来的洪水暂时蓄存起来，待下游洪峰过后，再加大水库泄量，把上游来的洪水放出来，从而大大减低下游的洪峰和洪水灾害，例如1998年8月长江中下游发生近百年一遇的特大洪水，由于及时准确的洪水预报，对葛洲坝水库、隔河岩水库和漳河水库科学调度，使

3、三峡以上来的洪水和清江、沮漳河洪水的洪峰互相错开，大大降低了荆江河段的洪峰水位，避免了荆江分洪损失，为战胜该年发生的特大洪水做出了巨大贡献。表为1998葛洲坝水库、隔河岩水库在错峰、调峰中，降低沙市水位发挥作用的分析结果。表葛洲坝水库、隔河岩水库在错峰调度对沙市水位的影响另外，洪水预报还可较好地解决水库防洪与兴利的矛盾，在预报的洪水未进库之前，先打开泄洪闸门腾空一部分库容，以便洪水来临时能蓄存更多的水量；当洪水即将结束时，预知近期没有很大的洪水入库，则可超蓄洪水尾部的一些水量，用于多发电、多灌溉，使现有工程发挥更多的效益。水文预报的分类1.按预报的项目，水文预报可分为l 径流预报：预报的要素主

4、要是水位和流量，水位预报指的是水位高程及其出现时间；流量预报则是流量的大小、涨落时间及其过程。径流预报又可分洪水预报和枯水预报。l 冰情预报：冰情预报是利用影响河流冰情的前期气象因子,预报流凌开始、封冻与开冻日期,冰厚、冰坝及凌汛最高水位等。l 沙情预报：沙情预报则是根据河流的水沙相关关系,结合流域下垫面因素,预报年、月和一次洪水的含沙量及其过程。l 水质预报：预测河流中污染物迁移转化的时空变化过程。l 施工水文预报：在工程施工期间要进行的特殊预报项目。2. 按预见期的长短，水文预报可分为l 短期水文预报：主要由水文要素作出的预报l 中长期水文预报：包括气象预报性质在内的水文预报预报的预见期是

5、指发布预报与预报要素出现的时间间距。在水文预报中，预见期的长与短并没有明确的时间界限。9.1.3水文预报工作的基本程序水文预报工作大体上分为两大步骤。(1) 制定预报方案：分析预报要素的形成规律，建立由过去的观测资料推算水文预报要素大小和出现时间的一整套计算方法，即水文预报方案。制定的方案按国家水文情报预报规范要求的允许误差进行评定和检验。只有质量优良和合格的方案才能付诸应用。图为我国2000年6月30日开始实施的水文情报预报规范（SL250-2000）。图 我国现行的水文情报预报规范（SL250-2000）(2) 进行作业预报:将现时发生的水文气象信息，经过预报方案算出即将发生的水文预报要素

6、大小和出现时间，及时将信息发布出去，这个过程称为作业预报。若现时水文气象信息是通过自动化采集、自动传送到预报中心的计算机内,由计算机直接按存储的水文预报模型程序计算出预报结果。这样的作业预报称为联机作业实时水文预报。复习思考题1.水文预报方案的评定和检验，是d a、都用制定预报方案时的全部实测资料进行b、都用作业预报过程中的时实测资料进行c、用以上两种资料进行d、用制定预报方案中使用的资料来进行评定，而用没有参加方案编制的预留资料进行检验。2. 洪水预报的预见期是c a、洪峰出现的时间b、洪水从开始到终止的时间c、从发布预报时刻到预报洪水出现时刻所隔的时间d、从收到报汛站资料开始到预报的洪水出

7、现所隔的时间第9.2节 短期洪水预报内容提要1. 河段中的洪水波运动；2. 相应水位(流量)法；3. 合成流量法；4. 流量演算法；5. 降雨径流预报。学习要求掌握相应水位(流量)法、合成流量法和马斯京根流量演算方法。河段中的洪水波运动1. 洪水波流域上大量降水后,产生的净雨沿坡地迅速汇集,注入河槽, 由于降雨量时空分布不均匀、河网干支流和分布形状的不同，以及水流汇集速度的快慢，河道接纳的水量沿程不同，使河道沿程水面发生高低起伏的一种波动,称为洪水波。i0ii图9.2.1 河道洪水波水面比降示意图2. 附加比降附加比降是洪水波的主要特征之一。附加比降是指洪水波水面比降与同水位稳定流水面比降之差

8、，即。图为河道洪水波水面比降示意图。l 当涨洪时： ；l 当落洪时： ；l 当水流稳定时： 。3.河段洪水波的传播与变形由于河槽的调蓄作用，洪水波向下游传播过程中,不断发生变形，如图（动画）所示。在沿棱柱形河槽运动中其变形有两种形态：图河道洪水波传播与变形过程示意图l 洪水波展开：洪水波在传播过程中,波长不断加大,波高不断减小的现象称为洪水波的展开,即A2C2>A1C1,h2< h1。l 洪水波扭曲：洪水波在传播过程中,波前水量不断向波后转移的现象称为洪水波的扭曲。在自然河道中，河道断面边界的差异和河段区间入流等条件变化，都对洪水波变形有显著的影响。 相应水位(流量)法根据河段洪水

9、波运动和变形规律,利用河段上断面的实测水位(流量),预报河段下断面未来水位(流量)的方法,称为相应水位(流量)法。1.基本原理相应水位（流量）相应水位（流量）是指在河段同次洪水过程线上，处于同一位相点上、下站的水位（流量）称为相应水位（流量）。如图9-2-3所示某次洪水过程线上的各个特征点,例如上游2点洪峰水位经过河段传播时间,在下游站2点的洪峰水位,就是同位相的水位。图9-2-3 上、下游站相应水位过程线图相应水位(流量)法的基本方程l 河段无区间入流：设河段上游站流量为Q上,t，经过时间的传播,下游站的相应流量为Q下,t+,两者的关系为：Q下,t+=Q上,t-QL (9-2-1)式中QL-

10、洪水波展开量,与附加比降有关。l 河段有区间入流q两者的关系为：Q下,t+=Q上,t-QL + q (9-2-2)2.无支流河段的相应水位预报在制定相应水位法的预报方案时,一般采取水位过程线上的特征点,如洪峰等,作出该特征点的相应水位关系曲线与传播时间曲线。简单的相应水位法在无支流汇入的河段上,若影响洪水波传播的因素比较单纯,可用简单的相应水位法。l 预报方案：根据上游站和下游站的实测水位过程线,摘录相应的特征点即洪峰水位值及其出现时间(见表9-2-1),并绘制相应洪峰水位相关曲线及其传播时间曲线(图9-2-4)作为预报方案。l 作业预报：按t时上游出现的洪峰水位Z上,t，在Z上,tZ上,t+

11、曲线上查得Z下,t+,在Z上,t曲线上查得,从而预报出t+时下游将要出现的洪峰水位Z上,t+。表9-2-1 长江万县站三斗坪站相应洪峰水位及传播时间 图9-2-4 万县站三斗坪站相应洪峰水位及传播时间关系曲线图以下游站同时水位为参数的相应水位法下游站同时水位Z下,t就是上游站水位Z上,t出现时刻的下游水位,它与Z上,t一起能反映t时刻的水面比降变化,同时也间接地反映区间入流和断面冲淤以及回水顶托等因素的影响。l 预报方案：制作预报方案时，以下游站同时水位Z下,t为参数作等值线,分别绘制Z上,t Z下,t Z下,t+和Z上,t Z下,t 相关曲线,如图9-2-5所示。l 作业预报：按t时刻的水位

12、Z上,t及Z下,t，按图9-2-5的箭头方向查得Z下,t+和，从而预报出t+时下游将要出现的洪峰水位Z上,t+。图9-2-5 以下游站同时水位为参数的水位及传播时间关系曲线图以上游站涨差为参数的水位相关法上述各种洪峰水位预报方案,可近似地用来预报下游站的洪水过程。但由于它们没有反映洪水过程中附加比降的变化等因素,使预报的洪水过程常常有比较大的系统误差。为克服这种缺点,可用以上游站水位涨差为参数的水位相关法。图9-2-6 长江万县水文站宜昌水文站以上游站涨差为参数的水位预报方案洪水波通过某一断面时,波前的附加比降为正,使涨水过程的涨率Z上/t为正;波后的附加比降为负,使落水过程的涨率为负。因此,

13、 水位(流量)过程线的涨(落)率在很大程度上反映了附加比降和水面比降的大小。l 预报方案：图9-2-6是长江万县水文站宜昌水文站河段以Q上为参数的水位预报方案。l 作业预报：已知t时刻的Z上,t(或Z下,t)、Z上(或Q上)，在图上查出预报的下游水位Z下,t+和预见期为。合成流量法在有支流汇入的河段，按照上游干、支流各站的传播时间，把各站同时刻到达下游站的流量叠加起来得合成流量，然后建立合成流量与下游站相应流量的关系曲线，进行预报的方法称为合成流量法。该法预报下游站流量的关系式为: (9-2-3)式中，为上游干、支流各站相应流量；i为上游干、支流各站到下游站的洪水传播时间；n-上游干、支流的测

14、站数目。根据式(9-2-3)的关系该法的预见期取决于上游各站中传播时间最短的一个。一般情况下,上游各站中以干流站的流量为最大,从预报精度的要求出发,常常用它的传播时间作为预报方案的预见期。流量演算法流量演算法是在圣维南方程组简化的基础上,利用河槽的水量平衡方程替代连续性方程,用河段的蓄泄关系替代动力方程, 然后联立求解,将河段的入流过程演算为出流过程的方法。1.基本原理河段流量演算是由以下两个基本公式组成：l 河槽时段水量平衡方程 （9-2-4） 若当河段内有区间入流量q，将q值并入到上断面的入流量中进行演算,即: l 河段蓄水量与泄流量方程 S=f(Q) （9-2-5） 式中S为河段内某一流

15、量所对应的蓄水量。 图9-2-7 河段时段水量平衡示意图2.马斯京根法及其槽蓄曲线方程该方法由G.T.麦卡锡(G-T-McCarthy) 1938年应用于美国马斯京根河流域而得名。该法主要是建立马斯京根槽蓄曲线方程,并与水量平衡方程联立求解,进行河段洪水演算。图9-2-8 河段槽蓄量示意图在马斯京根槽蓄曲线方程中,河段槽蓄量由柱蓄和楔蓄两部分组成，如图9-2-8。令x为流量比重因素,SQ.上、SQ.下分别为上下断面在稳定流情况下的蓄量,它们与河段内总蓄量S的关系为：（9-2-6）若将河道中的断面流量与相应的槽蓄量近似地按稳定流处理时（SQ =KQ）,由式（9-2-6）可得到马斯京根槽蓄曲线方程

16、：（9-2-7）式中,K为稳定流情况下的河段传播时间，Q称为示储流量：。3.马斯京根流量演算方程联解水量平衡方程式（9-2-4）和马斯京根槽蓄曲线方程式（9-2-7）,可得马斯京根流量演算方程:（9-2-8）式中： （9-2-9）4. 试算法确定x、K值对某一次洪水,假定不同的x值计算Q,作出Sf(Q)关系曲线,其中能使二者关系成为单一直线的x值即为此次洪水所求的x值,而该直线的斜率即为所求的K值。取多次洪水作相同的计算和分析,就可以确定该河段的x、K值。【例9-2-1】根据某河段一次实测洪水资料(表9-2-2),用马斯京根法进行河段洪水演算。表9-2-2 马斯京根法S与Q值计算表根据河段和资

17、料情况,取时段长t=12h。将河段实测洪水资料列于表中的第(2)(4)栏。因区间无实测值,将河段入流总量与出流总量差值作为区间入流总量,其流量过程近似地按入流过程的比值分配到各时段中去。按水量平衡方程式,分别计算各时段槽蓄量S(表中第(7)栏),然后逐时段累加S得槽蓄量S(表中第(8)栏)。假定x值,按Q=xQ上+(1-x)Q下计算Q值。本例分别假设x=0.2和x=0.3,计算结果列于表中第(9)、(10)栏。按第(9)、(10)两栏的数据,分别点绘两条SQ关系线,其中以x=0.2的关系线近似于直线(见图9-2-9),该x值即为所求。该直线的斜率K=S/Q=800×12/800=12

18、h。图9-2-9 马斯京根法SQ关系曲线图假定本河段的K值已选定,即x=0.2,K=12h以及t=12h,代入式（9-2-9）得: Co=0.231,C1=0.538,C2=0.231，而且Co+C1+C2=1.0,计算无误。因此,该河段的洪水演算方程为:Q下,2=0.231Q上,2+0.538Q上,1+0.231Q下,1根据本河段另一场洪水的上断面流量资料(见表9-2-3中第(2)栏),用上述洪水演算方程,可算出河段下断面的流量(见表9-2-3第(6)栏)。4.马斯京根法中几个问题的分析(1) K值的综合l K具有时间的因次,它基本上反映河道稳定流时河段的传播时间。l 在不稳定流情况下,按流

19、量分级,根据不同的流量取不同的K值。(2) x值的综合l 流量比重因素x反映楔蓄对流量的作用外,还反映河段的调蓄能力。l 对于一定的河段,x在洪水涨落过程中基本稳定,但也有随流量增加而减小的趋势。l 流量比重因素x，一般从上游向下游逐渐减小,介于0.20.45之间,特殊情况下也有小于零的。(3) 计算时段t的选择l t不能取得太长,以保证流量过程线在t内近于直线;l 为在计算中不漏掉洪峰,选取的t最好等于河段传播时间。(4) 预见期问题l 马斯京根法用于预报时是没有预见期的。l 如果取t=2kx，则Co=0,有一个时段的预见期。表9-2-3 马斯京根洪水演算表降雨径流预报降雨径流预报是利用流域

20、降雨量经过产流计算和汇流计算,预报出流域出口断面的径流过程。主要包括两方面的内容:l 由降雨量推求净雨量;l 由净雨过程推求流域出口断面的径流过程。有关这两个问题的计算原理和计算方法,见第七章的内容。在降雨径流预报中，有以下主要工作内容：1.编制降雨径流方案根据流域自然地理特征和实测资料条件,建立流域产、汇流计算方案,如降雨径流相关图、单位线等。并对方案的预报精度进行评定和检验。2.作业预报在作业预报中,如图9-2-10所示 to时刻以前的降雨量是实测的； to之后的降雨量是预报的。预报雨量与预报方案的误差,两者都影响预报的精度。因此,在作业预报时,应根据实测时段降雨量或实测流量对预报的径流进

21、行逐时段修正。图9-2-10降雨径流预报法预报洪水过程示意图（动画）复习思考题1.洪水波的波速时洪水波上b a、某一断面的平均流速b、某一位相点的运动速度c、最大断面平均流速d、某一水质点的运动速度2.对于一定的河段，附加比降影响洪水波变形，就洪峰而言，下述哪种说法是对的a a、附加比降愈大，洪水波变形愈大，传播时间愈短b、附加比降愈大，洪水波变形愈大，传播时间愈长c、附加比降愈大，洪水波变形愈大，对时间没什么影响d、附加比降愈大，洪水波变形愈小，对时间没什么影响3.引起河槽洪水波变形的原因是d a、附加比降 b、区间径流 c、河网调蓄及断面冲淤 d、以上因素都有影响4.存在附加比降是洪水波的

22、主要特征，在涨洪时ba、=0 b、>0 c、<0 d、不能肯定5.利用相应水位法作预报方案，加入下游站同时水位作参数的目的主要是c a、考虑洪水到来前的底水作用b、考虑上游降雨的影响c、考虑水面比降和底水的作用d、考虑河道坡降的作用6.相应水位法由上断面洪峰水位预报下断面洪峰水位的预见期为两断面间的b a、洪水上涨历时 b、洪水传播时间 c、整个洪水历时 d、洪水退落历时7.用马斯京根法作洪水预报，当计算时段=（，为该法的两上参数）时预见期为aa、 b、 c、零 d、+8.马斯京根流量演算方程中的系数da、大于零 b、小于零 c、等于零 d、三种情况都可能存在9.马斯京根法中的值，

23、从理论上说，应b a、随流量的增大而增大b、随流量的增大而减小c、不随流量变化d、以上三种情况都有可能10.在有支流汇入的河段，采用合成流量法预报下游站时刻的流量，它的计算表达式（设为上游第个站的流量，为上游第个站到下游站的洪水传播时间）为a a、= b、= c、= d、= 11.河段涨洪时，同时刻上游站洪峰流量与下游站洪峰流量（假设无区间入流）的关系是a a、> b、= c、< d、不能肯定第9.3节 枯水预报内容提要1. 基本原理；2. 退水曲线法；3. 前后期径流量相关法。学习要求掌握退水曲线法和前后期径流量相关法。9.3.1基本原理枯水期河中的流量主要是由滞留在流域中的蓄水

24、量消退形成，其次是来源于枯季 的降雨。流域蓄水量包括地面、地下蓄水量两部分：l 地面蓄水量存在于地面洼地、河网、水库、湖泊和沼泽之中；l 地下蓄水量存在于土壤孔隙、岩石裂隙、溶隙和层间含水带之中。由于地下蓄水量的消退比地面蓄水量慢得多，故长期无雨后河中水量几乎全由地下水补给，流域的水量平衡方程式和蓄量方程式分别为：（9-3-1）（9-3-2）式中Q(t)为枯水期中t时刻的流域出流量（m3/s）；W（t）为枯水期中t时刻的流域蓄水量（m3）；K为流域退水参数（s）。联立求解式（9-3-1）和（9-3-2），得枯水期流量消退规律的表达式为（9-3-3）式中Q(0)为枯水期某一初始时刻的流域出流量（

25、m3/s）；为消退系数。因此，只要分析出流域Kr值，就可掌握该流域的退水规律。9.3.2 退水曲线法在枯水预报中常把退水曲线表示成相邻时间t的流量的相关关系，由式（9-3-3），有（9-3-4）Kr随所采用的时段长而变，当相邻时间t取一固定时段时，相邻流量的比值则为： （9-3-5）式中消退系数Kr，可通过建立相邻时间流量的相关图（如图9-3-1）推求。图9-3-1 清江搬鱼咀站枯水期Qt+tQt关系图【例9-3-1】表9-3-1为清江搬鱼咀站1972年1月枯水期的流量摘录，相邻时间t为2d。点绘相邻流量关系如图9-3-1。由相关分析，得相关方程中a=0,则有表9-3-1 清江搬鱼咀站1972

26、年1月枯水期流量资料9.3.3 前后期径流相关法前后期径流相关法是根据流域前期径流量来预报未来的后期径流量的。该法的预见期一般较长，如10天、1个月等，可以直接预报出时段枯水径流量。图9-3-2为滏阳河东武仕站11月平均流量与10月平均流量的关系线（一年一个点子），关系很好。图9-3-2滏阳河东武仕站关系线复习思考题1.枯水径流变化相当稳定，是因为它主要来源于ba、地表径流 b、地下潜水 c、河网蓄水 d、融雪径流2. 枯水预报中的退水曲线法与前后期径流量相关法，预报的流量值是有差别的，预见期越长其差别越明显。T第9.4节 施工水文预报内容提要1. 坝址处流量预报；2. 围堰上、下游水位预报；

27、3. 截流期水情预报学习要求掌握围堰水情预报和截流期水情预报方法。施工水文预报是指对水电工程施工期，受到施工回水影响河段的水文预报。按施工阶段，施工水文预报主要分为围堰水情预报和截流期水情预报。围堰水情预报在修筑围堰（图9-4-1）及导流建筑物阶段，要求预报围堰前的水位或流量，以防止河水漫入施工区。图9-4-1 施工围堰平面示意图 1预报坝址处的流量 采用马斯京根流量演算法计算坝址处的流量，但此时要注意的是围堰修建以后，天然河道情况下的马斯京根槽蓄曲线已不适用了。具体计算方法步骤如下： 采用水力学方法计算各级稳定流量Qi相应的水面曲线； 计算出上游为入库站，下游为坝址各级稳定流量Qi的槽蓄量W

28、i； 假定修筑围堰后，原马斯京根参数x值不变（修围堰后x值最好以实测资料分析而得），计算出示储流量Q，点绘WQ关系，推求出K值。求得修筑围堰后的演算公式，由于围堰上、下游两端距离很短，推算的流量可作为围堰上、下游的流量。 2. 围堰上、下游水位预报 修筑围堰后，围堰上游天然情况下的水位流量关系发生了变化，此时应重新建立上游水位流量关系曲线。根据坝址流量，可用下列公式近似推求束窄河段水位的壅高值Z，（9-4-1）（9-4-2）（9-4-3）式中，Z上、Z下分别为上、下游断面水位；V上、Vc分别为上游及束窄断面平均流速；A上、Ac分别为上游及束窄处断面面积；Q为稳定流量，为动能修正系数，取1.0

29、1.1。g为重力加速度。采用试算法计算Z，在计算时，要求具备有下游断面的水位流量关系Qf(Z下),上游及束窄断面的水位面积曲线A上f1(Z上)，Ac=f2(Z下)。具体计算方法步骤如下：拟定过水流量Q，查Qf(Z下)曲线得Z下；由Z下值查Ac=f2(Z下)曲线得Ac，由此计算出Vc，并算出；假定壅水高度Z，则得上游水位Z上=Z下+Z，由A上f1(Z上)曲线查得A上，计算出V上，并计算出； 按式（9-4-1）计算壅水高度Z，若计算出的Z与假定的Z相符，则试算完毕，否则重新计算。计算出各级流量的壅水高度，即可建立上游壅高后的水位流量关系曲线Qf(Z下+Z)= f(Z上)，见图（9-4-2）。围堰下

30、游的水位流量关系仍是天然情况下的，即Qf(Z下)。有了围堰上、下游水位流量关系，便可利用前面预报的流量Q，推求出上游水位Z上、和Z下，完成围堰上、下游的水位预报。图9-4-2 围堰上、下游的水位流量关系曲线 截流期水情预报水利水电工程施工的截流一般是在枯季进行，枯季河水流量小，流速慢给截流施工创造了有利条件。只有预先掌握了截流期河道流量的大小，采取相应措施，施工截流才能顺利的进行 。因此，截流期水情预报是施工截流中不可缺少的工作。截流期水情预报方法可参见第9.3节的枯水预报。复习思考题1. 施工区河道流量演算中的马斯京根参数x值与河道天然情况相同，不必以实测资料分析得到。F2. 施工区河道流量

31、演算中的马斯京根参数k值与河道天然情况相同，可以不必重新分析得到。F3. 施工截流期河道流量预报，可以采用天然河道情况下枯水预报模型。 T第9.5节 水文实时预报方法内容提要1. 水文预报误差的来源；2. 实时预报校正方法；3. 递推最小二乘估计算法。学习要求了解水文预报误差的来源及最小二乘估计的递推算法。水文实时预报的意义水文现象受到自然界中众多因素的影响,人们采用的各种方法或模型都不可能将复杂的水文现象模拟得十分确切,水文预报估计值与实际出现值的偏离,即预报误差是不可避免的。实时预报就是利用在作业预报过程中，不断得到的预报误差信息,及时地校正、改善预报估计值或水文预报模型中的参数,使以后阶

32、段的预报误差尽可能减小。 水文实时预报方法分类1.水文预报误差的来源水文预报误差的来源大致有如下几个方面： 模型结构误差在对水文循环过程的模拟中，采用了不同程度简化的模型或不完善的处理方法，由此引起的误差称为模型结构误差。模型参数估计误差水文模型中估计的模型参数对其真值来讲,总是存在着误差的。根据各场洪水优选的模型参数,它是综合各场洪水的最优值,而对某一特定场次的洪水,它并非就是最合适的。模型的输入误差进行水文预报所输入的资料通常是降雨、流量和流域蒸散发,这些资料或由实测获得,或根据天气预报估算得到。前者存在着测验和时段统计平均误差,后者则存在着相当大的预报误差。2. 实时预报校正模型洪水实时

33、预报校正方法包含两方面的内容：一是实时预报校正模型，二是实时校正方法。实时预报校正模型在很大程度上取决于水文模型的结构。“显式”结构的水文模型当水文模型相对于模型的待定参数是线性关系模型时，称为为“显式”结构。一种处理方法是将水文系统视为动态系统，模型的动态参数“在线”识别和实时预报是关键。另一种处理方法是将水文预报模型改造成系统状态方程和系统观测方程，利用滤波的方法进行实时校正。“隐式”结构水文模型一般来讲，流域概念性水文模型是较复杂的“隐式”结构。目前在处理这类模型时，一种方法是对模型进行“显式”化处理；另一种方法是基于确定性流域水文模型的预报流量与实测流量的误差序列，建立流量误差预报的“

34、显式”模型（如AR或ARMA模型），流域洪水预报即用预报的流量残差叠加到模型的计算流量上，从而完成洪水实时校正预报。3.实时预报校正方法的分类根据不同预报误差的来源,实时预报校正方法可分为以下三种。对模型参数实时校正若认为水文预报方法或模型的结构是有效的,只是由于存在数据的观测误差,导致率定的模型参数不准确,或是率定的模型参数对具体场次洪水并非最优，可以在实际作业预报过程中,根据实际的预报误差不断地修正模型参数。对模型参数进行实时校正的方法有最小二乘估计等方法。对模型预报误差进行预测 对已出现的预报误差时序过程,建立合适的预报误差的模型。通过预报未来的误差值以校正尚未出现的预报值,从而提高水文

35、预报的精度。对状态变量进行估计一个预报模型中能控制当前及以后时刻系统状态和行为的变量,称为状态变量。对状态变量的估计是认为预报误差来源于状态估计的偏差和实际观测的误差,通过实时修正状态变量来提高预报的精度。卡尔曼滤波就是对状态变量进行实时校正的一种算法。 水文实时预报的最小二乘方法在水文预报模型参数的估计中，最小二乘法是一种常用的估计方法。通过最小二乘估计可以获得一个在最小方差意义上与实测数据拟合最好的模型。1.最小二乘估计的基本算法设y(i)和为在i时刻（i=1,2,，m）所观测得的数据，可以用m个方程表示这些数据的关系（9-5-1）式中y为系统观测输出值，X为系统观测输入值，V为误差向量，是一个待定的参数向量。上式用矩阵形式表示为（9-5-2）式中： （9-5-3）在水文模型中，经常遇到测量次数m超过方程组所需的定解条件数n。最小二乘原理指出，最可信赖的参数值应在使残余误差平方和最小的条件下求出。 目标函数为： （9-5-4）将J对求偏微分，并令其等于零，则可求得使J趋于最小的估计值，有（9-5-5）2. 序贯递推最小二乘算法它是基于最小二乘推导出的、利用新息来改进参数的递推估计算法，使参数实现在线识别。参数向量的递推算式： （9-5-6） （9-5-7） （9-5-8） 在递推过程中，当计及数据量不多时，新