付兵彬，賈春蓉，楊昌海，張中丹

（甘肅省電力公司電力經(jīng)濟(jì)技術(shù)研究院，蘭州 730050）

光伏并網(wǎng)發(fā)電系統(tǒng)仿真模型的參數(shù)辨識(shí)

付兵彬，賈春蓉，楊昌海，張中丹

（甘肅省電力公司電力經(jīng)濟(jì)技術(shù)研究院，蘭州 730050）

為保證光伏發(fā)電系統(tǒng)輸出結(jié)果的準(zhǔn)確性，開(kāi)展了光伏發(fā)電參數(shù)的辨識(shí)研究。針對(duì)光伏發(fā)電系統(tǒng)，從光伏陣列及并網(wǎng)逆變器兩個(gè)方面進(jìn)行分析，找出影響光伏模型輸出的關(guān)鍵參數(shù)，應(yīng)用遞推最小二乘法制定了光伏并網(wǎng)發(fā)電系統(tǒng)參數(shù)辨識(shí)方法，并利用光伏實(shí)測(cè)數(shù)據(jù)進(jìn)行了光伏電池及光伏逆變器的參數(shù)辨識(shí)，最終得到一組較為準(zhǔn)確的光伏發(fā)電系統(tǒng)參數(shù)。最后將解析計(jì)算、系統(tǒng)辨識(shí)下的仿真輸出與光伏實(shí)測(cè)輸出進(jìn)行了對(duì)比，驗(yàn)證了光伏發(fā)電系統(tǒng)參數(shù)方法及結(jié)果的可行性，為進(jìn)一步研究配電網(wǎng)分布式光伏接入配電網(wǎng)問(wèn)題奠定基礎(chǔ)。

光伏發(fā)電系統(tǒng)；仿真建模；參數(shù)辨識(shí)；最小二乘法

據(jù)預(yù)測(cè)太陽(yáng)能光伏發(fā)電在21世紀(jì)會(huì)占據(jù)世界能源消費(fèi)的重要席位，并將成為世界能源供應(yīng)的主體。隨著電力政策的放寬，光伏發(fā)電接入配電網(wǎng)成為必然趨勢(shì)。為研究光伏大規(guī)模接入對(duì)配電網(wǎng)的影響，國(guó)內(nèi)外科研機(jī)構(gòu)均將電力系統(tǒng)數(shù)字仿真作為研究的重要手段[1]。

光伏發(fā)電數(shù)字仿真技術(shù)的關(guān)鍵是保證光伏發(fā)電模型的輸出特性符合系統(tǒng)實(shí)測(cè)特性，為此需要較為精確合理的光伏發(fā)電系統(tǒng)參數(shù)。在光伏發(fā)電系統(tǒng)控制參數(shù)設(shè)置方面，常用的方法有理論解析法和系統(tǒng)辨識(shí)法，理論解析法能最大限度再現(xiàn)光伏發(fā)電系統(tǒng)內(nèi)部過(guò)程，因此被研究人員廣泛應(yīng)用。但由于光伏系統(tǒng)的復(fù)雜性，理論解析法在求解控制參數(shù)時(shí)遇到非線性微積分方程，且部分不可測(cè)變量會(huì)影響參數(shù)計(jì)算結(jié)果，致使求得的參數(shù)不能直接用于光伏控制，往往需要多次人工調(diào)整；系統(tǒng)辨識(shí)法是利用模型實(shí)測(cè)的輸入和輸出來(lái)求解未知參數(shù)，因此模型最終輸出與系統(tǒng)實(shí)測(cè)的擬合度較高，其在航空航天、經(jīng)濟(jì)系統(tǒng)和電機(jī)控制等領(lǐng)域已得到成功應(yīng)用。但在光伏并網(wǎng)發(fā)電領(lǐng)域，目前僅針對(duì)光伏陣列本身開(kāi)展了辨識(shí)研究工作，對(duì)光伏并網(wǎng)發(fā)電系統(tǒng)的整體尚無(wú)研究成果[2，3]。

