吳堅(jiān)，項(xiàng)頌，閻誠，吳曉丹，馬繼濤，劉福鎖

（1.國網(wǎng)內(nèi)蒙古東部電力有限公司，內(nèi)蒙古呼和浩特 010000；2.東南大學(xué)電氣工程學(xué)院，江蘇南京210096；3.國電南瑞科技股份有限公司，江蘇南京 211106）

0 引言

風(fēng)能作為一種高效的清潔能源，具有強(qiáng)烈的隨機(jī)性、間歇性和逆負(fù)荷性等特點(diǎn)[1]。據(jù)統(tǒng)計(jì)，2019年全球風(fēng)電新增和累積裝機(jī)容量較2018年同比增長(zhǎng)19%和10%[2]。為了降低高滲透率風(fēng)電對(duì)電網(wǎng)帶來的不確定性，提高風(fēng)電消納能力，實(shí)現(xiàn)系統(tǒng)安全穩(wěn)定運(yùn)行，準(zhǔn)確的風(fēng)電超短期功率預(yù)測(cè)具有重要的意義[3]。

風(fēng)電功率時(shí)間序列具有很強(qiáng)的時(shí)序特征。風(fēng)電功率序列預(yù)測(cè)方法以統(tǒng)計(jì)學(xué)方法和人工智能算法為主。統(tǒng)計(jì)學(xué)方法僅依賴時(shí)間序列方法進(jìn)行建模，能夠有效捕捉時(shí)間序列中的非線性關(guān)系，但存在模型定階和參數(shù)估計(jì)問題[4]。人工智能算法能夠有效捕捉風(fēng)電功率時(shí)間序列中的非線性關(guān)系，但對(duì)樣本數(shù)據(jù)質(zhì)量依賴較高[5]，其中，最小二乘支持向量機(jī)（Least-Squares Support Vector Machine，LSSVM）通過映射關(guān)系將原支持向量機(jī)模型參數(shù)訓(xùn)練中的非線性問題替換成線性問題，能有效增強(qiáng)模型計(jì)算效率，具有較強(qiáng)的泛化能力[6]。對(duì)于風(fēng)電功率序列復(fù)雜的波動(dòng)特性，單一預(yù)測(cè)模型因受制于樣本維度和數(shù)量等因素而具有局限性[7]。

基于數(shù)據(jù)預(yù)處理過程的組合預(yù)測(cè)模型是提升單一預(yù)測(cè)方法適應(yīng)能力和模型精度的重要途徑[8]。通過對(duì)原始高頻非線性時(shí)間序列進(jìn)行分解，能夠有效降低序列的復(fù)雜程度，采用分解后的子序列分別建立一系列適應(yīng)其序列特征的預(yù)測(cè)模型，顯著提升了預(yù)測(cè)性能。傳統(tǒng)數(shù)據(jù)預(yù)處理方法包括以經(jīng)驗(yàn)?zāi)B(tài)分解（Empirical Mode Decomposition，EMD）[9]和小波分析（Wavelet Decomposition，WD）[10]為基礎(chǔ)衍生出的一系列序列分解方法，但是分解序列結(jié)果存在顯著的端點(diǎn)效應(yīng)和模態(tài)混疊現(xiàn)象，極大地增強(qiáng)了預(yù)測(cè)過程中的難度。而奇異譜分析（Singular Spectrum Analysis，SSA）作為一種非線性時(shí)間序列分解重構(gòu)的分析方法，能夠有效處理復(fù)雜序列，無須選取先驗(yàn)基函數(shù)和復(fù)雜的運(yùn)算過程，具有較好的客觀性和自適應(yīng)性[11]，使其廣泛地應(yīng)用于各工程學(xué)科領(lǐng)域。文獻(xiàn)[12]采用SSA提取風(fēng)速序列的趨勢(shì)成分、振蕩成分和噪聲成分，分別建立預(yù)測(cè)模型，實(shí)現(xiàn)了準(zhǔn)確可靠的風(fēng)速區(qū)間預(yù)測(cè)。但是SSA存在嵌入維度選取困難的問題，過高的嵌入維度導(dǎo)致關(guān)鍵序列分量被過度分解為噪聲成分，而較小的嵌入維度無法有效提取序列的關(guān)鍵變化特征[13]。

