一种微电网储能经济运行控制方法和装置与流程

本发明涉及微电网技术领域，具体涉及一种微电网储能经济运行控制方法和装置。

背景技术：

随着经济的快速发展和能耗日益增加，可再生能源的开发及可持续能源系统的构建成为电力工业的必然趋势，由于微电网具有能源利用率高和对环境友好等优点，所以微电网受到广泛的关注。微电网作为一种新的电力网络系统，其由分布式电源、储能单元和用电负荷等组成。储能单元作为微电网重要的组成部分，对储能单元的智能调度关系着微电网的大力发展。

储能单元在微电网中的作用主要体现在：(1)削峰填谷；(2)提高电能质量；(3)提高微电网经济性；(4)稳定可再生能源输出。综上所述，如何配置合理的各微电源及储能单元的功率出力，以达到满足调度要求的系统输出特性和经济性，是微电网智能调度的主要内容。目前对微电网中储能单元调度通常采用以下两种方法：方法1)人为制定策略，如分布式电源的出力大于负荷时，储存多余的电能，分布式电源出力不足以满足负荷时，储能单元则开始放电；方法2)利用智能算法或者经典算法计算储能单元的充放电功率。上述方法1)未考虑各时间段的关联，没有进行统筹安排调度，灵活性低且通常不能达到最优目标，上述方法2)在含有多个局部最优解问题时无法保证收敛到全局最优解，参数选取困难。

技术实现要素：

为了克服上述现有技术的不足，本发明提供一种微电网储能经济运行控制方法和装置，先确定储能单元电量的有效状态，然后根据储能单元电量的有效状态计算储能单元的充放电功率，并根据储能单元的充放电功率确定阶段指标函数，最后根据阶段指标函数确定储能单元的最优运行成本，并根据储能单元的最优运行成本最终确定储能单元充放电功率的最优调度策略，实现对微电网中储能单元的智能调度，确保储能单元的经济稳定运行。

为了实现上述发明目的，本发明采取如下技术方案：

本发明提供一种微电网储能经济运行控制方法，包括：

确定储能单元电量的有效状态；

根据储能单元电量的有效状态计算储能单元的充放电功率，并根据储能单元的充放电功率确定阶段指标函数；

根据阶段指标函数确定储能单元的最优运行成本，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略。

所述确定储能单元电量的有效状态包括：

判断i阶段储能单元的电量s(i)和n-i阶段储能单元的电量s(n-i)是否各自满足若是，s(i)和s(n-i)为有效状态，否则s(i)和s(n-i)为无效状态；

其中，s(i)表示i阶段储能单元的电量，s(n-i)表示n-i阶段储能单元的电量，s(i-1)min和s(i-1)max分别表示i-1阶段储能单元的电量的最小值和最大值，s(n-i+1)min和s(n-i+1)max分别表示n-i+1阶段储能单元的电量的最小值和最大值，pbcmax表示储能单元的最大充电功率，pbdmax表示储能单元的最大放电功率，ηd表示储能单元的放电效率，ηc表示储能单元的充电效率，δt表示相邻阶段的时间间隔。

所述根据储能单元电量的有效状态计算储能单元的充放电功率包括：

根据下式计算i阶段储能单元的充放电功率：

其中，pb(i)表示i阶段储能单元的充放电功率；sup(i)表示i阶段储能单元电量变化量，且sup(i)＝(1-τ)s(i-1)-s(i)，s(i-1)表示i-1阶段储能单元的电量，τ表示储能单元的自放电率；

根据下式计算n-i+1阶段储能单元的充放电功率：

其中，pb(n-i+1)表示n-i+1阶段储能单元的充放电功率；sup(n-i+1)表示n-i+1阶段储能单元电量变化量，sup(n-i+1)＝(1-τ)s(n-i)-s(n-i+1)，s(n-i)表示n-i阶段储能单元的电量。

所述根据储能单元的充放电功率确定阶段指标函数包括：

根据pb(i)确定如下式的i阶段储能单元的指标函数：

v(i)＝fb(i)+fgrid(i)

其中，v(i)表示i阶段指标函数，fb(i)表示i阶段储能单元的运行成本，fgrid(i)表示i阶段电网购售电成本，fb(i)按下式计算：

fb(i)＝μb×|pb(i)|

其中，μb表示储能单元单位充放电成本；

fgrid(i)按下式计算：

