一种基于分散式ESN模型的风力发电机故障检测方法技术

本发明专利技术公开一种基于分散式ESN模型的风力发电机故障检测方法，通过利用ESN来建立风力发电机各个测量变量之间的量化模型，从而利用数据驱动模型的误差来反映风力发电机运行过程中出现的故障。具体来讲，本发明专利技术方法通过结合使用ESN与分散式建模策略的优势，为风力发电机各个测量变量建立分散式的ESN模型，再通过监测分散式ESN模型的估计误差来实现风力发电机故障检测目的。一方面，本发明专利技术方法利用ESN模型量化描述了风力发电机各测量变量之间体现在时序特征上的非线性关系。另一方面，本发明专利技术方法通过监测分散式ESN模型的估计误差的异常变化情况来反映是否出现故障，借鉴参考了利用机理模型生成误差的思想的优势。用机理模型生成误差的思想的优势。用机理模型生成误差的思想的优势。

【技术实现步骤摘要】
一种基于分散式ESN模型的风力发电机故障检测方法


[0001]本专利技术涉及一种风力发电机故障检测方法，特别涉及一种基于分散式ESN模型的风力发电机故障检测方法。

技术介绍

[0002]随着自然资源的日益短缺，充分利用可再生能源领域受到了越来越多的关注。其中，风力发电占据着举足轻重的地位。在风力持续的情况下，风力发电机能持续不断的提供电能，因此风力发电机的应用规模与日俱增。与此同时，及时的检测出风力发电机运行过程中出现的故障对于保证风力发电高效的电力输出具有重要的研究意义。由于风力发电机是一个大型且复杂的系统，通过物理原理建立相应的动力学模型和电路模型是非常困难的。因此，利用机理模型实施故障检测的方案是不可取的。
[0003]一般而言，风力发电机这个系统中通常会安装有多个传感器，会实时反馈诸如发电机转速，生成的电功率，加速度等数据信息。这些传感器以及相应配套的数据存储设备为实施数据驱动的风力发电机故障检测提供了坚实的数据基础。在当前智能制造与大数据的风潮下，利用这些采样数据实施风力发电机故障检测的方案是非常合时宜的，已有多种按照模式分类思想实施风力发电机故障检测与诊断的数据驱动的技术。然而，风力发电机的工作状态会受到外部环境风速的直接影响，会随着风速的变化而不断发生变化的。因为风力的间歇特性、非线性特性、时序变化特性并非人为可精准预测或可控制的，风力发电机的工作状态直接受风速影响的这种工作特性给数据驱动的故障检测方案增加了挑战。
[0004]从风力发电机运行机理上讲，每个时刻的风力发电机自身数据都是会受到外部环境风速的直接影响，而且风力发电机自身的测量数据之间也是相互关联的。从这个角度上讲，如何挖掘这些测量变量之间的关联性，对于实施数据驱动的风力发电机故障检测具有重要的作用。假设能建立精确的风力发电机的机理模型，则该模型就能描述出环境风速与风力发电机自身测量变量(如：风叶转速，功率等)之间的相关关系，相应的机理模型误差的异常变化即可反映出风力发电机运行出现故障问题。由于精确的机理模型无法建立，可行的解决思路只能通过数据驱动的方法来建立这些测量变量之间的关系模型，从而使用数据驱动模型的误差实现故障检测。
[0005]近年来，在时间序列分析相关的研究领域，回声状态网络(Echo State Networks，缩写：ESN)得到了广泛的研究与应用。然而，ESN是一类有监督的建模策略，需要多故障类别的风力发电机数据才可实施故障诊断，无法直接应用ESN解决只有正常工况数据下的风力发电机故障检测问题。虽然，ESN从表面上看应用于风力发电机故障检测看似困难重重，但从ESN在解决时序建模问题的优势来看，是可以用来解决风速的时序变化特性建模问题。从这个角度出发，如何将ESN用于风力发电机故障检测是值得研究并丞待解决的一个问题。

技术实现思路

