一种混合故障预警模型的建模方法及混合故障预警模型

摘要

本发明实施例提供一种混合故障预警模型的建模方法及混合故障预警模型，该方法包括：基于HAZOP或FMEA，生成功能分析模块；基于FMEA分析结果和随机过程理论，生成退化分析模块；根据状态监测数据和维修活动信息，生成事件分析模块；根据功能分析模块和退化分析模块的输出结果，结合动态贝叶斯网络DBN理论，生成行为分析模块；以事件分析模块的输出为推理证据，利用DBN推理算法在同一时间片内进行前后向推理，生成评价模块，以输出系统故障因素和后果；以评价模块和事件分析模块的输出结果为推理证据，利用DBN推理算法在不同的时间片之间进行前后向推理，生成预测模块，输出系统各个部件未来的退化趋势。该模型可追踪系统故障因素、推理可能故障后果及概率。

说明书

技术领域

本发明涉及安全工程技术领域，具体地涉及一种混合故障预警模型的建模方法及混合 故障预警模型。

背景技术

随着状态监测技术的不断发展，越来越多的专家学者致力于研究系统故障的检测与分 类，并开发了不少成熟的状态监测与诊断软硬件系统。但“百病不如一防”，为了提高系 统的本质安全，必须通过故障预警从根本上避免故障产生的条件。现有的方法往往关注于 单个部件或一个独立子系统的退化机制和剩余寿命的研究。然而，复杂系统大多不是在实 验室那样稳定的环境中运行，而是受到各种内部和外部随机事件的影响，系统之间的交互、 环境的变化、人为失误、机器老化以及其他不确定因素都会导致故障耦合作用的发生及演 化，并使得系统的行为具有随机性。简单的将其拆分为独立的单元虽然降低了问题的难度， 却得不到最优解，也不符合实际工程情况。尤其当系统存在共因失效(common cause  failure)、相关失效或部件具有多种失效模式的情况时，现有预测模型忽略了系统变量之 间以及故障失效模式之间的相互作用及影响程度，推理与决策的合理性较低，并容易产生 误报或漏报现象。

复杂系统中各个部件故障模式自身以及相互影响的发展变化过程体现为一条故障因 果链，其通过因果关系将由最初故障模式引起的一系列故障模式关联起来。复杂系统的故 障耦合作用实际上是一个以故障因果链为基础的网络，一个或少数几个节点发生的故障 (可能是随机发生，也可能是人为因素造成的)会通过网络节点之间的耦合关系引起其他 节点发生故障，最终导致相当一部分节点甚至整个网络的崩溃。其中，故障因果链是指： 复杂系统中各个部件故障模式自身以及相互影响的发展变化过程体现为一条故障因果链， 其通过因果关系将由最初故障模式引起的一系列故障模式关联起来。

因此，为了避免故障产生的条件，提高故障根源性因素辨识的合理性，亟待研究出一 种可实现如下功能的基于故障因果链的混合故障预警模型：

(1)希望能在发生故障之前，通过混合故障预警模型进行根源性危险的辨识，采取 预控措施使系统保持在规定状态；

(2)当发生早期单点故障时，希望通过混合故障预警模型能预测已发生故障的发展 趋势以及对其他单点状态的影响，有利于采用故障隔离主动切断故障传播路径，防止故障 扩散并导致其他单点状态的破坏甚至系统的崩溃。

发明内容

本发明的目的在于，弥补国内现有的复杂系统故障预警模型的不足，提供一种准确、 合理、有效的混合故障预警模型，通过故障因果链的定量建模与推理，突破传统方法所需 的部件失效独立性等假设条件，有效的实现故障预警分析中对多部件、多危险源系统故障 根源性因素的准确辨识，以及对部件未来退化趋势及其剩余使用寿命的合理预测。

一方面，本发明实施例提供了一种混合故障预警模型的建模方法，所述方法包括：基 于危险与可操作性分析HAZOP或失效模式与效果分析FMEA，生成功能分析模块；基于 FMEA分析结果和离散时间马尔科夫随机过程理论，生成退化分析模块；根据状态监测数 据和维修活动信息，生成事件分析模块；根据所述功能分析模块和所述退化分析模块的输 出结果结合动态贝叶斯网络理论，生成行为分析模块；以事件分析模块输出的可观测状态 变量实时监测值为推理证据，利用动态贝叶斯网络推理算法在所述行为分析模块的同一时 间片内进行前后向推理，生成评价模块，输出系统故障因素和故障引发的可能后果；以所 述评价模块的输出结果和所述事件分析模块输出的可观测状态变量实时监测值为推理证 据，利用动态贝叶斯网络推理算法不同的时间片之间进行前后向推理，生成预测模块，输 出系统各个部件未来的退化趋势。

