煤矿瓦斯和煤尘的监测和控制案例.docx

《煤矿瓦斯和煤尘的监测和控制案例.docx》由会员分享，可在线阅读，更多相关《煤矿瓦斯和煤尘的监测和控制案例.docx（25页珍藏版）》请在三一文库上搜索。

1、 最新资料推荐最新资料推荐煤矿瓦斯和煤尘的监测和控制案例暑期数学建模培训之煤矿瓦斯和煤尘的监测和控制案例组员：机电工程学院自动化083何亮亮理学院数学聂翠白信息学院通信郑文明指导老师：焉化彪摘要；本文研究的是有关煤矿瓦斯和煤尘的监测和控制的问题。煤矿安全生产时我国亟待解决的问题之一，做好井下瓦斯和煤尘的监测和控制是实现安全生产的关键环节。因此，本文就所给的信息和相应的数据，按照煤矿开采的实际情况解决了以下三个难题：针对问 题一，根据相关的数据信息， 我们建立了 基本的数 学模型，并应用MATLAB软件进行对附表2中总回风巷的风速及瓦 斯进行求日平均值， 再根据煤矿安全规章 第一百三十三条的 分

2、类标准， 求出相对瓦斯量和绝对瓦斯量浓度。因此，我们利用 软件月 平均相 对瓦 斯量和 月 平均 绝对瓦 斯量分别为：9.75603m /min，23.20193m /min.根据判断高瓦斯矿井还是低瓦斯矿井的标准得到：该煤矿属于高瓦斯矿井针对问题二，首先，我们对附表一进行分析， 引入降低浓度 系数k.在不同空气中瓦斯浓度情况下， 建立一般情况下煤尘爆炸下 限浓度与不同瓦斯对应的煤尘爆炸下限浓度的关系：k=m ；再次，利用MATLAB软件，对改进数据进行二次拟合，得到瓦斯浓度降低浓度系数曲线方程。最后，利用分别对各个监测点进行绘制的图形结合 煤矿规程 第一百三十八条规定进行分析， 从而，我们可

3、以得到该煤矿不安全 度几乎上是接近0。针对问题三，根据对问题的分析，我们建立了线性目标规划模型。利用Lin go软件或是MATLAB软件求的最优结果即：该煤矿所需要的最优总通风量为：zQ =25.14785s /m3;与及两个采煤工作面所需要的风量分别为：1Q=41V=2.257498*4=9.029992s/m3 ； 2Q =4606762.242V10.427048s /m3 ；同时，我们还得到局部通风机的额定风量和主 巷道所需要风量分别为：=59674375.05Q77V3.86975s/m3；=5138162.15Q33V5.69081S /m3；另外，根据计算结果我们得到附图1各个

4、巷道的风速值如下表所示：各巷道进风巷1( V1) 进风巷2 (V2)回风巷1 (V4) 回 风 巷2 (V5)主巷道(V3) 待 掘 进巷(V6 )道局 部 通风 机 巷道(乏风道)(V7)风(m/s) 速 2.257498 2.606762 1.8059982.985410 1.138162 0.9105294 0.9674375 最后， 本文就模型的 评价和推广展开了的叙述。关键字：MATLAB软件Lingo 软件降低浓度系数相对瓦斯浓度绝对瓦斯浓度 线性规划模型一.问题的重述 二. 问题的分析 现今，煤矿安全生产是我国目前亟待解决的问题之一，因此，做好井下瓦斯和煤尘的监测和控制是实现安全

5、生产的关键所在。其中，瓦斯爆炸要满足三个条件： .空气中瓦斯要达到一定的浓度。 .足够的氧气。 .一定温度的引火源。而煤尘爆炸的条件： .煤尘本身具有爆炸性。 .煤尘悬浮在空气中并达到一定的浓度。 .存在引爆的咼温热源。另外，国家煤矿安全规程 给出了煤矿预防瓦斯爆炸的措施 和操作过程， 与及相应的专业指标标准。根据要求煤矿必须安装完善的通风系统和瓦斯自动控制系统，并且，所有的采煤工作面、掘进面、 和回风巷都必须安装甲烷传感器。当井下浓度超标时，控制中心将自动切断电源， 停止采煤工作，工作人员都必须撤离采煤现场。因此，根据这几种相关的情况。要对煤矿瓦斯和煤尘的监测和控制有所措施和解决问题的相应

