一种场馆供电系统可靠性分析方法、装置、设备及介质与流程

1.本发明属于场馆供电技术领域，具体涉及一种场馆供电系统可靠性分析方法、装置、设备及介质。


背景技术：

2.随着经济的快速发展和综合国力的不断提升，大型场馆、博览中心、会场和公共区域作为社会活动的主体，承担着确保活动圆满成功的重任，其供电可靠性直接关系到社会影响，且活动期间发生停电事故会产生不可估量的经济损失。因此，亟需研究科学合理的考虑多种供电保障措施的场馆供电复杂冗余系统可靠性评价方法，以确保供电系统的安全稳定运行以及终端负荷的可靠供电。
3.场馆供电系统的供电保障措施常涉及重要低压负荷处ups、备用电源、元件备用等系统内部复杂的逻辑结构、动态顺序特性与部件依赖关系问题。工程应用中普遍采用markov方法来获取系统状态以定量评价系统可靠性，但复杂系统包含大量元件，使得马尔科夫链在多阶故障下其状态维数将急剧上升。此外需重点研究系统的动态变化过程，以及维修、故障等行为对系统故障特性的影响。传统静态故障树分析方法采用布尔逻辑门来表示故障事件间的组合关系，没有考虑事件间的动态依赖关系。考虑多种供电保障措施的场馆供电系统存在动态时间历程、部件间相互作用和系统状态划分等求解难题，动态故障树在传统故障树分析方法基础上考虑了故障事件间的动态特性，因此增强了对系统可靠性分析的适用度。


技术实现要素：

