发明名称 一种丙烯聚合生产中熔融指数检测的故障诊断系统及方法
摘要 一种丙烯聚合生产中熔融指数检测的故障诊断系统,包括与丙烯聚合生产过程连接的现场智能仪表、DCS系统以及上位机,所述的DCS系统由数据接口、控制站、数据库构成;智能仪表、DCS系统、上位机依次相连,所述的上位机包括标准化处理模块、独立成分分析模块、支持向量机分类器功能模块、信号采集模块、待诊断数据确定模块以及故障判断模块。以及提出了一种故障诊断方法。本发明提供一种求解相对方便、能够得到良好的诊断效果、有效减小误报率的丙烯聚合生产中熔融指数检测的故障诊断系统及方法。
申请公布号 CN101021723A 申请公布日期 2007.08.22
申请号 CN200610155415.1 申请日期 2006.12.22
申请人 浙江大学 发明人 刘兴高
分类号 G05B19/048(2006.01);G05B19/418(2006.01);G05B13/02(2006.01);G06F17/00(2006.01);G01N25/04(2006.01) 主分类号 G05B19/048(2006.01)
代理机构 杭州天正专利事务所有限公司 代理人 王兵;黄美娟
主权项 1、一种丙烯聚合生产中熔融指数检测的故障诊断系统,包括与丙烯聚合生产过程连接的现场智能仪表、DCS系统以及上位机,所述的DCS系统由数据接口、控制站、数据库构成;智能仪表、DCS系统、上位机依次相连,其特征在于:所述的上位机包括:标准化处理模块,用于对数据库中采集系统正常时关键变量的数据进行标准化处理,各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:1)计算均值:<![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>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>2)计算方差:<![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>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>3)标准化:<![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>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,TX为训练样本,N为训练样本数,<img file="A2006101554150002C4.GIF" wi="63" he="48" />为训练样本的均值;独立成分分析模块,用于依照独立成分个数,采用基于定点迭代的快速ICA算法计算解混矩阵W,具体的步骤如下:①随机选取范数为1的初始权重向量w<sub>i</sub>,若i≥2,则<![CDATA[ <math><mrow><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><msub><mi>w</mi><mi>i</mi></msub><mo>-</mo><msub><mi>W</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msubsup><mi>W</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup><msub><mi>w</mi><mi>i</mi></msub><mo>,</mo></mrow></math>]]></maths>其中W<sub>i-1</sub>=[w<sub>1</sub>w<sub>2</sub>L w<sub>i-1</sub>],i=1,...,m;②对w<sub>i</sub>进行迭代更新:<![CDATA[ <math><mrow><msubsup><mi>w</mi><mi>i</mi><mo>+</mo></msubsup><mo>=</mo><mi>E</mi><mo>{</mo><mi>xg</mi><mrow><mo>(</mo><msubsup><mi>w</mi><mi>i</mi><mi>T</mi></msubsup><mi>x</mi><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mi>E</mi><mo>{</mo><msup><mi>g</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msubsup><mi>w</mi><mi>i</mi><mi>T</mi></msubsup><mi>x</mi><mo>)</mo></mrow><mo>}</mo><msub><mi>w</mi><mi>i</mi></msub><mo>,</mo></mrow></math>]]></maths>其中w<sub>i</sub><sup>+</sup>表示更新后的权重向量,E是数学期望,g代表形式为g(x)=xexp(-x<sup>2</sup>/2)的函数,g’为g的倒数;③标准化处理<![CDATA[ <math><mrow><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>w</mi><mi>i</mi><mo>+</mo></msubsup><mo>/</mo><mo>|</mo><mo>|</mo><msubsup><mi>w</mi><mi>i</mi><mo>+</mo></msubsup><mo>|</mo><mo>|</mo><mo>,</mo></mrow></math>]]></maths>其中‖w<sub>i</sub><sup>+</sup>‖表示w<sub>i</sub><sup>+</sup>的范数;④若不收敛,返回②,否则一直迭代至i=m;⑤当更新的w<sub>i</sub>与原w<sub>i</sub>点积为1时,判别为收敛;⑥计算独立成分:S=WX;其中,S是独立成分矩阵,W是解混矩阵,X是输入矩阵;支持向量机分类器功能模块,用于依照支持向量机核参数和置信概率,计算核函数,采用径向基函数K(x<sub>i</sub>,x)=exp(-‖x-x<sub>i</sub>‖/σ<sup>2</sup>),将训练过程化为二次规划求解问题:<![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>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>从而得到分类函数,即如下函数的符号函数f(x):<![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>5</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中,α<sub>i</sub>(i=1,…,N)是拉各朗日乘子,x<sub>i</sub>(i=1,…,N)是输入向量,y是输出变量,ω是支持向量机超平面的法向量,决定超平面的的方向,b为决定超平面位置的参数,δ为核参数;信号采集模块,用于设定每次采样的时间间隙,采集现场智能仪表的信号;待诊断数据确定模块,用于将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;故障诊断模块,用于对待检测数据VX用训练时得到的<img file="A2006101554150003C3.GIF" wi="63" he="47" />和σ<sub>x</sub><sup>2</sup>进行标准化处理,并将标准化处理后的数据作为独立成分分析模块的输入,用训练时得到的解混矩阵W对输入进行变换,变换后矩阵输入到支持向量机分类器功能模块,将输入数据代入训练得到的判别函数f(x),计算判别函数值,当f(x)>=0,数据样本处于正常状态;当f(x)<0时,处于异常状态。
地址 310027浙江省杭州市西湖区浙大路38号