水利工程论文-库群联合调度的轮库寻优程序.doc_第1页
水利工程论文-库群联合调度的轮库寻优程序.doc_第2页
水利工程论文-库群联合调度的轮库寻优程序.doc_第3页
水利工程论文-库群联合调度的轮库寻优程序.doc_第4页
水利工程论文-库群联合调度的轮库寻优程序.doc_第5页
已阅读5页,还剩4页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

水利工程论文-库群联合调度的轮库寻优程序摘要:按一定目标对库群进行联合调度,具有较大的实用价值。州河地处津冀交界处,流域水资源对河北省地方经济发展及天津市用水十分重要。运用轮库寻优算法编写程序,对州河流域三库进行联合调度,可有效地解决流域水资源优化利用问题。关键词:联合调度库寻优算法程序1前言随着地区经济的迅速发展,区域水资源短缺程度日趋严重,各地区在水资源分配与利用上的矛盾越发突出。如何挖掘区域内水资源潜力,充分利用现有水利工程,实现库群水资源的联合优化调度,具有较高的科研与实用价值。这里,以州河流域内三库水资源联合优化调度为例,编写轮库寻优计算程序,使现实问题得以解决,在库群联合调度方面进行了有益的探索。2流域状况与资料年限2.1流域状况州河流域北部和东部与滦河为邻,南部和西部分别与泃河水系接壤。流域呈扇形,东西长约56km,南北宽近50km,总面积2114km2。流域内岩土以第四纪冲击层为主,北部山区多片麻岩,土壤为淋溶褐土,森林茂密,植被率达50%60%;中部为冲积平原,土壤多为砂质粘土,农业发达。流域属大陆性季风气候,多年平均降水量750mm。降水的年内、年际分配极不均匀。年内降水主要集中在69月份,占全年的80%以上。年际变化很大,丰水年多达1256mm,枯水年仅有387mm。州河为蓟运河左支,上游干游为魏进河。自东向西有黎河、沙河、淋河等3条主要支流,集水面积分别为560km2、887km2、252km2。三支流及干流经国河汇入于桥水库,出库后始称州河,向南汇流右支泃河后称蓟运河。三支流均为山溪性河流,地处燕山迎风坡暴雨中心区,每遇汛期暴雨,洪水较大。流域内主要水利工程为3座大、中型水库,即上关、般若院及于桥水库,详见图1。图1州河流域水利工程示意图2.2资料年限所采用的水文资料及成果,为于桥水库19601999年、上关水库19742000年、般若院水库19732000年逐年各月入库径流量资料及其计算成果,流域内各水文站19601999年的观测资料。资料年限最长达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)KiEi/10+QO3ii=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个月份库蓄水量,均为起调蓄量。3.3数学模型先确定三库联合调度的决策变量为:起调蓄量wj1、入库水量QIwji、需供水量Uwji、蒸发系数Ki、蒸发量Ei,则综合目标函数及约束方程,可获得州河流域三库联合调度的优化数学模型如下:MINEQ=Min(S1i+S2i+S3i)KiEi/10+QO3is.t.wj1=wj13Wji+UjiwjiWji式中,j=1,2,3;i=1,2,13。4轮库寻优算法与程序流程示意图4.1轮库寻优算法由州河流域三库数学模型求取目标函数,可采用轮库寻优算法,其思路如下:(1)根据一般经验、分析判断或用其它简便方法,先给三水库定出起调蓄量Wj1。由起调蓄量,可对三水库分别确定一条满足约束条件且各月不超过最高蓄量的初始调度线:Wj1,Wj2,Wj12,计算目标函数值。(2)固定般若院水库的初始调度线,再将一定步长的蓄量变化W1作为上关水库向于桥水库的放水量QO1,重新确定上关及于桥水库相应的蓄量调度线,并计算其目标函数值。之后,逐次对上关水库进行减量优化调度,比较各次目标函数值,记录最小水量损失及相应三库调度线,直至不满足上关水库的约束条件为止。(3)将一定步长的蓄量变化W2作为般若院水库向于桥水库的放水量QO2,并固定般若院水库减量后的蓄量调度线。再重复进行(2)中对上关水库的减量优化调度步骤。这样反复轮换优化于桥水库上游的上关及般若院两水库,直至不满足般若院水库的约束条件为止。(4)按上述步骤,可计算出相应步长三库所有蓄量调度线的目标函数值,进而寻得最小目标函数值(即最少损失水量)及三库相应的蓄量调度线、放水量。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_Q(MINQ)分别表示三库最小水量损失、最小蒸发水量、于桥水库最小出库水量,用WLij、QOLij表示i水库第j个月份的最优蓄水量、出库水量。对轮库寻优算法的求解步骤,则可用C语言编程如下:inti,j;/初始化寻优起始状态for(i=1;i4;i+)for(j=1;j14;j+)WLij=0;MINEQ=99999999;MINE=99999999;MINQ=99999999;W11=?,W21=?,W31=?;/轮库寻优算法inthh=0,tt=0;for(intfp=0;fp+)/按一定步长要求生成般若院水库蓄水量调度线for(intg=2;g14;g+)W2g=0;for(intgg=1;gg13;gg+)QO2gg=0;for(inth=2;h2800)QO2h-1=QO2h-1+W2h-2800;W2h=2800;if(W213W21)QO212=W213-W21+QO212;W213=W21;for(i=1;i13;i+)/判断并结束般若院水库过程线if(W2iDW2+UWi)break;if(W213for(intfs=0;fs+)/按一定步长要求生成上

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论