CN104701871A - 一种含风光水多源互补微网混合储能容量最佳配比方法 - Google Patents

Abstract

一种含风光水多源互补微网混合储能容量最佳配比方法，该方法根据当地风、光、水等自然资源的分布情况，模拟风力发电、光伏发电、以及水力发电的年输出功率曲线，结合微网的年负荷曲线，以系统成本和功率波动为目标函数，以蓄电池容量和超级电容容量为优化变量，同时确定功率平衡约束、最大瞬时功率约束、供电可靠性约束、超级电容充放电电流及电压约束、蓄电池SOC约束等约束条件，建立含风、光、水的微网混合储能优化配置模型；采用含模糊决策的多目标规划GA-PSO算法对目标函数进行优化求解，得到混合储能容量的最佳配比。本发明与传统GA算法和PSO算法相比，收敛速度更快，且较好地回避了多目标优化算法中目标函数相互冲突的问题。

Description

一种含风光水多源互补微网混合储能容量最佳配比方法

技术领域

本发明涉及一种含风光水多源互补微网混合储能容量最佳配比方法，属微网混合储能技术领域。

背景技术

与煤炭和石油等传统化石能源相比，风能、太阳能以及水能等能源具有污染少、储量大、可再生等特点，加大这些清洁能源的开发力度，实现其大规模的利用，对于缓解当今世界严重的环境污染和资源枯竭问题具有重要的意义，已成为当前电力行业的一个工作重点。

然而，风能、太阳能和水能等可再生清洁能源均具有波动性的特点，其中风能和太阳能的短期波动较大，而水能的短期波动较小，季节性波动较大。这种波动的差异性，使得风力发电、光伏发电以及水力发电之间具有一定的互补性，但由风、光、水构成的多源互补微网输出功率仍不平稳，给大电网的安全稳定运行带来了新的挑战。配置储能装置是平抑含风光水多源互补微网输出功率波动的一种有效措施。蓄电池具有能量密度高、功率密度低，充放电速度慢的特点，而超级电容器具备功率密度高、充放电速度快的特点。构建含蓄电池和超级电容器的混合储能系统，充分发挥两种储能装置的各自优势，不仅可以降低储能系统的成本，而且可以实现含风光水多源互补微网与大电网的友好互动。

发明内容

本发明的目的是，针对可再生清洁能源风、光、水构成的多源互补微网输出功率波动的问题，提出一种含风光水多源互补微网混合储能容量最佳配比方法。

本发明的技术方案是，根据当地风、光、水等自然资源的分布情况，模拟风力发电、光伏发电、以及水力发电的年输出功率曲线；结合微网的年负荷曲线，建立以系统成本和功率波动最小的优化目标函数，同时确定功率平衡约束、最大瞬时功率约束、供电可靠性约束、超级电容充放电电流及电压约束、蓄电池SOC约束等约束条件，采用含模糊决策的多目标规划GA-PSO算法对目标函数进行优化求解，得到最优的混合储能容量配比。

所述方法的步骤如下：

(1)获取当地风、光、水等自然资源的分布情况，模拟风力发电、光伏发电、以及水力发电的年输出功率曲线；

(2)获取微网的年负荷曲线；

(3)以系统成本和功率波动最小为目标函数，以蓄电池容量和超级电容容量为优化变量，同时确定功率平衡约束、最大瞬时功率约束、供电可靠性约束、超级电容充放电电流及电压约束、蓄电池SOC约束等约束条件，建立含风、光、水的互补微网混合储能优化配置模型；

(4)采用含模糊决策的多目标规划GA-PSO算法对目标函数进行优化求解；

(5)统计不同目标函数非负权重系数下各种蓄电池容量与超级电容容量配比，找到最佳配比。

所述目标函数包括：

(1)成本目标函数

对混合储能系统的生命周期费用进行分析，建立了包含购置、运行、维护以及处理费用在内的成本目标函数：

minF₁＝C_b+C_o+C_m+C_d

其中，C_b为购置费用；C_o为运行费用；C_m为维护费用；C_d为处理费用；

(2)功率波动目标函数

取混合储能参与平滑功率波动的时间段为[T₁,T₂]，建立功率波动目标函数：

