CN109950920B - 一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法 - Google Patents

一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法 Download PDF

Info

Publication number
CN109950920B
CN109950920B CN201910313548.4A CN201910313548A CN109950920B CN 109950920 B CN109950920 B CN 109950920B CN 201910313548 A CN201910313548 A CN 201910313548A CN 109950920 B CN109950920 B CN 109950920B
Authority
CN
China
Prior art keywords
energy storage
electrochemical energy
photovoltaic
power generation
photovoltaic 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
CN201910313548.4A
Other languages
English (en)
Other versions
CN109950920A (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 Jiangsu Electric Power Co ltd Yangzhou Power Supply Branch
Southeast University
State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Jiangsu Electric Power Co ltd Yangzhou Power Supply Branch
Southeast University
State Grid Jiangsu 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 Jiangsu Electric Power Co ltd Yangzhou Power Supply Branch, Southeast University, State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Jiangsu Electric Power Co ltd Yangzhou Power Supply Branch
Priority to CN201910313548.4A priority Critical patent/CN109950920B/zh
Publication of CN109950920A publication Critical patent/CN109950920A/zh
Application granted granted Critical
Publication of CN109950920B publication Critical patent/CN109950920B/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
    • 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
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin

Abstract

一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法,涉及平滑光伏发电技术领域。提供了一种减少光伏发电波动对电网的影响,提高储能资源利用率的方法。通过获取日前光伏发电预测功率和光伏场站侧电化学储能的基本信息,采用动态规划方法根据以上信息建立光伏场站侧电化学储能充放电计划的数学模型,以最小化光伏发电曲线方差为目标确定目标函数,以此计算出光伏站场侧电化学储能充放电功率的最优解。该方法通过改变光伏站场侧电化学储能充放电功率的方式,减少光伏发电曲线波动对电网的影响,提高储能资源利用率,有利于维护电网安全稳定和光伏电站经济运营。

Description

一种基于动态规划算法的光伏场站侧电化学储能充放电计划 制定方法
技术领域
本发明涉及电网调峰和光伏发电技术领域,尤其涉及一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法。
背景技术
光伏发电受太阳辐照度、环境温湿度、通风条件等因素的影响,具有随机性、波动性和间歇性的特点。大规模光伏系统并网会对电网的稳定性和安全性造成冲击。为了减少光伏发电波动对电网的影响,提高储能资源利用率,需要对光伏发电进行合理预测和调度。对光伏发电功率进行预测可作为电网调度的参考,有利于维护电网安全稳定和光伏电站经济运营。
同时,储能电源具有响应速度快、控制精度高的优点,能够在很大程度上改善光伏发电的波动。因此,在合理预测光伏发电功率的基础上,利用电化学储能在光伏发电波动性大时充电或放电,达到平滑光伏发电曲线的效果。利用电化学储能平滑光伏发电功率有利于维护光伏发电的稳定,提高储能设备利用率。
发明内容
本发明专利目的是针对上述背景技术的不足,提供了一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法。
本发明为实现上述发明目的采用如下技术方案:
获取日前光伏发电预测功率和光伏场站侧电化学储能的基本信息,采用动态规划方法根据以上信息建立光伏场站侧电化学储能充放电计划的数学模型,以最小化光伏发电曲线方差为目标确定目标函数,以此计算出光伏场站侧电化学储能充放电功率的最优解。
作为基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法,获取目前光伏发电预测功率的方法为:利用K均值算法(Kmeans)和支持向量机(SVM)的方法预测短期光伏发电功率。根据短期光伏发电特性和光伏发电季节特性,组织预测模型的训练样本集。通过K均值算法对训练样本集进行聚类分析,在聚类得到的各类别数据上分别训练支持向量机。预测时根据预测样本的类别使用相应的支持向量机进行发电功率预测。
作为基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法,建立光伏场站侧电化学储能充放电计划的数学模型的具体方法为:
首先,基于动态规划方法定义电化学储能充放电计划的状态变量为:St=(Ct,Dt),Ct为t时刻电化学储能的容量(MWh),Dt为t时刻光伏发电预测功率(MW)。其次,定义电化学储能充放电计划的决策变量为bt,代表t时刻电化学储能输出功率(MW),bt>0代表电化学储能充电,bt<0代表电化学储能放电。最后,以最小化光伏发电曲线方差为目标初步确定的目标函数为:
Figure BDA0002032329240000021
其中,f(b)为T时段内光伏发电曲线方差的最小值,T时段表示若干小时;Δt为电化学储能充放电的时间间隔,表示电化学储能充放电功率每隔Δt变化一次;ηc为电化学储能充电效率。
当电化学储能荷电状态初值SoCinit和终值SoCfinal确定后,公式因子
Figure BDA0002032329240000022
为常数,用常数M表示,则目标函数可以简化为:
Figure BDA0002032329240000023
作为基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法,求解光伏场站侧电化学储能充放电决策最优解的方法为:将T时段划分为N个阶段,时间间隔为Δt,每一阶段进行一次决策。则最优策略可以看成从初始状态到终值状态的最优路径。
将简化后的目标函数划分为N个阶段,则t时刻的光伏发电曲线方差的阶段指标函数为:
Figure BDA0002032329240000024
其中,ft(St,bt)代表t时刻电化学储能状态变量为St,确定电化学储能决策变量bt时光伏发电曲线方差的阶段指标函数。
为了求解每一时刻以及全局的最小化光伏发电曲线方差的最优解,得出t时刻至t+Δt时刻电化学储能容量状态转移方程为:
Figure BDA0002032329240000025
表示t时刻与t+Δt时刻电化学储能容量的关系。利用贝尔曼方程并基于状态转移方程定义新的最小化光伏发电曲线方差的目标函数为:
f′t(St,bt)=min{ft(St,bt)+f′t+Δt(St+Δt,bt+Δt)};
其中,f′t(St,bt)表示t时刻至最后一个时刻光伏发电曲线方差最小化的目标函数;St+Δt=(Ct+Δt,Dt+Δt),ft+Δt(St+Δt,bt+Δt)表示t+Δt时刻至最后一个时刻光伏发电曲线方差最小化的目标函数。通过递归算法最终得到T时段内每一时刻发电曲线方差指标函数的最小值,从而确定每一时刻电化学储能决策变量的最优解bt
本发明采用上述技术方案,具有以下有益效果:
1)本发明所提出的基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法,根据获取的光伏发电预测功率和光伏场站侧电化学储能的基本信息建立光伏场站侧电化学储能充放电计划的数学模型,以最小化光伏发电曲线方差为目标动态规划出光伏站场侧电化学储能充放电功率的最优决策。不仅考虑了光伏发电的季节特性,还考虑了电化学储能容量和功率等约束,在工作周期内光伏发电曲线方差不断缩小的过程中,能够客观反应电化学储能对平滑光伏发电曲线的有益效果,确定的光伏场站侧电化学储能的充放电计划更加合理;
2)本发明所提出的基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法,使用动态规划算法求解光伏场站侧电化学储能的充放电功率,将一个长时间段的问题划分成多个短时间段的问题,算法灵活度高,不仅降低了问题的求解难度,还反映了电化学储能平滑光伏发电曲线的动态演变过程,提高了计算结果的精确度和准确度,提高了该发明的可信度。
附图说明
图1为本发明光伏场站电化学储能充放电计划制定方法的流程示意图;
图2(a)、图2(b)分别为电化学储能平滑光伏发电功率曲线和电化学储能系统出力曲线。
具体实施方式
下面结合附图对发明的技术方案进行详细说明。
本发明公开的基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法如图1所示,包含以下四个步骤:
步骤1:获取日前光伏发电预测功率
本发明获取T时段内基于K均值算法(Kmeans)和支持向量机(SVM)的短期光伏发电预测功率,T时段表示若干小时。光伏发电受当前气象条件的影响,不同天气,不同时段的发电量差异很大。根据短期光伏发电特性和光伏发电季节特性,组织预测模型的训练样本集。通过K均值算法对训练样本集进行聚类分析,在聚类得到的各类别数据上分别训练支持向量机。预测时根据预测样本的类别使用相应的支持向量机进行发电功率预测。
步骤2:获取电化学储能的基本信息
充放电计划制定方法需要提前获取电化学储能的基本信息,用于建立最小化光伏发电曲线方差的目标函数,电化学储能的基本信息如下:
1)电化学储能额定功率Prated
2)电化学储能额定容量Crated
3)电化学储能充放电效率ηc·ηd
4)电化学储能荷电状态初值和终值SoCinit、SoCfinal
步骤3:基于动态规划建立光伏场站侧电化学储能充放电计划的数学模型
(3-1)定义光伏场站侧电化学储能充放电计划的状态变量St
定义光伏场站侧电化学储能充放电计划的状态变量St=(Ct,Dt)。其中,Ct为t时刻电化学储能的容量(MWh);Dt为t时刻光伏发电预测功率(MW)。
(3-2)定义光伏场站侧电化学储能充放电计划的决策变量bt
bt为t时刻电化学储能输出功率(MW),bt>0代表电化学储能充电,bt<0代表电化学储能放电。决策变量bt需要满足的电化学储能额定容量和额定功率约束如下:
Figure BDA0002032329240000041
其中,Δt为电化学储能充放电的时间间隔,表示电化学储能充放电功率每隔Δt变化一次。C0为电化学储能起始时刻容量,CT为电化学储能结束时刻容量,均满足电化学储能起始和结束时刻荷电状态的约束。
(3-3)确定最小化光伏发电曲线方差的目标函数
为了减少光伏发电对电网的影响,本发明的目标为光伏发电曲线尽量平坦。在数学上,方差可以反映随机变量偏离其均值的程度。曲线的方差可以反映发电曲线的平坦程度。本发明选取最小化光伏发电曲线的方差作为目标函数:
Figure BDA0002032329240000042
其中,f(b)为T时段内光伏发电曲线方差的最小值。
当电化学储能荷电状态初值SoCinit和终值SoCfinal确定后,
Figure BDA0002032329240000043
为常数,用常数M表示,式(1)可以简化为:
Figure BDA0002032329240000051
由式(2)可以看出,目标函数具有分离性,满足动态规划方法对指标函数的要求。
步骤4:求解光伏场站侧电化学储能充放电功率bt
采用动态规划方法求解电化学储能充放电决策的方法如下:将T时段划分为N个阶段,时间间隔为Δt。则最优策略可以看成从初始状态到终值状态的最优路径。则t时刻的光伏发电曲线方差的阶段指标函数为:
Figure BDA0002032329240000052
其中,ft(St,bt)代表t时刻电化学储能状态变量为St,确定电化学储能决策变量bt时光伏发电曲线方差的指标函数。
为了将每一阶段联系起来,可得t时刻至t+Δt时刻电化学储能容量的状态转移方程为;
Figure BDA0002032329240000053
为了求解每一时刻以及全局的最小化光伏发电曲线方差的最优解,基于式(3)和(4),利用贝尔曼最优方程定义新的最小化光伏发电曲线方差的目标函数为:
f′t(St,bt)=min{ft(St,bt)+f′t+Δt(St+Δt,bt+Δt)} (5)
其中,f′t(St,bt)表示t时刻至最后一个时刻光伏发电曲线方差最小化的目标函数,f′t+Δt(St+Δt:bt+Δt)表示t+Δt时刻至最后一个时刻光伏发电曲线方差最小化的目标函数。公式(4)表明了当前时刻与下一时刻目标函数之间的关系,通过递归算法最终得到T时段内每一时刻发电曲线方差指标函数的最小值,从而确定每一时刻电化学储能决策变量的最优解bt
下面结合一实施例对本发明技术方案的有益效果进行说明。
假设电化学储能额定容量为20MWh,额定放电功率为5MW,充放电效率分别为0.9/0.95,电化学储能起始荷电状态SoCinit和结束荷电状态SoCfinal均为0。获取浙江绍兴地区某日24小时光伏电站发电功率,曲线在上午、晚上各有一次高峰,在凌晨、中午各有一次低谷,Δt=30min。
利用电化学储能平滑光伏发电曲线的优化结果如图2所示。其中,图2(a)虚线为无储能参与的原始光伏发电功率曲线,实线为有储能参与的优化后的发电功率曲线;图2(b)为电化学储能系统出力,原始发电曲线方差为14.36,优化后曲线方差为2.01。由图可知,电化学储能在光伏发电功率高峰期放电,在低谷期充电,以达到平滑光伏发电曲线的目的。因此,有电化学储能参与的发电功率曲线更加平滑,储能参与后的方差比原始光伏发电曲线方差明显减小。

