一种基于主成分追踪的工业故障监测方法及应用技术

本发明专利技术公开了一种基于主成分追踪的工业故障监测方法及应用，属于工业过程监控与诊断技术领域。首先，对工业采集数据进行主成分追踪方法分解得到包含过程全部运行信息的低秩矩阵以及包含传感器噪声和过程故障的稀疏矩阵。其次，在低秩矩阵中利用T2统计量进行故障检测；在稀疏矩阵中利用均值相关系数的统计量进行故障检测。本发明专利技术根据主成分追踪方法对矩阵分解的特点，针对不同特点的矩阵采用不同的统计量，充分利用了数据中包含的有效信息，因此与其它现有的方法相比，本发明专利技术方法对于工业故障检测具有更高的准确率。

【技术实现步骤摘要】
一种基于主成分追踪的工业故障监测方法及应用
本专利技术属于工业过程监控与故障诊断领域，特别涉及一种基于主成分追踪的在线工业故障检测，利用T2统计量和均值相关系数统计量。
技术介绍
工业过程生产是国家经济发展的支柱产业，因此保证生产过程的高效性和稳定性十分重要。而故障检测是达到这个目标关键的一步。传统的故障检测方法有很多，如主元分析(PCA)，偏最小二乘(PLS)，支持向量机(SVM)，人工神将网络(ANN)以及在它们基础上面研究出来的改进方法。主成分追踪方法具有对异常值不敏感，对变量微小变化引起的故障比较敏感，可以解决非线性过程这些优点，并且主成分追踪方法本质为解决一个简单的凸优化问题。相比于其他方法，更加简单方便。主成分追踪将数据矩阵分解成为低秩矩阵、稀疏矩阵两个部分，其中低秩矩阵为不包含传感器噪声和过程故障的数据矩阵，稀疏矩阵为包含传感器噪声和过程故障的数据。由于在实际的工业生产中，各个变量是具有相关性的。因此只有在传感器测量过程中没有噪声干扰时，才能得到具有相关性的数据，也就是低秩矩阵，所以传感器噪声和过程故障被放置在稀疏矩阵中。理论上，可以通过观察稀疏矩阵中的元素幅值来进行故障检测，但是由于噪声的干扰，这么做会导致比较高的误报率。
技术实现思路
为了克服现有技术的不足，本专利技术的目的在于针对主成分追踪算法的特点，提供一种基于主成分追踪的工业故障监测方法及应用。一种基于主成分追踪的工业故障检测方法，主要采用基于T2统计量和基于均值相关系数统计量的主成分追踪算法，步骤如下：步骤一：离线建模，对工业过程采集的离线数据，进行数据归一化运算；利用主成分追踪的方法对归一化之后的数据进行分解得到低秩矩阵和稀疏矩阵；在低秩矩阵中使用奇异值分解得到负荷向量和奇异值矩阵，并且计算得到正常阈值，在稀疏矩阵中计算变量的均值和相关系数，得到正常阈值；步骤二：在线监测，对在线采集的数据首先按照离线建模数据的主元方向投影，构建T2统计量和离线建模得到的阈值比较进行故障检测；对在线监测到的数据和原来的离线建模数据组合成新的矩阵进行主成分追踪分解，计算均值，利用离线建模得到的相关系数进行在线故障检测。步骤一所述的离线建模过程如下：1)利用工业过程采集的离线数据构成X＝[x1,x2,…,xn]T∈Rn×m，其中m表示变量的个数，n表示样本个数，xi∈Rm,i＝1,…,n表示第i个样本；2)对采集的数据进行归一化处理得到X*，归一化采用均值和方法的方法，得到均值为0，方差为1的新矩阵；3)对归一化处理之后的数据X*进行主成分追踪方法分解，得到低秩矩阵A1和稀疏矩阵E1，分别在A1和E1中计算参数和阈值，3.1)低秩矩阵，3.1.1)对低秩矩阵A1进行奇异值分解，得到奇异值矩阵Σ和负荷向量P；3.1.2)根据公式(1)所示计算在离线建模过程中T2统计量的阈值，其中Fα(m,n-m)是指自由度为m和n-m的F-分布的上100α％临界点；3.2)稀疏矩阵，3.2.1)计算稀疏矩阵E1中每个变量和第一个变量的相关系数ci,i＝1,2,...,m-1；3.2.2)计算每个变量与第一个变量在相关系数中所占的权重pi,i＝1,2,...,m-1，如公式(2)中所示：3.2.3)计算离线建模的正常情况下的阈值M，将权重pi和每个变量的均值di,i＝1,2,...,m相乘得到阈值M，如公式(3)所示：M＝pi×mi(3)。步骤二所述的在线故障检测过程如下：1)在线采集数据Y∈Rn×m，其中n是在线监测样本的样本数，m是变量的个数，yi指的是在线监测数据集的第i行，也就是第i个采样时间的样本；2)根据公式(4)所示，计算在线监测的每个采样点的T2统计量；T2＝yiTPΣ-2PTyi,(4)3)将离线建模阶段的数据X*最后(n-i)×m,i＝1,2,...n行和在线采集到的数据集Y的前i×m行组成一个新的矩阵D；4)对矩阵D进行主成分追踪矩阵分解得到稀疏矩阵E2；5)将权重和稀疏矩阵E2相乘得到监测统计量G，如公式(5)所示：G＝pi×E2(5)；6)将矩阵G中最后一个变量的值存入统计量MG；7)重复上述3)～6)步骤，重复次数与在线监测数据的采样样本数相同，得到在线监测统计量；8)如果T2统计量超出离线建模阶段的正常阈值或者MG超出正常范围阈值M，则表示出现了故障。所述的工业故障为高炉冶炼过程故障。一种所述的方法用于高炉冶炼过程故障诊断。本专利技术具有以下有益效果：1.本专利技术首次提出一种应用在主成分追踪在线故障检测方法中的T2统计量和均值相关系数统计量，实现对复杂过程的故障检测；2.本专利技术能够针对主成分追踪将数据矩阵分解为低秩矩阵和稀疏矩阵的原理，充分利用两个矩阵中的数据的特点，分别构建出合适的统计量。能够将数据中的信息完全表达出来，提高故障检测的效果。附图说明图1是本专利技术方法的一种流程框图。具体实施方式本专利技术首先，对工业采集数据进行主成分追踪方法分解得到包含过程全部运行信息的低秩矩阵和包含传感器噪声和过程故障的稀疏矩阵。其次，在低秩矩阵中利用T2统计量进行故障检测；在稀疏矩阵中利用均值相关系数的统计量进行故障检测。本专利技术根据主成分追踪方法对矩阵分解的特点，针对包含不同信息的矩阵采用不同的统计量，充分利用了数据中包含的有效信息。因此，利用这种统计量，可以实现基于主成分追踪的在线故障检测。基于主成分追踪的故障检测具有对异常值不敏感，对变量微小变化引起的故障比较敏感的优点，并且仅仅进行了矩阵分解步骤，简单易行。本专利技术提出的一种基于T2统计量和均值相关系数统计量的主成分追踪的故障检测方法，其流程框图如图1所示，包括以下各步骤：步骤一：离线建模1)假设传感器采集m个变量的离线数据，样本个数为n，则得到的离线建模的数据构成X＝[x1,x2,…,xn]T∈Rn×m；2)对采集的数据进行归一化处理得到X*，归一化采用均值和方法的方法，得到均值为0，方差为1的新矩阵；3)对归一化处理之后的数据X*进行主成分追踪方法分解，主成分追踪算法的本质是解决一个凸优化函数问题，该凸优化函数为在限制条件X＝A+E下，最小化||A||*+λ||E||1，得到低秩矩阵A1和稀疏矩阵E1。分别在A1和E1中计算参数和阈值。3.1)低秩矩阵，3.1.1)对低秩矩阵A1进行奇异值分解，得到奇异值矩阵Σ和负荷向量P；3.1.2)根据公式(1)所示计算在离线建模过程中T2统计量的阈值，其中Fα(m,n-m)是指自由度为m和n-m的F-分布的上100α％临界点；3.2)稀疏矩阵，3.2.1)计算稀疏矩阵E1中每个变量和第一个变量的相关系数ci,i＝1,2,...,m-1；3.2.2)计算每个变量与第一个变量在相关系数中所占的权重pi,i＝1,2,...,m-1，如公式(2)中所示：3.2.3)计算离线建模的正常情况下的阈值M，将权重pi和每个变量的均值di,i＝1,2,...,m相乘得到阈值M，如公式(3)所示：M＝pi×mi(3)。步骤二：在线检测1)在线采集数据Y∈Rn×m，其中n是在线监测样本的样本数，m是变量的个数，yi指的是在线监测数据集的第i行，也就是第i个采样时间的样本；2)根据公式(4)所示，计算在线监测的每个采样点的T2统计量；T2＝yiTPΣ-2PTyi,(4)3)将离线建模本文档来自技高网...

