CN104979850B - 一种储能参与备用的含风电的电力***调度方法 - Google Patents

一种储能参与备用的含风电的电力***调度方法 Download PDF

Info

Publication number
CN104979850B
CN104979850B CN201510381736.2A CN201510381736A CN104979850B CN 104979850 B CN104979850 B CN 104979850B CN 201510381736 A CN201510381736 A CN 201510381736A CN 104979850 B CN104979850 B CN 104979850B
Authority
CN
China
Prior art keywords
mrow
msubsup
power
msub
energy storage
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
CN201510381736.2A
Other languages
English (en)
Other versions
CN104979850A (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.)
State Grid Corp of China SGCC
Economic and Technological Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Economic and Technological Research Institute of State Grid Shandong Electric Power Co Ltd
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 State Grid Corp of China SGCC, Economic and Technological Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510381736.2A priority Critical patent/CN104979850B/zh
Publication of CN104979850A publication Critical patent/CN104979850A/zh
Application granted granted Critical
Publication of CN104979850B publication Critical patent/CN104979850B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects
    • 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
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种储能参与备用的含风电的电力***调度方法,本方法包括以下步骤:实时监测火电机组功率、爬坡速率和主动负荷的功率;对风电历史预测误差进行统计,生成风电场景;依据***的约束条件,建立调度模型,求解该模型;统计调度结果。在备用配置中考虑了储能***作用,储能***能够通过改变充/放电功率为应对风电不确定性提供备用,在二次、一次调整中起到消除风电、负荷不确定性的作用,同时考虑了频率调节效应,最终能够减轻常规机组备用配置的压力,提高电网消纳风电的能力。

Description

