污染物扩散计算模式汇总.doc

... 大气稳定度分级 常用的大气稳定度分类方法有帕斯奎尔(Pasquill)法和国标原子能机构IAEA推荐的方法。这里介绍的是中国现有法规中推荐的修订帕斯奎尔分类法（简记P·S），分为强不稳定、不稳定、弱不稳定、中性、较稳定和稳定六级。它们分别表示为A、B、C、D、E、F。确定等级时首先计算出太阳高度角按表B1查出太阳辐射等级数，再由太阳辐射等级数与地面风速按表B2查找稳定等级。 注：云量（全天空十分制）观测规则与现国家气象局编定的《地面气象观测规范》相同。 注：地面风速（m/s）系指距地面10m高度处10min平均风速，如使用气象台（站）资料，其观测规则与国家气象局编定的《地面气象观测规范》相同。 太阳高度角ho使用下式计算： .....................(B1) 式中： ho----太阳高度角，deg ； ----当地纬度，deg.； λ----当地经度；deg ； t----进行观测时的北京时间； σ----太阳倾角，deg，可按下式计算： 式中： θo----360dn/365，deg； dn---- 一年中日期序数，0、1、2、······364。 地形对烟羽的影响 此前的扩散模式都假设地面是完全平整的（烟囱底部是一个无限大的水平面，其高程为0），因此在扩散过程中烟羽的中心线可保持水平不变。但如果在预测点（x,y,z）处，地面有一定的高程hT(0hTz)，则在对（x,y,z）式应用以上模式时，应对有效烟羽高度进行一些修正。 假定烟羽路径始终与起伏的地形保持平行，或者假设烟羽轴线保持固定的海拔高度，并与高于烟羽的地形相交，都是不正确的，实际情况应该是介于上述二者之间。具体的修正方法如下。 （1）中性和不稳定天气条件 令：hT为凸出的地形高度；He为烟轴高度（即有效高度）；T为烟轴高度修正系数（或地形系数），修正后的烟囱有效高度应该是THe oT则应按下式取值： （2）稳定天气条件（D-E、E、F） 在稳定天气条件下，当烟羽逼近孤立山体时，烟羽以临界高度Hc为界分成两部分，临界高度以上的烟羽有足够的动能爬越山体，而临界高度以下的烟羽，只能被迫绕着山体过去。临界高度H c 可由下式确定： Hc=Hm--u[θ/(gdθ/dz)]1/2 式中 Hm----孤立山体高度，m； Hc----临界高度，m； θ----z高度处大气位温，K； dθ/dz----z高度处位温梯度，K/m； u----平均风速，m/s； g----重力加速度，m/s2。 例如，θ=300K，dθ/dz=0.01K/m，u=2m/s，Hm=200m；则Hc=Hm-111=89m，烟囱有效高度大于89m时，烟羽将有足够的动能爬越山体。 对于山体高度Hm已定的情况，大气越稳定，则Hc越小。所以一般只需计算在F稳定度下的Hc，如果烟羽有效高度HeHc(F)，则可认为烟羽能够爬越山体。 点源扩散模式 持续排放源 有风模式（U10≥1.5m/s） 1.自由空间中的连续点源 实际上绝大多数污染源都是连续的，对于连续排放源，可理解为在时间上依次连续释放无穷多个烟团。因此，连续排放源的扩散模式可以通过将瞬时单烟团模式对to从—∞到t积分后求得。以烟团初始空间坐标为原点，下风方为x轴，烟羽轴线与x轴一直保持重合，都是x的函数，将对to的积分变换为对（x—uT）/σx的积分，可得最基本的烟羽扩散模式： 适用条件为：自由空间；风速要比较大（u10≥1.5m/s）；当大气不稳定状态时，可能带来一定的误差。 2.地面反射 用像源法，假想地平线为一镜面，在其下方有一与真实源完全对称的虚源，则这两个源叠加后的效果和真实源考虑到地面反射的结果是等价的。以烟囱地面位置的中心点为坐标原点，在考虑到地面反射后，污染源下风方任一点小于24小时取样时间的污染物浓度C（x，y，z）由下式给出： z=0时的地面浓度C（x,y,0），可简化为； 下风方X轴线上（y=0）的地面浓度C（x,0,0）为： 对于较低的排放源（例如He＜50m，具体限值由地面粗糙度、混合层高度等因素决定），一般可直接应以上式子计算。 3.混合层反射 对于高架源，当超过一定的下风距离时，需对烟羽在混合层顶的反射进行修正。同考虑地面反射类似，用像源法修正后，污染源下风方任一点小于24小时取样时间的污染物浓度C（x,y,z）可表示为： 式中 h——混合层高度； k——反射次数，一、二级项目取k=4已足够。 4.侧面反射 详见狭长山谷扩散模式。 小风静风模式（U101.5m/s） 小风静风时，污染物地面浓度C(x,y,0)可用下式计算： 式中η和G按下式计算： ， ， 和分别是横向和铅直向扩散参数的回归系数(σy=σz=T,σz=T)