一种基于变量重要性投影分块的过程监测方法

1.本发明属于故障检测技术领域，具体涉及一种基于变量重要性投影分块的多元指数加权移动平均核主成分分析过程监测方法。


背景技术：

2.现代工业过程中都增加了更多子流程或处理单元以追求更高的产品质量，提高生产效率。同时，为了采集更多过程相关的数据，每个单元内部都配备了许多高精度的工业传感器，这也导致数据维度成倍地增长。对于这类大规模复杂工业过程，直接采用所有变量建立单个模型进行监测这样的方式可能不再适合，数据高维性不仅会带来“维度灾难”，而且无关变量的引入还可能会降低模型的检测效果。针对这一问题，有研究人员提出采用变量分块或分布式监控框架，即将整个过程划分为几个有意义的块，然后构建子块模型进行监测。这种方式不仅考虑到过程中的局部信息，而且还有效简化了模型结构，降低了计算复杂度。近年来，关于分块线性监测模型相关的方法已有很多研究，然而，在实际大规模工业过程中，变量之间普遍存在较强的非线性关系，这些复杂的非线性结构以及高阶信息应在分块时给予更多的考虑。互信息源自于信息论，它可以从熵的角度定量地描述两个随机变量之间相互依赖的程度。互信息值越大，表示这两个变量的关系越密切，反之，则表示两个变量近似独立。虽然互信息对数据中的非线性结构和高阶信息描述的有效性已经得到证明，但忽略了由自变量组合所表示的潜变量对质量变量的影响。
3.主成分分析作为一个经典的多元统计方法在众多工业过程中得到广泛地应用，但难以处理数据中的非线性特征。因此，许多基于主成分分析的非线性方法被提出，其中核主成分分析由于其避免了非线性优化问题，使模型训练变得更容易。在非线性过程中，系统的运行状况可能会因为设备老化、设定点偏移或天气等原因发生微小的变化，这些故障随时间变化，通常在前期难以捕捉到。即使经过一段时间累积效应后在后期可以检测到该故障，但已经产生较大的检测时延。因此，及时准确地检测到这些微小故障是保证安全生产，减少企业经济损失的关键。为了解决这个问题，有学者提出将核主成分分析提取的核主成分得分和残差统计量作为多元指数加权移动平均的输入，以提高微小变化所导致故障的检测效果。因此，这种类似的混合过程监测方法应给予更多关注。


技术实现要素：

