CN109072875B - 使用实时增益计算的风力涡轮机的控制 - Google Patents

使用实时增益计算的风力涡轮机的控制 Download PDF

Info

Publication number
CN109072875B
CN109072875B CN201780026695.2A CN201780026695A CN109072875B CN 109072875 B CN109072875 B CN 109072875B CN 201780026695 A CN201780026695 A CN 201780026695A CN 109072875 B CN109072875 B CN 109072875B
Authority
CN
China
Prior art keywords
sensitivity
wind turbine
pitch
control
blade
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
CN201780026695.2A
Other languages
English (en)
Other versions
CN109072875A (zh
Inventor
J·D·格林内特
T·克吕格尔
A·扎伊博
J·S·汤姆森
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.)
Vestas Wind Systems AS
Original Assignee
Vestas Wind Systems AS
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 Vestas Wind Systems AS filed Critical Vestas Wind Systems AS
Publication of CN109072875A publication Critical patent/CN109072875A/zh
Application granted granted Critical
Publication of CN109072875B publication Critical patent/CN109072875B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors 
    • F03D7/02Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/04Automatic control; Regulation
    • F03D7/042Automatic control; Regulation by means of an electrical or electronic controller
    • F03D7/043Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
    • F03D7/045Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with model-based controls
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors 
    • F03D7/02Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/022Adjusting aerodynamic properties of the blades
    • F03D7/0224Adjusting blade pitch
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors 
    • F03D7/02Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/04Automatic control; Regulation
    • F03D7/042Automatic control; Regulation by means of an electrical or electronic controller
    • F03D7/043Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
    • F03D7/046Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with learning or adaptive control, e.g. self-tuning, fuzzy logic or neural network
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2260/00Function
    • F05B2260/82Forecasts
    • F05B2260/821Parameter estimation or prediction
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2260/00Function
    • F05B2260/84Modelling or simulation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/30Control parameters, e.g. input parameters
    • F05B2270/32Wind speeds
    • 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

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Wind Motors (AREA)

Abstract

一种风力涡轮机控制系统,包括:配置成输出用于风力涡轮机的控制机构的控制信号的至少一个控制模块;以及用于计算与所述控制模块相关联的增益参数的增益计算器,其中所述增益参数基于计算机化的实时叶片模型计算,所述计算机化的实时叶片模型使用确定的风力涡轮机操作点作为输入。所述叶片模型可以是叶片元件动量模型。另一方面,本发明可以被表述为一种控制风力涡轮机的控制机构的方法。

Description

