宋景慧，　馬繼帥，　湛志鋼，　代彥軍

(1.廣東電網(wǎng)有限責(zé)任公司電力科學(xué)研究院，廣州 510080；2.上海交通大學(xué) 制冷與低溫工程研究所，上海 200240)

?

線性菲涅爾集熱器光學(xué)特性實(shí)例分析與模擬

宋景慧1，馬繼帥2，湛志鋼1，代彥軍2

(1.廣東電網(wǎng)有限責(zé)任公司電力科學(xué)研究院，廣州 510080；2.上海交通大學(xué) 制冷與低溫工程研究所，上海 200240)

摘要：對(duì)線性菲涅爾集熱器的聚光性能和光學(xué)效率進(jìn)行了模擬與計(jì)算.根據(jù)幾何光學(xué)原理，對(duì)集熱器鏡場(chǎng)各項(xiàng)光學(xué)損失(如余弦損失、陰影與遮擋損失)建立數(shù)學(xué)模型，計(jì)算出每項(xiàng)光學(xué)損失對(duì)應(yīng)的光學(xué)效率.再用TracePro光學(xué)軟件建立集熱器的幾何模型，利用光線追蹤的方法，模擬入射到鏡場(chǎng)的光線在模型空間的傳播.光線在模型表面發(fā)生吸收、反射和散射等過(guò)程，追蹤每束光線的光通量，計(jì)算得到集熱器的光學(xué)效率和聚光比等性能參數(shù).結(jié)果表明：通過(guò)數(shù)學(xué)模型和光學(xué)軟件模擬得出的集熱器光學(xué)效率一致，2種方法分別從細(xì)節(jié)與整體上剖析了影響集熱器光學(xué)性能的因素，在集熱器設(shè)計(jì)中可以結(jié)合使用，互相補(bǔ)充.

關(guān)鍵詞：線性菲涅爾集熱器； 鏡場(chǎng)參數(shù)； 光學(xué)效率； 能流密度分布； 聚光比

在太陽(yáng)能中溫?zé)崂妙I(lǐng)域，大多使用槽式集熱器和線性菲涅爾集熱器對(duì)太陽(yáng)輻射能量進(jìn)行采集.后者因其易生產(chǎn)、易安裝和易運(yùn)行維護(hù)等優(yōu)點(diǎn)，應(yīng)用越來(lái)越廣泛[1].線性菲涅爾集熱器一般由鏡場(chǎng)(Linear Fresnel Reflector, LFR)、吸收器、跟蹤裝置和支架組成.其中，鏡場(chǎng)由一系列離散的條狀平面鏡組成，鏡面均跟蹤太陽(yáng)位置，將光線準(zhǔn)確反射到吸收器，實(shí)現(xiàn)聚光.

鏡場(chǎng)的光學(xué)特性可以反映集熱器整體的光熱轉(zhuǎn)化性能，光學(xué)效率則集中體現(xiàn)了集熱器的光學(xué)性能.在進(jìn)行鏡場(chǎng)光學(xué)特性分析時(shí)，通常有2種方法：一是建立數(shù)學(xué)理論模型，計(jì)算出給定集熱器的光學(xué)效率，并分析影響光學(xué)效率的每一個(gè)因素，如各項(xiàng)光學(xué)損失；二是利用軟件建立幾何模型進(jìn)行光學(xué)模擬，可以得到集熱器的吸收器吸收面上的能流密度分布，進(jìn)而計(jì)算出光學(xué)效率和聚光比等性能參數(shù).

吸收器吸收面的能流密度一般不是均勻分布的，會(huì)沿著某一坐標(biāo)方向發(fā)生較大變化，比如集熱器吸熱管表面的熱流密度由于遮蔽效應(yīng)，會(huì)分為遮蔽效應(yīng)區(qū)、熱流密度遞增區(qū)和熱流密度衰減區(qū)等[2].合理的熱流密度分布有助于提高集熱器的光熱轉(zhuǎn)化效率，并且能指導(dǎo)吸收器的優(yōu)化設(shè)計(jì).

筆者首先建立線性菲涅爾集熱器鏡場(chǎng)的數(shù)學(xué)理論模型，對(duì)影響集熱器光學(xué)效率的因素進(jìn)行分析，重點(diǎn)介紹余弦損失、末端損失及陰影與遮擋損失的數(shù)學(xué)模型，并計(jì)算出光學(xué)效率.然后利用TracePro光學(xué)模擬軟件，建立不同太陽(yáng)光線入射傾角下的線性菲涅爾集熱器的幾何模型，模擬得出集熱器的能流密度分布，再計(jì)算出光學(xué)效率.最后，將2種方法計(jì)算得到的結(jié)果進(jìn)行對(duì)比分析.

