CN107846045B - 一种电池储能系统跟踪风光计划出力的控制方法 - Google Patents

一种电池储能系统跟踪风光计划出力的控制方法 Download PDF

Info

Publication number
CN107846045B
CN107846045B CN201711211856.3A CN201711211856A CN107846045B CN 107846045 B CN107846045 B CN 107846045B CN 201711211856 A CN201711211856 A CN 201711211856A CN 107846045 B CN107846045 B CN 107846045B
Authority
CN
China
Prior art keywords
energy storage
storage system
battery energy
output
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.)
Active
Application number
CN201711211856.3A
Other languages
English (en)
Other versions
CN107846045A (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 GASU ELECTRIC POWER RESEARCH INSTITUTE
State Grid Corp of China SGCC
State Grid Gansu Electric Power Co Ltd
Lanzhou Jiaotong University
Original Assignee
STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
State Grid Corp of China SGCC
State Grid Gansu Electric Power Co Ltd
Lanzhou Jiaotong 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 STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE, State Grid Corp of China SGCC, State Grid Gansu Electric Power Co Ltd, Lanzhou Jiaotong University filed Critical STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
Priority to CN201711211856.3A priority Critical patent/CN107846045B/zh
Publication of CN107846045A publication Critical patent/CN107846045A/zh
Application granted granted Critical
Publication of CN107846045B publication Critical patent/CN107846045B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/383
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • 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
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/10Flexible AC transmission systems [FACTS]
    • 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

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Secondary Cells (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

本发明对风电场、光伏电站所具有的随机性和间歇性的特点,电场实际发电功率不能跟踪计划发电功率,从而导致对电网安全性、稳定性以及电能质量等造成巨大的冲击和影响,以及出现弃风弃光现象的问题。发明提出了一种基于模型预测控制的电池储能系统跟踪风光计划出力的控制方法,提高风光储联合发电系统中储能系统跟踪计划出力能力,平滑风光联合发电输出、削峰填谷,提高风光发电利用率。该方法根据采集风电场、光伏电站日前功率计划出力数据、建立储能系统功率出力模型,制定出参考模型轨迹,建立储能系统出力闭环预测模型以及控制算法的滚动优化,求出满足跟踪期望轨迹优化目标函数最小的控制量,实现对储能系统跟踪风光计划出力精确的控制。

Description

一种电池储能系统跟踪风光计划出力的控制方法
技术领域
本发明属于新能源发电控制技术领域,涉及一种基于模型预测控制的电池储能系统跟踪风光计划出力的控制方法。
背景技术
近年来,由于在风光储联合发电中应用储能系统跟踪计划出力以提高风光发电效率的应用价值极高,不但可以平滑风光联合发电输出、削峰填谷,而且可以减少弃风弃光的现象,提高风光发电利用率,因此风光储联合发电中应用储能系统跟踪计划出力以提高风光发电效率的控制问题近几年来成为一个热点问题。
目前,国内外均已有了多项关于风光储联合应用的示范项目工程,但对储能系统的控制与分析更多的是关于平滑风电输出、削峰填谷以及系统调频几方面,基于风光储联合应用以提高跟踪风光发电计划出力能力的研究成果还相对较少。由于蓄电池储能系统有着充放电功率及储能容量的约束条件,并且过大的充、放电深度均会对储能电池的寿命造成影响,特别是在储能系统运行过程中,需要控制其荷电状态保持在适宜的范围内,荷电状态过高或者过低都会减少储能系统的剩余可用容量,降低储能系统的充放电功率最大允许值,从而降低储能系统的控制效果;同时又由于风光发电日前预测误差经常会出现在一个时间段内,出现持续不满足预测误差要求的情况,对储能系统在某一时段的持续出力有了更严格的要求,故控制好电池储能系统的荷电状态,实时优化控制储能系统的出力对于改善风电跟踪计划出力具有重要的意义。
发明内容
针对风电场、光伏电站所具有的随机性和间歇性的特点,电场实际发电功率不能跟踪计划发电功率,从而导致对电网安全性、稳定性以及电能质量等造成巨大的冲击和影响,以及出现弃风弃光现象的问题。为实现对储能系统出力精确的控制,提出了一种基于模型预测控制的电池储能系统跟踪风光计划出力的控制方法,进而提高风光储联合发电系统中储能系统跟踪计划出力能力,平滑风光联合发电输出、削峰填谷,也提高了风光发电利用率。该方法根据采集风电场、光伏电站日前功率计划出力数据、建立储能系统功率出力模型,制定出参考模型轨迹,通过建立储能系统出力闭环预测模型以及控制算法的滚动优化,求出满足跟踪期望轨迹优化目标函数最小的控制量,从而实现对储能系统跟踪风光计划出力精确的控制。
具体步骤为:
步骤一,采集风电场、光伏电站的功率计划出力数据;
步骤二,基于第一步采集风电场、光伏电站的功率计划出力数据,建立电池储能系统的功率出力模型;
步骤三,基于第二步建立电池储能系统的功率出力模型,建立参考轨迹模型;
步骤四,基于第三步建立参考轨迹模型,建立闭环预测模型;
步骤五,基于第四步建立闭环预测模型,根据模型预测控制算法计算得出最优控制量。
进一步的,步骤一中采集数据,求出系统设定值
ys(k)=yf(k)-yw(k) (3)
式中,yf(k)为k时刻系统风电场、光伏电站的计划功率出力值;yw(k)为k时刻系统风电场、光伏电站的实际功率出力值;ys(k)为系统k时刻系统的设定值,它是由电池系统的计划功率出力值与实际功率出力值做差所得。
进一步的,步骤二中电池储能系统的功率出力模型采用对电池储能系统的电量进行逐时刻积累的方式,进而来表示电池储能系统的充、放电的过程,用电池储能系统的充、放电过程中的效率近似的等效充、放电过程中的电能损耗;电池储能系统不断变化的剩余电量与任意时刻的充放电功率、充放电效率和自放电率相关,剩余电量递推关系如下:
充电过程
E(k)=(1-σsdr)E(k-1)-y(k)ΔtηC (4)
放电过程
E(k)=(1-σsdr)E(k-1)-y(k)Δt/ηD (5)
式中,E(k)为k时刻末该电池储能系统的剩余电池容量,MW·min;y(k)为k时刻电池储能系统的充、放电功率值;σsdr为电池储能系统的自放电率,min;ηC为该电池储能系统的充电效率;ηD为该电池储能系统的放电效率;Δt为计算窗口时长,min;
在电池储能系统充、放电过程中,其允许的最大充电功率和最大放电功率以及剩余容量的边界分别如下:
充电过程、放电过程对于电池储能系统最大出力的限制分别为:
Figure BDA0001484745290000031
Figure BDA0001484745290000032
电池储能系统剩余容量的限制为:
Emin≤E(k)≤Emax (8)
电池储能系统k时刻充放电功率限制由该时刻电池储能系统的剩余电量和电池储能系统自身特性共同决定,具体表达式如下:
电池储能系统k时刻,充电过程对于充电功率限制值为:
Figure BDA0001484745290000033
电池储能系统k时刻,放电过程对于放电功率限制值为:
Figure BDA0001484745290000034
式中,
Figure BDA0001484745290000035
为电池储能系统最大允许充电功率值(MW),由电池储能系统自身特性决定;/>
Figure BDA0001484745290000036
为电池储能系统最大允许放电功率值(MW),由电池储能系统自身特性决定;ymin为电池储能系统最小容量限制(MW·h);ymax为电池储能最大;容量限制(MW·h);ylimit(k)为电池储能介质充、放电功率限制值。
进一步的,步骤三中在线校正建立电池储能系统的闭环预测模型方法为:
将第k步的实际对象的输出测量值与预测模型输出之间的误差附加到模型的预测输出ym(k+i)上,得到闭环预测模型,用yp(k+1)表示,即:
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (11)
式中:
yp(k+1)=[yp(k+1),yp(k+2),…,yp(k+P)]T,h0=[1,1,…,1]T (12)
如果以i时刻为预测模型的起点,可以列出其闭环预测模型为:
Figure BDA0001484745290000041
已知脉冲响应和阶跃响应系数之间有如下关系:
Figure BDA0001484745290000042
Figure BDA0001484745290000043
展开式稍加整理,然后将后两式代入,可以得到闭环预测模型为:
yp(k+1)=h0y(k)+p+AΔu(k+1) (16)
式中,ys为采集到的风电场、光伏电站出力计划值作为系统的设定值;P为预测步长;ym(k)为电池储能系统的预测模型输出;
p=[P1,P2,…,PP]T
h0=[1,1,…,1]T
yp(k+1)=[yp(k+1),yp(k+2),…,yp(k+P)]T
Δu(k+1)=[Δu(k),Δu(k+1),…,Δu(k+P-1)]T
Figure BDA0001484745290000044
进一步的,步骤四中参考轨迹模型预测控制的目标是使电池储能系统的输出y(t)沿着一条预先规定的轨迹逐渐达到设定值ysp;这条预先规定的轨迹称为参考轨迹yr;参考轨迹采用从现在时刻实际值出发的一阶指数形式;它在未来P个时刻的值为:
Figure BDA0001484745290000051
式中,yr(k+i)为参考模型轨迹函数;T为采用周期;α为柔化系数,且0<α=exp(-T/τ)<1;τ为参考轨迹的时间常数。
进一步的,步骤五中滚动优化步骤为:
首先,假设被控对象实际脉冲相应及其预测模型的脉冲响应分别为:
h=[h1,h2,…,hN]T (18)
Figure BDA0001484745290000052
可以得到闭环预测模型为:
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (20)
可以得到闭环预测控制模型的单步控制输出为:
Figure BDA0001484745290000053
已知电池储能系统的预测模型和闭环校正预测模型分别为:
Figure BDA0001484745290000054
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (23)
并且已知电池储能系统输出的参考轨迹为yr(k+1),设系统误差方程为:
e(k+1)=yr(k+1)-yp(k+1) (24)
优化性能函数为:
Figure BDA0001484745290000055
式中,Q1为非负定加权对称矩阵;Q2为正定控制加权对称矩阵;
Figure BDA0001484745290000061
使得上式优化性能目标函数最小,可以求得最优控制量Δu2为:
Figure BDA0001484745290000062
式中,e'为参考轨迹与零输入响应下闭环预测输入之差,即为:
Figure BDA0001484745290000063
本发明的优点在于,采用模型预测控制的电池储能系统跟踪风光计划出力的控制方法时,提高风光储联合发电系统中储能系统跟踪计划出力能力,并且平滑风光联合发电输出、削峰填谷,提高了风光发电利用率,安全性、稳定性和电能质量等都很大的得到了提高;同时该方法还具有相应速度快、风光联合发电的预测性高,出力偏差小的优点。
附图说明
图1预测控制原理示意图
图2预测控制算法原理图
图3预测控制系统结构框图
图4滚动优化原理图
图5在线误差校正原理图
图6电池储能系统模型预测控制设计框图
图7基于模型预测控制的电池储能系统跟踪风光计划出力的流程图
具体实施方式
1.理论基础
预测控制原理示意图如图1所示。
如图1中ysp为设定值,即系统要弥补的差值;y(t)为控制对象的输出值;0时刻代表当前的时刻,0时刻左右两边的曲线分别代表过去、将来的输出与控制;预测控制的算法就是要根据预测输出yp(k)与期望输出yr(k)的差值e(k),计算出当前和将来L时刻使得e(k)最小的最优控制量u(k)。
预测控制是是通过达到某一性能指标最优来确定未来的控制作用。传统意义的最优控制一般是指全局最优,而预测控制中的优化是一种有限时间段的滚动优化。即是说,在预测控制中,优化不是一次离线进行,而是不同时间段内反复在线进行的。
预测控制算法原理图如图2所示。
如图2中,u(k)为对被控对象的控制作用;y(k+1)为被控对象k+1时刻的输出值;y(k)为被控对象k时刻的输出值;ym(k+1)为被控对象预测模型k+1时刻的预测输出值;ym(k)为被控对象预测模型k时刻的预测输出值。根据算法实现的具体形式,可以得到不同的预测控制衍生算法。
预测控制系统结构框图如图3所示。
如图3中,预测控制系统由:参考轨迹,预测模型,滚动优化,在线校正等四个部分组成,如图3所示。预测控制是以不同的预测模型为基础,采用在线滚动优化指标和反馈自校正策略,克服被控对象的不确定性、不稳定性等因素的影响,达到预期的控制目标——参考轨迹输入,使系统具有很好的鲁棒性和稳定性。
(1)参考轨迹
在预测控制过程中,考虑到过程的动态特性,为了避免控制过程出现输入和输出的剧烈变化,一般情况下会要求控制过程的轨迹由当前时刻输出值y(k)沿着一条期望的平缓曲线达到设定值,这条由设定值经过在线柔化的曲线被称为参考轨迹,一般表示为yr(k),最广泛采用的参考轨迹为一阶指数变化形式:
yr(k+i)=αiy(k)+(1-αi)yr (1)
其中
Figure BDA0001484745290000071
式中,Ts为采用周期;T为参考轨迹的时间常数;下标r指示参考值(ReferenceValue);yr为设定值。显然T越小,α越小,参考轨迹也就能够很快地达到设定值yr
(2)预测模型
预测模型的作用即是根据被控对象的历史信息和未来的输入,预测被控对象未来的输出,具有展示系统未来动态行为的能力。从方法的角度讲,只要是具有预测功能的信息集合,均可作为预测模型。
(3)滚动优化
滚动优化原理图如图4所示。
如图4所示,1—为参考轨迹;2—为最优预测输出;3—为最优控制作用。
在预测控制中的在线优化方法是通过某一性能指标的最优来确定未来的控制作用,性能指标中涉及系统未来的行为,是根据预测模型由未来的控制动作决定的,在MPC中的优化不是采用一个不变的全局优化指标,而是采用滚动式的、通常是有限时域的优化策略。如图4所示,在每一采样时刻,优化性能指标通常只涉及未来的有限时间,而到下一采样时刻,这一优化时域向前推移。因此MPC在每一时刻有一个相对于该时刻的优化性能指标。不同的时刻优化性能指标的相对形式是相同的,但其绝对形式则是不同的。通常优化不是一次离线进行、而是反复在线运行,这就是滚动优化的含义所在。
(4)在线校正
在线误差校正原理图如图5所示。
如图5中,uk-1为k-1时刻的控制作用;uk为k时刻的控制作用;
Figure BDA0001484745290000081
为k时刻模型预测的k+1时刻的输出;/>
Figure BDA0001484745290000082
为k+1时刻的实际输出;y1为k时刻模型预测的输出;y2为k+1时刻模型预测的输出。在线校正的形式是多种多样的,可以在保持模型不变的基础上,对未来的误差做出预测并加以补偿,也可以采用在线辨识的方法直接修改预测模型。
2.电池储能系统跟踪风光计划出力的预测控制
电池储能系统模型预测控制设计框图如图6所示。
如图6中,ysp为系统设定值,即系统要弥补的差值,y(k)为电池储能系统的功率出力,yr(k+1)为对象输出的参考轨迹,e(k+1)为系统的误差,ym(k)为电池储能系统k时刻的预测模型输出,ym(k+1)为电池储能系统k+1时刻的预测模型输出,yp(k+1)为电池储能系统k+1时刻在线校正的预测模型输出,Δu(k)为系统计算得出的最优控制量。
模型预测控制的电池储能系统跟踪风光计划出力的流程图如图7所示。
如图7所示,第一步采集风电场、光伏电站的功率计划出力数据;第二步基于第一步采集风电场、光伏电站的功率计划出力数据,建立电池储能系统的功率出力模型;第三步基于第二步建立电池储能系统的功率出力模型,建立参考轨迹模型;第四步基于第三步建立参考轨迹模型,建立闭环预测模型;第五步基于第四步建立闭环预测模型,根据模型预测控制算法计算得出最优控制量。
2.1采集数据,求出系统设定值
ys(k)=yf(k)-yw(k) (3)
式中,yf(k)为k时刻系统风电场、光伏电站的计划功率出力值;yw(k)为k时刻系统风电场、光伏电站的实际功率出力值;ys(k)为系统k时刻系统的设定值,它是由电池系统的计划功率出力值与实际功率出力值做差所得。
2.2、电池储能系统的功率出力模型
本办法采用对电池储能系统的电量进行逐时刻积累的方式,进而来表示电池储能系统的充、放电的过程,用电池储能系统的充、放电过程中的效率近似的等效充、放电过程中的电能损耗。电池储能系统不断变化的剩余电量与任意时刻的充放电功率、充放电效率和自放电率相关,剩余电量递推关系如下:
充电过程
E(k)=(1-σsdr)E(k-1)-y(k)ΔtηC (4)
放电过程
E(k)=(1-σsdr)E(k-1)-y(k)Δt/ηD (5)
式中,E(k)为k时刻末该电池储能系统的剩余电池容量,MW·min;y(k)为k时刻电池储能系统的充、放电功率值;σsdr为电池储能系统的自放电率,min;ηC为该电池储能系统的充电效率;ηD为该电池储能系统的放电效率;Δt为计算窗口时长,min。
在电池储能系统充、放电过程中,其允许的最大充电功率和最大放电功率以及剩余容量的边界分别如下:
充电过程、放电过程对于电池储能系统最大出力的限制分别为:
Figure BDA0001484745290000091
Figure BDA0001484745290000092
电池储能系统剩余容量的限制为:
Emin≤E(k)≤Emax (8)
电池储能系统k时刻充放电功率限制由该时刻电池储能系统的剩余电量和电池储能系统自身特性共同决定,具体表达式如下:
电池储能系统k时刻,充电过程对于充电功率限制值为:
Figure BDA0001484745290000101
电池储能系统k时刻,放电过程对于放电功率限制值为:
Figure BDA0001484745290000102
式中,
Figure BDA0001484745290000103
为电池储能系统最大允许充电功率值(MW),由电池储能系统自身特性决定;/>
Figure BDA0001484745290000104
为电池储能系统最大允许放电功率值(MW),由电池储能系统自身特性决定;ymin为电池储能系统最小容量限制(MW·h);ymax为电池储能最大;容量限制(MW·h);ylimit(k)为电池储能介质充、放电功率限制值。
在电池储能系统的实际运行过程中,为了使得电池储能系统具有比较强的充电能力和放电能力,剩余电量最佳的状态是保持在额定容量的60%左右,剩余电量过高和过低都会降低储能系统允许的最大充、放电功率值,从而降低了电池储能系统对风电功率的控制效果。
2.3、在线校正建立电池储能系统的闭环预测模型
本预测控制中的优化不仅是基于模型,而且利用了反馈信息构成了闭环优化。通过优化方法确定了一系列未来的控制作用后,为防止模型误差或环境干扰引起控制效果偏离理想输出,并不是把所有控制作用逐一全部施加,是将下一时刻的控制作用,要通过此时刻的实际输出与模型预测输出之间的偏差对预测模型进行修正后,再通过新的优化产生新的控制作用,作用到电池储能系统上去。即第k步的实际对象的输出测量值与预测模型输出之间的误差附加到模型的预测输出ym(k+i)上,得到闭环预测模型,用yp(k+1)表示,即:
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (11)
式中:
yp(k+1)=[yp(k+1),yp(k+2),…,yp(k+P)]T,h0=[1,1,…,1]T (12)
如果以i时刻为预测模型的起点,可以列出其闭环预测模型为:
Figure BDA0001484745290000111
已知脉冲响应和阶跃响应系数之间有如下关系。
Figure BDA0001484745290000112
Figure BDA0001484745290000113
展开式稍加整理,然后将后两式代入,可以得到闭环预测模型为:
yp(k+1)=h0y(k)+p+AΔu(k+1) (16)
式中,ys为采集到的风电场、光伏电站出力计划值作为系统的设定值;P为预测步长;ym(k)为电池储能系统的预测模型输出;
p=[P1,P2,…,PP]T
h0=[1,1,…,1]T
yp(k+1)=[yp(k+1),yp(k+2),…,yp(k+P)]T
Δu(k+1)=[Δu(k),Δu(k+1),…,Δu(k+P-1)]T
Figure BDA0001484745290000114
从上式可以看出,由于每个预测时刻都引入了当时实际对象输出和模型输出的偏差,使闭环预测模型不断得到修正,显然这种方法可以有效地克服模型的不精确性和系统中存在的不确定性。因此,反馈修正是预测控制的特点之一。
2.4、参考轨迹模型
预测控制的目标是使电池储能系统的输出y(t)沿着一条预先规定的轨迹逐渐达到设定值ysp。这条预先规定的轨迹称为参考轨迹yr。通常参考轨迹采用从现在时刻实际值出发的一阶指数形式。它在未来P个时刻的值为:
Figure BDA0001484745290000121
式中,yr(k+i)为参考模型轨迹函数;T为采用周期;α为柔化系数,且0<α=exp(-T/τ)<1;τ为参考轨迹的时间常数。
从式子中可以看到,采用这种形式的参考轨迹,将减小过量的控制作用,使系统的输出能平滑地达到设定值。并且,柔化系数α改变设定值参考轨迹的形状,在其他参数已经确定不变的情况下,增加α会使参考轨迹曲线变化更平缓,从而减小控制量,使被控过程更加平稳,系统的稳定性和鲁棒性也得以提高。相反,减小α则会改善系统的动态响应,但稳定性和鲁棒性将变差。因此,应在两者兼顾的原则下预先设计和在线调整α值。
2.5、滚动优化
首先,假设被控对象实际脉冲相应及其预测模型的脉冲响应分别为:
h=[h1,h2,…,hN]T (18)
Figure BDA0001484745290000122
可以得到闭环预测模型为:
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (20)
可以得到闭环预测控制模型的单步控制输出为:
Figure BDA0001484745290000123
已知电池储能系统的预测模型和闭环校正预测模型分别为:
Figure BDA0001484745290000131
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (23)
并且已知电池储能系统输出的参考轨迹为yr(k+1),设系统误差方程为:
e(k+1)=yr(k+1)-yp(k+1) (24)
优化性能函数为:
Figure BDA0001484745290000132
式中,Q1为非负定加权对称矩阵;Q2为正定控制加权对称矩阵。
Figure BDA0001484745290000133
使得上式优化性能目标函数最小,可以求得最优控制量Δu2为:
Figure BDA0001484745290000134
式中,e'为参考轨迹与零输入响应下闭环预测输入之差,即为:
Figure BDA0001484745290000135
以上所述仅为本发明的较佳实例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。

Claims (4)

1.一种电池储能系统跟踪风光计划出力的控制方法,其特征在于,步骤为:
步骤一,采集风电场、光伏电站的功率计划出力数据;
步骤二,基于第一步采集风电场、光伏电站的功率计划出力数据,建立电池储能系统的功率出力模型;
步骤三,基于第二步建立电池储能系统的功率出力模型,建立参考轨迹模型;
步骤四,基于第三步建立参考轨迹模型,建立闭环预测模型;
步骤五,基于第四步建立闭环预测模型,根据模型预测控制算法计算得出最优控制量;
步骤三中参考轨迹模型建立方法为:
参考轨迹模型预测控制的目标是使电池储能系统的输出y(t)沿着一条预先规定的轨迹逐渐达到设定值ysp;这条预先规定的轨迹称为参考轨迹yr;参考轨迹采用从现在时刻实际值出发的一阶指数形式;它在未来P个时刻的值为:
Figure QLYQS_1
式中,yr(k+i)为参考模型轨迹函数;α为柔化系数,且0<α=exp(-T/τ)<1;τ为参考轨迹的时间常数;
步骤四建立中闭环预测模型方法为:
将第k时刻的实际对象的输出测量值与预测模型输出之间的误差附加到模型的预测输出ym(k+1)上,得到闭环预测模型,用yp(k+1)表示,即:
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (11)
式中:
yp(k+1)=[yp(k+1),yp(k+2),···,yp(k+P)]T,h0=[1,1,…,1]T (12)
如果以i时刻为预测模型的起点,可以列出其闭环预测模型为:
Figure QLYQS_2
已知脉冲响应和阶跃响应系数之间有如下关系:
Figure QLYQS_3
/>
Figure QLYQS_4
展开式稍加整理,然后将
p=[P1,P2,···,PP]T
h0=[1,1,···,1]T
yp(k+1)=[yp(k+1),yp(k+2),···,yp(k+P)]T
Δu(k+1)=[Δu(k),Δu(k+1),···,Δu(k+P-1)]T
Figure QLYQS_5
代入,可以得到闭环预测模型为:
yp(k+1)=h0y(k)+p+AΔu(k+1) (16)
式中,光伏电站出力计划值作为系统的设定值;p为预测步长;Δu为控制变量序列增量。
2.根据权利要求1所述的一种电池储能系统跟踪风光计划出力的控制方法,其特征在于,步骤一中采集数据,求出系统设定值
ys(k)=yf(k)-yw(k) (3)
式中,yf(k)为k时刻系统风电场、光伏电站的计划功率出力值;yw(k)为k时刻系统风电场、光伏电站的实际功率出力值;ys(k)为系统k时刻系统的设定值,它是由电池系统的计划功率出力值与实际功率出力值做差所得。
3.根据权利要求1所述的一种电池储能系统跟踪风光计划出力的控制方法,其特征在于,步骤二中电池储能系统的功率出力模型采用对电池储能系统的电量进行逐时刻积累的方式,进而来表示电池储能系统的充、放电的过程,用电池储能系统的充、放电过程中的效率近似的等效充、放电过程中的电能损耗;电池储能系统不断变化的剩余电量与任意时刻的充放电功率、充放电效率和自放电率相关,剩余电量递推关系如下:
充电过程
E(k)=(1-σsdr)E(k-1)-y(k)ΔtηC (4)
放电过程
E(k)=(1-σsdr)E(k-1)-y(k)Δt/ηD (5)
式中,E(k)为k时刻末该电池储能系统的剩余电池容量,MW·min;y(k)为k时刻电池储能系统的充、放电功率值;σsdr为电池储能系统自放电率;ηC为该电池储能系统的充电效率;ηD为该电池储能系统的放电效率;Δt为计算窗口时长,min;
在电池储能系统充、放电过程中,其允许的最大充电功率和最大放电功率以及剩余容量的边界分别如下:
充电过程、放电过程对于电池储能系统最大出力的限制分别为:
Figure QLYQS_6
Figure QLYQS_7
电池储能系统剩余容量的限制为:
Emin≤E(k)≤Emax (8)
电池储能系统k时刻充放电功率限制由该时刻电池储能系统的剩余电量和电池储能系统自身特性共同决定,具体表达式如下:
电池储能系统k时刻,充电过程对于充电功率限制值为:
Figure QLYQS_8
电池储能系统k时刻,放电过程对于放电功率限制值为:
Figure QLYQS_9
式中,
Figure QLYQS_10
为电池储能系统最大允许充电功率值(MW),由电池储能系统自身特性决定;/>
Figure QLYQS_11
为电池储能系统最大允许放电功率值(MW),由电池储能系统自身特性决定;Emin电池储能系统最小电池储能系统剩余容量限制(MW·h)、Emax电池储能系统最大电池储能系统剩余容量限制(MW·h);ylimit(k)为电池储能介质充、放电功率限制值。
4.根据权利要求1所述的一种电池储能系统跟踪风光计划出力的控制方法,其特征在于,步骤五中滚动优化步骤为:
首先,假设被控对象实际脉冲响应及其预测模型的脉冲响应分别为:
h=[h1,h2,···,hN]T (18)
Figure QLYQS_12
可以得到闭环预测模型为:
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (20)
可以得到闭环预测控制模型的单步控制输出为:
Figure QLYQS_13
已知电池储能系统的预测模型和闭环校正预测模型分别为:
Figure QLYQS_14
yp(k+1)=ym(k+1)+h0[y(k)-ym(k)] (23)
并且已知电池储能系统输出的参考轨迹为yr(k+1),设系统误差方程为:
e(k+1)=yr(k+1)-yp(k+1) (24)
优化性能函数为:
Figure QLYQS_15
式中,Q1为非负定加权对称矩阵;Q2为正定控制加权对称矩阵;
Figure QLYQS_16
使得上式优化性能目标函数最小,可以求得最优控制量Δu2为:
Figure QLYQS_17
式中,e'为参考轨迹与零输入响应下闭环预测输入之差,即为:
Figure QLYQS_18
/>
CN201711211856.3A 2017-11-28 2017-11-28 一种电池储能系统跟踪风光计划出力的控制方法 Active CN107846045B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711211856.3A CN107846045B (zh) 2017-11-28 2017-11-28 一种电池储能系统跟踪风光计划出力的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711211856.3A CN107846045B (zh) 2017-11-28 2017-11-28 一种电池储能系统跟踪风光计划出力的控制方法

Publications (2)

Publication Number Publication Date
CN107846045A CN107846045A (zh) 2018-03-27
CN107846045B true CN107846045B (zh) 2023-06-02

Family

ID=61680583

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711211856.3A Active CN107846045B (zh) 2017-11-28 2017-11-28 一种电池储能系统跟踪风光计划出力的控制方法

Country Status (1)

Country Link
CN (1) CN107846045B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108631368B (zh) * 2018-05-17 2021-09-07 东南大学 计及储能运行损耗下风储系统联合调度的储能配置方法
CN109120003A (zh) * 2018-09-07 2019-01-01 中国南方电网有限责任公司 一种基于mpc算法的分布式光伏储能系统优化控制方法
CN110071512B (zh) * 2019-03-13 2022-07-05 国网浙江平阳县供电有限责任公司 储能装置平抑光伏输出功率波动的控制方法
CN110970913B (zh) * 2019-11-13 2021-09-28 国网内蒙古东部电力有限公司经济技术研究院 一种基于云计算的电池储能系统动态控制方法
CN111478375B (zh) * 2020-05-09 2022-10-14 国网江苏省电力有限公司 计及相关性的风电场、光伏电站群出力非机理建模方法
CN114844027B (zh) * 2022-05-17 2024-04-12 兰州交通大学 一种基于有限控制集模型预测的光伏系统直流母线电压控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102003349A (zh) * 2010-11-05 2011-04-06 甘肃省电力公司 风电受限状态下的大规模风电有功优化方法
CN102522776A (zh) * 2011-12-23 2012-06-27 中国电力科学研究院 储能系统提升风电跟踪计划出力能力的方法
CN105162149A (zh) * 2015-09-07 2015-12-16 中国电力科学研究院 基于模糊自适应调节的光储系统跟踪发电计划出力方法
CA2874269A1 (en) * 2014-12-11 2016-06-11 Meaghan Ann Charest-Finn Model predictive controller and method with correction parameter to compensate for time lag
CN106505635A (zh) * 2016-09-20 2017-03-15 北京恒泰实达科技股份有限公司 弃风最小的有功调度模型及调度系统
CN106786702A (zh) * 2017-01-19 2017-05-31 上海电力学院 用于混合储能系统的全状态模型预测能量调度方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102003349A (zh) * 2010-11-05 2011-04-06 甘肃省电力公司 风电受限状态下的大规模风电有功优化方法
CN102522776A (zh) * 2011-12-23 2012-06-27 中国电力科学研究院 储能系统提升风电跟踪计划出力能力的方法
CA2874269A1 (en) * 2014-12-11 2016-06-11 Meaghan Ann Charest-Finn Model predictive controller and method with correction parameter to compensate for time lag
CN105162149A (zh) * 2015-09-07 2015-12-16 中国电力科学研究院 基于模糊自适应调节的光储系统跟踪发电计划出力方法
CN106505635A (zh) * 2016-09-20 2017-03-15 北京恒泰实达科技股份有限公司 弃风最小的有功调度模型及调度系统
CN106786702A (zh) * 2017-01-19 2017-05-31 上海电力学院 用于混合储能系统的全状态模型预测能量调度方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《一种适用于混合储能系统的控制策略》;丁明 等;《中国电机工程学报》;20120305;全文 *
《基于超短期风电预测功率的储能系统跟踪风电计划出力控制方法》;闫鹤鸣 等;《电网技术》;20150205;全文 *

Also Published As

Publication number Publication date
CN107846045A (zh) 2018-03-27

Similar Documents

Publication Publication Date Title
CN107846045B (zh) 一种电池储能系统跟踪风光计划出力的控制方法
CN105162149B (zh) 基于模糊自适应调节的光储系统跟踪发电计划出力方法
US20200191117A1 (en) Adaptive dynamic planning control method and system for energy storage station, and storage medium
AU2019202751B2 (en) Method for Controlling Power Fluctuation Ramps having Energy Storage Systems in Plants for Intermittent Energy Generation
CN103023074B (zh) 一种基于模型预测控制的大电网有功实时调度方法
CN106253315B (zh) 一种考虑电动汽车充电站可调度性的储能容量配置方法
CN107732897B (zh) 融合虚拟储能系统的楼宇微网模型预测调控方法
CN104993522A (zh) 一种基于mpc的主动配电网多时间尺度协调优化调度方法
CN104779631A (zh) 基于风电预测功率的储能跟踪风电出力计划方法及其系统
CN110165991B (zh) 一种光伏组件跟踪系统控制方法及装置
WO2015154318A1 (zh) 梯级水电站群日发电计划编制方法
CN110581571A (zh) 一种主动配电网动态优化调度方法
JP2006304402A (ja) 分散型エネルギーシステムの制御装置、制御方法、プログラム、および記録媒体
CN109787251A (zh) 一种集群温控负荷聚合模型、系统参数辨识及反推控制方法
CN104103020B (zh) 考虑可变寿命特性的确定电力系统调峰用电池储能容量的方法
CN110739726A (zh) 一种考虑海上风电接入的多类型电源容量长期规划方法
CN104871388A (zh) 控制装置及方法和程序、以及具备该控制装置的自然能发电装置
CN113659620A (zh) 基于动态频率约束的水风互补发电系统日前调度方法
CN104537577A (zh) 一种基于参考线的水电站确定性优化调度方法
CN107958306A (zh) 一种基于参考线的水电站随机优化调度方法
CN111798060A (zh) 一种基于机组爬坡率估计的功率指令优化分配方法
CN112886566A (zh) 一种电网调峰需求计算方法和系统
CN109102110B (zh) 一种径流式小水电短期出力预测方法和装置
CN108631368B (zh) 计及储能运行损耗下风储系统联合调度的储能配置方法
CN110311424B (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