“白屋顶计划”对降低城市热岛效应作用的研究.doc

1、参赛队号 # 22161“白屋顶计划”对降低城市热岛效应作用的研究1 引言在近地面等温线图上，郊区气温相对较低，而市区则形成一个明显的高温区，如同露出水面的岛屿，被形象地称之为“城市热岛”。城市热岛中心，气温一般比周围郊区高 左右，最高可达 以上。1C 6C图 1 城 市 热 岛 效 应 示 意 图城市热岛效应的形成主要有以下原因：(1) 城市下垫面特性的影响。城市内有大量的人工构筑物，如混凝土、柏油路面和各种建筑墙面等，改变了下垫面的热力属性。这些人工构筑物吸热快而热容量小，在相同的太阳辐射条件下，它们比自然下垫面如绿地、水面等升温快，因而其表面温度明显高于自然下垫面。(2) 人工热源的影响

2、。工厂生产、交通运输以及居民生活都需要燃烧各种燃料，每天都在向外排放大量的热量。(3) 城市里中绿地、林木和水体的减少。随着城市化的发展，城市人口的增加，城市中的建筑、广场和道路等大量增加，绿地、水体等却相应减少，缓解热岛效应的能力被削弱。(4) 城市中的大气污染。城市中的机动车、工业生产以及居民生活，产生了大量的氮氧化物、二氧化碳和粉尘等排放物。这些物质会吸收下垫面热辐射，产生温室效应，从而引起大气进一步升温。预防、降低城市热岛效应的方法有多种，建筑物淡色化以增加热量的反射就是一种简单、有效的方法，其中以“白屋顶”计划最为著名。早期的屋顶采用焦油涂层材料制成，目的在于防水。在冬季，黑色屋顶可

3、以增加热量，对建筑物来说不是坏事。但在其它季节特别是夏季，黑色的屋顶也加重了城市的热岛效应。参赛队号 # 22162为此，纽约哥伦比亚大学的斯图尔特 加芬(Stuart Gaffin)开展了一项关于屋顶颜色对温度影响的研究。他发现，白色屋顶可以降低城市的热岛效应。2011 年 7 月 22 日，遭遇热浪袭击的纽约城市用电量突破纪录。研究显示：在这最热的一天中，纽约市一些深色的、具有吸光表面的传统屋顶的测量温度达到了 170 华氏度(约 76.7 摄氏度)；而一种白色屋面材料，其测量温度相对要低大约 42 华氏度(23.3 摄氏度)。图 2 “白 屋 顶 ”计 划 示 意 图加芬的研究表明，在

4、2011 年整个夏季，这种白色屋顶的温度峰值比典型的黑色屋顶平均降低了 43 华氏度(23.9 摄氏度)，这是第一个白色屋顶可以减少城市热岛效应的科学结论。相关研究论文已发表在 2012 年 3 月 7 日环境研究快报网络版上。但是，对于“白屋顶”计划也有一些反对意见。斯坦福大学土木与环境工程学教授马克.雅各布森声称通过计算机模拟显示：白屋顶实际上并不能冷却城市地表，相反会导致全球变暖。原因在于：白屋顶提高了空气稳定性，使云量减少，而云量的减少会导致更多的阳光照射到地面。导致热岛效应的原因众多，人们也先后提出了研究热岛效应的多种方法，大致可分为气象资料分析法、布点观测法、遥感法和数学模型模拟法

5、 等。14根据题意，本文选用建筑群热时间常数(Cluster Thermal Time Constant，CTTC)模型 ，从黑白屋顶对阳光的吸收和对空辐射的差异角度，57对“白屋顶”计划对降低城市热岛效应方面的作用进行了一些研究。2 数学模型的建立以色列环境学家 H. Swaid 和 M. E. Hoffman 在研究太阳辐射所引起的温升参赛队号 # 22163效果时，将建筑群简化成为周期性起伏的“城市峡谷” ，提出了所谓 CTTC 模型，即把特定地点的温度视为几个单独因素作用的叠加5(1)absollwTttTt其中， 为所研究的建筑群空气温度； 为基准背景温度，通常选为郊区的aTt b日

