理想条件下空大气污染物扩散参数Briggs参数

q(x,y,z;H)~Q,u,y,z,H,……扩散参数怎样给出？3.3.1早期的大气扩散参数处理3.3.1.1萨顿模式给出RL()的具体形式,根据Taylor公式求出扩散参数，再代入Gauss公式得到萨顿公式萨顿扩散参数Cx,Cy时间间隔宏观粘滞度，，湍流特征量0222zuNwvu萨顿模式RL()nLLLvNNRRR2)(0)(,1)0(2)(vNfRLmzzuummn)(,1211nzznyyxCxC2222222121nnnznnnyuwunnNCuvunnNC12221222)2)(1(4)2)(1(43.3.1.2直接测量湍流量的方法1)、风向脉动标准差(H.E.Cramer,1957)由Taylor公式，在近距离，运行时间不长，有当水平风向脉动摆角不大（20o）,有2222221)(uxvTvRyL22222tantanAAuvuuvuuvAA222222xxEzAyCramer模式的基本表达式表3.1,p67采用水平风向摆动角的范围作为特征量，导出扩散参数，并以高架源（108m）扩散试验为依据，得出了相应的经验表达式和扩散曲线。见表3.2,p682)BNL(M.E.Smith,1951)方法，或ASME方法Taylor公式的谱函数形式经过时间T后，在x轴向距离uT的位置上，y向扩散范围不仅与横向湍流强度有关，也与拉氏湍流谱有关。粒子扩散是由不同频率的脉动速度引起的，高频脉动是小湍涡活动的结果，低频脉动是大湍涡活动的结果。扩散是各种不同尺度湍涡活动的结果。拉格朗日参量欧拉参量02cos)()(dnnnFRLLdnnTnTnFTvdtddnnnFvyLLtLTLy22022000222sin)(2cos)(23)J.S.Hay和F.Pasquill方法（1）T足够小，脉动速度的频率n不可能很高，T小到一定程度后，总有：各种频率的速度分量都按其所具有的权重作用于粒子扩散。（2）任意T：频率高于1/T的湍涡（小湍涡）对扩散几乎没有贡献，低于0.1/T的湍涡（大湍涡）不受影响，按其自身的权重作用于扩散。介于两者之间的频率分量虽有贡献，但受到一定程度的削弱。1)()(sin22nTnT22220222sin)(TvdnnTnTnFTvLyTnTnnTnT1.0110)()(sin22(3)T很大时，1/T很小，在构成v’的全部分量中，大都满足n1/T，它们对扩散都不起作用，只有那些频率非常低的湍涡才对扩散有贡献上面的参量都是拉氏形式，按下列假定转化成欧拉形式的参量假设1只要观测v’的时间（取样时间）大于粒子运行时间的5倍，就可以用定点观测得到脉动速度方差的代替跟踪观测得到的脉动速度方差。225ELvvT时，当假定2：拉氏自相关系数与欧拉自相关系数形式完全相同，仅在时间尺度上差了一个倍数。）（）（时，当tRRtEL)(2cos)(42cos)(4)(00nFtdtntRdnRnFELLL欧拉变数表示的湍流时间尺度较小，它的频率就高。由图上可以看出：欧拉谱向高频端有一个位移。将上式代入Taylor公式的谱形式：dnnTnTnFTvyEy2202222/)/(sin)(/,222)(TvTy下标表示取样时间，T/表示先对v’取T/时间的平均以后在计算方差。利用Cramer公式可以写出xxxxxxuxETEvzuxATApy,,,,p、v分别表示以角度（弧度）表示的横风向和垂直向浓度分布标准差。当扩散物质的运行时间T足够小，有：时间尺度是唯一的待定参量，起初=0.44/i，值改变时，扩散参数不灵敏。给定值以后，就可以根据水平风向脉动记录计算各个距离的p，其步骤如下：1）根据取样时间确定平均时间，最大平均时间sn=τ/5。2）用不同的平均时间s1、s2,……sn等对原始的风向脉动记录进行滑动平均，得到经过平滑的n个、风向变化序列s1、s2，……sn。3）对上述风向变化序列，分别计算标准差，得到(A)s1、(A)s2,……(A)sn。4）根据结果做(A)s-s曲线图。5）根据图求出任意y。的确定方法：1）根据定义2）由扩散参数的比求原则上有3）利用经验关系求值：值与大气稳定度有关，并与湍流强度成反比，值范围1~10。大气越稳定，湍流强度越小，值越大。中性层结取4。200)()()()()(vtvtvRdRdRTTvELELTTEvLy])[()(TTEvEy/])[()(2.1()()yLyE欧拉时间尺度。）式的拉氏时间尺度和分别代表满足（和1.2,,ELELTTTTvvyucucicC为实验常数，0.35~0.80，取0.60。适用于开阔平坦地形上的小尺度扩散问题，利用A和E计算扩散参数，精确度高。但对观测要求较高，计算量大。3.3.2稳定度级别与扩散曲线法利用常规气象资料把大气的扩散能力划分为六个稳定度等级，从A到F（极不稳定—稳定）用实验曲线给出每个等级的y和z随距离的变化，这些实验曲线称为扩散曲线。稳定度的划分方法多种，有差别。3.3.2.1P-G-T方法Pasquill把风速和辐射状况作为划分稳定度等级的指标。辐射状况由白天日照和夜间云量来估计。表3.3表3.3Pasquill稳定度分级方法地面风速（m/s）日间日射强度夜间天空状况强中等弱薄云遮天或云量≤4/8云量≤3/82AA-BB2-3A-BBCEF3-5BB-CCDE5-6CC-DDDD6CDDDD不论何种天空状况，夜晚头1及最后1h稳定度为中性。（日落前1h至日出后1h为夜晚）。仲夏晴天中午为强日照，寒冬晴天中午为弱日照。Turner稳定度划分首先，依据研究地域的纬度和赤纬、时角定出太阳高度角hsinh=sinsin+coscoscos其中时角=(t-12)15o，t为地方时。其次，按太阳高度角和天空状况确定日射等级，见表3.5，3.6最后，结合地面风速给出稳定度，表3.7。地面风速日射等级43210-1-22AA-BBCDEF2~3A-BBCDDD-EE3~5BB-CCDDDD-E5~6CCDDDDD-E6CDDDDDD表3.7Turner的稳定度分类方法时间天空状况日射等级无论日间或夜间总云量10/10，而且云高2000m0总云量≤4/10-2总云量4/10-1日间h15o115oh35o235oh60o3h60o4云高为2000m的低云量6~9，且：h60o1h≤60o0云高2000米m的低云量9，不论h0夜间表3.6日射等级确定规则Turner的稳定度划分方法较为确切。只要有地面风速、云量、云高常规观测资料，就可以客观地定出稳定度级别。所以这种方法是实际工作中最常应用的方法。此方法原为机场每小时的气象观测资料而设计的，机场多在城郊，下垫面比较平坦，在应用于城市下垫面时应做稳定度等级的订正。上述Pasquill和Turner稳定度计算方法，都是用宏观的气象资料间接地给出微观湍流特征。大气扩散速率最终取决于流场的湍流特性，用上述宏观的气象条件划分稳定度等级是一种间接方法。有时，大气宏观条件大致相同，但流场湍流特性会差别很大。因此用湍流场的特征量表征大气稳定度更为有效。a.常规气象资料表征稳定度太粗糙，没有考虑各种大、中尺度过程的影响，也没有考虑下垫面的影响。与其它稳定度分类法相比，差别较大。b.稳定度等级代表不同扩散速率的等级，扩散试验表明，使用定量的稳定度指标有时不能说明扩散速率的快慢。如风向发生缓慢偏转，显示有大尺度湍涡的作用，即使是逆温的情况也不例外。稳定度分类扩散曲线a.没有考虑地面粗糙度的影响，扩散曲线是依据十分平坦的观测资料得出的，对城市和其它粗糙地形的扩散速率估计过低。需进行订正。b.扩散曲线的取样时间为3分钟，一般短时间平均浓度需要的取样时间为30~60分钟，必须订正到相同的时间才可以应用。c.垂直扩散参数的曲线适用于地面源，不一定适用于高架源。P-G适用性p75p763.3.3.扩散参数法的修改完善3.3.3.1国家标准20世纪70年代以来，在我国环保研究实践中积累了很多经验，对扩散曲线不断修改和完善，并结合实际情况制定了国家标准GB/T13201-91。稳定度、扩散参数确定：云量+太阳高度角太阳辐射等级太阳辐射等级+地面风速大气稳定度由稳定度等级查算扩散参数有风时（风速u1.5m/s）扩散参数小风(0.5u1.5m/s)和静风(u0.5m/s)扩散参数注意：根据下垫面条件调整确定稳定度等级；取样时间订正。p79表3.10-a太阳辐射等级云量，1/10太阳辐射等级数总云量/低云量夜间h15o15oh35o35h65oh65o4/4-2-1+1+2+35~7/4-10+1+2+38/4-100+1+15/5~70000+18/800000表3.10-b大气稳定度等级地面风速m/s太阳辐射等级+3+2+10-1-21.9AA~BBDEF2~2.9A~BBCDEF3~4.9BB~CCDDE5~5.9CC~DDDDD6DDDDDD表3.11-a横向扩散参数幂函数表达式数据扩散参数稳定度等级11下风距离y=1x1A0.9010740.8509340.4258090.6020520~10001000B0.9143700.8650140.2818460.3963530~10001000B~C0.9193250.8750860.2295000.3142380~10001000C0.9242790.8851570.1771540.2321230~10001000C~D0.9268490.8869400.1439400.1893960~10001000D0.9294180.8887230.1107260.146690~10001000D~E0.9251180.8927940.09856310.1243080~10001000E0.9208180.8968640.08640010.1019470~10001000F0.9294180.8887230.05536340.07333480~10001000表3.11-c小风和静风时扩散参数的系数稳定度0102U100.5m/s1.5m/sU100.5m/sU100.5m/s1.5m/sU100.5m/sA0.930.761.57B0.760.560.47C0.550.350.21D0.470.270.12E0.440.240.07F0.440.240.05•取样时间的订正取样时间指连续取得观测数据并用以计算平均值的时间长度。对污染物浓度计算，即指大气采样的时间。假设大气湍流场是均匀、平稳的，这是Gauss公式的前提，所以在Guass公式中没有指明平均风速和扩散参数的采样时间及观测高度。但在应用中，为了使各种观测资料和计算结果具有可比性，必须考虑这些因素。污染物扩散在铅直方向受地面限制，风向倾角变化范围较小，当取样时间超过几分钟后，铅直风向脉动范围变化不大，所以对铅直向扩散参数可以不做取样时间的订正。水平风向可在360o范围内改变，而表3.11中的扩散参数适用于采样时间为0.5小时的情形，如果采样时间大于或小于这个时段需进行订正。qyy)(1212适用时间范围，hq11000.30.510.2