丁 文 荣

(云南师范大学 旅游与地理科学学院, 云南 昆明 650500)

变化环境下滇中地区典型流域水资源演变特征

丁 文 荣

(云南师范大学 旅游与地理科学学院, 云南 昆明 650500)

[目的] 以典型流域为研究对象，揭示变化环境下滇中地区水资源的演变特征，为区域水资源综合管理服务。 [方法] 以滇中地区蜻蛉河流域1963—2012年的月值径流实测资料为基础，采用线性倾向估计、Mann-Kendall非参数检验、小波分析和R/S分析等作为主要研究方法。 [结果] (1) 蜻蛉河流域径流量年内分配极不均匀，主要集中在夏秋两季，而冬春两季占全年的百分比均较小； (2) 自20世纪60年代以来，蜻蛉河流域径流量总体呈现出了减少的趋势，速率为4.80×107m3/10 a，但其变化趋势的尚未达到0.05的显著性水平； (3) 研究时段内蜻蛉河流域径流量的变化经历了“多—少—多—少”的演变过程，2010年为一个典型突变年份； (4) 从多时间尺度特征来看，蜻蛉河流域径流量变化局域性特点较为突出，没有贯穿整个研究时段的周期性波动，相对较为明显的周期性变化为15 a； (5) 蜻蛉河流域年径流量的变化存在比较显著的赫斯特现象，未来持续减少的概率较大。 [结论] 蜻蛉河流域水资源季节性分布不均，年际波动幅度大，未来一段时间内将继续维持减少趋势，需提前做好合理规划。

变化环境; 滇中地区; 水资源; 蜻蛉河

文献参数： 丁文荣.变化环境下滇中地区典型流域水资源演变特征[J].水土保持通报，2017,37(2)：274-277.DOI:10.13961/j.cnki.stbctb.2017.02.041； Ding Wenrong. Evolutions of Water Resources in Typical Basin of Central Yunnan Province Under Changing Environment[J]. Bulletin of Soil and Water Conservation, 2017,37(2)：274-277.DOI:10.13961/j.cnki.stbctb.2017.02.041

近百年来，地球系统正经历着以气候变暖为主要特征的显著变化，加之剧烈的人类活动影响，不同程度改变了水文循环过程，严重威胁着人类生存和社会发展。由此，变化环境下的水资源效应受到各国学者的广泛关注，成为当今国际水科学研究的前沿课题[1]。自20世纪80年代以来，关于变化环境对水文过程影响的研究已经在国内外积累了丰富的研究成果[2]，为缓解乃至解决当今人类社会所面临的水资源、水环境、水生态和水灾害问题提供了途径。滇中地区是云南省的政治、经济和文化中心与交通枢纽。2015年的统计结果表明，该区域人口1 702.8万人，占该省总人口的37.4%；GDP占该省的57.8%，人均GDP约为19 400元，是该省人均GDP的1.55倍，已成为云南的名副其实经济核心区。滇中地处金沙江、南盘江、红河3大水系分水岭地带，多年来在气候变化和土地利用/覆被变化(LUCC)等人类活动的双重影响下，滇中地区流域水资源时空分布发生了深刻的变化，资源型、工程型和水质型缺水日益凸显，不仅影响区内人民生活水平的提高，也严重制约着全省社会经济的可持续发展[3-4]。鉴此，本研究以滇中地区典型区域蜻蛉河流域为研究对象，采用多种时间序列分析方法，系统探讨该流域水资源的时间变化规律，以期为区内水资源的综合管理提供决策参考。

1 研究区概况

