本發(fā)明涉及風(fēng)電場(chǎng)功率預(yù)測(cè)技術(shù)領(lǐng)域，特別涉及一種應(yīng)用于新建風(fēng)電場(chǎng)的短期功率預(yù)測(cè)方法。

背景技術(shù)：

近年來(lái)，隨著我國(guó)風(fēng)力發(fā)電迅速發(fā)展，短時(shí)間內(nèi)大量新建風(fēng)電場(chǎng)并網(wǎng)投產(chǎn)，這些風(fēng)電場(chǎng)缺乏準(zhǔn)確可靠的歷史運(yùn)行數(shù)據(jù)，無(wú)法通過(guò)傳統(tǒng)的統(tǒng)計(jì)建模方法建立高精度短期功率預(yù)測(cè)模型，極大影響了新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)精度。為確保大規(guī)模風(fēng)電安全、可靠、高效的并網(wǎng)運(yùn)行，開(kāi)發(fā)出滿足工程應(yīng)用的風(fēng)電功率預(yù)測(cè)系統(tǒng)非常關(guān)鍵。

技術(shù)實(shí)現(xiàn)要素：

有鑒于此，本發(fā)明的主要目的在于，提供一種應(yīng)用于新建風(fēng)電場(chǎng)的短期功率預(yù)測(cè)方法，包括步驟：

A、利用BP神經(jīng)網(wǎng)絡(luò)，建立新建風(fēng)電場(chǎng)區(qū)域中尺度數(shù)值天氣預(yù)報(bào)預(yù)測(cè)因子與風(fēng)電場(chǎng)前期測(cè)風(fēng)塔實(shí)測(cè)風(fēng)速之間的歷史數(shù)據(jù)關(guān)系，生成測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度的風(fēng)速統(tǒng)計(jì)降尺度模型；

B、依據(jù)風(fēng)電場(chǎng)區(qū)域中尺度數(shù)值天氣預(yù)報(bào)的預(yù)測(cè)因子，以及步驟A中所述測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度統(tǒng)計(jì)降尺度模型，生成測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度的短期預(yù)測(cè)風(fēng)速；

C、依據(jù)風(fēng)速-功率特性曲線解析數(shù)據(jù)以及風(fēng)電場(chǎng)綜合廠用電耗損比例，建立新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型，計(jì)算出新建風(fēng)電場(chǎng)短期預(yù)測(cè)功率。

由上，通過(guò)對(duì)中尺度模式數(shù)值天氣預(yù)報(bào)進(jìn)行統(tǒng)計(jì)降尺度，得出測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度的短期預(yù)測(cè)風(fēng)速，進(jìn)而建立新建風(fēng)電場(chǎng)短期功 率預(yù)測(cè)模型，實(shí)現(xiàn)對(duì)新建風(fēng)電場(chǎng)整體出力短期預(yù)測(cè)。有效減少了中尺度分辨率不足帶來(lái)的不確定性，顯著提高了新建風(fēng)電場(chǎng)短期風(fēng)速預(yù)測(cè)精度，有效解決了新建風(fēng)電場(chǎng)因缺乏風(fēng)機(jī)運(yùn)行數(shù)據(jù)無(wú)法建立高精度功率預(yù)測(cè)模型的難題。

可選的，步驟A中，所述中尺度數(shù)值天氣預(yù)報(bào)的預(yù)測(cè)因子至少包括：500hPa位勢(shì)高度、850hPa位勢(shì)高度、風(fēng)機(jī)輪轂高度的風(fēng)速、風(fēng)向、壓強(qiáng)和相對(duì)濕度。

由上，由于上述預(yù)報(bào)因子可以被較為準(zhǔn)確的預(yù)測(cè)，且不同預(yù)測(cè)因子間是弱相關(guān)或無(wú)關(guān)的，因此將上述預(yù)報(bào)因子作為BP神經(jīng)網(wǎng)絡(luò)的輸入量，可以較為精確的對(duì)BP神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練。

