发明名称 一种低误报率的改进主元分析过程监测方法及其监测系统
摘要 一种低误报率的改进主元分析过程监测方法及其监测系统,属于工业过程监视方法技术领域。收集能体现工业过程正常运行时的工作状态的数据;对收集的过程数据进行标准化处理,建立改进主元分析模型;计算统计指标PVR和CVR的低误报率控制限;过程监测系统在线监测过程时,计算改进主元分析的统计指标PVR和CVR,与计算得到的低误报率控制限进行比对;当统计指标超过对应的控制限时,系统报警。本发明还提供一种使用如上述的低误报率的改进主元分析过程监测方法进行过程监测的过程监测系统。与现有技术相比,本发明能有效地降低基于改进主元分析方法的监测系统的误报率,并能保持改进主元分析方法对微小故障的检测能力;提高了监测系统的可用性。
申请公布号 CN103207567A 申请公布日期 2013.07.17
申请号 CN201310075192.8 申请日期 2013.03.08
申请人 华北电力大学 发明人 邱天;杨婷婷;牛玉广;曾德良;刘吉臻
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 薄观玖
主权项 1.一种低误报率的改进主元分析过程监测方法,在过程监测系统中使用,其特征在于,该过程监测方法步骤如下:步骤一:收集能体现工业过程正常运行时的工作状态的数据;步骤二:建立改进主元分析模型对收集的过程数据进行标准化处理,对某一过程变量x的数据进行标准化处理,得到的处理结果为y,<maths num="0001"><![CDATA[<math><mrow><mi>y</mi><mo>=</mo><mfrac><mrow><mi>x</mi><mo>-</mo><mi>&mu;</mi></mrow><mi>&sigma;</mi></mfrac></mrow></math>]]></maths>其中,μ为过程变量x的均值,σ为过程变量x的标准差;将标准化后的过程数据放入二维矩阵Y<sub>m×n</sub>中,其中,m是过程数据的组数,n是过程变量个数;然后对矩阵Y进行主元分析,提取主元方向,并计算对应的特征值,得到全部特征向量构成的矩阵U=[u<sub>1</sub> u<sub>2</sub> …u<sub>n</sub>]和对应的特征值向量λ=[λ<sub>1</sub> λ<sub>2</sub> … λ<sub>n</sub>],特征值矩阵Λ=diag(λ)。设主元个数为k,载荷矩阵P=[u<sub>1</sub> u<sub>2</sub> … u<sub>k</sub>];计算各变量x与主元的复相关系数;根据变量与主元之间的复相关系数的强弱,将变量分为两组,设相关性强的变量个数为s,该组记为PV,另一组记为CV;步骤三:计算统计指标PVR和CVR的低误报率控制限定义两个矩阵:<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>&Phi;</mi><mo>=</mo><mi>U</mi><msup><mi>&Lambda;</mi><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><msup><mi>U</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>I</mi><mo>-</mo><mi>PP</mi><mo>)</mo></mrow><msub><mi>D</mi><mi>PV</mi></msub><mrow><mo>(</mo><mi>I</mi><mo>-</mo><mi>PP</mi><mo>)</mo></mrow><mi>U</mi><msup><mi>&Lambda;</mi><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><msup><mi>U</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><mi>&Theta;</mi><mo>=</mo><mi>U</mi><msup><mi>&Lambda;</mi><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><msup><mi>U</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>I</mi><mo>-</mo><mi>PP</mi><mo>)</mo></mrow><msub><mi>D</mi><mi>CV</mi></msub><mrow><mo>(</mo><mi>I</mi><mo>-</mo><mi>PP</mi><mo>)</mo></mrow><mi>U</mi><msup><mi>&Lambda;</mi><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><msup><mi>U</mi><mo>&prime;</mo></msup></mtd></mtr></mtable></mfenced></math>]]></maths>其中,<img file="FDA00002898933000013.GIF" wi="294" he="73" /><img file="FDA00002898933000014.GIF" wi="352" he="76" />I为单位阵;选取置信度α,通常选择95%或99%,计算PVR和CVR对应的控制限为:<maths num="0003"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&delta;</mi><mi>PVR</mi><mn>2</mn></msubsup><mo>=</mo><msub><mi>&eta;</mi><mn>1</mn></msub><msup><mrow><mo>[</mo><mfrac><mrow><msub><mi>c</mi><mi>&alpha;</mi></msub><msqrt><mn>2</mn><msub><mi>&eta;</mi><mn>2</mn></msub><msubsup><mi>h</mi><mn>1</mn><mn>2</mn></msubsup></msqrt></mrow><msub><mi>&eta;</mi><mn>1</mn></msub></mfrac><mo>+</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>&eta;</mi><mn>2</mn></msub><msub><mi>h</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>h</mi><mn>1</mn></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><msubsup><mi>&eta;</mi><mn>1</mn><mn>2</mn></msubsup></mfrac><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><msub><mi>h</mi><mn>1</mn></msub></mrow></msup></mtd></mtr><mtr><mtd><msubsup><mi>&delta;</mi><mi>CVR</mi><mn>2</mn></msubsup><mo>=</mo><msub><mi>&mu;</mi><mn>1</mn></msub><msup><mrow><mo>[</mo><mfrac><mrow><msub><mi>c</mi><mi>&alpha;</mi></msub><msqrt><mn>2</mn><msub><mi>&mu;</mi><mn>2</mn></msub><msubsup><mi>h</mi><mn>1</mn><mn>2</mn></msubsup></msqrt></mrow><msub><mi>&mu;</mi><mn>1</mn></msub></mfrac><mo>+</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>2</mn></msub><msub><mi>h</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>h</mi><mn>1</mn></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><msubsup><mi>&mu;</mi><mn>1</mn><mn>2</mn></msubsup></mfrac><mo>]</mo></mrow><mrow><mn>1</mn><mo>/</mo><msub><mi>h</mi><mn>2</mn></msub></mrow></msup></mtd></mtr></mtable></mfenced></math>]]></maths>其中,<maths num="0004"><![CDATA[<math><mrow><msub><mi>&eta;</mi><mi>j</mi></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></msubsup><msubsup><mi>&sigma;</mi><mi>i</mi><mi>j</mi></msubsup><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mi>j</mi></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>q</mi></msubsup><msubsup><mi>&delta;</mi><mi>i</mi><mi>j</mi></msubsup><mo>,</mo></mrow></math>]]></maths>j=1,2,3,<maths num="0006"><![CDATA[<math><mrow><msub><mi>h</mi><mn>1</mn></msub><mo>=</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mi>&eta;</mi><mn>1</mn></msub><msub><mi>&eta;</mi><mn>3</mn></msub><mo>/</mo><mrow><mo>(</mo><mn>3</mn><msubsup><mi>&eta;</mi><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mi>h</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mi>&mu;</mi><mn>1</mn></msub><msub><mi>&mu;</mi><mn>3</mn></msub><mo>/</mo><mrow><mo>(</mo><mn>3</mn><msubsup><mi>&mu;</mi><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>σ<sub>i</sub>为矩阵Φ的非零特征值,p为σ<sub>i</sub>的个数,σ<sub>i</sub>为矩阵Θ的非零特征值,q为δ<sub>i</sub>的个数,c<sub>α</sub>为标准正态分布在置信度为α下的临界值;步骤四:使用改进主元分析模型进行过程监视过程监测系统在线监测过程时,使用收集的过程数据计算改进主元分析的统计指标PVR和CVR并与步骤三中计算得到的低误报率控制限进行比对;当统计指标超过对应的控制限时,系统报警。
地址 102206 北京市昌平区朱辛庄北农路2号