一种基于KPCA混合模型的多工况过程监控方法技术

本发明专利技术公开了一种基于KPCA混合模型的多工况过程监控方法，属于工业过程监控与诊断技术领域。本发明专利技术结合高斯混合模型和核主元分析模型，利用核主元分析在处理工业过程非线性、降低数据维度方面的优势以及高斯混合模型在处理非高斯、多工况等问题的优越性能，对工业过程进行监控。相比于现有的其他方法，由于充分考虑了工业过程的非高斯、非线性、多工况的复杂特性，本发明专利技术方法能更精确地估计各工况的统计特性，从而更准确及时地检测出多工况过程的各种故障。

【技术实现步骤摘要】
一种基于KPCA混合模型的多工况过程监控方法
本专利技术属于工业过程监控与故障诊断领域，特别涉及一种基于KPCA混合模型的多工况过程监控方法。
技术介绍
随着工业过程复杂性的增长，工业过程监测和诊断的有效性对于保障生产过程安全、维持产品质量和优化产品利益变得日益重要。对于过程监控和故障诊断问题，传统的方法大多采用多元统计过程监控技术(MultivariableStatisticalProcessMonitoring,MSPM)，其中以主元分析(PrincipalComponentAnalysis,PCA)和偏最小二乘(PartialLeastSquares,PLS)为代表等方法已在工业过程监控中得到了成功的应用。传统的MSPM方法均假设过程数据服从高斯分布，变量之间是线性关系且数据来自单一的操作工况下，但是实际中测量数据难以满足这些假设条件，常呈现非高斯、非线性和多工况等特性。尽管，一些改进的方法，如针对非高斯ICA，核PCA(KPCA)等也被提出。但是，当上述非高斯、非线性和多工况等特性同时存在的时，这些方法仍然无法很好地解决。近年来，基于PCA混合模型的多工况过程监控方法被提出用于解决上述问题。将混合高斯模型和PCA相结合,用EM算法估计模型的工况数以及工况的分布参数和主元数,并对每个子模型构建T2、SPE统计量实现时多工况过程的监控。此种方法对每一个高斯成分模型建立了一个PCA模型，但是传统PCA只能处理变量间的线性关系，并不能提取变量间的非线性信息，而在工业过程，变量间的非线性关系是普遍存在的。由此，在工况存在强非线性情况下，可能会引起故障检测发生错误。
技术实现思路
本专利技术的目的在针对现有技术的不足，提供一种基于KPCA混合模型的多工况过程监控方法，利用核主元分析在处理工业过程非线性、降低数据维度方面的优势以及高斯混合模型在处理非高斯、多工况等问题的优越性能，对工业过程进行监控，从而更准确及时地检测出多工况过程的各种故障。一种基于KPCA混合模型的多工况过程监控方法，该方法的步骤如下：步骤一：离线建模，收集多工况过程正常运行的数据，构建高斯混合模型，在高斯混合模型基础上，对每个高斯元空间进行KPCA变换和降维，建立了基于KPCA的混合模型。采用EM算法来估计模型参数，并按照传统的KPCA方法构建各工况的控制限；步骤二：在线检测，采集在线运行数据，利用步骤一所述的混合模型将监控样本按其后验概率大小分类到相应的工况中，并按照传统的KPCA方法计算统计量，如果该统计量超出步骤一所建立的相对应的控制限，则判断故障发生。步骤一所述的离线建模过程如下：1)利用工业过程采集的多工况监测数据构成X＝[x1,x2,…,xn]T∈Rn×m，其中m表示监控变量的个数，n表示样本个数，xi∈Rm,i＝1,…,n表示第i个样本；2)它在有限GMM模型下的概率密度函数表示为其中K表示GMM中混合高斯成分数目，wi表示第i个单高斯成分的混合系数，且满足和θ＝{θ1,…,θK}分别表示局部和全局高斯模型参数集，即均值向量μi和协方差矩阵Σi。相应的第i个分量的多元高斯密度函数可表示为3)第i个单高斯成分样本子集表示为θi＝{μi,Σi}。对数据集Xi进行KPCA投影。3.1)引入核函数φ将数据集Xi投影到高维特征空间F，表示为Φ:Rm→F(1)3.2)计算核矩阵KKij＝〈Φ(xi),Φ(xj)〉＝K(xi,xj)(2)其中使用径向基核函数σ＝rm，r为常数。3.3)对核矩阵K进行中心化处理其中3.4)计算主成分tk4)采用EM算法来估计模型参数4.1)E-step其中p(s)(Ci|xj)表示第s次迭代后第j个训练样本属于第i个高斯成分的后验概率。4.2)M-step其中，和αi,j分别表示第(s+1)次迭代后，第i个高斯成分的均值，协方差，先验概率和KPCA特征向量。5)确定T2和SPE统计限步骤二所述的在线检测过程如下：1)在线采集测量数据y∈Rm；2)判断采集样本属于的工况类别。计算则xt所属工况为3)计算采集样本的T2和SPE统计量。T2＝[t1,…,tp]Λ-1[t1,…,tp]T(12)4)比较T2、SPE与式(10)、(11)中所建立的检测控制限T2lim、SPElim之间的大小，如果统计量超出控制限，则判断故障发生；如果统计量低于控制限，则说明过程正常运行。一种根据所述的方法用于高炉冶炼过程故障诊断。本专利技术有以下优势：1.本专利技术首次提出一种基于KPCA混合模型的多工况过程监控方法，实现对复杂过程的监测；2.本专利技术能够解决过程数据存在的非高斯、非线性和多模态等问题，从而能够更有效地监控。具体实施方式本专利技术提出的一种基于KPCA混合模型的多工况过程监控方法，包括以下各步骤：步骤一：离线建模1)利用工业过程采集的多工况监测数据构成X＝[x1,x2,…,xn]T∈Rn×m，其中m表示监控变量的个数，n表示样本个数，xi∈Rm,i＝1,…,n表示第i个样本；2)它在有限GMM模型下的概率密度函数表示为其中K表示GMM中混合高斯成分数目，wi表示第i个单高斯成分的混合系数，且满足和θ＝{θ1,…,θK}分别表示局部和全局高斯模型参数集，即均值向量μi和协方差矩阵Σi。相应的第i个分量的多元高斯密度函数可表示为3)第i个单高斯成分样本子集表示为θi＝{μi,Σi}。对数据集Xi进行KPCA投影。3.1)引入核函数φ将数据集Xi投影到高维特征空间F，表示为Φ:Rm→F(1)3.2)计算核矩阵KKij＝〈Φ(xi),Φ(xj)〉＝K(xi,xj)(2)其中使用径向基核函数σ＝rm，r为常数。3.3)对核矩阵K进行中心化处理其中3.4)计算主成分tk4)采用EM算法来估计模型参数4.1)E-step其中p(s)(Ci|xj)表示第s次迭代后第j个训练样本属于第i个高斯成分的后验概率。4.2)M-step其中，和αi,j分别表示第(s+1)次迭代后，第i个高斯成分的均值，协方差，先验概率和KPCA特征向量。5)确定T2和SPE统计限步骤二：在线监测1)在线采集测量数据y∈Rm；2)判断采集样本属于的工况类别。计算则xt所属工况为3)计算采集样本的T2和SPE统计量。T2＝[t1,…,tp]Λ-1[t1,…,tp]T(12)4)比较T2、SPE与式(10)、(11)中所建立的检测控制限T2lim、SPElim之间的大小，如果统计量超出控制限，则判断故障发生；如果统计量低于控制限，则说明过程正常运行。实施例钢铁冶炼作为国民经济中最重要的基础产业之一，是衡量一个国家的经济水平和综合国力的重要指标。而高炉炼铁是钢铁工业生产流程中最重要的环节，所以对大型高炉非正常工况诊断与安全运行方法进行研究具有重要意义。高炉是一个巨大的密闭反应容器，其内部冶炼过程是在高温、高压条件下，经过一系列复杂的物理化学和传热反应，是一个典型的“黑箱”操作。正是由于高炉内部的复杂性，使得其监测过程具有非线性、非高斯性以及多模态等特性。因此，我们提出的方法对高炉故障监测具有适应性。下面结合柳钢2号高炉来说明本专利技术方法的有效性。成立于1958年的柳钢炼铁厂，是一个有着56年辉煌历史的设备先进、装备水平较高的大型冶炼企业，主要产品为生铁，副产品有炉尘、炉渣、高炉本文档来自技高网...