6、平均温度； 为吸收太阳辐射所引起的温升； 为对天空的长波soltlwt辐射引起的温降。CTTC 模型提出后，许多学者对其进行了一系列的改进和应用研究 。6,72.1 吸收太阳辐射引起的温升的计算水平地面所接受的太阳辐射是引起建筑群温度变化的主要因素。假定城市覆盖层所吸收的太阳辐射热释放给空气，并服从指数衰减规律，Swaid 和Hoffman 给出的温升变化公式为 5(2)01exptsol penkmtkTtIhCT其中， 为下垫面对太阳辐射的吸收率，如沥青为 0.800.95，土壤为m0.600.90； 为综合导热系数； 为单位面积上太阳直辐射强度的平均步长hpenI变化；CTTC 为建筑群

7、的热时间常数。2.2 CCTC 参数的确定在现代城市的建设中，由于使用了一些建筑材料，如混凝土和沥青等，而大大改变了城市地区的热平衡。它们不仅改变了原有的自然景观，而且能在较短的时间内储存比同体积土壤更多的能量。由于这些密封性很强的材料代替了原来的植物带，下雨天不容易留住水分，因此水分蒸发吸热及植物蒸腾作用吸热均减小了。CTTC 参数反映了结构的蓄热能力和透热能力，它具有时间的量纲，称为“热时间常数” 。因此，预测模型也常称为 CTTC 模型，它的值可通过将所有活动面的热时间常数加权平均获得 ：5(3)1groundwalFAWACTCTTCSS其中， 为相对于小区面积的建筑面积； 为相对于小

8、区面积的外墙面积；FAS为地板的热时间常数，约为 8 小时； 为墙壁的热时间常数，groundCT wal约为 6 小时。2.3 单位面积上太阳直射平均辐射强度的确定参赛队号 # 22164CTTC 模型认为建筑群受到的太阳辐射是造成建筑群温度波动的主要因素。由于在建筑小区中，建筑物投下阴影，对太阳光形成了遮挡作用，并且这种遮挡效果在一天中随着太阳高度角的改变而发生变化。因此有必要对太阳辐射强度进行修正，求出平均的太阳辐射强度。公式如下所示 ：5(4)1penItIPSAt其中， 为时间 时无遮挡的太阳直射辐射强度； 为时间 时的遮挡面Itt tt积率。2.4 导热系数的确定综合导热系数包含了

9、辐射和对流的共同作用。对流导热系数主要与沿下垫面的风速有关。城市覆盖面下表面综合导热系数可由下式计算：(5)rch(6)5.841cUt其中， 为低温表面的发射率，约为 0.9； 为辐射导热系数，平均表面温度为 r时， 约为 0.57； 为对流导热系数； 为下垫面附近的风速。20Crhcht2.5 对天空长波辐射引起的温降的计算地表对天空的长波辐射包括由于地表的高温自身发射的长波辐射和地表反射的部分环境辐射，主要是大气逆辐射。计算长波辐射的方法有多种，常用的一个经验公式为 8(7)41lussldRTR其中， 为总长波辐射， 为大气逆辐射， 为地表比辐射率， 为地表温luRld sT度， 为转

10、换系数。在长波波段，一般地表的比辐射率均在 0.90 以上，而地表的大气逆辐射相对很小，即长波辐射主要由地表自身发射的长波辐射构成。可以通过测量地表温度、地表比辐射率和转换系数近似计算得到地表的长波辐射，但前两者的测量具有很多不确定性而难于精确测量 。8有学者应用统计分析的方法研究了长波辐射与地表温度的相互关系，结果显示长波辐射与地表温度呈较严格的线性关系 。8基于上述结论，本文提出计算屋顶对天空长波辐射引起的温降的方法如下：(1) 根据前面计算出的 获得地表近似温度数据；solTt(2) 利用文献 8中长波辐射和地表温度的线性关系计算出长波辐射数据；参赛队号 # 22165(3) 与 Swa

