CN109164794A - 基于偏f值selm的多变量工业过程故障分类方法 - Google Patents

Abstract

本发明涉及一种基于偏F值SELM的堆栈极限学习机故障分类方法，首先对训练数据集进行标准化处理；其次，求出所有故障类型对于每一个变量的偏F值，根据所有变量的偏F值得出所有变量的权值，进一步对训练以及测试数据进行加权；通过训练数据集利用堆栈极限学习机建模方法构建FSELM分类模型；在此基础上，对测试数据集进行标准化处理；最后通过测试数据集验证该分类模型的好坏。本发明能够明显提升故障诊断性能，提高故障分类的准确率。

Description

基于偏F值SELM的多变量工业过程故障分类方法

技术领域

本发明属于工业过程故障诊断技术领域，涉及一种利用基于偏F值堆栈极限学习机(简称：FSELM)对多变量工业过程故障进行分类的方法。

背景技术

分布式计算机控制技术的进步极大地推动了工业过程生产规模的扩大。对于复杂庞大的工业生产过程而言，高效的故障诊断技术不仅有助于提高装置生产的稳定性，而且也有助于保障工作人员的人身安全。故障诊断技术的核心在于故障出现后准确的识别故障类型，因此故障分类技术对于现场工程师排除故障具有重要意义。

极限学习机(Extreme Learning Machine,简称：ELM)是由Guang-Bin Huang等人于2006年提出的一种单隐含层前馈神经网络学习算法，该算法训练过程中无需多次迭代调整，只需合理选择隐含层神经元的个数、选择合适的激活函数，便可获得唯一的最优解。因其具有训练速度快、分类精度高、泛化能力强等优点，目前已经在时间序列预测、图像处理、故障分类等多个领域得到广泛应用。为处理复杂多变量工业数据，研究人员进一步提出了堆栈极限学习机(简称：SELM)，在堆栈学习框架下构建多层ELM，以实现对数据特征信息的深度挖掘和提取。然而，传统的SELM在构建模式分类器时将所有的过程变量同等对待，赋予相同的权值，使得部分有助于故障分类的重要变量信息被整体变量信息掩盖，从而降低了分类效果。因此，如何判断过程变量的重要性差异并赋予不同权值，成为SELM故障分类器构建中的挑战课题。

发明内容

本发明针对SELM故障分类器未充分挖掘局部重要变量信息导致分类效果降低的问题，提供一种基于偏F值SELM的多变量工业过程故障分类方法。该方法能够准确提取过程数据中有助于故障分类的部分重要变量，降低冗余变量对故障分类的影响，提高故障分类准确度。

为了达到上述目的，本发明提供了一种基于偏F值SELM的多变量工业过程故障分类方法，含有以下步骤：

(一)采集工业过程正常操作工况数据集X_o和C类故障工况数据集{X₁,X₂,…X_C}作为训练数据集，并使用正常操作工况数据集的均值mX_o和标准差sX_o对训练数据集进行标准化处理得到标准化后的正常操作工况数据集和C类故障工况数据集

(二)针对每个故障工况数据集，计算过程变量的偏F值F(i,j),其中F(i,j)表示第i类故障工况数据集中第j个变量的偏F值；

(三)对所有故障工况数据集得到的偏F值求和，得到第j个变量在所有故障工况情形下的总体偏F值Fs(j)，并根据各个变量偏F值在所有变量中的占比p(j)计算每个变量在故障分类过程中的重要性权值w(j)；

(四)利用重要性权值w(j)对故障工况数据集进行加权处理；

(五)利用堆栈极限学习机建模方法构建FSELM分类模型，FSELM分类模型表示为：y＝f(x,W⁽¹⁾,…,W^(P),V⁽¹⁾,…,V^(P),B⁽¹⁾,…,B^(P))，其中，x表示输入向量，即待分类的过程样本，y表示输出向量，为样本的类别编码，f(·)表示一个多层ELM堆栈网络函数，W⁽ⁱ⁾,i＝1,2,...,P表示输入层权重参数，V⁽ⁱ⁾,i＝1,2,...,P表示线性变换阵，B⁽ⁱ⁾,i＝1,2,...,P表示输出权重参数；