【技术保护点】
一种基于KPCA混合模型的多工况过程监控方法，其特征在于，该方法的步骤如下：步骤一：离线建模，收集多工况过程正常运行的数据，构建高斯混合模型，在高斯混合模型基础上，对每个高斯元空间进行KPCA变换和降维，建立了基于KPCA的混合模型，采用EM算法来估计模型参数，并按照KPCA方法构建各工况的控制限；步骤二：在线检测，采集在线运行数据，利用步骤一所述的混合模型将监控样本按其后验概率大小分类到相应的工况中，并按照KPCA方法计算统计量，如果该统计量超出步骤一所建立的相对应的控制限，则判断故障发生。

【技术特征摘要】
1.一种基于KPCA混合模型的多工况过程监控方法，其特征在于，该方法的步骤如下：步骤一：离线建模，收集多工况过程正常运行的数据，构建高斯混合模型，在高斯混合模型基础上，对每个高斯元空间进行KPCA变换和降维，建立了基于KPCA的混合模型，采用EM算法来估计模型参数，并按照KPCA方法构建各工况的控制限；步骤二：在线检测，采集在线运行数据，利用步骤一所述的混合模型将监控样本按其后验概率大小分类到相应的工况中，并按照KPCA方法计算统计量，如果该统计量超出步骤一所建立的相对应的控制限，则判断故障发生；步骤一所述的离线建模过程如下：1)利用工业过程采集的多工况监测数据构成X＝[x1,x2,…,xn]T∈Rn×m，其中m表示监控变量的个数，n表示样本个数，xi∈Rm,i＝1,…,n表示第i个样本；2)在有限GMM模型下的概率密度函数表示为其中K表示GMM中混合高斯成分数目，wi表示第i个单高斯成分的混合系数，且满足θi＝{μi,∑i}和θ＝{θ1,…,θK}分别表示局部和全局高斯模型参数集，即均值向量μi和协方差矩阵∑i，相应的第i个分量的多元高斯密度函数可表示为3)第i个单高斯成分样本子集表示为θi＝{μi,∑i}，对数据集Xi进行KPCA投影；3.1)引入核函数φ将数据集Xi投影到高维特征空间F，表示为Φ:Rm→F(1)；3.2)计算核矩阵KKij＝<Φ(xi),Φ(xj)>＝K(xi,xj)(2)；其中使用径向基核函数σ＝rm，r为常数；3.3)对核矩阵K进行中心化处理其中3.4)计算主成分tk4)采用EM算法来估计模型参数4.1)E-step

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