另一方面，本发明实施例提供了一种混合故障预警模型，包括：功能分析模块，用于 确定混合故障预警模型的动态节点与静态节点，动态节点对应的隐含状态变量以及隐含状 态变量的状态空间，静态节点对应的观测变量以及观测变量的状态空间，与各个节点关联 的故障原因与后果，以及故障原因与后果之间的关联关系；退化分析模块，用于根据历史 数据库中的部件失效数据与先验知识，确定混合故障预警模型各个动态节点对应的隐含状 态变量的状态转移规律和失效概率密度函数；事件分析模块，用于将监测到的实时数据和 长期失效数据存储进所述历史数据库，为评价模块和预测模块的推理过程提供推理证据； 行为分析模块，用于建立混合故障预警模型的网络结构和参数，同时利用所述历史数据库 中的状态监测数据、失效数据对所述混合故障预警模型的参数进行估计和更新；评价模块， 用于经过推理输出当前系统各个部件的隐含状态、系统故障根源性因素、各级危险原因、 危险后果以及相应的安全措施、检修计划中的至少一个；预测模块，用于根据所述各个部 件的隐含状态输出系统各个观测参数未来的变量值、部件未来的退化趋势、剩余使用寿命、 预测维修策略中的至少一个。

本发明实施例提供的上述技术方案，通过构建一系列“功能分析模块”、“退化分析 模块”、“行为分析模块”、“事件分析模块”并将其融合起来构建混合故障预警模型， 对设备系统故障因果链进行定量化建模。构建好的混合故障预警模型可以根据“事件分析 模块”实时存储的观测变量值，在同一时间片内(“空间维度”)以及不同的时间片之间 (“时间维度”)进行前后向推理，实现对系统故障根源性因素的诊断、故障可能后果的 预测、系统退化趋势的预测、以及预防性维修方案的制定等，保障系统的本质安全。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技 术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图仅仅是本发 明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可 以根据这些附图获得其他的附图。

图1为本发明实施例的混合故障预警模型构建及应用架构图；

图2为本发明实施例的混合预警模型的建模流程图；

图3a为本发明实施例的青海湖东压气站燃压机组系统整体示意图；

图3b为本发明实施例的青海湖东压气站燃压机组系统的轴流式压气机；

图3c为本发明实施例的青海湖东压气站燃压机组系统的燃气轮机燃烧室；

图3d为本发明实施例的青海湖东压气站燃压机组系统的燃气轮机辅助系统；

图4为本发明实施例的燃压机组“事件分析模块”应用实例的界面图；

图5为本发明实施例的滑油系统的混合故障预警模型网络结构示意图；

图6为本发明实施例的发生退化的滑油系统未来的可靠性发展趋势示意图；

图7为本发明实施例的在退化的滑油系统的影响下轴承的退化趋势示意图；

图8为本发明实施例的维修后滑油系统未来可靠性发展趋势示意图；

图9a为本发明实施例的空气系统退化和隐含状态变量未来变化趋势示意图；

图9b为本发明实施例的空气系统退化和观测参数变量未来变化趋势示意图；

图10a为本发明实施例的给定预防性维修方案下的空气系统退化和隐含状态变量未来 变化趋势示意图；

图10b为本发明实施例的给定预防性维修方案下的空气系统退化和观测参数变量未来 变化趋势示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的 附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本 发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例以动态贝叶斯网络为理论基础，提供一种混合故障预警模型及其建模方 法，对故障因果链进行定量化建模和推理。从多因多果的因果链中找出产生故障征兆的最 初始原因；并进一步寻找故障传播路径，预测已发生故障的未来发展趋势以及对其他单点 状态的负面影响，从而增强对故障根源性因素辨识的合理性，避免故障产生的条件，实现 故障预警的根本目标。本发明实施例在融合了系统结构/硬件、过程功能、历史的失效数 据和状态监测数据的基础上，执行系统及部件自身的内在演变过程和彼此之间交互过程的 抽象与建模。其包含6个建模过程，图1为本发明实施例的混合故障预警模型建模与应用 架构示意图，图2为本发明实施例的混合故障预警模型的建模流程图。结合参阅图1和图 2，具体实现过程包括如下步骤：

步骤100、基于HAZOP(Hazard and Operability，危险与可操作性)分析与FMEA分 析(Failure Mode and Effect Analysis，失效模式和效果分析)生成“功能分析模块”。

如图1所示，本发明的一个实施例在功能分析模块的构建过程中，需要P&ID图(工 艺和仪表流程图，Process & Instrument Diagram)和专家知识。HAZOP分析是一种现成并 已在石油化工安全分析领域普遍应用的一种方法。HAZOP分析采用标准化“引导词”对 系统中间变量设定“偏离”，沿“偏离”在系统中反向查找非正常“原因”，以及沿“偏 离”在系统中正向查找不利“后果”。查找系统中针对各重要“原因-后果”对偶已有的安 全措施，提出整改或新的安全措施建议。其中，“原因-后果”对偶是指引发危险的原因及 危险导致的后果。本实施例中，装置系统的生产过程是通过单元操作和物理化学反应使物 料发生预定物理和化学变化的过程。为了防止发生反应失控和生产事故，这些变化(多以 工艺参数表示)必须被严格控制，不能偏离预定的范围和形式，将这类工艺参数偏离预定 的范围和形式称为“偏离”(包括状态量偏离、观测量偏离)。

