可修复系统可靠性分析：马尔可夫状态空间建模与FD参数计算

简介：《电力系统规划与可靠性：5 可修复系统的可靠性(马氏FD串并)》是一份面向电力系统规划、可靠性工程及相关专业学生与从业者的教学PPT，重点讲解可修复系统可靠性的基本概念与分析方法。内容涵盖预防性维修与故障后维修两类策略，并通过可靠度、可用度与不可用度等指标，系统介绍状态空间/时间法、频率平衡法及框图法在单元件和多元件串并联系统中的应用；其中明晰了MTTF、MTTR、MTBF等参数关系，给出可用度A=MTTF/(MTTF+MTTR)、不可用度U=MTTR/(MTTF+MTTR)的单元件指标表达，以及串联系统等效故障率与修复时间的推导，配有发电机、变压器、线路等具体算例，便于读者对照工程实际理解。资源包为单个PPT文件，大小仅801KB，已有168人学习，轻量便携适合课内外反复研读。通过其中状态转移图、公式推导和示例计算，读者能较快掌握串并联可修复系统的可靠性建模与计算思路，也可作为课件复习、考研梳理或工程入门的参考，为后续电力系统可靠性评估和检修策略制定打下基础。

1. 可修复系统的可靠性为什么要用马尔可夫方程

电力设备有一个电子元器件没有的特性——坏了还能修。传统可靠性数学把元件当成一次性消耗品，只关心首次失效时间，得到的可靠度R(t)曲线放到电网规划里基本不参与决策。规划评审实际问的是三个数：长期运行时间占比多少、一年发生几次故障停运、每次停运平均持续多久。这三个数全部来自马尔可夫状态空间里的状态概率和状态转移频率。可修复系统这个前提一旦成立，串并联的直觉就失灵：两台变压器并联后的等效故障率并不是单台的一半，而是状态占比和转移频率联合作用的结果。下面顺着马氏FD串并这条主线，从两状态模型开始，把四状态枚举、FD参数计算和工程上的易错点讲透。

2. 马尔可夫FD法：两状态模型的微分方程与指标

2.1 可修复系统的马尔可夫假设与状态转移图

马尔可夫模型在可靠性现场能落地，依赖两条强假设：第一，元件寿命和修复时间都服从指数分布，失效率和修复率是常数λ和μ，不随设备投运时间增长而变化；第二，未来状态只与当前状态有关，与过去几个月的故障历史无关。这两条假设让状态数有限、转移率恒定，微分方程组有闭式解。工程上有人质疑指数分布过于理想，但变电站年度检修和定期巡检实际上把失效率和修复率都摊平成“年度平均”，指数假设造成的偏差在规划精度范围内完全可以接受。

拿一个最典型的可修复元件举例，状态转移图只有两个圈：状态0表示运行，状态1表示故障。0→1的转移率是λ，1→0的转移率是μ。这里的μ不是检修效率的笼统说法，而等于1/MTTR。比如一台变压器平均故障修复时间MTTR为365小时，则μ=8760/365=24次/年；若失效率λ为0.5次/年，那么这台变压器的平均无故障工作时间MTTF就是1/λ=2年。注意λ和μ的量纲都是“次/年”，计算时不能一个用年、一个用小时。

2.2 两状态微分方程与稳态可用度公式

设P0(t)、P1(t)分别是t时刻元件处于运行态和故障态的概率，可以写出科尔莫戈罗夫微分方程：

dP0(t)/dt = -λ·P0(t) + μ·P1(t) dP1(t)/dt = λ·P0(t) - μ·P1(t)

同时满足P0(t)+P1(t)=1。稳态时令两个导数同时为0，得到：

P0 = μ/(λ+μ)，P1 = λ/(λ+μ)

P0就是稳态可用度A，P1是稳态不可用度U。这个公式很便宜，却经常被用错。常见误区是把λ和μ的单位混在一起，或者只用λ/μ来估计U而忽略归一化条件。当λ远小于μ时，U≈λ/μ是有效近似；当λ和μ同量级时，必须用完整公式。