一种储能参与备用的含风电的电力***调度方法
技术领域
本发明涉及一种储能参与备用的含风电的电力***调度方法。
背景技术
生态环境日益恶化,威胁着人类生存,能源清洁化变革势在必行,电力***必须向绿色、低碳、环保和节能的智能化方向发展。在此背景下,风光等可再生能源发电得到了迅猛发展,使电力***运行过程中必须面对强不确定性接纳能力的考验。
专利号为201110001574.7的中国专利:“消纳风电接入的大电网实时调度方法”,给出了一种消纳风电接入的大电网实时调度方法,该专利通过将全网机组分类,然后从***获得机组的日前计划、实时出力、联络线计划及数值天气预报等相关信息,进行超短期负荷预测及风电出力预测,求得下一时刻实时调度机组的出力调节量并构建弃风损失最小的有功实时调度模型,采用单纯形法求得包括风电机组的实时调度机组出力调整量。
专利号为201210371334.0的中国专利:“基于风电出力不确定性集合的电力***调度方法”,公开了一种基于风电出力不确定性集合的电力***调度方法,解决了含大规模风电的电力***鲁棒调度问题中,风电出力不确定性建模的难题,保证电网机组启停及调度计划的安全可靠性,有效提高风电场风电出力的利用率。
专利号为201210176522.8的中国专利:“一种电动汽车与风电协同实时调度优化方法”,根据实际电网的电网模型建立以有功总电力最大为目标的优化模型;将总有功出力曲线变动关系的惩罚量加入到优化目标中,获得考虑发电曲线平滑的优化调度模型;将模型中的非线性因素线性化,采用对偶单纯形法求解,得出风光储联合发电***的有功发电曲线,上报至上级调度中心,并得到储能装置的充放电计划,下发子***执行。
以上专利所述的调度方法均是在调度中锁定储能***充/放电功率,完全由常规机组承担应对不确定性的备用。在常规机组主导地位逐渐削弱的态势下,这一研究思路将难以奏效。
发明内容
本发明为了解决上述问题,提出了一种储能参与备用的含风电的电力***调度方法,本方法在备用配置中考虑了储能***作用,储能***能够通过改变充/放电功率为应对风电不确定性提供备用,在二次、一次调整中起到消除风电、负荷不确定性的作用,同时考虑了频率调节效应,最终能够减轻常规机组备用配置的压力,提高电网消纳风电的能力。
为了实现上述目的,本发明采用如下技术方案:
一种储能参与备用的含风电的电力***调度方法,包括以下步骤:
(1)实时监测火电机组功率、爬坡速率和主动负荷的功率;
(2)对接入电网的风电场历史预测误差进行统计,生成风电场景;
(3)依据***的约束条件,建立电力***调度模型,求解该模型;
所述步骤(1)中,实时监测的数据包括:机组最大输出功率、最小输出功率、机组爬坡速率和主动负荷最大充/放电功率。
所述步骤(2)中,采用自回归滑动平均(ARMA)模型来对风电场输出功率预测误差进行估计,表达式为,
式中,p、q分别为ARMA模型中的阶数;αi、βj为模型参数,通过估计得到;ε为服从均值为0,方差为σ2的高斯分布的白噪声;为预测误差。
通过对风电历史预测误差进行统计,采用最小二乘方法获得式(1)中参数的取值。由于ε服从正态分布,随机模拟产生ε,可生成需要的风电场景。
所述步骤(3)中,电力***调度模型为:
考虑到风电不消耗燃料,不计储能***的成本费用和运行费用,目标函数表示成如下形式:
式中,PW,t为t时段风电场输出功率预测值;PL,t为t时段负荷预测值;PESS,t为t时段储能***有功功率,PESS,t>0表示储能***充电,PESS,t<0表示储能***放电。
所述步骤(3)中,需要考虑约束条件,约束条件包括:
功率平衡约束:
式中,PW,t为t时段风电场输出功率预测值;PL,t为t时段负荷预测值;PESS,t为t时段储能***有功功率,PESS,t>0表示储能***充电,PESS,t<0表示储能***放电。
所述步骤(3)中,约束条件包括:
机组输出功率上下限约束:
式中,分别为机组i输出功率的下限和上限。
所述步骤(3)中,约束条件包括:
机组爬坡速率约束:
-rGi·Δt≤PGi,t+1-PGi,t≤rGi·Δt (5)
式中,rGi为机组i输出功率最大调整速率,Δt为时段长度。
所述步骤(3)中,约束条件包括:储能***逆变器充放电功率约束:
式中,为逆变器允许的充放电功率最大值。
所述步骤(3)中,约束条件包括:各时段末储能***存储能量约束,
Et=Et-1+PESS,tΔt
所述步骤(3)中,约束条件包括:储能***存储能量约束,
Emin≤Et≤Emax (7)
式中,Emin和Emax分别为储能***储存能量的最小值和最大值。
所述步骤(3)中,约束条件包括,基于各个场景的约束:设为场景s中风电在时段t的输出功率,为场景s中t时段负荷值,当PESS,t>0,即储能***充电时,储能***可通过增加充电功率应对风电输出功率高于计划值或负荷低于预测值的场景,储能***可锁定的波动范围受逆变器允许的充电功率和最大存储能量的限制,具体可表示为
当PESS,t>0,即储能***充电时,储能***可通过减少充电功率或放电应对风电输出功率小于计划值或负荷高于预测值的场景,可锁定的范围受逆变器允许的放电功率和最小存储能量的限制,具体可表示为
对于场景s,需火电机组提供的向上和向下调节备用容量可分别表示为,
火电机组i在t时段向上和向下调节备用容量分别表示为考虑到***频率偏差为时负荷的频率调节效应,在风电和负荷实际功率偏离期望值时,若允许频率在一定范围内变动,应满足,
式中,D为单位频率变化所引起负荷有功功率增加或减少的系数;火电机组向上和向下调节备用容量应满足,
式中,RGi单位频率变化所引起机组i输出功率增加或减少的系数;
同时,火电机组输出功率应满足上下限约束,即
***频率偏差应满足,
式中,Δfmin、Δfmax分别为***允许的最小和最大频率偏差。
所述步骤(3)中,依据目标函数中成本函数取二次或线性表达,可分别利用二次规划或线性规划算法求解。
本发明的有益效果为:
(1)本发明在备用配置中考虑了储能***作用,储能***能够通过改变充/放电功率为应对风电不确定性提供备用,在二次、一次调整中起到消除风电、负荷不确定性的作用;
(2)本发明考虑了频率调节效应,最终能够减轻常规机组备用配置的压力,提高电网消纳风电的能力。
附图说明
图1为负荷、风电功率曲线图。
具体实施方式:
下面结合附图与实施例对本发明作进一步说明。
一种储能***参与备用的含风电的电力***调度方法,包括以下步骤:
1)获得火电机组、主动负荷相关参数和信息。火电机组最大调节速率为每分钟机组容量1%,火电机组调差系数为4%,即***频率变化4%引起机组输出功率100%的变化。储能***参数见表1。
表1储能***参数
2)根据风电历史预测误差的统计,生成风电场景。风电与负荷在前瞻时段内的期望值如图1所示。
3)建立调度模型,并对模型进行求解;
调度模型为:
考虑到风电不消耗燃料,不计储能***的成本费用和运行费用,目标函数表示成如下形式:
式中,PW,t为t时段风电场输出功率预测值;PL,t为t时段负荷预测值;PESS,t为t时段储能***有功功率,PESS,t>0表示储能***充电,PESS,t<0表示储能***放电。
约束条件包括:
1)功率平衡约束
式中,PW,t为t时段风电场输出功率预测值;PL,t为t时段负荷预测值;PESS,t为t时段储能***有功功率,PESS,t>0表示储能***充电,PESS,t<0表示储能***放电。
2)机组输出功率上下限约束
式中,分别为机组i输出功率的下限和上限。
3)机组爬坡速率约束
-rGi·Δt≤PGi,t+1-PGi,t≤rGi·Δt (4)
式中,rGi为机组i输出功率最大调整速率,Δt为时段长度。
4)储能***逆变器充放电功率约束
式中,为逆变器允许的充放电功率最大值。
5)各时段末储能***存储能量约束
Et=Et-1+PESS,tΔt (6)
6)储能***存储能量约束
Emin≤Et≤Emax (7)
式中,Emin和Emax分别为储能***储存能量的最小值和最大值。
7)基于各个场景的约束
通过可能出现的场景模拟风电和负荷的不确定性。假设为场景s中风电在时段t的输出功率,为场景s中t时段负荷值。实际运行中,若风电场输出功率偏离计划值或负荷偏离预测值,储能***可在二次和一次调整中改变输出功率,将风电和负荷尽量锁定为计划值。
当PESS,t>0,即储能***充电时,储能***可通过增加充电功率应对风电输出功率高于计划值或负荷低于预测值的场景,储能***可锁定的波动范围受逆变器允许的充电功率和最大存储能量的限制,具体可表示为
当PESS,t>0,即储能***充电时,储能***可通过减少充电功率或放电应对风电输出功率小于计划值或负荷高于预测值的场景,可锁定的范围受逆变器允许的放电功率和最小存储能量的限制,具体可表示为
对于场景s,需火电机组提供的向上和向下调节备用容量可分别表示为,
火电机组i在t时段向上和向下调节备用容量分别表示为考虑到***频率偏差为时负荷的频率调节效应,在风电和负荷实际功率偏离期望值时,若允许频率在一定范围内变动,应满足,
式中,D为单位频率变化所引起负荷有功功率增加或减少的系数。
火电机组向上和向下调节备用容量应满足,
式中,RGi单位频率变化所引起机组i输出功率增加或减少的系数。
火电机组输出功率应满足上下限约束,即
***频率偏差应满足,
式中,Δfmin、Δfmax分别为***允许的最小和最大频率偏差。
式(1)—(14)构成了基本的调度模型,依据目标函数中成本函数取二次或线性表达,可分别利用二次规划或线性规划算法求解。
4)统计调度结果。
假设风电预测误差以固定步长10%逐步增加,对应场景中风电波动性增强,下表给出了本文方法和传统调度方法的结果对比情况
表2风电不确定性增强时两种方法对比情况
由表2可知,随着风电不确定性增强,两种方法对应的成本都会增加。风电预测误差增加相同的量,本文方法成本的增加量要小于传统方法。并且,当对应的预测误差达到40%,传统方法无解,本文方法由于考虑了储能参与备用,可以应对风电更大范围的波动,增加了电网消纳可再生能源发电的能力。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (8)