fgrid(i)＝μgrid×[pl(i)-pb(i)-σpg(i)]

其中，μgrid表示电网单位购售电成本，pl(i)表示i阶段负荷消耗的电量，∑pg(i)表示i阶段分布式电源的总发电量；

根据pb(n-i+1)确定如下式的n-i+1阶段指标函数：

v(n-i+1)＝fb(n-i+1)+fgrid(n-i+1)

其中，v(n-i+1)表示n-i+1阶段指标函数，fb(n-i+1)表示n-i+1阶段储能单元的运行成本，fgrid(n-i+1)表示n-i+1阶段电网购售电成本，fb(n-i+1)按下式计算：

fb(n-i+1)＝μb×|pb(n-i+1)|

fgrid(n-i+1)按下式计算：

fgrid(n-i+1)＝μgrid×[pl(n-i+1)-pb(n-i+1)-σpg(n-i+1)]

其中，pl(n-i+1)表示n-i+1阶段负荷消耗的电量，σpg(n-i+1)为n-i+1阶段分布式电源的总发电量。

所述根据阶段指标函数确定储能单元的最优运行成本包括：

根据v(i)确定如下式的i阶段储能单元的最优运行成本：

f(s(i,m))＝min{l∈e(i,m)|f(s(i-1,l))+v(i)}

根据v(n-i+1)确定如下式的n-i+1阶段储能单元的最优运行成本：

f(s(n-i,k))＝min{u∈e(n-i,k)|f(s(n-i+1,u))+v(n-i+1)}

其中，f(s(i,m))表示i阶段中有效状态m下储能单元的最优运行成本，e(i,m)表示i阶段中有效状态的集合，f(s(i-1,l))表示i-1阶段中有效状态l下储能单元的最优运行成本；f(s(n-i,k))表示n-i阶段中有效状态k下储能单元的最优运行成本，e(n-i,k)表示n-i阶段中有效状态的集合，f(s(n-i+1,u))表示n-i+1阶段中有效状态u下储能单元的最优运行成本。

所述根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略包括：

判断是否满足i＝n-i，若是，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略；否则使i＝i+1，返回确定储能单元电量的有效状态，并执行后续操作。

本发明还提供一种微电网储能经济运行控制装置，包括：

第一确定模块，用于确定储能单元电量的有效状态；

计算模块，用于根据储能单元电量的有效状态计算储能单元的充放电功率，并根据储能单元的充放电功率确定阶段指标函数；

第二确定模块，根据阶段指标函数确定储能单元的最优运行成本，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略。

所述第一确定单元具体用于：

判断i阶段储能单元的电量s(i)和n-i阶段储能单元的电量s(n-i)是否各自满足若是，s(i)和s(n-i)为有效状态，否则s(i)和s(n-i)为无效状态；

其中，s(i)表示i阶段储能单元的电量，s(n-i)表示n-i阶段储能单元的电量，s(i-1)min和s(i-1)max分别表示i-1阶段储能单元的电量的最小值和最大值，s(n-i+1)min和s(n-i+1)max分别表示n-i+1阶段储能单元的电量的最小值和最大值，pbcmax表示储能单元的最大充电功率，pbdmax表示储能单元的最大放电功率，ηd表示储能单元的放电效率，ηc表示储能单元的充电效率，δt表示相邻阶段的时间间隔。

所述计算模块具体用于：

根据下式计算i阶段储能单元的充放电功率：

其中，pb(i)表示i阶段储能单元的充放电功率；sup(i)表示i阶段储能单元电量变化量，且sup(i)＝(1-τ)s(i-1)-s(i)，s(i-1)表示i-1阶段储能单元的电量，τ表示储能单元的自放电率；

根据下式计算n-i+1阶段储能单元的充放电功率：

其中，pb(n-i+1)表示n-i+1阶段储能单元的充放电功率；sup(n-i+1)表示n-i+1阶段储能单元电量变化量，sup(n-i+1)＝(1-τ)s(n-i)-s(n-i+1)，s(n-i)表示n-i阶段储能单元的电量。

所述计算模块还具体用于：

根据pb(i)计算如下式的i阶段储能单元的指标函数：

v(i)＝fb(i)+fgrid(i)

其中，v(i)表示i阶段指标函数，fb(i)表示i阶段储能单元的运行成本，fgrid(i)表示i阶段电网购售电成本，fb(i)按下式计算：

