一种基于量纲可变型独立元分析模型的故障检测方法技术

本发明专利技术公开一种基于量纲可变型独立元分析模型的故障检测方法，该方法在离线建模阶段首先按照传统独立元分析(ICA)模型中分离矩阵各列向量元素数值的差异，对应为各变量赋予不同的权重以体现量纲的差异。然后，由于分离矩阵中每个列向量都体现了测量变量在该投影方向上的差异，因此可有多种量纲可变形式，对应可建立多个ICA故障检测模型。在线监测时，则调用这多个ICA模型计算相应的监测统计量，并利用贝叶斯推理得到最终的概率型监测指标以方便故障决策。与传统方法相比，本发明专利技术方法不仅将测量变量的不同等重要性考虑进建模中，而且还利用多模型实施故障检测。该方法对正常数据特征的描述就更加全面，利用取得更优越的故障检测效果。

【技术实现步骤摘要】
一种基于量纲可变型独立元分析模型的故障检测方法
本专利技术涉及一种工业过程故障检测方法，尤其涉及一种基于量纲可变型独立元分析模型的故障检测方法。
技术介绍
随着计算机技术在工业系统中的广泛应用，现代化工业生产过程正在迈向“大数据时代”，海量的生产过程数据资源为数据驱动的监测、软测量、控制、优化等提供坚实的数据基础。其中，数据驱动的故障检测方法于近年来得到了空前的发展，各种故障检测模型层出不穷。以主成分分析(PrincipalComponentAnalysis，PCA)与独立成分分析(IndependentComponentAnalysis，ICA)算法为代表的多变量统计过程监测方法得到了学术界和工业界的广泛关注，产生了适合不同工况类型的故障检测实施方案。与基于PCA的故障检测方法相比较而言，ICA算法在建模的过程中旨在提取在高阶统计意义上相互独立的成分信息，可应对非高斯的过程数据对象。正因为如此，基于ICA的故障检测方法在监测现代化复杂工业过程对象时，其所取得故障检测效果通常优越于传统的PCA方法。值得指出的是，无论是基于PCA，还是基于ICA的故障检测模型，它们在离线建模阶段都需要对正常工况下的采样数据实施标准化预处理以消除测量变量量纲间的差异。由于任何一个测量变量对应的采样数据出现问题，都可以说过程进入非正常的运行状态，因此每个测量变量的重要性是一致的。在建模前为了消除各个测量变量变化范围(量纲)对建模结果的影响，标准化处理将所有测量变量的量纲统一化是一个势在必行的过程。可是，从故障检测的角度出发，若能将测量变量的量纲差异体现出来，就相当于为不同变量赋予了不同的权重系数，相应的故障检测模型就会对某些故障种类更为敏感，从而可以有效地提升故障检测效果。由于故障检测模型在离线建模阶段时只利用了正常工况下的采样数据集，可以看成是一种无监督型的建模思路。因此，只有正常数据而缺乏故障数据的条件下，无法有针对性的将量纲差异体现出来，这时标准化处理貌似就成了“不得已而为之”的预处理方案了。为了提高传统基于ICA的故障检测模型的监测效果，量纲可变是一种可行的实施方案。唯一的难点在于如何在只有正常数据可用来建模的前提下，有指导性的挖掘量纲差异，从而在此基础上建立量纲可变的故障检测模型。
技术实现思路
本专利技术所要解决的主要技术问题是：从正常数据统计特性角度出发，挖掘测量变量间的量纲差异，并建立基于量纲可变型的ICA故障检测模型以实施在线故障检测。具体来讲，本专利技术方法在离线建模阶段首先按照传统ICA建模的思路得到模型的分离矩阵，然后依据分离矩阵各列向量上元素数值的差异，对应地为过程测量变量赋予不同的权重以体现量纲的差异，并再次建立ICA故障检测模型。由于分离矩阵有多个列向量，每个列向量都体现了测量变量在该投影方向上的差异，对应有多种不同的量纲可变形式，可建立多个ICA故障检测模型。因此，在线实施故障检测时，则调用这多个ICA模型计算相应的监测统计量，并利用贝叶斯推理得到最终的概率型监测指标，从而简化最后是否触发故障报警的决策。本专利技术解决上述技术问题所采用的技术方案为：一种基于量纲可变型独立元分析模型的故障检测方法，包括以下步骤：(1)采集生产过程正常运行状态下的数据样本，组成训练数据集X∈Rn×m，并对每个变量进行标准化处理，得到均值为0，标准差为1的新数据矩阵其中，n为训练样本数，m为过程测量变量数，R为实数集，Rn×m表示n×m维的实数矩阵。(2)利用PCA方法对数据矩阵进行白化处理得到数据矩阵Z∈Rn×M。其中，M≤m为矩阵Z中变量个数。(3)利用ICA算法为训练数据求解得到分离矩阵W∈Rm×p，其中p为分离矩阵中的列向量的个数。(4)根据分离矩阵W＝[w1，w2，…，wp]中各个列向量wj∈Rm×1元素的差异，对应为矩阵中各列赋予不同的权重，相应的得到数据矩阵Xj，其中，下标号j＝1，2，…，p。(5)再次利用ICA算法分别为量纲改变后的训练数据矩阵X1，X2，…，Xp求解分离矩阵(W1，W2，…，Wp)与混合矩阵(A1，A2，…，Ap)，建立对应的p个ICA模型。(6)确定每个ICA模型故障监测统计量所对应的控制限，保留模型参数集Θ1，Θ2，…，Θp以备在线监测调用。(7)采集过程对象新采样时刻的数据xnew∈Rm×1，对其实施与训练数据集X相同的标准化处理得到(8)根据分离矩阵W中各个列向量wj，对中各测量变量实施相同的加权处理，对应得到量纲改变后的数据向量y1，y2，…，yp。(9)调用模型参数集Θ1，Θ2，…，Θp实施在线故障检测，并利用贝叶斯推理将多个监测统计量融合为一个概率型指标方便决策故障发生与否。与传统方法相比，本专利技术方法的优势主要在于：首先，本专利技术方法通过对不同变量赋予不同的权重，使各测量变量的量纲不一致，从而实现了将测量变量的不同等重要性考虑进建模中。此外，本专利技术方法实施量纲可变型建模时，采用了多组不同的量纲可变方式以建立多个ICA故障检测模型。相比于传统单个ICA模型的故障检测方法而言，本专利技术方法是用多模型实施故障检测的，对正常数据特征的描述就更加全面。因此，本专利技术方法可以提升传统基于ICA的故障检测方法的过程监测性能。值得指出的是，本专利技术方利用ICA算法建立故障检测模型主要考虑了过程数据的非高斯特性，本专利技术方法当然也可以借鉴其他非线性的建模算法，以实现对非线性过程对象的故障检测。附图说明图1为本专利技术方法的实施流程图。图2为ICA算法求解过程示意图。具体实施方式下面结合附图对本专利技术方法进行详细的说明。如图1所示，本专利技术公开一种基于量纲可变型独立元分析模型的故障检测方法。具体的实施步骤如下所示：步骤1：采集生产过程正常运行状态下的数据样本，组成训练数据集X∈Rn×m，并对每个变量进行标准化处理，得到均值为0，标准差为1的新数据矩阵其中，n为训练样本数，m为过程测量变量数，R为实数集，Rn×m表示n×m维的实数矩阵。步骤2：利用PCA方法对数据矩阵进行白化处理得到数据矩阵Z∈Rn×M。其中，M≤m为矩阵中变量个数，具体的实施过程如下所示：①计算的协方差矩阵其中S∈Rm×m；②计算矩阵S的所有特征值和特征向量，并剔除小于0.0001的特征值及其对应的特征向量，得到特征向量组成的矩阵B＝[b1，b2，…，bM]∈Rm×M以及特征值组成的对角矩阵D∈RM×M；③按照公式对进行白化处理以得到矩阵Z；步骤3：利用ICA算法为训练数据求解得到分离矩阵W∈Rm×p，其中p为分离矩阵中的列向量的个数，相应的实施过程如图2所示，具体包括以下流程：①设定需要提取的独立成分个数p，并初始化k＝1；②选取单位矩阵I∈RM×M中的第k列做为向量ck的初始值；③根据公式ck←E{Zg(ckTZ)}-E{h(ckTZ)}ck更新向量ck，其中，g(ckTZ)与h(ckTZ)的计算方式如下所示：g(ckTZ)＝4·(ckTZ)3(1)h(ckTZ)＝12·(ckTZ)2(2)E{}表示求取向量的平均值；④对更新后的向量ck依次按照下式进行正交归一化处理：ck←ck/||ck||(4)⑤重复③～④直至向量ck中各元素不再发生变化，然后保存向量ck⑥判断k＜p？若是，置k＝k+1后返回②；若否，将保存的d个向量c1，c2，…本文档来自技高网...

