一种并网风光发电中混合储能容量优化配置方法

一种并网风光发电中混合储能容量优化配置方法
【专利摘要】本发明公开了一种并网风光发电中混合储能容量优化配置方法，确定以优化蓄电池工作状态为原则，以提高储能系统整体经济性为目标的能量管理策略，基于该能量管理策略分析了并网风光发电系统能量损失率和能量缺失率的计算流程，根据全生命周期费用理论，建立了储能装置的年均费用函数表达式，并建立了以该函数值最小为目标，以系统能量损失率及能量缺失率等运行指标为约束的储能容量优化配置模型，最后运用改进混沌优化算法求解优化配置模型。本发明提出的改进混沌优化算法利用混沌运动具有的遍历性、随机性、规律性，可以有效完成该复杂非线性优化配置模型的计算。
【专利说明】
一种并网风光发电中混合储能容量优化配置方法
技术领域
[0001] 本发明涉及一种并网风光发电中混合储能容量优化配置方法，属于微网储能系统 配置优化技术领域。
【背景技术】
[0002] 随着风光发电规模的迅速扩大，并网型风光发电系统将成为一个极具应用前景的 发展方向。然而，风光发电功率具有的不稳定性和间歇性使其不能像传统发电方式那样满 足电网要求的调度灵活性，为解决这一问题需要开发和建设与之相配套的储能系统。
[0003] 储能技术应用到风光发电系统中将使系统的能量管理更加灵活，但对于目前以蓄 电池为常用储能元件的风光储项目，储能成本接近甚至大于购买单独的光伏发电组件或风 力发电机的成本。因此如何最大限度地延长蓄电池组的使用寿命是新能源发电系统设计人 员亟须解决的重大问题。在目前可应用的储能技术中，蓄电池和超级电容器的储能特性可 以达到很好的优势互补，因此该类混合储能技术将在风光系统中得到广泛应用。

【发明内容】