可選的，步驟A中，還包括對(duì)中尺度數(shù)值天氣預(yù)報(bào)的預(yù)測(cè)因子進(jìn)行歸一化處理。

可選的，步驟C中所述風(fēng)速-功率特性曲線解析數(shù)據(jù)包括：

采用多階高斯函數(shù)擬合風(fēng)速-功率特性曲線，得到風(fēng)速-功率特性曲線的解析函數(shù)P(v)，

P(v)＝a1*exp(-((v-b1)/c1)^2)+a2*exp(-((v-b2)/c2)^2)+a3*exp(-((v-b3)/c3)^2)+a4*exp(-((v-b4)/c4)^2)，

式中v表示風(fēng)速值，ai，bi，ci(i＝1～4)表示多階高斯擬合曲線的常值系數(shù)。

可選的，步驟C中所述風(fēng)電場(chǎng)綜合廠用電耗損比例＝(發(fā)電量—上網(wǎng)電量)/發(fā)電量＝綜合廠用電量/發(fā)電量；

綜合廠用電量＝風(fēng)電場(chǎng)內(nèi)部用電量+主變及場(chǎng)內(nèi)線路損耗電量+電抗器用電量。

可選的，所述新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型為N*P(v)*C；

其中N表示風(fēng)電場(chǎng)風(fēng)機(jī)數(shù)量，P(v)表示風(fēng)速-功率特性曲線的解析數(shù)據(jù)，C表示風(fēng)電場(chǎng)綜合廠用電耗損比例。

可選的，步驟A中，所述數(shù)值天氣預(yù)報(bào)預(yù)測(cè)因子的歷史預(yù)測(cè)數(shù)據(jù)、與中尺度數(shù)值天氣預(yù)報(bào)預(yù)測(cè)因子時(shí)間匹配的風(fēng)電場(chǎng)前期測(cè)風(fēng)塔實(shí)測(cè)風(fēng)速的時(shí)間長(zhǎng)度不少于6個(gè)自然月。

可選的，步驟C后還包括：結(jié)合實(shí)際功率對(duì)新建風(fēng)電場(chǎng)短期預(yù)測(cè)功率進(jìn)行校驗(yàn)，依據(jù)校驗(yàn)結(jié)果對(duì)所述新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型進(jìn)行修正的步驟。

附圖說(shuō)明

圖1為本發(fā)明流程圖；

圖2為本發(fā)明BP神經(jīng)網(wǎng)絡(luò)的原理示意圖；

圖3風(fēng)機(jī)的原廠風(fēng)速-功率特性曲線示意圖。

具體實(shí)施方式

為克服現(xiàn)有技術(shù)存在的缺陷，本發(fā)明提供一種應(yīng)用于新建風(fēng)電場(chǎng)的短期功率預(yù)測(cè)方法，通過(guò)對(duì)中尺度模式預(yù)測(cè)輸出進(jìn)行降尺度運(yùn)算，得出風(fēng)電場(chǎng)前期測(cè)風(fēng)塔所在位置輪轂高度的短期預(yù)測(cè)風(fēng)速，進(jìn)而依據(jù)風(fēng)速-功率特性曲線解析函數(shù)和風(fēng)電場(chǎng)綜合廠用電耗損比例建立新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型，實(shí)現(xiàn)對(duì)風(fēng)電場(chǎng)整體出力的短期預(yù)測(cè)。上述預(yù)測(cè)方法有效減少了中尺度分辨率不足帶來(lái)的不確定性，以及新建風(fēng)電場(chǎng)缺少大量歷史運(yùn)行數(shù)據(jù)導(dǎo)致無(wú)法建立高精度統(tǒng)計(jì)模型的缺點(diǎn)，顯著提高了新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)精度。

如圖1所示，一種應(yīng)用于新建風(fēng)電場(chǎng)的短期功率預(yù)測(cè)方法包括：

