第4章20大气污染扩散估算模式s

第四章大气污染物扩散模式1、湍流扩散的基本理论2、高斯扩散模式3、污染物浓度的估算方法4、特殊气象条件下的扩散模式5、城市及山区的扩散模式6、区域大气环境质量模型7、烟囱高度设计8、厂址选择第一节湍流扩散的基本理论扩散的要素风：平流输送为主，风大则湍流大湍流：扩散比分子扩散快105～106倍风、湍流是决定污染物在大气中稀释扩散的最直接因素。湍流的基本概念湍流——大气的无规则运动风速的脉动（上、下）风向的摆动（左、右）起因与两种形式热力：温度垂直分布不均（不稳定）机械：垂直方向风速分布不均匀及地面粗糙度湍流扩散理论主要阐述湍流与烟流传播及湍流与物质浓度衰减的关系1.梯度输送理论类比于分子扩散，污染物的扩散速率与负浓度梯度成正比2.湍流统计理论泰勒－图4-1，正态分布萨顿实用模式高斯模式（应用最为广泛）第二节高斯扩散模式坐标系右手坐标系（食指—x轴；中指—y轴；拇指—z轴），原点：为无界点源或地面源的排放点，或者高架源排放点在地面上的投影点；x为主风向；y为横风向；z为垂直向高斯模式的四点假设a．污染物浓度在y、z风向上分布为正态分布b．全部高度风速均匀稳定c．源强是连续均匀稳定的d．扩散中污染物是守恒的（不考虑转化）一、高斯模式的有关假定高斯扩散模式高斯扩散模式的坐标系二、无界空间连续点源扩散模式2222(,,)exp[()]222πyzyzqyzcxyzuū—平均风速，m/s；q—源强，g/s；σy—侧向扩散参数，污染物在y方向分布的标准偏差，m；σz—竖向扩散参数，污染物在z方向分布的标准偏差，m；三、高架连续点源扩散模式高架源须考虑到地面对扩散的影响。根据假设④可认为地面就象镜子一样对污染物起全反射作用，按全反射原理，可用“像源法”处理——把P点污染物浓度看成为两部分（实源和像源）作用之和。建立三个坐标系：1、以实源在地面的投影点为原点；P点坐标为（x，y，z）；2、以实源为原点；3、以像源为原点。（1）实源贡献：P点在以实源为原点的坐标系中的垂直坐标为（z-H）。不考虑地面的影响，实源在P点形成的污染物浓度为：2222()(,,,)exp[()]222πyyyzqyzHcxyzHu实源的贡献2222()(,,,)exp[()]222πyyyzqyzHcxyzHu实源的贡献ρ1（2）像源贡献：P点在以像源为原点的坐标系中的垂直坐标为（z+H），像源在P点形成的污染物浓度为：2222()(,,,)exp[()]222πyzyzqyzHcxyzHu像源的贡献2222()(,,,)exp[()]222πyzyzqyzHcxyzHu像源的贡献ρ2222222()()(,,,)exp(){exp[]exp[]}2222πyyzyzqyzHzHcxyzHu实际浓度222222()()(,,,)exp(){exp[]exp[]}2222πyyzyzqyzHzHcxyzHu实际浓度ρ（x,y,z,H）(3)高斯烟流的形态2222(,,)exp[()]222πyzyzqyzcxyzu高架连续点源扩散模式2222(,,0,)exp()exp()22πyzyzqyHcxyHu地面浓度模式：取z＝0代入上式，得2222(,,0,)exp()exp()22πyzyzqyHcxyHu地面浓度模式：取z＝0代入上式，得22(,0,0,)exp()2πzyzqHcxHu地面轴线浓度模式：再取y=0代入上式22(,0,0,)exp()2πzyzqHcxHu地面轴线浓度模式：再取y=0代入上式22(,0,0,)exp()2πzyzqHcxHuyz上式，x增大，则、增大，第一项减小，第二项增大，必然在某x处有最大值地面最大浓度模式：考虑地面轴线浓度模式22(,0,0,)exp()2πzyzqHcxHuyz上式，x增大，则、增大，第一项减小，第二项增大，必然在某x处有最大值yz上式，x增大，则、增大，第一项减小，第二项增大，必然在某x处有最大值地面最大浓度模式：考虑地面轴线浓度模式高架连续点源扩散模式yzconstd(,0,0,)0dzcxHmax22πzyqcuHemax|2czxxH地面最大浓度模式（续）：设（实际中成立）由此求得yzconstd(,0,0,)0dzcxHmax22πzyqcuHemax|2czxxH地面最大浓度模式（续）：设（实际中成立）由此求得2222(,,,0)exp[()]22πyzyzqyzcxyzu地面源高斯模式（令H＝0）：相当于无界源的2倍（镜像垂直于地面，源强加倍）2222(,,,0)exp[()]22πyzyzqyzcxyzu地面源高斯模式（令H＝0）：相当于无界源的2倍（镜像垂直于地面，源强加倍）kzy设zyk则)2exp(A)2exp(Q222222zzzzHHku＝上式变为：’’‘’‘）（）（）（＝）（）（＝22222222222222212)2exp(1)2exp(1A)2exp(1)2exp(1AzzzzzzzzzzHHHHH令ρ’=0211022Hz＝2Hz＝22max2QQ)1exp(QHeueuuyzzyzzy＝例题：一工厂在源高H=60m处以80g/s的速度排放SO2，风速为6m/s，在下风向距离600m处，扩散系数分别取σy=30m，σz=20m。计算烟流中心线上SO2的浓度；中心线以左60m、以下20m处SO2的浓度。颗粒物扩散模式粒径小于15μm的颗粒物可按气体扩散计算大于15μm的颗粒物：倾斜烟流模式地面反射系数2222(1)(/)(,,0,)exp()exp[]222πtyzyzaqyHvxucxyHu2pp18tdgv第三节污染物浓度的估算q源强计算或实测平均风速多年的风速资料H有效烟囱高度、扩散参数uyz烟云抬升的原因有两个：①是烟囱出口处的烟流具有一初始动量（使它们继续垂直上升）；②是因烟流温度高于环境温度产生的静浮力。这两种动力引起的烟气浮力运动称烟云抬升，烟云抬升有利于降低地面的污染物浓度。1.烟气抬升高度的计算初始动量：速度、内径浮力：烟温度烟气抬升sHHHsHH――烟囱几何高度――抬升高度有效源高sHHHsHH――烟囱几何高度――抬升高度有效源高烟气抬升高度的计算（1）Holland公式：当大气稳定度为中性，计算烟气抬升高度时，经常使用霍兰徳公式：Holland公式比较保守，特别在烟囱高、热释放率比较强的情况下。在实际计算中，不稳定条件(A、B稳定度)，ΔH需增加10％～20％；稳定条件(D、E、F稳定度)，ΔH需减少10％～20％。3ssaHs1(1.52.7)(1.59.610)svDTTHDvDQTuukW抬升高度计算式H11/32/3sH11/32/3sH21000kW10=0.36210=1.55当时sQxHHQxuxHHQHuH11/31/3H3/52/5Hs6/53/53/5Hs21000kW3*=0.3623*=0.332*=0.33当时QxxHQxuxxHQHxQHu（2）Briggs公式：适用不稳定及中性大气条件不稳定或中性大气下，布里格斯公式用来确定不同的热释放率和下风向距离条件下的烟气抬升高度：（3）我国“环境影响评价技术导则——大气环境”(HJ/T12.2-93)中的公式——在没有特别要求时，应优先使用国家标准规定的方法。12Hsa1nn0HsHaVasHH121sH12100kW()35K=0.351700kW2100kW1700=()4002(1.50.01)0.04=sQTTHnQHuTQPQTTTTQQHHHHvDQHu(1)当和时(2)当时HHsH1/43/8aH8(1700)1700kW35K2(1.50.01)=10m1.5m/sd=5.5(0.0098)dQuQTvDQHuTHQz(3)当或时(4)当高处的年平均风速小于或等于时例：某市远郊区电厂烟囱高160m，烟囱排出口内径5m，排烟速度12m/s。烟气温度135℃，周围大气温度15℃。大气稳定度C级，源高处风速612m/s。试分别用霍兰德、布里格斯、国家标准公式计算烟气抬升高度（假设下风向距离x=2km）大气扩散参数（σy，σz）的确定P－G曲线法帕斯奎尔在1961年推荐一种仅需要常规气象观测资料就能估算σy，σz的方法，吉福德(Gifford)进一步将它制成应用更方便的图表。应用观测到的风速、云量、云状和日照等天气资料，将大气扩散稀释能力分为6个等级：A—极不稳定，B—不稳定，C—弱不稳定，D—中性，E—弱稳定，F—稳定。若稳定级别为A～B，则表示按A、B级的数据内插。该法的要点：首先根据帕斯奎尔划分大气稳定度的方法来确定大气稳定度级别；然后从图4-4和图4-5中查得（或表4-4用内插法求出）对应的扩散参数σy和σz；最后将σy、σz代入前面介绍的一系列扩散模式中，就可估计出各种情况下的浓度值。扩散参数是表征湍流扩散剧烈程度的物理量，是影响污染物浓度的重要参数。扩散参数的确定－P－G曲线法P－G曲线的应用根据常规资料确定稳定度级别扩散参数的确定－P－G曲线法P－G曲线的应用利用扩散曲线确定和yz扩散参数的确定－P－G曲线法P－G曲线的应用地面最大浓度估算Hmax|2czxxzH~zxmaxcx~yxmaxC由和由曲线（图4-5）反查出由曲线（图4-4）查由式（4-10）求出yHmax|2czxxzH~zxmaxcx~yxmaxC由和由曲线（图4-5）反查出由曲线（图4-4）查由式（4-10）求出y扩散参数的确定－中国国家标准规定的方法我国在修订P-T法基础上产生了国家标准法(HJ/T2.2-93)。太阳高度角（式4-29，地理纬度，倾角）辐射等级稳定度云量（加地面风速）该方法的技术路线是：根据时间、地理位置确定日倾角、太阳高度角，利用天气条件确定辐射等级，然后利用辐射等级和风速确定大气稳定度，最后查扩散参数幂函数表，确定扩散参数。扩散参数的确定－中国国家标准规定的方法扩散参数的选取扩散参数的表达式为（取样时间0.5h，按表4-8查算）平原地区和城市远郊区，D、E、F向不稳定方向提半级工业区和城市中心区，C提至B级，D、E、F向不稳定方向提一级丘陵山区的农村或城市，同工业区取样时间大于0.5h，不变，1221,aayzxx2121()qyyz例：某冶炼厂烟囱高150m，烟气抬升高度75m，SO2排放量1000g/s。估算风速3m/s，大气稳定度C级时地面最大浓度是多少？发生在什么位置？（分别用P-G法和国家标准方法计算）第二步：确定出现地面最大浓度的下风向距离。第一步：确定出现地面最大浓度的Z向扩散参数。第三步：确定出现地面最大浓度的y向扩散参数。第四步：计算地面最大浓度。第四节特殊气象条件下的扩散模式在整层大气都具有同一稳定度（即温度层结构均一，实际中难以实现）、平坦地形的条件下应用高斯模式计算污染物浓度。如果整个大气层不均匀，污染物扩散所涉及的温度层结不止一个，或者地表粗糙度高，地势起伏大就需要特殊处理。下面讨论两种特殊情况：封闭型扩散和熏烟型扩散。一、封闭型扩散模式——有上部逆温层的扩散模式所谓封闭型扩散就是指在上部存在逆温层的气象条件下，污染物受到上部逆温层限制，只能在地面和逆温层之间扩散的情况第四节特殊气