6、方案。因此，对所给问题一一展开分析。2.1 问题一的分析要判断一个煤矿是低瓦斯矿井还是高瓦斯 矿井，就必须要寻找一个指标来进行比较。根据煤矿安全规程 规定的第一百三十三条分类标准。目标就转向求解该矿井的相对瓦斯浓度和绝对瓦斯浓度。考虑到对象是这个矿井，所以，我们只需要使用附表2中总回风巷的有关数据根据所建立的基本数学模型即可求解出相对瓦斯 浓度和绝对瓦斯浓度，然后根据标准即可判断该矿井属于何种瓦斯 矿井。2.2 问题二的分析 对于问题二，要求我们根据相关数据信息， 求解并断定矿井不安全的程度。根据附表一和结合题目 开始所说：当矿井空气中瓦斯浓度增加时，煤尘爆炸浓度下限降低。于是，我们引进降低浓

7、度系数k,建立起一般煤尘爆炸下限浓度 和在空气中不同瓦斯浓度所对应的煤尘爆炸下限浓度之间的关系。这样，我们就可以通过MATLAB软件对其进行拟合得到降低浓 度系数空气中阿瓦斯浓度曲线方程。由此，在分别就各个监测点的数据信息进行绘制图形，按照煤5 / 24 最新资料推荐最新资料推荐矿安全 规程的第一百六十八条的标准，进行分析。从而解决问题二。2.3 问题三的分析 为了保障安全生产煤矿采用两个可控风门 调节采煤工作面的风量，通过一个局部通风机和风筒实现掘进巷的 通风。首先，根据附图的示意图和煤矿安全规程的第一百零一条规定各巷道的风速约束， 再者，考虑到瓦斯和煤尘因素的影响。针对问题三，我们应建立一

8、个线性规划模型，就以最佳总的通风量为目标，其他的约束和因素都是约束条件。另外，我们分析示意图可和，就以总的进风巷为准，它由进风巷1和进风巷2与及主巷道三部分组成。还有，要理解局部通风机所在的巷道中至少有 15%的新鲜风。这句话的意思是还有85%的胀风即乏风将被其吸入并送至道掘 进工作面， 从而，进入主巷道。根据这点，我们又能建立一个表达关系式作为约束条件。在实际的矿井生产煤的情况下， 我们又要考虑在各个工作面的 巷道中不均衡的通风系数、 在各个掘进工作面的巷道中不均匀的备 用通风系数与及矿井通风系数等。通过这些来建立模型的约束条件。最后，根据线性规划模型的格式应用 Lingo软件或是MATLA

9、B 软件进行求解，得到最优解结果。三.模型的假设1.假设瓦斯的主要成分是甲烷，在矿井它通常从煤岩裂缝中涌出。2 . 假设国家煤矿安全生产给出的有关规程以及相应专业标准的信息是准确无误的；3 . 假设一般情况下， 煤尘爆炸下限浓度为：3050g/m3; 4.假设掘进巷需要安装局部通风机，其额定风量一般为：150400 m3/mi n; 5.假设在每个采煤工作面的巷道中不均衡的备用通风系数一般为：1.151.2 ;在此，我们取1.15; 6 .假设在每个采煤掘进工作面的巷道中不均匀的备用通风系数一般为：1.52 ;在此，我们取1.5; 7 .假设矿井通风系数为：1.25; 这样，进风巷和回风巷就能

