曹丁元，金秀章，孫小林，謝澤坤

(華北電力大學(xué) 控制與計算機(jī)工程學(xué)院，河北 保定071003)

0 引言

近些年環(huán)境污染、資源短缺一直是人們討論的熱門話題，具有無污染、可再生、經(jīng)濟(jì)、清潔等優(yōu)點(diǎn)的太陽能成為專家學(xué)者的研究對象。西安交通大學(xué)的魏進(jìn)家教授提出了復(fù)合拋物面聚光器(CPC)型太陽能光伏/光熱聯(lián)合利用(PV/T)系統(tǒng)［1］。采用CPC 聚光器的光電和光熱相結(jié)合的PV/T 系統(tǒng)，當(dāng)聚光比確定后，影響系統(tǒng)性能的主要參數(shù)是太陽電池的工作溫度。因此，確定太陽電池的工作溫度對于系統(tǒng)設(shè)計具有指導(dǎo)意義。目前為止，國內(nèi)外很多學(xué)者對CPC-PV/T 系統(tǒng)的性能分析進(jìn)行了不少研究［2～6］，但是涉及CPC-PV/T 綜合系統(tǒng)中太陽電池溫度的建模報道尚少。本文將對系統(tǒng)中太陽電池工作溫度的機(jī)理模型建立進(jìn)行討論。機(jī)理建模［7］是在充分了解了受控對象、執(zhí)行機(jī)構(gòu)及系統(tǒng)內(nèi)一切元件的運(yùn)動規(guī)律，將這些運(yùn)動規(guī)律用數(shù)學(xué)模型表示出來的過程。機(jī)理建模一直以來廣泛應(yīng)用于各個領(lǐng)域［8，9］，具有理論依據(jù)嚴(yán)密，在任何狀態(tài)下使用都不會引起定性錯誤的優(yōu)點(diǎn)。

1 數(shù)學(xué)模型分析

CPC-PV/T 系統(tǒng)簡圖如圖1 所示，它主要由兩部分構(gòu)成:太陽能光電系統(tǒng)以及太陽能換熱系統(tǒng)。太陽能光電系統(tǒng)是PV 電池板接收到來自聚光器的光照而產(chǎn)生電能的系統(tǒng);太陽能換熱系統(tǒng)是下層流道中的流體吸收PV 電池板的廢熱，從而使電池板溫度處在最佳工作點(diǎn)，同時產(chǎn)生可利用的熱能的系統(tǒng)。本文主要對太陽能換熱系統(tǒng)的模型進(jìn)行分析。

圖1 CPC-PV/T 系統(tǒng)簡圖

電池板的溫度主要受太陽輻射和下層流體流量的影響。太陽輻射照射到光伏電池上，會使電池板的溫度升高，從而影響到系統(tǒng)的光電效率。為了降低電池板的溫度，使其處在最佳的工作溫度下，本課題采用了平板PV/T 集熱器，當(dāng)下層流體流過電池板，帶走電池板的廢熱，即降低了電池板的溫度，又提高了流體的溫度，實現(xiàn)了熱量的充分利用。

對于下層流體與電池板的熱交換部分，建模方法分為集總參數(shù)建模法［10，11］和分布參數(shù)建模法［11，12］。在流體與電池板導(dǎo)熱過程中，流體與電池板的狀態(tài)參數(shù)不僅是時間而且還是空間的函數(shù)，具有明顯的分布參數(shù)特點(diǎn)。這種分布參數(shù)的動態(tài)模型都是非線性偏微分方程組，這些偏微分方程組含有長度和時間兩個自變量。這樣一個分布參數(shù)系統(tǒng)的非線性模型是十分復(fù)雜的，計算也有一定難度。為了簡化問題，在建立數(shù)學(xué)模型時，將下層流道內(nèi)的介質(zhì)狀態(tài)參數(shù)看成是均勻一致的，并在空間位置上選定一個有代表性的點(diǎn)，利用這一點(diǎn)介質(zhì)的參數(shù)作為環(huán)節(jié)的集總參數(shù)。在這種簡化假定下建立的模型就是集總參數(shù)模型。

為了便于分析，需要對換熱管道作幾個合理的簡化假設(shè):

(1)電池板無軸向?qū)?

(2)下層流道內(nèi)液相介質(zhì)為不可壓縮流體;

(3)忽略流體的軸向傳熱;

(4)電池板與下層流道的換熱系數(shù)為常數(shù);

(5)同一橫截面的流體的參數(shù)相同;

(6)電池板沿下層流道長度方向的放熱量均勻。

1.1 分布參數(shù)建模

電池板與下層流道簡圖如圖2 所示，在上述簡化假設(shè)條件下，可列出距入口x 處的微分流道段dx 中在單位時間內(nèi)的熱平衡方程如下:

圖2 電池板與下層流道簡圖

下層水熱平衡方程為:

式中:ρw，cw，t1分別為下層水的密度、比熱容、入口溫度;tc為電池板的溫度;Aw為流體流道換熱面積;f 為流體流道內(nèi)橫截面積;h 為電池板導(dǎo)熱系數(shù);L 為電池板長度;D 為下層流體的質(zhì)量流量。上式整理后得