考慮到風(fēng)電功率時(shí)間序列包含大量的時(shí)序性和自相關(guān)性特征，本文提出了一種基于SSA和LSSVM的組合預(yù)測(cè)模型。針對(duì)SSA嵌入維度存在不確定性的問題，使用Cao方法判別嵌入維度。對(duì)于分解得到的子序列分別使用LSSVM進(jìn)行訓(xùn)練，通過粒子群優(yōu)化算法（Particle Swarm Optimization，PSO）確定各子序列的LSSVM模型參數(shù)，實(shí)現(xiàn)對(duì)風(fēng)電功率的多步預(yù)測(cè)。引入樣本熵作為分解模型效果度量指標(biāo)，結(jié)合實(shí)測(cè)數(shù)據(jù)，驗(yàn)證了基于SSA的風(fēng)電功率超短期預(yù)測(cè)模型在多步預(yù)測(cè)中的優(yōu)越性能。

1 整體研究思路

通過對(duì)風(fēng)電功率時(shí)間序列進(jìn)行分解，根據(jù)所得到的不同分量的波動(dòng)特性，分別建立預(yù)測(cè)模型，以降低預(yù)測(cè)誤差。采用組合模型進(jìn)行預(yù)測(cè)的整體研究思路如圖1所示。

圖1 組合預(yù)測(cè)模型Fig.1 Diagram of combined forecasting model

首先，基于復(fù)雜系統(tǒng)混沌理論，采用Cao方法確定嵌入維度的選取，即分解算法所得到的分量數(shù)量，通過對(duì)原始訓(xùn)練樣本進(jìn)行SSA分解，得到一系列不同時(shí)間尺度的分量序列；然后，采用LSSVM分別對(duì)其進(jìn)行訓(xùn)練，采用PSO優(yōu)化預(yù)測(cè)模型的基本參數(shù)。類似的，在功率預(yù)測(cè)階段，將實(shí)時(shí)采集到的功率數(shù)據(jù)與歷史數(shù)據(jù)進(jìn)行拼接，根據(jù)選取的嵌入維度重新進(jìn)行分解，并采用LSSVM對(duì)各分量序列進(jìn)行預(yù)測(cè)；最后，疊加不同分量的風(fēng)電功率序列預(yù)測(cè)輸出。

2 基于SSA的風(fēng)電功率序列分解

2.1 SSA原理

SSA是基于奇異值分解的一種序列成分分析方法，與時(shí)序預(yù)測(cè)模型結(jié)合能有效提升預(yù)測(cè)性能。對(duì)于具有非線性特征的風(fēng)電功率序列樣本，通過構(gòu)建軌跡矩陣、奇異值分解、分組和對(duì)角平均化等流程實(shí)現(xiàn)風(fēng)電功率序列的分解[11]。具體流程如下。

提取樣本數(shù)量為N的風(fēng)電歷史功率序列z=[z1，z2，…，zN]，設(shè)定SSA的模型嵌入維度為L(zhǎng)，將其重構(gòu)為跡矩陣X：

在此基礎(chǔ)上，可根據(jù)奇異值對(duì)各子序列進(jìn)行分組，篩選分解后的各類型分量，剔除諧波分量。在本文中，待預(yù)測(cè)風(fēng)電功率序列已經(jīng)進(jìn)行了預(yù)處理剔除異常數(shù)據(jù)，因此在SSA階段中不考慮分組的情況。

對(duì)角平均化還原時(shí)間序列，當(dāng)X是一個(gè)L×K的矩陣，矩陣中元素為xij，令L*=min（L，K），K*=max（L，K），N=L+K-1，且滿足當(dāng)L＜K時(shí)，=xij；否則=xji，對(duì)子矩陣Xi進(jìn)行重構(gòu)得到的序列[y1，y2，…，yN]可表示為

