事故树分析PPT.ppt

例:以如图事故树为例，用最小径集法计算顶事件的发生概率。设各基本事件的发生概率q1=q2=q3=q4=q5=0.1. + · X2 X5 M2 T M1 + X1 M3 · X3 X5 + X4 M4 · X3 M5 解： 事故树有四个最小径集： · T + X1 X5 K2 + X1 X3 K1 + X2 K4 X5 X4 + X4 K3 X3 化简，代入具体数值,得: p(T)=0.02071 和前种方法计算结果一致 练习:以如图事故树为例，试用最小径集法写出顶事件的发生概率的表达式。设各基本事件的发生概率为q1 , q2 ,q3 ,q4 · T + X1 X4 K2 + X1 X3 K1 + X2 K4 X4 + X4 K3 X3 解：该事故树有四个最小径集： 由最小径集公式得顶事件发生概率： · T + X1 X4 K2 + X1 X3 K1 + X2 K4 X4 + X4 K3 X3 ?小结 在上述三种顶事件发生概率的精确算法中，状态值表法计算工作量一般较大，花费时间较长，所以一般不适合手工运算，但由于其规律性强，适用于计算机编制程序计算。 后两种方法相对较简单。虽然最小径集法计算公式较最小割集法要复杂一些，但一般情况下，由于事故树的最小径集的数量往往少于最小割集，所以最小径集法更为实用。 5.化相交集合为不相交集合理论 从最小割集法和最小径集法公式的计算项数看，两式的和差项分别为(2k-1)与2k项。当k足够大时，就会产生“组合爆炸”问题。 如K=40，则计算P(T)的公式和差项共有 240-1=1.1×1012 每一项又是许多数的连乘积，即使计算机也难以胜任。解决的办法就是首先求解出事故树的最小割集，再化相交和为不交和，求顶事件发生概率的精确解。 T Er · X3 + Es X4 · X1 X2 P(T) = P(Er∪Es )=P( Er )+P(Es )-P(Er∩ Es ) 而一般, 即公式第三项 P(Er∩ Es ) ≠0 Er、Es相容 Er● Es ≠Φ 思考：如何计算如图事故树顶上事件T的发生概率。已知x1，x2，x3，x4的发生概率为q1,q2,q3,q4。 逻辑加的基本公式适用吗？ 独立事件 证明：Er+Es=(Er+Er)(Er+Es) =ErEr+ErEs+ErEr+ErEs=Er+ErEs ? Er+Es=Er+ErEs P(T)=P(Er+ Es )=P(Er+ErEs ) Er∩Er Es = Φ P(Er∩Er Es)=0 为简化计算，将Er+Es转化为 Er与Er/ Es互不相容是关键 =P(Er )+P(Er Es )-P(Er∩Er Es ) = P(Er )+P(Er Es ) 某事故树有k个最小割集：E1，E2，…，Er，…， Ek，一般情况下它们是相容的，即两个最小割集可能同时发生。由文氏图可以看出， Er与Es相容，Er与ErEs不相容，如下图所示： 上式中，任意两项都是互不相容的,如前两项: P(T)= P(E1+E2+E3+···+EK ) =P(E1+E1/E2+E1/E2/E3+···+E1/E2/E3/···Ek-1/Ek) = P( E1)+P(E1/E2)+P(E1/E2/E3) +··· +P (E1/E2/E3/···Ek-1/Ek) E1∩E1E2 =Φ P( E1∩E1E2 )=0 P( E1∪E1E2)= P( E1)+P(E1E2)-P( E1∩E1E2) = P( E1)+P(E1E2) 所以得到 由上式可以推广到一般式： T=E1+E2+E3+···+EK =E1+E1E2+E1E2E3+···+E1E2E3···Ek-1Ek 命题2 若集合Er和Es包含共同元素，则： Er/Es= E/r←sEs Er←s表示Er中有的而Es中没有的元素的布尔积. 如Er={X1,X2,X3}, Es={X1,X2,X4},则: Er←s=X3 Er/Es=(X1X2X3)/(X1X2X4) E/r←sEs=X3/(X1X2X4) P(T)= P(E1+E2+E3+···+EK ) = P(E1)+P(E1/E2)+P(E1/E2/E3) +···+P(E1/E2/E3/···Ek-1/Ek) 为了简化上式，提出下面的四个命题： 命题1 集合Er和Es如不包含共同元素，则Er/Es可用不交化规则直接展开。 命题3 若集合Er 和Et包含共同元素，Es和Et 也包含共同元素，则： Er/Es/Et= E/r←tE/s←tEt 命题4 若集合Er和Et包含共同要