基于定量结构模型的故障诊断系统传感器优化配置方法
【技术领域】
[0001] 本发明属于故障诊断技术领域，涉及一种基于定量结构模型的故障诊断系统传感 器优化配置方法，特别是面向基于解析冗余关系的故障诊断系统。
【背景技术】
[0002] 在系统的运行过程中如果出现故障，有可能会在经济、产品质量和安全上造成灾 难性的后果，造成严重的损失。为了监测、隔离出这些不正常的系统状态，故障诊断系统变 得越来越重要。基于模型的故障诊断系统通常需要传感器测量系统的重要参数来提取有用 的故障信息。因此，故障诊断系统的性能很大程度上依赖于传感器的配置情况。
[0003] 传感器优化配置问题是故障诊断系统的基础，对故障诊断系统功能的实现至关重 要。传感器的优化配置主要考虑如何在满足故障诊断系统性能需求的前提下选择出最优的 传感器配置方案。本发明提出的传感器优化配置方法是面向基于解析冗余关系故障诊断系 统的。基于解析冗余关系故障诊断方法的基本思想是根据解析冗余关系残差值建立观测特 征，通过观测特征与故障特征的比较进行故障检测和隔离。这种故障诊断系统由于其优良 的故障检测、分离与辨识性能已经受到越来越多的关注。
[0004] "Optimal sensor placement for model-based fault detection and isolation"（Proceedings of the IEEE Conference on Decision and Control,2008: 2584-2589.)-文的传感器优化配置方法研究针对的是基于解析冗余关系的故障诊断系 统。虽然此公开文献给出了一种采用结构模型对故障诊断系统进行传感器优化配置的方 法，但这种方法只是定性地分析系统模型所含公式和变量之间的逻辑关系，假设传感器能 够以100%的概率检测到相应被测量。实际上，传感器自身的故障率以及系统检测故障的能 力会影响测量结果的准确性。因此，由这种方法建立的故障诊断系统模型是一种不考虑检 测不确定性的数学模型，在一定程度上偏离工程实际。虽然《考虑传感器故障检测能力的 PHM系统传感器优化配置方法》(宇航学报，2013,34(09): 1253-1258.)-文给出了一种基于 定量结构模型的故障诊断系统传感器优化配置方法。但这篇文献针对的并不是基于解析冗 余关系的故障诊断系统。目前，在基于解析冗余关系故障诊断系统领域尚未出现考虑检测 不确定性的传感器优化配置方法的发明报道。
【发明内容】
[0005] 要解决的技术问题
[0006] 为了避免现有技术的不足之处，本发明提出一种基于定量结构模型的故障诊断系 统传感器优化配置方法，对故障诊断系统进行性能评估时将传感器故障率和故障检测能力 考虑在内，并根据评估值得到更加准确和可靠的传感器优化配置结果。
[0007] 技术方案
[0008] -种基于定量结构模型的故障诊断系统传感器优化配置方法，其特征在于步骤如 下：
[0009] 步骤一建立由解析冗余关系矩阵B和故障特征矩阵Μ描述的故障诊断系统的结构 模型。
[0010] 1.确定系统的解析冗余关系集合A= {ai，a2，. . . am}、完备传感器(可测变量)集合S ={si，S2, . . .，sn}和故障集合F= {fi，f2, . . .，fp}，ai表示第i个解析冗余关系，Sj表示第j个 传感器配置位置，fk表示第k类故障。
[0011] 2.建立解析冗余关系矩阵Β=[1^]ωΧη，矩阵的列表示可测变量或传感器，行表示 系统的解析冗余关系，b^ = l表示解析冗余关系式&1中含有可测变量~，Μ = 0则表示解析 冗余关系式ai中不包含可测变量建立故障特征矩阵M=[mjk] mXp表示任一故障fk的故障 特征，mik= 1表示故障fk的出现会影响到解析冗余关系ai，此时ai不是满足的，mik = 0表示故 障fk的出现不会影响到解析冗余关系m。
[0012]步骤二求解考虑检测不确定性的故障诊断系统关键性能指标 故障检测率FDR 和故障隔离率FIR。
[0013] 1.确定故障诊断系统的可检测故障集合Fd和可隔离故障集合F:。故障诊断系统性 能指标--故障检测率FDR和故障隔离率FIR的计算公式：
[0014]
[0015] Ik/
[0016] 其中，λ为故障率，g卩单位时间内发生的故障数;Qk为故障fk总的可检测率，计算公 式为：
[0017] …
[0018] 其中，心表示第i个解析冗余关系的检测可靠性;dlk为解析冗余关系&1对故障f k的 可测性;m是解析冗余关系的个数。
[0019] 2.解析冗余关系检测可靠性心和故障可测性dlk是两个用于描述故障检测的不确 定性的参数，其取值可以通过计算得到。这两个参数可以定量地描述故障诊断系统。
[0020] 解析冗会关系检测可靠性Ri的计算公式如下：
[0021]
[0022] 其中，心表示第i个解析冗余关系的检测可靠性;^表示第j个传感器的失效率; 是第j个传感器的配置个数;η是可测变量的个数;是解析冗余关系矩阵中的元素。根据计 算公式，可以得到与解析冗余关系集合A= {ai，a2, . . .am}相对应的检测可靠性矢量DR = [Ri,R2 ,. ..Rm]
[0023] 故障可测性dik的计算公式如下：
[0024]
[0025] 其中，dlk为解析冗余关系&1对故障fk的可测性;TTD lk为故障fk初始时刻到解析冗 余关系ai检测到该故障的时间跨度；SyDlk为解析冗余关系&1检测到故障f k的持续时间； TTFlk为故障fk初始时刻到系统失效时的时间跨度；
=两个比值描述了解析冗余 关系ai对故障fk的故障检测时效性和故障可跟踪性;D=[djk]mxP为计算得到故障可测性矩 阵。
[0026] 3.综上所述，故障诊断系统的定量结构模型分成两部分表达，第一部分是结构模 型表达，即解析冗余关系矩阵B和故障特征矩阵M。第二部分是定量表达，由结构模型和相关 数据得到检测可靠性矢量DR和故障可测性矩阵D定量地描述故障诊断系统。最终，通过这些 定量信息计算得到故障诊断系统性能指标故障检测率FDR和故障隔离率FIR。
[0027] 步骤三基于故障诊断系统的定量结构模型，将故障检测率FDR和故障隔离率FIR作 为约束条件，经济成本作为优化目标，建立考虑检测不确定性时的传感器优化模型：
[0028]
[0029] 其中，S*为已配置传感器集合JH〇bj(S*)是以S*为变量的目标函数;Q为加 权值;W为第j个传感器的成本;FDR'，FIir是故障诊断系统必须满足的指标要求；s为最优 的传感器配置集合;该模型的意义在于将故障检测率FDR和故障隔离率FIR作为约束条件， 经济成本作为优化目标，得到约束条件下成本最低的传感器布局方案，即目标函数〇b j (S， 取最大值时的解?。
[0030] 步骤四根据遗传算法求解上述传感器优化模型，得到最优传感器配置集合，具体 过程如下：
[0031 ] 1.参数初始化，随机产生初始种群；
[0032] 2.根据传感器优化模型，建立适应度函数：
[0033]
[0034] 其中，Cost (i )为第i个传感器配置方案的成本;max (0，FDIT -FDR(i ))和max (0， Flf-FIRG))为惩罚函数，即当第i个传感器配置方案的故障检测率H)R(i)和故障隔离率 FIR( i)小于指标要求FDf和FIIT时，该方案的适应度值会减少;&、C2、C3为常数；
[0035] 3.计算每个染色体所对应的传感器集合的适应度值；
[0036] 4.根据适应度采用轮盘赌方法选择个体，并进行交叉操作，得到当前种群Pop ' ；
[0037] 对当前种群Pop'中的个体进行变异操作，得到下一代种群Pop"，判断算法是否到 达最大迭代次数，是则结束，运算结果即为最优传感器配置方案，否则返回步骤3。
[0038] 有益效果
[0039] 本发明提出的一种基于定量结构模型的故障诊断系统传感器优化配置方法，第一 步，运用结构分析方法建立故障诊断系统的结构模型。该结构模型的表现形式为两个矩阵： 解析冗余关系矩阵和故障特征矩阵。该结构模型能够定性地描述故障诊断系统。第二步，结 合故障诊断系统结构模型和传感器检测不确定性，得到故障诊断系统的定量描述，达到定 量分析故障诊断系统的目的，进而建立考虑检测不确定性的故障诊断系统关键性能指标的 计算公式。第三步建立故障诊断系统的传感器优化模型。最后一步，采用遗传算法对优化模 型进行求解。本发明提出的面向故障诊断系统的传感器优化配置方法可以有效节约故障诊 断系统的传感器成本且配置结果更加符合工程实际情况。
【附图说明】
[0040] 图1为本发明基于定量结构模型的故障诊断系统传感器优化配置方法流程图；
[0041] 图2为本发明采用的遗传算法各代故障检测率平均值；
[0042] 图3为本发明采用的遗传算法各代故障隔离率平均值；
[0043] 图4为本发明采用的遗传算法各代总成本平均值。
[0044] 图5为本发明采用的遗传算法流程图。
【具体实施方式】
[0045] 现结合实施例、附图对本发明作进一步描述：
[0046] 下面将某型号机电作动器故障诊断系统作为实例并结合附图对本发明进一步说 明。
[0047] 步骤一建立由解析冗余关系矩阵B和故障特征矩阵Μ描述的故障诊断系统的结构 模型。
[0048] 1.确定解析冗余关系集合A={ai，a2, . . .am}和系统的完备传感器(可测变量)集合 S= {si，S2, . . .，sn}，传感器信息如表1所示;确定系统中的P个故障的集合F= {fl，f2,...， fp}，故障信息如表2所示；
[0049] 2.建立解析冗余关系矩阵B=[bu]mXn，矩阵的列表示可测变量或传感器，行表示 系统的解析冗余关系，b^ = l表示解析冗余关系式&1中含有可测变量~，Μ = 0则表示解析 冗余关系式ai中不包含可测变量建立故障特征矩阵M=[mjk] mXp表示任一故障fk的故障 特征，mik= 1表示故障fk的出现会影响到解析冗余关系ai，此时ai不是满足的，mik = 0表示故 障fk的出现不会影响到解析冗余关系m，故障特征矩阵如表3所示。
[0050]步骤二求