CN104929863A - 一种基于改进爬山算法的风力发电最大功率点跟踪方法 - Google Patents
一种基于改进爬山算法的风力发电最大功率点跟踪方法 Download PDFInfo
- Publication number
- CN104929863A CN104929863A CN201510209324.0A CN201510209324A CN104929863A CN 104929863 A CN104929863 A CN 104929863A CN 201510209324 A CN201510209324 A CN 201510209324A CN 104929863 A CN104929863 A CN 104929863A
- Authority
- CN
- China
- Prior art keywords
- search
- influence coefficient
- size
- performance number
- maximum power
- 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
- 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
本发明公开了一种基于改进爬山算法的风力发电最大功率点跟踪方法,包括:测量风机当前的功率值;对风机的转速增加一个搜索步长;测量风机转速变化后的功率变化值;根据功率值,计算功率值影响系数;根据搜索步长和功率变化值,计算得出斜率影响系数;根据功率值影响系数、斜率影响系数、本次的搜索步长和功率变化值,计算下一次的搜索步长 ;根据预设的条件判断是否结束搜索。本发明通过引入功率值影响系数和斜率影响系数,能根据当前功率值变化的实际情况改变下一次搜索步长,从而实现加快低风速时搜索最大功率点搜索的速度的目的,并且在高风速时,避免了风机在最大功率点的振荡情况。本发明可广泛应用于风力发电控制器中。
Description
技术领域
本发明涉及风力发电领域,尤其涉及一种基于改进爬山算法的风力发电最大功率点跟踪方法。
背景技术
随着传统化石能源造成的环境压力越来越大,风力发电作为一种新能源得到了广泛的重视。风力发电受风速影响比较大,不同的风速条件下,风力发电机有不同的最大功率点,对应着不同的运行转速。风速具有随机性,要求风力发电机能够实现最大功率点跟踪(MPPT),自动调整运行转速,保证风力发电机始终工作在最大功率点,最大程度上利用风能。
风力发电最大功率点跟踪的方法,最容易实现的是爬山算法,因其无需测量风速和风轮转速,成本低,特别适用于小型风力发电设备。爬山算法利用了风力发电机功率-转速曲线的单峰值特性来搜索最大功率点。该方法首先给转速一个扰动,接着观察输出功率的变化,判断功率-转速曲线在当前转速位置的斜率。如果斜率为正,则风力发电机目前工作在最大功率点的左侧,即风机运行在功率-转速曲线的左半段上,当前转速小于最佳转速,然后沿着转速增加的方向继续对风机的转速添加扰动,进行下一步的搜索,直至斜率接近零,此时风机位于功率-转速曲线的峰值的位置(即功率-转速曲线的最大功率点)。反之亦然。这种方法使风机由功率-转速曲线的两侧爬向曲线中间的峰值位置,因此形象的成为爬山算法。
专利《基于Boost类变换器的风力发电系统最大功率点跟踪方法》(专利号:201410107195.X)公开了一种改进的爬山算法,该方法在给风力发电机转速添加扰动的时候,扰动的大小——也就是搜索步长(即扰动前后转速变化的差值)不是固定的,它是功率-转速曲线相应的点的斜率的函数,斜率的绝对值越大,搜索步长越大。由于基本上风机功率-转速曲线上,离峰值越远的点,斜率绝对值越大,因此该方法可以使风机在功率-转速曲线上远离最大功率点时可以快速改变转速,加快爬向峰值的速度,而在峰值附近,斜率绝对值变小,搜索步长变小,避免了风机在最大功率点处振荡。此种方法亦称为变步长爬山算法。
但是,当在不同风速条件下,风机的功率-转速曲线上,远离峰值的点的斜率并不一样。如果使用上述专利公开的方法,由不同风速的功率-转速曲线可知,当风速较大时,远离最大功率点的位置斜率较大,容易获得较大的搜索步长。而当风速较低时,远离最大功率点的位置斜率较小,导致搜索步长小,搜索速度慢。单靠增大斜率对搜索步长的影响系数,使同样的斜率值下搜索步长更大,虽然可以加快低风速时的爬山速度,但由于在高风速时,最大功率点附近的斜率依然比较大,容易导致风机在最大功率点附近振荡,例如风机工作在最大功率点左侧附近时,由于计算得到的搜索步长过大,导致下一步搜索时转速过大,使风机到达了最大功率点的右侧,然后又由于同样的原因,再下一次搜索时,风机又回到了最大功率点的左侧,如此往复,导致转速忽大忽小,形成在最大功率点附近振荡,很难到达最大功率点。
发明内容
为了解决上述技术问题,本发明的目的是提供一种能提高低风速时最大功率点的搜索速度,且在高风速接近最大功率点时避免振荡的一种基于改进爬山算法的风力发电最大功率点跟踪方法。
本发明所采用的技术方案是:
一种基于改进爬山算法的风力发电最大功率点跟踪方法,包括以下步骤:
A、测量风机当前的功率值P;
B、对风机的转速n增加一个搜索步长Δn(i),使得风机转速变为n+Δn(i);
C、测量风机转速变化后的功率变化值ΔP;
D、根据功率值P,计算功率值影响系数K1;
E、根据搜索步长Δn(i)和功率变化值ΔP,计算得出斜率影响系数K2;
F、根据功率值影响系数K1、斜率影响系数K2、本次的搜索步长Δn(i)和功率变化值ΔP,计算下一次的搜索步长Δn(i+1);
G、根据预设的条件判断是否结束搜索,若是,则结束;反之,则根据计算得到的下一次的搜索步长Δn(i+1)返回执行步骤A。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述功率值影响系数K1的计算公式为:
其中,K1_max和K1_min分别为最大和最小功率值影响系数,P为功率值,A、B分别为K1-P曲线的斜率和截距,P1为对应最大功率值影响系数的转折功率,P2为对应最小功率值影响系数的转折功率。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述斜率影响系数K2的计算公式为:
K2=C|△P/△n(i)|;
其中,Δn(i)表示本次的搜索步长,ΔP表示功率变化值,C为预设的常数项。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述下一次的搜索步长Δn(i+1)的计算公式为:
△n(i+1)=sign(△P)*sign(△n(i))*K1*K2;
其中,K1表示功率值影响系数,K2表示斜率影响系数,Δn(i)表示本次的搜索步长,ΔP表示功率变化值。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述下一次的搜索步长Δn(i+1)的大小与斜率的绝对值|△P/△n(i)|正相关,与功率值影响系数K1正相关,与功率值P负相关,当功率值P越小时,计算得出的下一次搜索步长Δn(i+1)会越大,向最大功率点接近的速度越快;当功率值P越接近最大功率点时,计算得出的下一次搜索步长Δn(i+1)会越小,功率值P向最大功率点逐步趋近。
本发明的有益效果是:
本发明一种基于改进爬山算法的风力发电最大功率点跟踪方法通过引入功率值影响系数和斜率影响系数,能根据当前功率值变化的实际情况改变下一次搜索步长,从而实现加快低风速时搜索最大功率点搜索的速度的目的,并且在高风速时,避免了风机在最大功率点的振荡情况。
附图说明
下面结合附图对本发明的具体实施方式作进一步说明:
图1是本发明一种基于改进爬山算法的风力发电最大功率点跟踪方法的步骤流程图。
具体实施方式
参考图1,本发明一种基于改进爬山算法的风力发电最大功率点跟踪方法,包括以下步骤:
A、测量风机当前的功率值P;
B、对风机的转速n增加一个搜索步长Δn(i),使得风机转速变为n+Δn(i);
C、测量风机转速变化后的功率变化值ΔP;
D、根据功率值P,计算功率值影响系数K1;
E、根据搜索步长Δn(i)和功率变化值ΔP,计算得出斜率影响系数K2;
F、根据功率值影响系数K1、斜率影响系数K2、本次的搜索步长Δn(i)和功率变化值ΔP,计算下一次的搜索步长Δn(i+1);
G、根据预设的条件判断是否结束搜索,若是,则结束;反之,则根据计算得到的下一次的搜索步长Δn(i+1)返回执行步骤A。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述功率值影响系数K1的计算公式为:
其中,K1_max和K1_min分别为最大和最小功率值影响系数,P为功率值,A、B分别为K1-P曲线的斜率和截距,P1为对应最大功率值影响系数的转折功率,P2为对应最小功率值影响系数的转折功率。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述斜率影响系数K2的计算公式为:
K2=C|△P/△n(i)|;
其中,Δn(i)表示本次的搜索步长,ΔP表示功率变化值,C为预设的常数项。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述下一次的搜索步长Δn(i+1)的计算公式为:
△n(i+1)=sign(△P)*sign(△n(i))*K1*K2;
其中,K1表示功率值影响系数,K2表示斜率影响系数,Δn(i)表示本次的搜索步长,ΔP表示功率变化值,sign为符号运算符。
作为所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法的进一步改进,所述下一次的搜索步长Δn(i+1)的大小与斜率的绝对值|△P/△n(i)|正相关,与功率值影响系数K1正相关,与功率值P负相关,当功率值P越小时,计算得出的下一次搜索步长Δn(i+1)会越大,向最大功率点接近的速度越快;当功率值P越接近最大功率点时,计算得出的下一次搜索步长Δn(i+1)会越小,功率值P向最大功率点逐步趋近。
本发明实施例中,当风力发电控制器初次进入步骤A时,Δn将被赋予一个初值Δn(0),而从第二次搜索开始,Δn将用上一次搜索过程的步骤F计算所得的结果进行迭代。
在功率值影响系数K1的计算公式中,功率值影响系数K1与功率值P负相关,即P越小,K1越大。P值小代表着两种情况,一是风速高,且风机工作在离功率-转速曲线的峰值点比较远的位置,二是风速较低时。根据斜率影响系数K2的计算公式和下一次搜索步长Δn(i+1)的计算公式,下一次的搜索步长的大小与斜率的绝对值|△P/△n(i)|正相关,与功率值影响系数K1正相关,即与功率值P负相关。因此,无论风机工作在远离功率-转速曲线峰值点的点时,即|△P/△n(i)|较大,P较小的情况,还是风速较低时,即P较小的情况,搜索步长都会变大,从而加快了这两种情况下向最大功率点接近的速度。反之,当风速较高,又接近最大功率点时,搜索步长就会随之变小,避免了高风速时风机在最大功率点附近振荡。
在下一次的搜索步长Δn(i+1)的计算公式中,sign(△P)*sign(△n(i))则决定了下一次搜索步长的正负,即决定了下一次搜索的方向。当风机工作在最大功率点的左侧时,sign(△P)*sign(△n(i))>0,下一次搜索步长为正,即下一次执行步骤B时,风机转速增加,从而进一步接近最大功率点。当风机工作在最大功率点的右侧时,sign(△P)*sign(△n(i))<0,下一次搜索步长为负,即下一次执行步骤B时,风机转速减少,从而接近最大功率点。
其中,本发明针对的仅是加快风机从功率-转速曲线的两侧向峰值位置搜索的速度,并不包括搜索停止机制。因此可以根据自定义设定的条件判断是否结束搜索。本发明具体实施例中,可以设定当斜率的绝对值小于一个阈值ε时(即斜率很小,接近于零,此时风机位于功率-转速曲线的峰值点附近,相当于|△P/△n(i)|接近零),搜索停止。本发明亦不包括搜索停止后的搜索重启机制。本发明具体实施例中,可以设定当功率值的变化大于一个阈值β时(功率发生变化,即表明风速发生了变化,需要重新搜索当前风速下的最大功率点,在本发明中,相当于ΔP大于一个阈值β),搜索将重启。无论使用何种搜索停止和重启机制,只要搜索最大功率点的过程使用了本方法,都应当属于本发明需要保护的范围。
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (5)
1.一种基于改进爬山算法的风力发电最大功率点跟踪方法,其特征在于,包括以下步骤:
A、测量风机当前的功率值P;
B、对风机的转速n增加一个搜索步长Δn(i),使得风机转速变为n+Δn(i);
C、测量风机转速变化后的功率变化值ΔP;
D、根据功率值P,计算功率值影响系数K1;
E、根据搜索步长Δn(i)和功率变化值ΔP,计算得出斜率影响系数K2;
F、根据功率值影响系数K1、斜率影响系数K2、本次的搜索步长Δn(i)和功率变化值ΔP,计算下一次的搜索步长Δn(i+1);
G、根据预设的条件判断是否结束搜索,若是,则结束;反之,则根据计算得到的下一次的搜索步长Δn(i+1)返回执行步骤A。
2.根据权利要求1所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法,其特征在于:所述功率值影响系数K1的计算公式为:
其中,K1_max和K1_min分别为最大和最小功率值影响系数,P为功率值,A、B分别为K1-P曲线的斜率和截距,P1为对应最大功率值影响系数的转折功率,P2为对应最小功率值影响系数的转折功率。
3.根据权利要求1所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法,其特征在于:所述斜率影响系数K2的计算公式为:
K2=C|ΔP/Δn(i)|;
其中,Δn(i)表示本次的搜索步长,ΔP表示功率变化值,C为预设的常数项。
4.根据权利要求1所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法,其特征在于:所述下一次的搜索步长Δn(i+1)的计算公式为:
Δn(i+1)=sign(ΔP)*sign(Δn(i))*K1*K2;
其中,K1表示功率值影响系数,K2表示斜率影响系数,Δn(i)表示本次的搜索步长,ΔP表示功率变化值。
5.根据权利要求1所述的一种基于改进爬山算法的风力发电最大功率点跟踪方法,其特征在于:所述下一次的搜索步长Δn(i+1)的大小与斜率的绝对值|ΔP/Δn(i)|正相关,与功率值影响系数K1正相关,与功率值P负相关,当功率值P越小时,计算得出的下一次搜索步长Δn(i+1)会越大,向最大功率点接近的速度越快;当功率值P越接近最大功率点时,计算得出的下一次搜索步长Δn(i+1)会越小,功率值P向最大功率点逐步趋近。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510209324.0A CN104929863B (zh) | 2015-04-28 | 2015-04-28 | 一种基于改进爬山算法的风力发电最大功率点跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510209324.0A CN104929863B (zh) | 2015-04-28 | 2015-04-28 | 一种基于改进爬山算法的风力发电最大功率点跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104929863A true CN104929863A (zh) | 2015-09-23 |
CN104929863B CN104929863B (zh) | 2017-10-31 |
Family
ID=54117219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510209324.0A Active CN104929863B (zh) | 2015-04-28 | 2015-04-28 | 一种基于改进爬山算法的风力发电最大功率点跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104929863B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108390406A (zh) * | 2018-02-05 | 2018-08-10 | 易事特集团股份有限公司 | 基于无刷双馈电机的风力发电系统及其控制方法 |
CN109944822A (zh) * | 2018-11-11 | 2019-06-28 | 南京理工大学 | 基于限制参考转速变化率的转速跟踪目标优化方法 |
CN110397554A (zh) * | 2019-09-05 | 2019-11-01 | 国电联合动力技术有限公司 | 智能寻优的风电机组偏航控制方法、装置及风电机组 |
CN110645145A (zh) * | 2018-06-27 | 2020-01-03 | 新疆金风科技股份有限公司 | 风力发电机组的控制方法和控制设备 |
CN112302865A (zh) * | 2019-07-31 | 2021-02-02 | 北京金风科创风电设备有限公司 | 风力发电机组的最优增益追踪方法和设备 |
CN114548611A (zh) * | 2022-04-27 | 2022-05-27 | 东方电气风电股份有限公司 | 一种风力发电机组最优增益参数的搜索方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101603502A (zh) * | 2008-06-11 | 2009-12-16 | 武汉事达电气股份有限公司 | 一种基于人工智能的风能控制方法 |
CN102242689A (zh) * | 2011-06-24 | 2011-11-16 | 南京理工大学 | 基于风力发电的最大功率点跟踪控制的改进爬山算法 |
CN102352814A (zh) * | 2011-06-21 | 2012-02-15 | 北京交通大学 | 一种大型直驱风电机组的最大功率跟踪方法 |
CN202228271U (zh) * | 2011-09-13 | 2012-05-23 | 辽宁力迅风电控制系统有限公司 | 离网型风力发电系统的最大功率跟踪装置 |
-
2015
- 2015-04-28 CN CN201510209324.0A patent/CN104929863B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101603502A (zh) * | 2008-06-11 | 2009-12-16 | 武汉事达电气股份有限公司 | 一种基于人工智能的风能控制方法 |
CN102352814A (zh) * | 2011-06-21 | 2012-02-15 | 北京交通大学 | 一种大型直驱风电机组的最大功率跟踪方法 |
CN102242689A (zh) * | 2011-06-24 | 2011-11-16 | 南京理工大学 | 基于风力发电的最大功率点跟踪控制的改进爬山算法 |
CN202228271U (zh) * | 2011-09-13 | 2012-05-23 | 辽宁力迅风电控制系统有限公司 | 离网型风力发电系统的最大功率跟踪装置 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108390406A (zh) * | 2018-02-05 | 2018-08-10 | 易事特集团股份有限公司 | 基于无刷双馈电机的风力发电系统及其控制方法 |
CN110645145A (zh) * | 2018-06-27 | 2020-01-03 | 新疆金风科技股份有限公司 | 风力发电机组的控制方法和控制设备 |
CN109944822A (zh) * | 2018-11-11 | 2019-06-28 | 南京理工大学 | 基于限制参考转速变化率的转速跟踪目标优化方法 |
CN109944822B (zh) * | 2018-11-11 | 2020-07-07 | 南京理工大学 | 基于限制参考转速变化率的转速跟踪目标优化方法 |
CN112302865A (zh) * | 2019-07-31 | 2021-02-02 | 北京金风科创风电设备有限公司 | 风力发电机组的最优增益追踪方法和设备 |
CN110397554A (zh) * | 2019-09-05 | 2019-11-01 | 国电联合动力技术有限公司 | 智能寻优的风电机组偏航控制方法、装置及风电机组 |
CN114548611A (zh) * | 2022-04-27 | 2022-05-27 | 东方电气风电股份有限公司 | 一种风力发电机组最优增益参数的搜索方法 |
CN114548611B (zh) * | 2022-04-27 | 2022-07-19 | 东方电气风电股份有限公司 | 一种风力发电机组最优增益参数的搜索方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104929863B (zh) | 2017-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104929863A (zh) | 一种基于改进爬山算法的风力发电最大功率点跟踪方法 | |
CA2755024C (en) | Hybrid intelligent control method and system for power generating apparatuses | |
Mesemanolis et al. | High-efficiency control for a wind energy conversion system with induction generator | |
Gloe et al. | Continuous provision of synthetic inertia with wind turbines: implications for the wind turbine and for the grid | |
CN103441529A (zh) | 一种变速风电机组模拟惯量响应控制方法 | |
CN102072083B (zh) | 一种双馈风力发电系统最大风能跟踪方法 | |
Thongam et al. | Neural network based wind speed sensorless MPPT controller for variable speed wind energy conversion systems | |
CN110454328A (zh) | 一种无风速传感器下的风力发电系统功率追踪方法 | |
CN104948385A (zh) | 风电机组恒额定转速区域提高发电量的变桨控制方法 | |
CN107394817B (zh) | 一种风电参与电力系统调频的方法及系统 | |
CN105649875A (zh) | 风力发电机组的变桨控制方法及装置 | |
CN105673324A (zh) | 一种实现风电机组最大功率点跟踪的方法 | |
Thongam et al. | An optimum speed MPPT controller for variable speed PMSG wind energy conversion systems | |
CN105332853A (zh) | 在额定风速以上协调风力发电机转速与功率的控制方法 | |
Liao et al. | Hybrid control of DFIGs for short‐term and long‐term frequency regulation support in power systems | |
Gajewski et al. | Analysis of Sliding Mode Control of variable speed wind turbine system with PMSG | |
Raza et al. | A novel algorithm for fast and efficient maximum power point tracking of wind energy conversion systems | |
CN103291541B (zh) | 风力发电机组控制方法及控制系统 | |
Merabet et al. | Adaptive sliding mode speed control for wind turbine systems | |
CN105484940B (zh) | 最大功率点跟踪控制方法及风力发电系统 | |
Abdelhameed et al. | Adaptive maximum power tracking control technique for wind energy conversion systems | |
CN107524563A (zh) | 一种基于滑模极值搜索的控制方法 | |
CN104500338B (zh) | 一种风力发电主动偏航变转速失速控制系统 | |
Raza et al. | A novel speed-sensorless adaptive hill climbing algorithm for fast and efficient maximum power point tracking of wind energy conversion systems | |
Jia et al. | Fluctuating characteristic and power smoothing strategies of WECS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |