CN102545214B - 一种含电池储能设备的风电场可靠性评估方法 - Google Patents

Abstract

本发明提供一种含电池储能设备的风电场可靠性评估方法：采集风电场的原始数据、建立风电机组功率时序模型、建立含电池储能设备的输出功率序列模型、建立含电池储能设备的风电场可靠性评估模型、建立计及运行特性参数的可靠性评估模型、将原始采集的数据输入计算机系统出处理，获得计及电池储能设备运行特性参数的风电场可靠性参数；在本发明中，考虑了不同类型电池储能设备运行特性参数对风电场可靠性的影响，提高了对含电池储能设备的风电场可靠性评估的准确性。

Description

一种含电池储能设备的风电场可靠性评估方法

技术领域

本发明涉及风电场可靠性评估方法，具体涉及一种计及电池储能设备运行特性的风电场可靠性评估方法，属于电力工程技术领域。 

背景技术

风能具有间歇性和随机性的特点，其输出功率的波动性是限制大规模风电场接入电网的重要因素之一。为了抑制波动、平滑风电场输出功率，常将储能设备接入以提高系统稳定性和可靠性，电池储能设备有功率型和能量型两种。 

现有技术中，针对电池储能设备风电场可靠性评估有记下三种方法：(1)计及电池储能的风电场可靠性概率评估模型与离散的风速威布尔分布结合，提出系统可靠性指标的计算分析方法(Reliability Evaluation of a Wind-Diesel Hybrid Power System with Battery Bank Using Discrete Wind Speed Frame Analysis[C].International Conference on Probabilistic Methods Applied to Power System，Stockholm，Sweden，2006)，该方法未考虑储能自身特性对系统运行状态的限制，提出的模型较为简单；(2)使用解析法评估包含风电机组和储能设备的发电系统可靠性(Reliability assessment of a wind-power system with integrated energy storage[J].Renewable Power Generation,IET，2010，4(3)：211-219，)，该方法初步分析了储能设备对风电场可靠性的提升效果，计及了储能设备自身故障率以 及容量限制、充放电速率等因素的影响，但未分析不同类型储能设备充放电过程差异对系统可靠性的影响；(3)计算包含储能设备大型风电系统可靠性指标的通用方法和模型是建立在前面两种方法的基础上的(Reliability evaluation of generating systems containing wind power and energy storage[J].IET Generation,Transmission&Distribution，2009，3(8)：783-79)，该方法在宏观上计及了储能设备不同能量存储释放控制策略对系统可靠性造成的差异，但仍未考虑设备自身充放电特性对可靠性的影响，不同类型储能设备充放电过程差异。 

综上所述，已有分风电场可靠性评估的方法和模型中，没有考虑不同类型储能设备在运行状态中各项参数对风电场可靠性的影响，做出的可靠性分析准确性不高。 

发明内容

针对现有技术存在的上述问题，本发明提供一种计及电池储能设备运行特性的风电场可靠性评估方法，能够对含电池储能设备的风电场可靠性做出准确地评估。 

本发明是这样实现的： 

一种含电池储能设备的风电场可靠性评估方法，具体步骤包括： 

步骤1:采集风电场原始数据，确定ARMA模型中风电场历史风速的每小时平均值μt和标准差σ，风电机组的额定功率、切入风速、额定风速、切出风速，负荷有功功率、切负荷功率；风电场中各风电机组、电池储能设备和开关设备之间的连接关系； 

步骤2：建立风电机组功率时序模型 

ARMA(n,m)模型数学表达式为 

y_t＝φ₁y_t-1+φ₂y_t-2+…+φ_ny_t-n+α_t-θ₁α_t-1-θ₂α_t-2-…-θ_mα_t-m

                                           (1) 

式中，y_t为时刻t的风速序列值，Φ_i(i＝1,2,…,n)和θ_j(j＝1,2,…,m)分别为ARMA模型的自回归参数和滑动平均参数； 

利用ARMA(n,m)模型获得风电场在任意时刻t的风速时间序列模拟值后，通过历史风速的每小时平均值μ_t和标准差σ_tσ_t求取时刻t的风电场风速SW_t： 

SW_t＝μ_t+σ_ty_t   (2) 

