一种用于深井底板突水风险程度的判别方法
【专利摘要】本发明公开了一种深井底板突水风险程度的判别方法，属于矿山安全分析领域。其首先采集影响底板突水的多个主控因素，然后根据主控因素权重选择评价指标，确定评价指标的隶属度函数，再结合实际矿井突水情况，选择底板风险程度明显的案例，并对其样本数据做归一化处理，运用MATLAB程序对案例进行分类，构建底板突水危险程度案列数据库；最后输入待测底板危险程度的评价指标数据来确定待测底板的突水危险等级。本发明方法可以对多个主控因素进行评价，为提高深井开采的安全性提供了技术支持。
【专利说明】
一种用于深井底板突水风险程度的判别方法
技术领域
[0001] 本发明属于矿山安全分析领域，具体涉及一种用于深井底板突水风险程度的判别 方法。
【背景技术】
[0002] 煤层底板突水是一个比较复杂的过程，是人为采掘工程活动导致煤层围岩体应力 场能量释放、煤层底板隔水岩体破坏、矿井局部充水水文地质条件突变的一种岩体失稳现 象。目前，深井底板突水风险程度的安全评价方法一般从数理统计、水力学、模糊数学、灰色 理论、神经网络等方面对底板突水危险程度进行评价，但这些方法还不成熟，且多以单一运 算为主，不能形成深井底板突水风险案例数据库。
[0003] 目前现有技术中有关矿井底板突水风险程度判别方法主要有"突水系数法"、"下 三带"理论和"下四带理论"，然而，这些底板突水评价还存在以下主要问题：
[0004] (1) "下三带"理论和"下四带"理论在采动条件下的矿井底板隔水层破坏情况和实 际情况比较吻合，但是由于现场实测费用较高、操作程序较复杂、数据获取较困难，因此没 有被现场广泛应用。
[0005] (2) "突水系数法"因其物理概念简单，计算方便，因而一直是我国煤层底板突水评 价的方法，但此方法仅考虑含水层水压和底板隔水层总厚度两个影响因素，未能体现煤层 底板突水受控于多因素影响且具有非常复杂机理的非线性动力现象。
[0006] 以模糊数学、灰色关联度、神经网络等理论或交叉理论为基础的煤层底板突水评 价方法，多以单一运算为主，没有形成底板突水风险案例库，不能将现有的突水案例加以利 用，造成资源的浪费。
【发明内容】
[0007] 针对上述现有技术中存在的技术问题，本发明提供了一种深井底板突水风险程度 的判别方法，该方法可以对多个主控因素进行评价，为提高深井开采的安全性提供了技术 支持。
[0008] 其技术解决方案包括：
[0009] -种深井底板突水风险程度的判别方法，依次包括以下步骤：
[0010] a采集影响底板突水的多个主控因素数据，所述主控因素包括含水层水压、断层导 水性、隔水层厚度、工作面斜长、强水源补给度、开采深度和煤层倾角；
[0011] b根据主控因素权重选择几个评价指标，并确定每个评价指标的隶属度函数；
[0012] c根据实际矿井突水情况，选取多个底板突水危险程度明显的样本数据，运用 MATLAB模糊C-均值综合分析法对其进行分类，将底板突水危险程度分为不突水、一般突水 和严重突水三类，其中，MATLAB模糊C-均值综合分析法计算方法如下：
[0013] (^初始化:取模糊加权指数m=2,聚类的类别数C(2<(Xn)，n为数据样本点的个 数，迭代停止阈值ε，初始的聚类中心值P (())，以及迭代次数1 = 0;
[0014] (32计算由隶属度的值所组成的划分矩阵U(1);
[0015] 对于任意的i，k，如果dik(1)>0，则
[0017] 对于任意的i，r，如果dik(1) = 0,贝丨J
[0018] jiir⑴=1，且 j辛r时，yij(1) = 0;
[0019] C3更新聚类中心值：
[0021] C4若| | PUP11 | <ε，则算法停止，否则转到步骤c2;
[0022] d根据步骤c的分类情况，将得到的数据存入底板突水危险程度案例数据库；
[0023] e将待测底板危险状况的评价指标数据输入到数据库中，运用步骤c所述的MATLAB 模糊C-均值综合分析法计算，根据所述底板突水危险程度案例数据即可确定底板突水危险 情况。
[0024] 作为本发明的一个优选方案，步骤b中，根据主控因素权重选择4个评价指标。
[0025] 作为本发明的另一个优选方案，步骤c中，选取6个底板突水危险程度明显的样本 数据。
[0026]本发明所带来的有益技术效果：
[0027] 本发明提供的一种用于深井底板突水风险程度的判别方法，对多个主控因素确 定，根据主控因素的重要度选取评价指标，对实际突水风险情况明显的矿井进行分类，根据 分类情况构建底板突水风险程度数据库，将待评价矿井底板突水评价指标输入数据库，确 定矿井底板突水风险级别，进而为提尚深井开米的安全性提供技术支持。
【附图说明】
[0028] 下面结合附图对本发明做进一步说明：
[0029] 图1为本发明判别方法的流程图。
【具体实施方式】
[0030] 本发明提出了一种用于深井底板突水风险程度的判别方法，为了使本发明的优 点、技术方案更加清楚、明确，下面结合具体实施例对本发明做详细说明。
[0031] 实施例1:
[0032] 用于深井底板突水风险程度的判别方法，如图1所示的，包括以下步骤：
[0033] 步骤一、采集影响底板突水的多个主控因素数据，分别包括:水文地质条件、导水 构造状态、底板隔水层状态、开采状态和深井状态；
[0034] 步骤二、根据专家打分方式，计算分析各因素的权重大小，选择等4个因素作为此 次底板突水危险程度的评价指标；
[0035] 步骤三、根据专家打分计算各因素的权重大小，含水层水压0.24、隔水层厚度 0.20、断层导水性0.16、工作面斜长0.14、强水源补给度0.10、开采深度0.09、煤层倾角 〇. 07，重要度程度顺序如下:含水层水压〉隔水层厚度〉断层导水性〉工作面斜长〉强水源补 给度〉开采深度〉煤层倾角，选择含水层水压、隔水层厚度、断层导水性、工作面斜长等4个主 控因素作为此次底板突水危险程度的评价指标；
[0036] 步骤四、上述含水层水压、隔水层厚度、断层导水性、工作面斜长这4个评价指标的 隶属度函数分别为：
[0037]隔水层厚度的隶属度函数为：
[0039]上式中，h-底板隔水层厚度；
[0040]含水层水压的隶属度函数为：
[0042]上式中，j 一底板含水层水压；
[0043]工作面斜长的隶属度函数为：
[0045] 上式中，k 一工作面斜长；
[0046]步骤五、选取现实中6个底板突水危险程度明显的样本数据，运用MATLAB模糊C-均 值综合分析法对其进行分类，将底板突水危险程度分为不突水、一般突水和严重突水三类； [0047]上述6个底板突水危险程度明显的样本数据如表1所示：
[0048] 表1底板突水样本数据
[0049]
[0050] 6个底板突水危险程度明显的样本数据归一化处理后的数据如表2所示：
[00511表2底板突水归一化处理数据 [0052]
[0053]上述MATLAB模糊C-均值综合分析法的程序语言为：
[0054] [center,U,obj_fcn] = fcm(x,3)；
[0055] maxU=max(U)；
[0056] indexl =f ind(U(l ,:) = =maxU)
[0057] index2 = find(U(2, :) = =maxU)
[0058] index3 = find(U(3, :) = =maxU) 〇
[0059] 计算方法如下：
[0060] 第一步、初始化:取模糊加权指数m = 2,聚类的类别数C(2<C<n)，n为数据样本点 的个数，迭代停止阈值ε，初始的聚类中心值P(())，以及迭代次数1 = 0;
[0061] 第二步、计算由隶属度的值所组成的划分矩阵U(1);
[0062] 对于任意的i，k，如果dik(1)>0,则
[0064] 对于任意的i ,r，如果dik(1) = 0,则
[0065] yir⑴=1，且 j辛r时，μυ(1) = 0;
[0066] 第三步、更新聚类中心值：
[0068]第四步、若| ΙΡ^Ρ" |〈ε，则算法停止，否则转到第二步；
[0069] 根据上述MATLAB运行结果并结合矿井实际突水情况得出，样本5和样本6为不突 水，样本3和样本4为一般突水，样本1和样本2为严重突水。
[0070] 步骤六、根据分类情况，将6个底板突水危险案例构建成底板突水危险程度数据 库。
[0071] 步骤七、将待测底板危险状况的评价指标数据输入到数据库中，运用MATLAB模糊 C-均值综合分析法，确定底板突水危险情况。
[0072] 上述步骤七中MATLAB模糊C-均值综合分析法的程序语言为：
[0073] e = ones(3, l)*xl;
[0074] f=(center-e)'
[0075] ff = sum(f. '2)；
[0076] [mini, index] =min(ff)；
[0077] disp(['此底板危险的类别为'，num2str(index), '类'])
[0078]计算方法与步骤五相同。
[0079]检测结果:待测样本的数据是含水层水压7MPa，断层导水性0.5，隔水层厚度30m， 工作面斜长ll〇m，将待测样本底板数据输入数据库中，得出此样本底板突水情况为一般突 水。
[0080] 实施例2:
[0081] 与实施例1不同之处在于:步骤五中选取现实中9个底板突水危险程度明显的样本 数据，9个底板突水危险程度明显的样本数据如表3所示，后续的归一化处理方法等与实施 例1类似。
[0082]表3底板突水样本数据
[0083]
[0084] 需要说明的是，本发明所例举的多个主控因素不仅仅局限于含水层水压、断层导 水性、隔水层厚度、工作面斜长、强水源补给度、开采深度和煤层倾角，在本说明书的教导下 本领域技术人员所做出的任何等同方式，或明显变型方式均应在本发明的保护范围内。
【主权项】
1. 一种深井底板突水风险程度的判别方法，其特征在于，依次包括W下步骤： a采集影响底板突水的多个主控因素数据，所述主控因素包括含水层水压、断层导水 性、隔水层厚度、工作面斜长、强水源补给度、开采深度和煤层倾角； 时良据主控因素权重选择几个评价指标，并确定每个评价指标的隶属度函数； C根据实际矿井突水情况，选取多个底板突水危险程度明显的样本数据，运用MATLAB模 糊C-均值综合分析法对其进行分类，将底板突水危险程度分为不突水、一般突水和严重突 水Ξ类，其中，MATLAB模糊C-均值综合分析法计算方法如下： C1初始化:取模糊加权指数m = 2,聚类的类别数0(2《0《11)，11为数据样本点的个数，迭 代停止阔值ε，初始的聚类中屯、值PW，W及迭代次数1 = 0; C2计算由隶属度的值所组成的划分矩阵UW ; 对于任意的i，k，如果dikW〉0,则对于任意的i，r，如果dikW = 0,则 yir(" = l，且 j声r时，μυ山=0; C3更新聚类中屯、值：C4若II pi+1-pi II <ε，则算法停止，否则转到步骤C2; d根据步骤C的分类情况，将得到的数据存入底板突水危险程度案例数据库； e将待测底板危险状况的评价指标数据输入到数据库中，运用步骤C所述的MATLAB模糊 C-均值综合分析法计算，根据所述底板突水危险程度案例数据即可确定底板突水危险情 况。2. 根据权利要求1所述的深井底板突水风险程度的判别方法，其特征在于:步骤b中，根 据主控因素权重选择4个评价指标。3. 根据权利要求1所述的深井底板突水风险程度的判别方法，其特征在于:步骤C中，选 取6个底板突水危险程度明显的样本数据。
【文档编号】G06Q10/06GK106095990SQ201610453481
【公开日】2016年11月9日
【申请日】2016年6月21日
【发明人】刘伟韬, 孟祥喜, 申建军, 赵吉园, 宋文成, 李光, 李强, 王建宁, 王东辉, 张勤, 范龙, 张茂鹏, 穆殿瑞, 谢祥祥, 董文程, 宰慧, 刘士亮, 廖尚辉
【申请人】山东科技大学