國電聯(lián)合動力技術(shù)有限公司 王果毅 牟金磊 李潤祥

風(fēng)資源評估及建模過程是決定宏觀區(qū)域選址、機型選擇、電網(wǎng)運行調(diào)度、容量可信度計算等工作的決定性因素。在風(fēng)資源評估過程中完整年風(fēng)速是各項評估的基礎(chǔ)依據(jù)。完整年風(fēng)速數(shù)據(jù)通常通過建立測風(fēng)塔進行采集。而測風(fēng)塔的精密部件大多長期暴露于自然環(huán)境中，測風(fēng)塔的損壞會導(dǎo)致風(fēng)速數(shù)據(jù)的缺失。

基于歷史實測數(shù)據(jù)提取信息所擬合的風(fēng)速概率模型為基礎(chǔ)，通過對該概率模型進行抽樣，對風(fēng)速缺失部分進行補全是處理風(fēng)速缺失問題的常用方法。在眾多的概率密度模型中，威布爾分布以其形式簡單、描述靈活等優(yōu)點為業(yè)內(nèi)普遍接受。威布爾分布中包含兩個關(guān)鍵參數(shù)，通常以K 和C 表示，兩參數(shù)將會在威布爾分布的仿真過程中共同影響風(fēng)速分布模型的分布情況及峰值位置。在實際操作中，威布爾分布的兩參數(shù)有多種求法，常用方法有最大風(fēng)速法、累積分布函數(shù)擬合、平均風(fēng)速法、平均風(fēng)速與標(biāo)準(zhǔn)差等。

縱觀上述兩參數(shù)的求法，都是通過以歷史實測風(fēng)速的年平均值及風(fēng)速標(biāo)準(zhǔn)差為基礎(chǔ)計算得到的，因此威布爾分布在對區(qū)域風(fēng)場的風(fēng)速分布概率的描述中有很好的表達效果。但是傳統(tǒng)的兩參數(shù)威布爾分布在描繪風(fēng)速分布時，并沒有考慮風(fēng)速序列相鄰數(shù)值間的相互影響關(guān)系，在仿真過程中僅能參考風(fēng)速分布密度得到無序的概率序列，與實際工況中存在變化規(guī)律的風(fēng)速序列存在較大差異。因此傳統(tǒng)兩參數(shù)威布爾分布在對風(fēng)資源的評估中依舊存在一定局限性。

若將風(fēng)速視為變量，變量間的相互影響關(guān)系通常以相關(guān)性進行描述。在眾多的相關(guān)性描述方法中Copula 函數(shù)形式靈活、應(yīng)用廣泛，在不同領(lǐng)域內(nèi)如水文、氣候、經(jīng)濟等方面均有良好的應(yīng)用效果，因此在描述風(fēng)速序列相鄰數(shù)值間的相互影響關(guān)系上同樣具備參考價值。

本文提供一種基于時間相關(guān)性的風(fēng)速描述方法，在傳統(tǒng)兩參數(shù)風(fēng)速威布爾分布的基礎(chǔ)上，結(jié)合基于Copula 函數(shù)對風(fēng)資源信息中的時間相關(guān)性信息進行提取，并基于時間相關(guān)性在基于威布爾分布對風(fēng)速抽樣過程進行約束。仿真證明，該方法可為風(fēng)資源評估過程中的風(fēng)速缺失部分提供更符合自然工況的仿真數(shù)據(jù)，對威布爾分布的建模過程有顯著的優(yōu)化效果，在實際工程中具備良好的應(yīng)用價值。

威布爾分布及Copula 函數(shù)

風(fēng)速威布爾分布

威布爾分布是在風(fēng)電行業(yè)內(nèi)被廣泛認(rèn)可的風(fēng)速密度函數(shù)，其主要形式如式（1）所示。

式中變量x 為風(fēng)速；c 為描述風(fēng)速的尺度參數(shù)；k 為描述風(fēng)速分布的形狀參數(shù)。C 值越大，威布爾分布所描述的風(fēng)速均值越大；k 值越大，威布爾分布描述的風(fēng)速分布越集中。關(guān)鍵參數(shù)有多種求法，其中較為常用的一種如式（2）所示。

式中α 為0 到1 之間均勻分布的隨機變量，x 為由α 及關(guān)鍵參數(shù)k 值及c 值所得到的模擬風(fēng)速數(shù)值。在仿真過程中通過對α 進行隨機抽樣來建立風(fēng)速序列。如圖1 所示，該圖為基于某風(fēng)場威布爾分布所模擬的間隔為10 分鐘、時長為一周的仿真風(fēng)速，該序列所包含的風(fēng)速數(shù)據(jù)有1008 個。該風(fēng)速的關(guān)鍵參數(shù)K=2.6、C=7.45。

如圖2 所示，該圖為該風(fēng)場的歷史實測風(fēng)速。

通過對比可以明顯看出，相比于具有明顯連續(xù)性的實測風(fēng)速，基于威布爾分布所擬合的風(fēng)速呈現(xiàn)明顯的排列混亂且波動劇烈。因此，僅基于威布爾分布對風(fēng)速進行模擬來補充風(fēng)速的缺失部分存在明顯的局限性。

