小波分析属于时频分析,
它不同于Fourier分析,
Fourier分析使用的是一种全局的变换,
无法表述信号的时频局域性质.
小波变换具有多分辨率分析的特点,
在时频域都具有表征信号局域特征的能力[1].
小波分析理论已广泛地应用到各个学科领域,
本文应用小波分析理论,
对黄河上、中、下游河川径流长系列进行深入细致的分析,
以揭示黄河全流域河川径流的变化特性,
使我们对黄河径流变化有更加深刻的认识,
为黄河流域规划、水资源优化配置等提供更加全面的技术支持.
2 黄河上中下游径流变化规律
选取兰州站、三门峡站和花园口站为黄河上中下游的代表站,
所用资料为1919~1997年共79年天然河川径流,
3站多年平均径流量分别为3.2755×1010,5.0153×1010和5.5852×1010m3.
为减小资料两端的影响,
采用数据反摺的方法,
外延资料两端的数据.
采用Morlet复值小波对3站径流序列进行连续小波变换,
可以得到小波变换系数的实部、虚部、模、模平方、相位等信息(具体见图1~3),
通过分析这些信息,
能够揭示河川径流变化的多时间尺度结构.
2.1 黄河上游小波变换结果及分析
(1)小波系数实部
小波系数的实部包含给定时间和尺度下,
相对于其他时间和尺度,
信号的强度和位相两方面的信息[4,5].
小波系数实部为正时,表示径流量偏多,
图中用实线绘出,
“+”表示正值中心;
为负时表示径流量偏少,
图中用虚线绘出,
“-”表示负值中心.
从小波系数的实部可以看出不同尺度下的丰枯位相结构,
表明不同的时间尺度所对应的径流丰枯变化是不同的.
小尺度的丰枯变化则表现为嵌套在较大尺度下的较为复杂的丰枯结构.
从图1(a)中可分析出径流存在明显的年际变化和年代际变化,
从上至下分析得出径流存在30年以上,10~30,7~17及10年以下4类尺度的周期变化规律.
从较大尺度30~40年分析,
径流变化出现丰枯交替的准2次振荡,
对应于这种大尺度的丰枯交替,
黄河上游径流变化表现出了明显的突变特征,
具体表现为1929年以前偏枯,
1930~1949年偏丰,
1950~1969年偏枯,
1970~1988年偏丰,
1988~1997年偏枯,
直到1997年等值线仍未闭合,1997年后一段时间仍处于偏枯期.
从10~30年尺度分析,径流出现准3次振荡,具体表现为1924年以前偏丰,1925~1935年偏枯,1936~1947年偏丰,1948~1958年偏枯,1959~1969年偏丰,1970~1980年偏枯,1981~1991年偏丰,1992年以后偏枯,而且这一偏枯曲线直到1997年仍未闭合,与30~40年尺度的变化一样,1997年以后一段时间仍将处于偏枯期.
30~40和10~30年尺度的周期变化在整个分析时段表现得非常稳定,具有全域性.
07~17年尺度的周期变化,主要在30~80年代表现活跃,存在负、正、负的周期振荡.
大尺度下的一个丰期(或枯期),包含小尺度下的若干个丰枯期.
对应于10年以下较小尺度的径流变化,
较为明显的周期为1~5年尺度的年际变化,
主要发生于20世纪40~70年代中期,
70年代中期之后,
周期中心上移,
变为5~10年左右尺度的周期变化,
到80年代中期,
周期中心又下移到1~5年尺度,
其他时间尺度的周期变化不明显.
对于小尺度而言,径流突变点增多.
不同时间尺度下,径流突变点个数、时间位置都有所不同,径流突变应针对具体时间尺度来讨论.
图1 三州站径流小波变换等值线图(a)实部;(b)模;(c)模平方
从图的上部可看到存在30年以上的周期变化,随着径流资料序列的延长,可以分析出更大尺度的周期变化.以22年尺度为例,各丰枯期的径流值与其小波系数的对照见表1,表中数据进一步表明小波系数的实部可以表征实际径流的丰枯位相变化,两者具有对应关系.
表1 丰枯期径流与小波系数实部对应表(单位:108m3)
注:表中时段19~24代表1919~1924年,以此类推
(2)小波系数模
Morlet小波变换系数的模值表示能量密度,
模值图把各种时间尺度的周期变化在时间域中的分布情况展示出来,小波变换系数的模值越大,表明其所对应的时段和尺度的周期性越明显.从图1(b)中可看出30~40和10~30年尺度的小波变换模值较大,说明这两个周期变化最明显;其次是7~17和10年以下尺度的周期变化.
(3)小波系数模平方
小波系数的模平方相当于小波能量谱,
所以从小波系数模平方图可分析出不同周期的振荡能量,
从图1(c)中可看到,10~30和30~40年尺度的周期非常显著,占据整个时段,能量最强,10年以下尺度的周期也比较显著,但10年尺度的周期变化具有局部化的特征.其他尺度的周期变化都较弱,能量较低.
2.2 黄河中游径流变化规律
(1)小波系数实部图2 (a)从上至下可分析得出代表黄河中游的三门峡站河川径流序列的变化具有30~40,10~30,7~17及10年尺度以下的周期变化.周期从大到小呈现嵌套结构,各尺度下的周期变化特征与上游相似.
(2)小波系数模从图2 (b)分析,模值最大的是10~30年尺度的周期,表明这一尺度的周期变化最明显,其次是30~40和10年以下尺度的周期.
(3)小波系数模平方从图2(c)可知,能量最强的是全域性的10~30年尺度的周期变化,呈条状分布在整个分析时段;30~40年尺度的周期变化,在80年代以后能量增强,与10~30年尺度的周期能量相当
2.3 黄河下游径流变化规律
(1)小波系数实部图3(a)从上至下,代表黄河下游的花园口站河川径流序列的变化具有3 0~40,10~30,7~17及10年尺度以下的周期变化,各尺度下的周期变化特征与上中游相似.
图2 三门峡站径流小波变换等值线图(a)实部;(b)模;(c)模平方
(2)小波系数模分析图3(b)可知,30~40,10~30和10年以下尺度的周期均具有全域性,其中10~30年尺度的模值最大,周期变化最明显;30~40年尺度的周期在80年代以前表现明显,10年以下尺度的模值与30~40年的相同.
(3)小波系数模平方从图3(c)中看到,能量最强的是全域性的10~30年的周期,呈条状占据整个分析时段,它对原径流时间序列的方差贡献最大.
图3 花园口站径流小波变换等值线图(a)实部;(b)模;(c)模平方
3 小波分析检验以上分析出黄河上中下游径流序列存在几个周期范围,如何判断哪个周期占据主要地位,即对径流序列的变化起主要作用呢?下面通过小波方差进行检验[5],小波方差的计算公式为
3 小波分析检验
以上分析出黄河上中下游径流序列存在几个周期范围,如何判断哪个周期占据主要地位,即对径流序列的变化起主要作用呢?下面通过小波方差进行检验[5],小波方差的计算公式为0
式中Wf(a)为小波方差,Wf(a,b)为小波系数.0利用(5)式,计算各时间尺度对应的小波方差,据此可以确定径流时间序列中存在的主要周期,上中下游小波方差图分别见图4.
图4 小波方差图4(a)中有5个峰值,分别对应上游4,9,15,22和38年的时间尺度,第一峰值是22年尺度,说明22年左右的周期振荡最强,为上游径流变化的第1主周期,第2~5主周期依次为38,15,4和9年.
图4(b)显示中游小波方差有4个明显的波峰,根据峰值大小依次是22,16,38,4和9年,中游径流变化第1主周期与兰州站一致是22年,第2~5主周期与兰州站稍有不同.
图4(c)表明下游小波方差有4个峰值,依据方差大小判断,下游径流变化的第1主周期是22年,其次是38,9和4年.
4 黄河径流变化趋势分析
根据小波方差检验出的主要周期,可绘出主要周期的小波系数图,
图5为22年尺度周期的小波系数图,
因为22尺度的周期变化能量最大,
对原始径流序列的方差贡献也最大,
所以黄河上中下游1997年后5年左右,
低频波仍将发展,
即径流大约到2002年前后仍将于偏枯阶段,
大约2003~2012年上游径流将处于偏丰期.
图5 22年尺度径流变化趋势
5 结论
通过以上分析,得出以下结论:
(1)黄河上、中、下游河川径流变化存在多时间尺度的特征,具有30年以上尺度,10~30,7~17及10年以下尺度的周期变化,大尺度的周期变化嵌套着小尺度的周期变化.
(3)根据主要周期可以分析径流变化趋势,黄河上、中、下游径流推测大约在2003至2012年将处于偏丰期.
(4)径流变化在时域中存在多层次时间结构的局部化特征,径流变化周期具有局部化的特征,某一周期变化在某一时段表现明显,而在其他时段表现可能不明显,径流变化周期只在一定的时段内有意义.