近60年黃河上游干流水沙變化及其關(guān)系

許文龍，趙廣舉,?，穆興民,，高 鵬,，孫文義,

(1.中國(guó)科學(xué)院水利部水土保持研究所 黃土高原土壤侵蝕與旱地農(nóng)業(yè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，712100，陜西楊凌； 2.中國(guó)科學(xué)院大學(xué)，100049，北京; 3.西北農(nóng)林科技大學(xué) 黃土高原土壤侵蝕與旱地農(nóng)業(yè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，712100，陜西楊凌)

黃河水少沙多、水沙異源、水沙關(guān)系不協(xié)調(diào)，黃河的徑流量主要來(lái)源于上游的天然降水、冰雪融水和地下水。黃河水資源的時(shí)空分布對(duì)流域工農(nóng)業(yè)生產(chǎn)與人民生活有著極其重要的影響，但20世紀(jì)90年代以來(lái)，黃河干流實(shí)測(cè)徑流量顯著減少，已經(jīng)影響到該區(qū)農(nóng)業(yè)生產(chǎn)、人民生活及經(jīng)濟(jì)的可持續(xù)發(fā)展[1]。黃河的泥沙主要來(lái)源于黃土高原及風(fēng)沙區(qū)，大量泥沙隨風(fēng)、水進(jìn)入黃河，河道湖庫(kù)淤積萎縮，河床不斷抬升，導(dǎo)致主河槽行洪輸沙能力大幅降低，上游寧蒙河段淤積嚴(yán)重，下游河道泥沙淤積，“地上懸河”等現(xiàn)象加劇，嚴(yán)重影響黃河的防凌防汛工作，威脅沿途群眾的生命財(cái)產(chǎn)安全。進(jìn)入21世紀(jì)以來(lái)，黃河輸沙量從16億t銳減至3億t，泥沙銳減引起廣泛的關(guān)注，很多學(xué)者對(duì)此展開大量的研究工作[2]。

已有研究[3-6]表明，大量的人類活動(dòng)是黃河徑流量和輸沙量銳減的主要原因。穆興民等[7]指出不同時(shí)期黃河水沙變化的主要因素有所差異，上世紀(jì)80年代，壩庫(kù)工程攔蓄占據(jù)主導(dǎo)地位，而2000年以后，國(guó)家實(shí)施大規(guī)模的“退耕還林(草)”政策，植被恢復(fù)成為這一時(shí)期黃河水沙減少的主要原因。然而，大量研究多集中在黃河中游地區(qū)，尤其是多沙粗沙區(qū)的水沙時(shí)空演變規(guī)律與驅(qū)動(dòng)因素的相關(guān)研究。受氣候變化、農(nóng)業(yè)灌溉、大規(guī)模的水土保持措施和水庫(kù)建設(shè)等影響，黃河上游徑流量和輸沙量及水沙關(guān)系亦發(fā)生了顯著變化。筆者通過對(duì)唐乃亥至頭道拐區(qū)間10個(gè)干流水文站的徑流量和輸沙量變化分析，探究黃河上游水沙變化及水沙關(guān)系規(guī)律以及沿程變化特征，為流域水資源合理利用和河道綜合治理提供依據(jù)。

1 研究區(qū)概況

黃河上游位于E 95.0°～111.1°、N 32.51°～41.81°，流域面積為38.6萬(wàn)km2，占黃河流域總面積的51.3%；干流河長(zhǎng)為3 472 km，占黃河干流總長(zhǎng)度的61.3%；上游河段落差3 496 m，平均比降為1‰，年均降水量為426.2 mm[8]。

筆者選取唐乃亥至頭道拐區(qū)間作為研究區(qū)域(圖1)。根據(jù)景觀特征差異，以下河沿為分界點(diǎn)將該區(qū)間分為2段。上段為唐乃亥至下河沿，為峽谷地段，此區(qū)間貴德至蘭州段是上游支流集中區(qū)段，洮河、湟水等重要支流在此區(qū)段匯入黃河，是黃河上游徑流的重要補(bǔ)給區(qū)；同時(shí)唐乃亥至青銅峽區(qū)間又是我國(guó)重要的水電開發(fā)建設(shè)基地，已建成包括龍羊峽、劉家峽和青銅峽在內(nèi)的近20座水電站。下段為下河沿至頭道拐區(qū)間，屬于沖積平原段，區(qū)間大部為荒漠和荒漠草原，無(wú)大的支流匯入，干流河床平緩，水流緩慢，兩岸有大片沖積平原，其中包括銀川平原與河套平原等。

圖1 黃河上游干流地理位置、水系及水文站點(diǎn)分布Fig.1 Geographical location, river network and hydrological stations in the main stream of the upper reaches of the Yellow River

2 數(shù)據(jù)與方法

根據(jù)黃河上游干流水文站點(diǎn)的分布，兼顧研究區(qū)域分布省份和流域景觀差異，筆者選取黃河干流的唐乃亥站、貴德站、循化站、小川站、下河沿站、青銅峽站、石嘴山站、巴彥高勒站、三湖河口站及頭道拐站等10個(gè)沿程水文站的徑流量和輸沙量數(shù)據(jù)，研究黃河上游1950—2012年時(shí)間范圍內(nèi)的水沙關(guān)系、水沙變化及沿程變化特征。筆者所采用的水沙數(shù)據(jù)來(lái)自黃河上游干支流水文泥沙資料數(shù)據(jù)集(1950—2012年)。

水文序列可看作為非周期成分(趨勢(shì)項(xiàng)、突變項(xiàng)、跳躍項(xiàng)等)、周期成分和隨機(jī)成分的累加[9]，故可通過統(tǒng)計(jì)特征分析、趨勢(shì)性分析、突變性分析和周期性分析研究黃河上游水沙變化特征。

筆者將選用均值、方差、變異系數(shù)(CV，coefficient of variation)等[10]表征水文時(shí)間序列數(shù)據(jù)的基本統(tǒng)計(jì)特征，采用非參數(shù)Mann-Kendall檢驗(yàn)法(下稱MK方法)和Spearman秩相關(guān)系數(shù)法(下稱Spearman方法)[11-12]來(lái)分析各站點(diǎn)水沙變化的趨勢(shì)特征，采用累計(jì)距平法[13]識(shí)別水文序列的突變點(diǎn)，采用小波分析方法[14]判斷黃河上游徑流泥沙近60年的周期特征。由于本文采用的方法均為水文序列常用的統(tǒng)計(jì)分析方法，受篇幅限制，在此不再贅述。

3 結(jié)果與分析

3.1 徑流量、輸沙量年際變化的統(tǒng)計(jì)特征

根據(jù)黃河上游干流水文站點(diǎn)1950年至2012年之間的年徑流量、年輸沙量數(shù)據(jù)進(jìn)行統(tǒng)計(jì)(表1)，統(tǒng)計(jì)特征結(jié)果如下：多年平均徑流量和輸沙量最小值出現(xiàn)在唐乃亥站, 分別為201.33億m3和0.12億t；最大值出現(xiàn)在下河沿站，分別為300.71億m3和0.24 億t。變異系數(shù)描述數(shù)據(jù)與均值的偏移程度，當(dāng)CV≤0.1時(shí)，數(shù)據(jù)具有弱變異型，表明水文時(shí)間序列相對(duì)穩(wěn)定；當(dāng)0.1≤CV≤1時(shí)，數(shù)據(jù)達(dá)到中度變異程度；當(dāng)CV≥1時(shí)，水文時(shí)間序列具有強(qiáng)變異性[15]。各站點(diǎn)年徑流量的CV值在0.2～0.4的范圍內(nèi)，屬于中等變異性；年輸沙量的CV值較徑流量大，且小川站多年平均輸沙量的CV>1，屬于強(qiáng)變異性，其余屬于中等變異性。輸沙量變化變異程度較高，其中作為水庫(kù)出口控制的貴德、小川、青銅峽3站變化尤為明顯，顯然是受到龍羊峽(1986年蓄水，1989年運(yùn)行)、劉家峽(1968年蓄水)[16]、青銅峽(1968年運(yùn)行)3座水庫(kù)(圖1所示)蓄水調(diào)沙以及沿途大量工農(nóng)業(yè)生產(chǎn)生活用水的影響。采用Pearson相關(guān)分析和一元線性回歸分析描述黃河上游干流各站點(diǎn)的水沙關(guān)系，均呈現(xiàn)明顯的正相關(guān)性，除小川站相關(guān)性不顯著外，其余各站的相關(guān)系數(shù)均>0.50且在0.01的置信水平上顯著。

表1 黃河上游干流水文站點(diǎn)多年平均徑流量、輸沙量統(tǒng)計(jì)特征Tab.1 Statistical characteristics of multi-year runoff and sediment load at hydrological stations in the main stream of the upper reaches of the Yellow River

注：** 表示相關(guān)系數(shù)在0.01的置信水平上顯著；CV為變異系數(shù)。Notes: ** indicates that the correlation coefficient is significant at the confidence level of 0.01. CV is coefficient of variation.

圖2 黃河上游干流水文站點(diǎn)不同年代平均年徑流量、年輸沙量變化Fig.2 Variations of average annual runoff and sediment load at hydrological stations in the main stream of the upper reaches of the Yellow River during different time period

以10年尺度劃分年代進(jìn)行斷代分析，黃河上游干流各站不同年代平均徑流量和輸沙量均表現(xiàn)為下降趨勢(shì)。如圖2所示，20世紀(jì)60年代各站的平均徑流量最大，90年代的平均徑流量最小，上游徑流量在90年代顯著地減少；小川站及以下站點(diǎn)輸沙量代際變化呈現(xiàn)明顯下降趨勢(shì)，小川站以上各站變化幅度不大，較為平穩(wěn)，其中青銅峽站50年代到60年代輸沙量下降幅度較大。小川站是劉家峽水庫(kù)出口控制站，劉家峽水庫(kù)于1968年蓄水運(yùn)行，故70年代小川站及以下站點(diǎn)的輸沙量下降幅度較大。青銅峽水庫(kù)是衛(wèi)寧灌區(qū)和青銅峽灌區(qū)的分界點(diǎn)，因青銅峽水庫(kù)淤積和寧蒙灌區(qū)引水引沙使得60年代青銅峽站輸沙量大幅減少。表明水庫(kù)對(duì)黃河上游泥沙的攔蓄起到重要的作用。

3.2 徑流量、輸沙量年際變化趨勢(shì)性

采用非參數(shù)MK方法和Spearman方法判斷水文時(shí)間序列數(shù)據(jù)的趨勢(shì)性，并通過2種方法的對(duì)比以期得到更加可信的結(jié)果。表2為黃河上游干流沿程各水文站點(diǎn)水沙趨勢(shì)檢驗(yàn)統(tǒng)計(jì)量表。沿程各站點(diǎn)年徑流量和年輸沙量在近60年間均呈減少趨勢(shì)，且2種方法所得結(jié)果一致。唐乃亥站的年徑流量和年輸沙量減少趨勢(shì)不顯著，這與已有的研究結(jié)果一致，河源至唐乃亥區(qū)間無(wú)大型水利工程，故唐乃亥站是黃河上游受人類活動(dòng)影響最小的站點(diǎn)；除唐乃亥站之外各站年輸沙量減少趨勢(shì)均達(dá)到0.01的置信水平，黃河上游干流年輸沙量整體呈現(xiàn)下降趨勢(shì)；下河沿站至頭道拐區(qū)間各站的年徑流量和年輸沙量呈現(xiàn)極其顯著的減少趨勢(shì)，該段乃黃河上游沖積平原，黃河灌區(qū)取用水和河道泥沙淤積是影響該段水沙的重要因素；貴德和循化站年徑流量減少趨勢(shì)分別到達(dá)0.10和0.05的置信水平，小川站年徑流量減少趨勢(shì)不顯著。

表2 黃河上游干流水文站點(diǎn)年徑流量與年輸沙量變化趨勢(shì)檢驗(yàn)Tab.2 Trend test of average annual runoff and sediment load at hydrological stations in the main stream of the upper reaches of the Yellow River

注：*** 表示在0.01的置信水平上顯著；** 表示在0.05的置信水平上顯著；*表示在0.1的置信水平上顯著。Notes: ***, ** and * indicate that the correlation coefficient is significant at the confidence level of 0.01, 0.05 and 0.1.

3.3 徑流量、輸沙量年際變化突變特征

圖3為上游各站點(diǎn)年徑流量和年輸沙量的累積距平曲線。可以看出，徑流和輸沙累積距平曲線在1968年、1985年、1986年和1989年前后有明顯的升降趨勢(shì)，表明時(shí)間序列在該時(shí)間點(diǎn)發(fā)生突變。唐乃亥站年徑流量和輸沙量發(fā)生突變的年份都是1989年，貴德站和循化站年徑流量和輸沙量突變時(shí)間點(diǎn)分別是1985年和1986年，小川站的突變時(shí)間則是1986年和1968年，下河沿站至三湖河口站發(fā)生突變的年份是1985年和1968年，頭道拐站發(fā)生突變的年份都是1985年。

1968年和1986年分別是劉家峽水庫(kù)和龍羊峽水庫(kù)開始蓄水的時(shí)間，并且1986年后龍劉水庫(kù)開始聯(lián)合調(diào)蓄黃河上游水沙[16-17]，1989年11月后龍羊峽水庫(kù)開始正常運(yùn)行。這與各站點(diǎn)徑流量與輸沙量的突變時(shí)間相契合，進(jìn)一步證明龍劉水庫(kù)對(duì)黃河上游水沙的調(diào)控。

圖3 黃河上游干流水文站點(diǎn)年徑流量和年輸沙量變化的突變特征Fig.3 Abrupt changes of annual runoff and annual sediment load at hydrological stations in the main stream of the upper reaches of the Yellow River

3.4 月徑流量、月輸沙量變化周期性

圖4 黃河上游干流水文站點(diǎn)月徑流量和月輸沙量小波變換Fig.4 Wavelet analysis for monthly runoff and sediment load at hydrological stations in the main stream of the upper reaches of the Yellow River

對(duì)上游各站點(diǎn)的月徑流量和月輸沙量數(shù)據(jù)進(jìn)行小波分析。小波分析圖可以顯示同一時(shí)間點(diǎn)上不同周期的振蕩強(qiáng)度和同一周期內(nèi)不同時(shí)間點(diǎn)的分布，反映水文時(shí)間序列的豐枯變化特征[18]。由于所選站點(diǎn)數(shù)量較多，本文在分析10個(gè)站點(diǎn)月徑流量和月輸沙量的小波分析結(jié)果的基礎(chǔ)上，篩選4個(gè)具有代表性的站點(diǎn)(唐乃亥站、小川站、青銅峽站和頭道拐站，如圖4)。全部10個(gè)站點(diǎn)的徑流量和輸沙量具有顯著的1年的振蕩周期，巴彥高勒站、三湖河口站和頭道拐站的徑流量以及除貴德站外的其他9個(gè)站的輸沙量具有0.5年的顯著周期，且達(dá)到了95%的置信度。水沙0.5～1年的周期變化主要受季節(jié)和年際降水變化的影響，三者之間的變化密切相關(guān)[19]。此外各站水沙的小波功率譜圖也顯示出3～5和7～9年的變化周期，這與副熱帶高壓脊線位置的準(zhǔn)3年周期和地極移動(dòng)振幅變化、ENSO的7年左右的周期特征相一致[18]，而副熱帶高壓和ENSO都是影響我國(guó)西北地區(qū)天氣系統(tǒng)的重要因素。唐乃亥站水沙的周期變化至今未出現(xiàn)間斷，其他9個(gè)站點(diǎn)月徑流量的顯著周期自90年代中期出現(xiàn)了間斷性變化。貴德月輸沙量的周期性自90年代初開始消失，循化月輸沙量的周期性于90年代初開始減弱并在2000年后消失，其余7個(gè)站點(diǎn)月輸沙量在70年代初開始出現(xiàn)間斷性變化并逐漸減弱，同樣在2000年后基本消散，表明上游水沙在這些階段受到龍劉水庫(kù)運(yùn)行調(diào)節(jié)和退耕還林(草)等水土保持措施的影響較大。

3.5 沿程水沙變化關(guān)聯(lián)性

圖5 黃河上游干流水文站點(diǎn)多年平均徑流量、輸沙量沿程變化Fig.5 Variations for average annual runoff and sediment load at hydrological stations in the main stream of the upper reaches of the Yellow River

黃河上游干流自上而下10個(gè)水文站點(diǎn)多年平均和年代際平均徑流量、輸沙量沿程分布情況如圖5所示。除20世紀(jì)50年代外，上游徑流量自上游至下游表現(xiàn)為先增大再減少的態(tài)勢(shì)，各年代最大值均出現(xiàn)在下河沿站，青銅峽站的平均徑流量明顯低于其前后的下河沿站和石嘴山站。輸沙量自上而下表現(xiàn)為逐漸增大的趨勢(shì)，且在下河沿站處明顯增加，下河沿之后各站變化幅度不大，較為平穩(wěn)。20世紀(jì)50年代，青銅峽站的徑流量與下河沿站相近且高于石嘴山站，當(dāng)時(shí)只有青銅峽灌區(qū)在大量引水灌溉，而衛(wèi)寧灌區(qū)并未投入使用。青銅峽站的輸沙量高于下河沿站和石嘴山站，而石嘴山站的輸沙量又明顯低于巴彥高勒站，輸沙量波動(dòng)的原因是因?yàn)樽鎱柡?、清水河和紅柳溝泥沙的匯入、青銅峽灌區(qū)秦渠、唐徠渠和漢渠引水引沙以及烏蘭布和沙漠風(fēng)沙入黃。60年代以后，青銅峽站的徑流量均低于下河沿站和石嘴山站，表明衛(wèi)寧灌區(qū)和青銅峽灌區(qū)消耗黃河水資源巨大，且青銅峽至石嘴山區(qū)間有明顯的灌區(qū)退水和地下水補(bǔ)給現(xiàn)象。20世紀(jì)90年代以后的平均徑流量、平均輸沙量均低于多年平均，表明近些年黃河上游的徑流量和輸沙量減少顯著。

4 討論

黃河徑流量與輸沙量變化是由氣候變化和人類活動(dòng)導(dǎo)致的下墊面變化等諸多因素影響的，目前研究發(fā)現(xiàn)影響黃河水沙變化的原因多集中在降水量減少、蒸發(fā)增加、水利工程調(diào)節(jié)、灌區(qū)引用水、風(fēng)沙入黃和水土保持措施等方面。

近些年黃河流域水沙量均呈現(xiàn)顯著減小趨勢(shì)，水沙變化的周期性受到降水的年際變化以及氣候系統(tǒng)的影響具有相似的周期(0.5～1、3～5和7～9年)[14]。黃河上游與中游水沙變化的突變點(diǎn)多集中在1968和1986年，這主要是受上游大型水庫(kù)建設(shè)的影響。對(duì)比而言，中游的部分水文站點(diǎn)顯示，徑流量和輸沙量突變年份集中在1979與1996年前后，這與黃河中游的水土保持措施實(shí)施關(guān)系更為密切。

水庫(kù)調(diào)節(jié)使得年內(nèi)徑流量重新分配，在汛期具有攔洪錯(cuò)峰的作用，能夠有效削減洪峰流量，延滯洪水下泄，在攔蓄入庫(kù)泥沙的同時(shí)也大大減小出庫(kù)水流的挾沙能力[17]。上游徑流量和輸沙量發(fā)生突變的時(shí)間點(diǎn)與劉家峽水庫(kù)、龍羊峽水庫(kù)和青銅峽水庫(kù)建設(shè)運(yùn)行的時(shí)間相契合，這也進(jìn)一步驗(yàn)證了水庫(kù)對(duì)水沙的調(diào)控作用。

黃河上中游途徑干旱半干旱地區(qū)，沿途農(nóng)業(yè)灌溉、工業(yè)用水、居民生活用水多依賴于黃河，灌區(qū)引水引沙是黃河水沙情勢(shì)變化的又一重要因素。寧蒙灌區(qū)是黃河最主要的灌溉引水區(qū)，包括衛(wèi)寧灌區(qū)、青銅峽灌區(qū)和河套灌區(qū)等多片區(qū)域。2000—2012年寧蒙灌區(qū)年均引水132.67億m3，凈引水91.37億m3，年均灌溉引沙量達(dá)到0.3億t[20]，年均凈引水量占同期下河沿站和頭道拐站年均徑流量的35.4%和56.4%，年均引沙量占同期下河沿站和頭道拐站年均輸沙量的70.9%和68.5%。

黃土高原的水土保持建設(shè)對(duì)上游部分流域及中游流域水沙影響顯著。各時(shí)期大規(guī)模的水土保持措施建設(shè)，如植被恢復(fù)、退耕還林(草)、梯田、淤地壩和風(fēng)沙治理等極大地改變了地表產(chǎn)匯水、輸沙過程，減沙阻沙效應(yīng)明顯。上世紀(jì)70、80年代，黃土高原修建了大量的淤地壩有效的攔蓄了溝道來(lái)水來(lái)沙；而自1999年始，以退耕還林草為主的植被恢復(fù)顯著改善了區(qū)域的植被條件，有效的控制了坡面水土流失。

5 結(jié)論

筆者采用MK方法、Spearman方法、累積距平法和小波變換等方法，分析了黃河上游干流唐乃亥至頭道拐區(qū)間10個(gè)水文站的徑流量、輸沙量在月、年、年代間的變化特征，主要結(jié)論如下。

1)近60年來(lái)，黃河上游干流大部分水文站的徑流和輸沙量均呈現(xiàn)減少趨勢(shì)，且大部分站點(diǎn)減少趨勢(shì)顯著，唐乃亥站水沙減少趨勢(shì)不顯著。黃河上游干流各站點(diǎn)年徑流量和年輸沙量突變年份多發(fā)生在1968、1985、1986和1989年，與上游大型水利工程(龍羊峽、劉家峽、青銅峽等)開始運(yùn)行的時(shí)間相一致，表明干流水利工程對(duì)河流水沙具有極大的調(diào)控作用。

2)黃河上游干流各站點(diǎn)月徑流量和月輸沙量具有0.5～1、3～5和7～9年的振蕩周期，主要是由降水的季節(jié)變化和影響西北地區(qū)天氣的氣候系統(tǒng)引起的，在20世紀(jì)70年代之后，干流水利水保工程和退耕還林(草)等顯著改變了徑流輸沙的周期特征。

3)黃河上游支流匯水使得徑流量沿程逐漸增加且在下河沿站達(dá)到最大值，而寧蒙灌區(qū)引用水使得徑流量在下河沿站之后出現(xiàn)減少的狀況；受入黃泥沙貢獻(xiàn)的影響，上游輸沙量沿程逐漸增加且在下河沿站劇增，但在寧蒙河段輸沙量變化不大。

4)大型水利工程、灌區(qū)引用水、部分支流的水土保持措施等是黃河上游唐乃亥至頭道拐區(qū)間水沙變化的主要原因。