水库多目标优化调度理论及其应用研究

1、水库多目的优化调度理论及其应用研究摘要：本文提出了综合利用水库的多目的优化调度的理论，并将该理论应用在综合利用水库优化调度过程中，在此应用中用马尔可夫单链弹性相关理论处理径流，并在引入“有效雨量的根底上，将供水量作为决策条件，以满足用水保证率条件下供水量最大为目的函数，建立了相应的数学模型和编制了相应的计算程序，绘出了综合利用水库三维优化调度图，利用三维优化调度图进展综合调节计算，计算结果理想、效益显著，且大大增加了调度过程的灵敏性。经沐浴水库等多个综合利用水库的理论证明，本方法是可靠有效的。关键词：优化调度弹性相关径流动态规划综合利用水库的优化调度受多因素影响，如径流，水库特性、用水特性以及

2、电站的机电特性等，其中径流的影响较大。本文采用马尔可夫单链弹性相关理论处理径流，以供水流量为决策变量，在考虑有效雨量的根底上建立了动态规划数学模型，编制了构造简明，功能完善，便于操作使用的大型优化调度计算程序，自动绘制出三维优化调度图，利用优化调度图进展综合利用水库调节计算，在几乎不增加投资的条件下，产生了宏大的经济效益。经理论证明，本方法准确可靠，合适于大、中、小型水库，也合适于平原水库、地下水库；更合适于我国北方水资源紧缺地区使用。1采用离散的马尔可夫随机过程描绘径流1.1用马尔可夫过程描绘径流为了计算和应用的方便，将时间序列离散化(即分为假设干时段：月)，相邻时段存在着依赖关系，以水库来

3、水的3个相邻时段t1、t2、t3间径流关系进展分析。用X1、X2、X3表示3个时段的径流，三者之间的相关情况可分为2种情况：(1)直接相关。即不管X2取值怎样(或不计X2取值的影响)的条件下，X1与X3相关，称为偏相关，其相关程度用相关系数表征，可用数量表示为13。(2)间接相关。即因存在着X1和X2、X2和X3之间的相邻时段相关关系，故X1的大小影响着X2的大小，从而又影响着X3的大校这种相关是由中间量X2传递的，不是直接的，因此叫间接相关。1.2计算相应条件概率当一年分成K个时段(月)，每个时段的径流以平均值来表示，记作QK(K=1,2,3,K)。应用相关理论分析，可以确定相邻时段径流QK

4、,QK-1(如图1所示)的条件概率分布函QK,QK-1的条件概率分布函数示意数F(QK/QK-1)。其条件概率分布是一个二维分布，用概率理论及水文统计原理来推求径流的条件概率计算式。图1相邻时段径流研究相邻时段的径流相关联络时，应用相关系数R及回归方程式求得(1)隔时段相关系数那么为：(2)式中：Q1i,Q2i,Q3i为第i年相邻时段的实测径流值；为平均值；n为径流实测系列年数。本时段径流的相关关系，应用相关中的直线相关，以自回归线性公式来表示：(3)式中：K,K-1分别为时段tk,tk-1的径流均方差；R1为相邻时段径流之间的相关系数。(4)式中：条件变差系数，其中vk为变差系数。一年划分为

5、K个时段，每个时段的径流划分为级(即个状态)，那么相邻时段的转移概率：Pkij(k=1,2,3,k;i,j=1,2,3,)表示的含义是tk-1时段径流为状态i时，tk时段径流为状态j时的概率而矩阵(5)那么表示tk-1时段到tk时段状态的转移概率矩阵，显然，这个矩阵的每行各非负元素之和为1，即：(6)为了计算Pkij转移概率的方便，取等分的10个概率5%，15%，95%，这样转移概率的值都为0.1，那么相应的条件概率的流量Qpi由式(4)即可求得。2动态规划动态规划法是美国数学家贝尔曼提出的，是一种研究多阶段决策过程的数学方法。近年来广泛应用于水资源规划管理领域中2.1动态规划数学模型把径流当

6、作随机过程的水库优化调度图的计算是一个多阶段的随机决策过程。它的计算模型如下。(1)阶段：将水库调度图按月(或者旬)划分成12个互相关连的阶段(时段)，以便求解(2)状态：因相邻两个阶段的入库平均流量Qt和Qt+1之间有相关关系，以面临时段初的库水位和本时段预报径流量Qt为状态变量St(Zt-1，Qt)(3)决策：在时段状态确定后，作一个相应的决定，即面临时段的供水量qt，同时确定了时段末水位，进展状态转移。水库水位分级，故有个状态转移，按0.618法在决策域内优选，对每一个状态变量St要选择一最优供水量qt，Stqt关系曲线为时段t的调度线，决策域为(QDin,t;Qxax,t)对决策变量供