利用(3)式计算得到任意t时刻下风电机组的输出功率： 

$P\left({\mathrm{SW}}_t\right)=\begin{array}{c}\left\{\begin{array}{cc}0,& 0\≤{\mathrm{SW}}_t\≤V_{\mathrm{ci}}\\ (A+B\×{\mathrm{SW}}_t+C\×& \\ {\mathrm{SW}}_t^2)\×P_r,& V_{\mathrm{ci}}\≤{\mathrm{SW}}_t\≤V_r\\ P_r,& V_r\≤{\mathrm{SW}}_t\≤V_{\mathrm{co}}\\ 0,& {\mathrm{SW}}_t\≥V_{\mathrm{co}}\end{array}\right.\end{array}---\left(3\right)$

式中，P_r、V_ci、V_r和V_co分别为风电机组的额定功率、切入风速、额定风速和切出风速，参数A、B、C可由V_ci、V_r和V_co进行计算； 

步骤3：建立含电池储能设备的输出功率序列模型 

风电场输出功率高于当前负荷的某一比例P％，常规发电机组输出功率低于当前负荷的1-P％，储能设备开始释放能量参与供电，基于此种策略建立储能设备充放电模型： 

1)利用(1)-(3)式获取在某一时间段T内的风电机组功率输出时间序列W_i,； 

2)根据时间段T内的常规发电机组功率时间序列{G_i,i＝1,2,.......T}以及利用负荷曲线得到的负荷时间序列{L_i,i＝1,2,......T}，P％表示风电场输出功率低于或高于当前负荷的比例，分别利用(4)式和(5)式求取风电机组剩余功率序列{P_Wi,i＝1,2,......T}和常规机组剩余功率序列{P_Gi,i＝1,2,......T}： 

P_Wi＝W_i-P％×L_i   (4) 

P_Gi＝G_i-(1-P％)×L_i   (5) 

3)求取电池储能设备的输出功率序列{P_Bi,i＝1,2,......T}： 

$P_{\mathrm{Bi}+1}=\left\{\begin{array}{cc}{P}_{\mathrm{Bi}}+P_{\mathrm{Wi}},& P_{\mathrm{Wi}}\&GreaterEqual;0,P_{\mathrm{Gi}}\&GreaterEqual;0\\ P_{\mathrm{Bi}}+P_{\mathrm{Gi}},& P_{\mathrm{Wi}}\&GreaterEqual;0,P_{\mathrm{Gi}}<0\\ P_{\mathrm{Bi}},& P_{\mathrm{Wi}}<0,(P_{\mathrm{Wi}}+P_{\mathrm{Gi}})\&GreaterEqual;0\\ P_{\mathrm{Bi}}+P_{\mathrm{Wi}}+P_{\mathrm{Gi}},& P_{\mathrm{Wi}}<0,(P_{\mathrm{Wi}}+P_{\mathrm{Gi}})<0\end{array}\right.---\left(6\right)$

步骤4：建立含电池储能设备的风电场可靠性评估模型 

将风电机组和储能设备看作负荷，对原始负荷曲线与输出或输入功率进行叠加，得到等效负荷曲线： 

L_Ni＝L_i+L_Di-ΣP_Wi-ΣP_Bi   (7) 

其中，P_Wi、P_Gi、P_Bi、L_i和L_Di分别为T时刻下的风电机组剩余功率、常规机组剩余功率、电池储能设备输出功率、负荷有功功率以及切负荷功率； 

由此可得到每一抽样年中发生的缺电次数L_LTi、缺电持续时间L_LDi和缺供电量ENS_i，对于N年的抽样时间，缺电时间期望LOLE电量不足期望LOEE的计算式可由(9)和(10)得到： 

ENS_i＝ΣL_LDiL_Ni   (8) 

$\mathrm{LOLE}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{L}_{\mathrm{LDi}}---\left(9\right)$

$\mathrm{LOLE}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{\mathrm{ENS}}_i---\left(10\right)$

步骤5：计及电池储能设备运行特性参数的可靠性评估模型 

计及电池储能设备运行参数设备容量、充放电速率和充放电电压稳定性的影响下，风电场可靠性评估模型为： 

1)在一个时间段内对常规机组进行抽样，得到每小时的机组功率序列P_Gi； 

