晋普山煤矿通风系统阻力测定与分析.pdf

晋普山煤矿通风系统阻力测定与分析 * 程根银1， 2教授朱锴1教授石建军1李楠1宋晓燕1 （1 华北科技学院安全工程学院2 北京科技大学土木与环境学院） 学科分类与代码 620.4030中图分类号 X924.3 文献标识码 A 【摘要】矿井通风系统阻力测定是煤矿通风技术管理的重要内容， 阻力的大小直接影响矿井的通风效 果， 笔者介绍了矿井通风阻力测定常用的几种方法， 分别评述了其优缺点及适用条件， 并根据晋普山煤矿通风系统 实际情况， 选用基点法对其通风阻力进行了测定， 全面、 准确地掌握了全矿井通风阻力分布情况， 为今后进行通风 系统调整， 提供了有参考价值的基础数据资料。同时依据测定结果， 对晋普山煤矿通风系统进行了合理分析。 【关键词】通风系统； 阻力； 基点法； 等积孔； 风量； 回风 Measurement and Analysis of Resistance in Ventilation System of Jinpushan Mine CHENG Gen-yin1， 2， Prof.ZHU Kai2， Prof.SHI Jian-jun2LI Nan2SONG Xiao-yan2 （1 Safety Engineering College， North China Institute of Science Technology 2 Civil and Environment Engineering School， Beijing University of Science Technology） Classification and code of disciplines620.4030 AbstractThe resistance measurement of ventilation system in mine is an important content in ventilation management. The magnitude of ventilation resistance directly affects its efficacy. Some s commonly used in ventilation resistance measurement are introduced， their advantages and disadvantage as well as their feasibilities are discussed respectively. Basic point is se- lected to measure the ventilation resistance according to actual situation of Jinpushan mine. Through measurement the resistance distribution in whole ventilation system is known providing valuable data from which rational analysis could be made. Key wordsventilation system；resistance；basic point ；equivalent orifice；air volume；air return 1引言 晋普山煤矿位于山西省晋城市郊区， 属国家司法部下属 的地方国有劳改煤矿， 年产量 90 万吨。目前， 开采的煤层为 3 煤层、 9 煤层， 煤层倾角较小， 属于缓倾斜煤层。该矿的 主、 副两平硐为主进风井， 另排矸立井也进一定风量 （约 400 m3/min） ； 1 号风井的 1 号风机型号是 70B2 - 21 - 28， 2 号 风机型号是 ZK2 - GW28， 3 号风井的 1 号风机和 3 号风机型 号均为 ZKX0 - 4W28 ［1］ 。 矿井通风阻力是伴随着矿井通风而产生的， 其对矿井通 风、 矿井安全及经济效益有着重要影响 ［23］ 。因此， 在矿井 设计及实际生产中， 要合理选择巷道参数， 采取有效措施， 降 低通风阻力， 提高经济效益。矿井通风阻力测定是煤矿通风 技术管理工作的重要内容之一 ［4］ 。其目的有二 其一， 了解矿井通风系统中通风阻力大小和分布情况， 为改善矿井通风状况、 降低阻力提供依据； 其二， 通过实测数据测算矿井各类巷道的井巷摩擦阻力 系数和风阻值， 为矿井的通风设计、 系统改造、 调节风压法控 制火灾提供可靠的基础资料。 煤矿安全规程 第 119 条规定“新井投产前应必须进 行一次矿井通风阻力测定， 以后每 3 年至少进行一次。矿井 转入新水平生产或改变一翼通风系统后， 必须重新进行矿井 通风阻力的测定” 。为准确掌握晋普山煤矿井下风网通风阻 力分布情况， 评价未来安全生产布局， 提高矿井防灾和救灾 的科学决策能力， 该课题组于 2004 年 7 月下旬对晋普山煤 矿通风系统进行了通风阻力测定。由于该次测定线路长， 测 量范围大， 要求参与测定的人员多， 所以选用省时省力， 数据 处理工作简单的基点法进行测定， 并对测定误差进行控制， 最后通过测算结果对晋普山煤矿通风系统存在的问题进行 分析并提出了有效的整改方案。 2通风阻力测定方法选择 在矿井通风阻力测定过程中， 因多种干扰因素的影响， 使测定值存在一定误差， 如果测定方法选择得当， 可以在满 第 1 Y 卷 第 9 期 2 0 0 Y 年 9 月 中国安全科学学报 ChinaSafetyScienceJournal Vol . 1Y No . 9 ══ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ Sep . 2 0 0 Y *文章编号 1003 - 3033 （200Y） 0001 - 04 - 27 - 01 足精度的情况下提高测定速度， 并能给测定工作带来方便。 常用的几种测定方法如下 比值校正法、 同步法、 基点法、 逐 点测量法。 2.1比值校正法 在地面井口附近设监测点， 测定地面气压在一定时间 Δt 内的变化量Δp0， 由此对井下测点 i 的静压在同一时间 内的变化值进行校正。此法需对风流流动的热力过程进行 校正， 计算误差较大。适用于浅井测定。 2.2同步法 用两台等精度和型号相同的气压计在风道的起点和末 点上约定时间同时读数， 使一台的读数作为另一台的修正。 该方法的测定精度较高， 适用于井下局部区域阻力测定。 2.3基点法 用两台气压计， 一台井上， 一台井下。井上的一台固定 在入风井口的基点上， 监测地面大气压的变化， 井下的一台 沿测定路线巡回测定。这种方法的优点是测定方便， 省时省 力， 数据处理工作简单， 缺点是误差较大。其测定结果能够 满足一般性要求。 2.4逐点测量法 用 3 台气压计， 一台固定在入风井口的基点上， 作为大 气压力变化校正用， 将另外两台气压计携至井下测点， 沿着 预定的测定路线逐点测定。这种方法的优点是测定精度高， 适用于测定时间长、 范围大的矿井。 3晋普山煤矿通风阻力测定 由于该项目意在掌握全矿井的通风阻力分布情况， 须对 全矿井各个采区的通风阻力进行测定， 且部分区域由于积水 不能通过， 故采用气压计法。同时， 由于测定线路长， 测量范 围大， 要求参与测定的人员多， 所以选用基点法进行测定。 本次测定程序如图 1 所示。 图 1 矿井通风系统阻力测定程序 3.1矿井通风阻力测定路线的选择 矿井通风阻力测定系统如图 2 所示， 共选择了 105 采 区、 903 采区和 901 采区 3 条主测路线， 另又选择了 5 条辅测 路线。选择的路线包含了采区和工作面等主要用风地点， 满 足全矿通风阻力测定的要求。 图 2 矿井通风系统阻力测定测点布置系统图 86 中国安全科学学报 ChinaSafetyScienceJournal 第15卷 2005 ══ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ 年 1）105 采区副平硐口副平硐105 北大巷105 皮 带巷10514 主巷10514 工作面10514 副巷105 回风 北回风1 号风井总回1 号风井。 2）903 采区 副平硐口副平硐西大巷903 北大巷 联络巷903 皮带巷90303 主巷90303 工作面90303 副巷90303 回风903 回风902 回风3 号风井总回 3 号风井。 3）901 采区 副平硐口副平硐南 740 大巷901 轨 道巷106 回风3 号风井总回3 号风井。 3.2主要计算参数 3.2.1空气密度计算 由于井下各测点的相对湿度均在 60以上， 且密度的 精确度对阻力测定的精确度影响不是很大， 故按简化的近 似公式进行计算 （地面的密度仍然按精确计算公式计算） 。 ρ 0.003 458 P / T 式中， ρ 测点的空气密度， kg/m 3； P 测点的绝对压力， Pa； T 测点空气的绝对温度， K。 T 273.15 td ， t d为测点空气的干球温度， ℃。 3.2.2测段风量计算 在同一测段中， 若风路中间无分岔， 测段风量按两个测 点的平均风量计算， 即 Q12 1 2 （Q 1 Q2） 若风路中间有分岔， 则按测段中绝大部分井巷中的风量 计算； 这时 Q12 Q1 或 Q12 Q2； 若中间分岔较多， 则只进行 阻力计算而不进行风阻、 百米风阻及摩擦阻力系数计算。 3.2.3测段平均断面积及周长的计算 本矿井井下大部分巷道断面按三心拱形计算， 部分巷 道按半圆拱形和梯形计算， 立井按圆形断面计算。 一般测段平均断面积是按两端点的断面积用算术平均 法计算。 若两端点断面积相差较大， 则按测段中绝大部分井 巷中的断面积计算， 即取其中的某一端点的断面积作为测 段平均断面积， 特殊情况下， 在中间专门实测。 测段周长是 根据测段断面形状， 由下式计算 U ﾍC S， 其中 U为测段平 均周长， S 为测段平均断面积， C 为断面类型系数， 对于三心 拱，C 3.85；对 于 半 圆 拱 形，C 3.90；对 于 梯 形， C 4.16。 3.2.4通风阻力计算 hr1-2 K1 （p 1- p2）- K2 （P 01- P02） 1 2 （ρ1v2 1- ρ2v 2 2）ρ1-2g （Z1- Z2） 式中， hr1-2 测点 1 与 2 间的通风阻力， Pa； P1 ， P 2 井下测量气压计在巷道内沿风流方向的 始、 末测点压力读值， Pa； P01 ， P 02 读取 P1及 P2压力值时， 基点校正气压计 的相应读值， Pa； K1 ， K 2 井下测量用气压计与基点校正气压计的 仪器校正系数； v1 ， v 2 始、 末测点处的平均风速， m/s； ρ1，ρ2 始、 末测点的空气密度， kg/m3； Z1 ， Z 2 始末测点的标高， m； ρ1-2 始末测点的空气平均密度， kg/m 3； 1 2ρ 1v21， 1 2ρ 2v22 风流断面1， 2处断面的平均速压。 3.2.测段风阻的计算 R1-2 hr1-2/ Q21-2 式中， R1-2 断面 1 至 2 间井巷的风阻， kg/m 7， Q1-2 断面 1 至 2 间井巷的平均风量。 3.2.’测段百米风阻的计算 R1-2 （100） 100R1-2/ L 式中，L 测段长度， m 3.2.摩擦阻力系数α值的计算 α1-2 R1-2S31-2/ L1-2U1-2 3.2.自然风压的计算 矿井进回风两侧空气柱作用在底部单位面积上的重力 差就是自然风压。特点是冷而重的空气向下流动， 热而轻 的空气向上流动。矿井自然风压 HN的计算公式如下 HN∮ ρ进风gdz -∮ρ回风gdz 4结论 4.1测算结果 晋普山矿井 1 号风井通风阻力汇总如表 1、 表 2 所示。 表 1 晋普山矿井 1 号风井通风阻力汇总表 节点 测段 节点 静压差 （Pa） 标高差 （m） 平均密度 （kg/m3） 断面积 （m 2） 风量 （m 3 /s） 阻力 （Pa） 累计阻力 （Pa） 风阻 （kg/m7） 长度 （m） 阻力系数 （kg/m3） 121.80.01.061 811.2047.601818.0000.007 94150.057 750 232.50.01.073 211.2547.6025.00043.0000.011 031800.006 759 347.3- 1.51.080 011.3048.5856.80099.8000.024 073700.007 252 452.70.01.081 510.4548.5827.000126.8000.011 44600.017 483 5623.8- 1.51.083 510.1576.47221.748 0348.5480.037 927370.004 386 6123.0- 1.11.086 69.4019.6318.047 4366.5950.046 842650.012 436 1213- 0.60.61.081 78.704.050.490 2367.0850.029 891230.014 089 96 ══ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ 第九期程根 银 等 晋 普 山 煤 矿 通 风 系 统 阻 力 测 定 与 分 析 节点 测段 节点 静压差 （Pa） 标高差 （m） 平均密度 （kg/m3） 断面积 （m 2） 风量 （m 3 /s） 阻力 （Pa） 累计阻力 （Pa） 风阻 （kg/m7） 长度 （m） 阻力系数 （kg/m3） 1317- 1.72.51.077 48.8310.759.935 0377.0200.085 973680.014 059 17185.53.71.088 88.4529.9095.285 6472.3060.106 584690.012 252 1820- 1.71.21.092 38.1031.57- 3.892 4468.413- 0.003 90151- 0.001 250 20211.2- 0.21.092 39.7030.889.815 4478.2290.010 291920.004 081 2122- 3.30.81.092 312.0523.43- 24.262453.967- 0.044 2080- 0.072 330 2223- 1.11.81.092 911.0522.878.672 2462.6390.016 581450.012 055 23252.5- 0.31.093 27.806.2521.720 4484.3600.556 04360.681 681 252624.84.11.092 84.2012.22292.805 0777.1651.960 816460.026 377 26273.5- 15.71.091 34.500.97- 136.330 0640.830- 144.900 048- 31.171 500 2728- 6.67.91.091 54.759.2020.228 5661.0590.238 991360.020 773 282928.9- 20.81.091 34.259.3262.009 6723.0690.713 881400.045 643 293016.9- 7.51.086 83.859.7587.490 0810.5590.920 341640.039 234 3031- 20.720.21.083 94.059.7511.947 8822.5060.125 681200.008 311 31320.54.81.085 63.9010.6357.108 8879.6150.505 40400.091 231 3233- 5.49.71.088 63.6010.7351.594 2931.2090.448 135340.004 960 33348.4- 5.51.089 97.6037.8724.055 5955.2650.016 771690.003 799 343514.3- 2.11.087 98.3039.00120.154 01 075.4200.079 002810.013 412 35368.6- 0.71.086 16.0046.9778.397 31 153.8200.035 541380.005 898 363727.0- 1.41.083 15.6551.47254.837 01 408.6500.096 204700.004 034 373818.7- 7.01.081 47.0053.28111.302 01 519.9500.039 214980.002 651 383990.4- 71.11.078 68.1582.48137.115 01 657.0700.020 161650.006 016 394040.5- 19.01.074 48.2085.68200.864 01 857.9300.027 36650.021 053 40地面1.070 68.2085.68243.0002 100.9300.033 10 表 2 晋普山矿井 3 号风井通风阻力汇总表 节点 测段 节点 静压差 （Pa） 标高差 （m） 平均密度 （kg/m3） 断面积 （m 2） 风量 （m 3 /s） 阻力 （Pa） 累计阻力 （Pa） 风阻 （kg/m7） 长度 （m） 阻力系数 （kg/m3） 019.60369.100 0369.1 123.7- 0.41.088 909.1040.9532.644 4401.7440.019 471720.007 344 230.0- 0.11.089 708.8038.70- 1.089 7400.655- 0.000 7038- 0.001 140 343.2- 0.31.089 808.5534.8028.730 6429.3850.023 721400.009 408 460.7- 1.41.090 208.9541.70- 8.262 8421.123- 0.004 80288- 0.001 030 670.701.090 509.7041.327.000 0428.1230.004 10670.004 658 797.5- 11.090 509.5538.0364.095 0492.2180.044 323500.009 269 9103.7- 0.61.090 4012.2058.4730.457 6522.6750.008 911880.006 399 10113.3- 0.71.089 8012.4047.1225.371 4548.0470.011 431570.010 236 111241.1- 36.81.087 658.4020.7010.744 8558.7910.025 081780.006 925 1213- 9.06.41.085 458.4018.67- 20.531 0538.260- 0.058 90100- 0.028 960 1314- 7.68.41.085 958.4018.6215.219 8553.4800.043 90670.032 209 14152.7- 31.087 305.305.92- 5.619 0547.861- 0.160 30110- 0.022 660 1516- 5.791.088 505.505.1340.965 0588.8261.556 603000.088 485 161711.1- 9.21.090 755.305.1310.651 0599.4770.404 723150.019 973 171816.0- 141.090 605.204.857.316 0606.7930.311 023000.015 367 181955.601.084 555.204.85556.000 01 162.79023.636 90 1920- 11.613.11.081 007.3029.7525.611 01 188.4000.028 941650.006 070 20211.75.71.082 057.7031.1278.676 91 267.0800.081 242350.013 672 21222.6- 0.91.081 857.9035.0216.263 41 283.3400.013 26290.019 282 222314.1- 6.61.080 958.0042.1369.657 31 353.0000.039 243150.005 421 232413.41.61.081 457.7043.63151.303 01 504.3000.079 483800.008 938 242515.85.21.082 957.2075.80214.313 01 718.6200.037 30 252871.1- 62.51.076 1575.8038.406 31 757.0200.006 68 28地面1.067 75111.32821.6002 578.6200.066 30 由于测定方法、 仪器精度、 自然环境及人为因素的影响， 测定结果难免存在一定误差， 但只要把误差控制在一定范围之内 （即 07 中国安全科学学报 ChinaSafetyScienceJournal 第15卷 2005 ══ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ 年 相对误差不超过 5） 即能满足精度要求。 将 1 号风井各通风系统误差分析、 3 号风井各通风系统误差分析和井下通风系统阻力分布情况汇总如表 3～ 表 4 所示。 表 3 1 号风井各通风系统误差分析 系统 名称 系统阻力 （Pa） 阻力平均值 （Pa） 绝对误差 （Pa） 相对误差 （） 105 采区2 100.93 2 113.06- 12.13- 0.57 105 采区2 101.51 2 113.06- 11.55- 0.55 903 采区2 138.82 2 113.0625.761.22 903 采区2 110.99 2 113.06- 2.07- 0.10 表 4 3 号风井各通风系统误差分析 系统 名称 系统阻力 （Pa） 阻力平均值 （Pa） 绝对误差 （Pa） 相对误差 （） 903 采区2 487.11 2 543.83- 56.72- 2.23 903 采区2 459.29 2 543.83- 84.54- 3.32 901 采区2 578.62 2 543.8334.791.37 901 采区2 578.63 2 543.8334.81.37 901 采区2 579.66 2 543.8335.831.41 901 采区2 579.65 2 543.8335.821.41 表 5 井下通风系统阻力分布表 系统名称总阻力 （Pa） 进风段用风段回风段 阻力值 （Pa）占百分比 （）阻力值 （Pa）占百分比 （）阻力值 （Pa）占百分比 （） 105 采区 （1 号风井）2 100.93 472.3122.48458.9021.841 169.7255.68 105 采区 （1 号风井）2 101.51 371.5517.68560.2326.661 169.7355.66 903 采区 （1 号风井）2 138.82 370.0517.3540.3325.261 228.4457.44 903 采区 （1 号风井）2 110.99 370.0517.53512.5024.281 228.4458.19 903 采区 （3 号风井）2 487.11 370.0514.88540.3321.731 576.7363.39 903 采区 （3 号风井）2 459.29 370.0515.05512.5020.841 576.7464.11 901 采区 （3 号风井）2 578.62 522.6820.27640.1124.821 415.8354.91 901 采区 （3 号风井）2 578.63 522.6820.27640.1224.821 415.8354.91 将矿井通风系统中的 3 条主测路线各段风路按其所在 位置及作用， 分为进风段、 用风段和回风段 3 部分， 各段阻力 分布情况如图 3、 图 4 和图 5 所示。 图 3 105 采区测定路线 图 4 903 采区测定路线 图 5 901 采区测定路线 4.2结论分析及对策 根据测算结果可以看出 1）矿井通风系统总体布局合理， 风量分配均匀， 其中 1 号风井回风5142m3/min， 系统总阻力2100.93Pa， 等积孔为 2.225 m2， 井巷风阻为 0. 286 055 N. S2/m7。3 号风井回风 6 678 m3/min， 系统总阻力为 2 578.62 Pa， 等积孔为2.608 m2， 井巷风阻为 0.208 160 N. S2/m7。两个系统的等积孔较大， 通 风阻力较小， 各系统通风容易。 2）矿井大部分巷道断面规整， 维护良好， 保证了正常的 生产用风。 3）从各测定路线阻力分布表可以看出， 井下各系统回 风段阻力基本都超过了总阻力的 55， 所占比例相对偏高。 17 ══ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ 第九期程根 银 等 晋 普 山 煤 矿 通 风 系 统 阻 力 测 定 与 分 析 其原因主要有 一是回风路线较长， 以 903 采区 3 号回风井系 统为例， 系统总长度为 7449 m， 而回风段长度为 3197 m， 占总 长度的 42.92。二是回风路线上设置了大量风门和调节风 门， 消耗了大量的能量。三是回风井内安装了铁栅栏， 增大 了矿井的通风阻力， 从 3 号风井的测定情况看， 井底 3 道铁 栅栏两侧的阻力达到了 215 Pa， 106 回风 2 道铁栅栏两侧的阻 力达到了97Pa。四是工作面顺槽和局部回风巷道断面较小， 例如， 1 号风井系统 105 回风巷的 35 ～ 37 测段 （10 502 掘进回 风 - 炸药库回风） ， 平均断面积不足 6 m2； 105 回风巷的 30 ～ 33 测段 （105 回风 I- 北回 1 号点） ， 平均断面积不足 4.3 m2。 4）从矿井 90 万吨/年的产量看， 两个通风系统的阻力分 别为 2 100 Pa 和 2 450 Pa， 阻力值偏大； 但从系统的实际情况 看， 通过扩大局部巷道断面， 调整调节风门和铁栅栏的设置， 完全有可能把系统的阻力降下来。 （收稿 2005 年 4 月； 作者地址 河北省三河燕郊开发区； 华北科技学院安全工程学院； 邮编 065201） 参 考 文 献 1程根银， 吴怀俊， 金龙哲.晋普山煤矿 3风井主要通风机安全技术参数测定与分析 ［J］ .中国安全科学学报， 2002， 12 （4） 48 ～ 54 2吴中立.矿井通风与安全 ［M］ .徐州 中国矿业大学出版社， 1989 3张国枢.通风安全学 ［M］ . 徐州 中国矿业大学出版社， 2000 4淮南煤炭学院通风安全教研室.矿井通风技术测定及其应用 ［M］ .北京 煤炭工业出版社， 1980 （上接第 56 页 “孟亦飞等 含氧氢气钢瓶释放过程危险性分析及处置” ） 参 考 文 献 1北川彻三 ［日］ .爆炸事故的分析 ［M］ .黄九华译.化学工业出版社， 1984 39 2汤广官.气体静电火灾的预防 ［J］ .安徽消防， 1996 （7） 27 3崔克清.安全工程大辞典 ［M］ .北京 化学工业出版社， 1995 58 4王凯全， 邵辉等.事故理论与分析技术 ［M］ .北京 化学工业出版社， 2004 127， 27 5W.Bartknecht ［德］ .爆炸过程和防护措施 ［M］ .何宏达译.北京 化学工业出版社， 1985 3， 18， 57 6周乐.化工厂静电火灾的防治 ［J］ .化工设计通讯， 2001， 27 （2） 48 ～ 49， 52 7白光弼， 张朝临.加油站静电灾害的形成与防护措施 ［J］ .陕西气象， 2001 （3） 30 ～ 32 8王军.油库、 加油站静电引燃爆炸危险及防范 ［J］ .电气防爆， 2003 （2） 1 ～ 7 9李荣华， 曲春香.静电危害的防治 ［J］ .科技情报开发与经济， 2002 （12） 118 ～ 119 10沈维道等.工程热力学 ［M］ .北京 高等教育出版社， 1998 286 11曾鉴荣， 于泊洋.氢气爆炸能量工程估算方法 ［J］ .爆轰波与冲击波， 1999 （3） 113 ～ 116 12国家安全生产监督管理局.安全评价 ［M］ .北京 煤炭工业出版社， 2004 399 27 中国安全科学学报 ChinaSafetyScienceJournal 第15卷 2005 ══ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ 年