7、水量qt进展所有状态优选计算时，还要进展库水位限制的检查判别，假设时段末蓄水量V2大于允许的最高蓄水位或限制水位，那么在水库蓄满前供水量仍按qt放水计算，当水库蓄满后那么按入库水量供水。当入库水量大于电厂最大过水才能时，超过部分作为弃水qt=q(Z初,Qt,tk)(7)tk为时段数，每一个决策就有一个相应的时段末库水位，水库进展了状态转移，假设将水库的水位划分为Z级，径流划分为级。一个时段的水库面临状态有Z种，全年水库运行状态有KZ种，水库优化调度图就是对全年各种运行状态作出相应决策变量的关系图。由式(7)可知，当时段tk的初始库水位和径流量已定时，时段的最优决策供水量是一个定值，因此下一时段

8、tK+1的初始库水位(即时段tk末的水位)也就是一个确定值。由于下一时段tK+1的径流不是一个确定值，而是依时段tK的径流Qt变化的随机值，其值由条件概率分布函数(弹性相关)决策。因此，水库在时段tK处于状态i，而时段tK+1处于状态j的状态转移概率为Pkij，那么有，而矩阵Pk=(Pkij)那么表示从时段tK到时段tK+1的水库状态转移概率矩阵，Pk完全由时段tK的调度方式和径流状态转移矩阵决定。经过多年运行后，水库的运行状态到达一个稳定的概率分布(5)效益函数：水库进展状态转移，伴随着产生了效益函数(包括了工业用水、生活用水、灌溉用水、发电用水及三个保证率)其中灌溉用水：因灌溉需水量每年、

