城市地下浅埋管沟可燃气体爆炸的灾害效应（Ⅱ）：影响因素分析及后果评估*

2023-03-02杨石刚蔡炯炜孙文盛门敬敏

爆炸与冲击 2023年1期

杨石刚，蔡炯炜,2，杨亚，孙文盛，门敬敏

（1.陆军工程大学爆炸冲击防灾减灾国家重点实验室，江苏南京 210007；2.73021 部队，浙江杭州 310012）

近年来，随着产业升级的不断深入，我国对天然气的需求量在稳步提高，天然气在我国能源系统中的主体地位也得到了稳固。然而天然气作为一种可燃气体，无论是工业生产还是日常使用极易发生泄漏，当遇到合适的点火源时，会发生严重的爆炸事故。如2021 年6 月，湖北十堰的一条煤气管道被腐蚀，导致煤气泄漏发生爆炸，事故造成了25 人死亡，138 人受伤。2020 年6 月，浙江温岭的一辆油罐车阀门遭到严重的侧撞，导致可燃气体泄漏发生爆炸，事故造成了20 人死亡，172 人受伤。而当泄漏的可燃气体进入到城市地下浅埋管沟中积聚爆炸时，会造成更严重的灾害后果，如2013 年11 月山东青岛的排水暗渠可燃气体爆炸，2014 年7 月台湾高雄的地下管沟可燃气体爆炸。由于爆炸事故会造成严重的人员伤亡和财产损伤，也会给社会带来极其恶劣的影响，因此对城市地下浅埋管沟可燃气体爆炸灾害后果进行分析和评估具有重要意义。

目前，对于爆炸灾害的后果评估，运用最多的是超压准则，超压准则认为超压值的大小是评估人员和建筑物是否受到伤害和破坏的唯一评判标准，当冲击波超压超过某一阈值时，会对目标造成一定的伤害[1-4]。而在某些情况下，超压峰值往往并不是很大，但是由于超压作用时间较长，也会造成较严重的影响，目前美国国防部在统一设施标准[5]中根据冲量准则给出了人员的伤害标准。超压-冲量准则综合考虑了目标性质、破坏等级、冲击波超压和冲量这4 个参数作为参考值，认为对目标的灾害效应评估应该综合考虑超压和冲量，只有当两者的组合达到某一临界值时才会对目标造成相应的伤害，不少学者也得到了一些研究成果[6-11]。此外，李峰[12]结合了超压准则和超压-冲量准则，首次提出了人体伤害当量的概念，用于预测地下交通爆炸对人员的伤害，并针对隧道内发生爆炸的场景给出了人员伤亡的分布区域。也有一些学者[13-15]利用概率方程，即根据不同的伤害类型，利用超压和冲量的组合得到伤害百分比，实现了对冲击波灾害效应的定量评估，以此可以确定冲击波对目标的伤害等级及危险距离。另外，Assael 等[2]列出了人员的鼓膜破裂、肺出血、头部及全身撞击死亡的概率方程，Aloson 等[16]则将概率方程应用到特征曲线中。目前，关于城市地下浅埋管沟可燃气体爆炸的灾害后果评估尚未见有公开报道，现有的灾害后果评估研究成果主要是借鉴固体炸药爆炸事故，对城市地下浅埋管沟可燃气体爆炸的适用性不强，功能性也不完善，需要进行更具体和深入地研究。

本文中，利用计算流体动力学软件FLACS，对城市地下浅埋管沟可燃气体爆炸荷载进行数值模拟，选用超压准则评估传播到地面的爆炸冲击波对建筑物破坏和人员伤害的危险距离，并对影响荷载分布和危险距离的因素进行分析，以期为建筑物的安全设计以及爆炸事故的灾害预防提供一定参考。

1 管沟可燃气体爆炸荷载数值模拟

1.1 模拟工况设置