1.一种储能参与备用的含风电的电力***调度方法,其特征是:包括以下步骤:
(1)实时监测火电机组功率、爬坡速率和主动负荷的功率;
(2)对接入电网的风电场历史预测误差进行统计,生成风电场景;
(3)依据***的约束条件,建立电力***调度模型,求解该模型;
所述步骤(2)中,采用自回归滑动平均(ARMA)模型来对风电场输出功率预测误差进行估计,表达式为,
<mrow> <msubsup> <mi>&amp;Delta;p</mi> <mrow> <mi>g</mi> <mi>w</mi> </mrow> <mi>t</mi> </msubsup> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>p</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>i</mi> </msub> <mo>&amp;CenterDot;</mo> <msubsup> <mi>&amp;Delta;p</mi> <mrow> <mi>g</mi> <mi>w</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mi>i</mi> </mrow> </msubsup> <mo>+</mo> <msub> <mi>&amp;epsiv;</mi> <mi>t</mi> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>q</mi> </munderover> <msub> <mi>&amp;beta;</mi> <mi>j</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>t</mi> <mo>-</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式中,p、q分别为ARMA模型中的阶数;αi、βj为模型参数,通过估计得到;ε为服从均值为0,方差为σ2的高斯分布的白噪声;为预测误差;
通过对风电历史预测误差进行统计,采用最小二乘方法获得式(1)中参数的取值;由于ε服从正态分布,随机模拟产生ε,可生成需要的风电场景;
所述步骤(3)中,约束条件包括,基于各个场景的约束:设为场景s中风电在时段t的输出功率,为场景s中t时段负荷值,当PESS,t>0,即储能***充电时,储能***可通过增加充电功率应对风电输出功率高于计划值或负荷低于预测值的场景,储能***可锁定的波动范围受逆变器允许的充电功率和最大存储能量的限制,具体可表示为
<mrow> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>d</mi> <mi>n</mi> </mrow> </msubsup> <mo>=</mo> <mi>min</mi> <mo>{</mo> <mfrac> <mrow> <msup> <mi>E</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msup> <mo>-</mo> <msub> <mi>E</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> <mo>,</mo> <msubsup> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> </mrow> <mi>max</mi> </msubsup> <mo>}</mo> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式中,Emin和Emax分别为储能***储存能量的最小值和最大值;Et-1为t-1时段存储***存储能量;为逆变器允许的充放电功率最大值,PESS,t为t时段储能***有功功率;
当PESS,t>0,即储能***充电时,储能***可通过减少充电功率或放电应对风电输出功率小于计划值或负荷高于预测值的场景,可锁定的范围受逆变器允许的放电功率和最小存储能量的限制,具体可表示为
<mrow> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> </mrow> </msubsup> <mo>=</mo> <mi>min</mi> <mo>{</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msup> <mi>E</mi> <mi>min</mi> </msup> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> <mo>,</mo> <msubsup> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> </mrow> <mi>max</mi> </msubsup> <mo>}</mo> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中,Emin和Emax分别为储能***储存能量的最小值和最大值;Et-1为t-1时段存储***存储能量;为逆变器允许的充放电功率最大值,PESS,t为t时段储能***有功功率;
对于场景s,需火电机组提供的向上和向下调节备用容量可分别表示为,
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>RES</mi> <mi>t</mi> <mrow> <mi>u</mi> <mi>p</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <mi>max</mi> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>W</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>W</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>s</mi> </msubsup> <mo>+</mo> <msubsup> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> </mrow> </msubsup> <mo>}</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>RES</mi> <mi>t</mi> <mrow> <mi>d</mi> <mi>n</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <mi>max</mi> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mrow> <mi>W</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>W</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>s</mi> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>d</mi> <mi>n</mi> </mrow> </msubsup> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
火电机组i在t时段向上和向下调节备用容量分别表示为PW,t为t时段风电场输出功率预测值;PL,t为t时段负荷预测值;rGi为机组i输出功率最大调整速率;Δt为时段长度;考虑到***频率偏差为Δft s时负荷的频率调节效应,在风电和负荷实际功率偏离期望值时,若允许频率在一定范围内变动,应满足,
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mi>g</mi> </mrow> </munderover> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>D&amp;Delta;f</mi> <mi>t</mi> <mi>s</mi> </msubsup> <mo>&amp;GreaterEqual;</mo> <msubsup> <mi>RES</mi> <mi>t</mi> <mrow> <mi>u</mi> <mi>p</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mi>g</mi> </mrow> </munderover> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>d</mi> <mi>n</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>D&amp;Delta;f</mi> <mi>t</mi> <mi>s</mi> </msubsup> <mo>&amp;GreaterEqual;</mo> <msubsup> <mi>RES</mi> <mi>t</mi> <mrow> <mi>d</mi> <mi>n</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
式中,D为单位频率变化所引起负荷有功功率增加或减少的系数;Ng为发电机的总数量;火电机组向上和向下调节备用容量应满足,
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msub> <mi>r</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>+</mo> <msub> <mi>R</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> </msub> <msubsup> <mi>&amp;Delta;f</mi> <mi>t</mi> <mi>s</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>d</mi> <mi>n</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msub> <mi>r</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>+</mo> <msub> <mi>R</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> </msub> <msubsup> <mi>&amp;Delta;f</mi> <mi>t</mi> <mi>s</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
式中,RGi单位频率变化所引起机组i输出功率增加或减少的系数;
同时,火电机组输出功率应满足上下限约束,即
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>max</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>d</mi> <mi>n</mi> <mo>,</mo> <mi>s</mi> </mrow> </msubsup> <mo>&amp;GreaterEqual;</mo> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>min</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
其中,分别为机组i输出功率的下限和上限,PGi,t为火电机组i在t时刻输出功率;
***频率偏差应满足,
Δfmin≤Δft s≤Δfmax (8)
式中,Δfmin、Δfmax分别为***允许的最小和最大频率偏差。
2.如权利要求1所述的一种储能参与备用的含风电的电力***调度方法,其特征是:所述步骤(1)中,实时监测的数据包括:机组最大输出功率、最小输出功率、机组爬坡速率和主动负荷最大充/放电功率。
3.如权利要求1所述的一种储能参与备用的含风电的电力***调度方法,其特征是:所述步骤(3)中,电力***调度模型为:
考虑到风电不消耗燃料,不计储能***的成本费用和运行费用,目标函数表示成如下形式:
<mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>T</mi> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mi>g</mi> </mrow> </munderover> <msub> <mi>C</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>C</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mrow> <mi>r</mi> <mi>e</mi> <mi>s</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;Delta;P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>d</mi> <mi>n</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式中,PW,t为t时段风电场输出功率预测值;PL,t为t时段负荷预测值;PESS,t为t时段储能***有功功率,PESS,t>0表示储能***充电,PESS,t<0表示储能***放电;PGi,t为火电机组i在t时刻输出功率,火电机组i在t时段向上和向下调节备用容量;CGi为火电机组i的发电量成本,为火电机组i的备用容量成本;Ng为发电机的总数量。
4.如权利要求1所述的一种储能参与备用的含风电的电力***调度方法,其特征是:所述步骤(3)中,需要考虑约束条件,约束条件包括:
功率平衡约束:
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mi>g</mi> </mrow> </munderover> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mi>W</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
式中,PW,t为t时段风电场输出功率预测值;PL,t为t时段负荷预测值;PESS,t为t时段储能***有功功率,PESS,t>0表示储能***充电,PESS,t<0表示储能***放电;PGi,t为火电机组i在t时刻输出功率;Ng为发电机的总数量;
机组输出功率上下限约束:
<mrow> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>max</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
式中,分别为机组i输出功率的下限和上限,PGi,t为火电机组i在t时刻输出功率。
5.如权利要求1所述的一种储能参与备用的含风电的电力***调度方法,其特征是:所述步骤(3)中,约束条件包括:
机组爬坡速率约束:
-rGi·Δt≤PGi,t+1-PGi,t≤rGi·Δt (12)
式中,rGi为火电机组i输出功率最大调整速率,Δt为时段长度;PGi,t为火电机组i在t时刻输出功率,PGi,t+1为火电机组i在t+1时刻输出功率。
6.如权利要求1所述的一种储能参与备用的含风电的电力***调度方法,其特征是:所述步骤(3)中,约束条件包括:储能***逆变器充放电功率约束:
<mrow> <mo>-</mo> <msubsup> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> </mrow> <mi>max</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>E</mi> <mi>S</mi> <mi>S</mi> </mrow> <mi>max</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
式中,为逆变器允许的充放电功率最大值。
7.如权利要求1所述的一种储能参与备用的含风电的电力***调度方法,其特征是:所述步骤(3)中,约束条件包括:各时段末储能***存储能量约束,
Et=Et-1+PESS,tΔt
式中,PESS,t为t时段储能***有功功率,Δt为时段长度;
所述步骤(3)中,约束条件包括:储能***存储能量约束,
Emin≤Et≤Emax (14)
式中,Emin和Emax分别为储能***储存能量的最小值和最大值。
8.如权利要求1所述的一种储能参与备用的含风电的电力***调度方法,其特征是:所述步骤(3)中,依据目标函数中成本函数取二次或线性表达,分别利用二次规划或线性规划算法求解。
CN201510381736.2A 2015-07-01 2015-07-01 一种储能参与备用的含风电的电力***调度方法 Active CN104979850B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510381736.2A CN104979850B (zh) 2015-07-01 2015-07-01 一种储能参与备用的含风电的电力***调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510381736.2A CN104979850B (zh) 2015-07-01 2015-07-01 一种储能参与备用的含风电的电力***调度方法

