楊立濱　張海寧　李春來(lái)　楊軍　王平

摘 要： 光伏逆變器是光伏并網(wǎng)系統(tǒng)的核心部件，將基于Hammerstein模型的非線性系統(tǒng)辨識(shí)方法引入到光伏并網(wǎng)逆變器的建模中，把單相光伏并網(wǎng)逆變器視為雙輸入單輸出的非線性黑箱系統(tǒng)。在Hammerstein模型的靜態(tài)非線性環(huán)節(jié)采用B樣條神經(jīng)網(wǎng)絡(luò)，動(dòng)態(tài)線性環(huán)節(jié)采用ARX模型，同時(shí)采用基于誤差學(xué)習(xí)準(zhǔn)則和最小二乘遞歸準(zhǔn)則的自適應(yīng)學(xué)習(xí)方法。實(shí)驗(yàn)測(cè)試結(jié)果表明，提出的BSNN?ARX光伏逆變器模型辨識(shí)方法可以對(duì)不同天氣條件下的逆變器輸出功率進(jìn)行高精度的辨識(shí)，從而為并網(wǎng)逆變器的建模提供一種有效途徑。

關(guān)鍵詞： 光伏逆變器； B樣條神經(jīng)網(wǎng)絡(luò)； ARX模型； 系統(tǒng)辨識(shí)

中圖分類號(hào)： TN711?34； TM615 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1004?373X（2017）07?0167?04

Photovoltaic inverter model identification based on BSNN?ARX

YANG Libin1， ZHANG Haining1， LI Chunlai1， YANG Jun1， WANG Ping2

（1. Electric Power Research Institute， State Grid Qinghai Electric Power Company， Xining 810008， China；

2. State Key Laboratory of Power Transmission Equipment & System Security and New Technology， Chongqing University， Chongqing 400044， China）

Abstract： The photovoltaic inverter is the core component of the photovoltaic grid?connected system. A nonlinear system identification method based on Hammerstein model is introduced into the modeling of the photovoltaic grid?connected inverter， in which the single?phase photovoltaic grid?connected inverter is considered as a nonlinear black?box system with dual inputs and single output. The B?spline neural network （BSNN） is adopted in the static nonlinear link of the Hammerstein model， and the ARX model is adopted in the dynamic linear link. The adaptive learning method based on error learning criterion and least square recursion criterion is employed. The experimental measuring results show that the photovoltaic inverter model identification method based on BSNN?ARX can identify the inverter′s output power with high accuracy under different weather conditions， and provide an effective way to model the grid?connected inverter.

Keywords： photovoltaic inverter； B?spline neural network； ARX model； system identification

0 引 言

光伏發(fā)電作為一種清潔的可再生能源，正處在快速發(fā)展和大規(guī)模應(yīng)用的階段[1?4]。模型參數(shù)的準(zhǔn)確性對(duì)電力系統(tǒng)的分析尤為重要。文獻(xiàn)[5]利用非線性自回歸模型對(duì)光伏發(fā)電系統(tǒng)進(jìn)行建模，所得模型包含了光伏電池和逆變器，可以反應(yīng)不同天氣條件下光伏發(fā)電系統(tǒng)的動(dòng)態(tài)行為，但天氣波動(dòng)的劇烈程度對(duì)辨識(shí)效果影響較大。文獻(xiàn)[6]提出了一種基于果蠅優(yōu)化FOV?Elman神經(jīng)網(wǎng)絡(luò)的光伏電站出力短期預(yù)測(cè)模型，取得了較滿意的預(yù)測(cè)效果。文獻(xiàn)[7]采用改進(jìn)小波神經(jīng)網(wǎng)絡(luò)對(duì)光伏發(fā)電系統(tǒng)進(jìn)行非線性模型辨識(shí)，但其辨識(shí)結(jié)果均受不同天氣狀況的影響較大。

本文針對(duì)現(xiàn)有光伏系統(tǒng)辨識(shí)方法對(duì)不同天氣狀況辨識(shí)能力魯棒性不強(qiáng)的缺點(diǎn)提出一種基于Hammerstein模型[8?9]（簡(jiǎn)稱H模型）的B樣條神經(jīng)網(wǎng)絡(luò)光伏逆變器系統(tǒng)模型辨識(shí)方法。該方法在H模型靜態(tài)非線性環(huán)節(jié)采用B樣條神經(jīng)網(wǎng)絡(luò)（B?Spline Neural Network，BSNN）[10?12]。動(dòng)態(tài)線性環(huán)節(jié)采用自回歸模型（Auto Regressive eXogenous System，ARX），并采用最小二乘法進(jìn)行參數(shù)識(shí)別。為了驗(yàn)證BSNN?ARX模型的光伏逆變器模型辨識(shí)方法的有效性，采用實(shí)測(cè)數(shù)據(jù)對(duì)不同天氣條件下的光伏逆變器進(jìn)行模型辨識(shí)與數(shù)據(jù)結(jié)果分析。

