地源热泵竖直埋管的有限长线热源模型.pdf

《地源热泵竖直埋管的有限长线热源模型.pdf》由会员分享，可在线阅读，更多相关《地源热泵竖直埋管的有限长线热源模型.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、收稿日期:2002-04-20;修订日期:2002-07-10作者简介:曾和义(1965-)男,四川广汉人,山东建筑工程学院高级工程师,硕士研究生.文章编号:1001-2060(2003)02-0166-04地源热泵竖直埋管的有限长线热源模型曾和义,刁乃仁,方肇洪(山东建筑工程学院 地源热泵研究所,山东 济南 250014)摘 要:对地热换热器竖直埋管的非稳态传热模型进行了分析讨论。采用虚拟热源和格林函数法给出了半无限大介质中有限长线热源产生的非稳态温度场的解析解表达式。与稳态温度场的解进行比较,讨论了温度场达到名义上的“稳态”所需的时间,同时对于达到稳态时的温度场也进行了分析,指出了现行教科

2、书中关于该问题的错误,提出了稳态时两个地热换热器孔壁代表性温度的定义,并对两者进行了比较,进而给出了可供工程应用的简化计算公式。基于以上分析,进一步讨论了全年冷热负荷不平衡对地热换热器长期性能的影响。关键词:地源热泵;地热换热器;导热;传热模型中图分类号:TK17 文献标识码:A1 引 言在地源热泵空调系统的设计中,确定地热换热器的换热模型并进行可靠的计算是非常重要的。在地热换热器诸多布置形式中,竖直埋管布置形式因其较水平埋管布置形式有许多的优点而得到广泛的应用1。竖直埋管布置的地热换热器也就是把管子埋在地表以下的竖直钻孔中并用封井材料填实,流体在管内流动并与土壤进行换热。现在工程上常见的地热

3、换热器的钻孔直径为100200 mm,钻孔深度40200 m。与所埋管子的深度相比较,钻孔直径是很小的。同时,土壤可以看作是一半无限大介质。因此,埋有管子并与土壤进行着热交换的钻孔,就可以被近似地看作是置于半无限大介质中的线热源而进行传热分析。这种管内流体与周围土壤的换热是非稳态的,对于地热换热器运行工况的分析,常采用无限大介质中的线热源模型2。这种模型比较简单,但由于忽略了地面作为一个边界的影响,当时间趋于无穷大时,温度场不会趋于稳定。因此不能用来讨论长时间的问题。在很多情况下,地热换热器全年的冷热负荷是不平衡的。例如在北方建筑物冬季的供热负荷和时间远大于夏季的空调负荷和时间;而在南方情况则

4、相反。如果在一年中冬季从地热换热器中抽取的热量与夏季向地热换热器输入的热量平衡,则地热换热器在数年的长时间运行后,地下的年平均温度没有变化,对地热换热器的性能没有影响3。但是,如果地热换热器的吸热和放热不平衡,多余的热量(或冷量)就会在地下积累,引起地下年平均温度的变化。由于地表作为一个边界的影响,当时间足够长以后,这种不平衡热量的效应将会达到一个基本稳定的状态。这一过程的长短取决于钻孔的深度,通常需要十余年。这种换热器周围土壤年平均温度的变化会影响地热换热器长期的换热性能,这是在设计时所必须考虑的。本文拟对此非稳态的传热过程进行分析,并与达到稳定时的结果进行对比分析。根据线性迭加原理、格林函

5、数方法及虚拟热源方法,可以得到有限长线热源在半无限大介质中温度分布的解析解。这个解比无限长线热源引起的温度场的解更好地反映了地热换热器的长期运行的实际工况,可以作为讨论和计算地热换热器实际的非稳态工况的一个基础,但本文的重点是讨论不平衡热量的效应对地源热泵地热换热器长期性能的影响。2 数学模型如图1所示,在半无限大介质中初始温度均匀为t0。如果在某一时刻开始,垂直于边界表面的、强度为ql(Wm)的有限长匀强线热源开始放热(或吸热)。由对称性可知,这一温度分布在柱坐标系中是二维的。选取介质表面温度t0(也就是初始时刻介质的温度)为过余温度的零点,即设=t-t0。利用虚拟热源法原理,在与线热源关于