蜻蛉河为龙川江左岸最大的一级支流，发源于楚雄彝族自治州姚安县太平乡海拔2 525 m的黎梅山，进入大姚县境后折向东北方向，最后在元谋县黑泥坡村汇入龙川江。河流全长132 km，多克(二)站以上流域面积3 546 km2，位于101°04′—101°47′E和25°20′—26°10′N，地处北亚热带季风气候区域，属典型亚热带干燥气候，具有气候温和，日照充足，干湿季明显，无霜期长等气候特点。年均气温15.6 ℃，日照2 526 h，无霜期年均309 d，相对湿度65%。多年平均降水796.3 mm，雨日113 d。研究区基岩主要是三叠纪页岩、砂岩和少量的花岗岩、石灰岩和第四纪沉积物，地貌表现为中山、丘陵、盆地相间分布的格局。土壤类型以北亚热带山原红壤为主，在海拔1 700 m以下的干热河谷区发育有燥红土。植被主要为亚热带常绿阔叶林，自20世纪50年代以来，受人类活动的影响，原生植被已消失殆尽，目前云南松人工次生林广泛分布，是区内荒山绿化和迹地更新的主要林种。

2 数据与方法

所用径流量资料为蜻蛉河下游国家基本水文站多克(2)站1963—2012年的月值径流实测序列，数据来源于云南省水文水资源局，经过了可靠性、一致性与代表性检验。其中1988年10—12月缺测值以多年平均值进行插补，分析时间序列长度为50 a。水资源变化趋势采用线性倾向方程来表示，即：

Y=α0+α1t

(1)

式中：Y——待分析资料序列；t——时间；α1——线性趋势项，若该值为正，表明为增加趋势，反之为减少趋势，α1×10表示水资源每10 a的变化量。

突变过程检测采用Mann-Kendall非参数检验方法，该法因不需要确定样本服从的分布，也不受外界异常因素的影响而被广泛应用[5]。

水资源时间序列变化的周期性分析采用小波分析方法，该方法通过伸缩平移运算对信号逐步进行多尺度细化，自动适应时频信号分析的要求，从而可聚焦到信号的任意尺度，最终达到高频处时间细分，低频处频率细分的目的[6]。

水资源的未来演变趋势采用了R/S分析方法，该方法假设一个时间序列ξ(t)，t=1，2，3，…，对于任意整数τ>1，先指定时均序列<ξ>τ，t=1，2，3，… ，然后计算累计离差X(t,τ)、级差R(τ)和标准差S(τ)。赫斯特曾证明了R(τ)/S(τ)随τ的增长呈1/2的指数率关系[7]，在比值R(τ)/S(τ)=R/S中，R随τ的增长呈发散状态，即：

(2)

式中：H——赫斯特指数，可以揭示水资源时间序列的未来演变趋势。

3 结果与分析

3.1 年内变化特征

地处云南高原的北部边缘，蜻蛉河流域地表径流的主要补给源为大气降水。冬春季水汽源为中纬度干旱大陆，空气干燥降水稀少，而夏秋季则有孟加拉湾和北部湾的暖湿气流，雨量充沛[8]。受到降水这种年内变化特征的控制，蜻蛉河流域6月上旬进入汛期而10月下旬汛期基本结束，河川径流量主要集中在每年的6—10月，该时段径流量占全年的81.06%，而每年的1—5月与11—12月的径流量只仅全年的18.94%。年内径流量的变差系数Cv值为1.12，表明径流年内分配极不均匀。按滇中地区气候季节划分的时间节点，将蜻蛉河流域径流量对应地划分为春夏秋冬四季。分别是春季3—5月，夏季6—8月，秋季9—11月和冬季12月至翌年2月。径流量的年内季节变化计算结果表明，蜻蛉河流域径流量主要集中在夏秋两季，分别占全年径流量的49.62%和42.25%，而冬春两季占全年的百分比均较小。尤其是春季，由于此时正值农业大春播种，土壤缺墒，加之气温缓慢升高蒸发加剧，而该季节地表径流量仅占全年的2.69%，成为流域内一年中缺水最为突出的季节。

3.2 年际变化趋势

自20世纪60年代以来，蜻蛉河流域径流量总体呈现出了减少的趋势，变化速率为4.80×107m3/10 a(图1)。特别是进入21世纪以来减少更为明显，尤其是2011年径流量仅为5.00×107m3，相当于丰水年份1966年1.39×109m3的3.60%。然而，采用Mann-Kendall检验对年径流量的变化趋势进行显著性检验后发现，其变化趋势的Z值为-1.88，尚未达到0.05的显著性水平。

