劉大龍，楊 柳，霍旭杰，劉加平

(1.西安建筑科技大學 建筑學院，西安 710055；2.陜西省西部綠色建筑協(xié)同創(chuàng)新中心，西安 710055)

?

建筑節(jié)能分析太陽總輻射模型研究綜述

劉大龍1,2，楊 柳1,2，霍旭杰1,2，劉加平1,2

(1.西安建筑科技大學 建筑學院，西安 710055；2.陜西省西部綠色建筑協(xié)同創(chuàng)新中心，西安 710055)

太陽輻射是建筑節(jié)能分析的重要基礎(chǔ)氣象參數(shù)，實測數(shù)據(jù)遠遠不能滿足需求，理論計算是目前獲取輻射數(shù)據(jù)的主要途徑。將常用水平面太陽總輻射模型歸納為氣象參數(shù)、空間插值和基于DEM三類，詳述了各自的原理和計算方法。對三類模型在建筑節(jié)能分析中的適用性進行了分析，展望了建筑節(jié)能分析用太陽輻射模型的發(fā)展趨勢：氣象參數(shù)模型與DEM模型的融合。

水平面太陽總輻射；氣象參數(shù)模型；空間插值模型；DEM輻射模型

太陽輻射數(shù)據(jù)是建筑節(jié)能設(shè)計、暖通空調(diào)設(shè)計、建筑采光設(shè)計、城市規(guī)劃和景觀設(shè)計等領(lǐng)域的重要氣象參數(shù)。實測數(shù)據(jù)是獲取太陽輻射最準確的途徑，但中國對太陽輻射數(shù)據(jù)的測量臺站少，測量歷史短，測量項目少。中國幅員遼闊，地形復(fù)雜，下墊面條件多樣，使得各地輻射條件相差懸殊，導(dǎo)致各地的太陽輻射強度相差很大。現(xiàn)有的測量太陽輻射數(shù)據(jù)遠遠不能滿足工程和科研需要。

太陽輻射模型是根據(jù)相關(guān)原理構(gòu)建數(shù)學模型來獲取太陽輻射數(shù)據(jù)的理論計算方法。該方法不受測量條件和時空限制是彌補太陽輻射數(shù)據(jù)不足的有效途徑[1]。太陽輻射模型包括水平面總輻射模型、直散分離模型、逐時輻射模型等，主要對水平面日總輻射模型進行歸納整理。

1 水平面太陽總輻射計算模型

水平面太陽總輻射模型種類眾多，根據(jù)與太陽輻射的關(guān)聯(lián)方式和計算原理的不同，將其歸納為氣象參數(shù)模型、空間插值模型和基于DEM的輻射模型三類，下面對其分別進行論述。

1.1 氣象參數(shù)模型

1.1.1 計算原理 組成氣候系統(tǒng)的各氣象要素之間相互關(guān)聯(lián)，太陽輻射是氣候的主動因素，它對其他氣象參數(shù)產(chǎn)生影響，而這些參數(shù)反過來也反映了太陽輻射的特征。因此，選擇與太陽輻射關(guān)聯(lián)密切且便于測試的氣象要素，構(gòu)建其與太陽輻射之間的函數(shù)關(guān)系，就可計算出太陽輻射值，這是氣象參數(shù)輻射模型的計算原理。用于構(gòu)建太陽輻射模型的主要氣象參數(shù)有日照時數(shù)、溫差、云量。此外，使用相對濕度、降雨量、露點溫度等要素也能夠建立水平面太陽總輻射模型，但這類氣象要素與太陽輻射的關(guān)系較弱，不能單獨完成輻射的計算，必須與前3個要素中的一個或多個共同構(gòu)建總輻射模型。云量與輻射具有重要相關(guān)性，但云量的單獨輻射模型很少，且應(yīng)用也少，多數(shù)情況是和其他氣象參數(shù)共同構(gòu)成總輻射模型，因此，文中沒有列出云量的單獨輻射模型，詳細內(nèi)容可參考文獻[2-3]。

1.1.2 日照時數(shù)模型 日照時數(shù)模型是氣象參數(shù)模型，是所有水平面太陽總輻射模型中使用最廣泛、計算結(jié)果最準確且計算參數(shù)最容易獲得的一類模型。在該類模型中，很多情況以日照百分率(S/S0)為參數(shù)進行計算，日照百分率是實際日照時數(shù)與日最大日照時數(shù)的比值。最早提出日照時數(shù)模型的是?ngstr?m[4]，該模型如式(1)所示，直觀簡潔地給出了月均日總輻射量與晴天日總輻射量的比值同日照百分比之間的線性關(guān)系。孫治安等[5]指出?ngstr?m 模型在晴天條件下太陽總輻射計算值的誤差最小。

(1)