fb(i)＝μb×|pb(i)|

其中，μb表示储能单元单位充放电成本；

fgrid(i)按下式计算：

fgrid(i)＝μgrid×[pl(i)-pb(i)-σpg(i)]

其中，μgrid表示电网单位购售电成本，pl(i)表示i阶段负荷消耗的电量，∑pg(i)表示i阶段分布式电源的总发电量；

根据pb(n-i+1)确定如下式的n-i+1阶段指标函数：

v(n-i+1)＝fb(n-i+1)+fgrid(n-i+1)

其中，v(n-i+1)表示n-i+1阶段指标函数，fb(n-i+1)表示n-i+1阶段储能单元的运行成本，fgrid(n-i+1)表示n-i+1阶段电网购售电成本，fb(n-i+1)按下式计算：

fb(n-i+1)＝μb×|pb(n-i+1)|

fgrid(n-i+1)按下式计算：

fgrid(n-i+1)＝μgrid×[pl(n-i+1)-pb(n-i+1)-σpg(n-i+1)]

其中，pl(n-i+1)表示n-i+1阶段负荷消耗的电量，σpg(n-i+1)为n-i+1阶段分布式电源的总发电量。

所述第二确定模块具体用于：

根据v(i)确定如下式的i阶段储能单元的最优运行成本：

f(s(i,m))＝min{l∈e(i,m)|f(s(i-1,l))+v(i)}

根据v(n-i+1)确定如下式的n-i+1阶段储能单元的最优运行成本：

f(s(n-i,k))＝min{u∈e(n-i,k)|f(s(n-i+1,u))+v(n-i+1)}

其中，f(s(i,m))表示i阶段中有效状态m下储能单元的最优运行成本，e(i,m)表示i阶段中有效状态的集合，f(s(i-1,l))表示i-1阶段中有效状态l下储能单元的最优运行成本；f(s(n-i,k))表示n-i阶段中有效状态k下储能单元的最优运行成本，e(n-i,k)表示n-i阶段中有效状态的集合，f(s(n-i+1,u))表示n-i+1阶段中有效状态u下储能单元的最优运行成本。

所述第二确定模块还具体用于：

判断是否满足i＝n-i，若是，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略；否则使i＝i+1，返回确定储能单元电量的有效状态，并执行后续操作。

与最接近的现有技术相比，本发明提供的技术方案具有以下有益效果：

本发明提供了一种微电网储能经济运行控制方法和装置，先确定储能单元电量的有效状态，然后根据储能单元电量的有效状态计算储能单元的充放电功率，并根据储能单元的充放电功率确定阶段指标函数，最后根据阶段指标函数确定储能单元的最优运行成本，并根据储能单元的最优运行成本最终确定储能单元充放电功率的最优调度策略，实现对微电网中储能单元的智能调度，确保储能单元的经济稳定运行；

本发明提供的技术方案不仅考虑了各时间阶段的关联，能够进行统筹安排调度，灵活性较高且能够达到经济最优的控制目标，而且本发明提供的技术方案在含有多个局部最优解问题时能够保证收敛到全局最优解，参数选取较为简单；

本发明提供的技术方案在确定储能单元电量的有效状态的过程中考虑了储能单元的充放电功率和电量，且在计算储能单元电量变化量的过程中考虑了储能单元的自放电率；

本发明提供的技术方案基于双向动态规划算法从顺序和逆序两个方向向中间进行搜索，通过改进各个阶段的状态数，大大减少了动态规划算法的计算量，提高了时间效率，且当计算时间要求不高时，该双向动态规划算法可以给出更符合优化目标的储能单元经济运行调度策略。

附图说明

图1是本发明实施例中双向动态规划算法示意图；

图2是本发明实施例中微电网储能经济运行控制方法流程图。

具体实施方式

下面结合附图对本发明作进一步详细说明。

本发明实施例提供的微电网储能经济运行控制方法基于动态规划算法确定储能单元充放电功率的最优调度策略，实现对微电网中储能单元的智能调度，确保储能单元的经济稳定运行。

动态规划算法的时间复杂度主要以下三个因素决定：

因素1)：问题中状态的总个数，

因素2)：每次状态转移所涉及的状态数，

因素3)：每次状态转移的时间。

若对上述三个因素中的一个或多个进行优化，就可以提高动态规划算法的时间效率。

