石油与天然气化工  2025, Vol. 54 Issue (1): 149-157
变压吸附氢气纯化装置爆炸事故定量风险研究
刘迪1,2,3 , 王雅真1,2,3 , 邝辰1,2 , 顾蒙1,2     
1. 化学品安全全国重点实验室;
2. 中石化安全工程研究院有限公司;
3. 中石化国家石化项目风险评估技术中心有限公司
摘要目的 为合理确定变压吸附(PSA)氢气纯化装置爆炸风险及其周边建筑物的抗爆设防荷载,有效进行抗爆设计及防护。方法 以某变压吸附氢气纯化装置为例,采用计算流体动力学(CFD)模拟构建40个爆炸计算场景,获取氢气纯化装置爆炸风险三维分布,及其对人员集中建筑物的影响。结果 根据10−4 次/年的风险可接受准则和风险可接受范围内爆炸超压最大化准则,距离氢气纯化装置吸附塔橇85 m范围内人员集中建筑物应进行抗爆改造或设计;典型PSA装置北侧的人员集中建筑物墙的爆炸冲击波峰值入射超压为9.4 kPa,最大爆炸冲量为318 Pa·s,需要进行抗爆改造。结论 基于CFD的爆炸风险定量分析,考虑了装置实际布局、阻碍物等对氢气泄漏扩散、爆炸超压的影响作用,且滤除高后果低概率场景对最终爆炸设防荷载的影响,为PSA装置安全间距设置及周边人员集中建筑物抗爆改造或设计提供了依据。
关键词变压吸附装置    氢气爆炸    定量风险评估    抗爆载荷    
Quantitative risk study on explosion accident of PSA hydrogen purification unit
LIU Di1,2,3 , WANG Yazhen1,2,3 , KUANG Chen1,2 , GU Meng1,2     
1. State Key Laboratory of Chemical Safety, Qingdao, Shandong, China;
2. Sinopec Research Institute of Safety Engineering Co., Ltd., Qingdao, Shandong, China;
3. Sinopec National Petrochemical Project Risk Assessment Technical Center Co., Ltd., Qingdao, Shandong, China
Abstract: Objective The aim is to reasonably determine the explosion risk of pressure swing adsorption(PSA) hydrogen purification unit and the anti-explosion fortification load of surrounding buildings so as to support the anti-explosion design and protection. Method Taking PSA hydrogen purification unit as an example, CFD simulation was used to construct 40 explosion calculation scenarios to identify the three-dimensional distribution of explosion risk of hydrogen purification unit and its influence on the people gathering buildings. Result According to the acceptable risk criterion of 10−4 times per year and the maximum explosion overpressure criterion within the acceptable risk range, anti-explosion transformation or design should be carried out on the people gathering building within 85 m distance from the adsorption tower skids of the hydrogen purification unit. Besides, this study also showed that the peak incident overpressure of the explosion shock wave of the people gathering building wall on the north side of the typical PSA unit is 9.4 kPa, and the maximum explosive impulse was 318 Pa·s, which requires anti-explosion transformation. Conclusion With the consideration of the effects of the actual layout of the unit and obstacles on hydrogen leakage as well as diffusion and explosion overpressure, while eliminating the effects of scenarios with high consequences and low probability on the explosion protection value proof load, the quantitative analysis of explosion risk based on CFD in this study could provide an effective basis for the safe distance setting of PSA unit and the anti-explosion renovation or design of people gathering buildings.
Key words: pressure swing adsorption(PSA) unit    hydrogen explosion    quantitative risk assessment    anti-explosion load    

我国工业副产氢资源丰富,可为氢能产业提供稳定、低成本的氢源。然而副产氢气纯度较低、成分复杂,需要通过纯化技术才能将工业副产氢气转化为可使用的氢能,氢气纯化对于氢能在汽车领域的应用至关重要。变压吸附(PSA)是常见的气体提纯方法,其具有流程简单、自动化程度高,操作维修费用低,产品纯度可调性强等特点[1-2],成为我国石化企业普遍采用的氢气纯化方法。氢气具有易燃易爆,燃烧范围宽(体积分数4%~75%)、点火能量低、扩散系数大等特点[3-6]。同时,石化企业生产、储存过程通常涉及大量易燃易爆危险介质,且石化装置区阻塞程度高,一旦氢气纯化装置发生泄漏爆炸等严重事故,将导致周边人员集中建筑物及生产装置损毁。因此,须开展变压吸附氢气纯化装置泄漏爆炸风险研究,明确氢气泄漏燃爆危害范围,针对高风险场景采取有效的防护措施,以减缓事故后果。

