CN115758290A - 一种基于lstm的风机齿轮箱高速轴温度趋势预警方法 - Google Patents

一种基于lstm的风机齿轮箱高速轴温度趋势预警方法 Download PDF

Info

Publication number
CN115758290A
CN115758290A CN202211237677.8A CN202211237677A CN115758290A CN 115758290 A CN115758290 A CN 115758290A CN 202211237677 A CN202211237677 A CN 202211237677A CN 115758290 A CN115758290 A CN 115758290A
Authority
CN
China
Prior art keywords
speed shaft
gearbox
data
temperature
lstm
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
CN202211237677.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.)
Chongqing Datang International Pengshui Hydropower Development Co ltd
Original Assignee
Chongqing Datang International Pengshui Hydropower Development Co ltd
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 Chongqing Datang International Pengshui Hydropower Development Co ltd filed Critical Chongqing Datang International Pengshui Hydropower Development Co ltd
Priority to CN202211237677.8A priority Critical patent/CN115758290A/zh
Publication of CN115758290A publication Critical patent/CN115758290A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明涉及一种基于LSTM的风机齿轮箱高速轴温度趋势预警方法,包括:获取风力发电机齿轮箱高速轴温度同一时间序列的其他相关参数历史数据,对数据预处理;进行皮尔逊相关性分析筛选出与齿轮箱高速轴温度相关性强的参数作为样本数据集;提取训练数据集与测试数据集,搭建LSTM神经网络;配置LSTM神经网络训练参数,对网络进行训练,使用LSTM神经网络高速轴温度趋势预测模型进行预测,通过均方根误差衡量模型的预测误差;计算高速轴预测值与实际值的残差,求得残差均值,设置预警逻辑得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。本发明能够解决现有齿轮箱高速轴温度监测手段无法实现高精度提前预测和温度趋势预警问题。

Description