Claims (1)

1.一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法,其特征在于,获取目前光伏发电预测功率和光伏场站侧电化学储能的基本信息,采用动态规划方法根据以上信息建立光伏场站侧电化学储能充放电计划的数学模型,以此计算出光伏站场侧电化学储能充放电功率的最优解;
以最小化光伏发电曲线方差为目标动态规划出光伏站场侧电化学储能的充放电功率的方法为:
1)、建立光伏场站侧电化学储能充放电计划的数学模型,以最小化光伏发电曲线方差为目标确定目标函数;
2)、将运行周期划分为若干个阶段,确定电化学储能运行周期内每一时刻光伏发电曲线方差的贝尔曼最优方程,根据电化学储能运行周期内每一时刻的贝尔曼最优方程并在各时刻电化学储能容量和功率的约束下计算出对应时刻电化学储能充放电功率的最优值;
建立光伏场站侧电化学储能充放电计划的数学模型,以最小化光伏发电曲线方差为目标确定的目标函数为:
Figure FDA0002549503500000011
其中,f(b)为T时段内光伏发电曲线方差的最小值,T时段表示若干小时,b表示电化学储能;Δt为电化学储能充放电的时间间隔;Dt为t时刻光伏发电预测功率;bt为t时刻电化学储能输出功率,bt>0代表电化学储能充电,bt<0代表电化学储能放电;ηc为电化学储能充电效率;
当电化学储能起始和结束的荷电状态确定后,
Figure FDA0002549503500000012
为常数,用常数M表示,则上述最小化光伏发电曲线方差的目标函数简化为:
Figure FDA0002549503500000013
将T时段划分为N个阶段,得到t时刻光伏发电曲线方差的阶段指标函数为:
Figure FDA0002549503500000014
ft(St,bt)代表电化学储能在t时刻的状态St下电化学储能充放电功率为bt时光伏发电曲线方差的阶段指标函数,St=(Ct,Dt),Ct为t时刻电化学储能的容量;
t时刻到t+Δt时刻电化学储能容量状态转移方程为:
Figure FDA0002549503500000021
表示t时刻与t+Δt时刻电化学储能容量的关系;
电化学储能运行周期内每一时刻光伏发电曲线方差的贝尔曼最优方程为:ft′(St,bt)=min{ft(St,bt)+f′t+Δt(St+Δt,bt+Δt)}
f′t(St,bt)表示t时刻至最后一个时刻光伏发电曲线方差最小化的目标函数;St+Δt=(Ct+Δt,Dt+Δt),f′t+Δt(St+Δt,bt+Δt)表示t+Δt时刻至最后一个时刻光伏发电曲线方差最小化的目标函数;
各时刻电化学储能充放电功率取值约束有:0≤Ct+bt×Δt≤Crated,|bt|≤Prated
Figure FDA0002549503500000022
Crated为电化学储能的额定容量,Prated为电化学储能的额定功率,C0为电化学储能起始时刻容量,CT为电化学储能结束时刻容量,SoCinit为电化学储能荷电状态初值,SoCfinal为电化学储能荷电状态终值。
CN201910313548.4A 2019-04-18 2019-04-18 一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法 Active CN109950920B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910313548.4A CN109950920B (zh) 2019-04-18 2019-04-18 一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910313548.4A CN109950920B (zh) 2019-04-18 2019-04-18 一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法

