尧图网页设计 
简介《电力系统规划与可靠性5 可修复系统的可靠性(马氏FD串并)》是一份面向电力系统规划、可靠性工程及相关专业学生与从业者的教学PPT重点讲解可修复系统可靠性的基本概念与分析方法。内容涵盖预防性维修与故障后维修两类策略并通过可靠度、可用度与不可用度等指标系统介绍状态空间/时间法、频率平衡法及框图法在单元件和多元件串并联系统中的应用其中明晰了MTTF、MTTR、MTBF等参数关系给出可用度AMTTF/(MTTFMTTR)、不可用度UMTTR/(MTTFMTTR)的单元件指标表达以及串联系统等效故障率与修复时间的推导配有发电机、变压器、线路等具体算例便于读者对照工程实际理解。资源包为单个PPT文件大小仅801KB已有168人学习轻量便携适合课内外反复研读。通过其中状态转移图、公式推导和示例计算读者能较快掌握串并联可修复系统的可靠性建模与计算思路也可作为课件复习、考研梳理或工程入门的参考为后续电力系统可靠性评估和检修策略制定打下基础。1. 可修复系统的可靠性为什么要用马尔可夫方程电力设备有一个电子元器件没有的特性——坏了还能修。传统可靠性数学把元件当成一次性消耗品只关心首次失效时间得到的可靠度R(t)曲线放到电网规划里基本不参与决策。规划评审实际问的是三个数长期运行时间占比多少、一年发生几次故障停运、每次停运平均持续多久。这三个数全部来自马尔可夫状态空间里的状态概率和状态转移频率。可修复系统这个前提一旦成立串并联的直觉就失灵两台变压器并联后的等效故障率并不是单台的一半而是状态占比和转移频率联合作用的结果。下面顺着马氏FD串并这条主线从两状态模型开始把四状态枚举、FD参数计算和工程上的易错点讲透。2. 马尔可夫FD法两状态模型的微分方程与指标2.1 可修复系统的马尔可夫假设与状态转移图马尔可夫模型在可靠性现场能落地依赖两条强假设第一元件寿命和修复时间都服从指数分布失效率和修复率是常数λ和μ不随设备投运时间增长而变化第二未来状态只与当前状态有关与过去几个月的故障历史无关。这两条假设让状态数有限、转移率恒定微分方程组有闭式解。工程上有人质疑指数分布过于理想但变电站年度检修和定期巡检实际上把失效率和修复率都摊平成“年度平均”指数假设造成的偏差在规划精度范围内完全可以接受。拿一个最典型的可修复元件举例状态转移图只有两个圈状态0表示运行状态1表示故障。0→1的转移率是λ1→0的转移率是μ。这里的μ不是检修效率的笼统说法而等于1/MTTR。比如一台变压器平均故障修复时间MTTR为365小时则μ8760/36524次/年若失效率λ为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就是稳态可用度AP1是稳态不可用度U。这个公式很便宜却经常被用错。常见误区是把λ和μ的单位混在一起或者只用λ/μ来估计U而忽略归一化条件。当λ远小于μ时U≈λ/μ是有效近似当λ和μ同量级时必须用完整公式。可以顺手验证一下物理意义λ0.5次/年、μ24次/年时A24/24.5≈0.9796U≈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.0204f≈0.4898次/年T_d≈0.0417年换算回来恰好等于365小时和MTTR一致。这说明FD法把可用度、频率、持续时间三个指标封闭地串在一起知道任意两个就能推出第三个。指标符号含义与λ/μ的关系可用度A长期运行时间占比μ/(λμ)不可用度U长期停运时间占比λ/(λμ)故障频率f每年进入故障态的次数A·λ U·μ平均持续时间T_d每次故障的平均小时数U/f多状态系统里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λB0.1次/年、μAμB1次/年为例四个稳态概率分别为P(1,1) 0.8264P(1,0) P(0,1) 0.0826P(0,0) 0.0083串联系统只有(1,1)一个正常状态所以系统可用度A_sys0.8264不可用度U_sys0.1736。从正常状态集到故障状态集的单向转移频率为f_sys P(1,1)·(λAλB) 0.8264×0.2 0.1653次/年于是系统等效失效率λ_sysf_sys/A_sys0.2次/年正好等于λAλB。这个结论有一般性串联系统的等效失效率等于各元件失效率之和前提是元件失效率远小于修复率忽略两个元件同时故障的联合项。等效修复率μ_sysf_sys/U_sys0.9524次/年对应的MTTR≈1.05年。这个MTTR值得多说一句它比单个元件的MTTR1年还长。原因是系统停运状态里包含了(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.10.0826×0.1 0.0165次/年等效失效率λ_sysf_sys/A_sys≈0.0167次/年等效修复率μ_sysf_sys/U_sys2次/年。并联系统的MTTF从串联的5年拉长到60年MTTR则从1.05年缩短到0.5年。这个0.5年也很有规律它就是λ/(2μ)推导出来的组合修复率等于两个元件修复率的和μAμB。在双元件并联、多修理工的假设下“等效修复率相加”是马尔可夫模型给出的漂亮对偶结论和串联时“失效率相加”形成镜像。3.3 精确解与工程近似公式的边界工程手册里经常出现两个快速估算公式串联系统U_sys≈U1U2并联系统U_sys≈U1×U2。拿上面的数字验证串联近似给0.1818马尔可夫精确解是0.1736并联近似给0.00826精确解是0.00826。串联近似偏大6%并联近似几乎完全重合。系统参数串联精确解串联近似并联精确解并联近似不可用度U0.17360.18180.008260.00826等效λ(次/年)0.20000.20000.01670.0167等效μ(次/年)0.9524—2.0000—MTTF(年)5.05.060.060.0MTTR(年)1.05—0.50—近似公式在λ/μ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。直接解齐次方程会有无穷解需要把最后一个方程替换成归一化条件ΣP1变成非齐次线性方程组求解。# 稳态方程: 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(fMTTF {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)两项。稳态下这个方向的总频率与反方向相等所以取单向即可。跑出来的结果与第三节手动计算完全一致串联MTTF5年、MTTR1.05年并联MTTF60年、MTTR0.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这一行补上重新跑一遍求解脚本才能拿到运维评审认可的停运概率。本文还有配套的精品资源点击获取 可修复系统可靠性分析：马尔可夫状态空间建模与FD参数计算
简介《电力系统规划与可靠性5 可修复系统的可靠性(马氏FD串并)》是一份面向电力系统规划、可靠性工程及相关专业学生与从业者的教学PPT重点讲解可修复系统可靠性的基本概念与分析方法。内容涵盖预防性维修与故障后维修两类策略并通过可靠度、可用度与不可用度等指标系统介绍状态空间/时间法、频率平衡法及框图法在单元件和多元件串并联系统中的应用其中明晰了MTTF、MTTR、MTBF等参数关系给出可用度AMTTF/(MTTFMTTR)、不可用度UMTTR/(MTTFMTTR)的单元件指标表达以及串联系统等效故障率与修复时间的推导配有发电机、变压器、线路等具体算例便于读者对照工程实际理解。资源包为单个PPT文件大小仅801KB已有168人学习轻量便携适合课内外反复研读。通过其中状态转移图、公式推导和示例计算读者能较快掌握串并联可修复系统的可靠性建模与计算思路也可作为课件复习、考研梳理或工程入门的参考为后续电力系统可靠性评估和检修策略制定打下基础。1. 可修复系统的可靠性为什么要用马尔可夫方程电力设备有一个电子元器件没有的特性——坏了还能修。传统可靠性数学把元件当成一次性消耗品只关心首次失效时间得到的可靠度R(t)曲线放到电网规划里基本不参与决策。规划评审实际问的是三个数长期运行时间占比多少、一年发生几次故障停运、每次停运平均持续多久。这三个数全部来自马尔可夫状态空间里的状态概率和状态转移频率。可修复系统这个前提一旦成立串并联的直觉就失灵两台变压器并联后的等效故障率并不是单台的一半而是状态占比和转移频率联合作用的结果。下面顺着马氏FD串并这条主线从两状态模型开始把四状态枚举、FD参数计算和工程上的易错点讲透。2. 马尔可夫FD法两状态模型的微分方程与指标2.1 可修复系统的马尔可夫假设与状态转移图马尔可夫模型在可靠性现场能落地依赖两条强假设第一元件寿命和修复时间都服从指数分布失效率和修复率是常数λ和μ不随设备投运时间增长而变化第二未来状态只与当前状态有关与过去几个月的故障历史无关。这两条假设让状态数有限、转移率恒定微分方程组有闭式解。工程上有人质疑指数分布过于理想但变电站年度检修和定期巡检实际上把失效率和修复率都摊平成“年度平均”指数假设造成的偏差在规划精度范围内完全可以接受。拿一个最典型的可修复元件举例状态转移图只有两个圈状态0表示运行状态1表示故障。0→1的转移率是λ1→0的转移率是μ。这里的μ不是检修效率的笼统说法而等于1/MTTR。比如一台变压器平均故障修复时间MTTR为365小时则μ8760/36524次/年若失效率λ为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就是稳态可用度AP1是稳态不可用度U。这个公式很便宜却经常被用错。常见误区是把λ和μ的单位混在一起或者只用λ/μ来估计U而忽略归一化条件。当λ远小于μ时U≈λ/μ是有效近似当λ和μ同量级时必须用完整公式。可以顺手验证一下物理意义λ0.5次/年、μ24次/年时A24/24.5≈0.9796U≈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.0204f≈0.4898次/年T_d≈0.0417年换算回来恰好等于365小时和MTTR一致。这说明FD法把可用度、频率、持续时间三个指标封闭地串在一起知道任意两个就能推出第三个。指标符号含义与λ/μ的关系可用度A长期运行时间占比μ/(λμ)不可用度U长期停运时间占比λ/(λμ)故障频率f每年进入故障态的次数A·λ U·μ平均持续时间T_d每次故障的平均小时数U/f多状态系统里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λB0.1次/年、μAμB1次/年为例四个稳态概率分别为P(1,1) 0.8264P(1,0) P(0,1) 0.0826P(0,0) 0.0083串联系统只有(1,1)一个正常状态所以系统可用度A_sys0.8264不可用度U_sys0.1736。从正常状态集到故障状态集的单向转移频率为f_sys P(1,1)·(λAλB) 0.8264×0.2 0.1653次/年于是系统等效失效率λ_sysf_sys/A_sys0.2次/年正好等于λAλB。这个结论有一般性串联系统的等效失效率等于各元件失效率之和前提是元件失效率远小于修复率忽略两个元件同时故障的联合项。等效修复率μ_sysf_sys/U_sys0.9524次/年对应的MTTR≈1.05年。这个MTTR值得多说一句它比单个元件的MTTR1年还长。原因是系统停运状态里包含了(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.10.0826×0.1 0.0165次/年等效失效率λ_sysf_sys/A_sys≈0.0167次/年等效修复率μ_sysf_sys/U_sys2次/年。并联系统的MTTF从串联的5年拉长到60年MTTR则从1.05年缩短到0.5年。这个0.5年也很有规律它就是λ/(2μ)推导出来的组合修复率等于两个元件修复率的和μAμB。在双元件并联、多修理工的假设下“等效修复率相加”是马尔可夫模型给出的漂亮对偶结论和串联时“失效率相加”形成镜像。3.3 精确解与工程近似公式的边界工程手册里经常出现两个快速估算公式串联系统U_sys≈U1U2并联系统U_sys≈U1×U2。拿上面的数字验证串联近似给0.1818马尔可夫精确解是0.1736并联近似给0.00826精确解是0.00826。串联近似偏大6%并联近似几乎完全重合。系统参数串联精确解串联近似并联精确解并联近似不可用度U0.17360.18180.008260.00826等效λ(次/年)0.20000.20000.01670.0167等效μ(次/年)0.9524—2.0000—MTTF(年)5.05.060.060.0MTTR(年)1.05—0.50—近似公式在λ/μ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。直接解齐次方程会有无穷解需要把最后一个方程替换成归一化条件ΣP1变成非齐次线性方程组求解。# 稳态方程: 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(fMTTF {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)两项。稳态下这个方向的总频率与反方向相等所以取单向即可。跑出来的结果与第三节手动计算完全一致串联MTTF5年、MTTR1.05年并联MTTF60年、MTTR0.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这一行补上重新跑一遍求解脚本才能拿到运维评审认可的停运概率。本文还有配套的精品资源点击获取