水文序列突变点识别方法的研究.pdfVIP

水文序列突变点识别方法研究 荆、贵山 (黄河水利委员会上游水文水资源局 兰州 730030) 摘要：当前在水文序列变异点识别中常采用的几种统计方法都对数据有较多假设，当假设不满 足时，识别结果通常并不理想。通过探求多方法联合检验的方式，应用滑动T检验法、Men— 确的结果，并采用该方法对黄河兰州站70多年径流序列进行了变异(或突变)点识别。 关键词：R／S法；Mann—Kendall法；滑动T检验；变异点；水文时间序列；兰州站 水文序列作为水文观测及分析应用的基本形式，运用中必须满足准确性、代表性和一 致性。水文循环、自然条件的剧烈变化和流域内人类活动的影响通常会引起流域水文系统 运行方式的改变，进而使径流序列发生变异、呈现出阶段性特征(使直接观测得到的资料 出现前后不一致性)。了解和诊断径流序列的变异规律，对水文分析、模拟和预测中统计 方法的选择以及制定防洪减灾对策和水环境治理措施都具有重要意义。当前在水文序列变 异(或突变)点识别中常采用的几种统计方法都对数据有较多假设，当假设不满足时，识 别结果通常并不理想[1]。在用多方法联合判别序列变异点方面，文献2中给出了一种较为 完整的多方法联合检验识别变异点的方法，但认为工作量较大，同时存在方法重复(滑动 T检验法与有序聚类法检验原理相似，结果相同)的可能；本文考虑检验效果及方便应 结合(滑动平均)过程线及成因验证，可以获得较为满意的结果，在很大程度上避免了单 一方法判别上的不足。 1基本原理及方法 1．1 Men—Kendall法 在时间序列分析时，Men—Kendall法是世界气象组织推荐的非参数检验方法，并已 广泛地用来分析降水、径流和气温等要素地时间序列的变化情况。改检验不需要样本服从 一定的分布，也不受少数异常值的干扰，适合水文、气象等非正态分布的数据，计算简 便‘2|。具体计算方法如下‘2。4]。 对于具有咒个样本量的时间序列z，构成一秩序列： k S女一∑mk一1，2，…，咒 (1) f=1 其中 作者简介：孙贵山(1965一)，男，甘肃会宁县人，工程师，学士，主要从事水文测验及水文水资源分析研究。 水文序列突变点识别方法研究口 旷J1五i置，产1’2，…，i 0 n—1 z。≤zi。刮_’…“ 可见，秩序列S。是i时刻值大于第歹时刻值个数的累计数。 在时间序列随机独立的假定下，定义统计量： UFk--鐾型，k一1，2，…，n (2) √ynr(Sk) 互独立且有相同连续分布时，可由下式算出： fVa椰。)=坐；乌型丑 M j (3) I E(S。)一掣 UF。为标准正态分布，它是按时间序列z的顺序X。，z。，…，z。计算出的统计量序 列，给定显著性水平a，查正态分布表，若f晒。I阢，则表明序列存在明显的趋势 变化。 按时间序列z的逆序z。，z。一。，…，z，，再重复上述过程，同时使UBe(序列的逆 UB；，如果两条曲线出现交点，且交点在临界直线直接，那么交点对应的时刻就是突变开 始的时刻。 1．2滑动T检验法 传统的T检验法：服从正态分布的两序列总体的分布函数为Ft(z)、F2(z)，今从总 体F，(z)和F。(z)中分别抽取容量为咒。和他的两个样本，要求检验原假设：F，(z)一 F2(z)。对于滑动点r前后，n—n。+竹。定义统计量为[a-4]：