可燃性气体在工业和生活中得到了广泛应用，但由于生产、储存、输配、使用等环节中设备老化、维护滞后及人为操作等因素，气体爆炸事故时有发生，且多发生在厂房、仓库及民用住宅等建筑屋室内部，造成的后果十分严重[1]。泄爆是此类气体内爆炸事故的有效防护措施，在管道、容器等方面研究较多。例如：胡俊等[2-3]进行了大长径比柱形容器开口泄爆过程的实验研究，探究了不同开口面积以及不同泄爆压力对压力与火焰速度的影响；Alexiou等[4-5]和Kamani等[6]对圆形管道内的燃爆泄爆过程进行了相应的实验研究,研究了不同点火位置及不同泄爆位置对压力的影响；周灿等[7]、师喜林等[8]对球形容器泄爆过程开展了大量的研究工作。对于大空间泄爆过程，Chao等[9]分别在63.7和2.32 m3空间内进行泄爆实验，分析了不同类型可燃气体、不同泄压面积、不同体积以及障碍物等因素对最大压力的影响。与圆柱形管道或球形容器不同，常见建筑屋室多为容积较大的方形空间，内部存在障碍物，爆炸时门、窗等构件可破裂形成泄压口，其位置及泄压面积任意，加上气体类型、浓度、点火位置等影响，泄爆过程十分复杂，研究时应当根据实际需求，分析其中主要的一种或几种因素对空间内压力和火焰发展的影响。当前在此方面的实验研究不多，如文献[9-11]仅得到了泄爆口外部火焰传播特性和压力特性，内部火焰传播特性难以从实验手段获取，而空间内部火焰传播与压力变化的相互作用关系也不容忽视[2]。常用的解析模型未考虑动力学问题，假设整个空间内部压力任何时刻都均匀分布，与实际工况相差较大[12]。现阶段CFD研究取得了显著进展，肖华华[13]、李润之等[14]借助Fluent软件研究火焰形状，得到了与实验非常相近的结果。数值模拟研究可以详细揭示火焰传播特性，有助于理解泄爆过程。

本文中，通过实验研究，分析1.21 m3方形空间内泄压面积和气体体积分数对泄爆压力的影响，在此基础上，针对7%体积分数乙烯-空气混合气体泄爆过程进行数值模拟，分析泄爆过程中火焰阵面发展及火焰传播速度对泄爆过程空间内部压力变化的影响，为方形空间的工程泄爆防护提供参考。

式中，ρave为销蚀碎片形成的圆管的平均密度；Rm为弹体蘑菇头半径；r0为弹体初始半径。Rm与r0之比为[22]

1 实 验

1.1 实验装置与方法

图1 实验装置图 Fig.1 Schematic of experimental setup

图1为自制的方形泄爆实验装置，由方形腔体、混气系统、点火系统、高速摄像系统、压力检测与数据采集系统等组成。方形腔体一端封闭、一端完全开口，内部空间尺寸为2 000 mm×1 100 mm×500 mm，容积为1.21 m3。腔体一侧沿轴线每隔300 mm均匀开设6个的进气孔，安装进气管道，中间开设1个抽气孔。采用气体分压原理，参照文献[15]的方法，混气时根据所需配置乙烯体积分数先抽取部分空气，然后通过开设多个气孔的进气管道充入乙烯气体至标准大气压，进气过程的局部湍流将加速直管气孔周边空间内的气体混合，进气完成后静止10 min。点火装置设置于封闭端，采用直流电源加热电阻丝的方式点火。压力传感器选用PCB113B24，采样频率为200 kHz，分别位于腔体顶部与点火端相距600 mm和靠近泄爆口100 mm处，即图1中的P1和P2处；高速摄像机记录泄爆过程，拍摄频度为1 000 s-1。

在腔体开口端安装两种泄压面积不同的泄压装置：第1类泄压装置为聚乙烯薄膜，尺寸为300 mm×600 mm，面积为0.18 m2；第2类为硅酸钙板，尺寸为1 200 mm×600 mm，泄压面积为0.55 m2。实验时，分别充入体积分数为3%、4%、5%、6%、7%、8%的乙烯气体。

