基于在線支持向量機(jī)的風(fēng)速預(yù)測(cè)方法研究

云南能投新能源開(kāi)發(fā)有限公司 ■ 王飛 段家華

0 引言

在全球變暖和能源日益枯竭的今天，對(duì)可再生能源的研究和利用具有十分重要的意義，風(fēng)力發(fā)電已成為近幾年發(fā)展最快的新能源利用技術(shù)之一[1]，而風(fēng)速模型和預(yù)測(cè)是風(fēng)力發(fā)電技術(shù)的研究熱點(diǎn)之一。對(duì)風(fēng)速的預(yù)測(cè)，有助于提前了解即將入網(wǎng)的風(fēng)電功率，從而使電網(wǎng)部門能進(jìn)行合理調(diào)度，保證供電質(zhì)量[2]。對(duì)風(fēng)能資源可利用性的分析很大程度上依賴合適的風(fēng)速模型，并且風(fēng)速模型的準(zhǔn)確性直接影響風(fēng)速模擬結(jié)果數(shù)據(jù)的可靠性。目前存在的一些主流風(fēng)速預(yù)測(cè)方法大致可分為：持續(xù)法[3，4]、卡爾漫濾波法[5，6]、隨機(jī)時(shí)間序列法[7]、模糊邏輯算法[8]、人工神經(jīng)網(wǎng)絡(luò)法[9]和空間相關(guān)性法[10，11]。在線支持向量機(jī)與結(jié)構(gòu)復(fù)雜的神經(jīng)網(wǎng)絡(luò)方法相比，識(shí)別過(guò)程簡(jiǎn)潔快速，準(zhǔn)確率高，能有效降低相對(duì)誤差。因此本文通過(guò)建立風(fēng)速模型獲取實(shí)驗(yàn)數(shù)據(jù)樣本，應(yīng)用在線支持向量機(jī)建模并對(duì)風(fēng)速進(jìn)行預(yù)測(cè)，通過(guò)對(duì)風(fēng)速數(shù)據(jù)的預(yù)測(cè)值和真實(shí)值的比較驗(yàn)證該方法的可行性。

1 建模原理和方法

1.1 風(fēng)速數(shù)學(xué)模型

為了更準(zhǔn)確描述風(fēng)速的隨機(jī)性、間歇性和突變性等特點(diǎn)，構(gòu)建由基本風(fēng)速(vwb)、時(shí)陣風(fēng)速(vwg)、噪聲風(fēng)速(vwn)和漸變風(fēng)速(vwr)4個(gè)部分[12，13]組成的風(fēng)速模型。

1.1.1 基本風(fēng)速模型

在風(fēng)力發(fā)電機(jī)組的整個(gè)工作周期中存在著基本風(fēng)速。基本風(fēng)速是作用于葉輪上的一個(gè)平均風(fēng)速，決定了風(fēng)機(jī)組的出力大小。其表達(dá)式為：

由式(1)可知，基本風(fēng)速不隨時(shí)間變化，可用一個(gè)常數(shù)kb代替，即：

1.1.2 噪聲風(fēng)速

噪聲風(fēng)速用于描述風(fēng)電場(chǎng)葉輪高度某處風(fēng)速的隨機(jī)性[14]。表達(dá)式如下：

式(3)中，vnmax為隨機(jī)風(fēng)速的峰值，m/s ；Ram(-1，1)為-1～1之間均勻分布的隨機(jī)數(shù)；wn為風(fēng)速波動(dòng)的平均距離，一般取0～2π rad/s；φn為0～2π間均勻分布的隨機(jī)量。

1.1.3 時(shí)陣風(fēng)速

時(shí)陣風(fēng)速用于表示模擬風(fēng)速的突然變化的特性[14]。表達(dá)式如下：

式中，Gmax為陣風(fēng)峰值，m/s；t1g為陣風(fēng)開(kāi)始時(shí)間，s；Tg為陣風(fēng)周期，s；t為時(shí)間，s。

1.1.4 漸變風(fēng)速

漸變風(fēng)速用于表示模擬風(fēng)速隨時(shí)間逐漸緩慢變化的特征，表示為：

式中，Rmax為漸變風(fēng)速峰值，m/s；t1r為漸變風(fēng)速開(kāi)始時(shí)間，s；t2r為漸變風(fēng)速結(jié)束時(shí)間，s；t3r為風(fēng)速持續(xù)時(shí)間，s。

1.1.5 實(shí)際風(fēng)速模型

根據(jù)以上4種風(fēng)速的子模型可建立并模擬出實(shí)際風(fēng)速模型為：

