CN101918710B - 桨距和负载缺陷的诊断 - Google Patents

桨距和负载缺陷的诊断 Download PDF

Info

Publication number
CN101918710B
CN101918710B CN2008801189827A CN200880118982A CN101918710B CN 101918710 B CN101918710 B CN 101918710B CN 2008801189827 A CN2008801189827 A CN 2008801189827A CN 200880118982 A CN200880118982 A CN 200880118982A CN 101918710 B CN101918710 B CN 101918710B
Authority
CN
China
Prior art keywords
pitch
load
fault
blade
angle
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
CN2008801189827A
Other languages
English (en)
Other versions
CN101918710A (zh
Inventor
P·布拉特
S·达尔斯高
M·布兰克
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 CN101918710A publication Critical patent/CN101918710A/zh
Application granted granted Critical
Publication of CN101918710B publication Critical patent/CN101918710B/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/0204Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor for orientation in relation to wind direction
    • 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
    • F03D17/00Monitoring or testing of wind motors, e.g. diagnostics
    • 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
    • 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/83Testing, e.g. methods, components or tools therefor
    • 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/10Purpose of the control system
    • F05B2270/103Purpose of the control system to affect the output of the engine
    • 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/328Blade pitch angle
    • 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/331Mechanical loads
    • 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/40Type of control system
    • F05B2270/404Type of control system active, predictive, or anticipative
    • 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)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Wind Motors (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明涉及一种用于例如风轮机桨距和/或负载缺陷的诊断的方法、系统以及计算机可读代码,以及涉及一种使用所述诊断方法和/或包括所述诊断系统的风轮机。

Description

桨距和负载缺陷的诊断
技术领域
本发明涉及用于例如风轮机桨距和/或负载缺陷的诊断的方法、系统以及计算机可读代码,以及涉及一种使用所述诊断方法和/或包括所述诊断系统的风轮机。
背景技术
近年来,风能成为一种相当吸引人的可再生能量来源,并且全世界范围内已经建成了许多风轮机。这增加了对于风轮机部件进行更好维修和监测的需要。当前,部件维修有三种方案:失效维修;基于平均部件寿命的预防性维修;以及基于部件健康监测的预报性维修。然而,关于预报性维修策略,存在一些部件故障(例如涉及到叶片和相关传感器)难以通过已知方法和技术诊断出来,这些故障会导致恶化的效率或者使叶片过载以及缩短其寿命。这种故障不一定导致快速的失效,但是通常在一段时间之后才被注意到,例如当风轮机的平均年发电量与预期发电量不同时。
此外,一些缺陷只能通过现有方法而被检测到,但是如果没有人为干预的话,引起缺陷的部件不能够被隔离。
EP0995904公开了一种WT系统,该系统具有转子,转子具有可调几何攻角的至少两个叶片(B)和传感器(KB,KR,KT),传感器提供测量参数,它给出作用在该系统结构的元件上的当前负载的测量值。叶片角根据测量参数调节,这代表了结构元件的加速度或变形并且给出了力或转矩的测量值。传感器被安装在转子叶片上和/或转子叶片可单独地被调节。
由此需要一种新的方法和系统用于叶片以及相关部件的诊断。特别是用于风轮机叶片以及相关传感器的情况。
关于本发明,在此给出了该专利申请中大多数使用术语的定义。然而,这些术语定义不应限于本发明,并且能够使用任意其它适当的和/或适合的术语定义。
桨距角θ被定义为叶片叶根弦与转子旋转平面之间的角度。桨距旋转是指叶片沿着其纵向轴线的旋转。
方位角δ通过转子的角度位置被定义为叶片(例如第一叶片,叶片No.1,B1)当前角度位置与垂直向下方向之间的顺时针或者逆时针角度(从前面看)。
叶片负载是指叶片由于作用在叶片整体或者部分上的不同的力而承受的转矩。
桨距基准角θref是指桨距角θ的期望值以及方位角基准角δref是指方位角δ的期望值。
发明内容
这里描述的方法是一种预测性维修策略,但是由于故障种类以及本发明,可以在没有人力干预的情况下进行故障处理。
本发明允许叶片和/或传感器的桨距和/或负载故障以足够准确且及时的方式被诊断,从而防止故障发展成失效。
具有被控叶片角(桨距)的风轮机具有用于叶片角的传感器,该传感器被称为桨距传感器。如果该传感器具有故障,那么相关风轮机叶片的桨距将会是不正确的。桨距传感器故障(例如桨距增益故障或者桨距偏离)将会导致上述问题和缺点。
负载传感器可被用于风轮机中,用于测量单个叶片上的负载,如公知的那样风轮机叶片上的负载对于叶片寿命以及对于风轮机效率而言是决定性的。负载传感器故障(例如负载增益故障或者负载偏离)将同样由于使得测量结果几乎无用而导致问题和缺点。
本发明利用负载和桨距测量结果。
可以看到,本发明一个目的是提供一种方法和系统,用于根据桨距和负载测量结果来进行可能故障的诊断,以及提供一种计算机可读或可用介质上的计算机可读代码,该代码适合于执行所述方法。
可以看到,本发明的另一个的是提供一种风轮机,其中使用了所述方法和/或系统。
优选地,本发明以单独的方式或者以任意组合的方式减轻、缓和或者消除了上面提到的一种或多种缺点或者其它缺点。
特别地,可以看到,本发明的目的是提供一种对所述缺点的解决方案。
在本发明的一个方面,提供了一种用于叶片和/或相关传感器中至少一个可能故障的诊断的方法,该方法包括以下步骤:
-在一组方位角下和/或在一组桨距基准角下进行负载和桨距的测量;
-推导负载测量值与桨距角之间的关系;以及
-将所述负载-桨距或者桨距-负载关系与预期的负载-桨距或者桨距-负载关系进行比较,从而确定叶片和/或相关传感器是否出现至少一个可能故障。
在本发明的另一个方面,提供了一种用于诊断与叶片和/或相关传感器中桨距和/或负载相联系的至少一个可能故障的方法,该方法包括以下步骤:
a)通过使用下面两个程序中的任意一个或者两者进行负载和桨距测量:
i)运行桨距循环并且测量负载,以及随后在一个或多个方位角下重复桨距循环测量步骤,和/或
ii)运行方位角循环并且测量负载,以及随后在一个或多个桨距基准角重复方位角循环测量步骤;
b)推导负载测量值与桨距角之间的关系;以及
c)将所述负载-桨距或者桨距-负载关系与预期的负载-桨距或者桨距-负载关系进行比较,从而确定出叶片和/或叶片传感器是否出现至少一个可能故障。
根据本发明的诊断程序由多种活动组成,如故障的检测和隔离,以及可能地还有故障大小的估计。
本发明的主要特征在独立权利要求中描述。
发明的附加特征在从属权利要求中提出。
一般地,在撰写出通过本发明“是有利的”以及涉及到优点时,必须理解的是,该优点可以被看作是由本发明提供的可能优点,并且它还可以被理解为对于获得所述优点而言,本发明以特别的方式、但不是以排他的方式是有利的。
基本上,本发明的各个方面和优点可以在本发明范围内以任意可能方式组合及结合。
通过参考下文中描述的实施例,本发明的这些及其它方面、特征和/或优点将会明显并且得以阐明。
附图说明
参考附图通过示例对本发明实施例进行描述,其中:
图1a和1b显示了具有待诊断的叶片和传感器的风轮机以及一个叶片的横截面,还显示了桨距和方位角;
图2显示了根据本发明在无风静止条件下的叶片负载结构;
图3a、3b、3c和3d显示了当系统经历桨距偏离(3a)、负载偏离(3b)、负载增益故障(3c)或者桨距增益故障(3d)时测量出的桨距和负载如何相对于预期的桨距和负载进行变化;
图4a、4b和4c显示了通过主动隔离算法的故障估计。
具体实施方式
本发明提出了一种方法和系统,用于对叶片桨距和负载测量中增益和/或偏离故障的诊断。本发明不仅允许检测到传感器中的故障,而且能够检测到在安装叶片时的缺陷。它仅仅使用拍打方向负载测量以及使用桨距(叶片角)和方位角(转子位置)基准值的预定序列来诊断故障。故障可以被诊断出来,甚至是在桨距和负载同时出现增益和偏离故障时。主动诊断是在低风速或者零风速下进行的,此时轮机没有生产风能。数据序列的分析没有必要在时域内进行,但是可以作为桨距和方位角的函数,这样能够根据主导条件而在不同时刻完成序列收集数据。
风轮机叶片上的桨距测量偏离将会导致稳定叶片位置的控制器处于错误的桨距。这改变了特定叶片的空气动力学性能并且将会导致下降的年发电量以及可以预见到的负载行为的不同偏差。该测量故障的特征在于,故障并不严重并且在直观上难以诊断出来。如果重新配置闭环,那么不会经历性能下降。
为了能够实现全部的重新配置,必须使用轮机测量值而得到精确的桨距及负载故障估计。故障应当被诊断出来,使得不需要任何人工检查或者人工测量。
图1a显示了根据本发明的风轮机WT,其包括(应当注意到一些部件没有显示出来):
-转子,
-至少一个风轮机叶片B1,B2,B3,其中所述至少一个风轮机叶片B1,B2,B3从转子开始延伸,从而在受到风力时和/或转子通过具有动力的转子旋转装置而旋转时使转子旋转,以呈现出不同方位角δ,并且所述至少一个风轮机叶片B1,B2,B3的至少一部分设置有叶片旋转装置,用于使所述至少一个风轮机叶片B1,B2,B3的至少一部分旋转,从而呈现出不同桨距角θ,
-负载测量装置,用于测量由于重力和/或由于所述至少一个风轮机叶片B1,B2,B3受到的风而导致的所述至少一个风轮机叶片B1,B2,B3上的负载,
-桨距测量装置,用于测量所述至少一个风轮机叶片的至少一部分的桨距角θ,
-用于控制所述风轮机叶片的所述至少一部分的桨距角θ的调整的装置。
可以使用用于故障大小估计的计算装置和/或用于故障调节的计算装置。
此外,在另一个可能实施例中,转子可以被调整,从而呈现出一个或多个方位角δ。方位角可被用于借助测量的计算方法或者其它方法。
在另一个可能实施例中,转子不能够调整,以及方位角可被用于借助测量的计算方法或者其它方法。
在另一个可能的实施例中,转子旋转装置可包括轴承,其中旋转由风力提供并且旋转角度可通过制动器而停止,和/或可以包括电力或者液压马达和/或两者之间的传动装置。
在另一个可能的实施例中,叶片旋转装置可包括设置在转子与风轮机叶片之间的轴承,和/或可以包括电力或者液压马达和/或两者之间的传动装置。
负载测量装置可以设置在风轮机叶片中或者设置在转子中风轮机叶片与转子相连接的位置上。
在另一个可能实施例中,负载测量装置可以是应变仪和/或可以通过发射以及接收的光所提供,并且由此风轮机叶片的偏转和叶片上的负载可以被计算出来。然而,可以使用任意其它适当的负载测量装置或类型,并且本发明不应局限于上面提到的。
桨距角测量装置可通过接近传感器、角度传感装置或者任意其它适当的测量装置所提供。然而,可以使用任意其它适当的桨距角测量装置或类型,并且本发明不应局限于上面提到的。
诊断系统可包括至少一个处理器,用于各种计算等等,例如与用于所述至少一个可能故障的诊断的数学方法和/或算法相结合,和/或与所述期望的或基准负载-桨距或桨距-负载关系的计算相结合,和/或与负载和/或桨距故障的隔离和/或大小估计和/或调节相结合。
图1a中所示的三个叶片B1,B2,B3和/或叶片传感器使用本发明的方法进行诊断。叶片的虚线位置指示出等于零度的方位角(δ=0°)。风轮机的结构不局限于三个叶片B1,B2,B3。也可使用更少或更多数目的叶片。图1b显示了叶片B1的叶片横截面,横截面具有桨距θ,拍打方向负载转矩N1以及拍打方向运动(以虚线/箭头显示)。
根据本发明的诊断方法和系统基于叶片模型,当满足以下条件时,该模型是最准确的:
-没有显著的空气动力学力作用于叶片B1,B2,B3上,由此低风速或者零风速;和/或
-静止转子状态。
当目的在于隔离桨距和/或负载缺陷以及估计桨距和/或负载故障大小时,包括风速与风轮机产生的功率之间关系的模型在实际当中被认为是太不确定而不能给出有用的结果。
第二个需要通常是低风速的结果。当空气动力学负载是可忽略时,重力负载是影响叶片B1,B2,B3的唯一负载。这里给出了描述如图2所示变量之间正常行为的等式集合(约束)的结构,其中黑色方框代表约束以及圆圈代表变量。已知变量,例如测量值和基准值,被显示成灰色背景。当i=1,2或3时,图2显示出叶片上的负载转矩N1i通过约束C1i而取决于叶片桨距θ1以及方位角δ,约束C1i被称为负载模型。负载测量值mN1i通过约束M1i而与物理负载N1i相关,约束M1i被称为测量约束(constraint)。系统的输入是通过约束Cθi控制桨距的桨距基准θref,i。方位角δ通过闭环(被描述成约束Cδ)被控制成与设定点δref匹配。注意的是,关于如何测量桨距和负载(传感器类型等)没有做出任何假定。对于如何控制桨距和方位角没有做出任何假定,仅仅是沿循着假定没有故障存在时稳定状态下的基准值。
叶片和相关传感器的数学模型的分析显示出负载和/或桨距故障是可以检测到的,但是不能够通过源于故障检测和隔离理论的传统方法而隔离。
I表示可隔离集合以及D表示可检测故障的集合(违背约束),
I={Cδ}
D1={ml1,Cθ1,Cl1}(1)
D2={ml2,Cθ2,Cl2}
D3={ml3,Cθ3,Cl3}
这些可隔离与可检测故障的组合是被动诊断系统能够实现的。方位致动器Cδ的故障可以被隔离,而其余集合是可检测的,即任何关系{m1i,Cθi,C1i}中的故障将只能够被检测到,其中i=1,2或者3。术语“被动诊断系统”意味着没有系统输入被诊断系统所改变。被动诊断系统能够仅仅证实闭环的一致性,而不会改变行为和/或基准值。然而,如果诊断允许轻微地改变或者摄动输入,则能够将一些可检测组合划分成更小的组合(或者甚至使约束可以隔离)。这是主动故障隔离的概念。所描述的模型(零风速以及静止风轮机)特别地关注主动故障隔离。当模型有效时,能够实现主动故障隔离而不会影响风轮机发电量。
在图2的结构化模型中,可以看到从输入到输出的路径,从θref以及δrefmN1是不同的。路径θrefmN1是{Cθ,C1,m1}。路径δrefmN1是{Cδ,C1,m1}。由于这些不同,如果输入被给出了不同的基准信号或者任一基准信号被给出摄动,则存在使约束中故障隔离的可能性。这被称为主动故障隔离。
为了通过详细分析来评估该结果,通过下式给出负载测量:
N l m = g l ( δ , θ ) + f N l = f A g l ( δ ref , θ ref - f θ ) + f N l - - - ( 2 )
其中负载模型gl取决于方位角和桨距角,fA是负载测量的增益故障,
Figure GPA00001146419900082
是负载偏离或者负载传感器的偏离,fθ是桨距角或者桨距传感器的偏离。扩展到包括桨距增益故障是简单的,但是为了简洁它没有包括在此。测量到的负载相对于等式2中三个故障中任一个的灵敏度给出如下:
∂ m N l ∂ f θ = f A ∂ ∂ f θ g l ( δ ref , θ ref + f θ ) - - - ( 3 a )
∂ m N l ∂ f N l = 1 - - - ( 3 b )
∂ m N l ∂ f A = g l ( δ ref , θ ref + f θ ) - - - ( 3 c )
通过上面等式3a,3b和3c,可以清楚可以在同一输出中看到所有三个故障(没有任何一个灵敏度是零)。更重要地,桨距故障的灵敏度取决于能够实现主动故障隔离的系统输入(δref和θref)。负载增益故障fA灵敏度也取决于输入,但是与fθ的分析表达式不同。全部三个故障可能以不同灵敏度(Eq.3a,3b和3c)影响同一输出。这使得能够进行主动故障隔离。
关于桨距故障的灵敏度取决于负载行为(模型),而负载测量偏离在测量负载中具有一致的灵敏度。假定准确的模型,g1,那么如果能够确定出适当的输入信号摄动序列,则全部三个故障都能以良好的准确度而被隔离。
如果方位角为零(δ=0),那么重力信号的模型不取决于桨距增益及偏离故障以及负载增益故障。由此,该信号中只有负载偏离而且由此能够在该位置轻易地估计出来。然而,负载偏离还可在不等于零的方位角位置估计出来。由此,没有必要在算法上将本发明限制在该位置(δ=0),但是该位置不应当被排除作为一个可能性。
使用基于残差的主动隔离可以得出相同的结论。这可以从残差表达式中看到:
R=mNl-glref,θref)(4)
由于项glref,θref)是预期行为,因此它不取决于故障。该测量不取决于以等式2描述方式的缺陷。因此,R与mN1的灵敏度是相同的,并且可以得出相同的结论。
概括地,主动故障方案能够在静止等于零风速条件下诊断出桨距和负载故障,例如桨距和负载测量偏离以及增益故障。该故障可以使用以下程序加以区别:
程序1:在不同方位角下循环桨距。
1)对将被诊断的风轮机叶片设定期望的(例如非零)方位角δref
2)运行桨距θref循环(例如从极值到极值)以及测量重力负载,从而推导出桨距-负载或者负载-桨距关系;
3)在不同方位角下运行步骤2),从而增加故障估计的鲁棒性;
4)将桨距-负载或者负载-桨距关系与预期的桨距-负载或者负载-桨距关系作比较;以及
5)从桨距-负载或者负载-桨距数据集合推导出桨距和/或负载测量故障(例如桨距增益或/和桨距偏离故障和/或负载增益或/和负载偏离故障)。
程序2:在不同桨距角下的循环方位角。
1)对将被诊断的风轮机叶片设定预期的桨距θref
2)运行方位角δref循环(例如从极值到极值)以及测量重力负载,从而推导出桨距-负载或者负载-桨距关系;
3)在不同桨距角下运行步骤2),从而增加故障估计的鲁棒性;
4)将桨距-负载或者负载-桨距关系与预期的桨距-负载或者负载-桨距关系作比较;以及
5)从桨距-负载或者负载-桨距数据集合推导出桨距和/或负载测量故障(例如桨距增益或/和桨距偏离故障和/或负载增益或/和负载偏离故障)。
程序3:在方位角旋转与桨距改变同时发生的情况下获得桨距及负载测量的数据集合,从而获得用于分析的数据集合。
这些程序中的任意一个都可以被看作是本发明用于找到桨距-负载或负载-桨距响应的目标。
桨距和方位角可以通过反馈作用(闭环控制)或者通过直接致动没有反馈(开环)而引起。
可能的故障可以是以下中的至少一种:桨距增益故障,桨距偏离故障,负载偏离故障,负载增益故障,和/或与叶片和/或传感器正常行为不同的其它偏差。
图3a,3b,3c和3d分别显示了桨距偏离、负载偏离、桨距增益故障和负载增益故障的示例。这四个曲线图中,轻(细)曲线(线条)显示出通过桨距和负载测量装置所测量出的,以及重(粗)曲线(线条)显示出所期望的。
下面显示如何通过刚性体模型以及所描述的程序来估计桨距偏离故障以及负载增益和偏离。然而,本发明并不局限于刚性体模型,该模型仅仅作为示例给出。任意其它适当的模型以及由此任意适当的算法都可以与本发明一起使用。
通过刚性体重力负载来解释主动故障隔离的可用性。拍打方向重力负载行为是:
N l m ( δ ref , θ ref ) = m b l b g sin ( δ ) sin ( θ ) + f N l
= f A m ‾ b l ‾ b g sin ( δ ref ) sin ( θ ref + f θ ) + f N l - - - ( 5 )
其中mb和lb是实际叶片质量和质心,θ和δ分别是实际的桨距和方位角。当使用标准数值(
Figure GPA00001146419900103
Figure GPA00001146419900104
)时,倍增故障fA被包括。增益故障指的是
Figure GPA00001146419900105
项中的故障,其可被解释成负载传感器增益故障。此外,负载和桨距测量偏离故障被包括。
测量信号具有以下关于桨距测量偏离的灵敏度:
∂ m N l ∂ f θ = f A m ‾ b l ‾ b g sin ( δ ref ) cos ( θ ref - f θ ) - - - ( 6 )
其中当方位角基准值被设定为零(δref=0)时,桨距偏离故障的相关性明显地消失。这能够直接地读取负载测量偏离以及应用上面所述的诊断程序。在诊断程序运行之后,下面讨论估计故障的任务(该方案假设没有桨距增益故障)。
桨距基准值以从极值到极值的方式改变,其中故障应当被估计。重力负载是周期性的(相对于桨距而言),具有360度的周期,但是振幅、相位和偏离都是未知的。使用通过精确线性化的正弦拟合,能够估计出相位、偏离和振幅。负载测量被假定为:
n=1…N    (7)
这是具有未知振幅(A),偏离(N0)以及相位
Figure GPA00001146419900111
但是具有已知频率(f)的正弦函数。测量噪声为(v)。等式直接展开为:
Figure GPA00001146419900112
其相对于
Figure GPA00001146419900113
Figure GPA00001146419900114
以及λ3=N0呈线性。这些参数可以使用线性回归而被估计出来:
N l m ( 1 ) N l m ( 2 ) · · · N l m ( N ) = sin ( 2 πf 1 T ) cos ( 2 πf 1 T ) 1 sin ( 2 πf 2 T ) cos ( 2 πf 2 T ) 1 · · · · · · · · · sin ( 2 πfNT ) cos ( 2 πfNT ) 1 λ 1 λ 2 λ 3 + v 1 v 2 · · · v N - - - ( 9 )
或者mNl=Φλ+v
最小均方差参数求解是一种对线性最小方差问题的解法并且具有闭式解法:
λ ^ = ( Φ T Φ ) - 1 Φ T m N l - - - ( 10 )
相位可以通过下式估计:
Figure GPA00001146419900117
桨距偏离故障可以直接地通过相位而被估计出来,不需要计算出模型响应相位,原因是期望信号的模型相位为零。负载测量偏离故障同样可以通过正弦拟合而被估计出来。
N ^ 0 = f ^ N l = λ ^ 3 - - - ( 12 )
可选择地,在δref=0下负载偏离估计可以通过该正弦拟合而被证实。信号幅度同样能够被估计出来。
A ^ = f ^ A m ‾ b l ‾ b g sin ( δ ref ) = λ ^ 1 2 + λ ^ 2 2 - - - ( 13 )
认识到测量增益误差fA可以从振幅估计中发现。当测量包含白噪声时,(10)的解法是相对于的最优最小均方差估计。
使用基于残差的诊断(代替上面讨论的输出行为)获得了残差结构:
R = m b l b sin ( δ ) sin ( θ ) - m ‾ b l ‾ b sin ( δ ref ) sin ( θ ref ) - - - ( 14 )
= m ‾ b l ‾ b sin ( δ ref ) ( f A cos ( f θ ) sin ( θ ref ) + f A sin ( f θ ) cos ( θ ref ) -sin ( θ ref ) ) + f N l - - - ( 15 )
其中可以在变换响应上进行相同的回归分析:
R + sin ( θ ref ) m ‾ b l ‾ b sin ( δ ref ) - - - ( 16 )
算法的功能性可以通过运行仿真而得到证实,该仿真具有下面的、从基准到负载的静态响应:
N l m = f A m ‾ bg l ‾ b sin ( δ ref ) sin ( θ ref + f θ + f N l + v ) - - - ( 17 )
v∈N(μ,σ)
其中
Figure GPA00001146419900122
fA=1.05,fθ=4°,以及测量噪声具有零均值并且方差σ2=(10kNm)2。主动隔离算法生成了连续方位角下的桨距-负载或者负载-桨距曲线。每个曲线都显示了在方位角从0°到90°跨度内(间隔为10°)的响应。
在方位角为零时,算法可以估计出负载测量的负载偏离作为平均值(当改变桨距时)。随后,算法使在方位角集合产生的桨距-负载或负载-桨距曲线进行正弦拟合,从而获得了取决于方位角的相位、振幅和偏离估计。对于故障诊断仅仅需要一个方位角,但是使用不同角度的集合将会增加方法的鲁棒性。
样本的故障估计曲线在图4a,4b和4c中显示。图4b显示了负载偏离估计,其估计出25kNm的偏离。测量噪声影响了负载偏离估计,但是显示出没有显著的方位角相关性。图4a显示了桨距偏离估计,其在方位角大于10°时估计出4°偏离。正弦拟合的振幅与项
Figure GPA00001146419900125
相对应并且通过
Figure GPA00001146419900126
对估计进行缩放,从而获得如图4c所示的精确故障估计。图线是方位角的函数并且证实了该估计中没有方位角相关性。
知道的是,本发明的特定实施例使用了δref=δ并且获得了与上面所述相同的结果。
负载行为的分析因素证实了桨距偏离以及负载增益及偏离故障可以通过主动隔离而被隔离。通过用于刚性体的重力负载,已经进行仿真证实了这些结论。已经作出最小均方差参数估计,用于估计桨距和负载偏离以及负载测量增益故障。它们全部可以通过最小均方差估计和主动隔离而诊断出来。该诊断程序仅仅使用来自于一个叶片的测量,由此同一方案可以独立地用于其它轮机叶片上。
然而本发明不局限于上述仿真或算法,而是建立在其原理之上。
根据本发明的程序已经得出,能够仅仅通过被负载测量偏离和增益故障所破坏的负载测量而诊断出桨距偏离。通过变频正弦拟合估计量算法,扩展到诊断桨距增益故障的算法容易完成。
尽管本发明已经结合优选实施例进行了描述,但是并不意味着被限定为在此提出的特定形式。相反地,本发明的范围仅仅被附加权利要求所限定。
在该部分,所公开实施例的某一具体细节是出于解释的目的而非限制目的提出,从而提供了本发明清楚且全面的理解。然而,本领域技术人员应当容易理解到,本发明可以在没有完全符合在此提出细节的其它实施例中实施,而不会显著地脱离该公开的精神与范围。此外,在行文中,以及出于简要清楚的目的,对公知设备、电路和方法学的详细描述进行了省略,从而避免不必要的细节和可能的混淆。
在权利要求中,术语“包括”没有排除其它元件或步骤的存在。此外,尽管单个特征可以包含在不同权利要求中,但是这些特征可以有利地进行组合,并且包含在不同权利要求中并不意味着特征的组合是不可行的和/或不利的。此外,单数标记没有排除多个。由此,标记“一个”、“第一”、“第二”等没有排除多个。附图标记包含在权利要求中,然而,附图标记的包含仅仅为了清楚的缘故,并且不应当被认为是权利要求范围的限制。
术语表
拉丁字母
Figure GPA00001146419900141
希腊字母
Figure GPA00001146419900142
Figure GPA00001146419900151

