陶佳輝，卞錦宇，敖天其，呂良華，耿雷華，王 歡

(1.南京水利科學(xué)研究院水文水資源與水利工程科學(xué)國家重點實驗室，江蘇 南京 210029；2.四川大學(xué)水利水電學(xué)院，四川 成都 610065)

2002年8月，國務(wù)院審議通過了《南水北調(diào)工程總體規(guī)劃》(以下簡稱《總體規(guī)劃》)，自此南水北調(diào)工程由規(guī)劃階段進入實施階段。南水北調(diào)工程是實現(xiàn)水資源南北調(diào)配、東西互濟，提高水資源與人口、經(jīng)濟時空匹配性的國家戰(zhàn)略性工程。全球氣候變化改變區(qū)域水循環(huán)條件[1]，使南水北調(diào)工程所涉長江和黃淮海流域來水條件發(fā)生了動態(tài)調(diào)整，改變了調(diào)水區(qū)可調(diào)水量和受水區(qū)水資源需求，從而影響調(diào)水規(guī)模和工程運行方式[2-4]。目前，學(xué)者們對南水北調(diào)工程調(diào)水區(qū)及受水區(qū)水文情勢變化特征開展了一系列研究工作[5-7]。夏軍等[8]基于Mann-Kendall(M-K)檢驗、彈性系數(shù)法和水文模擬法等方法，發(fā)現(xiàn)1961—2013年漢江上游流域徑流呈明顯下降趨勢，降水及潛在蒸發(fā)變化趨勢不明顯；左其亭等[9]從地理學(xué)的角度，對南水北調(diào)中線調(diào)水區(qū)主要特征參數(shù)進行系統(tǒng)梳理，分析了水資源特征、水資源開發(fā)利用演變過程以及存在的問題；張曉松等[10]根據(jù)東線受水區(qū)1951—2012年降水資料，分析得出江蘇省沂南水利分區(qū)降水呈不顯著減少趨勢，變化周期為27 a；方思達等[11]基于中線工程水源區(qū)與受水區(qū)歷史旱澇變化特征，分析了近500 a來的旱澇遭遇特征及調(diào)水保障概率；楊鵬鵬等[12]通過Spearman、小波分析等方法，分析得出西線調(diào)水區(qū)年徑流量變化趨勢總體不明顯，認為氣候變化是流域徑流量變化的主要影響因素。相關(guān)研究對南水北調(diào)工程的規(guī)劃與運行提供了一定的科學(xué)依據(jù)[13-15]，但《總體規(guī)劃》以2000年作為規(guī)劃的現(xiàn)狀水平年，論證可調(diào)水量和工程規(guī)模所采用的水文資料系列均截至2000年甚至更短。而長江和黃淮海四大流域降水、徑流系列已延長近20 a，并且2000年以來長江和黃淮海流域降水、徑流等氣象水文要素已檢測到一定的變異信號[16-17]。鑒于此，本文以1956—2016年為研究期，以南水北調(diào)東線工程調(diào)水區(qū)、受水區(qū)為研究范圍，以1956—2000年為基準期，以2001—2016年為比較期，基于東線工程調(diào)水區(qū)、受水區(qū)長系列降水和徑流資料，采用M-K趨勢檢驗、小波分析等定量分析方法，通過比較期與基準期降水、徑流序列統(tǒng)計特征值，分析數(shù)據(jù)延長后統(tǒng)計特性上的差異，闡明數(shù)據(jù)系列變化對南水北調(diào)工程水資源條件的影響，為確保南水北調(diào)東線工程科學(xué)建設(shè)和安全經(jīng)濟運行提供參考。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

南水北調(diào)東線工程從長江下游江都樞紐和寶應(yīng)抽水站取水，基本沿京杭大運河向北輸水。主要供水目標是解決調(diào)水線路沿線和山東半島城市及工業(yè)用水，并在北方需要時，提供農(nóng)業(yè)和生態(tài)環(huán)境用水。根據(jù)《南水北調(diào)東線工程規(guī)劃》(2001年修訂)，東線工程一期調(diào)江水量89.37億m3，二期規(guī)模148.17億m3(含一期規(guī)模)。東線一期工程已建成并于2013年11月正式通水。