本文對(duì)光伏陣列和并網(wǎng)逆變器進(jìn)行了系統(tǒng)性地參數(shù)辨識(shí)。首先分析光伏發(fā)電系統(tǒng)模型，得到影響光伏并網(wǎng)發(fā)電系統(tǒng)模型輸出的關(guān)鍵參數(shù)；然后利用遞推最小二乘法確定了光伏并網(wǎng)發(fā)電系統(tǒng)參數(shù)辨識(shí)方法，并利用光伏實(shí)測(cè)數(shù)據(jù)與光伏并網(wǎng)發(fā)電系統(tǒng)模型進(jìn)行多次辨識(shí)，得到一組確切的控制參數(shù)；最后將系統(tǒng)辨識(shí)與理論解析參數(shù)下的光伏運(yùn)行特性進(jìn)行對(duì)比，驗(yàn)證了光伏發(fā)電控制參數(shù)辨識(shí)的可行性，增強(qiáng)了光伏發(fā)電系統(tǒng)仿真的可信度。

1 光伏發(fā)電系統(tǒng)結(jié)構(gòu)

光伏發(fā)電系統(tǒng)是由光伏陣列將光能轉(zhuǎn)化為直流電能，并通過(guò)逆變器將直流電能變換成交流電能，實(shí)現(xiàn)并網(wǎng)發(fā)電。典型的光伏發(fā)電結(jié)構(gòu)如圖1所示。

圖1 單級(jí)式光伏并網(wǎng)發(fā)電系統(tǒng)結(jié)構(gòu)Fig.1 System structure of single-stage PV grid

圖1中光伏發(fā)電并網(wǎng)主要包括兩個(gè)部分：光伏電池和光伏逆變器。光伏電池的理想等效電路如圖2所示，由光生電流源和一個(gè)二極管并聯(lián)得到。

圖2 光伏電池理想電路模型Fig.2 Idealcircuitmodelof PV cells

其伏安關(guān)系為

式中：V為光伏電池輸出電壓；I為光伏電池輸出電流；Iph為光生電流源電流；Is為二極管飽和電流；q為電子電量常量，值為1.602×10-19C；k為玻耳茲曼常數(shù)，為1.381×10-23J/K。對(duì)于該等值電路，當(dāng)其開(kāi)路時(shí)輸出V對(duì)應(yīng)的開(kāi)路電壓為Voc，ref；其短路時(shí)輸出電流I對(duì)應(yīng)的短路電流為Isc，ref。

為了最大限度的利用太陽(yáng)能，光伏電池中需要使用最大功率點(diǎn)跟蹤MPPT（maximum power point tracking）控制，每隔一定的時(shí)間增加或減小光伏方陣輸出電壓，并觀測(cè)光伏方陣輸出功率的變化，以此決定下一步的調(diào)整方向，此時(shí)等值電路的輸出電壓V和輸出電流I對(duì)應(yīng)的最大功率點(diǎn)電壓為Vmp，ref和最大功率點(diǎn)電流為Imp，ref。綜上，對(duì)于光伏發(fā)電系統(tǒng)需要對(duì)Voc，ref、Isc，ref和Vmp，ref、Imp，ref這4個(gè)參數(shù)進(jìn)行辨識(shí)和控制[4～6]。

光伏系統(tǒng)的逆變器采用電壓功率外環(huán)控制和電流內(nèi)環(huán)控制。其中電壓功率外環(huán)是將光伏陣列經(jīng)過(guò)MPPT控制后的直流電壓和逆變器輸出的無(wú)功功率分別與參考信號(hào)進(jìn)行比較，并對(duì)誤差進(jìn)行PI控制，從而得到內(nèi)環(huán)控制器的參考信號(hào)Idref與Iqref。電流內(nèi)環(huán)控制器采取“dq旋轉(zhuǎn)坐標(biāo)系下的內(nèi)環(huán)控制”，即利用正交派克變換，將三相電壓電流變換到以頻率ω旋轉(zhuǎn)的dq0坐標(biāo)系下，此時(shí)的三相對(duì)稱平衡分量將成為直流量。派克坐標(biāo)變換公式為

在逆變器控制系統(tǒng)中關(guān)鍵參數(shù)有外環(huán)PI參數(shù)有功比例系數(shù)KP、有功積分系數(shù)kiP、無(wú)功比例系數(shù)KQ、無(wú)功積分系數(shù)kiQ和內(nèi)環(huán)PI參數(shù)直軸比例系數(shù)Kd、直軸時(shí)間常數(shù)Td、交軸比例系數(shù)Kq、交軸時(shí)間常數(shù)Tq[7，8]。

2 光伏并網(wǎng)仿真模型的參數(shù)辨識(shí)

參數(shù)辨識(shí)是根據(jù)系統(tǒng)的輸入輸出信息，在某種準(zhǔn)則下，估計(jì)出模型的未知參數(shù)，基本原理如圖3所示。h（k）和z（k）是系統(tǒng)輸入輸出變量，θ為未知的模型參數(shù)。