式中：G為月均日總輻射量；Gc為月均日晴天總輻射量；S為測量的月均日日照時數(shù)；S0為月均日最大可能日照時數(shù)；a、b為回歸系數(shù)。

系數(shù)a和b是使用?ngstr?m 模型的關(guān)鍵問題。可在已知輻射和日照時數(shù)的情況下通過回歸獲得系數(shù)a和b，然后將其用于氣候相近地區(qū)計算當?shù)氐奈粗椛渲怠２煌赜虻南禂?shù)a、b不同。高國棟等[6]以該模型計算了我國不同地區(qū)的a、b值，指出兩系數(shù)的分布與地理條件和氣候狀況有密切。兩個系數(shù)反映了輻射與日照率關(guān)系模型具有較強的地域性。系數(shù)a、b不僅具有地域特性，而且還具有季節(jié)特性，Soler[7]根據(jù)歐洲100個氣象站的輻射數(shù)據(jù)，通過回歸給出了每個月不同的系數(shù)a、b，見表1。鞠曉慧等[8]根據(jù)我國建站30年以上的輻射資料研究也表明，需按不同月份確定系數(shù)a、b。

表1 Soler模型中各月份的系數(shù)a、b

有學者發(fā)現(xiàn)了?ngstr?m 模型中系數(shù)a、b與日照百分率之間存在函數(shù)關(guān)系，Rietveld[9]給出了如下關(guān)系。

(2a)

(2b)

Bahel[10]給出了式3所示的關(guān)系。

(3a)

(3b)

該模型在使用中月均日晴天總輻射量較難獲得。Prescott等[11]對該模型進行了修正，如式(4)所示，將日晴天總輻射量用天文輻射替換。天文輻射根據(jù)緯度、赤緯角等信息便于計算。王炳忠等[1]提出采用理想大氣日總輻射量代替天文輻射量，原因是理想大氣輻射量的計算中考慮了海拔和緯度的因素，而海拔因素是影響輻射的重要因素。

(4)

式中：G0為月平均日天文總輻射量

有學者根據(jù)當?shù)貧夂蛱卣?，將日照時數(shù)模型發(fā)展為非線性關(guān)系。Newland[12]在模型中引入了對數(shù)關(guān)系，如式(5)所示。Bakirci[13]提出了指數(shù)關(guān)系的日照時數(shù)模型，如式6所示。

(5)

(6)

Ogelman將日照時數(shù)模型發(fā)展成了二次完全非線性關(guān)系[14], 如式7所示。

(7)

Bahel等[15]在Ogelman模型基礎(chǔ)上將日照時數(shù)模型發(fā)展成為了三次非線性模型。日照時數(shù)模型變得越來越復(fù)雜，隨著模型復(fù)雜性的提高，其地域的適用性比計算準確性的改善更為顯著，即高次非線性的日照時數(shù)模型能夠在更廣泛的地區(qū)適用[3]。

1.1.3 溫差模型 日照時數(shù)模型雖然準確度較高，但是日照時數(shù)并不是常用的氣象參數(shù)，其數(shù)據(jù)獲取有一定的局限性，這一點限制了該模型的廣泛應(yīng)用。氣溫是最常見、也最方便測量的氣象參數(shù)，但是研究表明，最容易獲取的平均氣溫與水平面日總輻射之間并無有效的函數(shù)關(guān)系[16]，而日最高與最低氣溫之差與總輻射之間具有函數(shù)關(guān)系。

Hargreaves等人[17]提出了一個溫差的非線性模型，如式8所示。式中系數(shù)a體現(xiàn)了地域性差異，內(nèi)陸地區(qū)a取值0.16，沿海地區(qū)取值0.19。Allen[18]發(fā)展了Hargreaves模型，在模型中考慮了大氣壓的影響，如式9所示。

(8)

(9)

式中：Kra為經(jīng)驗系數(shù);PS為當?shù)卮髿鈮海琸Pa;P0為標準大氣壓，101.3 kPa。

Annandale等[19]引入了海拔參數(shù)對Hargreaves模型進行了修改，模型如式(10)所示。

(10)

Bristow等[20]提出了指數(shù)形式溫差輻射模型，如式(11)所示。Meza等[21]將公式中的系數(shù)a設(shè)為0.75，c設(shè)為2，系數(shù)b仍為經(jīng)驗系數(shù)，對Bristow模型進行了具體化，這樣可以降低計算誤差。

(11)

陳仁生等[22]提出對數(shù)形式的溫差輻射模型，如式12所示。

(12)

1.1.4 多參數(shù)模型 除了日照時數(shù)、溫差等單氣象參數(shù)以外，還有多參數(shù)構(gòu)成的日總輻射模型，這類模型是以日照時數(shù)或者溫差為主要參數(shù)，綜合了云量、大氣壓、相對濕度等參數(shù)對太陽輻射的影響。

