CN105205569A - 风机齿轮箱状态在线评估模型建立方法及在线评估方法 - Google Patents

风机齿轮箱状态在线评估模型建立方法及在线评估方法 Download PDF

Info

Publication number
CN105205569A
CN105205569A CN201510666344.0A CN201510666344A CN105205569A CN 105205569 A CN105205569 A CN 105205569A CN 201510666344 A CN201510666344 A CN 201510666344A CN 105205569 A CN105205569 A CN 105205569A
Authority
CN
China
Prior art keywords
model
state
degree
normal
normal cloud
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.)
Granted
Application number
CN201510666344.0A
Other languages
English (en)
Other versions
CN105205569B (zh
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.)
Huaqiao University
Original Assignee
Huaqiao University
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 Huaqiao University filed Critical Huaqiao University
Priority to CN201510666344.0A priority Critical patent/CN105205569B/zh
Publication of CN105205569A publication Critical patent/CN105205569A/zh
Application granted granted Critical
Publication of CN105205569B publication Critical patent/CN105205569B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

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

Abstract

本发明涉及一种基于趋势预测-正态云模型的风机齿轮箱状态在线评估模型建立方法,以及一种风机齿轮箱状态在线评估方法。趋势状态分析法在实际应用时,由于齿轮箱状态转换边界具有模糊性和不确定性的特征,如何判断趋势状态的转换一直是制约该方法应用的难题,本发明解决了趋势状态分析法所存在的问题。结果表明,本发明所述的模型能够对齿轮箱的早期缺陷及时预警,达到在线状态监测的目的,有助于防范齿轮箱严重故障的发生,提高风电机组安全性、可靠性和经济性。

Description

