水文時間序列趨勢與突變分析系統開發與應用

王毓森

（甘肅省水文水資源局，甘肅 蘭州 730000）

水文時間序列趨勢與突變分析系統開發與應用

王毓森

（甘肅省水文水資源局，甘肅 蘭州 730000）

在借鑒趨勢檢驗和突變分析各種統計應用軟件使用方法的基礎上，結合水文工作實際，對數據輸入格式進行統一約定，對計算分析過程進行模塊化編程，對數據輸出結果進行規范化處理，基于Net Framework4.5、Visual Basic. Net和SharpDevelop 4環境開發完成了水文時間序列趨勢檢驗和突變分析系統。介紹了水文時間序列趨勢與突變分析系統的理論基礎、總體結構、設計功能及其應用。系統設計的趨勢分析和突變分析各項功能達到預期目標，較好地實現了對水文資料的一致性分析，分析結果科學、合理，對水文水資源工作具有一定的參考意義。

水文時間序列；趨勢分析；突變分析；系統開發

隨著全球氣候變化和人類活動的雙重影響，極端天氣和水突發事件日益頻繁，已嚴重危害到經濟社會發展和人類生存安全。為了有效避免和控制這些破壞性的全球環境變化事件，對水文氣象時間序列的趨勢變化、突變分析及預測研究是當前行業十分關注的問題之一。

目前，時間序列趨勢與突變分析各種方法在實際應用過程中散見于 MATLAB、DPS、SPSS、SAS、Excel等各種統計計算分析軟件中[1]。工作人員有時需要安裝、掌握多種應用軟件，有時甚至需要學習、編制程序來實現相關功能，致使原本簡單的工作復雜化、工作效率得不到提高，也使相關軟件在行業內部、單位內部得不到廣泛應用。為此，在借鑒趨勢檢驗和突變分析各種統計應用軟件使用方法的基礎上，對數據輸入格式進行統一約定，對計算分析過程進行模塊化編程，對數據輸出結果進行規范化處理，開發適合水文工作實際的趨勢檢驗和突變分析系統很有必要。

1 系統的理論基礎

水文時間序列是一個重要的自然要素時態數據，從大量水文基礎觀測中獲得，是一個觀測值序列。一般來說，一個時間序列是由一系列局部極點（最大點或最小點）組成的，這些局部極點通常被稱為跳躍點、臨界點、轉折點，因為這些極點顯現出了這個時間序列的趨勢變化及突變情況。在日常水文計算分析中，水文時間序列資料必須滿足可靠性、一致性、代表性要求，而趨勢檢驗和突變分析則是一致性分析的重要內容[2]。一致性分析是指獲得的水文時間序列是否來自同一總體，序列中的每個數據是否具有同一成因。通過一致性分析，即可判斷水文時間序列變化趨勢及是否發生突變，為進一步還原資料、分析氣候變化和人類活動對其影響提供依據。

水文氣象要素趨勢檢驗和突變分析方法大多基于現代數學和統計學理論，普遍采用參數統計、非參數檢驗等分析方法，并在其廣泛研究的基礎上，凝練出水文氣象序列趨勢與突變分析的理論與方法體系，這為水文時間序列趨勢與突變分析奠定了良好的理論基礎[3]。謝平等人提出的“水文變異診斷系統”，為今后理論方法的進一步改進及應用發展提供了參考[4]。

2 系統設計

2.1 系統框架

經過對水文時間序列資料的整理分析，根據分析計算數據流向，可將系統總體結構概況為：數據準備、數據預處理、數據分析、結果保存4個部分，如圖1所示。

圖1 系統開發思路及總體結構圖

“數據準備”即對初始觀測資料進行收集、整理，了解數據質量、數量，若采用經過整匯編的水文資料則更為可靠。這部分數據格式常見的有DAT、XLS、TXT等。

“數據預處理”即對收集到的水文資料按系統規定的文本格式（TXT）進行處理，主要信息包括資料名稱、資料類別、起訖時間、數據序列等。

“數據分析”即選擇不同檢驗方法，進行趨勢和跳躍分析，得出結論。

“結果保存”即將原始數據、中間計算過程、檢驗結果按系統規定格式（TXT、JPG等）進行保存。

2.2 系統功能

根據前述分析可知系統主要功能是趨勢分析、突變分析及結果保存，每個功能又根據具體的應用方法分為若干子功能，如圖2所示。

圖2 系統功能分布圖

“趨勢分析”功能實現了定性分析、定量檢驗，主要有累計濾波器法、肯德爾秩次相關法、斯波曼秩次相關法、線性趨勢相關法及非參數檢驗 （Z檢驗）等方法。