Garg等[23]采用氣溫和降雨對?ngstr?m 模型的經(jīng)驗系數(shù)a、b進行了擬合，如式(13)所示。

a=0.379 1-0.000 4T-0.017 6P

(13a)

b=0.481 0+0.004 3T+0.009 7P

(13b)

式中：T為氣溫，℃;P為降雨量，cm。

陳仁生等[22]提出了溫差和日照時數(shù)的非線性輻射模型，如式(14)所示。曹雯等[24]將該模型中參數(shù)c設(shè)定為1。

(14)

Swartman等[25]提出了日照百分率和相對濕度的輻射模型，如式(15)所示。

(15)

式中：a、b、c為經(jīng)驗系數(shù);RH為相對濕度。

Jong等[26]提出了溫差和降雨量兩參數(shù)組合的輻射模型，如式(16)所示。

(16)

Supit等[27]提出了溫差和云量兩參數(shù)組合的輻射模型，如式(17)所示。

(17)

式中：C為云量

Abdalla等[28]提出了日照時數(shù)，平均氣溫和相對濕度三類參數(shù)的輻射模型，如式(18)所示。

(18)

Ojosu等[29]提出了日照時數(shù)，氣溫極值和相對濕度三類參數(shù)的輻射模型，如式(19)所示。

(19)

1.2 空間插值模型

1.2.1 計算原理 空間插值輻射模型是無輻射測量地域獲取輻射數(shù)據(jù)的另一類重要方法。在一定區(qū)域內(nèi)當氣候具有較好的相似性，而獲得氣象參數(shù)較為困難時，空間插值模型是計算太陽輻射數(shù)據(jù)的較好途徑?？臻g插值模型對于觀測臺站十分稀少而臺站分布又非常不合理的地區(qū)具有十分重要的實際意義[30]。

空間位置上越靠近的點，越可能具有相似的特征值;而距離越遠的點，其特征值相似的可能性越小[31]，這是“地理學第一定律的假設(shè)”，是最早的幾何空間插值技術(shù)基本原理, 距離權(quán)重法(Distance Weighting)屬于幾何空間插值法。空間統(tǒng)計學被引入了空間插值方法，用統(tǒng)計的概念去研究空間中的相近性問題，提出空間相似的程度是通過點對的平均方差度量的[32]?？肆⒏穹?Kriging)屬于空間統(tǒng)計法的空間插值。樣條插值法(Spline methods)屬于函數(shù)類空間插值方法，通過構(gòu)造平滑的函數(shù)曲線來進行插值，不需要對空間結(jié)構(gòu)進行預(yù)估計，也不需要做統(tǒng)計假設(shè)[30]?？臻g插值法多種多樣，但將任何一種插值技術(shù)應(yīng)用于太陽輻射的計算，必須充分考慮其輻射資源的相似性，插值技術(shù)理論假設(shè)和應(yīng)用條件等因素。

1.2.2 距離權(quán)重法 距離權(quán)重法較為簡便，只以兩地距離為依據(jù)進行插值，如式(20)所示[33]。該方法的實質(zhì)是以插值點與采樣點間距離為權(quán)重的一種加權(quán)平均法，其權(quán)重賦予離插值點越近的采樣點賦予估值權(quán)重越大。這對于與緯度、海拔等多種因素相關(guān)的太陽輻射不太合適。

(20)

式中：Z為計算站點的太陽總輻射；Zi為第i個站點的太陽總輻射。

Nalder等[33]提出了距離權(quán)重法的改進方法-梯度距離平方反比法(Gradient Plus Inverse Distance Squared)。在距離權(quán)重的基礎(chǔ)上，本方法考慮了氣象要素隨海拔和經(jīng)緯向的剃度變化。

(21)

式中：Z為計算站點的氣象要素；X、Xi為計算站點與參考氣象站點的經(jīng)度；Y、Yi為氣象站點的Y軸緯度；E、Ei為氣象站點的海拔高度；Cx、Cy和Ce為經(jīng)緯度與海拔高度對應(yīng)的回歸系數(shù)；di為計算站點到第I站點的大地球面距離；n為用于插值的氣象站點的數(shù)目；Zi為氣象站點的氣象要素測量值。

1.2.3 普通克立格法

普通克立格法來源于地統(tǒng)計學中[34]，以區(qū)域化變量理論為基礎(chǔ)，半變異函數(shù)為分析工具，能提供最佳線性無偏估計而逐漸被廣泛運用于需要空間插值的諸多領(lǐng)域，但是計算復(fù)雜且計算量大。其插值公式如式22所示。

(22)

式中：λi為氣象要素的Z(xi)的權(quán)重；Z(xi)為測試值

權(quán)重系數(shù)由“克里格方程組”決定，如式(23)所示。

(23)

式中：C(Xi,Xj)為采樣點間的協(xié)方差；C(Xi,X′)為采樣點與插值點間的協(xié)方差；μ為極小化處理時的拉格朗日乘子

