CN108958226A - 基于生存信息势—主成分分析算法的te过程故障检测方法 - Google Patents
基于生存信息势—主成分分析算法的te过程故障检测方法 Download PDFInfo
- Publication number
- CN108958226A CN108958226A CN201810899623.5A CN201810899623A CN108958226A CN 108958226 A CN108958226 A CN 108958226A CN 201810899623 A CN201810899623 A CN 201810899623A CN 108958226 A CN108958226 A CN 108958226A
- Authority
- CN
- China
- Prior art keywords
- data
- fault
- spe
- sip
- pca
- 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
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 29
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 11
- 238000012360 testing method Methods 0.000 title abstract description 5
- 230000008569 process Effects 0.000 claims abstract description 25
- 238000001514 detection method Methods 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000006870 function Effects 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 13
- 230000004083 survival effect Effects 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000002790 cross-validation Methods 0.000 claims description 4
- 238000011478 gradient descent method Methods 0.000 claims description 4
- 238000013499 data model Methods 0.000 claims description 3
- 238000013179 statistical model Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 abstract description 8
- 238000009826 distribution Methods 0.000 description 14
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 5
- 239000000047 product Substances 0.000 description 5
- 238000001311 chemical methods and process Methods 0.000 description 3
- 239000003054 catalyst Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000000376 reactant Substances 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 239000006227 byproduct Substances 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000012824 chemical production Methods 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010978 in-process monitoring Methods 0.000 description 1
- 239000011261 inert gas Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000012263 liquid product Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000011112 process operation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric 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/0224—Process history based detection method, e.g. whereby history implies the availability of large amounts of data
- G05B23/024—Quantitative 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Debugging And Monitoring (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及非高斯随机系统故障诊断领域,具体为基于生存信息势—主成分分析算法的TE过程故障检测方法,解决TE过程不服从高斯分布且变量个数较多,传统的方法结果不精确的问题,步骤:一、采集正常数据归一化处理;二、计算平均重构误差;三、构建SIP‑PCA算法;四、估计SPE及T2统计量的控制限;五、采集故障数据与控制限对比,检测是否发生故障。优点:计算复杂度低,迭代不涉及指数运算,收敛速度快;无分布假设,更具一般性;控制限的计算过程中采用核密度估计法,充分考虑非高斯噪声的影响,提高了故障检测的精度;使用SIP使变量的随机性变小,在非高斯噪声干扰下的估计结果具有更好的平滑性。
Description
技术领域
本发明涉及非高斯随机系统故障诊断领域,具体为基于生存信息势(SurvivalInformation Potential,SIP)—主成分分析(Principal Component Analysis,PCA)算法的田纳西—伊斯曼(Tennessee Eastman,TE)过程故障检测方法。
背景技术
化工过程具有技术密集、工艺复杂、生产过程联系紧密等特点。化工过程如果轻视安全的重要性,就可能造成重大事故,这样不仅会影响生产而且会造成巨大的经济损失甚至是人员的伤亡,如果大量的有毒有害的物质发生泄漏将还会对环境造成重大的污染。为了确保生产过程运行正常,降低设备的维修率,保证产品的生产质量,需要对生产过程进行全程监测,所以故障诊断成为有实际意义的研究课题。
随着故障检测和诊断在内的多变量统计方法得到了迅速发展。许多于从海量数据集中提取过程信息并解释这些信息的统计技术已在各个领域得到发展。这种方法不需要准确的经验知识也不需要复杂数学模型的建立,而且釆用的数据是工业过程中的第一手资料,因此更接近真实情况。附加成本低,假定前提少,且极易维护,因此适用于复杂的工业过程。已证明其中特别有效的一种方法是使用PCA结合T2图表和SPE图表。
但是传统的基于PCA的监测都假设数据服从高斯分布,实际上,并非所有的过程变量都是高斯分布的,其中一些可能遵循不同类型的非高斯分布,例如TE过程,该过程是依据实际化工联合反应过程开发的模型仿真平台,这个过程模型包括42个测量变量和12个控制变量,其产生的数据具体时变、强耦合、非线性特征,且数据分布不服从高斯分布。在这种情况下,传统的基于PCA的监控方法可能无法很好地运行,可能会导致过程监控的不准确。因此,为了解决以上这些问题,本发明提出一种基于生存信息势—主成分分析算法的TE过程故障检测方法。
发明内容
本发明旨在解决TE过程不服从高斯分布且变量个数较多,采取传统的主成分分析方法会出现监测结果不精确的问题,提出了一种基于生存信息势—主成分分析算法的TE过程故障检测方法。
本发明是通过以下操作步骤实现的:基于生存信息势—主成分分析算法的TE过程故障检测方法,包括以下操作步骤:
一、采集正常工况下的数据X∈RN×m(x1,x2,...,xm)作为正常数据,并且对其进行归一化处理,其中N表示样本数,m表示变量个数;
二、从重构误差的角度考虑,PCA就是寻找一个投影矩阵,使得平均重构误差最小,平均重构误差的计算公式如下:
即
其中xi为任一样本点,为样本的投影点;ti=PTxi是m×1维的得分向量,即为主成分向量;P为加载矩阵,即为主成分的投影方向;
三、基于SIP—PCA的算法:
1.构建SIP—PCA算法:SIP的本质是在传统的信息势表达式中,用生存函数S(x)=P{X>x}=1-F(x)取代概率密度函数;由此,对于一个向量其α-阶(α>0)生存信息势的定义为:
其中,是随机向量|X|的多元生存函数;|X|代表随机向量X的绝对值,是一个包含元素|X1|,|X2|,K,|Xm|的m维随机向量;
为了计算简便,引入了经验SIP,假设有N个样本X(1),X(2),...,X(N),在每个样本点的经验生存函数为:
因此,经验SIP可以表示为:
当样本数据为标量数据时,即m=1,对于N个样本,0≤X(1)≤X(2)≤…≤X(N),可以得到:
假设X(0)=0,可以得到一种更一般的形式:
其中:
由于式7在x(i)=0处不平滑,故使用有序样本的平方0≤X2(1)≤X2(2)≤…≤X2(N)来计算经验SIP,即:
然后,我们将式2带入式9,得到SIP-PCA的性能指标:
其中,P要求为正交矩阵,最终得出:
其中,tj=PTxj,P=[p1,p2,...,pm];
2.采用梯度下降法对性能指标进行优化,得到所有方向的方向矩阵P:
PCA是一种统计投影技术,其本质是坐标变换,由上述建立SIP—PCA算法,当重构误差最小时说明已经在保证数据不失真的情况下,将原始数据投影到了新的坐标系中;
SIP是衡量估计结果随机性的度量,即当SIP最小时,原始数据X与重构数据之间的误差最小;
因此,通过优化性能指标即可得到最优的投影矩阵P;由于性能指标式11是带约束的优化问题,本发明使用拉格朗日因子法,将上述带约束问题转换成无约束优化问题,如式12所示:
然后,采用梯度下降法对该性能指标进行优化,主要分为以下两步:
b1、对加载矩阵P和拉格朗日因子求偏导:
b2、根据递推公式求得未知参数P和
其中μ和η为学习率,且μ∈[0,1],η∈[0,1];其存在是为了均衡偏导部分在幅值上对所求最小值的影响;
3.使用交叉验证法确定PCA的主元个数:
使用PCA构造主元模型时,必须确定主元个数,而主元个数的确定需要考虑两个方面的因素:原始数据维数的降低和原始测量数据信息的丢失。主元个数的选择会直接影响故障诊断的结果,可见,主元个数的选择十分重要。常用的主元个数的选择方法有累计方差贡献率法、未重构方差法、交叉检验法等。因为前两种方法都需要人为设置一个阈值来确定要选取得主元个数,所以本发明采用交叉验证方法求取主元,具体步骤如下:
c1、将步骤一采集的原始数据X∈RN×m分为两部分,去掉第i行得到X(-i)∈R(N-1)×m;以及只有第i行的x(i)∈R1×m;
c2、使用数据X(-i)∈R(N-1)×m执行SIP-PCA算法,得到加载矩阵P(-i);
c3、将数据x(i)∈R1×m投影到加载矩阵P(-i)上求得其得分:
t(i)=x(i)P(-i) 式15
c4、确定x(i)∈R1×m的残差:
重复执行c1-c4,直至得到所有的残差;
c5、计算预测残差平方和:
PRESS最小值所对应的j即为所求的最佳主元个数;
四、将置信度设置为α,使用核密度估计法来估计SPE统计量以及T2统计量的概率密度函数,通过对两个统计量的概率密度函数积分得到故障检测控制限,包括如下步骤:
1.T2统计量表示得分向量的平方和,其具体计算为:
T2=xTPPTx 式18
2.SPE统计量表示每个采样数据与统计模型之间的误差,具体的计算为:
SPE=||(I-PPT)x||2 式19
3.由于数据是非高斯的,仍使用传统基于高斯假设的控制限的求法,可能会导致结果不准确,因此使用核密度估计法估计SPE和T2的概率密度函数:
4.根据置信度α,对T2和SPE统计量的概率密度函数分别求积分,得到对应的控制限和CLSPE:
五、采集故障工况下的数据作为故障数据,对其进行归一化处理,随后将故障数据投影到加载矩阵上,计算得到故障数据模型的SPE统计量以及T2统计量,通过与式22和式23计算出的控制限作对比,检测系统是否发生了故障。
与现有的技术相比本发明具有以下几个优势:实际化工过程中,并非所有的过程变量都是高斯分布的,其中一些可能遵循不同类型的非高斯分布。本方法使用了一种新的衡量估计结果随机性的性能指标—SIP,其本质是在传统的信息势表达式中,用生存函数取代概率密度函数。SIP不涉及核计算,使得其计算复杂度低,每次迭代不涉及指数运算,收敛速度快。SIP没有任何的分布假设,更具一般性;与传统的求取控制限的方法不同,本发明充分考虑到非高斯噪声的影响,在计算控制限的过程中采用核密度估计法,有效解决了传统PCA算法基于高斯分布的假设,提高了故障检测的精度。此外,使用SIP不仅可以使变量的随机性变小,还能使误差接近于零,且在非高斯噪声干扰下的估计结果具有相对较好的平滑性。
附图说明
图1为本发明SIP—PCA算法的故障诊断结构框图;
图2a为其中一个变量的正常数据分布图
图2b为图2a变量加入故障后的数据分布图;由图2b与图2a对比可以看出,未出现故障时,数据的分布范围为[-4,4],当给定一个阶跃信号作为故障数据时,数据的分布范围变成了[0,8],偏离了正常范围;
图3a为T2统计量的概率密度函数图;
图3b为SPE统计量的概率密度函数图;已知曲线与X轴所围成的面积等于1,这里通过设置统计阈置信度α来具体确定置信限:例如,此时置信度设为0.99,那么对应的控制限就是一个T2与SPE的值,在该值右边曲线与X轴围城的面积占总面积的0.01,即占总面积的α;此时,便求出了SIP-PCA的控制限与CLSPE;
图4a为传统PCA算法的主元分析统计量变化图(故障检测结果图);
图4b为SIP-PCA算法的主元分析统计量变化图(故障检测结果图);由图4a与图4b对比可知,在250秒时人为加入故障,传统的PCA算法不能及时检测出故障,而本申请提出的SIP-PCA算法则在检测故障的快速性和准确性方面均显著优于传统PCA算法。
具体实施方式
以下结合具体实施例对本发明作进一步说明:在某化工生产厂的整个TE过程主要有五个操作单元组成:反应器、产品冷凝器、气液分离器、循环压缩机和汽提塔。气态的反应物进入到反应器中,生成液态产品气相的反应是在一种不挥发的气相催化剂的作用下进行的。反应器内置有冷凝包用来移除反应产生的热量。产品以气态的形式出来,并夹杂有一些未反应物。催化剂仍然滞留在反应器中。
该TE生产过程主要由四种气态物料参与反应,分别为A、C、D和E,生产出两种产品G、H,并伴有一种副产品F,此外在产品的进料中含有少量的惰性气体B。整个过程主要由四种反应组成,反应方程式如下:
此TE过程主要有12个控制变量和41个测量变量,在确定过程监测变量时,通常选取常用的与过程运行状况紧密相关的16个过程变量,如表1所示;作为正常数据对系统进行故障诊断,数据集分为两部分:训练数据和测试数据,其中采样时间设置为1s,训练样本个数为500,测试数据在250s引入故障。
表1 TE过程检测变量列表
表2所示为TE过程的故障列表,其列举了多种类型的故障,可以体现实际应用中可能会出现的各种故障。
表2 TE过程故障列表
该实施例采用基于生存信息势—主成分分析算法的TE过程故障检测方法,包括以下操作步骤:
一、采集正常工况下的数据X∈RN×m(x1,x2,...,xm)作为正常数据,并且对其进行归一化处理,其中N表示的样本数在本实施例中为16,m表示的变量个数为500;
二、计算平均重构误差:
即
其中xi为任一样本点,为样本的投影点;ti=PTxi是m×1维的得分向量,即为主成分向量;P为加载矩阵,即为主成分的投影方向;
三、基于SIP—PCA的算法:
1.构建SIP—PCA算法:
2.采用梯度下降法对性能指标进行优化,得到所有方向的方向矩阵P:
b1、对加载矩阵P和拉格朗日因子求偏导:
b2、根据递推公式求得未知参数P和
其中μ和η为学习率,且μ∈[0,1],η∈[0,1];
3.使用交叉验证法确定PCA的主元个数,num(num<16),则原始数据的维度降低X∈R500×16→X∈R500×num,具体步骤如下:
c1、将步骤一采集的原始数据X∈RN×m分为两部分,去掉第i行得到X(-i)∈R(N-1)×m;以及只有第i行的x(i)∈R1×m;
c2、使用数据X(-i)∈R(N-1)×m执行SIP-PCA算法,得到加载矩阵P(-i);
c3、将数据x(i)∈R1×m投影到加载矩阵P(-i)上求得其得分:t(i)=x(i)P(-i) 式15
c4、确定x(i)∈R1×m的残差:
重复执行c1-c4,直至得到所有的残差;
c5、计算预测残差平方和:
PRESS最小值所对应的j即为所求的最佳主元个数;至此,完成了SIP-PCA算法的构建;
四、将置信度设置为α,使用核密度估计法来估计SPE统计量以及T2统计量的概率密度函数,通过对两个统计量的概率密度函数积分得到故障检测控制限,包括如下步骤:
1.T2统计量表示得分向量的平方和,其具体计算为:
T2=xTPPTx 式18
2.SPE统计量表示每个采样数据与统计模型之间的误差,具体的计算为:
SPE=||(I-PPT)x||2 式19
3.使用核密度估计法估计SPE和T2的概率密度函数:
4.根据置信度α,对T2和SPE统计量的概率密度函数分别求积分,得到对应的控制限和CLSPE:
五、采集故障工况下的数据作为故障数据,对其进行归一化处理,随后将故障数据投影到加载矩阵上,计算得到故障数据模型的SPE统计量以及T2统计量,通过与式22和式23计算出的控制限作对比,检测系统是否发生了故障。
例如,给定一个数据表达式如下:
X=UY+2E 式25
其中X∈R500×7;Y是利用高斯混合模型得到的非高斯分布的数据;E是服从Beta分布的噪声,在250s时刻,给数据X的第二个变量一个阶跃信号作为故障数据,其故障诊断结果如图4b所示,从结果可以看到,在250s时,T2和SPE统计量的值均超过了控制限,说明系统出现了故障。
Claims (1)
1.一种基于生存信息势-主成分分析算法的TE过程故障检测方法,其特征在于:包括以下操作步骤:
一、采集正常工况下的数据X∈RN×m(x1,x2,...,xm)作为正常数据,并且对其进行归一化处理,其中N表示样本数,m表示变量个数;
二、计算平均重构误差:
即
其中xi为任一样本点,为样本的投影点;ti=PTxi是m×1维的分向量,即为主成分向量;P为加载矩阵,即为主成分的投影方向;
三、基于SIP-PCA的算法:
1.构建SIP-PCA算法:
2.采用梯度下降法对性能指标进行优化,得到所有方向的方向矩阵P:
b1、对加载矩阵P和拉格朗日因子求偏导:
b2、根据递推公式求得未知参数P和
其中μ和η为学习率,且μ∈[0,1],η∈[0,1];
3.使用交叉验证法确定PCA的主元个数,具体步骤如下:
c1、将步骤一采集的原始数据X∈RN×m分为两部分,去掉第i行得到X(-i)∈R(N-1)×m;以及只有第i行的x(i)∈R1×m:
c2、使用数据X(-i)∈R(N-1)×m执行SIP-PCA算法,得到加载矩阵P(-i);
c3、将数据x(i)∈R1×m投影到加载矩阵P(-i)上求得其得分:
t(i)=x(i)P(-i)
c4、确定x(i)∈R1×m的残差:
重复执行c1-c4,直至得到所有的残差;
c5、计算预测残差平方和:
PRESS最小值所对应的j即为所求的最佳主元个数;至此,完成了SIP-PCA算法的构建;
四、将置信度设置为α,使用核密度估计法来估计SPE统计量以及T2统计量的概率密度函数,通过对两个统计量的概率密度函数积分得到故障检测控制限,包括如下步骤:
1.T2统计量表示得分向量的平方和,其具体计算为:
T2=xTPPTx
2.SPE统计量表示每个采样数据与统计模型之间的误差,具体的计算为:
SPE=||(I-PPT)x||2
3.使用核密度估计法估计SPE和T2的概率密度函数:
4.根据置信度α,对T2和SPE统计量的概率密度函数分别求积分,得到对应的控制限和CLSPE:
五、采集故障工况下的数据作为故障数据,对其进行归一化处理,随后将故障数据投影到加载矩阵上,计算得到故障数据模型的SPE统计量以及T2统计量,通过与步骤四、4计算出的控制限作对比,检测系统是否发生了故障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810899623.5A CN108958226B (zh) | 2018-08-08 | 2018-08-08 | 基于生存信息势—主成分分析算法的te过程故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810899623.5A CN108958226B (zh) | 2018-08-08 | 2018-08-08 | 基于生存信息势—主成分分析算法的te过程故障检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108958226A true CN108958226A (zh) | 2018-12-07 |
CN108958226B CN108958226B (zh) | 2021-03-19 |
Family
ID=64468859
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810899623.5A Active CN108958226B (zh) | 2018-08-08 | 2018-08-08 | 基于生存信息势—主成分分析算法的te过程故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108958226B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507972A (zh) * | 2018-12-19 | 2019-03-22 | 中国计量大学 | 基于分层式非高斯监测算法的工业生产过程故障监测方法 |
CN110009126A (zh) * | 2019-01-23 | 2019-07-12 | 北京化工大学 | 基于pls模型与pca贡献度融合的在线报警分析方法 |
CN110598370A (zh) * | 2019-10-18 | 2019-12-20 | 太原理工大学 | 基于sip和ekf融合的多旋翼无人机鲁棒姿态估计 |
CN110942258A (zh) * | 2019-12-10 | 2020-03-31 | 山东科技大学 | 一种性能驱动的工业过程异常监测方法 |
CN112560165A (zh) * | 2020-06-11 | 2021-03-26 | 中车青岛四方机车车辆股份有限公司 | 一种城轨车辆及其客室车门故障诊断方法 |
CN112836431A (zh) * | 2021-02-05 | 2021-05-25 | 河北工业大学 | 基于pso-lstm的青霉素发酵过程故障预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001123335A (ja) * | 1999-10-21 | 2001-05-08 | Nippon Ester Co Ltd | 分割型ポリエステル複合繊維 |
CN105182955A (zh) * | 2015-05-15 | 2015-12-23 | 中国石油大学(华东) | 一种多变量工业过程故障识别方法 |
CN105700518A (zh) * | 2016-03-10 | 2016-06-22 | 华中科技大学 | 一种工业过程故障诊断方法 |
CN107092242A (zh) * | 2017-06-02 | 2017-08-25 | 宁波大学 | 一种基于缺失变量pca模型的工业过程监测方法 |
CN108153267A (zh) * | 2017-12-15 | 2018-06-12 | 宁波大学 | 一种基于误差主元分析模型的工业过程监测方法 |
-
2018
- 2018-08-08 CN CN201810899623.5A patent/CN108958226B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001123335A (ja) * | 1999-10-21 | 2001-05-08 | Nippon Ester Co Ltd | 分割型ポリエステル複合繊維 |
CN105182955A (zh) * | 2015-05-15 | 2015-12-23 | 中国石油大学(华东) | 一种多变量工业过程故障识别方法 |
CN105700518A (zh) * | 2016-03-10 | 2016-06-22 | 华中科技大学 | 一种工业过程故障诊断方法 |
CN107092242A (zh) * | 2017-06-02 | 2017-08-25 | 宁波大学 | 一种基于缺失变量pca模型的工业过程监测方法 |
CN108153267A (zh) * | 2017-12-15 | 2018-06-12 | 宁波大学 | 一种基于误差主元分析模型的工业过程监测方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507972A (zh) * | 2018-12-19 | 2019-03-22 | 中国计量大学 | 基于分层式非高斯监测算法的工业生产过程故障监测方法 |
CN110009126A (zh) * | 2019-01-23 | 2019-07-12 | 北京化工大学 | 基于pls模型与pca贡献度融合的在线报警分析方法 |
CN110598370A (zh) * | 2019-10-18 | 2019-12-20 | 太原理工大学 | 基于sip和ekf融合的多旋翼无人机鲁棒姿态估计 |
CN110598370B (zh) * | 2019-10-18 | 2023-04-14 | 太原理工大学 | 基于sip和ekf融合的多旋翼无人机鲁棒姿态估计 |
CN110942258A (zh) * | 2019-12-10 | 2020-03-31 | 山东科技大学 | 一种性能驱动的工业过程异常监测方法 |
CN110942258B (zh) * | 2019-12-10 | 2022-02-25 | 山东科技大学 | 一种性能驱动的工业过程异常监测方法 |
CN112560165A (zh) * | 2020-06-11 | 2021-03-26 | 中车青岛四方机车车辆股份有限公司 | 一种城轨车辆及其客室车门故障诊断方法 |
CN112836431A (zh) * | 2021-02-05 | 2021-05-25 | 河北工业大学 | 基于pso-lstm的青霉素发酵过程故障预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108958226B (zh) | 2021-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108958226B (zh) | 基于生存信息势—主成分分析算法的te过程故障检测方法 | |
Jia et al. | Quality-related fault detection approach based on dynamic kernel partial least squares | |
CN108062565B (zh) | 基于化工te过程的双主元-动态核主元分析故障诊断方法 | |
CN107632592B (zh) | 基于高效递推核主元分析的非线性时变过程故障监测方法 | |
KR101360790B1 (ko) | 발전소 계측기 성능 감시 예측 방법 | |
Van Sprang et al. | Critical evaluation of approaches for on-line batch process monitoring | |
CN108549908B (zh) | 基于多采样概率核主成分模型的化工过程故障检测方法 | |
KR101178235B1 (ko) | 에프에스브이알과 지엘알티를 이용한 발전소 계측기 성능감시 방법 및 시스템 | |
Chen et al. | Probabilistic contribution analysis for statistical process monitoring: A missing variable approach | |
CN109739214A (zh) | 工业过程间歇故障的检测方法 | |
CN105676833A (zh) | 发电过程控制系统故障检测方法 | |
US20090076773A1 (en) | Method for identifying unmeasured disturbances in process control test data | |
CN111368428B (zh) | 一种基于监控二阶统计量的传感器精度下降故障检测方法 | |
CN110308713A (zh) | 一种基于k近邻重构的工业过程故障变量识别方法 | |
CN111639304B (zh) | 基于Xgboost回归模型的CSTR故障定位方法 | |
CN110751217A (zh) | 基于主元分析的设备能耗占比预警分析方法 | |
CN111983994B (zh) | 一种基于复杂工业化工过程的v-pca故障诊断方法 | |
Cong et al. | Multirate partial least squares for process monitoring | |
Hu et al. | Direct projection to latent variable space for fault detection | |
Alvarez et al. | Batch process monitoring in the original measurement's space | |
CN113253682B (zh) | 非线性化工过程故障检测方法 | |
CN114995338A (zh) | 一种基于规范变量分析与js散度融合的工业过程微小故障检测方法 | |
CN110794814B (zh) | 一种基于广义主成分的故障确定方法及系统 | |
CN115047853A (zh) | 基于递推规范变量残差和核主元分析的微小故障检测方法 | |
Babamoradi et al. | Comparison of bootstrap and asymptotic confidence limits for control charts in batch MSPC strategies |
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 |