Publications (2)

Publication Number Publication Date
CN104979850A CN104979850A (zh) 2015-10-14
CN104979850B true CN104979850B (zh) 2017-09-26

Family

ID=54276086

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510381736.2A Active CN104979850B (zh) 2015-07-01 2015-07-01 一种储能参与备用的含风电的电力***调度方法

Country Status (1)

Country Link
CN (1) CN104979850B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107203855B (zh) * 2017-08-03 2020-07-07 国网江苏省电力公司宿迁供电公司 含风电***实时调度的鲁棒两层优化模型及转化等效方法
CN109742779B (zh) * 2018-12-14 2022-09-02 中国电力科学研究院有限公司 一种分布式储能调度方法和装置
CN110348637B (zh) * 2019-07-12 2022-10-11 哈尔滨工业大学 一种考虑场-网因素的风电爬坡事件预警方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103927695B (zh) * 2014-04-22 2017-11-24 国家电网公司 基于自学习复合数据源的风电功率超短期预测方法
CN104166946B (zh) * 2014-08-15 2018-01-19 国家电网公司 促进新能源并网消纳的备用及调峰辅助服务成本分摊方法
CN104167762B (zh) * 2014-08-18 2016-08-03 国家电网公司 一种含化学储能***和间歇式电源的机组组合优化方法
CN104300585B (zh) * 2014-10-13 2016-12-07 江苏大学 一种基于大规模风电并网的经济调度优化方法
CN104734200B (zh) * 2015-03-26 2017-03-01 国家电网公司 一种基于虚拟发电的主动配电网优化调度方法

