郭宏偉，徐海量，趙新風，凌紅波

(1. 中國科學院新疆生態(tài)與地理研究所，新疆 烏魯木齊 830011; 2. 中國科學院大學，北京 100049 )

塔里木河流域最大灌溉面積與超載情況探討*

郭宏偉1,2，徐海量1，趙新風1，凌紅波1

(1. 中國科學院新疆生態(tài)與地理研究所，新疆 烏魯木齊 830011; 2. 中國科學院大學，北京 100049 )

基于塔里木河流域水資源量及利用水平資料，結(jié)合《塔里木河流域工程與非工程措施五年實施方案》，通過Mann-kendall單調(diào)趨勢檢驗及未來趨勢變化Hurst指數(shù)檢驗分析了近20 a來塔里木河流域耕地面積變化特征，計算探討了該流域現(xiàn)狀水平年(2010年)以及規(guī)劃年(2020年)可承載的最大灌溉面積、現(xiàn)狀水平年灌溉面積超載情況、不同來水頻率下規(guī)劃年規(guī)劃灌溉面積超載和生態(tài)水保證情況。結(jié)果表明：① 1990-2010年流域耕地面積顯著增加，且在未來一段時間將持續(xù)增加。② 現(xiàn)狀水平年流域可承載最大灌溉面積為129.06萬hm2，除和田河外其它均超載，超載面積共計43.33萬hm2。③ 規(guī)劃年在25%來水頻率下，除了開都-孔雀河流域其他均不超載，但是相等規(guī)劃灌溉面積在枯(特枯)水期會出現(xiàn)超載現(xiàn)象，規(guī)劃灌溉面積的制定還要參照枯(特枯)水期流域能承載的最大灌溉面積；在50%、75%及90%來水頻率下，流域規(guī)劃灌溉面積超載嚴重，管理部門應認識到灌溉面積規(guī)劃中存在的問題，提出更加科學的規(guī)劃方案。④ 在25%、50%、75%及90%等4種來水頻率下2020年流域生態(tài)水分別可保證95.7%、81.2%、61.4%、44.2%，表明在規(guī)劃年的枯水期生態(tài)水供應有一定壓力。本研究為遏制塔里木河流域自然生態(tài)和環(huán)境惡化趨勢提供了基礎依據(jù)，對流域以水資源量來定人工綠洲規(guī)模、優(yōu)化水資源利用規(guī)劃與配水制度意義重大。

塔里木河流域；最大灌溉面積；超載面積；不同來水頻率

塔里木河流域是中國乃至世界生態(tài)最為脆弱的地區(qū)之一。流域氣候干旱，降水稀少，屬于典型的暖溫帶大陸性氣候，年均氣溫10.6～11.5 ℃，年均降水量大約為116.8 mm[1-2]。近年來降水則在流域源流區(qū)及干流都有增加趨勢。流域水資源的形成以冰雪融水補給為主，并有降雨徑流加入，各河流均以河流出山口為界，出山口以上為徑流形成區(qū)，徑流量沿程遞增；河流出山以后，沿程滲漏、蒸發(fā)，用于灌溉、流入湖泊或盆地，徑流量沿程遞減，最后消失于湖泊、灌區(qū)或沙漠中[3-4]。流域多年平均水資源總量為430.2億m3。自20世紀80年代末尤其進入21世紀以來，流域各源流普遍進入豐水期，阿克蘇河、和田河和葉爾羌河近10 a平均來水量比多年平均來水量增加了27.45億m3，對干流下游生態(tài)輸水影響較大的開都河近10年的水量也比多年平均高出0.74億m3。然而，由于源流區(qū)灌溉面積增加了一倍多，年用水量平均約增加50億m3，導致補給干流的總水量并未有明顯變化[5-6]。因此，在現(xiàn)有的水土資源擴張開發(fā)模式下，降水與源流來水的增多對干流來水增加的影響并不大。

從20世紀50年代開始，在源流阿克蘇河、和田河、葉爾羌河，甚至在塔里木河干流，開發(fā)了大規(guī)模的灌溉農(nóng)業(yè)。過度的水土資源開發(fā)利用，特別是1972年下游大西海子攔河水庫建成以后，直接導致河道斷流。為了保護流域的生態(tài)環(huán)境，拯救斷流近30 a瀕于毀滅的下游“綠色走廊”，中國政府投資107億元，于2001年啟動了塔里木河流域近期綜合治理工程。如今，流域生態(tài)環(huán)境得到初步改善，近期綜合治理雖取得了明顯成效，但這與近10 a和田河、葉爾羌河、阿克蘇河和開都-孔雀河(以下簡稱“開-孔河”)4條源流天然來水持續(xù)偏豐的客觀原因是密切相關的；同時，前期治理從2000-2009年實施后，源流通過常規(guī)和高科技節(jié)水增加11.4億m3，開發(fā)地下水1.4億m3，合計12.8億 m3，但是這些水量并沒有輸往干流，而是全部被開荒占用[7]，耕地的大規(guī)模、無序、快速擴張毀壞了大面積的天然植被并擠占了大量的生態(tài)用水，從20世紀70年代初開始，下游300余km的河道斷流近30 a，過度無序的水土開發(fā)、嚴重惡化的生態(tài)環(huán)境和由此引發(fā)的各種社會矛盾，引起了國內(nèi)外的廣泛關注[8]，厘清近年來各流域耕地面積變化情況、以水資源為主線計算并分析規(guī)劃年各流域能承載的最大灌溉面積及超載情況等相關研究的開展勢在必行。然而，由于流域面積巨大，各流域耕地耗水情況各不相同且十分復雜，國內(nèi)相關研究的開展甚少。

