CN116401577A - 一种基于mcf-occa的质量相关故障检测方法 - Google Patents
一种基于mcf-occa的质量相关故障检测方法 Download PDFInfo
- Publication number
- CN116401577A CN116401577A CN202310323940.3A CN202310323940A CN116401577A CN 116401577 A CN116401577 A CN 116401577A CN 202310323940 A CN202310323940 A CN 202310323940A CN 116401577 A CN116401577 A CN 116401577A
- Authority
- CN
- China
- Prior art keywords
- quality
- train
- variable
- fault detection
- occa
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 69
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims description 26
- 238000012549 training Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000009776 industrial production Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000010219 correlation analysis Methods 0.000 abstract description 7
- 238000010187 selection method Methods 0.000 abstract 1
- 238000004519 manufacturing process Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000001311 chemical methods and process Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开尤其涉及一种基于MCF‑OCCA的质量相关故障检测方法,并在此基础上建立了质量相关的故障检测模型。具体来讲,首先,为了衡量过程变量对质量变量的影响,本发明提出了一种基于互信息和典型相关分析的特征选择方法,计算过程变量对质量变量的影响,选择对质量变量影响较大的过程变量。然后,基于挑选出的变量利用正交典型相关分析建立质量相关故障检测模型。相比于传统的方法,本发明方法剔除了与质量不相关的过程变量,增加了建模的准确性;并构建选择出的过程变量和质量变量的系数矩阵,对其进行SVD分解,得到与质量变量正交和与质量变量相关的子空间,解决了质量变量在线不可测的问题,是一种更优的质量相关故障检测方法。
Description
技术领域
本发明涉及一种数据驱动的故障检测,尤其涉及一种基于MCF-OCCA的质量相关故障检测方法。
背景技术
过程监控和故障检测是维持高质量产品和操作安全的关键。两种常见的故障检测技术是基于模型的和数据驱动的。基于模型的方法通过了解工业过程的工作原理,从而建立机理模型来模拟实际的工业过程。但是工业过程的规模越来越大,过程的非线性,高耦合,非平稳,多模态等特性使得建立机理模型的准确度不高而且很困难。因此人们提出了基于数据驱动的故障检测方法。数据驱动的故障检测技术形式简单,对设计工作量要求低,而且传感器技术的发展,使其在工业应用中得到了广泛的应用。多元统计分析方法是一类典型的数据驱动方法,典型多元统计分析方法主要有主成分分析(Principal ComponentAnalysis,PCA),偏最小二乘法(Projection to Latent Structures,PLS),典型相关分析(Canonical Correlation Analysis,CCA)。PCA是一种无监督的故障检测方法,只能检测过程中发生的故障,无法判别故障是否会影响产品质量。PLS通过最大化过程变量和质量变量的协方差来构建残差进行故障检测,PLS考虑了质量变量和过程变量之间的关系。基于CCA的故障检测方法可以视为基于PCA和基于PLS故障检测方法的一种扩展。CCA提取两组变量之间的最大相关系数,从而构建残差进行检测故障。
质量相关的故障是指过程中发生的能够影响产品质量的故障。在实际的工业过程中,人们往往更关心能够对质量产生影响的故障。利用CCA进行故障检测的原理是利用CCA获得过程变量和质量变量的最大相关系数和投影矩阵,然后构建残差来获得统计量。但是在实际工业过程中质量变量通常是在线不可测的,而且并不是所有的过程变量都会对质量变量产生影响,与质量变量无关的过程变量会降低建模的准确性,增加建模的计算量。因此为了解决以上两个问题,本文提出了一种基于最大相关特征-正交典型相关分析(MaximumCorrelation Feature-Orthogonal Canonical Correlation Analysis,MCF-OCCA)的质量相关故障检测方法。针对质量变量在线不可测的问题,在利用CCA获得过程变量和质量变量的最大相关系数和投影矩阵后,构建过程变量和质量变量的系数矩阵,并对系数矩阵进行SVD分解,将过程变量分为与质量变量正交的子空间和与质量变量相关的子空间,建立统计量对每个子空间进行监测。针对变量冗余问题,在建模前首先利用互信息和CCA对过程变量进行特征提取,选取与质量变量相关的过程变量。
发明内容
本发明所要解决的主要技术问题是:第一是质量变量在线不可测的问题,在利用CCA获得过程变量和质量变量的最大相关系数和投影矩阵后,构建过程变量和质量变量的系数矩阵,并对其进行SVD分解,将过程变量分为与质量变量正交的子空间和与质量变量相关的子空间,对每个子空间进行故障监测。第二是变量冗余问题,在建模前首先利用互信息和CCA对过程变量进行特征提取,选取与质量变量相关的过程变量。
本发明为解决上述问题所采用的技术方案为:一种基于MCF-OCCA的质量相关故障检测方法,包括以下步骤:
(1)采集工业生产对象在正常运行下的样本数据,组成训练数据集,其中过程变量集为质量变量为/>训练数据集的每一列是一个测量点,即特征,每一行是一个样本数据,计算训练数据集/>的均值μ和标准差σ,并按照公式(1)对/>进行标准化,得到标准化后的数据集X∈Rn×m,以同样的方式对/>进行标准化,得到标准化后的数据y∈Rn×1;
其中x为标准化后数据集X的一行数据,x∈R1×m,μ=[μ1,μ2,…,μm],
σ=diag(σ1,σ2,…,σm);
(2)利用CCA和互信息,提取X中质量相关的信息,将其作为训练数据Xtrain∈Rn×v'进行建模,其具体实施过程如下所示:
(3)利用OCCA对Xtrain和y建立模型,得到质量相关投影矩阵U1和质量无关投影矩阵U2的具体实施过程如下所示:
(4)构造统计量T2=xtrainU1inv(cov(XtrainU1))U1 Txtrain T,D2=xtrainU2inv(cov(XtrainU2))U2 Txtrain T。其中inv()为求逆函数,cov()为协方差函数;
上述步骤(1)~(5)为本发明方法的离线建模阶段,如下所示步骤(6)~(9)为本发明方法的在线动态过程监测实施过程;
'
(7)选择与离线建模中Xtrain相同的特征,得到xtrain∈R1×v;
(8)计算在线样本的统计量Tn 2=xntrainU1inv(cov(XtrainU1))U1 Txntrain T,
与传统的方法相比,本发明方法的优势在于:
利用CCA获得过程变量和质量变量的最大相关系数和投影矩阵后,构建过程变量和质量变量的系数矩阵,并对其进行SVD分解,将过程变量分为与质量变量正交的子空间和与质量变量相关的子空间,对每个子空间进行故障检测。解决了质量变量在线不可测的问题。利用互信息和CCA对过程变量进行特征提取,选取与质量变量相关的过程变量。利用一种新的特征提取方法解决了变量冗余的问题。
图1为本发明方法实施流程图。
图2为TE过程工业结构
具体实施方式
下面结合附图于具体的实施案例对本发明方法进行详情对比图。
如图1所示,一种基于MCF-OCCA的质量相关故障检测方法。下面结合一个具体的工业过程的例子来说明本发明的具体实施过程,以及相对于传统基于典型相关分析的故障检测方法的优越性。
应用对象是来自于美国田纳西-伊斯曼(TE)化工过程实验。其具体的工业过程如图2所示,包含了反应器,冷凝器,压缩机,分离器和汽提塔五个主要单元。整个TE过程包含22个测量变量,12个操作变量(其中一个为固定值)和19个成分变量。本案例选择表1所示的22个测量变量和11个操作变量作为过程变量,成分变量36作为质量变量进行建模,接下来结合改TE过程对本发明具体实施步骤进行详细的阐述。
表1:TE过程监控变量。
首先利用TE过程在正常运行状态下采集的960个样本进行离线训练建立故障检测模型,包括以下步骤:
(2)利用CCA和互信息,提取X中质量相关的信息,将其作为训练数据Xtrain∈Rn×v'进行建模,其具体步骤如下:
⑦在向量mi中的33个互信息值中,最大的为MI19,它位于第30列。因此选择的前30列,选择出的特征为7,27,2,3,9,23,15,11,14,33,24,6,1,17,12,25,20,21,29,13,22,32,18,30,10,8,5,26,4,19,所以Xtrain∈R960×30。
(3)利用OCCA对Xtrain和y建立模型,得到质量相关投影矩阵U1和质量无关投影矩阵U2,其具体步骤如下:
(4)构造统计量T2=xtrainU1inv(cov(XtrainU1))U1 Txtrain T,D2=xtrainU2inv(cov(XtrainU2))U2Txtrain T。其中inv()为求逆函数,cov()为协方差函数。
(7)选择7,27,2,3,9,23,15,11,14,33,24,6,1,17,12,25,20,21,29,13,22,32,18,30,10,8,5,26,4,19特征,得到新的训练数据xntrain∈R1×30。
将本方法与传统CCA方法进行对比。
两种方法对TE中质量相关的故障,故障8的检测结果,如表2所示。
表2:对质量相关的故障的检测结果
CCA | 本方法 | |
检测率(%) | 71 | 77 |
两种方法对TE中质量无关的故障,故障14的检测结果,如表3所示。
表3:对质量无关的故障的检测结果
CCA | 本方法 | |
误报率(%) | 15.21 | 0.42 |
通过表2和表3可以看出与传统的CCA相比较,本发明方法能够提升对质量相关故障的故障检测率,降低质量无关故障的误报率。证明了本发明的有效性。
上述实施案例只用来解释说明本发明的具体实施,而不是对本发明进行限制。故凡依本发明之形状,原理所作的变化,均应涵盖在本发明的包含范围内。
Claims (3)
1.一种基于MCF-OCCA的质量相关故障检测方法,其特征在于:所述方法包括以下步骤:
离线建模阶段的实施过程如下所示:
步骤(1)采集工业生产对象在正常运行下的样本数据,组成训练数据集,其中过程变量集为质量变量为/>训练数据集的每一列是一个测量点,即特征,每一行是一个样本数据,计算训练数据集/>的均值μ和标准差σ,并按照公式(1)对/>进行标准化,得到标准化后的数据集X∈Rn×m,以同样的方式对/>进行标准化,得到标准化后的数据y∈Rn×1;
其中x为标准化后数据集X的一行数据,x∈R1×m,μ=[μ1,μ2,…,μm],σ=diag(σ1,σ2,…,σm);
步骤(2)利用CCA和互信息,提取X中质量相关的信息,将其作为训练数据Xtrain∈Rn×v'进行建模;
步骤(3)利用OCCA对Xtrain和y建立模型,得到质量相关投影矩阵U1和质量无关投影矩阵U2;
步骤(4)构造统计量T2=xtrainU1inv(cov(XtrainU1))U1 Txtrain T,D2=xtrainU2inv(cov(XtrainU2))U2 Txtrain T,其中inv()为求逆函数,cov()为协方差函数;
在线监测阶段的实施过程如下所示:
步骤(7)选择与离线建模中Xtrain相同的变量,得到xntrain∈R1×v′;
2.根据权利要求1所述的一种基于MCF-OCCA的质量相关故障检测方法,其特征在于,所述步骤(2)中利用CCA和互信息,提取X中质量相关的信息,将其作为训练数据Xtrain∈Rn×v'进行建模的具体实施过程如下所示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310323940.3A CN116401577A (zh) | 2023-03-30 | 2023-03-30 | 一种基于mcf-occa的质量相关故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310323940.3A CN116401577A (zh) | 2023-03-30 | 2023-03-30 | 一种基于mcf-occa的质量相关故障检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116401577A true CN116401577A (zh) | 2023-07-07 |
Family
ID=87006838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310323940.3A Pending CN116401577A (zh) | 2023-03-30 | 2023-03-30 | 一种基于mcf-occa的质量相关故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116401577A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070112715A1 (en) * | 2005-11-07 | 2007-05-17 | Nec Laboratories America | System failure detection employing supervised and unsupervised monitoring |
CN110501986A (zh) * | 2019-09-03 | 2019-11-26 | 山东科技大学 | 一种基于加权部分互信息的质量相关过程监控方法 |
CN112348088A (zh) * | 2020-11-09 | 2021-02-09 | 北京航空航天大学 | 一种基于互信息典型相关分析的跨模态数据配准方法 |
CN114200914A (zh) * | 2021-12-09 | 2022-03-18 | 华东理工大学 | 一种基于mw-occa的质量相关早期故障检测方法 |
-
2023
- 2023-03-30 CN CN202310323940.3A patent/CN116401577A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070112715A1 (en) * | 2005-11-07 | 2007-05-17 | Nec Laboratories America | System failure detection employing supervised and unsupervised monitoring |
CN110501986A (zh) * | 2019-09-03 | 2019-11-26 | 山东科技大学 | 一种基于加权部分互信息的质量相关过程监控方法 |
CN112348088A (zh) * | 2020-11-09 | 2021-02-09 | 北京航空航天大学 | 一种基于互信息典型相关分析的跨模态数据配准方法 |
CN114200914A (zh) * | 2021-12-09 | 2022-03-18 | 华东理工大学 | 一种基于mw-occa的质量相关早期故障检测方法 |
Non-Patent Citations (2)
Title |
---|
BING SONG等: "Multisubspace Orthogonal Canonical Correlation Analysis for Quality-Related Plant-Wide Process Monitoring", 《IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS》, 7 August 2020 (2020-08-07) * |
金雨婷等: "基于KVAE-OCCA的质量相关故障检测方法及应用", 《控制工程》, 20 February 2022 (2022-02-20), pages 349 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108803520B (zh) | 一种基于变量非线性自相关性剔除的动态过程监测方法 | |
CN108062565B (zh) | 基于化工te过程的双主元-动态核主元分析故障诊断方法 | |
CN112101480B (zh) | 一种多变量聚类与融合的时间序列组合预测方法 | |
CN109522972B (zh) | 一种基于潜变量自回归模型的动态过程监测方法 | |
Zhang et al. | Fault detection and diagnosis of chemical process using enhanced KECA | |
CN110262450B (zh) | 面向汽轮机的多种故障特性协同分析的故障预测方法 | |
CN105955219A (zh) | 基于互信息的分布式动态过程故障检测方法 | |
CN104793606B (zh) | 基于改进的kpca和隐马尔科夫模型的工业故障诊断方法 | |
CN108345284B (zh) | 一种基于两变量块的质量相关故障检测方法 | |
CN109409425B (zh) | 一种基于近邻成分分析的故障类型识别方法 | |
CN108445867B (zh) | 一种基于分散式icr模型的非高斯过程监测方法 | |
CN112904810B (zh) | 基于有效特征选择的流程工业非线性过程监测方法 | |
CN107065842A (zh) | 一种基于粒子群优化核独立元分析模型的故障检测方法 | |
CN108549908A (zh) | 基于多采样概率核主成分模型的化工过程故障检测方法 | |
WO2021114320A1 (zh) | 一种oica和rnn融合模型的污水处理过程故障监测方法 | |
CN108830006B (zh) | 基于线性评价因子的线性-非线性工业过程故障检测方法 | |
CN108181893B (zh) | 一种基于pca-kdr的故障检测方法 | |
CN114757269A (zh) | 一种基于局部子空间-邻域保持嵌入的复杂过程精细化故障检测方法 | |
CN110209145B (zh) | 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法 | |
CN112000081A (zh) | 基于多块信息提取和马氏距离的故障监测方法及系统 | |
CN108572639B (zh) | 一种基于主成分自相关性剔除的动态过程监测方法 | |
CN110705129B (zh) | 一种基于全相关动态kpls的故障诊断方法 | |
CN111915121B (zh) | 一种基于广义典型变量分析的化工过程故障检测方法 | |
CN112947649A (zh) | 一种基于互信息矩阵投影的多变量过程监控方法 | |
CN116305733B (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 |