Claims (8)

1.一种用于叶片和相关传感器的诊断的方法,所述方法的特征在于,它是用于与桨距和/或负载相关的至少一个可能故障的诊断的方法,该方法包括以下步骤:
a)通过使用下面两个程序中的任意一个或者两者来进行负载和桨距测量:
i)运行桨距(θ)循环并且测量负载,以及随后在一个或多个方位角(δ)下重复桨距(θ)循环测量步骤,和/或
ii)运行方位角(δ)循环并且测量负载,以及随后在一个或多个桨距基准角(θref)重复方位角(δ)循环测量步骤;
b)推导负载测量值与桨距角之间的关系;以及
c)将所述负载-桨距或者桨距-负载关系与预期的负载-桨距或者桨距-负载关系进行比较,从而确定出叶片(B1;B2;B3…)和/或叶片传感器(s)是否存在至少一个可能故障。
2.如权利要求1所述的方法,其特征在于,程序i)还包括在一个或多个方位角下测量在一组桨距角(θ)下的负载。
3.如权利要求1所述的方法,其特征在于,测量是在低风速或者零风速下进行。
4.如权利要求1或2所述的方法,其特征在于,测量是在静止转子状态下进行。
5.如权利要求1或2所述的方法,其特征在于,所述可能故障是以下中的至少一个:桨距增益故障、桨距偏离故障、负载偏离故障、和负载增益故障。
6.如权利要求1或2所述的方法,其特征在于,至少一个数学方法被用于所述至少一个可能故障的诊断。
7.如权利要求1或2所述的方法,其特征在于,所述诊断方法的结果被用于负载和/或桨距故障的诊断、和/或用于故障大小的估计、和/或用于故障的调节。
8.如权利要求1或2所述的方法,其特征在于,至少一个数学算法被用于所述至少一个可能故障的诊断。
CN2008801189827A 2007-11-07 2008-11-07 桨距和负载缺陷的诊断 Active CN101918710B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
EP07120177.6 2007-11-07
EP07120177 2007-11-07
US99085807P 2007-11-28 2007-11-28
US60/990,858 2007-11-28
PCT/DK2008/000395 WO2009059606A2 (en) 2007-11-07 2008-11-07 Diagnosis of pitch and load defects