為此，本文基于塔里木河流域水資源量及水資源利用水平資料，結(jié)合《塔里木河流域工程與非工程措施五年實施方案》(2002年)(以下簡稱《五年實施方案》)[9]，通過Mann-kendall單調(diào)趨勢檢驗及未來趨勢變化Hurst指數(shù)檢驗分析了近20 a來塔里木河流域耕地面積變化特征，計算探討了流域現(xiàn)狀水平年(2010年)以及規(guī)劃年(2020年)可承載的最大灌溉面積、現(xiàn)狀水平年灌溉面積超載情況、不同來水頻率下規(guī)劃年規(guī)劃灌溉面積超載和生態(tài)水保證情況。其目的是為遏制流域自然生態(tài)和環(huán)境惡化趨勢提供基礎依據(jù)，這對塔里木河流域以水資源量來定人工綠洲規(guī)模、優(yōu)化水資源利用規(guī)劃與配水制度意義重大。

1 研究區(qū)概況

塔里木河流域位于新疆南部，是環(huán)塔里木盆地的阿克蘇河、渭干河、迪那河、喀什噶爾河、葉爾羌河、和田河、開都河—孔雀河、克里雅河、車爾臣河等9大水系共144條河流的總稱，流域總面積102萬km2，多年平均徑流量約401億m3。流域分為干流區(qū)和源流區(qū)，葉爾羌河和阿克蘇河在肖夾克會合后始稱塔里木河，干流自身不產(chǎn)流，從肖夾克至臺特瑪湖河道全長1 321 km，其中上游河段肖夾克至英巴扎，長495 km；中游河段英巴扎至恰拉，長398 km；下游河段恰拉至臺特瑪湖，長428 km。受人類活動與氣候變化的影響，許多源流早已中斷了向干流的水量補給，目前9大向心聚流水系中，阿克蘇河、和田河、葉爾羌河和開都河-孔雀河4條源流與干流有地表水力聯(lián)系，本文以下簡稱為“四源一干”(圖1)。

圖1 塔里木河流域水系示意圖Fig.1 The water system of Tarim River Basin

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源

塔里木河流域“四源一干”不同來水頻率下的水量分配方案由《塔河流域近期綜合治理規(guī)劃》(2001)[10]和《塔里木河流域“四源一干”地表水水量分配方案》(2010)(以下簡稱《分配方案》)[11]獲得。流域綜合毛灌溉定額和2010及2020年生活、工業(yè)、牲畜規(guī)劃用水量參照了《塔里木河流域綜合治理規(guī)劃》(2012)。塔里木河“四源一干”2010年實際灌溉面積和2020年規(guī)劃面積參照了《五年實施方案》。在塔里木河“四源一干”各流域行政區(qū)劃范圍的基礎上，以1991-2011年《新疆統(tǒng)計年鑒》[12]與《新疆生產(chǎn)建設兵團統(tǒng)計年鑒》[13]為依據(jù)，統(tǒng)計出1990-2010年塔里木河流域“四源一干”耕地面積。

各流域生態(tài)需水量前人研究結(jié)果有差異，郭斌[14]分析了開-孔河流域供需水平衡，并對該流域未來水量進行模擬預測，將該流域生態(tài)需水定為20.01億m3；雷志棟等[15]認為葉爾羌河與阿克蘇河流域生態(tài)需水量分別為23.87億、13.11億m3；關洪濤等[16]在對和田河綠色走廊天然生態(tài)需水量進行估算時得出和田河流域生態(tài)需水量為7.05 m3。有關部門指出“打通沿線、穿越沙漠另需水10億m3”。本文采用以上兩項之和17.05億m3作為和田河流域的生態(tài)需水量；2011年徐海量[17]在完成《塔里木河生態(tài)需水及對策研究》研究報告時，將干流流域生態(tài)需水量定為37.4億m3。

2.2 研究方法

2.2.1 耕地面積變化趨勢檢驗 基于秩的Mann-Kendall統(tǒng)計檢驗方法常用來預測時間序列資料的長期變化趨勢，使用該方法進行趨勢檢驗時，無需事先假定樣本的統(tǒng)計分布。本文在以往研究成果[18]的基礎上通過Mann-Kendall單調(diào)趨勢檢驗近20 a來塔里木河流域的耕地面積變化趨勢；Hurst指數(shù)是定量描述時間序列長程依賴性的有效方法[19]，在水文學、經(jīng)濟學、氣候?qū)W、地質(zhì)和地球化學等領域有著廣泛的應用，近期在10年以上時間序列的植被覆蓋變化研究中應用起來[20-21]。本文采用最常用的R/S分析法[22]計算Hurst指數(shù)，以檢驗流域耕地未來趨勢變化。其過程為：