风机齿轮箱状态在线评估模型建立方法及在线评估方法
技术领域
本发明涉及电力系统技术领域,更具体地说,涉及一种基于趋势预测-正态云模型的风机齿轮箱状态在线评估模型建立方法,以及一种风机齿轮箱状态在线评估方法。
背景技术
为提高风电机组可用性,减少停机时间,现有技术已经做了大量工作以提高风电机组的可靠性。根据风电机组故障和停机时间统计数据,齿轮箱的平均停机时间为12天/年,虽然齿轮箱的故障率并不高,但是其维修、更换和停机损失的总费用却很高昂,因而对齿轮箱进行实时状态监测,及时发现异常,对于降低齿轮箱运维成本和提高机组可靠性具有重大意义。
现有的风电机组齿轮箱状态监测方法中,应用较为广泛的是在线监测系统(CMS)中振动信号的频谱分析,该方法能够给出明确的齿轮箱异常预警信号,达到了较好的诊断效果,但是用于振动信号测量所需的加速度传感器非常昂贵,导致该方法成本过高。
现代风电机组大多建立了联接机组、气象站和远程计算机的数据采集与监控(SCADA)系统,其特点是无需额外安装传感器,符合成本效益,并且由于SCADA系统能够提供全面的设备监测数据,因而被越来越多地应用于实现齿轮箱的状态监测。其中,趋势状态分析法是齿轮箱状态监测的一种有效方法,其基本思想是基于历史数据建立齿轮箱正常运行时相关指标的预测模型,并利用该模型对齿轮箱的运行参数进行预测,当预测趋势背离实际演化趋势时则认为齿轮箱状态出现了异常,由于其概念清晰、实现方便,因而得到了广泛应用。
趋势状态分析法虽然取得了良好的应用效果,但是在实际应用时,由于齿轮箱状态转换边界具有模糊性和不确定性的特征,因此,如何判断趋势状态的转换一直是制约该方法应用的难题。
通过分析残差的统计特性分析进行诊断,但是其残差阀值仍需运行人员根据经验确定;采用统计过程控制技术分析齿轮箱油温和轴承温度的实际值与估计值的残差,该方法建立在数理统计方法的基础之上,需要大量的统计数据支撑。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于趋势预测-正态云模型的风机齿轮箱状态在线评估模型建立方法,以及一种风机齿轮箱状态在线评估方法。
本发明的技术方案如下:
一种风机齿轮箱状态在线评估模型建立方法,包括如下步骤:
1)基于SCADA历史数据,采用SVM建立齿轮箱各温度指标正常状态下的预测模型,分别将正常时的监测数据、异常时的监测数据作为模型输入,进行温度预测,并计算相对误差序列eijk,其中,i=1,2;j=1,…,q;k=1,…,t,当i=1表示正常状态,i=2表示异常状态,j=1,…,q表示齿轮箱各温度指标,k表示预测样本序号;
2)采用改进无确定度逆向正态云发生器模型,利用所求取的相对误差序列eijk提取正常云和异常云的数字特征(Exij,Enij,Heij),利用所求得的数字特征(Exij,Enij,Heij),构建齿轮箱状态正常云模型和异常云模型;其中,期望值Ex是最能够代表定性概念的点,是这个概念量化的最典型样本;熵En反映定性概念的不确定性,表示在论域空间可以被定性概念接受的取值范围大小,即模糊度;超熵He反映代表定性概念值的样本出现的随机性,揭示了模糊性和随机性的关联;
3)计算正态云模型间关联度kj和贡献度wj,完成在线评估模型建立;其中,贡献度wj表示状态云之间的关联度kj对区分不同定性概念的贡献程度,关联度越大,则其对区分定性概念的贡献越小,反之则越大。
作为优选,步骤1)中,进行温度预测时,当齿轮箱状态正常,通过趋势状态分析法得出的参数预测值与真实值的预测误差,小于当齿轮箱状态异常的参数预测值与真实值的预测误差,进而采用预测误差的变化序列描述齿轮箱趋势状态的演化过程。
作为优选,步骤1)中,进行温度预测时,先分别选取与齿轮箱各温度指标相关性符合预设值的监测指标,以选取的监测指标的当前值及齿轮箱各温度指标的上一时刻值作为SVM模型输入,以齿轮箱各温度指标当前时刻值为SVM模型输出,建立齿轮箱输入轴温度预测模型。
作为优选,步骤1)中,在SVM模型的建模过程中,核函数、惩罚系数C和核宽参数σ决定着模型的性能,采用径向基函数RBF作为核函数,惩罚系数C和核宽参数σ由网格法和十折交叉验证法进行寻优,确定最优参数组合。
作为优选,步骤2)中,正态云模型用数字特征(Exij,Enij,Heij)进行表征,设U是一个用精确数值表示的定量论域,C是论域U的定性概念,若x∈X是定性概念C上的一次随机出现,若同时满足以下两个条件:
其中
x对C的关联度 y ( x ) = exp ( - ( x - E x ) 2 / ( 2 E n ′ 2 ) ) ;
则x在论域U上的分布称为正态云。
作为优选,正态云发生器包括正向正态云发生器、逆向正态云发生器和条件正态云发生器;
正向正态云发生器完成从定性概念到定量值的映射,根据正态云模型数字特征(Exij,Enij,Heij)产生云滴;
逆向正态云发生器是实现从定量值到定性概念的转换,将一定数量的精确数据转换为以数字特征(Exij,Enij,Heij)表示的定性概念;
条件正态云发生器包括X条件正态云发生器、Y条件正态云发生器,其中,X条件正态云发生器中,当已知正态云模型数字特征(Exij,Enij,Heij)后,给定x=x0条件,求取y(x)=y(x0)。
作为优选,通过如下算法过程得到超熵He:
2.1)根据数据样本xi,分别得到样本均值一阶样本绝对中心距A、样本方差S2,公式如下:
X ‾ = 1 n Σ i = 1 n x i ;
A = 1 n Σ i = 1 n | x i - X ‾ | ;
S 2 = 1 n - 1 Σ i = 1 n ( x i - X ‾ ) 2 :
2.2) E x = X ‾ ;
2.3) E n = π 2 × A ;
2.4)若S2-En2<0,则删除m=n×1%个离期望值Ex最近的样本点,在剩余样本中重新计算方差S2,直到S2-En2>0,再跳转步骤2.5);若S2-En2>0直接跳转步骤2.5);
2.5) H e = S 2 - En 2 .
作为优选,步骤3)中的关联度kj的计算过程如下:
若把数x作为一个云滴,y表示数x对正态云模型的关联程度,则称y为云滴关联度,由X条件正态云发生器求取,公式如下:
E n ′ = n o r m r n d ( E n , H e 2 ) ;
y ′ = exp ( - ( x 0 - E x ) 2 / ( 2 E n ′ 2 ) ) ;
y = Σ i = 1 k y i ′ / k ;
在定性概念的正态云模型中,正态云模型C1和C2之间的关联度kj定义为:其中,N和M分别表示区间交集和并集:
N = { ( E x 1 - 3 E n 1 , E x 1 + 3 E n 1 ) ∩ ( E x 2 - 3 E n 2 , E x 2 + 3 E n 2 ) } ;
M = { ( E x 1 - 3 E n 1 , E x 1 + 3 E n 1 ) ∪ ( E x 2 - 3 E n 2 , E x 2 + 3 E n 2 ) } ;
贡献度wj的计算公式为:
一种利用所述的风机齿轮箱状态在线评估模型建立方法进行在线评估的方法,包括如下步骤:
A)获取实时SCADA数据,作为SVM模型的输入,求取齿轮箱各温度指标的相对误差序列,进而计算得到云滴关联度,结合正态云模型间关联度kj和贡献度wj,计算待测试样本对正常正态云模型和异常正态云模型的贴近程度ρ1和ρ2
B)按照最大贴近度的原则确定齿轮箱状态;若ρ1>ρ2,则判定齿轮箱为正常状态,若ρ1≤ρ2,则判定齿轮箱为异常状态。
作为优选,贴近度表示样本对不同定性概念正态云模型的贴近程度,计算公式如下:
ρi=yij×wj,i=1,2;j=1,…,q;
其中,ρi表示定性概念正态云模型贴近度,yij表示状态i的第j个指标的云滴关联度,wj表示第j个指标的贡献度。
本发明的有益效果如下:
本发明所述的基于趋势预测-正态云模型的风机齿轮箱状态在线评估模型建立方法,应用支持向量机法建立齿轮箱运行温度预测模型,并利用该模型分析齿轮箱不同状态时的趋势状态特征,求取齿轮箱运行温度预测相对误差序列,以获取描述齿轮箱趋势状态的量化数据;采用改进无确定度逆向正态云发生器模型,利用所求取的相对误差序列提取正常云和异常云的数字特征,构建齿轮箱状态监测云模型。
计算齿轮箱当前状态对正常状态和异常状态的贴近度,将反映齿轮箱趋势状态的量化数据转换为齿轮箱状态的定性评估。利用某风机齿轮箱的实测SCADA系统数据对所提模型进行了验证,结果表明,该模型能够对齿轮箱的早期缺陷及时预警,达到在线状态监测的目的,有助于防范齿轮箱严重故障的发生,提高风电机组安全性、可靠性和经济性。
趋势状态分析法在实际应用时,由于齿轮箱状态转换边界具有模糊性和不确定性的特征,如何判断趋势状态的转换一直是制约该方法应用的难题,本发明解决了趋势状态分析法所存在的问题。
附图说明
图1是本发明的整体流程图,图中,左半部分为离线训练过程(即在线评估模型建立方法流程),右半部分为在线状态监测过程(即在线评估方法流程);
图2是该风机齿轮箱输入轴温度正常状态时预测模型的预测结果曲线图;
图3是该状态时预测相对误差波动曲线图;
图4是预测模型对系统报警前70分钟的预测结果曲线图;
图5是异常状态时预测相对误差波动曲线图;
图6是正态云模型贴近度变化曲线图;
图7是齿轮箱在线状态评估结果曲线图。
具体实施方式
以下结合附图及实施例对本发明进行进一步的详细说明。
本发明为了解决趋势状态分析法所存在的问题,提供一种基于趋势预测-正态云模型的风机齿轮箱状态在线评估模型建立方法,应用支持向量机法建立齿轮箱运行温度预测模型,并利用该模型分析齿轮箱不同状态时的趋势状态特征,求取齿轮箱运行温度预测相对误差序列,以获取描述齿轮箱趋势状态的量化数据;采用改进无确定度逆向正态云发生器模型,利用所求取的相对误差序列提取正常云和异常云的数字特征,构建齿轮箱状态监测云模型,计算齿轮箱当前状态对正常状态和异常状态的贴近度,将反映齿轮箱趋势状态的量化数据转换为齿轮箱状态的定性评估。
所述的基于趋势预测-正态云模型的风机齿轮箱状态在线评估模型建立方法,包括如下步骤:
1)基于SCADA历史数据,采用SVM建立齿轮箱各温度指标正常状态下的预测模型,分别将正常时的监测数据、异常时的监测数据作为模型输入,进行温度预测,并计算相对误差序列eijk,其中,i=1,2;j=1,2,3,4;k=1,…,t,当i=1表示正常状态,i=2表示异常状态,j=1,2,3,4表示齿轮箱各温度指标,k表示预测样本序号。具体包括以下几个分步骤:
1.1)预测误差状态分析。
基于齿轮箱正常状态时的历史数据建立齿轮箱运行温度的预测模型后,将其用于齿轮箱运行温度预测。当齿轮箱状态正常时,趋势状态分析法能够准确挖掘指标参数的非线性关系,得出的参数预测值可较好拟合真实值,预测误差很小;而若齿轮箱处于异常状态时,其参数内部关系发生变化,已由当前关系转移至其他不确定的关系,导致参数趋势的变化并不在预测算法的“记忆”中,此时参数预测效果必然较差,导致误差值很大,因此可以采用预测误差的变化序列描述齿轮箱趋势状态的演化过程。
常用误差指标包括了相对误差,均方根误差,绝对值误差等。其中,相对误差指标将单个样本的预测值与实际值的绝对误差除以实际值,表示绝对误差占实际值的百分比,从而将误差进行规范化。相对误差比起绝对误差更能反映预测的可信程度[19],比起方根误差计算更简单,因而本发明选取相对误差表征预测误差,其公式如下:
e R E = | y p i - y c i y c i | × 100 % , 1 ≤ i ≤ n ;
其中,ypi表示预测值,yci表示实际值,n表示实际值或真实值的个数。
1.2)输入量和输出量的确定。
本发明以风机齿轮箱为状态评估对象,其SCADA系统采样频率为1次/min;该SCADA系统包含了47项连续监测指标,齿轮箱的监测指标包括齿轮箱输入轴温度,齿轮箱输出轴温度,齿轮箱油温和主轴齿轮箱侧轴承温度;为对这四个指标参数值进行准确预测,需先选取与其相关性较强的监测指标,本发明采用SPSS19.0软件所带的变量相关性分析功能[20],对齿轮箱指标进行相关性分析;以齿轮箱输入轴温度为例,分析结果表明其与齿轮箱输出轴温度、齿轮箱油温、发电机转速、叶轮转速以及主轴齿轮箱侧轴承温度高度相关,Pearson相关系数在[0.8,1.0]范围内;因而以上述相关变量的当前值及其齿轮箱输入轴温度上一时刻值作为SVM模型输入,以齿轮箱输入轴温度作为SVM模型输出,建立齿轮箱输入轴温度预测模型。
1.3)SVM模型性能参数选取和训练。
在SVM模型的建模过程中,核函数、惩罚系数C和核宽参数σ决定着模型的性能;由于齿轮箱各温度指标的预测模型是高度非线性模型,因此本发明采用径向基函数(RBF)作为核函数,惩罚系数C和核宽参数σ由网格法和十折交叉验证法进行寻优,确定最优参数组合。
2)采用改进无确定度逆向正态云发生器模型,利用所求取的相对误差序列eijk提取正常云和异常云的数字特征(Exij,Enij,Heij),利用所求得的数字特征(Exij,Enij,Heij),构建齿轮箱状态正常云模型和异常云模型;其中,期望值Ex是最能够代表定性概念的点,是这个概念量化的最典型样本;熵En反映定性概念的不确定性,表示在论域空间可以被定性概念接受的取值范围大小,即模糊度;超熵He反映代表定性概念值的样本出现的随机性,揭示了模糊性和随机性的关联。
正态云模型理论:设U是一个用精确数值表示的定量论域,C是论域U的定性概念,若x∈X是定性概念C上的一次随机出现,若同时满足以下两个条件:
其中
x对C的关联度 y ( x ) = exp ( - ( x - E x ) 2 / ( 2 E n ′ 2 ) ) ;
则x在论域U上的分布称为正态云。
正态云模型通过数字特征(Exij,Enij,Heij)进行表征。
正态云发生器包括正向正态云发生器、逆向正态云发生器和条件正态云发生器;
正向正态云发生器完成从定性概念到定量值的映射,根据正态云模型数字特征(Exij,Enij,Heij)产生云滴;
逆向正态云发生器是实现从定量值到定性概念的转换,将一定数量的精确数据转换为以数字特征(Exij,Enij,Heij)表示的定性概念;
条件正态云发生器包括X条件正态云发生器、Y条件正态云发生器,其中,X条件正态云发生器中,当已知正态云模型数字特征(Exij,Enij,Heij)后,给定x=x0条件,求取y(x)=y(x0)。
改进无确定度逆向正态云发生器求取云数字特征。
在利用逆向正态云发生器生成齿轮箱状态评估云模型时,预测相对误差序列不包括属于定性概念的确定度信息,因此本发明采用无确定度的逆向正态云发生器建模。传统无确定度的逆向正态云发生器生成评估云数字特征时,由于超熵计算过程中如果出现方差S2小于的情况,则会导致超熵He的结果为虚数而无法计算。考虑到熵En已是全局最优解[24],因而采用如下算法改进:
2.1)根据数据样本xi,分别得到样本均值一阶样本绝对中心距A、样本方差S2,公式如下:
X ‾ = 1 n Σ i = 1 n x i ;
A = 1 n Σ i = 1 n | x i - X ‾ | ;
S 2 = 1 n - 1 Σ i = 1 n ( x i - X ‾ ) 2 :
2.2) E x = X ‾ ;
2.3) E n = π 2 × A ;
2.4)若S2-En2<0,则删除m=n×1%个离期望值Ex最近的样本点,在剩余样本中重新计算方差S2,直到S2-En2>0,再跳转步骤2.5);若S2-En2>0直接跳转步骤2.5);
2.5) H e = S 2 - En 2 .
齿轮箱定性概念评估策略:若把数x作为一个云滴,y表示数x对正态云模型的关联程度,则称y为云滴关联度,其可由X条件正态云发生器求取,具体步骤如下:
公式如下:
E n ′ = n o r m r n d ( E n , H e 2 ) ;
y ′ = exp ( - ( x 0 - E x ) 2 / ( 2 E n ′ 2 ) ) ;
所求云滴关联度y′不是一个数值,而是小范围内的随机数集合;本文计算随机数的平均值,即最终确定云滴关联度y。
在定性概念的正态云模型中,99.74%的云滴位于(Ex-3En,Ex+3En)区间内,若将这个区间视为一个集合,则正态云模型C1和C2之间的关联度kj定义为:其中,N和M分别表示区间交集和并集:
N = { ( E x 1 - 3 E n 1 , E x 1 + 3 E n 1 ) ∩ ( E x 2 - 3 E n 2 , E x 2 + 3 E n 2 ) } ;
M = { ( E x 1 - 3 E n 1 , E x 1 + 3 E n 1 ) ∪ ( E x 2 - 3 E n 2 , E x 2 + 3 E n 2 ) } ;
贡献度wj表示正态云模型间的关联度对区分不同定性概念的贡献程度,关联度越大,则其对区分定性概念的贡献越小,反之则越大;贡献度wj的计算公式为:
基于上述的风机齿轮箱状态在线评估模型建立方法,本发明还提供一种在线评估的方法,如图1所示,包括如下步骤:
A)获取实时SCADA数据,作为SVM模型的输入,求取齿轮箱各温度指标的相对误差序列,进而计算得到云滴关联度,结合正态云模型间关联度kj和贡献度wj,计算待测试样本对正常正态云模型和异常正态云模型的贴近程度ρ1和ρ2
B)按照最大贴近度的原则确定齿轮箱状态;若ρ1>ρ2,则判定齿轮箱为正常状态,若ρ1≤ρ2,则判定齿轮箱为异常状态。
至此,完成了对风机齿轮箱状态评估过程。
贴近度表示样本对不同定性概念正态云模型的贴近程度,计算公式如下:
ρi=yij×wj,i=1,2;j=1,2,3,4;
其中,ρi表示定性概念正态云模型贴近度,yij表示状态i的第j个指标的云滴关联度,wj表示第j个指标的贡献度。
本实例中,采用MATLAB对各具体步骤进行仿真。
(1)正常和异常状态下齿轮箱相关温度指标预测。
获取前文所述风电机组齿轮箱的SCADA系统监测数据,以齿轮箱输入轴温度为例,建立正常状态SVM温度模型,进行温度预测。选取该风机齿轮箱在半个月内的正常运行数据,采用其中400个连续样本作为训练数据,以间隔10个小时以上的200个连续样本作为测试数据,保证了模型的科学性和适用性。
预测结果如图2、图3所示,平均相对误差为0.45%,最大相对误差为1.19%。此时预测平均相对误差为0.45%,最大相对误差为1.19%,模型的预测精度较高,能够准确反映齿轮箱输入轴温度正常状态时的趋势变化。
该风电机组在2012年1月22日2时41分发出齿轮箱故障预警,由于齿轮箱故障的发生是一个量变到质变的缓慢发展过程,SCADA系统能够记录该时间段内的异常监测数据,因而本发明获取此次故障发生前的70个连续样本作为测试数据进行预测,预测结果如图4、图5所示。此时平均相对误差为18.07%,最大相对误差为24.44%。
显然,由于齿轮箱预测模型是根据正常状态的数据所建立,因而当齿轮箱状态发生异常时,其趋势状态变化特征已经与正常状态背离,难以准确描述。
得到一定数量正常和异常状态时的预测样本之后,计算各预测样本的相对误差指标。同样地,根据上述过程能够求得齿轮箱输出轴温度,齿轮箱油温和主轴齿轮箱侧轴承温度在正常和异常状态下的相对误差序列。
(2)计算正态云模型的数字特征。
仍以齿轮箱输入轴温度为例,将正常预测样本的相对误差作为改进无确定度逆向正态云发生器的输入,得到正常正态云模型的数字特性为(Ex1,En1,He1)=(0.45,0.28,0.02)。在超熵He1的计算过程中,若采用传统发生器,得到的结果为虚数0.26i,采用改进发生器剔除22个样本点后,在剩余178个样本点中重新计算得到He=0.02。同理,异常正态云模型的数字特征为(Ex2,En2,He2)=(14.51,5.17,8.53)。表1所示为齿轮箱各指标状态的正态云模型数字特征。
表1
指标 正常云 异常云
齿轮箱输入轴温度 (0.45,0.28,0.02) (13.97,5.36,0.33)
齿轮箱输出轴温度 (0.82,0.62,0.17) (23.92,8.32,1.55)
齿轮箱油温 (0.60,0.48,0.17) (17.39,6.35,0.62)
主轴齿轮箱侧轴承温度 (0.62,0.44,0.56) (18.67,7.42,0.56)
(3)齿轮箱在线状态评估可行性验证。
该风电机组于2012年3月21日14:27发出齿轮箱故障警报,本发明获取了该故障发生前3个小时内的SCADA系统数据,共180组监测数据,应用本发明模型对其进行在线状态评估,
由图6可直观看出,在前半部分,正常的贴近度曲线存在波动,但其值明显大于异常的贴近度。在后半部分,正常贴近度快速下降为0,而异常贴近度迅速上升至较大数值。从图7的在线状态评估结果来看,在第105个样本点处,即故障前76分钟,本发明所提模型评估齿轮箱为异常状态,可提醒机组监控人员及时排查齿轮箱问题。
实例分析的结果验证了本发明所提模型能够根据齿轮箱实测数据在线评估其状态,评估结果符合实际情况,能够对出现的早期缺陷进行及时预警,具有一定的故障早报能力。
上述实施例仅是用来说明本发明,而并非用作对本发明的限定。只要是依据本发明的技术实质,对上述实施例进行变化、变型等都将落在本发明的权利要求的范围内。

