mapinfo-大气环境容量测算模型

附件二：

大气环境容量测算模型简介

说明：本部分内容是“重点城市大气环境容量核定工作方案”中提到的各推荐模型的简介，主要目的是为了使各城市了解各模型的功能和基本原理，同时，了解如选用该模型，都需要准备哪些输入数据，以便各城市根据本市的实际情况，提前准备。

第一部分大气扩散烟团轨迹模型

1大气扩散烟团轨迹模型简介

该模型由国家环境保护总局环境规划院开发。

烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。

1.1三维风场的计算

首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放烟团的时间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下：

V i =v（t i） V（t2）-V（t i）l丄

n

n二花-鮎?:t

式中：V(t i )、V(t 2)—分别为第1和第2个观测时刻的风场值;

=t —烟团释放时间步长；

n —为t i 、t 2间隔内的时间步长数目；

V i —表示t i 、t 2间隔内第i 个时间步长上的风场值

1.2烟团轨迹的计算

位于源点的某污染源，在t o 时刻释放出第1个烟团，此烟团按t o 时刻源点处 的风向风速运行，经一个时间步长

t 后在t i 时刻到达P ii ，经过的距离为D ii ，

从t i 开始，第一个烟团按P ii 处t i 时刻的风向风速走一个时间步长，在

t 2时刻到

达P i2,其间经过距离D i2，与此同时，在t i 时刻从源点释放出第2个烟团，按 源点处t i 时刻的风向风速运行，在t 2时刻到达P 22，其经过的距离为D 22，以此 类推，从t o 时刻经过

j 个氏，到t j 时刻共释放出了 j 个烟团，这时，这j 个烟团 的中心分别位于Pij ，i=i , 2,…j ,设源的坐标为(Xs , Ys ，Zs(t))，Zs(t)为t 时 刻烟团的有效抬升高度，Pij 的坐标为

(Xij , Yij , Zij ), u 、v 分别为风速在X 、 丫方向的分量，则有如下计算公式：

t i 时刻：

Xu 二 X s U[t °,X s ,Y s ,Z s (t 。)]凤 Y ii =Y s V[t °,X s ,Y s ,Z s (t 。)] Zu =Z s W[t °,X s ,Y s ,Z s (t 。)] D i =D ii 「(Xu 二X s )2

—(丫ii 二Y s )2

t 2时刻：

X i2 =X ii U[t i ,X ii ,Y i , Z ii ]厶 t Y i2 =Y ii ?V[t i ,X ii ,Y ii ,Z ii ]

Z

i2 - Z i i

W[t

i

, X ii ,Y ii , Z ii ] -A

D i 2 =D ii +D i2

+J(X i2 —Xu)2 +(Y i2 —Y ii )2

X 22 二X s U[t i ,X s ,Y s ,Z s (t i )] ：t Y 22 二Y s V[t i ,X s ,Y s ,Z s (t i )] ：t

Z 22 二Z s W[t i ,X s ,Y s ,Z s (t i )]

D ；二 D 22 = ..(X 22 —X s )2 (Y 22 —Y s )2

以此类推，至y tj 时刻，共释放出j 个烟团，这些烟团最后的中心位置分别在

Pij ，Xij ，Yij ，Zij ，i=1，2,…j ,对于第 i 个烟团有：

=X

i(jJ)

' U [t j J ，X i (j d)，Y i(j d)，Z i( j d)K'：

t

Y ij =Y i(jd)

' V [t

j J ，X

i(j J)，

Y

i( j J)

，

Z

i(j J)

K-'：t

Z

ij - Z i( j J)

' W[t jJ ，X i(jJ)，Y i(jd)，Z i(jJ)]川

j

_____________________________________________

D/ 八 D ik 二 D i jJ

■ ,(X j -X i (j d

))2 - (Y j -丫心二))2

k 4

D i j 为i 个烟团从源点释放后到tj 时刻所经过的距离。

1.3浓度公式

由前一个小节的计算，已找到由S 点(Xs ，Ys)的污染源释放出来的所有烟团 在第j 个时刻所处的位置，这样S 处的污染源在第j 个时刻在地面某接受点R(X 、 丫、0)处造成的浓度就是所有i 个烟团的浓度贡献之和。考虑中心位于Pij 的烟团 对R 点的浓度贡献，则有：

C i 二

2

Qs

C X C Y C Z C b C d

（2兀）6巧6

式中：Qs —源强，mg/s ；

二y 、二z : — X 方向、丫方向、Z 方向的大气扩散参数，m ；

