西江干流徑流年際及年內(nèi)變化趨勢(shì)分析

摘要：為了分析西江干流徑流年際變化趨勢(shì)及年內(nèi)分配規(guī)律，根據(jù)天峨、遷江、梧州三站60多年實(shí)測(cè)徑流資料，采用MannKendall趨勢(shì)檢驗(yàn)法、R/S法、有序聚類法、滑動(dòng)T檢驗(yàn)法等多方法分別對(duì)徑流年際變化的趨勢(shì)性和突變性進(jìn)行分析；同時(shí)，對(duì)表征徑流年內(nèi)分配的各指標(biāo)進(jìn)行趨勢(shì)性檢驗(yàn)和未來(lái)變化趨勢(shì)預(yù)測(cè)。結(jié)果表明：三站的MK統(tǒng)計(jì)量Z值分別為247、-238和-147，Hurst指數(shù)分別為072、083和080，均大于05；天峨站和遷江站的突變點(diǎn)均為2002年，梧州站的突變點(diǎn)為1954年和2002年；西江徑流60多年來(lái)顯著減少，且具較強(qiáng)的持續(xù)性；徑流年內(nèi)分配不均，多集中在4月-9月，占年徑流的75%左右；2002年后徑流年內(nèi)分配顯著性地趨于均勻，集中程度和變幅亦顯著性地變小，集中期前移12日左右。

關(guān)鍵詞：西江；徑流；MK趨勢(shì)檢驗(yàn)法；R/S法；年際變化趨勢(shì)；年內(nèi)變化趨勢(shì)

中圖分類號(hào)：TV121文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：

16721683（2018）04007408

Study on interannual and annual variation characteristics of runoff in the main stream of Xijiang River

CHEN Lihua1，2，3，LIU Weifu1，LENG Gang1，WANG Yan1

（

1.College of Civil engineering and Architecture，Guangxi University，Nanning 530004，China；

2.Key Laboratory of Disaster Prevention and Structural Safety of Ministry of Education，Guangxi University，Nanning 530004，China；

3.Guangxi Key Laboratory of Disaster Prevention and Engineering Safety，Guangxi University，Nanning 530004，China）

Abstract：

In order to study the interannual variation trend and annual distribution of runoff in Xijiang River，we used MannKendall trend test，R/S method，ordered cluster method，and sliding T test to analyze the variation trend and abrupt changes of interannual runoff，and we conducted trend tests and prediction of the future change trend for each index of annual runoff distribution.The analysis was based on the measured runoff data of the main stream of Xijiang River for over 60 years.The results showed that MK statistic values in the three stations were -247，-238 and -147 respectively.The Hurst indexes of the three stations were 072，083 and 080 respectively，all greater than 05.The abrupt change points at Tian'e and Qianjiang stations were both in 2002，while the abrupt change points at Wuzhou station were in 1954 and 2002.Xijiang River runoff has significantly declined over the past 60 years with a strong continuity.The annual distribution of runoff was uneven.The runoff concentrated during April to September，which accounted for about 75% of the annual runoff.After 2002，the annual distribution of runoff showed a remarkable tendency towards uniformity； the degree of concentration and variation amplitude both significantly decreased，and the concentration period moved about 12 days forward.

Key words：

Xijiang River； runoff； MannKendall trend test； R/S method； interannual variation trend； annual variation trend

