第三章理想条件下空气污染物散布-扩散参数

1、q(x, y, z; H) Q, u, y , z , H, 扩散参数怎样给出？3.3.1 早期的大气扩散参数处理早期的大气扩散参数处理3.3.1.1 萨顿模式给出RL ()的具体形式,根据Taylor公式求出扩散参数，再代入Gauss公式得到萨顿公式 萨顿扩散参数Cx,Cy时间间隔宏观粘滞度，湍流特征量0222zuNwvu萨顿模式RL ()nLLLvNNRRR2)(0)(, 1)0(2)(vNfRLmzzuummn)(,1211nzznyyxCxC2222222121nnnznnnyuwunnNCuvunnNC12221222)2)(1 (4)2)(1 (43.3.1.2 直接测量湍流量的方

2、法直接测量湍流量的方法1)、风向脉动标准差(H. E. Cramer,1957)由Taylor公式，在近距离，运行时间不长，有当水平风向脉动摆角不大（1/T，它们对扩散都不起作用，只有那些频率非常低的湍涡才对扩散有贡献w 上面的参量都是拉氏形式，按下列假定转化成欧拉形式的参量w 假设1 只要观测v的时间（取样时间）大于粒子运行时间的5倍，就可以用定点观测得到脉动速度方差的代替跟踪观测得到的脉动速度方差。225ELvvT 时，当假定2：拉氏自相关系数与欧拉自相关系数形式完全相同，仅在时间尺度上差了一个倍数。）（）（时，当tRRtEL)(2cos)(42cos)(4)(00nFtdtntRdnRn

3、FELLL欧拉变数表示的湍流时间尺度较小，它的频率就高。由图上可以看出：欧拉谱向高频端有一个位移。 将上式代入Taylor公式的谱形式：dnnTnTnFTvyEy2202222/)/(sin)(/,222)(TvTy下标表示取样时间，T/表示先对v取T/时间的平均以后在计算方差。利用Cramer公式可以写出xxxxxxuxETEvzuxATApy,p、 v分别表示以角度（弧度）表示的横风向和垂直向浓度分布标准差。当扩散物质的运行时间T足够小，有： 时间尺度是唯一的待定参量，起初=0.44/i，值改变时，扩散参数不灵敏。 给定值以后，就可以根据水平风向脉动记录计算各个距离的p，其步骤如下： 1）

4、根据取样时间确定平均时间，最大平均时间sn =/5。 2）用不同的平均时间s1、s2, sn 等对原始的风向脉动记录进行滑动平均，得到经过平滑的n个、风向变化序列s1、s2， sn 。 3）对上述风向变化序列，分别计算标准差，得到(A)s1、(A)s2, (A)sn 。 4）根据结果做(A)s -s曲线图。 5）根据图求出任意y。w 的确定方法：1）根据定义2）由扩散参数的比求原则上有3）利用经验关系求值：值与大气稳定度有关，并与湍流强度成反比，值范围110。大气越稳定，湍流强度越小，值越大。中性层结取4。200)()()()()(vtvtvRdRdRTTvELELTTEvLy)()(TTEv

5、Ey/)()(2.1()()yLyE欧拉时间尺度。）式的拉氏时间尺度和分别代表满足（和1 . 2,ELELTTTTvvyucucicC为实验常数，0.350.80，取0.60。w 适用于开阔平坦地形上的小尺度扩散问题，利用A 和E 计算扩散参数，精确度高。但对观测要求较高，计算量大。3.3.2 稳定度级别与扩散曲线法稳定度级别与扩散曲线法 利用常规气象资料把大气的扩散能力划分为六个稳定度等级，从A到F（极不稳定稳定） 用实验曲线给出每个等级的y和z随距离的变化，这些实验曲线称为扩散曲线。 稳定度的划分方法多种，有差别。3.3.2.1 P-G-T方法 Pasquill把风速和辐射状况作为划分稳定

6、度等级的指标。 辐射状况由白天日照和夜间云量来估计。表3.3表3.3 Pasquill稳定度分级方法 地面风速（m/s）日间日射强度夜间天空状况强中等弱薄云遮天或云量4/8云量3/86CDDDD 不论何种天空状况，夜晚头不论何种天空状况，夜晚头1及最后及最后1h稳定度为中性。稳定度为中性。（日落前1h至日出后1h为夜晚）。 仲夏晴天中午为强日照，寒冬晴天中午为弱日照仲夏晴天中午为强日照，寒冬晴天中午为弱日照。 Turner稳定度划分 首先，依据研究地域的纬度和赤纬、时角定出太阳高度角h sinh=sinsin+coscoscos 其中时角=(t-12)15o ，t为地方时。 其次，按太阳高度角

7、和天空状况确定日射等级，见表3.5，3.6 最后，结合地面风速给出稳定度，表3.7。地面风速 日射等级 43210-1-26CDDDDDD表3.7 Turner的稳定度分类方法时间天空状况日射等级无论日间或夜间 总云量10/10，而且云高4/10-1日间 h15o115oh35o235oh60o4云高为60o1h60o0云高9，不论h0夜间表3.6 日射等级确定规则w Turner的稳定度划分方法较为确切。只要有地面风速、云量、云高常规观测资料，就可以客观地定出稳定度级别。所以这种方法是实际工作中最常应用的方法。此方法原为机场每小时的气象观测资料而设计的，机场多在城郊，下垫面比较平坦，在应用于

8、城市下垫面时应做稳定度等级的订正。w 上述Pasquill和Turner稳定度计算方法，都是用宏观的气象资料间接地给出微观湍流特征。大气扩散速率最终取决于流场的湍流特性，用上述宏观的气象条件划分稳定度等级是一种间接方法。有时，大气宏观条件大致相同，但流场湍流特性会差别很大。因此用湍流场的特征量表征大气稳定度更为有效。a.常规气象资料表征稳定度太粗糙，没有考虑各种大、中尺度过程的影响，也没有考虑下垫面的影响。与其它稳定度分类法相比，差别较大。b.稳定度等级代表不同扩散速率的等级，扩散试验表明，使用定量的稳定度指标有时不能说明扩散速率的快慢。如风向发生缓慢偏转，显示有大尺度湍涡的作用，即使是逆温的

9、情况也不例外。稳定度分类扩散曲线 a.没有考虑地面粗糙度的影响，扩散曲线是依据十分平坦的观测资料得出的，对城市和其它粗糙地形的扩散速率估计过低。需进行订正。b.扩散曲线的取样时间为3分钟，一般短时间平均浓度需要的取样时间为3060分钟，必须订正到相同的时间才可以应用。c.垂直扩散参数的曲线适用于地面源，不一定适用于高架源。P-G适用性p75p763.3.3. 扩散参数法的修改完善扩散参数法的修改完善3.3.3.1 国家标准国家标准 20世纪70年代以来，在我国环保研究实践中积累了很多经验，对扩散曲线不断修改和完善，并结合实际情况制定了国家标准GB/T13201-91。 稳定度、扩散参数确定：w

10、 云量 + 太阳高度角 太阳辐射等级w 太阳辐射等级 + 地面风速 大气稳定度w 由稳定度等级查算扩散参数 有风时（风速u 1.5m/s）扩散参数 小风(0.5 u 1.5m/s)和静风(u 0.5m/s)扩散参数w 注意：根据下垫面条件调整确定稳定度等级；取样时间订正 。p79表3.10-a 太阳辐射等级 云量，1/10太阳辐射等级数总云量/低云量夜间h15o15oh35o3565o4/4-2-1+1+2+357/4-10+1+2+38/4-100+1+15/570000+18/800000表3.10-b 大气稳定度等级地面风速 m/s太阳辐射等级+3+2+10-1-21.9AABBDEF2

11、2.9ABBCDEF34.9BBCCDDE55.9CCDDDDD6DDDDDD表3.11-a 横向扩散参数 幂函数表达式数据扩散参数稳定度等级11下风距离y=1x1A0.9010740.8509340.4258090.602052010001000B0.9143700.8650140.2818460.396353010001000BC0.9193250.8750860.2295000.314238010001000C0.9242790.8851570.1771540.232123010001000CD0.9268490.8869400.1439400.189396010001000D0.929

12、4180.8887230.1107260.14669010001000DE0.9251180.8927940.09856310.124308010001000E0.9208180.8968640.08640010.101947010001000F0.9294180.8887230.05536340.0733348010001000表3.11-c 小风和静风时扩散参数的系数稳定度0102U100.5m/s1.5m/sU100.5m/sU100.5m/s 1.5m/sU100.5m/sA0.930.761.57B0.760.560.47C0.550.350.21D0.470.270.12E0.44

13、0.240.07F0.440.240.05 取样时间的订正取样时间的订正w 取样时间指连续取得观测数据并用以计算平均值的时间长度。对污染物浓度计算，即指大气采样的时间。假设大气湍流场是均匀、平稳的，这是Gauss公式的前提，所以在Guass公式中没有指明平均风速和扩散参数的采样时间及观测高度。但在应用中，为了使各种观测资料和计算结果具有可比性，必须考虑这些因素。w 污染物扩散在铅直方向受地面限制，风向倾角变化范围较小，当取样时间超过几分钟后，铅直风向脉动范围变化不大，所以对铅直向扩散参数可以不做取样时间的订正。水平风向可在360o范围内改变，而表3.11中的扩散参数适用于采样时间为0.5小时的

14、情形，如果采样时间大于或小于这个时段需进行订正。qyy)(1212适用时间范围，hq11000.30.510.2y2、y1分别为取样时间为2、1时的横风向扩散参数。等级A（度）A22.5B17.522.5C12.517.5D7.512.5E3.87.5F3.8(2)温度梯度法 p81(3)PBL湍流参量法823.3.3.3 不同源高、不同下垫面不同源高、不同下垫面 稳定度的高度变化 扩散曲线依据的源高、源特性 不同 z0的影响，p87 表3.233.3.3.4 扩散曲线的内插完善扩散曲线的内插完善 综合，表3.253.3.3.2 其它稳定度分类方法其它稳定度分类方法 (1)风向脉动标准差美国环

15、保局(EPA, 1990)，推荐了用表征湍流强度的风向标准差A， E 划分稳定度的方法。p79表3.25 Briggs扩散参数 (单位：m) (a) 开阔乡间条件稳定度级别yz A0.22x(1+0.0001x)-1/20.20 xB0.16x(1+0.0001x)-1/20.12xC0.11x(1+0.0001x)-1/20.08x(1+0.0002x)-1/2D0.08x(1+0.0001x)-1/20.06x(1+0.0015x)-1/2E0.06x(1+0.0001x)-1/20.03x(1+0.0003x)-1F0.04x(1+0.0001x)-1/20.016x(1+0.0003x

16、)-1(b)城市条件yz A-B0.32x(1+0.0004x)-1/20.24x(1+0.001x)-1/2C0.22x(1+0.0004x)-1/20.20 xD0.16x(1+0.0004x)-1/20.14x(1+0.0003x)-1/2E-F0.11x(1+0.0004x)-1/20.08x(1+0.0015x)-1/23.3.4 风向脉动与扩散函数法将扩散参数y和z与风向脉动标准差A和E联系起来,由此形成了一种由风向脉动与扩散函数确定扩散系数的方法，这是扩散参数研究的主要进展之一。由Taylor公式： tL拉格朗日时间尺度，fy、fz扩散函数当T 很小时，有Taylor公式可得：

17、)(LyvytTfT)(LzwztTfTTxvATxwE A、E可以从风向脉动资料算出，只要知道扩散函数fy、fz和估计tL的方法，就可算出扩散参数。许多外场扩散试验都有扩散参数和A、E的同步观测资料，据此可以反推得出扩散函数fy、fz的形式。w P. R. Draxler汇总了11次试验资料，将扩散资料按源高、稳定度及铅直向或水平向进行分类，分别求出扩散函数。 其中T为扩散时间，Ti为fy或fz等于0.5所需的时间。这种方法对于不稳定时地面源的铅直扩散和稳定时高架源的铅直扩散是不适合的。w Pasquill 5 . 0)(90. 011izyTTff)()()(LeyALeyevLyvytu

18、xxftuxfuxtTTf)(xxfAy汇总扩散函数表达式，表汇总扩散函数表达式，表3.27 p953.3.5 扩散参数研究状况(1) 实用体系(2) 观测试验(3) 施放高度对扩散参数的影响(4) PBL结构及湍流参量与扩散(1) 实用体系a. 稳定度扩散级别与扩散曲线法 b. 风向脉动与扩散函数法c. 对上述两种确定扩散参数的实用体系有如下实施意见wy 1）10km以内，研究透彻，资料多，数据可靠，不同扩散参数结果彼此接近。 2）最佳公式： 3）y可用P-G曲线做粗略估计，但应根据粗糙度和取样时间进行订正)(xfxAywz 1) 研究少，测量数据少，可靠性差，与y比，更多地依赖于温度层结。 2) P-G曲线的z根据近地面源1km以内的扩散试验资料得出，其它都是外推的结果，故P-G曲线适用于近地面源的近距离扩散。 wz 3)一般情况下应考虑z0的影响，只