9、每月、每天都不一样，因此是随机变量，极难编制计算机程序计算，故首次引入?农田水利学?的“有效雨量概念，使整个优化计算大大简化，完全解决了水量平衡问题，整个优化计算，水量平衡到达100%有效雨量的计算：从水库灌区试验站获取资料ij即从19521999年历年(i=19521999，j为第i年各月(或旬)的灌溉定额(是由历年灌溉试验站实测作物需水量采用通用电算程序计算出的)，而ax是48年中最枯水年的灌溉定额。ax-ij=P0ij,i=1952,1999,j=1,12,逐一列表进展计算。把每年每月的有效雨量加到每年每月的来水量Qt中，因ax是常数，所以仅有随机变量ij。其数学表达式如下：ixj=Ai

10、xj-Bixj，即：(8)式中ij为i年系列j时段(月)的有效雨量，aij为i年系列j时段农作物需水量(j可按日计算后归纳成各农作物生长期所需水量，再换算成月)。bij为i年系列j时段各类农作物综合灌溉水量。(6)目的函数：根据水库水资源缺乏的详细情况，拟定在满足生活用水和工业用水保证率的条件下，尽量满足农业用水。目的函数可表示为：满足用水量保证率条件下供水量最大。目的函数计算可用以下分段线性函数求得：f(st,qt)=qtQxaxqtQxin(9)f(st,qt)=qt+A(qt-Qxin)QxinqtQDinf(st,qt)=Qxax+E(qt-Qxax)QDaxqtQxax式中：qt为水

11、库供水量，QDin为系统供水下限，即保证城市生活用水和工业用水的下限；Qxin为农业保证供水量与QDin之和；QDax为电厂的最大过水才能；Qxax为农业供水量上限与QDin之和；E为发电专用水量小于Qxin时的折算系数，A为供水量小于Qxin时的折算系数，在计算中，可先任意假设A、E,A、E与Qxin的保证率成正比。给定一个A、E就可递推得出一张优化调度图，用水库多年入库流量资料按调度图进展历时操作计算，假设计算结果所得保证率低于要求的保证率，那么修改A、E重新递推计算(一般递推23次即可)，求得另一优化调度图，再进展历时操作，直至所得保证率符合要求为止。即经过试算选择满足保证率要求的A、E

12、值。2.2动态规划递推方程以qt为t阶段的决策变量，St(Zt-1,Qt)为t阶段的状态变量，那么其逆时序动态规划最优递推方程为：Ft(St,qt)=axft(St,qt)+Ft+1(St+1)qtQtt=1,2,N(10)式中：Ft(St,qt)代表水库从时刻t处于状态St出发至水库运行终了时刻N(计算周期末)的目的函数值；ft(St,qt)代表时刻t水库处于状态St取供水量qt时面临时段效益期望值；Ft+1(St+1)代表水库从时刻t+1处于St+1(j状态)出发至时刻期间各时段均采用最优决策时所得的效益期望值；Qt表示计算中t时段所用的入库径流序列；pi,j为t时刻采取qt决策，系统由第

13、t阶段的第i种状态St转移为第t+1阶段的第j种状态St+1时的条件概率，Ft+1相应St+1状态最优决策的效益。递推方程的约束条件如下：库水位约束Vin,tVtVax,t，即各时段的库水位不低于死水位Vin,t，也不能超过该时段允许的最高蓄水位Vax,t。水量平衡约束Vt+1=Vt+(Qt-qt)t-yt-Et,式中Vt+1、Vt代表时段t末、初的蓄水量；Qt、qt代表t时段平均入库径流量和供水量；yt为弃水量，Et为水库蒸发渗漏损失。供水约束和输水才能约束QDax,tqtQDin,t。t时段内供水量不能超过水轮机的最大过水才能QDax,t，也不能小于下限QDin,t2.3动态规划递推计算采

14、取逆时序逐时段动态规划递推计算，即每时段对所有状态逐一地优选对应的最优决策。对时段的多个入库流量代表值所产生的效益期望值。优选方法采用0.618法，规定搜索点为20个|Ft(Si)n-Ft(Si)(n+1)|/Ft(Si)(n+1)(11)式中：Ft(Si)n代表第n年时段t递推线上相应于状态Si的将来效益值；Ft(Si)(n+1)那么是第n+1年时段t递推线上同一状态Si相应的将来效益值，取0.001。一般最多递推两年就可以收敛，即可得出12时段或36个时段(旬)的最优调度线。这时各时段的最优决策构成一个最优策略，即为优化调度图。显然，因考虑月(或旬)、相隔月(旬)的相关，即多用了一项概率预

15、报，那么相应增加了经济效益。由于采用了马尔可夫单链弹性相关理论对径流进展处理，使水库调度图从二维坐标变成三维坐标，形成空间水库优化调度图，再由调度图换成一组以Qt为参数的方程，递推线也由一条变成一组，即优化调度线由一条线变成一组，形成一族调度曲线图，为便于实际调度时使用。2.5动态规划计算程序动态规划的计算是一个非常复杂的过程，不同的规划问题，要用不同的计算程序。我们根据最优化(pt)问题的数学模型2，用VISUL编制了计算程序，用递推方程找出最优解。该程序在P微机上调试成功，经理论证明其具有功能强大，使用方便，运行速度快等特点，并能自动绘出三维空间水库优化调度图及带有一组参数的调度曲线图。3

16、应用例如本方法已应用于山东沐寓跋山和黄前等几个大中型水库，都获得理想效果。仅以沐浴水库多目的优化调度的应用情况来说明。沐浴水库位于山东省烟台地区莱阳市，控制流域面积455k2，总库容1.87亿3。兴利库容1.07亿3，年平均来水量6900万3。水库每年向莱阳市供水180.0多万3，灌溉面积0.93万h2，水电站分东西电厂，装机容量共为1800k，是一座具有灌溉、防洪、城市工业、生活供水、发电、养殖等综合利用的大型水利工程。如图2所示。在沐浴水库优化调度过程中，我们用马尔可夫单链弹性相关理论对径流进展处理，将供水流量作为决策条件，在引入有效雨量的根底上，采用优选迭代试算来满足3个保证率(生活用水

17、保证率、工业用水保证率和灌溉用水保证率)的动态规划算法，协调了生活、工业、灌溉和发电之间的关系。图2沐浴水库运用系统示意应用满足用水保证率条件下供水量最大为目的函数合理地解决3个保证率的计算问题；建立了动态规划数学模型5，利用其优化调度程序计算，计算结果理想，输出了大量的表格，(如表1所示，限于篇幅，仅列一小部分)，并自动绘出了水库优化调度空间图及多族调度曲线图(如图3、4所示)。利用优化调度图进展综合调节计算，在几乎不增加投资的情况下，增加了宏大的经济效益。表1沐浴水库优化调度年序：1月份：8(单位：亿3)水位来水量(Q)0.63960.43680.32520.25910.21080.167

18、10.12690.09380.06160.0295最优决策水量(qt)63.0064.0065.00.81.0082.00.0.029500.046500.06650.0.122620.13155.0.029290.046170.06603.0.130630.05824.0.029090.045850.06557.0.129710.05784.0.028880.045530.06511.0.128800.05743.0.028680.045210.06466.0.127900.05703.0.028480.044900.06420.0.127010.05663.0.028280.044580.

19、06376.0.126120.05663.0.028080.044270.06331.0.125230.05584.0.027890.043960.06287.0.124360.05546.0.027690.043650.06243.0.123490.05506.年序：48月份：12(单位：亿3)水位来水量(Q)0.02230.01700.01340.01160.01070.00890.00630.00540.00450.0027最优决策水量(qt)63.0064.00.81.0082.000.002700.01545.0.014410.015450.002680.01535.0.015350.015350.002660.01524.0.015240.015240.002640.0153.0.015530.015530.002630.01503.0.015030.015030.002610.01492.0.014920.014920.002590.01482.0.014820.014820.002570.01471.0.014710.014710.002550.01461.0.014610.014610.002530.01451.0.014510.01451根据制定的水库优化调度图即马尔可夫调度图，对19521999年共48年水文年度的径流资料进展长系列操作计算，计算结果说明，