首先，對時間序列(X1，X2，X3，…，Xn)依次比較，結(jié)果記為sgnθ：

(1)

接著便可用如下公式計算出Mann-Kendall統(tǒng)計值：

(2)

式中xk、xi為要進行檢驗的隨機變量，n為所選數(shù)據(jù)序列的長度。

則與此相關的檢驗統(tǒng)計量為：

(3)

隨著n的逐漸增加,Zc很快收斂于標準化正態(tài)分布，當-Z1-a/2≤Zc≤Z1-a/2時，接受原假設，表明樣本沒有明顯變化趨勢，其中，±Z1-a/2是標準正態(tài)分布中值為1-a/2時對應的顯著性水平a下的統(tǒng)計值。當統(tǒng)計量Zc為正值，說明序列有上升趨勢；Zc為負值，則表示有下降趨勢。

在Mann-Kendall檢驗中，常用β表示時間序列變化率，其公式為：

?j

(4)

式中，1

再在以上基礎上進行基于R/S分析的Hurst指數(shù)計算，其基本步驟如下:

對于一個時間序列x(t),t=1,2…。均值序列為:

(5)

累積離差為

(6)

極差為

τ=1,2…

(7)

標準差為

τ=1,2…

(8)

Mandelbrotetal證實了Hurst的研究,并得出了更廣泛的指數(shù)律,即

(9)

式中H為Hurst指數(shù)(00.5時，意味著未來的變化狀況與過去一致，這個過程具有持續(xù)性。

2.2.2 最大農(nóng)業(yè)水量及最大灌溉面積 根據(jù)流域水系組成及其徑流特性、人工綠洲與天然綠洲分布格局以及生態(tài)環(huán)境的地域特征綜合分析，流域生態(tài)調(diào)度應通過“初始水權(quán)管理、源流總量控制，豐平水期確保生態(tài)流量、枯水期保證最小灌溉面積”實現(xiàn)流域水資源配置方案目標，其技術路線如圖2，以此為依據(jù)得出不同來水頻率下最大農(nóng)業(yè)水量的計算方法：

1)豐、平水年份，最大農(nóng)業(yè)水量是用區(qū)域多年平均水資源總量扣除人類生活、牲畜及工業(yè)用水量[10]等后的水量，還要留出滿足下游天然植被保護和恢復所需的水量，即：

W農(nóng)業(yè)=W總-W生活-W牲畜-W工業(yè)-W生態(tài)

(10)

2)枯水年(75%來水頻率)農(nóng)業(yè)最大水量的計算方法是：如果流域的水資源量在滿足生活、牲畜、工業(yè)用水、給干流下泄，以及生態(tài)需水量后，留給農(nóng)業(yè)灌溉的水量大于《五年實施方案》規(guī)劃面積所需水量，此時流域最大農(nóng)業(yè)水量的計算同豐、平水期的計算方法一樣；如果流域水資源總量減去以上各用水量后，剩余的水量小于《五年實施方案》規(guī)劃面積用水時，將《五年實施方案》規(guī)劃面積用水量的90%作為農(nóng)業(yè)用水。

3)特枯水年(90%來水頻率)農(nóng)業(yè)最大水量的計算方法是：如果流域的水資源量在滿足生活、牲畜、工業(yè)用水、給干流下泄，以及生態(tài)需水量后，留給農(nóng)業(yè)灌溉的水量大于《五年實施方案》[9]規(guī)劃面積所需水量，此時流域最大農(nóng)業(yè)水量的計算同豐、平水期的計算方法一樣；如果流域水資源總量減去以上各用水量后，剩余的水量小于《五年實施方案》規(guī)劃面積用水時，將《五年實施方案》[9]規(guī)劃面積用水量的75%作為農(nóng)業(yè)用水。

圖2 塔里木河流域水資源配置技術路線圖Fig.2 Technical route of water resources allocation in Tarim River Basin

流域可承載的最大灌溉面積在最大農(nóng)業(yè)水量的基礎上通過綜合毛灌溉定額計算。其計算公式為：最大灌溉面積=最大農(nóng)業(yè)水量/綜合毛灌溉定額。

2.2.3 生態(tài)水保證率計算 在干旱區(qū)，降水極少，因此生態(tài)用水補給主要通過河流，如果流域灌溉面積超載，一定的超載面積必定會擠占一定的生態(tài)水量，導致天然植被得不到補給而衰敗甚至消失。為了準確反映塔里木河流域“四源一干”生態(tài)水的保證情況，本文計算了流域生態(tài)水保證率，計算方法為：流域能夠用到生態(tài)的水資源量除以該流域總的生態(tài)需水量[23]。

3 結(jié)果與分析

3.1 流域耕地面積變化趨勢

地域等因素造成塔里木河流域水資源分布不均，而水資源是綠洲面積的決定性因素，那么耕地的面積變化也會隨水分條件差異呈現(xiàn)不同的特點。統(tǒng)計資料顯示，塔里木河“四源一干”中，葉爾羌河、阿克蘇河流域是河川徑流量最大的流域，同時也是擁有耕地面積最大的兩大流域。圖3顯示了近20 a“四源一干”耕地面積變化情況。

