洪洪水水预预报报与与调调度度 2012.10.22 2012.10.22授课人：俞授课人：俞 淞淞uu 姓名：俞淞姓名：俞淞uu 教育背景：教育背景：uu 籍贯：安徽省滁州市定远县籍贯：安徽省滁州市定远县 自我介绍自我介绍1999.92003.7 1999.92003.7 武汉大学水文与水资源工程专业武汉大学水文与水资源工程专业 本科本科2003.92005.7 2003.92005.7 武汉大学水文学及水资源专业武汉大学水文学及水资源专业 硕士硕士2005.92008.7 2005.92008.7 武汉大学水文学及水资源专业武汉大学水文学及水资源专业 博士（直博）博士（直博）uu 出生年月：出生年月：19821982年年3 3月月樱花烂漫樱花时节春意盎然珞珈夏韵层林尽染樱园雪景五彩缤纷日落时分老图的天空珞珈夜色倦鸟归巢雨后樱花路 樱花烂漫 Email:yusongbnu.edu.cnQQ：85283180联系方式联系方式本本次次课课程程及及应应用用实实例例洪洪洪洪水水水水预预预预报报报报实实实实时时时时校校校校正正正正方方方方法法法法介介介介绍绍绍绍：残残残残差差差差回回回回归归归归、自自自自适适适适应反馈、卡尔曼滤波等方法应反馈、卡尔曼滤波等方法应反馈、卡尔曼滤波等方法应反馈、卡尔曼滤波等方法水文预报，流域产汇流计算水文预报，流域产汇流计算水文预报，流域产汇流计算水文预报，流域产汇流计算主主 要要 内内 容容水水水水库库库库洪洪洪洪水水水水调调调调节节节节计计计计算算算算：主主主主要要要要介介介介绍绍绍绍水水水水库库库库的的的的调调调调洪洪洪洪作作作作用用用用，设计标准，调洪计算基本原理和方法。设计标准，调洪计算基本原理和方法。设计标准，调洪计算基本原理和方法。设计标准，调洪计算基本原理和方法。水水水水库库库库汛汛汛汛期期期期分分分分期期期期方方方方法法法法介介介介绍绍绍绍：数数数数理理理理统统统统计计计计、变变变变点点点点分分分分析析析析、模糊集合分析、矢量统计等方法的比较及应用。模糊集合分析、矢量统计等方法的比较及应用。模糊集合分析、矢量统计等方法的比较及应用。模糊集合分析、矢量统计等方法的比较及应用。引引言言u水孕育了生命，然而从各种自然灾害发生的次数来看，水孕育了生命，然而从各种自然灾害发生的次数来看，每年有每年有6060以上以上是是洪水灾害洪水灾害，特别是亚洲、北美发生的，特别是亚洲、北美发生的次数更多。次数更多。u我国我国是世界上受洪涝灾害影响最大的国家之一，频繁的是世界上受洪涝灾害影响最大的国家之一，频繁的洪水灾害每年都给社会经济和人民的生命财产造成洪水灾害每年都给社会经济和人民的生命财产造成巨大巨大损失损失。据历史记载，从公元前。据历史记载，从公元前206206年到年到19491949年中华人民共年中华人民共和国成立，共和国成立，共21552155年间，中国就发生过较大的水灾年间，中国就发生过较大的水灾10921092次，较大的旱灾次，较大的旱灾10561056次，平均次，平均几乎每年几乎每年都有较大水旱灾都有较大水旱灾害。害。 u 研究的主要意义研究的主要意义 随着社会经济的高速发展，水库在防洪中随着社会经济的高速发展，水库在防洪中的地位和作用越来越的地位和作用越来越突出突出。水库不仅要确保。水库不仅要确保自身的自身的安全安全，还要担负发电等多种，还要担负发电等多种综合利用综合利用任务，任务，洪水预报和调度洪水预报和调度任务越来越任务越来越艰巨而复艰巨而复杂杂。 概概 述述 （一）洪水预报与调度在水库防洪中的作用（一）洪水预报与调度在水库防洪中的作用水库特征水位与特征库容：死水位Z死与死库容V死 正常蓄水位Z蓄与兴利库容V蓄 防洪限制水位Z限与结合库容V结防洪高水位Z防与分洪库容V防设计洪水位Z设与拦洪库容V拦校核洪水位Z校与调洪库容V调 水库特征水位及其相应库容示意图 Z死、Z蓄、Z防、Z设、 Z校建设时已经确定， Z限设计时一般对于正常蓄水位Z蓄，即 Z限=Z蓄。由于后来的工程质量、水文情况变化等，常使Z限Z蓄 1 1保证自身安全保证自身安全 设计洪水来临时，水库Z= Z限，最高洪水位Z设 校核洪水来临时，水库Z= Z限，最高洪水位Z校 考虑预报时考虑预报时 预泄：在保证兴利的条件下，适当腾空部分兴利库容，提高水库的抗洪标准 超蓄：保证安全的条件下，在洪水尾部适当超蓄临时多拦蓄部分洪水（超蓄的洪水能在预见期内腾空），提高水库的兴利作用2. 2. 保证下游安全保证下游安全 下游洪水防洪标准时，水库ZZ防,下游防洪点QQ安 考虑预报时 通过防洪调度，与下游防洪点错峰，提高下游的防洪标准 （二）洪水复核与汛限水位设计（二）洪水复核与汛限水位设计洪水复核：洪水复核：新情况下的设计洪水、校核洪水、调洪规则洪水分期洪水分期工程任务与工程质量（三）（三）防洪预案研究防洪预案研究对于不同标准的洪水对于不同情况的实际特大洪水（四）（四）基本步骤基本步骤洪水预报 离线预报：流域降雨径流预报 河段洪水预报 在线预报实时洪水预报：最小二乘方法 卡尔曼滤波方法水库调度流流域域降降雨雨径径流流预预报报水文循环示意图水文循环示意图水文循环示意图水文循环示意图 流域降雨径流预报流域降雨径流预报一、产流计算一、产流计算 净雨过程计算净雨过程计算 一个流域的Wm、b一定，可按上式由降雨、蒸发连续预报W 、R R进一步划分为二水源：地面净雨RS、地下净雨Rg三水源：地面净雨、壤中流净雨、地下净雨二、汇流计算二、汇流计算 1 1地面汇流（或河网汇流）地面汇流（或河网汇流）时段单位线（经验单位线）瞬时单位线J.E.Nash单位线 纳希瞬时单位线仅有两个参数n和K，便于进行单站和地区综合，故在洪水预报和设计洪水计算中应用较为普遍2 2地下汇流（包括攘中流净雨汇流）地下汇流（包括攘中流净雨汇流） 按地下水库进行演算预报方案评定预报方案评定 我国颁布的SD1382000水文情报预报规范采用以下误差评定方法 1 1评定标准评定标准（1）确定性系数dy（2）合格率 m预报值误差许可误差的场数； n预报的总场数，一般不少于40次 许可误差：对于降雨径流预报，净雨深的许可误差采用实测值的20%，许可误差大于20mm时，以20mm为上限；许可误差小于3mm时，以3mm为下限。 洪峰流量的许可误差取实测值的20%，并以流量测验误差为下限。 洪峰流量出现时间的许可误差，取预报根据时间至实测峰现时间间距的30%，并以3h或一个计算时段为下限。 a)确定性系数 b)洪峰合格率 c)峰现时差合格率 2 2预报方案的评定标准预报方案的评定标准（1）按确定性系数评定有效性 评定（检验）方案的有效性时按下列标准进行。（2）按许可误差评定合格率 预报方案合格率是评定（检验）中计算值与实测值之差不超过许可误差的次数（m）占全部次数（n）的百分率（m/n100）。按其合格率将预报方案划分为以下3个等级。 预报方案经评定达到上述甲、乙两个等级者，即可用于作业预报；达到丙等的方案用于参考性预报；丙等以下的方案不能用于作业预报，只能作参考性估报。3.3.作业预报的评定作业预报的评定三、河段洪水预报三、河段洪水预报河槽蓄水量与泄流量关系方程蓄泄方程 W=f（Q上，Q下）=f(Q) (3-2) 马斯京根法将该方程表达为Q上，Q下的线性组合联解式（3-1）、（3-2），即可求得Q下t2.马斯京根法 (1)马斯京根流量演算方程该法将槽蓄量S分为二部分柱蓄和锲蓄组成：柱蓄量=河段通过稳定流量Q下的槽蓄量SQ下=KQ下， K为稳定流情况下的河段流量传播时间锲蓄量=x(SQ上-SQ下)=x (KQ上- KQ下) x楔蓄系数，或称流量比重因数则槽蓄量 S=SQ下+ x(SQ上-SQ下)= KQ下+ x (KQ上- KQ下) 亦 S=Kx Q上+(1- x) Q下 =KQ (3-3) Q= x Q上+(1- x) Q下= Q下+ x(Q上-Q下) (3-4) Q示储流量（表示槽蓄量S相应的稳定流流量）式（3-1）式、（3-3）联解，即得马斯京根流量演算方程 （3-5）(2) K和x的确定 马斯京根蓄泄方程的S和Q呈线性关系，表明Q是S下的稳定流流量，所以，可假设x，作S和Q关系线，当为直线时，如图，x即为所求，此时的直线斜率即为K值。(3) 马斯京根法的几个问题与处理方法（1）K值的综合分析多次洪水的K值，如果变化不大，可取平均值如果K值变化较大，一般可按流量分级，应用时，随预报流量的变化取相应的K值（2）x值的综合 x值一般变化不大，约0.20.45，可取多次洪水分析的平均值（3） 的选取 一般 ，最好取等于河段的洪水传播时间K。如果难以满足，则可将整个河段按n=K/分为n个河段，采用分段连续演算法计算（4）马斯京根法本身无预见期，一般利用降雨径流法预报出上断面Q上t，再用马法计算Q下t。但取=2Kx时，C0=0，则有 的预见期洪洪水水预预报报实实时时校校正正方法方法洪水预报实时校正方法洪水预报实时校正方法水文预报是人类对自然界水文现象的模拟和预测，影响因素水文预报是人类对自然界水文现象的模拟和预测，影响因素众多。一般来讲，洪水预报方案是根据以往的实测资料编众多。一般来讲，洪水预报方案是根据以往的实测资料编制的，其中洪水预报模型的参数，反映的都是以往资料情制的，其中洪水预报模型的参数，反映的都是以往资料情况下况下平均而言平均而言的最优值，但在作业预报时，当实际情况偏的最优值，但在作业预报时，当实际情况偏离方案的状态时，就会使预报结果发生偏差。如果我们能离方案的状态时，就会使预报结果发生偏差。如果我们能够利用作业预报过程中随时得到的够利用作业预报过程中随时得到的误差信息（新息）误差信息（新息），恰，恰当的调整下一步用产汇流预报方案推算所得的数值，将会当的调整下一步用产汇流预报方案推算所得的数值，将会提高预报精度，这种利用新息进行预报校正的方法称为提高预报精度，这种利用新息进行预报校正的方法称为实实时校正时校正预报。预报。 洪水预报校正模型主要采用了两种类型，一是采洪水预报校正模型主要采用了两种类型，一是采用用相关分析相关分析校正模型，二是采用校正模型，二是采用自回归自回归校正模型。校正模型。采用的校正技术主要有，利用采用的校正技术主要有，利用最小二乘最小二乘分析技术分析技术建立的建立的定常参数定常参数校正模型方法；利用递推最小二校正模型方法；利用递推最小二乘分析技术建立的乘分析技术建立的自适应反馈自适应反馈洪水实时预报校正洪水实时预报校正模型方法和利用卡尔曼滤波技术建立的洪水流量模型方法和利用卡尔曼滤波技术建立的洪水流量实时校正方法，经比较后，选择一种适合水库情实时校正方法，经比较后，选择一种适合水库情况的校正模型和方法。况的校正模型和方法。模型预报流量与实测流量相关校正方法模型预报流量与实测流量相关校正方法选取湖北省漳河水库选取湖北省漳河水库 利用利用1990199020022002年汛期反推的入库洪水过程和产年汛期反推的入库洪水过程和产汇流预报模型模拟计算的流量过程资料，汇流预报模型模拟计算的流量过程资料，模型预报模型预报流量流量 与实时校正后的预报流量与实时校正后的预报流量 相关分析相关分析校正方程为：校正方程为： 其相关系数其相关系数 r = 0.91745 r = 0.91745 。 多阶自回归洪水实时校正方法多阶自回归洪水实时校正方法 概念性流域水文预报模型使得水文模型呈现为较复杂的概念性流域水文预报模型使得水文模型呈现为较复杂的“隐式隐式”结构，模型参数识别存在着计算的复杂性，很难用实结构，模型参数识别存在着计算的复杂性，很难用实测数据实时修正模型参数，成为定常参数的流域水文模型。测数据实时修正模型参数，成为定常参数的流域水文模型。目前在利用定常参数的流域水文模型进行洪水实时预报时，目前在利用定常参数的流域水文模型进行洪水实时预报时，采用的有效方法是用确定性流域水文模型加上实时校正处理采用的有效方法是用确定性流域水文模型加上实时校正处理算法进行洪水实时校正预报。即利用流域水文模型的算法进行洪水实时校正预报。即利用流域水文模型的预报流预报流量序列量序列与与实时观测的流量序列实时观测的流量序列的的残差序列残差序列dQ(j),j=1,.,tdQ(j),j=1,.,t，建立残差预报模型，将预报的残差，建立残差预报模型，将预报的残差dQ(j),j=t+1,.dQ(j),j=t+1,.叠加叠加到预报流量上，完成流域洪水预报校正，从而提高洪水预报到预报流量上，完成流域洪水预报校正，从而提高洪水预报精度。精度。 对于漳河水库对于漳河水库 一阶残差自回归洪水预报校正模型为：一阶残差自回归洪水预报校正模型为： 二阶残差自回归洪水预报校正模型为：二阶残差自回归洪水预报校正模型为：自适应反馈洪水实时校正方法自适应反馈洪水实时校正方法 1. 1. 洪水实时预报校正方法与技术洪水实时预报校正方法与技术 在在洪洪水水实实时时预预报报中中，目目前前阶阶段段采采用用的的自自适适应应反反馈馈洪洪水水信信息息实实时时校校正正的的方方法法主主要要有有：递递推推最最小小二二乘乘法法和和卡卡尔尔曼曼滤滤波波方方法法。这这些些实实时时校校正正的的方方法法共共同同特特点点是是，能能实实时时地地处处理理水水文文系系统统最最新新出出现现的的预预报报误误差差，以以此此作作为为修修正正预预报报模模型型参参数数、状状态态或或预预报报输输出出值值的的依据，从而使预报系统迅速适应现时状况。依据，从而使预报系统迅速适应现时状况。 2. 2. 自适应洪水实时校正模型自适应洪水实时校正模型 将将一一阶阶自自回回归归洪洪水水预预报报校校正正模模型型改改写写为为时时变变参数的模型，其结构形式为：参数的模型，其结构形式为： 将将上上述述时时变变参参数数的的模模型型，结结合合改改进进的的可可变变遗遗忘忘因因子子最最小小二二乘乘参参数数估估计计自自适适应应递递推推算算法法，从从而而得得到到一阶自适应洪水实时校正模型具体结构形式为：一阶自适应洪水实时校正模型具体结构形式为： （1 1）一阶自适应洪水实时校正模型递推算法为：）一阶自适应洪水实时校正模型递推算法为：预报：预报：参数估计：参数估计： 增益因子：增益因子： 协方差阵：协方差阵： 可变遗忘因子：可变遗忘因子： 上述表达式中上述表达式中: : （2 2）二阶自适应洪水实时校正模型递推算法为：）二阶自适应洪水实时校正模型递推算法为：预报：预报：参数估计：参数估计： 增益因子：增益因子： 协方差阵：协方差阵： 可变遗忘因子：可变遗忘因子： 上述表达式中上述表达式中: : 卡尔曼滤波洪水实时校正方法卡尔曼滤波洪水实时校正方法卡尔曼滤波是卡尔曼和布西在卡尔曼滤波是卡尔曼和布西在6060年代初提出的，它是具有完年代初提出的，它是具有完整的严密理论体系的实时预报方法，在自动控制等领域应用整的严密理论体系的实时预报方法，在自动控制等领域应用取得很好的效果。自从取得很好的效果。自从7070年代被引入到洪水实时预报以来，年代被引入到洪水实时预报以来，人们进行了许多探索性的研究工作。人们进行了许多探索性的研究工作。卡尔曼滤波是根据对水文系统建立的状态方程和观测方程，卡尔曼滤波是根据对水文系统建立的状态方程和观测方程，采用线性递推的算法进行实时预报的。即采用线性递推的算法进行实时预报的。即由现时段的预报值由现时段的预报值和获得的观测值作滤波计算，然后预报下一时段的系统状态，和获得的观测值作滤波计算，然后预报下一时段的系统状态，待下一时段的观测值出现后，再滤波，再预报下一个时段的待下一时段的观测值出现后，再滤波，再预报下一个时段的状态，如此连续循环和预报状态，如此连续循环和预报。 卡尔曼滤波洪水实时预报校正模型卡尔曼滤波洪水实时预报校正模型 本课中的卡尔曼滤波洪水实时预报研究，是以自回归残差本课中的卡尔曼滤波洪水实时预报研究，是以自回归残差预报模型为基础，进一步导出卡尔曼滤波洪水实时预报模型预报模型为基础，进一步导出卡尔曼滤波洪水实时预报模型和计算方法。和计算方法。 其中自回归残差预报模型，可表示为：其中自回归残差预报模型，可表示为： 将自回归残差预报模型，改写成状态方程和观测方程的表将自回归残差预报模型，改写成状态方程和观测方程的表达形式，有达形式，有 （1 1）状态方程）状态方程可进一步简写为：可进一步简写为： 其中状态变量、噪声分配矩阵、状态转移矩阵、噪声方差其中状态变量、噪声分配矩阵、状态转移矩阵、噪声方差矩阵为：矩阵为：，（2）观测方程）观测方程在在t+1时刻流量是可观测的，此时可推求得到流量预报误差时刻流量是可观测的，此时可推求得到流量预报误差，假设测量的流量值受到误差影响，从而有：，假设测量的流量值受到误差影响，从而有：可简写为：可简写为：式中，式中，的方差取的方差取R=147。实时校正方法比较与选用实时校正方法比较与选用漳河水库洪水预报校正模型及算法结果综合表漳河水库洪水预报校正模型及算法结果综合表模型方法类别确定性系数(%)洪峰合格率(%)峰现时差合格率(%)产汇流预报模型80.1283.3383.33流量相关校正法84.1741.6783.33一阶残差自回归校正模型89.6310083.33二阶残差自回归校正模型89.6710083.33一阶自适应实时校正模型89.6410083.33二阶自适应实时校正模型85.5683.3391.67卡尔曼滤波实时校正算法89.5710083.33水水库库调调度度电网调度管理条例实施办法电网调度管理条例实施办法（1994年10月11日中华人民共和国电力工业部令第3号发布自发布之日起施行）第四条电网调度机构一般应当进行下列主要工作：（一）、组织编制和执行电网的调度计划（运行方式）；（二）、指挥调度管辖范围内的设备的操作；（三）、指挥电网的频率调整和电压调整；（四）、指挥电网事故的处理，负责进行电网事故分析，制订并组织实施提高电网安全运行水平的措施；（五）、编制调度管辖范围内的设备的检修进度表，批准其按计划进行检修；（六）、负责本调度机构管辖的继电保护和安全自动装置以及电力通信和电网调度自动化设备的运行管理；负责对下级调度机构管辖的上述设备和装置的配置和运行进行技术指导；（七）、组织电力通信和电网调度自动化规划的编制工作，组织继电保护及安全自动装置规划的编制工作；（八）、参与电网规划编制工作，参与电网工程设计审查工作；（九）、参加编制发电、供电计划，监督发电、供电计划执行情况，严格控制按计划指标发电、用电；（十）、负责指挥全电网的经济运行；（十一）、组织调度系统有关人员的业务培训；（十二）、（十二）、统一协调水电厂水库的合理运用；统一协调水电厂水库的合理运用；（十三）、协调有关所辖电网运行的其他关系。第十八条调度机构应当按年、月、日编制并下达发电调度计划。凡由调度机构统一调度并纳入电网进行电力、电量平衡的发电设备，不论其产权归属和管理形式，均必须纳入发电调度计划的范围。月度发电调度计划，须在年度发电预期计划的基础上，综合考虑用电负荷需求、月度水情、燃料供应、核燃料的燃耗、供热机组供热等情况和电网设备能力、设备检修情况等因素进行编制。日发电调度计划在月发电调度计划的基础上，综合考虑日用电负荷需求、近期内水情、燃料供应情况和电网设备能力、设备检修情况等因素后，编制出日发电（有功、无功功率或者电压）曲线，下达各发电厂执行。第二十二条跨省电网管理部门和省级电网管理部门编制发电、供电计划以及调度机构编制发电、供电调度计划时，对具有综合效益的水电厂（站）的水库，不论其产权归属和管理形式，均应当根据批准的水电厂（站）的设计文件，合理运用水库蓄水，保证其正常运用，不允许发生水库长期处于降低出力区运行的情况。多年调节水库在蓄至正常蓄水位后，供水期末水位一般应当控制在年消落水位附近。遭遇连续枯水年时，需降至死水位运行，必须按规定报批。年调节水库一般在每年汛末应当蓄至正常蓄水位。多年调节水库和年调节水库的最低运行水位，均不允许低于死水位，不得破坏水库的调节能力。调节性能差的水库一般也不能低于死水位运行。第二十四条跨省电网管理部门和省级电网管理部门遇有下列情形之一，而且不能在短期内解决，需要调整年、月度发电、供电计划指标时，可以适当调整。调整后，应当立即通报有关地方人民政府的有关部门，并与上述有关部门共同研究协调，采取有效措施：、大、中型水电厂（站）水库实际来水与编制发电计划依据的来水预计相差较大；、火电厂的燃料库存低于规定的火电厂最低燃料库存量警戒线；、其他需调整发电、供电计划的情形。水调自动化系统概述水调自动化系统概述一、建立水调自动化系统的目的与意义一、建立水调自动化系统的目的与意义 社会主义市场经济体制的健康发展，加速了电力系统的建设步伐社会主义市场经济体制的健康发展，加速了电力系统的建设步伐, ,我国电力工业正在形成投资、产权、利益主体多元化的格局，如何建我国电力工业正在形成投资、产权、利益主体多元化的格局，如何建立高效的电力市场机制，立高效的电力市场机制，优化优化网内资源配置、网内资源配置、降低电网运行成本降低电网运行成本、保、保证电力系统经济可靠性、提高水电站水库调度自动化程度，是各水电证电力系统经济可靠性、提高水电站水库调度自动化程度，是各水电厂、网局乃至国家电力公司普遍关注和亟需解决的重大课题。因此，厂、网局乃至国家电力公司普遍关注和亟需解决的重大课题。因此，实现水电站水库及水电系统调度高度自动化，开发实用型高级应用软实现水电站水库及水电系统调度高度自动化，开发实用型高级应用软件，建立专家决策支持系统，确保电厂、电网安全、稳定、优质、经件，建立专家决策支持系统，确保电厂、电网安全、稳定、优质、经济运行，达到济运行，达到无人或少人值班无人或少人值班的高度集成自动化调度，对于提高科学的高度集成自动化调度，对于提高科学管理水平、提高全网水电厂的水能利用率，增加发电效益，获取整体管理水平、提高全网水电厂的水能利用率，增加发电效益，获取整体最佳的国民经济效益，促进工农业发展，为市场经济条件下最佳的国民经济效益，促进工农业发展，为市场经济条件下水库（群）水库（群）优化调度提供切实可行的途径，无疑具有重大的科技意义和广阔的应优化调度提供切实可行的途径，无疑具有重大的科技意义和广阔的应用前景用前景。 水电水电作为一种洁净的作为一种洁净的可再生能源可再生能源，具有，具有运行费用低廉运行费用低廉、机组启停机组启停灵活、调峰能力强灵活、调峰能力强的优势，早已成为国际国内首选开发的能源，正因的优势，早已成为国际国内首选开发的能源，正因为如此，水电在整个电网中的作用与地位越来越重要和突出。可见尽为如此，水电在整个电网中的作用与地位越来越重要和突出。可见尽快实现水电系统调度自动化及专家决策支持系统势在必行，其社会和快实现水电系统调度自动化及专家决策支持系统势在必行，其社会和经济效益不可低估。经济效益不可低估。二、水调自动化系统现状二、水调自动化系统现状 水水电电系系统统调调度度自自动动化化及及专专家家决决策策支支持持系系统统是是一一个个集集水水库库调调度度专专业业知知识识，自自动动化化硬硬件件设设备备与与接接口口、计计算算机机及及网网络络通通讯讯技技术术、决决策策支支持持理理论论为为一一体体的的大大规规模模集集成成系系统统，其其研研究究需需要要花花费费大大量量的的人人力力、物物力力和和财财力力。国国外外由由于于自自动动化化技技术术和和计计算算机机通通信信技技术术水水平平较较高高, ,其其水水调调自自动动化化系系统统研研究究起起步步早早, ,实实用用化化程程度度高高。从从我我们们收收集集到到国国外外1010余余个个水水调调自自动动化化系系统统研研究究报报告告看看, ,主主要要依依据据数数据据库库、模模型型库库和和计计算算机机交交互互式式图图形形界界面面来来确确定定水水库库长长、中中、短短期期运运行行计计划划。如如美美国国加加利利福福尼尼亚亚中中央央河河谷谷工工程程（C C）长长、中中、短短期期优优化化调调度度、美美国国田田纳纳西西流流域域（）的的水水资资源源系系统统周周和和日日的的调调度度，加加拿拿大大马马尼尼托托巴巴大大学学开开发发的的水水库库调调度度决决策策支支持持系系统统等等。但但是是由由于于体体制制、管管理理模模式式等等差差异异，以以及及实实际际客客观观条条件件和和经经费费等等限限制制，近近年年来来国国内内一一些些水水电电站站引引进进的的水水调调自自动动化化系系统统主主要要集集中中在在硬硬件件，引引进进的的软软件件因因汉汉化化等等原原因因使使用用的的效效果果并并不不理理想想。显显然然，水水调调自自动动化化系系统统的的核核心心-应应用用软软件件的的开开发发研研制制，必必须须在在汲汲取取世世界界先先进进经经验验的的基基础础上上立立足于国内科技力量开发足于国内科技力量开发。 国内，自年代末开展水库优化调度以来，在理论和应用方面都取国内，自年代末开展水库优化调度以来，在理论和应用方面都取得了长足发展，尤其是应用确定性和随机性优化模型编制水电站水库（群）得了长足发展，尤其是应用确定性和随机性优化模型编制水电站水库（群）中长期调度方案，结合水文预报开展短期预报调度（防洪及发电）和实时中长期调度方案，结合水文预报开展短期预报调度（防洪及发电）和实时调度，创造了巨大的经济效益（例如柘溪水库优化调度、白山与丰满水库调度，创造了巨大的经济效益（例如柘溪水库优化调度、白山与丰满水库联合优化调度、闽江流域梯级水电站联合经济调度等）。大系统理论、多联合优化调度、闽江流域梯级水电站联合经济调度等）。大系统理论、多目标理论及神经网络理论等新的理论方法也不断用来研究和解决水电站水目标理论及神经网络理论等新的理论方法也不断用来研究和解决水电站水库（群）优化调度问题，但主要是模型构造和求解方法的研究，成功的、库（群）优化调度问题，但主要是模型构造和求解方法的研究，成功的、成熟的、实用的库群优化调度软件较少。成熟的、实用的库群优化调度软件较少。“八五八五”期间，以北京水利水电期间，以北京水利水电科学研究院自动化研究所和南京自动化研究院为主研制的科学研究院自动化研究所和南京自动化研究院为主研制的水情自动测报水情自动测报系系统取得了可喜的成就，在全国电力行业已建和在建的水情自动测报系统达统取得了可喜的成就，在全国电力行业已建和在建的水情自动测报系统达6060多个；在此期间，多个；在此期间，水调自动化工作水调自动化工作也取得了一定的进展，一些网局，省也取得了一定的进展，一些网局，省局初步实现了水调信息的自动采集和监视功能，局初步实现了水调信息的自动采集和监视功能，19961996年由武汉水利电力大年由武汉水利电力大学和东北电力集团公司联合进行了东北电网水库调度自动化系统总体规划，学和东北电力集团公司联合进行了东北电网水库调度自动化系统总体规划，19971997年由武汉水利电力大学和中国水科院联合研制的西北电网水调自动化年由武汉水利电力大学和中国水科院联合研制的西北电网水调自动化系统，已投入试运行。系统，已投入试运行。 19991999年东北电力集团公司与武汉水利电力大学签年东北电力集团公司与武汉水利电力大学签署协议，联合开展东北电网水库调度自动化系统软件开发工作。署协议，联合开展东北电网水库调度自动化系统软件开发工作。 河海大学、大连理工大学、河海大学、大连理工大学、清华大学、华中理工大学清华大学、华中理工大学等在水电站综合等在水电站综合自动化系统及有关模型理论方法等方面也取得了可喜的成果。自动化系统及有关模型理论方法等方面也取得了可喜的成果。水库防洪调度水库防洪调度 （一）（一） 任务任务 根据调洪规则由入库洪水过程推求出库洪水过程、库水位变化过程、最高库水位等 （二）（二） 水库调洪计算原理水库调洪计算原理 1泄洪建筑物的泄流能力闸门全开情况 q=f(Z) 或 q=f(V)溢洪道堰流公式泄洪洞孔口出流公式（二）（二）水库调洪计算基本原理水库调洪计算基本原理根据初始、边界条件，逐时段联介式水库水量平衡方程(14-1)和蓄泄方程(14-2):（14-1）式中-分别为计算时段始,末入库流量-分别为计算时段始,末出库流量-分别为计算时段始,末的水库蓄水量-计算时段,洪水陡涨陡落时,取短些. （14-2） 通过逐时段联解,就可以求得水库下泄流量过程线,最大下泄流量,这场洪水所需的防洪库容和水库最高洪水位.如图14-3所示. 3. 调洪计算（四）水库防洪调度1.控制泄流与自由泄流相结合调度无下游防洪任务，入库流量Q泄流能力 q，使q=Q 入库流量Q泄流能力 q，敞开泄流2.固定泄流或分级控制泄流调度下游有防洪任务，按下游安全泄量Q安的情况，或固定泄流，或分若干级控制泄流3.错峰调度 根据下游防洪点的安全泄量Q安的要求，应Qc,t-tca= Q安t - QB,t-tBA - QABC,t-t区 或 Qc,t= Q安t+tca - QB,t-tBA+ tca - QABC,t-t区+ tca = Q安 - QB,t+tD - QABC,t+tD区 4.防洪预报调度 水水库库汛汛期期分分期期方方法法研研究究研究的背景及意义研究的背景及意义uu我国是世界上严重缺水的国家之一我国是世界上严重缺水的国家之一，人均水资源占有量仅，人均水资源占有量仅为世界平均水平的为世界平均水平的1/41/4，是世界上，是世界上1313个贫水国之一。个贫水国之一。uu我国又是世界上用水量最多的国家之一我国又是世界上用水量最多的国家之一，据统计，仅，据统计，仅20022002年，全国淡水取用量达到年，全国淡水取用量达到54975497亿亿m m3 3，约占世界年取用，约占世界年取用量的量的13%13%。我国从。我国从2020世纪世纪7070年代以来就开始出现水资源供年代以来就开始出现水资源供需矛盾，并由局部逐渐蔓延至全国大部分地区。需矛盾，并由局部逐渐蔓延至全国大部分地区。uu随着人口增长和社会经济的发展，污水排放量也在不断增随着人口增长和社会经济的发展，污水排放量也在不断增加，水质恶化加剧了水资源短缺矛盾。加，水质恶化加剧了水资源短缺矛盾。已成为国民经济已成为国民经济可持续发展的制约因素。可持续发展的制约因素。 静静静静态态态态控控控控制制制制法法法法设计理念设计理念动动动动态态态态控控控控制制制制法法法法汛期分期汛期分期固固固固定定定定汛汛汛汛限限限限水水水水位位位位分分分分期期期期汛汛汛汛限限限限水水水水位位位位气候成因分析法气候成因分析法：对流域气候分析。经验统计法经验统计法：利用多年各月、旬或日的实测雨 量、洪峰、洪量等资料，选择统计指标，通过数理统计方法得出汛期的变化规律，从而确定汛期的划分点。有多年降雨；多年径流；超定量洪峰；年最大旬径流分析等。模糊分析法模糊分析法：经验隶属度法和理论隶属度法。汛期分期研究汛期分期研究结论：结论：结合以上各种方法作分析比较，前汛期46月，5月和6月为主汛期，4月为前期过渡期（归为前汛期）；后汛期79月，7月为后期过渡期；相邻月份间有明显变化。 u一般地说，变点就是“模型中的某个或某些量起突然变化之点”。这种突然变化往往反映事物的某种质的变化，在自然界、社会及各种领域中很常见且具有重要性。变点分析方法采用如下符号表示整个水文序列和前、后两部分。采用如下符号表示整个水文序列和前、后两部分。变点分析的任务就是检验变点最有可能发生的时间位置和所研究的统计特征值的变化幅度。 均值变点分析方法均值变点分析方法 逐逐月分期及其汛限水位计算由报告中的调洪规则一和规则二分别调洪演算，其中7月和8月为后汛期，其调洪规则与前汛期采用相同调洪规则，更趋于安全。调洪规则二控制条件相对于调洪规则一较少，则计算的最高水位应高于调洪规则一计算的最高水位，但由于规则二调洪演算中可能比规则一提前敞泄，部分最高水位会低于规则一最高水位。逐日滑动汛限水位逐日滑动动态汛限水位u在多年资料基础上，通过统计汛期中每日为起点的最大来水量，求出固定频率的设计来水量，选择合适的典型洪水过程线进行调洪计算后得到每天的汛限水位 优点：有效地节约了水资源不足：工作人员每天控制一个水位，操作过程中比较麻烦 水库汛期分期方法比较u统计旬月一日或三日最大值序列，采用峰量综合控制法经过调洪计算求得相应的各旬和各月的汛限水位优点：一定程度上考虑了季节变化的物理成因不足：较逐日法不够细 旬月划分 u根据防洪安全和兴利蓄水要求、河道及水库工程状态，利用水文气象和统计规律分析流域“汛期”或“主汛期”的确切涵义优点：结果比较合理，有较高的可靠性而被普遍采用不足：需要对成灾天气进行大量的分析，尤其对于大流域，成灾天气有很多组合方式，因而工作量大，同时分期有一定的主观性，也难于将汛期分到较细的时段（如日） 成因分析u利用实测历史流量（雨量）资料，选择统计指标，分析指标在年内（或汛期）的变化规律，最后通过数理统计理论得出汛期的变化规律优点：简单实用不足：在分析过程中，对是不是洪水这一临界值带有一定的主观性；同时也很难将汛期划分得较细 数理统计法 u采用模糊集理论，使用融成因分析、数理统计、模糊统计为一体的模糊集合综合分析法对水库的汛期划分进行研究优点：考虑了汛期在时间上的模糊性，在理论上有了较大的发展，具有先进性不足：分析结果对所选用的指标阈值比较敏感，而指标阈值在取值上任意性较大：同时该方法主要用于主汛期划分，对于有多个分期要求则尚嫌不足 模糊分析法 u洪峰散点序列在一定尺度范围内表现出了自相似性，认为洪峰散点序列是一种分形，用分维作指标划分洪水分期优点：比较客观，受经验和人为影响较小不足：分析计算的工作量较大。如何判断对象是分形，至今仍是一个需要进一步讨论的问题 分形分析法 u基于统计理论，用于检测时间序列突变，同时可以进行假设检验的划分时间序列的方法优点：结论合理，划分汛期可精确到日，并且在一定程度上更为客观、可靠不足：调度缺乏连续性 变点分析法 u选用描述流域降雨洪水等特性的多个因子，计算论域元素之间的相似系数以构成模糊相似矩阵，通过系统聚类分析进行汛期分期优点：避免了采用单因子进行分析带来的片面性，更具合理性不足：该方法对于汛期分期这样具有时序分布的样本聚类问题，仍具有一定的局限性，可能会产生将非连时序的时段归为一类的不合理现象 系统聚类法 u把每场洪水的发生日期看作一个矢量，根据各个矢量之间的方向相似性来判断分割点，即作为汛期分期点根据实际需要优点：比较直观，结论合理，划分汛期可精确到日不足：其应用具有一定的限制性，对相似矢量聚集的情况比较适用，而对于相同矢量累计的情况则分期效果不明显 矢量统计法 基于AM取样的矢量统计法分期 u 展望展望 l水文现象是一种水文现象是一种复杂的、不断变化复杂的、不断变化着的现象，它具有着的现象，它具有确定性与非确定性确定性与非确定性（不确定性）两个基本方面。我们（不确定性）两个基本方面。我们对水文规律的认识还很有限，对水文规律的认识还很有限，还有许多问题至今未为还有许多问题至今未为人们所认识和理解人们所认识和理解。l洪水预报结果总是存在或大或小的洪水预报结果总是存在或大或小的误差误差，采取一定的，采取一定的措施措施方法，提高精度，减少误差，这是水文工作者方法，提高精度，减少误差，这是水文工作者最基本的任务，也是洪水预报一直追求的目标。但是最基本的任务，也是洪水预报一直追求的目标。但是由于影响预报误差的因素非常多由于影响预报误差的因素非常多非常复杂，至今还非常复杂，至今还不可能不可能完全认识和采取尽善尽美的方法完全消除它。完全认识和采取尽善尽美的方法完全消除它。虽然做了一些工作，但不可能面面俱到，还有很多问虽然做了一些工作，但不可能面面俱到，还有很多问题有待解决。题有待解决。l人类活动对水文过程的影响越来越强烈，人类活动对水文过程的影响越来越强烈，生态环境生态环境日日趋恶化，且趋恶化，且气候变化气候变化异常，水文预报异常，水文预报复杂性和不确定复杂性和不确定性比较突出，给水文工作者带来了诸多挑战。如何做性比较突出，给水文工作者带来了诸多挑战。如何做好好水文中长期预报水文中长期预报；从；从可持续发展可持续发展的角度，使的角度，使水量与水量与环境保护相结合环境保护相结合，做好水质水量管理等方面还有很多，做好水质水量管理等方面还有很多工作要做。工作要做。科研之外科研之外u个人成长个人成长历史预测历史预测宇宙演化宇宙演化u随着研究的深入，更加赞叹自然的奇妙随着研究的深入，更加赞叹自然的奇妙和生命的可贵（身体本身就是一个复杂和生命的可贵（身体本身就是一个复杂系统）。珍惜现在，好好生活和观察身系统）。珍惜现在，好好生活和观察身边的一切。边的一切。u上善若水：水善利万物而不争上善若水：水善利万物而不争 谢谢！谢谢！