【技术保护点】
一种基于主成分追踪的工业故障检测方法，其特征在于，主要采用基于T2统计量和基于均值相关系数统计量的主成分追踪算法，步骤如下：步骤一：离线建模，对工业过程采集的离线数据，进行数据归一化运算；利用主成分追踪的方法对归一化之后的数据进行分解得到低秩矩阵和稀疏矩阵；在低秩矩阵中使用奇异值分解得到负荷向量和奇异值矩阵，并且计算得到正常阈值，在稀疏矩阵中计算变量的均值和相关系数，得到正常阈值；步骤二：在线监测，对在线采集的数据首先按照离线建模数据的主元方向投影，构建T2统计量和离线建模得到的阈值比较进行故障检测；对在线监测到的数据和原来的离线建模数据组合成新的矩阵进行主成分追踪分解，计算均值，利用离线建模得到的相关系数进行在线故障检测。

【技术特征摘要】
1.一种基于主成分追踪的工业故障检测方法，其特征在于，主要采用基于T2统计量和基于变量的均值和相关系数统计量的主成分追踪算法，步骤如下：步骤一：离线建模，对工业过程采集的离线数据，进行数据归一化运算；利用主成分追踪的方法对归一化之后的数据进行分解得到低秩矩阵和稀疏矩阵；在低秩矩阵中使用奇异值分解得到负荷向量和奇异值矩阵，并且计算得到正常阈值，在稀疏矩阵中计算变量的均值和相关系数，得到正常阈值；步骤二：在线监测，对在线采集的数据首先按照离线建模数据的主元方向投影，构建T2统计量和离线建模得到的阈值比较进行故障检测；对在线采集到的数据和原来的离线建模数据组合成新的矩阵进行主成分追踪分解，计算均值，利用离线建模得到的相关系数进行在线故障检测。2.根据权利要求1所述的方法，其特征在于，步骤一所述的离线建模过程如下：1)利用工业过程采集的离线数据构成X＝[x1,x2,…,xn]T∈Rn×m，其中m表示变量的个数，n表示样本个数，xi∈Rm,i＝1,…,n表示第i个样本；2)对采集的数据进行归一化处理得到X*，归一化采用均值和方差的方法，得到均值为0，方差为1的新矩阵；3)对归一化处理之后的数据X*进行主成分追踪方法分解，得到低秩矩阵A1和稀疏矩阵E1，分别在A1和E1中计算参数和阈值，3.1)低秩矩阵，3.1.1)对低秩矩阵A1进行奇异值分解，得到奇异值矩阵Σ和负荷向量P；3.1.2)根据公式(1)所示计算在离线建模过程中T2统计量的阈值，其中Fα...

【专利技术属性】
技术研发人员：杨春节，潘怡君，王琳，孙优贤，安汝峤，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33