2)结合当前负荷大小，确定电池储能系统的初始充放电功率P_ini

$P_{\mathrm{ini}}=\left\{\begin{array}{cc}-P_{\mathrm{disch}-\max }& P_{\mathrm{Gi}}-P_{\mathrm{Li}}\&le;-P_{\mathrm{disch}-\max }\\ P_{\mathrm{Gi}}-P_{\mathrm{Li}}& -P_{\mathrm{disch}-\max }\&le;P_{\mathrm{Gi}}-P_{\mathrm{Li}}\&le;P_{\mathrm{ch}-\max }\\ P_{\mathrm{ch}-\max }& P_{\mathrm{Gi}}-P_{\mathrm{Li}}>P_{\mathrm{ch}-\max }\end{array}\right.---\left(11\right)$

其中，P_disch-max和P_ch-max分别为设备最大放电功率和设备最大充电功率； 

3)当前时段内设备每小时释放或充入的能量，能量型和功率型设备的E_Bi计算表达式分别如(12)和(13)所示： 

$E_{\mathrm{Bi}+1}=\left\{\begin{array}{cc}{E}_{B\min }& E_{\mathrm{Bi}}+T\&times;P_{\mathrm{ini}}\&le;E_{B\min }\\ E_{\mathrm{Bi}}+T\&times;P_{\mathrm{ini}}& E_{B\min }\&le;E_{\mathrm{Bi}}+T\&times;P_{\mathrm{ini}}\&le;E_{B\max }\\ E_{B\max }& E_{\mathrm{Bi}}+T_{\mathrm{Pini}}>E_{B\max }\end{array}\right.---\left(12\right)$

$E_{\mathrm{Bi}+1}=\left\{\begin{array}{cc}{E}_{B\min }& E_{\mathrm{Bi}}+E_N\&le;E_{B\min }\\ E_{\mathrm{Bi}}+E_N& E_{B\min }\&le;E_{\mathrm{Bi}}+E_N\&le;E_{B\max }\\ B_{\max }& E_{\mathrm{Bi}}+E_N>E_{B\max }\end{array}\right.---\left(13\right)$

其中， 

$E_N={\&Integral;}_0^T{f}_P\left(t\right)\mathrm{dt}---\left(14\right)$

由于E_N与充放电初始功率P_ini有关，因此f_p(t)实际上是功率型电池设备的功率时间函数曲线簇，即： 

f_p(t)＝f_pb(t)+P_ini,0≤t≤T   (15) 

4)获取时间段T内设备的可用充放电功率P_avai，可用放电功率P_avai-out和可用充电功率P_avai-in，计算公式如下： 

$P_{\mathrm{avai}-\mathrm{out}}=\left\{\begin{array}{cc}(E_{\mathrm{Bi}}-E_{B\min })/T& (E_{\mathrm{Bi}}-E_{B\min })/T\&le;P_{\mathrm{disch}-\max }\\ P_{\mathrm{disch}-\max }& (E_{\mathrm{Bi}}-E_{B\min })/T>P_{\mathrm{disch}-\max }\end{array}\right.---\left(16\right)$

$P_{\mathrm{avai}-\mathrm{in}}=\left\{\begin{array}{cc}(E_{B\max }-E_{\mathrm{Bi}})/T& (E_{B\max }-E_{\mathrm{Bi}})/T\&le;P_{\mathrm{ch}-\max }\\ P_{\mathrm{ch}-\max }& (E_{B\max }-E_{\mathrm{Bi}})/T>P_{\mathrm{ch}-\max }\end{array}\right.---\left(17\right)$

步骤6：计算计及电池储能设备运行特性参数的风电场可靠性参数 

根据(6)式电池储能设备的输出功率序列以及时间段T内的可用充放电功率，结合时间段T内设备的可用充放电功率P_avai，获取在n个模拟年内储能设备的完整充放电功率序列{P_Bt(i)，i＝1,2,…,8736×n}，根据式(4)、(5)和(6)分别求得P_Wi、P_Gi、和P_Bi、以及利用式(7)表示的等效负荷曲线，使用序贯蒙特卡洛法求取相应的可靠性指标。 

附图说明

图1-含风电机组及电池储能设备的发电系统模型。 

图2-两类设备在风电场VL下的可靠性指标。 

图3-两类设备在风电场WI下的可靠性指标。 

图4-LOLE指标在不同设备容量下的变化。 

图5-LOLE指标在不同放电限制下的变化。 

具体实施方式

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

参见图1，一种含电池储能设备的风电场可靠性评估方法，首先采集风电场的原始数据、然后建立风电机组功率时序模型和含电池储能设备的输出功率序列模型、再建立含电池储能设备的风电场可靠性评估模型和建立计及运行特性参数的可靠性评估模型、最后将采集的数据输入计算机系统出处理，获得计及电池储能设备运行特性参数的风电场可靠性参数；具体步骤包括： 

步骤1:采集风电场原始数据，确定ARMA模型中风电场历史风速的每小时平均值μ_t和标准差σ，风电机组的额定功率、切入风速、额定风速、切出风速，负荷有功功率、切负荷功率；风电场中各风电机组、电池储能设备和开关设备之间的连接关系； 

步骤2：建立风电机组功率时序模型 

ARMA(n,m)模型数学表达式为 

y_t＝φ₁y_t-1+φ₂y_t-2+…+φ_ny_t-n+α_t-θ₁α_t-1-θ₂α_t-2-…-θ_mα_t-m

                                           (1) 

式中，y_t为时刻t的风速序列值，Φ_i(i＝1,2,…,n)和θ_j(j＝1,2,…,m) 分别为ARMA模型的自回归参数和滑动平均参数； 

利用ARMA(n,m)模型获得风电场在任意时刻t的风速时间序列模拟值后，通过历史风速的每小时平均值μ_t和标准差σ_tσ_t求取时刻t的风电场风速SW_t： 

SW_t＝μ_t+σ_ty_t   (2) 

利用(3)式计算得到任意t时刻下风电机组的输出功率： 

$P\left({\mathrm{SW}}_t\right)=\begin{array}{c}\left\{\begin{array}{cc}0,& 0\&le;{\mathrm{SW}}_t\&le;V_{\mathrm{ci}}\\ (A+B\&times;{\mathrm{SW}}_t+C\&times;& \\ {\mathrm{SW}}_t^2)\&times;P_r,& V_{\mathrm{ci}}\&le;{\mathrm{SW}}_t\&le;V_r\\ P_r,& V_r\&le;{\mathrm{SW}}_t\&le;V_{\mathrm{co}}\\ 0,& {\mathrm{SW}}_t\&GreaterEqual;V_{\mathrm{co}}\end{array}\right.\end{array}---\left(3\right)$

式中，P_r、V_ci、V_r和V_co分别为风电机组的额定功率、切入风速、额定风速和切出风速，参数A、B、C可由V_ci、V_r和V_co进行计算； 

步骤3：建立含电池储能设备的输出功率序列模型 

风电场输出功率高于当前负荷的某一比例P％，常规发电机组输出功率低于当前负荷的1-P％，储能设备开始释放能量参与供电，基于此种策略建立储能设备充放电模型： 

1)利用(1)-(3)式获取在某一时间段T内的风电机组功率输出时间序列W_i,； 

2)根据时间段T内的常规发电机组功率时间序列{G_i,i＝1,2,.......T}以及利用负荷曲线得到的负荷时间序列{L_i,i＝1,2,......T}，P％表示风电场输出功率低于或高于当前负荷的比例，分别利用(4)式和(5)式求取风电机组剩余功 率序列{P_Wi,i＝1,2,......T}和常规机组剩余功率序列{P_Gi,i＝1,2,......T}： 

P_Wi＝W_i-P％×L_i   (4) 

P_Gi＝G_i-(1-P％)×L_i   (5) 

3)求取电池储能设备的输出功率序列{P_Bi,i＝1,2,......T}： 

$P_{\mathrm{Bi}+1}=\left\{\begin{array}{cc}{P}_{\mathrm{Bi}}+P_{\mathrm{Wi}},& P_{\mathrm{Wi}}\&GreaterEqual;0,P_{\mathrm{Gi}}\&GreaterEqual;0\\ P_{\mathrm{Bi}}+P_{\mathrm{Gi}},& P_{\mathrm{Wi}}\&GreaterEqual;0,P_{\mathrm{Gi}}<0\\ P_{\mathrm{Bi}},& P_{\mathrm{Wi}}<0,(P_{\mathrm{Wi}}+P_{\mathrm{Gi}})\&GreaterEqual;0\\ P_{\mathrm{Bi}}+P_{\mathrm{Wi}}+P_{\mathrm{Gi}},& P_{\mathrm{Wi}}<0,(P_{\mathrm{Wi}}+P_{\mathrm{Gi}})<0\end{array}\right.---\left(6\right)$

步骤4：建立含电池储能设备的风电场可靠性评估模型 

将风电机组和储能设备看作负荷，对原始负荷曲线与输出或输入功率进行叠加，得到等效负荷曲线： 

L_Ni＝L_i+L_Di-ΣP_Wi-ΣP_Bi   (7) 

其中，P_Wi、P_Gi、P_Bi、L_i和L_Di分别为T时刻下的风电机组剩余功率、常规机组剩余功率、电池储能设备输出功率、负荷有功功率以及切负荷功率； 

由此可得到每一抽样年中发生的缺电次数L_LTi、缺电持续时间L_LDi和缺供电量ENSi，对于N年的抽样时间，缺电时间期望LOLE电量不足期望LOEE的计算式可由(9)和(10)得到： 

ENS_i＝ΣL_LDiL_Ni   (8) 

$\mathrm{LOLE}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{L}_{\mathrm{LDi}}---\left(9\right)$

$\mathrm{LOLE}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{\mathrm{ENS}}_i---\left(10\right)$

步骤5：计及电池储能设备运行特性参数的可靠性评估模型 

计及电池储能设备运行参数设备容量、充放电速率和充放电电压稳定性的影响下，风电场可靠性评估模型为： 

1)在一个时间段内对常规机组进行抽样，得到每小时的机组功率序列P_Gi； 

2)结合当前负荷大小，确定电池储能系统的初始充放电功率P_ini