一种基于LSTM的风机齿轮箱高速轴温度趋势预警方法
技术领域
本发明属于风力发电技术领域,尤其涉及一种基于LSTM的风机齿轮箱高速轴温度趋势预警方法。
背景技术
风能是一种清洁可再生能源,取之不尽用之不竭,风力发电将迎来高速发展,而对风力发电机组的运行可靠性提出了更高的要求。
目前在风力发电机各设备温度监测预警方式上最普遍的是设置定值阈值报警,触发报警条件单一,当设备温度超过设定温度阈值时触发报警,这样的方式只考虑了设备最极限的温度值,无法根据不同工况条件进行动态设置阈值,无法反映设备缓慢劣化趋势。国内外的研究中使用较多的是ARIMA差分整合滑动平均自回归模型,但该方法在变量输入上只考虑自身历史数据,无法感知在风速、发电机转速、机舱温度等变化情况,造成预测准确度较低。
随着人工神经网络的逐渐成熟,将其与各领域传统经验、方法相融合,实现模拟人类认知,在很多领域已经取得了成功应用,近几年也有研究人员使用 BP神经网络用于风力发电机设备温度预测,其做法主要是确定输入因子及其对应的输出,建立BP神经网络包括输入层、隐藏层、输出层,通过正向传播计算损失函数,再反向传播,利用梯度下降和链式求导法则最小化损失函数,完成 BP神经网络训练。该方法可以有效的提取各输入因子的特征,建立输入与输出的深层关系,实现温度预测。但其不足之处在于没有考虑到温度数据的时间序列特性,对历史温度数据不存在记忆性,无法有效感知风力发电机组运行时风资源不稳定引起发电机转速、有功功率突变,造成预测精度较低。
发明内容
本发明的目的是提供一种基于LSTM的风机齿轮箱高速轴温度趋势预警方法,针对齿轮箱高速轴温度因风资源不可靠造成瞬变性、波动大和非线性的特点,设计基于LSTM神经网络的复合模型,解决现有齿轮箱高速轴温度监测手段无法实现高精度提前预测和温度趋势预警问题。
本发明提供了一种基于LSTM的风机齿轮箱高速轴温度趋势预警方法,包括如下步骤:
步骤1,获取风力发电机齿轮箱高速轴温度同一时间序列的其他相关参数历史数据,对数据进行预处理,通过插值法填补缺失值,处理异常值,归一化处理;
步骤2,进行皮尔逊相关性分析筛选出与齿轮箱高速轴温度相关性强的参数作为样本数据集;
步骤3,提取训练数据集与测试数据集,搭建LSTM神经网络;
步骤4,配置LSTM神经网络训练参数,对网络进行训练,使用LSTM神经网络高速轴温度趋势预测模型进行预测,通过均方根误差衡量模型的预测误差;
步骤5,计算高速轴预测值与实际值的残差,求得残差均值,设置预警逻辑得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。
进一步地,所述步骤1包括:
从风力发电机组本身配备的SCADA系统获取风力发电机实际运行的同一时间序列的风速、发电机有功功率、发电机转速、机舱温度、发电机绕组温度、偏航压力和齿轮箱高速轴温度7组历史数据;所采集的历史数据根据同一时间序列排序,时间间隔为60秒;
采用滑动窗口方法,设定窗口大小为n,对样本数据进行检验,判断连续n个时间步数据中出现缺失值的个数<m是否成立,若是,将使用拉格朗日插值法对缺失值进行插补,否则,直接删除缺失数据的时间步;
拉格朗日插值计算公式如下:
Figure BDA0003881727510000021
式中,坐标点(x,L(x))为窗口n内的缺失数据,x为缺失数据的时间坐标,L(x) 为经过拉格朗日插值法计算后的插入值,坐标点(xi,yi)为窗口n内的非缺失数据;
使用极差法对数据进行归一化处理,计算公式如下:
Figure BDA0003881727510000031
式中,x为数据归一化前的初始值,xmin为所述数据中的最小值,xmax为所述数据中的最大值。
进一步地,所述步骤2包括:
对步骤1中归一化处理后的风速、发电机有功功率、发电机转速、机舱温度、发电机绕组温度、偏航压力分别与齿轮箱高速轴温度进行皮尔逊相关性分析;
皮尔逊相关性系数计算如下:
Figure BDA0003881727510000032
Figure BDA0003881727510000033
Figure BDA0003881727510000034
Figure BDA0003881727510000035
Figure BDA0003881727510000036
Figure BDA0003881727510000037
式中,
Figure BDA0003881727510000038
分别表示风速集F和齿轮箱高速轴温度集A的均值;Cov(F,A)表示风速集F与齿轮箱高速轴温度集A的协方差;SF、SA分别表示风速集F和齿轮箱高速轴温度集A的标准差;rFA表示皮尔逊相关系数;使用相同的方法分别得到其余参数数据与齿轮箱高速轴温度的相关性系数值;
通过比较相关性系数绝对值的大小,筛选出风速、发电机有功功率、发电机转速和机舱温度和齿轮箱高速轴温度共5组历史数据作为样本数据,整合如下:
Figure BDA0003881727510000041
式中,IN为5组历史数据整合后的样本数据矩阵,F为风速集,P为发电机有功功率集,R为发电机转速集,E为机舱温度集,A为齿轮箱高速轴温度集,s 为数据集第s个时间步。
进一步地,所述步骤3包括:
将IN矩阵样本数据中前10%数据作为测试样本集O,后90%数据作为训练样本集k,在Q中从i=1个时间步开始单时间步滑动,每滑动一步提取长度为k的连续数据作为输入训练集;在齿轮箱高速轴温度集A中,提取第k个值作为输出标签集得到:
Figure BDA0003881727510000042
Ytrain=(Ak+u,Ak+u+1,…,As+u)T
式中,Xtrain作为LSTM神经网络的训练集输入;Ytrain作为LSTM神经网络的训练集标签;u∈[1,2,…,10]表示预测时间步长,即使用前k个时间步的训练集输入,预测第k+u个时间步的齿轮箱高速轴温度;使用相同的方法从 O中提取LSTM神经网络测试集输入Xtest其对应测试集标签Ytest
搭建LSTM神经网络,包括输入层、LSTM层、Dropout层、全连接层和回归输出层;其中输入层特征维度为5;LSTM层设为单层,设置LSTM层输入维度为128,输入样本个数为s-k+1,沿时间步扩展神经元个数为k,输出模式为序列最后一个;Dropout层,设置随机丢弃概率为0.2;全连接层输出层神经元个数为1;输出层为线性回归层得到预测标签;
所述LSTM神经元内部包含输入门it、遗忘门ft、输出门ot和内部记忆单元 ct,LSTM神经元输入含t时刻的序列输入xt,t-1时刻的记忆细胞ct-1、细胞状态ht-1;输出包含t时刻记忆细胞ct、细胞状态ht,其中记忆细胞ct代表长期记忆信息,细胞状态ht代表短期记忆信息,得到ct与ht具体计算公式如下:
it=σ(Wi·[ht-1,xt]+bi)
ft=σ(Wf·[ht-1,xt]+bf)
ot=σ(Wo·[ht-1,xt]+bo)
Figure BDA0003881727510000051
Figure BDA0003881727510000052
ht=ot⊙tanh(ct)
式中,σ、tanh分别表示sigmoid和双曲正切激活函数;it、ft、ot
Figure BDA0003881727510000053
ct、 ht分别表示输入门、遗忘门、输出门、记忆细胞输入状态的候选状态、记忆细胞和细胞状态;Wi、Wf、Wo、Wc分别表示对应的权重矩阵;bi、bf、bo、bc分别表示对应的偏置矩阵;⊙表示各元素按位相乘。
进一步地,所述步骤4包括:
配置LSTM神经网络训练参数;随初始化LSTM神经网络参数;选用adam 优化器,梯度阈值为1;选择CPU进行训练;设置最大训练轮数为100;设置最小批量为2000;设置最初学习率为0.005,学习率计划为分段学习,学习率下降周期为125,学习率下降系数为0.2;
将训练集输入Xtrain和训练集标签Ytrain放入在网络进行训练,过程中采用时间反向传播算法,根据输入的时间序列不断计算每个时刻输出值与真实标签值的误差值,再逐级将误差项传向上一层,重新更新相关权重矩阵和偏置矩阵的值,训练完成后得到高速轴温度趋势预测模型,将测试集输入Xtest输入到训练好的模型中,得到归一化的齿轮箱高速轴温度预测值Ypre,再对其进行反归一化得到最终齿轮箱高速轴温度预测值;
选用RMSE均方根误差作为误差的衡量标准,公式如下:
Figure BDA0003881727510000054
式中,yi为实际高速轴温度;
Figure BDA0003881727510000061
为预测的高速轴温度。
进一步地,所述步骤5包括:
利用步骤4得到的高速轴温度趋势预测模型进行预测,计算齿轮箱高速轴温度预测值与实际值之间的残差,求出残差均值,设置预警逻辑得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。
借由上述方案,通过基于LSTM的风机齿轮箱高速轴温度趋势预警方法,针对齿轮箱高速轴温度因风资源不可靠造成瞬变性、波动大和非线性的特点,通过基于LSTM神经网络的复合模型,解决了现有齿轮箱高速轴温度监测手段无法实现高精度提前预测和温度趋势预警问题。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附图详细说明如后。
附图说明
图1是本发明的实现流程图;
图2是本发明搭建的LSTM神经网络模型结构;
图3是本发明LSTM神经网络模型预测效果;
图4是本发明LSTM神经网络模型预测效果(局部);
图5是本发明预测值与实际值的残差及残差阈值图;
图6是本发明设置的高速轴温度预警逻辑。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
术语解释:
LSTM(Long Short-Term Memory)长短期记忆神经网络是一种特殊RNN循环神经网络,LSTM的神经元通过增加输入门it、遗忘门ft、输出门ot和内部记忆单元ct,具有选择性记忆功能,可以选择记忆重要信息,过滤掉噪声信息,有效解决了梯度消失和梯度爆炸问题。
参图1所示,本实施例提供了一种基于LSTM的风机齿轮箱高速轴温度趋势预警方法,包括如下步骤:
S01、获取风力发电机齿轮箱高速轴温度同一时间序列的其他相关参数历史数据,对数据预处理,通过插值法填补缺失值,处理异常值,归一化处理;
S02、进行皮尔逊相关性分析筛选出与齿轮箱高速轴温度相关性强的参数作为样本数据集;
S03、提取训练数据集与测试数据集,搭建LSTM神经网络;
S04、配置LSTM神经网络训练参数,对网络进行训练,使用LSTM神经网络高速轴温度趋势预测模型进行预测,通过均方根误差衡量模型的预测误差;
S05、计算高速轴预测值与实际值的残差,求得残差均值,设置预警逻辑得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。
对于步骤S01:获取风力发电机齿轮箱高速轴温度同一时间序列的其他相关参数历史数据,对数据预处理,通过插值法填补缺失值,处理异常值,归一化处理。包括以下步骤:
具体的,从重庆市丰都县某风场11号风力发电机组本身配备的SCADA系统(数据采集与监视控制系统)获取2021年8月24日至2021年11月30日机组正常运行期间共129305条数据,采样间隔为60S。数据源包括同一时间序列的风速、发电机有功功率、发电机转速、机舱温度、发电机绕组温度、偏航压力和齿轮箱高速轴温度7组参数。
具体的,采用滑动窗口方法,设定窗口大小为n,对样本数据进行检验,判断连续n个时间步数据中出现缺失值的个数<m是否成立,若是,将使用拉格朗日插值法对缺失值进行插补,否则,直接删除缺失数据的时间步。
拉格朗日插值计算公式如下:
Figure BDA0003881727510000071
式中,坐标点(x,L(x))为窗口n内的缺失数据,x为缺失数据的时间坐标,L(x) 为经过拉格朗日插值法计算后的插入值,坐标点(xi,yi)为窗口n内的非缺失数据。本实例中设置n=10,m=2。
具体的,使用极差法对数据进行归一化处理,计算公式如下:
Figure BDA0003881727510000081
式中,x为数据归一化前的初始值,xmin为所述数据中的最小值,xmax为所述数据中的最大值。
在本实施例中,步骤S02包括以下步骤:
对步骤01中归一化后的风速、发电机有功功率、发电机转速、机舱温度、发电机绕组温度、偏航压力分别与齿轮箱高速轴温度进行皮尔逊相关性分析。
皮尔逊相关性系数计算如下:
Figure BDA0003881727510000082
Figure BDA0003881727510000083
Figure BDA0003881727510000084
Figure BDA0003881727510000085
Figure BDA0003881727510000086
Figure BDA0003881727510000087
式中,
Figure BDA0003881727510000088
分别表示风速集F和齿轮箱高速轴温度集A的均值;Cov(F,A)表示风速集F与齿轮箱高速轴温度集A的协方差;SF、SA分别表示风速集F和齿轮箱高速轴温度集A的标准差;rFA表示皮尔逊相关系数。使用相同的方法可分别得到其余参数数据与齿轮箱高速轴温度的相关性系数值。
通过比较相关性系数绝对值的大小,筛选出风速、发电机有功功率、发电机转速和机舱温度和齿轮箱高速轴温度共5组历史数据作为样本数据,整合如下:
Figure BDA0003881727510000091
式中,IN为5组历史数据整合后的样本数据矩阵,F为风速集,P为发电机有功功率集,R为发电机转速集,E为机舱温度集,A为齿轮箱高速轴温度集,s 为数据集第s个时间步。
在本实施例中,步骤S03包括以下步骤:
提取训练数据集与测试数据集,搭建LSTM神经网络。
具体的,从预处理好的IN矩阵样本数据中将前10%共11520条序列数据作为测试样本集O,后90%共117682条序列数据作为训练样本集Q,在Q中从i=1 个时间步开始单时间步滑动,每滑动一步提取长度为k的连续数据作为输入训练集。在齿轮箱高速轴温度集A中,提取第k个值作为输出标签集得到:
Figure BDA0003881727510000092
Ytrain=(Ak+u,Ak+u+1,…,As+u)T
式中,Xtrain作为LSTM神经网络的训练集输入;Ytrain作为LSTM神经网络的训练集标签;u∈[1,2,…,10]表示预测时间步长,即使用前k个时间步的输入,预测第k+u个时间步的齿轮箱高速轴温度。使用相同的方法从O中提取LSTM神经网络测试集输入Xtest其对应测试集标签Ytest。本实例中令k=30、 u=1,即使用前30个时间步的输入,预测第31个时间步的齿轮箱高速轴温度。
具体的,如图2所示,搭建LSTM神经网络,包括输入层、LSTM层、Dropout 层、全连接层和回归输出层。其中输入层特征维度为5;LSTM层设为单层,设置LSTM层输入维度为128,输入样本个数为s-k+1,沿时间步扩展神经元个数为k,输出模式为序列最后一个;Dropout层,设置随机丢弃概率为0.2;全连接层输出层神经元个数为1;输出层为线性回归层得到预测标签。
所述LSTM神经元内部包含输入门it、遗忘门ft、输出门ot和内部记忆单元ct,LSTM神经元输入含t时刻的序列输入xi,t-1时刻的记忆细胞ct-1、细胞状态ht-1;输出包含t时刻记忆细胞ct、细胞状态ht,其中记忆细胞ct代表长期记忆信息,细胞状态ht代表短期记忆信息,得到ct与ht具体计算公式如下:
it=σ(Wi·[ht-1,xt]+bi)
ft=σ(Wf·[ht-1,xt]+bf)
ot=σ(Wo·[ht-1,xt]+bo)
Figure BDA0003881727510000101
Figure BDA0003881727510000102
ht=ot⊙tanh(ct)
式中,σ、tanh分别表示sigmoid和双曲正切激活函数;it、ft、ot
Figure BDA0003881727510000103
ct、 ht分别表示输入门、遗忘门、输出门、记忆细胞输入状态的候选状态、记忆细胞和细胞状态;Wi、Wf、Wo、Wc分别表示对应的权重矩阵;bi、bf、bo、bc分别表示对应的偏置矩阵;⊙表示各元素按位相乘。
在本实施例中,步骤S04包括以下步骤:
配置LSTM神经网络训练参数,对网络进行训练,使用LSTM神经网络高速轴温度趋势预测模型进行预测,通过均方根衡量模型的预测误差。
具体的,随初始化LSTM神经网络参数;选用adam优化器,梯度阈值为1;选择CPU进行训练;设置最大训练轮数为100;设置最小批量为2000;设置最初学习率为0.005,学习率计划为分段学习,学习率下降周期为125,学习率下降系数为0.2。
具体的,对网络进行训练。将训练集输入Xtrain和训练集标签Ytrain放入在网络进行训练,过程中采用时间反向传播算法,根据输入的时间序列不断计算每个时刻输出值与真实标签值的误差值,再逐级将误差项传向上一层,重新更新相关权重矩阵和偏置矩阵的值,训练完成后得到高速轴温度趋势预测模型。
具体的,如图3,图4,带入测试集数据,通过均方根衡量模型的预测误差。将测试集输入Xtest输入到训练好的模型中,得到归一化的齿轮箱高速轴温度预测值Ypre,再对其进行反归一化得到最终齿轮箱高速轴温度预测值。
选用RMSE均方根误差作为误差的衡量标准,公式如下:
Figure BDA0003881727510000111
式中,yi为实际高速轴温度;
Figure BDA0003881727510000112
为预测的高速轴温度。本实例中RMSE= 0.541。
在本实施例中,步骤S05包括以下步骤:
计算高速轴预测值与实际值的残差,求得残差均值,设置预警逻辑得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。
具体的,如图5所示,利用步骤S04得到的高速轴温度趋势预测模型进行预测,计算齿轮箱高速轴温度预测值与实际值之间的残差,求出残差均值。设置预警逻辑,如图6,设置一个滑动窗口长度为g=10,沿预测时间步逐步滑动,如果残差与残差均值相减绝对值大于4倍均方根误差时,滑动窗口内的当前时间步数值计“1”,为超过则计“0”。当窗口内出现计“1”的位数连续且大于等于3个或者大于4个时发出红色预警,否则判断窗口内出现计“1”的位数是否满足大于等于2个且小于等于4个,满足时发出黄色预警,否则不发出报警,得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。发出红色预警、发出黄色预警或不发出报警后重新开始,进入计算齿轮箱高速轴温度预测值与实际值之间的残差,求出残差均值步骤。
该基于LSTM的风机齿轮箱高速轴温度趋势预警方法,针对齿轮箱高速轴温度因风资源不可靠造成瞬变性、波动大和非线性的特点,通过基于LSTM神经网络的复合模型,解决了现有齿轮箱高速轴温度监测手段无法实现高精度提前预测和温度趋势预警问题。
以上所述仅是本发明的优选实施方式,并不用于限制本发明,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。

