CN112924173A - 一种风力发电机组变桨轴承故障诊断方法 - Google Patents

一种风力发电机组变桨轴承故障诊断方法 Download PDF

Info

Publication number
CN112924173A
CN112924173A CN202110224981.8A CN202110224981A CN112924173A CN 112924173 A CN112924173 A CN 112924173A CN 202110224981 A CN202110224981 A CN 202110224981A CN 112924173 A CN112924173 A CN 112924173A
Authority
CN
China
Prior art keywords
pitch bearing
generating set
wind generating
fault diagnosis
matrix
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.)
Pending
Application number
CN202110224981.8A
Other languages
English (en)
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.)
Guangzhou Special Pressure Equipment Inspection and Research Institute
Original Assignee
Guangzhou Special Pressure Equipment Inspection and Research Institute
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 Guangzhou Special Pressure Equipment Inspection and Research Institute filed Critical Guangzhou Special Pressure Equipment Inspection and Research Institute
Priority to CN202110224981.8A priority Critical patent/CN112924173A/zh
Publication of CN112924173A publication Critical patent/CN112924173A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Wind Motors (AREA)

Abstract

本发明涉及一种风力发电机组变桨轴承故障诊断方法。本发明所述的风力发电机组变桨轴承故障诊断方法包括如下步骤:获取风力发电机组的运行状态参数;通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的变桨轴承的运行状态敏感的P种状态参数;将所述P种状态参数输入训练好的所述变桨轴承故障诊断网络模型,得到所述风力发电机组的变桨轴承故障诊断结果,其中,所述变桨轴承故障诊断网络模型为单隐含层前馈神经网络模型。本发明所述的一种风力发电机组变桨轴承故障诊断方法具有快速、在线、无损、实时、高效监测和诊断故障的优点。

Description