“突變分析”功能實現了累計距平、有序聚類、M-K等檢驗方法。在推廣應用中，又根據用戶的建議陸續實現了滑動T、滑動F、滑動游程、滑動秩和和MWP檢驗等方法。

“結果保存”功能實現了水文時間序列原始數據、每個子功能中間計算關鍵結果、常規檢驗結果的保存及子功能計算結果直觀展示。特別是對中間計算關鍵結果和常規檢驗結果的處理極大的方便了用戶在Excel等數據處理軟件中進行二次編輯，如圖3所示，以便引用。

圖3 經Excel辦公軟件二次編輯的分析成果

2.3 開發技術

水文時間序列趨勢與突變分析系統基于.Net Framework4.5、Visual Basic.Net和Sharp Develop4環境開發。NET Framework是微軟用于Windows托管代碼編程模型的組件，它以通用語言運行庫（CLR）為基礎，實現多種語言 （C#、VB、C++、F、Python、ASP、HTML等）的開發。Visual Basic.Net是基于微軟.NET Framework的面向對象編程語言，是目前最流行的Windows平臺應用程序開發語言之一[5]。Sharp Develop4是一款完全兼容visual Studio平臺的所有語言的輕量級開源IDE，支持插件擴展，具有安裝體積小、編譯封裝快捷、跨平臺運行等特點[6]。

3 系統應用實例

系統編制完成后，在各類水文水資源專項研究工作中，特別是涉及到基礎研究數據“一致性”處理方面得到了廣泛應用，并在工作實踐中不斷加以完善。在《氣候變化對祁連山區水文循環的影響研究》中通過對57個降水、氣溫、蒸發代表站、21個徑流代表站資料分析，得到了祁連山區三大流域氣候變化趨勢及發生突變的事實[7]。在《黑河流域東部子水系近60年來泥沙演變規律分析》中通過對3個代表水文站泥沙、徑流資料的分析，得到了黑河干流泥沙變化的趨勢及發生突變的年份[8]。在《基于變異診斷分析的大通河流量預報模型研究》中通過對享堂水文站年平均流量資料分析，結果表明大通河流域自20世紀90年代以來發生了顯著突變，在此基礎上進行了序列還原，并建立了預報模型[9]。上述計算成果結合物理成因分析，對進一步分析氣候變化和人類活動對水文時間序列的影響提供了依據。

4 結論與建議

1）系統設計的功能達到了預期目標，較好的實現了對水文資料的一致性分析，分析結果科學、合理，對水文水資源工作具有一定的參考意義。

2）系統可借鑒“水文變異診斷系統”綜合診斷的理念，對各種計算方法進行綜合評價，給出權重系數，從而確定最終變異形式[10]。

3）系統應對比分析各種方法計算結果的差異性，對各突變點相異的可能原因進行分析，對突變點代表水文時間序列變化發生的具體階段有待進一步探索。

[1]唐啟義.DPS數據處理系統∶實驗設計、統計分析及數據挖掘[M].北京∶科學出版社,2010∶990-995.

[2]SL322-2013,建設項目水資源論證導則[S].北京∶中國水利水電出版社,2013∶12-13.

[3]周園園,師長興,范小黎,等.國內水文序列變異點分析方法與在各流域應用研究進展[J].地理科學進展,2011,30 (11)∶1361-1369.

[4]謝平,陳廣才,雷紅富,等.水文變異診斷系統[J].水力發電學報,2010,29(2)∶85-92.

[5]謝爾頓,等.NET4.5與Visual Basic 2012高級編程(王凈譯)[M].北京∶清華大學出版社,2014∶1-3.

[6]侯穆,等.C#軟件項目開發全程剖析：全面透視Sharp Develop軟件的開發內幕(薛興濤等譯)[M].北京∶清華大學出版社,2009∶2-5.

[7]牛最榮.氣候變化對祁連山區水文循環的影響研究[M].蘭州∶甘肅人民出版社,2013∶36-41.

[8]王毓森,黃維東,崔亮,等.黑河流域東部子水系近60年來泥沙演變規律分析[J].甘肅水利水電技術,2016,52(1)∶5-8.

[9]王毓森,黃維東.基于變異診斷分析的大通河流量預報模型研究[J].人民黃河,2016,38(2)∶19-23..

[10]謝平,陳廣才,李德,等.水文變異綜合診斷方法及其應用研究[J].水電能源科學,2005,23(4)∶11-14.

p333.6