目前,爆炸计算方法包括TNT 当量法、计算流体动力学(CFD)方法、TNO-ME(TNO多能法)和Baker-Strehlow-Tang(BST)模型等[7],部分学者采用上述计算方法开展了相关研究。刘奎[8]等通过研究TNO多能法,推导了爆炸后果与距离的关系式,用于评估石化企业爆炸对建筑物的爆炸影响。翟良云等[9]提出基于CFD的气体爆炸风险评估方法,结合密集管道与复杂装置布局、气云大小等因素影响,实现超压波的近场预测,并确定不同空间位置处的爆炸超压,更适用于石化行业的爆炸风险评估。夏兰生[10]对爆炸的特点、模型及其计算方法等进行研究和比较,同时结合爆炸力计算方面的相关实例,分析了适合工程设计使用的爆炸计算模型及计算方法。其中,TNT 当量法、TNO多能法、BST方法等均为基于经验公式的计算模型[11],无法考虑装置实际布局、阻碍物等对氢气泄漏扩散、爆炸超压的影响作用,导致计算结果与实际过程具有一定偏差。本研究采用CFD方法,利用FLACS Risk软件考虑复杂几何区域的通风、气体泄漏扩散和爆炸的过程,同时结合天气数据、泄漏频率、点火概率和事故后果减缓措施,对典型变压吸附氢气纯化装置爆炸风险进行模拟分析,建立了基于风险的爆炸荷载确定方法,为PSA装置安全间距及周边人员集中建筑物抗爆改造或设计提供了依据。

1 爆炸场景及发生频率
1.1 典型PSA工艺

典型变压吸附氢气纯化工艺流程分为原料缓冲、变压吸附、解吸气升压和氢气升压4个部分。某石化企业变压吸附氢气纯化工艺如图1所示,自原料气缓冲罐来的原料氢进入橇装设备的变压吸附塔,吸附脱除杂质组分,制取满足燃料电池用氢和高纯氢标准的产品氢气。

图 1     典型变压吸附氢气纯化工艺流程图 XV10101—进料切断阀;XV30201—氢气压缩机入口缓冲罐前切断阀;XV60101—解析气缓冲罐出口切断阀;D-101—原料气缓冲罐;RF-301A/B—旋转阀;C-301A~C-301L—变压吸附塔;D-601—解析气缓冲罐;D-901—氢气压缩机入口缓冲罐;K-901A/B—氢气压缩机;K-701A/B—解析气压缩机。

1.2 泄漏单元划分

PSA氢气纯化装置区有PSA橇装设备区、压缩机棚区,主要包括原料气缓冲罐、变压吸附塔橇、氢气压缩机入口缓冲罐、氢气压缩机、解析气缓冲罐、解析气压缩机等,相关工艺管线、工艺设备及其相关附件均可能发生失效泄漏。以发生泄漏时工艺前后所能关断或紧急切断的阀门位置为节点,参照设施的区域位置布置及操作参数(压力和温度等)的变化,将该装置整个工艺系统划分为不同的过程泄漏单元[12],划分结果见图1

1.3 爆炸场景发生频率

爆炸事故场景的发生需要考虑初始泄漏、风向条件、泄漏方向、点火源情况等,定义每个爆炸场景的发生频率计算公式如式(1)所示。

$ \mathit{P} _{ \mathrm{s}} \mathrm= \mathit{P} _{ \mathrm{w}} \mathit{P} _{ \mathrm{lf}} \mathit{P} _{ \mathrm{ld}} \mathit{P} _{ \mathrm{i}} $ (1)

式中:Ps为爆炸场景发生频率,次/年;Pw为该场景风向条件的概率;Plf为初始泄漏频率,次/年;Pld为泄漏方向的概率;Pi为延迟点火概率。

1.3.1 泄漏频率

