國家發(fā)展和改革委員會(huì)能源研究所 ■ 王斯成
對(duì)于傾斜光伏方陣面輻射量的計(jì)算,國內(nèi)外已有很多研究,但是絕大多數(shù)文獻(xiàn)僅限于簡單的固定方陣[1,2],雖然也有對(duì)于太陽跟蹤方陣的研究[3](被RetScreen采用),但并未系統(tǒng)地對(duì)不同坐標(biāo)系的太陽跟蹤器進(jìn)行定義和區(qū)分,也沒有對(duì)于各種光伏方陣和太陽跟蹤器給出完整的方陣面輻照度計(jì)算的數(shù)學(xué)模型。本文在文獻(xiàn)[4]的基礎(chǔ)上,對(duì)地平坐標(biāo)系和赤道坐標(biāo)系太陽跟蹤器進(jìn)行了定義和分類,從太陽和地球的運(yùn)行規(guī)律及球面三角的基本原理推導(dǎo)出了不同坐標(biāo)系下方陣面太陽入射角、方陣傾角和方陣面輻照度的完整計(jì)算方法和計(jì)算公式,為正確計(jì)算光伏方陣面(包括背表面)輻照度和輻射量提供了理論依據(jù)。
為了推導(dǎo)光伏方陣面輻照度的數(shù)學(xué)模型,必須先了解下面幾個(gè)基本定律:球面三角邊的余弦定律、太陽輻射的直散分離原理、太陽高度角的計(jì)算公式和太陽入射光的余弦定律。
球面三角邊的余弦公式[5]可用文字表達(dá)為:球面三角形任意邊的余弦等于其他兩邊余弦的乘積加上這兩邊的正弦及其夾角余弦的連乘積。
圖1為球面三角形的邊角示意圖,其中,a、b、c代表邊,A、B、C代表角。
圖1 球面三角形的邊角示意圖
邊a的余弦定律為:
余弦定律為:任意傾斜面的輻照度同該表面法線與入射線方向之間夾角的余弦成正比。圖2為太陽直接輻照度的關(guān)系圖,其中,SH′為水平面的太陽直接輻照度;SD′為太陽直接輻照度;ST′為傾斜方陣面的太陽直接輻照度;θ為任意時(shí)刻太陽直射光的入射角;θZ為太陽天頂角,是太陽射線與天頂軸的夾角;α為太陽高度角,是太陽射線與地平面的夾角,范圍在0°~90°;Z為光伏方陣傾角。
圖2 太陽直接輻照度的關(guān)系圖
由圖2可推導(dǎo):
所以,
太陽高度角α的表達(dá)式可由球面三角原理導(dǎo)出[3]。圖3為太陽高度角和地平面的關(guān)系圖,其中,與水平面垂直的軸是天頂軸,天頂軸和太陽射線在天球上可以組成一個(gè)球面三角形。根據(jù)球面三角邊的余弦定律(見式1),可以導(dǎo)出太陽高度角α的正弦公式(見式5)[3]。
圖3 太陽高度角和地平面的關(guān)系圖
由于:
a:太陽天頂角θZ,
所以:
因此,代入式(1),可以導(dǎo)出太陽高度角α的正弦表達(dá)式為:
式中,φ為當(dāng)?shù)鼐暥?0?~90?),是所在地法線與地心的連線與赤道平面的夾角,有北緯、南緯之分;ω為時(shí)角(地球自轉(zhuǎn)一周360?,24 h),即15?/h或1?/4min,正午為零,上午為正,下午為負(fù);δ為太陽赤緯角,是陽光射線與赤道平面的夾角,太陽照射到北半球時(shí)為正,照射到南半球時(shí)為負(fù),春秋分時(shí)為零,一年中太陽赤緯角的變化規(guī)律如圖4所示。
著名的計(jì)算太陽赤緯角的Cooper方程為:
式中,N為從1月1日起算的天數(shù)。
圖4 一年中太陽赤緯角的變化圖
水平面和傾斜面上所接收到的輻射量均符合直散分離原理,即總輻射等于直接輻射與散射輻射之合,只不過水平面所接收到的輻射量沒有地面反射分量,而傾斜方陣面上所接收到的輻射量包括地面反射分量。
光伏方陣面上所接收到的全部輻射量的公式為:
式中,QH為水平面接收到的日總輻射量;SH為水平面接收到的日直接輻射量;DH為水平面接收到的日散射輻射量;QT為傾斜方陣面接收到的日總輻射量;ST為傾斜方陣面接收到的日直接輻射量;DT為傾斜方陣面接收到的日散射輻射量;RT為傾斜方陣面接收到的日地面反射輻射量。
同理,傾斜方陣面上的總輻照度QT′可表達(dá)為[1]:
式中,DT′為傾斜方陣面上的散射輻照度;RT′為傾斜方陣面上的地面反射輻照度。
由式(2)和式(4)可得:
傾斜方陣面上的散射輻照度DT′ 可表達(dá)為:
式中,DH′為水平面的散射輻照度;Z′為任意時(shí)刻的光伏方陣傾角(0?~90?)。
傾斜方陣面上的地面反射輻照度RT′可表示為:
式中,ρ為地面反射率;QH′為水平面的總輻照度。
由上面的公式可知,光伏方陣面上所接收到的散射輻射和地面反射被認(rèn)為是各向同性的,與太陽光的入射角度無關(guān)。而光伏方陣所接收到的天空散射與光伏方陣面的視天空比例有關(guān):光伏方陣水平向上,光伏方陣傾角Z為0°,視天空為100%,DT′=DH′;光伏方陣水平向下,Z為180°,cosZ=-1,視天空為零,DT′=0;方陣垂直安放,則Z為 90°,視天空為 50%,DT′=1/2DH′。同樣,光伏方陣所接收到的地面反射與光伏方陣面的視地表比例有關(guān),不再贅述。
由式(9)~式(12),光伏方陣面總輻照度可描述為[1]:
由此可知,光伏方陣面總輻照度的計(jì)算需要5個(gè)變量:水平面的總輻照度QH′;水平面的太陽直接輻照度SH′;水平面的散射輻照度DH′;任意時(shí)刻太陽直射光的入射角θ;任意時(shí)刻的光伏方陣傾角Z′。
上述散射輻射模型假設(shè)天空散射各向同性,但經(jīng)過實(shí)際測(cè)量和驗(yàn)證,該假設(shè)僅適用于低輻照、多云和陰天情況,而對(duì)于晴朗天氣,應(yīng)當(dāng)采用各向異性的散射輻射模型。而對(duì)于各向異性的散射輻射模型,主要有Hay模型[6]、Klucher模型[7]和Perez模型[8-9](被PVSyst采用),楊金煥也有相關(guān)論文述及[10-11]。
到底采用各向同性還是各向異性模型,以及到底采用哪一種各向異性模型,并不是本論文的研究重點(diǎn),設(shè)計(jì)者可以根據(jù)自己的情況來確定需要選擇采用哪種模型。
即使采用各向異性散射模型,仍然需要獲得本文1.4節(jié)所列出的相關(guān)變量,主要散射輻射各向異性模型的公式如下:
1) Hay 模型[6,8,10]:
式中,K為地面日總輻射量與大氣層上界日輻射量之比(復(fù)合大氣透明度)。K可用式 (15)求得:
式中,Q0為大氣層上界日總輻射量。
式中,S0為太陽常數(shù) 1367 W/m2;λ為日-地距離變化引起大氣層上界太陽輻照度的修正值。
式中,n為從1月1日算起的天數(shù)。
2) Klucher模型[7]:
3) Perez模型(被PVSyst所采用):
1986年模型[8]為:
或1990年的簡化模型[9]為:
由各模型公式可以看出,這些各向異性模型傾斜面散射輻照度DT′ 也都是DH′ 、θ和Z′ 的函數(shù),即DT′ =f(DH′ ,θ,Z′ ),只要有了DH′、θ和Z′ ,即使采用散射輻射各向異性模型,也可以順利得出傾斜面散射輻照度。
各向異性散射輻射的詳細(xì)計(jì)算過程請(qǐng)參考相關(guān)文獻(xiàn)。(注:式(19)、式(20)中的F1、F2、a、b、c、d等都是原著中的經(jīng)驗(yàn)參數(shù),并非本文中使用的參數(shù),詳細(xì)的參數(shù)定義請(qǐng)見原著。)
從上節(jié)可知,水平面總輻照度QH′ 、水平面太陽直接輻照度SH′ 和水平面散射輻照度DH′ 是計(jì)算水平面太陽輻照度的基礎(chǔ)。水平面太陽輻照度根據(jù)不同的數(shù)據(jù)來源,可通過下列幾種方法獲得。
多年平均的水平面小時(shí)總輻射量和小時(shí)散射輻射量可以從當(dāng)?shù)貧庀笳精@得,如果當(dāng)?shù)貧庀笳静荒芴峁┒嗄昶骄r(shí)輻射量的數(shù)據(jù),也可以從NASA、NREL、Mateonorm或PVSyst的數(shù)據(jù)庫中得到;小時(shí)直接輻射量是小時(shí)總輻射量和小時(shí)散射輻射量的差值。
小時(shí)輻射量也可以看成是該小時(shí)的平均輻照度,對(duì)于固定光伏方陣的計(jì)算已足夠,但對(duì)于太陽跟蹤器方陣面輻照度的計(jì)算來說就遠(yuǎn)遠(yuǎn)不夠了。如果需要更細(xì)的積分步長,可以通過插值的方法得到以分鐘為步長的輻照度。
很多氣象站不能提供多年平均小時(shí)輻射量的實(shí)測(cè)數(shù)據(jù),只能提供多年平均日總輻射量和日散射輻射量的數(shù)據(jù)。對(duì)于這種情況,我們可以采用文獻(xiàn)[4]或文獻(xiàn)[12]的方法獲得小時(shí)輻射量(文獻(xiàn)[12]的方法被RetScreen軟件[13]所采用),即以小時(shí)為步長的輻照度。
文獻(xiàn)[12]的計(jì)算模型為:
式中,QH′為水平面小時(shí)平均總輻照度;DH′為水平面小時(shí)散射輻照度;SH′為水平面小時(shí)直接輻照度;rt為水平面小時(shí)總輻射與日總輻射的比值;rd為水平面小時(shí)散射量與日散射量的比值;ωr為日出時(shí)角;a、b均為原著中的經(jīng)驗(yàn)系數(shù),
大多數(shù)氣象臺(tái)站只能提供多年平均水平面月總輻射量,多年平均水平面日總輻射量QH可以通過各月平均水平面日輻射量采用插值的方法取得[4]。有了多年平均水平面日總輻射量QH,就可以通過文獻(xiàn) [14]的方法(被RetScreen所采用)獲得多年平均水平面日散射輻射量DH,模型如下:
有了QH和DH,就可以通過式(21)~式(23)求得QH′、DH′和SH′。
為了求出任意時(shí)刻太陽直射光的入射角的余弦公式,必須先要了解2個(gè)不同的坐標(biāo)系——地平坐標(biāo)系和赤道坐標(biāo)系,由此,太陽跟蹤器可以分為地平坐標(biāo)太陽跟蹤器和赤道坐標(biāo)太陽跟蹤器。任意時(shí)刻太陽直射光的入射角在不同的坐標(biāo)系中有著不同的描述,下面將分別進(jìn)行介紹。
在地平坐標(biāo)系中,光伏方陣將地平面作為參照系。太陽的位置由太陽高度角和太陽方位角來確定。光伏方陣和太陽射線的關(guān)系如圖5所示[3],其中,β為太陽方位角,是太陽射線在地面上的投影與正南方向的夾角,正南方向?yàn)榱?,東為正,西為負(fù);γ為光伏方陣任一時(shí)刻方位角,是方陣法線在水平面上的投影與正南方向的夾角,正南方向?yàn)榱?,東為正,西為負(fù)。
圖5 光伏方陣與太陽射線在地平坐標(biāo)中的關(guān)系
由于:
所以:
因此,代入式(1),得到:
地平坐標(biāo)系光伏方陣的安裝方式有3種,分別為:固定式光伏方陣(含傾角可以按照季節(jié)手動(dòng)調(diào)節(jié)的方陣)、太陽方位角跟蹤光伏方陣和雙軸跟蹤光伏方陣(同時(shí)跟蹤太陽的高度角和方位角)。圖6為3種安裝方式的示意圖。
圖6 地平坐標(biāo)系光伏方陣的運(yùn)行方式
1)對(duì)于固定式光伏方陣,方陣傾角全年固定或按照季節(jié)可調(diào)(在一天當(dāng)中也是固定不變化),方陣方位角一般朝向正南,則有Z′=Z,γ=0;代入式(27),可得:
2)東西跟蹤時(shí),光伏方陣的傾角不變,只跟蹤太陽的方位角,則有Z′=Z,γ=β;代入式 (27),可得:
3)全跟蹤時(shí),光伏方陣始終跟蹤太陽的高度角和方位角,則有Z′= 90°-α,γ=β;代入式 (27),可得:
即入射角θ=0°,始終準(zhǔn)確跟蹤太陽。
赤道坐標(biāo)系是以赤道平面和南北地軸作為參照系。在赤道坐標(biāo)系中,太陽的位置由太陽赤緯角和太陽時(shí)角決定。按照赤道坐標(biāo)系設(shè)計(jì)的用來跟蹤太陽赤緯角和太陽時(shí)角的太陽跟蹤器分為4種類型:平單軸東西向太陽跟蹤器、平單軸帶傾角東西向太陽跟蹤器、極軸太陽跟蹤器和雙軸太陽跟蹤器(同時(shí)跟蹤太陽赤緯角和太陽時(shí)角)。圖7為4種赤道坐標(biāo)太陽跟蹤器的跟蹤示意圖,圖中,Ω為赤道坐標(biāo)系中光伏方陣主軸的旋轉(zhuǎn)角,正南為零,東旋為正,西旋為負(fù)。
圖7 赤道坐標(biāo)系光伏方陣的運(yùn)行方式
平單軸東西向跟蹤可以看成將光伏組件平鋪在地面,方陣法線與天頂軸平行,平單軸光伏方陣和太陽射線的關(guān)系類似于太陽射線和地平面的關(guān)系(見前文圖3和式(5)),如圖8所示。
由于:
A=ω-Ω(光伏方陣需要東西向旋轉(zhuǎn)) ;
所以:
圖8 平單軸光伏方陣與太陽射線的關(guān)系
代入式(1),可以得到平單軸東西向太陽跟蹤器的θ的表達(dá)式為:
赤道坐標(biāo)系其他類型的太陽跟蹤器與平單軸東西向太陽跟蹤器的差別僅僅是主軸的傾角不同而已。為了更清楚地推導(dǎo)θ的公式,圖9給出了光伏方陣在球面三角坐標(biāo)中的圖示。圖中,Z在赤道坐標(biāo)系中為光伏方陣主軸傾角,極軸跟蹤時(shí)有Z=φ;z在赤道坐標(biāo)系中為光伏組件與主軸的夾角,與主軸平行時(shí)為零,南傾為正,北傾為負(fù)。
圖9 赤道坐標(biāo)系光伏方陣與太陽射線的關(guān)系
由于:
所以:
代入式(1),可以得到赤道坐標(biāo)太陽跟蹤器的θ的通用表達(dá)式為:
1)平單軸東西向跟蹤時(shí):
理想條件下,Ω=ω,A=0°,Z=0°,z=0°;于是:
所以,θ=φ-δ。
2)平單軸帶傾角東西向跟蹤時(shí):
理想條件下,Ω=ω,A=0°,Z=0°;于是:
所以,θ=φ-z-δ。
3)極軸跟蹤時(shí)有:
理想條件下:Ω=ω,A=0°,Z=φ,z=0°;于是:
所以,θ=δ。
4)雙軸跟蹤時(shí)有:
理想條件下:Ω=ω,A=0°,Z=φ,z=-δ(春分之后δ為正值,但傾角減小;秋分之后δ為負(fù)值,但傾角增加);于是:
所以,θ=0°,意味著光伏方陣準(zhǔn)確跟蹤太陽。
如果太陽跟蹤器固定不跟蹤時(shí),θ的表達(dá)式應(yīng)當(dāng)與地平坐標(biāo)系θ的表達(dá)式一致,以下來證明這一點(diǎn)。
1)平單軸東西向太陽跟蹤器固定平放時(shí):
地平坐標(biāo)系有:Z′=0°,γ=0°;代入式 (27),得到 cosθ= sinα。
赤道坐標(biāo)系有:Z=0°,z=0°,Ω=0°;代入式(31),得到 cosθ=sinφsinδ+cosφcosδcosω;
根據(jù)地平坐標(biāo)和赤道坐標(biāo)的參數(shù)變換式[4,15]:
可得 cosθ= sinα。
證實(shí)了赤道坐標(biāo)系任意時(shí)刻太陽直射光的入射角與地平坐標(biāo)系任意時(shí)刻太陽直射光的入射角的表達(dá)式一致的結(jié)論。
2)平單軸帶傾角東西向太陽跟蹤器固定安放時(shí):
地平坐標(biāo)系有:Z′=Z,γ=0°;代入式 (27),得到 cosθ=cosZsinα+sinZcosαcosβ。
根據(jù)參數(shù)變換式(33)~式(35)可得:
赤道坐標(biāo)系有:Z=0°,z=Z,Ω=0°;代入式 (32),得到:
證實(shí)了赤道坐標(biāo)系任意時(shí)刻太陽直射光的入射角與地平坐標(biāo)系任意時(shí)刻太陽直射光的入射角的表達(dá)式一致的結(jié)論。
3)極軸太陽跟蹤器固定安放時(shí):
地平坐標(biāo)系有:Z′=φ,γ=0°;代入式 (27),得到 cosθ=cosφsinα+sinφcosαcosβ。
根據(jù)參數(shù)變換式(33)~式(35)可得:
赤道坐標(biāo)系有:Z=φ,z=0°,Ω=0°;代入式(32),得到 cosθ=cos90°sinδ+sin90°cosδcosω=cosδcosω。
證實(shí)了赤道坐標(biāo)系任意時(shí)刻太陽直射光的入射角與地平坐標(biāo)系任意時(shí)刻太陽直射光的的表達(dá)式一致的結(jié)論。
上述結(jié)果可說明這2個(gè)坐標(biāo)系任意時(shí)刻太陽直射光的入射角θ的余弦公式完全正確。
從式(13)可知,要想準(zhǔn)確計(jì)算光伏方陣面的總輻照度,必須知道任意時(shí)刻的光伏方陣傾角Z′,下面分別就2個(gè)坐標(biāo)系進(jìn)行推導(dǎo)。
1)對(duì)于固定式光伏方陣,方陣傾角始終不變,總有Z′=Z;
2)對(duì)于太陽方位角跟蹤光伏方陣,方陣傾角也始終不變,總有Z′ =Z;
3)對(duì)于雙軸跟蹤光伏方陣,方陣傾角始終等于太陽高度角的余角,有Z′ = 90°-α。
為了推導(dǎo)赤道坐標(biāo)系任意時(shí)刻的光伏方陣傾角,首先要給出光伏方陣在天球上的定位三角形[4],如圖10所示。
由于:
圖10 赤道坐標(biāo)系光伏方陣在天球上的定位三角形
A=Ω;
在赤道坐標(biāo)系中,Z為光伏方陣主軸傾角,z為光伏組件在主軸上的傾角。
所以:
代入式(1),可以得到赤道坐標(biāo)系任意時(shí)刻的光伏方陣傾角的表達(dá)式為:
1)平單軸東西向跟蹤時(shí),Z=0°,z=0°;所以cosZ′= cosΩ,即Z′=Ω。
2)平單軸帶傾角東西向跟蹤時(shí),Z=0?;所以cosZ′= coszcosΩ。
3)極軸跟蹤時(shí),Z=φ,z=0?;所以 cosZ′=cosφcosΩ。
4)雙軸跟蹤時(shí),Z=Ω,z=-δ,Ω=ω;所以cosZ′= sin(-δ)sinφ+cos(-δ)cosφcosω。
有了水平面總輻照度QH′、水平面散射輻照度DH′和水平面直接輻照度SH′,就可以根據(jù)式(10)~式 (12)求得傾斜面輻照度QT′、DT′和ST′;然后再對(duì)時(shí)角積分,即可得到每日的輻射量。由于本論文約定采用太陽時(shí),上午下午對(duì)稱分布,積分只需要從日出時(shí)角積分到正午(正午時(shí)角為0°),然后乘2即可。公式為:
如果采用的水平面輻照度是小時(shí)平均值,則積分步長為15°(或1 h);如果采用的是分鐘級(jí)的輻照度,則積分步長需要相應(yīng)變更。值得注意的是,積分步長應(yīng)當(dāng)與太陽跟蹤器的動(dòng)作步長相一致。
將式(37)~式(39)的結(jié)果代入式(8),可得到日總輻射量QT。
各月和全年的方陣面總輻射量,可以先計(jì)算出一年365天的日總輻射量,然后將結(jié)果相加即可。
近年來,雙面光伏組件的應(yīng)用越來越普遍,在計(jì)算雙面組件背面發(fā)電量時(shí),需要計(jì)算出背表面所接收到的輻射量,而背表面輻射量的計(jì)算仍然需要知道任意時(shí)刻背表面的輻照度,其計(jì)算需要用到同樣的5個(gè)變量:QH′、DH′、SH′、θ和Z′,計(jì)算公式與正表面輻射量的計(jì)算公式(見式(13))類似。
背表面輻射量QBT′ 的計(jì)算公式為:
式中,QBT′、SBT′、DBT′、RBT′分別為傾斜方陣背表面的總輻照度、直接輻照度、散射輻照度和反射輻照度。
說明:
1)春分之后,日出太陽方位角會(huì)大于90?,對(duì)于固定方陣來說,在早晚會(huì)有一段時(shí)間的陽光是照在光伏陣列背面的,所以,背表面是會(huì)接收到直接輻射的;
2)由于方陣正面和背面的法線是同一根,所以背面直接輻照度的計(jì)算公式與正面的一致;
3)對(duì)于散射輻照度和反射輻照度,背面輻照度計(jì)算時(shí)只需要將正面計(jì)算公式括號(hào)中的加號(hào)、減號(hào)對(duì)調(diào)即可;
4)采用雙面組件時(shí),根據(jù)經(jīng)驗(yàn),組件的安裝高度至少距離地面1 m;
5)背面屬于低輻照度,不考慮散射輻射的各向異性;
6)計(jì)算暫時(shí)沒有考慮不均勻度的影響。
本文從太陽和地球的運(yùn)行規(guī)律和相互關(guān)系出發(fā),運(yùn)用球面三角的幾何原理,系統(tǒng)地推導(dǎo)出了不同坐標(biāo)系中各種光伏方陣面上實(shí)時(shí)太陽輻照度的計(jì)算方法和公式,為光伏發(fā)電系統(tǒng)的設(shè)計(jì)提供了計(jì)算理論和工具。方陣面上實(shí)時(shí)太陽輻照度可用于以下幾方面:
1)為各類太陽跟蹤器跟蹤太陽的原理提供了理論依據(jù);
2)用于計(jì)算各類不同光伏系統(tǒng)在任何工況條件下日、月、年方陣面(包括背表面)的輻射量,進(jìn)一步可以計(jì)算光伏系統(tǒng)的發(fā)電量;
3)可以得到太陽輻照度在光伏方陣面上的日分布曲線,在知道負(fù)載日負(fù)荷曲線的情況下,可以預(yù)測(cè)“自發(fā)自用系統(tǒng)”自發(fā)自用電量的比例和項(xiàng)目的經(jīng)濟(jì)性;
4)當(dāng)光伏系統(tǒng)在微電網(wǎng)中需要削峰填谷或者錯(cuò)峰運(yùn)行時(shí),方陣面太陽輻照度的日分布曲線可以用于測(cè)算所需配備蓄電池的容量;
5)當(dāng)需要優(yōu)化光伏-逆變器容配比時(shí),需要預(yù)測(cè)在特定光伏-逆變器容配比條件下的限功率運(yùn)行所損失的發(fā)電量,方陣面上實(shí)時(shí)太陽輻照度日分布曲線可以測(cè)算限功率運(yùn)行的時(shí)段,從而計(jì)算出限功率運(yùn)行所損失的發(fā)電量,進(jìn)而優(yōu)化設(shè)計(jì)光伏-逆變器容配比;
6)在將來與其他類型發(fā)電單元、儲(chǔ)能系統(tǒng)組成的混合發(fā)電系統(tǒng)或微電網(wǎng)時(shí),準(zhǔn)確的光伏系統(tǒng)發(fā)電曲線將為各電源、負(fù)載和儲(chǔ)能系統(tǒng)間的平衡、優(yōu)化提供設(shè)計(jì)依據(jù)。
[1]Liu B Y H, Jordan R C. A Rational Procedure for Predicting The Long-Term Average Performance of Flat-Plate Solar-Energy Collectors[J]. Solar Energy, 1963, 7: 53-74.
[2]Klein S A. Calculation of Monthly Average Insolation on Tilted Surfaces[J]. Solar Energy, 1977, 19(4): 325-329.
[3]Braun J E, Mitchell J C. Solar Geometry for Fixed and Tracking Surfaces[J]. Solar Energy, 1983, 31(5): 439-444.
[4]王斯成. 光伏電源系統(tǒng)的計(jì)算機(jī)輔助設(shè)計(jì)[J]. 太陽能學(xué)報(bào),1986, 7(3): 251-264.
[5]張楚斌. 球面三角學(xué)[M]. 北京: 人民教育出版社, 1959.
[6]Hay J E. Calculation of Monthly Mean Solar Radiation for Horizontal and Inclined Surfaces[J]. Solar Energy, 1979, 23(4): 301-307.
[7]Klucher T M. Evaluation of Models to Predict Insolation on Tilted Surfaces[J]. Solar Energy, 1979, 23(2): 111-114.
[8]Perez R, Stewart R, Arbogast C, et al. An Anisotropic Hourly Diffuse Radiation Model for Sloping Surfaces: Description,performance validation, site dependency evaluation [J]. Solar Energy, 1986, 36(6): 481-497.
[9]Perez R, Seals R, Ineichen P, et al. A New Simplif i ed Version of the Perez Diffuse Irradiance Model for Tilted Surfaces[J]. Solar Energy, 1987, 39(3): 221-231.
[10]楊金煥. 固定式光伏方陣最佳傾角的分析[J]. 太陽能學(xué)報(bào),1992, 13(1): 86-92.
[11]楊金煥. 不同方位傾斜面上太陽輻射量及最佳傾角的計(jì)算[J]. 上海交通大學(xué)學(xué)報(bào), 2002, 36(7): 1032-1036.
[12]Klein S A, Theilacker J C. An Algorithm for Calculating Monthly-Average Radiation on Inclined Surfaces[J]. Journal Solar Energy Engineering, 1981, 103(1): 29-33.
[13]Natural Resources Canada (NASA, UNEP, GEF).RETSCREEN?ENGINEERING & CASES TEXTBOOK[EB/OL]. http://www.retscreen.net, 2004.
[14]Duffie J A, Beckman W Z. Solar Engineering of Thermal Process(3rdedition) [M]. John Wiley & Sons, 2005.
[15]日本太陽能學(xué)會(huì)[編]. 太陽能的基礎(chǔ)和應(yīng)用[M]. 上海:上??茖W(xué)技術(shù)出版社, 1982.