一种风力发电机组变桨轴承故障诊断方法
技术领域
本发明涉及风电领域,特别是涉及一种风力发电机组变桨轴承故障诊断方法。
背景技术
变桨轴承作为水平轴风力发电机组的重要组成部分,主要功能是通过调节桨距角进行功率调节。其中,变桨轴承又是风力发电机组变桨轴承的关键部件,在工作中同时承受轴向载荷、径向载荷及倾覆力矩的影响,工况载荷情况复杂、工况环境差及整体刚度较低。因此变桨轴承的使用性能,决定着风力发电机组整体的稳定性和安全性。
风电场调研分析表明,随着机组运行时间的增加,变桨轴承故障概率呈增长趋势。例如,某风电场的各台风力发电机组陆续运行时间已有7-10年,变桨轴承出现多起严重事故:发生21起导致叶片坠落的事故中,其中9起出现变桨轴承外圈断裂;在运维检查中发现近10台次轴承外圈开裂,与叶片坠落时轴承断裂情况类似;发生1起轴承卡涩导致叶轮无法回桨、高速轴过速制动引发机组着火事故;同时还发现50多例较为严重的磨损卡涩问题。该型机组变桨轴承的损坏主要表现为轴承外圈断裂和轴承轨道磨损卡涩问题,给风电场的安全生产带来极大的隐患。
因此,有必要不断研发和改进风力发电机组变桨轴承在线状态监测与故障诊断方法,以便及时准确地发现变桨轴承故障,并根据故障诊断结果安排对应的检修计划,减少设备事故率,降低维修费用,大大提高风力机运行安全性。
发明内容
基于此,本发明的目的在于,提供一种风力发电机组变桨轴承故障诊断方法,其具有快速、在线、无损、实时、高效监测和诊断故障的优点。
一种风力发电机组变桨轴承故障诊断方法,其特征在于,包括如下步骤:
获取风力发电机组的运行状态参数;
通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的运行状态敏感的P种状态参数,其中,P为所述运行状态参数的输入维数;
将所述P种状态参数输入训练好的所述变桨轴承故障诊断网络模型,得到所述风力发电机组的变桨轴承故障诊断结果,其中,所述变桨轴承故障诊断网络模型为单隐含层前馈神经网络模型。
本发明所述的风力发电机组变桨轴承故障诊断方法,通过实时获取风电机组变桨轴承的运行状态参数及实时诊断,可及时准确地发现变桨轴承故障,并根据故障诊断结果安排对应的检修计划,减少设备事故率,降低维修费用,大大提高风力机运行安全性。
进一步地,所述通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的运行状态敏感的P种状态参数,包括:
1)标准化处理:
构建状态参数数据矩阵X:
Figure BDA0002956927670000021
通过如下公式计算协方差矩阵xi'j
Figure BDA0002956927670000022
(i=1,2,...,n;j=1,2,...,m)
其中,xnm表示第m组参数中的第n个状态参数;
标准化后的数据矩阵X'为:
Figure BDA0002956927670000023
X'协方差矩阵为:
Figure BDA0002956927670000024
2)利用协方差矩阵S对风力发电机组的运行状态参数进行降维处理:
将协方差矩阵S进行特征值分解,求出特征值及其特征向量:
S=PΛPT
式中,Λ是对角矩阵,包含递减的非负实特征值λ1≥λ2≥...≥λm≥0,P为特征值降序排列后前m行组成的向量矩阵;
3)计算主成分贡献率:
将最大特征值及其对应的特征向量,定义为第一主成分的方差和方向,以此类推,直至确定最后一个主成分的方差和方向;
每个主成分的方差在样本总方差中所占比值,即为该主成分对于该样本的贡献率,第k个主成分的贡献率Vk定义为:
Figure BDA0002956927670000031
将特征值λm从大到小排列,前k个主成分的累计贡献率CPV(k)定义为:
Figure BDA0002956927670000032
累计贡献率达到CPV阈值时主成分k的个数,即为降维后风力发电机组的运行状态参数的输入维数P。
通过上述公式处理,可从多种所述风力发电机组的变桨轴承运行状态参数中,获得所述运行状态参数的输入维数P。
进一步地,所述CPV阈值为90%,保证降维后风力发电机组的运行状态参数能包括原始状态参数的绝大多数信息。
进一步地,所述变桨轴承故障诊断网络模型包括一个输入层、一个隐含层和一个输出层;
输入层节点数为n个,将所述P种状态参数作为输入参数,即n=P;
隐含层节点数为i;
输出层节点数为3,分别表示变桨轴承的三个典型状态:正常状态、部件磨损状态、部件断裂状态。
进一步地,所述变桨轴承故障诊断网络模型的训练过程包括如下步骤:
获取风力发电机组的运行状态参数;
通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的运行状态敏感的P种状态参数,其中,P为所述运行状态参数的输入维数;
通过OS-W-ELM训练算法来训练所述变桨轴承故障诊断网络模型,包括以下阶段:
1)OS-W-ELM初始化学习阶段:
初始阶段选取N0组数据
Figure BDA0002956927670000033
其中
Figure BDA0002956927670000034
其中,xi是一个数组,其包含了一组特征值,指同一时刻,所有的风力机运行状态参数;
ti为这一组特征值所对应的时刻;
1-1)随机产生输出权重矩阵wi和偏置矩阵bi,其中,
Figure BDA0002956927670000041
1-2)通过如下公式计算初始的隐含层输出矩阵:
Figure BDA0002956927670000042
其中,g(x)为激活函数,g(w1x1+b1)是隐含层输出矩阵中的值,根据随机的输出权重矩阵wi和偏置矩阵bi获得,然后在神经网络自我学习过程中自动修正输出权重矩阵wi和偏置矩阵bi
1-3)通过如下公式计算初始的输出权重矩阵β0
Figure BDA0002956927670000043
其中,
Figure BDA0002956927670000044
其中,M0
Figure BDA0002956927670000045
计算得到,T0是单独提出的时间矩阵,对应
Figure BDA0002956927670000046
中的ti
1-4)设k=0,k为块数,表示初始化学习阶段;
2)OS-W-ELM在线序列学习阶段:
2-1)设第k+1块样本集为:
Figure BDA0002956927670000047
其中,Nj是第j块样本集的样本数量;
通过如下公式计算隐含层输出矩阵Hk+1
Figure BDA0002956927670000048
2-2)通过如下公式计算输出权重矩阵βk+1
Figure BDA0002956927670000051
其中,
Figure BDA0002956927670000052
Mk+1通过上述公式计算得到,用于不断地修正权重矩阵βk+1,使神经网络能达到自我学习的目的;
2-3)令k=k+1,转至步骤2-1)直至结束。
与传统的训练方法相比,所述OS-W-ELM训练算法具有学习速度快、泛化性能好等优点。通过学习的诊断网络,就可以用于变桨轴承的故障诊断。在网络诊断使用过程中,若发现有新的故障数据样本,采用在线学习方法,快速学习新的样本,无需重新开始学习,大大地节约了时间成本,提高了系统运行效率。
进一步地,所述风力发电机组的运行状态参数包括以下至少一项:
风力发电机机舱外平均风速、风向、风力发电机第一桨叶至第三桨叶的桨距角、低速轴转速、发电机有功功率、变桨电机的电压和电流、变桨电机温度、桨叶变流器温度、变桨电池电压、变桨轴承油脂泵出口油压和变桨齿轮油脂泵出口油压。
以上参数基本覆盖风力机电动变桨系统常见故障所需的运行状态参数,可通过常用传感器获取,不额外增加成本。
进一步地,所述风力发电机组的运行状态参数通过SCADA系统获得,所述SCADA系统包括信号传感器、数据采集器、SCADA网络和故障诊断终端;
所述信号传感器设置于风力发电机组上,并与所述数据采集器连接;所述数据采集器和所述故障诊断终端分别接入所述SCADA网络;
所述信号传感器设置于风力发电机组上,用于感应所述风力发电机组的运行状态,并将感应到的参数信号传输至所述数据采集器;
所述数据采集器用于将所述参数信号转换为风力发电机组的运行状态参数信号,并通过所述SCADA网络发送给所述故障诊断终端;
所述故障诊断终端搭载信号处理软件和检测诊断软件系统,用于通过所述SCADA网络获取所述变桨轴承的运行状态参数,并输出风力发电机组变桨轴承故障结果。
所述SCADA系统可实现所述风力发电机组的运行状态参数的实时获取、传输和处理计算,有效提高风力发电机组运行的安全性和可靠性。
进一步地,所述SCADA系统还包括信号放大器,所述信号放大器分别与所述传感器和所述数据采集器连接,所述信号放大器用于接收所述传感器的参数信号,对所述参数信号进行滤波、除噪和放大处理后,传输给所述数据采集器。
所述信号放大器可将所述传感器感应到的参数信号进行标准化处理,便于数据采集器进行下一步的处理与转换。
进一步地,所述SCADA系统还包括数据库服务器和web服务器,所述数据库服务器与所述web服务器接入所述SCADA网络,用于接收并存储所述风力发电机组的运行状态参数。
所述数据库服务器和所述web服务器可实现所述风力发电机组的运行状态参数的实时存储,并进一步便于与故障诊断终端的数据交互。
进一步地,所述SCADA系统还包括监测终端;
所述监测终端接入所述SCADA网络,
或者,
所述监测终端通过Internet与所述web服务器连接,并通过所述web服务器获取所述风力发电机组的运行状态。
用户可灵活设置监测终端的位置,即时在线监视风力发电机组的变桨系统的运行情况。
为了更好地理解和实施,下面结合附图详细说明本发明。
附图说明
图1为风力发电机组变桨轴承故障诊断方法步骤图;
图2为风力发电机组的运行状态参数的示意图;
图3为基于PCA的数据处理流程图;
图4为基于极限学习机的变桨轴承故障诊断网络模型的结构示意图;
图5为基于OS-W-ELM训练算法训练变桨轴承故障诊断网络模型的流程图;
图6为变桨轴承故障诊断软件功能结构示意图;
图7为本发明一个实施例中的SCADA系统结构示意图;
图8为本发明的另一个实施例中的SCADA系统结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容。
在本发明的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶部”、“底部”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
需要说明的是,当元件被称为“固定于”另一个元件,它可以是直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件,它可以是直接连接到另一个元件或者可能同时存在居中元件。
下面给出几个具体的实施例,用于详细介绍本申请的技术方案。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。
如图1所示,在一个具体的实施例中,本发明所提供的一种风力发电机组变桨轴承故障诊断方法包括如下步骤:
S101:获取风力发电机组的运行状态参数。
其中,如图2所示,在一个实施例中,风力发电机组的运行状态参数包括风力发电机机舱外平均风速、风向、风力发电机第一桨叶至第三桨叶的桨距角、低速轴转速、发电机有功功率、变桨电机的电压和电流、变桨电机温度、桨叶变流器温度、变桨电池电压、变桨轴承油脂泵出口油压和变桨齿轮油脂泵出口油压中的至少一项。
在一个优选的实施例中,本发明提供的一种风力发电机组变桨轴承故障诊断方法通过SCADA系统获取以上参数。具体过程为:
1)将风速传感器、风向传感器、叶片桨距角传感器、低速轴转速传感器、发电机电功率传感器、变桨电机的电压和电流传感器、变桨电机温度传感器、桨叶变流器箱温度传感器、变桨电池电压传感器、变桨轴承油脂泵出口油压传感器、变桨齿轮油脂泵出口油压传感器接入SCADA系统。
2)SCADA系统实时采集各个传感器信号,并通过SCADA系统自身的信号转换系统,将传感器信号(电流或电压信号)转换成风力发电机组的运行状态参数(包括但不限于风速值、风向、叶片桨距角度、低速轴转速、发电机电功率、变桨电机的电压和电流值、变桨电机温度、桨叶变流器箱温度、变桨电池电压值、变桨轴承油脂泵出口油压值、变桨齿轮油脂泵出口油压值),并将风力发电机组的运行状态参数存储在数据库服务器中。
3)利用专门的数据读取程序读取风电机组SCADA系统数据库服务器中的实时风力机运行状态参数数据,从中提取风速、风向、叶片桨距角、低速轴转速、发电机电功率、变桨电机的电压和电流、变桨电机温度、桨叶变流器箱温度、变桨电池电压、变桨轴承油脂泵出口油压、变桨齿轮油脂泵出口油压等数据。
S102:通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的变桨轴承的运行状态敏感的P种状态参数,如图3所示,具体包括如下步骤:
1)标准化处理:
构建状态参数数据矩阵X:
Figure BDA0002956927670000081
通过如下公式计算协方差矩阵xi'j
Figure BDA0002956927670000082
(i=1,2,...,n;j=1,2,...,m)
其中,xnm表示第m组参数或在m时刻的第n个状态参数;
标准化后的数据矩阵X'为:
Figure BDA0002956927670000083
X'协方差矩阵为:
Figure BDA0002956927670000084
2)利用协方差矩阵S对风力发电机组的运行状态参数进行降维处理:
将协方差矩阵S进行特征值分解,求出特征值及其特征向量:
S=PΛPT
式中,Λ是对角矩阵,包含递减的非负实特征值λ1≥λ2≥...≥λm≥0,P为特征值降序排列后前m行组成的向量矩阵;
3)计算主成分贡献率:
将最大特征值及其对应的特征向量,定义为第一主成分的方差和方向,以此类推,直至确定最后一个主成分的方差和方向;
每个主成分的方差在样本总方差中所占比值,即为该主成分对于该样本的贡献率,第k个主成分的贡献率Vk定义为:
Figure BDA0002956927670000091
将特征值λm从大到小排列,前k个主成分的累计贡献率CPV(k)定义为:
Figure BDA0002956927670000092
累计贡献率达到CPV阈值时主成分k的个数,即为降维后风力发电机组的运行状态参数的输入维数P。
在一个优选的实施例中,CPV阈值为90%,保证降维后的风力发电机组的运行状态参数能包括原始状态参数的绝大多数信息。
S103:将所述P种状态参数输入训练好的所述变桨轴承故障诊断网络模型,得到所述风力发电机组的变桨轴承故障诊断结果,其中,所述变桨轴承故障诊断网络模型为单隐含层前馈神经网络模型。
优选地,变桨轴承故障诊断网络模型为单隐含层前馈神经网络模型,如图4所示,变桨轴承故障诊断网络模型的具体结构为:
包括一个输入层、一个隐含层和一个输出层;
所述输入层的节点数为n个,将所述P种状态参数作为输入参数,即n=P;
所述隐含层的节点数为i;
输出层节点数为3,分别表示变桨轴承的三个典型状态:正常状态、部件磨损状态、部件断裂状态。
其中,隐含层又称隐藏层或中间层。
优选地,采用OS-W-ELM训练算法训练变桨轴承故障诊断网络模型,如图5所示,具体包括以下阶段:
1)OS-W-ELM初始化学习阶段:
初始阶段选取N0组数据
Figure BDA0002956927670000093
其中
Figure BDA0002956927670000095
其中,xi是一个数组,其包含了一组特征值,指同一时刻,所有的风力机运行状态参数;
ti为这一组特征值所对应的时刻;
1-1)随机产生输出权重矩阵wi和偏置矩阵bi,其中,
Figure BDA0002956927670000094
1-2)通过如下公式计算初始的隐含层输出矩阵:
Figure BDA0002956927670000101
其中,g(x)为激活函数,g(w1x1+b1)是隐含层输出矩阵中的值,根据随机的输出权重矩阵wi和偏置矩阵bi获得,然后在神经网络自我学习过程中自动修正输出权重矩阵wi和偏置矩阵bi
常用的激活有函数Sigmoid函数、Tanh函数、ReLU函数等。
1-3)通过如下公式计算初始的输出权重矩阵β0
Figure BDA0002956927670000102
其中,
Figure BDA0002956927670000103
其中,M0
Figure BDA0002956927670000104
计算得到,T0是单独提出的时间矩阵,对应
Figure BDA0002956927670000105
中的ti
1-4)设k=0,k为块数,表示初始化学习阶段;
2)OS-W-ELM在线序列学习阶段:
2-1)设第k+1块样本集为:
Figure BDA0002956927670000106
其中,Nj是第j块样本集的样本数量;
通过如下公式计算隐含层输出矩阵Hk+1
Figure BDA0002956927670000107
2-2)通过如下公式计算输出权重矩阵βk+1
Figure BDA0002956927670000111
其中,
Figure BDA0002956927670000112
Mk+1通过上述公式计算得到,用于不断地修正权重矩阵βk+1,使神经网络能达到自我学习的目的;
2-3)令k=k+1,转至步骤2-1)直至结束。
与传统的训练方法相比,所述OS-W-ELM训练算法具有学习速度快、泛化性能好等优点。通过学习的诊断网络,就可以用于变桨轴承的故障诊断。在网络诊断使用过程中,若发现有新的故障数据样本,采用在线学习方法,快速学习新的样本,无需重新开始学习,大大地节约了时间成本,提高了系统运行效率。
在一个优选的实施例中,SCADA系统设置有计算机作为故障诊断终端,用于将P种状态参数输入训练好的变桨轴承故障诊断网络模型,得到风力发电机组的变桨轴承故障诊断结果的功能。
优选地,计算机内部安装专用的信号处理软件和轴承故障诊断软件,用于对风力发电机组的变桨轴承运行状态数字信号进行一系列处理与变换,对变桨系统轴承的运行状态进行监测诊断,对数据进行存储与查询、显示。
优选地,如图6所示,计算机搭载的轴承故障诊断软件的基本功能有:
(1)用户登录:用户通过设定密码登入系统。
(2)参数设置:系统中关于数据采集中各通道的采样参数该如何设置进行描述。
(3)系统维护:通过该模块,可以实现对系统的维护。
(4)帮助:对系统使用过程可能出现的问题及其处理方法进行解释;对系统进行维护。
(5)信号曲线:对系统采集到的原始信号,可绘制信号对时间的关系曲线,两个原始信号之间的相互关系曲线。
(6)特征曲线:对系统提取的特征参数,可绘制特征参数对时间的关系曲线,两个(或多个)特征参数之间的相互关系曲线。
(7)PCA分析:对采集到的原始数据进行主成分分析,提取与轴承状态密切相关的P各状态参数。
(8)特征提取:从敏感状态参数中提取故障特征,并将特征数据存入数据库。
(9)网络训练:利用上述(7)-(8)功能模块中提取的数据样本,对诊断网络进行离线和在线训练。
(10)自动诊断:系统一旦启动,就进入自动诊断模式,诊断的时间间隔可以设定。
(11)报告生成:当系统诊断到变桨系统轴承发生故障,就自动记录下故障时间和故障类型、故障严重程度,并自动生成诊断报告。
(12)数据保存:系统可诊断报告进行保存,存入数据库。
如图7和图8所示,在一个具体的实施例中,本发明所提供的一种风力发电机组变桨轴承故障诊断方法通过SCADA系统进行数据的获取、传输、保存和输出功能。SCADA系统包括信号传感器、数据采集器、SCADA网络和故障诊断终端。
信号传感器设置于风力发电机组上,并与数据采集器连接;数据采集器和故障诊断终端分别接入SCADA网络;
信号传感器设置于风力发电机组上,用于感应风力发电机组的运行状态,并将感应到的参数信号传输至数据采集器;数据采集器用于将参数信号转换为风力发电机组的运行状态信号,并通过SCADA网络发送给故障诊断终端;故障诊断终端搭载信号处理软件和检测诊断软件系统,用于通过SCADA网络获取风力发电机组的运行状态参数,并输出风力发电机组变桨轴承故障结果。
在一个优选的实施例中,SCADA系统还包括信号放大器,信号放大器分别与传感器和数据采集器连接,用于接收传感器的参数信号,对参数信号进行滤波、除噪和放大处理后,传输给数据采集器。
信号放大器可将传感器感应到的参数信号进行标准化处理,便于数据采集器进行下一步的处理与转换。
优选地,SCADA系统还包括数据库服务器和web服务器,数据库服务器与web服务器接入SCADA网络,用于接收并存储风力发电机组的运行状态参数。
数据库服务器和web服务器可实现风力发电机组的运行状态参数的实时存储,并进一步便于与故障诊断终端的数据交互。
在一个优选的实施例中,SCADA系统还包括监测终端。监测终端接入SCADA网络,或者,监测终端通过Internet与web服务器连接,并通过web服务器获取风力发电机组的运行状态。
用户可灵活设置监测终端的位置,即时在线监视风力发电机组的变桨系统的运行情况。
本发明采用基于主成分分析法(Principal Component Analysis,简称PCA)的数据处理方法,对风力发电机组的运行状态参数降维后,提取了用于诊断变桨轴承故障的主要状态参数,在确保诊断精度的前提下简化了诊断网络的结构。
同时,本发明采用改进的在线序列加权极限学习机(Online Sequential ExtremeWeighted Learning Machine,简称OS-W-ELM)的轴承故障诊断网络学习算法,学习速度快、泛化性能好在网络诊断使用过程中,若发现有新的故障数据样本,采用在线学习方法,快速学习新的样本,无需重新开始学习,大大地节约了时间成本。
结合现有的风电机组SCADA系统,本发明提供的一种风力发电机组变桨轴承故障诊断方法能够快速、在线、无损、实时、高效监测和诊断风电机组变桨系统轴承故障,提高了系统运行效率,为保障风电机组运行的安全性和可靠性提供了有力技术支持。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。