(六)采集测试数据集X_t，并利用正常操作工况数据集的均值和方差进行标准化处理，得到标准化后的测试数据集

(七)利用重要性权值w(j)对测试数据集进行加权处理；

(八)将加权后的测试数据集输入FSELM分类模型，根据FSELM分类模型输出确定故障类型。

进一步的，利用正常操作工况数据集的均值mX_o和标准差sX_o通过公式(1)对正常操作工况数据集X_o和C类故障工况数据集{X₁,X₂,…X_C}进行标准化处理，公式(1)的表达式为：

正常操作工况数据集X_o和C类故障工况数据集{X₁,X₂,…X_C}由上述公式(1)标准化处理后，即得到标准化后的正常操作工况数据集和C类故障工况数据集

进一步的，步骤(二)中，针对故障工况数据集通过计算公式(2)计算第j变量的偏F值F(i,j)，计算公式(2)表示为：

式中，v＝n₀+n_i-2，n₀为正常操作工况数据集中的样本数目，n_i为故障工况数据集中的样本数目，m为正常操作工况数据集和故障工况数据集两个数据集的变量个数，T_i ²为考虑所有过程变量时的T²统计量，表示为：

式中，v＝n₀+n_i-2，n₀为正常操作工况数据集中的样本数目，n_i为故障工况数据集中的样本数目，m为正常操作工况数据集和故障工况数据集两个数据集的变量个数，T_i ²为考虑所有过程变量时的T²统计量，表示为：

式中，为b₀去掉第j个变量后的向量，为b_i去掉第j个变量后的向量；为S_i去掉第j行和第j列后的矩阵。

进一步的，步骤(三)中，对所有故障数据集得到的偏F值求和，得到第j个变量在所有故障情形下的总体偏F值Fs(j)，总体偏F值Fs(j)表示为：

通过公式(6)计算第j个变量的总体偏F值Fs(j)在所有变量中的占比p(j)，公式(6)表示为：

根据各个变量偏F值占比p(j)通过公式(7)计算每个变量在故障分类过程中的重要性权值w(j)，公式(7)表示为：

式中，m正常操作工况数据集和故障工况数据集两个数据集的变量个数。

进一步的，步骤(四)中，利用重要性权值w(j)通过公式(8)对故障工况数据集进行加权处理，公式(8)表达为：

式中，为加权处理后的故障工况数据集。

进一步的，步骤(五)中，利用堆栈极限学习机建模方法构建FSELM分类模型的步骤为：

(1)构建训练数据集并根据数据类别编制相应的输出矩阵Y；设置堆栈极限学习机的超参数P、L、L1；

(2)随机初始化输入层权重参数W⁽ⁱ⁾,i＝1,2,...,P，其中，W⁽¹⁾为L×m维矩阵，W⁽²⁾，…W^(P)为(L-L1)×m维矩阵；初始化隐含层序号p＝1；

(3)计算第p个隐含层输出矩阵H^(p)＝g(X(W^(p))^T)，其中，g()表示sigmoid激活函数；

(4)计算输出权重矩阵B^(p)＝(H^(p))^*Y，其中，B^(p)为L×C维矩阵,^*表示MP广义逆运算；

(5)对输出权重矩阵B^(p)进行主元分析，得到L×L1维的线性变换阵V^(p)；

(6)求出第p+1个隐含层矩阵H^(p+1)＝[H^(p)V^(p)g(X(W^(p+1))^T)]；

(7)令p＝p+1，如果p<P，将新的隐含层矩阵H^(p+1)代替H^(p)，重复步骤(4)-(6)，进一步计算B^(p+1)、V^(p+1)；否则达到堆栈层数，FSELM模型构建完成；构建的FSELM分类模型表示为：