$P_{\mathrm{ini}}=\left\{\begin{array}{cc}-P_{\mathrm{disch}-\max }& P_{\mathrm{Gi}}-P_{\mathrm{Li}}\&le;-P_{\mathrm{disch}-\max }\\ P_{\mathrm{Gi}}-P_{\mathrm{Li}}& -P_{\mathrm{disch}-\max }\&le;P_{\mathrm{Gi}}-P_{\mathrm{Li}}\&le;P_{\mathrm{ch}-\max }\\ P_{\mathrm{ch}-\max }& P_{\mathrm{Gi}}-P_{\mathrm{Li}}>P_{\mathrm{ch}-\max }\end{array}\right.---\left(11\right)$

其中，P_disch-max和P_ch-max分别为设备最大放电功率和设备最大充电功率； 

3)当前时段内设备每小时释放或充入的能量，能量型和功率型设备的E_Bi计算表达式分别如(12)和(13)所示： 

$E_{\mathrm{Bi}+1}=\left\{\begin{array}{cc}{E}_{B\min }& E_{\mathrm{Bi}}+T\&times;P_{\mathrm{ini}}\&le;E_{B\min }\\ E_{\mathrm{Bi}}+T\&times;P_{\mathrm{ini}}& E_{B\min }\&le;E_{\mathrm{Bi}}+T\&times;P_{\mathrm{ini}}\&le;E_{B\max }\\ E_{B\max }& E_{\mathrm{Bi}}+T_{\mathrm{Pini}}>E_{B\max }\end{array}\right.---\left(12\right)$

$E_{\mathrm{Bi}+1}=\left\{\begin{array}{cc}{E}_{B\min }& E_{\mathrm{Bi}}+E_N\&le;E_{B\min }\\ E_{\mathrm{Bi}}+E_N& E_{B\min }\&le;E_{\mathrm{Bi}}+E_N\&le;E_{B\max }\\ E_{B\max }& E_{\mathrm{Bi}}+E_N>E_{B\max }\end{array}\right.---\left(13\right)$

其中， 

$E_N={\&Integral;}_0^T{f}_P\left(t\right)\mathrm{dt}---\left(14\right)$

由于E_N与充放电初始功率P_ini有关，因此f_p(t)实际上是功率型电池设备的功率时间函数曲线簇，即： 

f_p(t)＝f_pb(t)+P_ini,0≤t≤T   (15) 

4)获取时间段T内设备的可用充放电功率P_avai，可用放电功率P_avai-out和可用充电功率P_avai-in，计算公式如下： 

$P_{\mathrm{avai}-\mathrm{out}}=\left\{\begin{array}{cc}(E_{\mathrm{Bi}}-E_{B\min })/T& (E_{\mathrm{Bi}}-E_{B\min })/T\&le;P_{\mathrm{disch}-\max }\\ P_{\mathrm{disch}-\max }& (E_{\mathrm{Bi}}-E_{B\min })/T>P_{\mathrm{disch}-\max }\end{array}\right.---\left(16\right)$

$P_{\mathrm{avai}-\mathrm{in}}=\left\{\begin{array}{cc}(E_{B\max }-E_{\mathrm{Bi}})/T& (E_{B\max }-E_{\mathrm{Bi}})/T\&le;P_{\mathrm{ch}-\max }\\ P_{\mathrm{ch}-\max }& (E_{B\max }-E_{\mathrm{Bi}})/T>P_{\mathrm{ch}-\max }\end{array}\right.---\left(17\right)$

步骤6：计算计及电池储能设备运行特性参数的风电场可靠性参数 

根据(6)式电池储能设备的输出功率序列以及时间段T内的可用充放电功率，结合时间段T内设备的可用充放电功率P_avai，获取在n个模拟年内储能设备的完整充放电功率序列{P_Bt(i)，i＝1,2,…,8736×n}，根据式(4)、(5)和(6)分别求得P_Wi、P_Gi、和P_Bi、以及利用式(7)表示的等效负荷曲线，使用序贯蒙特卡洛法求取相应的可靠性指标。 

实施例：使用RBTS系统进行分析，对其分别接入两个风电场作为待比较的两个测试系统，其风速数据来源于简称为VL的荷兰Vlieland和简称为WI的Wijdenes两个风电场过去10年的历史统计，风电场VL中，平均风速为5.1m/s，机组的V_ci、V_r和V_co分别为3.5m/s、13m/s、26m/s；风电场WI中，平均风速为4.3m/s，机组的V_ci、V_r和V_co分别为3m/s、15m/s、27.5m/s。 

两个风电场的风速ARMA模型分别如式(21)和式(22)所示： 

y_t＝3.482y_t-1-4.615y_t-2+2.77y_t-3+α_t-2.544α_t-1+2.238α_t-2

                                                  (21) 

α_t∈NID(0,0.31039²) 

y_t＝1.832y_t-1-0.7834y_t-2+0.0504y_t-3+α_t-0.9676α_t-1+0.0065α_t-2

                                                      (22) 

α_t∈NID(0,0.37593²) 

RBTS系统总装机容量为240MW，其中5MW机组2台，10MW机组1台，20MW机组5台，40MW机组3台。两个风电场各自的装机总容量均为18.0MW，渗透率P％为9.73％，风电机组单机容量为2MW，每个风电场中的机组数量为9台。 

选用的功率型电池储能设备和能量型电池储能设备总容量均为10MWh，对于功率型设备，通过统计数据进行二次拟合后得到的功率时间函数曲线簇表达式为： 

f_p(t)＝0.3873t²-4.7775t+12.7202+P_ini   (23) 

风速数据每1小时取一次，因此电池设备充放电时间段T为1小时。 

(1)两类电池储能设备对可靠性影响的比较 

电池充电设备输出功率序列模型的建立决定于电池储能设备的充放策略，本申请计及电池储能设备运行特性，在步骤3建立含电池储能设备的输出功率序列模型时，使用的策略为：风电场输出功率高于当前负荷的某一比例P％，常规发电机组输出功率低于当前负荷的1-P％，储能设备开始释放能量参与供电。 

现有技术中计及电池储能设备运行特性的策略有：策略2：风电场输出功率低于当前负荷的某一比例P％，储能设备就开始释放能量参与供电；策略3：风电场输出功率高于当前负荷的某一比例P％，常规发电机组输出功率低于当前负荷的1-P％，储能设备开始释放能量参与供电。 

(2)基于策略2和策略3建立含电池储能设备的输出功率序列模型 

1)利用(1)-(3)式获取在某一时间段T内的风电机组功率输出时间序列W_i,； 

2)根据此时间段内的常规发电机组功率时间序列{G_i,i＝1,2,...T}以及利用负荷曲线得到的负荷时间序列{L_i,i＝1,2,...T}，P％表示风电场输出功率低于或高于当前负荷的比例，分别利用(24)式和(25)式求取风电机组剩余功率序列{P_Wi,i＝1,2,...T}和常规机组剩余功率序列{P_Gi,i＝1,2,...T}： 

P_Wi＝W_i-P％×L_i   (24) 

P_Gi＝G_i-(1-P％)×L_i   (25) 

3)求取电池储能设备的输出功率序列{P_Bi,i＝1,2,...T}： 

