近60年來渭河甘肅段徑流演變及驅動因素

田晉華, 李正勤, 柯浩成, 梁 靚, 田孟涵, 高金芳

(1.甘肅省水利廳 蘭州水土保持科學試驗站, 蘭州 730020; 2.蘭州理工大學, 蘭州 730050)

水資源是影響流域環境的重要因素。近年來，全球性的氣候變化改變了流域水文循環態勢[1]。另一方面，人類活動諸如水利工程、水土保持工程等直接或間接地影響到流域產匯流過程，其中以黃土高原渭河流域等最為典型。一般認為，徑流過程受氣候變化及人類活動的綜合影響[2]，因此明確區分徑流變化原因有助于科學認識水文過程，為指導地區水資源利用與管理提供科學依據。

徑流變化驅動因素響應一直是行業研究的熱門問題，國內外學者關于氣候變化及人類活動對徑流過程的影響進行了大量研究。張建云等[3]通過模型場景預測，分析了徑流對氣候變化的響應，發現氣溫升高和降水減少導致黃河中游徑流減少。劉廷璽等[4]分析1957—2010年遼河流域徑流變化，定量研究徑流驅動因素響應，發現不同站點徑流變化對氣候變化與人類活動影響的響應程度差距較大。郭愛軍等[5]研究表明近50 a人類活動對渭河徑流減少的貢獻率接近70%。岳永杰等[6]研究發現降水是根河徑流變化的主要原因，但不同時期徑流對降水量的響應不同。龔珺夫等[7]研究延河流域徑流過程響應過程，發現人類活動對徑流減少的貢獻率達到56%，地區暖干化趨勢明顯。Li等[8]對漢江的研究表明氣候變化是各季徑流變化的主導因素，氣候變化與人類活動對徑流的影響呈現明顯的季節動態性。Mwangi等[9]對肯尼亞馬拉河流域的研究表明土地利用變化對徑流影響貢獻率達到97.5%。Yan等[10]通過VIC模型研究灤河流域，發現人類活動主要影響汛期徑流，耕地轉為建設用地對徑流變化有巨大影響。由此可見，氣候變化與人類活動對不同地區徑流的影響差距較大。

現有渭河徑流變化及驅動因素影響的研究一般都集中于渭河關中段，對于渭河甘肅段鮮有涉及。研究渭河甘肅段水文氣象要素的演變及驅動因素的響應，對甘肅省東南部水土保持與水資源高效利用具有十分重要的意義。

本文以渭河甘肅段為研究區域，分析地區徑流變化特征規律，在此基礎上計算氣候變化與人類活動對渭河徑流變化的貢獻率，為地區水土保持工作及水資源高效可持續利用提供理論依據。

1 材料與方法

1.1 研究區概況

渭河發源于甘肅省渭源縣，是黃河最大支流，其干支流經陜甘寧三省，于陜西潼關匯入黃河。渭河干流全長818 km，流域總面積13.48萬km2，其中甘肅省境內干流長316 km，北道站以上254 km，控制面積2.49萬km2。流域屬于溫帶季風氣候區，大陸性氣候特征顯著，降水年內分布差異大，主要集中在7—9月。渭河甘肅段區域年降水量400～600 mm，年平均氣溫4～10℃[11]。研究區位置見圖1。

圖1 研究區概況示意圖

1.2 數據來源

本文選取渭河干流北道水文站1953—2010年的月徑流資料數據，資料來自于甘肅省水資源公報。北道站是渭河甘肅段的控制站點，可以較好的反映渭河甘肅段的徑流特征。氣象數據資料選擇1953—2010年天水、華家嶺、崆峒3個站點月氣象數據資料，數據來自于中國氣象數據網(http:∥data.cma.cn/)。土地利用現狀遙感監測數據選用1980年與2010年30 m土地利用數據，數據來自于地理國情監測云平臺(http:∥www.dsac.cn/)。

1.3 研究方法

1.3.1 Mann-Kendall法 Mann-Kendall法[11]是一種非參數統計檢驗方法，其應用不受個別異常值的干擾，計算過程簡便，因此適用于水文氣象等非正態分布的數據，可用來分析水文氣象時間序列趨勢性與變異性。

構造時間序列x的秩序列Sk

(1)

式中:Sk是第i時刻數值大于第j時刻數值個數的累計數,i≥j。

定義統計量

(2)

式中：

(3)

(4)