Claims (10)

1.一种风机齿轮箱状态在线评估模型建立方法,其特征在于,包括如下步骤:
1)基于SCADA历史数据,采用SVM建立齿轮箱各温度指标正常状态下的预测模型,分别将正常时的监测数据、异常时的监测数据作为模型输入,进行温度预测,并计算相对误差序列eijk,其中,i=1,2;j=1,…,q;k=1,…,t,当i=1表示正常状态,i=2表示异常状态,j=1,…,q表示齿轮箱各温度指标,k表示预测样本序号;
2)采用改进无确定度逆向正态云发生器模型,利用所求取的相对误差序列eijk提取正常云和异常云的数字特征(Exij,Enij,Heij),利用所求得的数字特征(Exij,Enij,Heij),构建齿轮箱状态正常云模型和异常云模型;其中,期望值Ex是最能够代表定性概念的点,是这个概念量化的最典型样本;熵En反映定性概念的不确定性,表示在论域空间可以被定性概念接受的取值范围大小,即模糊度;超熵He反映代表定性概念值的样本出现的随机性,揭示了模糊性和随机性的关联;
3)计算正态云模型间关联度kj和贡献度wj,完成在线评估模型建立;其中,贡献度wj表示状态云之间的关联度kj对区分不同定性概念的贡献程度,关联度越大,则其对区分定性概念的贡献越小,反之则越大。
2.根据权利要求1所述的风机齿轮箱状态在线评估模型建立方法,其特征在于,步骤1)中,进行温度预测时,当齿轮箱状态正常,通过趋势状态分析法得出的参数预测值与真实值的预测误差,小于当齿轮箱状态异常的参数预测值与真实值的预测误差,进而采用预测误差的变化序列描述齿轮箱趋势状态的演化过程。
3.根据权利要求2所述的风机齿轮箱状态在线评估模型建立方法,其特征在于,步骤1)中,进行温度预测时,先分别选取与齿轮箱各温度指标相关性符合预设值的监测指标,以选取的监测指标的当前值及齿轮箱各温度指标的上一时刻值作为SVM模型输入,以齿轮箱各温度指标当前时刻值为SVM模型输出,建立齿轮箱输入轴温度预测模型。
4.根据权利要求3所述的风机齿轮箱状态在线评估模型建立方法,其特征在于,步骤1)中,在SVM模型的建模过程中,核函数、惩罚系数C和核宽参数σ决定着模型的性能,采用径向基函数RBF作为核函数,惩罚系数C和核宽参数σ由网格法和十折交叉验证法进行寻优,确定最优参数组合。
5.根据权利要求1所述的风机齿轮箱状态在线评估模型建立方法,其特征在于,步骤2)中,正态云模型用数字特征(Exij,Enij,Heij)进行表征,设U是一个用精确数值表示的定量论域,C是论域U的定性概念,若x∈X是定性概念C上的一次随机出现,若同时满足以下两个条件:
x ~ N ( E x , E n ′ 2 ) , 其中 E n ′ ~ N ( E x , H e 2 ) ;
x对C的关联度 y ( x ) = exp ( - ( x - E x ) 2 / ( 2 E n ′ 2 ) ) ;
则x在论域U上的分布称为正态云。
6.根据权利要求5所述的风机齿轮箱状态在线评估模型建立方法,其特征在于,正态云发生器包括正向正态云发生器、逆向正态云发生器和条件正态云发生器;
正向正态云发生器完成从定性概念到定量值的映射,根据正态云模型数字特征(Exij,Enij,Heij)产生云滴;
逆向正态云发生器是实现从定量值到定性概念的转换,将一定数量的精确数据转换为以数字特征(Exij,Enij,Heij)表示的定性概念;
条件正态云发生器包括X条件正态云发生器、Y条件正态云发生器,其中,X条件正态云发生器中,当已知正态云模型数字特征(Exij,Enij,Heij)后,给定x=x0条件,求取y(x)=y(x0)。
7.根据权利要求5所述的风机齿轮箱状态在线评估模型建立方法,其特征在于,通过如下算法过程得到超熵He:
2.1)根据数据样本xi,分别得到样本均值一阶样本绝对中心距A、样本方差S2,公式如下:
X ‾ = 1 n Σ i = 1 n x i ;
A = 1 n Σ i = 1 n | x i - X ‾ | ;
S 2 = 1 n - 1 Σ i = 1 n ( x i - X ‾ ) 2 :
2.2) E x = X ‾ ;
2.3) E n = π 2 × A ;
2.4)若S2-En2<0,则删除m=n×1%个离期望值Ex最近的样本点,在剩余样本中重新计算方差S2,直到S2-En2>0,再跳转步骤2.5);若S2-En2>0直接跳转步骤2.5);
2.5) H e = S 2 - En 2 .
8.根据权利要求7所述的风机齿轮箱状态在线评估模型建立方法,其特征在于,步骤3)中的关联度kj的计算过程如下:
若把数x作为一个云滴,y表示数x对正态云模型的关联程度,则称y为云滴关联度,由X条件正态云发生器求取,公式如下:
E n ′ = n o r m r n d ( E n , H e 2 ) ;
y ′ = exp ( - ( x 0 - E x ) 2 / ( 2 E n ′ 2 ) ) ;
y = Σ i = 1 k y i ′ / k ;
在定性概念的正态云模型中,正态云模型C1和C2之间的关联度kj定义为:其中,N和M分别表示区间交集和并集:
N = { ( E x 1 - 3 E n 1 , E x 1 + 3 E n 1 ) ∩ ( E x 2 - 3 E n 2 , E x 2 + 3 E n 2 ) } ;
M = { ( E x 1 - 3 E n 1 , E x 1 + 3 E n 1 ) ∪ ( E x 2 - 3 E n 2 , E x 2 + 3 E n 2 ) } ;
贡献度wj的计算公式为:
9.一种利用权利要求1至8任一项所述的风机齿轮箱状态在线评估模型建立方法进行在线评估的方法,其特征在于,包括如下步骤:
A)获取实时SCADA数据,作为SVM模型的输入,求取齿轮箱各温度指标的相对误差序列,进而计算得到云滴关联度,结合正态云模型间关联度kj和贡献度wj,计算待测试样本对正常正态云模型和异常正态云模型的贴近程度ρ1和ρ2
B)按照最大贴近度的原则确定齿轮箱状态;若ρ1>ρ2,则判定齿轮箱为正常状态,若ρ1≤ρ2,则判定齿轮箱为异常状态。
10.根据权利要求9所述的进行在线评估的方法,其特征在于,贴近度表示样本对不同定性概念正态云模型的贴近程度,计算公式如下:
ρi=yij×wj,i=1,2;j=1,…,q;
其中,ρi表示定性概念正态云模型贴近度,yij表示状态i的第j个指标的云滴关联度,wj表示第j个指标的贡献度。
CN201510666344.0A 2015-10-15 2015-10-15 风机齿轮箱状态在线评估模型建立方法及在线评估方法 Active CN105205569B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510666344.0A CN105205569B (zh) 2015-10-15 2015-10-15 风机齿轮箱状态在线评估模型建立方法及在线评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510666344.0A CN105205569B (zh) 2015-10-15 2015-10-15 风机齿轮箱状态在线评估模型建立方法及在线评估方法

