蔡 樂，趙學麗，王亞娟

（北京市水文總站，100089，北京）

降水的形成機制和預報屬于氣象工作的研究內(nèi)容，水文工作研究降水量時空分布特征和降水資料的整理及應用。降水量合理性檢驗是判斷降水量數(shù)據(jù)異常的有效方法，也是發(fā)布實時水情信息和降水量資料整編的重要前提。傳統(tǒng)的降水量合理性檢驗完全由人工完成，由于數(shù)據(jù)量很大，檢驗方法規(guī)則不統(tǒng)一，插補方式因人而異，從而導致效率低，數(shù)據(jù)的可靠性和準確性受到影響。因此，結(jié)合計算機的使用，尋求一種高效的降水量數(shù)據(jù)處理平臺，對水情發(fā)報和資料整編都具有重要意義。

2007年引進的水文水環(huán)境基礎(chǔ)信息管理系統(tǒng)（Water Management Information System Kisters，簡稱 WISKI）成功應用于北京水文、地下水、水環(huán)境業(yè)務系統(tǒng)。WISKI由德國Kisters公司研發(fā)，經(jīng)過20多年的開發(fā)與應用，已成為一個比較完善的軟件產(chǎn)品。該產(chǎn)品具有強大的業(yè)務擴展功能，可根據(jù)具體工作需要進行二次開發(fā)，進而配置出符合實際要求的各類統(tǒng)計、分析和糾錯功能。WISKI平臺的成功搭建為雨量合理性檢驗和插補工作提供了充分條件，同時也為尋求一種全新的水文數(shù)據(jù)整理模式提供了探索空間。

一、基本思路

WISKI系統(tǒng)在時序數(shù)據(jù)的基本運算和數(shù)據(jù)生成方面有獨特之處，具有較強的整體功能和完整性。WISKI本身沒有提供降水量合理性檢測的解決方案，但它提供了時序計算平臺和一系列通用的時序計算方法，可根據(jù)需求進行系統(tǒng)的二次開發(fā)，使得降水量合理性檢驗的適應性和靈活性更強。

根據(jù)合理性檢驗的需求，以獨立的雨量站降水量時序為基本單位，選取相關(guān)性最好的雨量站作為關(guān)聯(lián)站，進行每場次降水量數(shù)據(jù)的橫向比較，通過數(shù)據(jù)值之間的定量關(guān)系判別異常數(shù)據(jù)類型，從而完成降水量合理性檢驗，完整實施流程見圖1。

1.數(shù)據(jù)采集與錯誤數(shù)據(jù)類型

圖1 WISKI進行降水量合理性分析流程圖

表1 40a長時間序列雨量站之間相關(guān)系數(shù)所占總體百分比分布表

目前，北京市水文總站在全市范圍內(nèi)建設了121個遙測雨量站，雨量傳感器的分辨率為0.5mm，并且在雨水情系統(tǒng)建設中首次使用了CDMA—1X網(wǎng)絡作為主信道、GSM為備用信道的傳輸方案，針對目前的硬件設施情況和近年的使用經(jīng)驗，遙測降水量數(shù)據(jù)的錯誤主要有以下幾種類型：

①遷站引起數(shù)據(jù)缺失。由于山區(qū)房屋改造，遙測雨量站被迫遷站而引起數(shù)據(jù)缺失。

②儀器堵塞引起數(shù)據(jù)異常。遙測儀器出現(xiàn)堵塞，會造成一次降水后遙測數(shù)據(jù)連續(xù)多日維持在一個較小值范圍內(nèi)，或者無遙測數(shù)據(jù)。

③設備故障引起數(shù)據(jù)異常。由于設備故障，在較長時間內(nèi)沒有遙測數(shù)據(jù)。

④人為注水調(diào)試引起數(shù)據(jù)異常。人為注水調(diào)試的數(shù)據(jù)一般會人工剔除，但也不免遺漏未剔除的人工注水數(shù)據(jù)。

⑤RTU硬件原因?qū)е聰?shù)據(jù)異常。由RTU的硬件原因?qū)е碌漠惓?，主要表現(xiàn)為“跳數(shù)”，即數(shù)據(jù)出現(xiàn)較大的波動。

降水量合理性檢驗的核心工作就是對有以上問題而導致的異常數(shù)據(jù)值進行排查。

2.關(guān)聯(lián)雨量站的選取

（1）相關(guān)系數(shù)的計算

傳統(tǒng)人工進行雨量合理性檢驗是將同流域或者鄰近雨量站之間的數(shù)據(jù)進行對比，這種比較多數(shù)情況下是可行的，但是缺乏有力的科學根據(jù)和足夠的量化數(shù)據(jù)作支撐。雨量站之間距離的大?。ǖ匦螚l件）是降水是否同步的必要條件，而非充分條件，即關(guān)聯(lián)性強的雨量站距離一定很近，但距離很近的雨量站關(guān)聯(lián)性不一定強。我們需要引入一個具體的量化指標來判斷雨量站之間關(guān)聯(lián)性的強弱，因此，相關(guān)系數(shù)的計算是必不可少的。

選取長時間序列（1979—2009年）的日雨量值進行兩站之間相關(guān)系數(shù)計算，同時間的降水量值分別為序列 x1，x2，x3…xn-1，xn和 y1，y2，y3…yn-1，yn，n的數(shù)值不取決于兩個時間序列中的任何一個時序數(shù)，某時間序列的降水量值不為0，與其對應另一組時序便參與相關(guān)系數(shù)計算，無雨量取0值。采用EXCEL求相關(guān)系數(shù)方便快捷，CORREL函數(shù)自動將同為空白的信息濾掉，滿足分析的要求。雨量站之間的相關(guān)系數(shù)分布情況見表1。通過計算可以得出，平原區(qū)雨量站之間的相關(guān)系數(shù)普遍較大，山區(qū)雨量站之間相關(guān)性不及平原區(qū)明顯。通常相關(guān)系數(shù)大于0.8時，兩組變量有很強的線性相關(guān)性，為了考慮山區(qū)雨量站相關(guān)性略低的影響，選取相關(guān)系數(shù)大于0.75的雨量站為候選關(guān)聯(lián)雨量站。