步驟10：收集風(fēng)電場(chǎng)建立前測(cè)風(fēng)塔的實(shí)測(cè)風(fēng)速數(shù)據(jù)，以及與之時(shí)間匹配的風(fēng)電場(chǎng)區(qū)域中尺度數(shù)值天氣預(yù)報(bào)預(yù)測(cè)因子的歷史預(yù)測(cè)數(shù)據(jù)。

中尺度數(shù)值天氣預(yù)報(bào)預(yù)測(cè)因子的選擇是應(yīng)用統(tǒng)計(jì)降尺度方法過(guò)程中重要的環(huán)節(jié),預(yù)測(cè)因子的選擇很大程度上決定了風(fēng)電場(chǎng)局地氣象條件的預(yù)報(bào)特征。選取對(duì)風(fēng)速有顯著影響的預(yù)測(cè)因子，減少其數(shù)量，可以有效降低統(tǒng)計(jì)降尺度模型的復(fù)雜程度，減少模型計(jì)算量，避免引入額外的干擾信息。選取預(yù)測(cè)因子遵循以下方法：

1、預(yù)測(cè)因子必須能夠被中尺度數(shù)值天氣預(yù)報(bào)模式較準(zhǔn)確地預(yù)測(cè)，并且預(yù)測(cè)因子與風(fēng)電場(chǎng)內(nèi)部站點(diǎn)氣象要素之間具有顯著的非線性統(tǒng)計(jì)關(guān)系，并且這種統(tǒng)計(jì)關(guān)系是穩(wěn)定而有效的；

2、預(yù)測(cè)因子必須能夠反映重要的中尺度大氣物理變化過(guò)程；

3、不同預(yù)測(cè)因子間是弱相關(guān)或無(wú)關(guān)的。

基于上述幾點(diǎn)，所選取的預(yù)測(cè)因子分別為：中尺度數(shù)值天氣預(yù)報(bào)所提供的500hPa位勢(shì)高度、850hPa位勢(shì)高度、風(fēng)機(jī)輪轂高度的風(fēng)速、風(fēng)向、壓強(qiáng)和相對(duì)濕度。

所述風(fēng)電場(chǎng)前期測(cè)風(fēng)塔歷史數(shù)據(jù)時(shí)間長(zhǎng)度至少為半年，較佳的，為一年。

步驟20：使用BP神經(jīng)網(wǎng)絡(luò)建立風(fēng)電場(chǎng)前期測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度的統(tǒng)計(jì)降尺度模型。

使用BP神經(jīng)網(wǎng)絡(luò)建立統(tǒng)計(jì)降尺度模型，具體來(lái)說(shuō)，即中尺度數(shù)值天氣預(yù)報(bào)預(yù)測(cè)因子與測(cè)風(fēng)塔位置輪轂高度風(fēng)速之間的統(tǒng)計(jì)關(guān)系，將此統(tǒng)計(jì)關(guān)系作為測(cè)風(fēng)塔的統(tǒng)計(jì)降尺度模型應(yīng)用于產(chǎn)生該測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度預(yù)測(cè)風(fēng)速。

BP神經(jīng)網(wǎng)絡(luò)是一種按誤差逆向傳播算法訓(xùn)練的多層前饋網(wǎng)絡(luò)，能以任意精度逼近任何非線性映射，其拓?fù)浣Y(jié)構(gòu)由輸入層、隱含層和輸出層組成，該模型輸出為：式中，y_k(t)為模型輸出的t時(shí)刻第k臺(tái)風(fēng)機(jī)輪轂高度風(fēng)速預(yù)測(cè)值；f為傳遞函數(shù)，本實(shí)施例中，傳遞函數(shù)為正切雙曲函數(shù)；w_j為連接隱含層和輸出層的權(quán)重系數(shù)；x_i(t)為模型輸入t時(shí)刻第i個(gè)預(yù)測(cè)因子值；v_ij為連接輸入層和隱含層的權(quán)重系數(shù)；n為輸入層維數(shù)；m為隱含層的維數(shù)；θ_j為隱含層閾值；θ₀為輸出層閾值。