由此可见，空间内部压力发展与火焰传播速度变化密切相关。火焰传播速度越快，化学反应越剧烈，单位时间产生的已燃气体体积越大，内部压力上升速率越快，将产生更大的压力峰值。泄压口开启后，气流泄放对空间内部的扰动引发爆燃火焰速度的加快，随着火焰前锋接近泄爆口，泄压口泄放效应对火焰传播影响更显著，导致火焰阵面变形更大，火焰传播速度也更大，因而产生了两次升压过程和的两个压力峰值。

1.2 实验结果

图2～3分别为不同气体体积分数时、两类泄爆条件下靠近泄爆口处测点的压力曲线。由图可见，泄爆压力曲线均出现两个峰值。点火后，空间内部压力上升使得泄压装置作用，泄放气体介质，产生第1个压力峰值，泄压面积不同时，第1个压力峰值均在10 kPa左右，且随体积分数变化不明显，而空间内部升压速率随乙烯体积分数增大而增大。随后，压力经过了先下降后上升的过程，并形成第2个峰值：在小泄压面积时，压力下降幅度随着气体体积分数增大而减小，第2个压力峰值则随着气体体积分数的逐渐增大而增大，在6%～8%体积分数时，泄爆压力第2峰值超过第1峰值，体积分数为8%时，第2个压力峰值达到26 kPa；在大泄压面积时，压力下降幅度随着气体体积分数变化不大，均下降至-4kPa左右，相比小泄压面积时低，第2峰值压力均较小且随体积分数变化不明显。

图2 0.18 m2泄压面积时不同体积分数压力曲线 Fig.2 Overpressure curves for different volume fractions with 0.18 m2 vent area

图3 0.55 m2泄压面积时不同体积分数压力曲线 Fig.3 Overpressure curves for different volume fractions with 0.55 m2 vent area

2 数值模拟

2.1 数值模拟方法和初始条件

图4 计算网格 Fig.4 Numerical mesh

对泄压面积为0.18 m2、体积分数7%乙烯-空气混合气体爆燃泄爆过程进行数值模拟研究。使用GAMBIT软件对实验装置进行三维建模，为简化计算，模拟过程中忽略容器中进气管路的影响，根据对称性取二分之一模型并划分六面体网格，如图4所示。其中：方形空间为12.5 mm的三维均一结构网格，泄爆口外部流场为大小呈梯度变化的结构网格，总网格数为629 760个；壁面边界条件为绝热边界条件，泄爆口初始时设置为壁面边界类型，假设内部压力到达开启压力10 kPa后泄爆口瞬间完全开启，设置为内部边界类型。

RNG k-ε湍流模型的方程为：

计算选用的CFD软件为Fluent 17.0，采用压力基瞬态求解器，基于RNG k-ε湍流模型和EDC燃烧模型，压力-速度耦合求解算法采用SIMPLEC算法。质量方程、动量方程、组分方程和能量方程的控制方程为：

(1)

(2)

(3)

(4)

（2）提高任课教师的实践水平，具体可以通过“走出去和引进来”方法实现。走出去就是鼓励教师参加社会实践，让任课教师到银行国际结算业务部或比较大型的外贸企业进行参观、实践。将自身具备的专业知识与实际业务相结合，获取最新的业务做法，将前沿实践操作带进课堂。为了鼓励教师积极进行调研、培训、挂职锻炼，建议高校对教师参加实践和课堂授课给予同等待遇，同时将实践能力作为教师职称评定、岗位竞聘、评优评先等晋升的条件之一，从而达到调动教师参与实践的积极性，全面提高专业教师实践水平的目的。

(5)

(6)

(7)

初始压力为标准大气压，容器中完全静止，各方向速度分量为零，湍流动能为1 m2/s2，湍流耗散率为1 m2/s2，容器中混合气体体积分数分别为乙烯0.07，氧气0.024，其余均为氮气，初始温度为300 K。点火区域设置于容器一端中心处，点火半径为5 cm，初始温度为1 200 K，假设点火时乙烯质量分数的一半发生反应。压力监测点位于容器内靠近泄爆口的上壁面，与实验中的P2压力传感器位置相同。