6、边界面对称的位置上设一虚拟线热汇,其强度为-ql,长度同样为H。第18卷(总第104期)热能动力工程2003年3月 1995-2005 Tsinghua Tongfang Optical Disc Co.,Ltd.All rights reserved.图1线热源和线热汇及其几何关系由于问题的线性性质,时刻点M(,z,)处的温度就是线热源与线热汇上微元段在此点产生的过余温度的迭加,我们导得:=ql4kH0erfc2+(z-h)22 2+(z-h)2-erfc2+(z+h)22 2+(z+h)2dh(1)k和分别为介质的导热系数和热扩散系数。对式(1)进行无量纲化:令Z=zH,H=hH,R=H,

7、F0=H2,=4k(t-t0)ql可得:=10erfcR2+(Z-H)22F0R2+(Z-H)2-erfcR2+(Z+H)22F0R2+(Z+H)2dH(2)可见,无量纲温度 是无量纲变量Z、R、F0的函数,即=f(Z,R,F0)。如果取Z=0.5,R取定为孔壁处的相对半径Rb=rbH,则孔壁中点的无量纲温度b是Rb和F0的函数。对于不同的钻孔,由式(2)可得孔壁中点无量纲温度随F0(随时间)的变化关系图(如图2)。从图2可以看出,随着时间的推移(即随着F0的增大),无量纲温度 刚开始是急剧增大,然后逐渐趋于平缓。理论上讲,只有当F0 时,温度场才能到达稳定。如果定义:当某处非稳态的温度达到稳

