一种水文时间序列小波互相关分析方法.pdf

1、10申请公布号CN102033851A43申请公布日20110427CN102033851ACN102033851A21申请号201010297456022申请日20100930G06F17/0020060171申请人南京大学地址210093江苏省南京市汉口路22号72发明人王栋吴吉春桑燕芳祝晓彬74专利代理机构江苏圣典律师事务所32237代理人贺翔54发明名称一种水文时间序列小波互相关分析方法57摘要本发明公开了一种水文时间序列小波互相关分析方法，首先选择小波函数和时间尺度范围，然后对待分析的水文时间序列进行连续小波变换分析；计算水文时间序列在不同时间尺度上及不同时滞下的小波互协方差；求解两

2、序列在不同时间尺度上及不同时滞下的小波互相关系数，计算小波互相关度以描述两序列在整体时间域上的互相关程度，绘制小波互相关系数等值线图，实现对时间序列之间互相关关系进行时频综合分析。实例分析结果显示了小波互相关分析方法的有效性和优越性，该方法能够分析和定量描述非平稳时间序列在特定时间尺度和指定时滞下的互相关关系，可克服传统互相关分析方法的局限，具有更好的灵活性和适用性。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书7页附图3页CN102033854A1/2页21一种水文时间序列小波互相关分析方法，其特征在于包括以下步骤1选择小波函数和时间尺度范围，然后对待分析

3、的水文时间序列进行连续小波变换分析；2计算水文时间序列在不同时间尺度上及不同时滞下的小波互协方差；3根据小波互协方差计算结果，求解两序列在不同时间尺度上及不同时滞下的小波互相关系数；4求得不同时间尺度上及不同时滞下的小波互相关系数之后，计算小波互相关度；以描述两序列在整体时间域上的互相关程度；5求得小波互相关系数和小波互相关度之后，绘制小波互相关系数等值线图，并通过详细分析小波互相关系数等值线图，掌握所研究序列之间由整体到局部的互相关关系；实现对时间序列之间互相关关系进行时频综合分析。2根据权利要求1所述的水文时间序列小波互相关分析方法，其特征在于步骤2根据连续小波变换分析结果存在有实部和模两

4、种不同的情况，分别定义了小波互协方差的求解方法或式中WCOVXYA，KEWXA，BWYA，BK3上式中，WXA，B和WYA，B分别表示分析序列X和Y时得到的连续小波变换系数，A表示时间尺度因子，B表示时间位置因子，K表示时滞，R和I分别表示求解括弧内复数的实部和虚部，E表示求解均值，WCXYA，K表示在时间尺度A和时滞K下求得的小波互协方差。3根据权利要求1或2所述的水文时间序列小波互相关分析方法，其特征在于步骤3根据小波互协方差的两种不同的情况，分别定义了小波互相关系数的求解方法或上式中，WRXYA，K表示在时间尺度A和时滞K下求得的小波互相关系数，表示求解绝对值，其余公式符号同上。4根据权

5、利要求1或2所述的水文时间序列小波互相关分析方法，其特征在于步骤4根据小波互相关系数求解结果，定义了小波互相关度的求解方法，以描述两序列在整体时间域上的互相关程度，步骤如下41在求得两时间序列在尺度A和时滞K下小波互相关系数WRXYA，K的基础上，通过积分求得两时间序列在时滞K下对应整体时间域上的小波互相关程度的总和WRXYKWRXYKWRXYA，K2DA642然后，求解各时间尺度A下的小波互相关系数WRXYA，K的权重系数FWRXYA，KWRXYA，K2/WRXYK7权利要求书CN102033851ACN102033854A2/2页343求解两时间序列在时滞K下的小波互相关度为WCCDXYK

