METHOD FOR BUILDING URBAN CANOPY MODEL BASED ON TROPICAL ISLAND CLIMATE CHARACTERISTICS

Abstract

Provided is a method for building an urban canopy model based on tropical island climate characteristics. Adjacent regions are linked together, multiple streets of finite lengths within the regions affecting each other. Net radiation heat flux Q*=long-wave radiation+short-wave radiation. The tropical island urban canopy model considers the strong solar radiation and high temperature and high humidity climate characteristics of tropical cities and the influence of perennial monsoons on island cities, improves the methods of processing long-wave radiation flux, short-wave radiation flux, sensible and latent heat flux, street canyon wind velocity, heat storage flux, anthropogenic heat flux and horizontal heat flux on the basis of an urban canopy model, and has higher adaptability to the studies on the tropical island-type urban heat island effect.

Claims

1. A method for building an urban canopy model based on tropical island climate characteristics, wherein the method comprises: linking adjacent regions together, multiple streets of finite lengths within the regions affecting each other; wherein an energy balance equation in a street is as follows:
Q*.sub.s+Q.sub.F,s=Q.sub.H,s+Q.sub.E,s+ΔQ.sub.S,s+ΔQ.sub.A,s  Equation 1 wherein Q* is a net radiation heat flux, in unit of W/m.sup.2; Q.sub.F is an anthropogenic heat production, in unit of W/m.sup.2; Q.sub.H is a sensible heat flux, in unit of W/m.sup.2; Q.sub.E is a latent heat flux, in unit of W/m.sup.2; ΔQ.sub.S is a net heat storage flux, in unit of W/m.sup.2; ΔQ.sub.A is a net convective heat flux, in unit of W/m.sup.2; s is the street; wherein net radiation heat flux Q*=a long-wave radiation+a short-wave radiation.