Copula 函數(shù)

Copula 函數(shù)的主要功能，是反向求取多元隨機變量的相關(guān)性關(guān)系，其表達形式如式（4）所示。

F（x1，x2，…，xn）＝C（F1（x1），F(xiàn)2（x2），…，F(xiàn)n（xn）） （4）

式中F（x1，x2，…，xn）為隨機變量的聯(lián)合分布；Fn（xn）為聯(lián)合分布中第n 個變量的邊緣分布；C（F1（x1），F(xiàn)2（x2），…，F(xiàn)n（xn））為用以表達變量間相關(guān)性的Copula 函數(shù)。Copula 函數(shù)的實際原理是通過擬合隨機變量的分布得到不同變量0 到1 之間均勻分布的特征值，以二元Copula 函數(shù)為例，式（5）為兩個隨機變量。

式中X1、X2為需要尋求相關(guān)性的隨機變量，其中將x1、x2為[0.1]均勻分布的隨機數(shù)，分別是X1、X2的特征值。當(dāng)X1、X2表示風(fēng)速時，結(jié)合式（3）可得，x1、x2分別是不同風(fēng)速的α 值，f（x）、g（x）分別為不同關(guān)鍵參數(shù)威布爾分布的反函數(shù)。Copula函數(shù)相關(guān)性建模原理如式（6）所示。

式中K 為x1、x2基于Copula 函數(shù)相關(guān)性特性系數(shù)，x3為[0.1]均勻分布的隨機數(shù)。通過對隨機變量的特征值間建立函數(shù)關(guān)系，再由具備函數(shù)關(guān)系的特征值可以得到具備相關(guān)性的變量值。

由此可見，應(yīng)用Copula 函數(shù)對隨機變量相關(guān)性建模時充分地考慮了兩變量之間的相關(guān)性及變量各自的隨機性。

基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法

根據(jù)上一節(jié)中的分析，基于威布爾分布的風(fēng)速建模方法的局限性，主要在于對變量特征值α 的采樣為隨機抽樣，忽略了風(fēng)速序列相鄰數(shù)值間的相互影響關(guān)系。

由Copula 函數(shù)對數(shù)值相關(guān)性的靈活且良好表達效果，提出基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法。通過Copula 函數(shù)對風(fēng)速數(shù)值間的相關(guān)性進行發(fā)掘，并結(jié)合Copula函數(shù)對風(fēng)速建模過程進行優(yōu)化。

圖1 基于威布爾分布的模擬風(fēng)速

圖2 歷史實測風(fēng)速

假設(shè)有一組風(fēng)速，序列中包含（v1-vn）n 個風(fēng)速數(shù)據(jù)，將風(fēng)速序列需包含的數(shù)據(jù)如式（7）所示進行組合。

將V1 與V2 中的元素一一對應(yīng)，就可以構(gòu)成如式（8）所示的對應(yīng)關(guān)系。

由式（8）所示的對應(yīng)關(guān)系所得到的相關(guān)性關(guān)系可以認(rèn)為是風(fēng)速序列相鄰點數(shù)值之間的相關(guān)性，即風(fēng)速的時間相關(guān)性。對于相關(guān)性強弱的評估，通常以相關(guān)性系數(shù)作為評估標(biāo)準(zhǔn)，常用的相關(guān)性系數(shù)有kendall 相關(guān)性系數(shù)、spearman 相關(guān)性系數(shù)、pearson 相關(guān)性系數(shù)。由表1 所示，該表為實測數(shù)據(jù)及威布爾分布擬合的風(fēng)速序列在不同類型相關(guān)性系數(shù)。

表1 不同類型相關(guān)性對比

根據(jù)表1 可以看出，實測風(fēng)速具有極強的時間相關(guān)性，不同的相關(guān)性系數(shù)均達到0.8 以上，而基于威布爾分布的模擬風(fēng)速相鄰數(shù)據(jù)之間毫無相關(guān)性可言。因此提出基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法，應(yīng)用實測風(fēng)速的強相關(guān)性對基于威布爾分布的風(fēng)速建模過程進行約束。如式（9）～（11）所示。

式（9）對應(yīng)式（1），vn為風(fēng)速序列中最后一個數(shù)據(jù)，基于威布爾分布求取該風(fēng)速所對應(yīng)的特征值αn。

式（10）對應(yīng)式（6），基于αn及風(fēng)速的時間相關(guān)性得到vn下一點需要擬合的風(fēng)速vn+1所對應(yīng)的特征值αn+1。

式（11）對應(yīng)式（3），應(yīng)用威布爾分布的反函數(shù)及由式（10）所得到的vn+1所對應(yīng)的特征值αn+1共同得到風(fēng)速vn+1的實際數(shù)值。

此時得到的vn+1將參考風(fēng)速，根據(jù)式（9）～（11）求取vn+2。以此步驟進行循環(huán)，直至將風(fēng)速的缺失部分填充完整。