其計算步驟為(1) 分別構造順序與逆序時間序列的秩序列Sk，計算對應UFk與UBk。(2) 給定置信度α，若UFk>Uα/2，證明通過置信度為α的顯著性檢驗，時間序列趨勢顯著。(3) 繪制UFk與UBk兩條曲線，臨界區間內的曲線交點為潛在突變發生時刻。

1.3.2 R/S分析法 Mann-Kendall檢驗只能對已有數據進行趨勢性判斷，如對未來趨勢進行判斷，則需要利用Hurst系數。R/S分析法[12]定義了定義極差R與標準差S。其比值即

R(t)/S(t)=(Ct)H

(5)

式中:H為Hurst指數,通過線性回歸分析即可求取H。

H>0.5，未來變化表現為與原有趨勢一致的正持續性，H越大，正持續性越強。H=0.5，則代表未來序列隨機獨立。如果0≤H<0.5，則代表反持續性。

1.3.3 累積距平法 累積距平法[13]可根據圖像曲線變動直觀判斷數據變化。對于序列x，其t時刻的累積距平值表示為

(6)

1.3.4 變差系數Cv變差系數Cv[14]是評價數據離散程度的相對指標，一般用來分析數據波動情況。Cv值越大，代表數據變化幅度越大。其計算公式為

(7)

式中:Ki為模比系數,即實測序列與均值的比值。

1.3.5 小波分析 小波分析[15]可以通過變換突出時間序列中的多種變化周期，因此用來分析水文氣象要素的周期性變化。在水文氣象要素分析中選擇小波函數時，一般選擇復Morlet小波函數做為基函數。為降低噪聲影響需要先對水文要素作標準化處理，其標準化方法如下[16]:

(8)

1.3.6 累積量斜率變化率法 累積量斜率變化率法[17]是一種定量分離出氣候變化和人類活動對徑流量變化的貢獻率的方法，該方法通過計算基準期與人類活動影響期各累積量斜率變化率與徑流累積量斜率變化率的比值，從而計算各因素變化對徑流變化影響的貢獻率。

將基準期定為a時期，影響期定為b時期，則前后兩個時期的徑流量為SRa與SRb，氣溫分別為STa與STb，降水分別為SPa與SPb。設降水對徑流變化的貢獻率為CP，氣溫對徑流變化貢獻率為CT，人類活動對徑流變化貢獻率Ch=1-CP-CT。其中

CP=(SPa-SPb)/|SPb|/((SRa-SRb)/|SRb|)×100%

(9)

CT=-(STa-STb)/|STb|/((SRa-SRb)|SRb|)×100%

(10)

2 結果與分析

2.1 渭河水文氣象要素變化

2.1.1 渭河甘肅段徑流變化特征 對渭河甘肅段水文氣象要素年際變化特征值進行計算處理，其詳細結果見表1，表2。渭河近58年年徑流量平均值為11.06億m3，最大年徑流量為30.34億m3(1967年)，最小年徑流量僅1.287億m3(1997年)。通過線性回歸分析對渭河甘肅段年徑流趨勢進行了判斷，其結果見圖2，由圖2可知，渭河甘肅段徑流量年際整體呈現下降趨勢，年徑流量以0.2億m3/a的速度減少。徑流序列Mann-Kendall檢驗通過置信水平α=0.01的顯著性檢驗，Hurst指數H=0.87，表明未來一段時間內徑流仍會呈現下降趨勢，且下降趨勢十分強烈。

表1 渭河甘肅段徑流降水氣溫年際變化特征值

表2 水文氣象要素趨勢分析相關參數

2.1.2 渭河甘肅段降水變化特征 渭河甘肅段近58年平均年降水量499.84 mm，其年際變化情況見圖3。由圖3可知，渭河甘肅段年降水量整體呈下降趨勢，平均以1.52 mm/a的趨勢減少。降水序列Mann-Kendall檢驗通過置信水平α=0.05的顯著性檢驗，Hurst指數H=0.72，證明降水量未來下降趨勢依然較為顯著。

圖2 徑流量年際趨勢

2.1.3 渭河甘肅段氣溫變化特征 近58 a來，渭河流域甘肅段年平均氣溫為7.88℃，其年際變化情況見圖4，由圖4可知流域平均氣溫呈上升趨勢，其H=0.76，表明趨勢性較強。年平均最低氣溫序列H=0.88，趨勢性十分強烈，年平均最高氣溫序列H=0.67，趨勢性較強，3個特征值Mann-Kendall檢驗均通過置信水平α=0.01的顯著性檢驗，其中最低氣溫上升對于氣溫升高的貢獻更大，未來氣溫整體依然會延續顯著上升趨勢，與全球變暖背景相對應。

