安 東，吳寶國，宋 倍，李 琦

（黃河水利委員會(huì) 信息中心，河南鄭州450003）

1 背 景

水資源具有時(shí)空變化特性，隨著人類活動(dòng)和全球氣候變化的影響，極端水事件頻發(fā)，水資源短缺問題日益嚴(yán)重，水資源開發(fā)利用管理面臨嚴(yán)峻考驗(yàn)，必須實(shí)行最嚴(yán)格水資源管理制度[1]。與傳統(tǒng)的水資源管理相比，嚴(yán)格化、精細(xì)化和系統(tǒng)化是最嚴(yán)格水資源管理的三大基本特征，必須從時(shí)間和空間上進(jìn)行水資源精細(xì)化管理[2]，根據(jù)水資源情勢變化及時(shí)響應(yīng)、動(dòng)態(tài)調(diào)度，從自然、社會(huì)、市場多方面獲取信息，這就要求實(shí)時(shí)、動(dòng)態(tài)、精準(zhǔn)、全面、細(xì)致地掌握水資源情勢變化。

對水資源的時(shí)空變化趨勢進(jìn)行分析，是水資源情勢分析和預(yù)測的關(guān)鍵手段，在水資源趨勢分析學(xué)術(shù)研究中被廣泛使用。但是其在水資源管理部門實(shí)際業(yè)務(wù)工作中的應(yīng)用與普及遇到了很多難題，主要包括：①數(shù)據(jù)資料收集整理不便，趨勢分析需要大量、多專業(yè)的歷史數(shù)據(jù)和不斷更新的實(shí)時(shí)數(shù)據(jù)以保證分析的時(shí)效性和準(zhǔn)確性；②算法、工具的專業(yè)性強(qiáng)，且相對獨(dú)立和分散，應(yīng)用普及難度大；③缺乏全面、系統(tǒng)的趨勢分析業(yè)務(wù)梳理。

在塔里木河流域水資源信息化整合項(xiàng)目的建設(shè)過程中，為了解決全流域綜合水資源情勢分析的需求，基于信息化技術(shù)和水資源管理要求，建設(shè)了一套獨(dú)立的通用趨勢分析應(yīng)用系統(tǒng)，能夠方便地融合到水資源業(yè)務(wù)應(yīng)用平臺(tái)，易于推廣應(yīng)用。

2 需求分析

2.1 塔里木河流域水資源管理概況

塔里木河是我國最大的內(nèi)陸河，承擔(dān)著新疆5個(gè)地州的供水重任。塔里木河流域由九大水系組成，其中地表直接匯入塔里木河干流的源流包括阿克蘇河、葉爾羌河、和田河、開都-孔雀河[3]，稱為“四源一干”。目前塔里木河流域“四源一干”已基本建成采集、傳輸、存儲(chǔ)、分析、調(diào)度為一體的水聯(lián)網(wǎng)體系，以及服務(wù)于全流域水量統(tǒng)一調(diào)度的塔里木河流域管理局（簡稱塔管局）綜合業(yè)務(wù)應(yīng)用平臺(tái)。塔管局從管理體制和信息化手段上均已實(shí)現(xiàn)全流域統(tǒng)一調(diào)度和管理。

2.2 趨勢分析需求

水資源管理包括對水資源供、用、耗、排各個(gè)環(huán)節(jié)的統(tǒng)一管理，在水資源合理開發(fā)利用的同時(shí)，還要保障水資源可持續(xù)利用，降低對生態(tài)環(huán)境的影響。在進(jìn)行流域水資源統(tǒng)一管理與調(diào)度時(shí)，需要掌握水資源各環(huán)節(jié)的現(xiàn)狀特征及變化趨勢。從發(fā)生變化的水文水資源情勢中分析自然水循環(huán)水文特征的自然規(guī)律，準(zhǔn)確把握人工水文循環(huán)供、用、耗、排機(jī)理，揭示水資源供需變化和生態(tài)環(huán)境變化，是水資源合理配置和生態(tài)環(huán)境保護(hù)的依據(jù)，也是對實(shí)行最嚴(yán)格水資源管理、水生態(tài)文明建設(shè)效果的直觀反映[4-6]。

2.3 趨勢分析分類

針對塔管局的水資源管理需求，具體到水資源管理部門，可以將趨勢分析需求劃分為來水和供水、需水和取用水、用水效益、生態(tài)環(huán)境4類（見表1）。

表1 趨勢分析內(nèi)容分類