图5为小泄压面积(0.18 m2)、7%体积分数时数值模拟结果与实验中测得的压力曲线的对比。

式中：为密度，p为压力，ui、uj为速度分量，t为时间，μ为动力黏性系数，Ys、ωs、hs、Js分别为组分s的质量分数、反应率、焓和扩散流，Sh为反应热，keff为热导率，(τij)eff为应力张量，E为总能，k为湍流动能，ε为湍流耗散率，Gk是由层流速度梯度而产生的湍流动能，Gb是由浮力产生的湍流动能，YM是由在可压缩湍流中过渡的扩散产生的波动，各常数取值Cμ=0.084 5，η0=4.38，β=0.012，αk=αε=1.393，C1ε=1.42，C2ε=1.68。

2.2 数值模拟结果

采用的化学反应为总包反应，化学反应方程为：C2H4+3O2→2CO2+2H2O。

图5 实验与数值模拟的压力曲线 Fig.5 Overpressure curves between experimental and numerical results

图7～8分别为方形空间气体泄爆过程中火焰速度与火焰前锋位置和爆燃压力的关系，其中压力为滤波后的曲线。在密闭爆燃阶段(t<78 ms)，气体燃烧缓慢，火焰速度小于10 m/s且振荡上升，压力上升速率也较慢且呈现一定的波动，火焰速度的震荡可能是火焰阵面受到空间内部压力波作用的影响。随着火焰阵面的发展，火焰传播速度不断增大，当t=78 ms时泄压口开启，气体泄放引起的扰动开始影响火焰前锋，尽管图6表明泄爆前后火焰阵面形状并没有明显变化，但计算表明火焰速度在3～4 ms内先急剧增大到40 m/s后再突然减小，形成一个峰值，该时段内火焰前锋位置也经历一个先变化较快后减缓的过程(见图7)，而压力同样产生一个较小峰值(见图8)。产生该现象的原因可能为，泄压口开启后空间内气体瞬间泄放使得火焰加速向前传播，遇到冷的预混气体后，燃烧反应速率变慢，火焰传播速度减慢。随后，受泄放影响火焰传播速度又振荡上升，压力也产生了较大幅度振荡；随着火焰阵面逐渐接近泄压口，它受到的拉伸作用越加显著，火焰传播速度不断增大，火焰阵面前端受拉伸变形呈现水滴型，在火焰前锋到达泄压口时达到最大速度90 m/s，此时空间内部压力也达到第2个峰值；此后，火焰传出泄压口，压力逐渐回落。

随着我国经济快速发展，能源危机日益临近，而传统能源对环境的危害也日益显著，所以急需开发利用新能源来降低对传统能源依赖，并减少环境污染。地热能是由地壳抽取的天然热能，具有清洁无污染、能量储存量大等特点，受到越来越多国家的关注。国家能源局也提出，“十三五”时期，我国将全方位推进地热能发展。

3 结果与分析

3.1 泄爆过程火焰面形态变化

图6为方形空间内部不同时刻火焰发展情况。由图可见：在初始点火时刻到t=40 ms时，火焰呈半球形向外扩展；当t=78 ms时，火焰阵面接近空间上下壁面，火焰形状呈扁平状，尽管已燃气体体积较小，但空间内部的压力已达到泄爆开启压力10 kPa，火焰形状主要由于壁面约束而呈扁平状并向前方和两侧发展；当t=100 ms时，火焰阵面接近空间两侧壁面，其后由于四周壁面限制仅向泄压口方向传播，火焰阵面前端呈扁球形，此时火焰阵面距泄压口较远，气流泄放对火焰形状的影响不大；当t=130 ms时，火焰阵面前端接近泄爆口，由泄爆产生的扰动效应使火焰阵面前端开始出现拉伸变形，并逐渐呈水滴形；当t=140 ms，火焰阵面前端变形加大并被拉伸至泄爆口，其后从泄爆口中心位置传出，在泄爆口外逐渐形成蘑菇状的火焰前端，这与实验中观察到外部火焰形态基本一致。上述分析表明，本文数值模拟方法可合理地描述方形空间内部气体爆燃火焰的发展过程。