SSA具有很強(qiáng)的非參數(shù)性，因此其嵌入維度L決定了模型分解性能。對(duì)于復(fù)雜波動(dòng)的風(fēng)電功率序列，L越大會(huì)削弱不同分量之間特征值差異，利于周期規(guī)律變化的風(fēng)電功率的成分提取，但會(huì)產(chǎn)生多余分量；L較小能有效反映風(fēng)電功率序列的波動(dòng)態(tài)勢(shì)，但無法充分挖掘分量中蘊(yùn)含的關(guān)鍵信息。L一般通過已知數(shù)據(jù)進(jìn)行經(jīng)驗(yàn)選取，不能大于給定序列長(zhǎng)度的1/2。而風(fēng)電功率序列兼具短時(shí)間尺度的高頻波動(dòng)特性和長(zhǎng)時(shí)間尺度變化趨勢(shì)，因此合理的L選取具有重要的意義。SSA中構(gòu)建跡矩陣與混沌系統(tǒng)理論中的相空間重構(gòu)一致，可借鑒混沌序列相空間重構(gòu)嵌入維度的選取方法。

2.2 Cao方法

Cao方法具有較強(qiáng)的適應(yīng)性和計(jì)算效率[14]，采用Cao方法計(jì)算風(fēng)電功率序列分解的L步驟如下。

計(jì)算不同L下最近鄰點(diǎn)的距離變化值：

式中：i=1，2，…，N-L；||·||為∞-范數(shù)；Xi（L+1）為L(zhǎng)+1的第i個(gè)重構(gòu)得到的相空間向量。

對(duì)于L+1，Xn（i，M）（L+1）為∞-范數(shù)下與Xi（L+1）距離最近的向量，n（i，M）為滿足1≤n（i，M）≤（N-L）的整數(shù)。

計(jì)算相同L下空間距離波動(dòng)值的均值：

采用E1（L）監(jiān)測(cè)不同L下E（L）的變化情況：

當(dāng)滿足L＞L0時(shí)，對(duì)應(yīng)E1數(shù)值停止變化或波動(dòng)小于一定閾值，得到最小L為L(zhǎng)0+1。為了保證L選取的客觀性和合理性，引入E1波動(dòng)率ei。

設(shè)定閾值e0和波動(dòng)邊界k，滿足1＜k＜N-1，k為當(dāng)?shù)谝淮螡M足ei＜e0時(shí)的下標(biāo)。

在此基礎(chǔ)上，定義e：

當(dāng)i+1滿足ei+2＜ei+1＜ei，且ei+1＜e時(shí)，得到最佳L=i+1。

2.3 樣本熵

為對(duì)比不同分解方法下風(fēng)電功率序列的分解效果，引入樣本熵描述其復(fù)雜程度。樣本熵越大的序列，復(fù)雜程度越大。作為一種復(fù)雜性度量方法，樣本熵對(duì)數(shù)據(jù)長(zhǎng)度不敏感，具有較好的穩(wěn)定性。對(duì)于長(zhǎng)度為N的時(shí)間序列X=[X（1），X（2），…，X（N）]，其樣本熵計(jì)算步驟如下[15]。

基于重構(gòu)維數(shù)m對(duì)原風(fēng)電功率分解子序列進(jìn)行重構(gòu)，得到一組m維的空間向量Xm：

Xm（i）表示為

定義運(yùn)算Xm（i）和Xm（j）之間距離的絕對(duì)值為

計(jì)算Dm{Xm（i），Xm（j）}小于設(shè)定相似度閾值r的數(shù)量，并記為，Bm。

樣本熵計(jì)算式為

基于一致性原理，通常取樣本熵m為1，r為0.15倍風(fēng)電功率序列統(tǒng)計(jì)標(biāo)準(zhǔn)差[16]。

3 SSA-LSSVM風(fēng)電功率組合預(yù)測(cè)模型

3.1 LSSVM原理