使用实时增益计算的风力涡轮机的控制
技术领域
本发明涉及一种用于控制风力涡轮机的机构或致动系统的方法、一种相关联的控制系统并且还涉及一种配备有这种控制系统的风力涡轮机。
背景技术
现代公用事业规模的风力涡轮机被设计用于在各种风况下操作以便最大化从风的能量提取。为了实现该目的,风力涡轮机通常配备有叶片,叶片的桨距角可以通过叶片桨距控制系统来控制。通过使每个叶片围绕其纵向轴线旋转来控制叶片桨距,该旋转改变叶片相对于迎面而来的风流的迎角。
桨距控制是高于额定风速的主要控制机制,其中风力涡轮机将倾向于满载荷策略操作,该策略的目的是保持恒定的功率输出。这可以通过控制叶片桨距角来实现,以确保由转子产生的转矩不会增加超过可允许的极限,例如由功率参考设定点确定的极限。
叶片桨距控制可以集体地和单独地应用。例如,集体部件是其中所有叶片的桨距角同时改变的桨距调节,并且是用于调节由转子产生的转矩的主要控制动作。相反,单独的桨距部件可以彼此独立地调节所选叶片的桨距并且可以用于平衡施加在转子轴上的载荷。
桨距控制系统的配置,特别是在控制系统中实施的各种控制器的环路增益,很大程度上由风力涡轮机的操作范围确定。通常,在设计阶段期间执行建模以确定确保桨距控制系统能够跨越风力涡轮机的大致上整个预测的操作范围实现其期望的控制目的所需的增益。然而,该方法具有局限性,因为仅针对风力涡轮机的操作范围的选定点模拟了环路增益。这可能降低系统在更不寻常或极端的操作条件期间控制转子载荷的能力。
正是在这种背景下设计了本发明的各实施方式。
发明内容
根据本发明的第一方面,提供了一种风力涡轮机控制系统,其包括:配置成输出用于风力涡轮机的控制机构的控制信号的至少一个控制模块;以及用于计算与所述控制模块相关联的增益参数的增益计算器,其中基于计算机化的实时叶片模型来计算增益参数,所述计算机化的实时叶片模型使用确定的风力涡轮机操作点作为输入。叶片模型可以是叶片元件动量模型。
另一方面,本发明可以被表述为一种控制风力涡轮机的控制机构的方法,该方法包括:确定至少一个增益参数,其中所述增益参数基于来自计算机化的实时叶片模型的输出,所述计算机化的实时叶片模型使用至少一个测量的风力涡轮机操作点作为输入;以及基于所确定的至少一个增益参数确定用于控制机构的控制信号。
本发明还可以被认为存在于:一种风力涡轮机,其包括如上定义的控制系统;一种包括处理器、存储器模块和输入/输出系统的控制器,并且其中所述存储器包括一组程序代码指令,当由处理器执行该程序代码指令时,实施如上定义的方法;以及一种从通信网络可下载和/或存储在机器可读介质上的计算机程序制品,该计算机程序制品包括用于实施如上定义的方法的程序代码指令。
在一个实施例中,增益参数是反馈控制模块的一部分,并且所述增益参数是在反馈控制模块中实施的控制律的控制项。尽管可以计算单个增益参数,但是本发明还可应用于包括与一个或多个控制模块相关联的多个增益模块的实施例。
有利地,本发明提供了一种更灵活和响应更快的控制系统,即使那些风况超出通常预测的范围,该控制系统也能够适应变化的风况。由于控制器增益参数是实时计算的并且基于建模的叶片特征,因此叶片响应于控制变量(例如叶片桨距角)的变化的方式可以被调整为考虑风况的严重性。
所述增益参数或每个增益参数可以基于由计算机化的实时叶片模型计算的相应空气动力学灵敏度值,其中该灵敏度值是力响应于控制变量的变化的灵敏度的指示。在这样的实施例中,灵敏度值是以下组中的一个:
桨距对推力灵敏度,其提供由叶片产生的沿着转子轴线的推力依据叶片桨距角变化的灵敏度的指示;
桨距对转矩灵敏度,其提供转子转矩依据叶片桨距角变化的灵敏度的指示;
桨距对翼载荷灵敏度,其提供叶片翼载荷响应于桨距角变化的灵敏度的指示;以及
桨距对弦向弯矩灵敏度,其提供弦向弯矩响应于桨距角变化的灵敏度的指示。
在以上定义的各实施方式中,计算机化的实时叶片模型可以接收估算的风速值作为输入。估算的转子平面风速值可以是风力涡轮机的旋转叶片的平面扫掠区域的风速。可以依据实际产生的功率值和估算产生的功率值之间的差异来计算估算的风速值。此外,估算产生的功率值可以从计算机化的实时叶片模型获得。
在通常的风力涡轮机控制中,基于由风力涡轮机的转子经历的风域来选择风力涡轮机的期望操作点。操作点被定义为多维参数空间中的点,所述多维参数空间跨越发电机速度、桨距角、电功率、电转矩、风速以及用于控制风力涡轮机的进一步参数中的两个或多个参数。在一实施方式中,基于转子的估算的/经历的风域来选择风力涡轮机的操作点。
在本申请的范围内,明确地意图在前面的段落中、在权利要求书中和/或在以下描述和附图中阐述的各种方面、实施方式、示例和替代方案,并且特别是其各个特征,可以独立地或以任何组合采取。也就是说,任何实施方式的所有实施方式和/或特征可以以任何方式和/或组合进行组合,除非这些特征是不兼容的。申请人保留更改任何最初提交的权利要求或相应地提交任何新权利要求的权利,包括修改任何最初提交的权利要求以从属于和/或包括任何其它权利要求的任何特征的权利,尽管最初并未那样声明。
附图说明
为了更全面地理解本发明,现在将仅通过示例的方式参考以下附图来描述本发明,其中:
图1是具有桨距控制系统的风力涡轮机的视图;
图2是更详细地示出风力涡轮机的桨距控制系统的功能框图;
图3是更详细地示出图2中的桨距控制系统的一部分的功能框图;
图4是更详细地示出图2中的桨距控制系统的另一部分的功能框图;以及
图5是如图3中所示的风力估算器的功能框图。
具体实施方式
首先参考图1,风力涡轮机10包括支撑机舱14的塔架12,转子16安装到该机舱14。转子16包括联接到轮毂20的一组叶片18。该示例中的风力涡轮机10是包括三个叶片的水平轴风力涡轮机(HAWT)。然而,技术人员将意识到替代构造。
如在构造成用于变速操作的风力涡轮机设计中常见的那样,风力涡轮机10配备有桨距控制系统22,叶片18能够通过该桨距控制系统22被控制,使得叶片18围绕其纵向轴线角度上可调节。这可以通过桨距控制器24来实现,桨距控制器24可操作于命令一组桨距致动器26到相应的桨距位置。
叶片桨距控制系统22的实施例在图2中更详细地示出。总体上,叶片桨距控制器24包括多个桨距控制模块30a-d,所述桨距控制模块可集体操作以向相应的桨距致动器32a-c提供三个桨距控制信号。更特别地,第一叶片的桨距由第一桨距致动器32a响应于第一桨距控制信号PA设定,第二桨叶的桨距由第二桨距致动器32b响应于第二桨距控制信号PB设定,并且第三叶片的桨距由第三桨距控制致动器30c响应于第三桨距控制信号PC设定。
三个桨距致动器32a-c在此示出为单个单元,但是应当理解,在实践中各桨距致动器会特定于各叶片18中的每一个。
在该实施方式中,桨距控制模块30a-d中的每一个实施合适的控制律,该控制律被配置成控制过程变量(例如测量的桨距位置)使得所述过程变量满足参考或控制变量(例如目标桨距位置)。这种控制律可以是基于PID的控制律,但是如果情况允许可以是简化的控制律。
在图2中,桨距控制器被示为单独的实体,然而通常上桨距控制器的功能可以实施为风力涡轮机的通用或主要控制器的功能元件。
四个桨距控制模块30a-d一起工作以提供复合的桨距信号输入,所述桨距信号输入被组合以形成独立的桨距控制信号PA-PC。因此,控制模块30a-30d中的每一个提供作为对于每个叶片的最终桨距信号PA-PC的组成部分的相应桨距信号。
总的来说,四个桨距控制模块30a-30d是:推力控制模块30a、速度控制模块30b、独立桨距控制(IPC)模块30c以及侧向塔架阻尼(SSTD)控制模块30d。
推力控制模块30a负责调节在转子上产生的轴向推力。所述推力倾向于使机舱倾斜并且因此在风力涡轮机的塔架上施加弯矩,所述弯矩期望被控制在可接受的限度内。
速度控制模块30b负责将由转子产生的转矩控制为等于要求的或参考转矩值。该参考值可以由风力涡轮机的更高级别的速度/功率控制器设定。
IPC控制模块30c负责控制叶片的独立桨距以便满足每个叶片的预定的翼载荷限制。这种控制信号分量本质上可以是周期的,以便当每个叶片上的载荷随着转子旋转而变化时(例如响应于转子盘中不同高度处的阵风),将翼载荷控制在可接受的限度内。
SSTD控制模块30d负责为机舱的侧向运动(也就是说横向于转子轴线的运动)提供阻尼。SSTD控制模块30d通过提供关于每个叶片的周期桨距调节以便调整作用在转子上的水平力并且从而调整风力涡轮机的支撑结构来完成这一点。
推力控制模块30a和速度控制模块30b的输出信号(在此标记为P1和P2)分别组合成单个集体桨距控制信号PCOLL。然后将该信号分成三个分量并输入相应的求和点34a-c。然后将集体叶片桨距信号PCOLL与由IPC控制模块30c和SSTD控制模块30d输出的独立叶片桨距控制信号组合,如下所述。
IPC控制模块30c和SSTD控制模块30d各自输出三个桨距控制信号,在此分别示出为P3-P5和P6-P8。由于IPC控制模块30c和SSTD控制模块30d输出特定于每个叶片之一的控制信号,所以控制信号对在相应求和点36a-c处组合。更特别地,控制信号P3和P6在求和点36a处组合以得到与桨距控制致动器32a相关联的控制信号P10;控制信号P4和P7在求和点36b处组合以得到与桨距致动器32b相关联的控制信号P11;并且控制信号P5和P8在求和点36c处组合以得到与桨距致动器32c相关联的控制信号P12。
最后,叶片特定的桨距控制信号P10-P12在求和点34a-c处与集体桨距信号PCOLL组合,并且被分别输出为叶片桨距信号PA、PB和PC
在这一点上应该理解的是,尽管在此提供了四个桨距控制模块,但是可替代和稍微简化的系统可以包括更少的控制模块。例如,如果通过控制风力涡轮机的速度可以实现足够的桨距控制,则可能仅需要一个控制模块。然而,实际上,桨距控制方法通常以集体桨距控制元件和独立桨距控制元件为特征,特别是在商业规模的风力涡轮机系统中。
控制模块30a-30d中的每一个可以包括至少一个控制器增益值。在已知的控制方法中,控制器增益值会在风力涡轮机的安装之前或安装期间被优化。可替代地,已知基于预设调度依据风力涡轮机操作所处的操作条件(风速、发电机功率、转子速度)改变控制器增益值。这样的增益调度在离线环境中确定并且被配置成一旦达到风力涡轮机的预定操作点就更新一个或多个控制模块的增益值。
在本发明的所示实施方式中,控制器增益值依据计算机化的在线或“实时”叶片元件动量(BEM)模型50可变化。BEM模型实时运行并且向叶片控制模块30a-d输出各种空气动力学灵敏度参数。反过来,叶片控制模块30a-d可操作于计算在其中的控制律中实施的更新的增益值。
总的来说,BEM模型50被实施为从适合的来源接收操作数据输入的算法,如将描述的,并且该BEM模型50可操作用于输出:
桨距对推力灵敏度信号52,其提供由叶片产生的(沿着转子轴线的)推力依据桨距角变化的灵敏度的指示;
桨距对转矩灵敏度信号54,其提供转子转矩依据桨距角变化的灵敏度的指示;
桨距对翼载荷弯矩灵敏度信号56,其提供叶片翼载荷响应于桨距角变化的灵敏度的指示;以及
桨距对弦向弯矩灵敏度信号58,其提供弦向弯矩响应于桨距角变化的灵敏度的指示。
为了导出上述空气动力学灵敏度信号52-58,BEM模型50接收以下输入:发电机功率信号60、转子速度信号62、测量的叶片桨距角信号64以及空气密度信号66。这些信号可以从与风力涡轮机相关联的适当传感系统获得,并且例如被直接传输到控制器或者源自系统数据总线。
现在还将参考图3,其示出了BEM模型50的可能的结构实施例。
BEM模块50包括实施叶片元件动量模型的空气动力学灵敏度计算模块70,以及风力估算模块72。如技术人员所知,BEM模型是叶片元件理论和动量理论的融合,其用于分析风力涡轮机转子的性能。在此将不提供完整的描述,因为这些理论会被风力涡轮机叶片设计领域中工作的工程师很好地理解,并且此外在当代教科书中记录,例如,《风力涡轮机的空气动力学》,作者M Hansen:ISBN no.978-1-84407-438-9,以及《风能手册》,作者TBurton:ISBN no.978-0-470-69975-1。
空气动力学灵敏度计算模块70基于叶片18的升力和阻力系数实施叶片载荷计算,该升力和阻力系数是基于叶片18设计的预定值。该计算基于风力涡轮机的操作点,所述操作点可以由转子速度信号62、测量的叶片桨距角信号64以及空气密度信号66定义。除了这些信号之外,空气动力学灵敏度计算模块70还接收指示出估算的通过转子区域的风速的信号,该信号由风速估算模块72确定。
风速估算模块72提供由风力涡轮机的叶片18限定的转子平面处或附近的估算的风速值。该值不可通过通常安装到风力涡轮机的机舱的标准风速传感器来测量。重要的是注意到,通过常规风速测量和估算技术(例如风速计和LIDAR传感器(光照方向和测距))测量的“自由风速”不是在此讨论的在线叶片模型计算中使用的合适的参数。以下讨论将使这点清楚。
为了能够计算沿着叶片产生的力,有必要知道相对风速(Vrel)和叶片的迎角。如技术人员会理解的,相对风速取决于两个主要参数:转子平面处的风速(Va)和由转子的旋转引起的风速(Vrot)。如本领域技术人员所理解的,动量理论指出自由风速(Vo)和转子平面风速(Va)通过转子的轴向感应(a)相关。所述轴向感应取决于转子平面风速(Va)、叶片的桨距角以及转子速度。更特别地,该表达式可以被认为是Va=(1-a)Vo。基于上述理论对转子平面风速的计算是计算密集型的并且不适合关键机器控制环境中的实时实施。
然而,在该实施方式中,风速估算模块72基于其作为来自空气动力学灵敏度计算模块70的输入接收的参数即预测的空气动力功率(PBEM)和作为直接输入接收的实际产生的功率(PGEN)的比较来计算估算的转子平面风速(VEST)。因此,避免了进行计算密集型实时计算的需要,同时还实现了转子的平面处风速(也就是说在转子旋转时通过由转子的叶片限定的平面扫琼区域的风速)的准确估算。
风速估算模块72的一个实施例在图5中示出。在该实施例中,基于风力涡轮机产生的实际产生的功率与如通过叶片载荷计算器模块估算的估算产生的功率之间的差异来计算转子盘处的风速的估算值。
风速估算模块72的一个实施例在图5中示出。在该实施例中,基于风力涡轮机的实际产生的功率与如通过空气动力学灵敏度计算模块70估算的估算产生的空气动力学功率之间的差异来计算转子盘处的风速的估算值。如从以下讨论中将理解的,这种计算转子风速的方法比仅基于涉及转子的轴向引入因子的实时迭代计算的动量理论的估算方法更不计算密集。
在图5中,通过从测量的电网功率信号PGEN中减去估算的电网功率信号PGEST在求和点73处确定功率误差信号PERR。然而,信号PGEST不是直接输入,而是基于空气动力学功率信号PBEM计算的,如由空气动力学灵敏度计算模块70所提供的。更特别地,空气动力学功率信号PBEM被馈送到损耗计算器74,该损耗计算器计算与风力涡轮机传动系相关联的力学损耗和与风力涡轮机转换器设备相关联的电损耗,由此提供功率损耗信号PLOSS。然后,在求和点75处从空气动力学功率信号PBEM中减去功率损耗信号PLOSS以得到估算的电网功率信号PGEST
基于功率误差信号PERR,通过使功率误差信号PERR通过增益模块75确定了校正信号PCORR。校正信号PCORR从先前样本被添加到转子平面风力估算信号VEST(通过单元延迟76馈送)并且将得到的信号VEST_0供应到稳定性限制模块77。稳定性限制模块77的功能实际上是验证对到来的VEST_0信号做出的校正。为此,稳定性限制模块77还接收桨距角信号64和转子速度信号62作为输入。使用这两个信号,稳定性限制模块77确定转子风速的有效上限和下限并确保到来的转子风速信号VEST_0在这些限制内。稳定性限制模块77可以以技术人员理解的各种方式实施该功能。作为示例,一个实施例可以基于查找表,该查找表基于桨距角和转子速度将可接受的上转子速度值和下转子速度值与特定操作点相关联。
从上面的讨论应当理解,风力估算模块72不需要基于叶片模型的计算量重的计算,而是将估算的转子风速基于由风力涡轮机传递的测量的功率与由在空气动力学灵敏度计算模块70中实施的叶片模型确定的估算的功率之间的收敛。
返回图3,由空气动力学灵敏度计算模块70接收的输入信号转子速度信号62、桨距角信号64、空气密度信号66和转子平面风力估算信号VEST向其提供与风力涡轮机的当前操作点有关的信息;即当前叶片桨距、转子风速、转子速度和空气密度。基于当前操作点的该定义,空气动力学灵敏度计算模块70能够计算各种空气动力学灵敏度。
例如,转子推力对桨距角变化的灵敏度可以通过以下等式确定(例如,0阶估算):
Figure GDA0003049222120000091
以上推导涉及如上所述的桨距对推力灵敏度信号52。
进一步的示例是叶片翼载荷弯矩对桨距角变化的灵敏度,其可以通过以下等式检测:
Figure GDA0003049222120000092
以上推导涉及如上所述的桨距对翼载荷弯矩灵敏度信号56。
此外,下面的等式提供了转子转矩对桨距角变化的灵敏度,其涉及如上所述的桨距对转矩灵敏度信号54:
Figure GDA0003049222120000093
而下面的等式提供了弦向方向或“平面中”弯矩对桨距角变化的灵敏度,其对应于如上所述的桨距对弦向弯矩灵敏度信号58。
Figure GDA0003049222120000094
在以上各等式中,θ0、V0、Ω0(桨距角、风速、转子速度)是计算的操作点,ΔV、Δθ是完成0阶近似的风速和桨距范围,并且Mx、Fy分别是相应的估算的叶片根部平面中弯矩和平面外力。总体上,将理解BEM模型50能够确定以下灵敏度测量:
·推力Fy对桨距、转子速度和风速的灵敏度。
·驱动弯矩My对桨距、转子速度和风速的灵敏度。
·叶片弯矩Mx对桨距、转子速度和风速的灵敏度。
·横向力Fx对桨距、转子速度和风速的灵敏度。
基于这些基础灵敏度可以计算其它灵敏度,例如:
·倾斜和横摆弯矩对周期桨距、风速和转子速度的灵敏度。
·功率和转矩对集体桨距、风速和转子速度的灵敏度。
·侧向力和上下力对周期桨距、风速和转子速度的灵敏度。
·转子推力对集体桨距、风速和转子速度的灵敏度。
从BEM模型50输出的空气动力学灵敏度变量被输入到控制模块52-58并且此后用于实时计算在相应的控制模块30a-d内使用的更新的增益值。以下关于可以如何更新合适的控制增益的示例将提供对该实施例的充分理解。
推力控制模块30a在图4中更详细地示出。可以看出,推力控制模块30a包括控制律80、增益块82和增益计算器84。如已经提到的,控制律80可以适当的是PID控制律或其简化版的形式。
控制律80作用于修正的误差输入E'并输出集体桨距信号PCOLL_THRUST,以便将修正的误差信号减小到零。修正的误差信号E'由增益块82确定,该增益块82作用于初始误差信号E,该初始误差信号E是从目标推力值86与在求和点90处计算的实际或测量的推力值88之间的差异导出的。注意到,目标推力值86可以从更高级别的控制器接收,并且测量的推力值88可以从风力涡轮机的数据总线接收或者直接从合适的测量或控制系统接收。
误差信号E被输入到增益块82,然后该增益块82计算修正的误差信号E'。如上所述,推力控制模块30a从实时BEM模型50接收连续更新的空气动力学灵敏度信号,更具体地说是桨距对推力灵敏度信号52。这基于风力涡轮机10的主要操作点为推力控制模块30a提供了同步值,该同步值具有由桨距变化所得的推力变化的灵敏度或响应性。作为示例,为了计算合适的增益值,增益计算器84可以使用以下公式:
Figure GDA0003049222120000101
在以上公式中:
·
Figure GDA0003049222120000102
是由在线BEM模型50计算的灵敏度。
·
Figure GDA0003049222120000103
是表示“标称”灵敏度的参数,以及
·KTL是标称灵敏度下的期望控制器增益。
在此应该注意的是,标称灵敏度可以是在转子推力对额定风速下桨距变化的灵敏度,这将是推力最高的位置。因此,在该操作点处可以计算控制器增益KTL,并且增益调度将针对其它操作点进行校正。
将要理解,与用于图4中所示的推力控制模块30a相同的原理和相同的结构也可以应用于其它控制模块30b、30c和30d。
例如,侧向塔架阻尼控制模块30d可以具有共同的结构,尽管可以利用由增益计算器84计算的不同增益值来更新增益块。例如,增益计算器可以实施以下公式来计算实时所需的增益:
Figure GDA0003049222120000111
在以上公式中:
·KSS是标称增益。
·
Figure GDA0003049222120000112
是标称侧向力对周期桨距振幅的灵敏度,以及
·
Figure GDA0003049222120000113
是如从在线BEM推导的平面中的力对桨距的灵敏度。
从上面的讨论将要理解,本发明的各实施方式提供了一种控制系统,其中至少部分地基于风力涡轮机的主要操作点来连续更新该控制系统的选定增益项,如通过计算机化的实时叶片模型所实现的。这提供了对控制下的控制机构(在这种情况下是桨距控制机构)的更准确的调节,因为相关控制器的增益项能够响应于风力涡轮机的操作点的变化而被调节。特别的优点是即使风力涡轮机的操作点超出通常预期的操作点,控制器性能仍然很高,使得风力涡轮机在异常条件(例如极端温度操作、在强风条件期间启动和关闭等等)下表现良好。
技术人员将理解,在不脱离由权利要求限定的本发明构思的情况下,可以对上面讨论的各特定实施方式做出修改。
例如,在示出的实施方式中,增益块82作用于误差信号E,并且然后将修改的误差信号传递给控制律80。然而,应该注意,控制模块30a可以被配置成使得增益计算器84更新包含在控制律内的增益项,例如任何比例、积分或微分增益项。
在示出的各实施方式中,在系统中实施的计算机化的叶片模型是叶片元件动量模型。可以使用具有类似目的的其它叶片模型,例如基于计算流体动力学理论的模型,尽管这样的模型可能不是如此计算高效的。
注意,在上面的讨论中,各种控制模块和功能块被描述为分离的功能单元。然而,这将被理解为公认的惯例,并且不应认为这种结构对如何在软件、固件或硬件中实施该功能施加任何限制。

