发明名称 丙烯聚合生产数据检测及故障诊断系统及方法
摘要 一种丙烯聚合生产数据检测及故障诊断系统,包括与丙烯聚合生产过程连接的现场智能仪表、DCS系统以及上位机,所述的DCS系统由数据接口、控制站、数据库构成;智能仪表、DCS系统、上位机依次相连,所述的上位机包括标准化处理模块、主元分析模块、残差分析模块、信号采集模块、待诊断数据确定模块以及故障判断模块。以及提出了一种故障诊断方法。本发明提供一种能够实现在线测量熔融指数、能够得到良好的诊断效果、有效减小误报率的丙烯聚合生产数据检测及故障诊断系统及方法。
申请公布号 CN101038485A 申请公布日期 2007.09.19
申请号 CN200610155416.6 申请日期 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)计算均值:<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> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math> 2)计算方差:<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> 3)标准化:<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> 其中,TX为训练样本,N为训练样本数,TX为训练样本的均值;主元分析模块,用于进行主元分析提取主成分,依照设置的主元方差提取率,采用协方差奇异值分解的方法,采用以下步骤来实现:①计算X的协方差阵,记为∑X;②对∑X进行奇异值分解,得到特征根λ1,λ2,...,λp,其中λ1≥λ2≥...≥λp,对应的特征向量矩阵为U;③计算总方差和每个特征值对应的方差贡献率,按各个特征值的方差贡献率从大到小累加直到总的方差贡献率达到给定值,记选取个数为k;④选取特征向量矩阵U的前k列,作为变换矩阵T;⑤计算主元,通过式子F=T×X计算得到主元F;残差分析模块,用于计算残差分析控制限,设置检验水平为α,控制限按下式(4)计算:<math> <mrow> <msub> <mi>Q</mi> <mi>&alpha;</mi> </msub> <mo>=</mo> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> <msup> <mrow> <mo>[</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mi>&alpha;</mi> </msub> <msub> <mi>h</mi> <mn>0</mn> </msub> <msqrt> <mn>2</mn> <msub> <mi>&theta;</mi> <mn>2</mn> </msub> </msqrt> </mrow> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <msub> <mi>&theta;</mi> <mn>2</mn> </msub> <msub> <mi>h</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <msub> <mi>&theta;</mi> <mn>2</mn> </msub> </mfrac> <mo>]</mo> </mrow> <mfrac> <mn>1</mn> <msub> <mi>h</mi> <mn>0</mn> </msub> </mfrac> </msup> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math> 式中:<math> <mrow> <msub> <mi>&theta;</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msubsup> <mi>&lambda;</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1,2,3</mn> </mrow> </math> <math> <mrow> <msub> <mi>h</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> <msub> <mi>&theta;</mi> <mn>3</mn> </msub> </mrow> <mrow> <mn>3</mn> <msubsup> <mi>&theta;</mi> <mn>2</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow> </math> λi是X协方差矩阵的特征值,Cα是正态分布置信度为α的统计;信号采集模块,用于设定每次采样的时间间隙,采集现场智能仪表的信号;待诊断数据确定模块,用于将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待诊断数据VX;故障诊断模块,用于对待检测数据VX用训练时得到的TX和σx 2进行标准化处理,并将标准化处理后的数据作为主元分析模块的输入,用训练时得到的变换矩阵T对输入进行变换,变换后矩阵输入到残差分析模块,计算输入数据的Q统计量,若Q<Qα,判定样本点Q统计正常,过程对象正常;若Q>Qα,判定样本点Q统计异常,过程对象出现故障。
地址 310027浙江省杭州市西湖区浙大路38号