圖3 1990-2010年塔里木河“四源一干”耕地面積變化Fig.3 The area change of cultivated land of four source streams and the mainstream of Tarim river in 1990-2010

經(jīng)Mann-Kendall單調(diào)趨勢檢驗，近20 a來，塔里木河流域“四源一干”耕地面積呈顯著增加的趨勢(表1)。其中，增加最快的流域是葉爾羌河流域，以平均每年0.96萬hm2的速度增加，開-孔河和阿克蘇河流域增加速度分別為每年0.93萬hm2和0.50萬hm2。對塔里木河流域耕地面積計算時間序列的Hurst指數(shù)。結(jié)果顯示，塔里木河“四源一干”各流域耕地面積在未來一段時間將呈持續(xù)增加的態(tài)勢。按目前的趨勢發(fā)展下去，2020年塔里木河流域耕地面積將遠遠超出流域承載力。

表1 塔里木河流域“四源一干”耕地面積變化Mann-Kendall趨勢檢驗及未來趨勢檢驗1)Table1 The Mann-Kendall trend test and future trend test of area changes of cultivated land of four source streams and the mainstream of Tarim river basin

1)Sig.為顯著性，**表示極顯著；R表示接受原假設

3.2 現(xiàn)狀水平年最大灌溉面積及超載情況

在現(xiàn)狀水平年，流域處于平水期，最大農(nóng)業(yè)水量按照豐、平水期最大農(nóng)業(yè)水量計算方法計算。通過計算，塔里木河流域“四源一干”在可承載能力水平下各流域最大的灌溉面積如表2所示。

比較 “四源一干”2010年實際灌溉面積與其可承載的最大灌溉面積(表3)，2010年塔里木河流域“四源一干”中除和田河流域外，其余流域灌溉面積均處于超載水平。其中開-孔河灌溉面積超載最為嚴重，超載百分比為85.3%。

表2 塔里木河“四源一干”2010年可承載的最大灌溉面積Table 2 The maximum irrigation area of four source streams and the mainstream of Tarim river basin in 2010

表3 塔里木河“四源一干”2010年耕地超載情況Table 3 Cultivated land overload situations of four source streams and the mainstream of Tarim river basin in 2010

近20 a(1990-2010年)塔里木河“四源一干”耕地面積不斷增加，按這一趨勢發(fā)展下去，規(guī)劃年(2020年)塔里木河流域耕地面積將遠遠超出流域承載力，這無疑會導致生產(chǎn)用水搶占生態(tài)用水, 上游用水搶占下游用水破壞下游水過程的完整性被，使得綠洲退化，對社會經(jīng)濟和生態(tài)環(huán)境保護可持續(xù)協(xié)調(diào)發(fā)展產(chǎn)生嚴重不利影響[24-26]。水資源是塔里木河流域農(nóng)業(yè)發(fā)展和生態(tài)環(huán)境建設最主要的制約因素, 有關部門必須從戰(zhàn)略高度來認識塔里木河流域現(xiàn)階段水資源短缺問題與耕地面積不合理擴張之間存在的矛盾，嚴格控制灌溉面積增長,建立高效節(jié)水農(nóng)業(yè)體系,發(fā)展高效生態(tài)農(nóng)業(yè), 以扭轉(zhuǎn)塔里木河流域灌溉面積持續(xù)超載的嚴重危機[27]。

3.3 規(guī)劃年最大灌溉面積、超載及生態(tài)水保證情況

基于國家《水文基本術語和符號標準》(GB/T50095-98)[28]中河川徑流豐枯年型劃分方法，根據(jù)塔里木河多年來水量的皮爾遜Ⅲ型頻率分布曲線，當某年來水頻率小于25%時，則河流處于特豐水年；當來水頻率介于25%～50%時，河流處于偏豐水年；當來水頻率介于50%～75%時，河流處于平水年；當來水頻率大于75%時，則河流水系進入枯水年；當某年來水頻率介于75%～90%時，定義河流水系處于偏枯水年；當來水頻率大于90%時，則河流水系進入特枯水年。本文分別計算了各流域四種來水頻率下規(guī)劃年的最大灌溉面積、超載及生態(tài)水保證情況。

3.3.1 25%來水頻率 當來水頻率為25%時，流域處于是豐水期，最大農(nóng)業(yè)水量按照豐、平水期最大農(nóng)業(yè)水量計算方法計算。通過計算，2020年塔里木河流域“四源一干”除開-孔河流域規(guī)劃灌溉面積超出最大灌溉面積14.6%外，其它流域均不超載(表4)。

3.3.2 50%來水頻率 當來水頻率為50%時，流域處于平水期，最大灌溉面積的計算與豐水期一致。通過計算，2020年塔里木河流域“四源一干”除和田河流域不超載外，其他流域均超載，超載主要發(fā)生在開-孔河和塔里木河干流，分別超出其最大灌溉面積的31.6%和34.1%(表5)。

