甘肅省渭河流域天然徑流還原計算研究

孫棟元，王啟優(yōu)，牛最榮，牛瓊燕，賈 玲，趙文琦，周 敏

(1.甘肅農(nóng)業(yè)大學，甘肅 蘭州 730070；2.甘肅省水文站，甘肅 蘭州 730000)

正確評價水資源是合理開發(fā)利用、保護和統(tǒng)一管理水資源的基礎[1- 3]。精準的水資源評價和水資源數(shù)量計算是實現(xiàn)水資源可持續(xù)發(fā)展和利用的重要前提和關鍵內(nèi)容[4]。徑流還原計算為水資源評價、分配、調(diào)度和生產(chǎn)服務等方面提供關鍵支撐[4- 5]。隨著水利工程的不斷修建和人類社會加速發(fā)展，取用水量不斷增加是導致河川徑流量銳減的主要原因之一，同時下墊面條件受人類建設活動不斷影響，造成流域產(chǎn)匯流情況變化，從而影響徑流狀況，使得徑流觀測資料不能代表河流的天然狀況，導致資料代表性和一致性出現(xiàn)一定偏差，因此，需要將觀測徑流還原為天然狀況下的徑流量[3- 6]。通過對徑流的準確還原計算，摸清流域水資源狀況，使得還原結(jié)果滿足水資源評價、水利工程設計和規(guī)劃要求，為區(qū)域水資源合理配置、高效利用提供技術保障。

眾多學者針對不同區(qū)域、采用不同方法與模型開展了諸多徑流還原方面的計算與分析，涉及有分項調(diào)查、降水徑流關系等方法，各個方法各有其優(yōu)缺點。喬云峰等將投影尋蹤法應用于徑流還原計算研究[2]；陳佳蕾等采用SWAT模型進行徑流還原計算[3]；范輝等利用VIC模型研究了天然徑流還原[7]；劉強等利用BP神經(jīng)網(wǎng)絡模型進行了徑流量還原計算[8]；阿不力克木和宋星原對西北干旱地區(qū)徑流量還原計算有關問題進行了探討[4]；張洪波等探討了地下水強擾動地區(qū)的徑流還原計算方法[5]；夏傳清對雅魯藏布江中游徑流還原計算方法進行了探析[6]；薛青和張亮還原計算了渭干河徑流量[9]；吳英超等還原計算了撫河徑流[10]，牛最榮等基于徑流還現(xiàn)對洮河流域徑流變化特征進行了研究[11]。盡管徑流還原已有相關成果，但針對渭河流域等特定河流仍需進一步開展徑流還原計算研究。牛最榮、黃維東、王海榮、聶文晶和盧路等分別對渭河流域水資源、地表徑流、水沙演變、河流健康等方面進行了相關研究[12- 18]，而關于渭河流域徑流還原計算依然薄弱，需要進一步加強該方面研究。因此，本文以甘肅省渭河流域為研究區(qū)域，開展流域天然徑流還原計算研究，為流域科學調(diào)控和利用水資源提供技術支持。本文選取渭河流域武山、北道2個典型水文站，借助多年觀測資料，利用分項調(diào)查和降水徑流雙累積方法，對渭河天然徑流進一步進行還原研究，從而為流域水資源科學管理和優(yōu)化配置提供科學依據(jù)。

1 研究區(qū)概況

渭河位于甘肅省中東部，省內(nèi)渭河流域總面積5.91×104km2，其中渭河干流流域面積為2.56×104km2。

渭河源頭在甘肅渭源縣，干流由西向東，省內(nèi)干流長297km，南北兩岸有榜沙河、秦祁河等14條河流，水量南岸較豐，北岸較貧。榜沙河、大南河等5條支流多年平均徑流量超過1.0×108m3。

2 資料與方法

2.1 資料來源與監(jiān)測站點情況

數(shù)據(jù)資料來源于甘肅省水文部門。甘肅省渭河流域涉及到渭源、武山、甘谷、北道等10個水文站，本文主要選取武山和北道站，渭河水文站信息見表1。

表1 渭河水文站信息表

選用渭河及相鄰流域社棠、天水等79個雨量站資料進行面雨量計算。資料系列長度不夠或缺測的流量和降水量可以通過插補延長得到，并經(jīng)過合理性檢驗后確定相應數(shù)值。

2.2 研究方法

利用徑流還原完成對典型站徑流還原計算。借助分項還原原理和過程[11，19- 20]，完成徑流還原計算和合理性檢驗，其還原耗損量項目有農(nóng)業(yè)、工業(yè)和生活用水等，跨區(qū)域引入(出)水量，水庫蓄水變量等。采用降水與徑流量雙累積關系圖進行徑流修正，確定明確拐點年份，作為時間段劃分依據(jù)，同時修正相應時間段徑流量[11]。通過年徑流深值和修正系數(shù)等指標，計算修正后天然徑流結(jié)果。

3 徑流還原計算與分析

3.1 徑流還原計算

采用統(tǒng)計分析和分項還原方法，得出1956—2016年渭河流域武山、北道2站多年平均實測徑流量、多年平均還原水量和還原后天然徑流量分別為5.357×108、11.238×108m3，0.404×108、1.892×108m3和5.761×108、13.130×108m3。2站對比結(jié)果如圖1所示，還原后徑流量也呈現(xiàn)增加趨勢，武山和北道站均呈現(xiàn)先增加后減小的徑流量變化趨勢，說明人類活動影響下徑流量呈現(xiàn)減小的變化趨勢。

圖1 武山與北道水文站天然與實測多年徑流量變化

3.2 徑流一致性分析

通過對武山、北道站控制面積內(nèi)的29個和79個雨量站分別進行合理性和一致性分析，繪制武山、北道2站面平均年降水量與徑流量雙累積圖，發(fā)現(xiàn)兩者關系突變拐點在1993年，以此將1993年作為分割點，將1956—1993年和1994—2016年作為2個劃分時段，并完成一致性修正如圖2—3所示。

圖2 武山與北道水文站降水量-徑流量雙累積及修正曲線圖

2站的徑流修正系數(shù)根據(jù)降水量大小而確定，同時將降水與修正系數(shù)作為流域徑流還原修正參考依據(jù)見表2。

圖3 武山與北道站前后2個時段徑流深-降水量關系圖

表2 武山和北道站降水量與徑流修正系數(shù)關系表

依據(jù)修正后2站降水量，借助確定的相對應修正系數(shù)，計算得到修正后兩站徑流量分別為3.893×108、8.403×108m3，其結(jié)果與還原后結(jié)果相比分別減少1.868×108、4.727×108m3。武山、北道站天然徑流量與還原徑流量對比結(jié)果如圖4所示。

圖4 武山與北道站天然與修正后多年徑流量變化

4 結(jié)語

通過對渭河流域武山、北道徑流量修正還原計算，客觀反映流域水資源開發(fā)利用情況和新的變化規(guī)律。研究結(jié)果為流域水資源綜合評價、科學管理與有效調(diào)配以及水利工程管理規(guī)劃提供科學依據(jù)。徑流還原結(jié)果受多方面因素影響，強人類活動使得區(qū)域下墊面狀況更加趨于復雜，同時受人類活動引用水量的影響，其還原水量也受到一定影響，再加之還原方法的有限性，使得徑流還原呈現(xiàn)更多多變性和復雜性。因此，為了提高還原計算的準確可靠，應進一步加強變化環(huán)境下多方法、多因素徑流的還原計算研究，使得計算結(jié)果更精確地反映區(qū)域水資源動態(tài)變化規(guī)律。