張嵐婷，王文圣，劉淺奎，鄭芳芳

(四川大學(xué) 水利水電學(xué)院，四川 成都 610065)

流域水資源變化量、水循環(huán)過程皆和氣候變化及人類活動(dòng)有著顯著關(guān)聯(lián)[1-3]。確定徑流量的變化特征，量化氣候因素和人類活動(dòng)對(duì)徑流變化的作用，既具理論意義又有實(shí)用價(jià)值。在徑流變化特征研究中，廣泛應(yīng)用了統(tǒng)計(jì)方法，并取得了大量研究成果。對(duì)于大渡河流域上的徑流變化特征，國內(nèi)學(xué)者已開展了一些工作，例如陳媛利用沙坪站1937—2008 年徑流資料，對(duì)大渡河流域徑流的變化特征進(jìn)行了定性識(shí)別和定量分析[4]。定量分析氣候變化和人類活動(dòng)對(duì)徑流影響的研究方法種類主要有水文模型模擬法類型和以氣候彈性系數(shù)法、降水-徑流雙累積曲線法、敏感性分析法等為代表的定量評(píng)估法類型[5-8]，前者具備較好的物理基礎(chǔ)，但參數(shù)存在一定的不確定性[9]；后者計(jì)算簡(jiǎn)單，但需要較長(zhǎng)的數(shù)據(jù)序列且序列中的噪聲會(huì)對(duì)評(píng)估結(jié)果造成干擾[10]。由王隨繼等[3]提出的累積量斜率變化率比較法，可以有效剔除噪聲，能方便地分離出氣候變化和人類活動(dòng)對(duì)徑流的影響，并被廣泛應(yīng)用于黃河流域、松花江流域和黑河流域。

本文收集1951—2012 年大渡河流域內(nèi)氣象站年降水量、年平均氣溫和銅街子站年徑流量資料，采用多種方法、多種途徑系統(tǒng)分析大渡河流域徑流量的變化特征，同時(shí)探討氣候變化和人類活動(dòng)對(duì)流域徑流量變化的影響程度。

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

1.1 研究區(qū)概況

大渡河是岷江的最大支流，干流長(zhǎng)1 062 km，流域面積77 400 km2(不含青衣江)。瀘定以上為上游，集水面積為58 943 km2，占全流域集水面積的76.2%；瀘定至銅街子為中游，集水面積為17 440 km2，占全流域集水面積的22.5%；銅街子以下為下游，集水面積為1 017 km2，占全流域集水面積的1.3%[11]。銅街子水文站為大渡河流域徑流的控制站，上、中、下游共選了14 個(gè)氣象站。水文站及氣象站分布見圖1。

分別收集了1951—2012 年銅街子站的年徑流量資料和14 個(gè)氣象站年降水量、年平均氣溫資料。14 個(gè)氣象站較均勻地分布在流域上，運(yùn)用算術(shù)平均法計(jì)算流域面降水量和流域平均氣溫。

1.2 研究方法

1.2.1 徑流變化特征分析方法 本文基于有序聚類法[12]、累積距平法及Pettitt 檢驗(yàn)[13]對(duì)銅街子站年徑流量序列進(jìn)行突變點(diǎn)識(shí)別，3 種方法互為驗(yàn)證，拐點(diǎn)可信度更高；采用Mann-Kendall 秩次相關(guān)檢驗(yàn)法[12]檢驗(yàn)銅街子站年徑流量的變化趨勢(shì)；利用Morlet 小波分析法[14]識(shí)別出年徑流量的多尺度時(shí)間變化和主要的周期成分。限于篇幅，這里不再贅述相關(guān)分析方法。

圖1 大渡河流域氣象站及水文站位置Fig. 1 Location of meteorological stations and hydrological stations of the Dadu River basin

1.2.2 徑流變化歸因分析方法 本文應(yīng)用王隨繼等提出的累積量斜率變化率比較法（SCRCQ）[3]定量分析大渡河流域氣候因素和人類活動(dòng)對(duì)徑流量變化的影響。累積量斜率變化率比較法的基本原理和步驟如下。

先確定突變年份，以突變年份為拐點(diǎn)，點(diǎn)繪拐點(diǎn)前、后時(shí)段不同累積量與年份的關(guān)系曲線，然后擬合拐點(diǎn)前、后時(shí)段相應(yīng)累積量與年份的線性回歸方程，突變前、后斜率分別計(jì)為Ka和Kb。

對(duì)于年徑流量，突變前、后的斜率為KRa和KRb，累積徑流量斜率變化率定義為：

