行
业
产
品
功
能
北极星电力网技术频道 作者:佚名 2008/1/16 18:20:58
摘要：本文提出了综合利用水库的多目标优化调度的理论，并将该理论应用在综合利用水库优化调度过程中，在此应用中用马尔可夫单链弹性相关理论处理径流，并在引入“有效雨量”的基础上，将供水量作为决策条件，以满足用水保证率条件下供水量最大为目标函数，建立了相应的数学模型和编制了相应的计算程序，绘出了综合利用水库三维优化调度图，利用三维优化调度图进行综合调节计算，计算结果理想、效益显著，且大大增加了调度过程的灵活性。经沐浴水库等多个综合利用水库的实践证明，本方法是可靠有效的。
关键词：优化调度弹性相关径流动态规划
综合利用水库的优化调度受多因素影响，如径流，水库特性、用水特性以及电站的机电特性等，其中径流的影响较大。本文采用马尔可夫单链弹性相关理论处理径流，以供水流量为决策变量，在考虑有效雨量的基础上建立了动态规划数学模型，编制了结构简明，功能完善，便于操作使用的大型优化调度计算程序，自动绘制出三维优化调度图，利用优化调度图进行综合利用水库调节计算，在几乎不增加投资的条件下，产生了巨大的经济效益。经实践证明，本方法准确可靠，适合于大、中、小型水库，也适合于平原水库、地下水库；更适合于我国北方水资源紧缺地区使用。
1.1用马尔可夫过程描述径流
为了计算和应用的方便，将时间序列离散化(即分为若干时段：月)，相邻时段存在着依赖关系，以水库来水的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,qk-1(如图1所示)的条件概率分布函qk,qk-1的条件概率分布函数示意数f(qk/qk-1)。其条件概率分布是一个二维分布，用概率理论及水文统计原理来推求径流的条件概率计算式。
相邻时段径流之间应用自回归线性相关时，其间隔时段的径流对回归线的偏离值即误差的分布，经刚性和弹性相关比较后，采用了弹性相关处理方法即偏态分布，按皮尔逊ⅲ型曲线分布。相应于条件概率的流量qpk可由下式求得：
而矩阵
动态规划法是美国数学家贝尔曼提出的，是一种研究多阶段决策过程的数学方法。近年来广泛应用于水资源规划管理领域中
2.1动态规划数学模型
(1)阶段：将水库调度图按月(或者旬)划分成12个相互关连的阶段(时段)，以便求解
(2)状态：因相邻两个阶段的入库平均流量qt和qt 1之间有相关关系，以面临时段初的库水位和本时段预报径流量qt为状态变量st(zt-1，qt)
(3)决策：在时段状态确定后，作一个相应的决定，即面临时段的供水量qt，同时确定了时段末水位，进行状态转移。水库水位分m级，故有m个状态转移，按0.618法在决策域内优选，对每一个状态变量st要选择一最优供水量qt，st～qt关系曲线为时段t的调度线，决策域为(qdmin,t;qxmax,t)
对决策变量供水量qt进行所有状态优选计算时，还要进行库水位限制的检查判别，若时段末蓄水量v2大于允许的最高蓄水位或限制水位，则在水库蓄满前供水量仍按qt放水计算，当水库蓄满后则按入库水量供水。当入库水量大于电厂最大过水能力时，超过部分作为弃水
(4)状态转移：水库状态和调度图形式有关，因考虑当时入库径流和短期径流因素，水库调度中将一年划分为k个时段，每个时段由时段初库水位z初和时段流量qt组成水库的运行状态，而每一种状态有一个相应的决策变量供水流量qt，用函数关系表示为：
由式(7)可知，当时段tk的初始库水位和径流量已定时，时段的最优决策供水量是一个定值，因而下一时段tk 1的初始库水位(即时段tk末的水位)也就是一个确定值。由于下一时段tk 1的径流不是一个确定值，而是依时段tk的径流qt变化的随机值，其值由条件概率分布函数(弹性相关)决策。因此，水库在时段tk处于状态i，而时 段tk 1处于状态j的状态转移概率为pkij，则有，而矩阵pk=(pkij)则表示从时段tk到时段tk 1的水库状态转移概率矩阵，pk完全由时段tk的调度方式和径流状态转移矩阵决定。经过多年运行后，水库的运行状态达到一个稳定的概率分布
(5)效益函数：水库进行状态转移，伴随着产生了效益函数(包括了工业用水、生活用水、灌溉用水、发电用水及三个保证率)
其中灌溉用水：因灌溉需水量每年、每月、每天都不相同，因此是随机变量，极难编制计算机程序计算，故首次引入《农田水利学》的“有效雨量”概念，使整个优化计算大大简化，完全解决了水量平衡问题，整个优化计算，水量平衡达到100
有效雨量的计算：从水库灌区试验站获取资料mij即从1952～1999年历年(i=1952～1999，j为第i年各月(或旬))的灌溉定额(是由历年灌溉试验站实测作物需水量采用通用电算程序计算出的)，而mmax是48年中最枯水年的灌溉定额。mmax-mij=p0ij,i=1952,…,1999,j=1,…,12,逐一列表进行计算。把每年每月的有效雨量加到每年每月的来水量qt中，因mmax是常数，所以仅有随机变量mij。其数学表达式如下：cixj=aixj-bixj，即：
(6)目标函数：根据水库水资源不足的具体情况，拟定在满足生活用水和工业用水保证率的条件下，尽量满足农业用水。目标函数可表示为：满足用水量保证率条件下供水量最大。目标函数计算可用下列分段线性函数求得：
2.2动态规划递推方程以qt为t阶段的决策变量，st(zt-1,qt)为t阶段的状态变量，则其逆时序动态规划最优递推方程为：
递推方程的约束条件如下：①库水位约束vmin,t≤vt≤vmax,t，即各时段的库水位不低于死水位vmin,t，也不能超过该时段允许的最高蓄水位vmax,t。②水量平衡约束vt 1=vt (qt-qt)·δt-yt-et,式中vt 1、vt代表时段t末、初的蓄水量；qt、qt代表t时段平均入库径流量和供水量；yt为弃水量，et为水库蒸发渗漏损失。③供水约束和输水能力约束qdmax,t≥qt≥qdmin,t。t时段内供水量不能超过水轮机的最大过水能力qdmax,t，也不能小于下限qdmin,t
2.3动态规划递推计算采取逆时序逐时段动态规划递推计算，即每时段对所有状态逐一地优选对应的最优决策。对时段的多个入库流量代表值所产生的效益期望值。优选方法采用0.618法，规定搜索点为20个
2.4优化调度图howard用z变换方法证明式(10)随年数t增加计算是收敛的，进行递推计算采取逆时序递推，即从n时段开始递推到1时段，只要知道fn(sn)即可按式(10)递推计算。开始可取库水位(库容)～蓄水量关系曲线作为初始递推线fn(sn)。当对第一个时段的所有状态优选出最优决策后，即可往前递推一个时段。当第一年逐个时段全部递推计算完毕后，还要进行第二年周期的递推计算，是因为初始递推fn(sn)是任意假设的，故第一年周期递推所得的策略并非稳定的最优策略，必需继续递推至各时段的递推线均收敛为止，这时所得的策略才是稳定的最优策略。递推线收敛的准则是：前后两年周期中同一时段的递推线相差小于规定的相对误差ε即：
2.5动态规划计算程序动态规划的计算是一个非常复杂的过程，不同的规划问题，要用不同的计算程序。我们根据最优化(opt)问题的数学模型［2］，用visulc编制了计算程序，用递推方程找出最优解。该程序在pⅱ微机上调试成功，经实践证明其具有功能强大，使用方便，运行速度快等特点，并能自动绘出三维空间水库优化调度图及带有一组参数的调度曲线图。
参考文献：
［1］张勇传.水电站优化调度［m］.北京：水利电力出版社，1983.
［2］魏权，等.数学规划与优化调度［m］.北京：水利电力出版社，1984.
［3］廖昭懋，杨文礼.概率论与数理统计［m］.北京：师范大学出版社，1988.
［4］loucksdp,stedingerjr,haithda.waterresourceresourcesystemsplanningandanalysis［m］.cornelluniversity,1981.
［5］［美］n.伯拉斯著，戴国瑞，冯尚友译.水资源科学分配［m］.北京：水利电力出版社，1983
本方法已应用于山东沐浴、跋山和黄前等几个大中型水库，都取得理想效果。仅以沐浴水库多目标优化调度的应用情况来说明。
在沐浴水库优化调度过程中，我们用马尔可夫单链弹性相关理论对径流进行处理，将供水流量作为决策条件，在引入有效雨量的基础上，采用优选迭代试算来满足3个保证率(生活用水保证率、工业用水保证率和灌溉用水保证率)的动态规划算法，协调了生活、工业、灌溉和发电之间的关系。
64.00
65.00
...
81.00
82.00
...
0.04650
0.06650
...
0.12262
0.13155
...
0.04617
0.06603
...
0.13063
0.05824
...
0.04585
0.06557
...
0.12971
0.05784
...
0.04553
0.06511
...
0.12880
0.05743
...
0.04521
0.06466
...
0.12790
0.05703
...
0.04490
0.06420
...
0.12701
0.05663
...
0.04458
0.06376
...
0.12612
0.05663
...
0.04427
0.06331
...
0.12523
0.05584
...
0.04396
0.06287
...
0.12436
0.05546
...
0.04365
0.06243
...
0.12349
0.05506
...
64.00
...
81.00
82.00
0.01545
...
0.01441
0.01545
0.01535
...
0.01535
0.01535
0.01524
...
0.01524
0.01524
0.0153
...
0.01553
0.01553
0.01503
...
0.01503
0.01503
0.01492
...
0.01492
0.01492
0.01482
...
0.01482
0.01482
0.01471
...
0.01471
0.01471
0.01461
...
0.01461
0.01461
0.01451
...
0.01451
0.01451
广告直拨： 媒体合作/投稿：陈女士 13693626116
关于北极星 | 广告服务 | 会员服务 | 媒体报道 | 营销方案 | 成功案例 | 招聘服务 | 加入我们 | 网站地图 | 联系我们 | 排行
京ICP证080169号京ICP备09003304号-2
京公网安备11010502034458号电子公告服务专项备案
网络文化经营许可证 [2019] 5229-579号广播电视节目制作经营许可证 (京) 字第13229号出版物经营许可证新出发京批字第直200384号人力资源服务许可证1101052014340号
Copyright © 2022 Bjx.com.cn All Rights Reserved. 北京火山动力网络技术有限公司 版权所有