1理論分析

1.1集熱器參數(shù)

本實(shí)例的線性菲涅爾集熱器的腔體吸收器為一半圓倒扣槽式腔體吸收器，長(zhǎng)度方向沿南北子午線.鏡場(chǎng)由12塊反射鏡構(gòu)成，鏡面寬度相同，且鏡面均勻排布于腔體吸收器正下方，吸收器鉛垂面每側(cè)鏡面數(shù)N=6，如圖1所示.鏡場(chǎng)設(shè)計(jì)參數(shù)如表1所示.

為保證經(jīng)鏡場(chǎng)反射到吸收器的光斑寬度小于吸收器開(kāi)口寬度，反射鏡采用內(nèi)凹的微弧度柱面鏡[3]，不同位置的鏡面弧度有細(xì)微的不同.從左到右(從西向東)將反射鏡依次編號(hào)為-6,-5,…,-1,1,…,5,6，由于鏡場(chǎng)鏡面位置的對(duì)稱性，表2只給出鏡場(chǎng)右邊一側(cè)反射鏡的位置與結(jié)構(gòu)參數(shù).表2中，Qn為第n塊反射鏡到鏡場(chǎng)中心線的距離；rn為第n塊反射鏡的弧度圓半徑；dn為第n塊反射鏡的弧高.

圖1　線性菲涅爾集熱器示意圖

表1　集熱器鏡場(chǎng)設(shè)計(jì)參數(shù)

表2　反射鏡位置與結(jié)構(gòu)參數(shù)

集熱器鏡場(chǎng)所在平面與水平面平行，鏡場(chǎng)的太陽(yáng)光線入射角可通過(guò)太陽(yáng)高度角和太陽(yáng)方位角來(lái)描述，如圖2所示.太陽(yáng)入射光線與水平面的夾角稱為太陽(yáng)高度角，記為αs；太陽(yáng)入射光線在水平面的投影與水平面南北水平線之間的夾角稱為太陽(yáng)方位角，記為γs.

圖2　太陽(yáng)入射光線分解示意圖

如圖3所示，太陽(yáng)入射光線在東西立面的投影與向西的水平射線的夾角為α，稱為太陽(yáng)光線入射傾角.鏡面反射光線與水平向西的射線的夾角記為γ；鏡面相對(duì)于水平位置轉(zhuǎn)過(guò)的角度稱為鏡面傾角β，逆時(shí)針為正；鏡面光線入射角記為θ，入射光線在法線左側(cè)為正，右側(cè)為負(fù).

圖3　鏡場(chǎng)光學(xué)幾何參數(shù)

鏡面傾角β、鏡面光線入射角θ與太陽(yáng)光線入射傾角α的關(guān)系如下：

(1)

(2)

表3給出了太陽(yáng)光線入射傾角α分別為15°、30°、45°、60°、75°和90°時(shí)，各個(gè)鏡面傾角β的分布.

表3　不同太陽(yáng)光線入射傾角下鏡面傾角的分布

從表3可以看出，在一定的太陽(yáng)光線入射傾角α下，鏡面傾角從西向東依次增大，角度變化均接近線性；隨著α增大，每個(gè)鏡面朝向均向右旋轉(zhuǎn)；α>60°時(shí)，西側(cè)部分鏡面傾角小于0°，表示這部分鏡面在跟蹤裝置的帶動(dòng)下已朝向東側(cè)；α=90°時(shí)，光線垂直入射鏡場(chǎng)平面，鏡面傾角東西兩側(cè)對(duì)稱分布.α>90°時(shí)，鏡場(chǎng)光學(xué)特性與α<90°時(shí)成對(duì)稱關(guān)系，故只需分析太陽(yáng)光線入射傾角0<α<90°的情況.

1.2光學(xué)效率

光學(xué)效率是聚焦型集熱器一個(gè)非常重要的參數(shù)，其定義[4]如下：

(3)

式中：Qu為腔體吸收器吸收的太陽(yáng)輻射能量，W；A為反射鏡總面積，m2；Ib為太陽(yáng)直射輻射強(qiáng)度，W/m2.

影響集熱器光學(xué)效率的因素眾多，主要包括余弦損失、陰影與遮擋損失、末端損失、吸收器與支架陰影損失及集熱器各表面光學(xué)性能等.每一項(xiàng)光學(xué)損失均對(duì)應(yīng)著一個(gè)效率，集熱器的光學(xué)效率即為各項(xiàng)效率的乘積.

