CN105893758A - 一种wmnt瓦斯爆炸灾害风险评价方法 - Google Patents

Abstract

本发明公开了一种WMNT瓦斯爆炸灾害风险评价方法，属于煤矿安全技术领域。本发明方法包括如下步骤：一、确定瓦斯爆炸灾害风险指标体系，对瓦斯爆炸灾害风险划分风险等级；二、确定所述指标体系中的各指标的权重向量；三、确定成本型隶属度函数或效益型隶属度函数，并得到指标隶属度矩阵；四、计算所述风险等级的目标隶属度向量；五、判断瓦斯爆炸灾害的风险等级。本发明提供的方法计算简单、操作性强、评价结果可靠，是一种非线性评价瓦斯爆炸灾害风险等级的模型。

Description

一种WMNT瓦斯爆炸灾害风险评价方法

技术领域

本发明涉及一种用于煤矿安全生产中对于瓦斯爆炸灾害评价的方法，特别是涉及一种瓦斯爆炸灾害风险评价方法，属于煤矿安全生产技术领域。

背景技术

煤矿生产系统是一个典型的非线性系统，具有复杂性、动态性和模糊性等特点。当前，煤矿安全成为制约煤炭工业发展的主要问题之一。瓦斯爆炸灾害事故的频频发生严重影响着煤矿职工的生命安全和煤炭工业的可持续发展。由于我国存在大量的高瓦斯矿井，且随着开采深度的增加及机械化程度的加强，瓦斯涌出量逐渐增大；而瓦斯爆炸危险因素具有非线性动力学特征，因此对瓦斯爆炸灾害进行风险评价是一个重要的课题。

目前，众多专家学者对瓦斯爆炸灾害风险评价方法进行了研究。如：随着现代应用数学在安全评价理论中的应用和发展，早期出现了模糊数学综合评价方法、神经网络评价方法、AHP-GT危险性评价模型、G-K评价与粗糙集模型以及属性模式识别模型。这些研究促进了煤矿安全技术的发展，但上述方法均为线性评价方法，不能很好的适应非线性煤矿生产系统的发展。

期刊“基于IAHP-ECM的瓦斯爆炸灾害风险评估”使用区间层次分析法(IAHP)建立层次结构体系，用区间数特征根法求解指标权向量，设置指标风险等级分界点，确立单项指标分段功效函数改进功效系数法(ECM)，从而构建了基于IAHP-ECM的灾害风险评估模型。该模型虽然操作性强，但在计算指标权重时采用的IAHP方法过于繁琐。期刊“基于网络分析和联系熵的煤与瓦斯突出预测研究”将网络分析法和联系熵理论相耦合，建立了煤与瓦斯突出预测的ANP-CE模型。该模型使用网络分析法建立了煤与瓦斯突出预测指标的网络模型，并且计算出指标的权重、划分危险性等级，最终由联系熵累加确定了各风险等级的联系熵范围。该模型在煤与瓦斯突出危险性预测上具有可行性和合理性，但利用联系熵的累加使得计算复杂度有所增加。

通过上述分析，目前瓦斯爆炸灾害风险评价存在的问题主要有：(1)常用的评价方法仍为线性评价方法；(2)模型计算复杂，操作性相对较弱。

发明内容

为解决上述技术问题，本发明提供了一种WMNT瓦斯爆炸灾害风险评价方法。首先，从“人-机-环境-管理”4个方面考虑，确定瓦斯灾害风险的指标体系，对瓦斯爆炸灾害风险划分风险等级；针对AHP定权方法存在的问题，本发明采用改进的AHP单准则排序方法确定所述指标体系中的各指标的权重向量(Weight,W)；本发明根据所述各指标匹配的各风险等级的临界值确定隶属度函数并得到指标隶属度矩阵(Membership Function,M)；针对当前用加权平均线性模型实现隶属度转换的不足，本发明使用非线性隶属度转换的方法和所述指标体系中各指标的权重向量计算所述风险等级的目标隶属度向量(Nonlinear Membership Transform,NT)，最终根据置信度识别准则和所述目标隶属度向量判断瓦斯爆炸灾害的风险等级。本发明为瓦斯爆炸灾害风险评价提供了新的方法。

