附件二：大气环境容量测算模型简介说明：本部分内容是“重点城市大气环境容量核定工作方案”中提到的各推荐模型的简介，主要目的是为了使各城市了解各模型的功能和基本原理，同时，了解如选用该模型，都需要准备哪些输入数据，以便各城市根据本市的实际情况，提前准备。第一部分大气扩散烟团轨迹模型1大气扩散烟团轨迹模型简介该模型由国家环境保护总局环境规划院开发。烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。1.1三维风场的计算首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放烟团的时间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下：Vi=V(t1)+[V(t2)−V(t1)]⋅inn=(t2−t1)/Δt式中：V(t1)、V(t2)—分别为第1和第2个观测时刻的风场值；Δt—烟团释放时间步长；n—为t1、t2间隔内的时间步长数目；Vi—表示t1、t2间隔内第i个时间步长上的风场值。1.2烟团轨迹的计算位于源点的某污染源，在t0时刻释放出第1个烟团，此烟团按t0时刻源点处的风向风速运行，经一个时间步长Δt后在t1时刻到达P11，经过的距离为D11，从t1开始，第一个烟团按P11处t1时刻的风向风速走一个时间步长，在t2时刻到达P12，其间经过距离D12，与此同时，在t1时刻从源点释放出第2个烟团，按源点处t1时刻的风向风速运行，在t2时刻到达P22，其经过的距离为D22，以此类推，从t0时刻经过j个Δt，到tj时刻共释放出了j个烟团，这时，这j个烟团的中心分别位于Pij，i=1，2，…j，设源的坐标为（Xs，Ys，Zs(t)），Zs(t)为t时刻烟团的有效抬升高度，Pij的坐标为（Xij，Yij，Zij），u、v分别为风速在X、Y方向的分量，则有如下计算公式：t1时刻：X11=Xs+U[t0,Xs,Ys,Zs(t0)]⋅ΔtY11=Ys+V[t0,Xs,Ys,Zs(t0)]⋅ΔtZ11=Zs+W[t0,Xs,Ys,Zs(t0)]⋅ΔtD1=D11=√(X11−Xs)2+(Y11−Ys)2t2时刻：X12=X11+U[t1,X11,Y11,Z11]⋅ΔtY12=Y11+V[t1,X11,Y11,Z11]⋅ΔtZ12=Z11+W[t1,X11,Y11,Z11]⋅ΔtD12=D11+D12=D11+√(X12−X11)2+(Y12−Y11)2X22=Xs+U[t1,Xs,Ys,Zs(t1)]⋅ΔtY22=Ys+V[t1,Xs,Ys,Zs(t1)]⋅ΔtZ22=Zs+W[t1,Xs,Ys,Zs(t1)]⋅ΔtD22=D22=√(X22−Xs)2+(Y22−Ys)2以此类推，到tj时刻，共释放出j个烟团，这些烟团最后的中心位置分别在Pij，Xij，Yij，Zij，i=1，2，…j，对于第i个烟团有：Xij=Xi(j−1)+U[tj−1,Xi(j−1),Yi(j−1),Zi(j−1)]⋅ΔtYij=Yi(j−1)+V[tj−1,Xi(j−1),Yi(j−1),Zi(j−1)]⋅ΔtZij=Zi(j−1)+W[tj−1,Xi(j−1),Yi(j−1),Zi(j−1)]⋅ΔtDij=∑k=1jDik=Dij−1+√(Xij−Xi(j−1))2+(Yij−Yi(j−1))2Dij为i个烟团从源点释放后到tj时刻所经过的距离。1.3浓度公式由前一个小节的计算，已找到由S点(Xs，Ys)的污染源释放出来的所有烟团在第j个时刻所处的位置，这样S处的污染源在第j个时刻在地面某接受点R(X、Y、0)处造成的浓度就是所有i个烟团的浓度贡献之和。考虑中心位于Pij的烟团对R点的浓度贡献，则有：Ci=Qs(2π)2/3σxσyσzCX⋅CY⋅CZ⋅Cb⋅CdCX=EXP[−(X−Xij)22σx2]⋅CY=EXP[−(Y−Yij)22σy2]Cb=EXP(−b⋅jΔt)Cd=EXP[−(Vd⋅jΔt)22σz2]式中：Qs—源强，mg/s；σx、σy、σz：—X方向、Y方向、Z方向的大气扩散参数，m；Cx、Cy、Cz：—X、Y、Z方向扩散项，Cz在后面给出算式；Cb为污染物转化项，b为转化率，1/s；Cd为污染物沉降项，Vd为沉降速率，m/s。由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，垂直扩散项分以下几种情况讨论：当混合层高为零时（即无混合层时）有：Cz=EXP[−(Z+Zij)22σz2]+EXP[−(Z−Zij)22σz2]计算地面浓度时，Z=0，则有：Cz=EXP(−Zij22σz2)当混合层高度Zi不为零时，垂直扩散项分以下几种情况计算。设排放源几何高度为hs，混合层高度为Zi，令Zi'=Zi-hs，设烟气抬升高为Δh（烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气穿透率：p=1.5−Zi'Δh，按不同的P值，分别...
VIP