Claims (10)

1.一种风力发电机组变桨轴承故障诊断方法,其特征在于,包括如下步骤:
获取风力发电机组的运行状态参数;
通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的变桨轴承的运行状态敏感的P种状态参数;
将所述P种状态参数输入训练好的变桨轴承故障诊断网络模型,得到所述风力发电机组的变桨轴承故障诊断结果,其中,所述变桨轴承故障诊断网络模型为单隐含层前馈神经网络模型。
2.根据权利要求1所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于,所述通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的变桨轴承的运行状态敏感的P种状态参数,包括如下步骤:
1)标准化处理:
构建状态参数数据矩阵X:
Figure FDA0002956927660000011
通过如下公式计算协方差矩阵x′ij
Figure FDA0002956927660000012
其中,xnm表示第m组参数或在m时刻的第n个状态参数;
标准化后的数据矩阵X'为:
Figure FDA0002956927660000013
X'协方差矩阵为:
Figure FDA0002956927660000014
2)利用协方差矩阵S对所述运行状态参数进行降维处理:
将协方差矩阵S进行特征值分解,求出特征值及其特征向量:
S=PΛPT
式中,Λ是对角矩阵,包含递减的非负实特征值λ1≥λ2≥...≥λm≥0,P为特征值降序排列后前m行组成的向量矩阵;
3)计算主成分贡献率:
将最大特征值及其对应的特征向量,定义为第一主成分的方差和方向,以此类推,直至确定最后一个主成分的方差和方向;
每个主成分的方差在样本总方差中所占比值,即为该主成分对于该样本的贡献率,第k个主成分的贡献率Vk定义为:
Figure FDA0002956927660000021
将特征值λm从大到小排列,前k个主成分的累计贡献率CPV(k)定义为:
Figure FDA0002956927660000022
累计贡献率达到CPV阈值时主成分k的个数,即为降维后所述运行状态参数的输入维数P。
3.根据权利要求2所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于:
所述CPV阈值为90%。
4.根据权利要求1所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于:
所述变桨轴承故障诊断网络模型包括输入层、隐含层和输出层;
所述输入层的节点数为n个,将所述P种状态参数作为输入参数,即n=P;
所述隐含层的节点数为i;
所述输出层的节点数为3,分别表示所述变桨轴承的三个典型状态:正常状态、部件磨损状态、部件断裂状态。
5.根据权利要求4所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于,所述变桨轴承故障诊断网络模型的训练过程包括如下步骤:
获取所述风力发电机组的运行状态参数;
通过PCA算法对所述运行状态参数进行处理,获得对所述风力发电机组的变桨轴承的运行状态敏感的P种状态参数;
通过OS-W-ELM训练算法训练所述变桨轴承故障诊断网络模型,包括以下阶段:
1)OS-W-ELM初始化学习阶段:
初始阶段选取N0组数据
Figure FDA0002956927660000031
其中
Figure FDA0002956927660000032
其中,xi是一个数组,其包含了一组特征值,指同一时刻,所有的所述运行状态参数;
ti为这一组特征值所对应的时刻;
1-1)随机产生输出权重矩阵wi和偏置矩阵bi,其中,
Figure FDA0002956927660000033
1-2)通过如下公式计算初始的隐含层输出矩阵H0
Figure FDA0002956927660000034
其中,g(x)为激活函数,g(w1x1+b1)是隐含层输出矩阵中的值,根据随机的输出权重矩阵wi和偏置矩阵bi获得,然后在神经网络自我学习过程中自动修正输出权重矩阵wi和偏置矩阵bi
1-3)通过如下公式计算初始的输出权重矩阵β0
Figure FDA0002956927660000035
其中,
Figure FDA0002956927660000036
其中,M0
Figure FDA0002956927660000037
计算得到,T0是单独提出的时间矩阵,对应
Figure FDA0002956927660000038
中的ti
1-4)设k=0,k为块数,表示初始化学习阶段;
2)OS-W-ELM在线序列学习阶段:
2-1)设第k+1块样本集为:
Figure FDA0002956927660000039
其中,Nj是第j块样本集的样本数量;
通过如下公式计算隐含层输出矩阵Hk+1
Figure FDA0002956927660000041
2-2)通过如下公式计算输出权重矩阵βk+1
Figure FDA0002956927660000042
其中,
Figure FDA0002956927660000043
Mk+1通过上述公式计算得到,用于不断地修正权重矩阵βk+1,使神经网络能达到自我学习的目的;
2-3)令k=k+1,转至步骤2-1)直至结束。
6.根据权利要求1所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于:
所述风力发电机组的运行状态参数包括以下至少一项:
风力发电机机舱外平均风速、风向、风力发电机第一桨叶至第三桨叶的桨距角、低速轴转速、发电机有功功率、变桨电机的电压和电流、变桨电机温度、桨叶变流器温度、变桨电池电压、变桨轴承油脂泵出口油压和变桨齿轮油脂泵出口油压。
7.根据权利要求1所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于:
所述风力发电机组的运行状态参数通过SCADA系统获得,所述SCADA系统包括信号传感器、数据采集器、SCADA网络和故障诊断终端;
所述信号传感器设置于风力发电机组上,并与所述数据采集器连接;所述数据采集器和所述故障诊断终端分别接入所述SCADA网络;
所述信号传感器设置于风力发电机组上,用于感应所述风力机电动变桨轴承的运行状态,并将感应到的参数信号传输至所述数据采集器;
所述数据采集器用于将所述参数信号转换为变桨轴承运行的状态信号,并通过所述SCADA网络发送给所述故障诊断终端;
所述故障诊断终端搭载信号处理软件和检测诊断软件系统,用于通过所述SCADA网络获取所述变桨轴承的运行状态参数,并输出风力发电机组变桨轴承故障结果。
8.根据权利要求7所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于:
所述SCADA系统还包括信号放大器,所述信号放大器分别与所述传感器和所述数据采集器连接,所述信号放大器用于接收所述传感器的参数信号,对所述参数信号进行滤波、除噪和放大处理后,传输给所述数据采集器。
9.根据权利要求7所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于:
所述SCADA系统还包括数据库服务器和web服务器,所述数据库服务器与所述web服务器接入所述SCADA网络,用于接收并存储所述变桨轴承运行的状态信号。
10.根据权利要求7所述的一种风力发电机组变桨轴承故障诊断方法,其特征在于:
所述SCADA系统还包括监测终端;
所述监测终端接入所述SCADA网络,
或者,
所述监测终端通过Internet与所述web服务器连接,并通过所述web服务器获取所述变桨轴承运行的状态信号。
CN202110224981.8A 2021-03-01 2021-03-01 一种风力发电机组变桨轴承故障诊断方法 Pending CN112924173A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110224981.8A CN112924173A (zh) 2021-03-01 2021-03-01 一种风力发电机组变桨轴承故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110224981.8A CN112924173A (zh) 2021-03-01 2021-03-01 一种风力发电机组变桨轴承故障诊断方法

