CN109886473B - 一种考虑下游生态的流域风光水***多目标优化调度方法 - Google Patents
一种考虑下游生态的流域风光水***多目标优化调度方法 Download PDFInfo
- Publication number
- CN109886473B CN109886473B CN201910067807.XA CN201910067807A CN109886473B CN 109886473 B CN109886473 B CN 109886473B CN 201910067807 A CN201910067807 A CN 201910067807A CN 109886473 B CN109886473 B CN 109886473B
- Authority
- CN
- China
- Prior art keywords
- wind
- output
- ith
- reservoir
- period
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 121
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005457 optimization Methods 0.000 title claims abstract description 29
- 238000010248 power generation Methods 0.000 claims abstract description 32
- 238000013178 mathematical model Methods 0.000 claims abstract description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 16
- 230000005855 radiation Effects 0.000 claims description 16
- 230000035772 mutation Effects 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 230000005611 electricity Effects 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 10
- 230000013011 mating Effects 0.000 claims description 8
- 239000011717 all-trans-retinol Substances 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000009827 uniform distribution Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- FPIPGXGPPPQFEQ-OVSJKPMPSA-N all-trans-retinol Chemical compound OC\C=C(/C)\C=C\C=C(/C)\C=C\C1=C(C)CCCC1(C)C FPIPGXGPPPQFEQ-OVSJKPMPSA-N 0.000 claims description 3
- 238000013486 operation strategy Methods 0.000 claims description 3
- 238000010187 selection method Methods 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 12
- 238000004458 analytical method Methods 0.000 abstract description 7
- 230000008901 benefit Effects 0.000 description 23
- 230000002860 competitive effect Effects 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 4
- 230000002441 reversible effect Effects 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 241001048891 Jatropha curcas Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 235000010678 Paulownia tomentosa Nutrition 0.000 description 1
- 240000002834 Paulownia tomentosa Species 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑下游生态的流域风光水***多目标优化调度方法,包括获取梯级水电站的基本信息;获取流域风光电站的基本信息;根据流域风光水电站基本信息,建立考虑下游生态的流域风光水***目标以及硬性约束的多目标发电优化调度数学模型;采用改进的AR‑MOEA对建立的多目标发电优化调度数学模型进行求解;采用SMAA‑2对得到的非支配方案进行多属性决策。本发明将风光水电***进行解耦,避免可行解在在进化过程中遭到破坏;随后得到完整Pareto前沿,应用能考虑决策者主观意愿和属性信息不确定性的随机多准则可接受性分析(SMAA‑2)的多属性决策模型,给决策者提供了更多方案排序的信息,能更科学合理地进行决策。
Description
技术领域
本发明涉及新能源和水利水电调度方法,特别涉及一种考虑下游生态的流域风光水***多目标优化调度方法。
背景技术
风电和光伏发电出力具有随机性、波动性、和间歇性等特点,直接上网给电力***的运行带来困扰,而水电具有可储存、调节速度快,能有效缓解间歇性能源出力波动给电力***带来的影响。因此,研究水电补偿风光出力具有重要意义。
目前流域风光水电***联合调度的研究大多集中在***效益最大、***出力稳定性和最小化***运行费用。但是,水电站一般都承担着多个任务,如生态、发电和防洪等。水电站根据风光的不稳定出力适应性调整自身运行方式补偿风光出力,不仅使得自身的运行状态发生了改变,进而对下游河道生态造成影响。因此,建立能满足风光水电站稳定运行、风光水***效益和下游河道生态健康的多目标优化调度模型是流域风光水***管理发展的必然趋势。
目前,处理多目标问题的方法通常可以分为三类,一是约束法,将目标变为约束条件,使多目标问题转化为单目标问题进行优化求解。二是权重法,通过一组权重值将多目标组合成单目标问题进行优化求解。三是近年来出现的多目标进化算法,包括动态权重进化算法、NSGA-II、(NSGA-Ⅲ)和AR-MOEA等。其中AR-MOEA可以在无偏好的情况下运行一次得到完整的Pareto前沿,且对于Pareto前沿不连续、不可微、非凸等情况均具有较好的鲁棒性。
风光水电***多目标优化调度涉及多部门(电力、水利、生态等部门)的利益,不同利益部门对于各个目标之间的偏好是不一样的,给出的主观权重往往不一致。客观权重可以反应方案集数据自身的属性信息,不同赋权方法得到的权重差异较大。因此,在多目标优化得到Pareto最优解之后,建立考虑指标权重不确定性的随机多属性决策是至关重要的,它能给决策者提供更多方案排序的信息,从而能够科学合理地进行决策。
发明内容
发明目的:本发明提供一种考虑下游生态的流域风光水***多目标优化调度方法,以解决现有技术存在的上述问题。
技术方案:为实现上述发明目的,本发明采用以下技术方案:
一种考虑下游生态的流域风光水***多目标优化调度方法,包括以下步骤:
(1)获取梯级水电站的基本信息,该基本信息包括各个水电站的入库径流资料、水位库容曲线、尾水流量关系曲线和水轮机机组综合曲线资料;
(2)获取流域风光电站的基本信息,该基本信息包括各个风光电站的地理位置、风速、太阳辐射、气温和装机容量资料;
(3)根据流域风光水电站基本信息,建立考虑下游生态的流域风光水***目标以及硬性约束的多目标发电优化调度数学模型;
(4)采用改进的AR-MOEA对所述步骤(3)建立的多目标发电优化调度数学模型进行求解;
(5)采用SMAA-2对步骤(4)得到的非支配方案进行多属性决策。
进一步的,步骤(3)中多目标优化调度数学模型的目标函数为发电量、最大调度期内最小时段出力最大和对下游生态影响最小;
发电量最大目标函数为:
其中,f1为风光水电***总发电量;T为调度期时段数;PPVt、PWt和PHt分别是光伏电站、风力发电站和水电站第t时段出力;Δt为时段小时数;
最小时段出力最大目标函数为:
其中,f2为风光水电***最小出力;Pt=PPVt+PWt+PHt;
水电站出库流量与天然流量偏差最小目标函数为:
其中,f3为水电站的出库流量与天然流量的偏差;Oi',t为第i水库t时段的出库流量,QNi,t为第i库t时段下游断面天然流量;I为水库数目。
进一步的,步骤(3)中硬性约束具体为:
(a)各库的水量平衡约束
Vi,t=Vi,t-1+(Qi,t-O'i,t)Δt;
(b)各库的上、下限水位约束
(c)流量约束
(d)调度期末水位约束
Zi,T=Zi,end,i∈[1,I];
其中,Zi,end为第i库调度期的期末水位;
(e)光伏电站出力约束
(f)风力发电站出力约束
(g)水电站出力约束
(h)外送断面约束
进一步的,步骤(3)中风光水电***出力计算具体为:
(a)风力发电站出力计算方式如下:
其中,SA是风力发电机轮毂的面积,ρ是空气密度,Nk是第风力发电站的风力发电机的台数,uk,t是风力发电机轮毂处的风速,具体见下式:
(b)光伏电站出力计算:建立光伏发电***与太阳辐射和温度的关系;
其中,Pstc是标准条件下光伏板的出力,Gstc是标准条件下的太阳辐射强度,是第d光伏电站第t时段的光伏板的温度,Tref是标准条件下的温度,是第d光伏电站t时段的实际太阳辐射强度,β是光伏板的温度系数,γ是光伏板的辐射系数,是第d光伏电站光伏板的面积;
(c)水电站的出力如下式所示:
PHi,t=g(Oi,t,ΔHi,t);
其中,PHi,t是第i水电站第t时段出力,g(·)函数为水电站出力特性函数,Oi,t是发电流量,ΔHi,t是发电水头。
进一步的,步骤(4)中改进的AR-MOEA算法具体为:
(41)读取风速、太阳辐射和气温资料,根据根据风力发电站出力计算公式和光伏电站出力计算公式计算风、光出力;
(42)随机生成M个初始个体;选择水库水位作为决策变量进行实数编码,具体的生成方式如下:
(43)将初始种群复制到外部档案AC,并且生成一个大小为NR的初始参考点集;
(44)将步骤(41)计算得到的风、光***的出力输入到梯级水电站***中,计算种群的适应度函数和IGD-NS贡献度,并采用锦标赛选择法选择贡献度大的个体到交配池,IGD-NS指标具体为:
(45)由于梯级水电站上、下游水库存在水量联系,各库运行策略相互影响,存在着多变量耦合约束,直接对个体进行交叉和变异不能保证子代个体满足约束条件;因此,在交叉、变异前基于梯级水电站复杂约束建立动态可行域,具体为:
满足所有约束的动态可行域为:
其中,FRj为第j个约束条件的可行域,CountC是可能被破坏的约束条件个数;
确定完动态可行域之后,在动态可行域内进行交叉和变异,具体如下所示:
改进的变异算子为:
(46)基于步骤(45)在动态可行域产生的子代个体,进行更新外部档案AC,然后利用更新后的外部档案AC进行自适应调整参考点R',使得参考点R'与Pareto前沿的形状相似;
(47)合并父代与子代种群,并对其进行非支配排序,然后将基于IGD-NS指标的环境选择策略作用于第kM个前沿面的解上,直到k满足|Front1∪...∪FrontkM|≥M;
(48)判断是否达到最大迭代次数,若是,退出循环;若否,重复步骤(44)到(47),直至达到最大迭代次数。
更进一步的,步骤(46)中通过选择交配池种群P'并经过交叉和变异得到子代后,利用新产生的解来更新外部档案AC和自适应参考点集合R',具体步骤如下:
1)基于当前种群P各个目标中的理想点和最差点N为目标数,外部档案AC、种群P和参考点R各个目标分别同时减去并同时乘以使得转换之后种群P、外部档案AC和参考点R标准化到相同的范围使得均匀分布的参考点R能够在不同的目标空间内生成均匀分布的方案;
2)基于标准化后的参考点R,计算外部档案AC的IGD-NS贡献值,删除外部档案AC里互相支配的和冗余的方案,并将剩余的方案复制到外部档案Acon;
3)将外部档案Acon复制到新的外部档案AC',基于新的外部档案AC'计算参考点R的IGD-NS值,并将有贡献的参考点保留,然后从新的外部档案AC'剩余的解选择最不拥挤的解,直到参考点R'的大小与新的外部档案AC'大小一致。
进一步的,步骤(5)中SMAA-2模型具体为:
首先,指标权重采用均匀分布来描述,如下所示:
指标值决策矩阵X=[xmn]M×N的计算公式为:
式中,xmn为方案m在属性n上的值;
假设决策矩阵X=[xmn]M×N,属性权重W={w1,w2,...,wN},SMAA-2通过线性效用函数对每个属性的效用值进行加权求和得到每个方案的综合效用um=u(xm,w);通过其加权值计算各方案的优劣排序,并从中选出满足决策要求的均衡方案;
式中,αr为二级权重,r越小,相应的二级权重越大,表明越看重排名靠前时的可接受度。
有益效果:与现有技术相比,本发明将多目标优化与多属性决策有机结合起来,提出了一个风光水电***多目标优化-决策的基本框架。首先考虑了传统的***效益和***出力的稳定性,同时还考虑了下游河道的生态健康,提出了一个风光水电***多目标优化调度模型;由于风光水电***设计电源多、拓扑关系复杂,直接求解难度较大,因此,本文将风光水电***进行解耦,解耦成风光被补偿子***和水电补偿子***;水电站子***约束较为复杂,因此在水电站子***构造决策变量交叉、变异的动态可行域,避免可行解在在进化过程中遭到破坏;随后得到完整Pareto前沿,应用能考虑决策者主观意愿和属性信息不确定性的随机多准则可接受性分析(SMAA-2)的多属性决策模型,给决策者提供了更多方案排序的信息,能更科学合理地进行决策。
附图说明
图1是本发明方法的流程图;
图2是单库动态可行域;
图3是库群动态可行域;
图4是雅砻江流域概化图;
图5是AR-MORA和NSGA-III的超体积指数随着迭代次数变化情况;
图6是风光水电***非劣解集;
图7是风光水电***非劣解集的二维投影图;
图8是桐子林的出流过程;
图9是方案A、B各库的水位过程;
图10是风光水电***Pareto曲面;
图11是不同风光水***调度方案的排名可接受指标;
图12是不同风光水***调度方案的中心权重;
图13是不同风光水***调度方案的全局可接受性指标。
具体实施方式
下面结合附图和实施案例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示,本发明提出的一种考虑下游生态的流域风光水***多目标优化调度方法,主要包括以下步骤:
(1)获取梯级水电站的基本信息,该基本信息包括各个水电站的入库径流资料、水位库容曲线、尾水流量关系曲线和水轮机机组综合曲线资料。
(2)获取流域风光电站的基本信息,该基本信息包括各个风光电站的地理位置、风速、太阳辐射、气温和装机容量等资料。
(3)根据流域风光水电站基本信息建立考虑下游生态的流域风光水***目标以及硬性约束的多目标发电优化调度数学模型;
根据流域风光水电站基本信息建立考虑水量平衡约束、水库上下限水位约束、流量约束、负荷约束和调度期末水位约束的多目标优化调度数学模型,其中目标函数为调度期内最小时段出力最大、发电量最大和对下游生态影响最小;
(31)发电量最大目标函数为:
其中,f1为流域风光水***总发电量;T为调度期时段数;PPVt、PWt和PHt分别是光伏电站、风力发电站和水电站第t时段出力;Δt为时段小时数;
最小时段出力最大目标函数为:
其中,Pt=PPVt+PWt+PHt;
水电站出库流量与天然流量偏差最小:
式中,O'i,t为第i水库t时段的出库流量,QNi,t为第i库t时段下游断面天然流量;I为水库数目。
(32)水量平衡以及硬性约束具体为:
(a)各库的水量平衡约束
Vi,t=Vi,t-1+(Qi,t-O'i,t)Δt(4);
(b)各库的上、下限水位约束
(c)流量约束
(d)调度期末水位约束
Zi,T=Zi,end,i∈[1,I](7);
其中,Zi,end为第i库调度期的期末水位。
(e)光伏电站出力约束
(f)风力发电站出力约束
(g)水电站出力约束
(h)外送断面约束
(33)风光水电***出力计算具体为:
(a)风力发电站出力计算:风力发电站出力计算方式如下:
式中,SA是风力发电机轮毂的面积,ρ是空气密度,Nk是第风力发电站的风力发电机的台数,uk,t是风力发电机轮毂处的风速,将气象站实测风速(10m)通过风速转换关系换算成轮毂处高程(80m高空)的风速,具体见下式:
(b)光伏电站处理计算
本实施例研究不具储能***的大型光伏发电***,建立光伏发电***与太阳辐射和温度的关系。
其中,Pstc是标准条件下(对应太阳辐射强度Gstc=1000W/m2,温度Tref=25℃)光伏板的出力,是第d光伏电站t时段的实际太阳辐射强度,β是光伏板的温度系数,γ是光伏板的辐射系数,是第d光伏电站光伏板的面积。
(c)水电站的出力:水电站的出力如下式所示:
PHi,t=g(Oi,t,ΔHi,t)(18);
其中,PHi,t是第i水电站第t时段出力,g(·)函数为水电站出力特性函数,Oi,t是发电流量,ΔHi,t是发电水头。
(4)采用改进的AR-MOEA对步骤(3)建立的多目标发电优化调度数学模型进行求解。
AR-MOEA(an adaptive reference point-based multi-objectiveevolutionary algorithm)算法是Tian于2017年提出的一种基于增强的反转欧式距离(IGD-NS)指标的多目标进化算法。AR-MOEA算法的基本框架与常用的基于性能指标的多目标进化算法框架类似。他们都是先随机初始化一个大小为N的种群P并将它复制到外部档案A,然后利用与基于目标分解的多目标进化算法相同的方法产生一个大小为NR的初始参考点集R。不同的地方在于AR-MORA算法基于IGD-NS指标的环境选择策略选出交配池P’,并且AR-MOEA算法采用了一组自适应的参考点作为计算IGD-NS指标的参考点集,这能够改善最终得到的种群的多样性和Pareto前沿的分布情况。
改进的AR-MOEA算法具体为:
(41)读取风速、太阳辐射、气温等资料,根据公式(15)-(17)计算风、光出力。
(42)随机生成M个初始个体。选择水库水位作为决策变量进行实数编码,具体的生成方式如下:
(43)将初始种群复制到外部档案AC,并且生成一个大小为NR的初始参考点集。
(44)将步骤(41)计算得到的风、光***的出力输入到梯级水电站***中,计算种群的适应度函数和IGD-NS贡献度,并采用锦标赛选择法选择贡献度大的个体到交配池。
IGD-NS指标具体为:
假设Q是种群P中无贡献解的集合。定义,即无贡献解表示对于任意的参考点集R,Q中的点离参考点集R的距离均大于待评价种群P与参考点R的距离。后半部分额外考虑了种群Q中无贡献解到参考点R中的距离。因此,在IGD值相同的情况下,含有更少且离R更近的无贡献解的种群能拥有更好的IGD-NS值。
(45)由于梯级水电站上、下游水库存在水量联系,各库运行策略相互影响,存在着多变量耦合约束,直接对个体进行交叉和变异不能保证子代个体满足约束条件,甚至可能破坏可行解。因此,本文在交叉、变异前基于梯级水电站复杂约束建立动态可行域。
对于单库来说,假设交叉、变异点为tp+1时刻,对于第j个约束,由于tp时刻水位固定,按照约束条件在tp时段不破坏的原则,可以确定tp+1时刻的正向可行域Forwardfeasible region(FFRj)(j=1~CountC,CountC是可能被破坏的约束条件个数);同理,tp+2时刻的水位固定,按照约束条件在tp+1时段不破坏原则进行反向计算,可以确定tp+1时刻的反向可行域Reverse Feasible Region(RFRj),此时对于第j个约束条件的可行域FRj是正向可行域和逆向可行域的交集:
FRj=FFRj∩RFRj(21);
满足所有约束的动态可行域为:
如图2所示。
1)出库流量约束
入库流量Qs,tp已知,根据水量平衡原理,Vs,tp+1的可行域变化范围:
合并入库流量tp+1时刻的入库流量正向可行域和反向可行域,得到tp+1入库流量可行域FR1(Vs,tp+1),如下所示:
FR1(Vs,tp+1)=FFR1(Vs,tp+1)∩RFR1(Vs,tp+1)(25);
2)出力约束
本发明采用耗水率函数描述出力,即:
gi(Oi,t,ΔHi,t)=min{Oi,t,OM(ΔHi,t)}/ξ(ΔHi,t)(26);
gi(Oi,t,ΔHi,t)=min{Oi,t,OM(ΔHi,t)}/ξ(ΔHi,t)(27);
其中,gi(Oi,t,ΔHi,t)为出力函数;OM(·)为机组在不同水头下的满发流量;ξ(·)为耗水率函数,根据水电站出力约束条件可知:
根据水量平衡得知,出力变化的可行域为:
合并tp+1时刻的出力正向可行域和反向可行域,得到tp+1出力可行域FR2(Vs,tp+1),如下所示:
FR2(Vs,tp+1)=FFR2(Vs,tp+1)∩RFR2(Vs,tp+1)(31);
由于梯级水库存在水量联系,上库放水策略可能会影响到下库的可行域。因此,采用轮库迭代的思想对梯级水电站群进行解耦,当对第s库求解动态可行域时,假定其他库水位不变,即将梯级水电站群的耦合约束转变成单库约束,如图3所示。
确定完动态可行域之后,在动态可行域内进行交叉和变异,具体如下所示:
改进的变异算子为:
(46)基于步骤(45)在动态可行域产生的子代个体,进行更新外部档案AC,然后利用更新后的外部档案AC进行自适应调整参考点R',使得参考点R'与Pareto前沿的形状相似。具体为:
通过选择交配池种群P'并经过交叉和变异得到子代后,利用新产生的解来更新外部档案AC和自适应参考点集合R',具体步骤如下:
1)基于当前种群P各个目标中的理想点和最差点N为目标数,外部档案AC、种群P和参考点R各个目标分别同时减去并同时乘以可以使得转换之后种群P、外部档案AC和参考点R标准化到相同的范围可以使得均匀分布的参考点R能够在不同的目标空间内生成均匀分布的方案。
2)基于标准化后的参考点R,计算外部档案AC的IGD-NS贡献值,删除外部档案AC里互相支配的和冗余的方案,并将剩余的方案复制到外部档案Acon。
3)将外部档案Acon复制到新的外部档案AC',基于新的外部档案AC'计算参考点R的IGD-NS值,并将有贡献的参考点保留,然后从新的外部档案AC'剩余的解选择最不拥挤的解,直到参考点R'的大小与新的外部档案AC'大小一致。
(47)合并父代与子代种群,并对其进行非支配排序,然后将基于IGD-NS指标的环境选择策略作用于第kM个前沿面的解上,直到k满足|Front1∪...∪FrontkM|≥M。
(48)判断是否达到最大迭代次数,若是,退出循环;若否,重复(44)到(47),直至达到最大迭代次数。
(5)采用SMAA-2对步骤(4)得到的非支配方案进行多属性决策。
使用多目标进化算法求解风光水电***多目标问题,得到Pareto最优解集,但是Pareto最优解集之间不存在优劣。在流域风光水电***复杂多变的决策环境下,决策属性信息伴随着决策者的主观意愿和自身的属性存在着一定的不确定性,决策者难以给出准确的偏好信息给决策者科学决策带来了困难。因此,本发明使用多属性决策方法,从Pareto非支配方案中进行评估和选择最优的解决方案,以供决策者使用。
SMAA-2模型具体为:
SMAA-2(Stochastic multicriteria acceptability analysis)模型在决策时不需要专家提前给定权重,而是通过反权重空间分析的决策流程,利用决策模型搜索在整个可行权重空间内随机分布的权重组合来揭示何种权重组合能够使某方案达到最优或某一排序,并计算在整个可行权重空间中的比例,并以此度量某方案达到最优或某一排序的概率,从而实现随机环境下的方案排序与优选。因此,本发明应用SMAA-2对上述风光水***的优化调度方案(风光水***发电效益、供电可靠程度和生态保障程度),考虑权重的不确定性条件下,进行随机多属性决策。指标权重采用均匀分布来描述,如下所示。
方案集合Α={Am|m=1,2,…,M}由上述改进的AR-MOEA算法求解得到,指标值决策矩阵为X=[xmn]M×N:
式中,xmn为方案m在属性n上的值:
假设决策矩阵X=[xmn]M×N,属性权重W={w1,w2,...,wN},SMAA-2通过线性效用函数对每个属性的效用值进行加权求和得到每个方案的综合效用um=u(xm,w)。通过其加权值计算各方案的优劣排序,并从中选出满足决策要求的均衡方案。
式中,αr为二级权重,r越小,相应的二级权重越大,表明越看重排名靠前时的可接受度。
下面以雅砻江流域风光水电站发电多目标优化调度为例,说明方法的有效性与合理性。
雅砻江流域位于中国青藏高原东部,具有丰富的风、光、水能资源,目前规划风、光、水电总装机容量6000万kW,其中水电规划装机容量2885万kW,风、光电规划装机容量3077万kW,是世界上目前规划的最大风光水互补清洁能源示范基地,具体概化图如图4所示。
表2显示了水电站的具体信息。此外,为了满足梯级水库群的实际需求,本次以1年为调度期,计算时段长为月,给定各水库与区间的来水过程和梯级负荷约束条件,各水库的约束参数见下表。
表1雅砻江风光水电站***装机容量
本次数值实验水电站径流数据采用的是水文站2016年实测数据,风速、太阳辐射和气温是在国家气象局网下载的逐日数据。经过多次试验,本实施例得到AR-MOEA求解风光水电***的最优参数。具体如下,种群规模取1000,交叉概率取1,变异概率取0.1,最大迭代次数为500。本次数值实验是在Matlab下运行,电脑配置是内存16GB,CPU是i7。分别采用AR-MOEA和NSGA-Ⅲ对雅砻江流域风光水电***进行求解。
为了进一步评估这两个多目标进化算法的性能,引入了超体积指标,超体积指标越大,表明算法得到的Pareto前沿的收敛性和分布性越好。图5为两个算法20次独立实验的超体积迭代过程的平均值。
从图5可以看出,在迭代初期,AR-MOEA的超体积指数大于NSGA-Ⅲ。虽然AR-MOEA和NSGA-Ⅲ算法的超体积指标均在迭代100次附近趋于稳定,但后者在优化结束前的超体积指标值较前者相比更少。在迭代400次之后,两种算法的超体积指标值均保持不变,这意味着Pareto前沿收敛到最终的状态。通过上述超体积指标随着迭代次数的变化情况可以看出,AR-MOEA算法在得到Pareto前沿方面表现要优于NSGA-Ⅲ,并且也不容易陷入局部最优解,因此,应用AR-MOEA算法求解效率会更高。
由于AR-MOEA算法在初始化具有随机性,为了消除随机性的影响,本发明采用改进的AR-MOEA算法独立地进行了20次实验,选取超体积指数最大的那次作为分析。下图显示了应用AR-MOEA算法求解风光水电***多目标问题得到的非劣解集,从图6可以看出,若要增大三个目标中一个目标的效益,则需要以牺牲其余两个目标作为代价,体现出风光水***发电效益、发电稳定性和下游生态三者之间的竞争博弈关系。
为了进一步分析三个目标之间的关系,对三维坐标下的点距进行二维投影,获得目标的投影,如图7所示。
由图7a可知,当保证出力固定时,随着发电量的增大,APFD值也随之增大,表现出生态效益随着发电量的增大而减小。当保证出力减小时,发电量和生态效益两个目标之间的关系右移,即随着保证出力减小,表现出生态效益和发电量分别呈现出增加的趋势。由图7b可以看出,当发电量固定时,APFD随着最小出力的增加而增大,即生态效益随着出力稳定性的增大而减小。由图7c可以看出,当APFD值较小时,即生态效益较优,表现出发电量与保证出力的竞争关系不太明显,随着APFD值增大,发电量与保证出力呈现出一个正比的关系,即保证出力增大,发电量随之增大。
由图6可以得出:在风光水电***多目标优化调度中,三个目标之间存在着竞争关系,其中发电量目标和生态目标的竞争关系最强,其次是生态效益目标和出力稳定性目标的竞争关系。
为了再进一步分析风光水***多目标优化调度过程特征,选择A和B两个方案,其保证出力均为10011.87万kW,但是方案A发电量比方案B多了198.94亿kWh。取出方案A、B和天然情况下在下游断面的出库流量,如图8所示,方案A、B两河口、锦屏一级和二滩的水位过程,如图9所示。
由图8和图9可知,A和B方案桐子林下游断面的出流过程在蓄水期相差不大,这主要是因为汛期过后,水库需要在年末回蓄正常高水位,水库在蓄水期尽量蓄水,水库的水位过程相似,因此,水库的出流过程相差不大。而在消落期和汛期这两个阶段,A方案较B方案弃水多了472.83亿m3/s,这是因为A方案为了出库流量尽可能的与天然流量过程一致,前期尽量减少出库流量,导致后期弃水较多,因此发电量较少。而B方案为了尽可能的增加风光水***的发电量,在前期增加了出库流量,为汛期的到来腾空了库容,减少了弃水,从而增加了发电量。
由上述分析可知,风光水***发电效益、出力稳定性和下游河道生态效益之间存在激烈的竞争关系。同时,不同决策者对于不同目标的偏好不一样,且不同目标的自身属性也不一样,因此存在着一定的不确定性。因此,本发明采用第3节建立的随机多属性决策模型进行决策。从上述Pareto曲面均匀选出12个方案,如图10所示。
由于上述决策问题包含12个方案和3个指标,若直接进行求解涉及36维积分,问题的复杂度较大,因此,本发明采用蒙特卡罗的方法对上述随机多属性决策模型进行求解,为了保证计算精度,随机模拟次数取10000次。得到的排名可接受指标如图11所示。
图11直观地展示了排序可接受性指标的三维柱状图。由图11可知,上述12个方案均具有一定概率获得排序第一,其中,方案2获得排序第一具有最大的可接受指标,但其也具有一定概率获得排序中间和靠后。方案4获得排序第一具有第二大的可能性,但是其获得排序最后一名的可接受性指标也很大。方案12具有最大的概率获得排序最后,但是其也有一定概率获得排序靠前和居中。
图12给出了各个方案排序最优时的中心权重向量,各个方案排序最优时的中心权重向量均不一样。具体地,方案4较其他方案发电量效益占优,下游河道生态效益和出力稳定性均呈现弱势,由于发电效益和出力稳定性是正向指标,生态效益是负向指标,因此对于其获得最优排名时的中心权重向量,赋予了发电效益这个指标较大权重,而下游河道生态和出力稳定性较小指标。对于方案12,发电量和出力稳定性较差,而下游河道生态最优,因此,赋予了下游河道生态较大权重,发电量和出力稳定性赋予了较小的权重。
为了从整体上衡量各个方案的综合排序,采用公式(40)计算全局可接受性指标,结果如图13所示。其中,二级权重αr采用重心形式,排序越靠前分配的权重越大,由图10可知,方案2、3、4和1排序靠前的概率较大,因此,方案2、3、4和1的全局可接受性指标较大。从排序结果来看,这种排序结果是合理的。
Claims (1)
1.一种考虑下游生态的流域风光水***多目标优化调度方法,其特征在于,包括以下步骤:
(1)获取梯级水电站的基本信息,该基本信息包括各个水电站的入库径流资料、水位库容曲线、尾水流量关系曲线和水轮机机组综合曲线资料;
(2)获取流域风光电站的基本信息,该基本信息包括各个风光电站的地理位置、风速、太阳辐射、气温和装机容量资料;
(3)根据流域风光水电站基本信息,建立考虑下游生态的流域风光水***目标以及硬性约束的多目标发电优化调度数学模型;具体为:
多目标优化调度数学模型的目标函数为发电量、最大调度期内最小时段出力最大和对下游生态影响最小;
发电量最大目标函数为:
其中,f1为风光水电***总发电量;T为调度期时段数;PPVt、PWt和PHt分别是光伏电站、风力发电站和水电站第t时段出力;△t为时段小时数;
最小时段出力最大目标函数为:
其中,f2为风光水电***最小出力;Pt=PPVt+PWt+PHt;
水电站出库流量与天然流量偏差最小目标函数为:
其中,f3为水电站的出库流量与天然流量的偏差;O'i,t为第i水库t时段的出库流量,QNi,t为第i库t时段下游断面天然流量;I为水库数目;
硬性约束具体为:
(a)各库的水量平衡约束
Vi,t=Vi,t-1+(Qi,t-O'i,t)△t ;
(b)各库的上、下限水位约束
(c)流量约束
(d)调度期末水位约束
Zi,T=Zi,end,i∈[1,I];
其中,Zi,end为第i库调度期的期末水位;
(e)光伏电站出力约束
(f)风力发电站出力约束
(g)水电站出力约束
(h)外送断面约束
风光水电***出力计算具体为:
(a)风力发电站出力计算方式如下:
其中,SA是风力发电机轮毂的面积,ρ是空气密度,Nk是第风力发电站的风力发电机的台数,uk,t是风力发电机轮毂处的风速,具体见下式:
(b)光伏电站出力计算:建立光伏发电***与太阳辐射和温度的关系;
其中,Pstc是标准条件下光伏板的出力,Gstc是标准条件下的太阳辐射强度,是第d光伏电站第t时段的光伏板的温度,Tref是标准条件下的温度,是第d光伏电站t时段的实际太阳辐射强度,β是光伏板的温度系数,γ是光伏板的辐射系数,是第d光伏电站光伏板的面积;
(c)水电站的出力如下式所示:
PHi,t=g(Oi,t,△Hi,t);
(4)采用改进的AR-MOEA对所述步骤(3)建立的多目标发电优化调度数学模型进行求解;改进的AR-MOEA算法具体为:
(41)读取风速、太阳辐射和气温资料,根据风力发电站出力计算公式和光伏电站出力计算公式计算风、光出力;
(42)随机生成M个初始个体;选择水库水位作为决策变量进行实数编码,具体的生成方式如下:
(43)将初始种群复制到外部档案AC,并且生成一个大小为NR的初始参考点集;
(44)将步骤(41)计算得到的风、光***的出力输入到梯级水电站***中,计算种群的适应度函数和IGD-NS贡献度,并采用锦标赛选择法选择贡献度大的个体到交配池,IGD-NS指标具体为:
(45)由于梯级水电站上、下游水库存在水量联系,各库运行策略相互影响,存在着多变量耦合约束,直接对个体进行交叉和变异不能保证子代个体满足约束条件;因此,在交叉、变异前基于梯级水电站复杂约束建立动态可行域,具体为:
满足所有约束的动态可行域为:
其中,FRj为第j个约束条件的可行域,CountC是可能被破坏的约束条件个数;
确定完动态可行域之后,在动态可行域内进行交叉和变异,具体如下所示:
改进的变异算子为:
(46)基于步骤(45)在动态可行域产生的子代个体,进行更新外部档案AC,然后利用更新后的外部档案AC进行自适应调整参考点R',使得参考点R'与Pareto前沿的形状相似;具体为:
通过选择交配池种群P'并经过交叉和变异得到子代后,利用新产生的解来更新外部档案AC和自适应参考点集合R',具体步骤如下:
1)基于当前种群P各个目标中的理想点和最差点N为目标数,外部档案AC、种群P和参考点R各个目标分别同时减去并同时乘以使得转换之后种群P、外部档案AC和参考点R标准化到相同的范围使得均匀分布的参考点R能够在不同的目标空间内生成均匀分布的方案;
2)基于标准化后的参考点R,计算外部档案AC的IGD-NS贡献值,删除外部档案AC里互相支配的和冗余的方案,并将剩余的方案复制到外部档案Acon;
3)将外部档案Acon复制到新的外部档案AC',基于新的外部档案AC'计算参考点R的IGD-NS值,并将有贡献的参考点保留,然后从新的外部档案AC'剩余的解选择最不拥挤的解,直到参考点R'的大小与新的外部档案AC'大小一致;
(47)合并父代与子代种群,并对其进行非支配排序,然后将基于IGD-NS指标的环境选择策略作用于第kM个前沿面的解上,直到k满足|Front1∪...∪FrontkM|≥M;
(48)判断是否达到最大迭代次数,若是,退出循环;若否,重复步骤(44)到(47),直至达到最大迭代次数;
(5)采用SMAA-2对步骤(4)得到的非支配方案进行多属性决策;
SMAA-2模型具体为:
首先,指标权重采用均匀分布来描述,如下所示:
指标值决策矩阵X=[xmn]M×N的计算公式为:
式中,xmn为方案m在属性n上的值;
假设决策矩阵X=[xmn]M×N,属性权重W={w1,w2,…,wN},SMAA-2通过线性效用函数对每个属性的效用值进行加权求和得到每个方案的综合效用um=u(xm,w);通过其加权值计算各方案的优劣排序,并从中选出满足决策要求的均衡方案;
式中,αr为二级权重,r越小,相应的二级权重越大,表明越看重排名靠前时的可接受度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910067807.XA CN109886473B (zh) | 2019-01-24 | 2019-01-24 | 一种考虑下游生态的流域风光水***多目标优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910067807.XA CN109886473B (zh) | 2019-01-24 | 2019-01-24 | 一种考虑下游生态的流域风光水***多目标优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109886473A CN109886473A (zh) | 2019-06-14 |
CN109886473B true CN109886473B (zh) | 2020-05-05 |
Family
ID=66926735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910067807.XA Expired - Fee Related CN109886473B (zh) | 2019-01-24 | 2019-01-24 | 一种考虑下游生态的流域风光水***多目标优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109886473B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110880048B (zh) * | 2019-11-06 | 2022-06-21 | 国网湖北省电力有限公司宜昌供电公司 | 一种梯级水库生态随机优化调度模型及求解方法 |
CN110991703B (zh) * | 2019-11-14 | 2024-06-07 | 华中科技大学 | 一种基于改进的nsga-ii算法的风光水联合发电调度方法 |
CN111342486B (zh) * | 2020-02-24 | 2022-10-14 | 河海大学 | 一种含梯级水电的风光水互补发电***的优化调度方法 |
CN111626512B (zh) * | 2020-05-28 | 2023-11-14 | 国家电网公司西南分部 | 一种考虑电网断面约束的梯级水电日前调度方法 |
CN111898253B (zh) * | 2020-07-15 | 2022-04-15 | 武汉大学 | 水库调度及其下游河流生态环境保护的合作价值评估方法 |
CN112183839B (zh) * | 2020-09-22 | 2022-08-30 | 河海大学 | 一种含梯级水电风光水互补发电***中长期优化调度方法 |
CN112184070B (zh) * | 2020-10-28 | 2021-12-17 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种协同生态流量需求的梯级水电站多目标优化调度方法及*** |
CN113205250B (zh) * | 2021-04-28 | 2022-08-19 | 河海大学 | 一种考虑水库调度风险的风光水互补发电计划编制方法 |
CN113610638B (zh) * | 2021-08-20 | 2023-11-21 | 东北大学秦皇岛分校 | 基于smaa-ds的信用等级与违约损失率相匹配的评级***及方法 |
CN113572168B (zh) * | 2021-09-26 | 2022-02-08 | 江西江投电力技术与试验研究有限公司 | 一种考虑风光水能源多重不确定性风险的优化调度方法 |
CN114154404B (zh) * | 2021-11-22 | 2024-05-21 | 大连理工大学 | 一种利用观测数据推断相邻水电站运行状态及参数的方法 |
CN115809718B (zh) * | 2022-10-09 | 2024-03-15 | 当涂南科院科研试验有限公司 | 基于多目标竞争关系量化的梯级电站发电与生态协同优化方法和*** |
CN115765044B (zh) * | 2022-11-26 | 2023-05-09 | 水利部水利水电规划设计总院 | 风光水电***联合运行及风险分析方法和*** |
CN115600822B (zh) * | 2022-11-28 | 2023-02-28 | 水利部水利水电规划设计总院 | 跨省江河流域初始水权分配方法和*** |
CN116415742B (zh) * | 2023-06-12 | 2023-08-22 | 水利部水利水电规划设计总院 | 风光水电***两阶段随机多属性决策方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104538992A (zh) * | 2015-01-14 | 2015-04-22 | 云南电力调度控制中心 | 一种大小水电及风电协调优化调度方法 |
CN105048516A (zh) * | 2015-08-18 | 2015-11-11 | 四川大学 | 一种风光水火多源互补优化调度方法 |
CN107506873A (zh) * | 2017-09-21 | 2017-12-22 | 武汉大学 | 基于调节性能的水风光电站群协同组合划分方法及*** |
CN107527119A (zh) * | 2017-09-06 | 2017-12-29 | 河海大学 | 基于改进的多目标量子遗传算法的水资源优化调度方法 |
CN108805434A (zh) * | 2018-05-25 | 2018-11-13 | 河海大学 | 一种基于改进nsga-ⅲ的梯级水电站多目标优化调度方法 |
-
2019
- 2019-01-24 CN CN201910067807.XA patent/CN109886473B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104538992A (zh) * | 2015-01-14 | 2015-04-22 | 云南电力调度控制中心 | 一种大小水电及风电协调优化调度方法 |
CN105048516A (zh) * | 2015-08-18 | 2015-11-11 | 四川大学 | 一种风光水火多源互补优化调度方法 |
CN107527119A (zh) * | 2017-09-06 | 2017-12-29 | 河海大学 | 基于改进的多目标量子遗传算法的水资源优化调度方法 |
CN107506873A (zh) * | 2017-09-21 | 2017-12-22 | 武汉大学 | 基于调节性能的水风光电站群协同组合划分方法及*** |
CN108805434A (zh) * | 2018-05-25 | 2018-11-13 | 河海大学 | 一种基于改进nsga-ⅲ的梯级水电站多目标优化调度方法 |
Non-Patent Citations (1)
Title |
---|
阿里地区能源利用方案及多能互补***研究;万久春;《中国优秀博硕士学位论文全文数据库(硕士)社会科学I辑(经济政治与法律)》;20040315(第1期);第4-5章 * |
Also Published As
Publication number | Publication date |
---|---|
CN109886473A (zh) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109886473B (zh) | 一种考虑下游生态的流域风光水***多目标优化调度方法 | |
Motevasel et al. | Expert energy management of a micro-grid considering wind energy uncertainty | |
Moghaddam et al. | Multi-objective operation management of a renewable MG (micro-grid) with back-up micro-turbine/fuel cell/battery hybrid power source | |
Liu et al. | Multi-objective optimization scheduling of wind–photovoltaic–hydropower systems considering riverine ecosystem | |
CN112508221A (zh) | 一种有限储能下考虑源-荷不确定性的日前调度决策方法 | |
CN108205717A (zh) | 一种光伏发电功率多时间尺度预测方法 | |
CN104484833A (zh) | 基于遗传学算法改进的rbf-bp神经网络的光伏发电输出功率追踪算法 | |
Li et al. | Hierarchical multi-reservoir optimization modeling for real-world complexity with application to the Three Gorges system | |
CN108320062A (zh) | 一种基于多目标种群群搜索算法的联合调度方法及*** | |
CN103440370A (zh) | 一种输变电工程造价评估方法及装置 | |
CN111786420A (zh) | 一种考虑再生能源混合发电的电力***经济调度的方法 | |
Bandopadhyay et al. | Application of hybrid multi-objective moth flame optimization technique for optimal performance of hybrid micro-grid system | |
CN103580020A (zh) | 一种基于NSGA-II和Look-ahead的含风电场电力***多目标动态优化调度方法 | |
CN110264012A (zh) | 基于经验模态分解的可再生能源功率组合预测方法及*** | |
CN116667325B (zh) | 一种基于改进布谷鸟算法的微电网并网运行优化调度方法 | |
CN111724064B (zh) | 一种基于改进免疫算法的含储能配电网规划方法 | |
CN117578409A (zh) | 一种电力市场环境下的多能互补优化调度方法及*** | |
CN116757446A (zh) | 基于改进粒子群算法的梯级水电站调度方法及*** | |
Kyomugisha et al. | Performance of various voltage stability indices in a stochastic multiobjective optimal power flow using mayfly algorithm | |
CN112668186B (zh) | 基于elm的输配一体储能***选址定容协同优化方法 | |
CN113807566A (zh) | 基于可再生能源和内部负荷不确定的社区能源的调度方法 | |
Su et al. | Optimal placement and capacity sizing of energy storage systems via NSGA-II in active distribution network | |
Pandya et al. | Stochastic wind-thermal power plants integrated multi-objective optimal power flow | |
Tiwari et al. | Economic dispatch in renewable energy based microgrid using manta ray foraging optimization | |
CN114580298A (zh) | 一种基于rbf-isbo的微电网优化调度方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200505 |