CN109139363B - 一种提升多机型风力机性能的最大功率点跟踪控制方法 - Google Patents
一种提升多机型风力机性能的最大功率点跟踪控制方法 Download PDFInfo
- Publication number
- CN109139363B CN109139363B CN201710449967.1A CN201710449967A CN109139363B CN 109139363 B CN109139363 B CN 109139363B CN 201710449967 A CN201710449967 A CN 201710449967A CN 109139363 B CN109139363 B CN 109139363B
- Authority
- CN
- China
- Prior art keywords
- wind turbine
- electromagnetic torque
- formula
- maximum power
- power point
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000011217 control strategy Methods 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
-
- 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
- F05B2260/00—Function
- F05B2260/84—Modelling or simulation
-
- 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
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Wind Motors (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种提升多机型风力机性能的最大功率点跟踪控制方法,该方法应用减小电磁转矩增益控制方法来实现最大功率点跟踪控制,根据风力机参数中的转动惯量来设置电磁转矩增益系数,其所用公式为:Kd=1‑α×J,且根据电磁转矩增益系数Kd和转动惯量J的上限值、下限值来确定电磁转矩增益系数Kd的调整系数α。本发明的改进方法根据风力机的转动惯量动态设定电磁转矩增益系数,能够改善多机型风力机的风能捕获效率,具有良好的适应性;同时将电磁转矩增益系数限定在最佳电磁转矩增益系数附近,从而使风力机获得更高的风能捕获效率;且该控制方法仅依赖风力机转动惯量这一风力机参数,无需复杂的迭代计算,简便易行。
Description
技术领域
本发明属于风力发电技术领域,具体地说是一种提升多机型风力机性能的最大功率点跟踪控制方法。
背景技术
为了提高低于额定风速区间的风能捕获效率,变速恒频风力发电机组一般采用最大功率点跟踪(Maximum Power Point Tracking,MPPT)控制策略。为了改善风力机的动态性能,美国国家可再生能源实验室的Johnson K.E.等人提出了减小电磁转矩增益(Decreased Torque Gain,DTG)控制,该控制方法通过设置发电机电磁转矩增益系数,以减小发电机电磁转矩来提高风力机在跟踪渐强阵风时的加速性能。
电磁转矩增益系数的设置将决定风力机的风能捕获效率,然而该系数与风力机转动惯量关系密切。研究表明,对于某一固定的风力机转动惯量,存在最佳电磁转矩增益系数使其风能捕获效率最大,而且不同的转动惯量具有不同的最佳电磁转矩增益系数,换言之,电磁转矩增益系数应随风力机转动惯量的变化而变化,因而采用固定的电磁转矩增益系数是不合理的。此外,某一地区或某一风电场内,实际采用的风力机型号较多,因此采用固定的电磁转矩增益系数并不能使所有风力机均达到最优的控制效果。
综上所述,根据不同机型动态设置电磁转矩增益系数,对于提升多机型风力机风能捕获效率十分必要。目前对于该方面的研究尚未见报道。
发明内容
本发明针对固定电磁转矩增益系数的控制方法难以适应多机型风力机最大功率点跟踪控制需求的问题,提供一种提升多机型风力机性能的最大功率点跟踪控制方法;具体而言,当某一地区采用不同型号的风力机时,该方法均能很好的适应各种机型以提升风力机的风能捕获效率。
本发明所要解决的技术问题是通过以下技术方案实现的:
一种提升多机型风力机性能的最大功率点跟踪控制方法,其特征在于:应用减小电磁转矩增益控制方法来实现最大功率点跟踪控制,该方法所用的公式为:
式(1)和式(2)为风力机数学模型,式(3)和式(4)为最大功率点跟踪控制策略;式中:J为转动惯量,Tm(v,ω)为风力机的机械转矩,Te(ω)为发电机电磁转矩,v为风速,ω为风力机转速,为风力机角加速度,ρ为空气密度,R为风力机半径,Cp(λ)为风能利用系数,λ=ωR/v是叶尖速比,ωbgn为起始发电转速,Kd为电磁转矩增益系数,Topt(ω)为风力机最优转矩,λopt为最佳叶尖速比,为最大风能利用系数;
根据风力机参数中的转动惯量来设置电磁转矩增益系数,其所用公式为:
Kd=1-α×J (5)
式(5)为最大功率点跟踪控制策略,式(5)中α为Kd的调整系数,该调整系数α的取值范围由转动惯量J的上限值、下限值和电磁转矩增益系数Kd的上限值、下限值确定,按式(3)、式(4)和式(5)调整发电机电磁转矩,即可实现最大功率点跟踪。
应用该方法实现最大功率点跟踪的具体步骤如下:
S1、输入风力机参数:风力机半径、转动惯量、最大风能利用系数、最佳叶尖速比,输入风电场环境参数:风电场空气密度;设定电磁转矩增益系数Kd的调整系数α;
S3、测量发电机实际功率Pe,并根据公式Te=Pe×ω计算发电机实际电磁转矩;
步骤S1中对调整系数α的设置方法为:
S11、确定应用该方法的风力机的转动惯量的取值范围[Jmin,Jmax],其中Jmin为多机型风力机的最小转动惯量,Jmax为最大转动惯量;
S12、设定电磁转矩增益系数Kd的取值范围,即Kdmin≤Kd≤Kdmax;
S13、根据Kd=1-α×J,结合J的取值范围[Jmin,Jmax],可得α的取值范围为
步骤S12中的电磁转矩增益系数Kd的下限值范围为[0.75,0.85]、上限值范围为[0.9,1.0]。
步骤S4中的控制器为PI控制器。
本发明相比现有技术有如下优点:
1)本发明是基于减小电磁转矩增益控制的改进方法,能够通过设置电磁转矩增益系数很好的提升风力机的跟踪性能;由于转动惯量是制约风力机跟踪性能的关键因素,该改进控制方法根据风力机的转动惯量动态设定电磁转矩增益系数,能够改善多机型风力机的风能捕获效率,具有良好的适应性。
2)本发明的控制方法通过设置电磁转矩增益系数的取值范围,将电磁转矩增益系数限定在最佳电磁转矩增益系数附近,从而使风力机获得更高的风能捕获效率。
3)本发明的控制方法仅依赖风力机转动惯量这一风力机参数,无需复杂的迭代计算,简便易行。
附图说明
附图1为一种提升多机型风力机性能的最大功率点跟踪控制方法原理图;
附图2为J=5.602×105kgm2时改进算法与DTG控制的风能利用系数的对比;
附图3为J=1.1204×106kgm2时改进算法与DTG控制的风能利用系数的对比;
附图4为J=2.2408×106kgm2时改进算法与DTG控制的风能利用系数的对比。
具体实施方式
为了进一步描述本发明的技术特点和效果,以下结合附图和具体实施方式对本发明做进一步描述。
本发明的实施例首先根据风力机参数确定电磁转矩增益系数Kd的调整系数α。然后通过仿真算例验证本发明的优越性。
一、实施例的仿真模型
(1)简化风力机模型的参数
在matlab/simulink中建立风力机仿真模型。三种风力机的转动惯量如表1所示。
表1三种风力机的转动惯量
仿真采用的其他主要参数如表2所示。
表2仿真采用的主要参数
(2)风速模型
本文利用matlab建立中长期风速模型,其中表征长期风速特征的平均风速根据Van der Hoven谱随机产生,短期湍流风速则采用卡尔曼滤波随机产生具有von Karman功率谱的短期风速时间序列。本文设置风速湍流等级为IEC-614000-1标准规定的A级,并随机产生50组风速时间序列,供仿真分析使用。
二、本发明方法的实现
(1)具体的实现方法
一种提升多机型风力机性能的最大功率点跟踪控制方法,该控制方法应用减小电磁转矩增益控制方法来实现最大功率点跟踪控制,该方法所用的公式为:
式(1)和式(2)为风力机数学模型,式(3)和式(4)为最大功率点跟踪控制策略;式中:J为转动惯量,Tm(v,ω)为风力机的机械转矩,Te(ω)为发电机电磁转矩,v为风速,ω为风力机转速,为风力机角加速度,ρ为空气密度,R为风力机半径,Cp(λ)为风能利用系数,λ=ωR/v是叶尖速比,ωbgn为起始发电转速,Kd为电磁转矩增益系数,Topt(ω)为风力机最优转矩,λopt为最佳叶尖速比,为最大风能利用系数;
根据风力机参数中的转动惯量来设置电磁转矩增益系数,其所用公式为:
Kd=1-α×J (5)
式(5)为最大功率点跟踪控制策略,式(5)中α为Kd的调整系数,该调整系数α的取值范围由转动惯量J的上限值、下限值和电磁转矩增益系数Kd的上限值、下限值确定,按式(3)、式(4)和式(5)调整发电机电磁转矩,即可实现最大功率点跟踪。
结合附图1,应用该方法实现最大功率点跟踪的具体步骤如下:
S1、输入风力机参数:风力机半径、转动惯量、最大风能利用系数、最佳叶尖速比,输入风电场环境参数:风电场空气密度;设定电磁转矩增益系数Kd的调整系数α;
S2、测量风力机转速ω,并根据式(3)、式(4)和式(5)确定的最大功率点跟踪控制策略计算发电机电磁转矩参考值
S3、测量发电机实际功率Pe,并根据公式Te=Pe×ω计算发电机实际电磁转矩;
上述步骤S1中对调整系数α的设置方法为:
S11、确定应用该方法的风力机的转动惯量的取值范围[Jmin,Jmax],其中Jmin为多机型风力机的最小转动惯量,Jmax为最大转动惯量;本实施例中Jmin=5.602×105kgm2,Jmax=2.2408×106kgm2;
S12、设定电磁转矩增益系数Kd的取值范围,即Kdmin≤Kd≤Kdmax;本实施例中Kdmin=0.8,Kdmax=1.0;
三、实施例的结果分析
采用50组模拟风速序列对本发明提出方法的优越性进行分析。具体地,将本发明提出的一种提升多机型风力机性能的最大功率点跟踪控制方法与DTG控制进行比较,以验证本发明提供的控制方法的优越性。
针对50组时长为10h的随机风速序列,分别应用DTG控制和本发明提出的方法,可获得每次迭代对应的平均风能利用率ηfavg,其表达式如式(9)所示;那么,整个风速序列的ηfavg的平均值记为如式(10)所示。
式(9)中,Pcap为实际功率,Pwy为最优功率,ψ为偏航误差角,本文设为0度,nc为一个迭代周期内的采样次数。
结合附图2、附图3和附图4,三种风力机应用DTG控制和本发明的方法的风能利用系数Cp的对比情况如附图2、附图3和附图4所示。附图2、附图3和附图4选取了某一条10h风速序列中600-700s的仿真结果。由附图2、附图3和附图4可知,本发明提出的方法具有比DTG控制更高的风能利用系数。
每种方法采用的电磁转矩增益系数Kd以及50组仿真算例对应的风能捕获效率如表3所示。由表3可知,DTG控制采用固定的电磁转矩增益系数(固定为0.8),而本发明根据风力机转动惯量动态调整电磁转矩增益系数,从而使得本发明提出的方法在风能捕获效率方面优于DTG控制。
本发明的控制方法是基于减小电磁转矩增益控制的改进方法,能够通过设置电磁转矩增益系数很好的提升风力机的跟踪性能;由于转动惯量是制约风力机跟踪性能的关键因素,该改进控制方法根据风力机的转动惯量动态设定电磁转矩增益系数,能够改善多机型风力机的风能捕获效率,具有良好的适应性;通过设置电磁转矩增益系数的取值范围,将电磁转矩增益系数限定在最佳电磁转矩增益系数附近,从而使风力机获得更高的风能捕获效率;该控制方法仅依赖风力机转动惯量这一风力机参数,无需复杂的迭代计算,简便易行。
上述实施例不以任何形式限定本发明,凡采取等同替换或等效变换的形式所获得的技术方案,均落在本发明的保护范围之内;本发明未涉及的技术均可通过现有技术加以实现。
Claims (5)
1.一种提升多机型风力机性能的最大功率点跟踪控制方法,其特征在于:应用减小电磁转矩增益控制方法来实现最大功率点跟踪控制,该方法所用的公式为:
式(1)和式(2)为风力机数学模型,式(3)和式(4)为最大功率点跟踪控制策略;式中:J为转动惯量,Tm(v,ω)为风力机的机械转矩,Te(ω)为发电机电磁转矩,v为风速,ω为风力机转速,为风力机角加速度,ρ为空气密度,R为风力机半径,Cp(λ)为风能利用系数,λ=ωR/v是叶尖速比,ωbgn为起始发电转速,Kd为电磁转矩增益系数,Topt(ω)为风力机最优转矩,λopt为最佳叶尖速比,为最大风能利用系数;
根据风力机参数中的转动惯量来设置电磁转矩增益系数,其所用公式为:
Kd=1-α×J (5)
式(5)为最大功率点跟踪控制策略,式(5)中α为Kd的调整系数,该调整系数α的取值范围由转动惯量J的上限值、下限值和电磁转矩增益系数Kd的上限值、下限值确定,按式(3)、式(4)和式(5)调整发电机电磁转矩,即可实现最大功率点跟踪。
2.根据权利要求1所述的一种提升多机型风力机性能的最大功率点跟踪控制方法,其特征在于:应用该方法实现最大功率点跟踪的具体步骤如下:
S1、输入风力机参数:风力机半径、转动惯量、最大风能利用系数、最佳叶尖速比,输入风电场环境参数:风电场空气密度;设定电磁转矩增益系数Kd的调整系数α;
S3、测量发电机实际功率Pe,并根据公式Te=Pe×ω计算发电机实际电磁转矩;
4.根据权利要求3所述的一种提升多机型风力机性能的最大功率点跟踪控制方法,其特征在于:步骤S12中的电磁转矩增益系数Kd的下限值范围为[0.75,0.85]、上限值范围为[0.9,1.0]。
5.根据权利要求2所述的一种提升多机型风力机性能的最大功率点跟踪控制方法,其特征在于:步骤S4中的控制器为PI控制器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710449967.1A CN109139363B (zh) | 2017-06-15 | 2017-06-15 | 一种提升多机型风力机性能的最大功率点跟踪控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710449967.1A CN109139363B (zh) | 2017-06-15 | 2017-06-15 | 一种提升多机型风力机性能的最大功率点跟踪控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109139363A CN109139363A (zh) | 2019-01-04 |
CN109139363B true CN109139363B (zh) | 2020-02-04 |
Family
ID=64829667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710449967.1A Active CN109139363B (zh) | 2017-06-15 | 2017-06-15 | 一种提升多机型风力机性能的最大功率点跟踪控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109139363B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850714B (zh) * | 2019-10-31 | 2023-01-24 | 全球能源互联网研究院有限公司 | 发电机有功功率及风电机组参数计算方法、风电机组模型 |
CN110889781B (zh) * | 2019-12-04 | 2022-05-27 | 浙江大学 | 基于滑模控制的风电机组保性能最大功率跟踪方法 |
CN111162721A (zh) * | 2020-01-08 | 2020-05-15 | 深圳易能电气技术股份有限公司 | 负载参数辨识方法、控制系统、装置及可读存储介质 |
CN113323818B (zh) * | 2021-06-10 | 2024-01-19 | 北京国电思达科技有限公司 | 多类型风机的偏航误差测量方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102434391B (zh) * | 2011-12-27 | 2013-07-17 | 南京理工大学 | 一种基于起始转速调整的改进最大功率点跟踪控制方法 |
KR20130099479A (ko) * | 2012-02-29 | 2013-09-06 | 목포대학교산학협력단 | 풍력발전 제어계통에 대한 센서리스 온라인 신경망 제어방법 |
CN103835878B (zh) * | 2013-04-07 | 2017-05-17 | 南京理工大学 | 基于神经网络优化起始转速的最大功率点跟踪控制方法 |
CN104141591B (zh) * | 2014-07-16 | 2017-01-25 | 南京工程学院 | 风力发电最大功率点跟踪的改进自适应转矩控制方法 |
CN105844544B (zh) * | 2016-04-11 | 2018-02-09 | 南京工程学院 | 基于变系数转矩控制的风力机最大功率点跟踪控制方法 |
-
2017
- 2017-06-15 CN CN201710449967.1A patent/CN109139363B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109139363A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104141591B (zh) | 风力发电最大功率点跟踪的改进自适应转矩控制方法 | |
CN109139363B (zh) | 一种提升多机型风力机性能的最大功率点跟踪控制方法 | |
CN103835878B (zh) | 基于神经网络优化起始转速的最大功率点跟踪控制方法 | |
CN105649875B (zh) | 风力发电机组的变桨控制方法及装置 | |
CN108488035B (zh) | 永磁直驱风力发电机组失速和变桨混合控制方法 | |
CN107947228B (zh) | 基于Markov理论的含风电电力系统随机稳定性分析方法 | |
CN107045574B (zh) | 基于svr的风力发电机组低风速段有效风速估计方法 | |
CN106777499B (zh) | 一种双馈异步风力发电机组的整机动态建模方法 | |
CN109298747A (zh) | 基于iiwo优化的smesc风电系统mppt方法 | |
CN111209638B (zh) | 基于运行攻角偏差修正的低风速风力机叶片气动设计方法 | |
CN106777525B (zh) | 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法 | |
CN109995078A (zh) | 光伏并网机电暂态仿真装置 | |
CN108757306A (zh) | 一种风力发电机控制方法、装置以及风力发电机 | |
CN117028141A (zh) | 一种风电机组的变桨协调优化控制方法 | |
CN112211782B (zh) | 基于起始转速自适应搜索的收缩跟踪区间控制方法 | |
CN105844544B (zh) | 基于变系数转矩控制的风力机最大功率点跟踪控制方法 | |
CN109359426A (zh) | 一种风力机叶片气动参数与控制器参数的联合优化方法 | |
CN101252334A (zh) | 变速恒频风力发电机动态最优能量捕获方法 | |
CN113464378A (zh) | 一种基于深度强化学习的提升风能捕获的转速跟踪目标优化方法 | |
CN104863793A (zh) | 一种根据平均值触发风力发电机变桨动作指令的控制方法 | |
CN112682258A (zh) | 基于Backstepping的大型风力机最大功率点跟踪控制方法 | |
CN105986962B (zh) | 一种风力发电机组的最大风能捕获方法 | |
Nilsson | Numerical computations of wind turbine wakes and wake interaction: Optimization and control | |
Hosseinkhani et al. | Performance Prediction of a SANDIA 17-m Vertical Axis Wind Turbine Using Improved Double Multiple Streamtube | |
CN112884262A (zh) | 一种风电机组载荷适应性确定方法和系统 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |