基于鲁棒萤火虫-粒子群混合优化的风-光-燃料电池微电网频率控制方法与流程

本发明涉及一种微电网频率控制方法。

背景技术：

电力系统通常向位于不同位置的不同类型的负载供电。由于在一些偏远和孤立的地方无法获得这种常规电力，由于安装额外的输电和配电线路的地理限制，以及能源需求的日益增长，导致将可再生能源引入常规电力系统，如太阳能光伏、风力、微型水力、等。以满足日益增长的需求，并最大限度地减少环境污染、传输损耗。然而，这些可再生能源具有高度不确定性，气发电取决于天气情况。该现象可能会导致电力系统频率发生较大波动，威胁系统的稳定运行。为此，本专利针对含有风机、光伏、燃料电池以及各类储能系统的微网，提出一种基于萤火虫和粒子群混合优化技术用于调节pid控制器参数，使得微网在不同运行条件下(如风速变化和负荷需求变化)频率偏差最小。

首先，之前大多数研究主要针对电力系统内负荷不确定性造成的频率波动提出一些控制方法，但当光伏等可再生能源普及后，其可能会造成能量管理的难题。有研究针对风力-燃料电池-柴油发电机和储能元件组成的微网提出了一种小信号分析方法，但他们没有应用任何控制器来获得更好的频率控制曲线。

尽管在基于启发式优化技术的文献中有很多关于频率控制的工作，但是它们主要在传统的火电、水电系统中进行研究，在这些系统中，唯一的干扰是负载变化的形式。然而，在当前的研究中，频率控制问题是在微电网中研究的，该微电网存在间歇性的可再生能源，例如风能和太阳能光伏。这里，除了负载需求的变化之外，频率控制器还受到风速和太阳强度变化的影响。此外，在以负载需求、风速和太阳强度变化为形式的多输入扰动下，基于风/光伏/风力发电/燃料电池/储能系统微电网的能量管理是一个至关重要的问题。

技术实现要素：

为了克服已有技术的不足，本发明提供了一种基于鲁棒萤火虫-粒子群混合优化的风-光-燃料电池微电网频率控制方法，使得微网在不同运行条件下(如风速变化和负荷需求变化)频率偏差最小。

本发明解决其技术问题所采用的技术方案是：

一种基于鲁棒萤火虫-粒子群混合优化的风-光-燃料电池微电网频率控制方法，所述方法包括以下步骤：

步骤1、分布式电源建模，过程如下：

1.1风/光/燃料电池/柴油发电机模型

电力由风力、光伏、燃料电池和柴油发电机产生，以满足负荷需求；其中风能和光伏所产生的电功率10％-15％被水电解槽用于产生氢气，然后由燃料电池根据负载要求发电；柴油发电机可视为备用电源，在风力、光伏等其他电源不可用的情况下，它可以自动向所连接的负载供电，风电和光伏的传递函数可以忽略其非线性，用简单的一阶线性传递函数表示：

式中，kwpg和kpv是增益常数；twpg和tpv分别为风和光伏的时间常数；和ppvg为系统内第k各风机输出功率；pw为风机的机械功率；为太阳辐照度；

水电解槽、燃料电池和柴油发电机的一阶传递函数为:

式中，kae，kfc，kdeg是增益常数；tae、tfc、tdeg分别为电解槽、燃料电池、柴油发电机的时间常数；pfc为燃料电池输出功率；pae为电解槽输出功率；pdeg为柴油发电机输出功率；δf为系统频率偏差；

1.2储能系统模型

在非高峰时段或在太阳能强度和风速较高时产生的额外电能可以储能在飞轮、电池或超级电容等储能设备中，储存的电力可以在高峰负荷期间或风能和光伏发电不够用期间重新利用，飞轮、电池或超级电容储能系统的线性传递函数如下所示：

式中，kfes，kbes，kuc是增益常数，tfes，tbes，tuc分别是飞轮储能/电池储能/超级电容的时间常数；pfes为飞轮储能输出功率；pbes为电储能输出功率；puc为超级电容输出功率；

步骤2、构建用于频率控制的微电网模型，过程为：

微网总发电量pmg表示为:

pmg＝pt+pdeg+pfcg+ppvg±pfes±(pbesorpuc)(7)

其中

pt＝pwpg+ppvg-pae(8)

式中，pt为风机和光伏的净功率；

供需之间的电力平衡通过控制各发电单元和储能来实现，表示为：

δpe＝pmg-pl(9)

式中，pmg是微电网总发电量；pl是总电力需求；

频率偏差δf计算如下：

式中，ksc为微网的特征常数，此时，系统的传递函数gsys表示为：

式中，m＝ksys和d＝ksystsys分别为系统的等效惯性常数和阻尼常数；

步骤3，ff-pso混合优化pid控制，过程为：

