CN110504701A - 一种考虑动态充放电效率的电池储能***调度方法 - Google Patents
一种考虑动态充放电效率的电池储能***调度方法 Download PDFInfo
- Publication number
- CN110504701A CN110504701A CN201910778785.8A CN201910778785A CN110504701A CN 110504701 A CN110504701 A CN 110504701A CN 201910778785 A CN201910778785 A CN 201910778785A CN 110504701 A CN110504701 A CN 110504701A
- Authority
- CN
- China
- Prior art keywords
- energy storage
- storage system
- battery energy
- electricity
- formula
- 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.)
- Granted
Links
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/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- 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
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/10—Flexible AC transmission systems [FACTS]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
一种考虑动态充放电效率的电池储能***调度方法,包括以下步骤:A、获取电池储能***基本信息,对充放电效率特性曲线进行分段线性化;B、获取电池储能***运行限制、初末电量、预测分时电价或价格乘子等调度信息;C、建立调度模型;D、计算电池储能***各时段存储电量的动态上下限;E、动态规划反向递推,获得各时段需要反复求解的核心子问题的优化模型;F、利用核心子问题优化模型的几何特征,从下到上依次求解各水平子区域中最优路径;G、求解水平子区域内的最优路径;H、处理最优路径跳变;I、动态规划正向递推,计算各时段的最优充放电电量及功率决策;J、随时间发展进入下一时段,更新步骤B中的调度信息,执行步骤C。本方法为电池储能***的优化运行提供了一种有效途径,具有工程实用价值。
Description
技术领域
本发明涉及电池储能技术领域,尤其涉及一种考虑动态充放电效率的电池储能***调度方法。
背景技术
可再生能源的不断快速接入电网,其不确定性对电力***的控制和管理提出了新的挑战。近年来,储能已成为新的行业热点。通过储能***的广泛应用,不仅能够平抑功率波动,提高电能质量,增强电网***运行的可靠性。
一般情况下,储能的充放电效率随着充放电功率动态变化。例如铅酸电池和锂电池储能***在放电时荷电状态(SOC)处于20%到95%,在这个范围内,放电时电压基本不变,效率随着放电电流的增大会有一定程度的降低。充电过程视为放电的逆过程,效率变化和放电类似,在偏离额定工况时,损耗增大,效率降低。
目前已有电池储能***调度方法通常是根据未来一天的预测数据,制定储能的运行计划,在实际运行偏离运行计划时根据设定量进行修正,优化电池储能***的运行效益。但是目前在求解电池储能***的调度模型时一般都采用通用的商业优化软件进行求解,缺少有针对性、适应性强的专用求解算法。
因此,目前迫切需要有一种能够针对调度模型自身结构特性的专用算法能够实现调度问题的高效求解。
发明内容
本发明的目的是提供一种考虑动态充放电效率的电池储能***调度方法,本发明方法为电池储能***的优化运行提供了一种有效途径,具有工程实用价值。
为实现上述目的,本发明采用的技术方案包括以下步骤:
一种考虑动态充放电效率的电池储能***调度方法,运用该方法能够制定电池储能***的优化运行决策,提高电池储能***的运行经济性,该方法包括以下步骤:
A、确定电池储能***的容量、动态充放电功率和效率关系、充放电功率限制和电池储能***容量上下限,根据调度的精度要求选择合适分段数对电池储能***的充放电功率—效率关系进行分段线性化;
B、获取电池储能***的初始存储电量、调度结束时的目标存储电量、各时段上网电量约束以及分时预测电价或者等效价值乘子;
C、建立电池储能***的调度模型,调度模型的优化目标是制定电池储能***各时段充放电功率决策,使电池储能***的总发电收益最大化,调度模型的目标函数为:
式中,T为调度时段总数,t为调度时段编号,λt为第t时段的预测电价或者等效价值乘子,qt为第t时段的电池储能***的充放电电量,放电为正值,充电为负值,p(qt)第t时段的电池储能***的充放电功率,为qt的一元分段线性函数,δ为调度时段长度;
调度模型的约束条件包含:
1)电池储能***的电量平衡约束
式中,et为第t时段末电池储能***的存储电量,Lt为已知的第t时段电池储能***的固定损失电量,为负值;
2)上网电量约束
式中,和分别为第t时段电池储能***的上网电量上下限;
3)电池储能***的充放电功率约束
式中,Pmax和Pmin分别为电池储能***放电功率和充电功率的物理上限,函数p是单调的分段线性函数,充放电功率约束(4)能够转换为充放电量的约束:
式中,函数p'为函数p的逆函数,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限;
4)电量SOC约束
式中,Emax和Emin分别为电池储能***的存储电量的上下限;
5)电池储能***的初末电量约束
e0=E0,eT=ET (7)
式中,E0和ET分别是电池储能***初始时刻电量和调度结束时的目标电量,e0和eT分别为调度开始时刻和第T时段末的电池储能***的存储电量;
由(1),(2),(3),(5),(6),(7)式构成了电池储能***的调度模型;
D、计算电池储能***各时段存储电量的动态上下限约束,根据调度末电池储能***调度结束时目标存储电量ET和电池储能***各时段存储电量上下限约束(2)和(5),通过从T时段开始反向递推,计算出电池储能***各时段存储电量的动态上下限为:
式中,j为时间变量,和分别为反向推导获得的第t时段末电池储能***各时段存储电量的动态上下限,ET为调度结束时电池储能***的目标存储电量,Lj为已知的第j时段电池储能***的固定损失电量,为负值,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限,和分别为第j时段电池储能***上网电量的上下限;
根据电池储能***的初始电量E0和电池储能***各时段存储电量上下限约束(2)和(5),从调度开始时段进行正向递推,计算出电池储能***各时段存储电量的动态上下限为:
式中,j为时间变量,和分别为正向推导获得的第t时段末的电池储能***各时段存储电量的动态上下限,E0为调度开始时电池储能***的存储电量,Lj为已知的第j时段电池储能***的固定损失电量,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限,和分别为第j时段电池储能***上网电量的上下限;
将正向推导和反向推导结果结合,计算电池储能***的存储电量的动态上下限为:
式中,和et分别为计算获得的第t时段的电池储能***各时段存储电量的动态上下限,结合(6)式,得到新的电池储能***各时段存储电量的上下限为:
式中,et为第t时段末电池储能***的存储电量,Emax和Emin分别为电池储能***的存储电量的物理上下限;
经过等效变换,由(1),(2),(3),(5),(7),(13)式构成了电池储能***的调度模型;
E、动态规划反向递推阶段,根据动态规划原理,得到
式中,ft(et)为从第t个调度时段末到调度结束时的最大收益,置fT(ET)=0,将电池储能***的电量平衡约束(2)带入(14),得到需要在各个时段反复求解的动态规划cost-to-go函数:
(15)式在求解过程中受到(2),(3),(5),(7)和(13)式的约束,将(15)式定义为动态规划反向递推阶段中核心子问题,该核心子问题的优化模型由(2),(3),(5),(7),(13)和(15)式构成;
F、对(15)式进行变量代换,将决策量qt+1表示为x,et表示为y,最大收益ft(et)表示为h(x),同时将ft+1(et-Lt+1-qt+1)平移Lt+1得到函数g(y-x),动态规划反向递推阶段中核心子问题的优化模型等效变换为;
式中,a和b分别表示由(3)和(5)式共同确定的第t+1时段充放电量qt+1的上下限,c和d分别表示由(7)和(13)式共同确定的第t时段末电池储能***各时段存储电量et的上下限,根据(2)式,α和β实际上分别表示电池储能***在第t+1时段末存储电量的上下限;
根据(16)式的几何意义,函数h(x)和g(y-x)的分段转折线将定义域剖分成若干子区域,在每个平行四边形或类似于平行四边形的子区域,能够证明最优值一定是在边界处取得,函数值相对于参变量x都是线性变化的,即边界所在直线关于函数值的分段转折点只能在平行四边形子区域的端点取得,根据各端点y值将各平行四边形子区域的端点进行排序,相邻端点所在的水平线将问题的可行域分成了N个水平子区域,从下n=1到上n=N依次寻找各水平子区域中最优路径;
G、水平子区域最优路径求解,从n=1开始,水平子区域内有L条90°垂直线和45°斜线,分别对应的是函数h(x)和g(y-x)的分段转折线,求解这这些分段转折线和水平子区域上下边界交点函数值分别为zl(yu)和zl(yd),找到下边界函数值最大点zl'(yd)及对应的路径l',判断这些路径和上边界交点的函数值zl'(yu)是否为所有上边界交点函数值的最大点,如果不是,则说明在该区域内出现了最优路径的跳变,执行步骤H;如果是,则该条路径是这个水平子区域内的最大路径,判断是否n>N,如果是,则第t时段核心子问题求解流程结束,执行步骤I,如果不是,则置n=n+1,继续求解该水平子区域上方的临近子区域的最优路径;
H、最优路径跳变是指下边界的函数值最大的点通过最初的最优路径上升,而在某转折点处跳变到另外一条最优路径继续上升而后到达了上边界函数值最大的点,将l'作为基准路径,初始化l'=l's,s=1,计算最早与l'相交路径的交点及斜率,获得其相交路径l's+1的转折点坐标(ys+1,zl's+1(ys+1))以及l's+1与上边界交点的函数值zl's+1(yu),判断zl's+1(yu)是否为上边界最大值,如果不是,令s=s+1,将将l's+1作为新的基准路径,重新执行H;如果是,则输出各转折点坐标(ys+1,zl's+1(ys+1))及zl'(yd),zl's(yu)值,得到最优路径l's,s=1,2,...,S,该水平子区域最优路径求解结束,判断是否n>N,如果是,则第t时段核心子问题求解流程结束,执行步骤I,如果不是,则置n=n+1,继续求解该水平子区域上方的临近子区域的最优路径;
I、置t=t-1,判断如果t>0,返回执行步骤E,继续求解t时段的核心子问题,如果t=0,执行步骤J;
J、动态规划正向递推阶段,令对t=1,2,...,T依次计算:
式中,qt *为电池储能***第t时段的最优充放电电量,为第t时段的电池储能***的最优存储电量,Lt为已知的第t时段电池储能***的固定损失电量;求解(17)式时受(3)-(6)式约束,根据动态规划原理,求解(3)-(6)和(17)式定义的极大化问题获得的从第1时段到第T时段的电池储能***的充放电量决策及对应的充放电功率就是电池储能***调度问题的最优解;
J、随时间发展进入下一时段,更新步骤B中电池储能***的初始存储电量、调度结束时的目标存储电量、各时段上网电量约束以及分时预测电价或者等效价值乘子,继续执行步骤C。
和现有技术相比较,本发明具备如下优点:
本发明在建立模型时对电池储能***的非线性的充放电效率特性曲线进行了分段线性化,建立了混合整数线性规划模型,比传统的模型更加精确;此外,在求解混合整数线性规划模型时一般都只能采用CPLEX等商业软件包才能进行高效求解,本发明首选通过步骤D计算电池储能***各时段存储电量的动态上下限约束,缩小了解的可行域范围,减少求解时间;其次,在步骤F中针对动态规划需要反复求解的核心子问题的几何结构特性,在步骤G和步骤H中提出了能够高效求解,同时能够充分保证计算过程的数值稳定性的精细化算法,从而在求解这类电池储能***调度问题时能够替代商业软件。
附图说明
图1为整个调度算法流程图。
图2为电池放电效率-功率曲线。
图3为动态规划算法中核心子问题的几何结构特征。
图4为处理“最优路径跳变”子方法示意图。
图5为处理“最优路径跳变”子方法算法流程图。
具体实施方式
本发明旨在提供一种考虑动态充放电效率特性的电池储能***调度方法,能根据电池储能自身特性、运行约束及预测的市场电价制定未来一段时间内电池储能***运行的最优决策,从而指导电池储能***运行,提高***运行效率及运行收益。
下面结合附图1对发明中的调度方法的实施方式及整个流程做进一步的说明。
一种考虑动态充放电效率的电池储能***调度方法,运用该方法能够制定储能电池***的优化运行决策,提高电池储能***的运行经济性,该方法包括以下步骤:
A、确定电池储能***的容量、动态充放电效率、充放电功率限制和电池储能***容量上下限,参见附图2,对电池储能***的充放电功率—效率关系图进行分段线性化,根据精度要求选择合适分段数。
B、获取电池储能***的初始存储电量、调度结束时的目标存储电量、各时段上网电量约束和分时预测电价或者等效价值乘子。
C、建立电池储能***的调度模型,调度模型的优化目标是制定电池储能***各时段充放电功率决策,使电池储能***的总发电收益最大化,模型的目标函数为:
式中,T为调度时段总数,t为调度时段编号,λt为第t时段的预测电价或等效价值乘子,qt为第t时段的电池储能***的充放电电量,放电为正值,充电为负值,p(qt)第t时段的电池充放电功率,为qt的一元分段线性函数,δ为调度时段长度。
调度模型约束条件包含:
1)电池储能***的电量平衡约束
式中et为第t时段末电池储能***的存储电量,Lt为已知的第t时段电池储能***的存储的固定损失电量,为负值。
2)上网电量约束
式中和分别为第t时段电池储能***的上网电量上下限。
3)电池储能***的充放电功率约束
式中Pmax和Pmin分别为电池储能***放电功率和充电功率的物理上限,函数p是单调的分段线性函数,充放电功率约束(4)能够转换为充放电量的约束:
式中函数p'为函数p的逆函数,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限。
4)电量(SOC)约束
式中Emax和Emin分别为电池储能***的存储电量的上下限。
5)电池储能***的初末电量约束
e0=E0,eT=ET (7)
式中E0和ET分别是电池储能***初始时刻电量和调度结束时的目标电量,e0和eT分别为调度开始时刻和第T时段末的电池储能***的存储电量。
由(1),(2),(3),(5),(6),(7)式构成了电池储能***的调度模型。
D、计算电池储能***各时段存储电量的动态上下限约束,根据调度末电池储能***调度结束时目标存储电量ET和电池储能***各时段存储电量上下限约束(2)和(5),通过从T时段开始反向递推,计算出电池储能***各时段存储电量的动态上下限为:
式中j为时间变量,和分别为反向推导获得的第t时段末电池储能***各时段存储电量的动态上下限,ET为调度结束时电池储能***的目标存储电量,Lj为已知的第j时段电池储能***的固定损失电量,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限,和分别为第j时段电池储能***上网电量的上下限。
根据电池储能***的初始电量E0和电池储能***各时段存储电量上下限约束(2)和(5),从调度开始时段进行正向递推,计算出电池储能***各时段存储电量的动态上下限为:
式中j为时间变量,和分别为正向推导获得的第t时段末的电池储能***各时段存储电量的动态上下限,E0为调度开始时电池储能***的存储电量,Lj为已知的第j时段电池储能***的固定损失电量,为负值,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限,和分别为第j时段电池储能***上网电量的上下限。
将正向推导和反向推导结果结合,计算电池储能***的存储电量的动态上下限为:
式中和e t分别为计算获得的第t时段的电池储能***各时段存储电量的动态上下限,结合(6)式,得到新的电池储能***各时段存储电量的上下限为:
式中,et为第t时段末电池储能***的存储电量,Emax和Emin分别为电池储能***的存储电量的物理上下限;
经过等效变换,由(1),(2),(3),(5),(7),(13)式构成了电池储能***的调度模型。
E、动态规划反向递推阶段,根据动态规划原理,得到
式中ft(et)为从第t个调度时段末到调度结束时的最大收益。置fT(ET)=0,将电池储能***的电量平衡约束(2)带入(14),得到需要在各个时段反复求解的动态规划cost-to-go函数:
(15)式在求解过程中受到(2),(3),(5),(7)和(13)式的约束,将(15)式定义为动态规划反向递推阶段中核心子问题,该核心子问题的优化模型由(2),(3),(5),(7),(13)和(15)式构成。
F、对(15)式进行变量代换,将决策量qt+1表示为x,et表示为y,最大收益ft(et)表示为h(x),同时将ft+1(et-Lt+1-qt+1)平移Lt+1得到函数g(y-x),动态规划反向递推阶段中核心子问题的优化模型等效变换为:
式中,a和b分别表示由(3)和(5)式共同确定的第t+1时段充放电量qt+1的上下限,c和d分别表示由(7)和(13)式共同确定的第t时段末电池储能***各时段存储电量et的上下限,根据(2)式,α和β实际上分别表示电池储能***在第t+1时段末存储电量的上下限。
参见附图3,根据核心子问题优化模型的几何意义,函数h(x)的分段转折线为x=xθ,θ=1,2,3...Φ,Φ为函数h(x)为分段数,对应图中的90°数值点划线;g(y-x)的分段转折线为y-x=ρν,v=1,2,3.....V,V为函数g(y-x)的分段数,对应图中的45°虚斜线。函数h(x)和g(y-x)的分段转折线将定义域剖分成若干子区域,在每个平行四边形或类似于平行四边形的子区域,能够证明最优值一定是在边界处取得,函数值相对于参变量x都是线性变化的,即边界所在直线关于函数值的分段转折点只能在平行四边形子区域的端点取得,根据各端点y值将各平行四边形子区域的端点进行排序,相邻端点所在的水平线将问题的可行域分成了N个水平子区域,从下n=1到上n=N依次寻找各水平子区域中最优路径。
G、水平子区域最优路径求解,参见附图4,从n=1开始,水平子区域内共有L条90°垂直线和45°斜线,分别对应的是函数h(x)和g(y-x)的分段转折线,求解这这些分段转折线和水平子区域上下边界交点函数值分别zl(yu)和zl(yd),找到下边界函数值最大点zl'(yd)及对应的路径l',判断这些路径和上边界交点的函数值zl'(yu)是否为所有上边界交点函数值的最大点,如果不是,则说明在该区域内出现了最优路径的跳变,执行步骤H;如果是,则该条路径是这个水平子区域内的最大路径,判断是否n>N,如果是,则第t时段核心子问题求解流程结束,执行步骤I,如果不是,则置n=n+1,继续求解该水平子区域上方的临近子区域的最优路径。
H、最优路径跳变是指下边界的函数值最大的点通过最初的最优路径上升,而在某转折点处跳变到另外一条最优路径继续上升而后到达了上边界函数值最大的点;参见附图5,通过一个简单例子说明处理最优路径“跳变”方法的思想:分别记录水平子区域中上下边界和所有l条路径交点的函数值分别为zl(yu)和zl(yd),能够得到在zoy平面上路径l对应的线性函数zl(y),从下边界函数值最大的a点出发,在与其它路径相交之前l1为最优路径,而当其和路径l3相交在b点之后,由于l3的上升速率大于原最优路径l1,因此从b点开始最优路径变为l3,同理在c点最优路径又变为l4直至到达上边界函数值最大的d点;详细计算流程参见附图4,将l'作为基准路径,初始化l'=l's,s=1,计算最早与l'相交路径的交点及斜率,获得其相交路径l's+1的转折点坐标(ys+1,zl's+1(ys+1))以及l's+1与上边界交点的函数值zl's+1(yu),判断zl's+1(yu)是否为上边界最大值,如果不是,令s=s+1,将将l's+1作为新的基准路径,重新执行H;如果是,则输出各转折点坐标(ys+1,zl's+1(ys+1))及zl'(yd),zl's(yu)值,得到最优路径l's,s=1,2,...,S,该水平子区域最优路径求解结束,判断是否n>N,如果是,则第t时段核心子问题求解流程结束,执行步骤I,如果不是,则置n=n+1,继续求解该水平子区域上方的临近子区域的最优路径。
I、置t=t-1,判断如果t>0,返回执行步骤E,继续求解t阶段的核心子问题,如果t=0,执行步骤J;
J、动态规划正向递推阶段,令对t=1,2,...,T依次计算:
式中,qt *为电池储能***第t时段的最优充放电电量,为第t时段的电池储能***的最优存储电量,Lt为已知的第t时段电池储能***的固定损失电量;求解(17)式时受(3)-(6)式约束,根据动态规划原理,求解(3)-(6)和(17)式定义的极大化问题获得的从第1时段到第T时段的电池储能***的充放电量决策及对应的充放电功率就是电池储能***调度问题的最优解。
J、随时间发展进入下一时段,更新步骤B中电池储能***的初始存储电量、调度结束时的目标存储电量、各时段上网电量约束以及分时预测电价或者等效价值乘子,继续执行步骤C。
Claims (1)
1.一种考虑动态充放电效率的电池储能***调度方法,其特征在于:运用该方法能够制定电池储能***的优化运行决策,提高电池储能***的运行经济性,该方法包括以下步骤:
A、确定电池储能***的容量、动态充放电效率、充放电功率限制和电池储能***容量上下限,根据调度的精度要求选择合适分段数对电池储能***的充放电功率—效率关系进行分段线性化;
B、获取电池储能***的初始存储电量、调度结束时的目标存储电量、各时段上网电量约束以及分时预测电价或者等效价值乘子;
C、建立电池储能***的调度模型,调度模型的优化目标是制定电池储能***各时段充放电功率决策,使电池储能***的总发电收益最大化,调度模型的目标函数为:
式中,T为调度时段总数,t为调度时段编号,λt为第t时段的预测电价或者等效价值乘子,qt为第t时段的电池储能***的充放电电量,放电为正值,充电为负值,p(qt)第t时段的电池储能***的充放电功率,为qt的一元分段线性函数,δ为调度时段长度;
调度模型的约束条件包含:
1)电池储能***的电量平衡约束
式中,et为第t时段末电池储能***的存储电量,Lt为已知的第t时段电池储能***的固定损失电量,为负值;
2)上网电量约束
式中,和分别为第t时段电池储能***的上网电量上下限;
3)电池储能***的充放电功率约束
式中,Pmax和Pmin分别为电池储能***放电功率和充电功率的物理上限,函数p是单调的分段线性函数,充放电功率约束(4)能够转换为充放电量的约束:
式中,函数p'为函数p的逆函数,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限;
4)电量SOC约束
式中,Emax和Emin分别为电池储能***的存储电量的上下限;
5)电池储能***的初末电量约束
e0=E0,eT=ET (7)
式中,E0和ET分别是电池储能***初始时刻电量和调度结束时的目标电量,e0和eT分别为调度开始时刻和第T时段末的电池储能***的存储电量;
由(1),(2),(3),(5),(6),(7)式构成了电池储能***的调度模型;
D、计算电池储能***各时段存储电量的动态上下限约束,根据调度末电池储能***调度结束时目标存储电量ET和电池储能***各时段存储电量上下限约束(2)和(5),通过从T时段开始反向递推,计算出电池储能***各时段存储电量的动态上下限为:
式中,j为时间变量,和分别为反向推导获得的第t时段末电池储能***各时段存储电量的动态上下限,ET为调度结束时电池储能***的目标存储电量,Lj为已知的第j时段电池储能***的固定损失电量,为负值,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限,和分别为第j时段电池储能***上网电量的上下限;
根据电池储能***的初始电量E0和电池储能***各时段存储电量上下限约束(2)和(5),从调度开始时段进行正向递推,计算出电池储能***各时段存储电量的动态上下限为:
式中,j为时间变量,和分别为正向推导获得的第t时段末的电池储能***各时段存储电量的动态上下限,E0为调度开始时电池储能***的存储电量,Lj为已知的第j时段电池储能***的固定损失电量,p'(Pmax)和p'(Pmin)分别为电池储能***放电电量和充电电量的物理上限,和分别为第j时段电池储能***上网电量的上下限;
将正向推导和反向推导结果结合,计算电池储能***的存储电量的动态上下限为:
式中,和e t分别为计算获得的第t时段的电池储能***各时段存储电量的动态上下限,结合(6)式,得到新的电池储能***各时段存储电量的上下限为:
式中,et为第t时段末电池储能***的存储电量,Emax和Emin分别为电池储能***的存储电量的物理上下限;
经过等效变换,由(1),(2),(3),(5),(7),(13)式构成了电池储能***的调度模型;
E、动态规划反向递推阶段,根据动态规划原理,得到
式中,ft(et)为从第t个调度时段末到调度结束时的最大收益,置fT(ET)=0,将电池储能***的电量平衡约束(2)带入(14),得到需要在各个时段反复求解的动态规划cost-to-go函数:
(15)式在求解过程中受到(2),(3),(5),(7)和(13)式的约束,将(15)式定义为动态规划反向递推阶段中核心子问题,该核心子问题的优化模型由(2),(3),(5),(7),(13)和(15)式构成;
F、对(15)式进行变量代换,将决策量qt+1表示为x,et表示为y,最大收益ft(et)表示为h(x),同时将ft+1(et-Lt+1-qt+1)平移Lt+1得到函数g(y-x),动态规划反向递推阶段中核心子问题的优化模型等效变换为;
式中,a和b分别表示由(3)和(5)式共同确定的第t+1时段充放电量qt+1的上下限,c和d分别表示由(7)和(13)式共同确定的第t时段末电池储能***各时段存储电量et的上下限,根据(2)式,α和β实际上分别表示电池储能***在第t+1时段末存储电量的上下限;
根据(16)式的几何意义,函数h(x)和g(y-x)的分段转折线将定义域剖分成若干子区域,在每个平行四边形或类似于平行四边形的子区域,能够证明最优值一定是在边界处取得,函数值相对于参变量x都是线性变化的,即边界所在直线关于函数值的分段转折点只能在平行四边形子区域的端点取得,根据各端点y值将各平行四边形子区域的端点进行排序,相邻端点所在的水平线将问题的可行域分成了N个水平子区域,从下n=1到上n=N依次寻找各水平子区域中最优路径;
G、水平子区域最优路径求解,从n=1开始,水平子区域内有L条90°垂直线和45°斜线,分别对应的是函数h(x)和g(y-x)的分段转折线,求解这这些分段转折线和水平子区域上下边界交点函数值分别为zl(yu)和zl(yd),找到下边界函数值最大点zl'(yd)及对应的路径l',判断这些路径和上边界交点的函数值zl'(yu)是否为所有上边界交点函数值的最大点,如果不是,则说明在该区域内出现了最优路径的跳变,执行步骤H;如果是,则该条路径是这个水平子区域内的最大路径,判断是否n>N,如果是,则第t时段核心子问题求解流程结束,执行步骤I,如果不是,则置n=n+1,继续求解该水平子区域上方的临近子区域的最优路径;
H、最优路径跳变是指下边界的函数值最大的点通过最初的最优路径上升,而在某转折点处跳变到另外一条最优路径继续上升而后到达了上边界函数值最大的点,将l'作为基准路径,初始化l'=l's,s=1,计算最早与l'相交路径的交点及斜率,获得其相交路径l's+1的转折点坐标(ys+1,zl's+1(ys+1))以及l's+1与上边界交点的函数值zl's+1(yu),判断zl's+1(yu)是否为上边界最大值,如果不是,令s=s+1,将将l's+1作为新的基准路径,重新执行H;如果是,则输出各转折点坐标(ys+1,zl's+1(ys+1))及zl'(yd),zl's(yu)值,得到最优路径l's,s=1,2,...,S,S为最优路径总数,该水平子区域最优路径求解结束,判断是否n>N,如果是,则第t时段核心子问题求解流程结束,执行步骤I,如果不是,则置n=n+1,继续求解该水平子区域上方的临近子区域的最优路径;
I、置t=t-1,判断如果t>0,返回执行步骤E,继续求解t时段的核心子问题,如果t=0,执行步骤J;
J、动态规划正向递推阶段,令对t=1,2,...,T依次计算:
式中,qt *为电池储能***第t时段的最优充放电电量,为第t时段的电池储能***的最优存储电量,Lt为已知的第t时段电池储能***的固定损失电量;求解(17)式时受(3)-(6)式约束,根据动态规划原理,求解(3)-(6)和(17)式定义的极大化问题获得的从第1时段到第T时段的电池储能***的充放电量及对应的充放电功率就是电池储能***调度问题的最优解;
J、随时间发展进入下一时段,更新步骤B中电池储能***的初始存储电量、调度结束时的目标存储电量、各时段上网电量约束以及分时预测电价或者等效价值乘子,继续执行步骤C。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910778785.8A CN110504701B (zh) | 2019-08-22 | 2019-08-22 | 一种考虑动态充放电效率的电池储能***调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910778785.8A CN110504701B (zh) | 2019-08-22 | 2019-08-22 | 一种考虑动态充放电效率的电池储能***调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110504701A true CN110504701A (zh) | 2019-11-26 |
CN110504701B CN110504701B (zh) | 2023-01-20 |
Family
ID=68588764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910778785.8A Active CN110504701B (zh) | 2019-08-22 | 2019-08-22 | 一种考虑动态充放电效率的电池储能***调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110504701B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110880777A (zh) * | 2019-12-20 | 2020-03-13 | 国电南瑞科技股份有限公司 | 一种电网侧储能应急调峰备用能力的评估方法和装置 |
CN114565164A (zh) * | 2022-03-03 | 2022-05-31 | 国网江苏综合能源服务有限公司 | 一种储能日前调度混合整数规划的三阶段优化方法及*** |
CN116388334A (zh) * | 2023-04-14 | 2023-07-04 | 云储新能源科技有限公司 | 一种可重构电池储能*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968098B (zh) * | 2012-11-05 | 2014-12-10 | 清华大学 | 一种对集群内电动汽车充电功率的分布式优化方法 |
CN106684898A (zh) * | 2016-10-31 | 2017-05-17 | 国电南瑞科技股份有限公司 | 一种基于价值网络的储能***调度优化方法 |
CN106684916A (zh) * | 2017-02-16 | 2017-05-17 | 上海电力学院 | 一种带蓄电池的并网光伏***运行优化方法 |
CN107482679A (zh) * | 2017-08-25 | 2017-12-15 | 河海大学 | 考虑储能***充放电次数的主动配电网日前优化调度方法 |
US20180358812A1 (en) * | 2017-05-05 | 2018-12-13 | The Board Of Trustees Of The Leland Stanford Junior University | Systems and Methods for Coordinating Distributed Energy Storage |
-
2019
- 2019-08-22 CN CN201910778785.8A patent/CN110504701B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968098B (zh) * | 2012-11-05 | 2014-12-10 | 清华大学 | 一种对集群内电动汽车充电功率的分布式优化方法 |
CN106684898A (zh) * | 2016-10-31 | 2017-05-17 | 国电南瑞科技股份有限公司 | 一种基于价值网络的储能***调度优化方法 |
CN106684916A (zh) * | 2017-02-16 | 2017-05-17 | 上海电力学院 | 一种带蓄电池的并网光伏***运行优化方法 |
US20180358812A1 (en) * | 2017-05-05 | 2018-12-13 | The Board Of Trustees Of The Leland Stanford Junior University | Systems and Methods for Coordinating Distributed Energy Storage |
CN107482679A (zh) * | 2017-08-25 | 2017-12-15 | 河海大学 | 考虑储能***充放电次数的主动配电网日前优化调度方法 |
Non-Patent Citations (2)
Title |
---|
JIE YANG: "Optimal Dispatching Strategy for Shared Battery Station of Electric Vehicle by Divisional Battery Contro", 《IEEE ACCESS》 * |
文乐: "并网型多能互补***优化规划方法", 《热力发电》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110880777A (zh) * | 2019-12-20 | 2020-03-13 | 国电南瑞科技股份有限公司 | 一种电网侧储能应急调峰备用能力的评估方法和装置 |
CN114565164A (zh) * | 2022-03-03 | 2022-05-31 | 国网江苏综合能源服务有限公司 | 一种储能日前调度混合整数规划的三阶段优化方法及*** |
CN114565164B (zh) * | 2022-03-03 | 2024-02-06 | 国网江苏综合能源服务有限公司 | 一种储能日前调度混合整数规划的三阶段优化方法及*** |
CN116388334A (zh) * | 2023-04-14 | 2023-07-04 | 云储新能源科技有限公司 | 一种可重构电池储能*** |
CN116388334B (zh) * | 2023-04-14 | 2024-06-07 | 云储新能源科技有限公司 | 一种可重构电池储能*** |
Also Published As
Publication number | Publication date |
---|---|
CN110504701B (zh) | 2023-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113176511B (zh) | 考虑健康状态的储能充放电优化方法及*** | |
CN110504701A (zh) | 一种考虑动态充放电效率的电池储能***调度方法 | |
CN107565607B (zh) | 一种基于实时电价机制的微电网多时间尺度能量调度方法 | |
CN105846423A (zh) | 一种计及需求响应的光伏微电网储能多目标容量配置方法 | |
CN107394802A (zh) | 分布式储能参与自动发电控制的协调控制方法 | |
CN108710992A (zh) | 一种需求响应参与电网削峰填谷效果的评估方法 | |
CN109217290A (zh) | 计及电动汽车充放电的微网能量优化管理方法 | |
CN112101637B (zh) | 基于分时电价和能耗控制的动态充电优化方法及其应用 | |
CN110176765A (zh) | 一种峰谷电价驱动的储能调峰日前优化调度方法 | |
CN111064214A (zh) | 基于电动汽车两阶段滚动式策略的配电网优化调度方法 | |
CN103336998A (zh) | 一种风电场功率波动平抑目标值的优化计算方法 | |
CN114899856A (zh) | 一种电动汽车充电桩功率调节方法、***、设备及介质 | |
CN116470543A (zh) | 一种虚拟电厂的运行控制方法、装置、设备及介质 | |
CN109768560A (zh) | 提供多服务的蓄电池的充放电控制装置及方法 | |
CN110649638B (zh) | 一种弥补风电功率预测误差的储能***的优化方法 | |
CN117254464A (zh) | 一种储能***的控制方法及*** | |
CN105659464B (zh) | 蓄电池控制装置 | |
CN113437757B (zh) | 基于前景理论的风-蓄联合***电量分解方法 | |
CN110110929A (zh) | 一种充换电站的资源配置及充换电调度联合优化方法 | |
CN107332262B (zh) | 一种多类型混合储能的能量优化管理方法 | |
CN114330999A (zh) | 一种基于云储能***的虚拟电网资源匹配方法 | |
CN114194060A (zh) | 确定充电站的方法及装置 | |
CN109598856B (zh) | 一种储能计费方法及装置 | |
CN108599248B (zh) | 基于调控措施代价性能比指标的开机计划迭代寻优方法 | |
CN113320413A (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 |