FMEA分析是指：根据系统可以划分为子系统、设备和元件的特点，按实际需要将系 统进行分割，然后分析各自可能发生的故障类型及其产生的影响，以便采取相应的对策。

进一步地，关于HAZOP和FMEA分析，本领域技术人员可参考下列文献(其中有详 细的可能的危险原因与后果的分析步骤)：[1]邵辉，《系统安全工程》，石油工程出版社， 2008.5。

其中，上述基于HAZOP和FMEA分析而生成的“功能分析模块”，具体可以用于： 确定混合故障预警模型中的关键参数，包括混合故障预警模型中的动态节点和静态节点； 以及，确定各个节点对应的变量及其参数空间，例如动态节点对应的隐含状态变量及其状 态空间，静态节点对应的观测变量及其状态空间；以及，各个部件之间及各部件与环境之 间的交互影响关系。其中，状态空间是指所有可能状态的集合。参数空间包括状态空间及 各状态发生的概率。

步骤102、基于FMEA的分析结果和离散时间马尔可夫随机过程理论，生成“退化分 析模块”。

本实施例中，该退化分析模块，可以用于对FMEA分析中的各个失效模式定量化地抽 象为混合故障预警模型的动态节点对应的隐含状态变量，并根据失效数据与先验知识，确 定各个隐含状态变量的状态转移规律、先验失效概率等模型参数。

其中，上述离散时间马尔可夫随机过程理论本领域技术人员可参阅下列参考文献：周 萌清，《随机过程理论》，电子工业出版社，2009.4。

步骤104、根据状态监测数据与维修活动信息生成“事件分析模块”。

具体地，上述状态监测数据可以包括从现场传感器、PLC(Programmable logic  Controller，可编程控制器)或SCADA(Supervisory Control And Data Acquisition，数据采 集与监视控制)系统获得的实时监测数据以及各种操作事件记录，这些信息被存储进历史 数据库。上述维修活动信息可以包括维修记录或者各种预防性维修计划。

本实施例中，该事件分析模块可将监测到的数据存储进历史数据库，以便为后述的行 为分析模块所执行的参数估计过程提供不断更新与完善的数据基础，该过程是保障故障预 警分析准确、可靠的基础。

步骤106、根据上述功能分析模块和退化分析模块的输出结果，并基于动态贝叶斯网 络(Dynamic Bayesian Networks，DBN)生成“行为分析模块”。

具体地，步骤106的具体处理过程可以包括：将“功能分析模块”与“退化分析模块” 输出的混合故障预警模型的隐含状态变量、观测参数变量、上述各变量自身内在的随时间 演变过程以及变量之间的交互过程融合起来(具体融合过程后文有详述)构成“行为分析 模块”；再利用历史数据库中的数据(例如失效数据、历史观测参数变量值)，采用DBN 参数估计算法对“行为分析模块”的参数进行估计。

其中，上述DBN参数估计算法本领域技术人员可参阅下列参考文献：肖秦琨，高嵩， 高晓光，《动态贝叶斯网络推理学习理论及应用》，北京：国防工业出版社，2007。

步骤108、以上述“事件分析模块”输出的可观测状态变量实时监测值为推理证据， 利用动态贝叶斯网络推理算法在行为分析模块的同一时间片内进行前后向推理，生成“评 价模块”，以输出系统故障根源性因素和故障引发的可能后果。

在本实施例中，混合故障预警模型的“评价模块”利用动态贝叶斯网络推理算法在行 为分析模块的同一时间片内进行前后向推理(Forward and backward reasoning algorithm)， 以“事件分析模块”输出的可观测状态变量实时监测值为推理证据，从“空间维度”上推 导混合故障预警模型各个动态节点对应的隐含状态变量以及相应的发生概率，从而向终端 输出系统当前各个部件的隐含状态、故障根源性因素、各级危险原因、可能的危险后果、 安全措施与检修计划中的至少一个。

步骤110、以上述评价模块的输出结果和上述事件分析模块输出的可观测状态变量实 时监测值为推理证据，利用动态贝叶斯网络推理算法在不同的时间片之间进行前后向推 理，生成预测模块，以输出系统各个部件未来的退化趋势。

在本实施例中，混合故障预警模型的“预测模块”结合“评价模块”输出的结果(如， 隐含状态变量)以及“事件分析模块”输出的可观测状态变量实时监测值为推理证据，利 用动态贝叶斯网络进行时间片段延拓(“时间维度”)，在混合故障预警模型的行为分析 模块的不同的时间片之间进行前后向推理，预测得到部件退化直至失效的时间演化过程。 从而向终端输出系统未来的观测参数变量值、隐含状态变量值、未来的可靠性变化趋势等， 为制定准确、合理的预防维修最佳时间或间隔周期提供依据。

其中，上述时间片段延拓相关概念本领域技术人员可参阅下列参考文献：肖秦琨，高 嵩，高晓光，《动态贝叶斯网络推理学习理论及应用》，北京：国防工业出版社，2007。

下面对本发明实施例图2所示方法作进一步详细说明。