针对微电网中储能单元的特性，本发明实施例根据储能单元的最大放电功率和储能单元的最大充电功率，舍掉每个阶段中无效状态，减少了状态总数即可降低了动态规划算法的空间复杂度，此时可从整体上提高算法的时间效率。

另外，由于储能单元在一个调度周期内的初末电量状态要求一致以满足循环使用，即初末状态都已经确定，故可采用双向动态规划算法即从顺序和逆序两个方向向中间进行搜索(如图1分别从阶段0和阶段n为起点向中间阶段进行搜索)，搜索沿着两个方向同时进行，一是正向搜索即初始状态向目标状态进行搜索，另一个是逆向搜索即从目标状态向初始状态进行搜索，当两个方向的搜索汇合时终止搜索并进行判断，得到问题的答案，这将大大减少计算量，提高时间效率。

本发明实施例提供的微电网储能经济运行控制方法具体流程图如图2所示，其具体过程如下：

s101：确定储能单元电量的有效状态；

s102：根据s101确定的储能单元电量的有效状态计算储能单元的充放电功率，并根据储能单元的充放电功率确定阶段指标函数；

s103：根据s102确定的阶段指标函数确定储能单元的最优运行成本，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略。

上述s101中，确定储能单元电量的有效状态具体过程如下：

判断i阶段储能单元的电量s(i)和n-i阶段储能单元的电量s(n-i)是否各自满足若是，s(i)和s(n-i)为有效状态，否则s(i)和s(n-i)为无效状态；

其中，s(i)表示i阶段储能单元的电量，s(n-i)表示n-i阶段储能单元的电量，s(i-1)min和s(i-1)max分别表示i-1阶段储能单元的电量的最小值和最大值，s(n-i+1)min和s(n-i+1)max分别表示n-i+1阶段储能单元的电量的最小值和最大值，pbcmax表示储能单元的最大充电功率，pbdmax表示储能单元的最大放电功率，ηd表示储能单元的放电效率，ηc表示储能单元的充电效率，δt表示相邻阶段的时间间隔。

上述s102中，根据储能单元电量的有效状态计算储能单元的充放电功率主要过程如下：

1)根据下式计算i阶段储能单元的充放电功率：

其中，pb(i)表示i阶段储能单元的充放电功率；sup(i)表示i阶段储能单元电量变化量，且sup(i)＝(1-τ)s(i-1)-s(i)，s(i-1)表示i-1阶段储能单元的电量，τ表示储能单元的自放电率；

2)根据下式计算n-i+1阶段储能单元的充放电功率：

其中，pb(n-i+1)表示n-i+1阶段储能单元的充放电功率；sup(n-i+1)表示n-i+1阶段储能单元电量变化量，sup(n-i+1)＝(1-τ)s(n-i)-s(n-i+1)，s(n-i)表示n-i阶段储能单元的电量。

上述s102中进一步根据储能单元的充放电功率确定阶段指标函数具体过程如下：

1)根据pb(i)确定如下式的i阶段储能单元的指标函数：

v(i)＝fb(i)+fgrid(i)

其中，v(i)表示i阶段指标函数，fb(i)表示i阶段储能单元的运行成本，fgrid(i)表示i阶段电网购售电成本，fb(i)按下式计算：

fb(i)＝μb×|pb(i)|

其中，μb表示储能单元单位充放电成本；

fgrid(i)按下式计算：

fgrid(i)＝μgrid×[pl(i)-pb(i)-∑pg(i)]

其中，μgrid表示电网单位购售电成本，pl(i)表示i阶段负荷消耗的电量，σpg(i)表示i阶段分布式电源的总发电量；

2)根据pb(n-i+1)确定如下式的n-i+1阶段指标函数：

v(n-i+1)＝fb(n-i+1)+fgrid(n-i+1)

其中，v(n-i+1)表示n-i+1阶段指标函数，fb(n-i+1)表示n-i+1阶段储能单元的运行成本，fgrid(n-i+1)表示n-i+1阶段电网购售电成本，fb(n-i+1)按下式计算：

fb(n-i+1)＝μb×|pb(n-i+1)|

fgrid(n-i+1)按下式计算：

fgrid(n-i+1)＝μgrid×[pl(n-i+1)-pb(n-i+1)-σpg(n-i+1)]

其中，pl(n-i+1)表示n-i+1阶段负荷消耗的电量，σpg(n-i+1)为n-i+1阶段分布式电源的总发电量。

