劉麗珺，梁友嘉

(1.甘肅省氣象服務(wù)中心，甘肅　蘭州　730020；2.武漢理工大學(xué)資源與環(huán)境工程學(xué)院，湖北　武漢　430070)

?

BJ-RUC模式的低層風(fēng)場(chǎng)預(yù)報(bào)產(chǎn)品質(zhì)量檢驗(yàn)
——以甘肅省當(dāng)金山風(fēng)電場(chǎng)為例

劉麗珺1，梁友嘉2

(1.甘肅省氣象服務(wù)中心，甘肅蘭州730020；2.武漢理工大學(xué)資源與環(huán)境工程學(xué)院，湖北武漢430070)

利用2013年甘肅華電阿克塞當(dāng)金山49.5 mW風(fēng)電場(chǎng)風(fēng)機(jī)監(jiān)控?cái)?shù)據(jù)采集系統(tǒng)的觀測(cè)數(shù)據(jù)，對(duì)BJ-RUC模式的低層風(fēng)場(chǎng)預(yù)報(bào)產(chǎn)品進(jìn)行質(zhì)量檢驗(yàn)。結(jié)果表明：(1)總體上有效風(fēng)速區(qū)在各季節(jié)預(yù)測(cè)效果較滿發(fā)風(fēng)速區(qū)好，且這2個(gè)風(fēng)速區(qū)秋季預(yù)報(bào)效果最好，冬季預(yù)報(bào)效果最差；(2)年尺度上，觀測(cè)值的Weibull分布眾值偏右，形態(tài)相對(duì)接近正態(tài)分布，而預(yù)測(cè)值的Weibull分布眾值偏左，形態(tài)相對(duì)偏離正態(tài)分布。月尺度上，除1月外，其余各月實(shí)況和預(yù)測(cè)的Weibull分布與全年的基本一致，且各月預(yù)測(cè)風(fēng)速形狀參數(shù)k值接近2，為中性風(fēng)，說明BJ-RUC模式模擬強(qiáng)陣性風(fēng)速有局限性；(3)模擬的主導(dǎo)風(fēng)向與實(shí)測(cè)的風(fēng)向有較大偏差，僅對(duì)接近滿發(fā)風(fēng)速區(qū)的風(fēng)向模擬得較準(zhǔn)確，表明BJ-RUC模式對(duì)風(fēng)向模擬有待提高。

風(fēng)速；風(fēng)向；BJ-RUC模式；誤差分析

引　言

隨機(jī)變化的低層風(fēng)場(chǎng)通常會(huì)導(dǎo)致電廠風(fēng)電輸出功率呈現(xiàn)出波動(dòng)性和間歇性的特點(diǎn)[1]。風(fēng)電場(chǎng)輸出功率的預(yù)測(cè)研究已成為提高電網(wǎng)中風(fēng)電并網(wǎng)能力的重要技術(shù)支撐和研究熱點(diǎn)[2]。風(fēng)電功率預(yù)測(cè)通常以數(shù)值天氣預(yù)報(bào)(Numerical Weather Prediction，NWP)模式為基礎(chǔ)，并與物理、統(tǒng)計(jì)和耦合預(yù)測(cè)算法集成，最終預(yù)測(cè)風(fēng)電場(chǎng)輸出功率。NWP模式的數(shù)據(jù)質(zhì)量和不確定性顯著影響并制約風(fēng)電功率的預(yù)測(cè)精度，亟需對(duì)NWP模式數(shù)據(jù)進(jìn)行科學(xué)評(píng)估和檢驗(yàn)，這也是提高不同時(shí)間尺度風(fēng)電功率預(yù)報(bào)精度的難點(diǎn)[3]。