（2）雙累計分析

通過長時間序列求得的相關(guān)系數(shù)可有效反映日降水量數(shù)據(jù)相關(guān)性，關(guān)聯(lián)雨量站的選取還需要在相關(guān)性較強的雨量站之間選擇數(shù)據(jù)一致性好的。WISKI自帶的雙累計分析為數(shù)據(jù)一致性檢查提供了有效方式。選取5分鐘時段的雨量值進行雙累計曲線分析，通過分析雙累計曲線可知，相關(guān)系數(shù)高的雨量站之間雙累積曲線形態(tài)一般較好，但也存在極個別特殊情況。浦洼站與十渡站40年降水量資料相關(guān)系數(shù)為0.92，但是雙累計曲線形態(tài)不如與黃塔站的好，浦洼站與黃塔站40年降水量資料相關(guān)系數(shù)為0.88。

通過相關(guān)系數(shù)的計算和雙累積曲線分析，每個雨量站選取相關(guān)系數(shù)高、雙累積曲線形態(tài)較好的3個雨量站為其關(guān)聯(lián)站，在WISKI中建立日雨量檢查時序，配置該時序的源信息，時間序列配置為待檢驗雨量站及其關(guān)聯(lián)站的日降水量值。

3.規(guī)則的制定

降水量合理性檢驗規(guī)則的制定是一個突破傳統(tǒng)，由定性分析向定量分析轉(zhuǎn)變的一個步驟?？蓞⒖棘F(xiàn)成的降水量合理性檢查原則極少，降水量合理性檢驗規(guī)則的制定參考以往的工作經(jīng)驗進行。針對降水量數(shù)據(jù)可能出現(xiàn)的錯誤情況，共制定了如下幾種異常值類型，并且統(tǒng)一編碼，用異常代碼的形式展示。將該雨量站的錯誤類型總體分為降水量大于0與等于0兩大類，再細化分具體的異常類型（見表2）。

表2 雨量站遙測數(shù)據(jù)異常代碼統(tǒng)計表

表3 5月降水量合理性檢查成果表（部分）（異常值成立:true，不成立:false）

4.成果輸出

WISKI將所有的測量值和計算值都以時序數(shù)據(jù)方式保存，每個時序數(shù)據(jù)對應于一個站點的指定參數(shù)，可以配置多種數(shù)據(jù)輸出模式，根據(jù)參數(shù)與時段生成相關(guān)的報表，并且可以設置在指定的時間自動生成報表，為成果自動輸出提供了一個合理的解決方案。由于降水量合理性檢查沒有具體的定量指標來判別異常值一定成立，所以最終輸出的成果還需要進行人工校核，但如果在檢驗規(guī)則較為完善的前提下，人工校核的工作量就會很小。首先經(jīng)WISKI平臺自動將疑似異常值的數(shù)據(jù)從海量的降水量數(shù)據(jù)庫中挑出，再根據(jù)人工的經(jīng)驗對WISKI挑選的異常值進行排查，可有效減少數(shù)據(jù)的遺漏，進而達到提高效率、減少錯誤的目的。

二、應用實例

表3是運用WISKI生成的5月降水量合理性檢查成果表的部分信息，經(jīng)人工校核，WISKI自動生成的成果表1/2以上的異常值均成立，并且無一異常值被遺漏。

三、結(jié) 語

①二次開發(fā)WISKI使其應用于雨量合理性檢驗是可行的。這種全新的數(shù)據(jù)處理方式將以往大量繁雜的人工工作通過計算機平臺實現(xiàn)，減輕了工作強度，提高了時效性，同時也提高了水文資料處理的先進性、科學性，具有積極的意義。

②由于WISKI平臺二次開發(fā)具有很強的靈活性，雨量合理性檢驗規(guī)則可根據(jù)實際需求具體制定，并可以根據(jù)工作經(jīng)驗不斷完善，使其更接近于真實的降水情況。當然，使用計算機平臺進行水文資料的處理代替了大量人工工作量，但其中人的作用還是不可或缺的，還需要經(jīng)驗豐富的水文工作人員對從海量數(shù)據(jù)中提取的異常值進行排查。這樣人工經(jīng)驗和計算機運用得到了有機結(jié)合，提高了水文數(shù)據(jù)處理的可靠性。

③WISKI應用于雨量合理性檢驗僅是二次開發(fā)WISKI平臺的探索之一。根據(jù)水文工作的實際需求對WISKI進行有目的的二次開發(fā)，可實現(xiàn)很多水文數(shù)據(jù)的分析計算功能。除WISKI之外的一些平臺也可經(jīng)過開發(fā)代替部分人工勞作，可應用到水文數(shù)據(jù)的處理工作中，促進水文行業(yè)信息化技術(shù)的發(fā)展。

[1]降水量測驗規(guī)范（SL 21—2006）[M].北京：中國水利水電出版社，2006.

[2]水文資料整編規(guī)范（SL 247—1999）[M].北京：中國水利水電出版社,2000.

[3]趙學麗，池宸星，吳海山.WISKI在北京水文數(shù)據(jù)管理中的應用[J].北京水務，2008（6）.

[4]吳金塔.水文遙測數(shù)據(jù)智能糾錯及插補技術(shù)在洪水預報中的應用研究[J].水利水電科技，2005（11）.