基于GAMLSS 模型的宜昌站年径流序列趋势分析.pdfVIP

第67 卷第11期 地 理 学 报 Vol.67, No.11 2012 年11 月 ACTA GEOGRAPHICA SINICA Nov., 2012 基于GAMLSS 模型的宜昌站年径流序列趋势分析 江 聪, 熊立华 (武汉大学水资源与水电工程科学国家重点实验室, 武汉430072) 摘要：在变化环境下，研究非一致性水文序列的变化趋势具有重要意义。常规的趋势分析方 法一般只能分析水文序列均值的线性趋势。本文引入GAMLSS 模型对宜昌站1882-2009 年间的 年平均流量序列和年最小月流量序列分别进行趋势分析，将序列的趋势分析从均值扩展至均 方差 (或变差系数) 、偏态系数等其他统计参数。研究发现宜昌站年平均流量序列的均值有明 显线性减少的趋势，而宜昌站年最小月流量序列线性趋势不明显。在此基础上，建立基于多 项式回归的GAMLSS 模型，结果表明宜昌站年最小月流量序列并非平稳序列，其均值表现为 非线性的趋势变化，偏态系数呈现线性的趋势变化。 关键词：趋势分析；年径流量；GAMLSS 模型；时变统计参数；长江 1 研究背景 在变化环境下，水文序列由于受到气候变化和人类活动影响两方面的扰动，不再满足 一致性的要求，导致无法直接运用基于一致性假设的方法进行水文频率分析计算。因此， 在变化环境下，进行水文序列的非一致性分析对保证水文频率分析计算结果的合理性有重 要意义。水文序列的趋势分析是非一致性分析中的重要内容，常用的水文序列趋势分析方 法主要有：滑动平均法、Mann-Kendall 趋势检验法、Spearman 秩次相关法、线性趋势回归 [1] [2-4] 检验法等 。这些方法已经在实际中得到了广泛的应用 ，比如，国内水文学者已经运用 [5] 上述方法对于宜昌水文站年径流序列进行了广泛的研究。熊立华等 应用Mann-Kendall 、 Spearman 秩次相关检验法和线性趋势回归检验法对长江宜昌站年径流的趋势进行了检 验，结果发现年平均和年最小日径流序列有明显减少的趋势，而年最大日流量序列的趋势 [6] 不明显。许继军等 运用Mann-Kendall 法对宜昌站1950-2000 年间的径流序列进行了趋势 [7] 分析，发现年平均径流有减少的趋势。秦年秀 选取宜昌站，分别应用1882-2000 年、 1870-2000 年和1950-2000 年的月平均流量资料，对年代际、月径流、季节性径流的变化以 及径流的变化趋势及突变进行了分析研究，并使用M-K 法来检验径流的趋势变化。 但是，上述提到的常用趋势分析方法一般只考虑了水文序列均值的线性趋势变化，这 显然与实际情况是不相符的。事实上，水文序列的趋势会存在于水文序列的任何统计参 [8] 数 ，不仅仅是均值，还包括方差、偏态系数等。而且，水文序列的趋势不仅会出现在序 列的全过程，也会出现在序列的某一段时期，或在不同时期内出现不同的趋势。另外，水 文序列的趋势可能是线性的，也可能是非线性的。因此本文中将引入位置、尺度、形状的 广义可加模型 (Generalized Additive Models for Location, Scale and Shape ，简称 GAMLSS)[9] ，基于该模型对宜昌站年平均流量序列和年最小月流量序列的所有统计参数进 行线性以及非线性趋势变化进行分析，以弥补常规方法在水文序列趋势分析中的不足。 收稿日期：2012-07-24; 修订日期：2012-08-15 基金项目：国家自然科学基金重大项目; 国家自然科学基金项目 [Foundation: Major Program of the National Natural Science