CN110018670B - 一种基于动态关联规则挖掘的工业过程异常工况预测方法 - Google Patents

一种基于动态关联规则挖掘的工业过程异常工况预测方法 Download PDF

Info

Publication number
CN110018670B
CN110018670B CN201910244877.8A CN201910244877A CN110018670B CN 110018670 B CN110018670 B CN 110018670B CN 201910244877 A CN201910244877 A CN 201910244877A CN 110018670 B CN110018670 B CN 110018670B
Authority
CN
China
Prior art keywords
association rule
data
window
item set
training
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.)
Active
Application number
CN201910244877.8A
Other languages
English (en)
Other versions
CN110018670A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201910244877.8A priority Critical patent/CN110018670B/zh
Publication of CN110018670A publication Critical patent/CN110018670A/zh
Application granted granted Critical
Publication of CN110018670B publication Critical patent/CN110018670B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/406—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by monitoring or safety
    • G05B19/4065—Monitoring tool breakage, life or condition
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/37—Measurements
    • G05B2219/37616—Use same monitoring tools to monitor tool and workpiece

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明公开了一种基于动态关联规则挖掘的工业过程异常工况预测方法，可应用于工业过程的故障预测与健康管理。本发明利用滑动窗口的方式挖掘工业设备运行参数动态关联规则，并将其引入工业过程异常工况预测中。本发明考虑关联规则的时间特性，利用滑动窗口限制数据长度，提出一种用于两两运行参数动态关联规则挖掘的关联规则挖掘算法，随后将关联规则挖掘结果引入小波神经网络预测中，以动态关联规则不断更新网络，以获得更为准确的预测结果。对于工程上的故障预测和健康管理有重大应用价值。

Description