2. The method according to claim 1, wherein the long-wave radiation comprises: a net long-wave radiation of a pavement L.sub.r* is:
L.sub.r*=ϵ.sub.rΨ.sub.rL.sup.↓−ϵ.sub.rσT.sub.r.sup.4+ϵ.sub.rϵ.sub.w(1−Ψ.sub.r)σT.sub.w.sup.4+ϵ.sub.r(1−ϵ.sub.w)(1−Ψ.sub.r)Ψ.sub.wL.sup.↓+ϵ.sub.rϵ.sub.w(1−ϵ.sub.w)(1−‥.sub.r)(1−2Ψ.sub.w)σT.sub.w.sup.4+ϵ.sub.r(1−ϵ.sub.w(1−Ψ.sub.r)Ψ.sub.wσϵ.sub.rT.sub.r.sup.4 and a net long-wave radiation of a wall L.sub.w* is:
L.sub.w*=ϵ.sub.wΨ.sub.wL.sup.↓−ϵ.sub.wσT.sub.w.sup.4+ϵ.sub.wΨ.sub.wσϵ.sub.rT.sub.r.sup.4+ϵ.sub.w.sup.2(1−2Ψ.sub.w)σT.sub.w.sup.4+ϵ.sub.w(1−ϵ.sub.r)Ψ.sub.rΨ.sub.wL.sup.↓+ϵ.sub.w(1−ϵ.sub.w)Ψ.sub.w(1−2Ψ.sub.w)L.sup.↓+ϵ.sub.w.sup.2(1−ϵ.sub.w)(1−2Ψ.sub.w).sup.2σT.sub.w.sup.4+ϵ.sub.w.sup.2(1−ϵ.sub.r)Ψ.sub.w(1−Ψ.sub.r)σT.sub.w.sup.4+ϵ.sub.w(1−ϵ.sub.w)Ψ.sub.w(1−2Ψ.sub.w)σϵ.sub.rT.sub.r.sup.4 where L↓ is an amount of solar radiation, σ is a standard deviation, T.sub.r and T.sub.w are temperatures of the pavement and the wall, and ϵ.sub.r and ϵ.sub.w are emissivities of the pavement and the wall; for the pavement, if Ψ.sub.r is a sky viewing angle coefficient of the pavement to the sky, amount of solar radiation, (1−Ψ.sub.r) is a sky viewing angle coefficient of the pavement to the walls on both sides; a sky viewing angle coefficient of the wall to the sky is Ψ.sub.w, a sky viewing angle coefficient to the pavement is Ψ.sub.w, then a sky viewing angle coefficient to the opposite wall is (1−2Ψ.sub.w), and the sky viewing angle coefficient is 1.0 for a roof; the sky viewing angle coefficient is calculated using a plane angle, the sky viewing angle coefficient at a w/2 position of the pavement is ${\Psi }_r\left(y=w/2\right)=1+\frac{1}{\pi }{\tan }^{-1}\left(\frac{h/w}{{\left(h/w\right)}^2-1/4}\right)\mathrm{for}h/w\le 1/2$ ${\Psi }_r\left(y=w/2\right)=1-\frac{1}{\pi }{\tan }^{-1}\left(\frac{h/w}{{\left(h/w\right)}^2-1/4}\right)\mathrm{for}h/w\ge 1/2$ the sky viewing angle coefficient at an intersection of the wall and the pavement is ${\Psi }_w\left(z=0\right)=\frac{1}{\pi }{\tan }^{-1}\left(\frac{1}{h/w}\right),$ where h represents a height of a street canyon, and w represents a width of the street canyon.

3. The method according to claim 1, wherein the short-wave radiation comprises: average direct solar radiant fluxes of the pavement , a west wall , an east wall and the roof are calculated according to a perpendicular angle of the street to a sun direction: $S_r^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\{\begin{array}{cc}0& \mathrm{for}.Math.\lambda .Math.>{\lambda }_0\\ \left[1-\frac{h}{w}\tan \left(\lambda \right)\right]S^{\Downarrow }& \mathrm{for}.Math.\lambda .Math.<{\lambda }_0\end{array}$ $S_{ww}^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\{\begin{array}{cc}\frac{1}{2}\frac{w}{h}{S}^{\Downarrow }& \mathrm{for}-\frac{\pi }{2}\le \lambda \le -{\lambda }_0\\ \frac{1}{2}\tan \left(\lambda \right)S^{\Downarrow }& \mathrm{for}-{\lambda }_0\le \lambda \le 0\\ 0& \mathrm{for}0\le \lambda \le \frac{\pi }{2}\end{array}$ $S_{we}^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\{\begin{array}{cc}0& \mathrm{for}-\frac{\pi }{2}\le \lambda \le -{\lambda }_0\\ \frac{1}{2}\tan \left(\lambda \right)S^{\Downarrow }& \mathrm{for}-{\lambda }_0\le \lambda \le 0\\ \frac{1}{2}\frac{w}{h}{S}^{\Downarrow }& \mathrm{for}0\le \lambda \le \frac{\pi }{2}\end{array}$ $S_R^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\chi S^{\Downarrow }$ where is direct solar radiation on a horizontal surface, θ is an angle between a sun angle and an axial direction of the canyon, λ is an angle between a sun direction and a normal direction of the wall, χ is a ratio of a direct radiation to a total radiation at a top of the street canyon, h is a height of the street canyon, and w is a width of the street canyon; according to an orientation change of the street canyon, the width w of the street canyon is corrected as w/sin θ; after a heat flux of the wall is obtained, the heat flux of the wall is multiplied by sin θ for correction, where θ.sub.0 is an orientation of the street canyon where the pavement does not receive direct sunlight at all ${\theta }_0=\arcsin \left(\min \left(\frac{w}{h}\frac{1}{\tan \lambda };1\right)\right)$ all direct radiant fluxes obtained by the street canyon are averaged according to all possible changes in the direction of the street canyon, wherein two integrals are used, one between θ=0 and θ=θ.sub.0, and the other between θ=θ.sub.0 and $\theta =\frac{\pi }{2};$ wherein average direct solar fluxes of the wall, the pavement and the roof are: $S_r^{\Downarrow }=S^{\Downarrow }\left[\frac{2{\theta }_0}{\pi }-\frac{2}{\pi }\frac{h}{w}\tan \left(\lambda \right)\left(1-\cos \left({\theta }_0\right)\right)\right]$ $S_w^{\Downarrow }=S^{\Downarrow }\left[\frac{w}{h}\left(\frac{1}{2}-\frac{{\theta }_0}{\pi }\right)+\frac{1}{\pi }\tan \left(\lambda \right)\left(1-\cos \left({\theta }_0\right)\right)\right]$ $S_R^{\Downarrow }=S^{\Downarrow }$ S.sup.↓ is a scattered solar radiation available on the horizontal surface, and an amount of scattered solar radiation received by a surface in the street canyon is directly obtained from the sky viewing angle coefficient; due to an influence of a shape of the street canyon and high-reflectivity building surface materials, the short-wave radiation is calculated to solve a geometric system with an infinite number of reflecting surfaces, a reflecting processes of which are assumed to be isentropic processes; when direct and diffuse reflectivities of each surface are the same, an energy stored by the pavement A.sub.r and an energy stored by the wall A.sub.w when a first reflection occurs are:
A.sub.r(0)=(1−α.sub.r)(+S.sub.r.sup.↓)
A.sub.w(0)=(1−α.sub.w)(+S.sub.w.sup.↓) where α.sub.r and α.sub.w represent the reflectivities of the pavement and the wall, respectively; the energies of the reflected parts from the pavement R.sub.r and from the wall R.sub.w are:
R.sub.r(0)=α.sub.r(+S.sub.r.sup.↓)
R.sub.w(0)=α.sub.w(+S.sub.w.sup.↓) after n reflections occur,
A.sub.r(n+1)=A.sub.r(n)+(1−α.sub.r)(1−Ψ.sub.r)R.sub.w(n)
A.sub.w(n+1)=A.sub.w(n)+(1−α.sub.r)‥.sub.wR.sub.r(n)+(1−α.sub.w)(1−2Ψ.sub.w)R.sub.w(n)
R.sub.r(n+1)=α.sub.r(1−Ψ.sub.r)R.sub.w(n)
R.sub.w(n+1)=α.sub.w‥.sub.wR.sub.r(n)+α.sub.w(1−2Ψ.sub.w)R.sub.w(n) the following is obtained according to recursive formulas, $A_r\left(n+1\right)=A_r\left(0\right)+\left(1-{\alpha }_r\right)\left(1-{\Psi }_r\right)\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)$ $A_w\left(n+1\right)=A_w\left(0\right)+{\Psi }_w\left(1-{\alpha }_r\right)\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)+\left(1-2{\Psi }_w\right)\left(1-{\alpha }_w\right)\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)$ $\mathrm{and}$ $\underset{k=0}{\overset{n}{.Math.}}{R}_r\left(k\right)=\left(1-{\Psi }_r\right){\alpha }_r\underset{k=0}{\overset{n-1}{.Math.}}{R}_w\left(k\right)+R_r\left(0\right)$ $\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)={\alpha }_w{\Psi }_w\underset{k=0}{\overset{n-1}{.Math.}}{R}_r\left(k\right)+{\alpha }_w\left(1-2{\Psi }_w\right)\underset{k=0}{\overset{n-1}{.Math.}}{R}_w\left(k\right)+R_w\left(0\right)$ for the infinite reflections, the following is obtained by solving a geometric system, $\underset{k=0}{\overset{\infty }{.Math.}}{R}_r\left(k\right)=\frac{R_r\left(0\right)+\left(1-{\Psi }_r\right){\alpha }_r\left(R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)\right)}{1-\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}\underset{k=0}{\overset{\infty }{.Math.}}{R}_r\left(k\right)=\frac{R_r\left(0\right)+\left(1-{\Psi }_r\right){\alpha }_r\left(R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)\right)}{1-\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}$ M.sub.r is assumed to be a sum of a pavement reflection and M.sub.w is assumed to be a sum of a wall reflection, $M_r=\frac{R_r\left(0\right)+\left(1-{\Psi }_r\right){\alpha }_r\left(R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)\right)}{1-\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}{M}_w=\frac{R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)}{1\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}$ where
R.sub.r(0)=α.sub.r+α.sub.rS.sub.r.sup.↓
R.sup.w(0)=α.sub.r+α.sub.rS.sub.w.sup.↓ a total solar radiation absorbed by the pavement is S.sub.r*, a total solar radiation absorbed by the wall is S.sub.w.sup.*, a total solar radiation absorbed by the roof is S.sub.R*:
S.sub.r*=(1−α.sub.r)+(1−α.sub.r)S.sub.r.sup.↓+(1−α.sub.r)(1−Ψ.sub.r)M.sub.w
S.sub.w*=(1−α.sub.w)+(1−α.sub.w)S.sub.w.sup.↓+(1−α.sub.w)(1−2Ψ.sub.w)M.sub.w+(1−α.sub.w)Ψ.sub.wM.sub.r
S.sub.R*=(1−α.sub.R)+(1−α.sub.R)S.sub.R.sup.↓.

4. The method according to claim 1, wherein the anthropogenic heat production specifically comprises: a current anthropogenic heat flux in the street canyon is Q.sub.F=Q.sub.FV+Q.sub.FH+Q.sub.FM; where Q.sub.FV, Q.sub.FH and Q.sub.FM are heat generated by vehicles, fixed heat sources and biological metabolism, respectively.

5. The method according to claim 1, wherein the sensible heat flux Q.sub.H comprises:
Q.sub.H,r,ww,we=ρC.sub.pC.sub.H1U.sub.can(T.sub.r,ww,we−T.sub.can)
Q.sub.H,R=ρC.sub.pC.sub.H2U.sub.top(T.sub.R−T.sub.air)
Q.sub.H,can=ρC.sub.pC.sub.H2U.sub.air(T.sub.can−T.sub.air) where r, ww, we, R, and can refer to the pavement, the west wall, the east wall, the roof, and the street canyon, respectively; ρ is an air density; C.sub.p is a specific heat under a constant pressure; T.sub.can is a temperature in a center of the street canyon (w/2, h/2); U.sub.can and U.sub.top are a wind velocity in the center of the street canyon (w/2, h/2) and a wind velocity above the street canyon; U.sub.air and T.sub.air are an input wind velocity and an input temperature at a reference height of a turbulence model, and C.sub.H1 and C.sub.H2 are dimensionless velocity transfer coefficients; differences between the C.sub.H1 and C.sub.H2 are only a height and a roughness of a reference layer; the same zero plane layer and roughness are used, and the values of the two are equal, and are calculated as follows: $C_{H*}{u}_{*}=\frac{{\mathrm{ku}}_{*}}{{\Psi }_h}$ where k is a Von Karman constant, u.sub.* is a friction velocity of the reference layer, and Ψ.sub.h is a general integral function, ${\Psi }_h={\int }_{{\zeta }^T}^{{\zeta }^{\prime }}\frac{{\varphi }_h}{\zeta }d\zeta$ where ζ′=(z.sub.a−d)/L; ζ.sup.T=z.sub.T/L, z.sub.T is a roughness length of a heat flow; L is an Obukhov stability length, $L=-\frac{\rho C_p{\mathrm{Tu}}_{*}^3}{{\mathrm{kgH}}_a}$ where T is an average temperature of this layer, H.sub.a is an air flux between the street canyon and an atmosphere, and L is an implicit function, which is solved by simplified iteration; when a specific heat of air in an urban canopy is ignored, H.sub.a is a weighted average of a wall flux and a road flux in the street canyon, that is,
H.sub.a=2(h/w)Q.sub.w+Q.sub.R in a TEB model proposed by Masson, C.sub.H*u.sub.* is a reciprocal of aerodynamic resistance, i.e. 1/RES.sub.*, which is determined by the wind velocities in the street canyon and at a top of the street canyon; if a surface covered by plants such as green space is not considered, an average sensible heat flow of the street canyon depends on a weighted average area of the roof, the wall and the pavement, $Q_H=\frac{{\mathrm{bQ}}_{H,R}+{\mathrm{wQ}}_{H,r}+h\left(Q_{H,\mathrm{ww}}+Q_{H,\mathrm{we}}\right)}{w+b},$ wherein h is a height of the street canyon, w is a width of the street canyon, and b is an average width of the buildings.

6. The method according to claim 1, wherein the latent heat flux Q.sub.E comprises: a direct latent heat flow between a building roof and the atmosphere
Q.sub.E,R=l.sub.vB.sub.RρC.sub.H2U.sub.top(q.sub.R−q.sub.air) where l.sub.v is a latent heat of evaporation, B.sub.R is a humidity parameter of the roof, between 0 and 1, 0 is completely dry, 1 is completely wet, a value of B depends on the plant and a water conditions of the surface, ρ is a density of the air, C.sub.H2 is a dimensionless velocity transfer coefficient, and q.sub.R is a humidity of the roof; q.sub.air is a humidity at the reference height, a latent heat flow is calculated using a similarity law for the air on the pavement and the wall and in the street canyon
Q.sub.E,r=l.sub.vB.sub.rρC.sub.H1U.sub.can(q.sub.r−q.sub.can)
Q.sub.E,w=0 wherein C.sub.H2 is a dimensionless velocity transfer coefficient, and q.sub.r is a humidity of the pavement; q.sub.can is a humidity at the street canyon, a latent heat flow between an interior of the street canyon and a top atmosphere is
Q.sub.E,can=l.sub.vρC.sub.H2U.sub.air(q.sub.can−q.sub.air).

7. The method according to claim 1, wherein the net heat storage flux ΔQ.sub.S comprises: because there is a temperature gradient inside the building or the pavement, the roof, the wall and the pavement are assumed to be of at least three-layer structures; for an outermost layer structure, heat transfer equations of the three planes are written as, $C_{R1}\frac{\partial T_{R1}}{\partial t}=\frac{{\lambda }_{R1}}{d_{R1}}\left(S_R^{*}+L_R^{*}-H_R-{\mathrm{LE}}_R-G_{R1-2}\right)C_{w1}\frac{\partial T_{w1}}{\partial t}=\frac{{\lambda }_{w1}}{d_{w1}}\left(S_w^{*}+L_w^{*}-H_w-G_{w1-2}\right)C_{r1}\frac{\partial T_{r1}}{\partial t}=\frac{{\lambda }_{r1}}{d_{r1}}\left(S_r^{*}+L_r^{*}-H_r-{\mathrm{LE}}_r-G_{r1-2}\right)$ where T.sub.*i is a temperature of the i-th layer; C.sub.*i is a specific heat capacity of the air; d.sub.*i is a layer thickness, fluxes S.sub.**, L.sub.**, LE.sub.*, and G.sub.*1-2 are net solar radiation, net infrared radiation, sensible heat, latent heat, and thermal conductivity between the surface layer and the next layer, and the thermal conductivity is calculated using a Fourier heat conduction equation, $G_{*1-2}=\frac{{\lambda }_{*1-2}}{0.5\left(d_{*1}+d_{*2}\right)}\left(T_{*1}-T_{*2}\right)$ an average thermal conductivity between two adjacent layers λ.sub.*1-2 is calculated using a geometric average method: ${\lambda }_{*1-2}=\frac{d_{*1}+d_{*2}}{\left(d_{*1}/{\lambda }_{*1}\right)+\left(d_{*2}/{\lambda }_{*2}\right)}$ where λ.sub.*i is a thermal conductivity of the i-th layer; an inner first layer of the surface is assumed to be a very thin surface, and a temperature of the first layer is simplified to an outer surface temperature; for other inner i-th layers, a thermal conductivity between adjacent layers is calculated; for an innermost layer, such as the n-th layer, an internal temperature of the building is used for the roof and the wall surface, and the 0 flux is used for the pavement; $G_{\mathrm{Rn}-n+1}=\frac{{\lambda }_{\mathrm{Rn}}}{0.5d_{\mathrm{Rn}}}\left(T_{\mathrm{Rn}}-T_{in}\right)G_{\mathrm{wn}-n+1}=\frac{{\lambda }_{\mathrm{wn}}}{0.5d_{\mathrm{wn}}}\left(T_{\mathrm{wn}}-T_{in}\right)G_{\mathrm{rn}-n+1}=0$ the internal temperature of the building T.sub.in and the temperature of the external street canyon are assumed in a quasi-steady equilibrium state, then, if it is assumed that the internal temperature T.sub.in of the building under air conditioning or natural ventilation is substantially constant in a tropical island climate, an average temperature in a center of the interior of the building T.sub.in is $T_{in}=\frac{{\left(h/b\right)}^2\left(T_{\mathrm{ww}}+T_{\mathrm{ww}}\right)+T_R}{2{\left(h/b\right)}^2+1}$ where b is an average width of the building.

8. The method according to claim 1, wherein the wind velocity comprises: in the street canyon, the wind velocity is decomposed into a vertical velocity W.sub.can along the wall and a horizontal velocity U.sub.can along the length of the street; the horizontal velocity along the width of the street is ignored; according to an observation, in a part close to the top of the street canyon, regardless of an air stability and a wind direction above the street canyon, a standard deviation σ.sub.w of a vertical wind velocity is equal to a friction velocity u.sub.*; the part σ.sub.w/u.sub.* close to the roof is 1.15, which is the same order of magnitude as an observed result; for an inertial boundary layer, a deviation of u.sub.* is not more than 10%; therefore, for any aspect ratio of the street canyon, the vertical velocity is assumed to be
W.sub.can=u.sub.*=√{square root over (C.sub.d)}┌U.sub.air┐ where U.sub.air is a wind velocity of the first layer of an atmospheric model, and C.sub.d is a drag coefficient, which is calculated from the temperature/humidity in and above the street canyon, a roughness Z.sub.0, and a stability effect; the horizontal wind velocity at the top of the street canyon U.sub.can is obtained by means of a Log approximate curve, a processing range of the Log curve being from h/3 of a lower part of the roof to a height of the first layer of the atmospheric model, wherein h is the height of the street canyon; when all street canyon orientations are considered, 360° integral processing is performed, then the velocity at the top of the street canyon U.sub.top is $U_{\mathrm{top}}=\frac{2}{\pi }\frac{\ln \left(\frac{h/3}{z_0}\right)}{\ln \left(\frac{\Delta z+h/3}{z_0}\right)}.Math.U_{\mathrm{air}}.Math.$ where Δz is a height from the roof to the first layer of the atmospheric model; the horizontal wind velocity U.sub.can is determined according to the wind velocity at ½ height of the street canyon; in order to calculate U.sub.can, a reasonable change law of U.sub.can in the vertical direction needs to be assumed; according to a continuity assumption of the wind velocity, a change curve of U.sub.can in the vertical direction has the following form
U.sub.can=U.sub.top exp(−N/2) where a value of N is slightly different; according to an aspect ratio of the street canyon (h/w=1-4), the value of U.sub.can varies from 0.75 U.sub.top to 0.4 U.sub.top; N=0.5(h/w), the horizontal wind velocity in the street canyon U.sub.can is $U_{\mathrm{can}}=\frac{2}{\pi }\exp \left(-0.25\frac{h}{w}\right)\frac{\ln \left(\frac{h/3}{z_0}\right)}{\ln \left(\frac{\Delta z+h/3}{z_0}\right)}.Math.U_{\mathrm{air}}.Math.$ calculations of aerodynamic roughness of the pavement and the wall in the street canyon are simplified, and the two are considered to have equal aerodynamic roughness, which is unrelated to the stability inside and outside the street canyon,
RES.sub.w=RES.sub.r=(11.8+4.2√{square root over (U.sub.can.sup.2+W.sub.can.sup.2)}).sup.−1 where the parameters RES.sub.w and RES.sub.r are inverses of C.sub.pC.sub.H1 and C.sub.pC.sub.H2, and are used for calculating sensible and latent heat flows.

9. The method according to claim 1, wherein the net convective heat flux is: under quasi-steady state conditions, the wind flow in an x-axis direction has been stable along the length of the street canyon; if an influence of pedestrians and vehicles inside the street is not considered, laws of mass conservation and momentum conservation are used in the x-axis direction, so as to obtain the horizontal movement of air inside the street canyon; if the air density and the horizontal velocity in the street canyon are processed as quasi-steady state variables, then in the case where an entrance velocity and an exit velocity satisfy outflow conditions, the laws of mass conservation and momentum conservation are written as follows according to a one-dimensional flow equation in the x-axis direction: $\{\begin{array}{c}\frac{\partial \left(p\stackrel{\_}{u}\right)}{\partial x}={\stackrel{.}{m}}_c\\ \frac{\partial \left(p{\stackrel{\_}{u}}^2\right)}{\partial x}+\frac{\partial p}{\partial x}+\frac{\partial {\tau }_w}{\partial x}=0\\ \stackrel{\_}{u}{.Math.}_{x=x_0}={\stackrel{\_}{u}}_{in}\\ \stackrel{\_}{u}{.Math.}_{x=x_0+L}={\stackrel{\_}{u}}_{\mathrm{out}}\\ \frac{d\stackrel{\_}{u}}{\mathrm{dx}}{.Math.}_{x=x_0+L}=0\end{array}$ where ρ is the air density; ū is an average velocity in the x-axis direction; {dot over (m)}.sub.c is a specific volumetric mass flow of air entering or exiting a control body, a specific volume being a ratio of the mass flow of the air entering or exiting the control body to the volume of the control body; p is an average pressure of a cross section of the street canyon; τ.sub.w is an average frictional stress of the wall surface and the street surface to the air flow; ū.sub.in and ū.sub.out are average air velocities at the entrance and exit of the street canyon, respectively; x.sub.0 is an entrance position of the street canyon, and the flow velocity at the entrance of the street canyon is measured by an instrument; from a perspective of regional scale, streets are usually connected to form a street network; by studying a road network, a horizontal flux of the urban heat island phenomenon can be replaced by a street canyon formed by only one independent street; when the horizontal flux of a crossroad is calculated, a Kirchhoff's principle for calculating a fluid network can be used; according to a topological structure and plan theory of the street network, a street network can be represented by a corresponding adjacency matrix, and the horizontal air flow of each branch can be solved; for a scenario at a crossroad, it is assumed that the horizontal flux at the exit of street m is Q.sub.m,out and the horizontal flux at the exit of street j is Q.sub.j,out, the horizontal flux at the entrance of street i is Q.sub.i,out and the horizontal flux at the entrance of street n is Q.sub.n,out; then, according to a law of energy conservation, a horizontal flux at a node is:
(ρ{dot over (V)}.sub.mQ.sub.m,out)+(ρ{dot over (V)}.sub.jQ.sub.j,out)=ρQ.sub.min({dot over (V)}.sub.n+{dot over (V)}.sub.i) where {dot over (V)} is an air volume flow of each street canyon; Q.sub.mix is a mixed horizontal flux; the mixed horizontal flux Q.sub.mix is the horizontal flux flowing into the entrances of street canyons n and i at the node, $Q_{\mathrm{mix}}=\frac{\rho {.Math.}_{\mathrm{All}\mathrm{outflowing}\mathrm{branch}}{\left(\stackrel{.}{V}\right)}^k}{{.Math.}_{\mathrm{All}\mathrm{inflowing}\mathrm{branch}}{\left(\rho \stackrel{.}{V}{Q}_{\mathrm{out}}\right)}^k}.$

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0082] In order to illustrate the technical solutions of the embodiments of the present invention more clearly, the following briefly introduces the accompanying drawings used in the description of the embodiments. Apparently, the drawings in the following description are merely some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained from these drawings without any creative effort.

[0083] FIG. 1 is a cross-sectional view of a street canyon canopy model according to an embodiment of the present invention.

[0084] FIG. 2 is an air flow analysis diagram along the length of a street according to an embodiment of the present invention.

[0085] FIG. 3 shows horizontal flux of a street canyon network at a node according to an embodiment of the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0086] Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the drawings show the exemplary embodiments of the present disclosure, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments illustrated herein. On the contrary, these embodiments are provided so that the present disclosure can be understood more thoroughly, and can fully convey the scope of the present disclosure to those skilled in the art.

[0087] The terms “include” and “have” and any variations thereof in the description embodiments, claims and drawings of the present invention are intended to cover non-exclusive inclusion, for example, including a series of steps or elements.

[0088] The technical solution of the present invention will be described in further detail below with reference to the accompanying drawings and embodiments.

[0089] Main assumptions of a tropical island urban canopy model are as follows:

[0090] The space formed by a road and opposite sides of buildings on both sides is referred to as a street canyon.

[0091] The length of the road is much greater than its width. The aspect ratio of a street changes with the actual situation of the road.

[0092] The angle between the direction of the road and the sun exposure is arbitrary.

[0093] The buildings on both sides of the road have the same height and width, and the roof height is the surface layer height of an atmospheric model.