本文所研究的風(fēng)電超短期功率預(yù)測(cè)，要求風(fēng)電場(chǎng)內(nèi)輸出功率序列進(jìn)行就地和實(shí)時(shí)的準(zhǔn)確預(yù)測(cè)，因此對(duì)于預(yù)測(cè)模型精度和復(fù)雜程度提出了一定的要求。本文通過LSSVM構(gòu)建風(fēng)電功率分解后子序列的預(yù)測(cè)模型，而LSSVM作為支持向量機(jī)（Support Vector Machine，SVM）的改進(jìn)模型，將SVM模型訓(xùn)練中復(fù)雜的二次規(guī)劃問題轉(zhuǎn)變成線性方程組問題[17]，更有利于在工程實(shí)際中的應(yīng)用?？杀硎緸?/p>

式中：φ（x）為核空間映射函數(shù)；w為權(quán)重；b為偏差量。

根據(jù)結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則，LSSVM的參數(shù)優(yōu)化模型可表示為

式中：ζ為松弛因子；樣本數(shù)量滿足i=1，2，…，N。

構(gòu)造拉格朗日函數(shù)如下：

式中：αi為拉格朗日乘子。

根據(jù)Karush-Kuhn-Tucker最優(yōu)條件，得到?L/?ω=0，?L/?b=0，?L/?ξ=0和?L/?αi=0，模型中權(quán)值系數(shù)w和誤差變量ξ可進(jìn)行化簡(jiǎn)消去，最終得到方程組：

式中：l=[1，1，…，1]T；K（xi，xj）為根據(jù)Mercer條件的核函數(shù)，文中LSSVM的核函數(shù)選取徑向基核函數(shù)（Radial Basis Function，RBF）。

求解得到LSSVM模型為

LSSVM模型的性能很大程度取決于懲罰系數(shù)γ以及核函數(shù)的參數(shù)σ的取值，文中采用的PSO為優(yōu)化算法，以LSSVM輸出和實(shí)測(cè)值的均方根誤差作為PSO算法的適應(yīng)度函數(shù)，對(duì)超參數(shù)γ和σ進(jìn)行尋優(yōu)，具體原理可見文獻(xiàn)[18]。

3.2 SSA-LSSVM組合預(yù)測(cè)模型

風(fēng)電功率序列在短時(shí)間尺度存在高非線性和波動(dòng)性的特征。單一預(yù)測(cè)方法無法在預(yù)測(cè)過程中對(duì)這些特征進(jìn)行準(zhǔn)確捕捉，因此預(yù)測(cè)結(jié)果普遍存在一定的時(shí)延性。本文結(jié)合SSA處理非平穩(wěn)序列的優(yōu)勢(shì)和LSSVM在線高效計(jì)算的特點(diǎn)，提出了一種SSA-LSSVM的風(fēng)電超短期功率組合預(yù)測(cè)模型。首先采用Cao方法確定功率時(shí)間序列的最優(yōu)L，對(duì)原始序列進(jìn)行SSA分解，得到一系列多時(shí)間尺度的固有模態(tài)分量（Intrinsic Mode Function，IMF），實(shí)現(xiàn)序列的關(guān)鍵特征提??；然后分別采用LSSVM對(duì)各分量進(jìn)行訓(xùn)練，通過PSO優(yōu)化模型參數(shù)；最后通過測(cè)試樣本的預(yù)測(cè)值疊加得到實(shí)際風(fēng)電功率的預(yù)測(cè)結(jié)果。具體步驟如下：

①將功率序列進(jìn)行訓(xùn)練集和測(cè)試集劃分，以驗(yàn)證預(yù)測(cè)模型的泛化能力。初始化模型參數(shù)：預(yù)測(cè)步長(zhǎng)T和PSO中的尋優(yōu)參數(shù)，包括粒子種群數(shù)量、粒子的位置和速度等；

②在離線模型訓(xùn)練階段，采取Cao方法確定L，對(duì)序列進(jìn)行SSA分解，采用PSO選取各分量的最佳參數(shù)，訓(xùn)練LSSVM模型；