河川徑流變化通常包含徑流的量值變化和過程的結(jié)構(gòu)變化[13]，前者包括徑流總量、流量等數(shù)值上的變化，而后者側(cè)重于年內(nèi)不同時(shí)段徑流組分分配情況。河川徑流過程是氣候條件和流域下墊面綜合作用的產(chǎn)物[4]，氣候變化及人類活動(dòng)引起全球水文循環(huán)變化加劇，導(dǎo)致各種極端水文事件[56]頻繁發(fā)生。近年來(lái)珠江流域[7 8]由氣候變化引起的極端水文事件時(shí)有發(fā)生，西江作為珠江流域的主要水系干流，具有流域面積大、水文地質(zhì)、地貌和覆蓋條件復(fù)雜等特點(diǎn)，在氣候變化及人類活動(dòng)雙重影響下，其河川徑流發(fā)生不同程度的改變[911]，進(jìn)而影響河床及水生態(tài)系統(tǒng)演變，增大流域水資源開發(fā)利用難度[1213]，其中下游珠江三角洲河口地區(qū)供水安全問題亟待解決。隨著珠江三角洲地區(qū)經(jīng)濟(jì)和人口快速增長(zhǎng)，區(qū)域水資源短缺、水生態(tài)惡化、流域水資源調(diào)配能力低等問題日趨突出，特別是冬春季咸潮入侵對(duì)珠江三角洲地區(qū)1500萬(wàn)居民飲水安全構(gòu)成極大威脅，而抑咸流量的80%來(lái)源于西江。因此，西江流域徑流量變化與地表水文過程的研究備受關(guān)注，在水文要素量值變化方面已有一些研究成果[1415]，但缺少對(duì)量值突變前后對(duì)比研究及對(duì)地表水文過程結(jié)構(gòu)特征的研究。

近年來(lái)，徑流趨勢(shì)變化及突變分析成為徑流變化規(guī)律分析的研究熱點(diǎn)，目前對(duì)徑流序列趨勢(shì)分析的方法主要有MK趨勢(shì)檢驗(yàn)法、小波分析法、R/S分析法等。MK趨勢(shì)檢驗(yàn)法[16]和R/S分析法[17]相較于其他方法，具有原理簡(jiǎn)單、可操作性強(qiáng)、表征效果好等特點(diǎn)，故采用上述兩種方法進(jìn)行徑流趨勢(shì)變化分析。對(duì)徑流序列突變分析的方法有滑動(dòng)平均法、有序聚類法、滑動(dòng)T檢驗(yàn)法、MannKendall法等。基于有序聚類法和滑動(dòng)T檢驗(yàn)法計(jì)算簡(jiǎn)單、資料要求低、結(jié)果直觀等特點(diǎn)，在渭河流域[18]、洞庭湖[19]、新村流域[20]等地區(qū)得到較好應(yīng)用，故采用上述兩種方法對(duì)徑流進(jìn)行綜合突變?cè)\斷。對(duì)徑流年內(nèi)分配的研究方法，燕華云等[21]將表示降水量年內(nèi)分配的向量法引入徑流集中度和集中期并結(jié)合不均勻系數(shù)、變化幅度等定量指標(biāo)進(jìn)行表示，以突變點(diǎn)劃分徑流階段，對(duì)表征徑流年內(nèi)分配各指標(biāo)進(jìn)行了趨勢(shì)性檢驗(yàn)和未來(lái)變化趨勢(shì)預(yù)測(cè)，以期為西江流域水資源合理開發(fā)利用及防洪減災(zāi)提供參考。

1流域概況與分析方法

1.1西江流域概況與研究數(shù)據(jù)

西江是珠江流域的主干流，自上而下依次為南盤江、紅水河、黔江、潯江和西江五個(gè)河段，全長(zhǎng)2 075 km，流域面積達(dá)3531萬(wàn)km2。根據(jù)西江干流水文站點(diǎn)分布，河段徑流代表性，及水文資料完整性、準(zhǔn)確性和觀測(cè)時(shí)間的持續(xù)性，選取西江干流天峨、遷江、梧州三站60多年的月平均流量資料，分析西江干流年際及年內(nèi)變化趨勢(shì)。三站分別位于天峨縣六排鎮(zhèn)、來(lái)賓市興賓區(qū)遷江鎮(zhèn)和梧州市，多年平均徑流量分別為1 8004 m3/s、2 3769 m3/s、7 5672 m3/s[14]。西江干流水文站點(diǎn)及資料情況見表1。

1.2徑流分析方法

MannKendall法是由Mann[22]1945年首次提出，而后Kendall[23]1975年進(jìn)行改進(jìn)的一種被廣泛應(yīng)用于氣象與水文時(shí)間序列的趨勢(shì)分析[24]的非參數(shù)檢驗(yàn)方法。當(dāng)MK統(tǒng)計(jì)值Z的絕對(duì)值大于等于128、164、232 時(shí)，表示對(duì)應(yīng)的徑流序列變化趨勢(shì)分別通過了信度為90%、95%和99%的顯著性檢驗(yàn)。

