《安全事故现状与趋势分析方法研的究毕业论文.doc》由会员分享，可在线阅读，更多相关《安全事故现状与趋势分析方法研的究毕业论文.doc（16页珍藏版）》请在三一文库上搜索。
1、安徽审计职业学院毕业论文（设计）安全事故现状与趋势分析方法研究摘 要：包括安全生产事故预防、控制在内的安全生产监督管理工作需要客观把握安全生产的现状并准确地判断其总体发展趋势。本文基于我国安全生产事故快报数据，对近年安全生产现状做出季节性等分析，并为短期安全生产形势变化的预测与判断不仅提出了三个可行的ARIMA-BP、ARIMA-RBF以及ARIMA-GRNN非线性组合模型而且进一步基于RBF，对前述三个模型再次进行非线性组合，给出了一种新的双重非线性组合趋势分析方法。实证结果表明，双重非线性组合能够较为精确地预测安全生产事故的发展趋势，可以为安全生产事故的预防、控制和应对提供管理和决策支持。
2、关键词：安全生产；ARIMA模型；BP模型；RBF模型；GRNN模型The Analysis and Forecast of National Accident Trends in Work SafetyAbstract: Based on the data query system of accidents in work provided by the State Administration of Work Safety, not only does this paper analyze the statistics of work safety accidents but also e
3、xplore the development trend of work safety accidents from 2001 to 2008. (accidents in work safety are classified by the number of deaths).At the same time, the ARIMA-BP ，ARIMA-RBF and ARIMA-GRNN nonlinear combining forecast model are constructed. At last, this paper constructs a double nonlinear co
4、mbining forecast model by combining three single nonlinear combining forecast models with RBF network. The results indicate that the accuracy of the double nonlinear combining forecasting models are higher than the accuracy of the single ARIMA model and three single nonlinear combining forecast mode
5、ls. This study provides a new method for promoting the development and improvement of strategic planning of work safety. Key words: Work safety；ARIMA model；BP model；RBF model；GRNN model1 引言目前，安全生产仍然是世界各国都必须面对和应对的重大问题。对此，我国坚持“安全第一，预防为主，综合治理”的方针。显然，通过对各类安全生产事故的现状和特征分析有助于认识安全生产事故的规律，通过对安全生产事故发展趋势的准确判断有
6、助于对安全生产事故的预防、控制和应对做出有效的决策。国内外关于安全生产事故研究的报告1和文献很多，主要集中在以下方面:(1)某一类安全事故的预测预警，如李书全和窦艳杰(2008) 2针对建筑安全事故，提出了RSSVM预测模型；钟铭(2009) 3采用模糊数据融合算法，对煤矿瓦斯进行预警；(2)安全事故发生的规律分析，主要包括：安全事故的成因分析，安全事故与相关经济指标之间的关系。如吕海燕(2004)4对我国工矿企业、消防火灾、交通等进行总量统计分析，提出了相应的事故统计分析指标体系框架；Kopits和Cropper(2005)5专门对交通死亡事故与经济增长的关系作了深入研究；任荣森(2009)
7、 6等人建立了油库的多层递阶解释结构模型(ISM)，将油库安全事故成因分为根本原因、直接原因和间接原因三个层次；(3)对安全生产管理的研究，主要包括：对事故管理绩效的评估、城市生产安全事故应急管理能力评价、安全事故损失的计算、安全事故中的风险管理。如牛跃林等人(2008)7建立了城市生产安全事故应急管理能力评价指标体系，为城市的生产安全事故应急体系建设提供决策依据；Karen(2009) 8和 Shi(2009) 9 对煤矿安全事故风险管理及造成的损失从政府、行业规模和运营等方面作了探讨。上述已有研究主要围绕特定行业领域展开，而且主要集中在煤矿、建筑工程、交通安全生产事故。固然，安全事故规律需
8、要分不同行业领域进行研究，但也需要整体认识和把握，这样才能制定更科学更完善的国家安全生产战略规划。基于此，本文重点针对安全生产事故的起数及其造成的死亡人数开展方法建模和实证研究。首先把安全生产事故按其造成的死亡人数进行分类，具体地将安全生产事故分成10类：死亡个人的事故为第类(=1,2,9)，死亡10人及以上的事故为第10类。这样一种分类方式符合国家安全生产按不同伤亡和损失规模做整体统计的原则。因此，相应的分析结果有助于国家安全生产事故预防、控制和应对的监管工作。显然，对不同类型的事故所采用的处置和应对方式以及相应的资源配置也不同。本文对从2001年1月到2008年12月国家安全生产快报包含的
9、1万7千6百多条数据进行了分类型、分时段、分区域的整理分析和比较,特别结合各类安全生产月平均事故量序列的季节指数图，指出了需要重点关注的月份。由于安全生产管理决策离不开预测，而预测的精度又决定了预测结果的可用性和有效性。所以基于时间序列，建立了10类安全生产事故的数量以及安全生产事故总量的ARIMA-BP、 ARIMA-RBF以及ARIMA-GRNN非线性组合预测模型，并结合RBF，构建了基于前述三个模型的双重非线性组合预测方法，提高了预测精度。2 安全生产事故现状分析根据国家安全生产监督管理总局官方网站发布的安全生产事故快报数据，对2001-2008年的安全生产状况进行了分类型、分时段、分区
10、域的整理分析和比较,并结合各类安全生产月平均事故量序列季节指数图，指出了需重点关注的月份。2.1 年度状况为了分析全国安全生产事故近年发展趋势及其变化，对2001-2008年间的安全事故进行了整理和10个类别的分类统计。整体年度状况见表21。基于表21，把2001-2008年各类安全生产事故的发生变化趋势制成如下的图21。表2-1 2001-2008年各类安全生产事故总量统计 事故起数 年份死亡1人死亡2人死亡3人死亡4人死亡5人死亡6人死亡7人死亡8人死亡9人死亡10人及以上合计2001年53311120912398523733259013112002年1445230412200121874
11、0452710927162003年17932655472271418151463211332962004年13821976402841599659341811729862005年13391697202961829745342212630302006年134227312581521195549199416332007年34762283142915546327814962008年4456721512366552925881176合计66331002458818861118689397316200815图2-1 2001-2008年各类安全生产事故每年事故量趋势变化图2.2 月度状况采用季节指数对2
12、001-2008年的安全生产事故进行月度趋势分析。假设观测值序列的数据结构以期为一个周期，共有个周期，本文中=12,=8。用表示第类生产事故年第月的事故起数,则第类事故在第月的季节指数 。 如果季节指数值大于1，那么第类事故该月的事故平均值高于其总的事故平均值；如果季节指数值小于1，那么第类事故该月的事故平均值低于其总的事故平均值。各类事故的季节指数图汇总在下面的图22。15图2-2 2001-2008年各类安全事故月平均事故量序列季节指数图从季节指数上看，第1类至第5类事故每年的8月份是事故高发期，需要加强监管。特别地，第1类和第5类每年的7月份也需要高度关注。第6类至第10类事故呈现与前面
13、5类不同的特点，第6类需要加强监管的月份是每年的2月、5月和10月；第7类需要加强监管的月份是每年的1月、2月和8月；第8类需要加强监管的月份是每年的1月、4月和7月；第9类需要加强监管的月份是每年的4月、5月和10月；第10类需要加强监管的月份是每年的3月、8月和11月。2.3 区域安全生产事故状况为了掌握不同地区安全生产事故的情况，通过数据分析整理，给出了2001-2008年各类安全生产事故在全国各地区的分布情况，见表22:表2-2 我国31个省市区2001-2008年各季度安全生产事故平均总量统计表 省市区第一季度第二季度第三季度第四季度01北京市4.0006.3756.0004.750
14、02天津市0.7501.2501.3751.25003上海市1.1252.3752.1252.37504重庆市42.87544.87551.12543.50005四川省71.62574.87584.87583.25006黑龙江省13.87515.00025.75019.12507吉林省12.87515.75018.25018.75008辽宁省9.75014.00020.62514.37509河北省16.50025.50020.5019.62510山西省11.62514.25017.12514.50011山东省7.87511.37512.87512.75012陕西省24.37526.87532.
15、87529.50013河南省10.12513.75014.5007.62514湖北省10.50010.1258.3757.62515湖南省41.87548.50061.00047.00016江西省4.5008.0009.5006.37517安徽省19.25021.87521.75021.62518江苏省6.6257.3755.8757.87519浙江省12.75012.00010.50017.00020福建省21.00019.87523.12517.50021广东省22.75022.12516.00015.75022海南省2.37502.0001.6251.87523贵州省36.62542.3
16、7559.75049.87524云南省24.25025.75022.00021.00025甘肃省15.25018.50021.87521.00026青海省2.2502.0001.3751.87527广西壮族自治区14.8758.87514.37518.12528内蒙古自治区7.8758.87510.759.37529宁夏回族自治区2.8755.5004.0005.37530新疆维吾尔自治区5.8756.87513.00010.87531西藏自治区4.6253.6258.1257.1253 各类安全生产事故总量预测第2节从不同角度分析和考察了2001-2008年的安全生产事故现状，本节基于时间序
17、列模型及神经网络，通过三阶段方式给出了双重非线性组合预测方法。 安全生产事故量随着时间的变化而变化，由于ARIMA模型在各类安全生产事故量预测过程中既考虑了这些数量在时间序列上的相关性，又考虑了随机变化的干扰性，对短期趋势预测具有较高的准确率。所以，第一阶段对各类安全生产事故量采用ARIMA建模。因为安全生产事故涉及多种复杂的影响因素，应用单一的预测方法难于满足预测的有效性，所以有必要采用多个预测方法进行组合预测，综合利用各种预测方法提供的信息，以便减少预测误差而获得更好的预测结果。Bates10和Clemen 11都曾提出，可以将各种单一的预测结果进行组合得到一种组合预测结果，以达到改善预测
18、效果的目的。20世纪90年代初，Granger 和Terasvirta12等人研究发现，借助非线性函数关系的非线性组合预测比线性组合预测所得的结果误差要小。又由于非线性函数关系确定非常困难，而用神经网络来实现这种非线性函数映射相对容易并且神经网络具有自组织、自学习能力、非线性逼近能力强的特点。所以，第二阶段在第一阶段结论的基础上应用神经网络技术构建了三种可行的安全生产事故量非线性组合预测模型。最后，第三阶段在第二阶段基础上给出了双重非线性组合预测方法。对以下各阶段得出的预测模型，我们都利用2001年到2009年9月的安全生产事故快报数据作为已知信息，对2009年10月和11月的情况做出预测，并
19、与实际值做出比较以验证模型的精度。3.1 第一阶段预测3.1.1 ARIMA模型13ARIMA模型又叫BoxJenkins法，是以美国统计学家George E.P.Box和英国统计学家Gwilym M.Jenkins命名的时间序列分析方法。其结构如下： 简记为ARIMA()。当在阶差分后又作了步长为的周期差分，则此时模型结构为：。式中，为后移算子，；其中，是需要估计的自回归系数和滑动平均系数，；为零均值白噪声序列。3.1.2 ARIMA模型预测ARIMA建模过程包括数据检验、模型识别、参数估计、模型适应性诊断等。选择2001-2009年9月之间的数据拟合ARIMA模型。记（=1,2,9）代表死
20、亡个人的事故总量，代表死亡10人及以上的事故总量，代表总的事故量。(=1,2,9)代表第期死亡个人的事故总量，代表第期死亡10人及以上的事故总量。代表第期总的事故量。得到如下11个预测模型。 (1A)(1+0.41189-0.17536+0.39474+0.24147-0.20759)=， (0,477.19)；(2A) (1+0.54802+0.21288+0.18633+0.25675+0.20701)=， (0,14.84)；(3A) (1 +0.56396)=(1-0.38228)， (0,108.02)；(4A) (1 - 0.26609- 0.33063)(-19.06625)=
21、， (0,39.22)；(5A)=， (0,19.68)；(6A)=， (0,13.62)；(7A)=，(0,10.28)；(8A) (1+0.81153+0.61162+0.51758+0.20181)=， (0,3.65)；(9A)=， (0,3.93)；(10A) (1 +0.30446)=(1-0.79891)， (0,12.49) ；(11A) (1-0.4416-0.2876+0.18249-0.21841)=(1-0.77191)， (0,937.366)根据上述模型，得到2009年10-11月份的各类安全生产事故量的点预测值，见表35。同时也给出点预测值的置信水平为95%的置信
22、限以及安全生产事故总量时间序列拟合图，见表36和图31。 其中置信限计算公式如下：由3.1.1知，记，又令，根据的同幂次系数相等的原理，可以确定，即其中，。 在均方误差最小原则下，容易求出,，所以， ，的置信水平为的置信限为：其中，为标准正态分布的分位数，本文取5%。（注：为时期的预测误差）3.1.3 安全生产事故总量ARIMA拟合预测图以下是安全生产事故总量ARIMA拟合图。图中星号为实际值，中间实曲线为依据3.1.2中(11A)得到的预测值，上下两条虚曲线为预测值置信水平为95%的置信上限和置信下限。图31 安全生产事故总量时间序列拟合预测图3.2 第二阶段预测3.2.1 BP模型14BP
23、神经网络由Rumelhart等人于1986年建立，是一种多层前馈型神经网络，由输入层、隐层和输出层构成。隐层神经元的传递函数通常用log-sigmoid型函数、tan-sigmoid型函数、purelin纯线性型函数。输出层的函数常用purelin纯线性型函数。它是前向网络的核心部分，体现了人工神经网络的精华。本文采用BP建立了第一种非线性组合预测方法，记作：ARIMA-BP。BP模型表达式如下：,其中为节点输出；为节点连接权值；为节点输入；为神经元阈值；为输入/输出的传递函数。3.2.2 ARIMA-BP模型预测 依据3.1.2的11个模型，对10类事故以及事故总量进行第二阶段建模。ARIM
24、A-BP模型的实现过程由模型输入、网络训练、网络测试、模型输出等四个部分组成。确定神经网络为3层，输入层、隐层及输出层，隐层采用S型函数作为传递函数，输出层取线性函数。网络训练采用基于数值最优化理论的Levenberg-Marquardt15-16方法进行，然后调用这些训练好的权值与阈值，经输出层输出即可得到预测值。依据3.1.2，确定了安全生产事故量随时间变化的规律，因此构建的输入输出样本对以及网络结构各层神经元个数见下面的表31： 表31 ARIMA-BP模型结构表事故类型网络输入期望输出网络结构11-10-110-10-13-6-14-5-15-5-111-10-114-12-16-6-
25、19-10-15-5-111-10-1应用ARIMA-BP模型关于2009年10月和11月事故起数的预测结果见表35。3.2.3 RBF模型17-18径向基函数（RBF）是多维空间插值的传统技术，由Powell于1985年提出。1988年，Broomhead和Lowe根据生物神经元具有局部响应这一特点，将RBF引入神经网络设计中，产生了RBF神经网络。它属于前馈型神经网络，是一种三层的前向网络。第一层为输入层；第二层为隐层（径向基层）；第三层为线性输出层。1989年，Jackson论证了RBF神经网络对非线性函数的一致逼近性能。因而将它应用于复杂的时间序列预测会取得较好的效果。本文采用RBF建
26、立了第二种非线性组合预测方法，记作：ARIMA-RBF。RBF模型表达式如下：径向基层输入为，径向基层输出为，线性输出层输出为。其中，为输入向量；和分别为径向基层和线性输出层的节点连接权值；和分别为径向基层和线性输出层的神经元阈值。 3.2.4 ARIMA-RBF模型预测ARIMA-RBF模型构建的输入输出样本对与ARIMA-BP完全相同。隐层的径向基函数设计为高斯函数，输入为输入向量与权值向量的距离乘以阈值，输出即为神经元的输出：。其中为输入向量；为径向基层的节点连接权值；为径向基层的神经元阈值。另外，还需要给出径向基函数的散布常数（SPREAD）。本文综合以下三个原则确定：(1)通常，散布
27、常数的选取取决于输入向量之间的距离，要求是大于最小距离，小于最大距离；(2)网络的预测误差和逼近误差； (3)散布常数等于1作为参照。首先由输入向量间距离矩阵确定SPREAD的范围，然后在此范围内结合预测误差、逼近误差以及参考散布常数1决定最终的SPREAD值。ARIMA-RBF模型结构见下面的表32：表32 ARIMA-RBF模型结构表事故类型网络输入期望输出SPREAD0.100.100.040.151.801.501.500.701.700.502.00应用ARIMA-RBF模型关于2009年10月和11月事故起数的预测结果见表35。3.2.5 GRNN模型19-20 广义回归神经网络（
28、GRNN）是美国学者Donald F.Specht在1991年提出的。在结构上与RBF网络较为相似。它由输入层、模式层、求和层和输出层四层构成。GRNN也具有很强的非线性映射能力，适用于解决非线性问题。GRNN网络最后收敛于样本量积聚较多的优化回归面，并且在样本数据较少时，预测效果也比较好。另外它还可以处理不稳定的数据。本文采用GRNN建立了第三种非线性组合预测方法，记作：ARIMA-GRNN。GRNN模型表达式如下：，其中是在输入向量为的条件下，输出向量的期望值。是和的联合概率密度函数，可由Parzen非参数估计出来。3.2.6 ARIMA-GRNN模型预测ARIMA-GRNN模型构建的输入
29、输出样本对也与ARIMA-BP相同。通过GRNN模型不断学习归纳出输入和输出变量之间的关系。类似的，GRNN的SPREAD值的确定原则与3.2.4一样，ARIMA-GRNN模型结构见下面的表33：表33 ARIMA-GRNN模型结构表事故类型网络输入期望输出SPREAD0.100.100.06 0.151.801.501.500.701.700.501.00应用ARIMA-GRNN模型关于2009年10月和11月事故起数的预测结果见表35。3.3第三阶段预测3.3.1 双重非线性组合模型的构建将第二阶段三种预测方法关于训练样本的预测结果作为RBF的输入向量，输入节点数为3，而实际值作为期望输出
30、向量训练网络，输出节点数为1。SPREAD值的确定原则与3.2.4一样。最后，把上述三种预测方法关于2009年10月和11月事故起数的预测值作为输入向量来预测10月和11月的事故起数。上述基于RBF的ARIMA-BP、ARIMA-RBF以及ARIMA-GRNN的双重非线性组合模型记作： 3.3.2 双重非线性组合模型预测由3.3.1得到的双重非线性组合模型结构见下面的表34：表34 模型结构表事故类型网络输入期望输出SPREAD0.0110.1200.3300.7901.0000.0401.0000.1001.0001.0001.000应用双重非线性组合模型关于2009年10月和11月事故起数
31、的预测结果见表35。3.4模型预测效果的比较我们把各种模型对2009年10-11月各类安全生产事故的数量和安全生产事故总量的预测结果总结在下面的表35。同时用表9的误差评价表，即平均绝对误差(MAE)和均方根误差(RMSE)来比较各种模型。表35 各模型预测结果与实际值对比表事故类型月份ARIMAARIMA-BPARIMA-RBFARIMA-GRNN双重非线性组合实际值X1Oct-090.34580.23631.00001.00001.00001Nov-090.59811.18461.00001.00001.00001X2Oct-090.54011.26311.00001.00001.8623
32、1Nov-090.00000.74731.00001.00002.29603X3Oct-0958.166970.548262.383161.684061.898468Nov-0959.568360.909260.375145.000048.362148X4Oct-0919.300720.590218.720714.999319.041619Nov-0923.065820.227920.253320.291320.253321X5Oct-0910.670414.775512.757111.106813.526617Nov-0914.32788.63217.273012.04817.49497X6
33、Oct-096.61297.89515.18086.21226.00006Nov-097.60945.99254.69836.45396.00006X7Oct-092.27731.90032.28363.83342.88892Nov-095.00133.04372.55963.72572.55052X8Oct-092.00342.47682.66782.90143.00233Nov-092.24423.62872.90893.01993.00353X9Oct-090.63592.01661.14371.73111.15721Nov-091.04353.91982.76781.90442.775
34、23X10Oct-095.58899.16117.71479.90217.88158Nov-095.48907.48498.041810.80108.37026X11Oct-09131.7766130.4441133.0933116.1790132.3547132Nov-09118.5314115.6359113.3715110.2932113.3678116表36 2009年10-11月ARIMA预测值置信水平为95%的置信区间置信限X1X2X3X4X5X6X7X8X9X10X11Oct-09(0,43.1606)(0,8.0898)(37.7969,78.5369)(7.0259,31.5
35、755)(1.9745,19.3663)(0,13.8451)(0,8.5616)(0,5.7497)(0,4.5212)(0,12.5155)(71.7695,191.7837)Nov-09(0,50.2683)(0,8.0448)(37.3461,81.7906)(7.7273,33.1312)(4.9910,23.6646)(0.3773,14.8416)(0,11.2856)(0,6.0564)(0,5.1012)(0,12.4156)(46.3108,190.7520)事故类型评价标准ARIMAARIMA-BPARIMA-RBFARIMA-GRNN双重非线性组合X1MAE0.5281
36、0.47420.00000.00000.0000RMSE0.54300.55560.00000.00000.0000X2MAE1.73001.25791.00001.00000.7832RMSE2.14611.60371.41421.41420.7871X3MAE10.70077.72878.99604.65803.2319RMSE10.73589.30439.60974.94434.3221X4MAE1.18331.18150.51302.35470.3942RMSE1.47611.25000.56372.87300.5288X5MAE6.82871.92832.25805.47061.98
37、42RMSE6.84691.95093.00645.48692.4809X6MAE1.11110.95131.06050.33310.0000RMSE1.21771.34001.08760.35430.0000X7MAE1.63930.57170.71161.77960.7197RMSE2.13130.74140.72771.78040.7393X8MAE0.87620.57600.21170.05930.0029RMSE0.88440.57840.24350.07140.0028X9MAE1.16030.96820.18800.91340.1910RMSE1.40720.96940.1931
38、0.93130.1941X10MAE1.46111.32301.16363.35161.2444RMSE1.74281.33291.45783.65151.6781X11MAE1.37740.96021.860910.76391.4935RMSE1.79691.12992.846811.89211.8781总体MAE2.59971.62921.63302.78950.9132总体RMSE2.81171.88701.92283.03631.1465表37 各模型误差评价表从表35和表37看，ARIMA-BP 、ARIMA-RBF 、ARIMA-GRNN和均优于ARIMA。三个单一的非线性组合预测
39、模型中，最优的是ARIMA-RBF，其次是ARIMA-BP，最后是ARIMA-GRNN。而本文最后给出的不仅优于单一的ARIMA预测模型而且也优于三个单一的非线性组合预测模型。最小MAE为0，最大MAE为3.2319。最小RMSE为0，最大RMSE为4.3221。的最大MAE分别仅为 ARIMA、ARIMA-BP、ARIMA-RBF、ARIMA-GRNN最大MAE的30.2%、41.82%、35.93%和30.03%。的最大RMSE分别仅为 ARIMA、ARIMA-BP 、ARIMA-RBF、ARIMA-GRNN最大RMSE的40.26%、46.45%、44.98%和36.34%。由此可见，双
40、重非线性组合模型不仅提高了预测精度而且还能满足非线性的预测要求，这为解决因安全生产事故量随机变动而造成的预测困难及精度问题提供了新的思路和方法，是安全生产中可行的预测模型。4 结语本文按生产事故造成伤亡和损失规模的原则重新对安全生产事故进行分类，并对各类安全生产事故进行了现状分析，借助季节指数分析法对各类安全生产事故指出了应重点预防的时段建议。最后为了精确地预测安全生产事故的发展趋势，对各类安全生产事故量以及安全生产事故总量先建立了三个可行的非线性组合预测模型，即ARIMA-BP、 ARIMA-RBF以及ARIMA-GRNN，然后在此基础上构建了一种新的双重非线性组合趋势分析方法。为提高预测的
41、准确性提供了新的方法和思路，为管理和决策部门制定安全生产战略规划提供了理论依据。通过本文建立的模型及其实证研究结果，一定程度上说明就我国目前所处的经济发展阶段而言，各类安全生产事故发生的最严重程度及经过努力可以实现的改进程度是可以做出统计分析和预测的。显然，各类事故发生的频度也有很大差别，进一步弄清其原因是今后值得特别关注的问题。各类事故在一天24小时内的不同时间段发生的频度同样值得进一步关注，所有这些对控制和改进安全生产状况都具有重要意义。参考文献:1 王显政. 安全生产与经济社会发展报告R.北京：煤炭工业出版社,2006.2 李书全,窦艳杰. 基于RSSVM模型的建筑安全事故预测模型J.统
42、计与决策, 2008, (19)：56-58.3 钟铭.基于模糊数据融合的煤矿瓦斯预警研究D.西安科技大学博士论文,2009.4 吕海燕.生产安全事故统计分析及预测理论方法研究D.北京林业大学博士论文,2004. 5 Kopits,E., Cropper,M.Traffic Fatalities and Economic GrowthJ.Accident Analysis and Prevention,2005,37(1) ：169-178.6 任荣森,雍歧东,张宁. 基于ISM的油库安全事故成因分析J.后勤工程学院学报,2009,25(5)：58-61. 7 牛跃林,王花平,谢金亮.基于模糊数学的城市生产安全事故应急管理能力综合评价J.江西理工大学学报,2008,29(1)：53-56.8 Karen,P. Blood o