CN107579515B - 一种改善电厂调频性能的容量配置方法 - Google Patents

一种改善电厂调频性能的容量配置方法 Download PDF

Info

Publication number
CN107579515B
CN107579515B CN201710636185.9A CN201710636185A CN107579515B CN 107579515 B CN107579515 B CN 107579515B CN 201710636185 A CN201710636185 A CN 201710636185A CN 107579515 B CN107579515 B CN 107579515B
Authority
CN
China
Prior art keywords
energy storage
frequency modulation
value
unit
time
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.)
Active
Application number
CN201710636185.9A
Other languages
English (en)
Other versions
CN107579515A (zh
Inventor
陈丽娟
姜宇轩
许晓慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201710636185.9A priority Critical patent/CN107579515B/zh
Publication of CN107579515A publication Critical patent/CN107579515A/zh
Application granted granted Critical
Publication of CN107579515B publication Critical patent/CN107579515B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明公开了一种改善电厂调频性能的容量配置方法,包括根据《两个细则》中AGC的考核指标,提出调节速率、调节精度,响应时间三项指标的具体计算方法;提出不同阶段提升调频各项指标的储能充放电策略;计算储能参与电厂调频的收益;计算储能参与电厂调频的成本,同时计及储能频繁动作对其寿命的影响;建立储能调频的经济性模型;以净收益最大为目标函数,采用差分进化算法计算保证净收益最大的储能容量。本发明的方法可以较为准确地计算出Kp值,为后续充放电策略制定以及有效性的计算奠定基础;可以显著提升电厂的AGC性能;并将修正后的寿命纳入储能成本项,采用差分进化算法计算出系统配置储能的最优容量和最大净收益。

Description

一种改善电厂调频性能的容量配置方法
技术领域
本发明涉及储能调频方法,特别是涉及一种改善电厂调频性能的容量配置方法。
背景技术
大规模新能源的并网缓解了电力紧张,但其本身的间歇性和不确定性给电网的电能质量带来了严重的影响。其中,频率作为电能质量的一项重要指标,对电力系统的稳定运行有着关键的影响。
在我国,传统机组响应速度和精度有限,不能完全满足当下调频需要。而储能作为一种新兴的产业,在世界电力调频市场得到越来越广泛的认可。
考虑到储能参与电厂调频,经济性最优是首要目标。在现有的储能参与电网调频中,很少与现有政策结合,并且实用性较低。
发明内容
发明目的:针对现有技术的缺陷,提供一种改善电厂调频性能的容量配置方法。
技术方案:一种改善电厂调频性能的容量配置方法,该方法包括如下步骤:
(1)根据《两个细则》中AGC的考核指标,提出调节速率、调节精度,响应时间三项指标的具体计算方法;
(2)提出不同阶段提升调频各项指标的储能充放电策略;
(3)计算储能参与电厂调频的收益;
(4)计算储能参与电厂调频的成本,同时计及储能频繁动作对其寿命的影响;
(5)建立储能调频的经济性模型;
(6)以净收益最大为目标函数,采用差分进化算法计算保证净收益最大的储能容量。
进一步的,所述步骤(1)包括:
(11)计算调节速率指标K1
判断每段有效区间的起始时刻和结束时刻,当机组出力跨过死区即判断为有效计算区间的开始;而当AGC指令和机组的出力绝对值之差小于预定限值,即可算作有效区间的结束;即满足下式:
Figure BDA0001364918230000021
其中,t1为计算响应速率区间的开始时刻,t2为计算响应速率区间的结束时刻,Pd1和Pd2分别为判断计算响应速率区间的阈值;Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
调节速率指标K1的计算公式为:
Figure BDA0001364918230000022
Figure BDA0001364918230000023
Figure BDA0001364918230000024
其中,PSi和PEi分别是第i次调节过程开始和结束时刻的机组出力,TSi和TEi分别是第i次调节过程的开始和结束时刻,K1表示调节速率指标,vi表示第i次调节的速率,v是n次调节的平均速率,vN是机组标准调节速率,与机组类型有关;
(12)计算调节精度指标K2
其有效区间选取原则为AGC指令保持不变,实际出力和设点出力之间的差值始终小于预定阈值,如下式:
Figure BDA0001364918230000025
其中,const为常数,Pa(t)表示t时刻AGC的指令值,Pd3为判断计算响应精度区间的阈值,Pg(t)表示t时刻机组的出力值;
调节精度指标K2的计算公式如下:
Figure BDA0001364918230000031
Figure BDA0001364918230000032
其中,△Pi表示平稳运行阶段的功率偏差均值,Pi(t)为机组平稳运行阶段的出力,Pj为平稳运行阶段AGC的指令值,TSi和TEi分别为平稳运行阶段的开始和结束时刻,K2表示调节精度指标,PN表示平稳运行阶段机组的额定功率;
(13)计算响应时间指标K3
响应时间是指EMS系统发出指令之后,机组出力在原出力点的基础上,可靠地跨出与调节方向一致的调节死区所用的时间;K3计算公式如下式:
Figure BDA0001364918230000033
Figure BDA0001364918230000034
其中,△Ti是第i次的响应时间,△T是响应时间的平均值,TN是响应时间基准值;
计算响应时间指标需要判断有效区间的开始和结束时刻,如下式所示:
Figure BDA0001364918230000035
其中,Pd4和Pd5分别为判断计算响应时间指标区间的阈值值,t3为判断计算响应时间区间的开始时刻,t4为判断计算响应时间区间的结束时刻。
进一步的,所述步骤(2)包括:
(21)在机组响应AGC指令阶段,提升机组的响应速度
在AGC指令下调阶段,当Pa(t)-Pg(t)<0,储能为充电状态,充电功率如下式:
Pc(t)=max(-Pmax,Pa(t)-Pg(t))
其中,Pc(t)表示t时刻储能的充电功率,Pmax表示储能的额定功率,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
在AGC指令上调阶段,当Pa(t)-Pg(t)>0,储能为放电状态,放电功率如下式:
Pd(t)=min(Pmax,Pa(t)-Pg(t))
其中,Pd(t)表示t时刻储能的放电功率,Pmax表示储能的额定功率,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
(22)在机组进入稳定出力阶段,提升机组的响应精度
判断机组进入稳定出力区间后,利用储能进行充放电补偿,平滑综合出力,此时储能出力如下式:
P(t)=Pa(t)-Pg(t)
其中,P(t)表示储能的出力,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
(23)在系统下发AGC指令时,降低机组的响应时间
当系统下发AGC指令,机组并不能及时动作,当(Pa(t)-Pa(tn))(Pg(t)-Pg(tn))<0时,说明机组出现反方向调节,此时储能应及时动作,降低响应时间,功率如下式:
P(t)=Pg(tn)-Pg(t)
其中,P(t)表示储能的出力,tn是响应时间的开始时刻,t是当前时刻,Pg(tn)表示tn时刻机组的出力值,Pg(t)表示t时刻机组的出力值。
进一步的,所述步骤(3)包括:
根据AGC性能指标Kp值确定机组参与AGC调频服务的经济补偿
Figure BDA0001364918230000041
如下式:
Figure BDA0001364918230000042
Di=|Pg(ti)-Pg(ti-1)|
其中,YAGC表示AGC调节性能补偿标准,Di表示第i次AGC调节深度,ti为第i次调节的时刻,Pg(ti)表示第i次调节的时刻机组的出力值,Pg(ti-1)表示第i-1次调节的时刻机组的出力值。
进一步的,所述步骤(4)包括:
(41)计算储能的成本
储能的成本由初始建设投资成本和运行维护成本两部分组成;
储能的初始建设投资成本Cin的计算公式如下式:
Cin=Cb+Cpcs+Caf
其中,Cb是蓄电池成本,Cpcs是能量转换装置成本,Caf是辅助设施成本;
初始投资成本的等年值Cin,year的计算公式如下式:
Figure BDA0001364918230000051
其中,r是贴现率,a是储能的使用寿命;
储能的运行维护成本COM的计算公式如下式:
COM=Cfix+Cvar
其中,Cfix是储能固定运行成本,Cvar是可变运行维护成本;
(42)计算储能参与调频后的实际寿命
储能的频繁动作会降低其使用寿命,需计算储能参与调频后的实际寿命N,其中,储能一天的等效循环次数的计算公式如下:
m=∑(soc(i)-soc(i+1)) soc(i)>soc(i+1)
其中,m是储能一天的等效循环次数,soc(i)表示时刻i的储能荷电状态;
储能参与调频后的实际寿命计算公式如下:
Figure BDA0001364918230000052
其中,a是蓄电池的出厂寿命,q是蓄电池的最大循环次数,N是储能参与调频后的实际寿命;
然后对储能投资成本项中的蓄电池的出厂寿命a进行修正得到实际的成本项C'in,year
Figure BDA0001364918230000061
更进一步的,所述步骤(41)包括:
所述的储能各项建设投资成本计算公式如下:
Figure BDA0001364918230000062
Cpcs=CpPb
Caf=CAEb
其中,CE是储能本体的单位能量价格,Cp是PCS的单位功率价格,CA是辅助设施的单位能量价格,Eb是储能的额定容量,Pb是储能的额定功率,η是储能的充放电效率;
所述的储能的固定运行成本Cfix和可变运行维护成本Cvar计算公式如下式:
Cfix=CFPb
Figure BDA0001364918230000063
其中,CF是单位功率固定运行成本,Pb表示储能的额定功率,Cc是当前的充电电价,P是储能的充电功率,η表示储能的充放电效率。
进一步的,所述步骤(5)包括:
储能调频的经济性模型为储能参与调频的净收益:
Figure BDA0001364918230000064
其中,C表示储能参与调频的净收益,
Figure BDA0001364918230000071
是机组参与AGC调频服务的经济补偿,Cin是储能的初始建设投资成本,COM是储能的运行维护成本。
进一步的,所述步骤(6)中差分进化算法的步骤如下:
(1)生成初始种群
在n1维空间里随机生成满足约束条件的h个染色体,其中,h为种群规模,n1为变量个数,
Figure BDA0001364918230000072
分别为求解变量的上下限值,xi(g)表示第i条染色体或个体的g代值,第xi,j(g)表示第i个染色体上的第j个基因的g代值,rand(0,1)表示0和1之间的随机数;在实际算例中
Figure BDA0001364918230000073
分别为取为储能容量的上下限值;初始种群xi,j(0)计算公式为:
Figure BDA0001364918230000074
(2)计算不同储能容量下储能参与调频净收益的相反数,且差分进化算法默认为求最小值,因此求最大值问题需要将目标函数前面加上负号;
(3)判断是否满足终止条件
若满足终止条件,则直接输出最优结果;若不满足终止条件,则执行步骤(4);
(4)变异操作
从群体中随机选择3个染色体
Figure BDA0001364918230000076
且(p1≠p2≠p3),其中
Figure BDA0001364918230000077
为差异化向量,η1为缩放因子,vi(g+1)为变异后的第g+1代个体计算公式为:
Figure BDA0001364918230000078
(5)交叉操作
具体操作如下:
Figure BDA0001364918230000075
其中,ui,j(g+1)为交叉操作后第i条染色体第j个基因g+1代的值,CR为交叉概率,jrand为[1,2,...,n1]中的随机数;
(6)选择操作
Figure BDA0001364918230000081
其中,i=1,2,...,h,xi(g+1)为染色体i的g+1代值,f(·)为问题的目标函数(默认最小值问题,最大值加负号),实际算例中即为储能参与电厂调频净收益的相反数;
(7)令g=g+1,返回步骤(2)继续操作,直至达到最大的迭代次数Tmax
进一步的,所述的储能参与调频的约束条件如下:
1)储能充放电功率约束
Figure BDA0001364918230000082
其中,Pd(t)是储能t时刻的放电功率,Pd,max是储能放电功率的最大值,Pc(t)是储能t时刻的充电功率,Pc,max是储能充电功率的最大值;
2)储能SOC约束
S1≤SOC(t)≤S2
其中,S1是时刻t的储能荷电状态的下限,S2是时刻t的储能荷电状态的上限。
有益效果:与现有技术相比,本发明的方法提出了AGC性能计算方法,可以较为准确地计算出Kp值,为后续充放电策略制定以及有效性的计算奠定基础。其次,本发明从提升调节速率,调节精度,降低响应时间三方面出发,制定储能的充放电策略,可以显著提升电厂的AGC性能。再次,本发明定量分析了储能频繁动作对寿命的影响,并将修正后的寿命纳入储能成本项,采用差分进化算法计算出系统配置储能的最优容量和最大净收益。
附图说明
图1是储能参与电厂调频的流程图;
图2是典型AGC机组设点控制过程示意图;
图3是差分进化算法流程图。
具体实施方式
下面结合附图对本发明的技术方案进行详细的说明。
如图1所示,本发明提供的一种改善电厂调频性能的容量配置方法,该方法给出了储能参与电厂调频的充放电策略,在此基础上,建立了经济性模型,计算出电厂配置的储能最优容量。具体包含以下步骤:
步骤1,根据《两个细则》中AGC(Automatic Generation Control自动发电控制)的考核指标,提出调节速率、调节精度和响应时间三项指标的具体计算方法。
步骤1.1,计算机组的调节速率指标K1
如图2中T1-T2和T4-T5时段,机组的调节速率指标K1的计算公式如下:
Figure BDA0001364918230000091
Figure BDA0001364918230000092
Figure BDA0001364918230000093
其中,PSi和PEi分别是第i次调节过程开始和结束时刻的机组出力,TSi和TEi分别是第i次调节过程的开始和结束时刻,K1表示调节速率指标,vi表示第i次调节的速率,v是n次调节的平均速率,vN是机组标准调节速率,与机组类型有关,单位是MW/min;
计算调节速率指标K1的关键是判断每段有效区间的起始时刻和结束时刻,当机组出力跨过死区即判断为有效计算区间的开始;而当AGC指令和机组的出力绝对值之差小于预定限值,即可算作有效区间的结束。即满足下式:
Figure BDA0001364918230000094
其中,t1为计算响应速率区间的开始时刻,t2为计算响应速率区间的结束时刻,Pd1和Pd2分别为判断计算响应速率区间的阈值;Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
步骤1.2,计算机组的调节精度指标K2
如图2中的T2-T3时段,在机组平稳运行阶段,机组出力围绕P2轻微波动,类似这样的时段,对实际出力与设点指令之差的绝对值进行积分,再除以积分时间,即为该时段的调节偏差量。计算公式如下:
Figure BDA0001364918230000101
Figure BDA0001364918230000102
其中,△Pi表示平稳运行阶段的功率偏差均值,Pi(t)为机组平稳运行阶段的出力,Pj为平稳运行阶段AGC的指令值,TSi和TEi分别为平稳运行阶段的开始和结束时刻,K2表示调节精度指标,PN表示平稳运行阶段机组的额定功率。
调节精度指标K2的有效区间选取原则为AGC指令保持不变,实际出力和设点出力之间的差值始终小于预定阈值,如下式:
Figure BDA0001364918230000103
其中,Pa表示有效区间内AGC的指令值,const为常数,Pd3为判断计算响应精度区间的阈值,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
步骤1.3,计算机组的响应时间指标K3
响应时间是指EMS系统发出指令之后,机组出力在原出力点的基础上,可靠地跨出与调节方向一致的调节死区所用的时间。如图2中,△Ti=T1-T0。K3计算公式如下式:
Figure BDA0001364918230000104
Figure BDA0001364918230000111
其中,△Ti是第i次的响应时间,△T是响应时间的平均值,TN是响应时间基准值,与机组类型相关。
计算响应时间指标需要判断有效区间的开始和结束时刻,如下式所示:
Figure BDA0001364918230000112
其中,Pd4和Pd5分别为判断计算响应时间指标区间的阈值值,t3为判断计算响应时间区间的开始时刻,t4为判断计算响应时间区间的结束时刻。
步骤2,提出不同阶段提升调频各项指标的储能充放电策略
步骤2.1,在机组响应AGC指令阶段,需提升机组的响应速度
通常在AGC指令下调阶段,如图2中T4到T5阶段,当Pa(t)-Pg(t)<0,储能应为充电状态,充电功率如下式:
Pc(t)=max(-Pmax,Pa(t)-Pg(t))
其中,Pc(t)表示t时刻储能的充电功率,Pmax表示储能的额定功率,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值。
通常在AGC指令上调阶段,如图2中T1到T2阶段,当Pa(t)-Pg(t)>0,储能应为放电状态,放电功率如下式:
Pd(t)=min(Pmax,Pa(t)-Pg(t))
其中,Pd(t)表示t时刻储能的放电功率,Pmax表示储能的额定功率,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值。
储能参与系统调频的充放电功率约束条件如下:
Figure BDA0001364918230000113
其中,Pd(t)是储能t时刻的放电功率,Pd,max是储能放电功率的最大值,Pc(t)是储能t时刻的充电功率,Pc,max是储能充电功率的最大值。
步骤2.2,在机组进入稳定出力阶段,需提升机组的响应精度
判断机组进入稳定出力区间后,如图2中T2到T3阶段,利用储能进行充放电补偿,平滑综合出力,此时储能出力如下式:
P(t)=Pa(t)-Pg(t)
其中,P(t)表示储能出力,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值。
步骤2.3,在系统下发AGC指令时,需降低机组的响应时间
当系统下发AGC指令,机组并不能及时动作,如图2中T0到T1阶段,当(Pa(t)-Pa(tn))(Pg(t)-Pg(tn))<0时,其中,tn是响应时间的开始时刻,t是当前时刻。说明机组出现反方向调节,此时储能应及时动作,降低响应时间,功率应如下式:
P(t)=Pg(tn)-Pg(t)
其中,P(t)表示储能的出力,tn是响应时间的开始时刻,t是当前时刻,Pg(tn)表示tn时刻机组的出力值,Pg(t)表示t时刻机组的出力值。
步骤3,计算储能参与电厂调频的收益
根据AGC性能指标Kp值可以确定机组参与AGC调频服务的经济补偿
Figure BDA0001364918230000121
如下式:
Figure BDA0001364918230000122
Di=|Pg(ti)-Pg(ti-1)|
其中,YAGC表示AGC调节性能补偿标准,Kp表示AGC性能指标,Di表示第i次AGC调节深度,ti为第i次调节的时刻,Pg(ti)表示第i次调节的时刻机组的出力值,Pg(ti-1)表示第i-1次调节的时刻机组的出力值。
步骤4,计及储能寿命计算其参与调频的成本
步骤4.1,计算储能的成本
储能的成本通常由初始建设投资成本和运行维护成本两部分组成。储能的建设投资成本Cin的计算公式如下式:
Cin=Cb+Cpcs+Caf
其中,Cb是蓄电池成本,Cpcs是能量转换装置成本,Caf是辅助设施成本。
初始投资成本的等年值Cin,year的计算公式如下式:
Figure BDA0001364918230000131
其中,r是贴现率,a是储能的使用寿命;
储能的运行维护成本COM的计算公式如下式:
COM=Cfix+Cvar
其中,Cfix是蓄电池固定运行成本,Cvar是可变运行维护成本。
本发明中储能各项建设投资成本及投资成本的等年值计算公式如下:
Figure BDA0001364918230000132
Cpcs=CpPb
Caf=CAEb
其中,CE是储能本体的单位能量价格,Cp是PCS的单位功率价格,CA是辅助设施的单位能量价格,Eb是储能的额定容量,Pb是储能的额定功率,η是储能的充放电效率。
储能的固定运行成本Cfix和可变运行维护成本Cvar计算公式如下式:
Cfix=CFPb
Figure BDA0001364918230000133
其中,CF是单位功率固定运行成本,Pb表示储能的额定功率,Cc是当前的充电电价,P是储能的充电功率,η表示储能的充放电效率。
步骤4.2,计算储能参与调频后的实际寿命
储能的频繁的动作将会降低其使用寿命,因此需计算储能参与调频后的实际寿命N。
为评估储能频繁动作对寿命的影响,首先需要统计出电池平均每天的等效充放电循环次数m。由于SOC的回归策略使得一天中储能的充电放电电量近似相等,因此可以将储能一天中所有放电或者充电的荷电数进行累加,得到等效充放电循环次数m,其计算公式如下:
m=∑(soc(i)-soc(i+1)) soc(i)>soc(i+1)
其中,soc(i)表示表示时刻i的储能荷电状态。
再通过与电池的额定循环次数进行比较,得到参与系统调频后的实际寿命N,如下式:
Figure BDA0001364918230000141
其中,a是蓄电池的出厂寿命,q是蓄电池的最大循环次数,N是储能参与调频后的实际寿命。
然后对储能成本项中的蓄电池的出厂寿命a进行修正得到实际的成本项C'in,year
Figure BDA0001364918230000142
储能参与电厂调频的SOC约束条件如下:
S1≤SOC(t)≤S2
其中,S1是时刻t的储能荷电状态的下限,S2是时刻t的储能荷电状态的上限。
步骤5,建立储能调频经济性模型
该模型为储能参与电厂调频的净收益,计算公式为:
Figure BDA0001364918230000151
其中,C表示储能参与电厂调频的净收益。
步骤6,以净收益最大为目标函数,采用差分进化算法计算保证净收益最大的储能容量
储能的最优容量配置属于单变量单目标优化问题,优化目标为净收益最大,优化变量为储能容量。采用差分进化算法对经济性模型进行求解,计算电厂配置储能的最优容量。
差分进化算法的步骤如下:
(1)生成初始种群
在n1维空间里随机生成满足约束条件的h个染色体。在本发明的实施例中,h为种群规模,n1为变量个数,
Figure BDA0001364918230000152
分别为求解变量的上下限值,xi(g)表示第i条染色体(或个体)的g代值,第xi,j(g)表示第i个染色体上的第j个基因的g代值,rand(0,1)表示0和1之间的随机数。在实际算例中
Figure BDA0001364918230000153
分别为取为储能容量的上下限值。初始种群xi,j(0)计算公式为:
Figure BDA0001364918230000154
(2)计算种群中每个个体的适应值,在实际算例中即为计算不同储能容量下储能参与调频净收益的相反数。且差分进化算法默认为求最小值,因此求最大值问题需要将目标函数前面加上负号。
(3)判断是否满足终止条件
若满足终止条件,则直接输出最优结果;若不满足终止条件,则执行步骤(4),实际算例中终止条件即为最大迭代次数。
(4)变异操作
从群体中随机选择3个染色体
Figure BDA0001364918230000155
且(p1≠p2≠p3),其中
Figure BDA0001364918230000156
为差异化向量,η1为缩放因子。vi(g+1)为变异后的第g+1代个体计算公式为:
Figure BDA0001364918230000163
(5)交叉操作
交叉操作是为了增加群体的多样性,同时为了确保变异体vi(g+1)的每个染色体至少有一个基因遗传给下一代,具体操作如下:
Figure BDA0001364918230000161
其中,ui,j(g+1)为交叉操作后第i条染色体第j个基因g+1代的值,CR为交叉概率,jrand为[1,2,...,n1]中的随机数。
(6)选择操作
Figure BDA0001364918230000162
其中,i=1,2,...,h,xi(g+1)为染色体i的g+1代值,f(·)为问题的目标函数(默认最小值问题,最大值加负号),实际算例中即为储能参与电厂调频净收益的相反数。
(7)令g=g+1,返回步骤(2)继续操作,直至达到最大的迭代次数Tmax