管沟的数值模型如图1 所示，其中L 为泄爆口之间的距离，l 为气云的长度，D 为正方形泄爆口的边长。分别考虑点火点位置、泄爆口大小、气云长度和管沟横截面面积这4 个因素的影响，工况记录如表1 所示。工况1～4 分别表示点火点位于管沟内部(4/8)L、(5/8)L、(6/8)L、(7/8)L 处，用以研究点火位置对危险距离的影响。在实际管沟中，雨水井是日常检修预留的孔口，将工况1、5 和6 的泄爆口边长D 分别设置为 1.0、0.8 和0.6 m，研究泄爆口大小对危险距离的影响。气体爆炸的强度与气体量的大小紧密相关，工况1、7～10 分别表示气云长度l 为90、60、40、20 和10 m 的可燃气体爆炸，研究气云长度对危险距离的影响。参照相关规范[17]，将工况1、11 和12 分别表示管沟横截面面积为1、2 和3 m2的可燃气体爆炸，研究管沟横截面面积对危险距离的影响。

图1 管沟的数值模型Fig.1 The numerical model established for a pipe trench

表1 管沟可燃气体爆炸数值模拟工况Table1 Numerical simulation conditions on combustible gas explosion in a pipe trench

1.2 影响因素分析

由已有研究结果可知，管沟内可燃气体爆炸冲击波通过泄爆口传播到地面的过程可分为稳定段、超压峰值Δp1段、超压峰值Δp2段[18]。由于超压峰值Δp1较小，对地面的影响较有限，而超压峰值Δp2大、危险性高，因此本节仅选取超压峰值Δp2作为研究内容，在空气域中沿X、Y、Z 方向上布置的测点作为研究对象，用于记录超压峰值Δp2的变化，分析其影响因素。

1.2.1 点火位置的影响

图2(a)和(b)分别为在不同的点火位置点火时，在X 和Y 方向上的超压峰值分布，其中红色虚线表示泄爆口所在位置。可以发现，当点火点位于(4/8)L 处时超压峰值最大，而当点火点位于(7/8)L 处时超压峰值最小；在泄爆口两侧沿Y 方向上超压峰值大致呈对称衰减的趋势，但是在泄爆口两侧X 方向上超压峰值呈不对称衰减的特征。

图2 不同点火位置时的超压峰值分布Fig.2 Peak overpressure distribution at different ignition positions

在Z 方向上，超压峰值随高度的上升而逐渐减小（见图2(c)）。当Z＜5 m 时，工况1 中超压峰值最大；当Z＞5 m 时，工况2 中超压峰值最大，这是由于空气流场的带动在泄爆口上方有大量的残余燃料，当火焰锋面传播到泄爆口上方时，残余燃料与火焰波阵面接触后会继续燃烧，此时部分火焰已传递到泄爆口上方约5 m 处，导致此处的超压峰值衰减较慢。而工况1 中，燃料受空气流场的影响，被大量排到泄爆口上方区域，而此时火焰波阵面仅传播到管沟内部，因此泄爆口上方的残余燃料并没有参与燃烧而是仅被气流推动到空气域中，在泄爆口上方测点的超压峰值呈现快速衰减的趋势。因此在Z=5 m 的分界点处，工况2 的超压峰值出现反超的情况。而工况4 中，在Z 方向上超压峰值始终最小，可能的原因是：当处于该点火位置时，燃料参与燃烧反应的量较小，导致火焰波传播到泄爆口处的能量也较低，超压峰值偏小。

1.2.2 泄爆口大小的影响

图3 为在不同的泄爆口面积下，测点的超压峰值分布情况。在X 负方向上，各工况之间的超压峰值相差较小，不超过10%；而在X 正方向上，当D=0.8 m 时，测点的超压峰值最大。在Y 方向上，当D=0.8 m 时，泄爆口处的超压峰值达到7.8 kPa，分别比D=0.6 m 和D=1.0 m 时的大12.2%和29.3%。这是由于X、Y 方向的超压峰值主要受内部爆炸荷载大小、爆炸波从泄爆口传出角度2 个因素的影响。随着泄爆口尺寸的增大，在泄爆作用下，一方面内部爆炸荷载随之减小，两者呈负反馈机制，另一方面爆炸波传播方向更偏向地面，对于X、Y 方向的超压是正反馈机制。所以在两者的共同作用下，X、Y 方向的超压峰值与泄爆口大小之间并不呈规律性变化。而在Z 方向上，D=0.6 m 时的超压峰值最大，这是由于在Z 方向的超压峰值主要受内部爆炸荷载的影响，随着泄爆口尺寸的减小，对应的超压峰值增大，同时受泄爆口的约束作用，其向空气域传播的方向也更偏向于沿Z 方向传播。

