库群联合调度的轮库寻优程序(1)

1、    库群联合调度的轮库寻优程序(1)    按一定目标对库群进行联合调度，具有较大的实用价值。州河地处津冀交界处，流域水资源对河北省地方经济发展及天津市用水十分重。运用轮库寻优算法编写程序，对州河流域三库进行联合调度，可有效地解决流域水资源优化利用问题。 关键词：联合调度 库寻优算法 程序 1 前言随着地区经济的迅速发展，区域水资源短缺程度日趋严重，各地区在水资源分配与利用上的矛盾越发突出。如何挖掘区域内水资源潜力，充分利用现有水利工程，实现库群水资源的联合优化调度，具有较高的科研与实用价值。这里，以州河流域内

2、三库水资源联合优化调度为例，编写轮库寻优计算程序，使现实问题得以解决，在库群联合调度方面进行了有益的探索。2 流域状况与资料年限2.1 流域状况州河流域北部和东部与滦河为邻，南部和西部分别与泃河水系接壤。流域呈扇形，东西长约56 km，南北宽近50 km，总面积2 114 km2。流域内岩土以第四纪冲击层为主，北部山区多片麻岩，土壤为淋溶褐土，森林茂密，植被率达50%60%；中部为冲积平原，土壤多为砂质粘土，农业发达。流域属大陆性季风气候，多年平均降水量750 mm。降水的年内、年际分配极不均匀。年内降水主集中在69月份，占全年的80%以上。年际变化很大，丰水年多达1 256 mm ，枯水年仅

3、有387 mm。州河为蓟运河左支，上游干游为魏进河。自东向西有黎河、沙河、淋河等3条主支流，集水面积分别为560 km2、887 km2、252 km2。三支流及干流经国河汇入于桥水库，出库后始称州河，向南汇流右支泃河后称蓟运河。三支流均为山溪性河流，地处燕山迎风坡暴雨中心区，每遇汛期暴雨，洪水较大。流域内主水利工程为3座大、中型水库，即上关、般若院及于桥水库，详见图1。图1 州河流域水利工程示意图2.2 资料年限所采用的水文资料及成果，为于桥水库19601999年、上关水库19742000年、般若院水库19732000年逐年各月入库径流量资料及其计算成果，流域内各水文站19601999年的观

4、测资料。资料年限最长达40年，有利于对州河流域三水库水资源进行联合调度分析。3 联合调度的数学模型3.1 最优准则与目标函数研究库群联合优化调度必先确定最优准则。这里根据州河流域三库水文特性及水资源利用现状，确定最优准则为：三水库控制并利用水资源量最大；三水库内水面蒸发量最小；于桥水库向下游排放水量最小。上关、般若院及于桥水库自起始调度月起，水库各月平均水面面积分别记作：S11 ，S12，S112 ；S21 ，S22，S212 ；S31 ，S32，S312 。各月平均水面蒸发系数与蒸发量分别记作：K1 ，K2，K12 。E1 ，E2，E12 ；于桥水库各月出库水量分别记作：QO31 ，QO32

5、，QO312 。则由最优准则可写出求解三库水量损失的数学表达式，亦即目标函数如下：Min(S1i S2i S3i)·Ki·Ei/10 QO3i i=1，2，123.2 约束方程（1）蓄水量约束WjiwjiWjij=1，2，3； i=1，2，13式中，wji为j水库第i个月份的库蓄水量；Wji、Wji分别为j水库第i个月份的允许最小、最大库蓄水量。（2）需供水量约束Wji Ujiwjij=1，2，3； i=1，2，13式中，Uji为j水库第i个月份需供水量，其它同上。（3）起调蓄量约束wj1=wj13j=1，2，3式中，wj1、wj13分别为j水库第1个月份、第13个月份库蓄