美国桑迪亚实验室结合类似企业(炼化厂、核电站、海上石油工业等)设备泄漏历史数据,利用贝叶斯模型修正获取了临氢9类关键设备的失效频率[13]。每类设备包括5种泄漏类型(微小泄漏、较小泄漏、中等泄漏、较大泄漏、完全破裂),泄漏面积分别为总流通面积的0.01%、0.10%、1.00%、10.00%、100.00%。本研究基于桑迪亚实验室临氢设施失效数据,根据工艺流程图统计不同过程泄漏单元中包含9类设施(压缩机、储氢罐、管道、过滤器、法兰、阀门、软管、连接、仪表)的数量,得到对应典型泄漏孔径下各泄漏单元的泄漏频率,汇总见表1

表 1    各单元泄漏频率计算结果

1.3.2 点火概率

点火概率模型采用作为氢泄漏速率函数的点火概率值[14],具体见表2。该点火概率被美国桑迪亚实验室采用,作为HYRAM软件中点火概率,氢气泄漏直接点火概率在0.008~0.230之间,延迟点火概率在0.004~0.120之间。本研究依据不同场景下的泄漏速率获取对应爆炸场景的延迟点火概率。

表 2    高压氢气泄漏点火概率

2 典型爆炸事故后果影响分析
2.1 纯化装置物理模型

根据该PSA装置平面布局图及设备尺寸图,构建的装置整体模型如图2所示,由西向东依次为解析气压缩机、氢气压缩机、PSA橇装设备,PSA橇装设备南侧由东向西依次为原料气缓冲罐D-101、氢气压缩机入口缓冲罐D-901、解析气缓冲罐601等,利用FLACS RISK软件在三维笛卡尔坐标系中建立纯化装置1∶1三维模型。纯化装置,Y轴正方向为装置正北方向。重点关注氢气纯化装置周边北侧人员集中建筑物受到的爆炸危害后果,在该建筑物北墙、南墙、东墙、西墙表面设置监测点,监测PSA装置爆炸事故作用在该建筑物上的爆炸超压及冲量。

图 2     PSA装置及周边三维物理模型图

2.2 典型爆炸事故场景后果分析

表1中原料气缓冲罐D-101、变压吸附塔组C-301、解析气缓冲罐D-601、氢气压缩机入口缓冲罐D-901本体发生0.01%~100.00%泄漏后,延迟点火形成可燃气云爆炸的后果虽然很严重,但是发生频率非常低,以上几个储罐本体发生0.01%泄漏爆炸的频率在10−9~10−8 次/年,储罐本体发生100%泄漏爆炸的频率在10−8~10−7 次/年。因此,此处重点说明非储罐本体的泄漏爆炸事故后果。

PSA装置附近,平时人员暴露概率较低,爆炸超压直接作用导致人员伤亡风险较小。有统计表明,石油化工厂事故中,75%的人员死伤与蒸气云爆炸对建筑物的毁坏相关[15],多数情况下人员严重伤亡是由于建筑物坍塌或碎片飞溅造成的,因此,重点分析PSA装置爆炸事故对周边人员集中建筑的影响。

2.2.1 吸附塔橇装设备区泄漏爆炸

根据图1表1可知,泄漏单元2涉及的设施主要集中在PSA橇装设备区,相关管道、阀门典型尺寸为80 mm,当管道或连接处发生完全破裂即泄漏压力为2.1 MPa、温度为40 ℃、泄漏孔径为80 mm时,可燃气云体积最大为518 m3,可燃气云遇到点火源发生爆炸的危害后果如图3所示。

图 3     吸附塔橇装设备区泄漏爆炸后果云图

根据上述计算结果可知,吸附塔橇装设备区管道或连接处发生完全破裂泄漏爆炸后,超压由点燃位置处向四周传播,当遇到阻塞度较高的管道或设备(如吸附塔橇装设备、压缩机等)时,超压明显增强。这是由于障碍区域氢气可燃气云内混合气体的湍流程度较强,燃烧波面与未燃气体的接触面积较大,从而加快了气体化学反应速率,即燃烧速率[16],气体瞬间吸收大量燃烧产热后体积急剧膨胀,超压迅速上升。冲击波传播过程中,随着远离点燃处,超压逐渐减弱,直到最后完全消失。取PSA装置北侧人员集中建筑4面墙体最大爆炸超压对应监测点处的超压曲线,以及最大爆炸冲量监测点处的冲量曲线,具体如图4所示。PSA装置北侧人员集中建筑物承受最大爆炸超压为2.8 kPa(南墙),最大爆炸冲量为144 Pa·s(南墙),该建筑物受到该爆炸场景的危害后果较小。

