CN104976139A - 一种基于高斯模型的机械设备状态诊断方法 - Google Patents

一种基于高斯模型的机械设备状态诊断方法 Download PDF

Info

Publication number
CN104976139A
CN104976139A CN201510312396.8A CN201510312396A CN104976139A CN 104976139 A CN104976139 A CN 104976139A CN 201510312396 A CN201510312396 A CN 201510312396A CN 104976139 A CN104976139 A CN 104976139A
Authority
CN
China
Prior art keywords
gauss model
state
current
normal
equipment
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
CN201510312396.8A
Other languages
English (en)
Other versions
CN104976139B (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.)
Beijing University of Chemical Technology
Original Assignee
Beijing University of Chemical Technology
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 Beijing University of Chemical Technology filed Critical Beijing University of Chemical Technology
Priority to CN201510312396.8A priority Critical patent/CN104976139B/zh
Publication of CN104976139A publication Critical patent/CN104976139A/zh
Application granted granted Critical
Publication of CN104976139B publication Critical patent/CN104976139B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing And Monitoring For Control Systems (AREA)

Abstract

一种基于高斯模型的机械设备状态诊断方法,属于设备故障诊断技术领域的。具体步骤为:(1)根据设备类型确定需要获取的特征参数的种类;(2)获取设备正常状态的特征参数数据并构建对应的高斯模型库;(3)获取设备当前的特征参数数据并构建高斯模型;(4)根据公式计算当前设备状态的高斯模型与正常状态高斯模型的偏离度;(5)根据计算得到的偏离度确定设备当前运行状态。本发明用于对机械设备进行实时的状态诊断,使得设备管理人员掌握设备运行状态。

Description

