CN107632592A - 基于高效递推核主元分析的非线性时变过程故障监测方法 - Google Patents
基于高效递推核主元分析的非线性时变过程故障监测方法 Download PDFInfo
- Publication number
- CN107632592A CN107632592A CN201710777870.3A CN201710777870A CN107632592A CN 107632592 A CN107632592 A CN 107632592A CN 201710777870 A CN201710777870 A CN 201710777870A CN 107632592 A CN107632592 A CN 107632592A
- Authority
- CN
- China
- Prior art keywords
- data
- fault
- monitoring
- vector
- varying
- 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
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于高效递推核主元分析的非线性时变过程故障监测方法,属于故障检测与诊断技术领域。首先,在田纳西伊斯曼过程仿真器中采集具有非线性和缓慢时变特性且包含故障的数据,将采集的正常数据利用高斯核函数将其投影到高维特征空间并中心化,建立初始的离线监测模型并采用核密度估计函数来确定控制限。然后,当采集到新的过程数据时,通过引入一阶干扰理论方法在离线模型获取的特征值和特征矢量基础上直接更新模型,并将新数据投影到更新的核空间和残差空间计算T2和SPE统计量。当超出其对应的控制限时,认为监测的故障发生,反之,整个过程正常运行。本发明方法主要解决两个问题:1)核主元分析进行非线性时变过程的故障监测会产生较高的误报率;2)基于特征值分解的递推算法计算负荷较高的问题。
Description
技术领域
本发明属于故障监测与诊断技术领域,提出了一种基于高效递推核主元分析的非线性时变过程故障监测方法。
背景技术
随着工业过程复杂性的增长,如何提高工业过程系统的安全性、可靠性,防止和避免过程故障的发生已成为一个亟待解决的问题。过程监测就是为解决这一类问题而发展起来的技术。
对于过程监控和故障诊断问题,传统的以主元分析(Principal ComponentAnalysis,PCA)和偏最小二乘(Partial Least Squares,PLS)为代表的多变量统计方法已在工业过程监测领域得到了成功应用。传统的多变量统计方法均假设过程数据来自单一工况且服从高斯分布,变量间是线性关系,但是实际测量数据往往难以满足这些假设条件,常呈现非高斯、非线性、时变和多模态等特性。尽管,一些改进的方法,如针对数据非线性特性的,核主元分析(KPCA);针对过程时变特性的递推主元分析(RPCA)等已被提出;但是,当上述非高斯、非线性和时变等特性同时存在时,已提出的方法仍然无法有效地解决复杂工业过程的故障监测问题。
近年来,基于递推主元分析(RPCA)的非线性时变过程监测方法被提出用于解决上述问题。该方法通过引入一阶干扰理论实现了递推特征值分解,并且降低了递推主元分析(RPCA)的计算负荷。当T2和Q统计量的控制限依然采用基于数据符合高斯分布的假设。此种方法对通过遗忘因子策略对初始PCA模型进行更新,但是传统PCA仅能有效处理线性特性的数据,并不能有效的提取变量间的非线性特征。然而在工业过程中变量间的非线性关系是普遍存在的。因此在存在强非线性情况下,可能会产生较高的故障误报率。
发明内容
发明目的:本发明的目的是为了解决现有技术中的不足,提供一种基于高效递推核主元分析(ERKPCA)的非线性时变过程故障监测方法,以达到降低误报率的目的。利用递推核主成分分能同时处理工业过程非线性和缓慢变化的优势,以及通过引入一阶干扰理论来实现是较低计算负荷的递推奇异值分解。对数据不符合高斯分布的问题,采用核密度估计的方法来确定统计量的控制限,降低了故障误报率。该方法用于工业过程故障监测,不仅能适应过程正常变化,而且也能准确及时地检测出过程的各种故障。
技术方案:本发明所述的基于高效递推核主元分析的非线性时变过程故障监测方法,包括如下步骤:
步骤1:离线建模
采集田纳西伊斯曼化工过程正常运行数据,并将其标准化后用于建立离线的核主元分析模型;通过径向基核函数将训练数据映射到高维线性特征空间,然后利用特征值分解求得样本协方差矩阵的特征值和特征向量,建立田纳西伊斯曼化工过程的离线模型;同时,借助核密度估计算法来确定统计量的控制限;
步骤2:在线监测
当采集到的田纳西伊斯曼化工过程数据xt是正常数据时,采用递推核主元分析方法在步骤2中田纳西伊斯曼化工过程初始离线监测模型提供的特征值和特征向量的基础上,通过引入一阶干扰理论来实现高效的递推特征值分解,从而直接更新特征值和特征向量,并重新计算更新后的田纳西伊斯曼化工过程的负荷向量、残差向量;通过计算Hotelling的T2统计和平方预测误差Q统计量对田纳西伊斯曼化工过程进行故障监测,即可判断出田纳西伊斯曼化工过程是否发生故障,当T2统计和Q统计超出各自的控制限时,认为有故障发生,反之,整个过程正常。
进一步的,步骤1所述的建立田纳西伊斯曼化工过程的初始离线监测模型,方法如下:
1)利用田纳西伊斯曼化工过程采集具有非线性、缓缓时变特性的训练数据集:其中n和m分别是采样个数和变量数,xi∈Rm,i表示第i组采样数据;
2)将X中m个变量分别进行去均值除以标准差的标准化处理,通过高斯核函数K(i,j)=exp(-||x(i)-x(j)||2/σ),将数据映射到高维特征空间,其中σ核函数的宽度参数,通常根据经验设定;
3)对核矩阵K进行中心化处理K=K-ZK-KZ+ZKZ,
其中,
4)在高维线性特征空间求得样本协方差矩阵并通过特征值分解求得特征值λi及其对应的特征向量vi,其中1≤i≤n;
5)根据累计贡献率确定主元个数p:
6)将测试故障数据xt投影到高维特征空间,并计算得分向量:
其中,kt=Φ(Xt)Φ(xt);
7)计算监控统计量T2=[t1,t2…,tp]Λ-1[t1,t2…,tp]T,
其中,Λ-1是p个特征值的逆矩阵,是所有特征值中非零特征值的个数;
8)确定统计量T2和Q的核密度估计:
其中,x和xi是被估计数据向量和向量中对应的第i个数值,n是数据向量包含的总个数,h是核函数的宽度,Kp采用高斯核函数;
9)根据核密度估计方法确定统计量T2和Q的控制限:
其中,α是指定的置信度。
进一步的,步骤2所述高效递推核主元分析方法通过引入一阶干扰理论,在田纳西伊斯曼化工过程初始监测模型的基础上,直接递推更新特征值及其特征向量,高效递推核主元分析监测方法如下:
1)在线采集测量数据xt∈Rm,并利用径向基核函数K(i,j)=exp(-||x(i)-x(j)||2/σ)投影到高维特征空间并进行中心化:
kt=kt-1tkt-kt1t+1tkt1t,其中,1t=1/n[1 … 1]∈R1×n;
2)取得离线模型求得的k-1时刻特征值λk-1及其对应的特征向量vk-1,引入一阶干扰理论实现递推特征值分解直接求得k时刻特征值λk及其对应的特征向量vk;递推特征值分解到步骤如下:
a)计算xt=(ε)1/2xt,其中ε是一个接近于零的小正数;
b)初始化然后计算
c)更新第i个特征值λk,i=(1-ε)λk,i-1+fi 2;
d)更新第i个特征向量
e)并进行标准化vk,i=(vk,i)/(||vk,i||);
3)依据离线模型确定的主元个数来划分负荷向量和残差向量,并计算监控统计量T2和Q,
4)将监控统计量T2和Q与离线模型获取的控制限进行对比,如果或Q>QUCL(α),则表示监测到故障的发生;否则,跳转到步骤1)继续循环执行。
进一步的,所述的工业过程具有复杂特性,包括非线性和缓慢时变特性。
进一步的,所述的工业故障为田纳西伊斯曼化工过程典型故障,其中包含的两类缓慢变化均在正常工况范围内。
有益效果:本发明的有益效果如下:
1.本发明首次提出一种具有较低计算负荷的基于高效递推核主元分析的监测方法,实现对具有非线性和缓慢时变特性的复杂工业过程的故障监测,从而在有效降低故障监测的误报率;
2.本发明能够解决工业过程数据存在的非线性和时变特性等问题,从而能够更有效地适应过程的缓慢变化并且及时准确的检测到过程的突发故障。
具体实施方式
下面结合具体实施例对本发明的技术方案作进一步详细说明。
本发明提出的一种基于高效递推核主元分析的非线性时变过程故障监测方法,包括以下各步骤:
步骤1:离线建模。采集田纳西伊斯曼化工过程正常运行数据,并将其标准化后用于建立离线的核主元分析模型。通过径向基核函数将训练数据映射到高维线性特征空间,然后利用特征值分解求得样本协方差矩阵的特征值和特征向量,建立田纳西伊斯曼化工过程的离线模型;同时,借助核密度估计算法来确定统计量的控制限。
步骤2:在线监测。当采集到的田纳西伊斯曼化工过程数据xt是正常数据时,采用递推核主元分析方法在步骤2中田纳西伊斯曼化工过程初始离线监测模型提供的特征值和特征向量的基础上,通过引入一阶干扰理论来实现高效的递推特征值分解,从而直接更新特征值和特征向量,并重新计算更新后的田纳西伊斯曼化工过程的负荷向量、残差向量。通过计算Hotelling的T2统计和平方预测误差Q统计量对田纳西伊斯曼化工过程进行故障监测,即可判断出田纳西伊斯曼化工过程是否发生故障,当T2统计和Q统计超出各自的控制限时,认为有故障发生,反之,整个过程正常。
步骤一所述的离线建模过程如下:
1)利用田纳西伊斯曼化工过程采集具有非线性、缓缓时变特性的训练数据集其中n和m分别是采样个数和变量数,xi∈Rm,i表示第i组采样数据。
2)将X中m个变量分别进行去均值除以标准差的标准化处理,通过高斯核函数K(i,j)=exp(-||x(i)-x(j)||2/σ),将数据映射到高维特征空间,其中σ核函数的宽度参数,通常根据经验设定;
3)对核矩阵K进行中心化处理K=K-ZK-KZ+ZKZ,
其中,
4)在高维线性特征空间求得样本协方差矩阵并通过特征值分解求得特征值λi及其对应的特征向量vi,其中1≤i≤n;
5)根据累计贡献率确定主元个数p,
6)将测试故障数据xt投影到高维特征空间,并计算得分向量
其中,kt=Φ(Xt)Φ(xt);
7)计算监控统计量T2=[t1,t2…,tp]Λ-1[t1,t2…,tp]T,
其中,Λ-1是p个特征值的逆矩阵,是所有特征值中非零特征值的个数;
8)确定统计量T2和Q的核密度估计
其中,x和xi是被估计数据向量和向量中对应的第i个数值,n是数据向量包含的总个数,h是核函数的宽度,Kp采用高斯核函数;
9)根据核密度估计方法确定统计量T2和Q的控制限
其中,α是指定的置信度。
步骤二所述的在线监测方法如下:
1)在线采集测量数据xt∈Rm,并利用径向基核函数K(i,j)=exp(-||x(i)-x(j)||2/σ)
投影到高维特征空间并进行中心化;
kt=kt-1tkt-kt1t+1tkt1t
其中,1t=1/n[1 … 1]∈R1×n;
2)取得离线模型求得的k-1时刻特征值λk-1及其对应的特征向量vk-1,引入一阶干扰理论实现递推特征值分解直接求得k时刻特征值λk及其对应的特征向量vk;递推特征值分解到步骤如下:
a)计算xt=(ε)1/2xt,其中ε是一个接近于零的小正数;
b)初始化然后计算
c)更新第i个特征值λk,i=(1-ε)λk,i-1+fi 2;
d)更新第i个特征向量
e)并进行标准化vk,i=(vk,i)/(||vk,i||);
3)依据离线模型确定的主元个数来划分负荷向量和残差向量,并计算监控统计量T2和Q,
4)将监控统计量T2和Q与离线模型获取的控制限进行对比。如果或Q>QUCL(α),则表示监测到故障的发生;否则,跳转到步骤1)继续循环执行。
TE过程的实例
化工过程作为国民经济中最重要的支柱产业之一,是衡量一个国家的综合国力和经济水平的重要指标。由于化工过程具有复杂性、非线性、非高斯以及时变和多模态等特性,因此,对大型化工过程进行故障监测保障其安全运行具有重要的实际意义。
田纳西伊斯曼化学公司开发的Tennessee Eastman(TE)过程模拟器,为工业和科研界提供了一个实际工业过程的仿真平台。近年来,被广泛用于验证各种过程建模与控制、故障监测与诊断方法等,成为一种国际上领域内公认的标准仿真模型。TE过程包括五个主要操作单元:反应器、冷凝器、压缩机、分离器和汽提塔;而且还包括八种成分4种气体进料、2个气液放热反应生成的2种主产品、两个衍生放热反应生成的2种副产品等,过程机理复杂,变量较多,包括阶跃、随机变化、慢漂移、粘滞和恒定位置等多种数据故障类型。由于开环不稳定,采用一种分布式控制策略,其中包括9个PI控制环和21种设定值。该过程包含41个测量变量和12个控制变量。每种测量变量都叠加加性噪声来模拟实际工业过程中的噪声。每个数据集有52个变量,大部分变量的采样时间是3分钟;其中有14个变量的采样时间是6分钟,5个变量的采样时间是15分钟。TE过程的Simulink仿真代码可从http://depts.washington.edu/control/LARRY/TE/download.html下载。利用TE过程仿真程序,采集了两组数据集,前200个采样时刻为正常运行数据,第201-500个时刻是过程发生缓慢变化斜坡变化和微小阶跃变化的数据,在第501采样时刻分别引入故障2(B成分,A/C进料比不变)和故障7(C存在压力损失—可用性降低(流4))。
接下来结合具体的TE过程对本发明的实施步骤进行详细地阐述:
离线建模:
1)利用田纳西伊斯曼化工过程采集具有非线性、缓缓时变特性的训练数据集其中n=200和m=52分别是采样个数和变量数,xi∈Rm,i表示第i组采样数据;
2)将X中m个变量分别进行去均值除以标准差的标准化处理,通过高斯核函数K(i,j)=exp(-||x(i)-x(j)||2/σ),将数据映射到高维特征空间,其中σ核函数的宽度参数,通常根据经验设定;
3)对核矩阵K进行中心化处理K=K-ZK-KZ+ZKZ,
其中,
4)在高维线性特征空间求得样本协方差矩阵并通过特征值分解求得特征值λi及其对应的特征向量vi,其中1≤i≤n;
5)根据累计贡献率确定主元个数p,
6)将测试故障数据xt投影到高维特征空间,并计算得分向量
其中,kt=Φ(Xt)Φ(xt);
7)计算监控统计量T2=[t1,t2…,tp]Λ-1[t1,t2…,tp]T,
其中,Λ-1是p个特征值的逆矩阵,是所有特征值中非零特征值的个数;
8)确定统计量T2和Q的核密度估计
其中,x和xi是被估计数据向量和向量中对应的第i个数值,n是数据向量包含的总个数,h是核函数的宽度,Kp采用高斯核函数;
9)根据核密度估计方法确定统计量T2和Q的控制限
其中,α是指定的置信度。
步骤二所述的在线监测方法如下:
1)在线采集测量数据xt∈Rm,并利用径向基核函数K(i,j)=exp(-||x(i)-x(j)||2/σ)投影到高维特征空间并进行中心化;
kt=kt-1tkt-kt1t+1tkt1t
其中,1t=1/n[1 … 1]∈R1×n;
2)取得离线模型求得的k-1时刻特征值λk-1及其对应的特征向量vk-1,引入一阶干扰理论实现递推特征值分解直接求得k时刻特征值λk及其对应的特征向量vk;递推特征值分解到步骤如下:
a)计算xt=(ε)1/2xt,其中ε是一个接近于零的小正数;
b)初始化然后计算
c)更新第i个特征值λk,i=(1-ε)λk,i-1+fi 2;
d)更新第i个特征向量
e)并进行标准化vk,i=(vk,i)/(||vk,i||);
3)依据离线模型确定的主元个数来划分负荷向量和残差向量,并计算监控统计量T2和Q,
4)将监控统计量T2和Q与离线模型获取的控制限进行对比。如果或Q>QUCL(α),则表示监测到故障的发生;否则,跳转到步骤1)继续循环执行。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (5)
1.基于高效递推核主元分析的非线性时变过程故障监测方法,其特征在于:包括如下步骤:
步骤1:离线建模
采集田纳西伊斯曼化工过程正常运行数据,并将其标准化后用于建立离线的核主元分析模型;通过径向基核函数将训练数据映射到高维线性特征空间,然后利用特征值分解求得样本协方差矩阵的特征值和特征向量,建立田纳西伊斯曼化工过程的离线模型;同时,借助核密度估计算法来确定统计量的控制限;
步骤2:在线监测
当采集到的田纳西伊斯曼化工过程数据xt是正常数据时,采用递推核主元分析方法在步骤2中田纳西伊斯曼化工过程初始离线监测模型提供的特征值和特征向量的基础上,通过引入一阶干扰理论来实现高效的递推特征值分解,从而直接更新特征值和特征向量,并重新计算更新后的田纳西伊斯曼化工过程的负荷向量、残差向量;通过计算Hotelling的T2统计和平方预测误差Q统计量对田纳西伊斯曼化工过程进行故障监测,即可判断出田纳西伊斯曼化工过程是否发生故障,当T2统计和Q统计超出各自的控制限时,认为有故障发生,反之,整个过程正常。
2.根据权利要求1所述的基于高效递推核主元分析的非线性时变过程故障监测方法,其特征在于:步骤1所述的建立田纳西伊斯曼化工过程的初始离线监测模型,方法如下:
1)利用田纳西伊斯曼化工过程采集具有非线性、缓缓时变特性的训练数据集:其中n和m分别是采样个数和变量数,xi∈Rm,i表示第i组采样数据;
2)将X中m个变量分别进行去均值除以标准差的标准化处理,通过高斯核函数K(i,j)=exp(-||x(i)-x(j)||2/σ),将数据映射到高维特征空间,其中σ核函数的宽度参数,通常根据经验设定;
3)对核矩阵K进行中心化处理K=K-ZK-KZ+ZKZ,
其中,
4)在高维线性特征空间求得样本协方差矩阵并通过特征值分解求得特征值λi及其对应的特征向量vi,其中1≤i≤n;
5)根据累计贡献率确定主元个数p:
<mrow>
<mi>C</mi>
<mi>P</mi>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
<mo>;</mo>
</mrow>
6)将测试故障数据xt投影到高维特征空间,并计算得分向量:
其中,kt=Φ(Xt)Φ(xt);
7)计算监控统计量T2=[t1,t2…,tp]Λ-1[t1,t2…,tp]T,
其中,Λ-1是p个特征值的逆矩阵,是所有特征值中非零特征值的个数;
8)确定统计量T2和Q的核密度估计:
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>n</mi>
<mi>h</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>h</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,x和xi是被估计数据向量和向量中对应的第i个数值,n是数据向量包含的总个数,h是核函数的宽度,Kp采用高斯核函数;
9)根据核密度估计方法确定统计量T2和Q的控制限:
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>U</mi>
<mi>C</mi>
<mi>L</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mi>dT</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mi>&alpha;</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>U</mi>
<mi>C</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>Q</mi>
<mo>=</mo>
<mi>&alpha;</mi>
</mrow>
其中,α是指定的置信度。
3.根据权利要求1所述的基于高效递推核主元分析的非线性时变过程故障监测方法,其特征在于:步骤2所述高效递推核主元分析方法通过引入一阶干扰理论,在田纳西伊斯曼化工过程初始监测模型的基础上,直接递推更新特征值及其特征向量,高效递推核主元分析监测方法如下:
1)在线采集测量数据xt∈Rm,并利用径向基核函数K(i,j)=exp(-||x(i)-x(j)||2/σ)投影到高维特征空间并进行中心化:
kt=kt-1tkt-kt1t+1tkt1t,其中,1t=1/n[1 … 1]∈R1×n;
2)取得离线模型求得的k-1时刻特征值λk-1及其对应的特征向量vk-1,引入一阶干扰理论实现递推特征值分解直接求得k时刻特征值λk及其对应的特征向量vk;递推特征值分解到步骤如下:
a)计算xt=(ε)1/2xt,其中ε是一个接近于零的小正数;
b)初始化然后计算
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>v</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>v</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msubsup>
<mi>q</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>v</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>q</mi>
<mi>i</mi>
<mi>v</mi>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
c)更新第i个特征值λk,i=(1-ε)λk,i-1+fi 2;
d)更新第i个特征向量
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>/</mo>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>v</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>q</mi>
<mi>i</mi>
<mi>v</mi>
</msubsup>
<mo>;</mo>
</mrow>
e)并进行标准化vk,i=(vk,i)/(||vk,i||);
3)依据离线模型确定的主元个数来划分负荷向量和残差向量,并计算监控统计量T2和Q,
<mrow>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>p</mi>
</msub>
<mo>&rsqb;</mo>
<msup>
<mi>&Lambda;</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>p</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>,</mo>
<mi>Q</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mover>
<mi>n</mi>
<mo>&OverBar;</mo>
</mover>
</munderover>
<msubsup>
<mi>t</mi>
<mi>j</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<msubsup>
<mi>t</mi>
<mi>j</mi>
<mn>2</mn>
</msubsup>
<mo>;</mo>
</mrow>
4)将监控统计量T2和Q与离线模型获取的控制限进行对比,如果或Q>QUCL(α),则表示监测到故障的发生;否则,跳转到11)继续循环执行。
4.根据权利要求1所述的基于高效递推核主元分析的非线性时变过程故障监测方法,其特征在于:所述的工业过程具有复杂特性,包括非线性和缓慢时变特性。
5.根据权利要求1所述的基于高效递推核主元分析的非线性时变过程故障监测方法,其特征在于:所述的工业故障为田纳西伊斯曼化工过程典型故障,其中包含的两类缓慢变化均在正常工况范围内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710777870.3A CN107632592B (zh) | 2017-09-01 | 2017-09-01 | 基于高效递推核主元分析的非线性时变过程故障监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710777870.3A CN107632592B (zh) | 2017-09-01 | 2017-09-01 | 基于高效递推核主元分析的非线性时变过程故障监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107632592A true CN107632592A (zh) | 2018-01-26 |
CN107632592B CN107632592B (zh) | 2021-03-16 |
Family
ID=61099795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710777870.3A Active CN107632592B (zh) | 2017-09-01 | 2017-09-01 | 基于高效递推核主元分析的非线性时变过程故障监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107632592B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108764305A (zh) * | 2018-05-14 | 2018-11-06 | 浙江大学 | 一种改进的群智能机器学习故障诊断系统 |
CN108803337A (zh) * | 2018-06-28 | 2018-11-13 | 杭州电子科技大学 | 一种新型的化工过程数据监测方法 |
CN109145256A (zh) * | 2018-11-14 | 2019-01-04 | 南通大学 | 基于规范变量非线性主成分分析的非线性动态过程监测方法 |
CN109242046A (zh) * | 2018-10-10 | 2019-01-18 | 中国工程物理研究院计算机应用研究所 | 基于非线性复杂系统过程数据的在线故障诊断方法 |
CN109885022A (zh) * | 2019-02-21 | 2019-06-14 | 山东科技大学 | 一种基于潜故障敏感子空间的故障检测方法 |
CN110009126A (zh) * | 2019-01-23 | 2019-07-12 | 北京化工大学 | 基于pls模型与pca贡献度融合的在线报警分析方法 |
CN110224160A (zh) * | 2019-06-12 | 2019-09-10 | 北京亿华通科技股份有限公司 | 一种燃料电池系统故障诊断方法 |
CN110687895A (zh) * | 2019-10-24 | 2020-01-14 | 上海工程技术大学 | 一种基于自适应核主成分分析的化工过程故障检测方法 |
CN110795695A (zh) * | 2019-10-28 | 2020-02-14 | 华北理工大学 | 消去时序相关性的多模型动态pca故障监测方法 |
CN111025222A (zh) * | 2019-12-31 | 2020-04-17 | 三峡大学 | 基于递推主元分析的电子式电压互感器状态监测方法 |
CN111209973A (zh) * | 2020-01-09 | 2020-05-29 | 浙江理工大学 | 基于混合核pca-cca及核密度估计的过程监测方法 |
CN111852837A (zh) * | 2020-07-14 | 2020-10-30 | 北京博华信智科技股份有限公司 | 基于聚类的往复式压缩机自适应故障监测方法 |
CN111914886A (zh) * | 2020-06-13 | 2020-11-10 | 宁波大学 | 一种基于在线简略核学习的非线性化工过程监测方法 |
CN112329191A (zh) * | 2019-08-05 | 2021-02-05 | 华中科技大学 | 一种多模态过程中模态辨识方法和过渡模态故障检测方法 |
CN113253682A (zh) * | 2021-05-17 | 2021-08-13 | 中国石油大学(华东) | 非线性化工过程故障检测方法 |
CN113673884A (zh) * | 2021-08-25 | 2021-11-19 | 国网吉林省电力有限公司长春供电公司 | 基于数字孪生映射模型进行台区数据监测的方法 |
CN114167826A (zh) * | 2021-11-26 | 2022-03-11 | 华中科技大学 | 一种用于铸造生产过程混合多变量的监控方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1655082A (zh) * | 2005-01-27 | 2005-08-17 | 上海交通大学 | 基于核主元分析的非线性故障诊断的方法 |
CN103853152A (zh) * | 2014-03-21 | 2014-06-11 | 北京工业大学 | 一种基于ar-pca的间歇过程故障监测方法 |
CN105700518A (zh) * | 2016-03-10 | 2016-06-22 | 华中科技大学 | 一种工业过程故障诊断方法 |
CN107065843A (zh) * | 2017-06-09 | 2017-08-18 | 东北大学 | 基于独立子空间的多方向kica间歇过程故障监测方法 |
-
2017
- 2017-09-01 CN CN201710777870.3A patent/CN107632592B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1655082A (zh) * | 2005-01-27 | 2005-08-17 | 上海交通大学 | 基于核主元分析的非线性故障诊断的方法 |
CN103853152A (zh) * | 2014-03-21 | 2014-06-11 | 北京工业大学 | 一种基于ar-pca的间歇过程故障监测方法 |
CN105700518A (zh) * | 2016-03-10 | 2016-06-22 | 华中科技大学 | 一种工业过程故障诊断方法 |
CN107065843A (zh) * | 2017-06-09 | 2017-08-18 | 东北大学 | 基于独立子空间的多方向kica间歇过程故障监测方法 |
Non-Patent Citations (2)
Title |
---|
商亮亮 等: ""基于自适应高效递推规范变量分析的多模过程软传感器建模"", 《控制理论与应用》 * |
商亮亮 等: ""基于递推规范变量分析的时变过程故障检测"", 《东北大学学报(自然科学版)》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108764305A (zh) * | 2018-05-14 | 2018-11-06 | 浙江大学 | 一种改进的群智能机器学习故障诊断系统 |
CN108803337A (zh) * | 2018-06-28 | 2018-11-13 | 杭州电子科技大学 | 一种新型的化工过程数据监测方法 |
CN109242046A (zh) * | 2018-10-10 | 2019-01-18 | 中国工程物理研究院计算机应用研究所 | 基于非线性复杂系统过程数据的在线故障诊断方法 |
CN109242046B (zh) * | 2018-10-10 | 2021-11-23 | 中国工程物理研究院计算机应用研究所 | 基于非线性复杂系统过程数据的在线故障诊断方法 |
CN109145256A (zh) * | 2018-11-14 | 2019-01-04 | 南通大学 | 基于规范变量非线性主成分分析的非线性动态过程监测方法 |
CN109145256B (zh) * | 2018-11-14 | 2022-09-16 | 保控(南通)物联科技有限公司 | 一种基于规范变量非线性主成分分析的过程监测方法 |
CN110009126A (zh) * | 2019-01-23 | 2019-07-12 | 北京化工大学 | 基于pls模型与pca贡献度融合的在线报警分析方法 |
CN109885022A (zh) * | 2019-02-21 | 2019-06-14 | 山东科技大学 | 一种基于潜故障敏感子空间的故障检测方法 |
CN109885022B (zh) * | 2019-02-21 | 2021-07-09 | 山东科技大学 | 一种基于潜故障敏感子空间的故障检测方法 |
CN110224160A (zh) * | 2019-06-12 | 2019-09-10 | 北京亿华通科技股份有限公司 | 一种燃料电池系统故障诊断方法 |
CN112329191A (zh) * | 2019-08-05 | 2021-02-05 | 华中科技大学 | 一种多模态过程中模态辨识方法和过渡模态故障检测方法 |
CN110687895A (zh) * | 2019-10-24 | 2020-01-14 | 上海工程技术大学 | 一种基于自适应核主成分分析的化工过程故障检测方法 |
CN110687895B (zh) * | 2019-10-24 | 2022-11-18 | 上海工程技术大学 | 一种基于自适应核主成分分析的化工过程故障检测方法 |
CN110795695A (zh) * | 2019-10-28 | 2020-02-14 | 华北理工大学 | 消去时序相关性的多模型动态pca故障监测方法 |
CN111025222A (zh) * | 2019-12-31 | 2020-04-17 | 三峡大学 | 基于递推主元分析的电子式电压互感器状态监测方法 |
CN111209973A (zh) * | 2020-01-09 | 2020-05-29 | 浙江理工大学 | 基于混合核pca-cca及核密度估计的过程监测方法 |
CN111209973B (zh) * | 2020-01-09 | 2023-04-07 | 浙江理工大学 | 基于混合核pca-cca及核密度估计的过程监测方法 |
CN111914886A (zh) * | 2020-06-13 | 2020-11-10 | 宁波大学 | 一种基于在线简略核学习的非线性化工过程监测方法 |
CN111914886B (zh) * | 2020-06-13 | 2022-07-26 | 宁波大学 | 一种基于在线简略核学习的非线性化工过程监测方法 |
CN111852837B (zh) * | 2020-07-14 | 2022-02-08 | 北京博华信智科技股份有限公司 | 基于聚类的往复式压缩机自适应故障监测方法 |
CN111852837A (zh) * | 2020-07-14 | 2020-10-30 | 北京博华信智科技股份有限公司 | 基于聚类的往复式压缩机自适应故障监测方法 |
CN113253682A (zh) * | 2021-05-17 | 2021-08-13 | 中国石油大学(华东) | 非线性化工过程故障检测方法 |
CN113673884A (zh) * | 2021-08-25 | 2021-11-19 | 国网吉林省电力有限公司长春供电公司 | 基于数字孪生映射模型进行台区数据监测的方法 |
CN114167826A (zh) * | 2021-11-26 | 2022-03-11 | 华中科技大学 | 一种用于铸造生产过程混合多变量的监控方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107632592B (zh) | 2021-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107632592A (zh) | 基于高效递推核主元分析的非线性时变过程故障监测方法 | |
CN108062565B (zh) | 基于化工te过程的双主元-动态核主元分析故障诊断方法 | |
Zhao et al. | Fault-relevant principal component analysis (FPCA) method for multivariate statistical modeling and process monitoring | |
Cai et al. | A new fault detection method for non-Gaussian process based on robust independent component analysis | |
CN109739214B (zh) | 工业过程间歇故障的检测方法 | |
Deng et al. | Modified kernel principal component analysis using double-weighted local outlier factor and its application to nonlinear process monitoring | |
Jiang et al. | Weighted kernel principal component analysis based on probability density estimation and moving window and its application in nonlinear chemical process monitoring | |
CN108958226B (zh) | 基于生存信息势—主成分分析算法的te过程故障检测方法 | |
Zhang et al. | Fault detection in the Tennessee Eastman benchmark process using principal component difference based on k-nearest neighbors | |
Rotem et al. | Ethylene compressor monitoring using model‐based PCA | |
CN101169623A (zh) | 基于核主元分析贡献图的非线性过程故障辨识方法 | |
CN109145256B (zh) | 一种基于规范变量非线性主成分分析的过程监测方法 | |
CN113642754A (zh) | 一种基于rf降噪自编码信息重构和时间卷积网络的复杂工业过程故障预测方法 | |
CN111639304B (zh) | 基于Xgboost回归模型的CSTR故障定位方法 | |
CN110119579B (zh) | 一种基于oica的复杂工业过程故障监测方法 | |
Zhang et al. | Two‐step support vector data description for dynamic, non‐linear, and non‐Gaussian processes monitoring | |
CN112214006A (zh) | 考虑两维动态特性的间歇过程故障检测方法及系统 | |
CN110362063B (zh) | 基于全局保持无监督核极限学习机的故障检测方法及系统 | |
Li et al. | Novel adaptive fault detection method based on kernel entropy component analysis integrating moving window of dissimilarity for nonlinear dynamic processes | |
Linzhe et al. | A nonlinear quality-relevant process monitoring method with kernel input-output canonical variate analysis | |
CN114995338A (zh) | 一种基于规范变量分析与js散度融合的工业过程微小故障检测方法 | |
CN116304823A (zh) | 一种垃圾焚烧过程在线诊断方法 | |
CN116048036A (zh) | 基于动态非平稳投影结构的故障检测方法 | |
CN111983994B (zh) | 一种基于复杂工业化工过程的v-pca故障诊断方法 | |
Cao et al. | Quality-relevant batch process fault detection using a multiway multi-subspace CVA method |
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 |