一种永磁风机PI参数优化方法及装置与流程

本发明属于风力发电领域，具体涉及一种适用于不同风速下基于自适应多维萤火虫算法的永磁风机的pi参数优化的方法。

背景技术：

能源是人类社会进步和经济发展的重要物质基础，能源工业的发展规模和速度影响着社会经济的发展，70年代初期，由于“石油危机”，出现了资源紧张的问题，人们意识到常规化石能源的不稳定性和有限性及其大量消耗并造成的环境污染和资源短缺等问题，使人们开始寻找清洁的可再生资源。相对于传统的石油、天然气等化石能源，新能源具有污染少的特点，对于解决当今世界资源短缺和改善环境污染，推动技术进步，实现人与自然可持续发展。风能是清洁的可再生资源，与传统化石能源相比，具有无环境污染、能源充足等优点。风力发电是将风能转换为机械能，进而转换为电能，其基本的工作原理是利用风吹动风轮，通过风轮的机械转动驱动发电机转子旋转，从而产生电能。在当今世界的新能源开发技术中，风力发电是最有利用价值的发电方式，全球风力发电量所占比也在逐步增加，并且发电技术也较为成熟。风力发电机一般包括叶片、发电机、尾翼、塔架、储能装置等构件组成。

永磁直驱风机是国内常用的风力发电机，其优点是实用低风速的能力强，耗能较低，并且维护成本也比较低。根据我国低风速的三类风区占全部风能资源的一半以上，更适合实用永磁直驱风机，也是我国未来的发展趋势。考虑到实际工作中的风电场，研究风力发电机最大风能捕获问题具有实际意义，也可以降低成本。风力发电机组机侧变流器是功率控制的关键环节。目前应用十分广泛的算法是定子磁链定向的矢量控制，可以实现高动态性能变频调速系统的设计，将电机定子电流分解成转矩电流分量和励磁电流分量去实现解耦，进而按照线性系统的控制理论方法来设计转速和磁链的pi控制器，从而可以实现连续控制，十分适用于风力发电机组此类宽范围调速系统和伺服系统。矢量控制一般采用pi控制器作为实现策略，具有结构简单，稳定性能好，可靠性高的优点。pi控制算法效果取决于控制器的环节参数kp、ki的取值组合，而风电场风功率是非线性的、还需要考虑一些因天气变化、频率变化、电网电压等不确定因素，而且风速通常以高度随机且快速的时变模式变化，发电机转速变化的离散型非常大，所以风力发电机很难实现最大风功率。传统的参数优化方法通常是操作者在实际系统中经过多次反复调试得出目标参数，这种方法效率低下，严重依赖于生产经验；现代控制方法如适应控制以实现全局一致的控制性能，实际难以实现。

技术实现要素：

针对现有技术的不足，本发明提出一种永磁风机pi参数优化方法及装置。具体而言，本发明提供了以下的技术方案：

一方面，本发明提供了一种永磁风机pi参数优化方法，所述方法包括：

s1、构造pi参数控制的目标函数，并检测当下风速；

s2、初始化萤火虫算法参数，所述算法参数包括种群维数、萤火虫数目、每组萤火虫种群数、最大吸引度、光强吸收系数、最大迭代次数、随机萤火虫位置；

s3、计算各组萤火虫吸引度、相对亮度、自适应移动步长，并在随后开始移动；

s4、更新萤火虫位置，并重新计算吸引度、相对亮度；

s5、计算搜索精度，当搜索精度达到预设要求或者迭代次数达到最大阈值时，执行s6，否则检测当下风速，并执行s2；

s6、若每组萤火虫满足搜索精度要求或迭代次数均达到最大阈值，则得到所述目标函数的全局最优解；

s7、输出上一阶段风速对应的最优值，若当下风速与所述上一阶段风速差值大于预设值，则获取当下风速，并返回s2，若小于等于预设值，则结束。