一种基于动态关联规则挖掘的工业过程异常工况预测方法
技术领域
本发明属于可靠性维护工程技术领域，涉及一种基于动态关联规则挖掘的工业过程异常工况预测方法。
背景技术
随着复杂***的不断出现以及工业过程实时监测的需求不断增加，现代工业设备在运行过程中往往配备多个传感器对其运行状态进行监测。同时，设备运行过程中可能会出现多种故障模式，某一故障可能对应若干征兆，在此情况下，单传感器信息已无法完全体现设备运行状态，基于多传感器信息的故障预测应运而生。基于多传感器信息的故障预测旨在利用全面的传感器信息分析设备的运行状态，从而进行更可靠的设备诊断和预测。随着传感技术的持续发展，使用多个传感器进行设备的状态监测、故障诊断和预测已经成为发展趋势。
设备运行过程中其运行参数之间存在一定关联性，在故障预测领域，目前仍鲜有将关联规则挖掘与故障预测相结合的工作。传统的关联规则挖掘算法认为，对于某一对象***所发现的关联规则是不变的即在分析该***其他数据集时也是有效的。如此挖掘出来的关联规则并没有考虑到时间因素，由此得到的是一种静态的关联规则。而实际上，随着数据的更新，可能伴随着***运行状态的变化，随之而来的是数据特性的变化，从而关联规则也可能发生变化。在工业设备运行过程中，运行状态并不是固定不变的。随着数据的更新或当设备运行状态发生变化，其运行参数的特性也会随之变化，从而运行参数的关联性也可能发生了变化。在此情况下，若一直使用从最初的数据集中挖掘得到的关联规则来指导后续决策显然是不合理的。
发明内容
针对现有技术的现状，本发明的目的是解决现有数据驱动的预测技术中鲜有考虑传感器数据动态关联的情况，提出一种基于运行参数动态关联规则的设备异常工况预测方法，利用动态关联规则更新小波神经网络并进行异常工况预测(故障预测)。
现将本发明的构思阐述如下：
本发明针对工业设备运行过程中运行参数之间的关联性可能随设备运行状态的变化而发生变化的问题，引入了动态关联规则挖掘。本发明利用滑动窗口的方式，对每一个窗口内的数据进行关联规则挖掘，从而随着窗口的滑动得到动态变化的关联规则。为得到参数关联性的完整变化情况，本发明所提方法在参数层面的关联规则挖掘时并没有设置最小支持度和置信度阈值，而是保存所有两两参数的关联规则，并以支持度和置信度描述其强度。随后，将动态关联规则挖掘结果引入小波神经网络的训练和更新中，通过由动态关联规则挖掘获得的关联参数进行异常工况预测。
根据以上发明构思，本发明提出了一种基于动态关联规则挖掘的工业过程异常工况预测方法，具体步骤如下：
步骤1：基于滑动窗口对工业过程中传感器所采集的测量值时间序列进行数据预处理，对时间序列进行分段线性化表示、聚类、符号化，生成适用于关联规则挖掘的事务集；
步骤2：采用两阶段的频繁项集产生策略生成频繁项集，并挖掘两两参数的关联规则；
步骤3：利用初始数据集进行关联规则挖掘，基于初始关联规则挖掘结果，训练初始小波神经网络模型；
步骤4：基于数据更新进行数据窗口的滑动，对于新数据窗口进行关联规则更新和小波神经网络模型更新；
步骤5：基于更新的小波网络模型进行异常工况预测，并在预测到异常工况发生前不断更新关联规则和小波神经网络模型。
基于上述方案，各步骤可具体采用如下实现方式：
作为优选，所述的步骤1具体通过如下子步骤实现：
步骤1.1：记传感器测量序列为
Figure BDA0002010777970000021
N为传感器数量，K为序列长度,k＝1,2,…,K；记滑动窗口长度为L，每次滑动距离设为l，将滑动窗口记为Wk，窗口中包含的数据为
Figure BDA0002010777970000022
步骤1.2：对于每一数据窗口进行时间序列分段线性化表示、聚类、符号化处理，具体过程如下：
对每一个数据窗口Wk执行1.2.1～1.2.5过程：
1.2.1.设置初始拟合起点为
Figure BDA0002010777970000031
终点为
Figure BDA0002010777970000032
h＝2；拟合误差阈值为ωE；初始化分段点计数值count＝1；
1.2.2.依次对每个拟合起点
Figure BDA0002010777970000033
执行步骤1.1)-步骤1.4)：
1.1)首先计算end＝start+h；
1.2)对于数据
Figure BDA0002010777970000034
利用最小二乘法进行拟合，计算拟合误差ERR；
1.3)若拟合误差ERR不大于拟合误差阈值ωE，则h＝h+1，重新跳转至步骤1)；
1.4)若拟合误差ERR大于拟合误差阈值ωE，则保存分段点
Figure BDA0002010777970000035
Figure BDA0002010777970000036
至Pi；保存拟合线段
Figure BDA0002010777970000037
采用三元组的方式表示该拟合线段
Figure BDA0002010777970000038
并保存至
Figure BDA0002010777970000039
ki表示该线段的斜率，
Figure BDA00020107779700000310
表示该线段在时间轴上的跨度，ri表示该线段数据的增长率；重置拟合起点start＝start+h，重置h＝ 2；更新count＝count+1；
1.2.3.循环执行1.2.2直到end大于k+L-1结束，得到拟合后的线段化时间序列
Figure BDA00020107779700000311
和分段点
Figure BDA00020107779700000312
组成的分段点序列Pi；
1.2.4.以Wk-1窗口的聚类中心为初始聚类中心，采用K-means聚类算法对三元组序列
Figure BDA00020107779700000313
进行聚类，并为不同类的线段分配不同符号得到符号化序列
Figure BDA00020107779700000314
其中聚类过程中采用基于欧氏距离的指标描述两条线段si和sj之间的相似度：
Figure BDA00020107779700000315
其中，dij即表示线段si和sj的相似度，dij越小，则表示两条线段有更相似的变化形态，ωk和ωr表示权重；
1.2.5.对于两个运行参数Vi和Vj，分别得到符号化后的序列
Figure BDA00020107779700000316
和
Figure BDA00020107779700000317
合并两个运行参数的分段点序列Pi和Pj，并按合并后的分段点对其符号化序列
Figure BDA0002010777970000041
和
Figure BDA0002010777970000042
进行分割重构，得到重构后的符号化序列
Figure BDA0002010777970000043
和
Figure BDA0002010777970000044
由其构成事务集，nij-1为Pi和Pj合并后的分段点个数。
作为优选，所述的步骤2具体通过如下子步骤实现：
步骤2.1：由步骤1获得的事务集的每一事务记为
Figure BDA0002010777970000045
另将
Figure BDA0002010777970000046
和
Figure BDA0002010777970000047
中所包含的线段类别符号分别记为
Figure BDA0002010777970000048
和
Figure BDA0002010777970000049
mi和mj分别为
Figure BDA00020107779700000410
和
Figure BDA00020107779700000411
中所包含的线段类别数；设频繁项集产生过程的最小支持度阈值为minsup1；
步骤2.2：对于初始数据窗口W1进行关联规则挖掘，其具体过程如下：
2.2.1)对所有传感器的传感器测量序列数据集进行扫描以计算每个项的支持度，记σ(·)表示某项或项集的支持度计数，初始为0；假设
Figure BDA00020107779700000412
的类别符号为tk， e表示i或j，t表示c或d；
对每一个
Figure BDA00020107779700000413
计算σ(tk)＝σ(tk)+1；
对每一个tk，判断
Figure BDA00020107779700000414
是否成立，若成立则认为tk为频繁1- 项集，保留tk并记录相应的支持度计数σ(tk)；
2.2.2)使用已得到的频繁1-项集tk构成2-项集，并发现频繁2-项集，分别记cp和dq分别为经过步骤2.2.1)从原线段类别符号
Figure BDA00020107779700000415
和
Figure BDA00020107779700000416
中保留的项；
对于每一个{cp,dq}执行操作：对每一个存在于
Figure BDA00020107779700000417
中的{cp,dq}，计算σ({cp,dq})＝σ({cp,dq})+1；如果
Figure BDA00020107779700000418
不小于minsup1，则认为{cp,dq}为频繁2-项集，保留{cp,dq}并记录相应的支持度计数；
2.2.3)使用步骤2.2.2中获得的频繁2-项集{cp,dq}计算每两个运行参数Vi和 Vj在整个数据集中的支持度，计算过程如下：对每两个运行参数Vi和Vj构成的项集{Vi,Vj}，计算σ({Vi,Vj})＝sum(σ({cp,dq}))，两个参数的关联规则支持度
Figure BDA0002010777970000051
并计算σ(Vi)＝sum(σ(cp))；σ(Vj)＝ sum(σ(dq))；
2.2.4)对于每组{Vi,Vj}，产生如下关联规则：Vj→Vi和Vi→Vj；对于每一个关联规则，计算其置信度为
Figure BDA0002010777970000052
步骤2.3：对于每一个新的数据窗口，采用以下过程进行频繁项集产生并生成关联规则：
对每一个k>1的数据窗口Wk执行步骤2.3.1)～2.3.3)：
2.3.1)对每一滑出窗口的事务
Figure BDA0002010777970000053
记为{ca1,da2}，判断{ca1,da2}是否为频繁项集，若是则更新σ({ca1,da2})＝σ({ca1,da2})-1；
2.3.2)对每一滑入窗口的事务
Figure BDA0002010777970000054
记为{cb1,db2}，判断{cb1,db2}是否为频繁项集，若是，则更新σ({cb1,db2})＝σ({cb1,db2})+1，若否，则在数据窗口Wk内计算该{cb1,db2}的支持度计数为σ({cc1,dc2})；
2.3.3)对每一个{ca1,da2}，若
Figure BDA0002010777970000055
则更新{ca1,da2}相应项集的支持度计数，若不满足则剔除相应项集；对每一个{cb1,db2}，若
Figure BDA0002010777970000056
则更新{cb1,db2}相应项集的支持度计数，若不满足则剔除相应项集；对每一个{cc1,dc2}，若
Figure BDA0002010777970000057
则更新{cc1,dc2}相应项集的支持度计数，若不满足则剔除相应项集；
2.3.4)利用步骤2.2.3)和2.2.4)更新每一关联规则的支持度和置信度。
作为优选，所述的步骤3具体通过如下子步骤实现：
记预设的预测步长为lp，由训练数据集挖掘到的关联规则中提取的一组关联参数为{V1,V2,…,Vu}，u为提取出的一组关联参数个数；对于每一参数Vi，其测量值序列为
Figure BDA0002010777970000058
构造以下矩阵Itrain为初始的神经网络训练输入：
Figure BDA0002010777970000059
其中，Itrain中每一列为一个训练输入样本；构造训练输出Otrain为，
Figure BDA0002010777970000061
利用上述公式构造的训练样本，对小波神经网络进行训练；在网络初始化时，利用关联规则Vi→Vu对应的初始置信度ωi设置网络输入层和隐含层之间的初始权值，i＝1,2,…u-1。
作为优选，所述的步骤4具体通过如下子步骤实现：
步骤4.1：关联规则更新做法为：记在设备运行过程中传感器已采集到的c 个新数据为
Figure BDA0002010777970000062
c>L；利用新采集到的数据中最近的L个数据
Figure BDA0002010777970000063
通过步骤2.3更新关联规则，从而更新关联规则Vi→Vu所对应的置信度记为
Figure BDA0002010777970000064
步骤4.2：模型更新做法为：结合训练数据以及新采集到的数据构造新的模型训练样本，即数据集为
Figure BDA0002010777970000065
构造矩阵Itest为神经网络更新的输入，构造Otest为输出，
Figure BDA0002010777970000066
Figure BDA0002010777970000067
利用上述公式构造新的训练样本，对小波神经网络进行训练；在网络初始化时，利用新的置信度
Figure BDA0002010777970000068
设置网络输入层和隐含层之间的初始权值。
作为优选，所述的步骤5具体通过如下子步骤实现：
记预设的异常工况(失效)发生阈值为ωp，对于最新采集到的数据，利用步骤4中更新的模型进行lp步预测，若得到的目标参数预测值相对于初始正常的漂移量超过所设阈值，则认为异常工况发生。
作为优选，在预测过程结束前，每更新预定数量的数据，就回到步骤4进行关联规则和模型的更新。
本发明利用滑动窗口的方式挖掘工业设备运行参数动态关联规则，并将其引入工业过程异常工况预测中。本发明考虑关联规则的时间特性，利用滑动窗口限制数据长度，提出一种用于两两运行参数动态关联规则挖掘的关联规则挖掘算法，随后将关联规则挖掘结果引入小波神经网络预测中，以动态关联规则不断更新网络，以获得更为准确的预测结果。对于工程上的故障预测和健康管理有重大应用价值。
附图说明
图1为实施例中IDV(13)关联规则(v13,v16,v36→v7)置信度变化曲线；
图2为实施例中IDV(13)关联规则(v35,v36→v11)置信度变化曲线；
图3为实施例中IDV(13)变量7预测结果；
图4为实施例中IDV(13)变量7预测误差；
图5为实施例中IDV(13)变量11预测结果；
图6为实施例中IDV(13)变量11预测误差。
具体实施方式
现结合附图对本发明的具体实施方式作进一步的说明。
下面本实施例通过田纳西-伊斯曼(TE)过程仿真数据来具体阐述具体操作步骤以及验证方法的效果。
该数据集的采样间隔为3分钟，数据集记录了该采样间隔下各传感器采集到的变量测量值。在每一种运行条件下(正常运行状态以及21种预设故障下的故障运行状态)，仿真过程的测量数据都将生成两类数据集，即训练集和测试集。其中，对于训练集的采集过程，是在仿真过程运行25个小时的情况下获得的所有52个变量的测量值，其中，除正常运行状态采集到的训练集外，其余21个训练集数据的采集是在仿真过程运行1小时候引入故障，并只记录随后24个小时的测量数据。也即是说，正常运行状态的训练集有500个观测样本，其余21个故障状态下采集的训练集均为480个观测样本。另外，对于22个测试集，其数据是仿真过程运行48个小时所采集到的所有变量测量值，也就是说，每一个测试集中包含960个样本数据。需要注意的是，在对21种过程故障进行仿真时，相应的故障是在仿真运行8小时后引入的。因此，对于21个故障运行状态下的测试集，其前160个观测样本为正常数据，后800个观察样本为故障数据。在TE过程仿真模型中，只有IDV(13)是一个缓变故障，因此，在本例中，我们使用IDV(13)的相关数据进行实验。
工业过程异常工况预测方法具体执行过程如下：
步骤1：基于滑动窗口对工业过程中传感器所采集的测量值时间序列进行数据预处理，对时间序列进行分段线性化表示、聚类、符号化，生成适用于关联规则挖掘的事务集。本步骤具体通过如下子步骤实现：
步骤1.1：记传感器测量序列为
Figure BDA0002010777970000081
N为传感器数量，K为序列长度,k＝1,2,…,K；记滑动窗口长度为L，每次滑动距离设为l，将滑动窗口记为Wk，窗口中包含的数据为
Figure BDA0002010777970000082
需要注意的是，在本发明中，i、j作为上标是表示传感器的编号，作为下标是仅表示序数，与传感器编号无关。
步骤1.2：对于每一数据窗口进行时间序列分段线性化表示、聚类、符号化处理，具体过程如下：
对每一个数据窗口Wk执行1.2.1～1.2.5过程：
1.2.1.设置初始拟合起点为
Figure BDA0002010777970000083
终点为
Figure BDA0002010777970000084
h＝2；拟合误差阈值为ωE；初始化分段点计数值count＝1；
1.2.2.依次对每个拟合起点
Figure BDA0002010777970000085
执行步骤1.1)-步骤1.4)：
1.1)首先计算end＝start+h；
1.2)对于数据
Figure BDA0002010777970000086
利用最小二乘法进行拟合，计算拟合误差ERR；
1.3)若拟合误差ERR不大于拟合误差阈值ωE，则h＝h+1，重新跳转至步骤1)；
1.4)若拟合误差ERR大于拟合误差阈值ωE，则保存分段点
Figure BDA0002010777970000087
Figure BDA0002010777970000088
至Pi；保存拟合线段
Figure BDA0002010777970000089
采用三元组的方式表示该拟合线段
Figure BDA00020107779700000810
并保存至
Figure BDA00020107779700000811
ki表示该线段的斜率，
Figure BDA00020107779700000812
表示该线段在时间轴上的跨度，ri表示该线段数据的增长率；重置拟合起点start＝start+h，重置h＝ 2(重置拟合起点和h)；更新count＝count+1；
1.2.3.循环执行1.2.2直到end大于k+L-1结束，得到最小二乘法拟合后的线段化时间序列
Figure BDA0002010777970000091
(其中具有多条拟合线段)和分段点
Figure BDA0002010777970000092
组成的分段点序列Pi；
1.2.4.以Wk-1窗口的聚类中心为初始聚类中心，采用K-means聚类算法对三元组序列
Figure BDA0002010777970000093
进行聚类，并为不同类的线段分配不同符号得到符号化序列
Figure BDA0002010777970000094
其中聚类过程中采用基于欧氏距离的指标描述两条线段si和sj之间的相似度：
Figure BDA0002010777970000095
其中，dij即表示线段si和sj的相似度，dij越小，则表示两条线段有更相似的变化形态，ωk和ωr表示权重；
1.2.5.对于两个运行参数Vi和Vj，分别得到符号化后的序列
Figure BDA0002010777970000096
和
Figure BDA0002010777970000097
合并两个运行参数的分段点序列Pi和Pj，并按合并后的分段点对其符号化序列
Figure BDA0002010777970000098
和
Figure BDA0002010777970000099
进行分割重构，得到重构后的符号化序列
Figure BDA00020107779700000910
和
Figure BDA00020107779700000911
由其构成事务集，nij-1为Pi和Pj合并后的分段点个数。
步骤2：采用两阶段的频繁项集产生策略生成频繁项集，并挖掘两两参数的关联规则。本步骤具体通过如下子步骤实现：
步骤2.1：由步骤1获得的事务集的每一事务记为
Figure BDA00020107779700000912
另将
Figure BDA00020107779700000913
和
Figure BDA00020107779700000914
中所包含的线段类别符号分别记为
Figure BDA00020107779700000915
和
Figure BDA00020107779700000916
mi和mj分别为
Figure BDA00020107779700000917
和
Figure BDA00020107779700000918
中所包含的线段类别数；设频繁项集产生过程的最小支持度阈值为minsup1；在本例中，设置最小支持度阈值为0.2；
步骤2.2：对于初始数据窗口W1进行关联规则挖掘，其具体过程如下：
2.2.1)对所有传感器的传感器测量序列数据集进行扫描以计算每个项的支持度，记σ(·)表示某项或项集的支持度计数，初始为0；假设
Figure BDA00020107779700000919
的类别符号为tk， e表示i或j，t表示c或d(即ck或dk)；
对每一个
Figure BDA0002010777970000101
更新计算σ(tk)＝σ(tk)+1；
对每一个tk，判断
Figure BDA0002010777970000102
是否成立，若成立则认为tk为频繁1- 项集，保留tk并记录相应的支持度计数σ(tk)；
2.2.2)使用已得到的频繁1-项集tk构成2-项集，并发现频繁2-项集，分别记cp和dq分别为经过步骤2.2.1)从原线段类别符号
Figure BDA0002010777970000103
和
Figure BDA0002010777970000104
中保留的项；
对于每一个{cp,dq}执行操作：对每一个存在于
Figure BDA0002010777970000105
中的{cp,dq}，计算σ({cp,dq})＝σ({cp,dq})+1；如果
Figure BDA0002010777970000106
不小于minsup1，则认为{cp,dq}为频繁2-项集，保留{cp,dq}并记录相应的支持度计数；
2.2.3)使用步骤2.2.2中获得的频繁2-项集{cp,dq}计算每两个运行参数Vi和 Vj在整个数据集中的支持度，计算过程如下：对每两个运行参数Vi和Vj构成的项集{Vi,Vj}，计算σ({Vi,Vj})＝sum(σ({cp,dq}))，两个参数的关联规则支持度
Figure BDA0002010777970000107
并计算σ(Vi)＝sum(σ(cp))；σ(Vj)＝ sum(σ(dq))；
2.2.4)对于每组{Vi,Vj}，产生如下关联规则：Vj→Vi和Vi→Vj；对于每一个关联规则，计算其置信度为
Figure BDA0002010777970000108
在本例中，对于IDV(13)，通过训练集挖掘初始的关联规则。设置支持度阈值为0.2，置信度阈值0.7，得到几组不同的关联规则，其结果可见表1，其中包含两组关联参数，即(变量13,变量16,变量36→变量7)和(变量35,变量36→变量11)。
步骤2.3：对于每一个新的数据窗口，采用以下过程进行频繁项集产生并生成关联规则。在本例中，设置动态关联规则挖掘的滑动窗口长度为100，且每次滑动1个数据点。
对每一个k>1的数据窗口Wk执行步骤2.3.1)～2.3.3)：
2.3.1)对每一滑出窗口的事务
Figure BDA0002010777970000109
记为{ca1,da2}，判断{ca1,da2}是否为频繁项集，若是则更新σ({ca1,da2})＝σ({ca1,da2})-1；
2.3.2)对每一滑入窗口的事务
Figure BDA0002010777970000111
记为{cb1,db2}，判断{cb1,db2}是否为频繁项集，若是，则更新σ({cb1,db2})＝σ({cb1,db2})+1，若否，则在数据窗口Wk内计算该{cb1,db2}的支持度计数为σ({cc1,dc2})；
2.3.3)对每一个{ca1,da2}，若
Figure BDA0002010777970000112
则更新{ca1,da2}相应项集的支持度计数，若不满足则剔除相应项集；对每一个{cb1,db2}，若
Figure BDA0002010777970000113
则更新{cb1,db2}相应项集的支持度计数，若不满足则剔除相应项集；对每一个{cc1,dc2}，若
Figure BDA0002010777970000114
则更新{cc1,dc2}相应项集的支持度计数，若不满足则剔除相应项集；
2.3.4)利用步骤2.2.3)和2.2.4)更新每一关联规则的支持度和置信度。
利用IDV(13)的测试集，从第100个点开始，每更新1个采样数据，本例就用包括这个点的前100个点来挖掘关联规则，得到的置信度变化情况如图1和图2所示。由图可知，关联规则(置信度)随着数据的更新在不断变化，但对于IDV(13) 的这两组参数，其数值均保持在一个较高的水平。
步骤3：利用初始数据集进行关联规则挖掘，基于初始关联规则挖掘结果，训练初始小波神经网络模型。本步骤具体通过如下子步骤实现：
记预设的预测步长为lp，在本例中为10。由训练数据集挖掘到的关联规则中提取的一组关联参数为{V1,V2,…,Vu}，u为提取出的一组关联参数个数；对于每一参数Vi，其测量值序列为
Figure BDA0002010777970000115
构造以下矩阵 Itrain为初始的神经网络训练输入：
Figure BDA0002010777970000116
其中，Itrain中每一列为一个训练输入样本；构造训练输出Otrain为，
Figure BDA0002010777970000117
利用上述公式构造的训练样本，对小波神经网络进行训练；在网络初始化时，利用关联规则Vi→Vu对应的初始置信度ωi设置网络输入层和隐含层之间的初始权值，i＝1,2,…u-1。
步骤4：基于数据更新进行数据窗口的滑动，对于新数据窗口进行关联规则更新和小波神经网络模型更新。本步骤具体通过如下子步骤实现：
步骤4.1：关联规则更新做法为：记在设备运行过程中传感器已采集到的c 个新数据为
Figure BDA0002010777970000121
c>L；利用新采集到的数据中最近的L个数据
Figure BDA0002010777970000122
通过步骤2.3更新关联规则，从而更新关联规则Vi→Vu所对应的置信度记为
Figure BDA0002010777970000123
需要注意的是，其中c需大于预设的滑动窗口长度L。也就是说，仅当设备开始运行且采集到的数据长度大于L时，才开始进行关联规则的更新。
步骤4.2：模型更新做法为：结合训练数据以及新采集到的数据构造新的模型训练样本，即数据集为
Figure BDA0002010777970000124
构造矩阵Itest为神经网络更新的输入，构造Otest为输出，
Figure BDA0002010777970000125
Figure BDA0002010777970000126
利用上述公式构造新的训练样本，对小波神经网络进行训练；在网络初始化时，利用新的置信度
Figure BDA0002010777970000127
设置网络输入层和隐含层之间的初始权值。在本例中，利用测试集的前300个数据，每更新10个数据，便对关联规则和模型进行更新。
步骤5：基于更新的小波网络模型进行异常工况预测，并在预测到异常工况发生前不断更新关联规则和小波神经网络模型。本步骤具体通过如下子步骤实现：
记预设的异常工况(失效)发生阈值为ωp，本例中异常工况(失效)发生的阈值为ωp＝0.015。对于最新采集到的数据，利用步骤4中更新的模型进行lp步预测，若得到的目标参数预测值相对于初始正常的漂移量超过所设阈值，则认为异常工况发生。在预测过程结束前，每更新U个数据则回到步骤4进行关联规则和模型的更新，U的数值可以调整，本例中U为10。
表1关联规则
规则前项 规则后项 置信度
变量13 变量7 0.7527
变量16 变量7 0.7446
变量36 变量7 0.7017
变量35 变量11 0.7513
变量36 变量11 0.7390
表2总预测误差率
引入动态关联规则 引入关联规则 未引入关联规则
变量7 0.7602 1.0247 1.7423
变量11 0.8075 0.8309 1.1884
本例为验证本章所用方法的有效性，将引入动态关联规则的预测结果与未引入关联规则以及引入静态关联规则的情况进行对比，预测结果和预测误差的对比见图3-图6。在图3和图5中，竖实线为在预设的失效阈值下的实际失效时刻，竖点划线为引入动态关联规则的失效时刻预测值，竖虚线为引入静态关联规则的失效时刻预测值，竖点线为未引入关联规则的失效时刻预测值。另外，为了更好地量化预测误差，本例计算两个参数预测在不同条件下的总误差率，如表2所示。
从预测结果(图3和图5)来看，与另外两种情况相比，引入动态关联规则可以获得更准确的异常工况预测结果，尤其是与未引入关联规则的情况相比，从预测曲线中便可看出引入动态关联规则具有明显优势。需要注意的是，在图5 中，引入静态关联规则以及未引入关联规则的情况下的预测值与真实值出现了较大的偏差，这其实是由于该参数的失效实际值正好处于某一极值点，故若在此处存在细微的预测误差便会导致异常工况预测结果出现较大偏差。实际上，可以设置多个不同的越限阈值对参数进行监测，以避免这样的情况带来的风险。另外，从预测误差的角度(图4和图6)来看，与引入静态关联规则的情况相比，引入动态关联规则可以在一定程度上减小预测误差，从表2可以看出，对于变量7，引入动态关联规则的提升效果较为明显，而对于变量11，由于关联规则变化并不大，所以其效果并不显著；与未引入关联规则的情况相比，从图4和图6中可以看出引入动态关联规则可以明显降低预测误差，由表2也可以更直观得看出引入动态关联规则相较于未引入关联规则具有明显优势。

