4、故障树定量计算：顶事件概率计算、最小割集与最小径集求法、重要度分析

故障树分析（FTA）是安全评价中最重要的演绎分析方法之一。在完成故障树的定性分析（即找出所有最小割集）之后，定量计算成为进一步量化系统风险的关键步骤。定量计算的核心任务包括三个层面：一是计算顶事件的发生概率，以衡量系统整体的风险水平；二是通过最小割集与最小径集的求取，识别系统的薄弱环节；三是通过重要度分析，为维修策略、安全投入和系统优化提供量化依据。本章将围绕这三项核心任务，结合工程实例，系统讲解计算方法与实操要点。

4.1 顶事件概率计算：从底事件到系统风险

4.1.1 核心原理

顶事件概率计算是故障树定量分析的基础。其基本原理是：在已知各底事件发生概率的前提下，通过故障树的逻辑结构（与门、或门），利用概率论方法逐层向上推算，最终得到顶事件的发生概率。对于由与门连接的输入事件，其输出概率为各输入事件概率的乘积；对于由或门连接的输入事件，其输出概率为1减去各输入事件不发生概率的乘积。

4.1.2 计算方法与操作流程

顶事件概率计算通常采用以下三种方法：

操作流程如下：

1. 数据准备：收集各底事件的失效概率（可从设备可靠性数据库、历史故障统计或专家判断获得）。
2. 逻辑简化：检查故障树是否存在冗余逻辑，必要时先进行布尔化简。
3. 逐层计算：从最底层开始，按门类型逐层向上计算中间事件概率。
4. 结果验证：将计算结果与历史故障数据或行业基准进行对比，验证合理性。

4.1.3 案例拆解：化工储罐超压保护系统

某化工储罐设有压力安全阀（A）和紧急泄压系统（B），顶事件为“储罐超压破裂”。故障树逻辑为：顶事件 = A失效 ∩ B失效（即两个保护同时失效）。已知：P(A失效)=0.01（安全阀拒动概率），P(B失效)=0.02（泄压系统故障概率）。

计算过程：

P(顶事件) = P(A失效) × P(B失效)
           = 0.01 × 0.02
           = 0.0002（即2×10⁻⁴）

若该储罐年运行时间为8760小时，则每年因超压破裂的风险概率约为0.0002次/年，即平均5000年发生一次。该结果可用于判断是否需要增加第三重保护（如爆破片）。

4.2 最小割集与最小径集求法

4.2.1 核心原理

最小割集是导致顶事件发生的最少底事件组合。只要最小割集中的所有底事件同时发生，顶事件必然发生。最小割集揭示了系统的“最薄弱环节”——一阶割集（单事件割集）意味着该单一事件失效即可导致系统失效，是系统设计中最需要关注的。

最小径集是使顶事件不发生的最少底事件组合。只要最小径集中的所有底事件都不发生，顶事件就不会发生。最小径集揭示了系统的“冗余路径”——径集越多，系统越安全。

4.2.2 求取方法

4.2.3 案例拆解：双发动机飞机推力丧失

某双发飞机，顶事件为“推力完全丧失”。故障树逻辑：顶事件 = 发动机1失效 ∩ 发动机2失效。每台发动机失效又分为“燃油系统故障”和“机械故障”的或门组合。

布尔表达式：

T = E1 ∩ E2
E1 = F1 ∪ M1
E2 = F2 ∪ M2

展开：T = (F1 ∪ M1) ∩ (F2 ∪ M2)
     = (F1∩F2) ∪ (F1∩M2) ∪ (M1∩F2) ∪ (M1∩M2)

最小割集：{F1,F2}、{F1,M2}、{M1,F2}、{M1,M2}，共4个二阶割集。没有一阶割集，说明系统具有冗余设计，单一部件失效不会导致顶事件发生。

最小径集：将故障树转化为成功树，求成功树的最小割集后取对偶，得到最小径集为：{F1,M1}和{F2,M2}。这意味着只要保证任一发动机的燃油系统和机械系统都正常，飞机就不会完全丧失推力。

4.3 重要度分析：量化底事件对顶事件的贡献

4.3.1 核心原理

重要度分析用于衡量各底事件（或最小割集）对顶事件发生概率的贡献程度。三种重要度各有侧重：结构重要度仅考虑故障树的结构关系，不考虑底事件概率差异；概率重要度考虑底事件概率变化对顶事件概率的影响程度；关键重要度综合了概率重要度和底事件自身概率，反映底事件概率与顶事件概率的相对变化率。

4.3.2 计算方法

4.3.3 案例拆解：化工反应釜温度失控

某反应釜温度控制系统，顶事件为“温度失控导致爆炸”。底事件：A=温度传感器故障（P=0.01），B=控制器故障（P=0.005），C=冷却阀卡滞（P=0.02）。故障树逻辑：T = A ∪ (B ∩ C)。

计算步骤：

1. 顶事件概率：P(T) = P(A) + P(B)×P(C) − P(A)×P(B)×P(C) = 0.01 + 0.005×0.02 − 0.01×0.005×0.02 ≈ 0.010099
2. 概率重要度：
   I_g(A) = ∂P(T)/∂P(A) = 1 − P(B)×P(C) = 1 − 0.0001 = 0.9999
   I_g(B) = P(C) − P(A)×P(C) = 0.02 − 0.0002 = 0.0198
   I_g(C) = P(B) − P(A)×P(B) = 0.005 − 0.00005 = 0.00495
3. 关键重要度：
   I_c(A) = (0.01/0.010099) × 0.9999 ≈ 0.9901
   I_c(B) = (0.005/0.010099) × 0.0198 ≈ 0.0098
   I_c(C) = (0.02/0.010099) × 0.00495 ≈ 0.0098

结果分析：温度传感器（A）的关键重要度高达0.99，远超其他底事件。这意味着要降低顶事件风险，最有效的措施是提高温度传感器的可靠性（如采用冗余传感器或提高检修频率），而不是优先改进冷却阀或控制器。

4.3.4 常见误区与总结复盘

• 误区一：只计算结构重要度就下结论。结构重要度未考虑底事件实际概率差异，可能导致排序失真。例如，某底事件结构上很关键，但其失效概率极低（如10⁻⁶），实际风险贡献可能很小。
• 误区二：忽略底事件之间的相关性。当两个底事件存在共因失效（如同一型号传感器、同一电源供电）时，独立假设会导致概率计算偏低。
• 误区三：将概率重要度与关键重要度混用。概率重要度反映“敏感度”，关键重要度反映“实际贡献”，两者排序可能不同，需根据决策目标选择。