Claims (8)

1.一种改善电厂调频性能的容量配置方法,其特征在于,该方法包括如下步骤:
(1)根据《两个细则》中AGC的考核指标,提出调节速率、调节精度,响应时间三项指标的具体计算方法;
(2)提出不同阶段提升调频各项指标的储能充放电策略;具体为:
(21)在机组响应AGC指令阶段,提升机组的响应速度
在AGC指令下调阶段,当Pa(t)-Pg(t)<0,储能为充电状态,充电功率如下式:
Pc(t)=max(-Pmax,Pa(t)-Pg(t))
其中,Pc(t)表示t时刻储能的充电功率,Pmax表示储能的额定功率,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
在AGC指令上调阶段,当Pa(t)-Pg(t)>0,储能为放电状态,放电功率如下式:
Pd(t)=min(Pmax,Pa(t)-Pg(t))
其中,Pd(t)表示t时刻储能的放电功率,Pmax表示储能的额定功率,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
(22)在机组进入稳定出力阶段,提升机组的响应精度
判断机组进入稳定出力区间后,利用储能进行充放电补偿,平滑综合出力,此时储能出力如下式:
P(t)=Pa(t)-Pg(t)
其中,P(t)表示储能的出力,Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
(23)在系统下发AGC指令时,降低机组的响应时间
当系统下发AGC指令,机组并不能及时动作,当(Pa(t)-Pa(tn))(Pg(t)-Pg(tn))<0时,说明机组出现反方向调节,此时储能应及时动作,降低响应时间,功率如下式:
P(t)=Pg(tn)-Pg(t)
其中,P(t)表示储能的出力,tn是响应时间的开始时刻,t是当前时刻,Pg(tn)表示tn时刻机组的出力值,Pg(t)表示t时刻机组的出力值;
(3)计算储能参与电厂调频的收益;
(4)计算储能参与电厂调频的成本,同时计及储能频繁动作对其寿命的影响;
(5)建立储能调频的经济性模型;
(6)以净收益最大为目标函数,采用差分进化算法计算保证净收益最大的储能容量。
2.根据权利要求1所述的改善电厂调频性能的容量配置方法,其特征在于,所述步骤(1)包括:
(11)计算调节速率指标K1
判断每段有效区间的起始时刻和结束时刻,当机组出力跨过死区即判断为有效计算区间的开始;而当AGC指令和机组的出力绝对值之差小于预定限值,即可算作有效区间的结束;即满足下式:
Figure FDA0002361907770000021
其中,t1为计算响应速率区间的开始时刻,t2为计算响应速率区间的结束时刻,Pd1和Pd2分别为判断计算响应速率区间的阈值;Pa(t)表示t时刻AGC的指令值,Pg(t)表示t时刻机组的出力值;
调节速率指标K1的计算公式为:
Figure FDA0002361907770000022
Figure FDA0002361907770000023
Figure FDA0002361907770000024
其中,PSi和PEi分别是第i次调节过程开始和结束时刻的机组出力,TSi和TEi分别是第i次调节过程的开始和结束时刻,K1表示调节速率指标,vi表示第i次调节的速率,v是n次调节的平均速率,vN是机组标准调节速率,与机组类型有关;
(12)计算调节精度指标K2
其有效区间选取原则为AGC指令保持不变,实际出力和设点出力之间的差值始终小于预定阈值,如下式:
Figure FDA0002361907770000031
其中,const为常数,Pa(t)表示t时刻AGC的指令值,Pd3为判断计算响应精度区间的阈值,Pg(t)表示t时刻机组的出力值;
调节精度指标K2的计算公式如下:
Figure FDA0002361907770000032
Figure FDA0002361907770000033
其中,△Pi表示平稳运行阶段的功率偏差均值,Pi(t)为机组平稳运行阶段的出力,Pj为平稳运行阶段AGC的指令值,TSi和TEi分别为平稳运行阶段的开始和结束时刻,K2表示调节精度指标,PN表示平稳运行阶段机组的额定功率;
(13)计算响应时间指标K3
响应时间是指EMS系统发出指令之后,机组出力在原出力点的基础上,可靠地跨出与调节方向一致的调节死区所用的时间;K3计算公式如下式:
Figure FDA0002361907770000034
Figure FDA0002361907770000035
其中,△Ti是第i次的响应时间,△T是响应时间的平均值,TN是响应时间基准值;
计算响应时间指标需要判断有效区间的开始和结束时刻,如下式所示:
Figure FDA0002361907770000041
其中,Pd4和Pd5分别为判断计算响应时间指标区间的阈值值,t3为判断计算响应时间区间的开始时刻,t4为判断计算响应时间区间的结束时刻。
3.根据权利要求1所述的改善电厂调频性能的容量配置方法,其特征在于,所述步骤(3)包括:
根据AGC性能指标Kp值确定机组参与AGC调频服务的经济补偿
Figure FDA0002361907770000042
如下式:
Figure FDA0002361907770000043
Di=|Pg(ti)-Pg(ti-1)|
其中,YAGC表示AGC调节性能补偿标准,Di表示第i次AGC调节深度,ti为第i次调节的时刻,Pg(ti)表示第i次调节的时刻机组的出力值,Pg(ti-1)表示第i-1次调节的时刻机组的出力值。
4.根据权利要求1所述的改善电厂调频性能的容量配置方法,其特征在于,所述步骤(4)包括:
(41)计算储能的成本
储能的成本由初始建设投资成本和运行维护成本两部分组成;
储能的初始建设投资成本Cin的计算公式如下式:
Cin=Cb+Cpcs+Caf
其中,Cb是蓄电池成本,Cpcs是能量转换装置成本,Caf是辅助设施成本;
初始投资成本的等年值Cin,year的计算公式如下式:
Figure FDA0002361907770000044
其中,r是贴现率,a是储能的使用寿命;
储能的运行维护成本COM的计算公式如下式:
COM=Cfix+Cvar
其中,Cfix是储能固定运行成本,Cvar是可变运行维护成本;
(42)计算储能参与调频后的实际寿命
储能的频繁动作会降低其使用寿命,需计算储能参与调频后的实际寿命N,其中,储能一天的等效循环次数的计算公式如下:
m=∑(soc(i)-soc(i+1)) soc(i)>soc(i+1)
其中,m是储能一天的等效循环次数,soc(i)表示时刻i的储能荷电状态;
储能参与调频后的实际寿命计算公式如下:
Figure FDA0002361907770000051
其中,a是蓄电池的出厂寿命,q是蓄电池的最大循环次数,N是储能参与调频后的实际寿命;
然后对储能投资成本项中的蓄电池的出厂寿命a进行修正得到实际的成本项C'in,year
Figure FDA0002361907770000052
5.根据权利要求4所述的改善电厂调频性能的容量配置方法,其特征在于,所述步骤(41)包括:
所述的储能各项建设投资成本计算公式如下:
Figure FDA0002361907770000053
Cpcs=CpPb
Caf=CAEb
其中,CE是储能本体的单位能量价格,Cp是PCS的单位功率价格,CA是辅助设施的单位能量价格,Eb是储能的额定容量,Pb是储能的额定功率,η是储能的充放电效率;
所述的储能的固定运行成本Cfix和可变运行维护成本Cvar计算公式如下式:
Cfix=CFPb
Figure FDA0002361907770000061
其中,CF是单位功率固定运行成本,Pb表示储能的额定功率,Cc是当前的充电电价,P是储能的充电功率,η表示储能的充放电效率。
6.根据权利要求1所述的改善电厂调频性能的容量配置方法,其特征在于,所述步骤(5)包括:
储能调频的经济性模型为储能参与调频的净收益:
Figure FDA0002361907770000062
其中,C表示储能参与调频的净收益,
Figure FDA0002361907770000063
是机组参与AGC调频服务的经济补偿,Cin是储能的初始建设投资成本,COM是储能的运行维护成本。
7.根据权利要求1所述的改善电厂调频性能的容量配置方法,其特征在于,所述步骤(6)中差分进化算法的步骤如下:
(1)生成初始种群
在n1维空间里随机生成满足约束条件的h个染色体,其中,h为种群规模,n1为变量个数,
Figure FDA0002361907770000064
分别为求解变量的上下限值,xi(g)表示第i条染色体或个体的g代值,第xi,j(g)表示第i个染色体上的第j个基因的g代值,rand(0,1)表示0和1之间的随机数;在实际算例中
Figure FDA0002361907770000065
分别为取为储能容量的上下限值;初始种群xi,j(0)计算公式为:
Figure FDA0002361907770000071
(2)计算不同储能容量下储能参与调频净收益的相反数,且差分进化算法默认为求最小值,因此求最大值问题需要将目标函数前面加上负号;
(3)判断是否满足终止条件
若满足终止条件,则直接输出最优结果;若不满足终止条件,则执行步骤(4);
(4)变异操作
从群体中随机选择3个染色体
Figure FDA0002361907770000072
且(p1≠p2≠p3),其中
Figure FDA0002361907770000073
为差异化向量,η1为缩放因子,vi(g+1)为变异后的第g+1代个体计算公式为:
Figure FDA0002361907770000074
(5)交叉操作
具体操作如下:
Figure FDA0002361907770000075
其中,ui,j(g+1)为交叉操作后第i条染色体第j个基因g+1代的值,CR为交叉概率,jrand为[1,2,...,n1]中的随机数;
(6)选择操作
Figure FDA0002361907770000076
其中,i=1,2,...,h,xi(g+1)为染色体i的g+1代值,f(·)为问题的目标函数(默认最小值问题,最大值加负号),实际算例中即为储能参与电厂调频净收益的相反数;
(7)令g=g+1,返回步骤(2)继续操作,直至达到最大的迭代次数Tmax
8.根据权利要求1所述的改善电厂调频性能的容量配置方法,其特征在于,所述的储能参与调频的约束条件如下:
1)储能充放电功率约束
Figure FDA0002361907770000081
其中,Pd(t)是储能t时刻的放电功率,Pd,max是储能放电功率的最大值,Pc(t)是储能t时刻的充电功率,Pc,max是储能充电功率的最大值;
2)储能SOC约束
S1≤SOC(t)≤S2
其中,S1是时刻t的储能荷电状态的下限,S2是时刻t的储能荷电状态的上限。
CN201710636185.9A 2017-07-31 2017-07-31 一种改善电厂调频性能的容量配置方法 Active CN107579515B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710636185.9A CN107579515B (zh) 2017-07-31 2017-07-31 一种改善电厂调频性能的容量配置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710636185.9A CN107579515B (zh) 2017-07-31 2017-07-31 一种改善电厂调频性能的容量配置方法

