泄爆和氮气惰化耦合作用对氢-空气爆炸影响的实验研究*

张 凯，杜赛枫，陈 昊，郭 进，王金贵，洪溢都

（福州大学环境与安全工程学院，福建 福州 350116）

与石油和煤炭等化石燃料相比，氢是一种清洁高效的能源载体。目前，氢能的开发利用已呈可持续发展的趋势。但是，氢气具有可燃范围广、点火能量低、燃烧速度快等特点，当它在密闭空间内积聚，遇到意外火源时极易爆炸并产生超压，会对临近人员的安全和建筑结构的完整性构成威胁。如2011 年3 月11 日日本福岛核事故和2018 年12 月26 日我国北京交通大学实验室爆炸事故，主要就是由氢气爆燃引起的。因此，在以氢能为基础的工业中，控制氢气爆炸已经成为一个急需解决的重要安全问题。

泄爆和抑爆是两种常用的防爆减灾措施[1-2]。泄爆是通过预设低强度材料为泄爆口，在爆炸早期提前泄压；抑爆则是通过加入惰性物质来降低受限空间内爆炸性混合物中的氧浓度，抑制其反应性，减弱爆炸超压。目前，学者们已通过大量实验和数值模拟研究了泄放参数[3-7]（如燃料浓度、泄爆面积、点火位置、破膜压力、障碍物和泄放容器等）对泄爆过程的影响，并阐明了爆炸容器内不同压力峰值、亥姆霍兹振荡和声学振荡及外部爆炸的形成机制[8-11]。近年来，学者们也大量研究了各种惰化剂的性质和浓度对气体爆炸的抑制效果[1,12-14]。这些研究表明，单独采用泄爆或惰化都可有效降低最大爆炸超压。那么，这两种方式的结合能否进一步降低爆炸超压，值得深入研究。

崔益虎等[2]和Zhang 等[15]分别在不同的容器中研究了泄爆和CO2惰化共同作用对甲烷-空气混合物爆燃的影响，证实与单一方法相比，二者的耦合作用可以显著降低爆炸超压。这个结论也在泄爆与细水雾耦合[16-17]或同时使用泄爆与细水雾、氮气[18-19]对氢-空气爆炸特性影响的实验及数值模拟中得到验证。此外，Wen 等[20]在有障碍物的泄放管道中研究了泄爆和超细水雾协同作用对氢-甲烷-空气混合物爆炸特性的影响；Pei 等[21]则指出，泄爆与氮气和超细水雾的联合应用也可以显著提高对管内瓦斯爆炸超压的抑制效果；Ingram 等[22]还联合使用氢氧化钠细水雾和泄爆的方式来抑制H2/O2/N2预混气体的爆炸过程。刘洋等[23]利用FLACS 软件，模拟了在末端点火时长30 m 管道内不同浓度CO2、N2和水蒸气分别与泄爆耦合作用对甲烷-空气混合物的抑爆效果，发现当这3 种气体的体积浓度分别为25%、26%和30%时能完全抑爆和抑制超压振动。张迎新等[24]分别讨论了泄爆与3 种不同体积分数的N2和CO2耦合对瓦斯爆炸超压的影响，实验指出瓦斯爆炸超压会随着惰性气体添加量的增加而减小。Lu 等[25]在泄爆口附近及下游位置，通过喷射氮气抑制甲烷-空气爆炸发现，氮气的喷出压力决定能否防止泄爆过程中火焰的传播。对泄爆和其他惰性剂如七氟丙烷[26]和碳酸氢钠粉末[27-28]等耦合在抑爆方面的效果，也已有大量的研究。

以上表明，关于泄爆和惰化对气体爆燃耦合效应的研究主要集中在碳氢化合物-空气混合物，而对反应性和扩散性更强的氢气的研究[16-19]却十分有限；而且，之前研究人员注重于分析二者耦合作用对最大超压峰值及其上升速度的影响，而对抑爆过程中火焰的演化、内部超压不同峰值和振荡的产生机制、外部爆炸的形成及其对内部压力的影响等却较少涉及。此外，使用水雾或同时水雾和氮气作为惰性剂来抑制氢-空气的爆燃，在一些特殊场所中是不允许的。因此，有必要进一步研究惰性气体如氮气，是如何缓解氢-空气爆炸强度的。

本文中，在顶部泄放的矩形容器内，进行泄爆和氮气惰化对内部火焰行为、外部爆炸形成和内外压力曲线耦合效应的实验，以期阐明氮气添加量对恒定当量比氢-空气混合物泄爆特性的影响。

1 实 验

实验是在高1 m、截面0.3 m×0.3 m 的竖直爆炸容器中开展的，如图1 所示。泄爆口位于容器顶部中央，其面积为0.25 m×0.25 m；用法兰固定的铝膜作为泄爆封口，其静态破膜压力为(8±1) kPa。爆炸室中央有长0.75 m、宽0.3 m 的透明观察窗，高速摄像机可通过该窗口以频率2 kHz 记录火焰图像。内外压力分别由4 个量程为0～150 kPa 的压阻式压力传感器PS1～PS4 测量。其中，PS1～PS3 分别安装在爆炸室内部距离泄爆口875、500 和125 mm 处；PS4 则安装在爆炸室外，距离泄爆口350 mm。预混气体在爆炸室中心被能量约500 mJ 的电火花点火器点燃。

图1 矩形泄放容器Fig. 1 Rectangular vented vessel

在对爆炸容器配气前，先用盲板将泄爆口密封，再用真空泵将容器抽至真空状态，确保容器内压没有波动，通过顶部装有压力表的配气柱，将氢气、氮气和空气分别按道尔顿分压定律的特定分压充入容器内。配气结束后，将盲板沿着容器顶端水平移开，同时沿着盲板挪开的方向密封一片铝膜，然后在铝膜上方用螺栓固定一块方形法兰。

混合气体制备完成后，为确保其均匀性，在点火前保持静止15 min。最后，由同步控制器输出信号，同时触发点火系统，高速摄像机和数据采集系统记录数据。所有实验均在初始压力101 kPa 和初始温度285 K 下进行。在实验中，混合物的当量比均为1.6，而氮气体积分数不同。

2 结果和讨论

2.1 泄爆过程中的典型火焰

在实验中，不同氮气体积分数φ下的火焰演化行为相似，但亮度会随着φ的增加而变得黯淡。以φ=20%为例，说明泄爆过程中爆炸室内外典型的火焰行为，如图2 所示。

图2 当φ=20%时容器内外火焰演化过程Fig. 2 Flame evolution inside and outside the vessel for φ=20%

图2(a)为点火11 ms 后泄爆封口破裂时的火焰图像，此时火焰的上锋面仍远离泄爆口，因此未燃烧的氢-氮-空气混合物先排出容器，并在外部形成可燃气云。火焰也在泄放气流作用下被拖向泄爆口，如图2(b)所示。此外，图2(b)中火焰的下锋面较图2(a)的更高，这个现象在φ更高时更加明显。当内部火焰从泄爆口喷射出来后，外部的可燃气云随即被点燃，并由此触发外部爆炸，该过程已经被通过纹影系统和示踪粒子方式的研究[10-11,29]验证。外部爆炸过程的独特现象是，在很短时间内泄爆口附近形成了一个明亮的火球。

外部爆炸发生后，内部火焰的下锋面开始上下剧烈地跳动，见图2(d)～(f)。根据文献[8,30]，这是亥姆霍兹振荡现象。当火焰向下跳跃，即火焰内部低密度燃烧产物向容器底部高密度未燃烧氢-氮-空气混合物加速时，就会在火焰下锋面触发形成瑞利-泰勒不稳定现象[30-31]，火焰下锋面的泰勒不稳定现象将会导致火焰在容器内部的燃烧速率和燃烧面积增大[30]，如图2(f)所示。随着火焰的传播，容器内外压差逐渐减小，在阻力的作用下火焰下锋面上下跳动的振幅随时间逐渐减小，当火焰下锋面接近容器底部时，几乎无法分辨。与此同时，容器外的火球也逐渐变形为射流结构火焰。

图3 为火焰上下锋面在点火后的位置和传播速度。位置是指火焰锋面与点火电极之间的距离；火焰的速度是指火焰锋面位移与时间之比，负值为朝向容器底部的运动。火焰位置和速度可以通过高速摄像机的火焰图像计算。考虑火焰亮度对数据精度的影响，分析φ=0%时火焰上下锋面的速度和位置在泄爆过程中的演变。

图3 当φ=0%时点火后火焰锋面的位置和速度Fig. 3 Flame front locations and velocities after ignition for φ=0%

在破膜前，由于浮力的作用[32-33]，火焰向上膨胀的速度较快。需要说明：由于火焰速度随着爆炸容器内φ的增高而减小，因而浮力对内部火焰的影响在φ较高时会较突出；破膜后，火焰上锋面在泄放气流作用下向泄爆口明显加速，同时下锋面也朝着泄爆口的方向运动。当爆炸室外火球开始形成时，火焰上锋面有最大速度约253 m/s，随后迅速下降。此外，在外部爆炸发生后不久，火焰下锋面开始剧烈振荡，其最大速度至171 m/s、最小速度至-147 m/s。这种振荡与内部超压的亥姆霍兹型振荡同时发生，将会在下节详细讨论。

2.2 氮气体积分数对内部超压的影响

图4 为氮气体积分数从0%增加到50%的内部压力曲线，φ显著影响内部压力演变过程。当φ＜20%时，内部压力曲线与φ=10%的类似，可见到两个压力峰值p1、p2和亥姆霍兹振荡。p1由泄爆膜破裂而产生[7,34]，p2则是外部爆炸和火焰下锋面泰勒不稳定耦合作用的结果[5,35]。外部爆炸降低了泄爆口内外压力梯度，火焰下锋面出现的泰勒不稳定则增加了火焰面积，这些因素的耦合将使内部压力快速上升。

有学者认为，燃烧气体泄放会造成爆炸容器内的体积流出率急剧增加，导致内部压力迅速下降而形成压力峰值p2，其出现的时间通常与内部火焰到达泄爆口时相对应[3,6]。但火焰传播速度很快，使该峰值在很多情况下无法与外部爆炸导致的内部超压正确区分，尤其是小尺度实验容器。在本实验中，当φ=10%, 25%时，结合同步的火焰图像可发现，p2出现的时间始终滞后于火焰到达泄爆口处的时间；当φ=25%时，火焰到达泄爆口时并未导致内部压力峰值，且外部压力峰值甚至早于p2出现（下文将展示）。这些现象表明，p2在火焰冲出爆炸容器并触发外部爆炸之后出现。

随后，在压力曲线出现了亥姆霍兹振荡，它由外部爆炸发生后、内部火焰下锋面持续上下跳动所导致。它如同活塞，以很高速度周期性地压缩容器底部的未燃气体。在所有实验中，都可以观察到亥姆霍兹振荡，当φ＞20%时更易辨识。此外，这种活塞式运动会随着时间逐渐衰弱，因此亥姆霍兹振荡的振幅会逐渐减小。

当φ=25%, 30%时，内部压力曲线上出现了频率约1.25 kHz 的声学振荡和第2 个压力峰值p3。由文献[3,5,9]，这种类型的振荡由混合气体在容器底角的声学增强燃烧引起。在有些情形中，p3甚至成为最大内部压力峰值[9,30,36]，但在本文中，p3始终小于p2。当φ＞30% 时，声学振荡消失。声学振荡只在φ=25%, 30%时才出现，说明燃烧过程和声波之间微妙的耦合作用只出现在特定情形。

由图4 还可以发现，φ能显著影响亥姆霍兹振荡的频率。这是因为，爆炸室声速会随着φ的增加而减小[8]，导致叠加在火焰下锋面活塞式运动上的压力振荡，其频率会从φ=0% 时的220 Hz 减小到φ=50%时的130 Hz。

图4 容器内部压力曲线Fig. 4 Internal pressure curves of the vessel

此外，当φ＜50%时，p2始终是最大的内部压力峰值，而当φ=50%时，p1成为最大压力峰值。由此可以推断，当φ≥50%时，内部最大爆炸超压将与φ无关，因为p1近似等于破膜压力，而后者在燃烧速率低的混合气体泄爆中几乎不会变化[30]。

2.3 氮气体积分数对最大内部超压的影响

图5 为爆炸容器内部最大超压和氮气体积分数的关系。显然，对于给定的φ，3 个压力监测点的pmax存在较大差异。除φ=0%，在其他实验中，pmax都随着压力传感器与泄爆口距离的增加而增加。而对φ=0%时的情形，可以由图6 解释：点火约13 ms 时，PS3 的压力曲线出现了一个持续时间很短的压力尖峰，几乎与外部压力峰值同时出现，这说明该压力尖峰是由外部爆炸产生的冲击波引起的。由于爆炸容器内部的压力取决于燃烧的体积生成速率和泄放的体积流出速率。显然，当φ=0%时，外部爆炸最剧烈，在此情形下外部爆炸产生的冲击波最大限度地限制了爆炸室内部燃烧产物的体积流出速率，使PS3 的最大超压大于PS2 的。

图6 当φ=0%时容器内外压力曲线Fig. 6 Internal and external pressure curves of the vessel for φ=0

由图5 中还发现，PS1～PS3 的pmax的差值随着φ的增加而减小，且当φ增加到50%时，压差几乎减小为零。在相同的压力监测点，pmax总体上随着φ的增加而减小。如，PS1 的pmax先从φ=0%时的145 kPa迅速减小到φ=35%时的19 kPa，后又缓慢减小到φ=40%时的12 kPa，最后随着φ进一步增加到50%，pmax几乎不再减小。pmax随着φ的增大而减小的原因是，φ增加会降低内部氢-氮-空气混合物的燃烧速率[37]。

图5 容器最大内部超压和氮气体积分数的关系Fig. 5 Relations between maximum internal overpressures of the vessel and nitrogen addition ratios

2.4 氮气体积分数对外部爆炸的影响

如上所述，外部爆炸由爆炸室外可燃云爆燃导致。外部爆炸会在容器外产生一个压力峰值，典型的外部压力曲线如图7 所示。当φ＜40%时，外部压力曲线与φ=20%的类似。其中，泄爆封口破裂会产生一个很小的压力峰值。之后，由于外部爆炸，外部压力曲线上会出现另外一个压力峰值pext。而在本文实验中，当φ更高时，可在外部压力曲线中观察到两个压力峰值pext,1和pext,2，外部压力曲线与φ=40%的类似。这可通过容器外火焰的行为解释，而后者又与可燃云团的形状密切相关，下段作具体解释。

图7 当φ=20%, 40%时容器外部压力曲线Fig. 7 External pressure curves of the vessel for φ=20%, 40%

Proust 等[11]通过在预混气体中添加示踪粒子（细二氧化硅颗粒）的方式，揭示了泄爆口外可燃云的形状会从蘑菇状变成射流结构。当φ＜40%时，未燃气体会在泄爆口附近形成蘑菇状可燃云。当蘑菇状可燃云被点燃时，会形成单个火球并产生单个压力峰值。然而，由于容器内混合气体的燃烧速度会随着φ的增加而降低，这会使火焰在容器内的传播时间延长，为可燃气云从蘑菇状转变为射流结构提供更长的时间。因此，当射流结构的可燃云被点燃时，会在几毫秒内在外部流场中形成两个明显的火球，同步火焰图像如图8 所示，由此产生的两个外部压力峰值也能在外部压力曲线中区分。

图8 当φ=40%时容器外部火焰演变过程Fig. 8 External flame evolution of the vessel for φ=40%

图9 为不同氮气体积分数的最大外部超压峰值。由文献[6,10]，pext取决外部可燃云团的体积、形状、湍流水平和反应性等。在本文实验中，pext随着φ的增加而单调减小的趋势，主要由外部可燃云团活性的改变而导致。后者随着φ的增加而降低，这同样可通过外部爆炸发生时火球亮度的变化中得到证实。如图10 所示，火球在φ=0% 时非常明亮，但随着φ的增加，火球的亮度明显变暗，当φ增加到50%时，外部火球甚至无法辨别出来。

图9 容器最大外部超压与氮气体积分数的关系Fig. 9 Relations between maximum external overpressures of the vessel and nitrogen addition ratios

图10 不同氮气体积分数下最大外部超压峰值时的火球Fig. 10 Fireballs at maximum external overpressure under different nitrogen addition ratios

由文献[5,35,38]，外部爆炸会以降低泄爆口内外压力梯度的方式来影响爆炸室内的泄放过程。在本文实验中也存在这种现象，且与外部爆炸的强度无关。如图6 所示，当φ=0%时外部爆炸会使容器内部出现一个振幅很大的压力峰值。虽然φ的增加能够显著削弱外部爆炸的强度，如图11 所示，当φ=25%,40%时，容器内的压力峰值仍几乎与外部压力峰值同时出现。即使当φ=40%时，虽然两个外部压力峰值只有几千帕，也能清晰地从PS3 的压力曲线中区分相应的内部压力峰值（见图11(b)）。

图11 容器外部爆炸对内部压力的影响Fig. 11 Effects of external explosion on internal pressure of the vessel

3 结 论

(1) 泄爆容器内的最大压力峰值和压力振荡取决于氮气体积分数：当φ≤40%时，由外部爆炸导致的内部压力峰值最大，而当φ＞40%时，则由泄爆膜破裂产生的压力峰值最大；在所有实验中，都观测到了亥姆霍兹振荡，且其频率随着φ的增加而降低；声学振荡只在φ=25%, 30% 时才出现，振荡频率约为1. 25 kHz。

(2) 内部最大爆炸超压随着氮气体积分数和压力监测点到泄爆口之间的距离而变化：对于给定的φ，pmax随着压力监测点与泄爆口之间距离的增加而增加，仅在φ=0%时，由于剧烈的外部爆炸，靠近泄爆口的pmax才略高于泄放容器中心的pmax。容器内部3 个测点间的压差随着φ的增加而减小，当φ≥40%时，压差几乎可以忽略不计。总体上，pmax随着φ从0%增加、到50%降低。

(3) 在一些实验中，外部压力曲线中出现了由外部爆炸产生的两个压力峰值。外部爆炸的强度随着φ的增加而降低。此外，外部爆炸显著地影响内部压力曲线，即使在外部爆炸超压很低的情况下，也能观测到由其导致的内部超压峰值。

（4）本文实验在小尺度的爆炸容器内开展，火焰图像和压力数据可为氢气爆炸抑制机理的研究提供基础数据。实验结果可为氢能源在工业生产、储运和利用过程中的防爆设计，开发和验证气体爆炸相关的CFD 软件（如FLACS）等，提供参考。