风电全消纳下的配电网储能可调鲁棒优化配置

Adjustable Robust Optimization for Energy Storage System in Distribution Network Based on Wind Power Full Accommodation

朱嘉远, 刘洋, 许立雄, 蒋卓臻, 林潇

四川大学 电气信息学院,四川省 成都市 610065

ZHU Jiayuan, LIU Yang, XU Lixiong, JIANG Zhuozhen, LIN Xiao

College of Electrical Engineering and Information Technology, Sichuan University, Chengdu 610065, Sichuan Province, China

  • 朱嘉远(1993),男,硕士研究生,研究方向为储能在含配电网中的应用,E-mail:zhujy@stu.scu. edu.cn;

  • 刘洋(1982),男,博士,副教授,研究方向为电力系统高性能计算及算法优化,E-mail:yang.liu @scu.edu.cn;

  • 许立雄(1982),男,博士,讲师,通信作者,研究方向为电力系统稳定与控制,E-mail:xulixiong @163.com。

文章编号: 1000-3673(2018)06-1875-08 中图分类号: TM73

摘要

随着现代电力系统中分布式风电的渗透率不断提高,其出力的不确定性和波动性对电网规划运行的影响愈发显著。侧重于对配电网中储能系统进行优化配置以应对分布式风电的不确定性,提出了综合考虑网损成本、储能投资运行成本以及主网购电成本的可调鲁棒规划模型,最终获得满足风电功率全部消纳要求且总成本最少的储能布点以及容量配置方案。在规划求解过程中,引入二阶锥松弛技术对原模型进行凸化松弛,使非凸非线性优化问题转化成为可有效求解的二阶锥优化问题,并通过对凸锥模型进行拉格朗日对偶将原双层规划问题转化成为含双线性项的单层规划问题。由于鲁棒调整参数的引入,整个模型在提高规划结果鲁棒性的同时,防止规划结果过于保守。最后通过算例证明了所提出模型的有效性。

关键词 : 风电不确定性; 储能选址定容; 二阶锥规划; 可调鲁棒规划;

DOI:10.13335/j.1000-3673.pst.2017.2621

ABSTRACT

With increasing penetration of distributed wind generation (DWG) in modern power system, uncertainty and volatility of wind power have a significant effect on power grid planning. This paper focuses on configuration of energy storage system (ESS) in distribution network to accommodate uncertainty of wind power. An adjustable robust optimization planning method is proposed to obtain minimum capacity and appropriate distribution of ESS. The costs of network loss, ESS and main network purchase are all taken into account. In order to simplify optimization calculations, the second-order cone relaxation technique is implemented to transform non-convex nonlinear optimization problem into a second-order cone programming. Then, Lagrange duality theory is employed to convert the bi-level convex programming into a monolayer programming problem with bi-linear terms. In addition to improving robustness of the planning results, the model can prevent the planning results from being overly conservative because of introduction of robust adjustment parameter. Finally, experiment results indicate effectiveness of the proposed method.

KEY WORDS : uncertainty of wind power; energy storage system planning; second-order cone programming; adjustable robust planning;

0 引言

近年来,随着能源需求不断增加以及环境问题日益严峻,电力系统中分布式风电(distributed wind generation,DWG)的渗透率逐渐提高。储能系统(energy storage system,ESS)具有快速响应以及对功率进行双向调节的能力,可以实现电能的时空平移,而被广泛地用于应对DWG出力不确定性及波动性问题。但由于ESS造价昂贵,如何对其进行合理的选址定容从而获得较高经济收益成为现阶段的研究重点。

目前,对电力系统中ESS进行优化配置实现功率调控从而提高风电消纳能力的研究主要可以分为两方面。一方面,在DWG侧配置ESS[1-7],将两者视为一个整体,根据要求对ESS进行充放电控制,对风电预测出力进行平滑或跟踪,使其满足一定输出要求。此种方法主要针对风电侧出力的调整,而对系统的整体考虑较少。本文侧重于另一方面,即在含DWG的电力系统模型中加入ESS约束,寻找其最优的布点及功率容量配置方案,在提高风电消纳能力的同时,使整个系统运行更具经济性及充裕性[8-14]。在ESS的优化配置过程中,对于风电场所存在的不确定出力,文献[3]将其用模糊变量表示,提出了基于模糊相关机会规划的储能优化控制方法,该方法采用软约束形式,且规划中模糊隶属度函数常根据有限的样本及决策者经验来确定,使得规划结果不够精确。文献[7-9]采用机会约束优化方法,该方法也允许所作决策在一定程度上不满足约束条件,决策者可以根据实际需求调控风险,相比确定模型提高了灵活性,但当系统规模增大时,模型的计算时间普遍较长[15]。文献[10-13]采用确定性场景规划法,将不确定性因素可能的取值通过某种规则选取出来并进行枚举,由此不确定性问题转变成为确定性问题,然而构成具有代表性的场景且保证模型的精度往往较为困难。文献[14]采用两阶段随机优化,对风电柴油系统中ESS的容量进行合理配置,其在求解过程中依赖于精确的概率密度函数,但由于测量技术、经济性或实际条件所限,决策者往往做出一定程度的近似和假设,使得规划结果易产生较大误差。

