降水與人類活動(dòng)對小清河上游徑流量的影響

1、<p>　　降水與人類活動(dòng)對小清河上游徑流量的影響</p><p>　　摘要：根據(jù)小清河流域黃臺(tái)橋水文站1965～2013年徑流序列資料以及黃臺(tái)橋雨量站的降水量數(shù)據(jù)，采用5a滑動(dòng)平均法和Mann-Kendall非參數(shù)檢驗(yàn)法，分析了小清河上游徑流量和降水量的變化特征；在此基礎(chǔ)上，利用徑流量與降水量的雙累積曲線對降水-徑流序列進(jìn)行了突變點(diǎn)分析，并對突變點(diǎn)前后徑流量與降水量的Pearson相關(guān)系數(shù)進(jìn)行了分析與

2、檢驗(yàn)。結(jié)果表明：小清河上游徑流量與降水量均呈現(xiàn)增加趨勢，但徑流量較降水量的增加趨勢顯著。將小清河上游徑流演變過程分為1965～1974年，1975～2003年和2004～2013年3個(gè)階段，1965～1974年為天然徑流階段，1975-2013年為受人類活動(dòng)影響顯著階段。在徑流變化過程中，徑流量與降水量的Peason相關(guān)系數(shù)逐漸減小，而人類通過排污、引水等活動(dòng)加大了對徑流量的影響程度，成為影響徑流變化的決定性因素?？茖W(xué)地區(qū)分降水與人類活

3、動(dòng)對徑流變化的影響，可以為研究流域水文循環(huán)過程提供借鑒，對實(shí)現(xiàn)水資源可持續(xù)利用也具有現(xiàn)實(shí)意義。 </p><p>　　關(guān)鍵詞：徑流變化；降水；人類活動(dòng)；小清河 </p><p>　　中圖分類號(hào)：TV121.1文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：16749944（2016）18014706 </p><p><b>　　1引言 </b></p>

4、<p>　　隨著全球氣候變化和人類活動(dòng)的加劇，河川徑流發(fā)生了顯著變化。徑流不僅對氣候波動(dòng)響應(yīng)敏感，而且受人類活動(dòng)強(qiáng)度影響顯著。實(shí)測徑流不僅包括天然徑流，還包括人類活動(dòng)影響下的徑流。眾多學(xué)者已經(jīng)對黃河、長江等大江大河及其支流徑流量的變化規(guī)律以及影響因素做了大量的研究＼[2～7＼]。目前，評估氣候變化對徑流影響的方法主要有流域?qū)Ρ确?、時(shí)間序列對比法以及流域水文模型3種[8]。流域?qū)Ρ确椒ㄊ沁x擇氣候相似的研究流域，對比分析氣候變化

5、或人類活動(dòng)的干擾。張永勇[9]等對比分析了氣候變化對三江源地區(qū)徑流變化的影響。時(shí)間序列對比法應(yīng)用統(tǒng)計(jì)方法分析水文氣象數(shù)據(jù)，即建立降水徑流相關(guān)曲線，分析其相關(guān)關(guān)系。孫衛(wèi)國等＼[10～12＼]采用交叉小波變換方法，討論了黃河源區(qū)徑流與區(qū)域氣候變化之間的多時(shí)間尺度相關(guān)。曹建廷＼[13，14＼]等采用Mann-Kendall趨勢檢驗(yàn)法演示了長江、黃河源區(qū)的徑流變化過程。流域水文模型既對流域內(nèi)發(fā)生的水文過程進(jìn)行模擬計(jì)算建立起的數(shù)學(xué)模型。劉潔等＼[

6、15～18＼]通過建立流域降雨-徑流模型分析了降雨和人類活動(dòng)對流域徑流的影響。隨著濟(jì)南市的發(fā)展，小清河流域社會(huì)經(jīng)濟(jì)活動(dòng)越發(fā)頻繁。在人類活動(dòng)干擾下</p><p><b>　　2研究區(qū)域 </b></p><p>　　小清河源于濟(jì)南西郊睦里閘，干流流經(jīng)濟(jì)南市槐蔭區(qū)、天橋區(qū)、歷城區(qū)，于章丘出市境，流向?yàn)I州市的鄒平，又經(jīng)高青、桓臺(tái)、博興、廣饒，由壽光市的羊角溝注入渤海。小清

7、河濟(jì)南市河段長70.5 km，流域面積2792 km2，是濟(jì)南及沿岸地區(qū)的唯一排洪河道，擔(dān)負(fù)著濟(jì)南市承瀉洪水的重要任務(wù)[19]。黃臺(tái)橋水文站位于小清河上游，控制流域面積351 km2，流域內(nèi)的主要支流有南太平河、北太平河、興濟(jì)河、工商河、西濼河、東濼河、柳行河、全福河。支流主要分布于主河道南岸，多屬雨源型山溪河流。小清河流域?qū)贉貛Ъ撅L(fēng)氣候，降水變率大，約70%的降水集中在6～9月[20]。均在黃臺(tái)橋水文站以上流入小清河主河道。本文主要研