y＝f(x,W⁽¹⁾,…,W^(P),V⁽¹⁾,...,V^(P),B⁽¹⁾,...,B^(P)) (9)。

进一步的，步骤(六)中，利用正常操作工况数据集的均值mX_o和标准差sX_o通过公式

(10)对测试数据集X_t进行标准化处理，公式(10)表示为：

测试数据集X_t由上述公式(10)标准化处理后，即得到标准化后的测试数据集

进一步的，步骤(七)中，利用重要性权值w(j)通过公式(11)对测试数据集进行加权处理，公式(11)表示为：

式中，为加权处理后的测试数据集。

进一步的，步骤(八)中，将加权后的测试数据集输入FSELM分类模型，根据FSELM分类模型输出得出测试数据集的故障类型矩阵故障类型矩阵表示为：

与现有技术相比，本发明的有益效果在于：

本发明提供的基于偏F值SELM的多变量工业过程故障分类方法，对采集的训练数据集进行标准化处理，求出所有故障类型对于每一个变量的偏F值，根据所有变量的偏F值得出所有变量的权值，进一步对训练以及测试数据进行加权，通过训练数据集利用堆栈极限学习机建模方法构建FSELM分类模型，通过测试数据集验证该分类模型的好坏，并根据FSELM分类模型对故障类型进行分类。本发明挖掘有助于故障分类的过程变量并给予较大权重，准确提取过程数据中有助于故障分类的部分重要变量，降低冗余变量对故障分类的影响，能够有效地提高故障分类的准确率。

附图说明

图1为本发明所述基于偏F值SELM的多变量工业过程故障分类方法的流程图；

图2为本发明实施例所述田纳西-伊斯曼(简称：TE)过程的结构图；

图3a为本发明实施例所述故障1各个变量的偏F值；

图3b为本发明实施例所述故障2各个变量的偏F值；

图3c为本发明实施例所述故障3各个变量的偏F值；

图3d为本发明实施例所述故障4各个变量的偏F值。

具体实施方式

下面，通过示例性的实施方式对本发明进行具体描述。然而应当理解，在没有进一步叙述的情况下，一个实施方式中的元件、结构和特征也可以有益地结合到其他实施方式中。

参见图1，本发明揭示了一种基于偏F值SELM的多变量工业过程故障分类方法，含有以下步骤：

(一)采集工业过程正常操作工况数据集X_o和C类故障工况数据集{X₁,X₂,…X_C}作为训练数据集，并使用正常操作工况数据集的均值mX_o和标准差sX_o通过公式(1)对训练数据集进行标准化处理，公式(1)的表达式为：

训练数据集由上述公式(1)标准化处理后，即得到标准化后的正常操作工况数据集和C类故障工况数据集

(二)针对故障工况数据集通过计算公式(2)计算第j变量的偏F值F(i,j)，其中F(i,j)表示第i类故障工况数据集中第j个变量的偏F值。计算公式(2)表示为：

式中，v＝n₀+n_i-2，n₀为正常操作工况数据集中的样本数目，n_i为故障工况数据集中的样本数目，m为正常操作工况数据集和故障工况数据集两个数据集的变量个数，T_i ²为考虑所有过程变量时的T²统计量，表示为：

式中，v＝n₀+n_i-2，n₀为正常操作工况数据集中的样本数目，n_i为故障工况数据集中的样本数目，m为正常操作工况数据集和故障工况数据集两个数据集的变量个数，T_i ²为考虑所有过程变量时的T²统计量，表示为：

式中，为b₀去掉第j个变量后的向量，为b_i去掉第j个变量后的向量；为S_i去掉第j行和第j列后的矩阵。

若针对某类故障数据的偏F值越大，说明此变量对该类故障进行分类时的贡献越大；反之，针对某类故障数据的偏F值越小，说明此变量对该类故障进行分类的贡献越小。

(三)对所有故障数据集得到的偏F值求和，得到第j个变量在所有故障情形下的总体偏F值Fs(j)，总体偏F值Fs(j)表示为：