以建立測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度統(tǒng)計(jì)降尺度模型為例具體說(shuō)明，如圖2所示，建模數(shù)據(jù)為2012年1月至2012年6月的中尺度數(shù)值天氣預(yù)報(bào)模式預(yù)測(cè)因子和測(cè)風(fēng)塔實(shí)測(cè)風(fēng)速數(shù)據(jù)，數(shù)據(jù)時(shí)間分辨率為15min，數(shù)據(jù)長(zhǎng)度選取次日0～24h，將模式輸出的500hPa位勢(shì)高度、850hPa位勢(shì)高度、風(fēng)速V、風(fēng)向D正弦值、風(fēng)向D余弦值、大氣壓強(qiáng)P、濕度RH作為模型訓(xùn)練輸入量，風(fēng)機(jī)輪轂高度風(fēng)速實(shí)測(cè)值作為訓(xùn)練輸出量，建立三層的BP神經(jīng)網(wǎng)絡(luò)模型。該BP神經(jīng)網(wǎng)絡(luò)輸入層神經(jīng)元數(shù)為6；隱含層神經(jīng)元數(shù)經(jīng)過(guò)測(cè)試擇優(yōu)確定；網(wǎng)絡(luò)輸出層神經(jīng) 元數(shù)為1。訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)時(shí)，對(duì)輸入和輸出層數(shù)據(jù)進(jìn)行歸一化處理，這里采用如下方法進(jìn)行歸一化處理：以壓強(qiáng)為例進(jìn)行說(shuō)明，X表示預(yù)測(cè)壓強(qiáng)值，X_min表示預(yù)設(shè)的最小壓強(qiáng)值，X_max表示預(yù)設(shè)的最大壓強(qiáng)值，表示歸一化的壓強(qiáng)值。

經(jīng)過(guò)篩選分析，確定隱含層神經(jīng)元數(shù)為21時(shí)，訓(xùn)練樣本誤差最小(此步驟為BP正常流程)，此時(shí)各權(quán)重系數(shù)矩陣和閾值矩陣也隨之確定。權(quán)重值和閾值矩陣系數(shù)確定后的BP神經(jīng)網(wǎng)絡(luò)模型作為測(cè)風(fēng)塔的統(tǒng)計(jì)降尺度預(yù)測(cè)模型。

步驟30：依據(jù)測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度統(tǒng)計(jì)降尺度模型和中尺度數(shù)值天氣預(yù)報(bào)提供的風(fēng)電場(chǎng)區(qū)域短期預(yù)測(cè)風(fēng)速，生成測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度短期預(yù)測(cè)風(fēng)速。

完成測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度統(tǒng)計(jì)降尺度建模后，將中尺度數(shù)值天氣預(yù)報(bào)預(yù)測(cè)因子(500hPa位勢(shì)高度、850hPa位勢(shì)高度、風(fēng)機(jī)輪轂高度風(fēng)速、風(fēng)向、壓強(qiáng)和相對(duì)濕度)輸入測(cè)風(fēng)塔統(tǒng)計(jì)降尺度模型，便可輸出測(cè)風(fēng)塔位置風(fēng)機(jī)輪轂高度的短期預(yù)測(cè)風(fēng)速。

中尺度數(shù)值天氣預(yù)報(bào)統(tǒng)計(jì)降尺度方法，使用上述6個(gè)預(yù)測(cè)因子便可計(jì)算出測(cè)風(fēng)塔位置輪轂高度的高分辨率短期風(fēng)速預(yù)測(cè)值，計(jì)算量小且計(jì)算速度快，完全可以滿足風(fēng)電場(chǎng)短期功率預(yù)測(cè)工程要求，該方法相比于現(xiàn)有動(dòng)力降尺度技術(shù)，大大提高計(jì)算效率。

步驟40：建立風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型。