鲁棒规划由于仅需不确定性参数的边界,避免了对精确的概率密度函数抑或模糊函数的需求,在电力系统规划运行中得到了广泛应用[16-18]。文献[16]采用不确定参数均值以及波动区间描述不确定性,基于满足风电功率全部接纳的要求,引入鲁棒线性优化方法,建立了包含储能容量、功率以及位置的模型。鲁棒规划方法通过将随机变量在确定性范围内的最差场景纳入模型中得出规划方案,保证了最优解在整个随机变量确定性变化范围内的鲁棒性。

本文基于风电出力完全消纳,提出了配电网储能系统可调鲁棒规划模型,通过对配电网(distribution network,DN)中ESS进行合理选址定容,达到减小网损,降低主网购电成本的目的,并且相比于确定性规划拥有更强的抵抗不确定风险的能力。针对DN的非凸非线性规划问题,本文采用二阶锥松弛技术(second-order cone relaxation,SOCR)将其凸化松弛为一个可有效求解的混合整数二阶锥优化问题,然后利用拉格朗日对偶将其进一步转化为等价的单层凸规划问题,从而便于利用商业优化软件求得最优解。最后通过算例仿真验证了本文方法的正确性。

1 风电全消纳下的储能可调鲁棒规划模型

目前较多关于含DWG电力系统中的ESS优化配置研究均涉及提高可再生能源的接纳能力,并多以经济效益最优为规划目标。本文基于风机出力完全消纳,以系统投资运行综合成本最小为目标,构建模型以获得储能系统的最佳配置位置及容量。在优化配置过程中,为保证规划结果的鲁棒性,通常制定最恶劣场景下的规划方案并予以实施,但容易造成结果过于保守,因此本文引入可调参数对不确定因素进行限制,从而得出更具经济效益的储能选址定容方案。可调鲁棒模型表达如下:

$\left\{ \begin{align} \underset{x}{\mathop{\min }}\,\underset{u}{\mathop{\max }}\,C(x,u,\Gamma ) \\ \text{s}\text{.t}\text{.}H(x,u,\Gamma )=0 \\ G(x,u,\Gamma )\le 0 \\\end{align} \right.$ (1)

式中:x表示优化变量;u表示不确定参数;Γ表示可调参数,用于对不确定参数进行限制。目标函数C为系统经济效益的综合表达式;HG分别为模型的等式约束及不等式约束,主要包括配电网中功率平衡约束,电流电压约束以及储能运行约束等。

1.1 不确定参数

根据文献[19],上述模型中不确定参数u的主要体现为与不确定风速相关的风电出力Pwind,i,t。本文应用文献[19-20]的不确定集合表征方法,将风速不确定集表示为包含预测值vpre i,t及波动量vwidth i,t的区间形式如下所示:

${{v}_{i,t}}\in [v_{i,t}^{\text{pre}}-v_{i,t}^{\text{width}},v_{i,t}^{\text{pre}}+v_{i,t}^{\text{width}}\text{ }\!\!]\!\!\text{ }$ (2)

由于风电输出功率与风速之间满足式(3)的分段函数关系,风电出力不确定区间可表示如式(4)。

式中:Pwind,i,t表示节点i处风机t时刻的输出功率;vi,t为节点it时刻风速;vcut in和vcut out分别为风机的切入风速和切出风速;vrated为风机的额定风速;φi为节点i处风电的装机容量;Pwind,i(vi,t)maxPwindi(vi,t)min分别为风电输出功率的上下限。

由于引入可调鲁棒参数Γ对不确定区间做出部分妥协来保证规划结果的经济性,风电场各时段输出功率还需满足如下约束:

式中:Ppre wind,i,t为风电t时刻预测出力;PtP- t分别为风电出力上下波动范围;ρt以及ρ- t为0-1变量,若风机出力大于预测出力ρt=1,ρ- t=0,反之ρt=0,ρ- t=1。

1.2 目标函数

目标函数主要包括一个运行周期内储能投资和运行维护成本Cess、网损成本Cpowerloss以及主网购电成本Cgrid。以上各部分具体表述如下:

式中:第1项为ESS一个运行周期的投资运行成本,由年投资运行成本折算而来:T为运行周期的小时数;累和项是ESS的年投资运行成本,ΩE表示所有配置ESS的节点;αess为ESS资金等年值系数,具体表达如式(7)所示;ress为ESS贴现率;hess为ESS使用寿命;IC ess为ESS单位容量投资成本;Eess,i为节点i处ESS装机容量;Gc ess为ESS单位容量年运行维护成本;第2项为网损成本;Kloss为网损电价;Iij,t为支路ijt时刻流过的电流;rij为支路ij的电阻;c(i)表示所有以i为首端节点的末端节点集合;N为总节点数;第3项为主网购电成本;ηt为分时电价;Pgrid,t为根节点在t时刻输入的有功功率。

1.3 约束条件
1.3.1 网络约束

\(\begin{align} {{P}_{\text{ess},i,t}}+{{P}_{\text{wind},i,t}}+{{P}_{\text{grid},i,t}}-{{P}_{\text{load},i,t}}\text{=} \\ \sum\limits_{i=1}^{N}{\sum\limits_{j\in c(i)}{U_{i,t}^{2}{{g}_{ij}}-{{U}_{i,t}}{{U}_{j,t}}({{g}_{ij}}\cos {{\theta }_{ij,t}}}}+{{b}_{ij}}\sin {{\theta }_{ij,t}}) \\\end{align}\) (8)

\(\begin{align} {{Q}_{\text{ess},i,t}}+{{Q}_{\text{grid},i,t}}-{{Q}_{\text{load},i,t}}\text{=} \\ \sum\limits_{i=1}^{N}{\sum\limits_{j\in c(i)}{-U_{i,t}^{2}b{}_{ij}+{{U}_{i,t}}{{U}_{j,t}}}({{b}_{ij}}\cos {{\theta }_{ij,t}}}-{{g}_{ij}}\sin {{\theta }_{ij,t}}) \\\end{align}\) (9)

$U_{i}^{\min }\le {{U}_{i,t}}\le U_{i}^{\max }$ (10)

\(\begin{align} I_{ij,t}^{2}=(U_{i,t}^{2}+U_{j,t}^{2}-2{{U}_{i,t}}{{U}_{j,t}}\cos {{\theta }_{ij,t}})\cdot \\ \text{ }(g_{ij}^{2}+b_{ij}^{2})\le I{{_{ij}^{\max }}^{2}} \\\end{align}\) (11)

式中:Ui,tUj,t分别为节点ijt时刻的电压幅值;θij,t为在t时刻节点i到节点 j 间的电压相角差;θij,t=θi,t-θj,t;gijbij分别为节点i到节点j间的电导和电纳;Pwind,i,t为节点i处的风机在t时刻输出的有功功率;Pess,i,tQess,i,t分别为节点i处的储能在t时刻输出的有功功率和无功功率;Pgrid,i,tQgrid,i,t分别为根节点在t时刻输出的有功功率和无功功率;Pload,i,tQload,i,t分别为节点i处的负荷在t时刻消耗的有功功率和无功功率;Umin iUmax i分别为节点i处最小和最大电压允许值;Iij,t为线路ijt时刻的电流;Imax ij为线路ij的最大允许电流。

1.3.2 储能约束

本文规划所用ESS为文献[21]所述蓄电池储能系统,由于换流器的存在,整个储能系统可提供一定无功辅助。同时假定ESS接入点换流器已提供,即换流器容量取为固定值。其中,ESS向配电网输电定为其功率正方向。根据文献[21],ESS需满足功率约束、荷电状态(state of charge,SOC)连续性约束以及容量约束如式(12)—(16),且为保证连续运行,ESS须遵循自身始末电量平衡约束,如式(17)。

\(\sqrt{Q_{\text{ess},i,t}^{2}+P_{\text{ess},i,t}^{2}}\le S_{\text{ess},i}^{\text{rate}}{{\varepsilon }_{i}}\) (12)

\(-P_{\text{ess},i}^{\text{charge},\max }{{\varepsilon }_{i}}\le {{P}_{\text{ess},i,t}}\le P_{\text{ess},i}^{\text{dis}ch\arg e,\max }{{\varepsilon }_{i}}\) (13)

${{e}_{\text{ess},i,t+\Delta t}}={{e}_{\text{ess},i,t}}-{{P}_{\text{ess},i,t}}\Delta t$ (14)

${{S}_{\text{ess},i,t}}=\frac{{{e}_{\text{ess},i,t}}}{E_{\text{ess},i}^{{}}}$ (15)

$S_{\text{ess},i}^{\min }\le {{S}_{\text{ess},i,t}}\le S_{\text{ess},i}^{\max }$ (16)

${{S}_{\text{ess},i,0}}={{S}_{\text{ess},i,T}}$ (17)

式中:εi为0-1变量,εi=1表示节点i处连有ESS,εi=0表示节点i处未连有ESS;Srate ess,i为换流器的额定容量;Pmax,charge ess,iPmax,discharge ess,i分别为ESS充、放电功率上限;eess,i,t为节点i处ESS在t时刻的容量;Sess,i,t为节点i处ESS在t时刻的荷电状态;Smin ess,iSmax ess,i分别为SOC上下限;T为某一固定的运行周期。

综上,式(12)—(17)构成了ESS出力数学模型。

2 基于二阶锥松弛及拉格朗日对偶的模型转化

上节所述问题已是一个复杂的混合整数非凸非线性规划问题,若在其基础上进一步考虑不确定参数并进行鲁棒优化会使模型求解更为困难。因此本文先采用二阶锥松弛技术将原模型凸化,然后利用鲁棒对偶转化方法将原双层凸规划问题转化为易于求解的单层凸规划问题。具体转化方法如下。

2.1 基于二阶锥松弛的模型凸化
2.1.1 目标函数及网络约束转化

采用文献[22]的方法,通过如下变量替代和松弛,使原规划问题转变为二阶锥规划问题。替换变量如下所示:

$\left\{ \begin{align} {{X}_{i,t}}=U_{i,t}^{2} \\ {{Y}_{ij,t}}={{U}_{i,t}}{{U}_{j,t}}\cos {{\theta }_{ij,t}} \\ {{Z}_{ij,t}}={{U}_{i,t}}{{U}_{j.t}}\sin {{\theta }_{ij,t}} \\\end{align} \right.$ (18)

将上述替换变量代入后,原数学模型的目标函数(6)中的网损成本以及网络约束条件(8)—(11)可作如下转变:

\({{C}^{\text{powerloss}}}=\sum\limits_{t=1}^{T}{[{{C}^{\text{loss}}}\sum\limits_{i,j\in {{\Omega }_{\text{ N}}}}^{{}}{{{\overset{\text{ }\!\!\tilde{\ }\!\!\text{ }}{\mathop{I}}\,}_{ij,t}}{{r}_{ij}}}]}\) (19)

$\begin{align} {{P}_{\text{ess},i,t}}+{{P}_{\text{wind},i,t}}+{{P}_{\text{grid},i,t}}-{{P}_{\text{load},i,t}}\text{=} \\ \sum\limits_{j\in c(i)}{({{X}_{i,t}}{{g}_{ij}}-{{Y}_{ij,t}}{{g}_{ij}}-{{Z}_{ij,t}}{{b}_{ij}})} \\\end{align}$ (20)

$\begin{align} {{Q}_{\text{ess},i,t}}+{{Q}_{\text{grid},i,t}}-{{Q}_{\text{load},i,t}}\text{=} \\ \sum\limits_{j\in c(i)}{(-{{X}_{i,t}}{{b}_{ij}}+{{Y}_{ij,t}}{{b}_{ij}}-{{Z}_{ij,t}}{{g}_{ij}})} \\\end{align}$ (21)

$U{{_{i}^{\min }}^{2}}\le {{X}_{i,t}}\le U{{_{i}^{\max }}^{2}}$ (22)

\({{\overset{\text{ }\!\!\tilde{\ }\!\!\text{ }}{\mathop{I}}\,}_{ij,t}}=({{X}_{i,t}}+{{X}_{j,t}}-2{{Y}_{ij,t}})(g_{ij}^{2}+b_{ij}^{2})\le I{{_{ij}^{\max }}^{2}}\) (23)

此外,替换变量Xi,tYij,tZij,t还需满足自身约束式(24)。至此原模型已经转变为线性方程式和二次等式的结合。为满足凸规划要求,将二次等式(24)进一步松弛并转化为标准二阶锥形式,如式(25)所示。文献[23]指出模型中的松弛是严格的,且根据文献[24-26]所述要求,本文含网损成本的模型严格符合转化条件。

${{X}_{i,t}}{{X}_{j,t}}=Y_{ij,t}^{2}+Z_{ij,t}^{2}$ (24)

2.1.2 储能约束转化

储能约束中式(14)—(17)为线性约束,符合二阶锥规划要求。非线性约束式(12)(13)可经如下转化变为锥约束及线性约束的结合,其中含0-1变量的非线性表达式可用Big-M法进行处理。

式中:M1M2M3为较大的常数;H1H2H3为新增辅助变量。

εi=1时,模型中换流器额定容量、储能充放电功率上限值分别取到Srate ess,i、Pmax,charge ess,i和Pmax,discharge ess,i,由于M为较大常数,约束(26)(27)将不会影响H的取值。当εi=0时,约束(26)(27)分别将换流器额定容量、储能充放电功率上限值约束至0。

经过上述处理,原混合整数非凸非线性规划问题转化为混合整数二阶锥规划问题。简化表达如下:

式中:x表示二阶锥模型中的所有优化变量,包括Pgrid,i,tQgrid,i,tPess,i,tQess,i,tEess,iεiXi,tYij,tZij,t等;uwind=[umin, umax],umaxumin分别为风电出力上下限,且不确定参数u受可调鲁棒参数Γ的约束;xk对应于优化变量中的替换变量以及储能功率变量,此两者均受到二阶锥约束,二阶锥简化模型kni如下所示:

$\sum\limits_{j=1}^{{{n}_{i}}-1}{x_{ij}^{2}}-x_{i0}^{2}\le 0,{{x}_{i0}}\ge 0$ (29)

2.2 基于拉格朗日对偶的单层模型转化
2.2.1 鲁棒模型对偶转化

根据文献[27],二阶锥模型(28)中的目标函数等价于下式:

$\underset{u}{\mathop{\max }}\,\underset{x}{\mathop{\min }}\,{{c}^{\text{T}}}x$ (30)

再根据锥优化的强对偶理论[28],原二阶锥模型可转化为

式中: α是一个自由变量;δ为对偶二阶锥,其推导如式(32)所示。其中αβτ为对偶过程中引入的辅助变量,无实际意义。

\(\begin{align} & \tau (\sum\limits_{j=1}^{{{n}_{i}}-1}{x_{ij}^{2}-x_{i0}^{2}})\text{=}-\text{ }\!\![\!\!\text{ (}\tau {{x}_{i0}}\text{)}{{x}_{i0}}+\sum\limits_{j=1}^{{{n}_{i}}-1}{(-\tau {{x}_{ij}}){{x}_{ij}}}]= \\ & \text{ }-(\tau {{x}_{i0}},-\tau {{x}_{i1}},...,-\tau {{x}_{i({{n}_{i}}-1)}}){{x}_{i}}\text{=}-\delta x \\ \end{align}\) (32)

通过上述处理方法,原鲁棒规划问题转变为单层双线性规划问题。根据文献[20],对于式(31)中存在的双线性项uTα可用Big-M法进行处理如下:

式中:M是一个较大的常数;$\widetilde{\alpha }$为辅助变量;μ为0-1变量。

由上式可见:当μi=1时,${{\widetilde{\alpha }}_{i}}$=αi,此时风电出力取到最大值,其对应的αi为正;当μi=0时,${{\widetilde{\alpha }}_{i}}$=0,此时风电出力取到最小值,其对应的αi为负。

2.2.2 可调参数的代入

为防止结果过于保守,根据式(5)代入可调鲁棒参数Γ,改变不确定参数取值如下:

式中:μiμ- i用于对各时间断面的不确定区间进行调整,为简化起见,均定为0-1变量;Γ为整个时段内的鲁棒总调节参数。

因此,模型(33)进一步转化如下:

通过上述二阶锥转化及鲁棒对偶转化,原问题已成为一个易求解的单层混合整数凸规划问题,可利用成熟的商业软件包进行计算。

3 算例分析

3.1 算例介绍

为验证本文提出的可调鲁棒规划模型的正确性,在MATLAB R2014 b下进行建模,并利用算法包CPLEX12.6.3进行计算,系统硬件环境为Intel Core I5 CPU,3.30 GHz,8 GB内存,操作系统为Win10 64bit。

本文利用改进的IEEE-33节点标准算例进行验证分析,算例拓扑图及风机接入节点见附录图A1。IEEE-33节点系统最大有功负荷为7.43 MW,负荷波动见附录图A2,基准电压为12.66 KV,除根节点电压恒定外,其他节点最大电压偏差范围为±10%。风机额定容量均为0.5 MW,总渗透率约为40%。风机的切入风速为5 m/s, 切出风速为25 m/s, 额定风速为15 m/s。鉴于风机往往在夜间出力较大,本文取风速波动及不确定区间如附录图A3所示。

为保证模型的精确性,网损电价取为0.1 USD/ (kW•h)。分时电价参考文献[29]。单个蓄电池储能单元的参数如附录表A1所示。

3.2 算例结果及分析比较
3.2.1 可调鲁棒模型分析

1)规划结果经济性分析。

风速的预测误差取为±20%,不确定区间总调节度Γ取为12,连接储能的节点总数取为3。规划结果如表1所示。从表1中的规划结果可以看出,由于考虑了风电出力不确定性,可调鲁棒规划中的储能总投资比确定性规划更多。为了对两种规划结果进行比较。利用蒙特卡洛法随机生成500个场景,对鲁棒规划结果进行经济性分析。运行成本结果对比见表2。

表1 系统规划结果 Tab. 1 Planning results

表2 可调鲁棒规划方案与传统规划方案对比 Tab. 2 Comparison for operation costs between adjustable robust method and traditional method

表中总运行成本的最大值为总成本最多场景的取值,并非各项成本最大值的简单累加。从表中可以看出:①在考虑不确定性的情况下,系统中节点2处的储能容量配置有一定的增加,说明节点2处配置储能相比其他节点具有更强的应对风电出力不确定性的能力。②可调鲁棒规划不仅能够减少主网购电成本,还能有效降低网损。这是由于更多的储能投入使得配电网对电能具有更强的时空平移能力,尽可能选择在电价较低时段进行购电,并且储能降低等效负荷所带来网损减少量大于增加等效负荷所引起的网损增量,从而使总体网损成本有所降低。③当风电出力随机波动时,可调鲁棒规划方案的每日总运行成本均值和最大值均小于确定性规划方案。系统平均每日总运行成本降低约210.93 USD,在不计折现率情况下6~7年即可收回额外增加的投资成本。综上,本文可调鲁棒规划模型相比确定性规划模型具有更高的经济价值。

2)储能系统对配电网运行成本影响。