可以顺手验证一下物理意义：λ=0.5次/年、μ=24次/年时，A=24/24.5≈0.9796，U≈0.0204。也就是说这台变压器长期看每年平均停运约0.0204×8760≈179小时。这个数字和“每年故障多少次”是完全不同的两个维度，但规划评审两样都要。

2.3 频率-持续时间（FD）法：从状态概率到工程指标

状态概率本身不回答“一年停几次”。FD法的核心操作，是把状态概率乘上离开该状态的转移率，得到状态转移频率。以两状态系统为例：

f_0→1 = P0 × λ

这就是系统从运行态进入故障态的期望频率，也就是“每年故障次数”。稳态下单方向的0→1频率和1→0频率相等，所以取哪个方向计算都行。每次故障的平均持续时间用FD法的另一个公式：

T_d = U / f

代入上面例子：U≈0.0204，f≈0.4898次/年，T_d≈0.0417年，换算回来恰好等于365小时，和MTTR一致。这说明FD法把可用度、频率、持续时间三个指标封闭地串在一起，知道任意两个就能推出第三个。

多状态系统里FD法的价值更明显：状态数量多了以后，直接看状态概率很难判断系统整体表现，但把状态按“正常/故障”分组后，用组间转移频率算出的f和T_d，就是可以直接写进可靠性评估报告的等效参数。

3. 串并联可修复系统的四状态马尔可夫建模

3.1 串联系统：状态枚举与系统等效参数

电力系统里最常见的可修复串联结构，是一台主变带一条出线的供电单元：任何一台设备故障，系统就停运。以两个元件A、B串联为例，每个元件有运行/故障两个状态，整个系统共有四种状态。设计数1表示运行、0表示故障，四种状态分别是：

• (1,1)：A运行、B运行，系统正常
• (1,0)：A运行、B故障，系统故障
• (0,1)：A故障、B运行，系统故障
• (0,0)：A故障、B故障，系统故障

设两个元件的失效率和修复率分别为λA、λB、μA、μB，按状态转移图列写稳态方程。以λA=λB=0.1次/年、μA=μB=1次/年为例，四个稳态概率分别为：

P(1,1) = 0.8264，P(1,0) = P(0,1) = 0.0826，P(0,0) = 0.0083

串联系统只有(1,1)一个正常状态，所以系统可用度A_sys=0.8264，不可用度U_sys=0.1736。从正常状态集到故障状态集的单向转移频率为：

f_sys = P(1,1)·(λA+λB) = 0.8264×0.2 = 0.1653次/年

于是系统等效失效率λ_sys=f_sys/A_sys=0.2次/年，正好等于λA+λB。这个结论有一般性：串联系统的等效失效率等于各元件失效率之和，前提是元件失效率远小于修复率，忽略两个元件同时故障的联合项。等效修复率μ_sys=f_sys/U_sys=0.9524次/年，对应的MTTR≈1.05年。

这个MTTR值得多说一句：它比单个元件的MTTR（1年）还长。原因是系统停运状态里包含了(1,0)、(0,1)、(0,0)三种情况，其中(0,0)需要把两个元件都修复才能恢复，拉高了整体修复时间。如果直接用单台MTTR代替系统MTTR做备用容量分析，会把恢复时间算短。

3.2 并联系统：状态分组与等效修复率

并联结构对应双回路供电、双主变互为备用的场景，系统停运条件是“所有元件同时故障”。同样四个状态，但正常状态集合变成了三个：(1,1)、(1,0)、(0,1)，只有(0,0)是系统故障。

沿用相同的稳态概率，并联系统可用度A_sys = P(1,1)+P(1,0)+P(0,1) = 0.9917，不可用度U_sys = P(0,0) = 0.0083。从正常状态集到故障状态集的转移频率是：