Cx 、Cy 、Cz ：— X 、丫、Z 方向扩散项，Cz 在后面给出算式； C b 为污染物转化项，b 为转化率，1/s ； C d 为污染物沉降项，V d 为沉降速率，m/s 。

由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，

垂直扩散

C X =EXP

「

（X-X ij ）2

〕 2cr 2

J X

2 EXP 一（Y_Y j ）

C b

C d 2CT 2 厶y

-EXP -b j 4

,（v d iN ）2 =

项分以下几种情况讨论：

当混合层高为零时（即无混合层时）有:

计算地面浓度时，Z=0，则有: Z 2 j

C z 二 EXP （ 4）

2f

当混合层高度Zi 不为零时，垂直扩散项分以下几种情况计算。

设排放源几何高度为hs ,混合层高度为Zi ，令Z ；= Z j -hs ，设烟气抬升高 为丄h （烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气

I

穿透率：卩=1.5-色，按不同的P 值，分别计算Cz 0

Ah

2 '

当P=0，即< 2Z i 时，认为污染物全在混合层内，按封闭性扩散式计算,

3

即污染物在混合层与地面间多次反射。

式中：N —为反射次数，一般取为 N=4即可。

当P > 1时，即-：h 2Z ；时，认为污染物完全穿透混合层，并在混合层以 上的稳定层中扩散，由混合层的阻挡而不能到达地面，这时令

Cz=0o

2

当0

C ,EXP [(^L EXP

? 疋 一

-(Z-Z ij )2 31

「加 z 2

N

C z

EXP[

n — -N

(Z j -2nZ i )2

乙2

C z1 = P EXP(打)

2b z

而(1-P )部分的烟团在Zi 处按封闭扩散:

(Z j - 2n Z i )2]

1.4大气扩散参数

(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

横向扩散参数幕函数表达式数据

N

C z2=(1—P)' EXP[-

n

1.4.1有风时扩散参数c

y 、

(T z 的确定(0.5h 取样时间)

表垂直扩散参数幕函数表达式数据

扩散参数

稳定度等级 (P ? S) a 2

L

下风距离，m

1.12154 0.0799904 0~300

A

1.5260 0.00854771 300~500

2.10881 0.000211545 >500

B 0.941015 0.127190 0~500

1.09356 0.0570251 >500

B~C

0.941015 0.114682 0~500

1.00770 0.0757182 >500

C 0.917595 0.106803 0

0.838628 0.126152 0~2000

C~D

0.756410 0.235667 2000~10000

0.815575 0.136659 >10000

0.826212 0.104634 1~1000

D

0.632023 0.400167 1000~10000

0.555360 0.810763 >10000

O.776864 0.104634 0~2000

D~E

0.572347 0.400167 2000~10000

0.499149 1.03810 >10000

0.788370 0.0927529 0~1000

E

0.565188 0.433384 1000~10000

0.414743 1.73241 >10000

0.78440 0.0620765 0~1000

F

0.525969 0.370015 1000~10000

0.322659

2.40691

>10000

(3) 丘陵山区的农村或城市，其扩散参数选取方法同工业区

1.4.2小风和静风(U io V 1.5m/s)时，0.5h 取样时间的扩散参数按表1.4-3选取

表1.4-3小风和静风扩散参数的系数 01、

02

稳定度(P ? S )

%

%2

U 10 v 0.5m/s

1.5m/s > U 10》

0.5m/s U 10v 0.5m/s

1.5m/s > U 10》

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

x

y

=1

01

> 二

z

1.5烟气抬升公式

1.5.1有风时，中性和不稳定条件的烟气抬升高度厶H (m)

(1)当烟气热释放率Q h大于或等于是2100KJ/S,且烟气温度与环境温度的差值△ T 大于或等于35K时，△ H采用下式计算：

H n o Q“H n2U J

Q h = 0.35P：.Q.—

T讥-「.

式中：no——烟气热状况及地表系数，见表1.5-1；

n1――烟气热释放率指数，见表1.5-1；

n2 ------- 排气筒高度指数，见表1.5-1 ；

Qh——烟气热释放率，KJ/s；

H――排气筒距地面几何高度，m，超过去240m时，取H=240m；P a

----------- 大气压力，hP a；

Q v ------ 实际排烟率，m%；

△ T――烟气出口温度与环境温度差，K ；

T s――烟气出口温度，K ；

T a――环境大气温度，K ；

U ――排气筒出口处平均风速，m/s。

。、i、 2 的选取

表

(2)当1700 kJ/s v Q h v 2100KJ/s 时，

Qh 1700

■V. - ?氷2Y门

400

■ ：^2 1.5V s D 0.01Q h /U -0.048 Q h T700 /U