三分仓回转式空预器的详细热力计算方法.pdf
- 文档编号:14653750
- 上传时间:2023-06-25
- 格式:PDF
- 页数:5
- 大小:289.92KB
三分仓回转式空预器的详细热力计算方法.pdf
《三分仓回转式空预器的详细热力计算方法.pdf》由会员分享,可在线阅读,更多相关《三分仓回转式空预器的详细热力计算方法.pdf(5页珍藏版)》请在冰点文库上搜索。
江苏电机工程三分仓回转式空气预热器是一种旋转换热的装置,广泛应用于大型锅炉。
锅炉的三分仓空气预热器与普通的二分仓回转式空气预热器相比,多出一个通道,它将空气流通区域分成一次风和二次风两个通道,满足不同空气出口压力和温度的要求。
三分仓空气预热器温度分布的精度对火电节能及换热元件腐蚀具有重要意义。
目前,二分仓空气预热器的热力计算方法相对比较成熟1,而三分仓空气预热器的热力计算到目前为止并没有完整成熟的方法。
文献2基于二分仓空气预热器热力计算推广到三分仓空气预热器热力计算。
文献3利用数值方法研究了由偏微分方程得到的二维传热方程。
文献4建立复杂的微分方程利用数值方法对内部的温度场进行了模拟。
文献5考虑金属蓄热板竖直切面平均温度沿旋转方向的非线性特性得到出口空气温度的解析表达式。
文献6以大量的试验和统计数据为基础,通过计算直接得到烟气空气传热系数的变化量简化计算。
文献7基于文献5中的方法推广到了四分仓的计算。
上述文献介绍的方法存在一些不足,文献2,3没有考虑到蓄热板平均温度沿旋转方向的非线性特性;文献4,5计算量大,不适用于工程应用;文献5,7虽然避免了上述不足,但是均需将一、二次风合并处理,无法分别计算三分仓的一、二次风区的漏风。
文中在前期研究的基础上5,考虑三个仓之间的漏风差异,建立了三分仓微元体模型,同时考虑金属蓄热板竖直切面平均温度沿旋转方向的非线性及流体温度的非线性特性,对三分仓各个仓进行计算,无需将一、二次风合并处理,从而得到了三分仓空气预热器的详细热力计算方法,且可划分数量不等的计算区域进行计算,以满足不同的工程设计和分析需求。
1数学模型如图1所示,三分仓回转式空气预热器转子转动时,金属受热面依次经过烟气、二次风、一次风,高温烟气对金属受热面加热,金属再加热冷空气。
可以利用微元体建立微元方程,然后在适当简化的条件下得到三分仓空气预热器的计算模型。
1.1微元方程根据空气预热器的工作原理、稳态工况下的流体以及受热面的能量平衡与传热关系可以建立空气预热器转子的微元方程。
如图2所示,微元体的半径与空气预热器转子半径相同,角度为d。
在忽略金属和流体轴向及切向导热的条件下,流体放热或吸热的微元体方程表示为:
cffwf(ad)(tfldl+tf)-tf=f(hddl)(tm-tf)
(1)受热面蓄热或放热的微元体方程表示为:
摘要:
为提高锅炉三分仓空预器的热力计算精度,根据三分仓空气预热器的传热机理及温度分布特点,结合各分仓之间的漏风,建立三分仓微元体数学模型,得到各分仓的微分解析式,采用迭代计算的方法得到各分仓流体平均温度及金属蓄热板的平均温度。
某1000MW三分仓回转式空预器算例表明,该算法具有较高的计算精度,能够考虑各分仓的漏风影响,且可划分数量不等的计算区域进行计算以满足不同的工程需求,可以得到金属蓄热板的平均温度分布,为低温腐蚀和积灰提供温度参考。
关键词:
三分仓;微元模型;传热;对流;漏风中图分类号:
TK222文献标志码:
A文章编号:
1009-0665(2015)03-0010-04陈欢,周克毅,黄军林(东南大学能源与环境学院,江苏南京210096)三分仓回转式空预器的详细热力计算方法收稿日期:
2015-02-12;修回日期:
2015-03-18基金项目:
国家自然科学基金项目(51176031)图1三分仓回转式空气预热器传热示意图二次风一次风烟气图2三分仓空预器微元体示意图ddltmtf+tfldletm+tmdetfeeeeJiangsuElectricalEngineering2015年5月江苏电机工程第34卷第3期10cmmm2n60dl(tmd+tm)-tm=f(hddl)(tf-tm)
(2)式中:
cf为流体比热容,kJ/(kg);pf为流体密度kJ/m;wf为流速,m/s;a为单位弧度流体通流截面积,m2/rad;为通道所占角度,rad;tf为流体温度,;l为受热面高度,m;f为流体与金属的换热系数,kW/(m2);h为单位高度、单位弧度受热面的换热面积,m2/(radm);tm为金属受热面温度,;cm为金属比热容,kJ/(kg);mm为单位高度、单位弧度金属质量,kg/(mrad);n为转子转速,r/min。
由式
(1)和式
(2)可得到如下的微分方程:
tfl=fhcffwfa(tm-tf)(3)tm=60fh2cmmmn(tf-tm)(4)1.2计算区划分及漏风式(3)、(4)的准确性与换热系数、流体的物性参数、金属的物性参数有密切的关系,均与温度相关,即与三分仓空预器内部温度场有关。
理论上当空预器沿高度方向被分成的段数越多、沿圆周方向被分成的区域越多时,由式(3)和式(4)获得的温度场也更精确。
但是,随分成区段增加,其计算量也会增加。
考虑到空气预热器一般都分为热段、冷段,两者材料和板型均不同,为求解的简便,这里可将空预器分成冷热2段。
整个空气预热器被分成6个区分别进行计算,即热段烟气区、热段二次风区、热段一次风区、冷段烟气区、冷段二次风区以及冷段一次风区。
在实际运行中,由于密封问题将不可避免造成各个区之间的漏风,包括一次风区向二次风区漏风、一次风区向烟气区的漏风,以及二次风区向烟气区的漏风。
漏风会对区域的传热和温度分布造成影响,所以在计算各个区域时,要将各漏风量加入对应区域的流量中,并考虑漏风温度的不同,然后确定各区域的流量、流速、平均温度、换热系数等。
按常规处理方法,本文也将各区域的进口作为漏风点1。
1.3各区模型首先考虑漏风影响对进口焓值的改变:
I=I0+1I1+2I2(5)式中:
I0,I分别为考虑漏风前后进口温度对应的焓值;1、2为相邻两区对应的漏风系数,漏进为正,漏出为负;I1,I2为相邻两区对应的平均焓值。
然后由新的焓值求出新的进口温度,再用差分替代式(3)中的微分,并用考虑漏风影响后各区流体进出口的平均温度作为流体传热温度,可得各区传热和蓄热关系式:
tj,n-tj,n-1=fj,nHj,n2cfj,nfj,nWfj,naj,n(tmj,n-tj,n+tj,n-12)(6)tmj,n=60fj,nHj,n2cmMj,nn(tj,n+tj,n-12-tmj,n)(7)式中:
下标n=1,2,表示沿流体方向将空预器分成2段,靠近进口段n为1,靠近出口段n为2;下标j=1,2,3,分别表示烟气区、二次风区、一次风区;fj,n,cfj,n,fj,n,aj,n为第j区第n段考虑漏风后的参数;Hj,n为第j区第n段的换热面积;tj,n-1,tj,n为第j区第n段流体考虑漏风后的进口平均温度和出口平均温度,;tmj,n为第j区第n段的金属受热面竖直切面的平均温度,;Mj,n为第j区第n段金属转子总的质量,kg。
Wfj,n为第j区第n段在(+1+2)过量空气系数下求得的流速,其中为该区域的过量空气系数,m/s。
由于漏风影响导致质量流量变化,流速将会改变,流速的计算为:
wfj,n=Bca1V(t+273)273F(8)式中:
Bca1为计算燃料消耗量,kg/s;V为(+1+2)过量空气系数下的体积,m3/kg;t为考虑漏风后区域的平均温度,;F为通流截面积,m2。
记:
Aj,n=fj,nHj,n2cfj,npfj,nWfj,naj,nBj,n=60fj,nHj,n2cmMj,nnnnnnnnnnnnnnnnn(9)代入式(6)、(7)得各个分区所满足的关系式:
tj,n-tj,n-1=Aj,n(tmj,n-tj,n+tj,n-12)(10)tmj,n=Bj,n(tj,n+tj,n-12-tmj,n)(11)由式(10)可得到:
tmj,n=2+Aj,n2Aj,ntj,n-2-Aj,n2Aj,ntj,n-1(12)tj,n=2Aj,n2+Aj,ntmj,n+2-Aj,n2+Aj,ntj,n-1(13)将式(13)代入式(11),可得到:
tj,n=-2Bj,nAj,n+2(tj,n-tj,n-1)(14)将式(14)积分得到弧度为处的流体出口温度:
tj,n()=tj,n-1+Cj,neRj,n(15)式中:
Rj,n=-2Bj,n/(Aj,n+2);Cj,n为常数。
将式(15)代入式(12)得到金属蓄热板过转子圆心的竖直切面平均温度:
tmj,n()=tj,n-1+2+Aj,n2Aj,nCj,neRj,n(16)eeeeeeeeeeee陈欢等:
三分仓回转式空预器的详细热力计算方法11江苏电机工程由式(15)积分,可得各个分仓流体的平均出口温度为:
t軃j,n=1jj0乙tj,n()d=1jj0乙(tj,n-1+Cj,neRj,n)d(17)式中:
j为第j区在空气预热器中所占的弧度,rad。
若已知流体进口平均温度和出口平均温度,可将已知条件代入式(17)求解积分方程式得到常数,从而由式(16)求得金属蓄热板竖直切面进口平均温度:
tmj,n(0)=tj,n-1+2+Aj,n2Aj,nCj,n(18)和出口平均温度:
tmj,n(j)=tj,n-1+2+Aj,n2Aj,nCj,neRj,n(19)1.4求解的唯一性整个空气预热器共分成6个区,每个区都有3个方程,即式(14)、(15)和(16),共有18个方程。
另外每个相邻区金属蓄热板温度相等,增加6个方程。
这样整个空气预热器合计有24个方程。
每个区都有4个未知数,即t軃j,n,Cj,n,tmj,n(0)和tmj,n(j),共有24个未知数,与24个方程相对应,可解得唯一解。
2计算2.1计算步骤根据文中计算模型,可分别对冷热段各个仓进行计算,由式(17)、(18)和(19)迭代求解出各个区流体的出口平均温度。
每段的计算流程如图3所示,其总体计算步骤为:
(1)已知热段烟气进口平均温度、冷段一次风及二次风进口平均温度,假定热段、冷段烟气出口平均温度,热段、冷段二次风出口平均温度,热段、冷段一次风出口平均温度,考虑漏风影响得到新的进口温度、流速等参数,由式(17)得到常数C1,1,C1,2,最后由式(15)和式(18)得到热段、冷段烟气区金属蓄热板的进出口平均温度tm1,n,tm1,n;
(2)按照烟气区求解的方法,得出二次风区金属蓄热板的进出口平均温度tm2,n,tm2,n;(3)按照烟气区求解的方法,得出一次风区金属蓄热板的进出口平均温度tm3,n,tm3,n;(4)最后将冷热段的金属蓄热板进出口平均温度进行比较,若满足:
tm1,1-tm2,2,tm2,n-tm3,n,tm3,2-tm1,1,tm1,2-tm2,1,tm3,1-tm1,2,则迭代结束,否则重新假定,转
(1)。
2.2参数选取
(1)考虑漏风时,一次风区向烟气区和二次风区的漏风、二次风区向烟气区的漏风取设计值。
在计算漏风影响时,先考虑漏风对进口流体温度、流量的影响,再按照模型计算流体出口温度及金属蓄热板温度。
(2)考虑到烟气侧及空气侧的散热,文献采用与计算标准一致的方法,计算A1,n时,在其分母上乘以保温系数,即A1,n=f1,nH1,n/2cf1,nf1,nWf1,na1,n。
(3)考虑到受热面污染、未完全冲刷等因素的影响,在计算换热系数fj,n时,考虑热利用系数的影响。
(4)在计算Aj,n,Bj,n时,涉及到计算各参数,有些参数的计算在标准中没有涉及。
因为烟气是由各种气体混合而成,所以烟气的部分参数计算如下:
cf=hy(1+L0)T(20)式中:
hy为平均烟焓,kJ/kg;L0为完全燃烧1kg燃料理论上需要的空气质量,kg。
fwf=GyBj2S(21)式中:
Gy为1kg燃料烟气质量,kg/kg;Bj为计算燃料消耗量,kg/s。
空气的部分参数如下:
f=v=Pr/cfv(22)其余参数可参考文献8计算获得。
3计算结果及分析文献8详细介绍了三分仓空气预热器的热力计算,其结果具有较好的参考性,而文献5在文献8的基础上提出了改进,所得的结果更加精确。
为验证文献改进方法的可行性及计算结果的合理性,分别按照文献8、文献5和文中提出的方法,编制了三分仓热力计算程序。
图3每段的计算流程开始计算考虑漏风影响后各个区的进口流体温度假设烟气出口温度二次风出口温度一次风出口温度计算得出各个区金属板进出口温度判断金属板烟气区出口与二次风区进口温度二次风区出口与一次风区进口温度一次风区出口与烟气区进口温度是否相等结束NY12陈欢等:
三分仓回转式空预器的详细热力计算方法(下转第16页)项目数值项目数值转子直径/m16.4二次风进口温度/18热段受热面高度/m1.042一次风通道所占份额0.139冷段受热面高度/m1二次风通道所占份额0.278转速/(rmin-1)1.2烟气通道所占份额0.458热段总受热面积/m293250一次风区向烟气区漏风0.04冷段总受热面积/m262120二次风区向烟气区漏风0.04热段烟气进口温度/374一次风区向二次风区漏风0.06一次风进口温度/20表1某1000MW超超临界三分仓空气预热器原始数据表2算例计算结果参数设计值文献8方法文献5方法文中两段文中四段出口烟气温度124125.6125.6123.6124.0一次风出口温度323307.5317.4319.0319.5二次风出口温度336332.6332.5334.5333.0图4出口空气温度分布050100150文中空气温度;150一二次风出口温度/300310320330340350360370周向位置/()文献5空气温度;文献5平均温度;文中平均温度图5金属蓄热板竖直切面平均温度分布密封区密封区密封区060120180240300360金属蓄热板竖直方向平均温度/周向位置/()100140180220260300热段冷段烟气区二次风区一次风区算例为某1000MW超超临界锅炉2台对称分布的三分仓空气预热器,其计算所需原始数据及3种方法的计算结果分别列于表1和表2。
3.1计算结果分析由表2的计算结果及计算模型,可知文中详细热力计算方法有以下特点。
(1)具有较高计算精度。
文中热力计算方法中排烟温度、一次风温度和二次风温度与设计值的差值分别为0.4、4和1.5,均小于文献5和文献8。
因为文中方法不但考虑了金属蓄热板竖直切面平均温度沿旋转方向的非线性及流体温度的非线性特性,而且分别考虑各个漏风的影响,并且不涉及二分仓回转式空气预热器的计算方法,故文中计算方法更为合理。
(2)可以对空预器进行更加详细的划分,以满足不同的分析及设计要求。
表2中可以看到:
将空预器沿高度分成冷热两段共6个区域或将冷热两段再分两段得到四段共12个区域的结果差值比较小,但分成四段后可得到更多的金属和空气温度,可以根据需求划分区域得到空气预热器不同位置的温度分布。
3.2空气温度变化和金属温度分布根据式(16)可以得到出口空气温度沿旋转方向的变化,将文中国得到的曲线与文献5进行比较。
由图4可见:
文献5一、二次风区的曲线斜率基本一致,即认为空气出口温度在一次风区、二次风区变化快慢一致。
而文中计算方法所获得的曲线在二次风区变化快,一次风区变化较慢。
因为金属蓄热板经过烟气区加热后依次经过二次风区、一次风区冷却,其在二次风区的温度较一次风区要高,二次风区的传热温差大,且二次风区的换热系数大于一次风区,所以二次风区换热更强烈,二次风区空气出口温度变化应该比一次风区要快,故文中方法的计算结果更合理。
在计算烟气、空气出口温度的同时,还能得到金属蓄热板竖直切面平均温度沿旋转方向的温度分布,图5为算例金属蓄热板竖直切面平均温度随旋转方向的变化。
图中是以烟气区金属蓄热板进口为起点,温度曲线上升、水平、下降所对应的范围分别为烟气区、密封区、空气区。
在烟气、空气区内,温度变化接近线性特性。
图中实线代表冷段金属蓄热板竖直切面的平均温度,为判断受热面低温腐蚀和积灰提供了理论依据。
若想获得更准确的温度分布,还可以对低温段进行更详细的划分。
4结束语对三分仓回转式空气预热器微元体进行数学建模推导,得到了一种详细的热力计算方法,该方法根据推导的各分仓微分解析式,采用迭代计算的方法得到各分仓出口流体平均温度及金属蓄热板的平均温度。
与以前介绍的热力计算方法比较,文中的热力计算方法精度较高,考虑了各分仓的漏风影响,同时能得到金属13江苏电机工程蓄热板竖直切面平均温度的分布。
实际应用中,还可以划分更多计算区,得到沿烟气流动方向上不同位置处的金属和流体温度,以满足不同的分析或设计要求。
参考文献:
1北京锅炉厂译.锅炉热力计算标准方法M.北京:
机械工业出版社,1976:
86-88.2周俊虎,杨卫娟,靳彦涛,等.三分仓空气预热器热力计算的研究J.动力工程,2003,23(6):
2810-2813.3冷伟,王渡.一种改进的回转式空气预热器热力计算方法J.动力工程,2005,25(3):
392-395.4WANGHY,ZHAOLL,XUZG,etal.TheStudyonHeatTransferModelofTri-sectionalRotaryAirPreheaterBasedontheSemi-analyticalMethodJ.AppliedThermalEngineering,2008,28:
1882-1888.5郑凯,周克毅,黄军林.一种改进的三分仓空气预热器热力计算方法J.动力工程学报,2012,32
(1):
36-40.6李杨,周元祥,谢天,等.三分仓空气预热器热力计算的简化方法J.动力工程学报.2013,33
(1):
37-417陈昌贤,孙奉仲,李飞,等.四分仓回转式空气预热器热力计算方法J.动力工程学报,2014,44(4):
59-638周强泰,周克毅,冷伟,等.锅炉原理M.2版.北京:
中国电力出版社,2009:
201-204.作者简介:
陈欢(1989),男,江苏无锡人,硕士研究生,从事电站锅炉特性方面研究;周克毅(1955),男,江苏泰州人,博士生导师,从事电站锅炉方面研究;黄军林(1986),男,山东潍坊人,博士研究生,从事电站锅炉方面研究。
ADetailedThermalCalculationMethodForTri-sectionalRotaryAirPreheatersCHENHuan,ZHOUKeyi,HUANGJunlin(SchoolofEnergyandEnvironment,SoutheastUniversity,Nanjing210096,China)Abstract:
Toimprovetheaccuracyofthermalcalculationfortri-sectionalrotaryairheaters,adifferentialelementmodelhasbeensetupbasedontheheattransferprincipleofrotaryairheatersandthecharacteristicsoftemperaturedistribution,consideringtheinfluenceofairleakage.Andthentheanalyticalexpressionsofeachsectorhavebeenobtained.Atlast,thefluidtemperatureandthemeantemperatureonverticalsectionofmetallicheataccumulationplatecanbecalculatedbymeansofiterationbasedontheanalyticalexpressionsofeachsector.When1000MWtri-sectionalregeneratorisselectedascomputationalobject,theresultsindicatethatthenewmethodhashighprecision,andthecalculationareacanbedividedasrequiretomeetdifferentdemandsofproject.Also,thismethodcanobtainthemeantemperaturefieldoftheverticalsectionofmetallicheataccumulationplate,whichcouldprovidetemperaturereferenceforanalyzinglow-temperaturecorrosionandashdeposition.Keywords:
tri-sectionalairpreheater;differentialelementmodel;heattransfer;convection;leakage参考文献:
1刘振亚.特高压交直流电网M.北京:
中国电力出版社,2013:
90-91.2国家发展和改革委员会.DL/T10892008直流换流站与线路合成场强离子流密度测量测量方法S.北京:
中国电力出版社,2008.3国家发展和改革委员会.DL/T10882008800kV特高压直流线路电磁环境参数限值S.北京:
中国电力出版社,2008.作者简介:
李国奇(1985),男,河北唐山人,工程师,从事电网环境监测工作;华伟(1979),男,江苏扬州人,高级工程师,从事电网环境监测工作;傅高健(1980),男,江苏镇江人,工程师,从事电网环境监测工作;祁建民(1980),男,江苏盐城人,工程师,从事电网环境监测工作。
MonitoringandanalysisofthetotalelectricfieldofDClinesinJiangsuprovinceLIGuoqi,HUAWei,FUGaojian,QIJianmin(JiangsuFrontierElectricTechnologiesCo.Ltd.,Nanjing211102,China)Abstract:
ThetotalelectricfieldofDClinesatgroundlevelandsensitivenearbygroundinJiangsuprovincehasbeenmeasured.Accordingtothemeasurementresults,severalconclusionshavebeenobtained:
thetotalelectricfieldatgroundlevelisweakestatthemiddlepositionofDCline;thetotalelectricfieldatsensitivenearbygroundhasatrendthatitincreasesfirstandthenreduceswiththeincreaseofthedistance;theheightofDClinehasanobviousinfluenceontotalelectricfieldofDClinesatgroundlevel;andalargerdistanceoftheindicatesaweaktotalelectricfieldofDClinesatsensitivenearbyground.Keywords:
DCline;totalelectricfieldatgroundlevel;motoring:
(上接第13页)16
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 三分 回转 式空预器 详细 热力 计算方法