（气象学专业论文）沿海地区大气质量模式体系性能研究.pdf

摘要 本论文运用模式性能评价方法，对新一代平原地区大气质量模式和沿海熏烟 模式从参数到模式设置的各个方面进行通盘的考虑，并以定量量度表征，同时再 通过实测资料检验和工程验算等多方面的工作，针对沿海地区大气扩散条件的特 点，提出适合于我国沿海地区的大气质量模式体系的建议。 对国外有代表性的新一代平原地区大气质量模式系统( l p d m 、p p s p 、a l l s ) 及 我国目前广泛采用的e i a 模式对比研究表明：h p 蹦模式在扩散模式、扩散参数、 边界层参数及其它参数的选取上理论更为合理，物理框架更为完善。利用福鼎气 象台一整年( 每天2 4 次) 常规气象观测资料进行模式验证表明，h p d m 模式计算结 果更合理。工程验算时选取3 4 种源条件组合和2 0 0 种气象条件组合，得到了不同 组合条件下的地面浓度分布，h p d m 模式不论对低矮源还是对高架源都有较好的模 拟性能，验算结果的比较与合理性分析是一致的。本论文还指出h p d 模式尚需要 进一步完善之处，即对混合层高度的计算仍有不合理性的地方，在处理烟气穿透 时具有相当程度的保守性。并提出一整套改进建议。 对五个沿海熏烟模式g l u 艘模式、v a nd o p 模式、m i s r a 模式、d e a r d o r f f 修 正模式( m i s r a 修正模式i ) 、s t u n d e r 修正模式( m i s r a 修正模式i i ) 对比研究表 明：m i s r a 模式在理论推导上合理，采用统一的浓度计算公式模拟整个扩散过程， 保证了整个浓度场的连续性，方法较为简便。其中扩散参数采用参数化方法，使 得该模式又具有更大的适用性。利用国电环境保护研究院在杭州湾北仑港电厂进 行的双元示踪试验的结果进行检验，并采用w i l l m o t c ( 1 9 8 2 ) 推荐的度量空气质量模 式性能的方法，对各模式进行评价和比较表明，m i s r a 模式的实测资料检验及模式 评价中的各项指标都明显优于其它模式，这与牧集到的国外已有的检验结果是一 致的。 通过对沿海熏烟模式和新一代平原地区大气质量模式的性能评价，认为m i s r a 模式和h p d m 模式性能较优，能够较好地模拟沿海地区大气污染物输送扩散的情形 和大气污染物的分布，最后，提出适合于我国沿海地区的大气质量模式体系的初 步意见。 关键词；铪海地区空气质量模式检验性能评价对比研究 a b s t r a c t i nt h i sp a p e r ，t h en e wg e n e r a t i o np l a i na i rq u a l i t ym o d e la n dt h ec o a s t a l d i s p e r s i o nm o d e lw e r ec o m p r e h e n s i v e l yc o n s i d e r e di na s p e c t so fm o d e l p a r a m e t e r sa n dm o d e ls e t t i n g su s i n gm o d e lq u a l i t ye v a l u a t i o nm e t h o d ，a n d t h e nt h e yw e r ei n d i c a t e dw i t ht h eq u a n t i f i e dm e a s u r e a tt h em e a nt i m e ， e n g i n e e r i n gc h e c k i n gc o m p u t a t i o n sw h i c hw e r ec a l c u l a t e dw i t hc o n v e m i o n a i o b s e r v a t i o na n do b s e r v e dm e t e o r o l o g i c a ld a t aw e r es t a t i s t i c a l l ya n a l y z e d ； a f t e ra l l ，w ep r o p o s e dt h ea d v i c eo na i rq u a l i t ym o d e ls y s t e ma p p r o p r i a t et o c h i n a sc o s t a la r e a - c o m p a r i s o n sw e r ed o n ea m o n g f o u rp l a i na i rq u a l i t ym o d e l s ，i n c l u d i n g h p d m ，p p s p ，a m sa n de 队i ts h o w st h a th p d mm o d e li sm u c hb e t t e r0 1 1 d i s p e r s i o ns c h e m e a n dp a r a m e t e r s ，b o u n d a r yp a r a m e t e ra n do t h e rp a r a m e t e r s w ea l s ou s eo b s e r v e dd a t at ov e e r yt h ef o u rm o d e l s ，r e s u l t ss h o wt h a t h p d mi sm o r em o d e r a t e 砧ec o n c e n t r a t i o nd i s p u t a t i o n sw e r eg o tu s i n g e n g i n e e r i n gc h e c k i n gc o m p u t a t i o n si n3 4d i f f e r e n ts o u r c ec o m p o u n df o r m s a n d2 0 0d i f f e r e n tm e t e o r o l o g yc o n d i t i o nc o m p o u n df o r m s t h eh p d m m o d e ls h o w sab e t t e rp e r f o r m a n c eo ns i m u l a t i o nd i s a b i l i t y t h ev e i l f i c a t i o n r e s u l ti sc o n s i s t e n tw i t hm o d e r a t ea n a l y s i s m e a n w h i l e h p d mm o d e la l s o n e e d ss o m ef u r t h e ri m p r o v e m e n t s s u c ha sam u c hm o r em o d e r a t e c o m p u t a t i o ns c h e m eo fm i x e dl a y e rh e i g h ta n dam o r ea c t i v es c h e m eo f s m o k ep e n e t r a t i o n w ea l s og i v es o m es c h e m e sf o rm o d e li m p r o v e m e n t s t h ec o m p a r i s o no ft h ef i v ec o a s t a ld i s p e r s i o nm o d e l s ( g l u m p , v a n d o p , m i s r a , d e a r d o r f fa n ds t u n d e r ) s h o w st h a tm i s r am o d a li sm o r e r e a s o n a b l eo nt h e o r e t i c a ld e r i v a t i o n , w h i c hs i m u l a t e st h ew h o l ed i s p e r s i o n p r o c e s s ，i t l lt h eu n i f o r mf o r m u l a r yt og u a r a n t yt h ec o n t i n u i t yo ft h ew h o l e c o n c e n t r a t i o nf i e l d ，t h em o d e lb e o o m em o r ea p p l i c a b l ea f t e rp a r a m e t e r i z e i t sd i s p e r s i o np a r a m e t e r t h e nw er l x a k eu s eo ft h er e s u l t so ft h ef u m i g a t i o n e x p e r i m e n t , w h i c hw a sd o n eb ys t a t ep o w e re n v i r o n m e n t a lp r o t e c t i o n r e s e a r c hi n s t i t u t ei nb e i l u np o w e rp l a n tn e a r b yb a yo fh a n g z h o ut o 2 e v a l u a t et h em o d e l s w ec h o s ea q u a n t i f i e de v a l u a t i o no fa i rq u a l i t ym o d e l c a p a b i l i t yc o m m e n d e db yw i l l m o t t ( 1 9 8 2 ) i ts h o w st h a tt h em i s r am o d e lh a s ab e t t e rp e r f o r m a n c eo no b s e r v e dd a t av e r i f i c a t i o na n dm o d e le v a l u a t i o n t h a no t h e rm o d e l s ，w h i c hi sc o n s i s t e n tw i t hf o r m e rr e s e a r c h e s a c c o r d i n gt ot h ec a p a b i l i t ye v a l u a t i o no ft h en e wg e n e r a t i o np l a i na i r q u a l i t ym o d e l sa n dc o a s t a ld i s p e r s i o nm o d e l s ，m i s r am o d e la n dh p d m m o d e lh a v eb e t t e rp e r f o r m a n c e t h e yc a nd oab e t t e rs i m u l a t i o no ft h e c o a s t a la i rp o l l u t i o nd i s p e r s i o n ，t r a n s p o r t a t i o na n dd i s t r i b u t i o n l a s t l y , w e p r e s e n t e dt h ep r i m a r ya d v i c eo fa i rq u a l i t ym o d e ls y s t e ma p p r o p r i a t et ot h e c o s t a la l e ai nc h i 帆 k e yw o r d s ：c o a s t a la r e a ，a i rq u a l i t ym o d a l ，v e r i f i c a t i o n ，q u a l i t y e v a l u a t i o n ，c o m p a r a t i v er e s e a r c h 3 第一章绪论 1 1 研究的背景目的和意义 随着我国改革开放的不断深入，沿海地区不仅人口越来越密集，而 且工业发展的步伐大大加快，新建、扩建工业企业也越来越多，由此带 来的大气环境的污染问题也越来越引起人们的重视。仅从煤电发展预测， 到2 0 2 0 年，中国沿海地区火电装机容量将达到3 0 一4 0 万娜，年耗煤量 将达到l o 亿吨以上，二氧化硫年排放量将达到2 0 0 0 万吨，扣除脱硫减 排的二氧化硫排放量，年排放量也在2 0 0 万吨以上，如此大量的污染物 排入大气，所带来的污染不容忽视。由于沿海地区处于水、陆交界，水、 陆热力和动力特性的巨大差异产生了一系列特殊的污染气象过程，如海 陆风、热力内边界层，熏烟等，这些与一般的平原地区有明显区别，而 出现上述特殊气象过程，也需要具备一定的条件，在其他情况下，一般 平原地区的大气扩散模式也适用于陆地地形较为平缓的沿海地区。 目前，国内正在考虑对环境影响评价技术导则大气环境 ( h u t 2 2 - 1 9 9 3 ) 进行修订，其中环境空气预测模式计划引进国外新一 代的模式系统。在众多的第二代模式中如何选择更为科学合理且在我国 有应用价值和应用条件的模式系统正是本文探讨的问题。 1 2 国内外研究进展 1 2 1 局地大气模式系统研究进展 我国目前广泛采用的局地空气质量模式和扩散参数体系是国标 g b 3 8 9 0 1 9 9 1 所推荐的高斯烟流模式和扩散参数体系，也是环境影响 评价技术导则大气环境( h u t 2 2 - 9 3 ) 所采用的模式系统。该模式 实际上是以美国e p a 的c r s t e r 模式为蓝本编制而成的，国内外大量的研 4 究和实际工作的经验表明，该模式存在着一些缺陷：( 1 ) c r s t e r 模式对 于热浮力高架点源不完全合适，在某些特殊气象条件下，该模式所预测 的地面浓度偏高，另一些情况又明显偏低。( 2 ) 该模式所提出的p 吒一t 稳定度分类和扩散参数体系，采用的是跃变式的稳定度分级方法，即把 大气环境简单地分为六类不同的稳定度( 也有细分为l o 类的) 。显然不 能细致地反映大气环境的状况。 自上世纪九十年代以来，大气边界层湍流与扩散的研究取得了较快 进展，国外一些发达国家陆续发展和建立了新一代的空气质量模式系统， 如加拿大的a 挞系统，丹麦的o p i l 系统，英国的a d 姻系统以及美国的 h p d m 系统等阻”。相对于第一代空气质量模式，新一代( 也可称为第 二代) 空气质量模式的主要特点是摒弃了传统的p - g _ t 稳定度分类以及 扩散参数体系，在对流条件下，多数采用了非高斯的扩散模式，对e i a 模式进行了全面的发展和升级，摒弃了跃变式的p a s q u i l l 曲线或幂指数 形式的扩散参数体系，采用连续性普适函数或无量纲表达式的形式棘： 烟气达到混合层顶后，考虑到烟气仍保持一定的能量，具有部分穿透的 能力；不稳定条件下( 除近中性外) 基本上摒弃了高斯模式体系。 在沿海地区，熏：疆扩散是造成沿海地区较大污染的主要形式之一， 因此，沿海熏烟扩散模式是沿海地区空气质量模式体系中的一个重要组 成部分，预测向岸流条件下的熏烟浓度的模式一直受到人们的重视。美 国和加拿大曾在大湖区组织过多次大规模的现场试验，并陆续发展了一 批处理熏烟扩散的大气质量模式，主要有g l u m p 模式i s l 咖、v a nd o p 模式1 ”、m i s r a 模式1 味1 1 1 |o e a r d o r f f 修正模式 7 1下沉涡旋修正 模式n 2 。1 3 1 1 1 等。 1 2 2 模式性能研究进展 当出现一批同一类型的大气质量模式时，如何甄别其优劣，上世纪 八十年代，w i l l m o t t ( 1 9 8 2 ) “1 率先提出了模式性能评价方法，在上世 纪九十年代，国内外也出现了些按此方法对同类型模式进行比较研究 的先例“，主要针对模式中的部分参数进行性能评价，系统性较差。本 论文旨在运用模式性能评价方法对各模式系统从参数到模式设置的各个 方面进行通盘的考虑，并以定量量度表征，同时再通过实测资料检验和 工程验算等多方面的工作，提出适合于我国沿海地区的大气质量模式体 系的建议。 l - 3 主要研究内容和思路 鉴于沿海地区大气扩散条件的一般性和特殊性，本论文首先对平原 地区扩散模式的性能进行评价，推荐合适的扩散模式作为我国沿海地区 大气质量模式体系的一部分；然后在针对沿海地区的特殊污染气象过程， 对多种沿海熏烟模式进行性链评价，并推荐其中合适的熏烟模式。 本论文主要研究内容包括： 1 通过谪研和牧资，熟悉各模式及其参数，在计算机上建模。 2 对各模式参数进行对比研究，分析其优缺点。 3 运用模式性能评价方法给出模式性能定量量度。 4 运用常规气象资料和实际观测资料，开展工程验算，并对验算结 果进行统计分析。 5 在上述工作的基础上，提出模式选择的建议。 6 提出适合于我国沿海地区的大气质量模式体系的初步意见。 本论文采用比对、统计、计算机建模、模式验算等方法，对模式的 性能进行全面的分析。技术路线为：调研与收资一模式参数进行比对 运用模式性箍评价方法进行统计分析- - 7 程验算一提出模式选择的建议 一提出适合于我国沿海地区的大气质量模式体系的初步意见。 6 第二章平原地区大气质量模式性能研究 2 1 基本考虑 自上世纪九十年代以来，大气边界层湍流与扩散的研究取得了较快 进展，一些发达国家己逐步发展和建立起新一代的空气质量模式系统， 如加拿大的a m s 系统，丹麦的o m l 系统，英国的a d m s 系统以及美国的 p p s p 系统和h p d m 系统。而我国环境空气影响评价广泛采用的局地空气 质量模式和扩散参数体系是环境影响评价技术导则大气环境 ( h j t 2 2 - 9 3 ) 所推荐的高斯烟流模式和扩散参数体系，该模式实际上 是以美国e p a 的c r s t e r 模式为蓝本编制而成的。 我国目前也正在考虑对现有的空气质量模式系统进行改进和发展， 本论文试图选取国外有代表性的新一代空气质量模式系统( h p d m 、p p s p 、 a m s ) 和我国目前广泛采用的模式( 称为e i a 模式) 进行对比分析，指 出所推荐的模式系统尚需要进一步完善之处，并提出改进建议。 模式的计算和比较主要从以下几个方面进行： 1 扩散模式所需要的边界层参数求取方法的比较； 2 不同气象条件下地面浓度计算方法的比较： 3 模式验证计算。 2 2 边界层参数的求取 以下从稳定度分类、表面热通量q l 、混合层高度z 。等边界层参数求 取方法对h p d m 、p p s p 、a m s 和e i a 四个模式系统进行比较。 7 2 2 1 稳定度分类 上述四个模式系统所采用的稳定度分类方法各不相同( 表2 1 ) ， 其中h p d m 、蠲s 模式完全摒弃了p a s q u i l l 稳定度分类法，将大气边界 层分为稳定、近中性和不稳定三大类，这是稳定度分类方法上的重大突 破。利用江苏省一些主要气象台站某年份1 月和4 月的常规观测资料定 量分析不同模式求取稳定度分类的差异( 表2 2 ) ，其结果表明，h p d m 、 p p s p 和e i a 模式所计算的稳定度分类结果大体上一致，且这种分布趋势 基本上与实际相符，而删s 模式所计算的不稳定类明显偏高。 表2 一l 稳定度分类方法 h p d mp p s p删se i a 白天报据u w 分 类： 根据平均风速与对流速度尺度 甜 之比u w 和m o n i n - o b u k h o v 长 a 3 5 w l 度l 划分三类： 稳定；o 4 w 稳定：q a 5类 一 1 3 9 不稳定时风速超过 w 饥 6 w 归为中性 l o 是根据g o l d ( 1 9 7 2 ) 方法确 d 一1 3 , 9 定的中性时的l w 夜同按t u r n e r 稳定 度分类法分类 8 表2 2 四种模式计算的不同稳定度所占的百分比范围( 单位：) 月份稳定度 呻d m。a 粥p p s pe i a 不穆定5 9约5 0约2 5 5 月 中性4 5 5 5 5 、2 55 0 6 55 0 7 5 稳定3 5 4 - z 5 4 51 0 2 02 3 4 5 不稳定1 2 1 8约5 0约3 5约1 0 四 月 中性4 0 5 05 2 04 0 5 53 5 5 5 稳定3 5 4 03 0 4 51 0 、2 32 5 4 5 结合四个模式所采取的稳定度分类方法和定量计算结果分析，a s 模式分类方法由于仅仅考虑了热力因素的作用，导致不稳定比例明显偏 大，中性比例明显偏小；e i a 模式采用的是卜州稳定度分类法，粗 略地考虑了热力和动力的因素，不稳定比例有些偏小：p p s p 模式夜间情 况下采用t u r n e r 分类法，同样具有一定的局限性，稳定比例有些偏小； 只有h p d m 模式比较充分地考虑了热力因素和动力因素对大气稳定度的 影响，理论上更加合理，分类结果总体上比较理想。 2 2 2 表面热通量q h h p d m 、p p s p 、a m s 模式计算表面热通量q 。的方法有所区别( 表2 3 ， e i a 模式不需要进行q 。的计算) ，测算出的各稳定度下的q h 也有所区别 ( 表2 4 ) 。从计算结果来看，中性条件下，h p d m 模式计算的q h 僮是 几个模式中最小的，一月份! q h i 5 w m 2 ，四月份i q h l 1 0 w m 2 ，这与中 性稳定度的理论意义相一致。 9 表2 3 表面热通量如计算方法 h p d mp p s pa h se i a 白天， 白天： q 铲寸一q e q g q ：净辐射通量 h = o 4 ( 0 6 + o 2s i n 伊) 电；潜热通量 s r s i 行妒) ( 1 q k ) q h = a 凡+ q l q g ：土壤热通量 妒t 太阳高度角 r ，：太阳入射辐射 夜间； 龟一长波损失 q f p c p u e s f = 太阳常数 u 。；摩擦速度 u c ：云反射，吸收因子 9 ；温度尺度 夜间；q h = 一p c p u o 表2 4 四种模式计算的平均值的范围( 单位；i m 2 ) 稳 l p d ha 惦p p s pe i a 定 度 一月四月一月四月一月四月一月四月 a3 0 “7 05 0 7 0 b2 0 8 0 3 0 幸0 5 0 1 0 01 0 0 2 0 04 2 01 5 2 0 c2 52 5 d f 岛 ( 5j 岛f 1 0 0 0 l i 1 2 0 0i l l i s 0 0 4 0 0 6 0 02 0 0 - 3 0 0 0 e5 0 2 0 09 0 - 2 0 0 1 5 - 1 2 01 4 1 1 05 0 一1 0 05 0 一1 0 0 f 约1 4 约1 4 表2 8 四种模式计算的t l 。平均值的范围( 单位：m s ) 稳 弼d ma m sp i s pe l a 定 度 一月四月一月四月一月四月一月四月 a约0 30 3 0 5 b0 1 5 o 3 00 1 5 加4 0 o 3 0 4 约0 40 3 - - 0 40 4 0 5 c0 4 - 060 如o 5 d约0 , 40 4 0 5 约0 4 o 5 加60 2 0 3o 2 0 3 eo 2 0 40 3 , 0 5 0 1 加3o 1 加3约0 2约0 2 f0 0 4 0 0 5 约0 0 5 表2 9 四种模式计算的e 平均值的范围( 单位；k ) 稳 h p d ma m spo s pe 1 a 定 度 一月四月一月四月一月四月月四月 a0 1 - 一0 2一o 1 0 一一0 1 5 b- 0 1 一一0 3 5- 0 1 5 - - 0 3 0一o 1 4 一一0 2 l一0 2 5 一一0 3 5约一o 0 2 一o 0 2 - - 0 0 3 c 约一o ，0 4约- 0 0 0 5 d 0 i o0 7 10 t 0 1 5 1 日i 0 0 2 1o i 5 5 5 5 ) 1 幽：“3 ( 争s ( 旦翳：” “ g 彰阿 觞踽= o 二一n + = 三) 2 l n 蝣助 对流：矗= 3 z 。刀7 5 中性： 不稳定、中性： 风= ( 1 + 号) 2 ” 脯= n o 甜h ”z u 4 w朋 挑”，( 1 + 轰3 u ：烟囱出口处风速 rf 同 醯=警 稳定： 眠： 2 6 ( 云) “3 ( “o ) h p d m u u t 【4 f “4 s 3 ” = 0 )靛：幽：2 6 芒) m 幽；2 6 ( 曼) 啪：h - - - m i n ( h i a h 2 , f ：无量纲浮力通量 h 3 ，h 4 ) l o ：h = m i n ( ah i ，h 3 ，a h 5 ) h ：烟气抬升高度 f = 浮力通量 s ：稳定度参数 表2 一1 2 烟气穿透方案 i p d mp p s pa 惦e i a 当三塾i 5 时，p = 0 ； “幽一一” b r l g g s ( 1 9 7 5 ) 方法 对流： 当 等d s p = 1 5 一垃 m a r t i n s ( 1 9 7 9 ) 不考虑部分穿 a h , 方法稳定及中透，要么完全穿 时一”一簪； p 为穿透系数 性： 透，要么没有穿 a h 。是稳定层结下的 b r i g g s ( 1 9 7 5 ) 透 抬升高度方法 当譬如s 吼间 表2 一1 3 扩散参数 ih p d mp p s pa m se i a i 采用普适函数或无量纲表 采用b r i g g s 扩散 采用普适函数的 采用幂指数形式 ( p g 扩散曲 l 达式( 不稳定) 的形式曲线 形式 线) 表2 1 4 扩散模式 h p d mp p s pa m se i a 1 不稳定低浮力烟羽一 p d f 模式对流：v j c a t r a m 2 不稳定高浮力烟羽一小( 1 9 8 0 ) 模式 风对流模式中性：g a u s s 模 3 ，中性高浮力烟羽一l o f i g a u s s 模式g a u s s 模式 式 模式 稳定：g a u s s 模 4 中性低浮力烟羽和稳定式 g a t l s $ 模式 2 。4 模式验证计算 2 4 1h p d m 模式已有的验证和对比情况 h p d m 模式在建立和完善的过程中，经过了长期的大量的外场观测试 验资料的验证，并与美国e p a 的同类模式进行过分析和对比( 有关情况 见表2 1 5 ) ，表明该模式具有明显的改进，成为目前世界上比较全面和 完善的适用于火电等行业大气环境影响预测的模式。 2 4 2 地面浓度的计算 为了考察四种模式求取地面浓度的差异，利用福鼎气象台( 站号 5 8 7 5 4 ) 一整年( 每天2 4 次) 常规气象观测资料进行模式计算，结果歹 于表2 一1 6 。计算时源条件取为：烟囱高度l l s = 1 5 0 m ，烟气热释放率 q h = 4 6 0 0 0 k j s ，源强q = 2 0 5 6 5 k g h r 。 1 6 表2 一1 5h p d m 的资料验证和模式对比情况 h p d m 3 0 h p d m 4 0 、h p d m 4 2 k m c a i d ( i l l i n o i s ) 电厂 i n d i a n a p o l i s ( i n d i a n a ) 电厂 烟囱高度：1 8 7 m烟囱高度：8 3 m s 0 2 资料：3 0 个测站，5 分钟采样一s f 6 示踪资料：8 9 个样本( 小时平 次，连续疆测一年均浓度) ，离源o 5 1 2 k i n 范围内布 s f 6 示踪资料：2 0 0 个测站，离源置了1 6 0 个测站 0 5 5 0 k m 呈弧形分布，观测分三个阶 段，每阶段3 周时间 b u l lr u n ( t e n n e s s e e ) 电厂 b a l d w i n ( 1 l l i n o s s ) 电厂 烟囱高度：2 4 3 m熘囱高度：1 8 3 m ( 3 个) s f 6 示踪资料：2 0 0 个测站，离源s 0 2 资料：8 7 6 0 样本( 小时平均浓 0 5 5 0 k i n 呈弧形分布，观测分两个阶度) ，离源1 3 - 9 8 k m 内布置了l o 段，每阶段5 周时间 个测站 b a l d w i n ( i l l i n o i s ) 电厂 s u m m i tc o u n t y ( o h ) 工业区 资 烟囱高度：1 8 4 m ( 3 个)烟囱高度； 况 s 0 2 资料：一个固定测站，观测时段 1 9 7 4 - 1 9 7 7 d i c k e r s o n ( m a n y l a n d ) 电厂 烟囱高度： l ，则烟羽高度上的位温梯度采用： 堂：o t o 哆7 伪 ( 三锄) ( 1 5 ) d z