一种基于高斯模型的机械设备状态诊断方法
技术领域
本发明属于机械设备故障诊断领域,涉及一种基于高斯模型的机械设备状态诊断方法。
背景技术
随着科学技术和现代工业的不断发展和进步,机械设备正朝着大型化、复杂化、高速化等方向发展。这些机械设备对人们日常生活和生产所产生的影响越来越大、越来越明显。设备一旦出现故障,不仅会带来严重的经济损失,还会给人身安全造成威胁。对设备进行实时评估可以较早地发现故障,使得设备管理人员掌握设备运行的状态,这样不仅可以预防事故和减少事故的发生,保证人身和机械设备的安全,还为维修决策提供依据,提高了企业的经济效益和社会效益。
近年来,在设备故障诊断领域,针对不同设备研究了不同的机械设备的状态评估方法。王均华等研究了变电设备状态评估及故障诊断系统,李晋等针对大型机电设备研究了模糊-灰度-变权求和的分层设备健康状态综合评估模型,王孟莲等使用模糊神经网络和支持向量机构建了船舶电力推进系统的状态评估方法,孙志斌等研究了针对六氟化硫封闭式组合电器设备的状态评估方法和检修策略。目前的故障诊断方法研究均是针对不同设备的特点,构建专用的状态评估方法。并未见有文献或专利设计本发明中提到的基于高斯模型的机械设备状态诊断方法。
本发明提出的一种基于高斯模型的机械设备状态诊断方法,获取设备开机正常时间段的特征参数数据,并根据公式构建设备正常状态的高斯模型。在进行设备状态评估的过程中,获取设备当前的特征参数数据构建高斯模型,使用匹配公式计算当前设备状态特征参数高斯模型和设备正常状态高斯模型的匹配度,根据匹配度判断设备当前运行状态正常或异常。
发明内容
本发明的目的是在于实现一种机械设备状态诊断方法,该方法将已知状态评估等级的设备特征参数数据构建成高斯模型,使用模型匹配算法对当前设备的状态进行诊断,最终实现一套能够准确评估设备状态的基于高斯模型的机械设备状态诊断方法。
一种基于高斯模型的机械设备状态诊断方法,其特征在于,包括:
1)获取机械设备开机正常的特征参数数据,生成机械设备正常状态特征参数数据的高斯模型,并将正常状态高斯模型保存到高斯模型库中;
2)实时采集机械设备特征参数数据,生成实时特征参数数据的高斯模型;
3)将实时特征参数数据的高斯模型与状态评估等级及高斯模型库中各个等级特征参数的高斯模型进行匹配计算,得到实时特征参数数据的高斯模型与正常状态高斯模型的偏离度;
4)根据计算得到的高斯模型的偏离度,判定设备运行状态。若计算得到的偏离度超过设定的阈值,则判定设备异常,若不超过设定的阈值,则判定设备运行正常。
根据获取的机械设备一段时间的特征参数数据xi(t),其中i=1,2,…,N,i代表不同种类的特征。不同的设备获取不同数量不同种类的特征参数数据。根据公式(1)(2)(3)构建不同特征参数数据的高斯模型f(xi(t)),
f ( x i ( t ) ) = 1 2 π σ exp ( - ( x i ( t ) - μ ) 2 2 σ 2 ) - - - ( 1 )
其中,
μ = x i ( t 1 ) + x i ( t 2 ) + . . . + x i ( t n ) n - - - ( 2 )
σ = ( x i ( t 1 ) - μ ) 2 + ( x i ( t 2 ) - μ ) 2 + . . . + ( x i ( t n ) - μ ) 2 n - - - ( 3 )
n表示i组特征参数个数,xi(t1)、xi(t2)...xi(tn)表示i组特征参数数据具体数值。
将实时采集得到的当前机械设备的特征参数,根据公式(1)(2)(3)构建当前特征参数数据的高斯模型fcurrent(xi(t));
将当前特征参数数据的高斯模型fcurrent(xi(t))与设备正常状态的高斯模型库fnormal(xi(t))进行匹配识别。模型匹配算法:
P = p 1 + p 2 + . . . + p n n - - - ( 4 )
其中,pi i=1,2,…,n为根据公式(5)计算得到的第i个特征值的偏差度。
p i = 1 - ∫ f current ( x i ( t ) ) f normal ( x i ( t ) ) dx ∫ f current ( x i ( t ) ) dx ∫ f normal ( x i ( t ) ) dx - - - ( 5 )
其中,fcurrent(xi(t))为当前运行状态第i个特征值的高斯模型,fnormal(xi(t))为正常状态第i个特征值的高斯模型。
根据公式(4)(5)计算当前状态的高斯模型fcurrent(xi(t))与设备正常状态高斯模型库fnormal(xi(t))的偏差度P。若计算的到的偏离度p≤Thre,则判定设备运行状态正常;若计算得到的偏离度p>Thre,则判定设备发生故障;其中,Thre为偏离度阈值。偏离度阈值设定为0.05-0.1之间的任意值,建议取0.05。
附图说明
图1一种基于高斯模型的机械设备状态评估方法流程图;
图2第一组多级离心泵工艺特征参数数据趋势
图3第一组多级离心泵振动特征参数数据趋势
图4第二组多级离心泵工艺特征参数数据趋势
图5第二组多级离心泵振动特征参数数据趋势
具体实施方式
下面结合附图和实例对本发明进行具体说明。
一种基于高斯模型的机械设备状态诊断方法流程图如图1所示,主要包括:
1)根据设备的实际情况,由设备管理人员设定需要采集的特征参数数据类型;
2)获取设备开机时间段正常状态下的特征参数数据,并根据公式(1)(2)(3)构建设备正常状态高斯模型;
3)获取设备当前运行状态的特征参数数据,将当前状态的特征参数数据根据公式(1)(2)(3)构建当前状态特征参数数据的高斯模型;
4)根据公式(4)(5)计算当前状态特征参数数据的高斯模型和设备正常状态特征参数数据高斯模型的偏离度p;
5)若计算的到的偏离度p≤Thre,则判定设备运行状态正常;若计算得到的偏离度p>Thre,则判定设备发生故障;
本发明以某多级离心泵试验台为对象对基于高斯模型的机械设备状态诊断方法进行校验。
设定多级离心泵状态评估需要采集的特征参数包括:流量特征参数、入口压力特征参数、出口压力特征参数、多级离心泵驱动端轴承振动特征参数、多级离心泵非驱动端轴承振动特征参数、多级离心泵壳体振动特征参数;
获取多级离心泵4个状态等级的流量、入口压力、出口压力、驱动端轴承振动、非驱动端轴承振动、离心泵壳体振动特征参数,根据公式(1)(2)(3)构建离心泵正常状态的特征参数数据的高斯模型。
分别获取2组离心泵运行一段时间的多级离心泵的6个特征参数数据;
第一组特征参数数据如图2和图3所示;第二组特征参数数据如图4和图5所示。
根据根据公式(1)(2)(3)计算当前状态特征参数数据的高斯模型,根据(4)(5)计算当前特征参数数据高斯模型的偏离度,设定多级离心泵偏离度阈值Thre=0.05。
根据上述公式计算得到第一组特征参数数据的高斯模型偏离度为P=0.0126,判定设备运行状态正常。
根据上述公式计算得到第二组特征参数数据的高斯模型偏离度为:P=0.6231,判定设备运行状态异常。
根据多级离心泵记录,第一组数据是来自离心泵正常运行状态的数据,第二组数据是来离心泵自气蚀故障的数据,该方法诊断准确率100%。