單純考慮鏡場(chǎng)余弦損失的光學(xué)效率ηθ為：

(4)

式中：Ap為鏡場(chǎng)等效面積，m2.

圖4為相鄰鏡面的陰影與遮擋示意圖，其中，Ci、Di和C表示水平距離.單純考慮鏡面間陰影與遮擋損失的光學(xué)效率ηE為：

(5)

單純考慮吸收器末端損失的光學(xué)效率ηend為：

(6)

單純考慮吸收器與支架陰影損失的光學(xué)效率ηsupport為：

(7)

式中：Areceiver和Asupport分別為吸收器與支架的陰影面積，m2.

最后，可得線性菲涅爾集熱器光學(xué)效率的表達(dá)式[5]為：

(8)

式中：ρm為反射材料的反射率，取0.935；τr為吸收器外透明覆蓋物的透過(guò)率，取0.91；αr為吸收腔體對(duì)太陽(yáng)光的吸收率，取0.9；ηtrace為考慮跟蹤誤差所產(chǎn)生的輻射損失后的效率.

(a) 陰影

(b) 遮擋

1.3計(jì)算結(jié)果及分析

本實(shí)例中線性菲涅爾集熱器的安裝地點(diǎn)在廣州市，屬低緯度地區(qū)，所使用的腔體吸收器的長(zhǎng)度大于反射鏡長(zhǎng)度，能減小或消除末端損失.在計(jì)算分析時(shí)，可認(rèn)為末端損失為零，同時(shí)假定鏡面跟蹤精確，即跟蹤誤差為零.

當(dāng)太陽(yáng)光線入射傾角α為15°、30°、45°、60°、75°和90°時(shí)，分別計(jì)算集熱器的光學(xué)效率，結(jié)果如表4.

表4所示為考慮了鏡場(chǎng)余弦效應(yīng)、陰影與遮擋、支架與吸收器陰影及材料光學(xué)性能等因素后，集熱器對(duì)應(yīng)的光學(xué)效率隨著太陽(yáng)光線入射傾角α變化而變化的情況.α增大時(shí)，ηθ增大，表明鏡場(chǎng)余弦損失減小，減小趨勢(shì)較為平緩；ηE增大，表明鏡面間陰影與遮擋損失快速減小，α>45°時(shí)已經(jīng)為零；ηsupport先維持在一個(gè)較大的值，然后減小，再緩慢增大，這是因?yàn)棣?45°時(shí)，吸收器的陰影落在鏡場(chǎng)之外，未造成光學(xué)損失，α>45°時(shí)吸收器陰影落在鏡場(chǎng)內(nèi)，陰影面積減小，但變化很小，所以造成的光學(xué)損失基本維持不變.

整體而言，余弦效應(yīng)和鏡面間的陰影與遮擋在太陽(yáng)光線入射傾角α<45°時(shí)，對(duì)集熱器光學(xué)效率有很大影響，后者尤甚.鏡場(chǎng)余弦損失始終存在，但鏡面間陰影與遮擋造成的損失可以避免.

表4　不同太陽(yáng)光線入射傾角下集熱器的光學(xué)效率

由表4可知，太陽(yáng)光線入射傾角α增大，集熱器的光學(xué)效率ηoptical隨之增大，α=90°時(shí)達(dá)到最大值66.5%，該值為集熱器的純光學(xué)效率.

2光學(xué)軟件模擬

2.1TracePro簡(jiǎn)介

光學(xué)模擬軟件TracePro運(yùn)用“普適光線追跡”技術(shù)來(lái)追跡光線和光通量，用離散分布的光線試樣傳播來(lái)代替連續(xù)的光線分布傳播.在建立幾何模型時(shí),運(yùn)用MonteCarlo算法模擬光線的反射、散射、折射和衍射，其中散射和衍射被認(rèn)為是一種隨機(jī)過(guò)程.在幾何模型物理表面的每個(gè)點(diǎn)上，光線都遵從以上定律.TracePro可以考慮復(fù)雜的表面特性，如鏡面反射、各向異性反射和各向異性折射等，只需針對(duì)不同的截面輻射特性建立相應(yīng)的物理模型，同時(shí)可以考慮到模型之間重疊、遮擋和交叉等問(wèn)題.

2.2模型建立