(1)混合故障预警模型中“功能分析模块”的建模过程：

“功能分析模块”基于HAZOP与FMEA分析，确定混合故障预警模型中的关键参数， 包括混合故障预警模型中的动态节点，静态节点，各个节点对应的变量及其参数空间，以 及各个部件之间及其与环境之间的交互影响关系。

“功能分析模块”的主要目的是对系统故障之间的定性因果关系进行研究分析，建立 系统部件之间的相互影响关系模型(交互过程)。较佳地，本发明实施例采用一种双向功 能-失效分析机制(dual functioning-malfunctioning reasoning)。

(2)混合故障预警模型中“退化分析模块”的建模过程：

“退化分析模块”确定混合故障预警模型中动态节点对应的隐含状态变量及其状态转 移规律、先验失效概率等模型参数。

首先从“功能分析模块”的输出中选择关键的失效模式建立退化过程，然后辨识退化 过程之间的相关性。其中关键的失效模式，通常意义上指对部件安全有严重影响的失效模 式。可由离散时间马尔可夫随机过程(DTMP)来对系统关键的退化过程进行建模，具体 算法如下：

一个独立的退化过程{X_k}可由其离散状态空间χ^X和对应的转移矩阵P_X来表示。对于 随机相关(非独立)退化过程的建模，本发明实施例采用一种简化方法，将有相互作用的 几个过程合并入一个单一模型中。例如，将两个相关的退化过程{A_t}、{B_t}合并为一个“宏 观过程”{AB_t}，这种合并方式是基于DTMP，并有如下定义：

1、合并状态空间χ^A和χ^B得到状态空间χ^AB(创建一个“宏观”失效状态A_f^B_f，并 去除集合中所有系统不可达状态)；

2、“宏观过程”参数p_ih，jl由参数p_ij和p_hl的乘积转化而来，即根据式(1)和式(2) 将过程A和过程B的状态转移进行合并。其中，p_ij为状态i向状态j转化的概率，p_hl为状 态h向状态l转化的概率，p_ih，jl为同时具有状态i和状态h向同时具有状态j和状态l转化 的概率。

$p_{\mathrm{ih},\mathrm{il}}=(1-\underset{m\ne i}{\Sigma }{p}_{\mathrm{im}})p_{\mathrm{hl}}---\left(1\right)$

$p_{\mathrm{ih},\mathrm{jh}}=(1-\underset{n\ne h}{\Sigma }{p}_{\mathrm{hn}})p_{\mathrm{ij}}---\left(2\right)$

(3)混合故障预警模型中“行为分析模块”的建模过程：

“行为分析模块”，将“功能分析模块”与“退化分析模块”的输出融合起来，并进 一步地采用DBN参数估计算法对混合故障预警模型的参数进行定量估计，从而建立以 DBN为框架的更为精确的定量的混合故障预警模型，该模型能够进行定量计算并给出评价 和预测的具体数值，而不是定性的分析和描述，具体融合机制如下：

基于动态贝叶斯网络理论，将“退化分析模块”输出的各个退化过程转换为混合故障 预警模型动态变量节点；将“功能分析模块”输出的各个传输流属性转换为混合故障预警 模型静态变量节点。同时，以动态贝叶斯网络有向图模型为框架，根据“功能分析模块” 输出的因果关系，从原因节点(父节点)出发做一条有向箭头指向其直接后果节点(子节 点)，直至将所有节点都连接在DBN网络上。即，将系统结构信息(由组件与传输流描述)、 定性因果关系以及动态退化过程综合起来，对系统故障因果链进行基于DBN的定量化建 模，生成行为分析模块这种包含动态、静态节点与各个节点之间的有向连线的网络即为混 合预警模型的结构。而各个节点对应的变量信息(如各个节点隐含状态变量及其条件概率、 先验失效概率等)以及节点之间连线的定量关系(如，一组节点的联合概率、不同时间片 段之间变量的转移概率等)定义为混合预警模型的参数。这些参数可基于大量历史数据通 过样本统计学习法进行估计。

其中，上述混合预警模型参数学习的相关算法本领域技术人员可参阅下列参考文献： Baum L E，Petrie G S，Weiss N.A maximization technique occurring in the statistical analysis of  probabilistic functions of Markov chains.The Annals of Mathematical Statistics，1970，41(1)： 164-171.

(4)混合故障预警模型中“事件分析模块”的建模过程：

生成“事件分析模块”的目的在于，提供在上述混合故障预警模型中定义的参数(隐 含状态变量和观测参数变量等)的最新数据，即对混合故障预警模型中的节点变量进行动 态更新，并为故障预警分析提供推理依据。因此，该“事件分析模块”为混合故障预警模 型提供以下两种信息：

①提供状态监测系统采集的实时观测参数变量值，即对混合故障预警模型静态节点对 应的观测参数变量的数据值(如压力值、流量值等)进行周期性更新(具体周期视状态监 测系统，如SCADA系统的数据采集频率而定)，称“监测事件分析子模块”；