圖3 參數(shù)辨識(shí)原理Fig.3 Schematic diagram of parameter identification

常用的參數(shù)辨識(shí)方法是最小二乘法，即利用系統(tǒng)的輸入輸出數(shù)據(jù)以最小方差為目標(biāo)對(duì)參數(shù)估計(jì)值進(jìn)行不斷修正，以取得更為準(zhǔn)確的參數(shù)估計(jì)值。對(duì)于從t1到tn時(shí)刻輸入和輸出的觀測(cè)值xi、y以及辨識(shí)參數(shù)θi，構(gòu)建m個(gè)線性方程組為

針對(duì)該系統(tǒng)的誤差ε=y-Xθ，則對(duì)應(yīng)的方差J為

此時(shí)，存在一組參數(shù)θ?滿足

為降低最小二乘法在矩陣維數(shù)較大情況下計(jì)算復(fù)雜或無(wú)法求解的情況，需要加以改進(jìn)和完善。其基本思想是：新的估計(jì)值=舊的估計(jì)值+修正項(xiàng)，即估計(jì)值是利用本次的觀測(cè)數(shù)據(jù)對(duì)預(yù)測(cè)值進(jìn)行修正得到的，隨著新的觀測(cè)數(shù)據(jù)逐次引入，逐步地進(jìn)行參數(shù)估計(jì)，直到估計(jì)值達(dá)到滿意的精度為止[9～11]。

根據(jù)相關(guān)規(guī)定，配電網(wǎng)光伏單點(diǎn)接入容量不大于6MW，本文研究選取1MW光伏發(fā)電系統(tǒng)，并利用數(shù)字仿真軟件DigSILENT進(jìn)行參數(shù)辨識(shí)前后的對(duì)比。其中光伏逆變器和出口變壓器容量取1MW，逆變器交流側(cè)電壓取0.328 kV，逆變器控制采用UdcQ控制，控制功率因數(shù)為0.995，并通過(guò)升壓變壓器以10 kV接入配電網(wǎng)。

圖4 光伏參數(shù)辨識(shí)流程Fig.4 Flow chartof PV parameter identification

光伏發(fā)電控制參數(shù)辨識(shí)的流程如圖4所示。針對(duì)光伏發(fā)電系統(tǒng)的關(guān)鍵控制參數(shù)Voc，ref、Isc，ref、Vmp，ref、Imp，ref，以及逆變器外環(huán)控制參數(shù)KP、TiP、KQ、TiQ和內(nèi)環(huán)控制參數(shù)Kd、Td、Kq、Tq，其目標(biāo)函數(shù)是實(shí)測(cè)數(shù)據(jù)和仿真結(jié)果的差值平方和，即

式中：V1i和I1i分別各次實(shí)測(cè)和仿真中光伏電池的電壓和電流，θi為需要光伏陣列需要辨識(shí)的參數(shù)（Voc，ref，Isc，ref，Vmp，ref，Imp，ref）T，V2i和I2i分別為各次實(shí)測(cè)和仿真輸出的逆變器輸出電壓和電流，ωi為光伏逆變器需要辨識(shí)的外環(huán)控制參數(shù)和內(nèi)環(huán)控制參數(shù)（KP，TiP，KQ，TiQ，Kd，Td，Kq，Tq）T。

3 光伏系統(tǒng)辨識(shí)結(jié)果及驗(yàn)證

首先通過(guò)解析方法對(duì)光伏系統(tǒng)全部控制參數(shù)進(jìn)行計(jì)算，然后設(shè)定光伏發(fā)電系統(tǒng)中逆變器控制參數(shù)保持不變，從而對(duì)光伏電池控制參數(shù)進(jìn)行辨識(shí)，其解析參數(shù)與辨識(shí)后系統(tǒng)參數(shù)如表1所示。

表1 光伏陣列參數(shù)辨識(shí)結(jié)果Tab.1 Parameter identification of PV array

針對(duì)已有光伏電池模型，分別設(shè)定其控制參數(shù)為解析法所得參數(shù)和辨識(shí)法所得參數(shù)，并設(shè)定光照變化在1 000～600W/m2范圍內(nèi)突變，記錄此時(shí)光伏輸出電壓和輸出電流波形，同時(shí)與光伏實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，仿真結(jié)果如圖5和圖6所示。可以看出，當(dāng)光伏系統(tǒng)發(fā)生光照擾動(dòng)時(shí)，光伏并網(wǎng)點(diǎn)電壓、電流均會(huì)有所波動(dòng)，然后恢復(fù)穩(wěn)定狀態(tài)，其中應(yīng)用系統(tǒng)辨識(shí)后參數(shù)的光伏發(fā)電模型輸出與實(shí)測(cè)情況在穩(wěn)態(tài)運(yùn)行和波動(dòng)趨勢(shì)上的擬合度均高于解析法仿真結(jié)果。

