CN104806450A - 一种基于万有引力神经网络的风电***mppt控制方法 - Google Patents
一种基于万有引力神经网络的风电***mppt控制方法 Download PDFInfo
- Publication number
- CN104806450A CN104806450A CN201510133361.8A CN201510133361A CN104806450A CN 104806450 A CN104806450 A CN 104806450A CN 201510133361 A CN201510133361 A CN 201510133361A CN 104806450 A CN104806450 A CN 104806450A
- Authority
- CN
- China
- Prior art keywords
- speed
- wind
- particle
- power
- fan
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Wind Motors (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
一种基于万有引力神经网络的风电***MPPT控制方法,所述方法在构造大量功率-转速-风速样本的基础上,建立风速的万有引力神经网络预测模型并利用该模型进行风速估计,继而通过最佳叶尖速比法预测出最大功率点所对应的最优风机转速,然后将风机的转速调节到预测的最优风机转速,并以该转速为初始值,采用占空比扰动观察法以设定的扰动步长跟踪风机的最大功率。本发明采用估计的方法获取风速,无需装设风速传感器,可有效节约***的控制成本,提高***的可靠性;该方法利用万有引力搜索算法优化神经网络模型,可有效提高风速估计的精度;此外,本发明还具有跟踪速度快的优点,可提高风机的发电效率。
Description
技术领域
本发明涉及一种基于万有引力神经网络的风力发电***最大功率跟踪控制方法,属于风力发电技术领域。
背景技术
随着社会经济的不断发展和能源危机的日益加剧,风能作为可再生能源越来越受到当今世界的重视。
风力发电机的最大功率跟踪是指在额定风速以下,通过调节风轮转速,使风能利用系数保持在最大值,获得最大功率输出。基于最大功率跟踪的基本控制原理,风力发电机的控制方法大致可以分为三类:叶尖速比法,功率反馈法以及爬山搜索法。
叶尖速比法:在某一固定的风速下,通过调节风机转速,使叶尖速比保持最佳值,就可以实现最大功率跟踪。该方法的优点在于思路比较简单,缺点在于需要通过风速传感器进行有效风速的精确测量,实际工程难度较大。另外,风速传感器的引入不仅会增加***成本,而且还会降低***的可靠性。
功率反馈法:通过转速传感器测量出风力发电机的转速ω,利用发电机电磁转矩直接控制发电机的有功功率来跟踪已知的最佳转速—功率特性曲线,从而实现***的最大功率跟踪控制。此方法的难点 在于需要掌握风机的最佳转速—功率特性曲线,而获得该曲线需要进行大量的计算和实验,过程十分复杂。
爬山搜索法:对风机转速控制的指令值以一定的步长Δω进行扰动,然后观察风车的功率变化。如果功率增加,那么风机转速的扰动方向不变,如果功率减小,则将风车转速扰动反向。爬山搜索法的优点在于不需要进行风速测量,难点在于扰动周期以及扰动步长的选择。另外,该方法要求转速对风速变化具有良好的瞬时响应特性,当风机容量较大时,由于***存在惯性,转速扰动并不能得到及时的功率响应,致使爬山搜索法的实际控制效果并不理想。
综上所述,现有的风力发电***控制方法存在成本高、跟踪准确性差等缺点,限制了风力发电的发展,因此有必要加以改进。
发明内容
本发明的目的在于针对现有技术之弊端,提供一种控制成本低、跟踪准确性高的基于万有引力神经网络的风电***MPPT控制方法。
本发明所述问题是以下述技术方案实现的:
一种基于万有引力神经网络的风电***MPPT控制方法,所述方法在构造大量功率-转速-风速样本的基础上,建立风速的万有引力神经网络预测模型并利用该模型进行风速估计,继而通过最佳叶尖速比法预测出最大功率点所对应的最优风机转速,然后将风机的转速调节到预测的最优风机转速,并以该转速为初始值,采用占空比扰动观察法以设定的扰动步长跟踪风机的最大功率。
上述基于万有引力神经网络的风电***MPPT控制方法,所述方 法包括以下步骤:
a.构造功率-转速-风速样本:
设定风速Vw以ΔV作为步长从启动风速V0,逐步增加至额定风速Vt,并在每一风速条件下,令风机转速ω以Δω为步长从起始转速ω0升至额定转速ωt,在每一风速-风机转速(Vw(i),ω(i))环境下,利用下式计算出风机机械功率Pm(i):
式中:R为风轮半径,ρ为空气密度,Cp为风能利用系数,λ为叶尖速比,β为叶片桨距角,
由此得到大量风速-风机转速-功率样本{ω(i),Pm(i),Vw(i)};
b.建立风速的万有引力神经网络预测模型:
利用BP神经网络建立风速预测模型,并利用万有引力搜索算法对BP神经网络模型进行优化,具体步骤如下:
①确定BP神经网络结构并设置初始的权值及阈值;
②设置粒子群体规模N以及各粒子初始位置 D为粒子维数,并采用随机数生成的方式初始化各粒子位置;
③计算各粒子的适应度函数值:
定义粒子的适应度函数fi为BP神经网络在训练样本上的均方差:
其中:m为输出节点个数;P为训练样本的个数;为网络期望输出值;xpj为网络实际输出值;
④更新群体中的fbest与fworst(fbest=minfj,fworst=maxfj),按下式计算各粒子的质量Mi:
⑤按照下式计算各粒子的引力合力Fi与加速度ai:
式中:t为迭代次数,为第i个粒子受到的来自其它粒子引力合力Fi的第d维分量,为第i个粒子的第d维的加速度,randi为[0,1]之间的随机数,G(t)为引力时间常数;Mpi(t)和Maj(t)分别为第i个粒子的被动惯性质量和第j个粒子的主动惯性质量,和 分别为第i个粒子与第j个粒子的第d维的位置;
⑥按照下式更新各粒子的速度vi与位置Pi:
式中:为第i个粒子的第d维的位置;为第i个粒子的第d维的速度;
⑦返回步骤③,达到最大迭代次数后停止迭代,此时fbest所对应的粒子位置即为优化后的BP神经网络初始权值与阈值;
⑧利用步骤a所构造的功率-转速-风速样本开展模型训练:
将(Pm,ω)作为万有引力神经网络模型的输入,输出为估计风速Vpre,每次迭代后,判断是否满足迭代终止条件,若满足则算法终止,否则继续迭代;
c.采集风机的功率信号与转速信号,利用训练好的万有引力神经网络预测模型进行风速估计;
d.利用估计风速,采用最佳叶尖速比法,得出所估计风速下的风机最佳转速ωref:
其中,Vpre为估计风速;
e.通过PI控制方法,将风机转速调节到最佳转速ωref;
f.以最佳转速ωref为初始值,采用占空比扰动观察法进行最大功率跟踪,当前后两次扰动的功率差值大于设定的阈值φ时,说明风速发生了突变,重复步骤c~f。
上述基于万有引力神经网络的风电***MPPT控制方法,考虑实际工程中风机转速信号与功率信号的测量存在噪声,对每个风速-风 机转速-功率样本{ω(i),Pm(i),Vw(i)}样本中的ω(i)与Pm(i)均叠加一定的高斯噪声,叠加高斯噪声后的转速ω′(i)和功率Pm′(i)如下式所示:
ω′(i)=ω(i)+δω
Pm′(i)=Pm(i)+δPm
其中,δω为转速信号噪声,服从N(0,σω 2)的高斯分布,δPm为功率信号噪声,服从N(0,σP 2)的高斯分布。
上述基于万有引力神经网络的风电***MPPT控制方法,风机功率信号的采集方法是:测量风力发电机A相与B相之间的线电压VAB,并检测其过零点的频率f,设电机的极对数为n,则风机转速信号为ω=πf/n;风机转速信号的采集方法是:测量整流器输出电压Vz和输出电流Iz,则功率信号为Pm=Iz*Vz。
本发明的有益效果在于:
1)采用估计的方法获取风速,无需装设风速传感器,可有效节约***的控制成本,提高***的可靠性;
2)利用万有引力搜索算法优化神经网络模型,可有效提高风速估计的精度;
3)当外界环境发生变化时,借助于风速估计模型可以直接将工作转速调节至最佳风机转速的附近,避免了传统爬山搜索算法逐步试探的过程,具有跟踪速度快的优点,从而提高了风机的发电效率。
附图说明
下面结合附图对本发明作进一步详述。
图1为本发明所采用的硬件结构图;
图2为万有引力神经网络预测模型估计风速的示意图;
图3为万有引力神经网络预测模型训练流程图;
图4为叶尖速比与风能利用系数关系图;
图5为占空比扰动观察法算法流程图。
图中各标号为:1-永磁直驱风力发电机;2-第一电压传感器;3-整流器;4-第二电压传感器;5-第一电容;6-第二电容;7-电流传感器。
图中和文中各符号为:为第i个粒子受到的来自其它粒子引力合力Fi的第d维分量,Vw为设定风速,V0为启动风速,Vt为额定风速,ω为风机转速,ω0为起始转速,ωt为额定转速,Pm(i)为风机机械功率,R为风轮半径,ρ为空气密度,Cp为风能利用系数,λ为叶尖速比,β为叶片桨距角,wlh为输入层神经元与隐含层神经元之间的连接权值集合,whm为隐含层神经元与输出层神经元之间的连接权值集合,bh和bm分别表示隐含层神经元与输出层神经元的阈值向量,fi为粒子的适应度函数,m为输出节点个数;P为训练样本的个数;为网络期望输出值;xpj为网络实际输出值;Mi为粒子的质量,t为迭代次数,为第i个粒子的第d维的加速度,randi为[0,1]之间的随机数,G(t)为引力时间常数;Mpi(t)和Maj(t)分别为第i个粒子的被动惯性质量和第j个粒子的主动惯性质量;为第i个粒子的第d维的位置;为第i个粒子的第d维的速度,ωref为风机最佳转速,Vpre为估计风速。
具体实施方式
参考图1,本发明所采用的硬件结构包括MPPT控制器、整流器、第一电压传感器、第二电压传感器、电流传感器、DC-DC变换器、驱动模块、第一电容和第二电容。该装置结合风机与负载即可构成永磁直驱风力发电***。其中,整流器的三相输入端与风机的三相输出端连接,单相输出正端与第一电容的正极连接,单相输出负端接地;第一电压传感器的测量正端与风机的A相输出端连接,测量负端与风机的B相输出端连接,测量输出端与MPPT控制器连接;第一电容负极接地;第二电压传感器的测量正端与第一电容正极连接,测量负端接地,测量输出端与MPPT控制器连接;电流传感器的待测电流输入端与电压传感器正极连接,待测电流输出端与DC-DC变换器的输入端连接,测量信号输出端与MPPT控制器连接;DC-DC变换器的输入正端与电流传感器的待测电流输出端连接,输入负端接地,脉宽调制信号输入端与驱动模块一端连接,输出正端与第二电容的正极连接,输出负端接地;第二电容负极接地;驱动模块另外一端与MPPT控制器连接;MPPT控制器由DSP与FPGA两部分构成,其中DSP的端口1至端口11与FPGA的端口1至端口11对应连接。负载一端与第二电容正极连接,另外一端接地。
本发明在构造大量功率-转速-风速样本的基础上,建立风速的万有引力神经网络预测模型并利用该模型进行风速估计,继而通过最佳叶尖速比法预测出最大功率点所对应的最优风机转速,然后将风机的转速调节到预测的最优风机转速,并以该转速为初始值,采用占空比 扰动观察法以设定的扰动步长跟踪风机的最大功率。
本方法主要包括以下步骤:
1.构造风速估计模型训练样本;
根据风力机的空气动力学特性,可以得出风力机输出机械功率为:
式中:R为风轮半径,m;ρ为空气密度,Vw为风速,m/s;Cp为风能利用系数,其大小与叶尖速比λ以及叶片桨距角β有关,可近似表示为:
其中,叶尖速比λ为风机叶片的叶尖速度与风速的比值,有
可知,在风轮叶片大小、风速和空气密度一定时,影响功率输出的唯一因素是风能利用系数Cp。对于永磁直驱风力机,风速在额定风速以下时,一般采用不变桨距方式,即β=0,则可得出风能利用系数Cp与叶尖速比λ之间的关系曲线,此时,存在一个最佳的叶尖速比λopt(约为6.3),对应最大的风能利用系数Cpmax。在不同的风速下,通过调节风机转速ω使叶尖速比λ始终在工作在λopt,即可实现最大功率跟踪。
本发明的训练样本构造过程如下:
设定风速Vw以ΔV作为步长从启动风速V0,逐步增加至额定风速Vt,并在每一风速条件下,令风机转速ω以Δω为步长从起始转速ω0升至额定转速ωt,并根据式(1)(2)(3)计算出每一风速-风机转速(Vw(i),ω(i))环境下的风机机械功率Pm(i),从而可构成大量风速-风机转速-功率样本{ω(i),Pm(i),Vw(i)},即在已知风机转速ω(i)和功率Pm(i)的情况下,可得出风速值为Vw(i)。考虑实际工程中风机转速信号与功率信号的测量存在噪声,对每个样本中的ω(i)与Pm(i)均叠加一定的高斯噪声,如下式所示:
ω′(i)=ω(i)+δω (4)
Pm′(i)=Pm(i)+δPm (5)
其中,δω服从N(0,σω 2)的高斯分布,δPm服从N(0,σP 2)的高斯分布。
2.基于万有引力神经网络建立风速估计模型;
利用BP神经网络建立风速估计模型,利用步骤1所收集的训练样本开展模型训练。传统的BP训练算法采用随机数生成的方式设置初始权值以及阈值,具有容易陷入局部最优的缺点。本发明提出利用万有引力搜索算法对BP神经网络的初始权值以及阈值进行优化,从而保证训练结果的全局最优,得到高精度的风速估计模型。
万有引力搜索算法是Rashedi等人首先提出的一种基于万有引力定律进行寻优的智能优化方法。该方法将优化问题的解视为一组在空间运行的粒子,粒子之间通过万有引力作用相互吸引,粒子运动遵 循动力学规律,万有引力的作用使得粒子朝着质量最大的粒子移动,而质量最大的粒子占据最优位置,从而可求出优化问题的最优解。算法通过个体间的万有引力相互作用实现优化信息的共享,引导群体向最优解区域搜索。
该算法的原理为:
在D维的搜索空间中,设有N个粒子,定义第i个粒子位置为 第t次迭代时,粒子j对i的万有引力定义为
式中:G(t)为引力时间常数;Mpi(t)和Maj(t)分别为粒子i的被动惯性质量和粒子j的主动惯性质量。
假设引力质量与惯性质量相等,根据每个粒子的适应度函数值计算粒子惯性质量:
Mai=Mpi=Mi (7)
对于极小化问题,fbest=minfj,fworst=maxfj。
对于第i个粒子,受到来自其它粒子引力的合力可用引力的随机加权和表示为:
基于牛顿第二定律,粒子i产生的加速度为:
则粒子的速度和位置可以更新为:
式中:为第i个粒子的第d维的位置;为第i个粒子的第d维的速度;第i个粒子的第d维的加速度;randi为[0,1]之间的随机数。
本发明利用万有引力搜索算法对BP神经网络模型进行优化的具体步骤如下:
a)确定BP神经网络结构并设置初始的权值及阈值;
设BP神经网络为三层结构,输入层有l个神经元,隐含层有h个神经元,输出层有m个神经元,wlh={wji|j=1,…,l,i=1,…,h}为输入层神经元与隐含层神经元之间的连接权值集合,whm={wkj|k=1,…,m,j=1,…,h}为隐含层神经元与输出层神经元之间的连接权值集合,bh和bm分别表示隐含层神经元与输出层神经元的阈值向量。
b)设置粒子群体规模以及各粒子初始位置;
设置群体规模为N,各粒子位置为Pi=[wlh,bh,whm,bm],并采用随机数生成的方式初始化各粒子位置。
c)计算各粒子适应度函数值fi;
定义粒子的适应度函数fi为BP神经网络在训练样本上的均方差:
其中:m-输出节点个数;P-训练样本的个数;-网络期望输出值;xpj-网络实际输出值。
d)更新群体中的fbest与fworst,按照式(9)计算各粒子的质量Mi;
e)按照公式(10)(11)计算各粒子的引力合力Fi与加速度ai;
f)按照公式(12)(13)更新各粒子的速度vi与位置Pi;
g)返回步骤c),达到最大迭代次数后停止迭代,此时fbest所对应的粒子位置即为优化后的BP神经网络初始权值与阈值;
h)应用BP训练算法对网络进行训练,即迭代计算隐含层与输出层输出,以及隐含层与输出层的误差,并根据误差对权值与阈值更新;每次迭代后,判断是否满足迭代终止条件,若满足则算法终止,否则继续迭代。
3.通过采集功率信号与转速信号,利用万有引力神经网络模型进行风速估计;
利用第一电压传感器采集风力发电机A相与B相的线电压VAB,并检测其过零点的频率f,设电机的极对数为n,则可获取风机转速信号ω=πf/n;利用第二电压传感器采集整流器输出电压Vz,利用电流传感器采集整流器输出电流Iz,则近似有功率信号Pm=Iz*Vz,将(Pm,ω)作为万有引力神经网络模型输入,输出即为估计风速Vpre。
4.利用最佳叶尖速比法,得出所估计风速下的风机最佳转速ωref:
其中,Vpre为估计风速,λ为叶尖速比;R为叶轮半径。
5.通过PI控制方法,调节风机转速到ωref;
6.以ωref为初始值,采用占空比扰动观察法进行最大功率跟踪。当前后两次扰动的功率差值大于设定的阈值φ时,说明风速发生了突变,重复步骤3~6。
下面结合附图,对优选实施例作详细说明。应该强调的是下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
图1为本发明硬件结构图。实施例中,小型永磁直驱风力发电机2的主要参数为:风轮直径为1.3m,额定功率为300W,额定电压为24V,额定转速800r/min,启动风速1m/s,额定风速10m/s,安全风速25m/s;MPPT控制器中DSP采用dsPIC33FJ06GS101单片机,其主要作用为进行信号采集与控制信号的输出;FPGA采用cycloneIV EP4CE22F17,其主要作用为进行神经网络风速估计模型的运算,并通过端口1至端口11与dsPIC33FJ06GS101单片机的端口1至端口11对应连接;DC-DC变换器采用CUK电路,驱动模块选用MCP14E3,第一电压传感器和第二电压传感器选用LV28-P,电流传感器选用LA25-NP,第一电容C1=10uF,第二电容C2=100uF。
图2为万有引力神经网络预测模型估计风速的示意图;
图3为万有引力神经网络预测模型训练流程图;
所述永磁直所述驱风力发电***最大功率跟踪控制方法,主要包 括以下步骤:
1.构造风速估计模型训练样本;
实施例中,取风轮半径R=0.65m,空气密度ρ=1.205kg/m3,桨距角β=0,设定风速Vw以ΔV=0.1m/s作为步长从启动风速V0=1m/s逐步增加至额定风速Vt=10m/s,并在每一风速条件下,令风机转速ω以Δω=20rad/s为步长从ω0=1rad/s升至额定转速ωt=84rad/s,并根据式(1)(2)(3)计算出每一风速-风机转速(Vw(i),ω(i))环境下的风机机械功率Pm(i),从而可构成P=5915个风速-风机转速-功率样本{ω(i),Pm(i),Vw(i)},并分别按照式(4)和(5)对每个样本中的ω(i)与Pm(i)均叠加一定的高斯噪声,δω服从N(0,σω 2)的高斯分布,δPm服从N(0,σP 2)的高斯分布,实施例中,取σω=5,σP=20。
2.基于万有引力神经网络建立风速估计模型;
本发明利用万有引力搜索算法对BP神经网络模型进行优化的具体步骤如下:
a)确定BP神经网络结构并设置初始的权值及阈值;
实施例中,采用三层BP神经网络,其中输入层神经元节点数为2,隐含层神经元节点数为7,输出层节点数为1,神经元转移函数采用tansig函数。由于所建立的BP网络为2-7-1结构,故l=2,h=7,m=1,则权值和阈值总数为29个,对应优化问题的维数为29;
b)初始化群体。设置粒子群体规模以及各粒子初始位置;
实施例中,设置粒子群体的规模为30个,并采用在[0,1]之间 生成随机数的方式初始化各粒子的位置,取最大迭代次数为400次。
c)依据公式(14)计算各粒子适应度函数值fi;
d)更新群体中的fbest与fworst,按照式(9)计算各粒子的质量Mi;
e)按照公式(10)(11)计算各粒子的引力合力Fi与加速度ai;
f)按照公式(12)(13)更新各粒子的速度vi与位置Pi;
g)返回c,达到最大迭代次数后停止迭代,此时fbest所对应的粒子位置即为优化后的BP神经网络初始权值与阈值;
实施例中,优化后输入层到隐含层的权值和阈值为:
优化后隐含层到输出层的权值和阈值为:
ω2=(0.0318,0.2769,0.0462,0.0971,0.8235,0.6948,0.3171),
b2=0.8502。
h)应用BP训练算法对网络进行训练,即迭代计算隐含层与输出层输出,以及隐含层与输出层的误差,并根据误差对权值与阈值更新;每次迭代后,判断是否满足迭代终止条件,若满足则算法终止,否则继续迭代。
实施例中,设置迭代终止条件为在训练样本上的均方误差MSE<0.03,可得到全局最优下的输入层到隐含层的权值和阈值为:
全局最优下的隐含层到输出层的权值和阈值为:
ω2=(0.8687,0.0844,0.3998,0.2599,0.8001,0.4314,0.9106),
b2=0.8693
3.通过采集功率信号与转速信号,利用万有引力神经网络模型进行风速估计;
实施例中,利用第一电压传感器、第二电压传感器和电流传感器所采集信号计算出功率信号Pm与转速信号ω,电机的极对数n=6;某一工况下得出ω=77.5rad/s,Pm=55.1W,将(P,ω)=(55.1,77.5)作为万有引力神经网络模型输入,可估计出风速为Vpre=5.80m/s;
4.利用最佳叶尖速比法,得出所估计风速下的风机最佳转速ωref:
实施例中,取λ=6.7,
5.通过PI控制方法(设定Kp=0.02,KI=0.1),调节风机转速到ωref;
6.以ωref为初始值,占空比扰动观察法进行最大功率跟踪(扰动周期T=2s,步长Δd=0.02),经过连续2拍的负方向扰动后,进入稳态过程,跟踪时间为4s,Pm平均值为69.8W。当前后两次扰动的功率差值大于设定的阈值φ时,说明风速发生了突变,重复步骤3~6。实施例中,取φ=15W。
为验证本发明的功效,在上述同样工况下(ω=77.5rad/s,Pm=55.1W),直接采用占空比扰动观察法进行最大功率跟踪(扰动周期T=2s,扰动步长Δd=0.02),则需经过连续17拍的负方向扰动后进入稳态过程,跟踪时间为34s,可知,本发明所述方法可以明显的缩短跟踪时间;若取扰动步长Δd=0.05,则需经过连续7拍的负方向扰动后进入稳态过程,跟踪时间为14s,但此时Pm存在较大波动,其平均值为66.5W,与本发明所述方法相比,存在一定的功率损失。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可推想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (4)
1.一种基于万有引力神经网络的风电***MPPT控制方法,其特征是,所述方法在构造大量功率-转速-风速样本的基础上,建立风速的万有引力神经网络预测模型并利用该模型进行风速估计,继而通过最佳叶尖速比法预测出最大功率点所对应的最优风机转速,然后将风机的转速调节到预测的最优风机转速,并以该转速为初始值,采用占空比扰动观察法以设定的扰动步长跟踪风机的最大功率。
2.根据权利要求1所述的基于万有引力神经网络的风电***MPPT控制方法,其特征是,所述方法包括以下步骤:
a.构造功率-转速-风速样本:
设定风速Vw以ΔV作为步长从启动风速V0,逐步增加至额定风速Vt,并在每一风速条件下,令风机转速ω以Δω为步长从起始转速ω0升至额定转速ωt,在每一风速-风机转速(Vw(i),ω(i))环境下,利用下式计算出风机机械功率Pm(i):
式中:R为风轮半径,ρ为空气密度,Cp为风能利用系数,λ为叶尖速比,β为叶片桨距角,
从而得到大量风速-风机转速-功率样本{ω(i),Pm(i),Vw(i)};
b.建立风速的万有引力神经网络预测模型:
利用BP神经网络建立风速预测模型,并利用万有引力搜索算法对BP神经网络模型进行优化,具体步骤如下:
①确定BP神经网络结构并设置初始的权值及阈值;
②设置粒子群体规模N以及各粒子初始位置,i=1,2,3,…,N,D为粒子维数,并采用随机数生成的方式初始化各粒子位置;
③计算各粒子的适应度函数值:
定义粒子的适应度函数fi为BP神经网络在训练样本上的均方差:
其中:m为输出节点个数;P为训练样本的个数;为网络期望输出值;xpj为网络实际输出值;
④更新群体中的fbest与fworst(fbest=minfj,fworst=maxfj),按下式计算各粒子的质量Mi:
⑤按照下式计算各粒子的引力合力Fi与加速度ai:
式中:t为迭代次数,为第i个粒子受到的来自其它粒子引力合力Fi的第d维分量,为第i个粒子的第d维的加速度,randi为[0,1]之间的随机数,G(t)为引力时间常数;Mpi(t)和Maj(t)分别为第i个粒子的被动惯性质量和第j个粒子的主动惯性质量,和分别为第i个粒子与第j个粒子的第d维的位置;
⑥按照下式更新各粒子的速度vi与位置Pi:
式中:为第i个粒子的第d维的位置;为第i个粒子的第d维的速度;
⑦返回步骤③,达到最大迭代次数后停止迭代,此时fbest所对应的粒子位置即为优化后的BP神经网络初始权值与阈值;
⑧利用步骤a所构造的功率-转速-风速样本开展模型训练:
将(Pm,ω)作为万有引力神经网络模型的输入,输出为估计风速Vpre,每次迭代后,判断是否满足迭代终止条件,若满足则算法终止,否则继续迭代;
c.采集风机的功率信号与转速信号,利用训练好的万有引力神经网络预测模型进行风速估计;
d.利用估计风速,采用最佳叶尖速比法,得出所估计风速下的风机最佳转速ωref:
其中,Vpre为估计风速;
e.通过PI控制方法,将风机转速调节到最佳转速ωref;
f.以最佳转速ωref为初始值,采用占空比扰动观察法进行最大功率跟踪,当前后两次扰动的功率差值大于设定的阈值φ时,说明风速发生了突变,重复步骤c~f。
3.根据权利要求2所述的基于万有引力神经网络的风电***MPPT控制方法,其特征是,考虑实际工程中风机转速信号与功率信号的测量存在噪声,对每个风速-风机转速-功率样本{ω(i),Pm(i),Vw(i)}样本中的ω(i)与Pm(i)均叠加一定的高斯噪声,叠加高斯噪声后的转速ω′(i)和功率Pm′(i)如下式所示:
ω′(i)=ω(i)+δω
Pm′(i)=Pm(i)+δPm
其中,δω为转速信号噪声,服从N(0,σω 2)的高斯分布,δPm为功率信号噪声,服从N(0,σP 2)的高斯分布。
4.根据权利要求3所述的基于万有引力神经网络的风电***MPPT控制方法,其特征是,风机功率信号的采集方法是:测量风力发电机A相与B相之间的线电压VAB,并检测其过零点的频率f,设电机的极对数为n,则风机转速信号为ω=πf/n;风机转速信号的采集方法是:测量整流器输出电压Vz和输出电流Iz,则功率信号为Pm=Iz*Vz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510133361.8A CN104806450B (zh) | 2015-03-25 | 2015-03-25 | 一种基于万有引力神经网络的风电***mppt控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510133361.8A CN104806450B (zh) | 2015-03-25 | 2015-03-25 | 一种基于万有引力神经网络的风电***mppt控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104806450A true CN104806450A (zh) | 2015-07-29 |
CN104806450B CN104806450B (zh) | 2017-07-14 |
Family
ID=53691579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510133361.8A Expired - Fee Related CN104806450B (zh) | 2015-03-25 | 2015-03-25 | 一种基于万有引力神经网络的风电***mppt控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104806450B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105715452A (zh) * | 2014-12-18 | 2016-06-29 | 通用电气公司 | 用于增加风力涡轮机功率输出的***和方法 |
CN106096246A (zh) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 |
CN106194582A (zh) * | 2016-09-19 | 2016-12-07 | 华能新能源股份有限公司辽宁分公司 | 基于风速测量与估计的风电***mppt控制装置及方法 |
CN106372718A (zh) * | 2016-08-29 | 2017-02-01 | 华北电力大学(保定) | 雾霾积灰条件下的光伏发电功率输出减少率估计方法 |
CN107272815A (zh) * | 2017-07-24 | 2017-10-20 | 广东电网有限责任公司电力调度控制中心 | 一种基于sqp信赖域算法的光伏单峰mppt方法 |
CN107923368A (zh) * | 2015-12-23 | 2018-04-17 | Vdm-泰克尼卡有限公司 | 调整风力涡轮机取力器的方法 |
CN108347072A (zh) * | 2017-01-24 | 2018-07-31 | 罗伯特·博世有限公司 | 用于运行混合能量存储器***的方法 |
CN108547736A (zh) * | 2018-03-27 | 2018-09-18 | 中南大学 | 风速风向预测方法及风力发电机组的偏航控制方法 |
CN108594895A (zh) * | 2018-05-31 | 2018-09-28 | 江苏大学 | 基于温差发电的汽车尾气能量回收装置电控***及方法 |
CN108843495A (zh) * | 2018-05-30 | 2018-11-20 | 上海电机学院 | 小型风力发电***mppt控制方法 |
CN110108457A (zh) * | 2019-04-18 | 2019-08-09 | 大唐东北电力试验研究院有限公司 | 一种基于万有引力神经网络的一次风机轴温诊断方法 |
CN110985290A (zh) * | 2019-12-04 | 2020-04-10 | 浙江大学 | 一种基于支持向量回归的最优转矩控制方法 |
CN111305999A (zh) * | 2020-04-01 | 2020-06-19 | 宁夏三聚阳光新能源科技有限公司 | 一种智能调整风力发电*** |
CN112083753A (zh) * | 2020-09-28 | 2020-12-15 | 东莞市钜大电子有限公司 | 一种光伏并网逆变器的最大功率点跟踪控制方法 |
CN114386718A (zh) * | 2022-03-16 | 2022-04-22 | 广州兆和电力技术有限公司 | 一种结合粒子群神经网络的风电场输出功率短时预测算法 |
CN115268559A (zh) * | 2022-05-23 | 2022-11-01 | 北京华能新锐控制技术有限公司 | 用于永磁同步风力发电机的最大功率点跟踪鲁棒控制方法 |
CN116345967A (zh) * | 2022-10-31 | 2023-06-27 | 广东电网有限责任公司江门供电局 | 一种波浪发电***的功率控制方法、***和设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009091923A (ja) * | 2007-10-04 | 2009-04-30 | Univ Of Ryukyus | 風力発電装置の最大電力追従制御装置 |
CN102352814A (zh) * | 2011-06-21 | 2012-02-15 | 北京交通大学 | 一种大型直驱风电机组的最大功率跟踪方法 |
CN102522937A (zh) * | 2011-12-09 | 2012-06-27 | 北京交通大学 | 一种风力发电最大功率跟踪控制方法 |
US20120205913A1 (en) * | 2009-10-23 | 2012-08-16 | Juan Carlos Garcia Andujar | Wind turbine control methods for improving the production of energy |
CN103410662A (zh) * | 2013-08-06 | 2013-11-27 | 江苏科技大学 | 风电***最大风能捕获神经网络补偿控制方法 |
CN103437955A (zh) * | 2013-08-13 | 2013-12-11 | 华北电力大学(保定) | 小型永磁直驱风力发电***最大功率跟踪装置及控制方法 |
CN103835878A (zh) * | 2013-04-07 | 2014-06-04 | 南京理工大学 | 基于神经网络优化起始转速的最大功率点跟踪控制方法 |
-
2015
- 2015-03-25 CN CN201510133361.8A patent/CN104806450B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009091923A (ja) * | 2007-10-04 | 2009-04-30 | Univ Of Ryukyus | 風力発電装置の最大電力追従制御装置 |
US20120205913A1 (en) * | 2009-10-23 | 2012-08-16 | Juan Carlos Garcia Andujar | Wind turbine control methods for improving the production of energy |
CN102352814A (zh) * | 2011-06-21 | 2012-02-15 | 北京交通大学 | 一种大型直驱风电机组的最大功率跟踪方法 |
CN102522937A (zh) * | 2011-12-09 | 2012-06-27 | 北京交通大学 | 一种风力发电最大功率跟踪控制方法 |
CN103835878A (zh) * | 2013-04-07 | 2014-06-04 | 南京理工大学 | 基于神经网络优化起始转速的最大功率点跟踪控制方法 |
CN103410662A (zh) * | 2013-08-06 | 2013-11-27 | 江苏科技大学 | 风电***最大风能捕获神经网络补偿控制方法 |
CN103437955A (zh) * | 2013-08-13 | 2013-12-11 | 华北电力大学(保定) | 小型永磁直驱风力发电***最大功率跟踪装置及控制方法 |
Non-Patent Citations (1)
Title |
---|
李满天,王春林,王鹏飞,查富生,郭伟: "一种基于万有引力的自生长神经网络算法", 《机械与电子》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105715452A (zh) * | 2014-12-18 | 2016-06-29 | 通用电气公司 | 用于增加风力涡轮机功率输出的***和方法 |
CN107923368A (zh) * | 2015-12-23 | 2018-04-17 | Vdm-泰克尼卡有限公司 | 调整风力涡轮机取力器的方法 |
CN107923368B (zh) * | 2015-12-23 | 2020-03-13 | Vdm-泰克尼卡有限公司 | 调整风力涡轮机取力器的方法 |
CN106096246A (zh) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 |
CN106096246B (zh) * | 2016-06-06 | 2018-08-31 | 华北电力大学(保定) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 |
CN106372718B (zh) * | 2016-08-29 | 2018-07-31 | 华北电力大学(保定) | 雾霾积灰条件下的光伏发电功率输出减少率估计方法 |
CN106372718A (zh) * | 2016-08-29 | 2017-02-01 | 华北电力大学(保定) | 雾霾积灰条件下的光伏发电功率输出减少率估计方法 |
CN106194582B (zh) * | 2016-09-19 | 2018-09-04 | 华能新能源股份有限公司辽宁分公司 | 基于风速测量与估计的风电***mppt控制装置及方法 |
CN106194582A (zh) * | 2016-09-19 | 2016-12-07 | 华能新能源股份有限公司辽宁分公司 | 基于风速测量与估计的风电***mppt控制装置及方法 |
CN108347072A (zh) * | 2017-01-24 | 2018-07-31 | 罗伯特·博世有限公司 | 用于运行混合能量存储器***的方法 |
CN108347072B (zh) * | 2017-01-24 | 2023-05-26 | 罗伯特·博世有限公司 | 用于运行混合能量存储器***的方法 |
CN107272815A (zh) * | 2017-07-24 | 2017-10-20 | 广东电网有限责任公司电力调度控制中心 | 一种基于sqp信赖域算法的光伏单峰mppt方法 |
CN108547736A (zh) * | 2018-03-27 | 2018-09-18 | 中南大学 | 风速风向预测方法及风力发电机组的偏航控制方法 |
CN108843495A (zh) * | 2018-05-30 | 2018-11-20 | 上海电机学院 | 小型风力发电***mppt控制方法 |
CN108594895B (zh) * | 2018-05-31 | 2023-09-26 | 江苏大学 | 基于温差发电的汽车尾气能量回收装置电控***及方法 |
CN108594895A (zh) * | 2018-05-31 | 2018-09-28 | 江苏大学 | 基于温差发电的汽车尾气能量回收装置电控***及方法 |
CN110108457B (zh) * | 2019-04-18 | 2021-11-09 | 大唐东北电力试验研究院有限公司 | 一种基于万有引力神经网络的一次风机轴温诊断方法 |
CN110108457A (zh) * | 2019-04-18 | 2019-08-09 | 大唐东北电力试验研究院有限公司 | 一种基于万有引力神经网络的一次风机轴温诊断方法 |
CN110985290B (zh) * | 2019-12-04 | 2022-02-11 | 浙江大学 | 一种基于支持向量回归的最优转矩控制方法 |
CN110985290A (zh) * | 2019-12-04 | 2020-04-10 | 浙江大学 | 一种基于支持向量回归的最优转矩控制方法 |
CN111305999B (zh) * | 2020-04-01 | 2020-12-01 | 宁夏三聚阳光新能源科技有限公司 | 一种智能调整风力发电*** |
CN111305999A (zh) * | 2020-04-01 | 2020-06-19 | 宁夏三聚阳光新能源科技有限公司 | 一种智能调整风力发电*** |
CN112083753A (zh) * | 2020-09-28 | 2020-12-15 | 东莞市钜大电子有限公司 | 一种光伏并网逆变器的最大功率点跟踪控制方法 |
CN114386718A (zh) * | 2022-03-16 | 2022-04-22 | 广州兆和电力技术有限公司 | 一种结合粒子群神经网络的风电场输出功率短时预测算法 |
CN115268559A (zh) * | 2022-05-23 | 2022-11-01 | 北京华能新锐控制技术有限公司 | 用于永磁同步风力发电机的最大功率点跟踪鲁棒控制方法 |
CN116345967A (zh) * | 2022-10-31 | 2023-06-27 | 广东电网有限责任公司江门供电局 | 一种波浪发电***的功率控制方法、***和设备 |
CN116345967B (zh) * | 2022-10-31 | 2023-08-04 | 广东电网有限责任公司江门供电局 | 一种波浪发电***的功率控制方法、***和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN104806450B (zh) | 2017-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104806450A (zh) | 一种基于万有引力神经网络的风电***mppt控制方法 | |
Tafticht et al. | Output power maximization of a permanent magnet synchronous generator based stand-alone wind turbine | |
Mousavi G | An autonomous hybrid energy system of wind/tidal/microturbine/battery storage | |
CN103437955B (zh) | 小型永磁直驱风力发电***最大功率跟踪装置及控制方法 | |
CN109642536A (zh) | 控制或处理***和方法 | |
CN101272121B (zh) | 一种风力发电机组的最大功率点跟踪方法 | |
CN102419394B (zh) | 一种预测分辨率可变的风光功率预测方法 | |
Mozafarpoor-Khoshrodi et al. | Improvement of perturb and observe method for maximum power point tracking in wind energy conversion system using fuzzy controller | |
Mohammadi et al. | Using a new wind turbine emulator to analyze tower shadow and yaw error effects | |
CN102013702B (zh) | 外部电网故障时并网风电场的动态等值方法 | |
CN106194582A (zh) | 基于风速测量与估计的风电***mppt控制装置及方法 | |
CN104454346A (zh) | 一种小型永磁直驱风力发电***最大功率跟踪控制方法 | |
Wang et al. | A hybrid, self-adapting drag-lift conversion wind energy harvesting system for railway turnout monitoring on the Tibetan Plateau | |
CN103410662A (zh) | 风电***最大风能捕获神经网络补偿控制方法 | |
CN110206686A (zh) | 一种用于风力发电机组的自适应最大功率跟踪控制方法 | |
CN103410659A (zh) | 基于High-Gain观测器的风力发电机组有效风速获取方法 | |
Ren et al. | Control strategy of maximum wind energy capture of direct-drive wind turbine generator based on neural-network | |
CN103986189B (zh) | 一种储能型双馈风电场简化模型建模方法 | |
Hossain | Application of wind energy into the transportation sector | |
Lu et al. | A dynamic wind turbine simulator of the wind turbine generator system | |
Liang et al. | The modeling and numerical simulations of wind turbine generation system with free vortex method and simulink | |
Iqbal et al. | A novel vertical axis wind turbine for energy harvesting on the highways | |
CN102156044A (zh) | 适用于直驱型风电机组测试的风力机模拟器及其选型方法 | |
CN109653946A (zh) | 风力发电***最大功率的捕获方法、***及其实现方法 | |
CN103375332B (zh) | 变速变桨风力发电机组最优阻力矩动态优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170714 Termination date: 20190325 |
|
CF01 | Termination of patent right due to non-payment of annual fee |