CN110071530B - 基于lstm的含储能电力系统风电爬坡协调调度方法 - Google Patents
基于lstm的含储能电力系统风电爬坡协调调度方法 Download PDFInfo
- Publication number
- CN110071530B CN110071530B CN201910418323.5A CN201910418323A CN110071530B CN 110071530 B CN110071530 B CN 110071530B CN 201910418323 A CN201910418323 A CN 201910418323A CN 110071530 B CN110071530 B CN 110071530B
- Authority
- CN
- China
- Prior art keywords
- wind power
- climbing
- unit
- energy storage
- power
- 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
- 230000009194 climbing Effects 0.000 title claims abstract description 169
- 238000000034 method Methods 0.000 title claims abstract description 88
- 238000004146 energy storage Methods 0.000 title claims abstract description 70
- 238000013528 artificial neural network Methods 0.000 claims abstract description 49
- 230000033228 biological regulation Effects 0.000 claims abstract description 22
- 239000003245 coal Substances 0.000 claims description 9
- 238000007599 discharging Methods 0.000 claims description 8
- 230000003828 downregulation Effects 0.000 claims description 7
- 230000003827 upregulation Effects 0.000 claims description 7
- 230000000295 complement effect Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000011056 performance test Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 230000007774 longterm Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000009472 formulation Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 230000010485 coping Effects 0.000 description 2
- 238000005034 decoration Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000013076 uncertainty analysis Methods 0.000 description 1
Images
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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H02J3/386—
-
- 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
本发明的公开了一种基于LSTM的含储能电力系统风电爬坡协调调度方法,首先,利用LSTM神经网络对风功率爬坡事件进行提前多步预测;然后,提出一种应对风电爬坡事件的提前‑滞后调节方法,利用储能系统进行提前滞后两步补偿,降低风功率爬坡时的速率,同时降低了对储能功率和容量的需求;基于上述预测和调节方法,建立了考虑风电爬坡事件的火储协调调度模型,提前协调机组,使机组在爬坡时有足够的可调备用容量,通过对不同的风电爬坡条件的判断协调机组和储能,实现更经济和安全的爬坡调度方案。通过对预测方法的性能测试和爬坡调度仿真,验证了预测方法、调节方法和调度模型的有效性。
Description
技术领域
本发明涉及一种基于LSTM的含储能电力系统风电爬坡协调调度方法,属风电的不确定性分析和新能源并网调度领域。
背景技术
随着能源短缺和环境污染的加剧,以风能为代表的清洁能源得到广泛关注,风电并网容量也在日益提高。然而,风电具有较强的随机性和波动性,并网后会给电力系统的安全运行带来风险。特别是发生风电爬坡事件,即风功率幅值在短时间内大幅度波动现象。因其发生的概率小而难以预测,一旦发生,会使电力系统无法制定及时有效的调整措施,给电力系统的稳定运行和安全供电带来严重影响。因此,在含风电并网的电力系统运行调度中,考虑风电爬坡事件的影响,对其进行提前预测和制定调整策略至关重要。
关于含风电的电力系统优化调度问题,国内外已有了很多研究,如概率模型、场景法、模糊理论、随机规划等,但针对风电爬坡事件的调度方法研究较少,目前的方法大多由上述理论发展而来,主要有以下几类。一是基于预测值的随机规划模型;二是构建极限场景集,建立应对最坏爬坡场景的鲁棒优化模型;三是基于博弈论;四是协调储能系统参与协调调度。上述方法能够在一定程度上应对风电爬坡事件,但均有局限性。风电爬坡事件为小概率事件,不确定性参数难以用概率方法准确描述;场景理论依赖于预测值的准确性,否则会造成结果的较大偏离;博弈论仅依靠对风电场之间的协调调节应对风电爬坡事件,没有考虑机组和储能系统的协调,故调节能力有限;利用储能系统的调整依赖于合理调节策略的制定和对风功率爬坡的准确预测,若调节不当会增加调度成本。
对风电爬坡事件的预测,要兼顾时效性和准确性。例如基于自回归对数模型对风概率爬坡事件进行概率预测,同时估计各种阈值下的风电爬坡事件的概率;例如采用NWP模型,基于风电场景集进行爬坡预测。上述两种方法只能实现1~2步较准确的预测,随着预测步长的增加,误差增大,无法满足爬坡预测要求。基于此,相关学者提出基于EEMD去噪和改进RSPA的风电功率预测模型,能够对风功率进行4h的多步预测,但在预测风功率爬坡事件时误差仍较大。
发明内容
发明目的:针对上述技术,提出一种基于LSTM的含储能电力系统风电爬坡协调调度方法,提高风电爬坡的准确性,并减少或避免弃风现象。
技术方案:基于LSTM的含储能电力系统风电爬坡协调调度方法,包括以下步骤:
步骤1:利用LSTM神经网络对风功率爬坡事件进行提前多步预测;
步骤2:利用储能系统进行提前-滞后两步补偿,降低风功率爬坡时的速率,同时降低对储能功率和容量的需求;
步骤3:基于步骤1的预测结果和步骤2的调节方法,使用考虑风电爬坡事件的火储协调调度模型,提前协调机组。
进一步的,所述步骤2中包括如下具体步骤:设风电爬坡事件的开始时间为t0,持续n=4个时段,持续时间为ΔT=4×Δt,Δt为每个时段的时长;当发生风电上爬坡事件,在时段(t0-Δt,t0+2Δt)内,提前风电爬坡开始时间Δt使机组下爬坡,使机组以最大爬坡率调节,所述时段(t0-Δt,t0+2Δt)内由储能系统放电补足功率差额;在时段(t0+2Δt,t0+5Δt)内,滞后风电爬坡结束时间Δt使机组停止爬坡,所述时段(t0+2Δt,t0+5Δt)内由储能系统充电补足功率差额;当风电发生下爬坡事件时的调节方法与发生风电上爬坡事件时的调节方向相反。
进一步的,所述步骤3中考虑风电爬坡事件的火储协调调度模型建立方法包括如下步骤:
考虑机组煤耗成本和储能使用成本,得到调度时段总成本f最小为目标:
minf=CG+CB+CWC+CLS
其中
式中:CG、CB、CWC、CLS为机组煤耗成本、储能运行成本、弃风成本和切负荷成本;T为调度总时间;N为火电机组总数;ai、bi、ci为机组i煤耗成本系数;Pi,t为机组i在时段t的出力;Pt B,ch、Pt B,dc为储能在时段t的充、放电功率;kB为储能单位容量的使用成本系数;Pt WC、kWC为弃风功率和单位弃风功率惩罚系数;Pt LS、kLS为切负荷功率和单位切负荷功率惩罚系数;
约束条件包括:
1)机组和储能出力约束:
式中:Pi为机组i瞬时出力,Pi max、Pi min为机组i最大、最小技术出力;PB,max为储能最大充放电功率;
2)功率平衡约束:
式中:Pt f、Pt L为风电功率和负荷功率预测值;
3)机组爬坡约束:
式中:ri u、ri d为机组i单位时间的爬坡率;Δt为调度单时段持续时间;Pi,t-1为机组i在t-1时段的出力;
4)旋转备用约束:
5)储能容量约束:
6)爬坡对策判断条件:
式中:ΔPt g,u、ΔPt g,d为所有机组在时段t可提供的上调和下调备用;ΔPt w,u、ΔPt w ,d为风电在t时刻的上、下爬坡率;ΔPt L,u、ΔPt L,d为负荷在t时刻的上、下爬坡率;
若机组最大爬坡率大于风电和负荷爬坡率之和时,通过提前协调机组应对风电爬坡,此时无需使用储能系统,模型中Pt B,ch=Pt B,dc=0;反之,则需要进行提前弃风或使用储能系统参与协调调度,按步骤2提前-滞后调节法和所述火储协调调度模型进行优化求解。
有益效果:本发明的一种基于LSTM的含储能电力系统风电爬坡协调调度方法,首先,利用LSTM神经网络对风功率爬坡事件进行提前多步预测;然后,提出一种应对风电爬坡事件的提前-滞后调节方法,利用储能系统进行提前滞后两步补偿,降低风功率爬坡时的速率,同时降低了对储能功率和容量的需求;基于上述预测和调节方法,建立了考虑风电爬坡事件的火储协调调度模型,提前协调机组,使机组在爬坡时有足够的可调备用容量,通过对不同的风电爬坡条件的判断协调机组和储能,实现更经济和安全的爬坡调度方案。通过对预测方法的性能测试和爬坡调度仿真,验证了预测方法、调节方法和调度模型的有效性。
附图说明
图1是LSTM神经网络子模块结构图;
图2是风电上爬坡事件简图;
图3是风电上爬坡事件时的“弃风调节法”和“补偿弃风法”示意图;
图4是风电上爬坡事件时的“提前-滞后”调节方法示意图;
图5是算例2的风功率预测值曲线和实际值曲线;
图6是弃风调节法和补偿弃风法的调节过程示意图;
图7是提前-滞后调节法的调节过程示意图;
图8是算例3的预测值曲线和实际值曲线;
图9是GR神经网络风功率调度示意图;
图10是LSTM神经网络风功率调度示意图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
一种基于LSTM神经网络的含储能电力系统风电爬坡协调调度方法,包括以下步骤:
步骤1:利用LSTM神经网络对风功率爬坡事件进行提前多步预测,具体原理如下:
长短期记忆(Long-Short Term Memory,LSTM)神经网络是一种特殊的时间递归神经网络(RNN),它解决了传统神经网络无法做到信息的时序关联问题,且具有独特的记忆和遗忘模式,对具有时序相关性特征的风功率预测效果有很大的提高;同时,LSTM神经网络解决了RNN在BPTT训练中出现的长期依赖(Long-Term Dependencies)问题,能够关联到间隔时间较长的信息,从而充分利用历史信息对风功率做到更准确的预测,尤其在风功率多步提前预测方面表现明显。
LSTM神经网络子模块由输入层、隐藏层和输出层构成。与RNN不同的是,LSTM的隐藏层不是由单一的神经网络层构成,而是由4个层以一种特殊的方式交互。LSTM子模块的关键是元胞状态(cell state),即图1中顶端的水平直线。它是LSTM的记忆单元,通过门(gates)的控制有选择性的储存信息。LSTM子模块的工作步骤如下:第一,输入信息ht-1和xt合并后通过σ函数层后输出,称此函数层为遗忘门(forget gate),它决定哪些信息需要被丢弃,如式(1);第二,输入信息通过被称为输入门(input gate)的σ层和一个tanh层,决定哪些信息通过并产生一个候选值向量,如式(2)、(3);第三,旧的元胞状态Ct-1通过遗忘门和输入门更新其状态Ct,如式(4);第四,由被称为输出门(output gate)的σ层决定元胞状态的哪一部分需被输出,并通过一个tanh层进行规范化,得到输出ht,如式(5)、(6)。
ft=σ(Wf·[ht-1,xt]+bf) (1)
it=σ(Wi·[ht-1,xt]+bi) (2)
ot=σ(Wo[ht-1,xt]+bo) (5)
ht=ot*tanh(Ct) (6)
式中:ft、it、ot分别为遗忘门、输入门、输出门的状态结果;Wf、Wi、WC、Wo分别为遗忘门、输入门、输入元胞状态、输出门的权重矩阵;bf、bi、bo分别为遗忘门、输入门、输出门的偏置项;为可被加入元胞状态的候选值向量;Ct-1、Ct为上一时刻和当前时刻LSTM的元胞状态;ht-1、ht为上一时刻和当前时刻LSTM的隐藏状态;σ为sigmoid激活函数;*表示逐点相乘。
统计研究表明,95%的风电爬坡事件持续时间低于4.04h。因此,本发明的预测时间窗为未来4h,以采样分辨率15min算,共16步预测。假设当前时刻为t,LSTM的输入为{t,t-z},输出为{t,t+16}。通过当前时刻之前z个样本一次多步预测未来16个时刻的风电功率输出。一旦出现风电爬坡事件,这代表着未来4小时内、16个点的预测值曲线将会即时表现出来。
根据国家电网公司规定的技术指标,风功率超短期预测(未来4h)的月平均绝对误差要求小于20%,月均方根误差要求小于15%。采用标准化均方根误差(NRMSE)和标准化平均绝对误差(NMAE)对LSTM神经网络预测结果进行评估。计算公式如下:
步骤2:建立应对风电爬坡事件的提前-滞后调节方法,即利用储能系统进行提前-滞后两步补偿,降低风功率爬坡时的速率,同时降低对储能功率和容量的需求。其中,应对风电爬坡事件的提前-滞后调节方法的制定过程如下:
风电具有随机性和波动性的特征,风电爬坡事件是短时间内风功率发生大幅变化的事件,是风电波动性的体现。风电爬坡事件分为上爬坡和下爬坡,风电上爬坡时,当机组的调节速度无法跟随风电爬坡速度时,必须采取弃风措施确保电力系统运行安全;同理,风电下爬坡时,需采取必要的切负荷措施。
本发明将风功率表示为Pw,Δt为风功率采样分辨率,也是调度单时段持续时间。设风电爬坡事件的开始时间为t0,持续n个时段,持续时间为ΔT=n×Δt,幅值变化为ΔP。以上爬坡为例,风电爬坡事件的简图如图2所示。
图2中,斜线段为风电正在上爬坡,在t0时刻发生上爬坡事件,持续n=4个时段,持续时间为ΔT=4×Δt。
本发明将时段(t,t+Δt)内首末时刻风功率之差定义为t时段风电爬坡时的爬坡率ΔPw(t),如式(9)所示:
ΔPw(t)=Pw(t+Δt)-Pw(t) (9)
若t时刻的爬坡率ΔPw(t)大于给定阈值Pthreshold,则认为风电在t时段处于爬坡状态,即:
|ΔPw(t)|>Pthreshold (10)
常规调度计划制定时,由于对风功率的多步预测误差较大,难以预测爬坡事件,因而难以对爬坡事件做出准确对策,可能致使某些机组在风电爬坡时刻处于满载或最小技术出力,即无法提供上调或下调的爬坡备用,致使机组爬坡率下降,应对风电爬坡事件的能力降低。而通过LSTM神经网络,可以对风电功率的爬坡事件提前准确预测,这为提前协调各机组应对风电爬坡事件提供了可能。
将通过提前协调各机组使其能够在爬坡时段提供的最大爬坡率称为ΔPg。当预测到风电爬坡事件时,首先要判断风电在爬坡时段各时刻的爬坡率ΔPw是否大于机组系统的最大爬坡率ΔPg。若判定条件ΔPw<ΔPg,则通过合理分配爬坡事件发生之前时段各机组的出力,使其在风电发生爬坡的时刻有充足的爬坡备用,并使所有时段机组发电的总成本最小。若判定条件ΔPw>ΔPg,此时已无法通过机组的协同调节应对风电爬坡,为保证供电安全,避免切负荷现象,常规的做法是采取弃风手段,降低风电爬坡率,如图3所示。
与图2中的爬坡事件相同,图3中斜率较大的线段表示风电正在上爬坡,斜率较小的线段为机组在各时刻可提供的最大下爬坡备用,并已达到其最大爬坡率,则横线阴影部分为必须采取弃风的功率。弃风时段为t0到t0+ΔT+mΔt,弃风调节滞后风电爬坡结束时间mΔt。图3中是m=2的情况,当风电的爬坡率不同时,需要的调节时间也不同。同理,应对下爬坡事件时需提前mΔt采取弃风,且此时的提前预测显得更加重要。若无法设置准确的弃风量则仍会发生切负荷现象,而且这种传统的调节方法造成大量的风力资源浪费。本发明称这种采用弃风的调节方法为“弃风调节法”。
为减少或避免弃风现象,加入储能系统进行协调调度,在风电上爬坡需要弃风时进行充电,图3中需要弃风的部分即为需要充电的功率。同理,在风电下爬坡时提供放电。本发明称这种采用储能系统补偿弃风功率的方法为“补偿弃风法”。这样,弃风现象可以得到很大改善。但是,单位容量的储能系统配置成本较高,当风电爬坡持续时间较长时需要配置较大容量的储能系统以确保安全。因此,本发明采取“提前-滞后”的调节方法,可降低对储能容量的需要,如图4所示。
由图4所示,当发生风电上爬坡事件,在时段(t0-Δt,t0+2Δt)内,提前风电爬坡开始时间Δt使机组下爬坡,使其以最大爬坡率调节,此时段内由储能系统放电补足功率差额;在时段(t0+2Δt,t0+5Δt)内,滞后风电爬坡结束时间Δt使机组停止爬坡,此时段内由储能系统充电补足功率差额。风电发生下爬坡事件的调节方法与此类似,调节方向相反。与图3相比,此调节方法对储能容量的需求最大可降低50%(图4恰好为50%),且在发生风电下爬坡时,需要提前调节的时间由图3中的2Δt减少至Δt,即对预测的精度需求有所降低。图2至图4结合风电上爬坡的一种情况对本发明的调节方法作出说明,在具体调节中,提前和滞后的调节时间需根据风电的爬坡率优化调整。
步骤3:基于步骤1的预测结果和步骤2的调节方法,使用考虑风电爬坡事件的火储协调调度模型,提前协调机组,使机组在爬坡时有足够的可调备用容量,通过对不同的风电爬坡条件的判断协调机组和储能,实现更经济和安全的爬坡调度方案。其中,考虑风电爬坡事件的火储协调调度模型建立方法包括:
考虑机组煤耗成本和储能使用成本,得到调度时段总成本f最小为目标:
minf=CG+CB+CWC+CLS (11)
其中
式中:CG、CB、CWC、CLS为机组煤耗成本、储能运行成本、弃风成本和切负荷成本;T为调度总时间;N为火电机组总数;ai、bi、ci为机组i煤耗成本系数;Pi,t为机组i在时段t的出力;Pt B,ch、Pt B,dc为储能在时段t的充、放电功率;kB为储能单位容量的使用成本系数;Pt WC、kWC为弃风功率和单位弃风功率惩罚系数;Pt LS、kLS为切负荷功率和单位切负荷功率惩罚系数;
约束条件包括:
1)机组和储能出力约束:
式中:Pi为机组i瞬时出力,Pi max、Pi min为机组i最大、最小技术出力;PB,max为储能最大充放电功率;
2)功率平衡约束:
式中:Pt f、Pt L为风电功率和负荷功率预测值,用风功率预测值Pt f代替步骤2中的风电功率Pw作为调节策略制定的依据;
3)机组爬坡约束:
式中:ri u、ri d为机组i单位时间的爬坡率;Δt为调度单时段持续时间;Pi,t-1为机组i在t-1时段的出力;
4)旋转备用约束:
5)储能容量约束:
6)爬坡对策判断条件:
式中:ΔPt w,u、ΔPt w,d为风电在t时刻的上、下爬坡率;ΔPt L,u、ΔPt L,d为负荷在t时刻的上、下爬坡率;ΔPt g,u、ΔPt g,d为所有机组在时段t可提供的上调和下调备用,二者统称为步骤2判定条件ΔPw<ΔPg中的ΔPg,并且考虑负荷的爬坡,将判断条件具体为式(18);
当所有机组均可调时,称ΔPt g为机组在时段t的最大爬坡率。若机组最大爬坡率大于机组和负荷爬坡率之和时,即满足此约束条件时,通过提前协调机组应对风电爬坡,此时无需使用储能系统,模型中Pt B,ch=Pt B,dc=0;反之,若不满足此约束条件,则需要进行提前弃风或使用储能系统参与协调调度,按步骤2提前-滞后调节法和的所述火储协调调度模型进行优化求解。
本实施例采用IEEE39节点系统进行仿真分析,共有10台火电机组。风电场数据为ELIA比利时电力运营商公开的运行数据。储能系统总容量为100MW,其它参数见表1。
表1储能系统参数设置
算例1:基于LSTM神经网络的风功率预测效果评估
为验证LSTM神经网络的预测效果,采用ELIA比利时电力运营商公开的2014年1月~5月的数据分别对BP神经网络、GR神经网络、Elman神经网络和LSTM神经网络进行训练,并分别对2014年6月共30天的数据进行预测。式(1)、(2)中样本数s=2880。表2为不同预测方法预测结果的评估指标比较,预测提前的时间从提前15min到提前4h。其中底部方框中表示超出国家电网公司规定指标,右侧方框表示四种预测方法的最优值。
表2不同预测方法和不同提前预测时间的NRMSE(%)和NMAE(%)比较
由表2可见,LSTM神经网络在提前不同时间情况下所得NRMSE和NMAE指标与其它三种神经网络相比均是最低的。如提前4h预测所得NMAE,BP神经网络为20.09%,Elman神经网络为19.68%,GR神经网络为17.28%,均超出了国家电网公司公司规定的指标(15%),而LSTM神经网络为3.47%,与其它三种神经网络相比分别降低了16.62%、16.21%和13.81%,远低于规定指标。又如提前1h预测所得NRMSE,BP神经网络为10.58%,Elman神经网络为10.76%,GR神经网络为6.05%,而LSTM神经网络为5.49%,与其它三种神经网络相比分别降低了5.09%、5.27%和0.56%,远低于规定的指标(20%)。说明LSTM神经网络的预测精度在四种预测方法里是最高的。
比较提前不同时间的预测结果,可见,LSTM神经网络在预测提前时间增大的过程中两项性能指标变化不大,依旧能够实现高精度预测,而其它三种神经网络性能指标变化明显,预测精度迅速衰减。在预测提前时间从15min变化到4h的过程中,LSTM神经网络的NMAE从2.48%变化至3.47%,仅下降了0.99%的精度,而在此过程中,BP神经网络的NMAE从3.30%变化至20.09%,精度下降了16.79%,同样,Elman神经网络和GR神经网络精度分别下降了16.49%和14.25%。由此可见,LSTM神经网络能够实现较长期的高精度预测,得益于其独特的记忆和遗忘设计,能够充分关联较长期的历史信息。这对提前制定和调节火储系统发电计划以应对风电爬坡事件具有非常重要的指导意义。
算例2:风电上爬坡事件调度算例
将风电爬坡事件的阈值Pthreshold设定为15min内变化量为装机容量的10%,在2014年6月1日至2015年5月31日一年内的风功率数据中选取风电爬坡事件样本进行预测和调度,并按照装机容量500MW进行归一化折算。取LSTM神经网络得到的多步预测结果,即预测未来4h内的风功率。调度总时段T=16,单时段持续时间Δt=15min,调度周期为4h。为促进风电消纳,尽量避免切负荷,设置单位功率弃风成本系数kWC=100$/MW,单位功率切负荷成本系数kLS=1000$/MW。
取2014年10月4日15:30~19:30发生的风电上爬坡事件。图5左纵轴为分别采用BP神经网络和LSTM神经网络预测所得的风功率预测值曲线和实际值曲线,右纵轴为负荷曲线。
由图5预测结果可见,BP神经网络所得预测值在t=6时刻开始误差明显增大,说明其对多步预测,特别是发生风电爬坡事件时的预测误差较大,其16个时段的NMAE为23.41%,超出了规定指标。与此相比,LSTM网络所得结果在大多数时段误差均小于BP神经网络,更接近实际值,其16个时段的NMAE为3.68%,比前者降低了19.73%,能实现风电爬坡事件的多步准确预测。
为说明本发明所提调节方法和调度模型的有效性,本例暂不考虑其它预测值,不考虑预测误差,按图5中LSTM网络得到的预测值和负荷曲线进行分析和调度。可见,风功率在t=5时刻发生上爬坡事件,而负荷发生下爬坡,出现了极端情景。由公式(16)和(18),通过提前调节能使所有机组在15min内可提供的最大备用为127.5MW,而在t=[8,9]的15min内风功率和负荷的爬坡功率之和达到163.8MW,超出机组可用备用36.3MW。此时,需采取步骤2中所述的三种调节方式。采取弃风调节法和补偿弃风法的调节过程如图6所示,采用本发明提出的提前-滞后调节方法的调节过程如图7所示(只画出有调节的附近几个时段)。
图6所示为总弃风功率最小的调节策略。可见,弃风调节法需在t=9时刻和t=10时刻分别弃风46.6MW和7.6MW才能满足机组的爬坡能力范围,共造成54.2MW的弃风。若采用补偿弃风法,利用储能平抑弃风,则需单时段充电功率最少为46.6MW,在本例设置储能约束下,储能在t=9时刻最多可充电25MW,造成弃风功率21.6MW。若要避免弃风,需配置更大输出功率和容量的储能系统才能满足需要,储能配置成本会显著增加。且若出现连续的风电上爬坡,储能系统将一直处于充电状态,易造成过度充电,影响储能系统寿命。相比之下,采用如图7所示的提前-滞后调节法,提前1个时刻进行调节,使储能系统在t=8时刻放电25.0MW,在t=9时刻充电21.6MW,降低风功率的爬坡速率;同时,使机组在其爬坡能力范围内提前开始提供下爬坡备用,降低出力,接纳风电。利用火储协调调度,使整个调度过程的弃风总功率降至0。在调节过程中,储能的单时段放电功率最大为25.0MW,所设容量为100MW储能系统满足要求。同时,所需的储能总输出功率由54.2MW降低至46.6MW,降低了14%。可见,提前-滞后调节法不仅避免了弃风,且对储能输出功率和配置容量的需求降低,从而节省了储能配置和运行成本。
在不计预测误差的前提下,三种调节方法的调度成本比较如表3所示。
表3各调节方法的调度成本比较(单位:$)
由表3可见,提前-滞后法所需总成本在三种方法中是最低的,在本例的风电爬坡事件下,能够在保证不切负荷、不弃风的情况下实现最经济的调度计划。
算例3:风电下爬坡事件调度算例
为说明LSTM神经网络多步预测方法在爬坡调度中的指导意义,本例考虑预测误差,取2015年5月18日17:30~21:30发生的风电下爬坡事件。图8左纵轴为分别采用GR神经网络和LSTM神经网络预测所得的风功率预测值曲线和实际值曲线,右纵轴为负荷曲线,对两种方法得到的预测值分别进行分析和调度。
由图8可见,GR神经网络在前2个时刻预测较准确,而在第3时刻至第16时刻误差逐渐增大,说明GR网络在进行多步预测时预测精度明显下降,其16个时段的NRMSE为19.03%。与此相比,LSTM神经网络所得结果在大部分时段预测误差均小于GR神经网络,更接近实际值曲线,其16个时段的NRMSE为5.31%,比前者降低了13.72%,能够实现更多时间步的精确预测。
由图8,风功率在t=6时刻发生急剧下爬坡事件,此时负荷处于上爬坡状态,出现了极端情景。利用火储系统、采用提前-滞后的调节策略分别对两种预测方法得到的风功率爬坡预测值进行爬坡调节,并将优化前风功率的实际值与优化后风功率的调度值和进行比较。
由图9所示,由于GR神经网络对爬坡时段的风功率预测误差较大,若按此预测值制定调度计划,会致使大部分时段的调度计划很大程度的偏离实际值,从而使火储系统提前调节的功率不足。若得不到及时的修正,会在4h的调度周期内共造成554.7MW的切负荷和20.4MW的弃风,给电力系统的安全供电带来严重影响。若GR神经网络的预测值与实际值相比不是偏大而是偏小同等的误差,则会使火储系统调节过度,还会造成554.7MW的弃风,造成能源的浪费和电力系统运行经济性的严重下降。与此相比,由图10可见,由于LSTM神经网络能够保证较长时段预测的准确性,从而为应对风电爬坡调节策略的提前制定提供很大的依据。在4h调度周期内的总切负荷量为35.6MW,比前者降低了519.1MW,降低了93.6%。根据两种预测值进行调度的成本比较如表4所示。
表4各调节方法的调度成本比较(单位:$)
由表4可见,根据LSTM神经网络进行调度时,由于避免了大规模的切负荷风险,总成本比GRNN神经网络调度降低了51.2%。说明根据LSTM的多步精确预测,使火储系统能够提前制定更准确的调节策略,应对风电爬坡事件。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (2)
1.基于LSTM的含储能电力系统风电爬坡协调调度方法,其特征在于,包括以下步骤:
步骤1:利用LSTM神经网络对风功率爬坡事件进行提前多步预测;
步骤2:利用储能系统进行提前-滞后两步补偿,降低风功率爬坡时的速率,同时降低对储能功率和容量的需求;
步骤3:基于步骤1的预测结果和步骤2的调节方法,使用考虑风电爬坡事件的火储协调调度模型,提前协调机组;
所述步骤2中包括如下具体步骤:设风电爬坡事件的开始时间为t0,持续n=4个时段,持续时间为ΔT=4×Δt,Δt为每个时段的时长;当发生风电上爬坡事件,在时段(t0-Δt,t0+2Δt)内,提前风电爬坡开始时间Δt使机组下爬坡,使机组以最大爬坡率调节,所述时段(t0-Δt,t0+2Δt)内由储能系统放电补足功率差额;在时段(t0+2Δt,t0+5Δt)内,滞后风电爬坡结束时间Δt使机组停止爬坡,所述时段(t0+2Δt,t0+5Δt)内由储能系统充电补足功率差额;当风电发生下爬坡事件时的调节方法与发生风电上爬坡事件时的调节方向相反。
2.根据权利要求1所述的基于LSTM的含储能电力系统风电爬坡协调调度方法,其特征在于:所述步骤3中考虑风电爬坡事件的火储协调调度模型建立方法包括如下步骤:
考虑机组煤耗成本和储能使用成本,得到调度时段总成本f最小为目标:
minf=CG+CB+CWC+CLS
其中
式中:CG、CB、CWC、CLS为机组煤耗成本、储能运行成本、弃风成本和切负荷成本;T为调度总时间;N为火电机组总数;ai、bi、ci为机组i煤耗成本系数;Pi,t为机组i在时段t的出力;Pt B ,ch、Pt B,dc为储能在时段t的充、放电功率;kB为储能单位容量的使用成本系数;Pt WC、kWC为弃风功率和单位弃风功率惩罚系数;Pt LS、kLS为切负荷功率和单位切负荷功率惩罚系数;
约束条件包括:
1)机组和储能出力约束:
式中:Pi为机组i瞬时出力,Pi max、Pi min为机组i最大、最小技术出力;PB,max为储能最大充放电功率;
2)功率平衡约束:
式中:Pt f、Pt L为风电功率和负荷功率预测值;
3)机组爬坡约束:
式中:ri u、ri d为机组i单位时间的爬坡率;Δt为调度单时段持续时间;Pi,t-1为机组i在t-1时段的出力;
4)旋转备用约束:
5)储能容量约束:
6)爬坡对策判断条件:
式中:ΔPt g,u、ΔPt g,d为所有机组在时段t可提供的上调和下调备用;ΔPt w,u、ΔPt w,d为风电在t时刻的上、下爬坡率;ΔPt L,u、ΔPt L,d为负荷在t时刻的上、下爬坡率;
若机组最大爬坡率大于风电和负荷爬坡率之和时,通过提前协调机组应对风电爬坡,此时无需使用储能系统,模型中Pt B,ch=Pt B,dc=0;反之,则需要进行提前弃风或使用储能系统参与协调调度,按步骤2提前-滞后调节法和所述火储协调调度模型进行优化求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418323.5A CN110071530B (zh) | 2019-05-20 | 2019-05-20 | 基于lstm的含储能电力系统风电爬坡协调调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418323.5A CN110071530B (zh) | 2019-05-20 | 2019-05-20 | 基于lstm的含储能电力系统风电爬坡协调调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110071530A CN110071530A (zh) | 2019-07-30 |
CN110071530B true CN110071530B (zh) | 2020-10-02 |
Family
ID=67371072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910418323.5A Expired - Fee Related CN110071530B (zh) | 2019-05-20 | 2019-05-20 | 基于lstm的含储能电力系统风电爬坡协调调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110071530B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110707745B (zh) * | 2019-10-16 | 2023-05-16 | 中国矿业大学 | 基于改进vmd的电热综合系统多时间尺度经济调度方法 |
CN110751346B (zh) * | 2019-11-04 | 2023-06-13 | 重庆中涪科瑞工业技术研究院有限公司 | 一种基于驾驶速度预测和博弈论的分布式能源管理方法 |
CN111429020B (zh) * | 2020-04-07 | 2023-12-29 | 中国矿业大学 | 考虑热网储热特性的电热系统多时间尺度经济调度方法 |
CN111582557A (zh) * | 2020-04-20 | 2020-08-25 | 哈尔滨工业大学 | 基于变差函数的风电爬坡事件多级预警方法 |
CN112036658B (zh) * | 2020-09-09 | 2022-06-14 | 华东交通大学 | 一种基于随差遗忘长短期记忆的风电功率实时预测方法 |
CN112134304B (zh) * | 2020-09-22 | 2022-04-01 | 南方电网数字电网研究院有限公司 | 一种基于深度学习的微电网全自动导航方法、系统与装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107909227B (zh) * | 2017-12-20 | 2022-07-15 | 北京金风慧能技术有限公司 | 超短期预测风电场功率的方法、装置及风力发电机组 |
CN108615097A (zh) * | 2018-05-10 | 2018-10-02 | 广东工业大学 | 一种风速预测方法、系统、设备及计算机可读存储介质 |
-
2019
- 2019-05-20 CN CN201910418323.5A patent/CN110071530B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110071530A (zh) | 2019-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110071530B (zh) | 基于lstm的含储能电力系统风电爬坡协调调度方法 | |
CN109687530B (zh) | 一种考虑阻塞和储能分时电价的电网混合滚动调度方法 | |
CN105846461B (zh) | 一种大规模储能电站自适应动态规划的控制方法和系统 | |
CN107818385B (zh) | 一种梯级水电站群实时运行趋势预测方法 | |
CN109599861B (zh) | 考虑本地负荷调峰能力的送端电网电源结构规划方法 | |
CN109347151B (zh) | 一种新能源参与调峰的送端电网电源结构优化方法 | |
CN110581571A (zh) | 一种主动配电网动态优化调度方法 | |
WO2017161785A1 (zh) | 基于储能运行状态的光伏功率稳定输出控制方法 | |
CN109390973B (zh) | 一种考虑通道约束的送端电网电源结构优化方法 | |
CN109765787B (zh) | 一种基于日内-实时滚动控制的配电网源荷快速跟踪方法 | |
WO2023065113A1 (zh) | 风光水多能互补系统灵活性需求量化及协调优化方法 | |
CN105162149A (zh) | 基于模糊自适应调节的光储系统跟踪发电计划出力方法 | |
CN112186811B (zh) | 一种基于深度强化学习的agc机组动态优化方法 | |
CN110348599B (zh) | 一种考虑弃水风险的跨流域水电站群调峰优化调度方法 | |
CN104181895A (zh) | 适应新能源接入的短期与超短期协调滚动调度优化策略 | |
CN111476407B (zh) | 联合风电光伏电站的梯级水电站中长期隐随机调度方法 | |
CN108054771A (zh) | 一种储能系统充放电控制方法及系统 | |
CN110137952B (zh) | 一种实现源-荷-储协调运行的优化模型与方法 | |
Li et al. | Direct control strategy of real-time tracking power generation plan for wind power and battery energy storage combined system | |
CN113408962A (zh) | 一种电网多时间尺度多目标能量优化调度方法 | |
CN104795830A (zh) | 一种利用多类型储能系统跟踪发电计划出力的控制方法 | |
CN114819373A (zh) | 一种基于合作博弈的共享混合储能电站的储能规划方法 | |
CN114243794A (zh) | 风光水多能互补系统灵活性需求量化及协调优化方法 | |
CN115566741A (zh) | 基于lstm预测校正的储能多时间尺度优化控制方法 | |
CN108009672A (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 | ||
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: 20201002 Termination date: 20210520 |