③在線預(yù)測(cè)階段，根據(jù)T選取測(cè)試樣本，對(duì)其進(jìn)行SSA分解，以此輸入LSSVM中生成預(yù)測(cè)結(jié)果。根據(jù)T調(diào)整LSSVM輸入值，滾動(dòng)迭代完成該步長(zhǎng)下所有分量的預(yù)測(cè)后，將各分量預(yù)測(cè)值累加得到真實(shí)功率預(yù)測(cè)值；

④重復(fù)步驟③，直到完成整個(gè)預(yù)測(cè)周期的功率預(yù)測(cè)。

具體功率預(yù)測(cè)模型流程如圖2所示。

圖2 組合預(yù)測(cè)模型流程圖Fig.2 Flow chart of combined forecasting model

值得注意的是，采用組合模型通過序列分解的方式進(jìn)行預(yù)測(cè)會(huì)存在一個(gè)誤區(qū)，就是對(duì)包含訓(xùn)練樣本和測(cè)試樣本的原始風(fēng)速序列直接進(jìn)行分解，然后將各分量帶入相應(yīng)的模型進(jìn)行預(yù)測(cè)。因?yàn)閷?duì)含有未來數(shù)據(jù)的序列進(jìn)行分解，得到的分量結(jié)果本身包含了未來真實(shí)值波動(dòng)的一個(gè)趨勢(shì)，所以上述方法在實(shí)際工程應(yīng)用中是不可能實(shí)現(xiàn)的。采用組合模型進(jìn)行預(yù)測(cè)時(shí)，必須在每一步預(yù)測(cè)前對(duì)含有當(dāng)前功率真實(shí)值的序列重新進(jìn)行分解。同時(shí)由于分解算法原理上的問題，當(dāng)待分解的時(shí)間序列邊緣的真實(shí)值增加到一定數(shù)量時(shí)，模型輸出的各模態(tài)分量會(huì)與原始測(cè)試樣本產(chǎn)生的分量存在一定程度的畸變。采用EMD方法進(jìn)行預(yù)測(cè)時(shí)，由于其模態(tài)分解數(shù)量不可控性，待分解的序列不斷更新，EMD輸出模態(tài)分量的數(shù)量也會(huì)發(fā)生改變，而SSA可以克服上述缺點(diǎn)。基于上述考慮，由于每一步預(yù)測(cè)均須要進(jìn)行SSA分解以更新各分量的真實(shí)值且計(jì)算量大，因此組合模型不適用于單步預(yù)測(cè)和中長(zhǎng)期的預(yù)測(cè)。

為了克服上述問題，文中采用的多步預(yù)測(cè)模型采用迭代法，即每次只向前預(yù)測(cè)一步，且每步都采用新得到的預(yù)測(cè)值而不是實(shí)測(cè)值來繼續(xù)下一步的預(yù)測(cè)。相較于直接向前預(yù)測(cè)T步的直接法，迭代法多步預(yù)測(cè)具有更高的精度[19]。文中通過虛假臨近點(diǎn)法確定LSSVM輸入維度[20]。

4 算例及結(jié)果分析

算例采用中國某風(fēng)電場(chǎng)1.5 MW風(fēng)機(jī)在2015年6月19-26日SCADA采集的功率序列為實(shí)驗(yàn)樣本，輸入風(fēng)速采樣間隔為10 min，共計(jì)1 152個(gè)功率采樣點(diǎn)，選取前1 008個(gè)采樣點(diǎn)為模型實(shí)驗(yàn)樣本，后144個(gè)點(diǎn)為測(cè)試樣本。測(cè)試前已進(jìn)行數(shù)據(jù)清洗，故在SSA分析過程中不再考慮對(duì)樣本重構(gòu)分量的篩選。實(shí)驗(yàn)在MATLAB 2019a平臺(tái)上進(jìn)行。選取Cao方法維度篩選閾值e0=0.01，Cao方法中所定義E1的波動(dòng)情況和不同維度下最大最小特征值的比值λmax/λmin如圖3所示。

圖3 Cao方法選取嵌入維度Fig.3 Select embedding dimension using Cao method

