一种电站风机状态预警方法、系统及其应用与流程

1.本发明涉及火电机组风机的状态监测技术领域，具体涉及一种电站风机状态预警方法、系统及其应用。


背景技术：

2.在当前“双碳”目标的背景下，随着电力市场体制改革的不断深入，传统火电企业面临着日渐加剧的行业竞争，对机组的运行灵活性和调峰能力也提出了更高的要求，保持机组安全、经济和环保的长期稳定运行愈加重要，如何合理地安排机组检修周期，提高各主辅设备利用率成为关键的环节。然而实际中对于电站辅机的状态监测主要包括dcs的限幅报警与现场人员的点检定修，这类方法动态性与实时性较低，属于事后维护。随着运行年限的增加以及机组负荷变化的频率增加，磨损、老化程度增高，严重影响机组的经济性与安全性。
3.目前现有设备监控技术由于故障数据缺乏、检测手段单一等原因，容易导致误报、故障不能识别等问题；由于监控过程数据传输量通常较大，上位机处理负担较重，造成预警的实时性较差。并且深度调峰过程中，设备将出现的一系列安全性问题(由于小流量引起的风机失速、喘振、振动大)，使得现有技术的不足进一步得到放大。因此，必须开发高实时性的辅机状态定量评估技术，完成由“事后维护”到“事前预知”的转变。


技术实现要素：