如圖3所示，風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型基于原廠風(fēng)速-功率特性曲線數(shù)據(jù)P(v)和風(fēng)電場(chǎng)綜合廠用電耗損比例C構(gòu)建。采用高階高斯函數(shù)對(duì)原廠風(fēng)速-功率特性曲線數(shù)據(jù)進(jìn)行非線性擬合，獲取原廠風(fēng)速-功率特性曲線的解析函數(shù)P(v)；風(fēng)電場(chǎng)綜合廠用電耗損因素包含風(fēng)電場(chǎng)內(nèi)部線路耗損、變壓器耗損和風(fēng)電場(chǎng)自用電等。

具體的，采用多階高斯函數(shù)擬合原廠風(fēng)速-功率特性曲線，得到原廠風(fēng)速-功率特性曲線的解析函數(shù)，該解析函數(shù)具有如下的表達(dá)式：

P(v)＝a1*exp(-((v-b1)/c1)^2)+a2*exp(-((v-b2)/c2)^2)+a3*exp(-((v-b3)/c3)^2)+a4*exp(-((v-b4)/c4)^2)

式中v表示風(fēng)速值，ai，bi，ci(i＝1～4)表示多階高斯擬合曲線的常值系數(shù)。

另外，風(fēng)電場(chǎng)綜合廠用電耗損比例C＝(發(fā)電量—上網(wǎng)電量)/發(fā)電量＝綜合廠用電量/發(fā)電量；

綜合廠用電量＝風(fēng)電場(chǎng)內(nèi)部用電量+主變及場(chǎng)內(nèi)線路損耗電量+電抗器用電量。

通過(guò)上述計(jì)算式，可以得出發(fā)生在風(fēng)電場(chǎng)內(nèi)部的電量耗損比列。

新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型采用如下表達(dá)，風(fēng)電場(chǎng)短期預(yù)測(cè)功率＝N*P(v)*C，式中N為風(fēng)電場(chǎng)風(fēng)機(jī)數(shù)量。

步驟50：對(duì)風(fēng)電場(chǎng)短期功率進(jìn)行預(yù)測(cè)。

將經(jīng)過(guò)統(tǒng)計(jì)降尺度處理的風(fēng)電場(chǎng)前期測(cè)風(fēng)塔位置的短期預(yù)測(cè)風(fēng)速輸入到新建風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型，獲得風(fēng)電場(chǎng)短期預(yù)測(cè)功率。

步驟60：結(jié)合實(shí)際功率，對(duì)步驟S50中的預(yù)測(cè)數(shù)據(jù)進(jìn)行校驗(yàn)，并依據(jù)校驗(yàn)結(jié)果對(duì)步驟40中的風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型進(jìn)行修正。

對(duì)于新建風(fēng)電場(chǎng)，其實(shí)際功率數(shù)據(jù)較少，則需要依據(jù)較少的實(shí)際功率數(shù)據(jù)對(duì)于預(yù)測(cè)值進(jìn)行準(zhǔn)確度判斷。舉例來(lái)說(shuō)，當(dāng)預(yù)測(cè)功率小于實(shí)際功率時(shí)，則對(duì)應(yīng)減小風(fēng)電場(chǎng)內(nèi)部的電量耗損比列C，重新計(jì)算預(yù)測(cè)功率。反之，當(dāng)預(yù)測(cè)功率大于實(shí)際功率時(shí)，則對(duì)應(yīng)增加風(fēng)電場(chǎng)內(nèi)部的電量耗損比列C，直至所計(jì)算出的預(yù)測(cè)功率與實(shí)際功率相近，確定出最接近真實(shí)值的風(fēng)電場(chǎng)內(nèi)部的電量耗損比列C。從而實(shí)現(xiàn)對(duì)于步驟40中風(fēng)電場(chǎng)短期功率預(yù)測(cè)模型的修正。

以上所述僅為本發(fā)明的較佳實(shí)施例而已，并不用以限制本發(fā)明?？傊?，凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。