图 4     吸附塔橇装设备区泄漏爆炸后果曲线图

2.2.2 氢气压缩机区泄漏爆炸

图1表1可知,泄漏单元3涉及的设施主要集中在氢气压缩机区,相关管道、阀门典型尺寸为50 mm,当管道或连接处发生完全破裂即泄漏压力为22 MPa、温度为40 ℃、泄漏孔径为50 mm时,可燃气云体积最大为916 m3,遇到点火源发生爆炸的危害后果如图5所示。

图 5     氢气压缩机区泄漏爆炸后果云图

取PSA装置北侧人员集中建筑4面墙体最大爆炸超压对应监测点处的超压曲线,以及最大爆炸冲量监测点处的冲量曲线,具体如图6所示。该建筑物承受最大爆炸超压为9.4 kPa(东墙),最大爆炸冲量为318 Pa·s(南墙),该建筑物受到该场景下的爆炸危害较大,未经抗爆设计时建筑物在该超压及冲量作用下结构受损,存在坍塌风险。

图 6     氢气压缩机区泄漏爆炸后果曲线图

2.2.3 解析气压缩机区泄漏爆炸

泄漏单元4即解析气缓冲罐出口切断阀XV60101—氢气管网涉及的设施主要集中在解析气压缩机区,相关管道、阀门典型尺寸为50 mm。当管道或连接处发生完全破裂即泄漏压力为2.2 MPa、温度为40 ℃、泄漏孔径为50 mm时,可燃气云体积最大为196 m3,遇到点火源发生爆炸的危害后果如图7所示。

图 7     解析气压缩机区泄漏爆炸后果云图

取PSA装置北侧人员集中建筑4面墙体最大爆炸超压对应监测点处的超压曲线,以及最大爆炸冲量监测点处的冲量曲线,具体如图8所示。PSA装置北侧人员集中建筑物承受最大爆炸超压为2.7 kPa(南墙),最大爆炸冲量为92 Pa·s(南墙),该建筑物受到该场景下的爆炸危害较小。

图 8     解析气压缩机区泄漏爆炸后果曲线图

3 爆炸超压风险分析

根据《中国石化既有建筑物抗爆治理指导意见(试行)》:当建筑物受到的爆炸冲击波超压≥6.9 kPa或者爆炸冲量≥207 kPa·ms时,建筑物需要进行抗爆设计。基于不同泄漏单元不同泄漏孔径对应的泄漏频率、结合不同泄漏速率下的延迟点火概率,同时考虑气象条件,计算得到所有爆炸场景的叠加风险结果,如图9图10所示。

图 9     1.0×10−4 次/年条件下爆炸超压6.9~50.0 kPa的影响范围

图 10     1.0×10−4 次/年条件下爆炸冲量207~2000 Pa·s的影响范围

PSA氢气纯化装置涉及到的危险物料为氢气,高压氢气泄漏后产生的可燃气体在高度阻塞区可能形成爆炸性云团,1.0×10−4 次/年频率条件下爆炸超压6.9 kPa最大影响范围约距离吸附塔橇85 m,爆炸冲量207 Pa·s最大影响范围约距离吸附塔橇68 m。根据人员集中建筑物的爆炸风险控制标准可知,距离吸附塔橇85 m范围内的人员集中建筑物应进行抗爆改造或设计,PSA装置北侧人员集中建筑物承受超压大于6.9 kPa,需进行抗爆改造治理。

图11可知,爆炸导致的死亡风险变压吸附由大到小为:塔橇装设备区>氢气压缩机棚区 > 解析气压缩机棚区。其中,变压吸附塔橇装设备区爆炸死亡风险为4×10−4~1×10−3 次/年,这是由于泄漏单元2中进料切断XV10101—氢气压缩机入口缓冲罐前切断阀XV30201/解析气缓冲罐出口切断阀XV60101之间的主要设施涉及阀门、连接、法兰等数量较多,导致该泄漏单元的初始泄漏频率较高,故爆炸场景发生频率高。氢气压缩机棚区由于氢气压力较高、存量较大,导致氢气泄漏速率较大、可燃气云较大,因而爆炸后果相对严重。

图 11     爆炸事故导致的死亡风险10−7~10−3次/年的分布

