删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于递推规范变量分析的时变过程故障检测

本站小编 Free考研考试/2020-03-23

商亮亮1,2, 刘建昌1, 谭树彬1, 王国柱1
1.东北大学 信息科学与工程学院,辽宁 沈阳 110819;
2.南通大学 电气工程学院,江苏 南通 226019
收稿日期: 2015-08-04
基金项目: 国家自然科学基金资助项目(61374137);流程工业综合自动化国家重点实验室基础科研业务费资助项目(2013ZCX02-03)。
作者简介: 商亮亮(1981-),男,山东淄博人,东北大学博士研究生;
刘建昌(1960-),男,辽宁黑山人,东北大学教授,博士生导师。

摘要: 由于规范变量分析(CVA)不适应过程的时变特性,容易把正常的过程改变识别为故障.因此,针对时变过程提出一种故障检测方法是十分必要的.采用指数权重滑动平均来更新过去观测矢量的协方差矩阵.递推CVA有较高的计算负荷是需要解决的关键问题.通过引入一阶干扰理论来递推更新Hankel矩阵的奇异值分解(SVD).与普通奇异值分解相比,显著降低了递推算法的计算负荷.将提出的基于一阶干扰理论的递推规范变量分析(RCVA-FOP)应用于田纳西伊斯曼化工过程中.仿真结果表明,所提出方法不仅能有效适应过程的时变特性,而且可以有效检测到两种类型的故障.
关键词:一阶干扰理论规范变量分析时变过程故障检测
Recursive Canonical Variate Analysis for Fault Detection of Time-Varying Processes
SHANG Liang-liang1,2, LIU Jian-chang1, TAN Shu-bin1, WANG Guo-zhu1
1.School of Information Science & Engineering, Northeastern University, Shenyang 110819, China;
2.School of Electrical Engineering, Nantong University, Nantong 226019, China
Corresponding author: SHANG Liang-liang, E-mail: sll-15@163.com
Abstract: Because CVA (canonical variate analysis) is unable to adapt the characteristics of time-varying processes, by which the normal changes of the process is easily identified as faults, it is very necessary to propose a monitoring approach for time-varying processes. The exponential weighted moving average approach was adopted to update the covariance of the past observation vectors. The most critical problem faced by recursive CVA algorithm is the high computation cost. To reduce the computation cost, the first order perturbation theory was introduced to update recursively the singular value decomposition (SVD) of the Hankel matrix. The computation cost of recursive SVD based on the first order perturbation theory is significantly less compared to the SVD. Recursive canonical variate analysis based on the first order perturbation (RCVA-FOP) was applied in the Tennessee Eastman chemical process. Simulation results indicate that the proposed method not only can effectively adapt to the normal change of time-varying processes, but also can detect two types of faults.
Key Words: first order perturbation theoryCVA (canonical variate analysis)time-varying processesfault detection
为了提高现代工业过程的生产效率,近年来故障检测越来越重要.故障检测有三种主要方法:基于知识的方法、基于模型的方法、基于数据驱动的方法[1].基于知识的方法的缺点是需要大量不容易获取的过程知识.尽管数据驱动方法已被广泛应用,但在过去的二十年中基于子空间辨识模型的方法吸引了过程建模和监控领域的学者和专家的关注.常规的子空间辨识方法包括规范变量分析(CVA),子空间状态空间系统辨识的数值算法(N4SID),多变量输出误差状态空间(MOESP)[2]. CVA对线性过程静态操作条件表现出了非常好的性能,但是,许多工业过程需要进行预定的正常改变,比如设定值改变、化学反应器进料比的改变等.因此,针对时变过程开发在线故障检测方法是非常必要的.
递推辨识方法存在的关键问题是较高的计算负荷.根据获取状态变量和可观测矩阵的不同,递推子空间辨识方法可以分为两类:一类是避免使用奇异值分解.Gustafuson[3]首次提出了基于投影近似子空间跟踪(PAST)的递推子空间辨识方法.Oku等[4]提出了使用梯度类型子空间跟踪的递推4SID算法.Houtzager等[5]提出了基于predictor的递推辨识方法.另一类是直接使用奇异值分解或降低奇异值的计算负荷.Choi等[6]提出了基于CVA更新状态空间模型的自适应故障检测方法.Ding等[7]提出了基于等价空间和观测器的故障检测和隔离系统.
为了降低Hankel矩阵奇异值分解的计算负荷,引入一阶干扰理论实现Hankel矩阵的递推奇异值分解更新.与常规奇异值分解相比,基于一阶干扰理论的递推奇异值分解显著降低了计算负荷.本文提出了基于一阶干扰理论的递推规范变量分析(RCVA-FOP)时变过程故障检测方法,并将该方法用于田纳西伊斯曼化工过程的在线故障检测.
1 规范变量分析多变量状态空间模型为
其中:x(k)Rny(k)Rmy是采样时刻k的状态和观测矢量;ARn×nBRn×nCRmy×nDRn×n是系统矩阵;w(k)v(k)是随机干扰,假定是零均值正态分布白噪声.
在CVA方法中,测量矢量由过去和未来的测量值扩展而成,分别形成过去、未来矢量zp (r)yf (r),其中r表示一类序数.
其中,p代表过去和未来观测值的长度.一般情况,n/myp.
设置r=p+1, p+2, …, p+N,过去和未来Hankel矩阵YpRmp×NYfRmp×N定义如下:
一组由l个观测值组成的变量,yp (p+1)最后的两个元素是y(1)u(1),然而yf (p+N)最后的元素是y(l).因此,这些Hankel矩阵的最大列数是N=l-p-f+1.
过去和未来观测值的协方差和互方差矩阵可以使用如下公式计算:
CVA试图找到aTyf (r)bTyp (r)之间的最优线性组合,以便于过去和未来观测值之间的相关关系最大化.该最优化问题的解决方案可由Hankel矩阵H的奇异值分解求得.
对上式进行SVD,规范变量z(r)可基于过去测量值求得:
其中J=VTΣpp-1/2是转换过去测量值到规范变量空间的转换矩阵.规范变量空间可以被分成两个部分:主空间和残差空间,它们的维数由奇异值的维数确定.前n个主要奇异值确定主空间x(r)Rn,并且剩余的奇异值确定残差空间r(r)Rmp-n,其中n < mp,奇异值个数n应事先给出.因此,规范变量可以写成:
其中,Jx=VxTΣpp-1/2Vx包含V的前n列.状态变量不仅是估计的规范变量的子集,也是过去观测矢量的线性组合.因此,像规范变量一样状态变量定义为过去观测矢量的线性组合.在系统状态被确定之后,矩阵A, B, C, D可以通过线性最小二乘递推来估计.最后,可以重新计算系统的输出.
2 基于递推规范变量分析的故障检测对多模过程,每个模式有不同的特性,例如变量之间均值、方差和相关关系结构的改变导致主要状态变量个数的改变.每次采集到新的测量变量时,均值和协方差随着模型结构改变程度进行更新.本文中,采用指数权重滑动平均而不需要回调过去的训练数据.
过去观测矢量的协方差矩阵可以通过指数权重滑动平均来估计.
其中β是遗忘因子.
2.1 基于一阶干扰理论更新状态矢量基于干扰理论更新奇异值和奇异值矢量的根源可以参考文献[8]中提到的结果.在例子中,奇异值和奇异值矢量的初始集来自离线辨识.在线更新状态矢量开始于当前测量值.zk表示为zk=zf (r)zf (r)Tzf (r),其中,zf (r)是列矢量包含当前和未来测量值.相似地,增补变量矢量可以利用以前的测量值zτ=zp (r)zp (r)Tzp (r).
定义如下过渡变量pu, 0=qu, 0=[0 … 0].
首先,更新左奇异值矢量执行如下操作:
下标1≤iq代表矩阵列或矢量的第i元素.
然后,更新右奇异值矢量,定义两个过渡变量:pv, 0=qv, 0=[0 … 0].
右奇异值矢量可计算如下:
奇异值矢量更新如下:
最后标准化后如下:
给定更新的奇异值矢量,新的规范变量可以获得.应该事先给出递推模型的阶数,计算状态和残差矢量的方程如下:
其中:x(r)是在时间点r的状态矢量; Vx包含Vr在时间点r的前n列; zp (r)是递推建模过程中新的过去观测矢量.
仅对ΦzRq×q的完整更新进行复杂度和计算成本计算.如果表示想要估计的子空间的秩,通常l < < q.许多可供选择的要求较少操作的方案被开发来避免极度高的计算复杂度.几种常见算法SVD、逆迭代、秩一更新、投影近似子空间跟踪(PAST)和RSVD-FOP的复杂度分别为O(q3),O(q2l),O(q3),O(q2),O(q2).与SVD相比,基于一阶干扰理论的递推SVD的计算成本大大减小.
2.2 确定主空间的维数确定CVA模型阶数的方法有很多.累积方差百分比(CPV)、Akaike信息评价准则(AIC)、最大似然比检验和交叉验证等.然而,不是所有方法都适合确定递推CVA的阶数,比如交叉验证法.本文采用了累积方差百分比来确定主要的奇异值个数:
2.3 故障检测统计量及其控制限为了监控过程故障,本文采用两个分别来自状态空间和噪声空间的统计量.基于状态空间定义T2统计量为T2=x(r)Tx(r).
其置信水平α的控制限计算如下:
其中,Fn, N-n(α)在置信水平为α时自由度为nN-nF分布的临界值.
总的过程噪声统计量计算公式如下:
基于正态分布假设,Q统计量的控制限计算如下:
其中:θi=hr-nhr是Hankel矩阵的秩;h0=1-(2θ1θ3/3θ22);cα是对应1-α百分位的正态偏差.
3 应用于田纳西伊斯曼化工过程采用田纳西伊斯曼(TE)化工过程仿真平台产生的数据进行仿真研究和分析.测试数据集使用TE过程的simulink仿真代码产生.在第300采样时刻将产品设定值从22.89改变为24.89,并持续到第800采样时刻.为了验证提出的方法,假定设定值改变属于正常操作.两种类型的过程变量故障F4,F10在第800采样时刻发生并持续到仿真时间停止.首先使用测试数据中前200个采样时刻的数据来建立CVA模型,并取得k时刻的过去观测矢量的协方差矩阵.然后,再使用测试数据进行模型更新.因为计算时间和p2成正比,在本次仿真研究中将p的值设置为3.
首先,对故障4进行仿真对比研究.图 1图 2分别给出了基于CVA和基于RCVA-FOP方法的在线故障检测结果.
图 1(Fig. 1)
图 1 TE过程故障4的CVA监控图Fig.1 Monitoring charts of CVA for fault 4 of TE process

