CN108667077A - 一种风蓄联合***优化调度方法 - Google Patents
一种风蓄联合***优化调度方法 Download PDFInfo
- Publication number
- CN108667077A CN108667077A CN201810313720.1A CN201810313720A CN108667077A CN 108667077 A CN108667077 A CN 108667077A CN 201810313720 A CN201810313720 A CN 201810313720A CN 108667077 A CN108667077 A CN 108667077A
- Authority
- CN
- China
- Prior art keywords
- wind
- power
- formula
- hours
- value
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种风蓄联合***优化调度研究方法,其内容包括:定义抽水蓄能机组运行状态,建立抽水蓄能状态模型;考虑风电机组出力特性,对风力发电机组进行建模;基于可信性理论求解机会约束规划;建立基于可信性理论的机会约束规划的风蓄联合优化调度模型;结合改进的粒子群算法和随机模拟技术对所述风蓄联合优化调度模型进行求解。本发明提出的优化策略相比现有技术更加简捷方便、实用性高,而且可以更有效地兼顾减少风电的间歇性和波动性、提高风电消纳量和增加风电上网经济效益的效果。
Description
技术领域
本发明涉及新能源发电运行领域,尤其涉及一种风蓄联合***优化调度方法。
背景技术
风能作为一种可再生绿色能源得到人们的广泛重视,随着风力发电技术迅猛发展,其发电成本已得到大幅度下降,具有较好的社会经济意义和可持续发展意义,但是受到风速不确定因素的影响,风力发电具有很强的间歇性和随机波动性。不仅导致大量的风能成为弃风,大大降低了风电的消纳;而且大规模风电的接入会大大增加电力***的不确定性,对电力***运行过程中的优化调度、***稳定性和电能质量等方面造成影响。针对以上问题,采取合理的方式平抑风电功率输出,减少弃风量,实现电力***的低碳经济调度具有重要的社会意义。
近年来,随着我国越来越多的大中型风电场相继投入运行,为保证供电的均衡和连续,储能装置已经成为风力发电***的关键配套部件。其中,最成熟也是世界应用最普遍的物理储能--抽水储能,作为存储电能的一种方式,它的优点不仅是释放时间可以从几小时到几天,且能量转换率高(70%-85%),抽水储能可以间接储存电能,在风能过剩的时候利用多余的风能驱动水泵,将水从下水库抽到上水库储存起来,然后在风能缺乏的时候将水放出发电,并流入下水库.
现阶段,针对风电和抽水蓄能联合运行优化问题,国内外专家学者进行了一些探索性的工作,这些研究多集中于风电和抽水蓄能联合运行策略、调度模型优化求解以及***中风电预测出力不确定性问题的处理等方面。虽然都针对风电的消纳和经济效益等问题进行了处理,但整体考虑的都不够全面。
综上所述,有必要发明一种风蓄联合***优化调度研究策略,即增加了风电的消纳、又增加了风电的上网经济效益、同时使得协调后的整体的风力发电间歇和波动性相比协调前有了明显减少,使整个***的运行较协调前平稳了很多。
发明内容
本发明目的在于提供一种风蓄联合***优化调度方法,该方法具有简便、考虑因素全面、以及实用性高的特点。
为实现上述目的,本发明是通过以下技术方案实现的:
一种风蓄联合***优化调度方法,所述方法内容包括如下步骤:
步骤1定义抽水蓄能机组运行状态,建立抽水蓄能状态模型;
步骤2考虑风电机组出力特性,对风力发电机组进行建模;
步骤3基于可信性理论求解机会约束规划;
步骤4建立基于可信性理论的机会约束规划的风蓄联合优化调度模型;
步骤5结合改进的粒子群算法和随机模拟技术对步骤4的模型进行求解。
进一步的,在步骤1中,所述定义抽水蓄能机组运行状态,建立抽水蓄能状态模型,其具体过程如下:
定义抽水蓄能机组运行状态,建立抽水蓄能机组状态模型;
抽水蓄能机组与风电机组联合后组成新的优化调度,可以充分发挥抽水蓄能机组的优势,对风电机组出力进行有效平滑;本发明方法是根据日前风电出力预测值来确定抽水蓄能机组的运行状态,但是由于风电预测出力存在误差,因此,根据实时的风电出力变化,建立风蓄联合协调出力的新型节能调度方式;
本发明方法引入0-1整型变量来考虑机组的抽水和发电两个阶段,建立抽水蓄能机组的状态模型如下:
X(t,k)+Y(t,k)=1 (1)
(1)式中:X(t,k)=1时表示机组在第t小时第k(协调周期为24小时,每小时分四个时段,k=1,2,3,4)时段为发电阶段,Y(t,k)=1时机组为抽水阶段。
进一步的,在步骤2中,所述考虑风电机组出力特性,其具体过程如下:
1.单台风电机出力特性
风速为服从Weibull分布的随机变量,风电机组的功率特性曲线一般由风机制造商提供,也可通过实测得到。在计算中,根据风电机组的功率特性曲线,单台风电机的输出功率PIF与风速关系v可近似用分段函数表示为:
(2)式中:vr和PWr分别表示风机在第t小时第k时段的额定风速和额定功率;vin为切入风速;vout为切出风速;
2.描述风电机组总出力曲线特性
总的风电机组的发电功率PW可表示为:
Pw(t,k)=efNW(t,k)PIF(t,k) (3)
0≤PIF(t,k)≤PWR (4)
(3)和(4)式中:PIF(t,k),Pw(t,k),分别为第t小时第k时段单台风机和总的风电机组发电功率(MW);e,f分别为风机的传动效率和发电效率(%);NW(t,k)为第t小时第k时段风电场中正常运行的风机台数;PWR(t,k)为每台机组的额定功率(MW)。
进一步的,在步骤3中,所述基于可信性理论求解机会约束规划,其具体求解过程如下:
机会约束规划是随机规划理论的重要分支,由Charnes和Cooper提出,目前,有两种方法用来求解机会约束规划,第一种是利用清晰类等价算法将其转化确定性模型;本发明方法采用第二种方法即通过随机模拟技术来求解机会约束规划:
机会约束规划是一种用于解决在给定置信水平下含有随机变量的优化问题的随机规划方法,且考虑到决策变量在随机变量满足约束条件后才能作出决策。本发明方法提出如下模型来求解机会约束规划,具体分析如式(1)~式(3)所示。
Cr{gj(x,ξ)≤0,j=1,2,…,p}≥β (7)
在(5)、(6)和(7)式中:f(x,ξ)为目标函数;为目标函数的乐观值;gj(x,ξ)为约束函数;x为决策向量;ξ为一个随机参数向量;α、β分别为决策者预先给定的置信水平;Cr{·}为{·}中事件的可信性。
可信性测度是描述机会约束规划不确定性情况的一种重要工具。假设(Θ,P(Θ),Pos)为可能性空间,A为幂集P(Θ)中的一个元素,则其可信性测度可表示为:
(8)式中:Pos{A}为事件A的可能性测度;Ac为A的对立事件。
进一步的,在步骤4中,所述建立基于可信性理论的机会约束规划的风蓄联合优化调度模型,其具体实现过程如下:
确定目标函数及约束条件
本发明方法考虑了风电场的上网经济和抽水蓄能机组抽水的花费,通过约束风电场输出有功功率使其趋于平滑,同时满足风电场的功率输出限制,最大程度地利用风能,以风电***的上网经济效益最大为目标。
此外,由于风电场出力的随机性,导致风蓄联合出力的取值也具有随机性,因此整个周期内通过抽水机组获取的抽水功率和发电功率均具有随机性。若仍采用传统经济调度约束条件中等式或不等式的形式,则结果将过于保守,且没有实际意义,因此本发明方法以概率的形式来描述约束条件,即约束条件含有置信水平,而置信水平的取值来自于***实际的要求。因此,基于可信性理论的机会约束规划的风蓄联合运行优化调度的数学模型如下:
风蓄联合运行优化目标函数:
X(t,k)+Y(t,k)=1 (14)
式中:为风电上网经济效益;Pws(t,k)为第t小时第k时段风蓄联合出力值;Pt Y,Pt X分别为风蓄联合***在抽水状态和发电状态下的联合出力值;Pw(t,k)为第t小时第k时段风电机组实际出力;Pw.t为第t小时的风电机组实际出力;Pp(t,k)和Pg(t,k)分别为机组在第t小时第k时段的抽水功率和发电功率;C(t,k)为第t小时第k时刻的上网电价;Ci(t,k)为第t小时第k时刻的抽水电价;T为一个周期即24小时;K为每小时内时段数即4个时段。
约束条件:
①风蓄联合***确保小时内各个时段的协调总出力保持稳定,满足风电场输出有功功率相对平滑:
Pws(t,k)=Pws(t,k+1) (17)
②目标利润的约束条件:
式中:α0为满足目标利润最大的置信水平。
③风蓄联合***出力在抽水状态下约束:
式中:φ1为风蓄联合机组在第t小时第k时刻抽水状态的联合出力约束置信水平;表示机组在第t小时第k时段风电机组实际出力最小值。
④风蓄联合***出力在发电状态下约束:
式中:φ2为风蓄联合机组在第t小时第k时刻的发电状态下的联合出力约束置信水平;表示机组在第t小时第k时段风电机组实际出力最大值;Pmax(t,k)为第t小时第k时段风电场输送到电网的功率最大值。
⑤风蓄联合***出力约束:
Pmin(t,k)≤Pws(t,k)≤Pmax(t,k) (21)
式中:Pmin(t,k)、Pmax(t,k)为第t小时第k时段风电场输送到电网的功率最小值和最大值;
⑥风蓄联合***抽水功率和发电功率约束:
式中:ηg为水泵发电效率;其余见上文解释。
⑦风电机组实际出力约束:
Pdmin(t,k)≤Pw(t,k)≤Pdmax(t,k) (23)
式中:Pdmin(t,k)、Pdmax(t,k)为第t小时第k时段风电机组装机容量的最大值和最小值。
⑧抽水蓄能电站容量约束:
0≤Q(t,k)≤PpgR (27)
0≤Pp(t,k)≤Pgmax(t,k) (29)
式中:PpgR为机组额定容量;c为机组循环效率;Q(t,k)为第t小时第k时段同步机抽水功率(MW);Ppmax(t,k)和Pgmax(t,k)分别为抽水功率和发电功率上限;E(t,k)为抽水蓄能电站存储量;相邻两个时段的间隔为ΔT。
⑨抽水蓄能电站水库约束:
本发明方法将上下水库看成一个整体,并将水库储水量转化为存储电量,不考虑水位对***的影响,本发明方法将存储电量和水量之间看作是线性关系,储水量受到上下限约束,且在***运行一个调度周期后,储水量未发生变化:
Emin(t,k)≤E(t,k)≤Emax(t,k)= (30)
E(T,k)=T(t,k) (31)
式中:Emin(t,k)、Emax(t,k)分别表示第t小时第k时刻风蓄联合***储水量的最大和最小值。
⑩在调度周期内,风蓄联合***的储水量变化公式满足:
式中:ηp为水泵抽水效率;ηg为水泵发电效率;ηp取0.85;ηg取0.8。
进一步的,在步骤5中,所述结合改进的粒子群算法和随机模拟技术对模型进行求解,其具体实现过程如下:
粒子群算法是Kennedy博士和Eberhart教授提出的一种智能寻优方法。其对求解含有大量约束条件的非线性规划问题具有很快的收敛速度及较好的全局寻优能力。粒子群基本公式为:
式中:m表示粒子的编号;t分别表示算法运行的迭代步数;d表示搜索空间维度D的第d维分量;h表示全局最优;Pm为粒子的个体极值,xmd,vmd,pmd,vmax.d,vmin.d,hd分别表示Xm,Vm,Pm,vmax,vmin,h的第d维分量;c1,c2通常取值为2;r1,r2是介于[0,1]之间的随机数。
由于粒子群算法中惯性权重对寻优搜索影响较大,本发明方法参考基于自适应度调整粒子群算法思想对其进行了改进。搜索初始阶段,设定惯性权重取较大值,可以提高寻优搜索的效率,在搜索中后期,随着搜寻结果的不断优化,采用改进后的惯性权重取值公式使惯性权重随着搜索步骤的进行不断减小,从而进一步提高寻优搜索的准确性,获取更多的经济效益。改进后的惯性权重ω可表示为:
式中:ωmax为初始权重;ωmin为最终权重;nmax为最大迭代次数;n0为当前迭代次数。ωmax取0.8。
随机模拟方法,也称为蒙特·卡罗(Monte Carlo)方法。它的基本思想是:首先建立一个概率模型或随机过程求各方面的问题,然后对模型或过程进行分析计算,最后给出所求解的近似值。本发明方法将随机模拟技术结合改进后的粒子群算法来求解机会约束规划模型,提高寻优搜索的精度和收敛速度。具体求解思路如下:
(1)输入***数据,置信水平以及算法控制参数;
(2)以各个小时内机组在抽水状态下和发电状态下的联合出力作为随机变量,在可行域范围内由概率密度函数生成N个随机变量ξN即产生一组初始方案作为粒子群算法的初始种群,并初始化粒子群位置和速度;
(3)任意给定决策变量x,设其初始值m=0,将x和ξ值代入式(3),若满足使gj(x,ξ)≤0成立,则令m=m+1,当N取值较大,且满足m/N≥α时,式(3)成立,若不满足,返回步骤(2);
(4)若满足,此时设初始值为n,判断N次随机序列中若n/N≥β,则开始排序置fi=f(x,ξi),其中i=1,2,…,N,得到序列{f1,f2,…,fN},置K为βiN的整数部分,则序列中第K个最大元素即为目标函数的最优值;
(5)结合步骤(4)及适应度优劣程度确定粒子群算法的个体极值和全局极值;
(6)更新粒子速度和位置,利用式(34)计算惯性权重,并采用式(32)、和式(33)对种群进行迭代计算,得到新一代的种群;
(7)重复步骤(3)至(6),当满足最大允许迭代次数时,停止搜索,并输出结果即为风蓄联合机组在抽水状态下和发电状态下的最优联合出力值;否则返回步骤(3)继续搜索。
由于采用上述技术方案,本发明提供的与现有技术相比具有这样的有益效果:
本发明提出的优化策略相比现有技术更加简捷方便、实用性高,而且可以更有效地兼顾减少风电的间歇性和波动性、提高风电消纳量和增加风电上网经济效益的效果。且本发明提出的将改进惯性权重后的粒子群算法和随机模拟技术结合对模型求解,提高算法的寻优精度和收敛速度。
附图说明
图1是本发明方法***结构图;
图2是本发明方法中风能与抽水蓄能的协调方式图;
图3是本发明方法的总体流程图。
具体实施方式
本发明方法是在考虑大规模分布式电源接入主动配电网情况下提出的。
下面结合图1、图2和图3,详细说明本发明所述的一种风蓄联合***优化调度方法,如图3所示,该方法具体步骤如下:
(1)定义抽水蓄能机组运行状态,建立抽水蓄能状态模型
抽水蓄能机组与风电机组联合后组成新的优化调度,可以充分发挥抽水蓄能机组的优势,对风电机组出力进行有效平滑。本发明方法是根据日前风电出力预测值来确定抽水蓄能机组的运行状态,但是由于风电预测出力存在误差,因此,根据实时地风电出力变化,建立风蓄联合协调出力的新型节能调度方式。
本发明方法引入0-1整型变量来考虑机组的抽水和发电两个阶段,建立抽水蓄能机组的状态模型如下:
X(t,k)+Y(t,k)=1 (1)
式中:X(t,k)=1时表示机组在第t小时第k(协调周期为24小时,每小时分四个时段,k=1,2,3,4)时段为发电阶段,Y(t,k)=1时机组为抽水阶段。
(2)考虑风电机组出力特性,对风力发电机组进行建模
描述风电机组出力特性;
(2-1)单台风电机出力特性
风速为服从Weibull分布的随机变量,风电机组的功率特性曲线一般由风机制造商提供,也可通过实测得到。在计算中,根据风电机组的功率特性曲线,单台风电机的输出功率PIF与风速关系v可近似用分段函数表示:
(2)式中:vr和PWr分别表示是风机在第t小时第k时段的额定风速和额定功率;vin为切入风速;vout为切出风速;
(2-2)描述风电机组总出力特性
总的风电机组的发电功率PW可表示为:
Pw(t,k)=efNW(t,k)PIF(t,k) (3)
0≤PIF(t,k)≤PWR (4)
(3)和(4)式中:PIF(t,k),Pw(t,k),分别为第t小时第k时段单台风机和总的风电机组发电功率(MW);e,f分别为风机的传动效率和发电效率(%);NW(t,k)为第t小时第k时段风电场中正常运行的风机台数;PWR(t,k)为每台机组的额定功率(MW)。
(3)基于可信性理论求解机会约束规划
机会约束规划是随机规划理论的重要分支,由Charnes和Cooper提出,目前,有两种方法来求解机会约束规划,第一种是利用清晰类等价算法将其转化确定性模型;本发明方法采用第二种方法即通过随机模拟技术来求解。
机会约束规划是一种用于解决在给定置信水平下含有随机变量的优化问题的随机规划方法,且考虑到决策变量在随机变量满足约束条件后才能作出决策。本发明方法提出如下模型来求解机会约束规划,具体分析如式(1)~式(3)所示。
Cr{gj(x,ξ)≤0,j=1,2,…,p}≥β (7)
在(5)、(6)和(7)式中:f(x,ξ)为目标函数;为目标函数的乐观值;gj(x,ξ)为约束函数;x为决策向量;ξ为一个随机参数向量;α、β分别为决策者预先给定的置信水平;Cr{·}为{·}中事件的可信性。
可信性测度是描述机会约束规划不确定性情况的一种重要工具。假设(Θ,P(Θ),Pos)为可能性空间,A为幂集P(Θ)中的一个元素,则其可信性测度可表示为:
(8)式中:Pos{A}为事件A的可能性测度;Ac为A的对立事件。
(4)建立基于可信性理论和机会约束规划的风蓄联合优化调度模型
(4-1)风蓄联合***结构
我国响应大力开发可再生能源的号召,以及跟据国际上提出的节能减排的要求,本发明方法建立了包含风电机组与抽水蓄能机组协调的风蓄联合机组、常规火电机组的优化调度模型,并以优先调度风电为原则,尽可能地多利用风电,其***结构如图1所示。
(4-2)风蓄联合***的协调方式
以一天24小时为周期,当某一小时内的总的风力发电值大于周期内平均每小时发电值时,则该小时被定义为抽水蓄能机组的抽水状态,此阶段通过控制抽水功率使得此阶段的每时段向***提供的风力发电量相同,多余的风力发电用于抽水蓄能抽水出力。当某一小时内的总的风力发电值小于周期内平均每小时发电值时,则该小时被定义为抽水抽水蓄能机组的发电状态,此阶段每时段的风力发电较少,此时将水从上水库放出发电,通过控制发电功率使得每时段抽水机组发出的电量和风能的发电量协调后向***提供电能相同,这样协调后的整体的风力发电间歇和波动性相比协调前有了明显减少。协调方式如图2。
(4-3)确定目标函数及约束条件
本发明方法考虑了风电场的上网经济和抽水蓄能机组抽水的花费,通过约束风电场输出有功功率使其趋于平滑,同时满足风电场的功率输出限制,最大程度地利用风能,以风电***的上网经济效益最大为目标。
此外,由于风电场出力的随机性,导致风蓄联合出力的取值也具有随机性,因此整个周期内通过抽水机组获取的抽水功率和发电功率均具有随机性。若仍采用传统经济调度约束条件中等式或不等式的形式,则结果将过于保守,且没有实际意义,因此本发明方法以概率的形式来描述约束条件,即约束条件含有置信水平,而置信水平的取值来自于***实际的要求。因此,基于可信性理论的机会约束规划的风蓄联合运行优化调度的数学模型如下:
风蓄联合运行优化目标函数:
Pws(t,k)=Pt XX(t,k)+Pt YY(t,k) (11)
Pp(t.k)=Pw(t,k)Y(t,k)-Pt Y (12)
Pg(t,k)=Pt X-Pw(t,k)X(t,k) (13)
X(t,k)+Y(t,k)=1 (14)
式中:为风电上网经济效益;Pws(t,k)为第t小时第k时段风蓄联合出力值;Pt Y,Pt X分别为风蓄联合***在抽水状态和发电状态下的联合出力值;Pw(t,k)为第t小时第k时段风电机组实际出力;Pw.t为第t小时的风电机组实际出力;Pp(t,k)和Pg(t,k)分别为机组在第t小时第k时段的抽水功率和发电功率;C(t,k)为第t小时第k时刻的上网电价;Ci(t,k)为第t小时第k时刻的抽水电价;T为一个周期即24小时;K为每小时内时段数即4个时段。
约束条件:
①风蓄联合***确保小时内各个时段的协调总出力保持稳定,满足风电场输出有功功率相对平滑:
Pws(t,k)=Pws(t,k+1) (17)
②目标利润的约束条件:
式中:α0为满足目标利润最大的置信水平。
③风蓄联合***出力在抽水状态下约束:
式中:φ1为风蓄联合机组在第t小时第k时刻抽水状态的联合出力约束置信水平;表示机组在第t小时第k时段风电机组实际出力最小值。
④风蓄联合***出力在发电状态下约束:
式中:φ2为风蓄联合机组在第t小时第k时刻的发电状态下的联合出力约束置信水平;表示机组在第t小时第k时段风电机组实际出力最大值;Pmax(t,k)为第t小时第k时段风电场输送到电网的功率最大值。
⑤风蓄联合***出力约束:
Pmin(t,k)≤Pws(t,k)≤Pmax(t,k) (21)
式中:Pmin(t,k)、Pmax(t,k)为第t小时第k时段风电场输送到电网的功率最小值和最大值;
⑥风蓄联合***抽水功率和发电功率约束:
式中:ηg为水泵发电效率;其余见上文解释。
⑦风电机组实际出力约束:
Pdmin(t,k)≤Pw(t,k)≤Pdmax(t,k) (23)
式中:Pdmin(t,k)、Pdmax(t,k)为第t小时第k时段风电机组装机容量的最大值和最小值。
⑧抽水蓄能电站容量约束:
0≤Q(t,k)≤PpgR (27)
0≤Pp(t,k)≤Pgmax(t,k) (29)
式中:PpgR为机组额定容量;c为机组循环效率;Q(t,k)为第t小时第k时段同步机抽水功率(MW);Ppmax(t,k)和Pgmax(t,k)分别为抽水功率和发电功率上限;E(t,k)为抽水蓄能电站存储量;相邻两个时段的间隔为ΔT。
⑨抽水蓄能电站水库约束:
本发明方法将上下水库看成一个整体,并将水库储水量转化为存储电量,不考虑水位对***的影响,本发明方法将存储电量和水量之间看作是线性关系,储水量受到上下限约束,且在***运行一个调度周期后,储水量未发生变化:
Emin(t,k)≤E(t,k)≤Emax(t,k) (30)
E(T,k)=T(t,k) (31)
式中:Emin(t,k)、Emax(t,k)分别表示第t小时第k时刻风蓄联合***储水量的最大和最小值。
⑩在调度周期内,风蓄联合***的储水量变化公式满足:
式中:ηp为水泵抽水效率;ηg为水泵发电效率;ηp取0.85;ηg取0.8。
(5)结合改进的粒子群算法和随机模拟技术对步骤4的模型进行求解
粒子群算法是Kennedy博士和Eberhart教授提出的一种智能寻优方法。其对求解含有大量约束条件的非线性规划问题具有很快的收敛速度及较好的全局寻优能力。粒子群基本公式为:
式中:m表示粒子的编号;t分别表示算法运行的迭代步数;d表示搜索空间维度D的第d维分量;h表示全局最优;Pm为粒子的个体极值,xmd,vmd,pmd,vmax.d,vmin.d,hd分别表示Xm,Vm,Pm,vmax,vmin,h的第d维分量;c1,c2通常取值为2;r1,r2是介于[0,1]之间的随机数。
由于粒子群算法中惯性权重对寻优搜索影响较大,本发明方法参考基于自适应度调整粒子群算法思想对其进行了改进。搜索初始阶段,设定惯性权重取较大值,可以提高寻优搜索的效率,在搜索中后期,随着搜寻结果的不断优化,采用改进后的惯性权重取值公式使惯性权重随着搜索步骤的进行不断减小,从而进一步提高寻优搜索的准确性,获取更多的经济效益。改进后的惯性权重ω可表示为:
式中:ωmax为初始权重;ωmin为最终权重;nmax为最大迭代次数;n0为当前迭代次数。ωmax取0.8。
随机模拟方法,也称为蒙特·卡罗(MonteCarlo)方法。它的基本思想是:首先建立一个概率模型或随机过程求各方面的问题,然后对模型或过程进行分析计算,最后给出所求解的近似值。本发明方法将随机模拟技术结合改进后的粒子群算法来求解机会约束规划模型,提高寻优搜索的精度和收敛速度。具体求解思路如下:
(1)输入***数据,置信水平以及算法控制参数;
(2)以各个小时内机组在抽水状态下和发电状态下的联合出力作为随机变量,在可行域范围内由概率密度函数生成N个随机变量ξN即产生一组初始方案作为粒子群算法的初始种群,并初始化粒子群位置和速度;
(3)任意给定决策变量x,设其初始值m=0,将x和ξ值代入式(3),若满足使gj(x,ξ)≤0成立,则令m=m+1,当N取值较大,且满足m/N≥α时,式(3)成立,若不满足,返回步骤(2);
(4)若满足,此时设初始值为n,判断N次随机序列中若n/N≥β,则开始排序置fi=f(x,ξi),其中i=1,2,…,N,得到序列{f1,f2,…,fN},置K为βiN的整数部分,则序列中第K个最大元素即为目标函数的最优值;
(5)结合步骤(4)及适应度优劣程度确定粒子群算法的个体极值和全局极值;
(6)更新粒子速度和位置,利用式(34)计算惯性权重,并采用式(32)、和式(33)对种群进行迭代计算,得到新一代的种群;
(7)重复步骤(3)至(6),当满足最大允许迭代次数时,停止搜索,并输出结果即为风蓄联合机组在抽水状态下和发电状态下的最优联合出力值。否则返回步骤(3)继续搜索。
Claims (6)
1.一种风蓄联合***优化调度方法,其特征在于:所述方法内容包括如下步骤:
步骤1 定义抽水蓄能机组运行状态,建立抽水蓄能状态模型;
步骤2 考虑风电机组出力特性,对风力发电机组进行建模;
步骤3 基于可信性理论求解机会约束规划;
步骤4 建立基于可信性理论的机会约束规划的风蓄联合优化调度模型;
步骤5 结合改进的粒子群算法和随机模拟技术对步骤4的模型进行求解。
2.根据权利要求1所述的一种风蓄联合***优化调度方法,其特征在于:在步骤1中,所述定义抽水蓄能机组运行状态,建立抽水蓄能状态模型,其具体过程如下:
定义抽水蓄能机组运行状态,建立抽水蓄能机组状态模型;
抽水蓄能机组与风电机组联合后组成新的优化调度,可以充分发挥抽水蓄能机组的优势,对风电机组出力进行有效平滑;根据日前风电出力预测值来确定抽水蓄能机组的运行状态,但是由于风电预测出力存在误差,因此,根据实时的风电出力变化,建立风蓄联合协调出力的新型节能调度方式;
引入0-1整型变量来考虑机组的抽水和发电两个阶段,建立抽水蓄能机组的状态模型如下:
X(t,k)+Y(t,k)=1 (1)
(1)式中:X(t,k)=1时表示机组在第t小时第k(协调周期为24小时,每小时分四个时段,k=1,2,3,4)时段为发电阶段,Y(t,k)=1时机组为抽水阶段。
3.根据权利要求1所述的一种风蓄联合***优化调度方法,其特征在于:在步骤2中,所述考虑风电机组出力特性,其具体过程如下:
①单台风电机出力特性
风速为服从Weibull分布的随机变量,风电机组的功率特性曲线一般由风机制造商提供,也可通过实测得到;在计算中,根据风电机组的功率特性曲线,单台风电机的输出功率PIF与风速关系v可近似用分段函数表示为:
(2)式中:vr和PWr分别表示风机在第t小时第k时段的额定风速和额定功率;vin为切入风速;vout为切出风速;
②描述风电机组总出力曲线特性
总的风电机组的发电功率PW可表示为:
Pw(t,k)=efNW(t,k)PIF(t,k) (3)
0≤PIF(t,k)≤PWR (4)
(3)和(4)式中:PIF(t,k),Pw(t,k),分别为第t小时第k时段单台风机和总的风电机组发电功率(MW);e,f分别为风机的传动效率和发电效率(%);NW(t,k)为第t小时第k时段风电场中正常运行的风机台数;PWR(t,k)为每台机组的额定功率(MW)。
4.根据权利要求1所述的一种风蓄联合***优化调度方法,其特征在于:在步骤3中,所述基于可信性理论求解机会约束规划,其具体求解过程如下:
采用随机模拟技术来求解机会约束规划:
机会约束规划是一种用于解决在给定置信水平下含有随机变量的优化问题的随机规划方法,且考虑到决策变量在随机变量满足约束条件后才能作出决策。本发明方法提出如下模型来求解机会约束规划,具体分析如式(1)~式(3)所示。
Cr{gj(x,ξ)≤0,j=1,2,…,p}≥β (7)
在(5)、(6)和(7)式中:f(x,ξ)为目标函数;为目标函数的乐观值;gj(x,ξ)为约束函数;x为决策向量;ξ为一个随机参数向量;α、β分别为决策者预先给定的置信水平;Cr{·}为{·}中事件的可信性;
假设(Θ,P(Θ),Pos)为可能性空间,A为幂集P(Θ)中的一个元素,则其可信性测度可表示为:
(8)式中:Pos{A}为事件A的可能性测度;Ac为A的对立事件。
5.根据权利要求1所述的一种风蓄联合***优化调度方法,其特征在于:在步骤4中,所述建立基于可信性理论的机会约束规划的风蓄联合优化调度模型,其具体实现过程如下:
确定目标函数及约束条件
以概率的形式来描述约束条件,即约束条件含有置信水平;而置信水平的取值来自于***实际的要求;因此,基于可信性理论的机会约束规划的风蓄联合运行优化调度的数学模型如下:
风蓄联合运行优化目标函数:
Pws(t,k)=Pt XX(t,k)+Pt YY(t,k) (11)
Pp(t.k)=Pw(t,k)Y(t,k)-Pt Y (12)
Pg(t,k)=Pt X-Pw(t,k)X(t,k) (13)
X(t,k)+Y(t,k)=1 (14)
式中:为风电上网经济效益;Pws(t,k)为第t小时第k时段风蓄联合出力值;Pt Y,Pt X分别为风蓄联合***在抽水状态和发电状态下的联合出力值;Pw(t,k)为第t小时第k时段风电机组实际出力;Pw.t为第t小时的风电机组实际出力;Pp(t,k)和Pg(t,k)分别为机组在第t小时第k时段的抽水功率和发电功率;C(t,k)为第t小时第k时刻的上网电价;Ci(t,k)为第t小时第k时刻的抽水电价;T为一个周期即24小时;K为每小时内时段数即4个时段;
约束条件:
①风蓄联合***确保小时内各个时段的协调总出力保持稳定,满足风电场输出有功功率相对平滑:
Pws(t,k)=Pws(t,k+1) (17)
②目标利润的约束条件:
式中:α0为满足目标利润最大的置信水平;
③风蓄联合***出力在抽水状态下约束:
式中:φ1为风蓄联合机组在第t小时第k时刻抽水状态的联合出力约束置信水平;表示机组在第t小时第k时段风电机组实际出力最小值;
④风蓄联合***出力在发电状态下约束:
式中:φ2为风蓄联合机组在第t小时第k时刻的发电状态下的联合出力约束置信水平;表示机组在第t小时第k时段风电机组实际出力最大值;Pmax(t,k)为第t小时第k时段风电场输送到电网的功率最大值;
⑤风蓄联合***出力约束:
Pmin(t,k)≤Pws(t,k)≤Pmax(t,k) (21)
式中:Pmin(t,k)、Pmax(t,k)为第t小时第k时段风电场输送到电网的功率最小值和最大值;
⑥风蓄联合***抽水功率和发电功率约束:
式中:ηg为水泵发电效率;其余见上文解释;
⑦风电机组实际出力约束:
Pdmin(t,k)≤Pw(t,k)≤Pdmax(t,k) (23)
式中:Pdmin(t,k)、Pdmax(t,k)为第t小时第k时段风电机组装机容量的最大值和最小值;
⑧抽水蓄能电站容量约束:
[Pp(t,k)-Q(t,k)]×Y(t,k)=0 (26)
0≤Q(t,k)≤PpgR (27)
0≤Pp(t,k)≤Pgmax(t,k) (29)
式中:PpgR为机组额定容量;c为机组循环效率;Q(t,k)为第t小时第k时段同步机抽水功率(MW);Ppmax(t,k)和Pgmax(t,k)分别为抽水功率和发电功率上限;E(t,k)为抽水蓄能电站存储量;相邻两个时段的间隔为ΔT;
⑨抽水蓄能电站水库约束:
将上下水库看成一个整体,并将水库储水量转化为存储电量,不考虑水位对***的影响,将存储电量和水量之间看作是线性关系,储水量受到上下限约束,且在***运行一个调度周期后,储水量未发生变化:
Emin(t,k)≤E(t,k)≤Emax(t,k) (30)
E(T,k)=T(t,k) (31)
式中:Emin(t,k)、Emax(t,k)分别表示第t小时第k时刻风蓄联合***储水量的最大和最小值;
⑩在调度周期内,风蓄联合***的储水量变化公式满足:
式中:ηp为水泵抽水效率;ηg为水泵发电效率;ηp取0.85;ηg取0.8。
6.根据权利要求1所述的一种风蓄联合***优化调度方法,其特征在于:在步骤5中,所述结合改进的粒子群算法和随机模拟技术对模型进行求解,其具体实现过程如下:
粒子群基本公式为:
式中:m表示粒子的编号;t分别表示算法运行的迭代步数;d表示搜索空间维度D的第d维分量;h表示全局最优;Pm为粒子的个体极值,xmd,vmd,pmd,vmax.d,vmin.d,hd分别表示Xm,Vm,Pm,vmax,vmin,h的第d维分量;c1,c2通常取值为2;r1,r2是介于[0,1]之间的随机数;
采用改进后的惯性权重取值公式使惯性权重随着搜索步骤的进行不断减小,从而进一步提高寻优搜索的准确性,获取更多的经济效益;改进后的惯性权重ω可表示为:
式中:ωmax为初始权重;ωmin为最终权重;nmax为最大迭代次数;n0为当前迭代次数。ωmax取0.8;
将随机模拟技术结合改进后的粒子群算法来求解机会约束规划模型,提高寻优搜索的收敛速度和精度,具体求解思路如下:
(1)输入***数据,置信水平以及算法控制参数;
(2)以各个小时内机组在抽水状态下和发电状态下的联合出力作为随机变量,在可行域范围内由概率密度函数生成N个随机变量ξN即产生一组初始方案作为粒子群算法的初始种群,并初始化粒子群位置和速度;
(3)任意给定决策变量x,设其初始值m=0,将x和ξ值代入式(3),若满足使gj(x,ξ)≤0成立,则令m=m+1,当N取值较大,且满足m/N≥α时,式(3)成立,若不满足,返回步骤(2);
(4)若满足,此时设初始值为n,判断N次随机序列中若n/N≥β,则开始排序置fi=f(x,ξi),其中i=1,2,…,N,得到序列{f1,f2,…,fN},置K为βiN的整数部分,则序列中第K个最大元素即为目标函数的最优值;
(5)结合步骤(4)及适应度优劣程度确定粒子群算法的个体极值和全局极值;
(6)更新粒子速度和位置,利用式(34)计算惯性权重,并采用式(32)、和式(33)对种群进行迭代计算,得到新一代的种群;
(7)重复步骤(3)至(6),当满足最大允许迭代次数时,停止搜索,并输出结果即为风蓄联合机组在抽水状态下和发电状态下的最优联合出力值;否则返回步骤(3)继续搜索。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810313720.1A CN108667077A (zh) | 2018-04-10 | 2018-04-10 | 一种风蓄联合***优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810313720.1A CN108667077A (zh) | 2018-04-10 | 2018-04-10 | 一种风蓄联合***优化调度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108667077A true CN108667077A (zh) | 2018-10-16 |
Family
ID=63783187
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810313720.1A Pending CN108667077A (zh) | 2018-04-10 | 2018-04-10 | 一种风蓄联合***优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108667077A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109449980A (zh) * | 2018-11-28 | 2019-03-08 | 武汉大学 | 一种含海水抽蓄电站的主动配电网优化调度方法 |
CN110135631A (zh) * | 2019-04-26 | 2019-08-16 | 燕山大学 | 基于信息间隙决策理论的电气综合能源***调度方法 |
CN111882131A (zh) * | 2020-07-31 | 2020-11-03 | 南方电网调峰调频发电有限公司 | 一种基于双工况最优效率跟踪的变速海水抽蓄电站容量优化方法 |
CN113809773A (zh) * | 2021-09-22 | 2021-12-17 | 南方电网科学研究院有限责任公司 | 基于海上风电机组的电力***生产模拟方法、装置及设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2495839A1 (de) * | 2011-03-02 | 2012-09-05 | ABB Technology AG | Energiesystem |
US20130024044A1 (en) * | 2011-07-18 | 2013-01-24 | Nec Laboratories America, Inc. | Optimal storage sizing for integrating wind and load forecast uncertainties |
CN104124708A (zh) * | 2014-07-25 | 2014-10-29 | 国家电网公司 | 一种能源协调优化机组组合方法 |
CN104362681A (zh) * | 2014-11-18 | 2015-02-18 | 湖北省电力勘测设计院 | 一种考虑随机性的孤岛微电网容量优化配置方法 |
CN105281320A (zh) * | 2015-09-16 | 2016-01-27 | 河海大学 | 协调经济性及切负荷率公平性的紧急负荷控制优化方法 |
CN105576709A (zh) * | 2016-01-06 | 2016-05-11 | 南京工程学院 | 一种基于混合算法的风蓄火联合运行的优化方法 |
JP2017084360A (ja) * | 2015-10-23 | 2017-05-18 | 富士通株式会社 | 太陽光発電システムの作動方法 |
CN107171367A (zh) * | 2017-05-19 | 2017-09-15 | 山东大学 | 风电与储能***互补下的火电机组组合优化决策方法 |
CN107171365A (zh) * | 2017-06-23 | 2017-09-15 | 国网山东省电力公司经济技术研究院 | 基于场景解耦和异步迭代的多目标随机动态经济调度方法 |
CN107679679A (zh) * | 2017-11-20 | 2018-02-09 | 国网辽宁省电力有限公司电力科学研究院 | 联合发电机组调度运行方法 |
-
2018
- 2018-04-10 CN CN201810313720.1A patent/CN108667077A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2495839A1 (de) * | 2011-03-02 | 2012-09-05 | ABB Technology AG | Energiesystem |
US20130024044A1 (en) * | 2011-07-18 | 2013-01-24 | Nec Laboratories America, Inc. | Optimal storage sizing for integrating wind and load forecast uncertainties |
CN104124708A (zh) * | 2014-07-25 | 2014-10-29 | 国家电网公司 | 一种能源协调优化机组组合方法 |
CN104362681A (zh) * | 2014-11-18 | 2015-02-18 | 湖北省电力勘测设计院 | 一种考虑随机性的孤岛微电网容量优化配置方法 |
CN105281320A (zh) * | 2015-09-16 | 2016-01-27 | 河海大学 | 协调经济性及切负荷率公平性的紧急负荷控制优化方法 |
JP2017084360A (ja) * | 2015-10-23 | 2017-05-18 | 富士通株式会社 | 太陽光発電システムの作動方法 |
CN105576709A (zh) * | 2016-01-06 | 2016-05-11 | 南京工程学院 | 一种基于混合算法的风蓄火联合运行的优化方法 |
CN107171367A (zh) * | 2017-05-19 | 2017-09-15 | 山东大学 | 风电与储能***互补下的火电机组组合优化决策方法 |
CN107171365A (zh) * | 2017-06-23 | 2017-09-15 | 国网山东省电力公司经济技术研究院 | 基于场景解耦和异步迭代的多目标随机动态经济调度方法 |
CN107679679A (zh) * | 2017-11-20 | 2018-02-09 | 国网辽宁省电力有限公司电力科学研究院 | 联合发电机组调度运行方法 |
Non-Patent Citations (4)
Title |
---|
HUAJIE DING,等: ""Stochastic optimization of the daily operation of wind farm and pumped-hydro-storage plant"", 《RENEWABLE ENERGY》 * |
于佳,等: ""基于机会约束规划的风–蓄联合动态经济调度"", 《电网技术》 * |
张刘冬,等: ""基于成本效益分析的风电–抽水蓄能联合运行优化调度模型"", 《电网技术》 * |
沈洲,等: ""基于机会约束规划和随机模拟技术的含风电场电力***多目标优化调度 "", 《电力学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109449980A (zh) * | 2018-11-28 | 2019-03-08 | 武汉大学 | 一种含海水抽蓄电站的主动配电网优化调度方法 |
CN110135631A (zh) * | 2019-04-26 | 2019-08-16 | 燕山大学 | 基于信息间隙决策理论的电气综合能源***调度方法 |
CN110135631B (zh) * | 2019-04-26 | 2022-02-22 | 燕山大学 | 基于信息间隙决策理论的电气综合能源***调度方法 |
CN111882131A (zh) * | 2020-07-31 | 2020-11-03 | 南方电网调峰调频发电有限公司 | 一种基于双工况最优效率跟踪的变速海水抽蓄电站容量优化方法 |
CN113809773A (zh) * | 2021-09-22 | 2021-12-17 | 南方电网科学研究院有限责任公司 | 基于海上风电机组的电力***生产模拟方法、装置及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107276127B (zh) | 考虑联络线电力交易计划的多区域互联电力***的风电消纳优化方法 | |
JP7261507B2 (ja) | 電気ヒートポンプ-熱電併給システムを最適化する調整方法及びシステム | |
CN106099993B (zh) | 一种适应新能源大规模接入的电源规划方法 | |
CN113822496B (zh) | 一种多机组热电厂供热模式及参数在线寻优方法 | |
CN109508499A (zh) | 多时段多场景分布式电源最优接入位置与容量研究方法 | |
CN108009693A (zh) | 基于两级需求响应的并网微电网双层优化方法 | |
CN108667077A (zh) | 一种风蓄联合***优化调度方法 | |
CN108039737B (zh) | 一种源网荷协调运行模拟*** | |
CN111064229A (zh) | 基于q学习的风-光-气-蓄联合动态经济调度优化方法 | |
CN103699941A (zh) | 一种电力***调度运行年方案的制定方法 | |
CN103762589A (zh) | 一种电网中新能源容量配比分层优化方法 | |
CN109245180A (zh) | 一种风光火储协调优化运行方法 | |
CN109103929A (zh) | 基于改进动态克里丝金模型的配电网经济优化调度方法 | |
CN110957717A (zh) | 一种多电源电力***多目标日前优化调度方法 | |
CN114707403B (zh) | 基于抽水蓄能调节的地区配电网多能源协调优化调度方法 | |
CN109740808A (zh) | 一种风光水互补发电计划计算方法及*** | |
CN114519543B (zh) | 一种乡村多能源***边缘自治运行方法及*** | |
CN104809545B (zh) | 一种虚拟电厂运行建模方法 | |
Li et al. | Optimized operation of hybrid system integrated with MHP, PV and PHS considering generation/load similarity | |
Liu et al. | Two-stage optimal economic scheduling for commercial building multi-energy system through internet of things | |
Shi et al. | Research on energy management of hydrogen electric coupling system based on deep reinforcement learning | |
CN114154744A (zh) | 综合能源***的扩容规划方法、装置及电子设备 | |
CN114066204A (zh) | 一种综合能源***的一体化优化规划与运行方法及装置 | |
CN105633950B (zh) | 一种考虑风电注入不确定性的多目标随机模糊动态最优潮流求解方法 | |
Ma et al. | Integrated strategy of the output planning and economic operation of the combined system of wind turbines-pumped-storage-thermal power units |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181016 |