安全系统工程事故树分析(3)事故树分析进一步讨论曹庆贵教授/博士/博导2008年10月习题课:事故树分析一、锅炉超压事故分析(课本第72~76页)1事故树事故树:图3-18。基本事件发生概率:表3-16。锅炉超压C2C3D1D2E1X1X2X3X9X7X8X10X11X12X13X14X15X6X4X5安全阀未卸压不排气排起不够安全阀锈死未定期进行手动试验调压过高规格小起跳高度不够压力超过安全值压力上升未减弱燃烧压力视值偏低工作失误未监视脱岗读数误差安装位置不当表盘直径小表盘看不清光线不足未标红线压力表损坏压力表失灵超期未校未定期冲洗三通关闭TA1A2B1B2B3C1图3-18锅炉超压事故分析图代号基本事件名称qi1-qiX1未定期做手动试验10-20.99X2调压过高10-40.9999X3安全阀规格选小10-30.999X4安全阀起跳高度不够10-30.999X5压力上升5×10-20.95X6压力表损坏10-50.99999X7脱岗5×10-20.95X8未监视10-20.99X9安装位置不当10-30.999X10表盘直径小10-40.9999X11光线不足10-30.999X12未标红线显示5×10-20.95X13超期未校10-60.999999X14没定期冲洗10-30.999X15三通关闭10-40.9999表3-16基本事件发生概率取值表2事故树定性分析(1)最小割集与最小径集①最小割集与最小径集的数目公式式中i-门的编号或代码;λi-第i个门输入事件的数量;Xi,j-第i个门的第j个输入变量(j=1,2…λi),当输入变量是基本事件时,Xi,j=1;当输入变量是门k时,Xi,j=Xk;Xi=门i的变量,如果门i是紧接顶事件T的门,则Xi=XT即为矩阵CMm×n的行数m。经计算,割集为40个,径集为3个。)(X)(XXXii,2,1,ii,2i,1i,为或门为与门iXXiXiii做出锅炉超压事故树的成功树:图3-19。T'A1'X1'X2'X3'X4'A2'X5'B3'C2'C3'D1'D2'X12'X11'X10'X9'X6'X14'X15'X13'X7'X8'图3-19图3-18的成功树②求取最小径集写出成功树的结构式,并化简,求取其最小割集:结构函数为:T’=A1’+A2’=X1’X2’X3’X4’+(X5’+B3’)=X1’X2’X3’X4’+(X5’+C2’C3’)=X1’X2’X3’X4’+[X5’+(D1’X6’D2’)(X7’X8’)]=X1’X2’X3’X4’+[X5’+(X9’X10’X11’x12’)X6’(X13’X14’x15’)(X7’X8’)]=X1’X2’X3’X4’+X5’+X6’X7’X8’X9’X10’X11’x12’X13’X14’x15’从而得到事故树的最小径集为:151413121110987634321251,,,,,,,,,,,,,,xxxxxxxxxxPxxxxPxP(2)结构重要度分析由于3个最小径集中均不含共同元素,所以得到:IΦ(5)IΦ(1)=IΦ(2)=IΦ(3)=IΦ(4)IΦ(6)=IΦ(7)=IΦ(8)=IΦ(9)=IΦ(10)=IΦ(11)=IΦ(12)=IΦ(13)=IΦ(14)=IΦ(15)3.事故树定量分析(1)顶上事件发生概率根据公式(课本式3-15),(课本式3-16),得到P(T)=[1-(1-q5)][1-(1-q1)(1-q2)(1-q3)(1-q4)][1-(1-q6)(1-q7)(1-q8)(1-q9)(1-q10)(1-q11)(1-q12)(1-q13)(1-q14)(1-q15)]将表3-16所列个基本事件发生概率值代入,得:P(T)=6.61X10-5(2)概率重要度分析与临界重要度分析——选作4.结论总结分析结果,提出预防措施。ripxiprqTP1)(ripxiprqTP11)(二、木工平刨伤手事故分析(冯肇瑞书第167~173页)1.事故树事故树:图6-58。图6-58木工平刨伤手事故树分析图X11X10X9X8X7X6X5X4X3X2X1EBD1D2AT两人配合推接不好料掌握不稳动作不协调用力过大脚下滑绊加工小薄木料料有硬节振动脱手手的动作异常身体失去平衡误动作使手伸向刀口直接用手推加工木料开机用手清木屑或异物开机量尺寸开机修理加工时手指接近旋转刀口手异常接近旋转刀口安全装置故障失灵木工平刨伤手C2.事故树定性分析(1)最小割集与最小径集①最小割集与最小径集的数目公式式中i-门的编号或代码;λi-第i个门输入事件的数量;Xi,j-第i个门的第j个输入变量(j=1,2…λi),当输入变量是基本事件时,Xi,j=1;当输入变量是门k时,Xi,j=Xk;Xi=门i的变量,如果门i是紧接顶事件T的门,则Xi=XT即为矩阵CMm×n的行数m。经计算,割集为9个;径集为3个。)(X)(XXXii,2,1,ii,2i,1i,为或门为与门iXXiXiii②求取最小径集做出原事故树的成功树:图6-59。写出成功树的结构式,并化简,求取其最小割集:T’=A1’+X11’=B1’X8’X9’X10’+X11’=(C’+X1’)X8’X9’X10’+X11’=(C’+X1’)X8’X9’X10’+X11’=……=X1’X8’X9’X10’+X2’X3’x4’X5’X6’X7’X8’X9’X10’+X11’从而得到事故树的最小径集为:11310987654322109811,,,,,,,,,,,,,xPxxxxxxxxxPxxxxPX6'X7'E'X1'X9'X10'X8'X11'T'A'B'C0'D1'D2'X5'X4'X3'X2'图6-59木工平刨伤手事故树成功树(2)结构重要度分析IΦ(11)IΦ(8)=IΦ(9)=IΦ(10)IΦ(1)IΦ(2)=IΦ(3)=IΦ(4)=IΦ(5)=IΦ(6)=IΦ(7)3.事故树定量分析(1)基本事件发生概率估计值为了计算,最重要的是确定故障率数据。而现在只能凭经验估计。从理论上讲,事故发生概率应为任—瞬间发生的可能性,是一无量纲值。但从工程实践出发,许多文献皆采用计算频率的办法代替概率的计算,即计算单位时间事故发生的次数。表6—14中的数据是从这一点出发给出的。表6-14基本事件发生概率估计值代号基本事件发生概率qi(1/小时)X1X2X3X4X5X6X7X8X9X10X11直接用手推加工木料脚下滑绊料掌握不稳二人配合推接不好用力过大料有硬节振动脱手加工小薄木料开机修理开机量尺寸开机用手清木屑或异物安全装置故障失灵0.15×10-35×10-210-410-310-510-22.5×10-610-510-34×10-4(2)顶上事件发生概率g=0.000003009/h(3)概率重要度分析与临界重要度分析——选作4.结论总结分析结果,提出预防措施。作业(课本中思考题)A.课本第76页题1。B.课本第76页题3。C.课本第77页题4(增加内容):(1)利用状态枚举法计算图3-22事故树中顶上事件的发生概率;(2)求各基本事件的结构重要系数,进行结构重要度分析。D.课本第77页题5(减少内容):已知图3-23所示事故树,试用最小割集法、最小径集法分别计算顶上事件的发生概率。