11、id 和 Hoffman 给出的温升变化公式(2)原理类似，假定长波辐射引起的温降近似服从下列指数衰减规律(8)01exptlwspenk tkTtICT其中， 为地表比辐射率，通常取为 0.90； 为单位面积上长波辐射强度的s penI平均步长变化；CTTC 为建筑群的热时间常数。因为吸收太阳辐射引起的温升与对天空长波辐射引起的温降机理完全一致，所以我们认为做上述假设是合理的。3 模型的模拟计算下面根据专业手册查阅到的参数和相关文献提供的某城市的具体数据 ，6对 CTTC 模型进行模拟计算。基准温度 ，沥青或混凝土屋顶对太阳辐射的吸收率 ，26.3bTC 0.9m白屋顶对太阳辐射的吸收率 ，

12、建筑面积比 ，屋顶面积比0.65m0.4FAS，地板的热时间常数 ，墙壁的热时间常数0.4WAS 8groundT，遮挡面积率 ，表面发射率 ，辐射导热系数walC.1PSAt .9，下垫面附近的风速 ，各时段太阳辐射强度如下：.57rh4U表 1 各 时 段 太 阳 辐 射 模 拟 数 据时刻 6:00 7:00 8:00 9:00 10:00 11:00辐射强度 213.6 417.7 607.9 771.3 896.7 975.512:00 13:00 14:00 15:00 16:00 17:00 18:001002.4 975.5 896.7 771.3 607.9 417.7 21

13、3.63.1 屋顶吸收太阳辐射引起的温升的计算根据上述数据，利用 Maple 编程，可以求出黒( )和白( )屋0.9m0.65顶吸收太阳辐射所引起的温升。Maple 程序见附录 1，相应结果图见下( 红色为黒屋顶，蓝色为白屋顶 )：参赛队号 # 22166图 3 黒 白 屋 顶 吸 收 太 阳 辐 射 所 引 起 的 温 升 结 果 图需要注意的是，计算结果是屋顶表面的温升，而不是屋顶附近空气的温升。从图 3 中可以看出：不同吸收率的黑白屋顶的温升效果差异比较明显；这种差异从早晨 6:00 开始逐渐变大，下午 3:00 左右达到最大，然后逐渐变小；差异的最大值高达 50%。3.2 屋顶对天空

14、的长波辐射引起的温降的计算根据 2.5 中给出的计算方法，对长波辐射引起的温降进行模拟计算。取屋顶比辐射率 ，其余参数同上，长波辐射和地表温度的线性关系0.9s采用文献8中根据实测数据拟合出的结果： 。6.329lwRTMaple 程序见附录 2，相应结果图见下：图 4 屋 顶 对 空 长 波 辐 射 所 引 起 的 温 降 结 果 图从图 4 中可以看出：在相同的参数条件下，屋顶对天空长波辐射引起的温降幅度明显低于屋顶吸收太阳辐射引起的温升。3.3 太阳辐射对屋顶温度的综合影响参赛队号 # 22167为了准确反映同时考虑吸收太阳辐射引起的温升和长波辐射引起的温降后屋顶温度的变化情况，下面给出

15、太阳辐射对屋顶温度的综合影响效果图：图 5 太阳辐射对屋顶温度的综合影响结 果 图图中红、蓝色分别为黒、白屋顶的结果。图 5 与图 3 非常相似，即长波辐射引起的温降对屋顶温度的变化影响不大，我们可以得出与 3.1 中同样的结论。 4 结论CTTC 模型模拟计算结果清楚地表明：尽管由于白屋顶对天空的长波辐射率远高于黒屋顶，从而会导致较小的温降，但白屋顶和黒屋顶的温降差与两者的温升差相比是微小的，换句话说，所研究的建筑群空气温度 主要取决于aTt吸收太阳辐射所引起的温升 ，对天空的长波辐射引起的温降 对solTtlwt其贡献甚微，即“白屋顶”计划对降低城市热岛效应有一定程度的作用。需要指出的是，