1 B樣條神經(jīng)網(wǎng)絡(luò)ARX模型

1.1 靜態(tài)非線性環(huán)節(jié)

BSNN屬于單層前饋神經(jīng)網(wǎng)絡(luò)，通過(guò)B樣條激勵(lì)函數(shù)實(shí)現(xiàn)輸入到輸出的映射，其基函數(shù)由一些局部多項(xiàng)式組成，本文采用的BSNN結(jié)構(gòu)如圖1所示。

設(shè)輸入信號(hào)[uj]在區(qū)間[[uminj umaxj]]內(nèi)，節(jié)點(diǎn)設(shè)置為[Kj={kj，1，kj，2，…，kj，p}，]區(qū)間劃分為[uminj=kj，1

[Ni，j（uj）=（uj-kj，i）22（Δkj）2， ki，j≤uj

式中：[Δkj=（umaxj-uminj）p，][j=1，2，…，q；i=1，2，…，n；][n]是隱含層基函數(shù)的個(gè)數(shù)，[p]是節(jié)點(diǎn)數(shù)，且[n=p-3；][q]是輸入信號(hào)的維數(shù)。

多輸入單輸出BSNN的輸出信號(hào)[x（u）]可表示為：

[x（u）=j=1qw2j?i=1p-3w1i，j?Ni，j（uj）] （2）

多輸入單輸出BSNN網(wǎng)絡(luò)權(quán)重從輸入到輸出分為兩個(gè)部分：第一部分是每個(gè)BSNN網(wǎng)絡(luò)中神經(jīng)元到輸出的權(quán)重[w1i，j（t）]；第二部分是[q]個(gè)BSNN網(wǎng)絡(luò)輸出求和過(guò)程中的權(quán)重[w2j（t）]。式（2）中第二層的權(quán)值更新公式為：

[Δw2j=η2?i=1p-3w1i，j?Ni，j（uj）?（x-x）?1-w2j（t）， j=1，2，…，q] （3）

第一層的權(quán)值更新變?yōu)椋?/p>

[Δw1i，j=η1，j?Ni，j（uj）?w2j?（x-x）?1-w1i，j（t）， i=1，2，…，n， j=1，2，…，q] （4）

式中：[η]是學(xué)習(xí)因子，為了使BSNN?ARX模型收斂，取[0<η<2m，][m]是模型中插值樣條的次數(shù)；[η2]表示從BSNN輸出到總體輸出的學(xué)習(xí)因子；[η1， j]表示第[j]個(gè)信號(hào)從BSNN神經(jīng)元到輸出的學(xué)習(xí)因子；[x]和[x]分別是BSNN網(wǎng)絡(luò)的目標(biāo)值和實(shí)際輸出值。

1.2 動(dòng)態(tài)線性部分

在H模型的動(dòng)態(tài)線性環(huán)節(jié)采用ARX模型可用式（5）表示：

[y（t）=B（z-1）A（z-1）x（t）] （5）

信號(hào)經(jīng)過(guò)BSNN網(wǎng)絡(luò)后由多輸入變成了單輸出，對(duì)單輸入單輸出線性系統(tǒng)式（5）可以寫成如下表達(dá)式：

[y（t）=-a1y（t-1）-…-any（t-n）+b1x（t-1）+…-bmx（t-m）] （6）

式中：[y（t）]是輸出；[x（t）]是輸入；[n]是[A]的階數(shù)；[m]是[B]的階數(shù)。模型中待估計(jì)的參數(shù)[θT（t-1）]可以寫成：

[θT（t-1）=a1（t-1），…，an（t-1），b1（t-1），…，bm（t-1）] （7）

可以采用遞推算法對(duì)ARX模型進(jìn)行建模，則遞歸向量表示為：

[?T（t-1）=-y（t-1），…，y（t-n），x（t-1），…，x（t-m）]

則式（6）可以進(jìn)一步寫成：

[y0（t+1）=θT（t）?（t）] （8）

式中：[θT（t）]表示[θ（t）]的轉(zhuǎn)置運(yùn)算，ARX部分的遞推最小二乘法參數(shù)估計(jì)如下：

[θ（t+1）=θ（t）+K（t+1）y（t+1）-?T（t+1）θ（t）K（t+1）=P（t）?T（t+1）1+?（t+1）P（t）?T（t+1）P（t+1）=1-K（t+1）?（t+1）P（t）] （9）

式中：[P（t）]是自適應(yīng)增益矩陣，為了啟動(dòng)遞推公式，由式（10）確定初始值[P（0）]和[θ（0）]：

[P（0）=α， α∈（104～1010）θ（0）=ε， ε?1] （10）

2 BSNN?ARX算法實(shí)現(xiàn)

BSNN?ARX模型由單層B樣條神經(jīng)網(wǎng)絡(luò)級(jí)聯(lián)ARX模型組成。每個(gè)神經(jīng)元的激勵(lì)函數(shù)都由式（1）所示的函數(shù)表達(dá)式構(gòu)成。為了對(duì)BSNN及ARX的模型參數(shù)進(jìn)行識(shí)別，將1.1節(jié)和1.2節(jié)所述算法進(jìn)行結(jié)合，自適應(yīng)調(diào)節(jié)部分的具體算法實(shí)現(xiàn)如圖2所示。

自適應(yīng)調(diào)節(jié)部分程序運(yùn)行結(jié)束后，只執(zhí)行前向的計(jì)算，調(diào)整后計(jì)算BSNN的輸出YBSNN作為ARX模型的輸入，再根據(jù)[θ（tL）]計(jì)算ARX模型的輸出，即可得到光伏電站模型的輸出功率。

本文以光伏逆變器直流電壓和電流作為輸入[uii=1，2，]有功功率作為輸出[ytt=1]?；贐SNN?ARX模型的光伏電站模型辨識(shí)流程圖如圖3所示。

3 仿真驗(yàn)證和結(jié)果分析

3.1 實(shí)驗(yàn)數(shù)據(jù)采集與處理

以青海某光伏電站系統(tǒng)為研究對(duì)象，分別采集晴天、多云、雨天三種不同天氣條件下的直流電壓、直流電流及輸出功率數(shù)據(jù)。采樣時(shí)間間隔為15 min，每種天氣類型各采集20天數(shù)據(jù)，將20天的數(shù)據(jù)分成兩部分，前15天作為訓(xùn)練組，實(shí)現(xiàn)模型相關(guān)參數(shù)的估計(jì)，后5天作為模型的驗(yàn)證數(shù)據(jù)。另外，還需要對(duì)數(shù)據(jù)進(jìn)行預(yù)處理，去掉日出前和日落后的無(wú)效數(shù)據(jù)。經(jīng)過(guò)處理的晴天、多云和雨天三種典型天氣條件下的光伏逆變器輸入電壓電流及輸出功率信號(hào)如圖4所示。

3.2 性能評(píng)價(jià)指標(biāo)

單相光伏并網(wǎng)逆變器的辨識(shí)模型與參數(shù)確定后，需要對(duì)其進(jìn)行比較驗(yàn)證，以確保模型的有效性，驗(yàn)證采用如下指標(biāo)：

（1） 平均絕對(duì)誤差：用來(lái)反映預(yù)測(cè)值偏離真實(shí)值誤差的大小。

[eMAE=1ni=1nyi-y*i] （11）

（2） 計(jì)算模型輸出與實(shí)測(cè)數(shù)據(jù)的平均相對(duì)誤差，誤差越小，模型越可靠。

[eMRE=1ni=1nyi-y*iyi] （12）

（3） 計(jì)算模型輸出與實(shí)測(cè)輸出的最佳擬合度，所得數(shù)值最大者的精度最高，即：

[fBestFit=1-y-y*y-y×100%] （13）

式中：[y]為實(shí)測(cè)輸出功率；[y*]為BSNN?ARX模型輸出功率；[y]是實(shí)測(cè)輸出功率[y]的平均值。

3.3 實(shí)驗(yàn)測(cè)試

為了驗(yàn)證提出的BSNN?ARX模型對(duì)光伏并網(wǎng)逆變器建模的有效性，將BSNN?ARX模型與BSNN模型，ARX模型進(jìn)行比較。BSNN神經(jīng)元個(gè)數(shù)取5，激勵(lì)函數(shù)為二次B樣條函數(shù)，表達(dá)式為式（1），學(xué)習(xí)常數(shù)取為0.005。對(duì)晴天、多云和雨天的數(shù)據(jù)進(jìn)行辨識(shí)，輸出功率的辨識(shí)結(jié)果如圖5所示。

從圖5可見(jiàn)，BSNN?ARX模型在不同天氣情況下都可以對(duì)輸出功率進(jìn)行很好的預(yù)測(cè)。表1給出了不同天氣條件下，訓(xùn)練樣本和測(cè)試樣本的預(yù)測(cè)誤差。預(yù)測(cè)誤差采用平均絕對(duì)誤差[eMAE]和平均相對(duì)誤差[eMRE。]平均絕對(duì)誤差反映預(yù)測(cè)值偏離真實(shí)值的誤差大小，平均相對(duì)誤差反映該預(yù)測(cè)誤差的可靠性。

從表1可見(jiàn)，對(duì)于訓(xùn)練樣本，BSNN?ARX模型與BSNN模型相比，平均絕對(duì)誤差在晴天、多云和雨天三種不同天氣情況下分別降低了33.02%，53.83%和56.8%，而對(duì)訓(xùn)練數(shù)據(jù)則分別降低了43.98%，56.79%和66.03%。BSNN?ARX模型與ARX模型相比，對(duì)于訓(xùn)練樣本和測(cè)試樣本，平均絕對(duì)誤差分別下降了62.23%和57.75%。BSNN?ARX模型的平均相對(duì)誤差在三種模型中是最優(yōu)的。

從表2可見(jiàn)，本文提出BSNN?ARX模型的輸出功率在不同天氣條件下均具有較高的最佳擬合度，平均最佳擬合度為94.74%，并且不同天氣條件下最佳擬合度差異很小，對(duì)于天氣條件具有更好的適應(yīng)性。

4 結(jié) 論

本文提出一種基于BSNN?ARX的光伏逆變器模型辨識(shí)方法。該方法采用B樣條神經(jīng)網(wǎng)絡(luò)描述H模型的非線性部分，線性自回歸模型描述H模型的線性部分，采用青海某光伏電站的實(shí)則數(shù)據(jù)對(duì)晴天、多云和雨天三種天氣進(jìn)行BSNN?ARX模型辨識(shí)與實(shí)驗(yàn)測(cè)試。與其他常用的參數(shù)辨識(shí)估計(jì)算法相比，基于BSNN?ARX方法的辨識(shí)模型能夠自適應(yīng)不同天氣條件下的動(dòng)態(tài)變化，并且具有較高的辨識(shí)精度，在該方法下構(gòu)建的辨識(shí)模型對(duì)于天氣條件具有更廣泛的適應(yīng)性。

參考文獻(xiàn)

[1] 周德佳，趙爭(zhēng)鳴，袁立強(qiáng).具有改進(jìn)最大功率跟蹤算法的光伏并網(wǎng)控制系統(tǒng)及其實(shí)現(xiàn)[J].中國(guó)電機(jī)工程學(xué)報(bào)，2008，28（31）：94?100.

[2] HU Bo， NONAKA Yuto， YOKOYAMA Ryuichi. Influence of large scale grid?connected photovoltaic system on distribution networks [J]. Automation of electric power system， 2012， 36（3）： 34?38.

[3] 巢睿祺，王杰.太陽(yáng)能光伏并網(wǎng)系統(tǒng)的控制與仿真[J].電網(wǎng)與清潔能源，2014，30（2）：81?87.

[4] 劉東冉，陳樹勇，馬敏，等.光伏發(fā)電系統(tǒng)模型綜述[J].電網(wǎng)技術(shù)，2011，35（8）：47?52.

[5] XIONG Xiaofu， ZHENG Wei， CHEN Kang， et al. Identification approach in modeling a photovoltaic grid?tie inverter [J]. International review of electrical engineering， 2012， 7（3）： 4575?4584.

[6] 韓偉，王宏華，杜煒.基于FOA?Elman神經(jīng)網(wǎng)絡(luò)的光伏電站短期出力預(yù)測(cè)模型[J].電測(cè)與儀表，2014，51（12）：120?124.

[7] 鄭凌蔚，劉士榮，謝小高.基于改進(jìn)小波神經(jīng)網(wǎng)絡(luò)的光伏發(fā)電系統(tǒng)非線性模型辨識(shí)[J].電網(wǎng)技術(shù)，2011，35（10）：159?164.

[8] 王東風(fēng)，任燕燕，劉長(zhǎng)良，等.基于量子遺傳算法的循環(huán)流化床鍋爐Hammerstein模型辨識(shí)[J].動(dòng)力工程學(xué)報(bào)，2013，33（10）：782?788.

[9] HONG X， MITCHELL R， CHEN S. Modelling and control of Hammerstein system using B?spline approximation and the inverse of De Boor algorithm [J]. International journal of system science， 2012， 43（10）： 1976?1984.

[10] 吳忠強(qiáng)，賈文靜，吳昌韓，等.基于PSO?BSNN的短期風(fēng)速預(yù)測(cè)[J].電力系統(tǒng)保護(hù)與控制，2015，43（15）：36?41.

[11] BEYHAN S， ALCI M. Fuzzy functions based ARX model and new fuzzy basis function models for nonlinear system identification [J]. Applied soft computing， 2010， 10（2）： 439?444.

[12] HONG X， CHEN S. The system identification and control of Hammerstein system using non?uniform rational b?spline neural network and particle swarm optimization [J]. Neurocomputing， 2012， 82（4）： 216?223.