目前，國(guó)際上可用于風(fēng)電場(chǎng)功率預(yù)測(cè)的經(jīng)典NWP模式主要有：MetUM(Met Office Unified Model)、GEM(Global Environmental Multiscale Model)、ECMWF(European Centre for Medium-Range Weather Forecasts)、WRF(Weather Research Forecasting)、MM5(Fifth-Generation Penn State/NCAR Mesoscale Model)、RUC(Rapid Update Cycle)和COSMO(Consortium for Small Scale Modeling)等[4]。近年來，國(guó)內(nèi)相關(guān)研究也發(fā)展迅速，其中，北京市氣象局在RUC模式基礎(chǔ)上開發(fā)的北京快速循環(huán)更新同化系統(tǒng)(Beijing Rapid Update Cycle，BJ-RUC)已廣泛應(yīng)用于國(guó)內(nèi)低層風(fēng)場(chǎng)的風(fēng)電功率預(yù)測(cè)[5-6]。相應(yīng)地，針對(duì)多種低層風(fēng)場(chǎng)數(shù)值模擬產(chǎn)品的檢驗(yàn)方法也開始迅速發(fā)展，如：Powers[7]對(duì)6種主流中尺度模式的模擬值與觀測(cè)值分別進(jìn)行風(fēng)速偏差(Mean Error，ME)、平均絕對(duì)誤差(Mean Absolute Error，MAE)和均方根誤差(Root Mean Square Error，RMSE)分析，系統(tǒng)評(píng)價(jià)了各模式數(shù)據(jù)預(yù)測(cè)質(zhì)量；Jiménez等[8]對(duì)復(fù)雜地形下的低層風(fēng)場(chǎng)進(jìn)行2 km高分辨率的數(shù)值模擬，并使用標(biāo)準(zhǔn)化Taylor方法對(duì)各站點(diǎn)RMSE及相關(guān)系數(shù)進(jìn)行檢驗(yàn)；Storm等[9]通過擬合的威布爾(Weibull)風(fēng)速概率密度函數(shù)對(duì)比分析了低空風(fēng)切變條件下的風(fēng)速模擬精度；李軍等[10]利用雙參數(shù)Weibull分布模型分析內(nèi)蒙古草原6層高度近1 a的風(fēng)速序列資料，并利用最大似然法估計(jì)多時(shí)段、多高度Weibull分布尺度參數(shù)c和形狀參數(shù)k；Wang等[11]基于WRF模式預(yù)報(bào)的甘肅酒泉電場(chǎng)區(qū)域風(fēng)速和風(fēng)向預(yù)報(bào)產(chǎn)品，分析了西北地區(qū)風(fēng)場(chǎng)預(yù)報(bào)中的誤差特征及其影響因素。本文利用2013年甘肅阿克塞當(dāng)金山風(fēng)電場(chǎng)風(fēng)機(jī)實(shí)況數(shù)據(jù)，對(duì)BJ-RUC模式的低層風(fēng)場(chǎng)數(shù)值預(yù)報(bào)產(chǎn)品進(jìn)行評(píng)估檢驗(yàn)，并探討B(tài)J-RUC模式在干旱區(qū)開展風(fēng)預(yù)報(bào)應(yīng)用的優(yōu)勢(shì)和局限性，可為電廠提高風(fēng)電功率產(chǎn)出和掌握區(qū)域風(fēng)速特征和風(fēng)能資源狀況提供科學(xué)參考，并為BJ-RUC模式區(qū)域性應(yīng)用和改進(jìn)提供案例參考。

1　數(shù)據(jù)與方法

1.1數(shù)據(jù)

實(shí)況數(shù)據(jù)源于2013年甘肅華電阿克塞當(dāng)金山49.5 mW風(fēng)電場(chǎng)風(fēng)機(jī)監(jiān)控?cái)?shù)據(jù)采集系統(tǒng)(Supervisory Control and Data Acquisition，SCADA)，該系統(tǒng)由中國(guó)華電集團(tuán)福新能源股份有限公司甘肅分公司提供。風(fēng)電場(chǎng)位于甘肅省酒泉市阿克塞哈薩克族自治縣(簡(jiǎn)稱阿克塞縣)南部約70 km處的高原戈壁，場(chǎng)址北側(cè)為當(dāng)金山，三面環(huán)山。場(chǎng)區(qū)海拔在3 050～3 200 m之間，地形平坦、地勢(shì)開闊，附近有215國(guó)道，施工交通便利。風(fēng)功率模擬數(shù)據(jù)由BJ-RUC模式獲取，BJ-RUC主模式為WRF V2.2，前處理模塊為WPS V2.2和WRF-VAR V2.1，模擬時(shí)間為2013年全年。模式同化和預(yù)報(bào)區(qū)域設(shè)置為27 km、9 km、3 km三重嵌套，其中西北地區(qū)選用9 km分辨率數(shù)據(jù)，模式模擬范圍為96.45°E—127.53°E、37.20°N—54.48°N，垂直方向38層，模式層頂為50 hPa。模式參數(shù)化方案為：WSM6顯式微物理方案、Kain-Fritseh(new Eta)對(duì)流參數(shù)化方案、YSU邊界層方案、RRTM長(zhǎng)波輻射方案、Goddard短波輻射方案、Noah陸面模式[12]。模式采用全球環(huán)流模式背景場(chǎng)資料(T639的GFS資料)進(jìn)行再分析獲得模式初值，并將常規(guī)及加密探空觀測(cè)、常規(guī)及加密地面觀測(cè)、船舶/浮標(biāo)觀測(cè)、飛機(jī)觀測(cè)等全球觀測(cè)資料作為業(yè)務(wù)運(yùn)行同化資料[13]。