式中：RSR為累積徑流量斜率變化率。

同理，對(duì)于年降水，

式中：RSP為累積降水量的斜率變化率；KPa和KPb分別為突變前、后的斜率。

對(duì)于年氣溫，

式中：RST為累積氣溫的斜率變化率；KTa和KTb分別為突變前、后的斜率。

降水和氣溫對(duì)徑流量的貢獻(xiàn)率CP、CT由斜率變化率計(jì)算得出，

氣候要素對(duì)徑流量變化的貢獻(xiàn)率（CC）為降水和氣溫因素貢獻(xiàn)率的總和，即：CC=CP+CT。

人類活動(dòng)變化和氣候要素共同影響徑流量變化，則人類活動(dòng)對(duì)徑流量的貢獻(xiàn)率（CH）為：CH=1?CP?CT。

2 結(jié)果與討論

2.1 年徑流量變化特征

2.1.1 年徑流量的趨勢(shì)性 根據(jù)年徑流量突變點(diǎn)識(shí)別中的突變年份把銅街子水文站年徑流序列分為3 段：1951—1968 年、1969—1987 年和1988—2012 年。各時(shí)段的徑流量均值顯著不同且呈減少趨勢(shì)。采用Mann-Kendall 趨勢(shì)檢驗(yàn)法分別對(duì)銅街子站不同時(shí)段的年徑流量序列變化趨勢(shì)進(jìn)行識(shí)別，結(jié)果見表1。

由表1 可知，1951—2012 年期間，銅街子站年徑流量序列表現(xiàn)為不顯著的減少趨勢(shì)；1951—1968 年期間和1988—2012 年期間，年徑流量序列存在不顯著的減少趨勢(shì)；1969—1987 年期間，年徑流序列表現(xiàn)出不顯著的增加趨勢(shì)。

表1 銅街子水文站年徑流趨勢(shì)Mann-Kendall 檢驗(yàn)結(jié)果Tab. 1 Mann-Kendall test results of annual runoff trend at Tongjiezi Station

2.1.2 年徑流量的周期性 采用小波分析法識(shí)別年徑流量的周期成分，圖2 繪制了年徑流量序列小波方差，圖3 給出了年徑流量序列小波變換系數(shù)實(shí)部時(shí)頻分布。由圖3 能夠看到年徑流量演化過程中明顯存在著多時(shí)間尺度特征，其中6～8 年、10～12 年、15～20 年、22～30 年尺度的周期震蕩尤為突出；由圖2 可知，年徑流量的主要周期成分為7、11、17、28 年，其中28 年的時(shí)間尺度與最大峰值相對(duì)應(yīng)，為流域年徑流量變化的第一主周期，即28 年左右的周期震蕩最強(qiáng)。

圖2 銅街子站年徑流量序列小波方差Fig. 2 Wavelet variance of annual runoff series at Tongjiezi station

圖3 銅街子站年徑流序列小波變換系數(shù)實(shí)部時(shí)頻分布Fig. 3 Real part change process of wavelet analysis for annual runoff series at Tongjiezi Station

2.1.3 年徑流量的突變性 本文采用有序聚類法、累積距平法及Pettitt 檢驗(yàn)等3 種方法，對(duì)銅街子水文站1951—2012 年期間年徑流量序列的突變點(diǎn)進(jìn)行綜合判別。

有序聚類法識(shí)別成果及累積距平法診斷成果見圖4，識(shí)別出的年徑流量序列突變年份均為1968 年。由Pettitt 檢驗(yàn)（給定顯著水平α=0.05）表明年徑流量序列突變亦發(fā)生在1968 年。3 種檢驗(yàn)方法可以推斷出大渡河流域年徑流量在1968 年存在顯著性突變。

圖4 1951—2012 年徑流量序列突變點(diǎn)識(shí)別Fig. 4 Diagnosis of change point based on annual runoff from 1951 to 2012

同樣采用3 種檢驗(yàn)方法對(duì)銅街子水文站1969—2012 年期間年徑流量序列進(jìn)行突變點(diǎn)識(shí)別。由圖5(a)的總離差平方和（Sn）變化量可知，有序聚類法無法尋求突變點(diǎn)。從圖5(b)可以發(fā)現(xiàn)，該時(shí)期年徑流量累積距平表現(xiàn)出先減少后增加的變化趨勢(shì)，其突變點(diǎn)為1987 年。由Pettitt 檢驗(yàn)亦得出1969—2012 年徑流量序列突變年份為1987 年，故可推斷大渡河流域年徑流量在1987 年亦發(fā)生顯著性突變。

