基于灰色系统矿井涌水量预测分析.docVIP

基于灰色系统矿井涌水量预测分析 　　[摘要]灰色预测模型已被用于矿井涌水量的预测，但灰色数据列中X（0）（1）作为一个灰数，带有一定程度的随机性，且X（1）（1）也没有被完全白化，故模型的精度对X（0）（1）的可靠度具有很大依赖性。本文尝试采用重复建模求英平均值的方法，提高X（0）（1）白化程度，以削弱其随机性提高GM模型的精度。并以济宁三号井为例进行涌水量预测。 　　[关键字]涌水量 灰色预测 模型精度 　　[中图分类号] TD742 [文献码] B [文章编号] 1000-405X（2013）-5-137-2 　　在采矿生产中，矿井涌水是由许多因素共同决定的，如水文地质条件、地质构造、降水量及开采阶段等。人们能在一定程度上对矿井涌水量的影响因素（矿区地质、水文地质条件等）有所了解，但不可能了解的很透彻。故在矿井涌水系统中，既有已知信息又有未知信息，实际上是一个灰色系统。 　　前人已应用灰色系统对矿井涌水量进行了预测.利用离散数据建立微分方程的动态模型（Greymodel），一般采用一阶单序列GM（1，1）。 　　X（1）（k+1）=（X（0）（1）-u/a）e-ak+u/a其中X（0）称为X（1）的初始值，模型曲线必过该定点，但是由灰色理论可知，在灰色数列X（0）中，X（0）（1）为一灰数，带有很大的随机性，即使是X（1）也是一阶弱随机序列，其中的X（1）（1）也没有被完全白化。因此，X（1）（1）的取值有可能与X（1）所蕴含的规律相悖，这时模型的精度就必然受到影响。且经过对一些文献的分析，发现初值法的确存在不稳定性。为了提高模型精度，本文尝试采用依次去掉老数据分别重复建模的方法，得到一预测区间，并对预测区间求其平均值。其优点在于：一是削弱了X（0）（1）的随机性。二是随时间的变化，涌水量变化趋势也在改变，维数的减少从一定程度上更能反映涌水量新的变化趋势。 　　l基本原理 　　1.1GM（1，1）建模 　　首先对原始数据序列x（0）={x（0）（i），i=1，2，…，N}据式 求出相应的一阶累加序列X（1）={X（i）（i），i=1，2，....，N} 　　求微分方程系数向量： 　　得微分方程动态模型及时间响应函数： 　　据上式求生成数的回代计算值 　　求原始数据还原值： ???　1.2精度检验 　　采用后验检验法，将精度分为四级（表1）。 　　检验精度是否达到要求，如未达到，可进行残差修正。 　　1.3残差修正 　　求原始序列和模型计算值之间的残差： 　　若残差序列中的最小值ε2min小于0，则将横坐标下移2×|ε（0）min|，得非负残差序列ε（0），即 。 　　将ε（0）（i）作一阶累加生成，建立GM（1，1）模型求出回代计算值ε（1）（i）.并还原得ε（0）（i）。 　　将横坐标上移2×|ε（0）min|得ε（0）（i），即 。 　　用ε（0）（i）修正同时刻的 。即： 。 　　1.4重复建模 　　将原始数列中的x（0）（1）、x（0）（2）、x（0）（3）…依次作为x（0）（1）分别建模，得到多个时间响直函数，同时进行多次预测，输出一个预测区间，并对区间取其平均值.该值即为预测结果。 　　2以济宁三号煤矿涌水量预测为例 　　济宁三号煤矿目前主要开采3煤层，直接充水含水层为山西组砂岩和太原组第三层石灰岩，间接充水含水层为上侏罗统砂岩含水层。山西组砂岩漏水孔率为23.8%，单位涌水量0.1301～0.4755L（s·m），富水性中等；太原组第三层石灰岩漏水孔率为1.2%—14.5%，单位涌水量0.00306～0.00701L/（s·m），富水性较弱。至2012年底，矿井年平均涌水量最大为351.6m3/h，最小为48.9m3/h，平均为280.0m3/h，单月最大涌水量为552.5m3/h。 　　近年来，矿井涌水量平均值和最大值均呈现出上升的趋势；2001年以后矿井涌水量平均值则处于相对稳定状态，而涌水量最大值则有逐年下降的趋势。主要原因可能在于；投产初期，随着开采面积的逐步增加以及开采上侏罗统砂砾岩水下泄，矿井涌水量处于不稳定状态，呈现出上升的趋势；经过几年连续较大强度的排水，上侏罗统水位大幅度下降，上侏罗统对其他采区基本上没有直接的影响，因此，其对开采产生的影响逐渐减弱；3煤顶板砂岩的富水性和补给条件一般，初次揭露时有较大的出水，但随着开采，水位下降并逐渐被疏放，最大涌水量逐渐减小。本文采用2006年以后数据进行建模（表如下）。 　　2.1年最大涌水量预测 　　（1）2006-2012年数据建模 　　取2006-2012年数据建立GM（1，1）模型，得到a=0.0901，u=612.2217，经过后验差检验，c0=0.4569，p=0.8571.符合预