矿井瓦斯涌出量预测

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,矿井瓦斯涌出量预测,1,基本概念,矿井瓦斯涌出量,在矿井建设和生产过程中涌出到采掘空间的瓦斯数量。,绝对瓦斯涌出量,单位时间内涌出的瓦斯量，,m,3,/min, m,3,/d,。,相对瓦斯涌出量,矿井在正常生产条件下，平均日产一吨煤所涌出的,瓦斯量，,m,3,/t,。,式中,q ,相对瓦斯涌出量，,m,3,/t,；,Q ,绝对瓦斯涌出量，,m,3,/d,；,A ,平均日产煤量，,t/d,；,2,矿井瓦斯等级,根据矿井相对瓦斯涌出量、矿井绝对瓦斯涌出量和瓦斯涌出形式划分为三类：,(1) 低瓦斯矿井：矿井相对瓦斯涌出量小于或等于10,m,3,/t,且矿井绝对瓦斯涌出量小于或等于40,m,3,/min。,(2),高瓦斯矿井：矿井相对瓦斯涌出量大于10,m,3,/t,或矿井绝对瓦斯涌出量大于40,m,3,/min。,(3),煤（岩）与瓦斯（二氧化碳）突出矿井。,3,矿井瓦斯涌出量预测方法研究现状,从国内外研究现状来看，矿井瓦斯涌出量预测方法可分为两类：一类是建立在数理统计基础上的,矿山,统计法，这种方法依据矿井瓦斯涌出量随开采深度变化的统计规律，外推到预测的新区；另一类是以煤层瓦斯含量为基本预测参数的瓦斯含量法。,瓦斯地质数学模型法是近年来提出的一种新的预测方法，这种方法通过建立包括开采深度在内的多变量数学模型，对矿井未采区域的瓦斯涌出量进行预测。,4,1,矿山统计法,1.1,瓦斯梯度法,所谓瓦斯梯度，是指甲烷带内相对瓦斯涌出量平均每增加,1m,3,/t,时深度的增加量 。,式中,H,1,、H,2,瓦斯风化带以下两次测定涌出量的深度，,m，H,2,H,1,；,q,1,、q,2,对应于,H,1,、H,2,的相对瓦斯涌出量，,m,3,/t。,5,q,0,q,1,q,2,q,H,H,2,H,1,H,0,6,利用求得的瓦斯梯度，可对深部的瓦斯涌出量进行预测：,式中,q,待求深度的相对瓦斯涌出量，,m,3,/t；,H,对应于,q,的深度，,m。,7,1.2,一元回归法,如果已采区域内有比较多的瓦斯涌出量实测数据，采用回归分析方法建立瓦斯涌出量对深度的一元线性回归方程，可得到更高的预测精度。,瓦斯涌出量实测数据,单元,编号,开采深度,H,（,m,）,相对瓦斯涌出量,q,(m,3,/t),1,180,5.8,2,320,8.6,n,420,11.9,8,H,q,b,0,q=b,0,+b,1,H,9,统计检验,相关系数：,q,H,10,曲线回归,令,:,则,:,q,H,11,1.3,采用矿山统计法预测时需注意的问题,如果某一井田范围内存在较明显的构造分区，而且受地质构造等因素的影响，不同构造分区的瓦斯涌出量随开采深度的变化梯度有明显差异，井田内采用一个瓦斯涌出量梯度或回归方程进行预测，会造成预测误差增大，达不到生产要求。在这种情况下，,应在各个构造分区分别计算瓦斯涌出量梯度或回归方程，分区进行瓦斯涌出量预测。,12,2,瓦斯含量法,这种方法以煤层瓦斯含量为基本的预测参数，通过计算井下各涌出源的瓦斯涌出量 ，得到矿井或某一预测范围的瓦斯涌出量预测值。,2.1,矿井瓦斯涌出来源,矿井瓦斯的涌出来源如图所示，包括7个基本涌出源：,开采煤层瓦斯涌出,矿井瓦斯涌出,生产采区瓦斯涌出,已采采区采空区瓦斯涌出,回采工作面瓦斯涌出,生产采区采空区瓦斯涌出,掘进巷道瓦斯涌出,邻近煤层瓦斯涌出,围岩瓦斯涌出,巷道煤壁瓦斯涌出,掘进落煤瓦斯涌出,矿井瓦斯涌出来源,13,（,1,）,开采煤层瓦斯涌出,开采煤层瓦斯涌出量的大小主要受以下一些因素的影响：,煤层瓦斯含量,煤层瓦斯含量越大，涌出的瓦斯量越多。,煤的残存瓦斯含量,煤中的瓦斯全部解吸需要较长的时间。煤炭运至地表后仍然含有部分瓦斯。运至地表后单位重量的煤中所含的瓦斯量称为煤的残存瓦斯含量。残存瓦斯含量的大小主要与煤的变质程度有关。通常，煤的变质程度越高，其残存瓦斯含量越大，这是由于不同变质程度的煤其孔隙结构不同。,煤层厚度与采高的比值,相对瓦斯涌出量是根据绝对瓦斯涌出量和采出的煤量计算的。如果煤层厚度大，不能一次采全高，没有采出的煤也会向采掘空间涌出瓦斯。所以，煤层厚度与采高的比值越大，相对瓦斯涌出量也越大。,14,掘进巷道排放瓦斯的影响,工作面回采之前，已掘进的回风巷和运输巷会使工作面内沿巷道方向形成一定宽度的瓦斯排放带。瓦斯排放带越宽，回采过程中涌出的瓦斯量将越小。,工作面丢煤的影响,回采工作面的煤并不能全部开采出来，都用一定的回采率。丢失在工作面内的煤也会涌出瓦斯。而相对瓦斯涌出量是根据采出的煤量计算的，所以工作面丢煤越多，相对瓦斯涌出量越大。,15,（,2,）邻近煤层瓦斯涌出,煤层开采过程中，围岩原有的应力状态遭到破坏，煤层顶板会依次出现冒落带、裂隙带、整体移动带，煤层底板会依次出现裂隙带、整体膨胀变形带。位于顶底板冒落带和裂隙带内的邻近煤层中的瓦斯，在瓦斯压力梯度作用下，会不断涌入开采煤层的采掘空间。,邻近煤层的瓦斯涌出量主要与下面一些因素有关：,邻近层与开采层的层间距,层间距越小，邻近层越容易向开采煤层涌出瓦斯。,邻近层的瓦斯含量和残存瓦斯含量,邻近层的瓦斯含量越大、残存瓦斯含量越小，涌入开采煤层采掘空间的瓦斯量就会越大。,邻近层数目,要计算邻近煤层的瓦斯涌出量，需要将若干临近煤层的瓦斯涌出量累加起来,16,邻近层厚度与采高的比值,邻近层厚度相对于开采层采高的比值越大，以开采煤量计算的相对瓦斯涌出量应该越大。,（,3,）,围岩的瓦斯涌出,煤层的围岩中也含有一定数量的瓦斯，它们以游离状态存在于围岩的孔隙和裂隙之中。煤层开采过程中，围岩中的瓦斯会很快从岩层中涌出。但由于围岩的瓦斯含量目前还难以准确测定，因而煤层围岩的瓦斯涌出量通常采用比例系数来估算，一般按开采层涌出量的5%-25%估算煤层围岩的瓦斯涌出量。,（,4,）掘进巷道煤壁瓦斯涌出,讨论巷道煤壁的瓦斯涌出量，涉及到巷道煤壁瓦斯涌出强度（单位时间、单位面积煤壁的瓦斯涌出量）。,17,根据前苏联的研究，煤壁瓦斯涌出强度符合下列经验公式：,式中,Q,t,t,时刻煤壁瓦斯涌出强度，,m,3,/(m,2,min,)；,Q,0,煤壁瓦斯涌出初始强度，,m,3,/(m,2,min,) ；,t,煤壁暴露时间，,min,。,初始强度,Q,0,可以在巷道中实测。从公式可以看出，煤壁瓦斯涌出强度,Q,t,与暴露时间,t,呈幂函数关系。,设,L,巷道长度，,m;,n,煤壁暴露面数；,m,0,巷道中的煤层厚度，,m；,18,V,掘进速度，,m/min,；,Q,1,巷道煤壁绝对瓦斯涌出量，,m,3,/min,。,则一分钟的掘进面积为,nm,0,V,。,于是，整条巷道的绝对瓦斯涌出量（,m,3,/min,）,相当于一分钟的掘进面积连续涌出瓦斯,L/V,分钟。因此，巷道煤壁绝对瓦斯涌出量可表示为：,根据,Q,0,的物理意义，当,t=0,时，,Q,1,= nm,0,VQ,0,，,故,C=-nm,0,VQ,0,，,于是,19,（,5,）掘进落煤瓦斯涌出,掘进落煤绝对瓦斯涌出量,Q,2,（m,3,/min）,可表示为：,Q,2,=,每分钟落煤量(煤层瓦斯含量煤的残存瓦斯含量),（,6,）采空区瓦斯涌出,采空区可分为生产采区采空区和已采采区采空区。生产采区采空区包括生产采区内的已采工作面、残留煤柱、已掘巷道、报废巷道、采区内各种峒室等，已采采区采空区包括老采区、开拓巷道、井底车场内各种峒室等。,从目前的研究现状来看，影响采空区瓦斯涌出的因素很多，采空区瓦斯涌出量很难准确计算。各国都是采用采空区瓦斯涌出系数来估算采空区的瓦斯涌出量。采空区瓦斯涌出系数可分为生产采区采空区瓦斯涌出系数和已采采区采空区瓦斯涌出系数。,20,2.2,前苏联提出的预测公式,（,1,）,开采煤层（包括围岩）相对瓦斯涌出量,式中,q,k,开采煤层（包括围岩）相对瓦斯涌出量，,m,3,/t；,k,1,围岩瓦斯涌出系数。全部陷落法管理顶板，,k,1,=1.25,；,局部充填法，,k,1,=1.20,；,全部充填法，,k,1,=1.10,；,m,0,煤层厚度，,m,；,m,煤层采高，,m,；,k,2,掘进巷道瓦斯排放系数，,k,2,=(L-2h)/L。,前进式开采时，如上,部相邻工作面已采，,k,2,1；,如上部相邻工作面未采，,k,2,=(L+2h+2b)/(L+2b),。,L,采煤工作面长度，,m,；,h,瓦斯排放带宽度，无,烟煤及贫煤,h=10m,，,瘦煤及焦煤,h=14m,，,其它煤种,h=18m,；,21,k,3,煤柱瓦斯涌出系数，,k,3,=,l/L,；,l,煤柱沿倾斜方向的宽度，,m,；,x,0,煤层瓦斯含量，,m,3,/t,；,x,1,煤的残存瓦斯含量，,m,3,/t,。,不同变质程度煤的残存瓦斯含量见下表。,运到地表的煤中残存瓦斯含量,煤的挥发分,V,r,(%),残存瓦斯含量,x,1,(m,3,/t,可燃质,),28,128,812,87,1218,76,1826,65,2635,54,3542,43,4250,32,22,表中的残存瓦斯含量,x,1,，,应用时应按下式换算为煤层原始瓦斯含量,x,1,：,式中,W,f,、,A,f,煤的水分和灰分，,%,。,（,2,）邻近煤层相对瓦斯涌出量,式中,q,l,邻近煤层相对瓦斯涌出量，,m,3,/t；,n,邻近层数目；,m,i,第,i,邻近层厚度，,m；,m,开采层采高，,m；,23,x,0i,第,i,邻近层瓦斯含量，,m,3,/t；,x,1i,第,i,邻近层残存瓦斯含量，,m,3,/t；,h,i,第,i,邻近层与开采层的层间距，,m；,h,p,采动后煤层顶底板岩石受到影响的范围，,m。,顶板的影响范围按下式计算：,h,p,=Z,k,m(1.2+cos,),式中,煤层倾角，度；,m,开采层采高，,m；,Z,k,与顶板管理方法有关的系数。采高2.5,m，,全部陷落法管理,顶板时，,Z,k,=60；,局部充填法，,Z,k,=45；,全部充填法:,Z,k,=25,。,底板的影响范围：,h,p,=35m；,急倾斜煤层，,h,p,=Z,k,m(1.2-cos,)。,当某一邻近煤层已先期开采时，开采煤层已作为邻近层排放了瓦斯。此时，开采煤层的瓦斯含量,x,0,已降低为,x,0,：,24,（,3,）掘进巷道煤壁绝对瓦斯涌出量,式中,Q,1,掘进巷道煤壁绝对瓦斯涌出量，,m,3,/min,；,Q,0,煤壁瓦斯涌出初始强度，,m,3,/(m,2,min,)，可用经验公式推算：,（,4,）掘进落煤绝对瓦斯涌出量,25,式中,Q,2,掘进落煤绝对瓦斯涌出量，,m,3,/min,；,m,0,巷道中的煤层厚度，,m；,b,巷道宽度，,m；,V,掘进速度，,m/min；,d,煤的容重，,t/m,3,。,（,5,）矿井相对瓦斯涌出量,式中,q,矿,矿井相对瓦斯涌出量，,m,3,/t,；,K,采空区瓦斯涌出系数，1.151.25；,n,矿井内采煤区个数；,q,ki,第,i,采煤区开采层相对瓦斯涌出量，,m,3,/t；,q,li,第,i,采煤区邻近层相对瓦斯涌出量，,m,3,/t,；,26,A,i,第,i,采煤区平均日产煤量，,t；,m,矿井内掘进巷道条数；,Q,1j,第,j,掘进巷道煤壁绝对瓦斯涌出量，,m,3,/min；,Q,2j,第,j,掘进巷道落煤绝对瓦斯涌出量，,m,3,/min；,A,矿井平均日产煤量，,t。,2.3,分源预测法抚顺分院的改进公式,（,1,）开采煤层（包括围岩）相对瓦斯涌出量,薄及中厚煤层不分层开采时按下式计算：,27,式中,q,k,开采煤层（包括围岩）相对瓦斯涌出量，,m,3,/t；,k,1,围岩瓦斯涌出系数。其值取决于回采工作面顶板管理方法：,全部陷落法管理顶板时，,k,1,1.20,；,局部充填法管理顶板时，,k,1,=,1.15,；,全部充填法管理顶板时，,k,1,1.10；,k,2,工作面丢煤瓦斯涌出系数,，其值为工作面回采率的倒数；,k,3,准备巷道排瓦斯对工作面煤体瓦斯涌出影响系数；计算公式,同前苏联公式，,瓦斯排放带宽度,h,见后表。,m,0,煤层厚度（,夹矸层按层厚1/2计算,），,m；,m,煤层开采厚度，,m；,x,0,煤层原始瓦斯含量，,m,3,/t；,x,1,煤的残存瓦斯含量，,m,3,/t，,按下表选取。,28,运至地表时煤中残存瓦斯含量,煤的挥发份含量,V,daf,(%),6,8,8,12,12,8,18,26,26,35,35,42,42,50,纯煤残存瓦斯含量,X,1,( m,3,/t),9,6,6,4,4,3,3,2,2,2,2,表中残存瓦斯含量的单位为每一吨纯煤(即无灰干燥煤)的瓦斯体积，在应用开采煤层相对瓦斯涌出量计算公式时，应按下式换算为原煤残存瓦斯含量：,式中,x,1,表中查出的纯煤残存瓦斯含量，,m,3,/t,；,A,ad,原煤中灰份含量，,%,；,M,ad,原煤中水份含量，%。,29,巷道预排瓦斯等值宽度,h,巷道煤壁,暴露时间,(,d),不同煤种巷道预排瓦斯等值宽度（,m,）,无,烟,煤,瘦煤,焦煤,肥煤,气煤,长,焰,煤,25,6.5,9.0,9.0,11.5,11.5,11.5,50,7.4,10.5,10.5,13.0,13.0,13.0,100,9.0,12.4,12.4,16.0,16.0,16.0,160,10.5,14.2,14.2,18.0,18.0,18.0,200,11.0,15.4,15.4,19.7,19.7,19.7,250,12.0,16.9,16.9,21.5,21.5,21.5,300,13.0,18.0,18.0,23.0,23.0,23.0,30,厚煤层分层开采时按下式计算：,式中,k,fi,取决于煤层分层数量和顺序的分层开采瓦斯涌出系数，,k,fi,可按下表选取。,厚煤层分层开采瓦斯涌出系数,k,f,两,分,层,开,采,三,分,层,开,采,k,f1,k,f2,k,f1,k,f2,k,f3,1.504,0.496,1.820,0.692,0.488,31,（,2,）邻近煤层相对瓦斯涌出量,邻近煤层相对瓦斯涌出量,q,l,的计算公式与前苏联公式相同。,（,3,）掘进巷道煤壁绝对瓦斯涌出量,掘进巷道煤壁绝对瓦斯涌出量,Q,1,的计算公式与前苏联公式相同。,（,4,）掘进落煤绝对瓦斯涌出量,掘进落煤绝对瓦斯涌出量,Q,2,的计算公式与前苏联公式相同。,（,5,）生产采区相对瓦斯涌出量,32,式中,q,k,生产采区相对瓦斯涌出量，,m,3,/t,；,k,生产采区内采空区瓦斯涌出系数，,1.15,1.25,；,n,采区内回采工作面个数；,q,ki,第,i,回采面开采层相对瓦斯涌出量，,m,3,/t,；,q,li,第,i,回采面邻近层相对瓦斯涌出量，,m,3,/t,；,A,i,第,i,回采面平均日产煤量，,t,；,m,采区内掘进巷道条数；,Q,1j,第,j,掘进巷道煤壁绝对瓦斯涌出量，,m,3,/min,；,Q,2j,第,j,掘进巷道落煤绝对瓦斯涌出量，,m,3,/min,；,A,c,采区平均日产煤量，,t,。,33,式中,q,矿,矿井相对瓦斯涌出量，,m,3,/t；,k,已采采区采空区瓦斯涌出系数，1.101.25；,n,生产采区个数；,q,ci,第,i,生产采区相对瓦斯涌出量，,m,3,/t；,A,ci,第,i,生产采区平均日产煤量，,t；,（,6,）矿井相对瓦斯涌出量,34,3,瓦斯地质数学模型预测方法,3.1,对以往预测方法的讨论,瓦斯含量法,经过近二、三十年不断的研究发展，目前已基本上达到实用化阶段。在实际应用中还存在的问题是：,（,1,）由于瓦斯含量法以煤层瓦斯含量作为预测的基础依据，因而对煤层瓦斯含量测定值的可靠性和含量点的分布及密度有较高的要求。如果预测区内只有很少或没有的瓦斯含量测定点（这种情况在地方煤矿是比较普遍的），那么第一步的瓦斯含量预测就不可靠，由此而进行的瓦斯涌出量预测，其精度将难以保证。,35,（,2,）国内外这类预测方法有近十种，由于各种方法所考虑的影响因素不同，或采用的物理模型不一样，因而形成了不同的计算公式。对同一矿井采用不同的瓦斯含量法，其预测结果相差较大。,矿山统计法,由于矿山统计法仅考虑瓦斯涌出量与开采深度一个因素之间的关系，故其适用范围受到一定的限制。对于地质条件简单的矿井，瓦斯涌出量的变化主要受开采深度的影响，预测结果可以满足生产要求。而在很多生产矿井，由于矿井地质条件的变化，瓦斯涌出量除了与开采深度有关以外，与其它地质因素也存在较密切的关系。在这种情况下，只考虑开采深度的预测方法将难以达到生产要求的预测精度。,36,3.2,基本原理,瓦斯地质数学模型法的基本原理是：通过瓦斯地质规律研究，分析瓦斯涌出量的变化规律，筛选影响瓦斯涌出量变化的主要地质因素；在此基础上，根据矿井已采地区的瓦斯涌出量实测资料和相关的地质资料，综合考虑包括开采深度在内的多种影响因素，采用一定的数学方法，建立预测瓦斯涌出量的多变量数学模型（预测方程）；利用所建立的数学模型，对矿井未采区域的瓦斯涌出量进行预测。,37,瓦斯地质数学模型法采用数量化理论,作为建模工具。数理化理论,是数量化理论的方法之一，用于解决从定性的或兼有定量的自变量出发对因变量的预测问题。,在瓦斯地质相关因素定量分析中，某些地质因素难以定量化，如煤层的顶、底板岩性，只是某种属性的描述，而没有量的概念，这类变量称为定性变量。某些定性变量有时是影响瓦斯涌出量变化的主要因素。另外，在实际应用中如果某些定量变量对瓦斯涌出的影响是趋势性的，将其转化为定性变量参加建立数学模型可能会得到更好的预测效果。,3.3,预测方法及步骤,（,1,）分析瓦斯涌出的影响因素,38,瓦斯涌出的影响因素,地质因素,开采因素,煤 层 围 岩,煤层组合特征,煤的变质程度,煤厚及其变化,地 质 构 造,水文地质条件,岩 浆 侵 入,煤层埋藏深度,开 采 顺 序,采 煤 方 法,顶板管理方式,瓦斯涌出的影响因素,39,（,2,）,建立预测瓦斯涌出量的数学模型,确定统计单元,选择变量,建立数学模型,式中，,x,u,(u=1,2,h),为第,u,个定量变量的数据，,(j,k),（,j =1,，,2,，,，,m;k=1,，,2,，,r,j,）为第,j,个项目第,k,各类目上的反应。,统计检验,（3）,未采区瓦斯涌出量预测,未采区域有完整的设计图,未采区只有规划设计，无完整的设计图,40,3.4,建模工具数量化理论,数量化理论,(,Theory of Quantification),是多元分析的一个分支，适合于同时处理定性数据和定量数据。数量化理论始于本世纪五十年代，最初应用于社会学研究，六十年代以后，随着电子计算机的广泛使用，才逐渐应用到自然科学领域中来，且应用日益增多。,研究数量化理论较早的数学家之一是日本的林知已夫，他最先采用了数量化理论这个名称，并先后提出了数理化理论、直至等一系列方法。,3.4.1,数量化理论的方法原理,数量化理论用于根据定性变量,(,可以兼有定量变量,),对定量的基准变量,(,因变量,),进行预测的问题。在数量化理论中，定性变量称为项目，定性变量的各种不同的取“值”称为类目。,41,（,1,）反应及反应矩阵,设因变量受,m,个项目,x,1,，,x,2,，,，,x,m,的影响，第一个项目,x,1,有,r,1,个类目,c,11,，,c,12,，,，,c,1 r1,，,第二个项目,x,2,有,r,2,个类目,c,21,，,c,22,，,，,c,2 r2,，,，,样 品,号,基准变量,x,1,x,2,x,m,c,11,c,12,c,1 r1,c,21,c,22,c,2 r2,c,m1,c,m2,c,m rm,1,2,3,n,y,1,y,2,y,3,y,n,1,(1,1),1,(1,2),1,(1,r,1,),2,(1,1),2,(1,2),2,(1,r,1,),3,(1,1),3,(1,2),3,(1,r,1,),n,(1,1),n,(1,2),n,(1,r,1,),1,(2,1),1,(2,2),1,(2,r,2,),2,(2,1),2,(2,2),2,(2,r,2,),3,(2,1),3,(2,2),3,(2,r,2,),n,(2,1),n,(2,2),n,(2,r,2,),1,(m,1),1,(m,2),1,(m,r,m,),2,(m,1),2,(m,2),2,(m,r,m,),3,(m,1),3,(m,2),3,(m,r,m,),n,(m,1),n,(m,2),n,(m,r,m,),个类目。如果有,n,个样品，其测定结果可表示为下列项目、类目反应表。,项目、类目反应表,第,m,个项目,x,m,有,r,m,个类目,c,m1,，,c,m2,，,，,c,m rm,，,m,个项目共有,42,表中,y,i,是因变量,y,在第,i,个样品中的测定值，,i,(j,k) (i=1,2,，,n,；,j=1,2,m,；,k=1,2,r,j,),称为,j,项目之,k,类目在第,i,个样品中的反应，其值按下式确定:,由元素,i,(j,k),构成的,n,p,阶矩阵,X,称为反应矩阵,43,这是由于任一样品在每个项目中只有一个类目的反应是1，其余类目的 反应皆为0之故。,（,2,）数量化理论的数学模型,假定因变量与各项目、类目的反应间遵从下列线性模型：,反应,i,(j,k),有个重要性质，即对每个固定的,i,和,j，,有,44,为了使方程的预测值尽量接近实测值，可以根据最小二乘法原理寻求系数,b,jk,，,即使,达到最小值。为此，求,q,对,b,uv,的偏导数，并令偏导数等于0，得到,设,是使,q,达到最小值的,b,jk,，,则,应满足上式，即,45,如果用矩阵形式来表示，上式可写成,式中,以上方程称为正规方程(组),从正规方程解出,后，便得到下面的预测,式中,(j,k),表示任一样品在,j,项目,k,类目上的反应，当取得一样品时，,便可由其反应,(j,k)，,利用上式算出,就是因变量,y,的预测值。,方程：,46,（,3,）正规方程的特点和解法,可以证明，上述正规方程具有如下特点：,正规方程的系数矩阵,X,X,是对称的。,正规方程系数矩阵的各列及右端常数项的列中，对应各项目的元素之和皆相等。, 正规方程中最多有,个方程是线性无关的，因此,其系数矩阵,X,X,是不满秩的，其秩,最多是,，方,程有无穷多组解。, 正规方程的任意一组解皆使,达到同一最小值。,47,对于正规方程的任一组解,因此，解正规方程时，,可以对每个,j=2,m,，,删去第,j,项目第一,预测值,都是一样的。,类目的方程，并取,，,以使删除后的方程组成为满秩的，故,可唯一地解出其余的,，这样得到的解不失一般性，且确使,q,达,到最小。,（,4,）预测方程的显著性检验及预测精度,预测方程的显著性检验：,48,或复相关系数,估计预测精度：,（,5,）项目对预测的贡献,衡量各项目对预测贡献大小的指标有三个：偏相关系数、方差比和范围。,49,偏相关系数,为了考虑每个项目单独对因变量的贡献，将,看作是第,i,个样品中第,j,个项目的定量数据。,令,式中,50,这样，,m+1,个变量（包括因变量）之间的样本相关矩阵为：,其中,以,r,uv,表示,R,的逆矩阵,R,-1,中的元素，则因变量,y,与第,u,个项目间的样本偏相关系数为：,51,方差比,范围,（,6,）选择项目的方法, 增加项目法,该方法每一次增加一个项目，选取能使剩余均方,为最小的项目。式中,u,为进入预测方程的项目数。,52, 减少项目法,该方法先取全部,m,个项目，使用数量化理论，然后逐次减少一个项目，减少的项目按偏相关系数不显著的标准来确定。,为此，计算偏相关系数,y,j,构成的统计量,3.4.2,兼有定性和定量自变量的数理化理论,设自变量中有,h,个是定量变量，它们在第,i,个样品中的数据为,x,i,(u) (u=1,2,h; i=1,2,n),。,有,m,个定性变量，即,m,个项目，其中第,j,个项目有,r,j,个类目,它们在第,j,个样品中的反应是,(j,k) (j=1,2, m ; k=1,2,r,j,; i=1,2,n),。,因变量的数据为,y,i,(i=1,2,n),。,53,假定因变量与各定量自变量及项目、类目的反应间遵从如下线性模型：,式中,b,u,(u=1,2,h)、b,jk,(j=1,2,m；k=1,2, r,j,),是未知系数，,i,(i=1,2, n),是随机误差。类似只有定性自变量的推导,过程，可以得到,b,u,和,b,jk,的最小二乘估计,满足正规方程,54,式中,正规方程可在,的条件下求解，而得到预测方程,55,可按前述同样的方法计算,F,统计量、复相关系数和偏相关系数，用来衡量预测效果和各变量在预测中的贡献。只是在计算样本相关矩阵,R,时，对于定性变量用,对于定量变量则用,x,i,(u) (u=1,2,h ; i=1,2, n ),即可。,56,3.4,应用实例,我们对某矿采用瓦斯地质数学模型法进行了瓦斯涌出量预测。选定,17,个回采工作面作为已知统计单元。根据瓦斯地质分析得出的规律并经过统计检验，选择了煤层底板标高、顶板砂岩比、煤层厚度、地质构造等,4,个与瓦斯涌出有密切关系的因素作为自变量，建立了如下瓦斯地质数学模型：,57,式中,瓦斯涌出量预测值，,m,3,/t,；,煤层底板标高，,m,（,定量变量）；,砂岩比项目之,0.4,类目；,煤层厚度项目之,35,m,类目；,煤层厚度项目之,5,m,类目；,地质构造项目之地垒类目；,地质构造项目之地堑类目。,58,经统计检验，预测方程的,F,统计量为,114.28,，在,0.01,水平下显著；计算了因变量对每个自变量的偏相关系数，其,t,统计量均在,0.01,水平下显著，说明所建立的数学模型有实际意义。,利用所建立的数学模型，对深部,28,个设计工作面进行了瓦斯涌出量预测，预测结果从略。预测工作完成后，预测区,25031,、,24061,二个工作面已先后回采完毕，利用其实测的瓦斯涌出量值对预测结果进行了实践性检验。检验结果，二个工作面的平均预测误差为,13.75%,，说明预测精度较高。,59,3.5,瓦斯地质数学模型预测软件,采用瓦斯地质数学模型法进行矿井瓦斯涌出量预测，需要进行大量的计算工作， 为此，我们采用,VC+6.0,语言开发了瓦斯地质数学模型法预测软件。,软件的主对话框由两个组合框和,6,个按钮组成。两个组合框用来初始化变量和输出运行结果。,6,个按钮的主要功能如下如下：,(1),提示按钮,(4),详细资料按钮,(2),输入或修改数据文件按钮,(5),预测位置单元按钮,(3),运行数据文件按钮,(6),退出按钮,60,瓦斯地质数学模型软件主对话框,3.5,瓦斯地质数学模型预测软件,61,按下预测位置单元按钮，程序将进入预测对话框，在预测对话框里可以对未知样品的因变量进行预测。,预 测 对 话 框,3.5,瓦斯地质数学模型预测软件,62,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,预测系统的研究目标,（,1,）实现数据提取、数学建模、预测过程和预测结果的数字化、自动化和可视化，提高瓦斯涌出量预测工作的质量和效率。,（,2,）将基于,GIS,的空间数据管理与常规的属性数据管理相结合，实现钻孔、瓦斯含量、瓦斯压力和瓦斯涌出量等瓦斯地质数据的实时查询和图形显示。,（,3,）将瓦斯涌出量预测同瓦斯地质图的编制结合起来，为矿井瓦斯地质信息的管理和分析提供数字化工作平台。,63,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,开发平台,以,.net2003,为开发平台 以,MapObjects,为开发组件,基于,GIS,的二次开发,64,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,数据的导入和导出,65,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,地物的添加和属性的编辑,鼠标拾取添加,对话框添加,批量添加,66,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,变量的选择与设置,67,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,统计单元数据的提取,68,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,建立瓦斯涌出量预测模型,69,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,未采区域瓦斯涌出量预测,70,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,绘制瓦斯涌出量预测图,71,3.6,基于,GIS,的瓦斯涌出量可视化预测系统研究,绘制瓦斯涌出量预测图,72,