优选地，所述目标函数为：

其中，

ωm机械转动速度，ωm^*追踪最佳参考值，α1和α2为权重，pmax为风力发电机组的最大风能捕捉，vd、vq分别为定子电压的d-q分量，为瞬时定子电流d轴分量基准值，isd为瞬时定子电流d轴分量。

优选地，所述s2中，所述种群维数为7。

优选地，所述s3中，各组萤火虫i的吸引度为：

其中，β0为最大吸引度；γ为光强吸收系数；为当前时刻萤火虫i与j之间的距离

优选地，所述s3中，各组萤火虫i的相对亮度为：

其中，为萤火虫i当前时刻的最大荧光亮度，γ为光强吸收系数；为当前时刻萤火虫i与j之间的距离

优选地，所述s3中，萤火虫i自适应移动步长如下：

式中：表示第i只萤火虫的步长,riz表示第i只萤火虫中亮度与最高亮度萤火虫z之间的距离。

优选地，所述s4中，更新萤火虫位置采用如下方式：

其中，代表当前时刻萤火虫i所处位置，代表前一时刻萤火虫i所处位置；代表前一时刻萤火虫j所处位置；rand为[0,1]上服从均匀分布的随机因子。

优选地，所述搜索精度的预设要求为：搜索精度小于等于0.2。

优选地，所述s7中的预设值为3m/s。

另一方面，本发明还提供了一种永磁风机pi参数优化装置，所述装置包括：

风速检测模块，用于检测当下风速；

参数优化计算模块，基于目标函数以及当下风速，通过萤火虫算法计算所述目标函数的全局最优解，输出上一阶段风速对应的最优值，从而得到优化的pi参数；

风速判断模块，获取所述上一阶段风速，并基于所述当下风速以及上一阶段风速的差值，控制所述参数优化计算模块进行计算或停止。

优选地，所述风速判断模块中，若当下风速与所述上一阶段风速差值大于预设值，则获取当下风速，控制所述参数优化计算模块进行计算，若小于等于预设值，则结束所述参数优化计算模块的计算。

优选地，所述目标函数为：

其中，

ωm机械转动速度，ωm^*追踪最佳参考值，α1和α2为权重，pmax为风力发电机组的最大风能捕捉，vd、vq分别为定子电压的d-q分量，为瞬时定子电流d轴分量基准值，isd为瞬时定子电流d轴分量。

基于上述技术方案可见，与现有技术相比，本发明提高了pi参数优化过程中的寻优速度，有效的避免了陷入局部最优，并应用于变风速下这一复杂且在实际中常见的情况，解决了在变风速下传统手动测量数据的效率低下问题，结合本发明的技术方案即可实现高效率并保证了较高的精度。

附图说明

图1为本发明实施例的方法流程图；

图2为本发明实施例的pi控制装置结构示意图。

具体实施例

下面将结合本发明实施例中的图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1

在一个具体的实施例中，本发明通过改进的萤火虫算法实现永磁风机pi控制的优化。萤火虫算法(firelyalgorithm,fa)是模拟自然界中成虫发光的生物学特性发展而来的，是一种启发式算法也是基于群体的随机优化算法在萤火虫的群聚活动中，每只萤火虫通过散发荧光素与同伴进行寻觅食物以及求偶等信息交流。一般来说，荧光素越亮的萤火虫其号召力也就越强，最终会出现很多萤火虫聚集在一些荧光素较亮的萤火虫周围。萤火虫不分性别，它将会被吸引到所有其他比它更亮的萤火虫那去；萤火虫的吸引力和相对亮度成正比，对于任何两只萤火虫，其中一只会向着比它更亮的另一只移动，然而，相对亮度是随着距离的增加而减少的；相对比于其他算法，自适应步长的萤火虫算法可以提高收敛速度和精度，并且避免进入局部最优，可以达到理想的后期效果，实现最大风功率和最低控制成本的最优解。本发明基于风电机组应用的特殊要求，对萤火虫算法进行了改进，以更快地达到求解目标函数的最优解的目的。

