一种基于网格化毒气扩散模拟的方法
【专利摘要】本发明公开了一种基于网格化毒气扩散模拟的方法，其特征在于，基于区域网格划分技术及高斯烟团模型，综合储罐泄露速度、有效源高（泄露高度），以及风向、风速、环境稳定级多因素变化，研究区域内一个或多个储罐中有毒有害物在连续时间内浓度分布规律。主要包括首先，以储罐为中心，将研究区域网格化为若干单元作为定量研究和叠加的基础；然后，运用高斯烟团理论构建有毒有害物的泄露扩散浓度分布模型；最后，结合单泄露源和多泄露源条件，给出有毒有害物连续性扩散算法，即根据任影响因素的变化来划分时间段，按时间顺序，结合浓度叠加效应，基于Matlab模拟区域内有毒有害物浓度过程。本发明可有效地为工业化工园区复杂条件下，有毒有害物突发性泄漏事故应急决策提供理论支持。可广泛用于化工园区规划和应急事故处理。
【专利说明】一种基于网格化毒气扩散模拟的方法【技术领域】
[0001]本发够涉及化工园区毒气扩散研究,録;g/是涉及基于网格化毒气扩散模拟的方法。
【背景技术】
[0002]城市化工园区生产、储存、运输着大量易燃、易爆、有毒的危险化学品。园区内化学储罐多，种类繁杂，安全性脆弱，一旦发生泄漏事故，大量的有毒、易挥发、易燃性物质迅速释放到大气中，会形成气云并逐渐扩散，将对周围人员和环境造成巨大危害。更严重的是，有毒有害物质可能在很低的浓度下就会对人体造成伤害，从而酿成重大事故，据统计2000年到2008年间，公安消防部分参与处置的化学灾害事故35438起，平均每周处置75起、每天处置11起。
[0003]目前，对化学工业园区有毒有害物扩散造成重特大事故的研究，主要集中在环境条件相对固定的单一危险源有害物质扩散区域浓度分析，而对由于单个储罐火灾爆炸事故诱发的多米诺连锁事故研究相对较少。基于以上考虑，重点研究多危险源事故期间，在泄漏量、环境平均风速、泄漏时间、有效源高、环境稳定等级及风向等影响因素耦合作用条件下，区域内有害物质连续扩散浓度分布情况。运用Matlab技术，研究多因素多储罐高斯烟团模型的修正问题，通过模型修正设计、分析、叠加等技术手段，找出灾变时期储罐区域有毒有害物质浓度分布规律，实现对整个区域进行合理规划设计，达成有效的事故应急救援方案。
【发明内容】
[0004]针对化工园区毒气泄漏问题，对于其模拟方面存在的问题，本发明提出一种基于网格化毒气扩散模拟的方法。
[0005]1.网格划分理论及高斯烟团模型改造
化工园区有害物扩散浓度模拟的目的是给出区域内各位置有害物浓度大小。为计算需要，将研究区域进行网格化，即将评价区域划分成等步长正方形网格，以mXn 二矩阵表示这个笛卡尔坐标系下二维空间。首先计算区域内各单元的有害物浓度，然后进行叠加，其余网格坐标值[x，y]通过离散与内插计算得到。对于步长的选择可根据区域大小、制图精度和实际需要等进行调整。一般步长选择25~300m较为合适，步长的大小直接影响评价结果的精度，步长越小则精度越高，但计算相对耗时。
[0006]基础模型为高斯烟团模型为瞬时泄漏模型，其数学模型为:
【权利要求】
1.一种基于网格化毒气扩散模拟的方法，差藍您杰基于区域网格划分技术及高斯烟团模型改造，综合储罐泄露速度、有效源高(泄露高度)，以及风向、风速、环境稳定级多因素变化，确定区域内一个或多个储罐中有毒有害物在连续时间内浓度分布规律，其包蕾勿下步骤:售先,以储罐为中心，将研究区域网格化为若干单元作为定量研究和叠加的基础；然后，运用高斯烟团理论构建有毒有害物的泄露扩散浓度分布模型；最后，结合单泄露源和多泄露源条件，给出有毒有害物连续性扩散算法，即根据任影响因素的变化来划分时间段，按时间顺序，结合浓度叠加效应，基于Matlab模拟区域内有毒有害物浓度过程，本I够可有效地为工业化工园区复杂条件下，有毒有害物突发性泄漏事故应急决策提供理论支持。
2.根据权利要求1所述的区域网格划分技术，其特征在于，化工园区有害物扩散浓度模拟的目的是给出区域内各位置有害物浓度大小，为计算需要，将研究区域进行网格化，即将评价区域划分成等步长正方形网格，以mXn且m=n 二矩阵表示这个笛卡尔坐标系下二维空间，首先计算区域内各单元的有害物浓度，然后进行叠加，其余网格坐标值[X，y]通过离散与内插计算得到；对于步长的选择可根据区域大小、制图精度和实际需要等进行调整；一般步长选择25~300m。
3.根据权利要求1所述的高斯烟团模型改造，其特征在于，使用的计算瞬时泄漏的公式是对高斯烟团模型的改造，具体方法通过编程实现，公式说明如下:
4.根据权利要求1所述的单泄露源模拟，其特征在于，对于单一危险源泄漏的解决方案:设某一时段t1;设不变，那么对于单一泄漏源A运用公式(2)可得到经过时段h后的,在下一时段t2时，如果其中之一有变化，那么对于继续泄漏的单一危险源A重新根据公式(2)计算，得到，但是在t2时间段内，中每个单元内的有害气体在KEC, IA,"条件下都产生扩散，在这种情况下设:
【文档编号】G06F19/00GK103793591SQ201310029240
【公开日】2014年5月14日 申请日期:2013年1月27日 优先权日:2013年1月27日
【发明者】刘涛, 李颖, 张洋 申请人:辽宁工程技术大学