第五章：毒物泄漏及扩散模型-第四次

1化工安全工程Chemicalengineeringsafetyfundamentalandapplication邹海魁Email:[e-mail]:64443134第四章泄漏源模型一、源模型简介依据描述物质释放时所表现出的物理化学过程的理论，或传递过程理论及经验方程而建立。作用：源模型给出了流出速率、流出总量(总时间)和流出状态。二、释放机理大孔和有限孔第四章泄漏源模型三、7个基本模型液体经孔洞流出液体经贮罐上的孔洞流出液体经管道流出蒸汽经孔洞流出，塞压气体经管道流出闪蒸液体液池蒸发或沸腾第五章：毒物泄漏及扩散模型主要内容：•有毒物质泄漏及扩散模型的意义、步骤、作用。•扩散影响参数。•中性浮力扩散模型、14个例子；•Pasquill-Gifford模型•重气扩散•毒性作用标准•释放动量和浮力的影响•释放缓解前言1、有毒物质泄漏扩散模型的意义•事故期间，过程单元或设备释放出大量有毒物质，形成的危险气云可能传播到整个工厂区域和当地社区。•重大化学品事故（如Bhopal）的发生及后果，使人们意识到应急计划的重要性，将工厂设计成毒物释放事故发生最少、事故后果最小化的重要性。•需要毒物释放模型，描述其迁移转化过程规律。2、作用毒物释放和扩散模型是后果模拟步骤的重要部分；是应急管理、应急响应、应急决策的基础。前言3、步骤–确定泄漏事件、泄漏源（第4章介绍了7种泄漏源模型）–建立源模型，对事故释放、释放速率、释放量等进行描述。–应用扩散模型估算下方向有毒有害物质的浓度，并由此根据一些准则来评估释放的后果及影响。5.1扩散影响参数•扩散模型用于描述事故释放后由释放源开始向其他地点及大范围区域的大气输送过程。释放发生后，大气中的毒害物质在风作用下以烟羽方式（连续源，图5-1）、或云团方式（瞬时源，图5-2）带走。1.风图5-1物质连续泄漏形成烟羽物质连续泄漏形成烟羽连续泄漏处形状与风速有关、与气云性质有关、地形，释放源位置及性质有关烟羽：通过与新鲜空气混合而消散风1.风图5-2物质瞬时泄漏形成烟团物质瞬时泄漏形成烟团物质瞬时泄漏生成烟团t1时的烟团t2时的烟团风三个烟团表面浓度相等：c环烟团向下风移动并通过与新鲜空气混合而消散5.1扩散影响参数多个因素影响着毒害物质在大气中的扩散：风速；大气稳定度；地形条件（建筑物、地面类型、地面构筑物）；释放源离地面的高度；释放的初始动量和浮力。2.大气稳定度•大气稳定度与空气的垂直混合有关。白天，空气温度随高度增加而下降，地面处空气受热密度变小（蒸腾），向上运动；夜晚，空气温度随高度增加下降不多，垂直运动较少。此外，由于土壤、地表导热系数大，夜间无太阳辐射，地面附近空气温度要降低，空气密度变大，垂直方向运动较小。图5-3，昼间和夜间空气温度随高度的变化，温度梯度影响空气的垂直运动1119753500400300200100-10高度M温度摄氏度夜间昼间2.大气稳定度稳定度划分：不稳定、中性和稳定划分标准：对地面加热速度与地面散热速度相对快慢•1加热速度＞地面散热速度地面附近的空气温度比高处的空气温度高，地表附近空气的密度小，上层空气密度大，密度小的空气在这种浮力作用下上升，导致大气不稳定。[晴天上午9、10点后，肉眼会观测到地表升腾；春秋早晨水雾消散]。F浮＞F重•2加热速度=散热速度。热量对大气扰动很小，但很难长久保持。F浮=F重•3加热速度＜散热速度。地面附近的温度比高处空气的温度低，地表附近空气密度大于高处空气的密度。F浮＜F重。重力影响抑制了大气机械湍流。3.地面条件•地面条件通过对大气施加曳力（摩擦力）改变风速分布及大小。•平坦的地面、水面等曳力小。•高层建筑、密集建筑群等曳力、起作用的范围均比较大。图5-3地面情况对垂直风速梯度的影响城市光滑水面郊区5004003002001000风速梯度高度M地面条件：影响地表的机械混合和随高度变化的风速4.释放高度•随着泄漏高度的增加，烟羽达到地面的时间和距离均增加；此外，烟羽到达地面时的浓度也会更低，原因是烟羽被更多的稀释。图5-5泄漏高度增加地面浓度降低对地面浓度的影响很大，释放高度越大，地面浓度越低。L2L1H1H2。5.释放物质的浮力和动量释放物质的浮力和动量改变了泄漏的有效高度；释放物质的物化性质决定了扩散影响的区域和程度。内部浮力占支配地位转变区周围环境湍流占支配地位初始加速和稀释风源图5-6泄漏物质的初始动量和浮力影响烟羽的特性,,D•喷射气体动量将气体带到释放点上方高处，使有效释放高度更高。•气体喷出后气体物化性质（ρ，μ，D）等导致气体是（上浮还是下沉）及其快慢。•经过一定时间在空气里传播后，释放的气体被充分稀释混合，与局部大气环境物化性质一致，变为中性浮力气体。5.释放物质的浮力和动量释放物质的浮力和动量改变了泄漏的有效高度；释放物质的物化性质决定了扩散影响的区域和程度。喷泉5.2中性浮力扩散模型•中性浮力扩散模型用来估算释放后所释放出的气体与空气混合，导致混合气体具有中性浮力后下风向各处的浓度；适用于低浓度气体。•经常有两种中性浮力蒸气云扩散模型：烟羽（连续源释放）和烟团模型（瞬时释放）。烟团模型可用来描述烟羽；对涉及动态烟羽的研究（如风向变化），须使用烟团模型。（5-1）•考虑固定质量Qm的物质瞬时泄漏到无限膨胀扩张的空气中，坐标系固定在释放源处。假设不发生反应，或不存在分子扩散，释放所导致的物质的浓度C可由式(5-1)。•考虑到湍流的影响，用平均值和随机量来代替速度；5.2中性浮力扩散模型5.2中性浮力扩散模型jtjjjccucKtxxx（5-7）Kt：湍流扩散系数，取决于大气湍流程度、地形条件等对式（5-7），可给定适当的初始条件（t=0的条件）和边界条件，即可对各种情况进行求解。释放源0xzy风向，风速u图5-7：风作用下连续点源泄露（烟羽）x：风向；y：横风向；z：垂直风向用于扩散模型的坐标系5.2中性浮力扩散模型图5-8风作用下烟团随风的移动x：风向；y：横风向；z：垂直风向5.2中性浮力扩散模型0xzyt2时的烟团t1时的烟团，以速度u向下风处移动风向，风速u不同时刻烟团表面浓度相等5.2中性浮力扩散模型：处理CASE-1：无风情况下、稳态、连续点源释放*jK=K●湍流扩散系数不变（所有方向上）●边界条件：rr,c0;222222ccc0xyzmx,y,z*222Q1C4Kx+y+zQm=常数●质量释放速率不变：ju0●无风：rr,c0;A：采用湍流扩散系数Kj的情形jtjjjccucKtxxxc0t●稳态：经坐标变换和积分：（5-16）rr,c0;CASE-2：无风烟团扩散*jK=K●湍流扩散系数不变（所有方向上）●边界条件：●一定量的物质瞬间释放*mQ*222mx,y,z,t3**2Qx+y+zCexp4Kt8(Kt)22222221cccccrK*txyzrrddddju0●无风：A：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理jtjjjccucKtxxx（5-20）x,y,zt0,C0CASE-3：无风情况下的非稳态连续点源释放(从开始释放到达到稳态前)Qm=常数ju0*jK=K●质量释放速率不变：●湍流扩散系数不变（所有方向上）●无风：●初始条件：●边界条件：rr,c0;2222221ccccK*txyz222mx,y,z,t**x+y+zQCerfc4Kr2KtA：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理（5-22）（5-7）化简为：●湍流扩散系数不变（所有方向上）●质量释放速率不变：CASE-4：有风情况下的稳态连续点源释放Qm=常数*jK=K●风沿一个方向如x方向：●边界条件：rr,c0;jxuuu常数222*222uccccKxxyz222m(x,y,z)**222QuCexpx+y+zx2K4Kx+y+z沿烟羽中心线yz0mx,0,0*QC4KxA：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理（5-27）（5-24）任意点处平均浓度：ju0rr,c0;constQm*CASE-5：无风时的烟团,湍流扩散系数是各向异性●烟团释放，●湍流扩散系数不同●边界条件：●无风：xyzK,K,K222xyz222ccccKKKtxyzzyxzyxtzyxKzKyKxtKKKtQmC2222/3,,,41exp8*A：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理jtjjjccucKtxxx（5-7）（5-29）化简为：CASE-6：有风情况下稳态连续点源释放，湍流扩散系数各向异性Qm=常数●连续释放，●●风沿一个方向如x方向：●边界条件：rr,c0;jxuuu常数与CASE-4相同，差别在于扩散系数xyzKKK222xyz222ccccuKKKtxyz沿烟羽中心，yz0mxyzQC4xKK22,,exp44mxyzyzyzQuyzCxKKxKKA：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理（5-31）（5-32）化简为：CASE-7：有风下的烟团●烟团释放，●各方向湍流扩散系数不同，●风沿x方向：●边界条件：与CASE-5相同，但有风rr,c0;jxuuu常数xyzKKKConstQm*222mx,y,z,t32xyzxyzx-utQ1yzCexp4tKKK8tKKKA：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理（5-33）通过坐标移动，经处理后可以得到：CASE-8：释放源在地面上的无风时的烟团与CASE-5相同，但释放源在地面。释放源在地面代表不能渗透的条件，因此浓度是CASE-5的两倍222mx,y,z,t32xyzxyzQ1xyzCexp4tKKK4tKKK（CASE-5）222mx,y,z,t32xyzxyzQ1xyzCexp4tKKK8tKKKA：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理（5-34）（5-29）CASE-9：释放源在地面上的稳态烟羽与CASE-6相同，但释放源位于地面，如图5-9所示。地面不能渗透。结果是浓度应是CASE-6的浓度的2倍22mx,y,zyzyzQuyzCexp4xKK2xKK22mx,y,zyzyzQuyzCexp4xKK4xKK（CASE-6）A：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理（5-35）（5-31）CASE-10：连续的稳态源。释放源在地面上方Qm=常数●连续释放，●●风沿一个方向x方向：●边界条件：rr,c0;jxuuu常数xyzKKK对此种情形，地面起着距源H处的不能渗透的边界作用。A：采用湍流扩散系数Kj的情形5.2中性浮力扩散模型：处理222mx,y,zyzzyzQuyuuCexpexpzHrexpzHr4Kx4Kx4Kx4xKK若rH0，简化为释放源在地面上的情况（5-36）5.2中性浮力扩散模型：处理B：简便有效的处理方法-Pasquill-Gifford模型对于上面10种情形都依赖于K。一般K随位