②针对各个过程中的支持组件，基于各种备选维修计划(引起混合故障预警模型节点 状态改变的未来的动作)，提供混合故障预警模型动态节点对应的隐含状态变量的更新条 件与更新的数据，称“维修事件分析子模块”。

(5)混合故障预警模型“评价模块”的预警分析过程：

“评价模块”的功能是通过“事件分析模块”提供的最新证据利用DBN滤波推理机 制，确定混合故障预警模型中各个节点的变量值，尤其是隐含状态变量值，并进而在混合 故障预警模型的行为分析模块的同一时间片内(“空间维度”)进行前后向推理得到故障 的根源性因素、各级故障原因、故障的可能后果以及对应的安全控制措施和检测方案。

在整个评价模块的执行过程中不断追踪隐含状态变量的当前值，即已知观测参数变量 Y_t值(由“事件分析模块”提供)，通过下列滤波推理算法不间断估计出隐含状态变量X_t：

P(X_t|y_1:t)∝P(y_t|X_t，y_1:t-1)P(X_t|y_1:t-1)＝P(y_t|X_t)[∑x_t-1P(X_t|x_t-1)P(x_t-1|y_1:t-1)] (3)

(6)混合故障预警模型“预测模块”的预警分析过程：

“预测模块”的功能是基于混合故障预警模型“评价模块”给出的部件隐含状态变量 值、“事件分析模块”提供的可观测状态变量实时监测值，通过DBN预测推理机制，在混 合故障预警模型的行为分析模块的不同的时间片之间(“时间维度”)进行推理，得到隐含 状态变量未来(τ＜t)各种可能取值的概率，从而实现预测系统或部件未来的可靠度发展 趋势(以预测的隐含状态取值概率为指标)。

“预测模块”的预测推理机制进行一步预测计算隐含状态变量未来能取值的概率 或计算观测参数变量未来能取值的概率的公式如下：

$P\left(x_{t+1}\right|Y_0^T)=\frac{{\text{Σ}}_{x_t}P\left(x_{t+1}\right|x_t){\alpha }_t\left(x_t\right)}{{\Sigma }_{x_t}{\alpha }_t\left(x_t\right)}---\left(4\right)$

或$P\left(y_{t+1}\right|Y_0^t)=\frac{{\Sigma }_{x_{t+1}}{\alpha }_{t+1}\left(x_{t+1}\right)}{{\Sigma }_{x_t}{\alpha }_t\left(x_t\right)}---\left(5\right)$

其中α_t(i)＝P(X_t＝i|y_1:t)，表示已知隐含状态变量y在时间段[1，t]的所有取值时，观 测变量X在时刻t取值为状态i的概率；P(x_t+1|x_t)表示已知任一变量当前时刻状态时，未 来下一时刻状态发生的概率。

“预测模块”的输出结果可用于维修决策：如果对未来状态的预测值符合安全标准并 可视为能够“安全”或“成功”完成系统既定目标，则不采取任何维修或检测动作。否则 必须立即执行(基于安全考虑)维修任务，或在一个计划最优时间内(预防性维修)执行。

因此，通过一系列“功能分析模块”、“退化分析模块”、“行为分析模块”、“事 件分析模块”的执行，并将其融合起来对设备系统故障因果链进行定量化建模，从而构建 了混合故障预警模型，此为混合故障预警模型的建模阶段。在故障预警分析阶段，该模型 通过“评价模块”和“预测模块”分别在“空间维度”和“时间维度”上进行推理，实现 对系统故障根源性因素的追踪、故障可能后果的推理、系统退化趋势的预测，并为预防性 维修方案的制定提供依据，从而能够保障系统的本质安全。

以下通过实际的例子来更进一步的说明本发明实施例的上述方案。

“涩-宁-兰”输气管线湖东增压站天然气燃压机组系统(GTCS)由索拉大力神130型 燃气轮机和曼透平RV050/04型离心压缩机组成。由于其组成结构复杂、工况调节频繁、 自然环境恶劣，为了提高设备本质安全水平，避免发生非正常停机、天然气泄漏等事故， 须对燃压机组系统的故障因果链进行建模，动态评价、预测其安全状态，并对其剩余使用 寿命进行评估。

调研并采集燃压机组系统的状态监测数据、历史维修记录、运行状态日报表和相关故 障统计资料等(燃压机组系统及关键部件结构见附图3a-图3d所示)，获得充分的关于燃 压机组系统以及类似的旋转机械设备的先验知识和专家经验，对混合预警模型参数与结构 的修正提供合理的依据。其中，图3d绘示的关键部件包括：主油泵、排气管、气动马达 和进气管。

根据本发明实施例中描述的上述建模方法与详细实例，针对燃压机组系统故障因果链 建立混合故障预警模型，包括如下步骤：

第一步，基于HAZOP和FMEA分析，生成“功能分析模块”，确定燃压机组系统各 个关键子系统或部件的主要失效模式、过程参量偏差、可能的危险原因和后果等。

具体地，HAZOP和FMEA分析是系统安全工程学科中的常用方法，本领域技术人员 可参考下列文献(其中有详细的可能的危险原因与后果的分析步骤)：[1]邵辉，《系统安全 工程》，石油工程出版社，2008.5。

