- Home
- 產品說明
預報系統 (DCA/CWB/GFS) 風能
本產品分為風速、風能密度、風能發電量。利用本局全球預報系統(CWB/GFS)所產製100米高度的氣象預報資料[東西風(m/s)、南北風(m/s)、氣壓(hpa)、溫度(°C)、比溼(kg/kg)],經由統計分析小組將模式產出的東西風及南北風透過誤差衰退平均法校正後轉換成風速,再透過簡易運算式將風速轉換成風能密度及風能發電量。
風能密度 : \(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公尺高度風場內插而得、氣壓場由σ層及地面氣壓內插而得。
什麼是DCA
誤差衰退平均法(Decaying Average:DCA)是簡化形式的卡門濾波(Kalman Filter),用意是讓最近時間的訊息擁有更多的權重,就好比我們會更重視一個人最近一份工作的表現狀態,而不是第一份工作到最近一份工作同等平均的表現。
因此,我們可以保留一段時間,數值模式預報與真實值之間的誤差,且透過DCA的方式令最近的誤差擁有較多的權重,時間越久遠的誤差影響越少,這麼作能夠保留些許過去的誤差訊息,也能讓最近的誤差表現出來,再利用此誤差來修正最新產出的數值預報模式。
誤差衰退平均法 (Decaying Average)
- 原理: 將數值模式近期模擬的平均誤差視為模式系統性偏差,藉由簡化形式卡門濾波(Kalman Filter)動態估計後予以移除。
- 作法:
- 誤差計算 以下式估計t時刻格點i,j的τ時間預報場誤差值\(b_i_^τ_{,j}(t)\)。其中,\(f_i_^τ_{,j}(t)\)為t時刻格點i,j的τ時間預報場預報值,而\(a_{i,j}(t)\)為t時刻格點i,j的分析值/真實值。
- 更新系統性偏差 以下式與t-1時刻校正結果更新系統性偏差\(B_i_^τ_{,j}(t)\)。其中權重係數w代表該更新值考量幾次近期計算的偏差資訊(圖1)。若以每日00Z模擬為例,權重係數0.02代表掌握最近100-150次(權重係數倒數的2-3倍)模式模擬與分析場的系統偏差。
- 以系統性偏差校正預報值 將t時刻格點i,j的τ時間預報場預報值減去更新後的系統性偏差值\(B_i_^τ_{,j}(t)\),\(F_i_^τ_{,j}(t)\)得到t時刻格點i,j的τ時間預報場校正值 ,即:
- 參考文獻: Cui, B., Z. Toth, Y. Zhu, and D. Hou, 2012: Bias correction for global ensemble forecast. Wea. Forecasting, 27, 396-410.
\(b_i_^τ_{,j}(t)=f_i_^τ_{,j}(t)-a_{i,j}(t)\)(1)
\(B_i_^τ_{,j}(t)=(1-w)⋅B_i_^τ_{,j}(t-1)+w⋅b_i_^τ_{,j}(t)\)(2)
\(F_i_^τ_{,j}(t)=f_i_^τ_{,j}(t)-B_i_^τ_{,j}(t)\)(3)
圖1逐次(橫軸)預報誤差佔系統性偏差比例(縱軸)隨不同衰退權重係數(w)的變化情形
(曲線下面積總和為一)
DCA的作業方式
DCA會於每日產製00Z、06Z、12Z、18Z逐六時預報資料(共192小時),以00Z為例,每日下午2點前會將DCA校正後的逐六時預報資料上傳至下游端。 DCA取模擬範圍逐網格點之東西向百米風(100m U)、南北向百米風(100m V)及地表短波輻射(Radiation),統計前一天的CWB GFS百米風數值預報模式與CWB WRF MDAS百米風分析場以及CWB GFS地表短波輻射數值預報模式與本局預報中心修正後的地面格點輻射分析場之間的誤差,再利用誤差來修正今天產出的數值預報。
風能計算方法
風能即為空氣流動所產生的動能,當空氣流動越快,其蘊含的動能越高。風力發電機的原理是利用風的動能轉動發電機來產生電力,而在目前發電機的有效利用風速範圍大約是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 理想功率曲線