圖5 不同參數(shù)下光伏陣列電壓輸出結(jié)果Fig.5 Resultsof PV voltagew ith differentparameters

圖6 不同參數(shù)下光伏陣列電流輸出結(jié)果Fig.6 Resultsof PV currentw ith differentparameters

應(yīng)用辨識(shí)后的光伏陣列控制參數(shù)，對(duì)并網(wǎng)逆變器兩環(huán)控制參數(shù)進(jìn)行參數(shù)辨識(shí)，其解析參數(shù)與辨識(shí)后系統(tǒng)參數(shù)如表2所示。

表2 光伏逆變器參數(shù)辨識(shí)結(jié)果Tab.2 Resultof parameters identification of PV inverter

分別設(shè)定控制參數(shù)為解析法所得參數(shù)和辨識(shí)法所得參數(shù)，并設(shè)定逆變器直流側(cè)電壓在0.5 s時(shí)發(fā)生擾動(dòng)參考值變?yōu)?.95，記錄光伏發(fā)電系統(tǒng)的功率輸出波形，同時(shí)與光伏實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，仿真結(jié)果如圖7所示?？梢钥闯觯捎脜?shù)辨識(shí)后光伏逆變器的功率特性較采用解析法的功率特性能夠更快地恢復(fù)穩(wěn)定狀態(tài)，與實(shí)測(cè)系統(tǒng)變化趨勢(shì)更為接近，能更好地反映光伏發(fā)電系統(tǒng)的實(shí)際情況。

圖7 不同參數(shù)下光伏逆變器的輸出功率波形Fig.7 Outputpower waveform sof PV inverter w ith differentparameters

4 結(jié)論

（1）分析光伏發(fā)電系統(tǒng)的結(jié)構(gòu)及原理得到了影響光伏系統(tǒng)模型輸出的陣列與逆變器的關(guān)鍵參數(shù)。

（2）應(yīng)用遞推最小二乘法研究分布式光伏的參數(shù)辨識(shí)方法，并得到了光伏發(fā)電模型控制參數(shù)。

（3）光伏系統(tǒng)在解析參數(shù)、辨識(shí)參數(shù)運(yùn)行仿真數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的對(duì)比，驗(yàn)證了本文辨識(shí)方法及仿真結(jié)果的準(zhǔn)確性。

值得指出的是，參數(shù)辨識(shí)前的初值設(shè)定與辨識(shí)結(jié)果成功與否有一定制約因素，有待后續(xù)進(jìn)一步完善。但本次研究成果仍有助于掌握光伏發(fā)電系統(tǒng)的運(yùn)行特性，為后續(xù)光伏接入配電網(wǎng)的影響研究提供必要的技術(shù)支撐。

[1]王飛，余世杰，蘇建徽，等（Wang Fei，Yu Shijie，Su Jianhui，etal）．太陽(yáng)能光伏并網(wǎng)發(fā)電系統(tǒng)的研究（Research on photovoltaic grid-connected power system）[J]．電工技術(shù)學(xué)報(bào)（Transactions of China Electrotechical Society），2005，20（5）：72-74，91．

[2]朱豫才．過(guò)程控制的多變量系統(tǒng)辨識(shí)[M]．長(zhǎng)沙：國(guó)防科技大學(xué)出版社，2005．

[3]周建良，王冰，張一鳴，等（Zhou Jianliang，Wang Bing，Zhang Yiming，etal）．基于實(shí)測(cè)數(shù)據(jù)的光伏陣列參數(shù)辨識(shí)與輸出功率預(yù)測(cè)（Parameter identification and output power prediction of photovoltaic array based on themeasured data）[J]．可再生能源（Renewable Energy Resources），2012，30（7）：1-4．

[4]劉邦銀，段善旭，劉飛，等（Liu Bangyin，Duan Shanxu，Liu Fei，etal）．基于改進(jìn)擾動(dòng)觀察法的光伏陣列最大功率點(diǎn)跟蹤（Photovoltaic array maximum power point tracking based on improved perturbation and observation method）[J]．電工技術(shù)學(xué)報(bào)（Transactionsof China Electrotechical Society），2009，24（6）：91-94．