4 结论

1) 装置爆炸场景发生频率、可燃气云体积、阻塞度和爆炸源距离是影响爆炸设防荷载的关键变量。基于CFD模拟开展PSA装置爆炸风险定量计算,并结合风险可接受准则和风险可接受范围内爆炸超压最大化准则确定爆炸设防载荷,可以滤除高后果低概率场景对最终爆炸设防荷载的影响,避免因场景选择离散造成爆炸设防荷载计算结果的随机性。

2) 氢气纯化装置发生氢气泄漏爆炸事故后,1.0×10−4次/年频率条件下爆炸超压6.9 kPa最大影响范围约距离吸附塔橇85 m,爆炸冲量207 Pa·s最大影响范围约距离吸附塔橇68 m。根据人员集中建筑物的爆炸风险控制标准可知,距离吸附塔橇85 m范围内人员集中建筑物应进行抗爆改造或设计,PSA装置北侧人员集中建筑物承受超压大于6.9 kPa,需进行抗爆改造治理。

3) 高压氢气泄漏爆炸导致的死亡风险由大到小为:变压吸附塔橇装设备区>氢气压缩机棚区>解析气压缩机棚区。变压吸附塔橇装设备区涉及设施种类数量多,导致该区域初始泄漏频率较高,故爆炸场景泄漏频率高;氢气压缩机棚区氢气压力较高、存量较大,导致氢气泄漏速率较大、可燃气云较大,因而爆炸后果相对严重。

参考文献
[1]
氢气纯化变压吸附循环的热效应[J]. 工程热物理学报, 2018, 39(5): 1104-1111.
[2]
变压吸附提氢过程数值模拟研究[J]. 石油与天然气化工, 2017, 46(2): 27. DOI:10.3969/j.issn.1007-3426.2017.02.006
[3]
油氢合建站典型泄漏火灾热辐射后果模拟分析[J]. 消防科学与技术, 2021, 40(6): 851-854. DOI:10.3969/j.issn.1009-0029.2021.06.017
[4]
基于深度学习的管线泄漏仿真与测试平台[J]. 天然气与石油, 2022, 40(5): 10-16.
[5]
基于LMD-STA/LTA模型的油气管道泄漏检测方法[J]. 世界石油工业, 2022, 29(4): 71-76.
[6]
掺氢天然气长输管道泄漏扩散规律数值模拟[J]. 石油与天然气化工, 2023, 52(6): 44-52.
[7]
基于风险的液化烃罐区建筑物抗爆荷载评估方法[J]. 中国安全生产科学技术, 2020, 16(7): 30-35.
[8]
应用TNO多能法估算石油化工企业建筑物抗爆需求[J]. 炼油技术与工程, 2021, 51(3): 69-72. DOI:10.3969/j.issn.1002-106X.2021.03.018
[9]
石化行业控制室承爆风险评估方法研究[J]. 中国安全科学学报, 2009, 19(6): 129-134.
[10]
中心控制室抗爆模型选择与计算[J]. 石油化工自动化, 2015, 51(2): 7-11.
[11]
HU Q C, ZHANG X H, HAO H. A review of hydrogen-air cloud explosions: the fundamentals, overpressure prediction methods, and influencing factors[J]. International Journal of Hydrogen Energy, 2023, 48(36): 13705-13730. DOI:10.1016/j.ijhydene.2022.11.302
[12]
CPR18E. Purple book,guidelines for quantitative risk assessment[M]. The Hague,The Netherlands:Sdu uitgevers,1999.
[13]
LACHANCE J L,HOUF W G,FLUER L,et al. Analyses to support development of risk-informed separation distances for hydrogen codes and standards[R]. Albuquerque:Sandia National Laboratories,2009.
[14]
SUZUKI T, SHIOTA K, IZATO Y I, et al. Quantitative risk assessment using a Japanese hydrogen refueling station model[J]. International Journal of Hydrogen Energy, 2021, 46(11): 8329-8343. DOI:10.1016/j.ijhydene.2020.12.035
[15]
爆炸风险分析在建筑物抗爆设计中的应用[J]. 化工安全与环境, 2024, 37(7): 12-17.
[16]
加氢站高压氢气泄漏爆炸事故模拟及分析[J]. 浙江大学学报(工学版), 2015, 49(7): 1389-1394.