6、水量，均为起调蓄量。3.3 数学模型先确定三库联合调度的决策变量为：起调蓄量wj1、入库水量QIwji、需供水量Uwji、蒸发系数Ki、蒸发量Ei，则综合目标函数及约束方程，可获得州河流域三库联合调度的优化数学模型如下：MINEQ=Min(S1i S2i S3i)·Ki·Ei /10 QO3is.t. wj1=wj13 Wji UjiwjiWji 式中，j=1，2，3；i=1，2，13。4 轮库寻优算法与程序流程示意图4.1 轮库寻优算法由州河流域三库数学模型求取目标函数，可采用轮库寻优算法，其思路如下：      

7、;      摘按一定目标对库群进行联合调度，具有较大的实用价值。州河地处津冀交界处，流域水资源对河北省地方经济发         本篇论文是由3COME文档频道的网友为您在网络上收集整理饼投稿至本站的，论文版权属原作者，请不用于商业用途或者抄袭，仅供参考学习之用，否者后果自负，如果此文侵犯您的合法权益，请联系我们。（1）根据一般经验、分析判断或用其它简便方法，先给三水库定出起调蓄量Wj1。由起调蓄量，可对三水库分别确定一条满足约束条件且各月不超过最高蓄量

8、的初始调度线：Wj1 ，Wj2，Wj12，计算目标函数值。（2）固定般若院水库的初始调度线，再将一定步长的蓄量变化W1作为上关水库向于桥水库的放水量QO1，重新确定上关及于桥水库相应的蓄量调度线，并计算其目标函数值。之后，逐次对上关水库进行减量优化调度，比较各次目标函数值，记录最小水量损失及相应三库调度线，直至不满足上关水库的约束条件为止。（3）将一定步长的蓄量变化W2作为般若院水库向于桥水库的放水量QO2，并固定般若院水库减量后的蓄量调度线。再重复进行（2）中对上关水库的减量优化调度步骤。这样反复轮换优化于桥水库上游的上关及般若院两水库，直至不满足般若院水库的约束条件为止。（4）按上述步骤，

9、可计算出相应步长三库所有蓄量调度线的目标函数值，进而寻得最小目标函数值（即最少损失水量）及三库相应的蓄量调度线、放水量。4.2 程序流程示意图对上述轮库寻优算法的求解步骤，设计其计算机程序流程示意图，见图2。图2 州河流域三库优化调度程序流程示意图    4.3 程序清单分别用Wij、Sij、QIij 、QOij、UWij表示i水库第j个月份的蓄水量、库水面面积、入库水量、出库水量、需供水量，用DWi 、steps表示i水库的死库容、寻优步长，用Kj、Ej分别表示第j个月份的蒸发系数、蒸发量，用Min_EQ （MINEQ）、Min_E（MINE）、Min