Publications (2)

Publication Number Publication Date
CN107579515A CN107579515A (zh) 2018-01-12
CN107579515B true CN107579515B (zh) 2020-04-24

Family

ID=61034345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710636185.9A Active CN107579515B (zh) 2017-07-31 2017-07-31 一种改善电厂调频性能的容量配置方法

Country Status (1)

Country Link
CN (1) CN107579515B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108711891B (zh) * 2018-04-12 2020-10-23 广东电网有限责任公司 一种考虑机组调频性能的调频辅助服务优化系统及方法
CN108521133B (zh) * 2018-04-18 2022-02-01 国网山西省电力公司电力科学研究院 基于集合经验模态分解的火-储联合调度调频功率分配方法
CN110070387A (zh) * 2019-03-19 2019-07-30 武汉华中思能科技有限公司 一种基于边际成本的电力调频报价方法及系统
CN110112761B (zh) * 2019-04-16 2022-10-04 中国能源建设集团广东省电力设计研究院有限公司 适用于电厂储能辅助调频的储能定容方法、装置及设备
CN110266031B (zh) * 2019-07-03 2021-01-26 北京五洲驭新科技有限公司 储能并网充放电量控制方法、装置、服务器及存储介质
CN110994639B (zh) * 2019-10-24 2021-05-14 中国能源建设集团广东省电力设计研究院有限公司 一种电厂储能辅助调频的仿真定容方法、装置及设备
CN110994646B (zh) * 2019-12-11 2024-02-09 国网河南省电力公司电力科学研究院 电网agc调节的运行效果评估方法、系统及存储介质
CN111106628B (zh) * 2020-01-04 2021-08-24 深圳市中业智能系统控制有限公司 基于发电机和储能系统的调频控制方法、装置及终端设备
CN111398726B (zh) * 2020-04-30 2021-12-28 南方电网电力科技股份有限公司 一种储能调频系统充放电响应时间测试方法、系统及设备
CN112350344B (zh) * 2020-05-25 2022-03-25 清华大学 储能系统-火电机组联合调频控制方法
CN111817339B (zh) * 2020-06-04 2022-06-10 中国大唐集团科学技术研究院有限公司火力发电技术研究院 一种量化电化学储能提升agc调节性能效果的方法
CN112086975B (zh) * 2020-09-01 2021-11-26 东南大学 一种协调多储能单元参与二次调频的优化调度方法
CN112465665A (zh) * 2020-11-17 2021-03-09 西安热工研究院有限公司 改善电厂调频收益的燃气机组电储能系统容量配置的方法
CN112736954B (zh) * 2021-02-01 2023-05-02 吉林省电力科学研究院有限公司 一种提高agc调节性能的储能控制方法
CN113178883B (zh) * 2021-03-17 2022-11-08 北京天启鸿源新能源科技有限公司 一种联合储能系统的储能配置方法
CN113629768B (zh) * 2021-08-16 2023-06-20 广西大学 一种差分进化变参数矢量情感深度强化学习发电控制方法
CN115525041B (zh) * 2022-09-21 2023-07-25 华能威海发电有限责任公司 电厂性能指标考核的故障诊断方法、系统、设备及介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105207242A (zh) * 2015-09-17 2015-12-30 山东大学 储能装置参与机组调频的优化控制与容量规划系统及方法
CN106487038A (zh) * 2016-11-28 2017-03-08 北京交通大学 一种用于风电‑火电‑储能联合调频系统的控制优化方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105207242A (zh) * 2015-09-17 2015-12-30 山东大学 储能装置参与机组调频的优化控制与容量规划系统及方法
CN106487038A (zh) * 2016-11-28 2017-03-08 北京交通大学 一种用于风电‑火电‑储能联合调频系统的控制优化方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
基于风电预测误差随机性的火储联合相关机会规划调度;赵书强等;《中国电机工程学报》;20141120;第34卷;全文 *
混合储能辅助机组参与调频的优化控制与容量配置;牛阳;《中国优秀硕士学位论文全文库》;20170215;C042-2021 *
牛阳.混合储能辅助机组参与调频的优化控制与容量配置.《中国优秀硕士学位论文全文库》.2017, *
考虑储能设备的AGC竞价市场研究;任晓朦;《中国优秀硕士学位论文全文库》;20160215;C042-1228 *
采用全寿命周期成本模型的用户侧电池储能经济可行性研究;薛金花等;《电网技术》;20160831;第40卷(第8期);第2471-2476页 *