Publications (2)

Publication Number Publication Date
CN105205569A true CN105205569A (zh) 2015-12-30
CN105205569B CN105205569B (zh) 2018-11-27

Family

ID=54953238

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510666344.0A Active CN105205569B (zh) 2015-10-15 2015-10-15 风机齿轮箱状态在线评估模型建立方法及在线评估方法

Country Status (1)

Country Link
CN (1) CN105205569B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105718740A (zh) * 2016-01-24 2016-06-29 辽宁工程技术大学 一种确定系统可靠性运行环境的方法
CN106446540A (zh) * 2016-09-20 2017-02-22 华北电力大学(保定) 一种风电机组健康状态实时评估方法
CN106597574A (zh) * 2016-12-30 2017-04-26 重庆邮电大学 一种基于时变云模型的天气温度预测方法及装置
CN107247230A (zh) * 2017-07-03 2017-10-13 浙江大学 一种基于支持向量机和数据驱动的旋转电机状态监测方法
CN107643178A (zh) * 2017-09-15 2018-01-30 苏州大成有方数据科技有限公司 一种进行齿轮分析的试验台数据管理系统及其工作方法
CN109632293A (zh) * 2018-12-29 2019-04-16 河南理工大学 一种基于云平台的齿轮云测量系统及齿轮云测量方法
CN110032981A (zh) * 2019-04-19 2019-07-19 电子科技大学 基于改进支持向量机的旋转机械故障识别方法
CN111060323A (zh) * 2019-12-05 2020-04-24 山东超越数控电子股份有限公司 一种基于rbf神经网络的航空发动机过振的自动识别方法
CN111075661A (zh) * 2019-12-25 2020-04-28 明阳智慧能源集团股份公司 基于温度变化趋势判断风电机组主轴轴承健康状况的方法
CN111144009A (zh) * 2019-12-27 2020-05-12 广东电科院能源技术有限责任公司 一种风机的运行状态评估方法、装置、设备和存储介质
CN112001632A (zh) * 2020-08-25 2020-11-27 中国船舶重工集团海装风电股份有限公司 一种风电机组齿轮箱轴承性能退化状态评估方法
CN112994047A (zh) * 2021-04-23 2021-06-18 华侨大学 基于二维云模型的自适应分数阶pid负荷频率控制方法
CN113076694A (zh) * 2019-06-20 2021-07-06 北京实力伟业环保科技有限公司 基于物联网的风机故障评价设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103679273A (zh) * 2013-12-20 2014-03-26 南京邮电大学 一种基于隶属云理论的不确定性推理方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103679273A (zh) * 2013-12-20 2014-03-26 南京邮电大学 一种基于隶属云理论的不确定性推理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
郭鹏等: "风电机组齿轮箱温度趋势状态监测及分析方法", 《中国电机工程学报》 *
黄必清: "基于模糊综合评价的海上直驱风电机组运行状态评估", 《清华大学学报(自然科学版)》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105718740A (zh) * 2016-01-24 2016-06-29 辽宁工程技术大学 一种确定系统可靠性运行环境的方法
CN106446540B (zh) * 2016-09-20 2019-05-03 华北电力大学(保定) 一种风电机组健康状态实时评估方法
CN106446540A (zh) * 2016-09-20 2017-02-22 华北电力大学(保定) 一种风电机组健康状态实时评估方法
CN106597574A (zh) * 2016-12-30 2017-04-26 重庆邮电大学 一种基于时变云模型的天气温度预测方法及装置
CN107247230A (zh) * 2017-07-03 2017-10-13 浙江大学 一种基于支持向量机和数据驱动的旋转电机状态监测方法
CN107247230B (zh) * 2017-07-03 2020-02-07 浙江大学 一种基于支持向量机和数据驱动的旋转电机状态监测方法
CN107643178A (zh) * 2017-09-15 2018-01-30 苏州大成有方数据科技有限公司 一种进行齿轮分析的试验台数据管理系统及其工作方法
CN109632293A (zh) * 2018-12-29 2019-04-16 河南理工大学 一种基于云平台的齿轮云测量系统及齿轮云测量方法
CN110032981A (zh) * 2019-04-19 2019-07-19 电子科技大学 基于改进支持向量机的旋转机械故障识别方法
CN110032981B (zh) * 2019-04-19 2022-07-26 电子科技大学 基于改进支持向量机的旋转机械故障识别方法
CN113076694A (zh) * 2019-06-20 2021-07-06 北京实力伟业环保科技有限公司 基于物联网的风机故障评价设备
CN113076694B (zh) * 2019-06-20 2023-10-20 北京实力伟业环保科技有限公司 基于物联网的风机故障评价设备
CN111060323A (zh) * 2019-12-05 2020-04-24 山东超越数控电子股份有限公司 一种基于rbf神经网络的航空发动机过振的自动识别方法
CN111075661A (zh) * 2019-12-25 2020-04-28 明阳智慧能源集团股份公司 基于温度变化趋势判断风电机组主轴轴承健康状况的方法
CN111075661B (zh) * 2019-12-25 2021-11-09 明阳智慧能源集团股份公司 基于温度变化趋势判断风电机组主轴轴承健康状况的方法
CN111144009A (zh) * 2019-12-27 2020-05-12 广东电科院能源技术有限责任公司 一种风机的运行状态评估方法、装置、设备和存储介质
CN112001632A (zh) * 2020-08-25 2020-11-27 中国船舶重工集团海装风电股份有限公司 一种风电机组齿轮箱轴承性能退化状态评估方法
CN112001632B (zh) * 2020-08-25 2022-07-19 中国船舶重工集团海装风电股份有限公司 一种风电机组齿轮箱轴承性能退化状态评估方法
CN112994047A (zh) * 2021-04-23 2021-06-18 华侨大学 基于二维云模型的自适应分数阶pid负荷频率控制方法

