周期图方法在径流序列分析中的应用.pdfVIP

周期图方法在径流序列分析中的应用 李 万 绪 (陕西省石 泉水 电站 ) 摘要 对水文序列 中多年周期的变化问题 ，常用功率谱方法、最大熵谱 方 法 厦 周 期 圈方法采进行分析。本文应用周期 图分析原理，建立 了一种适用于小 样 本 容 量 (30~50年 )水文序列的多车周期识别方法。该方法还改进 了周期 圈方法 中对周期 的识剐所产生 的不连 续性。 主题 词 周期 径泣 圈 径流计算 水文气象要素的周期现象已为许多研究所肯定。但是，除了年、季节性变化 以外，一般 周翔变化都是不确定的，其周期的长度也在一定的范围内波动。虽然一般的周期变化不具严 格的确定性，然而分析水文气象序列的频率特性，仍对生产实践具有指导意义，特别在长期 水文气象预报领域。 水文气象序列属时闯上的离敲序列，而离散痔列周期结构分析的尉l奶图方法，需使被分 析序列具有正态分布特性，以便利用某些统计分布函数来检验分析结论。对水文气象要素， 一 般认为多雨地区的年降水量，多雨季节长时段 (如一月以上 )降水量的概率分布属正态分 布 “】。而流域某出口断面上长时问段平均径流是该 区域 内长时段降雨的直接反映，故也属 正态分布。 对于具有正态分布的径流序列，一般也认为其具有平稳性，即序列的期望值不随时问变 化，并且可 以用时间平均代替总体平均。对于某些平稳性不十分确定的序列或人类活动影响 较大的流域，其径流序列也可用模拟方法检验其是否具有平稳性，如辩程检验方法 “’等。具 有正态分布的平稳时间序列，可用周期图方法判断其周期结构。 1 分析原理 1．1 周期豳 对具有纯随机噪声的离散序列孔 (t—l，…，N)，设其数学模型是 ． xt= ∑ cicos(2n，i+ )+8t， 1，… ． k I t 上式中 、 (c}、 {，t}为常数，t j}是在 (一 ． )内均匀分布的 独 立 随机 变 鹫。 {e)是独立于 { }的纯随机过程，且 E (￡．}=0， E{￡t }=a (1)式也可 以写成 25 I= ∑ (0icos2~!it+扫jsin2~!jt)+8I (2) |。 I 其 中dj—cICOS J，bJ掌sin j。 在给定频率，时，d和bi的最小二乘估计式 ㈨ 是 一 n N 。 os2~!i ， b tsln2，i (4) ， 那末，N个观察序列 t，…， 的周期匿定义为 了c，垒I e一，：ml。，一≤，≤。 上式 中，为虚数单位， “垒”表示定义于。 实际计算中Ix(，)的频率取离散值0，]／N，2／N，…，故记 1P垒 IN(，P) ，=P／N，P=0，1，…， [N／23。 以，代 替 (3)， (4)两式中的，t得周期图估计式 JP=N／2(d +b ) ． 州应频率，，下的周期 『口， Tp=N／P 振幅C， √ + 。 1．2 周期捡验 对获得 的 (N／2] 个 J值，计算Fisher统计疑g， s垒 ㈩ j，