【技术实现步骤摘要】
本专利技术属于工业过程监测领域，具体涉及一种针对非平稳工业过程的故障监测方法和装置。
技术介绍
1、有效的过程监测是保障工业系统高效稳定运行的关键。一个完整的过程监测主要包括故障检测和故障诊断两个环节。通过及时发现故障并确定故障发生的位置，维护人员可以快速对其处理，降低发生重大事故的风险。然而，在实际的生产中，由于原料的变化、生产目标的调整以及未知的外部干扰等，工业过程往往呈现出非平稳特性，即采集到的观测变量的均值和协方差会随时间发生变化。
2、非平稳工业过程为传统的多元统计过程监测(mspm)带来两个方面的挑战，如图1所示。其一，传统的mpsm方法如主成分分析(pca),独立成分分析(ica)等仅对平稳过程有较好的故障检测效果。而对于非平稳过程，他们无法提取到正常训练数据的非平稳变化趋势，导致在线检测阶段，正常的数据常被误检为故障。而真正的故障效应很容易被非平稳趋势所掩盖，从而出现漏检。其二，在变量的非平稳变化以及故障的涂抹效应的双重影响下，很难准确评价各过程变量对故障的贡献，无法定位到故障变量。
3、针对非平稳过程，近些年已经有了一系列方法去解决故障检测时的误检、漏检以及变量冗余的干扰问题。在研究之初，一些方法通过计算非平稳时间序列的差值来消除非平稳特征的影响。然而，作差处理之后可能会对原始数据的动态信息和故障特征造成破坏，从而影响检测性能。后来，针对非平稳过程的故障检测方法主要有基于自适应策略，基于慢特征分析(sfa)，基于基于协整分析(ca)以及基于平稳子空间分析(ssa)等发展方向。基于自适应策
4、虽然上述方法一定程度解决了非平稳过程故障检测的系列问题，然而实现精准的故障诊断依旧是很大的挑战。故障检测和故障诊断是过程监测中密切相关且互为补充的两个重要阶段，但是关于非平稳过程的故障诊断的研究还并未广泛展开，目前只有一些初步的成果。文献[5]通过将最小绝对收缩和选择算子应用于协整分析，利用稀疏重构的策略实现在线识别故障。然而该方法依据每个变量被选择的频率确定的非平稳变量范围较为宽泛，受缺失重建比的阈值影响较大，可能导致诊断结果的不正确。文献[6]提出了基于主导趋势的逻辑回归的故障诊断方法来监测非平稳过程。该方法使用滑动窗口来收集过程数据和提取主要的趋势特征，用于区分故障并确定故障类型。虽然该方法可以实现较为准确的故障诊断效果，但是受到滑动窗口的大小和频率的影响，可能出现故障检测的延迟，而且高频的在线更新也会造成计算复杂度的提高。综上所述，目前非平稳过程依旧缺乏有效的故障诊断策略。
技术实现思路
1、针对非平稳工业过程的故障检测和诊断难题，本专利技术提供一种针对非平稳工业过程的故障监测方法和装置，实现高故障检测率以及精准确定故障位置。
2、为实现上述技术目的，本专利技术采用如下技术方案：
3、一种针对非平稳工业过程的故障监测方法，包括：
4、利用稀疏平稳投影矩阵对工业过程的多维观测变量进行稀疏线性组合构建多维度平稳源，并构建平稳源的监测统计指标；
5、其中，所述稀疏平稳投影矩阵，采用平稳子空间分析方法并在优化问题中引入l2,p范数作为稀疏约束，并利用由多维观测变量构成的训练数据求解得到；
6、根据每条训练数据，均计算平稳源的监测统计指标，并根据核密度估计计算平稳源的控制限；
7、在线监测阶段，根据实时观测数据计算平稳源的监测统计指标，并与控制限比较以完成故障检测，即判断工业过程是否发生异常；
8、再采用关键信息动态重构的方法，通过将每个观测变量对于监测统计指标的贡献值进行动态重构，通过比较重构贡献值的大小来确定故障变量，完成故障诊断。
9、进一步地，利用稀疏平稳投影矩阵对工业过程的多维观测变量进行稀疏线性组合构建平稳源的表达式为：
10、
11、其中，为t时刻的m个平稳源，为t时刻的d-m个非平稳源，x(t)为t时刻的d维观测变量，x(t)∈rd；bs和bn分别为平稳投影矩阵和非平稳投影矩阵，bs∈rm×d，bn∈r(d-m)×d。
12、进一步地，所述采用平稳子空间分析方法并在优化问题中引入l2,p范数作为稀疏约束，得到优化问题为：
13、
14、其中，bs为平稳投影矩阵；tr()表示括号内矩阵的迹；s为中间量，ui和∑i分别为第i个时期的观测训练数据的均值和协方差，和分别为所有n个时期的观测训练数据的平均均值和平均协方差，λ＞0是正则化系数，用以调节平稳性和稀疏性之间的平衡关系，p表示矩阵范数参数，右角标t表示转置运算，||·||2,p表示l2,p范数；im为m行m列的单位矩阵，m为平稳源的数目。
15、进一步地，求解优化问题以得到平稳投影矩阵具体为：
16、a1、将式(2)所示的优化问题转换为迹的形式：
17、
18、其中，d为中间量：
19、
20、其中，bj代表平稳投影矩阵bs的第j列向量；
21、a2、令各列向量：
22、
23、其中，是一个极小值；
24、a3、利用拉格朗日乘子方法来求解式(3)所示的优化问题，其中拉格朗日函数为：
25、l＝tr[bss(bs)t]+λtr[bsd(bs)t]-tr[λ(bs(bs)t-im)] (6)
26、其中，λ∈rm×m是拉格朗日乘子的矩阵；
27、a4、对bs求导得：
28、
29、a5、先保持d不变，更新平稳投影矩阵bs：
30、bs(s+λd)＝λbs (8)
31、式中，这是一个关于s+λd的特征值分解问题，bs由s+λd的m个最小特征值对应的m个特征向量给出平稳投影矩阵
32、然后根据更新的bs并按照式子(4)计算d；
33、a6、重复迭代优化步骤a5，直至收敛。
34、进一步地，采用马氏距离构建所有平稳源的整体监测统计指标：...
【技术保护点】
1.一种针对非平稳工业过程的故障监测方法，其特征在于，包括：
2.根据权利要求1所述的针对非平稳工业过程的故障监测方法，其特征在于，利用稀疏平稳投影矩阵对工业过程的多维观测变量进行稀疏线性组合构建平稳源的表达式为：
3.根据权利要求1所述的针对非平稳工业过程的故障监测方法，其特征在于，所述采用平稳子空间分析方法并在优化问题中引入l2,p范数作为稀疏约束，得到优化问题为：
4.根据权利要求3所述的针对非平稳工业过程的故障监测方法，其特征在于，求解优化问题以得到平稳投影矩阵具体为：
5.根据权利要求3所述的针对非平稳工业过程的故障监测方法，其特征在于，采用马氏距离构建所有平稳源的整体监测统计指标：
6.根据权利要求3所述的针对非平稳工业过程的故障监测方法，其特征在于，故障诊断具体包括：
7.一种基于权利要求1所述故障监测方法的故障监测装置，其特征在于，包括：
8.根据权利要求7所述的故障监测装置，其特征在于，所述采用平稳子空间分析方法并在优化问题中引入l2,p范数作为稀疏约束，得到优化问题为：
...【技术特征摘要】
1.一种针对非平稳工业过程的故障监测方法，其特征在于，包括：
2.根据权利要求1所述的针对非平稳工业过程的故障监测方法，其特征在于，利用稀疏平稳投影矩阵对工业过程的多维观测变量进行稀疏线性组合构建平稳源的表达式为：
3.根据权利要求1所述的针对非平稳工业过程的故障监测方法，其特征在于，所述采用平稳子空间分析方法并在优化问题中引入l2,p范数作为稀疏约束，得到优化问题为：
4.根据权利要求3所述的针对非平稳工业过程的故障监测方法，其特征在于，求解优化问题以得到平稳投影矩阵具体为：
5.根据权利要求3所述的针对非平稳工业过程的故障监测方法，其特征在于，采用马氏距离构建所有平稳源的整体...
【专利技术属性】
技术研发人员：黄科科，李俊贤，阳春华，吴德浩，刘一顺，周灿，桂卫华，
申请(专利权)人：中南大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。