8、究黃臺(tái)橋以上流域。 </p><p>　　3研究方法和數(shù)據(jù)來源 </p><p><b>　　3.1研究方法 </b></p><p>　　3.1.1滑動(dòng)平均法 </p><p>　　滑動(dòng)平均法是趨勢擬合技術(shù)最基礎(chǔ)的方法，相當(dāng)于低波過濾器[21]。用確定時(shí)間序列的平滑值來表示變化趨勢。樣本量為n的時(shí)間序列x的滑動(dòng)平均序列

9、可表示為： </p><p>　　j=1k∑ki=1xi+j-1j=1，2，…，n-k+1 （1） </p><p>　　式中k為滑動(dòng)長度，選擇合適的k可以較好地反應(yīng)時(shí)間序列的變化趨勢。 </p><p>　　3.1.2Mann-Kendall非參數(shù)趨勢檢驗(yàn) </p><p>　　用Mann-Kendall[16]非參數(shù)趨勢檢驗(yàn)法檢測徑流序列

10、變化趨勢，該統(tǒng)計(jì)檢驗(yàn)方法的優(yōu)點(diǎn)是不需樣本服從一定的分布，也不受少數(shù)異常值的干擾，被廣泛應(yīng)用于氣象、水文要素等時(shí)間序列的趨勢和突變分析。 對于具有n個(gè)樣本量的時(shí)間序列x1，x2，…，xn構(gòu)造一秩序列： </p><p>　　sk=∑ki=1rik=2，3，…，n（2） </p><p><b>　　其中 </b></p><p>　　ri=+1，

11、當(dāng)xi>xj </p><p>　　0，當(dāng)xi<xjj=1，2，…，i </p><p>　　sk是時(shí)間序列中第i時(shí)刻的數(shù)值大于j時(shí)刻數(shù)值個(gè)數(shù)的累計(jì)數(shù)。 </p><p>　　在時(shí)間序列隨機(jī)獨(dú)立的假定下，定義統(tǒng)計(jì)量 </p><p>　　UFk=＼[sk-E（sk）＼]var（sk）k=1，2，…，n（3） </p>

12、<p>　　其中，E（sk）、var（sk）是累計(jì)數(shù)sk的均值和方差，可由下式算出： </p><p>　　E（sk）=k（k-1）4 </p><p>　　var（sk）=k（k-1）（2k+5）72k=2，3，…，n （4） </p><p>　　按時(shí)間序列x逆序xn，xn-1，…，x1，再重復(fù)以上過程，可得UBk。 </p><

13、p>　　若UFk或者UBk的值大于0，則表明序列呈上升趨勢，小于0則表明呈下降趨勢。當(dāng)它們超過臨界直線時(shí)，表明上升或下降趨勢顯著[21]。 </p><p>　　3.1.3雙累積曲線和累積距平法 </p><p>　　雙積累曲線方法是目前用于水文氣象要素一致性或長期演變趨勢分析中最簡單、最直觀、最廣泛的方法[17]。當(dāng)兩個(gè)變量累積值之間直線的斜率發(fā)生改變，那么斜率發(fā)生突變的年份就是兩

14、個(gè)變量累積關(guān)系發(fā)生突變的時(shí)間[22]。 </p><p>　　累積距平也是一種常用的、由曲線直觀判斷變化趨勢的方法。對于序列x，其某一時(shí)刻t的累積距平表示為： </p><p>　　i=∑ti=1（xi-）t=1，2，…，n （5） 　　累積距平曲線呈上升趨勢，表示距平值增加，曲線呈下降趨勢則表示距平值減小。從曲線的上下起伏，可以判斷序列長期顯著的演變趨勢及持續(xù)性變化[21]。 <

15、/p><p><b>　　3.2資料來源 </b></p><p>　　小清河流域徑流量選用山東省水文局提供的黃臺(tái)橋水文站1965～2013年逐年實(shí)測徑流量資料，降水量選用山東省水文局提供的黃臺(tái)橋雨量站1965～2013年的逐年降水量。 </p><p>　　朱琳，等：降水與人類活動(dòng)對小清河上游徑流量的影響地理與規(guī)劃 </p><

16、;p><b>　　4結(jié)果分析 </b></p><p>　　3.1徑流量變化特征 </p><p>　　根據(jù)1965～2013年小清河黃臺(tái)橋水文站的實(shí)測年徑流量資料，選取5a線性滑動(dòng)平均法對小清河年徑流量變化情況進(jìn)行統(tǒng)計(jì)，結(jié)果如圖2，小清河徑流量呈現(xiàn)出不規(guī)則的周期性波動(dòng)。20世紀(jì)70年代和90年代前期徑流量呈增長趨勢，80年代和90年代后期徑流量呈減小趨勢。進(jìn)入

