一種光伏系統(tǒng)輸出功率預(yù)測(cè)方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及光伏發(fā)電領(lǐng)域，特別是涉及一種光伏系統(tǒng)輸出功率預(yù)測(cè)方法。
【背景技術(shù)】
[0002] 光伏發(fā)電是一個(gè)非線性的隨機(jī)過程，具有顯著的間歇性和隨機(jī)性，易受天氣變化、 日夜交替、季節(jié)推移等因素的影響，光伏發(fā)電一旦大規(guī)模入網(wǎng)，必將導(dǎo)致電網(wǎng)的運(yùn)行方式、 配網(wǎng)結(jié)構(gòu)和潮流方向等發(fā)生重大變化，因此我們?cè)谘芯亢凸夥l(fā)電的同時(shí)，如何去保證其 大規(guī)模入網(wǎng)后，電網(wǎng)依然能夠穩(wěn)定、健康的運(yùn)行，將成為一項(xiàng)非常重要也非常值得去研究的 課題。隨著光伏發(fā)電系統(tǒng)容量的不斷擴(kuò)大，光伏發(fā)電預(yù)測(cè)技術(shù)對(duì)于減輕光伏發(fā)電系統(tǒng)輸出 電能的隨機(jī)性對(duì)電力系統(tǒng)的影響具有重要意義。目前常見的光伏系統(tǒng)輸出功率預(yù)測(cè)方法有 BP神經(jīng)預(yù)測(cè)，灰色預(yù)測(cè)、粒子群優(yōu)化算法和遺傳算法組合預(yù)測(cè)等，這類方法雖然避免了預(yù)測(cè) 模型的復(fù)雜建模，但依賴于輸出功率的歷史數(shù)據(jù)，對(duì)原始數(shù)據(jù)樣本的選擇具有更高的要求。
[0003] 高斯過程（Gaussianprocess，GP)是近年來基于統(tǒng)計(jì)學(xué)理論發(fā)展而來的一種全新 學(xué)習(xí)機(jī)，它對(duì)處理高維數(shù)、小樣本、非線性等復(fù)雜分類和回歸問題具有很好的適應(yīng)性，且泛 化能力強(qiáng)，在不犧牲性能的條件下，與人工神經(jīng)網(wǎng)絡(luò)和支持向量機(jī)相比有著容易實(shí)現(xiàn)的特 點(diǎn)；同時(shí)，其超參數(shù)可通過求取訓(xùn)練樣本的對(duì)數(shù)似然函數(shù)的極大值自適應(yīng)獲取，有著靈活的 非參數(shù)推斷和預(yù)測(cè)輸出的概率解釋，是一個(gè)具有概率意義的核學(xué)習(xí)機(jī)。

【發(fā)明內(nèi)容】

