煤矿瓦斯和煤尘的监测与控制模型.doc

煤矿瓦斯和煤尘的监测与控制模型 【摘要】如何划分煤矿瓦斯等级、判断煤矿不安全程度和确定煤矿所需要的最佳通风量，对煤矿安全生产具有非常重要的意义。本文通过建立三种不同的模型并通过MATLAB、EXCLE计算分别对矿井分类问题、矿井不安全程度问题和矿井所需最佳通风量问题进行了合理解决。在问题1中，利用EXCLE作数据分析，根据煤矿安全规程第一百三十三条的分类标准鉴别出该矿是属于 “高瓦斯矿井”；在问题2中，我们分三种情况讨论，根据试验表明，一般情况下煤尘的爆炸浓度是得出煤尘单独影响的不安全度为0，再利用matlab对数据进行拟合，求出该矿井瓦斯和煤尘的综合影响的不安全度为0，我们又根据煤矿安全规程第一百六十八条的规定用excle对数据进行处理得到该矿井的不安全程度为4.63% 。在问题3中，首先根据瓦斯浓度、煤尘浓度与风速的关系，找出安全生产前提下采煤工作面实际需要的最小风量。通过简单计算得到，采煤工作面一和采煤工作面二实际所需最小风量分别为424.908m3/min和554.09m3/min。根据各井巷风量的分流情况以及对各井巷中风速的具体要求，建立了煤矿实际所需最佳通风量的单目标最优化模型。得出该矿井局部通风机的额定风量和煤矿所需要的最佳(总)风量分别为150m3/min和1566.56m3/min。最后，通过模型求解和分析知，本文对不同问题的建模方法合理，对煤矿公司等高风险的企业的管理具有较强的借鉴作用。关键词：matlab 瓦斯含量 煤尘浓度 爆炸极限一、问题重述煤矿安全生产是我国目前亟待解决的问题之一，做好井下瓦斯和煤尘的监测与控制是实现安全生产的关键环节。瓦斯是一种无毒、无色、无味的可燃气体，其主要成分是甲烷，在矿井中它通常从煤岩裂缝中涌出。瓦斯爆炸需要三个条件：空气中瓦斯达到一定的浓度；足够的氧气；一定温度的引火源。煤尘是在煤炭开采过程中产生的可燃性粉尘。煤尘爆炸必须具备三个条件:煤尘本身具有爆炸性；煤尘悬浮于空气中并达到一定的浓度；存在引爆的高温热源。试验表明，一般情况下煤尘的爆炸浓度是30 2000g/m3，而当矿井空气中瓦斯浓度增加时，会使煤尘爆炸下限降低。国家煤矿安全规程给出了煤矿预防瓦斯爆炸的措施和操作规程，以及相应的专业标准。规程要求煤矿必须安装完善的通风系统和瓦斯自动监控系统，所有的采煤工作面、掘进面和回风巷都要安装甲烷传感器，每个传感器都与地面控制中心相连，当井下瓦斯浓度超标时，控制中心将自动切断电源，停止采煤作业，人员撤离采煤现场。现按照煤矿开采的实际情况考虑下列问题： 问题一：根据煤矿安全规程第一百三十三条的分类标准(见附件2)，鉴别该矿是属于“低瓦斯矿井”还是“高瓦斯矿井”。问题二：根据煤矿安全规程规定，并参照附表1，判断该煤矿不安全的程度(即发生爆炸事故的可能性)有多大？ 问题三：如果为了保障安全生产，利用两个可控风门调节各采煤工作面的风量，通过一个局部通风机和风筒实现掘进巷的通风。根据附图1所示各井巷风量的分流情况、对各井巷中风速的要求(见煤矿安全规程第一百零一条)，以及瓦斯和煤尘等因素的影响，确定该煤矿所需要的最佳（总）通风量，以及两个采煤工作面所需要的风量和局部通风机的额定风量(实际中，井巷可能会出现漏风现象)。二、问题的背景煤矿安全生产是目前社会重点关注的热点问题之一，尤其是在能源紧张，对煤碳的需求量不断增加的情况下，煤矿的安全生产问题更是值得我们关注，这也是建设平安和谐社会的重要组成部分。据不完全统计，2005年全国煤矿发生事故405起，死亡2157人。瓦斯事故死亡人数占全国煤矿事故总死亡人数的36.0%。在58起一次死亡10人以上的特大事故中，瓦斯事故40起，占69 %。一次死亡百人以上的事故5起。这些事故所造成的经济损失是重大的，给社会和伤亡人员的家庭所造成影响与损失是无法估量的。我们注意到，大部分事故的罪魁祸首都是瓦斯或煤尘爆炸，瓦斯在煤矿的开采中是不可避免的。因此，矿井下的瓦斯和煤尘对煤矿的安全生产构成了重大威胁，做好煤矿井下瓦斯和煤尘的监测与控制是保证煤矿安全生产的关键所在。3、 模型假设、符号说明、名词解释3.1模型假设1.采煤工作面1与采煤工作面2的每天的采煤量相同，早、中、晚班的采煤速率相同；2.假设附表2中的监测数据能准确反映监测点的实际情况，数据是真实可靠的；3.煤层的瓦斯涌出速率和煤尘的生成速率基本稳定；4假定工作区域、掘井区域、回风巷区域以及总回风巷区域的瓦斯浓度、煤尘浓度是均匀的；5.假设早班、中班、晚班的工作时间相同，均为八小时；6.假设漏风很小，但为了安全考虑我们设安全系数为1.2；7.假设局部通风机不会在此过程中发生机械故障；3.2符号说明见（表（1）表（1）符号的说明符号符号定义单位风速风流中瓦斯平均浓度绝对瓦斯涌出量相对瓦斯涌出量日产量工作面横截面积该煤矿需要的通风量局部通风机的额定风量主巷道断面面积其它各采煤区的进风巷和回风巷以及掘进巷的断面面积 局部通风机的风筒的断面面积 3.3名词解释1采煤工作面：矿井中进行开采的煤壁 (采煤现场)。2掘进巷：用爆破或机械等方法开凿出的地下巷道，用以准备新的采煤区和采煤工作面。3掘进工作面：掘进巷尽头的开掘现场。4新鲜风：不含瓦斯和煤尘等有害物质的风。5乏风：含有一定浓度的瓦斯和煤尘等有害物质的风流。四、模型建立4.1问题（1）一、绝对瓦斯涌出量：单位时间涌出的瓦斯体积。（1-1）二、 相对瓦斯涌出量：平均日产一吨煤同期所涌出的瓦斯量。 （1-2）根据式子（1-1）与（2-2）我们得到模型一如下；（1-3）我们用excle对数据进行运算和筛选（见附录1）根据煤矿安全规程第一百三十三条由于低瓦斯矿井必须要同时满足两个条件：矿井相对瓦斯涌出量小于或等于10；矿井绝对瓦斯涌出量小于或等于40。而高瓦斯矿井只须满足（矿井相对瓦斯涌出量大于10；矿井绝对瓦斯涌出量大于40）两者中的一个即可。我们根据附录（1）中的数据不难判断该矿井为高瓦斯矿井。4.2问题（2）我们从三个方面考虑该矿井的不安全度：4.2.1煤尘产生的不安全度:试验表明，一般煤尘的爆炸浓度是我们用excle对煤尘数据进行筛选得出煤尘的最大值为，该煤矿是不会发生煤尘爆炸。4.2.2瓦斯和煤尘的综合因素的不安全度：我们用matlab对表（2）中的数据进行拟合，出于安全考虑我们考虑煤尘爆炸下限。表（2）瓦斯浓度和煤尘爆炸下限浓度关系空气中瓦斯浓度(%)00.51.01.52.02.53.03.54煤尘爆炸下限浓度(g/m3)305022.537.5152510.517.56.512.54.57.5352.53.51.52.5我们通过点的分布考虑它为对数函数，用matlab处理（见附录（2）得到的函数如下： （21）我们建立 （2-2）当表中的煤尘浓度超过拟合函数计算出来的煤尘浓度我们就记一次数，依此累加，我们用matlab编程运算得出煤尘浓度超标的次数为0（见附录（3），我们可以判断瓦斯和煤尘的综合因素的不安全度为0。4.2.3瓦斯产生的不安全度： 我们根据煤矿安全规程第一百六十八条的规定，各个考察点在任意天任意工作时段报警、断电、复电等事故的发生与否，我们用excle对数据进行筛选得到的结果见表（3）表（3） 各个工作面的报警、断电、复电次数表 监测点次数工作面工作面掘进工作面回风巷回风巷报 警050020断 电000020复 电000012通过对各个工作面的考察，可以得出在一个月中所有统计的次数中有25次报警，而这25次报警都有可能造成安全事故，则我们可以运用简单的统计概率计算公式综合上面三种因素我们得出该煤矿的不安全度为0.04634.3问题（3）由煤矿安全规程第一百零三条，矿井所需风量应符合本规定对风流中有关瓦斯浓度、煤尘浓度的要求。而由煤矿开采的实际情况，不难得知在满足各种约束条件的前提下，该煤矿所需要的最佳（总）通风量，为两个采煤工作面、一个掘进工作面和主巷道所需最小风量之和。因此我们首先求出两个采煤工作面所需要的实际风量。 ( 1 ) 根据煤矿安全规程规定，按采煤工作面回风巷中瓦斯浓度不超过1%的要求，并按照瓦斯绝对涌出量计算采煤工作面所需要的最小风量为： （3-1）式中：表示采煤工作面一、采煤工作面二实际需要风量；表示主巷道的风量；表示采煤工作面回风巷风流中日平均瓦斯绝对涌出量；表示采煤工作面瓦斯不均衡系数。100采煤工作面回风流中瓦斯浓度不超过1%所换算的常数。(2 )对于采煤工作面一工作面回风巷风流中日平均瓦斯绝对涌出量可以利用excel电子表格算出正常生产条件下，连续观测1 个月月平均日瓦斯绝对涌出量=3.7945。最大瓦斯绝对涌出量为4.24944对于采煤工作面二工作面回风巷风流中日平均瓦斯绝对涌出量可以利用excel电子表格算出连续观测1 个月月平均日瓦斯绝对涌出量=4.4021。最大瓦斯绝对涌出来为5.54112。 （3-2）计算得采煤工作面一的=1.1198，计算其风量=采煤工作面二的=1.2587，计算其分量=(3)考虑风量对煤尘浓度的影响由于采煤工作面回风巷中瓦斯浓度不超过1%的规定，我们可以根据附表1（瓦斯浓度与煤尘爆炸下限浓度的关系表）可知，当空气中瓦斯浓度不超过1%时，煤尘爆炸下限浓度的最小值为15，然而我们从附表2（各监测点的风速，瓦斯和煤尘的监测数据表）知在采煤工作面的煤尘浓度最大值为8.42。远小于瓦斯浓度不超标时的煤尘爆炸下限浓度。所以，在正常生产的情况下，风量对煤尘浓度影响不大，即是只要瓦斯浓度不超标时，煤尘浓度绝对不会达到爆炸下限。故我们不考虑风量对煤尘浓度的影响。为了保证安全生产和符合经济要求，确定该煤矿所需要的最佳（总）通风量，除了采煤工作面一和采煤工作面二所需风量外，还要考虑掘进工作面和主巷道所需的风量。因此我们根据煤矿安全规程第一百零一条建立了相关的优化模型四： （4）说明：根据实际瓦斯绝对涌出量可以得到采煤工作面一和采煤工作面二的所需风量，我们可以知道进入采煤工作面所需风量不得低于采煤工作面的所需风量。对于约束条件5是由于至少要有15%的余裕风量在巷道中，否则可能会出现负压乏风逆流。按照瓦斯绝对涌出量可以计算采煤工作面所需要的风量，其结果为采煤工作面一和采煤工作面二所需风量分别424.908、554.09主巷道的风量为。模型三是一个多约束条件的最优化问题模型，解得到最佳通风量为。局部通风机的额定风量为150。七 模型的结果分析7.1模型一的结果分析根据煤矿瓦斯等级鉴定方法和计算公式，我们运用了(EXCEL)处理数据和绝对瓦斯涌出量以及相对瓦斯涌出量的计算方法求解得出了工作面一、工作面二以及掘进工作面的绝对瓦斯涌出量和相对瓦斯涌出量的值，可以发现采煤工作面二的相对瓦斯涌出量超出了分类标准(大于10),所以,我们可以根据煤矿安全规程第一百三十三条的分类标准得知该矿井为高瓦斯矿井。7.2模型二的结果分析由于煤尘爆炸的可能性为0。于是我们只需根据所求出的报警次数与可能发生爆炸点总和之比来计算出爆炸的可能性。其可能性即是该煤矿的不安全的程度。7.3模型三的结果分析我们利用不均衡系数和假设（6）计算最佳通风量，不仅可以保证矿井内瓦斯浓度小于它的上限值1%，煤尘含量低于它的爆炸下限，既保证了安全生产，而且还合乎经济要求。因此，得出我们计算出的结果是合理的。八 模型的评价与推广8.1模型的评价我们从各个问题的出发点，分析了应该考虑的各种情况，建立了三个不同模型，并进行了理论论证和实践论证。从而证明，我们建立的数学模型能较好的解决问题。此模型的最大好处在于模型极为广泛的使用性。 (一) 模型的优点1. 建立的模型能与实际紧密联系，结合实际情况对所提出的问题进行求解，使模型更贴近实际，通用性、推广性较强。2. 模型一中，我们利用Excel软件对数据进行了验证，从而确保了模型建立的合理性。3.在模型二中，我们利用了模糊统计方法，统计出了所有监测点的可能发生爆炸次数，具有更为准确的估计了煤矿的不安全程度4. 通过利用数学工具和matlab编程的方法，严格的对模型求解，具有科学性。 5. 本文对模型的合理性以及算法进行了讨论，为模型的推广和解决同类型问题提供有价值的参考。(二) 模型的缺点1.在模型建立中，我们只考虑了影响矿井安全生产的几个主要因素，忽略了其他因素的影响，使得结果可能存在一定的误差。2.在处理巷道内的风量问题时，没有考虑到风通过巷道的阻力问题，还有工人工作是所需要的呼吸风量，所以此模型有待改进，求出与时间更符合的总需风量。8.2模型的推广本文所建立的三个模型可以分别去判定一个煤矿是低瓦斯煤矿还是高瓦斯煤矿、安全程度以及最经济和安全条件下的最佳通风量问题。我们建立模型的方法和思想对其他类似的煤矿问题也很实用，基本囊括了几种优秀的数据挖掘方法。特别是本问题中的煤矿安全问题是我国目前亟待解决的问题之一，所以我们做的安全评价和最经济生产条件下的各项监测与控制可以推广到煤矿企业中。进一步建立安全指标体系并应用改进式层次分析法对各指标权重处理,运用模糊评判模型,对煤矿企业进行安全评价。据此再给出的系统整体的安全状况评价结果而对安全管理具有很大的意义。参考文献1 姜启源，谢金星，叶俊.数学模型(第三版) 高等教育出版社 20032 苏金明、王永利，MATLAB实用指南，北京：电子工业出版社，2004年3 萧树铁，数学实验，浙江：高等教育出版社，1999年 附录（1）工作面绝对瓦斯浓度相对瓦斯浓度4.6699211.264134.243210.234854.243210.149854.7030411.24984.485610.729674.485610.108394.6147210.399374.7150410.625444.8153611.256694.384810.250184.2537610.041664.9420811.666554.888811.972574.9627212.15365.08812.460414.2134410.425014.6699211.554444.6699211.115184.9934412.291544.586411.28964.4265610.281044.345210.092084.4140810.45444.8052811.306544.4642410.504095.5411213.167024.263610.131334.263610.370924.4642410.555844.939211.756115.1170412.17944.7601611.329974.586411.367334.2537610.542884.6648811.561844.7318411.061444.7827211.253464.2931210.101464.2134410.095434.8614411.648044.586410.989044.6588810.890894.737611.074914.7095211.533524.9586412.143614.838411.849144.6147210.632314.6648810.747884.2537610.243174.737611.408275.0133612.072314.384810.217014.9862411.637254.5283210.778154.5645610.864414.96812.208054.939211.47169总回风巷绝对瓦斯涌出量相对瓦斯涌出量10.351524.968449.64823.2715610.035924.20729.25222.1310310.348824.75469.895223.6695810.587623.8593810.25723.114379.544521.5087310.546824.654868.943620.907129.326421.801979.534322.50729.343822.05759.808523.1544910.058424.6328210.220425.0295510.342825.329319.34223.114239.634823.8386810.156825.130239.451822.496859.318622.179819.695723.077379.461123.288869.979224.564189.030622.2291710.10123.460399.727822.59369.37821.781169.48622.175069.461122.116869.753622.800629.10821.571588.926221.1419.028821.38410.526424.7689.969623.4578810.