1.2 在線支持向量機(jī)

根據(jù)統(tǒng)計(jì)學(xué)理論，最優(yōu)滿足Karush Kuhn Tucker(KKT)條件，得到傳統(tǒng)支持向量機(jī)[15]的對(duì)偶最優(yōu)化問(wèn)題。

式中，Qij=φ(xi)Tφ(xj)T=K(xi·xj)；K(xi·xj)為核函數(shù)。

根據(jù)式(7)得到SVM回歸函數(shù)如下：

由W的一階條件引出KKT條件，因此訓(xùn)練集X中的樣本可分為以下幾類：

1) Set E：誤差支持向量集，E(gi＜0，θi=C)和E(gi*＜0，θi= -C)。

2) Set S：邊界支持向量集，S(gi=0，0＜θi＜C)和S(gi*=0， -C＜θi＜0)。

3) Set R：剩余支持向量集，R(gi＞0，gi*＞0，θi=0)。

在增加新的訓(xùn)練樣本(xc， yc)時(shí)，首先要初始化θc=0且保持訓(xùn)練樣本的參數(shù)θi不變；然后逐步修改θc；并保持在其他訓(xùn)練樣本的KKT條件不變的前提下，調(diào)整對(duì)應(yīng)參數(shù)θi，直到原樣本和新樣本都符合KKT條件。具體算法可參見(jiàn)文獻(xiàn)[16]。

2 實(shí)驗(yàn)分析

2.1 風(fēng)速模型

通過(guò)對(duì)前面風(fēng)速模型的簡(jiǎn)介，文獻(xiàn)[17]已驗(yàn)證了該風(fēng)速模型的正確性，以及與實(shí)際風(fēng)速模型高度相似。在實(shí)驗(yàn)過(guò)程中通過(guò)Matlab/Simulink建立該風(fēng)速模型來(lái)獲取風(fēng)速預(yù)測(cè)的實(shí)驗(yàn)樣本。圖1a為陣風(fēng)與基本風(fēng)的模型圖，圖1b為仿真結(jié)果圖，參數(shù)設(shè)置為vwb=6 m/s，Gmax=3 m/s，t1g=1 s，Tg=20 s。

圖1 陣風(fēng)與基本風(fēng)速模型仿真圖

噪聲風(fēng)速仿真圖如圖2所示，參數(shù)設(shè)置為vnmax=1 m/s，wn=2 π rad/s。

圖2 噪聲風(fēng)速仿真圖

圖3為漸變風(fēng)速仿真圖，參數(shù)設(shè)置為Rmax=3 m/s，t1r=0 s，t2r=10 s，t3r=10 s。

圖3 漸變風(fēng)速仿真圖

通過(guò)以上4部分的合成得到合成風(fēng)速，其結(jié)果如圖4所示。

圖4 風(fēng)速合成圖形

2.2 風(fēng)速預(yù)測(cè)

風(fēng)速模型的建立是為了有效獲得在線支持向量機(jī)的實(shí)驗(yàn)樣本；風(fēng)速預(yù)測(cè)中，模型的預(yù)測(cè)精度很大程度上取決于訓(xùn)練樣本的大小，數(shù)據(jù)太少不足以反映風(fēng)速的全部信息，一般樣本數(shù)至少選擇120個(gè)，預(yù)報(bào)時(shí)效至少可達(dá)24 h，樣本越多，預(yù)測(cè)時(shí)效越長(zhǎng)。根據(jù)風(fēng)速模型實(shí)驗(yàn)結(jié)果，選取200個(gè)數(shù)據(jù)作為訓(xùn)練樣本，100個(gè)數(shù)據(jù)作為測(cè)試樣本；基于風(fēng)速模型在線支持向量機(jī)的數(shù)學(xué)原理在Matlab進(jìn)行編程實(shí)現(xiàn)，建立在線支持向量機(jī)模型，預(yù)測(cè)結(jié)果如圖5所示。

圖5 實(shí)際值與預(yù)測(cè)值對(duì)比圖

由圖5可知，預(yù)測(cè)風(fēng)速與樣本風(fēng)速最大相對(duì)誤差為11.45%，平均誤差為3.02%，很大程度降低了實(shí)驗(yàn)誤差，驗(yàn)證了該方法的可行性。實(shí)驗(yàn)中由于風(fēng)速模型中構(gòu)建包含隨機(jī)風(fēng)速，存在不確定性，導(dǎo)致個(gè)別點(diǎn)預(yù)測(cè)結(jié)果偏差較大，但該方法的大多數(shù)數(shù)據(jù)都能取得較好的預(yù)測(cè)結(jié)果，具有較好的實(shí)用價(jià)值。

