看到這個(gè)標(biāo)題，你大概會(huì)聯(lián)想到氣象局的天氣預(yù)報(bào)吧！但事實(shí)上，這里說(shuō)的風(fēng)功率預(yù)報(bào)是指對(duì)未來(lái)24小時(shí)乃至未來(lái)3天的風(fēng)電場(chǎng)輸出功率做出預(yù)測(cè)的短期風(fēng)功率預(yù)報(bào)，其預(yù)測(cè)的準(zhǔn)確性不僅可以為電力部門制定調(diào)整發(fā)電計(jì)劃提供參考，還可降低風(fēng)電輸出功率大幅波動(dòng)給輸電網(wǎng)帶來(lái)的風(fēng)險(xiǎn)。
　　
　　目前，市場(chǎng)上主流的短期風(fēng)功率預(yù)報(bào)產(chǎn)品都是依賴于數(shù)值天氣預(yù)報(bào)結(jié)果開(kāi)發(fā)的。這是因?yàn)榻乇盹L(fēng)機(jī)輪轂高度風(fēng)的強(qiáng)度及風(fēng)向，受天氣過(guò)程、低空急流、局地地形強(qiáng)迫、晴空湍流等不同時(shí)間、空間、尺度的物理過(guò)程共同作用，形成了變化復(fù)雜且缺乏顯著的周期性特征。因此，當(dāng)風(fēng)速及功率預(yù)報(bào)時(shí)長(zhǎng)超過(guò)數(shù)小時(shí)之后，再優(yōu)秀的純粹統(tǒng)計(jì)算法也束手無(wú)策，這就需要數(shù)值天氣預(yù)報(bào)發(fā)揮威力了。也就是說(shuō)，需要先獲得風(fēng)電場(chǎng)所在數(shù)值天氣預(yù)報(bào)模式格點(diǎn)處的風(fēng)速、風(fēng)向及其他要素的預(yù)報(bào)結(jié)果，然后傳入到后續(xù)運(yùn)算模塊中，最終輸出功率預(yù)報(bào)。
　　
　　說(shuō)到這里，可能你還不了解數(shù)值天氣預(yù)報(bào)的原理。那么，本文就將對(duì)此做個(gè)簡(jiǎn)要介紹，以增進(jìn)你對(duì)數(shù)值天氣預(yù)報(bào)及其結(jié)果的理解。
　　
　　數(shù)值天氣預(yù)報(bào)的前世今生
　　
　　從史實(shí)來(lái)看，古人對(duì)天氣預(yù)報(bào)的認(rèn)識(shí)僅僅停留在看云識(shí)天氣以及出海前去龍王廟里拜一拜的階段。十九世紀(jì)和二十世紀(jì)之交，Abbe和Bjerknes提出可以基于大氣運(yùn)動(dòng)物理定律預(yù)報(bào)天氣，將天氣預(yù)報(bào)過(guò)程抽象為求解一組偏微分方程，而氣象要素觀測(cè)值可以作為偏微分方程的初值。到了二十年代，英國(guó)科學(xué)家Richardson嘗試著預(yù)報(bào)未來(lái)6小時(shí)歐洲地區(qū)地面氣壓場(chǎng)的變化，這也是有記載的人類第一次發(fā)布數(shù)值天氣預(yù)報(bào)。遺憾的是，這一次的預(yù)報(bào)結(jié)果與真實(shí)數(shù)據(jù)偏差太大，此后關(guān)于數(shù)值天氣預(yù)報(bào)的研究一度陷于沉寂。
　　
　　到了第二次世界大戰(zhàn)，電子計(jì)算機(jī)在此期間誕生，其強(qiáng)大的計(jì)算能力使得偏微分方程的求解擁有了時(shí)效性。氣象觀測(cè)網(wǎng)的逐步建立、特別是高空探測(cè)的發(fā)展，為數(shù)值天氣預(yù)報(bào)提供了越來(lái)越豐富的初值信息。與此同時(shí)，科學(xué)家對(duì)于數(shù)值天氣預(yù)報(bào)的理論研究也在不斷深入。他們認(rèn)識(shí)到Richardson之所以失敗，是因?yàn)樗蠼獾姆匠探M中包括了極易發(fā)散的高頻物理信息。1950年，Charney等美國(guó)氣象學(xué)家使用計(jì)算機(jī)求解，做了濾波近似后的方程組，首次成功發(fā)布了北美地區(qū)500Pa高度場(chǎng)的數(shù)值天氣預(yù)報(bào)結(jié)果。
　　
　　伴隨著氣象及科技的飛速發(fā)展，數(shù)值天氣預(yù)報(bào)能力在過(guò)去的40年里取得了長(zhǎng)足進(jìn)步。高性能超級(jí)計(jì)算機(jī)的運(yùn)算速度以指數(shù)級(jí)增長(zhǎng)，為數(shù)值模式時(shí)空分辨率的不斷提高提供了可能；觀測(cè)技術(shù)的發(fā)展，特別是氣象衛(wèi)星的發(fā)射，使得數(shù)據(jù)采集時(shí)不會(huì)遺漏全球任何一個(gè)角落。氣象模型以及針對(duì)模式無(wú)法直接求解的物理過(guò)程參數(shù)化不斷優(yōu)化，提高了模式描述真實(shí)物理世界的能力。下圖展示了數(shù)值天氣預(yù)報(bào)在南北半球預(yù)報(bào)能力的進(jìn)步（以500hPa位勢(shì)高度為比較基準(zhǔn)）。
　　
　　令人欣慰的是，如今對(duì)5天后氣壓場(chǎng)的預(yù)報(bào)精度已經(jīng)可以達(dá)到20年前對(duì)3天后氣壓場(chǎng)的預(yù)報(bào)，差不多每10年預(yù)報(bào)能力便提升1天！這是緩慢的進(jìn)步，卻也是巨大的進(jìn)步。 　　
　　主流天氣預(yù)報(bào)產(chǎn)品
　　
　　要強(qiáng)調(diào)的是，全球數(shù)值天氣預(yù)報(bào)的研發(fā)需要有強(qiáng)大的科研能力，且要投入大量的資源，世界上能夠玩得起這場(chǎng)燒錢競(jìng)賽的國(guó)家與組織屈指可數(shù)。
　　
　　要提的第一家便是美國(guó)大氣與海洋局，這是全球大氣科學(xué)學(xué)界的圣地，其開(kāi)發(fā)的全球預(yù)報(bào)系統(tǒng)(Global Forecast System, GFS)大概是這個(gè)世界上使用最為廣泛的數(shù)值天氣預(yù)報(bào)模式了。一方面GFS模式的預(yù)報(bào)質(zhì)量的確很值得信賴，另一方面也得益于其免費(fèi)政策。換句話說(shuō)，任何組織與個(gè)人不受國(guó)籍限制，都能自由獲取GFS模式數(shù)據(jù)并在此基礎(chǔ)上做研究。GFS每日發(fā)布4次，每次發(fā)布當(dāng)前時(shí)刻至未來(lái)15天的數(shù)值天氣預(yù)報(bào)，預(yù)報(bào)時(shí)長(zhǎng)分辨率為3小時(shí)。GFS預(yù)報(bào)數(shù)據(jù)的水平分辨率為0.25°，在中國(guó)區(qū)域相當(dāng)于近30km。
　　
　　另一家要提到的是歐洲中期天氣預(yù)報(bào)中心(European Center for Medium-Range Weather Forecasts, ECMWF)。歐洲中心的數(shù)值天氣預(yù)報(bào)質(zhì)量堪比GFS數(shù)據(jù)，但其對(duì)成員國(guó)外的用戶收取的服務(wù)費(fèi)用較為高昂，故其用戶比GFS用戶少很多。值得一提的是，為了獲得更高質(zhì)量的天氣預(yù)報(bào)數(shù)據(jù)，為客戶提供更為準(zhǔn)確的風(fēng)功率預(yù)報(bào)產(chǎn)品，遠(yuǎn)景也采購(gòu)了這一套預(yù)報(bào)數(shù)據(jù)。根據(jù)實(shí)測(cè)，這套數(shù)據(jù)的風(fēng)速預(yù)報(bào)精度的確比GFS預(yù)報(bào)結(jié)果要好。歐洲中心數(shù)據(jù)每日發(fā)布兩次，每次發(fā)布未來(lái)10天的數(shù)值預(yù)報(bào)，預(yù)報(bào)時(shí)長(zhǎng)分辨率亦為3小時(shí)。歐洲中心預(yù)報(bào)數(shù)據(jù)的水平分辨率為0.1°，相當(dāng)于10km，遠(yuǎn)超出GFS數(shù)據(jù)。
　　
　　我國(guó)對(duì)于數(shù)值天氣預(yù)報(bào)也非常重視，投入了大量人力、物力用于開(kāi)發(fā)自己的模式。中國(guó)氣象局自主研發(fā)的GRAPES全球預(yù)報(bào)系統(tǒng)2.0版本正式開(kāi)始業(yè)務(wù)化運(yùn)行并下發(fā)產(chǎn)品。該產(chǎn)品同化了大量衛(wèi)星資料，且在質(zhì)量控制和偏差訂正技術(shù)方面優(yōu)勢(shì)明顯。有報(bào)道稱本產(chǎn)品在中國(guó)區(qū)域的短期降水預(yù)報(bào)能力已接近歐洲中心，因此對(duì)于其在近地面處風(fēng)的預(yù)報(bào)能力表現(xiàn)也值得期待。
　　
　　數(shù)值天氣預(yù)報(bào)中的物理過(guò)程描述
　　
　　值得注意的是，大氣中不同物理過(guò)程在時(shí)間尺度和空間尺度上跨越了數(shù)個(gè)量級(jí)。大氣長(zhǎng)波的波長(zhǎng)長(zhǎng)達(dá)幾千千米，三四個(gè)波就能夠在緯向環(huán)繞地球一圈。臺(tái)風(fēng)、溫帶氣旋等過(guò)程的空間尺度較大氣長(zhǎng)波小一個(gè)量級(jí)，而海陸風(fēng)、山谷風(fēng)、強(qiáng)對(duì)流天氣等過(guò)程的尺度又要小一個(gè)量級(jí)。晴空積云的空間尺度只有千米量級(jí)，存在時(shí)間不過(guò)數(shù)小時(shí)，而湍流的空間尺度最小僅有數(shù)十厘米，生命時(shí)間為秒級(jí)。下圖展示了大氣中不同物理過(guò)程的空間尺度及特征時(shí)間。 　　
　　數(shù)值天氣預(yù)報(bào)的水平分辨率為十千米量級(jí)，簡(jiǎn)單理解就是在模式中十千米范圍內(nèi)就一個(gè)格點(diǎn)，這一個(gè)格點(diǎn)要代表十千米范圍內(nèi)的總體狀態(tài)。對(duì)于空間尺度遠(yuǎn)大于十千米的物理過(guò)程，在模式中可以用物理方程描述。
　　
　　問(wèn)題是那些尺度小于十千米的物理過(guò)程又該如何處理呢？這就需要引入?yún)?shù)化概念了，對(duì)于那些無(wú)法直接用方程顯示求解的過(guò)程，可以用一些經(jīng)驗(yàn)關(guān)系表達(dá)出來(lái)。比如對(duì)于尺度非常小的湍流項(xiàng)，模式通過(guò)格點(diǎn)的平均風(fēng)速的梯度進(jìn)行計(jì)算，但由于表達(dá)的只是湍流對(duì)模式格點(diǎn)風(fēng)速影響的統(tǒng)計(jì)特征，因此在具體的不同時(shí)刻可能會(huì)帶來(lái)不同的誤差引入，最終使得模式模擬的結(jié)果因時(shí)而異。