潘娅英, 骆月珍， 王亚男, 张 青, 朱占云

(浙江省气象服务中心, 杭州 310017)

气候变化和人类的活动是影响河川径流变化的最主要的两大因素，由于气候因素会造成降水的改变，从而影响地表、地下水资源的多寡，进而引起径流量的丰枯变化。同时人类活动也会通过改变下垫面而影响河川的产汇流机制，径流量也随之改变。

目前，在水文和气象专业的很多学者对降水和径流的变化都做了较多的有意义的研究，比如孙洋洋[1]、潘彬[2]、张建兴[3]等对黄河上、中、下游的径流变化特征做了研究，得出黄河各段的径流量均有明显减少的趋势，且减少趋势显著。宋小燕等[4]认为松花江的径流在100年时间尺度上并未表现出显著的增加或减小的趋势，但径流量的年际变化具有明显的阶段性特征。潘扎荣等[5]的分析表明淮河干流径流量总体表现出递减的变化趋势。郭爱军[6]、刘伟[7]、古丽孜巴·艾尼瓦尔[8]等利用小波分析、Copula函数等方法探讨了泾河、清水河等流域径流量的周期特征以及与径流、降水的演变规律。从已有的研究来看，径流量的变化与降雨的关系非常密切。

本文选取新安江流域近57 a的径流和降水资料，诊断分析降雨、径流的关系特征，该研究成果有助于新安江流域合理进行水资源的优化开发和高效利用，对区域的经济发展和人民生命财产安全的保障均有重要的实际意义。

1 研究区概况

新安江发源于安徽省休宁县与江西省交界处的五股尖山，东入浙江省西部，经淳安至建德与兰江汇合后为钱塘江干流桐江段、富春江段，东北流入钱塘江，是钱塘江正源。干流长373 km，流域面积11 720 km2。全区面积大，地貌类型多样，平原、岗地、丘陵和山地俱全，各个副区间差异显著。新安江流域属于亚热带季风气候，冬夏季风交替明显，由于季风强度出现的时间年际变化较大，流域内易出现旱、涝、高温、低温等自然灾害。

新安江流域内的新安江水库建立于20世纪50年代，装机容量达216.26亿m3，为多年调节水库。流域范围117.64°—119.34°E，29.19°—30.23°N，控制流域面积10 480 km2，占新安江流域面积的89.4%。是一座以发电为主，兼防洪、旅游、饮用水源、养殖以及工农业用水等诸多功能为一体的大型水库。

2 资料和方法

2.1 资 料

本文选用新安江水库1960—2016年共57 a的实测入库流量资料以及流域内4个气象站点资料和10个水文站的降雨资料作为分析资料。流域面雨量通过各站点的降雨资料经泰森多边形法计算得到。为降雨和径流比较方便，将入库流量换算成径流深，单位为mm。

2.2 分析方法

采用Mann-Kendall检验法分析降雨、径流序列的突变检验分析。Mann-Kendall检验是提取序列变化趋势的有效工具，被广泛应用于气候和水文序列的分析[9-10]。该方法是一种非参数统计检验方法，其优点是不需要样本遵从一定的分布，也不受少数异常值的干扰，更适用于类型变量和顺序变量，计算也比较简便。

另外，利用小波分析法讨论降雨、径流的周期变化特征。小波分析是比较成熟的数学分析方法，广泛地应用于大气、水文等领域，取得了显著成果[11]。采用Morlet小波可以得到Morlet小波系数实部的时频结构图和小波方差，由此可确定序列的主要周期。

3 结果与分析

3.1 降雨、径流年际变化特征分析

新安江流域1960—2016年降雨、径流变化过程如图1所示：(1) 降雨、径流的年际分配不均匀，但两者演变趋势具有一致性，均呈现波动增多的趋势，其中降水增多趋势略微大于径流。(2) 降雨、径流的极大值、极小值均同期出现，分别是在1999年和1978年。其中，1999年新安江流域面雨量达2 410 mm，径流1 715.94 mm，为历年极大值。1978年新安江流域面雨量仅1 040 mm，径流478 mm，为历年极小值。