【技术保护点】
一种基于量纲可变型独立元分析模型的故障检测方法，其特征在于，包括以下步骤：(1)：采集生产过程正常运行状态下的数据样本，组成训练数据集X∈R

【技术特征摘要】
1.一种基于量纲可变型独立元分析模型的故障检测方法，其特征在于，包括以下步骤：(1)：采集生产过程正常运行状态下的数据样本，组成训练数据集X∈Rn×m，并对每个变量进行标准化处理，得到均值为0，标准差为1的新数据矩阵其中，n为训练样本数，m为过程测量变量数，R为实数集，Rn×m表示n×m维的实数矩阵；(2)：利用主成分分析(PrincipalComponentAnalysis，PCA)方法对数据矩阵进行白化处理得到数据矩阵Z∈Rn×M。其中，M≤m为矩阵中变量个数；(3)：利用独立元分析(IndependentComponentAnalysis，ICA)算法为训练数据求解得到分离矩阵W∈Rm×p，其中p为分离矩阵中的列向量的个数；(4)：根据分离矩阵W＝[w1，w2，…，wp]中各个列向量wj∈Rm×1元素的差异，对应为矩阵中各列赋予不同的权重，按照如下所示公式得到数据矩阵Xj：其中，下标号j＝1，2，…，p，对角线上的元素wj，1，wj，2，…，wj，m分别是向量wj中的第1，2，…，m个元素；(5)：再次利用ICA算法分别为量纲改变后的训练数据矩阵X1，X2，…，Xp求解分离矩阵(W1，W2，…，Wp)与混合矩阵(A1，A2，…，Ap)，具体的实施过程与步骤3类似，这里不再赘述；(6)：确定每个ICA模型故障监测统计量所对应的控制限，保留模型参数集Θ1，Θ2，…，Θp以备在线监测调用；(7)：采集过程对象新采样时刻的数据xnew∈Rm×1，对其实施与训练数据集X相同的标准化处理得到(8)：根据分离矩阵W中各个列向量wj，按照如下所示公式对中各测量变量实施相同的加权处理：对应得到量纲改变后的数据向量yj∈R1×m，其中，上标号T表示矩阵或向量的转置；(9)：调用模型参数集Θ1，Θ2，…，Θp实施在线故障检测，并利用贝叶斯推理将多个监测统计量融合为概率指标从而方便决策故障发生与否，具体的实施过程如下所示：①利用参数集Θj中的分离矩阵Wj和Aj，按照如下所示公式计算监测统计量Ij2与Qj，hm：Ij2＝yjWjWjTyjT(3)Qj＝||yj-yjWjAjT||2(4)上式中，符号||||表示求取向量的长度；②按照如下所示公式计算条件概率与上式中，N和F分别表示正常和故障工况；③按照如下所示公式计算概率上式中，β为置信限，通常取β＝99％；④按照如下所示公式计算新数据属于故障的概率

【专利技术属性】
技术研发人员：史旭华，蓝艇，童楚东，
申请(专利权)人：宁波大学，
类型：发明
国别省市：浙江,33