通过公式(6)计算第j个变量的总体偏F值Fs(j)在所有变量中的占比p(j)，公式(6)表示为：

根据各个变量偏F值占比p(j)通过公式(7)计算每个变量在故障分类过程中的重要性权值w(j)，公式(7)表示为：

(四)利用重要性权值w(j)通过公式(8)对故障工况数据集进行加权处理，公式(8)表达为：

式中，为加权处理后的故障工况数据集。

(五)利用堆栈极限学习机建模方法构建FSELM分类模型，FSELM分类模型表示为：y＝f(x,W⁽¹⁾,...,W^(P),V⁽¹⁾,…,V^(P),B⁽¹⁾,…,B^(P))，其中，x表示输入向量，即待分类的过程样本，y表示输出向量，为样本的类别编码，f(·)表示一个多层ELM堆栈网络函数，W⁽ⁱ⁾,i＝1,2,...,P表示输入层权重参数，V⁽ⁱ⁾,i＝1,2,...,P表示线性变换阵，B⁽ⁱ⁾,i＝1,2,...,P表示输出权重参数。构建FSELM分类模型的具体步骤为：

(1)构建训练数据集并根据数据类别编制相应的输出矩阵Y；设置堆栈极限学习机的超参数P、L、L1；

(2)随机初始化输入层权重参数W⁽ⁱ⁾,i＝1,2,...,P，其中，W⁽¹⁾为L×m维矩阵，W⁽²⁾，…W^(P)为(L-L1)×m维矩阵；初始化隐含层序号p＝1；

(3)计算第p个隐含层输出矩阵H^(p)＝g(X(W^(p))^T)，其中，g()表示sigmoid激活函数；

(4)计算输出权重矩阵B^(p)＝(H^(p))^*Y，其中，B^(p)为L×C维矩阵,^*表示MP广义逆运算；

(5)对输出权重矩阵B^(p)进行主元分析，得到L×L1维的线性变换阵V^(p)；

(6)求出第p+1个隐含层矩阵H^(p+1)＝[H^(p)V^(p)g(X(W^(p+1))^T)]；

(7)令p＝p+1，如果p<P，将新的隐含层矩阵H^(p+1)代替H^(p)，重复步骤(4)-(6)，进一步计算B^(p+1)、V^(p+1)；否则达到堆栈层数，FSELM模型构建完成；构建的FSELM分类模型表示为：

y＝f(x,W⁽¹⁾,…,W^(P),V⁽¹⁾,…,V^(P),B⁽¹⁾,…,B^(P)) (9)。

(六)采集测试数据集X_t，利用正常操作工况数据集的均值mX_o和标准差sX_o通过公式(10)对测试数据集X_t进行标准化处理，公式(10)表示为：

测试数据集X_t由上述公式(10)标准化处理后，即得到标准化后的测试数据集

(七)利用重要性权值w(j)通过公式(11)对测试数据集进行加权处理，公式(11)表示为：

式中，为加权处理后的测试数据集。

(八)将加权后的测试数据集输入FSELM分类模型，根据FSELM分类模型输出得出测试数据集的故障类型矩阵故障类型矩阵表示为：

上述方法中，步骤(一)至(五)为离线建模阶段，步骤(六)至(八)为在线测试阶段。

本发明上述故障分类方法对训练数据集进行标准化处理；求出所有故障类型对于每一个变量的偏F值，根据所有变量的偏F值得出所有变量的权值，进一步对训练以及测试数据进行加权；在此基础上，对测试数据集进行标准化处理；通过训练数据建立基于堆栈极限学习机的故障分类模型。本发明挖掘有助于故障分类的过程变量并给予较大权重，准确提取过程数据中有助于故障分类的部分重要变量，降低冗余变量对故障分类的影响，能够有效地提高故障分类的准确率。

为了能更清楚地说明本发明上述故障分类方法的有益效果，以下结合具体实施例对本发明上述故障分类方法做出进一步说明。