[0004] 本发明所要解决的技术问题是克服现有技术的缺陷，提供一种并网风光发电中混 合储能容量优化配置方法，以蓄电池和超级电容器构成的混合储能系统为基础，对并网风 光发电中储能容量优化配置问题进行研究。
[0005] 为解决上述技术问题，本发明提供一种并网风光发电中混合储能容量优化配置方 法，具体包括如下步骤：
[0006] 1)确定能量管理策略；
[0007] 2)计算并网风光发电系统运行指标；
[0008] 3)以混合储能装置的年均费用最小为目标函数；
[0009] 4)确定约束条件；
[0010] 5)采用改进混沌优化算法求解优化配置模型。
[0011] 前述的步骤1)确定能量管理策略，是指超级电容器承担由发电功率与调度功率不 相等产生的不平衡功率Pub中的频繁波动功率;蓄电池只负责承担不平衡功率P ub中的基本 功率。
[0012] 前述的步骤2)中，采用能量损失率和能量缺失率作为并网风光发电系统运行指 标;所述能量损失率Rlppp是指在系统运行过程中，能量损失量Elpp与电源发电总量E g的比值； 所述能量缺失率Rlpsp是指在系统运行过程中，能量缺失量Elps与调度总需求量En的比值;计 算公式如下：
[0013]
(!)
[0014] 在kt时刻，若光伏阵列输出功率Ppv(kt)和风机输出功之和大于电网调度 功率Pn(kt)，则存在剩余功率卩1(^，此时，如果储能系统剩余能量E。nr未处于储能系统能量最 大值Ecnmax，则电源为储能系统充电；若储能系统剩余能量处于满额状态，将产生能量损失， 即：
[0015] ELPp(kt) = [(Ppv(kt)+Pw(kt)) · nc-Pn(kt)]t,
[0016] 其中，Ew(kt)表示kt时刻的能量损失量;η。表示逆变器的转换效率;t为单位时间 间隔；
[0017] 当光伏阵列和风机输出功率之和小于电网调度功率要求时，出现缺失功率plac;k， 若储能系统剩余能量Ec^大于储能系统能量最小值Ε_ ιη，则储能系统进行放电;若储能系统 剩余能量不大于储能系统能量最小值，则由电网来弥补能量缺失：
[0018] ELPs(kt) = {Pn(kt)-[Ppv(kt)+Pw(kt) ]nc}t (2)
[0019] 其中，ELPs(kt)表示kt时刻的能量缺失量。
[0020] 前述的步骤3)中，以混合储能装置的年均费用最小为目标函数，具体表达式如下：
[0021]
(14)
[0022]其中，C表示混合储能装置的年均费用，Clv指混合储能装置的购置成本，fba、f uc分 别代表蓄电池组和超级电容器组的单价，m、η分别表示蓄电池组和超级电容器组的个数， kdeba、kdeuj别表示蓄电池组和超级电容器组的年折旧值，C。。指混合储能装置的运行成本， 分别表示蓄电池和超级电容器的运行成本系数，C mc指混合储能装置的维护成本， 1^^表示蓄电池的维护成本系数，Cdc指混合储能装置的处置成本，kdcb a、kdcuc分别表示蓄电 池和超级电容器的处置成本系数。
[0023]前述的步骤4)中，约束条件如下：
[0024]
(15)
[0025] 其中，Pba(kt，m)、Puc(kt，n)分别表示蓄电池和超级电容器的充电功率，m、n分别表 示蓄电池组和超级电容器组的个数;？0(1^)^1^)、？11(1^)表示1^时刻光伏阵列输出功 率、风机输出功率和电网调度功率;表示逆变器的转换效率;Pban(HI)、Eban(m)分别表示蓄 电池的额定功率和额定容量;α表示不平衡功率P ub= 中基本 功率所占比例;RLPPP(m，n)、RLPSP(m，n)分别表不能量损失率和能量缺失率，RLPPPmax、RLPSPmax分 别表示能量损失率最大值和能量缺失率最大值;Pu?ax(n)表示超级电容器最大充电功率； Ebamin(m)表示蓄电池的最小容量;Eba(kt，m)表示蓄电池经过充放电后的剩余能量;E_i n (!^^^^"!^分别表示超级电容器组能量的最小值和最大值^^仏扒"表示超级电容器组 经过充放电后的剩余能量。
[0026] 前述的步骤5)采用改进混沌优化算法求解优化配置模型，具体步骤如下：
[0027] (5-1)算法初始化，令k = 0,给出N个具有微小差异的初值Xi,〇,且Xi,〇辛1且Xi,〇辛 0.75且Xi, ο辛0.5且Xi, ο辛0.25且Xi, ο辛0,然后按照Logi sti c映射表达式得至IjN个轨迹不同的 混纯变量{Xi,k};
[0028] (5-2)将生成的混沌变量{Xl,k}按下式进行第一次载波，得到优化变量范围[x min， Xmax]内的混纯变量W i,k}:
[0029]
(17)；
[0030] (5-3)对生成的混沌变量Yi, !^进行可行性检验，重复执行步骤(5-1)和步骤(5-2) 进行迭代搜索得到满足约束条件的可行解{Vik}，并计算相应的目标函数值f(k) = f ({x'k});
[0031]
;否贝 IJ，令 k = k+ I，进入下一循环，直到产变化小于0.01或达到最大迭代次数Nmaxl;
[0032] (5-4)按下式对混沌变量Xi,k+1进行第二次载波，得到混沌变量{V i,k}:
[0033] Z7 i,k = x*+Pxi,k+i (18)
[0034] 其中，β为调节因子；
[0035] (5-5)对混沌变量乂 i,行可行性检验，当其满足约束条件时，若f (乂 i,k) <f'则 = (z\,k);否则令k = k+l，进入下一循环，直到产变化小于0.01或达到最大迭 代次数Nmax2，最终的/即为最优解，产为最优值。
[0036]前述的步骤(5-1)中，Logistic映射表达式为：
[0037] xi,k+i = yxi,k(l-xi,k) ,Xi.kG [0,1] (16)
[0038] 其中，μ为控制参数;Xi,k为混纯变量。
[0039]本发明所达到的有益效果：
[0040] 相较于传统方法，在容量配置优化计算中，本发明充分考虑蓄电池和超级电容器 的互补运行特性，以优化蓄电池工作状态为原则提出的能量管理策略能够有效提高储能系 统的整体经济性；以全生命周期费用理论为基础建立的储能容量优化目标函数，比只考虑 初始购置成本的优化目标更符合实际情况，具有现实意义。本发明提出的改进混沌优化算 法利用混沌运动具有的遍历性、随机性、规律性，可以有效完成该复杂非线性优化配置模型 的计算。
【附图说明】
[0041] 图1是本发明并网风光发电中混合储能容量优化配置方法流程图；
[0042]图2是能量损失率和能量缺失率计算流程图；
[0043]图3是储能装置充电时能量损失量的计算流程图；
[0044]图4是储能装置放电时能量缺失量的计算流程图；
[0045] 图5是实施例中风力发电功率、光伏功率和调度功率0.5h发电预测值；
[0046] 图6是实施例中单一储能时蓄电池组输出功率曲线；
[0047] 图7是实施例中混合储能时蓄电池组和超级电容器组输出功率曲线。
【具体实施方式】
[0048] 下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明 的技术方案，而不能以此来限制本发明的保护范围。
[0049] 本发明的并网风光发电中混合储能容量优化配置方法，如图1所示，包括如下步 骤：
[0050] 步骤一:确定能量管理策略；
[0051] 步骤二:计算并网风光发电系统运行指标；
[0052] 步骤三:以混合储能装置的年均费用最小为目标函数；
[0053]步骤四：确定约束条件；
[0054]步骤五:采用改进混沌优化算法求解优化配置模型。
[0055] 具体实施过程如下：
[0056]以优化蓄电池工作状态为原则，以提高储能系统整体经济性为目标的能量管理策 略。基于该能量管理策略分析了并网风光发电系统能量损失率和能量缺失率的计算流程， 根据全生命周期费用理论，建立了储能装置的年均费用函数表达式，并建立了以该函数值 最小为目标，以系统能量损失率及能量缺失率等运行指标为约束的储能容量优化配置模 型，最后运用改进混沌优化算法求解优化配置模型。下文以某国家示范工程一期建设项目 为例进行分析计算，验证所提方法的技术合理性和经济实用性。
[0057] 1)确定能量管理策略
[0058]考虑到间歇式电源功率的波动性，并结合蓄电池和超级电容器互补的储能特性， 提出如下适合该混合储能系统的能量管理策略:为充分发挥超级电容器功率密度大、循环 寿命长、响应速度快的优势，同时有效回避其能量密度低的缺陷，超级电容器主要承担由发 电功率与调度功率不相等产生的不平衡功率P ub中的频繁波动功率;蓄电池只负责承担不平 衡功率Pub中的基本功率，这样可以有效避免蓄电池频繁充放电造成的小循环充放电现象， 且蓄电池均以额定功率充放电，保证蓄电池处于优化的工作状态，延长蓄电池的使用寿命， 提高储能系统的整体经济性。
[0059] 2)计算并网风光发电系统运行指标
[0060] 本发明采用能量损失率(LPPP)和能量缺失率(LPSP)作为并网风光发电系统运行 指标。定义能量损失率Rlppp为在系统运行过程中，能量损失量Elpp与电源发电总量E g的比值； 定义能量缺失率Rlpsp为在系统运行过程中，能量缺失量Elps与调度总需求量En的比值。其满 足：
[0061] Cl)
[0062]基于上述所提出的能量管理策略，能量损失率Ruw和能量缺失率Rlpsp的计算流程 如图2所示。图中：t为单位时间间隔；T为总计算次数;Pw(kt)、Ppv(ktWPP n(ktb^m*kt 时刻风机输出功率、光伏阵列输出功率和电网调度功率;Pirft和Pid分别表示剩余功率和缺 失功率;E cnr、EbaJPEucr*别为储能系统、蓄电池组和超级电容器组充放电之前的剩余能量； Eba和Euc表示蓄电池组和超级电容器组经过充放电后的剩余能量;Ec_x、E cnmin分别表示储能 系统能量最大值和最小值;η。表示逆变器的转换效率。
[0063] 在kt时刻，若光伏阵列输出功率Ppv(kt)和风机输出功率?,(!^)之和大于电网调度 功率P n(kt)，则存在剩余功率?1(^，此时，如果储能系统剩余能量Ecn^处于满额状态E cnmax， 则电源为储能系统充电;若储能系统剩余能量处于满额状态，将产生能量损失，
[0064] gpELPP(kt) = [(Ppv(kt)+Pw(kt)) · nc_Pn(kt)]t，
[0065] 其中，ELPP(kt)表示kt时刻的能量损失量;Ppv(kt)、Pw(kt)、P n(kt)表示kt时刻风 机、光伏阵列输出功率和电网调度功率;η。表示逆变器的转换效率。当光伏阵列输出功率和 风机输出功率之和小于电网调度功率要求时，出现功率缺失Pid，若储能系统剩余能量Emr 大于最小剩余能量Emmin储能系统进行放电，以弥补功率缺失。若储能系统剩余能量不大于 最小剩余能量，则由电网来弥补能量缺失：
[0066] ELPs(kt) = {Pn(kt)-[Ppv(kt)+Pw(kt) ]nc}t (2)
[0067] 其中：ELPs(kt)表示kt时刻的能量缺失量;Ppv(kt)、Pw(kt)、P n(kt)表示kt时刻风 机、光伏阵列输出功率和电网调度功率;η。表示逆变器的转换效率。
[0068]从图3中，根据不同的充电情况，能量损失量Elpp的计算流程共有7种情况，即S&~ SC7,在每一种充电情况下储能装置的充电功率和剩余能量情况、系统的能量损失量计算过 程如下：
[0069] (I)SC1。系统的剩余功率大到足以使蓄电池以额定功率Pban充电，超级电容器以最 大功率Pucmax充电，而且蓄电池经过该单位时间段的充电达到了额定容量Eban。在这种情况下 的能量损失量Elpp包括直接废弃的电源功率、充电效率产生的损失以及因蓄电池满额产生 的浪费，具体计算如下式所示：
[0070] ElPP 一 [ Plef t/π。- ( Pban+Pucmax )]t+[Pban(l -Hbac ) +Pucmax ( I -Hucc )] t+ ( Eba-Eban ) ( 3 )
[0071] 式中：Ρ?Α=[ρρν(?α)+ρν(?α)]η「Ρ η(?α)表示剩余功率;n。表示逆变器的转换效 率;Pucmax表示超级电容器的最大充电功率;t为单位时间间隔;Hb ac、η%。分别为蓄电池和超级 电容器的充电效率;Eba表示蓄电池经过充放电后的剩余能量。
[0072] (2) SC2。该状态类似于状态SC1，唯一不同的是:蓄电池并没有充电至额定容量Eban， 因此这种情况下能量损失量Elpp的计算如下式所示：
[0073] ElPP 一 [ Plef t/π。- ( Pban+Pucmax )]t+[Pban(l -Hbac ) +Pucmax ( I -Hucc )] t ( 4 )
[0074] 式中：Pleft表示剩余功率;η。表示逆变器的转换效率;pu_x表示超级电容器的最大 充电功率;t为单位时间间隔;ri bac；、rw。分别为蓄电池和超级电容器的充电效率。
[0075] (3)SC3。起初蓄电池以额定功率P-优先充电，超级电容器充电功率P uc = PirftAb-Pban。经过tl = (Eban-Ebar)/(Pbannba)c后，其中tl<t，蓄电池达到额定容量Eban，因此在之后的 七士时间段仅有超级电容器充电，且P uc = P_ax。这种情况下能量损失量Elpp的计算如下式所 示：
[0076] ElPP - Pban( l-TUiacOtl+l! (Pleft/nc_Pban)tl+Pucmax(t_tl) ] ( I-flucc)+(Pleft/nc_Pucmax) (t_tl) ( 5 )
[0077] 式中：Plrft表示剩余功率;n。表示逆变器的转换效率;Pirnmax表示超级电容器的最大 充电功率;t为单位时间间隔;ri bac；、rw。分别为蓄电池和超级电容器的充电效率。
[0078] (4)SC4。该状态与SC3类似，唯一不同的是:在t-ti时间段，超级电容器充电功率P uc = PirftAb。这种情况下能量损失量Elpp的计算如下式所示：
[0079] Elpp = Pban( I-Tlbac )tl+[ (Pleft Alc-Pban) tl+Pleft Alc (t-tl) ]( 1-riucc) (6)
[0080] 式中：？化11表示蓄电池额定功率;Pieft表示剩余功率;n。表示逆变器的转换效率； Hbac、Hu。。分别为蓄电池和超级电容器的充电效率。
[00811 (5) SC5。整个t时间段内，蓄电池的充电功率Pba = Pban、超级电容器的充电功率Pu。= PirftAb-Pban。可见该状态比SC3和SC4简单得多，原因在于蓄电池经过时间t的充电并未达到 额定容量。这种情况下系统能量损失量Elpp的计算如下式所示：
[0082] Elpp= [Pban( l_rlbac) + (Pleft/rlc_Pban) ( I-rIucc) ]t (7)
[0083] 式中：？^"表示蓄电池额定功率;Pleft表示剩余功率;η。表示逆变器的转换效率； nbac、rwj别为蓄电池和超级电容器的充电效率;t为单位时间间隔。
[0084] (6) SC6。系统的剩余功率小于蓄电池额定功率，因此只有超级电容器充电，且充电 功率为Pirc = Puanax能量损失量Elpp的计算如下式所示：
[0085] Elpp= (Pief t/nc-Pu cmax )t+P u cmax (l-n〇c〇)]t (8)
[0086] 式中：Pleft表示剩余功率;n。表示逆变器的转换效率;pu_x表示超级电容器的最大 充电功率;rw。为超级电容器的充电效率;t为单位时间间隔。
[0087] (7) SC7。系统的剩余功率小于蓄电池额定功率，因此只有超级电容器充电，且充电 功率为puc; = ple3ft/ri。，能量损失量Elpp的计算如下式所示：
[0088] Elpp = PieftAc (l-nucc)t (9)
[0089] 式中：Pleft表示剩余功率;η。表示逆变器的转换效率;rw。为超级电容器的充电效 率;t为单位时间间隔。
[0090] 类似于储能装置充电时对系统能量损失量Em的计算分析，放电过程中系统能量 缺失量Em的计算流程如图4所示。从图4中可以看出，根据不同的放电情况Em的计算共有5 种情况，即SDi~SD 5,在每一种放电状态下储能装置的放电功率和剩余能量、系统的能量缺 失量的计算如下：
[0091] (I)SDjPSD2这两种状态相似。系统的缺失功率大于蓄电池和超级电容器所能提供 的最大放电功率之和P ban+Pu_x，因此蓄电池和超级电容器放电功率分别是Pbar^PP uraax。唯 一不同的是，在SD2中，经过单位时间t的放电后蓄电池未放电至最小容量Ebamin，而SDl却正 好相反。这两种状态的能量缺失量Ews的计算分别下式所示：
[0092] ElPS 一 [Plack- ( Pban-I-Pucmax ) He ] (Ebamin-Eba ) TlbacTlc (10)
[0093] Elps= [Plack-(Pban+Pucmax) He] t (11)
[0094] 式中：Piack=pn(kt)-[ppv(kt)+p w(kt)]nc表示缺失功率;η。表示逆变器的转换效 率;nba。为蓄电池的充电效率;t为单位时间间隔。
[0095] (2)SD3。系统缺失功率小于储能系统最大放电功率，但是大于蓄电池额定功率，因 此根据上述提出的能量管理策略，该缺失功率由蓄电池以额定功率P ban放电，超级电容器以 功率Puc = PlackAc-Pban放电，二者共同进行功率弥补。经过tl = (Ebar-Ebamin )/( Pbannbac )，其中 tl<t，蓄电池放电至最小容量Ebamin，之后的t-tl时间段仅有超级电容器放电，但当Plack彡 pu_xq。时，存在能量缺失，计算如下式所示：
[0096] Elps= (Piack-Pucnc) (t-ti) (12)
[0097] 式中:Plack表示缺失功率;η。表示逆变器的转换效率;t为单位时间间隔。
[0098] (3)SD4。该状态与SD3类似，系统缺失功率介于储能系统最大放电功率和蓄电池额 定功率之间，因此由蓄电池以额定功率Pban放电、超级电容器以功率pu。= plac;k/nc;-pban放电， 二者共同弥补功率缺失。相比SD3,该情况下蓄电池没有放电至最小容量，所以不存在两阶 段放电情况，也没有能量缺失。
[0099] (4) SD5。系统缺失功率小于蓄电池额定功率，因此这种情况下只有超级电容器放 电，在Piad^PucWl。情况下系统能量缺失Elps为：
[0100] Elps= (Piack-Pucnc )t (13)
[0101] 式中:Pkk表示缺失功率;η。表示逆变器的转换效率;t为单位时间间隔。
[0102] 3)以混合储能装置的年均费用最小为目标函数
[0103] 储能系统优化目标模型定义为四大成本之和，即购置成本、运行成本、维护成本以 及处置成本。为尽量简化目标函数，以免优化模型过于复杂，本发明引用不同的比例系数表 示相应的成本。鉴于计算的成本包含两种储能元件，因此为体现合理性，将目标函数折算为 储能系统的年均费用，折算系数用储能元件每年各自的折旧值表示。具体表达式如下：
[01 04] niinC - Civ+Coc+Cmc+Cdc-(1+kocba+kmcba+kdcba) · kdebanrfba+(l+kocuc+kdcuc)kdeucnfuc ( 14 )
[0105] 式中：C表示储能装置的年均费用，Clv指储能装置的购置成本，包括储能装置及其 附加设备的购置成本;f ba、fu。分别代表蓄电池组和超级电容器组的单价;m、n分别表示蓄电 池组和超级电容器组的个数;kd eba、kdeuc分别表示蓄电池组和超级电容器组的年折旧值;C。。 指储能装置的运行成本，包括储能装置的试验、安装、损耗、停运成本、人工费用等;kocba、 k_c分别表示蓄电池和超级电容器的运行成本系数;Cmc指储能装置的维护成本，包括故障 前的检修成本和故障后的维修成本 ;kmcba表示蓄电池的维护成本系数，超级电容器可免维 护，因此不需考虑它的维护成本;C d。指储能装置的处置成本，包括储能装置的报废成本和残 值;1^^、1^。11。分别表示蓄电池和超级电容器的处置成本系数。
[0106] 由于蓄电池和超级电容器的储能特性有很大区别，它们的年折旧值需按不同的方 法计算。（1)蓄电池的折旧进程可以通过其容量损失程度来表示。一些研究表明，蓄电池的 容量损失与放电深度1_呈线性关系，所以具体计算其折旧值时，可先将不同放电深度的充 放电循环次数Ndod折算为等效的1 DQD，1DQD=100%充放电循环次数Ne。折算公式可表示为Ne = NdodIDQD。若蓄电池进行IDQD = 100 %充放电循环寿命次数为Nba，则蓄电池的折旧值kdeba = Ne/ Nba。（2)超级电容器的折旧值可以直接用循环次数来计算。若超级电容器进行从Eucmax到 Eucmin的充放电循环次数为Nuc，则进行N次充放电循环后，它的折旧值可表示为kd_ = N/Nuc。
[0107] 上述优化目标模型克服了传统经济性优化问题中只注重降低初始购置成本的局 限性，将储能装置成本延伸到安装、运行、维护、失效报废等整个寿命期过程中，使优化结果 更符合实际。
[0108] 4)确定约束条件
[0109] 根据上述能量管理策略，并考虑并网风光发电系统的运行特性，建立满足并网风 光发电系统安全可靠运行的约束条件如下：
[0110]
.(15)
[0111] 式中：Pba、PUc^v别表示蓄电池和超级电容器的充电功率;m、n分别表示蓄电池组和 超级电容器组的个数;？以1^)^(1^)几(1^)表示1^时刻光伏阵列、风机输出功率和电网 调度功率;η。表示逆变器的转换效率;p ban、Eban分别表示蓄电池的额定功率和额定容量； EUCmin(n)、Euaiiax(n)分别表示超级电容器组能量的最小值和最大值;α表示不平衡功率P ub = Pn(kt)-[Ppv(kt)+Pw(kt) ]n。I中基本功率所占比例;RlPPP、RlPSP分别表示能量损失率和能量 缺失率。
[0112] 上述公式中：第1个公式属于系统合理性的要求;第2、3、4、5个公式是保证蓄电池 和超级电容器正常运行的基本约束，可以看出为延长蓄电池寿命，蓄电池功率值或为零，或 为额定功率;第6个公式体现了所提的能量管理策略，即蓄电池的功率输出维持在基本负荷 以内；第7个公式表示系统的能量损失率必须小于给定的最大值RLPPPm ax，以确保风光发电功 率利用率;第8个公式表示系统的能量缺失率必须小于给定的最大值RLPSPmax。
[0113] 5)采用改进混沌优化算法求解优化配置模型
[0114] 该优化方法直接采用混沌变量进行搜索，搜索顺序按照下式所示的Logistic映射 表达式产生的点列依次进行。
[0115] xi,k+i = yxi,k(l-xi,k) ,Xi.kG [0,1] (16)
[0116] 式中:μ为控制参数;Xi,k为混纯变量。
[0117] Logistic映射是[0，1 ]上的不可逆映射。当μ取值4.0时，系统处于混沌状态，任意 取初始点，可以得到在[0，1 ]范围内不重复的遍历所有点。由于混纯运动具有遍历性、随机 性、规律性，因而混沌优化算法较其他随机性的智能优化方法获得最优解的可能性更强，是 一种极有前途的优化手段。本发明所用混沌优化算法主要分为两个阶段：（1)在整个优化变 量范围内按Logistic映射表达式依次检验产生的各点是否满足约束条件，在筛选出的符合 要求的点中找出当前最优点；（2)达到一定迭代次数后，认为当前最优点已在实际最优点附 近，然后以当前最优点为初值，仍按照（1)的方法进行小范围搜索最终找到最优点。
[0118] 混沌优化方法的基本求解步骤如下：
[0119] (1)算法初始化。令k = 0,给出N个具有微小差异的初值Xi,〇,且Xi,〇辛1且Xi,〇辛0.75 且Xi, Q辛0 · 5且Xi, Q辛0 · 25且Xi, Q辛0，然后按照Logi Stic映射表达式得到N个轨迹不同的混纯 变量{Xi,k}。
[0120] (2)将生成的混沌变量{Xl,k}按下式进行第一次载波，得到优化变量范围[Xmin， Xmax]内的混纯变量{X7 i,k}。
[0121] X i, k - Xmin+( Xmax-Xmin) Xi, k ( 17 )
[0122] (3)对生成的混沌变量Y i, k，进行可行性检验，重复执行步骤（1)和步骤(2)进行迭 代搜索得到满足约束条件的可行解{x"^}，并计算相应的目标函数值。令 父*=&〃1,()}，产=汽0)，如果以1〇<产，那么)(*={^ 1,1{}，产=以1〇;否则，令1^ = 1^1，进入下一 循环，直至Ijf变化小于0.01或达到最大迭代次数Nmaxl。
[0123] (4)按下式对混沌变量Xl,k+1进行第二次载波，β为调节因子，一般取值在优化变量 所属区间长度(Xmax-x_)的1 %和10 %之间。对于所属区间长度较小的优化变量，辟交大;反 之较小。
[0124] Z7 i,k = x*+Pxi,k+i (18)
[0125] (5)对生成的混沌变量k，进行可行性检验，当其满足约束条件时，若k)< f*，则x* = z\,k，f* = f (z\,k);否则令k = k+l，进入下一循环，直到f*变化小于0.01或达到最 大迭代次数Nmax2,最终的/即为最优解，产为最优值。
[0126] 为了进一步说明本发明方法的准确性和可靠性，以某国家示范工程一期建设项目 为例进行分析计算，该建设项目包含IOOMff的风电系统和40MW的太阳能光伏发电系统。假设 电网调度功率额定值为96丽，逆变器效率为0.95，系统供电指标RLPPPma#RLPSPmax分别为0.3 和0.05。现在为其配置储能容量。详细参数如表1所示。
[0127] 本发明采取50个传统固定型蓄电池单体、100个超级电容器单体串并联，分别组成 蓄电池组和超级电容器组，根据图5所示的风机年发电预测数据、光伏年发电预测数据、年 调度功率预测数据，运用改进混沌优化算法对混合储能系统的配置容量进行了优化计算。
[0128] 表1储能系统参数
[0130]为了对比蓄电池-超级电容器混合储能方式与蓄电池单一储能方式的经济优越 性，采用上述给出的相同原始数据，利用改进混沌优化算法进行分析计算，优化计算结果如 表2所示。
[0132]通过表2的数据可以看出，在相同的系统运行指标约束下，采用蓄电池单一储能方 式下需要配置蓄电池组的数目为5002组，储能系统的年均最小费用为3003余万元，而采用 蓄电池和超级电容器混合储能时，需配置5691组蓄电池和9875组超级电容器，储能系统年 均最小费用仅977余万元，约是前一种储能方式的1/3。由此说明，基于本发明提出的能量管 理策略的混合储能方式能够充分发挥两种储能元件的互补优势，可以有效降低储能系统成 本。风力发电功率、光伏功率和调度功率〇.5h发电预测值见图5。根据表2优化结果数据，绘 制两种储能方式下储能装置在4300-4400h时间段的输出功率变化曲线，分别如图6和7所 示。图7中：P b^示蓄电池输出功率;Puc表示超级电容器输出功率。将图6和图7对比可看出， 由于超级电容器组的参与，混合储能方式下，图7中蓄电池组的输出功率变化比图6平坦得 多，原因在于峰值不平衡功率和频变不平衡功率主要由超级电容器组提供，这样可以有效 地优化蓄电池组的工作状态，延长其使用寿命，最终减少储能系统的运行维护等费用。 [0133]以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形 也应视为本发明的保护范围。
【主权项】
1. 一种并网风光发电中混合储能容量优化配置方法，其特征在于，具体包括如下步骤： 1) 确定能量管理策略； 2) 计算并网风光发电系统运行指标； 3) 以混合储能装置的年均费用最小为目标函数； 4) 确定约束条件； 5) 采用改进混沌优化算法求解优化配置模型。2. 根据权利要求1所述的一种并网风光发电中混合储能容量优化配置方法，其特征在 于:所述步骤1)确定能量管理策略，是指超级电容器承担由发电功率与调度功率不相等产 生的不平衡功率P ub中的频繁波动功率;蓄电池只负责承担不平衡功率Pub中的基本功率。3. 根据权利要求1所述的一种并网风光发电中混合储能容量优化配置方法，其特征在 于:所述步骤2)中，采用能量损失率和能量缺失率作为并网风光发电系统运行指标;所述能 量损失率Rlppp是指在系统运行过程中，能量损失量Elpp与电源发电总量Eg的比值;所述能量 缺失率Rlpsp是指在系统运行过程中，能量缺失量Elps与调度总需求量En的比值;计算公式如 下：在kt时刻，若光伏阵列输出功率Ppv(kt)和风机输出功率?|(1^)之和大于电网调度功率 Pn(kt)，则存在剩余功率卩1(^，此时，如果储能系统剩余能量未处于储能系统能量最大值 Ecnmax，则电源为储能系统充电;若储能系统剩余能量处于满额状态，将产生能量损失，即： ELPp(kt) = [(Ppv(kt)+Pw(kt)) · nc-Pn(kt)]t, 其中，ELPP(kt)表示kt时刻的能量损失量;η。表示逆变器的转换效率;t为单位时间间隔； 当光伏阵列和风机输出功率之和小于电网调度功率要求时，出现缺失功率Pkk，若储能 系统剩余能量大于储能系统能量最小值Ε_ιη，则储能系统进行放电;若储能系统剩余能 量不大于储能系统能量最小值，则由电网来弥补能量缺失： ELPs(kt) = {Pn(kt)-[Ppv(kt)+Pw(kt)]nc}t (2) 其中，ELPs(kt)表示kt时刻的能量缺失量。4. 根据权利要求1所述的一种并网风光发电中混合储能容量优化配置方法，其特征在 于:所述步骤3)中，以混合储能装置的年均费用最小为目标函数，具体表达式如下： IllinC - Civ+C〇c+Cmc+Cdc -( l+k〇cba+kmcba+kdcba) * kdebafflfba+( l+k〇cuc+kdcuc)kdeucnfuc ( 14) 其中，c表示混合储能装置的年均费用，Clv指混合储能装置的购置成本，fba、fu。分别代 表蓄电池组和超级电容器组的单价，m、η分别表示蓄电池组和超级电容器组的个数，kdeba、 kdeuc分别表示蓄电池组和超级电容器组的年折旧值，(：。。指混合储能装置的运行成本，kQcba、 k。?。分别表示蓄电池和超级电容器的运行成本系数，Cm。指混合储能装置的维护成本，kmcba 表示蓄电池的维护成本系数，Cdc指混合储能装置的处置成本，kdcba、kdcuc分别表示蓄电池和 超级电容器的处置成本系数。5. 根据权利要求1所述的一种并网风光发电中混合储能容量优化配置方法，其特征在 于:所述步骤4)中，约束条件如下：其中，？^(1^，!11)、？11。(1^，11)分别表示蓄电池和超级电容器的充电功率，111、11分别表示蓄 电池组和超级电容器组的个数;？以1^)^(1^)几(1^)表示1^时刻光伏阵列输出功率、风 机输出功率和电网调度功率;η。表示逆变器的转换效率;P ban(m)、Eban(m)分别表示蓄电池的 额定功率和额定容量;α表示不平衡功率P ub= 中基本功率所 占比例;RlPPP (m，η )、RlPSP (m，η )分别表不能量损失率和能量缺失率，RLPPPmax、RLPSPmax分别表不 能量损失率最大值和能量缺失率最大值;PUMax(n)表示超级电容器最大充电功率;Ebamin(m) 表示蓄电池的最小容量;Eb a(kt，m)表示蓄电池经过充放电后的剩余能量;Euc；min(n)、E uc；max (η)分别表示超级电容器组能量的最小值和最大值;Euc(kt，n)表示超级电容器组经过充放 电后的剩余能量。6. 根据权利要求1所述的一种并网风光发电中混合储能容量优化配置方法，其特征在 于:所述步骤5)采用改进混沌优化算法求解优化配置模型，具体步骤如下： (5-1)算法初始化，令1^ = 0，给出~个具有微小差异的初值^()，且&,()辛1且&,()辛0.75且 Xi,o辛0.5且Xi,Q辛0.25且Xi,Q辛0,然后按照Logistic映射表达式得到N个轨迹不同的混纯变 量{xi，k}; (5-2)将生成的混沌变量{Xl,k}按下式进行第一次载波，得到优化变量范围[Xmin， Xmax] 内的混纯变量{χ7 i,k}: X i, k - Xmin+( Xmax_Xmin) Xi，k (17); (5-3)对生成的混沌变量进行可行性检验，重复执行步骤(5-1)和步骤(5-2)进行 迭代搜索得到满足约束条件的可行解Ix^.k}，并计算相应的目标函数值f (k) =f ({x"^}); 令)(*={^1,()}，产=以0)，如果以1〇<产，那么)(*={^1,1{}，产=以1〇;否则，令1^=1^1，进 入下一循环，直到产变化小于〇. 01或达到最大迭代次数Nmaxl; (5-4)按下式对混沌变量Xl,k+1进行第二次载波，得到混沌变量{V 1>k}: z7 i,k=x*+Pxi,k+i (18) 其中，β为调节因子； (5-5)对混沌变量进行可行性检验，当其满足约束条件时，若则f = z' = fXz' i,k);否则令k = k+l，进入下一循环，直到产变化小于0.01或达到最大迭代次 数Nmax2，最终的X#即为最优解，f $为最优值。7. 根据权利要求6所述的一种并网风光发电中混合储能容量优化配置方法，其特征在 于:所述步骤(5-1)中，Logistic映射表达式为： Xi，k+l = yXi,k(l_Xi,k)，Xi,kG[〇，l] (16) 其中，μ为控制参数;xi,k为混沌变量。
【文档编号】H02J3/28GK106058900SQ201610581661
【公开日】2016年10月26日
【申请日】2016年7月21日 公开号201610581661.7, CN 106058900 A, CN 106058900A, CN 201610581661, CN-A-106058900, CN106058900 A, CN106058900A, CN201610581661, CN201610581661.7
【发明人】杨志超, 陆文伟, 葛乐, 马寿虎, 陆文涛, 顾佳易, 王蒙
【申请人】南京工程学院