獨立型太陽能發電系統在負載變動下的動態特性分析

(1)

獨立型太陽能發電系統在負載變動下的動態特性分析

廖立智1_、陳盟仁1_、吳有基2_、林建文3 1_{國立高雄應用科技大學電機工程系} 2_{聯合大學電機系} 3_{金陽機電工程有限公司} E-mail: [email protected]

摘要

本論文主要在探討太陽能發電系統在獨立運轉情況下負載變動時的動態特性。該系統主要包括太陽能發電系統、功率調節系統、三相變壓器、感應電動機、以及靜態負載等。研究方法是先推導系統組件的數學模型，接著利用 SimPowerSystems 模組開發並依系統架構連接，最後進行動態特性模擬與分析。研究結果顯示在此種運轉模式下系統變數的變動情形是可以接受的並且應可和實際運轉情形一致。本研究最大的價值是可做為太陽能系統規劃、運轉、以及系統擴充的重要參考。關鍵詞：太陽能發電系統、動態特性分析、功率調節系統、MATLAB/Simulink、SimPowerSystems

1. 前言

台灣地區自產能源缺乏，99%以上能源需仰賴進口。原油價格日益高漲，連帶天然氣、煤炭價格都受到影響，新能源或綠色再生能源的開發熱潮將進一步掀起。面對能源使用量增加、化石能源枯竭的情況，唯有提高再生能源的使用比例，才能面對嚴峻能源挑戰。根據經濟部能源局的規劃，再生能源發電的總裝置容量，將在2025 年達到全國總裝置容量 15%。要達到這個目標，必須增設再生能源發電容量約 8450MW。台電積極配合政府的「再生能源發展方案」，進行各類型再生能源的發展評估。具發展潛力的小水力、風力、太陽光電、海洋溫差以及波浪發電等低碳能源均列為發展重點。到2008 年底為止，台電再生能源（含慣常水力）裝置容量為2184MW，占系統裝置容量 5.7%，淨尖峰能力為 1326MW，占系統淨尖峰能力 3.4%[1][2]。經濟部於2009 年啟動 5 項節能減碳政策，其中之一是公共建築引進太陽光電及打造台灣最大太陽光電發電廠。台電公司目前已完工運轉中的太陽光電發電系統有樹林綜合研究所、台北市區營業處、大林電廠、南投區處、高訓中心、金門區處、南部展示館、台中電廠、花蓮及台東區處等共 233.5kW。2010 年 2 月在高雄縣永安鹽灘的自有土地上建置4MW 的太陽能發電系統，完成後與旁邊的台電興達發電廠併聯發電[3]。此外，台電也制定「外界提供處所合作設置太陽能發電系統作業要點」以配合其他重要電源開發。實施方式是由公司以外的機構提供適宜場所設置太陽光電系統，擴大太陽光電的普及性。目前較大型的太陽能光電系統有高雄世運會場、台南柳營科技園區、北台光電遊憩城[4]、[5]。太陽能模組的輸出為一個直流源，需經過轉換器提供輸出，因此在太陽能發電系統控制方面更顯得重要。Shataf 等人提出一個新的最大功率追蹤器控制法則，其使用四象限的 PWM 轉換器，藉此獲得良好的最大功率追蹤速度[6]。Abdulaziz 等人使用 ANFIS 模型，提出線性相關性分析實驗數據。該研究使用短路電流和開路電壓作為輸入因素，模糊控制器利用ANFIS 的輸出電壓進行最大功率追蹤因而達到高效率與低漣波[7]。Xiaojin 等人提出一種新型微太陽能系統結構與最大功率追蹤方法。該研究採用 Fibonacci 序列可改

(2)