为了使风机在变风速下达到最大风功率并且控制成本最低，本实施例中，选出需要优化的参数，结合算法具体包括以下步骤：

步骤一，构造出目的函数fp，使其实现最大风功率和最低控制成本；

步骤二，测当下风速，开始算法；

步骤三，初始化：萤火虫种群维数7维，萤火虫数目χ，每一组萤火虫种群数w，最大吸引度β0，光强吸收系数γ，最大迭代次数m，随机初始化萤火虫位置；

步骤四，计算各组萤火虫吸引度βi和相对亮度ii；

步骤五，计算萤火虫自适应移动步长αi，从而确定萤火虫移动步长，萤火虫开始移动；

步骤六，更新萤火虫位置，并重新计算吸引度，相对亮度；

步骤七，计算搜索精度，当精度达到要求或达到最大迭代次数时进行下一步，否则再次进行搜索，跳到步骤二；

步骤八，若7组萤火虫满足精度要求或最大迭代次数时，得到全局最优解，即找到目标函数fp最小，即成本最低和最大风功率；

步骤九，输出上一阶段风速最优值后，若此时风速与差值大于3m/s，则跳到步骤二，若小于，则结束算法，直至时，重启算法。

具体而言，本实施例的详细实现方法，可以通过如下方式进行：

步骤1：构造目标函数fp实现最大功率和最低成本：

永磁风机典型的解耦d-q电流pi控制方程可基于风电的最优求解设置如下：

vsd^*＝vsd'-ωelqisq

vsq^*＝vsq'+ωeldisd+ωeψf

其中rs是定子绕组的电阻，vsd、vsq、isd、isq分别是瞬时定子电压、电流的d-q分量，ωe是电气转速，ψf是由永磁体产生的磁链，ld、lq是定子的漏电感与定子和转子之间dq轴互感之和。d-q电压基准vsd^*和vsq^*是d-q电压的总和，vsd′和vsq′表示补偿项目。

其中，cp风能利用系数，β：风力发电机叶片的桨距角；λ：叶尖速比，ω：风轮旋转角速度，r：风轮半径；

风力发电机组的最大风能捕捉与发电机转速控制密切相关

kω＝0.5ρscρmax

ωt风力机风轮转速，n：齿轮箱变比；kω机械转速下的风力在风力涡轮机上产生的机械动力。

功率系数cp(λ，β)应在工作范围内的各种风速下均能够保持在最大点cp^*处，

cp^*＝cp(λ^*)

其中，λ^*：叶尖速比最佳值，ωm机械转动速度，ωm^*追踪最佳参考值：，桨距角β＝0°，vf风速。

控制成本通过正弦脉宽调制(spwm)技术调制，用g(x)表示

当fp取得最小值时，即为实现最大风功率和最低控制成本。

其中t为风速下的时间，kpi、kii(i＝3)为比例积分系数，由于有3个pi控制器，所以总共有6个参数，其余各项均在其最大值最小值之间变化，由于功率与成本的数量级相差甚远，所以引入权重α1和α2用来分权等，使控制成本不被过度优化。

步骤2：捕获风速，开始算法

步骤3：初始化

在本实施例中，优选地，设置萤火虫种群维数为7，7维萤火虫位置分别代表kp1、ki1、kp2、ki2、kp3、ki3、ωr；

初始化7维萤火虫数目χ，每组萤火虫种群数目w，最大吸引度β0，光强吸收系数γ，最大迭代次数m，初始化萤火虫随机位置，所有萤火虫的最大亮度位置表示各个参数值最优；萤火虫最大亮度表示能达到最大风功率和最低控制成本，目标函数为fp，当fp最小时，满足要求。

由于萤火虫会相对更亮的而吸引，当每组萤火虫都聚集在一起时，即为最大亮度。即表示符合要求。