本发明的目的通过以下技术方法来实现：一种WMNT瓦斯爆炸灾害风险评价方法，包括以下步骤：

步骤1：从“人-机-环境-管理”4个方面考虑，确定瓦斯爆炸灾害风险指标体系，对瓦斯爆炸灾害风险划分风险等级；

步骤2：采用改进的AHP单准则排序方法确定所述指标体系中的各指标的权重向量；

步骤3：根据所述各指标匹配的各风险等级的临界值，将隶属度函数分为成本型隶属度函数或效益型隶属度函数；依据各风险等级隶属度函数的确定公式来确定所述成本型隶属度函数或效益型隶属度函数，并根据所述成本型隶属度函数或效益型隶属度函数得到相应的指标隶属度矩阵μ；

步骤4：根据非线性隶属度转换算法和所述指标体系中各指标的权重向量，将所述的指标隶属度矩阵μ转换为目标隶属度向量；

步骤5：根据置信度识别准则和所述目标隶属度向量，判断瓦斯爆炸灾害的风险等级。

本发明的有益效果在于：

(1)本发明利用改进的AHP单准则排序方法确定所述指标体系中的各指标的权重向量，具有一定的理论基础。并且实现了因素层各因素的无量纲归一化，满足风险评价的分级要求；

(2)本发明采用的隶属度转换算法是非线性转换，既不增加新的先验知识，也不造成已有信息损失。该模型具有可行性和正确性，为瓦斯爆炸灾害风险评价提供了新的思路；

(3)本发明提出的WMNT瓦斯爆炸灾害风险评价模型的计算简单、操作性强、评价结果可靠，实现了非线性评价瓦斯爆炸灾害风险等级的方法。该模型在瓦斯爆炸灾害风险评价上具有科学性和准确性。

附图说明

图1是本发明的WMNT瓦斯爆炸灾害风险评价方法流程图；

图2是本发明的成本型隶属度函数示意图；

图3是本发明的效益型隶属度函数示意图。

具体实施方式

下面结合附图和具体实施方式对本发明做进一步详细的说明。

图1是本发明的一种WMNT瓦斯爆炸灾害风险评价方法流程图，该方法包括以下步骤：

步骤1：从“人-机-环境-管理”4个方面考虑，确定瓦斯爆炸灾害风险的指标体系，对瓦斯爆炸灾害风险划分风险等级；

确定的瓦斯爆炸灾害风险指标体系，如表1所示。

根据煤矿生产系统的特点、专家经验及常用《国家突发公共事件总体应急预案》中规定的分级标准，将瓦斯爆炸灾害风险等级划分为5个等级：V＝{v₁-安全，v₂-较安全，v₃-一般安全，v₄-较不安全，v₅-不安全}。

表1瓦斯爆炸灾害风险指标体系

步骤2：由于AHP给出的基于判断矩阵一致性检验的特征根定权法中，临界值的确定缺乏必要的理论基础，采用改进的AHP单准则排序方法确定所述指标体系中的各指标的权重向量；

步骤21：当所述指标体系中指标i与j进行重要性比较时，用s_ij表示所述指标i相对于所述指标j的重要性得分，用s_ji表示所述指标j相对于所述指标i的重要性得分，且s_ij，s_ji满足：

