CN103362735B - 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法 - Google Patents

变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法 Download PDF

Info

Publication number
CN103362735B
CN103362735B CN201210096910.5A CN201210096910A CN103362735B CN 103362735 B CN103362735 B CN 103362735B CN 201210096910 A CN201210096910 A CN 201210096910A CN 103362735 B CN103362735 B CN 103362735B
Authority
CN
China
Prior art keywords
optimum
speed
resisting moment
changing
wind
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
Application number
CN201210096910.5A
Other languages
English (en)
Other versions
CN103362735A (zh
Inventor
贾利民
刘展
雷涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing nenggaopukang measurement and Control Technology Co., Ltd
Original Assignee
BEIJING NEGO AUTOMATION TECHNOLOGY Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BEIJING NEGO AUTOMATION TECHNOLOGY Co Ltd filed Critical BEIJING NEGO AUTOMATION TECHNOLOGY Co Ltd
Priority to CN201210096910.5A priority Critical patent/CN103362735B/zh
Publication of CN103362735A publication Critical patent/CN103362735A/zh
Application granted granted Critical
Publication of CN103362735B publication Critical patent/CN103362735B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Landscapes

  • Wind Motors (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开一种变速变桨风力发电机组基于传动系统传动特性补偿控制的最大功率追踪控制方法。该方法通过在主轴低速端与轮毂连接处安装转矩传感器并将检测的阻力矩反馈实现最优阻力矩闭环控制,有效的使轮毂与主轴低速轴连接处阻力矩跟踪最优阻力矩变化曲线。基于该方法实现的最大功率追踪控制方法相对传统的自寻优控制算法有效的对传动系统的非理想传动特性进行了补偿控制,同时有效的消除了由于风速的随机性和不确定性对最优阻力矩的动态影响,成功的解决了传统的自寻优控制算法进行最大功率追踪时难以准确跟踪最优阻力矩跟踪的情况,有效的提升了风轮的最大功率追踪能力。

Description

变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法
技术领域
本发明属于风力发电技术领域,涉及变速变桨风力发电机组在低于额定风速以下的最大功率追踪控制方法。
背景技术
风能是目前最具有经济价值的可再生能源,根据中国国家发改委能源研究所发布《中国风电发展路线图2050》,中国风电未来40年的发展目标:到2020年、2030年和2050年,风电装机容量将分别达到2亿、4亿和10亿千瓦,到2050年,风电将满足17%的国内电力需求。
目前,风力发电机组已经成为风能利用的主要设备,而变速变桨风力发电机组是主流的风力发电机型。
变速变桨风力发电机组通过风轮吸收风能,利用传动系统将吸收风能传递到发电机轴端,通过发电机将机械能转化为电能,而传动系统的摩擦阻尼损失是传动效率损失的主要来源。
对于大功率变速变桨风力发电机组,传动系统一般情况下可以按照理想刚性环节考虑进行简化建模,若要实现传动系统相对准确的动态分析及控制,通常可以按照柔性的多质量块环节进行建模。
大功率变速变桨风力发电机组在低于额定风速以下的一定风速段内做最大功率追踪控制,风力发电机组通过控制电磁力矩利用风轮的自寻优气动特性实现最大功率追踪。
对于传统的最大功率追踪控制,控制策略中对传动系统考虑的相对比较理想,没有考虑传动系统的传动特性影响,因此在实际控制过程中相对于最优控制力矩轨迹有一定的偏差,不能实现最优最大功率追踪控制。
变速变桨风力发电机组传动系统输入环节中,发电机电磁力矩为可控环节,风轮气动力矩由于风的随机性和不确定性为不可控环节。
发明内容
本发明提出一种变速变桨风力发电机组在额定风速以下的最大功率追踪控制方法,该方法基于对风力发电机组传动系统的传动特性进行补偿控制,通过闭环控制方式实现低于额定风速以下轮毂与主轴低速轴连接处最优阻力矩跟踪,实现风轮风能吸收最大功率追踪控制。
本发明的技术方案是变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法。
如图1所示,在理想认为风力发电机组传动系统是严格的刚性环节且不考虑阻尼的情况下,风力发电机组在低于额定风速以下一定的风速段内通过控制风力机的电磁力矩符合最优力矩变化曲线如图1中的BC段,实现最大功率追踪。
而对于实际的变速变桨风力发电机组,由于其传动系统是柔性环节且存在阻尼,按照传统的控制方法进行控制后,由于传动系统传动特性影响,实际作用于轮毂与主轴低速轴连接处与旋转方向相反的阻力矩难以准确跟踪最优力矩致使风轮未能实现最大功率追踪。
针对变速变桨风力发电机组传动系统特性,在认为变速变桨风力发电机组传动系统为柔性环节且存在摩擦的情况下,如图2所示,变速变桨风力发电机组传动系统可以等效为低速轴和高速轴两质量块环节,按照以下方法建模:
                                                
式中
表示低速轴(风轮端)转速;
表示风轮气动力矩;
表示低速轴(风轮端)转动惯量;
表示轮毂与主轴低速轴连接处与旋转方向相反的阻力矩;
表示高速轴(发电机端)转速;
表示高速轴(发电机端)转动惯量;
表示发电机电磁力矩;
表示传动系统增速比;
表示传动系统传动效率;
表示传动系统等效到低速轴刚性系数;
表示传动系统等效到低速轴阻尼系数。
设状态变量为,输入变量为,输出为,则传动系统可以表示为一个三阶模型:
 
其中
 
表示低速轴等效到高速轴后的转动惯量。
则系统传递函数可以表示为:
系统传递函数框图如图3所示,可以表示为下式:
 
一般情况下,低速轴等效到高速轴后的转动惯量也明显小于低速轴(风轮端)转动惯量,因此从稳态增益的角度而言,轮毂与主轴低速轴连接处与旋转方向相反的阻力矩主要由产生,由于系统中风轮气动力矩由风轮气动力矩产生,为不可控变量,可以看做为外部输入扰动。
针对变速变桨风力发电机组电磁力矩响应环节建模,通常电磁力矩的响应速度很快,可近似认为瞬时响应,当需要精确控制时,电磁力矩响应模型可近似用一阶惯性环节表示如下:
基于以上模型可以设计出单回路闭环负反馈控制系统如图4所示,在该控制系统中,控制目标为最优阻力矩,最优阻力矩可以按照下式计算:
式中
表示最优阻力矩系数;
表示空气密度;
表示风轮扫风半径;
表示风轮在叶片低于额定风速工作时时(桨距角最小)的最大风能利用系数;
表示在叶片低于额定风速工作时(桨距角最小)的最佳叶尖速比。
单位负反馈量为通过传感器检测的轮毂与主轴低速端连接处的转矩测量值。
控制目标偏差可以按照下式计算:
式中
表示传感器检测的轮毂与主轴低速端连接处的转矩测量值。
控制器输入为控制目标偏差,控制器输出为电磁力矩设定值,控制器输出信号至电磁力矩响应环节
控制器可以采用常用的PID控制器,也可以采用其它的控制器形式,当采用PID控制器时,控制器传递函数可以表示为:
控制率为:
式中
表示PID控制器比例环节系数;
表示PID控制器积分环节系数;
表示PID控制器微分环节系数。
电磁力矩响应环节的输入为电磁力矩设定值,输出为真实电磁力矩
传动环节的输入变量为,其中风轮气动力矩为外部不可控输入,可以看作为扰动,发电机电磁力矩为可控变量,输出为轮毂与主轴低速轴连接处与旋转方向相反的阻力矩可以通过传感器测量,并将测量值单位负反馈回输入端。
附图说明
图1变速变桨风力发电机组最优力矩控制曲线示意图
图2 传动系统两质量块建模示意图
图3传动系统传递函数框图
图4基于最优阻力矩跟踪的最大功率追踪闭环控制框图
具体实施方式
根据变速变桨风力发电机组叶片及轮毂质量距离旋转中心分布计算传动系统低速轴总的转动惯量;根据变速变桨风力发电机组发电机转子质量距离旋转中心分布数据计算传动系统高速轴总的转动惯量,一般而言传动轴或齿轮箱相对转动惯量很小可以近似忽略。
获取变速变桨风力发电机组叶片翼型数据,根据动量-叶素定理、有限元计算方法或者依据相关商业计算软件计算风轮在低于额定风速以下的最优阻力矩参数
传动系统建模相关数据可以根据传动系统的传动方式查询相关机组资料获取传动效率、等效于低速轴的刚性系数和等效于低速轴的摩擦系数数据,基于数据进行机理建模;也可以通过测量相关辨识数据用模型拟合参数辨识的方式直接获取传动系统辨识模型。
由于变速变桨风力发电机组在不同工作点的传动效率、等效于低速轴的刚性系数和等效于低速轴的摩擦系数等数据可能有一定差异,因此当系统参数偏差较大时可以考虑在不同工作区间建立分段线性化近似模型,在不同的工作区间基于不同的分段线性化模型设计不同参数的控制器进行控制从而获取更好的控制效果。
当电磁力矩响应环节时间常数与传动模型输出相应时间常数相差不大时,必须考虑电磁力矩动态响应特性,可以通过监测采集相关辨识数据用模型拟合参数辨识的方式直接获取力矩响应辨识模型。
轮毂与主轴低速轴连接处与旋转方向相反的阻力矩可以通过在轮毂与主轴低速轴连接处安装动态转矩传感器精确测量,通常为了确保机组运行可靠性可以安装两套动态转矩传感器通过冗余提高系统可靠性。
PID参数调节方法可以采用工业上常用的经验法、Z-N法或者通过离线相关智能搜索算法如遗传算法、蚁群算法等优化进行整定,也可以通过经典控制理论中的零极配置方法进行参数整定。
本发明针对变速变桨风力发电机组由于其传动系统的非理想传动特性导致在低于额定风速以下机组按照传统的自寻优控制算法进行最大功率追踪时难以准确跟踪最优阻力矩跟踪的情况,提出了一种轮毂与主轴低速轴连接处最优阻力矩闭环控制方法,该方法可以有效的使轮毂与主轴低速轴连接处阻力矩跟踪最优阻力矩变化曲线,基于该方法实现的最大功率追踪控制方法相对传统的自寻优控制算法有效的对传动系统的非理想传动特性进行了补偿控制,同时有效的消除了由于风速的随机性和不确定性对最优阻力矩的动态影响,有效的提升了风轮的最大功率追踪能力。

Claims (7)

1.变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法,其特征是该方法基于对风力发电机组传动系统的传动特性进行补偿控制,通过闭环控制方式实现低于额定风速以下轮毂与主轴低速轴连接处最优阻力矩跟踪,实现风轮风能吸收最大功率追踪控制;传动系统可以等效为低速轴和高速轴两质量块环节,按照以下方法建模:
dω b dt = T b J b - T J b dω g dt = ηT NJ g - T e J g dT dt = K ( ω b - ω g N ) + B ( dω b dt - 1 N dω g dt )
式中
ωb表示低速轴转速;
Tb表示风轮气动力矩;
Jb表示低速轴转动惯量;
T表示轮毂与主轴低速轴连接处与旋转方向相反的阻力矩;
ωg表示高速轴转速;
Jg表示高速轴转动惯量;
Te表示发电机电磁力矩;
N表示传动系统增速比;
η表示传动系统传动效率;
K表示传动系统等效到低速轴刚性系数;
B表示传动系统等效到低速轴阻尼系数;
设状态变量为 x = ω b ω g T , 输入变量为 u = T b T e , 输出为y=T,则传动系表示为一个三阶模型:
其中
J t = N 2 J g η
Jt表示低速轴等效到高速轴后的转动惯量;
则系统传递函数表示为:
y ( s ) = [ G 1 ( s ) , G 2 ( s ) ] u ( s ) G 1 ( s ) = BJ t s + KJ t J b J t s 2 + B ( J b + J t ) s + K ( J b + J t ) G 2 ( s ) = BJ b J t s + KJ b J t NJ g J b J t s 2 + BNJ g ( J b + J t ) s + KNJ g ( J b + J t )
2.根据权利要求1所述变速变桨风力发电机组最优阻力矩跟踪的最大功率追踪控制方法,其特征是电磁力矩响应模型可近似用一阶惯性环节表示如下:
G 3 ( s ) = 1 τs + 1
τ表示电磁力矩响应时间常数。
3.根据权利要求1所述变速变桨风力发电机组最优阻力矩跟踪的最大功率追踪控制方法,其特征是最大功率追踪控制系统中控制目标为最优阻力矩Topt,最优阻力矩按照下式计算:
T opt = K opt ω b 2 K opt = 1 2 ρπ R 5 C P max λ C P max 3
式中
Kopt表示最优阻力矩系数;
ρ表示空气密度;
R表示风轮扫风半径;
CPmax表示风轮在叶片低于额定风速工作时的最大风能利用系数;
表示在叶片低于额定风速工作时的最佳叶尖速比。
4.根据权利要求1所述变速变桨风力发电机组最优阻力矩跟踪的最大功率追踪控制方法,其特征是控制目标偏差e(t)按照下式计算:
( t ) = T opt ( t ) - T ‾ ( t )
式中
Topt表示最优阻力矩;
表示传感器检测的轮毂与主轴低速端连接处的转矩测量值。
5.根据权利要求1所述变速变桨风力发电机组最优阻力矩跟踪的最大功率追踪控制方法,其特征是控制器输入为控制目标偏差e(t),控制器输出为电磁力矩设定值控制器输出信号至电磁力矩响应环节G3(s);控制器可以采用常用的PI D控制器,也可以采用其它的控制器形式,当采用PI D控制器时,控制器传递函数表示为:
c ( s ) = K P + K I + 1 s + K D s
控制率为:
T ‾ e ( t ) = K P e ( t ) + K I ∫ 0 t e ( t ) dt + K D ed ( t ) dt
式中
KP表示PID控制器比例环节系数;
KI表示PID控制器积分环节系数;
KD表示PID控制器微分环节系数。
6.根据权利要求1所述变速变桨风力发电机组最优阻力矩跟踪的最大功率追踪控制方法,其特征是电磁力矩响应环节G3(s)的输入为电磁力矩设定值输出为真实电磁力矩Te
7.根据权利要求1所述变速变桨风力发电机组最优阻力矩跟踪的最大功率追踪控制方法,其特征是传动环节的输入变量为 T b T e , 其中风轮气动力矩Tb为外部不可控输入,可以看作为扰动,发电机电磁力矩Te为可控变量,输出为轮毂与主轴低速轴连接处与旋转方向相反的阻力矩T可以通过传感器测量,并将测量值单位负反馈回输入端。
CN201210096910.5A 2012-04-05 2012-04-05 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法 Active CN103362735B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210096910.5A CN103362735B (zh) 2012-04-05 2012-04-05 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210096910.5A CN103362735B (zh) 2012-04-05 2012-04-05 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法

Publications (2)

Publication Number Publication Date
CN103362735A CN103362735A (zh) 2013-10-23
CN103362735B true CN103362735B (zh) 2015-10-28

Family

ID=49364823

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210096910.5A Active CN103362735B (zh) 2012-04-05 2012-04-05 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法

Country Status (1)

Country Link
CN (1) CN103362735B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104141591B (zh) * 2014-07-16 2017-01-25 南京工程学院 风力发电最大功率点跟踪的改进自适应转矩控制方法
CN104329224B (zh) * 2014-11-13 2017-02-15 湖南世优电气股份有限公司 一种直驱型永磁风力发电机组变速率变桨系统及变桨方法
CN105649877B (zh) * 2016-01-13 2018-11-06 湖南世优电气股份有限公司 一种大型风电机组的蚁群pid独立变桨控制方法
CN106286151B (zh) * 2016-08-26 2019-03-22 中车株洲电力机车研究所有限公司 一种风电机组低速轴扭转载荷监测方法及载荷分析方法
CN109751185B (zh) * 2019-01-04 2020-01-07 上海交通大学 一种变速风电机组的分频控制方法及控制系统
CN111188732B (zh) * 2020-01-17 2022-05-13 湖南工业大学 一种风力发电变桨鲁棒容错控制方法
CN111219295B (zh) * 2020-03-02 2020-12-08 南京信息工程大学 一种风力发电机风功率曲线补偿控制装置及方法
CN112462798B (zh) * 2020-12-04 2021-05-28 三生万物(北京)人工智能技术有限公司 一种无人机及提高无人机航线飞行性能的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1538262A (zh) * 2003-03-11 2004-10-20 欧姆龙株式会社 最大电力跟踪控制装置
CN101012809A (zh) * 2007-02-08 2007-08-08 上海交通大学 基于风向标和输出功率的风力机偏航控制方法
WO2010044163A1 (ja) * 2008-10-16 2010-04-22 三菱重工業株式会社 風力発電システム、及びその制御方法
CN101793235A (zh) * 2010-04-15 2010-08-04 哈尔滨工业大学 具有能量预测功能的最大功率跟踪型风力发电装置及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005025026A1 (de) * 2003-09-03 2005-03-17 Repower Systems Ag Verfahren zum betrieb bzw. regelung einer windenergieanlage sowie verfahren zur bereitstellung von primärrefelleistung mit windenergieanlagen

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1538262A (zh) * 2003-03-11 2004-10-20 欧姆龙株式会社 最大电力跟踪控制装置
CN101012809A (zh) * 2007-02-08 2007-08-08 上海交通大学 基于风向标和输出功率的风力机偏航控制方法
WO2010044163A1 (ja) * 2008-10-16 2010-04-22 三菱重工業株式会社 風力発電システム、及びその制御方法
CN101793235A (zh) * 2010-04-15 2010-08-04 哈尔滨工业大学 具有能量预测功能的最大功率跟踪型风力发电装置及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
风机变桨距非线性预测控制算法;雷涛等;《计算机仿真》;20110131;第28卷(第1期);第319-322页 *

Also Published As

Publication number Publication date
CN103362735A (zh) 2013-10-23

Similar Documents

Publication Publication Date Title
CN103362735B (zh) 变速变桨风力发电机组基于最优阻力矩跟踪的最大功率追踪控制方法
CN103362738B (zh) 变速变桨风力发电机组基于前馈解耦控制的最大功率追踪控制方法
CN102797629B (zh) 一种风电机组的控制方法、控制器及其控制系统
CN103850876B (zh) 一种适用于无载荷测量的风电机组独立变桨控制方法
CN101776043B (zh) 基于误差补偿模型的风电机组最大风能捕获控制方法
CN104314757A (zh) 一种风力发电机组偏航控制方法及系统
CN102072083B (zh) 一种双馈风力发电系统最大风能跟踪方法
CN103244350A (zh) 一种风力发电机组最佳叶尖速比跟踪控制方法
KR20090083371A (ko) 풍력 발전 장치, 풍력 발전 시스템 및 풍력 발전 장치의 발전 제어 방법
CN103441537A (zh) 配有储能电站的分散式风电场有功优化调控方法
CN103061980A (zh) 基于激光测风雷达的风力发电机组的前馈控制系统及其控制方法
CN103807096A (zh) 风力涡轮机及其控制方法
CN107061158B (zh) 一种低风速下风力发电机的预测与跟踪控制方法
CN103758697B (zh) 一种基于风机有效跟踪区间的改进最大功率跟踪控制方法
CN104074679A (zh) 一种变速变桨距风电机组全风速限功率优化控制方法
CN103939286A (zh) 变速恒频风电机组变速变桨距联合控制方法
CN104797813A (zh) 使用入射风速估算的风轮机控制方法
CN101592126B (zh) 直驱永磁同步风力发电机组风能捕获跟踪控制方法
CN104454347A (zh) 一种独立变桨距风力发电机桨距角的控制方法
CN103362736B (zh) 变速变桨风力发电机组基于内模控制的最大功率追踪控制方法
CN103590975A (zh) 一种环境友好型的风电机组控制系统及方法
CN105134490A (zh) 提高风力发电机组经济性的方法
CN114033617B (zh) 一种控制参量自适应调整的可控风力发电方法及系统
Suryanarayanan et al. Control of large wind turbines: Review and suggested approach to multivariable design
CN103375332B (zh) 变速变桨风力发电机组最优阻力矩动态优化方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20200402

Address after: 100071 207, 2 / F, building 3, yard 128, South Fourth Ring West Road, Fengtai District, Beijing

Patentee after: Beijing nenggaopukang measurement and Control Technology Co., Ltd

Address before: 100044 Beijing city Haidian District Xizhimen, Beijing Jiaotong University Simon Jiaotong University for six storey building

Patentee before: BEIJING NEGO AUTOMATION TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right