由圖3可以看出，當(dāng)L達(dá)到13時(shí)，E1停止波動(dòng)，同時(shí)λmax/λmin與L呈線性遞增關(guān)系上漲。L過大會(huì)造成計(jì)算量顯著增加和關(guān)鍵特征序列被淹沒，而Cao方法不須要依賴經(jīng)驗(yàn)，可確定最佳L?；诖?，對(duì)訓(xùn)練樣本進(jìn)行SSA分解，構(gòu)成13個(gè)特征值遞減的模態(tài)分量，如圖4所示。由圖可以看出，SSA分量為1的奇異值顯著高于其他分量，說明該分量對(duì)原風(fēng)電功率軌跡矩陣貢獻(xiàn)值最大。

圖4 SSA分解序列奇異值Fig.4 Comparison of Singular value of SSA decomposition sequence

為了驗(yàn)證所提模型的預(yù)測(cè)性能，本文采用PSO-LSSVM[18]和以EMD方法、完備經(jīng)驗(yàn)?zāi)B(tài)分解（Ensemble EMD，EEMD）以及改進(jìn)總體完備經(jīng)驗(yàn)?zāi)B(tài)分解（Improved Complete Ensemble EMD，ICEEMD）[21]為基礎(chǔ)，與LSSVM結(jié)合的預(yù)測(cè)方法進(jìn)行對(duì)比分析，以下分別稱為EMD-LSSVM，EEMD-LSSVM和ICEEMD-LSSVM方法。EEMD在EMD的基礎(chǔ)上通過增加白噪聲降低EMD的模態(tài)混疊現(xiàn)象。進(jìn)一步，ICEEMD通過利用添加自適應(yīng)白噪聲和選取各模態(tài)分量當(dāng)前的局部均值以重新定義真實(shí)模態(tài)的方式對(duì)EMD方法進(jìn)行改進(jìn)，克服了虛假模態(tài)和分量存在延遲特性等問題，是當(dāng)前EMD一類算法應(yīng)用于預(yù)測(cè)模型中的典型代表，具體算法流程可參見文獻(xiàn)[22]。

樣本熵在度量低維數(shù)據(jù)復(fù)雜程度上具有巨大的優(yōu)勢(shì)，通過計(jì)算可得原始時(shí)間序列樣本熵為1.606 9，各方法樣本熵對(duì)比結(jié)果如表1所示。

表1 不同分解方法樣本熵計(jì)算結(jié)果Table 1 Calculation results of sample entropy for different decomposition methods

表1顯示經(jīng)SSA分解后的各分量信號(hào)仍保持較高的復(fù)雜性特征，而只通過EMD類方法分解得到的分量平均樣本熵遠(yuǎn)低于原始序列樣本熵。為了量化分析不同預(yù)測(cè)方法和步長(zhǎng)下多步預(yù)測(cè)時(shí)的預(yù)測(cè)精度，本文采用歸一化均方根誤差（Normalized Root-Mean-Squared Error，NRMSE）、歸一化絕對(duì)平均誤差（Normalized Mean Absolute Error，NMAE）和歸一化絕對(duì)誤差（Normalized Absolute Error，NAE），定義如下：

式中：y0（i）為真實(shí)值；y（i）為預(yù)測(cè)值；N為預(yù)測(cè)序列總長(zhǎng)度；PN為風(fēng)機(jī)裝機(jī)容量。

預(yù)測(cè)模型的四步預(yù)測(cè)結(jié)果及其對(duì)應(yīng)的不同預(yù)測(cè)方式下的NAE如圖5所示。

圖5 四步預(yù)測(cè)結(jié)果及誤差Fig.5 Curves of four-step wind power prediction results and error