采用萤火虫算法进行全局搜索，采用粒子群算法进行局部搜索，首先，利用萤火虫算法识别搜索空间的有效区域，然后采用粒子群算法进行下一阶段的挖掘；为了获得所提出的ff-pso方法的优越性，对萤火虫和粒子群的特性进行了协调，以获得最佳的控制效果。

进一步，所述ff-pso方法包括以下步骤：

3.1、首先初始化萤火虫随机种群，其中包括对萤火虫数量、荧光强度、迭代次数iter的设定；

3.2、随机初始化第i各萤火虫在目标函数搜索范围内的位置，并计算萤火虫发光亮度；

3.3、寻找光强最大的萤火虫，并更新萤火虫的位置；

3.4、判断是否达到最大迭代数或者要求精度，如果达到则进行下一步，否则转向3.2；

3.5、选择10只发光亮度最大的萤火虫作为最为pso的输入；

3.6、计算出每个萤火虫的适应度，并对每个萤火虫，用它的适应度值和个体极值比较，如果适应度值大于个体极值，则个体极值用适应度值替换；

3.7、判断是否达到最大迭代数或者要求精度，如果达到则停止，否则更新萤火虫位置和速度并产生新的种群。

本发明中，风电、光伏发电以及负荷的变化都具有不确定性，该不确定性将导致微网中发电量和需求量之间不匹配，造成电力系统频率偏差，从而威胁系统的稳定运行。为此，本专利从两方面去解决该问题。首先，在系统结构上除了微网本身自带的风机、光伏以外增加了燃料电池和各类储能系统。其次，提出一种萤火虫和粒子群混合优化技术用于调节pid控制器参数，使得微网在不同运行条件下(如风速变化和负荷需求变化)频率偏差最小。

本发明的有益效果主要表现在：使得微网在不同运行条件下(如风速变化和负荷需求变化)频率偏差最小。

附图说明

图1为基于风/光/燃料电池的微网框图。

图2为微电网能量管理策略流程图。

图3为微网线性模型。

图4为ff-pso算法的流程图。

具体实施方式

下面结合附图对本发明作进一步描述。

参照图1～图3，一种基于鲁棒萤火虫-粒子群混合优化的风-光-燃料电池微电网频率控制方法，包括以下步骤：

步骤1、分布式电源建模，过程如下：

各种分布式电源与储能系统集成在一起，以提高负荷供电的质量和可靠性，然而，由于风速和光照强度的变化，风力和光伏等发电资源具有随机性和不确定性，该特性可能会使得系统中有功功率和无功功率的不匹配度增加从而导致系统不稳定；因此，在系统中引入燃料电池，通过与风能、光伏以及电池、飞轮、超级电容等储能系统相结合，形成可靠性更好的混合动力系统；线性化电源和储能系统模型，忽略其非线性，并以简化的一阶传递函数形式表示，以降低其复杂性；虽然简化模型可能会丢失系统的一些动态信息，但由于该研究为小信号分析，负荷/风速变化的百分比较小，因此上述线性化方法对系统行为和响应的影响较小；

1.1风/光/燃料电池/柴油发电机模型

电力由风力、光伏、燃料电池和柴油发电机产生，以满足负荷需求；其中风能和光伏所产生的电功率10％-15％被水电解槽(ae)用于产生氢气，然后由燃料电池根据负载要求发电；柴油发电机可视为备用电源，在风力、光伏等其他电源不可用的情况下，它可以自动向所连接的负载供电，风电和光伏的传递函数可以忽略其非线性，用简单的一阶线性传递函数表示：

式中，kwpg和kpv是增益常数；twpg和tpv分别为风和光伏的时间常数；pwpgk和ppvg为系统内第k各风机输出功率；pw为风机的机械功率；为太阳辐照度；

水电解槽、燃料电池和柴油发电机一阶传递函数为:

式中，kae，kfc，kdeg是增益常数；tae、tfc、tdeg分别为电解槽、燃料电池、柴油发电机的时间常数；pfc为燃料电池输出功率；pae为电解槽输出功率；pdeg为柴油发电机输出功率；δf为系统频率偏差；

1.2储能系统模型

在非高峰时段或在太阳能强度和风速较高时产生的额外电能可以储能在飞轮、电池或超级电容等储能设备中。储存的电力可以在高峰负荷期间或风能和光伏发电不够用期间重新利用，飞轮、电池或超级电容储能系统的线性传递函数如下所示：

式中，kfes，kbes，kuc是增益常数，tfes，tbes，tuc分别是飞轮储能/电池储能超级电容的时间常数；pfes为飞轮储能输出功率；pbes为电储能输出功率；puc为超级电容输出功率；

步骤2、构建用于频率控制的微电网模型，过程为：

