CN110454329A - 一种风电机组桨距角控制方法 - Google Patents
一种风电机组桨距角控制方法 Download PDFInfo
- Publication number
- CN110454329A CN110454329A CN201910830549.6A CN201910830549A CN110454329A CN 110454329 A CN110454329 A CN 110454329A CN 201910830549 A CN201910830549 A CN 201910830549A CN 110454329 A CN110454329 A CN 110454329A
- Authority
- CN
- China
- Prior art keywords
- wind speed
- pitch angle
- wind turbine
- control signal
- model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000012549 training Methods 0.000 claims abstract description 37
- 238000013178 mathematical model Methods 0.000 claims abstract description 26
- 230000004044 response Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 3
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 5
- 238000011217 control strategy Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 238000012938 design process Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- OVBPIULPVIDEAO-LBPRGKRZSA-N folic acid Chemical compound C=1N=C2NC(N)=NC(=O)C2=NC=1CNC1=CC=C(C(=O)N[C@@H](CCC(O)=O)C(O)=O)C=C1 OVBPIULPVIDEAO-LBPRGKRZSA-N 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 229940064302 folacin Drugs 0.000 description 1
- 235000019152 folic acid Nutrition 0.000 description 1
- 239000011724 folic acid Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012827 research and development Methods 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
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/022—Adjusting aerodynamic properties of the blades
- F03D7/0224—Adjusting blade pitch
-
- 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)
- Physics & Mathematics (AREA)
- Fluid Mechanics (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)
Abstract
本发明涉及一种风电机组桨距角控制方法,建立基于ELM的有效风速估计模型:获取一段时间内有效风速信息和相应输出数据训练有效风速估计模型;获取反馈控制信号:建立风电机组数字模型;定义调节误差,求动态特性,求反馈控制信号的初步表达式;对风电机组数学模型的不确定及干扰项进行估计;对初步表达式变换处理获得反馈控制信号;获取前馈控制信号:获取一段时间内有效风速信息和机组相应输出数据中的桨距角信号,训练风电机组动态非线性风速‑桨距角关系模型,该模型能够在线实时给出前馈控制信号。本发明针对传统桨距角控制器存在的由于前馈桨距角信息不准确导致的载荷优化效果不佳和由于需要机组线性模型而带来的计算复杂的问题。
Description
技术领域
本发明属于风力发电机组控制技术领域,特别涉及一种风电机组桨距角控制方法。
背景技术
在过去的几十年里,风力发电技术在全球范围内得到了迅速发展。根据全球风能协会发布的全球风电展望报告,如果能保持当前的增长速度,到2030年全球能源消耗的30%将由风电提供。然而,风电度电成本的居高不下以及风能的间隙性、不确定性仍然是限制风电快速发展和大规模并网的主要因素。作为风电技术领域的关键技术,智能控制是解决上述问题的有效途径,而风电机组的智能控制技术一直以来都被德国、丹麦等西方发达国家所垄断,因此,研发具有自主知识产权的智能风电控制技术具有重要的战略意义和实践价值。
与定速风电机组相比,变速风电机组具有更高的风能捕获效率和更低的机械载荷,因此,变速风电机组已经成为当前风电产业中的主流机组。变速风电机组主要有两个运行区域:低风速运行区域(风速大于切入风速小于额定风速)和高风速运行区域(风速大于额定风速小于切出风速)。最大风能捕获是低风速运行区域的主要控制目标,这一目标可以通过调节电磁转矩从而使得机组运行在最佳转速上来实现。当机组运行在高风速工作区域时,通常将电磁转矩维持在额定值,将桨距角视为控制信号,同时,将电磁转矩设置为额定值,而控制目标包括两个方面:首先,为了保证整个风电系统的安全运行,转速和发电功率要维持在额定值附近;其次,为了延长机组的服役寿命,机组关键部件(包括叶片、轮毂、塔架、变桨轴承和变桨系统等)的载荷需要进行优化。值得注意的是,因为机组大部分的疲劳载荷是在高风速运行区域产生的,所以风电机组的载荷优化对于高风速运行区域意义更大,以叶片为例,高风速运行区域产生的疲劳载荷占其全部疲劳载荷的95%。
为了在不同的运行区域实现不同的控制目标,风电机组在不同的风速区间采用不同的控制策略。对于高风速运行区域,目前工业上普遍采用PI控制方法。然而,仅采用负反馈的PI控制策略,只能在系统调节误差产生后进行作用,同时,现代大型机组的惯量较大,导致系统响应较慢,因此,PI控制方法的转速/功率调节效果并不理想,且不具备载荷优化的功能。为了解决这一问题,学者们在现有反馈环的基础上加入了前馈控制策略,通过利用有效风速信息提高调节性能和载荷优化效果,目前已有的桨距角前馈控制方法主要包括两类:查表法和传递函数法。查表法主要思想是使用有效风速信息查询风速和桨距角之间的静态关系数据表,从而给出前馈桨距角的大小,然而,风速和桨距角之间的静态关系表不能准确描述机组动态运行过程中的有效风速和桨距角之间的对应关系;传递函数法的主要思想是利用桨距角与转速调节误差、风速与转速调节误差之间的开环传递函数计算出前馈补偿函数,该补偿函数能够根据风速信息计算出前馈桨距角的值,然而,这种方法需要对机组的数学模型在某个工作点处进行线性化,不仅计算复杂,而且会丢失机组动态的非线性部分,进而影响其控制性能。
通过上述分析可以看出,有效风速信息对于机组前馈控制来讲是必备的。为了克服这一困难,风电技术研究人员提出了很多有效风速估计方法,这些方法可以分为两类:基于模型的方法和基于数据的方法,前者依赖于精确的机组数学模型或者机械功率表达式,尽管在仿真中能够取得不错的结果,但是在实际场测中效果却大打折扣;后者由于没有充分考虑机组不同的运行状态特征,导致其在线使用效果较差。
发明内容
本发明的目的在于提供一种风电机组桨距角控制方法,本发明针对传统前馈-反馈桨距角控制器存在的由于前馈桨距角信息不准确导致的载荷优化效果不佳和由于需要机组线性模型而带来的计算复杂的问题。
为实现上述目的,本发明所采取的技术方案是:
一种风电机组桨距角控制方法,其特征在于包括:
建立基于ELM的有效风速估计模型:获取一段时间内风电机组的有效风速信息和相应输出数据,训练基于ELM的有效风速估计模型,有效风速估计模型能够在线实时给出有效风速估计值;
获取反馈控制信号:建立风电机组数字模型;定义调节误差,求得调节误差的动态特性,根据调节误差的动态特性求得反馈控制信号的初步表达式;根据UDE原理对风电机组数学模型的不确定及干扰项进行估计,获得不确定及干扰项估计值;先将不确定及干扰项估计值带入初步表达式,再对初步表达式进行拉普拉斯变换和拉普拉斯反变换获得反馈控制信号的最终表达式;
获取前馈控制信号:获取一段时间内机组有效风速信息和机组相应输出数据中的桨距角信号,训练基于ELM的风电机组动态非线性风速-桨距角关系模型,动态非线性风速-桨距角关系模型能够在线实时给出前馈控制信号。
进一步地,所述建立基于ELM的有效风速估计模型包括:
(1)分别通过LIDAR测风系统和SCADA系统获得一段时间内风电机组的有效风速信息Y和相应输出数据X;
(2)将相应输出数据X进行归一化处理,作为有效风速估计模型的训练特征集,将有效风速信息Y作为有效风速估计模型的训练目标值,将训练特征集和训练目标值作为有效风速估计模型的训练集T;
(3)通过训练集T训练得到基于ELM的有效风速估计模型。
进一步地,所述建立风电机组数字模型包括:
(1)对风电机组的结构进行分析,建立风电机组的初步数字模型;
其中,ωr是风轮角速度,Jt是系统总惯量,Ta是气动转矩,Kt是系统总阻尼,Tg是等效发电机电磁转矩,dt是外部干扰;
(2)对所述风电机组的初步数字模型进行改写,得到风电机组数字模型:
其中,是风电机组数学模型的已知项;D(t)=Δ(ωr,β)+d(t)风电机组数学模型的不确定及干扰项, 是有效风速估计值。
进一步地,所述定义调节误差包括:
(1)定义调节误差e=ωr-ωd,其中ωd表示额定风轮转速;
(2)根据风电机组数学模型,求得调节误差的动态特性表示如下:
e的动态特性满足k>0;由此,可以求得反馈控制信号的初步表达式如下:
进一步地,所述获得不确定及干扰项估计值包括:
使用UDE原理,通过低通滤波器Gf(s)对风电机组数学模型中的不确定及干扰项D(t)进行逼近,记录低通滤波器Gf(s)的冲击响应为gf(t),得到D(t)的不确定及干扰项估计值 的表达式如下:
其中,“*”表示卷积。
进一步地,所述获得反馈控制信号的最终表达式包括:
将不确定及干扰项估计值代替反馈控制信号的初步表达式中的D(t),得到:
将式①进行拉普拉斯变换可得:
β(s)=f(s)+(sωr(s)-f(s)+β(s))Gf(s)+ke(s) ②
其中β(s)、f(s)、ωr(s)和e(s)分别表示β、ωr和e的拉普拉斯变换,继续整理式②,可得:
将式③进行拉普拉斯反变换可得反馈控制信号的最终表达式如下:
进一步地,所述Gf(s)选取为一阶低通滤波器τ是滤波器的时间常数。
进一步地,所述获取前馈控制信号包括:获取一段时间内机组有效风速信息和机组相应输出数据中的桨距角信号,使用当前时刻以及之前连续三个时刻点的有效风速值作为训练特征,将当前时刻的桨距角值作为训练目标训练基于ELM的风电机组动态非线性风速-桨距角关系模型fELMFF,从动态非线性风速-桨距角关系模型fELMFF中获取在线实时前馈控制信号:
其中是当前时刻的有效风速估计值,和是从当前时刻连续向前数三个时刻点的有效风速估计值,其中两个连续时刻点之间的时间长度是控制周期。
进一步地,还包括获取风电机组桨距角控制信号:
风电机组桨距角控制信号=反馈控制信号+前馈控制信号。
本发明的有益效果是:本发明针对传统前馈-反馈桨距角控制器存在的由于前馈桨距角信息不准确导致的载荷优化效果不佳和由于需要机组线性模型而带来的计算复杂的问题,设计基于ELM的动态非线性风速-桨距角前馈控制机制,该机制能够较准确地给出前馈桨距角的数值,在提高转速\发电功率的恒值调节性能的同时,减少反馈环带来的变桨系统的频繁动作,同时降低机组关键部件的载荷;针对目前现有风速估计方法由于依赖于系统精确数学模型导致的实际场测效果差、在线效果差的问题,设计基于ELM的有效风速估计方法,合理选择机组相关输出,摆脱对系统模型的依赖,训练过程简单,得到的风速估计模型能够以较高的精度实时地给出风速估计值;针对现有桨距角控制器设计过程复杂、鲁棒性不佳的问题,设计基于有效风速信息和UDE(Uncertainty and DisturbanceEstimator)的反馈控制策略,该反馈控制算法设计过程简单,需要调试的参数少。该风电机组的桨距角控制算法,能够实现在提高机组转速\发电功率的恒值调节性能的同时降低机组关键部件的疲劳载荷,进而降低机组故障率,延长风电机组的服役寿命,对模型不确定性和环境干扰具有较好鲁棒性,控制效果好,简单易实施,提高风电场的经济效益,降低风电场的运维成本。
附图说明
图1为风电机组桨距角示意图。
图2为基于有效风速估计与UDE的桨距角控制方法设计流程图。
图3为基于有效风速估计与UDE的桨距角器框图。
图4为风速真实值、风速估计值对比图。
图5为风轮转速对比图。
图6为发电功率对比图。
图7为桨距角对比图。
图8为桨距角变化率对比图。
具体实施方式
为了更好地理解本发明,下面结合实施例和附图对本发明的技术方案做进一步的说明。
本发明提供的一种基于有效风速估计和UDE的新型风电机组桨距角控制方法,包括:
一、建立基于ELM的有效风速估计模型:
S101、分别通过LIDAR测风系统和SCADA系统获得一段时间内机组的有效风速信息和相应输出数据,
得到的机组有效风速信息用Y表示,Y=[vti],i=1,2,3,...,l,
得到的机组相应输出数据用X'表示,X'=[x'(i,j)],i=1,...,l,j=1,2,3,
其中l是样本个数,x'(i,:)是机组相应输出数据的一次采样,x'(i,:)表达式为:
其中,β是桨距角,是桨距角加速度,afa是塔架前后加速度。
S102、将上述X'进行归一化,为了避免由异常值带来的数据整体偏斜问题,本发明采用均值标准差归一化方式,具体计算公式如下:
其中,x'(:,j)是中X'的列分量,mean(x'(:,j))和std(x'(:,j))分别是x'(:,j)的均值和方差,x(:,j)是有效风速估计模型的训练特征集X(X=[x(i,j)],i=1,...,l,j=1,2,3)中的列分量,将有效风速估计模型的训练目标集选择为Y,则有效风速估计模型的训练集可以表示为:T=[X,Y]。
S103、运用上述训练集T训练基于ELM的有效风速估计模型,此处的ELM模型指的是普遍使用的极限学习机模型,该ELM包括输入层,隐含层和输出层三层,输入层有三个节点,隐含层的节点个数建议设置为9,输出层的个数为1,隐含层的激活函数选择为sigmoid函数。训练完成后的ELM模型记为fELM,该模型能够在线实时给出有效风速估计值:
二、获取反馈控制信号:
S201、对风电机组的结构进行分析,根据叶素理论和物理学基本定律,得到如下的风电机组的初步数学模型:
其中,ωr是风轮角速度,Jt是系统总惯量,Ta是气动转矩,Kt是系统总阻尼,Tg是等效发电机电磁转矩,dt是外部干扰。
Ta的表达式为:
其中,ρ是空气密度,R是风轮半径,是叶尖速比,v是有效风速,β是桨距角,其定义为风轮扫略面和弦线的夹角,如图1所示。Cp为功率系数,能够衡量风能捕获能力,其表达式可以写为:
S202、为了简化控制器设计过程,将所述风电机组的初步数学模型进行改写,得到改写后的风电机组数学模型:
其中,
是风电机组数学模型的已知项,
D(t)=Δ(ωr,β)+d(t)风电机组数学模型的不确定及干扰项,
是有效风速估计值。
S203、为了设计反馈控制器,定义调节误差e=ωr-ωd,其中ωd表示额定风轮转速。根据改写后的风电机组数学模型,求得调节误差的动态特性表示如下:
为了使得调节误差收敛到零(系统稳定),e的动态特性必须满足k>0,由此,可以求得反馈控制信号的初步表达式如下:
由于D(t)未知的,控制信号的表达式并不能在实际中使用。
S204、使用UDE(不确定和干扰估计因子)原理,选择一个带宽合适的严格正则的低通滤波器Gf(s)对改写后的风电机组数学模型中的信号D(t)进行逼近,记低通滤波器Gf(s)的冲击响应为gf(t),使用改写后的风电机组数学模型,可以得到D(t)的估计值的表达式如下:
其中,“*”表示卷积。
S205、使用步骤S204中的表达式代替步骤S203中反馈控制信号的初步表达式中的D(t),则可以得到:
将式①进行拉普拉斯变换可得:
β(s)=f(s)+(sωr(s)-f(s)+β(s))Gf(s)+ke(s)②
其中β(s)、f(s)、ωr(s)和e(s)分别表示β、ωr和e的拉普拉斯变换,继续整理式②,可得:
将式③进行拉普拉斯反变换可得反馈控制信号的最终表达式如下:
为了清晰起见,本发明中将反馈控制信号记为βfb,而将前馈控制信号记为βff。其中Gf(s)选取为一阶低通滤波器τ是滤波器的时间常数。
三、获取前馈控制信号:
为了获得前馈控制信号βff,使用步骤S101中获得的机组有效风速信息和机组相应输出数据中的桨距角信号,训练基于ELM的风电机组动态非线性风速-桨距角关系模型,此处的ELM模型值的是普遍使用的极限学习机模型,该ELM包括输入层,隐含层和输出层三层,输入层有4个节点,隐含层的节点个数建议设置为12,输出层的个数为1,隐含层的激活函数选择为sigmoid函数,即本方案使用当前时刻以及之前连续三个时刻点的风速值作为训练特征,而将当前时刻的桨距角值作为训练目标。训练完成后的动态非线性风速-桨距角关系模型记为fELMFF,该模型能够在线实时给出前馈控制信号:
其中是当前时刻的有效风速估计值,和是从当前时刻连续向前数三个时刻点的有效风速估计值,其中两个连续时刻点之间的时间长度是控制周期。
四、获取风电机组桨距角控制信号
求得的最终表达式如下:β=βff+βfb
本发明的原理:图2是基于有效风速估计与UDE的桨距角控制方法设计流程图。使用机组历史输出数据训练基于ELM的有效风速估计模型,该模型能够在线给出有效风速估计值序列。对风电机组的结构进行分析,得到风电机组的数学模型,为了反馈控制器设计的方便,对该数学模型进行改写,得到改写后的风电机组数学模型。使用改写后的风电机组数学模型和有效风速估计模型,设计基于UDE原理的反馈桨距角控制信号。使用历史风速和桨距角数据训练得到风速-桨距角动态关系模型,根据有效风速信息得到前馈桨距角控制信号;将反馈桨距角控制信号和前馈桨距角控制信号相加,得到最终的桨距角控制信号。图3为基于有效风速估计与UDE的桨距角器框图。机组的实时输出经过归一化后,输入到有效风速估计模型中,风速估计模型实时给出有效风速估计值,基于ELM的风速-桨距角动态关系模型根据有效风速估计值给出前馈桨距角的大小,而UDE控制器利用反馈误差和风速估计值的大小给出反馈桨距角的大小,反馈桨距角和前馈桨距角相加得到最终的桨距角控制信号。
实施例
为量化比较本发明提出的基于有效风速估计的前馈-反馈控制器和传统的(基于查表法)前馈-反馈控制器的载荷优化效果,需要对机组关键部件的载荷评价指标进行说明。对于桨距系统而言,本发明使用如下的指标(PF,Pitch Fatigue)进行评价:
其中,β(k)是风电机组的桨距角,上述指标的单位是deg/s,它反映了在机组运行时间段Tβ秒内,平均每秒桨距角的变化量,能够比较好的评价桨距系统的载荷。对于机组的其他关键部件,包括叶片(BR,Blade Root,)、轮毂(FH,Fixed Hub,)、塔架(TB,TowerBottom)、变桨轴承(PB,Pitch Bearing),则使用DEL(Damage Equivalent Loads,损伤等小载荷)来评价。
本实施例使用工业上普遍采用的风电技术开发软件GH Bladed仿真平台,将本发明的方法与目前工业上普遍采用的基于查表法的前馈-反馈桨距角控制器(以下简称传统方法)进行对比,说明本方法的有效性和创新性。
实施例中使用1.5MW三叶片水平轴变速风力发电机组模型,其主要参数如下表所示:
图1为风电机组桨距角示意图,根据叶素理论,风电机组的桨距角定义为风轮扫掠面与叶片弦线的夹角。
图2是基于有效风速估计与UDE的桨距角控制方法设计流程图。
图3为基于有效风速估计与UDE的桨距角器框图。
图4为风速真实值、风速估计值对比图。
图5为风轮转速对比图。经计算,传统方法的风轮转速超调量是8.57%,本方法的风轮转速超调量是5.85%,降低31.74%。
图6为发电功率对比图。经计算,传统方法的发电功率超调量是8.73%,前馈-反馈控制的发电功率超调量是5.96%,降低31.73%。
图7为桨距角对比图。
图8为桨距角变化率对比图。经计算,传统方法的桨距系统载荷评价指标PF=0.8586,本方法的桨距系统载荷评价指标PF=0.7683,降低10.52%。机组的其他关键部件,包括叶片(BR)、轮毂(FH)、塔架(TB)、变桨轴承(PB)的DEL对比如下表所示,其中表中的DEL数值均以传统方法为基准进行了正规化,可见机组各关键部件的疲劳载荷都得到了不同程度的降低。
以上说明仅为本发明的应用实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明申请专利范围所作的等效变化,仍属本发明的保护范围。
Claims (9)
1.一种风电机组桨距角控制方法,其特征在于包括:
建立基于ELM的有效风速估计模型:获取一段时间内风电机组的有效风速信息和相应输出数据,训练基于ELM的有效风速估计模型,有效风速估计模型能够在线实时给出有效风速估计值;
获取反馈控制信号:建立风电机组数字模型;定义调节误差,求得调节误差的动态特性,根据调节误差的动态特性求得反馈控制信号的初步表达式;根据UDE原理对风电机组数学模型的不确定及干扰项进行估计,获得不确定及干扰项估计值;先将不确定及干扰项估计值带入初步表达式,再对初步表达式进行拉普拉斯变换和拉普拉斯反变换获得反馈控制信号的最终表达式;
获取前馈控制信号:获取一段时间内机组有效风速信息和机组相应输出数据中的桨距角信号,训练基于ELM的风电机组动态非线性风速-桨距角关系模型,动态非线性风速-桨距角关系模型能够在线实时给出前馈控制信号。
2.根据权利要求1所述的一种风电机组桨距角控制方法,其特征在于,所述建立基于ELM的有效风速估计模型包括:
(1)分别通过LIDAR测风系统和SCADA系统获得一段时间内风电机组的有效风速信息Y和相应输出数据X;
(2)将相应输出数据X进行归一化处理,作为有效风速估计模型的训练特征集,将有效风速信息Y作为有效风速估计模型的训练目标值,将训练特征集和训练目标值作为有效风速估计模型的训练集T;
(3)通过训练集T训练得到基于ELM的有效风速估计模型。
3.根据权利要求1所述的一种风电机组桨距角控制方法,其特征在于,所述建立风电机组数字模型包括:
(1)对风电机组的结构进行分析,建立风电机组的初步数字模型;
其中,ωr是风轮角速度,Jt是系统总惯量,Ta是气动转矩,Kt是系统总阻尼,Tg是等效发电机电磁转矩,dt是外部干扰;
(2)对所述风电机组的初步数字模型进行改写,得到风电机组数字模型:
其中,是风电机组数学模型的已知项;D(t)=Δ(ωr,β)+d(t)风电机组数学模型的不确定及干扰项, 是有效风速估计值。
4.根据权利要求3所述的一种风电机组桨距角控制方法,其特征在于,所述定义调节误差包括:
(1)定义调节误差e=ωr-ωd,其中ωd表示额定风轮转速;
(2)根据风电机组数学模型,求得调节误差的动态特性表示如下:
e的动态特性满足k>0;由此,可以求得反馈控制信号的初步表达式如下:
5.根据权利要求4所述的一种风电机组桨距角控制方法,其特征在于,所述获得不确定及干扰项估计值包括:
使用UDE原理,通过低通滤波器Gf(s)对风电机组数学模型中的不确定及干扰项D(t)进行逼近,记录低通滤波器Gf(s)的冲击响应为gf(t),得到D(t)的不确定及干扰项估计值 的表达式如下:
其中,“*”表示卷积。
6.根据权利要求5所述的一种风电机组桨距角控制方法,其特征在于,所述获得反馈控制信号的最终表达式包括:
将不确定及干扰项估计值代替反馈控制信号的初步表达式中的D(t),得到:
将式①进行拉普拉斯变换可得:
β(s)=f(s)+(sωr(s)-f(s)+β(s))Gf(s)+ke(s) ②
其中β(s)、f(s)、ωr(s)和e(s)分别表示β、ωr和e的拉普拉斯变换,继续整理式②,可得:
将式③进行拉普拉斯反变换可得反馈控制信号的最终表达式如下:
7.根据权利要求6所述的一种风电机组桨距角控制方法,其特征在于,所述Gf(s)选取为一阶低通滤波器τ是滤波器的时间常数。
8.根据权利要求1所述的一种风电机组桨距角控制方法,其特征在于,所述获取前馈控制信号包括:获取一段时间内机组有效风速信息和机组相应输出数据中的桨距角信号,使用当前时刻以及之前连续三个时刻点的有效风速值作为训练特征,将当前时刻的桨距角值作为训练目标训练基于ELM的风电机组动态非线性风速-桨距角关系模型fELMFF,从动态非线性风速-桨距角关系模型fELMFF中获取在线实时前馈控制信号:
其中是当前时刻的有效风速估计值,和是从当前时刻连续向前数三个时刻点的有效风速估计值,其中两个连续时刻点之间的时间长度是控制周期。
9.根据权利要求1所述的一种风电机组桨距角控制方法,其特征在于,还包括获取风电机组桨距角控制信号:
风电机组桨距角控制信号=反馈控制信号+前馈控制信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910830549.6A CN110454329B (zh) | 2019-09-04 | 2019-09-04 | 一种风电机组桨距角控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910830549.6A CN110454329B (zh) | 2019-09-04 | 2019-09-04 | 一种风电机组桨距角控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110454329A true CN110454329A (zh) | 2019-11-15 |
CN110454329B CN110454329B (zh) | 2021-03-30 |
Family
ID=68490521
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910830549.6A Active CN110454329B (zh) | 2019-09-04 | 2019-09-04 | 一种风电机组桨距角控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110454329B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110985286A (zh) * | 2019-12-04 | 2020-04-10 | 浙江大学 | 一种基于elm的新型风电机组桨距角控制方法 |
CN110985288A (zh) * | 2019-12-04 | 2020-04-10 | 浙江大学 | 一种基于ude的风电机组保性能控制方法 |
CN111075647A (zh) * | 2019-12-04 | 2020-04-28 | 浙江大学 | 一种基于elm的变速风电机组最大风能捕获方法 |
CN111219293A (zh) * | 2020-02-18 | 2020-06-02 | 南京理工大学 | 基于线性自抗扰控制的变桨控制器设计方法 |
CN112648140A (zh) * | 2020-12-21 | 2021-04-13 | 北京华能新锐控制技术有限公司 | 基于信号重构的风电机组桨距角编码器故障容错方法 |
CN114320770A (zh) * | 2021-11-11 | 2022-04-12 | 华能新能源股份有限公司 | 风电机组的有效风速测量方法及装置 |
CN115370537A (zh) * | 2022-10-25 | 2022-11-22 | 山东大学 | 基于雷达前馈控制的风电调频载荷优化控制方法及系统 |
CN116378898A (zh) * | 2023-06-05 | 2023-07-04 | 南京信息工程大学 | 一种风电系统桨距角控制方法、装置、系统及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105673325A (zh) * | 2016-01-13 | 2016-06-15 | 湖南世优电气股份有限公司 | 基于rbf神经网络pid的风电机组独立变桨控制方法 |
EP3269974A1 (en) * | 2016-07-12 | 2018-01-17 | Siemens Aktiengesellschaft | Rotor blade pitch arrangement |
CN108167120A (zh) * | 2017-12-11 | 2018-06-15 | 浙江大学 | 一种变速风力发电机组的变桨变矩联合控制方法 |
-
2019
- 2019-09-04 CN CN201910830549.6A patent/CN110454329B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105673325A (zh) * | 2016-01-13 | 2016-06-15 | 湖南世优电气股份有限公司 | 基于rbf神经网络pid的风电机组独立变桨控制方法 |
EP3269974A1 (en) * | 2016-07-12 | 2018-01-17 | Siemens Aktiengesellschaft | Rotor blade pitch arrangement |
CN108167120A (zh) * | 2017-12-11 | 2018-06-15 | 浙江大学 | 一种变速风力发电机组的变桨变矩联合控制方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110985286A (zh) * | 2019-12-04 | 2020-04-10 | 浙江大学 | 一种基于elm的新型风电机组桨距角控制方法 |
CN110985288A (zh) * | 2019-12-04 | 2020-04-10 | 浙江大学 | 一种基于ude的风电机组保性能控制方法 |
CN111075647A (zh) * | 2019-12-04 | 2020-04-28 | 浙江大学 | 一种基于elm的变速风电机组最大风能捕获方法 |
CN111075647B (zh) * | 2019-12-04 | 2021-06-22 | 浙江大学 | 一种基于elm的变速风电机组最大风能捕获方法 |
CN110985288B (zh) * | 2019-12-04 | 2021-07-02 | 浙江大学 | 一种基于ude的风电机组保性能控制方法 |
CN111219293A (zh) * | 2020-02-18 | 2020-06-02 | 南京理工大学 | 基于线性自抗扰控制的变桨控制器设计方法 |
CN111219293B (zh) * | 2020-02-18 | 2021-12-17 | 南京理工大学 | 基于线性自抗扰控制的变桨控制器设计方法 |
CN112648140A (zh) * | 2020-12-21 | 2021-04-13 | 北京华能新锐控制技术有限公司 | 基于信号重构的风电机组桨距角编码器故障容错方法 |
CN114320770A (zh) * | 2021-11-11 | 2022-04-12 | 华能新能源股份有限公司 | 风电机组的有效风速测量方法及装置 |
CN115370537A (zh) * | 2022-10-25 | 2022-11-22 | 山东大学 | 基于雷达前馈控制的风电调频载荷优化控制方法及系统 |
CN116378898A (zh) * | 2023-06-05 | 2023-07-04 | 南京信息工程大学 | 一种风电系统桨距角控制方法、装置、系统及存储介质 |
CN116378898B (zh) * | 2023-06-05 | 2023-09-12 | 南京信息工程大学 | 一种风电系统桨距角控制方法、装置、系统及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110454329B (zh) | 2021-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110454329B (zh) | 一种风电机组桨距角控制方法 | |
CN103850876B (zh) | 一种适用于无载荷测量的风电机组独立变桨控制方法 | |
CN108334672A (zh) | 基于有效风速估计的变速风力发电机组最大风能捕获方法 | |
CN105179164B (zh) | 基于t-s模糊模型的风能转换系统滑模控制方法及装置 | |
CN106979126B (zh) | 基于svr的风力发电机组高风速段有效风速估计方法 | |
CN106774276B (zh) | 风电场自动发电控制系统测试平台 | |
CN102588211A (zh) | 一种风力发电机组全工况模型预测控制方法及系统 | |
CN107045574B (zh) | 基于svr的风力发电机组低风速段有效风速估计方法 | |
CN110345006A (zh) | 一种风力发电机组低风速区最大功率跟踪优化控制方法 | |
CN110500239A (zh) | 一种基于高阶滑模观测器的风电机组自适应反演滑模优化控制方法 | |
CN111425347B (zh) | 基于转矩增益系数优化的风电机组最大功率点跟踪控制方法 | |
CN112784509A (zh) | 基于改进高斯分布尾流模型的风力发电机功率计算方法 | |
CN111502915A (zh) | 基于激光雷达测风的风电机组智能控制系统及方法 | |
CN111188732B (zh) | 一种风力发电变桨鲁棒容错控制方法 | |
CN111075647B (zh) | 一种基于elm的变速风电机组最大风能捕获方法 | |
CN110714881B (zh) | 一种风力机性能预测控制方法及装置 | |
CN115807734A (zh) | 一种基于尾流跟踪的海上风电场场级协同控制策略 | |
Wang et al. | Composite model-free adaptive predictive control for wind power generation based on full wind speed | |
CN111502914A (zh) | 基于线性变参数系统的风机变桨控制器设计方法 | |
CN110263477B (zh) | 一种风力发电机组叶尖速比获取方法 | |
CN115788773A (zh) | 一种漂浮式风电机组多目标变桨控制、系统及设备 | |
CN110985290B (zh) | 一种基于支持向量回归的最优转矩控制方法 | |
CN110985286B (zh) | 一种基于elm的风电机组桨距角控制方法 | |
Mahdizadeh | LIDAR-Assisted Exact Output Regulation for Multi-Megawatt Wind Turbines | |
Shang et al. | A genetic algorithm for optimizing yaw operation control in wind power plants |
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 | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A pitch angle control method for wind turbine unit Effective date of registration: 20211119 Granted publication date: 20210330 Pledgee: Wuhan area branch of Hubei pilot free trade zone of Bank of China Ltd. Pledgor: WINDMAGICS (WUHAN) CO.,LTD. Registration number: Y2021420000126 |