大气环境容量计算模型-A-P值法在大气环境容量测算
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
0.124308
0~1000
>1000
F
0.929481
0.888723
0.0553634
0.073348
0~1000
>1000
表1.4-2垂直扩散参数幂函数表达式数据
扩散参数
稳定度等级
(P·S)
α2
下风距离，m
A
1.12154
1.5260
2.10881
0.0799904
0.00854771
0.000211545
2．按功能分区的控制浓度（标准年平均浓度限值）Ci
3．确定各个功能区总量控制系数Ai值
4．确定各个功能区允许排放总量：
5.根据总量控制区所在地区，按GB/T13201-91表1查取低源分担率值,确定各个功能区低矮源(面源)允许排放总量：
6．计算总量控制区允许排放总量 和低矮面源允许排放总量
，
7．如果计算出的 值小于上级部门的指令允许排放总量，则在总量控制区内就使用该 值可以继续采用A-P值法确定总量控制区内各个功能分区内的点源允许排放量，也可以在该市的辖区内适当增加控制区面积（即增加新的开发区）以使A-P值法计算的 值与指令总量接近，但是不得超过指令值。
0.886940
0.143940
0.189396
0~1000
>1000
D
0.929481
0.888723
0.110726
0.146669
0~1000
>1000
D~E
0.925118
0.892794
0.0985631
0.124308
0~1000
>1000
E
0.920818
0.896864
0.086001
1.03810
0~2000
2000~10000
>10000
E
0.788370
0.565188
0.414743
0.0927529
0.433384
1.73241
0~1000
1000~10000
>10000
F
0.78440
0.525969
0.322659
0.0620765
0.370015
2.40691
0~1000
t1时刻：
t2时刻：
以此类推，到tj时刻，共释放出j个烟团，这些烟团最后的中心位置分别在Pij，Xij，Yij，Zij，i=1，2，…j，对于第i个烟团有：
为i个烟团从源点释放后到tj时刻所经过的距离。
由前一个小节的计算，已找到由S点(Xs，Ys)的污染源释放出来的所有烟团在第j个时刻所处的位置，这样S处的污染源在第j个时刻在地面某接受点R(X、Y、0)处造成的浓度就是所有i个烟团的浓度贡献之和。考虑中心位于Pij的烟团对R点的浓度贡献，则有：
(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
0~300
300~500
>500
B
0.941015
1.09356
0.127190
0.0570251
0~500
>500
B~C
0.941015
1.00770
0.114682
0.0757182
0~500
>500
C
0.917595
0.106803
0
C~D
0.838628
0.756410
0.815575
0.126152
当P=0，即 时，认为污染物全在混合层内，按封闭性扩散式计算，即污染物在混合层与地面间多次反射。
式中：N－为反射次数，一般取为N=4即可。
当P＞1时，即 时，认为污染物完全穿透混合层，并在混合层以上的稳定层中扩散，由混合层的阻挡而不能到达地面，这时令Cz=0。
当0<P<1，即 时认为是部分穿透情形，这时有部分污染物抬升到混合层以上，而（1-P）部分被封闭在混合层以内，Cz按下式计算：
许多文献认为穿透到混合层以上的污染物被阻挡后不能向地面扩散，当地区大气层结处于中性偏稳定结构时，混合层对污染物的阻挡作用并不是很强，这时可设计成让这部分烟团在 高度上向下扩散，则有：
而（1-P）部分的烟团在Zi处按封闭扩散：
(1)平原地区农村及城市远郊区的扩散参数选取方法如下：
A、B、C级稳定度直接由表1.4-1和表1.4-2查算，D、E、F级稳定度则需向不稳定方向提半级后由表1.4-1和表1.4-2查算。
城市及近郊区
0.292
3/5
2/5
（2）当1700kJ/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时，
式中：Qs—源强，mg/s；
：—X方向、Y方向、Z方向的大气扩散参数，m；
Cx、Cy、Cz：—X、Y、Z方向扩散项，Cz在后面给出算式；
Cb为污染物转化项，b为转化率，1/s；
Cd为污染物沉降项，Vd为沉降速率，m/s。
由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，垂直扩散项分以下几种情况讨论：
0~1000
>1000
B
0.914370
0.865014
0.281846
0.396353
0~1000
>1000
B~C
0.919325
0.875086
0.229500
0.314238
0~1000
>1000
C
0.924279
0.885157
0.177154
0.232123
0~1000
>1000
C~D
0.926849
如果 ，则取
8．再按以下公式计算各个功能区内所有点源的最终允许排放量（吨/小时）
四、列出总量控制区和各个功能分区的允许排放总量、低矮面源允许排放总量和各个点源的允许排放量清单
五、将城市大气污染物排放总量控制A-P值法试用于控制pm10时，低源分担率值在长江以北可适当放宽到0.4,在长江以南可适当放宽到0.5
当混合层高为零时（即无混合层时）有：
计算地面浓度时，Z=0，则有：
当混合层高度Zi不为零时，垂直扩散项分以下几种情况计算。
设排放源几何高度为hs，混合层高度为Zi，令 ，设烟气抬升高为 （烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气穿透率： ，按不同的P值，分别计算Cz。
城市大气污染物（以SO2为例）排放总量控制A-P值法简介
一、采用A-P值法确定总量控制区允许排放总量时所需的资料
1．总量控制区面积S
2．总量控制区内的功能分区的面积Si
3．功能分区的控制浓度（标准浓度限值）Ci
二、采用A-P值法确定总量控制区允许排放总量（万吨/年）的步骤
1．根据总量控制区所在地区，按GB/T13201-91表1查取总量控制系数A值（取中值）
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；
2．按以下公式计算各个功能区内所有点源的初始允许排放量（吨/小时）
-点源的有效高度（烟囱的实体高度加上抬升高度）,这里Ci用标准日平均浓度限值
3．按点源的实体高度分类为低架源（排气筒高度小于30米）、中架源（排气筒高度大于或等于30米但小于100米）、高架源（排气筒高度大于或等于100米）
4．在功能分区内，将属于中架源的点源初始允许排放量相加，并乘以8760小时得到中架源的年初始允许排放总量Qmi，并用（万吨/年）表示
1000~10000
>10000
(3)丘陵山区的农村或城市，其扩散参数选取方法同工业区。
表1.4-3小风和静风扩散参数的系数 、
稳定度（P·S）
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
第一部分
1
该模型由国家环境保护总局环境规划院开发。
烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。
首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放烟团的时间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下：
但 取值不宜小于0.01K/m。
2
共四行：
第1行：
X方向网格点的数目（MX），Y方向网格点的数目（MY），Z方向风的观测数据层数目（MZ），最大有效烟团数（NT，默认110），污染源数目（MSC），气象观测小时数目（NTimes），稳定度数目（NeleTa，默认24）；
第2行：
烟团的时间步长（分，默认30.0）
8．如果计算出的 值大于上级部门的指令允许排放总量，则在总量控制区内用下式计算出A值后，再从本节第3条向下继续计算。
用该 值可以继续采用A-P值法确定总量控制区内各个功能分区内的点源允许排放量，
三、采用A-P值法确定总量控制区内各个功能分区内的点源允许排放量的步骤
1．根据总量控制区所在地区，按GB/T13201-91表1查取总量控制系数P值
0.235667
0.136659
0~2000
2000~10000
>10000
D
0.826212
0.632023
0.555360
0.104634
0.400167
0.810763
1~1000
1000~10000
>10000
D~E
o.776864
0.572347
0.499149
0.104634
0.400167
城市大气污染物（以SO2为例）排放总量控制A-P值法简介结束。
附件二：
大气环境容量测算模型简介
说明：本部分内容是“重点城市大气环境容量核定工作方案”中提到的各推荐模型的简介，主要目的是为了使各城市了解各模型的功能和基本原理，同时，了解如选用该模型，都需要准备哪些输入数据，以便各城市根据本市的实际情况，提前准备。
第3行：
X方向步长（m，默认1000.0），Y方向步长（m，默认1000.0），大气压力(hPa，默认1013.25)，规划区类型（1农村，2城市）；
第4行：
化学转化率1/s，沉降速率m/s；
网格点上的背景浓度值，单位mg/m3；
Ta——环境大气温度，K；
U——排气筒出口处平均风速，m/s。
表1.5-1no、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.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）当烟气热释放率Qh大于或等于是2100KJ/s，且烟气温度与环境温度的差值△T大于或等于35K时，△H采用下式计算：
式中：no——烟气热状况及地表系数，见表1.5-1；
5．计算各个功能分区内的点源调整系数i
如果 ，则取
6．在总量控制区内，将属于中架源的点源初始排放量相加，并乘以8760小时得到中架源的年初始允许排放总量Qm；将属于高架源的点源初始排放量相加，并乘以8760小时得到高架源的年初始允许排放总量Qc,二者都用（万吨/年）表示。
7．计算总量控制区内的点源调整系数
式中：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方向的分量，则有如下计算公式：
0~1000
>1000
F
0.929481
0.888723
0.0553634
0.073348
0~1000
>1000
表1.4-2垂直扩散参数幂函数表达式数据
扩散参数
稳定度等级
(P·S)
α2
下风距离，m
A
1.12154
1.5260
2.10881
0.0799904
0.00854771
0.000211545
2．按功能分区的控制浓度（标准年平均浓度限值）Ci
3．确定各个功能区总量控制系数Ai值
4．确定各个功能区允许排放总量：
5.根据总量控制区所在地区，按GB/T13201-91表1查取低源分担率值,确定各个功能区低矮源(面源)允许排放总量：
6．计算总量控制区允许排放总量 和低矮面源允许排放总量
，
7．如果计算出的 值小于上级部门的指令允许排放总量，则在总量控制区内就使用该 值可以继续采用A-P值法确定总量控制区内各个功能分区内的点源允许排放量，也可以在该市的辖区内适当增加控制区面积（即增加新的开发区）以使A-P值法计算的 值与指令总量接近，但是不得超过指令值。
0.886940
0.143940
0.189396
0~1000
>1000
D
0.929481
0.888723
0.110726
0.146669
0~1000
>1000
D~E
0.925118
0.892794
0.0985631
0.124308
0~1000
>1000
E
0.920818
0.896864
0.086001
1.03810
0~2000
2000~10000
>10000
E
0.788370
0.565188
0.414743
0.0927529
0.433384
1.73241
0~1000
1000~10000
>10000
F
0.78440
0.525969
0.322659
0.0620765
0.370015
2.40691
0~1000
t1时刻：
t2时刻：
以此类推，到tj时刻，共释放出j个烟团，这些烟团最后的中心位置分别在Pij，Xij，Yij，Zij，i=1，2，…j，对于第i个烟团有：
为i个烟团从源点释放后到tj时刻所经过的距离。
由前一个小节的计算，已找到由S点(Xs，Ys)的污染源释放出来的所有烟团在第j个时刻所处的位置，这样S处的污染源在第j个时刻在地面某接受点R(X、Y、0)处造成的浓度就是所有i个烟团的浓度贡献之和。考虑中心位于Pij的烟团对R点的浓度贡献，则有：
(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
0~300
300~500
>500
B
0.941015
1.09356
0.127190
0.0570251
0~500
>500
B~C
0.941015
1.00770
0.114682
0.0757182
0~500
>500
C
0.917595
0.106803
0
C~D
0.838628
0.756410
0.815575
0.126152
当P=0，即 时，认为污染物全在混合层内，按封闭性扩散式计算，即污染物在混合层与地面间多次反射。
式中：N－为反射次数，一般取为N=4即可。
当P＞1时，即 时，认为污染物完全穿透混合层，并在混合层以上的稳定层中扩散，由混合层的阻挡而不能到达地面，这时令Cz=0。
当0<P<1，即 时认为是部分穿透情形，这时有部分污染物抬升到混合层以上，而（1-P）部分被封闭在混合层以内，Cz按下式计算：
许多文献认为穿透到混合层以上的污染物被阻挡后不能向地面扩散，当地区大气层结处于中性偏稳定结构时，混合层对污染物的阻挡作用并不是很强，这时可设计成让这部分烟团在 高度上向下扩散，则有：
而（1-P）部分的烟团在Zi处按封闭扩散：
(1)平原地区农村及城市远郊区的扩散参数选取方法如下：
A、B、C级稳定度直接由表1.4-1和表1.4-2查算，D、E、F级稳定度则需向不稳定方向提半级后由表1.4-1和表1.4-2查算。
城市及近郊区
0.292
3/5
2/5
（2）当1700kJ/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时，
式中：Qs—源强，mg/s；
：—X方向、Y方向、Z方向的大气扩散参数，m；
Cx、Cy、Cz：—X、Y、Z方向扩散项，Cz在后面给出算式；
Cb为污染物转化项，b为转化率，1/s；
Cd为污染物沉降项，Vd为沉降速率，m/s。
由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，垂直扩散项分以下几种情况讨论：
0~1000
>1000
B
0.914370
0.865014
0.281846
0.396353
0~1000
>1000
B~C
0.919325
0.875086
0.229500
0.314238
0~1000
>1000
C
0.924279
0.885157
0.177154
0.232123
0~1000
>1000
C~D
0.926849
如果 ，则取
8．再按以下公式计算各个功能区内所有点源的最终允许排放量（吨/小时）
四、列出总量控制区和各个功能分区的允许排放总量、低矮面源允许排放总量和各个点源的允许排放量清单
五、将城市大气污染物排放总量控制A-P值法试用于控制pm10时，低源分担率值在长江以北可适当放宽到0.4,在长江以南可适当放宽到0.5
当混合层高为零时（即无混合层时）有：
计算地面浓度时，Z=0，则有：
当混合层高度Zi不为零时，垂直扩散项分以下几种情况计算。
设排放源几何高度为hs，混合层高度为Zi，令 ，设烟气抬升高为 （烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气穿透率： ，按不同的P值，分别计算Cz。
城市大气污染物（以SO2为例）排放总量控制A-P值法简介
一、采用A-P值法确定总量控制区允许排放总量时所需的资料
1．总量控制区面积S
2．总量控制区内的功能分区的面积Si
3．功能分区的控制浓度（标准浓度限值）Ci
二、采用A-P值法确定总量控制区允许排放总量（万吨/年）的步骤
1．根据总量控制区所在地区，按GB/T13201-91表1查取总量控制系数A值（取中值）
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；
2．按以下公式计算各个功能区内所有点源的初始允许排放量（吨/小时）
-点源的有效高度（烟囱的实体高度加上抬升高度）,这里Ci用标准日平均浓度限值
3．按点源的实体高度分类为低架源（排气筒高度小于30米）、中架源（排气筒高度大于或等于30米但小于100米）、高架源（排气筒高度大于或等于100米）
4．在功能分区内，将属于中架源的点源初始允许排放量相加，并乘以8760小时得到中架源的年初始允许排放总量Qmi，并用（万吨/年）表示
1000~10000
>10000
(3)丘陵山区的农村或城市，其扩散参数选取方法同工业区。
表1.4-3小风和静风扩散参数的系数 、
稳定度（P·S）
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
第一部分
1
该模型由国家环境保护总局环境规划院开发。
烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。
首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放烟团的时间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下：
但 取值不宜小于0.01K/m。
2
共四行：
第1行：
X方向网格点的数目（MX），Y方向网格点的数目（MY），Z方向风的观测数据层数目（MZ），最大有效烟团数（NT，默认110），污染源数目（MSC），气象观测小时数目（NTimes），稳定度数目（NeleTa，默认24）；
第2行：
烟团的时间步长（分，默认30.0）
8．如果计算出的 值大于上级部门的指令允许排放总量，则在总量控制区内用下式计算出A值后，再从本节第3条向下继续计算。
用该 值可以继续采用A-P值法确定总量控制区内各个功能分区内的点源允许排放量，
三、采用A-P值法确定总量控制区内各个功能分区内的点源允许排放量的步骤
1．根据总量控制区所在地区，按GB/T13201-91表1查取总量控制系数P值
0.235667
0.136659
0~2000
2000~10000
>10000
D
0.826212
0.632023
0.555360
0.104634
0.400167
0.810763
1~1000
1000~10000
>10000
D~E
o.776864
0.572347
0.499149
0.104634
0.400167
城市大气污染物（以SO2为例）排放总量控制A-P值法简介结束。
附件二：
大气环境容量测算模型简介
说明：本部分内容是“重点城市大气环境容量核定工作方案”中提到的各推荐模型的简介，主要目的是为了使各城市了解各模型的功能和基本原理，同时，了解如选用该模型，都需要准备哪些输入数据，以便各城市根据本市的实际情况，提前准备。
第3行：
X方向步长（m，默认1000.0），Y方向步长（m，默认1000.0），大气压力(hPa，默认1013.25)，规划区类型（1农村，2城市）；
第4行：
化学转化率1/s，沉降速率m/s；
网格点上的背景浓度值，单位mg/m3；
Ta——环境大气温度，K；
U——排气筒出口处平均风速，m/s。
表1.5-1no、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.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）当烟气热释放率Qh大于或等于是2100KJ/s，且烟气温度与环境温度的差值△T大于或等于35K时，△H采用下式计算：
式中：no——烟气热状况及地表系数，见表1.5-1；
5．计算各个功能分区内的点源调整系数i
如果 ，则取
6．在总量控制区内，将属于中架源的点源初始排放量相加，并乘以8760小时得到中架源的年初始允许排放总量Qm；将属于高架源的点源初始排放量相加，并乘以8760小时得到高架源的年初始允许排放总量Qc,二者都用（万吨/年）表示。
7．计算总量控制区内的点源调整系数
式中：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方向的分量，则有如下计算公式：