通過BJ-RUC將WRF數(shù)值預(yù)報(bào)模式產(chǎn)品9 km分辨率提高至3 km，并進(jìn)一步利用CALMET(California Meteorological Model)動(dòng)力降尺度模式將分辨率再次提高到1 km，然后插值到每個(gè)風(fēng)機(jī)的輪轂高度(80 m)，并輸出未來72 h內(nèi)每15 min格點(diǎn)的風(fēng)速、風(fēng)向、垂直風(fēng)速、溫度、相對(duì)濕度和壓強(qiáng)等氣象要素資料(圖1)。

圖1　技術(shù)路線

1.2基礎(chǔ)數(shù)據(jù)檢驗(yàn)

參考《風(fēng)電場(chǎng)風(fēng)能資源評(píng)估方法》(GB/T 18710-2002)，數(shù)據(jù)檢驗(yàn)主要關(guān)注數(shù)據(jù)的完整性和合理性。數(shù)據(jù)完整性是指其數(shù)量和時(shí)間順序是否完整。其中，數(shù)據(jù)數(shù)量應(yīng)等于預(yù)期記錄的數(shù)量，而時(shí)間順序應(yīng)符合預(yù)期的開始和結(jié)束時(shí)間，且時(shí)間區(qū)間應(yīng)連續(xù)。利用測(cè)風(fēng)有效數(shù)據(jù)完整率指標(biāo)檢驗(yàn)，風(fēng)電廠所需有效數(shù)據(jù)完整率一般需>90%，具體公式如下：

(1)

式(1)中，Rvi為有效數(shù)據(jù)完整率，No為應(yīng)測(cè)數(shù)目(即15 min測(cè)量間隔數(shù))，Nm為缺測(cè)數(shù)目(即未記錄的15 min間隔平均值數(shù)目)，Ni為無效數(shù)據(jù)數(shù)目(即不合理的15 min間隔平均值數(shù)目)。合理性檢驗(yàn)主要是根據(jù)數(shù)據(jù)的合理性范圍來檢驗(yàn)數(shù)據(jù)。通常，風(fēng)速范圍為0～40 m·s-1，風(fēng)向范圍為0°～360°。數(shù)據(jù)初始檢驗(yàn)后，對(duì)缺測(cè)和不合理的數(shù)據(jù)進(jìn)行二次判別，進(jìn)一步甄別符合實(shí)際情況的有效數(shù)據(jù)，并將其加入原始數(shù)據(jù)組。并分析備用和參考傳感器的同期記錄數(shù)據(jù)，用于填補(bǔ)無效或缺測(cè)數(shù)據(jù)。

表1是當(dāng)金山風(fēng)電場(chǎng)2013年逐月測(cè)風(fēng)有效數(shù)據(jù)的完整率，通過檢驗(yàn)其各月的實(shí)況數(shù)據(jù)，剔除并替換無效數(shù)據(jù)，最終獲取研究區(qū)實(shí)測(cè)的逐15 min風(fēng)速風(fēng)向數(shù)據(jù)。

1.3誤差分析

在數(shù)據(jù)檢驗(yàn)基礎(chǔ)上，進(jìn)一步提高BJ-RUC模式模擬的NWP數(shù)據(jù)的準(zhǔn)確性和可靠性，通過多尺度數(shù)據(jù)檢驗(yàn)方法實(shí)現(xiàn)：(1)誤差指數(shù)：包括標(biāo)準(zhǔn)化平均絕對(duì)誤差(NMAE)、標(biāo)準(zhǔn)化均方根誤差(NRMSE)、平均絕對(duì)百分比誤差(MAPE)和相關(guān)系數(shù)(R)[14]，其中NMAE、NRMSE和MAPE越低則表明預(yù)報(bào)效果越好，R越高表明模擬值與實(shí)況值的相關(guān)性越好；(2)誤差分布：包括風(fēng)速頻率威布爾分布和風(fēng)向風(fēng)速玫瑰圖。根據(jù)上述指標(biāo)分析風(fēng)速區(qū)間和年、月尺度下模擬與實(shí)測(cè)數(shù)據(jù)的誤差分布特征。

表1　2013年各月當(dāng)金山風(fēng)電場(chǎng)測(cè)風(fēng)有效數(shù)據(jù)的完整性

1.3.1誤差指數(shù)