（1）來水和供水趨勢分析，主要對降水、徑流、水庫蓄水等因素的長序列和動(dòng)態(tài)實(shí)時(shí)變化趨勢進(jìn)行分析。

（2）需水和取用水趨勢分析，從需用水用途、水源兩方面進(jìn)行分析，從而嚴(yán)格用水總量控制，合理安排生產(chǎn)布局，與水資源條件相適應(yīng)。

（3）用水效益趨勢分析，包括萬元國內(nèi)生產(chǎn)總值用水量、萬元工業(yè)增加值用水量、農(nóng)田灌溉有效利用系數(shù)3個(gè)方面，是最嚴(yán)格水資源管理的重要考核指標(biāo)。

（4）生態(tài)環(huán)境趨勢分析，從地下水情、生態(tài)輸水、水域面積、植被覆蓋率等方面進(jìn)行綜合分析。

3 通用趨勢分析系統(tǒng)設(shè)計(jì)

3.1 總體設(shè)計(jì)

該系統(tǒng)依托水聯(lián)網(wǎng)體系提供強(qiáng)大的數(shù)據(jù)資源，利用成熟的水文序列趨勢分析方法，搭建數(shù)據(jù)庫層通用趨勢分析算法庫，提供多對象、多要素、多尺度的水資源綜合趨勢分析服務(wù)，建立趨勢分析專用圖庫實(shí)現(xiàn)水資源情勢分析的可視化表達(dá)，結(jié)合實(shí)際的水資源業(yè)務(wù)應(yīng)用系統(tǒng)，實(shí)現(xiàn)簡潔、實(shí)用、靈活的趨勢分析功能。通用趨勢分析系統(tǒng)的總體結(jié)構(gòu)包括業(yè)務(wù)應(yīng)用、分析服務(wù)、外部數(shù)據(jù)源和數(shù)據(jù)庫4部分（見圖1）。

圖1 通用趨勢分析系統(tǒng)總體結(jié)構(gòu)

（1）業(yè)務(wù)應(yīng)用。通用趨勢分析系統(tǒng)以Web方式，提供自定義趨勢分析、預(yù)定義趨勢分析、可視化趨勢分析組件3種標(biāo)準(zhǔn)化應(yīng)用。

（2）分析服務(wù)。通用趨勢分析系統(tǒng)提供趨勢分析處理服務(wù)，進(jìn)行定時(shí)或人工趨勢分析處理。該服務(wù)分別調(diào)用趨勢分析定義服務(wù)（定義趨勢分析的數(shù)據(jù)范圍、分析算法）、數(shù)據(jù)抓取服務(wù)（從外部數(shù)據(jù)源中獲取數(shù)據(jù)）和算法調(diào)用服務(wù)來實(shí)現(xiàn)功能處理。

（3）外部數(shù)據(jù)源。通用趨勢分析系統(tǒng)進(jìn)行趨勢分析的數(shù)據(jù)來源于所嵌入系統(tǒng)的水資源業(yè)務(wù)數(shù)據(jù)源。通過趨勢分析服務(wù)定義的數(shù)據(jù)從外部數(shù)據(jù)源中獲取。對于水調(diào)管理業(yè)務(wù)，主要包括來水水情、需水水情、取用水水情、下泄水情、水費(fèi)及社會(huì)效益、用水效益等水調(diào)業(yè)務(wù)數(shù)據(jù)庫。

（4）數(shù)據(jù)庫。通用趨勢分析系統(tǒng)數(shù)據(jù)庫包括算法庫和成果庫。其中：算法庫包括各類算法函數(shù)、算法參數(shù)基礎(chǔ)表等，成果庫存儲(chǔ)趨勢分析結(jié)果公式及系數(shù)、計(jì)算結(jié)果時(shí)序數(shù)組等。系統(tǒng)趨勢分析的原始數(shù)據(jù)來源于外部數(shù)據(jù)源。

3.2 業(yè)務(wù)應(yīng)用功能

針對水資源趨勢分析業(yè)務(wù)需求，充分考慮系統(tǒng)的易操作性、靈活性和可擴(kuò)展性，設(shè)計(jì)了3個(gè)主要功能模塊。

（1）自定義趨勢分析方案模塊。提供可自定義時(shí)間尺度、時(shí)間范圍和統(tǒng)計(jì)對象及要素的趨勢分析服務(wù)，以滿足水資源情勢分析涉及多對象、多要素、多尺度的分析要求。自定義設(shè)置的趨勢分析方案可作為模板保存，以方便后續(xù)使用。

