☰
发输电组合系统可靠性评估:状态解析法算例详解与LOLP指标计算
2026/9/30 8:08:16 网站建设 项目流程

简介:这份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最终给出EENS=4767.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负荷,而不是简单靠G1+G2把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], repeat=len(comps)): # mask=0表示正常,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(f"LOLP = {lolp:.6f}") print(f"EENS = {eens:.2f} MWh/年")

这段脚本的枚举逻辑很简单:itertools.product生成每个元件的0/1状态组合,0代表可用、1代表故障;每个组合的概率由可用率和不可用率连乘得出;然后通过失负荷量字典判断这个状态是否失效,并取对应的削减量。最后汇总LOLP和EENS。

有两点要说明。第一,失负荷量字典在PPT里是现成的,因为网络约束分析已经做完了;如果你要评估一个新系统,这一列必须由潮流计算或直流潮流分析产生,不能拍脑袋填。第二,脚本里未列出的故障组合我认为不失负荷,这对应PPT表中失负荷量为0的那些行,实际使用时要逐行核对。

3.3 表格里的系统可靠性指标:LOLP、LOLF、EENS怎么验算

PPT算例1最终给出系统可靠性指标:LOLP=0.013405,LOLF=2.896923次/年,EENS=4767.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回输电线。设每台发电机和每回线路的故障率、修复率分别相等:发电机λG=2次/年、μG=36.5次/年;线路λL=2次/年、μL=365次/年。发电机的可用率A=μ/(λ+μ)=36.5/38.5≈0.94805,不可用率U≈0.05195;线路的可用率A=365/367≈0.99727,不可用率U≈0.00273。

这里要注意量纲。μG=36.5次/年对应的平均修复时间是8760/36.5=240小时,也就是10天;μL=365次/年对应的平均修复时间是24小时。发电机的修复时间远长于线路,这是符合实际的——发电机故障后检修周期以天计,线路故障后巡线、抢修通常能在一天内完成。所以在同一个系统里,不同元件的λ、μ必须保持同一单位再做除法,否则可用率会算错一个数量级。

全部元件正常工作的概率是0.94805的三次方乘以0.99727的平方,约等于0.8475。这个值看起来很“安全”,但真正决定可靠性指标的是那些低概率高损失的事件组合——比如两台发电机同时故障,概率只有0.002545,却要削减200MW负荷,对EENS的贡献远大于它的概率数字。

4.2 二阶事件表解读:哪些组合在拖指标

PPT按二阶事件列出了系统状态分析表。我挑几个关键行来看:

事件故障概率频率(次/年)负荷削减(MW)失电量(MWh)
G1、G2故障0.0025450.1959382004458.227
G1、L2故障0.0001270.051675400445.437
G2、G3故障0.0025450.1959382004458.227
G2、L2故障0.0001270.051675200222.723
L1、L2故障0.0000060.00467440022.250
L2单回线故障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,同时增装一回线路L3,L3的故障率和修复率与L1、L2相同。增强前后的系统可靠性指标对比非常直观:

指标改善前改善后变化幅度
LOLP0.0105960.007626下降约28%
LOLE(小时/年)92.768466.804下降约28%
LOLF(次/年)1.7138660.603159下降约65%
EENS(MWh/年)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故障明明只剩G1=100MW,高于负荷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连乘得出,再确认频率用了离开率求和,最后用EENS=8760×Σ(概率×失负荷量)做闭环验算。三个数能互相咬合,才敢把这份算例的结论拿到规划方案里去用。希望帮到你。

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

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询