图 2(Fig. 2)
图 2 TE过程故障4的RCVA-FOP监控图Fig.2 Monitoring charts of RCVA-FOP for fault 4 of TE process

因为基于常规CVA的检测方法没有及时更新模型,所以产生了误报.如图 1所示,其监控统计量从第300采样时刻以后开始逐渐超过控制限,产生了较高的误报率.与基于CVA的检测方法对比,如图 2所示基于RCVA-FOP的在线检测方法不仅能够通过模型更新适应过程的正常改变,而且还能检测到故障4 (反应器冷却水入口温度阶跃变化),降低了误报率,提高了过程的在线监控性能.
然后,对故障10进行仿真对比研究.图 3图 4分别给出了基于CVA和基于RCVA-FOP方法的在线故障检测结果.
图 3(Fig. 3)
图 3 TE过程故障10的CVA监控图Fig.3 Monitoring charts of CVA for fault 10 of TE process

图 4(Fig. 4)
图 4 TE过程故障10的RCVA-FOP监控图Fig.4 Monitoring charts of RCVA-FOP for fault 10 of

因为基于CVA的检测方法没有根据过程的正常变化及时更新模型,如图 3所示其监控统计量从第300采样时刻后开始逐渐超过控制限,增大了误报率.与基于CVA的检测方法对比,如图 4所示提出的RCVA-FOP能够通过模型更新适应过程的正常改变,同时还能检测到故障10 (流2中C的进料温度随机变化),大大降低了误报率,提高了过程的在线监控性能.由于故障产生初期变化比较缓慢,所以存在一定的检测延时.
4 结语将提出的递推规范变量分析应用于田纳西伊斯曼化工过程中进行仿真研究.基于一阶干扰理论的递推奇异值分解降低了算法的计算负荷.仿真结果表明,提出的方法不仅可以有效适应过程的正常改变,而且可以检测到两种类型的故障.
参考文献
[1]Sammaknejad N, Huang B, Fatehi A, et al. Adaptive monitoring of the process operation based on symbolic episode representation and hidden Markov models with application toward an oil sand primary separation[J].Computers & Chemical Engineering, 2014, 71 : 281–297.
[2]Qin S J. An overview of subspace identification[J].Computer & Chemical Engineering, 2006, 30(10/11/12) : 1502–1513.
[3] Gustafsson T.Recursive system identification using instrumental variable subspace tracking[C]// Proceedings of the 11th IFAC Symposium on System Identification.Fukuoka, 1997.
[4]Oku H, Kimura H. Recursive 4SID algorithms using gradient type subspace tracking[J].Automatica, 2002, 38(6) : 1035–1043.DOI:10.1016/S0005-1098(01)00286-2
[5]Houtzager I, van Wingerden J W, Verhaegen M. Recursive predictor-based subspace identification with application to the real-time closed-loop tracking of flutter[J].IEEE Transactions on Control Systems Technology, 2012, 20(4) : 934–949.DOI:10.1109/TCST.2011.2157694
[6]Choi S W, Martin E B, Morris A J, et al. Adaptive multivariate statistical process control for monitoring time-varying processes[J].Industrial & Engineering Chemistry Research, 2006, 45(9) : 3108–3118.
[7]Ding S X, Zhang P, Naik A, et al. Subspace method aided data-driven design of fault detection and isolation systems[J].Journal of Process Control, 2009, 19(9) : 1496–1510.DOI:10.1016/j.jprocont.2009.07.005
[8]Naik A S, Yin S, Ding S X, et al. Recursive identification algorithms to design fault detection systems[J].Journal of Process Control, 2010, 20(8) : 957–965.DOI:10.1016/j.jprocont.2010.06.018

