电气工程学报, 2022, 17(3): 210-218 doi: 10.11985/2022.03.025

新能源发电与电能存储

考虑电池储能寿命损耗的风光储微网随机日前调度策略*

赵诣,1, 方陈,2, 冯冬涵,1, 李恒杰1,3, 时珊珊2, 周云,1

1.电力传输与功率变换控制教育部重点实验室(上海交通大学) 上海 200240

2.华东电力试验研究院有限公司 上海 200437

3.兰州理工大学电气工程与信息工程学院 兰州 730050

Stochastic Day-ahead Scheduling Optimization for Wind-PV-ES Microgrid Considering Battery Life Loss

ZHAO Yi,1, FANG Chen,2, FENG Donghan,1, LI Hengjie1,3, SHI Shanshan2, ZHOU Yun,1

1. Key Laboratory of Control of Power Transmission and Conversion, Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240

2. East China Electric Power Test and Research Institute Co., Ltd., Shanghai 200437

3. School of Electrical Engineering and Information Engineering, Lanzhou University of Technology, Lanzhou 730050

通讯作者: 周云,男,1990年生,博士,讲师。主要研究方向为电力/综合能源系统应急及用电能效管理等。E-mail:yun.zhou@sjtu.edu.cn

收稿日期: 2021-08-22   修回日期: 2022-03-24  

基金资助: *国家自然科学基金.  51907122
国家自然科学基金.  52077139
上海市科委计划.  21DZ2204800
国家电网有限公司科技.  52094021000F

Received: 2021-08-22   Revised: 2022-03-24  

作者简介 About authors

赵诣,男,1998年生,硕士研究生。主要研究方向为用电能效管理。E-mail: Yi.Zhao.sjtu@sjtu.edu.cn

方陈,男,1983年生,高级工程师。主要研究方向为能源互联网多能互补规划与运行,分布式能源、储能、电动汽车、微电网规划与运行等。E-mail: fangc02@gmail.com

冯冬涵,男,1981年生,博士,教授。主要研究方向为智能电网中的策略与风险及电力市场理论与设计等。E-mail: seed@sjtu.edu.cn

摘要

随着我国“碳达峰”和“碳中和”目标的提出,以新能源为主的电力系统已经受到了广泛的关注,风光储微网就是其中的典型场景。研究电池储能的寿命损耗特性,建立精细的分段线性化电池储能寿命损耗模型。将建立的电池储能寿命损耗模型加入到风光储微网场景中,考虑储能技术在提高风光资源消纳率以及减少碳排放量方面的作用,以最小化微网和上级电网交互费用、弃风弃光成本、削负荷成本、碳排放成本以及电池储能寿命损耗成本为目标函数,对整个风光储微网的随机日前调度模型进行线性化求解。对风光储微网的调度模型算例分析表明,所提模型可以通过规范电池储能的充放电动作,有效提高电池的使用寿命以及提高风光储微网的经济性。

关键词: 风光储微网 ; 风光消纳 ; 电池储能寿命损耗 ; 随机调度 ; 线性化

Abstract

With the ‘carbon peak’ and ‘carbon neutrality’ goals proposed, the power system based on renewable resources has received extensive attention, and the wind-PV-ES(Energy storage) microgrid is the typical scenario. The piecewise linear life loss model is established by studying the life loss characteristics of batteries. Combined with the role of energy storage technology in improving the wind and solar power consumption and reducing carbon emissions, the economic effect of battery on the microgrid is added to the wind-PV-ES microgrid scenarios. Regarding the interaction cost of the grid, wind and solar abandonment cost, shaving load cost, carbon emissions cost, and battery life loss cost as the objective function, the day-ahead scheduling model of entire wind-PV-ES microgrid is linearized. The case is used to analyze that the introduction of batteries can effectively increase the battery life and improve the economy of the microgrid by regulating the charging and discharging actions of batteries.

Keywords: Wind-PV-ES microgrid ; wind-PV accommodation ; battery life loss ; stochastic scheduling ; linearization