6、WRXYA，KFWRXYA，KDA85根据权利要求1或2所述的水文时间序列小波互相关分析方法，其特征在于步骤5绘制小波互相关系数等值线图时，以横轴表示时滞K的取值，纵轴表示时间尺度A的取值，图中的某点数值表征了对应尺度A和时滞K下两序列的互相关系数的大小。6根据权利要求1或2所述的水文时间序列小波互相关分析方法，其特征在于步骤5根据小波互相关系数等值线图，对两时间序列之间的互相关关系进行详细分析，主要步骤如下61通过对小波互相关系数等值线图进行垂向截取，分析在固定时滞下，两序列在各时间尺度上互相关程度大小的变化情况；62通过对等值线图进行横向截取，分析在固定时间尺度上，两序列在各时滞下互相关程

7、度的变化情况；63通过分析各时间尺度上小波互相关系数值的正负性，掌握两序列在各时间尺度上互相关性的正负变化情况；64通过对比分析各时间尺度上小波互相关系数绝对值的大小，识别并提取出对应某个或若干个互相关性明显的时间尺度范围；65通过对比分析在各时滞下小波互相关系数值的大小，识别出两时间序列之间最显著的时间延迟关系。权利要求书CN102033851ACN102033854A1/7页4一种水文时间序列小波互相关分析方法技术领域0001本发明涉及一种水文时间序列分析方法，具体是一种水文时间序列小波互相关分析方法。背景技术0002小波分析方法WAVELETANALYSIS，WA具有对非平稳时间序列进行

8、时频综合分析的能力崔锦泰小波分析导论西安西安交通大学出版社，1995，因此适合于研究具有多时间尺度变化特性的复杂水文水资源系统王文圣，丁晶，李跃清水文小波分析北京化学工业出版社，2005；LABATDRECENTADVANCESINWAVELETANALYSESPART1AREVIEWOFCONCEPTSJOURNALOFHYDROLOGY，2005，314275288。随着理论研究的深入和解决实际问题的需要，WA在水文水资源学中的应用日益增多王文圣，丁晶，向红莲小波分析在水文学中的应用研究及展望水科学进展，2002，134515520。综合分析可以看出，目前WA主要用于揭示和描述水文系统自身

9、的内部结构和变化特性。而揭示各水文要素之间相互关系如降雨与径流序列、气象因子与降雨和径流序列、水位与流量等是认识水文循环过程和揭示水文演变机制的另一重要途径。传统互相关分析方法包括互谱分析方法由于存在以下主要缺陷丁晶，邓育仁随机水文学成都成都科技大学出版社，19881仅适用于平稳各态历经序列；2无法揭示序列在不同时间尺度范围内的互相关关系，使其在实际应用中有较大局限。小波互相关分析方法WAVELETCROSSCORRELATION，WCC能够实现对两非平稳时间序列在特定时间尺度和指定时滞下互相关关系的定量描述，具有更大的适用性和优越性，因此可以很好地用于研究和揭示各水文要素之间的相互关系。然而

10、，由于针对WCC的系统研究偏少，目前缺乏统一的求解公式和分析方法，且本身理论体系还不甚完善，仅在经济学、信号处理、临床医学等方面得到一定的应用，而在水文水资源学中的研究和应用非常少，且国内目前未见有相关报道。0003为此，本发明旨在探讨适合于研究水文水资源学问题的小波互相关分析理论和方法体系。首先经分析和整理，系统地介绍用于水文序列分析的基于连续小波变换CONTINUOUSWAVELETTRANSFORM，CWT的小波互相关分析方法；同时定义了基于CWT的小波互协方差和小波互相关度两个定量指标，用于描述两时间序列在整体时间域上的互相关关系；然后提出绘制小波互相关系数等值线图的方法，通过该等值线

11、图可达到对两时间序列互相关关系进行“时频综合分析”的目的。最后结合具体实例加以简要分析，以显示小波互相关分析方法用于研究水文时间序列的适用性和优越性。0004分析实测时间序列XT和YT之间的互相关关系时，实际中常用下式求解互相关系数。00050006式中，X、Y分别表示序列XT和YT的均方差，分别表示XT和YT说明书CN102033851ACN102033854A2/7页5的均值，CXYK表示两时间序列在时滞K下的互协方差，RXYK为两时间序列在时滞也称时移K下的互相关系数，N表示时间序列的长度。0007实际中，当需要定量描述各非平稳时间序列之间在某特定时间尺度上的互相关关系时，该式无能为力。

