空气污染学-第四章-理想条件下空气污染物散布的模式处理课件

第四章第四章 理想条件下空气污染物理想条件下空气污染物散布的模式处理散布的模式处理2主要内容主要内容l连续点源高斯扩散计算公式连续点源高斯扩散计算公式l连续线、面源和体源扩散计算公式连续线、面源和体源扩散计算公式l大气扩散参数大气扩散参数l烟流抬升高度烟流抬升高度l非扩散过程的处理非扩散过程的处理l特殊条件下的扩散特殊条件下的扩散3l欧拉扩散方程假设K为常数即斐克扩散，可得到正态分布形式的解l从统计理论出发，在平稳、均匀湍流的假定下，也可以证明粒子扩散位移的概率分布符合正态分布形式l对连续点源发出的烟流的大量试验研究和观测事实说明，尤其是对于平均烟流的情形，其浓度分布是符合正态也称高斯分布的44.1 连续点源高斯扩散公式连续点源高斯扩散公式 5一一 无界情形公式及物理意义无界情形公式及物理意义l湍流均匀定常，设源位于无界空间，取湍流均匀定常，设源位于无界空间，取X轴与平均风向一致，那么污染物浓度轴与平均风向一致，那么污染物浓度在在y和和z方向符合高斯分布，可得：方向符合高斯分布，可得：6物理意义物理意义lQ：源强，点、面、线、体源，影响直接、明显，影响大源强，点、面、线、体源，影响直接、明显，影响大l大气稀释因子：大气稀释因子：代表了不同气象条件和地形条件下物质散布的程度及其随空代表了不同气象条件和地形条件下物质散布的程度及其随空间距离的变化间距离的变化l正态分布形式项：在正态分布情况下，分布形式的影响不敏正态分布形式项：在正态分布情况下，分布形式的影响不敏感感7l镜像全反射镜像全反射-像源法像源法l实源：实源：l像源像源：二二 有界情形掌握有界情形掌握8 8总奉献：总奉献：源强源强平均风速平均风速扩散参数扩散参数有效源高有效源高9烟流有效源高：烟流有效源高：H=hs+h 归一化浓度：归一化浓度：10三三 地面源地面源l取取H0，有界情形是无界情形地面浓度两倍有界情形是无界情形地面浓度两倍11四四 地面浓度和地面最大浓度地面浓度和地面最大浓度1 地面浓度地面浓度令令 z=0，可得高架源的地面浓度，可得高架源的地面浓度122 地面轴线浓度地面轴线浓度令令y=0,z=0 可得高架源地面轴线浓度可得高架源地面轴线浓度13高斯烟流的浓度分布高斯烟流的浓度分布14源高和稳定度的影响源高和稳定度的影响151 与与 之比为常数之比为常数3 地面最大浓度的估算地面最大浓度的估算假设稳定度不变，增加有效源高假设稳定度不变，增加有效源高H H，那么会在，那么会在更远处出现到达最大浓度更远处出现到达最大浓度qm qm 所需的扩散参数所需的扩散参数 。162 假设 与 之比是变化的17一一 线源扩散公式线源扩散公式线源定义为呈线状分布的污染物排放线源定义为呈线状分布的污染物排放源。如繁忙的公路和城市的街道通源。如繁忙的公路和城市的街道通常被看作是线源常被看作是线源连续线源等价于连续点源沿着线源长连续线源等价于连续点源沿着线源长度范围的积分，其浓度场是线上无度范围的积分，其浓度场是线上无数点源浓度奉献之和数点源浓度奉献之和4.2 连续线、面、体源扩散公式连续线、面、体源扩散公式18l对于直线型的线源，可直接积分求出；对于很不对于直线型的线源，可直接积分求出；对于很不规那么的线源，只能用数值求和的方法解决规那么的线源，只能用数值求和的方法解决l点源计算一般取点源计算一般取x x轴与风向一致，线源计算时需考轴与风向一致，线源计算时需考虑风向与其交角以及线源的长度虑风向与其交角以及线源的长度191 无限长线源无限长线源 风向与其正交风向与其正交 为线源源强，为线源源强，mg/(sm)20一般一般 不适用不适用风向与线源成交角风向与线源成交角 时时21 风向与其平行，只有上风向有奉献，浓度与顺风位置无关。222 有限长线源有限长线源 设线源长度为范围为设线源长度为范围为 ，根据不，根据不同情况取积分有：同情况取积分有：无界情形有限线源：无界情形有限线源：有界情形高架线源：有界情形高架线源：23二二 面源扩散公式面源扩散公式l在水平方向呈面块状散布的污染物排放源在水平方向呈面块状散布的污染物排放源,称面源。如称面源。如散布很集中而高度低源强小的居民区污染源，可视为面散布很集中而高度低源强小的居民区污染源，可视为面源。源。l面源扩散公式，原那么上可由点源扩散公式沿面源扩散公式，原那么上可由点源扩散公式沿x x和和y y方向方向积分而得。积分而得。l假设面源源强为假设面源源强为 ，自整个上风方的半平面对，自整个上风方的半平面对x=0 x=0和和 y=0 y=0点造成的浓度奉献点造成的浓度奉献可分为两种：可分为两种：l 1.1.地面面源地面面源(即即H=0)H=0)情形情形l 2.2.近地层面源近地层面源(即即H )H )情形情形24 由点源沿由点源沿x x和和y y向积分给出，自上风向半平面向积分给出，自上风向半平面对对x=0 x=0，y=0y=0造成的浓度奉献造成的浓度奉献实际运用时，常处理积分并作源的编目和模实际运用时，常处理积分并作源的编目和模式化处理。式化处理。地面面源地面面源25l这里有一个面源如何合理处理的技术问题。经验这里有一个面源如何合理处理的技术问题。经验上采用上采用“虚点源法，此法大致分两步：虚点源法，此法大致分两步：1 1把大把大块面源划分为假设干较小的面源块情况下，可将块面源划分为假设干较小的面源块情况下，可将面源化为点源来处理，即将每一面源块又称面面源化为点源来处理，即将每一面源块又称面源单元简化为一个等效的点源源单元简化为一个等效的点源假定整个面源假定整个面源块的污染物排放集中在该点。这样就可用点源扩块的污染物排放集中在该点。这样就可用点源扩散公式来计算该面源所造成的污染物浓度。散公式来计算该面源所造成的污染物浓度。2 2由烟流半宽定义有由烟流半宽定义有26l此时，我们能在上游方向确定出虚拟点源的位置。即此时，我们能在上游方向确定出虚拟点源的位置。即在这一点上，使虚拟点源的扩散参数恰等于该面源的在这一点上，使虚拟点源的扩散参数恰等于该面源的初始扩散参数初始扩散参数 。这样，可得高架这样，可得高架源地面浓度为源地面浓度为27l在应用上式时在应用上式时,常采用经验方法给出初始扩散参数常采用经验方法给出初始扩散参数 ，如下图，如下图。假设。假设:图中面源方块边长为图中面源方块边长为a,虚拟点源出发虚拟点源出发的烟流在抵达面源中心位置时的横风向宽度的烟流在抵达面源中心位置时的横风向宽度2y。=a，横向初始扩散参数为横向初始扩散参数为28于是上式变成：于是上式变成：29 面积较小的面积较小的虚点源法虚点源法 30三三 体源扩散公式自学体源扩散公式自学l与面源类似与面源类似31重点重点l理解记忆掌握点源高斯扩散公式理解记忆掌握点源高斯扩散公式l理解掌握线、面源高斯公式理解掌握线、面源高斯公式32l早期大气扩散参数处理早期大气扩散参数处理l稳定度扩散级别与扩散曲线法稳定度扩散级别与扩散曲线法l扩散曲线讨论扩散曲线讨论l风向脉动与扩散函数法风向脉动与扩散函数法l扩散参数的研究现状扩散参数的研究现状4.3 4.3 大气扩散参数大气扩散参数33一一 早期的扩散参数模式早期的扩散参数模式 格雷厄姆萨顿，英国气象学家。1903年2月4日生于克温坎。毕业于威尔士大学、阿伯里斯威恩大学和牛津大学。19261928年在威尔士大学、阿伯里斯威恩大学任讲师，19281941年任助理教授。第二次世界大战期间从事国防科研工作。19421943年任英国国防部防化实验所所长。19431945年任坦克实验所所长。19451947年任英国雷达研究开展中心主任。19501955年任英国大气污染研究委员会主席。1951年任英国陆军部科学参谋。19521953年任英国皇家军事科学院教务长。1953年任皇家气象学会主席。19531956年任世界气象组织常务理事。19601966年任英国大地测量及地球物理学全国委员会主席。19651971年任英国自然环境研究委员会主席。在自然环境和气象研究方面取得了许多成果。曾获世界气象组织颁发的奖金。著作：?大气湍流?Atmos-pheric turbulence，1948；?微气象学?Micrometeoro-logy，195334具体步骤：具体步骤：1 找出泰勒公式中的拉格朗日相关系数的具体形式，找出泰勒公式中的拉格朗日相关系数的具体形式，即寻找它与某些可测气象参量的关系，代入泰勒即寻找它与某些可测气象参量的关系，代入泰勒公式求扩散参数。公式求扩散参数。2 将扩散参数代入根本高斯扩散，得到萨顿扩散公将扩散参数代入根本高斯扩散，得到萨顿扩散公式。式。1 萨顿模式萨顿模式353 基于简单物理考虑，认为拉格朗日相关系数与湍流基于简单物理考虑，认为拉格朗日相关系数与湍流特征量，宏观黏滞度，时间间隔相关，并通过量纲特征量，宏观黏滞度，时间间隔相关，并通过量纲分析得到所有量的组合。分析得到所有量的组合。以以y向为例向为例：n为由风速梯度观测确定的实验常数，为由风速梯度观测确定的实验常数，m为风速廓线幂指数。为风速廓线幂指数。36萨顿参数萨顿参数最早，但有局限性最早，但有局限性37l萨顿基于拉氏相关函数和泰勒公式导出萨顿萨顿基于拉氏相关函数和泰勒公式导出萨顿广义扩散参数广义扩散参数;他的功绩在于将大气扩散参数与他的功绩在于将大气扩散参数与可测量的气象参量联系起来；可测量的气象参量联系起来；l经验系数经验系数N,n,m 需要通过风速梯度观测才能确需要通过风速梯度观测才能确定定;萨顿曾在平坦地形做了小尺度扩散观测实验萨顿曾在平坦地形做了小尺度扩散观测实验,直接利用这些数据计算高架连续点源的扩散浓度直接利用这些数据计算高架连续点源的扩散浓度将比实际观测值偏高。数理推导也不严密，应用将比实际观测值偏高。数理推导也不严密，应用也很不方便。故目前已不大使用。也很不方便。故目前已不大使用。38和和 由双向风标测量，反映大气湍流扩由双向风标测量，反映大气湍流扩散能力。散能力。2 直接测量湍流特征量的方法直接测量湍流特征量的方法H.E.Cramer(1957)提出提出脉动风方位脉动风方位角标准差角标准差脉动风高度脉动风高度角标准差角标准差p，g与稳定度、下风向距离及地表粗糙度相关与稳定度、下风向距离及地表粗糙度相关393 BNL模式模式M.E.Smith,1951l特征量：水平风向摆动角的范围特征量：水平风向摆动角的范围l高架源高架源108m高塔施放油雾扩散试验高塔施放油雾扩散试验l简便、合理、实用，美国机械工程师协会沿用简便、合理、实用，美国机械工程师协会沿用至今至今404 J.S.Hay&F.Pasquill(1959)出发点：出发点：统计理论，泰勒公式统计理论，泰勒公式方法：方法：利用相关函数和湍流能谱关系，由湍流观利用相关函数和湍流能谱关系，由湍流观测资料做谱分析，计算扩散参数。测资料做谱分析，计算扩散参数。总结：总结：模型合理可取，反映湍流场本质，而且准模型合理可取，反映湍流场本质，而且准确度较高，其探讨有一定理论意义，但应用尚确度较高，其探讨有一定理论意义，但应用尚不普遍，观测要求高，计算工作量大。不普遍，观测要求高，计算工作量大。41 由大量扩散试验含气象观测和示踪物浓度观测资料分析及理论分析得出扩散参数随下风方距离x的变化曲线 P-G法，或者P-G-T法英国气象学家帕斯奎尔(F.Pasquill，1961)基于大量扩散实验资料的分析,建立了一套扩散参数计算方案,后经美国核气象学家杰富德(F.Giffod,1961)和美国公共卫生局的气象学家特纳尔(D.B.,Turner,1967)的改进与完善,构成了现今广为引用的P-G-T 扩散曲线系统。二二 稳定度扩散分级与扩散曲线法稳定度扩散分级与扩散曲线法421.P-G曲线法曲线法方法要点方法要点大气分成大气分成A-F共六个稳定度等级共六个稳定度等级 云，日照，风速。云，日照，风速。xy曲线六条曲线六条 对应对应A、B.F稳定度级稳定度级 43l帕斯奎尔帕斯奎尔(1961)按表给出划分按表给出划分6个稳定度扩个稳定度扩散级别散级别:A级级极不稳定极不稳定,B级级中等不稳定中等不稳定,C级级弱不稳定弱不稳定,D级级中性中性,E级级弱稳定弱稳定,F级级 中等稳定。中等稳定。44lP-G曲线的应用曲线的应用l根据常规资料确定稳定度级别根据常规资料确定稳定度级别Pasquill1961表表45太阳高度角太阳高度角云量云量Turner(1961)引入太阳高度角判定日射强弱的引入太阳高度角判定日射强弱的定量方法，确定稳定度级别。定量方法，确定稳定度级别。日射等级日射等级风速风速稳定度稳定度46表表3.6 日射等级确定规那日射等级确定规那么么47表表3.7 Turner的稳定度分级方法的稳定度分级方法48lPG曲线的应用曲线的应用10分钟平均图分钟平均图l利用扩散曲线确定利用扩散曲线确定 和和水平扩散参数水平扩散参数垂直扩散参数垂直扩散参数49l用用P-G扩散曲线方法确定扩散参数的步骤是：扩散曲线方法确定扩散参数的步骤是：l根据太阳高度角和云量确定日射等级根据太阳高度角和云量确定日射等级(见表见表3.6)；l根据日射等级和地面风速确定稳定度等级见根据日射等级和地面风速确定稳定度等级见表表3.7)；l根据稳定度等级和距离根据稳定度等级和距离x可查出扩散参数见可查出扩散参数见图图3.10)。50lP-G曲线的应用曲线的应用l地面最大浓度估算地面最大浓度估算HzxsyxsmaxC由由和和由由曲线（图曲线（图3.10b）反查出反查出由由曲线（图曲线（图3.10a）读出）读出由式由式求出求出HzxsmaxcxyxsmaxC由由和和由由由由由式由式ys3.10按照表按照表3.3或或3.7定出所处稳定度扩散级别定出所处稳定度扩散级别51三三 扩散曲线法的完善扩散曲线法的完善1.国家标准中的修改应用国家标准中的修改应用GB/T13201-91修正太阳高度角的计算方法修正太阳高度角的计算方法适应我国大量地面观测无云高观测适应我国大量地面观测无云高观测的情况的情况5253平原地区和城市远郊区，平原地区和城市远郊区，D、E、F向不稳定向不稳定方向提半级方向提半级工业区和城市中心区，工业区和城市中心区，C提至提至B级，级，D、E、F向不稳定方向提一级向不稳定方向提一级丘陵山区的农村或城市，同工业区丘陵山区的农村或城市，同工业区中国国家标准规定的方法中国国家标准规定的方法542 不同稳定度分类方法不同稳定度分类方法1风向脉动标准差风向脉动标准差EPA，19905556l以风速做细致调整，观测数据在粗糙度以风速做细致调整，观测数据在粗糙度z0=15cm，10m高度处测量得到。采样时间为高度处测量得到。采样时间为15min。l如果风向发生转折，为了尽量减小风向转折的影响，如果风向发生转折，为了尽量减小风向转折的影响，应该将长时间段分成小段进行计算，例如将应该将长时间段分成小段进行计算，例如将60min的的时间划分为时间划分为15min一段，最终小时量值：一段，最终小时量值：572与温度递减率有关方法与温度递减率有关方法 以温度递减率，即以两层以温度递减率，即以两层10m、60m的垂直温度梯度来表征水平和垂直向的垂直温度梯度来表征水平和垂直向的湍流状况。国际原子能机构的湍流状况。国际原子能机构(1980)推推荐荐58 以温度递减率来表征大气稳定度的方法对于以温度递减率来表征大气稳定度的方法对于稳定大气的情形比较可靠，不稳定或者高架源稳定大气的情形比较可靠，不稳定或者高架源的情形适宜用水平风向脉动标准差方法。的情形适宜用水平风向脉动标准差方法。59 以温度递减率和风速相结合同时考虑支配以温度递减率和风速相结合同时考虑支配湍流活动的机械和热力因子湍流活动的机械和热力因子60 分别以温度梯度和分别以温度梯度和 表征湍流的垂表征湍流的垂直和水平运动直和水平运动 61EPA(1990)推荐在缺乏云量和云高资料时，推荐在缺乏云量和云高资料时，采用表替换原采用表替换原P-G-T方法方法70035035050623边界层湍流参量法边界层湍流参量法6364注意：注意：l除了使用公认的已有统一标准的方案，例如除了使用公认的已有统一标准的方案，例如P-P-G-TG-T方案，国家标准给出的修改方案，采用其方案，国家标准给出的修改方案，采用其他任何方案都应当验证其可行性，提供充分的他任何方案都应当验证其可行性，提供充分的实验依据和例证，必要时还应做专门的论证实验依据和例证，必要时还应做专门的论证653 不同源高和不同下垫面的应用不同源高和不同下垫面的应用lP-G扩散曲线实验依据：平坦理想条件，大量扩散曲线实验依据：平坦理想条件，大量低矮源扩散试验。低矮源扩散试验。l不同源高不同源高l不同下垫面不同下垫面66lBriggs，1974内插完善内插完善4 扩散曲线法的内插完善扩散曲线法的内插完善67四 风向脉动和扩散函数法l扩散曲线法试验根底存在较大的经验性，扩散曲线法试验根底存在较大的经验性，方案结果有许多不确定性方案结果有许多不确定性l仅以宏观气象状况为判据，在稳定度级仅以宏观气象状况为判据，在稳定度级别和湍流特性之间缺少清晰关系别和湍流特性之间缺少清晰关系l目前重要开展由风向脉动与扩散函数确目前重要开展由风向脉动与扩散函数确定扩散参数的方法定扩散参数的方法68(1)方法原理方法原理按照湍流扩散理论，在均匀定常条件下，粒子位按照湍流扩散理论，在均匀定常条件下，粒子位移的总体平均由泰勒公式表示。移的总体平均由泰勒公式表示。69由泰勒公式可得由泰勒公式可得 为拉格朗日时间尺度为拉格朗日时间尺度f为普适函数，扩散参数，函数形式随为普适函数，扩散参数，函数形式随源高源高和和稳定度稳定度变化变化70方法原理与湍流统计理论根底一致方法原理与湍流统计理论根底一致舍弃别离的稳定度级别，采用连续性稳定度，舍弃别离的稳定度级别，采用连续性稳定度，接近实际接近实际考虑源高影响，认为考虑源高影响，认为f是稳定度状况函数是稳定度状况函数使用方便，可用于多种情况使用方便，可用于多种情况(2)特点特点71Pasquill(1976)给出试验资料所得的给出试验资料所得的f(x)数据数据由泰勒公式积分可得由泰勒公式积分可得(3)扩散函数扩散函数f 确实定确实定72两者中间范围一致，近范围，理论值稍高，远距离两者中间范围一致，近范围，理论值稍高，远距离相反。相反。73l由试验资料分析求取扩散函数的方由试验资料分析求取扩散函数的方法法-自学自学P93744几种扩散函数表达式几种扩散函数表达式75五五 扩散参数的现状和开展扩散参数的现状和开展自学完成自学完成76重点重点l大气稀释因子、烟流有效源高、虚拟点大气稀释因子、烟流有效源高、虚拟点源、源、lP-G-T方法方法l国家标准国家标准77一一 烟流抬升烟流抬升(掌握掌握)有效源高有效源高4.4 烟流抬升高度烟流抬升高度78烟气抬升取决于烟气抬升取决于：动力因子：烟气排放的出口速度动力因子：烟气排放的出口速度热力因子：烟气温度比环境气温高热力因子：烟气温度比环境气温高791烟流抬升的物理模型及影响因子烟流抬升的物理模型及影响因子连续点源的排放模型连续点源的排放模型垂直烟流无风垂直烟流无风弯曲烟流有风弯曲烟流有风80 喷出阶段喷出阶段 浮升阶段浮升阶段 瓦解阶段瓦解阶段 变平阶段变平阶段抬升过程分为如下几个阶段抬升过程分为如下几个阶段(掌握掌握)81影响热烟流抬升的根本因子：影响热烟流抬升的根本因子：排放源及排放烟气的性质源排放排放源及排放烟气的性质源排放烟气的初始动量和浮力烟气的初始动量和浮力 环境大气的性质环境大气的性质下垫面性质下垫面性质822 烟流抬升轨迹与终极抬升烟流抬升轨迹与终极抬升烟气排出烟气排出准垂直形式准垂直形式弯曲并扩大弯曲并扩大瓦解变平瓦解变平 终极抬升高度终极抬升高度 烟流抬升轨迹烟流抬升轨迹:z(x)终极抬升距离：终极抬升距离：xT83 烟流是由运动中的流体构成的，于是，烟流是由运动中的流体构成的，于是，大多数烟流抬升模式基于流体力学的根本大多数烟流抬升模式基于流体力学的根本定律，即质量守恒，浮力守恒和动量守恒定律，即质量守恒，浮力守恒和动量守恒来考察烟流抬升问题，通过积分这些守恒来考察烟流抬升问题，通过积分这些守恒表达式表达式(在烟流的截面上积分在烟流的截面上积分)，以方便的，以方便的形式来描述烟流轨迹。形式来描述烟流轨迹。84二二 烟气抬升高度的计算烟气抬升高度的计算1 Holland公式公式1953：适用于中：适用于中性大气条件稳定时减小，不稳时增加性大气条件稳定时减小，不稳时增加1020 计算结果明显低估计算结果明显低估2-3倍倍852Briggs公式：适用不稳定及中性大气条件公式：适用不稳定及中性大气条件 QH:热释放率，单位时间，单位质量烟气升高热释放率，单位时间，单位质量烟气升高T度所释度所释放的热量，单位放的热量，单位cal/s。烟流抬升高度与烟囱高度有关，烟囱高情况下烟流抬升高度与烟囱高度有关，烟囱高情况下计算保守，相反情况计算不平安计算保守，相反情况计算不平安863国家标准国家标准(GB/T13201-91)中的公式中的公式 87 因为对烟流抬升和地面浓度的计算，因为对烟流抬升和地面浓度的计算，平均风速都具有重要作用，所以在高斯扩散平均风速都具有重要作用，所以在高斯扩散公式和烟流抬升高度计算公式中都有公式和烟流抬升高度计算公式中都有平均风平均风速速。一般认为这个量值应该是整个烟流厚度。一般认为这个量值应该是整个烟流厚度范围里各层风速的平均值。实际工作中常采范围里各层风速的平均值。实际工作中常采用源高处的平均风速，或者也有采用有效源用源高处的平均风速，或者也有采用有效源高处的平均风速。高处的平均风速。88常用地面风速资料推算要求高度常用地面风速资料推算要求高度处的风速处的风速89l烟流抬升过程经过的几个阶段烟流抬升过程经过的几个阶段l终极抬升距离终极抬升距离l烟流热释放率烟流热释放率l浮力参数浮力参数l烟流抬升高度实用计算公式烟流抬升高度实用计算公式l国家标准中的计算国家标准中的计算l利用风速幂指数率的计算利用风速幂指数率的计算904.5 其他非扩散过程其他非扩散过程91非扩散过程非扩散过程:1.烟流抬升作用，增加烟源有效高度烟流抬升作用，增加烟源有效高度2.污染物干沉积、湿沉积和由化学变化形成的迁污染物干沉积、湿沉积和由化学变化形成的迁移转化移转化3.固体下边界下垫面的作用固体下边界下垫面的作用92一一 大气去除过程的一般表述大气去除过程的一般表述1.根本机制：根本机制：(1)干沉积干沉积由地面的土壤、水、植物、建筑物由地面的土壤、水、植物、建筑物等通过污染物质的重力沉降、碰撞与捕获、吸收等通过污染物质的重力沉降、碰撞与捕获、吸收与吸附、光合作用或其它生物、化学、物理过程与吸附、光合作用或其它生物、化学、物理过程实现；实现；93(2)湿沉积湿沉积由云、雾和降水由云、雾和降水(雨，雪等形式雨，雪等形式)等通过污染等通过污染物质被吸收进入水滴或随水滴被去除。物质被吸收进入水滴或随水滴被去除。此外，化学转化亦是一种去除机制，但它是由一种物质此外，化学转化亦是一种去除机制，但它是由一种物质转换成另一种物质，因此，它只是对原生污染物的去除，转换成另一种物质，因此，它只是对原生污染物的去除，同时却又滋生新的次生污染物。同时却又滋生新的次生污染物。94湿沉积湿沉积：令气体污染物在令气体污染物在(x,y,z)位置，位置，t时刻的平时刻的平均浓度为均浓度为q(x,y,z;t)，而平均的气溶胶尺度分布函，而平均的气溶胶尺度分布函数为数为F(Dp；x,y,z,t)。函数函数:冲洗系数冲洗系数,与与离地高度离地高度和和时间时间有关有关污染气体去除率污染气体去除率气溶胶去除率气溶胶去除率95干沉积干沉积 定义为那些使污染物质在地面被消除的过程。定义为那些使污染物质在地面被消除的过程。物质垂直向下的通量可用一个称之为沉积速度物质垂直向下的通量可用一个称之为沉积速度的经验参数的经验参数 乘上某高度乘上某高度 的物质浓度。的物质浓度。96沉积速度沉积速度 与以下诸因素有关：与以下诸因素有关：(1)被去除的粒子种类或物质；被去除的粒子种类或物质；(2)表征近地面层状态的特性的气象参数；表征近地面层状态的特性的气象参数；(3)地面本身的性质。地面本身的性质。97 对地面的湿通量是自空中所有体积单元湿去除之和。98根据此定义湿沉积速度根据此定义湿沉积速度如果正被去除的物质在厚度为如果正被去除的物质在厚度为H的一层里垂的一层里垂直分布是均匀的直分布是均匀的99定义冲洗比定义冲洗比(washout ratio)湿沉积速度与冲洗比之间有关系湿沉积速度与冲洗比之间有关系？100二二 高斯烟流扩散公式的修正形式高斯烟流扩散公式的修正形式经修正的高架连续点源高斯烟流扩散公式经修正的高架连续点源高斯烟流扩散公式:大粒子重大粒子重力沉降力沉降地面反射系数地面反射系数半衰期半衰期随距离变化随距离变化101大粒子重力沉降的影响。大粒子重力沉降的影响。下垫面反射系数。下垫面反射系数。污染物半衰期时间。污染物半衰期时间。源项随距离源项随距离x的变化的变化102三三 干沉积干沉积-下垫面去除下垫面去除1.大粒子的重力沉降大粒子的重力沉降 大气中各种处于气溶胶状态的粒子污染物有大气中各种处于气溶胶状态的粒子污染物有一定的尺度分布，对于粒径大于一定的尺度分布，对于粒径大于10m的大粒子的大粒子在大气中运动时会产生重力沉降。在大气中运动时会产生重力沉降。103 对于粒径对于粒径2060 m范围的粒子，粘滞性不大范围的粒子，粘滞性不大的情况下，可采用斯托克斯公式计算它们的末速的情况下，可采用斯托克斯公式计算它们的末速度度 ，其大小与粒子的密度，其大小与粒子的密度 有关，即有有关，即有粒子半径粒子半径空气动力空气动力粘滞系数粘滞系数空气空气密度密度104对于粒径更大的粒子，须对斯托克斯公式修正对于粒径更大的粒子，须对斯托克斯公式修正对粒径对粒径200 m的粒子，其沉降速度的粒子，其沉降速度100cm/s，可，可按弹道计算。按弹道计算。对粒子沉降速度在对粒子沉降速度在100cm/s，粒径，粒径200 m，可采，可采用偏斜烟流模式处理，修改高斯扩散公式，考虑重用偏斜烟流模式处理，修改高斯扩散公式，考虑重力沉降的影响。力沉降的影响。105106当考虑以下斜烟流轴线为基准的扩散时，有当考虑以下斜烟流轴线为基准的扩散时，有107考虑地面反射作用考虑地面反射作用1082 气体和小粒子的干沉降气体和小粒子的干沉降 气体和很小的粒子(粒径小于10-20 m)上述沉降作用可忽略不考虑。但由于湍流扩散和布朗运动沉积到各种外表。吸收、碰撞、光合作用和其它生物学、化学和物理学过程会使物质沉积到地面。109 一般都假定下垫面对这类小粒子的去除作用不影响大气中污染物的浓度分布形式，而只影响大气中污染物的总量，或者说影响有效源强。最普遍采用的一种方法是所谓“源损耗模式或称有效源强法。源项随距离源项随距离x的变化的变化110 干沉积主要是优先去除近地面处的污染，而并非干沉积主要是优先去除近地面处的污染，而并非在整个烟流厚度层范围均匀起作用的。在整个烟流厚度层范围均匀起作用的。故此，源损耗模式尚有必要进一步改善。故此，源损耗模式尚有必要进一步改善。111 重力沉降、湍流运动、布朗运动、惯性作用和静电作用等，是形成干沉积的主要物理过程。化学反响、溶解、解吸等那么是形成干沉积的主要化学过程。对于植被那样的生物沉积外表，植被子生长的形态特征、生命过程以及静电性质等那么是影响干沉积的主要生物学特征。这几方面过程都会受气象条件、污染物性质和沉积外表的特性的影响。1123 干沉积速度的理论计算干沉积速度的理论计算干沉积过程可视为三步：干沉积过程可视为三步：由湍流扩散支配物质向贴地层输送；由湍流扩散支配物质向贴地层输送；物质通过紧贴地面的片流子层向吸收体扩散，物质通过紧贴地面的片流子层向吸收体扩散，称为地面输送；称为地面输送；由在地面的物质可溶性或吸收率确定通过片流由在地面的物质可溶性或吸收率确定通过片流子层扩散的物质实际上有多少被消除，这个最后子层扩散的物质实际上有多少被消除，这个最后的过程称为转移过程。的过程称为转移过程。113 根据干沉积速度的定义，它是垂直通量与实根据干沉积速度的定义，它是垂直通量与实测浓度之间的比例常数。整个沉积的三个步骤测浓度之间的比例常数。整个沉积的三个步骤可类同于电流或热量传送的情形。可类同于电流或热量传送的情形。气体气体粒子粒子114四四 湿沉积湿沉积大气中的雨、雪等降水形式和其它水汽凝结物，大气中的雨、雪等降水形式和其它水汽凝结物，如云、雾、霜等都能对空气污染物，包括气体和如云、雾、霜等都能对空气污染物，包括气体和粒子起到去除作用，称为降水去除或湿沉积。粒子起到去除作用，称为降水去除或湿沉积。通常，把由降水造成的污染物去除过程称为雨洗通常，把由降水造成的污染物去除过程称为雨洗(或雪洗或雪洗)，这种过程将空气污染带到地面。按照，这种过程将空气污染带到地面。按照降水去除过程发生所在高度分成云下清洗降水去除过程发生所在高度分成云下清洗(washout)和云中清洗和云中清洗(rainout)。115116 为模拟上述降水去除过程对空气污染物散布的影响，一般采用两种方法，即定义去除系数和清洗比的处理方法。1：去除系数的处理方法2：清洗比的处理方法117欧拉方法污染物平流扩散方程欧拉方法污染物平流扩散方程平流输送项平流输送项湍能扩散项湍能扩散项源源排排放放项项干干湿湿沉沉降降化化学学反反应应118五五 化学转换和空气污染物的滞留与迁移化学转换和空气污染物的滞留与迁移 排放进入大气层的烟流中，各种化学物质和大气成分混合输送，在适当的气象条件下，如辐射、温度、湿度和降水等，会发生化学反响以致变性或生成二次污染物，经过在大气层中一定时间的滞留而后迁移。119对于空气质量影响的模拟，最普通并具典型意对于空气质量影响的模拟，最普通并具典型意义的此类变化与迁移的例子有如：义的此类变化与迁移的例子有如：(1)由煤的燃烧和其它工业源排放由煤的燃烧和其它工业源排放SO2和和NOx初级初级空气污染物氧化生成硫酸盐、硝酸盐和臭氧等空气污染物氧化生成硫酸盐、硝酸盐和臭氧等二次污染物，对区域性在大气酸性和降水酸化二次污染物，对区域性在大气酸性和降水酸化以及臭氧成分的变化有重大作用以及臭氧成分的变化有重大作用120(2)城市汽车废排放和其它工业源排放的城市汽车废排放和其它工业源排放的HC化合物和化合物和NOx等初级空气污染物在太阳辐射作用下生成光化学烟等初级空气污染物在太阳辐射作用下生成光化学烟雾的过程，其中主要组成成分有臭氧和过氧乙醛硝酸盐雾的过程，其中主要组成成分有臭氧和过氧乙醛硝酸盐(PAN)等二次污染物，对城市区城空气质量和人体健康等二次污染物，对城市区城空气质量和人体健康有重要影响；有重要影响；(3)放射性污染物质的衰变与迁移，如燃料中含有一定放放射性污染物质的衰变与迁移，如燃料中含有一定放射性物质，核反响设施的裂变产物的泄漏或废弃物中的射性物质，核反响设施的裂变产物的泄漏或废弃物中的微量残留物等；微量残留物等；1214粒子和气体污染物在干湿去除过程中的化粒子和气体污染物在干湿去除过程中的化学反响和各种途径的迁移，例如，上节所述气学反响和各种途径的迁移，例如，上节所述气体向水滴输送以及在水滴中可能发生的一些化体向水滴输送以及在水滴中可能发生的一些化学反响，生成的二次污染物质被去除到地，影学反响，生成的二次污染物质被去除到地，影响地面空气污染物浓度的分布并造成危害。响地面空气污染物浓度的分布并造成危害。122123通常，采用最为简单的处理方式，即假设一个通常，采用最为简单的处理方式，即假设一个随时间常数随时间常数 变化的污染物浓度指数化学衰减变化的污染物浓度指数化学衰减率，即有率，即有124放射性衰变放射性衰变放射性衰变亦是能影响污染物放射性衰变亦是能影响污染物(其中含有放射性其中含有放射性物质物质)迁移的一种简单的去除过程。迁移的一种简单的去除过程。例如，燃料煤和其它工业矿物原料的燃烧都排放例如，燃料煤和其它工业矿物原料的燃烧都排放含有一定量的放射性物质，如镭含有一定量的放射性物质，如镭(Ra)，它们被引，它们被引入城市大气层等等。显然，它们的浓度通常都很入城市大气层等等。显然，它们的浓度通常都很低，而且都遵循日渐降低简单的衰变规律低，而且都遵循日渐降低简单的衰变规律125污染物滞留时间污染物滞留时间空气污染物滞留时间亦称居留时间空气污染物滞留时间亦称居留时间(residence time)是表是表征空气污染物在一定时期，一定空间范围里的寿命的一征空气污染物在一定时期，一定空间范围里的寿命的一个特征量。个特征量。它定义为这样一个时间，即假定没有污染物补充，没有它定义为这样一个时间，即假定没有污染物补充，没有污染物穿过区域边界，该区域大气里污染随时间减少至污染物穿过区域边界，该区域大气里污染随时间减少至原来的原来的1/e倍所需的时间。倍所需的时间。126污染物在大气中的物理污染物在大气中的物理化学衰减化学衰减综合考虑上述所有导致空气污染物可能被去除或综合考虑上述所有导致空气污染物可能被去除或迁移的物理和化学过程，滞留时间或居留期反映迁移的物理和化学过程，滞留时间或居留期反映了空气污染物衰变的快慢。实际的空气污染物浓了空气污染物衰变的快慢。实际的空气污染物浓度，在综合了各种去除机制后，仍假定它是随时度，在综合了各种去除机制后，仍假定它是随时间呈指数形式变化：间呈指数形式变化：127为了表征空气污染物的被去除的快慢，可以定义为了表征空气污染物的被去除的快慢，可以定义一个半衰期或称半生命期一个半衰期或称半生命期()的特征时间，意的特征时间，意即由于去除或迁移作用使空气污染物浓度衰减为即由于去除或迁移作用使空气污染物浓度衰减为原来浓度的一半所需的时间。原来浓度的一半所需的时间。由此可得：由此可得：1284.6 各种条件下的扩散计算各种条件下的扩散计算1294.6.1 4.6.1 封闭型扩散封闭型扩散130扩散被局限在地面和扩散被局限在地面和上部逆温层底之间，上部逆温层底之间，具有两个反射面的特具有两个反射面的特征。应用像源法可以征。应用像源法可以导出封闭型扩散公式，导出封闭型扩散公式，实源在两个反射面的实源在两个反射面的另一侧等距离处都能另一侧等距离处都能成像成像131132l实际工作中按三个阶段作简化处理实际工作中按三个阶段作简化处理l近距离近距离(xxu)：污染物经过两个界面的屡次反射在垂直方向为均匀分布，但在横风向仍属于Gauss分布。此时污染物散布的公式l离源距离x在xDxxu之间，地面浓度采用内差的方法 1344.6.2 熏烟型扩散 l 晴夜，地面辐射冷却形成自地面向上的逆温层，其上部仍为递减状态。l日出后，地面受太阳照射增温，夜间形成的逆温层自下而上逐渐消失，底部变成中性或不稳定层结。l原来积聚于稳定层结中的污染物，一方面由于上部仍有未消散完全的逆温层，污染物的扩散受抑制；135 另一方面由于下部逆温层消失一定高度后变为在中性或不稳定层结垂直混合到地，从而形成熏烟型熏烟型污染。l发生污染时，污染的范围较小，时间较短，但可以造成很高的地面浓度。136137l熏烟扩散的三种形式：l1逆温破坏，出现热对流，形成的混合层l2城市热岛效应生成城市混合层l3水陆交接地区岸上热力内边界层138l熏烟型扩散与封闭型扩散的不同1)熏烟型扩散的烟流原先处于稳定层结中，随后污染物由积聚到熏蒸到地；而封闭型扩散烟流原先处于不稳定层结内，污染物并无积聚，对地面影响较轻。2)熏烟型扩散过程持续的时间很短暂，而封闭型扩散过程往往持续相当长的时间。3)熏烟型扩散过程造成的地面污染物高浓度落地点范围窄小，但离源距离变化很大，视逆温强度、平均风速等而定；封闭型扩散过程造成的范围较大，影响相对稳定。1394.6.3 微风扩散 l微风逆温条件下的扩散有两个重要特性：1).垂直和水平扩散有明显差异，不能划分在同一稳定度类别，应分别处理。2).微风时风向不稳定，常伴有缓慢的大幅度的气流偏转，可能造成较大的水平扩散。1404.6.4 准静风烟团模式l该情况下扩散处理需考虑lx方向的扩散l扩散参数随时间变化而且水平扩散增大l风场随时间变化，还需考虑大湍涡的作用1414.6.5 长期平均浓度的计算l运用高斯模式计算的都是准定常条件下，假定风向风速和稳定度状况不随时间变化的。因此只适合于短时间，如1小时平均浓度估算。l计算长时间的浓度，如月、季、年平均浓度，不能再把风向风速和稳定度视为不变的量了。l假设在计算时段内污染源为定常，只考虑气象条件的变化既可。142l大气去除过程、干沉积、湿沉积大气去除过程、干沉积、湿沉积l高斯烟流扩散公式的修正形式高斯烟流扩散公式的修正形式l斯托克斯公式斯托克斯公式l偏斜烟流模式偏斜烟流模式l有效源强法有效源强法l干沉积速度的理论计算干沉积速度的理论计算l去除系数、清洗比、污染物滞留时间去除系数、清洗比、污染物滞留时间l封闭型扩散和熏烟型扩散封闭型扩散和熏烟型扩散