将各发电装置和储能系统连接形成的微电网如图1～图3所示。其中，图1为基于风/光/燃料电池的微网框图，图2为微电网能量管理策略流程图，图3为微网线性模型，将风能、光伏、燃料电池和柴油发电机等电力资源与电解槽结合起来，电解槽将风能或光伏产生的部分电能转化为氢气，作为燃料电池的输入；光伏、燃料电池和储能系统通常通过dc/dc和dc/ac连接到负载；为了简化模型本文未考虑这些变流器模型，通过协调风机-光伏-燃料电池-柴油发电机各自的优点和缺点，大大减少了能量管理问题；此外，为了提高供电的可靠性，在微网内添加了各类储能系统，有助于在高峰负荷需求期间减少电力不足；柴油发电机作为备用发电机补偿功率不足，这证明了孤立微电网配置的自主特性。

微网总发电量pmg表示为:

pmg＝pt+pdeg+pfcg+ppvg±pfes±(pbesorpuc)(7)

其中

pt＝pwpg+ppvg-pae(8)

式中，pt为风机和光伏的净功率。

供需之间的电力平衡通过控制各发电单元和储能来实现，可表示为:

δpe＝pmg-pl(9)

式中，pmg是微电网总发电量；pl是总电力需求。

频率偏差δf计算如下：

式中，ksc为微网的特征常数。此时，系统的传递函数gsys表示为:

式中，m＝ksys和d＝ksystsys分别为系统的等效惯性常数和阻尼常数。在建模中考虑了这些参数，因为旋转机械的重质量转子会引入惯性和随频率变化的负载分量。

步骤3，ff-pso混合优化pid控制，过程为：

pid控制器是用途最广、应用最广泛的反馈控制器，因其优良的动态性能和鲁棒特性而受到人们的青睐，它依赖于三个基本控制器参数即比例增益(kp),积分增益(ki)和微分增益(kd)。这些增益通过萤火虫-粒子群(fireflyandparticleswarm，ff-pso)混合优化技术进行优化调整从而减小系统频率偏差以维持微电网的稳定；

粒子群优化是一种由eberhart和kennedy提出的基于鸟群的社会行为而设计的演变算法技术。粒子群优化算法是一种基于种群的优化工具，每个粒子根据自己和邻居的经验更新自己的位置和速度，该算法因其简单、收敛速度快、程序易读等特点而受到广泛欢迎；类似地，萤火虫算法是一种基于种群的算法，其中萤火虫的特征在于它们通过生化过程或生物发光产生闪光，闪光作为交配的主要求偶信号，基于以下三种行为：(a)所有的萤火虫都是雌雄同体的，不管其他萤火虫的性别如何，它们都会被吸引；(b)萤火虫的吸引力大小与其发光亮度成正比。它们的吸引力与它们的光强成正比。因此，对于任何两个闪烁的萤火虫，不太亮的萤火虫会向亮的那一个移动。亮度与距离成正比，亮度越高说明两只萤火虫之间的距离越小。如果两只闪烁的萤火虫亮度相同，那么它们就会随机移动；(c)萤火虫的亮度由目标函数确定。在系统不确定性和参数变化的情况下，粒子群算法或ff算法用于调节pid控制器参数有时会出现问题。

因此，为了进一步搜索最优多峰空间和更好的优化性能，设计了一种称为混合萤火虫和粒子群算法的混合方法来改进微电网的频率控制，该混合算法采用萤火虫算法进行全局搜索，采用粒子群算法进行局部搜索，首先，利用萤火虫算法识别搜索空间的有效区域，然后采用粒子群算法进行下一阶段的挖掘；为了获得所提出的ff-pso方法的优越性，对萤火虫和粒子群的特性进行了协调，以获得最佳的控制效果。ff-pso算法具体步骤如图4所示，包括以下步骤。

3.1、首先初始化萤火虫随机种群，其中包括对萤火虫数量、荧光强度、迭代次数iter的设定；

3.2、随机初始化第i各萤火虫在目标函数搜索范围内的位置，并计算萤火虫发光亮度；

3.3、寻找光强最大的萤火虫，并更新萤火虫的位置；

3.4、判断是否达到最大迭代数或者要求精度，如果达到则进行下一步，否则转向3.2；

3.5、选择10只发光亮度最大的萤火虫作为最为pso的输入；

3.6、计算出每个萤火虫的适应度，并对每个萤火虫，用它的适应度值和个体极值比较，如果适应度值大于个体极值，则个体极值用适应度值替换；

3.7、判断是否达到最大迭代数或者要求精度，如果达到则停止，否则更新萤火虫位置和速度并产生新的种群。

本说明书的实施例所述的内容仅仅是对发明构思的实现形式的列举，仅作说明用途。本发明的保护范围不应当被视为仅限于本实施例所陈述的具体形式，本发明的保护范围也及于本领域的普通技术人员根据本发明构思所能想到的等同技术手段。