3.3.3 75%來水頻率 當來水頻率為75%時，屬于枯水年份，在計算2020年塔里木河“四源一干”能承載的最大灌溉面積時，由于和田河與阿克蘇河流域這兩個流域的水資源量在滿足生活、牲畜、工業(yè)用水、給干流下泄，以及生態(tài)需水量后，留給農(nóng)業(yè)灌溉的水量大于《五年實施方案》規(guī)劃面積所需水量，此時這兩個流域最大農(nóng)業(yè)水量的計算同豐、平水期的計算方法一樣，得出和田河與阿克蘇河流域最大農(nóng)業(yè)用水分別為18.6億、42.7億m3。在計算開-孔河、葉爾羌河與塔里木河干流最大灌溉水量時，由于這3個流域水資源總量減去生活、工業(yè)、牲畜用水、給干流下泄水量，以及生態(tài)需水量后，留給農(nóng)業(yè)灌溉的水量小于《五年實施方案》規(guī)劃面積用水，依據(jù)枯水期調(diào)水原則，這三個流域75%來水頻率下的最大農(nóng)業(yè)水量為《五年實施方案》規(guī)劃面積用水量的90%，得出開-孔河、葉爾羌河與干流的最大灌溉水量分別為14.64億、35.03億、8.47億m3(表6)。

表4 塔里木河“四源一干”2020年25%來水頻率時最大灌溉面積及規(guī)劃灌溉面積超載情況Table 4 The maximum irrigation area and planning irrigation area overload situations of four source streams and the mainstream of Tarim river basin in 25% runoff frequency in 2020

表5 塔里木河“四源一干”2020年50%來水頻率時最大灌溉面積及規(guī)劃灌溉面積超載情況Table 5 The maximum irrigation area and planning irrigation area overload situations of four source streams and the mainstream of Tarim river basin in 50% runoff frequency in 2020

在75%來水頻率下，2020年塔里木河“四源一干”規(guī)劃灌溉面積與最大灌溉面積相比均超載，超載主要發(fā)生在塔里木河干流、葉爾羌河和開-孔河，分別超出其最大灌溉面積的54.8%、36.4%和34.2%。

表6 2020年75%來水頻率時塔里木河“四源一干”最大灌溉面積、規(guī)劃灌溉面積超載及生態(tài)水保證情況Table 6 The maximum irrigation area, planning irrigation area overload situations and ecological water guarantee situations of four source streams and the mainstream of Tarim river basin in 75% runoff frequency in 2020

在75%來水頻率下，生態(tài)可能受影響，因為耕地超載引起生態(tài)水擠占現(xiàn)象，造成生態(tài)水虧缺。和田河流域超載面積最少，能保證生態(tài)水的92.9%，葉爾羌河流域超載面積最多，僅能保證生態(tài)水的46.5%。經(jīng)合計，2020年在75%來水頻率下塔里木河流域“四源一干”26.78億m3的生態(tài)水將被擠占。

3.3.4 90%來水頻率 當來水頻率為90%時，屬于枯水年份，在計算2020年塔里木河“四源一干”能承載的最大灌溉面積時，由于和田河與阿克蘇河流域的水資源量在滿足生活、工業(yè)、牲畜用水、給干流下泄水量，以及生態(tài)需水量后，留給農(nóng)業(yè)灌溉的水量大于《五年實施方案》規(guī)劃面積用水(表12)，在計算這兩個流域90%來水頻率下最大農(nóng)業(yè)水量時，同豐、平水期的計算方法一樣，得出和田河與阿克蘇河流域最大農(nóng)業(yè)用水分別為17.4億m3、35.77億m3。在計算開-孔河、葉爾羌河與干流最大灌溉水量時，由于這3個流域水資源總量在扣除生活、牲畜、工業(yè)用水、給干流下泄水量，以及生態(tài)需水量后，留給農(nóng)業(yè)灌溉的水量小于《五年實施方案》規(guī)劃面積用水量，據(jù)枯水期調(diào)水原則，這3個流域90%來水頻率下的最大灌溉水量為《五年實施方案》規(guī)劃面積用水量的75%，得出開-孔河、葉爾羌河與干流的最大灌溉水量分別為12.2 億m3,29.19億m3,7.06億m3(表7)。

在90%來水頻率下，2020年塔里木河流域“四源一干”灌溉面積與規(guī)劃面積相比均超載，干流超載最為嚴重，超載比例為85.8%，開-孔河和葉爾羌河超載較為嚴重，超載百分比分別為61.0%和63.7%。

表7 塔里木河“四源一干”2020年90%來水頻率時最大灌溉面積、規(guī)劃灌溉面積超載及生態(tài)水保證情況Table 7 The maximum irrigation area, planning irrigation area overload situations and ecological water guarantee situations of four source streams and the mainstream of Tarim river basin in 90% runoff frequency in 2020

在90%來水頻率下，和田河流域超載面積最少，能保證生態(tài)水的85.9%，葉爾羌河流域超載面積多，僅能保證生態(tài)水的22.1%。經(jīng)合計，2020年在90%來水頻率下塔里木河流域“四源一干”44.59億m3的生態(tài)水將被擠占。