風(fēng)機(jī)運(yùn)轉(zhuǎn)對(duì)風(fēng)速有嚴(yán)格要求。風(fēng)速<3.0 m·s-1時(shí)，風(fēng)機(jī)難以運(yùn)轉(zhuǎn)；風(fēng)速在3.0～11.5 m·s-1區(qū)間內(nèi)，風(fēng)機(jī)輸出功率隨風(fēng)速增大而相應(yīng)增大；風(fēng)速>11.5 m·s-1后，風(fēng)機(jī)處于滿發(fā)狀態(tài)，其輸出功率不再隨風(fēng)速增加而增大；風(fēng)速≥25.0 m·s-1時(shí)，往往易造成風(fēng)機(jī)結(jié)構(gòu)損壞[15]。故參考上述標(biāo)準(zhǔn)，將風(fēng)速劃分為4個(gè)等級(jí)區(qū)間：低風(fēng)速區(qū)(0～3.5 m·s-1)、有效風(fēng)速區(qū)(3.5～11.5 m·s-1)、滿發(fā)風(fēng)速區(qū)(11.5～25.0 m·s-1)和高風(fēng)速區(qū)(≥25.0 m·s-1)。分析不同風(fēng)速區(qū)的預(yù)報(bào)特性及誤差分布規(guī)律，有針對(duì)性地提高風(fēng)功率預(yù)測(cè)精度。

為更好地描述風(fēng)速預(yù)報(bào)效果，首先，分別計(jì)算MAE、RMSE和R。其中，對(duì)MAE和RMSE進(jìn)行標(biāo)準(zhǔn)化處理以消除量綱影響，從而得到NMAE和NRMSE。其次，按照《中華人民共和國(guó)能源行業(yè)標(biāo)準(zhǔn)(NB/T 31046-2013)：風(fēng)電功率預(yù)測(cè)系統(tǒng)功能規(guī)范》要求：“單個(gè)風(fēng)電場(chǎng)短期預(yù)測(cè)月均方根誤差率應(yīng)<20%，超短期預(yù)測(cè)第4 h預(yù)測(cè)值月均方根誤差率應(yīng)<15%”進(jìn)行檢驗(yàn)和分析。為準(zhǔn)確評(píng)價(jià)預(yù)報(bào)在不同時(shí)段的效果，還需計(jì)算平均絕對(duì)百分比誤差(MAPE)。

1.3.2誤差分布

雙參數(shù)威布爾(Weibull)分布函數(shù)f(v)是一單峰的正偏態(tài)分布函數(shù)，更接近風(fēng)速實(shí)際分布[16]，其風(fēng)速v的函數(shù)式(3)。式(3)中：f(v)為風(fēng)速分布函數(shù)；k為無量綱形狀參數(shù)，反映風(fēng)速分布情況。k值越大，說明風(fēng)速分布越集中，表示風(fēng)的變化越小，陣性越弱；k值越小，說明風(fēng)速分布越分散，表示風(fēng)的變化越大，陣性更強(qiáng)；c為尺度因子(m·s-1)，其與平均風(fēng)速相關(guān)，平均風(fēng)速越大，c值越大；μ為估算的風(fēng)速平均值(m·s-1)；σ為估算的風(fēng)速標(biāo)準(zhǔn)差(m·s-1)；Γ(1+1/k)為伽馬函數(shù)。

(2)

(3)

此外，風(fēng)向風(fēng)速玫瑰圖是一種重要的圖形分析工具，可用于與風(fēng)速因子密切相關(guān)的氣象預(yù)報(bào)、城市規(guī)劃、建筑規(guī)劃和風(fēng)資源評(píng)估等研究中[17]。風(fēng)向風(fēng)速玫瑰圖基于網(wǎng)格化極坐標(biāo)系統(tǒng)制圖，通過不同方位匯集風(fēng)向及其頻數(shù)，頻數(shù)最大方位表示該風(fēng)向出現(xiàn)次數(shù)最多，用不同顏色帶區(qū)別風(fēng)速大小。

2　結(jié)果與分析

2.1不同等級(jí)風(fēng)速區(qū)的誤差

由四季不同等級(jí)風(fēng)速區(qū)間的頻率分布(圖2)發(fā)現(xiàn)：各季中有效風(fēng)速區(qū)間所占比例最大，且夏、秋季較高，分別占全季的83.43%和89.89%；低風(fēng)速區(qū)間次之，春、冬季較高，分別占全季的32.74%和33.85%；滿發(fā)和高風(fēng)速區(qū)間所占比例整體較小。其中，滿發(fā)風(fēng)速夏季占比最高為1.09%，秋季最低為0.28%；高風(fēng)速秋季占比最高為1.28%，冬季最低為0.94%。

圖2　2013年當(dāng)金山風(fēng)電場(chǎng)四季不同等級(jí)風(fēng)速頻率(單位：%)

