[发明专利]一种瓦斯抽采钻孔自燃危险预测方法

权利要求书

1.一种瓦斯抽采钻孔自燃危险预测方法，其特征在于：包括以下步骤：

步骤1：将测试点自燃危险单项指标煤碳化变质程度、煤氧化还原燃点温差、煤的吸氧能力、CO单位温升率、地质构造、煤层埋深、煤层厚度、煤层倾角、漏风强度、漏风时长、围岩温度、封孔方式、封孔长度、封孔深度、封孔材质作为评价指标，将评价指标分为主观指标和客观指标，利用改进CRITIC法对客观评价指标权重进行确定；

步骤1.1：对各评价指标包含的改进CRITIC信息量进行计算，公式如下：

式中：C_k为第k个评价指标的改进CRITIC信息量，σ_k为第k个评价指标的标准差，u_k为第k个评价指标的均值，为第k个指标和其他指标互相影响程度的量化值，t_ik为评价指标i和k的相关系数；

步骤1.2：计算各评价指标的权重w_ck，公式如下：

步骤2：利用G2赋权法对主观指标权重进行确定；

步骤2.1：根据评价目标的影响力权重，即偏好系数，对步骤1.1中包含评价指标信息量的原始指标集{u_i}中m个指标按重要性从大到小进行重新排序，排序后指标集记为{u_i1，…，u_ik，…，u_im}，其中，u_i1为最重要指标，u_im为最不重要指标，对评价指标u_ik与u_im重要程度之比a_k进行计算：

a_k＝u_ik/u_im，k＝1，2，…，m-1；

步骤2.2：当某个评价指标的信息量不足，不能计算出a_k准确值，此时赋予a_k一个取值范围，记为D_k，D_k＝[d_1k，d_2k]，通过区间长度e(D_k)以及区间中点n(D_k)的计算得出区间映射函数φε(D_k)，其中：

e(D_k)＝d_2k-d_1k

φε(D_k)＝n(D_k)+εe(D_k)

式中：ε为风险态度因子(|ε|≤0.5)，当-0.5≤ε≤0时，为保守型；当ε＝0时，为中立型；当0≤ε≤0.5时，为风险型；

若{D_k}赋值准确，则计算第k个指标的G2法权重w_Gk，如下：

步骤3：利用改进CRITIC法对G2法进行修正，确定评价指标的综合权重；

将各评价指标改进CRITIC信息量C_k的置信度区间记为[C_1k，C_2k]，其中k＝1，2，…，m，利用改进CRITIC信息量的置信度区间上下限之比计算两指标重要性程度之比的取值区间，代替G2法中的理性赋值区间，即：

式中：R_km为置信度区间；

计算第k个评价指标改进CRITIC修正G2区间判指标的赋值权重w为：

步骤4：构建初始评判矩阵；

步骤5：确定正负理想解计算相似度；

步骤6：构建G2-TOPSIS综合评判模型，根据组合赋权原理，确定准则层指标的加权综合判据矩阵Z_ij与基于改进CRITIC改进G2法重新排序后对准则层指标打分所构成的评判矩阵V结合得出评判对象综合评判结果L；

L＝Z_ijV

步骤7：对抽采钻孔自燃危险等级进行建立，将抽采钻孔自燃危险分为5个级别，I级-基本不可能发生、II级-较不可能发生、III级-可能发生、IV级-较可能发生、V级-极可能发生，各等级赋值均在[0,100]，各指标赋值越高说明该指标引起抽采钻孔自燃危险性越小，抽采钻孔发生自燃的概率越小。