綜上所述，在規(guī)劃年，當來水頻率為25%時，塔里木河流域“四源一干”中除開-孔河流域規(guī)劃面積超出了最大灌溉面積的15.0%外，其它均不超載，在50%、75%和90%來水頻率下，規(guī)劃灌溉面積分別超載13.49萬 hm2，28.05萬 hm2和46.98萬 hm2。這表明豐水期流域水量基本可以滿足農(nóng)業(yè)灌溉，規(guī)劃灌溉面積在合理范圍之內(nèi)，但是相等規(guī)劃灌溉面積在枯(特枯)水期會出現(xiàn)超載現(xiàn)象，所以規(guī)劃灌溉面積的制定還要參照其他來水頻率下流域所能承載的最大灌溉面積；規(guī)劃灌溉面積的嚴重超載勢必會導致流域內(nèi)耕地的大規(guī)模、快速擴張，使得大面積的天然植被被毀壞，不利于流域水土資源合理配置，管理部分應該充分認識到現(xiàn)階段規(guī)劃灌溉面積制定中存在的問題，在遵守自然規(guī)律的基礎上提出更加科學、準確的流域規(guī)劃方案，進而指導流域良性發(fā)展[29]。在75%和90%來水頻率下，即枯水期，流域生態(tài)水保證率過低，表明在規(guī)劃年枯水期生態(tài)水的供應有一定壓力，而持續(xù)枯水期引發(fā)輸往下游水量減少，勢必會導致下游天然植被得不到灌溉而死亡[30]，所以要在規(guī)劃灌溉面積的制定過程中參照枯(特枯)水期流域所能承載的最大灌溉面積，同時對現(xiàn)階段枯(特枯)水期生態(tài)水供應就已經(jīng)困難的區(qū)域進行強行退地，以遏制流域生態(tài)環(huán)境持續(xù)惡化。

4 結(jié) 論

通過以上分析與研究，可以得出以下幾點結(jié)論：

1)近20 a(1990-2010年)塔里木河“四源一干”耕地面積呈顯著增加的趨勢：耕地面積增加最快的流域是葉爾羌河流域，以平均每年0.96萬 hm2的速度增加，開-孔河和阿克蘇河流域增加速度分別為每年0.93萬hm2和0.50萬hm2。

2)現(xiàn)狀水平年塔里木河“四源一干”能承載的最大灌溉面積為129.06萬 hm2，除了和田河流域灌溉面積不超載外，其它流域均超載，總計超載面積為43.33萬 hm2，其中葉爾羌河、阿克蘇河、開-孔河、干流河流域分別超載8.77、14.04、16.59、3.93萬 hm2，超載比例分別為21.1%、35.0%、85.3%、64.6%。

3)規(guī)劃年在25%來水頻率下，塔里木河流域“四源一干”能承載的最大的灌溉面積是167.80萬 hm2，除了開-孔河流域規(guī)劃面積超出了最大灌溉面積的15.0%外，其它均不超載。在50%來水頻率下，可承載的最大灌溉面積為144.78萬 hm2，開-孔河、干流、葉爾羌河流域規(guī)劃灌溉面積分別超出其最大灌溉面積的31.6%、34.1%、10.7%；在75%來水頻率下，可承載的最大灌溉面積為126.94萬 hm2，開-孔河、干流、葉爾羌河流域規(guī)劃灌溉面積分別超出最大灌溉面積的34.2%、54.8%、36.4%，阿克蘇河與和田河流域超載比例較小，分別超出最大灌溉面積的7.4%、6.5%；在90%來水頻率下，可承載的最大灌溉面積為108.01萬 hm2，開-孔河、干流、葉爾羌河、阿克蘇河、和田河流域灌溉面積分別超出最大灌溉面積61.0%、85.8%、63.7%、28.2%、13.8%。

4)隨來水量的減少，生態(tài)水被擠占現(xiàn)象更加嚴重，在25%、50%、75%、90%四種來水頻率下，2020年生態(tài)水分別可保證95.7%、81.2%、61.4%、44.2%。

[1] CHEN Y N, LI W H, XU C C, et al. Effects of climate change on water resources in Tarim River Basin, Northwest China[J]. Journal of Environmental Sciences, 2007, 19(4):488-493.

[2] CHEN Y N, XU Z X. Plausible impact of global climate change on water resources in the Tarim River Basin[J]. Science China(Earth Sciences), 2005, 48(1):127-133.

[3] HARTMANN H, SNOW J A, STEIN S, et al. Predictors of precipitation for improved water resources management in the Tarim River basin[J]. Journal of Arid Environments, 2016,125:31-42.

[4] 李紅軍，姜志紅，魏文壽.近40年來塔里木河流域旱澇的氣候變化[J].地球科學,2007, 27(6):1216-1222. LI H J, JIANG Z H, WEI W T. Drought and flood change of tarim river basin in recent 40 years [J]. Scientia Geographica Sinica,2007, 27(6):1216-1222.

[5] XU J H, CHEN Y N, LI W H, et al. The complex nonlinear systems with fractal as well as chaotic dynamics of annual runoff processes in the three headwaters of the Tarim River[J]. Journal of Geographical Sciences, 2009, 19(1):1227-1239.

[6] CUI Z, CHEN Y, ZHAO R. Responses of surface runoff to climate change and human activities in the arid region of central asia: a case study in the Tarim River Basin, China[J]. Environmental Management, 2013, 51(4):327-335.