12、发明内容0008本发明针对常规互相关分析方法存在的缺陷，并基于小波分析方法，提出了一种水文时间序列小波互相关分析方法，并绘制了小波互相关系数等值线图，通过该等值线图可达到对两时间序列互相关关系进行“时频综合分析”的目的。0009本发明所述的一种水文时间序列小波互相关分析方法，其特征在于包括以下步骤00101选择小波函数和时间尺度范围，然后对待分析的水文时间序列进行连续小波变换分析CONTINUOUSWAVELETTRANSFORM，CWT；00112计算水文时间序列在不同时间尺度上及不同时滞下的小波互协方差；00123根据小波互协方差计算结果，求解两序列在不同时间尺度上及不同时滞下的小波互相关

13、系数；00134求得不同时间尺度上及不同时滞下的小波互相关系数之后，计算小波互相关度；以描述两序列在整体时间域上的互相关程度；00145求得小波互相关系数和小波互相关度之后，绘制小波互相关系数等值线图，并通过详细分析小波互相关系数等值线图，掌握所研究序列之间由整体到局部的互相关关系；实现对时间序列之间互相关关系进行时频综合分析。0015上述步骤2根据连续小波变换分析结果存在有实部和模两种不同的情况，分别定义了小波互协方差的求解方法00160017或0018式中WCOVXYA，KEWXA，BWYA，BK30019上式中，WXA，B和WYA，B分别表示分析序列X和Y时得到的连续小波变换系数，A表示

14、时间尺度因子，B表示时间位置因子，K表示时滞。R和I分别表示求解括弧内复数的实部和虚部，E表示求解均值。WCXYA，K表示在时间尺度A和时滞K下求得的小波互协方差。0020上述步骤3根据小波互协方差的两种不同的情况，分别定义了小波互相关系数的求解方法00210022或0023上式中，WRXYA，K表示在时间尺度A和时滞K下求得的小波互相关系数，表示说明书CN102033851ACN102033854A3/7页6求解绝对值，其余公式符号同上。0024上述步骤4根据小波互相关系数求解结果，定义了小波互相关度的求解方法，以描述两序列在整体时间域上的互相关程度，步骤如下002541在求得两时间序列在尺

15、度A和时滞K下小波互相关系数WRXYA，K的基础上，通过积分求得两时间序列在时滞K下对应整体时间域上的小波互相关程度的总和0026WRXYKWRXYA，K2DA6002742然后，求解各时间尺度A下的小波互相关系数WRXYA，K的权重系数0028FWRXYA，KWRXYA，K2/WRXYK7002943求解两时间序列在时滞K下的小波互相关度WAVELETCROSSCORRELATIONDEGREE，WCCD为0030WCCDXYKWRXYA，KFWRXYA，KDA80031步骤5绘制小波互相关系数等值线图时，以横轴表示时滞K的取值，纵轴表示时间尺度A的取值，图中的某点数值表征了对应尺度A和时滞

16、K下两序列的互相关系数的大小。0032步骤5绘制小波互相关系数等值线图之后，对两时间序列之间的互相关关系进行详细分析，主要步骤如下003351通过对小波互相关系数等值线图进行垂向截取，分析在固定时滞下，两序列在各时间尺度上互相关程度大小的变化情况；003452通过对等值线图进行横向截取，分析在固定时间尺度上，两序列在各时滞下互相关程度的变化情况；003553通过分析各时间尺度上小波互相关系数值的正负性，掌握两序列在各时间尺度上互相关性的正负变化情况；003654通过对比分析各时间尺度上小波互相关系数绝对值的大小，识别并提取出对应某个或若干个互相关性明显的时间尺度范围；003755通过对比分析在