R/S法是由Hurst在分析尼羅河水文資料時(shí)提出的一種分析方法[25]，可通過R/S法估算Hurst指數(shù)（0

以有序聚類法來(lái)推求最優(yōu)化分割點(diǎn)，使同類之間離差平方和較小，而類與類之間的離差平方和較大。該法原理[26]清晰簡(jiǎn)單，但對(duì)于多個(gè)突變點(diǎn)需要進(jìn)一步分割序列進(jìn)行檢驗(yàn)。

1.3年內(nèi)分配特征指標(biāo)

采用徑流不均勻系數(shù)C[WTB1X]v[WTBX]、完全調(diào)節(jié)系數(shù)Cr表征徑流年內(nèi)分配均勻程度；通過最大月平均流量Rmax、最小月平均流量Rmin與年平均流量之比，分別稱為極大比Cmax、極小比Cmin及極值比Cm，表征徑流相對(duì)變化幅度；引入集中度Cn和集中期D表征徑流年內(nèi)分配集中的程度和集中的重心出現(xiàn)月份。集中度就是將各月流量按一定角度以向量方式累加，其各分量之和的合成向量的模與年徑流量的比值；集中期則是合成向量的方向，以角度表示，可以換算到月份，本研究從1月至12月每月的方位角θ分別為0°，30°，60°，…，360°，并將每個(gè)月平均流量分解到x和y兩個(gè)方向上。各指標(biāo)計(jì)算公式見文獻(xiàn)[21]。

[BT2-*6]2徑流年際變化特征

2.1徑流年際變化趨勢(shì)分析

采用MK法及線性回歸對(duì)西江干流天峨、遷江、梧州三個(gè)水文站年平均流量進(jìn)行趨勢(shì)分析，結(jié)果見圖1及表2。

由表2可知，20世紀(jì)60年代至今，三站流量序列ZC均小于0，表明西江干流年平均流量總體呈下降趨勢(shì)，其中天峨、遷江站流量序列|ZC|>ZCα/2（顯著性α=005的臨界值為164），分別以-626 m3/（10a·s）和-637 m3/（10a·s）的趨勢(shì)顯著下降；梧州站的|ZC|=147

采用R/S法估算三站年平均流量Hurst指數(shù)（見表3），可知三站徑流序列R/S擬合曲線擬合度R2值均超過098，表明擬合效果良好。三站Hurst指數(shù)分別為072、083和080，均明顯大于05，表明未來(lái)年平均流量與過去60多年具相同的變化趨勢(shì)。結(jié)合三站年平均流量ZC統(tǒng)計(jì)值和Hurst指數(shù)可預(yù)測(cè)：未來(lái)西江干流年平均流量仍將呈持續(xù)減少趨勢(shì)，其中中上游流域流量呈顯著下降趨勢(shì)，而下游流量減少趨勢(shì)不顯著。

2.2徑流年際突變特征

采用累積距平和有聚序類法對(duì)天峨、遷江、梧州三站年平均流量突變特征進(jìn)行綜合檢驗(yàn)（見圖2）。圖2（a）至圖2（c）為累計(jì)距平曲線，三站均存在階段變化，階段上升代表豐水年組，階段下降代表枯水年組。三站均在2002年累計(jì)距平值達(dá)到極大值，為豐轉(zhuǎn)枯節(jié)點(diǎn)。圖2（d）至圖2（f））為有序聚類法檢驗(yàn)結(jié)果，天峨站和遷江站的突變點(diǎn)均為2002年。梧州站在檢驗(yàn)出1954年為突變點(diǎn)后，考慮其后序列較長(zhǎng)，徑流變化強(qiáng)烈，可能存在第二個(gè)突變點(diǎn)，對(duì)1955-2015年徑流序列進(jìn)行檢測(cè)，發(fā)現(xiàn)2002年為潛在突變點(diǎn)。

因累積距平和有序聚類法突變點(diǎn)檢測(cè)結(jié)果存在不一致的缺點(diǎn)，因此輔以滑動(dòng)T檢驗(yàn)法進(jìn)一步對(duì)可能存在的突變點(diǎn)進(jìn)行顯著性檢驗(yàn)，結(jié)果見圖3。

經(jīng)滑動(dòng)T檢驗(yàn)（見圖3）發(fā)現(xiàn)三站年徑流均在2002年發(fā)生突變，且均超過005置信水平（對(duì)應(yīng)臨界值均為20），為顯著突變點(diǎn)，但原因有所差異。文獻(xiàn)[9]采用可變范圍（RVA）法評(píng)估龍灘水庫(kù)建設(shè)與運(yùn)行下對(duì)天峨站徑流的影響發(fā)現(xiàn)水庫(kù)修建后對(duì)逐月徑流量有顯著影響而對(duì)年徑流量影響有限；文獻(xiàn)[11]通過比較分析2006-2008年龍灘水電站蓄水量與遷江站年徑流量的響應(yīng)關(guān)系發(fā)現(xiàn)其調(diào)蓄作用對(duì)

下游年徑流量變化同樣影響有限，因此2000年末投產(chǎn)運(yùn)行的天生橋一級(jí)電站蓄水調(diào)節(jié)影響并非造成中上游年徑流突變的主要原因。天峨站和遷江站來(lái)水單一，徑流一致性較好，綜合上述因素可知中上游徑流突變主要受氣候變化影響。下游梧州站年徑流突變?cè)蜉^為復(fù)雜，一方面，朱穎潔等[10]基于神經(jīng)網(wǎng)絡(luò)年徑流量模擬的氣候變化和人類活動(dòng)的定量分析認(rèn)為氣候變化是梧州站年徑流量變化的主要原因；另

一方面，遷江至梧州區(qū)間流域廣闊，降水差異較大，且三大支流匯入致使水文情勢(shì)復(fù)雜，加上2002年郁江百色水利樞紐大江截流的有限影響綜合導(dǎo)致了下游年徑流突變。因此，氣候變化是西江干流年徑流突變的主要原因，人類活動(dòng)如水利工程建設(shè)的影響有限。因三站年徑流均在2002年發(fā)生突變，以2002年為分界點(diǎn)將研究數(shù)據(jù)分為兩段，進(jìn)一步探討和分析西江干流年際及年內(nèi)變化特征。

2.3徑流年際變化特征值

以突變節(jié)點(diǎn)（2002年）為分界點(diǎn)將序列劃分為兩段，統(tǒng)計(jì)分析三站徑流突變前后特征值變化（見表4）。由表可知：（1）三站徑流年際極值比與C[WTB1X]v[WTBZ]值均較小，表明西江干流徑流年際變幅平緩；（2）2002年前后兩個(gè)時(shí)段，三站徑流的時(shí)段特征值存在較大差異。2002后三站徑流的時(shí)段均值、距平值、時(shí)段最大年徑流、極值比和C[WTB1X]v[WTBZ]值均低于前一時(shí)段的對(duì)應(yīng)值，即干流徑流年際變化趨于均勻，變幅減小，說明中上游梯級(jí)水庫(kù)發(fā)揮了一定的調(diào)控作用；對(duì)于時(shí)段最小年徑流，遷江站和梧州站增大，天峨站減小，一方面為2005年1月起通過天生橋一級(jí)、龍灘、巖灘、百色及長(zhǎng)洲水利樞紐進(jìn)行的枯季壓咸補(bǔ)淡應(yīng)急調(diào)水；另一方面流域干旱時(shí)有發(fā)生，如水文年2004-2005年和2009-2010年均為特枯年、2006-2007年和2007-2008年為偏枯年，且龍灘上游的降水量少于受多個(gè)暴雨中心影響的下游降水量[27]，這就造成了天峨站最小年徑流較少，而遷江站和梧州站增大。

3徑流年內(nèi)變化特征

3.1徑流年內(nèi)分配及變化

