基于小波分析的淮河干流上游年径流量周期分析研究

　　【摘 要】随着淮河流域区域经济的发展，水资源供需矛盾日益加剧。而淮河干流上游对整个淮河流域起关键性作用。因此，研究淮河干流上游的径流量变化规律，对整个淮河流域的水资源规划管理、开发利用具有一定的意义。

　　【关键词】小波分析 小波方差 时频变化 年径流序列

　　1 引言

　　淮河发源于河南省桐柏山，自西向东流经鄂、豫、皖、苏四省，主流在三江营入长江，全长约 1000km，总落差 200m。从河源到洪河口为上游，流域面积 3.06 万 km2 ，河长 360km，多年平均径流量为 92.48 亿 m(3 1956~2010 年平均值，下同)，径流量年内分配不均，年际变化较大。王家坝水文站是淮河干流上游重要的控制站，其 1956~2010 年径流量过程线见图 1。

　　河川径流序列是一个具有多时间尺度特征的复杂过程，径流的变化规律和丰枯变化趋势在不同的时间尺度下是不相同的，即多时间尺度。所谓多时间尺度变化，其含义就是指研究系统的变化并不存在一种真正传统意义上的周期性，即没有非常明显或者严格的交替循环过程，但是又存在时而以这种周期变化，时而以另一种周期变化，而且在相同的一个时段内，研究系统又同时包括各种时间尺度上的周期变化的现象，即这种肉眼无法观察或从直观数据中无法得到的多时间尺度的周期性和多层次时间尺度结构。所以在刻画径流序列中，水文要素的多时间尺度分析是无法避免的。而小波分析的优势就是在时域和频域上均具有良好的局部化特征和多分辨功能，对于多时间尺度的径流序列，不仅可以对其进行局部化分析，还可以分析其内部精细的结构特征，得到其在不同时间尺度下的周期性和演变情况。

　　2 小波分析

　　2.1 小波函数

　　小波函数指的是具有震荡特性，在有限的区域内能够迅速衰减到 0 的一类函数 (t)。

　　3 淮河干流上游年径流序列的小波分析

　　用上述小波分析的方法，对延伸后标准化的年径流量序列进行 Morlet 小波变换，用 Matlab 绘制出小波系数实部、模的等值线图和小波方差图。

　　3.1 年径流序列小波变化模时频变化分析

　　从图 2 中可以看出，研究地区不同时段各时间尺度的震动能量密度时频强弱分布。在 25a 左右的时间尺度上其能量密度比较强，在时域上分布也比较明显，主要发生在 1985～2010 年之间，震动中心在 2005 年左右;在 10～15a 的时间尺度上，出现了两个震动中心，一个是全局最强震动能量密度，发生在 1956～1970 年，震动中心在 1956 年，说明该时间尺度上变化周期最明显。另一个在发生在 1983～ 1996 年，震动中心在 1990 年左右;在 5a 左右的时间尺度上，也出现震动能量密度较强的中心，发生在 1957 年左右;在 1～5a 的时间尺度上，震动能量较强信号发生在 1998～ 2006 年，震动中心发生在 2003 年左右。 3.2 年径流序列小波变化实部时频变化分析。

　　图 3 是研究区域 1956～2010 年径流量序列小波系数实部的时频分布图。图中清晰地显示了研究区域径流量序列不同时间尺度变化、突变点分布及位相结构。在时间尺度 a 相同的情况下，正的小波系数与丰水期相对应;负的小波系数与枯水期相对应;小波系数为零与突变点相对应。可以看出该年径流量序列存在着 0～3a、4～8a、9～15a、16～32a 左右的四类时间尺度的变化规律。其中 9～15a 和 16～32a 左右的丰枯交替变化最为明显。对于 16～32a 的时间尺度来说，年径流经历了枯→丰→枯→丰→枯→丰→枯 7 个交替变化。具体枯水年时段为：1956～1965 年、1975～1982 年、 1991～1999 年、2007～2010 年，且 2010 年其实部等值线未闭合，说明在 2010 年以后可能还会存在一段枯水期。丰水年时段为：1966～1974 年、1983～1990 年、2000～2006 年。

　　而对于 9～15a 的时间尺度来说其年径流序列经历了丰→枯→丰→枯→丰→枯→丰→枯→丰→枯→丰→枯→ 丰→枯→丰→枯 16 个交替变化。枯水年时段具体表现为： 1958 ～1962 年 、1967 ～1968 年 、1972 ～1974 年 、1979 ～ 1981 年、1986～1989 年、1994～1996 年、2001～2004 年、 2008～2010 年。丰水年时段为：1956～1957 年、1963～1966 年、1669～1971 年、1975～1978 年、1982～1985 年、1990～ 1993 年、1997～2000 年、2005～2007 年。

　　对于 4～8a 的时间尺度，周期变化显示较弱，丰枯周期变化竟有 13 次之多，位相结构不稳定。由此得知该时间尺度下年径流变化有明显的突变特征。同时，对于 3a 左右的时间尺度，丰枯周期变化更是达到了 14 次，可知在 3a 左右的时间尺度下，其年径流突变性较强，周期性较弱。

　　3.3 年径流序列小波变化方差分析

　　图 4 为王家坝站 1956～2010 年径流量序列小波系数方差图。小波方差图能够反映年径流时间序列的波动幅度随时间尺度 a 的分布情况，可以用来辨识时间序列中各种尺度的扰动强弱和周期变化特征。由此可以确定年径流变化的主要周期。由图可知：小波方差出现 4 个较为明显的峰值，依次为 3a、6a、11a 和 25a 的时间尺度。其中，11a 的时间尺度对应着最大峰值，所以 11a 左右的时间尺度对应的周期变化最强烈，是其年径流变化的第一主周期。25a 左右的时间尺度对应着第二峰值，是其年径流变化的第二主周期; 6a 左右的时间尺度对应着第三峰值，是其年径流变化的第三主周期;而 3a 左右的时间尺度对应着最小的峰值，是其年径流变化的第四主周期，相对于前三个主周期第四主周期的周期变化较弱。

　　绘制第一、第二、第三和第四主周期所对应的小波变换系数实部过程线，如图 5 所示。

　　将年径流序列小波变化方差分析与年径流序列小波变化实部时频变化分析结合，可以看出：在王家坝站年径流变化的四个主周期和实部分析中，四类时间尺度的变化规律相对应，并且四个主周期分别是是其时间尺度中心。

　　4 结论

　　通过Morlet 小波分析法对王家坝站 1956～2010 年共计 55 年的年径流序列进行小波变化，结果表明：王家坝站存在 0～3a、4～8a、9～15a、16～32a 左右的四类时间尺度的变化规律，其时间尺度中心分别对应为 3a、6a、11a、25a。根据小波方差图可知，11a 的时间尺度为年径流变化的第一主周期，25 年的时间尺度为其第二主周期，6a 的时间尺度为其第三主周期，3a 的时间尺度为其第四主周期。

　　推荐阅读：长洲水利枢纽船闸上游引航道潜坝基础处理设计

　　摘要：本文主要介绍长洲水利枢纽船闸上游引航道潜坝基础处理设计，即在水库内堆填造陆，采用水泥土搅拌桩进行混凝土潜坝基础加固的方法。

《基于小波分析的淮河干流上游年径流量周期分析研究》

上一篇：法医昆虫标本的数字库建立与应用

下一篇：浅谈背景技术及技术效果的撰写