f_sys = P(1,0)·λA + P(0,1)·λB = 0.0826×0.1+0.0826×0.1 = 0.0165次/年

等效失效率λ_sys=f_sys/A_sys≈0.0167次/年，等效修复率μ_sys=f_sys/U_sys=2次/年。并联系统的MTTF从串联的5年拉长到60年，MTTR则从1.05年缩短到0.5年。这个0.5年也很有规律：它就是λ/(2μ)推导出来的组合修复率，等于两个元件修复率的和μA+μB。在双元件并联、多修理工的假设下，“等效修复率相加”是马尔可夫模型给出的漂亮对偶结论，和串联时“失效率相加”形成镜像。

3.3 精确解与工程近似公式的边界

工程手册里经常出现两个快速估算公式：串联系统U_sys≈U1+U2，并联系统U_sys≈U1×U2。拿上面的数字验证：串联近似给0.1818，马尔可夫精确解是0.1736；并联近似给0.00826，精确解是0.00826。串联近似偏大6%，并联近似几乎完全重合。

近似公式在λ/μ<0.1时误差可以接受，但当元件修复率很低、λ/μ大于0.3时，并联近似的误差会显著扩大。原因是U1×U2忽略了“两个元件同时故障期间修复过程叠加”的效应，而马尔可夫精确解天然包含这个耦合。所以做可修复系统评估时，手算用近似找量级，正式计算建议直接走状态空间法。

4. 用Python把马氏FD参数一次算完

4.1 转移率矩阵Q的构建规则

手工列四状态方程还算轻松，元件增加到5个以上时状态数按2的n次方爆炸，必须用程序化方式构建状态转移矩阵。核心思路：先枚举全部状态，再对每个状态逐项判断每个元件能发生什么转移，把转移率填入矩阵Q。

import numpy as np # 状态编码: 1=运行 0=故障，状态(元件A, 元件B) states = [(1,1),(1,0),(0,1),(0,0)] idx = {s:i for i,s in enumerate(states)} n = len(states) lam = 0.1 # 单台失效率(次/年) mu = 1.0 # 单台修复率(次/年) mode = "parallel" # 改成 series 可切换到串联系统 def sys_ok(s, mode): if mode == "series": return s == (1,1) # 串联必须双元件运行 return s != (0,0) # 并联允许单台运行 # 构建转移率矩阵 Q，行和为零 Q = np.zeros((n, n)) for s in states: i = idx[s] a_ok, b_ok = s if a_ok: # A 从运行转移为故障 j = idx[(0, b_ok)] Q[i, j] += lam else: # A 从故障修复为运行 j = idx[(1, b_ok)] Q[i, j] += mu if b_ok: j = idx[(a_ok, 0)] Q[i, j] += lam else: j = idx[(a_ok, 1)] Q[i, j] += mu Q[i, i] = -Q[i, :].sum() # 对角线保证行和为零

这段代码的转移方向是“行索引为当前状态，列索引为目标状态”。A元件从运行转故障时，B元件状态不变，所以目标状态是(0,b_ok)；A元件从故障修好转运行时，目标状态是(1,b_ok)。每个状态的λ和μ只填一次，Q矩阵自然满足行和为零，这是后面解稳态方程的前提。

4.2 稳态求解与串并联指标输出

稳态方程是dP/dt = Q^T·P = 0。直接解齐次方程会有无穷解，需要把最后一个方程替换成归一化条件ΣP=1，变成非齐次线性方程组求解。