隨著西江干流梯級(jí)水電站相繼建成運(yùn)行，其徑流年內(nèi)分配波動(dòng)明顯增大[9，11]，有必要針對(duì)變化環(huán)境下的徑流年內(nèi)分配特征進(jìn)行詳細(xì)分析。由表5可知，上游天峨站夏季（6月-8月）占年徑流量的512%，冬季（12月-2月）僅占96%，連續(xù)最大4個(gè)月為6月-9月，占全年徑流量的635%；中游遷江站夏季占年徑流量的541%，冬季僅占82%，連續(xù)最大4個(gè)月同樣為6月-9月，占全年徑流量的658%，兩站各季的徑流占比接近，具較高的一致性；而下游梧州站夏季占年徑流量的505%，冬季僅占82%，春季增長(zhǎng)明顯，達(dá)到201%，可見西江下游冬春交替之際補(bǔ)水形勢(shì)依然嚴(yán)峻。

基于各階段月平均徑流計(jì)算得表征三站年內(nèi)徑流的均勻度、變化幅度及集中程度特征指標(biāo)（見表6），可知：（1）三站C[WTB1X]v[WTBX] 、Cr值均較小，表明西江干流徑流年內(nèi)分配較均勻；（2）三站Cm、Cmax、Cmin值較接近，梧州站略小，表明西江下游年內(nèi)月徑流變化幅度較??；（3）就集中期D而言，天峨站徑流年內(nèi)集中期為7月下旬至8月上旬，遷江站為7月中旬至8月初，而梧州站主要集中在7月，基本與實(shí)測(cè)徑流系列的平均月徑流最大值出現(xiàn)時(shí)間相一致，表明集中度和集中期在量化徑流年內(nèi)分配方面很具代表性。（4）西江干流徑流的年內(nèi)分配在2002年前后存在較大差異，2002年后的C[WTB1X]v[WTBX] 、Cr、Cn、D值明顯變小，Cm、Cmax、Cmin值向中值靠攏，表明西江干流徑流年內(nèi)分配趨于均勻，徑流集中程度和變化幅度均減小，集中期前移約12日。

3.2變化趨勢(shì)分析

采用MK法對(duì)三站各階段年內(nèi)月徑流特征指標(biāo)的趨勢(shì)性進(jìn)行檢驗(yàn)（見表7）?？芍?）整個(gè)時(shí)域內(nèi)，天峨站和遷江站年內(nèi)月徑流變化趨勢(shì)一致性較高，兩站C[WTB1X]v[WTBX]、Cr、Cmax、Cm、Cn、D序列的變化趨勢(shì)均[JP+1]為顯著下降，Cmin序列為顯著上升，表明西江干流中上游徑流年內(nèi)分配顯著性地趨于均勻，集中程度[CM（22]與變幅也顯著性地變??；（2）下游梧州站年內(nèi)月徑流

3.3R/S分析

利用R/S法分析西江干流年內(nèi)月徑流特征指標(biāo)序列的Hurst指數(shù)（見表8）。結(jié)果表明：三站年內(nèi)月徑流的Hurst指數(shù)均大于05，表明西江干流徑流年內(nèi)分配變化趨勢(shì)與過去60多年徑流具有相同的變化趨勢(shì)。具體地，從C[WTB1X]v[WTBX]、Cr的Hurst值來(lái)看，徑流年內(nèi)分配趨于均勻的持續(xù)性：天峨>梧州>遷江；Cmax、Cmin、Cm的Hurst值表明，年內(nèi)月徑流變化幅度變小趨勢(shì)的持續(xù)性同樣為：天峨>梧州>遷江；表明西江干流徑流年內(nèi)分配趨于均勻和年內(nèi)變幅變小的持續(xù)性自上游至下游減弱，但梧州站因干支流徑流來(lái)水的復(fù)雜性，呈現(xiàn)較強(qiáng)的持續(xù)性；由Cn、D的Hurst值可知，三站徑流集中度變小和集中期前移的持續(xù)性為：遷江>天峨>梧州 ，分析認(rèn)為，天峨站經(jīng)上游水庫(kù)徑流調(diào)節(jié)年內(nèi)徑流集中度變小和集中期遷移具較強(qiáng)的持續(xù)性，遷江站則經(jīng)梯級(jí)水庫(kù)徑流調(diào)節(jié)其持續(xù)性最強(qiáng)，而梧州站雖然處于下游，同樣經(jīng)過梯級(jí)電站調(diào)節(jié)，但因來(lái)水多元化，不確定性較大，其持續(xù)性有所減弱。

