库群联合调度的轮库寻优程序 免费.docVIP

库群联合调度的轮库寻优程序 1?前言 随着地区经济的迅速发展，区域水资源短缺程度日趋严重，各地区在水资源分配与利用上的矛盾越发突出。如何挖掘区域内水资源潜力，充分利用现有水利工程，实现库群水资源的联合优化调度，具有较高的科研与实用价值。这里，以州河流域内三库水资源联合优化调度为例，编写轮库寻优计算程序，使现实问题得以解决，在库群联合调度方面进行了有益的探索。 2?流域状况与资料年限 2.1?流域状况 州河流域北部和东部与滦河为邻，南部和西部分别与泃河水系接壤。流域呈扇形，东西长约56 km，南北宽近50 km，总面积2 114 km2。流域内岩土以第四纪冲击层为主，北部山区多片麻岩，土壤为淋溶褐土，森林茂密，植被率达50%～60%；中部为冲积平原，土壤多为砂质粘土，农业发达。流域属大陆性季风气候，多年平均降水量750 mm。降水的年内、年际分配极不均匀。年内降水主要集中在6～9月份，占全年的80%以上。年际变化很大，丰水年多达1 256 mm?，枯水年仅有387 mm。 州河为蓟运河左支，上游干游为魏进河。自东向西有黎河、沙河、淋河等3条主要支流，集水面积分别为560 km2、887 km2、252 km2。三支流及干流经国河汇入于桥水库，出库后始称州河，向南汇流右支泃河后称蓟运河。三支流均为山溪性河流，地处燕山迎风坡暴雨中心区，每遇汛期暴雨，洪水较大。流域内主要水利工程为3座大、中型水库，即上关、般若院及于桥水库，详见图1。 图1?州河流域水利工程示意图 2.2?资料年限 所采用的水文资料及成果，为于桥水库1960—1999年、上关水库1974—2000年、般若院水库1973—2000年逐年各月入库径流量资料及其计算成果，流域内各水文站1960—1999年的观测资料。资料年限最长达40年，有利于对州河流域三水库水资源进行联合调度分析。 3?联合调度的数学模型 3.1?最优准则与目标函数 研究库群联合优化调度必先确定最优准则。这里根据州河流域三库水文特性及水资源利用现状，确定最优准则为：三水库控制并利用水资源量最大；三水库内水面蒸发量最小；于桥水库向下游排放水量最小。 上关、般若院及于桥水库自起始调度月起，水库各月平均水面面积分别记作：S11，S12，…，S112?；S21?，S22，…，S212?；S31?，S32，…，S312?。各月平均水面蒸发系数与蒸发量分别记作：K1?，K2，…，K12?。E1?，E2，…，E12?；于桥水库各月出库水量分别记作：QO31?，QO32，…，QO312?。则由最优准则可写出求解三库水量损失的数学表达式，亦即目标函数如下： Min∑[(S1i+S2i+S3i)·Ki·Ei/10+ QO3i?]　　　　i=1，2，…，12 3.2?约束方程 （1）蓄水量约束 Wji≤wji≤Wji　　j=1，2，3；?i=1，2，…，13 式中，wji为j水库第i个月份的库蓄水量；Wji、Wji分别为j水库第i个月份的允许最小、最大库蓄水量。 （2）需供水量约束 Wji+Uji≤wji　　j=1，2，3；?i=1，2，…，13 式中，Uji为j水库第i个月份需供水量，其它同上。 （3）起调蓄量约束 wj1=wj13　　j=1，2，3 式中，wj1、wj13分别为j水库第1个月份、第13个月份库蓄水量，均为起调蓄量。 3.3?数学模型 先确定三库联合调度的决策变量为：起调蓄量wj1、入库水量QIwji、需供水量Uwji、蒸发系数Ki、蒸发量Ei，则综合目标函数及约束方程，可获得州河流域三库联合调度的优化数学模型如下： MINEQ=Min∑[(S1i+S2i+S3i)·Ki·Ei?/10+QO3i] s.t. wj1=wj13 　　?Wji+Uji≤wji≤Wji 式中，j=1，2，3；i=1，2，…，13。 ? 4?轮库寻优算法与程序流程示意图 4.1?轮库寻优算法 由州河流域三库数学模型求取目标函数，可采用轮库寻优算法，其思路如下： （1）根据一般经验、分析判断或用其它简便方法，先给三水库定出起调蓄量Wj1。由起调蓄量，可对三水库分别确定一条满足约束条件且各月不超过最高蓄量的初始调度线：Wj1?，Wj2，…，Wj12，计算目标函数值。 （2）固定般若院水库的初始调度线，再将一定步长的蓄量变化△W1作为上关水库向于桥水库的放水量QO1，重新确定上关及于桥水库相应的蓄量调度线，并计算其目标函数值。之后，逐次对上关水库进行减量优化调度，比较各次目标函数值，记录最小水量损失及相应三库调度线，直至不满足上关水库的约束条件为止。 （3）将一定步长的蓄量变化△W2作为般若院水库向于桥水库的放水量QO2，并固定般若院水库减量后的蓄量调度线。再重复进行（2）中对上关水库的减量优化调度步骤。这样反复轮换优化于桥水库上游的上关及