1.2.4 樣條插值法

樣條插值是根據(jù)已知點值來擬合出平滑的樣條函數(shù)，然后使用樣條函數(shù)值作為插值結(jié)果。樣條函數(shù)易操作,計算量不大，多用于氣象要素的時間序列插值。它適合于已知點密度較大的情況，缺點是難以對誤差進行估計，點稀時效果不好。樣條插值是函數(shù)逼近的方法，3次樣條函數(shù)和薄盤光滑樣條函數(shù)是兩類常用的樣條函數(shù)。

3次樣條函數(shù)的定義是[35]:已知平面上n個點(xi,yi)(i=1,2,…,n),其中x1

1)S(xi)=yi,(i=1,2,…,n),

2)S(x)在每個子區(qū)間[xi,xi+1]上為3次多項式

S(x)=ci1(x-xi)3+ci2(x-xi)2+ci3(x-xi)+ci4

(24a)

3)S(x)在整個區(qū)間上有連續(xù)的一階及二階導(dǎo)數(shù)。

則三次樣條插值模型如式(24)所示[36]。

(24b)

薄盤光滑樣條函數(shù)是對樣條函數(shù)法的曲面擴展, 常用于不規(guī)則分布數(shù)據(jù)的多變量平滑插值。利用光滑參數(shù)來達到數(shù)據(jù)逼真度和擬合曲面光滑度之間的優(yōu)化平衡, 保證了插值曲面光滑連續(xù), 且精度可靠。它除通常的樣條自變量外, 允許引入線性協(xié)變量子模型。薄盤光滑樣條函數(shù)如式(25)所示[37]。

Z=f(x)+bTy+e

(25)

式中：Z為位于空間點的插值結(jié)果；x為樣條獨立變量矢量；f為要估算的關(guān)于x未知光滑函數(shù)；y為獨立協(xié)變量矢量；b為y的向量系數(shù)；T為轉(zhuǎn)置符號；e為自變量隨機誤差。

1.3 基于數(shù)字高程的輻射計算模型

1.3.1 計算原理

地形對太陽輻射具有重要影響，坡度、坡向以及周圍地形的遮蔽都會顯著影響水平地面接收到的總輻射，前面介紹的氣候模型和空間插值模型都不能解決復(fù)雜地形下的輻射計算問題。隨著地理信息系統(tǒng)技術(shù)的發(fā)展，數(shù)字高程技術(shù)被用于復(fù)雜地形條件下的輻射計算。

數(shù)字高程模型(Digital Elevation Model,簡稱DEM)是對地球表面地形屬性為高程時的一種離散的數(shù)字表達。通過DEM可以直接獲得地形的坡度、坡向等地形信息,用于計算地形遮擋狀態(tài)下的地面接收到的水平總輻射數(shù)據(jù)。數(shù)字高程的優(yōu)勢表現(xiàn)在坡度、坡向、地形遮蔽度的計算以及模擬結(jié)果可視化表達方面。采用DEM技術(shù)輻射模型主要是考慮地形對輻射的遮蔽作用，用地形遮蔽因子來體現(xiàn)，不同的DEM輻射計算模型主要是地形遮蔽因子的計算方法不同，遮蔽因子可作用于散射或者反射分量的計算。圖1是采用數(shù)字高程模型進行總輻射計算的流程圖。

圖1 采用數(shù)字高程模型進行總輻射計算的流程圖[43]Fig.1 The flow chart of global radiation calculation using DEM

1.3.2 DEM輻射模型

Dozier[38]最早提出了利用數(shù)字高程模型模擬太陽輻射的方法。中國基于DEM的輻射模型起源于對山地地形輻射計算模型的研究。傅抱璞[39]對于任意地形條件下太陽輻射進行了開創(chuàng)性研究。翁篤鳴、李占清[40-41]等發(fā)展了這一方法，之后李新[42]、楊昕[43]等開展了基于DEM技術(shù)的復(fù)雜地形輻射計算

模型研究。

李新等[42]提出了依據(jù)DEM技術(shù)計算我國任意地形條件下太陽輻射模型，模型中利用計算機圖形學的光線追蹤算法生成形狀因子計算地形對坡面的反射輻射。其模型如式(26)所示。

G=Gdir+Gdif+Gref

(26a)

式中：G為水平面總輻射；Gdir為直射輻射；Gdif為散射輻射；Gref為反射輻射。

其中某個坡元j的反射輻射計算模型為:

(26b)

(26c)

式中：Fij為坡元i到坡元j的形狀因子；Ai，Aj為坡元i,j的面積；r為坡元i,j間的距離；φi，φj為坡元i,j法線與它們連線的夾角；HID為取值0或1，取決于第i個坡元能否“看到”第j個坡元，采用光線追蹤法計算。