Publications (2)

Publication Number Publication Date
CN101918710A CN101918710A (zh) 2010-12-15
CN101918710B true CN101918710B (zh) 2013-10-16

Family

ID=39865724

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008801189827A Active CN101918710B (zh) 2007-11-07 2008-11-07 桨距和负载缺陷的诊断

Country Status (6)

Country Link
US (1) US7954372B2 (zh)
EP (1) EP2225460B1 (zh)
CN (1) CN101918710B (zh)
AT (1) ATE528504T1 (zh)
ES (1) ES2374860T3 (zh)
WO (1) WO2009059606A2 (zh)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DK2386751T3 (en) 2010-05-12 2016-11-28 Siemens Ag Windmill
DE102010024251A1 (de) * 2010-06-18 2011-12-22 Robert Bosch Gmbh Verfahren und Vorrichtung zur Ermittlung eines Schätzwerts für zumindest eine Messgröße einer Windkraftanlage
GB2481461A (en) * 2010-06-21 2011-12-28 Vestas Wind Sys As Control of a downstream wind turbine in a wind park by sensing the wake turbulence of an upstream turbine
EP2588755B1 (en) 2010-06-29 2019-12-11 Vestas Wind Systems A/S Calibration of wind turbine sensor
GB2482038B (en) 2010-07-14 2014-07-23 Vestas Wind Sys As Ice detection method and system for wind turbine blades
GB2482009B (en) * 2010-07-14 2014-07-23 Vestas Wind Sys As Ice detection and system for wind turbine blades
CN102644545B (zh) * 2011-02-18 2013-07-31 华锐风电科技(集团)股份有限公司 风电机组故障处理方法及系统
US8757003B1 (en) 2011-12-15 2014-06-24 Shaw Shahriar Makaremi Multi-frequency-band blade condition monitoring system
US8511177B1 (en) 2011-12-15 2013-08-20 Shaw Shahriar Makaremi Blade condition monitoring system
CN104246466B (zh) * 2011-12-30 2018-01-02 维斯塔斯风力系统集团公司 估计和控制结构中经受的负荷
US9606234B2 (en) 2013-10-18 2017-03-28 Tramontane Technologies, Inc. Amplified optical circuit
ES2612208T3 (es) * 2014-01-21 2017-05-12 Ssb Wind Systems Gmbh & Co. Kg Sistema de medición del ángulo de paso para turbina eólica
CN105332855B (zh) 2014-06-11 2019-06-28 通用电气公司 用于风力涡轮机的控制方法和控制系统
WO2015192852A1 (en) * 2014-06-20 2015-12-23 Mita-Teknik A/S System for dynamic pitch control
RU2582906C1 (ru) * 2015-04-06 2016-04-27 Фонд поддержки научной, научно-технической и инновационной деятельности "Энергия без границ" Система мониторинга крутильных колебаний вращающегося валопровода турбоагрегата
US11022100B2 (en) 2015-12-17 2021-06-01 General Electric Company System and method for controlling wind turbines
CN107327367B (zh) * 2017-06-30 2018-12-07 北京金风科创风电设备有限公司 风力发电机组变桨异常的识别方法及装置
CN107355342B (zh) * 2017-06-30 2019-04-23 北京金风科创风电设备有限公司 风力发电机组变桨异常的识别方法及装置
JP6476250B1 (ja) 2017-08-29 2019-02-27 三菱重工業株式会社 風力発電装置の診断方法及び診断システム
CN108150355B (zh) * 2017-12-29 2020-04-21 华润新能源(临武)风能有限公司 滤波器的控制方法和滤波器的控制装置
CN112523946B (zh) * 2019-09-17 2022-10-18 北京金风科创风电设备有限公司 风力发电机组的变桨控制方法、装置及介质
CN113565699B (zh) * 2021-08-13 2023-09-12 远景智能国际私人投资有限公司 风力发电机组的桨距角的检测方法、装置及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0995904A2 (de) * 1998-10-20 2000-04-26 Tacke Windenergie GmbH Windkraftanlage
CN101022260A (zh) * 2006-02-15 2007-08-22 中国科学院自动化研究所 多风轮风电机组装置及其控制方法

