专利名称：基于小波分析的工业生产过程故障诊断装置的制作方法
技术领域：
本实用新型涉及工业过程故障诊断领域，特别地，涉及一种基于小波分析的 工业生产过程故障诊断装置。(二) 背景技术由于产品质量、经济效益、安全以及环保的要求，工业过程和相关的控制系 统变的非常复杂，为了保证工业系统的正常运作，故障的诊断和检测在工业过程 中扮演着非常重要的角色。近年来，统计分析应用于过程监控和故障诊断得到了 广泛的研究。利用工业实测数据，采用统计的方法进行故障诊断，避开了复杂的机理分析， 求解相对方便。由于工业过程从本质上来说是具有多尺度特性的，比如微观、宏 观聚合动力学的空间多尺度特性，过程变量总是在不同采样率或尺度下得到的， 相应的控制或操作也是在不同的时间和空间尺度下发生的，多尺度特性往往对工 业过程产生重要有时甚至是本质上的影响。所以只有全面考虑到过程的复共线性、 非线性特性以及多尺度特性才能对过程进行准确、可靠的监控。但是目前的故障 诊断方法往往只考虑了工业过程的复共线性和非线性特性，而没有考虑到过程的 多尺度特性，对于多尺度特性影响严重的复杂工业过程的故障诊断往往难以得到 较好的诊断效果。(三) 发明内容为了克服已有的故障诊断系统的没有考虑过程的多尺度特性、难以得到较好 的诊断效果的不足，本实用新型提供一种同时考虑工业过程数据的复共线性、非 线性特性和多尺度特性、能够得到良好的诊断效果的的基于小波分析的工业生产 过程故障诊断装置。本实用新型解决其技术问题所釆用的技术方案是-一种基于小波分析的工业生产过程故障诊断装置，包括与工业过程对象连接 的现场智能仪表、DCS系统以及上位机，所述的DCS系统由数据接口、控制站、 数据库构成；智能仪表、DCS系统、上位机依次相连，所述的上位机包括标准化处理模块，用于对数据进行标准化处理，各变量的均值为0,方差为1， 得到输入矩阵X，采用以下过程来完成一 l Wi)计算均值- tx-:2]7^，A !'=12)计算方差《=^"^(巧—3)标准化X=CT —^(6)(7)(8)其中，TX为训练样本，N为训练样本数；小波分解功能模块，用于采用Mallat塔式分解算法将原始信号分解为一系列 近似信息和细节信息，采用以下步骤来实现① 将原始信号空间Fo可以分解为一系列逼近空间Fj与细节空间『;，其 中J是最粗的尺度，也称为分解的尺度；② 计算逼近空间^。空间^由尺度函数(化"/)jez〉张成，采用下式来计算-^ 、(K0 = 2-J/2 (9)③ 计算细节空间^。细节空间^由小波函数^,必，"l，…"Ae Z}张成，采用下式来计算R ,2力/2 (TO) 其中/是尺度因子，A:是平移因子； 从而得到原始信息的分解信息，采用下式进行计算/(o=2XA"o+i2化,"o(11 )其中，第一项代表近似信息，第二项为细节信息，逼近因子^,*与细节因子4,*采 用Mallat算法计算；所述的近似信息A/W与细节信息A7(0 (/ = 1, 2， ...，■/)，定义如下 4/(0 = 5>^^(0 (12)，)=2化,) (13) 主元分析功能模块，用于进行主元分析提取主成分，采用协方差奇异值分解的方法，采用以下步骤来实现① 计算X的协方差阵，记为Z^;② 对I^进行奇异值分解，得到特征根、，&， ...， V其中？i^^…Ap，
对应的特征向量矩阵为U; ◎计算总方差和每个特征值对应的方差贡献率，按各个特征值的方差贡 献率从大到小累加直到总的方差贡献率达到给定值；④ 选取特征向量矩阵U的前k列，作为变换系数矩阵T;⑤ 计算主元F = TxX;小波重构功能模块，用于进行小波重构，根据小波理论，将各尺度下得到的 主元相加，即得到总的主元；支持向量机分类器功能模块，用于核函数采用径向基函数 《(X,x) = exp(-jc, II/o"2)，将训练过程化为如下二次规划求解问题 ) = -^Za"局叫，x》 (14)! =1 ^ ,',/=1得到分类函数，即如下函数的符号函数/(x) = |>,",I(;c,，;c) + Z> (15)定义当f(X)〉K)，数据样本处于正常状态；当f(X)O时，处于异常状态； 信号采集模块，用于设定每次采样的时间间隙，采集现场智能仪表的信号； 待诊断数据确定模块，用于将采集的数据传送到DCS实时数据库中，在每个 定时周期从DCS数据库的实时数据库中，得到最新的变量数据作为待诊断数据VX;故障诊断模块，用于对待检测数据VX用训练时得到的^和o"〗进行标准化 处理，并将标准化处理后的数据作为小波分解模块的输入，用训练时相同的参数 对输入数据进行小波分解，得到的系数作为中主元分析模块的输入；用训练时得 到的变换矩阵T对输入进行变换，变换后矩阵输入到小波重构模块；将对应的数 据相加即得到原待测数据的主成分，并将所得到的主成分输入到支持向量机分类 器模块；将输入代入训练得到的判别函数，计算判别函数值，判别过程的状态。作为优选的一种方案所述的上位机还包括判别模型更新模块，用于定期 将过程状态正常的点添加到训练集VX中，输出到标准化处理模块、小波分解模 块、主元分析功能模块、小波重构功能模块，并更新支持向量机分类器的分类模型。作为优选的另一种方案所述的上位机还包括结果显示模块，用于将故障 诊断结果传给DCS系统，并在DCS的控制站显示过程状态，并通过DCS系统和 现场总线将过程状态信息传递到现场操作站进行显示。本实用新型同时考虑工业过程数据的复共线性、非线性特性和多尺度特性， 将主元分析、支持向量机和小波分析相结合，对工业过程进行故障诊断。其中， 主元分析用于处理工业过程生产数据的复相关性，支持向量机分类器用于解决非 线性分类问题，小波分析用于获取过程在不同尺度下的信息。本实用新型的有益效果主要表现在同时兼顾工业过程数据中的复相关性、 多尺度特性和非线性特性等特性，充分利用了主元分析，多尺度系统理论和支持 向量机各自的特点，将主元分析的解相关性能力，小波分析对信息不同尺度下的 强分解与重构能力以及支持向量机的多变量非线性映射能力很好地结合了起来， 发挥了各自的优势，使得故障诊断更加可靠有效，能更好的指导生产，提高生产 效益。
图1是本实用新型所提出的故障诊断系统的硬件结构图。图2是本实用新型所提出的故障诊断系统功能模块图。 图3是小波分析的Mallat算法分解图。 图4是本实用新型上位机的原理框图。 具体实施方式
以下结合附图对本实用新型作进一步描述。参照图l、图2、图3以及图4， 一种基于小波分析的工业生产过程故障诊断 系统，包括与工业过程对象1连接的现场智能仪表2、 DCS系统以及上位机6， 所述的DCS系统由数据接口3、控制站4、数据库5构成；智能仪表2、 DCS系 统、上位机6通过现场总线依次相连，所述的上位机6包括标准化处理模块7，用于对数据进行标准化处理，各变量的均值为O，方差为 1，得到输入矩阵X，采用以下过程来完成一i w1)计算均值<formula>formula see original document page 7</formula>(6)2) 计算方差^^^"7^(巧-^)， (7)<formula>formula see original document page 8</formula>3) 标准化1=汉—汉， (8) 其中，TX为训练样本，N为训练样本数；小波分解功能模块8，用于采用Mallat塔式分解算法将原始信号分解为一系 列近似信息和细节信息，采用以下步骤来实现① 将原始信号空间Fo可以分解为一系列逼近空间K/与细节空间％，其 中J是最粗的尺度，也称为分解的尺度；② 计算逼近空间K/。空间K/由尺度函数(^^)，AE々张成，采用下式 来计算K ，"K,) = 2-J/2 (9)③ 计算细节空间％■。细节空间^由小波函数{^,必，_/= 1， A e Z} 张成，采用下式来计算^={ 2力/2 (10) 其中/是尺度因子，6是平移因子； 从而得到原始信息的分解信息，采用下式进行计算AeZ 乂=1 AeZ其中，第一项代表近似信息，第二项为细节信息，逼近因子"^与细节因子《4采用Mallat算法计算；所述的近似信息A/(/)与细节信息"y/W (/ = 1， 2, ...， J)，定义如下4/")=2Xw) (12) A/w人2X,w(0 (13)主元分析功能模块9:"用于进行主元分析提取主成分，采用协方差奇异值分解的方法，采用以下步骤来实现① 计算X的协方差阵，记为Z^;② 对L"进行奇异值分解，得到特征根、，、，...，V其中、^21..Ap，对应的特征向量矩阵为U;③ 计算总方差和每个特征值对应的方差贡献率，按各个特征值的方差贡 献率从大到小累加直到总的方差贡献率达到给定值；
选取特征向量矩阵U的前k列，作为变换系数矩阵T; ⑤计算主元F = TxX; 小波重构功能模块IO，用于进行小波重构，根据小波理论，将各尺度下得到 的主元相加，即得到总的主元；支持向量机分类器功能模块11，用于核函数采用径向基函数 尺(X,，X) = exp(- IIJC-X, II /O"2)，将训练过程化为如下二次规划求解问题=— ZQr,. 少,力《(x,，x》 (14)' =1 ^ ',7=1得到分类函数，即如下函数的符号函数/C0 = ￡x",〖(w) + 6 (15)定义当《">=0，数据样本处于正常状态；当f(x)O时，处于异常状态； 信号采集模块12，用于设定每次采样的时间间隙，采集现场智能仪表的信号； 待诊断数据确定模块13，用于将采集的数据传送到DCS实时数据库中，在 每个定时周期从DCS数据库的实时数据库中，得到最新的变量数据作为待诊断数据VX;故障诊断模块14，用于对待检测数据VX用训练时得到的^和c^进行标准 化处理，并将标准化处理后的数据作为小波分解模块的输入，用训练时相同的参 数对输入数据进行小波分解，得到的系数作为中主元分析模块的输入；用训练时 得到的变换矩阵T对输入进行变换，变换后矩阵输入到小波重构模块；将对应的 数据相加即得到原待测数据的主成分，并将所得到的主成分输入到支持向量机分 类器模块；将输入代入训练得到的判别函数，计算判别函数值，判别过程的状态。所述的上位机还包括判别模型更新模块15,用于定期将过程状态正常的点 添加到训练集VX中，输出到标准化处理模块7、小波分解模块8、主元分析功能 模块9、小波重构功能模块IO，并更新支持向量机分类器模块ll的分类模型。所述的上位机还包括结果显示模块16，用于将故障诊断结果传给DCS系 统，并在DCS的控制站显示过程状态，并通过DCS系统和现场总线将过程状态 信息传递到现场操作站进行显示。
本实施例的工业过程故障诊断系统的硬件结构图如附图1所示，所述的故障 诊断系统核心由包括标准化模块7、小波分解模块8、主元分析模块9、小波重构 模块10、支持向量机分类器模块11等五大功能模块和人机界面的上位机6构成， 此外还包括现场智能仪表2， DCS系统和现场总线。所述的DCS系统由数据接 口3、控制站4、数据库5构成；工业过程对象l、智能仪表2、 DCS系统、上位 机6通过现场总线依次相连，实现信息流的上传和下达。故障诊断系统在上位机 6上运行，可以方便地与底层系统进行信息交换，及时应对系统故障。本实施例的故障诊断系统的功能模块图如附图2所示，主要包括标准化处理 模块7、小波分解模块8、主元分析模块9、小波重构模块IO、支持向量机分类器 模块ll等五大功能模块。本实用新型所述的故障诊断方法按照如下步骤进行实施1、 确定故障诊断所用的关键变量，从DCS数据库5的历史数据库中分别采集 系统正常和故障时这些变量的数据作为训练样本TX;2、 在上位机6的小波分解模块8、主元分析模块9和支持向量机分类器模块 11中，分别设置小波分解层数、主元分析方差提取率、支持向量机核参数 和置信概率等参数，设定DCS中的采样周期；3、 训练样本TX在上位机6中，依次经过标准化处理7、小波分解8、主元分 析9、小波重构10、支持向量机11等功能模块，采用以下步骤来完成诊断 系统的训练1)在上位机6的标准化处理功能模块7中，对数据进行标准化处理，使得各 变量的均值为0，方差为l，得到输入矩阵X。采用以下过程来完成-①计算均值<formula>formula see original document page 10</formula> (6)②计算方差<formula>formula see original document page 10</formula> (7)③标准化<formula>formula see original document page 10</formula>， 其中N为训练样本数。
上位机6的标准化处理功能模块7所进行的标准化处理能消除各变量因为量纲不同造成的影响。2) 在上位机6的小波分解功能模块8中，采用Mallat塔式分解算法将原始信号分解为一系列近似信息和细节信息。所述的上位机6中小波分解模块8的小 波分析采用db3小波，分解层数取为3 — 7。采用以下步骤来实现① 将原始信号空间Fo可以分解为一系列逼近空间J^与细节空间巧，其中J 是最粗的尺度，也称为分解的尺度；② 计算逼近空间F>空间^由尺度函数{似^)， A E Z》张成，采用下式来计算K ={ (,)| (0 = 2—"W2力-(9)③ 计算细节空间巧。细节空间巧由小波函数{牧必，7'= 1， ..." A: E Z〉张成，采用下式来计算-K = I W) = 2力7 V(2-々-(10)其中j'是尺度因子，A:是平移因子。
从而得到原始信息的分解信息，采用下式进行计算/") = Z "w (o +1 ￡ "" (o ( 11 )其中第一项代表近似信!l，第二项为fl等z信息。逼近因子aw与细节因子4^采用 Mallat算法计算。其中所采用的Mallat算法的塔式分解图，如图2所示。所述的近似信息^//(/)与细节信息A/(0 (/ = 1， 2，J)，定义如下 4州=1>為(0 (12)= IX,^,W) (13) 实际工业过程从本质12来说是具有多尺度特性的，各个尺度下的信息含量和所体现的系统特性是不同的，如果直接对系统进行建模，将忽略这种不同，从而导致会导致结果上的偏差。用小波分解提取各个尺度的信息，能更充分的挖掘过程数据所携带的信息，提升结果的精确度。3) 在上位机6的主元分析功能模块9中，进行主元分析提取主成分。所述的 主元分析总方差提取率大于80%，计算过程采用协方差奇异值分解的方法，采用 以下步骤来实现① 计算X的协方差阵，记为Z";② 对S"进行奇异值分解，得到特征根、，X2, ...， V其中、^^2^.2Xp， 对应的特征向量矩阵为U; ③计算总方差和每个特征值对应的方差贡献率，按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值； 选取特征向量矩阵U的前k列，作为变换系数矩阵T;⑤计算主元F = TxX。很显然，分析系统在一个低维空间要比在一个高维空间容易得多。主元分析 在力求数据信息丢失最少的原则下，对高维的变量空间降维，以得到变量体系的 少数几个线性组合，并且这几个线性组合所构成的综合变量将尽可能多地保留原 变量变异方面的信息。4) 在上位机6的小波重构功能模块10中，进行小波重构。 根据小波理论，将各尺度下得到的主元相加，即得到总的主元。5) 训练上位机6中的支持向量机分类器功能模块11的分类模型。 所述的上位机6中的支持向量机分类器功能模块9的核函数，采用径向基函数《0c,，x)-exp(-x, ||/o"2),将训练过程化为如下二次规划求解问题<formula>formula see original document page 12</formula> (14)从而得到分类函数，即如下函数的符号函数/(义)=Z 乂 , x) + ^ (15)定义当！^)>=0，数据样本处于正常状态；当f(x)O时，处于异常状态。 支持向量机基于统计学习理论，采用结构风险最小化准则，很好地解决了小样本、局部极小点、高维数等难题，用于分类问题能提高分类精度。4、 系统开始投运1) 用定时器，设置好每次采样的时间间隔；2) 现场智能仪表2检测过程数据并传送到DCS数据库5的实时数据库中；3) 上位机6在每个定时周期从DCS数据库5的实时数据库中，得到最新 的变量数据，作为待诊断数据VX;4) 上位机6的标准化处理功能模块7，用训练时得到的玩和^对待检测 数据VX进行标准化处理，并将标准化处理后的数据作为小波分解模 块8的输入；5) 上位机6的小波分解模块8，用训练时同样的参数对输入数据进行小波 分解，得到的系数作为上位机6中主元分析模块9的输入；12
6) 上位机6的主元分析模块9，用训练时得到的变换矩阵T对输入进行 变换，变换后矩阵输入到上位机6的小波重构模块10;7) 上位机6的小波重构模块10，将对应的数据相加即得到原待测数据的 主成分，并将所得到的主成分输入到上位机6的支持向量机分类器模 块ll;8) 上位机6的支持向量机分类器模块11，将输入代入训练得到的判别函 数，计算判别函数值，判别过程的状态，并在上位机6的人机界面上 显示过程的状态9) 上位机6将故障诊断结果传给DCS，并在DCS的控制站4显示过程状 态，同时通过DCS系统和现场总线将过程状态信息传递到现场操作站 进行显示，使得现场操作工可以及时应对。5、 分类器模型更新在系统投运过程中，定期将过程状态正常的点添加到训练集TX中，重复步 骤3的训练过程，以便及时更新上位机6的支持向量机分类器11的分类模型，使 分类器模型具有较好的分类效果。
权利要求1、一种基于小波分析的工业生产过程故障诊断装置，包括与工业过程对象连接的现场智能仪表、DCS系统以及上位机，所述的DCS系统由数据接口、控制站、数据库构成；智能仪表、DCS系统、上位机依次相连，其特征在于所述的上位机包括用于对数据进行标准化处理的标准化处理模块；用于采用Mallat塔式分解算法将原始信号分解为一系列近似信息和细节信息的小波分解功能模块；用于进行主元分析提取主成分，采用协方差奇异值分解的方法的主元分析功能模块；用于进行小波重构，根据小波理论，将各尺度下得到的主元相加，即得到总的主元的小波重构功能模块；用于核函数采用径向基函数K(xi，x)＝exp(-‖x-xi‖/σ2)，将训练过程化为二次规划求解问题的支持向量机分类器功能模块；用于设定每次采样的时间间隙，采集现场智能仪表的信号的信号采集模块；用于将采集的数据传送到DCS实时数据库中，在每个定时周期从DCS数据库的实时数据库中，得到最新的变量数据作为待诊断数据VX的待诊断数据确定模块；用于对待检测数据VX用训练时得到的和σx2进行标准化处理，并将标准化处理后的数据作为小波分解模块的输入，用训练时相同的参数对输入数据进行小波分解，得到的系数作为中主元分析模块的输入；用训练时得到的变换矩阵T对输入进行变换，变换后矩阵输入到小波重构模块；将对应的数据相加即得到原待测数据的主成分，并将所得到的组成分输入到支持向量机分类器模块；将输入代入训练得到的判别函数，计算判别函数值，判别过程的状态的故障诊断模块，所述现场智能仪表与信号采集模块数据连接，所述信号采集模块连接待诊断数据确定模块，所述的待诊断数据确定模块连接故障诊断模块，所述标准化处理模块与小波分解功能模块连接，所述小波分解功能模块与主元分析功能模块连接，所述主元分析功能模块与小波重构模块连接，所述小波重构模块与支持向量机分类器功能模块连接，所述支持向量机分类器功能模块与故障诊断模块连接。
2、 如权利要求1所述的基于小波分析的工业生产过程故障诊断装置，其特征在于 所述的上位机还包括用于定期将过程状态正常的点添加到训练集VX中，输出到标准化处理模块、 小波分解模块、主元分析功能模块、小波重构功能模块，并更新支持向量机分类 器的分类模型的判别模型更新模块，所述判别模型更新模块连接支持向量机分类 器功能模块。
3、 如权利要求1或2所述的的基于小波分析的工业生产过程故障诊断装置，其特 征在于所述的上位机还包括用于将故障诊断结果传给DCS系统，并在DCS的控制站显示过程状态，并通过DCS系统和现场总线将过程状态信息传递到现场操作站进行显示的结果显示模块，所述结果显示模块连接故障诊断模块。
专利摘要一种基于小波分析的工业生产过程故障诊断装置，包括与工业过程对象连接的现场智能仪表、DCS系统以及上位机，所述的DCS系统由数据接口、控制站、数据库构成；智能仪表、DCS系统、上位机依次相连，所述的上位机包括标准化处理模块、小波分解模块、主元分析功能模块、小波重构功能模块、支持向量机分类器功能模块以及故障判断模块。本实用新型提供一种同时考虑工业过程数据的复共线性、非线性特性和多尺度特性、能够得到良好的诊断效果的基于小波分析的工业生产过程故障诊断装置。
文档编号G05B19/418GK201017233SQ200620140289
公开日2008年2月6日 申请日期2006年11月23日 优先权日2006年11月23日
发明者刘兴高 申请人:浙江大学