Claims (12)

1.一种风力涡轮机控制系统,包括:
配置成输出用于风力涡轮机的控制机构的控制信号的至少一个控制模块,
用于计算与所述控制模块相关联的增益参数的增益计算器,以及
使用基于由风力涡轮机的转子经历的风域选择的风力涡轮机操作点作为输入的计算机化实时叶片模型,所述计算机化实时叶片模型包括:
风速估算模块,其配置成确定估算的转子平面风速值,以及
空气动力学灵敏度计算器模块,其设置成基于所述估算的转子平面风速值和已确定的所述风力涡轮机操作点来确定空气动力学灵敏度值;
其中基于所述空气动力学灵敏度值计算所述增益参数,
其中所述风力涡轮机操作点由转子速度信号、测量的叶片桨距角信号以及空气密度信号定义,
其中所述控制机构是桨距致动器,
其中所述估算的转子平面风速值是旋转叶片的平面扫掠区域的风速,
其中所述空气动力学灵敏度值是以下组中的一个:
桨距对推力灵敏度,其提供由叶片产生的沿着转子轴线的推力依据叶片桨距角变化的灵敏度的指示;
桨距对转矩灵敏度,其提供转子转矩依据叶片桨距角变化的灵敏度的指示;
桨距对翼载荷灵敏度,其提供叶片翼载荷响应于桨距角变化的灵敏度的指示;以及
桨距对弦向弯矩灵敏度,其提供所述弦向弯矩响应于桨距角变化的灵敏度的指示。
2.根据权利要求1所述的风力涡轮机控制系统,其中所述增益参数是反馈控制模块的一部分。
3.根据权利要求2所述的风力涡轮机控制系统,其中所述增益参数是在所述反馈控制模块中实施的控制律的控制项。
4.根据权利要求2或3所述的风力涡轮机控制系统,其中确定了多个增益参数,并且其中所述多个增益参数中的一个或多个与一个或多个控制模块中的每一个相关联。
5.根据权利要求4所述的风力涡轮机控制系统,其包括多个控制模块,其中每个控制模块提供所述控制信号的分量。
6.根据权利要求1至3中任一项所述的风力涡轮机控制系统,其中所述增益参数或每个增益参数基于由所述计算机化实时叶片模型计算的相应空气动力学灵敏度值,其中所述空气动力学灵敏度值指示出力响应于控制变量变化的灵敏度。
7.根据权利要求1至3中任一项所述的风力涡轮机控制系统,其中所述计算机化实时叶片模型是叶片元件动量模型。
8.根据权利要求1至3中任一项所述的风力涡轮机控制系统,其中依据实际产生的功率值和估算产生的功率值之间的差异来计算所述估算的转子平面风速值。
9.根据权利要求8所述的风力涡轮机控制系统,其中所述估算产生的功率值是从所述计算机化实时叶片模型获得的。
10.一种风力涡轮机,其包括根据权利要求1至9中任一项所述的控制系统。
11.一种控制风力涡轮机的控制机构的方法,所述方法包括:
使用基于由风力涡轮机的转子经历的风域选择的至少一个风力涡轮机操作点作为输入来实施计算机化实时叶片模型,以确定估算的转子平面风速值,并基于所述估算的转子平面风速值以及确定的风力涡轮机操作点来确定空气动力学灵敏度值;
基于所述空气动力学灵敏度值来确定与所述风力涡轮机相关联的至少一个增益参数;以及
基于确定的所述至少一个增益参数确定用于控制机构的控制信号,
其中风力涡轮机操作点由转子速度信号、测量的叶片桨距角信号以及空气密度信号定义,
其中控制机构是桨距致动器,
其中所述估算的转子平面风速值是旋转叶片的平面扫掠区域的风速,
其中所述空气动力学灵敏度值是以下组中的一个:
桨距对推力灵敏度,其提供由叶片产生的沿着转子轴线的推力依据叶片桨距角变化的灵敏度的指示;
桨距对转矩灵敏度,其提供转子转矩依据叶片桨距角变化的灵敏度的指示;
桨距对翼载荷灵敏度,其提供叶片翼载荷响应于桨距角变化的灵敏度的指示;以及
桨距对弦向弯矩灵敏度,其提供所述弦向弯矩响应于桨距角变化的灵敏度的指示。
12.一种机器可读介质,包括存储在其上的计算机程序,其中所述计算机程序包括用于实施根据权利要求11所述的方法的程序代码指令。
CN201780026695.2A 2016-03-30 2017-03-30 使用实时增益计算的风力涡轮机的控制 Active CN109072875B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
DKPA201670175 2016-03-30
DKPA201670175 2016-03-30
PCT/DK2017/050093 WO2017167343A1 (en) 2016-03-30 2017-03-30 Control of a wind turbine using real-time gain calculation