4結(jié)論

（1）根據(jù)西江三站60多年徑流資料分析可知，西江干流中上游天峨站與遷江站年平均流量均呈顯著下降趨勢(shì)，且均在2002年發(fā)生突變；下游梧州站年平均流量為不顯著下降趨勢(shì)，突變點(diǎn)為1954年和2002年。說明2002年后西江流域徑流受氣候變化和人類活動(dòng)干擾加劇，還值得深入研究氣候變化和人類活動(dòng)對(duì)西江徑流的影響。

（2）西江干流徑流多集中在4月-9月，冬春季占比低。2002年后西江干流中上游徑流年內(nèi)分配強(qiáng)顯著性地趨于均勻，集中程度和變幅強(qiáng)顯著性地變??；下游徑流年內(nèi)分配趨于均勻，但變化不顯著，說明遷江站至梧州站區(qū)間來(lái)水可控性較差。

（3）天峨站徑流集中期為7下旬至8月上旬，遷江站為7月中旬至8月初，而梧州站主要集中在7月，表明干流中上游來(lái)水對(duì)梧州流量影響有限，若深入分析不同水平年梧州流量地區(qū)組成將對(duì)流域防洪、生態(tài)調(diào)度研究具有重要意義。

參考文獻(xiàn)（References）：

[1]HANNAH D M，SMITH B P G，GURNELL A M，et al.An approach to hydrograph classification[J].Hydrological Processes，2000，14（2）：317338.DOI：10.1002/（SICI）10991085（20000215）14：2<317：：AIDHYP929>3.0.CO；2T.

[2]HARRIS N M，GURNELL A M，HANNAH D M，et al.Classification of river regimes：a context for hydroecology[J].Hydrological Processes，2000，14（1617）：28312848.DOI：10.1002/10991085（200011/12）14：16/17<2831：：AIDHYP122>3.3.CO；2F.

[3]鄭紅星，劉昌明.黃河源區(qū)徑流年內(nèi)分配變化規(guī)律分析[J].地理科學(xué)進(jìn)展，2003，22（6）：585590.（ZHENG H X，LIU C M.Changes of annual runoff distribution in the headwater of the Yellow River Basin[J].Progress in Geography，2003，22（6）：585590.（in Chinese）） DOI：10.11820/dlkxjz.2003.06.006.

[4]王云琦，齊實(shí)，孫閣，等.氣候與土地利用變化對(duì)流域水資源的影響——以美國(guó)北卡羅萊納州Trent流域?yàn)槔齕J].水科學(xué)進(jìn)展，2011，22（1）：5158.（WANG Y Q，QI S，SUN G，et al.Impacts of climate and landuse change on water resources in a watershed：a case study on the Trent River Basin in North Carolina，USA[J].Advances in Water Science，2011，22（1）：5158.（in Chinese）） DOI：10.14042/j.cnki.32.1309.2011.01.006.

[5]APURV T，MEHROTRA R，SHARMA A，et al.Impact of climate change on floods in the Brahmaputra basin using CMIP5 decadal predictions[J].Journal of Hydrology，2015，527：281291.

[6]TRENBERTH K E.Changes in precipitation with climate change[J].Climate Research，2011，47（1）：123138.DOI：10.3354/cr00953.issn：0936577x.

[7]DUAN L M，ZHANG J Y，LI W，et al.Multivariate properties of extreme precipitation events in the Pearl River Basin，China：Magnitude，frequency，timing，and related causes[J].Hydrological Processes，2017，31（21）：36623671.DOI：10.1002/hyp.11286.