17、21世紀(jì)以來，小清河徑流量增長趨勢明顯，年均增長率為10.61%。統(tǒng)計(jì)資料顯示，1965～2013年的多年平均徑流總量3.03×108 m3，最大年徑流量5.64×108 m3（2012年），最小年徑流量0.99×108 m3（1968年），極值比5.66。 </p><p>　　由Mann-Kendall趨勢檢驗(yàn)可以看出（圖3），1974年后小清河徑流量變化的整體趨勢是增加的，在2

18、0世紀(jì)70年代后期，UFk值超過顯著性水平0.05（u0.05=±1.96）趨勢線，表明在該時(shí)間段內(nèi)，徑流量增長趨勢顯著。2005年后，UFk值甚至超過了顯著性水平0.001（u0.001=±2.56）趨勢線，徑流量的增長趨勢十分顯著。 </p><p>　　圖3小清河上游徑流量Mann-Kendall統(tǒng)計(jì)量曲線 </p><p>　　4.2降水量對徑流量的影響 <

19、;/p><p>　　4.2.1降水量的變化趨勢 </p><p>　　根據(jù)1965～2013年小清河黃臺(tái)橋雨量站的降水量資料，選取5a滑動(dòng)平均對小清河降水量變化情況進(jìn)行統(tǒng)計(jì)，結(jié)果見圖4。在研究時(shí)間序列內(nèi)，降水量呈現(xiàn)不規(guī)則的周期性波動(dòng)，降水量的最小值為297.5 mm（1968年），最大值為1076.7 mm（2004年）。20世紀(jì)60、70年代，降水量增加趨勢顯著；80、90年代，降水量整體變

20、化趨勢比較平穩(wěn)。2004年后，降水量有下降趨勢。整體來講，小清河的降水量呈現(xiàn)不太顯著的增加趨勢。 </p><p>　　4.2.2降水量與徑流量相關(guān)性分析 </p><p>　　由徑流量與降水量的距平百分率曲線可以看出（圖5），二者的變化趨勢大體一致，但并非完全吻合。1974年以前，徑流量的年際波動(dòng)幅度較大，1974～2003年二者年際波動(dòng)幅度相差不大，2004年后徑流量的波動(dòng)幅度大于降水

21、量的波動(dòng)幅度?？梢?，降水并非是影響徑流量的唯一因素，非降水因素也影響著小清河徑流量的變化。2004年后，徑流量距平百分率超過40%，而降水量在平均值附近波動(dòng)，變化不大。表明在非降水因素的影響下，小清河上游的徑流量呈顯著的增加趨勢。 </p><p>　　采用Mann-Kendall非參數(shù)檢驗(yàn)和雙累積曲線法識(shí)別徑流序列發(fā)生突變的臨界年份。Mann-Kendall檢驗(yàn)結(jié)果如圖2，UFk與UFb兩條曲線的交點(diǎn)出現(xiàn)在19

22、74年，表明徑流序列在1974年發(fā)生突變。由圖6可見，雙累積曲線的斜率在1974年、2003年發(fā)生明顯變化，且3個(gè)階段線性擬合趨勢線的斜率呈增加趨勢。3條趨勢線R2均超過0.99，線性擬合程度較高。因此，可將降水-徑流變化過程分為3個(gè)階段：1965～1974年、1975～2003年和2004～2013年。第一階段為天然徑流階段，人類活動(dòng)對小清河的干預(yù)較少，其徑流量主要受降水等自然條件的影響。第二、三階段為受人類活動(dòng)影響顯著階段，雙累積曲

23、線表明，在人類活動(dòng)的干預(yù)下，小清河徑流量呈增加趨勢。 </p><p>　　不同時(shí)期徑流量與降雨量的相關(guān)關(guān)系發(fā)生變化，如圖7。隨著時(shí)間序列的推移，徑流量與降水量的關(guān)系點(diǎn)向右推移，相同量級(jí)降水量下的徑流量明顯增加。因此，將1965～1974年劃分為受人類活動(dòng)影響較小的天然徑流階段是有意義的。 </p><p>　　4.3人類活動(dòng)對徑流量的影響 </p><p>　　1

24、965-1974年為天然徑流階段，小清河徑流量主要來自于大氣降水，受人類活動(dòng)影響相對較小。第二、三階段為受人類活動(dòng)影響顯著階段，徑流量與降水量的相關(guān)性明顯減小。對小清河多年平均月徑流量和降水量進(jìn)行統(tǒng)計(jì)，如圖8、圖9所示。2004～2013年小清河月平均徑流量明顯增加，但月平均降水量的增加趨勢卻不顯著。且降水主要集中在夏季，徑流量的年內(nèi)變化較降水量緩和。這在一定程度上表明，人類活動(dòng)增加了小清河尤其是枯水期的徑流量，緩和了小清河徑流量的年內(nèi)