為了驗(yàn)證文中所提方法的有效性，本文將直接通過LSSVM預(yù)測(cè)和以EMD類分解方法為基礎(chǔ)的LSSVM預(yù)測(cè)結(jié)果與所提出的方法進(jìn)行對(duì)比。由圖5可以看出，SSA-LSSVM方法能夠在多步預(yù)測(cè)過程中很好地捕捉原始序列的變化趨勢(shì)，精度較未采用序列分解而直接預(yù)測(cè)的結(jié)果有顯著的提升。由于對(duì)時(shí)間序列的充分分解，原來具有高頻波動(dòng)特征的風(fēng)電功率序列的復(fù)雜程度被充分降低，有利于預(yù)測(cè)模型的擬合。EMD類方法無法主觀設(shè)定分解維數(shù)，因此在預(yù)測(cè)過程中會(huì)出現(xiàn)分解維度突變情況，適應(yīng)性較差。ICEEMD相較于EMD和EEMD方法，克服了數(shù)據(jù)預(yù)處理過程中的端點(diǎn)效應(yīng)，但是在風(fēng)電功率突變點(diǎn)，預(yù)測(cè)結(jié)果與SSALSSVM方法預(yù)測(cè)結(jié)果變化趨勢(shì)有明顯差別，說明單獨(dú)的ICEEMD模態(tài)分解方法無法對(duì)原序列的多時(shí)間尺度分量實(shí)現(xiàn)充分的分解，整體預(yù)測(cè)精度低于SSA-LSSVM方法。而只采用LSSVM的方法受訓(xùn)練樣本容量的限制，超短期預(yù)測(cè)結(jié)果中存在多處滯后現(xiàn)象，預(yù)測(cè)效果最差。在四步預(yù)測(cè)步長(zhǎng)下，SSA-LSSVM方法相對(duì)誤差始終維持在一個(gè)較低水平，說明組合模型能有效提升機(jī)器學(xué)習(xí)方法對(duì)時(shí)間序列的擬合水平。SSA方法能夠很好地挖掘原始時(shí)間序列中的關(guān)鍵特征，而LSSVM方法可以模擬序列分解得到非線性特性變化。

圖6給出了不同預(yù)測(cè)方法在各預(yù)測(cè)步長(zhǎng)下的誤差性能指標(biāo)，SSA-LSSVM預(yù)測(cè)模型在不同預(yù)測(cè)步長(zhǎng)下的預(yù)測(cè)誤差均小于其他方法，隨著預(yù)測(cè)步長(zhǎng)的提升，同一預(yù)測(cè)模型的預(yù)測(cè)性能也會(huì)逐漸下降。單一LSSVM預(yù)測(cè)性能最差，相較于傳統(tǒng)EMD方法，采用EEMD-LSSVM和ICEEMD-LSSVM方法有效降低了模型預(yù)測(cè)誤差，但是對(duì)預(yù)測(cè)性能的改善作用十分有限，隨著預(yù)測(cè)步長(zhǎng)的增大，這一改善作用會(huì)逐漸減弱。在各預(yù)測(cè)步長(zhǎng)下，基于SSA的預(yù)測(cè)方法的NRMSE和NMAE始終低于EMD類方法誤差性能指標(biāo)，說明與EMD類方法相比，SSA方法對(duì)序列關(guān)鍵特征分量提取能力更強(qiáng)，輸出分量質(zhì)量較高，有利于LSSVM學(xué)習(xí)序列特征，提升組合模型的預(yù)測(cè)精度。

圖6 預(yù)測(cè)誤差性能指標(biāo)Fig.6 Prediction error performance index

5 結(jié)論

針對(duì)風(fēng)電功率序列復(fù)雜的波動(dòng)特性，本文提出了一種基于SSA的超短期風(fēng)電功率組合多步預(yù)測(cè)方法。算例對(duì)比研究結(jié)果表明：

①通過Cao方法確定SSA的最佳嵌入維度，有效避免了人為因素的影響以及單純追求預(yù)測(cè)結(jié)果而造成的序列特征被淹沒和計(jì)算量暴增等問題；

②SSA具有直接提取風(fēng)電功率序列不同成分的能力，與其他數(shù)據(jù)預(yù)處理方法相比，分解后序列平均樣本熵最大，有效保留了原始序列關(guān)鍵特征；

③所提基于SSA的預(yù)測(cè)模型能夠獲取時(shí)間序列中的時(shí)序演變特征和波動(dòng)趨勢(shì)，有效提升了不同步長(zhǎng)下的預(yù)測(cè)精度，計(jì)算效率更高。