綜上所述，基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法的流程如圖3 所示。

仿真對比驗證

模擬風(fēng)序列的特性對比

圖3 基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法流程

圖4 基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法模擬風(fēng)速

表2 不同類型相關(guān)性對比

應(yīng)用基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法同樣擬合，同樣以10 分鐘為間隔，時間跨度為一周，數(shù)據(jù)量為1008。模擬風(fēng)速如圖4 所示。

如式（7），將應(yīng)用三種不同相關(guān)性系數(shù)對基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法的模擬風(fēng)速分為兩個相鄰風(fēng)速序列，同樣應(yīng)用三種不同的相關(guān)性對兩序列進行相關(guān)性評價，如表2 所示。

根據(jù)表2 所示的基于風(fēng)速時間相關(guān)性的模擬風(fēng)速時間相關(guān)性系數(shù)，對比表1 的兩種不同風(fēng)速的相關(guān)性可以看到，基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法所得到的模擬風(fēng)速的時間相關(guān)性更接近自然條件下實測風(fēng)速的時間相關(guān)性。同樣對比圖4 及圖1 可以看出，在時間相關(guān)性的約束下，所建立的風(fēng)速具備更好的連續(xù)性及更具規(guī)律的波動性。

圖5 基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法模擬風(fēng)速

如圖5 所示，為基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法的模擬風(fēng)速對分布概率進行擬合的分布曲線，及實測風(fēng)速對分布概率進行擬合的分布曲線進行比較。

圖中虛線為實測風(fēng)速的分布曲線，實線為基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法的模擬風(fēng)速的分布曲線。

通過圖5 的對比可以得到結(jié)論，基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法在對風(fēng)速模擬過程優(yōu)化的同時，依舊能夠良好地保持風(fēng)速的分布情況，該方法不僅對威布爾分布有顯著的優(yōu)化效果，同時對風(fēng)速分布有良好的還原效果。模擬風(fēng)序列的出力對比

風(fēng)機出力模型如式（12）所示。

表中，V 表示風(fēng)速；Vci為風(fēng)機切入風(fēng)速；Vco為風(fēng)機切出風(fēng)速；Vr表示額定風(fēng)速；Pr表示風(fēng)機額定功率；P（V）表示風(fēng)速為V時對應(yīng)的風(fēng)機輸出功率。風(fēng)機工作中，風(fēng)機功率輸出如圖6 所示。

圖6 風(fēng)機功率輸出曲線

由式（12）可知，風(fēng)機出力與風(fēng)速呈現(xiàn)三次方關(guān)系，風(fēng)速的劇烈波動將會被急劇放大。如圖7、圖8 所示，圖7 為基于威布爾分布模擬風(fēng)速的風(fēng)電場出力曲線，圖8 為基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法的模擬風(fēng)速的風(fēng)電場出力曲線。該風(fēng)場擬采用20 臺2.5 兆瓦風(fēng)機。風(fēng)機參數(shù)見表3。

圖7 基于威布爾分布模擬風(fēng)速的風(fēng)電場出力曲線

表3 風(fēng)機參數(shù)表

通過圖7、圖8 的對比可以得到結(jié)論，基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法模擬風(fēng)速在應(yīng)用于風(fēng)電場出力計算的工作中具備更好的連續(xù)性，可以更好地還原風(fēng)電場的出力模式，相比于基于威布爾分布的風(fēng)速模擬結(jié)果擁有更好的工程應(yīng)用價值。

結(jié)論

圖8 基于風(fēng)速時間相關(guān)性模擬風(fēng)速的風(fēng)電場出力曲線

風(fēng)資源評估是風(fēng)電工程的基礎(chǔ)，而完整的風(fēng)速數(shù)據(jù)是風(fēng)資源評估工作的保障。因此，對缺失風(fēng)速數(shù)據(jù)進行良好的補充，對風(fēng)電工程的順利開展具有重要意義。

基于威布爾分布的風(fēng)速模擬在工程應(yīng)用中依舊存在一定的缺陷，最明顯的就是在對風(fēng)速的模擬過程中未對風(fēng)速序列的時間相關(guān)性進行充分考慮，導(dǎo)致所得到的風(fēng)速序列波動過大、連續(xù)性差，而這個問題在對風(fēng)電場出力評估的過程中將被劇烈放大。

為改善威布爾分布在風(fēng)速模擬工作中的缺陷，本文提出一種基于風(fēng)速時間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法，主要原理為應(yīng)用Copula 函數(shù)對風(fēng)速的時間相關(guān)性進行發(fā)掘，并用風(fēng)速的時間相關(guān)性對風(fēng)速模擬過程進行約束。經(jīng)過仿真驗證表明，該方法對風(fēng)速模擬工作有巨大的改良效果，同時對風(fēng)速分布依舊有良好的還原效果，可以得到還原性更好的風(fēng)速仿真數(shù)據(jù)，更符合實際工程需求，在對風(fēng)速理論改進的同時擁有良好的工程應(yīng)用前景。