8、态时温度的0.98倍时,称在该处非稳态传热达到了“稳态”,则可得出非稳态传热达到稳态所需的时间,如图2。图中示出了在五条不同的相对半径条件下无量纲温度随傅里叶数变化的关系曲线图。从图中可以看出,温度达到稳态时所用时间与点的坐标有关系,如果距线热源距离越远(即R越大),按以上定义达到稳态所用时间也就越长,这是“半无限大介质”边界条件引起的比较特殊的情况。按图2所给出的曲线关系,可用回归方法得到坐标R与该处达到名义上的“稳态”所需的时间的关系式为:F0=3.8Rb+0.29(0.000 1Rb0.75(3)图2无量纲温度与傅里叶数间关系曲线对于一个典型的钻孔:rb=0.055 m,H=55 m,则

9、Rb=0.001,则孔壁非稳态温度达到名义上的稳态所需的时间为F0=0.293 8。如=310-6m2s,则=9.394年。当F0 时,以上非稳态的解趋于稳态的温度场,其解析表达式为4:s=q14kln(H-z)2+2+(H-z)(H+z)2+2+(H+z)2z2+2z2+z2+22(4)同样引入无量纲量对上式进行化简,则式(4)可无量纲化为:s=ln(1-Z)2+R2+1-Z(1+Z)2+R2+1+Z761第2期曾和义,等:地源热泵竖直埋管的有限长线热源模型 1995-2005 Tsinghua Tongfang Optical Disc Co.,Ltd.All rights reserve

10、d.2Z2+2ZR2+Z2R2R2(5)可知稳态时温度场中无量纲温度 s只是无量纲坐标Z和R函数,即s=f(Z,R)。3孔壁中点稳态温度简化计算式在地源热泵的地热换热器的设计中,最关心的是钻孔壁上的温度,即在=rb处的温度。对于给定的钻孔尺寸,孔壁的无量纲温度将是相对深度Z的函数。给定不同的钻孔直径,可得不同的孔壁稳态过余温度随相对深度的变化关系图。图3是相对半径Rb=rbH分别为0.000 5、0.002 75、0.005时孔壁的无量纲稳态过余温度随相对深度变化的关系曲线。从图中可看出,随着相对深度从零逐渐地增加,无量纲稳态过余温度先是急剧增大,在孔的中段经过一个相对平缓的变化段,在接近孔的

11、底部时又急剧下降。当Z超过1后,无量纲稳态过余温度渐渐趋于零。图3无量纲温度沿钻孔深度变化关系图在工程实际中,常常希望用一个代表性的温度来表征孔壁上的过余温度。一个常用的方案是取z=H2处孔壁的温度作为孔壁的代表温度。通常情况下,在地热换热器中,钻孔半径rb相对于钻孔深度H是一个很小的量,即rbH1。通过对式(5)进行化简,可以得到计算孔壁中点稳态过余温度的近似计算公式:b=rb,H2ql2klnH3rb,(rbH1)(6)值得指出的是,在很多传热手册和传热学教科书59中都讨论了这一问题,并给出了这一问题的解:3b=ql2kln2Hrb,(rbH1)(7)实际上这个解是错误的。式(7)所表示的

12、解适合于边界绝热的半无限大介质中有限长线热源产生的稳态温度场,但都被错误地用作等温边界条件的解。当Rb=rbH=0.001时,这两者之间的偏差可达19.5%;对于较大的相对半径Rb,两者的偏差更大。因此在使用这一公式时应注意纠正教科书中错误。4孔壁积分平均稳态温度用孔壁深度方向中点处的温度作为孔壁的代表温度虽然比较方便,但在定义上显得有些随意性。也可以以沿整个孔壁深度方向温度的积分平均值作为孔壁代表稳态温度,这样在传热计算中更为合理,也就是:b=1HH0(rb,z)dz(8)或 b=4kbql=H0(Rb,Z)dz(9)积分平均温度可以由式(5)通过数值积分的方法得到。从式(6)和式(9)可得

13、知,稳态时无量纲的孔壁中点温度和积分平均温度均只与相对半径大小有关。图4分别绘制出这两个温度值随相对半径变化的关系曲线图,从图4可以看出,积分平均值与管壁中点温度的值相比较,积分平均值要小于管壁中点温度的值。对于工程中常用的钻孔半径(相对半径为0.000 50.005),积分平均值与由简化公式计算出的管壁中点温度的值其相对误差大小在3.4%5%之间。但是,积分平均值计算过于繁琐,用于工程实际显然是不方便的。从图4可以看出,积分平均温度与相对半径的对数值近似于直线关系,由此可以用线性回归的方法确定孔壁积分平均温度与相对半径的近似关系式。这样,可以整理出形式如式(6)便于应用的计算公式:bql2k

14、lnH2.2rb(10)861热能动力工程2003年 1995-2005 Tsinghua Tongfang Optical Disc Co.,Ltd.All rights reserved.图4积分平均温度与管壁中点温度比较通过与积分平均温度的精确值进行比较可知,当0.0005 rb/H 0.005时,式(10)的相对误差均小于0.5%。5影响孔壁稳态过余温度的因素由式(4)可知,影响孔壁稳态过余温度的因素较多,包括线热源全年平均热强度ql、介质导热系数k、埋管深度H、钻孔的半径rb。经过无量纲处理后可知,影响稳态时孔壁无量纲过余温度的因素为相对埋管深度Z和相对钻孔半径Rb。对于钻孔孔径以及

15、埋管所需释放的热负荷(即线热源强度ql)一定的实际工程,为保证换热器的长期运行性能,通常须保证稳态时孔壁过余温度不得超过某一给定的值。对于一个典型的单个钻孔地热换热器,rb=0.055m,H=55 m,则Rb=10-3。由式(10)可得此时无量纲稳态温度为=12.23。对于一个由N个钻孔组成的地热换热器,在某一个孔壁处的在任一时刻的温升可由该组换热器的各钻孔在该处引起的温升迭加得到,即e=(rb,)+N-1i=1(ri,)(11)其中:ri表示第i个钻孔与所考虑的钻孔间的距离。当然上式也适用于稳态的情况。如果地热换热器是一个由六个钻孔组成的23矩形阵列,钻孔间的间隔都为6 m,则按照迭加原理,

16、可以计算得到稳态时最不利的钻孔的孔壁的无量纲温升为=12.23+32.85+22.16=25.10。如果地热换热器每个钻孔每米孔深所分担的全年不平衡的热负荷为3.15105kJm,即全年平均的稳态热负荷为ql=10 Wm,土壤的导热系数为k=2.5W/(m.K),则可得以上单孔换热器的孔壁稳态温升为b=3.9;六个钻孔组成地热换热器的孔壁稳态温升为b=8.0。如果钻孔的数量再增加或间距减小,相邻钻孔的影响还会进一步增加。以上计算结果表明,对于全年冷热负荷不平衡的地热换热器,例如北方寒冷地区的地源热泵系统,必须考虑这种由于冷热负荷的不平衡所带来的长期影响,否则将不能保证系统在若干年的运行以后仍能

17、承担正常的负荷。其次,对于孔径一定的钻孔来说,加大孔深H将不利于从地表的吸热(或散热),因而会增大不平衡热负荷所引起的稳态温差。不过,在工程通常所能调节的孔深范围内,单纯靠减小孔深来改善这种长期工况的特性,其作用是有限的。对于由多孔组成的地热换热器,加大孔间距有利于减轻各孔间的相互影响,但可能受占地面积的限制。因此,对于选择的地源热泵系统,关键还是在于减小冷热负荷的不平衡程度,亦即减小全年平均的稳态热负荷ql。为此可以考虑采用混合系统9,即当地热换热器全年的冷热负荷不平衡时,可以考虑按其中较小的负荷设计地热换热器。全年运行中地热换热器负担平衡的冷热负荷,而不平衡的部分由其它辅助的加热(或冷却)

18、设备来完成。6结 论(1)半无限大介质中有限长线热源产生的温度响应更好地描述了竖直埋管地热换热器在长期运行时的传热过程,可以作为分析和计算竖直埋管地热换热器的传热的基础。(2)对于地热换热器的工程实际应用,可以考虑采用孔壁中点温度或孔壁积分平均温度作为达到稳态时孔壁的代表性温度,但两者的差别可达5%。由简化计算公式(6)计算的孔壁中点温度作孔壁的代表性温度值是通常推荐的方法;而由孔壁积分平均温度拟合得到的计算公式(10)则是更合理、同时也没有增加计算难度的计算公式。(3)在设计地源热泵系统时,冷热负荷的不平衡程度对系统的设计有重要的影响。为保证地源热泵系统在长期运行当中能有足够的制热(或制冷)

19、能力,在系统设计选择时应力求减小冷热负荷的不平衡程度,即减小全年平均的稳态热负荷ql,最好能够使冷热不平衡度ql为零。否则,应增加地热换热器埋管的长度,或考虑采用混合系统。961第2期曾和义,等:地源热泵竖直埋管的有限长线热源模型 1995-2005 Tsinghua Tongfang Optical Disc Co.,Ltd.All rights reserved.收稿日期:2002-03-20;修订日期:2002-06-02作者简介:李运泽(1972-),男,河北邢台人,清华大学博士研究生.文章编号:1001-2060(2003)02-0170-03大型锅炉长期动态特性研究中的烟气计算模型

20、李运泽,杨献勇(清华大学 热能工程系,北京 100084)摘 要:为了简便、准确、快速地解决大容量锅炉汽轮机组长期动态特性研究中的烟气放热计算问题,本文针对现有仿真用烟气放热模型的主要不足,由大容量锅炉热力计算的杜卜斯基-卜劳赫公式推导出新的仿真与动态特性分析用炉膛烟气放热模型;并通过简要的分析、简化和推导,得出了动态过程中各种锅炉换热面烟气出口温度和放热量计算方法,这些烟气计算模型能够较为准确、简单地反映动态过程中高温烟气在大容量锅炉的炉膛和各烟道换热面的烟温及放热量的变化规律,应用十分方便。最后将本文的烟气计算模型应用于某600MW超临界机组的动态特性仿真,仿真结果合理、正确。关键词:大容

21、量锅炉;长期动态特性;仿真;烟气计算模型中图分类号:TM621 文献标识码:A1 引 言锅炉汽轮机组的动态特性分析与仿真离不开锅炉的烟气放热量计算13。目前烟气计算使用的炉膛模型有两种:一种是文献23中介绍的基于四次方定律的四次方模型,这种模型存在难以合理确定炉膛烟气代表温度,计算结果也不能较为精确地反映炉膛放热量随燃料量变化的变化规律;另一种是基于前苏联锅炉热力计算标准中古尔维奇公式3的非四次方模型,这一模型能较为简便、准确地反映中小型锅炉炉膛出口烟温和放热量的变化规律,但由于古尔维奇公式应用在大容量锅炉上会出现较大误差45,因而在大容量锅炉动态特性分析和计算中采用这一模型会产生较大的误差。

22、文献45中研究表明大容量锅炉的炉膛尺寸较大,炉膛内烟气温度分布的不均匀度加大,在其炉膛出口烟温的热力计算中采用充分考虑上述因素并拟合了大量实际机组统计数据的杜卜斯基-卜劳赫公式较古尔维奇公式和四次方定律更为准确合理。因而推导出基于杜卜斯基-卜劳赫公式的仿真用炉膛烟气模型会提高大容量锅炉动态特性计算的精度。本文在总结前人研究15的基础上,由大容量锅炉炉膛热力计算的新方法(杜卜斯基-卜劳赫公式45)推导出新的仿真与动态特性分析用炉膛烟气放热模型,其仿真实验和分析表明这一新的非四次方模型能较为全面、合理地反映动态过程中大容量锅炉炉膛出口烟温与放热量的变化规律,应用也十分简单。参考文献:1BOSEJ

23、E,PARKERJ D,MCQUISTON F C.Designdata manual forclosed2loop ground2coupled heat pump systemM.Atlanta:AmericanSociety of Heating,Refrigerating and Air2Conditioning Engineers,1985.2 柳晓雷,王德林,方肇洪.竖直埋管地热换热器的传热模型与计算J.建筑热能通风空调,2001,21(2):1-3.3 崔 萍,刁乃仁,方肇洪.地热换热器间歇运行工况分析J.山东建筑工程学院学报,2001,16(1):52-57.4 曾和义,刁乃仁

24、,方肇洪.竖直埋管地热换热器的稳态温度场分析J.山东建筑工程学院学报,2002,17(1):1-6.5ECKERT E R G,DRAKE R M.Analysis of heat and mass transferM.New Y ork:Mc Graw2Hill,1972.104-106.6ROHSENOW W M.Handbook of heat transferM.New Y ork:Mc Graw-Hill,1973.7 杨世铭,陶文铨.传热学M.(第三版).北京:高等教育出版社,1998.8 张洪济.热传导M.北京:高等教育出版社,1992.9 钱滨江.简明传热手册M.北京:高等教育

25、出版社,1984.10PHETTEPLACE G,SULLIVAN W.Performance of a hybrid ground2coupled heat pump systemJ.ASHRAE Transactions,1998,104(1B):763-770.(何静芳 编辑)第18卷(总第104期)热能动力工程2003年3月 1995-2005 Tsinghua Tongfang Optical Disc Co.,Ltd.All rights reserved.neering,Shanghai Jiaotong University,Shanghai,China,Post Code:2

26、00240)Journal of Engineering for Thermal En2ergy&Power.-2003,18(2).163-165A numerical simulation was conducted of the particle movement in a fluidized bed with two2dimensional non2uniformlydistributed air.In addition to dealing with a gas2phase field by the use of Euler method a Lagrangian method is

