《大气环境容量计算模型.doc》由会员分享，可在线阅读，更多相关《大气环境容量计算模型.doc（34页珍藏版）》请在金锄头文库上搜索。
1、附件二：大气环境容量测算模型简介说明：本部分内容是“重点城市大气环境容量核定工作方案”中提到的各推荐模型的简介，主要目的是为了使各城市了解各模型的功能和基本原理，同时，了解如选用该模型，都需要准备哪些输入数据，以便各城市根据本市的实际情况，提前准备。第一部分 大气扩散烟团轨迹模型1 大气扩散烟团轨迹模型简介该模型由国家环境保护总局环境规划院开发。烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。1.1 三维风场的计算首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放
2、烟团的时间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下：式中： 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，其间经
3、过距离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个烟团从源点释放后
4、到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。由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，垂直扩散项分以下几种情况讨论：当混合层高为
5、零时（即无混合层时）有：计算地面浓度时，Z=0，则有：当混合层高度Zi不为零时，垂直扩散项分以下几种情况计算。设排放源几何高度为hs，混合层高度为Zi，令，设烟气抬升高为（烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气穿透率：，按不同的P值，分别计算Cz。当P=0，即时，认为污染物全在混合层内，按封闭性扩散式计算，即污染物在混合层与地面间多次反射。式中：N为反射次数，一般取为N=4即可。当P1时，即时，认为污染物完全穿透混合层，并在混合层以上的稳定层中扩散，由混合层的阻挡而不能到达地面，这时令Cz=0。当0P1000B0.9143700.8650140.281846
6、0.396353010001000BC0.9193250.8750860.2295000.314238010001000C0.9242790.8851570.1771540.232123010001000CD0.9268490.8869400.1439400.189396010001000D0.9294810.8887230.1107260.146669010001000DE0.9251180.8927940.09856310.124308010001000E0.9208180.8968640.0860010.124308010001000F0.9294810.8887230.05536340
7、.073348010001000表 1.4-2 垂直扩散参数幂函数表达式数据扩散参数稳定度等级(PS)2下风距离，mA1.121541.52602.108810.07999040.008547710.0002115450300300500500B0.9410151.093560.1271900.05702510500500BC0.9410151.007700.1146820.07571820500500C0.9175950.1068030CD0.8386280.7564100.8155750.1261520.2356670.1366590200020001000010000D0.8262120
8、.6320230.5553600.1046340.4001670.8107631100010001000010000DEo.7768640.5723470.4991490.1046340.4001671.038100200020001000010000E0.7883700.5651880.4147430.09275290.4333841.732410100010001000010000F0.784400.5259690.3226590.06207650.3700152.406910100010001000010000(3) 丘陵山区的农村或城市，其扩散参数选取方法同工业区。1.4.2 小风和静
9、风（U101.5m/s) 时，0.5h取样时间的扩散参数按表1.4-3选取表1.4-3 小风和静风扩散参数的系数、稳定度（PS）U100.5m/s1.5m/sU100.5m/sU100.5m/s1.5m/sU100.5/sA0.930.760.151.57B0.760.560.470.47C0.550.350.210.21D0.470.270.120.12E0.440.240.070.07F0.440.240.050.051.5 烟气抬升公式1.5.1 有风时，中性和不稳定条件的烟气抬升高度H（m）（1）当烟气热释放率Qh大于或等于是2100KJ/s，且烟气温度与环境温度的差值T大于或等于35
10、K时，H采用下式计算：式中： no烟气热状况及地表系数，见表1.5-1；n1烟气热释放率指数，见表1.5-1；n2排气筒高度指数，见表1.5-1；Qh烟气热释放率，KJ/s；H排气筒距地面几何高度，m ，超过去240m时，取H=240m；Pa大气压力，hPa；Qv实际排烟率，m3/s；T烟气出口温度与环境温度差，K；Ts烟气出口温度，K；Ta环境大气温度，K；U排气筒出口处平均风速，m/s 。表1.5-1 no、n1、n2的选取Qh，KJ/s地表状况（平原）non1n2Qh，KJ/s农村或城市远郊区1.4271/32/3城市及近郊区1.3031/32/32100Qh21000且T35K农村或城
11、市远郊区0.3323/52/5城市及近郊区0.2923/52/5（2）当1700 kJ/sQh2100KJ/s时，式中：Vs排气筒出口处烟气排出速度，m/s；D排气筒出口直径，m；H2按（1）方法计算，no、n1、n2按表1.5-1中Qh值较小的一类选取；Qh ，U与（1）中的定义相同。（3）当Qh1700kJ/s或者T35K时，1.5.2 有风时,稳定条件按下式计算烟气抬升高度H(m)。1.5.3 静风和小风时，按下式计算烟气抬升高度H(m).。但取值不宜小于0.01K/m。2 模型运行所需数据数据文件1：共四行：第1行：X方向网格点的数目（MX），Y方向网格点的数目（MY），Z方向风的观测
12、数据层数目（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；(C1(I,J),I=1,MX),J=1,MY)数据文件3：网格点上的高程（地形值）， 单位m；数据文件4：共六部分。第1部分：污染源数据，共12列，MSC行；序号1 23456789101112
