Comprehensive Fault Diagnosis of Shaft Furnace Roasting Processes Using Simplified Concurrent Projection to Latent Structures
-
摘要: 竖炉焙烧过程因运行条件异常变化或操作不当会造成上火、冒火、过还原和欠还原等运行故障.这些故障直接影响过程运行安全和产品质量(比如,磁选管回收率),但难以采用基于模型和基于知识的方法建模故障与产品质量的关系,以及诊断故障变量.针对上述问题,本文提出数据驱动的基于并发潜结构映射(Concurrent projection to latent structures,CPLS)的竖炉焙烧过程综合故障诊断方法.首先,将并发潜结构映射分解的过程变量共有子空间与残差空间精简合并来建立磁选管回收率相关的过程变化空间,提出基于精简并发潜结构映射模型的竖炉焙烧过程综合监控方法;接下来,定义相应的重构贡献图并与竖炉焙烧过程相结合,提出CPLS精简重构贡献方法用于竖炉焙烧过程故障变量诊断;最后,利用竖炉焙烧过程半实物仿真平台采集的数据进行实验研究,结果表明所提方法不仅可以诊断出质量相关的故障,而且可诊断出回路设定值之外的故障变量.Abstract: Operational faults of shaft furnace roasting processes can appear when operational conditions change abnormally or operators do not react properly or timely. Typical operational faults, including fire-emitting, flame-out, under-reduction and over-reduction, are highly related to process safety and product quality, e.g., magnetic tube recovery rate (MTRR). Fault diagnosis of shaft furnace roasting processes deserves more attentions. However, it is difficult to apply model-based or knowledge-based fault diagnosis methods. In particular, it is difficult to model the relations between fault and product quality. In this paper data-driven concurrent projection to latent structures (CPLS) based fault diagnosis is developed for shaft furnace roasting processes. First, a CPLS based comprehensive monitoring method for shaft furnace roasting processes is proposed by combining co-variation and residual of process spaces of concurrent projection to latent structures into a simplified MTRR-relevant process-variation space. Secondly, a corresponding simplified reconstruction-based contribution method is proposed and used to pinpoint the faulty variable. Finally, the proposed methods are verified using the data collected from a hardware-in-loop simulation platform. The results demonstrate that the quality-relevant faults as well as faulty variables are successfully diagnosed.
-
竖炉焙烧过程是选矿生产的关键工序, 用于对原料赤铁矿高温焙烧, 将弱磁性矿物转化成强磁性矿物, 作为下游磁选工序的原料.竖炉焙烧原料矿石的大小、种类、成分变化频繁, 在运行条件异常变化或操作不当时, 可发生上火和冒火等与运行安全密切相关的故障, 也可发生过还原和欠还原等与质量指标(即磁选管回收率)下降密切相关的故障, 有必要有效监控与诊断.
竖炉焙烧过程生产条件变化频繁、物理和化学反应机理复杂, 一方面难以依据物理和化学规律建立过程的精确数学模型; 另一方面, 竖炉焙烧过程的运行故障不是传统的传感器和执行机构故障, 难以建立故障模型.因此, 以过程模型与故障模型为基础的基于模型方法难以用于竖炉焙烧过程运行故障诊断[1].
近年来, 学者提出基于知识的竖炉焙烧过程故障诊断方法.比如, Chai等[2]、吴峰华等[3]建立了基于规则推理的竖炉焙烧过程故障诊断专家系统, 对过程变量的观测值利用产生式规则进行判断与推理, 以规则推理的结论作为故障诊断结论; 严爱军等[4]提出基于案例推理的竖炉焙烧过程故障诊断方法, 将新故障与案例库中的历史故障进行比较, 得到以概率形式表示的故障诊断结论.但上述方法存在两个问题: 1)只能诊断出故障发生, 但不能诊断故障原因, 特别是无法诊断回路设定值之外的故障原因变量; 2)由于竖炉焙烧过程的过程变量相关关系复杂, 难以获得用于规则推理及案例推理所需的过程先验知识, 方法实施与维护难度大.
针对上述问题, Lu等[5]结合数据驱动故障诊断方法无需过程精确数学模型、降维可视化以及实施维护方便等优点[6-8], 利用历史正常数据, 采用主元分析(Principal component analysis, PCA)技术建立竖炉焙烧过程的多元统计模型与故障监控指标, 在此基础上进行竖炉焙烧过程故障诊断.但该方法属于单一层面的潜结构建模与故障诊断, 将质量变量等同于一般的过程变量, 无法诊断故障与产品质量磁选管回收率之间的关系.针对单一层面数据建模与故障诊断方法的问题, 多位学者研究了多层面数据建模与故障诊断[7-16].比如, MacGraegor等[8]提出基于偏最小二乘(Partial least squares, PLS)的故障诊断方法, Zhou等[9]、Li等[10-11]提出全潜结构映射的建模与监控方法, Qin等[12]、Liu等[13]提出并发潜结构映射(Concurrent projection to latent structures, CPLS)多层面建模与监控方法, 实现过程与质量的联合监控, 已成为过程工业数据建模与故障诊断的新方向[17].
然而, CPLS方法将与质量直接相关的过程变化以及与质量潜在相关的过程变化, 划分成两个不同的子空间, 即共有子空间与过程变化残差空间, 分别用$T^2$统计指标与$Q$统计指标监控.对于竖炉焙烧过程而言, 新发故障可能同时影响这两个子空间.为此, 本文提出将两个统计指标合并为一个统一指标, 利用竖炉焙烧过程历史正常数据(以磁选管回收率作为质量变量, 以燃烧室温度、还原煤气流量、炉顶废气温度、炉内负压、燃烧煤气热值和加热空气流量作为过程变量), 建立竖炉焙烧过程并发潜结构映射模型及精简CPLS模型用于综合监控.在此基础上, 考虑到CPLS重构贡献难以满足正确诊断的必要条件, 本文提出CPLS精简重构贡献方法, 用于诊断竖炉焙烧过程质量相关故障变量.
1. 竖炉焙烧过程工艺及故障描述
1.1 竖炉焙烧过程工艺描述
首先, 矿石经皮带输送到炉顶的储槽, 之后进入炉膛内通过预热带和炉顶废气热交换进行预热, 使温度升至150$\sim $200${^\circ}$C; 预热后的矿石进入加热带继续加热至1 045$\sim $1 160${^\circ}$C左右; 其后矿石进入还原带与还原煤气接触, 使其温度降低到570${^\circ}$C左右.期间, 与还原煤气中的CO和H$_{2}$进行还原反应, 生成强磁性铁矿物Fe$_{3}$O$_{4}$.
适宜条件下, 还原带中的矿石发生如下反应:
$ \begin{array}{l} 3\text{Fe}_2 \text{O}_3 +\text{CO}\buildrel {570{^\circ}\rm C} \over \longrightarrow 2\text{Fe}_3 \text{O}_4 +\text{CO}_2 \\ 3\text{Fe}_2 \text{O}_3 +H_2 \buildrel {570{^\circ}\rm C} \over \longrightarrow 2\text{Fe}_3 \text{O}_4 +\text{H}_2 \text{O} \\ \end{array} $
其中, 生成物Fe$_{3}$O$_{4}$为强磁性铁矿, 磁选管回收率较高, 利于磁选出高品位精矿.
还原反应不充分时, 生成物为Fe$_{3}$O$_{4}$和Fe$_{2}$O$_{3}$的混合物.还原反应过度时, 还会发生下面的反应:
$ \begin{array}{l} \text{Fe}_3 \text{O}_4 +\text{CO}\buildrel {>570{^\circ}\rm C} \over \longrightarrow 3\text{FeO}+\text{CO}_2 \\ \text{Fe}_3 \text{O}_4 +\text{H}_2 \buildrel {>570{^\circ}\rm C} \over \longrightarrow 3\text{FeO}+\text{H}_2 \text{O} \\ \end{array} $
其中, 生成物FeO是弱磁性物质, 不利于磁选工序提取强磁性的高品位精矿.由此看出, 还原不充分与过还原都会降低产品质量, 以磁选管回收率来衡量.磁选管回收率降低会直接影响下游球磨工序和磁选工序的工艺指标:金属回收率和精矿品位.
图 1中, 质量变量为磁选管回收率$\gamma$; 相关的过程变量主要包括燃烧室温度、还原煤气流量、炉顶废气温度、炉内负压、加热煤气热值和加热空气流量.正常工况下, 上述过程变量都将维持在合适的范围内, 保证还原反应正常运行.因物料能量耦合关系, 以及底层回路控制与上层运行控制作用使得过程变量在正常工况下具有特定的相关关系潜结构, 可利用该相关关系建立正常工况数据并发潜结构模型.
1.2 竖炉焙烧过程故障描述
竖炉焙烧实际生产过程中, 由于原料矿石的大小、成分与种类频繁变化, 在过程操作与运行条件不相符时, 可发生与运行安全密切相关的冒火和上火故障, 也可发生欠还原与过还原这类与磁选管回收率相关的运行故障, 故障具体描述为
1) 冒火:炉内负压较低, 使得炉膛内加热带中的火焰不能及时上升, 从火眼喷出炉外;
2) 上火:未完全燃烧的煤气进入烟道燃烧, 使得炉内的火焰从炉顶喷出;
3) 欠还原:矿石在竖炉焙烧过程没有充分的还原就被搬出炉体;
4) 过还原:还原成的磁铁矿(主要成分Fe$_{3}$O$_{4}$)进一步被还原为弱磁铁矿后, 才搬出炉体.
竖炉焙烧过程工况异常, 直接影响生产与人员安全; 磁选管回收率降低, 直接影响下游工序的工艺指标(金属回收率和精矿品位).长期以来, 竖炉焙烧过程的操作者主要依据经验知识通过观察炉体表面的现象进行判断与决策以避免工艺设备可能会发生的故障, 但由于竖炉焙烧过程的过程变量较多, 仅仅依靠操作员观察炉体外观来进行故障诊断, 无法诊断故障原因, 也无法分析故障是否影响质量指标磁选管回收率.而且, 故障的现有理解将上火故障与过还原等区分为不同的故障, 而实际生产中, 上火故障等异常工况也会影响磁选管回收率.
异常运行工况可导致过程变量的异常变化、过程变量间的相关关系改变以及过程变量与质量变量间相关关系的改变.为此, 本文下一节提出利用历史数据建立磁选管回收率与过程变量间的潜结构相关关系; 在此基础上, 针对新发故障可能影响历史已发故障残差空间异常变化, 提出并发潜结构映射竖炉焙烧过程综合监控方法和基于CPLS精简重构贡献的故障变量诊断方法.
2. 基于并发潜结构映射的竖炉焙烧过程综合故障诊断
2.1 竖炉焙烧过程综合监控
2.1.1 竖炉焙烧过程并发潜结构建模
首先, 以燃烧室温度$T_{O1}$、还原煤气流量$F_h$、炉顶废气温度$T_{O2}$、炉内负压$P_f $、加热煤气热值$h$和加热空气流量$F_k$作为过程变量, 构造样本向量, $\boldsymbol{x}(k)=\left[{T_{O1} (k), F_h (k), T_{O2} (k), P_f (k), h(k), F_k (k)} \right]^{\rm T}$ ($k$为采样时刻), 过程变量数为$m=6$; 以磁选管回收率$\gamma $作为质量变量, 构造样本向量${\boldsymbol{y}}^\ast (k)=\left[{{\gamma }^\ast (k)} \right]$, 质量变量数为$p=1$.
正常工况下, 根据磁选管回收率同步采集过程变量, 构建过程变量建模数据阵${\pmb X}^{\ast }$和质量变量建模数据阵${\pmb Y}^{\ast }$, 如式(1)和(2)所示:
$ \begin{equation} \label{eq1} {{\pmb X}}^\ast =\left[\boldsymbol{x}^\ast (1) \quad {\boldsymbol{x}^\ast (2) \quad \cdots \quad \boldsymbol{x}^\ast (n)} \right]^{\rm T} \end{equation} $
(1) $ \begin{equation} {{\pmb Y}}^\ast =\left[{{\boldsymbol{y}}^\ast (1) \quad {\boldsymbol{y}}^\ast (2) \quad \cdots \quad {\boldsymbol{y}}^\ast (n)} \right]^{\rm T} \end{equation} $
(2) 式中, $n$为样本数;
接下来, 对{${\pmb X}^{\ast }$, ${\pmb Y}^{\ast }$}进行标准化处理, 即对每一列各元素减去该列的均值后除以该列的标准差, 以消除变量单位和测量范围的影响, 得到${\{}{{\pmb X}, {\pmb Y}}{\}}$.
其后, 利用过程数据阵${\pmb X}$与质量数据阵${\pmb Y}$, 采用并发潜结构映射算法建立过程变量与磁选管回收率间的潜结构模型[12],
$ \begin{equation} \label{eq6} {\pmb X}={\pmb T}_c {{\pmb R}}_c^{{\pmb † }} {{ +\pmb T}}_x { {\pmb P}}_x^{\rm T} {{ +\tilde {\pmb X}}} \end{equation} $
(3) $ \begin{equation} \label{eq7} {\pmb Y}={\pmb T}_c {{\pmb Q}}^{\rm T}_c {{ +\pmb T}}_y {{\pmb P}}_y^{\rm T} {{ +\tilde{\pmb Y}}} \end{equation} $
(4) 其中, ${{\pmb R}}_c \in {\bf R}^{6\times l_c }$、${{\pmb P}}_x \in {\bf R}^{6\times l_x }$、${{\pmb Q}}_c \in {\bf R}^{1\times l_c }$和${ {\pmb P}}_y \in {\bf R}^{1\times l_y }$为负荷阵, $l_c$、$l_x$和$l_y$分别为共有变化、过程特有变化和磁选管回收率特有变化的潜变量个数; ${{\pmb T}}_c $、${ {\pmb T}}_x$和${{\pmb T}}_y$为得分矩阵; ${{\pmb T}}_c $为可用于预测磁选管回收率的过程共有变化, ${{\pmb T}}_x $为对预测磁选回收率无用的过程特有变化, ${{\pmb T}}_y $表示不能由过程变量预测的磁选管回收率变化, ${{\tilde{\pmb X}}}$为过程变量残差, ${\tilde{\pmb {Y}}}$为磁选管回收率残差.
对单一样本, CPLS模型表示为
$ \begin{equation} \label{eq8} \boldsymbol{x}=\hat {\boldsymbol{x}}_c +{\tilde{\pmb {x}}}_c ={{\pmb R}}_c^{{{\pmb †}}{\rm T}} {{\pmb t}}_c {+{\pmb P}}_x {{\pmb t}}_x +\tilde {\boldsymbol{x}} \end{equation} $
(5) $ \begin{equation} \label{eq9} {\boldsymbol{y}}={\hat{\pmb {y}}}_c +{\tilde{\pmb {y}}}_c ={{\pmb Q}}_c {{\pmb t}}_c {+{\pmb P}}_y {{\pmb t}}_y +\tilde {\pmb{y}} \end{equation} $
(6) 其中
$ \begin{equation} \label{eq10} {{\tilde {\boldsymbol{x}}}}_c ={\boldsymbol{x}-{\pmb R}}_c^{{{\pmb †}}{\rm T}} {{\pmb t}}_c {={\pmb P}}_x {{\pmb t}}_x + \tilde {\boldsymbol{x}} \end{equation} $
(7) $ \begin{equation} \label{eq11} {\tilde {\boldsymbol{y}}}_c ={\pmb y-{\pmb Q}}_c {{\pmb t}}_c ={{\pmb P}}_y {{\pmb t}}_y +\tilde {\boldsymbol{y}} \end{equation} $
(8) $ \begin{equation} \label{eq12} {{\pmb t}}_c {={\pmb R}}_c^{\rm T} { \boldsymbol{x}} \end{equation} $
(9) $ \begin{equation} \label{eq13} {{\pmb t}}_x {={\pmb P}}_x^{\rm T} {\tilde {\boldsymbol{x}}}_c ={{\pmb P}}_x^{\rm T} \boldsymbol{x} \end{equation} $
(10) $ \begin{equation} \label{eq14} {{\pmb t}}_y {={\pmb P}}_y^{\rm T} {{\tilde{\pmb y}}}_c \end{equation} $
(11) $ \begin{equation} \label{eq15} {{\tilde {\boldsymbol{x}}}}=\left( {{ {\pmb I}-{\pmb P}}_x {{\pmb P}}_x^{\rm T}} \right){{\tilde { \boldsymbol{x}}}}_c =\left( {{{\pmb I}-{\pmb P}}_x {{\pmb P}}_x^{\rm T} } \right)\boldsymbol{x} \end{equation} $
(12) $ \begin{equation} \label{eq16} {\tilde {\boldsymbol{y}}}=\left( {{{\pmb I}-{\pmb P}}_y {{\pmb P}}_y^{\rm T}} \right){\tilde {\boldsymbol{y}}}_c \end{equation} $
(13) 需要说明的是, 过程变化${{\pmb T}}_c $与${{\tilde {\pmb X}}}$ (或${ {\pmb t}}_c$与${{\tilde{\boldsymbol{x}}}}$)均可能与磁选管回收率相关, 有必要将二者合并.从而, 同时监控历史已发及未发的与磁选管回收率相关的过程故障; ${ {\pmb T}}_y $与${{\tilde {\pmb Y}}}$ (或${{\pmb t}}_y $与${\tilde{\pmb {y}}}$)均为不能由过程变化预测的磁选管回收率变化, 有必要将二者合并.从而, 同时监控历史已发及未发的不可由过程变化预测的磁选管回收率异常.然而, ${ {\pmb T}}_c $ (或${{\pmb t}}_c $)与${{\pmb T}}_y $ (或${{\pmb t}}_y $)为得分变化, 适于采用$T^2$统计量监控; ${{\tilde {\pmb X}}}$ (或${ {\tilde {\boldsymbol{x}}}}$)与${\tilde{\pmb {Y}}}$ (或${\tilde{\pmb {y}}}$)为残差变化, 适于采用Q统计量监控, 不能将二者简单合并.因此, 这里并未给出显式的精简并发潜结构模型; 在第3.1.2节中, 通过相应的综合指标来建立精简并发潜结构空间; 在此基础上, 实现竖炉焙烧过程的综合监控.
2.1.2 基于精简CPLS的竖炉焙烧过程综合监控
首先, 采集过程变量实时数据$\boldsymbol{x}^*$, 采用BP神经网络[18], 计算质量变量磁选管回收率预测值$\boldsymbol{y}^*$, 并根据${\{}{\pmb X}^*, {\pmb Y}^*{\}}$的均值和标准差对$\{\boldsymbol{x}^*, \boldsymbol{y}^*\} $归一化处理, 得到$\boldsymbol{x}, \boldsymbol{y}$.
接下来, 利用第3.1.1节建立的CPLS模型, 对给定的新数据样本$\{\boldsymbol{x}, \boldsymbol{y}\}$, 按照式(8)$\sim $(16)计算得分与残差, 对得分和残差分别建立二次型统计指标, 由$T^2$统计量监控得分的变化, 由$Q$统计量监控残差的变化.假设数据服从多元正态分布, $T^2$统计量服从$F$分布, $Q$统计量服从$\chi^2$分布, 从而根据Qin等[12]方法列出5个统计指标及其控制限如表 1所示.
表 1 竖炉焙烧过程统计指标及其控制限[12]Table 1 The statistics and control limits for shaft furnace roasting processes[12]统计指标 控制限 ${array}{l} T_c^2 ={{\pmb t}}_c^{\rm T} {{\pmb \Lambda }}_c^{-1} {{\pmb t}}_c= \\ \quad \, \boldsymbol{x}^{\rm T}{{\pmb R}}_c {{\pmb \Lambda }}_c^{-1} {{\pmb R}}_c^{\rm T} \boldsymbol{x} \\ {array}$ $\tau _c^2 =\frac{l_c (n^2-1)}{n(n-l_c)}F_{l_c, n-l_c, \alpha } $ ${array}{l} T_x^2 ={{\pmb t}}_x^{\rm T} {{\pmb \Lambda }}_x^{-1} {{\pmb t}}_x =\\ \quad \, \boldsymbol{x}^{\rm T}{{\pmb P}}_x {{\pmb \Lambda }}_x^{-1} {{\pmb P}}_x^{\rm T} \boldsymbol{x} \\ {array}$ $\tau _x^2 =\frac{l_x (n^2-1)}{n(n-l_x)}F_{l_x, n-l_x, \alpha } $ ${array}{l} T_y^2 ={{\pmb t}}_y^{\rm T} {{\pmb \Lambda }}_y^{-1} {{\pmb t}}_y =\\ \quad \, {{\tilde {\pmb y}}}_c^{\rm T} {{\pmb P}}_y {{\pmb \Lambda }}_y^{-1} {{\pmb P}}_y^{\rm T} {{\tilde{ \boldsymbol{y}}}}_c \\ {array}$ $\tau _y^2 =\frac{l_y (n^2-1)}{n(n-l_y)}F_{l_y, n-l_y, \alpha} $ ${array}{l} Q_x =\left\| {{{\tilde{ \boldsymbol{x}}}}} \right\|^{{2}} =\\ \quad \, \boldsymbol{x}^{\rm T}\left({{{\pmb I-P}}_x {{\pmb P}}_x^{\rm T} } \right)\boldsymbol{x} \\ {array}$ $\delta _x^2 =g_x \cdot \chi _{h_x, \alpha }^2 $ ${array}{l} Q_y =\left\| {{{\tilde {\boldsymbol{y}}}}} \right\|^{{2}} =\\ \quad \, {{\tilde {\boldsymbol{y}}}}_c^{{\rm T}} \left({{{\pmb I-P}}_y {{\pmb P}}_y^{{\rm T}} } \right){{\tilde {\boldsymbol{y}}}}_c \\ {array}$ $\delta _y^2 =g_y \cdot \chi _{h_y, \alpha }^2 $ 其中
$ \begin{equation} \label{eq17} {\pmb\Lambda} _c =\frac{1}{n-1}{{\pmb T}}_c^{\rm T} {{\pmb T}}_c \end{equation} $
(14) $ \begin{equation} \label{eq18} {\pmb\Lambda} _x =\frac{1}{n-1}{{\pmb T}}_x^{\rm T} {{\pmb T}}_x \end{equation} $
(15) $ \begin{equation} \label{eq19} {\pmb\Lambda} _y =\frac{1}{n-1}{{\pmb T}}_y^{\rm T} {{\pmb T}}_y \end{equation} $
(16) $g_x $、$h_x$为服从$\chi^2$分布的$Q_x $统计指标的系数和自由度; $g_y $、$h_y $为对于服从$\chi^2$分布的$Q_y$的系数和自由度, 可按Qin等[12]方法求得.
对于竖炉焙烧过程而言, 新发故障可能引起历史已发故障残差空间异常变化.因此, 根据表 1中求得的5个统计指标及其控制限, 将其精简为三个统计指标, 得到综合监控竖炉焙烧过程磁选管相关故障和过程相关故障的方案如下.
1) $T_c^2 $超出控制限, 检测出可预测的磁选管回收率相关故障; $Q_x $超出控制限, 检测出潜在的磁选管回收率相关故障.二者均为可能与磁选管回收率有关的过程变化.为此, 参照Yue等[19]的统计指标合并方法, 将二者合并为一个综合指标$\phi$, 定义如下:
$ \begin{equation} \label{eq20} \phi =\frac{T_c^2 }{\tau _c^2 }+\frac{Q_x }{\delta _x^2 }={\boldsymbol{x}}^{\rm T}{{\pmb \Phi }}_x \boldsymbol{x} \end{equation} $
(17) 其中
$ \begin{equation}{{\pmb \Phi }}_x =\frac{{{\pmb R}}_c {{\pmb \Lambda }}_c^{-1} {{\pmb R}}_c^{\rm T} }{\tau _c^2 }+\frac{{{\pmb I}-{\pmb P}}_x {{\pmb P}}_x^{\rm T} }{\delta _x^2 }\end{equation} $
(18) $\phi $的控制限$\varsigma _1 ^2 $近似服从$\chi^2$分布[12,20], 按照如下方式计算:
$ \begin{equation} \varsigma _1 ^2 =g_{\phi } \cdot \chi _{h_{\phi }, \alpha }^2 \end{equation} $
(19) 其中
$ \begin{equation} g_{\phi } =\frac{{\rm tr}({{\pmb S}}_x {{\pmb \Phi }}_x )^2}{{\rm tr}({ {\pmb S}}_x {{\pmb \Phi }}_x )} \end{equation} $
(20) $ \begin{equation} h_{\phi } =\frac{\left[{{\rm tr}({{\pmb S}}_x {{\pmb \Phi }}_x )} \right]^2}{{\rm tr}({{\pmb S}}_x {{\pmb \Phi }}_x )^2} \end{equation} $
(21) 且${{\pmb S}}_x ={{\pmb X}}^{\rm T}{{\pmb X}}/{(n-1)}$为$\boldsymbol{x}$样本协方差, $1-\alpha $为置信水平($\alpha =0.05$).
如果$\phi>\varsigma _1^2 $, 检测出与磁选管回收率相关的异常过程变化.
2) 当质量变量数较少时, 通常不再适合对${{\tilde {\pmb Y}}}_c $做进一步主元分解.因此, 当无法进行主元分解时, 直接将${{ \tilde {\pmb Y}}}_c $作为${\tilde{\pmb {Y}}}$建立相应的$Q$统计指标来进行监控; 否则, 由于$T_y^2 $或$Q_y $超出控制限, 均可检测出不可预测的质量故障, 按照Qin等[12]方法将二者合并为一个综合指标$\kappa$, 其定义如下:
$ \begin{equation} \kappa =\frac{T_y^2 }{\tau _y^2 }+\frac{Q_y }{\delta _y^2 }= {\tilde {\boldsymbol{y}}}_c^{\rm T} {{\pmb \Phi }}_y {\tilde{\pmb {y}}}_c \end{equation} $
(22) 其中
$ \begin{equation} {{\pmb \Phi }}_y =\frac{{{\pmb P}}_y {{\pmb \Lambda }}_c^{-1} { {\pmb P}}_c^{\rm T} }{\tau _y^2 }+\frac{{{\pmb I}-{\pmb P}}_y {{\pmb P}}_y^{\rm T} }{\delta _y^2 } \end{equation} $
(23) $\kappa$的控制限$\varsigma_2^2 $近似服从$\chi^2$分布, 按照如下方式计算[12]:
$ \begin{equation} \varsigma _2^2 =g_{\kappa} \cdot \chi _{h_{\kappa}, \alpha }^2 \end{equation} $
(24) 其中
$ \begin{equation} g_{\kappa} =\frac{{\rm tr}({{\pmb S}}_y {{\pmb \Phi }}_y )^2}{{\rm tr}({ {\pmb S}}_y {{\pmb \Phi }}_y )} \end{equation} $
(25) $ \begin{equation} h_{\kappa } =\frac{\left[{{\rm tr}({{\pmb S}}_y {{\pmb \Phi }}_y )} \right]^2}{{\rm tr}({{\pmb S}}_y {{\pmb \Phi }}_y )^2} \end{equation} $
(26) 且${{\pmb S}}_y ={{\tilde {\pmb Y}}}_c^{\rm T} {\tilde{\pmb {Y}}}_c/(n-1) $为${{\tilde {\boldsymbol{y}}}}_c $样本协方差.
如果$\kappa>\varsigma_2 ^2$, 检测出不可由过程变量预测的磁选管回收率相关故障.需要说明的是, 计算$\kappa$统计指标用到${\tilde {\boldsymbol{y}}}_c$, 根据式(8)实时计算${\tilde{\pmb {y}}}_c$需要与过程变量同步快速率采样的质量数据.为此, 通过计算质量预测值, 就可以在质量化验数据得到前实时监控质量特有变化, 诊断质量特有故障.在获得新的质量数据测量值时, 进行质量预测模型有效性的校验与修正.
3) 如果$T_x^2$超出控制限, 检测出过程特有故障; 如果只关心磁选管回收率相关故障, 可以忽略该类过程异常.
综合监控$\phi$、$T_x^2$和$\kappa$三个统计指标, 相应的精简并发潜结构映射(Simplified concurrent projection to latent structures, S-CPLS)模型说明如下.
1) 对于过程变化, $\phi$和$T_x^2$根据${\boldsymbol{x}}$重写为二次型的一般形式:
$ \begin{equation} {\rm Index}(\boldsymbol{x})=\boldsymbol{x}^{\rm T}{{\pmb M}\boldsymbol{x}}={{\pmb z}}_x^{\rm T} { {\pmb z}}_x \end{equation} $
(27) 其中, ${{\pmb z}}_x ={{\pmb M}}^{\frac{1}{2}}{\boldsymbol{x}}$, 对于$\phi$而言, ${\pmb M}=\Phi_x$, 对于$T_x^2$而言, ${\pmb M}={{\pmb P}}_x {{\pmb \Lambda }}_x^{-1} {{\pmb P}}_x^{\rm T}$.
从而, 将原始过程变化$\boldsymbol{x}$空间, 投影到新的精简过程变化${{\pmb z}}_x $空间.
2) 对于磁选管回收变化, $\kappa$根据${{\tilde {\boldsymbol{y}}}}_c $重写为二次型的一般形式:
$ \begin{equation} {\rm Index}({\boldsymbol{y}})={{\tilde {\boldsymbol{y}}}}_c ^{\rm T}{{\pmb M}}_y { \tilde{\pmb {y}}}_c ={{\pmb z}}_y^{\rm T} {{\pmb z}}_y \end{equation} $
(28) 其中, ${{\pmb z}}_y ={{\pmb M}}_y^{\frac{1}{2}} {{\tilde {\boldsymbol{y}}}}_c $, 对于$\kappa$而言, ${\pmb M}_{y}=\Phi_y$.
从而, 将原始${{\tilde {\boldsymbol{y}}}}_c $数据空间, 投影到新的精简磁选管回收率变化${{\pmb z}}_y$空间.
2.2 CPLS精简重构贡献用于诊断竖炉焙烧过程故障变量
在竖炉焙烧过程综合监控的基础上, 进行竖炉焙烧过程故障变量的诊断.虽然重构贡献比贡献图方法诊断准确率高[21], 但竖炉焙烧过程CPLS模型的$T_c^2$指标重构不能满足重构的必要条件.为此, 基于上节建立的S-CPLS模型, 提出CPLS精简重构贡献方法用于诊断竖炉焙烧过程故障变量.
重构贡献方法以沿着变量方向的故障监控指标的重构量作为变量对故障的贡献.
为了定义磁选管回收率相关故障与磁选管回收率无关故障的重构贡献, 本文参照Alcala等[21]的重构贡献定义, 利用$\phi$和$T_x^2$根据$\boldsymbol{x}$重写为二次型的一般形式${\rm Index}({\boldsymbol{x}})=\boldsymbol{x}^{\rm T}{\pmb {Mx}}$.沿着第$i$个过程变量的方向${{\pmb \xi }}_i $ ($i=1, 2, \cdots, m$)定义二次型重构指标为
$ \begin{equation} {\rm Index}(\boldsymbol{x})=\boldsymbol{x}^{\rm T}{{\pmb M}\boldsymbol{x}}={{\pmb z}}_x^{\rm T} { {\pmb z}}_x \end{equation} $
(29) 其中, 故障强度$f$最优值可根据${\rm Index}(\boldsymbol{x}_i^r)$对$f$求偏导为0求得, 为$f=\left({\pmb \xi }_i^{\rm T} {\pmb \xi }_i \right)^{-1}\left({\pmb \xi }_i^{\rm T} {\pmb {Mx}} \right)$, ${ {\pmb \xi }}_i $为单位阵的第$i$列.
从而, $\phi$和$T_x^2 $两个统计指标的重构贡献为
$ \begin{equation} \label{eq21} \begin{array}{l} \mbox{RBC}_i^{{\rm Index}\left( \boldsymbol{x} \right)} ={\rm Index}( \boldsymbol{x})-{\rm Index}(\boldsymbol{x}_i^r ) =\\ \quad \quad \quad \quad \quad \;\dfrac{\left( {\pmb \xi }_i^{\rm T} {{\pmb M}\boldsymbol{x}} \right)^2 }{{\pmb \xi }_i^{\rm T} {{\pmb M}{\pmb \xi }_i}} \\ \end{array} \end{equation} $
(30) 在此基础上, 根据$\mbox{RBC}_i^{{\rm Index}\left(\boldsymbol{x} \right)} $的相对大小诊断故障变量, 变量重构贡献的值越大, 对故障贡献越大.
3. 实验研究
利用竖炉焙烧过程运行控制仿真实验平台[22], 进行竖炉焙烧过程正常工况与异常工况运行实验.该仿真平台采用的回路控制策略、运行控制策略以及运行条件, 与实际现场一致.采集正常工况数据和异常工况数据, 对本文所提方法进行验证.
首先, 在历史正常运行工况下, 采集10小时时间段内的质量变量磁选管回收率测量值及与其同步采样的过程变量, 建立竖炉焙烧过程数据潜结构模型, 具体是: ${l_c} = 1,{l_x} = 4,{l_y} = 1;$${{\pmb R}}_c=[-1.7165, $ $ 0.3585, $ $0.6597, $ $-0.4038, $ $-0.6156, $ $0.9101]^{\rm T}$; $\varsigma _\gamma ^2 =1.9209$; $\varsigma _\rho ^2 =3.1502$; $\tau _x^2=9.6768$;
$ {{\pmb P}}_x =\left[{{\begin{array}{*{20}c} {0.1644} &{-0.6627} &{-0.4172}&{-0.0391} \\ {-0.1327} &{-0.7198} &{0.3628}&{0.1084} \\ {-0.3552} &{-0.0374} &{-0.6565}&{-0.5028} \\ {0.6137} &{0.1479}&{-0.1696}&{-0.1046} \\ {0.0549} &{-0.0467}& {0.4809} &{-0.8474} \\ {0.6705} &{-0.1314} &{-0.0578} &{-0.0701} \\ \end{array} }}\right] $
接下来, 正常工况采集的400个样本如图 2(a)所示, 上火故障工况下采集的200个样本如图 2(b)所示, 以第36个样本时刻引入上火故障, 102个样本时刻恢复正常.对于正常工况数据和故障工况数据, 分别将本文所提S-CPLS综合监控方法与PCA、PLS、CPLS三种方法进行比较分析.
正常工况下, PCA方法、PLS方法、CPLS方法以及S-CPLS方法, 过程监控的实验结果分别如图 3(a)$\sim$(d)所示, 经统计, 误报率分别为4%、3.25%、2%、1.33%, S-CPLS方法误报率最低.
上火故障工况下, 相比较图 2(a)而言, 图 2(b)中各过程变量单一变量均看不出明显异常变化, 因而难以采用基于规则及人工方法诊断. PCA、PLS、CPLS及S-CPLS四种方法的故障监控实验结果见图 4(a)$\sim $(d). 图 4(a)中, PCA方法以$T^2$指标和$Q$指标超限监控过程运行异常, $T^2$指标和$Q$指标从36个样本时刻开始超限, 从而检测到异常工况, 但无法诊断其是否与磁选管回收率异常有关. 图 4(b)中, PLS方法以$T^2$指标超限监控与质量相关的过程运行异常, $T^2$指标从36个样本时刻开始超限, 从而诊断出磁选管回收率相关的过程故障, 但PLS的$T^2$指标在第120个样本时刻后超限, 与如图 2(b)所示的磁选管回收率恢复正常不符. 图 4(c)中, CPLS方法以$T_c^2$和$Q_x$指标超限监控与质量相关的过程运行异常、以$Q_y$指标超限监控质量特有异常, $T_c^2$指标36个样本时刻开始超限、102个样本开始恢复正常, 有效诊断出该时段内发生磁选管回收率相关的过程故障; 但$Q_x$指标在130个样本时刻后仍超限, 这与102个样本时刻磁选管回收率已恢复正常不一致, 为误报; 另外, 由图 4(c)的$T_c^2$和$Q_x$, CPLS方法过程变量与质量变量的共有子空间以及残差子空间都发生异常变化.为此, S-CPLS方法将二者合并为如图 4(d)所示的$\phi$综合指标. 图 4(d)中, S-CPLS方法以$\phi$指标超限监控质量相关异常、以$T_x^2$指标超限监控过程特有异常, $\phi$指标第36个样本时刻开始超限、第102个样本时刻恢复到正常, 有效诊断出磁选管回收率相关的过程故障, 并由120个样本时刻后$T_x^2$指标超限检测到过程特有变化, 经分析是由运行控制引起的与磁选管回收率无关的正常变化. S-CPLS方法的上述故障诊断结果与实际的仿真实验条件一致, 表明了S-CPLS综合监控的有效性.统计故障工况下, PCA、PLS、CPLS以及S-CPLS的误报率分别为49.63%、52.59%、17.78%、14.07%, 漏报率分别为1.54%、44.62%、0%、0%, 采用本文所提的S-CPLS综合监控方法的诊断准确率更高.
在综合监控的基础上, 进行故障变量的诊断.对于故障发生时刻(第36个样本), CPLS和S-CPLS的重构贡献图如图 5(a)$\sim $(b)所示.由图 5(a)可以看出, 因质量变量数少, 沿各个变量方向对$T_c^2 $指标重构贡献相等, 无法有效诊断磁选管回收率相关故障变量. 图 5(b)诊断出与磁选管回收率下降相关的故障变量为燃烧室温度.根据过程知识, 与上火故障发生机理, 以及伴随其发生过还原发生机理一致, 说明诊断出故障原因变量; 由图 5(b)的$T_x^2$指标重构贡献图诊断出炉内负压与加热空气流量偏离正常工况范围, 为过程特有变化故障变量.
4. 结语
本文提出了数据驱动的基于精简并发潜结构映射的竖炉焙烧过程多层面综合故障诊断方法, 并利用竖炉焙烧过程半实物仿真平台进行实验研究, 结果表明:
1) 正常工况下, 相比PCA、PLS和CPLS方法, 本文所提出的S-CPLS综合监控方法误报率更低;
2) 故障工况下, 本文所提出的S-CPLS的综合监控方法不仅检测出竖炉焙烧过程故障工况, 且诊断出上火故障会降低磁选管回收率, 以及与之相应的故障变量;
3) 本文基于S-CPLS重构贡献诊断出回路设定值之外的故障变量.
本文所提方法不仅在竖炉焙烧过程实验平台成功应用, 还可以推广应用于实际的竖炉焙烧过程.在方法软件化的基础上, 利用现场实际数据建立模型, 给出可视化的监控图与重构贡献图, 可以实现实际竖炉焙烧过程的综合故障诊断.此外, 本文所提出的精简重构贡献图等研究成果还可应用于高炉炼铁等其他复杂系统的多层面监控与诊断.
-
表 1 竖炉焙烧过程统计指标及其控制限[12]
Table 1 The statistics and control limits for shaft furnace roasting processes[12]
统计指标 控制限 $\begin{array}{l} T_c^2 ={{\pmb t}}_c^{\rm T} {{\pmb \Lambda }}_c^{-1} {{\pmb t}}_c= \\ \quad \, \boldsymbol{x}^{\rm T}{{\pmb R}}_c {{\pmb \Lambda }}_c^{-1} {{\pmb R}}_c^{\rm T} \boldsymbol{x} \\ \end{array}$ $\tau _c^2 =\frac{l_c (n^2-1)}{n(n-l_c)}F_{l_c, n-l_c, \alpha } $ $\begin{array}{l} T_x^2 ={{\pmb t}}_x^{\rm T} {{\pmb \Lambda }}_x^{-1} {{\pmb t}}_x =\\ \quad \, \boldsymbol{x}^{\rm T}{{\pmb P}}_x {{\pmb \Lambda }}_x^{-1} {{\pmb P}}_x^{\rm T} \boldsymbol{x} \\ \end{array}$ $\tau _x^2 =\frac{l_x (n^2-1)}{n(n-l_x)}F_{l_x, n-l_x, \alpha } $ $\begin{array}{l} T_y^2 ={{\pmb t}}_y^{\rm T} {{\pmb \Lambda }}_y^{-1} {{\pmb t}}_y =\\ \quad \, {{\tilde {\pmb y}}}_c^{\rm T} {{\pmb P}}_y {{\pmb \Lambda }}_y^{-1} {{\pmb P}}_y^{\rm T} {{\tilde{ \boldsymbol{y}}}}_c \\ \end{array}$ $\tau _y^2 =\frac{l_y (n^2-1)}{n(n-l_y)}F_{l_y, n-l_y, \alpha} $ $\begin{array}{l} Q_x =\left\| {{{\tilde{ \boldsymbol{x}}}}} \right\|^{{2}} =\\ \quad \, \boldsymbol{x}^{\rm T}\left({{{\pmb I-P}}_x {{\pmb P}}_x^{\rm T} } \right)\boldsymbol{x} \\ \end{array}$ $\delta _x^2 =g_x \cdot \chi _{h_x, \alpha }^2 $ $\begin{array}{l} Q_y =\left\| {{{\tilde {\boldsymbol{y}}}}} \right\|^{{2}} =\\ \quad \, {{\tilde {\boldsymbol{y}}}}_c^{{\rm T}} \left({{{\pmb I-P}}_y {{\pmb P}}_y^{{\rm T}} } \right){{\tilde {\boldsymbol{y}}}}_c \\ \end{array}$ $\delta _y^2 =g_y \cdot \chi _{h_y, \alpha }^2 $ -
[1] 柴天佑, 丁进良, 王宏, 苏春翌.复杂工业过程运行的混合智能优化控制方法.自动化学报, 2008, 34(5):505-515 http://www.aas.net.cn/CN/abstract/abstract13476.shtmlChai Tian-You, Ding Jin-Liang, Wang Hong, Su Chun-Yi. Hybrid intelligent optimal control method for operation of complex industrial processes. Acta Automatica Sinica, 2008, 34 (5):505-515 http://www.aas.net.cn/CN/abstract/abstract13476.shtml [2] Chai T Y, Ding J L, Wu F H. Hybrid intelligent control for optimal operation of shaft furnace roasting process. Control Engineering Practice, 2011, 19 (3):264-275 doi: 10.1016/j.conengprac.2010.05.002 [3] 吴峰华, 丁进良, 岳恒, 柴天佑.竖炉焙烧过程智能故障诊断系统.南京航空航天大学学报, 2006, 38 (S1):91-94 http://d.wanfangdata.com.cn/Periodical/njhkht2006z1023Wu Feng-Hua, Ding Jin-Liang, Yue Heng, Chai Tian-You. Intelligent fault diagnosis system for roasting process of shaft furnace. Journal of Nanjing University of Aeronautics & Astronautics, 2006, 38 (S1):91-94 http://d.wanfangdata.com.cn/Periodical/njhkht2006z1023 [4] 严爱军, 王普, 曾宇.竖炉燃烧过程智能故障预报系统.化工学报, 2008, 59 (7):1768-1772 http://www.oalib.com/paper/4208380Yan Ai-Jun, Wang Pu, Zeng Yu. Intelligent fault prediction system of combustion process in shaft furnace. Journal of Chemical Industry and Engineering (China), 2008, 59 (7):1768-1772 http://www.oalib.com/paper/4208380 [5] Lu X L, Liu Q, Chai T Y, Qin S J. Data-driven fault diagnosis of shaft furnace roasting processes using reconstruction and reconstruction-based contribution approaches. IFAC Proceedings Volumes, 2014, 47 (3):8897-8902 doi: 10.3182/20140824-6-ZA-1003.01188 [6] Qin S J. Survey on data-driven industrial process monitoring and diagnosis. Annual Reviews in Control, 2012, 36 (2):220-234 doi: 10.1016/j.arcontrol.2012.09.004 [7] 文成林, 吕菲亚, 包哲静, 刘妹琴.基于数据驱动的微小故障诊断方法综述.自动化学报, 2016, 42 (9):1285-1299 http://www.aas.net.cn/CN/abstract/abstract18918.shtmlWen Cheng-Lin, Lv Fei-Ya, Bao Zhe-Jing, Liu Mei-Qin. A review of data driven-based incipient fault diagnosis. Acta Automatica Sinica, 2016, 42 (9):1285-1299 http://www.aas.net.cn/CN/abstract/abstract18918.shtml [8] MacGregor J F, Kourti T. Statistical process control of multivariate processes. Control Engineering Practice, 1995, 3 (3):403-414 doi: 10.1016/0967-0661(95)00014-L [9] Zhou D H, Li G, Qin S J. Total projection to latent structures for process monitoring. AIChE Journal, 2010, 56 (1):168-178 doi: 10.1002/aic.11977/abstract [10] Li G, Qin S J, Ji Y D, Zhou D H. Total PLS based contribution plots for fault diagnosis. Acta Automatica Sinica, 2009, 35 (6):759-765 http://or.nsfc.gov.cn/bitstream/00001903-5/71521/1/1000004259143.pdf [11] Li G, Qin S J, Zhou D H. Geometric properties of partial least squares for process monitoring. Automatica, 2010, 46 (1):204-210 doi: 10.1016/j.automatica.2009.10.030 [12] Qin S J, Zheng Y Y. Quality-relevant and process-relevant fault monitoring with concurrent projection to latent structures. AIChE Journal, 2013, 59 (2):496-504 doi: 10.1002/aic.v59.2 [13] Liu Q, Qin S J, Chai T Y. Multiblock concurrent PLS for decentralized monitoring of continuous annealing processes. IEEE Transactions on Industrial Electronics, 2014, 61 (11):6429-6437 doi: 10.1109/TIE.2014.2303781 [14] 曹玉苹, 黄琳哲, 田学民.一种基于DIOCVA的过程监控方法.自动化学报, 2015, 41 (12):2072-2080 http://www.aas.net.cn/CN/abstract/abstract18780.shtmlCao Yu-Ping, Huang Lin-Zhe, Tian Xue-Min. A process monitoring method using dynamic input-output canonical variate analysis. Acta Automatica Sinica, 2015, 41 (12):2072-2080 http://www.aas.net.cn/CN/abstract/abstract18780.shtml [15] Zhu Q Q, Liu Q, Qin S J. Concurrent canonical correlation analysis modeling for quality-relevant monitoring. IFAC-PapersOnLine, 2016, 49 (7):1044-1049 doi: 10.1016/j.ifacol.2016.07.340 [16] Dong Y N, Qin S J. Dynamic-inner partial least squares for dynamic data modeling. IFAC-PapersOnLine, 2015, 48 (8):117-122 doi: 10.1016/j.ifacol.2015.08.167 [17] 刘强, 秦泗钊.过程工业大数据建模研究展望.自动化学报, 2016, 42(2):161-171 http://www.aas.net.cn/CN/abstract/abstract18807.shtmlLiu Qiang, Qin S J. Perspectives on big data modeling of process industries. Acta Automatica Sinica, 2016, 42 (2):161-171 http://www.aas.net.cn/CN/abstract/abstract18807.shtml [18] 严爱军, 柴天佑.磁选管回收率智能混合预报方法.信息与控制, 2005, 34 (6):759-764 http://d.wanfangdata.com.cn/Periodical/xxykz200506024Yan Ai-Jun, Chai Tian-You. Intelligent hybrid prediction method of magnetic tube recovery rate. Information and Control, 2005, 34 (6):759-764 http://d.wanfangdata.com.cn/Periodical/xxykz200506024 [19] Yue H H, Qin S J. Reconstruction-based fault identification using a combined index. Industrial & Engineering Chemistry Research, 2001, 40 (20):4403-4414 doi: 10.1021/ie000141%2B [20] Box G E P. Some theorems on quadratic forms applied in the study of analysis of variance problems, I. effect of inequality of variance in the one-way classification. The Annals of Mathematical Statistics, 1954, 25 (2):290-302 doi: 10.1214/aoms/1177728786 [21] Alcala C F, Qin S J. Reconstruction-based contribution for process monitoring. Automatica, 2009, 45 (7):1593-1600 doi: 10.1016/j.automatica.2009.02.027 [22] 周平, 代伟, 柴天佑.竖炉焙烧过程运行优化控制系统的开发及实验研究.控制理论与应用, 2012, 29 (12):1565-1572 http://www.doc88.com/p-9979653061816.htmlZhou Ping, Dai Wei, Chai Tian-You. Exploration of operational optimization control system for shaft furnace roasting process and its experiment study. Control Theory & Applications, 2012, 29 (12):1565-1572 http://www.doc88.com/p-9979653061816.html 期刊类型引用(1)
1. 王光,杨晶惠,徐一平,焦建芳. 基于核直接分解的过热蒸汽系统故障诊断方法. 控制理论与应用. 2023(05): 825-832 . 百度学术
其他类型引用(2)
-