发明名称 一种基于故障树信息的设备状态预测方法
摘要 本发明公开了一种基于故障树信息的设备状态预测方法,首先搜索故障树模型,将模型中的各事件用不同变量表示,将模型中的逻辑门用变量的结构关系集合表示,依次计算故障预测模型中故障原因变量先验概率分布集合,计算故障预测模型中故障模式变量的条件概率分布和故障传递变量的条件概率分布,将设备实际故障检测现象用故障预测模型中对应的故障检测变量表示,根据故障检测变量结构关系集合中表述的关联关系计算故障预测模型中故障检测变量的条件概率分布,根据现场采集到的故障检测现象实时信息计算设备故障模式后验概率状态分布。本发明能够快速准确的预测设备实时状态、指导对设备的监控及维修,有效提高维修效率、降低维修成本。
申请公布号 CN101950327B 申请公布日期 2012.05.23
申请号 CN201010277907.4 申请日期 2010.09.09
申请人 西北工业大学 发明人 司书宾;蔡志强;孙树栋;王宁;兑红炎;李淑敏;张丽丽
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 1.一种基于故障树信息的设备状态预测方法,其特征在于包括下述步骤:1)确定需要预测的一个设备故障模式,并在故障树中搜索以该故障模式为顶事件的故障树模型F,F中所有事件集合为T={T<sub>1</sub>,T<sub>2</sub>,...,T<sub>i</sub>,...,T<sub>n</sub>},其中T<sub>i</sub>表示第i个事件,n表示故障树事件总数;故障树具体搜索方式如下:针对设备故障模式,首先在故障树信息中找到与该模式对应的故障树事件T<sub>1</sub>,然后以该事件为顶事件向下搜索与顶事件通过逻辑门关联的故障树中间事件,不断向下搜索与中间事件通过逻辑门关联的其它故障树中间事件,直到搜索到故障树底事件为止;其中,底事件是指没有其它故障树事件通过逻辑门与其关联的事件;2)将集合T中的顶事件用故障预测模型中的故障模式变量M表示,当顶事件不发生时,对应变量的取值为0,当顶事件发生时,对应变量的取值为1,且整个模型中只有一个故障模式变量;将集合T中的其它中间事件用故障预测模型中的故障传递变量集合D={D<sub>1</sub>,D<sub>2</sub>,...,D<sub>j</sub>,...,D<sub>m</sub>}表示,其中D<sub>j</sub>表示第j个中间事件对应的故障传递变量,<img file="FSA00000262447500011.GIF" wi="790" he="136" />j=1,...,m,m表示故障传递变量总数;将集合T中的底事件用故障预测模型中的故障原因变量集合C={C<sub>1</sub>,C<sub>2</sub>,...,C<sub>k</sub>,...,C<sub>l</sub>}表示,其中C<sub>k</sub>表示第k个底事件对应的故障原因变量,<img file="FSA00000262447500012.GIF" wi="740" he="137" />k=1,...,l,l表示故障原因变量总数,且1+m+l=n;3)将故障树模型F中表示的逻辑门集合L={L<sub>0</sub>,L<sub>1</sub>,L<sub>2</sub>,...,L<sub>j</sub>,...,L<sub>m</sub>}用故障预测模型中变量的结构关系集合J={J<sub>0</sub>,J<sub>1</sub>,J<sub>2</sub>,...,J<sub>j</sub>,...,J<sub>m</sub>}表示,其中L<sub>0</sub>表示顶事件对应的逻辑门,L<sub>j</sub>表示第j个中间事件对应的逻辑门,故障树中共有1+m个逻辑门;J<sub>0</sub>表示故障模式变量M对应的结构关系,J<sub>j</sub>表示第j个传递变量对应的结构关系;在<img file="FSA00000262447500013.GIF" wi="268" he="76" />中,π(D<sub>j</sub>)表示影响第j个中间事件状态的故障树事件相对应的变量,(·)表示与第j个逻辑门中逻辑关系相对应的关联关系;根据逻辑门建立变量结构关系的具体方式如下:<img file="FSA00000262447500021.GIF" wi="807" he="393" />其中j=0,...,m其中,<img file="FSA00000262447500022.GIF" wi="191" he="74" />表示π(D<sub>j</sub>)中的变量为并联结构关系;<img file="FSA00000262447500023.GIF" wi="193" he="76" />表示π(D<sub>j</sub>)中的变量为串联结构关系;用<img file="FSA00000262447500024.GIF" wi="185" he="76" />表示π(D<sub>j</sub>)中变量与D<sub>j</sub>的结构关系,命名为反联结构;用<img file="FSA00000262447500025.GIF" wi="197" he="76" />表示π(D<sub>j</sub>)中变量与D<sub>j</sub>的结构关系,命名为多连结构,且在r个影响D<sub>j</sub>状态的变量中需要至少有q个同时故障时,D<sub>j</sub>才会发生故障;4)根据故障树模型包含的各底事件可靠性函数集合R={R<sub>1</sub>,R<sub>2</sub>,...,R<sub>k</sub>,...,R<sub>l</sub>}计算故障预测模型中故障原因变量先验概率分布集合P={P(C<sub>1</sub>=0),P(C<sub>2</sub>=0),...,P(C<sub>k</sub>=0),...,P(C<sub>l</sub>=0)},P(C<sub>k</sub>=0)=R<sub>k</sub>,其中k=1,...,l,其中R<sub>k</sub>表示第k个底事件的可靠性函数,P(C<sub>k</sub>=0)表示第k个故障原因变量的状态概率分布;5)根据结构关系集合J={J<sub>0</sub>,J<sub>1</sub>,J<sub>2</sub>,...,J<sub>j</sub>,...,J<sub>m</sub>}中所包含的关联关系,计算故障预测模型中故障模式变量的条件概率分布P(M|π(M))和故障传递变量的条件概率分布CP<sub>D</sub>={P(D<sub>1</sub>|π(D<sub>1</sub>)),P(D<sub>2</sub>|π(D<sub>2</sub>)),...,P(D<sub>j</sub>|π(D<sub>j</sub>)),..,P(D<sub>m</sub>|π(D<sub>m</sub>))},其中P(D<sub>j</sub>|π(D<sub>j</sub>))表示第j个故障传递变量的条件概率分布,传递变量条件概率分布的总数也是l个;根据结构关系建立变量条件概率分布的具体方式如下:当π(D<sub>j</sub>)中变量为并联结构关系,<img file="FSA00000262447500026.GIF" wi="928" he="137" />当π(D<sub>j</sub>)中变量为串联结构关系,<img file="FSA00000262447500027.GIF" wi="921" he="129" />当π(D<sub>j</sub>)中变量与D<sub>j</sub>为反联结构关系,<img file="FSA00000262447500028.GIF" wi="819" he="137" />当π(D<sub>j</sub>)中变量为多联结构关系,<img file="FSA00000262447500029.GIF" wi="991" he="131" />6)将设备实际故障检测现象用故障预测模型中对应的故障检测变量E={E<sub>1</sub>,E<sub>2</sub>,...,E<sub>x</sub>,...,E<sub>p</sub>}表示,其中E<sub>x</sub>表示第x个故障检测现象对应的故障检测变量,<img file="FSA00000262447500031.GIF" wi="835" he="137" />x=1,...,p,p表示故障检测变量总数;按<img file="DEST_PATH_RE-FSB00000622515700012.GIF" wi="1077" he="161" />其中x=1,...,p的方式建立故障检测变量的结构关系集合G={G<sub>1</sub>,G<sub>2</sub>,...,G<sub>x</sub>,...,G<sub>p</sub>},其中G<sub>x</sub>表示第x个故障检测变量所包含的结构关系,结构关系的总数也为p个;在<img file="DEST_PATH_RE-FSB00000622515700013.GIF" wi="272" he="71" />中,π(E<sub>x</sub>)表示影响第x个故障检测现象状态的故障原因相对应的变量,(·)表示第x个故障检测现象及其对应故障原因间的关联关系;<img file="DEST_PATH_RE-FSB00000622515700014.GIF" wi="193" he="71" />表示π(E<sub>x</sub>)中的变量为并联结构关系;<img file="DEST_PATH_RE-FSB00000622515700015.GIF" wi="193" he="69" />表示π(E<sub>x</sub>)中的变量为串联结构关系;7)根据故障检测变量结构关系集合G中表述的关联关系,分别按π(D<sub>j</sub>)中变量为并联结构关系和串联结构关系的方式计算故障预测模型中故障检测变量的条件概率分布CP<sub>E</sub>={P(E<sub>1</sub>|π(E<sub>1</sub>)),P(E<sub>2</sub>|π(E<sub>2</sub>)),...,P(E<sub>x</sub>|π(E<sub>x</sub>)),...,P(E<sub>p</sub>|π(E<sub>p</sub>))},其中P(E<sub>x</sub>|π(E<sub>x</sub>))表示第x个故障检测变量的条件概率分布;8)根据现场采集到的故障检测现象实时信息e={e<sub>1</sub>,e<sub>2</sub>,...,e<sub>x</sub>,...,e<sub>p</sub>},其中e<sub>x</sub>表示第x个故障检测现象的实际状态信息,利用故障预测模型推理能力和条件概率分布,基于贝叶斯定理,计算设备故障模式后验概率状态分布P(M=0|E=e)=P(M=0|π(E))×P(π(E)|E=e)。 
地址 710072 陕西省西安市友谊西路127号