[0004] 針對(duì)現(xiàn)有技術(shù)上存在的不足，本發(fā)明的目的是基于高斯過程對(duì)光伏系統(tǒng)的輸出功 率進(jìn)行合理的預(yù)測(cè)，公開一種光伏系統(tǒng)輸出功率預(yù)測(cè)方法。
[0005] 為了實(shí)現(xiàn)上述目的，本發(fā)明是通過如下的技術(shù)方案來實(shí)現(xiàn)：
[0006] -種光伏系統(tǒng)輸出功率預(yù)測(cè)方法，該方法包括如下步驟：
[0007] (1)、選預(yù)測(cè)日前一天的輸出功率歷史數(shù)據(jù)、全天的平均溫度、平均輻照度、預(yù)測(cè)日 的平均溫度和輻照度組成原始的訓(xùn)練樣本集&，
[0008] (2)、構(gòu)建&的滿足高斯分布隨機(jī)過程的函數(shù)t，預(yù)測(cè)日相應(yīng)的實(shí)測(cè)數(shù)據(jù)作為測(cè)試 樣本I，預(yù)測(cè)日滿足高斯分布隨機(jī)過程的函數(shù)的實(shí)測(cè)數(shù)據(jù)作為I；，利用貝葉斯推理算法，建 立t和I；的分布函數(shù)關(guān)系，
[0009] (3)、選擇組合核函數(shù)作為高斯方法中的協(xié)方差函數(shù)，通過對(duì)訓(xùn)練樣本的對(duì)數(shù)似然 函數(shù)的極大化獲得協(xié)方差函數(shù)的最優(yōu)超參數(shù)，進(jìn)一步可確定預(yù)測(cè)函數(shù)的均值和方差，求得 輸出功率的預(yù)測(cè)函數(shù)，實(shí)現(xiàn)光伏系統(tǒng)發(fā)電量的預(yù)測(cè)。
[0010] 上述步驟（2)中建立t和T*的分布函數(shù)關(guān)系的步驟如下：
[0011] 將預(yù)測(cè)日前一天的輸出功率歷史數(shù)據(jù)、全天的平均溫度、平均輻照度、預(yù)測(cè)日的平 均溫度和輻照度，記尤=，…，，％；，xei ,&，_%j，設(shè)共獲得N個(gè)訓(xùn)練樣本，即該數(shù)據(jù) 集為
[0012]D={XoUXj|i= 1，2,…，N} (1)
[0013] 將輸出功率、溫度等數(shù)據(jù)X通過高斯函數(shù)映射為t，假設(shè)
[0014] tx=f(X,)+￡,,! =l,2,-,n (2)
[0015] 其中，X為N維輸入矢量，q為輸出標(biāo)量，f(X)服從高斯分布的隨機(jī)過程，即：
[0016] f(X)~GP(y(X)，C(X，X'）） (3)
[0017] y(X)為高斯過程的均值，變量X'GX為N維輸入矢量，C(X，X'）為對(duì)稱正定協(xié)方 差矩陣，其第i行，第j列的元素記為Clj，表示輸出功率訓(xùn)練樣本X和X'之間的相關(guān)性，其 值通過正定核函數(shù)k(X，X'）計(jì)算獲得，e為獨(dú)立的隨機(jī)變量，符合高斯分布， 〇 "為方差；
[0018] 在貝葉斯線性回歸f(X)=巾⑴Tw框架下，采用參數(shù)向量w的隨機(jī)分布w~ N(0,A)，再由公式（3)得公式（2)對(duì)應(yīng)的輸出功率數(shù)據(jù)集t滿足如下分布：
[0019]
(4；)
[0020]記預(yù)測(cè)日實(shí)測(cè)數(shù)據(jù)為待測(cè)試數(shù)據(jù)為f由公式（4)推出訓(xùn)練樣本t和測(cè)試樣本T$， 所形成的聯(lián)合高斯先驗(yàn)分布為
[0021]
[0022] C(X，X)為nXn階對(duì)稱正定的協(xié)方差矩陣，C(X，X，為測(cè)試f與訓(xùn)練的所有點(diǎn)輸入 點(diǎn)X的nX1階協(xié)方差矩陣，CU' 測(cè)試點(diǎn)f自身的協(xié)方差，I為單位矩陣。
[0023] 上述步驟（3)中輸出功率的預(yù)測(cè)函數(shù)如下：
[0024] 在給定測(cè)試點(diǎn)訓(xùn)練集D的條件下，貝葉斯概率預(yù)測(cè)的目標(biāo)是計(jì)算概率 IA/)，預(yù)測(cè)日的輸出功率f滿足如下分布：
[0025]
C6)
[0026]f的期望和方差分別為：
[0029] 確定了f的期望和方差，就得到預(yù)測(cè)日的輸出功率，T為矩陣的轉(zhuǎn)置，C(x'X)為 測(cè)試x#與訓(xùn)練的所有點(diǎn)輸入點(diǎn)X的協(xié)方差矩陣。
[0030] 上述步驟（3)中選擇組合核函數(shù)作為高斯方法中的協(xié)方差函數(shù)，通過對(duì)訓(xùn)練樣本 的對(duì)數(shù)似然函數(shù)的極大化獲得協(xié)方差函數(shù)的最優(yōu)超參數(shù)的計(jì)算方法如下：
[0031] 由于高斯方法中協(xié)方差函數(shù)在有限輸入點(diǎn)集上要求是正定的，且是一個(gè)滿足 Mercer條件的對(duì)稱函數(shù)，故協(xié)方差函數(shù)等價(jià)于核函數(shù)，公式（7)改寫成如下形式：
[0032]
(9)：
[0033] 預(yù)測(cè)值的均值是核函數(shù)C的線性組合，將非線性關(guān)系的數(shù)據(jù)映射到特征空間后轉(zhuǎn) 換為線性關(guān)系，從而使復(fù)雜非線性問題轉(zhuǎn)化為容易處理的線性問題；
[0034] 高斯過程選擇不同的協(xié)方差函數(shù)，3種單一協(xié)方差函數(shù)：
[0035]平方指數(shù)協(xié)方差函數(shù)
[0041 ] 令0=G/V/丨，〇_丨，〇:.:0為包含所有超參數(shù)的向量；_ =diagd2)為超參數(shù)的對(duì)稱矩 陣；h參數(shù)為核函數(shù)的信號(hào)方差，用來控制局部相關(guān)性的程度；s克洛內(nèi)克爾符號(hào)，其值 越大，表示輸入與輸出相關(guān)性越??；a為核函數(shù)的形狀參數(shù)；〇n為方差，通過對(duì)訓(xùn)練樣本 的對(duì)數(shù)似然函數(shù)的極大化獲得最優(yōu)超參數(shù)9。
[0042] 上述對(duì)數(shù)似然函數(shù)的極大化獲得最優(yōu)超參數(shù)0的計(jì)算方法如下：
[0043]
[0045] 其中，％麵為對(duì)數(shù)似然函數(shù)"為矩陣的轉(zhuǎn)置。
[0046] 上述步驟（3)中的組合核函數(shù)考慮到當(dāng)f\(X)，fJX)，…，fN⑴都是相互獨(dú)立的高 斯隨機(jī)過程時(shí)，隨機(jī)過程
也為一個(gè)高斯過程，故采用公式（15)的核函數(shù)建 立高斯過程回歸模型：
[0049] 式中：K為組合核函數(shù)；a=QH，Q(X，X'）為X和X'的協(xié)方差矩陣；
[0050] 組合核函數(shù)，由公式（11)和公式（12)組合而成，即：
[0051]CdXj) =CRQ+CM〇
[0052] 本發(fā)明的實(shí)施過程簡明，利用輸出功率的歷史數(shù)據(jù)對(duì)光伏系統(tǒng)的輸出功率進(jìn)行預(yù) 測(cè)，避免創(chuàng)建復(fù)雜的預(yù)測(cè)模型，降低了預(yù)測(cè)的難度，有助于電網(wǎng)電量的調(diào)度。
【附圖說明】
[0053] 下面結(jié)合附圖和【具體實(shí)施方式】來詳細(xì)說明本發(fā)明；
[0054] 圖1是一種光伏系統(tǒng)輸出功率預(yù)測(cè)方法流程圖。
【具體實(shí)施方式】
[0055] 為使本發(fā)明實(shí)現(xiàn)的技術(shù)手段、創(chuàng)作特征、達(dá)成目的與功效易于明白了解，下面結(jié)合
【具體實(shí)施方式】，進(jìn)一步闡述本發(fā)明。
[0056] 如圖1所示，一種光伏系統(tǒng)輸出功率預(yù)測(cè)方法，包括如下步驟：
[0057](1)、選擇預(yù)測(cè)日前一天的輸出功率歷史數(shù)據(jù)、全天的平均溫度、平均輻照度、預(yù)測(cè) 日的平均溫度和輻照度組成原始的訓(xùn)練樣本集&，
[0058] (2)、構(gòu)建&的滿足高斯分布隨機(jī)過程的函數(shù)t，預(yù)測(cè)日滿足高斯分布隨機(jī)過程的 函數(shù)的實(shí)測(cè)數(shù)據(jù)作為!；，利用貝葉斯推理算法，建立t和I；的分布函數(shù)關(guān)系；
[0059] 將預(yù)測(cè)日前一天的輸出功率歷史數(shù)據(jù)、全天的平均溫度、平均輻照度、預(yù)測(cè)日的平 均溫度和輻照度，記^，％2，? ，設(shè)共獲得N個(gè)訓(xùn)練樣本，即該數(shù)據(jù)集 為
[0060] D= ?"$)|i= 1，2,…，N} (1)
[0061] 將輸出功率、溫度等數(shù)據(jù)X通過高斯函數(shù)映射為t，假設(shè)
[0062] tx=f(X,)+￡,,! =l,2,-,n (2)
[0063] 其中，X為N維輸入矢量，q為輸出標(biāo)量，f(X)服從高斯分布的隨機(jī)過程，即： f(X)~GP(y(X)，C(X，X'）） (3)
[0064] y(X)為高斯過程的均值，變量X'GX為N維輸入矢量，C(X，X'）為對(duì)稱正定協(xié)方 差矩陣，其第i行，第j列的元素記為Clj，表示輸出功率訓(xùn)練樣本X和X'之間的相關(guān)性，其 值通過正定核函數(shù)k(X，X'）計(jì)算獲得。e為獨(dú)立的隨機(jī)變量，符合高斯分布，g~/V(0.大）， 〇 "為方差。
[0065] 在貝葉斯線性回歸f(X)=<HX)Tw框架下，采用參數(shù)向量w的隨機(jī)分布w~ N(0,A)，再由公式（3)得公式（2)對(duì)應(yīng)的輸出功率數(shù)據(jù)集t滿足如下分布：
[0066]
(4)
[