17、各时滞下小波互相关系数值的大小，识别出两时间序列之间最显著的时间延迟关系。0038本发明提供了一种适合于研究水文水资源学问题的小波互相关分析理论和方法体系。首先经分析和整理，系统地提出了用于水文序列分析的基于连续小波变换CWT的小波互相关分析方法；同时定义了基于CWT的小波互协方差和小波互相关度两个定量指标，用于描述两时间序列在整体时间域上的互相关关系；并提出绘制小波互相关系数等值线图的方法，通过该等值线图可达到对两时间序列互相关关系进行“时频综合分析”的目的。实例分析结果显示了小波互相关分析方法的有效性和优越性，应用其可实现由整体到局部认识时间序列互相关关系的目的。小波互相关分析方法能够分析

18、和定量描述非平稳时间序列在特定时间尺度和指定时滞下的互相关关系，可克服传统互相关分析方法的局限，具有更好的灵活性和适用性。因此，通过应用小波互相关分析方法对各水文要素进行分析和描述，更有助于深入认识复杂的水文循环过程和揭示水文演变机制，同时也可促进小波互相关分析自身理论和方法体系的不断完善。说明书CN102033851ACN102033854A4/7页7附图说明0039图1小波互相关系数等值线图绘制流程。0040图2利津站和花园口站年径流序列小波互相关系数等值线图图2A、不同时滞下图2B和不同时间尺度上图2C小波互相关系数曲线。0041图3利津站和花园口站年径流序列小波互相关度曲线和互相关系数

19、曲线。具体实施方式0042以下具体介绍基于连续小波变换的小波互相关分析方法0043令L2R表示定义在实轴上、可测的平方可积函数空间，则对于时间序列XTL2R，其连续小波变换CONTINUOUSWAVELETTRANSFORM，CWT公式为0044A，BR，A090045式中，T为T的复共轭函数，WXA，B为XT经连续小波变换后的小波系数。A为时间尺度因子，B为时间位置因子。A，BT为小波函数，是由一满足“容许性”条件式10的小波母函数T经尺度伸缩和平移后得到。式10中，表示复共轭函数T在频率处的FOURIER变换。00460047小波变换的实质是采用一种大小可变，位置可动的变窗对时间序列进行谱

20、分析，因此可满足序列时频局部化分析的要求。对水文序列进行小波变换后得到的小波系数是序列在不同时间尺度和不同时间位置上的投影，可用来刻画和描述水文序列的组成结构和多时间尺度变化特性。此外，通过研究两时间序列对应小波系数之间的关系，同样可以刻画和描述两序列之间的互相关关系。00481、基于CWT的小波互相关系数求解方法0049依据传统时间序列互协方差的定义式，在此处定义小波互协方差WAVELETCROSSCOVARIANCE。由于某些连续小波变换系数结果例如用“MORLET”小波函数进行小波变换有实部和模两个重要变量，实部表示不同特征时间尺度信号在不同时间位置上的分布和相位两方面的信息，模的大小主

21、要表示特征时间尺度信号的强弱。因此，此处分别定义对应的小波互协方差WCXYA，K00500051或0052式中WCOVXYA，KEWXA，BWYA，BK30053其中，WXA，B和WYA，B分别表示时间序列XT和YT在尺度A下对应的连续小波变换系数，K表示时滞，R和I分别代表括号内变量的实部REALPART和虚部IMAGINARYPART，E表示方括弧内结果的均值，WCXYA，K表示序列XT和YT在尺度A和时滞K下对应的小波互协方差。0054在求得小波互协方差的基础之上，首先介绍小波局部互相关系数WAVELETLOCALCORRELATIONCOEFFICIENT，WLCC，用于分析两序列在特

22、定时间尺度和时间位置点上的互说明书CN102033851ACN102033854A5/7页8相关关系。00550056然后对时间序列XT和YT进行小波互相关分析WAVELETCROSSCORRELATION，WCC。此处同样针对连续小波变换系数结果的实部和模两个变量，分别定义对应的小波互相关系数WRXYA，K00570058或0059WRXYA，K定量描述了时间序列XT和YT在时间尺度A上和时滞K下相应的互相关程度。00602、时间序列小波互相关度0061WCC方法主要用于分析两时间序列在特定时间尺度上和指定时滞下的互相关关系。为便于分析和描述时间序列之间在整体时间域上的互相关程度，笔者经分析