[5]崔巖，蔡炳煌，李大勇，等（Cui Yan，Cai Binghuang，Li Dayong，etal）．太陽(yáng)能光伏系統(tǒng)MPPT控制控制算法的對(duì)比研究（Comparative studieson theMPPT controlalgorithms of solar energy photovoltaic system）[J]．太陽(yáng)能學(xué)報(bào)（Acta Energiae SolarisSinica），2006，27（6）：535-539．

[6]王凌，李培強(qiáng)，李欣然，等（Wang Ling，Li Peiqiang，Li Xinran，etal）．微電源建模及其在微電網(wǎng)仿真中的應(yīng)用（Modeling ofmicro-powers and its application inmicrogrid simulation）[J]．電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)（Proceedingsof the CSU-EPSA），2010，22（3）：32-38．

[7]劉邦銀，段善旭，康勇（Liu Bangyin，Duan Shanxu，Kang Yong）．高能效直流模塊式光伏發(fā)電系統(tǒng)性能評(píng)估（Performance evaluation of high efficiency DC-modulebased photovoltaic power systems）[J]．太陽(yáng)能學(xué)報(bào)（Acta Energiae Solaris Sinica），2008，29（9）：1107-1111．

[8]李廣凱，李庚銀，梁海峰，等（LiGuangkai，LiGengyin，Liang Haifeng，etal）.基于電壓源換流器的光伏并網(wǎng)系統(tǒng)暫態(tài)特性研究（Research on dynamic characteristicsof VSC based photovotaic grid connected system）[J]．太陽(yáng)能學(xué)報(bào)（Acta Energiae Solaris Sinica），2007，28（7）：715-720．

[9]汪海寧，蘇建徽，丁明，等（Wang Haining，Su Jianhui，Ding Ming，etal）．光伏并網(wǎng)功率調(diào)節(jié)系統(tǒng)（Photovoltaic grid connected power conditioner system）[J]．中國(guó)電機(jī)工程學(xué)報(bào)（Proceedingsof the CSEE），2007，27（2）：75-79．

[10]張仰飛，袁越，郝思鵬，等（Zhang Yangfei，Yuan Yue，Hao Sipeng，etal）.數(shù)字PI控制器的參數(shù)辨識(shí)及實(shí)驗(yàn)驗(yàn)證（Parameter identification ofdigitalPIcontrollerand experimentvalidation）[J].電力自動(dòng)化設(shè)備（Electric Power Automation Equipment），2010，30（11）：40-43．

[11]薛凌云（Xue Lingyun）．一種基于過(guò)程辨識(shí)的PID參數(shù)自整定方法（A self-tuningmethod of PID parameter based on process identification）[J]．機(jī)電工程（Mechanical &Electrical EngineeringMagazine），2003，20（3）：34-37．

Parameter Identification of Distributed Photovoltaic Power Generation System

FU Bing-bin，JIA Chun-rong，YANGChang-hai，ZHANG Zhong-dan
（Gansu Electric PowerResearch Instituteof Economicsand Technology，Lanzhou 730050，China）

In order toguarantee the outputaccuracy of the photovoltaic powergeneration system，this paper conductsa research on photovoltaic parameters identification.In photovoltaic powergeneration system，two aspectsofphotovoltaic（PV）array and the grid inverterareanalyzed.The key parametersaffecting PVmodeloutputare found.The parameter identificationmethod for PV grid power system ismade by using recursive leastsquaresmethod.Then，the parameters of photovoltaic PV cellsand PV inverter are identified by using themeasurement data.Finally，a setofmore accurate parameters of photovoltaic power generation system can be obtained.By comparing the simulation outputof the parse calculation and system identification with themeasured outputof PV，the feasibility of parametermethods and results for photovoltaic power generation system is verified，which can lay the foundation for further research on the access problemsofdistribution network distributed photovoltaic.

photovoltaic power generation system；modeling and simulation；parameter identification；least square method

TM744

A

1003-8930（2013）05-0116-05

付兵彬（1972—），男，碩士，高級(jí)工程師，從事電網(wǎng)建設(shè)管理工作。Email：fubb@g s.s g cc.co m.cn

2013-03-16；

2013-06-07

賈春蓉（1973—），女，碩士，高級(jí)工程師，從事電網(wǎng)規(guī)劃工作。Email：jcr730200@163.co m

楊昌海（1982—），男，本科，助理工程師，從事電網(wǎng)規(guī)劃工作。Email：ych0608@126.co m