第二步，基于第一步的FMEA分析结果和离散时间马尔可夫随机过程理论，生成“退 化分析模块”，利用燃压机组已有的历史失效数据、故障数据、设备运转统计数据以及现 场专家经验，为各个子系统或关键部件的失效模式建立定量的退化过程，并确定其参数。

第三步，在动态贝叶斯网络框架下，生成“行为分析模块”，将“功能分析模块”与 “退化分析模块”所得结果进行融合，建立混合故障预警模型。其中将“退化分析模块” 所得各个单元状态的退化过程对应生成混合故障预警模型的动态节点及相关参数。该相关 参数包括状态空间的转移概率和先验失效概率。将“功能分析模块”所得各个单元的观测 参数变量对应生成混合故障预警模型的静态节点和相关参数。上述单元是指系统分解后的 子系统、部件等的集合。同时混合故障预警模型中父节点与子节点集合及其之间的连线显 示了混合故障预警模型的网络结构，如图5所示，即表征系统中各个单元之间状态演变和 参数变化的交互作用。在图5中，D5_1：油管退化过程(OTD)；D5_2：滑油冷却器退化 过程(OCD)；D5_3：油滤元件退化过程(OFD)；D5_4：油泵退化过程(OPD)；D5_5： 滑油加热器退化过程(OHD)；S5_1：供油压力；S5_2：供油温度：S5_3：油箱液位：S5_4： 油滤压差；S5_5：油箱温度；E1_1：环境温度。图5中白色圆圈代表动态节点；灰色圆圈 代表静态节点；其中，K代表当前时间单元，K-1代表上一时刻的时间单元。

燃压机组混合预警模型中的各个参数是在已有资料(2005-2009年度设备运行情况， 以及部分2005年之前的设备维修记录等)的基础上建立的。但由于系统运行环境的不断 变化，必需对上述模型参数提供一个实时或定期的更新机制，使得混合故障预警模型的结 果符合当前安全管理和操作的需求。因此，对燃压机组系统建立了“事件分析模块”，其 应用软件界面见图4所示(软件是基于C++Builder6.0和SQL Server 2005平台上开发研 制)。燃压机组“事件分析模块”应用软件利用OPC(OLE for Process Control，用于过程 控制的OLE)远程高速数据采集技术进行数据采集与状态监测，并将各种实时数据进行分 类(如起停机数据、工况调节数据、正常运行数据、故障状态数据等)，存储到“事件分 析模块”中与监测事件相关的安全数据库中，便于数据检索与处理。“事件分析模块”还 可以定期将设备维修记录存储到“事件分析模块”中与维修事件相关的安全数据库中，便 于定期对失效数据进行统计分析。随着“事件分析模块”不断采集并存储各种有效数据， 一方面为混合故障预警模型中的参数和网络结构的动态修正与更新提供合理的依据，使模 型更加准确并能够符合环境、工况操作调节的实际情况；另一方面为混合故障预警模型中 的“评价模块”和“预测模块”(见第四、五步)的推理过程提供软/硬证据，增强故障 根源性因素辨识以及故障趋势预测结果的合理性。

第四步，混合故障预警模型“评价模块”对系统故障根源性因素辨识过程验证。

2008年5月对青海湖东站天然气燃压机组进行实时在线状态监测，某天测得润滑油平 均温度为57℃，润滑油压力为0.195MPa，均偏离正常值范围(润滑油正常温度区间为[35， 55]℃，润滑油压力正常取值范围为[0.210，0.449]MPa)。传统HAZOP分析报告将显示表 1所示的10个可能的危险原因甚至出现冗余的分析结果，且各个危险原因不分主次并列输 出，给安全决策带来了困难。若模型节点超过10个，则传统HAZOP分析结果将十分庞大， 往往出现“组合爆炸”情形。系统安全员根据列表逐个检查，即浪费人力物力，又可能错 过了系统最佳维修时间，从而使系统运行状态恶化。而本发明实施例提供的故障混合预警 模型则能够很好的解决这一难题。

现根据本发明实施例中描述的“评价模块”推理算法，各节点状态概率计算分析过程 如下：

(1)将观测节点偏差进行模糊定量化作为混合故障预警模型中节点后验概率。其中， 关于模糊定量化的技术手段，本领域技术人员可参考如下文献：胡瑾秋，张来斌，梁伟， 王朝晖.基于模糊信息融合的燃压机组定量HAZOP分析.系统工程理论与实践，2009，29 (8)：153-159。机组运行时润滑油正常温度区间为[35，55]℃，现测得润滑油平均温度为 57℃，则观测节点S72的后验概率经计算得{0，0.1546，0.7939，0.0515}，分别对应于 可能的状态{偏低，正常，偏高，超高}。同理润滑油压力正常取值范围为[0.210，0.449]MPa， 则观测节点S71的后验概率经计算得{0.0593，0.8536，0.0871，0}，分别对应于可能的 状态{超低，偏低，正常，偏高}。其中，后验概率是根据模糊定量化处理计算得到，例 如，每个可能的状态(偏低、偏高、正常等)对应各自的模糊隶属度函数，然后将实际观 测值带入，计算各个状态的模糊隶属度，作为各个状态的后验概率。异常的观测值均发生 在滑油系统中，其混合故障预警模型的网络结构如图5所示。与系统组件相对应的动态节 点D5_1，D5_2，D5_3，D5_4，D5_5为正常状态的概率可作为该部件的可靠度估计值， 经“评价模块”推理，其分别为0.3140，0.6583，0.7244，0.5395以及0.9016。其中D5_1和 D5_4节点位于正常状态的概率值最低，表明滑油管路和油泵的可靠性较差，存在安全隐患， 需要检查或修理。

