(高教社杯数学建模论文)煤矿瓦斯和煤尘的监测与控制.doc

D题 煤矿瓦斯和煤尘的监测与控制 海南软件职业技术学院 09级数学建模 某参赛生 摘要 本文由绝对涌出量相对涌出量的定义，以及对原题附表附图的分析，并结合监测数据建立模型，又用数学关系式清楚的表达模型中各种数据之间的映射关系。模型包括如下几个主要部分1.绝对涌出量和相对涌出量的算法；2.风速与瓦斯浓度、风速和煤尘浓度之间的映射关系；3.串联在同一“U”形道里工作面和回风巷之间的速度的映射关系；4.危险程度的算法，危险级别的拟定。 在模型求解1里，结合工具软件“Matlab65”用计算机计算出每一次监测数据对应的绝对涌出量，30天来每一天的相对涌出量。把计算结果结合煤矿安全规程第一百三十三条规定得出该矿属于高瓦斯矿。 在模型求解2里，由危险程度算法模型，给出了判断该矿安全级别的两个评判方案方案1.从长期的监测数据中看该矿的危险度；方案2.从该矿最危险时（断电的临界瓦斯浓度）的危险度。并用计算机编程分别算出两种方法的计算结果均显示该矿长期处于安全生产的状态。 在模型求解3里，把最优风速定义为在满足安全生产、高效生产（不发生断电）的基础上对工人健康危险最小的风速定为最优风速。通过计算把最优风速定为工作面Ⅰ1.8m/s，工作面Ⅱ2.35m/s，局部通风机207.6m3/min，回风巷Ⅰ1.62m/s，回风巷Ⅱ2.35m/s，总回风巷4.1360m/s。同时考虑到在实际生产中可根据工人对新鲜风的需求会有所变动。模型的解给出最优风速可调节的范围为工作面Ⅰ1.84.0m/s，工作面Ⅱ2.354.0m/s，局部通风机207.6400m3/min，回风巷Ⅰ1.623.6m/s，回风巷Ⅱ2.354.0m/s，总回风巷8m/s。 关键词 瓦斯相对涌出量，瓦斯绝对涌出量，危险度商数（危商），危险级别，最优风速，最优风范围。 1.问题重述 煤矿安全生产是我国目前亟待解决的问题之一，做好井下瓦斯和煤尘的监测与控制是实现安全生产的关键环节（见附件1）。 瓦斯是一种无毒、无色、无味的可燃气体，其主要成分是甲烷，在矿井中它通常从煤岩裂缝中涌出。瓦斯爆炸需要三个条件空气中瓦斯达到一定的浓度；足够的氧气；一定温度的引火源。 煤尘是在煤炭开采过程中产生的可燃性粉尘。煤尘爆炸必须具备三个条件煤尘本身具有爆炸性；煤尘悬浮于空气中并达到一定的浓度；存在引爆的高温热源。试验表明，一般情况下煤尘的爆炸浓度是30～ 2000g/m3，而当矿井空气中瓦斯浓度增加时，会使煤尘爆炸下限降低，结果如附表1所示。 国家煤矿安全规程给出了煤矿预防瓦斯爆炸的措施和操作规程，以及相应的专业标准 见附件2。规程要求煤矿必须安装完善的通风系统和瓦斯自动监控系统，所有的采煤工作面、掘进面和回风巷都要安装甲烷传感器，每个传感器都与地面控制中心相连，当井下瓦斯浓度超标时，控制中心将自动切断电源，停止采煤作业，人员撤离采煤现场。具体内容见附件2的第二章和第三章。 附图1是有两个采煤工作面和一个掘进工作面的矿井通风系统示意图，请你结合附表2的监测数据，按照煤矿开采的实际情况研究下列问题 （1）根据煤矿安全规程第一百三十三条的分类标准 见附件2，鉴别该矿是属于“低瓦斯矿井”还是“高瓦斯矿井”。 （2）根据煤矿安全规程第一百六十八条的规定，并参照附表1，判断该煤矿不安全的程度（即发生爆炸事故的可能性）有多大 （3）为了保障安全生产，利用两个可控风门调节各采煤工作面的风量，通过一个局部通风机和风筒实现掘进巷的通风（见下面的注）。根据附图1所示各井巷风量的分流情况、对各井巷中风速的要求（见煤矿安全规程第一百零一条），以及瓦斯和煤尘等因素的影响，确定该煤矿所需要的最佳（总）通风量，以及两个采煤工作面所需要的风量和局部通风机的额定风量（实际中，井巷可能会出现漏风现象）。 注 掘进巷需要安装局部通风机，其额定风量一般为150400 m3/min。局部通风机所在的巷道中至少需要有15的余裕风量（新鲜风）才能保证风在巷道中的正常流动，否则可能会出现负压导致乏风逆流，即局部通风机将乏风吸入并送至掘进工作面。 名词解释 （1）采煤工作面矿井中进行开采的煤壁 采煤现场。 （2）掘进巷用爆破或机械等方法开凿出的地下巷道，用以准备新的采煤区和采煤工作面。 （3）掘进工作面掘进巷尽头的开掘现场。 （4）新鲜风不含瓦斯和煤尘等有害物质的风流。 （5）乏风含有一定浓度的瓦斯和煤尘等有害物质的风流。 2.问题分析 该问题是个生产优化问题，在求解过程中会用到大量计算，计算可以交给计算机。对于计算过程中用到的变量之间的映射关系要根据相关定义进行演绎，并结合计算机对相关数据的拟合来搞定。于是所建立的模型要搞清楚下列变量之间的映射关系 1.瓦斯绝对涌出量和风速、巷道截面、瓦斯浓度的映射关系。 2.瓦斯相对涌出量和绝对涌出量、矿煤量的映射关系。 3.由于漏风和瓦斯涌出等原因，串联在同一“U”形通风巷道里的工作面和回风巷的风速不一定相等。所以要搞定串联在同一“U”形通风巷道里的工作面和回风巷的风速映射，以及总回风巷的风速和支回风巷风速的映射。 4.煤尘浓度和风速的映射。 5.煤尘爆炸下限和空气中瓦斯浓度的映射。 6.危险度和爆炸物浓度、爆炸下限的映射。 原题中要求判断该矿属于高瓦斯还是低瓦斯矿，则可以由以上的映射关系结合煤矿安全规程第一百三十三条对低瓦斯矿井和高瓦斯矿井的判别条件（低瓦斯矿井矿井相对瓦斯涌出量小于或等于10m3 /t且矿井绝对瓦斯涌出量小于或等于40m3 /min；高瓦斯矿井矿井相对瓦斯涌出量大于10m3 /t或矿井绝对瓦斯涌出量大于40m3 /min）来确定。对于求解该矿的不安全程度以及最优风速等问题可以由以上的映射关系结合煤矿安全规程的相关规定来求解。 3. 模型假设及符号说明 3.1 模型假设 1.漏风的多少和风巷中风量成正比不会出现漏风率偶尔多偶尔少的现象。 2.忽略少量漏风对巷道中瓦斯浓度的影响。漏风会使计算出的瓦斯涌出量小于实际值，当漏风量相对于巷道风流量很小时可以忽略。 3.总回风巷里的瓦斯和煤尘完全来自其它支巷道（在模型建立时会对此假设的合理性给出证明）。 3.2 符号说明 n--------天数 n∈（1，2，330），多数情况下n做脚码用； i--------i为1、2、3、4、5、6分别代表工作面I、工作面II、掘进工作面、回风巷I、回风巷II、总回风巷,多数情况下i做脚码用； Q -------产煤量（t）； RJ ------绝对涌出量m3/min； RX-------相对涌出量； v -------风速； S -------巷道中的风筒的横截面积； Pg------ 空气中瓦斯浓度，单位体积百分比； Pc------ 空气中煤尘浓度g/m3； δ-------空气中有瓦斯时的煤尘爆炸下限g/m3； δg------瓦斯的爆炸下限，一般为5～16； δc------煤尘的爆炸下限，一般为30～50g/m3； K--------瓦斯参与使煤尘爆炸下限降低的降低系数； D--------整体的危险程度； Dg-------瓦斯的爆炸危险程度； Dc-------煤尘的爆炸危险程度； 注释声明 a. 本文所有的计算和编程代码全是工具软件“Matlab65”,不敢保证在其它版本的Matlab中能运行。 b. 下文用到的命令“bianliang”代码都是指输入D2006之附表2里的数据的Matlab程序。凡用到“bianliang”这个命令的程序代码均要有此“bianliang.m”文件为基础。“bianliang”的源代码见本文“附码一”。 c. 为了和原题中的符号统一，我们用δ表示爆炸下限，又由于所用工具软件不支持希腊字母做变量，故在程序代码中用英文字母“U”代替“δ”。 4.模型建立与必要的模型检验 4.1绝对涌出量建模。由瓦斯绝对涌出量的定义可知 瓦斯绝对涌出量0.60瓦斯浓度风速风巷断面 4.2瓦斯浓度与风速建模。由式1导出瓦斯浓度和风速的关系 瓦斯的绝对涌出量RJ在长期行为下是个常数。同一巷道里的截面S也是个常数，所以上试可以简化为 4.3相对涌出量建模。由相对涌出量的定义可知 某日瓦斯相对涌出量1440（一天的分钟数）该天绝对涌出量该天产煤量 4.4工作面风速和对应回风巷风速建模。考虑同一个U型通风巷里的工作面和回风巷的风速。由于他们串联在U型风巷里，而U型风巷在实际操作中可能会设置漏风洞，同时涌出瓦斯气体也会对采煤工作面和回风巷的风速造成影响，由于漏风率不变（假设1），所以串联在同一个U型通风巷里的风速会存在简单的正比关系。 即 1．计算r1。 代码命令及结果 bianliang; r1 meanv4./v1 r1 0.8964 为了方便计算把meanr1最近似值0.90。于是得到 v40.9v1 （4） 检验式4的拟合度，程序代码 bianliang; clc; CC11-abs1-v4./0.9*v1; string0[最大拟合度 ,num2str maxCC1];dispstring0; string0[最小拟合度 ,num2str minCC1];dispstring0; string0[平均拟合度 ,num2str meanCC1];dispstring0; 运行结果 最大拟合度 0.99952 最小拟合度 0.91775 平均拟合度 0.97862 平均拟合度为0.97862。证明式（4）合理可用。 二．计算r2。代码命令和运行结果 bianliang; r2 meanv5./v2 r2 1.0101 为了方便计算，在不会造成较大误差的情况下r2最近似值1，即得 v3v4 （5） 检验式5的拟合度代码 bianliang; CC11-abs1-v5./v2; string0[最大拟合度 ,num2str maxCC1];dispstring0; string0[最小拟合度 ,num2str minCC1];dispstring0; string0[平均拟合度 ,num2str meanCC1];dispstring0; 运行结果 最大拟合度 1 最小拟合度 0.98077 平均拟合度 0.98995 平均拟合度为0.98995。证明式5可用。 4.5总回风巷风速建模。单位时间内，总回风巷气流量回风巷Ⅰ，回风巷Ⅱ和掘进工作面气流量之和。即 检验式（6）的拟合度，代码命令 bianliang; CC11-abs5*v6./4*v3v4v5-1; string0[最大拟合度 ,num2str maxCC1];dispstring0; string0[最小拟合度 ,num2str minCC1];dispstring0; string0[平均拟合度 ,num2str meanCC1];dispstring0; 运行结果 最大拟合度 0.99962 最小拟合度 0.99123 平均拟合度 0.99581 平均拟合度为0.99581。证明式6可用。 4.6总回风巷瓦斯浓度建模。忽略总回风巷土层涌出的少量瓦斯气体，总回风巷瓦斯浓度等于回风巷Ⅰ、回风巷Ⅱ和掘进工作面气流混合在一起的瓦斯浓度。即是 检验式（7）的拟合度代码 bianliang; pLv3.*Pg3v4.*Pg4v5.*Pg5./v3v4v5; CC11-abs1-pL./Pg6; string0[最大拟合度 ,num2str maxCC1];dispstring0; string0[最小拟合度 ,num2str minCC1];dispstring0; string0[平均拟合度 ,num2str meanCC1];dispstring0; 运行结果 最大拟合度 0.99944 最小拟合度 0.95729 平均拟合度 0.98582 平均拟合度为0.98582。证明式7可用。同时也证明了假设3成立。 4.7总回风巷煤尘浓度建模。总回风巷里的煤尘来自回风巷Ⅰ，回风巷Ⅱ和掘进工作面，那么可建立求总回风巷煤尘浓度的模型 数学式 检验式（8）程序代码 bianliang; PLv3.*Pc3v4.*Pc4v5.*Pc5./v3v4v5; CC11-abs1-PL./Pc6; string0[最大拟合度 ,num2str maxCC1];dispstring0; string0[最小拟合度 ,num2str minCC1];dispstring0; string0[平均拟合度 ,num2str meanCC1];dispstring0; 运行结果 最大拟合度 0.98977 最小拟合度 0.92094 平均拟合度 0.95743 平均拟合度为0.95743。证明式8可用。 4.8危险度商数建模。用空气中爆炸物的浓度与该物爆炸下限的比值做为危险测量危险度商数的标准，并且用不同测量地点或不同爆炸物之间最大的爆炸危险度商数做为该矿的危险度商数。从而可建立下列关系式“危险度商数爆炸物浓度爆炸下限100” 。 于是 瓦斯危商瓦斯浓度瓦斯爆炸下限100 DgPg/δg （9） 煤尘危商煤尘浓度有瓦斯时煤尘的爆炸下限100 DcPc/δ （10） 危商Max瓦斯危商，煤尘危商 DmaxDg,Dc 由附表1分析空气中有瓦斯时煤尘的危险度远大于瓦斯的危险度，所以由上式可得 DPc/δ （11） 危险级别标准D≥100 时爆炸， 90≤D＜100定为极高危险， 80 ≤D＜90定为高危险，70 ≤D＜80为比较危险，60 ≤D＜70为比较安全，D＜60很安全。由此标准做如下“危险级别对照表” D的取值范围（） 小于60 6070 7080 8090 90100 大于等于100 危险度级别 很安全 比较安全 比较危险 高危险 极高危险 爆炸 4.9煤尘的爆炸下限与瓦斯浓度映射关系建模。由原题附表1，空气中煤尘的爆炸下限和空气中瓦斯浓度有个对应关系。 空气中瓦 浓度 0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4 煤尘爆炸下限浓度g/m3 3050 22.537.5 1525 10.517.5 6.512.5 4.57.5 35 2.53.5 1.52.5 参照煤矿安全规程第一百六十八条的规定，瓦斯浓度在达到1.5时巷道断电。所以为了模型的实用性，可以把上表中瓦斯浓度大于1.5的数据剔除，把剩下的有效数据用计算机进行拟合，出于人身生命安全的考虑，在拟合模型的时候爆炸下限的取值范围里取最小值（如3050时取30，22.537.5时取22.5）。经过多种拟合结果的比较得知，1次拟合的结果较为理想，即不会造成较大误差又很实用。用计算机拟合的代码和运行结果如下 Pg[0,0.5,1.0,1.5]; U[30,22.5,15,10.5]; polyfitPg,U,1 ans -13.2000 29.4000 为便于计算，小数点后保留两位有效数字，即得δ关于Pg的映射 δ-13.2Pg29.4（12） 把理论值和实际值进行分画图比较，代码如下 Pg[0,0.5,1.0,1.5]; U[30,22.5,15,10.5]; scatterPg,U,r*; hold on x[00.011.5]; y-13.2*x 29.4; plotx,y,. 运行结果得到图1，红色“*”点是附表1给出的实际爆炸下限。蓝线为理论值曲线。从图1可以看出理论值和实际值很吻合。式12可用。 4.10空气中煤尘浓度与风速映射关系建模。假设工作面I、工作面II、掘进工作面、回风巷I、回风巷II里的煤尘浓度和风速遵守同样的映射关系总回风巷里煤尘浓度遵守式（8），故在此不考虑。那么就从这5个监测点监测到的数据进行数据拟合得到模型。先对监测数据进行分析 以速度为自变量，经验表明，在同样误差情况下，自变量区间越小，数据模型越难拟合。甚至集中在小区间的数据点会让整个模型的建立起到负面作用。用计算机计算工作面I、工作面II、掘进工作面、回风巷I、回风巷II在30内检测到的风速各自的变化区间，命令代码如下 maxvi-minvi为监测点i的风速区间 bianliang; [maxv1-minv1,maxv2-minv2,maxv3-minv3,maxv4-minv4,maxv5-minv5] 运行结果 0.3900 0.4400 0.0800 0.0900 0.3500 由计算结果可以看出工作面II和回风巷II的变化区间分别只有0.08和0.09。这个风速区间相对于别的三个监测点太小。甚至可以认为这个变化区间是由风巷中人流绝对涌出量的波动和监测器误差造成的。但是这小小的风速区间里却有和别的监测点相同数目的数据点。这些挤在一起的数据点会使整体数据拟合出的模型误差更大。所以工作面2和回风巷2的数据不参与模型拟合，但遵守拟合模型。先画出i1,3,4时的vi与Pci的散点图，命令代码 bianlinag; x[v1,v3, v4]; y[Pc1, Pc3 ,Pc4]; scatterx,y,. 运行结果得图2。 由图2可以看出空气中煤尘浓度随巷道风速的增大而增大，在变化区间里增长趋势类似一条直线。于是可建立模型 计算α、β的值，命令代码 bianlinag; x[v1,v3, v4]; y[Pc1, Pc3 ,Pc4]; polyfitx,y,1 运行结果 ans 1.5980 3.9986 小数点后保留两位有效数字，取近似值即得α1.60，β4.00。于是有 Pc1.6v4 13 检验式（13）的拟合度，程序代码 bianliang; v[v1,v2,v3,v4,v5]; Pc[Pc1,Pc2,Pc3,Pc4,Pc5]; Pc01.6*v4; CC11-absPc./Pc0-1; string0[最大拟合度 ,num2str maxCC1];dispstring0; string0[最小拟合度 ,num2str minCC1];dispstring0; string0[平均拟合度 ,num2str meanCC1];dispstring0; 运行结果 最大拟合度 0.99973 最小拟合度 0.92896 平均拟合度 0.97399 从拟合度检验结果来看，式（13）拟合度很高。证明式（13）可用。 至此所需变量的映射关系已经搞清，并清楚地用数学式表达出来。同时也对于数据拟合出来的映射关系给出了拟合度证明。 5.模型求解 5.1 判断该矿属于低瓦斯矿还是高瓦斯矿。 5.1.1计算该矿的绝对涌出量。由式（8）总回风巷瓦斯浓度模型分析可以认为总回风巷的瓦斯的绝对涌出量为该矿的绝对涌出量。再依式（1）计算总回风巷的绝对涌出量，程序代码如下 bianliang; RJ660*5*v6.*Pg6*0.01; [aa,bb]maxRJ6; [cc,dd]minRJ6; eemeanRJ6; barRJ6; string0[最大绝对涌出量 ,num2straa,， 属于第,num2strbb,次检测];dispstring0; string0[最小绝对涌出量 ,num2strcc,， 属于第,num2strdd,次检测];dispstring0; string0[平均绝对涌出量 ,num2stree];dispstring0; string0[每次监测详细结果如图3];dispstring0; 运行结果 最大绝对涌出量 10.9752， 属于第69次检测 最小绝对涌出量 8.6832， 属于第86次检测 平均绝对涌出量 9.7536 每次监测详细结果如图3 5.1.2 计算相对涌出量。参照式（3）可得计算30天的平均相对涌出量公式为 求第n天的绝对涌出量，用该天早、中、晚绝对涌出量的平均值，再结合式（1）可得 依式（3）、（14）、（15）编写计算相对涌出量的程序代码 bianliang; for nn130 RJnnnv63*nn*Pg63*nnv63*nn-1*Pg63*nn-1v63*nn-2*Pg63*nn-2; 截面5乘以时间60s乘以百分比单位0.01除以监测次数3刚好等于1.故上式是第nn天的绝对涌出量 RXnnn 1440*RJnnn/Qnn; end RX1440*sumRJn/sumQ; [nn1,dat1]maxRXn;[nn2,dat2]minRXn; BarRJn; string0[本月最大相对涌出量 ,num2strnn1,,在第,num2strdat1,日;];dispstring0; string0[本月最小相对涌出量 ,num2strnn2,,在第,num2strdat2,日;];dispstring0; string0[本月平均相对涌出量 ,num2strRX,;];dispstring0; string0[详细结果如图4];dispstring0; 运行结果 本月最大相对涌出量 25.3991,在第23日; 本月最小相对涌出量 21.3655,在第12日; 本月平均相对涌出量 23.1832; 每天的详细结果如图4 平均相对涌出量为23.183210m3/t。依煤矿安全规程第一百三十三条高瓦斯矿井矿井相对瓦斯涌出量大于10m3/t或矿井绝对涌出量大于40m3/min的标准。该矿属于高瓦斯矿。 5.2评算该矿的危险程度。 方案一用每次监测到的瓦斯浓度结合式（12）算出监测时的煤尘爆炸下限，再结合式（11）用当时监测到的煤尘浓度比上当时的爆炸下限。以此计算结果来分析该矿长期工作中的危险度情况。用计算机编程运算，并输出最大危险度和平均危险度。 程序代码 bianliang; Pg[Pg1,Pg2,Pg3,Pg4,Pg5,Pg6]; Pc[Pc1,Pc2,Pc3,Pc4,Pc5,Pc6]; U-13.2*Pg 29.4; DPc./U; aa100*maxD; bb100*meanD; string0[本月最大危险度 ,num2straa,, 平均危险度 ,num2strbb,.];dispstring0; 运行结果 本月最大危险度 52.5174, 平均危险度 37.082. 参照危险级别对照表 D的取值范围（） 小于60 6070 7080 8090 90100 大于等于100 危险度级别 很安全 比较安全 比较危险 高危险 极高危险 爆炸 由方案一的计算结果可以看出该矿长期处于很安全的状态。 方案二把断电时的危险度做为该矿的危险度。瓦斯浓度高可以不怕，怕的是瓦斯浓度高的时候工作面和风巷不断电。依煤矿安全规程第一百六十八条的规定，高瓦斯矿井的采煤工作面断电浓度Pg1.5，复电浓度Pg1.0,复电浓度Pg1.96;v31.02。 由第一百零四条规定，依规定可得出如下约束条件 0.25≤v1,v2,v3≤4； v6≤8；结合式（6）可得4/5（v1v2v3）≤8； 综合上述解不等式的结果可得 1.80≤v1≤4；2.35≤v2≤4; 1.02≤v3≤4；v1v2v3≤10； 在安全生产和生产效率的基础上还要考虑工人健康。由于空气中煤尘对人体危害很大，而空气不煤尘的浓度是随风速的增大而增大的，所以在安全生产和不影响生产效率的情况下取最小风速为最优风速，于是最优风速可以定为 v11.8m/s；v22.35m/s；v31.02m/s；v40.9v11.62；v5v22.35m/s； v60.8v1v2v34.1360m/s； 由v3计算局部通风机的风量。由于局部通风机所在的巷道中至少需要有15的余裕风量（新鲜风）才能保证风在巷道中的正常流动，所以局部通风机的风量不大于所在巷道通风量的85，又由于v3已经是取的最小值了，所以就让局部通风机的风量等于所在巷道通风量的85。即 风量（局部）85v3S0.851.0243.46（m3/s）207.6m3/min。 207.6m3/min在局部通风机额定风量150400 m3/min范围内。 于是最终确定最优风速为工作面Ⅰ1.8m/s，工作面Ⅱ2.35m/s，局部通风机207.6m3/min，回风巷Ⅰ1.862m/s，回风巷Ⅱ2.35m/s，总回风巷4.14m/s。 考虑到在实际生产中可根据工人对新鲜风的需求来灵活调节最优风速，调动范围为工作面Ⅰ1.84.0m/s，工作面Ⅱ2.354.0m/s，局部通风机207.6400m3/min，回风巷Ⅰ1.623.6m/s，回风巷Ⅱ2.354.0m/s，总回风巷4.148m/s。 6.模型的进一步讨论 本模型有一定的灵活性，只要有相关统计数据，进行小量的改进工作就可以胜任相同煤矿类型的更大规模的生产优化问题。 在实际生产会遇到更复杂的煤矿生产优化问题。在不改变煤矿类型只扩大生产规模的情况下该模型的变量之间关系式的一些数值可能有所变动，但变量之间的映射模型不会变。在上文中建立了13个关系式，其中的变量映射模型有较广的实用性，如果换成另外一个类似的矿井，只需要将相关数据输入计算机用同样的方法编程计算，就可以方便的得出精确的关系式，进而可供计算使用。 本模型也有局限性，无法用于不同类型的煤矿生产的优化问题，例如对瓦斯突出矿井的生产优化问题就无能为了。 附码一 该附码为“bianling.m”的代码。只需在matlab命令中输入命令“bianliang”,附表2里的速度m/s、瓦斯浓度1、煤尘浓度分别以数组vi、Pgi、Pci存入计算机内存，然后便可随意调用相应的数据。 v1[2.4000 2.5600 2.2400 2.2700 2.4100 2.4000 2.2700 2.2200 2.2900 2.3100 2.3700 2.4000 2.3400 2.4300 2.2700 2.2900 2.3900 2.3100 2.3600 2.4000 2.5200 2.4900 2.2100 2.2700 2.3900 2.3700 2.2700 2.3600 2.3600 2.3600 2.2900 2.3100 2.3700 2.3700 2.3600 2.3600 2.2700 2.3700 2.2100 2.4300 2.3600 2.3100 2.5200 2.2700 2.4600 2.3400 2.2100 2.2700 2.3100 2.3700 2.3400 2.4300 2.4100 2.3600 2.3900 2.3100 2.3600 2.4000 2.3400 2.2700 2.3600 2.3400 2.3100 2.2400 2.2900 2.3700 2.5600 2.2700 2.3100 2.3600 2.3600 2.4500 2.3900 2.2700 2.5200 2.3400 2.3600 2.3400 2.3100 2.3600 2.3100 2.4100 2.2500 2.2200 2.4000 2.6000 2.3400 2.3100 2.3700 2.4000]; 工作面1的一个月内所测的风速 v2[ 2.0700 2.0800 2.0900 2.0800 2.1300 2.1000 2.1000 2.0900 2.0900 2.0900 2.0800 2.1000 2.1100 2.1400 2.0800 2.1000 2.1100 2.1200 2.1300 2.0900 2.0700 2.1400 2.0900 2.0700 2.0600 2.1000 2.0800 2.1200 2.1300 2.1100 2.0700 2.0900 2.09