实施例：田纳西-伊斯曼(以下简称：TE)过程是由美国伊斯曼化学公司的Downs和Vogel根据一个实际的化工过程建立的实验平台，现被广泛用于验证控制算法和过程监控方法的优劣。参见图2，TE过程主要由五个单元组成，包括反应器、产品冷凝器、气液分离器、循环压缩机和汽提塔组成。TE过程共53个变量，其中包括22个连续测量变量、19个成分变量和12个操作变量。其中34个变量包括22个连续测量变量和12个操作变量，如表1所示。

表1

实验中采用TE过程故障的4类故障数据，如表2所示，用来验证各方法的分类准确率。

表2

在本实施中，主要是对TE过程的故障1、故障2、故障3、故障4这4类故障类型进行分类，其中训练集中每类故障各300个样本，测试集中每类故障各200个样本。

图3a表示的是故障1中各个变量的偏F值，图3b表示的是故障2中各个变量的偏F值，图3c表示的是故障3中各个变量的偏F值，图3d表示的是故障4中各个变量的偏F值。根据各个变量的偏F值和各个变量的偏F值在所有变量中的占比得出其中的15个变量的权值为1，即变量1、44、50、40、9、17、52、21、38、18、35、41、34、29、10。

分别利用极限学习机(简称：ELM)、堆栈极限学习机(简称：SELM)、本发明基于偏F值SELM的多变量工业过程故障分类方法(简称：FSELM)对以上4类TE故障进行分类。故障分类准确率结果如表3所示。

表3

	故障1	故障2	故障3	故障4	平均值
						ELM	0.6662	0.5858	0.5458	0.9160	0.6785
SELM	0.8598	0.6110	0.5558	0.9432	0.7425
						FSELM	0.9925	0.7180	0.6535	0.9655	0.8324

由上述表3可知，本发明提出的基于偏F值SELM故障分类方法能够明显的提高分类的准确率。

以上所举实施例仅用为方便举例说明本发明，并非对本发明保护范围的限制，在本发明所述技术方案范畴，所属技术领域的技术人员所作各种简单变形与修饰，均应包含在以上申请专利范围中。

Claims (9)

1.一种基于偏F值SELM的多变量工业过程故障分类方法，含有以下步骤：

(一)采集工业过程正常操作工况数据集X_o和C类故障工况数据集{X₁,X₂,…X_C}作为训练数据集，并使用正常操作工况数据集的均值mX_o和标准差sX_o对训练数据集进行标准化处理，得到标准化后的正常操作工况数据集和C类故障工况数据集

(二)针对每个故障工况数据集，计算过程变量的偏F值F(i,j),其中F(i,j)表示第i类故障工况数据集中第j个变量的偏F值；

(三)对所有故障工况数据集得到的偏F值求和，得到第j个变量在所有故障工况情形下的总体偏F值Fs(j)，并根据各个变量偏F值在所有变量中的占比p(j)计算每个变量在故障分类过程中的重要性权值w(j)；

(四)利用重要性权值w(j)对故障工况数据集进行加权处理；

(五)利用堆栈极限学习机建模方法构建FSELM分类模型，FSELM分类模型表示为：y＝f(x,W⁽¹⁾,…,W^(P),V⁽¹⁾,…,V^(P),B⁽¹⁾,…,B^(P))，其中，x表示输入向量，即待分类的过程样本，y表示输出向量，为样本的类别编码，f(·)表示一个多层ELM堆栈网络函数，W⁽ⁱ⁾,i＝1,2,...,P表示输入层权重参数，V⁽ⁱ⁾,i＝1,2,...,P表示线性变换阵，B⁽ⁱ⁾,i＝1,2,...,P表示输出权重参数；

(六)采集测试数据集X_t，并利用正常操作工况数据集的均值和方差进行标准化处理，得到标准化后的测试数据集

(七)利用重要性权值w(j)对测试数据集进行加权处理；

(八)将加权后的测试数据集输入FSELM分类模型，根据FSELM分类模型输出确定故障类型。