图1 蜻蛉河流域年径流量的变化趋势

从不同年代的变化也可以看出(表1)，蜻蛉河流域径流量冬春两季占全年的比例均较少，春季尤为突出，而夏秋两季占全年的比例均超过90%。而从夏秋季来看，两个季节占全年的比例相差不大，夏季略多。此外，从20世纪60年代至今，春夏季两季径流量占全年的比例略有增加，秋冬两季略有减少，但总体而言不同季节径流量占年径流总量的比例并没有发生明显的变化。

表1 研究区不同季节径流量占年径流总量的比例 %

3.3 突变分析

采用Mann-Kendall非参数检验方法进行突变分析时，当统计量UF和UB构成的曲线超过了正负临界值直线，则表明原序列具有显著的上升或下降趋势，反之趋势不明显。若UF和UB在正负临界值之间出现交点，则表明原序列存在跳跃变异(突变)，且两曲线相交点的时刻即为突变开始的时刻。从图2中UF曲线的变化可以看出，自20世纪60年代中期至1983年，蜻蛉河流域径流量以减少趋势为主，1983—1997年为小幅波动而没有明显的趋势，1997年以来有增加的趋势并维持到2002年左右，然后持续减少至今。研究时段内径流量的变化经历了“多—少—多—少”的演变过程。此外，从图2还可以看出，UF和UB曲线在2010年有一个交点，2010年后逆序列UB曲线超过了UF曲线且UF曲线呈下降趋势，表明径流量呈现出减少趋势。此外，该交点在正负临界值之间，表明2010年为一个突变年份。经过计算也发现，2010年后径流量由多年平均值4.72×108m3骤减至1.18×108m3左右。这与2009年以来中国西南地区遭遇百年不遇的极端干旱是相对应的[9-10]。

图2 蜻蛉河流域年径流量的M-K检验

3.4 演变周期

研究表明，用复值Morlet小波进行时间序列分析时，小波系数实部的变化趋势与待分析信号的起伏基本一致，而Morlet小波系数模值的大小代表某一尺度成分的多少，反映了它所对应时段的周期性强弱[11]。由于小波系数模值都是大于或等于零的实数，因此周期性的强弱可通过模值的大小反映出来，模值越大则周期性越显著，反之不显著。由图3蜻蛉河流域年径流量的复值Morlet小波系数实部可以看出，在研究时段内，蜻蛉河流域径流量变化局域性特征较为突出，没有贯穿整个研究时段的周期性波动。较为明显的是15 a左右的周期性变化，存在2个明显的丰水中心和少水中心年份，丰水中心年份分别是1970和1985年，2个少水中心年份分别为1977和1992年。蜻蛉河流域径流量变化的复值Morlet小波变换系数模反映了径流量变化在各时间尺度的局部化特征(图3)。由图3可以看出，15 a时间尺度的振荡在研究时段内相对明显。而其它尺度则局部化特征很强，如20 a尺度的振荡最弱，5 a尺度的振荡在局部性特征尤为突出。

图3 蜻蛉河流域径流量的复值Morlet小波系数

3.5 演变趋势

Mandelbrot等[12]对赫斯特指数H的大小与未来发展趋势之间的关系进行了探讨，表明H=0.5是一个临界值，当0.5

4 结 论

(1) 蜻蛉河流域地表径流年内时间分配极不均匀，主要集中在夏秋两季，而冬春两季占全年的比例均较小。

(2) 自20世纪60年代以来，蜻蛉河流域径流量总体呈现出了减少的趋势，速率为4.80×107m3/10 a，但其变化趋势的尚未达到0.05的显著性水平。

(3) 自20世纪60年代中期至今，蜻蛉河流域径流量的变化经历了“多—少—多—少”的演变过程，2010年为一个明显的突变年份。

(4) 在研究时段内，蜻蛉河流域径流量变化局域性特征较为突出，没有贯穿整个研究时段的周期性波动，相对而言较为明显的是15 a左右的周期性变化。