27、 employedto treat a discrete particle field,directly keeping track of each particle in the particle field.The results of simulation in2dicate that there exists in the fluidized bed with non2uniformly distributed air an inner circulation movement of the parti2cles.Hence,the mixing characteristics of

28、the particles in the fluidized bed under discussion are superior to those in aconventional fluidized bed with uniformly distributed air.Key words:non2uniformly distributed air,fluidized bed,dis2crete element method地源热泵竖直埋管的有限长线热源模型=A Model of Finite2length Linear Heat2source for the Vertical Em2bedd

29、ed Pipe of a Ground2source Heat Pump 刊,汉ZENG He2yi,DIAO Nai2ren,FANG Zhao2hong(ResearchInstitute of Ground2source Heat Pumps under the Shangdong Institute of Architectural Engineering,Jinan,China,PostCode:250014)Journal of Engineering for Thermal Energy&Power.-2003,18(2).166-169Analyses and discussi

30、ons were conducted of a non2steady heat transfer model for the vertical embedded pipe of a geother2mal heat exchanger.With the use of a virtual heat source and Green function method obtained is an analytical solution ex2pression for the non2steady temperature field generated by a finite2length linea

31、r heat source in semi2infinite large media.By way of comparison with a steady2state temperature field solution discussed is the time required for the temperature fieldto attain a nominal”steady state”.Meanwhile,an analysis is performed of the temperature field when it has reached asteady state.In th