[0006]本专利技术所要解决的主要技术问题是：如何利用ESN来建立风力发电机各个测量变
量之间的量化模型，从而利用数据驱动模型的误差来反映风力发电机运行过程中出现的故障。具体来讲，本专利技术方法通过结合使用ESN与分散式建模策略的优势，为风力发电机各个测量变量建立分散式的ESN模型，再通过监测分散式ESN模型的估计误差的异常变化来实现风力发电机故障的实时检测目的。
[0007]本专利技术方法解决上述问题所采用的技术方案为：一种基于分散式ESN模型的风力发电机故障检测方法，包括以下所示步骤：
[0008]步骤(1)：在风力发电机正常运行状态下，每间隔1分钟采集一次样本数据，具体包括11个数据，并将每次采集的这11个数据组成一个列向量；其中，列向量中的11个数据依次为：风速，旋翼转速，发电机转速，机械转矩，发电功率，叶片螺距角，叶片方位角，叶片根部力矩，顶部水平轴加速度，顶部纵轴加速度，偏航误差。
[0009]步骤(2)：按照采集的时间先后顺序，将风速在5m/s至20m/s之间的N个列向量x1，x2，
…
，x
N
组成矩阵X＝[x1，x2，
…
，x
N
]，并对X∈R
11
×
N
中各个行向量实施归一化处理从而得到新矩阵其中，R
11
×
N
表示11
×
N维的实数矩阵，R
11
×1表示11
×
1维的实数向量，归一化处理的具体方式为：将行向量减去该行向量的最小值后，再除以该行向量的最大值与最小值之差。
[0010]步骤(3)：依次分别将中第2行至第11行的行向量用作ESN模型的输出，并将中其它10行的行向量组成输入矩阵后用作ESN模型的输入，建立由10个ESN模型组成的分散式ESN模型，并保留模型参数集，具体的实施过程如步骤(3.1)至步骤(3.9)所示。
[0011]步骤(3.1)：初始化j＝2。
[0012]步骤(3.2)：将新矩阵中第j行的行向量记做为输出向量z
j
∈R1×
N
，并将新矩阵中其余11行的行向量组成输入矩阵Z
j
∈R
10
×
N
；其中，R1×
N
表示1
×
N维的实数向量，R
10
×
N
表示10
×
N维的实数矩阵。
[0013]步骤(3.3)：搭建一个回声状态网络模型，具体包括：8个输入节点，L个储蓄池节点，和1个输出节点；其中，储蓄池节点的激活函数f(u)＝tanh(u)为双曲正切函数，u为函数自变量。
[0014]步骤(3.4)：初始化权值连接矩阵W
j
∈R
L
×
L
；其中，W
j
的各个元素都是按照均匀分布在区间[-1，1]范围内随机产生的。
[0015]步骤(3.5)：先对权值连接矩阵W
j
进行稀疏度修正，再对W
j
进行谱密度修正。
[0016]所谓的稀疏度修正，就是将W
j
中小于等于0.1且大于等于-0.1的非对角线上的元素修改成0；实施稀疏度修正时，不修改W
j
中对角线上的元素，只对除对角线以外的其他元素进行修正。
[0017]而谱密度修正需要首先计算权值连接矩阵W
j
的所有特征值，并对所有特征值计算绝对值后，再将最大的绝对值记做α；其次，判断α是否小于1；若否，则根据公式W
j
＝0.8
·
W
j
/α更新权值连接矩阵W
j
；若是，则不对权值连接矩阵W
j
进行任何处理。
[0018]步骤(3.6)：初始化输入权值矩阵初始化反馈权值向量V
j
∈R
L
×1；其中，输入权值矩阵和反馈权值向量中的元素都是按照均匀分布在区间[-1，1]范围内随机本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于分散式ESN模型的风力发电机故障检测方法，其特征在于，包括以下步骤：步骤(1)：在风力发电机正常运行状态下，每间隔1分钟采集一次样本数据，具体包括11个数据，并将每次采集的这11个数据组成一个列向量；其中，列向量中的11个数据依次为：风速，旋翼转速，发电机转速，机械转矩，发电功率，叶片螺距角，叶片方位角，叶片根部力矩，顶部水平轴加速度，顶部纵轴加速度，偏航误差；步骤(2)：按照采集的时间先后顺序，将风速在5m/s至20m/s之间的N个列向量x1，x2，
…
，x
N
组成矩阵X＝[x1，x2，
…
，x
N
]，并对X∈R
11
×
N
中各个行向量实施归一化处理从而得到新矩阵其中，R
11
×
N
表示11
×
N维的实数矩阵，R
11
×1表示11
×
1维的实数向量，归一化处理的具体方式为：将行向量减去该行向量的最小值后，再除以该行向量的最大值与最小值之差；步骤(3)：依次分别将中第2行至第11行的行向量用作ESN模型的输出，并将中其它10行的行向量组成输入矩阵后用作ESN模型的输入，建立由10个ESN模型组成的分散式ESN模型，并保留模型参数集，具体的实施过程如步骤(3.1)至步骤(3.9)所示：步骤(3.1)：初始化j＝2；步骤(3.2)：将新矩阵中第j行的行向量记做为输出向量z
j
∈R1×
N
，并将新矩阵中其余10行的行向量组成输入矩阵Z
j
∈R
10
×
N
；其中，R1×
N
表示1
×
N维的实数向量，R
10
×
N
表示10
×
N维的实数矩阵；步骤(3.3)：搭建一个回声状态网络模型，具体包括：8个输入节点，L个储蓄池节点，和1个输出节点；其中，储蓄池节点的激活函数f(u)＝tanh(u)为双曲正切函数，u为函数自变量；步骤(3.4)：初始化权值连接矩阵W
j
∈R
L
×
L
；其中，W
j
的各个元素都是按照均匀分布在区间[-1，1]范围内随机产生的；步骤(3.5)：先对权值连接矩阵W
j
进行稀疏度修正，再对W
j
进行谱密度修正；步骤(3.6)：初始化输入权值矩阵初始化反馈权值向量V
j
∈R
L
×1；其中，输入权值矩阵和反馈权值向量V
j
中的元素都是按照均匀分布在区间[-1，1]范围内随机产生的；步骤(3.7)：根据公式计算输出权值向量其中，表示z
j
中第2个至第N个元素组成的向量，上标号T表示矩阵或向量的转置符号，λ为正则系数，I为K
×
K维的单位矩阵，n＝N-1，K＝11+L，矩阵Z
j
的组成方式如下所示：上式中，表示X
j
中第2列至第N列的列向量组成的矩阵，表示z
j
中第1个至第N-1个元素组成的向量，表示状态矩阵S
j
∈R
L
×
N
中第2列至第N列的列向量组成的矩阵，状态矩阵S
j
中第一列向量为零向量，S
j
中第i列向量的计算方式如下所示：
上式中，表示X
j
中的第i列向量，i∈{2，3，
…
，N}，表示状态矩阵S
j
中第i-1列向量，表示z
j
中的第i-1个元素；步骤(3.8)：保留模型参数集并根据公式计算输出估计向量y
j
∈R1×
n
后，再计算估计误差向量步骤(3.9)：判断是否满足条件：j＜11；若是，则设置j＝j+1后，返回步骤(3.2)；若否，则得到模型参数集φ2，φ4，
…
，φ
11
，以及误差向量e2，e4，
…
，e
11
；步骤(4)：将10个误差向量e2，e3，
…
，e
11
组成误差矩阵后，再计算误差矩阵E中所有行向量的平均值向量
μ
∈R1×
10
，并根据公式Λ＝(E-U)
T
(E-U)/(N-1)计算协方差矩阵Λ；其中，均值矩阵U∈R...

【专利技术属性】
技术研发人员：ꢀ五一IntClF零三D一七零零，
申请(专利权)人：宁波大学，
类型：发明
国别省市：