10、够建立起相应的关系式；8假设在空气中只有瓦斯和煤尘的影响，其他的因素不用考虑。9.假设不考虑因突发性事件的发生而造成瓦斯或是煤尘爆炸的情况。四.模型的有关规定4.1名词解释1 .采煤工作面:矿井中进行开米的煤壁（米煤现场）；2 .掘进巷:用爆破或是机械等方法开凿出的地下巷道，用以准备新的采煤9 / 24区和采煤工作面；3. 掘进工作面:掘进巷尽头的开凿现场；4. 新鲜风:不含瓦斯和煤尘等有害物质的风流； 5 .乏风:含有一定浓度的瓦斯和煤尘等有害物质的风流；6. 绝对瓦斯浓度：单位内涌出瓦斯体积；7 .相对瓦斯浓度：平均日产量煤的过程中 涌出的瓦斯量；4.2 符号规定Qi J : 表示第i个监

11、测点的每天的绝对瓦斯浓度 QJ:表示绝对瓦斯浓 度；iXQ :表示第i个监测点每天的相对瓦斯浓度；XQ :表示相对瓦斯浓度；Gd:表示日产量煤；Vi:表示各监测点的 风速；Wi:表示各个监测点的瓦斯含量；K:表示降低浓度系数；K1:表示在每个采煤工作面中不均衡的通风系数；K2 :表示在每个采煤掘进工作面中不均匀的备用通风系数；K3:表示矿井备用通风系数；zQ :表示最佳总的通风量；j V :表示各个巷道中的风速；(j=1 2 3 4 5 6 7 ) ; m : 表示煤尘在空气中不同瓦斯，表示一般煤尘的一般情况Si :表示在各采煤各个巷道的断面积iQ :表示各进风巷的风量。(i=1 2 3 )

12、wiQ :表示第i个巷道因瓦斯或是煤尘影响所得到 的总风量。五.模型的建立5.1问题一的模型建立1 . 相关名词解释 相对瓦斯浓度：表示是平均日产量煤的过程中涌出的瓦斯量；绝对瓦斯浓度：表示单位时间涌出的瓦斯体积。2 .判断高低瓦斯矿井的标准低瓦斯矿井：相对瓦斯涌出量小于或是等于 10min/3m且绝对瓦斯涌出量小于 或是等于40min/3m ;高瓦斯矿井：相对瓦斯涌出量大于10min/3m ;或是，绝对瓦斯涌出量大于40min/3m;3 . 模型的建立如下：相对瓦斯浓度：XQ =(QJ/ Gd)6024 (tm /3)绝对瓦斯浓度：QJ二V66W60/100 ;其中：V6:表示总回风巷每天的

13、风速；W6:表示总回风巷每天的瓦斯含量；一般表示方法：iQX= Qi J/ Gd; Qi J=ViWi60/100; 5.2问题二模型的建立5.2.1模型数据准备 首先，针对附表一的有关数据，根据全面的分析，我们在次引入K,得到相对应的表格：空气中瓦斯量0 0.5 1.0 降低浓度系数1 0.75 0.5 此表所对 应的数学表达式为：1. 5 2.0 0.35 2.5 0.15 3.0 0.1 3.5 0.07 4 0.05 0.25 m=K ;5.2.2根据MATLAB软件求解出二次拟合曲线。最终我们建立数据数学模型，分别对各个监测点进行图形绘制分析，得出相应的结果。最新资料推荐具体相关解答

14、步骤在模型的求解时。会一一展现在我们面前。最主要一点就是编写程序。5.3 问题三的模型建立 531所用符号规定 最佳总的通风量：Qz ；进风巷的通风量：Qi ; 第i个巷道因瓦斯或煤尘所造成的通风量：wiQ各个巷道的风速：Vi;采煤工作区各个巷道断面面积：Si 5.3.2 各自的量纲表示 Qz二Q321Q+Q+; wiQ =100*QJ*K1 wiQ =100* QJ*K2 ;其中：Qi=ViSi ； 5.3.3由于漏风阻风的影响我们假设回风巷1、回风巷2、掘进巷的巷道的风速分别为：V4、V5、V6;而进风巷1、进风巷2、主巷道的各自巷道的风 速可分别为：V1 、V2、V3;由于在回风巷当中会