[8]ZHENG Y H，HE Y H，CHEN X H.Spatiotemporal pattern of precipitation concentration and its possible causes in the Pearl River Basin，China[J].Journal of Cleaner Production，2017，161（10）：10201031.DOI：10.1016/j.jclepro.2017.06.156.

[9]陳立華，賴河濤，葉明.龍灘水庫(kù)工程對(duì)紅水河天峨站水文情勢(shì)的影響分析[J].廣西大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，42（1）：394402.（CHEN L H，LAI H T，YE M.Influence of Longtan Reservior Project on hydrological regine of Hongshui River[J].Journal of Guangxi University（Nat Sci Ed），2017，42（1）：394402.（in Chinese）） DOI：10.13624 /j.cnki.issn.10017445.2017.0394.

[10]朱穎潔，郭純青，黃夏坤.氣候變化和人類活動(dòng)影響下西江梧州站降水徑流演變研究[J].水文，2010，30（3）：5055.（ZHU Y J，GUO C Q，HUANG X K.Change of precipitation and runoff at Wuzhou Station on Xijiang River under influence of climate changes and human activities[J].Journal of China Hydrology，2010，30（3）：5055.（in Chinese）） DOI：10.3969/j.issn.10000852.2010.03.011.

[11]陳立華，張利娜，劉為福.水電站運(yùn)行影響下紅水河下游水文情勢(shì)分析[J].廣西大學(xué)學(xué)報(bào)（自然科學(xué)版），2018（1）：397405.（CHEN L H，ZHANG L N，LIU W F.Hydrologic analysis of the lower reaches of Hongshui river under operation of hydropower stations[J].Journal of Guangxi University（Nat Sci Ed），2018（1）：397405.（in Chinese））.

[12]宋曉猛，張建云，占車生，等.氣候變化和人類活動(dòng)對(duì)水文循環(huán)影響研究進(jìn)展[J].水利學(xué)報(bào)，2013，44（7）：779790.（SONG X M，ZHANG J Y，ZHAN C S，el at.Review for impacts of climate change and human activities on water cycle[J].Journal of Hydraulic for Engineering，2013，44（7）：779790.（in Chinese）） DOI：10.13243/j.cnki.slxb.2013.07001.

[13]邱臨靜.氣候要素變化和人類活動(dòng)對(duì)延河流域徑流泥沙影響的評(píng)估[D].楊凌：西北農(nóng)林科技大學(xué)，2012.（QIU L J.Evaluation of impacts of the climaticvariable changes and human activity on the runoff and sediment in the Yanhe Watershed[D].Yangling：North A& F University，2012.（in Chinese））.

[14]陳立華，葉明，王佳佳，等.西江流域干流徑流演變特征與趨勢(shì)分析[J].廣西大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，39（4）：920927.（CHEN L H，YE M，WANG J J，el at.Analysis on characteristics and tendency of annual runoff changes of main stream of Xijiang River Basin.[J].Journal of Guangxi University（Nat Sci Ed），2014，39（4）：920927.（in Chinese）） DOI：10.13624/j.cnki.issn.10017445.2014.04.005.

[15]汪麗娜，陳曉宏，李粵安，等.西江流域徑流演變規(guī)律研究[J].水文，2009，29（4）：2225.（WANG L N，CHEN X H，LI Y A，el at.Rules of runoff variation of Xijiang River Basin[J].Journal of China Hydrology，2009，29（4）：2225.（in Chinese）） DOI：10.3969/j.issn.10000852.2009.04.005.

[16]王博，邱現(xiàn)訛，郜軍艷.灤河流域主要?dú)庀笠丶皬搅髯兓厔?shì)分析[J].人民黃河，2011，33（8）：4649.（WANG B，QIU X E，GAO J Y.LongTerm trends of major climate variables and runoff in Luanhe River Basin[J].Yellow River，2011，33（8）：4649.（in Chinese）） DOI：10.3969/j.issn.10001379.2011.08.016.

