全国大学生数学建模竞赛优秀论文Word格式文档下载.doc
- 文档编号:6872038
- 上传时间:2023-05-07
- 格式:DOC
- 页数:37
- 大小:1.46MB
全国大学生数学建模竞赛优秀论文Word格式文档下载.doc
《全国大学生数学建模竞赛优秀论文Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《全国大学生数学建模竞赛优秀论文Word格式文档下载.doc(37页珍藏版)》请在冰点文库上搜索。
煤尘本身具有爆炸性;
煤尘悬浮于空气中并达到一定的浓度;
存在引爆的高温热源。
试验表明,一般情况下煤尘的爆炸浓度是,而当矿井空气中瓦斯浓度增加时,会使煤尘爆炸下限降低,结果如附表1所示。
国家《煤矿安全规程》给出了煤矿预防瓦斯爆炸的措施和操作规程,以及相应的专业标准(见附件2)。
规程要求煤矿必须安装完善的通风系统和瓦斯自动监控系统,所有的采煤工作面、掘进面和回风巷都要安装甲烷传感器,每个传感器都与地面控制中心相连,当井下瓦斯浓度超标时,控制中心将自动切断电源,停止采煤作业,人员撤离采煤现场。
具体内容见附件2的第二章和第三章。
1.2问题提出
附图1是有两个采煤工作面和一个掘进工作面的矿井通风系统示意图,请你结合附表2的监测数据,按照煤矿开采的实际情况研究下列问题:
(1)根据《煤矿安全规程》第一百三十三条的分类标准(见附件2),鉴别该矿是属于“低瓦斯矿井”还是“高瓦斯矿井”。
(2)根据《煤矿安全规程》第一百六十八条的规定,并参照附表1,判断该煤矿不安全的程度(即发生爆炸事故的可能性)有多大?
(3)为了保障安全生产,利用两个可控风门调节各采煤工作面的风量,通过一个局部通风机和风筒实现掘进巷的通风(见下面的注)。
根据附图1所示各井巷风量的分流情况、对各井巷中风速的要求(见《煤矿安全规程》第一百零一条),以及瓦斯和煤尘等因素的影响,确定该煤矿所需要的最佳(总)通风量,以及两个采煤工作面所需要的风量和局部通风机的额定风量(实际中,井巷可能会出现漏风现象)。
二、问题的分析
2.1背景的分析
煤矿安全生产是目前社会重点关注的热点问题之一,尤其是在能源紧张,对煤碳的需求量不断增加的情况下,煤矿的安全生产问题更是值得我们关注,这也是建设平安和谐社会的重要组成部分。
根据统计资料,可知大部分煤矿事故的罪魁祸首都是瓦斯或煤尘爆炸。
因此,矿井下的瓦斯和煤尘对煤矿的安全生产构成了重大威胁,做好井下瓦斯和煤尘的监测与控制是实现煤矿安全生产的关键环节。
2.2基本预备知识
2.2.1《煤矿安全规程》第一百三十三条中,矿井瓦斯等级根据矿井相对瓦斯涌出量和矿井绝对瓦斯涌出量划分为:
(1)低瓦斯矿井:
矿井相对瓦斯涌出量小于或等于,且绝对瓦斯涌出量小于等于;
(2)高瓦斯矿井:
矿井相对瓦斯涌出量大于,或绝对瓦斯涌出量大于;
相对瓦斯量定义:
是指平均日产一吨煤所涌出的瓦斯量,单位为;
绝对瓦斯量定义:
是指矿井单位时间涌出的瓦斯体积,单位为或。
2.2.2根据《煤矿安全规程》第一百六十八条的规定,甲烷传感器报警浓度、断电浓度、复电浓度和断电范围必须符合表3规定(具体表3见附件2)。
2.2.3根据《煤矿安全规程》第一百零一条,井巷中的风流速度应符合表2要求(具体表2见附件2)。
2.3问题的分析
2.3.1问题1的分析
需根据《煤矿安全规程》第一百三十三条的分类标准,鉴别该矿是属于“低瓦斯矿井”还是“高瓦斯矿井”。
由分类标准可知,须考察出该矿的相对瓦斯涌出量和绝对瓦斯涌出量的值,与其分类标准值进行鉴别。
由附表2所给监测值,可根据绝对瓦斯涌出量与相对瓦斯涌出量的计算公式,算出各监测点的绝对瓦斯涌出量与相对瓦斯涌出量。
如果经考察出的监测点的相对瓦斯量有小于或等于且绝对瓦斯量小于等于,则鉴定该煤矿属于低瓦斯矿井。
而如果经考察出的监测点的相对瓦斯量有大于或绝对瓦斯量大于,则鉴定该煤矿属于高瓦斯矿井。
2.3.2问题2的分析
根据《煤矿安全规程》第一百六十八条的规定,并参照附表1,判断煤矿不安全的程度(即发生爆炸事故的可能性)有多大。
可知对煤矿不安全程度评价一般采用间接的方法,通过对影响事故发生可能性和后果严重程度的各内因和外因的分析与综合,可以得到不安全性的相对值。
目前一般的不安全性评价结果,都是相对不安全性,而非绝对不安全性。
即煤矿发生爆炸的可能性为相对的,不是绝对的。
因此,假设只考虑瓦斯爆炸与煤尘爆炸,需定义不同浓度瓦斯与煤尘发生爆炸事故的可能性,并综合瓦斯爆炸的可能性与煤尘爆炸的可能性,即为该煤矿发生爆炸事故的可能性.
2.3.3问题3的分析
满足各井巷中风速的要求及瓦斯和煤尘等因素的影响约束,确定煤矿所需要的最佳总通风量,以及两个采煤工作面所需要的风量和局部通风机的额定风量,这是一个有多约束条件的优化问题。
首先需根据各井巷风量的分流情况,确定总通风量为哪些巷道的风量之和。
再根据《煤矿安全规程》第一百零一条规定,其各巷道的风速就满足规定的风速要求,及满足瓦斯和煤尘浓度的要求的约束。
三、模型的假设
1、各监测站点的工作是相互独立的;
2、附表中的监测值均为有效值,忽略其测量误差,且每天各班次的监测数据为该班次内的平均监测值;
3、煤矿的生产是严格按照国家《煤矿安全规程》进行生产;
4、煤矿爆炸只考虑由瓦斯爆炸和煤尘爆炸,不考虑其他如矿井温度,机器摩擦及一些由人为失误造成的爆炸;
5、煤尘爆炸下限取其中位数,瓦斯爆炸下限取值为。
四、符号约定
:
监测点的风速(单位:
);
巷道横断面面积(单位:
监测点的风量(单位:
矿井的绝对瓦斯涌出量(单位:
风流中的平均瓦斯浓度,即体积百分比(单位:
%);
矿井中的煤尘(单位:
矿井的相对瓦斯涌出量(单位:
矿井的日产量(单位:
煤尘爆炸下限(单位:
瓦斯爆炸下限(单位:
在空气中有瓦斯时,煤尘降低系数;
在空气中有瓦斯时,煤尘发生爆炸的下限(单位:
煤尘爆炸对矿井的不安全性大小;
瓦斯爆炸对矿井的不安全性大小;
煤矿的不安全性大小。
五、模型的建立与求解
5.1问题1的分析与求解
5.1.1绝对瓦斯涌出量与相对瓦斯涌出量的计算公式
由问题的分析,鉴定矿井是属于“低瓦斯矿井”还是“高瓦斯矿井”,需算出该矿的绝对瓦斯量与相对瓦斯涌出量值,与分类标准值进行鉴别。
由绝对瓦斯涌出量与相对瓦斯涌出量的定义,结合相关的符号约定,可知
风量为风速在1分钟传播的距离乘以相应巷道横断面面积,公式为:
……
(1)
绝对瓦斯涌出量计算公式为:
……
(2)
一天24小时,且1小时60分钟,绝对瓦斯量的单位为,所以相对瓦斯涌出量的计算公式为:
……(3)
5.1.2矿井的绝对瓦斯涌出量与相对瓦斯涌出量的计算
由附表2所给各监测点的风速、瓦斯和日产量数据,及各监测点巷道的横断面面积,代入计算公式可得出矿井在各监测点各班次的绝对瓦斯量及相对瓦斯量。
对各监测点进行编号,有个监测点,,依次表示工作面,工作面,掘进工作面,回风巷,回风巷;
对矿井日生产班次进行编号,有个班次,,依次表示为早班,中班,晚班;
共有30天的监测数据,记天数为,。
则对应的有第天第个监测点第个班次的风速、瓦斯的监测值和日产量监测值,分别记为:
,,;
第天第个监测点巷道记为。
日产量取为30天的月平均日产量:
根据公式
(1),则可得第天第个监测点第个班次的风速为:
……(4)
根据公式
(2),第天第个监测点第个班次的绝对瓦斯涌出量为:
……(5)
第天第个监测点的平均绝对瓦斯涌出量为:
……(6)
根据公式(3),第天第个监测点第个班次的相对瓦斯涌出量为:
……(7)
第天第个监测点的平均相对瓦斯涌出量为:
……(8)
第个监测点30天的平均绝对瓦斯涌出量为:
……(9)
第个监测点30天的平均相对瓦斯涌出量为:
……(10)
则代入附表2所给的相应数据,可得各监测点30天的平均绝对瓦斯量和相均绝对瓦斯量。
5.1.3煤矿类型的鉴别
根据附图1(煤矿的通风系统示意图),认为井巷中出现漏风的可能性较小,并由数据可得,工作面、工作面、掘进工作面的瓦斯涌出量之和与回风巷、回风巷的瓦斯涌出量之和与总回风巷的瓦斯涌出量可认为是等价的。
这里取总回风巷的相对瓦斯涌出量和绝对瓦斯涌出量的值,依据矿井的分类标准来鉴别该矿井。
经MATLAB软件编程得出,30天回风巷的相对瓦斯涌出量如表1所示,均大于,且30天的平均相对瓦斯涌出量为,大于。
绝对瓦斯涌出量为,小于。
所以,根据矿井的分类标准,该矿井属于高瓦斯矿井。
表1回风巷30天的相对瓦斯涌出量(N:
天数;
A:
相对瓦斯涌出量,单位)
N
1
2
3
4
5
6
7
8
9
10
A
23.8
23.37
24.08
22.83
22.73
24.26
23.08
22.55
22.56
23.14
11
12
13
14
15
16
17
18
19
20
22.74
21.44
24.24
23.36
22.24
24.77
23.28
22.62
23.12
21
22
23
24
25
26
27
28
29
30
23.25
23.15
24.65
22.82
22.97
22.80
22.46
23.18
23.56
5.2问题2的分析与求解
5.2.1煤矿不安全程度(即发生爆炸事故的可能性)的定义
煤矿发生爆炸是随机、不确定的,所以煤矿发生爆炸的可能性是相对的,不是绝对的。
由问题的分析可知《煤矿安全规程》第一百六十八条规定,并参照附表1,判断煤矿不安全的程度有多大?
在此假设煤矿爆炸只考虑由瓦斯浓度引起的爆炸和煤尘浓度引起的爆炸,不考虑其他如矿井温度,机器摩擦及一些由人为失误造成的爆炸。
引用煤尘与在瓦斯浓度影响下煤尘的爆炸下限的偏离程度来恒量,由煤尘引起爆炸的可能性;
引用瓦斯浓度与瓦斯爆炸下限的偏离程度来恒量瓦斯爆炸的可能性。
若偏离值越大,煤矿的安全性越好;
若偏离值越小,煤矿的安全性越差。
在此采用了一个最大型心理函数计算其不安全的可能性。
通过给瓦斯不安全程度函数与煤尘不安全程度函数赋予不同权系数,平衡两者的不安全程度,综合两种发生爆炸的可能性即为该煤矿发生爆炸的可能性。
综合时,进行了对瓦斯的不安性与煤尘的不安全性赋权处理。
5.2.2煤矿安全性的计算
第监测点第班次的瓦斯浓度、煤尘、在有瓦斯时煤矿降低系数及相应的煤尘发生爆炸的下限分别为,,,。
煤尘爆炸下限一般为,取其中位值;
且瓦斯爆炸下限,取值为。
采用求解其偏离值的大小,即其不安全性的大小。
则第监测点第班次煤尘对矿井的不安全性为:
……(11)
其中,;
用MATLAB软件编程中的线性最小二乘法进行拟合可以将不同瓦斯浓度对应的煤尘降低系数算出。
共对30天进行监测,且监测出每天3个班次的数据,则对应的各监测点共有90个班次监测值。
第监测点平均每班次煤尘爆炸对矿井的不安全性大小为:
……(12)
煤尘爆炸对矿井的不安全性大小先对6个监测点赋权处理综合成一个点,为:
……(13)
……(14)
第监测点平均每班次瓦斯爆炸对矿井的不安全性大小为:
……(15)
瓦斯爆炸对矿井的不安全性大小同样先对6个监测点赋权处理综合成一个点,为:
……(16)
煤尘爆炸与瓦斯爆炸可以看为两个相互独立的事件,只要煤尘爆炸或瓦斯爆炸这两个事件任意有一个发生,则整个煤矿是不安全的。
根据概率统计知识中任意事件概率的加法公式,得煤矿的不安全性的大小为:
……(17)
经MATLAB软件编程得出取其对六个面影响最大的作为整个矿井的不安全程度,即煤矿发生爆炸的不安全性(煤矿发生爆炸事故的可能性)大小为。
在此根据以上模型把给出的煤尘浓度与瓦斯浓度对应的矿井不安全可能性的大小如列表2如下。
根据表2所得结果,人们可以从中看出煤尘浓度与瓦斯浓度对应的矿井不安全性的大小。
5.3问题3的分析与求解
5.3.1总通风量的定义及公式
根据附图1(煤矿的通风系统示意图)中各巷道的分布位置及各处风的流向(即分流情况),可把总通风量分为三大块,进风巷、进风巷及局部通风机所在巷(包括局
部通风机的风量)的风量(分别记为,,)。
表2:
矿井不安全性的大小(表示煤尘浓度;
表示瓦斯浓度)
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
6.0
0.01
0.02
0.03
0.04
0.05
0.06
0.08
0.10
0.12
0.15
6.5
0.07
0.09
0.11
0.14
0.16
7.0
0.13
0.18
7.5
0.17
0.20
8.0
0.19
0.22
8.5
0.21
0.25
9.0
0.123
0.24
0.28
9.5
0.90
0.32
0.23
0.27
0.31
0.36
局部通风机所在的巷道中至少需要有的余裕风量(新鲜风)才能保证风在巷道中的正常流动,否则可能会出现负压导致乏风逆流,即局部通风机将乏风吸入并送至掘进工作面。
记余裕通风量为,局部通风量为。
所以求最小总通风量的目标函数为:
……(18)
其中;
;
。
5.3.2风速的约束
由《煤矿安全规程》第一百零一条的规定,得各巷道的风速范围约束。
记风速为,,分别表示进风巷、进风巷、采煤工作面、回风巷、采煤工作面、回风巷、总回风巷、掘进工作面的风速。
相应的对于采煤工作面,采煤工作面,掘进工作面,其区域内有绝对瓦斯涌出量,则其处的风量应为进风巷的风量加绝对瓦斯涌出量,等于回风巷的风量。
所以,,,。
总回风巷的风量为进风巷、进风巷的风量之和与采煤工作面、采煤工作面、局部通风机所在巷的绝对瓦斯涌出量之和。
则总回风巷的风速值为:
由《煤矿安全规程》第一百零一条的规定,各巷道的风速范围约束为:
(),,,(),。
还需考虑各巷道中瓦斯和煤尘等因素的影响,首先通过附表2所给的数据用MATLAB软件编程,采用线性最小二乘法把风速与瓦斯及风速与煤尘的函数关系式插值拟合出来,分别记为,。
通过编程,得出风速对应各巷道瓦斯的函数关系式。
记,,,,,分别表示风速采煤工作面,回风巷,采煤工作面,回风巷,总回风巷,掘进工作面对应瓦斯的关系式,其为:
由上关系式,得对应风速,则有各巷道的瓦斯浓度,记为,,,,,。
由所给监测数据,取各监测点瓦斯浓度的最大值,作为最优瓦斯浓度的约束。
即
,,,,,。
相应的煤尘也有一安全浓度,即煤尘应小于其对应在有瓦斯时煤尘的爆炸下限。
对应值通过附表1,同样采用线性最小二乘法把瓦斯浓度与的函数关系插值拟合出来。
经MATLAB软件编程,得
煤尘爆炸下限仍取中位数,则
由不同巷道的瓦斯浓度,有其对应的煤尘浓度为:
则:
5.3.3最佳总风量的模型
综上所得,求得其最佳总风量的模型如下:
s.t.
()
(),
,()
其中,,,,,,,,,,,
5.3.4模型求解
经LINGO软件编程求解,求解最佳总通风量为,采煤工作面的风量为,采煤工作面的风量为,局部通风机的额定风量为。
六、误差分析
误差来源:
1、各监测站点在实际监测中,有观测误差,即存在监测数据与实际数据的误差。
2、在模型的建立中,有模型误差。
即对监测数据的处理大部分对其取平均值,由模型所得的解与实际问题的解之间存在一定的误差。
3、煤矿发生爆炸,在此只考虑瓦斯爆炸和煤尘爆炸,由此所得的煤矿发生爆炸的可能性与实际煤矿生产中有一定的误差。
4、模型的误差:
通过过建立的模型求出的解与实际的值间存在一定的误差,如取不安全的可能性大小,只是取其相对值。
而现实中是否发生不安全事故是随机的,不确定的。
5、舍入误差:
在计算时取的是小数点后两位,其数据有一定的误差。
七、模型的改进
问题2的改进:
7.2.1煤矿安全性的分析
7.2.2煤矿安全性的计算
采用最小二乘法求解其偏离值的大小,即其不安全性的大小。
用MATLAB软件编程可以将不同瓦斯浓度对应的煤尘降低系数可由三次样条插值法进行插值算出。
共对30天进行监测,且监测出每天3个班次的数据,则对应的各监测点共有90个监测值。
煤尘爆炸对矿井的不安全性大小取6个监测点的平均值,为:
……(14)
瓦斯爆炸对矿井的不安全性大小取6个监测点的平均值,为:
……(15)
……(16)
经MATLAB软件编程得出,即煤矿发生爆炸的不安全性(煤矿发生爆炸事故的可能性)为。
问题3的改进:
为了计算,我们进一步简化公式,即风速的约束作进一步简化,作为问题3的改进。
7.3.1总通风量的定义及公式
根据附图1(煤矿的通风系统示意图)中各巷道的分布位置及各处风的流向(即分流情况),可把总通风量分为三大块,进风巷、进风巷及局部通风机所在巷(包括局部通风机的风量)的风量(分别记为,,),和各巷道漏的风量。
由《煤矿安全规程》第一百一十条,各巷道的漏风率不超过,取漏风率最大值,则对应的各巷道通风量需加上对应的通风量的。
并知局部通风机所在的巷道中至少需要有的余裕风量(新鲜风)才能保证风在巷道中的正常流动,否则可能会出现负压导致乏风逆流,即局部通风机将乏风吸入并送至掘进工作面。
所以根据掘进巷道图,局部通风机所在巷的通风量为局部通风机的通风量除以(),记局部通风机的通风量为。
……(17)
,。
7.3.2总通风量的约束
……(18)
记,,分别表示风速对应进风巷、进风巷及局部通风机所在巷瓦斯的关系式,其为:
由上关系式,得对应风速,则有每一巷道的瓦斯浓度,记为,,。
由《煤矿安全规程》第一百六十八条的规定,其瓦斯浓度有一安全浓度约束,得
,,
对应值通过附表1,同样采用线性最小二乘法把瓦斯
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 全国大学生 数学 建模 竞赛 优秀论文