策略2：P_Bi+1＝P_Bi+P_Wi

(26) 

策略3： 

$P_{\mathrm{Bi}+1}=\left\{\begin{array}{cc}{P}_{\mathrm{Bi}}+P_{\mathrm{Wi}},& P_{\mathrm{Wi}}\&GreaterEqual;0,P_{\mathrm{Gi}}\&GreaterEqual;0\\ P_{\mathrm{Bi}}+P_{\mathrm{Gi}},& P_{\mathrm{Wi}}\&GreaterEqual;0,P_{\mathrm{Gi}}<0\\ P_{\mathrm{Bi}}+P_{\mathrm{Wi}},& P_{\mathrm{Wi}}<0\end{array}\right.---\left(27\right)$

(3)在三种储能策略下，分别将能量型设备和功率型设备接入RBTS系统中，计算得到的可靠性指标以及各指标。 

1)两类设备在风电场VL中的可靠性指标如表1和图2所示 

表1.风电场VL中不同设备和储能策略的可靠性指标 

由图表比较分析可知，在三种充放电策略下，两类设备对风电场可靠性均有所改善。但能量型电池储能设备对系统可靠性的改善优于功率型电池储能设备，功率型设备对可靠性改善的程度不如能量型设备。事实上，功率型设备等效后的功率与能量型设备稳定的功率水平相比，其差距并不大，因此图表中反映的两类指标数据也十分接近。 