15、有胀风所以，他们之间存在一定关系：Vi=K2*Vj ; 其中：i=1,2,3;j=4,5,6; 5.3.4在各个工作点的要求情况对于每一个进风巷的风量必须要大于或是等于由于瓦斯或是煤尘影响所产 生的通风量； 即:Vi*Si Qwi其内容:i=1,2,7; 5.3. 5 局部通风机的要求约束 V37V*5.80二536根据上述的数学数据准备，我们最终得到下面的线性规划模型：目标函数：min二 V1*4+V2*4+3V *5;约束条件：wiijiQSVVVVKVVVVVVVV=*85. 0*3425. 08625. 0625. 08625. 0;625. 0i677654321 六.模型的求解与分

16、析6.1问题一的求解根据我们建立的基本数学模型，我们应用MATLAB软件进行了计算， 对附表2的总回风巷的数据进行了统计分析， 我们得到了相对瓦斯 的日平均值和绝对瓦斯日平均值。如下表1所示：（程序代码见附录）天数相对量绝对量天数相对量绝对量天数相对量 1 2 3 4 5 6 7 8 9 24.1672 23.5064 22.820722.4671 22.6012 24.9985 24.0228 22.5982 23.3557 10.0190 9.827010.1267 9.6109 9.7541 10.2077 9.7092 9.4944 9.4883 10 11 12 1314 15 16

17、 17 18 22.6737 22.3671 21.3730 24.0119 23.3418 22.7717 22.9572 24.8190 24.2863 9.7365 9.5682 9.0242 10.2051 9.8230 9.3617 9.7090 10.4274 9.7988 19 20 21 22 23 24 25 26 27 22.2452 22.8800 23.449 22.7693 25.4057 22.1312 23.2615 23.3526 22.9423 绝对量天数相对量绝对量 通过上表：我们还可以求出月 的平均值分别为：9.7560 和 23.2019。所以，根据准则

18、，我们得知该矿井属于高瓦斯矿井。同样，如果是以最大值来进行比较我们也能得到同样的结果:由观察上表我们、。O、； 天数相对量绝对量天数相对量绝对量天数相对量绝对量 天数 相 11.3636 11.1546 11.40311.3539.51609.72409.7850 9.7402 10.3740 9.6056 9.6600 9.5930 9.8301 28 29 3022.4961 24.0229 23.02359.45159.7760 9.912912 3 4 5 6 7 8 11.5200 10. 7823 10.9499 11.0338 10.7470 11.697612.2484 11.

19、2257 4.7760 4.5076 4/8590。4.7200 4.5528 4.0265 4.9504 4.7164 9 10 11 12 13 14 15 169.0561 10.9078 11.0796 10.9149 11.7659 10.7976 11.7568 11.39583.6790 4.6964 4.7396 4.6085 5.0005 4.5440 4.8333 4.8195 17 1819 20 21 22 23 24 11.3063 11.6588 11.332710.7762 11.477310.9145 11 . 9486 10.5726 4.7502 4.704

20、0 4.8479 4.5799 4.79024.6690 4.8790 4.5888 25 26 27 28 29 30 12.14410.798对量 绝对量 3 7 8 1 4 . 7149 4.7872 4.8860 4.7701 4.9423 4.6492 其他表在此就不必再画表了，在附录当中都有，对应的每个监测点的情况；6.2 问题二的解答6.2.1首先，根据数学软件得到二 次拟合曲线图如下：程序代码见附录即可：622针对每个监测点我们都进行绘制图形进行分析研究:注明：所有程序代码都建见附录1 .工作面1 :2. 工作面2 :3. 掘进面：4. 回风巷1 :5. 回风巷2 :通过对这些