Claims (6)

1.一种基于动态关联规则挖掘的工业过程异常工况预测方法，其特征在于，具体步骤如下：
步骤1：基于滑动窗口对工业过程中传感器所采集的测量值时间序列进行数据预处理，对时间序列进行分段线性化表示、聚类、符号化，生成适用于关联规则挖掘的事务集；
步骤2：采用两阶段的频繁项集产生策略生成频繁项集，并挖掘两两参数的关联规则；
步骤3：利用初始数据集进行关联规则挖掘，基于初始关联规则挖掘结果，训练初始小波神经网络模型；
步骤4：基于数据更新进行数据窗口的滑动，对于新数据窗口进行关联规则更新和小波神经网络模型更新；
步骤5：基于更新的小波网络模型进行异常工况预测，并在预测到异常工况发生前不断更新关联规则和小波神经网络模型；
所述的步骤1具体通过如下子步骤实现：
步骤1.1：记传感器测量序列为
Figure FDA0002458389690000011
N为传感器数量，K为序列长度,k＝1,2,...,K；记滑动窗口长度为L，每次滑动距离设为l，将滑动窗口记为Wk，窗口中包含的数据为
Figure FDA0002458389690000012
步骤1.2：对于每一数据窗口进行时间序列分段线性化表示、聚类、符号化处理，具体过程如下：
对每一个数据窗口Wk执行1.2.1～1.2.5过程：
1.2.1.设置初始拟合起点为
Figure FDA0002458389690000013
终点为
Figure FDA0002458389690000014
拟合误差阈值为ωE；初始化分段点计数值count＝1；
1.2.2.依次对每个拟合起点
Figure FDA0002458389690000015
执行步骤1.1)-步骤1.4)：
1.1)首先计算end＝start+h；
1.2)对于数据
Figure FDA0002458389690000016
利用最小二乘法进行拟合，计算拟合误差ERR；
1.3)若拟合误差ERR不大于拟合误差阈值ωE，则h＝h+1，重新跳转至步骤1)；
1.4)若拟合误差ERR大于拟合误差阈值ωE，则保存分段点
Figure FDA0002458389690000021
Figure FDA0002458389690000022
至Pi；保存拟合线段
Figure FDA0002458389690000023
采用三元组的方式表示该拟合线段
Figure FDA0002458389690000024
并保存至
Figure FDA0002458389690000025
ki表示该线段的斜率，
Figure FDA0002458389690000026
表示该线段在时间轴上的跨度，ri表示该线段数据的增长率；重置拟合起点start＝start+h，重置h＝2；更新count＝count+1；
1.2.3.循环执行1.2.2直到end大于k+L-1结束，得到拟合后的线段化时间序列
Figure FDA0002458389690000027
和分段点
Figure FDA0002458389690000028
组成的分段点序列Pi；
1.2.4.以Wk-1窗口的聚类中心为初始聚类中心，采用K-means聚类算法对三元组序列
Figure FDA0002458389690000029
进行聚类，并为不同类的线段分配不同符号得到符号化序列
Figure FDA00024583896900000210
其中聚类过程中采用基于欧氏距离的指标描述两条线段si和sj之间的相似度：
Figure FDA00024583896900000211
其中，dij即表示线段si和sj的相似度，dij越小，则表示两条线段有更相似的变化形态，ωk和ωr表示权重；
1.2.5.对于两个运行参数Vi和Vj，分别得到符号化后的序列
Figure FDA00024583896900000212
和
Figure FDA00024583896900000213
合并两个运行参数的分段点序列Pi和Pj，并按合并后的分段点对其符号化序列
Figure FDA00024583896900000214
和
Figure FDA00024583896900000215
进行分割重构，得到重构后的符号化序列
Figure FDA00024583896900000216
和
Figure FDA00024583896900000217
由其构成事务集，nij-1为Pi和Pj合并后的分段点个数。
2.根据权利要求1所述的一种基于动态关联规则挖掘的工业过程异常工况预测方法，其特征在于：所述的步骤2具体通过如下子步骤实现：
步骤2.1：由步骤1获得的事务集的每一事务记为
Figure FDA00024583896900000218
另将
Figure FDA00024583896900000219
和
Figure FDA00024583896900000220
中所包含的线段类别符号分别记为
Figure FDA00024583896900000221
和
Figure FDA00024583896900000222
mi和mj分别为
Figure FDA0002458389690000031
和
Figure FDA0002458389690000032
中所包含的线段类别数；设频繁项集产生过程的最小支持度阈值为minsup1；
步骤2.2：对于初始数据窗口W1进行关联规则挖掘，其具体过程如下：
2.2.1)对所有传感器的传感器测量序列数据集进行扫描以计算每个项的支持度，记σ(·)表示某项或项集的支持度计数，初始为0；假设
Figure FDA0002458389690000033
的类别符号为tk，e表示i或j，t表示c或d；
对每一个
Figure FDA0002458389690000034
计算σ(tk)＝σ(tk)+1；
对每一个tk，判断
Figure FDA0002458389690000035
是否成立，若成立则认为tk为频繁1-项集，保留tk并记录相应的支持度计数σ(tk)；
2.2.2)使用已得到的频繁1-项集tk构成2-项集，并发现频繁2-项集，分别记cp和dq分别为经过步骤2.2.1)从原线段类别符号
Figure FDA0002458389690000036
和
Figure FDA0002458389690000037
中保留的项；
对于每一个{cp,dq}执行操作：对每一个存在于
Figure FDA0002458389690000038
中的{cp,dq}，计算σ({cp,dq})＝σ({cp,dq})+1；如果
Figure FDA0002458389690000039
不小于minsup1，则认为{cp,dq}为频繁2-项集，保留{cp,dq}并记录相应的支持度计数；
2.2.3)使用步骤2.2.2中获得的频繁2-项集{cp,dq}计算每两个运行参数Vi和Vj在整个数据集中的支持度，计算过程如下：对每两个运行参数Vi和Vj构成的项集{Vi,Vj}，计算σ({Vi,Vj})＝sum(σ({cp,dq}))，两个参数的关联规则支持度
Figure FDA00024583896900000310
并计算σ(Vi)＝sum(σ(cp))；σ(Vj)＝sum(σ(dq))；
2.2.4)对于每组{Vi,Vj}，产生如下关联规则：Vj→Vi和Vi→Vj；对于每一个关联规则，计算其置信度为
Figure FDA00024583896900000311
步骤2.3：对于每一个新的数据窗口，采用以下过程进行频繁项集产生并生成关联规则：
对每一个k>1的数据窗口Wk执行步骤2.3.1)～2.3.3)：
2.3.1)对每一滑出窗口的事务
Figure FDA0002458389690000041
记为{ca1,da2}，判断{ca1,da2}是否为频繁项集，若是则更新σ({ca1,da2})＝σ({ca1,da2})-1；
2.3.2)对每一滑入窗口的事务
Figure FDA0002458389690000042
记为{cb1,db2}，判断{cb1,db2}是否为频繁项集，若是，则更新σ({cb1,db2})＝σ({cb1,db2})+1，若否，则在数据窗口Wk内计算该{cb1,db2}的支持度计数为σ({cc1,dc2})；
2.3.3)对每一个{ca1,da2}，若
Figure FDA0002458389690000043
则更新{ca1,da2}相应项集的支持度计数，若不满足则剔除相应项集；对每一个{cb1,db2}，若
Figure FDA0002458389690000044
则更新{cb1,db2}相应项集的支持度计数，若不满足则剔除相应项集；对每一个{cc1,dc2}，若
Figure FDA0002458389690000045
则更新{cc1,dc2}相应项集的支持度计数，若不满足则剔除相应项集；
2.3.4)利用步骤2.2.3)和2.2.4)更新每一关联规则的支持度和置信度。
3.根据权利要求2所述的一种基于动态关联规则挖掘的工业过程异常工况预测方法，其特征在于：所述的步骤3具体通过如下子步骤实现：
记预设的预测步长为lp，由训练数据集挖掘到的关联规则中提取的一组关联参数为{V1,V2,...,Vu}，u为提取出的一组关联参数个数；对于每一参数Vi，其测量值序列为
Figure FDA0002458389690000046
构造以下矩阵Itrain为初始的神经网络训练输入：
Figure FDA0002458389690000047
其中，Itrain中每一列为一个训练输入样本；构造训练输出Otrain为，
Figure FDA0002458389690000048
利用公式(1)构造的训练输入样本，对小波神经网络进行训练；在网络初始化时，利用关联规则Vi→Vu对应的初始置信度ωi设置网络输入层和隐含层之间的初始权值，i＝1,2,...u-1。
4.根据权利要求3所述的一种基于动态关联规则挖掘的工业过程异常工况预测方法，其特征在于：所述的步骤4具体通过如下子步骤实现：
步骤4.1：关联规则更新做法为：记在设备运行过程中传感器已采集到的c个新数据为
Figure FDA0002458389690000051
利用新采集到的数据中最近的L个数据
Figure FDA0002458389690000052
通过步骤2.3更新关联规则，从而更新关联规则Vi→Vu所对应的置信度记为
Figure FDA0002458389690000053
步骤4.2：模型更新做法为：结合训练数据以及新采集到的数据构造新的模型训练样本，即数据集为
Figure FDA0002458389690000054
构造矩阵Itest为神经网络更新的输入，构造Otest为输出，
Figure FDA0002458389690000055
Figure FDA0002458389690000056
利用上述公式构造新的训练样本，对小波神经网络进行训练；在网络初始化时，利用新的置信度
Figure FDA0002458389690000057
设置网络输入层和隐含层之间的初始权值。
5.根据权利要求4所述的一种基于动态关联规则挖掘的工业过程异常工况预测方法，其特征在于：所述的步骤5具体通过如下子步骤实现：
记预设的异常工况发生阈值为ωp，对于最新采集到的数据，利用步骤4中更新的模型进行lp步预测，若得到的目标参数预测值相对于初始正常的漂移量超过所设阈值，则认为异常工况发生。
6.根据权利要求5所述的一种基于动态关联规则挖掘的工业过程异常工况预测方法，其特征在于：在预测过程结束前，每更新预定数量的数据，就回到步骤4进行关联规则和模型的更新。
CN201910244877.8A 2019-03-28 2019-03-28 一种基于动态关联规则挖掘的工业过程异常工况预测方法 Active CN110018670B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910244877.8A CN110018670B (zh) 2019-03-28 2019-03-28 一种基于动态关联规则挖掘的工业过程异常工况预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910244877.8A CN110018670B (zh) 2019-03-28 2019-03-28 一种基于动态关联规则挖掘的工业过程异常工况预测方法