23、探讨，在此处提出基于CWT的时间序列小波互相关度WAVELETCROSSCORRELATIONDEGREE，WCCD的概念，用于刻画时间序列之间整体时间域上的互相关程度。具体如下0062在求得两时间序列在尺度A和时滞K下小波互相关系数WRXYA，K的基础上，通过积分求得两时间序列在时滞K下对应整体时间域上的小波互相关程度的总和0063WRXYKWRXYA，K2DA60064然后，各时间尺度A下的小波互相关系数WRXYA，K的权重系数可定义为0065FWRXYA，KWRXYA，K2/WRXYK70066定义两时间序列的小波互相关度WAVELETCROSSCORRELATIONDEGREE，WCC

24、D为0067WCCDXYKWRXYA，KFWRXYA，KDA80068WCCDXYK表征了两时间序列在时滞K下整体时间域上互相关程度的大小，其实质是在相同时滞K下，对各时间尺度上的互相关程度求解其加权期望值，因此可综合反映两时间序列在各时间尺度上关于互相关性的大小和分布两方面信息。0069此外，通过绘制WCCDXYK随时滞K的变化曲线WCCDXYKK曲线，可以分析两序列在整体时间域上随时滞变化时互相关性的变化特点和规律。00703、小波互相关系数等值线图0071类似于通过绘制小波系数等值线图可以分析时间序列自身的结构和多时间尺度变化特性，在此处提出绘制小波互相关系数等值线图的方法，并利用小波互

25、相关系数等值线图由整体到局部定量分析时间序列之间在各时间尺度和各时滞下的互相关关系，即达到对时间序列之间互相关关系进行“时频综合分析”的目的。0072小波互相关系数等值线图的绘制方法图1如下1首先选择合理的小波函数和时间尺度范围，对所研究各时间序列进行连续小波变换，得到各序列对应的小波系数结果；2利用式4或式5求解两时间序列在某时间尺度A和时滞K下对应的小波互相说明书CN102033851ACN102033854A6/7页9关系数；3依次取不同的A和K值，分别对两序列进行小波互相关分析，最终得到两时间序列在各尺度和各时滞上对应的小波互相关系数；4绘制小波互相关系数等值线图，其中横轴表示时滞K的

26、取值，纵轴表示时间尺度A的取值，图中的某点数值表征了对应尺度A和时滞K下两序列的互相关系数的大小；5通过定量分析小波互相关系数等值线图，掌握所研究序列之间由整体到局部的互相关关系；6通过求解WCCD，掌握两序列不同时滞下在整体时间域上的互相关程度。0073小波互相关系数等值线图对描述和刻画两时间序列之间的相互关系具有重要的应用意义。主要如下1通过对小波互相关系数等值线图进行垂向截取，可以分析在固定时滞下，两序列在各时间尺度上互相关程度大小的变化情况；2通过对等值线图进行横向截取，可以分析在固定时间尺度上，两序列在各时滞下互相关程度的变化情况；3通过分析各时间尺度上小波互相关系数值的正负性，可以

27、掌握两序列在各时间尺度上互相关性的正负变化情况；4通过对比分析各时间尺度上小波互相关系数绝对值的大小，可以识别并提取出对应某个或若干个互相关性明显的时间尺度范围；5通过对比分析在各时滞下小波互相关系数值的大小，可以识别出两时间序列之间最显著的时间延迟关系等。0074实例分析0075选择黄河花园口站和利津站两水文站点实测的54年19502003年径流序列为例进行分析。研究两站年径流序列之间的相互关系，对认识黄河下游的径流过程和水文变化特性具有重要的意义。0076应用前述的基于CWT的小波互相关分析方法，深入分析在不同时间尺度上以及不同时滞下两年径流序列之间的互相关关系。此处选用式4求解小波互相关