Publications (1)

Publication Number Publication Date
CN112924173A true CN112924173A (zh) 2021-06-08

Family

ID=76172732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110224981.8A Pending CN112924173A (zh) 2021-03-01 2021-03-01 一种风力发电机组变桨轴承故障诊断方法

Country Status (1)

Country Link
CN (1) CN112924173A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113312731A (zh) * 2021-06-28 2021-08-27 北京南洋思源智能科技有限公司 一种变桨轴承故障检测方法及装置、存储介质
CN113326589A (zh) * 2021-06-28 2021-08-31 北京南洋思源智能科技有限公司 一种敏感特征筛选方法及装置、存储介质
CN113343579A (zh) * 2021-06-25 2021-09-03 海南小鲨鱼智能科技有限公司 一种电池输出功率的控制方法、装置及计算机存储介质
CN113898528A (zh) * 2021-09-30 2022-01-07 江苏金风软件技术有限公司 风机变桨轴承的异常检测方法、模型训练方法及相关装置
CN114235415A (zh) * 2022-02-28 2022-03-25 浙江大学 基于神经网络的风力发电机变桨轴承故障诊断方法及装置
CN116857118A (zh) * 2023-09-01 2023-10-10 江苏纳泉弘远新能源科技有限公司 一种基于物联网技术的风力发电变桨方法和系统
WO2023213366A1 (en) * 2022-05-02 2023-11-09 Vestas Wind Systems A/S Controlling a wind turbine based on wear to wind turbine rotor blade pitch bearings

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105740619A (zh) * 2016-01-28 2016-07-06 华南理工大学 基于核函数的加权极限学习机污水处理在线故障诊断方法
CA2965340A1 (en) * 2016-05-11 2017-11-11 Mahmoud Ismail An enhanced system and method for conducting pca analysis on data signals
CN110146293A (zh) * 2019-06-04 2019-08-20 昆明理工大学 一种基于pca与elm的滚动轴承故障诊断方法
CN111783531A (zh) * 2020-05-27 2020-10-16 福建亿华源能源管理有限公司 一种基于sdae-ielm的水轮机组故障诊断方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105740619A (zh) * 2016-01-28 2016-07-06 华南理工大学 基于核函数的加权极限学习机污水处理在线故障诊断方法
CA2965340A1 (en) * 2016-05-11 2017-11-11 Mahmoud Ismail An enhanced system and method for conducting pca analysis on data signals
CN110146293A (zh) * 2019-06-04 2019-08-20 昆明理工大学 一种基于pca与elm的滚动轴承故障诊断方法
CN111783531A (zh) * 2020-05-27 2020-10-16 福建亿华源能源管理有限公司 一种基于sdae-ielm的水轮机组故障诊断方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
何群 等: "基于相关主成分分析和极限学习机的风电机组主轴承状态监测研究", 《计量学报》 *
叶永伟 等: "基于PCA的汽车涂装线设备信号特征提取", 《仪器仪表学报》 *
孙毅刚 等: "加权在线贯序极限学习机算法及其应用", 《计算机工程与设计》 *
许玉格等: "基于核函数的加权极限学习机污水处理在线故障诊断", 《化工学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113343579A (zh) * 2021-06-25 2021-09-03 海南小鲨鱼智能科技有限公司 一种电池输出功率的控制方法、装置及计算机存储介质
CN113312731A (zh) * 2021-06-28 2021-08-27 北京南洋思源智能科技有限公司 一种变桨轴承故障检测方法及装置、存储介质
CN113326589A (zh) * 2021-06-28 2021-08-31 北京南洋思源智能科技有限公司 一种敏感特征筛选方法及装置、存储介质
CN113898528A (zh) * 2021-09-30 2022-01-07 江苏金风软件技术有限公司 风机变桨轴承的异常检测方法、模型训练方法及相关装置
CN114235415A (zh) * 2022-02-28 2022-03-25 浙江大学 基于神经网络的风力发电机变桨轴承故障诊断方法及装置
CN114235415B (zh) * 2022-02-28 2022-06-03 浙江大学 基于神经网络的风力发电机变桨轴承故障诊断方法及装置
WO2023213366A1 (en) * 2022-05-02 2023-11-09 Vestas Wind Systems A/S Controlling a wind turbine based on wear to wind turbine rotor blade pitch bearings
CN116857118A (zh) * 2023-09-01 2023-10-10 江苏纳泉弘远新能源科技有限公司 一种基于物联网技术的风力发电变桨方法和系统
CN116857118B (zh) * 2023-09-01 2023-11-10 江苏纳泉弘远新能源科技有限公司 一种基于物联网技术的风力发电变桨方法和系统