上述s103中，根据阶段指标函数确定储能单元的最优运行成本具体过程如下：

1)根据v(i)确定如下式的i阶段储能单元的最优运行成本：

f(s(i,m))＝min{l∈e(i,m)|f(s(i-1,l))+v(i)}

2)根据v(n-i+1)确定如下式的n-i+1阶段储能单元的最优运行成本：

f(s(n-i,k))＝min{u∈e(n-i,k)|f(s(n-i+1,u))+v(n-i+1)}

其中，f(s(i,m))表示i阶段中有效状态m下储能单元的最优运行成本，e(i,m)表示i阶段中有效状态的集合，f(s(i-1,l))表示i-1阶段中有效状态l下储能单元的最优运行成本；f(s(n-i,k))表示n-i阶段中有效状态k下储能单元的最优运行成本，e(n-i,k)表示n-i阶段中有效状态的集合，f(s(n-i+1,u))表示n-i+1阶段中有效状态u下储能单元的最优运行成本。

上述s103中进一步根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略具体过程如下：

判断是否满足i＝n-i，若是，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略；否则使i＝i+1，返回确定储能单元电量的有效状态，并执行后续操作。

实施例

假设图1中阶段0-阶段n对应为一天的24小时，即调度周期为24小时，储能单元从电量最满到电量为0(这里以电量多少为参数区分不同的状态)的这一过程包括有0～n个状态，图1中黑色圆圈代表储能单元为有效状态、白色圆圈代表储能单元为无效状态，下面对本发明实施例的方案进行详细描述：

如图1所示，取区间为1小时，则阶段分别为0，1，2，……，24，图1中的m表示储能单元电量的初末状态，假设阶段0和阶段24的电量状态为20kwh，即m取20kwh，储能单元的最大容量为50kwh，最大充放电功率为10kw，充放电效率为1，若把电池电量分为6份，即n＝5，则阶段i(i不等于0和24)有6个状态(即状态电量)，6个状态分别为0，10，20，30，40，50，因此阶段1的有效状态为10，20，30，同样此时阶段23的有效状态为10，20，30，依次计算得到所有阶段的有效状态。

下面确定各阶段的调度策略，如：

阶段1下有三个有效状态需要依次计算，当有效状态为10kwh时，对应阶段0的状态20kwh，此时的调度策略为(10-20＝-10kw)放电功率为10kw，根据这个调度策略和储能单元的最优运行成本函数得到一个成本，因为阶段0只有一个状态所以结束计算，跳到阶段1下一个状态20kwh，对应阶段0的状态20kwh，得到的调度策略为(20-20＝0)即为浮充状态，充放电功率为0，由这个调度策略和储能单元的最优运行成本函数得到一个成本，再计算阶段1状态为30kwh的成本。

这样得到了阶段1三个状态下的成本值，跳到阶段2，根据之前的保留的有效状态，阶段2的状态应该是0，10，20，30，40依次遍历，当状态为0时，对应的阶段1的状态有0和10，依次计算这两个状态下的成本，取成本低的那个状态，同时观察此成本值下阶段1的状态；

相同计算方法也在阶段24开始向前，当到阶段12时，两个方向的搜索汇合，此时判断阶段12各个有效状态的最低成本，并得到他们各自此成本值下前后一阶段下的状态，逐阶段搜索，得到24小时的所有调度策略，即得到了最优调度策略。

基于同一发明构思，本发明实施例还提供了微电网储能经济运行控制装置，这些设备解决问题的原理与微电网储能经济运行控制方法相似，该微电网储能经济运行控制装置可以包括第一确定模块、计算模块和第二确定模块，下面分别介绍上述3个模块的具体功能：

其中的第一确定模块，用于确定储能单元电量的有效状态；

其中的计算模块，用于根据第一确定模块确定的储能单元电量的有效状态计算储能单元的充放电功率，并根据储能单元的充放电功率确定阶段指标函数；

其中的第二确定模块，用于根据计算模块确定的阶段指标函数确定储能单元的最优运行成本，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略。

上述的第一确定单元确定储能单元电量的有效状态具体过程如下：

判断i阶段储能单元的电量s(i)和n-i阶段储能单元的电量s(n-i)是否各自满足若是，s(i)和s(n-i)为有效状态，否则s(i)和s(n-i)为无效状态；

