CN110705129B - 一种基于全相关动态kpls的故障诊断方法 - Google Patents
一种基于全相关动态kpls的故障诊断方法 Download PDFInfo
- Publication number
- CN110705129B CN110705129B CN201911048863.5A CN201911048863A CN110705129B CN 110705129 B CN110705129 B CN 110705129B CN 201911048863 A CN201911048863 A CN 201911048863A CN 110705129 B CN110705129 B CN 110705129B
- Authority
- CN
- China
- Prior art keywords
- matrix
- variables
- input
- kpls
- output
- 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 54
- 238000003745 diagnosis Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 230000002596 correlated effect Effects 0.000 claims abstract description 8
- 239000013598 vector Substances 0.000 claims description 28
- 230000003416 augmentation Effects 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 8
- 238000001311 chemical methods and process Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 description 20
- 238000004519 manufacturing process Methods 0.000 description 18
- 238000004422 calculation algorithm Methods 0.000 description 15
- 238000001514 detection method Methods 0.000 description 11
- 238000005086 pumping Methods 0.000 description 10
- 238000012360 testing method Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 230000000875 corresponding effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012880 independent component analysis Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明提出一种基于全相关动态KPLS(all correlated dynamic KPLS,ADKPLS)的非线性化工过程故障诊断方法。该方法首先把原始数据分为输入变量与输出变量两个部分,分别进行动态特征分析研究,使组成的数据矩阵能很好地反映变量间的动态关系。然后分析证明了KPLS中输出变量的变化会影响到输入残差空间,设计一个输出变量辅助矩阵,表征输入变量与输出变量的全相关性。最后,采用基于输入变量与输出变量之间的全相关信息构建贡献图以识别故障源变量。与KPLS相比,本发明在输入和输出变量之间建立了更直接的关系。
Description
技术领域
本发明涉及故障诊断领域,具体涉及一种基于全相关动态KPLS的故障诊断方法。
背景技术
由于现代工业过程在生产流程和自动化程度上都非常复杂,对可靠性和安全性的要求越来越高,系统任何部分的故障都可能导致整个系统故障。因此,早期检测潜在或发生的故障在工业过程中非常重要。在过去的几十年中,多变量统计过程监测(MSPM)技术被广泛用于复杂工业过程的故障诊断。MSPM主要包括主成分分析(PCA),独立成分分析(ICA)和偏最小二乘法(PLS)。通过使用统计分析方法从大量输入和输出过程数据中提取潜在特征,以诊断工业过程的故障。因此,称为数据驱动方法的MSPM通常用于工业系统中的过程监控和故障诊断。
在MSPM技术中,PLS是复杂工业过程中最常用的数据驱动方法之一。在PLS算法中,从过程变量中提取潜在变量。PLS的建模目标是使用潜在变量来描述输入和输出变量之间的关系,使用此模型可以通过输入变量预测输出变量,同时PLS方法使用不相关的潜变量进行建模和分析。因此,PLS具有处理缺陷数据和高度相关的过程数据的能力。针对较为复杂的非线性过程,标准PLS不能更好实现的监控和故障监测。为了解决这个问题,已经提出了许多非线性PLS方法来克服非线性输入和输出数据。在这些方法中核偏最小二乘法(KPLS)是最常用的方法,它是通过引入核函数来解决非线性问题,非线性输入数据被映射到高维特征空间。在高维空间中,将线性不可分割性问题转化为线性分离问题。因此,KPLS方法可以有效地解决非线性过程数据的线性化问题。
基于KPLS模型,已经提出了许多故障诊断方法并成功应用于工业系统。有些现有技术,提出了动态核偏最小二乘建模方法来建立测量和质量指标之间的动态关系。有些现有技术,开发了一个全核PLS(T-KPLS)用于非线性质量相关的过程监控,其中输入空间被分为四个部分。有些现有技术,提出了一种定向核偏最小二乘监测方法来提取输出相关变量,构建了输入和输出变量之间更精细的关系。有些现有技术,开发了一种基于最优偏好矩阵(OPM)的改进KPLS方法,以解决主成分误导性解释的问题,从而调整过程变量的分布和协方差矩阵的特征值。有些现有技术,提出了一种基于多核局部最小二乘法的分散故障检测方法,用于检测大规模工业过程。
尽管标准KPLS及其最近的修改方法被广泛用于监视和诊断工业故障,但在工业过程监控中仍然没有完全解决一些问题。在残余子空间中,残差局部变量的显著变化过程仍然与输出变量有关。所以,它使得KPLS方法不适合用潜变量来诊断输出相关的故障。因此,有效的故障诊断方法对于现代工业过程仍然是必不可少的。
发明内容
本发明的发明目的是提供一种适合用潜变量来诊断输出相关故障的故障诊断方法。
本发明提供一种基于全相关动态KPLS的故障诊断方法,包括以下步骤:
离线建立主元模型,获取统计量控制限;
获取实时运算数据并计算实时统计量;
将统计量控制限与实时统计量进行比对,若实时统计量超出统计量控制限则发出报警。
离线建立主元模型,获取统计量控制限包括以下步骤:
获取正常运行的数据;
对正常运行的数据进行预处理;
建立ADKPLS模型,
所述建立建立ADKPLS模型步骤包括以下步骤:
第一步:将原始数据集X分为输入变量A和输出变量B两个矩阵,分别确定动态阶次,分别构建增广矩阵Φ(x)和Y,其集合为增广矩阵Xk(x),
Xk(x)表示原始训练集X的增广矩阵;
此处为计算全相关辅助矩阵,替代原始数据中Φ(x)和Y的相关性,为辅助矩阵与Y进行主元提取得到主要的相关性变量提供依据。
k表示主元个数,uk表示全相关辅助矩阵;
第四步:计算输入向量得分矩阵tk:tk=Kuk,
tk表示输出变量的得分矩阵,K表示Φ(x)Φ(x)T的核运算;
第五步:标准化得分矩阵:tk=tk/||tk||;
pk表示输入向量负载矩阵;
第八步:重复第三部至第七步直到收敛;
Kk+1表示更新后的K值;
第十步:返回第四步。
本发明的有益效果是:本发明在输入和输出变量之间建立了更直接的关系。在抽油机生产过程上的仿真结果表明,本发明具有良好的故障检测与识别性能。
附图说明
图1KPLS残差子空间与输出变量之间的相关性示意图。
图2基于动态主元分析的故障检测与诊断流程示意图。
图3KPLS、ADKPLS检测故障1结果示意图。
图4故障1所有变量贡献率示意图。
图5KPLS、ADKPLS检测故障2结果示意图。
图6故障2所有变量贡献率示意图。
具体实施方式
传统KPLS(KPLS表示基于核函数的偏最小二乘算法)在构造输入变量与输出变量时未考虑样本间的动态特性,使故障样本信息易被其他样本掩盖;同时,KPLS未考虑所有变量对故障样本的影响,使变量间隐藏信息不能被完全表达。针对以上问题,本发明提出一种基于全相关动态KPLS(all correlated dynamic KPLS,ADKPLS)的非线性化工过程故障诊断方法。该方法首先对原始数据进行动态特征分析研究,使组成的数据矩阵能很好地反映变量间的动态关系。然后分析证明了KPLS中输出变量的变化会影响到输入残差空间,设计一个输出变量辅助矩阵,表征输入变量与输出变量的全相关性。最后,采用基于输入变量与输出变量之间的全相关信息构建贡献图以识别故障源变量。与KPLS相比,在输入和输出变量之间建立了更直接的关系。在抽油机生产过程上的仿真结果表明,所提方法具有良好的故障检测与识别性能。
下面对本发明所提出方法的推导过程和具体实施方式进行说明。
1动态回归模型
1.1 DKPLS模型
在抽油机生产过程中,获得的数据中存在着非线性和时序相关性等问题,采用核函数转换为线性问题。针对原始数据中存在的时序相关特性,采用多元回归模型AR(B)使数据间具有动态特征,可表示为:
x(k)=α1x(k-1)+α2x(k-2)+αBx(k-B)+…+αB+1y(k-1)+αB+2y(k-2)+αky(k-B)+εk(1)
其中,x(k)是动态变量,α1,α2,αD是动态回归系数,εk是动态误差,B是该模型的时滞值。用x(k)=(k-1),(k-2),…,(k-B)表示原始数据矩阵,能够良好地反映变量间动态特征。
1.2自回归模型阶次确定
DAKPLS实现自相关问题的关键就是要明确自回归阶次B,通常取B=1或2。如果需要更为精确的时滞值阶次,可采用统计假设检验方法。其基本思想是判断各变量之间是否具有自相关性,从而得到更为精确的阶次B。
由多元回归模型推导时滞值B阶的模型为xt=α1xt-1+α2xt-2+…+αpxt-B+εt,经计算得到该模型的残差平方和为SB;同理,B-1阶的模型为xt=α1xt-1+α2xt-2+…+αB-1xt-B+1+εt,残差平方和为SB-1。
假设B0:bp=0成立时,可作F分布统计量为F=SB-SB-1/2SB/(n-2B),选取显著水平α,以分子自由度1,分母自由度n-B,查表得Fα。若F>Fα,则说明B0不成立,B阶与B-1阶模型之间存在差异,取B阶。反之,取B-1阶。AR模型阶次B确定后,就可以得到原始数据对应的增广矩阵XB如下:
通过分别构建输入变量和输出变量的增广矩阵,使其分别体现数据之间的良好动态特性,挖掘出各自空间中存在的潜在变量,为后续构建全相关辅助矩阵进行相关性分析奠定基础。
2构建ADKPLS模型
在本节中,提出ADKPLS算法。根据KPLS算法推导输入残差子空间与输出变量的相关性,构建输出变量辅助矩阵,基于动态模型,提出ADKPLS算法。
2.1 KPLS算法
在KPLS中,输入与输出向量分别建立回归模型如下:
其中,i表示KPLS中保留的隐变量个数,通过交叉检验得到,tl为输入向量的得分矩阵,pl为输入向量的负载矩阵,ql为输出向量负载矩阵,φi(x)和分别为输入和输出向量残差矩阵。KPLS主要参数如下表所示:
2.2输入残差子空间与输出变量的相关性
KPLS中输入残差子空间φi(x)与输出向量Y存在着相关性,证明如下:
tl为KPLS中的主元得分向量,φi(x)为KPLS中的残差向量。
当i≥l时,φi T(x)tl=0 (5)
公式(4)可以化简为
由于在KPLS中输入、输出向量的负载向量pi和qi分别可以写为
pi=φi-1 T(x)ti/ti Tti (7)
将公式(7)和(8)带入公式(6)得到
图1所示为KPLS残差子空间与输出变量关系图,其中RA为KPLS中输入向量φ(x)变化范围,RB为KPLS中输出向量Y的变化范围,且φi(x)与为正交状态,若残差子空间与输出向量之间无关,即通过图2可以看出,与ti为正交状态,即第i次迭代时的F与t为正交状态,同样,第i+1次迭代时,由于迭代次数的不同,本发明并不能知道第i+1次迭代时的主元空间与第i次迭代时的残差空间之间的关系,从而无法保证与ti+1是否满足正交规则,实际上,因此,输入残差与输出变量之间存在着相关性。
2.3构造全相关辅助矩阵
因为Φ(x)与ξy完全无关,有
cov(Φ(x),ξy)=ε(Φ(x),ξy)=0 (12)
所以公式(11)可以简化为
为了使输入变量和输出变量之间的协方差达到最大,采用拉格朗日乘子法
其中w1和c1分别是输入变量E0和输出变量F0的轴向量(E0=X和F0=Y1),对s分别求关于w1,c1,τ,υ的偏导且置0,记θ=2τ=2υ,则θ是优化问题的目标函数且使θ达到最大必须有
在KPLS推导中,公式(16)可以写为
由上得出
φi-1(x)φi-1 T(x)F0ci=θφi-1(x)wi (18)
由公式(16)得到
φi-1(x)wi=φi-1(x)φi-1 T(x)F0ci/θ (20)
将公式(20)引入(19)结合F0ci=ti有
YYTφi-1(x)pi=ti (22)
φi(x)=Φ(x)-YYTΦ(x)/YTY (23)
其中,1/YTY是用于归一化输出相关变化的值,所以
因为ΦT(x)Φ(x)是一个满秩方阵,其中I为单位矩阵。所以,通过构造全相关辅助矩阵能够达到输出变量空间与输入残差空间之间不再存在相关变量,从而得到输入向量与输出向量之间的所有相关变量,再进行主元特征提取,能够更好的体现输入向量与输出变量之间的主要相关变量。
2.4 ADKPLS算法
基于动态回归模型理论上,推导了输入变量残差与输出变量存在着相关性,从而构造全相关辅助矩阵基于ADKPLS算法的故障检测与诊断步骤包括数据动态分析、变量相关性分析、构建全相关模型和在线故障检测与诊断,具体步骤如下:
ADKPLS具体算法如下表所示:
在模型实现主元提取之前,需对高维空间进行中心化运算:
3基于ADKPLS算法的故障检测
与PLS和KPLS类似,Hotelling’s T2统计量和SPE统计量仍是IDKPLS故障检测两个重要统计量。T2统计量反映变量在主元空间的变化:
T2=tnewΛ-1tnew (27)
其中,Λ-1是得分矩阵的协方差的逆。其控制线可通过F分布获得T2=i(n2-1)×F(i,n-i,α)/n(n-i),F(i,n-i,α)是对应于检验水平为α,自由度为i和n-i条件下的F分布临界点。
SPE统计量反映样本向量在残差空间投影的变化:
其中,n1是测试数据样本个数,xt=knewuk/(tk Tknewuk),knew是测试核矩阵。其控制限可通过其近似分布来计算:
4实验研究
在本节中,通过中国大港油田油井的实际生产数据用于验证本发明提出的算法ADKPLS的有效性。在实际工业生产中,抽油机的电压和电流参数反映了故障情况,电压和电流数据集不能直接用来建立模型,可通过傅里叶模型变幻得到。本发明在实验中使用抽油机的相应电压和电流数据来建立故障诊断模型,这里故障1是连抽带喷,故障2是气体影响,故障3是供液不足。本发明从2016年11月20日至2018年4月16日抽取了连续1000组数据集,训练集由500组样本数据构成,测试集由500组样本数据构成。在训练集中,前250个数据集是正常样本,后250个数据集是故障样本。在测试集中,前300个数据集是正常样本,后200个数据集是故障样本。
对于故障1,监测结果如图3所示。图3a可以看出,KPLS在T2和SPE均在第300个样本处发生跳变,但控制线不能很好的检测到故障,正常样本依旧属于故障,误报率结果高。图3b,ADKPLS在第300个样本成功检测到故障,同时正常样本数值超过控制限较少,效果相对比较明显;SPE统计量中正常样本比T2统计量超过控制限明显较少,且输出相对稳定变化,说明输入残差子空间的变化确实受到输出变量的影响。
贡献率通过分析每个变量对检测指标的敏感度来解释故障原因,可以成功定位故障组件。由图4可知,图4a中T2和SPE贡献率最大分别是变量5和变量37,图4b中T2和SPE贡献最大分别是变量41和变量28,故障1在实际生产中,T2统计量是由变量41引起的故障,即是由于电压、电流引起的故障,从而得到第300个采样点T2统计量贡献图法故障诊断与故障1所描述故障相符。证明通过动态回归模型使使数据间的相关性更为密切。
对于故障2,监测结果如图5所示。图5a可以看出,T2能在第300个样本点处成功检测到故障,同时正常样本数值超过控制限较少,SPE因控制限太小,不能正常进行故障检测。图5b中T2和SPE均能在第300个样本点处成功检测到故障,并且效果显著。SPE统计量中正常样本比T2统计量超过控制限明显较少,且输出相对稳定变化。与KPLS相比,ADKPLS方法在监测故障2方面更有效。
由图6可知,图6a中T2和SPE贡献率最大分别是变量23和变量28,图3b中T2和SPE贡献最大分别是变量36和变量28。故障2在实际生产中,T2统计量是由变量36引起的故障,SPE是由变量28引起的故障,从而得到第300个采样点SPE统计量贡献图法故障诊断与故障2所描述故障相符。证明通过ADKPLS模型对数据间隐藏变量的挖掘,使数据间的相关性更为密切;再构造输入变量与输出变量全相关特征,使数据间故障变量定位更准确。
通过对抽油机生产过程数据进行统计过程监控,利用KPLS和ADKPLS算法得到的各变量贡献率进行分析,最终用T2和SPE统计量的误报率和漏报率作为衡量监控性能的指标,结果如下表所示。
表1统计过程监控性能对比(%)
在多故障分类中,计算故障的误报率和漏报率时,其他故障样本被认为是“正常样本”。因此,误报警率是被识别为正常的“故障样本”的数量与正常样本的总数的比率;漏报警率是故障的“正常样本”的数量与故障样本的总数的比率。由表分析可知,所提出算法的误报警率和漏报警率在控制限基础上都比较低,说明基于ADKPLS的多元统计过程监控方法在抽油机生产过程故障检测与诊断方面具有优良的监控性能,适用于复杂的非线性和动态工业过程。
5结论
在本发明中,结合动态自回归模型,使原始数据间具有强动态性,再分析证明了KPLS输入残差变量与输出变量之间的相关性,通过构造全相关辅助矩阵,使输入变量与输出变量全相关。结果证明了ADKPLS监测方法比KPLS监测方法表现出更好的监测性能,表明所提出的ADKPLS监测方法对于监测质量相关过程更有效。
本发明提出了一种新的数据驱动的ADKPLS算法用于抽油机系统。首先,全相关性的KPLS与动态模型进行关联,构建后的数据能够良好地体现变量间隐藏的动态特性。其次,输出变量的变化会影响到输入残差子空间,通过构建一个输出变量全相关辅助矩阵,使其与所有输入变量具有相关性。最后,在抽油机系统中输入和输出变量之间有着更直接的关系,实验结果表明了所提出的ADKPLS算法的有效性。
Claims (1)
1.一种基于全相关动态KPLS的故障诊断方法,其特征在于,包括以下步骤:
步骤一、离线建立主元模型,获取统计量控制限,包括:第一步,获取正常运行的数据;第二步,对正常运行的数据进行预处理;第三步,建立ADKPLS模型,所述建立ADKPLS模型包括以下步骤:
第一步:将原始数据集X分为输入变量A和输出变量B两个矩阵,分别确定动态阶次,分别构建增广矩阵Φ(x)和Y,其集合为增广矩阵Xk(x),
Xk(x)表示原始训练集X的增广矩阵;
第四步:计算输入向量得分矩阵tk:tk=Kuk;
tk表示输入向量的得分矩阵,K表示Φ(x)Φ(x)T的核运算;
第五步:标准化得分矩阵:tk=tk/||tk||;
第六步:计算输入向量负载矩阵:Pk=YTtk
Pk表示输入向量负载矩阵;
第七步:计算迭代后的uk:uk=Ypk
第八步:重复第三步至第七步直到收敛;
第九步:计算残差矩阵:K=(I-tktk T)K(I-tktk T)
Y=Y-tktk TY
左边的K表示更新后的K值;
第十步:返回第四步;
步骤二、获取实时运算数据并计算实时统计量;
步骤三、将统计量控制限与实时统计量进行比对,若实时统计量超出统计量控制限则发出报警。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911048863.5A CN110705129B (zh) | 2019-10-31 | 2019-10-31 | 一种基于全相关动态kpls的故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911048863.5A CN110705129B (zh) | 2019-10-31 | 2019-10-31 | 一种基于全相关动态kpls的故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110705129A CN110705129A (zh) | 2020-01-17 |
CN110705129B true CN110705129B (zh) | 2023-03-14 |
Family
ID=69202969
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911048863.5A Active CN110705129B (zh) | 2019-10-31 | 2019-10-31 | 一种基于全相关动态kpls的故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110705129B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112229212A (zh) * | 2020-08-18 | 2021-01-15 | 广东工业大学 | 一种基于动态主元分析的辊道窑能耗异常检测方法 |
CN113282884B (zh) * | 2021-04-28 | 2023-09-26 | 沈阳航空航天大学 | 通用根因分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103472820A (zh) * | 2013-09-18 | 2013-12-25 | 哈尔滨工业大学 | 一种基于偏最小二乘算法的推进系统故障诊断方法 |
CN104331058A (zh) * | 2014-11-24 | 2015-02-04 | 重庆科技学院 | 一种故障诊断本体构建方法 |
CN104914847A (zh) * | 2015-04-09 | 2015-09-16 | 东北大学 | 基于方向核偏最小二乘的工业过程故障诊断方法 |
CN108803520A (zh) * | 2018-06-11 | 2018-11-13 | 宁波大学 | 一种基于变量非线性自相关性剔除的动态过程监测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050176057A1 (en) * | 2003-09-26 | 2005-08-11 | Troy Bremer | Diagnostic markers of mood disorders and methods of use thereof |
US10345046B2 (en) * | 2017-05-25 | 2019-07-09 | Northeastern University | Fault diagnosis device based on common information and special information of running video information for electric-arc furnace and method thereof |
-
2019
- 2019-10-31 CN CN201911048863.5A patent/CN110705129B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103472820A (zh) * | 2013-09-18 | 2013-12-25 | 哈尔滨工业大学 | 一种基于偏最小二乘算法的推进系统故障诊断方法 |
CN104331058A (zh) * | 2014-11-24 | 2015-02-04 | 重庆科技学院 | 一种故障诊断本体构建方法 |
CN104914847A (zh) * | 2015-04-09 | 2015-09-16 | 东北大学 | 基于方向核偏最小二乘的工业过程故障诊断方法 |
CN108803520A (zh) * | 2018-06-11 | 2018-11-13 | 宁波大学 | 一种基于变量非线性自相关性剔除的动态过程监测方法 |
Non-Patent Citations (3)
Title |
---|
A kernel least squares based approach for nonlinear quality-related fault detection;Wang G 等;《IEEE transactions on industrial electronics》;20170430;第64卷(第4期);全文 * |
偏最小二乘线性模型及其非线性动态扩展模型综述;孔祥玉等;《控制与决策》;20180113(第09期);全文 * |
基于核典型相关分析的故障检测方法;任嵬 等;《控制工程》;20190420;第26卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110705129A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yu | Local and global principal component analysis for process monitoring | |
Ge et al. | Improved kernel PCA-based monitoring approach for nonlinear processes | |
Ma et al. | Hierarchical monitoring and root-cause diagnosis framework for key performance indicator-related multiple faults in process industries | |
JP5431178B2 (ja) | 発電所計測器性能監視予測方法 | |
CN106644162B (zh) | 基于邻域保持嵌入回归算法的环网柜线芯温度软测量方法 | |
CN104714537B (zh) | 一种基于联合相对变化分析和自回归模型的故障预测方法 | |
Zhang et al. | Fault detection and diagnosis of chemical process using enhanced KECA | |
CN104914847B (zh) | 基于方向核偏最小二乘的工业过程故障诊断方法 | |
Wang et al. | Fault isolation based on residual evaluation and contribution analysis | |
Wang et al. | Data-driven sensor fault diagnosis systems for linear feedback control loops | |
CN111680725B (zh) | 基于重构贡献的气体传感器阵列多故障隔离算法 | |
CN110705129B (zh) | 一种基于全相关动态kpls的故障诊断方法 | |
Maleki et al. | Simultaneous monitoring of multivariate process mean and variability in the presence of measurement error with linearly increasing variance under additive covariate model (research note) | |
CN110942258B (zh) | 一种性能驱动的工业过程异常监测方法 | |
CN104503436B (zh) | 一种基于随机投影和k近邻的快速故障检测方法 | |
Xu et al. | Multivariate time delay analysis based local KPCA fault prognosis approach for nonlinear processes | |
CN107918381A (zh) | 一种基于组合核函数的类均值核主元故障诊断方法 | |
Hu et al. | Direct projection to latent variable space for fault detection | |
CN112149054B (zh) | 基于时序扩展的正交邻域保持嵌入模型的构建与应用 | |
CN117575027A (zh) | 一种利用知识增强时空因果的传感器故障诊断方法 | |
CN114371677B (zh) | 基于谱半径-区间主成分分析的工业过程状态监测方法 | |
CN114995338A (zh) | 一种基于规范变量分析与js散度融合的工业过程微小故障检测方法 | |
Ji et al. | Orthogonal projection based statistical feature extraction for continuous process monitoring | |
CN114611606A (zh) | 基于核混合空间投影的故障检测方法 | |
Sun et al. | A MPRM-based approach for fault diagnosis against outliers |
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 |