在已知储能系统选址定容方案的基础上,进一步分析储能系统对配电网运行成本的影响。考虑两种储能系统选址定容方案,一种是采用本文可调鲁棒规划得到的方案另一种是采用确定性规划得到的方案,在风电出力最小场景下,以主网购电成本与网损成本之和最小为目标,确定两种方案下储能系统的分时充放电功率及向主网购买的电量。

图1是两种方案对应的主网购电量,储能的加入使得配电网在满足负荷供给的情况下,在实时电价相对较低的时刻购电进行储存,在实时电价较高的时刻大量减少购电,从而降低总购电成本。采用可调鲁棒规划得到的储能配置方案,比采用确定性规划得到的储能配置方案,在用电高峰期时的购电量更少,相应地配电网的购电成本也减少了。

图1 系统主网供电曲线 Fig. 1 Active power supplied by the main network

图2是可调鲁棒规划得到的储能配置方案对应的储能系统分时充放电功率及SOC。在2:00—8:00,储能系统主要处于充电状态,是由于这一时间段实时电价较低,适当增加购电进行储存以便于应对之后的用电高峰期,实现电能的时空转移。用电高峰时段中,储能系统在10:00—16:00处于放电状态,对电力系统进行削峰。在16:00—20:00这一时段由于电价相对较低,储能系统又转至充电状态使系统在下一个实时电价较高时段减少购电。