图3 不同泄爆口面积时的超压峰值分布Fig.3 Peak overpressure distribution under different vent areas

1.2.3 气云长度的影响

由图4 可知，测点的超压峰值受管沟内气云长度l 的影响。当气云的长度为90 m 时超压峰值最大，而气云长度为10 m 时超压峰值最小。但同时也发现当气云长度为60 和40 m 时，超压峰值与气云长度为90 m 时的相差不大。这说明，当气云达到一定长度后，继续增大气云长度对超压的影响不大。可能的原因是，管沟内参与燃烧反应的气云量有限，未能参与反应的气云只能通过泄爆口被排放到空气域中稀释，而并未参与到燃烧反应中，因此对超压峰值的影响有限。

图4 不同气云长度下的超压峰值分布Fig.4 Peak overpressure distribution under different gas cloud lengths

1.2.4 截面面积的影响

由图5 可知，当管沟横截面面积S 增大2 倍时，X、Y、Z 方向上的最大超压峰值分别增大到3.3、4.6 和4.8倍；而当管沟横截面面积增大3 倍时，X、Y、Z 方向上的最大超压峰值分别增大到4.8、7.7 和6.9 倍。这是由于当管沟横截面面积增大时，管沟内参与燃烧反应的气云量也增加，加剧了燃气爆炸的反应程度，使得测点处的超压峰值也变大。管沟横截面面积变化时，Y、Z 方向的超压峰值变化更为敏感，这是由于冲击波在管沟内沿X 方向有限距离传播，而传播到泄爆口后进入到无限空气域中，管沟横截面积的影响会相对弱化。

图5 不同横截面面积的管沟内燃气爆炸超压峰值分布Fig.5 Peak distribution of gas explosion overpressure in trenches with different cross-sectional areas

2 管沟可燃气体爆炸灾害效应评估

选用超压准则，来评估管沟内的爆炸冲击波传播到地面后对人员伤害和建筑物破坏的危险距离。冲击波超压对建筑物的影响可以参照化工安全的行业标准[19]，如表2 所示。本文中分别以2.07、6.90、34.50 kPa 作为爆炸冲击波超压造成建筑物轻度破坏、中度破坏和重度破坏的临界值。

表2 冲击波超压对建筑物的影响[19]Table2 Impact of shock wave overpressure on buildings[19]

根据爆炸造成人员伤亡的不同概率，可以将爆炸危险源的中心从内向外依次可划分为死亡区、重伤区、轻伤区、安全区[8]，表3 为冲击波超压对人员伤害的区域划分[3]。本文中分别以10、45、75 kPa 作为爆炸冲击波超压造成人员轻伤、重伤和死亡的临界值。

表3 冲击波超压对人员的影响[3]Table3 Impact of shock wave overpressure on personnel[3]

2.1 冲击波对建筑物破坏的危险距离

图6 为当Y=0 时，沿XOZ 面爆炸冲击波对建筑物破坏区域的剖面图，蓝色表示轻度破坏区域，红色表示中度破坏区域，绿色表示重度破坏区域，但由于爆炸冲击波的能量有限，在地面上并没有出现重度破坏的区域。由图可知，轻度、中度破坏区域的空间形状并不规则，且轻度破坏区域比中度破坏区域要大得多。

图6 可燃气体爆炸对建筑物破坏区域的剖面图Fig.6 Sectional view of the damage area of the building caused by the explosion of combustible gas

为了能定量分析爆炸冲击波对建筑物造成破坏的最大危险距离，汇总了不同高度处建筑物破坏区域的二维分布图（见图7）并得到最大的危险距离，图中蓝色为轻度破坏区域，红色为中度破坏区域。中度及轻度破坏区域与高度Z 有关，当Z＞7 m 时，随着爆炸冲击波的衰减，并没有出现中度破坏（见图7(a)），随后将其各自的等值线投影到XOY 平面（见图7(b)）。