由于高、低風(fēng)速區(qū)不能發(fā)電，這里僅分析有效風(fēng)速區(qū)和滿發(fā)風(fēng)速區(qū)各季節(jié)風(fēng)預(yù)報(bào)效果(圖3)。從季節(jié)來看，實(shí)況與預(yù)報(bào)值之間的相關(guān)系數(shù)R在2個(gè)風(fēng)速區(qū)都表現(xiàn)為秋季>冬季>春季>夏季。其中，在有效風(fēng)速區(qū)實(shí)況與預(yù)報(bào)值四季均為正相關(guān)，而在滿發(fā)風(fēng)速區(qū)二者秋、冬季表現(xiàn)為正相關(guān)，春、夏季表現(xiàn)為負(fù)相關(guān)；除了冬季滿發(fā)風(fēng)速區(qū)誤差較高外，其余季節(jié)2個(gè)風(fēng)速區(qū)的NMAE和NRMSE都為0.5～1.0 m·s-1，且自秋、夏、春季依次降低；2個(gè)風(fēng)速區(qū)的MAPE季節(jié)分布略有不同。在有效風(fēng)速區(qū)，MAPE表現(xiàn)為冬季>夏季>秋季>春季，而在滿發(fā)風(fēng)速區(qū)則表現(xiàn)為冬季>夏季>春季>秋季。對(duì)比2個(gè)風(fēng)速區(qū)的誤差和相關(guān)系數(shù)發(fā)現(xiàn)，平均R，有效區(qū)>滿發(fā)區(qū)；MAPE，除夏季外，有效區(qū)<滿發(fā)區(qū)；NMAE和NRMSE，有效區(qū)<滿發(fā)區(qū)。綜上分析可見，有效風(fēng)速在各季節(jié)的預(yù)測(cè)效果較滿發(fā)風(fēng)速好，且秋季預(yù)報(bào)效果總體較好，冬季預(yù)報(bào)效果最差。

2.2不同時(shí)間尺度誤差

2.2.1年尺度

年尺度誤差分析能夠體現(xiàn)BJ-RUC模式綜合預(yù)報(bào)效果。由圖4可知，當(dāng)金山風(fēng)電廠2013年全風(fēng)速區(qū)BJ-RUC模式預(yù)報(bào)值與實(shí)況的NMAE為0.51 m·s-1，NRMSE為0.65 m·s-1，R為0.25，MAPE為39%，接近可發(fā)電風(fēng)速區(qū)的誤差結(jié)果，說明該模式產(chǎn)品總體質(zhì)量較優(yōu)。

圖3　BJ-RUC模式預(yù)報(bào)的當(dāng)金山風(fēng)電場(chǎng)2013年四季有效風(fēng)速區(qū)和滿發(fā)風(fēng)速區(qū)與實(shí)況之間的R、NMAE、NRMSE和MAPE

圖4　2013年BJ-RUC模式預(yù)報(bào)的當(dāng)金山風(fēng)電場(chǎng)風(fēng)速與實(shí)況的NMAE、NRMSE、R及MAPE

根據(jù)全年實(shí)況(Obs)和預(yù)測(cè)(SL)輪轂高度風(fēng)速估計(jì)的Weibull曲線(圖5)得到：觀測(cè)的全年輪轂高度c_Obs為5.78，k_Obs為2.84；預(yù)測(cè)的全年輪轂高度c_SL為6.33，k_SL為1.74。年實(shí)況Weibull分布眾值偏右，形態(tài)相對(duì)接近正態(tài)分布；年預(yù)測(cè)Weibull分布眾值偏左，形態(tài)相對(duì)偏離正態(tài)分布。

分析當(dāng)金山風(fēng)電場(chǎng)2013年風(fēng)機(jī)輪轂高度處的風(fēng)向玫瑰圖(圖6)發(fā)現(xiàn)：該電廠區(qū)觀測(cè)的主導(dǎo)風(fēng)向?yàn)槟?S)和南南東(SSE)，其頻率分別達(dá)26.5%和22%；模擬的主導(dǎo)風(fēng)向?yàn)槟夏蠔|(SSE)和北北西(NNW)，與實(shí)測(cè)有較大偏差，該模式對(duì)風(fēng)向預(yù)報(bào)效果需提高。

2.2.2月尺度

總體上，各月實(shí)況(Obs)與預(yù)測(cè)(SL)風(fēng)速數(shù)據(jù)的Weibull曲線(圖7)與全年大致相同。除1月外，各月實(shí)況Weibull分布眾值偏右，相對(duì)接近正態(tài)分布；預(yù)測(cè)值分布的眾值偏左，相對(duì)偏離正態(tài)分布。各月k_Obs在2.29～4.33之間，其中1—4月均<3，4月最小，5—12月均>3，9月最大；c_Obs在4.16～7.06之間，其中6月最大，12月最小?？傮w上，春、冬季c_Obs和k_Obs比夏、秋季小。各月c_SL在4.17～8.09之間，3月最大，12月最?。籯_SL在1.28～2.12之間，12月最小，8月最大，1—10月均接近2。根據(jù)《風(fēng)電場(chǎng)風(fēng)能資源評(píng)估方法》(GB/T 18710-2002)：當(dāng)k=2時(shí)，Weibull分布變?yōu)镽ayleigh分布。由于各月k_SL在2.0左右，故風(fēng)的變化一般，陣性強(qiáng)度中等，說明該模式未模擬出陣性較強(qiáng)的風(fēng)速序列，預(yù)測(cè)較實(shí)況風(fēng)速變化平緩。

