发明名称 一种汽轮机系统流量计故障监测方法
摘要 一种汽轮机系统流量计故障监测方法,属于汽轮机在线监测领域。本法利用数据协调方法,对故障流量计进行检测和识别。具体技术步骤包括:1)利用稳态工况测量数据,并根据测量参数测量值的不确定度及系统约束方程进行数据协调计算;2)比较数据协调最小目标函数值与故障检测阈值的大小,判断是否发生流量计故障;3)若检测出流量计发生故障,则逐个对流量计进行故障识别,将其对应的测量参数当作未测量参数,重新计算数据协调最小目标函数值;比较新的最小目标函数值和故障检测阈值,从而识别故障流量计。该方法可有效利用汽轮机系统中冗余的流量测量信息,进行流量计故障的检测与识别,方法简单易于实现,具有高效低成本的优点。
申请公布号 CN103674189B 申请公布日期 2016.06.15
申请号 CN201310596063.3 申请日期 2013.11.21
申请人 清华大学 发明人 李政;蒋晓隆;刘培;仲晓波;朱泓逻
分类号 G01F25/00(2006.01)I 主分类号 G01F25/00(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 一种汽轮机系统流量计故障监测方法,其特征在于该方法按如下步骤进行:1)从电厂在线数据库采集稳态工况下的压力、温度和流量测量参数,测量参数包括:高压给水加热器抽汽入口的压力P<sub>1</sub>,P<sub>2</sub>,…P<sub>N</sub>、给水泵出口的压力P<sub>N+1</sub>、除氧器抽汽入口的压力P<sub>N+2</sub>、除氧器凝结水入口的压力P<sub>N+3</sub>,其中N为高压给水加热器的个数;高压给水加热器抽汽入口的温度T<sub>1</sub>,T<sub>2</sub>,…T<sub>N</sub>、高压给水加热器给水出口的温度T<sub>N+1</sub>,T<sub>N+2</sub>,…T<sub>2N</sub>、高压给水加热器疏水出口的温度T<sub>2N+1</sub>,T<sub>2N+2</sub>,…T<sub>3N</sub>、给水泵出口的温度T<sub>3N+1</sub>、除氧器抽汽入口的温度T<sub>3N+2</sub>、除氧器凝结水入口的温度T<sub>3N+3</sub>;最终给水的流量、给水泵出口的流量和除氧器凝结水入口的流量D<sub>1</sub>,D<sub>2</sub>,D<sub>3</sub>;上述压力、温度和流量测量参数依次记为X<sub>1</sub>,X<sub>2</sub>,…X<sub>n</sub>,测量参数的测量值记为x<sub>1</sub>,x<sub>2</sub>,…x<sub>n</sub>,其中n为所有测量参数的个数;2)根据汽轮机系统中压力传感器、温度传感器和流量计的精度等级,评定测量参数测量值的不确定度,分别记为σ<sub>1</sub>,σ<sub>2</sub>,…σ<sub>n</sub>;3)系统中未测量参数包括:高压给水加热器抽汽的流量d<sub>1</sub>,d<sub>2</sub>,…d<sub>N</sub>、高压给水加热器疏水的流量d<sub>N+1</sub>,d<sub>N+2</sub>,…d<sub>2N</sub>、除氧器抽汽的流量d<sub>2N+1</sub>、除氧器给水出口的流量d<sub>2N+2</sub>;上述未测量参数依次记为Y<sub>1</sub>,Y<sub>2</sub>,…Y<sub>p</sub>,其中p为所有未测量参数的个数;4)利用汽轮机系统中高压给水加热器、给水泵和除氧器的质量平衡和能量平衡关系,并忽略高压给水加热器中给水侧和抽汽侧的压降,建立关于测量参数X<sub>1</sub>,X<sub>2</sub>,…X<sub>n</sub>和未测量参数Y<sub>1</sub>,Y<sub>2</sub>,…Y<sub>p</sub>的约束方程组f<sub>m</sub>,其中m为约束方程的个数:f<sub>m</sub>(X<sub>1</sub>,X<sub>2</sub>,…X<sub>n</sub>,Y<sub>1</sub>,Y<sub>2</sub>,…Y<sub>p</sub>)=0   (1)5)利用测量参数的测量值x<sub>1</sub>,x<sub>2</sub>,…x<sub>n</sub>、测量值的不确定度σ<sub>1</sub>,σ<sub>2</sub>,…σ<sub>n</sub>以及约束方程组f<sub>m</sub>进行数据协调计算,使得数据协调目标函数<img file="FDA0000939893600000011.GIF" wi="43" he="54" />最小,并且测量参数的协调值<img file="FDA0000939893600000012.GIF" wi="218" he="70" />和未测量参数的估计值y<sub>1</sub>,y<sub>2</sub>,…y<sub>p</sub>满足系统约束方程组:<img file="FDA0000939893600000013.GIF" wi="982" he="151" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><msub><mi>f</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mo>,</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>2</mn></msub><mo>,</mo><mo>...</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>n</mi></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>...</mo><msub><mi>y</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000939893600000014.GIF" wi="789" he="87" /></maths>6)将求解得到的数据协调最小目标函数值<img file="FDA0000939893600000015.GIF" wi="85" he="55" />与流量计故障检测阈值<img file="FDA0000939893600000016.GIF" wi="106" he="78" />进行比较;检测阈值<img file="FDA0000939893600000017.GIF" wi="107" he="77" />是自由度为r的卡方分布的1‑α分位数;r为系统冗余度,等于约束方程个数m减去未测量参数的个数p,r=m‑p;α为流量计故障检测的显著性水平;如果<img file="FDA0000939893600000021.GIF" wi="269" he="77" />表明流量计没有发生故障;7)如果<img file="FDA0000939893600000022.GIF" wi="270" he="79" />表明流量计发生故障,则逐个对流量计进行故障识别;当对第j个流量计进行故障识别时,将该流量计对应的测量参数X<sub>j</sub>当作未测量参数Y<sub>p+1</sub>,重新进行数据协调计算:<img file="FDA0000939893600000023.GIF" wi="1349" he="152" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><msub><mi>f</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mo>,</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>2</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>x</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>n</mi></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>...</mo><msub><mi>y</mi><mi>p</mi></msub><mo>,</mo><msub><mi>y</mi><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000939893600000024.GIF" wi="1149" he="91" /></maths>将求解得到的新的数据协调最小目标函数值<img file="FDA00009398936000000212.GIF" wi="104" he="53" />与新的流量计故障检测阈值<img file="FDA0000939893600000025.GIF" wi="137" he="78" />进行比较;如果<img file="FDA0000939893600000026.GIF" wi="318" he="79" />表明除第j个流量计以外的其它流量计没有发生故障,则识别出发生故障的流量计为第j个流量计,故障识别结束;如果<img file="FDA0000939893600000027.GIF" wi="323" he="78" />表明除第j个流量计以外的其它流量计仍然有故障,继续对下一个流量计进行故障识别;8)按步骤7)逐个对流量计进行故障识别时,如果所有识别结果均为<img file="FDA0000939893600000028.GIF" wi="318" he="79" />表明有多个流量计发生故障;则首先判定第一轮故障识别过程中<img file="FDA0000939893600000029.GIF" wi="121" he="63" />最小的流量计为故障流量计,并将该流量计对应的测量参数X<sub>j</sub>当作未测量参数Y<sub>p+1</sub>,再按步骤7)逐个对其余流量计进行下一轮故障识别,其中,第R轮流量计故障识别的故障检测阈值为<img file="FDA00009398936000000210.GIF" wi="167" he="76" />9)当完成的故障识别的轮数R=r‑1时,流量计故障识别结束。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室