简介这份PPT围绕电力系统发输电组合系统可靠性评估展开适合电力系统规划、运行及可靠性方向的高年级本科生、研究生或工程技术人员学习参考。内容以状态解析法为主线系统介绍了枚举系统状态、失效判断、失负荷量统计及可靠性指标累计的迭代流程并结合RBTS标准算例展示LOLP、LOLF、EENS等关键指标的计算过程同时指出状态解析法物理概念清晰、精度高但存在“维数灾”的局限。资源为单个PPT文件压缩包大小约2.19MB讲述逻辑由方法原理到算例验证便于课堂讲解或自学梳理。目前已有95人学习下载适合希望掌握发输电组合系统可靠性评估基本思路与指标计算方法的读者作为配套课件使用。1. 发输电组合系统可靠性评估一份能让你把LOLP指标算明白的课件发输电组合系统可靠性评估最劝退的不是公式本身而是那张动辄几十行的事件枚举表。拿到这份PPT的人多半是想把LOLP、LOLF、EENS这几个指标真正算明白的工程师或者研究生。这份资料用两个完整算例——一个5元件系统、一个含3台发电机和2回输电线的系统——把状态解析法的全部计算步骤走了一遍连增强措施前后的指标对比都直接给齐了。它适合刚开始接触可靠性评估的规划人员也适合想找一套可复现算例去校核自己程序的从业者。最关键的是PPT末尾还补了配电网可靠性的FMEA法等于一次给了两套评估思路。2. 状态解析法的骨架四步迭代、事件枚举与维数灾的边界2.1 状态解析法的四步流程从枚举到指标更新状态解析法的基本思路是把系统里每个元件“可用/不可用”的所有组合列出来逐个状态去判断系统是否失效。PPT里把这套过程拆成了四步枚举产生一个系统状态分析该状态是否为失效状态若是失效状态则记录各负荷点的失负荷量计算该失效状态的可靠性指标修改或更新累计指标。这四个步骤循环迭代直到所有枚举状态都被处理完。这个流程的物理概念很清晰。发输电组合系统和单纯发电系统可靠性评估最大的区别在于它不仅要看发电容量是否够还要看输电网络能不能把电送到负荷点。所以“分析系统状态”这一步内部其实要跑一次网络约束判断比如线路潮流是否越限、变压器是否过载然后才能得出“该状态削减了多少负荷”。PPT里算例1和算例2给出的失负荷量本质上就是这一步的输出结果。我一般会把“枚举”和“判定”拆开来看。枚举是纯组合问题N个元件就有2的N次方个状态判定则是网络分析问题可能要调潮流计算。PPT里两个算例的元件数都不多但已经把这两层逻辑展示得很完整——这也是为什么它适合作为入门资料。2.2 概率、频率与失电量的计算口径状态的概率计算核心是可用率和不可用率的连乘。某系统状态的概率P(s)等于所有可用元件的可用率A相乘再乘上所有故障元件的不可用率U。比如算例2里发电机G1、G2、G3和线路L1、L2全部正常工作的概率就是0.94805的三次方乘以0.99727的平方约等于0.8475。频率的计算稍微绕一点。PPT里给出的频率公式f(s)Σλ_k·P(s)在实际使用时要格外小心这里的λ_k应该是“该状态下所有可能转移出去的转移率之和”也就是说不仅要考虑故障元件修复回来的修复率还要考虑可用元件继续故障的故障率。如果只把元件故障率加一遍频率算出来会偏小。失电量的计算则是概率乘以失负荷量。某个失效状态对EENS的贡献是状态概率乘以该状态削减的负荷功率再乘以8760小时得到年化电量。PPT里的算例1最终给出EENS4767.36MWh/年这个数值就是把所有失负荷事件的概率、削减量与时间维度加权汇总后得到的结果。注意EENS和LOLP不是同一个维度的指标——LOLP是失负荷状态的总概率无量纲EENS是失电量有能量的量纲。2.3 “维数灾”的量化边界元件数到几就枚举不动了PPT里明确点出状态解析法存在“维数灾”问题。这很实在。5个元件只有32个状态手算表格勉强可行10个元件是1024个状态20个元件就已经是104万多个状态。如果每个状态都要跑一次网络分析计算量会迅速失控。我见过有人用状态解析法去评估一个40多个元件的区域电网枚举状态数量到了千亿级别最后只能靠截断处理和近似方法收场。实际工程中当元件数超过15到20个状态解析法通常会配合“截断”使用比如只枚举到二阶或三阶故障事件。PPT算例2里明确说“为简化只讨论二阶事件”这就是截断的思路。但截断带来的误差有多大取决于系统的冗余度——冗余越高低阶故障导致失负荷的概率越小截断就越安全。这一点第五章还会展开讲。提示判断一份可靠性评估结果是否可信第一步就是看它枚举到几阶事件以及被截断掉的概率量级有多大。高阶事件概率如果已经落到10的负6次方以下截断通常可以接受。3. 算例1手算复现RBTS五元件系统的LOLP、LOLF、EENS从哪来3.1 RBTS单线图的元件参数与负荷设置算例1用的是RBTS测试系统的一个简化局部系统里有5个元件发电机G1和G2、变压器T1和T2、线路L1系统负荷90MW。元件可靠性参数如下表元件故障率次/年修复时间小时/次容量MW发电机G12.050.0100发电机G22.050.050变压器T10.05120120变压器T20.05120120线路L10.816150注意容量这一列T1、T2、L1的“容量”在发输电组合系统里指的是输变电能力不是发电容量。所以这个系统实际上只有G1和G2两个电源G1容量100MW、G2容量50MW总装机150MW负荷90MW。从发电容量看冗余很高但变压器和线路的网络约束会把一部分可用容量卡住——这就是为什么T1故障时系统要削减23MW负荷而不是简单靠G1G2把90MW全部带起来。3.2 手动复现事件表概率、频率、失电量三步走PPT给出了这个系统所有事件的分析表从所有元件正常、一阶事件、二阶事件一直列到五阶全故障。我拿到这套表之后习惯用脚本把它重新算一遍既验证表里的数值也方便改参数做敏感性分析。下面是一个用Python枚举全部状态并汇总指标的骨架脚本。# 状态解析法枚举脚本验算PPT算例1的事件表 import itertools # 元件配置名称, 可用率, 不可用率 # 可用率和不可用率由故障率λ与修复率μ换算Aμ/(λμ)Uλ/(λμ) comps [ {name: G1, A: 0.98887223, U: 0.01112777}, {name: G2, A: 0.98887223, U: 0.01112777}, {name: T1, A: 0.99933233, U: 0.00066767}, {name: T2, A: 0.99933233, U: 0.00066767}, {name: L1, A: 0.99857565, U: 0.00142435}, ] # 失负荷量字典键为故障元件名值为该状态下的负荷削减量(MW) # 这里按PPT事件表的网络约束分析结果填入实际工程中应由潮流计算给出 loss_dict { (G1,): 40, (T1,): 23, (L1,): 49, (G1, G2): 100, # ... 其余状态按PPT表补全 } states [] for mask in itertools.product([0, 1], repeatlen(comps)): # mask0表示正常1表示故障 prob 1.0 failed [] for i, (comp, bit) in enumerate(zip(comps, mask)): if bit 1: prob * comp[U] failed.append(comp[name]) else: prob * comp[A] loss loss_dict.get(tuple(sorted(failed)), 0) if loss 0: states.append((tuple(failed), prob, loss)) # 汇总可靠性指标 lolp sum(p for _, p, _ in states) eens sum(p * loss * 8760 for _, p, loss in states) print(fLOLP {lolp:.6f}) print(fEENS {eens:.2f} MWh/年)这段脚本的枚举逻辑很简单itertools.product生成每个元件的0/1状态组合0代表可用、1代表故障每个组合的概率由可用率和不可用率连乘得出然后通过失负荷量字典判断这个状态是否失效并取对应的削减量。最后汇总LOLP和EENS。有两点要说明。第一失负荷量字典在PPT里是现成的因为网络约束分析已经做完了如果你要评估一个新系统这一列必须由潮流计算或直流潮流分析产生不能拍脑袋填。第二脚本里未列出的故障组合我认为不失负荷这对应PPT表中失负荷量为0的那些行实际使用时要逐行核对。3.3 表格里的系统可靠性指标LOLP、LOLF、EENS怎么验算PPT算例1最终给出系统可靠性指标LOLP0.013405LOLF2.896923次/年EENS4767.36MWh/年。这三个数怎么来的可以直接从事件表里验。LOLP是全部失负荷事件概率之和。把表里所有失负荷量大于0的事件概率加总就得到0.013405。这意味着该长期运行工况下系统处于失负荷状态的时间占比约为1.34%折算成时间大约是每年117个小时处于失负荷状态。LOLF是失负荷事件的发生频率单位是次/年。这个指标在状态解析法里不是简单把事件概率相加而是每个失效状态的概率乘以该状态的离开率再求和。离开率包括故障元件的修复率和正常元件的故障率所以在算例1里G1故障这个事件的频率不是直接拿G1的故障率2次/年而是要经过状态转移率的换算。EENS的验算我通常会拆成两步。先把每个失效状态的概率乘以失负荷量得到该状态的单次期望失电量再乘以8760小时换算成年化指标。把所有失效事件加总后和PPT给出的4767.36MWh/年对照误差应该在0.5%以内。如果对不上优先检查可用率是否由λ和μ正确换算其次检查失负荷量是否抄错行。提示验算时最常犯的错是把“累积概率”当成“LOLP”。PPT事件表的累积概率是按概率从大到小排序后的累加它逼近的是1.0而LOLP只是其中失负荷事件概率的加总两者含义不同。4. 算例2方案对比二阶事件枚举与增强措施的效果量化4.1 元件可用率计算从λ和μ到A值算例2的系统配置是3台发电机加2回输电线。设每台发电机和每回线路的故障率、修复率分别相等发电机λG2次/年、μG36.5次/年线路λL2次/年、μL365次/年。发电机的可用率Aμ/(λμ)36.5/38.5≈0.94805不可用率U≈0.05195线路的可用率A365/367≈0.99727不可用率U≈0.00273。这里要注意量纲。μG36.5次/年对应的平均修复时间是8760/36.5240小时也就是10天μL365次/年对应的平均修复时间是24小时。发电机的修复时间远长于线路这是符合实际的——发电机故障后检修周期以天计线路故障后巡线、抢修通常能在一天内完成。所以在同一个系统里不同元件的λ、μ必须保持同一单位再做除法否则可用率会算错一个数量级。全部元件正常工作的概率是0.94805的三次方乘以0.99727的平方约等于0.8475。这个值看起来很“安全”但真正决定可靠性指标的是那些低概率高损失的事件组合——比如两台发电机同时故障概率只有0.002545却要削减200MW负荷对EENS的贡献远大于它的概率数字。4.2 二阶事件表解读哪些组合在拖指标PPT按二阶事件列出了系统状态分析表。我挑几个关键行来看事件故障概率频率次/年负荷削减MW失电量MWhG1、G2故障0.0025450.1959382004458.227G1、L2故障0.0001270.051675400445.437G2、G3故障0.0025450.1959382004458.227G2、L2故障0.0001270.051675200222.723L1、L2故障0.0000060.00467440022.250L2单回线故障0.0023200.8630032004064.465注意L2单回线故障这行特别显眼概率只有0.00232但频率0.863003次/年几乎是全部二阶事件频率的一半失电量4064.465MWh排第二仅次于双发电机故障。这说明什么单回L2故障虽然只削减200MW不算最严重但它发生的频率高对系统可靠性的拖累反而比那些概率更低的高阶组合更大。再看G1、L2组合故障削减400MW而G2、L2组合故障只削减200MW。原因是G1是200MW机组L2故障后再叠加G1退出系统缺口的规模更大G2故障叠加L2故障时还有G1和G3两台机组撑着削减量相对小。这提示我们评估发输电组合系统时不能只看“几阶故障”还要看具体是哪几个元件组合——发电机容量差异和线路传输能力差异会直接改变削减量。4.3 增强措施怎么量化扩容与加线的指标对比PPT给出一组增强性措施把线路L1的容量从300MW提升到500MW同时增装一回线路L3L3的故障率和修复率与L1、L2相同。增强前后的系统可靠性指标对比非常直观指标改善前改善后变化幅度LOLP0.0105960.007626下降约28%LOLE小时/年92.768466.804下降约28%LOLF次/年1.7138660.603159下降约65%EENSMWh/年18797.7213371.44下降约29%这个对比表含金量很高。LOLP和EENS都下降了近三成但LOLF下降了近三分之二原因在于增强措施让很多原本会造成削减的单重故障直接变成了“无影响事件”。改善后的事件表里L1单回线故障、L2单回线故障、G1与任意单回线的组合故障失负荷量全都归零。单回线故障不再导致停电整个系统的失负荷频率自然大幅下降。我特别想强调PPT里那句话实际系统中采用增强措施的依据要综合考虑设备投资费用、运行维修费用、线损改变费用和停电损失费用并且要用同一年值口径去比。也就是说把容量从300MW提到500MW、增加一回线路这些举措的决策是否划算不能只看可靠性指标改善了多少还得把停电损失折算成钱和投资费用放同一个时间坐标上做等年值比较。这才是规划评估落地时真正要回答的问题。5. 状态解析法避坑指南五个我差点信了公式的高频翻车点5.1 频率算出来总是对不上表现象用PPT里的频率公式fΣλ·P去算结果和表里的LOLF总差一截差的不是小数而是成倍的偏差。原因频率公式里的λ是“离开该状态的所有转移率之和”不是一个元件故障率那么简单。比如算例2里G1、G2双故障状态离开这个状态有两个途径——G1修复或者G2修复两条路径的转移率都要算进去分别是μG1和μG2。如果只把λG1和λG2加进去方向就反了数值自然对不上。解决先把每个故障元件的修复率μ按8760/修复小时数换算好状态频率状态概率×(Σ故障元件μΣ正常元件λ)再按月汇总。这个口径在工程软件里也是通用的直接背下来用。5.2 把“不可用率”当成“故障率”去算概率现象算元件故障概率时直接拿λ2次/年当作概率代入连乘得出一个比正确结果大几十倍的LOLP。原因λ的量纲是次/年它是转移率不是概率。只有当λ数值极小且时间窗取1年时λ才近似等于概率。发输电系统里发电机λ2次/年对应不可用率0.05195差了十几倍线路λ2次/年对应不可用率只有0.00273差了三个数量级完全不能混用。解决所有状态概率必须用A和U计算而A、U由λ和μ共同决定。动手前先列一张元件参数换算表把每个元件的A、U算好再进枚举脚本。我自己的习惯是固定用“Aμ/(λμ)Uλ/(λμ)”不搞近似。5.3 失负荷量用“容量差”代替“网架约束”现象评估算例1时G2故障明明只剩G1100MW高于负荷90MW但把G2故障当作无影响事件结果算出来的LOLP比PPT少了快一半。原因算例1里T1故障要削减23MW、L1故障要削减49MW这些削减不是发电容量不够造成的而是变压器、线路的传输能力约束造成的。发输电组合系统评估的失负荷量必须来自网络分析不能简单用发电容量减负荷。解决把失负荷量作为枚举环节的外部输入由潮流计算或直流潮流分析提供。PPT里的事件表已经把网络约束结果固化下来了我用脚本复现时直接查表即可新场景下必须先跑网络分析再回填。注意状态解析法的枚举可以自动化但“判定是否失效、削减多少”永远是网络分析问题。把这两层解耦程序才不会越写越乱。5.4 累积概率排序列越看越糊涂现象把事件表按累积概率排序后发现LOLP直接等于1减去某一个累积概率不对怎么减都对不上。原因事件表的累积概率是按“单个状态概率从大到小”排序后的累加和系统是否失负荷没有对应关系。全正常状态概率最大排第一二阶、三阶事件按概率大小插入其中累积概率只是用来查“概率分布尾部”的参考列不是用来反推LOLP的。解决LOLP只统计失负荷事件把表中失负荷量大于0的行挑出来概率直接相加。不要试图从累积概率列反推那列是给人看概率分布的不是给公式用的。验算时我通常让脚本重新按失负荷事件单独求和和表尾给的LOLP值对照。5.5 二阶截断之后发现高阶事件漏掉了大额失电量现象某个系统按二阶截断算完EENS看起来不高但把三阶、四阶事件补进去后EENS涨了一倍。原因虽然高阶事件概率小但如果削减量大概率×削减量的乘积未必可以忽略。算例1里四阶事件概率已经到10的负8次方量级但失负荷100MW甚至90MW贡献的失电量还是微小的吗——到10的负8次方量级×90MW确实微小但如果系统冗余低、低阶故障概率被网络约束压掉高阶组合反而成为主要的失效模式。解决截断前先做一版无截断枚举观察EENS随阶数的收敛速度如果三阶事件对EENS的贡献超过5%就得把截断阶数提高或者改用解析网络法、蒙特卡洛法。PPT算例2明确只讨论二阶事件是因为这个系统二阶事件已经覆盖了全部主要失效模式这是个案不能直接推广。6. 从发输电到配电网用FMEA法反向验证枚举结果FMEA法故障模式与影响分析是配电网可靠性评估的经典方法和状态解析法看起来完全不同骨子里其实是同一件事。状态解析法枚举的是“元件状态组合”FMEA枚举的是“元件故障事件”两者的核心动作都是“故障事件→分析影响→累计指标”。配电网通常辐射状结构网络分析比输电网简单但开关设备的存在让故障隔离范围成为关键变量。断路器能切负荷电流和短路电流故障时与保护配合快速切除隔离开关只能在无负荷电流时分合提供检修断点负荷开关介于两者之间能切额定负荷电流但不能切短路电流。这三个设备在FMEA表里决定了一个元件故障后“停电范围有多大、多久能恢复”。我在做配网可靠性评估时会把每个开关的动作逻辑画成一张影响表再和状态解析法的枚举结果交叉验证——如果某个元件的FMEA影响分析和状态枚举给出的失负荷量矛盾那一定是开关配置或者枚举判定里有一处错了。从那以后我每次拿到可靠性评估结果都强制自己走一遍“先验概率、再验频率、最后对LOLP量级”的三板斧先确认所有概率都由A、U连乘得出再确认频率用了离开率求和最后用EENS8760×Σ(概率×失负荷量)做闭环验算。三个数能互相咬合才敢把这份算例的结论拿到规划方案里去用。希望帮到你。本文还有配套的精品资源点击获取 发输电组合系统可靠性评估：状态解析法算例详解与LOLP指标计算
简介这份PPT围绕电力系统发输电组合系统可靠性评估展开适合电力系统规划、运行及可靠性方向的高年级本科生、研究生或工程技术人员学习参考。内容以状态解析法为主线系统介绍了枚举系统状态、失效判断、失负荷量统计及可靠性指标累计的迭代流程并结合RBTS标准算例展示LOLP、LOLF、EENS等关键指标的计算过程同时指出状态解析法物理概念清晰、精度高但存在“维数灾”的局限。资源为单个PPT文件压缩包大小约2.19MB讲述逻辑由方法原理到算例验证便于课堂讲解或自学梳理。目前已有95人学习下载适合希望掌握发输电组合系统可靠性评估基本思路与指标计算方法的读者作为配套课件使用。1. 发输电组合系统可靠性评估一份能让你把LOLP指标算明白的课件发输电组合系统可靠性评估最劝退的不是公式本身而是那张动辄几十行的事件枚举表。拿到这份PPT的人多半是想把LOLP、LOLF、EENS这几个指标真正算明白的工程师或者研究生。这份资料用两个完整算例——一个5元件系统、一个含3台发电机和2回输电线的系统——把状态解析法的全部计算步骤走了一遍连增强措施前后的指标对比都直接给齐了。它适合刚开始接触可靠性评估的规划人员也适合想找一套可复现算例去校核自己程序的从业者。最关键的是PPT末尾还补了配电网可靠性的FMEA法等于一次给了两套评估思路。2. 状态解析法的骨架四步迭代、事件枚举与维数灾的边界2.1 状态解析法的四步流程从枚举到指标更新状态解析法的基本思路是把系统里每个元件“可用/不可用”的所有组合列出来逐个状态去判断系统是否失效。PPT里把这套过程拆成了四步枚举产生一个系统状态分析该状态是否为失效状态若是失效状态则记录各负荷点的失负荷量计算该失效状态的可靠性指标修改或更新累计指标。这四个步骤循环迭代直到所有枚举状态都被处理完。这个流程的物理概念很清晰。发输电组合系统和单纯发电系统可靠性评估最大的区别在于它不仅要看发电容量是否够还要看输电网络能不能把电送到负荷点。所以“分析系统状态”这一步内部其实要跑一次网络约束判断比如线路潮流是否越限、变压器是否过载然后才能得出“该状态削减了多少负荷”。PPT里算例1和算例2给出的失负荷量本质上就是这一步的输出结果。我一般会把“枚举”和“判定”拆开来看。枚举是纯组合问题N个元件就有2的N次方个状态判定则是网络分析问题可能要调潮流计算。PPT里两个算例的元件数都不多但已经把这两层逻辑展示得很完整——这也是为什么它适合作为入门资料。2.2 概率、频率与失电量的计算口径状态的概率计算核心是可用率和不可用率的连乘。某系统状态的概率P(s)等于所有可用元件的可用率A相乘再乘上所有故障元件的不可用率U。比如算例2里发电机G1、G2、G3和线路L1、L2全部正常工作的概率就是0.94805的三次方乘以0.99727的平方约等于0.8475。频率的计算稍微绕一点。PPT里给出的频率公式f(s)Σλ_k·P(s)在实际使用时要格外小心这里的λ_k应该是“该状态下所有可能转移出去的转移率之和”也就是说不仅要考虑故障元件修复回来的修复率还要考虑可用元件继续故障的故障率。如果只把元件故障率加一遍频率算出来会偏小。失电量的计算则是概率乘以失负荷量。某个失效状态对EENS的贡献是状态概率乘以该状态削减的负荷功率再乘以8760小时得到年化电量。PPT里的算例1最终给出EENS4767.36MWh/年这个数值就是把所有失负荷事件的概率、削减量与时间维度加权汇总后得到的结果。注意EENS和LOLP不是同一个维度的指标——LOLP是失负荷状态的总概率无量纲EENS是失电量有能量的量纲。2.3 “维数灾”的量化边界元件数到几就枚举不动了PPT里明确点出状态解析法存在“维数灾”问题。这很实在。5个元件只有32个状态手算表格勉强可行10个元件是1024个状态20个元件就已经是104万多个状态。如果每个状态都要跑一次网络分析计算量会迅速失控。我见过有人用状态解析法去评估一个40多个元件的区域电网枚举状态数量到了千亿级别最后只能靠截断处理和近似方法收场。实际工程中当元件数超过15到20个状态解析法通常会配合“截断”使用比如只枚举到二阶或三阶故障事件。PPT算例2里明确说“为简化只讨论二阶事件”这就是截断的思路。但截断带来的误差有多大取决于系统的冗余度——冗余越高低阶故障导致失负荷的概率越小截断就越安全。这一点第五章还会展开讲。提示判断一份可靠性评估结果是否可信第一步就是看它枚举到几阶事件以及被截断掉的概率量级有多大。高阶事件概率如果已经落到10的负6次方以下截断通常可以接受。3. 算例1手算复现RBTS五元件系统的LOLP、LOLF、EENS从哪来3.1 RBTS单线图的元件参数与负荷设置算例1用的是RBTS测试系统的一个简化局部系统里有5个元件发电机G1和G2、变压器T1和T2、线路L1系统负荷90MW。元件可靠性参数如下表元件故障率次/年修复时间小时/次容量MW发电机G12.050.0100发电机G22.050.050变压器T10.05120120变压器T20.05120120线路L10.816150注意容量这一列T1、T2、L1的“容量”在发输电组合系统里指的是输变电能力不是发电容量。所以这个系统实际上只有G1和G2两个电源G1容量100MW、G2容量50MW总装机150MW负荷90MW。从发电容量看冗余很高但变压器和线路的网络约束会把一部分可用容量卡住——这就是为什么T1故障时系统要削减23MW负荷而不是简单靠G1G2把90MW全部带起来。3.2 手动复现事件表概率、频率、失电量三步走PPT给出了这个系统所有事件的分析表从所有元件正常、一阶事件、二阶事件一直列到五阶全故障。我拿到这套表之后习惯用脚本把它重新算一遍既验证表里的数值也方便改参数做敏感性分析。下面是一个用Python枚举全部状态并汇总指标的骨架脚本。# 状态解析法枚举脚本验算PPT算例1的事件表 import itertools # 元件配置名称, 可用率, 不可用率 # 可用率和不可用率由故障率λ与修复率μ换算Aμ/(λμ)Uλ/(λμ) comps [ {name: G1, A: 0.98887223, U: 0.01112777}, {name: G2, A: 0.98887223, U: 0.01112777}, {name: T1, A: 0.99933233, U: 0.00066767}, {name: T2, A: 0.99933233, U: 0.00066767}, {name: L1, A: 0.99857565, U: 0.00142435}, ] # 失负荷量字典键为故障元件名值为该状态下的负荷削减量(MW) # 这里按PPT事件表的网络约束分析结果填入实际工程中应由潮流计算给出 loss_dict { (G1,): 40, (T1,): 23, (L1,): 49, (G1, G2): 100, # ... 其余状态按PPT表补全 } states [] for mask in itertools.product([0, 1], repeatlen(comps)): # mask0表示正常1表示故障 prob 1.0 failed [] for i, (comp, bit) in enumerate(zip(comps, mask)): if bit 1: prob * comp[U] failed.append(comp[name]) else: prob * comp[A] loss loss_dict.get(tuple(sorted(failed)), 0) if loss 0: states.append((tuple(failed), prob, loss)) # 汇总可靠性指标 lolp sum(p for _, p, _ in states) eens sum(p * loss * 8760 for _, p, loss in states) print(fLOLP {lolp:.6f}) print(fEENS {eens:.2f} MWh/年)这段脚本的枚举逻辑很简单itertools.product生成每个元件的0/1状态组合0代表可用、1代表故障每个组合的概率由可用率和不可用率连乘得出然后通过失负荷量字典判断这个状态是否失效并取对应的削减量。最后汇总LOLP和EENS。有两点要说明。第一失负荷量字典在PPT里是现成的因为网络约束分析已经做完了如果你要评估一个新系统这一列必须由潮流计算或直流潮流分析产生不能拍脑袋填。第二脚本里未列出的故障组合我认为不失负荷这对应PPT表中失负荷量为0的那些行实际使用时要逐行核对。3.3 表格里的系统可靠性指标LOLP、LOLF、EENS怎么验算PPT算例1最终给出系统可靠性指标LOLP0.013405LOLF2.896923次/年EENS4767.36MWh/年。这三个数怎么来的可以直接从事件表里验。LOLP是全部失负荷事件概率之和。把表里所有失负荷量大于0的事件概率加总就得到0.013405。这意味着该长期运行工况下系统处于失负荷状态的时间占比约为1.34%折算成时间大约是每年117个小时处于失负荷状态。LOLF是失负荷事件的发生频率单位是次/年。这个指标在状态解析法里不是简单把事件概率相加而是每个失效状态的概率乘以该状态的离开率再求和。离开率包括故障元件的修复率和正常元件的故障率所以在算例1里G1故障这个事件的频率不是直接拿G1的故障率2次/年而是要经过状态转移率的换算。EENS的验算我通常会拆成两步。先把每个失效状态的概率乘以失负荷量得到该状态的单次期望失电量再乘以8760小时换算成年化指标。把所有失效事件加总后和PPT给出的4767.36MWh/年对照误差应该在0.5%以内。如果对不上优先检查可用率是否由λ和μ正确换算其次检查失负荷量是否抄错行。提示验算时最常犯的错是把“累积概率”当成“LOLP”。PPT事件表的累积概率是按概率从大到小排序后的累加它逼近的是1.0而LOLP只是其中失负荷事件概率的加总两者含义不同。4. 算例2方案对比二阶事件枚举与增强措施的效果量化4.1 元件可用率计算从λ和μ到A值算例2的系统配置是3台发电机加2回输电线。设每台发电机和每回线路的故障率、修复率分别相等发电机λG2次/年、μG36.5次/年线路λL2次/年、μL365次/年。发电机的可用率Aμ/(λμ)36.5/38.5≈0.94805不可用率U≈0.05195线路的可用率A365/367≈0.99727不可用率U≈0.00273。这里要注意量纲。μG36.5次/年对应的平均修复时间是8760/36.5240小时也就是10天μL365次/年对应的平均修复时间是24小时。发电机的修复时间远长于线路这是符合实际的——发电机故障后检修周期以天计线路故障后巡线、抢修通常能在一天内完成。所以在同一个系统里不同元件的λ、μ必须保持同一单位再做除法否则可用率会算错一个数量级。全部元件正常工作的概率是0.94805的三次方乘以0.99727的平方约等于0.8475。这个值看起来很“安全”但真正决定可靠性指标的是那些低概率高损失的事件组合——比如两台发电机同时故障概率只有0.002545却要削减200MW负荷对EENS的贡献远大于它的概率数字。4.2 二阶事件表解读哪些组合在拖指标PPT按二阶事件列出了系统状态分析表。我挑几个关键行来看事件故障概率频率次/年负荷削减MW失电量MWhG1、G2故障0.0025450.1959382004458.227G1、L2故障0.0001270.051675400445.437G2、G3故障0.0025450.1959382004458.227G2、L2故障0.0001270.051675200222.723L1、L2故障0.0000060.00467440022.250L2单回线故障0.0023200.8630032004064.465注意L2单回线故障这行特别显眼概率只有0.00232但频率0.863003次/年几乎是全部二阶事件频率的一半失电量4064.465MWh排第二仅次于双发电机故障。这说明什么单回L2故障虽然只削减200MW不算最严重但它发生的频率高对系统可靠性的拖累反而比那些概率更低的高阶组合更大。再看G1、L2组合故障削减400MW而G2、L2组合故障只削减200MW。原因是G1是200MW机组L2故障后再叠加G1退出系统缺口的规模更大G2故障叠加L2故障时还有G1和G3两台机组撑着削减量相对小。这提示我们评估发输电组合系统时不能只看“几阶故障”还要看具体是哪几个元件组合——发电机容量差异和线路传输能力差异会直接改变削减量。4.3 增强措施怎么量化扩容与加线的指标对比PPT给出一组增强性措施把线路L1的容量从300MW提升到500MW同时增装一回线路L3L3的故障率和修复率与L1、L2相同。增强前后的系统可靠性指标对比非常直观指标改善前改善后变化幅度LOLP0.0105960.007626下降约28%LOLE小时/年92.768466.804下降约28%LOLF次/年1.7138660.603159下降约65%EENSMWh/年18797.7213371.44下降约29%这个对比表含金量很高。LOLP和EENS都下降了近三成但LOLF下降了近三分之二原因在于增强措施让很多原本会造成削减的单重故障直接变成了“无影响事件”。改善后的事件表里L1单回线故障、L2单回线故障、G1与任意单回线的组合故障失负荷量全都归零。单回线故障不再导致停电整个系统的失负荷频率自然大幅下降。我特别想强调PPT里那句话实际系统中采用增强措施的依据要综合考虑设备投资费用、运行维修费用、线损改变费用和停电损失费用并且要用同一年值口径去比。也就是说把容量从300MW提到500MW、增加一回线路这些举措的决策是否划算不能只看可靠性指标改善了多少还得把停电损失折算成钱和投资费用放同一个时间坐标上做等年值比较。这才是规划评估落地时真正要回答的问题。5. 状态解析法避坑指南五个我差点信了公式的高频翻车点5.1 频率算出来总是对不上表现象用PPT里的频率公式fΣλ·P去算结果和表里的LOLF总差一截差的不是小数而是成倍的偏差。原因频率公式里的λ是“离开该状态的所有转移率之和”不是一个元件故障率那么简单。比如算例2里G1、G2双故障状态离开这个状态有两个途径——G1修复或者G2修复两条路径的转移率都要算进去分别是μG1和μG2。如果只把λG1和λG2加进去方向就反了数值自然对不上。解决先把每个故障元件的修复率μ按8760/修复小时数换算好状态频率状态概率×(Σ故障元件μΣ正常元件λ)再按月汇总。这个口径在工程软件里也是通用的直接背下来用。5.2 把“不可用率”当成“故障率”去算概率现象算元件故障概率时直接拿λ2次/年当作概率代入连乘得出一个比正确结果大几十倍的LOLP。原因λ的量纲是次/年它是转移率不是概率。只有当λ数值极小且时间窗取1年时λ才近似等于概率。发输电系统里发电机λ2次/年对应不可用率0.05195差了十几倍线路λ2次/年对应不可用率只有0.00273差了三个数量级完全不能混用。解决所有状态概率必须用A和U计算而A、U由λ和μ共同决定。动手前先列一张元件参数换算表把每个元件的A、U算好再进枚举脚本。我自己的习惯是固定用“Aμ/(λμ)Uλ/(λμ)”不搞近似。5.3 失负荷量用“容量差”代替“网架约束”现象评估算例1时G2故障明明只剩G1100MW高于负荷90MW但把G2故障当作无影响事件结果算出来的LOLP比PPT少了快一半。原因算例1里T1故障要削减23MW、L1故障要削减49MW这些削减不是发电容量不够造成的而是变压器、线路的传输能力约束造成的。发输电组合系统评估的失负荷量必须来自网络分析不能简单用发电容量减负荷。解决把失负荷量作为枚举环节的外部输入由潮流计算或直流潮流分析提供。PPT里的事件表已经把网络约束结果固化下来了我用脚本复现时直接查表即可新场景下必须先跑网络分析再回填。注意状态解析法的枚举可以自动化但“判定是否失效、削减多少”永远是网络分析问题。把这两层解耦程序才不会越写越乱。5.4 累积概率排序列越看越糊涂现象把事件表按累积概率排序后发现LOLP直接等于1减去某一个累积概率不对怎么减都对不上。原因事件表的累积概率是按“单个状态概率从大到小”排序后的累加和系统是否失负荷没有对应关系。全正常状态概率最大排第一二阶、三阶事件按概率大小插入其中累积概率只是用来查“概率分布尾部”的参考列不是用来反推LOLP的。解决LOLP只统计失负荷事件把表中失负荷量大于0的行挑出来概率直接相加。不要试图从累积概率列反推那列是给人看概率分布的不是给公式用的。验算时我通常让脚本重新按失负荷事件单独求和和表尾给的LOLP值对照。5.5 二阶截断之后发现高阶事件漏掉了大额失电量现象某个系统按二阶截断算完EENS看起来不高但把三阶、四阶事件补进去后EENS涨了一倍。原因虽然高阶事件概率小但如果削减量大概率×削减量的乘积未必可以忽略。算例1里四阶事件概率已经到10的负8次方量级但失负荷100MW甚至90MW贡献的失电量还是微小的吗——到10的负8次方量级×90MW确实微小但如果系统冗余低、低阶故障概率被网络约束压掉高阶组合反而成为主要的失效模式。解决截断前先做一版无截断枚举观察EENS随阶数的收敛速度如果三阶事件对EENS的贡献超过5%就得把截断阶数提高或者改用解析网络法、蒙特卡洛法。PPT算例2明确只讨论二阶事件是因为这个系统二阶事件已经覆盖了全部主要失效模式这是个案不能直接推广。6. 从发输电到配电网用FMEA法反向验证枚举结果FMEA法故障模式与影响分析是配电网可靠性评估的经典方法和状态解析法看起来完全不同骨子里其实是同一件事。状态解析法枚举的是“元件状态组合”FMEA枚举的是“元件故障事件”两者的核心动作都是“故障事件→分析影响→累计指标”。配电网通常辐射状结构网络分析比输电网简单但开关设备的存在让故障隔离范围成为关键变量。断路器能切负荷电流和短路电流故障时与保护配合快速切除隔离开关只能在无负荷电流时分合提供检修断点负荷开关介于两者之间能切额定负荷电流但不能切短路电流。这三个设备在FMEA表里决定了一个元件故障后“停电范围有多大、多久能恢复”。我在做配网可靠性评估时会把每个开关的动作逻辑画成一张影响表再和状态解析法的枚举结果交叉验证——如果某个元件的FMEA影响分析和状态枚举给出的失负荷量矛盾那一定是开关配置或者枚举判定里有一处错了。从那以后我每次拿到可靠性评估结果都强制自己走一遍“先验概率、再验频率、最后对LOLP量级”的三板斧先确认所有概率都由A、U连乘得出再确认频率用了离开率求和最后用EENS8760×Σ(概率×失负荷量)做闭环验算。三个数能互相咬合才敢把这份算例的结论拿到规划方案里去用。希望帮到你。本文还有配套的精品资源点击获取