21、监测点的绘图情况，我们绘制了 有关空气甲烷浓度和煤尘浓度的曲线关系。并且，结合国家煤矿安全生产规程 的第一百六十八条规定， 从每个监测点的绘制图形来看， 每个监测点不会因为瓦斯或者说是 因为煤尘而爆炸， 根据附表1的帮助，才能够得以完成第二问。所以：我们可以下决定：该矿井的不安全程度几乎为0 ;6.3 问题三的解答过程及分析 首先，根据建立模型，我们对其进行简化求解：目标函数：min 二V5V4V4321+ ;约束条件：60/ )5479.15 . 1100(4V60/ )402.455.11100(460/)7113. 415.1100(4125.125.125.185.06625.0625

22、.08625.0625.07263524137654321二=VVVVVVVV VVVVVV根据简化的模型应用 Lingo软件，很快就可以求出最优解的答案。所以，根据计算，我们得到结果如下：总的通风量为：25.14785m3/s=1508.871 m3/mi n; 同时，两个采煤工作面的通 风量分别为：Q1=2.257498*4=9.029992m3/s=541.7992m3/mi n;Q2=2.606762*4=10.427048m3/s=625.62288m3/mi n;Q7=0,9674375*4=3.86975 m3/s=231.185 m3/min;还有， 经过计算， 其他巷道的风速

23、如下表所示：各巷道进风巷1 2进风巷回风巷1回风巷2主巷 道待掘进巷道局部通风机巷道（乏风道）0.9674375 风（m/s）速 2.257498 2.606762 1.805998 2.985410 1.1381620.9105294七. 模型的评价本文主要是建立了两个重要模型，但实际上就只有后面的线性规划模型，在第一个模型中。我们只考虑了 空气中只含有瓦斯和煤尘两种有害成分。然而在实际生后中，由于空气中还有大气压，而且压力和温度对矿井的日常产煤量有影响，更糟糕的是可能还会引起瓦斯爆炸或者 煤尘爆炸情况。这些都是没有进行考虑的问题。所以说，根据我们自己建立模型所计算出来的结果是有误差的， 但

24、在另一方， 本模型通俗易懂，实用性强。形成两个目标。然后分别计算， 得出所需要的 解答。而在模型二中我们用的线性规划模型， 简单实用，准确度较 强， 至于问题二主要体现编写程序上， 这也体现了数据模型的特 色。当然在问题二中不安全程度没有用一个确定的数值，应该讲这是文章唯一一个不好的地方。附录：一.问题一的程序代码：1.1求解相对和绝对瓦斯量代码程序fun ction helia ngOA二5.15 5.36 5.31;5.14 5.39 5.32;5.19 5.26 5.05;5.175.145.36;5.21 5.37 5.03;5.08 5.01 5.07;5.19 5.18 5.29;

25、5.345.015.13;5.175.285.19;5.185.23 5.21;5.15.175.08;5.06 5.13 5.28;5.16 5.36 5.1;5.13 5.22 5.16;5.325.115.26;5.25.065.07;5.1455.27;5.385.25.31;5.2 5.2 5.2;5.285.24 5.08;5.165.08 5.21;5.23 5.285.2;5.35 5.23 5.38;5.12 5.05 5.24;5.245.16 5.35;5.23 5.45.27;5.12 5.26 5.06;5.11 5.19 5.11;5.1 5.36 5.14;5.2

26、35.23 5.11;B=0.670.6 0.63;0.6 0.64 0.62;0.68 0.65 0.63;0.68 0.58 0.58;0.61 0.58 0.65;0.66 0.68 0.68;0.6 0.620.64;0.59 0.62 0.63;0.61 0.63 0.58;0.65 0.62 0.6; 0.620.610.64;0.60.58 0.57;0.68 0.62 0.66;0.570.690.64;0.55 0.62 0.62;0.640.65 0.61;0.68 0.7 0.65;0.6 0.630.62;0.68 0.57 0.58;0.660.59 0.62;0.6

27、20.66 0.62;0.64 0.61 0.61;0.59 0.68 0.68;0.67 0.62 0.58;0.57 0.620.65;0.6 0.61 0.6;0.610.64 0.66;0.63 0.60.61;0.710.540.63;0.65 0.6 0.66; C二597 602 639 616 610 588 582 605 585 620616 608 612 606 592 609 605 581 616 612 601 616 588 625 598 618617605586620;x=zeros(1,30);y=zeros(1,30);s1=zeros(1,30);s2

28、=zeros(1,30);for i=1:30for j=1:3s1(i)=A(i,j)+s1(i);s2(i)=B(i,j)+s2(i);ends1(i)=s1(i)/3;s2(i)=s2(i)/3;endfor i=1:3017 / 24x(i)=(s1(i)*s2(i)*5*60)/100)+x(i);y(i)=(x(i)*24*60)/C(i)+y(i);enda=sum(x)/30;b=sum(y)/30; disp('显示所有结果:')disp('如下所示:')disp('月平均绝对瓦斯涌出量和相对瓦斯涌出量:')a b disp(&

29、#39;日绝对瓦斯涌出量和日相对瓦斯涌出量分别为:')x y运行结果:1.2除总回风巷外其他的监测点相对和绝对的瓦斯含量情况代码程序1.2.1 工作面一 function heliang01 A=2.42.56 2.24;2.27 2.41 2.4;2.27 2.222.29;2.312.372.4;2.34 2.43 2.27;2.29 2.39 2.31;2.362.4 2.52;2.49 2.21最新资料推荐仃/ 242.27;2.39 2.37 2.27;2.36 2.36 2.36;2.29 2.31 2.37;2.372.36 2.36;2.27 2.37 2.21;2.4

30、3 2.36 2.31;2.52 2.27 2.46;2.34 2.21 2.27;2.31 2.37 2.34;2.43 2.41 2.36;2.39 2.312.36;2.4 2.34 2.27;2.36 2.34 2.31;2.24 2.29 2.37;2.562.27 2.31;2.362.36 2.45;2.39 2.27 2.52;2.342.362.34;2.312.36 2.31;2.412.25 2.22;2.4 2.6 2.34;2.312.372.4; B二0.71 0.62 0.66;0.62 0.63 0.66;0.78 0.71 0.66;0.67 0.71 0.6

31、2;0.61 0.64 0.69;0.72 0.64 0.69;0.7 0.690.65;0.64 0.73 0.66;0.59 0.68 0.3;0.650.67 0.67;0.7 0.650.69;0.64 0.65 0.66;0.720.71 0.76;0.6 0.69 0.63;0.63 0.70.67;0.710.710.7;0.660.68 0.69;0.650.66 0.65;0.71 0.66 0.69;0.64 0.66 0.66;0.69 0.7 0.66;0.68 0.70.65;0.6 0.75 0.7;0.66 0.62 0.64;0.63 0.70.64;0.66

32、0.680.7;0.7 0.68 0.72;0.7 0.7 0.68;0.68 0.61 0.73;0.62 0.630.72; C=597 602 639 616 610 588 582 605 585 620 616 608 612606 592 609 605 581 616 612 601 616 588 625 598 618 617 605 586 620;s1二zeros(1,30);s2二zeros(1,30);x=zeros(1,30);y=zeros(1,30);s1(i)=A(i,j)+s1(i);s1(i)=s1(i)/3;fori=1:30for i=1:30for

33、j=1:3s2(i)二B(i,j)+s2(i);ends2(i)=s2(i)/3;end x(i)=(s1(i)*s2(i)*5*60)/100)+x(i);y(i)=(x(i)*24*60)/C(i)+y(i);enda二sum(x)/30;b=sum(y)/30; disp(显示所有结果：)disp(如下所示：)disp(工作一月平均绝对瓦斯涌出量和相对瓦斯涌出量：)a b disp(工作面一日绝对瓦斯涌出量和日相对瓦斯涌出量分别为：)x y运行结果：1.2.2 工作面二、掘进面、回风一、回风二相关的程序代码和工作一得相似，在此不作详写。通过MATLAB的数学软件求得相应的结果分别如下：

34、工作面二：掘进面：回风1 :回风2 :二. 问题二的程序代码 2.1求解二次拟合曲线function heliang04 figure(1); hold on x=0 0.5 0.75 1.0 1.5 2.02.5 3.0 3.5 4.0;y=1 0.75 0.6 0.5 0.35 0.25 0.15 0.1 0.07 0.05;si=0:0.02:4;p=polyfit(x,y,2)ri=polyval(p,si)plot(x,y, ' ro' ,si,ri, ' -.k ' ) xlabel( 空气中瓦斯浓度：)ylabel( 降低浓度系数：)-0.5178

35、2 x +0.97621 ' ) hold off2.2根据各个监测点的数据情况，我们通过用MATLAB软件绘制出空气中甲烷的含量与煤尘浓度的曲线关系；工作面一 ；functionhelia ng05figure(1) hold on x=0:0.02:4;y3=0:0.01:50;x1仁 1.5;% 断电浓度 x12=1; % 报警浓度 x2=0.71 0.62 0.66 0.62 0.63 0.66 0.78 0.71 0.66 0.67 0.71 0.62 0.61 0.64 0.69 0.72 0.640.69 0.7 0.69 0.65 0.64 0.73 0.66.0.59

36、 0.68 0.73 0.650.67 0.67 0.7 0.65 0.69 0.64 0.65 0.66 0.720.72 0.76 0.6 0.690.69 0.63 0.7 0.67 0.71 0.71 0.7.0.660.68 0.69 0.65 0.660.65 0.71 0.66 0.69 0.64 0.66 0.66 0.69 0.7 0.66.0.680.7 0.65 0.6 0.75 0.7 0.66 0.62 0.64 0.63 0.7 0.64.0.660.68 0.66 0.7 0.68 0.72 0.7 0.7 0.68 0.68 0.61 0.73 0.62 0.6