其中，s(i)表示i阶段储能单元的电量，s(n-i)表示n-i阶段储能单元的电量，s(i-1)min和s(i-1)max分别表示i-1阶段储能单元的电量的最小值和最大值，s(n-i+1)min和s(n-i+1)max分别表示n-i+1阶段储能单元的电量的最小值和最大值，pbcmax表示储能单元的最大充电功率，pbdmax表示储能单元的最大放电功率，ηd表示储能单元的放电效率，ηc表示储能单元的充电效率，δt表示相邻阶段的时间间隔。

上述的计算模块根据第一确定模块确定的储能单元电量的有效状态计算储能单元的充放电功率具体过程如下：

1)根据下式计算i阶段储能单元的充放电功率：

其中，pb(i)表示i阶段储能单元的充放电功率；sup(i)表示i阶段储能单元电量变化量，且sup(i)＝(1-τ)s(i-1)-s(i)，s(i-1)表示i-1阶段储能单元的电量，τ表示储能单元的自放电率；

2)根据下式计算n-i+1阶段储能单元的充放电功率：

其中，pb(n-i+1)表示n-i+1阶段储能单元的充放电功率；sup(n-i+1)表示n-i+1阶段储能单元电量变化量，sup(n-i+1)＝(1-τ)s(n-i)-s(n-i+1)，s(n-i)表示n-i阶段储能单元的电量。

上述的计算模块根据储能单元的充放电功率确定阶段指标函数具体过程如下：

1)根据pb(i)计算如下式的i阶段储能单元的指标函数：

v(i)＝fb(i)+fgrid(i)

其中，v(i)表示i阶段指标函数，fb(i)表示i阶段储能单元的运行成本，fgrid(i)表示i阶段电网购售电成本，fb(i)按下式计算：

fb(i)＝μb×|pb(i)|

其中，μb表示储能单元单位充放电成本；

fgrid(i)按下式计算：

fgrid(i)＝μgrid×[pl(i)-pb(i)-σpg(i)]

其中，μgrid表示电网单位购售电成本，pl(i)表示i阶段负荷消耗的电量，σpg(i)表示i阶段分布式电源的总发电量；

2)根据pb(n-i+1)确定如下式的n-i+1阶段指标函数：

v(n-i+1)＝fb(n-i+1)+fgrid(n-i+1)

其中，v(n-i+1)表示n-i+1阶段指标函数，fb(n-i+1)表示n-i+1阶段储能单元的运行成本，fgrid(n-i+1)表示n-i+1阶段电网购售电成本，fb(n-i+1)按下式计算：

fb(n-i+1)＝μb×|pb(n-i+1)|

fgrid(n-i+1)按下式计算：

fgrid(n-i+1)＝μgrid×[pl(n-i+1)-pb(n-i+1)-σpg(n-i+1)]

其中，pl(n-i+1)表示n-i+1阶段负荷消耗的电量，∑pg(n-i+1)为n-i+1阶段分布式电源的总发电量。

上述的第二确定模块根据计算模块确定的阶段指标函数确定储能单元的最优运行成本具体过程如下：

1)根据v(i)确定如下式的i阶段储能单元的最优运行成本：

f(s(i,m))＝min{l∈e(i,m)|f(s(i-1,l))+v(i)}

2)根据v(n-i+1)确定如下式的n-i+1阶段储能单元的最优运行成本：

f(s(n-i,k))＝min{u∈e(n-i,k)|f(s(n-i+1,u))+v(n-i+1)}

其中，f(s(i,m))表示i阶段中有效状态m下储能单元的最优运行成本，e(i,m)表示i阶段中有效状态的集合，f(s(i-1,l))表示i-1阶段中有效状态l下储能单元的最优运行成本；f(s(n-i,k))表示n-i阶段中有效状态k下储能单元的最优运行成本，e(n-i,k)表示n-i阶段中有效状态的集合，f(s(n-i+1,u))表示n-i+1阶段中有效状态u下储能单元的最优运行成本。

上述第二确定模块根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略具体过程如下：

判断是否满足i＝n-i，若是，根据储能单元的最优运行成本确定储能单元充放电功率的最优调度策略；否则使i＝i+1，返回确定储能单元电量的有效状态，并执行后续操作。

为了描述的方便，以上所述装置的各部分以功能分为各种模块或单元分别描述。当然，在实施本申请时可以把各模块或单元的功能在同一个或多个软件或硬件中实现。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。