CN108536943A - 一种基于多生产单元变量交叉相关解耦策略的故障监测方法 - Google Patents
一种基于多生产单元变量交叉相关解耦策略的故障监测方法 Download PDFInfo
- Publication number
- CN108536943A CN108536943A CN201810280187.3A CN201810280187A CN108536943A CN 108536943 A CN108536943 A CN 108536943A CN 201810280187 A CN201810280187 A CN 201810280187A CN 108536943 A CN108536943 A CN 108536943A
- Authority
- CN
- China
- Prior art keywords
- matrix
- vector
- variable
- error
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Hardware Design (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于多生产单元变量交叉相关解耦策略的故障监测方法,旨在从数据角度出发,将工业对象中不同生产单元之间的交叉关系考虑进分布式建模与监测过程中,从而实施更加可靠而有效的分布式故障监测。具体来讲,首先根据各生产单元测量变量的归属,将所有测量变量划分成多个变量子块;其次,利用回归模型将各个变量子块中与其他变量子块之间的交叉相关信息挖掘出来;最后,利用交叉相关解耦后的误差实施建模与故障监测。相比于传统方法,本发明方法利用回归模型将不同生产单元变量子块之间的交叉关系考虑进来,并对能够反映出不同生产单元之间交叉相关关系是否发生变化的误差实施监测,理应具备更优越的故障监测性能。
Description
技术领域
本发明涉及一种数据驱动的故障监测方法,尤其涉及一种基于多生产单元变量交叉相关解耦策略的故障监测方法。
背景技术
保证持续正常的生产状态对于降低企业生产成本与保证生产安全具有重要性意义,通常采取的技术手段无外乎对过程运行状态实施实时监测,从而及时地甄别出***出现的异常状态。近年来,随着工业大规模与“大数据”建设的推进,生产过程中可以采集大量的实时数据而无法建立精确的机理模型,这为数据驱动的故障监测方法的兴起与广泛应用做好了铺垫。在这一研究领域里,多变量统计过程监测得到了最多的研究与关注。这其中,又当以主成分分析(Principal Component Analysis,PCA)算法为最主流的实施技术手段。通过挖掘训练数据的潜在特征,基于PCA的故障监测方法将故障监测任务转换成一个单分类的模式识别问题。从这点上看,以PCA为核心的故障监测方法体系的核心在于如何有效地描述采样数据间的关联关系。
此外,考虑到现代工业过程的大规模性,基于PCA的故障监测方法已经从最开始的单个故障监测模型发展成分布式的故障监测策略。一般来讲,分布式的故障监测方法能够在建模时简化对过程对象分析的难度,而且多模型的思路通常泛化性能优越于单个模型。由于现代工业过程的大规模趋势,整个生成单元通常由诸多个子生产单元组成。因此,可以将所有的测量变量按照生产单元的归属划分成多个子变量块,然后对每个子变量块分别实施监测就是分布式故障监测方法的核心本质。
然而,现代工业过程各生产子单元之间还存在交叉的关系,比如产品的回流进入其他生产单元,以及利用不同生产单元间的反馈信号实施控制等等。直接按照生产单元分成的变量子块的做法,不能将不同生产单元之间的交叉关系考虑进来。因此,传统的分布式故障监测方法还有待进一步的改进。若是按照生产单元实施分布式故障监测时,能将不同生产单元之间的交叉关联关系考虑进来,将会提升相应故障监测方法的性能。可是,若是直接根据生产过程的机理知识分析这种交叉关联关系,需要深入细致的分析过程对象的组成结构情况。这在实际过程中是不可取的,也违背了数据驱动故障监测方法从数据角度出发的理念。
发明内容
本发明所要解决的主要技术问题是:如何从数据角度出发,将工业对象中不同生产单元之间的交叉关系考虑进分布式建模与监测过程中,从而实施更加可靠而有效的分布式故障监测。具体来讲,本发明方法首先根据各生产单元测量变量的归属,将所有测量变量划分成多个变量子块;其次,利用回归模型将各个变量子块中与其他变量子块之间的交叉相关信息挖掘出来;最后,利用交叉相关解耦后的误差实施建模与故障监测。
本发明解决上述技术问题所采用的技术方案为:一种基于多生产单元变量交叉相关解耦策略的故障监测方法,包括以下步骤:
(1)采集生产过程对象中所有生产单元正常运行状态下的样本数据,组成训练数据集X∈Rn×m,并对其进行标准化处理,得到均值为0、标准差为1的新数据矩阵其中,n为训练样本总数,m为过程对象所有测量变量的个数,R为实数集,Rn×m表示n×m维的实数矩阵。
(2)根据各测量变量的生产单元归属,将m个测量变量分成C个不同的变量子块,其中C表示过程对象中生产单元的个数。
(3)依据C个不同的变量子块将矩阵对应地分成C个不同的子矩阵X1,X2,…,XC,其中c=1,2,…,C,mc为第c个变量子块中变量个数且满足条件∑mc=m。
(4)将第c个子矩阵Xc做为回归模型的输出,同时将Y=[X1,…,Xc-1,Xc+1,…XC]做为回归模型的输入,利用偏最小二乘算法(Partial Least Square,PLS)建立输入-输出之间的回归模型,具体的实施过程如下所示:
①置k=1与Z=Y后,初始化向量uk为矩阵Xc的第一列;
②根据公式wk=Zuk/||Zuk||、sk=Zwk、qk=Xc Tsk/(sk Tsk)分别计算得到系数向量wk、得分向量sk、和系数向量qk,其中||Zuk||表示计算向量Ziuk的长度,上标号T表示矩阵或者向量的转置;
③根据公式unew=Xcqk/qk 2计算向量unew;
④判断是否满足条件||uk-unew||<10-6?若否,则置uk=unew后返回步骤②;若是,则执行⑤;
⑤根据公式pk=ZTsk/(sk Tsk)计算得到第k个投影向量pk∈R(m-1)×1,并保留向量pk、向量wk、和系数向量qk;
⑥判断矩阵Yk=skpk T中的最大元素是否大于0.01?若是,根据公式Z=Z-skpk T更新矩阵Zi后执行步骤⑦;若否,则得到投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2,…,wk]、和向量Q=[q1,q2,…,qk],并执行步骤⑧;
⑦判断k<m-mc?若是,则置k=k+1后,返回步骤②;若否,则得到最终的投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2,…,wk]、和向量Q=[q1,q2,…,qk];
⑧根据公式br=Wr(Pr TWr)-1Qr T计算在保留r个投影向量前提下的回归系数向量br,其中r=1,2,…,k,Pr、Wr、和Qr分别为矩阵P、W、和Q中前1至r列向量组成的矩阵;
⑨根据公式计算输出Xc的预测值后,计算保留不同个数投影向量所对应的模型预测误差其中diag[]表示将矩阵对角线的元素单独组成向量的操作,sum{}表示求取向量中各元素之和;
⑩根据MSE1,MSE2,…,MSEk的数值变化情况找出预测误差不再发生显著变化的收敛点位置,并将其下标号对应的具体数值作为偏最小二乘模型需保留的投影变量个数K,其中K≤k;
根据K确定最终的输入Y与输出Xc之间的PLS模型,即:
Xc=YWK(PK TWK)-1QK T+Ec=YΘc+Ec (1)
上式中,Θc=Wk(PK TWk)-1QK T,为回归模型误差;
(5)重复步骤(4)直至得到C个PLS模型,并将模型误差E1,E2,…,EC合并成一个误差矩阵E=[E1,E2,…,EC]∈Rn×m。
(6)将估计误差矩阵E作为新的训练数据矩阵,对其中的每一列实施标准化处理,得到均值为0,标准差为1的新数据矩阵
(7)利用PCA算法为建立相应的PCA故障监测模型,并保留模型参数集Θ={B,Λ,Dlim,SPElim}以备调用,其中B为投影变量矩阵,Λ是由特征值组成的对角矩阵,Dlim与SPElim分别表示监测统计量D与SPE的控制上限,具体的实施过程如下所示:
①计算的协方差矩阵
②求解S所有特征值γ1≥γ2≥…≥γm所对应的特征向量β1,β2…,βm;
③设置保留的主成分个数d为满足如下所示条件的最小值,将对应的d个最大特征值组成对角矩阵Λ∈Rd×d,并将对应的d个特征向量组成投影变换矩阵B=[β1,β2…,βd]
④根据如下所示公式计算Dlim与SPElim:
上两式中,置信水平α=99%,Fα(d,n-d)表示自由度为d与n-d的F分布,表示权重为g=a/2b,自由度为h=2a2/b的χ2分布,a与b分别是对应的统计量SPE的估计均值和估计方差。
(8)收集新采样时刻的数据样本x∈R1×m,对其实施与步骤(1)中相同的标准化处理得到新数据向量
(9)根据步骤(2)中的C个不同的变量子块,对应地将样本向量分成C个不同的子向量x1,x2,…,xC。
(10)调用步骤(4)中建立的第c个PLS模型,根据如下所示公式计算出多生产单元变量交叉相关解耦后的误差ec
ec=xc-yΘc (5)
上式中,y=[x1,…,xc-1,xc+1,,xC]。
(11)重复步骤(10)直至得到C个误差e1,e2,…,eC,并将这些误差合并成一个误差向量e=[e1,e2,…,eC]。
(12)对误差e实施与步骤(6)中相同的标准化处理得到新向量
(13)根据如下所示公式计算监测统计指标D与SPE:
上式中,I为m×m维的单位矩阵。
(14)判断是否满足条件:D≤Dlim且SPE≤SPElim?若是,则当前监测时刻***处于正常运行状态,返回步骤(8)继续实施在线故障监测;若否,则当前采样数据来自故障工况。
与传统方法相比,本发明方法的优势在于:
本发明方法利用回归模型将不同生产单元变量子块之间的交叉关系考虑进来,而不是直接建立分布式多个故障监测模型。由于回归模型的误差包含的是与模型输入不相关的信息,本发明方法利用的误差信息实为多生产单元变量交叉相关解耦后的成分信息,能够反映出不同生产单元之间交叉相关关系是否发生变化。此外,本发明方法对误差进行基于PCA算法的建模与监测,又能挖掘出误差信息之间的相关性特征。因此,本发明方法理应具备更优秀的故障监测性能,是一种更为优选的数据驱动故障监测方法。
附图说明
图1为本发明方法实施流程图。
图2为TE过程结构流程图。
图3为交叉相关解耦示意图。
图4为TE过程物料C进口温度故障变化的故障监测详情对比图。
具体实施方式
下面结合附图与具体的实施案例对本发明方法进行详细的说明。
如图1所示,一种基于多生产单元变量交叉相关解耦策略的故障监测方法。下面结合一个具体的工业过程的例子来说明本发明方法的具体实施过程,以及相对于传统分布式PCA方法的优越性。
应用对象是来自于美国田纳西-伊斯曼(TE)化工过程实验,原型是伊斯曼化工生产车间的一个实际工艺流程。如图2所示,TE过程的生产流程比较复杂性,包含5个主要的生产单元:反应器、冷凝器、分离塔、汽提塔、和压缩机。TE过程已作为一个标准实验平台被广泛用于故障监测研究,整个TE过程包括22个测量变量、12个操作变量、和19个成分测量变量。本次选取如表1所示的33个过程变量实施离线建模与在线故障监测,接下来结合该TE过程对本发明具体实施步骤进行详细的阐述。
表1:TE过程监测变量。
序号 | 变量描述 | 序号 | 变量描述 | 序号 | 变量描述 |
1 | 物料A流量 | 12 | 分离器液位 | 23 | D进料阀门位置 |
2 | 物料D流量 | 13 | 分离器压力 | 24 | E进料阀门位置 |
3 | 物料E流量 | 14 | 分离器塔底流量 | 25 | A进料阀门位置 |
4 | 总进料流量 | 15 | 汽提塔等级 | 26 | A和C进料阀门位置 |
5 | 循环流量 | 16 | 汽提塔压力 | 27 | 压缩机循环阀门位置 |
6 | 反应器进料 | 17 | 汽提塔底部流量 | 28 | 排空阀门位置 |
7 | 反应器压力 | 18 | 汽提塔温度 | 29 | 分离器液相阀门位置 |
8 | 反应器等级 | 19 | 汽提塔上部蒸汽 | 30 | 汽提塔液相阀门位置 |
9 | 反应器温度 | 20 | 压缩机功率 | 31 | 汽提塔蒸汽阀门位置 |
10 | 排空速率 | 21 | 反应器冷却水出口温度 | 32 | 反应器冷凝水流量 |
11 | 分离器温度 | 22 | 分离器冷却水出口温度 | 33 | 冷凝器冷却水流量 |
首先,利用TE过程正常工况下的采样数据离线建立故障监测模型,包括以下步骤:
(1):收集生产过程正常运行状态下的数据样本,组成训练数据集X∈R960×33,并对每个变量进行标准化处理,得到均值为0,标准差为1的新数据矩阵
表2:多生产单元变量子块划分结果。
生产单元名称: | 变量子块中的变量编号: |
反应器 | 1,2,3,6,7,8,9,21,23,24,25,32 |
冷凝器 | 33 |
分离塔 | 4,15,16,17,18,19,26,30,31 |
汽提塔 | 10,11,12,13,14,22,28,28 |
压缩机 | 5,20,27 |
(2):根据各测量变量的生产单元归属,将33个测量变量分成5个不同的变量子块,相应的变量子块划分结果如表2所示。
(3):依据5个不同的变量子块将矩阵对应地分成5个不同的子矩阵X1,X2,…,X5。
(4):将第c个子矩阵Xc做为回归模型的输出,同时将Y=[X1,…,Xc-1,Xc+1,…X5]做为回归模型的输入,利用PLS算法建立输入-输出之间的回归模型:Xc=YΘc+Ec。
(5):重复步骤(4)直至得到5个PLS模型,并将模型误差E1,E2,…,E5合并成一个误差矩阵E=[E1,E2,…,E5]。
为验证交叉相关解耦特性,特将E5与输入Y=[X1,X2,X3,X4]之间的交叉相关解耦示意图显示于图3中。从图3中可以发现,压缩机生产单元对应的误差E5与其他生产单元的测量变量间的时序相关性几乎全部为0,即体现出了交叉相关解耦。
(6):将估计误差矩阵E作为新的训练数据矩阵,对其实施标准化处理。
(7):为建立相应的PCA故障监测模型,并保留模型参数集Θ={B,Λ,Dlim,SPElim}。
利用TE过程物料C进口温度故障变量工况的采样数据作为测试数据,实施在线故障监测。值得指出的是,该测试数据集前160个数据为正常工况,故障工况从161个数据采样时刻引入。
(8):收集新采样时刻的数据样本x∈R1×33,对其实施与步骤(1)中相同的标准化处理得到新数据向量后,初始化i=1。
(9):根据步骤(2)中的5个不同的变量子块,对应地将样本向量对应地分成5个不同的子向量x1,x2,…,x5;
(10):调用步骤(4)中建立的第c个PLS模型,根据如下所示公式计算出多生产单元变量交叉相关解耦后的误差ec。
(11):重复步骤(10)直至得到5个误差e1,e2,…,e5,并将这些误差合并成一个误差向量e=[e1,e2,…,e5]。
(12):对误差e实施与步骤(6)中相同的标准化处理得到新向量
(13):计算监测统计指标D与SPE。
(14):判断是否满足条件:D≤Dlim且SPE≤SPElim?若是,则当前监测时刻***处于正常运行状态,返回步骤(8)继续实施在线故障监测;若否,则当前采样数据来自故障工况。
如图4所示,本发明方法以及传统分布式PCA方法在该测试故障上的监测详情对比图,可以很明显地发现本发明方法的故障漏报率明显低于传统分布式PCA方法。因此,本发明方法的监测效果要优越于传统分布式PCA故障监测方法。
上述实施案例只用来解释说明本发明的具体实施,而不是对本发明进行限制。在本发明的精神和权利要求的保护范围内,对本发明做出的任何修改,都落入本发明的保护范围。
Claims (3)
1.一种基于多生产单元变量交叉相关解耦策略的故障监测方法,其特征在于,包括以下步骤:离线建模的实施过程如下所示:
步骤(1):采集生产过程对象中所有生产单元正常运行状态下的样本数据,组成训练数据集X∈Rn×m,并对其进行标准化处理,得到均值为0、标准差为1的新数据矩阵其中,n为训练样本总数,m为过程对象所有测量变量的个数,R为实数集,Rn×m表示n×m维的实数矩阵;
步骤(2):根据各测量变量的生产单元归属,将m个测量变量分成C个不同的变量子块,其中C表示过程对象中生产单元的个数;
步骤(3):依据C个不同的变量子块将矩阵对应地分成C个不同的子矩阵X1,X2,…,XC,其中c=1,2,…,C,mc为第c个变量子块中变量个数且满足条件∑mc=m;
步骤(4):将第c个子矩阵Xc做为回归模型的输出,同时将Y=[X1,…,Xc-1,Xc+1,…XC]做为回归模型的输入,利用偏最小二乘算法(Partial Least Square,PLS)建立如下所示的输入-输出之间的回归模型:
Xc=YΘc+Ec (1)
上式中,Θc为回归矩阵,为回归模型误差;
步骤(5):重复步骤(4)直至得到C个PLS模型,并将模型误差E1,E2,…,EC合并成一个误差矩阵E=[E1,E2,…,EC]∈Rn×m;
步骤(6):将估计误差矩阵E作为新的训练数据矩阵,对其中的每一列实施标准化处理,得到均值为0,标准差为1的新数据矩阵
步骤(7):利用主成分分析(Principal Component Analysis,PCA)算法为建立相应的PCA故障监测模型,并保留模型参数集Θ={B,Λ,Dlim,SPElim}以备调用,其中B为投影变量矩阵,Λ是由特征值组成的对角矩阵,Dlim与SPElim分别表示监测统计量D与SPE的控制上限;
在线故障监测的实施过程如下所示:
步骤(8):收集新采样时刻的数据样本x∈R1×m,对其实施与步骤(1)中相同的标准化处理得到新数据向量
步骤(9):根据步骤(2)中的C个不同的变量子块,对应地将样本向量分成C个不同的子向量x1,x2,…,xC;
步骤(10):调用步骤(4)中建立的第c个PLS模型,根据如下所示公式计算出多生产单元变量交叉相关解耦后的误差ec
ec=xc-yΘc (5)
上式中,y=[x1,…,xc-1,xc+1,,xC];
步骤(11):重复步骤(10)直至得到C个误差e1,e2,…,eC,并将这些误差合并成一个误差向量e=[e1,e2,…,eC];
步骤(12):对误差e实施与步骤(6)中相同的标准化处理得到新向量
步骤(13):根据如下所示公式计算监测统计指标D与SPE:
上式中,I为m×m维的单位矩阵;
步骤(14):判断是否满足条件:D≤Dlim且SPE≤SPElim?若是,则当前监测时刻***处于正常运行状态,返回步骤(8)继续实施在线故障监测;若否,则当前采样数据来自故障工况。
2.根据权利要求1所述的一种基于多生产单元变量交叉相关解耦策略的故障监测方法,其特征在于,所述步骤(4)中利用偏最小二乘算法建立由输入Y与输出Xc之间回归模型的实施步骤具体为:
①置k=1与Z=Y后,初始化向量uk为矩阵Xc的第一列;
②根据公式wk=Zuk/||Zuk||、sk=Zwk、qk=Xc Tsk/(sk Tsk)分别计算得到系数向量wk、得分向量sk、和系数向量qk,其中||Zuk||表示计算向量Ziuk的长度,上标号T表示矩阵或者向量的转置;
③根据公式unew=Xcqk/qk 2计算向量unew;
④判断是否满足条件||uk-unew||<10-6?若否,则置uk=unew后返回步骤②;若是,则执行⑤;
⑤根据公式pk=ZTsk/(sk Tsk)计算得到第k个投影向量pk∈R(m-1)×1,并保留向量pk、向量wk、和系数向量qk;
⑥判断矩阵Yk=skpk T中的最大元素是否大于0.01?若是,根据公式Z=Z-skpk T更新矩阵Zi后执行步骤⑦;若否,则得到投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2,…,wk]、和向量Q=[q1,q2,…,qk],并执行步骤⑧;
⑦判断k<m-mc?若是,则置k=k+1后,返回步骤②;若否,则得到最终的投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2,…,wk]、和向量Q=[q1,q2,…,qk];
⑧根据公式br=Wr(Pr TWr)-1Qr T计算在保留r个投影向量前提下的回归系数向量br,其中r=1,2,…,k,Pr、Wr、和Qr分别为矩阵P、W、和Q中前1至r列向量组成的矩阵;
⑨根据公式计算输出Xc的预测值后,计算保留不同个数投影向量所对应的模型预测误差其中diag[ ]表示将矩阵对角线的元素单独组成向量的操作,sum{ }表示求取向量中各元素之和;
⑩根据MSE1,MSE2,…,MSEk的数值变化情况找出预测误差不再发生显著变化的收敛点位置,并将其下标号对应的具体数值作为偏最小二乘模型需保留的投影变量个数K,其中K≤k;
根据K确定最终的输入Y与输出Xc之间的PLS模型,即:
Xc=YWK(Pk TWk)-1Qk T+Ec=YΘc+Ec (1)
上式中,回归矩阵Θc=Wk(Pk TWK)-1QK T,为回归模型误差。
3.根据权利要求1所述的一种基于多生产单元变量交叉相关解耦策略的故障监测方法,其特征在于,所述步骤(7)中建立PCA故障监测模型的详细实施过程具体为:
①计算的协方差矩阵
②求解S所有特征值γ1≥γ2≥…≥γm所对应的特征向量β1,β2…,βm;
③设置保留的主成分个数d为满足如下所示条件的最小值,将对应的d个最大特征值组成对角矩阵Λ∈Rd×d,并将对应的d个特征向量组成投影变换矩阵B=[β1,β2…,βd]
④根据如下所示公式计算Dlim与SPElim:
上两式中,置信水平α=99%,Fα(d,n-d)表示自由度为d与n-d的F分布,表示权重为g=a/2b,自由度为h=2a2/b的χ2分布,a与b分别是对应的统计量SPE的估计均值和估计方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810280187.3A CN108536943B (zh) | 2018-03-19 | 2018-03-19 | 一种基于多生产单元变量交叉相关解耦策略的故障监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810280187.3A CN108536943B (zh) | 2018-03-19 | 2018-03-19 | 一种基于多生产单元变量交叉相关解耦策略的故障监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108536943A true CN108536943A (zh) | 2018-09-14 |
CN108536943B CN108536943B (zh) | 2021-09-21 |
Family
ID=63482142
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810280187.3A Active CN108536943B (zh) | 2018-03-19 | 2018-03-19 | 一种基于多生产单元变量交叉相关解耦策略的故障监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108536943B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111695229A (zh) * | 2019-03-12 | 2020-09-22 | 宁波大学 | 一种基于ga-ica的新型分散式非高斯过程监测方法 |
CN113283315A (zh) * | 2021-05-12 | 2021-08-20 | 中国航空工业集团公司上海航空测控技术研究所 | 一种多飞行状态直升机故障动态监测告警方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104714537A (zh) * | 2015-01-10 | 2015-06-17 | 浙江大学 | 一种基于联合相对变化分析和自回归模型的故障预测方法 |
CN106649202A (zh) * | 2016-12-07 | 2017-05-10 | 宁波大学 | 基于多样性变量加权plsr模型的工业过程软测量方法 |
CN106940808A (zh) * | 2017-04-28 | 2017-07-11 | 宁波大学 | 一种基于改进型主元分析模型的故障检测方法 |
BR102015026814A2 (pt) * | 2015-10-22 | 2017-07-11 | Universidade Federal Do Paraná | Process of ph estimation in palmito |
CN107045283A (zh) * | 2011-09-19 | 2017-08-15 | 费希尔-罗斯蒙特***公司 | 使用多阶段数据分离的推理过程建模、质量预测及故障检测 |
-
2018
- 2018-03-19 CN CN201810280187.3A patent/CN108536943B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107045283A (zh) * | 2011-09-19 | 2017-08-15 | 费希尔-罗斯蒙特***公司 | 使用多阶段数据分离的推理过程建模、质量预测及故障检测 |
CN104714537A (zh) * | 2015-01-10 | 2015-06-17 | 浙江大学 | 一种基于联合相对变化分析和自回归模型的故障预测方法 |
BR102015026814A2 (pt) * | 2015-10-22 | 2017-07-11 | Universidade Federal Do Paraná | Process of ph estimation in palmito |
CN106649202A (zh) * | 2016-12-07 | 2017-05-10 | 宁波大学 | 基于多样性变量加权plsr模型的工业过程软测量方法 |
CN106940808A (zh) * | 2017-04-28 | 2017-07-11 | 宁波大学 | 一种基于改进型主元分析模型的故障检测方法 |
Non-Patent Citations (4)
Title |
---|
AHMAD M. AL-ZAIDY等: "Novel approach of satellite health monitoring, diagnosis and prediction via PLS batch modelling", 《IEEE XPLORE》 * |
M. ZIYAN SHERIFF等: "Fault detection of nonlinear systems using an improved KPCA method", 《IEEE XPLORE》 * |
彭开香等: "复杂工业过程质量相关的故障检测与诊断技术综述", 《自动化学报》 * |
童楚东等: "基于互信息的PCA方法及其在过程监测中的应用", 《化工学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111695229A (zh) * | 2019-03-12 | 2020-09-22 | 宁波大学 | 一种基于ga-ica的新型分散式非高斯过程监测方法 |
CN111695229B (zh) * | 2019-03-12 | 2023-10-17 | 宁波大学 | 一种基于ga-ica的新型分散式非高斯过程监测方法 |
CN113283315A (zh) * | 2021-05-12 | 2021-08-20 | 中国航空工业集团公司上海航空测控技术研究所 | 一种多飞行状态直升机故障动态监测告警方法 |
CN113283315B (zh) * | 2021-05-12 | 2022-04-05 | 中国航空工业集团公司上海航空测控技术研究所 | 一种多飞行状态直升机故障动态监测告警方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108536943B (zh) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lu et al. | Sub‐PCA modeling and on‐line monitoring strategy for batch processes | |
CN108375965A (zh) | 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 | |
CN108919755A (zh) | 一种基于多块非线***叉关系模型的分布式故障检测方法 | |
CN108803520A (zh) | 一种基于变量非线性自相关性剔除的动态过程监测方法 | |
CN108508865A (zh) | 一种基于分散式osc-pls回归模型的故障检测方法 | |
CN106092625B (zh) | 基于修正型独立元分析和贝叶斯概率融合的工业过程故障检测方法 | |
CN109189028A (zh) | 基于多块信息提取的pca故障诊断方法 | |
JP2010506256A (ja) | 多変量データの監視および解析のための単変量方法 | |
CN105955219A (zh) | 基于互信息的分布式动态过程故障检测方法 | |
CN108388234A (zh) | 一种基于相关性划分多变量块pca模型的故障监测方法 | |
CN109062189A (zh) | 一种用于复杂故障的工业过程故障诊断方法 | |
CN107153409A (zh) | 一种基于缺失变量建模思路的非高斯过程监测方法 | |
CN110009020A (zh) | 一种基于多核主元分析模型的非线性过程监测方法 | |
CN108445867A (zh) | 一种基于分散式icr模型的非高斯过程监测方法 | |
CN108345284A (zh) | 一种基于两变量块的质量相关故障检测方法 | |
CN108153267A (zh) | 一种基于误差主元分析模型的工业过程监测方法 | |
CN106647650A (zh) | 基于变量加权pca模型的分散式工业过程监测方法 | |
CN108345294A (zh) | 一种基于分散式主元回归模型的故障检测方法 | |
CN108536943A (zh) | 一种基于多生产单元变量交叉相关解耦策略的故障监测方法 | |
CN109799808A (zh) | 一种基于重构技术的动态过程故障预测方法 | |
Lindner et al. | Data-driven fault detection with process topology for fault identification | |
CN108762242A (zh) | 一种基于多块典型相关分析模型的分布式故障检测方法 | |
CN106354125A (zh) | 一种利用分块pca检测化工过程故障的方法 | |
CN109766281A (zh) | 一种故障检测率下降变化的不完美调试软件可靠性模型 | |
CN108181893A (zh) | 一种基于pca-kdr的故障检测方法 |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20180914 Assignee: Ningbo Science and Technology Innovation Association Assignor: Ningbo University Contract record no.: X2023980033633 Denomination of invention: A Fault Monitoring Method Based on Multiple Production Cell Variable Cross Correlation Decoupling Strategy Granted publication date: 20210921 License type: Common License Record date: 20230317 |