Claims (4)

1.一种基于高斯模型的机械设备状态诊断方法,其特征在于,包括:
1)获取机械设备开机正常的特征参数数据,生成机械设备正常状态特征参数数据的高斯模型,并将正常状态高斯模型保存到高斯模型库中;
2)实时采集机械设备特征参数数据,生成实时特征参数数据的高斯模型;
3)将实时特征参数数据的高斯模型与状态评估等级及高斯模型库中各个等级特征参数的高斯模型进行匹配计算,得到实时特征参数数据的高斯模型与正常状态高斯模型的偏离度;
4)根据计算得到的高斯模型的偏离度,判定设备运行状态;若计算得到的偏离度超过设定的阈值,则判定设备异常,若不超过设定的阈值,则判定设备运行正常。
2.根据权利要求1所述的一种基于高斯模型的机械设备状态诊断方法,其特征在于:所述的特征参数数据高斯模型生成,根据获取的机械设备一段时间的特征参数数据xi(t),其中i=1,2,…,N,i代表不同种类的特征;不同的设备获取不同数量不同种类的特征参数数据;根据公式(1)(2)(3)构建不同特征参数数据的高斯模型f(xi(t)),
f ( x i ( t ) ) = 1 2 π σ exp ( - ( x i ( t ) - μ ) 2 2 σ 2 ) - - - ( 1 )
其中,
μ = x i ( t 1 ) + x i ( t 2 ) + . . . + x i ( t n ) n - - - ( 2 )
σ = ( x i ( t 1 ) - μ ) 2 + ( x i ( t 2 ) - μ ) 2 + . . . + ( x i ( t n ) - μ ) 2 n - - - ( 3 )
n表示i组特征参数个数,xi(t1)、xi(t2)...xi(tn)表示i组特征参数数据具体数值。
3.根据权利要求1所述的一种基于高斯模型的机械设备状态诊断方法,其特征在于:所述的当前状态特征参数数据高斯模型生成,将实时采集得到的当前机械设备的特征参数,根据公式(1)(2)(3)构建当前特征参数数据的高斯模型fcurrent(xi(t))。
4.根据权利要求1所述的一种基于高斯模型的机械设备状态诊断方法,其特征在于:所述的运行状态高斯模型匹配,将当前特征参数数据的高斯模型fcurrent(xi(t))与设备正常状态的高斯模型库fnormal(xi(t))进行匹配识别;模型匹配算法:
P = p 1 + p 2 + . . . + p n n - - - ( 4 )
其中,pi i=1,2,…,n为根据公式(5)计算得到的第i个特征值的偏差度;
p i = 1 - ∫ f current ( x i ( t ) ) f normal ( x i ( t ) ) dx ∫ f current ( x i ( t ) ) dx ∫ f normal ( x i ( t ) ) dx - - - ( 5 )
其中,fcurrent(xi(t))为当前运行状态第i个特征值的高斯模型,fnormal(xi(t))为正常状态第i个特征值的高斯模型;
根据公式(4)(5)计算当前状态的高斯模型fcurrent(xi(t))与设备正常状态高斯模型库fnormal(xi(t))的偏差度P;若计算的到的偏离度p≤Thre,则判定设备运行状态正常;若计算得到的偏离度p>Thre,则判定设备状态异常;其中,Thre为偏离度阈值;偏离度阈值设定为0.05-0.1之间的任意值。
CN201510312396.8A 2015-06-08 2015-06-08 一种基于高斯模型的机械设备状态诊断方法 Active CN104976139B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510312396.8A CN104976139B (zh) 2015-06-08 2015-06-08 一种基于高斯模型的机械设备状态诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510312396.8A CN104976139B (zh) 2015-06-08 2015-06-08 一种基于高斯模型的机械设备状态诊断方法