$\min F_2=\underset{i=T_1}{\overset{T_2}{\mathrm{\Σ}}}{(P_{\mathrm{DG},i}-P_{\mathrm{DG},i-1})}^2$

P_DG,i＝P_dg,i+P_bat,i+P_uc,i

$P_{\mathrm{dg},i}=\frac{E_i\left(w\right)}{T_2-T_1}+P_{s,i}+\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\left(P_{\mathrm{h\; j},i}\right)$

其中，i表示某一段时间；P_bat,i、P_uc,i分别表示该段时间内蓄电池和超级电容的输出功率；P_dg,i、P_DG,i分别表示经过混合储能装置平滑功率波动前后的系统实际发电功率；表示该段时间内风机的输出功率；表示水电集群在该段时间内的总输出功率，n为水电站个数；

利用目标加权法，定义个体的适应度函数为各目标函数的加权和：minF＝ω₁F₁+ω₂F₂，其中ω₁，ω₂是目标函数的非负权重系数，且满足ω₁+ω₂＝1。

所述含风、光、水的互补微网混合储能优化配置模型包括：

(1)风力发电单元

假定当前风速为v_i、风机启动风速为v_m、额定风速为v_H、停机风速为v_T、额定功率为P_N，则风机在不同风速时对应的发电量分别为：

$\left\{\begin{array}{cc}{E}_{w1}=\mathrm{\&Sigma;}[P_N\&CenterDot;{(v_i/v_H)}^3\&CenterDot;h_i]& (v_m\&le;v_i<v_H)\\ E_{w2}=\mathrm{\&Sigma;}(P_N\&CenterDot;h_i)& (v_H\&le;v_i<v_T)\end{array}\right.$

根据上式可以得到风机任一时间段内的发电总量：E_w＝E_w1+E_w2。

其中，h_i为保持当前风速v_i的小时数，单位为h；

(2)光伏发电单元

假设光伏阵列的额定功率为P_N，转换效率为η，则光伏阵列的发电量为：

E_s＝P_s·T_m＝N·P_N·η·T_m

其中，P_s为光伏阵列的实际输出功率，单位为kW；N为当月天数，单位为d；T_m为光伏阵列倾斜面上当月平均峰值日照时数，单位为h/d；

(3)水力发电单元

假设水电站的水轮机引水发电流量为Q，单位为m³/s；水电站的工作水头为H，m；水电站的水能转换效率为η，可表示为η＝η₁η₂η₃，其中，η₁表示水轮机效率，η₂表示发电机效率，η₃表示机组传动效率，则水电机组的发电量为：

E_h＝P_h·T＝A·Q·H·T

其中，P_h为水电站的实际输出功率，单位为kW；A的值可根据水电站的规模进行判定，一般A＝6.0～8.5；T为其机组发电小时数，单位为h；

(4)蓄电池单元

忽略蓄电池的极化内阻与相间微分电容，只考虑其端电压u_b和荷电状态SOC，可将蓄电池等效为由直流电压源E与内阻R_b串联构成的电路；

蓄电池输出电压u_b与荷电状态SOC可表示为：

u_b＝E-R_b·i_b

$\mathrm{SOC}={\mathrm{SOC}}_0+\frac{E_{\mathrm{bat}}}{Q}$

其中，E为蓄电池空载电压，V；R_b为其内阻，单位为Ω；i_b为其工作电流，单位为A；Q为其容量，单位为Ah；SOC₀为其初始荷电状态；为蓄电池的储能容量；

(5)超级电容单元

超级电容器可以等效为理想电容器C与电阻R_S相串联的模型；其中，U_C(t)为超级电容器的电压；I_S1、I_S2分别为超级电容器的充、放电电流。

设U₀为电容器初始电压，U₁(t)和U₂(t)分别是电容器充电和放电时的电压值，则电容器在充电和放电时，其计算公式分别为：

$U_1\left(t\right)=U_C\left(t\right)+R_S\&CenterDot;I_{S1}=U_0+\frac{1}{C}{I}_{S1}\left(t\right)+R_S\&CenterDot;I_{S1}$

$U_2\left(t\right)=U_C\left(t\right)-R_S\&CenterDot;I_{S2}=U_0+\frac{1}{C}{I}_{S2}\left(t\right)+R_S\&CenterDot;I_{S2}$

超级电容器充、放电时间分别表示为：

C·dv＝I_S1·t

C·dv-I_S2·C·R_S＝I_S2·t

其中，dv为超级电容器的端电压变化；t为超级电容的充放电时间，单位为h；

设U_w、U_v分别为电容充电完成与放电完成时的电压值，则超级电容的储能容量E_uc为：

$E_{\mathrm{uc}}=\frac{1}{2}\&CenterDot;C\&CenterDot;(U_w^2-U_v^2).$

所述约束条件包括：

(1)功率平衡约束：

在任何时刻，微网系统都必须满足功率平衡：

P_bat,i+P_uc,i+P_dg,i≥P_L,i

其中，P_L,i为该时间段微网负荷功率；

(2)供电可靠性约束：

微网系统一段时间内提供的总能量不小于负荷的用电量。

$E_{\mathrm{bat},i-1}+E_{\mathrm{uc},i-1}+E_{w,i}+E_{w,i}+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}\left(E_{\mathrm{h\; j},i}\right)-E_{L,i}\&GreaterEqual;0$