[0094] The air flow and energy exchange in the street canyon change slowly, with hourly time resolution. The thermal and dynamic processes of the heat island effect can be regarded as quasi-steady state processes.

[0095] Feedbacks of urban heat islands on the external atmospheric effect are ignored.

[0096] A 2-layer model is used to describe the heat conduction among a pavement, a wall and a roof.

[0097] As shown in FIG. 1, a method for building an urban canopy model based on tropical island climate characteristics includes:

[0098] linking adjacent regions together, multiple streets of finite lengths within the regions affecting each other;

[0099] The energy balance equation in a street is as follows:

Q*.sub.s+Q.sub.F,s=Q.sub.H,s+Q.sub.E,s+ΔQ.sub.S,s+ΔQ.sub.A,s

[0100] where Q* is a net radiation heat flux, in unit of W/m.sup.2; Q.sub.F is anthropogenic heat production, in unit of W/m.sup.2; Q.sub.H is a sensible heat flux, in unit of W/m.sup.2; Q.sub.E is a latent heat flux, in unit of W/m.sup.2; ΔQ.sub.S is a net heat storage flux, in unit of W/m.sup.2; ΔQ.sub.A is a net convective heat flux, in unit of W/m.sup.2; s is the street;

[0101] Net radiation heat flux Q*=long-wave radiation+short-wave radiation.

[0102] The long-wave radiation includes:

[0103] net long-wave radiation of a pavement L.sub.r* is:

L.sub.r*=ϵ.sub.rΨ.sub.rL.sup.↓−ϵ.sub.rσT.sub.r.sup.4+ϵ.sub.rϵ.sub.w(1−Ψ.sub.r)σT.sub.w.sup.4+ϵ.sub.r(1−ϵ.sub.w)(1−Ψ.sub.r)Ψ.sub.wL.sup.↓+ϵ.sub.rϵ.sub.w(1−ϵ.sub.w)(1−‥.sub.r)(1−2Ψ.sub.w)σT.sub.w.sup.4+ϵ.sub.r(1−ϵ.sub.w(1−Ψ.sub.r)Ψ.sub.wσϵ.sub.rT.sub.r.sup.4

[0104] and net long-wave radiation of a wall L.sub.w* is:

L.sub.w*=ϵ.sub.wΨ.sub.wL.sup.↓−ϵ.sub.wσT.sub.w.sup.4+ϵ.sub.wΨ.sub.wσϵ.sub.rT.sub.r.sup.4+ϵ.sub.w.sup.2(1−2Ψ.sub.w)σT.sub.w.sup.4+ϵ.sub.w(1−ϵ.sub.r)Ψ.sub.rΨ.sub.wL.sup.↓+ϵ.sub.w(1−ϵ.sub.w)Ψ.sub.w(1−2Ψ.sub.w)L.sup.↓+ϵ.sub.w.sup.2(1−ϵ.sub.w)(1−2Ψ.sub.w).sup.2σT.sub.w.sup.4+ϵ.sub.w.sup.2(1−ϵ.sub.r)Ψ.sub.w(1−Ψ.sub.r)σT.sub.w.sup.4+ϵ.sub.w(1−ϵ.sub.w)Ψ.sub.w(1−2Ψ.sub.w)σϵ.sub.rT.sub.r.sup.4

[0105] where L↓ is the amount of solar radiation, σ is a standard deviation, T.sub.r and T.sub.w are temperatures of the pavement and the wall, and ϵ.sub.r and ϵ.sub.w are emissivities of the pavement and the wall; for the pavement, if Ψ.sub.r is a sky viewing angle coefficient of the pavement to the sky, (1−Ψ.sub.r) is a sky viewing angle coefficient of the pavement to the walls on both sides; a sky viewing angle coefficient of the wall to the sky is Ψ.sub.w, a sky viewing angle coefficient to the pavement is Ψ.sub.w, then a sky viewing angle coefficient to the opposite wall is (1−2Ψ.sub.w), and the sky viewing angle coefficient is 1.0 for a roof;

[0106] The sky viewing angle coefficient is calculated using a plane angle, the sky viewing angle coefficient at the w/2 position of the pavement is

[00023]${\Psi }_r\left(y=w/2\right)=1+\frac{1}{\pi }{\tan }^{-1}\left(\frac{h/w}{{\left(h/w\right)}^2-1/4}\right)\mathrm{for}h/w\le 1/2{\Psi }_r\left(y=w/2\right)=1-\frac{1}{\pi }{\tan }^{-1}\left(\frac{h/w}{{\left(h/w\right)}^2-1/4}\right)\mathrm{for}h/w\ge 1/2$

[0107] The sky viewing angle coefficient at the intersection of the wall and the pavement is

[00024]${\Psi }_w\left(z=0\right)=\frac{1}{\pi }{\tan }^{-1}\left(\frac{1}{h/w}\right),$

where h represents the height of a street canyon, and w represents the width of the street canyon.

[0108] The short-wave radiation includes:

[0109] Average direct solar radiant fluxes of the pavement , west wall , east wall and roof are calculated according to the perpendicular angle of the street to the sun direction:

[00025]$S_r^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\{\begin{array}{c}0\mathrm{for}.Math.\lambda .Math.>{\lambda }_0\\ \left[1-\frac{h}{w}\tan \left(\lambda \right)\right]S^{\Downarrow }\mathrm{for}.Math.\lambda .Math.<{\lambda }_0\end{array}{S}_{\mathrm{ww}}^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\{\begin{array}{c}\frac{1}{2}\frac{w}{h}{S}^{\Downarrow }\mathrm{for}-\frac{\pi }{2}\le \lambda \le -{\lambda }_0\\ \frac{1}{2}\tan \left(\lambda \right)S^{\Downarrow }\mathrm{for}-{\lambda }_0\le \lambda \le 0\\ 0\mathrm{for}0\le \lambda \le \frac{\pi }{2}\end{array}{S}_{\mathrm{we}}^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\{\begin{array}{c}0\mathrm{for}-\frac{\pi }{2}\le \lambda \le -{\lambda }_0\\ \frac{1}{2}\tan \left(\lambda \right)S^{\Downarrow }\mathrm{for}-{\lambda }_0\le \lambda \le 0\\ \frac{1}{2}\frac{w}{h}{S}^{\Downarrow }\mathrm{for}0\le \lambda \le \frac{\pi }{2}\end{array}{S}_R^{\Downarrow }\left(\theta =\frac{\pi }{2}\right)=\chi S^{\Downarrow }$

[0110] where is direct solar radiation on a horizontal surface, θ is an angle between a sun angle and an axial direction of the street canyon, λ is an angle between a sun direction and a normal direction of the wall, χ is a ratio of a direct radiation to a total radiation at a top of the street canyon, h is a height of the street canyon, and w is a width of the street canyon;

[0111] According to the orientation change of the street canyon, the width w of the street canyon is corrected as w/sin θ; after a heat flux of the wall is obtained, the heat flux of the wall is multiplied by sin θ for correction, where θ.sub.0 is an orientation of the street canyon where the pavement does not receive direct sunlight at all

[00026]${\theta }_0=\arcsin \left(\min \left(\frac{w}{h}\frac{1}{\tan \lambda };1\right)\right)$

[0112] All direct radiant fluxes obtained by the street canyon are averaged according to all possible changes in the direction of the street canyon, wherein two integrals are used, one between θ=0 and θ=θ.sub.0, and the other between θ=θ.sub.0 and

[00027]$\theta =\frac{\pi }{2};$

[0113] Wherein average direct solar fluxes of the wall, the pavement and the roof are:

[00028]$S_r^{\Downarrow }=S^{\Downarrow }\left[\frac{2{\theta }_0}{\pi }-\frac{2}{\pi }\frac{h}{w}\tan \left(\lambda \right)\left(1-\cos \left({\theta }_0\right)\right)\right]S_w^{\Downarrow }=S^{\Downarrow }\left[\frac{w}{h}\left(\frac{1}{2}-\frac{{\theta }_0}{\pi }\right)+\frac{1}{\pi }\tan \left(\lambda \right)\left(1-\cos \left({\theta }_0\right)\right)\right]S_R^{\Downarrow }=S^{\Downarrow }$

[0114] S.sup.↓ is a scattered solar radiation available on the horizontal surface, and an amount of scattered solar radiation received by a surface in the street canyon is directly obtained from the sky viewing angle coefficient. Due to an influence of a shape of the street canyon and high-reflectivity building surface materials, the short-wave radiation is calculated to solve a geometric system with an infinite number of reflecting surfaces, a reflecting processes of which are assumed to be isentropic processes;

[0115] When direct and diffuse reflectivities of each surface are the same, an energy stored by the pavement A.sub.r and an energy stored by the wall A.sub.w when a first reflection occurs are:

A.sub.r(0)=(1−α.sub.r)(+S.sub.r.sup.↓)

A.sub.w(0)=(1−α.sub.w)(+S.sub.w.sup.↓)

[0116] where α.sub.r and α.sub.w represent the reflectivities of the pavement and the wall, respectively;

[0117] The energies of the reflected parts from the pavement R.sub.r and from the wall R.sub.w are:

R.sub.r(0)=α.sub.r(+S.sub.r.sup.↓)

R.sub.w(0)=α.sub.w(+S.sub.w.sup.↓)

[0118] After n reflections occur,

A.sub.r(n+1)=A.sub.r(n)+(1−α.sub.r)(1−Ψ.sub.r)R.sub.w(n)

A.sub.w(n+1)=A.sub.w(n)+(1−α.sub.r)Ψ.sub.wR.sub.r(n)+(1−α.sub.w)(1−2Ψ.sub.w)R.sub.w(n)

R.sub.r(n+1)=α.sub.r(1−Ψ.sub.r)R.sub.w(n)

R.sub.w(n+1)=α.sub.wΨ.sub.wR.sub.r(n)+α.sub.w(1−2Ψ.sub.w)R.sub.w(n)

[0119] The following is obtained according to recursive formulas,

[00029]$A_r\left(n+1\right)=A_r\left(0\right)+\left(1-{\alpha }_r\right)\left(1-{\Psi }_r\right)\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)$$A_w\left(n+1\right)=A_w\left(0\right)+{\Psi }_w\left(1-{\alpha }_r\right)\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)+\left(1-2{\Psi }_w\right)\left(1-{\alpha }_w\right)\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)$$\mathrm{and}\underset{k=0}{\overset{n}{.Math.}}{R}_r\left(k\right)=\left(1-{\Psi }_r\right){\alpha }_r\underset{k=0}{\overset{n-1}{.Math.}}{R}_w\left(k\right)+R_r\left(0\right)\underset{k=0}{\overset{n}{.Math.}}{R}_w\left(k\right)={\alpha }_w{\Psi }_w\underset{k=0}{\overset{n-1}{.Math.}}{R}_r\left(k\right)+{\alpha }_w\left(1-2{\Psi }_w\right)\underset{k=0}{\overset{n-1}{.Math.}}{R}_w\left(k\right)+R_w\left(0\right)$

[0120] For the infinite reflections, the following is obtained by solving a geometric system,

[00030]$\underset{k=0}{\overset{\infty }{.Math.}}{R}_r\left(k\right)=\frac{R_r\left(0\right)+\left(1-{\Psi }_r\right){\alpha }_r\left(R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)\right)}{1-\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}\underset{k=0}{\overset{\infty }{.Math.}}{R}_r\left(k\right)=\frac{R_r\left(0\right)+\left(1-{\Psi }_r\right){\alpha }_r\left(R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)\right)}{1-\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}$