圖1為南水北調(diào)東線工程調(diào)水區(qū)與受水區(qū)范圍，受水區(qū)包括淮河、海河和黃河流域的4個省25個地級城市及天津市，其中江蘇省5個(揚州、淮安、宿遷、徐州、連云港)、安徽省3個(蚌埠、淮北、宿州)、山東省14個(棗莊、濟寧、菏澤、泰安、濟南、淄博、濰坊、青島、聊城、德州、濱州、東營、威海、煙臺)，河北省2個(衡水、滄州)；調(diào)水區(qū)為長江流域，以長江干流大通水文站(集水面積170.5萬km2，占長江流域總面積的95.6%)作為分析節(jié)點，結(jié)合長江全流域研究其降水量與天然徑流量；受水區(qū)針對沿線25市研究其降水量與徑流量。

圖1 南水北調(diào)東線工程調(diào)水區(qū)和受水區(qū)范圍Fig.1 Scope of water transfer area and water receiving area of the East Route of South to North Water Transfer Project

1.2 數(shù)據(jù)資料

收集了研究區(qū)4個省25個地級城市及天津市1956—2016年逐年降水量和逐年地表徑流量、長江大通水文站1956—2016年逐月天然徑流量，結(jié)合全國水資源公報、長江水資源公報、黃河水資源公報、相關(guān)地市水資源公報等資料，作為本文的數(shù)據(jù)基礎(chǔ)。

1.3 研究方法

水文氣象序列趨勢分析主要判斷水文序列是否存在顯著趨勢特征以及該序列趨勢是上升還是下降[18]。本文采用非參數(shù)M-K檢驗法[19-20]進行趨勢分析，采用Morlet小波作為基函數(shù)進行小波變換分析[21]開展水文周期性研究，從而揭示南水北調(diào)東線工程調(diào)水區(qū)及受水區(qū)水文序列變化特征。

2 調(diào)水區(qū)變化特征

2.1 調(diào)水區(qū)降水

a.總體變化特征。調(diào)水區(qū)基準期多年平均年降水量為1 084.7 mm，比較期多年平均年降水量為1 069.2 mm，較基準期下降約1.4%，總體變化不明顯(圖2)。年降水量最大(1998年)為1 239.1 mm，最小(1978年)為922.5 mm，降水年際間有明顯差異。對基準期與比較期調(diào)水區(qū)年降水量采用皮爾遜Ⅲ型曲線進行水文頻率分析，結(jié)果見表1，表中Cv為變差系數(shù)，Cs為偏態(tài)系數(shù)，對水文序列進行頻率分析時，一般要求序列長度大于30 a，比較期時段較短，因此僅計算了Cv值。

圖2 1956—2016年調(diào)水區(qū)年降水變化Fig.2 Annual precipitation change in water transfer area from 1956 to 2016

表1 調(diào)水區(qū)水文頻率分析結(jié)果Table 1 Hydrological frequency analysis results of water transfer area

b.M-K趨勢分析。對調(diào)水區(qū)降水進行M-K趨勢檢驗結(jié)果表明，經(jīng)95%置信水平的統(tǒng)計檢驗，基準期降水序列統(tǒng)計值Z為0.73，比較期降水序列Z值為1.22，1956—2016年調(diào)水區(qū)降水序列Z值為-0.02，可見基準期、比較期、研究期3個時段調(diào)水區(qū)年降水量變化趨勢均不顯著。

c.小波周期分析。圖3為1956—2016年調(diào)水區(qū)降水量小波變換系數(shù)實部等值線，圖3中綠色部分實線為正值等值線，代表豐水期，黃色部分虛線為負值等值線，代表枯水期。由圖3可見，在4 a、13 a左右的時間尺度豐枯交替變化比較明顯；在13 a的周期尺度上，降水量經(jīng)歷了“枯-豐-枯-豐”的循環(huán)交替過程，可推測未來一個周期調(diào)水區(qū)降水量將呈上升趨勢。圖4為1956—2016年調(diào)水區(qū)降水量小波變換方差過程線，同樣顯示出在4 a、13 a左右時間尺度的極值周期較明顯。綜合分析可知，在30 a左右的時間尺度上調(diào)水區(qū)降水存在一定的周期變化規(guī)律，但具體時間尺度不明顯，這可能是由于降水序列長度不夠，限于數(shù)據(jù)原因，無法分析出具體的主周期。