Also Published As

Publication number Publication date
CN105205569B (zh) 2018-11-27

Similar Documents

Publication Publication Date Title
CN105205569A (zh) 风机齿轮箱状态在线评估模型建立方法及在线评估方法
CN110926809B (zh) 一种基于大数据分析的风电机组传动链故障预警方法
CN104131950B (zh) 一种风电机组温度特征量的阈值分区确定方法
US20160223600A1 (en) Power generation performance evaluation method and apparatus for power generator set
CN104952000A (zh) 基于马尔科夫链的风电机组运行状态模糊综合评价方法
CN105160489A (zh) 一种水电机组变权劣化评估系统及评估方法
CN103150473A (zh) 一种风电机组发电效率实时监控诊断方法及装置
CN107728059B (zh) 一种变桨系统状态评估方法
CN115013261B (zh) 一种用于海上风电场的状态监测方法及系统
CN104951851A (zh) 一种基于灰关联回归支持向量机的风电机组状态预测模型建立方法
Butler et al. A feasibility study into prognostics for the main bearing of a wind turbine
CN109583075B (zh) 基于温度参数预测的永磁直驱风力机服役质量评价方法
CN111125895B (zh) 一种风力发电机整机服役功率特性在线评估方法
CN103912448A (zh) 一种区域风电场机组功率特性监测方法
CN104122086A (zh) 基于Kriging模型的风力机齿轮箱故障诊断方法
CN104156889A (zh) 一种基于wams数据的风电场性能评估系统及其评估方法
CN112632773B (zh) 一种风电机组可靠性预测方法
WO2023065580A1 (zh) 风电机组齿轮箱的故障诊断方法及装置
de Andrade Vieira et al. Wind turbine condition assessment based on changes observed in its power curve
CN114169718A (zh) 一种基于风电机组状态评估提高机组可靠性的方法
CN117649044A (zh) 一种海上风电机组多维度量化状态评估方法
Zhang et al. Probability warning for wind turbine gearbox incipient faults based on SCADA data
CN115577854A (zh) 一种基于eemd-rbf组合的分位数回归风速区间预测方法
CN114320773B (zh) 一种基于功率曲线分析与神经网络的风电机组故障预警方法
Souza et al. Evaluation of data based normal behavior models for fault detection in wind turbines

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant