土耳其kirazdere水库运行调度.doc

土耳其KIRAZDERE水库运行调度 摘要 Kirazdere水库是投资8.6亿美元为Lzmit地区和西部土耳其环海湾地区供水项目的一个组成部分，水库由DSI，国家水力委员会规划设计。但是，后来的调查表明：要满足水库在高可靠性上达到额定输出量有一定难度。虽然后来考虑在溢洪道设计水位以上蓄水，通过安装用于洪水控制的大型弧形门来增加蓄水量，但是，本文主要限于描述最终目标是最大可能的获得额定输出量的长期战略的过渡期，包括在假定溢洪道设计水位决定最大可控制蓄水量的前提下，确定供水限制的时间和水量。  1.概述 最近投资建成的Kirazdere水库拥有268km2集水区，是Lzmit地区和西部土耳其环海湾地区的主要供水源，为约一百二十万人口和正迅速发展的工业园提供水资源。Kirazdere水库是投资额为8.6亿美元的BOT工程的一部分，该工程包括一个供水处理厂和长140km，日供水量480,000m3的引水管道系统。土坝建造在一个地震多发区，坝高108m，宽400m，在溢洪道设计水位情况下，蓄水量约为45,700,000km3。但是，如果在溢洪道以上安装四个大型弧形门，可以使最大蓄水量提高到大约60,000,000km3。 八十年代初期，水库由DSI，国家水力委员会设计，但是由于资金短缺，部分工程没有完成。在和泰晤士河委员会签订协议以后，将工程转变成为一个BOT（built operate transfer）项目，开始了对水库的重新评价，包括：即使在最小补偿流和最大蓄水量60,000,000m3的情况下，整个系统也将努力实现在高可靠性上达到年额定输出量1.42亿km3。所以，为了保证供水的可靠性，考虑采用弧形门来蓄水。当然，这一切必须是在保证不增加下游形成坝的筑堤被淹没的危险性的条件下完成。然后，将这些目标系统化，形成了一个长期控制计划，即：每年不同月份的最大蓄水量依据出现洪水的危险性进行调整。最后，应用最优化模型来决定进行供水限制的时间和程度。 与其相应的短期控制策略是指具体实施长期控制策略的方法，包括为洪水防治提供策略。这需要一个检测控制和数据采集（SCADA）网络系统，将水文资料和控制设置传输到控制中心。这里，利用实时降雨量作为预测水库入流的主要指标，依次输入最优化程序，来确定弧形门的控制设置。采用这种方法，就有可能利用降水和实际入流之间的时间差，从水库的溢洪道排水，保证水量不超过下游最大流量限制，从而达到更大的蓄水量。 为了发展综合最优控制策略，在它得以实施以前，可以假定最大蓄水量限制为45,7000,000m3，以确定供水限制的时间和水量。这个阶段被称为长期控制策略的过渡期。（目前，从SCADA网络系统得到的数据资料还很不充分，不足以可靠的估测短期预报程序中所需要的水文参数。）本文主要介绍系统化长期控制策略过渡期过程中所采纳的方法。 2. 优化模型 为了尽可能的达到额定输出值，最优化程序的目标是使实际供水量和需水量的差额的平方值最小。表达式如下： Min (1) 其中：T－总时间长度；Rt －在t时间范围内的实际需水释放量；Dt－在t时间范围内的需水量。 约束条件： St+1 = St + Qt – Rt –Et (2) Rmin ≤ Rt ≤ Rmax (3) Smin ≤ St ≤ Smax (4) S0 = ST (5) 其中：St －时间范围t内水库最初的蓄水量；St+1 －时间范围t内水库最终的水库蓄水量；Qt －时间范围t内水库的入流量；Et －时间范围t内净蒸发量；Smax －最大可调控蓄水量；Smin －最小可调控蓄水量或死区； 为了应用动态过程程序，将水库蓄水量设定为状态变量，供水量设定为决策变量，阶段的变化由时间来代表。因此，回归关系式可以表示为： (6) 其中：－目前阶段t到可调控极限范围内总供水不足额的最小值；n－从目前阶段到可调控极限范围内剩余的时间步长数； 通过动态过程程序求解方程（1）到（5）获得最优运行策略。 3．通过特定模式确定最优运行策略 通过求解确定性动力模型获得的最优运行策略取决于具体的入流条件，所以基于某一年的具体条件而得到的策略的应用范围是有限制的。实际应用中，解决入流不确定性的常用方法是将可调控问题系统转化为一阶马尔可夫链。这样，在初始条件下，可以转化为特定模式分析的方法。在特定模式分析的方法中，入流的不确定性通过一系列特定模式分类，每一年的入流量被认为是一个特定的模式。确定性的优化过程通过很多特定的模式来完成，生成不同优化解的集合，每一个优化解代表与特定模式相匹配的优化策略。基于这些优化的运行策略的例子，通过模式识别技术获得良好的运行规律。 4．采用人工神经网络派生的运行规律 最近的研究表明：一个三层的前馈神经网