发明名称 基于定量结构模型的故障诊断系统传感器优化配置方法
摘要 本发明涉及一种基于定量结构模型的故障诊断系统传感器优化配置方法,第一步,运用结构分析方法建立故障诊断系统的结构模型。该结构模型的表现形式为两个矩阵:解析冗余关系矩阵和故障特征矩阵。该结构模型能够定性地描述故障诊断系统。第二步,结合故障诊断系统结构模型和传感器检测不确定性,得到故障诊断系统的定量描述,达到定量分析故障诊断系统的目的,进而建立考虑检测不确定性的故障诊断系统关键性能指标的计算公式。第三步建立故障诊断系统的传感器优化模型。最后一步,采用遗传算法对优化模型进行求解。本发明提出的面向故障诊断系统的传感器优化配置方法可以有效节约故障诊断系统的传感器成本且配置结果更加符合工程实际情况。
申请公布号 CN105678423A 申请公布日期 2016.06.15
申请号 CN201610059103.4 申请日期 2016.01.28
申请人 西北工业大学 发明人 谢蓉;孙无忌;李婷;王剑;王新民
分类号 G06Q10/04(2012.01)I;G06N3/12(2006.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于定量结构模型的故障诊断系统传感器优化配置方法,其特征在于步骤如下:步骤一 建立由解析冗余关系矩阵B和故障特征矩阵M描述的故障诊断系统的结构模型:1.系统的解析冗余关系集合为A={a<sub>1</sub>,a<sub>2</sub>,...a<sub>m</sub>}、完备传感器集合为S={s<sub>1</sub>,s<sub>2</sub>,...,s<sub>n</sub>}和故障集合为F={f<sub>1</sub>,f<sub>2</sub>,...,f<sub>p</sub>},a<sub>i</sub>表示第i个解析冗余关系,s<sub>j</sub>表示第j个传感器配置位置,f<sub>k</sub>表示第k类故障;2.建立解析冗余关系矩阵B=[b<sub>ij</sub>]<sub>m×n</sub>,矩阵的列表示可测变量或传感器,行表示系统的解析冗余关系,b<sub>ij</sub>=1表示解析冗余关系式a<sub>i</sub>中含有可测变量s<sub>j</sub>,b<sub>ij</sub>=0则表示解析冗余关系式a<sub>i</sub>中不包含可测变量s<sub>j</sub>;3.建立故障特征矩阵M=[m<sub>jk</sub>]<sub>m×p</sub>,表示任一故障f<sub>k</sub>的故障特征,m<sub>ik</sub>=1表示故障f<sub>k</sub>的出现会影响到解析冗余关系a<sub>i</sub>,此时a<sub>i</sub>不是满足的,m<sub>ik</sub>=0表示故障f<sub>k</sub>的出现不会影响到解析冗余关系a<sub>i</sub>;步骤二 求解考虑检测不确定性的故障诊断系统关键性能指标—故障检测率FDR和故障隔离率FIR:1.故障诊断系统的可检测故障集合F<sub>D</sub>和可隔离故障集合F<sub>I</sub>,故障诊断系统性能指标—故障检测率FDR和故障隔离率FIR的计算公式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>F</mi><mi>D</mi><mi>R</mi><mo>=</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mi>k</mi></msub><mo>&Element;</mo><msub><mi>F</mi><mi>D</mi></msub></mrow></munder><msub><mi>&lambda;</mi><mi>k</mi></msub><msub><mi>Q</mi><mi>k</mi></msub><mo>/</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mi>k</mi></msub><mo>&Element;</mo><mi>F</mi></mrow></munder><msub><mi>&lambda;</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0000916591690000011.GIF" wi="509" he="126" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>F</mi><mi>I</mi><mi>R</mi><mo>=</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mi>k</mi></msub><mo>&Element;</mo><msub><mi>F</mi><mi>I</mi></msub></mrow></munder><msub><mi>&lambda;</mi><mi>k</mi></msub><msub><mi>Q</mi><mi>k</mi></msub><mo>/</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mi>k</mi></msub><mo>&Element;</mo><msub><mi>F</mi><mi>D</mi></msub></mrow></munder><msub><mi>&lambda;</mi><mi>k</mi></msub><msub><mi>Q</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0000916591690000012.GIF" wi="549" he="119" /></maths>其中,λ为故障率,即单位时间内发生的故障数;Q<sub>k</sub>为故障f<sub>k</sub>总的可检测率,计算公式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>k</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><munderover><mo>&Pi;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>R</mi><mi>i</mi></msub><mo>&times;</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000916591690000013.GIF" wi="478" he="134" /></maths>其中,R<sub>i</sub>表示第i个解析冗余关系的检测可靠性;d<sub>ik</sub>为解析冗余关系a<sub>i</sub>对故障f<sub>k</sub>的可测性;m是解析冗余关系的个数;所述解析冗余关系检测可靠性R<sub>i</sub>的计算公式为:<img file="FDA0000916591690000021.GIF" wi="389" he="134" />其中,R<sub>i</sub>表示第i个解析冗余关系的检测可靠性;r<sub>j</sub>表示第j个传感器的失效率;x<sub>j</sub>是第j个传感器的配置个数;n是可测变量的个数;b<sub>ij</sub>是解析冗余关系矩阵中的元素。根据计算公式,可以得到与解析冗余关系集合A={a<sub>1</sub>,a<sub>2</sub>,...a<sub>m</sub>}相对应的检测可靠性矢量DR=[R<sub>1</sub>,R<sub>2</sub>,...R<sub>m</sub>];所述故障可测性d<sub>ik</sub>的计算公式为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>TTD</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow><mrow><msub><mi>TTF</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mn>0.5</mn></msup><mo>&times;</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>SyD</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow><mrow><msub><mi>TTF</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mn>0.2</mn></msup></mrow></mtd><mtd><mrow><msub><mi>TTD</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>&lt;</mo><msub><mi>TTF</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msub><mi>TTD</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>&GreaterEqual;</mo><msub><mi>TTF</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000916591690000022.GIF" wi="990" he="231" /></maths>其中,d<sub>ik</sub>为解析冗余关系a<sub>i</sub>对故障f<sub>k</sub>的可测性;TTD<sub>ik</sub>为故障f<sub>k</sub>初始时刻到解析冗余关系a<sub>i</sub>检测到该故障的时间跨度;SyD<sub>ik</sub>为解析冗余关系a<sub>i</sub>检测到故障f<sub>k</sub>的持续时间;TTF<sub>ik</sub>为故障f<sub>k</sub>初始时刻到系统失效时的时间跨度;<img file="FDA0000916591690000023.GIF" wi="137" he="142" />和<img file="FDA0000916591690000024.GIF" wi="134" he="135" />两个比值描述了解析冗余关系a<sub>i</sub>对故障f<sub>k</sub>的故障检测时效性和故障可跟踪性;D=[d<sub>jk</sub>]<sub>m×p</sub>为计算得到故障可测性矩阵;步骤三 基于故障诊断系统的定量结构模型:将故障检测率FDR和故障隔离率FIR作为约束条件,经济成本作为优化目标,建立考虑检测不确定性时的传感器优化模型:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open = "{" close = "}"><mtable><mtr><mtd><mrow><mover><mi>S</mi><mo>^</mo></mover><mo>=</mo><munder><mrow><mi>arg</mi><mi> </mi><mi>max</mi></mrow><mrow><msup><mi>S</mi><mo>*</mo></msup><mo>&SubsetEqual;</mo><mi>S</mi></mrow></munder><mi>O</mi><mi>b</mi><mi>j</mi><mrow><mo>(</mo><msup><mi>S</mi><mo>*</mo></msup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>O</mi><mi>b</mi><mi>j</mi><mrow><mo>(</mo><msup><mi>S</mi><mo>*</mo></msup><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mn>1</mn></msub><mo>&CenterDot;</mo><mn>1</mn><mo>/</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>&Element;</mo><msup><mi>S</mi><mo>*</mo></msup></mrow></munder><msub><mi>c</mi><mi>j</mi></msub><msub><mi>x</mi><mi>j</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>F</mi><mi>D</mi><mi>R</mi><mo>&GreaterEqual;</mo><msup><mi>FDR</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>F</mi><mi>I</mi><mi>R</mi><mo>&GreaterEqual;</mo><msup><mi>FIR</mi><mo>&prime;</mo></msup></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000916591690000025.GIF" wi="550" he="599" /></maths>其中,S*为已配置传感器集合;<img file="FDA0000916591690000026.GIF" wi="182" he="63" />Obj(S<sup>*</sup>)是以S*为变量的目标函数;C<sub>1</sub>为加权值;c<sub>j</sub>为第j个传感器的成本;FDR′,FIR′是故障诊断系统必须满足的指标要求;<img file="FDA0000916591690000031.GIF" wi="45" he="78" />为最优的传感器配置集合;将故障检测率FDR和故障隔离率FIR作为约束条件,经济成本作为优化目标,得到约束条件下成本最低的传感器布局方案,即目标函数Obj(S<sup>*</sup>)取最大值时的解<img file="FDA0000916591690000032.GIF" wi="63" he="87" />步骤四 根据遗传算法求解上述传感器优化模型,得到最优传感器配置集合,具体过程如下:1.参数初始化,随机产生初始种群;2.根据传感器优化模型,建立适应度函数:<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>F</mi><mi>i</mi><mi>t</mi><mi>F</mi><mi>u</mi><mi>n</mi><mo>=</mo><msub><mi>C</mi><mn>1</mn></msub><mo>&times;</mo><mfrac><mn>1</mn><mrow><mi>C</mi><mi>o</mi><mi>s</mi><mi> </mi><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><msub><mi>C</mi><mn>2</mn></msub><mo>&times;</mo><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mn>0</mn><mo>,</mo><msup><mi>FDR</mi><mo>&prime;</mo></msup><mo>-</mo><mi>F</mi><mi>D</mi><mi>R</mi><mo>(</mo><mi>i</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>C</mi><mn>3</mn></msub><mo>&times;</mo><mi>max</mi><mrow><mo>(</mo><mn>0</mn><mo>,</mo><msup><mi>FIR</mi><mo>&prime;</mo></msup><mo>-</mo><mi>F</mi><mi>I</mi><mi>R</mi><mo>(</mo><mi>i</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000916591690000033.GIF" wi="1053" he="214" /></maths>其中,Cost(i)为第i个传感器配置方案的成本;max(0,FDR′‑FDR(i))和max(0,FIR′‑FIR(i))为惩罚函数,即当第i个传感器配置方案的故障检测率FDR(i)和故障隔离率FIR(i)小于指标要求FDR′和FIR′时,该方案的适应度值会减少;C<sub>1</sub>、C<sub>2</sub>、C<sub>3</sub>为常数;3.计算每个染色体所对应的传感器集合的适应度值;4.根据适应度采用轮盘赌方法选择个体,并进行交叉操作,得到当前种群Pop';5.对当前种群Pop'中的个体进行变异操作,得到下一代种群Pop”,判断算法是否到达最大迭代次数,是则结束,运算结果即为最优传感器配置方案,否则返回步骤3。
地址 710072 陕西省西安市友谊西路127号