楊昕等[43]提出了基于DEM的山地總輻射模型，給出了地形遮蔽度因子的計算公式，將其用于散射輻射的計算。模型如式(27)所示。

G=S+D+R+r

(27a)

式中：G為水平面總輻射；S為直射輻射；D為天空散射輻射；R為周圍地形的短波反射輻射；r為研究點與遮蔽物間空氣散射輻射(程輻射)。

其中，天空散射輻射計算模型為(27b)：

(27b)

式中：Kd為地形遮蔽度因子；D0為水平面散射輻射通量密度；α為坡度；β為坡向；F(n)為云量函數(shù)；h′、A′為正午時刻太陽高度角和方位角。

(27c)

2 輻射模型的適用性分析

2.1 其他模型分析

除過文中介紹的氣象參數(shù)模型、空間插值模型和基于DEM的輻射模型之外，還有概率統(tǒng)計模型和衛(wèi)星遙感模型。概率統(tǒng)計模型以人工神經(jīng)網(wǎng)絡(luò)方法[44]和時間序列方法[45]為主要代表，其模型的計算需要大量實測輻射數(shù)據(jù)作為輸入，以統(tǒng)計學原理從海量數(shù)據(jù)中構(gòu)造出經(jīng)驗?zāi)Ｐ?。衛(wèi)星遙感模型是根據(jù)統(tǒng)計反演法或者物理反演法采用衛(wèi)星的可見光和紅外光波段測量數(shù)據(jù)建立的[46]，主要用于研究太陽輻射的垂直空間分布和大氣對輻射吸收等問題。上述兩類輻射模型尚處于發(fā)展階段，技術(shù)趨待完善，且因計算復(fù)雜計算結(jié)果存在較大的不確定性，因此，很少應(yīng)用于建筑節(jié)能分析。

2.2 輻射模型節(jié)能計算的適用性分析

太陽輻射數(shù)據(jù)在建筑節(jié)能領(lǐng)域應(yīng)用廣泛，而且具有重要作用，面對眾多輻射計算模型，對模型的選擇是個重要問題。從計算簡便性和準確性來看，氣象參數(shù)輻射模型最簡單，而且其計算準確度也高，劉大龍等[47]比較了日照時數(shù)和梯度距離平方反比法模型，研究表明日照時數(shù)模型在我國地域內(nèi)更為準確。該模型由左大康最早應(yīng)用于中國[48]，并得到廣泛應(yīng)用，已成為中國氣象行業(yè)的標準輻射計算方法[49]。鄧艷君等[50]采用中國實測輻射數(shù)據(jù)對三種氣象參數(shù)輻射模型進行了對比，日照百分率和溫差模型精度高，但經(jīng)典日照模型更為穩(wěn)定，而溫差模型誤差較大。如果測量的氣象參數(shù)具有較長的時間序列，則可得到長時間序列的輻射值，這是氣象參數(shù)模型的另一優(yōu)點。

日照時數(shù)等氣象參數(shù)類模型主要是用于晴天條件下輻射的計算，非晴天條件下該類模型誤差較大[5]。當需要計算有輻射值的相近站點的輻射值時，空間插值模型是較好的選擇?？臻g插值計算輻射有直接法和間接法兩種應(yīng)用方式，直接法就是對輻射數(shù)據(jù)進行插值；間接法可以對氣象參數(shù)模型中的相關(guān)參數(shù)進行插值，然后應(yīng)用再通過計算得到輻射數(shù)據(jù)。然而, 在眾多的氣象要素空間插值方法中, 并沒有一種適合每一個氣象要素的普適的最佳插值方法[51]。梯度距離平方反比法包含了經(jīng)緯度和海拔信息，比較適合輻射的直接插值；普通克立格法能夠準確控制計算誤差，但計算較為復(fù)雜，研究表明在溫度插值時其準確度和梯度距離平方反比法相當[51]；樣條插值比較適合氣溫和氣象要素的時間序列插值。

當缺乏用于計算輻射的氣象參數(shù)時，且地形對輻射具有重要影響時可采用基于DEM的輻射模型，這類模型根據(jù)天文輻射和地理要素可計算出不同季節(jié)、不同時刻的輻射值，而且還便于計算直射輻射和散射輻射。DEM輻射模型已成為我國研究復(fù)雜地形輻射數(shù)據(jù)主要方法[52-53]。DEM輻射模型具有一個非常適合建筑節(jié)能分析的特點，就是能夠計算建筑周圍微環(huán)境的輻射值。建筑節(jié)能分析絕大多數(shù)情況需要的是建筑周圍微環(huán)境的輻射數(shù)據(jù)。建筑更多的集中于城市，人為因素導(dǎo)致城市中出現(xiàn)了特殊的城市氣候，而城市氣候一個顯著特點就是差異性、下墊面、建筑布局、綠化等眾多因素導(dǎo)致城市不同區(qū)域具有明顯不同的微氣候環(huán)境，這就是城市氣候的差異性。因此，準確計算建筑周圍微環(huán)境對于建筑節(jié)能分析具有重要作用。

