[发明专利]一种基于公共子空间分离的电熔镁炉运行故障检测方法

摘要

一种基于公共子空间分离的电熔镁炉运行故障检测方法，该方法通过核偏最小二乘回归的方法，将电熔镁炉运行过程获取的过程变量和质量变量进行了相关性处理，运用处理后的质量相关的过程数据通过局部切空间排列算法进行公共子空间的提取，对原数据空间进行了进一步划分，分为了质量相关公共子空间和质量相关特殊子空间，并建立相应的核主元分析监测模型，通过电熔镁炉多个运行模式共享的质量相关公共子空间和各运行模式的质量相关特殊子空间组合监测，可以在运行模式切换时减少故障的误报率，同时也提高了故障检测的灵敏性。

主权项

一种基于公共子空间分离的电熔镁炉运行故障检测方法，其特征在于：按如下步骤进行：步骤1：采集离线历史数据，按其属性形成数据集合；在电熔镁炉运行过程中，设定有M个运行模式，当电熔镁炉稳定运行在不同模式时分别采集M个模式下三相电压值U_A、U_B、U_C、三相电流值I_A、I_B、I_C、三个电极位置P_A、P_B、P_C及炉温T；得到电熔镁炉运行过程离线历史数据集Χ_m＝[x₁^m,…,x_N^m]^T∈(N×J)及电熔镁炉产品质量离线历史数据集Y_m＝[y₁^m,…,y_N^m]^T∈(N×J)，其中，m＝1,2,…,M，上标m代表第m个运行模式；下标N为第m个运行模式下数据采样的个数；J为变量个数；步骤2：对数据进行处理，建立电熔镁炉运行过程的监测模型；对电熔镁炉M个运行模式下得到的电熔镁炉运行过程离线历史数据集和电熔镁炉产品质量离线历史数据集进行处理，提取电熔镁炉M个运行模式的质量相关的公共信息，得到电熔镁炉M个运行模式共享的质量相关公共子空间；将电熔镁炉每个运行模式区分为共享的公共子空间以及各自独享的特殊子空间；利用核主元分析法(KPCA)，分别为电熔镁炉M个运行模式共享的公共子空间和每个模式的特殊子空间建立监测模型；步骤3：利用电熔镁炉运行过程的监测模型，在线监测电熔镁炉运行过程以检测电熔镁炉运行故障；利用电熔镁炉运行过程的监测模型，分别对电熔镁炉M个运行模式共享的质量相关公共子空间的霍特林(Hotelling‑T²)统计量、以及电熔镁炉每个运行模式的质量相关特殊子空间的霍特林(Hotelling‑T²)统计量和SPE统计量进行计算和监测；步骤3.1：在线获取电熔镁炉运行过程的新采样数据；步骤3.2：对电熔镁炉M个运行模式共享的质量相关公共子空间的霍特林(Hotelling‑T²)统计量、以及电熔镁炉每个运行模式的质量相关特殊子空间的霍特林(Hotelling‑T²)统计量和SPE统计量进行在线计算和监测；新采样数据的质量相关公共子空间的得分向量及其对应的霍特林(Hotelling‑T²)统计量，如式(5)所示：$\left\{\begin{array}{c}{\hat{t}}_{\mathrm{new}}^C={\hat{P}}_g^T\mathrm{\Φ}\left(x_{\mathrm{new}}\right)\\ {\hat{T}}_c^2={\left({\hat{t}}_{\mathrm{new}}^C\right)}^T{\left({\mathrm{\Λ}}^C\right)}^{-1}\left({\hat{t}}_{\mathrm{new}}^C\right)\end{array}\right.---\left(5\right)$式中，为新采样数据的质量相关公共子空间的得分向量；由公式(3)求得；ΛC为新采样数据的质量相关公共子空间与保留的主元数有关的协方差矩阵；表示新采样数据的质量相关公共子空间的霍特林(Hotelling‑T²)统计量；新采样数据的质量相关特殊子空间的得分向量以及相应的霍特林(Hotelling‑T²)统计量的计算公式，如式(6)所示：$\left\{\begin{array}{c}{\hat{t}}_{\mathrm{new},m}^S={\hat{P}}_m^T\mathrm{\Φ}\left(x_{\mathrm{new}}\right)\\ {\hat{T}}_{S,m}^2={\left({\hat{t}}_{\mathrm{new},m}^S\right)}^T{\left({\mathrm{\Λ}}_m^s\right)}^{-1}\left({\hat{t}}_{\mathrm{new},m}^S\right)\end{array}\right.---\left(6\right)$式中，为新采样数据对应的电熔镁炉运行模式的质量相关特殊子空间的得分向量；由公式(4)计算得到；是与新采样数据对应的电熔镁炉运行模式的质量相关特殊子空间的保留主元数有关的协方差矩阵；表示新采样数据对应的电熔镁炉运行模式的质量相关特殊子空间的霍特林(Hotelling‑T²)统计量；电熔镁炉新采样数据的质量相关特殊子空间的SPE统计量的计算，按式(7)进行：${\mathrm{SPE}}_{\mathrm{new},m}={\left({\hat{e}}_{\mathrm{new},m}^S\right)}^T{\hat{e}}_{\mathrm{new},m}^S---\left(7\right)$式中，为电熔镁炉运行模式m的质量相关特殊子空间的残差，SPE_new,m为新采样数据对应的电熔镁炉运行模式m的SPE统计量；步骤3.3：根据步骤3.2得到的计算结果，及质量相关公共子空间的霍特林(Hotelling‑T²)统计量的置信限、以及质量相关特殊子空间的霍特林(Hotelling‑T²)统计量的置信限和SPE统计量的置信限，检测电熔镁炉运行故障；步骤3.3.1：在线监测过程中，判断新采样数据的质量相关公共子空间的霍特林(Hotelling‑T²)统计量是否超出质量相关公共子空间的霍特林(Hotelling‑T²)统计量的置信限，是，则认为整个电熔镁炉运行过程发生故障，并执行步骤3.3.2；否，整个电熔镁炉运行过程未发生故障；步骤3.3.2：判断新采样数据的质量相关特殊子空间的霍特林(Hotelling‑T²)统计量和SPE统计量是否分别超出质量相关特殊子空间的霍特林(Hotelling‑T²)的置信限和SPE统计量的置信限，若二者之一超出其置信限,则判定该新采样数据所在的电熔镁炉运行模式发生故障；若二者均未超出其各自的置信限，则确定该新采样数据所在的电熔镁炉运行模式未发生故障。