有毒物质泄漏及中性浮力扩散模型

1、有毒物质泄漏及中性浮力扩散模型有毒物质泄漏及中性浮力扩散模型中国环境科学研究院大气环境研究所 陈义珍第五章第五章 有毒物质泄漏及中性浮力扩散模型有毒物质泄漏及中性浮力扩散模型 5.1 扩散影响参数 5.2 中性浮力扩散模型 5.3 模型应用举例5.1 扩散影响参数 扩散模型描述了有毒物质远离事故发生池，并遍及整个工厂和社会的空中运输过程。释放发生后，空气中的有毒物质被风以烟羽的形式，或烟团的形式带走。有毒物质的最大浓度是在释放发生处（可能不在地面上）。由于有毒物质与空气的湍流混合和扩散，其在下风向的浓度较低。 众多因素影响着有毒物质在大气中的扩散：风速大气稳定度地面条件（建筑物、水、树）释放处

2、距离地面的高度物质释放的初始动量和浮力5.2 中性浮力扩散模型 中性浮力扩散模型，用于估算释放发生后释放气体与空气的混合，并导致混合气云具有中性浮力后下风向各处的浓度。因此，这些模型适用于低密度的气体，特别是浓度为10-6量级的。 经常用到两种类型的中性浮力蒸气云扩散模型：烟羽模型和烟团模型。烟羽模型描述来自连续源释放物质的稳态浓度。烟团模型描述一定量的单一物质释放后的暂时浓度。对于烟羽模型，典型例子是气体自烟窗的连续释放。稳态烟羽在烟窗下风向形成。对于烟团模型，典型例子是由于贮罐的破裂，一定量的物质突然泄漏，形成一个巨大的蒸气云团，并渐渐远离破裂处。 烟团模型能用来描述烟羽；烟羽只不过是连续

3、释放的烟团。然而，如果稳态烟羽信息是所需的所有信息，那么建议使用烟羽模型，因为它比较容易使用。对于涉及动态烟羽的研究（例如，风向的变化对烟羽的影响），必须使用烟团模型。f.a.gifford,”use of routine meteorological observations for estimating atmospheric dispersion,” nuclear safety(1961)，2(4):47. f.a.gifford,”turbulent diffusion-typing schemes；a review，” nuclear safety(1976)，17(1):68.

4、强烈的日光照射是指英国盛夏正午期间的充足的阳光。弱的日光照射是指严冬时期类似的情况。夜间是指日落前1h至破晓后1h这一段时间。这些数值是作者自己补充进去的，以使该表格完整。对于白天或夜晚的多云情况以及日落前或日出后数小时的任何天气情况，不管风速有多大，都应该使用中等稳定的等级d。注：稳定度等级 a-极度不稳定；b-中度不稳定；c-轻微不稳定；d-中性稳定；e-轻微稳定；f-中度稳定情况1: 地面瞬时点源的烟囱，坐标系固定在释放点，风速u恒定，风向仅沿x方向2*22223/21 x-utyzx,y, ,exp22mxyzxyzqcz t (5-1) 地面浓度可令z=0，求得2*223/21 x-

5、utyx,y,0,exp22mxyxyzqct (5-2) 地面上沿x轴的浓度可令y=z=0，求得2*3/21 x-utx,y,0,exp22mxxyzqct (5-3) 气云中心坐标在（ut，0，0）处。该移动气云中心的浓度为*3/2ut,0,0,2mxyzqct (5-4) 站在固定点（x，y，z）处的个体，所接受的全部剂量dtid是浓度的时间积分0 x,y,x,y, ,dtidzcz t dt(5-5) 地面的全部剂量，可依照式（5-5）对式（5-2）进行积分得到。结果为*221yx, y,0exp2myzyqdtidu (5-6) 地面上沿x轴的全部剂量为*x , 0, 0myzqd

