大气环境容量计算模型.doc

附件二： 大气环境容量测算模型简介 说明：本部分内容是“重点城市大气环境容量核定工作方案”中提到的各推荐模型的简介，主要目的是为了使各城市了解各模型的功能和基本原理，同时，了解如选用该模型，都需要准备哪些输入数据，以便各城市根据本市的实际情况，提前准备。 第一部分 大气扩散烟团轨迹模型 1 大气扩散烟团轨迹模型简介 该模型由国家环境保护总局环境规划院开发。 烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。 1.1 三维风场的计算 首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放烟团的时间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下： 式中： V(t1)、V(t2)—分别为第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，i=1，2，…j，设源的坐标为（Xs，Ys，Zs(t)），Zs(t)为t时刻烟团的有效抬升高度，Pij的坐标为（Xij，Yij，Zij），u、v分别为风速在X、Y方向的分量，则有如下计算公式： 　t1时刻： t2时刻： 以此类推，到tj时刻，共释放出j个烟团，这些烟团最后的中心位置分别在Pij，Xij，Yij，Zij，i=1，2，… j，对于第i个烟团有： 为i个烟团从源点释放后到tj时刻所经过的距离。 1.3 浓度公式 由前一个小节的计算，已找到由S点(Xs，Ys)的污染源释放出来的所有烟团在第j个时刻所处的位置，这样S处的污染源在第j个时刻在地面某接受点R(X、Y、0)处造成的浓度就是所有i个烟团的浓度贡献之和。考虑中心位于Pij的烟团对R点的浓度贡献，则有： 式中：Qs—源强，mg/s； ：—X方向、Y方向、Z方向的大气扩散参数，m； Cx、Cy、Cz：—X、Y、Z方向扩散项，Cz在后面给出算式； Cb为污染物转化项，b为转化率，1/s； Cd为污染物沉降项，Vd为沉降速率，m/s。 由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，垂直扩散项分以下几种情况讨论： 当混合层高为零时（即无混合层时）有： 计算地面浓度时，Z=0，则有： 当混合层高度Zi不为零时，垂直扩散项分以下几种情况计算。 设排放源几何高度为hs，混合层高度为Zi，令，设烟气抬升高为（烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气穿透率：，按不同的P值，分别计算Cz。 当P=0，即时，认为污染物全在混合层内，按封闭性扩散式计算，即污染物在混合层与地面间多次反射。 式中：N－为反射次数，一般取为N=4即可。 当P＞1时，即　时，认为污染物完全穿透混合层，并在混合层以上的稳定层中扩散，由混合层的阻挡而不能到达地面，这时令Cz=0。 当0P1，即时认为是部分穿透情形，这时有部分污染物抬升到混合层以上，而（1-P）部分被封闭在混合层以内，Cz按下式计算： 许多文献认为穿透到混合层以上的污染物被阻挡后不能向地面扩散，当地区大气层结处于中性偏稳定结构时，混合层对污染物的阻挡作用并不是很强，这时可设计成让这部分烟团在高度上向下扩散，则有： 而（1-P）部分的烟团在Zi处按封闭扩散： 1.4 大气扩散参数 1.4.1 有风时扩散参数σy、σz的确定（0.5h取样时间） (1) 平原地区农村及城市远郊区的扩散参数选取方法如下： A、B、C级稳定度直接由表1.4-1和表1.4-2查算，D、E、F级稳定度则需向不稳定方向提半级后由表1.4-1和表1.4-2查算。 (2) 工业区或城区中的点源，其扩散参数选取方法如下： A、B级不提级，C级提到B级，D、E、F级向不稳定方向提一级，再按表1.4-1和表1.4-2查算。 表 1.4-1 横向扩散参数幂函数表达式数据 扩散参数 稳定度等级（P·S） α1 下风距离，m A 0.901074 0.850934 0.425809 0.602052