其中，E_bat,i-1、E_uc,i-1分别为蓄电池与超级电容前一个工作时间段所剩电量，E_w,i为该时间段内风力总发电量，E_s,i为该时间段内光伏阵列总发电量，为该时间段内水电集群总发电量，E_L,i为该时间段负荷用电量；

(3)蓄电池SOC约束：

在蓄电池的充放电过程中，必须严格控制其荷电状态的上下限，以延蓄电池的使用寿命：

SOC_min＜SOC_i＜SOC_max

其中，SOC_min为蓄电池荷电状态的下限；SOC_max为蓄电池荷电状态的上限；

(4)超级电容充、放电电流及电压约束：

超级电容在任何时候，必须满足其充、放电电流及电压约束。

$\left\{\begin{array}{c}0<I_S\&le;I_{S\max }\\ 0<U_1\left(t\right)\&le;U_{\max }\\ U_{\min }<U_2\left(t\right)\&le;U_{\max }\end{array}\right.$

其中：U_min、U_max分别为超级电容最低与最高充、放电电压；I_min、I_max分别为超级电容最低与最高充、放电电流；

(5)最大瞬时功率约束：

在任何时刻，微网系统都必须要满足其最大瞬时功率的缺失：

$\left\{\begin{array}{c}{P}_{\mathrm{bat},i}+P_{\mathrm{uc},i}\&GreaterEqual;\mathrm{\&Delta;}{P}_{L\max }\\ \underset{o}{\overset{t}{\&Integral;}}{P}_{\mathrm{bat},i}\mathrm{dt}+\underset{0}{\overset{t}{\&Integral;}}{P}_{\mathrm{uc},i}\&GreaterEqual;\mathrm{\&Delta;}{P}_{L\max }\&CenterDot;t\end{array}\right.$

其中，ΔP_Lmax为该时间段内微网的最大瞬时功率缺失。

所述对目标函数进行优化求解过程为：

Step1初始化：确定种群规模pop_size，交叉概率p_c(一般取0.40-0.99)，变异概率p_m(一般取0.0001-0.1)，最大进化代数maxgen，产生初始种群pop；

Step2更新：求出pop中所有个体的适应度值，得出个体最优粒子gbest和全局最优粒子zbest，并按粒子群算法进行速度、位置更新，得到新的种群pop1；

Step 3交叉：按照交叉概率p_c，执行交叉算子；

Step4变异：按照变异概率p_m，执行交叉算子，得到新的种群pop2；

Step5选择：从种群pop1与pop 2中用精英选择算子选出pop_size个个体组成下一代种群pop，同时更新全局最优粒子zbest；

Step6判断：若满足迭代次数或满足最优个体所对应的网络输出误差精度要求，则停止；否则，转step2；

Step7模糊决策：改变目标函数的非负权重系数，多次重复计算，用以进行模糊决策。

本发明的有益效果是，本发明的一种含风光水多源互补微网混合储能容量最佳配比方法，具有如下优点：该方法考虑了风、光、水等多种清洁能源，具有较强的通用性；该方法与传统的GA算法和PSO算法相比，收敛速度更快；该方法采用含模糊决策的多目标规划GA-PSO算法较好地回避了多目标优化算法中目标函数相互冲突的问题。

本发明适用于含风光水多源互补微网混合储能容量最佳配比。

附图说明

图1为是本发明的一种含风光水多源互补微网混合储能容量最佳配比方法流程图；

图2是本发明的含模糊决策的多目标规划GA-PSO算法流程图；

图3是本发明的关于模糊决策的三角形隶属度函数。

具体实施方式

以下将结合图1和图2对本发明作进一步说明。

本发明的一种含风光水多源互补微网混合储能容量最佳配比方法具体实施方式如下：

(1)获取当地风、光、水等自然资源的分布情况，模拟风力发电、光伏发电、以及水力发电的年输出功率曲线；

(2)获取微网的年负荷曲线；

(3)以系统成本和功率波动最小为目标函数，以蓄电池容量和超级电容容量为优化变量，同时确定功率平衡约束、最大瞬时功率约束、供电可靠性约束、超级电容充放电电流及电压约束、蓄电池SOC约束等约束条件，建立含风、光、水的互补微网混合储能优化配置模型；

(4)采用含模糊决策的多目标规划GA-PSO算法对目标函数进行优化求解，其过程如下：

Step1初始化：确定种群规模pop_size，交叉概率p_c(一般取0.40-0.99)，变异概率p_m(一般取0.0001-0.1)，最大进化代数maxgen，产生初始种群pop。

Step2更新：求出pop中所有个体的适应度值，得出个体最优粒子gbest和全局最优粒子zbest，并按粒子群算法进行速度、位置更新，得到新的种群pop1；

Step 3交叉：按照交叉概率p_c，执行交叉算子；

Step4变异：按照变异概率p_m，执行交叉算子，得到新的种群pop2；

Step5选择：从种群pop1与pop 2中用精英选择算子选出pop_size个个体组成下一代种群pop，同时更新全局最优粒子zbest；

Step6判断：若满足迭代次数或满足最优个体所对应的网络输出误差精度要求，则停止；否则，转step2。

Step7模糊决策：改变目标函数的非负权重系数，多次重复计算，用以进行模糊决策。

(5)统计不同目标函数非负权重系数下各种蓄电池容量与超级电容容量配比，找到最佳配比。

下面针对本专利中涉及的含风光水互补微网混合储能系统模型、目标函数和约束条件作如下说明。

本实施例含风光水互补微网混合储能系统模型：

(1)风力发电单元

假定当前风速，风机启动风速、额定风速、停机风速与额定功率分别为v_i、v_m、v_H、v_T、P_N，则风机在不同风速时对应的发电量分别为：

$\left\{\begin{array}{cc}{E}_{w1}=\mathrm{\&Sigma;}[P_N\&CenterDot;{(v_i/v_H)}^3\&CenterDot;h_i]& (v_m\&le;v_i<v_H)\\ E_{w2}=\mathrm{\&Sigma;}(P_N\&CenterDot;h_i)& (v_H\&le;v_i<v_T)\end{array}\right.$

根据上式可以得到风机任一时间段内的发电总量：E_w＝E_w1+E_w2。

其中，h_i为保持当前风速v_i的小时数，单位为h。

(2)光伏发电单元

假设光伏阵列的额定功率为P_N，转换效率为η，则光伏阵列的发电量为：

E_s＝P_s·T_m＝N·P_N·η·T_m

其中，P_s为光伏阵列的实际输出功率，单位为kW；N为当月天数，单位为d；T_m为光伏阵列倾斜面上当月平均峰值日照时数，单位为h/d。

(3)水力发电单元

假设水电站的水轮机引水发电流量为Q，单位为m³/s；水电站的工作水头为H，单位为m；水电站的水能转换效率为η，可表示为η＝η₁η₂η₃，η₁、η₂、η₃分别表示水轮机效率、发电机效率、机组传动效率，则水电机组的发电量为：

E_h＝P_h·T＝A·Q·H·T

其中，P_h为水电站的实际输出功率，单位为kW；A的值可根据水电站的规模进行判定，一般A＝6.0～8.5；T为其机组发电小时数，单位为h。

(4)蓄电池单元

忽略蓄电池的极化内阻与相间微分电容，只考虑其端电压u_b和荷电状态SOC，可将蓄电池等效为由直流电压源E与内阻R_b串联构成的电路。

蓄电池输出电压u_b与荷电状态SOC可表示为：

u_b＝E-R_b·i_b

$\mathrm{SOC}={\mathrm{SOC}}_0+\frac{E_{\mathrm{bat}}}{Q}$

其中，E为蓄电池空载电压，单位为V；R_b为其内阻，单位为Ω；i_b为其工作电流，单位为A；Q为其容量，单位为Ah；SOC₀为其初始荷电状态；为蓄电池的储能容量。

(5)超级电容单元

超级电容器可以等效为理想电容器C与电阻R_S相串联的模型。其中，U_C(t)为超级电容器的电压；I_S1、I_S2分别为超级电容器的充、放电电流。

设U₀为电容器初始电压，U₁(t)和U₂(t)分别是电容器充电和放电时的电压值，则电容器在充电和放电时，其计算公式分别为：

$U_1\left(t\right)=U_C\left(t\right)+R_S\&CenterDot;I_{S1}=U_0+\frac{1}{C}{I}_{S1}\left(t\right)+R_S\&CenterDot;I_{S1}$

$U_2\left(t\right)=U_C\left(t\right)-R_S\&CenterDot;I_{S2}=U_0+\frac{1}{C}{I}_{S2}\left(t\right)+R_S\&CenterDot;I_{S2}$

超级电容器充、放电时间分别表示为：

C·dv＝I_S1·t

C·dv-I_S2·C·R_S＝I_S2·t

其中，dv为超级电容器的端电压变化；t为超级电容的充放电时间，单位为h。

设U_w、U_v分别为电容充电完成与放电完成时的电压值，则超级电容的储能容量E_uc为：

$E_{\mathrm{uc}}=\frac{1}{2}\&CenterDot;C\&CenterDot;(U_w^2-U_v^2).$

本实施例的目标函数：

(1)成本目标函数

对混合储能系统的生命周期费用进行分析，建立了包含购置、运行、维护以及处理费用在内的成本目标函数：

minF₁＝C_b+C_o+C_m+C_d

其中，C_b为购置费用；C_o为运行费用；C_m为维护费用；C_d为处理费用。

(2)功率波动目标函数

取混合储能参与平滑功率波动的时间段为[T₁,T₂]，建立功率波动目标函数：

$\min F_2=\underset{i=T_1}{\overset{T_2}{\mathrm{\&Sigma;}}}{(P_{\mathrm{DG},i}-P_{\mathrm{DG},i-1})}^2$

P_DG,i＝P_dg,i+P_bat,i+P_uc,i

$P_{\mathrm{dg},i}=\frac{E_i\left(w\right)}{T_2-T_1}+P_{s,i}+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}\left(P_{\mathrm{h\; j},i}\right)$

其中，i表示某一段时间；P_bat,i、P_uc,i分别表示该段时间内蓄电池和超级电容的输出功率；P_dg,i、P_DG,i分别表示经过混合储能装置平滑功率波动前后的系统实际发电功率；表示该段时间内风机的输出功率；表示水电集群在该段时间内的总输出功率，n为水电站个数。

利用目标加权法，定义个体的适应度函数为各目标函数的加权和：minF＝ω₁F₁+ω₂F₂，其中ω₁，ω₂是目标函数的非负权重系数，且满足ω₁+ω₂＝1。

本实施例的约束条件：

(1)功率平衡约束：

在任何时刻，微网系统都必须满足功率平衡。

P_bat,i+P_uc,i+P_dg,i≥P_L,i

其中，P_L,i为该时间段微网负荷功率。

(2)供电可靠性约束：

微网系统一段时间内提供的总能量不小于负荷的用电量。

$E_{\mathrm{bat},i-1}+E_{\mathrm{uc},i-1}+E_{w,i}+E_{w,i}+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}\left(E_{\mathrm{h\; j},i}\right)-E_{L,i}\&GreaterEqual;0$

其中，E_bat,i-1、E_uc,i-1分别为蓄电池与超级电容前一个工作时间段所剩电量，E_w,i为该时间段内风力总发电量，E_s,i为该时间段内光伏阵列总发电量，为该时间段内水电集群总发电量，E_L,i为该时间段负荷用电量。

(3)蓄电池SOC约束：

在蓄电池的充放电过程中，必须严格控制其荷电状态的上下限，以延蓄电池的使用寿命。

SOC_min＜SOC_i＜SOC_max

(4)超级电容充、放电电流及电压约束：

超级电容在任何时候，必须满足其充、放电电流及电压约束。

$\left\{\begin{array}{c}0<I_S\&le;I_{S\max }\\ 0<U_1\left(t\right)\&le;U_{\max }\\ U_{\min }<U_2\left(t\right)\&le;U_{\max }\end{array}\right.$

其中：U_min、U_max分别为超级电容最低与最高充、放电电压；I_min、I_max分别为超级电容最低与最高充、放电电流。

(5)最大瞬时功率约束：

在任何时刻，微网系统都必须要满足其最大瞬时功率的缺失。

$\left\{\begin{array}{c}{P}_{\mathrm{bat},i}+P_{\mathrm{uc},i}\&GreaterEqual;\mathrm{\&Delta;}{P}_{L\max }\\ \underset{o}{\overset{t}{\&Integral;}}{P}_{\mathrm{bat},i}\mathrm{dt}+\underset{0}{\overset{t}{\&Integral;}}{P}_{\mathrm{uc},i}\&GreaterEqual;\mathrm{\&Delta;}{P}_{L\max }\&CenterDot;t\end{array}\right.$

其中，ΔP_Lmax为该时间段内微网的最大瞬时功率缺失。

本实施例的模糊决策如下：

为了更直观、合理的描述目标函数的重要程度，定义目标函数的重要程度为语言变量，并定义模糊子集为：[很不重要，不重要，中等，重要，很重要]。三角形隶属度函数对各目标函数的重要程度做出相应的模糊描述，如图3所示。根据三角形隶属度函数，得出各模糊子集的隶属度函数。以目标函数的重要度作为论域，把个体对应的ω代入各段隶属度函数进行计算，得到ω对各模糊子集的隶属度，然后用隶属度最大的模糊子集对此个体进行描述。在取得Pareto最优解集后，通过模糊评价法对Pareto最优解集进行分类，得到9类结果。

Claims (7)

1.一种含风光水多源互补微网混合储能容量最佳配比方法，其特征在于，所述方法根据当地风、光、水自然资源的分布情况，模拟风力发电、光伏发电、以及水力发电的年输出功率曲线；结合微网的年负荷曲线，建立以系统成本和功率波动最小的优化目标函数，同时确定功率平衡约束条件、最大瞬时功率约束条件、供电可靠性约束条件、超级电容充放电电流及电压约束条件和蓄电池SOC约束条件，采用含模糊决策的多目标规划GA-PSO算法对目标函数进行优化求解，得到最优的混合储能容量配比。

2.根据权利要求1所述的一种含风光水多源互补微网混合储能容量最佳配比方法，其特征在于，所述方法的步骤如下：

(1)获取当地风、光、水自然资源的分布情况，模拟风力发电、光伏发电、以及水力发电的年输出功率曲线；

(2)获取微网的年负荷曲线；

(3)以系统成本和功率波动最小为目标函数，以蓄电池容量和超级电容容量为优化变量，同时确定包括功率平衡约束、最大瞬时功率约束、供电可靠性约束、超级电容充放电电流及电压约束、蓄电池SOC约束在内的约束条件；建立含风、光、水的微网混合储能优化配置模型；

(4)采用含模糊决策的多目标规划GA-PSO算法对目标函数进行优化求解；

(5)统计不同目标函数非负权重系数下各种蓄电池容量与超级电容容量配比，找到最佳配比。

3.根据权利要求2所述的一种含风光水多源互补微网混合储能容量最佳配比方法，其特征在于，所述对目标函数进行优化求解过程为：

Step1初始化：确定种群规模pop_size，交叉概率p_c，变异概率p_m，最大进化代数maxgen，产生初始种群pop；

Step2更新：求出pop中所有个体的适应度值，得出个体最优粒子gbest和全局最优粒子zbest，并按粒子群算法进行速度、位置更新，得到新的种群pop1；

Step3交叉：按照交叉概率p_c，执行交叉算子；

Step4变异：按照变异概率p_m，执行交叉算子，得到新的种群pop2；

Step5选择：从种群pop1与pop2中用精英选择算子选出pop_size个个体组成下一代种群pop，同时更新全局最优粒子zbest；

Step6判断：若满足迭代次数或满足最优个体所对应的网络输出误差精度要求，则停止；否则，转step2；

Step7模糊决策：改变目标函数的非负权重系数，多次重复计算，用以进行模糊决策。

4.根据权利要求2所述的一种含风光水多源互补微网混合储能容量最佳配比方法，其特征在于，所述目标函数包括：

(1)成本目标函数

对混合储能系统的生命周期费用进行分析，建立了包含购置、运行、维护以及处理费用在内的成本目标函数：

minF₁＝C_b+C_o+C_m+C_d

其中，C_b为购置费用；C_o为运行费用；C_m为维护费用；C_d为处理费用；

(2)功率波动目标函数

取混合储能参与平滑功率波动的时间段为[T₁,T₂]，建立功率波动目标函数：

$\min F_2=\underset{i=T_1}{\overset{T_2}{\mathrm{\&Sigma;}}}{(P_{\mathrm{DG},i}-P_{\mathrm{DG},i-1})}^2$

P_DG,i＝P_dg,i+P_bat,i+P_uc,i

$P_{\mathrm{dg},i}=\frac{E_i\left(w\right)}{T_2-T_1}+P_{s,i}+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}\left(P_{\mathrm{hj},i}\right)$

其中，i表示某一段时间；P_bat,i、P_uc,i分别表示该段时间内蓄电池和超级电容的输出功率；P_dg,i、P_DG,i分别表示经过混合储能装置平滑功率波动前后的系统实际发电功率；表示该段时间内风机的输出功率；表示水电集群在该段时间内的总输出功率，n为水电站个数；

利用目标加权法，定义个体的适应度函数为各目标函数的加权和：minF＝ω₁F₁+ω₂F₂，其中ω₁，ω₂是目标函数的非负权重系数，且满足ω₁+ω₂＝1。

5.根据权利要求2所述的一种含风光水多源互补微网混合储能容量最佳配比方法，其特征在于，所述含风、光、水的微网混合储能优化配置模型为：

(1)风力发电单元

假定当前风速为v_i、风机启动风速为v_m、额定风速为v_H、停机风速为v_T、额定功率为P_N，则风机在不同风速时对应的发电量分别为：

$\left\{\begin{array}{cc}{E}_{w1}=\mathrm{\&Sigma;}[P_N\&CenterDot;{(v_i/v_H)}^3\&CenterDot;h_i]& (v_m\&le;v_i<v_H)\\ E_{w2}=\mathrm{\&Sigma;}(P_N\&CenterDot;h_i)& (v_H\&le;v_i<v_T)\end{array}\right.$

根据上式可以得到风机任一时间段内的发电总量：E_w＝E_w1+E_w2；

其中，h_i为保持当前风速v_i的小时数，单位为h；

(2)光伏发电单元

假设光伏阵列的额定功率为P_N，转换效率为η，则光伏阵列的发电量为：

E_s＝P_s·T_m＝N·P_N·η·T_m

其中，P_s为光伏阵列的实际输出功率，单位为kW；N为当月天数，单位为d；T_m为光伏阵列倾斜面上当月平均峰值日照时数，单位为h/d；

(3)水力发电单元

假设水电站的水轮机引水发电流量为Q，单位为m³/s；水电站的工作水头为H，m；水电站的水能转换效率为η，可表示为η＝η₁η₂η₃，其中，η₁表示水轮机效率，η₂表示发电机效率，η₃表示机组传动效率，则水电机组的发电量为：

E_h＝P_h·T＝A·Q·H·T

其中，P_h为水电站的实际输出功率，单位为kW；A的值可根据水电站的规模进行判定，一般A＝6.0～8.5；T为其机组发电小时数，单位为h；

(4)蓄电池单元

忽略蓄电池的极化内阻与相间微分电容，只考虑其端电压u_b和荷电状态SOC，可将蓄电池等效为由直流电压源E与内阻R_b串联构成的电路；

蓄电池输出电压u_b与荷电状态SOC可表示为：

u_b＝E-R_b·i_b

$\mathrm{SOC}={\mathrm{SOC}}_0+\frac{E_{\mathrm{bat}}}{Q}$

其中，E为蓄电池空载电压，V；R_b为其内阻，单位为Ω；i_b为其工作电流，单位为A；Q为其容量，单位为Ah；SOC₀为其初始荷电状态；为蓄电池的储能容量；

(5)超级电容单元

超级电容器可以等效为理想电容器C与电阻R_S相串联的模型；其中，U_C(t)为超级电容器的电压；I_S1、I_S2分别为超级电容器的充、放电电流；

设U₀为电容器初始电压，U₁(t)和U₂(t)分别是电容器充电和放电时的电压值，则电容器在充电和放电时，其计算公式分别为：

$U_1\left(t\right)=U_C\left(t\right)+R_S\&CenterDot;I_{S1}=U_0+\frac{1}{C}{I}_{S1}\left(t\right)+R_S\&CenterDot;I_{S1}$

$U_2\left(t\right)=U_C\left(t\right)-R_S\&CenterDot;I_{S2}=U_0+\frac{1}{C}{I}_{S2}\left(t\right)+R_S\&CenterDot;I_{S2}$

超级电容器充、放电时间分别表示为：

C·dv＝I_S1·t

C·dv-I_S2·C·R_S＝I_S2·t

其中，dv为超级电容器的端电压变化；t为超级电容的充放电时间，单位为h；

设U_w、U_v分别为电容充电完成与放电完成时的电压值，则超级电容的储能容量E_uc为：

$E_{\mathrm{uc}}=\frac{1}{2}\&CenterDot;C\&CenterDot;(U_w^2-U_v^2).$

6.根据权利要求2所述的一种含风光水多源互补微网混合储能容量最佳配比方法，其特征在于，所述约束条件包括：

(1)功率平衡约束：

在任何时刻，微网系统都必须满足功率平衡：

P_bat,i+P_uc,i+P_dg,i≥P_L,i

其中，P_L,i为该时间段微网负荷功率；

(2)供电可靠性约束：

微网系统一段时间内提供的总能量不小于负荷的用电量；

$E_{\mathrm{bat},i-1}+E_{\mathrm{uc},i-1}+E_{w,i}+E_{w,i}+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}\left(E_{\mathrm{hj},i}\right)-E_{L,i}\&GreaterEqual;0$

其中，E_bat,i-1、E_uc,i-1分别为蓄电池与超级电容前一个工作时间段所剩电量，E_w,i为该时间段内风力总发电量，E_s,i为该时间段内光伏阵列总发电量，为该时间段内水电集群总发电量，E_L,i为该时间段负荷用电量；

(3)蓄电池SOC约束：

在蓄电池的充放电过程中，必须严格控制其荷电状态的上下限，以延蓄电池的使用寿命：

SOC_min＜SOC_i＜SOC_max

其中，SOC_min为蓄电池荷电状态的下限；SOC_max为蓄电池荷电状态的上限；

(4)超级电容充、放电电流及电压约束：

超级电容在任何时候，必须满足其充、放电电流及电压约束；

$\left\{\begin{array}{c}0<I_S\&le;I_{S\max }\\ 0<U_1\left(t\right)\&le;U_{\max }\\ U_{\min }<U_2\left(t\right)\&le;N_{\max }\end{array}\right.$

其中：U_min、U_max分别为超级电容最低与最高充、放电电压；I_min、I_max分别为超级电容最低与最高充、放电电流；

(5)最大瞬时功率约束：

在任何时刻，微网系统都必须要满足其最大瞬时功率的缺失：

$\left\{\begin{array}{c}{P}_{\mathrm{bat},i}+P_{\mathrm{uc},i}\&GreaterEqual;{\mathrm{\&Delta;P}}_{L\max }\\ \underset{o}{\overset{t}{\&Integral;}}{P}_{\mathrm{bat},i}\mathrm{dt}+\underset{0}{\overset{t}{\&Integral;}}{P}_{\mathrm{uc},i}\&GreaterEqual;{\mathrm{\&Delta;P}}_{L\max }\&CenterDot;t\end{array}\right.$

其中，ΔP_Lmax为该时间段内微网的最大瞬时功率缺失。

7.根据权利要求3所述的一种含风光水多源互补微网混合储能容量最佳配比方法，其特征在于，所述交叉概率p_c取值为0.40-0.99；变异概率p_m取值为0.0001-0.1。