CN114647231A - 一种数据驱动的工业过程故障检测方法 - Google Patents
一种数据驱动的工业过程故障检测方法 Download PDFInfo
- Publication number
- CN114647231A CN114647231A CN202210367561.XA CN202210367561A CN114647231A CN 114647231 A CN114647231 A CN 114647231A CN 202210367561 A CN202210367561 A CN 202210367561A CN 114647231 A CN114647231 A CN 114647231A
- Authority
- CN
- China
- Prior art keywords
- data
- matrix
- vector
- characteristic
- fault detection
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 69
- 238000001514 detection method Methods 0.000 title claims abstract description 60
- 239000011159 matrix material Substances 0.000 claims abstract description 214
- 238000000034 method Methods 0.000 claims abstract description 131
- 230000008569 process Effects 0.000 claims abstract description 65
- 238000012880 independent component analysis Methods 0.000 claims abstract description 43
- 238000010219 correlation analysis Methods 0.000 claims abstract description 30
- 238000000605 extraction Methods 0.000 claims abstract description 14
- 238000012847 principal component analysis method Methods 0.000 claims abstract description 13
- 238000009776 industrial production Methods 0.000 claims abstract description 5
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 115
- 238000012360 testing method Methods 0.000 claims description 66
- 238000012549 training Methods 0.000 claims description 37
- 238000000513 principal component analysis Methods 0.000 claims description 28
- 230000002159 abnormal effect Effects 0.000 claims description 25
- 230000006870 function Effects 0.000 claims description 18
- 238000013527 convolutional neural network Methods 0.000 claims description 12
- 230000004913 activation Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000013528 artificial neural network Methods 0.000 claims description 3
- 230000002087 whitening effect Effects 0.000 claims description 3
- 238000002372 labelling Methods 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 abstract description 4
- 238000013135 deep learning Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 12
- 230000006872 improvement Effects 0.000 description 8
- AOSFMYBATFLTAQ-UHFFFAOYSA-N 1-amino-3-(benzimidazol-1-yl)propan-2-ol Chemical compound C1=CC=C2N(CC(O)CN)C=NC2=C1 AOSFMYBATFLTAQ-UHFFFAOYSA-N 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 239000007789 gas Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 239000000498 cooling water Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000012263 liquid product Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000012495 reaction gas Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- 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/20—Pc systems
- G05B2219/24—Pc safety
- G05B2219/24065—Real time diagnostics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Test And Diagnosis Of Digital Computers (AREA)
Abstract
本发明公开了一种数据驱动的工业过程故障检测方法,包括获取工业生产的过程数据并进行中心化的预处理,然后将获得的中心化的矩阵分别经过主成分分析法、独立成分分析法和典型关联分析法提取特征,再将所提取的特征经过工业过程故障检测网络,推断得到过程数据的所对应的故障类型的分类标签。本发明通过数据驱动的方法来进行故障诊断,对比基于知识和模型的方法,通用性更强,建模过程通过手动特征提取后进行深度学习网络完成,能更好的适配规模较大和复杂的工业过程。
Description
技术领域
本发明涉及工业过程的控制与故障诊断领域,具体是一种数据驱动的工业过程故障检测方法。
背景技术
近年来,随着现代工业的不断发展,形成了许多规模巨大的工业生产制造系统和高度复杂的工业过程。在现代工业通过不断扩大和复杂化的工业过程为其带来便利的同时,也带来了一些问题和挑战。由于系统日趋复杂和大规模,生产制造过程中机器发生故障的概率大大增加,故障形成后容易在系统中扩散。若没有一个可靠的检测手段,可能会导致设备损坏,人员伤亡,产品不合格等严重后果。
针对大规模高复杂度的故障,目前检测方法一般有三类。第一类是基于知识的方法,优点为比较适用于当建立数学模型的过程比较容易的时候,这个方法被局限于那些过程变量较少的系统。第二类是基于模型的方法,特点是不需要精确的数学模型,只需要一些收集的基础数据,加入了一些固定的知识来进行检测但是通用性较差,局限性突出。第三类就是基于数据驱动的方法,特点是需要数据是运行的状态信息,通过基础数据,使用目前比较丰富的理论知识处理得到结果。
目前数据驱动的过程监测和故障检测是研究的热点,现在存在的方法主要存在训练网络复杂数据维度高,特征工程不够全面,导致检测单故障准确率高,多故障准确率低,计算量大等问题。因此需要改进现有技术从而能够准确快速地分析工业过程中的故障数据,保障工业生产稳定运行。
发明内容
本发明要解决的技术问题是提供一种数据驱动的工业过程故障检测方法,用以解决目前工业过程中检测单一故障准确率高,多故障检测准确率低,数据维度高计算量大的问题。
为了解决上述技术问题,本发明提供一种数据驱动的工业过程故障检测方法,包括的过程如下:
S01、获取工业生产的过程数据,过程数据的维度为53,总条数为m;然后构建一个m行53列的过程数据的原始矩阵Zsrc,将过程数据的原始矩阵Zsrc的每一列都减去训练好的列平均值矩阵Zmean,得到中心化的矩阵Zo;
S02、将中心化的矩阵Zo输入主成分分析法并利用训练好的特征向量矩阵A9计算特征数据Y1_online=A9ZO;将中心化的矩阵Zo输入独立成分分析法并利用训练好的矢量矩阵计算特征数据将中心化的矩阵Zo输入典型关联分析法并利用训练好的特征向量矩阵Σ9和训练好的右奇异向量正交基V9计算特征数据Y3_online:
其中,pk为中心化的矩阵Zo的第k行;k∈(1,m)为过程数据的当前时刻的条数,将k从1到m依次带入式(23)后获得特征数据Y3_online;
S03、工业过程故障检测网络包括PCA通道、ICA通道和CCA通道的三通道卷积神经网络,将特征数据Y1_online作为PCA通道的输入,特征数据Y2_online作为ICA通道的输入,特征数据Y3_online作为CCA通道的输入,分别输入到工业过程故障检测网络中,推断得到过程数据的分类标签,分类标签为1-6属于推断结果异常,将数据打上异常标签存入数据库;如果分类标签为0属于推断结果正常,则将数据打上正常标签存入数据库。
作为本发明的一种数据驱动的工业过程故障检测方法的改进:
步骤S02中所述主成分分析进行特征提取的过程为:
S0201、主成分分析法提取特征:
S020101、根据输入的中心化的矩阵求出协方差矩阵C:
其中,Z为输入的中心化的矩阵;
S020102、求出协方差矩阵C的特征值λi和对应的特征向量;
S020103、将协方差矩阵C的特征值λi从大到小排列作为主对角线的值,其它值全部为0,保存为特征值矩阵Λ;设定η>85%,根据式(3)来确定选取的特征值的个数k:
S020104:将特征向量按照对应的特征值λi大小,按从大到小的顺序从上到下按行排列成矩阵,保存为特征向量矩阵A,取特征向量矩阵A的前k行记为Ak;
S020105:所述主成分分析后的特征数据Y1:
Y1=AkZ (4)
经过主成分分析法的输出为特征数据Y1和特征向量矩阵Ak。
作为本发明的一种数据驱动的工业过程故障检测方法的进一步改进:
步骤S02中所述独立成分分析法进行特征提取的过程为:
S0202、独立成分分析法提取特征
S020201、获取白化矩阵W:
S020202:获取新的白化后的过程数据
Znew=WZ (6)
其中,Z为输入的中心化的矩阵;
S020203:令需要估计的分量个数为r,迭代次数为n,初始化n=0且选择随机的初始向量ω(0);
S020204:通过牛顿迭代法,按式(7)进行计算ω(n+1):
ω(n+1)=E{Znewg(ω(n)TZnew)}-E{g’(ωTZnew)}ω(n) (7)
其中,E表示均值,g代表非线性函数,a表示常数,g’(ωTZnew)表示函数g的导数,n为计算次数,ω(n)表示通过迭代法计算了n次的结果;
S020205、对得到的ω(n+1)按式(8)进行归一化,并按式(9)进行收敛判断:如果不满足式(9)的条件,那么n=n+1后,重复步骤S020204计算ω(n+1);
||ω(n+1)ω(n)T|-1|<ε (9)
其中,ε为误差;
S020206、所述独立成分分析提取的特征数据Y2为:
作为本发明的一种数据驱动的工业过程故障检测方法的进一步改进:
步骤S02中所述典型关联分析法进行特征提取的过程为:
S0203、典型关联分析法提取特征
S020301、构建历史向量pk与历史数据集Pk:
pk=[zk-1,zk-2,…,zk-h]T (12)
Pk=[pk-1,pk-2,…,pk+N-1]∈Rmh×N (13)
其中,h为历史向量长度,N为数据总条数,zk为输入的中心化的矩阵第k行数据,并记Pk的协方差矩阵为Σpp;
S020302:构建未来向量fk与未来数据集Fk:
fk=[zk,zk+1,…,zk+s]T (14)
Fk=[fk,fk+1,…,fk+N-1]∈Rsh×N (15)
其中,s为未来向量长度;
记未来数据集Fk的协方差矩阵为Σff,记未来向量fk与未来数据集Fk的协方差为Σpf;
S020303、使用典型关联分析法使历史数据集Pk和未来数据集Fk构成的数据集相关性最大,对相关系数进行奇异值分解:
ρpf=UΣVT (16)
根据奇异值分解公式可得:
其中,Σff为未来数据集Fk的协方差矩阵,Σpf是历史数据集Pk和未来数据集Fk的协方差,Σpp是历史数据集Pk的协方差矩阵,V为右奇异向量正交基,U为左奇异向量正交基,Σ为特征向量矩阵;可得:
Σ=Σpf (19)
S020304、所述典型关联分析法得到的特征数据Y3为:
其中,y(k)表示k时刻经过典型关联分析法得到的特征数据,k∈(1,m)表示输入的中心化的矩阵的当前时刻的条数,m为输入的中心化的矩阵的的总条数,将k从1到m带入式(21),得出经过所述典型关联分析得到的特征数据Y3;
经过典型关联分析法的输出为特征数据Y3、特征向量矩阵Σ和右奇异向量正交基V。
作为本发明的一种数据驱动的工业过程故障检测方法的进一步改进:
步骤S03中的所述工业过程故障检测网络包括三个等价通道的三通道卷积神经网络,分别为PCA通道、ICA通道和CCA通道;每个通道均包括四层结构:第一层为16个2x2的卷积核,步长为1,padding为1,激活函数为ReLU,输出为16x4x4;第二层为32个2x2的卷积核,步长为1,padding为0,激活函数也为ReLU,输出为32x3x3;第三层为16个2x2的卷积核,步长为1,padding为0,激活函数同为ReLU,输出为16x2x2;第四层为全连接层;
输入的特征数据经过PCA通道、ICA通道和CCA通道后均获得32维数据,将三通道的32维数据进行连接,并在通过一层全连接层,输出32维数据,最终通过输出层输出0到6的分类结果。
作为本发明的一种数据驱动的工业过程故障检测方法的进一步改进:
1)、田纳西伊斯曼的过程数据进行软件仿真获得样本数据并进行标签,正常数据的标签为0,异常数据的标签为1-6;将样本数据及其对应的标签按照8:2划分为训练集和测试集;
2)、对训练集进行中心化的预处理后获得中心化的矩阵Ztrain,同时保存训练集每一列的平均值为所述训练好的列平均值矩阵Zmean=zmean_1,zmean_2,…,zmean_53;
3)、中心化的矩阵Ztrain通过所述主成分分析法的获得特征向量矩阵A,取特征向量矩阵A的前9行并保存为所述训练好的特征向量矩阵A9,输出特征数据Y1=A9Ztrain;
5)、中心化的矩阵Ztrain通过所述经过典型关联分析法获得特征向量矩阵Σ和右奇异向量正交基V,然后维度取n=9,保存特征向量矩阵Σ的前9×9对角阵为所述训练好的特征向量矩阵Σ9和右奇异向量正交基V的前9×9对角阵为所述训练好的右奇异向量正交基V9,将k从训练集的第一条到最后一条依次代入按式(23),输出特征数据Y3:
其中,pk为中心化的矩阵Ztrain的第k行。
作为本发明的一种数据驱动的工业过程故障检测方法的进一步改进:
所述测试集构建一个m行53列的测试集的原始矩阵Zsrc,将原始矩阵Zsrc的每一列都减去训练好的列平均值矩阵Zmean,得到中心化的矩阵Ztest,然后进行特征提取:
将中心化的矩阵Ztest通过主成分分析法,并利用训练好的特征向量矩阵A9计算特征数据Y1_test=A9Ztest;将中心化的矩阵Ztest通过独立成分分析法,并利用训练好的矢量矩阵计算特征数据将中心化的矩阵Ztest通过典型关联分析法,并利用训练好的特征向量矩阵Σ9和训练好的右奇异向量正交基V9,将k从训练集的第一条到最后一条依次代入式(23)后输出特征数据Y3_test:
其中,pk为中心化的矩阵Ztest的第k行。
作为本发明的一种数据驱动的工业过程故障检测方法的进一步改进:
所述中心化预处理为:将输入的数据构建一个m行53列的原始矩阵为Zsrc=[z1,z2,z3,z4……,z53],对53个维度进行归一化,操作为:
其中:zi_new是中心化之后的数据,zi是第i列的列向量,包括同一个维度的数据集合,mean(zi)表示zi所在列的平均值;
作为本发明的一种数据驱动的工业过程故障检测方法的进一步改进:
所述工业过程故障检测网络的训练过程为:所述特征数据Y1输入PCA通道、特征数据Y2输入ICA通道、特征数据Y3输入CCA通道进行模型训练,学习率设置为0.001和损失函数为交叉熵损;共训练100个epoch病保存每个epoch训练好的模型,共保存100个训练好的工业过程故障检测网络模型,每个模型输出为训练集每一条样本数据对应的分类标签,其中标签0为正常数据,标签1-6代表故障数据;
所述工业过程故障检测网络的测试过程为:将特征数据Y1_test作为PCA通道的输入,特征数据Y2_test作为ICA通道的输入,特征数据Y3_test作为CCA通道的输入,将测试集分别输入所述100个训练好的工业过程故障检测网络模型中测试,得出每个模型的测试的F1分数,取出分数最高的工业过程故障检测网络模型为在线使用的所述工业过程故障检测网络。
本发明的有益效果主要体现在:
1.本发明针对复杂流程进行分析,通过分析发现在一个复杂过程中产生的数据变化规律不仅包含高斯性也包含非高斯性,且在时序数据在时间维度上存在一定关联,因此,本发明通过使用三种特征提取方法,分别能够将信号的高斯特征,非高斯特征,自相关特征加入到数据特征中。并结合三通道神经网络进行三通道输入,结果发现能在多种分类过程中有效提升准确率,能够更好的对生产过程中的故障进行溯源。
2.本发明通过数据驱动的方法来进行故障诊断,对比基于知识和模型的方法,通用性更强,建模过程通过手动特征提取后进行深度学习网络完成,能更好的适配规模较大和复杂的工业过程。
3.本发明通过上面特征提取方式也完成了数据降维,同时使用一个层数较少的三层卷积神经网络作为模型,在基本不降低异常数据判断准确的前提下,能够使用更少的计算来完成。
附图说明
图1为本发明的一种数据驱动的工业过程故障检测方法的整体流程图;
图2为三通道卷积神经网络结构图;
图3为实验中的单通道卷积神经网络结构图。
具体实施方式
下面结合具体实施例对本发明进行进一步描述,但本发明的保护范围并不仅限于此:
实施例1、本发明针对工业过程故障检测模型复杂,计算量大,多故障检测准确率低等问题提出了一种数据驱动的工业过程故障检测方法。
离线建模的数据为田纳西伊斯曼(Tennessee Eastman,TE)过程中的数据,将TE过程中的数据进行采集,完成建模过程。然后将建好的模型应用于实际生产亚硫酸氢铵过程中进行故障检测。如图1所示,包括以下步骤,步骤1-步骤4为离线建模阶段,步骤5为在线处理阶段,具体为:
步骤1:对历史生产过程中的正常数据导入并打上为正常标签,故障数据同理,并将数据保存成方便处理的csv格式;
数据集使用的是TE过程的数据。TE过程是依据实际化工反应过程,由美国Eastman化学公司开发了具有开放性和挑战性的化工模型仿真平台,其产生的数据具有时变、强耦合和非线性特征,广泛用于测试复杂工业过程的控制和故障诊断模型。TE过程主要包括如下5个操作单元:反应器(Reactor)、冷凝器(Condenser)、气/液分离器(Separator)、压缩机(Compressor)和汽提塔(Stripper)。TE过程的化学反应过程涉及A、B、C、D、E、F、G和H八种成分;其中,反应物包括气态物质A、C、D、E和惰性催化剂B,生成物包括液态产品G、H和副产品F,整体反应流程为,反应气体A、C、D、E进入反应器经过催化剂进行反应,产物通过冷凝器冷凝进行冷却,冷却后的形成气液混合流后进入分离器。分离器通过分离作用,将得到的气体重新循环进入反应仪器中,在进入循环压缩机前需要排出蒸汽,从分离器得到的冷凝成分会被放入气体塔进行分离,经过汽提塔的反应成分通过循环后在回到反应器中。反应过后在汽提塔底部获取生产得到生产产物G和H。过程中连续测量变量如表1:
表1:连续测量变量
成分分析变量如表2:
表2:成分分析变量
控制变量如表3:
表3:控制变量
变量描述 | 编号 |
D进料量(流2) | XMv(1) |
E进料量(流3) | XMv(2) |
A进料量(流l) | XMv(3) |
总进料量(流4) | XMv(4) |
压缩机再循环阀 | XMv(5) |
排放阀(流9) | XMv(6) |
分离器罐液流量(流10) | XMv(7) |
汽提塔液体产品流量(流II) | XMv(8) |
汽提塔水流阀 | XMv(9) |
反应器冷却水流量 | XMv(10) |
冷凝器冷却水流量 | XMv(11) |
搅拌速度 | XMv(12) |
故障类型如表4:
表4:故障类型
TE过程的数据维度一共53维,包含21种故障的异常情况。数据均通过Matlab仿真获得。
主要做了以下几种仿真:
1.针对正常情况进行73h仿真,获得了1460条正常数据。
2.针对21种异常情况进行仿真,每种异常情况均运行24h,获取到21种异常情况、每种异常情况为480条数据。
3.针对先正常情况后异常情况,对21种异常情况进行仿真,正常情况运行8h,异常情况运行40h,可以获得160条正常数据和800条对应的异常数据。如此,共进行了21次正常后异常的仿真,共获得正常数据是160*21条,21种异常数据每种各有800*21条。
将上述数据按表4的“分类”栏的值进行标签:正常数据的标签为0,异常数据按照表4“分类”栏进行分组合并,例如出现编号为8的异常数据打上的标签为2。样本总共1460+21*(480+960)=31700条数据,其中正常数据4820条,21种异常数据每种各1280条。
步骤2:对数据进行预处理,即进行中心化处理;
将输入的数据构建一个m行53列的原始矩阵为Zsrc=[z1,z2,z3,z4……,z53],原始矩阵Zsrc的数据行数为m对应于输入的数据的条数,代表时间,原始矩阵Zsrc的列和输入数据的维度相对应,代表着温度、含量、压力等信息;对53个维度进行归一化,操作为:
其中:zi_new是中心化之后的数据,zi是第i列的列向量,包括同一个维度的数据集合,mean(zi)表示zi所在列的平均值,i∈(1,53);
步骤3:采用主成分分析、独立成分分析和典型关联分析三种方法分别对中心化处理后的数据进行分析,获得三种不同特征矩阵作为下一步所建立的工业过程故障检测网络的输入;
步骤3-1:主成分分析法提取特征:
步骤3-1-1:根据步骤2中的中心化的矩阵Z,求出Z的协方差矩阵C:
步骤3-1-2:求出协方差矩阵C的特征值λi和对应的特征向量;
步骤3-1-3:将协方差矩阵C的特征值λi从大到小排列,并记为λ=[λ1,λ2,λ3,λ4,λ5……λm];将从大到小排列的特征值λi作为主对角线的值,其它值全部为0,保存为特征值矩阵记作Λ,选取贡献率η超过一个百分比,一般选取η>85%,根据式(3)来确定选取的特征值的个数k:
步骤3-1-4:将特征向量按照对应的特征值λi大小,按从大到小的顺序从上到下按行排列成矩阵,保存为特征向量矩阵A,取特征向量矩阵A的前k行记为Ak;
步骤3-1-5:计算主成分分析后的数据:
Y1=AkZ (5)
其中,Y1为主成分分析后获得的特征数据,Z为中心化的矩阵;
中心化的矩阵Z经过主成分分析法的输出为特征数据Y1和特征向量矩阵Ak;
步骤3-2:使用独立成分分析法提取特征,并保存特征矩阵;
步骤3-2-1:获取白化矩阵W,获取方法为:
其中Λ为步骤3-1-3中的特征值矩阵,A为步骤3-1-5中的特征向量矩阵;
步骤3-2-2:获取新的白化后的过程数据Znew,
Znew=WZ (7)
Z为中心化的矩阵;
步骤3-2-3:令需要估计的分量个数为r,迭代次数为n,初始化n=0且选择随机的初始向量ω(0);
步骤3-2-4:通过牛顿迭代法,按式(7)进行计算:
ω(n+1)=E{Znewg(ω(n)TZnew)}-E{g’(ωTZnew)}ω(n) (8)
其中,E表示均值,g代表非线性函数,可以取g1(x)=tanhax,g2(x)=x3等,a表示常数,g’(ωTZnew)表示函数g的导数,n为计算次数,ω(n)表示通过迭代法计算了n次的结果;
步骤3-2-5:对得到的ω(n+1)按式(8)进行归一化(第一次迭代是代入ω(0),n取0,得到的是ω(1)),并按式(9)进行收敛判断:
如果不满足式(9)的条件,那么n=n+1后,重复步骤3-2-4计算ω(n+1);
||ω(n+1)ω(n)T|-1|<ε (10)
其中,ε为误差;
步骤3-2-6:计算独立成分分析后的特征数据:
其中,Y2为独立成分分析后获得的特征数据,Z为步骤2中心化过后的矩阵
步骤3-3:使用典型关联分析法提取特征,并保存特征矩阵;
步骤3-3-1:选取历史向量长度为h,zk为步骤2中的中心化之后的数据zi_new的第k行,构建历史向量pk与历史矩阵Pk:
pk=[zk-1,zk-2,…,zk-h]T (13)
Pk=[pk-1,pk-2,…,pk+N-1]∈Rmh×N (14)
其中,N为数据总条数,并记Pk的协方差矩阵为Σpp;
步骤3-3-2:选取未来向量长度为s,构建未来向量fk与未来矩阵Fk:
fk=[zk,zk+1,…,zk+s]T (15)
Fk=[fk,fk+1,…,fk+N-1]∈Rsh×N (16)
记未来矩阵Fk的协方差矩阵为Σff,记历史矩阵Pk与未来矩阵Fk的协方差为Σpf;
步骤3-3-3:使用典型关联分析法使历史矩阵Pk和未来矩阵Fk构成的数据集相关性最大,对相关系数进行奇异值分解:
ρpf=UΣVT (17)
根据奇异值分解公式可得:
其中,Σff为未来矩阵Fk的协方差矩阵,Σpf是历史矩阵Pk和未来矩阵Fk的协方差,Σpp是历史数据集Pk的协方差矩阵,V为右奇异向量正交基,U为左奇异向量正交基,Σ为特征向量矩阵。
Σ=Σpf (20)
步骤3-3-4:典型关联分析法得到的数据
k∈(1,m)为输入数据的当前时刻的条数,m为输入数据的总条数,将k从1到m带入式(21),得出经过所述典型关联分析得到的特征数据Y3;
中心化的矩阵Z经过典型关联分析法的输出为特征数据Y3、特征向量矩阵Σ和右奇异向量正交基V。
步骤4:建立并训练工业过程故障检测网络
步骤4-1:构建工业过程故障检测网络
构建三通道卷积神经网络作为本发明的工业过程故障检测网络的模型,三通道卷积神经网络模型结构如图2所示,包括PCA通道、ICA通道和CCA通道。PCA通道、ICA通道和CCA通道的卷积核是相同的,所以三个通道是等价的。每个通道均包括四层结构:第一层为16个2x2的卷积核,步长为1,padding为1,激活函数为ReLU,输出为16x4x4;第二层为32个2x2的卷积核,步长为1,padding为0,激活函数也为ReLU,输出为32x3x3;第三层为16个2x2的卷积核,步长为1,padding为0,激活函数同为ReLU,输出为16x2x2;第四层为全连接层,将16x2x2的数据压平到32维。输入的特征数据分别经过PCA通道、ICA通道和CCA通道后均获得32维数据,将三通道的32维数据进行连接,并在通过一层全连接层,输出32维数据,最终通过输出层输出7种数据即0到6的分类结果,其中分类结果0为正常,其他故障如表4标签。
其中,激活函数:ReLU(x)=max(0,x) (23)。
步骤4-2:模型训练
将步骤1中所获得的31700条样本数据及其对应的标签按照8:2划分为训练集和测试集,并使得正常数据和异常数据均按8:2的比例分布于训练集和测试集,所得训练集有25360条样本数据,测试集一共6340条样本数据。
首先对训练集按步骤2进行中心化的预处理获得中心化的矩阵Ztrain,同时保存训练集每一列的平均值为训练好的列平均值矩阵Zmean=[zmean_1,zmean_1,…,zmean_53],然后提取特征:
1)、中心化的矩阵Ztrain通过步骤3-1主成分分析法的获得特征向量矩阵A,取特征向量矩阵A的前9行并保存为训练好的特征向量矩阵A9,输出特征数据Y1=A9Ztrain;
3)、中心化的矩阵Ztrain通过步骤3-3经过典型关联分析法获得特征向量矩阵Σ和右奇异向量正交基V,然后维度取n=9,保存特征向量矩阵Σ的前9×9对角阵Σ9(即训练好的特征向量矩阵Σ9)和右奇异向量正交基V的前9×9对角阵V9(即训练好的右奇异向量正交基V9),式(21)更新为:
其中,pk为中心化的矩阵Ztrain的第k行,将k从1到25360代入式(23)后输出特征数据Y3。
将上述中得到的训练集的特征数据Y1输入PCA通道、特征数据Y2输入ICA通道、特征数据Y3输入CCA通道进行模型训练,学习率设置为0.001和损失函数为交叉熵损;训练集所有数据都训练过一次为一个epoch,训练100个epoch,并保存每个epoch训练好的模型,共保存100个模型,模型输出为训练集每一条样本数据对应的分类标签,其中标签0为正常数据,标签1-6代表故障数据。
步骤4-3:模型测试
首先对测试集按步骤2进行中心化处理:将输入的测试集构建一个m行的原始矩阵Zsrc,将原始矩阵Zsrc的每一列都减去步骤4-2中保存的训练好的列平均值矩阵Zmean,即测试集的原始矩阵Zsrc第一列减去zmean_1,以此类推,得到中心化的矩阵Ztest,然后进行特征提取:
1):将中心化的矩阵Ztest通过主成分分析法,并利用训练好的特征向量矩阵A9计算特征数据Y1_test=A9Ztest;
3):将中心化的矩阵Ztest通过典型关联分析法,并利用训练好的特征向量矩阵Σ9和训练好的右奇异向量正交基V9计算特征数据Y3_test:
其中,pk为中心化的矩阵Ztest的第k行,将k从1到6340依次代入式(23)后输出特征数据Y3_test。
将上述提取的特征数据Y1_test作为PCA通道的输入,特征数据Y2_test作为ICA通道的输入,特征数据Y3_test作为CCA通道的输入,将测试集分别输入步骤4-2保存的100个模型中测试,得出每个模型的测试的F1分数,取出分数最高的模型保存为可以在线使用的工业过程故障检测网络。
其中F1分数的计算方法为:
其中,α为参数,取α=1,得到F1指标,P为准确率,R为召回率:
其中,TP代表将正类预测为正类数,TN代表将负类预测为负类数,FN代表将负类预测为正类数,FP代表将正类预测为负类数。
步骤5:在线处理工业过程故障检测及预警
在线处理与离线处理不同,离线处理是保存好的数据进行离线模型训练,是一个建模的过程,在线处理是实际运用的过程。通过在某化工厂生产亚硫酸氢铵中的过程来验证本模型的有效性。通过获取生产亚硫酸氢铵过程的数据,然后根据步骤4中训练得到的相关矩阵及模型,输出生产亚硫酸氢铵过程数据的分类标签,通过表4中的“分类”栏查询获得对应的故障“描述”栏。
步骤5-1:通过OPC UA协议获取生产亚硫酸氢铵中的实时过程数据,由于工厂已有采集传感设备,只需用OPC UA协议从中获取过程数据。采集的实时过程数据与步骤1中的TE过程的数据维度相同,也为53维,共采集5000条数据。
步骤5-2:在线处理方式与测试集一样,所有参与运算矩阵都来源于训练集计算。
首先对采集的实时过程数据按步骤2进行中心化处理:将输入的实时过程数据构建一个m行的原始矩阵Zsrc,将原始矩阵Zsrc的每一列都减去步骤4-2中保存的训练好的列平均值矩阵Zmean,即实时过程数据的原始矩阵Zsrc第一列减去列平均值矩阵zmean_1,以此类推,得到中心化的矩阵Zo,然后进行特征提取:
1)、将中心化的矩阵ZO通过主成分分析法,并利用训练好的特征向量矩阵A9计算特征数据Y1_online=A9ZO;
3)、将中心化的矩阵ZO通过典型关联分析法,并利用训练好的特征向量矩阵Σ9和训练好的右奇异向量正交基V9计算特征数据Y3_online:
其中,pk为中心化的矩阵Zo的第k行,将k从1到5000依次代入式(23)后输出特征数据Y3_online。
步骤5-3:将上述提取的特征数据Y1_online作为PCA通道的输入,特征数据Y2_online作为ICA通道的输入,特征数据Y3_online作为CCA通道的输入,分别输入到步骤4中获得的可在线使用的工业过程故障检测网络中,推断得到生产亚硫酸氢铵过程数据的分类标签;
步骤5-4:如果分类标签属于推断结果异常即判断标签不为0,则将当前数据,即原始收集到的生产亚硫酸氢铵过程数据53维数据打上异常标签存入数据库。如果分类标签属于推断结果正常即判断标签为0,则将数据打上正常标签存入数据库。
实验1:
为了验证本发明的有效性,为减少不同数据对最终分析结果的影响,统一采用实施例1中的实验数据来进行训练和测试,训练集有25360条数据,测试集有6340条数据,即与离线建模阶段所用的数据集和测试集一致。测试用机为macbook pro 2020m1,内存为16G。选取三种方案进行对比:
1)、基于主成分分析PCA的检测:
通过使用主成分分析法,将输入的训练集数据降维,根据T2统计量得出设置阈值,超出阈值的数据则判断为异常:
2)、使用常规卷积神经网络CNN
将训练集数据的中心化的矩阵Ztrain直接输入卷积神经网络CNN来进行网络训练,从而获得训练好的故障预警模型;
3)、单通道卷积神经网络
处理过程与本方法步骤完全一致,仅保留一个PCA通道,且删去另外两路通道同时去掉图2中的连接,最终模型图如图3所示。
4)、实施例1所述方法;
对四种方法使用测试集的测试结果分别统计F1指标、处理时间秒/万条,最终结果如表5:
表5:对比测试结果
根据结果不难看出,本发明在准确率方面是最高的。方法1通过PCA的方法来进行预测,虽然时间很短,但是准确率偏低。本方法的优越性为在构建三通道网络的过程中,对比传统网络降低了数据维度,减少了计算量,减少了运算时间,对比计算量小的方法1和方法3,本发明有效提高了多种故障检测的准确率。
最后,还需要注意的是,以上列举的仅是本发明的若干个具体实施例。显然,本发明不限于以上实施例,还可以有许多变形。本领域的普通技术人员能从本发明公开的内容直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (9)
1.一种数据驱动的工业过程故障检测方法,其特征在于:包括的过程如下:
S01、获取工业生产的过程数据,过程数据的维度为53,总条数为m;然后构建一个m行53列的过程数据的原始矩阵Zsrc,将过程数据的原始矩阵Zsrc的每一列都减去训练好的列平均值矩阵Zmean,得到中心化的矩阵Zo;
S02、将中心化的矩阵Zo输入主成分分析法并利用训练好的特征向量矩阵A9计算特征数据Y1_online=A9ZO;将中心化的矩阵Zo输入独立成分分析法并利用训练好的矢量矩阵计算特征数据将中心化的矩阵Zo输入典型关联分析法并利用训练好的特征向量矩阵∑9和训练好的右奇异向量正交基V9计算特征数据Y3_online:
其中,pk为中心化的矩阵Zo的第k行;k∈(1,m)为过程数据的当前时刻的条数,将k从1到m依次带入式(23)后获得特征数据Y3_online;
S03、工业过程故障检测网络包括PCA通道、ICA通道和CCA通道的三通道卷积神经网络,将特征数据Y1_online作为PCA通道的输入,特征数据Y2_online作为ICA通道的输入,特征数据Y3_online作为CCA通道的输入,分别输入到工业过程故障检测网络中,推断得到过程数据的分类标签,分类标签为1-6属于推断结果异常,将数据打上异常标签存入数据库;如果分类标签为0属于推断结果正常,则将数据打上正常标签存入数据库。
2.根据权利要求1所述的一种数据驱动的工业过程故障检测方法,其特征在于:
步骤S02中所述主成分分析进行特征提取的过程为:
S0201、主成分分析法提取特征:
S020101、根据输入的中心化的矩阵求出协方差矩阵C:
其中,Z为输入的中心化的矩阵;
S020102、求出协方差矩阵C的特征值λi和对应的特征向量;
S020103、将协方差矩阵C的特征值λi从大到小排列作为主对角线的值,其它值全部为0,保存为特征值矩阵Λ;设定η>85%,根据式(3)来确定选取的特征值的个数k:
S020104:将特征向量按照对应的特征值λi大小,按从大到小的顺序从上到下按行排列成矩阵,保存为特征向量矩阵A,取特征向量矩阵A的前k行记为Ak;
S020105:所述主成分分析后的特征数据Y1:
Y1=AkZ (4)
经过主成分分析法的输出为特征数据Y1和特征向量矩阵Ak。
3.根据权利要求2所述的一种数据驱动的工业过程故障检测方法,其特征在于:
步骤S02中所述独立成分分析法进行特征提取的过程为:
S0202、独立成分分析法提取特征
S020201、获取白化矩阵W:
S020202:获取新的白化后的过程数据
Znew=WZ (6)
其中,Z为输入的中心化的矩阵;
S020203:令需要估计的分量个数为r,迭代次数为n,初始化n=0且选择随机的初始向量ω(0);
S020204:通过牛顿迭代法,按式(7)进行计算ω(n+1):
ω(n+1)=E{Znewg(ω(n)TZnew)}-E{g’(ωTZnew)}ω(n) (7)
其中,E表示均值,g代表非线性函数,a表示常数,g’(ωTZnew)表示函数g的导数,n为计算次数,ω(n)表示通过迭代法计算了n次的结果;
S020205、对得到的ω(n+1)按式(8)进行归一化,并按式(9)进行收敛判断:如果不满足式(9)的条件,那么n=n+1后,重复步骤S020204计算ω(n+1);
||ω(n+1)ω(n)T|-1|<ε (9)
其中,ε为误差;
S020206、所述独立成分分析提取的特征数据Y2为:
4.根据权利要求3所述的一种数据驱动的工业过程故障检测方法,其特征在于:
步骤S02中所述典型关联分析法进行特征提取的过程为:
S0203、典型关联分析法提取特征
S020301、构建历史向量pk与历史数据集Pk:
pk=[zk-1,zk-2,...,zk-h]T (12)
Pk=[pk-1,pk-2,...,pk+N-1]∈Rmh×N (13)
其中,h为历史向量长度,N为数据总条数,zk为输入的中心化的矩阵第k行数据,并记Pk的协方差矩阵为∑pp;
S020302:构建未来向量fk与未来数据集Fk:
fk=[zk,zk+1,...,zk+s]T (14)
Fk=[fk,fk+1,...,fk+N-1]∈Rsh×N (15)
其中,s为未来向量长度;
记未来数据集Fk的协方差矩阵为∑ff,记未来向量fk与未来数据集Fk的协方差为∑pf;
S020303、使用典型关联分析法使历史数据集Pk和未来数据集Fk构成的数据集相关性最大,对相关系数进行奇异值分解:
ρpf=U∑VT (16)
根据奇异值分解公式可得:
其中,∑ff为未来数据集Fk的协方差矩阵,∑pf是历史数据集Pk和未来数据集Fk的协方差,∑pp是历史数据集Pk的协方差矩阵,V为右奇异向量正交基,U为左奇异向量正交基,∑为特征向量矩阵;可得:
Σ=∑pf (19)
S020304、所述典型关联分析法得到的特征数据Y3为:
其中,y(k)表示k时刻经过典型关联分析法得到的特征数据,k∈(1,m)表示输入的中心化的矩阵的当前时刻的条数,m为输入的中心化的矩阵的的总条数,将k从1到m带入式(21),得出经过所述典型关联分析得到的特征数据Y3;
经过典型关联分析法的输出为特征数据Y3、特征向量矩阵∑和右奇异向量正交基V。
5.根据权利要求4所述的一种数据驱动的工业过程故障检测方法,其特征在于:
步骤S03中的所述工业过程故障检测网络包括三个等价通道的三通道卷积神经网络,分别为PCA通道、ICA通道和CCA通道;每个通道均包括四层结构:第一层为16个2x2的卷积核,步长为1,padding为1,激活函数为ReLU,输出为16x4x4;第二层为32个2x2的卷积核,步长为1,padding为0,激活函数也为ReLU,输出为32x3x3;第三层为16个2x2的卷积核,步长为1,padding为0,激活函数同为ReLU,输出为16x2x2;第四层为全连接层;
输入的特征数据经过PCA通道、ICA通道和CCA通道后均获得32维数据,将三通道的32维数据进行连接,并在通过一层全连接层,输出32维数据,最终通过输出层输出0到6的分类结果。
6.根据权利要求5所述的一种数据驱动的工业过程故障检测方法,其特征在于:
1)、田纳西伊斯曼的过程数据进行软件仿真获得样本数据并进行标签,正常数据的标签为0,异常数据的标签为1-6;将样本数据及其对应的标签按照8∶2划分为训练集和测试集;
2)、对训练集进行中心化的预处理后获得中心化的矩阵Ztrain,同时保存训练集每一列的平均值为所述训练好的列平均值矩阵Zmean=[zmean_1,zmean_2,...,zmean_53];
3)、中心化的矩阵Ztrain通过所述主成分分析法的获得特征向量矩阵A,取特征向量矩阵A的前9行并保存为所述训练好的特征向量矩阵A9,输出特征数据Y1=A9Ztrain;
5)、中心化的矩阵Ztrain通过所述经过典型关联分析法获得特征向量矩阵∑和右奇异向量正交基V,然后维度取n=9,保存特征向量矩阵∑的前9×9对角阵为所述训练好的特征向量矩阵∑9和右奇异向量正交基V的前9×9对角阵为所述训练好的右奇异向量正交基V9,将k从训练集的第一条到最后一条依次代入按式(23),输出特征数据Y3:
其中,pk为中心化的矩阵Ztrain的第k行。
7.根据权利要求6所述的一种数据驱动的工业过程故障检测方法,其特征在于:
所述测试集构建一个m行53列的测试集的原始矩阵Zsrc,将原始矩阵Zsrc的每一列都减去训练好的列平均值矩阵Zmean,得到中心化的矩阵Ztest,然后进行特征提取:
将中心化的矩阵Ztest通过主成分分析法,并利用训练好的特征向量矩阵A9计算特征数据Y1_test=A9Ztest;将中心化的矩阵Ztest通过独立成分分析法,并利用训练好的矢量矩阵计算特征数据将中心化的矩阵Ztest通过典型关联分析法,并利用训练好的特征向量矩阵∑9和训练好的右奇异向量正交基V9,将k从训练集的第一条到最后一条依次代入式(23)后输出特征数据Y3_test:
其中,pk为中心化的矩阵Ztest的第k行。
9.根据权利要求8所述的一种数据驱动的工业过程故障检测方法,其特征在于:
所述工业过程故障检测网络的训练过程为:所述特征数据Y1输入PCA通道、特征数据Y2输入ICA通道、特征数据Y3输入CCA通道进行模型训练,学习率设置为0.001和损失函数为交叉熵损;共训练100个epoch病保存每个epoch训练好的模型,共保存100个训练好的工业过程故障检测网络模型,每个模型输出为训练集每一条样本数据对应的分类标签,其中标签0为正常数据,标签1-6代表故障数据;
所述工业过程故障检测网络的测试过程为:将特征数据Y1_test作为PCA通道的输入,特征数据Y2_test作为ICA通道的输入,特征数据Y3_test作为CCA通道的输入,将测试集分别输入所述100个训练好的工业过程故障检测网络模型中测试,得出每个模型的测试的F1分数,取出分数最高的工业过程故障检测网络模型为在线使用的所述工业过程故障检测网络。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210367561.XA CN114647231A (zh) | 2022-04-08 | 2022-04-08 | 一种数据驱动的工业过程故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210367561.XA CN114647231A (zh) | 2022-04-08 | 2022-04-08 | 一种数据驱动的工业过程故障检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114647231A true CN114647231A (zh) | 2022-06-21 |
Family
ID=81997595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210367561.XA Pending CN114647231A (zh) | 2022-04-08 | 2022-04-08 | 一种数据驱动的工业过程故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114647231A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116994709A (zh) * | 2023-06-01 | 2023-11-03 | 广州营悦营养健康咨询有限公司 | 一种个性化的饮食与运动推荐方法、系统及电子设备 |
-
2022
- 2022-04-08 CN CN202210367561.XA patent/CN114647231A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116994709A (zh) * | 2023-06-01 | 2023-11-03 | 广州营悦营养健康咨询有限公司 | 一种个性化的饮食与运动推荐方法、系统及电子设备 |
CN116994709B (zh) * | 2023-06-01 | 2024-04-05 | 精准营养医学科技研究院(深圳)有限公司 | 一种个性化的饮食与运动推荐方法、系统及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110132598B (zh) | 旋转设备滚动轴承故障噪声诊断算法 | |
CN108803520B (zh) | 一种基于变量非线性自相关性剔除的动态过程监测方法 | |
CN106843195B (zh) | 基于自适应集成半监督费舍尔判别的故障分类方法 | |
CN103914064B (zh) | 基于多分类器和d-s证据融合的工业过程故障诊断方法 | |
CN106649789B (zh) | 一种基于集成半监督费舍尔判别的工业过程故障分类方法 | |
CN108875772B (zh) | 一种基于堆叠稀疏高斯伯努利受限玻尔兹曼机和强化学习的故障分类模型及方法 | |
CN111580506A (zh) | 基于信息融合的工业过程故障诊断方法 | |
CN104699077A (zh) | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 | |
CN108345284A (zh) | 一种基于两变量块的质量相关故障检测方法 | |
CN111340110A (zh) | 一种基于工业过程运行状态趋势分析的故障预警方法 | |
CN111122811A (zh) | 一种oica和rnn融合模型的污水处理过程故障监测方法 | |
CN114330549A (zh) | 一种基于深度图网络的化工过程故障诊断方法 | |
CN114647231A (zh) | 一种数据驱动的工业过程故障检测方法 | |
CN114757269A (zh) | 一种基于局部子空间-邻域保持嵌入的复杂过程精细化故障检测方法 | |
Yao et al. | Fault diagnosis based on RseNet-LSTM for industrial process | |
CN115935277A (zh) | 基于小样本学习的增压锅炉故障诊断方法及故障诊断模型的训练方法和测试方法 | |
CN116415505A (zh) | 一种基于sbr-dbn模型的系统故障诊断与状态预测方法 | |
CN112947649B (zh) | 一种基于互信息矩阵投影的多变量过程监控方法 | |
CN108572639B (zh) | 一种基于主成分自相关性剔除的动态过程监测方法 | |
CN108491878B (zh) | 一种基于多个误差生成模型的故障分类诊断方法 | |
GB2616996A (en) | Method for recognizing type of vortex signal of evaporator of nuclear power plant on basis of LSTM-CNN | |
CN113157561A (zh) | 一种数控系统软件模块的缺陷预测方法 | |
CN110347579B (zh) | 基于神经元输出行为模式的深度学习测试用例的选择方法 | |
CN115017978A (zh) | 一种基于加权概率神经网络的故障分类方法 | |
CN111897310B (zh) | 基于一维多头卷积网络的工业过程故障分类方法和系统 |
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 |