圖5　2013年當(dāng)金山風(fēng)電場(chǎng)風(fēng)速實(shí)況與BJ-RUC模式模擬的威布爾風(fēng)頻圖

圖6　2013年當(dāng)金山風(fēng)電場(chǎng)的實(shí)況與BJ-RUC模式模擬的風(fēng)向風(fēng)速(單位：m·s-1)玫瑰圖

圖8是當(dāng)金山風(fēng)電廠2013年各月風(fēng)機(jī)輪轂高度處風(fēng)向風(fēng)速玫瑰圖，發(fā)現(xiàn)該廠區(qū)觀測(cè)的各月主導(dǎo)風(fēng)向均為偏南風(fēng)(S)，而預(yù)測(cè)的主導(dǎo)風(fēng)向除1月為東風(fēng)(E)外，其余月份均為南東(SE)、南西(SW)、北東(NE)和北西(NW)方向，說明該模式對(duì)月尺度風(fēng)向預(yù)報(bào)效果仍需提高。鑒于觀測(cè)的全年主導(dǎo)風(fēng)向均為偏南方向，應(yīng)調(diào)整模式模擬風(fēng)向。

2.2.3日尺度

波動(dòng)性和隨機(jī)性導(dǎo)致風(fēng)速很難準(zhǔn)確預(yù)報(bào)，在大尺度環(huán)流背景相對(duì)穩(wěn)定時(shí)，風(fēng)的局地性使其日變化可能較大。因此，進(jìn)一步計(jì)算日尺度模擬風(fēng)速與實(shí)況的NRMSE、MAPE和R。通過對(duì)當(dāng)金山風(fēng)電場(chǎng)2013年8月1—31日風(fēng)速模擬與實(shí)況間的誤差和相關(guān)系數(shù)(表2)分析發(fā)現(xiàn)：1個(gè)月中2/3的天數(shù)的NRMSE<20m·s-1，MAPE<25%，基本滿足規(guī)范中精度要求；相關(guān)系數(shù)R多為負(fù)值，表明BJ-RUC模式還需優(yōu)化參數(shù)化方案，以提高相關(guān)性，同時(shí)也說明此方法對(duì)模擬結(jié)果有較好的檢驗(yàn)效果。

圖7　2013年當(dāng)金山風(fēng)電場(chǎng)風(fēng)速逐月實(shí)況與BJ-RUC模式模擬的威布爾風(fēng)頻圖

日期NRMSE/m·s-1MAPE/%R日期NRMSE/m·s-1MAPE/%R12.2119.050.35173.0826.04-0.6821.6521.160.07182.4123.350.3331.9323.420.77192.7326.28-0.4643.4230.01-0.21202.8633.41-0.5555.7440.36-0.65213.1029.760.2763.1230.510.34222.5322.40-0.2071.9222.830.16232.2024.93-0.5183.5550.560.73243.4940.94-0.5992.1720.89-0.11253.0528.820.70102.9426.700.19262.8721.730.04111.2320.250.28274.5943.86-0.51122.0023.920.84281.7930.080.81132.0818.950.63292.7039.33-0.80143.7020.54-0.56302.7532.720.33151.9813.43-0.15311.7415.95-0.31162.9530.30-0.44

圖8　2013年當(dāng)金山風(fēng)電場(chǎng)的逐月實(shí)況(Obs.1-12)與BJ-RUC模式模擬(SL.1-12)的風(fēng)向風(fēng)速(單位：m·s-1)玫瑰圖

3　結(jié)論與討論

(1)當(dāng)金山電場(chǎng)2013年各季中有效風(fēng)速區(qū)所占比例最大，低風(fēng)速區(qū)次之，滿發(fā)和高風(fēng)速區(qū)所占比例整體較小。

(2)BJ-RUC模式模擬與觀測(cè)的有效風(fēng)速區(qū)和滿發(fā)風(fēng)速區(qū)的NMAE、NRMSE和MAPE值在秋季最小、冬季最大，且有效風(fēng)速區(qū)的誤差小于滿發(fā)風(fēng)速區(qū)，表明秋季2個(gè)風(fēng)速區(qū)的預(yù)報(bào)效果優(yōu)于其他季節(jié)，且各季節(jié)有效風(fēng)速區(qū)的預(yù)測(cè)效果整體較好。