Similar Documents

Publication Publication Date Title
CN112924173A (zh) 一种风力发电机组变桨轴承故障诊断方法
CN111275288B (zh) 基于XGBoost的多维数据异常检测方法与装置
CN106682814B (zh) 一种基于故障知识库的风电机组故障智能诊断方法
WO2022077901A1 (zh) 一种面向小样本数据集的轴承故障模式诊断方法及系统
CN109782603A (zh) 旋转机械耦合故障的检测方法及监测系统
CN109255333B (zh) 一种大型风电机组滚动轴承故障混合诊断方法
CN111103137A (zh) 基于深度神经网络的风力涡轮机齿轮箱故障诊断方法
CN111273125A (zh) 基于rst-cnn的电力电缆沟道故障诊断方法
CN113738595A (zh) 一种风力发电机组叶片状态监测方法及其系统
CN111198098A (zh) 一种基于人工神经网络的风力发电机轴承故障预测方法
CN110969185A (zh) 基于数据重构的设备异常状态检测方法
CN112651426A (zh) 一种风电机组滚动轴承故障诊断方法
CN111594392A (zh) 一种风力发电塔筒螺栓在线监测方法
CN114326486B (zh) 一种基于概率慢特征分析和弹性权重巩固的过程监控方法
Wang et al. Fault diagnosis of industrial robots based on multi-sensor information fusion and 1D convolutional neural network
CN111307455A (zh) 一种基于字典学习的列车转向架轴承故障监测方法及系统
CN113919102A (zh) 一种风电机组主传动系统载荷传递特性分析方法及系统
CN113240022A (zh) 多尺度单分类卷积网络的风电齿轮箱故障检测方法
CN112986821A (zh) 海上风电机组变桨电机转子断条故障诊断方法
CN117093938A (zh) 一种基于深度学习的风机轴承的故障检测方法和系统
CN114897292B (zh) 基于多尺度残差卷积变分网络故障诊断方法和诊断系统
CN115163426A (zh) 基于ai听诊的风机故障检测方法及系统、风机安全系统
CN115326393A (zh) 一种基于温度信息的风电机组轴承对故障诊断方法
Li et al. Deep learning based end-to-end rolling bearing fault diagnosis
CN116067489B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210608