10、_Q（MINQ）分别表示三库最小水量损失、最小蒸发水量、于桥水库最小出库水量，用WLij、QOLij表示i水库第j个月份的最优蓄水量、出库水量。对轮库寻优算法的求解步骤，则可用C语言编程如下：int i,j; /初始化寻优起始状态for(i=1;i<4;i )for(j=1;j<14;j )WLij=0;MINEQ= 99999999; MINE= 99999999; MINQ=99999999;W11=?，W21=?，W31=?; /轮库寻优算法int hh=0,tt=0;for(int fp=0;fp )/按一定步长求生成般若院水库蓄水量调度线for(int g=2;g<

11、14;g ) W2g=0;for(int gg=1;gg<13;gg ) QO2gg=0;for(int h=2;h<14;h )/过量蓄量调整QO2h-1=steps*fp;W2h=W2h-1 QI2h-1-steps*fp;if(W2h>2800)QO2h-1=QO2h-1 W2h-2800;W2h=2800;if(W213>W21)QO212=W213-W21 QO212;W213=W21;for(i=1;i<13;i )/判断并结束般若院水库过程线if(W2i< DW2 UWi) break;if(W213<> for(int fs=0;

12、fs )/按一定步长求生成上关水库蓄水量调度线for(g=2;g<14;g )W1g=0;for(gg=1;gg<13;gg )QO1gg=0;for(int t=2;t<14;t )QO1t-1=steps*fs;W1t=W1t-1 QI1t-1-steps*fs;if(W1t>2960)QO1t-1=QO1t-1 W1t-2960;W1t=2960; if(W113>W11)QO112=W113-W11 QO112;W113=W11;for(i=1;i<13;i )/判断并结束上关水库过程线if(W1i< DW1 UWi)break;if(W113

13、<> for(g=2;g<14;g )W3g=0;for(gg=1;gg<13;gg ) QO3gg=0;for(int r=2;r<14;r )/生成于桥水库蓄水量过程线W3r=W3r-1 app->QI3r-1 QO1r-1 QO2r-1;if(W3r>59791)QO3r-1=W3r-59791;W3r=59791;if(W313>W31)QO312=QO312 (W313-W31);W313=W31; /计算并判断最小水量损失，记录相应库状态for(j=1;j<13;j )/库水面面积与库容关系模拟S1j= 0.1*pow(10,-

14、9)*pow(W1j,3) - 0. 6*pow(10,-6)*pow(W1j,2) 0.0014*W1j 0.2583; S2j=0.6*pow(10,-10)*pow(W2j,3)-0.5*pow(10,-6)*pow(W2j,2) 0.0021*W2j 0.2345;            摘按一定目标对库群进行联合调度，具有较大的实用价值。州河地处津冀交界处，流域水资源对河北省地方经济发       

15、;  本篇论文是由3COME文档频道的网友为您在网络上收集整理饼投稿至本站的，论文版权属原作者，请不用于商业用途或者抄袭，仅供参考学习之用，否者后果自负，如果此文侵犯您的合法权益，请联系我们。S3j=0.5*pow(10,-12)*pow(W3j,3)-0.7*pow(10,-7)*pow(W3j,2) 0.0043*W3j 6.8841; Min_EQ=0，Min_E=0，Min_Q=0;for(int jj=1;jj<13;jj )/计算目标函数Min_EQ=Min_EQ (SO1jj SO2jj SO3jj)* Kjj*Ejj/10 QO3jj;Min_E=Min_E (

16、SO1jj SO2jj SO3jj)*Kjj*Ejj/10;Min_Q=Min_Q QO3jj;if(Min_EQB1< > MINEQ=Min_EQ; MINE=Min_E; MINQ=Min_Q;for(i=1;i<4;i )for(j=1;j<13;j )WLij=Wij;QOLij=QOij;/上关水库优化轮库结束/般若院水库优化轮库结束/输出寻优成果：MINEQ,MINE,MINQ,steps,WLij,QOLij;    5 三库联合优化调度成果由州河流域三库多年水文资料，可确定三水库在枯水年、平水年、丰水年的起调蓄量

17、、入库水量、需供水量、蒸发系数及蒸发量等参数。这里仅以丰水年的调度为例进行运算，并最终确定三水库在最优调度时的相应状态。根据已知资料，可确定上关、般若院及于桥水库在丰水年的起调蓄量分别为：1880、2455、29777（104m3），三库入库水量(表1)，需供水量(表.2)，蒸发系数及蒸发量(表.3)。运用轮库寻优算法，可算得实现联合调度时三库丰水年最小水量损失为81915.23(104m3)，最小蒸发量10960.25 (104m3)，于桥水库最小放水量70954.98(104m3)。三库优化调度方案中蓄量、水位及水库放水量(表.4)。6 结语州河流域尽管有较强的产流能力，却因流域面积以及上

18、关、般若院两水库库容相对较小，在防洪中同于桥水库进行联合调度的难度较大。这里，从径流调节及流域水资源联合调度角度出发，运用数学模型与轮库寻优算法，编写计算机程序对流域内三库进行寻优计算，使流域水资源损失减至最小，也为库群联合调度提供了实用的计算方法。参考文献：1 库群优化调度J.南京：河海大学，1997.PROGRAM SEEKING FOR THE OPTIMUM MANAGEMENT FOR A GROUP OF RESERVIORSAbstract: Zhou River locates the place where Tianjin and Hebei meet,and the basins water resources are important for Hebeis economic development and TianJins water supply. So, to gain some goals, united management for a group of reservoirs has great practical value. With this program to