1、第三章第三章 大气扩散估算大气扩散估算 3.1 大气扩散型 3.2 连续烟流高斯扩散模式3.3 模式计算参数 一、稳定度分类 二、大气扩散参数 三、风速 的计算 四、源强 的计算 五、烟气有效高度H的计算 1. 烟气抬升过程 2.实用抬升公式2222222)(exp2)(exp)2exp(2),(zzyzyHzHzyuQHzyxq;uQzy注意:计算公式中的风速,理论上指烟气有效高度处的风速,由于不方便给定,实用中取烟囱口处的风速。气象测风高度常规为10米:如何确定烟囱口处风速? u10uhsu2222222)(exp2)(exp)2exp(2),(zzyzyHzHzyuQHzyxq; 小球测
2、风小球测风 实际测量 2004-10-12风廓线拟合-平坦开阔地形 :中性层结 对数律非中性,偏离对数律, 符合幂指数律为应用方便,为应用方便, 统一用幂指数律统一用幂指数律 m为经验指数.0*lnzzuupzzzuu01mshuuhs1010风速等于零的高度为粗糙度GB/T 1320191 制定地方大气污染物排放标准的技术方法制定地方大气污染物排放标准的技术方法 下垫面特性:Z0 大,m大;稳定度状况: 越稳定,m越大 和 P140数值不同 例:张家口发电厂位于张家口和宣化之间。设计烟囱高240米。 ?不能用240, 国标规定为200 源强Q是单位时间内从污染源排放污染物的质量。 按时间分:
3、 瞬时源: g或mg 爆炸 间歇源: g或mg 工艺流程 连续点源: g/s或mg/s 烟囱、电厂 按几何形状分 点源: 工厂一定高度的烟囱排放 面源: 稠密居民区中的家用炉灶,城市 线源: 繁忙的公路干线,飞机 体源: 爆炸锅炉燃料不同,生成的污染物不一样。 燃煤锅炉 (SO2, 烟尘,NOx) 燃油锅炉(SO2,, NOx) 燃气锅炉(NOx) 1. 现场实测对有组织排放源 为实测浓度(mg/m3) 为废气体积排放率(m3/s)viviiQqQqQ2.物料衡算法 为投入物料总量 为所得产品量总和 为物料或产品流 失量总和2121GGGGGG3.经验估算法 根据生产过程中单位产品的经验排放系
4、数进行计算,求得污染物排放量的计算方法。只要取得准确的单位产品的经验排放系数,就可以使污染物排放量的计算工作大大简化。 SO2排放量估算 煤中可燃硫一般占总量的70%90%, 平均取80% 根据S氧化成SO2的化学关系: S+O2= SO2 即32克S氧化成64克SO2。因此: 燃煤生成SO2的排放量为: QSO2=2*0.8*B*S=1.6BS(kg/h) B为耗煤量kg/h, S为煤的含硫量%,范围 0.3-5%。 燃油生成SO2的排放量为: QSO2=2BS(kg/h) -全部燃烧 燃天然气 ,可忽略 燃煤产生NOx的估算 一般根据经验数据作估计: 1t煤约生成89kg NOx , 1t
5、石油约生成9.112.3kg NOx , 1t天然气约生成6.4kg NOx 。 燃煤烟尘排放量估算 B为耗煤量kg/hA为煤的尘分含量(%) 烟煤 无烟煤 煤矸石 12-16 34 50-60 d为烟气中烟尘占煤炭总尘分含量的百分数, 一般取值15%20%为除尘效率(%))1 (dABQ尘采取环保措施: 除尘 脱硫 脱硝)21 ()31 () 11 (五、烟气有效高度H 烟气-烟气抬升(烟气-烟流)2222222)(exp2)(exp)2exp(2),(zzyzyHzHzyuQHzyxq; 烟气从烟囱排放后,由于动力和热力作用的结果,烟气会继续上升,经过一段距离后逐渐变平,即烟流的有效高度比
6、烟囱实际高度要高,即: 有效高度H为烟囱几何高度hs与烟气抬升高度h之和。 Hs是一定的,关键是h的问题。 关于烟气抬升高度h的计算,已进行了大量的试验研究和理论探讨。得出过供各种条件下应用的许多烟气抬升高度计算公式,并已初步确定了烟气抬升的物理模型和基本理论体系。hhHs1、烟气抬升过程 造成抬升的原因(抬升因子) 烟气出口速度(初始动量)烟气比周围空气温度高(浮力)烟气自烟囱排出后形成烟流并可上升至相当高度,能够到达的高度与烟源自身条件和周围环境的气象条能够到达的高度与烟源自身条件和周围环境的气象条件有关件有关。一般烟气抬升能将烟源的实际排放高度提高到210倍的有效源高高度上,从而可能使地
7、面最大浓度降低3100倍。 烟气抬升几何模型 有风弯曲烟流模型,抬升大致经历4个阶段:喷出阶段、浮升阶段、瓦解阶段和变平阶段,最后达到烟流抬升高度h的终极抬升阶段。 抬升的四个阶段:抬升的四个阶段: 喷出阶段:烟气在自身具有的(由出口速度提供)作用下垂直向上喷射。 内部流动相对比较规则,边缘上烟气和周围空气的湍流交换尚未发展烟流轮廓清晰,内部基本维持原来状态 随着烟气上升,烟流内外空气开始逐渐发生湍流混合(由速度切边造成),烟流体扩大并获得水平动量(由源高处风速提供)烟流逐渐向水平下风方弯曲。 随着水平动量增大,因初始动量而具有的上升速度减小,主导地位消失,动力抬升转而由浮力作用取代为主导因子
8、。特点:内部流动均匀,外形清晰整齐 (一般维持至烟囱口径10倍左右水平距离的范围) 浮升阶段:由于烟气与环境造成的浮力加速度的作用,使烟气上升速度超过动力上升速度并使烟流继续上升,进入浮力抬升阶段。 烟流体积增大,烟流内外温差和浮力继续使烟流抬升,随着烟流内外的速度切变作用的加入,使更多空气参与混合,烟流边缘的卷夹过程加剧,产生边缘湍流活动带。起主导作用,由速度切变造成卷挟是导致烟气与周围空气混合的主要因素(自生湍流)。 环境湍流的作用还较弱。是热烟流抬升的主要阶段。 持续的混合是烟流体内外的温差不断减小,上升速度减缓烟流开始趋向变平而转入抬升的下一阶段。特点:烟流体增大。特点:烟流体增大。
9、瓦解阶段:瓦解阶段: 浮升阶段后期烟流上升速度逐渐减缓,由速度切变造成的自生湍流减弱。 随着烟流体的不断增大,环境湍流的作用明显增强并逐渐达到占主导。 当烟流逐渐增大到环境湍流含能湍涡尺度时,环境湍流的作用急剧增大,环境湍涡大量卷入烟流体,使其自身结构在短时间内瓦解,烟气原先的动力抬升和热力抬升的性质消失,烟流的抬升基本停止。 这个阶段通常较短。 变平阶段:变平阶段: 以大气中大湍涡为代表的环境湍流起主导作用,使烟流体继续胀大,环境湍流起主要作用,抬升完全停止烟流随风飘动,变平。 此时,烟流达终极抬升高度,并以此计算实际的烟流抬升高度h. 进行这样的阶段划分,可以比较清地知道各阶段的基本影响因
10、子,一方面有利于在理论处理中突出主要因子,另一方面可直接分析影响烟流抬升的各种因子及其作用。 抬升过程复杂多变,不可能完全按上述模型划出典型阶段。大量观测试验表明,影响热烟流抬升的基本因子可归纳成以下3类: 1)排放源及排放烟气的性质。烟气的初始动量和浮力是决定其上升高度的基本因素 。初始动量决定于烟气出口速度w0和源出口半径R0;浮力决定于烟气密度和周围环境的密度差。若不考虑烟气与空气成分的密度差异,压力相同情况下,密度差可用温度差表示,则浮力加速度关系为:该浮力项为烟流抬升的主要贡献项,涉及烟气的热释放量(率),是烟流抬升高度计算公式中的一个主要参量。aasssaTTTgg 2)环境大气的
11、性质。 烟气与周围空气混合的速率(平均风速和环境湍流强度 )对烟流抬升的影响也十分重要。 混合愈快相当于把烟气的初始动量和热量很快分散给周围空气,烟流上升速度很快减少,抬升高度则低,反之,烟流抬升就高。 与混合速率有关的因子主要有,平均风速 ,和环境湍流强度i。 平均风速越大,湍流强度越强,则混合越剧烈,抬升高度越低。烟气所在气层的温度层结表征大气稳定度状况,是影响烟气抬升的又一个重要因子。不稳定,烟流浮力抬升增强,反之,则会抑制烟流的抬升。 u3)下垫面性质。 地形,地面粗糙度 粗糙地面上空湍流活跃,不利于烟流抬升。 离地越高,地面粗糙度的影响减弱,有利于烟流抬升。 复杂地形,形成局部热力状
12、况的特殊分布,影响烟流抬升。 影响烟气抬升的因子有影响烟气抬升的因子有: 内因:初始动量、浮力 外因:大气层结、湍流状况 当烟气出口速度小,浮力小,而周围风速大 下泄,下洗(烟囱口变黑) “Top-Hat”模型模型在不同高度烟流横截面上烟流内部的物理属性,如速度、温度等分布均匀,只是到了烟流的边缘上才突变成与周围大气相应的属性,也就是说,这些属性在烟流内部有一定量值,在烟流外部则有不同的量值,在烟流半宽R处其量值不连续,分布成凸形式,称为“Top-Hat”模型。利用质量守恒、动量守恒、浮力守恒推导烟流抬升方程。利用质量守恒、动量守恒、浮力守恒推导烟流抬升方程。)()()()(2浮力守恒动量守恒
13、质量守恒SVdzdFwFdzVwdRvdzdVzze 为垂直方向的动量通量动量守恒静力稳定度参数浮力守恒,定义卷夹速度质量守恒VWWFdzVWdzgSSVdzdFwVRWRVdzdVzaazee,)(,.22)()()()(2浮力守恒动量守恒质量守恒SVdzdFwFdzVwdRvdzdVzze(A)(B)(C)基于Taylor假设:烟流里的湍流是由烟流和周围流体间的速度切变产生的,因此,所有湍流速度,包括夹卷速度,局部地以及平均烟流垂直速度 表征,也就是,烟流周围的空气被夹卷进入烟流的速率是与烟流和周围流体之间的速度切变成比例的,对垂直烟流 对于弯曲烟流wwvewve弯曲烟流弯曲烟流垂直烟流垂
14、直烟流SVdzdFwFdzVwdRvdzdVzze)(2质量守恒质量守恒动量守恒动量守恒浮力守恒浮力守恒eeRvdtdRorvdxdRu2R222质量守恒质量守恒zzFdtVwdoruFdxVwd)()(动量守恒动量守恒)( VwSdtdForSVdzdFzz浮力守恒浮力守恒有源高处平均风速有源高处平均风速u烟气上升速度烟气上升速度w排放出口半径排放出口半径R0烟流半宽烟流半宽R卷夹速度卷夹速度ve体积通量体积通量V22RwVRuV垂直烟流垂直烟流弯曲烟流弯曲烟流eeRvdtdRorvdxdRu2R222物理意义:单位时间物理意义:单位时间流入的空气体积随离流入的空气体积随离源距离的变化(左项
15、)源距离的变化(左项)与由于卷夹作用进入与由于卷夹作用进入烟流的空气量(右项)烟流的空气量(右项)有关,即维持质量守有关,即维持质量守恒关系恒关系zzFdtVwdoruFdxVwd)()(物理意义:烟气的垂物理意义:烟气的垂直动量增量(左项)直动量增量(左项)与浮力通量(与浮力通量(Fz)和)和风速相关,即维持动风速相关,即维持动量守恒关系。量守恒关系。)( VwSdtdForSVdzdFzz物理意义:烟气的浮力通量随物理意义:烟气的浮力通量随高度的变化(左项)与烟气体高度的变化(左项)与烟气体积通量和大气稳定度(右端相积通量和大气稳定度(右端相关)。关)。不稳定中性稳定0zgSaa大气稳定度
16、参量物理意义是单位质量空气在垂直方向移动单位距离,因浮力作功而引起的加速度变化,即该变量 和改变的成数在体积通量V一定的情况下,对应三种稳定度状况,分别有烟气活跃,不稳定发展浮力通量为常数烟气维持中性平衡烟气受抑制,0dzdFzzazaa1)s/m(3422szTTguRguRF 和T分别表示烟气排放并取得水平速度 后,与周围空气间的密度差和温度差。 物理意义为:单位时间单位质量烟气所受到的浮力贡献(增加动量),或者可以较方便地用排放口烟流参数表示, R0和w0分别为排放口半径和出口速度的初始值 Fz为浮力通量或浮力参数,定义为单位时间、单位质量烟气所受的浮力。)s/m(3422szTTguR
17、guRFusaszTTTgwRF020 考察t时间内排放质量为ms的一团烟气,作用在这团烟气上的力为 f=(ma-ms)g 式中烟气质量与周围空气质量分别为 净举力为 除以 可得浮力通量参数sstwRm2aatwRm2)(2satgwRfatasazgwRF2aassTTsaszTTTgwRF020 在烟流抬升高度的使用计算中常采用烟源的热释放率QH这一特征量,如 单位时间、单位质量烟气升高T度释放的热量aspsHTTcwRQ020cal/saapHzTcQgFTa=273 Ka=1.29 kg/m3g=9.8 m/scp=240 cal/(kgKh)HzQF5107 . 3 推导 通过AA截
18、面单位时间流入的质量 通过BB截面单位时间流出的质量 净增量 净增量由于湍流交换作用而从烟流边缘卷夹进入烟流的空气造成,卷夹速度为ve,则单位时间卷入的空气质量为uR2dxuRdxduR)(22dxuRdxd)(2dxRve2dxRvdxuRdxde2)(2eRvdxdRu22质量守恒质量守恒 通过AA截面单位时间垂直动量增量 式中 为烟流元的平均上升速度,这个动量的增加量是在dx距离(或dx/u时间内)浮力加速度作用的结果,所以有 右项表示单位时间内浮力引起的动量的增加量 消去 ,即得定义为弯曲烟流的浮力通量参数的表达式,可有 根据浮力通量参数的表达式,可有dxwuRdxd)(2wdxguR
19、dxwuRdxd)()(22uFdxwRduz)(2uFdxVwdz)(动量守恒动量守恒uRSguRdzddzdFz22)(22RwVRuV+zgSaaSVdzdFz烟气抬升理论 垂直抬升(无风大气) 定义:平均抬升速度: “体积”通量: 浮力通量: 2RwdsdswdswassszVgdswgF环境位温烟气位温22RwRwwdsV实验得:卷夹系数=即动量抬升时间很短秒设定义转换时间:得:由得:由(中性大气)浮力动量初始条件:转转转4t3mR,300KT,400KT, s/10mwRwTTTgRwFFttFFtFFVwRwdt)C(FF)(0S)(RwTTTgF)(RwVwF0as0200as
20、s20200z0m0z0m0z0m22t0wvwv0zz200ass0z202000m0B0.16(无浮力抬升)射流0.1250.155 (浮升烟云)为垂直方向的动量通量动量守恒浮力守恒VWWFdzVWdzgSSVdzdFzaaz,)(,(B)(C)t转换时间为动力抬升占优势转换为热力转换时间为动力抬升占优势转换为热力抬升占优势的时间抬升占优势的时间弯曲烟云抬升弯曲烟云抬升 定义:水平体积通量:定义:水平体积通量: 垂直动量的水平通量:垂直动量的水平通量: 因因 支配方程:支配方程: uRudsVB2wVwuRuwdsB2dxdzudtdzwt uxdtdzw,zBBzeeBFdtVWdwSV
21、dtdFdZdRWvRvudtdV,2卷夹速度(D)(E)(F) 实验得: 卷夹系数=0.4+1.2(W/u)-1 动量抬升0.6 (浮升)卷夹系数=0.16 动量0.1250.155 浮力由支配方程,可推得弯曲烟云抬升路径(中性大气) 中性层结条件下,S=0, Fz=F0Z=常数 tFFVWzm00积分公式(F)有:zBBzeBFdtVWdwSVdtdFRvudtdV2(D)(E)(F) 2,RuwwVdtdzw因tFFdtdzzuzm002)(可得zR3122020233tuFtuFzzm积分上式,最后可得wve抬升路径(中性大气): 只有动量作用 只有浮力作用实际问题,主要为浮力, 称为
22、中性大气弯曲烟气抬升Briggs2/3路径方程 中性条件下,动量抬升与抬升时间1/3成正比, 浮力抬升与抬升时间2/3成正比6 . 0,6 . 123233233321310321310232310231023120202xuFxuFztuFztuFztuFtuFzzzzmzmt ux 稳定时 S0, 为Brunt-Vaisala频率 利用无风垂直抬升方程中的动量守恒和浮力守恒方程,得到简谐振动方程 利用初始条件:t 0时, 初始动量 初始浮力 简谐振动方程解为: 由上式可知:只有动量作用,达到抬升最高时 只有浮力作用,达到抬升最高时 最大抬升高度:需要相似方法或半经验附加条件才能确定 实验得
23、到: 2122SNzgSaa vwdtvwd222zmFdtvwdFvw0000v200sincos抬升最高时wRwtFtFvwzmt2t8341max00.5SFzz稳定有风大气中的弯曲抬升 由弯曲烟云抬升方程,利用 可得出: 对浮升烟云(F0m=0)有:不稳定层结下的抬升 因 即浮力通量及上升速度均持续增大,抬升无止境。实际上不可能,浮力加强与湍流混合作用相互抵消。由于不稳定时湍流强,实验不易准确,观测资料离散大,实用时用中性公式,只是取比中性是大一些的实验系数。大气湍流对抬升的作用 烟云抬升不能到无穷,原因有二: 一是上部遇到稳定层结, 二是大气湍流作用,湍流越强,抬升越弱,呈反平方规律
24、2zuvB3103120max31003126 .26cos1sin3 suFsuFztFtFsuzzzzm实验0,0dtdFSz 环境湍流与烟气抬升 (不成熟) 环境湍流在烟流抬升的后期阶段有重要支配作用。为得到实用的合理的终极抬升高度,需要考虑环境湍流作用,作为一种闭合假设求解支配方程,求得终极抬升高度;或者根据试验观测,定出终极抬升距离xT,代入抬升高度公式以求得终极抬升高度。(1)引入环境湍流影响的闭合假设 分段作用模式分段作用模式, 假定前期抬升由自生湍流支配,后期则完全由环境湍流支配。 突然作用模式突然作用模式,假定环境湍流在完成抬升的主要阶段不起作用,整个抬升过程始终遵循三分之二
25、次律,一旦环境湍流发生作用,抬升就突然截止。 从抬升的前期就考虑环境湍流的作用,自生湍流和环境湍流的联合作用模型联合作用模型。 (2)受环境湍流制约的烟流抬升 溃散模式和触地模式 机械湍流制约时的烟气抬升 对流湍流制约时的烟气抬升 小结中性大气垂直抬升方程 很快过渡到2/3规律中性大气弯曲烟云抬升 Briggs(浮升烟云)稳定无风抬升稳定有风抬升 上述公式统称为Briggs抬升公式32310223tuFzz动量抬升 t1/3浮升 t2/38341max00 .5SFzz310max6 .2SuFzz不稳定层结下的抬升 实用时用中性公式,只是取比中性是大一些的实验系数。2.实用经验抬升公式烟流抬
26、升理论和实验研究的一个主要目的是建立可行的实用计算公式,供大气扩散模拟中确定抬升高度和有效源高公式带有较强经验性,必须予以观测检验理论研究发展与观测的局限性,不同公式的结果差别较大,缺少一致性以及可比性。实际应用时,应注意认真选择,一方面考虑它们的理论基础和应用条件,一方面注意它们实验依据和可靠性,以能满足扩散计算的准确度要求。 Holland公式(1953)U-出口处水平风速(m/s);w0-烟气出口速度(m/s);D-出口直径(m);QH-烟气的热释放率(卡/秒)。uQDwhH50101.45.132131312)23(xuFz理论公式 R0-排放口半径; W0-出口速度; s-烟气密度;
27、 Ts-烟气温度; Ta-大气温度; Cp-烟气的定压比热(kJ/kg.k)。评论:20世纪50-70年代很流行,但结果偏低很多。 我国:规定小烟囱可用,但需大2倍。)(020aspsHTTcwRQ Lucas公式(电厂) 60年代根据大电厂结果 结果对中、小烟囱偏高,大烟囱偏低。目前少用。 Concawe公式(西欧石油组织专家组) 大型烟囱结果偏低,中小烟囱适用。目前用得较多。uQhH41244321175. 0uQhh32131312)23(xuFz Briggs公式 中性.不稳定 Briggs认为,x=10hs时烟气抬升达到最大高度,故 (相应的下风向距离称为终极抬升距离) 结果偏高.3
28、213104.7szhuFh32131312)23(xuFz 稳定: F0z-浮力通量,m4/s2;u-平均风速,m/s;x-下风距离,m。 评论:Briggs公式理论很严密,国际上最流行,但中性偏 高一点,平坦地形、高烟囱好。 3106.2SuFhz有风:无风:834104SFhz不稳定中性稳定静力稳定度参数000szgSaa 我国标准推荐形式 在我国由国家环保局和国家技术监督局联合发布的国家标准制定地方大气污染物排放标准的技术方法GB/T13201-91,给出了抬升高度计算公式。该方案是以三分之二次律为基础,按我国专门的烟流抬升观测试验资料为依据,给出有关系数,结合国情给出计算公式和处理方
29、法 I.当 千焦耳/秒,且烟气温度与环境气温的差值 其中: 2100HQ时K35asTTT1021uhQnhnsnH;35. 0sVHTTQpQp为气压,取邻近气象站的年平均值,hPa;Qv为实际排烟率,米3/秒;KT35千焦耳/秒地表状况(平原地区)n0n1n2QH21000农村或城市远郊区1.4271/32/3城区及近郊区1.3031/32/3农村或城市远郊区0.3323/52/5城区近郊区0.2923/52/5210002100HQn0为烟气热状况及地表状况系数; n1为烟气热释放率指数; n2为烟囱高度指数;u用幂次律换算到烟囱口m值选择米米200,)10200(200,)10(101
30、0smsmshuuhhuu稳定度ABCD E F城市0.100.150.200.250.30乡村0.070.070.100.150.25 为源出口处(烟囱口高度处)平均风速,米/秒 。 uII. 当 千焦耳/秒时 式中 w0为烟气出口速度,米/秒, D为烟囱出口内径,米 。4001700)(121HQHHHhuQQDwHHH)1700(048. 0)01. 05 . 1 (20121001700HQ较小的类取表中计算按HQnH,2III. 当 千焦耳/秒或 时 IV. 对地面年平均风速 的地方 为排放源高度以上环境温度垂直变化率, 取值不小于0.01K/米。 1700HQKT35uQDwhh/
31、 )01. 05 . 1 (20sm/5 . 18341)0098. 0(50. 5dzdTQhahdzdTa 在抬升过程,出现下层大气中性或不稳定,上部有一定厚度的稳定层。抬升受上部稳定层限制:1)如稳定层底相当低,烟流可能完全伸进稳定层2)如稳定层底较高,烟流不能完全穿透,那么烟流可能部份穿透稳定层。 如果烟流能够穿透上部稳定层,即使只是部份穿透,污染物的地面浓度可以大为降低;否则烟流将封闭在稳定层以下,这会对地面造成较高的污染物浓度。 对大气层结垂直不均一情况,确定烟流抬升高度至关重要,通常,不可能以一个公式计算。 Briggs(1975)假定烟流的终极抬升 落在烟囱上方逆温层高度zi的
32、2倍范围内,只有部份烟(成数p)穿透逆温层。 他提出穿透部分(成数)的公式为hhzpi5 . 1烟流被逆温层底反射并向下散布为(1-p)。当式中: 时, 烟流在混合层顶以上, 全穿透; 时, 烟流在混合层内, 全封闭; 则为 部分穿透。 1pizh20p5 . 1izh 10 pZi2Zi 高架逆温层高度在烟囱顶以上zi处 小结 抬升实验技术较难 抬升问题复杂(地形,烟囱下洗,大气湍流) 推荐用Briggs公式 电厂烟源用concawe公式 复杂地形用Holland公式 一般- 用国标推荐公式GB/T13201-91 制定地方大气污染物排放标准的技术方法 课堂练习答:- 练习作业四 (实际烟囱
33、排放浓度预测计算) 某厂位于乡村,有一锅炉烟囱高30m,出口内径为1.0m,出口烟温140 ,出口烟速5m/s。年耗含硫量为3的重油15万吨,年生产300天。该厂所在地年平均风速为3m/s,年平均温度为22 。试求阴天下风向100、500、1000、2000m处的SO2地面轴线浓度。)()()/()m()2()(3.0/)04.05.1(220KTsKTasmwDwDTaTsQuQDwhoohh为烟温为环境温度为出口烟速为烟囱出口内径假定ABCDEFm0.050.100.150.200.250.30y0.41x0.860.24x0.880.14x0.890.11x0.89z0.05x1.10.12x0.910.12x0.910.82x0.48