图2 储能充放电曲线以及SOC曲线 Fig. 2 ESS charge (discharge) power curve and SOC curve

3)不同预测误差下规划结果比较。

为了分析不确定区间范围对规划结果造成的影响,分别对风速预测误差为0、10%、20%、30%、40%、50%、60%时进行规划,最终结果见表3。

表3可见,随着预测误差的增大,可调鲁棒规划下的储能投资成本及各项运行成本均不断升高。当风速预测误差达到40%以后,储能投资成本

表3 不同预测误差下的规划结果 Tab. 3 Comparison for the planning results under different prediction error

的增量开始减少,这是由于较多时间点的风速波动区间下界开始低于切入风速,使得风电出力下端点为0,无法继续减少,从而导致预测误差的增大对储能投资成本的影响逐渐减小。

4)鲁棒参数的调整。

为体现可调鲁棒模型在经济性与鲁棒性上的统一,将可调参数调整至最大值,得出最恶劣场景下的规划结果,即鲁棒规划结果。表4给出了Γ=12的可调鲁棒规划结果与鲁棒规划结果的经济性比对。表中增加的总投资成本及每日减少的运行成本均是相对于确定性规划结果的。鲁棒规划结果比可调鲁棒规划结果,虽然每日减少的运行成本更多,但增加的总投资成本也更多;若仅以每年减少的运行成本来回收增加的总投资成本,鲁棒规划结果约需7.1年,而可调鲁棒规划只需约6.7年。风电出力在预测区间内随机波动时,鲁棒规划模型在经济性上稍劣于总调节度Γ为12的可调鲁棒规划模型。