Family Cites Families (16)

* 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
WO1992018706A1 (fr) * 1991-04-12 1992-10-29 Komatsu Ltd. Niveleuse pour bulldozer
US5155375A (en) * 1991-09-19 1992-10-13 U.S. Windpower, Inc. Speed control system for a variable speed wind turbine
KR0157255B1 (ko) * 1994-01-31 1998-11-16 김무 블래이드 각도 조정가능한 냉각팬을 이용한 엔진 냉각 시스템
DE19948194C2 (de) * 1999-10-06 2001-11-08 Aloys Wobben Verfahren zur Überwachung von Windenergieanlagen
US7015595B2 (en) * 2002-02-11 2006-03-21 Vestas Wind Systems A/S Variable speed wind turbine having a passive grid side rectifier with scalar power control and dependent pitch control
CA2426711C (en) * 2002-05-02 2009-11-17 General Electric Company Wind power plant, control arrangement for a wind power plant, and method for operating a wind power plant
US7160083B2 (en) * 2003-02-03 2007-01-09 General Electric Company Method and apparatus for wind turbine rotor load control
US7445431B2 (en) * 2003-02-18 2008-11-04 Forskningscenter Riso Method of controlling aerodynamic load of a wind turbine based on local blade flow measurement
JP4738206B2 (ja) * 2006-02-28 2011-08-03 三菱重工業株式会社 風力発電システム、及びその制御方法
US8174136B2 (en) * 2006-04-26 2012-05-08 Alliance For Sustainable Energy, Llc Adaptive pitch control for variable speed wind turbines
DE112007001675A5 (de) * 2006-05-15 2009-04-16 Igus - Innovative Technische Systeme Gmbh Verfahren zur Überwachung der Beanspruchung von Rotorblättern von Windkraftanlagen
JP2010514978A (ja) * 2006-12-28 2010-05-06 クリッパー・ウィンドパワー・テクノロジー・インコーポレーテッド 風力タービンにおける、推定方法を使用するタワー共振運動および対称ブレード運動の減衰
US8334610B2 (en) * 2009-02-13 2012-12-18 Robert Migliori Gearless pitch control mechanism for starting, stopping and regulating the power output of wind turbines without the use of a brake
US20100226772A1 (en) * 2009-02-25 2010-09-09 Kenneth James Deering Blade control system
US7763989B2 (en) * 2009-07-07 2010-07-27 General Electric Company Method and apparatus for controlling the tip speed of a blade of a wind turbine

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0995904A2 (de) * 1998-10-20 2000-04-26 Tacke Windenergie GmbH Windkraftanlage
CN101022260A (zh) * 2006-02-15 2007-08-22 中国科学院自动化研究所 多风轮风电机组装置及其控制方法