[7] 鄧銘江．中國塔里木河治水理論與實踐[M]. 北京:科學出版社,2009. DENG M J．Chinese Tarim river flood control theory and practice [M]. Beijing: Science Press,2009.

[8] 饒瑞符. 塔里木河的變遷與整治[C]//毛德華. 塔里木河流域水資源、環(huán)境與管理. 北京:中國環(huán)境科學出版社,1998:184-193. RAO R F. Changes and regulation of the Tarim River [C]//MAO D H. Water resources, environment and management in Tarim River Basin. Beijing : China Environmental Science Press,1998:184-193.

[9] 新疆塔里木河流域管理局.塔里木河流域工程與非工程措施五年實施方案[R]. 2002. Xinjiang Tarim River Basin Management Bureau. Five year implementation plan of engineering and non engineering measures in Tarim River Basin [R]. 2002.

[10] 新疆塔里木河流域管理局.塔河流域近期綜合治理規(guī)劃[R]. 2001. Xinjiang Tarim River Basin Management Bureau. Comprehensive harnessing plan of Tarim River Basin in the near future [R]. 2001.

[11] 新疆塔里木河流域管理局.塔里木河流域“四源一干”地表水水量分配方案[R].2010. Xinjiang Tarim River Basin Management Bureau. Surface water allocation scheme of the four source streams and the main stream of the The Tarim River Basin [R]. 2010.

[12] 新疆維吾爾自治區(qū)統(tǒng)計局.新疆統(tǒng)計年鑒[M].北京：中國統(tǒng)計出版社,2013. Statistics bureau of xinjiang uygur autonomous region. Xinjiang statistical yearbook [M].Beijing：China Statistics Press,2013.

[13] 新疆生產(chǎn)建設兵團統(tǒng)計局，國家統(tǒng)計局兵團調(diào)查總隊.新疆生產(chǎn)建設兵團統(tǒng)計年鑒[M]. 北京：中國統(tǒng)計出版社,2013. Statistics bureau of xinjiang production and construction corps，National bureau of statistics corps investigation corps. Statistical yearbook of Xinjiang production and construction corps [M]. Beijing：China Statistics Press,2013.

[14] 郭斌.開都-孔雀河流域供需水平衡與模擬預測[D].烏魯木齊：中國科學院新疆生態(tài)與地理研究所,2013. GUO B. Water balance and simulation prediction of supply and demand in Kaidu-Kongque River Basin [D].Urumqi：Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences,2013.

[15] 雷志棟,胡和平,楊詩秀,等.塔里木盆地綠洲耗水分析[J].水利學報,2006,37(12):1470-1475. LEI Z D,HU H P,YANG S X, et al. Analysis on water consumption in oases of the Tarim Basin [J]. Journal of Hydraulic Engineering,2006,37(12):1470-1475.

[16] 關洪濤,楊澤慧,宋建鵬,等.和田河綠色走廊天然生態(tài)需水量估算[J].水資源與水工程學報,2005,16(1):71-75. GUAN H T, YANG Z H, SONG J P, et al. Estimation of the requiring water quantity for maintaining the natural vegetation in Hetianhe river green corridor [J]. Journal of Water Resources and Water Engineering,2005,16(1):71-75.

[17] 徐海量,白元,凌紅波,等.塔里木河生態(tài)需水及對策研究[M].烏魯木齊：新疆科學出版社,2011. XU H L, BAI Y, LING H B, et al. Study on ecological water demand of Tarim River and Its Countermeasures [M].Urumqi：Xinjiang Science Press,2011.

[18] HAMED K H, RAO A R. A modified Mann-Kendall trend test for autocorrelated data[J]. Journal of Hydrology, 1998, 204(1):182-196.

[19] 侯西勇，應蘭蘭，高猛，等.1998-2008年中國東部沿海植被覆蓋變化研究[J].地理科學, 2010,30(5):735-741. HOU X Y, YING L L, GAO J, et al. Study on vegetation cover change in eastern coastal areas of China during 1998-2008 [J]. Scientia Geographica Sinica, 2010,30(5):735-741.

[20] 袁麗華，蔣緯國，申文明，等.2000-2010年黃河流域植被覆蓋的時空變化[J].生態(tài)學報,2013,33(24):7798-7806. YUAN L H, JIANG W G, SHEN W M, et al. The spatio-temporal variations of vegetation cover in the Yellow River Basin from 2000 to 2010 [J]. Acta Ecologica Sinica,2013,33(24):7798-7806.

[21] 范娜，謝高地，張昌順，等.2001年至2010年瀾滄江流域植被覆蓋動態(tài)變化分析[J].資源科學,2012,34(7)：1222-1231. FAN N, XIE G D, ZHANG C S, et al. Spatial-temporal dynamic changes of vegetation cover in Lancang River Basin during 2001-2010 [J]. Resources Science,2012,34(7)：1222-1231.

[22] 王桂剛，周克法，孫莉，等.近10a新疆地區(qū)植被動態(tài)與R/S分析[J].遙感技術與應用,2010,25(1)：84-90. WANG G G, ZHOU K F, SUN L, et al. Vegetation dynamics and R/S analysis of 10 a in Xinjiang area [J]. Remote Sensing Technology and Application,2010,25(1)：84-90.

