发明名称 一种基于小波分析的工业生产过程故障诊断系统及方法
摘要 一种基于小波分析的工业生产过程故障诊断系统,包括与工业过程对象连接的现场智能仪表、DCS系统以及上位机,所述的DCS系统由数据接口、控制站、数据库构成;智能仪表、DCS系统、上位机依次相连,所述的上位机包括标准化处理模块、小波分解模块、主元分析功能模块、小波重构功能模块、支持向量机分类器功能模块以及故障判断模块。以及提出了一种故障诊断方法。本发明提供一种同时考虑工业过程数据的复共线性、非线性特性和多尺度特性、能够得到良好的诊断效果的基于小波分析的工业生产过程故障诊断系统及方法。
申请公布号 CN100480926C 申请公布日期 2009.04.22
申请号 CN200610154825.4 申请日期 2006.11.23
申请人 浙江大学 发明人 刘兴高;阎正兵
分类号 G05B19/418(2006.01)I 主分类号 G05B19/418(2006.01)I
代理机构 杭州天正专利事务所有限公司 代理人 王 兵;黄美娟
主权项 1、一种基于小波分析的工业生产过程故障诊断系统,包括与工业过程对象连接的现场智能仪表、DCS系统以及上位机,所述的DCS系统由数据接口、控制站、数据库构成;智能仪表、DCS系统、上位机依次相连,其特征在于:所述的上位机包括:标准化处理模块,用于对数据进行标准化处理,各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:1)计算均值:<maths num="0001"><![CDATA[<math><mrow><mover><mi>TX</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>T</mi><msub><mi>X</mi><mi>i</mi></msub><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>2)计算方差:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>TX</mi><mi>i</mi></msub><mo>-</mo><mover><mi>TX</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>3)标准化:<maths num="0003"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mi>TX</mi><mo>-</mo><mover><mi>TX</mi><mo>&OverBar;</mo></mover></mrow><msub><mi>&sigma;</mi><mi>x</mi></msub></mfrac><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,TX为训练样本,N为训练样本数;小波分解功能模块,用于采用Mallat塔式分解算法将原始信号分解为一系列近似信息和细节信息,采用以下步骤来实现:①将原始信号空间V<sub>0</sub>可以分解为一系列逼近空间V<sub>J</sub>与细节空间W<sub>j</sub>,其中J是最粗的尺度,也称为分解的尺度;②计算逼近空间V<sub>J</sub>。空间V<sub>J</sub>由尺度函数<img file="C200610154825C0002190746QIETU.GIF" wi="307" he="66" />张成,采用下式来计算:<img file="C200610154825C00024.GIF" wi="1055" he="64" />③计算细节空间W<sub>j</sub>。细节空间W<sub>j</sub>由小波函数{ψ<sub>j,k</sub>(t),j=1,...,J,k∈Z}张成,采用下式来计算:W<sub>j</sub>={ψ<sub>j,k</sub>(t)|ψ<sub>j,k</sub>(t)=2<sup>-j/2</sup>ψ(2<sup>-j</sup>t-k)}    (10)其中j是尺度因子,k是平移因子;④从而得到原始信息的分解信息,采用下式进行计算:<img file="C200610154825C00025.GIF" wi="1081" he="125" />其中,第一项代表近似信息,第二项为细节信息,逼近因子a<sub>J,k</sub>与细节因子d<sub>j,k</sub>采用Mallat算法计算;所述的近似信息A<sub>J</sub>f(t)与细节信息D<sub>j</sub>f(t)(j=1,2,...,J),定义如下:<img file="C200610154825C00031.GIF" wi="993" he="94" /><maths num="0004"><![CDATA[<math><mrow><msub><mi>D</mi><mi>j</mi></msub><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Element;</mo><mi>Z</mi></mrow></munder><msub><mi>d</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>&psi;</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>主元分析功能模块,用于进行主元分析提取主成分,采用协方差奇异值分解的方法,采用以下步骤来实现:①计算X的协方差阵,记为∑x;②对∑x进行奇异值分解,得到特征根λ<sub>1</sub>,λ<sub>2</sub>,...,λ<sub>p</sub>,其中λ<sub>1</sub>≥λ<sub>2</sub>≥...≥λ<sub>p</sub>,对应的特征向量矩阵为U;③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值;④选取特征向量矩阵U的前k列,作为变换系数矩阵T;⑤计算主元:F=T×X;小波重构功能模块,用于进行小波重构,根据小波理论,将各尺度下得到的主元相加,即得到总的主元;支持向量机分类器功能模块,用于核函数采用径向基函数K(x<sub>i</sub>,x)=exp(-‖x-x<sub>i</sub>‖/σ<sup>2</sup>),将训练过程化为如下二次规划求解问题:<maths num="0005"><![CDATA[<math><mrow><mi>&omega;</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><msub><mi>&alpha;</mi><mi>j</mi></msub><msub><mi>y</mi><mi>i</mi></msub><msub><mi>y</mi><mi>j</mi></msub><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths>得到分类函数,即如下函数的符号函数:<maths num="0006"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>y</mi><mi>i</mi></msub><msub><mi>&alpha;</mi><mi>i</mi></msub><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mi>b</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow></math>]]></maths>定义当f(x)&gt;=0,数据样本处于正常状态;当f(x)&lt;0时,处于异常状态;信号采集模块,用于设定每次采样的时间间隙,采集现场智能仪表的信号;待诊断数据确定模块,用于将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;故障诊断模块,用于对待诊断数据VX用训练时得到的TX和<img file="C200610154825C0003190941QIETU.GIF" wi="59" he="73" />进行标准化处理,并将标准化处理后的数据作为小波分解模块的输入,用训练时相同的参数对输入数据进行小波分解,得到的系数作为中主元分析模块的输入;用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到小波重构模块;将对应的数据相加即得到原待测数据的主成分,并将所得到的主成分输入到支持向量机分类器模块;将输入代入训练得到的判别函数,计算判别函数值,判别过程的状态。
地址 310027浙江省杭州市西湖区浙大路38号