CN112100577A - 基于长程相关性的设备运行稳定性在线监测方法及系统 - Google Patents
基于长程相关性的设备运行稳定性在线监测方法及系统 Download PDFInfo
- Publication number
- CN112100577A CN112100577A CN202010976074.4A CN202010976074A CN112100577A CN 112100577 A CN112100577 A CN 112100577A CN 202010976074 A CN202010976074 A CN 202010976074A CN 112100577 A CN112100577 A CN 112100577A
- Authority
- CN
- China
- Prior art keywords
- long
- range correlation
- data
- sequence
- operation stability
- 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
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 34
- 241001123248 Arma Species 0.000 claims abstract description 28
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 230000008859 change Effects 0.000 claims description 18
- 230000002159 abnormal effect Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 5
- 230000007774 longterm Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/005—Testing of complete machines, e.g. washing-machines or mobile phones
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明公开了基于长程相关性的设备运行稳定性在线监测方法及系统,其技术方案为:包括:获取监测数据,并对数据进行预处理;建立LRD时间序列模型,得到具有长程相关性的序列;将具有长程相关性的序列转换为具有短程相关性的序列,并确定ARMA模型的参数;利用已确定参数的ARMA模型预测数据;对预测数据进行异常决策。本发明将长程相关性引入机械设备监测领域,利用信号的长程相关性来检测机器的运行状态,解决了强环境噪声带来的机械设备在线监测设备运行稳定性的困难。
Description
技术领域
本发明涉及机械运行稳定性监测技术领域,尤其涉及一种基于长程相关性的设备运行稳定性在线监测方法及系统。
背景技术
现代机械设备向着高速化、大型化、自动化和智能化的方向发展,功能完善的同时也使得设备更加复杂。同时,现实中的工业系统几乎总是在随时间变化的不稳定条件下运行。因此,现代工业正越来越多地设法自动监测机器的动态运行状态,并在早期阶段识别出其状态变化点(例如异常、故障点)。机械设备运行稳定性在线监测可以在机械运行状态变化时发出警报从而及时采取措施,规避稳定性的进一步劣化,这不仅可以提高安全性,而且还可以提高设备和产品的可靠性、安全性、质量以及生产率等。
面向旋转机械稳定性监测存在的主要问题是由于设备的载荷波动、设备润滑程度、现场环境干扰等不可避免因素的存在导致机械设备的运行处于非平稳状态,这种时变工况对旋转机械运行稳定性的监测带来很大的困难,建立一个能够描述机器运行状态的动态模型是实现机械设备运行稳定性监测的关键步骤。然而,考虑到实际的工程场景,通过传感器收集到的状态信号通常是复杂和非平稳的。因此,设备运行中早期的变化信息很容易被强环境噪声所干扰和淹没,导致机器动态状态监测问题变得困难和复杂。
发明内容
针对现有技术存在的不足,本发明的目的是提供一种基于长程相关性的设备运行稳定性在线监测方法及系统,将长程相关性引入机械设备监测领域,利用信号的长程相关性来检测机器的运行状态,解决了强环境噪声带来的机械设备在线监测设备运行稳定性的困难。
为了实现上述目的,本发明是通过如下的技术方案来实现:
第一方面,本发明的实施例提供了一种基于长程相关性的设备运行稳定性在线监测方法,包括:
获取监测数据,并对数据进行预处理;
建立LRD时间序列模型,得到具有长程相关性的序列;
将具有长程相关性的序列转换为具有短程相关性的序列,并确定ARMA模型的参数;
利用已确定参数的ARMA模型预测数据;
对预测数据进行异常决策。
作为进一步的实现方式,计算长程相关的H参数,并计算用作长期依赖强度指标的参数d,其中,d表示微分阶。
作为进一步的实现方式,得到参数d后,利用分数阶差分算子对具有长程相关性的序列{Yt}序进行滤波,以得到具有短程相关性的序列{Zt}。
作为进一步的实现方式,用最小相似信息准则确定ARMA模型的参数。
作为进一步的实现方式,使用3σ原则进行假设检验;采用高斯分布,以检测在当前时间点k处是否发生变化。
第二方面,本发明实施例还提供了一种基于长程相关性的设备运行稳定性在线监测系统,包括:
数据预处理模块,用于获取监测数据,并对数据进行预处理;
LRD时间序列模型建立模块,用于建立LRD时间序列模型,得到具有长程相关性的序列;
ARMA模型参数确定模块,用于将具有长程相关性的序列转换为具有短程相关性的序列,并确定ARMA模型的参数;
数据预测模块,用于利用已确定参数的ARMA模型预测数据;
异常决策模块,用于对预测数据进行异常决策。
第三方面,本发明实施例还提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现所述的基于长程相关性的设备运行稳定性在线监测方法。
第四方面,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现所述的基于长程相关性的设备运行稳定性在线监测方法。
上述本发明的实施例的有益效果如下:
本发明的一个或多个实施方式将长程相关性引入机械设备监测领域,利用信号中的LRD现象来监测机器的运行状态,然后利用分数阶自回归积分滑动平均(FARIMA)模型定义LRD,有效解决了时变工况所带来的稳定性在线监测存在的困难,在实际工程应用中具有很好的潜力。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1是本发明根据一个或多个实施方式的长程相关性(LRD)预测过程示意图;
图2是本发明根据一个或多个实施方式使用FARIMA模型预测的数据与实际数据比较示例图;
图3是本发明根据一个或多个实施方式从50rpm到150rpm的速度变化检测示例。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
实施例一:
本发明的实施例提供了一种基于长程相关性的设备运行稳定性在线监测方法,包括:
获取监测数据,并对数据进行预处理;
建立LRD时间序列模型,得到具有长程相关性的序列;
将具有长程相关性的序列转换为具有短程相关性的序列,并确定ARMA模型的参数;
利用已确定参数的ARMA模型预测数据;
对预测数据进行异常决策。
具体的,如图1所示,包括以下步骤:
步骤1:通过安装在机械设备上的传感器如声传感器、振动传感器等,进行数据采集,采集的数据隐含了机器的运行状态信息,同时已证实所采集的信号,如声信号,具有长程相关性的特性。对采集到的数据流进行预处理,将给定的时间序列X(i),(i=1,2,…,N)用长度为m的非重叠窗口划分为若干个子序列,然后每个子序列的平均值计算如下:
子序列X(m)(k)的样本方差可由下式得到:
步骤2:计算长程相关的H参数。对于LRD,例如fGN(分数高斯模型)和FARIMA过程,方差还可以用(3)式表示:
步骤3:FARIMA(p,d,q)过程广泛应用于LRD时间序列的建模,其中p是自回归阶,d是微分阶,q是移动平均阶。p和q都是非负的中间值,而d是非整数值(d∈(-0.5,0.5)).一个FARIMA过程Xt:t=…,-1,0,1,…定义为:
Φ(B)(1-B)dXt=Θ(B)εt (4)
其中{∈t}为白噪声序列,B为后移运算符,定义为BXt=Xt-1,
Φ(B)=1-φ1B-φ2B2-…-φpBq (5)
Φ(B)=1+θ1B+θ2B2+…+θpBq (6)
其中:
Γ表示Gamma函数,由定义。显然,如果d=0,则FARIMA(p,d,q)模型变为本实施例所用的ARMA(p,q)模型。FARIMA(0,d,0)过程是FARIMA模型的最简单、最基本的形式。FARIMA(0,d,0)过程的性质可以描述长期依赖性,其中参数d是长期依赖强度的指标,就像FGN过程中的H参数一样,H=d+0.5。
步骤4:在得到参数d后,利用分数阶差分算子对{Yt}序列进行滤波,表达式如下:
其中:
根据公式的递归关系,可以获得以下公式:
对于d<0.5,上式为均方收敛的,因此算子是完整的。通过这种方式,具有长程相关性的序列{Yt}将转换为符合ARMA模型的具有短程相关性的序列{Zt}。
步骤5:用最小相似信息准则(AIC)确定ARMA模型的参数,Akaike信息准则(AIC)可以用来确定ARMA模型的最优参数(p,q)。在用时间序列数据拟合ARMA模型时,高斯似然可以看作是实似然函数,AIC可以写为:
以50rpm到150rpm的速度变化检测为例,原始数据、预测数据与残差变化如图3所示,变化点已在图中标示,其从上到下分别给出了原始信号数据、预测值和基于残差的检测结果。从原始信号数据可以看出,由于FARIMA模型的长相关性现象,预测值能够准确反映原始信号的变化趋势。同时,残差具有较高的预测精度,并且在变化时刻出现了显著的增长,从而成功的检测到了变化。
步骤8:在残差分析的基础上,使用3σ原则进行假设检验,该控方法采用高斯分布,以检测在当前时间点k处是否发生变化:
为了减少随机噪声的干扰,采用一种自适应检测方法,将待检测的r+1连续残差值检测为异常点,即实际变化点.决策过程如下:
1)当第一次检测到变化点时,设置r=0和变化点k为C0。
3)让r从1增加到r直至满足C0=Cr,并输出Cr作为实际变化点。随着r的增大,检测到的Cr变得更加稳定,这意味着随机噪声引起的状态变化得到了很好的抑制。
长程相关技术是从数据的角度对机械设备状态进行监测,并将机械设备数据中因机械转动与设备退化过程中监测数据得长相关性进行挖掘,建立相应得数据模型。相较于本实施例引用得长程相关性以及对应分析方法,研究较多得LSTM(长短期记忆网络)需要使用量得历史数据对模型进行训练与优化,同时对内部参数设置合理得收敛域,这需要花费大量得计算时间与内存,以致于在实现机械设备动态监测过程中会有较大得延时。而本实施例在保证精度得要求下,具有较快的计算速度和较小的计算负担。
本实施例将长程相关性引入机械设备监测领域,利用信号中的LRD现象来监测机器的运行状态,然后利用分数阶自回归积分滑动平均(FARIMA)模型定义LRD,有效解决了时变工况所带来的稳定性在线监测存在的困难,在实际工程应用中具有很好的潜力。
实施例二:
本实施例提供了一种基于长程相关性的设备运行稳定性在线监测系统,包括:
数据预处理模块,用于获取监测数据,并对数据进行预处理;
LRD时间序列模型建立模块,用于建立LRD时间序列模型,得到具有长程相关性的序列;
ARMA模型参数确定模块,用于将具有长程相关性的序列转换为具有短程相关性的序列,并确定ARMA模型的参数;
数据预测模块,用于利用已确定参数的ARMA模型预测数据;
异常决策模块,用于对预测数据进行异常决策。
实施例三:
本实施例提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现实施例一所述的基于长程相关性的设备运行稳定性在线监测方法。
实施例四:
本实施例提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现实施例一所述的基于长程相关性的设备运行稳定性在线监测方法。
以上实施例二-四中涉及的各步骤与方法实施例一相对应,具体实施方式可参见实施例一的相关说明部分。术语“计算机可读存储介质”应该理解为包括一个或多个指令集的单个介质或多个介质;还应当被理解为包括任何介质,所述任何介质能够存储、编码或承载用于由处理器执行的指令集并使处理器执行本发明中的任一方法。
本领域技术人员应该明白,上述本发明的各模块或各步骤可以用通用的计算机装置来实现,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。本发明不限制于任何特定的硬件和软件的结合。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (10)
1.基于长程相关性的设备运行稳定性在线监测方法,其特征在于,包括:
获取监测数据,并对数据进行预处理;
建立LRD时间序列模型,得到具有长程相关性的序列;
将具有长程相关性的序列转换为具有短程相关性的序列,并确定ARMA模型的参数;
利用已确定参数的ARMA模型预测数据;
对预测数据进行异常决策。
2.根据权利要求1所述的基于长程相关性的设备运行稳定性在线监测方法,其特征在于,计算长程相关的H参数,并计算用作长期依赖强度指标的参数d,其中,d表示微分阶。
3.根据权利要求2所述的基于长程相关性的设备运行稳定性在线监测方法,其特征在于,得到参数d后,利用分数阶差分算子对具有长程相关性的序列{Yt}序进行滤波,以得到具有短程相关性的序列{Zt}。
4.根据权利要求1所述的基于长程相关性的设备运行稳定性在线监测方法,其特征在于,用最小相似信息准则确定ARMA模型的参数。
7.根据权利要求6所述的基于长程相关性的设备运行稳定性在线监测方法,其特征在于,使用3σ原则进行假设检验;采用高斯分布,以检测在当前时间点k处是否发生变化。
8.基于长程相关性的设备运行稳定性在线监测系统,其特征在于,包括:
数据预处理模块,用于获取监测数据,并对数据进行预处理;
LRD时间序列模型建立模块,用于建立LRD时间序列模型,得到具有长程相关性的序列;
ARMA模型参数确定模块,用于将具有长程相关性的序列转换为具有短程相关性的序列,并确定ARMA模型的参数;
数据预测模块,用于利用已确定参数的ARMA模型预测数据;
异常决策模块,用于对预测数据进行异常决策。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-7任一项所述的基于长程相关性的设备运行稳定性在线监测方法。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-7任一项所述的基于长程相关性的设备运行稳定性在线监测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010976074.4A CN112100577A (zh) | 2020-09-16 | 2020-09-16 | 基于长程相关性的设备运行稳定性在线监测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010976074.4A CN112100577A (zh) | 2020-09-16 | 2020-09-16 | 基于长程相关性的设备运行稳定性在线监测方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112100577A true CN112100577A (zh) | 2020-12-18 |
Family
ID=73760215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010976074.4A Pending CN112100577A (zh) | 2020-09-16 | 2020-09-16 | 基于长程相关性的设备运行稳定性在线监测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100577A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112835883A (zh) * | 2021-02-01 | 2021-05-25 | 北京爱奇艺科技有限公司 | 一种滤波方法、装置及电子设备 |
CN113588259A (zh) * | 2021-08-03 | 2021-11-02 | 山东中科普锐检测技术有限公司 | 一种设备振动信号标度曲线转折点检测方法及工况监测装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107547269A (zh) * | 2017-08-14 | 2018-01-05 | 浙江大学 | 基于farima的智能变电站通信流量阈值模型的构建方法 |
CN108710757A (zh) * | 2018-05-18 | 2018-10-26 | 山东大学 | 基于时变参数预测模型的机械运行状态监测方法及装置 |
-
2020
- 2020-09-16 CN CN202010976074.4A patent/CN112100577A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107547269A (zh) * | 2017-08-14 | 2018-01-05 | 浙江大学 | 基于farima的智能变电站通信流量阈值模型的构建方法 |
CN108710757A (zh) * | 2018-05-18 | 2018-10-26 | 山东大学 | 基于时变参数预测模型的机械运行状态监测方法及装置 |
Non-Patent Citations (3)
Title |
---|
中国地震局监测预报司: "高速铁路光传送网(OTN)保护技术", 31 July 2002, pages: 121 - 124 * |
孙勇;白光伟;赵露;: "基于小波分形自回归整合滑动平均模型的网络流量预测", 计算机应用, no. 04 * |
安德洪 等: "一类基于FARIMA过程的电梯导轨振动模型", 机械设计, pages 1 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112835883A (zh) * | 2021-02-01 | 2021-05-25 | 北京爱奇艺科技有限公司 | 一种滤波方法、装置及电子设备 |
CN113588259A (zh) * | 2021-08-03 | 2021-11-02 | 山东中科普锐检测技术有限公司 | 一种设备振动信号标度曲线转折点检测方法及工况监测装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111189639B (zh) | 一种基于瞬时频率优化vmd的轴承故障诊断方法 | |
Shardt et al. | A new soft-sensor-based process monitoring scheme incorporating infrequent KPI measurements | |
Camci et al. | Feature evaluation for effective bearing prognostics | |
JP2018156151A (ja) | 異常検知装置及び機械学習装置 | |
CN112232447A (zh) | 一种电力设备状态监测数据的完整样本集的构建方法 | |
CN112100577A (zh) | 基于长程相关性的设备运行稳定性在线监测方法及系统 | |
JP2020060558A (ja) | 工場ロボットの状態ベースメンテナンスのための運動非感知性特徴 | |
US20190196458A1 (en) | Method for selecting leading associated parameter and method for combining critical parameter and leading associated parameter for equipment prognostics and health management | |
CN108257365B (zh) | 一种基于全局不确定性证据动态融合的工业报警器设计方法 | |
Wang et al. | Remaining useful life prediction considering joint dependency of degradation rate and variation on time-varying operating conditions | |
CN105675320B (zh) | 一种基于声学信号分析的机械系统运行状态实时监控方法 | |
Shao et al. | Gearbox vibration monitoring using extended Kalman filters and hypothesis tests | |
Heydarzadeh et al. | Gearbox fault diagnosis using power spectral analysis | |
CN116086537A (zh) | 一种设备状态监测方法、装置、设备及存储介质 | |
CN105718733B (zh) | 基于模糊贴近度和粒子滤波的故障预报方法 | |
CN113269327A (zh) | 一种基于机器学习的流量异常预测方法 | |
Dhini et al. | Review on some multivariate statistical process control methods for process monitoring | |
CN116360384A (zh) | 用于诊断和监控信息物理系统中的异常的系统和方法 | |
CN114818813A (zh) | 一种基于自注意力卷积神经网络的刀具磨损值预测算法 | |
EP3819608A1 (en) | Detecting rotor anomalies by determining vibration trends during transient speed operation | |
CN110045716B (zh) | 一种闭环控制系统早期故障检测和诊断方法和系统 | |
Zhao et al. | A unified framework for fault detection and diagnosis using particle filter | |
Raj | An Improved Accuracy and Efficiency Based Defect Detection Method for Industrial Signal Processing | |
Moraes et al. | Comparing the inertial effect of MEWMA and multivariate sliding window schemes with confidence control charts | |
Nishino et al. | Classification method based on Taguchi’s T-method for small sample sizes |
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 |