[23] 葉茂. 塔里木河流域環(huán)境承載力研究[D].烏魯木齊：新疆師范大學,2014. YE M. Study on environmental carrying capacity of Tarim River Basin [D].Urumqi：Xinjiang Normal University,2014.

[24] 滿蘇爾·沙比提,努爾卡木里·玉素甫.塔里木河流域綠洲耕地變化及其河流水文效應[J].地理研究,2010,29(12):2251-2260. MANSUR S, NURKAMIL Y. Oasis land use change and its hydrological response to Tarim River Basin [J]. Geographical Research,2010,29(12):2251-2260.

[25] JOSEPH S. Sustainable irrigation Agriculture in the Tarim River Basin of Western Xinjiang Province-A proposal for the establishment of a salinity laboratory of P. R. China[J].Arid Land Geography,2000,23(2):174-176.

[26] 徐海量,葉茂,宋郁東,等.塔里木河流域水文過程的特點初探[J].水土保持學報,2005,19(2):156-170. XU H L, YE M, SONG Y D, et al. Characteristic analysis on hydrological processes of Tarim River Basin [J]. Journal of Soil and Water Conservation,2005,19(2):156-170.

[27] 孟麗紅,陳亞寧,李衛(wèi)紅.新疆塔里木河流域水資源承載力評價研究[J].中國沙漠,2008,28(1):185-190. MENG L H, CHEN Y N, LI W H. Study of water resources carrying capacity in Tarim River Basin of Xinjiang [J]. Journal of Desert Research,2008,28(1):185-190.

[28] 中華人民共和國水利部.水文基本術語和符號標準(GB/T50095-98)[M].1999. Ministry of Water Resources of the PRC. Standard for essential technical terms and symbols in hydrology (GB/T50095-98)[M]. 1999.

[29] 高前兆,樊自立.塔里木河流域的環(huán)境治理與水土保持生態(tài)建設[J].水土保持學報,2002,16(1):11-15. GAO Q Z, FAN Z L. Environment rehabilitation and ecological construction of soil and water conservation in Tarim River Basin [J]. Journal of Soil and Water Conservation,2002,16(1):11-15.

[30] 儲國強,劉嘉麒,孫青,等.新疆克里雅河洪泛事件與樹輪記錄的初步研究[J].第四紀研究,2002,37(6):107-111. CHU G Q, LIU J Q, SUN Q, et al. Preliminary research on the flood events based on the studies of tree ring width (PopulusEuphratica) in the Keriya river, Xinjiang [J]. Quaternary Sciences,2002,37(6):107-111.

The maximum irrigation area and overload situations of the mainstream of Tarim river basin

GUOHongwei1,2,XUHailiang1,ZHAOXinfeng1,LINGHongbo1

(1.Xinjiang Institute of Ecology and Geography, Chinese Academy of Science, Urumqi 831400,China; 2. University of Chinese Academy of Sciences, Beijing 100049, China)

The changes of cultivated land area in the Tarim River Basin in recent years were analyzed with Mann-kendall monotonic trend test and future trend change Hurst index based on water resources quantity and water resources utilization level data in Tarim River Basin and combined with theMeasuresofTarimRiverBasinProjectandNon-projectImplementationPlanforFiveYear. The following factors were calculated for four source streams and mainstream of the Tarim River Basin: maximum irrigation area in 2010 and 2020, irrigation area overload situations in 2010, planning irrigation area overload situations and ecological water guarantee situations at different runoff frequencies in 2020. The results showed that: ① The cultivated land area increased significantly during 1990-2010, and it will keep increasing in the following years. ② The maximum irrigation area of the Tarim River basin in 2010 is 129.06×104hm2. All basins except for the Hotan River Basin are overloaded, and total overload area is 43.33×104hm2.③ In 2020, all basins except for Kai-Kong River Basin are not overloaded at 25% runoff frequency, but equal planning irrigation area will be overloaded in dry season. The formulation of planning irrigation area should refer to the maximum irrigation area that the basin can carry in dry season; planning irrigation area of the basin is seriously overloaded at 50%, 75% and 90% runoff frequency, respectively, so management department should recognize the problems in planning irrigation area and put forward more scientific planning schemes. ④ Guarantee rates of ecological water in the basin at water frequencies of 25%, 50%, 75% and 90% are 95.7%,81.2%,61.4%,44.2%, respectively, suggesting that in the planning year there is a certain pressure on the ecological water supply during the dry season. This study provides a basis for the containment of the ecological and environmental deterioration of the Tarim River Basin, and it is also significant for the determination of artificial oasis scale of the basin according to the water resources quantity, optimization of water resource utilization planning and water distribution.

Tarim River Basin; maximum irrigation area; overloading problem; different runoff frequencies

2016-03-31 基金項目：國家自然科學基金(31400466,30970549,41471099)；水利部公益性行業(yè)科研專項基金(201501060)

郭宏偉(1991年生),男;研究方向：恢復生態(tài)學；E-mail:ghwxzsy@163.com

趙新風(1981年生),女;研究方向：恢復生態(tài)學；E-mail:zxinfeng668@sina.com

10.13471/j.cnki.acta.snus.2017.02.021

S274

A

0529-6579(2017)02-0140-11