圖5 1969—2012 年徑流量序列突變點(diǎn)識(shí)別Fig. 5 Diagnosis of change point based on annual runoff from 1969 to 2012

2.2 年徑流量變化的歸因分析

年徑流量變化是諸多因素共同作用的結(jié)果，大渡河流域年徑流量變化主要受氣候因素和人類活動(dòng)影響。氣候變化主要通過水、熱兩方面來呈現(xiàn)影響，因此將氣候因子歸結(jié)為水因素和熱因素，分別以降水和氣溫表征。

2.2.1 基準(zhǔn)期和影響期下的年徑流量、年降水量和年平均氣溫過程 以大渡河流域年徑流量突變年份1968 年和1987 年為轉(zhuǎn)折點(diǎn)將研究時(shí)段分為3 個(gè)時(shí)期：A 時(shí)段為1951—1968 年，B 時(shí)段為1969—1987 年，C 時(shí)段為1988—2012 年。大渡河流域年徑流量在A 時(shí)段主要受氣候變化影響，人類活動(dòng)影響偏弱；B、C 時(shí)段人類活動(dòng)增加，徑流量受人類活動(dòng)和氣候變化共同作用。因此，將A 時(shí)段設(shè)為基準(zhǔn)期，B、C 時(shí)段為影響期。大渡河流域在基準(zhǔn)期及影響期年徑流量過程、年降水量過程及年平均氣溫過程見圖6?？梢?，大渡河流域年徑流量多年來總體呈不顯著下降趨勢(shì)，降水總體呈微弱下降趨勢(shì)，而氣溫總體呈上升趨勢(shì)。

圖6 大渡河流域年徑流、年降水、年平均氣溫過程及其均值變化Fig. 6 Annual series and its average for runoff, precipitation and average temperature in Dadu River basin

2.2.2 不同時(shí)期累積年徑流、降水量、氣溫與年份之間的相關(guān)關(guān)系 分別繪制累積年徑流量與年份、累積降水量與年份和累積氣溫與年份在A、B、C 時(shí)段的關(guān)系曲線，并擬合線性關(guān)系，結(jié)果見圖7。

以上各擬合關(guān)系的相關(guān)系數(shù)均大于0.999，表明年份與累積徑流、累積降水量和累積氣溫的線性關(guān)系很好。

2.2.3 氣候變化和人類活動(dòng)對(duì)徑流量的貢獻(xiàn)率分析 表2 給出了大渡河流域累積年徑流量、累積年降水量和累積年氣溫關(guān)系的斜率及其變化情況。

表2 不同時(shí)期年份與各累積量關(guān)系的斜率及其變化Tab. 2 Slope and slope change of relation between year and accumulated quantity in different periods

由表2 可知：

(1) 與基準(zhǔn)期A 相比，影響期B 累積年徑流量與年份線性關(guān)系的斜率減小了6.16 萬m3/a，斜率變化率為?11.13%；與基準(zhǔn)期A 相比，影響期C 的斜率減小了2.19 萬m3/a，斜率變化率為?3.96%。

(2) 與基準(zhǔn)期A 相比，影響期B 累積年降水量與年份線性關(guān)系的斜率減小了36.32 mm/a，變化率為?4.40%；與基準(zhǔn)期A 相比，影響期C 的斜率小了23.88 mm/a，變化率為?2.89%。斜率變化率為負(fù)，表明降水序列有下降趨勢(shì)，與圖6(b)是一致的。

(3) 與基準(zhǔn)期A 相比，累積年氣溫與年份線性關(guān)系的斜率增大了0.029 7 ℃/a，變化率為0.37%；與基準(zhǔn)期A 相比，影響期C 的斜率增大了0.124 7 ℃/a，變化率為1.56%。斜率變化率均為正，且C 時(shí)段比B 時(shí)段變化率更大，表明氣溫序列呈現(xiàn)持續(xù)上升趨勢(shì)，與圖6(c)是一致的。

對(duì)大渡河流域而言，降水量增加，徑流量亦增加，年徑流量與年降水量呈正相關(guān)；氣溫升高，徑流量減少，年徑流量與年平均氣溫呈負(fù)相關(guān)。

由計(jì)算結(jié)果可知：(1) 與基準(zhǔn)期A 相比，影響期B 的氣候變化對(duì)年徑流量減少的貢獻(xiàn)率為36.22%，人類活動(dòng)對(duì)年徑流量減少的貢獻(xiàn)率為63.78%。(2) 與基準(zhǔn)期A 相比，影響期C 的氣候變化對(duì)年徑流量減少的貢獻(xiàn)率為33.68%，人類活動(dòng)對(duì)年徑流量減少的貢獻(xiàn)率為66.32%。

