液氨泄漏扩散的模拟讲义.pptVIP

液氯泄漏扩散过程模拟研究 青岛科技大学 一、研究内容及意义 1.内容： 分析案例中相关条件建立数学模型； 用MATLAB对液氯泄漏扩散的具体案例进行模拟计算； 由计算结果绘出气体泄漏扩散浓度随距离的变化曲线图及危险区域分布图 2.意义： 通过结果分析熟悉液氯泄漏扩散的浓度变化规律； 为生产企业应急救援提供技术支持 二、非重气云团扩散模型简介 液氯泄漏扩散后，由于氯气的密度比空气大，属重气扩散，在扩散过程中，随着空气的不断卷入，氯气云团逐渐被稀释，当其密度与周围空气密度的差小于0.001Kg/m3时，认为氯气云团已转变为非重气云团。 这里只对液氯泄漏扩散后转变为中性气云后的扩散行为进行研究 ，可用于研究其点源泄漏扩散的模型有高斯烟羽模型和高斯烟团模型。 高斯烟羽模型：又称高架点连续点源扩散模型，适用于连续源的扩散，即连续源或泄放时间大于或等于扩散时间的扩散。 高斯烟团模型：适用于短时间泄漏的扩散，即泄放时间相对于扩散时间比较短的情形，如突发性泄放等。 三、案例计算 某化工厂室外高塔架上有一氯气储罐，该储罐最大可存储液氯50t，储罐直径4m,高7m,由于罐体破裂发生连续性泄漏，泄漏口当量孔径为3cm,泄漏口之上液位高度4m,当地大气压力为0.1MPa,平均风速为3m/s.试估算在不同大气稳定度下液氯泄漏可能的影响范围。 大气稳定度的分类 大气稳定度由10米高度以上的风速、白天的太阳辐射或夜间的云量等参数决定, Pasquill的分类方法如下（云量是指当地天空云层覆盖率）： 模型确定依据 （1）液氯泄漏后，首先形成氯气云团进行扩散，随着空气的不断卷吸，气云不断被稀释，最终转变为中性气云，此时便可采用高斯模型。 （2）该事故是高架储罐破裂孔洞发生的点源的连续性泄漏，故应选取高斯烟羽模型，即高架点的连续点源扩散模型。 高斯烟羽模型公式 其表达式为： 中 式中： X（x,y,z):下风向x米、横向y米、地面上方z米的扩散的气体浓度,kg/m3; Q:源强（即源释放速率 ),kg/s; u:平均风速,m/s; σy::水平扩散参数,m; σz：垂直扩散参数,m; t:泄漏后的时间,s; H:泄漏源有效高度,m; y:横向距离,m; z:垂直方向距离,m 扩散系数的选取 扩散系数的计算方法如下： 液氯泄漏后仅沿风向扩散 设风速方向为x轴方向，泄漏源中心在地面的投影为坐标原点，建立高斯烟羽模型坐标系。令y=0，采用高斯烟羽模型公式，运用MATLAB编程计算泄漏气体沿x轴（风向）扩散的气体浓度，输出结果如下所示： 大气稳定度为C： 质量浓度/kg/m3 大气稳定度为A 大气稳定度为B 大气稳定度为E 大气稳定度为F 质量浓度/Kg/m3 质量浓度/Kg/m3 液氯泄漏后沿风向和垂直风向扩散 仍选用高斯烟羽模型，用MATLAB软件编程计算输出结果如下（大气稳定度为C）: 质量浓度/Kg/m3 液氯泄漏扩散的危险区域分布 液氯泄漏扩散的危险区域划分为致死区、重伤区、轻伤区和吸入反应区，用MATLAB编程计算，输出结果如下（大气稳定度为B）： 大气稳定度为A 大气稳定度为C 大气稳定度为E 大气稳定度为F 四、结论 同一液氯泄漏源在不同大气稳定度下泄漏后扩散的浓度随距离变化曲线均呈正态分布，随着扩散的进行，氯气浓度值先增大后减小； 同一液氯泄漏源在不同大气稳定度下随着A至F大气稳定性的增强，其扩散浓度的峰值逐渐减小，而出现浓度峰值的点距泄漏源的距离逐渐增大； 同一液氯泄漏源泄漏后在较大尺度空间的扩撒浓度峰值比在较小尺度空间扩散的峰值小； 同一液氯泄漏源在不同大气稳定度下泄漏后扩散的危险区域随着A至F大气稳定性的增强呈逐渐扩大的趋势。