圖3 降水量年際趨勢

圖4 年平均氣溫年際趨勢

2.1.4 渭河甘肅段水文氣象要素趨勢預測 通過R/S分析法對渭河甘肅段徑流、降水及氣溫序列進行了分析，發現計算其Hurst指數均大于0.5，其中徑流序列H=0.87，降水序列H=0.72，平均氣溫序列H=0.88，3個特征均表現為強正持續性，即未來的水文氣象要素總體變化趨勢與過去趨勢相一致。結合Mann-Kendall趨勢檢驗結果，認為未來一段時期內，渭河流域甘肅段整體會沿著氣溫升高、徑流減少的趨勢發展。

2.2 渭河甘肅段水文氣象要素變異特征

2.2.1 渭河甘肅段徑流變異特征 繪制年徑流累積距平曲線判斷變異點見圖5。由圖5可以看出，1986年與1991年前后，累積距平值出現了趨勢變化，由增加轉為減少。

通過Mann-Kendall法對徑流序列進行變異點的識別，取置信水平α=0.05，U(0.05/2)=1.96，繪制Mann-Kendall檢驗曲線圖，其結果見圖6。由圖6可知，UF與UB兩條曲線在1991年前后相交，交點處于兩條臨界線之內，因此可以認為徑流序列突變發生在1991年。

圖5 徑流量累積距平

圖6 徑流Mann-Kendall檢驗

2.2.2 渭河甘肅段降水變異特征 繪制降水累積距平曲線見圖7。由圖7可知，渭河甘肅段降水量累積距平值在1991年前后發生趨勢改變，認為1991年為降水序列突變年。

通過Mann-Kendall法對徑流序列進行變異點的識別，取置信水平α=0.05，U(0.05/2)=1.96，繪制降水序列Mann-Kendall檢驗曲線圖，其結果見圖8。由圖8可知，渭河甘肅段降水序列在置信區間內出現了多個交點，表明研究期內存在多個潛在變異點。

圖7 降水累積距平

圖8 降水Mann-Kendall檢驗

2.2.3 渭河甘肅段徑流變異年份 有研究表明，進入20世紀90年代后，渭河流域出現干旱，且持續時間較長[18]。綜合徑流與累積距平曲線、Mann-Kendall檢驗曲線及實證分析，認為降水與徑流均于1991年發生了突變。因此將1953—1990年認為是受人類活動影響較小的基準期，1991—2010年為人類活動影響期。

2.3 渭河甘肅段徑流周期性變化

通過Matlab的小波分析工具箱[19]，計算處理序列的小波系數實部和小波方差，制作小波系數實部值等值線圖與時間尺度小波方差圖，分別見圖9—10。小波系數實部值的正負代表年徑流量豐枯，正值為豐水年，負值為枯水年，小波方差出現峰值的時間尺度即為徑流變化的周期。小波方差的峰值越大，則代表該尺度的周期性越明顯[20]。

分析小波系數實部等值線圖，研究發現徑流序列存在18～23 a，8～18 a，3～4 a等3個時間尺度的周期特征，小波方差圖中存在3 a，14 a以及23 a等3個明顯的峰值，與小波系數實部等值線圖相對應。其中3 a尺度的變化在六七十年代表現穩定，23 a尺度的變化在1975年后表現穩定。14 a尺度對應小波方差最大峰值，說明14 a尺度左右的周期性更強，徑流序列周期變化主要受14 a尺度周期控制，整個58 a在14 a時間尺度上經歷了約3個豐—枯轉換周期。渭河徑流序列在近58 a內存在明顯的豐枯變化，上述3個周期控制著徑流序列在整個研究時間區間內的的變化特征，周期性特征顯著。

圖9 徑流小波系數實部等值線

圖10 小波方差

2.4 渭河甘肅段徑流變化影響因素定量分析

利用高橋蒸發模型[21]對地區蒸發量進行估算，發現渭河甘肅段年平均蒸發量為363 mm/a，整體蒸發量呈現微弱下降態勢，線性回歸方程系數b=-0.312 4，蒸發量以每年0.3 mm的速度減少，相較于蒸發總量，可認為蒸發量變化不明顯。降水量的減少直接導致徑流量的減少，由于蒸發量相對變化不大，而徑流資料變動比氣象資料變動更劇烈，由此定性推測人類活動影響在其中發揮較大作用。