Publications (2)

Publication Number Publication Date
CN110018670A CN110018670A (zh) 2019-07-16
CN110018670B true CN110018670B (zh) 2020-07-10

Family

ID=67190170

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910244877.8A Active CN110018670B (zh) 2019-03-28 2019-03-28 一种基于动态关联规则挖掘的工业过程异常工况预测方法

Country Status (1)

Country Link
CN (1) CN110018670B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111060652A (zh) * 2019-11-21 2020-04-24 西南交通大学 基于长短期记忆网络的变压器油中溶解气体浓度预测方法
CN113515554A (zh) * 2020-04-09 2021-10-19 华晨宝马汽车有限公司 用于不规则采样的时间序列的异常检测方法和***
CN111858712A (zh) * 2020-07-20 2020-10-30 上海仪电（集团）有限公司中央研究院 原位水质巡检数据时空分析与异常检测方法和***
CN112016097B (zh) * 2020-08-28 2024-02-27 深圳泓越信息科技有限公司 一种预测网络安全漏洞被利用时间的方法
CN112329828A (zh) * 2020-10-26 2021-02-05 北京旋极信息技术股份有限公司 一种故障关联分析方法和装置
CN112380274B (zh) * 2020-11-16 2023-08-22 北京航空航天大学 一种面向控制过程的异常检测方法
CN112632151B (zh) * 2020-12-25 2023-02-10 南京邮电大学 一种基于二分过滤的移动对象传染行为挖掘方法
CN112862347A (zh) * 2021-03-02 2021-05-28 同济大学 基于联邦学习的设备异常监测方法、***、存储介质及终端
CN112800686A (zh) * 2021-03-29 2021-05-14 国网江西省电力有限公司电力科学研究院 一种变压器dga在线监测数据异常模式判定方法
CN112801426B (zh) * 2021-04-06 2021-06-22 浙江浙能技术研究院有限公司 一种基于关联参数挖掘的工业过程故障融合预测方法
CN112749370B (zh) * 2021-04-06 2021-07-02 广东际洲科技股份有限公司 一种基于物联网的故障跟踪定位方法和***
CN113297272B (zh) * 2021-05-30 2023-01-06 福建中锐网络股份有限公司 一种桥梁监测数据关联规则挖掘及健康预警方法及***
CN113591994B (zh) * 2021-08-03 2023-06-06 北京邮电大学 一种基于自动标注的终端行为预测方法
CN115689071B (zh) * 2023-01-03 2023-05-02 南京工大金泓能源科技有限公司 一种基于关联参数挖掘的设备故障融合预测方法及***

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101408754A (zh) * 2008-10-30 2009-04-15 中山大学 基于数据挖掘的智能家居优化***
CN101667197A (zh) * 2009-09-18 2010-03-10 浙江大学 基于滑动窗口的数据流关联规则挖掘方法
CN201898519U (zh) * 2010-09-01 2011-07-13 燕山大学 带有风险控制的设备维修预警装置
CN102818642A (zh) * 2012-07-18 2012-12-12 辽宁省海洋水产科学研究院 刺参病害预警***
CN103440248A (zh) * 2013-07-22 2013-12-11 西南交通大学 基于动态关联规则的网络视频事件挖掘框架
CN103676645A (zh) * 2013-12-11 2014-03-26 广东电网公司电力科学研究院 一种时间序列数据流中的关联规则的挖掘方法
CN104636479A (zh) * 2015-02-15 2015-05-20 西安电子科技大学 一种工业大数据驱动的总完工时间预测方法
CN108873859A (zh) * 2018-05-31 2018-11-23 浙江工业大学 基于改进关联规则的桥式抓斗卸船机故障预测模型方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7474987B2 (en) * 2004-02-10 2009-01-06 University Of Maryland Method and system for optimal data diagnosis

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101408754A (zh) * 2008-10-30 2009-04-15 中山大学 基于数据挖掘的智能家居优化***
CN101667197A (zh) * 2009-09-18 2010-03-10 浙江大学 基于滑动窗口的数据流关联规则挖掘方法
CN201898519U (zh) * 2010-09-01 2011-07-13 燕山大学 带有风险控制的设备维修预警装置
CN102818642A (zh) * 2012-07-18 2012-12-12 辽宁省海洋水产科学研究院 刺参病害预警***
CN103440248A (zh) * 2013-07-22 2013-12-11 西南交通大学 基于动态关联规则的网络视频事件挖掘框架
CN103676645A (zh) * 2013-12-11 2014-03-26 广东电网公司电力科学研究院 一种时间序列数据流中的关联规则的挖掘方法
CN104636479A (zh) * 2015-02-15 2015-05-20 西安电子科技大学 一种工业大数据驱动的总完工时间预测方法
CN108873859A (zh) * 2018-05-31 2018-11-23 浙江工业大学 基于改进关联规则的桥式抓斗卸船机故障预测模型方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
一种面向时空数据的关联规则更新算法;刘伯红等;《计算机与数字工程》;20151020(第10期);第1767-1770、1774页 *
基于关联规则的故障诊断方法及研究;张建明等;《化工自动化及仪表》;20031030(第05期);第11-14、18页 *
基于关联规则的火电厂传感器故障检测;鲍文;《中国电机工程学报》;20031231;第170-174页 *
基于滑动窗口的多元时间序列数据动态关联规则挖掘;熊腾飞;《中国优秀硕士学位论文全文数据库信息科技辑》;20170215;第16-38页 *
基于神经网络的铝电解混合控制模型研究;汪成亮等;《计算机应用研究》;20100715(第07期);第2536-2539、2559页 *