善日照度不均的情況[8]。Bellini 等人研究零電流轉換器，使太陽能發電系統在傳輸損失上有明顯的降低[9]。 Lee 等人開發一種新的高效率 DC/DC 轉換器，應用於實務上相當實惠[10]。Senjyu 等人提出使用最小階觀測控制方法。根據負載的變化控制太陽能功率得到最小的頻率誤差，進而有效地減少頻率的誤差並且達到最大功率輸出[11]。Hashimoto 等人提出一個使用半橋式整流變頻器加上新的控制電路新的變頻器，進而得到較高的效率以及較低的漣波電壓[12]。Kakimoto 等人利用雙層電容器控制太陽能發電的斜率，使得電容快速吸收太陽光電的漣波，進而改變輸出斜率[13] 本論文主要目的在探討獨立型太陽能發電系統在負載變動情況下的動態特性。研究方法是採用 MATLAB 內建的 Simulink 與架構在同一環境下的 SimPowerSystems 來建構模組，並且模擬負載變動時系統變數變化情形。此研究的價值在於可以作為獨立型太陽能發電系統規劃、運轉、以及擴充的重要參考

。

2. 獨立型太陽能發電系統

2.1 系統架構圖

圖1 為一個獨立型太陽能發電系統的架構，此系統的組件包括一組太陽能模組(Photovoltaic Array, PV)、一套功率調節系統（Power Conditioning System, PCS）、一部三相感應電動機(Induction Motor)、一台三相變壓器(Power Transformer)、以及集總的靜態負載(Lumped Static Load)。

2.2 太陽能模型 太陽能模組是由許多太陽能電池串並聯所組成。太陽能電池是由許多 P-N 接面的半導體組成，能夠直接將光能轉換成電能。如圖 2 所示，太陽能模組的等效電路可以用一個電流源表示，其中電流源 Iph 用來表示太陽能板經由光照射後所產生的電流，Dj 用來表示一個 P-N 接面的二極體，Rs 和 Rsh 則分別表示材料內部的等效串聯和並聯電阻，RL 表示外加負載電阻，Ipv 及 Vpv 則分別表示太陽能板的輸出電流及電壓。一般情況下，Rsh 的值很大，Rs 而的值很小。為了簡化分析可將 Rsh 和 Rs 忽略不計[14]。太陽能電池的輸出電流與輸出電壓的關係可以表示成 mp

I

−

+ ]

1 )

tt a st tt oc mp sc pv

T

t

E

t

E

t

E

V

C

V

C

I

t

I

(

)

=

{

1 −

[exp(

)

−

1 ]}

+

(

)

)[

(

)

+

0 .

002 (

2 1

α

(1)

))

(

02 .

E

_tt

t

0 )

(

)]

)