表4 规划结果经济性对比 Tab. 4 Economic comparison of planning results

3.2.2 算法正确性分析

1)松弛精确性验证。

为了验证式(25)中二阶锥松弛的精确性,即最优解处的替换变量能否满足等式要求,定义绝对误差指标如下:

${{\Delta }_{ij,t}}\text{=}{{X}_{i,t}}{{X}_{j,t}}-{{Y}_{ij,t}}^{2}-{{Z}_{ij,t}}^{2}$ (36)

以某极限场景(风电出力最小场景)为例计算各线路绝对误差指标。误差散点图如图3所示。

从图中可以看出,松弛后的绝对误差数量级为10-4,且有大部分为10-5,满足规划要求[30]

图3 松弛误差散点图 Fig. 3 Scatter for relaxation deviation

2)锥模型强对偶验证。

为了解决双层问题,文中用到了二阶锥规划的对偶转化。考虑到模型的复杂性,现以预测风速场景以及两个极端场景为例,对可调鲁棒模型对偶转化的精确性进行验证。原对偶问题的优化目标及相对误差总结如下,其中优化目标值为投资运行综合成本。

表5 对偶误差 Tab. 5 Deviation of dualization

由上表可见,对偶转化后模型的解与原模型的解相差不大误差较小,能满足规划需求。

4 结论

本文构建了基于风电全消纳的配电网储能系统可调鲁棒规划模型,然后进行适当处理使其转化为易于求解的等价模型。最终通过算例仿真得出以下结论:

1)由于考虑了风电出力的不确定性,鲁棒规划中储能的容量配置会适当提高,从而使系统在实际运行过程中风电出力出现不确定波动时,能有效降低主网购电成本以及网损成本。可调参数的引入使得规划人员能够根据实际情况设定总调节参数,从而获得更具经济性的规划方案。