2.如权利要求1所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(一)中，利用正常操作工况数据集的均值mX_o和标准差sX_o通过公式(1)对训练数据集进行标准化处理，公式(1)的表达式为：

训练数据集由上述公式(1)标准化处理后，即得到标准化后的正常操作工况数据集和C类故障工况数据集

3.如权利要求2所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(二)中，针对故障工况数据集通过计算公式(2)计算第j变量的偏F值F(i,j)，计算公式(2)表示为：

式中，v＝n₀+n_i-2，n₀为正常操作工况数据集中的样本数目，n_i为故障工况数据集中的样本数目，m为正常操作工况数据集和故障工况数据集两个数据集的变量个数，T_i ²为考虑所有过程变量时的T²统计量，表示为：

式中，b₀为正常类数据的均值向量；b_i是第i类故障数据的均值向量；S_i为组合矩阵的协方差矩阵；为剔除变量j之后的T_i ²值，即：

式中，为b₀去掉第j个变量后的向量，为b_i去掉第j个变量后的向量；为S_i去掉第j行和第j列后的矩阵。

4.如权利要求3所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(三)中，对所有故障数据集得到的偏F值求和，得到第j个变量在所有故障情形下的总体偏F值Fs(j)，总体偏F值Fs(j)表示为：

通过公式(6)计算第j个变量的总体偏F值Fs(j)在所有变量中的占比p(j)，公式(6)表示为：

根据各个变量偏F值占比p(j)通过公式(7)计算每个变量在故障分类过程中的重要性权值w(j)，公式(7)表示为：

5.如权利要求4所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(四)中，利用重要性权值w(j)通过公式(8)对故障工况数据集进行加权处理，公式(8)表达为：

式中，为加权处理后的故障工况数据集。

6.如权利要求5所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(五)中，利用堆栈极限学习机建模方法构建FSELM分类模型的步骤为：

(1)构建训练数据集并根据数据类别编制相应的输出矩阵Y；设置堆栈极限学习机的超参数P、L、L1；

(2)随机初始化输入层权重参数W⁽ⁱ⁾,i＝1,2,...,P，其中，W⁽¹⁾为L×m维矩阵，W⁽²⁾，…W^(P)为(L-L1)×m维矩阵；初始化隐含层序号p＝1；

(3)计算第p个隐含层输出矩阵H^(p)＝g(X(W^(p))^T)，其中，g()表示sigmoid激活函数；

(4)计算输出权重矩阵B^(p)＝(H^(p))^*Y，其中，B^(p)为L×C维矩阵,^*表示MP广义逆运算；

(5)对输出权重矩阵B^(p)进行主元分析，得到L×L1维的线性变换阵V^(p)；

(6)求出第p+1个隐含层矩阵H^(p+1)＝[H^(p)V^(p) g(X(W^(p+1))^T)]；

(7)令p＝p+1，如果p<P，将新的隐含层矩阵H^(p+1)代替H^(p)，重复步骤(4)-(6)，进一步计算B^(p+1)、V^(p+1)；否则达到堆栈层数，FSELM模型构建完成；构建的FSELM分类模型表示为：

y＝f(x,W⁽¹⁾,…,W^(P),V⁽¹⁾,…,V^(P),B⁽¹⁾,…,B^(P)) (9)。

7.如权利要求6所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(六)中，利用正常操作工况数据集的均值mX_o和标准差sX_o通过公式(10)对测试数据集X_t进行标准化处理，公式(10)表示为：

测试数据集X_t由上述公式(10)标准化处理后，即得到标准化后的测试数据集

8.如权利要求7所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(七)中，利用重要性权值w(j)通过公式(11)对测试数据集进行加权处理，公式(11)表示为：

式中，为加权处理后的测试数据集。

9.如权利要求8所述的基于偏F值SELM的多变量工业过程故障分类方法，其特征在于，步骤(八)中，将加权后的测试数据集输入FSELM分类模型，根据FSELM分类模型输出得出测试数据集的故障类型矩阵故障类型矩阵表示为：