32、is connection a mistake that appeared in current textbooks was indicated.Two representative steady2state borehole wall temperatures,i.e.,the temperature at the middle of the borehole and the integral mean temperaturealong the borehole,are defined.A comparison of the difference between these two temp

33、eratures has led to a simplifiedcalculation formula,suitable for engineering applications.On the basis of the above analyses discussed further is the im2pact of the annual imbalance between heating and cooling loads of geothermal heat exchangers on their long2term perfor2mance.Key words:ground2sourc

34、e heat pump,geothermal heat exchanger,heat conduction,heat transfer model大型锅炉长期动态特性研究中的烟气计算模型=Flue2gas Calculation Model Used in the Study of Long2termDynamic Characteristics of a Large2sized Boiler 刊,汉LI Yun2ze,Y ANGXian2yong(Department of Thermal En2ergy Engineering,Tsinghua University,Beijing,Chi

35、na,Post Code:100084)Journal of Engineeringfor Thermal Ener2gy&Power.-2003,18(2).170-172The studyof long2term dynamic characteristicsof a large2sized turbogenerator set requires a simplified,accurate and rapidcalculation of flue2gas heat release rate.T o eliminate the main deficiency of currently use

36、d flue2gas heat release simula2tion models the authors have deduced from Dybosky2Broch formula of large2sized boiler thermodynamic calculation a newflue2gas heat release model for the simulation and analysis of dynamic characteristics.Moreover,through a concise analy2sis,simplification and deduction

37、 obtained was a method for calculating flue2gas outlet temperatures and heat release ratesin a dynamic process for various boiler heat exchange surfaces.The recommended easy2to2use flue2gas calculation modeloffers a relatively accurate and simple expression of the variation mechanism of flue gas tem

38、perature and heat release rateof high2temperature gases in the furnace and various gas2pass heat exchange surfaces of a large2sized boiler during a dy2namic process.The flue2gas calculation model has been used to simulate the dynamic characteristics of a 600MW super2critical power generation unit wi

39、th rational and accurate simulation results being obtained.Key words:large2sized boil2er,long2term dynamic characteristics,simulation,flue2gas calculation model膜式水冷壁壁温影响因素的数值分析=Numerical Analysis of Factors H aving an Impact on the Wall Tempera2612热能动力工程2003年 1995-2005 Tsinghua Tongfang Optical Disc Co.,Ltd.All rights reserved.