Also Published As

Publication number Publication date
EP2225460B1 (en) 2011-10-12
ATE528504T1 (de) 2011-10-15
US20110001319A1 (en) 2011-01-06
US7954372B2 (en) 2011-06-07
WO2009059606A3 (en) 2009-12-23
WO2009059606A2 (en) 2009-05-14
ES2374860T3 (es) 2012-02-22
CN101918710A (zh) 2010-12-15
EP2225460A2 (en) 2010-09-08

Similar Documents

Publication Publication Date Title
CN101918710B (zh) 桨距和负载缺陷的诊断
US10072635B2 (en) Wind turbine and method for evaluating health state of blade thereof
US10883475B2 (en) Method for monitoring and assessing power performance changes of a wind turbine
US9240083B2 (en) Rotor system health monitoring using shaft load measurements and virtual monitoring of loads
CN103206342B (zh) 桨叶载荷传感器的标定
CN101821689B (zh) 用于为控制和/或监测应用提供至少一种输入传感器信号的方法和设备以及控制设备
EP2026160B1 (en) Event monitoring via combination of signals
US11519458B2 (en) Model-based method and system for monitoring the condition of a sliding bearing, particularly for wind turbines
CN108291527B (zh) 一种监测和评估风力涡轮机功率性能变化的方法
CN105891546B (zh) 基于大数据的风电机组偏航系统中风向标故障诊断的方法
US11313351B2 (en) Methods and systems of advanced yaw control of a wind turbine
JP7240527B2 (ja) 運転可能な風力タービンのコンポーネントの残存有効寿命を評価するための方法
Perišić et al. Cost‐effective shaft torque observer for condition monitoring of wind turbines
JP2017525891A (ja) 駆動システムの早期エラー検出方法、早期エラー検出システム、早期エラー検出システムを備える風力発電機、および早期エラー検出システムの使用
CN105527112A (zh) 一种基于使用与维修影响的旋转机械设备健康状态综合评估方法
CN113227926A (zh) 用于监测电机的状况监测装置和方法
CN111980855A (zh) 风力发电机组的偏航控制方法、装置、设备及存储介质
Badihi et al. Model-based active fault-tolerant cooperative control in an offshore wind farm
CN115135872A (zh) 确定风力涡轮机塔架的摆动运动的频率
Badihi et al. Model-free active fault-tolerant cooperative control in an offshore wind farm
BR102020006909A2 (pt) Aparato, sistema de monitoramento e método para avaliação da saúde de operação de máquinas movidas por motores elétricos
Ratkovac Rotor thrust force reconstruction for fatigue analysis of wind turbine support structures using modal identification based on long-term monitoring data
CN117684946A (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