基于模糊分区的风电机组变桨系统故障预警识别方法，基于梯形模糊分区理论，将风力等级分为四个模糊分区，也称模糊集合F，并计算模糊集合F的隶属度函数；将风电场SCADA得到的风速数据按照时间段进行风速等级模糊分区处理，确定各模糊分区处理，确定各模糊分区范围，通过步骤1所述的隶属度函数，得到各个分区的隶属度值，用ε1表示；进行长短型记忆网络深度学习算法分区预测，得到各个分区的预测值，用LSTM_out1表示，然后将各个分区的预测值LSTM_out1乘以该分区的隶属度值ε1，再将所有分区εi*LSTM_outi的值相加，得到本风机启动阶段的变桨轴承轴速预测值。本发明克服了现有技术的不足，能够实时预警变桨机械故障隐患状况，实现提前对设备故障隐患做出相应分析判断。
1.基于模糊分区的风电机组变桨系统故障预警识别方法，其特征在于，包括以下步骤：步骤1)基于梯形模糊分区理论，将风力等级分为小风、中风、大风，超大风四个模糊分区，也称模糊集合F，并计算模糊集合F的隶属度函数；
步骤2)将风电场SCADA得到的风速数据按照时间段进行风速等级模糊分区处理，确定各模糊分区处理，确定各模糊分区范围，通过步骤1所述的隶属度函数，得到各个分区的隶属度值，用ε1表示；
步骤3)进行长短型记忆网络深度学习算法分区预测，得到各个分区的预测值，用LSTM_out1表示，然后将各个分区的预测值LSTM_out1乘以该分区的隶属度值ε1，再将所有分区εi*LSTM_outi的值相加，得到本风机启动阶段的变桨轴承轴速预测值；
步骤1)中，所述模糊集合F是指对于任意的x∈X，都有一个确定的数εF(x)∈[0,1]与其对应，εF(x)表示x相对与F的隶属程度，映射：F的隶属度函数，所述隶属度函数选取梯形隶属度函数：其中a,b,c,d是梯形的4个顶点位置；
步骤21)确定无风模糊分区的范围为0～4m/s、小风模糊分区范围为2～8m/s、中风模糊分区范围为6～12m/s、大风模糊分区范围为11～16m/s、超大风模糊分区范围为14m/s以上；
步骤22)选取风机启动阶段连续时间内的风速数据V，取该时间段内的风速的算术平均值为该时间段的风速值，将该时间段的风速值V作为模糊分区的依据，即步骤1)所述的公式(1‑1)中的x值；各模糊分区的隶属度函数对应于公式(1‑2)(1‑3)中，a为该模糊分区的最小风速，b为该模糊分区的风速范围的期望值，c为该模糊分区的中位数，d为该模糊分区的最大风速，统计风机启动阶段这段时间内每秒钟的风速值，分别计算这段时间内风速序列的期望值b和中位数c,比较b和c的大小，如果b
2)，如果b≥c,风速模糊分区的隶属度函数选择公式(1‑3)；
有四层神经网络层以特殊的方式相互作用，而不是单个神经网络层；
步骤31)利用t‑1时刻风机的m个运行状态数据预测t时刻风机变桨轴速的预测值；输入t‑1时刻变桨系统的N组数据，得到的m个主成分量Xm(t‑1)，作为t‑1时刻网络输入层参数，得到t‑1时刻网络的变桨轴速输出值ht‑1，以及t‑1时刻的网络记忆状态Ct‑1，同时输入t时刻的m个风机主成分量Xm(t)，即ht‑1、Ct‑1、Xm(t)作为t时刻输入参数，预测计算得到t时刻网络的变桨轴速输出值ht 和当前时刻记忆状态C t ，以此类推；
步骤32)隐含层包含n个按前后时刻连接的同构LSTM单元，Ct‑1、ht‑1为前一个LSTM单元的状态和输出，计算公式如下所示：
ft＝σ(Wf·[ht‑1，xt]+bf) (2‑1)it＝σ(Wi·[ht‑1，xt]+bi) (2‑2)其中σ为sigmoid函数；
ot＝σ(Wo[ht‑1，xt]+bo) (2‑5)ht＝ot*tanh(Ct) (2‑6)式中，ot是输出门函数，Ct是单元的记忆值，输出层的输出值即为输出门函数与单元的记忆值Ct的正切函数的乘积；
步骤34)经模糊分区和长短型记忆网络预测后的变桨轴速预测值为：y＝εA·LSTM_outB+εB·LSTM_outA (2‑7)式中，εA、εB表示该时间段风速隶属于相应两个模糊分区的隶属度值，LSTM_outA、LSTM_outB表示在相应两个模糊分区下的长短型记忆网络预测模型的输出值。
基于模糊分区的风电机组变桨系统故障预警识别方法
技术领域
[0001] 本发明涉及风电机组故障预警方法领域，具体涉及基于模糊分区的风电机组变桨系统故障预警识别方法。
背景技术
[0002] 变桨直流电机的转矩只和电机电流大小有关，理想情况下，在叶片从90变到45的过程中：小风速阶段(3‑8)m/s变桨电机变桨速度是2°/s，变桨电机转矩是恒定值，变桨电机
电流也是恒定值,在中风速和大风速阶段，风机爬坡时，变桨速度≤2°/s，变桨电机的转矩
和电流也随之发生变化。变桨直流电机输入是电流，输出是扭矩，变桨速度与变桨受力(扭
矩)以及变桨自身健康状况有关，同样的力施加在健康的轴承和内部损坏的轴承，他们的变
桨速度是不一样的。
发明内容
[0003] 针对现有技术的不足，本发明提供了基于模糊分区的风电机组变桨系统故障预警识别方法，克服了现有技术的不足，能够实时预警变桨机械故障隐患状况，实现提前对设备
故障隐患做出相应分析判断。
[0004] 为实现以上目的，本发明通过以下技术方案予以实现：
[0005] 本发明公开了一种基于模糊分区的风电机组变桨系统故障预警识别方法，包括以下步骤：
[0006] 步骤1)基于梯形模糊分区理论，将风力等级分为小风、中风、大风，超大风四个模糊分区，也称模糊集合F，并计算模糊集合F的隶属度函数；
[0007] 步骤2)将风电场SCADA得到的风速数据按照时间段进行风速等级模糊分区处理，确定各模糊分区处理，确定各模糊分区范围，通过步骤1所述的隶属度函数，得到各个分区
的隶属度值，用ε1表示；
[0008] 步骤3)进行长短型记忆网络深度学习算法分区预测，得到各个分区的预测值，用LSTM_out1表示，然后将各个分区的预测值LSTM_out1乘以该分区的隶属度值ε1，再将所有
分区εi*LSTM_outi的值相加，得到本风机启动阶段的变桨轴承轴速预测值。
[0009] 优选地，所述步骤1)中，所述模糊集合F是指对于任意的x∈X，都有一个确定的数εF(x)∈[0,1]与其对应，εF(x)表示x相对与F的隶属程度，映射：
[0010]
[0011] F的隶属度函数，所述隶属度函数选取梯形形隶属度函数：
[0012]
[0013]
[0014] 其中a,b,c,d是梯形的4个顶点位置。
[0015] 优选地，所述步骤2)包括如下步骤：
[0016] 步骤21)确定无风模糊分区的范围为0～4m/s、小风模糊分区范围为2～8m/s、中风模糊分区范围为6～12m/s、大风模糊分区范围为11～16m/s、超大风模糊分区范围为14m/s
以上；
[0017] 步骤22)选取风机启动阶段连续时间内的风速数据V，取该时间段内的风速的算术平均值为该时间段的风速值，将该时间段的风速值V作为模糊分区的依据，即步骤1)所述的
公式(1‑1)中的x值；各模糊分区的隶属度函数对应于公式(1‑2)(1‑3)中，a为该模糊分区的
最小风速，b为该模糊分区的风速范围的期望值，c为该模糊分区的中位数，d为该模糊分区
的最大风速，统计风机启动阶段这段时间内每秒钟的风速值，分别计算这段时间内风速序
列的期望值b和中位数c,比较b和c的大小，如果b(1‑2)，如果b≥c,风速模糊分区的隶属度函数选择公式(1‑3)。
[0018] 优选地，所述步骤3包括如下步骤：
[0019] 有四层神经网络层以特殊的方式相互作用，而不是单个神经网络层；
[0020] 步骤31)利用t‑1时刻风机的m个运行状态数据预测t时刻风机变桨轴速的预测值。输入t‑1时刻变桨系统的N组数据，得到的m个主成分量Xm(t‑1)，作为t‑1时刻网络输入层参
数，得到t‑1时刻网络的变桨轴速输出值h(t‑1)，以及t‑1时刻的网络记忆状态c(t‑1)，同时
输入t时刻的m个风机主成分量Xm(t)，即h(t‑1)、c(t‑1)、Xm(t)作为t时刻输入参数，预测计
算得到t时刻网络的变桨轴速输出值h(t)和当前时刻记忆状态c(t)，以此类推；
[0021] 步骤32)隐含层包含n个按前后时刻连接的同构LSTM单元，ct‑1、ht‑1为前一个LSTM单元的状态和输出，计算公式为下面所示：
[0022] ft＝σ(Wf·[ht‑1，xt]+bf) (2‑1)
[0023] it＝σ(Wi·[ht‑1，xt]+bi) (2‑2)
[0024]
[0025]
[0026] 其中σ为sigmoid函数；
[0027] 步骤33)输出层的计算：
[0028] Ot＝σ(Wo[ht‑1，xt]+bo) (2‑5)
[0029] ht＝ot*tanh(Ct) (2‑6)
[0030] 式中，ot是输出门函数，Ct是单元的记忆值，输出层的输出值即为输出门函数与单元的记忆值Ct的正切函数的乘积。
[0031] 步骤34)经模糊分区和长短型记忆网络预测后的变桨轴速预测值为：
[0032] y＝εA·LSTM_outB+εB·LSTM_outA (2‑7)
[0033] 式中，εA、εB表示该时间段风速隶属于相应两个模糊分区的隶属度值，LSTM_outA、LSTM_outB表示在相应两个模糊分区下的长短型记忆网络预测模型的输出值。
[0034] 本发明提供了基于模糊分区的风电机组变桨系统故障预警识别方法。具备以下有益效果：通过分析风速变化的重要特征，以及风速与变桨轴速、变桨电机电流、变桨电机温
度之间的关系，采用模糊理论、深度学习智能算法相结合的模式。在进行风机变桨系统机械
部分故障预警时，先将风速按照时间段进行风力等级模糊分区处理，然后进行长短型记忆
网络分区预测变桨轴速，同时变桨轴速由梯形模糊数转化为区间值，分区变桨轴速预测值
乘以其隶属度值再将所有分区值相加即时是预测的变桨轴速值，变桨轴速与实际变桨轴速
相比，得到风机变桨轴承和齿轮箱的健康度，给出预警诊断结果。
附图说明
[0035] 为了更清楚地说明本发明或现有技术中的技术方案，下面将对现有技术描述中所需要使用的附图作简单地介绍。
[0036] 图1本发明通过长短记忆神经网络得到的变桨轴速预测结果图；
具体实施方式
[0037] 为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述。
[0038] 本发明公开了一种基于模糊分区的风电机组变桨系统故障预警识别方法，包括以下步骤：
[0039] 步骤1)基于梯形模糊分区理论，将风力等级分为小风、中风、大风，超大风四个模糊分区，也称模糊集合F，并计算模糊集合F的隶属度函数；
[0040] 步骤2)将风电场SCADA得到的风速数据按照时间段进行风速等级模糊分区处理，确定各模糊分区处理，确定各模糊分区范围，通过步骤1所述的隶属度函数，得到各个分区
的隶属度值，用ε1表示；
[0041] 步骤3)进行长短型记忆网络深度学习算法分区预测，得到各个分区的预测值，用LSTM_out1表示，然后将各个分区的预测值LSTM_out1乘以该分区的隶属度值ε1，再将所有
分区εi*LSTM_outi的值相加，得到本风机启动阶段的变桨轴承轴速预测值。
[0042] 优选地，所述步骤1)中，所述模糊集合F是指对于任意的x∈X，都有一个确定的数εF(x)∈[0,1]与其对应，εF(x)表示x相对与F的隶属程度，映射：
[0043]
[0044] F的隶属度函数，所述隶属度函数选取梯形形隶属度函数：
[0045]
[0046]
[0047] 其中a,b,c,d是梯形的4个顶点位置。
[0048] 优选地，所述步骤2)包括如下步骤：
[0049] 步骤21)确定无风模糊分区的范围为0～4m/s、小风模糊分区范围为2～8m/s、中风模糊分区范围为6～12m/s、大风模糊分区范围为11～16m/s、超大风模糊分区范围为14m/s
以上；
[0050] 步骤22)选取风机启动阶段连续时间内的风速数据V，取该时间段内的风速的算术平均值为该时间段的风速值，将该时间段的风速值V作为模糊分区的依据，即步骤1)所述的
公式(1‑1)中的x值；各模糊分区的隶属度函数对应于公式(1‑2)(1‑3)中，a为该模糊分区的
最小风速，b为该模糊分区的风速范围的期望值，c为该模糊分区的中位数，d为该模糊分区
的最大风速，统计风机启动阶段这段时间内每秒钟的风速值，分别计算这段时间内风速序
列的期望值b和中位数c,比较b和c的大小，如果b(1‑2)，如果b≥c,风速模糊分区的隶属度函数选择公式(1‑3)。
[0051] 优选地，所述步骤3包括如下步骤：
[0052] 有四层神经网络层以特殊的方式相互作用，而不是单个神经网络层；
[0053] 步骤31)利用t‑1时刻风机的m个运行状态数据预测t时刻风机变桨轴速的预测值。输入t‑1时刻变桨系统的N组数据，得到的m个主成分量Xm(t‑1)，作为t‑1时刻网络输入层参
数，得到t‑1时刻网络的变桨轴速输出值h(t‑1)，以及t‑1时刻的网络记忆状态c(t‑1)，同时
输入t时刻的m个风机主成分量Xm(t)，即h(t‑1)、c(t‑1)、Xm(t)作为t时刻输入参数，预测计
算得到t时刻网络的变桨轴速输出值h(t)和当前时刻记忆状态c(t)，以此类推；
[0054] 步骤32)隐含层包含n个按前后时刻连接的同构LSTM单元，ct‑1、ht‑1为前一个LSTM单元的状态和输出，计算公式为下面所示：
[0055] ft＝σ(Wf·[ht‑1，xt]+bf) (2‑1)
[0056] it＝σ(Wi·[ht‑1，xt]+bi) (2‑2)
[0057]
[0058]
[0059] 其中σ为sigmoid函数；
[0060] 步骤33)输出层的计算：
[0061] ot＝σ(Wo[ht‑1，xt]+bo) (2‑5)
[0062] ht＝ot*tanh(Ct) (2‑6)
[0063] 式中，ot是输出门函数，Ct是单元的记忆值，输出层的输出值即为输出门函数与单元的记忆值Ct的正切函数的乘积。
[0064] 步骤34)经模糊分区和长短型记忆网络预测后的变桨轴速预测值为：
[0065] y＝εA·LSTM_outB+εB·LSTM_outA (2‑7)
[0066] 式中，εA、εB表示该时间段风速隶属于相应两个模糊分区的隶属度值，LSTM_outA、LSTM_outB表示在相应两个模糊分区下的长短型记忆网络预测模型的输出值。
[0067] 图1所示为本发明经过长短记忆神经网络得到的变桨轴速预测结果。
[0068] 以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施
例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者
替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
