CN110209145A - 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法 - Google Patents

一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法 Download PDF

Info

Publication number
CN110209145A
CN110209145A CN201910410178.6A CN201910410178A CN110209145A CN 110209145 A CN110209145 A CN 110209145A CN 201910410178 A CN201910410178 A CN 201910410178A CN 110209145 A CN110209145 A CN 110209145A
Authority
CN
China
Prior art keywords
matrix
nuclear matrix
carbon dioxide
absorption tower
dioxide absorption
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
Application number
CN201910410178.6A
Other languages
English (en)
Other versions
CN110209145B (zh
Inventor
宋执环
余家鑫
叶凌箭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201910410178.6A priority Critical patent/CN110209145B/zh
Publication of CN110209145A publication Critical patent/CN110209145A/zh
Application granted granted Critical
Publication of CN110209145B publication Critical patent/CN110209145B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D53/00Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/14Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by absorption
    • B01D53/18Absorbing units; Liquid distributors therefor
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0224Process history based detection method, e.g. whereby history implies the availability of large amounts of data
    • G05B23/024Quantitative history assessment, e.g. mathematical relationships between available data; Functions therefor; Principal component analysis [PCA]; Partial least square [PLS]; Statistical classifiers, e.g. Bayesian networks, linear regression or correlation analysis; Neural networks
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D2257/00Components to be removed
    • B01D2257/50Carbon oxides
    • B01D2257/504Carbon dioxide
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P20/00Technologies relating to chemical industry
    • Y02P20/151Reduction of greenhouse gas [GHG] emissions, e.g. CO2

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Analytical Chemistry (AREA)
  • General Chemical & Material Sciences (AREA)
  • Oil, Petroleum & Natural Gas (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法,属于二氧化碳吸收塔的过程领域。本发明利用基于满秩分解的核矩阵近似方法结合典型相关分析的核学习模型,建立了一个有效的非线性过程二氧化碳吸收塔过程的故障诊断模型,并缓解了实际应用中核方法在数据样本维度较高时计算负荷大的问题,提高了二氧化碳吸收塔过程的故障检测与故障诊断的有效性和准确性,从而保证了现代流程工业的过程安全与可靠性、产品质量以及综合经济效益。

Description

一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法
技术领域
本发明属于二氧化碳吸收塔的过程领域,尤其涉及一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法。
背景技术
在大数据时代下,统计机器学习、数据挖掘等学科的热度高涨,并且这些理论目前在工业过程控制领域已取得了一些进展。在现代流程工业中,为了保证过程安全可靠、产品质量过关、经济效益最大化,故障诊断成为了工业控制系统中不可或缺的部分。对于一个实际的二氧化碳吸收塔,像温度、压力、流量等过程变量与质量变量天然具有较强非线性的数据特性,于是便需要针对测量数据的非线性进行建模。
为了达到上述的目的,二氧化碳吸收塔中常用的故障诊断方法有主成分分析(PCA)、偏最小二乘(PLS)以及它们的扩展方法等,但是这些方法对于过程数据具有的非线性特性的描述能力有限,无法对于具有强非线性的数据进行准确的建模。因此,为了对过程变量的非线性特性进行建模,引入了核方法(Kernel Method),产生了核主成分分析(KPCA)、核偏最小二乘(KPLS)以及它们的一些扩展模型来处理数据的非线性特性。
然而,以上基于核方法的模型具有一定局限性:其核矩阵的维度等于数据样本维度,其计算复杂度会以样本维度的三次方的速度增长,在模型参数的求解中需要对核矩阵进行特征值分解,在实际的应用中其计算复杂度高、计算耗时长,所以传统的核方法的使用收到了限制。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于核矩阵近似的非线性过程二氧化碳吸收塔的故障诊断方法。首先,对采样得到的过程变量数据集进行预处理后对其使用核方法并得到核矩阵;其次,对核矩阵使用典型相关分析进行非线性模型的建立,在求解模型参数时避免对原始样本维度的核矩阵进行特征值分解,而对核矩阵采取矩阵近似,结合矩阵满秩分解后求解近似核矩阵的特征向量,以此减小计算负荷、缩短了运行时间;然后,运用上述核矩阵近似进行典型相关分析建模得到的模型参数进行故障检测;在故障被检测出来后,通过模型参数进行故障诊断,从而实现非线性过程二氧化碳吸收塔的故障诊断。
本发明的目的是通过以下技术方案来实现的:
一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法,包括以下步骤:
(1)采集待检测二氧化碳吸收塔在正常运行时的数据作为训练样本,所述训练样本包括过程变量数据集XNOC,质量变量数据集YNOC,其中XNOC∈RD×M,YNOC∈RD×N,D为训练样本总数,M为过程变量总数,N为质量变量总数;
(2)对过程变量数据集XNOC,质量变量数据集YNOC进行归一化预处理,得到预处理后的训练样本集其中
(3)计算训练样本集的核矩阵KX和KY,其中KX∈RD×D,KY∈RD×D
(4)对核矩阵KX和KY利用典型相关分析进行建模,得到核矩阵K,其中K∈RD×D
(5)对核矩阵K进行矩阵近似得到近似核矩阵其中对近似核矩阵进行特征值分解得到KX和KY的投影矩阵P和Q;
(6)计算训练样本集的残差矩阵E,基于相似度构建SPE统计量,计算控制限SPElimit,得到核矩阵近似的典型相关分析模型;
(7)采集待检测的二氧化碳吸收塔中的数据作为测试样本,对测试样本进行归一化预处理后得到测试样本集
(8)将测试样本集输入步骤(6)所述的核矩阵近似的典型相关分析模型,计算测试样本集的SPE统计量SPEtest,判断是否超过步骤(6)所述控制限SPElimit,若超过,则判断为故障样本,若不超过,则判断为正常样本;
(9)取出所述故障样本,得到故障样本集Ff,从训练样本集中取出正常样本,记为正常样本集Fn,得到待诊断数据集其中Ff∈Rθ×(M+N)
(10)定义标签数据矩阵L作为待诊断数据集X的故障标签,其中
(11)对待诊断数据集X与标签数据矩阵L进行归一化预处理得到通过步骤(3-5)得到用于故障诊断的核矩阵KXd、KLd、近似核矩阵以及投影矩阵U,进一步得到核矩阵KXd的得分矩阵T,其中
(12)由得分矩阵T计算权重W,根据权重W计算待诊断数据集X中的过程变量和质量变量关于的故障关联度r,得到故障诊断结果,其中W∈R(M+N)×η,r∈R(M+N)
优选的,步骤(1)中,利用分布式控制系统采集的正常工况下的数据作为建模使用的训练样本集。对于二氧化碳吸收塔来讲,在一段时间内,采样总数与采样频率有关;变量总数一般由二氧化碳吸收塔的性质与主要影响因素有关。常见的过程变量包括了测量变量与质量变量,过程变量包括但不限于:温度、压力、流量、液位、产物浓度、重量、pH值等。
优选的,步骤(2)中,对测量变量数据集XNOC,质量变量数据集YNOC分别进行归一化预处理,即减去均值然后除以标准差,以消除其在水平和量纲上的差异,使得各个过程变量的均值为0,方差为1。
优选的,步骤(3)选用高斯径向基函数作为核函数计算核矩阵KX和KY
其中为过程变量归一化后的第i次采样,为过程变量归一化后的第j次采样,为质量变量归一化后的第i次采样,为质量变量归一化后的第j次采样,KX∈RD×D,KY∈RD×D,i=1,2,...,D,j=1,2,...,D,σ为核函数的参数。
优选的,步骤(4)具体为:
其中K∈RD×D
优选的,步骤(5)具体为:
(5-1)对核矩阵K采用分解得到:
K=CW+CT
其中C∈RD×Z,W∈RZ×Z,Z为分解中的参数;
(5-2)对W+进行特征值分解得到:
W+=SΛST
其中S∈RZ×V,Λ∈RV×V,W+为W的广义逆矩阵,V为KX的秩;
(5-3)对矩阵C与矩阵S的乘积进行满秩分解得到F,G,其中F∈RD×V,G∈RV×V
(5-4)由G与Λ计算矩阵M:
M=GΛGT
其中M∈RV×V
(5-5)对M进行特征值分解得到:
M=HΩHT
其中H∈RV×V,Ω∈RV×V
(5-6)计算基于核矩阵近似的特征向量矩阵J:
J=FH
其中J∈RD×V
(5-7)基于满秩分解的近似核矩阵的特征值分解为:
其中
(5-8)从J中得到KX和KY的投影矩阵P和Q:
其中P,Q∈RD×V
优选的,步骤(6)具体为:
计算训练样本集的残差矩阵E1,E2
E1=PTKX-ΩQKY
E2=QKYTPTKX
其中
构建SPE统计量:
计算SPE统计量控制限SPElimit
gh=mean(SPE)
2g2h=var(SPE)
其中g,h为卡方分布参数,mean(·)为均值运算符号,var(·)为方差运算符号。
优选的,步骤(8)具体为:
(8-1)选取高斯径向基函数作为核函数,计算测试样本集的核矩阵
其中为测试样本集过程变量归一化后的第i次采样,为训练样本集过程变量归一化后的第j次采样,为测试样本集质量变量归一化后的第i次采样,为训练样本集质量变量归一化后的第j次采样,
(8-2)根据核矩阵以及步骤(5)中所述的投影矩阵P和Q,计算测试样本集的残差:
其中
(8-3)构建SPEtest统计量:
判断SPEtest是否超过步骤(6)所述控制限SPElimit,若超过,则判断为故障样本,若不超过,则判断为正常样本;
优选的,步骤(10)中所述标签数据矩阵L如下所示:
其中其第一列中有个1与θ个0,其第二列中有个0与θ个1。
优选的,步骤(11)所述得分矩阵T的计算公式如下:
T=KXdU
其中η为隐变量的个数。
优选的,步骤(12)具体如下:
(12-1)由得分矩阵T计算权重矩阵W:
其中W∈R(M+N)×η
(12-2)计算待诊断数据集X中的过程变量和质量变量关于的故障关联度r:
其中KMI(·)表示核互信息的计算公式,tα表示得分矩阵T的第α列,w表示权重矩阵W的第i行,第α列的元素;
(12-3)根据故障关联度的大小,得到故障诊断结果。观察故障关联度的大小,得到故障诊断结果,其中故障关联度大的变量为故障变量的可能性更大。
本发明的有益效果是:本发明提出了一种基于满秩分解核矩阵近似的典型相关分析的模型。对于二氧化碳吸收塔过程中带有强非线性关系的过程变量数据,首先进行预处理后对其使用核方法并得到核矩阵,接着对核矩阵使用典型相关分析进行非线性模型的建立,在求解模型参数时避免对原始样本维度的核矩阵进行特征值分解,而对核矩阵采取矩阵近似,结合矩阵满秩分解后求解近似核矩阵的特征向量,然后,运用上述核矩阵近似得到的模型参数进行故障检测;在故障被检测出来后,通过模型参数计算得到故障关联度进行故障诊断,从而实现了非线性过程二氧化碳吸收塔的故障诊断。相比目前其它的二氧化碳吸收塔故障诊断方法,本发明不仅可以进行故障检测,并且能在故障被检测到后进行故障诊断以保证了过程监测的完整性,极大的提高非线性过程二氧化碳吸收塔过程故障诊断的效果,同时很大程度减小计算负荷、缩短了运行时间,提升了基于此发明的故障诊断方法的科学性和有效性。
附图说明
图1为本发明中基于核矩阵近似的二氧化碳吸收塔故障诊断方法的流程图。
具体实施方式
本方法是一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法,该方法针对二氧化碳吸收塔的故障诊断问题,首先利用分布式控制系统收集正常工况下的数据,建立基于满秩分解的核矩阵近似方法结合典型相关分析的核学习模型。核典型相关分析的模型参数由核矩阵的特征值分解得到,其中核矩阵的特征值分解通过其矩阵近似得到。在此基础上,基于该模型构建了故障检测统计量SPE及其对应的统计限SPElimit,对新采样的过程数据进行故障检测。在故障被检测到后,构建待诊断数据集,可再次利用以上模型构建方法,计算每个变量的故障关联度,并得到最终的故障检测结果。
第一步:利用分布式控制系统采集二氧化碳吸收塔正常工况下的数据,组成训练样本数据进行建模,其中测量变量、质量变量数据集分别为:
XNOC∈RD×M,XNOC=[x1,x2,...,xD]T
YNOC∈RD×N,YNOC=[y1,y2,...,yD]L
其中:D为采样总数,M为测量变量总数,N为质量变量总数;
第二步:对测量变量数据集XNOC,质量变量数据集YNOC分别进行数据预处理,以消除其在水平和量纲上的差异,使得各个过程变量的均值为0,方差为1,得到预处理后的训练样本集其中
第三步:选用高斯径向基函数作为核函数计算核矩阵:
其中为过程变量归一化后的第i次采样,为过程变量归一化后的第j次采样,为质量变量归一化后的第i次采样,为质量变量归一化后的第j次采样,KX∈RD×D,KY∈RD×D,i=1,2,...,D,j=1,2,...,D,σ为核函数的参数。
第四步:由核矩阵KX,KY建立核典型相关分析模型,计算核矩阵K:
其中K∈RD×D
第五步:核典型相关分析的模型参数可以通过对K矩阵进行特征值分解得到,但是该矩阵的维度等于样本维度,在样本很大时进行特征值分解的计算消耗量大、耗时长。因此,先对矩阵K采用分解:
K=CW+CT
其中C∈RD×Z,W∈RZ×Z,Z为分解的参数;
对W+进行特征值分解:
W+=SΛST
其中S∈RZ×V,Λ∈RV×V,W+为W的广义逆矩阵,V为KX的秩;
对矩阵C与矩阵S的乘积进行满秩分解得到F,G,其中F∈RD×V,G∈RV×V
计算矩阵M:
M=GΛGT
其中M∈RV×V
对M进行特征值分解:
M=HΩHT
其中H∈RV×V,Ω∈RV×V
计算近似特征向量矩阵:
J=FH
其中J∈RD×V
基于满秩分解的近似核矩阵的特征值分解为:
其中
从J中得到KX,KY的投影矩阵P,Q:
其中P,Q∈RD×V
第六步:通过上一步骤中计算得到的投影矩阵P,Q,计算训练样本集的残差矩阵E1,E2
E1=PTKX-ΩQKY
E2=QKYTPTKX
其中
基于模型的残差矩阵,构建SPE统计量进行故障检测:
SPE统计量控制限SPElimit服从卡方分布其估计方法为:
gh=mean(SPE)
2g2h=var(SPE)
其中g,h为卡方分布参数,mean(·)为均值运算符号,var(·)为方差运算符号。
第七步:利用分布式控制系统采集新的待检测的二氧化碳吸收塔的数据集作为测试样本,其中测量变量数据集Xtest,质量变量数据集Ytest
Xtest∈RD×M,Xtest=[x1,x2,…,xD]
Ytest∈RD×N,Ytest=[y1,y2,…,yD]
其中D为测试样本集的样本数,对其分别进行数据预处理,以消除其在水平和量纲上的差异,使得各个过程变量的均值为0,方差为1,得到预处理后的测试样本集
第八步:采用基于核矩阵近似的核典型相关分析模型对产生测试样本集的二氧化碳吸收塔进行在线故障检测,将测试样本集作为模型输入值,计算测试样本集的SPE统计量SPEtest,判断是否超过第六步所述的控制限SPElimit,若超过,则判断为故障样本,若不超过,则判断为正常样本;得到当前该二氧化碳吸收塔过程的故障检测结果,其具体过程如下:
同样按照第三步中的方法,选取高斯径向基函数作为核函数,由测试样本集与训练样本集,计算测试集的核矩阵:
其中为测试样本集过程变量归一化后的第i次采样,为训练样本集过程变量归一化后的第j次采样,为测试样本集质量变量归一化后的第i次采样,为训练样本集质量变量归一化后的第j次采样,
使用第五步中计算得到的投影矩阵P,Q以及核矩阵计算测试样本集的残差:
其中
基于测试样本集的残差矩阵,构建SPEtest统计量进行故障检测:
若SPEtest超过SPElimit则检测到故障。
第九步:取出所述故障样本,得到故障样本集Ff,将故障样本总数记为S,从训练样本集中取出与故障样本集中的样本数量相同的正常样本,记为正常样本集Fn,得到待诊断数据集其中Ff∈Rθ×(M+N)
第十步:定义标签数据矩阵L作为待诊断数据集X的故障标签:
其中其第一列中有个1与θ个0,其第二列中有个0与θ个1。
第十一步:对待诊断数据集X与标签数据矩阵L分别进行数据预处理,消除其在水平和量纲上的差异,使得各个过程变量的均值为0,方差为1,得到预处理后的待诊断数据集
按照第三步计算用于故障诊断的核矩阵KXd,KLd
其中
由核矩阵KXd,KLd建立核典型相关分析模型,计算核矩阵Kd
其中
核典型相关分析的模型参数可以通过对Kd矩阵进行特征值分解得到,但是该矩阵的维度等于样本维度,在样本很大时进行特征值分解的计算消耗量大、耗时长。选用第五步中的方法,求取基于满秩分解的近似核矩阵的特征值分解为:
其中Ωd∈RV×V
从Jd中得到KXd的投影矩阵U,其中η为隐变量的个数。
计算KXd的得分矩阵:
T=KXdU
其中
第十二步:由得分矩阵T计算权重W:
其中W∈R(M+N)×η
在权重W的指导下,计算待诊断数据集X中的过程变量和质量变量关于的故障关联度r:
其中KMI(·)表示核互信息的计算公式,tα表示得分矩阵T的第α列,w表示权重矩阵W的第i行,第α列的元素。观察故障关联度的大小,得到故障诊断结果,其中故障关联度大的变量为故障变量的可能性越大。
为了检验本发明的可行性,针对训练样本集和测试样本集均为5000个样本,在四种不同故障下,各重复试验5次,使用常规核矩阵特征值分解和使用核矩阵近似的特征值分解对于四种故障能都得到准确的故障诊断结果。为了验证基于满秩分解的核矩阵近似的特征值分解减少的计算损耗,表1记录了使用常规核矩阵特征值分解和使用核矩阵近似的特征值分解的运行时间。
表1基于核矩阵近似的二氧化碳吸收塔故障诊断方法特征值分解耗时比较

Claims (9)

1.一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于,包括以下步骤:
(1)采集待检测二氧化碳吸收塔在正常运行时的数据作为训练样本,所述训练样本包括过程变量数据集XNOC,质量变量数据集YNOC,其中XNOC∈RD×M,YNOC∈RD×N,D为训练样本总数,M为过程变量总数,N为质量变量总数;
(2)对过程变量数据集XNOC,质量变量数据集YNOC进行归一化预处理,得到预处理后的训练样本集
(3)计算训练样本集的核矩阵KX和KY
(4)对核矩阵KX和KY利用典型相关分析进行建模,得到核矩阵K;
(5)对核矩阵K进行矩阵近似得到近似核矩阵对近似核矩阵进行特征值分解得到KX和KY的投影矩阵P和Q;
(6)计算训练样本集的残差矩阵E,基于相似度构建SPE统计量,计算控制限SPElimit,得到核矩阵近似的典型相关分析模型;
(7)采集待检测的二氧化碳吸收塔中的数据作为测试样本,对测试样本进行归一化预处理后得到测试样本集
(8)将测试样本集输入步骤(6)所述的核矩阵近似的典型相关分析模型,计算测试样本集的SPE统计量SPEtest,判断是否超过步骤(6)所述控制限SPElimit,若超过,则判断为故障样本,若不超过,则判断为正常样本;
(9)取出所述故障样本,得到故障样本集Ff,从训练样本集中取出正常样本,记为正常样本集Fn,得到待诊断数据集
(10)定义标签数据矩阵L作为待诊断数据集X的故障标签;
(11)对待诊断数据集X与标签数据矩阵L进行归一化预处理得到通过步骤(3-5)得到用于故障诊断的核矩阵KXd、KLd以及投影矩阵U,进一步得到核矩阵KXd的得分矩阵T;
(12)由得分矩阵T计算权重W,根据权重W计算待诊断数据集X中的过程变量和质量变量关于的故障关联度r,得到故障诊断结果。
2.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于步骤(3)选用高斯径向基函数作为核函数计算核矩阵KX和KY
其中为过程变量归一化后的第i次采样,为过程变量归一化后的第j次采样,为质量变量归一化后的第i次采样,为质量变量归一化后的第j次采样,KX∈RD×D,KY∈RD×D,i=1,2,...,D,j=1,2,...,D,σ为核函数的参数。
3.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于步骤(4)具体为:
其中K∈RD×D
4.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于步骤(5)具体为:
(5-1)对核矩阵K采用分解得到:
K=CW+CT
其中C∈RD×Z,W∈RZ×Z,Z为分解中的参数;
(5-2)对W+进行特征值分解得到:
W+=SΛST
其中S∈RZ×V,Λ∈RV×V,W+为W的广义逆矩阵,V为KX的秩;
(5-3)对矩阵C与矩阵S的乘积进行满秩分解得到F,G,其中F∈RD×V,G∈RV×V
(5-4)由G与Λ计算矩阵M:
M=GΛGT
其中M∈RV×V
(5-5)对M进行特征值分解得到:
M=HΩHT
其中H∈RV×V,Ω∈RV×V
(5-6)计算基于核矩阵近似的特征向量矩阵J:
J=FH
其中J∈RD×V
(5-7)基于满秩分解的近似核矩阵的特征值分解为:
其中
(5-8)从J中得到KX和KY的投影矩阵P和Q:
其中P,Q∈RD×V
5.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于步骤(6)具体为:
计算训练样本集的残差矩阵E1,E2
E1=PTKX-ΩQKY
E2=QKYTPTKX
其中
构建SPE统计量:
计算SPE统计量控制限SPElimit
gh=mean(SPE)
2g2h=var(SPE)
其中g,h为卡方分布参数,mean(·)为均值运算符号,var(·)为方差运算符号。
6.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于步骤(8)具体为:
(8-1)选取高斯径向基函数作为核函数,计算测试样本集的核矩阵
其中为测试样本集过程变量归一化后的第i次采样,为训练样本集过程变量归一化后的第j次采样,为测试样本集质量变量归一化后的第i次采样,为训练样本集质量变量归一化后的第j次采样,
(8-2)根据核矩阵以及步骤(5)中所述的投影矩阵P和Q,计算测试样本集的残差:
其中
(8-3)构建SPEtest统计量:
判断SPEtest是否超过步骤(6)所述控制限SPElimit,若超过,则判断为故障样本,若不超过,则判断为正常样本。
7.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于,步骤(10)中所述标签数据矩阵L如下所示:
其中其第一列中有个1与θ个0,其第二列中有个0与θ个1。
8.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于步骤(11)所述得分矩阵T的计算公式如下:
T=KXdU
其中η为隐变量的个数。
9.根据权利要求1所述的基于核矩阵近似的二氧化碳吸收塔故障诊断方法,其特征在于,步骤(12)具体如下:
(12-1)由得分矩阵T计算权重矩阵W:
其中W∈R(M+N)×η
(12-2)计算待诊断数据集X中的过程变量和质量变量关于的故障关联度r:
其中KMI(·)表示核互信息的计算公式,tα表示得分矩阵T的第α列,w表示权重矩阵W的第i行,第α列的元素;
(12-3)根据故障关联度的大小,得到故障诊断结果。
CN201910410178.6A 2019-05-16 2019-05-16 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法 Active CN110209145B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910410178.6A CN110209145B (zh) 2019-05-16 2019-05-16 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910410178.6A CN110209145B (zh) 2019-05-16 2019-05-16 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法

Publications (2)

Publication Number Publication Date
CN110209145A true CN110209145A (zh) 2019-09-06
CN110209145B CN110209145B (zh) 2020-09-11

Family

ID=67787321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910410178.6A Active CN110209145B (zh) 2019-05-16 2019-05-16 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法

Country Status (1)

Country Link
CN (1) CN110209145B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112101731A (zh) * 2020-08-18 2020-12-18 华南理工大学 一种基于混合核典型相关分析的污水厂在线故障监测方法
CN112286169A (zh) * 2020-10-20 2021-01-29 浙江钱江机器人有限公司 一种工业机器人故障检测方法
CN116907764A (zh) * 2023-09-14 2023-10-20 国能龙源环保有限公司 脱硫设备气密性检测方法、装置、设备及存储介质

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101158873A (zh) * 2007-09-26 2008-04-09 东北大学 一种非线性过程故障诊断方法
CN101308385A (zh) * 2008-07-11 2008-11-19 东北大学 基于二维动态核主元分析的非线性过程故障检测方法
CN102426562A (zh) * 2011-08-15 2012-04-25 天津大学 一种基于支持向量机的核矩阵近似方法
JP2012181607A (ja) * 2011-02-28 2012-09-20 Toyota Central R&D Labs Inc 感覚データ識別装置及びプログラム
CN103903065A (zh) * 2014-04-01 2014-07-02 浙江大学 基于核学习技术的青霉素生产过程产品质量预测方法
CN106930751A (zh) * 2017-04-20 2017-07-07 东北大学 一种有杆泵抽油井故障分离方法
CN106951509A (zh) * 2017-03-17 2017-07-14 中国人民解放军国防科学技术大学 多标签核化典型相关分析检索方法
CN107153748A (zh) * 2017-06-07 2017-09-12 北京信息科技大学 基于加权核主元分析(wkpca)的回转窑故障诊断方法
CN107357275A (zh) * 2017-07-27 2017-11-17 中南大学 非高斯工业过程故障检测方法及系统
CN107463164A (zh) * 2017-07-26 2017-12-12 苏州嘉斯度智能装备有限公司 基于工业大数据多层核学习的连续加热釜快速故障检测法
CN107918379A (zh) * 2017-11-29 2018-04-17 东北大学 基于图半监督代价敏感的工业大数据早期故障检测方法
CN108960329A (zh) * 2018-07-06 2018-12-07 浙江科技学院 一种包含缺失数据的化工过程故障检测方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101158873A (zh) * 2007-09-26 2008-04-09 东北大学 一种非线性过程故障诊断方法
CN101308385A (zh) * 2008-07-11 2008-11-19 东北大学 基于二维动态核主元分析的非线性过程故障检测方法
JP2012181607A (ja) * 2011-02-28 2012-09-20 Toyota Central R&D Labs Inc 感覚データ識別装置及びプログラム
CN102426562A (zh) * 2011-08-15 2012-04-25 天津大学 一种基于支持向量机的核矩阵近似方法
CN103903065A (zh) * 2014-04-01 2014-07-02 浙江大学 基于核学习技术的青霉素生产过程产品质量预测方法
CN106951509A (zh) * 2017-03-17 2017-07-14 中国人民解放军国防科学技术大学 多标签核化典型相关分析检索方法
CN106930751A (zh) * 2017-04-20 2017-07-07 东北大学 一种有杆泵抽油井故障分离方法
CN107153748A (zh) * 2017-06-07 2017-09-12 北京信息科技大学 基于加权核主元分析(wkpca)的回转窑故障诊断方法
CN107463164A (zh) * 2017-07-26 2017-12-12 苏州嘉斯度智能装备有限公司 基于工业大数据多层核学习的连续加热釜快速故障检测法
CN107357275A (zh) * 2017-07-27 2017-11-17 中南大学 非高斯工业过程故障检测方法及系统
CN107918379A (zh) * 2017-11-29 2018-04-17 东北大学 基于图半监督代价敏感的工业大数据早期故障检测方法
CN108960329A (zh) * 2018-07-06 2018-12-07 浙江科技学院 一种包含缺失数据的化工过程故障检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
KAI WANG, JUNGHUI CHEN, ZHIHUAN SONG: "A new excitation scheme for closed-loop subspace identification using additional sampling outputs and its extension to instrumental variable method", 《JOURNAL OF THE FRANKLIN INSTITUTE》 *
付克昌等: "基于特征子空间的KPCA及其在故障检测与诊断中的应用", 《化工学报》 *
舒科等: "基于核典型相关分析的字符识别", 《光电子.激光》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112101731A (zh) * 2020-08-18 2020-12-18 华南理工大学 一种基于混合核典型相关分析的污水厂在线故障监测方法
CN112101731B (zh) * 2020-08-18 2024-03-12 华南理工大学 一种基于混合核典型相关分析的污水厂在线故障监测方法
CN112286169A (zh) * 2020-10-20 2021-01-29 浙江钱江机器人有限公司 一种工业机器人故障检测方法
CN116907764A (zh) * 2023-09-14 2023-10-20 国能龙源环保有限公司 脱硫设备气密性检测方法、装置、设备及存储介质
CN116907764B (zh) * 2023-09-14 2023-12-26 国能龙源环保有限公司 脱硫设备气密性检测方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN110209145B (zh) 2020-09-11

Similar Documents

Publication Publication Date Title
CN109146246B (zh) 一种基于自动编码器和贝叶斯网络的故障检测方法
CN110738274A (zh) 一种基于数据驱动的核动力装置故障诊断方法
CN110209145A (zh) 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法
CN108803520B (zh) 一种基于变量非线性自相关性剔除的动态过程监测方法
CN111780800B (zh) 一种传感器故障的监测、隔离以及重构方法和系统
CN110083860B (zh) 一种基于相关变量选择的工业故障诊断方法
CN107895224A (zh) 一种基于扩展核熵负载矩阵的mkeca发酵过程故障监测方法
CN111368428B (zh) 一种基于监控二阶统计量的传感器精度下降故障检测方法
CN108875108A (zh) 一种基于pca和brb的轨道车辆电路故障诊断方法
CN117290800B (zh) 一种基于超图注意力网络的时序异常检测方法及系统
CN105629109B (zh) 基于art1神经网络的油浸式变压器故障诊断方法
CN111122811A (zh) 一种oica和rnn融合模型的污水处理过程故障监测方法
CN108959498A (zh) 一种用于健康监测的大数据处理平台及其设计方法
CN115455746A (zh) 一种核动力装置运行监测数据异常检测与校正一体化方法
CN107918381A (zh) 一种基于组合核函数的类均值核主元故障诊断方法
CN111367253A (zh) 基于局部自适应标准化的化工系统多工况故障检测方法
CN111188761B (zh) 一种基于Fourier-CVA模型面向机泵设备的监测方法
CN107817784B (zh) 一种基于并发偏最小二乘的过程故障检测方法
CN106250937A (zh) 一种基于非相似度指标的故障分类诊断方法
CN116305733A (zh) 一种基于全局和局部特征提取的质量相关故障检测方法
CN114417704A (zh) 一种基于改进栈式自编码的风电机组发电机健康评估方法
CN111914886B (zh) 一种基于在线简略核学习的非线性化工过程监测方法
CN113298133A (zh) 超临界机组锅炉爆管故障诊断方法
CN115035962A (zh) 基于变分自编码器和生成对抗网络的虚拟样本生成及软测量建模方法
CN105699043A (zh) 一种提高风洞传感器测量稳定性和精准度的方法

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