400 likes | 582 Views
Progress report. 2009/03/19. 研究流程. 第二部份. 第一部份. 雨量輸入. 雷達估計降水 --- KDP-R ZH-R(300) ZDR-R DSD-R Zh-R(Steiner) 雨量內插 --- 距離平方反比權重 徐昇式多邊形法. 雨量站. 模擬區域 -- 石門水庫霞雲集水區. 集水區面積為 623 平方公里,土壤以壤土為主,山區多為闊葉林及混合林為主。 數值高程資料( DTM )為農林航測所製作,網格解析度為 40 公尺,本研究將資料內插到網格解析度為 200 公尺。. CASC-2D 敏感度分析. 第一部份.
E N D
Progress report 2009/03/19
研究流程 第二部份 第一部份
雨量輸入 雷達估計降水--- • KDP-R • ZH-R(300) • ZDR-R • DSD-R • Zh-R(Steiner) 雨量內插--- • 距離平方反比權重 • 徐昇式多邊形法 雨量站
模擬區域--石門水庫霞雲集水區 • 集水區面積為623平方公里,土壤以壤土為主,山區多為闊葉林及混合林為主。 • 數值高程資料(DTM)為農林航測所製作,網格解析度為40公尺,本研究將資料內插到網格解析度為200公尺。
第一部份 WRF模擬2005年龍王颱風降雨
WRF模式設定 • 初始場資料: NCEP-FNL資料(網格大小:1º*1 º) • 模式設定:
模擬結果 士杰學長的模擬結果 自己的模擬結果 摘至中央氣象局網站http://www.cwb.gov.tw/
第二部份 2001年納莉颱風五分山雷達資料推估降雨流量模擬結果
第二部份 2005年龍王颱風雷達推估降雨結果及CASC2D與FLO-2D的流量比較
RMSE EFFIC
雷達降雨估計 • Z-R relationship(張偉裕2002、 Javier et al.,2007) 層狀降水a=302.6、b=1.28 對流降水a=268.9、b=1.34 B=3.77 對流層狀判斷:Steiner法(Steiner et al.1995)
雷達降雨估計 • KDP-R relationship r1 r2
雨量內插方法 • 距離平方反比 雨量站4 雨量站1 D4 內插點 D1 D3 D2 雨量站3 雨量站2
雨量站內插 • 徐昇式 雨量站2 雨量站1 雨量站3
dry Hf溼潤邊緣壓力 wet
對流降水與層狀降水的判斷 • Steiner法(Steiner,1995) 對流降水中心 對流降水 判斷位置 網格降水資料點 位於對流中心影響半徑內 層狀降水 位於對流中心影響半徑外 Zg : 網格回波強度(dBZ) MBR:平均背景回波強度(Mean Background Reflectivity) (dBZ) 取自Steiner et al.,1995
雷達估計降雨公式 • 差異反射率(ZDR)與降雨關係式—R3 A=3*10-3 B=1.22 (Zrnic et al .,1999)
雷達估計降雨公式 • 雨滴粒徑譜與降雨關係式—R4(DSD-R) (In ms-1 for Deq in mm) (Doviak and Zrnic, 1993)
雷達估計降雨公式 • Λ:斜率 • μ:散佈形狀因子 • N0:常數 (Brandes et al. , 2002)
評估方式 • RMSE(root mean square error) Qis為第i筆模擬流量值(m3/hr)、 Qio為第i筆觀測流量值(m3/hr)、 N為觀測資料與模擬資料之個數。 RMSE 值越小代表模擬流量結果越好。
評估方式 • EFFIC(efficiency coefficient)(Kirkby et al. 1987) Qio為第i筆資料的觀測流量(m3/s) Qis為第i筆資料的模擬流量(m3/s) 為觀測值的平均流量(m3/s),n為觀測資料與模擬資料之個數 • EFFIC等於1則表示模擬值與觀測值完全符合,若等於0表示模擬結果有達到觀測流量平均值的水準,所以效率係數值越接近1模擬結果越佳。