模拟立式拱顶油罐内油气爆炸实验研究*

蔡运雄，蒋新生，王世茂，余彬彬，王子拓，王春辉，李玉玺

（1. 陆军勤务学院油料系，重庆 401331；2. 中国空气动力研究与发展中心空天技术研究所，四川 绵阳 621000）

石油是重要的战略物资，由于其易燃易爆特性，在石油储备过程中容易发生火灾，由此引发的沸溢、喷溅现象不仅难以扑救，也对周围事物造成重大危害，发生在储罐区内的火灾还容易导致多个油罐连续爆炸[1-2]。2005 年，英国邦斯菲尔德油库因设备故障引发油气爆炸、燃烧，损毁储罐20 座，40 多人受伤，直接经济损失2.5 亿英镑[3-4]。2016 年12 月，广州石化公司的一个5 000 m3汽油储罐因油气泄漏引发闪爆，造成8 人伤亡。2020 年6 月，沈海高速温岭段一辆液化石油气运输槽罐车因侧翻导致重大爆炸事故，造成195 人伤亡，直接经济损失9 477 万元。2021 年5 月，河北沧州的5 个2 000 m3储油罐由于防火堤内的切割作业发生火灾，喷溅的火焰引起了相邻储罐的闪爆起火，直接经济损失3 872 万元。发生于储罐中的爆炸通常伴随着油气泄爆现象，冲击波与火焰的传播相对复杂，并且两者之间的相互作用促使爆炸发展更为激烈，由此产生巨大的破坏力。因此，探究立式拱顶油罐中油气爆炸的发展规律对相关领域的安全防护工程有着十分重要的意义。

当前，有关可燃气体爆炸的实验研究多采用氢气、甲烷和丙烷等单质气体，关注于气体体积分数和点火位置等初始条件对爆炸超压等特性参数的影响。例如，Rui 等[5]在边长为0.5 m 的立方体容器中开展了一系列氢气-空气混合物的泄爆实验，揭示了点火位置和泄爆压力对外场超压峰值与火焰速度的影响规律。Li 等[6]在城市管道人孔结构中利用甲烷/空气混合物，开展了不同气体浓度和点火位置等条件下的爆炸实验，研究表明当甲烷浓度越接近化学当量比并且点火位置越低时，泄爆到外场的冲击波超压峰值越大。苏航等[7]在微小空间中探究了丙烷-空气混合物以及在加氢情况下的火焰传播情况，实验观察到了光滑、皱褶和断裂三种火焰锋面形态。相比于这些单质气体，油气混合物的成分较为复杂，但仍然可以借鉴其研究方法来开展油气爆炸实验。然而，当前有关油气爆炸的实验通常都是在圆形、方形管道或者小型容器中开展的研究。例如，吴松林等[8]在圆柱形管道中开展了电火花、高温热壁等不同点火方式下的油气爆炸研究，指出点火方式对火焰结构、火焰颜色等都有较大影响。Li 等[9-10]在半开口的方形管道中探究了障碍物位置和油气体积分数对爆炸的影响，分析了超压变化、火焰传播以及超压与火焰之间的耦合关系。Wang 等[11]在20 L 圆形容器中研究了油气体积分数和泄爆口尺寸对泄爆过程的影响，指出超压变化和火焰演变存在4 个阶段，超压峰值、火焰速度与火焰传播距离受到初始条件的影响较大。此外，还有部分学者在不同结构的管道中开展了油气泄爆研究[12-14]，然而这些研究多是聚焦于容器内部的爆炸发展，较少关注到外场空间的爆炸演变过程，尤其是涉及到外场冲击波传播特性的相关研究较为欠缺。而相比于这些管道容器，油料储罐的结构较为特殊，当罐内发生油气爆炸时，通常会从侧壁与拱顶连接处的一点开始泄爆，然后逐渐扩大直到整个拱顶被掀开。虽然有部分学者针对这种结构形式的容器开展了少量实验研究[15-16]，但并没有关注到油气体积分数对火焰形态，以及点火位置和液位变化对油气泄爆带来的影响。然而在实际的油罐爆炸中，这些影响因素都在较大程度上决定了爆炸的威力，因而有必要对此开展全面的研究。

针对目前研究的不足，本文中，以立式拱顶油罐内发生油气爆炸对周围事物造成重大危害为工程背景，利用可视化的模拟油罐，系统地开展不同初始油气体积分数、不同点火位置和不同液位的油气爆炸实验，记录油罐内外场空间超压参数与火焰形态变化，探索立式拱顶油罐内油气爆炸的发展规律，以期为相关领域的安全防护工程提供理论支撑。

1 实 验

1.1 实验系统

实验系统如图1 所示，主要包括油气爆炸发生的容器，即透明的模拟立式拱顶油罐；构建实验工况的仪器设备，即油气产生与循环系统、红外线分析器、高能无干扰点火器等；探测实验参数的仪器设备，即动态数据采集系统、高速摄影仪、同步触发控制器等。其中爆炸容器是根据相似原理，以5 000 m3立式拱顶油罐为原型同比例缩制而成，总高380 mm，内径500 mm，拱顶曲率半径600 mm。罐壁采用28 mm 厚PMMA 材质，具有良好透光性。罐底采用10 mm 厚钢材，罐体设计承压2 MPa。模拟油罐拱顶采用白铁皮冲压成型，与罐壁的连接处进行强度弱化处理，以模拟拱顶的弱连接结构，设定损坏压力为20 kPa。实验中采用动态信号测试分析系统，配合使用0～10 kPa 和0～50 kPa 两种量程的压力传感器来采集压力信号，使用光电型火焰传感器来采集光照强度信号。其他设备与文献[10]相同，其中爆炸介质为92 汽油蒸气与空气的混合物，点火能量为1.5 J，高速相机拍摄帧率为1 000 s-1。

1.2 实验内容

针对空罐、25%储油量、50%储油量、75%储油量和100%储油量等5 种液位情况，分别开展了0.9%、1.1%、1.4%、1.6%、1.7%、1.8%、2.0%、2.3%和2.6%等9 种初始油气体积分数下的爆炸实验。对于空罐情况，在最危险油气体积分数下，开展了4 种不同点火位置的油气爆炸实验。当液位变化时，为保证安全，先在油罐下层倒入一定量的水，再在水层上倒入0.05 m 厚的汽油层，控制总液位达到设定值来模拟不同液位下的油气爆炸。实验发现，爆炸过程中罐内液体并没有大量波动，因密度差汽油层始终位于水层之上。而爆炸中主要是上层汽油挥发后参入后续的燃烧，0.05 m 厚的汽油层已经足够，采用水代替部分汽油对爆炸超压等特性参数带来的影响可以忽略。实验中利用动态数据采集系统和高速摄影仪记录爆炸超压和火焰形态的发展过程。如图2 所示，测点I1～I4 用于放置点火杆，分别为侧壁上部、中部、下部和罐底中心位置等4 种点火位置。测点F 用于放置火焰传感器，测点P1～P3 和P4～P27 分别为油罐内外场压力测点。I1～I3 和P1～P3 分别位于罐壁两侧的上中下3 个位置，间隔距离为0.2 m。24 个外场测点与油罐侧壁同高，均匀分布在6 条夹角60°的直线上，与油罐中心的距离分别为0.65、1.05、1.45 和1.85 m。实验时周围环境温度约为298 K，压强约为 98.5 kPa，对于每种工况至少重复开展3 次以保证实验的可重复性。

2 实验结果与讨论

2.1 初始油气体积分数对油气爆炸的影响规律

2.1.1 油气体积分数对超压参数的影响

油气爆炸超压是重要的爆炸特征参数，超压变化可以直观地反映出油气爆炸的发展过程，并衡量爆炸的威力。由于爆炸发展迅速，实验所用的爆炸容器较小，所以3 个内场压力测点所测得的超压数据十分接近。图3 所示为点火杆位于I4，初始油气体积分数φCH=1.8%时，3 个内场压力测点的压力时序曲线。从图中可以看出，3 条曲线的变化历程一致，P1、P2 和P3 处最大值分别为27.61、27.44 和27.55 kPa，误差不超过0.6%，达到最大值所用的时间都接近于91 ms。由于容器相对较小，爆炸并未发展到爆轰状态，传感器所测压力更多是由于热量释放导致容器内压力的整体提升，另外3 个传感器的位置比较接近，所以最大值和相应时间都比较接近，即油罐气体空间内不同位置的压力变化几乎一致[11]。所以在下文中，仅选取测点P2 的数据对油罐内场的超压变化进行分析。

根据内场超压时序曲线，油罐油气爆炸的发展大致可以分为3 个阶段，分别为超压上升阶段（0～91 ms）、超压泄放阶段（91～101 ms）和超压振荡衰减阶段（101～180 ms）。在第1 阶段，超压发展由慢而快，火焰对称发展。在91 ms 时，罐内超压到达拱顶承受极限，开始泄爆并进入第2 阶段。泄爆从拱顶与罐壁的连接处产生，且方向具有随机性。因为罐内冲击波经过多次反射在拱顶和侧壁连接的角隅处产生了压力汇聚，并且施加在拱顶的压力传递到连接处时形成了较大的拉伸力。而连接处的强度相对较弱，不同位置的强度并不完全一致，爆炸发生后强度最弱的地方会首先被破坏，然后沿着连接处向两侧扩展，进而使得整个拱顶被掀起。在第2 阶段，泄流效应使得罐内超压迅速降至负值，在101 ms 时产生最大负超压。泄爆过程中未燃油气被抛向外场空间，继续发生强烈的燃烧反应。在第3 阶段，超压保持一定的频率振荡衰减，外场空间的油气混合物进一步燃烧消耗。

由实验可知，不同油气体积分数下的超压变化趋势相同，都可以用3 个阶段来简单表述，仅在数值和时间上有所区别，表1 所示为不同初始油气体积分数φCH所对应的油气爆炸内场最大超压峰值pin,max、内场最大超压峰值的形成时间tin,max和内场平均升压速率(dp/dt)in,ave。可以看出随着φCH的提升，pin,max先增后减，在φCH= 1.7%时取得最大值31.59 kPa。其中内场最大超压pin,max实际上就是泄爆压力，实验中设定拱顶连接处破坏压力为20 kPa，是指在升压相对较慢时的泄爆压力。但在油气体积分数不同时，化学反应的激烈程度有所不同，升压速率差异明显。而拱顶连接处的破坏属于动压破坏，按照塑性变形滞后理论，在爆炸剧烈时，泄爆压力会更大。

表1 不同初始体积分数下的油气爆炸内场超压参数Table 1 Internal field overpressure parameters of gasoline-air mixture explosion at different initial volume fractions

平均升压速率反映了爆炸发展的快慢，内外场平均升压速率都可用于衡量爆炸的威力，是重要的爆炸特性参数。其中外场平均升压速率(dp/dt)ext,ave选用24 个外场压力测点中最大的超压峰值pext,max除以发展时间text,max。分析实验数据发现，外场超压大小与测点的位置关系密切，在拱顶泄爆方向的超压数值整体大于其他方向，距离最近测点的数值大于其他数值。所以pext,max实际上就是泄爆方向最近测点的最大超压峰值，而发展时间text,max是从泄爆开始到形成外场最大超压峰值所用的时间，(dp/dt)ext,ave也反映了冲击波在空气中传播的快慢。如图4 所示为(dp/dt)in,ave和(dp/dt)ext,ave与φCH的关系曲线。由图4 可知，两者的大致变化趋势都是随着φCH的提升先增后减，在φCH=1.7%时取得最大值，分别为0.46 和0.05 MPa/s，分别是φCH=0.9%时的4.6 倍和8.5 倍，是φCH=2.6%时的11.5 倍和15.8 倍。由此可见初始油气体积分数对爆炸发展的快慢影响很大，φCH=1.7%时的油气爆炸最危险。而(dp/dt)in,ave比(dp/dt)ext,ave大很多，这与冲击波在空气中的快速衰减有关。

2.1.2 油气体积分数对火焰传播的影响

不同初始油气体积分数下爆炸的火焰形态会呈现一定的差异性，但同样具有一定的规律性，图5所示为不同初始体积分数下油气爆炸火焰传播图像，选取0.9%、1.4%、1.7%、2.0%和2.6%等5 种体积分数下的火焰传播图像进行对比，其中初始油气体积分数为0.9%的火焰图像颜色较浅，对图片亮度进行了调整以便于观察。

油气爆炸时，在剧烈的化学反应过程中会有中间产物CH、C2、OH 等自由基生成，这些自由基使得火焰呈现出相应的颜色[17]，从图5 可以看出，在火焰发展的初始阶段，中心呈现橙红色，是因为在火焰中心没有外界氧气的补充，而且温度较高，靠近中心的碳氢分子受热产生了C2基，伴随有碳颗粒的析出，使火焰中心呈现橙红色[18]。而火焰周围呈现蓝色，是因为反应过程中产生了大量的CH、OH 等自由基，并且不断有外界氧气加入，使反应较为充分，其中CH 基使火焰呈现出蓝色。在火焰冲出罐外之后，油罐内部靠近壁面处仍有部分未来得及燃烧的油气混合物继续反应，所以罐内火焰在短时间里仍然呈现蓝色，在氧气被逐渐消耗之后，自罐底向上，火焰逐渐变为橙红色。火焰经由破口处泄放到外部空间后，外界氧气的加入以及泄爆扰动使得反应突然变得剧烈，放出大量的热，所以火焰在开口处呈现亮白色。

由图5 中可以看出，当φCH较小时，反应速率较慢但反应充分，即使在黑暗环境中，火焰形态也难以分辨，火焰颜色呈现出较浅的蓝紫色。是因为氧的体积分数较高，在反应过程中产生的CH 自由基迅速被消耗，使得火焰锋面内CH 自由基的体积分数较低，所以火焰颜色较浅。随着油气体积分数的提升，罐内火焰颜色由较浅的蓝紫色向较深的蓝色转变，进而转变为蓝绿色，再向黄色转变。这种变化主要是由于罐内油气体积分数的提升，使得氧含量变少，反应过程中形成的CH、C2等自由基来不及消耗产生累积，自由基的累积使得火焰表现出相应的颜色。当φCH过大时，碳颗粒的析出使得火焰直接呈现出黄色，此时罐内的反应属于贫氧燃烧，反应缓慢且不充分。另外，当φCH较大时，罐外火焰传播的更远。这是因为泄爆时，在浮力和冲击波的作用下，油罐顶部聚集有大量的未燃油气混合物，当拱顶被破坏时，这些被压缩且体积分数较高的混合物随着冲击波泄放到罐外空间，与空气混合继续燃烧。而φCH越大，泄放到罐外空间的油气越多，所以罐外火焰传播得更远。

表2 所示为不同油气体积分数下，储罐内火焰强度的最大峰值Imax与形成时间tmax，为避免换算中引入误差，其中火焰强度采用火焰传感器所采集的光照强度信号值代替，用于对比数值大小和定性分析。由表2 可知，随着φCH的提升，Imax先升再降，然后再升。在爆炸较为激烈的1.7%、1.8%体积分数附近Imax较低，是因为火焰传感器检测的是光的照度。结合图5 可以看出，φCH=1.4% 时，火焰较亮，随着φCH的提升，在燃烧过程中产生的C2基和碳颗粒逐渐增多，使得火焰变得深暗，所以Imax降低。但在φCH较高时，此时火焰内部碳颗粒更多，使得火焰颜色逐渐变为亮黄色，所以φCH=2.3%和φCH=2.6%时Imax有所提升。Imax的形成时间tmax随着φCH的提升先减后增，在φCH=1.6% 时取得最小值184 ms，φCH=1.7%时的tmax与之相近。另外，出现火焰强度最大的时刻并不是在超压上升最快的罐内发展阶段，而是在泄爆产生之后，是因为爆炸前期的燃烧反应较为充分，各种自由基消耗较快，火焰颜色较淡。而在泄爆后，火焰锋面受到一定的扰动作用，锋面处反应速率加快，热量释放增多，因而火焰亮度变高，所以火焰强度最大的时间出现在泄爆后。

表2 不同初始体积分数下油气爆炸罐内火焰强度最大峰值和形成时间Table 2 Maximum peak flame intensity and time of formation in the tank at different initial volume fractions

2.2 点火位置对油气爆炸的影响规律

在工程实际中，油罐内出现电火花等点火源的位置并不明确。为了探究点火位置变化对油气爆炸的影响规律，在空罐情况下，开展了侧壁上部、中部、下部和罐底中心位置等4 种点火位置的油气爆炸实验，所用油气体积分数为上文中得出最危险体积分数，即φCH=1.7%。如表3 所示为相关爆炸特性超压参数，其中pin,max、(dp/dt)in,ave和pext,max都随着点火位置的下移而呈现增大趋势，在罐底中间点火时取得最大值，分别为31.59 kPa、0.47 MPa/s 和0.53 kPa，相比于侧壁上部点火分别增长了44.2%、38.2%和194.4%。由此可见，点火位置对油罐内场超压影响很大。其原因是点火位置靠近上面时，冲击波能够很快抵达弱连接处，及时泄爆使得超压上升阶段持续时间较短，密闭空间的爆炸转换为开敞空间的爆炸，爆炸威力无法因憋压而得到增强。对于外场超压参数，当点火位置靠近上部时，text,max较小，而冲击波相对较弱。当点火位置靠近下部时，text,max较大，而冲击波相对较强。所以(dp/dt)ext,ave并不是随着点火位置的下移单调变化，但同样在罐底中间点火时取得最大值0.05 MPa/s，相比于上中下3 个点火位置分别提升了25%、66.7%和66.7%。以pin,max、pext,max、(dp/dt)in,ave和(dp/dt)ext,ave作为衡量指标，油罐底部中心位置点火的油气爆炸威力最大。所以在实际应用中，当罐内油品用完之后，应及时排空油气，可采用惰气置换或者其他方式进行，同时避免罐底产生点火源。

表3 不同点火位置的油气爆炸超压参数Table 3 Overpressure parameters of gasoline-air mixture explosion at different ignition locations

图6 所示为不同点火位置的油气爆炸罐内火焰传播图像，从每种工况的最后一个火焰图像可以看出，点火位置越靠下，火焰泄放到罐外时罐内火焰锋面的面积越大，即泄爆时参入反应的油气混合物越多。结合图3 可以看出，在超压上升阶段，随着时间的增加，罐内超压的增长速度越来越快，也即平均泄压速率越来越大，这些从侧面解释了为什么pin,max和(dp/dt)in,ave均随点火位置的下移而增大。另外当点火位置为罐底中心时，拱顶连接处的破坏方向具有不可预测性，但当点火位置在侧壁时，拱顶一般从点火位置上方的连接处开始破坏，例如图6(a)中的68 ms、图6(b)中的70 ms 和图6(c)中的72 ms 时刻的火焰图像，连接处都以这种方式破坏。所以在工程实际中，更容易产生点火源方向的油罐安全距离应该更大。

2.3 液位对油气爆炸的影响规律

2.3.1 液位对最危险油气体积分数的影响

在实际应用中，油罐内液位不断变化，当罐内储存油料时，油气爆炸过程中会涉及到液面的变形、液滴的破碎与蒸发等较为复杂的相变过程，对爆炸的发展产生影响。另外，液位变化引起油罐上方油气混合物的空间发生变化，相当于爆炸容器发生改变。根据课题组的前期研究成果可知，受限空间的容积和形状变化都对油气爆炸的发展造成影响[19]。为探索液位不同对油气爆炸发展和爆炸威力的影响规律，对空罐、25%储油量、50%储油量、75%储油量和100%储油量共5 种情况，分别开展了9 种油气体积分数下的爆炸实验，分析相关参数，为控制无关变量，点火位置均选用图2(a)中I1 点火位置，压力参数选用P1 测点数据。

当液位较低时，罐内气体空间相对于空罐来说变化不大，因而选取75%和100%两种液位下的数据进行分析，其中100%液位是指到达油罐安全高度的液位，本文以罐壁高度的0.9 倍作为安全高度，在该液位时，油罐上部仍然具有较大的气体空间。表4 所示为不同液位下不同φCH的爆炸超压特性参数，由表中数据得到pin,max、pext,max与φCH的关系，分别如图7～8 所示。

表4 不同液位下不同初始体积分数油气爆炸内外场最大超压峰值Table 4 The maximum overpressure peak in the internal and external field of gasoline-air mixture explosion with different initial volume fraction at different liquid levels

由图7～8 可以看出，无论是75%液位还是100% 液位，随着φCH的增大，pin,max和pext,max都呈现出先增后减的趋势，在个别体积分数上有些波动，但不影响总体趋势，并且都在φCH=1.7%时取得最大值，分别为24.02 和0.22 kPa，22.47 和0.18 kPa。在75%液位下，pin,max和pext,max在φCH=1.7% 时，相比于φCH=0.9% 分别增大了21.9% 和100%，相比于φCH=2.6% 分别增大了20.4%和69.2%。在100%液位下，相比于φCH=0.9% 分别增大了13% 和260%，相比于φCH=2.6%分别增大了11.5%和500%，可见液位变化对pext,max的影响比pin,max要大。通过对其他液位下的数据进行分析，发现其超压数值表现出同样的规律。根据前面的分析可知，超压峰值和平均升压速率基本都在同一工况下出现。所以以pin,max和pext,max作为衡量指标，φCH=1.7%是油罐在任一液位下的最危险体积分数，即此体积分数下的油气爆炸威力最大，液位变化对最危险体积分数的取值影响不大。

2.3.2 液位对超压参数的影响

由表4 可以看出，对应于某一种油气体积分数，不同液位下油气爆炸的超压特性参数具有差异性，例如1.7%体积分数下，75%和100%液位的内场最大超压峰值pin,max分别为24.02 和22.47 kPa，说明在同一种体积分数下，液位变化对油罐内场超压具有影响。为了解其影响规律，选取φCH=1.7%时不同液位下油气爆炸参数进行对比，具体数据见表5。

表5 不同液位油气爆炸内场超压参数Table 5 Overpressure parameters inside tank of gasoline-air mixture explosion at different liquid levels

由表中数据可知，随着液位从0 到100%，tin,max逐渐减小，而pin,max和(dp/dt)in,ave都呈现先增后减的趋势，在50%液位时取得最大值，为25.40 kpa 和0.43 MPa/s。最大值并不是在空罐时产生，是因为点火位置在I1，发生油气爆炸时，拱顶是受到爆炸产生的入射冲击波和罐壁、罐底或液面反射的冲击波共同作用下才发生失稳的，当液位从0 升到50%时，液位的提升使得气体空间更小，冲击波的反射更加迅速，所以在拱顶失稳前的极短时间内，压力上升更快，平均升压速率的变化也一样。因而相比于0 和25%液位，50%液位下的pin,max分别提升了15.9%和9%，(dp/dt)in,ave分别提升了26.5%和16.2%。但是随着液位继续提升，从50%到100%，最大超压峰值逐渐降低，是因为油罐上方气体空间的减少使得油气爆炸的发展受限，在油气爆炸发展的初始阶段，受火焰不稳定性和压力冲击波的耦合作用，化学反应速率是逐渐加快的。但是当油罐上方气体空间太小时，油气爆炸的发展便受到体积限制，使得拱顶过早被破坏，泄压时间更早。因而相比于50% 液位，75% 和100% 液位的pin,max分别降低了5.4%和11.5%。

液位变化引起油罐上方气体空间体积产生变化，使得爆炸特性参数发生相应变化。为探究两者之间的对应关系首先需要计算不同液位下的气体空间体积，即油气混合物的体积V。根据模拟油罐的尺寸计算可得，0 至100%不同液位罐内油气混合物的体积分别为0.069、0.055、0.041、0.026 和0.012 m3，如图9 所示为φCH=1.7%时，pin,max、pext,max与V的关系。

由图可知，随着油气混合物体积V的增大，pin,max和pext,max均先增后减，在0.041 m3，即50%液位时取得最大值，分别为25.40 和0.27 kPa。对于泄压爆炸，当点火位置接近泄压口时，内场的爆炸发展受到空间体积的直接影响，最大超压峰值受到空间形状、可燃气体化学性质等因素的控制。当罐外没有发生猛烈的油料蒸汽云爆炸时，外场超压的产生主要是内场爆炸冲击波泄放到罐外空间造成的，内场超压越大外场超压也越大，因而两者呈现的变化趋势相同。以pin,max和pext,max作为衡量指标，50%液位是油罐的最危险液位，即此液位下的油气爆炸威力最大，而且在实际情况下，更可能伴随有油品的沸溅，引燃周围物体，进一步造成损害。但是对比表3 可知，当罐内没有油品时，若点火位置为油罐底部中心，内外场最大超压峰值分别为31.59 和0.53 kPa，即此时油气爆炸产生的毁伤力更大，造成两者不同的主要原因是点火位置的不同，在实际工程中，对于这两种情况都应该给予足够的重视。

2.3.3 液位对外场超压衰减的影响

探究外场冲击波的衰减对于认识储罐油气爆炸的威力十分重要，外场超压不但与油罐内部油气体积有关，还与测点位置有关。如图10 所示为空罐情况下，φCH=1.7%时外场P4～P27 压力测点的最大超压峰值分布，其中泄爆方向设为0°。由图可知，泄爆方向的超压数值大于其他角度的数值，与泄爆方向的夹角越大，数值通常越小。为了反映爆炸威力的传播距离，选取泄爆方向的超压数据分析不同液位外场冲击波的衰减规律，因此引入比例距离R:

式中：L为外场测点与模拟油罐中心的水平距离，V是模拟油罐中混合气体的体积。计算可得不同液位时泄爆方向上各测点的比例距离，空罐时，4 个测点的比例距离分别为1.58、2.56、3.53 和4.50；25%液位时，比例距离分别为1.71、2.76、3.81 和4.86；50%液位时，比例距离分别为1.89、3.05、4.22 和5.38；75%液位时，比例距离分别为2.19、3.53、4.88 和6.23；100%液位时，比例距离分别为2.85、4.61、6.36 和8.12，不同液位油气爆炸外场某一距离下的最大超压峰值与比例距离的关系如图11 所示。

由图可知，不同液位时外场超压的衰减规律相似，可以采用幂指数公式进行拟合[20,21]。

式中：pext表示油罐外场某一距离下的最大超压峰值，k为指前因子， α 为衰减系数。

式(3)～(7)分别为0、25%、50%、75%和100%液位的拟合公式，相关系数分别为0.990、0.976、0.976、0.995 和0.973。5 个公式的指前因子分别为0.32、0.37、0.70、0.71 和0.98，衰减系数分别为1.29、1.33、1.47、1.51 和1.65，即随着液位的升高，油气混合物体积的减小，两个参数都呈现增大趋势，两者与油气混合物体积V的对应关系都可以用二次多项式拟合：

拟合公式(8)和(9)的相关系数分别为0.938 和0.978。所以，在不同液位情况下油气爆炸外场冲击波超压峰值与距离和油气混合物体积的关系可以表示如下：

2.3.4 液下超压分析

通过对内场不同测点超压数据的分析发现，当存在液位时，放置于液面以下压力传感器测得的数值与液面以上压力传感器测得的数值存在一定的差异，接下来以50%液位时的油气爆炸内场超压数据进行分析。选用图2(a)所示P1 和P3 两个测点的数据，其中P1 所测数据为气体空间中的爆炸超压，P3 所测数据为液体空间中的爆炸超压。由于测点P3 位于液位以下，相比测点P1 多了液体施加的压强。为了便于对比分析，将测点P3 的数据减去液体压强，得到如图12 所示的不同测点的超压时序曲线，其中φCH=1.7%，点火杆位于油罐侧壁上部位置。

根据前文，将油气爆炸内场压力发展分为3 个阶段：超压上升阶段、超压泄放阶段和振荡衰减阶段。由图中曲线可知，在超压上升阶段，液位上、下2 个测点的超压数据都呈现上升趋势，但是测点P1 的最大超压峰值要大于测点P3，两者分别为25.40 和24.53 kPa，相差0.87 kPa，约为测点P1 超压峰值的3.4%。而测点P1 和测点P3 到达最大超压峰值的时刻分别为58.8 和62.1 ms，平均升压速率分别为0.43 和0.40 MPa/s，前者达到峰值的时间稍早，且平均升压速率稍大于后者。在超压泄放阶段，两个测点都在短时间内降为负值，但是测点P1 形成最大负超压的时间早于测点P3，而且两者的最大负压数值相差较大，分别为-3.51 和-15.11 kPa，测点P3 的最大负超压约为测点P1 的4.3 倍。在振荡衰减阶段，测点P1 的超压数据在第一个振荡周期里比较杂乱，这与冲击波在液面的反射和拱顶泄爆的扰动有关，之后超压缓慢振荡衰减，逐渐回到零点。而测点P3 的超压振幅和振荡频率明显大于前者，经过一段时间的振荡衰减之后回到零点，两者的振荡频率分别约为42.6 和53.8 Hz。

当压力冲击波在不同介质中传播时，介质物理属性对冲击波的影响较大[22]，当冲击波传播到介质交界面时，会产生扰动，影响冲击波的传播，并且不同介质内冲击波的传播存在明显差异[23]。发生于储罐内的油气爆炸属于爆燃，在初始阶段受到外界的扰动较少，化学反应速率较慢，火焰锋面要落后于前驱冲击波阵面，形成如图13 所示的双波三区结构。前驱冲击波的形成是化学反应区能量释放的一种方式，前驱冲击波和火焰锋面分别以速度Ds和Df传播，将罐内流场划分为三个区域，其中0 区的可燃气体未被扰动，处于初始状态，1 区的可燃气体分子经过冲击波的扰动后参与到化学反应中去，2 区为燃烧后的区域，即燃烧产物的聚集区。压力波自爆炸中心向四周传播，到达液面时产生反射波和透射波。由镜像爆心原理[24]，反射冲击波可看作是由与爆心关于液面对称的镜像爆心所产生的冲击波，但由于反射削弱作用，其强度要小于入射波。由图13 可知，点火位置I1 与测点P3 的距离比测点P1 远，更主要的是液相空间的压力是由气相空间传递而来的，所以在压力响应上测点P3 相比于测点P1 具有延后性，但由于模拟油罐的体积较小，两者时间间隔较短。冲击波在液面产生的反射，削弱了传入液体空间的透射冲击波强度，而且反射冲击波使得气体空间的压强更高。由于液体的黏性比空气大很多，阻尼系数更高，在冲击波传至液面以下超压测点时，已经消耗了部分能量，所以在最大超压峰值上，气体空间的数值大于液体空间。由相关研究可知，冲击波负压的形成与稀疏波和介质的惯性有直接关系[25]，当冲击波到达油罐壁面后会反射得到稀疏波，在介质惯性的共同作用下，该区域空间处于拉伸状态，即产生负压。由于气体和液体的密度、弹性模量等物质属性相差很大，所以在超压泄放阶段，液位以下测点P3 的最大负超压明显低于液位以上测点P1 的数值，负压的强度更大。在随后的超压振荡衰减阶段，由于液体的阻尼系数更高，测点P3 的超压衰减更加迅速，而且观察可知，液位以下超压的振荡频率更高，这是因为在液体中冲击波的传播速度更快。液相空间的超压变化在振荡衰减阶段振幅较大、频率较快，所以当罐内储存油品而发生爆炸时，罐底板在爆炸振荡中容易受损。