Publications (2)

Publication Number Publication Date
CN109950920A CN109950920A (zh) 2019-06-28
CN109950920B true CN109950920B (zh) 2020-08-14

Family

ID=67015453

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910313548.4A Active CN109950920B (zh) 2019-04-18 2019-04-18 一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法

Country Status (1)

Country Link
CN (1) CN109950920B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114925926B (zh) * 2022-05-31 2023-04-07 西安交通大学 一种具有多项式时间复杂度的直流通道规划与运行方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102624017A (zh) * 2012-03-22 2012-08-01 清华大学 一种基于负荷预测的电池储能系统削峰填谷实时控制方法
CN103248065A (zh) * 2013-05-22 2013-08-14 中国南方电网有限责任公司调峰调频发电公司 一种风电场中电池储能系统的充放电控制方法
CN108599206A (zh) * 2018-04-28 2018-09-28 国网湖南省电力有限公司 高比例不确定电源场景下的配电网混合储能配置方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016082617A (ja) * 2014-10-10 2016-05-16 富士通株式会社 需要調整計画生成装置、需要調整計画生成方法、及び需要調整計画生成プログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102624017A (zh) * 2012-03-22 2012-08-01 清华大学 一种基于负荷预测的电池储能系统削峰填谷实时控制方法
CN103248065A (zh) * 2013-05-22 2013-08-14 中国南方电网有限责任公司调峰调频发电公司 一种风电场中电池储能系统的充放电控制方法
CN108599206A (zh) * 2018-04-28 2018-09-28 国网湖南省电力有限公司 高比例不确定电源场景下的配电网混合储能配置方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"含光伏发电集中型充电站充放电策略的研究";田昊等;《可再生能源》;20131031;第31卷(第10期);第34-38页 *