(3)年尺度上，實(shí)況的Weibull分布眾值偏右，形態(tài)相對(duì)接近正態(tài)分布，而預(yù)測(cè)的Weibull分布眾值偏左，形態(tài)相對(duì)偏離正態(tài)分布；月尺度上，除1月外，其余各月觀測(cè)與模擬的Weibull分布與全年基本一致。風(fēng)速的c和k都有明顯的季節(jié)變化：總體上，春、冬季觀測(cè)的c_Obs值比夏、秋季小，而k_Obs值秋季最大，春季最??；預(yù)測(cè)的各月風(fēng)速估計(jì)k值接近2，較觀測(cè)值偏低，為中性風(fēng)，說明BJ-RUC模式模擬的風(fēng)速月變化較為平緩，觀測(cè)值的曲線較觀測(cè)更趨于對(duì)稱，對(duì)強(qiáng)陣性風(fēng)速預(yù)測(cè)有局限性。

(4)模式預(yù)測(cè)的風(fēng)機(jī)主導(dǎo)風(fēng)向與實(shí)測(cè)風(fēng)向有較大偏差，但接近滿發(fā)風(fēng)速區(qū)的風(fēng)向模擬較為準(zhǔn)確。

低層風(fēng)場(chǎng)既受天氣形勢(shì)變化影響，也受復(fù)雜下墊面影響，從而形成特殊的局地環(huán)流。BJ-RUC數(shù)值模式的邊界條件、物理參數(shù)化方案以及觀測(cè)數(shù)據(jù)等方面都存在不確定性。首先，BJ-RUC模式的預(yù)報(bào)中心定位為北京地區(qū)，西北干旱區(qū)與前者的地理位置、氣候條件和下墊面條件均不一致，在模式的邊界條件、物理參數(shù)化方案沒有進(jìn)一步本地化的情況下進(jìn)行預(yù)報(bào)勢(shì)必會(huì)影響預(yù)報(bào)精度；其次，當(dāng)金山風(fēng)電場(chǎng)無測(cè)風(fēng)塔，觀測(cè)數(shù)據(jù)來源于風(fēng)機(jī)輪轂高度處SCADA系統(tǒng)采集的數(shù)據(jù)，風(fēng)機(jī)數(shù)據(jù)受尾流影響較大，從而影響了模式的對(duì)比精度。BJ-RUC能夠快速更新循環(huán)同化最新的探測(cè)資料，其優(yōu)勢(shì)使其預(yù)報(bào)產(chǎn)品廣泛應(yīng)用于甘肅、陜西、內(nèi)蒙古等新能源服務(wù)地區(qū)，已成為一種不可替代的數(shù)值預(yù)報(bào)產(chǎn)品，對(duì)BJ-RUC產(chǎn)品的質(zhì)量評(píng)估與本地化應(yīng)用是未來新能源領(lǐng)域開展數(shù)值預(yù)報(bào)服務(wù)的研究重點(diǎn)。

致謝：感謝中國(guó)華電集團(tuán)福新能源股份有限公司甘肅分公司提供SCADA數(shù)據(jù)以及審稿專家和編輯的寶貴建議

[1] Alois Schaffarczyk. Understanding Wind Power Technology: Theory, Deployment and Optimisation[M]. John Wiley & Sons, Ltd., 2014,1-44.

[2] Liu Z Y. Electric Power and Energy in China[M]. John Wiley & Sons Singapore Pte. Ltd., 2013,73-136.

[3] 孫川永,陶樹旺,羅勇,等. 高分辨率中尺度數(shù)值模式在風(fēng)電場(chǎng)風(fēng)速預(yù)報(bào)中的應(yīng)用[J]. 太陽(yáng)能學(xué)報(bào), 2009,30(8):1097-1099.

[4] Ishak A M, Srivastava P K, Gupta M, et al. The Development of Numerical Weather Models-A review[J]. Bulletin of Environmental and Scientific Research, 2014,3(1):15-20.

[5] 江瀅,宋麗莉,程興宏. 風(fēng)電場(chǎng)風(fēng)速預(yù)報(bào)集合訂正方法的嘗試性研究[J]. 資源科學(xué), 2013,35(3):673-680.

[6] 盛春巖,劉詩(shī)軍,肖明靜. 十一運(yùn)會(huì)WRF-RUC系統(tǒng)應(yīng)用及典型個(gè)例評(píng)估[J]. 氣象科技, 2010,38(增刊):6-12.

[7] Powers J G. Numerical prediction of an Antarctic severe wind event with the Weather Research and Forecasting (WRF) model[J]. Monthly Weather Review, 2007,135(9):3134-3157.

[8] Jiménez P A, González-Rouco J F, García-Bustamante E, et al. Surface wind regionalization over complex terrain:evaluation and analysis of a high-resolution WRF simulation[J]. Journal of Applied Meteorology and Climatology, 2010,49(2):268-287.