图1 新安江流域降雨、径流年际变化

3.2 降雨、径流年内特征分析

新安江流域年内各月的降雨、径流如图2所示，新安江流域降雨各月分配极不均匀，从1月起逐渐增加，至6月达最大值319.6 mm，该月雨量占年总雨量的18.9%。7月起雨量迅速减小，至12月为全年最小值仅49.6，该月雨量占年总雨量的2.9%。新安江流域4—7月为年内降雨最为集中的时段，该时段总雨量达931.2 mm，占全年总雨量的55%。

新安江流域径流量从1月起逐渐增加，至6月达最大值216.2 mm，该月径流占年径流的21.7%。7月起径流量迅速减小，至12月达全年最小值21.3 mm，该月径流量仅占年径流量的2.1%。新安江流域4—7月是一年中径流量较为集中的时段，该时段的径流量627.2 mm,占年总径流量的63.1%。

因此，新安江流域降雨和径流的年内的月际分配具有相似特征，均呈一致的单峰型变化。4—5月新安江流域处于春季，天气系统活跃，易出现短时强降水及春季连阴雨天气，降水量较冬季明显增加，径流也明显增大。6—7月上中旬该流域处于雨量最大、最集中的梅汛期，易出现连续性阴雨天气和暴雨天气，径流迅速增大。7月中下旬—9月份新安江流域进入台风干旱期，当有台风影响的时段，该流域容易形成台风降水，多以过程性降水为主；当没有台风影响时该流域均为高温干旱天气，降水多为午后雷阵雨为主，雨量迅速减小，径流也迅速减小。10月—翌年3月是新安江流域的秋冬季节，降水明显减小，为一年之中雨量最少的时段，径流亦如此。4—7月是新安江流域降雨、径流是一年中最大、最为集中的时段，也是该流域洪水灾害多发时期。

图2 新安江流域降雨、径流月际变化

3.3 降雨、径流的突变性分析

采用Mann-Kendall检验法分别对新安江流域降雨、径流序列突变检验分析(图3)，在显著性水平0.05，临界值为±1.96的条件下，1960—1968年UF基本小于0，且在临界值之内，说明降雨(图3A)持续偏少，且有缓慢上升的趋势，但变化不显著。1969—1988年UF值大于0，且与UB交替上下，并小于临界值，说明趋势不显著。1989—2002年，这两段时间UF值呈上升趋势，说明降雨明显增多。2003—2012年UF呈减少趋势，说明降水量明显减少。2013年后UF值又上升，说明降水量增加。在1960—2016年的UF和UB曲线相交于1972年、1978年、1979年、1980年、1984年、1988年、2002年、2009年、2012年，结合UF，UB曲线的变化，可知新安江流域年雨量变化在1971年、1988年、2003年、2009年有4个明显的突变时间点，其中1971年、1988年、2009年是降雨由少转多的突变年份，2003年是降雨由多转少的突变年份。

如图3B所示，1960—1968年UF基本小于0，且在临界值之内，说明该时段处于径流偏少的枯水期，但有缓慢上升的趋势。1969—1978年、1980—2002年，这两段时间UF值先波动后上升，总体来看径流处于丰水期。2003—2009年UF呈减少的趋势，说明径流量逐渐减小的枯水期中。2010—2016年，UF值又上升，说明径流量呈上升趋势，由枯转丰。在1960—2016年的UF和UB曲线相交于1972年、1978年、1979年、2003年、2009年，结合UF，UB曲线的变化，可知新安江流域年径流变化在1971年、2003年、2009年是3个发生突变的时间点，其中1971年、2009年的年径流由枯转丰，2003年的年径流由丰转枯。可见，新安江流域径流与降雨的突变异常时间基本吻合，径流变化受降雨的影响显著。

图3 新安江流域降雨、径流突变检验

3.4 降雨、径流的周期分析

