基于Matlab采空区自燃“三带”数据处理研究.docVIP

基于Matlab采空区自燃“三带”数据处理研究 　　摘要：为了更好的划分采空区“三带”宽度，采用Matlab软件对束管检测数据进行处理，在基于氧浓度指标法前提下采用残差图剔除不可靠点后二维曲线的精确拟合计算方法以及采空区走向、倾向、氧浓度的三维等值线法直接观测“三带”范围。两种方法相互验证，最终通过应用Matlab软件的精确处理，得出采空区最大氧化带宽度为：30.4～158.95m；最小氧化带宽度为：32.17～158m。 　　关键词：采空区“三带”；氧浓度指标法；二维曲线；精确拟合；等值线 　　中图分类号：TD752.2文献标志码：A 　　[WT]文章编号：1672-1098（2012）04-0075-04 　　作者简介：张发亮（1985-），男，安徽马鞍山人，在读硕士，研究生方向：安全评价理论及技术。 　　随着煤矿开采的机械化程度的提高，综采放顶煤这种高产高效的工艺得到了广泛的应用[1]。但是，在煤矿产量提高的同时，采空区的遗煤、漏风等危险因素给工作面的安全开采带来了自燃发火的隐患。因此，对于采空区“三带”的准确预测和划分将对预防煤矿采空区的自燃发火和保障工作面的安全开采有着非常重要的意义。 　　“三带”的划分指标主要可分为3类，即以采空区内的氧气浓度、漏风风速和温度分布来划分[2-3]。目前，对于采空区“三带”的划分以氧浓度为标准的应用方式最为普遍[4]，根据煤炭科学总院重庆分院的试验研究[5-6]，不同氧化性的煤在供氧浓度大于等于5%～6%时，还可能在其着火温度前激烈氧化升温，有发火危险性，因此将氧化带划分的氧指标定为18%～6%。在此前提下，试验针对山西正株煤矿1511工作面进行研究，运用Matlab软件对束管检测的试验数据进行剔除甄别[7-8]，找出最佳的二维拟合曲线并进行三维立体成像，最终确定精确的“三带”范围。 　　1采空区自燃“三带”的现场试验 　　1.1工作面概况 　　1511工作面回采走向长度为990m，倾斜长150m；主采15煤的煤层倾角8°～10°，煤层灰分8.58%，挥发分26.05%，属于Ⅱ类易自燃煤层，煤厚最小5.12m，最大6.37m，平均5.74m。工作面采用长壁式综合机械化低位放顶煤全部垮落采煤法。工作面采高确定为2.3m，根据煤层实际厚度，考虑留设一定厚度底煤???利于拉架，确定顶煤厚度为2.9m，故平均采放比为1∶1.26。 　　1.2测点布置 　　沿工作面倾向布置5个测点，且测点由下顺槽进风巷向上顺槽回风巷依次为：1#测点，距进风巷外帮5m；2#测点，距1#测点35m；3#测点，距2#测点35m；4#测点，距3#测点35m；5#测点，距回风巷外帮5m（见表1）。在回风巷安置抽气泵，通过束管抽取气样，收集气体后将其进行色谱分析，用MATLAB处理所得的原始数据，从而研究并分析综放面采空区氧化“三带”的分布规律。 　　2基于Matlab的二维数据处理 　　2.1最小二乘法数据预处理 　　在所给出的5个检测点的数据中，随着工作面的前移，每个测点各得出19组氧浓度数据，共有95组样本点，在这些数据中可能存在有瑕疵的点，由于这些样本点会影响着实际拟合函数的结果，因此采用最小二乘法的算法方式将他们剔除[9]，实验运用多元线性回归函数regress（Y，X，alpha）和残差及其置信区间画图函数rcoplot（r，rint）进行数据处理分析，其中Y，X为所要处理的样本数据，alpha是显著性水平（取默认数值0.05），r，rint为残差及其置信区间，这2个值将由regress直接返回提供，这样，就能建立回归分析模型，主程序伪代码如下：alpha=0.05；X=[95个样本点对应的距工作面距离的数值]；Y=[95个样本点氧浓度数值]；X=[ones（size（Y）），自定义矩阵Q]；[b，bint，r，rint，stats]=regress（Y，X，0.05）；最终画出残差图（见图1）。 　　从图1中可以得出第33、56、75、81、94号这5个样本点的残差偏离了置信区间，且这5个结点对应的样本点为2#（125.8，10.2），3#（150.5，6.0），4#（150.5，5.9），5#（35.9，16.9），5#（150.5，6.1），为了更精确的分析数据，应该将这5个样本点去除以后再进行曲线回归拟合。 　　2.2拟合曲线的选取 　　在剔除偏离置信区间的样本点后，选用polytool多项式拟合工具拟合，拟合时分别选取3、4、5次多项式进行拟合，从而得到他们的拟合自由度R?分别为0.938，0.962，0.969，显然4次拟合的自由度要优于3次，由于高次多项式容易产生Runge现象[10]，亦即当次数变高时，插值多项式会逐渐偏离正常趋势线，使得拟合不精确，正常情