4.为克服现有技术存在的上述缺陷以满足电站风机安全可靠运行的基本需求，本发明提供了一种电站风机状态预警方法及系统，提高对风机状态的监测效率。首先利用传统及新型检测技术实现了电流、压力、流量、温度、振动、噪声等多维状态参数的在线检测，在前端提出了针对振动、噪声信号的特征提取方法，在终端提出了基于堆叠稀疏降噪自编码网络的故障状态评估及预警算法，并开发了基于边缘计算架构的重要辅机状态监控及预警系统，完成了辅机状态评估及预警系统的功能性验证。
5.为实现上述目的，本发明提供如下技术方案：
6.本发明提供了一种电站风机状态预警方法，具体包括以下步骤：
7.步骤1：采集风机正常运行振动、噪声信号；
8.步骤2：采集风机流量、电流、功率等过程参数时间序列；所述风机过程参数时间序列包括风机运行数据和所述风机运行数据对应的风机运行状态，所述风机状态包括正常运行状态和故障状态；
9.步骤3：采用基于小波包分解方法在所采集的风机正常运行状态下的振动数据时间序列中提取并构建小波能量谱；
10.步骤4：采用基于1/3倍频程分解方法在所采集的风机正常运行状态下的噪声数据时间序列中提取并构建频带能量谱；
11.步骤5：分别以所述过程参数、振动、噪声能量谱数据为输入，通过构建堆叠稀疏降
噪自编码网络进行模型训练，获得过程参数、风机振动、噪声预警模型；
12.步骤6：获取待评估时间段内的过程参数、振动、噪声能量谱数据xobs，分别输入过程参数、风机振动、噪声预警模型，得到模型重构值xrec；
13.步骤7：构建基于模型观测值x
obs
和重构值x
rec
之间欧氏距离和余弦距离的相似度指标；
14.步骤8：采用贝叶斯最大后验概率估计方法，对上述单一预警模型的相似度指标进行融合，构造一个基于多源信息融合的风机状态预警指标，设置状态预警阈值，并根据故障预警阈值判断是否出现故障。
15.本发明公开一种将上述的电站风机状态预警方法应用于电站风机状态预警系统。
16.本发明公开一种种火电机组风机的状态监测方法，其特征为：包括上述的电站风机状态预警方法。
17.有益效果：
18.一、本发明采用先进的声压传感器和加速度传感器，实现高精度的噪声和振动信号采集，弥补了当前电厂设备检测手段不足的缺点。
19.二、本发明采用了小波包分解的振动信号特征提取方法，不仅可以通过频率高低动态调整时间分辨率大小，还可以在提取低频信号的同时获取高频信号。
20.三、本发明选用了发电机功率、电机电流、电机绕组温度、驱动端轴承温度、自由端轴承温度、液压油压力、液压油温度、润滑油压力、润滑油温度、风机动叶开度、入口风温、出口风温和出口风压作为风机状态监测参数，在实际火电机组运行过程中，选择上述监测参数对风机进行监控是十分必要且有效的。
21.四、本发明采用堆叠稀疏降噪自编码网络算法，相比已有的机器学习或深度学习模型，增加了模型的鲁棒性。
22.五、本发明所采用的多源信息融合方法提升系统整体评估的可靠性。
23.六、本发明所设计的风机状态预警系统，能够实现风机状态的在线实时预警，在故障发生前发出预警信号，为现场人员争取处理时间，减少因故障引起的安全与经济性损失。
附图说明
24.图1风机振动原始信号及其频带能量；
25.其中a)风机水平振动原始信号；
26.b)水平振动信号特征提取结果；
27.c)风机垂直振动原始信号；
28.d)垂直振动信号特征提取结果；
29.图2风机噪声原始信号及1/3倍频程特征提取结果；
30.其中：a)为风机噪声原始信号
31.b)为：噪声信号特征提取结果
32.图3为堆叠稀疏降噪自编码网络模型训练流程图；
33.图4为堆叠稀疏降噪自编码网络模型训练误差
34.图5为风机在线状态预警流程图；
35.图6为风机状态预警系统结构示意图。
具体实施方式
36.一种电站风机状态预警方法，具体包括以下步骤：
37.步骤1：采集风机正常运行振动、噪声信号；
38.步骤2：采集风机流量、电流、功率等过程参数时间序列；所述风机过程参数时间序列包括风机运行数据和所述风机运行数据对应的风机运行状态，所述风机状态包括正常运行状态和故障状态；
39.步骤3：采用基于小波包分解方法在所采集的风机正常运行状态下的振动数据时间序列中提取并构建小波能量谱；
40.步骤4：采用基于1/3倍频程分解方法在所采集的风机正常运行状态下的噪声数据时间序列中提取并构建频带能量谱；
41.步骤5：分别以所述过程参数、振动、噪声能量谱数据为输入，通过构建堆叠稀疏降噪自编码网络进行模型训练，获得过程参数、风机振动、噪声预警模型；
42.步骤6：获取待评估时间段内的过程参数、振动、噪声能量谱数据xobs，分别输入过程参数、风机振动、噪声预警模型，得到模型解码输出即为重构值xrec；
43.步骤7：构建基于模型观测值x
obs
和重构值x
rec
之间欧氏距离和余弦距离的相似度指标；
44.步骤8：采用贝叶斯最大后验概率估计方法，对上述单一预警模型的相似度指标进行融合，构造一个基于多源信息融合的风机状态预警指标，设置状态预警阈值，并根据故障预警阈值判断是否出现故障。
45.进一步的，所述步骤1具体包括以下步骤：
46.步骤1.1，在距离风机电机和风机本体直线距离50cm处，以底盘固定式支架各加装一个声传感器，在风机轴承外壳x/y方向上以强力胶粘贴方式各加装一个振动加速度传感器；
47.步骤1.2，以25.6khz的采样频率采集风机两路噪声信号，以51.2khz的采样频率采集两路振动信号。
48.进一步的，步骤2包括以下步骤：
49.步骤2.1，通过与风机连接的分散控制系统或安全仪表系统采集发电机功率、电机电流、电机绕组温度、驱动端轴承温度、自由端轴承温度、液压油压力、液压油温度、润滑油压力、润滑油温度、风机动叶开度、入口风温、出口风温和出口风压数据；
50.步骤2.2，从历史数据中抽取20000组数据，计算各参数之间的皮尔逊相关系数绝对值，对于两两之间相关系数绝对值大于0.6的参数只保留其中一个参数作为建模变量。
51.进一步的，步骤3包括以下步骤：
52.步骤3.1，对采集的振动信号进行频谱分析，从而确定小波包分解时划分合适的频段；
53.步骤3.2，确定小波包分解的层数。在进行小波包分解时，分解的层数与时频域信号的分析精度有显著关联；
54.步骤3.3，对信号进行小波包分解与重构，并求取各子带信号的能量，得到重构后的信号能量谱；
55.步骤3.4，对重构后的信号能量谱进行划分得到相邻且相等的频段，计算出各个频
段的归一化信号能量比。
56.所述步骤3.3中，小波包函数表达式为：
[0057][0058]
式中：其中，l为振荡系数，k为位置指标，hk和gk分别为尺度函数导出的高通和低通正交滤波器系数；选用daubechies小波对振动信号进行三阶小波包分解。
[0059]
频带能量的计算公式为：
[0060][0061]
式中：ej为子带信号第j个节点的能量；sj(t)为原始信号，xj为原始信号离散点幅值，m/s2；n为采样点数。
[0062]
所述步骤3.4中，将经过三层小波包分解得到的第三阶4个低频分量与4个高频分量组成一个列为8的特征向量t，
[0063][0064]
式中，t为单个样本归一化后的特征向量；ej为第三层第j个重构信号。对特征向量进行归一化处理，每个向量为能量的百分比，便于去掉量纲。
[0065]
进一步的，步骤4包括以下步骤：
[0066]
步骤4.1，对采集的噪声信号进行快速傅里叶变换；
[0067]
步骤4.2，根据功率相等恢复系数添加汉宁窗得到噪声的初步频谱图；
[0068]
步骤4.3，提取中心频率标称值和1/3倍频程范围，并对1/3倍频程运用振动振级的方法将图形转变成未计权1/3倍频程声压级频谱图；
[0069]
步骤4.4，经过c计权声级修正，得到噪声的c计权1/3倍频程声压级频谱图。
[0070]
所述步骤4.2中，汉宁窗窗长归一化后的时域表达式为：
[0071][0072]
频域表达式为：
[0073][0074]
为使加窗后的功率谱或幅值谱不受窗函数的影响，必须根据一定原则推导出恢复函数，恢复系数计算公式为：
[0075]kt
＝a/a'a为谐波信号的幅值；a'为加窗后峰值频率处的幅值。
[0076]
可直接对窗函数作频谱分析，a'使窗函数幅值谱在零频率处的幅值，此时a等于1。根据此种方法计算得到汉宁窗的幅值恢复系数为2，所以功率相等恢复系数为1.633。
[0077]
所述步骤4.3中，声压级的定义如下：
[0078][0079]
式中:l
p
为声压级，db；p0＝2
×
10
－5
pa，为基准声压，基准声压本身的声压级为0db；pe为有效声压。
[0080]
所述步骤4.3中，噪声信号的频率提取范围为20hz-20khz，倍频程的各频带具有对应的上限频率与下限频率比值设定为2
1/3
，即：
[0081][0082]
其中fu为上限频率，fd为下限频率。
[0083]
所述步骤4.4中，c计权的计算公式为：
[0084][0085]
进一步的，步骤5包括以下步骤：
[0086]
步骤5.1，将历史数据作归一化处理，用以防止数据数量级的不同造成量级较大的属性占主导地位并缓解数量级的不同导致的迭代收敛速度变慢；
[0087]
步骤5.2，将已归一化数据划分成训练集数据和测试集数据；
[0088]
步骤5.3，对输入数据x进行增加噪声处理得新的输入生成堆叠降噪自编码器；
[0089]
步骤5.4，确定初始ssdae网络的隐含层数、各层的神经元数及激活函数；
[0090]
步骤5.5，堆叠降噪自编码器损失函数中加入约束生成深度稀疏降噪自编码器，即堆叠稀疏降噪自编码网络；
[0091]
步骤5.6，采用优化器对堆叠稀疏降噪自编码网络参数进行优化；
[0092]
步骤5.7，保存训练最优模型(即预警模型)及网络参数。
[0093]
步骤5.8，重复步骤5.1-5.7，分别以过程参数、风机振动、噪声为输入训练并保存预警模型。
[0094]
所述步骤5.1中，归一化处理采用minmax归一化处理方法。
[0095]
所述步骤5.2具体为：将归一化后的历史数据分为相同长度的多组数据形成矩阵，作为训练集和测试集的全部数据，采用留出法将该矩阵中的3/4组数据作为训练集，剩余1/4组数据作为测试集，并且训练集以及测试集数据为随机排列。
[0096]
所述步骤5.3中，数据增加噪声的方法为：
[0097]
步骤5.3.1，生成与x同维度的均匀分布噪声noise；
[0098]
步骤5.3.2，将x与noise做哈达玛乘积，具体公式为：
[0099][0100]
可选的，增加噪声还可采用椒盐噪声(salt-and-pepper noise)、高斯噪声(gaussian noise)、遮罩噪声(masking noise)、白噪声(white noise)。
[0101]
所述步骤5.4中，输出层采用sigmoid激活函数，其余层均采用leaky relu激活函数。
[0102]
所述步骤5.5中，堆叠稀疏降噪自编码器的优化目标为：
[0103][0104][0105]
式中，为隐藏层神经元的激活值；β为控制稀疏性惩罚因子的权重；x为预处理后的运行数据,为加噪后的运行数据；是分别以ρ和为均值的两个伯努利随机变量之间的相对熵，当时，kl随着ρ和之间的差异增加而单调递增，而当靠近0或者1的时候，相对熵则变得无穷大，因此最小化该惩罚因子，可以使趋近于ρ。
[0106]
所述步骤5.6中，优化器采用adam优化算法。
[0107]
进一步的，步骤7包括以下步骤：
[0108]
步骤7.1，分别计算观测值x
obs
和重构值x
rec
之间欧氏距离和余弦距离；
[0109]
步骤7.2，计算基于欧氏距离和余弦距离的模型相似度。
[0110]
所述步骤7.1中，欧氏距离的计算公式为：
[0111][0112]
余弦距离的计算公式为：
[0113][0114]
式中dot(
·
)表示内积运算，||
·
||表示向量的2范数。
[0115]
所述步骤7.2中，模型相似度的计算公式为：
[0116][0117]
式中，sim(x,y)表示向量x与y的相似度。本发明中构造的余弦距离值介于0与1之间，且越接近于0代表重构值与输入值的形态相似程度越低。
[0118]
所述步骤8中，贝叶斯最大后验概率估计计算公式为：
[0119][0120]
其中，xf表示风机状态指标，zi表示第i个信息源的相似度指标，σi表示第i个信息源的标准差。
[0121]
所述的预警阈值取值为0.8，当风机状态指标xf低于该阈值时，系统即发出风机状态异常警告。
[0122]
与现有技术相比，本发明具有以下优点：
[0123]
一、本发明采用先进的声压传感器和加速度传感器，实现高精度的噪声和振动信号采集，弥补了当前电厂设备检测手段不足的缺点。
[0124]
二、本发明采用了小波包分解的振动信号特征提取方法,不仅可以通过频率高低动态调整时间分辨率大小，还可以在提取低频信号的同时获取高频信号。有效弥补了传统信号分析方法的缺点，如短时傅里叶变换在不同频域内只有固定的分辨率；经验模态分解缺乏理论证明；小波变换只会提取低频信息，而高频细节信息在解析过程中易被丢失。所采用的1/3倍频程噪声信号特征提取方法，通过对噪声信号进行加窗快速傅里叶变换，成功避免了对信号截断处理而引起的频谱泄露和能量泄露，同时采用c计权算法对频带声压级进行修正，还原出人耳对声音的听觉特性。同时通过在设备边缘部署采集前端进行实时计算，一方面减少了数据传输量，减轻网络压力；另一方面将计算过程前移能有效缓解上位机的计算压力。
[0125]
三、本发明选用了发电机功率、电机电流、电机绕组温度、驱动端轴承温度、自由端轴承温度、液压油压力、液压油温度、润滑油压力、润滑油温度、风机动叶开度、入口风温、出口风温和出口风压作为风机状态监测参数，由于这些监测参数是风机监控系统中比较容易出现异常的现场运行参数，通过对这些参数的监控可以很好地判断风机的运行情况，比如：风机轴承的故障或润滑油质变差会使风机轴承温度、油温上升；叶片的磨损、转子不对中或旋转失速会导致风机异常振动。因此，选择上述监测参数对风机进行监控是十分必要且有效的。
[0126]
四、本发明采用堆叠稀疏降噪自编码网络算法，相比已有的机器学习或深度学习模型，采用深度稀疏降噪自编码器能有效提高模型的泛化能力；由于原始训练集数据的某一部分被破坏掉了，本发明能够在一定程度上减轻训练集和测试集之间的代沟，因此其在一定程度上更接近测试数据，从而增加了模型的鲁棒性。
[0127]
五、本发明所采用的多源信息融合方法，充分利用多种信息源，对多模型的输出相似度进行融合，提升系统整体评估的可靠性。以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是本发明的原理，在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明的范围内。本发明要求的保护范围由所附的权利要求书及其等同物界定。