基干MATLAB回采巷道围岩分类可视化系统开发与`应用.docVIP

基于MATLAB回采巷道围岩分类可视化系统开发及应用 　　[摘 要]随着煤层开采深度以每年10m左右的速度递增，巷道支护问题日益突出，必须对巷道围岩稳定性进行准确分类及顶底板移近量预测。针对这个问题，依模糊聚类理论为基础，采用巷道埋深、顶板岩层综合强度、护巷煤柱宽度等7个因素为分类指标，建立回采巷道围岩稳定性分类模型，将回采巷道围岩稳定性分为五类：非常稳定，稳定，较稳定，不稳定，极不稳定；并对巷道顶底板移近量进行预测。采用matlab语言编写，实现便于工程技术人员使用的可视化系统 [关键词]支护；围岩稳定性；模糊聚类；matlab；可视化 中图分类号：TD322 文献标识码：A 文章编号：1009-914X（2016）19-0339-03 0 引言 目前，煤矿巷道支护设计很大程度上依赖于工程技术人员的工程判断力和实践经验，具有很大的盲目性，导致有些巷道支护效果不佳，部分巷道出现前掘后修、前修后坏的状况。有些巷道支护设计参数偏于保守，制约着矿井安全高效水平的提高。科学地寻找支护参数设计在安全和经济两方面的最佳结合点，是开展巷道围岩稳定性分类研究的主要目的。据统计，相当一部分巷道失稳的原因是由于基本的支护方案不合理造成的，究其原因是不准确的巷道围岩分类。因此，准确地把握巷道围岩稳定性分类是巷道支护设计的基础，开发出便于现场工程技术人员能够直接使用的可视化围岩分类系统亟待解决。本文以模糊聚类理论基础，建立回采巷道围岩稳定性分类Fuzzy模型，利用matlab语言编写模型程序，实现了回采巷道围岩稳定性分类可视化系统 1 回采巷道围岩分类Fuzzy模型 影响回采巷道围岩稳定性的因素很多，依照指标选取的重要性原则、独立性原则、明确性原则、易获取原则[1]，选取巷道埋深、巷道底板岩层强度、巷道顶板8m范围内综合强度、巷道帮部岩（煤）层强度、直接顶厚度与采高比值N、直接顶初垮步距、巷煤柱宽度7个指标，以神东布尔台矿区为工程背景，收集回采巷道样本。其中顶板岩层的综合强度，采用顶板8m范围内各个岩层的单轴抗压强度的综合平均值。大量的实践经验表明，距离巷道顶部表面的岩层越近对巷道稳定性的影响越大[2]。如图1所示，顶板岩层的强度采用如下公式（1）计算 （1） 式中： A――巷道顶板上部三角形区域面积； a――巷道宽度；Ai为第i分层区域面积； ――第i分层的单轴抗压强度 1.1 数据无量纲化及单位化处理 样本分类的7个指标的量纲不尽相同，在数据使用必须对原始数据进行无量纲化处理，按照统计学原理，本文采用公式（2）进行无纲量化处理（Z-score） （2） （3） （4） 式中： ――第i个样本的的第j个指标； ――第j指标的平均值，计算公式（3）；――第j指标的标准差，计算公式（4） 为消除各指标绝对值大小对聚类分析的影响，要对去量纲化的样本数据进行单位化[2]，如公式（4）所示，把原始数据压缩在0～1之间 （5） 式中： {xij}min――第j个指标实测中最小值； {xij}max――第j个指标实测中最大值 1.2 回采巷道分类指标加权处理 数据无量纲化及单位化处理没有改变各指标对分类结果的影响，事实上，各个分类指标对巷道围岩稳定性的影响程度是不同的，有主次之分。如果把这些影响程度不同的指标平等的对待，无疑是要影响分类结果的准确性[3]。因此，在进行模糊聚类分析时，为区分这些指标对围岩稳定性的影响程度，需要对每一个指标进行加权处理。加权的具体实施方法，就是在各指标经标准化处理后的数据上乘以相应的权值。确定权值的途径很多，这是里采用多元回归分析法确定本分类中7个指标的权值，并用层次分析法加以检验[4]，最终得到各指标分类权值如表1所示 1.3 数据标定及聚类 标定就是计算出被分类对象间相似程度的统计量（i，j=1，2，……n。n为被分类对象的个数），从而确定论域上的模糊关系矩阵。常用的方法有：欧氏距离法、数量积法、相关系数法、夹角余弦法等[5]。本文采用欧氏距离法来进行标定。如公式（6）所示 （6） 式中： ――表示第条回采样本巷道的第k个标准化处理后的指标 ； ――表示第条回采样本巷道的第k个标准化处理后的指标 聚类就是在已经建立的模糊相似矩阵的基础上[6]，以不同的阈值进行截取，从而得到不同的分类。这里采用基于模糊等价关系的聚类法。具体作法是，将加权模糊相似矩阵进行改造，使之具有传递性，转化为加权模糊等价关系矩阵，给出不同的阈值进行聚类[7]。最终把寸草塔二矿回采巷道围岩稳定性分为分为5类：非常稳定；稳定；较稳定；不稳定，极不稳定。上述对围岩稳定性的评语用拉丁字母表示为