式中，Q1，ΔH1——分别为修正后的源强和抬升高度；
h——混合层高度。
2.风向切变和热浮力问题
高浮力烟羽在z方向受风向切变及热浮力的影响可按下式对σy和σz进行修正：
σ2yc=σ2y+（ΔH/3.5）2+（0.03Δα2x2）
σ2zc=σ2z+（ΔH/3.5）2
式中，σyc、σzc——分别为修正后y、z方向的扩散参数；
（ΔH/3.5）2——热浮力修正项；
（0.03Δα2x2）——风向切变修正项。
Δα按下表确定：
表风向切变修正参数Δα/ΔZ（度/m）
稳定度
A
B
C
D
E
F
Δα/ΔZ
0.005
0.010
0.015
0.020
0.025
0.035
1.2.4.3干沉积修正
1.干沉积的含义
在重力、湍流扩散、分子扩散、静电引力以及其他生物学、化学和物理学等因素的作用下，大气中的颗粒物或某些气体随时会被地表（土壤、植物、水体）滞留或吸收，使这些物质连续不断地从大气向地表作质量转移，从而减少其在空气中的浓度。
通常，把这一与降水作用无关的质量转移过程，称之为干沉积。
在对中距离大气污染物输送的研究中发现，差不多一半以上的质量转移是干沉积过程引起的。
为些，对一些物理、化学特性和空气相差较大的大气污染物（例如颗粒物），或者对于输送距离比较远的高架源的浓度预测，都应当考虑其干沉积的影响。
下面将介绍两种常用的干沉积扩散模式。
2.源亏损模式
源亏损模式主要用于粒径小于10μm易产生沉积的颗粒物或气体。
假定因各种机理造成的大气污染物向地面的沉积通量W[mg/（s·m2）]由下式表示：
W=VdC..........................................（4.9-79）
式中Vd——沉积速度，m/s；
C——大气污染物地面浓度，常取自距地面1m高度处。
大气污染物自烟囱出口排出后，其初始源强Q（0）因沉积作用将随下风距离逐渐减弱（亏损）。
根据式（4.9-79）可导出亏损后的源强Q（x）；因为将W从－∞到∞对y积分恰是—dQ（x）/dx，再对x积分则可得
由式（4.9-80）可见，问题的关键是给出沉积速度Vd。
长期以来，许多研究者为求得Vd，在实验或模式化方面都进行了大量的工作。
对于SO2，Hanna曾推荐如下Vd（cm/s）的野外实验值：
高约1m的庄稼地、城市、水面为0.7；高0.1m的草地为0.5；酸性土壤为0.4（干）或0.6（湿）；石灰质土壤为0.8；干积雪为0.1。
在模式化方面，通常是把沉积速度和电流类比，假定沉积速度与各种阻力的总和成反比。
上式可用数值积分计算。
采用1小时取样时间的国标扩散参数计算结果表明，D稳定度时亏损最明显，不稳定或更稳定时亏损依次减少。
设He=100m,Q（0）为100，Vd为0.007m/s，U为3m/s，则在10000m处Q为下表所示：
稳定度
A
B
B-C
C
C-D
D
D-E
E
F
Q（10km）
99.08
96.27
95.17
93.93
92.21
91.90
92.81
94.76
98.86
一般来说，3000m以内的扩散源亏损小于1%，可以忽略不计。
3.部分反射（斜烟羽）
部分反射模式主要用于粒径大于10μm的气载颗粒物。
因为这种污染物有明显的沉降作用，其烟羽重心是是逐渐降低的，此外，地面只能反射一部分的污染物。
对某一粒径的颗粒污染物，对He和Q进行以下调整后，再按相应气态物模式计算：
He（x）=He（0）-Vgx/u。
Q（x）=Q（0）（1+α）/2
式中α——颗粒物的地面反射系数（参阅下表）；
Vg——颗粒物沉降速度，用STOCKS公式计算：
；
d、ρ——分别为颗粒物的直径和密度；
g——重力加速度；
μ——空气粘性系数。
表地面反射系数α
粒径范围/μm
15~30
31~47
60
76~100
平均粒径/μm
22
38
50
85
反射系数α
0.8
0.5
0.3
0
1.2.4.4湿沉积修正
湿沉积系指大气污染物因降水而减少其在空气中浓度的过程。
严格地讲，湿沉积分云中和云下两种清除机制。
在工程应用中，常把这两种机制合起来考虑。
假设大气污染物的初始源强Q（0）因降水随下风距离x成指数衰减，则
Q（x）=Q（0）exp（—Λx/u）
式中x——接受点的下风距离；
u——烟囱出口处的风速；
Λ——清除系数，S－1。
令J为降水强度（mm/h），根据实验定测量结果，对于SO2
Λ=1.7×10－4J0.6
1.2.4.5化学迁移的修正
化学迁移的修正，在工程应用上与湿沉积类似，可令修正后的源强Q（x）等于初始源强Q（0）乘以修正因子fc。
式中Tc——大气污染物的时间常数；其他符号同前。
如果定义fc=1/2所对应的时间为该大气污染物的半衰期td，将其代入式（44.9-85），可得：
对于城市尺度，SO2的td≈4小时；长距离输送SO2的td典型值为4天。
1.3面源（体源）的点源修正算法
把面源的排放当作一个位于其几何中心的点源的排放，对扩散参数适当修正后，采用点源模式直接计算，用以近似代表该面源的扩散。
但这种方法仅适用于有风条件（U10≥1.5m/s）。
在小风静风的条件下，面源扩散如何计算，尚未有定论，这里暂按小风静风的点源方法计算。
如果测点位于面源（体源）之外，这种近似算法较为准确，但如果测点在面源内部，则计算结果就很不理想。
比如，在有风时，位于面源几何中心上风向各点（仍处于面源内部），按这种近似算法浓度应是0，实际不然。
此外，这种算法不能反映面源的形状。
这些都是点源修正算法的不足之外。
有风条件下，对于面源或体源内部的预测点，可采用下式计算：
式中
，x为测点离面源上风边界的距离。
这种方法使得在有风时，面源内任意一点均有浓度。
①面源。
采用直接修正法。
如果面源的面积较小（S≤1km2），面源外的Cs可按点源扩散模式计算，只是应附加一个初始扰动。
这一初始扰动使烟羽在x=0就有一个和面源横向宽度相等的横向尺度，以及和面源高度相等的垂直向尺度。
注意到烟羽的半宽度等于2.15σy或2.15σz，则修正后的σy和σz分别为：
式中x——自接受点至面源中心点的距离；
ay——面源在y方向的长度；
H——面源的平均排放高度。
②体源
当无组织排放源为体源时，地面浓度可按点源直接修正法计算；令αy、αz分别表示体源在y和z方向的边长，则修正后的σy和σz分别为：
x——自接受点至面源中心点的距离。
如果面源或体源为非正常排放源，可按上述修正后，采用相应的非正常排放点源模式。
但当风速U10<1.5m/s时，对于体源，可用在实际的时刻t中加一个初始时间t0的方法进行修正：
式中，ax、ay、az分别为体源在平均风向（x）、横向（y）、垂直方向（z）的边长。
1.4面源（体源）的数值积分算法
将面源（体源）分成一系列微小面源（dL×dW）,分别以点源代替这些微源，再用相应点源模式计算每一个微源对预测点的影响，再叠加所有微源在测点的浓度。
这种算法思维直接明了，算法准确，但无法用手工完成，以前计算机速度太慢也不现实，而现在的PC机的计算速度已经允许采用二次积分方法来实现这一算法，同时控制较高的精度。
设面源Ω是一封闭的区间，坐标系是以风向为正X轴，几何中心为原点，其几何中心坐标为（0,0）。
预测点在下风向X，横风向为Y，测点坐标（X,Y）。
对Ω内的任意一个点（x,y），设其代表一个面积为dx×dy的微源，则该微源对预测点的浓度贡献为：
C（X,Y）（x,y）=Qdxdyf（X-x,Y-y）
式中，Q是面源单位面积的源强[mg/（s.m2）]，Qdxdy为微元的源强（mg/s）。
f（）为点源计算公式，与风速、排放是否连续、稳定度等因素有关。
因此整个面源对（X,Y）的浓度等于Ω内所有点对（X,Y）浓度的叠加值，用积分式表示为：
f（X-x,Y-y）代表了面源中点x,y对测点（X,Y）的点源计算公式。
这个公式可以是有风模式，也可以是小风静风模式；可以是正常持续排放模式，也可是非正常排放的间断源模式。
还可以对源有效高和源强进行适当修正（如干沉、湿沉、化学迁移等等），总之，可以完全运用点源计算中的公式，详细见点源章节。
因Q实际上是与位置有关的（若要考虑源衰减），因此将Q归入f（）函数中，形成一个新的函数F（），则C（X,Y）的浓度为：
本软件中对上式采用变步长辛普生二重积分法进行积分，计算速度约为点源修正方法的百分之一。
如果采用PⅡ300以上处理器，这一速度还是可以接受的。
但若计算长期平均浓度时，一般PⅡ计算机显得太慢。
为简化求积过程和方便定义面源形状，本软件中所有面源（体源）均要求是矩形，但角度可作任意旋转。
优点：
可以准确描述面源与风向的关系，即使面源的边与风向不垂直时也能较好计算，面源内部的点也能算；可以直接使用点源公式，小风静风时也能计算。
从理论上来说，计算结果也比一般方法准确。
缺点：
计算量很大；要求准确定义面源位置和形状。
1.5多源叠加模式
如果需要评价的点源多于一个，计算浓度时，应将各个源对接受点浓度的贡献进行叠加。
在评价区内选一原点，以平均风的上风方为正x轴，各个源（坐标为xr,yr,0）对评价区内任一地面点（x,y）的浓度总贡献Cn可按下式计算：
式中Cr是第r个点源对（x,y,0）点的浓度贡献，其计算公式可根据不同条件选用本章给出的有关点源模式，但是注意坐标变换，（x,y,0）代以（x-xr，y-yr）。
1.6长期平均浓度计算方法
1.6.1孤立源长期平均浓度公式
当平均时间超过1小时之后，由于风向的摆动，任一风方位内的污染物浓度在横向都将趋于均匀分布。
为此，可将连续点源模式对y积分，并除以接受点所在位置的风方位宽度（或弧线长度）。
对于孤立排放源，以烟囱地面位置为原点，在某一稳定度（序号为j）和平均风速（序号为k）时，任意风向方位i的下风方x处的长期平均浓度（季、期或年均值）Cijk（x）（mg·m－3）为：
式中n——风向方位数，一般取16（2πx/16即为x处22.5°圆心角对应的弧长）；其他符号同前。
在可能出现的稳定度和平均风速条件下，任意风向位i的下风方x处的长期平均浓度Ci（x）（mg·m－3）为：
式中fijk——有风时风向方位、稳定度、风速联合频率；
Cijk——对应于该联合频率在下风方x处有风时的浓度值，由式（4.9-35）给出；
fLijk——静风或小风时，不同风方位和稳定度的出现频率（下标k只含有静风和小风两个风速段）；
CLijk——对应于fLijk的静风或小风时的地面浓度。
因为静风或小风时的风脉动角本来就比较大，CLijk可直接按小风静风模式计算。
式（4.9-36）中的j和k的加总数取决于所划分的稳定度和风速段数目，j的总数不宜少于3（稳定、中性、不稳定）；有风时k的总数一般也宜少于3。
在估算每个风速段的平均风速时，由于平均风速出现在公式的分母中，因而平均风速应等于单次风速倒数的平均值倒数。
其表达式为：
式中ui——第i个风速值；
N——总个数。
1.6.2多源长期平均浓度公式
如果评价区的烟囱多于一个，则任一接受点（x，y）的长期平均浓度为：
式中Crijk、CLrijk——分别是在接受点上风方对应于fijk和fLijk联合频率的第r个源对接收点的浓度贡献。
Crijk、CLrijk的公式形式分别和Cijk、CLijk相同（参阅4.9.6.1），但应注意坐标变换，将坐标转换到以接受点为原点，i风方位为正x轴的新坐标系后，再应用Cijk或CLijk公式。
1点源扩散模式
1.1持续排放源
1.1.1有风模式（U10≥1.5m/s）
1.1.2小风静风模式（U10<1.5m/s）
1.2非正常排放源
1.2.1有风情况（U10≥1.5m/s）
1.2.2小风静风（U10<1.5m/s）
1.3单源扩散的地面轴线最大浓度
1.4对源强和有效源高的修正
1.4.1地形对烟羽的影响
1.4.2热浮力对烟羽的修正
1.4.3干沉积修正
1.4.4湿沉积修正
1.4.5化学迁移的修正
2面源（体源）的点源修正算法
3面源（体源）的数值积分算法
4多源叠加模式
5长期平均浓度计算方法
5.1孤立源长期平均浓度公式
5.2多源长期平均浓度公式
2点源扩散模式
2.1持续排放源
2.1.1有风模式（U10≥1.5m/s）
1.自由空间中的连续点源
实际上绝大多数污染源都是连续的，对于连续排放源，可理解为在时间上依次连续释放无穷多个烟团。
因此，连续排放源的扩散模式可以通过将瞬时单烟团模式对to从--∞到t积分后求得。
以烟团初始空间坐标为原点，下风方为x轴，烟羽轴线与x轴一直保持重合，
都是x的函数，将对to的积分变换为对（x--uT）/σx的积分，可得最基本的烟羽扩散模式：
适用条件为：
自由空间；风速要比较大（u10≥1.5m/s）；当大气不稳定状态时，可能带来一定的误差。
2.地面反射
用像源法，假想地平线为一镜面，在其下方有一与真实源完
升级会员 