2)对于配电网,储能系统可在低电价时段存储电能留待高电价时段使用,能够减少高电价时段向主网的购电量,降低运行成本。

3)储能的配置容量与风速预测误差有关,预测误差越大,系统的储能需求越大。

4)通过对模型进行二阶锥松弛以及拉格朗日对偶,使得原问题转化为易于求解的单层凸规划问题,且转化误差能够满足规划要求。

附录A

表A1 单个储能单元参数 Tab. A1 Parameters of a single ESS unit

图A1 IEEE-33节点配电网系统 Fig. A1 Structure of IEEE 33-node distribution system

图A2 总有功负荷波动 Fig. A2 Active Load demand

图A3 风速波动区间 Fig. A3 Wind speed fluctuation interval

附录见本刊网络版(http://www.dwjs.com.cn/CN/volumn/current.shtml)。

参考文献

[1] Yuan Y,Zhang X S,Ju P,et al.Determination of economic dispatch of wind farm-battery energy storage system using genetic algorithm[J].International Transactions on Electrical Energy Systems,2014,24(2):264-280.

[2] 王成山,于波,肖峻,.平滑可再生能源发电系统输出波动的储能系统容量优化方法[J].中国电机工程学报,2012,32(16):1-8. Wang Chengshan,Yu Bo,Xiao Jun,et al.Hybrid energy storage system and control system design for wind power balancing[J].Proceedings of the CSEE,2012,32(16):1-8(in Chinese).

[3] 胡永强,刘晨亮,赵书强,.基于模糊相关机会规划的储能优化控制[J].电力系统自动化,2014,38(6):20-25. Hu Yongqiang,Liu Chenliang,Zhao Shuqiang,et al.Optimal control of energy storage based on fuzzy correlated-chance programming[J].Automation of Electric Power Systems,2014,38(6):20-25(in Chinese).

[4] Jiang Q,Gong Y,Wang H.A battery energy storage system dual-layer control strategy for mitigating windfarm fluctuations[J].IEEE Transactions on Power Systems,2013,28(3):3263-3273.

[5] 杨水丽,李建林,惠东,.用于跟踪风电场计划出力的电池储能系统容量优化配置[J].电网技术,2014,38(6):1485-1491. Yang Shuili,Li Jianlin,Hui Dong,et al.Optimal capacity configuration of battery energy storage system to track planned output of wind farm[J].Power System Technology,2014,38(6):1485-1491(in Chinese).

[6] Hartmann B,Dan A.Cooperation of a grid-connected wind farm and an energy storage unit-demonstration of a simulation tool[J].IEEE Trans on Sustainable Energy,2012,3(1):49-56.

[7] 兑潇玮,朱桂萍,刘艳章.考虑预测误差的风电场储能配置优化方法[J].电网技术,2017,41(2):434-439. Dui Xiaowei,Zhu Guiping,Liu Yanzhang.Research on battery storage sizing for wind farm considering forecast error[J].Power System Technology,2017,41(2):434-439(in Chinese).

[8] 赵书强,王扬,徐岩,.基于机会约束目标规划的高风电接入比例下大规模储能与火电协调调度[J].中国电机工程学报,2016,36(4):969-977. Zhao Shuqiang,Wang Yang,Xu Yan,et al.Coordinated dispatch of large scale energy storage system and thermal generation in high wind power penetration level system based on chance constrained goal programming[J].Proceedings of the CSEE,2016,36(4):969-977(in Chinese).

[9] 黎静华,文劲宇,程时杰,.基于p_有效点理论的含大规模风电电力系统最小储能功率配置方法[J].中国电机工程学报,2013,33(13):45-52. Li Jinghua,Wen Jinyu,Cheng Shijie,et al.Method of minimum energy storage power allocation for electric power systems with large-scale wind power based on p-efficient point theory[J].Proceedings of the CSEE,2013,33(13):45-52(in Chinese).

[10] 郑乐,胡伟,陆秋瑜,.储能系统用于提高风电接入的规划和运行综合优化模型[J].中国电机工程学报,2014,34(16):2533-2542. Zheng Le,Hu Wei,Lu Qiuyu,et al.Research on planning and operation model for energy storage system to optimize wind power integration[J].Proceedings of the CSEE,2014,34(16):2533-2542(in Chinese).

[11] Dutta S,Sharma R.Optimal storage sizing for integrating wind and load forecast uncertainties[C]//Innovative Smart Grid Technologies (ISGT).IEEE,2012:1-7.

[12] Zhang N,Kang C Q,Kirschen D S,et al.Planning pumped storage capacity for wind power integration[J].IEEE Transactions on Sustainable Energy,2013,4(2):393-401.

[13] 崔杨,朱加明,葛延峰,.计及限风特性与经济性的大规模集中充放电储能系统容量配置研究[J].电网技术,2016,40(2):484-490. Cui Yang,Zhu Jiaming,Ge Yanfeng,et al.Study on capacity configuration of large-scale energy storage system under concentrated charge and discharge based on wind curtailment characteristics and economics[J].Power System Technology,2016,40(2):484-490(in Chinese).

[14] Abbey C,Joós G.A stochastic optimization approach to rating of energy storage systems in wind-diesel isolated grids[J].IEEE Transactions on Power Systems,2009,24(1):418-426.

[15] 张沈习,程浩忠,邢海军,.配电网中考虑不确定性的分布式电源规划研究综述[J].电力自动化设备,2016,36(8):1-9. Zhang Shenxi,Cheng Haozhong,Xing Haijun,et al.Review of DG planning considering uncertainties for distribution network[J].Electric Power Automation Equipment,2016,36(8):1-9(in Chinese).

[16] 韩杏宁,黎嘉明,文劲宇,.含多风电场的电力系统储能鲁棒优化配置方法[J].中国电机工程学报,2015,35(9):2120-2127. Han Xingning,Li Jiaming,Wen Jinyu,et al.Optimization for robust energy storage allocation in power system with multiple wind farms integrated[J].Proceedings of the CSEE,2015,35(9):2120-2127(in Chinese).

[17] Ding Tao,Liu Shiyu,Yuan Wei,et al.A two-stage robust reactive power optimization considering uncertain wind power integration in active distribution networks[J].IEEE Transactions on Sustainable Energy,2016,7(1):301-311.

[18] Peng Xiong,Panida Jirutitijaroen,Chanan Singh.A distributionally robust optimization model for unit commitment considering uncertain wind power generation[J].IEEE Transactions on Power Systems,2017,32(1):39-49.

[19] 向月,刘俊勇,魏震波,.考虑可再生能源出力不确定性的微电网能量优化鲁棒模型[J].中国电机工程学报,2014,34(19):3063-3072. Xiang Yue,Liu Junyong,Wei Zhenbo,et al.Robust model of microgrid energy optimization with uncertain renewable energy sources[J].Proceedings of the CSEE,2014,34(19):3063-3072(in Chinese).

[20] 温俊强,曾博,张建华,.配电网中分布式风电可调鲁棒优化规划[J].电网技术,2016,40(1):227-233. Wen Junqiang,Zeng Bo,Zhang Jianhua,et al.Adjustable robust optimization for distributed wind power planning in distribution network[J].Power System Technology,2016,40(1):227-233(in Chinese).

[21] 赵金利,于莹莹,李鹏,.基于锥优化的储能系统参与配电网运行调节快速计算方法[J].电力系统自动化,2016,40(2):30-35. Zhao Jinli,Yu Yingying,Li Peng,et al.A fast calculation method of energy system in distribution network based on conic programming[J].Automation of Electric Power Systems,2016,40(2):30-35(in Chinese).

[22] 李洪美,崔翰韬,万秋兰.考虑电动汽车充电策略的配网重构二阶锥规划模型[J].中国电机工程学报,2015,35(18):4674-4681. Li Hongmei,Cui Hantao,Wan Qiulan.Distribution network reconfiguration based on second-order conic programming considering EV charging strategy[J].Proceedings of the CSEE,2015,35(18):4674-4681(in Chinese).

[23] Jabr R A.Radial distribution load flow using conic programming[J].IEEE Transactions on Power Systems,2006,21(3):1458-1459.

[24] Farivar M,Low S H.Branch flow model:relaxations and convexification (part I)[J].IEEE Transactions on Power Systems,2013,28(3):2554-2564.

[25] Taylor J A,Hover F S.Convex models of distribution system reconfiguration[J].IEEE Transactions on Power Systems,2012,27(3):1407-1413.

[26] Li N,Chen L,Low S H.Exact convex relaxation of OPF for radial networks using branch flow model[C]//2012 Smart Grid Communications(Smart Grid Comm).Tainan:IEEE,2012:7-12.

[27] Beck A,Bental A.Duality in robust optimization:primal worst equals dual best[J].Operations Research Letters,2009,37(1):1-6.

[28] 李静. 锥规划及其对偶锥规划的若干性质及应用[D].武汉:武汉理工大学,2005:20-33.

[29] 李秀磊,耿光飞,季玉琦,.主动配电网中储能和需求侧响应的联合优化规划[J].电网技术,2016,40(12):3803-3810. Li Xiulei,Geng Guangfei,Ji Yuqi,et al.Integrated optimal planning of energy storage and demand side response in active power distribution network[J].Power System Technology,2016,40(12):3803-3810(in Chinese).

[30] 高红均,刘俊勇,魏震波.主动配电网分层鲁棒规划模型及其求解方法[J].中国电机工程学报,2017,37(5):1389-1400. Gao Hongjun,Liu Junyong,Wei Zhenbo.A bi-level robust planning model of active distribution network and its solution method[J].Proceedings of the CSEE,2017,37(5):1389-1400(in Chinese).

  • 目录

    图1