2.3 節(jié)能分析用輻射模型展望

建筑節(jié)能技術(shù)主要應(yīng)用于城市建筑，因此，城市中微環(huán)境的輻射數(shù)據(jù)成為節(jié)能分析用輻射數(shù)據(jù)的重點。隨著我國氣象觀測能力的增強，輻射臺站密度有顯著增加，但其與城市的快速擴張速度相比相差甚遠，因此城市中微環(huán)境的輻射數(shù)據(jù)獲取依然離不開模型計算方法。節(jié)能分析用太陽輻射模型的發(fā)展趨勢是建筑周邊微氣候的輻射計算，輻射模型將會走氣象參數(shù)模型與DEM模型相結(jié)合的道路。通過氣象參數(shù)獲得大范圍、長時間序列的基本輻射數(shù)據(jù)，然后從DEM中依據(jù)地形、考慮建筑間的遮擋進行微環(huán)境、小尺度空間的精確化輻射計算。在采用DEM模型時，將坡元對輻射的影響用相鄰建筑來替代，這是城市環(huán)境輻射計算的新特征。

[1] 王炳忠.我國的太陽能資源及其計算[J].太陽能學報,1980,1(1):1-9.

Wang B Z.Solar Energy resources in China [J]. Acta Energiae Solaris Sinica,1980,1(1):1-9.(in Chinese)

[2] Black J N. The distribution of solar radiation over the earth’s surface [J]. Archiv fur Meteorologie，Geophysik，und Bioklimatologie Serie A Meteorologie und Geophysik 1956,7：165-189.

[3] Besharat F，Ali A,Dehghan A R,et al.Empirical models for estimating global solar radiation：A review and case study [J]. Renewable and Sustainable Energy Reviews，2013,21:798-821.

[4] ?ngstr?m A.,Solar and terrestrial radiation [J]. Quarterly Journal of Royal Meteorological Society，1924,50:121-125.

[5] 孫治安，施俊榮，翁篤鳴.中國太陽總輻射氣候計算方法的進一步研究[J].南京氣象學院學報，1992,15(2):21-28.

Sun Z A，Shi J R，Weng D M. A further research on the climatological calculation method of the global solar radiation over China [J].Journal of Nanjing Institute of Meteorology，1992,15(2):21-28.(in Chinese)

[6] 高國棟，陸瑜蓉.中國地表面輻射平衡與熱量平衡[M].北京：科學出版社，1982.

[7] Soler A.Monthly specific Rietveld’s correlations [J]. Solar and Wind Technology，1990，7：305-312.

[8] 鞠曉慧，屠其璞，李慶祥.我國太陽總輻射氣候?qū)W計算方法的再討論[J].南京氣象學院學報，2005,28(4):516-521.

Ju X H，Tu Q P，Li Q X.Discussion on the Climatological calculation of solar radiation [J]. Journal of Nanjing Institute of Meteorology，2005,28(4):516-521.(in Chinese)

[9] Rietveld M. A new method for estimating the regression coefficients in the formula relating solar radiation to sunshine [J]. Agricultural Meteorology,1978，19:243-252.

[10] Zabara K. Estimation of the global solar radiation in Greece [J]. Solar and Wind Technology，1986，3(4):267-272.

[11] Prescott J A. Evaporation from water surface in relation to solar radiation [J].Transactions of the Royal Society of Australia,1940,46:114-121.

[12] Newland F J. A study of solar radiation models for the coastal region of South China [J].Solar Energy,1988，31:227-235.

[13] Bakirci K. Correlations for estimation of daily global solar radiation with hours of bright sunshine in Turkey [J]. Energy，2009,34:485-501.

[14] Ogelman H，Ecevit A，Tasdemiroglu E.A new method for estimating solar radiation from bright sunshine data [J]. Solar Energy,1984,33:619-625.

[15] Bahel V, Bakhsh H，Srinivasan R.A correlation for estimation of global solar radiation [J]. Energy,1987，12:131-135.

[16] 劉大龍.區(qū)域氣候預(yù)測與建筑能耗演化規(guī)律研[D].西安：西安建筑科技大學，2010.

[17] Hargreaves G H，Samani Z A.Estimating potential evapotranspiration [J]. Journal of Irrigation and Drainage Engineering,1982,108(IR3):223-230.

[18] Allen R.Evaluation of procedures of estimating mean monthly solar radiation from air temperature [R]. Rome：FAO,1995.

[19] Annandale J G，Jovanic N Z，Benade N，et al.Software for missing data error analysis of penman-monteith reference evapotranspiration [J]. Irrigation Science,2002,21:57-67.