Publications (2)

Publication Number Publication Date
CN104976139A true CN104976139A (zh) 2015-10-14
CN104976139B CN104976139B (zh) 2017-04-05

Family

ID=54272975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510312396.8A Active CN104976139B (zh) 2015-06-08 2015-06-08 一种基于高斯模型的机械设备状态诊断方法

Country Status (1)

Country Link
CN (1) CN104976139B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105275833A (zh) * 2015-10-30 2016-01-27 北京航空航天大学 一种基于CEEMD-STFT时频信息熵和multi-SVM的离心泵故障诊断方法
CN106446339A (zh) * 2016-08-29 2017-02-22 北京化工大学 基于狄利克雷混合模型的转动机械运行状态异常检测方法
CN106523393A (zh) * 2016-11-18 2017-03-22 山东科技大学 一种用于井下排水系统的故障诊断方法
CN106571016A (zh) * 2016-11-03 2017-04-19 西安交通大学 一种基于报警次数跳变触发机制的机械早期故障判别方法
CN106855054A (zh) * 2015-12-08 2017-06-16 中国石油天然气股份有限公司 管道输油泵状态评价方法及装置
CN107066365A (zh) * 2017-02-20 2017-08-18 阿里巴巴集团控股有限公司 一种系统异常的监测方法及装置
CN107392401A (zh) * 2017-09-08 2017-11-24 上海理工大学 重型轧辊装载智能机器人的故障预测方法
CN109630450A (zh) * 2018-12-03 2019-04-16 中石化重庆天然气管道有限责任公司 一种天然气长输管道离心压缩机组停机保护方法
CN109779938A (zh) * 2018-12-29 2019-05-21 北京博华信智科技股份有限公司 一种离心压缩机组智能联锁保护方法
CN114790991A (zh) * 2022-04-14 2022-07-26 国家电投集团电站运营技术(北京)有限公司 给水泵汽蚀检测系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08189472A (ja) * 1994-10-31 1996-07-23 Meidensha Corp ポンプ吐出量演算方法および装置
CN103234753A (zh) * 2013-04-11 2013-08-07 华北电力大学 基于高斯过程建模的风电机组轴系状态监测方法
GB2500388A (en) * 2012-03-19 2013-09-25 Ge Aviat Systems Ltd Fault detection in system monitoring
CN104595170A (zh) * 2014-12-18 2015-05-06 中国矿业大学 一种自适应核高斯混合模型的空压机监控诊断系统及方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08189472A (ja) * 1994-10-31 1996-07-23 Meidensha Corp ポンプ吐出量演算方法および装置
GB2500388A (en) * 2012-03-19 2013-09-25 Ge Aviat Systems Ltd Fault detection in system monitoring
CN103234753A (zh) * 2013-04-11 2013-08-07 华北电力大学 基于高斯过程建模的风电机组轴系状态监测方法
CN104595170A (zh) * 2014-12-18 2015-05-06 中国矿业大学 一种自适应核高斯混合模型的空压机监控诊断系统及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张栋梁等: "基于状态监测的预警管理系统研究与应用", 《电子设计工程》 *
张龙等: "基于多域特征与高斯混合模型的滚动轴承性能退化评估", 《中国机械工程》 *
钟鑫: "基于逻辑回归和高斯混合模型的设备故障诊断技术研究与应用", 《中国优秀硕士学位论文全文数据库 ENGINEERING SCIENCE AND TECHNOLOGY Ⅱ》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105275833A (zh) * 2015-10-30 2016-01-27 北京航空航天大学 一种基于CEEMD-STFT时频信息熵和multi-SVM的离心泵故障诊断方法
CN105275833B (zh) * 2015-10-30 2017-01-11 北京航空航天大学 一种基于CEEMD-STFT时频信息熵和multi-SVM的离心泵故障诊断方法
CN106855054A (zh) * 2015-12-08 2017-06-16 中国石油天然气股份有限公司 管道输油泵状态评价方法及装置
CN106446339A (zh) * 2016-08-29 2017-02-22 北京化工大学 基于狄利克雷混合模型的转动机械运行状态异常检测方法
CN106446339B (zh) * 2016-08-29 2019-04-30 北京化工大学 基于狄利克雷混合模型的转动机械运行状态异常检测方法
CN106571016A (zh) * 2016-11-03 2017-04-19 西安交通大学 一种基于报警次数跳变触发机制的机械早期故障判别方法
CN106523393A (zh) * 2016-11-18 2017-03-22 山东科技大学 一种用于井下排水系统的故障诊断方法
CN107066365A (zh) * 2017-02-20 2017-08-18 阿里巴巴集团控股有限公司 一种系统异常的监测方法及装置
CN107392401A (zh) * 2017-09-08 2017-11-24 上海理工大学 重型轧辊装载智能机器人的故障预测方法
CN109630450A (zh) * 2018-12-03 2019-04-16 中石化重庆天然气管道有限责任公司 一种天然气长输管道离心压缩机组停机保护方法
CN109779938A (zh) * 2018-12-29 2019-05-21 北京博华信智科技股份有限公司 一种离心压缩机组智能联锁保护方法
CN114790991A (zh) * 2022-04-14 2022-07-26 国家电投集团电站运营技术(北京)有限公司 给水泵汽蚀检测系统及方法

