城市结构包括建筑墙体、屋顶及道路,道路与建筑两侧的空气形成类似于峡谷的地貌特征称为城市峡谷。城市峡谷内部吸收或反射太阳辐射不同于平坦地面,城市峡谷内部某一面的反射辐射被其他面阻拦不能直接“逃出”峡谷,形成多重反射。而多重反射极大地增加城市峡谷太阳辐射吸收量,进而加剧城市热岛效应。
AIDA等[1-2]通过实验测试不同季节不同峡谷走向的城市峡谷反射率,结果表明,城市峡谷反射率随时间变化并与峡谷走向有关。基于此实验,研究人员采用数值模型研究城市峡谷高反射率墙体材料的隔热性能[3]、城市峡谷反射率[4]、日平均加权反射率[5]以及城市峡谷反射率日变化[6]。研究发现,城市峡谷反射率明显低于路面和建筑墙体的反射率。提高路面反射率是提高城市峡谷反射率的一种有效方法[7-8]。而路面反射率是可控的,其大小取决于路面材料及其龄期,通常在0.20~0.80之间。研究还表明,采用浅色骨料[9-10]、高反射率的路面铺筑材料[11-13]或者在路面涂反射涂料[14]的反射沥青路面比传统沥青路面反射更多的太阳辐射;采用白色填充料[15]、浅色水泥基体[16]或者白色面层技术[17]的反射水泥混凝土路面也比传统水泥路面反射更多的太阳辐射。这些研究结果证实反射路面能够有效地降低路面温度,但这些研究均测量平坦路面的温度,而提高城市峡谷路面反射率能否有效提高城市峡谷反射率尚未证实。
笔者基于城市峡谷不同结构面之间的多重反射建立城市峡谷反射率数值模型,并用实验观测值验证模型的可靠性。通过改变模型参数,查证路面反射率对城市峡谷反射率的影响,同时探讨反射路面能否向邻近建筑物反射可观的额外散射辐射。
城市路面和建筑墙体是否受太阳直接辐射与太阳位置有关。而太阳位置随太阳赤纬角、太阳天顶角、太阳方位角及太阳时角而变化。其中太阳赤纬角δ由式(1)计算[18]
式中:N为日序,即1月1日为N=1;12月31日为N=365。
太阳天顶角θ由式(2)计算[19-20]
式中:φ为地球表面观测者所处的纬度;ω为太阳时角。
太阳时角ω可由式(3)计算[19]
其中,一天内太阳时角在负值与正值之间变化,日出前为π/2,正午时分为0,日落后为-π/2。
太阳方位角γ可由式(4)计算[5, 21]
当式(4)中sinγ>1或sinγ为无穷小时,太阳方位角可由式(5)给出[5]
为便于计算,取时间序列为日出至日落的偶数,则时间序列日出至日落为θ=0。正午时分(θ=0)的反射率由正午时分模拟的反射率采用内插法求出。
根据式(1)~(5)可估算城市峡谷内某个面为阴面或阳面的时间。图 1是典型的城市峡谷图,如果不存在阳面墙体,那么太阳直接辐射在x轴上的阴影长度x0可由式(6)计算
式中:h为建筑物高度;γc 为峡谷走向,东西方向时为0,南北方向时为π/2。
当太阳入射角较小时(x0≥w),图 1(a)向阳墙体仅有一部分受到太阳直接辐射,而路面和背阴墙体没有受到太阳直接辐射,此时向阳墙体上受太阳直接辐射的长度z为
当太阳入射角较大时(x0<w),图 1(b)向阳墙体完全受太阳直接辐射,部分路面受太阳直接辐射。路面受直接辐射的宽度x为
白天,一面墙体一直为阴面,另一面墙体完全或部分为阳面。为方便起见,建议计算z和x的绝对值,估算tanθ和cos(γ-γc) 的值确定墙体和路面受直接辐射的部分。对于一个对称的城市峡谷,用绝对值足可估算城市峡谷的吸收量和反射量。
峡谷内的路面与建筑墙体受太阳直接辐射和散射辐射作用。根据Orgill和Hollands 的模型[22],直接辐射强度 Ib 和散射辐射强度 Id 可按式(9)计算
式中:I为瞬时太阳辐射强度,单位W/m2,I=Id+Ib;kT为天空清晰度,晴朗无云时取kT=1.0,阴天时取kT=0。
一天内瞬时太阳辐射强度I呈余弦分布,其中,日出前与日落后为0,正午时分最大[23]
式中:I0为单日最大太阳辐射强度[24]
式中:S0为太阳辐射常数,其值为1 367 W/m2;τ为常数,0.61≤τ≤0.81,晴天时τ取0.81。
视角因子或辐射角系数,表示离开某一表面的辐射直接投射到另一表面的百分数,与另一表面的吸收能力无关。视角因子是散射强度的决定因素,峡谷内各个面到天空的视角因子共同决定城市峡谷太阳辐射吸收量。路面到墙体和天空的视角因子可由式(12)计算[21](下标w表示路面,h表示墙体,s表示天空,下同。为便于计算,假设两侧墙体等高):
墙体到天空、路面以及对面墙体的视角因子分别为[21]
根据式(12)和式(13)可计算城市峡谷内部某一表面到另一表面的视角因子。某一表面到天空的视角因子决定到达其表面的太阳散射辐射,路面和墙体接收到的太阳散射辐射量Idw和Idh (两面墙体的太阳散射辐射量相等)由式(14)计算[23]
当太阳入射角较小时,图 1(a),墙体上仅有高度为z的表面受到太阳直接辐射,其他面为阴面。z到天空、对面墙体以及路面的视角因子分别为[21](z到自身墙体的视角因子为0)
当太阳入射角较大时(图 1(b)),路面部分甚至全部受太阳直接辐射。假设宽度为x的路面受到太阳直接辐射,x到天空、邻近x的墙体以及远离x的墙体的视角因子分别为[21](受辐射部分到未受辐射部分的视角因子为0)
多种数值模型模拟太阳辐射在城市峡谷内的多重反射[4, 22, 25]。笔者提出一种较为简单的方法计算城市峡谷内部的多重反射。假设第i次多重反射左边墙体、右边墙体以及路面的反射量分别为Rlh,i、Rrh,i和Rw,i,这些反射除被自身以外的表面吸收,还被天空吸收,未被吸收的反射量变成对应的第i+1次反射量Rlh,i+1、Rrh,i+1和Rw,i+1,如式(17)
式中:ρh和ρw分别为墙体与路面的反射率。
第i+1次左边墙体、右边墙体、路面及天空的辐射吸收量分别为Alh,i+1、Arh,i+1、Aw,i+1及As,i+1,如式(18)
迭代过程中,若式(18)中的Alh,i+1、Arh,i+1、Aw,i+1及As,i+1的最大值分别是∑i=ni=1(Alh,i+1+Arh,i+1+Aw,i+1+As,i+1)的无穷小量,那么后续多重反射量可以忽略不计,迭代停止。
第一次反射循环时,墙体或部分路面为阳面,受太阳直接辐射的表面比未受太阳直接辐射的表面反射更多的太阳辐射。阳面到其它表面或天空的辐射包括2个部分:整个表面的散射辐射和受照射部分的直接辐射。这些反射量和相应的吸收量的初始值为Rlh,1、Rrh,1、Rw,1、Alh,1、Arh,1、Aw,1及As,1。由于直接辐射随太阳位置而变化,因此,计算时需明确太阳位置。
当太阳入射角较小时,部分墙体为阳面,因此,初始值计算如下:
式中:Iz=wIb√w2+z2为墙体z部分的直接辐射强度。右边第1项为不同表面对应的太阳直接辐射吸收或反射量,第2项为太阳散射辐射相应的吸收或反射量。式19(b)右边第1项为零,因为这部分反射量叠加到其他表面的吸收量与反射量以及天空的吸收量式(19c)~(19g)。
当太阳入射角较大时,部分路面为阳面,因此,初始值计算
式中:Iz=(w−x)Ib√h2+(w−x)2,Ix=hIb√h2+(w−x)2。式20(f)中直接辐射的反射量为零,因为这部分反射量叠加到其它表面的吸收量与反射量以及天空的吸收量式(20a)~(20d)、(20g)。
由于式(17)~(20)计算反射量与吸收量较为繁琐,本文采用求和法。计算每次迭代过程中反射量与吸收量,验证其与城市峡谷内部的太阳辐射量是否相等。
根据城市峡谷内的太阳总辐射量及峡谷向天空反射的辐射量,城市峡谷的反射率可由式(21)计算。
式中:n为迭代次数。式(21)表明城市峡谷反射率ρc与建筑物高度h、路面宽度w、墙体反射率ρh、路面反射率ρw、路基走向γc、太阳位置、时间、直接辐射以及散射辐射等参数有关。通过式(17)~(20)的关系式可估算参数的关系,进而计算出城市峡谷反射率。
城市峡谷反射率计算的流程图如图 2所示。
采用文献[1]中的试验观测值验证数值模型的可靠性,该试验测量人造城市峡谷反射率。人造城市峡谷模型用水泥砖做建筑物墙体与屋顶,水泥砖墙体之间的部分做城市路面。试验所用人造城市峡谷的纵横比为1.0,但峡谷走向不同。试验观测地纬度为35.28°,晴朗天气。
为验证模型,将墙体和屋顶反射率和城市峡谷路面反射率叠加即为观测反射率。鉴于墙体和路面同为水泥混凝土材料,假定墙体与屋顶反射率ρh及路面反射率ρw为测量平坦混凝土地面的反射率 (如图 3)。用线性内插法计算特定时间的平坦混凝土地面的反射率。图 3对比不同季节的试验观测值与模拟值,试验观测值与模拟值的偏差较小,表明数值模型具有可靠性,而该偏差可能是因为试验观测地不确定的气象条件与人为因素造成的。
城市峡谷反射率依赖于纵横比、墙体和路面的反射率、太阳位置及城市峡谷走向。为简化问题,本次模拟假设平坦混凝土地面的反射率为常数。以广西南宁(φ=22.8°)为例,除有特别说明外,本次模拟均采用以下参数值:h/w=1.0、ρh=0.4、ρw=0.4、kT=1.0及N=173(夏至)。
浅色墙体的实际反射率为0.40,灰黑色路面的反射率为0.10。不同走向和纵横比的城市峡谷反射率变化情况如图 4所示。根据图 4,城市峡谷反射率呈M型,日出后反射率迅速增大,之后达到最大值,正午时分反射率最小。
城市峡谷反射率随太阳位置而变化。当太阳入射角较小时,散射辐射占主导作用,且整个峡谷均受散射辐射作用,导致城市峡谷反射率较低。其根本原因是散射辐射到达峡谷后发生多重反射,与平坦地面相比峡谷反射到天空的辐射较少;随着太阳入射角增大,直接辐射占主导作用,太阳直接辐射照射到墙体上,一部分太阳辐射被反射到路面及对面墙体,另一部分被反射到天空。由于受直接辐射的墙体对天空的视角因子小于1/2,城市峡谷反射率的最大值小于0.2(0.4/2)(如图 4)。峡谷真实的最大反射率为0.04~0.06,小于0.2,这是因为散射辐射的低反射率抑制了城市峡谷反射率。最大反射率出现在太阳直接辐射照射到路面时,路面受辐射的部分也随着时间的推移而减少,路面的反射减少导致城市峡谷反射率降低。
除太阳位置外,时间与路基走向是影响城市峡谷反射率另两个重要因素。对于东西走向的城市峡谷,其反射率冬季较夏季大,原因在于冬季太阳直接辐射照射到墙体的时间较夏季长;而对于南北走向的城市峡谷,其反射率夏季较冬季大,原因在于太阳直接辐射首先照射到东面墙体,然后是路面,最后是西面墙体。
对东西走向的城市峡谷,最大辐射量出现在日出后不久,因为日出后路面受直接辐射(如图 4(a)和4(b));而对南北走向的峡谷,太阳直接辐射需要一定的时间才能照射到路面,因此,城市峡谷反射率最大值出现在日出后和日落前的1~2 h(如图 4(c)和4(d))。
纵横比是影响城市峡谷反射率的关键因素,纵横比大的城市峡谷反射率小。深峡谷内部某一表面反射的辐射被其他面吸收的可能性比浅峡谷大,从而导致深峡谷太阳辐射吸收量比浅峡谷大。根据图 4,正午时分纵横比为4.0时的城市峡谷反射率较纵横比为0.25或1.0时低,说明墙体吸收较多的路面反射辐射。
路面反射率的大小取决于路面材料。养护期的灰色水泥混凝土路面反射率为0.35~0.40,老化的水泥混凝土路面反射率为0.20~0.30,反射率减小的原因在于轮胎和路面的摩擦及路面掺杂的泥土使路面颜色加深。新铺筑的沥青混凝土路面反射率约为0.05,而老化的沥青混凝土路面反射率为0.10~0.15,其原因在于路面磨损骨料外露路面颜色变浅致使反射率变大。为研究不同路面反射率对城市峡 谷反射率的影响,本文模拟峡谷墙体反射率为0.40时,不同路面反射率(ρw=0.10、0.20、0.30、0.40)的城市峡谷反射率。根据图 5(b),冬季东西走向的城市峡谷反射率较稳定,路面反射率的变化对其影响较小,这是因为太阳直接辐射主要照射南面墙体,较少到达路面。 图 5表明,提高路面反射率使城市峡谷反射率曲线由M型转变为W型。对于低反射率的路面,正午时分,城市峡谷反射率降低到最小,原因在于正午时分整个路面受太阳直接辐射,较低的路面反射率使其吸收较多的太阳辐射,从而降低城市峡谷反射率,导致城市峡谷反射率日变化曲线呈M型;对于高反射率的路面,由于路面相对于墙体有较大的天空的视角因子,路面受太阳直接辐射时,更多的太阳辐射被反射到天空。因此,正午时分城市峡谷反射率有所增加,城市峡谷反射率日变化曲线呈W型。
城市峡谷反射率随着路面反射率的增大而增大。正午时分,路面的反射率每增加0.10,城市峡谷反射率增加大约0.05,其原因可能是纵横比为1.0时,路面对天空的视角因子约为0.5。这表明,天空视角因子或者纵横比是控制城市峡谷内部太阳辐射吸收量的主要因素。路面反射率对城市峡谷反射率的影响依赖于峡谷走向,东西走向的城市峡谷反射率日变化不同于南北走向。
使用冷却路面缓解城市热岛效应获得广泛关注。夏季反射路面向邻近建筑物反射太阳辐射,增加邻近建筑物的热存储。为估算存储热量,本文模拟不同路面反射率的城市峡谷内部建筑物表面与天空的太阳辐射吸收量。
图 6为夏季东西走向峡谷的路面与墙体的太阳辐射吸收量。图 6(c)表明,正午时分,路面的反射率每提高0.10,其相应的吸收量减少100 W/m2。其中,该减少量的一半被反射到天空,另一半分别被2个墙体吸收,吸收量均为25 W/m2(如图 5(a)和5(b)),这种分配比例是由于路面对天空的视角因子接近0.5(式(12)中h/w=1.0时, Fw→s=√2)。由于纵横比决定视角因子,这种太阳吸收量分配比例进一步说明纵横比对城市峡谷反射率起决定作用。
南北走向的峡谷具有同样的太阳吸收量分配模式(如图 7),但具体的能量分配过程不同。对南北走向的峡谷,随着时间变化太阳辐射先后照射东面墙体、路面以及西面墙体。尽管存在这种差别,但天空仍吸收路面辐射量的1/2,每面墙体各吸收1/4。
反射路面向邻近建筑物反射太阳辐射,促使建筑物存储热量。本文分别模拟峡谷东西走向与南北走向时峡谷内部不同面与天空的吸收量(图 8、9)。图 8表明,南面墙体比北面墙体和路面吸收较多的太阳辐射,因为太阳方位角远离路面方位角,导致南面墙体受太阳直接辐射的时间长,而路面受太阳直接辐射的时间短。路面的反射率每增加0.10,路面向天空和墙体反射的太阳辐射约增大15 W/m2,每面墙体总的额外吸收辐射大约为4 W/m2(15/2/2)。图 9反映峡谷南北走向时峡谷内部表面与天空吸收的额外辐射。南北走向的路面比东西走向的路面的照射时间更长,路面散射更多的热量。正午时分,路面的反射率每增加0.10,路面对天空与墙体的额外辐射量大约为8~10 W/m2。
提出城市峡谷反射率的理论计算模型,通过对比分析人造城市峡谷反射率试验观测值和理论计算值,验证了该数值模型的可靠性。分析模型得到如下结论:
1) 峡谷纵横比是影响城市峡谷反射率的决定性因素,而其它因素影响相对较小。当峡谷纵横比<1.0时,提高路面的反射率才能有效提高城市峡谷反射率。
2) 城市峡谷内部的反射路面能够向邻近建筑物反射可观的额外散射辐射。冬季时额外的太阳辐射反射到路面的反射量仅为夏季的1/5。东西走向峡谷的反射率与南北走向峡谷反射率没有明显区别。当峡谷纵横比<1.0时,建议峡谷路面采用反射路面缓解城市热岛。