本发明涉及专门适用于安全生产事故因素分析技术领域，特别是系统性分析多元事故因素层次间关系的一种方法。
背景技术：
安全生产事故通常不是由某个单一因素引起的，而是由一连串的多元因素综合造成的结果。同时，引起事故中的一个因素可能会影响另一个因素，这些因素之间可能互相会产生影响关系。在以往的安全生产处理中，仅仅强调单一或抽取主要因素对整个事件的影响，而忽略了事故中多个因素相互间的层次关系，会造成分析结果不全面、客观，安全事故仍旧发生的现象。在过去几十年中，对安全事故的分析过程中，已有系统的分析方法，例如人因分析与分类系统、功能性核磁共振分析法，这些方法在对事故发生的因素进行了有效地分析，但它们只是在设定好的框架中进行分析。目前对于安全事故常采用系统建模的方法对其进行分析，将事故分析的研究方向转向揭示系统的内在结构，使用较多的模型为解释结构模型(ism)，该方法利用数学方法推导出了结构中元素中客观的层级制并给出了两两间元素的关系，为系统工程中结构建模的流行算法，常被用于建立风险因素图以及专业问题的属性层次结构、评估模式。但该模型中的邻接矩阵依赖于经验判断，在对特殊的事件因素系统分析并不准确。因此须建立一种系统的、基于数据和客观计算的事故因素结构分析方法，减少人为干预的因素，有利于更加客观、准确地对安全事故进行系统分析，避免类似事故发生同时为安全生产部门提供科学的决策依据。
技术实现要素：
本发明的目的是为了克服目前对安全生产事故因素分析中仅凭借经验判断设立事故因素关联矩阵、事故因素层次分析不全面、不具体的缺点，提供了一种利用灰色关联确定事故因素层次结构的系统方法。
通过灰色关联分析，事故因素之间的关系根据执法检查清单记录中产生的时间序列来计算。通过矩阵运算，这些因子的灰色可达矩阵是由灰色关系组成的灰色邻接矩阵导出的。可达矩阵通过代数操作被划分成不同的层次，并且这些因子被安排在一个层次结构中。由此产生的层次结构提供了一个事故因素的整体情景，有助于有效追踪落后和发现相关事故因素的抽查情况，为进行综合治理并消除事故因素提供理论依据，为有效规避类似安全事故的发生提供可靠保障。
为了实现上述目的，本发明采用的分析方法可分为五个步骤：
步骤一、计算整体事故因素之间的灰色关联度；
步骤二、构建灰色邻接矩阵；
步骤三、建立可达性矩阵；
步骤四、将可达性矩阵划分为不同的级别；
步骤五、绘制分层结构图。
本发明的优点在于：
1.本发明建立了一套定量化方法，可以根据行为序列自主导出事故因素的图形层次结构，有利于更加客观、准确的分析事故因素层次间的关系，有利于减少安全事故的发生。
2.本发明将灰色关联方法应用于事故因素层次间的分析，建立事故因素间的灰色邻接矩阵以及可达化矩阵，为事故因素不同级别的划分提供了科学客观的依据。
3.本发明实现了为政府相关安全生产部门决策提供客观的参考和量化依据，依据执法部门检查清单数据计算分析结果，尽量减少人为因素以及经验判断对事故分析产生的影响，使得分析结果更贴近于现实情况。
附图说明
图1是本发明的整体技术路线示意图；
图2是元素可达性分析示意图；
具体实施方式
下面结合附图和实施例对本发明进行详细说明。
本发明提供的一种利用灰色关联确定安全事故因素层次结构的系统方法，系统框图如图1所示，具体包括步骤如下：
步骤一、计算整体事故因素之间的灰色关联度。通过从执法部门处取得的执法检查清单中获得安全事故因素并通过观察得到所有元素的行为序列，执法清单主要用作事故因素的数据收集，事故因素的统计是执法检查记录的总结，即不合格项目的频率。得到相关的安全事故因素时间序列后，首先要对这些数据进行归一化处理，其次，计算这些安全事故因素之间的灰色关联系数。得到了灰色关联系数后利用公式计算出安全事故因素间的灰色关联度。
步骤二、构建灰色邻接矩阵。通过计算事故因素之间的灰色关联度，事故因素两两之间的关系得以被量化，为了确定哪些因素与其他因素之间存在关系，将使用灰色邻接矩阵来表示。灰色关联度可以表明哪个相关行为序列更有可能导致另一个事故因素中的行为序列，选定安全事故因素中某个行为序列，参考对比其他的行为序列，从而选定的这个行为序列与其他因素之间的灰色关联可以被获得。
步骤三、建立可达性矩阵。可达性矩阵用于表示有向图中不同节点(顶点)通过某种特定路径可以达到的程度。要考虑的是传递性，表明如果安全事故因素与因素之间的关系。在安全事故的系统性分析中，传递性规则应该是真实的，元素ei与ej有关，且ej与ek有关，则ei与ek有关。即事故中a因素影响着b因素，b因素影响着c因素，那么a因素也影响了c因素，如果这个传递性条件满足且eirej≠ejrei，即满足不对称性，那么这些条件适用的矩阵将始终是一个层次结构。分析示意图如图2。
步骤四、将可达性矩阵划分为不同的级别。得到了事故因素之间的可达性矩阵后，可以定义级别。每一个因素ei都有一个可达性集合r(ei)，这个集合由所有ei路径上的来源于ei的元素和包含ei但不来源于ei的元素构成。根据定义，元素ei包含在两个集合中。元素ej将在ei的可达集中，当且仅当rij＝1。一个元素ej在ei的先行集中，当且仅当rji＝1。给定了可达性矩阵，矩阵中满足eir(ei)∩a(ei)＝r(ei)的元素被定义为有向图顶层的部分。
步骤五、绘制分层结构图。根据将可达矩阵进行分级的方式，可以确定结构图的所有层级，并在矩阵的参照下，可以通过从上到下放置不同的元素来绘制多级分层结构有向图。
上述五个步骤可以在计算机程序中实现。通过计算机的辅助，该方法可以将事故的因素行为序列自动转换为二元图。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，均同理包括在本发明的专利保护范围内。
实施例一
该实例实现对确定事故因素层次的技术包含四个部分：
(1)计算事故因素中的关联度。假设事故系统中有n个元素，集合为s＝{e1,e2,...,en}，通过对应观察执法清单中检查的项目，可以得到所有元素的行为序列，在这里是时间序列。x是这些时间序列的集合：x＝{xi|i＝0,1,2,...,n}。xi,yi∈x，符号γ(xi,yi)表示xi和yi之间的邻接度，即所谓的灰色关联系数。灰色关联的定理为：
1)可比。即x中的每一个序列都是无维度、等规模、均匀极性，且数量大于3。
2)规范间隔。0＜γ(x0,yi)≤1。(i)如果xi＝xj，γ(xi,yi)＝1。(ii)如果xi,yj∈φ，γ(xi,yi)＝0。
3)整体性。如果n≥3和i≠j，则γ(xi,yj)≠γ(yj,xi)
4)可接近。
只要满足定理1)-4)即可计算灰色关联度。关联系数γ(x0(k),xi(k))是通过领域本质和正常区间性质将灰色关联定理化计算得到。整体事故因素的灰色关联度的计算公式为：
(2)灰色关联矩阵与可达矩阵的构建
令x1为特征行为序列，可以得到e1与其他元素ei之间的灰色关联度，这将构成灰色关联矩阵的第一列。同理，令x2为特征行为序列，可以得到e2与其他元素ei之间的灰色关联度，这将构成灰色关联矩阵的第二列。重复此操作指导获得第n列。灰色可达矩阵可以通过将灰色邻接矩阵和单位矩阵加在一起并且多次应用操作获得。
(3)将可达矩阵划分为不同的级别
一旦确定了矩阵的顶层集合，顶层集合中的元素和与它们的连接的所有表示方式都可以从矩阵中移除，留下的子矩阵也将具有顶层集合。子矩阵的顶层集合将是矩阵的第二层集合。
(4)绘制分层结构图
得到了可达矩阵的不同级别划分，可以将元素从上到下置于层次结构图中，得到了安全事故的多级结构图。多级结构图的表层因素是事故直接因素，次深层因素是事故过渡因素，深层是事故根本因素。