首先，根據(jù)反射鏡設(shè)計(jì)參數(shù)、腔體吸收器的尺寸以及不同太陽(yáng)光線入射傾角下各鏡面傾角，在TracePro中建立鏡場(chǎng)和吸收器的幾何模型，如圖5(a)所示.定義好各表面屬性及光源，進(jìn)行光線追蹤模擬，可得到玻璃蓋板和腔體半圓形內(nèi)壁上的輻射能量分布，如圖5(b)所示.

2.3模擬結(jié)果

經(jīng)鏡場(chǎng)反射后的太陽(yáng)光線匯聚到吸收器下方的玻璃蓋板上，沿吸收器長(zhǎng)度方向形成一條高能量密度的光斑.聚焦后的光線被玻璃蓋板吸收或反射一小部分后，穿過(guò)玻璃，較為均勻地分散到半圓形腔體內(nèi)壁.

(a)

(b)

圖5線性菲涅爾集熱器幾何模型及光線追蹤模擬(α=45°)

Fig.5Geometry model of the collector and the ray tracing

2.3.1能流密度分布

在腔體吸收器中部取垂直于長(zhǎng)度方向的截面，在不同太陽(yáng)光線入射傾角的情況下，其上能流密度分布極坐標(biāo)圖和直角坐標(biāo)圖分別如圖6和圖7所示.

(a)α=45°

(b)α=90°

從圖6和圖7可以看出半圓形腔體吸收器內(nèi)壁能流密度分布的特點(diǎn)：(a)隨著太陽(yáng)光線入射傾角α的增大，吸收器內(nèi)壁的能流密度整體上增大；(b)α>45°時(shí)，在圓心角為180°-α的地方，吸收器內(nèi)壁的能流密度分布出現(xiàn)一個(gè)低谷，這是由腔體吸收器在鏡場(chǎng)的投影造成的；(c)不論α多大，吸收器內(nèi)壁的能量集中分布在30°～150°的圓心角范圍內(nèi).因此，吸收器吸熱管應(yīng)布置在這個(gè)角度范圍內(nèi).

圖7　吸收器內(nèi)壁能流密度分布直角坐標(biāo)圖

2.3.2玻璃蓋板光斑寬度

腔體吸收器玻璃蓋板上光斑寬度體現(xiàn)的是集熱器鏡場(chǎng)的聚光性能.圖 8為入射到玻璃蓋板的輻射能量圖.

圖8　入射到吸收器玻璃蓋板的輻射能量

從圖8可以看出，太陽(yáng)光線入射傾角α從0°變化到90°，玻璃蓋板上光斑寬度逐漸減小，以輻射通量300 W/m2為光斑邊緣，可得到光斑寬度的具體值，如表5所示.根據(jù)得到的光斑寬度，可以計(jì)算出集熱器的幾何聚光比.

2.3.3光學(xué)效率

利用數(shù)值積分，計(jì)算出吸收器內(nèi)壁吸收的輻射能量，據(jù)此可得到集熱器的光學(xué)效率，如表6所示.

從表6可以看出，根據(jù)數(shù)學(xué)模型計(jì)算得出的集熱器光學(xué)效率與根據(jù)TracePro模擬得出的光學(xué)效率變化趨勢(shì)一致，均隨著太陽(yáng)光線入射傾角α增大而迅速增大，后趨平緩.模擬得出的光學(xué)效率略低于計(jì)算得出的值，但最大差值不超過(guò)2.6%.其原因?yàn)椋号c數(shù)學(xué)計(jì)算相比，光學(xué)模擬考慮了光線散射，部分散射光線并沒(méi)有投射到吸收器，形成了額外能量損失.可以認(rèn)為，2種方法得出的光學(xué)效率能很好地吻合.

表5　玻璃蓋板光斑寬度與幾何聚光比

表6　計(jì)算與模擬得出的集熱器光學(xué)效率

注：1)差值=模擬光學(xué)效率-計(jì)算光學(xué)效率.

3結(jié)論

(1)余弦效應(yīng)和鏡面間的陰影與遮擋在太陽(yáng)光線入射傾角α<45°時(shí)，對(duì)集熱器光學(xué)效率有很大影響，后者尤甚.鏡場(chǎng)余弦損失始終存在，但鏡面間陰影與遮擋造成的損失可以避免.

(2)線性菲涅爾集熱器的最大光學(xué)效率為66.5%.