2)两类设备在风电场WI中的可靠性指标如表2和图3所示 

表2.接入风电场WI时不同设备和储能策略的可靠性指标 

对于不同的充放电策略，其可靠性指标优秀程度排序为：本申请策略>策略2>策略3。从较长的时间段内来看，电池设备处于满能量运行的时间较长，可从外界接受的能量一直处于较低的水平，一旦出现即将失负荷的情况，设备内的能量将用于填补这部分负荷缺额，因此对可靠性的提升最为明显； 策略2和策略3主要将储能设备作为平衡风电场渗透率的手段，目的在于平滑风电场的输出功率曲线，浪费了储能设备在负荷端的功率弥补能力，所以对系统可靠性的提升均不如本申请策略；而策略3面临削减负荷的概率更大，因此该策略的表现不如策略2。 

(4)储能设备容量对可靠性的影响 

参照图4可知，随着储能设备容量的增加，不论是何种储能策略，可靠性指标均得以不断提升，其中，在本申请策略下，可靠性指标提升较为明显，而另两种策略的提升则不大，尤其是策略3，当设备容量增加到30MW后，可靠性指标已几乎没有变化。 

(5)储能设备充放电速率对可靠性的影响 

参照图5可知，两类设备的P_disch-max越大，可靠性指标越好，但增加到一定水平后，指标的变化幅度并不大。事实上，放电限制的数值越高，设备在制造工艺上实现的难度越大，造价也越昂贵，因此在选择相应的电池储能设备时，应根据风电场容量、系统规模等实际情况进行优化选择。 

