CN110454328A - 一种无风速传感器下的风力发电系统功率追踪方法 - Google Patents
一种无风速传感器下的风力发电系统功率追踪方法 Download PDFInfo
- Publication number
- CN110454328A CN110454328A CN201910740476.1A CN201910740476A CN110454328A CN 110454328 A CN110454328 A CN 110454328A CN 201910740476 A CN201910740476 A CN 201910740476A CN 110454328 A CN110454328 A CN 110454328A
- Authority
- CN
- China
- Prior art keywords
- wind
- power
- generator system
- speed
- wind speed
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 238000005457 optimization Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 6
- 238000004804 winding Methods 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 5
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 230000009885 systemic effect Effects 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims 1
- 230000002411 adverse Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000013440 design planning Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P9/00—Arrangements for controlling electric generators for the purpose of obtaining a desired output
- H02P9/04—Control effected upon non-electric prime mover and dependent upon electric output value of the generator
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/10—Purpose of the control system
- F05B2270/101—Purpose of the control system to control rotational speed (n)
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/10—Purpose of the control system
- F05B2270/103—Purpose of the control system to affect the output of the engine
- F05B2270/1033—Power (if explicitly mentioned)
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/70—Type of control algorithm
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/70—Type of control algorithm
- F05B2270/707—Type of control algorithm fuzzy logic
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2101/00—Special adaptation of control arrangements for generators
- H02P2101/15—Special adaptation of control arrangements for generators for wind-driven turbines
-
- 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
Abstract
一种无风速传感器下的风力发电系统功率追踪方法,包括如下步骤:建立风力发电系统最大功率追踪控制的运动方程;采用变步长算法,对变风速下的输出功率和风轮转速进行迭代优化,进行最大功率点的跟踪,并计算得到最优叶尖速比;采用模糊控制器作为功率控制器,以叶尖速比的误差和误差变化率为输入,以定子侧调整电压的变化量作为输出,制定模糊规则,并将模糊控制器的输出进行解模糊运算,得到相应的准确变量;根据定子侧调整电压的变化量调整发电机转速,实现控制跟踪风力发电系统的风轮转速以及对风力发电系统的最大功率进行追踪;本发明避免了在风速测量不准确以及风速测量滞后的情况下对风电系统输出功率的不良影响。
Description
技术领域
本发明涉及风力发电领域,具体涉及一种无风速传感器下的风力发电系统功率追踪方法。
背景技术
煤炭、石油、天然气等是应用最为广泛的能源,它们都是不可再生的化石资源,而风能是一种经济、环境友好型、且储量巨大的可再生资源,因此世界各国都在加紧对风能的开发和利用。风力发电是目前世界上利用风能的主要形式,由于风力发电具有良好的前景,开发利用风力资源对于缓解能源短缺、保护生态环境具有重要意义。
在风里发电系统中,风速、风向以及风的强度都具有高度的随机性、瞬变性,因此,为了保证风力发电系统的持续最大功率发电,发明一种能够自动调整风轮机转速进而捕获或追踪风力发电的最大功率点,对风力发电系统具有重大意义。
目前变速风电系统一般采用最大功率点追踪控制策略,且主要有以下几种追踪方法:最优叶尖速比法、功率信号反馈法、直线扰动法以及爬山搜索法等,其中最优叶尖速比法是目前主要流行的方法之一,该方法需要建立在风速测量的基础上,再根据最优叶尖速比得到对应的最优转速,通过控制器改变发电机转速实现最大风能追踪,然而,由于风速的测量在实际应用中存在较多的问题,例如,风速计的安装位置、塔架、湍流和温度等因素会对风速的测量结果产生影响,而且风速的测量本身具有一定的滞后性,这些客观存在的问题都将影响风力发电系统功率追踪控制的性能表现;至于另外的几种主流方法,都需要对风力发电系统进行有针对性的设置,增加了系统设计的复杂度,实现较为繁琐。
发明内容
有鉴于此,为了解决现有技术中的风力发电系统在功率优化和控制过程中存有的技术问题,本发明经过理论研究与实验验证,提供一种无风速传感器下的风力发电系统功率追踪方法,应用在无风速传感器的风力发电系统中,在风力发电系统的功率特性未知的情况下,对风力发电系统进行最大功率点的追踪,不仅可以较为简易地得到实现和实际应用,具有较强的稳定性和鲁棒性,而且也节约了系统设计成本,尤其避免了现有的在风速测量不准确以及风速测量滞后的情况下对风电系统输出功率的不良影响。
基于上述目的,本发明提供一种无风速传感器下的风力发电系统功率追踪方法,在风力发电系统的功率特性未知的情况下,进行如下步骤追踪所述风力发电系统在变风速下的最大功率点:
步骤1:根据所述风力发电系统发电机的选型,建立风力发电系统最大功率追踪控制的运动方程;
步骤2:采用变步长算法,对变风速下的风力发电系统的输出功率P和风轮转速ωr进行迭代优化,进行风力发电系统最大功率点的跟踪,并计算得到最优叶尖速比λopt,并将其作为叶尖速比的参考值带入模糊控制器中;
步骤3:采用模糊控制器作为风力发电系统的功率控制器,所述模糊控制器采用二维控制结构,以风力发电系统的叶尖速比的误差e和误差变化率ec为输入,以风力发电系统的定子侧调整电压u的变化量Δu作为输出,制定出相应的模糊规则,并将模糊控制器的输出进行解模糊运算,得到相应的准确变量;
步骤4:根据定子侧调整电压u的变化量Δu调整发电机转速,控制发电机转速向最大功率输出下的额定转速靠近,间接实现控制跟踪风力发电系统的风轮转速ωr,opt以及对风力发电系统的最大功率进行追踪。
优选地,在所述步骤1中,首先建立电机电磁转矩表达式:
上式中,参数p是发电机极对数,m是发电机的相数,u是发电机输出电压,也即定子侧调整电压,c为发电机修正系数,ωg是发电机转速,ω1是发电机的同步转速,r1、x1是发电机定子绕组的电阻和漏抗,r2、x2是发电机转子绕组的电阻和漏抗;
然后,建立如下发电机运动方程:
其中,参数Tm是风轮的气动转矩,Te是电机的电磁转矩,μe是电机侧阻尼系数,Jg是发电机转动惯量,具体地,风轮的气动转矩Tm可通过如下表示式计算:
其中,参数P表示风力发电系统的输出功率,ρ是空气密度,R是风轮半径,υ是风速,ωr是风轮转速,Cp(λ,β)是风能利用系数,是桨距角β和叶尖速比λ的函数,可以按如下公式计算风能利用系数Cp(λ,β):
在额定风速以下,桨距角β通常为零,即β=0;叶尖速比λ的计算表达式为:
λ=Rωr/υ
优选地,在所述步骤2中,根据如下公式确定最大功率点对应的最优功率系数kopt:
其中,参数Cp,max表示风能利用系数的最大值,λopt分别表示最大功率点对应的最优叶尖速比;
所述变步长算法的核心即是由系统前一时刻的输出功率P(i-1)和转速ωr(i-1)得到当前时刻最优功率系数k(i),并迭代计算获得风轮转速变化量Δωr(i),结合当前时刻的输出功率P(i),得到下一时刻的风轮转速ωr(i+1),最终迭代得到最大功率点,具体运算过程是:
步骤2.1:根据风力发电系统的初始输出功率P(0)以及风轮转速ωr(0)计算获得初始功率系数k(0),计算公式如下:
步骤2.2:按照以下公式进行输出功率P和转速ωr的迭代:
ΔP(i)=P(i)-P(i-1)
如果当前迭代时间步的输出功率变化量|ΔP(i)|<ε1,表明两个时间步的输出功率差值处于较小范围,则加大风轮转速的扰动值,令
Δωr(i)=1.618(ωr(i)-ωr(i-1))
如果当前迭代时间步的输出功率变化量ε1<|ΔP(i)|<ε2,表明两个时间步的输出功率差值处于较大范围,则减小风轮转速的扰动值,令
Δωr(i)=0.618(ωr(i)-ωr(i-1))
其中,参数i表示迭代时间步数,ε1和ε2分别表示系统预设的阈值,常系数1.618和0.618分别是不同情况下的利用黄金分割点对迭代步长的优化系数;
步骤2.3:根据以下两式判断迭代是否终止:
ΔP(i)<0
|Δωr(i)|<ε3
上式中,参数ε3为系数预设的阈值,用来剔除由于风速减小造成的最大功率点的误判;
若有任一不等式不成立,利用公式ωr(i+1)=ωr(i)+Δωr(i)计算得到下一迭代步的风轮转速ωr(i+1),并进行下一迭代,直至当前时刻的输出功率变化量ΔP(i)和风轮转速变化量绝对值|Δωr(i)|满足终止条件;最终获得最大功率点对应的风轮转速ωr,opt以及最优功率系数kopt,通过公式计算即可获得最优叶尖速比λopt。
优选地,在所述步骤3中,所述风力发电系统的叶尖速比的误差e的计算公式是:
e=λopt-λt
其中,λt表示风力发电系统的实际叶尖速比。
相比于现有技术,本发明的技术效果是:
(1)本发明基于改进的叶尖速比控制策略,在不需要测量当前实际风速以及风电系统的功率特性曲线的基础上,通过变步长算法的迭代,计算获得风电系统的最大功率点对应的风轮转速ωr,opt以及最优功率系数kopt,进而获得最优叶尖速比λopt,由于传统的最优叶尖速比控制策略是通过仿真计算得到,需要事先获得风电系统的功率特性曲线,最终的功率追踪精度得不到良好保证,而本发明可以避免因为风电系统功率特性曲线的测量误差带来的最大功率点追踪误差的影响。
(2)本发明在追踪最大功率点时,采用变步长的迭代算法,在风力发电系统输出功率离最大功率点较近时,加大汽轮转速变化量的扰动值,在输出功率离最大功率点较远时,减小汽轮转速变化量的扰动值,从而提高了风力发电系统最大功率点的追踪速度,并减少了扰动功率的浪费,节约资源,从而解决了传统的最优叶尖速比法、功率信号反馈法等方法在实际应用中的实现困难和成本高的问题。
附图说明
图1是本发明的风力发电系统功率追踪方法流程示意图。
具体实施方式
以下描述用于揭露本发明以使本领域技术人员能够实现本发明。以下描述中的实施例只作为举例,本领域技术人员可以想到其他显而易见的变型。在以下描述中界定的本发明的基本原理可以应用于其他实施方案、变形方案、改进方案、等同方案以及没有背离本发明的精神和范围的其他技术方案。
参照附图1,本发明的原理如下:一种无风速传感器下的风力发电系统功率追踪方法,在风力发电系统的功率特性未知的情况下,进行如下步骤追踪所述风力发电系统在变风速下的最大功率点:
步骤1:根据所述风力发电系统发电机的选型,建立风力发电系统最大功率追踪控制的运动方程;
在优选的实施例中,可以采用双馈异步感应电机或者永磁同步电机作为风力发电系统的发电机,在建立风力发电系统的控制方程时,首先需要先建立电机电磁转矩表达式:
上式中,参数p是发电机极对数,m是发电机的相数,u是发电机输出电压,也即定子侧调整电压,c为发电机修正系数,ωg是发电机转速,ω1是发电机的同步转速,r1、x1是发电机定子绕组的电阻和漏抗,r2、x2是发电机转子绕组的电阻和漏抗;
然后,建立如下发电机运动方程:
其中,参数Tm是风轮的气动转矩,Te是电机的电磁转矩,μe是电机侧阻尼系数,Jg是发电机转动惯量,具体地,风轮的气动转矩Tm可通过如下表示式计算:
其中,参数P表示风力发电系统的输出功率,ρ是空气密度,R是风轮半径,υ是风速,ωr是风轮转速,Cp(λ,β)是风能利用系数,是桨距角β和叶尖速比λ的函数,可以按如下公式计算风能利用系数Cp(λ,β):
在额定风速以下,桨距角β通常为零,即β=0;叶尖速比λ的计算表达式为:
λ=Rωr/υ
需要指出的是,在风力发电系统运动过程中,诸如空气密度ρ、风轮半径R以及桨距角β等参数是固定不变的,即不受风力发电系统控制的影响,且默认为发电机的输出电功率即为风力发电系统的输出功率。
步骤2:采用变步长算法,对变风速下的风力发电系统的输出功率P和风轮转速ωr进行迭代优化,进行风力发电系统最大功率点的跟踪,并计算得到最优叶尖速比λopt,并将其作为叶尖速比的参考值带入模糊控制器中;
在优选的实施例中,根据如下公式确定最大功率点对应的最优功率系数kopt:
其中,参数Cp,max表示风能利用系数的最大值,λopt分别表示最大功率点对应的最优叶尖速比;
所述变步长算法的核心即是由系统前一时刻的输出功率P(i-1)和转速ωr(i-1)得到当前时刻最优功率系数k(i),并迭代计算获得风轮转速变化量Δωr(i),结合当前时刻的输出功率P(i),得到下一时刻的风轮转速ωr(i+1),最终迭代得到最大功率点,具体运算过程是:
步骤2.1:根据风力发电系统的初始输出功率P(0)以及风轮转速ωr(0)计算获得初始功率系数k(0),计算公式如下:
步骤2.2:按照以下公式进行输出功率P和转速ωr的迭代:
ΔP(i)=P(i)-P(i-1)
如果当前迭代时间步的输出功率变化量|ΔP(i)|<ε1,表明两个时间步的输出功率差值处于较小范围,则加大风轮转速的扰动值,令
Δωr(i)=1.618(ωr(i)-ωr(i-1))
如果当前迭代时间步的输出功率变化量ε1<|ΔP(i)|<ε2,表明两个时间步的输出功率差值处于较大范围,则减小风轮转速的扰动值,令
Δωr(i)=0.618(ωr(i)-ωr(i-1))
其中,参数i表示迭代时间步数,ε1和ε2分别表示系统预设的阈值,常系数1.618和0.618分别是不同情况下的利用黄金分割点对迭代步长的优化系数;
步骤2.3:根据以下两式判断迭代是否终止:
ΔP(i)<0
|Δωr(i)|<ε3
上式中,参数ε3为系数预设的阈值,用来剔除由于风速减小造成的最大功率点的误判;
若有任一不等式不成立,利用公式ωr(i+1)=ωr(i)+Δωr(i)计算得到下一迭代步的风轮转速ωr(i+1),并进行下一迭代,直至当前时刻的输出功率变化量ΔP(i)和风轮转速变化量绝对值|Δωr(i)|满足终止条件;最终获得最大功率点对应的风轮转速ωr,opt以及最优功率系数kopt,通过公式计算即可获得最优叶尖速比λopt,然后实施步骤3。
步骤3:采用模糊控制器作为风力发电系统的功率控制器,所述模糊控制器采用二维控制结构,以风力发电系统的叶尖速比的误差e和误差变化率ec为输入,以风力发电系统的定子侧调整电压u的变化量Δu作为输出,制定出相应的模糊规则,并将模糊控制器的输出进行解模糊运算,得到相应的准确变量;
在优选的实施例中,所述风力发电系统的叶尖速比的误差e的计算公式是:
e=λopt-λt
其中,λt表示风力发电系统的实际叶尖速比。
步骤4:根据定子侧调整电压u的变化量Δu调整发电机转速,控制发电机转速向最大功率输出下的额定转速靠近,间接实现控制跟踪风力发电系统的风轮转速ωr,opt以及对风力发电系统的最大功率进行追踪。
需要指出的是,可以选用现有的成熟的控制策略进行模糊控制,如采用模糊PID控制器,其具有常规PID控制器的基本特征,并采用模糊逻辑推理方法来调整PID控制算法的参数,经过模糊推理得到的结果不是直接作为模糊控制系统的输出,而是用该结果来整定PID的参数,再根据PID算法来决定模糊控制系统的最终输出;同时,基于模糊控制原理,控制器的维数越多,控制精度越高,但是控制器相应的设计规划难度也会以指数形式增长,实际使用起来更困难,从实际使用角度出发,本申请采用二维的模糊控制器,能够满足最终功率追踪精度的要求;本发明在具体设计模糊控制规则时,可参考现有的成熟控制算法,由此具体过程属于现有技术范畴,本发明在此不一一赘述。
本领域的技术人员应理解,上述描述及附图中的本发明的实施例只作为举例而不限制本发明。本发明的目的已经完整并有效地实现。本发明的功能及结构原理已在实施例中展示和说明,在没有背离原理下,本发明的实施方法可以有任何形式或修改。
Claims (4)
1.一种无风速传感器下的风力发电系统功率追踪方法,其特征在于,在风力发电系统的功率特性未知的情况下,进行如下步骤追踪所述风力发电系统在变风速下的最大功率点:
步骤1:根据所述风力发电系统发电机的选型,建立风力发电系统最大功率追踪控制的运动方程;
步骤2:采用变步长算法,对变风速下的风力发电系统的输出功率P和风轮转速ωr进行迭代优化,进行风力发电系统最大功率点的跟踪,并计算得到最优叶尖速比λopt,并将其作为叶尖速比的参考值带入模糊控制器中;
步骤3:采用模糊控制器作为风力发电系统的功率控制器,所述模糊控制器采用二维控制结构,以风力发电系统的叶尖速比的误差e和误差变化率ec为输入,以风力发电系统的定子侧调整电压u的变化量Δu作为输出,制定出相应的模糊规则,并将模糊控制器的输出进行解模糊运算,得到相应的准确变量;
步骤4:根据定子侧调整电压u的变化量Δu调整发电机转速,控制发电机转速向最大功率输出下的额定转速靠近,间接实现控制跟踪风力发电系统的风轮转速ωr,opt以及对风力发电系统的最大功率进行追踪。
2.如权利要求1所述的风力发电系统功率追踪方法,在所述步骤1中,首先建立电机电磁转矩表达式:
上式中,参数p是发电机极对数,m是发电机的相数,u是发电机输出电压,也即定子侧调整电压,c为发电机修正系数,ωg是发电机转速,ω1是发电机的同步转速,r1、x1是发电机定子绕组的电阻和漏抗,r2、x2是发电机转子绕组的电阻和漏抗;
然后,建立如下发电机运动方程:
其中,参数Tm是风轮的气动转矩,Te是电机的电磁转矩,μe是电机侧阻尼系数,Jg是发电机转动惯量,具体地,风轮的气动转矩Tm可通过如下表示式计算:
其中,参数P表示风力发电系统的输出功率,ρ是空气密度,R是风轮半径,υ是风速,ωr是风轮转速,Cp(λ,β)是风能利用系数,是桨距角β和叶尖速比λ的函数,可以按如下公式计算风能利用系数Cp(λ,β):
在额定风速以下,桨距角β通常为零,即β=0;叶尖速比λ的计算表达式为:
λ=Rωr/υ
3.如权利要求2所述的风力发电系统功率追踪方法,在所述步骤2中,根据如下公式确定最大功率点对应的最优功率系数kopt:
其中,参数Cp,max表示风能利用系数的最大值,λopt分别表示最大功率点对应的最优叶尖速比;
所述变步长算法的核心即是由系统前一时刻的输出功率P(i-1)和转速ωr(i-1)得到当前时刻最优功率系数k(i),并迭代计算获得风轮转速变化量Δωr(i),结合当前时刻的输出功率P(i),得到下一时刻的风轮转速ωr(i+1),最终迭代得到最大功率点,具体运算过程是:
步骤2.1:根据风力发电系统的初始输出功率P(0)以及风轮转速ωr(0)计算获得初始功率系数k(0),计算公式如下:
步骤2.2:按照以下公式进行输出功率P和转速ωr的迭代:
ΔP(i)=P(i)-P(i-1)
如果当前迭代时间步的输出功率变化量|ΔP(i)|<ε1,表明两个时间步的输出功率差值处于较小范围,则加大风轮转速的扰动值,令
Δωr(i)=1.618(ωr(i)-ωr(i-1))
如果当前迭代时间步的输出功率变化量ε1<|ΔP(i)|<ε2,表明两个时间步的输出功率差值处于较大范围,则减小风轮转速的扰动值,令
Δωr(i)=0.618(ωr(i)-ωr(i-1))
其中,参数i表示迭代时间步数,ε1和ε2分别表示系统预设的阈值,常系数1.618和0.618分别是不同情况下的利用黄金分割点对迭代步长的优化系数;
步骤2.3:根据以下两式判断迭代是否终止:
ΔP(i)<0
|Δωr(i)|<ε3
上式中,参数ε3为系数预设的阈值,用来剔除由于风速减小造成的最大功率点的误判;
若有任一不等式不成立,利用公式ωr(i+1)=ωr(i)+Δωr(i)计算得到下一迭代步的风轮转速ωr(i+1),并进行下一迭代,直至当前时刻的输出功率变化量ΔP(i)和风轮转速变化量绝对值|Δωr(i)|满足终止条件;最终获得最大功率点对应的风轮转速ωr,opt以及最优功率系数kopt,通过公式计算即可获得最优叶尖速比λopt。
4.如权利要求3所述的风力发电系统功率追踪方法,在所述步骤3中,所述风力发电系统的叶尖速比的误差e的计算公式是:
e=λopt-λt
其中,λt表示风力发电系统的实际叶尖速比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910740476.1A CN110454328A (zh) | 2019-08-12 | 2019-08-12 | 一种无风速传感器下的风力发电系统功率追踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910740476.1A CN110454328A (zh) | 2019-08-12 | 2019-08-12 | 一种无风速传感器下的风力发电系统功率追踪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110454328A true CN110454328A (zh) | 2019-11-15 |
Family
ID=68485934
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910740476.1A Pending CN110454328A (zh) | 2019-08-12 | 2019-08-12 | 一种无风速传感器下的风力发电系统功率追踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110454328A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111222718A (zh) * | 2020-02-27 | 2020-06-02 | 上海电力大学 | 一种风能转换系统的最大功率点跟踪方法与装置 |
CN111336062A (zh) * | 2020-03-05 | 2020-06-26 | 中国大唐集团科学技术研究院有限公司华中电力试验研究院 | 一种基于测量风速的风力发电机组最大风能捕获方法 |
CN111541238A (zh) * | 2020-04-10 | 2020-08-14 | 东方电气风电有限公司 | 一种风力发电机组调频验证模型建模方法 |
CN111894800A (zh) * | 2020-07-16 | 2020-11-06 | 湘电风能有限公司 | 一种提升环境适应性的风力发电机转速控制方法 |
CN112486019A (zh) * | 2020-11-17 | 2021-03-12 | 九江学院 | 不确定风力发电机系统的最大功率跟踪模糊控制方法 |
CN113090453A (zh) * | 2019-12-23 | 2021-07-09 | 新疆金风科技股份有限公司 | 风力发电机组的控制方法、装置和风力发电机组 |
CN113623126A (zh) * | 2021-06-23 | 2021-11-09 | 湖南大学 | 一种基于模糊控制的直驱永磁水力发电系统控制方法、系统、终端及可读存储介质 |
CN114673629A (zh) * | 2022-04-28 | 2022-06-28 | 中国华能集团清洁能源技术研究院有限公司 | 一种串列式双风轮风电机组的最大功率跟踪方法及系统 |
CN115268559A (zh) * | 2022-05-23 | 2022-11-01 | 北京华能新锐控制技术有限公司 | 用于永磁同步风力发电机的最大功率点跟踪鲁棒控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776043A (zh) * | 2010-02-25 | 2010-07-14 | 东南大学 | 基于误差补偿模型的风电机组最大风能捕获控制方法 |
CN102072083A (zh) * | 2010-12-27 | 2011-05-25 | 西安理工大学 | 一种双馈风力发电系统最大风能跟踪方法 |
CN108050014A (zh) * | 2017-12-08 | 2018-05-18 | 河海大学 | 一种无风速仪的水平轴风力机变速控制方法 |
-
2019
- 2019-08-12 CN CN201910740476.1A patent/CN110454328A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776043A (zh) * | 2010-02-25 | 2010-07-14 | 东南大学 | 基于误差补偿模型的风电机组最大风能捕获控制方法 |
CN102072083A (zh) * | 2010-12-27 | 2011-05-25 | 西安理工大学 | 一种双馈风力发电系统最大风能跟踪方法 |
CN108050014A (zh) * | 2017-12-08 | 2018-05-18 | 河海大学 | 一种无风速仪的水平轴风力机变速控制方法 |
Non-Patent Citations (1)
Title |
---|
李南: "大型风电机组功率优化控制研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113090453A (zh) * | 2019-12-23 | 2021-07-09 | 新疆金风科技股份有限公司 | 风力发电机组的控制方法、装置和风力发电机组 |
CN113090453B (zh) * | 2019-12-23 | 2023-03-03 | 新疆金风科技股份有限公司 | 风力发电机组的控制方法、装置和风力发电机组 |
CN111222718A (zh) * | 2020-02-27 | 2020-06-02 | 上海电力大学 | 一种风能转换系统的最大功率点跟踪方法与装置 |
CN111336062A (zh) * | 2020-03-05 | 2020-06-26 | 中国大唐集团科学技术研究院有限公司华中电力试验研究院 | 一种基于测量风速的风力发电机组最大风能捕获方法 |
CN111541238A (zh) * | 2020-04-10 | 2020-08-14 | 东方电气风电有限公司 | 一种风力发电机组调频验证模型建模方法 |
CN111541238B (zh) * | 2020-04-10 | 2022-05-24 | 东方电气风电股份有限公司 | 一种风力发电机组调频验证模型建模方法 |
CN111894800A (zh) * | 2020-07-16 | 2020-11-06 | 湘电风能有限公司 | 一种提升环境适应性的风力发电机转速控制方法 |
CN112486019A (zh) * | 2020-11-17 | 2021-03-12 | 九江学院 | 不确定风力发电机系统的最大功率跟踪模糊控制方法 |
CN112486019B (zh) * | 2020-11-17 | 2022-06-07 | 九江学院 | 不确定风力发电机系统的最大功率跟踪模糊控制方法 |
CN113623126A (zh) * | 2021-06-23 | 2021-11-09 | 湖南大学 | 一种基于模糊控制的直驱永磁水力发电系统控制方法、系统、终端及可读存储介质 |
CN114673629A (zh) * | 2022-04-28 | 2022-06-28 | 中国华能集团清洁能源技术研究院有限公司 | 一种串列式双风轮风电机组的最大功率跟踪方法及系统 |
CN115268559A (zh) * | 2022-05-23 | 2022-11-01 | 北京华能新锐控制技术有限公司 | 用于永磁同步风力发电机的最大功率点跟踪鲁棒控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110454328A (zh) | 一种无风速传感器下的风力发电系统功率追踪方法 | |
CN102072083B (zh) | 一种双馈风力发电系统最大风能跟踪方法 | |
CN104079228B (zh) | 双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 | |
CN105226723A (zh) | 一种双馈风电机组基于风功率跟踪自适应调节的虚拟惯性控制方法 | |
CN105545595A (zh) | 基于径向基神经网络的风力机反馈线性化功率控制方法 | |
CN104410107A (zh) | 一种双馈风电系统的无源积分滑模控制方法 | |
CN105591395A (zh) | 一种双馈风机虚拟惯量控制后的转速恢复方法 | |
CN104675629A (zh) | 一种变速风力发电机组的最大风能捕获方法 | |
CN104481803A (zh) | 一种风力发电系统追踪最大输出功率控制方法 | |
CN112542855A (zh) | 一种双馈风力发电系统相量模型建模及仿真方法 | |
Du et al. | MPPT control strategy of DFIG-based wind turbines using double steps hill climb searching algorithm | |
CN108494300A (zh) | 一种基于永磁同步发电机扰动观测器的最大功率跟踪滑模控制方法 | |
Hossain | Application of wind energy into the transportation sector | |
CN102661243B (zh) | 一种双馈感应风电机组的预报校正变桨控制方法 | |
CN111049178A (zh) | 一种直驱永磁风电机组经vsc-hvdc并网稳定控制分析方法 | |
Messaoud et al. | Modeling and optimization of wind turbine driving permanent magnet synchronous generator | |
CN110417047B (zh) | 基于复转矩系数分析双馈风机ssci阻尼特性的方法 | |
CN107453662A (zh) | 基于自适应反推控制的机械弹性储能用PMSG 闭环I/f 控制方法 | |
Yan et al. | Transient modelling of doubly‐fed induction generator based wind turbine on full operation condition and rapid starting period based on low voltage ride‐through testing | |
Moulay et al. | Application and Control of a Doubly Fed Induction Machine Integrated in Wind Energy System | |
CN112682258B (zh) | 基于Backstepping的大型风力机最大功率点跟踪控制方法 | |
Xiao et al. | VSCF wind turbine control strategy for maximum power generation | |
Le-peng et al. | Simulation for strategy of maximal wind energy capture of doubly fed induction generators | |
Pan et al. | MPPT of Wind Power System Based on Hybrid Perturbation Observation and Global Fast Integral Sliding Mode | |
Yunqi | A novel optimum power fuzzy control strategy for doubly-fed wind turbine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191115 |