总的来说，油气爆炸时，液位上下不同测点超压变化的总体趋势相近，但由于冲击波在气相和液相之间传播的差异性，液相空间测点的超压响应要慢于气相空间测点，在超压峰值的形成时间上具有延后性，而且气相空间的最大超压峰值稍大于液相空间的最大超压峰值，液相空间的最大负超压明显大于气相空间的最大负超压，在超压衰减阶段，液相空间内的超压振荡幅值和频率大于气相空间，而且衰减更快。在不同液位、不同初始油气体积分数或者不同点火位置时，液位上下测点的具体数值在大小上有所差别，但是两者之间的对应规律与之相似。

2.3.5 液位对火焰传播的影响

当罐内液位发生变化时，液面对火焰的传播产生影响，但液面的存在并不是简单地阻碍火焰锋面的传播，还涉及到热量的吸收、液体的相变等一系列复杂的相互作用，这使得存在液位时的火焰发展与空罐情况下的油气爆炸有所不同。为便于观察，选取25%和50%两种液位下的油气爆炸火焰传播图像进行分析，其中拱顶破坏前的火焰传播与空罐时类似，主要展示破口时和破口后的火焰传播，如图14 所示，初始油气体积分数都为1.7%。

由图可知，当液位较高时，火焰从破口处泄放的时间更早，这与表5 中液位越高到达最大超压峰值的时间越短是相对应的，这也显示了压力与火焰变化的同步性。如图14(a)中80 和84 ms 时刻的图像，在火焰泄放到外场空间时，有部分液体的相变参与反应，这些使火焰传播与超压变化变得更为复杂。在25%液位时，部分火焰传播到罐外空间后，罐内火焰以四分之一球形向前传播，如图14(a)中156 ms 时刻。当右侧火焰前端传到罐壁时，受拱顶和罐壁的共同限制，火焰前端开始向下转折，如图14(a) 中196 ms 时刻，直到受到液面的限制而停止。在50%液位时，罐内火焰的传播受到液面的限制更大，如图14(b)中76 ms 时刻，当火焰前锋即将接触液面时，受到马赫杆的作用[26]，罐内流场与罐壁平行，罐内火焰锋面逐渐变平。到128 ms 时，罐内火焰锋面已经变得平整，之后，受到拱顶反射冲击波的影响，右侧火焰锋面逐渐向下传播，直到受到液面限制而停止。从火焰锋面随时间变化的图像可以看出，罐内火焰传播时间越久，锋面越褶皱，越趋于形成胞格结构[27-28]，其原因是火焰固有的不稳定性以及冲击波对火焰前锋的扰动[29-30]。胞格结构的形成增大了锋面面积，使得燃烧速率大幅提升，让爆炸发展更为激烈。而更加激烈的爆炸会形成更强的冲击波和反射冲击波，在传播过程中继续扰动火焰的传播。当其穿过火焰锋面时，会在锋面内产生涡量，增大该区域的湍流度，增加活性自由基团相互碰撞的概率，从而提升化学反应速率，使爆炸发展更为激烈。而在泄爆时，冲击波的扰动作用还会促进未燃油气与新鲜空气的混合，进一步加快爆炸的发展，火焰与冲击波的相互作用、相互耦合机制贯穿整个油气爆炸的发展过程。

3 结 论

本文中利用可视化实验台架，开展了9 种初始油气体积分数、4 种点火位置和5 种液位下的油气爆炸研究，探究了立式拱顶油罐内油气爆炸的发展规律，主要结论如下。

(1)不同工况下的油气爆炸内场超压发展都可以分为超压上升、超压泄放和振荡衰减3 个阶段。空罐条件下，内外场平均升压速率都随着油气体积分数的提升先增后减，在1.7%时取得最大值。不同初始油气体积分数和不同阶段下的油气爆炸呈现出不同的火焰颜色变化，这与化学反应中间产物CH、C2、OH 等自由基的生成和空间分布有关。随着初始油气体积分数的提升，罐内火焰颜色由较浅的蓝紫色向较深的蓝色转变，进而转变为蓝绿色，再向黄色转变，火焰的光照强度先上升再下降，然后再上升，在油气体积分数较高并且在泄爆之后取得最大值。

(2)点火位置对油气爆炸的影响较大，点火位置越靠下，泄爆越晚，内外场最大超压峰值越大。罐底中心点火时，内外场平均升压速率取得最大值，分别为0.46 和0.05 MPa/s。空罐条件下，油罐底部中心位置点火时的油气爆炸威力最大。