RainfallQualityCheck.ppt

Rainfall Quality Check 李天浩 蔡孝忠 台灣大學土木研究所 研究目的 雨量觀測資料的正確與否，影響後續需要降雨資料的各項分析。由於測站數目龐大，由人工判別方式來檢核資料的正確性，無法有效率地即時剔除不合理的數據。故本計畫乃針對翡翠水庫之雨量站及相關聯之氣象局測站，進行歷史資料分析，由合理之檢核邏輯為基礎，建立自動化的檢核模式，提供作業人員做合理性的判斷依據。 研究方式概述 分析歷史降雨資料，找出與欲檢驗測站之降雨相關性高的測站，計算一個降雨量估計值及信心區間，用以判斷測站的觀測值是否合理。 颱風或暴雨期間常觀測到局部極大值，此種高離群值常是正確且不易由內插而得，故假設降雨觀測值是屬於對數常態分布。 分析方式 一、建立待測雨量站檢核之參考測站 1. 選擇和待雨量站距離小於30 km之雨量站為候 選參考測站 2. 選取同一時間待測站和參考測站記錄值皆大於 零的資料對 3. 計算待測站和所有候選參考測站資料對取對數 後的相關係數 4. 選擇相關係數大於0.4的候選參考測站為檢核之 參考測站 二、建立檢核係數資料庫 計算相關係數時是將資料數值為零的去除，但是實際檢核時，資料為零或大於零者皆須進行檢核，故將降雨量為零的資料以一個未知常數參數代表利用所有用來檢核的數據估計誤差變異係數最小的原則，決定未知常數參數，每個待測站的常數參數未必相同，若未來資料數據增加，常數參數也可能會改變。 a. 檢核邏輯 假設待測站A共有n個檢核參考測站，計算A站的觀測值 的估計值 (1) A站之估計值 A站大於零的歷史資料取對數後之平均值 檢核參考測站i同一時間資料取對數之數值 檢核參考測站i大於零的歷史資料取對數之數值 常數權重因子 檢核參考測站觀測值為零的代表常數參數值 b. 推估常數參數 常數參數 的求得方式是由待測站A的m筆估計值，和m筆觀測數值的差值平方和為最小的條件來估計，利用逼近法估計最佳的 數值。 三、評量檢核系數的適當性 研究發現某些測站的權重係數絕對數值過大時，估計值受到這幾個特定站觀測值的影響，會使估計值的可靠度降低，故結合克利金法得到另一組迴歸係數 做為檢核係數適當性之評量。當迴歸法檢核權重係數與克利金法權重係數偏離程度過大時，則由克利金法所求得的權重因子取而代之。 克利金法 克利金法起源於地質學家研究南非礦冶工程，主要以區域化變數理論探討自然資源在空間上分佈之相關性並應用於勘查及推估自然資源上。 是一種統計的內插方法，對於方程式假設為線性、最佳化、不偏差性的估計方法，考量變數的空間結構，以及提供克利金值標準偏差形式的內插誤差推估，可用推估模式校正前指定參數值的最小變異數。 克利金法之半變異圖 克利金法使用之半變異圖，由全省測站於1997年至2001年之時雨量觀測資料所求得。 檢核翡翠水庫所需之測站列表 時間序列測試 若連續三小時之降雨觀測值皆不通過檢驗式，且觀測值皆為零，則認為該測站故障，此後不再檢查該站，也不再利用該站檢核相鄰的測站，直到作業人員將故障的標誌移除;若檢驗時數未達三小時則假設該測站仍未確定已經故障，繼續利用雨量站空間協調性檢測方法，檢核該站之觀測值。 結論 本查核程式已於九十二年五月底安裝於中央氣象局，經測試對於測站觀測值為零，但周圍測站卻有明顯降雨的個案，可即時並且正確的自全省將近400個測站中挑出，檢驗結果正確性高。 對於局部最大降雨強度的雨量觀測，查核程式經常舉出可能發生誤差、需由作業人員檢核的測站，故再修正查核程式之檢驗標準，將誤差上限由原設計的2倍標準偏差放寬為3倍，以改進發生局部性大雨時，原檢核邏輯容易將觀測值列為可能發生錯誤測站的問題。 * * *