圖7 累積量與年份的關(guān)系曲線Fig. 7 Relationship curves of cumulative quantity and year

2.3 討 論

降水量對(duì)年徑流量的影響在C 時(shí)段遠(yuǎn)高于B 時(shí)段，平均氣溫對(duì)年徑流量的影響同樣在C 時(shí)段遠(yuǎn)高于B 時(shí)段。綜合考慮降水量和平均氣溫的共同影響，氣候變化對(duì)徑流量變化的貢獻(xiàn)率在B、C 時(shí)段分別為36.22%、33.68%，氣候變化的貢獻(xiàn)率在C 時(shí)段略低于B 時(shí)段，而人類活動(dòng)的貢獻(xiàn)率反而增大。

平均氣溫升高導(dǎo)致徑流量減少，一是江河湖泊水面蒸發(fā)增大致使徑流量減少；二是流域內(nèi)植物蒸騰作用及水分截留過程受到影響，隨之改變徑流量。

人類活動(dòng)是大渡河流域徑流量變化的主導(dǎo)因素。1970 年代以來，社會(huì)經(jīng)濟(jì)狀況的轉(zhuǎn)好、人口數(shù)量的增多使工農(nóng)業(yè)用水量和生活耗水量加大，從而直接導(dǎo)致徑流量減少；此外，大渡河流域上豐富的水力資源得到大量開發(fā)，諸多水利水電工程的修建與運(yùn)行皆對(duì)徑流量的自然變化過程產(chǎn)生影響；同時(shí)，工農(nóng)業(yè)開發(fā)及房地產(chǎn)業(yè)等的快速發(fā)展使得土壤利用類型和植被覆蓋率也發(fā)生了顯著變化，進(jìn)而使地表入滲、匯流及補(bǔ)給等水循環(huán)過程改變，間接影響到徑流量。因此人類活動(dòng)對(duì)年徑流量減少的貢獻(xiàn)率很大，且影響逐步加深。

累積量斜率變化率比較法在徑流量歸因分析的應(yīng)用，均是處于將氣候變化和人類活動(dòng)視作兩大完全獨(dú)立因素的情景中，本文將降水和氣溫也當(dāng)成互不干擾的兩個(gè)氣候因子來分別計(jì)算貢獻(xiàn)率。實(shí)際上，氣候變化與人類活動(dòng)相互影響且關(guān)系密切，降水和氣溫之間也存在著一定的相關(guān)關(guān)系。尋求新的研究方法，對(duì)氣候變化和人類活動(dòng)的影響進(jìn)行更科學(xué)地分離和判定，是歸因分析研究需進(jìn)一步開展的內(nèi)容。

3 結(jié) 語

利用大渡河流域徑流量及氣象資料，通過多種方法分析年徑流量變化特征，揭示了氣候變化與人類活動(dòng)對(duì)流域年徑流量的影響，主要結(jié)論如下：

（1）大渡河流域年徑流量在1951—2012 年、1951—1968 年和1988—2012 年期間呈不顯著的減少趨勢(shì)，1969—1987 年期間呈不顯著的增加趨勢(shì)；年徑流量變化普遍存在6～8 年、10～12 年、15～20 年、22～30 年尺度的周期震蕩，主要周期成分為7、11、17、28 年；流域年徑流量在1968 年、1987 年發(fā)生突變。

（2）以人類活動(dòng)輕微的A 時(shí)段（1951—1968 年）為基準(zhǔn)期，影響期B（1969—1987 年）氣候變化對(duì)年徑流量減少的貢獻(xiàn)率為36.22%，人類活動(dòng)對(duì)徑流減少的貢獻(xiàn)率達(dá)63.78%；影響期C（1988—2012 年）氣候變化對(duì)年徑流量減少的貢獻(xiàn)率為33.68%，人類活動(dòng)的貢獻(xiàn)率達(dá)66.32%。人類活動(dòng)是大渡河流域年徑流量減少的主導(dǎo)因素，氣候變化是次要因素。

（3）揭示了大渡河流域近60 年來年徑流量的變化趨勢(shì)及主要影響因素的貢獻(xiàn)率，為研究區(qū)徑流量變化歸因分析、水資源合理開發(fā)利用、生態(tài)環(huán)境修復(fù)與管理提供了參考。