相关话题/过程 变量

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 二元碱度对印尼钒钛矿烧结过程及烧结矿质量的影响
    高强健,魏国,刘常升,沈峰满东北大学冶金学院,辽宁沈阳110819收稿日期:2015-08-02基金项目:国家自然科学基金资助项目(51604069);中国博士后科学基金资助项目(2016M591445);东北大学博士后基金资助项目(20160302);国家自然科学基金-辽宁联合基金资助项目(U15 ...
    本站小编 Free考研考试 2020-03-23
  • 在现场确认过程中,若学生证丢失,可否用在读证明代替
    提问问题:在现场确认过程中,若学生证丢失,可否用在读证明代替学院:提问人:18***47时间:2018-09-1910:22提问内容:在现场确认过程中,若学生证丢失,可否用在读证明代替回复内容:请尽快到教务部门补办。 ...
    本站小编 东华大学 2019-11-25
  • 2019年化工过程机械复试录取比例
    提问问题:2019年化工过程机械复试录取比例学院:能源与动力工程学院提问人:17***28时间:2019-09-2312:44提问内容:你好老师,我想问一下2019年进入复试人数以及最终录取人数,录取考生成绩大概分布?谢谢回复内容:复试和录取人数比例不低于1.2:1 ...
    本站小编 南京理工大学 2019-11-25
  • 你校的化工学院是否招收080706化工过程机械专业
    提问问题:你校的化工学院是否招收080706化工过程机械专业学院:化工学院提问人:15***83时间:2017-09-2111:06提问内容:1、你学院是否招收“080706化工过程机械”专业(学术硕士,全日制)?2、如果招收,计划招收人数是多少人(不含推免生),报录比是多少?3、考试有哪些科目?( ...
    本站小编 南京理工大学 2019-11-25
  • 本科是针推的如果想考贵校的中内会不会在复试的过程中容易被刷?
    提问问题:中内学院:不区分院系所提问人:17***47时间:2019-09-1912:11提问内容:想请问一下老师1.我本科是针推的如果想考贵校的中内会不会在复试的过程中容易被刷?有没有报录比?2.确定导师是在复试之后分配还是复试之前就需要先联系导师?如果联系导师是在研究生官网上找么?3.中医内科学 ...
    本站小编 南京中医药大学 2019-11-25
  • 化工过程机械
    提问问题:化工过程机械学院:机械与动力工程学院提问人:18***49时间:2017-09-2011:51提问内容:老师您好,请问机械院的化机今年接收的推免名额是多少?本校推免到这个专业的多吗?会不会排斥二本院校的推免生啊?回复内容:参照往年推免人数情况。请去研究生院官网查看近几年录取情况,里面有往年 ...
    本站小编 南京工业大学 2019-11-25
  • 化工过程机械
    提问问题:化工过程机械学院:机械与动力工程学院提问人:18***80时间:2015-09-2309:31提问内容:请问一下哪里能看到往年录取线或录取名单还有就是麻烦告知今年化机打算统招多少人看到以前的是具体咨询学院派那怎么咨询麻烦一一告知谢谢回复内容:http://210.28.203.85/vie ...
    本站小编 南京工业大学 2019-11-25
  • 化工过程机械招生人数?
    提问问题:化工过程机械招生人数?学院:机械与动力工程学院提问人:21***cn时间:2014-09-2414:05提问内容:老师好,请问化工过程机械招生人数多少?谢谢了回复内容:具体咨询学院,谢谢! ...
    本站小编 南京工业大学 2019-11-25
  • 化工过程机械
    提问问题:化工过程机械学院:机械与动力工程学院提问人:21***cn时间:2014-09-2413:46提问内容:老师您好,请问化工过程机械历年复试分数线大概是多少?录取比呢?材料力学(机)是什么意思,指定课本使用哪个版本吗?提供真题吗?复试科目有哪些?谢谢了回复内容:具体咨询学院,谢谢! ...
    本站小编 南京工业大学 2019-11-25
  • 化工过程机械
    提问问题:化工过程机械学院:化工学院提问人:21***cn时间:2014-09-2613:16提问内容:老师好,请问化工过程机械专业对英语有没有特别的要求?四六级?谢谢老师了!回复内容:中国矿业大学2015年硕士研究生招生专业目录及招生简章http://yjsb.cumt.edu.cn/Infoma ...
    本站小编 中国矿业大学 2019-11-25