3 結(jié)論

本文對(duì)基于在線支持向量機(jī)的風(fēng)速預(yù)測(cè)模型進(jìn)行了討論，在線支持向量機(jī)所需樣本少，避免了維數(shù)災(zāi)難，有效減小了每次訓(xùn)練樣本的大小，提高了訓(xùn)練速度并減少了內(nèi)存的占用；在風(fēng)速預(yù)測(cè)中，很多程度上提高了預(yù)測(cè)精度，為風(fēng)電場(chǎng)風(fēng)速的研究提供了重要基礎(chǔ)。

[1] 任艷鋒. 基于最大風(fēng)能追蹤的風(fēng)速預(yù)測(cè)研究[D]. 蘭州： 蘭州理工大學(xué)， 2010.

[2] 王金翠. 基于實(shí)測(cè)數(shù)據(jù)的風(fēng)電場(chǎng)風(fēng)速和風(fēng)功率短期預(yù)測(cè)研究 [D]. 吉林 ： 東北電力大學(xué) ， 2010.

[3] Cameron W．Potter， Michael Negnevitsky．Very short term wind forecasting for tasmanian power generation[J]. IEEE Transactions on Power Systems， 2006， 21(2)： 965 － 972.

[4] Alexiadis M，Dokopoulos P，Sahsamanoglou H， et a1. Short term forecasting of wind speed and related electrical power[J].Solar Energy， 1998， 63(1)： 61 － 68．

[5] 吳國(guó)旸， 肖洋， 翁莎莎. 風(fēng)電場(chǎng)風(fēng)速預(yù)測(cè)探討[J]. 吉林電力，2005， 181(6)： 21 － 24.

[6] 潘迪夫， 劉輝， 李燕飛. 基于時(shí)間序列分析和卡爾曼濾波算法的風(fēng)電場(chǎng)風(fēng)速預(yù)測(cè)優(yōu)化模型[J]. 電網(wǎng)技術(shù)， 2008， 32(7)：82－86.

[7] Guoyang W， Yang X， Shasha W. Discussion about short-term forecast of wind speed on wind farm [J]. Jilin Electric Power，2005， 181(5)： 21 － 24.

[8] Damousis I G， Alexiadis M C， Theocharis J B， et al. A fuzzy model for wind speed prediction and power generation in wind parks using spatial correlation[J]. Energy Conversion， 2004， 19 (2)：352－361.

[9] Wang X， Sideratos G， Hatziargyriou N， et al. Wind speed forecasting for power system operational planning[A]. 2004 International Conference on Probablitistic Methods Applied to Power System[C]， Ames， IA， 2004.

[10] 杜穎， 盧繼平， 李青， 等. 基于最小二乘支持向量機(jī)的風(fēng)電場(chǎng)短期風(fēng)速預(yù)測(cè) [J]. 電網(wǎng)技術(shù) ， 2008， 32(15)： 62 － 66.

[11] 徐蓓蓓， 蔣鐵錚， 易宏. 三種基于支持向量機(jī)風(fēng)速預(yù)測(cè)方法對(duì)比研究 [J]. 電器技術(shù) ， 2012， (5)： 22 － 25.

[12] 吳學(xué)光 ， 張學(xué)成 ， 印永華 ， 等 . 異步風(fēng)力發(fā)電系統(tǒng)動(dòng)態(tài)穩(wěn)定性分析的數(shù)學(xué)模型及其應(yīng)用[J]. 電網(wǎng)技術(shù)， 1998， 22(6)：68-72.

[13] 郭永麗， 吳健， 溫步瀛， 等. 變速風(fēng)力機(jī)的建模與仿真[J].福建電力與電工 ， 2008， 28(3)： 1 － 4.

[14] 洪秀麗. 風(fēng)電場(chǎng)模型研究及其在DDRTS中的應(yīng)用[D]. 北京 ： 北京交通大學(xué) ， 2009.

[15] 孫鳳. 傳送帶給料生產(chǎn)加工站的神經(jīng)元在線優(yōu)化算法[D].安徽 ： 合肥工業(yè)大學(xué) ， 2010.

[16] Vapnik V N. The nature of statistical learning theory [M].New York： Springer-Verlag， 1995.

[17] 魏毅立， 韓素賢， 時(shí)盛志. 風(fēng)力發(fā)電系統(tǒng)中組合風(fēng)速的建模及仿真 [J]. 可再生資源 ， 2010， 28(2)： 18 － 20.