图7 可燃气体爆炸对建筑物破坏区域的二维分布图Fig.7 Two-dimensional distribution of the damage area of the building caused by the explosion of combustible gas

当Z=3 m 时，中度破坏区域的范围最大，因此为了得到该最大危险距离，画出该高度下中度破坏的等压力伤害线（见图8)。可以得到等压力伤害线上的点离泄爆口的距离互不相等，由于其沿Y 轴对称分布，因此当Y=0 时，等压力伤害线与X 轴正方向的交点的距离为3.4 m，为中度破坏的危险距离。由于在确定安全范围时，需把不确定因素考虑在内，因此危险范围可以扩大为以泄爆口为圆心、半径为3.4 m 的圆形区域。同理，当Z=3.5 m 时，轻度破坏区域的范围达到最大，该危险距离为7.5 m。

图8 中度破坏危险范围Fig.8 Dangerous range of moderate damage

图9(a)为点火位置与危险距离的关系，橙色部分为轻度破坏区域，绿色部分为中度破坏区域，红色部分为重度破坏区域，其余部分为安全区域。由图可知，对于任一点火位置，轻度破坏的危险距离比中度和重度破坏要大。当点火点位于管沟(4/8)L 处时，中度破坏和轻度破坏的危险距离最大，分别为7.0 和2.8 m。而当点火点位置位于管沟(7/8)L 处时，中度破坏和轻度破坏的危险距离为零。

图9(b)为管沟泄爆口大小与危险距离的关系。当泄爆口边长在0.6～1.0 m 的范围时，中度破坏和轻度破坏的危险距离分别在3～4、5～9 m 的范围之间，幅度波动不大。其中当泄爆口的边长为0.8 m 时，中度破坏及轻度破坏的危险距离最大，分别为8.6 和3.9 m；而当泄爆口边长为0.6 m 时，中度破坏及轻度破坏的危险距离最小，分别为5.2 和3.1 m。

图9(c)为气云长度对危险距离的影响。由图可知，气云的长度越长，危险距离就越大；但当气云达到一定的长度后，危险距离基本保持不变。

图9(d)为管沟截面面积对危险距离的影响。从图中可得，当管沟的截面面积越大时，轻度破坏和中度破坏的危险距离也越大。当截面面积为2 m2时，建筑物出现了重度破坏，其危险距离为1.5 m，而当管沟截面面积为3 m2时，重度破坏危险距离达到2.4 m。

图9 不同因素对建筑物破坏危险距离的影响Fig.9 Effects of different factors on the dangerous distance for building damage

2.2 冲击波对人员伤害的危险距离

汇总不同高度处人员危险距离的二维剖面图并得到最大危险距离，如图10 所示，图中红色为人员轻伤区域。可以得到在地面上并没有出现人员的重伤区和死亡区。

图10 人员轻伤区的二维剖面图Fig.10 Two-dimensional section view of lightly-injured area

由于在不同高度下的人员轻伤区具有不同的危险区域，因此可以得到其所对应的等压力伤害线，如图11 所示。当Z=2.9 m 时，危险距离达到最大值3.2 m。

图11 人员轻伤区危险距离Fig.11 Dangerous range of minor injury zone

图12 为点火位置、泄爆口大小、气云长度和管沟截面面积与人员危险距离的关系，橙色部分为轻伤区，红色部分为重伤区，其余为安全区。由图可知，当点火点位置位于管沟(4/8)L 时，人员受轻伤的危险距离最大；而当点火点位置位于管沟(6/8)L 和(7/8)L 处时，危险距离为零。当泄爆口大小变化时，人员轻伤的危险距离基本始终在3 m 左右小幅度波动。而当气云的长度越长时，人员轻伤的危险距离增加的幅值也越小，直至危险距离达到最大值。而当管沟的截面面积越大时，人员受伤的危险距离也越大。