# 稳态方程: Q^T @ P = 0，替换最后一行做归一化 A = Q.T.copy() A[-1, :] = 1.0 b = np.zeros(n) b[-1] = 1.0 P = np.linalg.solve(A, b) # 按系统正常/故障状态分组 n_idx = [idx[s] for s in states if sys_ok(s, mode)] f_idx = [idx[s] for s in states if not sys_ok(s, mode)] A_sys = P[n_idx].sum() U_sys = P[f_idx].sum() f_sys = sum(P[i] * Q[i, j] for i in n_idx for j in f_idx) lam_sys = f_sys / A_sys mu_sys = f_sys / U_sys print(f"系统可用度 A = {A_sys:.6f}") print(f"系统不可用度 U = {U_sys:.6f}") print(f"系统故障频率 f = {f_sys:.6f} 次/年") print(f"等效失效率 lambda = {lam_sys:.6f} 次/年") print(f"等效修复率 mu = {mu_sys:.6f} 次/年") print(f"MTTF = {1/lam_sys:.3f} 年, MTTR = {1/mu_sys:.4f} 年")

f_sys的计算方式是累加所有“从正常状态集到故障状态集”的单向转移频率。比如并联模式下，只有(1,0)→(0,0)和(0,1)→(0,0)两项；串联模式下则有(1,1)→(1,0)、(1,1)→(0,1)两项。稳态下这个方向的总频率与反方向相等，所以取单向即可。跑出来的结果与第三节手动计算完全一致：串联MTTF=5年、MTTR=1.05年，并联MTTF=60年、MTTR=0.5年。

4.3 参数改动后的工程解读

把lam改成0.5次/年、mu改成24次/年，模拟真实变压器参数，并联系统可用度会从0.9917变成0.9996，故障频率约0.00042次/年，约每2400年才发生一次双机同时停运。这个量级的变化说明：设备修复率μ对并联系统可用度的影响远大于失效率λ。规划里想提升供电可用度，与其压低单台故障率，不如缩短故障修复时间或者配置快速备用——这个结论不通过马氏模型很容易被直觉带偏。

修改代码时只需要改lam和mu的数值，或者把mode切换成"series"，就能在串联与并联之间快速对比。注意不要改动Q矩阵的构建逻辑。若是换成不同厂家的两台变压器，λ和μ也不再对称，可以在循环里分别用lam_a、lam_b、mu_a、mu_b，状态枚举和分组逻辑不用动。

5. 状态空间约简与三个工程易错点

5.1 同参数元件先合并状态再列方程

四状态模型手工可解，但同一个变电站里6台同型号变压器就有64个状态。常见做法是先把同参数、同工作方式的元件按“正常台数”合并：两台同参数并联变压器从四状态压缩成“2台运行、1台运行、0台运行”三状态；三机两用一备的机组组合也可以按“运行台数”建模，状态数从8降到4。合并后每个状态的含义从“具体哪台故障”变成“同时有几台故障”，转移率是单台λ的倍数和μ的倍数。这样做不丢精度，因为同参数元件之间没有可辨识性，系统指标只取决于故障台数。

5.2 修理工数量会改变并联等效修复率

第三节的并联模型默认两台故障时可以同时修，所以(0,0)→(1,0)和(0,0)→(0,1)两个方向的转移率各是μ，合计2μ。如果整个变电站只有一个检修班组，两台都故障时必须排队，此时(0,0)状态只能以μ的速率修复其中一台，等效修复率从2μ掉到μ，系统MTTR翻倍。同样参数下并联可用度会明显下降。写代码时把第二个μ改成0，对比一下就知道修理工配置在可靠性模型里有多关键。

5.3 共因失效必须单独画转移弧

双回路供电最怕同塔架设，两台主变最怕同一场洪涝。共因失效的意思是(1,1)状态直接跳到(0,0)，而不经过单台故障的中间状态。这时要在状态转移图里加一条(1,1)→(0,0)的转移弧，转移率记为λc。加入后并联系统等效故障频率从0.0165次/年直接跳升，可用度可能从0.9917掉到0.95以下。工程上忽略共因失效的并联计算，结果通常乐观一个数量级，N-1校核结论会被带偏。给模型加共因弧时，把Q矩阵里Q[idx[(1,1)], idx[(0,0)]] += λc这一行补上，重新跑一遍求解脚本，才能拿到运维评审认可的停运概率。

本文还有配套的精品资源，点击获取