Claims (6)

1.一种基于LSTM的风机齿轮箱高速轴温度趋势预警方法,其特征在于,包括如下步骤:
步骤1,获取风力发电机齿轮箱高速轴温度同一时间序列的其他相关参数历史数据,对数据进行预处理,通过插值法填补缺失值,处理异常值,归一化处理;
步骤2,进行皮尔逊相关性分析筛选出与齿轮箱高速轴温度相关性强的参数作为样本数据集;
步骤3,提取训练数据集与测试数据集,搭建LSTM神经网络;
步骤4,配置LSTM神经网络训练参数,对网络进行训练,使用LSTM神经网络高速轴温度趋势预测模型进行预测,通过均方根误差衡量模型的预测误差;
步骤5,计算高速轴预测值与实际值的残差,求得残差均值,设置预警逻辑得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。
2.根据权利要求1所述的基于LSTM的风机齿轮箱高速轴温度趋势预警方法,其特征在于,所述步骤1包括:
从风力发电机组本身配备的SCADA系统获取风力发电机实际运行的同一时间序列的风速、发电机有功功率、发电机转速、机舱温度、发电机绕组温度、偏航压力和齿轮箱高速轴温度7组历史数据;所采集的历史数据根据同一时间序列排序,时间间隔为60秒;
采用滑动窗口方法,设定窗口大小为n,对样本数据进行检验,判断连续n个时间步数据中出现缺失值的个数<m是否成立,若是,将使用拉格朗日插值法对缺失值进行插补,否则,直接删除缺失数据的时间步;
拉格朗日插值计算公式如下:
Figure FDA0003881727500000011
式中,坐标点(x,L(x))为窗口n内的缺失数据,x为缺失数据的时间坐标,L(x)为经过拉格朗日插值法计算后的插入值,坐标点(xi,yi)为窗口n内的非缺失数据;
使用极差法对数据进行归一化处理,计算公式如下:
Figure FDA0003881727500000021
式中,x为数据归一化前的初始值,xmin为所述数据中的最小值,xmax为所述数据中的最大值。
3.根据权利要求2所述的基于LSTM的风机齿轮箱高速轴温度趋势预警方法,其特征在于,所述步骤2包括:
对步骤1中归一化处理后的风速、发电机有功功率、发电机转速、机舱温度、发电机绕组温度、偏航压力分别与齿轮箱高速轴温度进行皮尔逊相关性分析;
皮尔逊相关性系数计算如下:
Figure FDA0003881727500000022
Figure FDA0003881727500000023
Figure FDA0003881727500000024
Figure FDA0003881727500000025
Figure FDA0003881727500000026
Figure FDA0003881727500000027
式中,
Figure FDA0003881727500000028
分别表示风速集F和齿轮箱高速轴温度集A的均值;Cov(F,A)表示风速集F与齿轮箱高速轴温度集A的协方差;SF、SA分别表示风速集F和齿轮箱高速轴温度集A的标准差;rFA表示皮尔逊相关系数;使用相同的方法分别得到其余参数数据与齿轮箱高速轴温度的相关性系数值;
通过比较相关性系数绝对值的大小,筛选出风速、发电机有功功率、发电机转速和机舱温度和齿轮箱高速轴温度共5组历史数据作为样本数据,整合如下:
Figure FDA0003881727500000031
式中,IN为5组历史数据整合后的样本数据矩阵,F为风速集,P为发电机有功功率集,R为发电机转速集,E为机舱温度集,A为齿轮箱高速轴温度集,s为数据集第s个时间步。
4.根据权利要求3所述的基于LSTM的风机齿轮箱高速轴温度趋势预警方法,其特征在于,所述步骤3包括:
将IN矩阵样本数据中前10%数据作为测试样本集O,后90%数据作为训练样本集Q,在Q中从i=1个时间步开始单时间步滑动,每滑动一步提取长度为k的连续数据作为输入训练集;在齿轮箱高速轴温度集A中,提取第k个值作为输出标签集得到:
Figure FDA0003881727500000032
Ytrain=(Ak+u,Ak+u+1,…,As+u)T
式中,Xtrain作为LSTM神经网络的训练集输入;Ytrain作为LSTM神经网络的训练集标签;u∈[1,2,…,10]表示预测时间步长,即使用前k个时间步的训练集输入,预测第k+u个时间步的齿轮箱高速轴温度;使用相同的方法从O中提取LSTM神经网络测试集输入Xtest其对应测试集标签Ytest
搭建LSTM神经网络,包括输入层、LSTM层、Dropout层、全连接层和回归输出层;其中输入层特征维度为5;LSTM层设为单层,设置LSTM层输入维度为128,输入样本个数为s-k+1,沿时间步扩展神经元个数为k,输出模式为序列最后一个;Dropout层,设置随机丢弃概率为0.2;全连接层输出层神经元个数为1;输出层为线性回归层得到预测标签;
所述LSTM神经元内部包含输入门it、遗忘门ft、输出门ot和内部记忆单元ct,LSTM神经元输入含t时刻的序列输入xt,t-1时刻的记忆细胞ct-1、细胞状态ht-1;输出包含t时刻记忆细胞ct、细胞状态ht,其中记忆细胞ct代表长期记忆信息,细胞状态ht代表短期记忆信息,得到ct与ht具体计算公式如下:
it=σ(Wi·[ht-1,xt]+bi)
ft=σ(Wf·[ht-1,xt]+bf)
ot=σ(Wo·[ht-1,xt]+bo)
Figure FDA0003881727500000041
Figure FDA0003881727500000042
ht=ot⊙tanh(ct)
式中,σ、tanh分别表示sigmoid和双曲正切激活函数;it、ft、ot
Figure FDA0003881727500000043
ct、ht分别表示输入门、遗忘门、输出门、记忆细胞输入状态的候选状态、记忆细胞和细胞状态;Wi、Wf、Wo、Wc分别表示对应的权重矩阵;bi、bf、bo、bc分别表示对应的偏置矩阵;⊙表示各元素按位相乘。
5.根据权利要求4所述的基于LSTM的风机齿轮箱高速轴温度趋势预警方法,其特征在于,所述步骤4包括:
配置LSTM神经网络训练参数;随初始化LSTM神经网络参数;选用adam优化器,梯度阈值为1;选择CPU进行训练;设置最大训练轮数为100;设置最小批量为2000;设置最初学习率为0.005,学习率计划为分段学习,学习率下降周期为125,学习率下降系数为0.2;
将训练集输入Xtrain和训练集标签Ytrain放入在网络进行训练,过程中采用时间反向传播算法,根据输入的时间序列不断计算每个时刻输出值与真实标签值的误差值,再逐级将误差项传向上一层,重新更新相关权重矩阵和偏置矩阵的值,训练完成后得到高速轴温度趋势预测模型,将测试集输入Xtest输入到训练好的模型中,得到归一化的齿轮箱高速轴温度预测值Ypre,再对其进行反归一化得到最终齿轮箱高速轴温度预测值;
选用RMSE均方根误差作为误差的衡量标准,公式如下:
Figure FDA0003881727500000051
式中,yi为实际高速轴温度;
Figure FDA0003881727500000052
为预测的高速轴温度。
6.根据权利要求5所述的基于LSTM的风机齿轮箱高速轴温度趋势预警方法,其特征在于,所述步骤5包括:
利用步骤4得到的高速轴温度趋势预测模型进行预测,计算齿轮箱高速轴温度预测值与实际值之间的残差,求出残差均值,设置预警逻辑得到最终的风力发电机齿轮箱高速轴温度趋势预警复合模型。
CN202211237677.8A 2022-10-10 2022-10-10 一种基于lstm的风机齿轮箱高速轴温度趋势预警方法 Pending CN115758290A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211237677.8A CN115758290A (zh) 2022-10-10 2022-10-10 一种基于lstm的风机齿轮箱高速轴温度趋势预警方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211237677.8A CN115758290A (zh) 2022-10-10 2022-10-10 一种基于lstm的风机齿轮箱高速轴温度趋势预警方法

