基于MATLAB的煤气泄漏扩散高斯模型影响范围的研究

工业煤气基于MATLAB仿真的泄漏扩散影响研究

摘要：本文的研究目的是研究企业范围空间煤气泄漏的扩散规律和影响范围。采用matlab模拟煤气泄漏后CO 的浓度分布和扩散距离规律。通过建立煤气泄漏扩散数学模型, 对其影响煤气扩散的主要因素进行了分析、探讨了煤气毒性范围的划分, 然后在对煤气泄漏造成的危害和泄漏原因的基础上，运用扩散模型，计算煤气泄露扩散影响范围，然后用MATLAB对此进行模拟，得出不同的距离下煤气的浓度，并对其进行分析。因为大气稳定度、风速对煤气泄漏扩散的浓度影响起着非常重要的作用。大气稳定度和风速会显著改变有害气体的扩散状态。在风速和泄漏增大时, 煤气在开放空间扩散距离大, 影响范围广, 应合理布置煤气监控点, 预防煤气中毒。本文还鉴于煤气泄漏的危害之大，根据CO 的特性,对于煤气柜这种重大危险源的管理和控制可以得出一些经验，为采取措施预防其危害提供一定的依据。

关键词： 煤气泄漏； MATLAB；数值模拟； 扩散

一、前言

煤气泄漏的研究的背景及意义

我国当代工业以煤炭为主要能源的结构特点,决定了我国大多数工业企业的生产性气源以焦炉煤气和高炉煤气等为主,而煤气具有易燃易爆性!易散发性!剧毒性的特点,随着煤气在石油!化工!冶金等行业的广泛应用,也随之增加了煤气在工业场所发生泄漏!扩散并且导致人员中毒!火灾甚至爆炸发生的危险性和可能性\例如,2002年12月4日,天津西青开发区某厂房发生一起一氧化碳泄漏事故,造成3人中毒死亡;2005年2月22日,湖北大冶市一公司发生煤气中毒事故,当班的4名工人因中毒相继坠入料仓死亡;2005年4月21日,内蒙古自治区乌海市同力冶炼有限责任公司发生高炉煤气泄漏事故,造成2人中毒死亡;2005年n月5日,包头市大安钢铁公司发生煤气泄漏事故,当场造成5人中毒死亡,1人受伤;2005年的10月26日,首钢动力厂发生一起煤气中毒事故,共有9人丧生;而时隔8个月,即2006年6月10日首钢动力厂再次煤气泄漏事故,至少有7人中毒,其中2人经抢救无效死亡\此类事故举不胜举\

近几年来市场对煤气及其相关产品的需求增大,企业不断扩大生产能力,同时煤气事故的次数也居高不下,鉴于以上事实,我们发现:工业场所煤气一旦发生事故性泄漏,往往会酿成人员中毒伤亡的严重后果,另外,若遇火源还可能导致火灾或爆炸等事故造成重大损失\因此,为减少因煤气事故泄漏事故带来的人员及财产损失,对工业场所煤气的泄漏!扩散进行数值模拟分析,加强对其微观规律的研究,为制定相应的煤气中毒预防及事故减灾策略有重要的理论意义\近年来我国工业煤气事故性泄漏屡有发生,尤其严重的是2005年和2006年首钢动力厂连续两次发生煤气泄漏事故,并造成重大人员伤亡,此事件发人深省\其重要原因之一就是人们对工业场所煤气泄漏扩散的规律不甚了解,尤其是煤气泄漏扩散后中毒伤害范围的变化,安全警戒撤离距离的确定等信息不能及时获得,从而延误了中毒区域内人员的救援时机,造成重大人员和财产损失\工业场所煤气泄漏扩散是一个综合而又复杂的过程,泄漏物质,泄源高度及面积!泄漏速度!泄漏时间!大气稳定度!地形等参数对扩散都有着重要的影响\因此,如何对工业场所煤气泄漏扩散的过程进行有效的模拟,以及时!准确!有效地获得各种参数,为煤气泄漏事故的应急救援提供科学依据就显得十分迫切\

国内外的研究现状

\国外在这方面的研究相对成熟,直到现在该领域的研究还比较活跃\国外学者提出了不少扩散的计算模型,同时也进行了许多大规模试验\主要的数值扩散模型有高斯(Gaussianplume/Puffmodel),BM(BritterandMeQuaid)模型Sutton模型 DEM(3一DFiniteElementModel)等等\高斯模型适用于点源的扩散,早在五六十年代就己经被应用\它从统计方法入手,考察扩散介质的浓度分布,适用于中等密度气团(非重气)扩散的模拟\烟羽模型(Plumemodel)适用于连续源的

扩散,烟团模型(Puffmodel)适用于短时间泄漏的扩散(即泄放时间相对于扩散时间比较短的情形,如突发性泄漏等)\高斯模型具有简单,易于理解,运算量小的特点,且由于提出的时间比较早,实验数据多,因而较为成熟\