圖3 1956—2016年調(diào)水區(qū)降水量小波變換系數(shù)實部等值線Fig.3 Real part isoline of wavelet transform coefficient of precipitation in water transfer area from 1956 to 2016

圖4 1956—2016年調(diào)水區(qū)降水量小波變換方差過程線Fig.4 Wavelet transform variance hydrograph of precipitation in water transfer area from 1956 to 2016

2.2 調(diào)水區(qū)徑流

a.總體變化特征。長江大通水文站1956—2016年平均天然徑流量為9 454.2億m3?；鶞势谀昃鶑搅髁繛? 501.1億m3，比較期平均徑流量為9 322.6億m3，較基準期下降1.9%(圖5)。大通站最豐年(1998年)徑流量達到13 132.6億m3，最枯年(1978年)徑流量僅為7 333.5億m3。表2為大通站水文頻率分析結(jié)果，基準期平水年徑流量為9 432.8m3，研究期平水年徑流量為9 368.7m3。

表2 大通站水文頻率分析結(jié)果Table 2 Hydrological frequency analysis results of Datong Station

圖5 1956—2016年大通站天然年徑流量變化Fig.5 Variation of natural annual runoff at Datong Station from 1956 to 2016

b.M-K趨勢分析。對調(diào)水區(qū)徑流進行M-K趨勢分析結(jié)果表明，經(jīng)95%置信水平的統(tǒng)計檢驗，基準期大通站年徑流序列Z值為2.07，比較期徑流序列Z值為0.69，研究期徑流序列Z值為1.18，可見大通站基準期年徑流量呈顯著上升趨勢，比較期年徑流量呈不顯著上升趨勢，研究期年徑流量呈一定上升趨勢但不顯著。

c.小波周期分析。圖6為1956—2016年大通站徑流量小波變換系數(shù)實部等值線，從圖6可見，大通站年徑流量序列主周期為28 a，次周期為12 a；28 a的周期尺度較為明顯，主次周期內(nèi)徑流量經(jīng)歷了“枯-豐-枯-豐”的循環(huán)交替過程，反映了大通站徑流具有明顯的周期性和波動性。同時，在28 a的時間尺度上，2000—2016年長江流域經(jīng)歷了一段徑流偏少期，降水變化周期類似，推測未來一個周期大通站將處于徑流偏多期，這為東線工程一期調(diào)水提供了穩(wěn)定的水源條件。圖7為1956—2016年大通站徑流量小波變換方差過程線，可見，在5 a、12 a、28 a左右時間尺度的極值周期較明顯。最大峰值對應(yīng)著28 a左右的時間尺度，說明28 a左右的周期震蕩最強，為大通站年徑流量的第一主周期；12 a左右時間尺度對應(yīng)著第二峰值，為第二主周期；5 a左右的時間尺度為第三主周期。

圖6 1956—2016年大通站徑流量小波變換系數(shù)實部等值線Fig.6 Real part contour of wavelet transform coefficient of runoff at Datong Station from 1956 to 2016

圖7 1956—2016年大通站徑流量小波變換方差過程線Fig.7 Wavelet transform variance hydrograph of runoff at Datong Station from 1956 to 2016

3 受水區(qū)變化特征

3.1 受水區(qū)降水

a.總體變化特征。按南水北調(diào)東線工程受水區(qū)涉及的25個地級行政區(qū)計算，受水區(qū)基準期平均降水量為699.0 mm，比較期平均降水量為703.7 mm，與基準期基本持平，總體變化不明顯(圖8)。表3為不同時段受水區(qū)降水量水文頻率分析結(jié)果，東線受水區(qū)基準期平水年降水量為691.6 mm，研究期平水年降水量為693.7 mm。根據(jù)1956—2016年受水區(qū)5個省相關(guān)地市的降水量資料，安徽省和江蘇省多年平均降水量均超過800 mm，山東省、河北省和天津市多年平均降水量為500～600 mm。與基準期相比，比較期江蘇省、安徽省、山東省受水區(qū)降水有一定增加，增長率分別達到3.7%、2.1%及3.7%；但比較期河北省和天津市受水區(qū)降水較基準期有所下降，兩市受水區(qū)合計減少60.1 mm，相對減少率分別為5.6%和5.2%。

(a) 全部受水區(qū)

(b) 江蘇受水區(qū)

(c) 安徽受水區(qū)

(d) 山東受水區(qū)

(e) 河北受水區(qū)

(f) 天津受水區(qū)圖8 1956—2016年東線受水區(qū)年降水量變化Fig.8 Variation of annual precipitation in the water-receiving area of Eastern Route from 1956 to 2016

表3 不同時段受水區(qū)降水量水文頻率分析結(jié)果Table 3 Hydrological frequency analysis results of precipitation in water receiving areas at different periods

b.M-K趨勢分析。表4為受水區(qū)降水量M-K趨勢分析結(jié)果，經(jīng)95%置信水平的統(tǒng)計檢驗，基準期受水區(qū)降水量整體呈現(xiàn)下降趨勢，其中山東省呈現(xiàn)顯著下降趨勢，其他省級行政區(qū)呈現(xiàn)不顯著下降趨勢。比較受水區(qū)各省級行政區(qū)基準期Z值大小，得到降水量下降趨勢顯著性由大到小排序為：山東省、河北省、江蘇省、天津市和安徽??；比較期受水區(qū)年降水量整體呈現(xiàn)下降趨勢，江蘇省、河北省和天津市呈不顯著上升趨勢，安徽省和山東省呈不顯著下降趨勢；研究期受水區(qū)年降水量整體呈現(xiàn)下降趨勢，各省市均呈現(xiàn)不顯著下降趨勢。

表4 受水區(qū)降水量M-K趨勢分析結(jié)果Table 4 M-K trend analysis results of precipitation in the catchment area

c.小波周期分析。圖9為受水區(qū)降水量小波轉(zhuǎn)換系數(shù)實部等值線，可見，天然年降水量表現(xiàn)為25 a、12 a左右時間尺度的周期，在該時間尺度上，年降水量的豐枯交替變化十分明顯。25 a的大周期與12 a的小周期呈現(xiàn)一定的嵌套關(guān)系，從25 a的大周期上看，降水量在該時間尺度上經(jīng)歷了3次“豐-枯”交替變化，目前東線受水區(qū)正處于年降水量偏少期。圖10為1956—2016年受水區(qū)降水量小波變換方差過程線，在4 a、12 a、17 a及25 a左右時間尺度的極值周期較明顯。最大峰值對應(yīng)著25 a左右的時間尺度，說明25 a左右的周期震蕩最強，為東線受水區(qū)天然年降水的第一主周期；17 a左右時間尺度對應(yīng)著第二峰值，為第二主周期；12 a、4 a左右的時間尺度分別為第三、第四主周期。

圖9 1956—2016年受水區(qū)降水量小波變換系數(shù)實部等值線Fig.9 Real part isoline of wavelet transform coefficient of precipitation in catchment area from 1956 to 2016

圖10 1956—2016年受水區(qū)降水量小波變換方差過程線Fig.10 Wavelet transform variance hydrograph of precipitation in watershed area from 1956 to 2016

3.2 受水區(qū)徑流

a.總體變化特征。1956—2016年受水區(qū)平均徑流量為284.4億m3?；鶞势谀昃鶑搅髁繛?81.5億m3，比較期平均徑流量為292.5億m3，較基準期略有上升，上升率為3.9%，總體變化不明顯(圖11)。最豐年(1964年)徑流量為731億m3，最枯年(2002年)徑流量為83.0億m3，極值比為8.8。表5為受水區(qū)徑流水文頻率分析結(jié)果，從表5可見，受水區(qū)研究期平水年徑流量為261.9m3，基準期平水年徑流量為259.2m3。河北省1956—2016年多年平均徑流量最小，僅為2.6億m3，山東省和江蘇省多年平均徑流量均大于100億m3，安徽省和天津市多年平均徑流量均大于10億m3。與基準期對比，比較期江蘇省、安徽省受水區(qū)徑流量有較大增加，增長率分別達到18.5%、10.0%，但山東省、河北省和天津市受水區(qū)天然徑流量較基準期有所下降，3省市受水區(qū)合計減少10.1億m3，相對減少率分別為5.8%、4.0%和9.6%。綜合分析受水區(qū)降水與徑流的均值變化情況，以山東省和天津市為例，比較期山東省降水較基準期增加3.5%，但天然徑流下降5.8%；天津市降水較基準期下降5.2%，但天然徑流量下降9.6%，要高于降水的下降率。分析原因可能是人類活動對下墊面產(chǎn)生了一定的影響，導(dǎo)致了以上兩種情況的出現(xiàn)，說明人類活動在一定程度上影響了區(qū)域水資源量[22-23]。

