CN104699077B - 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 - Google Patents
一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 Download PDFInfo
- Publication number
- CN104699077B CN104699077B CN201510077082.4A CN201510077082A CN104699077B CN 104699077 B CN104699077 B CN 104699077B CN 201510077082 A CN201510077082 A CN 201510077082A CN 104699077 B CN104699077 B CN 104699077B
- Authority
- CN
- China
- Prior art keywords
- fault
- data
- matrix
- variable
- component
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 91
- 238000004458 analytical method Methods 0.000 title claims abstract description 34
- 238000005192 partition Methods 0.000 title abstract 2
- 238000003745 diagnosis Methods 0.000 claims abstract description 67
- 230000008569 process Effects 0.000 claims abstract description 48
- 239000011159 matrix material Substances 0.000 claims description 136
- 239000013598 vector Substances 0.000 claims description 61
- 238000009826 distribution Methods 0.000 claims description 33
- 238000005259 measurement Methods 0.000 claims description 29
- 239000006185 dispersion Substances 0.000 claims description 11
- 238000002955 isolation Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000013144 data compression Methods 0.000 claims description 7
- 230000009467 reduction Effects 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000012824 chemical production Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 8
- 230000008439 repair process Effects 0.000 abstract description 4
- 239000004615 ingredient Substances 0.000 description 15
- 230000008859 change Effects 0.000 description 10
- 239000000498 cooling water Substances 0.000 description 10
- 238000001311 chemical methods and process Methods 0.000 description 9
- 239000007789 gas Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003292 glue Substances 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 239000010985 leather Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G05B23/0254—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 based on a quantitative model, e.g. mathematical relationships between inputs and outputs; functions: observer, Kalman filter, residual calculation, Neural Networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法,该方法充分发掘了过程数据所包含的潜在信息,能有效区分对故障起不同程度作用的变量,从而根据不同类别分别建立判别模型,改善了模型精度,大大提高了在线故障诊断的性能,增强了实际在线故障诊断的可靠性和可信度,有助于工业工程师对故障进行准确的修复,从而保证实际生产的安全可靠运行和产品的高质量追求。
Description
技术领域
本发明属于化工过程统计监测领域,特别是涉及一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法。
背景技术
作为工业生产中一种重要的生产方式,化工过程与人们的生活息息相关,已被广泛应用于冶金、炼油、造纸、制革等领域。如何保证化工过程生产安全,提高产品质量和经济效益是人们关注的焦点。随着化工过程日益复杂,在线故障检测及诊断也越来越重要。故障诊断是指在检测到故障发生后,进一步判断发生了哪种故障。而通过故障变量隔离则可以很好的识别发生的故障类别,保障生产的安全可靠运行以及产品的高质量,从而可以避免重大安全事故,减少人员伤亡以及提高经济效益。随着技术的发展,工业现场可以获得越来越多的数据,基于数据的在线故障诊断策略越来越受到研究人员和现场工程师的青睐。化工过程的统计建模、在线监测、故障诊断及质量预测已成为广泛的研究课题。
前人对此已经作了相应的研究与探讨,基于不同的角度提出了相应的在线故障诊断办法。归纳起来有如下几种:基于重构指标的方法、基于相似度计算的模式匹配方法、基于模型的改进隔离方法以及降维的方法。其中,降维的方法如主元分析、费舍尔判别分析能有效处理高维度、高相关性的数据,它们通过构造潜变量将高维度原始测量数据投影到低维度的监测空间,从而提高了故障诊断精度,被广泛应用于工业过程在线故障诊断。总体来说,上述几种基于降维的故障诊断方法各有各的适用场合与优缺点。相对而言,费舍尔判别分析方法侧重于区分具有不同特性的数据,在故障诊断方面更有优势。但是,基于传统费舍尔判别分析的在线故障诊断方法应用于实际化工过程时存在三方面的问题:首先,化工过程数据往往是高度耦合的,这可能导致类内散布矩阵是奇异的,从而无法进行奇异值分解提取过程数据潜在信息。其次,由于类间散布矩阵奇异,有可能导致判别成分的个数小于类别的个数,从而使得散布矩阵无法提供充足的过程信息。最后,在每个类中,所提取的判别成分是线性相关的,这就导致所提取的过程信息冗余。针对传统费舍尔判别分析方法的不足,研究人员提出了一系列的改进方法。总的来说,这些方法均采用两步法解决高耦合数据带来的散步矩阵奇异性问题,其关键在于如何在进行费舍尔判别分析之前进行数据降维。然而,前人的方法在解决奇异性问题时均存在一定程度的问题,如数据压缩不当而导致无法提取出过程数据的关键潜在信息,或者过程重要信息缺失等等,从而导致故障诊断精度欠缺。为了克服上述问题,研究人员提出一种迭代嵌套的费舍尔判别分析方法,提高了故障诊断性能。然而,该方法没有考虑到不同的变量对故障的影响程度不同,有些变量对故障影响很大而另一些变量与故障不甚相关,使用所有变量数据建立诊断模型,降低了模型精度,影响了在线故障诊断的性能。
本发明的内容深入考虑了化工过程的复杂性、数据的高维度高耦合性以及不同变量对故障作用的程度不同,提出了一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法。该方法充分发掘了故障相对于正常数据的变化情况,区分了对故障有重要影响的故障变量和没有影响的普通变量,并分别建立了故障诊断模型,大大提高了在线故障诊断性能。到目前为止,尚未见到与本发明相关的研究报道。
发明内容
本发明的目的在于针对现有针对化工生产过程的在线故障诊断技术的不足,提供一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法。
本发明的目的是通过以下技术方案来实现的:一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法,该方法包括以下步骤:
(1)获取过程分析数据:设一个化工生产过程具有J个测量变量和操作变量,则每一次采样可得到一个1×J的向量,采样K次后得到的数据表述为一个二维矩阵X(K×J),所述测量变量为运行过程中可被测量的状态参数,包括流量、温度、速率;所述操作变量包括给料量、阀门开度;分别获取正常数据二维矩阵Xn(K×J)和故障数据二维矩阵Xf,m(K×J),其中,下标n表示正常数据,下标f表示故障数据,m表示故障的类别;将正常数据和故障数据统一标示为Xi(K×J),其中下标i表示数据的类别;
(2)选取正常数据样本和一类故障数据样本作为总样本其中,由Xi(i=1,2)从上到下排列组成;
(3)数据准备:分别计算总样本均值向量每类样本均值向量总类内散布矩阵Sw和类间的散布矩阵Sb,计算公式如下:
其中,Si是每个类的散布矩阵;
(4)提取初始判别成分,该步骤由以下子步骤来实现:
(4.1)最大化类间离散度:求取使类间离散度最大的权重向量w,即相当于求取类间散布矩阵Sb的最大特征值所对应的特征向量w,所述类间离散度为wTSbw,获取w后,按公式(2)求取相应的总样本初始判别成分t;
其中,是减均值中心化后的总样本,那么对于每一类样本,其所对应的类判别成分为可知,t由ti从上到下依次排列构成;
(4.2)数据压缩:对减均值中心化后的总样本根据下式进行数据压缩:
其中:p表示总样本的负载向量,表示总样本中与t无关的残差;
同理,对于每类样本都可以通过公式(4)得到与ti无关的残差且由从上到下排列组成:
最后,用上述数据压缩关系wpT更新每一个类的信息,以保证判别成分的正交性:
Ei=Xi-XiwpT (5)
(4.3)迭代更新过程数据
(4.3.1)用步骤(4.2)中获得的Ei代替步骤(3)中的Xi,按步骤(3)重新计算总样本均值每类样本均值向量总类内散布矩阵Sw和类间的散布矩阵Sb,按步骤(4.1)、(4.2)再次提取初始判别成分;
(4.3.2)重复步骤(4.3.1)直到所提取的初始判别成分的个数等于Sw的阶数N;那么,同时可以得到由权重向量w组成的权重矩阵W(J×N)和相应的负载向量p组成的负载矩阵P(J×N)、总样本初始判别成分t组成的总样本的初始判别成分矩阵其中,T由Ti按从上至下排列构成,Ti是每个类的判别成分矩阵;最后,求取初始判别成分的系数矩阵R=W(PTW)-1,且T和Ti可直接由系数矩阵根据公式(6)求出:
(5)提取最终判别成分,该步骤通过以下子步骤来实现:
(5.1)过程数据预处理:使用XiR代替每类初始数据集合Xi,按步骤(3)重新计算每类样本均值总样本均值总类内散布矩阵Sw *以及类间散布矩阵Sb *;
(5.2)确定最终判别成分:最终判别成分通过以下步骤来确定:
(5.2.1)求取最优判别成分方向向量w*,使得类间散布矩阵与类内散布矩阵的比值J(θ)最大;其中,w*则可通过公式(7)求取矩阵Sw *-1Sb *最大特征值所对应的特征向量得到:
Sw *-1Sb *w*=λw* (7)
(5.2.2)求取每类的最终判别成分向量ti *:
(5.2.3)将ti *从上至下依次排列构成总样本的最终判别成分向量
(5.3)压缩过程数据:为了保证每类样本的判别成分之间是正交的,进行如下处理:
其中,是每类的负载向量,Ei *是与ti *无关的残差;
(6)迭代更新过程数据,该步骤包括以下子步骤:
(6.1)用步骤(5.3)中Ei *代替步骤(3)中的Xi,按步骤(3)重新计算每类样本均值、总样本均值,总类内散布矩阵Sw *以及类间散布矩阵Sb *,按步骤(4)和步骤(5)再次提取最终判别成分向量ti *;
(6.2)重复步骤(6.1)直至获得足够的最终判别成分ti *并构成最终的判别成分矩阵Ti *,Ti *所保留的最终判别成分个数为R,所述R通过交叉检验的方法确定;相应的,同时可以获得权重矩阵Θ(J×R)和负载矩阵Pi *(J×R);其中,Θ(J×R)和Pi *(J×R)分别由θ(J×1)和pi *(J×1)构成;
(6.3)求取最终系数矩阵Ri *(J×R):
Ri *=Θ(Pi *TΘ)-1 (10)
至此,步骤(2)中选取的该类故障最终系数矩阵被求取出来;
(7)选取正常数据和另一类故障数据作为总样本,重复步骤(4)-(6),获得该类故障样本的最终系数矩阵;
(8)重复步骤(7)直到M类故障的最终系数矩阵都被求取出来;
(9)故障数据衡量;该步骤由以下子步骤完成:
(9.1)求取每类故障Xf,m及相应正常数据Xn的类内成分矩阵Tf,m *和Tn,m *:
(9.2)求取每类故障Xf,m及相应正常数据Xn的D2指标:
其中,xn,k是Xn中一个样本,tn,m,k是xn,k的判别成分,是公式(11)中Tn,m *的均值向量,Σn则表示Tn,m *的协方差矩阵,如果方差过小而趋近于0,那么它们将被置为1;同理,xf,m是Xf,m中一个样本,tn,m,k是xf,m的判别成分;
(9.2)建立正常数据基于D2指标的控制限:由于过程数据服从多变量正态分布,那么可知D2统计量服从带权重的χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(10)比较每类故障数据的D2指标与相应的正常数据控制限若没有超限报警则说明该类故障与正常数据相似,无需进行以下步骤;若有超限报警,则执行以下步骤辨识和分离故障变量;
(11)变量贡献度衡量,该步骤包括以下子步骤;
(11.1)按公式(13)计算每类故障及相应正常数据的变量贡献度指标
其中,下标·表示正常数据或者故障数据,x·,k表示相应正常或者故障数据中的一个样本,t·,m,k表示相应样本x·,k的判别成分,x·,k,j是相应数据中的第k个样本的第j个变量,是正常数据中第j个变量的均值,rm,j是第j个变量的权重系数;
(11.2)因为正常数据服从一定的分布规律,根据分布规律确定正常数据贡献度指标的控制限;
(11.3)按公式(14)计算故障样本相对于正常样本的变量贡献度的比值
其中,表示正常数据的控制限;
(11.4)对于每个变量,将通过公式(14)所计算得到的M类故障对于正常样本的变量贡献度的比值构成向量并计算的均值
(12)故障变量选择;将步骤(11.4)所得到的排序,选择最大所对应的变量作为故障变量,将该变量从故障数据和正常数据中移除,得到新的故障数据和正常数据
(13)模型更新;对更新后的故障数据和正常数据执行步骤(3)~(8),求取新的最终系数矩阵,即代表故障方向;
(14)基于更新后的最终系数矩阵,重复执行步骤(9)~(12)直到所有故障变量都被选取出来;对于每类故障数据,故障变量构成矩阵其中,Jf,m表示故障变量的个数,Nf表示故障样本的个数;则每类故障被分为两部分数据:对故障有重大影响的以及对故障没有影响的其中,Jn,m=J-Jf,m;同样,对于正常数据也被分为两部分和其中,Nn表示正常数据的个数,Jn,m表示普通变量的个数;那么,对于每一对正常数据Xn和故障数据Xf,m,可以得到由故障变量构成的矩阵和以及由剩余的普通变量构成的矩阵和
(15)建立故障诊断模型,该步骤包括以下子步骤:
(15.1)用和分别代替步骤(2)中获取的正常样本数据和故障数据,重复步骤(3)至(8)提取出新的最终系数矩阵即表示有重要区分度的故障方向;
(15.2)计算每类故障子数据的判别成分和故障衡量统计指标
其中,表示判别成分矩阵,是的行向量的转置,是的均值向量,是的协方差矩阵;
(15.3)建立每类故障变量衡量指标的控制限;由于过程数据服从多维高斯分布,且样本数量足够大,故变量衡量指标服从加权χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(16)将步骤(14)中普通变量构成的矩阵和分别进行减均值除标准差的预处理后,构成大样本矩阵对进行PCA分解建模,得到负载矩阵表示所监测的波动方向;
(17)计算普通变量矩阵的故障衡量指标
其中,其中,表示判别成分矩阵,是的行向量的转置,是的均值向量,是的协方差矩阵;
(18)建立每类普通变量衡量指标的控制限;由于过程数据服从多维高斯分布,且样本数量足够大,故变量衡量指标服从加权χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(19)在线故障诊断;当发生故障后,需要检测发生故障的变量以及故障的类型,该步骤包括以下子步骤:
(19.1)按照步骤(1)获取新数据xnew(J×1),xnew(J×1)可以被分为两部分,故障变量部分和普通变量部分即
(19.2)调用故障诊断模型衡量故障变量:
其中,是新数据中故障变量部分,是的判别成分,和均由步骤(15.2)求得,为故障变量衡量指标;
(19.3)调用故障诊断模型衡量普通变量:
其中,是新数据中普通变量部分,是的判别成分,和均由步骤(17)求得,是普通变量衡量指标;
(19.4)依次比较和其控制限 和其控制限若均没有超限报警,则表示所调用的故障诊断模型能很好地诊断出故障,则判断该故障属于调用的故障诊断模型所表示的类别;若有一个以上指标超限,则继续调用该故障诊断模型对接下来的新故障数据进行诊断,按步骤(19.1)~(19.3)计算和若有连续新故障数据的和有一个以上指标超限,则表示故障不属于该故障诊断模型所表示的故障类别;考虑到误报率和漏报率,需要对连续多个新故障数据进行统计分析来判断故障的类别,记录新故障数据在调用每类故障诊断模型时两个D2指标的未超限个数,则新故障数据属于D2指标的未超限个数最大的故障类别。
本发明的有益效果是:该方法克服了传统方法存在的类内散布矩阵奇异性问题、判别成分的个数局限问题、判别成分线性相关问题,使得方法可以应用于高维度高耦合数据的复杂化工过程。该方法将传统的判别分析方法与变量选择方法有效结合,区分了故障变量和普通变量,建立了更准确的诊断模型,能有效地区分生产过程中故障类别,提高了在线故障诊断的性能,有助于工程师准确修复故障,确保生产的安全可靠运行以及产品的高质量追求。
附图说明
图1是本发明基于嵌套迭代费舍尔判别分析的故障变量隔离方法的流程图;
图2是本发明方法调用故障1的诊断模型针对故障2的在线故障诊断结果图,(a)为普通变量,(b)为故障变量;
图3是本发明方法调用故障2的诊断模型针对故障2的在线故障诊断结果图,(a)为普通变量,(b)为故障变量;
图4是本发明方法调用故障9的诊断模型针对故障2的在线故障诊断结果图(a)为普通变量,(b)为故障变量;
图5是本发明方法调用故障15的诊断模型针对故障2的在线故障诊断结果图,(a)为普通变量,(b)为故障变量。
具体实施方式
下面结合附图及具体实例,对本发明作进一步详细说明。
以田纳西-伊斯曼过程为例,田纳西-伊斯曼过程是一个典型的复杂化工生产过程,该过程由五个主要运行机构构成,分别是反应器、产品冷凝器、蒸汽液体分离器、循环压缩机以及产品剥离器。在整个过程中可以采集到两部分过程变量:41个测量变量以及11个操作变量。变量见表1和表2。
表1田纳西-伊斯曼过程测量变量表
序号 | 变量名称 | 序号 | 变量名称 | 序号 | 变量名称 |
1 | A组分进料流量 | 2 | D组分进料流量 | 3 | E组分进料流量 |
4 | A和C组分进料流量 | 5 | 循环流量 | 6 | 反应器进料流量 |
7 | 反应器压力 | 8 | 反应器液位 | 9 | 反应器温度 |
10 | 放空流量 | 11 | 产品分离温度 | 12 | 产品分离器液位 |
13 | 产品分离器压力 | 14 | 产品分离器下部出料 | 15 | 汽提塔液位 |
16 | 汽提塔压力 | 17 | 汽提塔下部出料 | 18 | 汽提塔温度 |
19 | 反应器冷却水出口温度 | 20 | 压缩机功率 | 21 | 汽提塔蒸汽流量 |
22 | 汽提塔冷却水出口温度 | 23 | 成分A(反应器进料) | 24 | 成分B(反应器进料) |
25 | 成分C(反应器进料) | 26 | 成分D(反应器进料) | 27 | 成分E(反应器进料) |
28 | 成分F(反应器进料) | 29 | 成分A(放空气体分析) | 30 | 成分B(放空气体分析) |
31 | 成分C(放空气体分析) | 32 | 成分D(放空气体分析) | 33 | 成分E(放空气体分析) |
34 | 成分F(放空气体分析) | 35 | 成分G(放空气体分析) | 36 | 成分H(放空气体分析) |
37 | 成分D(产品分析) | 38 | 成分E(产品分析) | 39 | 成分F(产品分析) |
40 | 成分G(产品分析) | 41 | 成分H(产品分析) |
表2田纳西-伊斯曼过程操作变量表
序号 | 变量名称 | 序号 | 变量名称 | 序号 | 变量名称 |
1 | A组分进料 | 2 | D组分进料 | 3 | E组分进料 |
4 | A和C组分进料 | 5 | 循环阀门 | 6 | 放空阀门 |
7 | 分离器液相流量 | 8 | 汽提塔液相流量 | 9 | 汽提塔蒸汽阀门 |
10 | 反应器冷却水流量 | 11 | 冷凝器冷却水流量 |
如图1所示,本发明一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法,包括以下步骤:
(1)获取过程分析数据:设一个化工生产过程具有J个测量变量和操作变量,则每一次采样可得到一个1×J的向量,采样K次后得到的数据表述为一个二维矩阵X(K×J),所述测量变量为运行过程中可被测量的状态参数,包括流量、温度、速率;所述操作变量包括给料量、阀门开度;分别获取正常数据二维矩阵Xn(K×J)和故障数据二维矩阵Xf,m(K×J),其中,下标n表示正常数据,下标f表示故障数据,m表示故障的类别;将正常数据和故障数据统一标示为Xi(K×J),其中下标i表示数据的类别;
表3田纳西-伊斯曼过程故障表
序号 | 故障变量 | 发生类型 | 序号 | 变量名称 | 发生类型 |
1 | APC进料变化 | 跃变 | 9 | 进料2温度变化 | 随机 |
2 | 成分B变化 | 跃变 | 10 | 进料4温度变化 | 随机 |
3 | 进料2温度变化 | 跃变 | 11 | 反应器冷却水温度变化 | 随机 |
4 | 反应器冷却水速度变化 | 跃变 | 12 | 冷凝器冷却水温度变化 | 随机 |
5 | 冷凝器冷却水速度变化 | 跃变 | 13 | 反应动力学特性变化 | 缓慢漂移 |
6 | 进料1损失 | 跃变 | 14 | 反应器冷却水阀门 | 粘滞 |
7 | 成分C进料压力下降 | 跃变 | 15 | 冷凝器冷却水阀门 | 粘滞 |
8 | 进料4中ABC组分变化 | 随机 |
(2)选取正常数据样本和一类故障数据样本作为总样本其中,由Xi(i=1,2)从上到下排列组成;
(3)数据准备:分别计算总样本均值向量每类样本均值向量总类内散布矩阵Sw和类间的散布矩阵Sb,计算公式如下:
其中,Si是每个类的散布矩阵;
(4)提取初始判别成分,该步骤由以下子步骤来实现:
(4.1)最大化类间离散度:求取使类间离散度最大的权重向量w,即相当于求取类间散布矩阵Sb的最大特征值所对应的特征向量w,所述类间离散度为wTSbw,获取w后,按公式(2)求取相应的总样本初始判别成分t;
其中,是减均值中心化后的总样本,那么对于每一类样本,其所对应的类判别成分为可知,t由ti从上到下依次排列构成;
(4.2)数据压缩:对减均值中心化后的总样本根据下式进行数据压缩:
其中:p表示总样本的负载向量,表示总样本中与t无关的残差;
同理,对于每类样本都可以通过公式(4)得到与ti无关的残差且由从上到下排列组成:
最后,用上述数据压缩关系wpT更新每一个类的信息,以保证判别成分的正交性:
Ei=Xi-XiwpT (5)
(4.3)迭代更新过程数据
(4.3.1)用步骤(4.2)中获得的Ei代替步骤(3)中的Xi,按步骤(3)重新计算总样本均值每类样本均值向量总类内散布矩阵Sw和类间的散布矩阵Sb,按步骤(4.1)、(4.2)再次提取初始判别成分;
(4.3.2)重复步骤(4.3.1)直到所提取的初始判别成分的个数等于Sw的阶数N;那么,同时可以得到由权重向量w组成的权重矩阵W(J×N)和相应的负载向量p组成的负载矩阵P(J×N)、总样本初始判别成分t组成的总样本的初始判别成分矩阵其中,T由Ti按从上至下排列构成,Ti是每个类的判别成分矩阵;最后,求取初始判别成分的系数矩阵R=W(PTW)-1,且T和Ti可直接由系数矩阵根据公式(6)求出:
(5)提取最终判别成分,该步骤通过以下子步骤来实现:
(5.1)过程数据预处理:使用XiR代替每类初始数据集合Xi,按步骤(3)重新计算每类样本均值总样本均值总类内散布矩阵Sw *以及类间散布矩阵Sb *;
(5.2)确定最终判别成分:最终判别成分通过以下步骤来确定:
(5.2.1)求取最优判别成分方向向量w*,使得类间散布矩阵与类内散布矩阵的比值J(θ)最大;其中,w*则可通过公式(7)求取矩阵Sw *-1Sb *最大特征值所对应的特征向量得到:
Sw *-1Sb *w*=λw* (7)
(5.2.2)求取每类的最终判别成分向量ti *:
(5.2.3)将ti *从上至下依次排列构成总样本的最终判别成分向量
(5.3)压缩过程数据:为了保证每类样本的判别成分之间是正交的,进行如下处理:
其中,是每类的负载向量,Ei *是与ti *无关的残差;
(6)迭代更新过程数据,该步骤包括以下子步骤:
(6.1)用步骤(5.3)中Ei *代替步骤(3)中的Xi,按步骤(3)重新计算每类样本均值、总样本均值,总类内散布矩阵Sw *以及类间散布矩阵Sb *,按步骤(4)和步骤(5)再次提取最终判别成分向量ti *;
(6.2)重复步骤(6.1)直至获得足够的最终判别成分ti *并构成最终的判别成分矩阵Ti *,Ti *所保留的最终判别成分个数为R,所述R通过交叉检验的方法确定;相应的,同时可以获得权重矩阵Θ(J×R)和负载矩阵Pi *(J×R);其中,Θ(J×R)和Pi *(J×R)分别由θ(J×1)和pi *(J×1)构成;
(6.3)求取最终系数矩阵Ri *(J×R):
Ri *=Θ(Pi *TΘ)-1 (10)
至此,步骤(2)中选取的该类故障最终系数矩阵被求取出来;
(7)选取正常数据和另一类故障数据作为总样本,重复步骤(4)-(6),获得该类故障样本的最终系数矩阵;
(8)重复步骤(7)直到M类故障的最终系数矩阵都被求取出来;
(9)故障数据衡量;该步骤由以下子步骤完成:
(9.1)求取每类故障Xf,m及相应正常数据Xn的类内成分矩阵Tf,m *和Tn,m *:
(9.2)求取每类故障Xf,m及相应正常数据Xn的D2指标:
其中,xn,k是Xn中一个样本,tn,m,k是xn,k的判别成分,是公式(11)中Tn,m *的均值向量,Σn则表示Tn,m *的协方差矩阵,如果方差过小而趋近于0,那么它们将被置为1;同理,xf,m是Xf,m中一个样本,tn,m,k是xf,m的判别成分;
(9.2)建立正常数据基于D2指标的控制限:由于过程数据服从多变量正态分布,那么可知D2统计量服从带权重的χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(10)比较每类故障数据的D2指标与相应的正常数据控制限若没有超限报警则说明该类故障与正常数据相似,无需进行以下步骤;若有超限报警,则执行以下步骤辨识和分离故障变量;
(11)变量贡献度衡量,该步骤包括以下子步骤;
(11.1)按公式(13)计算每类故障及相应正常数据的变量贡献度指标
其中,下标·表示正常数据或者故障数据,x·,k表示相应正常或者故障数据中的一个样本,t·,m,k表示相应样本x·,k的判别成分,x·,k,j是相应数据中的第k个样本的第j个变量,是正常数据中第j个变量的均值,rm,j是第j个变量的权重系数;
(11.2)因为正常数据服从一定的分布规律,根据分布规律确定正常数据贡献度指标的控制限;
(11.3)按公式(14)计算故障样本相对于正常样本的变量贡献度的比值
其中,表示正常数据的控制限;
(11.4)对于每个变量,将通过公式(14)所计算得到的M类故障对于正常样本的变量贡献度的比值构成向量并计算的均值
(12)故障变量选择;将步骤(11.4)所得到的排序,选择最大所对应的变量作为故障变量,将该变量从故障数据和正常数据中移除,得到新的故障数据和正常数据
(13)模型更新;对更新后的故障数据和正常数据执行步骤(3)~(8),求取新的最终系数矩阵,即代表故障方向;
(14)基于更新后的最终系数矩阵,重复执行步骤(9)~(12)直到所有故障变量都被选取出来;对于每类故障数据,故障变量构成矩阵其中,Jf,m表示故障变量的个数,Nf表示故障样本的个数;则每类故障被分为两部分数据:对故障有重大影响的以及对故障没有影响的其中,Jn,m=J-Jf,m;同样,对于正常数据也被分为两部分和其中,Nn表示正常数据的个数,Jn,m表示普通变量的个数;那么,对于每一对正常数据Xn和故障数据Xf,m,可以得到由故障变量构成的矩阵和以及由剩余的普通变量构成的矩阵和
(15)建立故障诊断模型,该步骤包括以下子步骤:
(15.1)用和分别代替步骤(2)中获取的正常样本数据和故障数据,重复步骤(3)至(8)提取出新的最终系数矩阵即表示有重要区分度的故障方向;
(15.2)计算每类故障子数据的判别成分和故障衡量统计指标
其中,表示判别成分矩阵,是的行向量的转置,是的均值向量,是的协方差矩阵;
(15.3)建立每类故障变量衡量指标的控制限;由于过程数据服从多维高斯分布,且样本数量足够大,故变量衡量指标服从加权χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(16)将步骤(14)中普通变量构成的矩阵和分别进行减均值除标准差的预处理后,构成大样本矩阵对进行PCA分解建模,得到负载矩阵表示所监测的波动方向;
(17)计算普通变量矩阵的故障衡量指标
其中,其中,表示判别成分矩阵,是的行向量的转置,是的均值向量,是的协方差矩阵;
(18)建立每类普通变量衡量指标的控制限;由于过程数据服从多维高斯分布,且样本数量足够大,故变量衡量指标服从加权χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(19)在线故障诊断;当发生故障后,需要检测发生故障的变量以及故障的类型,该步骤包括以下子步骤:
(19.1)按照步骤(1)获取新数据xnew(J×1),xnew(J×1)可以被分为两部分,故障变量部分和普通变量部分即
(19.2)调用故障诊断模型衡量故障变量:
其中,是新数据中故障变量部分,是的判别成分,和均由步骤(15.2)求得,为故障变量衡量指标;
(19.3)调用故障诊断模型衡量普通变量:
其中,是新数据中普通变量部分,是的判别成分,和均由步骤(17)求得,是普通变量衡量指标;
(19.4)依次比较和其控制限 和其控制限若均没有超限报警,则表示所调用的故障诊断模型能很好地诊断出故障,则判断该故障属于调用的故障诊断模型所表示的类别;若有一个以上指标超限,则继续调用该故障诊断模型对接下来的新故障数据进行诊断,按步骤(19.1)~(19.3)计算和若有连续新故障数据的和有一个以上指标超限,则表示故障不属于该故障诊断模型所表示的故障类别;考虑到误报率和漏报率,需要对连续多个新故障数据进行统计分析来判断故障的类别,记录新故障数据在调用每类故障诊断模型时两个D2指标的未超限个数,则新故障数据属于D2指标的未超限个数最大的故障类别。
本发明通过基于迭代嵌套费舍尔判别分析的故障变量隔离方法根据历史故障数据建立故障诊断模型,工程师可以实时获得新过程采样数据的在线故障诊断结果,判断故障的类别从而进行相应的修复措施。当基于历史故障数据建立的诊断模型能很好容纳新故障样本所求得的D2指标,即调用该类故障的诊断模型所得的新样本D2指标产生的报警信号显著少于调用其他类故障诊断模型所产生的报警信号,则说明新样本属于该类故障,那么当前故障被正确诊断出来。考虑到误报率和漏报率,统计了新样本D2在每类故障诊断模型下的未超限信号的个数,新样本则属于未超限个数最多的那类故障。为了描述故障诊断的性能,我们定义了一个指标—正确诊断指标,若故障被正确诊断,则该指标为1;否则,该指标值为0。表4总结了本发明所提出的方法与传统方法针对15种故障类型的正确诊断指标的比较结果,并且计算均值(Mean)以及标准差(STD)用于综合评估故障诊断的性能。
表4本方法与传统方法故障诊断性能对比(衡量指标:正确诊断指标)
从表中可以看出,基于本发明方法的故障诊断方法提高了实际在线故障诊断的可靠性和可信度,对于15种故障均能正确诊断出来,其性能优于传统方法。此外,图2至图5展示了本发明所提出方法针对于故障2的在线故障诊断结果。是实线表示每类故障基于D2指标的控制限,点线表示D2指标。可以看出,而本发明则能正确诊断出故障2,说明了本发明方法的有效性。总体来说,本发明的方法有优越的在线故障诊断性能,有助于工业工程师对故障进行准确修复,保证实际生产过程的安全可靠运行。
应该理解,本发明并不局限于上述具体实施例的田纳西-伊斯曼过程,凡是熟悉本领域的技术人员在不违背本发明精神的前提下还可做出等同变形或替换,这些等同的变型或替换均包含在本申请权利要求所限定的范围内。
Claims (1)
1.一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法,其特征在于,该方法包括以下步骤:
(1)获取过程分析数据:设一个化工生产过程具有J个测量变量和操作变量,则每一次采样可得到一个1×J的向量,采样K次后得到的数据表述为一个二维矩阵X(K×J),所述测量变量为运行过程中可被测量的状态参数,包括流量、温度、速率;所述操作变量包括给料量、阀门开度;分别获取正常数据二维矩阵Xn(K×J)和故障数据二维矩阵Xf,m(K×J),其中,下标n表示正常数据,下标f表示故障数据,m表示故障的类别;将正常数据和故障数据统一标示为Xi(K×J),其中下标i表示数据的类别;
(2)选取正常数据样本和一类故障数据样本作为总样本其中,由Xi(i=1,2)从上到下排列组成;
(3)数据准备:分别计算总样本均值向量每类样本均值向量总类内散布矩阵Sw和类间的散布矩阵Sb,计算公式如下:
其中,Si是每个类的散布矩阵;
(4)提取初始判别成分,该步骤由以下子步骤来实现:
(4.1)最大化类间离散度:求取使类间离散度最大的权重向量w,即相当于求取类间散布矩阵Sb的最大特征值所对应的特征向量w,所述类间离散度为wTSbw,获取w后,按公式(2)求取相应的总样本初始判别成分t;
其中,是减均值中心化后的总样本,那么对于每一类样本,其所对应的类判别成分为可知,t由ti从上到下依次排列构成;
(4.2)数据压缩:对减均值中心化后的总样本根据下式进行数据压缩:
其中:p表示总样本的负载向量,表示总样本中与t无关的残差;
同理,对于每类样本都可以通过公式(4)得到与ti无关的残差且由从上到下排列组成:
最后,用上述数据压缩关系wpT更新每一个类的信息,以保证判别成分的正交性:
Ei=Xi-XiwpT (5)
(4.3)迭代更新过程数据
(4.3.1)用步骤(4.2)中获得的Ei代替步骤(3)中的Xi,按步骤(3)重新计算总样本均值每类样本均值向量总类内散布矩阵Sw和类间的散布矩阵Sb,按步骤(4.1)、(4.2)再次提取初始判别成分;
(4.3.2)重复步骤(4.3.1)直到所提取的初始判别成分的个数等于Sw的阶数N;那么,同时可以得到由权重向量w组成的权重矩阵W(J×N)和相应的负载向量p组成的负载矩阵P(J×N)、总样本初始判别成分t组成的总样本的初始判别成 分矩阵其中,T由Ti按从上至下排列构成,Ti是每个类的判别成分矩阵;最后,求取初始判别成分的系数矩阵R=W(PTW)-1,且T和Ti可直接由系数矩阵根据公式(6)求出:
(5)提取最终判别成分,该步骤通过以下子步骤来实现:
(5.1)过程数据预处理:使用XiR代替每类初始数据集合Xi,按步骤(3)重新计算每类样本均值总样本均值总类内散布矩阵Sw *以及类间散布矩阵Sb *;
(5.2)确定最终判别成分:最终判别成分通过以下步骤来确定:
(5.2.1)求取最优判别成分方向向量w*,使得类间散布矩阵与类内散布矩阵的比值J(θ)最大;其中,w*则可通过公式(7)求取矩阵Sw *-1Sb *最大特征值所对应的特征向量得到:
Sw *-1Sb *w*=λw* (7)
(5.2.2)求取每类的最终判别成分向量ti *:
ti *=XiRw*=Xiθ (8)
θ=Rw*
(5.2.3)将ti *从上至下依次排列构成总样本的最终判别成分向量
(5.3)压缩过程数据:为了保证每类样本的判别成分之间是正交的,进行如下处理:
pi *T=(ti *Tti *)-1ti *TXi (9)
Ei *=Xi-ti *pi *T
其中,是每类的负载向量,Ei *是与ti *无关的残差;
(6)迭代更新过程数据,该步骤包括以下子步骤:
(6.1)用步骤(5.3)中Ei *代替步骤(3)中的Xi,按步骤(3)重新计算每类样本均值、总样本均值,总类内散布矩阵Sw *以及类间散布矩阵Sb *,按步骤(4)和步骤(5)再次提取最终判别成分向量ti *;
(6.2)重复步骤(6.1)直至获得足够的最终判别成分ti *并构成最终的判别成分矩阵Ti *,Ti *所保留的最终判别成分个数为R,所述R通过交叉检验的方法确定;相应的,同时可以获得权重矩阵Θ(J×R)和负载矩阵Pi *(J×R);其中,Θ(J×R)和Pi *(J×R)分别由θ(J×1)和pi *(J×1)构成;
(6.3)求取最终系数矩阵Ri *(J×R):
Ri *=Θ(Pi *TΘ)-1 (10)
至此,步骤(2)中选取的该类故障最终系数矩阵被求取出来;
(7)选取正常数据和另一类故障数据作为总样本,重复步骤(4)-(6),获得该类故障样本的最终系数矩阵;
(8)重复步骤(7)直到M类故障的最终系数矩阵都被求取出来;
(9)故障数据衡量;该步骤由以下子步骤完成:
(9.1)求取每类故障Xf,m及相应正常数据Xn的类内成分矩阵Tf,m *和Tn,m *:
Tf,m *=Xf,mRf,m * (11)
Tn,m *=XnRf,m *
(9.2)求取每类故障Xf,m及相应正常数据Xn的D2指标:
tn,m,k=xn,k TRf,m *
tf,m,k=xf,k TRf,m *
其中,xn,k是Xn中一个样本,tn,m,k是xn,k的判别成分,是公式(11)中Tn,m *的均值向量,Σn则表示Tn,m *的协方差矩阵,如果方差过小而趋近于0,那么它们将被置为1;同理,xf,m是Xf,m中一个样本,tn,m,k是xf,m的判别成分;
(9.3)建立正常数据基于D2指标的控制限:由于过程数据服从多变量正态分布,那么可知D2统计量服从带权重的χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(10)比较每类故障数据的D2指标与相应的正常数据控制限若没有超限报警则说明该类故障与正常数据相似,无需进行以下步骤;若有超限报警,则执行以下步骤辨识和分离故障变量;
(11)变量贡献度衡量,该步骤包括以下子步骤;
(11.1)按公式(13)计算每类故障及相应正常数据的变量贡献度指标
t·,m,k=x·,k TRf,m * (13)
其中,下标·表示正常数据或者故障数据,x·,k表示相应正常或者故障数据中的一个样本,t·,m,k表示相应样本x·,k的判别成分,x·,k,j是相应数据中的第k个样本的第j个变量,是正常数据中第j个变量的均值,rm,j是第j个变量的权重系数;
(11.2)因为正常数据服从一定的分布规律,根据分布规律确定正常数据贡献度指标的控制限;
(11.3)按公式(14)计算故障样本相对于正常样本的变量贡献度的比值
其中,表示正常数据的控制限;
(11.4)对于每个变量,将通过公式(14)所计算得到的M类故障对于正常样本的变量贡献度的比值构成向量并计算的均值
(12)故障变量选择;将步骤(11.4)所得到的排序,选择最大所对应的变量作为故障变量,将该变量从故障数据和正常数据中移除,得到新的故障数据和正常数据
(13)模型更新;对更新后的故障数据和正常数据执行步骤(3)~(8),求取新的最终系数矩阵,即代表故障方向;
(14)基于更新后的最终系数矩阵,重复执行步骤(9)~(12)直到所有故障变量都被选取出来;对于每类故障数据,故障变量构成矩阵其中,Jf,m表示故障变量的个数,Nf表示故障样本的个数;则每类故障被分为两部分数据:对故障有重大影响的以及对故障没有影响的其中,Jn,m=J-Jf,m;同样,对于正常数据也被分为两部分和 其中,Nn表示正常数据的个数,Jn,m表示普通变量的个数;那么,对于每一对正常数据Xn和故障数据Xf,m,可以得到由故障变量构成的矩阵 和以及由剩余的普通变量构成的矩阵和
(15)建立故障诊断模型,该步骤包括以下子步骤:
(15.1)用和分别代替步骤(2)中获取的正常样 本数据和故障数据,重复步骤(3)至(8)提取出新的最终系数矩阵即表示有重要区分度的故障方向;
(15.2)计算每类故障子数据的判别成分和故障衡量统计指标
其中,表示判别成分矩阵,是的行向量的转置,是的均值向量,是的协方差矩阵;
(15.3)建立每类故障变量衡量指标的控制限;由于过程数据服从多维高斯分布,且样本数量足够大,故变量衡量指标服从加权χ2分布,则可依据带权重的χ2分布建立D2统计量的控制限
(16)将步骤(14)中普通变量构成的矩阵和分别进行减均值除标准差的预处理后,构成大样本矩阵对进行PCA分解建模,得到负载矩阵表示所监测的波动方向;
(17)计算普通变量矩阵的故障衡量指标
其中,其中,表示判别成分矩阵,是的行向量的转置,是的均值向量,是的协方差矩阵;
(18)建立每类普通变量衡量指标的控制限;由于过程数据服从多维高斯分布,且样本数量足够大,故变量衡量指标服从加权χ2分布,则可依据带权重的χ2 分布建立D2统计量的控制限
(19)在线故障诊断;当发生故障后,需要检测发生故障的变量以及故障的类型,该步骤包括以下子步骤:
(19.1)按照步骤(1)获取新数据xnew(J×1),xnew(J×1)可以被分为两部分,故障变量部分和普通变量部分即
(19.2)调用故障诊断模型衡量故障变量:
其中,是新数据中故障变量部分,是的判别成分,和均由步骤(15.2)求得,为故障变量衡量指标;
(19.3)调用故障诊断模型衡量普通变量:
其中,是新数据中普通变量部分,是的判别成分,和均由步骤(17)求得,是普通变量衡量指标;
(19.4)依次比较和其控制限 和其控制限若均没有超限报警,则表示所调用的故障诊断模型能很好地诊断出故障,则判断该故障属于调用的故障诊断模型所表示的类别;若有一个以上指标超限,则继续调用该故障诊断模型对接下来的新故障数据进行诊断,按步骤(19.1)~(19.3)计算和若有连续新故障数据的和有一个以上指标超限,则表示故障不属于该故障诊断模型所表示的故障类别;考虑到误报率和漏报率,需要对连续多个新故障数据进行统计分析来判断故障的类别,记录新故障数据在调用每类 故障诊断模型时两个D2指标的未超限个数,则新故障数据属于D2指标的未超限个数最大的故障类别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510077082.4A CN104699077B (zh) | 2015-02-12 | 2015-02-12 | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510077082.4A CN104699077B (zh) | 2015-02-12 | 2015-02-12 | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104699077A CN104699077A (zh) | 2015-06-10 |
CN104699077B true CN104699077B (zh) | 2017-06-06 |
Family
ID=53346306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510077082.4A Active CN104699077B (zh) | 2015-02-12 | 2015-02-12 | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104699077B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105259890A (zh) * | 2015-08-18 | 2016-01-20 | 浙江中烟工业有限责任公司 | 超高速包装设备的统计监测与故障诊断方法 |
CN105573290B (zh) * | 2015-12-16 | 2017-12-29 | 浙江中烟工业有限责任公司 | 卷烟工厂超高速条盒包装机多工况过程在线监测与故障诊断方法 |
CN105487524B (zh) * | 2015-12-29 | 2017-12-08 | 浙江中烟工业有限责任公司 | 具有多工况特性的超高速小盒包装机状态监测与诊断方法 |
CN107220475A (zh) * | 2016-11-01 | 2017-09-29 | 重庆交通大学 | 一种基于线性判别分析的轴承特征数据分析方法 |
CN107103125A (zh) * | 2017-04-01 | 2017-08-29 | 宁波大学 | 一种基于两分类Fisher判别分析的故障诊断方法 |
CN106897542A (zh) * | 2017-04-18 | 2017-06-27 | 浙江中烟工业有限责任公司 | 基于显著故障变量提取的卷烟制叶丝段故障诊断方法 |
CN107153409B (zh) * | 2017-06-02 | 2019-08-16 | 宁波大学 | 一种基于缺失变量建模思路的非高斯过程监测方法 |
CN108492000B (zh) * | 2018-02-11 | 2021-07-02 | 浙江大学 | 面向百万千瓦超超临界机组非平稳特性的故障诊断方法 |
CN109407640B (zh) * | 2018-12-13 | 2021-03-09 | 宁波大学 | 一种基于动态正交成分分析的动态过程监测方法 |
CN111443686B (zh) * | 2020-03-23 | 2021-08-10 | 杭州电子科技大学 | 基于多目标优化与证据迭代更新的工业报警器设计方法 |
CN113079043B (zh) * | 2021-03-25 | 2021-12-10 | 浙江大学 | 一种基于码本的网络故障诊断方法 |
CN113484646A (zh) * | 2021-07-21 | 2021-10-08 | 国网河北省电力有限公司检修分公司 | 一种变电站二次回路复杂故障定位装置及诊断方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1996191A (zh) * | 2006-11-23 | 2007-07-11 | 浙江大学 | 一种基于fisher的工业过程非线性故障诊断系统及方法 |
WO2012040916A1 (zh) * | 2010-09-29 | 2012-04-05 | 东北大学 | 基于递归核主元分析的连续退火过程故障监测方法 |
CN103412557A (zh) * | 2013-07-17 | 2013-11-27 | 浙江工业大学 | 一种适于非线性过程在线监控的工业故障检测与诊断方法 |
CN103471865A (zh) * | 2013-09-12 | 2013-12-25 | 北京交通大学 | 基于线性判别法的列车悬挂系统故障分离方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7729789B2 (en) * | 2004-05-04 | 2010-06-01 | Fisher-Rosemount Systems, Inc. | Process plant monitoring based on multivariate statistical analysis and on-line process simulation |
US7349746B2 (en) * | 2004-09-10 | 2008-03-25 | Exxonmobil Research And Engineering Company | System and method for abnormal event detection in the operation of continuous industrial processes |
CN101169623B (zh) * | 2007-11-22 | 2010-07-07 | 东北大学 | 基于核主元分析贡献图的非线性过程故障辨识方法 |
CN101477372A (zh) * | 2009-01-07 | 2009-07-08 | 沈阳化工学院 | 化工生产非正常子域故障分离的技术方法 |
CN103777627B (zh) * | 2014-01-24 | 2016-03-30 | 浙江大学 | 一种基于少量批次的批次过程在线监测方法 |
-
2015
- 2015-02-12 CN CN201510077082.4A patent/CN104699077B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1996191A (zh) * | 2006-11-23 | 2007-07-11 | 浙江大学 | 一种基于fisher的工业过程非线性故障诊断系统及方法 |
WO2012040916A1 (zh) * | 2010-09-29 | 2012-04-05 | 东北大学 | 基于递归核主元分析的连续退火过程故障监测方法 |
CN103412557A (zh) * | 2013-07-17 | 2013-11-27 | 浙江工业大学 | 一种适于非线性过程在线监控的工业故障检测与诊断方法 |
CN103471865A (zh) * | 2013-09-12 | 2013-12-25 | 北京交通大学 | 基于线性判别法的列车悬挂系统故障分离方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104699077A (zh) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104699077B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 | |
CN104536439B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障诊断方法 | |
CN104714537B (zh) | 一种基于联合相对变化分析和自回归模型的故障预测方法 | |
CN107272667B (zh) | 一种基于平行的偏最小二乘法的工业过程故障检测方法 | |
CN109407652B (zh) | 基于主辅pca模型的多变量工业过程故障检测方法 | |
CN103926919B (zh) | 基于小波变换和Lasso函数的工业过程故障检测方法 | |
CN102736546A (zh) | 一种流程工业复杂机电系统的状态监测装置及方法 | |
CN110244692B (zh) | 化工过程微小故障检测方法 | |
CN105404280A (zh) | 基于自回归动态隐变量模型的工业过程故障检测方法 | |
CN105955219A (zh) | 基于互信息的分布式动态过程故障检测方法 | |
CN105259895A (zh) | 一种工业过程微小故障的检测和分离方法及其监测系统 | |
CN112904810B (zh) | 基于有效特征选择的流程工业非线性过程监测方法 | |
CN104914847A (zh) | 基于方向核偏最小二乘的工业过程故障诊断方法 | |
CN108388234A (zh) | 一种基于相关性划分多变量块pca模型的故障监测方法 | |
CN109298633A (zh) | 基于自适应分块非负矩阵分解的化工生产过程故障监测方法 | |
Zhong et al. | Multimode non‐Gaussian process monitoring based on local entropy independent component analysis | |
CN108830006B (zh) | 基于线性评价因子的线性-非线性工业过程故障检测方法 | |
CN106326915A (zh) | 一种基于改进核Fisher的化工过程故障诊断方法 | |
CN106354125A (zh) | 一种利用分块pca检测化工过程故障的方法 | |
CN110245460A (zh) | 一种基于多阶段oica的间歇过程故障监测方法 | |
CN103995985B (zh) | 基于Daubechies小波变换和弹性网的故障检测方法 | |
CN111983994B (zh) | 一种基于复杂工业化工过程的v-pca故障诊断方法 | |
CN114371677B (zh) | 基于谱半径-区间主成分分析的工业过程状态监测方法 | |
CN114879612A (zh) | 一种基于Local-DBKSSA的高炉炼铁过程监测方法 | |
CN114200914A (zh) | 一种基于mw-occa的质量相关早期故障检测方法 |
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 |