6、tidu (5-7) 通常情况下，需要用固定浓度定义气云边界。连接气云周围相等浓度的点的曲线称为等值线。对于指定的浓度 ,地面上的等值线通过用中心线浓度方程式（5-3）除以一般的地面浓度方程式（5-2）来确定。直接对y求解该方程*cx , 0, 0,2 lnx , y, 0,yctyct（5-8） 过程如下。情况2: 地面上的连续稳态源的烟囱，风向沿x轴，风速恒定为u22221yzx,y,exp2myzyzqczu （5-9） 地面浓度可令z=0求出21yx, y,0exp2myzyqcu （5-10） 下风向、沿烟羽中心线的浓度可令y=z=0求出x,0,0mxyqcu（5-11） 可使用类似

7、于情况1中使用的等值线求解过程来求得等值线。对于地面上的连续释放，最大浓度出现在释放处。情况3：位于地面hr高处的连续稳态源的烟羽，风向沿x轴，风速恒定为u2221y1 z-hr1 z+hrx,y,expexpexp2222myzyzzqczu （5-12） 地面浓度可令z=0求出221y1hrx,y,0exp22myzyzqcu （5-13） 地面中心线浓度可令y=z=0求得21hrx,0,0exp2myzzqcu （5-14） 地面上沿x轴的最大浓度 maxc由下式求得 2m axr2hmzyqceu（5-15） 下风向地面上的最大浓度出现的位置，可由下式求得rh2z（5-16） 求解最大

8、浓度和下风向距离的过程是：使用式（5-16）确定距离，然后使用式（5-15）计算最大浓度。 情况4：位于地面hr高处的瞬时点源的烟团，坐标系位于地面并随烟团移动222*3/21y1 z-hr1 z+hrx,y, ,expexpexp(2 )222mxyzyzzqcz t （5-17） 时间相关性通过扩散系数来完成，因为，随着烟团从释放处向下风向运动，它们的值也发生变化，如果没有风（u=0），式（5-17）预测的结果是不正确的。 在地面，即z=0，浓度可通过下述方程计算 22*r3/2h1y1x,y,0,exp222myzxyzqct （5-18） 沿地面中心线的浓度，可令y=z=0求出2*r3

9、/2h1x,0,0,exp22mzxyzqct （5-19） 通过应用式(5-5)(5-18)，即得到地面上的全部剂量。结果为22*rh1y1x, y,0exp22myzyzqdtidu （5-20） 情况5：位于地面hr高处的瞬时点源的烟团，坐标系位于地面的释放点21x-utx,y, ,-exp2xcz t使用移动坐标系的烟团方程式（5 54）式（5 56）（5-21） 式中，t为自从烟团释放后的时间。最坏事件情形 对于烟羽，最大浓度通常是在释放点处。如果释放是在高于地平面的地方发生，那么地面上的最大浓度出现在释放处的下风向上的某一点。 对于烟团，最大浓度通常在烟团的中心。对于释放发生在高于

10、地平面的地方，烟团中心将平行于地面移动，并且地面上的最大直接位于烟团中心的下方。对于烟团等值线，随着烟团向下风向的移动，等值线将接近于圆形。等值线的直径一开始随着烟团向下风向的移动而增加，然后达到最大，最后将逐渐减小。 如果天气条件未知或不确定，那么可进行某些假设，以得到一个最坏情形的结果，即估算一个最大浓度。 pasquill-gifford扩散方程中的天气条件可通过扩散系数和风速予以考虑。通过观察估算浓度用pasquill-gifford扩散方程，很明显扩散系数和风速在分母上。因此，通过选择导致最小值的扩散系数和风速的天气条件和风速，可使估算的浓度最大。我们知道，f类稳定度等级可以产生最小

11、的扩散系数。很明显，风速不能够为零，所以必须选择一个有限值。epa人为，当风速小到1.5m/s时，f类稳定度等级能够存在。一些风险分析专家使用2m/s的风速。在计算中所使用的假设，必须清楚地予以说明。5.3模型应用举例【例5-1】 在某个阴天，一个有效高度为60m的烟囱，正以80g/s的速度排放二氧化硫。风速为6m/s。烟囱位于农村。请确定：（1）下风向500m处地面上二氧化硫的平均浓度。（2）下风向500m、横风向50m处地面上二氧化硫的平均浓度。（3）径直风向地面上最大平均浓度的位置和数值。解 （1）属于连续排放。径直下风向地面浓度由式（5-14）给出21,0,0exp2rxyzhqmcx