（2）預(yù)定義趨勢分析方案模塊。該系統(tǒng)在分析水文水資源要素對流域水資源情勢影響的基礎(chǔ)上，預(yù)定義多種服務(wù)于水資源情勢分析的常用趨勢分析方案。方案中對趨勢分析的對象、要素、時(shí)間尺度、時(shí)間范圍均進(jìn)行了常規(guī)設(shè)定，選擇相應(yīng)方案即可查看來水趨勢分析、用水趨勢分析、下泄趨勢分析、生態(tài)趨勢分析、綜合趨勢分析等的結(jié)果。

（3）可視化趨勢分析模塊。將趨勢分析與水資源管理常用功能緊密結(jié)合，如在面向決策層的水資源綜合信息概覽功能中嵌入流域歷年來水變化趨勢圖，以及作為生態(tài)修復(fù)成效重要考核指標(biāo)的生態(tài)輸水量歷年變化趨勢圖和地下水位歷年變化趨勢圖等，為綜合決策提供直觀的水資源情勢分析。

3.3 處理流程

通用趨勢分析系統(tǒng)在處理流程設(shè)計(jì)上充分考慮用戶體驗(yàn)，提供操作簡單、結(jié)果明確的功能界面，其前后臺(tái)交互處理流程見圖2。

圖2 趨勢分析處理流程

（1）設(shè)置趨勢分析方案。提供簡單明確的功能界面，用戶根據(jù)需求選擇預(yù)定義的趨勢分析方案模板，或者自定義設(shè)置趨勢分析的對象要素類型、參與統(tǒng)計(jì)計(jì)算的時(shí)序數(shù)據(jù)的時(shí)間范圍以及時(shí)間粒度。

（2）獲取原始數(shù)據(jù)。根據(jù)用戶選擇的趨勢分析方案，從業(yè)務(wù)數(shù)據(jù)庫中獲取參與分析計(jì)算的原始時(shí)間序列數(shù)據(jù)。

（3）數(shù)據(jù)預(yù)處理。對原始時(shí)序數(shù)據(jù)進(jìn)行數(shù)據(jù)一致性和完整性預(yù)處理，包括異常數(shù)據(jù)檢查、缺漏數(shù)據(jù)處理、數(shù)據(jù)精度一致性處理等，降低因數(shù)據(jù)質(zhì)量參差不齊而導(dǎo)致的數(shù)據(jù)分析結(jié)果不可靠的概率。

（4）運(yùn)行算法進(jìn)行計(jì)算。將預(yù)處理后的中間數(shù)據(jù)輸入到選定的算法計(jì)算程序，運(yùn)行算法進(jìn)行趨勢分析計(jì)算。

（5）生成計(jì)算結(jié)果。趨勢分析算法的直接計(jì)算結(jié)果是趨勢線公式的系數(shù)。為實(shí)現(xiàn)計(jì)算結(jié)果的可視化展示，趨勢分析計(jì)算程序?qū)⑤敵鲆唤M時(shí)刻與數(shù)值對應(yīng)的時(shí)序數(shù)。

（6）可視化展示分析結(jié)果。在前端功能界面中嵌入趨勢圖庫中封裝的趨勢圖組件，將趨勢分析計(jì)算結(jié)果用專業(yè)的趨勢分析圖直觀地展示給用戶。

3.4 典型算法實(shí)現(xiàn)

趨勢分析的方法有很多，大體分為參數(shù)統(tǒng)計(jì)法和非參數(shù)統(tǒng)計(jì)法兩類。參數(shù)統(tǒng)計(jì)法是對原始數(shù)據(jù)直接統(tǒng)計(jì)分析，研究其整體變化趨勢，如線性滑動(dòng)平均法、累計(jì)距平法、線性回歸法等。非參數(shù)統(tǒng)計(jì)法是對基于原始數(shù)據(jù)獲得的秩進(jìn)行統(tǒng)計(jì)分析，其優(yōu)點(diǎn)是避免受少數(shù)異常值和缺失數(shù)據(jù)的影響，而且比參數(shù)統(tǒng)計(jì)法更適用于非正態(tài)分布，具體方法有Spearman秩相關(guān)系數(shù)法、Mann-Kendall趨勢分析法、小波分析法等[7-8]。根據(jù)塔管局綜合業(yè)務(wù)應(yīng)用平臺(tái)的數(shù)據(jù)情況和趨勢分析需求，在系統(tǒng)中采用了線性回歸最小二乘法。