[0121] M.sub.r is assumed to be a sum of a pavement reflection and M.sub.w is assumed to be a sum of a wall reflection,

[00031]$M_r=\frac{R_r\left(0\right)+\left(1-{\Psi }_r\right){\alpha }_r\left(R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)\right)}{1-\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}{M}_w=\frac{R_w\left(0\right)+{\Psi }_w{\alpha }_w{R}_r\left(0\right)}{1-\left(1-2{\Psi }_w\right){\alpha }_w+\left(1-{\Psi }_r\right){\Psi }_w{\alpha }_r{\alpha }_w}$

[0122] where

R.sub.r(0)=α.sub.r+α.sub.rS.sub.r.sup.↓

R.sub.w(0)=α.sub.r+α.sub.rS.sub.w.sup.↓

[0123] A total solar radiation absorbed by the pavement is S.sub.r*, a total solar radiation absorbed by the wall is S.sub.w* , a total solar radiation absorbed by the roof is S.sub.R*:

S.sub.r*=(1−α.sub.r)+(1−α.sub.r)S.sub.r.sup.↓+(1−α.sub.r)(1−Ψ.sub.r)M.sub.w

S.sub.w*=(1−α.sub.w)+(1−α.sub.w)S.sub.w.sup.↓+(1−α.sub.w)(1−2Ψ.sub.w)M.sub.w+(1−α.sub.w)Ψ.sub.wM.sub.r

S.sub.R*=(1−α.sub.R)+(1−α.sub.R)S.sub.R.sup.↓.

[0124] The anthropogenic heat production specifically includes:

[0125] current anthropogenic heat flux in the street canyon is Q.sub.F=Q.sub.FV+Q.sub.FH+Q.sub.FM;

[0126] where Q.sub.FV, Q.sub.FH and Q.sub.FM are heat generated by vehicles, fixed heat sources and biological metabolism, respectively.

[0127] The sensible heat flux Q.sub.H includes:

Q.sub.H,r,ww,we=ρC.sub.pC.sub.H1U.sub.can(T.sub.r,ww,we−T.sub.can)

Q.sub.H,R=ηC.sub.pC.sub.H2U.sub.top(T.sub.R−T.sub.air)

Q.sub.H,can=ρC.sub.pC.sub.H2U.sub.air(T.sub.can−T.sub.air)

[0128] where r, ww, we, R, and can refer to the pavement, the west wall, the east wall, the roof, and the street canyon, respectively; ρ is an air density; C.sub.p is a specific heat under a constant pressure; T.sub.can is a temperature in a center of the street canyon (w/2, h/2); U.sub.can and U.sub.top are a wind velocity in the center of the street canyon (w/2, h/2) and a wind velocity above the street canyon;

[0129] U.sub.air and T.sub.air are an input wind velocity and an input temperature at a reference height of a turbulence model, and C.sub.H1 and C.sub.H2 are dimensionless velocity transfer coefficients; differences between the C.sub.H1 and C.sub.H2 are only a height and a roughness of a reference layer; the same zero plane layer and roughness are used, and the values of the two are equal, and are calculated as follows:

[00032]$C_{H*}{u}_{*}=\frac{{\mathrm{ku}}_{*}}{{\Psi }_h}$

[0130] where k is a Von Karman constant, u.sub.* is a friction velocity of the reference layer, and Ψ.sub.h is a general integral function,

[00033]${\Psi }_h={\int }_{{\zeta }^T}^{{\zeta }^{\prime }}\frac{{\varphi }_h}{\zeta }d\zeta$

[0131] where ζ′=(z.sub.a−d)/L; ζ.sup.T=z.sub.T/L, z.sub.T is a roughness length of a heat flow; L is an Obukhov stability length,

[00034]$L=-\frac{\rho C_p{\mathrm{Tu}}_{*}^3}{{\mathrm{kgH}}_a}$

[0132] where T is an average temperature of this layer, H.sub.a is an air flux between the street canyon and an atmosphere, and L is an implicit function, which is solved by simplified iteration. When a specific heat of air in an urban canopy is ignored, H.sub.a is a weighted average of a wall flux and a road flux in the street canyon, that is,

H.sub.a=2(h/w)Q.sub.w+Q.sub.R

[0133] In a TEB model proposed by Masson, C.sub.H*u.sub.* is a reciprocal of aerodynamic resistance, i.e. 1/RES.sub.*, which is determined by the wind velocities in the street canyon and at a top of the street canyon;

[0134] If a surface covered by plants such as green space is not considered, an average sensible heat flow of the street canyon depends on a weighted average area of the roof, the wall and the pavement,

[00035]$Q_H=\frac{{\mathrm{bQ}}_{H,R}+{\mathrm{wQ}}_{H,r}+h\left(Q_{H,\mathrm{ww}}+Q_{H,\mathrm{we}}\right)}{w+b},$

wherein h is a height of the street canyon, w is a width of the street canyon, and b is an average width of the buildings.

[0135] Optionally, the latent heat flux Q.sub.E includes:

[0136] a direct latent heat flow between a building roof and the atmosphere

Q.sub.E,R=l.sub.vB.sub.RρC.sub.H2U.sub.top(q.sub.R−q.sub.air)

[0137] where l.sub.v is a latent heat of evaporation, B.sub.R is a humidity parameter of the roof, between 0 and 1, 0 is completely dry, 1 is completely wet, a value of B depends on the plant and water conditions of the surface, ρ is a density of the air, C.sub.H2 is a dimensionless velocity transfer coefficient, and q.sub.R is a humidity of the roof; q.sub.air is a humidity at the reference height,

[0138] A latent heat flow is calculated using a similarity law for the air on the pavement and the wall and in the street canyon

Q.sub.E,r=l.sub.vB.sub.rρC.sub.H1U.sub.can(q.sub.r−q.sub.can)

Q.sub.E,w=0

[0139] wherein C.sub.H2 is a dimensionless velocity transfer coefficient, and q.sub.r is a humidity of the pavement; a q.sub.can is a humidity at the street canyon,

[0140] A latent heat flow between an interior of the street canyon and a top atmosphere is

Q.sub.E,can=l.sub.vρC.sub.H2U.sub.air(q.sub.can−q.sub.air)

[0141] The net heat storage flux ΔQ.sub.S includes:

[0142] Because there is a temperature gradient inside the building or the pavement, the roof, the wall and the pavement are assumed to be of at least three-layer structures. For an outermost layer structure, heat transfer equations of the three planes are written as,

[00036]$C_{R1}\frac{\partial T_{R1}}{\partial t}=\frac{{\lambda }_{R1}}{d_{R1}}\left(S_R^{*}+L_R^{*}-H_R-{\mathrm{LE}}_R-G_{R1-2}\right)C_{w1}\frac{\partial T_{w1}}{\partial t}=\frac{{\lambda }_{w1}}{d_{w1}}\left(S_w^{*}+L_w^{*}-H_w-G_{w1-2}\right)C_{r1}\frac{\partial T_{r1}}{\partial t}=\frac{{\lambda }_{r1}}{d_{r1}}\left(S_r^{*}+L_r^{*}-H_r-{\mathrm{LE}}_r-G_{r1-2}\right)$

[0143] where T.sub.*i is temperature of the i-th layer; C.sub.*i is a specific heat capacity of the air; d.sub.*i is a layer thickness, fluxes S.sub.**, L.sub.**, H.sub.*, LE.sub.*, and G.sub.*1-2 are net solar radiation, net infrared radiation, sensible heat, latent heat, and thermal conductivity between the surface layer and the lower layer, and the thermal conductivity is calculated using a Fourier heat conduction equation,

[00037]$G_{*1-2}=\frac{{\lambda }_{*1-2}}{0\text{.5}\left(d_{*1}+d_{*2}\right)}\left(T_{*1}-T_{*2}\right)$

[0144] An average thermal conductivity between two adjacent layers λ.sub.*1-2 is calculated using a geometric average method:

[00038]${\lambda }_{*1-2}=\frac{d_{*1}+d_{*2}}{\left(d_{*1}/{\lambda }_{*1}\right)+\left(d_{*2}/{\lambda }_{*2}\right)}$

[0145] where λ.sub.*i is a thermal conductivity of the i-th layer;

[0146] An inner first layer of the surface is assumed to be a very thin surface, and a temperature of the first layer is simplified to an outer surface temperature; for other inner i-th layers, a thermal conductivity between adjacent layers is calculated; for an innermost layer, such as the n-th layer, an internal temperature of the building issued for the roof and the wall surface, and the 0 flux is used for the pavement;

[00039]$G_{Rn-n+1}=\frac{{\lambda }_{Rn}}{0.5d_{Rn}}\left(T_{Rn}-T_{in}\right)$$G_{wn-n+1}=\frac{{\lambda }_{wn}}{0.5d_{wn}}\left(T_{wn}-T_{in}\right)$$G_{rn-n+1}=0$

[0147] The internal temperature of the building T.sub.in and the temperature of the external street canyon are assumed in a quasi-steady equilibrium state, then, if it is assumed that the internal temperature T.sub.in of the building under air conditioning or natural ventilation is substantially constant in a tropical island climate, an average temperature in a center of the interior of the building T.sub.in is

[00040]$T_{in}=\frac{{\left(h/b\right)}^2\left(T_{ww}+T_{ww}\right)+T_R}{2{\left(h/b\right)}^2+1}$

[0148] where b is an average width of the building.

[0149] The advantage of simplifying the internal temperature of the building is that there is no need to assume a source term existing inside the building in the presence of a space heating or cooling system; moreover, the power consumption of the heating or cooling system inside the building is difficult to estimate. In this way, the heat flux storage inside the building can be uniformly treated with temperature boundary conditions.

[0150] The wind velocity includes: as shown in FIG. 2, in the street canyon, the wind velocity is decomposed into a vertical velocity W.sub.can along the wall and a horizontal velocity U.sub.can along the length of the street; the horizontal velocity along the width of the street is ignored;

[0151] According to an observation, in a part close to the top of the street canyon, regardless of an air stability and a wind direction above the street canyon, a standard deviation σ.sub.w of a vertical wind velocity is equal to a friction velocity u.sub.*;