表1：混合故障预警模型中“评价模块”推理结果-可能的危险原因与安全措施

(2)由混合故障预警模型的“评价模块”在“空间维度”前后向推理估计滑油系统 存在的可能的危险原因和后果(故障根源性因素)，见表1与表2所示。其显示具有最大 发生可能性(0.8035)的危险原因为“进油节流孔径偏小”；其次为“出口水温偏高”，可 能性为0.7142；具有最大发生可能性(0.9132)的危险后果为“压力油膜难以维持，甚至 造成轴瓦烧坏、密封磨损以及叶轮破坏等事故”。设备安全维护人员根据系统危险原因表 和相应的建议安全措施，按发生可能性从大到小依次检查维修。实际检查结果与“评价模 型”推理结果相符，发现系统运行潜在的安全隐患。同时根据最大可能性的危险后果制定 安全控制措施，对轴瓦、密封的温度、位移等相关参量进行密切监控。现场维修人员通过 增加冷却水量使节点观测值落入正常值区间，避免了故障的继续恶化；并在发现异常观测 值的两天后组织燃压机组清洗作业，重点对节流圈孔板进行清洗，以及检修润滑系统油泵、 油管，清洗滤油器。

表2：混合故障预警模型“评价模块”的推理结果-可能的危险后果

与传统单一HAZOP分析等安全评价、故障诊断方法分析相比，本发明实施例提供的 混合故障预警模型有利于快速、准确的发现故障根源性因素并制定安全措施，解决了传统 HAZOP定性分析的局限性与不确定性，不会输出庞大且冗余的推理结果。同时混合故障 预警模型增强了当系统存在多个节点偏差(即多个危险源)时进行辨识的准确性和决策的 合理性。现场维修人员只需依据故障原因发生可能性大小依次检查，往往检查2-3个因素 即可定位故障根源，从而主动采取事前的维修措施将这些诱发故障的因素控制在一个合理 的强度或水平以内，以防止诱发部件进一步的故障或失效。

第五步，混合故障预警模型“预测模块”对系统部件退化趋势和剩余使用寿命预测过 程验证。

