《基于MATLAB的煤气泄漏扩散高斯模型影响范围的研究.doc》由会员分享，可在线阅读，更多相关《基于MATLAB的煤气泄漏扩散高斯模型影响范围的研究.doc（14页珍藏版）》请在金锄头文库上搜索。
1、工业煤气基于MATLAB仿真的泄漏扩散影响研究摘要：本文的研究目的是研究企业范围空间煤气泄漏的扩散规律和影响范围。采用matlab模拟煤气泄漏后CO 的浓度分布和扩散距离规律。通过建立煤气泄漏扩散数学模型, 对其影响煤气扩散的主要因素进行了分析、探讨了煤气毒性范围的划分, 然后在对煤气泄漏造成的危害和泄漏原因的基础上，运用扩散模型，计算煤气泄露扩散影响范围，然后用MATLAB对此进行模拟，得出不同的距离下煤气的浓度，并对其进行分析。因为大气稳定度、风速对煤气泄漏扩散的浓度影响起着非常重要的作用。大气稳定度和风速会显著改变有害气体的扩散状态。在风速和泄漏增大时, 煤气在开放空间扩散距离大, 影响
2、范围广, 应合理布置煤气监控点, 预防煤气中毒。本文还鉴于煤气泄漏的危害之大，根据CO 的特性,对于煤气柜这种重大危险源的管理和控制可以得出一些经验，为采取措施预防其危害提供一定的依据。关键词： 煤气泄漏； MATLAB；数值模拟； 扩散一、前言煤气泄漏的研究的背景及意义我国当代工业以煤炭为主要能源的结构特点,决定了我国大多数工业企业的生产性气源以焦炉煤气和高炉煤气等为主,而煤气具有易燃易爆性!易散发性!剧毒性的特点,随着煤气在石油!化工!冶金等行业的广泛应用,也随之增加了煤气在工业场所发生泄漏!扩散并且导致人员中毒!火灾甚至爆炸发生的危险性和可能性例如,2002年12月4日,天津西青开发区某
3、厂房发生一起一氧化碳泄漏事故,造成3人中毒死亡;2005年2月22日,湖北大冶市一公司发生煤气中毒事故,当班的4名工人因中毒相继坠入料仓死亡;2005年4月21日,内蒙古自治区乌海市同力冶炼有限责任公司发生高炉煤气泄漏事故,造成2人中毒死亡;2005年n月5日,包头市大安钢铁公司发生煤气泄漏事故,当场造成5人中毒死亡,1人受伤;2005年的10月26日,首钢动力厂发生一起煤气中毒事故,共有9人丧生;而时隔8个月,即2006年6月10日首钢动力厂再次煤气泄漏事故,至少有7人中毒,其中2人经抢救无效死亡此类事故举不胜举近几年来市场对煤气及其相关产品的需求增大,企业不断扩大生产能力,同时煤气事故的次
4、数也居高不下,鉴于以上事实,我们发现:工业场所煤气一旦发生事故性泄漏,往往会酿成人员中毒伤亡的严重后果,另外,若遇火源还可能导致火灾或爆炸等事故造成重大损失因此,为减少因煤气事故泄漏事故带来的人员及财产损失,对工业场所煤气的泄漏!扩散进行数值模拟分析,加强对其微观规律的研究,为制定相应的煤气中毒预防及事故减灾策略有重要的理论意义近年来我国工业煤气事故性泄漏屡有发生,尤其严重的是2005年和2006年首钢动力厂连续两次发生煤气泄漏事故,并造成重大人员伤亡,此事件发人深省其重要原因之一就是人们对工业场所煤气泄漏扩散的规律不甚了解,尤其是煤气泄漏扩散后中毒伤害范围的变化,安全警戒撤离距离的确定等信息
5、不能及时获得,从而延误了中毒区域内人员的救援时机,造成重大人员和财产损失工业场所煤气泄漏扩散是一个综合而又复杂的过程,泄漏物质,泄源高度及面积!泄漏速度!泄漏时间!大气稳定度!地形等参数对扩散都有着重要的影响因此,如何对工业场所煤气泄漏扩散的过程进行有效的模拟,以及时!准确!有效地获得各种参数,为煤气泄漏事故的应急救援提供科学依据就显得十分迫切国内外的研究现状国外在这方面的研究相对成熟,直到现在该领域的研究还比较活跃国外学者提出了不少扩散的计算模型,同时也进行了许多大规模试验主要的数值扩散模型有高斯(Gaussianplume/Puffmodel),BM(BritterandMeQuaid)模
6、型Sutton模型 DEM(3一DFiniteElementModel)等等高斯模型适用于点源的扩散,早在五六十年代就己经被应用它从统计方法入手,考察扩散介质的浓度分布,适用于中等密度气团(非重气)扩散的模拟烟羽模型(Plumemodel)适用于连续源的扩散,烟团模型(Puffmodel)适用于短时间泄漏的扩散(即泄放时间相对于扩散时间比较短的情形,如突发性泄漏等)高斯模型具有简单,易于理解,运算量小的特点,且由于提出的时间比较早,实验数据多,因而较为成熟高斯(Gauss)模型属于非重气扩散模型,只适用于与空气密度相差不多的气体扩散但是,大多数危险性物质一旦泄漏到大气环境中就会由于较重的分子质
7、量(如C12)低温和化学变化(如HF)等原因形成比周围环境气体重的重气云,重气云的扩散机理与非重气云完全不同因此,重气云扩散机理的研究是国外众多学者竞相研究的热点课题国际上曾多次召开有关重气云扩散研究及其预防控制方面的系列学术会议,促进了重气云扩散的研究到目前为止,已提出大约200个重气云扩散模型重气云扩散模型可分为经验模型、箱模型、浅层模型以及三维流体力学模等等随着计算机的普及和计算能力的不断提高,加上近似计算方法,例如,有概述限差分法、有限元法、有限体积法等的发展,基于数值计算的计算流体力学(ComputationalFluidD0amics,CFD)方法形成并得到了迅速的发展正是Engl
8、and等(1978年)触发了采用CFD方法模拟重气扩散的三维非定常态湍流流动过程这种数值方法是通过建立各种条件下的基本守恒方程(包括质量、动量、能量及组分等),结合一些初始和边界条件,加上数值计算理论和方法,从而实现预报真实过程各种场的分布,例如,流场、温度场、浓度场等,以达到对扩散过程的详细描述用这种方法就克服了箱及相似模型中辨识和模拟重气的下沉、空气的卷吸、气云的受热等各种物理效应时所遇到的许多问题。我国关于有毒有害气体泄漏、扩散的研究起步较晚,始于九十年代初期,且投入力量尚不多八五期间,中国石油化工集团公司安全工程研究院(原化工部劳动保护研究所)对有毒物质的泄漏扩散进行了研究,在总结和建
9、立有毒物质泄漏模式及泄漏源模型的基础上提出了泄漏扩散的HLY模型,该模型与国外一些著名的泄漏扩散模型相比与试验数据的吻合程度更好。另外,北京市危险源控制技术研究中心自1997年成立以来,一直将泄漏计算机仿真技术研究作为其主攻方向之一,不断进行泄漏扩散模型和风洞模拟实验研究,先后完成了有毒物质泄漏扩散模型及监控技术研究有毒重气泄漏模拟实验与扩散模型研究、毒物泄漏扩散模型研究等科研课题,并将部分成果应用于一些环境卫生和安全与评价项目;大连理工大学在国家自然科学基金的资助下,通过对现有的扩散模型进行了研究分析,从气体动力学入手,通过对气体微元进行质量平衡、动量平衡、能量平衡的分析,采用平板模型并在风
10、洞实验的基础上对可燃及毒性气体的扩散过程进行了研究；北京化工大学对重气扩散的过程进行了分析,采用涡粘性模型中的k一双方程模型对模型中的扩散系数进行了修正,并对整个重气扩散过程进行了模拟,模拟结果与试验结果符合程度较好;南京工业大学在国家自然科学基金重点项目的资助下对泄漏物质在大气环境中的扩散过程也进行了较系统深入的研究。此外,北京市劳动保护科学研究所、东北大学、北京理工大学等单位也都开展了这方面的研究,也取得了一定的成果二、煤气扩散危险区域模拟的基本理论煤气的毒性易燃易爆：煤气中含有大量的一氧化碳等可燃气体，见明火或着火点以上的高温就会燃烧，而且易中毒。一氧化碳是一种无色无味无刺激性的气体，是
11、一种易燃易爆的气体么一氧化碳毒性很大，工业卫生标准规定空气中的一氧化碳最高允许浓度为30mg/m。当空气中一氧化碳浓度达到250mg/m是2-3个小时极有可能出现中毒症状，空气中得一氧化碳浓度达到1000mg/m时两小时可导致昏迷，如果浓度再高危险则越大，严重者当场死亡。如果煤气泄漏遇到助燃物加明火、静电、闪电或操作不当等都会引起爆炸、燃烧、火灾，在密闭空间会使人缺氧、窒息，甚至死亡，给公司的安全生产和国家及人民生命财产带来不可估量的损失。2、一般情况下,煤气柜区煤气泄漏的原因有以下几点:1）、煤气设备和管道发生破损、各处水封由于系统压力突然升高,可能会造成跑煤气现象从而导致操作人员煤气中毒。
12、2）、加压机管道连接处、阀门处密封不严密,煤气泄漏有可能导致中毒事故。3）、煤气柜放散系统及燃烧放散管,遇煤气压力波动及自然气候的低气压天气,煤气放散管烧嘴因故障熄灭,未能及时发现(排除) 隐患,可能引起煤气放散下风侧多人中毒的重大伤亡事故。4）、煤气柜和管道都用钢材建设,储存的气体对气柜和管道有一定的腐蚀作用,在炼厂冶炼废气较多的环境下,外部的SO2 ,H2S ,CO2 均对钢结构有较强的腐蚀作用,特别是SO2 腐蚀危害对钢结构有较大影响。除了大气和煤气中所带的腐蚀因素之外,煤气柜体焊缝等处由于金属与金属、金属与非金属因内外氧浓差不同形成浓差电池,缝隙内为阳极易被腐蚀产生泄漏。大部分冶炼厂中
13、的煤气柜系统均需隔一到二年即进行防腐蚀处理,否则将会因为腐蚀而出现煤气到处泄漏的情况3、煤气扩散过程中的主要影响因素( 1)泄漏源的情况。一般情况下, 可以根据泄漏面积的大小和泄漏持续时间的长短, 将泄漏源进行分类。泄漏包括: 泄漏位置, 裂孔的形状、尺寸; 气体压力和温度等; 事故发生地点及周围的地理环境; 泄漏形式是连续泄漏或是瞬时泄漏。( 2)泄漏物质的理化性质, 如分子量、沸点、闪点、比重、密度。( 3)环境条件, 包括大气稳定度, 风速、风向、气温以及泄漏源周围的地形地貌。大气稳定度表征湍流活动的强弱, 支配大气对污染物的稀释扩散能力,大气越稳定, 气云越不易向高空消散, 而贴近地表
14、扩散, 大气越不稳定, 空气垂直对流运动越强, 气云消散得越快; 风向决定泄漏气体的扩散方向, 风速决定泄漏气体的扩散速度。大气湿度对泄漏气体的扩散也有影响, 通常情况下, 湿度大则使气体不易扩散;气温或太阳辐射强弱主要是通过影响大气垂直对流运动而对气体的扩散产生影响。二、煤气扩散高斯模型从应用角度出发，目前国内外在实际操作中使用的大气扩散模型，多属描述中等密度云扩散的高斯模型及其变形模型。国家制定的制定地方大气污染物排放标准的技术方法)(GBT1320191)亦推荐该模型。由于在毒气扩散的全过程中，中等密度气云在大气中达到最大伤害面积。持续时间最长，通常将这个过程作为标绘重点。高斯模型又称高
15、架连续点源扩散模型，适用于连续点源扩散和中等密度毒气云团的伤害范围估算。当泄露毒气或毒气与空气混合后的密度接近空气密度时，重力下沉与浮力上升作用可以忽略，扩散主要是由空气湍流决定，在假设均匀湍流场的条件下，根据统计理论。有害物质在扩散截面的浓度呈高斯分布，所以称为高斯扩散。由于该模型建立比较早，数学运算简单，实验数据充分，因此得到广泛应用。在稠密气团(重气)扩散场合。由于重力作用时间非常短暂，也可以直接采用高斯模型。高斯模型适用条件如下：(1)云团下垫面近似水平，开阔、性质均匀；(2)扩散过程中污染物不发生化学反应，没有干沉降、降水清洗等衰减作用；(3)云团与空气相对运动可以忽略不计，即随大气
16、一起运动，平均流场平直稳定，平均风速和风向没有先住变化，平均风速不小于lms；(4)适用于尺度不超过20km的扩散。因此，设释放源投影点为坐标原点，z轴指向风向；考虑地面反射作用，可得到烟羽模型的浓度分布计。 式（1）式中， 下风向某点处的空气污染物的浓度；Q：源强即气体的流量或浓度，kgS；He：毒气源的有效高度，m；u：毒气源所在高度的平均风速，ms； Y：侧风方向离泄露源点的距离，m； z：垂直向上方的离泄露源点的距离，m。分别为水平方向和垂直方向的扩散参数，他们是下风距离x及大气稳定度的函数。三、MATLAB软件介绍20世纪80年代初期，Moler等一批数学家与软件专家组建了Mathworks软件开发公司，从事MATLAB的研究和开发， 1984年MathworkS公司推出第一个MATLAB商业版本，其核心是用C语言