高斯(Gauss)模型属于非重气扩散模型,只适用于与空气密度相差不多的气体扩散\但是,大多数危险性物质一旦泄漏到大气环境中就会由于较重的分子质量(如C12)低温和化学变化(如HF)等原因形成比周围环境气体重的重气云,重气云的扩散机理与非重气云完全不同\因此,重气云扩散机理的研究是国外众多学者竞相研究的热点课题\国际上曾多次召开有关重气云扩散研究及其预防控制方面的系列学术会议,促进了重气云扩散的研究\到目前为止,已提出大约200个重气云扩散模型\重气云扩散模型可分为经验模型、箱模型、浅层模型以及三维流体力学模等等\

随着计算机的普及和计算能力的不断提高,加上近似计算方法,例如,有概述限差分法、有限元法、有限体积法等的发展,基于数值计算的计算流体力学(ComputationalFluidD0amics,CFD)方法形成并得到了迅速的发展\正是England等(1978年)触发了采用CFD方法模拟重气扩散的三维非定常态湍流流动过程\这种数值方法是通过建立各种条件下的基本守恒方程(包括质量、动量、能量及组分等),结合一些初始和边界条件,加上数值计算理论和方法,从而实现预报真实过程各种场的分布,例如,流场、温度场、浓度场等,以达到对扩散过程的详细描述\用这种方法就克服了箱及相似模型中辨识和模拟重气的下沉、空气的卷吸、气云的受热等各种物理效应时所遇到的许多问题。

我国关于有毒有害气体泄漏、扩散的研究起步较晚,始于九十年代初期,且投入力量尚不多八五期间,中国石油化工集团公司安全工程研究院(原化工部劳动保护研究所)对有毒物质的泄漏扩散进行了研究,在总结和建立有毒物质泄漏模式及泄漏源模型的基础上提出了泄漏扩散的HLY模型,该模型与国外一些著名的泄漏扩散模型相比与试验数据的吻合程度更好。另外,北京市危险源控制技术研究中心自1997年成立以来,一直将泄漏计算机仿真技术研究作为其主攻方向之一,不断进行泄漏扩散模型和风洞模拟实验研究,先后完成了有毒物质泄漏扩散模型及监控技术研究有毒重气泄漏模拟实验与扩散模型研究、毒物泄漏扩散模型研究等科研课题,并将部分成果应用于一些环境卫生和安全与评价项目;大连理工大学在国家自然科学基金的资助下,通过对现有的扩散模型进行了研究分析,从气体动力学入手,通过对气体微元进行质量平衡、动量平衡、能量平衡的分析,采用平板模型并在风洞实验的基础上对可燃及毒性气体的扩散过程进行了研究；北京化工大学对重气扩散的过程进行了分析,采用涡粘性模型中的k一双方程模型对模型中的扩散系数进行了修正,并对整个重气扩散过程进行了模拟,模拟结果与试验结果符合程度较好;南京工业大学在国家自然科学基金重点项目的资助下对泄漏物质在大气环境中的扩散过程也进行了较系统深入的研究。此外,北京市劳动保护科学研究所、东北大学、北京理工大学等单位也都开展了这方面的研究,也取得了一定的成果\

二、煤气扩散危险区域模拟的基本理论

煤气的毒性

易燃易爆：煤气中含有大量的一氧化碳等可燃气体，见明火或着火点以上的高温就会燃烧，而且易中毒。一氧化碳是一种无色无味无刺激性的气体，是一种易燃易爆的气体么一氧化碳毒性很大，工业卫生标准规定空气中的一氧化碳最高允许浓度为30mg/m。当空气中一氧化碳浓度达到250mg/m是2-3个小时极有可能出现中毒症状，空气中得一氧化碳浓度达到1000mg/m时两小时可导致昏迷，如果浓度再高危险则越大，严重者当场死亡。

如果煤气泄漏遇到助燃物加明火、静电、闪电或操作不当等都会引起爆炸、燃烧、火灾，在密闭空间会使人缺氧、窒息，甚至死亡，给公司的安全生产和国家及人民生命财产带来不可估量的损失。

2、一般情况下,煤气柜区煤气泄漏的原因有以下几点:

1）、煤气设备和管道发生破损、各处水封由于系统压力突然升高,可能会造成跑煤气现象从而导致操作人员煤气中毒。

2）、加压机管道连接处、阀门处密封不严密,煤气泄漏有可能导致中毒事故。

3）、煤气柜放散系统及燃烧放散管,遇煤气压力波动及自然气候的低气压天气,煤气放散管烧嘴因故障熄灭,未能及时发现(排除) 隐患,可能引起煤气放散下风侧多人中毒的重大伤亡事故。

4）、煤气柜和管道都用钢材建设,储存的气体对气柜和管道有一定的腐蚀作用,在炼厂冶炼废气较多的环境下,外部的SO2 ,H2S ,CO2 均对钢结构有较强的腐蚀作用,特别是SO2 腐蚀危害对钢结构有较大影响。除了大气和煤气中所带的腐蚀因素之外,煤气柜体焊缝等处由于金属与金属、金属与非金属因内外氧浓差不同形成浓差电池,缝隙内为阳极易被腐蚀产生泄漏。大部分冶炼厂中的煤气柜系统均需隔一到二年即进行防腐蚀处理,否则将会因为腐蚀而出现煤气到处泄漏的情况

本文来源：https://www.bwwdw.com/article/k4o3.html