28、系数WRXYA，K。结合文献桑燕芳，王栋连续小波变换在黄河河口地区特性分析中的应用研究A第五届中国水论坛论文集C北京中国水利水电出版社，2007，766770中的相关分析结果，此处同样选用“MORLET”小波函数对花园口54年年径流序列进行连续小波变换。然后依据CWT分析结果，分别求解各时间尺度和各时滞下的小波互相关系数，并绘制小波互相关系数等值线图图2A。其中，最大时间尺度取50A，最大时滞取20A。为便于分析说明，分别绘制时滞K0、2和5时的WRXYA，KA曲线图2B，和时间尺度A3、7、11和20时的WRXYA，KK曲线图2C。最后绘制小波互相关度WCCDXYKK曲线和互相关系数曲线图3

29、。传统互相关分析方法计算结果显示图3两年径流序列之间的最大互相关系数为RXY009655。0077综合上述两径流序列互相关关系的分析结果，可得到如下主要结论00781对小波互相关系数等值线图进行分析，可以识别出两径流序列主要在四个时间尺度上存在较为明显的互相关关系A3、7、11和20。此四个时间尺度对应的也正好是各径流序列自身四个明显的周期值。00792对小波互相关系数等值线图进行垂向截取，可以了解两序列在固定时滞下各时间尺度上的互相关关系。图2B显示，在同一时滞下，不同时间尺度上两序列互相关性程度的大小不同；在不同时间尺度范围内两序列互相关的正负特性也不相同。以K5为例，在时间尺度4A、8A

30、、和50A上的小波互相关系数分别为058、095和088；随着时滞K由0增大到5，在时间尺度550A范围内两序列互相关关系有逐渐减弱的趋势，而在15A时间尺度范围内时滞5时的互相关关系要较时滞2的互相关关系更明显。说明书CN102033851ACN102033854A7/7页1000803对小波互相关系数等值线图进行横向截取，可以了解两序列在固定时间尺度上各时滞下的互相关关系。图2C显示，在相同时间尺度上，两径流序列互相关性的大小随时滞变化时存在波动，且存在相关性逐渐减弱的趋势。以A7为例，两序列在时滞K0、4、9、14和19时，互相关系数的波动变化为094、088、078、059、058；随

31、着时间尺度的逐渐增大，两径流序列互相关系数随时滞变化的波动周期也逐渐增大。00814小波互相关度分析结果图3显示，两序列整体时间域上在不同时滞下的互相关程度不同。小波互相关度在时滞K0时取得最大值096，与传统互相关分析方法的结果一致。之后随时滞增大时小波互相关度逐渐减小，在K13时小波互相关度又达到另一极大值080。与互相关系数曲线对比可以看出，小波互相关度曲线更加光滑、波动规律更加明显，这主要是由于小波互相关度曲线是上述各时间尺度上互相关关系的加权平均值，因此可以更有效地描述两序列在整体时间域上随时滞推移时的互相关关系变化情况。而传统互相关分析方法无法考虑和定量描述不同时间尺度上序列互相关

32、关系之间的差异，且曲线中得到的多个峰值点缺乏实际的物理意义。00825综合图2和图3的分析结果可以看出利津站和花园口站两个实测年径流序列在时滞为0时的互相关关系最明显，这与径流由花园口站至利津站的实际时间相符合；两年径流序列在时间尺度为3A、7A、11A和20A时的互相关关系最明显，即表明在序列自身的四个明显周期变化时间尺度范围内，两序列的互相关关系也最明显。6由实例分析可以看出，小波互相关分析方法可以揭示和描述非平稳时间序列在不同时间尺度上和不同时滞下的互相关关系，因此可克服传统互相关分析方法的局限，有助于对时间序列之间互相关关系进行全面细致的定量分析。说明书CN102033851ACN102033854A1/3页11图1图2A说明书附图CN102033851ACN102033854A2/3页12图2B图2C说明书附图CN102033851ACN102033854A3/3页13图3说明书附图CN102033851A