加氢站氢气泄漏事故模拟及后果分析.docxVIP

加氢站氢气泄漏事故模拟及后果分析 摘要：针对加氢站平安，通过理论模型分析和数值模拟两种方法，对 其开展事故模拟和后果分析。利用自行编制的MATLAB高斯扩散程序得到 爆炸危险区域的浓度曲线，分析环境风速对氢气扩散的影响，即风速越大, 危险区域越向泄漏口收缩；利用CFD软件Fluent建立加氢站氢气泄漏全 场景二维模型，模拟结果说明，无风情况下，氢气水平和垂直扩散速度很 快，容易富集并形成爆炸气团，而在风速10m/s情况下，泄漏氢气被带动、 吹散和稀释，难以富集，爆炸区域仅限于泄漏点附近。环境风不利于氢气 稳定扩散，对平安有利。 氢气被认为是一种可持续的、环境友好型的绿色能源，具有来源广泛、 燃烧热值高、可循环利用、储存方式多样等特点，被誉为21世纪的能源 之星。在氢能的众多前景中，氢燃料电池被认为是最有可能实现产业化的 应用之一。而伴随着各国氢燃料电池工程的开展实施，势必快速开展与之 配套的加氢站等根底设施。有专家预测，欧洲、美国和日本将在未来的五 到六年内实现氢燃料电池、氢能生产和加氢根底设施的商业化。我国燃料 电池技术虽刚刚起步，但在政府能源、环保战略的推进下，开展速度不断 加快。 国外学者研究说明，从能量利用的角度分析，高压储氢是最为经济合 理的选择，而由此也会带来较大的平安问题。高压氢气一旦发生泄漏，很 可能引起火灾和爆炸，造成巨大的人员伤亡和财产损失。因此，氢平安是 制约氢能开展和推广的主要瓶颈。 加氢站作为高压氢气储存较为集中的区域，其平安性必须受到高度的 重视。国际上常用的加氢站风险评价方法主要有快速风险评级和量化风险 评价，前者为定性评估，后者虽为定量评估，但在量化指标上存在较大争 议。国内浙江大学较早开展相关工作，研究了不同泄漏位置以及环境温度、 风速对高压储氢罐泄漏扩散的影响，但场景较为简单，不适合于真实、复 杂设施的事故模拟和后果分析。 考虑到氢气泄漏实验的危险性和高昂的费用，笔者选取某加氢站为工 程实例，采用理论模型分析和数值模拟(CFD)两种方法对其开展氢气泄漏 事故模拟和危险性分析，模拟结果可以为同类场所预防和抑制氢气扩散， 减小燃烧爆炸的可能性与破坏性提供工程依据和参考。 1原型描述 某高压氢气加氢站主要设备包括：大容量高压氢气管束车、无油隔膜 式氢压机、加气站高压储氢瓶组、高压氢气加注***、高压氢气专用阀件 和泄漏检测报警装置等。而其中，相比管束车、氢压机等流动或临时储氢 设施，高压储氢罐储氢时间长、储氢量大、压力高，瓶组的压力表、平安 阀一旦失灵或者超过使用压力，极易引起事故，造成燃烧爆炸。 加氢站采用分级充气的的传统压注方式，为了使车载储氢瓶内压力达 到35MPa,加氢站储氢罐的压力需要高达40~45MPa,才能获得快速加气。 笔者主要考察储氢罐发生泄漏的情况。具体参数见表1。这里不考虑 储氢罐内压力变化。 表1原型参数 氢气来源 工业副产氢气 氢气纯度 99.99% 泄漏位置 储氢罐顶部 泄漏压力 43 MPa 泄漏孔径 2 mm 2储氢罐泄漏和扩散理论模型 2. 1小孔泄漏模型 小孔泄漏模型适用于穿孔泄漏的情形，穿孔泄漏指管道或设备由于腐 蚀等原因形成小孔，燃气从小孔泄漏。穿孔直径在20nmi以下适用该模型。 小孔泄漏持续时间长，泄漏稳定且泄漏点多、不易发觉，潜在危险大。 在储存压力下，氢气以气相形式存在。气体从裂口泄漏的速率与其流 动状态有关，即首先判断气体流动属于音速还是亚音速流动。如果为音速 流动，那么有式(1)： P。c 2 P *《+ 1 式中：P为容器内介质压力，Pa;凡为环境压力，Pa;为 气体绝热指数，氢气取1.4o 气体泄漏量0(kg/s)计算公式见式(2): 式中：G为气体泄漏系数，泄漏口为圆形时取1 ,三角形 取0.95 ,长方形取0. 90 ;力为泄漏口面积jn2; M为相对 分子质量；R为气体常数,8.314 J/(mol ? K);『为气体 温度，K;*为氢气绝热系数。 如果为亚音速流动,那么有式(3): P。〉 2 P K+ 1 气体泄漏量0(kg/ s)计算公式见式(4): 0=EP 由］］ 式中:为泄漏口面积,m2 ;p为气体密度,kg/ n? ; Y为膨 胀因子，见式(5): Y=〔］1俨；：。「｛1.霹广｝ (5) 将原型参数代入上述公式得到，储氢罐顶部泄漏时, 气体属于音速流动，泄漏强度为0.082 8 kg/so 2. 2扩散模型 燃气泄漏后会在泄漏源附近形成气团，气团在大气中的扩散计算通常 采用高斯模型。高斯模型基于如下假设：(1)气体在扩散过程中没有沉降、 化合、分解以及地面吸收；(2)气体泄漏速率不随时间变化，即为稳定泄 漏；(3)扩散空间的风速、大气均匀稳定；(4)气云和环境无热交换；(5) 云团内部或云羽横截面上浓度、密度等参数服从正