最小二乘法回歸方程式為

式中：i為參與趨勢分析計(jì)算的原始數(shù)據(jù)個(gè)數(shù)；（Xi，Yi）為一組原始時(shí)序數(shù)據(jù)；、為方程的參數(shù)；ei為樣本（ Xi，Yi）的誤差。

在算法運(yùn)行中將原始數(shù)據(jù)的時(shí)間序列代入Xi、相應(yīng)時(shí)間對應(yīng)的原始數(shù)據(jù)代入Yi，即可求出趨勢線的標(biāo)準(zhǔn)方程。最后將時(shí)間序列代入趨勢線標(biāo)準(zhǔn)方程中，即可求得趨勢線上與時(shí)間序列對應(yīng)的趨勢線數(shù)值，該數(shù)據(jù)與時(shí)間序列對應(yīng)組成一組時(shí)間與數(shù)值對應(yīng)的時(shí)序數(shù)組，即繪制趨勢圖所需的結(jié)果數(shù)據(jù)。

4 塔管局綜合業(yè)務(wù)應(yīng)用平臺(tái)趨勢分析功能開發(fā)

4.1 系統(tǒng)總體結(jié)構(gòu)及趨勢分析功能

塔管局綜合業(yè)務(wù)應(yīng)用平臺(tái)（以下簡稱應(yīng)用平臺(tái)）主要由數(shù)據(jù)資源層、應(yīng)用支撐層、業(yè)務(wù)應(yīng)用層、數(shù)據(jù)交換平臺(tái)和應(yīng)用安全等部分組成（見圖3）。將趨勢分析系統(tǒng)的數(shù)據(jù)庫、服務(wù)層和應(yīng)用層分別部署到數(shù)據(jù)資源層、應(yīng)用支撐層、業(yè)務(wù)應(yīng)用層中。趨勢分析系統(tǒng)以內(nèi)嵌的方式與應(yīng)用平臺(tái)進(jìn)行整合，提供無縫的應(yīng)用集成。

圖3 塔管局綜合業(yè)務(wù)應(yīng)用平臺(tái)總體結(jié)構(gòu)

（1）數(shù)據(jù)資源層。趨勢分析系統(tǒng)的數(shù)據(jù)庫主要包括趨勢分析算法庫和分析結(jié)果庫，在應(yīng)用平臺(tái)的數(shù)據(jù)庫服務(wù)器中直接部署和初始化數(shù)據(jù)庫；趨勢分析系統(tǒng)所需要的水資源業(yè)務(wù)數(shù)據(jù)源直接使用應(yīng)用平臺(tái)的數(shù)據(jù)庫，通過數(shù)據(jù)訪問接口獲取數(shù)據(jù)。

（2）應(yīng)用支撐層。趨勢分析服務(wù)模塊部署在應(yīng)用支撐層的業(yè)務(wù)應(yīng)用服務(wù)中，使用與平臺(tái)兼容的服務(wù)協(xié)議。

（3）業(yè)務(wù)應(yīng)用層。趨勢分析應(yīng)用功能采用Web方式實(shí)現(xiàn)，根據(jù)應(yīng)用平臺(tái)的應(yīng)用需求，在應(yīng)用平臺(tái)中實(shí)現(xiàn)歷年來水下泄趨勢分析、歷年用水趨勢分析、生態(tài)斷面地下水情趨勢分析、歷次生態(tài)輸水趨勢分析、土地覆被變化趨勢分析、來水實(shí)時(shí)統(tǒng)計(jì)分析、用水實(shí)時(shí)統(tǒng)計(jì)分析等7個(gè)定制的數(shù)據(jù)統(tǒng)計(jì)和趨勢分析功能并嵌入到綜合監(jiān)控監(jiān)測和水資源管理模塊中。同時(shí)，采用Web組件方式的趨勢分析組件，內(nèi)嵌到綜合展示與服務(wù)信息系統(tǒng)的相關(guān)界面中。

4.2 趨勢分析功能設(shè)計(jì)

在應(yīng)用平臺(tái)中，基于趨勢分析系統(tǒng)的預(yù)定義趨勢分析模塊，以獨(dú)立功能的方式提供了以下7種趨勢分析功能。

（1）歷年來水下泄趨勢分析。對任意歷史年份各源流的上游來水情況及流域出口斷面的下泄情況進(jìn)行趨勢計(jì)算，呈現(xiàn)來水下泄的變化趨勢。

