CN102768115B - 一种风电机组齿轮箱健康状态实时动态监控方法 - Google Patents
一种风电机组齿轮箱健康状态实时动态监控方法 Download PDFInfo
- Publication number
- CN102768115B CN102768115B CN201210215628.4A CN201210215628A CN102768115B CN 102768115 B CN102768115 B CN 102768115B CN 201210215628 A CN201210215628 A CN 201210215628A CN 102768115 B CN102768115 B CN 102768115B
- Authority
- CN
- China
- Prior art keywords
- gear case
- health status
- health
- operating condition
- evaluation model
- 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.)
- Expired - Fee Related
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了旋转机械状态监测与故障诊断技术领域中的一种风电机组齿轮箱健康状态实时动态监控方法。本发明采集了齿轮箱振动信号和齿轮箱运行与监测数据;分析得到齿轮箱振动信号的时域特征、频域特征和包络信号的包络谱特征;建立了齿轮箱运行工况空间,并根据齿轮箱运行数据将齿轮箱运行工况空间划分为设定个齿轮箱运行工况子空间;得到运行监控历史数据子集,进而生成齿轮箱各部件的健康状态特征库;建立健康状态评价模型,得到各工况的健康状态实时评价模型和各部件的健康指数;最后求得齿轮箱的健康指数。本发明建立了健康状态评价模型,避开了复杂动态运行工况状态监测阈值确定的复杂性;避免了多特征参数简单递加造成线性失真的缺点。
Description
技术领域
本发明属于旋转机械状态监测与故障诊断技术领域,尤其涉及一种风电机组齿轮箱健康状态实时动态监控方法。
背景技术
齿轮箱是风电机组的关键部件,担负功率传递和变速功能。由于经常处于复杂恶劣的运行环境(交变载荷、潮湿、风沙、震动、高温、腐蚀等),风电机组齿轮箱故障率极高。据统计,运行5年后风电场齿轮箱损坏率高达40%~50%,极个别齿轮箱更换率几乎接近100%。近年来,齿轮箱的故障率逐年增高,已成为风电场机组可用率低、运行维护费用高的最主要原因。
目前,可以通过两方面的工作来降低风机齿轮箱的故障率,一方面是通过优化设计和提高制造、安装质量来提高齿轮箱可靠性;另一方面是利用状态监测技术,及时评估齿轮箱健康状态,并根据其健康状态和衰退趋势调整运行、合理安排维修。
当前主要通过在风机传动系统布置振动传感器和温度传感器来监测齿轮箱的运行状态,当测量值超过预先设定的阈值时报警或停机。然而,由于缺乏有效的特征提取和多元状态信息融合,很难实现齿轮箱整体健康状态的准确评估;另外,振动和温度状态阈值采用固定值,没有考虑运行工况(风速、转速、功率、扭矩、空气密度等)变化对阈值的影响,也影响了齿轮箱健康状态评估的准确性。
综上所述,适当配置齿轮箱健康状态监测测点,开展工况辨识基础上的基于流形学习和多元信息融合的复杂工况齿轮箱健康状态监控研究,实现齿轮箱健康状态的实时动态评估和故障预示,对提高风电场机组的可用率,降低运行与维护费用具有重要的现实意义。
发明内容
针对上述背景技术中提到当前风电机组齿轮箱健康状态评价在信息融合和阈值确定方面存在的不足,本发明提出了一种风电机组齿轮箱健康状态实时动态监控方法。
本发明的技术方案是,一种风电机组齿轮箱健康状态实时动态监控方法,其特征是该方法包括以下步骤:
步骤1:采集齿轮箱振动信号和齿轮箱运行与监测数据;
步骤2:对齿轮箱振动信号进行分析,得到齿轮箱振动信号的时域特征、频域特征和包络信号的包络谱特征;
步骤3:建立齿轮箱运行工况空间,根据齿轮箱运行与监测数据将齿轮箱运行工况空间划分为设定个齿轮箱运行工况子空间;
步骤4:在步骤2和步骤3的基础上,将齿轮箱时域特征、频域特征、运行与监测数据划分为运行监控历史数据子集;
步骤5:对运行监控历史数据子集进行特征提取,生成齿轮箱各部件的健康状态特征库;
步骤6:根据齿轮箱各部件的健康状态特征库建立健康状态评价模型,通过运行监控历史数据子集对健康状态评价模型进行训练,得到各工况的健康状态实时评价模型,通过健康状态实时评价模型对实时的齿轮箱运行与监测数据进行评价,得到各部件的健康指数;
步骤7:将各部件的健康指数带入变权综合公式,得到齿轮箱的健康指数。
所述齿轮振动信号时域特征为样本时间序列平均值、样本时间序列有效值、峭度指标、峰值指标、歪度指标和熵。
所述齿轮振动信号频域特征为滚动轴承故障特征频率、齿轮啮合频率处的能量和包络谱各故障特征频率对应的能量。
所述齿轮箱运行工况空间划分的方法是聚类分析方法。
所述对运行监控历史数据子集进行特征提取的方法是流形学习法。
所述健康状态评价模型为Logistics回归模型。
所述部件的健康指数的计算公式为:
其中:
为ti时刻对应第p个运行工况子空间的齿轮箱部件的健康指数;
为第p个工况子空间对应的健康状态评估函数;
yi为时刻ti齿轮箱相应部件特征提取后状态特征向量;
Θ(p)为第p个工况子空间对应健康状态评价模型的L+1个模型参数,
为健康状态评价模型的第一个模型参数;
为健康状态评价模型第L+1个模型参数;
yil为ti时刻第l个健康状态特征参数值;
ε为残差。
所述齿轮箱的健康指数为:
其中:
HI为ti时刻齿轮箱的健康指数;
为第j个部件健康状态评价时占的常权权重,j=1,2,…m;
为ti时刻第j个部件健康指数的γ次幂;
为第k个部件健康状态评价时占的常权权重,k=1,2,…m;
为ti时刻第k个部件健康指数的γ-1次幂;
γ为均衡系数。
本发明具有如下显著优点:本发明方法将风电机组复杂运行工况进行辨识基础上的分割,在各运行工况子空间建立健康状态评价模型,避开了复杂动态运行工况状态阈值确定的复杂性;采用基于统计模式识别的状态评价避免了现有状态评价算法中采用多特征参数线性叠加拟合可能造成的线性失真的缺点;另外通过模型的定期修正实现健康状态的动态监控,提高状态评价精度。本发明使得风电机组齿轮箱健康状态评价更加科学和全面,该方法设计合理,实用性强,为风电机组齿轮箱健康状态监控提供有效的技术支持。
附图说明
图1为风电机组齿轮箱健康状态实时动态监控总体流程图;
图2为风电机组齿轮箱结构及测点布置;
图3为风电机组齿轮箱信号处理、特征提取及运行工况辨识;
图4为风电机组齿轮箱健康状态实时动态评价流程;
图5为风电机组齿轮箱健康状态评估层次结构。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
本发明的步骤如下:
步骤1:采集齿轮箱振动信号和齿轮箱运行与监测数据;
步骤2:对齿轮箱振动信号进行分析,得到齿轮箱振动信号的时域特征、频域特征和包络信号的包络谱特征;
步骤3:建立齿轮箱运行工况空间,根据齿轮箱运行与监测数据将齿轮箱运行工况空间划分为设定个齿轮箱运行工况子空间;
步骤4:在步骤2和步骤3的基础上,将齿轮箱时域特征、频域特征、运行与监测数据划分为运行监控历史数据子集;
步骤5:对运行监控历史数据子集进行特征提取,生成齿轮箱各部件的健康状态特征库;
步骤6:根据齿轮箱各部件的健康状态特征库建立健康状态评价模型,通过运行监控历史数据子集对健康状态评价模型进行训练,得到各工况的健康状态实时评价模型,通过健康状态实时评价模型对实时的齿轮箱运行与监测数据进行评价,得到各部件的健康指数;
步骤7:将各部件的健康指数带入变权综合公式,得到齿轮箱的健康指数。
本发明通过合理配置齿轮箱振动传感器,采集齿轮箱振动信号,并进行信号处理,获得原始特征;结合运行工况历史数据,进行工况辨识,划分为若干工况子空间;开展基于流形学习的特征提取,并建立特征库;利用振动监测和运行历史数据训练齿轮箱多工况健康状态评估模型,并建立模型库;对运行的风机齿轮箱实时计算状态特征值,并利用相应工况的健康状态评价模型实时评价健康状态;最后,根据健康状态和衰退趋势调整运行并合理安排维修,最终实现风电机组齿轮箱健康状态的实时动态监控。
分别在风机齿轮箱的环形齿轮、行星齿轮支架、低速轴电机侧、中速轴电机侧、高速轴电机侧和高速轴轴承(迎风侧和背风侧)配置振动传感器。
本发明采用的技术方案如下:
1、数据采集及存储
根据齿轮箱的形式确定振动监测测点布置方案,根据各测点测量要求合理选用振动传感器。根据测量数据存储和分析显示要求设计并研制振动数据采集装置。通过连接机组数据采集与监控系统(SCADA)获取风电机组运行与监测数据,存储采集到的齿轮箱振动信号和齿轮箱运行与监测数据。
2、信号处理及分析
对齿轮箱各测点(如图2所示)振动历史数据进行时域分析并计算时域原始特征。计算出的时域原始特征包括样本时间序列平均值、样本时间序列有效值、峭度指标、峰值指标、歪度指标和熵,相应的计算公式如下:
①样本时间序列平均值
②样本时间序列有效值
③峭度指标
④峰值指标
⑤歪度指标
⑥熵
式中,N为时间序列样本个数,vi为样本时间序列。
快速傅里叶变换(FFT)可以将某一波形分解成一系列不同频率正弦波的和。有限区间快速离散傅里叶变换公式如下:
利用快速傅里叶变换(FFT)算法,将齿轮箱振动信号从时域转化为频率域。整个频谱可以分为一些子段,每一段由一些相近频率组成。这些子段的中心频率通常事先定义,比如机械系统滚动轴承故障特征频率(BPFI、BPFO、BSF)、齿轮啮合频率等。计算这些子段中心频率处的能量,作为齿轮、轴承健康状态评价的频域原始特征。
为进一步分析某一特征频率,利用希尔伯特变换(Hilbert)获得包络信号,对包络信号作谱分析,得到包络信号的包络谱,计算包络谱轴承故障特征频率(BPFI、BPFO、BSF)和齿轮啮合频率处能量,也将其作为健康状态评价的频域原始特征。Hilbert变换定义如下:
3、工况辨识、特征提取及健康评价模型训练
建立风电机组齿轮箱运行工况特征集u(包括风速、风向、空气密度、扭矩、转速、输出功率等)和齿轮箱运行工况空间Ω。假定齿轮箱运行工况空间Ω可利用某一聚类或空间分割算法fc(u)聚类成P个齿轮箱运行工况子空间,即Ω=(O1,O2,…,OP),则某一运行工况对各工况子空间的隶属度集C可表示为:
C=fc(u)=(C1,C2,…,CP)(9)
其中,Cp为某运行工况隶属于第p个齿轮箱运行工况子空间的隶属度,Cp=Member(u∈Op)。
上式中fc(u)可以是k-means或Fuzzyc-means等聚类算法。
4、运行监控历史数据划分
通过信号分析与处理,生成了包含齿轮箱运行与监测数据时间序列、振动信号时域、频域特征值时间序列的原始运行监控历史数据集X。
X=(x1,x2,…,xN)(10)
其中,x为M维原始状态特征向量;N为时间序列样本个数。
利用工况辨识确定各时刻样本所处的运行工况子空间,进而按照所处子空间将N个样本划分到P个齿轮箱运行工况子空间,得到P个运行监控历史数据子集:
X={X(1),X(2),…,X(P)}(11)
各运行监控历史数据子集用于相应运行工况子空间健康评价模型的训练。
5、特征提取
对齿轮箱各部件运行监控历史数据子集X(1),X(2),…,X(P)分别进行基于流形学习的特征提取,生成齿轮箱各部件对应运行工况子空间的特征提取变换矩阵A(具体过程见15页)。状态特征提取的变换函数如下:
xi→yi=ATxi,A=(a0,a1,…,aL-1)(12)
其中:
xi为时刻ti对应的M维原始状态特征值向量,即运行监控历史数据子集中的元素;
yi为时刻ti齿轮箱相应部件特征提取后状态特征向量;
A为变换矩阵;
a为变换矩阵的特征向量。
6、健康状态评价模型的建立
数据驱动的设备健康状态通用评价模型可表示如下:
(ti,ui,yi)→(ti,zi)(13)
其中:
ti为某一运行时刻;
ui为时刻ti运行工况特征向量;
zi为ti时刻设备的健康指数。
针对任一运行工况子空间Op,利用相应运行与监控历史数据子集X(p)进行健康状态评价模型的训练,得到对应工况的健康状态实时评价模型(模型参数为Θ(p)),则齿轮箱部件的健康指数可按下式计算:
其中:
为ti时刻对应第p个运行工况子空间的齿轮箱部件的健康指数;
为第p个工况子空间对应的健康状态评估函数;
Θ(p)为第p个工况子空间对应健康状态评价模型的L+1个模型参数,
为健康状态评价模型的第一个模型参数;
为健康状态评价模型第L+1个模型参数;
yil为ti时刻第l个健康状态特征参数值;
ε为残差。
整个齿轮箱的健康状态由其相应的主要组成部件的健康状态决定,而各组成部件的健康状态可以通过运行工况辨识、振动信号分析、特征提取及健康状态评估确定。建立相应层次结构后可以实现从底层到顶层的健康状态逐级评价。
7、针对齿轮箱某一运行时刻ti,通过健康状态实时评价模型对其运行工况数据、运行实时数据和状态监测数据进行信号处理和特征计算,得到该时刻运行工况向量和第j个部件状态特征向量(ui,yij),yij为时刻ti对应的第j个部件变换后状态特征向量。对其运行工况参数进行工况辨识,如果识别为某一已知工况,则按下式计算第j个部件的健康指数zij:
Cp=Member(ui∈Op)(16)
如果识别为一新工况,则建立新运行工况子空间,训练相应的新健康状态评估模型,并保存。
获得各部件健康指数后,利用下面的变权综合公式,最终得到整个齿轮箱的健康指数:
其中:
HI为ti时刻齿轮箱的健康指数;
为第j个部件健康状态评价时占的常权权重,j=1,2,…m;
为ti时刻第j个部件健康指数的γ次幂;
为第k个部件健康状态评价时占的常权权重,k=1,2,…m;
为ti时刻第k个部件健康指数的γ-1次幂;
γ为均衡系数,本发明取γ=0.1
风电机组通常设置在风力强的偏远地区,当前通常采用计划维修与事后维修相结合的维护策略。但这种维护策略存在以下问题:(1)计划维修是在停机状态下进行的检测维修,不能反映齿轮箱运行时的状况;计划维修时间间隔太长,即便运行中出现问题也不能及时发现。(2)事后维修一般都比较突然,缺乏应有的准备工作,需要较长的维修时间和较大的维修费用。
通过对齿轮箱健康状态进行实时动态评价并根据评价结果开展健康状态衰退趋势分析,可实现齿轮箱故障的早期预报。根据分析结果相应地进行运行调整或维修安排,最终达到提高机组可用率,降低运行与维护费用的目的。
下面结合附图说明具体实施方式。
本发明采用目前风电机组常用的一级行星两级圆柱齿轮箱为研究对象,如图2所示。整个齿轮箱分为低速级、中速级和高速级。其中,202为环形内齿圈,204为行星齿轮(共三个),206为太阳齿轮;208为中速级齿轮,210为中速级小齿轮;212为高速级齿轮,214为高速级小齿轮。当然本发明同样适用于其它形式的传动结构(如两级行星一级斜齿等)。
风电机组齿轮箱健康状态动态监控总体流程如图1所示,主要由离线健康状态评价模型训练和在线健康状态动态评估两部分组成。
离线部分:对齿轮箱振动和其它监测历史数据进行信号处理,获得原始状态特征序列;对风机运行工况历史数据进行基于聚类分析的工况辨识,并根据辨识结果将原始状态特征序列进行分割;对原始状态特征进行基于流形学习的特征提取,并作为用于健康状态评价的状态特征;针对不同工况子空间,利用相应的状态特征序列训练健康状态评估模型,并保存于模型库。
在线部分:对实时采集的振动和其它监测数据进行信号/数据处理,生成原始特征,调用离线特征提取结果,在线计算健康状态特征;利用运行工况实时数据,进行工况辨识,确定当前工况所属工况子空间;根据所属子空间,调用模型库相应的健康状态评估模型,进行健康状态的在线评估。
最后根据当前健康状态和历史健康状态衰退趋势,给出运行与维护建议。
1、数据采集
齿轮箱振动测点布置方案如图2所示,对应位置及测量值描述列于表1。
表1齿轮箱振动测点位置及相应测量参数
序号 | 位置 | 测量值描述 | 序号 | 位置 | 测量值描述 |
① | 环形齿轮12点 | 垂直向加速度 | ⑤ | 中速轴电机侧 | 垂直、水平加速度 |
② | 环形齿轮6点 | 垂直向加速度 | ⑥ | 高速轴电机侧 | 垂直、水平加速度 |
③ | 行星齿轮支架 | 垂直向加速度 | ⑦ | 高速轴轴承(迎风侧) | 垂直向加速度 |
④ | 低速轴电机侧 | 垂直、水平加速度 | ⑧ | 高速轴轴承(背风侧) | 垂直向加速度 |
通过振动数据采集装置实现振动数据的采集,为齿轮箱的健康状态实时动态监控提供数据。
运行工况数据和运行监测数据(如风速、转速、功率、发电机扭矩、空气密度、轴承温度、润滑油温等)通过开发与风电机组SCADA系统的接口程序来取得;风机设计数据由制造厂设计部门提供。
2、信号处理及分析
利用上文提到的时、频域信号分析技术,分别比较齿轮箱主要部件(轴、齿轮、轴承)正常状态和不正常状态下各对应特征参数值,并进而初步确定用于齿轮箱各部分健康状态评价的特征参数,形成原始特征库。本例采用原始特征包括:原始振动时间序列的RMS,振动谱峭度特征均值,信号时域同步平均分析峰峰值,信号时域同步平均分析能量算子,轴承外圈、内圈、滚动体故障特征频率频谱特征,齿轮啮合频率(1X、2X、3X)频谱特征,包络谱轴承外圈、内圈、滚动体频率频谱特征,轴承温度,润滑油温。
3、工况辨识、特征提取及健康评价模型训练
工况辨识、特征提取及健康评价模型训练总体流程如图3所示。一方面对振动和运行监测历史数据进行信号/数据处理,生成原始特征序列;另一方面对运行工况历史数据(风速、空气密度、输出功率、发电机扭矩、轴转速)开展基于聚类分析的工况辨识,将风电机组运行工况空间划分为若干子工况空间;根据工况划分将振动特征和运行监测历史数据划分入对应工况子空间;对各工况子空间,进行特征提取,并计算出相应状态特征值序列并保存;最后利用各子空间状态特征值序列训练相应健康状态评价模型并存入模型库。
针对运行工况历史数据,利用k-means聚类算法进行工况辨识,公式如下:
其中:
μi为聚类中心;
Ω=(O1,O2,…Oj,…OP)为工况空间;
ui为运行工况数据聚类样本,i=1,2,…N。
对各部件原始特征参数,进行基于保局投影(LPP)的特征提取。与PCA提取最具代表性的特征不同,LPP提取最具有判别性的特征。通过最小化下面的目标函数,实现从高维输入数据X=(x1,x2,…,xN)到低维子空间Y=(y1,y2,…,yN)的转变。
其中:
n为样本个数;
yj为时刻tj齿轮箱相应部件特征提取后状态特征向量,yj=ATxj;
S(i,j)为成近邻的第i个样本和第j个样本之间的相似性程度,可以取用k近邻或近邻的高斯权或均匀权作为相似性度量。
经过简单的数学推导,有:
其中:
X=(x1,x2,…,xN)为训练数据集;
L为拉普拉斯矩阵,L=D-S,D是对角权矩阵,它的元素就是S的行(或列,因为S是对称矩阵)元素的数值之和,Dii=∑Sji。
为去除上述最小化问题的一个缩放问题,并消除平移的任意性,可以加上下面的一个限制:
aTXDXTa=1(22)
最终的最小化问题转化为:
也即广义特征值问题:
XLXTa=λXDXTa(24)
假定a0,a1,…,al-1是上式的特征向量解,对应的特征值为λ0<λ1<…,λl-1,则最终得到的线性降维映射可以表示为:
xi→yi=ATxi(25)
其中,A=(a0,a1,…,al-1)是n×l维的变换矩阵。
利用运行监控历史数据训练健康状态评价模型。用于训练的数据按下面的方法选取:对该齿轮箱(或使用条件相似的同型齿轮箱)一个运行周期(从正常投入运行到故障)内的运行监控历史数据进行划分;取其前5%和后5%运行时间对应的运行监控历史数据作为模型训练的备用数据;将这部分备用数据再按辨识出的工况子空间进行划分,并作为对应工况健康状态评价模型的训练数据;计算各工况训练数据的特征值,利用回归算法确定模型参数。
本例的健康状态评价模型采用有监督学习的Logistics回归模型,对应的健康指数可按下式计算:
其中:
为ti时刻对应第p个运行工况子空间的齿轮箱部件的健康指数;
为第p个工况子空间对应的健康状态评估函数;
Θ(p)为第p个工况子空间对应健康状态评价模型的L+1个模型参数,
为健康状态评价模型的第一个模型参数;
为健康状态评价模型第L+1个模型参数;
yil为ti时刻第l个健康状态特征参数值;
ε为残差。
第p个工况子空间的健康评价模型训练规则如下:
(1)选取齿轮箱相应部件早期健康运行时的样本数据集并将输出设定为1,即zi=1;
(2)选取齿轮箱相应部件寿命末期故障运行时的样本数据集并将输出设定为0,即zi=0;(3)齿轮箱部件寿命中期运行时的样本数据集不作为模型训练数据。
本发明也可以采用神经网络模型或只需要健康数据进行模型训练的高斯混合模型(GMM)。
4、健康状态实时动态评估
齿轮箱健康状态实时评估流程如图4所示。
针对某一时刻运行的齿轮箱,对其运行和振动监测数据进行实时分析,计算对应状态特征值。
对其当前运行工况数据进行分析,辨识出当前的运行工况;如果该工况属于已有某一工况子空间,则调用相应的训练好的各部件状态评估模型评价当前健康状态;如果该工况为一新工况,则需要训练新的相应工况的健康状态评价模型,并保存入模型库。
在获得齿轮箱各部件(行星齿轮、低速齿轮轴、中速齿轮轴、高速齿轮轴)的健康指数后,按照图5建立的齿轮箱结构层次,利用式(17)获得整个齿轮箱的健康指数,为下一步的运行与维护决策提供依据。
健康状态评价模型需要动态更新,当验证偏差较大时,则需要进行模型修正或重新训练。
5、运行与维护建议
根据实时、动态评估出的齿轮箱(部件)健康状态和历史健康状态衰退趋势,进行综合分析,并给出运行与维修建议。
序号 | 健康状态指数 | 衰退趋势 | 运行与维修建议 |
1 | z≥0.75 | 平稳(渐变) | 健康状态良好,可以长时间持续运行 |
2 | z≥0.75 | 急剧(突变) | 受到冲击或某部件损坏,应切换为低产能模式观察 |
3 | 0.75>z>0.25 | 平稳(渐变) | 健康状态进入劣化阶段,应加强监视 |
4 | 0.75>z>0.25 | 急剧(突变) | 健康状态衰退加速,加强监视并提前准备检修 |
5 | z≤0.25 | 平稳(渐变) | 健康状态严重劣化,应提前安排检修 |
6 | z≤0.25 | 急剧(突变) | 健康状态严重劣化且衰退加剧,应立即停机检修 |
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (2)
1.一种风电机组齿轮箱健康状态实时动态监控方法,该方法首先采集齿轮箱振动信号和齿轮箱运行与监测数据;其特征是该方法包括以下步骤:
步骤1:对齿轮箱振动信号进行分析,得到齿轮箱振动信号的时域特征、频域特征和包络信号的包络谱特征;所述齿轮箱振动信号频域特征为滚动轴承故障特征频率、齿轮啮合频率处的能量和包络谱各故障特征频率对应的能量;
步骤2:建立齿轮箱运行工况空间,根据齿轮箱运行与监测数据将齿轮箱运行工况空间划分为设定个齿轮箱运行工况子空间;
步骤3:在步骤1和步骤2的基础上,将齿轮箱时域特征、频域特征、运行与监测数据划分为运行监控历史数据子集;
步骤4:对运行监控历史数据子集进行特征提取,生成齿轮箱各部件的健康状态特征库;
步骤5:根据齿轮箱各部件的健康状态特征库建立健康状态评价模型,通过运行监控历史数据子集对健康状态评价模型进行训练,得到各工况的健康状态实时评价模型,通过健康状态实时评价模型对实时的齿轮箱运行与监测数据进行评价,得到各部件的健康指数;
具体为:针对任一运行工况子空间Op,利用相应运行监控历史数据子集X(p)进行健康状态评价模型的训练,得到对应工况的健康状态实时评价模型,则齿轮箱部件的健康指数按下式计算:
其中:
为ti时刻对应第p个运行工况子空间的齿轮箱部件的健康指数;
为第p个工况子空间对应的健康状态评估函数;
yi为时刻ti齿轮箱相应部件特征提取后状态特征向量;
Θ(p)为第p个工况子空间对应健康状态评价模型的L+1个模型参
数,
为健康状态评价模型的第一个模型参数;
为健康状态评价模型第L+1个模型参数;
yil为ti时刻第l个健康状态特征参数值;
ε为残差;
整个齿轮箱的健康状态由其相应的主要组成部件的健康状态决定,而各组成部件的健康状态通过运行工况辨识、振动信号分析、特征提取及健康状态评估确定,建立相应层次结构后实现从底层到顶层的健康状态逐级评价,针对齿轮箱某一运行时刻ti,通过健康状态实时评价模型对其运行工况数据、运行实时数据和状态监测数据进行信号处理和特征计算,得到该时刻运行工况向量和第j个部件状态特征向量(ui,yij),ui为时刻ti运行工况特征向量;yij为时刻ti对应的第j个部件变换后状态特征向量;对其运行工况参数进行工况辨识,如果识别为某一已知工况,则按下式计算第j个部件的健康指数zij:
其中,Cp为某运行工况隶属于第p个齿轮箱运行工况子空间的隶属度,
Cp=Member(ui∈Op);
P’为齿轮箱运行工况子空间个数;
如果识别为一新工况,则建立新运行工况子空间,训练相应的新健康状态评价模型并保存;
步骤6:将各部件的健康指数带入下面的变权综合公式,得到整个齿轮箱的健康指数,则齿轮箱的健康指数按下式计算:
其中:
HI为ti时刻齿轮箱的健康指数;
为第j个部件健康状态评价时占的常权权重,j=1,2,…m;
为ti时刻第j个部件健康指数的γ次幂;
为第k个部件健康状态评价时占的常权权重,k=1,2,…m;
为ti时刻第k个部件健康指数的γ-1次幂;
γ为均衡系数。
2.根据权利要求1所述的一种风电机组齿轮箱健康状态实时动态监控方法,其特征是所述齿轮箱振动信号时域特征为样本时间序列平均值、样本时间序列有效值、峭度指标、峰值指标、歪度指标和熵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210215628.4A CN102768115B (zh) | 2012-06-27 | 2012-06-27 | 一种风电机组齿轮箱健康状态实时动态监控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210215628.4A CN102768115B (zh) | 2012-06-27 | 2012-06-27 | 一种风电机组齿轮箱健康状态实时动态监控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102768115A CN102768115A (zh) | 2012-11-07 |
CN102768115B true CN102768115B (zh) | 2016-04-20 |
Family
ID=47095596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210215628.4A Expired - Fee Related CN102768115B (zh) | 2012-06-27 | 2012-06-27 | 一种风电机组齿轮箱健康状态实时动态监控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102768115B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106934242A (zh) * | 2017-03-16 | 2017-07-07 | 杭州安脉盛智能技术有限公司 | 基于交叉熵法的多模式下设备的健康度评估方法及系统 |
CN110410282A (zh) * | 2019-07-24 | 2019-11-05 | 河北工业大学 | 基于som-mqe和sfcm的风电机组健康状态在线监测及故障诊断方法 |
Families Citing this family (69)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033359B (zh) * | 2012-12-19 | 2015-10-28 | 西安交通大学 | 一种多特征多级综合评判的风电机组主传动装置故障诊断方法 |
CN103217286A (zh) * | 2013-03-23 | 2013-07-24 | 中国水利电力物资有限公司 | 基于故障数据的风电机组传动系统故障识别方法和系统 |
CN103226651A (zh) * | 2013-03-23 | 2013-07-31 | 中国水利电力物资有限公司 | 基于相似度统计的风电机组状态评估预警方法和系统 |
CN103234753A (zh) * | 2013-04-11 | 2013-08-07 | 华北电力大学 | 基于高斯过程建模的风电机组轴系状态监测方法 |
CN103398846B (zh) * | 2013-08-16 | 2016-06-22 | 大连美恒时代科技有限公司 | 一种减速器健康分析方法及其分析平台系统 |
CN103983452B (zh) * | 2014-04-18 | 2016-09-07 | 中国人民解放军国防科学技术大学 | 利用混合域特征向量和灰色关联分析对行星齿轮箱进行故障模式识别的方法 |
CN104122086B (zh) * | 2014-08-11 | 2016-05-04 | 长沙理工大学 | 基于Kriging模型的风力机齿轮箱故障诊断方法 |
CN104408298B (zh) * | 2014-11-17 | 2017-10-03 | 同济大学 | 一种基于多重分组迭代区间分支定界技术的ddpc伺服系统故障预示方法 |
ES2947817T3 (es) | 2014-11-18 | 2023-08-21 | Hitachi Energy Switzerland Ag | Método y sistema de supervisión del estado de un aerogenerador |
CN104614179B (zh) * | 2015-01-28 | 2017-08-29 | 中国大唐集团新能源股份有限公司 | 一种风电机组齿轮箱状态监测方法 |
CN104807644B (zh) * | 2015-04-14 | 2019-04-05 | 北京中恒博瑞数字电力科技有限公司 | 一种风电机组变桨系统的故障预警方法及系统 |
CN105162413A (zh) * | 2015-09-08 | 2015-12-16 | 河海大学常州校区 | 一种基于工况辨识的光伏系统性能实时评估方法 |
CN105115594B (zh) * | 2015-10-09 | 2017-04-05 | 北京航空航天大学 | 基于小波熵和信息融合的齿轮箱振动信号故障特征提取方法 |
CN106228176A (zh) * | 2016-06-29 | 2016-12-14 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种基于工况数据检测设备状态的方法和系统 |
CN106446540B (zh) * | 2016-09-20 | 2019-05-03 | 华北电力大学(保定) | 一种风电机组健康状态实时评估方法 |
CN108009174B (zh) * | 2016-10-28 | 2021-06-01 | 沈阳高精数控智能技术股份有限公司 | 一种基于模式匹配的振动事件时间顺序区分方法 |
CN106682814B (zh) * | 2016-11-28 | 2021-01-01 | 华北电力大学 | 一种基于故障知识库的风电机组故障智能诊断方法 |
CN106681151B (zh) * | 2017-01-20 | 2019-05-10 | 华北电力大学 | 基于时序分析的变桨距风机多模态运行优化控制装置及运行方法 |
CN106952028A (zh) * | 2017-03-13 | 2017-07-14 | 杭州安脉盛智能技术有限公司 | 机电装备故障预诊与健康管理方法及系统 |
CN107016235B (zh) * | 2017-03-21 | 2020-06-19 | 西安交通大学 | 基于多特征自适应融合的设备运行状态健康度评估方法 |
EP3388910A1 (en) * | 2017-04-10 | 2018-10-17 | ABB Schweiz AG | Method and apparatus for monitoring the condition of subsystems within a renewable generation plant or microgrid |
CN107101827B (zh) * | 2017-06-19 | 2019-05-21 | 苏州微著设备诊断技术有限公司 | 一种低速重载齿轮裂纹故障在线检测方法 |
CN107831012B (zh) * | 2017-10-11 | 2019-09-03 | 温州大学 | 一种基于Walsh变换与Teager能量算子的轴承故障诊断方法 |
CN107832896B (zh) * | 2017-11-29 | 2021-03-12 | 广东电网有限责任公司电力科学研究院 | 一种电厂设备缓变故障预警方法及装置 |
CN108204895A (zh) * | 2017-11-29 | 2018-06-26 | 中国电力科学研究院有限公司 | 一种风电机组变速齿轮箱故障诊断方法及系统 |
CN108036940A (zh) * | 2017-12-04 | 2018-05-15 | 中国水利水电科学研究院 | 一种轴承故障诊断方法 |
CN108062603A (zh) * | 2017-12-29 | 2018-05-22 | 国网福建省电力有限公司 | 基于配电自动化终端设备全生命周期寿命预测方法及系统 |
CN108377209A (zh) * | 2018-01-17 | 2018-08-07 | 北京和利时智能技术有限公司 | 基于scada的设备故障检测系统和检测方法 |
CN108375729B (zh) * | 2018-02-02 | 2019-03-08 | 北京航空航天大学 | 基于Fisher判别的电机SOM聚类退化状态评估方法 |
CN108549220A (zh) * | 2018-03-29 | 2018-09-18 | 广东电网有限责任公司电力调度控制中心 | 燃煤机组运行状态实时评估方法及其系统 |
CN108595381B (zh) * | 2018-04-27 | 2022-03-22 | 厦门尚为科技股份有限公司 | 健康状态评价方法、装置及可读存储介质 |
CN109141881B (zh) * | 2018-07-06 | 2020-03-31 | 东南大学 | 一种深度自编码网络的旋转机械健康评估方法 |
CN109051618B (zh) * | 2018-09-06 | 2020-03-31 | 广东寰球智能科技有限公司 | 基于振动信号处理传送星轮故障的监测方法及监测系统 |
CN109357753B (zh) * | 2018-11-07 | 2021-01-19 | 深圳中广核工程设计有限公司 | 降低核电厂的松脱部件与振动监测误报警率的系统和方法 |
CN109376794A (zh) * | 2018-11-12 | 2019-02-22 | 武汉科技大学 | 一种铁水脱硫搅拌器的智能故障诊断方法及系统 |
CN111379624A (zh) * | 2018-12-29 | 2020-07-07 | 中国葛洲坝集团装备工业有限公司 | 一种多工况和时间深度并行诊断方法 |
CN110069814A (zh) * | 2019-03-13 | 2019-07-30 | 中交广州航道局有限公司 | 船机齿轮箱的状态预测方法、装置及计算机设备 |
CN110196157B (zh) * | 2019-04-08 | 2020-12-22 | 华电电力科学研究院有限公司 | 一种基于转机设备全息诊断远程分析系统的诊断方法 |
CN110044624A (zh) * | 2019-04-24 | 2019-07-23 | 华东交通大学 | 一种齿轮箱轴承故障检测方法及系统 |
CN110187275B (zh) * | 2019-06-06 | 2021-11-23 | 中车株洲电力机车研究所有限公司 | 一种永磁电机健康状态检测方法及系统 |
CN110428020A (zh) * | 2019-08-13 | 2019-11-08 | 树根互联技术有限公司 | 刀具健康评估模型训练方法、装置及电子设备 |
CN112541597A (zh) * | 2019-09-04 | 2021-03-23 | 上海杰之能软件科技有限公司 | 多设备维修方法及装置、存储介质、终端 |
CN110736926A (zh) * | 2019-10-14 | 2020-01-31 | 湖南银河电气有限公司 | 电机运行状态的特征参数提取方法和装置 |
CN110765699B (zh) * | 2019-10-17 | 2021-07-30 | 中国石油大学(北京) | 一种压裂装备作业时健康状态的评估方法及装置 |
TWI749416B (zh) * | 2019-11-29 | 2021-12-11 | 中國鋼鐵股份有限公司 | 變轉速設備異常監診方法 |
CN111061257B (zh) * | 2019-12-30 | 2021-02-19 | 杭州电子科技大学 | 一种基于动态全局lpp的工业过程监测方法 |
CN111308337A (zh) * | 2020-03-31 | 2020-06-19 | 张铭源 | 一种离心风机入口导叶调节阀门性能评价方法 |
CN111551361A (zh) * | 2020-03-31 | 2020-08-18 | 河北科技大学 | 行星齿轮传动系统故障状态识别方法 |
CN111307461A (zh) * | 2020-04-13 | 2020-06-19 | 郑州恩普特科技股份有限公司 | 一种基于特征向量基线法的滚动轴承故障检测方法 |
CN113723732A (zh) * | 2020-05-25 | 2021-11-30 | 中国石油化工股份有限公司 | 用于离心泵的状态判定方法及系统 |
WO2022036689A1 (en) * | 2020-08-21 | 2022-02-24 | Abb Schweiz Ag | Method for fault diagnosis using a fault diagnosis model trained by virtual data |
CN112177863A (zh) * | 2020-09-07 | 2021-01-05 | 上海电气风电集团股份有限公司 | 振动监控系统、风力发电系统及风电场 |
CN112149301B (zh) * | 2020-09-23 | 2023-05-16 | 创新奇智(上海)科技有限公司 | 模型训练方法、故障预测方法、装置及电子设备 |
CN112539827B (zh) * | 2020-12-04 | 2021-10-01 | 五凌电力有限公司 | 基于频率能量占比的机组机械振动评估方法和系统 |
CN112696481A (zh) * | 2020-12-11 | 2021-04-23 | 龙源(北京)风电工程技术有限公司 | 风电机组齿轮箱轴温异常智能诊断方法及装置 |
CN112729815A (zh) * | 2020-12-21 | 2021-04-30 | 云南迦南飞奇科技有限公司 | 基于无线网络的输送线健康状况在线故障大数据预警方法 |
CN113065097B (zh) * | 2021-03-18 | 2022-01-04 | 湖南南方宇航高精传动有限公司 | 一种行星齿轮箱均载测试载荷系数的计算方法 |
CN113051814B (zh) * | 2021-03-18 | 2023-11-24 | 王嘉寅 | 一种转动设备性能劣化评估方法、装置及存储介质 |
CN113052555A (zh) * | 2021-03-26 | 2021-06-29 | 浙江三一装备有限公司 | 损耗零部件管理方法及系统 |
CN113295419A (zh) * | 2021-05-26 | 2021-08-24 | 浙江运达风电股份有限公司 | 一种风电机组齿轮箱中速级轴承故障预警方法 |
CN113255777B (zh) * | 2021-05-28 | 2023-01-31 | 郑州轻工业大学 | 基于多模态敏感特征选取融合的设备故障预警方法及系统 |
CN113586505B (zh) * | 2021-07-01 | 2024-02-13 | 首钢京唐钢铁联合有限责任公司 | 一种应用于高炉鼓风机的停机控制方法及装置 |
CN113982850B (zh) * | 2021-09-10 | 2023-11-10 | 大唐国信滨海海上风力发电有限公司 | 融合高低频信号的风机综合健康分析方法及系统 |
CN113987713B (zh) * | 2021-11-04 | 2022-05-20 | 华北电力大学 | 一种风电机组齿轮箱螺旋角优化方法及系统 |
CN114033736B (zh) * | 2021-12-23 | 2024-05-31 | 国家能源集团泰州发电有限公司 | 一种基于气压脉动信号风机故障监测系统 |
CN114638280B (zh) * | 2022-01-30 | 2023-04-11 | 江苏东佳电气有限公司 | 基于本地设备网络的防爆电机轴承温度异常监测系统 |
CN115235762A (zh) * | 2022-09-26 | 2022-10-25 | 江苏东控自动化科技有限公司 | 冶金传动机构局域损伤振动包络信号获取及评估方法 |
CN115758260B (zh) * | 2023-01-05 | 2023-07-04 | 青岛明思为科技有限公司 | 基于高斯混合模型的机械设备状态检测的方法 |
CN117972547B (zh) * | 2024-03-26 | 2024-06-14 | 华电电力科学研究院有限公司 | 一种风电机组关键部件的故障预警方法、装置及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101634605A (zh) * | 2009-04-10 | 2010-01-27 | 北京工业大学 | 基于混合推理和神经网络的齿轮箱故障智能诊断方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012098213A (ja) * | 2010-11-04 | 2012-05-24 | Nsk Ltd | 転がり直動要素用稼働状態監視装置および転がり直動要素の稼働状態監視方法 |
-
2012
- 2012-06-27 CN CN201210215628.4A patent/CN102768115B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101634605A (zh) * | 2009-04-10 | 2010-01-27 | 北京工业大学 | 基于混合推理和神经网络的齿轮箱故障智能诊断方法 |
Non-Patent Citations (4)
Title |
---|
Logistic回归模型及其应用;常振海等;《延边大学学报(自然科学版)》;20120331;第38卷(第1期);第28-32页 * |
基于EMD的齿轮箱故障特征信息提取研究;王春林;《中国优秀硕士学位论文全文数据库工程科技II辑》;20120215(第2期);第19-21页,第43-66页 * |
基于支持向量机的风电机组运行工况分类方法;刘永前等;《太阳能学报》;20100930;第31卷(第9期);第1191页最后1段至第1195页最后1段 * |
大型汽轮机组变权综合状态评价模型研究;董玉亮等;《华北电力大学学报》;20050331;第32卷(第2期);第47页第1节至第48页第3节 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106934242A (zh) * | 2017-03-16 | 2017-07-07 | 杭州安脉盛智能技术有限公司 | 基于交叉熵法的多模式下设备的健康度评估方法及系统 |
CN110410282A (zh) * | 2019-07-24 | 2019-11-05 | 河北工业大学 | 基于som-mqe和sfcm的风电机组健康状态在线监测及故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102768115A (zh) | 2012-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102768115B (zh) | 一种风电机组齿轮箱健康状态实时动态监控方法 | |
CN103033359B (zh) | 一种多特征多级综合评判的风电机组主传动装置故障诊断方法 | |
CN106682814B (zh) | 一种基于故障知识库的风电机组故障智能诊断方法 | |
CN105787584B (zh) | 一种基于云平台的风电机群故障预警方法 | |
CN103631681B (zh) | 一种在线修复风电场异常数据的方法 | |
CN109255333B (zh) | 一种大型风电机组滚动轴承故障混合诊断方法 | |
CN105205569B (zh) | 风机齿轮箱状态在线评估模型建立方法及在线评估方法 | |
CN202814565U (zh) | 风机齿轮箱温度监控装置 | |
CN110766175A (zh) | 一种基于最优间隔分布机的变桨系统故障检测方法及装置 | |
Miao et al. | Condition multi-classification and evaluation of system degradation process using an improved support vector machine | |
Zhong et al. | Multi-fault rapid diagnosis for wind turbine gearbox using sparse Bayesian extreme learning machine | |
CN104819846A (zh) | 一种基于短时傅里叶变换和稀疏层叠自动编码器的滚动轴承声音信号故障诊断方法 | |
Zhan et al. | Abnormal vibration detection of wind turbine based on temporal convolution network and multivariate coefficient of variation | |
CN107907327A (zh) | 一种风电机组行星齿轮箱故障诊断方法 | |
CN111091236A (zh) | 一种按桨距角分类的多分类深度学习短期风功率预测方法 | |
CN107403015A (zh) | 一种基于时间序列相似性的短期光功率预测方法 | |
CN113255458A (zh) | 基于多视图关联特征学习的轴承故障诊断方法 | |
CN110595778A (zh) | 一种基于mmf与igra风电机组轴承故障诊断方法 | |
CN103852255A (zh) | 基于神经网络风力发电机组典型传动故障智能诊断方法 | |
CN103994199B (zh) | 基于状态监测的齿轮箱维护时间确定方法 | |
CN105046349A (zh) | 一种计及尾流效应的风电功率预测方法 | |
CN108832623B (zh) | 一种物理-统计混合两阶段风电功率预测方法 | |
Mutharasan et al. | Fault detection and power quality analysis of wind turbine system using integrated systems | |
CN114048767A (zh) | 一种风电主控系统故障监测及预警方法 | |
CN113821888B (zh) | 基于周期性冲击特征提取和回声状态网络的振动数据故障诊断方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160420 Termination date: 20190627 |
|
CF01 | Termination of patent right due to non-payment of annual fee |