[发明专利]一种水库群多目标作用机制和优化调度方案分析方法有效
申请号: | 201911240272.8 | 申请日: | 2019-12-06 |
公开(公告)号: | CN111005346B | 公开(公告)日: | 2021-03-19 |
发明(设计)人: | 董增川;陈牧风;姚弘祎;贾文豪;倪效宽;吴振天;陈雨菲;钟加星 | 申请(专利权)人: | 河海大学 |
主分类号: | G06F30/13 | 分类号: | G06F30/13;E02B1/02 |
代理公司: | 南京苏高专利商标事务所(普通合伙) 32204 | 代理人: | 柏尚春 |
地址: | 210024 *** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 水库 多目标 作用 机制 优化 调度 方案 分析 方法 | ||
1.一种水库群多目标作用机制和优化调度方案分析方法,其特征在于,包括以下步骤:
(1)水库群系统概化:对水库群系统主要过程和影响因素进行概化分析,实现系统模拟,建立从系统实际情况到抽象数学表达的映射关系;所述水库群系统概化为:将社会经济、生态环境和水资源三系统概化为取退水节点、计算单元水量传输系统以及流域单元水量传输系统三种元素;
(2)对水库效益功能进行识别,确定效益目标的目标函数表达形式及水库调度约束条件,构建水库群发电-供水-生态-航运-防洪五目标优化调度模型,其中防洪目标基于其对水库调度要求的优先性和强制性,转换为约束条件予以考虑,确定模型目标函数与约束条件;
(3)模型输入数据确定及模型计算:模型输入数据考虑根据计算时段和计算目的需求,采用实测洪水资料、设计流量资料及历史来水过程资料,以天、旬、月作为计算步长;
(4)基于智能优化算法求解模型:基于步骤(1)中的水库群系统概化,以步骤(2)识别的四个效益目标为目标函数,以及确定的约束条件为约束,以步骤(3)中确定的来水资料为输入,作为建模结果及运算准备;
(5)解集可视化并基于相关性计算和置换率的多目标关系分析:利用Pearson相关性系数计算两个目标函数之间的相关性系数,并据此判断量目标函数间是否有显著的相关性;对于具有显著相关性的两目标,引入置换率的概念加以定量分析,对于两目标非劣前沿,目标间置换率比欧式在其他目标保持值不变的情况下,当第m个目标值改变一个单位时,第n个目标值增加或减少Tmn个单位,以补偿由于m改变而对总体效益的影响,其计算方法为:
式中,fm,fn为参与置换率分析的目标函数,λmn为两目标函数之间的置换率;计算时对两目标二维投影结果进行函数拟合,再求其一阶导数,即可得两目标间的置换率函数;
(6)基于方案间比较的多目标关系分析:基于模糊优选的方法,以各目标优化程度相同为标准得到方案集中的均衡方案,以该方案中的各目标值为基础,考察各效益目标最优方案下其余目标的变化情况;
(7)目标间作用关系判别及水库调度方式获取:
基于步骤(5)和步骤(6)结果得出对所考虑的四个目标间的作用关系:步骤(5)中有强相关性的目标间,根据其置换率所表现的一个目标效益增加,另一个目标效益减少,则认为两者有强竞争关系;另一个目标效益增加,则两者有强协同关系;步骤(5)中未表现强相关性的目标通过步骤(6)进行判断,相比均衡方案中各目标效益值,若两目标均在偏另一目标方案中效益减小或一目标效益减小、另一目标效益未变化,则两目标存在弱竞争关系;若两目标在偏另一目标方案中效益值一增一减,则两目标存在弱竞争或弱协同关系;若两目标均在偏另一目标方案中效益值增加或一目标效益增加、另一目标效益未变化,则两目标存在弱协同关系;基于优化结果可得各目标最优方案,以及四个效益均衡时水库群调度方案。
2.根据权利要求1所述的一种水库群多目标作用机制和优化调度方案分析方法,其特征在于,步骤(2)所述目标函数为:
水库群全年总发电量最大作为目标函数:
式中,E为梯级水库群发电总量,为第m个水电站在第t时段的出力,kW;M表示梯级电站的总个数,T表示总时段长度,年;Δt表示单位计算时间步长,月;
需水缺水量最小作为供水目标:
式中,K为社会经济需水区数目;Qkt表示时段t需水区k的总社会经济需水量对应计算时段流量,m3/s;
适宜生态缺溢水量最小为目标:
式中,L为河段生态控制断面总数;为第l断面在时段t的适宜生态流量,m3/s;Rlt为第l断面在时段t的实际流量,m3/s;
通航缺溢水量最小为目标函数:
式中,Qst为各个计算时段下泄流量与航运适宜流量区间差值绝对值,m3/s;SUapp、SLapp分别表示航运适宜流量区间的上界和下界,m3/s。
3.根据权利要求1所述的一种水库群多目标作用机制和优化调度方案分析方法,其特征在于,步骤(2)所述约束条件:
水量平衡约束:
Vm,t-Vm,t-1=(Im,t-Qm.t)×Δt
Im,t=Qm-1,t+Inm-1,t-Em,t
式中:Vm,t、Vm,t-1为m库第t时段末、初水库的库容,m3;Im,t为m库第t时段内平均入库流量,m3/s;Inm,t为第t时段m库与m+1库之间的区间入流,m3/s;Qm,t为m库第t时段平均出库流量,m3/s;Em,t为m库第t时段损失流量,主要包括水量在两库间传输的过程中由于蒸发和下渗而产生的水量损失,m3/s;Δt表示单位计算时间步长,月;
下泄流量约束:
Qmt,min≤Qmt≤Qmt,max
式中,Qmt,min表示m库t时段最小允许出库流量,m3/s;Qmt,max表示m库t时段最大允许出库流量,考虑水轮机组最大过流能力和汛期下游防洪对流量的要求,m3/s;
水库水位约束:
Zmt,min≤Zmt≤Zmt,max
式中,Zmt表示m库t时段末水位,m;Zmt,max表示m库t时段末最大允许水位,在汛期设定为汛限水位以在目标中考虑防洪目标的要求,在非汛期为正常蓄水位,m;Zmt,min表示m库t时段末最小允许水位,即为死水位,m;
水库出力约束:
Nmt,min≤Nmt≤Nmt,max
式中,Nmt表示m库t时段出力,万kW;Nmt,min表示m库t时段最小出力,万kW;Nmt,max表示m库t时段最大出力,即装机容量,万kW。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于河海大学,未经河海大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201911240272.8/1.html,转载请声明来源钻瓜专利网。