PDF (382KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

赵诣, 方陈, 冯冬涵, 李恒杰, 时珊珊, 周云. 考虑电池储能寿命损耗的风光储微网随机日前调度策略*. 电气工程学报[J], 2022, 17(3): 210-218 doi:10.11985/2022.03.025

ZHAO Yi, FANG Chen, FENG Donghan, LI Hengjie, SHI Shanshan, ZHOU Yun. Stochastic Day-ahead Scheduling Optimization for Wind-PV-ES Microgrid Considering Battery Life Loss. Chinese Journal of Electrical Engineering[J], 2022, 17(3): 210-218 doi:10.11985/2022.03.025

1 引言

随着如今全球经济的快速发展,与之而来的是碳排放量不断上升的问题,节能减排已经成为各国的共识,我国对应提出了“碳达峰”和“碳中和”目标[1]。为尽快实现该目标,构建以风、光等可再生能源为主的电力系统受到了广泛关注。风、光等可再生能源具有零污染、分布广、蕴含量大等优点,在节能减排方面有着很大优势,实践证明,我国光伏资源和风力资源具有天然的互补特性。但由于风力发电和光伏发电都受到自然因素的影响,其对应出力具有波动性和间歇性,这种不确定性阻碍了微网的经济发展,所以国内外针对包括风光储微网在内的多种微网形式的经济配置与调度运行进行了研究。在独立型微网方面,文献[2]以独立型微网发电成本、切负荷补偿成本及环境效益最优为目标,实现多场景下独立微网的经济调度;文献[3]同样考虑独立型微网的经济性和环保性,针对各电源容量配置对算法进行创新;在并网型微网方面,文献[4]建立了考虑风光出力波动的并网型微网的最优潮流模型;文献[5]综合微网的经济性、可靠性和可再生能源利用率设计了一种综合能量调度策略,研究了并网型风光储微网容量配置问题。

储能设备不但可以用来克服可再生能源出力不稳定的情况,而且对于碳排放也起着重要作用[6],所以电池储能在提升风光储微网经济性方面扮演着重要角色。电池储能的使用寿命和整个微网运行过程中的充放电行为有着很大的关系,如果为了微网运行的经济性,让电池频繁充放电,会使其寿命下降迅速,从长远来看并不利于风光储微网运行的经济性。电池储能寿命损耗程度一般与其充放电倍率、工作环境温度和放电深度(Depth of discharge,DOD)有关[7]。如今国内外对电池储能寿命损耗成本也已经有了不少的研究,包括考虑总吞吐量的电池寿命损耗建模方式[8],利用电池全生命周期总吞吐量来拟合度电成本与荷电状态(State of charge,SOC)函数关系;还有基于试验和疲劳分析的建模方法[9],通过试验得到电池储能寿命循环次数和放电深度的关系曲线,再结合雨流计数法求取整个放电过程的电池储能寿命损耗,但这种方法难以嵌入优化模型中,且没有考虑到每次充放电初始SOC不同。为了进一步衡量电池储能对风光储微网经济性的作用,目前电池储能的寿命损耗模型已经被加入到了风光储微网的研究中,大部分文献均以储能固定的运维成本为寿命损耗[10]或通过类似雨流计数法思路,利用充放电次数计算储能寿命损耗加入到目标函数中[11],并未考虑储能寿命损耗和放电深度的关系,不符合实际的储能寿命损耗。

本文针对并网型风光储微网系统进行研究,考虑风光出力不确定性,利用电池储能寿命损耗的精细化模型,建立了以最小化电网交互费用、弃风弃光成本、削负荷成本、碳排放成本、电池储能寿命损耗成本为目标函数的风光储微网随机日前调度模型,并对其线性化求解,通过算例分析,文中所提模型可以通过规范储能的充放电动作,有效提高电池的使用寿命以及风光储微网的经济性。

2 风光储微网结构

2.1 系统模型

风光储微网系统结构如图1所示。主要由分布式电源、储能单元以及微网负荷这三部分组成,分布式电源包括光伏及风力装置,通过电力电子元件与交流母线相连,整个微网系统与主网并网运行。

图1

图1   风光储微网系统结构图


2.2 发电单元模型

2.2.1 风力发电模型

风力发电功率波动性比较大,其具体出力和实际风速相关,一般考虑出力表达式为[12]

${{P}_{\text{W}}}\text{=}\left\{ \begin{align} & \text{0}\ v<{{v}_{\text{c}}} \\ & {{P}_{\text{R}}}\cdot \frac{v-{{v}_{\text{c}}}}{{{v}_{\text{R}}}-{{v}_{\text{c}}}}\ {{v}_{\text{c}}}\le v\le {{v}_{\text{R}}} \\ & {{P}_{\text{R}}}\ {{v}_{\text{R}}}\le v\le {{v}_{\text{F}}} \\ & \text{0}\ v>{{v}_{\text{F}}} \\ \end{align} \right.$

式中,${{P}_{\text{W}}}$为风电出力;${{P}_{\text{R}}}$为风电额定功率;${{v}_{\text{c}}}$为切断风速;${{v}_{\text{R}}}$为额定风速;${{v}_{\text{F}}}$为截断风速。

2.2.2 光伏发电模型

光伏发电一般认为受到温度和光照强度的影响,具体出力表达式为[12]

${{P}_{\text{PV}}}={{P}_{\text{STC}}}\cdot {{G}_{\text{AC}}}\cdot \frac{1+k\cdot ({{T}_{\text{c}}}-{{T}_{\text{r}}})}{{{G}_{\text{STC}}}}$

式中,${{P}_{\text{PV}}}$为光伏出力;${{G}_{\text{AC}}}$为光照强度;${{P}_{\text{STC}}}$为标准测试条件下(太阳光入射强度1 000 $\text{W}/{{\text{m}}^{2}}$,环境温度为25 ℃)的最大测试功率;${{G}_{\text{STC}}}$为标准测试条件下的光照强度,值为1 000 $\text{W}/{{\text{m}}^{2}}$;$k$为功率温度系数;${{T}_{\text{c}}}$为电池板工作温度;${{T}_{\text{r}}}$为参考温度。

2.2.3 电池储能运行模型

电池储能在运行时主要关注的是其各个时刻的SOC。

充电时刻满足

${{S}_{\text{OC}}}(t)={{S}_{\text{OC}}}(t-1)+\frac{{{P}_{\text{c}}}\cdot {{\eta }_{\text{c}}}}{{{E}_{\text{rated}}}}\cdot \Delta t$

放电时刻满足

${{S}_{\text{OC}}}(t)={{S}_{\text{OC}}}(t-1)-\frac{{{P}_{\text{d}}}}{{{\eta }_{\text{d}}}\cdot {{E}_{\text{rated}}}}\cdot \Delta t$

式中,${{S}_{\text{OC}}}(t-1)$${{S}_{\text{OC}}}(t)$分别为$t-1$时刻和$t$时刻电池对应的荷电状态;$\Delta t$为离散时间间隔;${{P}_{\text{c}}}$${{P}_{\text{d}}}$分别为电池的充电和放电功率;${{\eta }_{\text{c}}}$${{\eta }_{\text{d}}}$分别为电池的充电和放电效率;${{E}_{\text{rated}}}$为电池的额定容量。

3 电池储能寿命损耗模型

3.1 电池储能寿命损耗原模型

电池储能的寿命和损耗与它的荷电状态和充放电功率密不可分。本文使用了一种精细化的电池储能寿命损耗模型,该模型建立在试验得到的DOD-循环寿命曲线上,该曲线是通过在不同DOD下,进行不断充放电试验得出的,比如DOD=80%,从SOC为100%放电至20%,再从20%充电至100%,记为一次充放电循环,重复该过程直至电池容量衰减至额定容量的80%,记录总计循环次数。一般使用自然指数拟合该函数曲线

${{N}_{\text{life}}}({{D}_{\text{OD}}})={{\mu }_{1}}\cdot \exp \left( -{{\mu }_{2}}\cdot {{D}_{\text{OD}}} \right)$

式中,${{N}_{\text{life}}}$表示拟合电池循环寿命与放电深度关系的自然指数函数;${{\mu }_{1}}$${{\mu }_{\text{2}}}$为拟合函数对应的系数。

假设电池每一次相同放电深度下对应的寿命损耗相同且充放电过程寿命损耗对称,定义电池SOC状态由${{S}_{\text{OC}}}$充电至100%的寿命损耗${{L}_{\text{loss}}}({{S}_{\text{OC}}},1)$

${{L}_{\text{loss}}}({{S}_{\text{OC}}},1)=\frac{1}{2{{N}_{\text{life}}}(1-{{S}_{\text{OC}}})}$

假设电池充放电状态在相邻时刻两个SOC状态之间保持不变,定义电池储能寿命损耗系数${{\lambda }_{\text{loss}}}({{S}_{\text{OC}}})$

$\begin{matrix} {{\lambda }_{\text{loss}}}({{S}_{\text{OC}}})=\underset{\Delta {{S}_{\text{OC}}}\to 0}{\mathop{\lim }}\,\frac{{{L}_{\text{loss}}}({{S}_{\text{OC}}},{{S}_{\text{OC}}}+\Delta {{S}_{\text{OC}}})}{\Delta {{S}_{\text{OC}}}\cdot {{E}_{\text{rated}}}}= \\ \frac{1}{{{E}_{\text{rated}}}}\cdot \frac{d{{L}_{\text{loss}}}({{S}_{\text{OC}}},1)}{d{{S}_{\text{OC}}}}=-\frac{1}{2{{E}_{\text{rated}}}}\cdot {{\left( \frac{1}{{{N}_{\text{life}}}(1-{{S}_{\text{OC}}})} \right)}^{\prime }} \\ \end{matrix}$

根据寿命损耗系数,电池荷电状态由${{S}_{\text{OC}}}(t-1)$变化至${{S}_{\text{OC}}}(t)$对应的寿命损耗可表示为

$\begin{align} & {{L}_{\text{loss}}}(t)=\left| \int_{{{S}_{\text{OC}}}(t-1)}^{{{S}_{\text{OC}}}(t)}{{{\lambda }_{\text{loss}}}({{S}_{\text{OC}}})d({{S}_{\text{OC}}}{{E}_{\text{rated}}})} \right|= \\ & \text{ }\left| F({{S}_{\text{OC}}}(t))-F({{S}_{\text{OC}}}(t-1)) \right| \\ \end{align}$

其中,$F({{S}_{\text{OC}}})$是电池储能寿命损耗系数的原函数,可表示为

$\begin{align} & F({{S}_{\text{OC}}})=\int_{0}^{{{S}_{\text{OC}}}}{{{\lambda }_{\text{loss}}}({{S}_{\text{OC}}})\cdot {{E}_{\text{rated}}}d{{S}_{\text{OC}}}}= \\ & \text{ }\frac{1}{2}\left( \frac{1}{{{N}_{\text{life}}}(1)}-\frac{1}{{{N}_{\text{life}}}(1-{{S}_{\text{OC}}})} \right)\text{ } \\ \end{align}$

3.2 电池储能寿命损耗分段线性化模型

通过精细化的电池储能寿命损耗模型,可以得到电池任意时刻SOC状态变化后对应的寿命损耗。但由于电池储能的DOD-循环寿命曲线较为复杂,难以写出显式表达式嵌入优化模型分析,所以利用分段线性化的思想得到线性化后的电池储能寿命损耗模型为

$F({{S}_{\text{OC}}})\approx F(S_{\text{OC}}^{\text{min}})+\sum\limits_{\lambda =1}^{\Lambda }{\varphi _{{{S}_{\text{OC}}}}^{\lambda }\cdot }\Delta _{{{S}_{\text{OC}}}}^{\lambda }$
$\sum\limits_{\lambda =1}^{N}{\Delta _{{{S}_{\text{OC}}}}^{\lambda }}={{S}_{\text{OC}}}-S_{\text{OC}}^{\text{min}}$
${{\bar{\Delta }}_{{{S}_{\text{OC}}}}}=(S_{\text{OC}}^{\text{max}}-S_{\text{OC}}^{\text{min}})/N$
$\varphi _{{{S}_{\text{OC}}}}^{\lambda }=\frac{F(S_{\text{OC}}^{\text{min}}+\lambda \cdot {{{\bar{\Delta }}}_{{{S}_{\text{OC}}}}})-F(S_{\text{OC}}^{\text{min}}+(\lambda -1)\cdot {{{\bar{\Delta }}}_{{{S}_{\text{OC}}}}})}{{{{\bar{\Delta }}}_{{{S}_{\text{OC}}}}}}$
$\Delta _{{{S}_{\text{OC}}}}^{\lambda }\text{=}\max \left\{ \min \{{{S}_{\text{OC}}}-S_{\text{OC}}^{\text{min}}-(\lambda -1)\cdot {{{\bar{\Delta }}}_{{{S}_{\text{OC}}}}},{{{\bar{\Delta }}}_{{{S}_{\text{OC}}}}}\},0 \right\}$

式中,N为分段数;$\lambda $为分段的序号;$\Delta _{{{S}_{\text{OC}}}}^{\lambda }$为对应分段的长度;${{\bar{\Delta }}_{{{S}_{\text{OC}}}}}$为分段的长度上限;$\varphi _{{{S}_{\text{OC}}}}^{\lambda }$为对应分段的斜率。通过式(14)在保证后一段取非零值的情况下,前面每一段长度都取为最大值,从而得到可以嵌入优化模型的电池储能寿命损耗模型。

为分析分段线性化后电池储能寿命损耗模型精确度,利用其和原函数之间的相对误差进行衡量,其具体含义如图2所示。

图2

图2   分段线性化后相对误差图


$\beta ({{S}_{\text{OC}}})=\left| F({{S}_{\text{OC}}})-F(S_{\text{OC}}^{\text{min}})-\sum\limits_{\lambda =1}^{\Lambda }{\varphi _{{{S}_{\text{OC}}}}^{\lambda }\cdot }\Delta _{{{S}_{\text{OC}}}}^{\lambda } \right|$
$\delta ({{S}_{\text{OC}}})=\frac{\beta ({{S}_{\text{OC}}})}{F({{S}_{\text{OC}}})}$
${{\delta }^{\text{avg}}}=\frac{\int_{S_{\text{OC}}^{\text{min}}}^{S_{\text{OC}}^{\text{max}}}{\beta ({{S}_{\text{OC}}})d({{S}_{\text{OC}}})}}{S_{\text{OC}}^{\text{max}}-S_{\text{OC}}^{\text{min}}}$
${{\delta }^{\max }}=\max \{\delta ({{S}_{\text{OC}}})\}$

式中,$\beta ({{S}_{\text{OC}}})$$\delta ({{S}_{\text{OC}}})$分别为分段线性化后模型和原函数之间的绝对误差和相对误差;${{\delta }^{\text{avg}}}$${{\delta }^{\max }}$分别对应相对误差的均值和最大值,作为衡量分段线性化模型精确度的指标。

4 风光储微网随机日前调度模型

4.1 风光出力场景

通过第2.2节的模型可以得到风电出力和光伏出力的预测值,但风电出力和光伏出力都具有不确定性,为了更准确地预测,利用叠加预测误差的方法对风电和光伏的出力进行建模

${{\tilde{P}}_{\text{W}}}(t)\text{=}{{P}_{\text{W}}}(t)+{{\varepsilon }_{\text{W,}t}}$
${{\tilde{P}}_{\text{PV}}}(t)\text{=}{{P}_{\text{PV}}}(t)+{{\varepsilon }_{\text{PV,}t}}$

式中,${{\varepsilon }_{\text{W,}t}}$${{\varepsilon }_{\text{PV,}t}}$分别为风电出力和光伏出力的预测误差,不确定性就体现在这一项中,一般认为风光出力的短期预测误差服从正态分布[12]

为了将风光出力的不确定性加入到风光储微网随机日前调度模型中,本文采取基于场景的随机优化方法,生成足够多的风光出力场景,使用K-means算法进行聚类,得到对应的风光出力场景和概率,聚类后场景的个数通过手肘法确定。

4.2 目标函数

并网型风光储微网在实际运行中会产生和电网交互的费用以及电池储能寿命损耗成本,为了提高新能源的消纳率以及减少碳排放,本文将弃风弃光成本和削负荷成本加入到微网的运行成本中。根据碳汇流的理论,从电网购得的电能在生产阶段有部分是火电机组产生的,会造成一定的$\text{C}{{\text{O}}_{\text{2}}}$排放,所以在本文的并网型风光储微网日前调度模型中需要考虑这部分的碳排放成本。

结合以上的分析,以并网型风光储微网的日运行成本最小化为目标,将一天按照15 min一个点进行离散分析,即单位时间间隔$\Delta t\text{=0}\text{.25}\ \text{h}$,可以得到对应的目标函数为

$\min C=\sum\limits_{s=1}^{S}{Pr(s)\cdot \sum\limits_{t=1}^{T}{{{C}^{s}}(t)}}$
${{C}^{s}}(t)=C_{\text{grid}}^{s}(t)+C_{\text{q}}^{s}(t)+C_{\text{x}}^{s}(t)+C_{\text{BESS}}^{s}(t)+C_{\text{c}{{\text{o}}_{\text{2}}}}^{s}(t)$

式中,T为优化的总时段数;$s$为场景编号;S为总场景数;$Pr(s)$为场景$s$发生的概率;$C$为微网的日运行总成本;${{C}^{s}}(t)$$C_{\text{grid}}^{s}(t)$$C_{\text{q}}^{s}(t)$$C_{\text{x}}^{s}(t)$$C_{\text{bess}}^{s}(t)$$C_{\text{c}{{\text{o}}_{\text{2}}}}^{s}(t)$分别为$t$时段场景$s$下微网的运行总成本、电网交互费用、弃风弃光成本、削负荷成本、电池储能寿命损耗成本和碳排放成本。

(1) 微网和电网交互费用

$C_{\text{grid}}^{s}(t)\text{=}P_{\text{grid}}^{s}(t)\cdot {{p}_{\text{grid}}}(t)\cdot \Delta t$

式中,$P_{\text{grid}}^{s}(t)$$t$时段场景$s$下微网和电网交互的功率;${{p}_{\text{grid}}}(t)$为对应的$t$时段微网与电网交互功率的电价。当$P_{\text{grid}}^{s}(t)\ge 0$时,$C_{\text{grid}}^{s}(t)$对应微网从电网购电的费用;当$P_{\text{grid}}^{s}(t)<0$时,$C_{\text{grid}}^{s}(t)$对应微网向电网售电的收益。

(2) 微网的弃风弃光成本

$C_{\text{q}}^{s}(t)\text{=}P_{\text{q}}^{s}(t)\cdot {{c}_{\text{q}}}\cdot \Delta t$

式中,$P_{\text{q}}^{s}(t)$$t$时段场景$s$下微网弃风弃光的功率; cq为对应的单位弃风弃光成本。

(3) 微网的削负荷成本

$C_{\text{x}}^{s}(t)\text{=}P_{\text{x}}^{s}(t)\cdot {{c}_{\text{x}}}\cdot \Delta t$

式中,$P_{\text{x}}^{s}(t)$$t$时段场景$s$下微网削负荷的功率;${{c}_{\text{x}}}$为对应的单位削负荷成本。

(4) 微网的电池储能寿命损耗成本

$C_{\text{BESS}}^{s}(t)\text{=}L_{\text{loss}}^{s}(t)\cdot {{C}_{\text{B}}}$
${{C}_{\text{B}}}\text{=}{{C}_{\text{bat}}}\text{+}{{C}_{\text{pcs}}}$
${{C}_{\text{bat}}}\text{=}{{c}_{\text{E}}}\cdot {{E}_{\text{rated}}}$
${{C}_{\text{pcs}}}\text{=}{{c}_{\text{p}}}\cdot {{P}_{\text{rated}}}$

式中,$L_{\text{loss}}^{s}(t)$$t$时段场景$s$下微网的电池储能寿命损耗,通过第2节模型求取;${{C}_{\text{B}}}$为电池及变流器总投资成本;${{C}_{\text{bat}}}$为电池投资成本;${{c}_{\text{E}}}$为单位容量电池的投资成本;${{C}_{\text{pcs}}}$为变流器投资成本;${{c}_{\text{p}}}$为单位功率变流器的投资成本;${{P}_{\text{rated}}}$为电池的额定功率。

(5) 微网的碳排放成本

$C_{\text{c}{{\text{o}}_{\text{2}}}}^{s}(t)\text{=}{{k}_{\text{H}}}\cdot P_{\text{equ}}^{s}(t)\cdot {{c}_{\text{c}{{\text{o}}_{\text{2}}}}}\cdot \Delta t$

式中,$P_{\text{equ}}^{s}(t)$$t$时段场景$s$下和电网交互功率中等效的火电机组出力;${{k}_{\text{H}}}$为火电机组的$\text{C}{{\text{O}}_{\text{2}}}$排放强度;${{c}_{\text{c}{{\text{o}}_{\text{2}}}}}$为对应的排放单位重量$\text{C}{{\text{O}}_{\text{2}}}$产生的成本。

4.3 约束条件

4.3.1 电池储能运行约束

(1) 电池储能容量约束

$S_{\text{OC}}^{\text{min}}\le S_{\text{OC}}^{s}(t)\le S_{\text{OC}}^{\text{max}}$

式中,$S_{\text{OC}}^{\text{min}}$$S_{\text{OC}}^{\text{max}}$分别对应整个研究阶段电池荷电状态的最小值和最大值。

(2) 电池储能功率约束

$P_{\text{B}}^{s}(t)=P_{\text{d}}^{s}(t)-P_{\text{c}}^{s}(t)$

式中,$P_{\text{B}}^{s}(t)$$P_{\text{c}}^{s}(t)$$P_{\text{d}}^{s}(t)$分别对应$t$时段场景电池的总功率和充、放电功率。

充电时刻满足功率约束

$\text{0}\le P_{\text{c}}^{s}(t)\le {{P}_{\text{rated}}}\ \ \ \ \ \ \ P_{\text{d}}^{s}(t)=0$

放电时刻满足功率约束

$\text{0}\le P_{\text{d}}^{s}(t)\le {{P}_{\text{rated}}}\ \ \ \ \ \ P_{\text{c}}^{s}(t)=0$

4.3.2 电网交互功率约束

微网和电网交互功率为

$-P_{\text{grid}}^{\text{max}}\le P_{\text{grid}}^{s}(t)\le P_{\text{grid}}^{\text{max}}$

式中,$P_{\text{grid}}^{\text{max}}$为微网和电网交互功率的上限。

4.3.3 功率平衡约束

$P_{\text{W}}^{s}(t)\text{+}P_{\text{PV}}^{s}(t)\text{+}P_{\text{B}}^{s}(t)\text{+}P_{\text{grid}}^{s}(t)\text{=}{{P}_{\text{load}}}(t)$

式中,${{P}_{\text{load}}}(t)$为对应的$t$时段的负荷;$P_{\text{W}}^{s}(t)$$P_{\text{PV}}^{s}(t)$分别为$t$时段场景$s$下叠加预测误差的风电和光伏出力。

当电源出力超过负荷时,需要进行弃风弃光

$P_{\text{q}}^{s}(t)\text{=}P_{\text{W}}^{s}(t)\text{+}P_{\text{PV}}^{s}(t)\text{+}P_{\text{B}}^{s}(t)\text{+}P_{\text{grid}}^{s}(t)-{{P}_{\text{load}}}(t)$

当电源出力小于负荷时,需要进行削负荷

$P_{\text{x}}^{s}(t)\text{=}{{P}_{\text{load}}}(t)-P_{\text{W}}^{s}(t)-P_{\text{PV}}^{s}(t)-P_{\text{B}}^{s}(t)-P_{\text{grid}}^{s}(t)$

4.3.4 等效出力约束

$P_{\text{equ}}^{s}(t)\text{=}\alpha \cdot P_{\text{grid}}^{s}(t)$

式中,$\alpha $为等效系数,由研究地区火电机组装机容量和所有发电机组装机容量的比例进行估算。

5 算例分析

5.1 算例设置

本文并网型风光储微网中配置风电100 kW、光伏200 kW,根据青海地区典型日的气象数据生成对应风光预测出力,选取同一时间对应的负荷曲线进行分析,如图3所示。通过第4.1节中的风光出力场景模型对预测值叠加预测误差,在本算例中生成1 500个随机风光出力场景,通过K-means算法进行聚类,经过手肘法验证四个聚类场景的效果比较好。

图3

图3   风电、光伏预测出力曲线以及负荷曲线


依据实际的风光储微网工程中各元件容量,对应配置该算例的电池储能参数,如表1所示。根据文献[13]中LI TEL48-170C磷酸铁锂电池的试验数据,处理得到对应电池储能的循环寿命拟合曲线为

${{N}_{\text{life}}}({{D}_{\text{OD}}})=21\ 870\times \exp \left( -1.957{{D}_{\text{OD}}} \right)$

通过式(15)~(18)计算分段线性化后模型的精确度,试验分析当电池储能寿命损耗模型分段数取5时,${{\delta }^{\text{avg}}}=0.42%$,${{\delta }^{\max }}=1.37%$,有好的求解精确度以及求解速度,所以在算例分析中分段数取为5。

表1   电池储能的参数

参数数值
额定容量${{E}_{\text{rated}}}$/(kW·h)400
额定功率${{P}_{\text{rated}}}$/kW100
充、放电效率${{\eta }_{\text{c}}}$,${{\eta }_{\text{d}}}$0.95
荷电状态最大值$S_{\text{OC}}^{\text{max}}$0.9
荷电状态最小值$S_{\text{OC}}^{\text{min}}$0.2
单位容量投资成本${{c}_{\text{E}}}$/[元/(kW·h)]1 000
单位功率变流器投资成本${{c}_{\text{p}}}$/(元/kW)1 500

新窗口打开| 下载CSV


根据算例中风光的配置容量,设置系统与上级电网间的交互功率上限为50 kW,系统和上级电网之间购售电使用青海地区的实际峰谷分时电价数据,取0.14元/(kW·h)(0:00—8:00)、0.35元/(kW·h) (8:00—9:00,12:00—18:00,23:00—24:00)及0.55元/(kW·h)(9:00—12:00,18:00—23:00);根据火电机组占总装机容量的比例,在本算例中取等效系数$\alpha =0.5$,根据文献[14]的算例参数可知,火电机组$\text{C}{{\text{O}}_{\text{2}}}$排放强度为0.89 kg/(kW·h),对应的碳排放成本为0.14元/kg$\text{C}{{\text{O}}_{\text{2}}}$;根据文献[15]的算例参数可知,系统的弃风弃光成本与削负荷成本分别取1.0 元/(kW·h)和8.0元/(kW·h)。

5.2 求解方法

为了解决风光储微网随机日前调度中决策变量较多的问题,一般采用对优化模型线性化后,利用商业软件CPLEX,进行有效快速求解。

(1) 电池储能运行模型线性化

$0\le P_{\text{c}}^{s}(t)\le v_{\text{c}}^{s}(t)\cdot {{P}_{\text{rated}}}$
$0\le P_{\text{d}}^{s}(t)\le v_{\text{d}}^{s}(t)\cdot {{P}_{\text{rated}}}$
$S_{\text{OC}}^{s}(t)=S_{\text{OC}}^{s}(t-1)+\frac{P_{\text{c}}^{s}(t)\cdot {{\eta }_{\text{c}}}}{{{E}_{\text{rated}}}}\cdot \Delta t-\frac{P_{\text{d}}^{s}(t)}{{{\eta }_{\text{d}}}\cdot {{E}_{\text{rated}}}}\cdot \Delta t$
$v_{\text{c}}^{s}(t)+v_{\text{d}}^{s}(t)\le 1$
$S_{\text{OC}}^{\text{min}}\le S_{\text{OC}}^{s}(t)\le S_{\text{OC}}^{\text{max}}$

式中,$v_{\text{c}}^{s}(t)$$v_{\text{d}}^{s}(t)$为0-1变量;$v_{\text{c}}^{s}(t)=1$$v_{\text{d}}^{s}(t)=1$分别表示电池处于充电和放电状态。

(2) 功率平衡约束线性化。对于弃风弃光成本$C_{\text{q}}^{s}(t)$以及削负荷成本$C_{\text{x}}^{s}(t)$,可以统一写为功率不平衡成本$C_{\text{unbal}}^{s}(t)$

$C_{\text{unbal}}^{s}(t)=\left\{ \begin{align} & C_{\text{q}}^{s}(t)P_{\text{q}}^{s}(t)\ge 0 \\ & C_{\text{x}}^{s}(t)P_{\text{x}}^{s}(t)\ge 0 \\ \end{align} \right.$

利用大M法将其线性化[16]

$M\cdot (z(t)-1)\le P_{\text{W}}^{s}(t)\text{+}P_{\text{PV}}^{s}(t)\text{+}P_{\text{B}}^{s}(t)\text{+}P_{\text{grid}}^{s}(t)-{{P}_{\text{load}}}(t)\le M\cdot z(t)$
$C_{\text{x}}^{s}(t)-M\cdot z(t)\le C_{\text{unbal}}^{s}(t)\le C_{\text{x}}^{s}(t)+M\cdot z(t)$
$C_{\text{q}}^{s}(t)-M\cdot (1-z(t))\le C_{\text{unbal}}^{s}(t)\le C_{\text{q}}^{s}(t)+M\cdot (1-z(t))$

式中,$z(t)$为辅助0-1变量;M为一个很大的正数。

结合式(10)~(14)、(21)~(23)、(26)~(30)和(39)~(49)分别对储能寿命损耗模型、储能运行模型以及目标函数中的弃风弃光成本和削负荷成本进行线性化,从而可以很好地嵌入整个优化模型中,得到整个线性化的风光储微网随机日前调度模型,使调度问题有着较好的求解速率和求解准确度。

5.3 算例结果

5.3.1 电池储能经济性分析

为了体现加入电池储能的经济性,将含电池储能微网和无电池储能微网的成本进行对比分析,如表2所示。

表2   有/无电池储能的微网运行成本对比结果

运行成本模式
含电池储能无电池储能
和电网交互费用/元38.4955.45
储能寿命损耗成本/元84.880
削负荷成本/元540.582 559.90
弃风弃光成本/元7.61302.63
碳排放成本/元30.9731.51
总运行成本/元702.532 949.49

新窗口打开| 下载CSV


在并网型风光微网中引入电池储能,可以有效减少整个微网的日运行成本,从2 949.49元下降到了702.53元,从而提高了风光储微网的经济性。

为了分析电池储能对风光储微网经济性方面的作用,结合微网的电池具体运行策略进行分析,如图4所示。

图4

图4   电池功率调度曲线


在弃风弃光成本和削负荷成本方面,4:30—8:00以及19:30—23:00阶段电池通过放电运行,有效缓解负荷远超风光出力的情况,降低了削负荷成本;在10:30—16:30阶段电池通过充电运行,提升了风光的消纳率,减少了风光出力远超负荷的情况,降低了弃风弃光成本。

在和电网交互费用方面,由于当地采用分时电价的缘故,电池储能可以实现一种变相的峰谷套利,在1:00—1:30负荷超过风光出力的阶段,通过自身的放电能力提供这一部分差额,在电网留有余力的情况下(即在该时刻不需要产生削负荷成本),电池在谷时电价下进行充电,在传统峰谷套利模式下,这部分充电量可以用在峰时电价阶段去向电网售电获利;但由于本文中削负荷单位成本较高,优先考虑降低削负荷成本实现变相的峰谷套利,这也是含电池储能微网和无电池储能微网在和电网交互费用上产生区别的缘故。

在碳排放成本方面,风光储微网模型中该成本只和购电量有关,而与售电量无关,根据上述分析在考虑电池储能峰谷套利的模式下会在一些阶段产生额外的购电量,这也就造成了两种模式下碳排放成本相差不大的原因。

5.3.2 电池储能寿命损耗分析

为了验证第3节储能寿命损耗模型在并网型风光储微网场景中的有效性,将不考虑电池储能寿命损耗的微网模型同考虑电池储能寿命损耗的模型进行对比,优化结果如表3所示。

表3   电池储能寿命损耗对微网运行成本影响结果

运行成本模式
计及电池寿命损耗不计及电池寿命损耗
和电网交互成本/元38.4934.91
储能寿命损耗成本/元84.880
削负荷成本/元540.58539.88
弃风弃光成本/元7.617.45
碳排放成本/元30.9735.08
总运行成本/元702.53617.32

新窗口打开| 下载CSV


由于单位电量的弃风弃光成本和削负荷成本都超过单位电量的电池储能寿命损耗和购电单价,所以在风光储微网中无论是否考虑电池储能寿命损耗,都会优先考虑提高风光消纳率及减少负荷超额量,体现在成本中就是这两种模式下削负荷成本及弃风弃光成本差异较小,说明在电池容量固定的情况下,风光消纳率的提升以及超额负荷的削减已经到了电池的最大优化程度,所以这两种模式的主要区别在于电池储能参与传统峰谷套利模式的程度。由于考虑了寿命损耗会减少电池的过充过放,所以从经济性看不考虑寿命损耗的电池会通过套利的模式增加对电网的售电收益,在该情况下电池储能电量会有更大的波动范围。

结合两种模式电池的SOC情况具体进行说明,从图5中可以直观体现考虑电池储能寿命损耗模型后,电池的SOC曲线会比不计及电池储能寿命损耗的模型的SOC曲线平缓,减少了电池充放电量波动大的出力方式,缓解了电池频繁充放电的情况,根据第2节的模型可以计算得出,按照考虑电池储能寿命损耗模型的微网调度运行方式,电池还可以在该算例的情况下工作大约6 480 d,而按照不考虑寿命损耗的运行模式,电池在该算例的情况下只能工作大约5 764 d,提升比例达到12.42%,体现了本文所提模型可以通过规范储能的充放电动作有效提高电池的使用寿命,所以在风光储微网的经济性分析中考虑电池储能的寿命损耗是具有必要性的。

图5

图5   电池储能寿命损耗对电池调度策略影响图


5.3.3 电池储能容量配置灵敏度分析

在并网型风光储微网中电池储能容量决定了其可以参与削负荷、消纳风光资源、峰谷套利以及影响碳排放量的程度,所以本节分析不同容量配置对微网各项运行成本的影响,优化结果如表4所示。

表4   电池储能不同容量配置下微网运行成本结果

电池储能配置
容量/(kW·h)
和电网交互
费用/元
储能寿命损耗
成本/元
削负荷成本/元弃风弃光成本/元碳排放成本/元总运行成本/元
30059.4369.43768.5881.3232.831 011.59
35049.0277.35654.7044.1831.90857.15
40038.4984.88540.587.6130.97702.53
45037.4092.79426.56030.06586.81
50040.5198.76312.48029.11480.86

新窗口打开| 下载CSV


随着电池配置容量的提高,尽管带来了电池储能寿命损耗成本的上升,但其参与整个并网型风光储微网运行的调节能力变强了,尤其是电池在削负荷以及消纳风光资源方面的巨大作用得以体现,使得整个微网的经济性运行水平有了显著的提高。但随着电池配置容量的提高,和电网交互费用有了一个先下降后上升的过程,这是因为在峰时电价时,电池通过更强的放电能力减少了购电量,而在谷时电价时,由于本算例中弃风弃光和削负荷的单位成本都要高于电池向电网售电单价,所以将更强的充电能力应用于消纳光伏以及减少负荷超额量,实现一种变相的削峰填谷。

进一步分析电池储能引入对微网经济效益的影响,研究电池配置容量和微网总运行成本的关系。从图6发现随着电池容量的上升,微网总运行成本为一个先下降后平稳的趋势,而电池寿命损耗成本占总运行成本的比重在不断上升,这说明电池储能引入可以有效地提升微网的经济性,但随着电池容量上升到一定阶段,实现了负荷超额量完全削减和风光完全消纳,电池在传统峰谷套利模式下的收入和寿命损耗成本逐渐平衡,使微网的经济性趋于稳定。

图6

图6   总运行成本及寿命损耗成本占比与电池容量关系


本文中电池储能容量配置灵敏度的算例分析,对风光储微网实际案例中的电池储能容量配置问题具有参考意义,如图6所示,在电池容量上升的后段,微网总运行成本趋于稳定,所以利用本文考虑电池储能寿命损耗的风光储微网随机日前调度策略可以为实际风光储系统配置最优电池储能容量提供可行思路。

6 结论

本文将精细化的电池储能寿命损耗模型加入到风光储微网随机日前调度模型中,将电池在提高风光消纳率、削峰填谷、减少碳排放量方面的作用转化为具体的经济性指标,综合各项经济指标和电池储能寿命损耗成本探究电池储能对风光储微网的经济作用,将储能寿命损耗模型、储能运行模型以及目标函数中的弃风弃光成本和削负荷成本进行线性化,从而使调度问题有着较好的求解速率和求解准确度。通过算例分析对文中所提的风光储微网日前调度模型进行验证,所得结论如下所述。

(1) 电池储能提升风光储微网的经济性主要体现在提高风光消纳率以及削减超额负荷方面,可以有效降低弃风弃光成本以及削负荷成本。

(2) 本文所提模型可以规范电池储能的充放电动作,避免电池的过度充电/放电,从而延长电池的使用寿命,经计算得到,本文所提模型相比于不考虑电池寿命损耗的模型,可延长其寿命12.42%。

(3) 随着电池容量的上升,微网的日运行成本有一个先下降后平稳的变化趋势,说明在考虑电池储能寿命损耗费用后,本文所提的优化调度方法仍具有较高的经济性,可以为风光储系统配置最优电池储能容量提供可行思路。

参考文献

米剑锋, 马晓芳.

中国CCUS技术发展趋势分析

[J]. 中国电机工程学报, 2019, 39(9):2537-2544.

[本文引用: 1]

MI Jianfeng, MA Xiaofang.

Development trend analysis of carbon capture,utilization and storage technology in China

[J]. Proceedings of the CSEE, 2019, 39(9):2537-2544.

[本文引用: 1]

苏适, 周立栋, 陆海, .

基于改进混沌粒子群算法的多源独立微网多目标优化方法

[J]. 电力系统保护与控制, 2017, 45(23):34-41.

[本文引用: 1]

SU Shi, ZHOU Lidong, LU Hai, et al.

Multi-objective optimization method of multi-source independent microgrid based on modified CPSO

[J]. Power System Protection and Control, 2017, 45(23):34-41.

[本文引用: 1]

谭颖, 吕智林, 李捷.

基于改进ELM的风/光/柴/储独立微网分布式电源多目标容量优化配置

[J]. 电力系统保护与控制, 2016, 44(8):63-70.

[本文引用: 1]

TAN Ying, Zhilin, LI Jie.

Multi-objective optimal sizing method for distributed power of wind-solar-diesel-battery independent microgrid based on improved electromagnetism-like mechanism

[J]. Power System Protection and Control, 2016, 44(8):63-70.

[本文引用: 1]

叶亮, 吕智林, 王蒙, .

基于最优潮流的含多微网的主动配电网双层优化调度

[J]. 电力系统保护与控制, 2020, 48(18):27-37.

[本文引用: 1]

YE Liang, Zhilin, WANG Meng, et al.

Bi-level programming optimal scheduling of ADN with a multi-microgrid based on optimal power flow

[J]. Power System Protection and Control, 2020, 48(18):27-37.

[本文引用: 1]

ZHANG Cuo, XU Yan, DONG Zhaoyang, et al.

Robust operation of microgrids via two-stage coordinated energy storage and direct load control

[J]. IEEE Transactions on Power Systems, 2017, 32(4):2858-2868.

DOI:10.1109/TPWRS.2016.2627583      URL     [本文引用: 1]

丁宇劼.

基于可再生能源与储能系统的微电网模型及其优化运行仿真研究

[J]. 电气工程学报, 2020, 15(1):36-40.

[本文引用: 1]

DING Yujie.

Microgrid model with renewable energy and energy storage system and its optimal operation simulation

[J]. Journal of Electrical Engineering, 2020, 15(1):36-40.

[本文引用: 1]

WANG D, GUAN X, WU J, et al.

Integrated energy exchange scheduling for multi-microgrid system with electric vehicles

[J]. IEEE Transactions on Smart Grid, 2016, 7(4):1762-1774.

DOI:10.1109/TSG.2015.2438852      URL     [本文引用: 1]

SHI Y, XU B, WANG D, et al.

Using battery storage for peak shaving and frequency regulation:Joint optimization for superlinear gains

[J]. IEEE Transactions on Power Systems, 2018, 33(3):2882-2894.

DOI:10.1109/TPWRS.2017.2749512      URL     [本文引用: 1]

WU K, ZHOU H, AN S, et al.

Optimal coordinate operation control for wind-photovoltaic-battery storage power-generation units

[J]. Energy Conversion and Management, 2015, 90(9):466-475.

DOI:10.1016/j.enconman.2014.11.038      URL     [本文引用: 1]

王开科, 南东亮, 李勇, .

虚拟电厂参与大规模新能源系统的储能侧双层优化运行策略

[J]. 电气工程学报, 2020, 15(2):24-33.

[本文引用: 1]

WANG Kaike, NAN Dongliang, LI Yong, et al.

Bi-level optimal operation strategy of energy storage with large-scale new energy system under the participation of virtual power plant

[J]. Journal of Electrical Engineering, 2020, 15(2):24-33.

[本文引用: 1]

孙运志, 蒋德玉, 张盛林, .

计及电池寿命损耗的多能源微网储能优化配置

[J]. 电力系统及其自动化学报, 2021, 33(5):128-133.

[本文引用: 1]

SUN Yunzhi, JIANG Deyu, ZHANG Shenglin, et al.

Optimal energy-storage configuration of multi-energy microgrid considering battery life degradation

[J]. Proceedings of the CSU-EPSA, 2021, 33(5):128-133.

[本文引用: 1]

吴江. 考虑风电光伏出力不确定性的微网经济优化调度[D]. 南京: 南京邮电大学, 2020.

[本文引用: 3]

WU Jiang. Economic optimal scheduling of microgrid considering uncertainty of wind power photovoltaic output[D]. Nanjing: Nanjing University of Posts and Telecommunications, 2020.

[本文引用: 3]

Featherweight lithium ion battery[EB/OL]. [2020-11-05]. http://www.cdtechno.com/pdf/lit/12_1093_0217.pdf.

URL     [本文引用: 1]

段建民, 王志新, 王承民, .

考虑碳减排效益的可再生电源规划

[J]. 电网技术, 2015, 39(1):11-15.

[本文引用: 1]

DUAN Jianmin, WANG Zhixin, WANG Chengmin, et al.

Renewable power planning considering carbon emission reduction benefits

[J]. Power System Technology, 2015, 39(1):11-15.

[本文引用: 1]

LIU G, TOMSOVIC K.

Quantifying spinning reserve in systems with significant wind power penetration

[J]. IEEE Transactions on Power Systems, 2012, 27(4):2385-2393.

DOI:10.1109/TPWRS.2012.2207465      URL     [本文引用: 1]

DING T, BO R, SUN H, et al.

Big-M based MIQP method for economic dispatch with disjoint prohibited zones

[J]. IEEE Transactions on Power System, 2014, 29(2):976-977.

DOI:10.1109/TPWRS.2013.2287993      URL     [本文引用: 1]

/