[9] Storm B, Basu S. The WRF model forecast-derived low-level wind shear climatology over the United States Great Plains[J]. Energies, 2010,3(2):258-276.

[10] 李軍,胡非,王丙蘭,等. 風(fēng)速的Weibull分布參數(shù)[J]. 太陽(yáng)能學(xué)報(bào), 2012,33(10):1667-1671.

[11] Wang C H, Jin S L. The Features of the Errors and their Possible Causes from a Low-level Wind Field Simulated by WRF[J]. Wind Energy. 2013, DOI:10.1002/we.1635.

[12] 魏東,尤鳳春,范水勇,等. 北京快速更新循環(huán)預(yù)報(bào)系統(tǒng)(BJ-RUC)模式探空質(zhì)量評(píng)估分析[J]. 氣象,2010,36(8):72-80.

[13] Benjamin S G, Dévényi D, Weygandt S S, et al. An hourly assimilation-forecast cycle:The RUC[J]. Monthly Weather Review, 2004,132(2):495-518.

[14] 劉永前,史潔,楊勇平,等. 基于預(yù)測(cè)誤差分布特性的風(fēng)電場(chǎng)短期功率預(yù)測(cè)不確定性研究[J]. 太陽(yáng)能學(xué)報(bào), 2012,33(12):2179-2184.

[15] Thomas Ackermann. Wind power in power systems[M]. John Wiley & Sons, Ltd., 2005,25-51.

[16] Bowden G, Barker P R, Shestopal V, et al. The Weibull distribution function and wind power statistics[J]. Wind Engineering, 1983,7(2):85-98.

[17] Mertens S. The energy yield of roof mounted wind turbines[J]. Wind Engineering, 2003,27(6):507-518.

Quality Assessment of Low-level Wind Prediction Product by BJ-RUC Model：A Case Study of Akesai Wind Farm in Gansu Province

LIU Lijun1, LIANG Youjia2

(1.GansuMeteorologicalServiceCentre,Lanzhou730020,China; 2.SchoolofResourcesandEnvironmentalEngineering,WuhanUniversityofTechnology,Wuhan430070,China)

Based on the observation data of Dangjinshan wind farm with 49.5 mW power in Akesai Kazak county of Gansu Province in 2013, the low-level wind prediction product of Beijing Rapid Update Cycle (BJ-RUC) with high-resolution system was assessed. The results are as follows: (1) The forecast effect for the effective wind speed (3.5 - 11.5 m·s-1) was better than the rated wind speed (11.5 - 25 m·s-1) in four seasons as a whole, and that was the best in autumn and the worst in winter. (2) The Weibull distribution of wind speed frequency from the observation in a year was close to the normal distribution, while that from the simulation deviated from the normal distribution. The Weibull distribution of wind speed frequency in other months was similar to annual wind speed frequency except January. Thekvalue from the simulated by BJ-RUC model in each month approximated 2.0, which indicated that the wind was neutral and simulated gustiness with strong wind speed by BJ-RUC model had limitations. (3) By comparison of the wind rose diagram, the simulations of dominant wind direction for four different levels wind in 2013 were greatly differ from the observation on the whole, but that for the approximative rated wind speed was accurate, which indicated the simulation of wind direction by BJ-RUC model need to be improved.

wind speed; wind direction; BJ-RUC model; error analysis

10.11755/j.issn.1006-7639(2016)-04-0743

2015-10-12；改回日期：2016-01-20

國(guó)家自然科學(xué)基金(41601184)和中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(2016IVA044)共同資助

劉麗珺(1988- )，女，甘肅定西人，工程師，主要從事CFD風(fēng)場(chǎng)模擬研究. E-mail: chenshi2124@126.com

梁友嘉(1985- )，男，甘肅慶陽(yáng)人，講師，主要從事生態(tài)系統(tǒng)評(píng)價(jià)和環(huán)境建模研究. Email:yjliang@whut.edu.cn

1006-7639(2016)-04-0743-09DOI:10.11755/j.issn.1006-7639(2016)-04-0743

P457.5

A

劉麗珺，梁友嘉.BJ-RUC模式的低層風(fēng)場(chǎng)預(yù)報(bào)產(chǎn)品質(zhì)量檢驗(yàn)——以甘肅省當(dāng)金山風(fēng)電場(chǎng)為例[J].干旱氣象，2016，34(4)：743-751, [LIU Lijun, LIANG Youjia. Quality Assessment of Low-level Wind Prediction Product by BJ-RUC Model：A Case Study of Akesai Wind Farm in Gansu Province[J]. Journal of Arid Meteorology, 2016, 34(4):743-751],