3.2 泄爆过程火焰传播速度特性

由图可见：在泄爆前密闭空间爆燃阶段，两者压力曲线较为一致；泄压口开启后，计算压力和实验压力均呈现先下降后上升再下降的趋势，但计算压力整体小于实验压力，且计算压力的下降速率要大于实验压力，而上升速率和第2个压力峰值小于实验压力。主要原因在于，数值模拟中假设泄爆口瞬时完全开启，未考虑泄爆口的开启过程，而在实验中，泄爆口的开启虽然短暂但仍需要一定的时间。该时段实验泄爆效果不如数值模拟结果。此外，实验装置内部含有6根进气管道，尽管阻塞率较小，但在火焰传播过程中仍会造成火焰速度加快，压力上升速率增大，在泄压面积一定时，将产生更大的压力峰值，数值模拟中未考虑其影响也是造成计算与实验压力误差的原因之一。上述分析表明，该数值模拟方法能够较合理地描述方形空间内气体爆燃泄放过程及其压力发展趋势。

图6 泄爆过程火焰传播图 Fig.6 Image of premixed flame propagation during explosion venting

图7 火焰阵面前端位置与速度 Fig.7 Flames speed and flame front position

图8 压力和火焰速度 Fig.8 Flames speed and overpressure

初始碳排放的探讨。初始碳排放就是二氧化碳自发排放数量，当碳交易年度总成交金额和地区生产总值均等于零时，碳排放的初始值应该是一个正的常数，代表居民自发碳排放。

4 结 论

通过对方形空间内不同体积分数乙烯-空气预混气体在两种不同泄压面积下的爆燃泄爆过程压力和火焰特性的实验和数值模拟研究，可以得出以下结论。

(1)在泄爆压力为10 kPa左右时，大泄压面积和小泄压面积条件下压力曲线均出现双峰。小泄压面积时，随着气体混合物中乙烯体积分数增大，升压速率越快，第2峰值越大，在6%～8%体积分数时，第2峰值超过第1峰值。而在大泄压面积时，泄爆后压力下降幅度较大，第2峰值较小，且第2峰值随体积分数增加变化不大。

(2)数值模拟中，假设泄爆口瞬时完全开启，未考虑泄爆口开启过程，与实验中泄爆口非瞬时开启的情况相比，泄爆后压力下降较快且幅度更大，第2超压峰值小。因此，尽量减小泄爆装置的开启时间，对于泄爆后的二次升压过程有一定的削弱作用。

2.历史因素——资源型城市对石油资源具有的依赖性。东营市是依靠油田开发建设而成立、发展、壮大起来的，东营市本埠产业结构比较单一，财政收入渠道狭窄。多年来，油区内的各级政府和村民存在“靠山吃山，靠油吃油”的思想，产业结构调整缓慢，地方经济发展滞后，“吃拿卡要”等成为部分地区的一种习惯。一些村民甚至认为，油田是在自己的地界，拿一点、用一点理所应当，促使某些人明目张胆地盗窃、破坏油田的生产物资，谋取私利，给油田正常的生产造成了极大的影响和破坏。

(3)火焰传播速度与压力发展存在对应关系，火焰传播速度的变化是压力产生波动的主要原因。泄爆口开启后所引发的湍流效应将导致空间内火焰速度明显增大，火焰前锋越接近泄爆口，湍流效应对火焰的影响越显著，表现为火焰变形增大，火焰阵面前端呈水滴形，传播速度加快。火焰传播速率的增大导致压力上升速率的增加，当火焰传播速率增至最大时，产生了第2个压力峰值。

除了“日月阴阳互藏其宅”的说法外，也有从日、月的“阴阳交感”来解释太阳在东，但“酉鸡”在西，月亮在西，但“卯兔”在东的现象。如明人郎瑛在《七修类稿》中有一段关于十二生肖的解释：

