二维haar小波matlab_洪泽湖入湖水沙序列的多时间尺度小波分析

摘要：

为了分析洪泽湖入湖水沙序列的多时间尺度变化过程及特征,预测洪泽湖入湖水沙变化趋势,利用洪泽湖1975—2015年的年来水、来沙量数据,在统计分析的基础上,对洪泽湖年来水、来沙序列进行延伸和距平处理,采用小波分析方法对处理后的洪泽湖入湖水沙距平序列进行Morlet连续复小波变换,得到了小波变换系数的实部等值线图和模方等值线图,小波方差图和主周期变化趋势图。结果表明:洪泽湖近41 a年来水、来沙序列存在着多时间尺度特征,在不同的时间尺度下,洪泽湖年来水、来沙量丰枯变化趋势不同;年来水、来沙序列存在着同样的第一主周期,为30a特征尺度,在该特征尺度下,预测2015年之后未来短时间内洪泽湖年来水量、来沙量将分别进入到枯水期和少沙期且都将处于相应的峰值期。

关键词：

洪泽湖; 入湖水沙; 小波分析; 多时间尺度; 趋势预测;

作者简介：

刘源(1995—),男,硕士研究生,研究方向为河流动力学及河床演变。E-mail:lyuantju@163.com;

徐国宾(1956—),男,教授,博士研究生导师,博士,研究方向为河流动力学及河床演变。

E-mail:xuguob@sina.com;

基金：

国家重点研发计划资助项目(2017YFC0405602);

引用：

刘源，徐国宾，段宇，等. 洪泽湖入湖水沙序列的多时间尺度小波分析[J]. 水利水电技术，2020,51(2)：128-135.

LIU Yuan, XU Guobin, DUAN Yu, et al. Multi-time scale wavelet analysis of streamflow and sediment sequences into Hongze Lake[J]. Water Resources and Hydropower Engineering，2020,51(2)：128-135.

---

0 引 言

全球气候变暖大背景下,极端降水事件发生的频率及强度不断增大,直接导致径流、泥沙等水文学研究对象更加具有不确定性和复杂性,常规的分析方法很难对它们的发展规律和变化特征有清晰的了解和准确的预测。20世纪80年代兴起的小波分析方法则能很好地分析出各种水文序列的变化特征,它是一种窗口大小固定但形状可变的时频局部化分析方法,而其中的Morlet复值小波能将水文时间序列分解成交织在一起的多尺度成分,并能够研究对象的微小细节,进而从不同时间尺度分析出水文研究对象的变化规律。

多位学者应用小波分析方法对不同河流及湖泊的水沙特性进行了大量研究,薛小杰等运用小波分析方法对黄河上中下游四个水文站的实测年径流序列进行了研究,结果表明多分辨分析用于水文时间序列趋势分析是合理可行的;桑燕芳等提出小波互相关度定量指标,并对黄河利津站和花园口站1950—2003年年径流序列进行分析,得出时滞为0时序列在时间尺度为3a、7a、11a和20a时具有良好的互相关关系;马海波等应用小波分析方法中良好的时频分析功能,分析了鄱阳湖湖口出湖径流序列的多时间尺度变换过程,得出鄱阳湖出湖流量在长时间尺度上丰枯变化较慢,短时间尺度上丰枯变化明显;胡光伟等应用小波分析方法研究了湘江主要控制站1951—2011年长时间序列径流与泥沙资料,研究表明径流量和输沙量均存在20a和7a的周期性变化规律;郭彦等应用小波分析方法及相关理论,分析了黄河上游唐乃亥站1957—2009年径流与输沙序列,结果表明黄河上游径流与输沙序列存在着年际尺度和年代际尺度的周期变化规律。

洪泽湖是淮河流域重要的大型水库和航运枢纽,在南水北调东线工程发挥着十分重要的调蓄功能。近年来多位学者对洪泽湖的水文序列特性进行了大量的研究,洪国喜等对洪泽湖入、出湖水沙特性进行了研究,研究表明洪泽湖出湖沙量小于入湖沙量,整体呈淤积趋势,入湖泥沙的多少主要受淮河控制;楚恩国分析了洪泽湖流域的径流及泥沙特性,分析结果表明洪泽湖流域径流年内分配不均,年际丰枯变化大,多年平均入湖河川径流量为325.1亿m3(其中6—9月占66%),洪泽湖区的最大年径流量达800亿m3/s(1954年),最小年径流量仅为28亿m3/s,洪泽湖泥沙的主要来源为流域表面随径流夹带入湖,泥沙含量6—7月最多,最大月平均含量为0.4～0.5 kg/m3,汛后则明显减少;虞邦义等对洪泽湖1983—2005年水沙特性进行了分析并总结出了洪泽湖冲淤的变化趋势,分析表明洪泽湖23年内洪泽湖多年平均来沙653万t,出沙326万t,多年平均淤积327万t,淤积率50. 1%,淮河干流来水来沙占洪泽湖总来水来沙的80%以上;王苓如等分析了在气候变化条件下洪泽湖水资源的变化程度及演变趋势,结果表明洪泽湖流域径流年内分配受降水变化影响较大;邓恒等对洪泽湖洪水调蓄能力进行了研究,并提出了表示湖泊对洪水蓄滞作用的参数——调节系数。而采用小波分析方法对洪泽湖水沙序列的多时间尺度变化特征分析的研究相对匮乏,几乎尚未涉及。近年来,研究和处理洪泽湖与淮河的河湖关系成为学术界研究的热点问题之一,而准河与洪泽湖的河湖关系随着水沙条件的变化而改变,近年来准河的径流量变化较小,但淮河输沙量却锐减,这种水沙特性的变化对干流河道发育及河湖关系演变产生了巨大的影响。因此,定量识别洪泽湖流域水沙序列的特征及变化规律是进一步明确两者河湖关系的研究关键。

本文根据洪泽湖1975—2015年的年来水、来沙量数据,在统计分析的基础上,采用小波分析方法对洪泽湖近41年的年来水、来沙序列进行研究,分析洪泽湖入湖水沙序列的时空变化规律及特征,找出水沙序列存在的主要尺度,分析水沙序列的丰枯变化过程,并预测未来短时间内的水沙变化趋势。这对预测洪泽湖的年际淤积趋势,建立准河与洪泽湖之间的水沙过程与河床湖盆演变的内在联系,揭示当前河湖关系特征及其形成机理具有重要的研究意义。

1 研究区域概况

洪泽湖是我国五大淡水湖之一,地处淮河流域中下游结合处的苏北平原中部偏西, 位于江苏省淮安、宿迁两市境内,地理位置在北纬33°06′—33°40′和东径118°10′—118°52′之间,水域面积2 069 km2,是淮河流域最大的湖泊型水库,属于过水性湖泊。洪泽湖上游入湖河流主要有淮河、怀洪新河、徐洪河、新汴河、濉河、老濉河和池河,下游泄洪出湖河道主要有淮河入江水道、入海水道、分淮入沂、苏北灌溉总渠、废黄河(见图1),承接着上游和中游15.8万km2的来水,淮河多年平均来沙量占洪泽湖总来沙量的83.0%,多年平均入湖水量占洪泽湖入湖水量的86.5%。