曹娥江低潮位序列突变成份分析.doc

《曹娥江低潮位序列突变成份分析.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、曹娥江低潮位序列突变成份分析摘要：本文提出了动态分维数的计算方法，并应用此方法对曹娥江低潮位序列中的突变成份进行了分析。结果表明，文中所用的方法对识别人为或自然因素引起的潮位时间序列中的突变成份是有效的，具有一定的实用价值。关键词：动态分维数时间序列河床演变1前言　　水文时间序列一般由确定成份和随机成份组成，暂态成份属确定的非周期成份，它包括趋势、跳跃、突变等成份，突变可看作是跟踪跳跃的一种特殊情况。潮位时间序列和其它水文要素时间序列一样，也往往是由两种或两种以上的成份合成的序列，如何识别序列中长时间的动态变化特征，识别其突变成份则可分析不同时期序列变

2、化程度，揭示其突变性，借以分析导致变化的过程和原因，是研究时间序列机制的重要内容之一。本文提出了动态分维数的计算方法，并将此方法应用于曹娥江低潮位序列分析，求得了序列的突变点，经合理性分析和成因分析表明，动态分维数对于进行序列分期，识别序列中的突变行为具有简单实用的特点，对于了解序列变化情况和分析导致该变化的过程、原因，进行预警分析具有实际意义。２突变成份分析　　随着时间的增长，对于潮位序列中各均值来说，或是增加或是减少，造成序列长期向上或向下缓慢地变动，这时序列的任何参数都将随着时间增长，呈现系统连续地递增或递减的变化，这种有一定规律的变化叫趋势。趋

3、势存在于序列的任何参数中，如均值、方差、自相关系数等，对于实际时间序列多在均值中寻找其趋势成份。而突变是指时间序列急剧变化的一种形式，即从某时间起序列前后时段总体水平发生急剧变化，从一种状态过渡到另一种状态（如陡升或陡降）。时间序列出现突变成份，其成因与趋势成份的出现基本一致，是由人为的或自然的原因引起。突变成份一般也出现在序列均值、方差等统计量中，实际上多在均值中寻找突变成份。识别序列中长时间的动态变化特征，识别其突变成份则可分析不同时期序列变化程度，揭示其突变性，借以分析导致变化的过程和原因，为河床演变分析提供依据。　　判断突变成份是否存在，传统的

4、作法是采用分割样本的方法进行检验，即先找出分割点，然后用有关方法检验分割点前后两个样本序列是否一致，以此判断时间序列中突变成份的存在与否。本文采用动态分维数来求序列的临界点（突变点），并以此来进行序列分期和统计分析。2.1动态分维数Dft10　　分形几何[1]指出,客体在单位时间内由于自身的发展而填充空间的能力称为空间维数,反过来,客体在单位空间内填充时间的能力可称为时间维数(至今尚无人研究此问题),因此,可在保留时空经典维数(拓扑维数)，即整数坐标的前提下,扩展运动客体的维数,如豪斯道夫维数、相似维数、关联维数、信息维数、容量维数、李雅普诺夫维数、复

5、维数、模糊维数等。目前所研究的这些整数维数和分数维数都是静态维数，即不考虑时间t的因素，如果在维数计算公式中引入了时间t因素，所求出的维数是随时间t而变化的，这样的维数称为动态维数，它描述了事物的发展变化规律。为此，本文以豪斯道夫分维数为基础，将时间因素引入其计算公式，求得一维动态豪斯道夫分维数Ｄft(1)式中x0、x1、x2...是时间序列Ｘt-1、Ｘt、Ｘt+1…是生成序列[2]，Ｄft为动态分维数，它描述了事物发展的否定之否定规律，其中，Ｋ为“肯定性”，Ｌ为“否定性”，利用Ｄft可以判断出系统的临界点，从而将序列分期，揭示其动态变化规律。2.2应

6、用实例2.2.1序列分期　　曹娥江（图１）是浙江省的八大水系之一，分析清楚曹娥江的河床演变对曹娥江的规划、整治具有重要的指导作用。在曹娥江河口，由于地形资料的缺乏而常通过分析低潮位的变化来分析其演变。曹娥江的潮位代表站为桑盆殿站，其潮位观测资料较详尽的记录起自1956年，至今已达40余年，潮位记录综合反映了系统各影响因素的综合作用和演变情况，其间由于人类活动、自然活动的影响使时间序列的一致性遭到破坏，需分期分段进行描述、统计、分析。经计算表明，常用的样本分割法难以较客观地进行序列分期，从而难以较客观地描述曹娥江的演变过程。为此本文采用动态分维数的方法，

7、对曹娥江低潮位序列进行突变分析，经计算发现分维数大于1的点共有7个（1964、1970、1977、1979、1985、1994、1996），其中与周围相邻点相比较差距较为明显的点有四个（1964、1970、1985、1996），即序列具有四个较显著的突变点（图２），对应的分维数分别为1.04，1.16，1.12，1.05。因此可根据这些临界点对桑盆殿低潮位分期统计如表１所示:由表可见在第一、三时期多年平均低潮位与多年平均低潮位较为接近，差距在0.1m左右，第二时期低潮位较低，低于多年平均水平0.3m，而第四时期有较大幅度增高，年均低潮位增高约0.5m，

8、最高低潮位则增高了1m。10图1曹娥江、钱塘江形势图MapoftheCaoeRiverandt