[0152] The part σ.sub.w/u.sub.* close to the roof is 1.15, which is the same order of magnitude as an observed result. For an inertial boundary layer, a deviation of u.sub.* is not more than 10%. Therefore, for any aspect ratio of the street canyon, the vertical velocity is assumed to be

W.sub.can=u.sub.*=√{square root over (C.sub.d)}┌U.sub.air┐

[0153] where U.sub.air is a wind velocity of the first layer of an atmospheric model, and C.sub.d is a drag coefficient, which is calculated from the temperature/humidity in and above the street canyon, a roughness Z.sub.0, and a stability effect;

[0154] The horizontal wind velocity at the top of the street canyon U.sub.can is obtained by means of a Log approximate curve, a processing range of the Log curve is from h/3 of a lower part of the roof to a height of the first layer of the atmospheric model, wherein h is the height of the street canyon. When all street canyon orientations are considered, 360° integral processing is performed, then the velocity at the top of the street canyon U.sub.top is

[00041]$U_{\mathrm{top}}=\frac{2}{\pi }\frac{\ln \left(\frac{h/3}{z_0}\right)}{\ln \left(\frac{\Delta z+h/3}{z_0}\right)}.Math.U_{air}.Math.$

[0155] where Δz is a height from the roof to the first layer of the atmospheric model;

[0156] The horizontal wind velocity U.sub.can is determined according to the wind velocity at ½ height of the street canyon;

[0157] In order to calculate U.sub.can, a reasonable change law of U.sub.can in the vertical direction needs to be assumed;

[0158] According to a continuity assumption of the wind velocity, a change curve of U.sub.can in the vertical direction has the following form

U.sub.can=U.sub.top exp(−N/2)

[0159] where a value of N is slightly different;

[0160] According to an aspect ratio of the street canyon (h/w=1-4), the value of U.sub.can varies from 0.75 U.sub.top to 0.4 U.sub.top;

[0161] N=0.5(h/w), the horizontal wind velocity in the street canyon U.sub.can is

[00042]$U_{can}=\frac{2}{\pi }\exp \left(-0.25\frac{h}{w}\right)\frac{\ln \left(\frac{h/3}{z_0}\right)}{\ln \left(\frac{\Delta z+h/3}{z_0}\right)}.Math.U_{air}.Math.$

[0162] Calculations of aerodynamic roughness of the pavement and the wall in the street canyon are simplified, and the two are considered to have equal aerodynamic roughness, which is unrelated to the stability inside and outside the street canyon,

RES.sub.w=RES.sub.r=(11.8+4.2√{square root over (U.sub.can.sup.2+W.sub.can.sup.2)}).sup.−1

[0163] where the parameters RES.sub.w and RES.sub.r are inverses of C.sub.pC.sub.H1 and C.sub.pC.sub.H2, and are used for calculating sensible and latent heat flows.

[0164] The net convective heat flux is:

[0165] Under quasi-steady state conditions, the wind flow in an x-axis direction has been stable along the length of the street canyon; if an influence of pedestrians and vehicles inside the street is not considered, laws of mass conservation and momentum conservation are used in the x-axis direction, so as to obtain the horizontal movement of air inside the street canyon;

[0166] If the air density and the horizontal velocity in the street canyon are processed as quasi-steady state variables, then in the case where an entrance velocity and an exit velocity satisfy outflow conditions, the laws of mass conservation and momentum conservation are written as follows according to a one-dimensional flow equation in the x-axis direction:

[00043]$\{\begin{array}{c}\frac{\partial \left(\rho \overline{u}\right)}{\partial x}={\dot{m}}_c\\ \frac{\partial \left(\rho {\overline{u}}^2\right)}{\partial x}+\frac{\partial p}{\partial x}+\frac{\partial {\tau }_w}{\partial x}=0\\ {\stackrel{\_}{u}.Math.}_{x=x_0}={\overline{u}}_{in}\\ {\stackrel{\_}{u}.Math.}_{x=x_0+L}={\overline{u}}_{\mathrm{out}}\\ {\frac{d\stackrel{\_}{u}}{\mathrm{dx}}.Math.}_{x=x_0+L}\end{array}$

[0167] where ρ is the air density; ū is an average velocity in the x-axis direction; {dot over (m)}.sub.c is a specific volumetric mass flow of air entering or exiting a control body, a specific volume being a ratio of the mass flow of the air entering or exiting the control body to the volume of the control body; p is an average pressure of a cross section of the street canyon; τ.sub.w is an average frictional stress of the wall surface and the street surface to the air flow; ū.sub.in and ū.sub.out are average air velocities at the entrance and exit of the street canyon, respectively; x.sub.0 is an entrance position of the street canyon, and the flow velocity at the entrance of the street canyon is measured by an instrument;

[0168] From a perspective of regional scale, streets are usually connected to form a street network; by studying a road network, a horizontal flux of the urban heat island phenomenon can be replaced by a street canyon formed by only one independent street; when the horizontal flux of a crossroad is calculated, a Kirchhoff's principle for calculating a fluid network can be used; according to a topological structure and plan theory of the street network, a street network can be represented by a corresponding adjacency matrix, and the horizontal air flow of each branch can be solved.

[0169] As shown in FIG. 3, for a scenario at a crossroad, it is assumed that the horizontal flux at the exit of street m is Q.sub.m,out and the horizontal flux at the exit of street/is Q.sub.j,out, the horizontal flux at the entrance of street i is Q.sub.i,out and the horizontal flux at the entrance of street n is Q.sub.n,out; then, according to aa law of energy conservation, a horizontal flux at a node is:

(ρ{dot over (V)}.sub.mQ.sub.m,out)+(ρ{dot over (V)}.sub.jQ.sub.j,out)=ρQ.sub.min({dot over (V)}.sub.n+{dot over (V)}.sub.i)

[0170] where {dot over (V)} is an air volume flow of each street canyon; Q.sub.mix is a mixed horizontal flux;

[0171] The mixed horizontal flux Q.sub.mix is the horizontal flux flowing into the entrances of street canyons n and i at the node,

[00044]$Q_{mix}=\frac{\rho {.Math.}_{\mathrm{All}\mathrm{outflowing}\mathrm{branch}}{\left(\stackrel{.}{V}\right)}^k}{{.Math.}_{\mathrm{All}\mathrm{inflowing}\mathrm{branch}}{\left(\rho \stackrel{.}{V}{Q}_{out}\right)}^k}.$

[0172] Beneficial Effects:

[0173] With the continuous rapid development of cities, the urban heat island effect is increasingly serious. The remission of urban heat islands can help suppress the spread of infectious diseases, reduce greenhouse gas emissions, and reduce building energy consumption. In order to study the heat island effect, the urban canopy model is required to model a region. The improvement on the adaptability of the model is the basis for improving the accuracy of the model for the study on the heat island effect.

[0174] The new tropical island urban canopy model proposed here considers the influence of horizontal convective flux on the heat island effect, which makes up for the deficiency of the classical urban canopy model in the study on the heat island effect of “tropical island” cities, and has higher adaptability to the “tropical island” cities. In addition, the new model, i.e. an urban canopy model of a three-dimensional space, is built based on the finite length of a street canopy. This improvement can greatly improve the calculation accuracy of the horizontal heat flux in the street canyon.

[0175] The above specific embodiments further describe the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above are only specific embodiments of the present invention, and are not intended to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention should be included within the protection scope of the present invention.