花园口站年最大洪峰流量时间序列的HHT分析 Analysis on Annual Maximum Peak Flow Series at Huayuankou Station Based on Hilbert-Huang Transform.pdfVIP

第28卷第3期 地理与地理信息科学 V01．28No．3 andGeo-InformationScience 2012 2012年5月 Geography May 花园口站年最大洪峰流量时间序列的HHT分析 任健1，史红玲2 (L中国水利水电科学研究院，北京100048；2．国际泥沙研究培训中心，北京100048) 摘要：Hilbert-Huang变换能够定量描述非线性、非平稳复杂时间序列的时频特性，较传统分析方法更具优势。通 过对时间序列进行EMD分解。得到变化过程的内在模态函数和趋势项函数，而后对各内在模态函数进行Hilbert- Huang变换，从而揭示出时间序列的多时间尺度特征。以黄河花园口站1952--2009年的年最大洪峰流量时间序 列为例，对其进行多时间尺度分析．得到不同波动周期的振荡分量及趋势分量，具体分析了各分量的变化特征。结 a 果表明，花园口年最大洪峰流量变化过程中存在准3．2a、准6．4a、准11．8a和准31，0a周期的波动，其中准3．2 和准6．4 a的周期波动是引起原序列波动的主要原因。近60年来花园口年最大洪峰流量变化呈递减趋势，由此揭 示了年最大洪峰流量变化过程的多时间尺度特征。在此基础上，探讨．『各波动分量变化的影响因素．其变化与大 气低频振荡、ENS()、太阳活动及气候变迁等因素有关。 关键词：Hilbert-Huang变换；经验模态分解；多时问尺度；年最大洪峰流量 中图分类号：TVl42文献标识码：A 文章编号：1672--0504(2012)03--0083--04 黄河花园口站是黄河下游重要的控制性监测断 稳数据序列的方法，它将时间信号，(￡)分解成一系 Mode 面。其控制流域面积达73万kmz，约占黄河流域面列本征模态函数(Intrinsic 积的97％，分析花园口站的水沙变化情势是研究黄 个IMF分量都具有如下特征：从全局特性看，极值 河下游河道演变的前提和基础。长期以来，关于黄 点数必须和过0点数一致或至多相差一个；在某个 河下游水沙变化特征的研究很多D-43，但多依据传局部点，极大值包络和极小值包络在该点的值的算 统的数理统计或简单的数据平滑等方法，不能揭示 术平均和是0。 水沙变化的多时问尺度特性[5]。而Fourier分析和 EMD分解过程为：找出序列，(f)所有的极大值 小波分析等方法虽然可以进行频谱分析，但其方法 和极小值点，分别用三次样条函数拟合成上下包络 本身存在一些不足和缺点【6J]。希尔伯特一黄变线，得到平均包络线m，，将原序列减去m，得到去掉 换[6’7】(Hilbert—HuangTransfoKin，HHT)技术是近低频序列的新序列h。。一般h-并不能满足IMF的 年来发展的一种最新的处理非线性、非平稳时间序 两个特征要求，仍是非平稳的，因而可以多次重复上 列的方法，它直接由时间序列本身自适应地构造出 述过程，使平均包络线趋近于o，得到第一个IMF分 基函数，进而分离得到不同尺度的特征分量，而不用 量c。，c-代表原始序列中最高频的分量。即： 预先设定基底函数再展开，因而该法可更准确反映 n(幻一∥f)一旬 (1) 序列的内在规律特征。本文运用HHT技术，对黄 对n(￡)继续进行上述分解，直到剩余部分为单 一信号或其值小于预先给定的值，分解结束。原始 河下游花园口站1952—2009年的实测年最大洪峰 流量序列进行多时间尺度分析，诊断花园口站年最 时间序列，(￡)可以表示为： 大洪峰变化过程中蕴含的多时间尺度振荡结构和特