根据本发明实施例中的上述“预测模块”推理算法，从当前退化状态出发在“时间维 度”上推理预测滑油系统未来的可靠性发展趋势，如图6所示。可见当前可靠度最低的D5_1 (油管)和D5_4(油泵)随部件使用年龄的增长其退化程度将更加严重，且在30天后其 预测的可靠度将低于0.5。若以工程上常用的可靠度50％阈值为设备失效判废标准，则其 剩余使用寿命预测为30天。同时，虽然D5_3(油过滤器)当前的可靠度位于安全水平， 但其退化速率很高，若不及时清洗或更换，将会影响其他部件的正常运行。考虑到相关失 效问题，退化的滑油系统对透平系统中轴承的影响如图7所示。其表明在润滑油温偏高， 润滑油压偏低的情况下，轴承(D5_1：GT 1^#2^#3^#轴承；D5_2：GT止推轴承)的退化趋势 (图7中实线D5_1和D5_2)比正常的退化趋势(分别对应于图7中虚线D5_1’和D5_2’) 更严重一些，即滑油系统的退化状态将缩短轴承系统的使用寿命。

根据本发明实施例的混合故障预警模型给出的评价与预测结果，制定维修方案，且将 维修操作定量化为“事件分析模块”中的备选维修计划，并对混合故障预警模型中的动态 节点对应的隐含状态产生影响，具体过程如下：

(1)更换滑油管路，并将该操作视为完全维修，则维修后节点D5_1的状态更新为[1 0 0 0]；

(2)清洗滑油过滤器，并将该操作视为不完全维修，则维修后节点D5_3的状态更新 为[0.9 0.10]；

(3)维修油泵，并将该操作视为不完全维修，则维修后节点D5_4的状态更新为[0.85 0.15 0]。

图8给出了在“事件分析模块”对混合故障预警模型进行修正后，即从维修后更新的 状态出发进行预测推理未来滑油系统可靠度的发展趋势。可见在未来的60天中，各个部 件工作性能良好，且可靠度基本都保持0.5以上(根据工程经验可靠度0.5以上视为可以 接受，并处于相对安全状态)。由于油过滤器自身的退化机制以及外部环境的作用，D5_3 节点对应的部件(油过滤器)的退化速率与其他部件相比较快，因此可以将油过滤器的维 修(一般为清洗或更换)周期制定为每60-90天，从而确保系统其他部件的正常工作。

第六步，预防性维修策略的制定。

混合故障预警模型的评价与预测结果可作为维修决策环节的输入：当系统(部件)可 靠性低于某一设定值(或标准)时，将在“事件分析模块”中触发相应的维修动作，从而 对混合故障预警模型的有关动态节点状态进行更新。以燃压机组空气系统为例，假定其中 的两个部件(D1_1：空气过滤器和D1_2：通风系统)均处于全新状态，在300天内未来 的系统退化趋势如图9a所示。当t＝168(天)时，空气过滤器(D1_1)的可靠度下降到 0.5(即，空气过滤器剩余使用寿命预测为168天)，而通风系统(D1_2)工作良好，无明 显退化趋势。分析代表系统运行状况的观测节点(S1_1，S1_2，S1_3，S1_4)的发展趋势 (其中S1_1为进口过滤器压差、S1_2为箱体压力、S1_3为箱体温度、S1_4为：燃气轮机入 口空气温度)，如图9b所示，S1_1(进口过滤器压差)将逐渐偏离正常值，并将对其相关 联的其他部件产生负面影响(如压气机系统、燃烧系统等)。

在这种情况下，空气过滤器的最大预防性维修周期可设定为168(天)(对应于50％安 全可靠度标准)。若为了确保充足的系统安全余量，可将维修周期设定为60-90(天)(当 可靠度低于0.7)。图10a显示了在60天维修周期条件下的空气系统退化和运行指标发展趋 势。其中维修动作假定为对空气过滤器的清洗，并视为不完全维修(通常情况下由于部件 材料的退化并存在一定的磨损与腐蚀，部件清洗后无法达到全新状态)，则维修后空气过 滤器D1_1部件状态更新为[0.9 0.10]。

如图10b所示，在维修计划的参与下，模型预测出空气系统整体的可靠度均在0.7以 上，且观测节点处于正常状态的概率也均为0.7以上，其表明该维修计划的制定使空气系 统能够长期保持良好的安全工作状态。同时，该预测结果也显示了虽然考虑了对空气系统 的定期预防性维修，但观测变量S1_2仍存在一定的退化趋势，其原因在于S1_2同时受空 气过滤器和通风系统的影响，而当时间跨度较长时，通风系统也存在缓慢的退化并影响节 点S1_2使其逐渐偏离正常值。

本实施例通过对燃压机组系统的现场应用实例，阐述了本发明的混合故障预警模型的 建立及模型参数的确定全过程，并根据现场采集数据对模型预警分析的正确性、有效性与 合理性给予了充分的验证。

本发明实施例的有益效果在于：

1、为了增强故障根源性因素辨识的合理性，本发明实施例中的混合故障预警模型在 处理复杂系统状态演变不确定性问题时具有以下突出的特点：

(1)能够处理不确定性和概率性的故障事件；

(2)能够用于学习故障因果关系；

(3)是一种将(专家)先验知识、实测数据以及物理机制进行综合的较佳表达模式；

(4)能够处理不完全信息(或部分数据丢失)的数据集；

(5)灵活性较好，能继续用专家知识和数据改进模型的性能，且模型结构和参数具 有明确的含义。

2、通过故障因果链的定量建模与推理，突破了传统方法所需的部件失效独立性等假 设条件，有效的实现了故障预警分析中对多部件、多危险源系统故障根源性因素的准确辨 识，以及故障预警控制中对剩余使用寿命的合理预测：

(1)本发明混合故障预警模型建模方法中“功能分析模块”通过功能分析和失效分 析建立了系统部件之间以及故障之间的相互影响关系模型(交互过程)，解决了如何定性 研究设备系统故障因果关系的难题。

(2)“退化分析模块”通过对所有影响系统支持组件的退化物理机制进行统计建模， 解决了如何表示故障预警模型中故障因果链节点的难题。

(3)“行为分析模块”将系统“功能分析模块”与“退化分析模块”所得结果进行 综合，以动态贝叶斯网络为系统故障因果链的有向图表达，解决了如何定量化构建故障因 果链的难题。

(4)“事件分析模块”通过实时跟踪监测事件和维修事件，提出了一种对混合故障 预警模型(结构与参数)进行动态修正与更新的机制，解决了由于外部环境的变化、工况 操作及维护活动的调节导致预警模型不适应的难题。

(5)混合故障预警模型中的“评价模块”通过“事件分析模块”提供的证据，利用 DBN滤波推理机制，从“空间维度”推理得出混合故障预警模型中各个动态节点对应的隐 含状态的概率，解决了如何基于故障因果链推理辨识故障的根源性因素、故障的可能后果 以及对应的预控措施难题，避免了故障产生的条件。

(6)混合故障预警模型中的“预测模块”通过DBN预测推理机制，从“时间维度” 预测系统或部件未来的退化演变过程，解决了如何基于故障因果链推理预测故障的发展趋 势难题。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任 何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都 应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为 准。