12、u （5-14） 由表5-1可知，大气稳定度等级为d。扩散系数可由表5-2得到。利用表5-2y=0.08x(1+0.0001x）-1/2=0.08500m(1+0.0001500m) -1/2 =39.0mz=0.06 x(1+0.0015x) -1/2=0.06500m(1+0.0015500m) -1/2 =22.7m代入式（5-14），得到280 /160m500m,0,0exp3.14 39.0mmm s2 22.7mg sc =1.4510-4g/m3（2）横风向50m处的平均浓度可通过式（5-13）及设y=50得到。直接应用来自（1）的结果21y500m,50m,0500m,0,0

13、 exp2ycc2-431 50m1.45 10 g/mexp2 39m（） = = -536.37 10 g/m（3）最大浓度的位置由式（5-16）确定zr60m=42.4m22hzy由pasquill-gifford烟羽模型扩散系数图可知，对于稳定度等级d，在下风向约1200m处=42.4m。由表5-2可知，=88m。最大浓度由式（5-15）确定2maxr22=e2 80 /42.4m =2.72 3.14m s60m88mzyqmcuhg s （）（5-15） =4.1810-4g/m3 【例5-2】在精细化工过程中使用氯。源模型研究表明，对于某个特殊的事故场景，1.0kg的氯将瞬时释放

14、出来。释放发生在地面。居民区距离释放源500m。请确定：（1）气云中心到达居民区所需的时间。假设风速为2m/s。（2）居民区氯的最大浓度。将结果与氯的erpg-1值1.010-6进行比较。怎样的大气稳定度和风速会产生最大浓度。（3）确定气云最大浓度低于erpg-1值所必须移动的距离。使用（2）中的条件。（4）根据erpg-1确定下风向地面5km处气云的大小。假设在所用的情况下氯气云释放都是中性浮力（这可能是不正确的假设）。解 （1）对于下风向500m距离和2m/s的风速，气云中心到达居民区所需要的时间为x500mt=250s=4.2 minu2m / s对于紧急警告，该时间很短。（2）最大浓度

15、出现在释放源下风向的气云中心。浓度由式（5-4）给出 *3/2,0,0,2mxyzqcutt （5-4） 选择稳定度条件，使式（5-4）中的c最大。这需要扩散系数的值最小。任意一个扩散系数的最小值都发生在f稳定度等级。这种情况对应着有薄云和少云的夜晚条件，同时风速小于3m/s。烟团中的最大浓度也出现在居民区中距离释放源最近的点，这发生在500m处。因此0.890.890.020.02 (500m=5.0mxyx）0.610.61=0.05x=0.05500m=2.2mz （）由式（5-4）有-3333/21.0kg=2.31 10 kg/m =2310mg/m2.m. mc 利用式(2-6)，

16、将其单位转变为10-6。假设压力为1atm，温度为298k，浓度为79810-6。这比高很多。如果是在户外和泄漏源的下风向，那么，工厂内的人和任何一个处在紧邻居民区的人都将过量暴露。（3）由表2-7可知，1.010-6的erpg-1值为3.0mg/m3或3.010-6kg/m3。气云中心的浓度由式（5-4）给出。代入已知数据，可以得到-633/221.0kg3.010 kg/m =2.yz243=4.24 10 myz 使用表5-3中提供的方程。求解下风向的距离。因此对于f稳定度等级： 20.89 20.6143= 0.020.05x=4.24 10 myzx （） （）通过试差法求解x，得到x=8.0km。 （4）下风向中心线的浓度由式（5-3）给出2*3/21x-utx,y,0,exp22mxxyzqct （5-3） 烟囱中心到达下风向5km处，所需时间为x5000 mt= 2500su2 m / s在下风向5km处，假设大气为f稳定度等级，计算0.890.610.0239.2m=0.05x=9.0mxyz