$\{\begin{array}{c}{s}_{ij}\≥0,s_{ji}\≥0\\ s_{ij}+s_{ji}=\left\{\begin{array}{cc}1& i\≠j\\ 0& i=j\end{array}\right.\end{array}---\left(1\right)$

显然，s_ij是一种标度，因为s_ij能表征所述指标体系中指标i与j重要性程度的差异，称s_ij是评分标度，同理s_ji也是评分标度。

步骤22：虽然s_ij∈[0,1]，但是s_ij不是[0,1]上的连续实数而是[0,1]区间上的离散值，通常可由比例标度a_i/a_j转换成评分标度。由比例标度a_i/a_j转换成评分标度的转换公式为：

这样，对区间[1/9,9]上任一比例标度a_i/a_j(i≠j)有唯一的评分标度s_ij∈[0.1,0.9]与之对应。

$s_{ij}=\left\{\begin{array}{cc}0& i=j\\ 0.5& i\≠j,\frac{a_i}{a_j}=1\\ 0.05\×\frac{a_i}{a_j}+0.45& i\≠j,1\≤\frac{a_i}{a_j}\≤9\\ 1-\[0.05\×\frac{a_j}{a_i}+0.45\]& i\≠j,\frac{1}{9}\≤\frac{a_i}{a_j}\≤1\end{array}\right.---\left(2\right)$

式中，a_i/a_j是比例标度，表示所述指标体系中的指标i与j重要性程度的比例差异；a_j/a_i表示所述指标体系中的指标j与i重要性程度的比例差异；s_ij表示得到的所述指标体系中指标i相对于所述指标体系中指标j的评分标度的值。

步骤23：通过步骤22所述评分标度，得到评分标度矩阵为：

B＝(s_ij)_m×m (3)

式中，s_ij表示得到的所述指标体系中指标i相对于所述指标体系中指标j的评分标度的值；m表示所述指标体系中指标的总个数；B表示得到的评分标度矩阵。

因B的主对角线元素均为0，所以B不是正矩阵，不存在“一致性”，更不存在“一致性检验问题”。注意到B的第i行是所述指标体系中指标i与m个指标(包含i)各进行一次比较后i的重要性得分序列，由于同一指标的重要性得分可以相加，用S(i)表示指标i的累积得分，则：

$S\left(i\right)=\underset{j=1}{\overset{m}{\Σ}}{s}_{ij}---\left(4\right)$

式中，s_ij表示得到的所述指标体系中指标i相对于所述指标体系中指标j的评分标度的值；

m个指标两两进行比较，共比较次，共得分，用r_i表示i的得分率，则：

$r_i=\frac{S\left(i\right)}{\frac{m\×(m-1)}{2}}=\frac{2\×\underset{j=1}{\overset{m}{\Σ}}{s}_{ij}}{m\×(m-1)},(i=1,2,\mathrm{...},m)---\left(5\right)$

显然r_i就是m项指标按重要度排序时所述指标体系中指标i的排序度量，俗称r_i是所述指标体系中指标i的重要性权重。

所以，m个指标的权重向量为：

r＝(r₁,r₂,…,r_m) (6)

步骤3：根据所述各指标匹配的各风险等级的临界值，将隶属度函数分为成本型隶属度函数或效益型隶属度函数；依据各风险等级隶属度函数的确定公式来确定所述成本型隶属度函数或效益型隶属度函数，并根据所述成本型隶属度函数或效益型隶属度函数得到相应的指标隶属度矩阵μ。

步骤31：依据对煤矿安全状况历史分析，科学合理的确定所述指标体系中各指标对应各风险等级的临界值，如表2所示。

表2各指标匹配的各风险等级临界值

从表2中可以看出，部分所述指标体系中的指标，如u_3-1、u_3-3、u_3-6以及u_3-7匹配的风险等级的临界值随着所述风险等级的增加而增大，符合成本型隶属度函数的特征；部分所述指标体系中指标匹配的风险等级临界值随着所述风险等级的增加而减小，符合效益型隶属度函数的特征。因此，根据所述指标体系中各指标匹配的各风险等级的临界值，将隶属度函数分为成本型隶属度函数或效益型隶属度函数。

步骤32：确定成本型隶属度函数的公式为：

$f_1\left(x\right)=\left\{\begin{array}{cc}1& x\&le;a\left(1\right)\\ \frac{a\left(2\right)-x}{a\left(2\right)-a\left(1\right)}& a\left(1\right)<x\&le;a\left(2\right)\\ 0& x>a\left(2\right)\end{array}\right.---\left(7\right)$

$f_2\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(1\right)\\ \frac{x-a\left(1\right)}{a\left(2\right)-a\left(1\right)}& a\left(1\right)<x\&le;a\left(2\right)\\ \frac{a\left(3\right)-x}{a\left(3\right)-a\left(2\right)}& a\left(2\right)<x\&le;a\left(3\right)\\ 0& x>a\left(3\right)\end{array}\right.---\left(8\right)$

$f_3\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(2\right)\\ \frac{x-a\left(2\right)}{a\left(3\right)-a\left(2\right)}& a\left(2\right)<x\&le;a\left(3\right)\\ \frac{a\left(4\right)-x}{a\left(4\right)-a\left(3\right)}& a\left(3\right)<x\&le;a\left(4\right)\\ 0& x>a\left(4\right)\end{array}\right.---\left(9\right)$

$f_4\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(3\right)\\ \frac{x-a\left(3\right)}{a\left(4\right)-a\left(3\right)}& a\left(3\right)<x\&le;a\left(4\right)\\ \frac{a\left(5\right)-x}{a\left(5\right)-a\left(4\right)}& a\left(4\right)<x\&le;a\left(5\right)\\ 0& x>a\left(5\right)\end{array}\right.---\left(10\right)$

$f_5\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(4\right)\\ \frac{x-a\left(4\right)}{a\left(5\right)-a\left(4\right)}& a\left(4\right)<x\&le;a\left(5\right)\\ 1& x>a\left(5\right)\end{array}\right.---\left(11\right)$

式(7)-(11)中，a(1),a(2),a(3),a(4),a(5)表示所述指标体系中指标匹配的各风险等级的临界值，之间的关系是：

a(1)＜a(2)＜a(3)＜a(4)＜a(5)；x表示瓦斯爆炸灾害中属于所述各指标的样本值；f₁(x),f₂(x),f₃(x),f₄(x),f₅(x)分别表示所述指标体系中各指标的样本值处于所述各风险等级的成本型隶属度函数。

效益型隶属度函数的确定公式与上述公式顺序相反。依据上述各风险等级隶属度函数的确定公式来确定所述成本型隶属度函数或效益型隶属度函数，如附图2和附图3所示。

步骤33：将所述指标体系中各指标的样本值与所述成本型隶属度函数或效益型隶属度函数相对应，并得到所述指标隶属度矩阵μ为：

式中，μ_jk表示瓦斯爆炸灾害中属于所述指标中第j个指标的样本值隶属于所述风险等级中第k个风险等级的隶属度；m表示所述指标体系中指标的总个数；p表示所述风险等级划分的总级数；μ表示瓦斯爆炸灾害中样本值的指标隶属度矩阵。

步骤4：根据非线性隶属度转换算法和所述指标体系中各指标的权重向量，将所述指标隶属度矩阵μ转换为目标隶属度向量；

设瓦斯爆炸灾害在m个指标属于所述风险等级中的第k个风险等级的目标隶属度为α_k(G)。而各种风险评估模型，几乎无例外地都用模糊综合判断的“加权平均”线性模型确定目标隶属度α_k(G)，即：

${\mathrm{\&alpha;}}_k\left(G\right)=\underset{j=1}{\overset{m}{\&Sigma;}}{r}_j\&times;{\mathrm{\&mu;}}_{jk}---\left(13\right)$

式中，r_j表示所述指标体系中指标j的权重向量；μ_jk表示瓦斯爆炸灾害中属于所述指标中第j个指标的样本值隶属于所述风险等级中第k个风险等级的隶属度。

实际上，隶属度转换是非线性转换。本发明采用非线性隶属度转换方法实现隶属度转换。由指标隶属度矩阵μ确定目标隶属度向量的非线性隶属度转换算法步骤如下：

步骤41：计算区分权w_j(G)(j＝1,2,…,m)

区分权w_j(G)的大小由所述的指标隶属度矩阵μ_jk取值集中与分散的程度决定。区分权表示的是所述指标体系中的指标j提供给所述瓦斯爆炸灾害并用隶属度向量表征的分类信息，能把所述瓦斯爆炸灾害所属风险等级“区分开”的程度。

区分权计算公式如下所示：

$H_j\left(G\right)=-\underset{k=1}{\overset{p}{\&Sigma;}}{\mathrm{\&mu;}}_{jk}\&times;{\mathrm{lg\&mu;}}_{jk}---\left(14\right)$

$V_j\left(G\right)=1-\frac{1}{lgp}\&times;H_j\left(G\right)---\left(15\right)$

$w_j\left(G\right)=\frac{V_j\left(G\right)}{\underset{j=1}{\overset{m}{\&Sigma;}}{V}_j\left(G\right)}---\left(16\right)$

式中，G表示瓦斯爆炸灾害；k表示所述风险等级中第k个风险等级；p表示所述风险等级划分的总级数；μ_jk表示瓦斯爆炸灾害中属于所述指标中第j个指标的样本值隶属于所述风险等级中第k个风险等级的隶属度；H_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的熵，刻画的是瓦斯爆炸灾害中属于所述指标体系中指标j的样本值对所述风险等级的划分所做的贡献值，熵的取值范围为[0,lgp]，熵越小，贡献越大；m表示所述指标体系中指标的总个数；V_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的峰值；w_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的区分权。

步骤42：计算可比和M_k(G)

可比和M_k(G)的物理意义是m种指标提供的分类信息使瓦斯爆炸灾害属于所述风险等级中的第k等级的“程度之和”。

$M_k\left(G\right)=\underset{j=1}{\overset{m}{\&Sigma;}}{r}_j\&times;w_j\left(G\right)\&times;{\mathrm{\&mu;}}_{jk},(k=1,2,...,p)---\left(17\right)$

式中，G表示瓦斯爆炸灾害；r_j表示所述指标体系中指标j的权重向量；w_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的区分权；w_j(G)×μ_jk表示使用区分权提取出的有效值；r_j×w_j(G)×μ_jk表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值隶属于所述风险等级中第k个风险等级的可比值；M_k(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值隶属于所述风险等级中第k个风险等级的可比和。

步骤43：计算目标隶属度α_k(G)

${\mathrm{\&alpha;}}_k\left(G\right)=\frac{M_k\left(G\right)}{\underset{j=1}{\overset{m}{\&Sigma;}}{M}_{jk}\left(G\right)}---\left(18\right)$

式中，M_k(G),M_jk(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值隶属于所述风险等级中第k个风险等级的可比和

输出的目标隶属度向量为：

α(G)＝(α₁(G),α₂(G),…,α_p(G)) (19)

此方法的特点是由状态转移矩阵实现的隶属度转换是非线性转换，并且转换过程既不增加新的先验知识，也不造成已有信息损失。

步骤5：根据置信度识别准则和所述目标隶属度向量，判断瓦斯爆炸灾害的风险等级。

置信度识别准则的定义为：

设{v₁,v₂,…,v_p}是一个有序的风险评价空间集，λ为置信度(λ＞0.5，通常取λ＝0.6或0.7)，令

$k=\min \left\{k\right|\underset{l=1}{\overset{k}{\&Sigma;}}{\mathrm{\&alpha;}}_l\left(G\right)\&GreaterEqual;\mathrm{\&lambda;},k=1,2,...,p\}$

式中，k表示所述风险等级中第k个风险等级；p表示所述风险等级划分的总级数；α_l(G)表示所述的风险等级中前k个风险等级中第l个风险等级的目标隶属度向量的值；λ表示置信度；

判断瓦斯爆炸灾害的风险等级为所述风险等级中第k个风险等级，并且有不低于λ的置信度。

Claims (5)

1.一种WMNT瓦斯爆炸灾害风险评价方法，其特征在于：包括以下步骤：

确定瓦斯爆炸灾害风险的指标体系，对瓦斯爆炸灾害风险划分风险等级；

采用改进的AHP单准则排序方法确定所述指标体系中的各指标的权重向量；

根据所述各指标匹配的各风险等级的临界值，将隶属度函数分为成本型隶属度函数或效益型隶属度函数；依据各风险等级隶属度函数的确定公式来确定所述成本型隶属度函数或效益型隶属度函数，并根据所述成本型隶属度函数或效益型隶属度函数得到相应的指标隶属度矩阵μ；

根据非线性隶属度转换算法和所述指标体系中各指标的权重向量，将所述的指标隶属度矩阵μ转换为目标隶属度向量；

根据置信度识别准则和所述目标隶属度向量，判断瓦斯爆炸灾害的风险等级。

2.根据权利要求1所述的一种WMNT瓦斯爆炸灾害风险评价方法，其特征在于：所述的采用改进的AHP单准则排序方法确定所述指标体系中的各指标的权重向量，包括：

当所述指标体系中指标i与j进行重要性比较时，用s_ij表示所述指标i相对于所述指标j的重要性得分，用s_ji表示所述指标j相对于所述指标i的重要性得分，且s_ij，s_ji满足：

$\left\{\begin{array}{c}{s}_{ij}\&GreaterEqual;0,s_{ji}\&GreaterEqual;0\\ s_{ij}+s_{ji}=\left\{\begin{array}{cc}1& i\&NotEqual;j\\ 0& i=j\end{array}\right.\end{array}\right.$

s_ij是一种标度，因为s_ij能表征所述指标体系中指标i与j重要性程度的差异，称s_ij是评分标度，同理s_ji也是评分标度。

由比例标度a_i/a_j转换成评分标度的转换公式为：

$s_{ij}=\left\{\begin{array}{cc}0& i=j\\ 0.5& i\&NotEqual;j,\frac{a_i}{a_j}=1\\ 0.05\&times;\frac{a_i}{a_j}+0.45& i\&NotEqual;j,1\&le;\frac{a_i}{a_j}\&le;9\\ 1-\&lsqb;0.05\&times;\frac{a_j}{a_i}+0.45\&rsqb;& i\&NotEqual;j,\frac{1}{9}\&le;\frac{a_i}{a_j}\&le;1\end{array}\right.$

式中，a_i/a_j是比例标度，表示所述指标体系中的指标i与j重要性程度的比例差异；a_j/a_i表示所述指标体系中的指标j与i重要性程度的比例差异；s_ij表示得到的所述指标体系中指标i相对于所述指标体系中指标j的评分标度的值；

通过所述评分标度，得到评分标度矩阵为：

B＝(si_j)_m×m

式中，m表示所述指标体系中指标的总个数；B表示得到的评分标度矩阵；

其中，B的第i行是所述指标体系中指标i与m个指标各进行一次比较后i的重要性得分序列，用S(i)表示所述指标体系中指标i的累积得分，则：

$S\left(i\right)=\underset{j=1}{\overset{m}{\&Sigma;}}{s}_{ij}$

m个指标两两进行比较，共比较次，共得分，用r_i表示i的得分率，则：

$r_i=\frac{S\left(i\right)}{\frac{m\&times;(m-1)}{2}}=\frac{2\&times;\underset{j=1}{\overset{m}{\mathrm{\&Sigma;}}}{s}_{ij}}{m\&times;(m-1)},(i=1,2,...,m)$

式中，S(i)表示所述指标体系中指标i的累积得分；表示m个指标两两进行比较后共得的分；

r_i就是m项指标按重要度排序时所述指标体系中指标i的排序度量，称r_i是所述指标体系中指标i的权重向量；

所以，m个指标的权重向量为：

r＝(r₁,r₂,…,r_m)。

3.根据权利要求1所述的一种WMNT瓦斯爆炸灾害风险评价方法，其特征在于：根据所述各指标匹配的各风险等级的临界值，将隶属度函数分为成本型隶属度函数或效益型隶属度函数；依据各风险等级隶属度函数的确定公式来确定所述成本型隶属度函数或效益型隶属度函数，并根据所述成本型隶属度函数或效益型隶属度函数得到相应的指标隶属度矩阵μ，包括：

依据对煤矿安全状况历史分析，确定所述各指标匹配的各风险等级的临界值；

部分所述指标体系中的指标匹配的风险等级的临界值随着所述风险等级的增加而增大，符合成本型隶属度函数的特征；部分所述指标体系中的指标匹配的风险等级临界值随着所述风险等级的增加而减小，符合效益型隶属度函数的特征；根据所述指标体系中各指标匹配的各风险等级的临界值，将隶属度函数分为成本型隶属度函数或效益型隶属度函数；

确定成本型隶属度函数的公式为：

$f_1\left(x\right)=\left\{\begin{array}{cc}1& x\&le;a\left(1\right)\\ \frac{a\left(2\right)-x}{a\left(2\right)-a\left(1\right)}& a\left(1\right)<x\&le;a\left(2\right)\\ 0& x>a\left(2\right)\end{array}\right.$

$f_2\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(1\right)\\ \frac{x-a\left(1\right)}{a\left(2\right)-a\left(1\right)}& a\left(1\right)<x\&le;a\left(2\right)\\ \frac{a\left(3\right)-x}{a\left(3\right)-a\left(2\right)}& a\left(2\right)<x\&le;a\left(3\right)\\ 0& x>a\left(3\right)\end{array}\right.$

$f_3\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(2\right)\\ \frac{x-a\left(2\right)}{a\left(3\right)-a\left(2\right)}& a\left(2\right)<x\&le;a\left(3\right)\\ \frac{a\left(4\right)-x}{a\left(4\right)-a\left(3\right)}& a\left(3\right)<x\&le;a\left(4\right)\\ 0& x>a\left(4\right)\end{array}\right.$

$f_4\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(3\right)\\ \frac{x-a\left(3\right)}{a\left(4\right)-a\left(3\right)}& a\left(3\right)<x\&le;a\left(4\right)\\ \frac{a\left(5\right)-x}{a\left(5\right)-a\left(4\right)}& a\left(4\right)<x\&le;a\left(5\right)\\ 0& x>a\left(5\right)\end{array}\right.$

$f_5\left(x\right)=\left\{\begin{array}{cc}0& x\&le;a\left(4\right)\\ \frac{x-a\left(4\right)}{a\left(5\right)-a\left(4\right)}& a\left(4\right)<x\&le;a\left(5\right)\\ 1& x>a\left(5\right)\end{array}\right.$

式中，a(1),a(2),a(3),a(4),a(5)表示所述指标体系中指标匹配的各风险等级的临界值，之间的关系是：a(1)＜a(2)＜a(3)＜a(4)＜a(5)；x表示瓦斯爆炸灾害中属于所述各指标的样本值；f₁(x),f₂(x),f₃(x),f₄(x),f₅(x)分别表示所述指标体系中各指标的样本值处于所述各风险等级的成本型隶属度函数；

所述效益型隶属度函数的确定公式与上述公式顺序相反；

依据上述各风险等级隶属度函数的确定公式来确定所述成本型隶属度函数或效益型隶属度函数；

将所述指标体系中各指标的样本值与所述成本型隶属度函数或效益型隶属度函数相对应，并得到所述指标隶属度矩阵μ

式中，μ_jk表示瓦斯爆炸灾害中属于所述指标中第j个指标的样本值隶属于所述风险等级中第k个风险等级的隶属度；m表示所述指标体系中指标的总个数；p表示所述风险等级划分的总级数；μ表示瓦斯爆炸灾害中样本值的指标隶属度矩阵。

4.根据权利要求3所述的一种WMNT瓦斯爆炸灾害风险评价方法，其特征在于：所述的根据非线性隶属度转换算法和所述指标体系中各指标的权重向量，将所述指标隶属度矩阵μ转换为目标隶属度向量，包括：

计算区分权w_j(G)(j＝1,2,…,m)

区分权w_j(G)的大小由所述的指标隶属度矩阵μ_jk取值集中与分散的程度决定；

区分权计算公式如下所示：

$H_j\left(G\right)=-\underset{k=1}{\overset{p}{\&Sigma;}}{\mathrm{\&mu;}}_{jk}\&times;{\mathrm{lg\&mu;}}_{jk}$

$V_j\left(G\right)=1-\frac{1}{\lg p}\&times;H_j\left(G\right)$

$w_j\left(G\right)=\frac{V_j\left(G\right)}{\underset{j=1}{\overset{m}{\&Sigma;}}{V}_j\left(G\right)}$

式中，G表示瓦斯爆炸灾害；k表示所述风险等级中第k个风险等级；p表示所述风险等级划分的总级数；μ_jk表示瓦斯爆炸灾害中属于所述指标中第j个指标的样本值隶属于所述风险等级中第k个风险等级的隶属度；H_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的熵，刻画的是瓦斯爆炸灾害中属于所述指标体系中指标j的样本值对所述风险等级的划分所做的贡献值，熵的取值范围为[0,lgp]，熵越小，贡献越大；m表示所述指标体系中指标的总个数；V_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的峰值；w_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的区分权；

计算可比和M_k(G)

$M_k\left(G\right)=\underset{j=1}{\overset{m}{\&Sigma;}}{r}_j\&times;w_j\left(G\right)\&times;{\mathrm{\&mu;}}_{jk},(k=1,2,...,p)$

式中，G表示瓦斯爆炸灾害；r_j表示所述指标体系中指标j的权重向量；w_j(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值的区分权；w_j(G)×μ_jk表示使用区分权提取出的有效值；r_j×w_j(G)×μ_jk表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值隶属于所述风险等级中第k个风险等级的可比值；M_k(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值隶属于所述风险等级中第k个风险等级的可比和；

计算目标隶属度α_k(G)

${\mathrm{\&alpha;}}_k\left(G\right)=\frac{M_k\left(G\right)}{\underset{j=1}{\overset{m}{\&Sigma;}}{M}_{jk}\left(G\right)}$

式中，M_k(G),M_jk(G)表示瓦斯爆炸灾害中属于所述指标体系中指标j的样本值隶属于所述风险等级中第k个风险等级的可比和；

输出的目标隶属度向量为：

α(G)＝(α₁(G),α₂(G),…,α_p(G))。

5.根据权利要求1所述的一种WMNT瓦斯爆炸灾害风险评价方法，其特征在于：根据置信度识别准则和所述目标隶属度向量，判断瓦斯爆炸灾害的风险等级，包括：

置信度识别准则的定义为：

设{v₁,v₂,…,v_p}是一个有序的风险评价空间集，λ为置信度(λ＞0.5，通常取λ＝0.6或0.7)，令

$k=\min \left\{k\right|\underset{l=1}{\overset{k}{\&Sigma;}}{\mathrm{\&alpha;}}_l\left(G\right)\&GreaterEqual;\mathrm{\&lambda;},k=1,2,...,p\}$

式中，k表示所述风险等级中第k个风险等级；p表示所述风险等级划分的总级数；α_l(G)表示所述的风险等级中前k个风险等级中第l个风险等级的目标隶属度向量的值；λ表示置信度；

判断瓦斯爆炸灾害的风险等级为所述风险等级中第k个风险等级，并且有不低于λ的置信度。