25、變化。人類干擾是造成小清河徑流量增大的主要原因，是影響徑流量的決定性因素。 </p><p>　　影響小清河上游徑流量的人類活動(dòng)主要包括：城市生活污水、工業(yè)廢水排入小清河和引水補(bǔ)源，即從其他水庫或湖泊向小清河調(diào)水。 </p><p>　　20世紀(jì)70年代以來，隨著濟(jì)南市小清河流域內(nèi)工農(nóng)業(yè)生產(chǎn)的迅速發(fā)展，工業(yè)廢水和生活污水的排放量逐年增加[23]，成為小清河水體的主要來源。以2000年為例，

26、濟(jì)南市生活污水排放量為1.3×108 m3，占全市廢水排放量的65.1%[23]，近90%的城市污水未經(jīng)處理直接進(jìn)入小清河。黃臺(tái)橋斷面接納了濟(jì)南市大部分工業(yè)廢水和生活污水，主要排放渠道是興濟(jì)河、東濼河、西濼河、工商河、全福河和柳行河六大排污系統(tǒng)。 </p><p>　　為保證小清河水質(zhì)要求，滿足其景觀功能，需要對小清河進(jìn)行補(bǔ)源：從臥虎山水庫放水，通過玉符河向小清河引水；引黃河水入玉符河，澄清后經(jīng)睦里閘直

27、接引入小清河；經(jīng)濟(jì)平干渠引東平湖水入小清河，2005年12月29日至2009年11月18日，東平湖共向小清河調(diào)水5次，調(diào)水量1.005×108 m3。通過對1965～2013年小清河流域月平均徑流量的分析發(fā)現(xiàn)，2004～2013年與1965～1974年相比，1月份徑流量平均增長率為1.91%，4月份平均增長率為1.62%，7月份平均增長率為0.89%，10月份平均增長率為1.61%。小清河枯水期徑流量增長趨勢較汛期顯著，因?yàn)闈M

28、足小清河的水質(zhì)和景觀需求，在枯水期通過引水工程增大了對小清河的引水量。 </p><p><b>　　5結(jié)論 </b></p><p>　　（1）1965～2013年小清河流域徑流量呈顯著的增加趨勢，多年平均增長率為2.07%。21世紀(jì)之后，徑流量增加趨勢更為明顯，年均增長率為10.61%。 </p><p>　?。?）利用Mann-Kenda

29、ll非參數(shù)檢驗(yàn)和雙累積曲線法對小清河徑流序列進(jìn)行突變點(diǎn)分析，發(fā)現(xiàn)小清河徑流序列出現(xiàn)1974、2003年兩個(gè)突變點(diǎn)，遂將小清河流域的徑流演變過程分為1965～1974年、1975～2003年、2004～2013年3個(gè)階段。1965～1974年為天然徑流階段，小清河平均年徑流量2.07×108 m3。后兩個(gè)階段為受人類活動(dòng)影響顯著階段，在人類活動(dòng)影響下，小清河徑流量增加趨勢顯著，1975～2003年平均年徑流量2.62×

30、108 m3，2004～2013年平均年徑流量4.88×108 m3。 </p><p>　?。?）小清河流域徑流量的變化主要是由氣候變化和人類活動(dòng)引起的。在徑流演變過程中，徑流量與降水量的相關(guān)性逐漸減小，人類活動(dòng)對徑流量的影響逐漸增大，并成為影響徑流量的主導(dǎo)因素。人類活動(dòng)主要通過向小清河排放工業(yè)廢水、生活污水、處理中水以及小清河引水工程來影響小清河徑流量。 </p><p>&

31、lt;b>　　參考文獻(xiàn)： </b></p><p>　　[1]王國慶，張建云，劉九夫，等.氣候變化和人類活動(dòng)對河川徑流影響的定量分析[J].中國水利，2008（2）：55～58. </p><p>　　[2]王彥君，王隨繼，蘇騰.降水和人類活動(dòng)對松花江徑流量變化的貢獻(xiàn)率[J].自然資源學(xué)報(bào)，2015，30（2）：304～314. </p><p>

32、　　[3]陳麗君.影響海河流域徑流變化的因素分析[J].綠色科技，2010（8）：150～151. </p><p>　　[4]王立紅，陳娟，張慧.濟(jì)南市錦繡川流域景觀變化對河川徑流的影響研究[J].綠色科技，2011（8）：154～157. </p><p>　　[5]劉蜀鄂，雷琳，高雁.牧羊河流域出水徑流年內(nèi)分配格局變化研究[J].綠色科技，2012（2）：203～204. </