参考文献：

[1] 叶宏.民用燃气爆炸及对建筑结构影响的分析与研究[D].北京:清华大学,1994.

[2] 胡俊,浦以康,万士昕,等.柱形容器开口泄爆过程中压力发展特性的实验研究[J].爆炸与冲击,2001,21(1):47-52.

HU Jun, PU Yikang, WAN Shixin, et al. Experimental investigations of pressure development during explosion venting from cylindrical vessels[J]. Explosion and Shock Waves, 2001,21(1):47-52.

[3] 胡俊,万士昕,浦以康,等.柱形容器开口泄爆过程中的火焰传播特性[J].爆炸与冲击,2004,24(4):330-336.

HU Jun, WAN Shixin, PU Yikang, et al. The characteristics of flame propagation during explosion venting from cylindrical vessel[J]. Explosion and Shock Waves, 2004,24(4):330-336.

[4] ALEXIOU A, ANDREWS G E, PHYLAKTOU H. Side-vented gas explosions in a long vessel: The effect of vent position[J]. Journal of Loss Prevention in the Process Industries, 1996,9(5):351-356.

[5] ALEXIOU A, ANDREWS G E, PHYLAKTOU H. A comparison between end-vented and side-vented gas explosions in large L/D vessels[J]. Process Safety and Environmental Protection, 1997,75(1):9-13.

[6] KASMANI R M, ANDREWS G E, PHYLAKTOU H N. Experimental study on vented gas explosion in a cylindrical vessel with a vent duct[J]. Process Safety and Environmental Protection, 2013,91(4):245-252.

[7] 周灿,王志荣,蒋军成.球形容器小口泄爆压力变化特性[J].实验流体力学, 2013,27(4):62-65.

ZHOU Can, WANG Zhirong, JIANG Juncheng. Change characteristics of gas explosion venting pressure inside spherical vessel with a small caliber[J]. Journal of Explosions in Fluid Mechanic, 2013,27(4):62-65.

[8] 师喜林,王志荣,蒋军成.球形容器内气体的泄爆过程[J].爆炸与冲击, 2009,29(4):390-394.

SHI Xilin, WANG Zhirong, JIANG Juncheng. Explosion-vented processes for methane-air premixed gas in spherical vessels with venting pipes[J]. Explosion and Shock Waves, 2009,29(4):390-394.

[9] CHAO J, BAUWENS C R, DOROFEEV S B. An analysis of peak overpressures in vented gaseous explosions[J]. Proceedings of the Combustion Institute, 2011,33(2):2367-2374.

[10] MERCX W P M, WINGERDEN C J, PASMAN H J. Venting of gaseous explosions[J]. Process Safety Progress, 1993,12(1):40-46.

[11] BAO Q, FANG Q, ZHANG Y, et al. Effects of gas concentration and venting pressure on overpressure transients during vented explosion of methane-air mixtures[J]. Fuel, 2016,175:40-48.

[12] HAN Yongli, CHEN Longzhu. Mechanical model of domestic gas explosion load[J]. Transactions of Tianjin University, 2008,14(6):434-440.

[13] 肖华华.管道中氢-空气预混火焰传播动力学实验与数值模拟研究[D].合肥:中国科学技术大学,2013.

[14] 李润之,司荣军.低温环境下甲烷爆炸流场特性模拟[J].爆炸与冲击, 2015,35(6):901-906.

LI Runzhi, SI Rongjun. Simulation study of flow field characteristics of gas explosion in low temperature environment[J]. Explosion and Shock Waves, 2015,36(6):901-906.

[15] 樊保龙,李斌,白春华,等.大尺寸密闭空间内甲烷-空气二元混合方法的改进与评定[J].高压物理学报,2016,30(3):242-248.

FAN Baolong, LI Bin, BAI Chunhua, et al. Improvement and evaluation on binary mixing method of methane and air in large-scale confined space[J]. Chinese Journal of High Pressure Physics, 2016,30(3):242-248.