对长江宜昌站月平均径流量突变的分析(.docVIP

对长江宜昌站月平均径流量突变的分析( 江剑民1 谷湘潜2 游性恬1 (1:中国气象局培训中心，北京100081;2:中国气象科学研究院，北京100081) 提 要 本文简要叙述了一阶矩突变扫描式t检验和二阶矩突变扫描式 F检验算法，并应用于长江宜昌水文站1882-2002年间历年月平均径流量的实测序列，划分出24个年际至十年际时间尺度的相对5级旱、涝时段及其变率特征，与有关长江上游洪水、干旱灾害的历史 记载相吻合. 关键词 突变 月径流 宜昌水文站 长江 分类号 P332.3 科学家们通常使用滑动平均或低通滤波算法来揭示时间序列的长期变化特征.然而, Karl Riebsame[1] 应用经典的学生氏t检验分析了美国气温和降水量10-20年际的振荡. Yamamoto, et al.[2] 采用修改的t检验来检测“气候跃变”. Goossen Berger[3] 采用Mann-Kendall秩和检验来鉴别气候突变. 这些方法每次计算都只能检测一个突变点. 在Foufoula-Georgiou Kumar[4] 编辑的《地球物理学中的小波分析》一书中，介绍了用小波变换来揭示大气湍流观测长时间序列中多次突变现象，以及分析两序列间突变相干性的方法.但是，小波变换不能给出判别突变点的统计显著性检验标准. 将小波变换技术移植到学生氏t检验和F检验[5]上, 分别拓展成一阶矩突变扫描式t检验[6，7] 和二阶矩突变扫描式 F检验算法.这样既可使学生氏t检验和F检验具有检测长时间序列中不同时间尺度上多次突变的现象，又能给出判别突变点的统计显著性检验标准.用扫描式t检验来检测长江宜昌水文站1882-2002年间历年月平均径流量的实测序列中不同时段相邻子序列之间平均值（一阶矩）的突变，从而划分出年际至十年际时间尺度的相对旱、涝期.又用扫描式F检验来检测该序列中不同时段相邻子序列之间均方差（二阶矩）的突变，从而划分出年际至十年际时间尺度上，长江宜昌水文站月平均径流量变化的相对稳定和不稳定阶段.然后综合上述扫描式t检验和扫描式 F检验的结果，概括得出24个相对稳定或不稳定的旱、涝期，与有关长江上游洪水、干旱灾害的历史记载相吻合. 1 方法 1.1 扫描式t 检验 定义一阶矩突变扫描式t检验的统计量t(n, j)，对于长序列中第j时间点前后两个相邻子序列样本量（n）相等时，为下列算式[6，7]： (1) 其中 ； 其中n表示两相邻子序列的长度,n=2, 3, ...， N/2, 也可以选择适当的间隔.j是检测突变现象的参考点, j=n+1, n+2, ..., N-n . t检验要求序列应是独立的和遵从正态分布的.在实际应用中需加注意, 或者进行适当的处理. 这里，采用一种称作“Table-Look-Up测试”[8] 来订正检验子序列的非独立性.该表是通过蒙托卡罗方法，按各子序列非独立程度下，一阶自相关系数与样本量大小计算得到的t检验统计显著性临界值.本文选用置信水平为95%的临界值作为显著性检验标准. 为使图表直观, 分析方便, 再取比值： (2) 作为一阶矩突变扫描式t检验的显著性指数. 显然, 当|tr(n,j)| 1.0时, 可判断该突变是通过了置信度为0.05的统计显著性标准的.其中tr(n,j) -1.0表示变化趋势为下降的突变事件, tr(n,j) 1.0代表上升的突变事件. 1.2 扫描式F检验 类似地,定义二阶矩突变扫描式F检验的统计量Fr(n,j)，对于长序列中第j时间点前后两个相邻子序列样本量（n）相等时，为下列算式： (Sj12 /Sj22 )/Fα, 当 Sj2Sj1 时 Fr(n,j)= 0, 当Sj2 = Sj1 , 或Sj1 = 0, 或Sj2 = 0时 (3) (Sj22/ Sj12 )/Fα, 当Sj2 Sj1 时 其中 Sj1 和 Sj2 的算法与上述 (1) 式相同，n =2, 3, ... N/2, j =n+1, n+2,..., N-n. Fα是经过对子序列非独立性进行订正之后的有效自由度下，置信水平在95%上F检验的显著性临界值.这样，Fr (n,j) -1.0表示在置信水平95%上子序列均方差（二阶矩）显著减小