1、 D 题: 煤矿瓦斯和煤尘的监测与控制(海南软件职业技术学院 09 级数学建模 某参赛生)摘要本文由绝对涌出量相对涌出量的定义,以及对原题附表附图的分析,并结合监测数据建立模型,又用数学关系式清楚的表达模型中各种数据之间的映射关系。模型包括如下几个主要部分:1.绝对涌出量和相对涌出量的算法;2.风速与瓦斯浓度、风速和煤尘浓度之间的映射关系;3.串联在同一“U”形道里工作面和回风巷之间的速度的映射关系;4.危险程度的算法,危险级别的拟定。在模型求解 1 里,结合工具软件“Matlab65”用计算机计算出每一次监测数据对应的绝对涌出量,30 天来每一天的相对涌出量。把计算结果结合煤矿安全规程第一百
2、三十三条规定得出该矿属于高瓦斯矿。在模型求解 2 里,由危险程度算法模型,给出了判断该矿安全级别的两个评判方案:方案 1.从长期的监测数据中看该矿的危险度;方案 2.从该矿最危险时(断电的临界瓦斯浓度)的危险度。并用计算机编程分别算出两种方法的计算结果均显示该矿长期处于安全生产的状态。在模型求解 3 里,把最优风速定义为:在满足安全生产、高效生产(不发生断电)的基础上对工人健康危险最小的风速定为最优风速。通过计算把最优风速定为:工作面1.8m/s,工作面2.35m/s,局部通风机 207.6m3/min,回风巷1.62m/s,回风巷2.35m/s,总回风巷 4.1360m/s。同时考虑到在实际
3、生产中可根据工人对新鲜风的需求会有所变动。模型的解给出最优风速可调节的范围为:工作面1.84.0m/s ,工作面2.354.0m/s,局部通风机 207.6400m3/min,回风巷 1.623.6m/s ,回风巷2.354.0m/s,总回风巷 8m/s。关键词:瓦斯相对涌出量,瓦斯绝对涌出量,危险度商数(危商) ,危险级别,最优风速,最优风范围。11.问题重述煤矿安全生产是我国目前亟待解决的问题之一,做好井下瓦斯和煤尘的监测与控制是实现安全生产的关键环节(见附件 1) 。瓦斯是一种无毒、无色、无味的可燃气体,其主要成分是甲烷,在矿井中它通常从煤岩裂缝中涌出。瓦斯爆炸需要三个条件:空气中瓦斯达
4、到一定的浓度;足够的氧气;一定温度的引火源。煤尘是在煤炭开采过程中产生的可燃性粉尘。煤尘爆炸必须具备三个条件:煤尘本身具有爆炸性;煤尘悬浮于空气中并达到一定的浓度;存在引爆的高温热源。试验表明,一般情况下煤尘的爆炸浓度是 30 2000g/m3,而当矿井空气中瓦斯浓度增加时,会使煤尘爆炸下限降低,结果如附表 1 所示。国家煤矿安全规程给出了煤矿预防瓦斯爆炸的措施和操作规程,以及相应的专业标准 (见附件 2)。规程要求煤矿必须安装完善的通风系统和瓦斯自动监控系统,所有的采煤工作面、掘进面和回风巷都要安装甲烷传感器,每个传感器都与地面控制中心相连,当井下瓦斯浓度超标时,控制中心将自动切断电源,停止
5、采煤作业,人员撤离采煤现场。具体内容见附件 2 的第二章和第三章。附图 1 是有两个采煤工作面和一个掘进工作面的矿井通风系统示意图,请你结合附表 2 的监测数据,按照煤矿开采的实际情况研究下列问题: (1)根据煤矿安全规程第一百三十三条的分类标准 (见附件 2),鉴别该矿是属于“低瓦斯矿井”还是“高瓦斯矿井” 。(2)根据煤矿安全规程第一百六十八条的规定,并参照附表 1,判断该煤矿不安全的程度(即发生爆炸事故的可能性)有多大? (3)为了保障安全生产,利用两个可控风门调节各采煤工作面的风量,通过一个局部通风机和风筒实现掘进巷的通风(见下面的注) 。根据附图 1 所示各井巷风量的分流情况、对各井
6、巷中风速的要求(见煤矿安全规程第一百零一条) ,以及瓦斯和煤尘等因素的影响,确定该煤矿所需要的最佳(总)通风量,以及两个采煤工作面所需要的风量和局部通风机的额定风量(实际中,井巷可能会出现漏风现象) 。注 掘进巷需要安装局部通风机,其额定风量一般为 150400 m3/min。局部通风机所在的巷道中至少需要有 15%的余裕风量(新鲜风)才能保证风在巷道中的正常流动,否则可能会出现负压导致乏风逆流,即局部通风机将乏风吸入并送至掘进工作面。名词解释(1)采煤工作面:矿井中进行开采的煤壁 (采煤现场)。(2)掘进巷:用爆破或机械等方法开凿出的地下巷道,用以准备新的采煤区和采煤工作面。(3)掘进工作面
7、:掘进巷尽头的开掘现场。(4)新鲜风:不含瓦斯和煤尘等有害物质的风流。2(5)乏风:含有一定浓度的瓦斯和煤尘等有害物质的风流。2.问题分析该问题是个生产优化问题,在求解过程中会用到大量计算,计算可以交给计算机。对于计算过程中用到的变量之间的映射关系要根据相关定义进行演绎,并结合计算机对相关数据的拟合来搞定。于是所建立的模型要搞清楚下列变量之间的映射关系:1.瓦斯绝对涌出量和风速、巷道截面、瓦斯浓度的映射关系。2.瓦斯相对涌出量和绝对涌出量、矿煤量的映射关系。3.由于漏风和瓦斯涌出等原因,串联在同一“U” 形通风巷道里的工作面和回风巷的风速不一定相等。所以要搞定串联在同一“U” 形通风巷道里的工
8、作面和回风巷的风速映射,以及总回风巷的风速和支回风巷风速的映射。4.煤尘浓度和风速的映射。5.煤尘爆炸下限和空气中瓦斯浓度的映射。6.危险度和爆炸物浓度、爆炸下限的映射。原题中要求判断该矿属于高瓦斯还是低瓦斯矿,则可以由以上的映射关系结合煤矿安全规程第一百三十三条对低瓦斯矿井和高瓦斯矿井的判别条件(低瓦斯矿井:矿井相对瓦斯涌出量小于或等于 10m3 /t 且矿井绝对瓦斯涌出量小于或等于 40m3 /min;高瓦斯矿井:矿井相对瓦斯涌出量大于 10m3 /t 或矿井绝对瓦斯涌出量大于40m3 /min)来确定。对于求解该矿的不安全程度以及最优风速等问题可以由以上的映射关系结合煤矿安全规程的相关规
9、定来求解。3. 模型假设及符号说明3.1 模型假设1.漏风的多少和风巷中风量成正比不会出现漏风率偶尔多偶尔少的现象。2.忽略少量漏风对巷道中瓦斯浓度的影响。漏风会使计算出的瓦斯涌出量小于实际值,当漏风量相对于巷道风流量很小时可以忽略。3.总回风巷里的瓦斯和煤尘完全来自其它支巷道(在模型建立时会对此假设的合理性给出证明) 。3.2 符号说明n-天数 n(1,2,330) ,多数情况下 n 做脚码用;3i-i 为 1、2、3、4、5、6 分别代表工作面 I、工作面 II、掘进工作面、回风巷 I、回风巷 II、总回风巷,多数情况下 i 做脚码用;Q -产煤量(t) ;RJ -绝对涌出量(m 3/mi
10、n);RX-相对涌出量;v -风速;S -巷道中的风筒的横截面积;Pg- 空气中瓦斯浓度,单位体积百分比;Pc- 空气中煤尘浓度(g/m 3);-空气中有瓦斯时的煤尘爆炸下限(g/m 3);g-瓦斯的爆炸下限,一般为 516%;c-煤尘的爆炸下限,一般为 3050(g/m 3);K-瓦斯参与使煤尘爆炸下限降低的降低系数;D-整体的危险程度;Dg-瓦斯的爆炸危险程度;Dc-煤尘的爆炸危险程度;注释声明:a. 本文所有的计算和编程代码全是工具软件“Matlab65”,不敢保证在其它版本的 Matlab 中能运行。b. 下文用到的命令“bianliang”代码都是指输入D2006之附表 2 里的数据
11、的 Matlab 程序。凡用到“bianliang”这个命令的程序代码均要有此“bianliang.m”文件为基础。“bianliang”的源代码见本文“附码一” 。c. 为了和原题中的符号统一,我们用 表示爆炸下限,又由于所用工具软件不支持希腊字母做变量,故在程序代码中用英文字母“U”代替“”。4.模型建立与必要的模型检验4.1 绝对涌出量建模。由瓦斯绝对涌出量的定义可知:4瓦斯绝对涌出量=0.60瓦斯浓度风速风巷断面 16.0iii SvPgRJ4.2 瓦斯浓度与风速建模。由式(1)导出瓦斯浓度和风速的关系 iii SvJg6.0瓦斯的绝对涌出量 RJ 在长期行为下是个常数。同一巷道里的截
12、面 S 也是个常数,所以上试可以简化为: 21iivPg4.3 相对涌出量建模。由相对涌出量的定义可知某日瓦斯相对涌出量=1440(一天的分钟数)该天绝对涌出量该天产煤量 3140nnQRJX4.4 工作面风速和对应回风巷风速建模。考虑同一个 U 型通风巷里的工作面和回风巷的风速。由于他们串联在 U 型风巷里,而 U 型风巷在实际操作中可能会设置漏风洞,同时涌出瓦斯气体也会对采煤工作面和回风巷的风速造成影响,由于漏风率不变(假设1) ,所以串联在同一个 U 型通风巷里的风速会存在简单的正比关系。即:1计算 r1。代码命令及结果:bianliang; r1= mean(v4./v1)r1 = 0
13、.8964为了方便计算把 mean(r1)最近似值 0.90。于是得到:v4=0.9v1 (4)检验式 4 的拟合度,程序代码:bianliang;clc;CC1=1-abs(1-v4./(0.9*v1);是 比 例 系 数rvrv142514255string0=最大拟合度 =,num2str( max(CC1);disp(string0);string0=最小拟合度 =,num2str( min(CC1);disp(string0);string0=平均拟合度 =,num2str( mean(CC1);disp(string0);运行结果:最大拟合度 =0.99952最小拟合度 =0.91
14、775平均拟合度 =0.97862平均拟合度为 0.97862。证明式(4)合理可用。二计算 r2。代码命令和运行结果:bianliang; r2= mean(v5./v2)r2 = 1.0101为了方便计算,在不会造成较大误差的情况下 r2 最近似值 1,即得:v3=v4 ( 5) 检验式 5 的拟合度代码:bianliang;CC1=1-abs(1-v5./v2);string0=最大拟合度 =,num2str( max(CC1);disp(string0);string0=最小拟合度 =,num2str( min(CC1);disp(string0);string0=平均拟合度 =,nu
15、m2str( mean(CC1);disp(string0);运行结果:最大拟合度 =1最小拟合度 =0.98077平均拟合度 =0.98995平均拟合度为 0.98995。证明式 5 可用。4.5 总回风巷风速建模。单位时间内,总回风巷气流量=回风巷,回风巷和掘进工作面气流量之和。即:665436 Svv检验式(6)的拟合度,代码命令:bianliang;CC1=1-abs(5*v6)./(4*(v3+v4+v5)-1);6string0=最大拟合度 =,num2str( max(CC1);disp(string0);string0=最小拟合度 =,num2str( min(CC1);dis
16、p(string0);string0=平均拟合度 =,num2str( mean(CC1);disp(string0);运行结果:最大拟合度 =0.99962最小拟合度 =0.99123平均拟合度 =0.99581平均拟合度为 0.99581。证明式 6 可用。4.6 总回风巷瓦斯浓度建模。忽略总回风巷土层涌出的少量瓦斯气体,总回风巷瓦斯浓度等于回风巷、回风巷和掘进工作面气流混合在一起的瓦斯浓度。即是: 气 流 量掘 进回 风回 风 瓦 斯 量掘 进回 风回 风总 回 风 巷 瓦 斯 浓 度 217543536 vPgPgv检验式(7)的拟合度代码:bianliang;pL=(v3.*Pg3+
17、v4.*Pg4+v5.*Pg5)./(v3+v4+v5);CC1=1-abs(1-pL./Pg6);string0=最大拟合度 =,num2str( max(CC1);disp(string0);string0=最小拟合度 =,num2str( min(CC1);disp(string0);string0=平均拟合度 =,num2str( mean(CC1);disp(string0);运行结果:最大拟合度 =0.99944最小拟合度 =0.95729平均拟合度 =0.98582平均拟合度为 0.98582。证明式 7 可用。同时也证明了假设 3 成立。4.7 总回风巷煤尘浓度建模。总回风巷里
18、的煤尘来自回风巷,回风巷和掘进工作面,那么可建立求总回风巷煤尘浓度的模型:气 流 量掘 进回 风回 风 煤 尘 量掘 进回 风回 风总 回 风 巷 煤 尘 浓 度 217数学式: 8543536 vPcPcv检验式(8)程序代码:bianliang;PL=(v3.*Pc3+v4.*Pc4+v5.*Pc5)./(v3+v4+v5);CC1=1-abs(1-PL./Pc6);string0=最大拟合度 =,num2str( max(CC1);disp(string0);string0=最小拟合度 =,num2str( min(CC1);disp(string0);string0=平均拟合度 =,n
19、um2str( mean(CC1);disp(string0);运行结果:最大拟合度 =0.98977最小拟合度 =0.92094平均拟合度 =0.95743平均拟合度为 0.95743。证明式 8 可用。4.8 危险度商数建模。用空气中爆炸物的浓度与该物爆炸下限的比值做为危险测量危险度商数的标准,并且用不同测量地点或不同爆炸物之间最大的爆炸危险度商数做为该矿的危险度商数。从而可建立下列关系式“危险度商数=爆炸物浓度爆炸下限100%” 。于是:瓦斯危商=瓦斯浓度 瓦斯爆炸下限100% Dg=Pg/g (9)煤尘危商=煤尘浓度 有瓦斯时煤尘的爆炸下限100% Dc=Pc/ (10)危商=Max(
20、瓦斯危商,煤尘危商) D=max(Dg,Dc) 由附表 1 分析空气中有瓦斯时煤尘的危险度远大于瓦斯的危险度,所以由上式可得:D=Pc/ (11)危险级别标准:D100% 时爆炸, 90%D100%定为极高危险, 80% D90% 定为高危险,70% D80%为比较危险, 60% D70%为比较安全,D60% 很安全。由此标准做如下“危险级别对照表 ” :D 的取值范围(%) 小于 60 6070 7080 8090 90100 大于等于 100危险度级别 很安全 比较安全 比较危险 高危险 极高危险 爆炸4.9 煤尘的爆炸下限与瓦斯浓度映射关系建模。由原题附表 1,空气中煤尘的爆炸下限8和空
21、气中瓦斯浓度有个对应关系。空气中瓦浓度 (%) 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 次
22、拟合的结果较为理想,即不会造成较大误差又很实用。用计算机拟合的代码和运行结果如下:Pg=0,0.5,1.0,1.5;U=30,22.5,15,10.5; polyfit(Pg,U,1)ans = -13.2000 29.4000为便于计算,小数点后保留两位有效数字,即得 关于 Pg 的映射:=-13.2Pg+29.4(12)把理论值和实际值进行分画图比较,代码如下:Pg=0,0.5,1.0,1.5;U=30,22.5,15,10.5; scatter(Pg,U,r*);hold onx=0:0.01:1.5;y=-13.2*x+ 29.4;plot(x,y,.)运行结果得到图 1,红色“*”点
23、是附表 1 给出的实际爆炸下限。蓝线为理论值曲线。从图 1 可以看出理论值和实际值很吻合。式 12 可用。94.10 空气中煤尘浓度与风速映射关系建模。假设工作面 I、工作面 II、掘进工作面、回风巷 I、回风巷 II 里的煤尘浓度和风速遵守同样的映射关系总回风巷里煤尘浓度遵守式(8) ,故在此不考虑。那么就从这 5 个监测点监测到的数据进行数据拟合得到模型。先对监测数据进行分析:以速度为自变量,经验表明,在同样误差情况下,自变量区间越小,数据模型越难拟合。甚至集中在小区间的数据点会让整个模型的建立起到负面作用。用计算机计算工作面 I、工作面 II、掘进工作面、回风巷 I、回风巷 II 在 3
24、0 内检测到的风速各自的变化区间,命令代码如下:%max(vi)-min(vi)为监测点 i 的风速区间bianliang;max(v1)-min(v1),max(v2)-min(v2),max(v3)-min(v3),max(v4)-min(v4),max(v5)-min(v5)运行结果:0.3900 0.4400 0.0800 0.0900 0.3500 由计算结果可以看出工作面 II 和回风巷 II 的变化区间分别只有 0.08 和 0.09。这个风速区间相对于别的三个监测点太小。甚至可以认为这个变化区间是由风巷中人流绝对涌出量的波动和监测器误差造成的。但是这小小的风速区间里却有和别的监
25、测点相同数目的数据点。这些挤在一起的数据点会使整体数据拟合出的模型误差更大。所以工作面 2 和回风巷 2 的数据不参与模型拟合,但遵守拟合模型。先画出 i=1,3,4 时的vi 与 Pci 的散点图,命令代码:bianlinag;x=v1,v3, v4;10y=Pc1, Pc3 ,Pc4;scatter(x,y,.)运行结果:得图 2。由图 2 可以看出空气中煤尘浓度随巷道风速的增大而增大,在变化区间里增长趋势类似一条直线。于是可建立模型: 为 常 数,vPc计算 、 的值,命令代码:bianlinag;x=v1,v3, v4;y=Pc1, Pc3 ,Pc4;polyfit(x,y,1)运行结
26、果:ans = 1.5980 3.9986小数点后保留两位有效数字,取近似值即得 =1.60,=4.00。于是有:Pc=1.6v+4 (13)检验式(13)的拟合度,程序代码:bianliang;v=v1,v2,v3,v4,v5;11Pc=Pc1,Pc2,Pc3,Pc4,Pc5;Pc0=1.6*v+4;CC1=1-abs(Pc./Pc0-1);string0=最大拟合度 =,num2str( max(CC1);disp(string0);string0=最小拟合度 =,num2str( min(CC1);disp(string0);string0=平均拟合度 =,num2str( mean(C
27、C1);disp(string0);运行结果:最大拟合度 =0.99973最小拟合度 =0.92896平均拟合度 =0.97399从拟合度检验结果来看,式(13)拟合度很高。证明式(13)可用。至此所需变量的映射关系已经搞清,并清楚地用数学式表达出来。同时也对于数据拟合出来的映射关系给出了拟合度证明。5.模型求解5.1 判断该矿属于低瓦斯矿还是高瓦斯矿。5.1.1 计算该矿的绝对涌出量。由式(8)总回风巷瓦斯浓度模型分析可以认为总回风巷的瓦斯的绝对涌出量为该矿的绝对涌出量。再依式(1)计算总回风巷的绝对涌出量,程序代码如下:bianliang;RJ6=60*5*v6.*Pg6*0.01; aa
28、,bb=max(RJ6); cc,dd=min(RJ6); ee=mean(RJ6);bar(RJ6);string0=最大绝对涌出量 =,num2str(aa),, 属于第,num2str(bb),次检测;disp(string0); string0=最小绝对涌出量 =,num2str(cc),, 属于第,num2str(dd),次检测;disp(string0);string0=平均绝对涌出量 =,num2str(ee);disp(string0); string0=每次监测详细结果如图 3;disp(string0); 运行结果:最大绝对涌出量 =10.9752, 属于第 69 次检测最
29、小绝对涌出量 =8.6832, 属于第 86 次检测平均绝对涌出量 =9.7536每次监测详细结果如图 3125.1.2 计算相对涌出量。参照式(3)可得计算 30 天的平均相对涌出量公式为:141401nQRJX求第 n 天的绝对涌出量,用该天早、中、晚绝对涌出量的平均值,再结合式(1)可得:1536.032 ni iiSvPgRJ依式(3) 、 (14) 、 (15)编写计算相对涌出量的程序代码:bianliang;for nn=1:30RJn(nn)=(v6(3*nn)*Pg6(3*nn)+v6(3*nn-1)*Pg6(3*nn-1)+v6(3*nn-2)*Pg6(3*nn-2);%截面
30、 5 乘以时间 60s 乘以百分比单位 0.01 除以监测次数 3 刚好等于 1.故上式是第 nn 天的绝对涌出量RXn(nn)= 1440*RJn(nn)/Q(nn);endRX=1440*sum(RJn)/sum(Q); nn1,dat1=max(RXn);nn2,dat2=min(RXn);Bar(RJn);string0=本月最大相对涌出量 =,num2str(nn1),在第,num2str(dat1),日;disp(string0);string0=本月最小相对涌出量 =,num2str(nn2),在第,num2str(dat2),日;disp(string0); 13string0
31、=本月平均相对涌出量 =,num2str(RX),;disp(string0);string0=详细结果如图 4;disp(string0); 运行结果:本月最大相对涌出量 =25.3991,在第 23 日;本月最小相对涌出量 =21.3655,在第 12 日;本月平均相对涌出量 =23.1832;每天的详细结果如图 4平均相对涌出量为 23.183210m3/t。依煤矿安全规程第一百三十三条(高瓦斯矿井:矿井相对瓦斯涌出量大于 10m3/t 或矿井绝对涌出量大于 40m3/min)的标准。该矿属于高瓦斯矿。5.2 评算该矿的危险程度。方案一:用每次监测到的瓦斯浓度结合式(12)算出监测时的煤
32、尘爆炸下限,再结合式(11)用当时监测到的煤尘浓度比上当时的爆炸下限。以此计算结果来分析该矿长期工作中的危险度情况。用计算机编程运算,并输出最大危险度和平均危险度。程序代码:bianliang;Pg=Pg1,Pg2,Pg3,Pg4,Pg5,Pg6;Pc=Pc1,Pc2,Pc3,Pc4,Pc5,Pc6;U=-13.2*Pg+ 29.4;D=Pc./U;aa=100*max(D); bb=100*mean(D);string0=本月最大危险度 =,num2str(aa),%, 平均危险度 =,num2str(bb),%.;disp(string0);运行结果:本月最大危险度 =52.5174%,
33、平均危险度 =37.082%.参照危险级别对照表:14D 的取值范围(%) 小于 60 6070 7080 8090 90100 大于等于 100危险度级别 很安全 比较安全 比较危险 高危险 极高危险 爆炸由方案一的计算结果可以看出该矿长期处于很安全的状态。方案二:把断电时的危险度做为该矿的危险度。瓦斯浓度高可以不怕,怕的是瓦斯浓度高的时候工作面和风巷不断电。依煤矿安全规程第一百六十八条的规定,高瓦斯矿井的采煤工作面断电浓度 Pg=1.5%,复电浓度 Pg=1.0%,复电浓度 Pg1.76;v21.96;v31.02。由第一百零四条规定,依规定可得出如下约束条件:0.25v 1,v2,v34
34、;v68;结合式(6)可得:4/5(v 1+v2+v3)8;综合上述解不等式的结果可得:1.80v 14;2.35v 24; 1.02v 34;v 1+v2+v310;在安全生产和生产效率的基础上还要考虑工人健康。由于空气中煤尘对人体危害很大,而空气不煤尘的浓度是随风速的增大而增大的,所以在安全生产和不影响生产效率的情况下取最小风速为最优风速,于是最优风速可以定为:v1=1.8m/s;v 2=2.35m/s; v3=1.02m/s;v 4=0.9v1=1.62;v 5=v2=2.35m/s;v6=0.8(v1+v2+v3)=4.1360m/s;由 v3 计算局部通风机的风量。由于局部通风机所在
35、的巷道中至少需要有 15%的余裕风量(新鲜风)才能保证风在巷道中的正常流动,所以局部通风机的风量不大于所在巷道通风量的 85%,又由于 v3 已经是取的最小值了,所以就让局部通风机的风量等于所在巷道通风量的 85%。即:19风量(局部)=85%v 3S=0.851.024=3.46( m3/s)=207.6m 3/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。考虑到在实际生产中可
36、根据工人对新鲜风的需求来灵活调节最优风速,调动范围为:工作面1.84.0m/s,工作面2.354.0m/s,局部通风机 207.6400m3/min,回风巷1.623.6m/s,回风巷2.354.0m/s,总回风巷 4.148m/s。6.模型的进一步讨论本模型有一定的灵活性,只要有相关统计数据,进行小量的改进工作就可以胜任相同煤矿类型的更大规模的生产优化问题。在实际生产会遇到更复杂的煤矿生产优化问题。在不改变煤矿类型只扩大生产规模的情况下该模型的变量之间关系式的一些数值可能有所变动,但变量之间的映射模型不会变。在上文中建立了 13 个关系式,其中的变量映射模型有较广的实用性,如果换成另外一个类
37、似的矿井,只需要将相关数据输入计算机用同样的方法编程计算,就可以方便的得出精确的关系式,进而可供计算使用。本模型也有局限性,无法用于不同类型的煤矿生产的优化问题,例如对瓦斯突出矿井的生产优化问题就无能为了。附码一%该附码为“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.3
38、700 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.
39、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 202.3100 2.3700 2.4000; %工作面 1 的一个月内所测的风速v2= 2.0700 2.0800 2.090
40、0 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.0900 2.0900 2.1000 2.1000 2.1300 2.0900 2.1200 2.0700 2.0800 2.0900 2.1000 2.0900 2.1100 2.12
41、00 2.0900 2.1000 2.1000 2.0700 2.1100 2.1000 2.1100 2.0900 2.1200 2.0700 2.1100 2.1200 2.0700 2.0800 2.0900 2.1100 2.1000 2.1100 2.1000 2.1000 2.1100 2.1300 2.1000 2.0900 2.0900 2.0800 2.1100 2.1000 2.1100 2.1200 2.1000 2.1000 2.1000 2.1200 2.1000 2.1200 2.0900 2.1100 2.0700 2.1200 2.1100 2.1000 2.1
42、100 2.1100; %工作面 2 的一个月内所测的风速v3= 2.2200 2.2900 2.4000 2.2100 2.2900 2.2900 2.2900 2.3100 2.1500 2.2100 2.1100 2.4400 2.2700 2.3600 2.2100 2.2700 2.1000 2.1800 2.2100 2.2900 2.2400 2.2700 2.2200 2.2200 2.2400 2.2900 2.2900 2.2700 2.2300 2.2400 2.1400 2.1800 2.0900 2.2200 2.1800 2.1800 2.1400 2.4400 2
43、.2100 2.1400 2.3100 2.2200 2.2900 2.2200 2.2200 2.1900 2.1500 2.1800 2.1500 2.1100 2.3700 2.3600 2.1400 2.3100 2.2500 2.2900 2.2700 2.2300 2.3400 2.2100 2.2100 2.1400 2.3100 2.3700 2.3600 2.2400 2.3600 2.3100 2.2900 2.1500 2.1000 2.2200 2.2700 2.3100 2.2300 2.3500 2.3500 2.3400 2.2200 2.2500 2.1400
44、2.0900 2.3400 2.2200 2.1400 2.1700 2.2700 2.2300 2.2400 2.1100; %掘进工作面的一个月内所测的风速v4= 2.1300 2.2700 2.0900 2.0800 2.2300 2.2100 2.0400 2.1100 2.0200 2.1100 2.1700 2.1100 2.0800 2.1700 1.9600 1.9300 2.0200 2.0000 2.0900 2.0200 2.2700 2.2100 1.9100 2.0800 2.1300 2.1500 2.0600 2.0400 2.1300 2.1000 2.1100
45、 2.1400 2.1100 1.9700 2.1000 2.2700 2.1400 2.1400 1.9700 2.1800 2.0900 2.1100 2.2100 2.0200 2.2100 2.1500 2.0000 2.0000 2.1300 2.0400 2.0600 2.2200 2.2100 2.1800 2.0900 2.0800 2.0800 2.2200 2.1100 2.0400 2.0900 2.0800 2.0800 2.0000 2.0900 2.1100 2.2100 2.0400 2.2500 2.1100 2.0900 2.1900 2.1100 2.000
46、0 2.2700 2.0000 2.2500 2.0900 2.0400 2.1500 2.0400 2.1400 2.0000 2.0000 2.1300 2.3500 2.0000 2.1700 2.1300 2.1500; %回风巷 1 的一个月内所测的风速v5= 2.0900 2.1000 2.1100 2.1100 2.1600 2.1200 2.1100 2.1100 2.1200 2.1200 2.1200 2.1200 2.1400 2.1600 2.1100 2.1200 2.1300 2.1300 2.1600 2.1200 2.0700 2.1600 2.1000 2.0900 2.0800 2.1200 2.1200 2.1300 2.1400 2.1400 2.0900 2.1100 2.1100 2.1100 2.1200 2.1200 2.1400 2.1100 2.1600 2.0900 2.0900 2.1100 2.1400 2.1200 2.1200 2.1400 2.1200 2.1300 2.1200 2.0700 2.