Publications (1)

Publication Number Publication Date
CN115758290A true CN115758290A (zh) 2023-03-07

Family

ID=85351175

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211237677.8A Pending CN115758290A (zh) 2022-10-10 2022-10-10 一种基于lstm的风机齿轮箱高速轴温度趋势预警方法

Country Status (1)

Country Link
CN (1) CN115758290A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116205375A (zh) * 2023-04-18 2023-06-02 北京市南水北调团城湖管理处 一种泵站机组运行状态预测方法和系统
CN116432542A (zh) * 2023-06-12 2023-07-14 国网江西省电力有限公司电力科学研究院 一种基于误差序列修正的开关柜母排温升预警方法及系统
CN116452067A (zh) * 2023-06-12 2023-07-18 昆山工业大数据创新中心有限公司 一种工业设备参数的劣化度评估方法、装置及设备
CN116995276A (zh) * 2023-09-27 2023-11-03 爱德曼氢能源装备有限公司 燃料电池发电系统冷却方法及系统
CN117272245A (zh) * 2023-11-21 2023-12-22 陕西金元新能源有限公司 一种风机齿轮箱温度预测方法、装置、设备和介质
CN117542169A (zh) * 2023-11-07 2024-02-09 国网江苏省电力有限公司镇江供电分公司 基于大数据分析的自动化设备温度异常预警方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116205375A (zh) * 2023-04-18 2023-06-02 北京市南水北调团城湖管理处 一种泵站机组运行状态预测方法和系统
CN116205375B (zh) * 2023-04-18 2023-09-01 北京市南水北调团城湖管理处 一种泵站机组运行状态预测方法和系统
CN116432542A (zh) * 2023-06-12 2023-07-14 国网江西省电力有限公司电力科学研究院 一种基于误差序列修正的开关柜母排温升预警方法及系统
CN116452067A (zh) * 2023-06-12 2023-07-18 昆山工业大数据创新中心有限公司 一种工业设备参数的劣化度评估方法、装置及设备
CN116452067B (zh) * 2023-06-12 2023-09-05 昆山工业大数据创新中心有限公司 一种工业设备参数的劣化度评估方法、装置及设备
CN116432542B (zh) * 2023-06-12 2023-10-20 国网江西省电力有限公司电力科学研究院 一种基于误差序列修正的开关柜母排温升预警方法及系统
CN116995276A (zh) * 2023-09-27 2023-11-03 爱德曼氢能源装备有限公司 燃料电池发电系统冷却方法及系统
CN116995276B (zh) * 2023-09-27 2023-12-29 爱德曼氢能源装备有限公司 燃料电池发电系统冷却方法及系统
CN117542169A (zh) * 2023-11-07 2024-02-09 国网江苏省电力有限公司镇江供电分公司 基于大数据分析的自动化设备温度异常预警方法
CN117272245A (zh) * 2023-11-21 2023-12-22 陕西金元新能源有限公司 一种风机齿轮箱温度预测方法、装置、设备和介质
CN117272245B (zh) * 2023-11-21 2024-03-12 陕西金元新能源有限公司 一种风机齿轮箱温度预测方法、装置、设备和介质