[20] Bristow K L，Campbell G S.On the relationship between incoming solar radiation and daily maximum and minimum temperature [J]. Agricultural and Forest Meteorology,1984,31:159-166.

[21] Meza F，Varas E.Estimation of mean monthly solar global radiation as a function of temperature [J]. Agric. For Meteorol.,2000,100:231-241.

[22] Chen R，Ersi K，Yang J，et al.Validation of five global radiation models with measured daily data in China [J].Energy Conversion and Management,2004，45:1759-69.

[23] Gariepy J.Estimation of global solar radiation [R]. International Report，Service of Meteorology，Governm-ent of Quebec，Canada; 1980.

[24] 曹雯，申雙和.我國太陽日總輻射計算方法的研究[J].南京氣象學院學報，2008,31(4):587-591.

Cao W，Shen S H. Estimation of daily solar radiation in China[J].Journal of Nanjing Institute of Meteorology，2008,31(4):587-591.(in Chinese)

[25] Swartman R K，Ogunlade O.Solar radiation estimates from common parameters [J].Solar Energy，1967,11:170-172.

[26] Jong D R，Stewart D W.Estimating global solar radiation from common meteorological observations in western Canada [J]. Canadian Journal of Plant Science，1993,73:509-518.

[27] Supit I，Kappel V R R.A simple method to estimate global radiation [J].Solar Energy，1998,63：147-160.

[28] Abdalla Y A G.New correlation of global solar radiation with meteorological parameters for Bahrain [J]. International Journal of Solar Energy,1994，16:111-120.

[29] Ojosu J O，Komolafe L K.Models for estimating solar radiation availability in south western Nigeria [J]. Nigerian Journal of Solar Energy,1987，6:69-77.

[30] 李新，程國棟，盧玲.空間內(nèi)插方法比較[J].地球科學進展，2000,15(3):260-264.

Li X，Chen G D，Lu L.Comparison of spatial interpolation methods [J].Advance in Earth Sciences，2000,15(3):260-264.(in Chinese)

[31] 烏仔倫，劉瑜，張晶，等.地理信息系統(tǒng)原理方法和應(yīng)用[M].北京:科學出版社，2001.

[31] Wu Z L，Liu Y，Zhang J，et al. Geographic Information System：Principle，Method and Application [M].Beijing：Science Press，2001.

[32] Goovaerts P.Geostatistics for natural resource evaluation [M]. New York：Oxford University Press，1997.

[33] Nalder I A，Wein R W.Spatial interpolation of climate normals:test of a new method in the Canadian boreal forest [J]. Agric. For Meteorol.,1998,92:211-225.

[34] Journel A G，Huijbregts C J.Mining Geostatistics [M]. London:Academic Press，1978.

[35] 樊天鎖,芮兵.樣條插值的MATLAB實現(xiàn)[J].佳木斯大學學報：自然科學版，2011,29(2)：238-210.

Fan T S，Rui B. The realization of spline interpolation based on MATLAB [J]. Journal of Jiamusi University：Natural Science Edition，2011,29(2):238-210.(in Chinese)

[36] Bi C X,Geng L,Zhang X Z.Cubic spline interpolation-based time-domain equivalent source method for modeling transient acoustic radiation [J].Journal of Sound and Vibration，2013,332：5939-5952.

[37] 鄭小波，羅宇翔，于飛.西南復(fù)雜山地農(nóng)業(yè)氣候要素空間插值方法比較[J].中國農(nóng)業(yè)氣象，2008，29(4)：458-462.

Zheng X B，Luo Y X，Yu F. Comparisons of spatial interpolation methods for agro-climate factors in complex mountain areas of southwest China [J]. Chinese Journal of Agrometeorology，2008，29(4)：458-462.(in Chinese)

[38] Dozier J，F(xiàn)rew J.Rapid calculation of terrain parameters for radiation modeling from digital elevati-on data [J].IEEE Transaction on Geo-science and Remote Sensing,1990,28(5):963-969.

[39] 傅抱璞，虞靜明，盧其堯.山地氣候資源與開發(fā)利用[M].南京:南京大學出版社,1996. 8-39.

[40] 翁篤鳴.中國輻射氣候[M]. 北京:氣象出版社，1997，250-278.

[41] 李占清，翁篤鳴.丘陵山地總輻射的計算模式[J].氣象學報,1988,46(4):461-468.

Li Z Q，Weng D M. A numerical approach toward global radiation over rugged areas [J]. Acta Mete-orologica Sinica，1988,46(4):461-468.(in Chinese)

[42] 李新，程國棟，陳賢章，等.任意地形條件下太陽輻射模型的改進[J].科學通報，1995,44(9)：993-998.

Li X，Chen G D，Chen X Z，et al. Impr-ovements of solar radiation models in terrain condition-s [J]. Chinese Science Bulletin，1995,44(9)：993-998.(in Chinese)