Also Published As

Publication number Publication date
CN107579515A (zh) 2018-01-12

Similar Documents

Publication Publication Date Title
CN107579515B (zh) 一种改善电厂调频性能的容量配置方法
CN105846461B (zh) 一种大规模储能电站自适应动态规划的控制方法和系统
CN110518643B (zh) 一种储能联合火电机组参与agc调频的控制方法及装置
CN112186811B (zh) 一种基于深度强化学习的agc机组动态优化方法
CN109873455B (zh) 一种储能辅助火电机组agc调频方法及系统
CN106887858B (zh) 一种接入新能源发电的储能系统跟踪计划出力方法及装置
CN108054771A (zh) 一种储能系统充放电控制方法及系统
CN110323768B (zh) 一种电化学储能电站功率分配方法及系统
CN108365637B (zh) 一种送电计划与抽水储能发电计划优化方法及系统
CN114256882B (zh) 一种基于风险控制的风储联合系统调度方法及系统
CN114820046A (zh) 一种区域电网混合储能辅助调频经济优化与补偿定价方法
CN115189367A (zh) 一种多时间尺度下的风光储蓄联合优化调度方法及装置
CN114358378A (zh) 一种用于考虑需量管理的用户侧储能优化配置系统及方法
CN113659623A (zh) 一种基于布林线理论的风储联合系统的优化方法及系统
CN110994639B (zh) 一种电厂储能辅助调频的仿真定容方法、装置及设备
CN112531735A (zh) 基于机器学习的自动发电控制系统的功率分配方法及装置
CN110890767A (zh) 一种多抽水蓄能电厂均衡调度方法
CN111404182B (zh) 基于磷酸铁锂介质的电化学储能电站的运行方法
CN113300388B (zh) 一种基于改进旋转门算法的风电功率波动平抑方法
CN110311424B (zh) 一种基于双时间尺度净负荷预测的储能调峰控制方法
CN112465665A (zh) 改善电厂调频收益的燃气机组电储能系统容量配置的方法
CN114725961A (zh) 平抑风电波动的制氢系统容量分层优化配置方法
CN109636096B (zh) 一种储能电站输出功率优化方法及装置
CN110233486B (zh) 基于主动预测soc的储能辅助模糊pid控制方法
CN111030149A (zh) 一种基于复合储能装置的微电网供电方法及电子设备

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