Also Published As

Publication number Publication date
CN104979850A (zh) 2015-10-14

Similar Documents

Publication Publication Date Title
WO2017000853A1 (zh) 主动配电网多时间尺度协调优化调度方法和存储介质
CN104362673B (zh) 基于调峰裕度的风电并网协调调度优化方法
CN108092324B (zh) 一种风电参与调峰调频的agc控制***和控制方法
CN112615399A (zh) 一种储能***参与电网调频优化控制方法、***及存储介质
CN105262117B (zh) 一种风电爬坡优化控制方法
CN105303267B (zh) 一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优化方法
CN104467029B (zh) 一种区域小水电、风电日内联合优化调度的方法
CN111737846B (zh) 一种水电厂参与电力现货市场出清计算的方法
CN104037805B (zh) 一种计及电网安全约束的光伏电站可发电裕度分配方法
CN104348190B (zh) 风电场群汇集***的有功功率调节方法
CN113644670B (zh) 一种储能容量优化配置的方法及***
CN112053035B (zh) 考虑经济性与灵活性的输电通道与储能联合规划方法
CN106026198B (zh) “风-水”打捆外送的agc协调控制***及控制方法
CN104979850B (zh) 一种储能参与备用的含风电的电力***调度方法
CN113241803A (zh) 一种基于新能源消纳的储能调度方法及计算机介质
CN105826946B (zh) 大规模光伏接入的配电网动态无功优化方法
CN114154790A (zh) 基于需量管理和柔性负荷的工业园区光储容量配置方法
CN115659096A (zh) 考虑源荷不确定性的微网多时间尺度能量调度方法及装置
CN105305485A (zh) 一种消纳大规模间歇性能源的安全约束经济调度方法
CN109687423B (zh) 一种直流电网的储能配置方法
CN109829624B (zh) 一种风电合作博弈爬坡控制方法和装置
CN107482657A (zh) 风电爬坡率实时平抑方法和***
CN110336308B (zh) 一种基于机会约束的主动配电网经济调度方法
CN117394404A (zh) 考虑碳效益和辅助调频的风电场储能容量配置方法
CN109617100B (zh) 一种数据驱动的风电场储能容量规划方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhang Jie

Inventor after: Mou Hong

Inventor after: Wang Yuan

Inventor after: Cao Xiangyang

Inventor after: Wang Nan

Inventor after: Shi Hongji

Inventor after: Yu Baijian

Inventor after: Yang Si

Inventor after: Zhao Linghan

Inventor after: Zhang Youquan

Inventor after: Liu Guojing

Inventor after: Sun Donglei

Inventor after: Wang Mingqiang

Inventor after: Han Xueshan

Inventor after: Liu Xiaoming

Inventor after: Wang Yan

Inventor after: Wang Fei

Inventor before: Zhang Jie

Inventor before: Wang Nan

Inventor before: Niu Xinsheng

Inventor before: Liu Guojing

Inventor before: Wang Mingqiang

Inventor before: Han Xueshan

Inventor before: Wang Yan

Inventor before: Wang Fei

Inventor before: Liu Xiaoming

Inventor before: Cao Xiangyang

GR01 Patent grant
GR01 Patent grant