4.本发明的目的在于提供一种场馆供电系统可靠性分析方法、装置、设备及介质，实现考虑系统顺序相关性、可修复性、元件备用情况和应急电源系统等供电保障措施及其动态特征的场馆供电系统可靠性分析评价。
5.本发明采取下述技术方案予以实现：
6.第一方面，一种场馆供电系统可靠性分析方法，包括以下步骤：
7.s1、依据场馆相关运行数据及负荷需求，建立保障不同场景情况下对负荷持续可靠供电的供电系统结构图；
8.s2、根据所述的供电系统结构图，判断内部逻辑关系，并引入动态逻辑门，建立以负荷供电中断故障为顶事项的动态故障树模型；
9.根据动态故障树模型内的10kv上级电源动态子树模块建立马尔科夫模型，并计算所述马尔科夫模型各状态的稳态概率；
10.s4、采用序贯蒙特卡罗仿真方法基于得到的马尔科夫模型各状态的稳态概率，对所述马尔科夫模型和所述动态故障树模型进行时序状态序列模拟，确定供电系统运行状态；
11.s5、根据所述供电系统运行状态计算可靠性评价指标，从而对供电系统可靠性进
行评价。
12.本发明的进一步改进在于：所述供电系统包括若干负载、10kv上级电源供电系统、0.4kv备用电源回路、0.4kv应急电源回路和不间断电源ups；
13.若干所述负载包括，负载a1、负载a2、负载a3、负载a4、负载a5和负载a6；
14.所述上级电源供电系统包括两路10kv市电主供电源s1、s2及一路10kv市电备供电源s3；
15.所述10kv市电备供电源s3与10kv母线b3通过馈线l3和断路器q3相连；
16.所述10kv市电主供电源s1与10kv母线b1通过馈线l1和断路器q1相连；
17.所述10kv市电主供电源s2与10kv母线b2通过馈线l2和断路器q2相连；
18.所述10kv母线b3与10kv母线b1通过断路器q4相连；
19.所述10kv母线b3与10kv母线b2通过断路器q5相连；
20.所述10kv母线b1与变流器t1相连；
21.所述10kv母线b2与变流器t2相连；
22.所述变流器t1与0.4kv母线b4通过馈线l4和断路器q6相连；
23.所述变流器t2与0.4kv母线b5通过馈线l5和断路器q7相连；
24.所述0.4kv母线b4与0.4kv母线b5通过断路器q8相连；
25.所述0.4kv母线b4与负载a2通过馈线l11和断路器q16相连；
26.所述0.4kv母线b4与双电源转换装置ats4通过馈线l10和断路器q15相连；
27.所述0.4kv母线b5与负载a3通过馈线l13和断路器q18相连；
28.所述0.4kv母线b5与双电源转换装置ats3通过馈线l12和断路器q17相连；
29.所述双电源转换装置ats4通过馈线l19与负载a6相连；
30.所述0.4kv备用电源回路包括两台备用电源g1和g2；
31.所述备用电源g1通过断路器q9与汇流排b6相连；
32.所述备用电源g2通过断路器q10与汇流排b6相连；
33.所述汇流排b6通过断路器q11和馈线l7与双电源转换装置ats1相连；
34.所述0.4kv母线b5通过馈线l6与双电源转换装置ats1相连；
35.所述双电源转换装置ats1与备用母线段b7相连；
36.所述备用母线段b7与双电源转换装置ats4通过馈线l14和断路器q19相连；
37.所述备用母线段b7与负载a4通过馈线l15和断路器q20相连；
38.所述0.4kv应急电源回路包括，两台应急电源g3、g4；
39.所述应急电源g3通过断路器q13与汇流排b8相连；
40.所述应急电源g4通过断路器q12与汇流排b8相连；
41.所述汇流排b8通过断路器q14和馈线l9与双电源转换装置ats2相连；
42.所述0.4kv母线b4通过馈线l8与双电源转换装置ats2相连；
43.所述双电源转换装置ats2与备用母线段b9相连；
44.所述备用母线段b9与双电源转换装置ats3通过馈线l16和断路器q21相连；
45.所述备用母线段b9与负载a1通过馈线l176和断路器q22相连；
46.所述不间断电源ups设置在双电源转换装置ats3与负载a5之间。
47.本发明的进一步改进在于：所述s2具体包括以下步骤：
48.s21、动态故障树模型顶事件选择，以负荷a1、a2、a3、a4、a5或a6任一个供电中断为顶事件k；
49.s22、确定10kv上级电源动态子树模块，以10kv市电主供电源s1、断路器q1和馈线故障l1作为底事件p1；10kv市电主供电源s2、断路器q2和馈线故障l2作为底事件p2；市电备供电源s3断路器q3和馈线故障l3作为底事件p3；
50.s23、确定0.4kv备用电源回路静态子树模块，以备用电源g1、备用电源g2、汇流排b6、双电源转换装置ats1、备用母线段b7及其相连的断路器、与馈线故障回路故障作为底事件p4；
51.确定0.4kv应急电源回路静态子树模块，以应急电源g3、应急电源g4、汇流排b8、双电源转换装置ats2、应急母线段b9及其相连的断路器与馈线故障回路故障作为底事件p5；
52.s24、将供电系统中除底事件p1、p2、p3、p4和p5外的故障事件作为触发底事件p6；
53.s25、将不间断电源ups、双电源转换装置ats3、双电源转换装置ats4及其相连的断路器与馈线故障作为底事件p7；
54.s26、基于底事件p1-p7引入动态逻辑门建立供电系统动态故障树模型。
55.本发明的进一步改进在于：所述动态逻辑门包括hsp热备件门与fdep功能相关门。
56.本发明的进一步改进在于：所述s3具体包括以下步骤：
57.s31、根据动态故障树模型内的10kv上级电源动态子树模块建立马尔科夫模型，并根据马尔科夫模型绘制马尔科夫状态空间图；
58.s32、根据马尔科夫状态空间图建立状态转移矩阵a：
[0059][0060]
式中，λ
12
、λ
23
、λ
24
、λ
35
和λ
45
为马尔科夫模型各状态的故障率；γ
12
、γ
23
、γ
24
、γ
35
和γ
45
为马尔科夫模型各状态的修复率；
[0061]
s33、根据状态转移矩阵a计算线性代数方程组：
[0062][0063]
式中，p为马尔科夫模型各状态的稳态概率，p＝{pi,i＝1,2,3,4,5}其中p5为马尔科夫模型的故障率。
[0064]
本发明的进一步改进在于：在s4中，具体包括以下步骤：
[0065]
s41、采用序贯蒙特卡罗仿真法模拟马尔科夫模型和底事件p3、p4、p5、p6和p7的抽样时间序列，底事件bi第j次抽样时间具体公式为
[0066][0067]
式中和分别为底事件bi在第j次的失效抽样时间和修复抽样时间；所述和为在区间[0,1]上的随机数；所述f
i,op-1
和f
i,fa-1
为第i个底事件的修复分布函数和失效分布函数；
[0068]
底事件bi的序列抽样时间t
i,j
：
[0069][0070]
式中，t
max为
最大仿真时间；
[0071]
底事件bi的抽样时间序列ti：
[0072]
ti＝{t
i,1
,t
i,2
,l,t
i,j
}。
[0073]
s42、计算各底事件在最大仿真时间t
max
内的时序状态序列，底事件bi运行状态函数：
[0074][0075]
式中，为底事件bi在抽样时间区间t
sim
∈[t
i,j-1
,t
i,j
)内的运行状态；所述t
sim
为仿真时间
[0076]
s43、根据底事件bi的抽样时间序列ti，获得最大仿真时间t
max
内的抽样时间序列t
′
，将所有底事件的抽样时间序列元素进行升序排列，得到
[0077]
t
′
＝{t
′1,t
′2,l,t
′m}；
[0078]
式中，0＝t
′0≤t
′1≤t
′2≤l≤t
′m≤t
max
；
[0079]
s44、计算在各抽样时间区间内的运行状态时间序列，并记录系统正常运行时间t
op,l
和次数a
l
，故障修复时间t
fa,l
和故障次数r
l
；
[0080]
动态故障树结构函数为：
[0081][0082]
所述φ
l
为第l个仿真时间区间[t
i,j-1
,t
i,j
)内系统的运行状态；
[0083]
运行状态时间序列：
[0084]
φ＝{φ1,φ2,l,φ
l
}
[0085]
第l个仿真时间区间t
sim
∈[t
i,j-1
,t
i,j
)；
[0086]
运行状态记录公式为
[0087][0088][0089]
式中，t
fa，l
表示故障修复时间，t
op，l
表示正常运行时间，r
l
表示故障修复次数，a
l
表示正产运行次数，l表示，供电系统在状态时间序列的单次仿真中，第l次状态改变的编号。
[0090]
本发明的进一步改进在于：所述s5具体包括以下步骤：
[0091]
s51、计算累计故障失效概率作为可靠性评价指标：
[0092][0093]
式中，n为仿真时段t
sim
≤t内的故障次数；t为变量，改变t的大小绘制得到累计故障率的时间变化曲线；
[0094]
s52、计算平均修复时间和平均无故障时间作为可靠性评价指标；
[0095]
所述平均修复时间t
mttr
：
[0096][0097]
式中，r为最大仿真时间t
max
内的故障修复次数；
[0098]
所述平均无故障时间t
mtbf
：
[0099][0100]
式中，a为最大仿真时间t
max
内的正常运行次数；
[0101]
s53、计算系统稳态可用度用于评价场馆供电系统长期运行的可靠性水平，取系统稳态可用度a(∞)：
[0102][0103]
第二方面，一种场馆供电系统可靠性分析装置，包括：
[0104]
供电系统建立模块：用于依据场馆相关运行数据及负荷需求，建立保障不同场景情况下对负荷持续可靠供电的供电系统结构图；
[0105]
动态故障树模型建立模块：用于根据所述的供电系统结构图，判断内部逻辑关系，并引入动态逻辑门，建立以负荷供电中断故障为顶事项的动态故障树模型；
[0106]
稳态概率计算模块：用于：根据马尔科夫理论和动态故障树模型计算马尔科夫模
型各状态的稳态概率；
[0107]
供电系统运行状态确认模块：用于通过序贯蒙特卡罗仿真方法对所述马尔科夫模型和所述动态故障树模型进行时序状态序列模拟，确定供电系统运行状态
[0108]
供电系统可靠性评价模块：用于根据所述供电系统运行状态计算可靠性评价指标，从而对供电系统可靠性进行评价。
[0109]
第三方面，一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述的一种场馆供电系统可靠性分析方法。
[0110]
第四方面，一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时上述的一种场馆供电系统可靠性分析方法。
[0111]
与现有技术相比，本发明至少包括以下有益效果：
[0112]
1、本发明提供一种综合故障树分析和马尔科夫模型的供电系统可靠性评价模型，考虑系统顺序相关性、可修复性、备件性等特征，并采用时序蒙特卡罗仿真方法解决马尔科夫模型在多阶故障下状态维数急剧上升的问题，具有强适应性和解决问题能力；
[0113]
2、本发明将可靠性评价对象拓展到终端用户，着重考虑场馆供电系统中元件备用情况、应急电源系统和不间断电源ups等低电压等级配电层面的供电保障措施，提高对场馆供电系统可靠性分析评价方法的针对性；
[0114]
3、本发明通过引入动态逻辑门以体现供电保障系统内部复杂的逻辑关系与动态特性，进行场馆供电系统相关动态分析，更清晰直观地描述了系统故障失效的动态行为，提高对场馆供电系统可靠性分析方法的适用性。
附图说明
[0115]
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
[0116]
图1为本发明一种场馆供电系统可靠性分析方法中供电系统结构示意图；
[0117]
图2为本发明一种场馆供电系统可靠性分析方法中动态故障树模型示意图；
[0118]
图3为本发明一种场馆供电系统可靠性分析方法中hsp热备件门示意图；
[0119]
图4为本发明一种场馆供电系统可靠性分析方法中hsp热备件门的马尔科夫状态空间图；
[0120]
图5为本发明一种场馆供电系统可靠性分析方法的流程图；
[0121]
图6为本发明一种场馆供电系统可靠性分析装置的结构框图。
具体实施方式
[0122]
下面将参考附图并结合实施例来详细说明本发明。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
[0123]
以下详细说明均是示例性的说明，旨在对本发明提供进一步的详细说明。除非另有指明，本发明所采用的所有技术术语与本发明所属领域的一般技术人员的通常理解的含义相同。本发明所使用的术语仅是为了描述具体实施方式，而并非意图限制根据本发明的示例性实施方式。
[0124]
实施例1
[0125]
如图5所示，一种场馆供电系统可靠性分析方法，包括以下步骤：
[0126]
s1、依据场馆相关运行数据及负荷需求，建立保障不同场景情况下对负荷持续可靠供电的供电系统结构图；
[0127]
s2、根据的供电系统结构图，判断内部逻辑关系，并引入动态逻辑门，建立以负荷供电中断故障为顶事项的动态故障树模型；
[0128]
s3、根据动态故障树模型内的10kv上级电源动态子树模块建立马尔科夫模型，并计算所述马尔科夫模型各状态的稳态概率；
[0129]
s4、采用序贯蒙特卡罗仿真方法基于得到的马尔科夫模型各状态的稳态概率，对马尔科夫模型和动态故障树模型进行时序状态序列模拟，确定供电系统运行状态；
[0130]
s5、根据供电系统运行状态计算可靠性评价指标，从而对供电系统可靠性进行评价。
[0131]
如图1所示，供电系统包括若干负载、10kv上级电源供电系统、0.4kv备用电源回路、0.4kv应急电源回路和不间断电源ups；
[0132]
若干负载包括，负载a1、负载a2、负载a3、负载a4、负载a5和负载a6；
[0133]
上级电源供电系统包括两路10kv市电主供电源s1、s2及一路10kv市电备供电源s3；
[0134]
10kv市电备供电源s3与10kv母线b3通过馈线l3和断路器q3相连；
[0135]
10kv市电主供电源s1与10kv母线b1通过馈线l1和断路器q1相连；
[0136]
10kv市电主供电源s2与10kv母线b2通过馈线l2和断路器q2相连；
[0137]
10kv母线b3与10kv母线b1通过断路器q4相连；
[0138]
10kv母线b3与10kv母线b2通过断路器q5相连；
[0139]
10kv母线b1与变流器t1相连；
[0140]
10kv母线b2与变流器t2相连；
[0141]
变流器t1与0.4kv母线b4通过馈线l4和断路器q6相连；
[0142]
变流器t2与0.4kv母线b5通过馈线l5和断路器q7相连；
[0143]
0.4kv母线b4与0.4kv母线b5通过断路器q8相连；
[0144]
0.4kv母线b4与负载a2通过馈线l11和断路器q16相连；
[0145]
0.4kv母线b4与双电源转换装置ats4通过馈线l10和断路器q15相连；
[0146]
0.4kv母线b5与负载a3通过馈线l13和断路器q18相连；
[0147]
0.4kv母线b5与双电源转换装置ats3通过馈线l12和断路器q17相连；
[0148]
双电源转换装置ats4通过馈线l19与负载a6相连；
[0149]
0.4kv备用电源回路包括两台备用电源g1和g2；
[0150]
备用电源g1通过断路器q9与汇流排b6相连；
[0151]
备用电源g2通过断路器q10与汇流排b6相连；
[0152]
汇流排b6通过断路器q11和馈线l7与双电源转换装置ats1相连；
[0153]
0.4kv母线b5通过馈线l6与双电源转换装置ats1相连；
[0154]
双电源转换装置ats1与备用母线段b7相连；
[0155]
备用母线段b7与双电源转换装置ats4通过馈线l14和断路器q19相连；
[0156]
备用母线段b7与负载a4通过馈线l15和断路器q20相连；
[0157]
0.4kv应急电源回路包括，两台应急电源g3、g4；
[0158]
应急电源g3通过断路器q13与汇流排b8相连；
[0159]
应急电源g4通过断路器q12与汇流排b8相连；
[0160]
汇流排b8通过断路器q14和馈线l9与双电源转换装置ats2相连；
[0161]
0.4kv母线b4通过馈线l8与双电源转换装置ats2相连；
[0162]
双电源转换装置ats2与备用母线段b9相连；
[0163]
备用母线段b9与双电源转换装置ats3通过馈线l16和断路器q21相连；
[0164]
备用母线段b9与负载a1通过馈线l176和断路器q22相连；
[0165]
不间断电源ups设置在双电源转换装置ats3与负载a5之间。
[0166]
其中10kv电源应至少由两个相对独立的上级变电站供电，并采用两路主供电源同时工作的运行方案；
[0167]
当10kv上级电源故障供电中断时，由配备的双路电源投切装置ats1切换供电来源为备用负荷供电；
[0168]
当10kv上级电源故障供电中断时，转换开关ats2为应急负荷供电；
[0169]
负荷采用不间断电源ups及双电源末端切换供电方式，实现电源转换零切换，满足负荷的高可靠性要求。
[0170]
动态故障树模型，如图2所示，s2具体包括以下步骤：
[0171]
s21、动态故障树模型顶事件选择，以负荷a1、a2、a3、a4、a5或a6任一个供电中断为顶事件k；
[0172]
s22、确定10kv上级电源动态子树模块，以10kv市电主供电源s1、断路器q1和馈线故障l1作为底事件p1；10kv市电主供电源s2、断路器q2和馈线故障l2作为底事件p2；市电备供电源s3断路器q3和馈线故障l3作为底事件p3；
[0173]
s23、确定0.4kv备用电源回路静态子树模块，以备用电源g1、备用电源g2、汇流排b6、双电源转换装置ats1、备用母线段b7及其相连的断路器、与馈线故障回路故障作为底事件p4；
[0174]
确定0.4kv应急电源回路静态子树模块，以应急电源g3、应急电源g4、汇流排b8、双电源转换装置ats2、应急母线段b9及其相连的断路器与馈线故障回路故障作为底事件p5；
[0175]
s24、将供电系统中除底事件p1、p2、p3、p4和p5外的故障事件作为触发底事件p6；
[0176]
s25、将不间断电源ups、双电源转换装置ats3、双电源转换装置ats4及其相连的断路器与馈线故障作为底事件p7；
[0177]
s26、基于底事件p1-p7引入动态逻辑门建立供电系统动态故障树模型。
[0178]
动态逻辑门包括hsp热备件门与fdep功能相关门。
[0179]
带公用备件的热备件门示意图如图3所示，两个带公用备件的热备件门中其备件可替换任意一个故障部件，将其转化为如图4所示的markov状态空间图；
[0180]
s3具体包括以下步骤：
[0181]
s31、根据动态故障树模型内的10kv上级电源动态子树模块建立马尔科夫模型，并根据马尔科夫模型绘制马尔科夫状态空间图，如图4所示；
[0182]
s32、根据马尔科夫状态空间图建立状态转移矩阵a：
[0183][0184]
式中，λ
12
、λ
23
、λ
24
、λ
35
和λ
45
为马尔科夫模型各状态的故障率；γ
12
、γ
23
、γ
24
、γ
35
和γ
45
为马尔科夫模型各状态的修复率；
[0185]
s33、根据状态转移矩阵a计算线性代数方程组：
[0186][0187]
式中，p为马尔科夫模型各状态的稳态概率，p＝{pi,i＝1,2,3,4,5}分别为马尔科夫模型各状态的稳态概率，其中p5为马尔科夫模型的故障率。
[0188]
在s4中，具体包括以下步骤：
[0189]
s41、采用序贯蒙特卡罗仿真法模拟马尔科夫模型和底事件p3、p4、p5、p6和p7的抽样时间序列，底事件bi第j次抽样时间具体公式为
[0190][0191]
式中和分别为底事件bi在第j次的失效抽样时间和修复抽样时间；和为在区间[0,1]上的随机数；f
i,op-1
和f
i,fa-1
为第i个底事件的修复分布函数和失效分布函数；
[0192]
底事件bi的序列抽样时间t
i,j
：
[0193][0194]
式中，t
max为
最大仿真时间；
[0195]
底事件bi的抽样时间序列ti：
[0196]
ti＝{t
i,1
,t
i,2
,l,t
i,j
}。
[0197]
s42、计算各底事件在最大仿真时间t
max
内的时序状态序列，底事件bi运行状态函数：
[0198][0199]
式中，为底事件bi在抽样时间区间t
sim
∈[t
i,j-1
,t
i,j
)内的运行状态；t
sim
为仿真时间
[0200]
s43、根据底事件bi的抽样时间序列ti，获得最大仿真时间t
max
内的抽样时间序列t
′
，将所有底事件的抽样时间序列元素进行升序排列，得到
[0201]
t
′
＝{t
′1,t
′2,l,t
′m}；
[0202]
式中，0＝t
′0≤t
′1≤t
′2≤l≤t
′m≤t
max
；
[0203]
s44、计算在各抽样时间区间内的运行状态时间序列，并记录系统正常运行时间t
op,l
和次数a
l
，故障修复时间t
fa,l
和故障次数r
l
；
[0204]
动态故障树结构函数为：
[0205][0206]
φ
l
为第l个仿真时间区间[t
i,j-1
,t
i,j
)内系统的运行状态；
[0207]
运行状态时间序列：
[0208]
φ＝{φ1,φ2,l,φ
l
}
[0209]
第l个仿真时间区间t
sim
∈[t
i,j-1
,t
i,j
)；
[0210]
运行状态记录公式为
[0211][0212][0213]
式中，t
fa，l
表示故障修复时间，t
op，l
表示正常运行时间，r
l
表示故障修复次数，a
l
表示正产运行次数，l表示，供电系统在状态时间序列的单次仿真中，第l次状态改变的编号。
[0214]
s5具体包括以下步骤：
[0215]
s51、计算累计故障失效概率作为可靠性评价指标：
[0216][0217]
式中，n为仿真时段t
sim
≤t内的故障次数；t为变量，改变t的大小绘制得到累计故障率的时间变化曲线；
[0218]
s52、计算平均修复时间和平均无故障时间作为可靠性评价指标；
[0219]
平均修复时间t
mttr
：
[0220][0221]
式中，r为最大仿真时间t
max
内的故障修复次数；
[0222]
平均无故障时间t
mtbf
：
[0223][0224]
式中，a为最大仿真时间t
max
内的正常运行次数；
[0225]
s53、计算系统稳态可用度用于评价场馆供电系统长期运行的可靠性水平，取系统稳态可用度a(∞)：
[0226][0227]
实施例2
[0228]
如图6所示，一种场馆供电系统可靠性分析装置，包括：
[0229]
供电系统建立模块：用于依据场馆相关运行数据及负荷需求，建立保障不同场景情况下对负荷持续可靠供电的供电系统结构图；
[0230]
动态故障树模型建立模块：用于根据的供电系统结构图，判断内部逻辑关系，并引入动态逻辑门，建立以负荷供电中断故障为顶事项的动态故障树模型；
[0231]
稳态概率计算模块：用于：根据马尔科夫理论和动态故障树模型计算马尔科夫模型各状态的稳态概率；
[0232]
供电系统运行状态确认模块：用于通过序贯蒙特卡罗仿真方法对马尔科夫模型和动态故障树模型进行时序状态序列模拟，确定供电系统运行状态
[0233]
供电系统可靠性评价模块：用于根据供电系统运行状态计算可靠性评价指标，从而对供电系统可靠性进行评价。
[0234]
实施例3
[0235]
一种计算机设备，包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述的一种场馆供电系统可靠性分析方法。
[0236]
实施例4
[0237]
一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时上述的一种场馆供电系统可靠性分析方法。
[0238]
由技术常识可知，本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此，上述公开的实施方案，就各方面而言，都只是举例说明，并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。
[0239]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序
产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0240]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0241]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0242]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0243]
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。