Similar Documents

Publication Publication Date Title
CN115758290A (zh) 一种基于lstm的风机齿轮箱高速轴温度趋势预警方法
CN112949945B (zh) 一种改进双向长短期记忆网络的风电功率超短期预测方法
CN109766583B (zh) 基于无标签、不均衡、初值不确定数据的航空发动机寿命预测方法
CN110738360B (zh) 一种设备剩余寿命预测方法及系统
CN111340292B (zh) 一种基于聚类的集成神经网络pm2.5预测方法
CN108596327B (zh) 一种基于深度学习的地震速度谱人工智能拾取方法
CN111814956B (zh) 一种基于多维度二次特征提取的多任务学习的空气质量预测方法
CN109583565A (zh) 基于注意力模型长短时记忆网络的洪水预测方法
CN114218872B (zh) 基于dbn-lstm半监督联合模型的剩余使用寿命预测方法
CN114282443B (zh) 基于mlp-lstm有监督联合模型的剩余使用寿命预测方法
CN111352977A (zh) 基于自注意力双向长短期记忆网络的时序数据监测方法
CN114781744A (zh) 基于编码解码器的深度学习多步长辐照度预测方法
CN115809405A (zh) 基于多特征融合的风机主轴齿轮箱温度异常检测方法
CN115293197A (zh) 一种基于长短期记忆网络的钻孔应变数据异常检测方法
CN114819340A (zh) 一种用于天然气日负荷的时序预测方法
CN115877483A (zh) 一种基于随机森林和gru的台风路径预报方法
Zhang et al. Remaining useful life prediction for rolling bearings with a novel entropy-based health indicator and improved particle filter algorithm
CN116793666A (zh) 一种基于lstm-mlp-lsgan模型的风电机组齿轮箱故障诊断方法
CN115330082A (zh) 一种基于注意力机制的lstm网络的pm2.5浓度预测方法
CN115577748A (zh) 一种集成挤压激励注意力机制的双通道风电预测方法
CN112598186B (zh) 一种基于改进的lstm-mlp的小型发电机故障预测方法
CN114841063A (zh) 一种基于深度学习的航空发动机剩余寿命预测方法
CN115809947A (zh) 一种水污染时序数据连续缺失值填补方法
CN114386672A (zh) 环境大数据物联网智能检测系统
CN112560252A (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