16、本文仅用 CTTC 模型从黒白屋顶吸收和辐射率不同这一角度研究了“白屋顶”计划对城市热岛效应影响，方法、模型及结论不可避免地带有一定的局限性。从近年的研究趋势看，随着遥感技术的发展，利用红外遥感技术进行城市热环境的研究将呈主流趋势。另一方面，遥感与气象方法相结合的城市热环境机理研究和定量研究也将成为重要研究方向。参考文献1 但尚铭 , 许辉熙, 叶强等 . 我国城市热岛效应研究方法综述J. 四川环境, 27(4): 88-91, 2008.2 桑建国. 城市热岛效应的分析解J. 气象学报, 44(2): 251-255, 1986.3 边海 , 铁学熊. 天津市夜间热岛的数值模拟 J. 地理学

17、报, 43(2): 150-158, 1988.参赛队号 # 221684 李小风. 热岛效应强迫下的中尺度环流的动力特征及极限风速的一种解析表达J. 地理学报, 48(3): 327-335, 1990.5 Swaid H, Hoffman M E. Prediction of urban temperature variations using the analytical CTTC model J. Energy and Building, 14: 313-324, 1990.6 Elnahls M M, Williamson J J. An improvement of the CTT

18、C model for prediction urban air temperature J. Energy and Building, 25: 41-49, 1997.7 王菲, 肖勇全 . 太阳辐射对不同建筑群产生温升效果的探讨J. 山东建筑工程学院学报, 19(1): 59-62, 2004.8 李超, 魏合理 , 徐青山等 . 一种基于地面气温的地表出射长波辐射的计算方案J. 激光与红外, 39(4): 399-402, 2009.附录1 屋顶吸收太阳辐射所引起的温升计算程序restart:Tb:=26.36:m:=0.9:Fa:=0.4:Wa:=0.4:CTTCg:=8:CTTCw:

19、=6:PSA:=0.1:epsilon:=0.9:hr:=0.57:U:=4:It:=213.6,417.7,607.9,771.3,896.7,975.5,1002.4,975.5,896.7,771.3,607.9,417.7,213.6:CTTC:=(1-Fa)*CTTCg+Wa*CTTCw:Ipen:=:for i from 2 to nops(It) doIpen:=op(Ipen),(Iti-Iti-1)*(1-PSA):od:h:=epsilon*hr+hc:hc:=5.8+4.1*U:参赛队号 # 22169T:=t-m/h*sum(Ipenk*(1-exp(-(t-k)/CT

20、TC),k=1t):L1:=:for i from 1 to nops(Ipen) doL1:=op(L1),i+5,T(i):od:L1;plot(L1,ytickmarks=4,labels=Time,Temperature,labeldirections=HORIZONTAL,VERTICAL);注：本程序计算、绘制的是 m=0.90 时的温升结果。2 屋顶对天空长波辐射所引起的温降计算程序restart:Tb:=26.36:epsilonr:=0.9:Fa:=0.4:Wa:=0.4:CTTCg:=8:CTTCw:=6:PSA:=0.1:epsilon:=0.9:hr:=0.57:U:=

21、4:T:=0,0.68,1.91,3.53,5.35,7.20,8.90,10.29,11.24,11.64,11.45,10.65:It:=:for i from 1 to nops(T) doIt:=op(It),6.33*(Ti+Tb)+296:od:CTTC:=(1-Fa)*CTTCg+Wa*CTTCw:Ipen:=:for i from 2 to nops(It) doIpen:=op(Ipen),(Iti-Iti-1)*(1-PSA):od:参赛队号 # 221610h:=epsilon*hr+hc:hc:=5.8+4.1*U:T:=t-epsilonr/h*sum(Ipenk*(1-exp(-(t-k)/CTTC),k=1t):L1:=:for i from 1 to nops(Ipen) doL1:=op(L1),i+5,T(i):od:L1;plot(L1,ytickmarks=4,labels=Time,Temperature,labeldirections=HORIZONTAL,VERTICAL);