(3)隨著太陽(yáng)光線入射傾角α的增大，吸收器玻璃蓋板上的光斑寬度逐漸減小，幾何聚光比增大，可達(dá)63.2；同時(shí)，吸收器內(nèi)壁的能流密度整體上增大，且能量始終集中分布在30°～150°的圓心角范圍內(nèi)，腔體吸收器的陰影雖會(huì)逐漸造成能量分布低谷帶，但對(duì)集熱器光學(xué)效率的影響較小.

(4)軟件模擬得出的集熱器光學(xué)效率與數(shù)學(xué)模型計(jì)算得出的光學(xué)效率能很好地吻合.用數(shù)學(xué)模型計(jì)算的方法能夠清晰地得到集熱器的各項(xiàng)光學(xué)損失，如余弦損失、陰影與遮擋損失等，以及對(duì)應(yīng)影響因素下的效率.光學(xué)軟件模擬的方法能夠得到鏡場(chǎng)吸收器各個(gè)表面的能流密度分布、光斑寬度與幾何聚光比，可以指導(dǎo)腔體吸收器的設(shè)計(jì).2種方法可分別從細(xì)節(jié)與整體上分析集熱器的光學(xué)性能，在集熱器設(shè)計(jì)中可以結(jié)合使用，互相補(bǔ)充.

參考文獻(xiàn)：

[1]張輝，王一平，朱麗，等. 條形平面鏡聚光器設(shè)計(jì)參數(shù)的分析及優(yōu)化[J]. 太陽(yáng)能學(xué)報(bào)，2013,34(11)：1882-1887.

ZHANG Hui, WANG Yiping, ZHU Li,etal. Analysis and optimization of the design parameters of linear flat mirror concentrator[J]. Acta Energiae Solaris Sinica, 2013,34(11): 1882-1887.

[2]韋彪，朱天宇. DSG太陽(yáng)能槽式集熱器聚光特性模擬[J]. 動(dòng)力工程學(xué)報(bào)，2011,31(10)：773-777.

WEI Biao, ZHU Tianyu. Simulation on concentrating characteristics of DSG catabolic trough collectors[J]. Journal of Chinese Society of Power Engineering, 2011,31(10): 773-777.

[4]DUFFIE J A, BECKMAN W A. Solar engineering of thermal processes[M]. Hoboken, New Jersey，USA:John Wiley & Sons, Inc, 2013.

[5]陳宇. 太陽(yáng)能菲涅爾線聚焦集熱器與中溫蓄熱裝置研究[D]. 上海：上海交通大學(xué)，2012：22-25.

Optical Analysis and Simulation of a Linear Fresnel Solar Collector

SONGJinghui1,MAJishuai2,ZHANZhigang1,DAIYanjun2

(1.Electric Power Research Institute of Guangdong Power Grid Co., Ltd., Guangzhou 510080, China;2. Institute of Refrigeration and Cryogenics, Shanghai Jiaotong University, Shanghai 200240, China)

Abstract：The optical efficiency and concentrating property of a linear Fresnel solar collector were calculated mathematically and simulated optically. Based on the principles of geometrical optics, mathematical models of relative optical losses (cosine loss, shading & blocking loss, etc.) in the collector mirror field were established so as to calculate related optical efficiency corresponding to each energy loss). Meanwhile, geometrical models of the collector were created using software TracePro to simulate the transmission of ray incident to the mirror field in model space by ray tracing method, including the absorption, reflecting and scattering of ray on the model surface, thus the optical efficiency and concentration ratio of the collector could be calculated by tracing the radiation flux of the ray. Results show that the optical efficiency obtained respectively by methamatical and geometrical model agrees well with one another. The detail and integral analysis on factors influencing the optical performance of linear Fresnel solar collectors may serve as a reference for the design of similar collectors, which may be used in combination or be mutually complementary with each other.

Key words：linear Fresnel collector; mirror field parameters; optical efficiency; energy flux distribution; concentration ratio

收稿日期：2015-07-29

修訂日期：2015-09-02

基金項(xiàng)目：南方電網(wǎng)有限責(zé)任公司科技基金資助項(xiàng)目(K-GD2013-0489)；國(guó)家科技支撐計(jì)劃課題資助項(xiàng)目(2012BAA05B04)

作者簡(jiǎn)介：宋景慧(1973-)，男，黑龍江密山人，教授級(jí)高工，碩士，主要從事能源高效清潔利用方面的研究.

文章編號(hào)：1674-7607(2016)07-0563-06中圖分類號(hào)：TK513

文獻(xiàn)標(biāo)志碼：A學(xué)科分類號(hào)：480.60

馬繼帥(通信作者)，男，碩士研究生，電話(Tel.)：13122183879；E-mail：biansaixingyin@163.com.