第一部分 大气扩散烟团轨迹模型 1 大气扩散烟团轨迹模型简介 该模型由国家环境保护总局环境规划院开发。
烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。
1.1 三维风场的计算 首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放烟团的时间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下 式中 Vt1、Vt2分别为第1和第2个观测时刻的风场值；
烟团释放时间步长；
n为t1、t2间隔内的时间步长数目；
Vi表示t1、t2间隔内第i个时间步长上的风场值。
1.2 烟团轨迹的计算 位于源点的某污染源，在t0时刻释放出第1个烟团，此烟团按t0时刻源点处的风向风速运行，经一个时间步长后在t1时刻到达P11，经过的距离为D11，从t1开始，第一个烟团按P11处t1时刻的风向风速走一个时间步长，在t2时刻到达P12，其间经过距离D12，与此同时，在t1时刻从源点释放出第2个烟团，按源点处t1时刻的风向风速运行，在t2时刻到达P22，其经过的距离为D22，以此类推，从t0时刻经过j个，到tj时刻共释放出了j个烟团，这时，这j个烟团的中心分别位于Pij，i1，2，j，设源的坐标为（Xs，Ys，Zst），Zst为t时刻烟团的有效抬升高度，Pij的坐标为（Xij，Yij，Zij），u、v分别为风速在X、Y方向的分量，则有如下计算公式 t1时刻 t2时刻 以此类推，到tj时刻，共释放出j个烟团，这些烟团最后的中心位置分别在Pij，Xij，Yij，Zij，i1，2， j，对于第i个烟团有 为i个烟团从源点释放后到tj时刻所经过的距离。
1.3 浓度公式 由前一个小节的计算，已找到由S点Xs，Ys的污染源释放出来的所有烟团在第j个时刻所处的位置，这样S处的污染源在第j个时刻在地面某接受点RX、Y、0处造成的浓度就是所有i个烟团的浓度贡献之和。考虑中心位于Pij的烟团对R点的浓度贡献，则有 式中Qs源强，mg/s；
X方向、Y方向、Z方向的大气扩散参数，m；
Cx、Cy、CzX、Y、Z方向扩散项，Cz在后面给出算式；
Cb为污染物转化项，b为转化率，1/s；
Cd为污染物沉降项，Vd为沉降速率，m/s。
由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，垂直扩散项分以下几种情况讨论 当混合层高为零时（即无混合层时）有 计算地面浓度时，Z0，则有 当混合层高度Zi不为零时，垂直扩散项分以下几种情况计算。
设排放源几何高度为hs，混合层高度为Zi，令，设烟气抬升高为（烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气穿透率，按不同的P值，分别计算Cz。
当P0，即时，认为污染物全在混合层内，按封闭性扩散式计算，即污染物在混合层与地面间多次反射。
式中N－为反射次数，一般取为N4即可。
当P＞1时，即 时，认为污染物完全穿透混合层，并在混合层以上的稳定层中扩散，由混合层的阻挡而不能到达地面，这时令Cz0。
当0 1000 BC 0.919325 0.875086 0.229500 0.314238 01000 1000 C 0.924279 0.885157 0.177154 0.232123 01000 1000 CD 0.926849 0.886940 0.143940 0.189396 01000 1000 D 0.929481 0.888723 0.110726 0.146669 01000 1000 DE 0.925118 0.892794 0.0985631 0.124308 01000 1000 E 0.920818 0.896864 0.086001 0.124308 01000 1000 F 0.929481 0.888723 0.0553634 0.073348 01000 1000 表 1.4-2 垂直扩散参数幂函数表达式数据 扩散参数 稳定度等级 PS α2 下风距离，m A 1.12154 1.5260 2.10881 0.0799904 0.00854771 0.000211545 0300 300500 500 B 0.941015 1.09356 0.127190 0.0570251 0500 500 BC 0.941015 1.00770 0.114682 0.0757182 0500 500 C 0.917595 0.106803 0 CD 0.838628 0.756410 0.815575 0.126152 0.235667 0.136659 02000 200010000 10000 D 0.826212 0.632023 0.555360 0.104634 0.400167 0.810763 11000 100010000 10000 DE o.776864 0.572347 0.499149 0.104634 0.400167 1.03810 02000 200010000 10000 E 0.788370 0.565188 0.414743 0.0927529 0.433384 1.73241 01000 100010000 10000 F 0.78440 0.525969 0.322659 0.0620765 0.370015 2.40691 01000 100010000 10000 3 丘陵山区的农村或城市，其扩散参数选取方法同工业区。
1.4.2 小风和静风（U10＜1.5m/s 时，0.5h取样时间的扩散参数按表1.4-3选取 表1.4-3 小风和静风扩散参数的系数、 稳定度（PS） U10＜0.5m/s 1.5m/s＞U10≥0.5m/s U10＜0.5m/s 1.5m/s＞U10≥0.5/s A 0.93 0.76 0.15 1.57 B 0.76 0.56 0.47 0.47 C 0.55 0.35 0.21 0.21 D 0.47 0.27 0.12 0.12 E 0.44 0.24 0.07 0.07 F 0.44 0.24 0.05 0.05 1.5 烟气抬升公式 1.5.1 有风时，中性和不稳定条件的烟气抬升高度△H（m） （1）当烟气热释放率Qh大于或等于是2100KJ/s，且烟气温度与环境温度的差值△T大于或等于35K时，△H采用下式计算 式中 no烟气热状况及地表系数，见表1.5-1；
n1烟气热释放率指数，见表1.5-1；
n2排气筒高度指数，见表1.5-1；
Qh烟气热释放率，KJ/s；
H排气筒距地面几何高度，m ，超过去240m时，取H240m；
Pa大气压力，hPa；
Qv实际排烟率，m3/s；
△T烟气出口温度与环境温度差，K；
Ts烟气出口温度，K；
Ta环境大气温度，K；
U排气筒出口处平均风速，m/s 。
表1.5-1 no、n1、n2的选取 Qh，KJ/s 地表状况（平原） no n1 n2 Qh，KJ/s 农村或城市远郊区 1.427 1/3 2/3 城市及近郊区 1.303 1/3 2/3 2100≤Qh＜21000 且△T≥35K 农村或城市远郊区 0.332 3/5 2/5 城市及近郊区 0.292 3/5 2/5 （2）当1700 kJ/s＜Qh＜2100KJ/s时， 式中Vs排气筒出口处烟气排出速度，m/s；
D排气筒出口直径，m；
△H2按（1）方法计算，no、n1、n2按表1.5-1中Qh值较小的一类选取；
Qh ，U与（1）中的定义相同。
（3）当Qh≤1700kJ/s或者△T＜35K时， 1.5.2 有风时,稳定条件按下式计算烟气抬升高度△Hm。
1.5.3 静风和小风时，按下式计算烟气抬升高度△Hm.。
但取值不宜小于0.01K/m。
2 模型运行所需数据 数据文件1 共四行 第1行 X方向网格点的数目（MX），Y方向网格点的数目（MY），Z方向风的观测数据层数目（MZ），最大有效烟团数（NT，默认110），污染源数目（MSC），气象观测小时数目（NTimes），稳定度数目（NeleTa，默认24）；
第2行 烟团的时间步长（分，默认30.0） 第3行 X方向步长（m，默认1000.0），Y方向步长（m，默认1000.0），大气压力hPa，默认1013.25，规划区类型（1农村，2城市）；
第4行 化学转化率1/s，沉降速率m/s；
数据文件2 网格点上的背景浓度值，单位mg/m3；
C1I,J,I1,MX,J1,MY 数据文件3 网格点上的高程（地形值）， 单位m；
数据文件4共六部分。
第1部分污染源数据，共12列，MSC行；
序号 1 2 3 4