发明名称 一种基于数据协调及假设检验的电站系统多故障诊断方法
摘要 一种基于数据协调及假设检验的电站系统多故障诊断方法,属于电厂性能监测及故障诊断领域。具体技术步骤包括:1)选取准静态测量数据,建立系统约束方程,进行数据协调计算;2)计算全局假设检验统计量,和允许的临界值对比,判断参数是否存在显著误差;3)若存在显著误差,识别包含误差的参数并在下一轮计算前剔除;4)计算降秩后的检验统计量,并和新的临界值对比,判断参数是否存在显著误差;5)若不存在,则结束运算;否则重复步骤3)及4),直至检测出所有的显著误差。该方法采用序列消去的思想,逐一实现快速、精准的故障定位。显著误差临界值可作为电站报警阈值,在故障早期及时提示操作人员,从而保护机组安全,降低生产成本。
申请公布号 CN106199241A 申请公布日期 2016.12.07
申请号 CN201610473991.4 申请日期 2016.06.24
申请人 清华大学 发明人 李政;刘培;郭思思
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 一种基于数据协调及假设检验的电站系统多故障诊断方法,其特征是该方法包括如下步骤:1)通过准静态数据分析,从电厂数据库选取准静态工况的测量数据;根据电厂系统的测点信息,选取参与热力系统平衡计算的测量变量x及未测量变量u,并将测量变量x的测量值记为y;<maths num="0001"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>y</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>y</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>u</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>p</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001029157940000011.GIF" wi="630" he="175" /></maths>其中n为所有测量变量的个数;p为所有未测量变量的个数;2)根据电厂热力系统中测量仪表的精度等级,按照下列方程(1)评定各个测量变量测量值的不确定度,分别记为σ<sub>1</sub>,σ<sub>2</sub>,…σ<sub>n</sub>;<maths num="0002"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>k&xi;</mi><mi>i</mi></msub></mrow><mrow><mn>1.96</mn><msqrt><msub><mi>N</mi><mi>s</mi></msub></msqrt></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001029157940000012.GIF" wi="1277" he="108" /></maths>其中,i代表第i个测量变量,i取1,2,3...n;ξ<sub>i</sub>代表置信区间为95%的仪表最大允许误差;N<sub>s</sub>代表测量同一变量的仪表个数;k为修正系数;3)根据电厂系统热力设备的平衡方程及特性方程,建立关于测量变量x及未测量变量u的约束方程组,如下所示:f(x,u)=0                                  (2)g(x,u)≤0                                  (3)其中,f代表等式约束方程组;g代表不等式约束方程组;4)结合测量变量的测量值、测量值的不确定度、约束方程组f和g,利用目标函数<img file="FDA0001029157940000019.GIF" wi="38" he="47" />计算得到测量变量x的协调值x<sup>*</sup>和未测量变量u的估计值u<sup>*</sup>;<img file="FDA0001029157940000013.GIF" wi="1286" he="63" />其中,<img file="FDA0001029157940000014.GIF" wi="511" he="183" />5)数据协调计算的最小目标函数值<img file="FDA0001029157940000015.GIF" wi="104" he="54" />等于全局假设检验统计量γ<sub>R</sub>;假设当测量变量中不存在显著误差时,统计分布的临界值为<img file="FDA0001029157940000016.GIF" wi="162" he="71" />将γ<sub>R</sub>与统计分布的临界值<img file="FDA0001029157940000017.GIF" wi="130" he="71" />进行比较;其中,R代表系统的冗余度,α为显著性水平;6)如果<img file="FDA0001029157940000018.GIF" wi="267" he="70" />则表明测量变量中不存在显著误差,即仪表或设备没有发生故障;相反,若<img file="FDA0001029157940000021.GIF" wi="297" he="71" />则表明测量变量中存在显著误差,需要进行故障识别;当对第j个测量变量进行故障识别时,将该测量变量当作未测量变量,则系统的冗余度减一,按照下式重新进行数据协调计算:<img file="FDA0001029157940000022.GIF" wi="1338" he="70" />其中,<img file="FDA0001029157940000023.GIF" wi="53" he="55" />为消去第j个测量变量后的目标函数;j的取值依次从1到n;此时,新的全局假设检验统计量γ<sub>R‑1</sub>等于数据协调计算得到的<img file="FDA0001029157940000024.GIF" wi="51" he="53" />的最小值;将γ<sub>R‑1</sub>与新的统计分布的临界值<img file="FDA0001029157940000025.GIF" wi="175" he="71" />进行比较,<img file="FDA0001029157940000026.GIF" wi="181" he="70" />代表系统冗余度降为R‑1后的统计分布的临界值;如果<img file="FDA0001029157940000027.GIF" wi="363" he="70" />表明除第j个测量变量以外的其它变量没有发生故障,则识别出第j个测量变量发生故障,故障识别结束;如果<img file="FDA0001029157940000028.GIF" wi="363" he="70" />表明除第j个测量变量以外的其它测量变量仍然有故障,则继续对下一个测量变量进行故障识别;7)按步骤6)逐个对测量变量进行故障识别时,如果所有识别结果均为<img file="FDA0001029157940000029.GIF" wi="155" he="48" /><img file="FDA00010291579400000210.GIF" wi="203" he="70" />则表明同时有多个测量变量发生故障;找出第一轮故障识别过程中检验统计量γ<sub>R‑1</sub>最小的情况,判定被消去的测量变量存在故障,并将该测量变量作为未测量变量,按步骤6)逐个对其余测量变量进行第二轮故障识别,此时系统的冗余度为R‑2,统计分布的临界值为<img file="FDA00010291579400000211.GIF" wi="195" he="70" />8)在第二轮故障识别过程中,找出第二个包含故障的测量变量;按步骤7)判断是否仍有其他故障存在;如果不存在,则故障识别结束;如果存在,则将第二个包含故障的测量变量作为未测量变量,进行第三轮故障识别;9)以此类推,重复步骤6)、7)、8),当第r轮故障识别过程中,出现<img file="FDA00010291579400000212.GIF" wi="162" he="54" /><img file="FDA00010291579400000213.GIF" wi="176" he="71" />或进行的轮数r=R时,故障识别结束,至此依次识别出测量变量中的全部故障。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室
您可能感兴趣的专利