對于電池板，其能量平衡方程為:

式中:mc，cc分別為電池板的質(zhì)量和比熱容。

電池板吸收的太陽輻射的熱量:

式中:ηth為光伏電池光熱轉(zhuǎn)換效率;S 為玻璃蓋板對太陽光的透射率;ηc為分頻玻璃的分頻效率;k 為聚光比;Ac為電池板的面積;G 為太陽輻射強(qiáng)度，可測量。

將(2)式代入上式，整理后得

電池板的溫度變化的原因可以來自于下層流道流體的流量(流體的初始溫度為恒定值)或者光照強(qiáng)度。在不同的擾動下，電池板溫度的動態(tài)特性是不同的。文獻(xiàn)［13～15］等分別針對不同的擾動因素，通過拉普拉斯(Laplace)變換，求出方程組的解析解。本文對不同擾動同時發(fā)生求出解析解。

對式(1)、(3)兩式化為增量形式，然后進(jìn)行二元拉氏變換，消去中間變量Δt1(x，s)得:

對(4)式進(jìn)行求解得到電池板初始溫度、流量和光照強(qiáng)度3 項擾動同時存在時的疊加響應(yīng)為:

式中:tc0(s)為電池板初始溫度。

為了得到各個擾動對系統(tǒng)的獨(dú)立作用，分別令其他兩個擾動為零。對于初始溫度擾動，令D(s)=0，G(s)=0，最后得到其傳遞函數(shù)為

同理，可以得到太陽光照強(qiáng)度ΔG(s)和水工質(zhì)流量ΔD(s)擾動下的傳遞函數(shù)分別為

由于式(6)～(8)均為超越型傳遞函數(shù)。計算起來很不方便，因此，采用低階逼近將它們進(jìn)行簡化［2，3］。具體方法是將上述傳遞函數(shù)在s=0 附近的領(lǐng)域內(nèi)展開為泰勒級數(shù)的形式，并與假設(shè)的傳遞函數(shù)相應(yīng)的展開式進(jìn)行比較，從而確定假設(shè)的傳遞函數(shù)的各項數(shù)。

假設(shè)各個擾動下的傳遞函數(shù)為

2017年，IPO審核加速帶動了創(chuàng)投基金的募集、投資和退出，政府引導(dǎo)基金的進(jìn)一步擴(kuò)容引導(dǎo)了社會資金更多地流進(jìn)創(chuàng)投基金和實體經(jīng)濟(jì)，港股、美股的窗口期為具有國際視野的創(chuàng)投基金帶來更大的機(jī)會，整個中國創(chuàng)投行業(yè)呈現(xiàn)出良好的發(fā)展態(tài)勢。無疑，創(chuàng)投行業(yè)的快速發(fā)展，已經(jīng)成為破解中小企業(yè)“融資難”問題、加快金融供給側(cè)改革的重要手段。

對于電池板初始溫度擾動，將式(6)與式(9)均在s=0 點(diǎn)處展開泰勒級數(shù)，取其前三項相等得到式(9)的系數(shù)分別為:，T。指數(shù)n 總是四舍五入取整數(shù)。

對于太陽光照強(qiáng)度G 擾動，將式(7)與式(9)均在s=0 點(diǎn)處展開泰勒級數(shù)，取其前三項相等得到(9)式的系數(shù)的系數(shù)分別為:K=。指數(shù)n 總是四舍五入取整數(shù)。

對于水工質(zhì)流量d 擾動，將式(8)與式(9)均在s=0 點(diǎn)處展開泰勒級數(shù)，取其前三項相等得到式(9)的系數(shù)分別為:，n=2，T。指數(shù)n 總是四舍五入取整數(shù)。

1.2 集總參數(shù)建模法

水的能量平衡:

式中:Qi為水與電池板之間的換熱量;t1，t2分別為水的進(jìn)出口溫度。

電池板對水的放熱量:

式中:tc為電池板的溫度。

電池板的能量平衡方程:

由式(2)、(11)、(12)可得:

對上述式(13)、(14)分別作拉氏變換，經(jīng)推導(dǎo)得溫度、流量和光照強(qiáng)度擾動的傳遞函數(shù)模型:

電池板初始溫度擾動:

水工質(zhì)流量D 擾動:

太陽光照強(qiáng)度G 擾動:

2 模型仿真與分析

根據(jù)表1 中的原始數(shù)據(jù)和簡化結(jié)果得到分布參數(shù)模型在各個擾動下的傳遞函數(shù)為:

表1 參數(shù)物理意義

電池板初始溫度擾動:

太陽光照強(qiáng)度擾動:

水工質(zhì)流量擾動:

根據(jù)表1 中的原始數(shù)據(jù)和式(15)～(17)得到集總參數(shù)模型在各個擾動下的傳遞函數(shù)為:

電池板初始溫度擾動:

太陽光照強(qiáng)度擾動:

水工質(zhì)流量擾動:

求取以上各傳遞函數(shù)的階躍響應(yīng)曲線，如圖3，4，5 所示，(a)圖均為分布參數(shù)模型，(b)圖均為集總參數(shù)模型。由圖3 可以看出，當(dāng)光照強(qiáng)度增加1 W/m2，電池板的溫度大約增加0．024 K;由圖4 可以看出，當(dāng)流體流量增加1 kg/s，其他擾動為零時，電池板溫度下降大約1．35 K;由圖5 可以看出，當(dāng)電池板初始溫度增加1 K，其他擾動為零時，電池板最終溫度在一定的時間范圍內(nèi)上升1 K。綜上，當(dāng)光照強(qiáng)度增加時電池板的溫度會升高。當(dāng)流體的質(zhì)量流量增加時，會帶走更多的熱量，從而電池板的溫度會有所下降，由此可判斷流量越快，電池板降溫越快。當(dāng)電池板的初始溫度發(fā)生變化的時候電池板的最終溫度會變化相同的幅度。對比分布參數(shù)模型和集總參數(shù)模型的曲線可以看出分布參數(shù)模型能更好地反映系統(tǒng)的動態(tài)性能，但是分布參數(shù)法計算復(fù)雜。

圖3 電池板溫度在光照強(qiáng)度擾動下的響應(yīng)

圖4 電池板溫度在流體流量擾動下的響應(yīng)

圖5 電池板溫度在初始溫度擾動下的響應(yīng)

3 結(jié)論

通過對復(fù)合拋物面聚光器(CPC)型太陽能光伏/光熱聯(lián)合利用(PV/T)系統(tǒng)生產(chǎn)工藝的研究，本文利用能量守恒定律得出了PV/T 系統(tǒng)的分布參數(shù)模型和集總參數(shù)模型，分別給出了在不同擾動下電池板溫度的傳遞函數(shù)，并且進(jìn)行了仿真。仿真結(jié)果表明，這兩種機(jī)理建模方法都能夠很好地反映系統(tǒng)的動態(tài)性能。由模型可以看出，流體流量是控制電池板溫度在最佳工作點(diǎn)的主要因素。機(jī)理模型的確定為CPC-PV/T 系統(tǒng)的進(jìn)一步優(yōu)化奠定了理論依據(jù)，為CPC-PV/T 系統(tǒng)控制算法的研究提供了技術(shù)支撐。

［1］魏進(jìn)家，謝胡凌，敬登偉．一種太陽能聚光分頻光伏光熱聯(lián)產(chǎn)裝置［P］．中國專利．CN201210224815．9，2012．

［2］劉亞雷，張紅，許輝，等．CPC 型聚光光伏光熱系統(tǒng)的性能分析［J］．可再生能源，2011，29(1):1-5．

［3］崔文智，于松強(qiáng)，廖全．聚光型混合光伏光熱系統(tǒng)熱電性能分析［J］．重慶大學(xué)學(xué)報，2009，32(1):86-90．

［4］孫健，王艷香，施明恒．復(fù)合拋物面聚光太陽能PV/T 系統(tǒng)的實驗研究［J］．太陽能學(xué)報，2012，33(1):86-91．

［5］Amrizal N，Chemisana D，Rosell J I．Hybrid photovoltaic-thermal solar collectors dynamic modeling［J］．Applied Energy，2013，101:797-807．

［6］Mohd Yusof Hj Othman，Baharudin Yatim，Kamaruzzaman Sopian，et al．Performance analysis of a double-pass photovoltaic/thermal(PV/T)solar collector with CPC and fins［J］．Renewable Energy，2005，30(13):2005-2017．

［7］韓璞，董澤，王東風(fēng)，等．智能控制理論及應(yīng)用［M］．北京:中國電力出版社，2012．

［8］李雅哲，塔式太陽能熱發(fā)電蒸汽系統(tǒng)建模與控制［D］．北京:華北電力大學(xué)，2011．

［9］馬進(jìn)，劉長良，李淑娜．穩(wěn)壓器壓力水位控制系統(tǒng)建模與仿真［J］．核科學(xué)與工程，2010，30(1):9-14．

［10］康英偉，薛陽，黃偉．電站鍋爐過熱器的集總參數(shù)動態(tài)建模與仿真［J］．計算機(jī)仿真，2012，29(9):332-334．

［11］阮剛，羅自學(xué)，周懷春．鍋爐過熱蒸汽溫度控制新策略動態(tài)仿真研究［J］．華中電力，2004，17(5):1-4．

［12］周少祥，胡三高，宋之平．管式換熱器分布參數(shù)模型的分段線性化方法研究［J］．中國電機(jī)工程學(xué)報，2002，22(6):123-125．

［13］李旭．蒸汽發(fā)生器的動態(tài)特性——模型、結(jié)構(gòu)參數(shù)的影響和控制［D］．北京:清華大學(xué)，1988．

［14］倪維斗，徐向東，李政，等．熱動力系統(tǒng)建模與控制的若干問題［M］．北京:科學(xué)出版社，1996．

［15］(美)E0．道別林，童鈞芳．系統(tǒng)的建模和響應(yīng)—理論及實驗方法［M］．關(guān)國樞譯．上海:科學(xué)技術(shù)文獻(xiàn)出版社，1984．