[17]張榮，段震，王文民，等.變化環(huán)境下大汶河水資源演變趨勢(shì)及影響分析[J].人民黃河，2015，37（1）：6769.（ZHANG R，DUAN Z，WANG W M，el at.Water resources evolution and its influence analysis in Dawen River Basin under the changing environment[J].Yellow River，2015，37（1）：6769.（in Chinese）） DOI：10.3969 /j.issn.10001379.2015.01.017.

[18]田甜，黃強(qiáng)，郭愛軍，等.基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的渭河流域干旱演變特征分析[J].水力發(fā)電學(xué)報(bào)，2016，35（2）：1627.（TIAN T，HUANG Q，GUO A J，el at.Drought evolution characteristics in Wei River Basin based on standardized precipitation evapotranspiration index[J].Journal of Hydroelectric Engineering，2016，35（2）：1627.（in Chinese）） DOI：10.11660/slfdxb.20160202.

[19]代穩(wěn)，呂殿青，李景保，等.19512014年洞庭湖水沙階段性演變特征及驅(qū)動(dòng)因素分析[J].水土保持學(xué)報(bào)，2017，31（2）：142150.（DAI W，LV D Q，LI J B，el at.Evolution characteristics and driving factors of runoff and sediment changes in Dongting Lake during 19512014[J].Journal of Soil and Water Conservation，2017，31（2）：142150.（in Chinese）） DOI：10.13870/j.cnki.stbcxb.2017.02.02.

[20]趙天力，雷德義，胡彩虹.新村流域徑流變異診斷及其影響因素定量分析[J].人民黃河，2016，38（9）：1923.（ZHAO T L，LEI D Y，HU C H.Diagnose and analysis of variation and its inflect factors in Xincun Watershed[J].Yellow River，2016，38（9）：1923.（in Chinese）） DOI：10.3969/j.issn.10001379.2016.09.006.

[21]燕華云，楊貴林，汪青春.長(zhǎng)江源區(qū)徑流年內(nèi)分配時(shí)程變化規(guī)律分析[J].冰川凍土，2006，28（4）：526529.（YAN H Y，YANG G L，WANG Q C，Change of annual runoff distribution in the headwaters of the Yangtze River[J].Journal of Glaciology an Geocryology，2006，28（4）：526 529.（in Chinese）） DOI：10.3969/j.issn.10000240.2006.04.010.

[22]MANN H B.Nonparametric tests against trend[J].Econometrica，1945，13（3）：245259.

[23]KENDALL M G.Rank Correlation Method[M].London：Charles Griffin，1975.

[24]李瑋，段利民，劉廷璽，等.19612015年內(nèi)蒙古高原內(nèi)陸河?xùn)|部流域極端降水時(shí)空變化特征分析[J].資源科學(xué)，2017，39（11）：21532165.（LI W，DUAN L M，LIU T X，el at.Spatiotemporal variations of extreme precipitation from 1961 to 2015 in the Eastern Inland River Basin of Inner Mongolian Plateau[J].ResourceS Science，2017，39（11）：21532165.（in Chinese）） DOI：10.18402/resci.2017.11.13.

[25]XU J H，CHEN Y N，LI W H，el at.Longterm trend and fractal of annual runoff process in mainstream of Tarim River[J].Chinese Geographical Science，2008，18（1）：7784.DOI：10.1007/s1176900800776.

[26]曹明亮，張弛，周惠成，等.豐滿上游流域人類活動(dòng)影響下的降雨徑流變化趨勢(shì)分析[J].水文，2008，28（5）：8689.（CAO M L，ZHANG C，ZHOU H C，el at.Analysis of changing trend of rainfall and runoff under the influence of human activities in upstream of Fengman Watershed[J].Journal of China Hydrology，2008，28（5）：8689.（in Chinese）） DOI：10.3969/j.issn.10000852.2008.05.021.

[27]姚章民，杜勇，張麗娜.珠江流域暴雨天氣系統(tǒng)與暴雨洪水特征分析[J].水文，2015（2）：8589.（YAO Z M，DU Y，ZHANG L N.Analysis of strom weather system and flood characteristic in Pearl River Basin[J].Journal of China Hydrology，2015（2）：8589.（in Chinese）） DOI：10.3969/j.issn.10000852.2015.02.018.