37、30.72; y4=8 8.4 7.61 7.69 7.99 7.87 7.75 7.71 7.68 7.67 7.897.99 7.83 8.1 7.64 7.81 7.96 7.73 7.81 7.95 8.23 8.19 7.597.67 .7.82 8.08 7.65 7.95 7.85 7.86 7.88 7.74 8.1 7.867.8 7.81 7.8 7.77 7.59 7.96 8 7.75 8.08 7.75 8.17 7.9 7.56 7.647.81 8 7.94 8 8.05 7.88 7.9 7.83 .7.95 7.88 7.85 7.627.59 7.85 7.

38、67 7.52 7.84 7.98 8.26 7.7 7.9 7.87 7.96 8.14 7.77.77 8.12 8 7.71 7.73 7.99 7.9 7.9 7.98 7.82 7.58 8.17 8.42 7.887.91 8.05 7.9;y1=30*(0.073818*x.-0.52339*x+0.97907);y2=50*(0.073818*x.-0.52339*x+0.97907);plot(x,y1, ' r* ' ,x,y2, ' sk' ,x11,y3,x12,y3,x2,y4,' m+ );最新资料推荐legend（ meic

39、henbaozanongduzuixiaozhi ' , ' meichenbaozazuida zhi ' , ' duandiannongdu' , ' baojinnongdu ' ）; xlabel（ 空气中 甲烷浓度）;ylabel（煤尘所含浓度）;title（ 工作面1的各指标的分布情况）;hold off运行结果：工作面二：functionheliang06 hold on x=0:0.02:4;y1=0:0.01:50;x11 = 1.5;% 断电浓度 x12=1; % 报警浓度 x2=0.94 0.85 0.82 0.8

40、5 0.92 0.89 0.89 0.92 0.94 0.96 0.72 0.87 0.84 0.8 0.99 0.97 0.98 1 0.77 0.84 0.94 0.78 0.8.0.941.01 0.91 0.69 0.87 0.850.84 0.84 0.81 0.83 0.88 0.81 0.67 0.94 0.89 0.82 0.78 1.11 0.85 0.79 0.85 0.76 0.8 .0.89 0.78 0.98 1.03 0.94 0.91 0.840.93 0.93 0.74 0.78 0.94 0.74 0.86 0.84 0.96 0.91 0.92 0.83 0

41、.94 0.93 0.97 0.96 .0.92 0.93 0.83 0.84 0.94 0.99 0.810.87 0.8 0.8 0.98 0.85 0.89 0.91 0.82 1 0.68 0.76 0.98 0.82 0.85; y2=7.67 7.51 7.51 7.7 7.74 7.85 7.67 7.56 7.51 7.55 7.58 7.67.77 7.67 7.65 7.75 7.54 7.46 7.66 7.63 7.8 7.78 7.7 7.68.7.66 7.57 7.51 7.54 7.85 7.57 7.51 7.85 7.61 7.65 7.64 7.68 7.

42、687.69 7.72 7.48 7.66 7.55 7.77 7.71 7.66 7.68 7.6 7.66 7.76 7.57.62 7.71.7.59 7.5 7.77 7.65 7.46 7.52 7.45 7.69 7.687.65 7.57 7.72 7.6 7.61 7.55 7.71 7.6 7.66 7.48 7.56 7.59 7.647.68 7.59 7.6 7.54 7.78 7.56 7.59 7.75 7.72 7.69 7.57 7.74 7.597.747.547.54;y3=30*（0.073818*x.-0.52339*x+0.97907）;y4=50*（

43、0.073818*x.-0.52339*x+0.97907）;plot（x,y 3, ' r* ' ,x,y4, ' sk' ,x11,y1,x12,y1,x2,y2,' m -'）;legend（ 煤尘爆炸浓度最小，煤尘爆炸浓度最大，断电浓 度，报警浓度）;xlabel（空气中甲烷浓度）;ylabel（ 煤尘所含浓度）;title（ 工作面2各指标的分布情况）;hold off运行结果：掘进面：functionheliang07clear ;clc; hold on x=0:0.02:4;y1=0:0.01:50; x11 = 1;%断电浓度

44、x12=1;% 报警浓度 x2=0.260.26 0.29 0.23 0.26 0.22 0.27 0.26 0.16 0.31 0.19 0.18 0.28 0.24 0.21 0.2 0.23 0.22 0.23 0.23 0.21 0.28 0.26.0.19 0.160.17 0.19 0.32 0.23 0.16 0.22 0.24 0.27 0.19 0.17 0.26 0.25 0.190.27 0.21 0.2 0.27 0.15 0.24 0.29 0.29 0.23 0.24 .0.240.27 0.24 0.16 0.25 0.24 0.29 0.21 0.2 0.27

45、0.3 0.24 0.26 0.230.2 0.25 0.24 0.13 0.17 0.23 0.27 0.33 0.2 0.19 0.2 0.22.0.19 0.23 0.18 0.23 0.24 0.18 0.28 0.2 0.14 0.19 0.29 0.21 0.30.2 0.24 0.26; y2=7.44 7.45 7.62 7.31 7.47 7.56 7.42 7.487.18 7.28 7.02 7.71 7.46 7.55 7.34 7.48 7.24 7.17 7.35 7.54 7.557.45 7.21 7.37 7.2 7.47 7.34 7.46 7.39 7.5

46、1 7.26 .7.246.99 7.21 7.33 7.2 7.07 7.74 7.337.33 7.73 7.34 7.59 7.26 7.167.16 7.28 6.98 .7.2 7.08 7.67 7.66 7.12 7.38 7.21 7.497.5 7.21 7.37 7.16 7.35 7.03 7.53 7.63 .7.7 7.48 7.557.58 7.45 7.09 7.14 7.29 7.2 7.44 7.3 7.34 7.68 7.48 7.15 7.24 .7.227.137.617.137.147.087.267.217.337.26;y3=30*（0.073818*x.-0.52339*x+0.97907）;y4=50*（0