Also Published As

Publication number Publication date
CN110018670A (zh) 2019-07-16

Similar Documents

Publication Publication Date Title
CN110018670B (zh) 一种基于动态关联规则挖掘的工业过程异常工况预测方法
CN110008253B (zh) 一种工业数据关联规则挖掘及异常工况预测方法
CN110263846B (zh) 基于故障数据深度挖掘及学习的故障诊断方法
CN111694879B (zh) 一种多元时间序列异常模式预测方法及数据采集监控装置
US9779361B2 (en) Method for learning exemplars for anomaly detection
CN102789545B (zh) 基于退化模型匹配的涡轮发动机剩余寿命的预测方法
CN110895526A (zh) 一种大气监测***中数据异常的修正方法
CN108985380B (zh) 一种基于聚类集成的转辙机故障识别方法
US20210216386A1 (en) Time-sequential data diagnosis device, additional learning method, and recording medium
KR20190072652A (ko) 정보 처리 장치 및 정보 처리 방법
CN114167838B (zh) 一种伺服***多尺度健康评估与故障预测方法
CN111931806A (zh) 一种多传感器数据融合的设备故障诊断方法及装置
CN113723010A (zh) 一种基于lstm温度-位移相关模型的桥梁损伤预警方法
CN109298633A (zh) 基于自适应分块非负矩阵分解的化工生产过程故障监测方法
CN112949735A (zh) 一种基于离群数据挖掘的液态危化品挥发浓度异常发现方法
CN115186762A (zh) 一种基于dtw-knn算法的发动机异常检测方法及***
Zhang et al. Remaining useful life prediction for rolling bearings with a novel entropy-based health indicator and improved particle filter algorithm
CN114036647A (zh) 一种基于实车数据的动力电池安全风险评估方法
CN111144018B (zh) 一种基于航后数据的航空发动机整机剩余性能提取方法
CN110956155B (zh) 基于co数据的综采工作面作业工序模糊聚类识别方法
CN110032585B (zh) 一种时间序列双层符号化方法及装置
CN115935285A (zh) 基于掩码图神经网络模型的多元时间序列异常检测方法和***
CN115659271A (zh) 传感器异常检测方法、模型训练方法、***、设备及介质
Li et al. A similarity-based and model-based fusion prognostics framework for remaining useful life prediction
CN112100919B (zh) 一种基于re-cf-ekf算法的滚动轴承剩余寿命预测方法

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
GR01 Patent grant
GR01 Patent grant