CN110501986B - 一种基于加权部分互信息的质量相关过程监控方法 - Google Patents
一种基于加权部分互信息的质量相关过程监控方法 Download PDFInfo
- Publication number
- CN110501986B CN110501986B CN201910825708.3A CN201910825708A CN110501986B CN 110501986 B CN110501986 B CN 110501986B CN 201910825708 A CN201910825708 A CN 201910825708A CN 110501986 B CN110501986 B CN 110501986B
- Authority
- CN
- China
- Prior art keywords
- quality
- matrix
- formula
- variable
- krc
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 230000008569 process Effects 0.000 title claims abstract description 50
- 238000012544 monitoring process Methods 0.000 title claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims description 48
- 239000013598 vector Substances 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 4
- 238000012549 training Methods 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 5
- 238000004519 manufacturing process Methods 0.000 abstract description 5
- 238000003745 diagnosis Methods 0.000 abstract description 2
- 230000004927 fusion Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41875—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by quality surveillance of production
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32368—Quality control
-
- 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)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开一种基于加权部分互信息的质量相关过程监控方法,属于故障诊断领域。本发明通过结合部分互信息和贝叶斯加权融合,在多质量指标的监督下,选择与质量密切相关的过程变量,通过相关元分析对所选变量进行建模,获得包含质量变主要信息的相关元。过程数据空间被分为两个子空间,并分别在子空间中构造了统计量用于故障检测,该方法可以提高质量相关故障的检测率,并显著减少误报警,提高实际工业过程的效率;本发明方法可以有效提取隐藏在过程变量和质量变量中的有意义的信息;通过分别构建两个正交子空间中的统计指标,检测与质量相关的故障,并且同时减少与质量无关的干扰报警。
Description
技术领域
本发明属于故障诊断领域,具体涉及一种基于加权部分互信息的质量相关过程监控方法。
背景技术
产品质量与实际工业过程中的生产效率密切相关。实际工业过程中发生的质量相关故障会严重影响产品质量。同时与质量无关的干扰报警的频繁发生,会使得用户对监控系统失去信任,影响正常工况运行。
发明内容
针对现有技术中存在的上述技术问题,本发明提出了一种基于加权部分互信息的质量相关过程监控方法,设计合理,克服了现有技术的不足,具有良好的效果。
为了实现上述目的,本发明采用如下技术方案:
一种基于加权部分互信息的质量相关过程监控方法,包括如下步骤:
步骤1:变量选择;具体包括如下步骤:
对原始过程数据X和质量数据Y分别标准化,如公式(1)和(2)所示:
其中,μx和对角阵Σx分别代表X的均值和标准差,μy和对角阵Σy分别代表Y的均值和标准差;
式中,p(x,y,s)是xi、yi和S的联合概率分布,p(x,y|s)是xi和yi在给定S下的条件概率分布,p(x)和p(y)为边缘概率密度函数,S为过程变量选择集,第一次迭代时过程变量选择集S为空;
步骤1.3:通过公式(4),根据贝叶斯推断,对于每个过程变量xi,计算加权融合后的综合指标BICPMI;
式中,ωj=p(xi|Rj),p(Rj|xi)为xi和yj相关的后验概率,p(xi|Rj)为条件概率;
步骤1.6:选择S中的前kx列以形成新的输入矩阵Xr;其中,kx通过累计质量相关信息百分比CPQI确定,如公式(5)所示:
式中,θxα决定了过程变量的选择比例;
步骤2:离线建模;具体包括如下步骤:
步骤2.1:根据公式(6),计算变量选择后得到的输入矩阵Xr的协方差矩阵C;
步骤2.2:根据公式(7),对协方差矩阵C进行特征值分解;
C=PΛPT (7);
其中,P为投影矩阵,T为得分矩阵,T=XrP;
步骤2.3:根据公式(8),计算得分矩阵T中每个得分向量和每个质量变量之间的部分互信息;
式中,p(t,y,v)是ti、yi和V的联合概率分布,p(t,y|v)是ti和yi在给定V下的条件概率分布,p(t)和p(y)为边缘概率密度函数,V为得分向量选择集,第一次迭代时,得分向量选择集V为空;
步骤2.4:针对每个得分向量,根据公式(9),计算得到加权融合后的综合指标BICPMI:
式中,ωj=p(ti|Rj),p(Rj|ti)为ti和yj相关的后验概率,p(ti|Rj)为条件概率;
步骤2.5:与变量选择相似,返回步骤2.3,继续进行对得分向量的选择和迭代过程,直至得分矩阵T为空;
步骤2.6:选择前krc个与质量最相关的得分向量,得到新的得分矩阵Tr和对应的负载矩阵Pr,krc同样通过累计质量相关信息百分比CPQI确定;
步骤2.7:将输入数据空间Xr分解成两个相互正交的子空间,如公式(10)所示:
式中,Tr为新的得分矩阵,Pr是对应的负载矩阵,Er为残差矩阵;
步骤2.8:根据公式(11)和(12),计算两个子空间相应的控制限;
步骤3:在线故障检测,具体包括如下步骤:
步骤3.1:对新采集的样本x进行标准化;
步骤3.2:根据公式(13)和(14),计算检测统计指标Tr 2和SPEr:
式中,对角阵Λr=diag{λ1,λ2,...,λkrc}代表了相关元的方差;
若:判断结果是检测统计量超过相应控制限,则认为有质量相关的故障发生。
本发明所带来的有益技术效果:
本发明提出一种称为基于加权部分互信息的相关成分分析(WPMI-RCA)的过程监控方法,通过结合部分互信息和贝叶斯加权融合,在多质量指标的监督下,选择与质量密切相关的过程变量,通过相关元分析对所选变量进行建模,获得包含质量变主要信息的相关元。过程数据空间被分为两个子空间,并分别在子空间中构造了统计量用于故障检测,该方法可以提高质量相关故障的检测率,并显著减少误报警,提高实际工业过程的效率;本发明方法可以有效提取隐藏在过程变量和质量变量中的有意义的信息;通过分别构建两个正交子空间中的统计指标,检测与质量相关的故障,并且同时减少与质量无关的干扰报警
附图说明
图1为仿真研究结果图;其中,图(a)为质量y1的变化趋势图;图(b)为质量y2的变化趋势图,图(c)为采用PLS+T2方法的故障检测结果示意图;图(d)为采用PLS+SPE方法的故障检测结果示意图;图(e)为采用WMIPCA+T2方法的故障检测结果示意图;图(f)为采用WMIPCA+SPE方法的故障检测结果示意图;图(g)为采用WPMI-RCA+T2方法的故障检测结果示意图;图(h)为采用WPMI-RCA+SPE方法的故障检测结果示意图。
图2为本发明基于加权部分互信息的质量相关过程监控方法的流程图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明:
一种基于加权部分互信息的质量相关过程监控方法,其流程如图2所示,包括如下步骤:
步骤1:变量选择;具体包括如下步骤:
对原始过程数据X和质量数据Y分别标准化,如公式(1)和(2)所示:
其中,μx和对角阵Σx分别代表X的均值和标准差,μy和对角阵Σy分别代表Y的均值和标准差;
式中,p(x,y,s)是xi、yi和S的联合概率分布,p(x,y|s)是xi和yi在给定S下的条件概率分布,p(x)和p(y)为边缘概率密度函数,S为过程变量选择集,第一次迭代时过程变量选择集S为空;
步骤1.3:通过公式(4),根据贝叶斯推断,对于每个过程变量xi,计算加权融合后的综合指标BICPMI;
式中,ωj=p(xi|Rj),p(Rj|xi)为xi和yj相关的后验概率,p(xi|Rj)为条件概率;
步骤1.6:选择S中的前kx列以形成新的输入矩阵Xr;其中,kx通过累计质量相关信息百分比CPQI确定,如公式(5)所示:
式中,θxα决定了过程变量的选择比例;
步骤2:离线建模;具体包括如下步骤:
步骤2.1:根据公式(6),计算变量选择后得到的输入矩阵Xr的协方差矩阵C;
步骤2.2:根据公式(7),对协方差矩阵C进行特征值分解;
C=PΛPT (7);
其中,P为投影矩阵,T为得分矩阵,T=XrP;
步骤2.3:根据公式(8),计算得分矩阵T中每个得分向量和每个质量变量之间的部分互信息;
式中,p(t,y,v)是ti、yi和V的联合概率分布,p(t,y|v)是ti和yi在给定V下的条件概率分布,p(t)和p(y)为边缘概率密度函数,V为得分向量选择集,第一次迭代时,得分向量选择集V为空;
步骤2.4:针对每个得分向量,根据公式(9),计算得到加权融合后的综合指标BICPMI:
式中,ωj=p(ti|Rj),p(Rj|ti)为ti和yj相关的后验概率,p(ti|Rj)为条件概率;
步骤2.5:与变量选择相似,返回步骤2.3,继续进行对得分向量的选择和迭代过程,直至得分矩阵T为空;
步骤2.6:选择前krc个与质量最相关的得分向量,得到新的得分矩阵Tr和对应的负载矩阵Pr,krc同样通过累计质量相关信息百分比CPQI确定;
步骤2.7:将输入数据空间Xr分解成两个相互正交的子空间,如公式(10)所示:
式中,Tr为新的得分矩阵,Pr是对应的负载矩阵,Er为残差矩阵;
步骤2.8:根据公式(11)和(12),计算两个子空间相应的控制限;
步骤3:在线故障检测,具体包括如下步骤:
步骤3.1:对新采集的样本x进行标准化;
步骤3.2:根据公式(13)和(14),计算检测统计指标Tr 2和SPEr:
式中,对角阵Λr=diag{λ1,λ2,...,λkrc}代表了相关元的方差;
若:判断结果是检测统计量超过相应控制限,则认为有质量相关的故障发生。
仿真研究
仿真实例模型选取如下:
其中,x1,...x8为过程变量,y1,y2为质量变量,ei,vj为标准差为0.001的高斯白噪声。x1,x3和x5是主要的质量相关变量。生成1000个正常样本以组成用于变量选择和建模的训练集,另外生成1000个样本作为在线监测的测试数据集。在第201到第500个样本之间,对质量相关变量x5引入幅值为0.1的阶跃故障和随机噪声;在第401到第700个样本之间,对质量相关变量x3引入幅值为0.1的阶跃故障;在第301到第1000个样本之间,对质量无关变量x4引入幅值为0.1的阶跃故障。仿真结果如图1所示。
从图1容易看出,图(g)和图(h)在第201到第700个采样时刻之间检测到异常,第700个采样时刻之后不再发出报警,这与质量变量的变化是一致的。
本发明所提出的方法在检测质量相关的故障的同时,减少了对质量无关的干扰和故障的报警。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (1)
1.一种基于加权部分互信息的质量相关过程监控方法,其特征在于:包括如下步骤:
步骤1:变量选择;具体包括如下步骤:
对原始过程数据X和质量数据Y分别标准化,如公式(1)和(2)所示:
其中,μx和对角阵Σx分别代表X的均值和标准差,μy和对角阵Σy分别代表Y的均值和标准差;
式中,p(x,y,s)是xi、yi和S的联合概率分布,p(x,y|s)是xi和yi在给定S下的条件概率分布,p(x)和p(y)为边缘概率密度函数,S为过程变量选择集,第一次迭代时过程变量选择集S为空;
步骤1.3:通过公式(4),根据贝叶斯推断,对于每个过程变量xi,计算加权融合后的综合指标BICPMI;
式中,ωj=p(xi|Rj),p(Rj|xi)为xi和yj相关的后验概率,p(xi|Rj)为条件概率;
步骤1.6:选择S中的前kx列以形成新的输入矩阵Xr;其中,kx通过累计质量相关信息百分比CPQI确定,如公式(5)所示:
式中,θxα决定了过程变量的选择比例;
步骤2:离线建模;具体包括如下步骤:
步骤2.1:根据公式(6),计算变量选择后得到的输入矩阵Xr的协方差矩阵C;
步骤2.2:根据公式(7),对协方差矩阵C进行特征值分解;
C=PΛPT (7);
其中,P为投影矩阵,T为得分矩阵,T=XrP;
步骤2.3:根据公式(8),计算得分矩阵T中每个得分向量和每个质量变量之间的部分互信息;
式中,p(t,y,v)是ti、yi和V的联合概率分布,p(t,y|v)是ti和yi在给定V下的条件概率分布,p(t)和p(y)为边缘概率密度函数,V为得分向量选择集,第一次迭代时,得分向量选择集V为空;
步骤2.4:针对每个得分向量,根据公式(9),计算得到加权融合后的综合指标BICPMI:
式中,ωj=p(ti|Rj),p(Rj|ti)为ti和yj相关的后验概率,p(ti|Rj)为条件概率;
步骤2.5:与变量选择相似,返回步骤2.3,继续进行对得分向量的选择和迭代过程,直至得分矩阵T为空;
步骤2.6:选择前krc个与质量最相关的得分向量,得到新的得分矩阵Tr和对应的负载矩阵Pr,krc同样通过累计质量相关信息百分比CPQI确定;
步骤2.7:将输入数据空间Xr分解成两个相互正交的子空间,如公式(10)所示:
式中,Tr为新的得分矩阵,Pr是对应的负载矩阵,Er为残差矩阵;
步骤2.8:根据公式(11)和(12),计算两个子空间相应的控制限;
步骤3:在线故障检测,具体包括如下步骤:
步骤3.1:对新采集的样本x进行标准化;
步骤3.2:根据公式(13)和(14),计算检测统计指标Tr 2和SPEr:
式中,对角阵Λr=diag{λ1,λ2,...,λkrc}代表了相关元的方差;
若:判断结果是检测统计量超过相应控制限,则认为有质量相关的故障发生。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910825708.3A CN110501986B (zh) | 2019-09-03 | 2019-09-03 | 一种基于加权部分互信息的质量相关过程监控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910825708.3A CN110501986B (zh) | 2019-09-03 | 2019-09-03 | 一种基于加权部分互信息的质量相关过程监控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110501986A CN110501986A (zh) | 2019-11-26 |
CN110501986B true CN110501986B (zh) | 2021-06-18 |
Family
ID=68591055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910825708.3A Active CN110501986B (zh) | 2019-09-03 | 2019-09-03 | 一种基于加权部分互信息的质量相关过程监控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110501986B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116401577B (zh) * | 2023-03-30 | 2024-06-21 | 华东理工大学 | 一种基于mcf-occa的质量相关故障检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104239536A (zh) * | 2014-09-22 | 2014-12-24 | 广西教育学院 | 基于互信息的完全加权课程正负关联模式挖掘方法及系统 |
CN105955219A (zh) * | 2016-05-30 | 2016-09-21 | 宁波大学 | 基于互信息的分布式动态过程故障检测方法 |
CN109407652A (zh) * | 2018-12-10 | 2019-03-01 | 中国石油大学(华东) | 基于主辅pca模型的多变量工业过程故障检测方法 |
-
2019
- 2019-09-03 CN CN201910825708.3A patent/CN110501986B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104239536A (zh) * | 2014-09-22 | 2014-12-24 | 广西教育学院 | 基于互信息的完全加权课程正负关联模式挖掘方法及系统 |
CN105955219A (zh) * | 2016-05-30 | 2016-09-21 | 宁波大学 | 基于互信息的分布式动态过程故障检测方法 |
CN109407652A (zh) * | 2018-12-10 | 2019-03-01 | 中国石油大学(华东) | 基于主辅pca模型的多变量工业过程故障检测方法 |
Non-Patent Citations (4)
Title |
---|
"Batch process monitoring based on just-in-time learning and multiple-subspace principal component analysis";Zhaomin Lv, etc.;《Chemometrics and Intelligent Laboratory Systems》;20141015(第137期);128-139 * |
"基于加权统计局部核主元分析的非线性化工过程微小故障诊断方法";邓佳伟 等;《化工学报》;20190731;第70卷(第7期);2594-2605 * |
"基于部分互信息和贝叶斯打分函数的基因调控网络构建算法";刘飞 等;《西北工业大学学报》;20171031;第35卷(第5期);876-883 * |
"面向互联网的局域底层软件运行故障检测仿真";高文莲 等;《计算机仿真》;20180731;第35卷(第7期);305-309 * |
Also Published As
Publication number | Publication date |
---|---|
CN110501986A (zh) | 2019-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104777830B (zh) | 一种基于kpca混合模型的多工况过程监控方法 | |
CN109407652B (zh) | 基于主辅pca模型的多变量工业过程故障检测方法 | |
CN105425779A (zh) | 基于局部邻域标准化和贝叶斯推断的ica-pca多工况故障诊断方法 | |
CN109189028A (zh) | 基于多块信息提取的pca故障诊断方法 | |
CN109507972B (zh) | 基于分层式非高斯监测算法的工业生产过程故障监测方法 | |
CN108536128A (zh) | 一种参数优化的机器学习故障诊断系统 | |
Wang et al. | Fault detection and identification using a Kullback-Leibler divergence based multi-block principal component analysis and Bayesian inference | |
Zheng et al. | Quality‐relevant dynamic process monitoring based on mutual information multiblock slow feature analysis | |
Huang et al. | Adaptive process monitoring via online dictionary learning and its industrial application | |
CN109298633A (zh) | 基于自适应分块非负矩阵分解的化工生产过程故障监测方法 | |
CN111736567B (zh) | 一种基于故障敏感慢特征的多块故障监测方法 | |
CN110501986B (zh) | 一种基于加权部分互信息的质量相关过程监控方法 | |
CN111639304A (zh) | 基于Xgboost回归模型的CSTR故障定位方法 | |
Qian et al. | Multichannel profile-based monitoring method and its application in the basic oxygen furnace steelmaking process | |
Kämpjärvi et al. | Fault detection and isolation of an on-line analyzer for an ethylene cracking process | |
CN112149054B (zh) | 基于时序扩展的正交邻域保持嵌入模型的构建与应用 | |
Zhang et al. | Decentralized dynamic monitoring based on multi-block reorganized subspace integrated with Bayesian inference for plant-wide process | |
CN106647274B (zh) | 一种连续生产过程中运行工况稳态判别方法 | |
Liu et al. | Structured sequential Gaussian graphical models for monitoring time-varying process | |
CN112631258A (zh) | 一种工业过程关键指标的故障预警方法 | |
CN116305733A (zh) | 一种基于全局和局部特征提取的质量相关故障检测方法 | |
CN110347138A (zh) | 自适应核主元分析的化工过程故障诊断方法 | |
CN114879612A (zh) | 一种基于Local-DBKSSA的高炉炼铁过程监测方法 | |
Zhong et al. | A Quality-related Fault Detection Method Based on Weighted Mutual Information | |
CN115407753A (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 |