Also Published As

Publication number Publication date
CN104976139B (zh) 2017-04-05

Similar Documents

Publication Publication Date Title
CN104976139A (zh) 一种基于高斯模型的机械设备状态诊断方法
CN107941537B (zh) 一种机械设备健康状态评估方法
CN104131950B (zh) 一种风电机组温度特征量的阈值分区确定方法
CN103033359B (zh) 一种多特征多级综合评判的风电机组主传动装置故障诊断方法
CN104200396B (zh) 一种风力发电机部件故障预警方法
CN102434388B (zh) 风力发电机组健康状态在线监测装置及其监测方法
Zhang et al. Fault diagnosis and prognosis using wavelet packet decomposition, Fourier transform and artificial neural network
CN111522858B (zh) 一种基于多维状态向量的抽蓄机组性能劣化预警方法
CN104792529A (zh) 基于状态空间模型的滚动轴承寿命预测方法
CN102352824B (zh) 一种基于电气信息的风力发电机健康状态监测系统和方法
CN103758696A (zh) 基于scada温度参数的风电机组安全状况评价方法
CN106650122B (zh) 一种设备变工况运行风险评估方法
CN103925155A (zh) 一种风电机组输出功率异常的自适应检测方法
CN105653851B (zh) 基于分阶段物理模型及粒子滤波的滚动轴承剩余寿命预测方法
CN103439091A (zh) 水轮机转轮叶片裂纹故障早期预警和诊断方法及系统
CN110080921B (zh) 一种抽水蓄能电站主进水阀在线监测评估方法和系统
CN103940611A (zh) 一种风电机组变工况下滚动轴承自适应异常检测方法
CN107701468A (zh) 一种混流泵在线综合监测方法及装置
CN107588937B (zh) 一种水力发电机组设备健康分析方法
CN105354768A (zh) 输变电设备状态的评估方法和系统
CN109595130A (zh) 一种风机叶片结冰故障预测方法及系统
CN106546918A (zh) 一种水电机组的故障诊断方法
CN105334472A (zh) 矿用本安电源在线剩余寿命预测方法
CN112983843B (zh) 一种潜水电泵的智能化控制系统及其控制方法
CN116050613B (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
GR01 Patent grant
GR01 Patent grant