- Home
- 產品說明
- 預報系統 (MDDS/CWB/GFS) 風能
- 1.1 什麼是GFS?
- 1.2 GFS的作業方式
- 1.3 GFS 100M風速的產出
- 1.4 什麼是MDDS?
- 1.5 MDDS的作業方式
- 1.6 風能計算方法
預報系統 (MDDS/CWB/GFS) 風能
本產品分為風速、風能密度、風能發電量。利用本局全球預報系統(CWB/GFS)所產製不同高度的氣象場資料,經過動力降尺度系統運算後,輸出100米高度的預報氣象資料[東西風(m/s)、南北風(m/s)、氣壓(hpa)、溫度(℃)、比溼(kg/kg)]。之後利用模式輸出的東西風及南北風計算出100米高度的風向及風速,再透過簡易運算式將風速轉換成風能密度及風能發電量。
風能密度 : \(P={C_P}×{ρ/2}×{A×V^3}\)
(\(C_p\)為風機功率係數,目前極限約為40%,ρ為密度,A為風機葉片掃掠面積,V為有效風速,3m/s~25m/s)
什麼是GFS?
本局從 1983 年起進行數值天氣預報的發展,並且於 1988 年上線第一代的全球數值天氣預報系統(簡稱 CWB/GFS),數值天氣預報包含資料同化系統與預報模式系統兩部份。現行作業版本為TCo639L72模式:水平解析度提升至約0.14個經緯度(約15公里),垂直座標更改為S−P hybrid座標,解析度也增加至72層,模式最上層頂為0.1 hPa。
GFS模式的運作,在前置資料處理的資料同化利用GSI (Gridpoint Statistical Interpolation),並採用三維系集變分混合法(Hybrid 3D-EnVAR)分析。觀測資料同時採用傳統(風場之Vor, Div,及T, Q, Terrain Pressure)及非傳統(輻射等)觀測資料在模式\(σ\)面上進行分析;預報的動力模組為全球波譜模式(TCo639L72),其動力核心由Eulerian(T511L60)升級至Semi-Lagrangian(TCo639L72),高斯網格點分布為2576 x 1280,垂直為S−P混合座標,共72層。預報變數包含渦度、散度、虛位溫、比濕、地表氣壓及雲水含量,時間積分採跳蛙和半隱式計算格式;物理參數化方面,在邊界層上採用Eddy Diffusivity/Mass Flux(EDMF) Monin (Siebesma et al. 2007),並加入尺度感知、深/淺對流Scale- and Aerosol- aware Simplified Arakawa-Schubert (Han et al. 2017 & Arakawa and Wu 2013);其它動力核心及參數化法,包含土壤模式、深/淺對流、網格尺度降水、地形/對流重力波拖曳及輻射參數化法則如表1所示。
表 1中央氣象局全球預報模式(TCo639L72)設定
變數 | TCo639L72 |
---|---|
資料同化 | Hybrid GSI |
水平解析度 | 15 公里 |
垂直解析度 | S-P Hybrid座標;72層 |
模式層頂 | 0.1 hpa |
網格尺度降雨 | Zhao and Carr (1997) |
動力核心 | Spectral method Semi-Lagrangian ( NDSL ) + Semi-implicit 3 time level |
深/淺對流 | Scale- and Aerosol- aware Simplified Arakawa-Schubert (Han et al. 2017 & Arakawa and Wu 2013) |
土壤模式 | Noah Land Surface model |
輻射 | RRTMG (Mlawar and Clough 1997; Hou 2002) |
PBL邊界層 | Eddy Diffusivity/Mass Flux(EDMF) Monin (Han et al. 2016) |
地形重力波拖曳 | Palmer et al. (1986) |
對流重力波拖曳 | stationary convectively forced gravity wave drag (Chun and Baik 1998) |
GFS的作業方式
目前CWB/GFS TCo639L72的模式作業流程為每日00Z、06Z、12Z、18Z產出,一共輸出16天的預報資料(大部分資料10天以內逐6小時、10天以上逐24小時)。於預報系統中僅使用8天內的預報資料,為配合綠能計劃需求,將100米資訊、2米溫度及地表短波輻射量資訊由1天增加輸出至16天逐3小時,並且利用模式輸出之100米高風場、溫度場、比濕、雲液態水混和比、以及離地表2 m高度之氣溫、地表太陽短波輻射量等氣象參數來做風能及太陽能資訊的監測及預報。
GFS 100M風速的產出
為了提供綠能開發計畫所需之100米高氣象資訊,利用CWB/GFS全球模式產出之σ層及近地表層之氣象場資料,以靜力方程及線性內插方式,計算離地表100公尺高之氣象場資料,並隨CWB/GFS全球模式每六小時預報一同輸出。
此模組提供離地表100公尺高之溫度、氣壓、風場、比濕以及雲液態水混和比。其中比濕以及雲液態水混和比由CWB/GFS全球模式之σ層資料內/外插而得、溫度場由σ層及2公尺高度氣溫內插而得、風場由σ層及10公尺高度風場內插而得、氣壓場由σ層及地面氣壓內插而得。
什麼是MDDS?
MDDS(Mesoscale Dynamic Downscaling System)為本局與美國國家海洋與大氣總署(NOAA)下之GSD (Global System Division)合作發展的複雜地形區的三維連續變分降尺度系統。降尺度是指利用統計法或動力法由尺度較大的模式結果得到較小尺度的現象。統計法是指利用統計的關係(通常為經驗關係),慎選較大尺度模式結果的參數(或稱之為預報因子),得到較小尺度的氣象特徵(如von Storch et al., 1993; Goodess and Palutikof, 1998; Wilby et al., 2004)。而動力法主要是利用較小尺度的模式,以較大尺度模式的結果為其邊界條件,而得到較小尺度的氣象特徵(如Zagar and Rakovec, 1999; Giorgi et al., 2001; Mearns et al., 2004),另外亦可廣義的利用物理關係為基礎,由尺度較大的模式結果得到較小尺度的現象。
在地形上的氣象場受地形的影響程度視地表的特性及所要解析的天氣系統尺度而決定。當地形越複雜且天氣系統尺度越小,氣象場受地形的影響越嚴重。在如此複雜的地形影響下,降尺度方法是否考慮地形的效應就決定了其結果的可信度。在複雜地形區的降尺度方法最簡便的即是利用區域數值模式的動力降尺度法,但此方法會面臨所有數值模式的極限,諸如初始場產生的模式調整問題(spin up problem)、模式動力及物理參數化問題等等。而在統計降尺度法中,Steincker et al. (2006) 利用其所謂的”fingerprint method”,將地形高度對氣象場的理想分布特性做為其方法的權重函數,可得到複雜地形區內較好的氣象場分佈。但此方法僅適用於單一的純量氣象變數,對向量氣象變數(風場)並不適用,且不同氣象變數間並無任何的物理關係存在。
本方法利用追隨地勢座標下的純控制方程三維變分法,試圖解出在靜態下因複雜地形產生的氣象場變化,且因為此方法利用較完整的中小尺度氣象物理方程式為變分法的控制方程,其所得結果在不同氣象變數間有相當大程度滿足其所使用的控制方程。考慮僅以動量方程、連續方程、熱力方程及水汽保守方程為控制方程的純三維變分,其價值函數可寫為
$$J={∭(α_1 P_1^2+α_2 P_2^2+α_3 P_3^2+α_4 P_4^2+α_5 P_5^2+α_6 P_6^2+α_7 P_7^2 )dΩ}$$
其中\(P_n\)為各項控制方程,包含三個方向的動量方程(\(P_1\)到\(P_3\))、連續方程(\(P_4\))、熱力方程(\(P_5\))、水汽保守方程(\(P_6\))及所有氣象變數的平滑懲罰項(\(P_7\)),\(α_n\)為各控制方程項的權重。
MDDS的作業方式
MDDS目標在於提供下游單位即時且不耗費大量計算資源的情況下,取得更高解析度的氣象場資訊,同時包含台灣複雜地形條件下氣象場更細微變化過程。由於接收到初始場為氣壓層場資料,為滿足MDDS之應用,會將其轉換成地勢追隨座標。
目前因應綠能計畫需求提供3種降尺度產品(MDDS-CWRF、MDDS-CGFS與MDDS-NCEP),水平與垂直解析度分別為2公里、250公尺,一共32層場資料且為追隨地勢座標下產品,降尺度範圍如圖一所示。即時作業流程每日4次(00、06、12、18 UTC),MDDS-CWRF初始條件使用目前本局區域預報模式系統之3公里水平解析度分析與預報場(CWBWRF-M04_WD02,DMS34key格式)進行動力降尺度,每次提供84小時逐時預報。MDDS-CGFS模式的初始條件使用目前本局全球預報模式系統之25公里水平解析度分析與預報場(CWBGFS-GH0G,DMS34key格式)進行動力降尺度,每次提供192小時逐6時預報。MDDS-NCEP模式的初始條件使用美國國家海洋和大氣管理局之全球預報模式之0.5度水平解析度分析與預報場(NCEP-GFS,grib2格式)進行動力降尺度,每次提供192小時逐6時預報。MDDS待上游分析與預報資料接收齊全之後立即啟動,平均計算時間花費約20~30分鐘,模式網格設定資訊及模式平均輸出時間列於附表一。本產品輸出格式為Bindary格式,提供給氣象資訊綠能營運中心使用,透過發電量評估模組介接系統將氣象資訊轉換風能發電量、風能密度等,為滿足離岸風機平均高度之需求,MDDS模式產品另外輸出高度100公尺層場資料。歷史備份資料為逐時(或逐六時)分析與預報場之100公尺高度之經緯度(deg)、東西風(m/s)、南北風(m/s)、氣壓(Pa)、相當位溫(K)、水氣含量(g/kg),預計儲存3年,規劃放置於氣象局綠能營運中心雲端資料庫。
表一、MDDS-CWRF模式網格設定資訊與輸出時間(LST)
CWBWRF(M04WD02) | MDDS/CWRF | |
投影法 | 藍伯特投影 | 藍伯特投影 |
座標系統 | 氣壓座標 | 追隨地勢座標 |
網格解析度 | 3公里 | 2 公里 |
X方向格點數 | 1158 | 257 |
Y方向格點數 | 673 | 257 |
垂直層數 | 24層 | 32層 |
參考經度 | 120° | 120.65° |
參考緯度 | 10°、40° | 30°、60° |
中心經度 | 122.27 | 120.63097 |
中心緯度 | 22.87 | 23.782600 |
模式頂高度 | 1000 Pa | 15500公尺 |
00Z(08LST) | 06Z(14LST) | 12Z(20LST) | 18Z(02LST) | |
CWBWRF(M04WD02) | 14:10 | 20:10 | 02:10 | 08:10 |
MDDS/CWRF | 14:40 | 20:40 | 02:40 | 08:40 |
表二、MDDS-CGFS模式網格設定資訊與輸出時間(LST)
CWBGFS(GH0G) | MDDS/CGFS | |
投影法 | 等經緯 | 藍伯特投影 |
座標系統 | 氣壓座標 | 追隨地勢座標 |
網格解析度 | 24公里(0.22度) | 2 公里 |
X方向格點數 | 1536 | 257 |
Y方向格點數 | 768 | 257 |
垂直層數 | 16層 | 32層 |
參考經度 | X | 120.65° |
參考緯度 | X | 30°、60° |
中心經度 | X | 120.63097 |
中心緯度 | X | 23.782600 |
模式頂高度 | 10000 Pa | 15500公尺 |
00Z(08LST) | 06Z(14LST) | 12Z(20LST) | 18Z(02LST) | |
CWBGFS(GH0G) | 13:40 | 19:40 | 01:40 | 07:40 |
MDDS/CGFS | 14:10 | 20:10 | 02:10 | 08:10 |
表三、MDDS-NCEP模式網格設定資訊與輸出時間(LST)
NCEPGFS | MDDS/NCEP | |
投影法 | 等經緯 | 藍伯特投影 |
座標系統 | 氣壓座標 | 追隨地勢座標 |
網格解析度 | 55公里(0.5度) | 2 公里 |
X方向格點數 | 720 | 257 |
Y方向格點數 | 361 | 257 |
垂直層數 | 26層 | 32層 |
參考經度 | X | 120.65° |
參考緯度 | X | 30°、60° |
中心經度 | X | 120.63097 |
中心緯度 | X | 23.782600 |
模式頂高度 | 1000Pa | 15500公尺 |
00Z(08LST) | 06Z(14LST) | 12Z(20LST) | 18Z(02LST) | |
NCEPGFS | 12:40 | 18:40 | 00:40 | 06:40 |
MDDS/NCEP | 13:10 | 19:10 | 00:10 | 07:10 |
圖一 、 MDDS降尺度範圍
風能計算方法
風能即為空氣流動所產生的動能,當空氣流動越快,其蘊含的動能越高。風力發電機的原理是利用風的動能轉動發電機來產生電力,而在目前發電機的有效利用風速範圍大約是3~25 m/s,因此在本計劃的風能計算方法,亦加入有效風速的限制條件,來計算有效風能密度。而根據貝茲極限理論,在物理上風力發電機是無法百分之百擷取風能,因此從風中取出的能量有一個極限值大約59%。
簡易的風能計算可由下式表示:
$$P={C_P×ρ/2×A×V^3}$$
(\(C_p\)為風機功率係數;本計劃目前設定為40% ;ρ為空氣密度;A為風機葉片掃掠面積,本計劃目前設定為1 \(m^2\);V為有效風速,3m/s~25m/s)
藉由上述方法取得風能密度,而風能發電量則從風力發電機的啟動風速至切出風速之功率曲線(圖1-1)去做計算。
圖1-1 理想功率曲線