(

log(

0539

.

0

1 [

)

(

T

t

E

t

E

V

t

V

_a st tt mp pv

=

+

β

+

(2)

1

(1

) e

mp 2

xp[

mp

]

sc oc

I

V

C

I

C V

−

= −

(3)

2

ln(

V

C

1

1 )

mp oc mp sc

V

I

−

=

(4)

(3)

其中Ipv 是太陽能電池輸出電流，Isc 是太陽能電池短路電流，Voc 是太陽能電池開路電壓，Vmp、Imp 是太陽能電池在最大功率點的電壓及電流，Ett 是太陽日照強度，Est 是太陽能電池參考日照強度，α 為太陽能電池電流溫度係數，β 為太陽能電池電壓溫度係數，Ta 是環境溫度，其參數設定參照 Siemens SP75 規格 [15]、[16]。圖 3 是太陽能電池在固定環境溫度下，當日照強度改變時，其輸出電流與輸出電壓的關係圖。由PV 模組的特性曲線得知，PV 模組既非電壓源亦非電流源。因此，與轉換器結合時，PV 模組可當作電壓源也可當作電流源，等效電路如圖4 所示。圖1 獨立型太陽能發電系統架構圖2 太陽能電池等效電路圖

(4)

圖3 環境溫度為 25℃，日照度變動時，太陽能電池輸出電流與輸出電壓關係圖圖4 PV 模組電流源轉換器與電壓源轉換器 2.3 功率調節系統模型 功率調節系統(PCS)的功能之一是將太陽能的直流電壓，轉換成市電的電壓及頻率。一般而言分散式電源與市電連接都須經過PCS 的轉換，如燃料電池、風能發電、太陽能發電及微輪機等。圖 5 顯示 PCS 的系統方塊圖，包含了昇壓轉換器、最大功率追蹤器、以及換流器三部份[17]。

(5)

圖5 功率調節系統方塊圖 2.4 昇壓轉換器模型 圖 6 為昇壓轉換器的架構圖。假設電容器 C 已被充電，則當功率開關 IGBT 導通時，VDS=0，電感電流逐漸增加；此時由於在二極體 D1 陽極之電位會小於輸出電壓 V0，所以，二極體 D1 就會因逆向偏壓而截止，負載電壓由電容釋放能量來維持[10]。昇壓轉換器的電壓方程式可以表示成

1

V

D

=

o

₋

(5)

其中V0 為昇壓後的直流電壓，V1 為昇壓前的直流電壓，D 為轉換器之責任週期。圖6 昇壓轉換器電路圖 2.5 最大功率追蹤器模型 目前最大功率追蹤有許多方法如電壓迴授法、功率迴授法、擾動與觀察法、增量電導法、直線近似法、實際量測法等。而擾動觀察法則是ㄧ般最多也最容易操作的方法，其追踨法則是在每間隔一段時間，加入一個擾動信號，來檢測目前是否在操作在最大功率點[18-20]。圖 7 所示為以 PWM 為擾動源之最大功率追蹤系統架構。此方法除可改善傳統擾動觀察法的缺點外，亦具有快速追踨最大功率的優點[20]。圖 7 中， Vpv 是太陽能電池端電壓，Ipv 是太陽能電池輸出電流，亦為電感器上的電流，p(t)是計算 Vpv 乘以 Ipv 而得的瞬時功率，Psample 是每一次 PWM 開始時取樣值，誤差量 e(t)為瞬時功率 p(t)減去 Psample，e(t)經由積分器輸出u 與三角波做比較，得到的訊號來調整 PWM 的責任週期，方程式可表示成

∫

−

=

K

p

t

(6)

圖7 最大功率追蹤器架構圖 2.6 換流器 換流器是可將輸入的直流電轉換為可與網路連接的交流電的裝置。利用閘極信號切換電力電子元件的開閉，經過濾波後可得到所需的交流電壓頻率。在本研究中使用弦波比較法改變 PWM 的脈波寬度作為閘極信號，將直流電壓轉換為60Hz 的交流電壓，便可與負載連接[21]。圖 8 為換流器的架構圖。換流器的 A 相電壓方程式可以表示成 o a

V

m

5 .

0

AN

V

=

₍₇₎

其中VAN 為 A 相的相電壓，V0 為昇壓後的直流電壓，ma 為調變指標。換流器的線電壓可表示成 o a

V

m

AN LL

V

(

)

0 .

612

2

3 ₌

=

)

1 (

m

a

≤

(8)

圖8 換流器電路圖 2.7 三相感應電動機模型 三相感應電動機的定子結構為三相對稱繞組，轉子結構則有繞線式與鼠籠式可供選擇。感應電動機的電壓方程式參考到靜止軸後可表示成

(7)

qs s ss m ds s ss qr m r m r rr r dr r m m r rr r

v

r

pL

0 pL

v

0 r

pL

0 v

pL

L

r

pL

v

L

pL

L

r

ω

+

⎡ ⎤ ⎡

⎢ ⎥ ⎢

₊

⎢ ⎥ ⎢

₌

⎢ ⎥ ⎢

−

+

⎢ ⎥ ⎢

₊

⎢ ⎥ ⎢

⎣ ⎦ ⎣

qs m ds rr qr rr dr

0 i

pL

i

L

i

pL

i

⎤ ⎡ ⎤

⎥ ⎢ ⎥

−

⎥ ⎢ ⎥

⎦ ⎣ ⎦

q1 m d 1 q 2 22 d 2 0 i pL i 0 i pL i ⎤ ⎡ ⎤ ⎥ ⎢ ⎥ ⎥ ⎢ ⎥ ⎥ ⎢ ⎥ ⎥ ⎢ ⎥ ⎥ ⎢ ⎥ ⎦ ⎣ ⎦ qk sk dk

i

r

pL

i

⎤ ⎡ ⎤

⎥ ⎢ ⎥

⎦ ⎣ ⎦

(8)

其中vqs 及 iqs 分別是 q 軸定子電壓及電流，vds 及 ids 分別是 d 軸定子電壓及電流，vqr 及 iqr 分別是 q 軸轉子電壓及電流，vdr 及 idr 分別是 d 軸轉子電壓及電流，Lm 是磁化電感，p 是微分運算子[22]、[23]。 2.8 三相變壓器模型 變壓器的主要目的是在相同頻率下，將能量從一個電壓、電流準位轉換至另一個電壓、電流準位。三相變壓器的電壓方程式可表示成 q1 1 11 m d 1 1 11 q 2 m 2 22 d 2 m 2 v r pL 0 pL v 0 r pL 0 v pL 0 r pL v 0 pL 0 r + ⎡ ⎤ ⎡ ⎢ ⎥ ⎢ ₊ ⎢ ⎥ ⎢₌ ⎢ ⎥ ⎢ + ⎢ ⎥ ⎢ ₊ ⎢ ⎥ ⎢ ⎣ ⎦ ⎣

₍₉₎

其中vd1 及 id1 是 d 軸一次側電壓及電流，vq1 及 iq1 是 q 軸一次側電壓及電流，vd2 及 id2 是 d 軸二次側電壓及電流，vq2 及 iq2 是 q 軸二次側電壓及電流，r1 及 r2 是一次側電阻及二次側電阻，L11 及 L22 是一次側自感及二次側自感，Lm 是磁化電感，p 是微分運算子。 2.9 靜態負載模型 集總的靜態負載可考慮成包含電阻性和電感性負載，其電壓方程式可表示成

0

qk sk sk dk sk

v

r

pL

v

+

⎡

⎤ ⎡

=

⎢

⎥ ⎢

₊

⎣

⎦ ⎣

₍₁₀₎

其中vdk 及 idk 是 d 軸電壓及電流，vqk 及 iqk 是 q 軸電壓及電流，rsk 及 Lsk 是電阻及電感，p 是微分運算子。

3. 動態特性模擬

3.1 SimPowerSystems 模組 圖 9 為獨立型太陽能發電系統的 SimPowerSystems 模組架構圖，主要包括四部份—A 部份是太陽能發電系統，B 部份是 PCS，C 部份是靜態負載，D 部份是電動機負載[24]、[25]。

(8)

圖9 獨立型太陽能系統架構圖 3.2 模擬順序 圖10 為模擬的時序圖。由圖中可看出負載順序加入，總模擬時間 12 秒。圖10 模擬時序圖 3.3 模擬結果 圖 11(a)、11(b)顯示太陽能系統輸出電壓與昇壓轉換器輸出電壓會隨著負載的併入而降低。圖 11(c)、 11(d)、11(e)顯示太陽能系統輸出電流、昇壓轉換器輸出電流、以及 PCS 的 A 相電流，會隨著負載的併入而增加。圖11(f)、11(g)顯示太陽能輸出功率會因為負載變動與日照度變化而不同。圖 11(h)顯示 PCS 的電壓在 10 秒後因為日照度較低，加上負載電流造成的壓降，導致明顯的壓降。圖 11(i)顯示 PCS 供應的實功率隨著負載的併入而增加。圖12(a)~12(f)顯示靜態負載的實功率與虛功率都與匯流排電壓變動呈現平方關係。圖12(g)~12(i)顯示電動機的轉速與轉矩都受到匯流排電壓變動的影響。

(9)

4. 結論

本論文主要探討太陽能發電系統獨立運轉時負載變動情形下的動態特性。模擬結果顯示日照度充足時，隨著負載順序併入，太陽能發電系統輸出功率會上升。當日照度不足時，端電壓會下降，特別是有電動機負載時，電壓會降得更低，併聯的靜態負載也會受到相當程度的影響。整體而言，太陽能發電系統在這種運轉模式下的動態特性是可以接受的並且應可和實際運轉情形一致。此研究最大的價值是可以做為太陽能發電系統規劃、運轉、以及擴充的重要依據。

參考文獻

[1] Website , http://www.un.org。 [2] 資策會產業情報研究所(MIC) , http://mic.iii.org.tw/intelligence。 [3] Website , http://www.digitimes.com.tw/。 [4] 台灣電力公司永續報告書 , http://www.taipower.com.tw, 2009。 [5] 大紀元新聞網 , http://www.epochtimes.com/b5/ncnews.htm。

[6] Sharaf A.M. and Yang L., “A Novel Maximum Power Trecking Controller for a Stand-alone Photovoltaic DC Motor Drive,” Canadian Conference on Electrical Engineering and Computer Engineering, pp.450-453, May 2006。

[7] Aldobhani M.S. and John R., “Maximum Power Point Tracking of PV System Using ANFIS Prediction and Fuzzy Logic Tracking,” International MultiConference of Engineers and Computer Scientists, Vol. II, pp978-988, March 2008。 [8] Wu X., Cheng Z., and Wei X., “Maximum Power Point Tracking of Micro PV Systems under Non-uniform Insolation,”

International Conference on Energy and Environment Technology, Vol.2, pp164-167, Oct. 2009。

[9] Bellini and Bifaretti S., “A Quasi-resonant ZCS Boost DC-DC Converter for Photovoltaic Applications,” IEEE International Symposium on Industrial Electronics, pp.815-820, June 2007。

[10] Lee J.P., Min B.D., Kim T.J., Yoo D.W., and Lee B.K., “A Novel Topology for Photovoltaic Series Connected DC/DC Converter with high Efficiency Under wide load Range,” European Conference on Power Electronics and Applications, pp.1-6, Sept. 2007。

[11] Senjyu T., Datta M., Yona A., and Kim C.H., “A Control Method for Small Utility Connected Large PV System to Reduce Frequency Deviation Using a Minimal-Order Observer,” IEEE Transactions on Energy Conversion, Vol. 24, No. 2, pp.520-528, June 2009。

[12] Hudson R.M., Behnke M.R., West R., Gonzalez S., and Ginn J., “Design Considerations for Three-phase grid Connected Photovoltaic Inverters,” IEEE Photovoltaic Specialists Conference,,pp.1396-1401, May 2002。

[13] Kakimoto N., Satoh H., Takayama S., and Nakamura K., “Ramp-Rate Control of Photovoltaic Generator With Electric Double-Layer Capacitor,” IEEE Transactions on Energy Conversion, Vol. 24, No.2, pp465-473, June 2009。

[14] 邱健琮，整合再生能源獨立發電系統之研究，國立勤益科技大學資訊與電能科技研究所，碩士論文，2007。 [15] Gow J. A. and Manning C. D. , “Development of a Photovoltaic Array Model for Use in Power-electronics Simulation

Studies,” IEE Proc. Electric Power Application, Vol.146, No.2, pp.193-200, March 1999。

[16] J.-S Lai, “A High-Performance V6 Converter for Fuel Cell Power Conditioning System,” IEEE Vehicle Power and Propulsion Conference, pp. 624-630, Sept 2005。

[17] 梁適安，交換式電源供給器之理論與實務設計，全華科技圖書公司，2004

[18] Esram T. and Chapman P. L., “Comparison of Photovoltaic Array Maximum Power Point Tracking Techniques”, IEEE Trans. Energy Convers., vol.22, no.2 , pp.439–449, June 2007。

[19] Yongji H. and Deheng L., “A new Method for Optimal Output of a Solar cell Array,” IEEE Proceedings of the Industrial Electronics International Symposium, vol.1, pp.456-459, May 1992。

[20] 楊嘉亨，新型太陽能最大功率追蹤技術之研究，國立高雄應用科技大學電機工程研究所，碩士論文，2007。 [21] Bimal K. Bose, Modern Power Electronics and AC Drives, Prentice-Hall, 2002

[22] Krause P.C., Analysis of Electric Machinery and Drive System, 2nd Ed, McGRAW-Hill Book Co., USA, 2001。 [23] Ong C.-M., Dynamic Simulation of Electric Machinery Using MATLAB/Simulink, McGRAW-Hill Book Co., USA,

1998。

[24] Using Simulink, The Mathworks Inc., 2005。

(10)

0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 時間(秒) PV 輸出電壓 (a) PV 輸出電壓 0 2 4 6 8 10 12 -1.5 -1 -0.5 0 0.5 1 1.5 2 時間(秒) PC S A 相電流 (e) PCS A 相電流 0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 時間(秒) PC S 實功率 (i) PCS 實功率 0 2 4 6 8 10 12 0 0.5 1 1.5 時間(秒) 昇壓轉換器輸出電壓 (b) 昇壓轉換器輸出電壓 0 2 4 6 8 10 12 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 時間(秒) 日照度 (f) PV 日照度 0 2 4 6 8 10 12 -0.2 0 0.2 0.4 0.6 0.8 1 1.2 時間(秒) PC S 虛功率 (j) PCS 虛功率 0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 時間(秒) PV 輸出電流 (c) PV 輸出電流 0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 時間(秒) PV 輸出功率 (g) PV 輸出功率 0 2 4 6 8 10 12 -1.5 -1 -0.5 0 0.5 1 1.5 時間(秒) SL1 A 相電流 (k) SL1 A 相電流 0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 時間(秒) 昇壓轉換器輸出電流 (d) 昇壓轉換器輸出電流 0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 時間(秒) PC S 平均電壓 (h) PCS 平均電壓 0 2 4 6 8 10 12 -1.5 -1 -0.5 0 0.5 1 1.5 時間(秒) SL2 A 相電流 (l) SL2 A 相電流圖11 PV 與 PCS 的變數變動情形

(11)

0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 時間(秒) S L1 平均電壓 (a) SL1 平均電壓 0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 1.2 1.4 3 時間(秒) S L2 平均電壓 (d) SL2 平均電壓 0 2 4 6 8 10 12 0 0.5 1 1.5 2 2.5 時間(秒) Mot o r 實功率 (g) Motor 實功率 1.2 0 2 4 6 8 10 12 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 時間(秒) SL1 實功率 (b) SL1 實功率 0 2 4 6 8 10 12 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 時間(秒) SL2 實功率 (e)SL2 實功率 0 2 4 6 8 10 12 0 0.2 0.4 0.6 0.8 1 時間(秒) Mot o r 轉速 (h) Motor 轉速 4 0 2 4 6 8 10 12 -0.05 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 時間(秒) SL1 虛功率 (c) SL1 虛功率 0 2 4 6 8 10 12 -0.05 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 時間(秒) SL2 虛功率 (f) SL2 虛功率 0 2 4 6 8 10 12 -3 -2 -1 0 1 2 3 時間(秒) Mot o r 轉矩 (i)Motor 轉矩圖12 SL1、SL2 與 Motor 的變數變動情形

獨立型太陽能發電系統在負載變動下的動態特性分析