[43] 楊昕，湯國安，王雷.基于DEM 的山地總輻射模型及實現(xiàn)[J].地理與地理信息科學，2004，20(5):41-44.

Yang X，Tang G A,Wang L，Modeling of global radiation over rugged areas based on DEM [J]. Geography and Geo-information Science，2004,20(5):41-44.(in Chinese)

[44] Yadav A K，Chande S S.Solar radiation prediction using Artificial neural network techniques：A review [J]. Renewable and Sustainable Energy Reviews,2014,33(5)：772-781.

[45] Voyant C，Paoli C，Muselli M，et al.，Multihorizon solar radiation forecasting for mediterranean locations using time series models [J]. Renewable and Sustainable Energy Reviews,2013，28(12)：44-52.

[46] 申彥波.近20 年衛(wèi)星遙感資料在我國太陽能資源評估中的應(yīng)用綜述[J].氣象，2010,36(9):111-115.

Shen Y B.Review of applications of satellite remote sensing data to solar energy resources assessment in China in recent 20 years [J]. Meteorological Monthly,2010,36(9):111-115.(in Chinese)

[47] Liu D L，Yang L，Liu J P.Compared Predict Model of Global Solar Radiation [R].Solari-S 2008，Hongkong,2008.12,.

[48] 左大康，王懿賢，陳建綏.中國地區(qū)太陽總輻射的空間分布特征[J].氣象學報,1963,33(1):78-96.

Zuo D K，Wang Y X，Chen J S.Characteristics of the distribution of total radiation in China [J]. Acta Meteorologica Sinica,1963,33(1):78-96.(in Chinese)

[49] 中華人民共和國氣象行業(yè)標準QX/T89-2008，太陽能資源評估方法[S].北京：氣象出版社，2008.

[50] 鄧艷君，邱新法，曾燕，等.幾種水平面太陽總輻射量計算模型的對比分析[J].氣象科學，2013，33(4):371-377.

Deng Y J，Qiu X F，Zeng Y，et al.Comparison of horizontal global solar radiation mod-els [J].Journal of Meteorological Sciences，2013，33(4):371-377.(in Chinese)

[51] 林忠輝，莫興國，李宏，等.中國陸地區(qū)域氣象要素的空間插值[J].地理學報，2002，57(1)：47-55.

Lin Z H，Mo X G，Li H，et al.Compar-ison of three Spatial interpolation methods for climate variables in China [J].Acta Geographica Sinica，2002，57(1)：47-55.

[52] 李凈，李新.基于DEM的坡地太陽總輻射估算[J].太陽能學報，2007,28(8):905-911.

Li J，Li X. Estimation solar radiation on slop-e using DEM [J]. Acta Energiae Solaris Sinica，2007,28(8):905-911.(in Chinese)

[53] 張海龍，劉高煥，姚玲，等.復(fù)雜地形任意天氣情形下太陽直射輻射量模擬[J].中國沙漠，2010,30(6)：1409-1407.

Zhang H L，Liu G H，Yao L，et al.Direct solar radiation simulation in all-weather and complex terrain conditions [J].Journal of Desert Research，2010，30(6)：1409-1407.(in Chinese)

(編輯 胡 玲)

Review of globe solar radiation model for building energy efficiency analysis

LiuDalong1，2，YangLiu1，2，HuoXujie1，2，LiuJiaping1，2
(1.Architecture school，Xi’AN University of Architecture and Technology，710055 Xi’an，P.R.China;2.Shaanxi provincial collaborative innovation centre of green building in western，710055 Xi’an，P.R.China)

The solar radiation is an important fundamental meteorological parameter for building energy efficiency analysis. Since the measured radiation data can not meet the demand for analysis，the theoretical calculation is the critical. The widely used horizontal solar radiation models were classified into three categories including meteorological parameter model，spatial interpolation model，DEM solar model. Principles and calculation methods of three type model were presented. The applicability of three models in building energy efficiency field was analyzed. The future of solar radiation for building energy efficiency analysis was predicted.

horizontal globe solar radiation；meteorological parameter model；spatial interpolation model；DEM solar model

10.11835/j.issn.1674-4764.2015.02.016

2014-10-25 基金項目：“十二五”國家科技計劃項目(2013BAJ03B04-01)；博士后基金面上項目(2014M552419)

劉大龍(1976-)，男，博士，副教授，主要從事建筑能耗模擬研究，(E-mail)coffeevc@xauat.edu.cn。

Foundation item：“Twelfth Five Year Plan” National Science and Technology Plan Project(No.2013BAJ03B04-01)；Post Doctoral Science Fund Projects(No.2014M552419)

TU111.3

A

1674-4764(2015)02-0101-08

Received：2014-10-25

Author brief：Liu Dalong(1976-),PhD.,associate professor,main research interest:building energy simulation，(E-mail)coffeevc@xauat.edu.cn.