采用Morlet小波对新安江流域年降雨、径流进行小波变换分析，小波实系数等值线图反映了降雨、径流包含的周期与振幅大小。如图4所示，降雨、径流偏多，小波实系数为正值用实线表示；反之，降雨、径流偏少，小波实系数等值线用虚线表示。根据1960—2016年新安江流域的年降雨量、径流的年际变化均存在着2～3 a，5～8 a，10～20 a的多重时间尺度上的周期。新安流域年降雨量具有10～20 a尺度的周期振荡规律，共经历了3次枯—丰循环交替过程。5～8 a的周期变化主要发生在20世纪70—80年代以及2008年以后，先后共经历了4次枯—丰循环交替过程。2～3 a的周期变化主要发生在1985年之前，以及1990—2005年，共经历了9个丰—枯循环交替变化过程。目前正处于大、小周期中的丰水年周期里。

新安江年降雨、径流序列的小波方差分析反映了其波动能量随着时间尺度的变化情况，以及小波变换系数在不同周期内方差贡献的相对极大值，据此，可以准确判断出降雨、径流变化过程中存在的主周期。新安江的年降雨、径流序列的周期振荡一致，第1主周期均为17 a，第2，3周期分别是7 a，3 a，可见，新安江流域径流的振荡周期是由降雨的振荡周期引起的。

图4 新安江流域降雨、径流小波分析实系数和小波方差

3.5 降雨、径流的相关分析

计算新安江流域年降雨、年径流序列的相关系数达0.962 4，月降雨、月径流序列的相关系数达0.951 3，新安江流域的径流变化与同期降雨高度相关。结合径流与降雨的时间序列变化、突变和周期等特征分析可知，新安江流域的径流与降雨变化有着极其吻合的同步变化特征，新安江流域属于雨水补给型河流，降雨是影响该流域径流最主要的因子。

将月降雨与月径流建立降雨径流的拟合分析，降雨与径流以线性形式拟合最好，回归方程如下：

Y=0.776X-26.672

(1)

式中：X为降雨量(mm)；Y为月径流(mm)，该回归方程的拟合优度即判定系数R2为0.905，说明该线性拟合模型的效果显著。

由此可以利用新安江流域月降雨量的预测来判断月径流量的丰枯变化，以2016年1—12月的月雨量作为试报变量，用该回归方程可计算出同期月径流的预测值(图5)，除7月份的绝对误差达27.6 mm，其余11个月的绝对误差均在10 mm以内；相对误差均20%以内；预测值与实况值较为一致。

在业务应用过程中，新安江流域特别需要有月尺度上的径流预测来作为生产和调度的决策依据，对于流域径流长期预报可以利用气象部门的长期预报中的流域月雨量预报，并通过该线性回归模型对流域月径流量做出较准确的预测。因此该线性回归模型在新安江流域长期径流量的预报业务中有一定的应用价值。

图5 月径流分月预测与实测对比

4 结 论

(1) 新安江流域降雨、径流的年际分配不均匀，但两者演变趋势具有一致性，均呈现波动增多的趋势。两序列的月际分配均呈一致的单峰型变化，春夏季的4—7月是该流域降雨和径流最为集中的时段，也是同期洪水灾害最为多发的时段。

(2) 新安江流域径流与降雨的突变异常时间基本吻合，径流变化受降雨的影响显著。新安江流域年雨量序列有4个明显的突变年份，其中1971年、1988年、2009年是降雨明显增多的突变年份，2003年是径流明显减小的突变年份。径流量序列有3个发生突变的年份，其中1971年、2009年的年径流由枯转丰，2003年的年径流由丰转枯。

(3) 新安江的年降雨、径流序列的周期振荡一致，第1主周期均为17 a，第2，3周期分别是7 a，3 a。目前正处于不同周期中的丰水年中。新安江流域径流的振荡周期是由降雨的振荡周期引起的。

(4) 新安江流域的径流与降雨变化有着极其吻合的同步变化特征，新安江流域属于雨水补给型河流，降雨是影响该流域径流最主要的因子。

(5) 新安江流域降雨与径流高度相关。利用这两组序列建立线性回归模型效果较好，在流域长期径流量预报业务中有一定的应用价值。