表5 受水區(qū)徑流量水文頻率分析結(jié)果Table 5 Hydrological frequency analysis results of runoff in the catchment area

(a) 全部受水區(qū)

(b) 江蘇受水區(qū)

(c) 山東受水區(qū)

(d) 安徽受水區(qū)

(e) 河北受水區(qū)

(f) 天津受水區(qū)圖11 東線受水區(qū)1956—2016年天然徑流量變化Fig.11 Variation of natural runoff in the water-receiving area of Eastern Route from 1956 to 2016

b.M-K趨勢分析。表6為受水區(qū)年徑流量M-K趨勢分析結(jié)果，經(jīng)95%置信水平統(tǒng)計檢驗，總受水區(qū)以及江蘇省、安徽省和天津市受水區(qū)年徑流量表現(xiàn)出一定下降態(tài)勢但不顯著，而山東省、河北省受水區(qū)年徑流量下降態(tài)勢較其他受水區(qū)明顯，這在一定程度上說明兩省受水區(qū)的水資源量呈下降趨勢，反映了南水北調(diào)工程東線調(diào)水的必要性。

表6 受水區(qū)年徑流量M-K趨勢分析結(jié)果Table6 M-K trend analysis results of annual runoff in the catchment area

c.小波周期分析。圖12為1956—2016年受水區(qū)徑流量小波變換系數(shù)實部等值線，可見徑流變化主周期為20 a，次周期為12 a，主次周期內(nèi)徑流量經(jīng)歷了“枯-豐-枯-豐”的循環(huán)交替過程。總體來說，在4 a、12 a、20 a左右時間尺度的豐枯交替變化特征比較明顯。在20 a的周期尺度上，共經(jīng)歷了4次較為完整的豐枯交替過程，推測目前受水區(qū)正處于徑流偏少期。圖13為1956—2016年受水區(qū)徑流量小波變換方差過程線，可見，在4 a、12 a、20 a左右時間尺度的極值周期較明顯。最大峰值對應(yīng)著20 a左右的時間尺度，說明20 a左右的周期震蕩最強，為東線受水區(qū)天然年徑流量的第一主周期；12 a左右時間尺度對應(yīng)著第二峰值，為第二主周期；4 a左右的時間尺度為第三主周期。

圖12 1956—2016年受水區(qū)徑流量小波變換系數(shù)實部等值線Fig.12 Real part contour of wavelet transform coefficient of runoff in catchment area from 1956 to 2016

圖13 1956—2016年受水區(qū)徑流量小波變換方差過程線Fig.13 Wavelet transform variance hydrograph of runoff in catchment area from 1956 to 2016

4 結(jié) 論

a.東線調(diào)水區(qū)水資源條件較為平穩(wěn)，其統(tǒng)計特征無顯著變化。相對于基準期，比較期東線調(diào)水區(qū)降水量變化不明顯，長江大通站天然徑流量較基準期略有下降。以長江大通控制斷面徑流計算，根據(jù)《總體規(guī)劃》，東線工程一、二期共調(diào)水195.23億m3，調(diào)水量僅占大通站多年平均年徑流量的2%左右。因此，可以認為長江流域天然水資源條件可穩(wěn)定支撐東線工程建設(shè)運行。

b.東線受水區(qū)部分地區(qū)本地水資源嚴重不足的局面在2000年之后無明顯變化。相對于基準期，比較期東線受水區(qū)年降水量較基準期基本無變化，天然徑流量總體較基準期上升3.9%，但山東省、河北省和天津市受水區(qū)天然徑流量較基準期有所下降，3省市受水區(qū)合計減少10.1億m3，相對減少率為6.6%。