基于B_G分割算法的河川年输沙量突变分析.pdfVIP

水 利 学 报 2010 年12月 SHUILI XUEBAO 第41 卷 第12期 文章编号：0559-9350（2010 ）12-1387-06 基于B-G分割算法的河川年输沙量突变分析 1 1 2 2 李海彬 ，张小峰 ，胡春宏 ，王延贵 （1. 武汉大学 水资源与水电工程科学国家重点实验室，湖北 武汉 430072 ；2. 国际泥沙研究培训中心，北京 100048） 摘要：引入一种适合非线性和非平稳时间序列突变检测的新方法——Bernaola-Galavan （B-G）分割算法，结合二 项系数滤波，对嘉陵江流域年输沙量变化进行分析。在B-G 算法的基础上，对嘉陵江北碚站近50 年的实测年输 沙量资料进行突变年检测，并从流域的泥沙来源、年内分配和年际变化、以及气候变化、水利水保工程对产输沙 的影响等方面来探讨输沙量发生突变的物理背景。研究结果表明：1984年和1993年为嘉陵江北碚站年输沙量时 间序列的两个突变点，北碚站的年输沙量明显减少，气候变化对流域产输沙的影响是主要的、直接的；20 世纪 90 年代以来修建的一批大中型水库的减沙作用也十分显著；北碚站年输沙量发生突变是气候变化和人类活动共 同作用的结果。 关键词: B-G 分割算法；年输沙量；突变；嘉陵江流域 中图分类号: TV143 文献标识码：A 20 世纪60 年代中期以来，以Thom 的工作为先导而逐步建立起来的突变理论，不仅广泛应用于 ［1-2 ］ 气候、地震等领域 ，在水沙循环异常现象的分析和预测方面，也起到重要的作用。传统的研究河 ［3-4 ］ 川水沙量突变的方法有双累积曲线法、有序聚类分析法、M-K 秩检验法等 ，这些方法各有特点， 能不同程度上表现水沙突变的特征，但上述方法在处理非线性、非平稳时间序列资料时，存在一定 ［5-6 ］ ［7 ］ 的缺陷 。B-G 分割算法是2001 年首次由Bernaola-Galavan 等 人在研究心电图序列中提出来检测非 线性、非平稳时间序列突变的方法。与传统方法相比，该方法基于t检验将非平稳序列分割成多个具 有不同均值的平稳子序列，各子序列表征了不同的物理背景，分解得到的各均值段的尺度具有可变 ［8］ 性，不受方法本身的限制 。文献［9 ］通过构建理想时间序列，用B-G 算法对序列进行了突变检验， 证明了该算法能够有效地检测突变。本文在介绍B-G 分割算法的基础上，将其运用于河流年输沙量 时间序列的突变年检测，并对突变前后时段影响泥沙输移主要因子的变化进行分析，探讨输沙量发 生突变的物理背景，为探索河川输沙波动规律及其演变趋势打下基础。 1 研究方法 对于一个观测次数为N的非线性时间序列 ｛x（t）｝｛t=1，2 ，….. ，N｝，B-G 分割算法的计算步骤 为： （1）计算第i点左右两段的平均值，分别记为U和U（i=2 ，3 ，……，N-1）。 L R （2 ）为了衡量U和U之间的差异，计算统计量 L R T（i）=（U-U ）/S （1） L R D