(5) 蜻蛉河流域年径流量的赫斯特指数大于0.5，表明存在比较显著的赫斯特现象，故将维持前期的演变趋势。自20世纪60年代以来蜻蛉河流域径流量为减少趋势，故未来持续减少的概率较大。

[1] Oki T, Kanae S. Global hydrological cycles and world water resources[J]. Science, 2006,313(5790):1068-1072.

[2] Boyer C, Chaumont D, Chartier I, et al. Impact of climate change on the hydrology of St. Lawrence tributaries[J]. Journal of Hydrology, 2010,384(1/2):65-83.

[3] 冯少辉,李靖,朱振峰,等.云南省滇中地区水资源脆弱性评价[J].水资源保护,2010,26(1):13-16.

[4] 伍立群,王超.滇中地区水资源特性与水安全问题分析[J].水文,2009,29(S1):28-30.

[5] Burn D H, Elnur M A H. Detection of hydrologic trends and variability[J]. Journal of Hydrology, 2002,255(1/4),107-122.

[6] 吴创收,杨世伦,黄世昌,等.1954—2011年珠江入海水沙通量变化的多尺度分析[J].地理学报,2014,69(3):422-432.

[7] Lloyd E H, Hurst H E, Black R P, et al. Long-term storage: An experimental study[J]. Journal of the Royal Statistical Society, 1966, 129(4):591.

[8] 曹杰,李华宏,姚平,等.北半球夏季印度洋和太平洋水汽交汇区及其空间分异规律研究[J].自然科学进展,2009,19(3):302-309.

[9] 贺晋云,张明军,王鹏,等.近50年西南地区极端干旱气候变化特征[J].地理学报,2011,66(9):1179-1190.

[10] 黄慧君,李庆红,高月忠.云南省2009/2010年秋冬季罕见干旱诊断分析[J].热带地理,2011,31(1):28-32.

[11] 尤卫红.气候变化的多尺度诊断分析和预测的多种技术方法研究[M].北京:气象出版社,1998.

[12] 王文圣,丁晶,李跃清.水文小波分析[M].北京:化学工业出版社,2005.

Evolutions of Water Resources in Typical Basin of Central Yunnan Province Under Changing Environment

DING Wenrong

(CollegeofTourismandGeographyScience,YunnanNormalUniversity,Kunming,Yunnan650500,China)

[Objective] A typical basin was taken as study area, we demonstrated the evolution of water resources in Central Yunnan Province, to supply scientific support for integrated water resources management. [Methods] Runoff data was collected in Qingling River basin of central Yunnan Province. Linear trend estimation, Mann-Kendall non parameter test, wavelet analysis method andR/Sanalysis were used. [Results] (1) Runoff is extremely non-uniform at the annual scale, and concentrating on the summer and winter. (2) Since 1960s, runoff of Qingling River basin showed a decreasing trend with a of 4.80×107m3/10 a, but the trend not yet reached the significant level of 0.05. (3) The periodical evolution of runoff experienced “rich—lack—rich—lack” processes, in 2010, an abrupt change occurred. (4) From the aspect of multi-time scale variation, it was locally prominent but not over all covered. Comparatively a 15 year-periodical variation was obvious. (5) A significant Hurst phenomenon existed for annual runoff of Qingling River basin, whereby the decrease probability in the future was great. [Conclusion] Water resources in Qingling River basin are unevenly distributed in seasonal and inter-annual time scale. A decrease trend will continue in the next period, hence rational planning is important.

changing environment; Central Yunnan Province; water resource; Qingling River

2016-07-18

2016-08-08

云南省教育厅重点项目“金沙江下游龙川江流域水沙输移过程及机制研究”(2015Z056); 国家自然科学基金项目(41101099； 41261044)

丁文荣(1979—),男(彝族),云南省昆明市人,博士,副教授 ,主要从事水文水资源方面的研究工作。E-mail:dingwenrong@163.com。

A

1000-288X(2017)02-0274-04

TV312.9, K903