每组萤火虫的位置范围由上述可知分别为[0,kp1]、[0,ki1]、[0,kp2]、[0,ki2]、[0,kp3]、[0,ki3]、[0,ωr]，计算fp作为各组萤火虫i的最大荧光亮度，由于萤火虫之间会向比它更亮的地方靠近，即吸引度和相对亮度

步骤4：

计算萤火虫i当前时刻的吸引度

其中：β0为最大吸引度；γ为光强吸收系数；为当前时刻萤火虫i与j之间的距离

计算萤火虫i当前时刻的相对亮度

其中：为萤火虫i当前时刻的最大荧光亮度，在数值上等于fp，萤火虫i由当前相对亮度确定移动方向，每一只萤火虫因为所处位置不同发出的荧光亮度也不同，通过比较各点的相对亮度，相对亮度高的萤火虫可以吸引相对亮度低的萤火虫向自己移动；

步骤5：

所述萤火虫i自适应移动步长如下：

式中：表示第i只萤火虫的步长,riz表示第i只萤火虫中亮度与最高亮度萤火虫z之间的距离，riz越大，步长越大，系统寻优速度就越快，反之若距离最亮的萤火虫越近，步长越小。

步骤6：更新萤火虫的位置，对处在最佳位置的萤火虫进行随机扰动，对于萤火虫的位置进行更新，所述位置更新公式如下：

其中，代表当前时刻萤火虫i所处位置，代表前一时刻萤火虫i所处位置；代表前一时刻萤火虫j所处位置；rand为[0,1]上服从均匀分布的随机因子。

根据更新后萤火虫的位置，重新计算萤火虫的最大荧光亮度；保存当前所有萤火虫个体中的最大荧光亮度及位置等信息。

步骤7：计算搜索精度当满足搜索精度要求或迭代次数t达到最大迭代次数m，转到上一步骤；否则，迭代次数加1，以更新后的萤火虫位置计算萤火虫当前时刻的吸引度和相对亮度，转步骤2，进行下一次搜索；

其中，表示第i只萤火虫当前所在位置；表示最大荧光亮度最高值的萤火虫当前所在位置；

步骤8，7组萤火虫满足精度或最大迭代次数，结束算法，给出最终求出的最优解，即找到目标函数fp最小，即成本最低和最大风功率

步骤9，输出上一阶段风速最优值后，若此时风速与差值大于3m/s，则跳到步骤二，若小于，则结束算法，直至时，重启算法。

实施例2

在又一个具体的实施例中，本发明还提供了一种永磁风机pi参数优化装置，这一装置可以执行如实施例1中所述的方法。在一个具体的实施方式中，该装置可以设置如下：

风速检测模块，用于检测当下风速；

参数优化计算模块，基于目标函数以及当下风速，通过萤火虫算法计算所述目标函数的全局最优解，输出上一阶段风速对应的最优值，从而得到优化的pi参数；

风速判断模块，获取所述上一阶段风速，并基于所述当下风速以及上一阶段风速的差值，控制所述参数优化计算模块进行计算或停止。

更为优选地，该参数优化计算模块包括正弦脉宽调制(spwm)单元，用于实现最大风功率和最低控制成本，从而得到全局最优解。

优选地，所述风速判断模块中，若当下风速与所述上一阶段风速差值大于预设值，则获取当下风速，控制所述参数优化计算模块进行计算，若小于等于预设值，则结束所述参数优化计算模块的计算。

优选地，所述目标函数为：

其中，

ωm机械转动速度，ωm^*追踪最佳参考值，α1和α2为权重，pmax为风力发电机组的最大风能捕捉，vd、vq分别为定子电压的d-q分量，为瞬时定子电流d轴分量基准值，isd为瞬时定子电流d轴分量。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read-onlymemory，rom)或随机存储记忆体(randomaccessmemory，ram)等。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。