4.本发明针对现有技术中存在的问题，提供一种基于改进变量重要性投影分块的多元指数加权移动平均核主成分分析过程监测方法，所述方法包括以下步骤：
5.s1.数据预处理：对输入数据x及输出数据y进行标准化，将其转化为均值为0，方差为1的无量纲数据；将训练集数据用于模型的构建及训练，测试集数据用于对模型的性能进行评估。
6.s2.构建改进的变量重要性投影指标：在变量重要性投影的基础上引入互信息对潜变量和质量变量之间的关系进行定量描述，并构建改进的变量重要性投影指标。
7.s3.构建动态模型：根据改进的变量重要性投影，将过程变量划分到质量相关空间和质量独立空间，同时对不同空间中的数据变量进行多元指数加权移动平均预处理，实现动态模型构建。
8.s4.构建多元指数加权移动平均核主成分分析模型：将通过多元指数加权移动平均预处理后的数据作为核主成分分析的输入，构建多元指数加权移动平均核主成分分析模型，以此来提高模型对数据中的动态信息的捕获能力及微小故障的敏感性。
9.s5.对数据进行反标准化，并完成模型性能的评估：将测试集数据带入训练后的模型进行测试，根据故障检测率这个指标完成模型性能的评估。
10.该方法的优势在于充分考虑了变量之间的非线性关系，在变量重要性投影基础上进一步引入互信息构建改进的变量重要性投影指标。基于此，将整个过程变量进行分块，适当的变量划分不仅能够降低模型计算复杂度，还能够挖掘过程中的局部信息。最后，分别在两个子空间内构建多元指数加权移动平均核主成分分析模型来提高对非线性过程中微小变化的敏感性。因此，在大规模复杂工业过程中，该模型可靠性高，能够充分地挖掘局部信息，具备较高的故障检测率。
11.采用上述方案后，本发明与现有技术相比具有如下效果：
12.与现有技术相比，基于改进的变量重要性投影分块的多元指数加权移动平均核主成分分析过程监测方法的提出，克服了大规模工业过程数据中非线性和时变性特征，解决了故障发生初期检测难的问题，有效降低了模型计算复杂度。利用改进的变量重要性投影指标对过程变量分块，使其不仅避免了传统分块方法在确定分块个数时的主观性，还增强了对监测结果的可解释性。因此，对工程师来说更容易对过程运行状态和工艺优化进行把控，有助于保障生产安全，减少企业经济损失。
附图说明
13.图1是田纳西-伊斯曼过程的工艺及控制流程图。
14.图2是基于改进变量重要性投影的变量分块结构图。
15.图3是基于改进变量重要性投影的变量分块结果。
16.图4是对于故障5基于改进变量重要性投影分块的多元指数加权移动平均核主成分分析的监测结果。
具体实施方式
17.下面将对本发明进行更清楚、完整的描述，显然，所描述的实例仅仅是本发明的一部分实例，而不是全部的实施例。
18.本发明提供的对大规模复杂工业过程监测所采取的技术方案是：
19.s1.数据预处理：根据公式(1)完成输入数据x及输出数据y的标准化处理；将其转化为均值为0，方差为1的无量纲数据；将训练集数据用于模型的构建及训练，测试集数据用于对模型的性能进行评估。
20.s2.构建改进的变量重要性投影指标：在变量重要性投影的基础上引入互信息对潜变量和质量变量之间的关系进行定量描述，并构建改进的变量重要性投影指标。
21.s3.构建动态模型：根据改进的变量重要性投影，将过程变量划分到质量相关空间
和质量独立空间，同时对不同空间中的数据变量进行多元指数加权移动平均预处理，实现动态模型构建。
22.s4.构建多元指数加权移动平均核主成分分析模型：将通过多元指数加权移动平均预处理后的数据作为核主成分分析的输入，构建多元指数加权移动平均核主成分分析模型，以此来提高模型对数据中的动态信息的捕获能力及微小故障的敏感性。
23.s5.对数据进行反标准化，并完成模型性能的评估：将测试集数据带入训练后的模型进行测试，根据故障检测率这个指标完成模型性能的评估。
24.步骤s1中，将数据转化为均值为0，方差为1的公式为：
[0025][0026]
式中，x
*
为原数据，x为标准化后的数据，μ和σ分别为所有训练集样本的均值和方差。
[0027]
步骤s2中，改进的变量重要性投影指标的构建如下：
[0028]
s21：使用偏最小二乘模型对标准化后的输入数据与输出数据作如下分解：
[0029][0030]
其中，是得分矩阵；和分别为x和y的负载矩阵；e和f分别为输入和输出矩阵的残差；l为潜变量的个数，通常由交叉验证来确定；
[0031]
s22：采用潜变量和质量变量之间的互信息值来替代原始通过计算每个成分所解释的方差的部分，即：
[0032][0033]
式中，为估计两个变量互信息值的函数；ti第i个主成分得分；yk为第k个质量变量。
[0034]
s23：计算第j个变量的改进的变量重要性投影(ivip)值：
[0035][0036]
步骤s3中，利用改进的变量重要性投影将过程变量进行分块，并构建动态模型：
[0037]
s31：通过计算每个过程变量的改进的变量重要性投影，将其划分到质量相关空间和质量独立空间：
[0038]
x＝[x
qr
，x
in
]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0039]
s32：分别对不同空间的数据变量进行预处理，构建动态模型：
[0040]
zi＝βxi (1-β)z
i-1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0041]
式中，0＜β≤1为平滑参数，控制着检测时延和灵敏度之间的平衡。
[0042]
步骤s4中，构建多元指数加权移动平均核主成分分析模型：
[0043]
s41：对处理后的数据利用核技巧构建gram矩阵
[0044]kij
＝k(zi，zj)＝《φ(zi)，φ(zj)》
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0045]
其中，k
ij
为gram矩阵k中的元素，核函数采用高斯核，定义为k(x，y)＝exp(-||x-y||/(2
×
σ2))；σ为可调节的核参数，通常由经验或者网格搜索确定。
[0046]
s42：对gram矩阵k进行中心化：
[0047][0048]
其中
[0049][0050]
s43：对进行特征分解：
[0051][0052]
式中，λ和α分别为的特征值和特征向量。
[0053]
s44：根据累计方差贡献率确定核主成分个数l；
[0054]
s45：计算非线性主成分得分k＝1，2，...，p，并构建t2和q统计量：
[0055][0056][0057]
s46：分别计算质量相关空间和质量独立空间的两个统计量控制限：
[0058][0059][0060]
g＝b/2a
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(15)
[0061]
h＝2a2/b
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(16)
[0062]
其中，f表示f分布；χ2表示卡方分布；a和b分别为q统计量的均值和方差。
[0063]
s47：统计量任何一个超出置信度的话，则认为发生了故障。否则认为该过程正常。
[0064]
步骤s5中，对数据进行反标准化，并完成模型性能的评估：
[0065]
将测试集数据带入训练后的模型进行测试，根据故障检测率(fdr)这个指标完成模型性能的评估。fdr的计算公式如下：
[0066][0067]
实施例1：
[0068]
以下结合田纳西-伊斯曼过程来说明本发明的有效性，它是美国伊斯曼化学品公司在1993年基于真实的化工过程所建立的一个仿真平台，旨在于对过程控制与监测方法的可行性和有效性的验证提供一个真实的案例，已获得学术界的高度认可，并广泛应用在控制、优化以及过程监测等领域。田纳西-伊斯曼过程的工艺及控制流程图如图1所示，主要由五个单元组成，包括反应器、冷凝器、压缩机、分离器和汽提塔。传感器记录有41个测量变量和12个操纵变量，均包含高斯噪声。41个测量变量中又可分为22个连续过程变量，和19个成
分变量。该仿真数据集预设了21种故障，其中有16种已知故障，5种未知故障。表1对这些故障进行详细描述。
[0069]
表1.田纳西-伊斯曼过程故障描述
[0070][0071]
下面对本发明做进一步详述：
[0072]
第一步：选用正常状态下的过程变量(1-22)和操纵变量(1-11)为输入变量，最终产物作为质量变量用于模型训练，并对其进行数据预处理。测试集则采用21种预设故障数据来评估模型的性能。
[0073]
第二步：在变量重要性投影基础上引入互信息，构建改进的变量重要性投影指标。通过式(4)计算每个过程变量的改进的变量重要性投影，如图2所示。
[0074]
第三步：根据改进的变量重要性投影，将大于阈值的过程变量划分到质量相关空间，其余则分到质量独立空间，如图3所示。图中的横线为设定的阈值。对不同空间中的数据变量进行多元指数加权移动平均预处理，实现动态模型构建。
[0075]
第四步：通过多元指数加权移动平均预处理后的数据作为核主成分分析的输入，构建多元指数加权移动平均核主成分分析模型，以此来提高模型对数据中的动态信息的捕获能力及微小故障的敏感性。
[0076]
第五步：将测试集数据进行反标准化，完成模型性能评估。本发明所提出的基于改进的变量重要性投影分块的多元指数加权移动平均核主成分分析方法不仅能够监测整个过程中的系统运行状况，而且对于工程师来说，这些与质量相关的变量被单独划分出来，加强了监测结果对质量变量的可解释性。图4(a)和(b)为所提出方法对于故障5分别在质量相关空间和质量独立空间的监测结果。相比于传统的全局模型，基于改进的变量重要性投影分块的多元指数加权移动平均核主成分分析过程监测方法能够有效检测到非线性过程中
的微小故障，并且工程师可以根据异常情况及时对过程进行调节和维修，避免安全事故和更大的损失。
[0077]
以上描述了本发明的基本原理、主要特征及本发明的优势所在。以上所述，仅为本发明较佳的具体实施方式，本发明的保护范围不局限于此，本行业的技术人员在本发明展示的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该由所附的权利要求书及其等效物界定。