Also Published As

Publication number Publication date
CN109950920A (zh) 2019-06-28

Similar Documents

Publication Publication Date Title
CN109103912B (zh) 考虑电网调峰需求的工业园区主动配电系统调度优化方法
CN111525627B (zh) 含抽蓄与新能源发电的柔性直流输电系统日前调度方法
GB2596240A (en) Coordinated control of renewable electric generation resource and charge storage device
CN105680474B (zh) 一种储能抑制光伏电站快速功率变化的控制方法
CN110544935A (zh) 一种电-氢多能互补直流微电网协调调度方法
CN106447131B (zh) 独立型微电网能量调控方法
CN107069776B (zh) 一种平滑微网联络线功率的储能前瞻分布式控制方法
CN110350589B (zh) 一种可再生能源和储能调度模型及调度方法
CN111934360B (zh) 基于模型预测控制的虚拟电厂-储能系统能量协同优化调控方法
CN111555356B (zh) 基于耦合Markov模型的区域综合能源系统随机变量分析方法
CN113824111A (zh) 一种光储能场景中储能容量配置和调度方法
CN114156951B (zh) 一种源网荷储系统的控制优化方法和装置
CN109950920B (zh) 一种基于动态规划算法的光伏场站侧电化学储能充放电计划制定方法
CN111654054A (zh) 一种储能基于自适应神经网络(ann)平抑短期风电波动的控制方法
CN114154790A (zh) 基于需量管理和柔性负荷的工业园区光储容量配置方法
CN115619153A (zh) 动态聚合下考虑分布式资源运行特性的响应方法和装置
CN114172180B (zh) 考虑退化成本的电池储能参与调峰双层滚动优化控制方法
CN111311032B (zh) 一种基于扇形雷达图模型的微网系统容量优化配置方法
CN109347120A (zh) 一种面向电网频率响应的决策方法和系统
CN110717694B (zh) 基于新能源消纳期望值的储能配置随机决策方法及装置
CN114723278A (zh) 一种考虑光伏储能的社区微电网调度方法及系统
CN114498690A (zh) 支持大规模可再生能源消纳的多元复合储能优化配置方法
CN107134773B (zh) 一种基于事件的智能微电网运行优化方法
CN110783935B (zh) 一种限定荷电状态下电池储能响应调频信号的方法
CN116388233B (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