风向变化下核污染扩散模型.docVIP

风向变化下核污染扩散模型 　　摘要本文基于高斯烟团原理建立了核泄漏中核素气体的扩散模型，通过分析了干沉积、核衰变对核污染的影响，对原高斯模型进行了修正，并进一步运用坐标变换解决了在风向变化的条件下气体扩散的问题。 　　关键词高斯扩散模型；核素气体；坐标变换；核衰变；干沉积 　　0引言 　　核反应堆失控导致夹杂着大量放射性物质的水蒸气排出，在热力和动力作用下上升，与大气混台，形成含有大量放射性核素的云团。如果能选择适当的气体扩散模型对核素气体浓度情况进行有效的预测，就可以为突发的核事故预警、为相关核辐射检测部门提供重要的数据支持，协助人们进行高效的物资调拨和人员撤离，最大程度地降低核事故所造成的伤害和损失。 　　1高斯浓度模型建立 　　本文假设污染物浓度符合正态分布，建立高斯浓度分布模型。由于核泄漏是突发事件，核反应堆一次性释放出大量的含有放射行为物质的气体，因此本文采用了烟团模式进行地面核素气体浓度的计算，设源强为：Q（kg/s），风速为：u（m/s），下风向、横风向和垂直风向的扩散参数（m）分别为：σx，σy，σz（计算公式由文献[4]给出），其取值与泄漏时间、大气稳定度和地面有效粗糙度等因素有关[4]，则某时刻下空间点（x,y,z）处的气体浓度为： 　　2高斯模型的修正 　　因为干沉积和核素衰变是影响污染粒子在大气中的衰减的重要因素，在核污染事故中碘-131的影响较大，故引入其浓度随时间的变化关系：C=C0*exp（-0.00358t） 　　干沉积量：D（计算公式由文献[2]给出） 　　则某时刻下空间点（x,y,z）处的气体浓度为： 　　（1） 　　（2） 　　3坐标变换 　　3.1坐标变换公式I 　　为了满足高斯模式中横轴与下风向保持一致的假设，我们需要在基准坐标系的基础上进行一次坐标变换。设核素气体泄漏源在IOJ下的坐标为D=（i0,j0），以该点为原点建立新的笛卡儿坐标系XDY，风向为中X轴正方向。设风向与I轴正方向的夹角为θ，则核素气体扩散区域内任意一点在两个坐标系下的坐标（x，y）和（i，j）满足关系式： 　　设t=0时发生核泄漏事故，并且风向在气体扩散的过程中保持不变。根据公式（2）（3），则t时刻下基准坐标系IOJ中（i，j）处的核素气体浓度为： 　　3.2坐标变换公式Ⅱ 　　核素气体在扩散过程中风向并不是一成不变的，当风向发生变化以后，原XDY坐标系X轴与新风向不一致，于是高斯模式的前提条件未能满足。如果想继续使用烟团模式计算核素气体扩散的情况，需要将旧坐标系XDY进行平移和旋转至X’D’Y’，使得X’轴与新风向相同。设风向在t=t1时发生一次变化，此时气团中心位于（x1,y1）=（ut1，0），X与X’的夹角为α。 　　通常情况下风向发生变化的幅度不大,所以我们不妨设风向变化在瞬间完成。在t1时刻内坐标系X’D’Y’中：设在[0，t1]时间内风向保持不变风速为u，那么气团从点D’开始扩散后，经过时间t1，其中心将位于（x1，y1）。且核素气体扩散区域内任意一点在两个坐标系下的坐标（x1,y1）和（x’,y’）满足： 　　设为t=0时发生核泄漏事故，在t=t1时发生了第1次风向变化后，风向不再发生变化。由（2）、（3）、（4）、（5）得，在t时刻下基准坐标系IOJ（i，j）处的气体浓度为： 　　3.3坐标变换公式Ⅲ 　　设为t=0时发生核泄漏事故，在整个扩散过程中共发生了N次风向变化，发生时间依次为：t1，t2···tN。则可在[0，t1]，[t1，t2]，···，[tN-1。，tN]，[tN。，∞）中反复使用坐标变换（6）来计算基准坐标系IOJ中任一点处（i，j）处的核素气体浓度。 　　4结论 　　本文利用高斯烟团模型，并引入干沉积、核素衰变等影响因素，求出了新型的核素气体浓度公式，进一步利用坐标变换的原理给出了，在有风向变化下，核素气体浓度变化公式。本文建立的核素气体扩散模型更适合于实际问题，可以为突发的核泄漏事故提供相应的理论与数据指导，尽可能核事故的危害降到最低。 　　参考文献 　　[1]张梦瑶，崔晋川.适用于风向变化条件的毒气扩散模型研究[J].运筹与管理，2008，17(3). 　　[2]董希琳，赵智键.核电站核事故核素污染评估模式[J].火灾科学，1999，8(2). 　　[3]柳盛.可变风向的空气污染扩散模拟系统的实现[J].科技信息，2010(27). 　　[4]化工部化工劳动保护研究所.重要有毒物质泄漏扩散模型研究，1996(3). 　　[5]姚仁太.核电厂年均大气弥散因子估算中有关参数和修正方法的影响分析. - 3 -