Claims (1)

1.一种含电池储能设备的风电场可靠性评估方法，其特征在于：具体步骤包括：

步骤1:采集风电场原始数据，确定 ARMA（n，m）模型中风电场历史风速的每小时平均值μ_t和标准差σ_t ，风电机组的额定功率、切入风速、额定风速、切出风速，负荷有功功率、切负荷功率以及风电场中各风电机组、电池储能设备和开关设备之间的连接关系；

步骤2：建立风电机组功率时序模型

ARMA(n,m)模型数学表达式为

 （1）

式中，为时刻t的风速序列值，Φ_i和θ_j分别为ARMA模型的自回归参数和滑动平均参数，其中，i=1, 2,…, n；j=1, 2,…, m；

利用ARMA(n,m)模型获得风电场在任意时刻t的风速时间序列模拟值后，通过历史风速的每小时平均值μ_t和标准差σ_t求取时刻t的风电场风速SW_t：

           （2）

利用（3）式计算得到任意t时刻下风电机组的输出功率：

   （3）

式中，P_r、V_ci、V_r和V_co分别为风电机组的额定功率、切入风速、额定风速和切出风速，参数A、B、C由V_ci、V_r和V_co进行计算；

步骤3：建立含电池储能设备的输出功率序列模型

风电场输出功率高于当前负荷的某一比例P%，常规发电机组输出功率低于当前负荷的1-P%，储能设备开始释放能量参与供电，基于此种策略建立储能设备充放电模型：