（2）歷年用水趨勢分析。對任意歷史年份地州、兵團(tuán)、市（縣）等各級用水單位的用水情況進(jìn)行趨勢分析。

（3）生態(tài)斷面地下水情趨勢分析。對任意時(shí)段生態(tài)斷面的地下水相關(guān)監(jiān)測指標(biāo)（如地下水位、埋深、礦化度、電導(dǎo)率等）進(jìn)行變化趨勢分析。

（4）歷次生態(tài)輸水趨勢分析。對任意歷史年份生態(tài)輸水情況進(jìn)行趨勢分析。

（5）土地覆被變化趨勢分析。對多年遙感影像解譯的土地覆被結(jié)果進(jìn)行趨勢分析計(jì)算，分析土地覆被變化情況。

（6）來水實(shí)時(shí)統(tǒng)計(jì)分析。根據(jù)每日實(shí)時(shí)更新的來水?dāng)?shù)據(jù)，進(jìn)行短期來水趨勢分析，及時(shí)掌握來水情勢，以便動(dòng)態(tài)響應(yīng)。

（7）用水實(shí)時(shí)統(tǒng)計(jì)分析。對實(shí)時(shí)采集的取水口取水?dāng)?shù)據(jù)、每日更新的用水單位用水量數(shù)據(jù)進(jìn)行趨勢分析，及時(shí)掌握用水情勢，以便開展實(shí)時(shí)調(diào)度。

4.3 典型功能介紹

4.3.1 來水下泄趨勢分析

流域來水變化趨勢是流域水資源情勢分析的關(guān)鍵因素，歷年來水變化趨勢分析是對流域當(dāng)前水資源量進(jìn)行評估以及對今后水資源量進(jìn)行預(yù)測的重要手段。來水下泄趨勢分析功能將流域來水變化趨勢和下泄變化趨勢分析結(jié)果呈現(xiàn)在同一個(gè)時(shí)間維度上，為流域水資源綜合情勢分析提供了更豐富的數(shù)據(jù)依據(jù)。根據(jù)用戶選擇的流域、時(shí)間范圍，從數(shù)據(jù)庫中獲取相應(yīng)的數(shù)據(jù)進(jìn)行趨勢分析計(jì)算，分析結(jié)果以統(tǒng)計(jì)表和折線圖的形式進(jìn)行展示（見圖4，圖中數(shù)據(jù)為測試數(shù)據(jù)）。

圖4 綜合趨勢分析結(jié)果展示

4.3.2 生態(tài)輸水歷年變化趨勢

塔里木河向下游生態(tài)輸水是流域生態(tài)文明建設(shè)的重要內(nèi)容之一，對恢復(fù)其尾閭臺(tái)特瑪湖生態(tài)環(huán)境至關(guān)重要。生態(tài)輸水量的逐年變化趨勢是生態(tài)修復(fù)工作成效的重要體現(xiàn)。生態(tài)輸水歷年變化趨勢分析功能根據(jù)用戶所選擇的年份時(shí)段查詢生態(tài)輸水?dāng)嗝娴闹鹉晗滦沽?，繼而對歷年生態(tài)輸水下泄量進(jìn)行趨勢分析計(jì)算，以統(tǒng)計(jì)表和折線圖的形式直觀地呈現(xiàn)趨勢分析結(jié)果（見圖5，圖中數(shù)據(jù)為測試數(shù)據(jù)）。

圖5 生態(tài)輸水歷年變化趨勢分析結(jié)果展示

5 結(jié) 語

基于本研究開發(fā)的趨勢分析系統(tǒng)在塔管局綜合業(yè)務(wù)應(yīng)用平臺(tái)中得到了實(shí)際驗(yàn)證，取得了良好的應(yīng)用效果。系統(tǒng)應(yīng)用將趨勢分析方法與水資源情勢分析業(yè)務(wù)相結(jié)合，增強(qiáng)了水資源管理的科學(xué)性，提升了對水資源管理過程及變化情況的監(jiān)管能力；實(shí)現(xiàn)了趨勢分析功能通用性和靈活性的有機(jī)統(tǒng)一，降低了應(yīng)用系統(tǒng)的設(shè)計(jì)開發(fā)難度；實(shí)現(xiàn)了與水聯(lián)網(wǎng)和大數(shù)據(jù)處理技術(shù)的結(jié)合，解決了趨勢分析應(yīng)用中數(shù)據(jù)獲取難、質(zhì)量差、總量小的問題，增強(qiáng)了水資源情勢趨勢分析的實(shí)時(shí)性和準(zhǔn)確性。