研究采用累積量斜率變化率法定量計算各種影響因素對徑流變化的貢獻率，該方法計算方便，已在黃河流域與黑河流域有較好應用[22]。

渭河流域甘肅段人類活動影響期較基準期相比，徑流量下降8.02億m3，下降率高達57.99%。降水量影響期較基準期下降61.97 mm，下降率為11.89%。年平均氣溫上升0.9℃，上升率為11.89%(表3)。

不同時期各影響要素累積線性擬合方程見表4，對影響因素貢獻率進行定量分析，發現氣候變化貢獻率為37.57%(降水貢獻15.99%，氣溫貢獻21.58%)，人類活動影響貢獻率高達62.43%。可以看出，在渭河甘肅段，人類活動的對徑流減少的影響極為顯著。

表3 基準期與人類活動影響期多年平均值變化

表4 不同時期累積量線性擬合方程

3 討 論

3.1 氣候變化對渭河甘肅段徑流的影響

氣候變化是改變徑流的主要因素。一般來說，徑流量與降水量成正相關，與氣溫呈負相關[23]。降水減少直接導致徑流減少，氣溫變化與潛在蒸發量成正相關關系[24]，氣溫升高，區域內的潛在蒸發量增加。夏軍等[25]研究表明渭河流域自1956—2000年干旱程度逐漸加深，氣溫增高導致徑流系數減小，且認為氣候的暖干化將人類活動的不利影響進一步放大。研究表明，渭河甘肅段近58年來年降水量下降78.3 mm，年平均氣溫上升1.44℃，兩個特征變化趨勢均對徑流減少有正向作用。

3.2 人類活動對渭河甘肅段徑流的影響

人類活動對徑流減少的影響主要體現在兩個方面，直接方面如人類大規模取水、調水工程等，導致徑流量的直接減少，間接方面在于人類活動例如城市地面硬化、水土保持工程等，改變了下墊面土地利用情況，影響到流域產匯流進而改變徑流量[26]。

本文選取基準期(1980年)與影響期(2010年)的土地利用現狀遙感監測數據，制作土地利用變化轉移矩陣分析基準期與影響期的土地利用變化情況，以分析人類活動對于徑流變化的影響(表5)。

表5 渭河甘肅段1980-2010年土地利用變化轉移矩陣 km2

從土地利用轉移矩陣中可以看出，從1980—2010年，渭河甘肅段建設土地面積顯著增加，增加幅度為25.37%，建設土地增加對應人類聚居區的擴大、人口增長與用水需求上升，有研究發現城市化背景下土地利用類型的變化是導致徑流變化的主要因素[10]。耕地總面積雖然減少，但是由于近年來渭河沿線諸如武山縣等大力發展蔬菜種植產業，加大了農業取水量，也在一定程度上導致渭河徑流量的減少。區域內林草總面積增加，植被通過增加降雨入滲、攔蓄徑流影響地表徑流[27]，導致地表徑流減少。此外，多年來渭河流域修建了眾多淤地壩，新修淤地壩具備攔蓄功能，淤地壩淤滿后則具備減洪作用[28]，綜合來看，黃土高原地區水土保持措施對徑流量的減少有促進作用[29]。

4 結 論

(1) 渭河甘肅段水文氣象要素特征分析表明，地區年徑流呈現極顯著下降趨勢，年降水呈現不顯著下降趨勢，年均氣溫呈現顯著上升趨勢。年徑流序列存在明顯突變，突變年份為1991年。氣溫、降水和徑流變化均表現為正持續性，未來一段時間內渭河流域甘肅段整體會沿著氣溫增加、徑流減少的趨勢發展。

(2) 渭河甘肅段徑流周期性分析表明，年徑流變化呈現明顯的周期性，存在3 a，14 a，23 a共3個時間尺度的周期特征，其中14 a的的時間尺度為第一主周期，徑流在過去58 a經歷了約3個豐—枯轉換期。

(3) 渭河甘肅段徑流驅動因素分析表明，人類活動影響期較基準期，多年平均徑流量減少8.02億m3，下降57.99%。其中氣候變化對徑流減少的貢獻率為37.57%(降水貢獻15.99%，氣溫貢獻21.58%)，人類活動影響占比62.43%。