1)利用(1)-(3)式获取在某一时间段T内的风电机组功率输出时间序列W_i；

2)根据时间段T内的常规发电机组功率时间序列以及利用负荷曲线得到的负荷时间序列，P%表示风电场输出功率高于当前负荷的比例，分别利用(4)式和(5)式求取风电机组剩余功率序列和常规机组剩余功率序列：

           (4)

        (5)

3)求取电池储能设备的输出功率序列：

(6)

是指第i+1个时段电池储能设备的输出功率； 

步骤4：建立含电池储能设备的风电场可靠性评估模型

将风电机组和储能设备看作负荷，对原始负荷曲线与输出或输入功率进行叠加，得到等效负荷曲线：

     (7)

其中，P_Wi、P_Gi、P_Bi、L_i和L_Di分别为T时刻下的风电机组剩余功率、常规机组剩余功率、电池储能设备输出功率、负荷有功功率以及切负荷功率； 

由此可得到每一抽样年中发生的缺电次数L_LTi、缺电持续时间L_LDi和缺供电量ENS_i，对于N年的抽样时间，缺电时间期望LOLE和电量不足期望LOEE的计算式由(9)和(10)得到：

         (8)

        (9)

        (10)

步骤5：计及电池储能设备运行特性参数的可靠性评估模型

计及电池储能设备运行参数设备容量、充放电速率和充放电电压稳定性的影响下，风电场可靠性评估模型为：

1）在一个时间段内对常规机组进行抽样，得到常规机组剩余功率序列； 

2）结合当前负荷大小，确定电池储能系统的初始充放电功率P_ini

 (11)

其中，P_disch-max和P_ch-max分别为设备最大放电功率和设备最大充电功率；

3）能量型和功率型设备的当前时段内设备每小时释放或充入的能量E_Bi计算表达式分别如(12)和(13)所示：(12)

其中，指的是能量型和功率型设备第i+1时段每小时释放或充入的能量；

    (13)

其中，

          (14)

由于E_N与充放电初始功率P_ini有关，因此f_p(t)实际上是功率型电池设备的功率时间函数曲线簇，即：

    (15)

其中，通过统计数据进行二次拟合获取； 

4）获取时间段T内设备的可用充放电功率P_avai，可用放电功率P_avai-out和可用充电功率P_avai-in，计算公式如下：

 (16)

 （17）

步骤6：计算计及电池储能设备运行特性参数的风电场可靠性参数根据(6)式电池储能设备的输出功率序列以及时间段T内的可用充放电功率，结合时间段T内设备的可用充放电功率P_avai，获取在n个模拟年内储能设备的完整充放电功率序列{P_Bt(i)，i=1, 2, …, 8736×n}，根据式（4）、（5）和（6）分别求得P_Wi、P_Gi、和P_Bi、以及利用式（7）表示的等效负荷曲线，使用序贯蒙特卡洛法求取相应的可靠性指标。