Publications (2)

Publication Number Publication Date
CN109072875A CN109072875A (zh) 2018-12-21
CN109072875B true CN109072875B (zh) 2021-11-30

Family

ID=59962633

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201780026695.2A Active CN109072875B (zh) 2016-03-30 2017-03-30 使用实时增益计算的风力涡轮机的控制

Country Status (5)

Country Link
US (1) US10907612B2 (zh)
EP (1) EP3436693B1 (zh)
CN (1) CN109072875B (zh)
ES (1) ES2870074T3 (zh)
WO (1) WO2017167343A1 (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2870074T3 (es) 2016-03-30 2021-10-26 Vestas Wind Sys As Control de una turbina eólica usando cálculo de ganancia en tiempo real
EP3667071A1 (en) * 2018-12-13 2020-06-17 Siemens Gamesa Renewable Energy A/S Wind turbine control method and system
CN109779836B (zh) * 2018-12-20 2020-02-21 明阳智慧能源集团股份公司 一种基于遗传算法寻优的风力发电机组发电量提升方法
EP3744973B1 (en) 2019-05-31 2023-10-04 Vestas Wind Systems A/S Controlling flap loading on a wind turbine blade based on predicted flap loading
US11629694B2 (en) * 2019-10-22 2023-04-18 General Electric Company Wind turbine model based control and estimation with accurate online models
WO2021143991A1 (en) 2020-01-16 2021-07-22 Vestas Wind Systems A/S Pitch control of wind turbine blades in a standby mode
CN111188732B (zh) * 2020-01-17 2022-05-13 湖南工业大学 一种风力发电变桨鲁棒容错控制方法
ES2949649T3 (es) 2020-03-02 2023-10-02 Vestas Wind Sys As Sistema de monitorización de control de paso de pala para una turbina eólica
EP4155539A1 (en) * 2021-09-28 2023-03-29 LM Wind Power A/S System and method for estimating energy production from a wind turbine
CN114934874A (zh) * 2022-05-31 2022-08-23 新疆金风科技股份有限公司 风电机组的控制方法和控制设备以及风电机组
US11802545B1 (en) * 2022-09-26 2023-10-31 General Electric Company Method and system for detection and mitigation of edge-wise vibrations in wind turbine blades
WO2024099527A1 (en) * 2022-11-09 2024-05-16 Vestas Wind Systems A/S Determining wind speed at a wind turbine

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4161658A (en) 1978-06-15 1979-07-17 United Technologies Corporation Wind turbine generator having integrator tracking
US4160170A (en) 1978-06-15 1979-07-03 United Technologies Corporation Wind turbine generator pitch control system
WO2005026537A1 (ja) 2003-09-10 2005-03-24 Mitsubishi Heavy Industries, Ltd. ブレードピッチ角度制御装置及び風力発電装置
US7342323B2 (en) * 2005-09-30 2008-03-11 General Electric Company System and method for upwind speed based control of a wind turbine
US8174136B2 (en) 2006-04-26 2012-05-08 Alliance For Sustainable Energy, Llc Adaptive pitch control for variable speed wind turbines
WO2010069323A2 (en) * 2008-12-15 2010-06-24 Vestas Wind Systems A/S Pitch control of a wind turbine
CN101660489B (zh) * 2009-09-23 2013-01-09 大唐南京自动化有限公司 一种兆瓦级风力发电机组联合控制策略
BRPI1000022A2 (pt) * 2010-02-08 2017-01-17 Mitsubishi Heavy Ind Ltd gerador de turbina eólica e método para controlar o ângulo de inclinação da lâmina do mesmo.
US20110229300A1 (en) 2010-03-16 2011-09-22 Stichting Energieonderzoek Centrum Nederland Apparatus and method for individual pitch control in wind turbines
ES2401857B1 (es) * 2011-01-31 2014-03-10 Gamesa Innovation & Technology S.L. Métodos y sistemas de control de aerogeneradores mejorados.
US20130259682A1 (en) * 2012-03-27 2013-10-03 General Electric Company Method of rotor-stall prevention in wind turbines
US8704393B2 (en) 2012-08-09 2014-04-22 General Electric Company System and method for controlling speed and torque of a wind turbine during post-rated wind speed conditions
CN102797631B (zh) * 2012-08-24 2014-05-14 国电联合动力技术有限公司 一种风电机组的最优增益在线自校正方法、系统及其装置
FR2996266B1 (fr) 2012-10-01 2014-09-12 IFP Energies Nouvelles Procede pour controler une eolienne au moyen d'une estimation de la vitesse du vent incident
DK2784303T3 (en) * 2013-03-27 2017-01-30 Alstom Renovables Espana Sl Method of operating a wind turbine
CN104214045B (zh) * 2013-05-30 2017-03-08 成都阜特科技股份有限公司 双馈式变速变桨风力发电机组的独立变桨距控制方法
CN105980703B (zh) * 2013-12-09 2019-03-01 维斯塔斯风力系统集团公司 用于风力涡轮机的操作方法
CN103698399B (zh) * 2013-12-19 2016-06-29 国家电网公司 一种汽轮机叶片枞树型根部超声成像检测方法
CN104533714B (zh) * 2014-11-03 2017-03-01 华北电力大学 一种风力发电机组控制器参数辨识与整定方法
ES2870074T3 (es) 2016-03-30 2021-10-26 Vestas Wind Sys As Control de una turbina eólica usando cálculo de ganancia en tiempo real

Also Published As

Publication number Publication date
CN109072875A (zh) 2018-12-21
EP3436693B1 (en) 2021-05-12
US20190113022A1 (en) 2019-04-18
WO2017167343A1 (en) 2017-10-05
ES2870074T3 (es) 2021-10-26
US10907612B2 (en) 2021-02-02
EP3436693A1 (en) 2019-02-06

Similar Documents

Publication Publication Date Title
CN109072875B (zh) 使用实时增益计算的风力涡轮机的控制
CN109154276B (zh) 使用实时叶片模型的风力涡轮机的控制
CN107532568B (zh) 使用中央控制器计算局部控制目标的多转子风力涡轮机的控制
CN106460793B (zh) 响应于风切变而对风力涡轮机的控制
EP2859223B1 (en) A wind turbine with a load controller
Bossanyi et al. Advanced controller research for multi‐MW wind turbines in the UPWIND project
CN107709766B (zh) 校准风力涡轮机的负载传感器的方法
CN203939626U (zh) 用于在偏航误差过程中减小风力涡轮机上负载的控制系统
US20130302161A1 (en) Controller of wind turbine and wind turbine
JP6494514B2 (ja) 入力風速の予測値を使用した風力タービン制御方法
TWI708893B (zh) 風力發電系統
Andersen et al. Load alleviation on wind turbine blades using variable airfoil geometry
EP4112919B1 (en) Control of a multi-rotor wind turbine system using local mpc controllers
EP3899250B1 (en) Control of side-side and fore-aft vibrational movement of a wind turbine
US10451034B2 (en) Methods and systems for alleviating loads in off-shore wind turbines
US20200072192A1 (en) Method for controlling a wind turbine
EP2499358A2 (en) Improved control of wind turbine blade lift regulating means
EP3987171B1 (en) Controlling noise emissions from individual blades of a wind turbine
EP4166782A1 (en) Wind turbine model predictive control (mpc) with thrust and/or power correction
WO2022179801A1 (en) Method of configuring a wind turbine pitch controller
Bergami et al. Preliminary Evaluation with Feed-Forward Cyclic Control

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