近40年黄河中游径流情势变化分析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
理等的影响。由图 ３ 可以看出， 流域径流量表现为夏季大、冬季 小，月径流量多集中在 ６～９ 月。３ 月份有一个小汛期， 这可能与冰 雪融水补给有关。１９７８ 年后， 平均月径流量除了 ６ 月份比 １９７８ 年 前 稍 有 增 加 ， 其 他 月 份 普 遍 减 少 ， 而 且 湿 季 （ ７～１０ 月 ） 径 流 量 的减少比较显著。
势不显著。
２．３ 突变检验方法
检验序列突变的方法很多 ， 本文使用了滑动 ｔ 检验法［６］来推
估年径流序列的突变点。
ｔ 检 验［７］可 用 来 检 验 两 个 序 列 样 本 平 均 值 之 差 是 否 显 著 ， 从
而判断该序列中是否有突变发生。当两个子序列样本单元数目
ｎ１≠ｎ２ 时，统计量为［４］：
根据各个假设突变点的前和后长度为 ｎ 的两个子序列计算出
多个统计量 ｔ ， 取其中的最大绝对值进行统计量显著性检验。对
给定的信度水平 α， 若计算所得的 ｔ 值大于 ｔ 分布表中查出的临
界值， 则该点可被认为是序列的突变点。
２．４ 流量历时曲线
流量历时曲 线（ＦＤＣ） 定 义 为 在 某 个 时 段 内 某 个 流 量 与 大 于
一个流域径流情势的变化是气候条件与下垫面共同作用的 结果， 充分了解径流情势、气候因素和人类作用间的关系对制定 可 持 续 的 流 域 管 理 规 划 起 到 十 分 重 要 的 作 用 ［２］。本 文 以 黄 河 中 游 无定河水系水土保持重点治理区— ——岔巴沟为例， 分析流域径 流的年内变化和年际变化， 推估水土保持措施对水文序列的显 著影响干扰点， 进而根据突变点前后两阶段水文序列特性的变 化， 来分析人类活动对该黄河中游径流情势的影响。
２ 研究对象和方法
２．１ 研究对象和资料 无定河是黄河中游河口镇—龙门区间右岸的最大支流，平
均每年约有 ２．２ 亿 ｔ 泥沙随 着 滚 滚 洪 流 输 入 黄 河 ， 占 黄 河 年 输 沙量 １６ 亿 ｔ 的 １４％，是多年来黄河中游水土 保 持 工 作 的 重 点 地 区。本文所研究的岔巴沟位于无定河流域西部的黄土丘陵沟壑 区， 见图 １， 集水面积 １８７ｋｍ２， 每逢暴雨干、支流均会出现较大的 洪水和高含沙水流。
月平均径流量（万 ｍ３） Ｑ５０（万 ｍ３） Ｃｖ Ｑ５（万 ｍ３） Ｑ９５（万 ｍ３）
９０．５ ４８．８６ １．３４ ３１０．６９ １３．９３
验是否存在趋势的统计量为：
Ｕ＝
!
１／２
（ １）
［Ｖａｒ（ !） ］
式中，
!＝ ４Ｐ － １
（ ２）
ｎ（ ｎ－ １）
对于随机序列有 Ｅ（ !） ＝０，则：
Ｖａｒ（ !） ＝ ２（ ２ｎ＋５）
（ ３）
９ｎ（ ｎ－ １）
当 ｎ 增加时， Ｕ 很快收敛于标准正态分布。当原假设为该序列无
趋势时， 一般采用双侧检验， 在给定信度水平 α下， 若 Ｕ ＞Ｕα／２， 则拒绝原假设， 即认为趋势是存在的， 否则接受原假设， 认为趋
表 ２ 年径流序列特性表
１９５９￣２０００
Ｑ（ 万 ｍ３）
Ｃｖ
８６３．４
０．４８
１９５９￣１９７７
Ｑ（ 万 ｍ３）
Ｃｖ
１０３２．５
０．４８
１９７８￣２０００
Ｑ（ 万 ｍ３）
Ｃｖ
７０９．６５
０．３６
３．２ 径流的年内变化 径 流 的 年 内 变 化 主 要 是 受 到 降 雨 季 节 性 变 化 、温 度 、流 域 管
面积随时间的变化。从图 ２ 可以看出， ７０ 年代坝地面积 显 著 增
加， １９７７ 年到达最大， １９７８ 年左右大量淤堤坝被水毁， 面积显著
减少。
２．２ 趋势检验
目前，在时间序列数理统计分析方法中，坎德尔秩次相关检
验 法 ［５］是 分 析 时 间 序 列 平 稳 性 、检 验 趋 势 性 的 有 效 手 段 。
图 １ 岔巴沟流域图 黄河中游大规模水土保持措施是人类活动导致黄河中游 径流情势变化的一个主要因素。本文选取了无定河水土保持重 点治理区岔巴沟为研究区域， 根据岔巴沟流域 １９５９～２０００ 年内 年径流、年降雨序列及各 项 水 保 措 施 的 占 地 面 积 随 时 间 的 变 化 情况来进行研究。岔巴沟流域水土保持治理工作始于 １９５６ 年， 主 要 措 施 主 要 为 梯 田 、造 林 、种 草 和 淤 地 坝 。淤 地 坝 是 工 程 措 施 中的主要部分， 对径流的影响相对比较迅速、明显。它的发展过
１ 引言
黄河是中 国 第 二 大 河 ， 流 经 ９ 省 ， 长 ５ ４６４ｋｍ， 流 域 面 积 ７５．２ 万 ｋｍ２， 从 西 到 东 横 跨 青 藏 高 原 、内 蒙 古 高 原 、黄 土 高 原 和 黄淮海平原四个地貌单元。流域地势西高东低， 西部河源地区由 一系列高山组成， 常年积雪， 冰川地貌发育； 东部主要由黄河冲 积平原组成； 中部地区流经世界上最大的黄土高原， 由于黄土土 质疏松、垂 直 节 理 发 育 ， 水 土 流 失 严 重 ， 导 致 了 严 重 的 生 态 恶 化 和经济损失。２０ 世纪 ６０ 年代以来黄河中游内进行了诸 如 修 建 水 库 、引 水 灌 溉 等 大 量 的 水 事 活 动 和 土 地 开 发 、水 土 保 持 等 一 系 列大规模生产活动， 不同程度地改变了流域下垫面条件， 使得它 的径流情势在过去的几十年里经历了较大的变化， 包括年径流 量 的 减 少 ，及 更 严 重 的 断 流 和 季 节 性 断 流 ［１－ ３］。
第２７卷 第 ５期 ２００７年１０月
水文 ＪＯＵＲＮＡＬ ＯＦ Ｃ水ＨＩＮ文Ａ ＨＹＤＲＯＬＯＧＹ
Ｖｏｌ．２７ Ｎｏ．５ Ｏｃｔ．第， ２２０７０卷７
近 ４０ 年黄河中游径流情势变化分析
王 蕊 １， 夏 军 ２
（ １． 武汉大学水资源与水电工程科学国家重点实验室， 湖北武汉 ４３００７２； ２．中国科学院地理科学与资源研究所， 北京 １００１０１）
对序列 ｘ１， ｘ２， …， ｘｎ， 先确定所有对偶值（ ｘｉ， ｘｊ， ｊ＞ｉ） 中 ｘｉ＜ｘｊ 的 出现个数（ 设为 Ｐ） 。如果按顺序前进的值全部大于前一个值， 则
是 一 种 上 升 趋 势， Ｐ＝（ ｎ－ １） ＋（ ｎ－ ２） ＋…１， 系 为 等 差 级 数 ， 其 总 和
为 Ｐ＝（ ｎ－ １） ｎ／２。如果序列全部倒过来， 则 Ｐ＝０， 即为下降趋势。检
显著干扰， 年径流量减少了 ３０％以上， 且径流年际、年内变化均有所减小。此外， 流域径流情势随淤地坝有起有落的发展过
程出现了一系列相应的变化。
关键词： 径流情势； 黄河中游； 流量历时曲线； 变化分析
中图分类号： Ｐ３３３．１
文献标识码： Ａ
文章编号： １０００－ ０８５２（２００７）０５－ ００７４－ ０４
摘 要： ２０ 世纪 ６０ 年代以来， 黄河中游地区一系列大规模人类生产活动如引水灌溉、水土保持等，使得 黄 河 中 下 游 的 径 流
情势在过去的几十年里经历了较大的变化， 造成流域水土流失、生态恶化和经济损失严重。本文以黄河中游的岔巴沟流域
为例，研究了 １９５９～２０００ 年黄河中游地区人类活动对径流情势的影响。结果表明： 流域径流在 １９７８ 年受到水土保持活动的
１９７８￣２０００
１９５９￣２０００
统计量值
－ ２．４６
－ ０．７８
－ ０．２
－ １．５７
坎德尔法
显著性
＊＊
ｎｓ
ｎｓ
ｎｓ
注： " 取＝０．１， ＊＊ 表示趋势非常显著， ｎｓ 表示趋势不显著。
年降雨量 １９５９￣１９７７
０．４１ ｎｓ
１９７８￣２０００
０．２３ ｎｓ
７６
水文
第２７卷
到了 ０．３６。可见， １９７８ 年以后， 年径流量和径流年际变化都比以 前有所减小。
表 ３ 岔巴沟流域 １９７８ 年前后月径流序列特性表
指标
Ｑ５０ （万 ｍ３） Ｃｖ
Ｑ５（万 ｍ３） Ｑ９５（万 ｍ３）
岔巴沟
１９５６ － １９７７ ４９．２５
１９７８ － ２０００ ３８．８８
１．４３
１．３７பைடு நூலகம்
３１０．６９
１７４．５７２
１３．９３
１６．１１３
但对于 １０％～９０％的中、小来水， 各年代径流量与 ６０ 年代比变化 很小。
收稿日期： ２００６－ １１－ １５ 基金项目： 国家自然科学基金项目（４０６７１０３５） 作者简介： 王蕊（１９８２－ ）， 女， 吉林通化人， 武汉大学硕士研究生， 主要从事水文及水资源研究。
第５期
王蕊等： 近 ４０ 年黄河中游径流情势变化分析
７５
程有起有落［４］， 图 ２ 显示了岔巴沟流域 １９５６～２０００ 年淤地坝占地
ｔ＝ ｘ" １－ ｘ"２
（ ４）
#Ｓ· １ ＋ １ ｎ１ ｎ２
且
２
２
#Ｓ＝ ｎ１ｓ１ ＋ｎ２ｓ２ ｎ１＋ｎ２－ ２
（ ５）
式中： ｘ" １、ｘ" ２ 和 ｓ１ 、ｓ２ 分别代表两个子序列样 本 的 平 均 值 和 均 方 差的无偏估计值。
由于突变点和它的位置是未知的， 滑动 ｔ 检验不是将整个
序列分为两个子序列进行 ｔ 检验， 而是在原序列上滑动突变点，
图 ３ 岔巴沟流域 １９７８ 年前后两阶段月径流过程图
３．４ 流量历时曲线的变化 图 ４ 是 岔 巴 沟 流 域 １９７８ 年 前 后 两 阶 段 的 月 流 量 历 时 曲
线。从图中可以看出， 同样频率下径流大多在减少， 频率在 １５％～９５％之 间 的 来 水 减 少 比 较 稳 定 。在 极 大 或 极 小 来 水 时 ， 径 流的变化幅度较大。表 ３ 对两个阶段月径流序列的一些特性 进 行了比较。１９７８ 年以后 ， 变 差 系 数 Ｃｖ 比 以 前 减 少 了 ４．１９％， Ｑ５０ 减少 ２１．０６％， Ｑ５ 减少 ４３．８１％， Ｑ９５ 增加了 １５．６７％。
等于该流量所对应的时间之间的相关关系， 流量历时曲线的形
状由降雨类型、流域大小和流域的 地 形 特 征 来 决 定 。ＦＤＣ 在 研
究径流情势中应用很广， 它综合地描述了某流域的径流从枯水
到洪水整个阶段的全部特征， 可以较好地反映流域的降雨径流
特性。它也可以研究某频率下径流的变化， 进而来评估由于人
１９８０ｓ
１９９０ｓ
图 ４ １９７８ 年前后流量历时曲线对比图
为了进一步探明径流变化规律， 画出了 ６０￣９０ 年代的 ４ 条 流量历时曲线， 并用曲线表示出各频率下 ７０～９０ 年代相对于 ６０ 年 代 的 所 减 少 的 径 流 量 ， 见 图 ５、图 ６。 图 形 表 明 ： 相 对 于 ６０ 年 代 ， ９０ 年 代 的 曲 线 变 化 最 大 ， ７０ 年 代 和 ８０ 年 代 的 曲 线 较 接 近 ，
由 表 ２ 中 计 算 结 果 可 知 ： １９７８～２０００ 年 的 年 平 均 径 流 量 比 １９５９～１９７８ 年减少了 ３０％以上， 年径 流 变 差 系 数 也 由 ０．４８ 减 少
序列名称 时段
表 １ １９５９￣２０００ 年年径流趋势检验结果表
年径流量
１９５９￣２０００
１９５９￣１９７７
类 活 动 和 气 候 变 化 导 致 的 径 流 情 势 的 改 变 ［８］。
３ 计算结果
３．１ 径流的年际变化 径流量的变化是流域径流情势变化的一个重要方面。采用
坎德尔秩次相关法对年径流序列和年降雨量进行趋势检验， 统 计 量 计 算 结 果 见 表 １。 当 取 显 著 性 水 平 α＝０．０５ 时 ， 临 界 值 为 １．６４５， 则 年 径 流 序 列 有 显 著 的 下 降 趋 势 ， 年 降 雨 量 序 列 无 显 著 趋势。然后用滑动 ｔ 检验法检验年径流序列的突变点， 计算可 得 ： 当 取 信 度 水 平 α＝０．１ 时 ， 年 径 流 序 列 在 １９７８ 发 生 突 变 。 若 １９７８ 是年径流序列的突 变 点，则 其 前 后 的 水 文 序 列 应 无 明 显 趋 势 性 。为 进 一 步 验 证 突 变 点 的 可 靠 性 ，对 突 变 点 前 后 两 个 时 段 的 水文序列进行趋势性检验， 表 １ 中的检验结果进一步说明以 １９７８ 年为年径流序列的突变点是合理的。
表 ４ 列出了各年代的一些径流特性， 数据表明： 与其他年 代相比， ９０ 年代月径流序列的变差系数 Ｃｖ 最大。月平均径流量 先逐渐减少到 ９０ 年代有所增加， 且 ９０ 年 代 的 Ｑ５ 和 Ｑ９５ 恢 复 到 与 ６０ 年代很接近。
表 ４ ６０￣９０ 年代年径流量特性表
年代
１９６０ｓ
１９７０ｓ
势不显著。
２．３ 突变检验方法
检验序列突变的方法很多 ， 本文使用了滑动 ｔ 检验法［６］来推
估年径流序列的突变点。
ｔ 检 验［７］可 用 来 检 验 两 个 序 列 样 本 平 均 值 之 差 是 否 显 著 ， 从
而判断该序列中是否有突变发生。当两个子序列样本单元数目
ｎ１≠ｎ２ 时，统计量为［４］：
根据各个假设突变点的前和后长度为 ｎ 的两个子序列计算出
多个统计量 ｔ ， 取其中的最大绝对值进行统计量显著性检验。对
给定的信度水平 α， 若计算所得的 ｔ 值大于 ｔ 分布表中查出的临
界值， 则该点可被认为是序列的突变点。
２．４ 流量历时曲线
流量历时曲 线（ＦＤＣ） 定 义 为 在 某 个 时 段 内 某 个 流 量 与 大 于
一个流域径流情势的变化是气候条件与下垫面共同作用的 结果， 充分了解径流情势、气候因素和人类作用间的关系对制定 可 持 续 的 流 域 管 理 规 划 起 到 十 分 重 要 的 作 用 ［２］。本 文 以 黄 河 中 游 无定河水系水土保持重点治理区— ——岔巴沟为例， 分析流域径 流的年内变化和年际变化， 推估水土保持措施对水文序列的显 著影响干扰点， 进而根据突变点前后两阶段水文序列特性的变 化， 来分析人类活动对该黄河中游径流情势的影响。
２ 研究对象和方法
２．１ 研究对象和资料 无定河是黄河中游河口镇—龙门区间右岸的最大支流，平
均每年约有 ２．２ 亿 ｔ 泥沙随 着 滚 滚 洪 流 输 入 黄 河 ， 占 黄 河 年 输 沙量 １６ 亿 ｔ 的 １４％，是多年来黄河中游水土 保 持 工 作 的 重 点 地 区。本文所研究的岔巴沟位于无定河流域西部的黄土丘陵沟壑 区， 见图 １， 集水面积 １８７ｋｍ２， 每逢暴雨干、支流均会出现较大的 洪水和高含沙水流。
月平均径流量（万 ｍ３） Ｑ５０（万 ｍ３） Ｃｖ Ｑ５（万 ｍ３） Ｑ９５（万 ｍ３）
９０．５ ４８．８６ １．３４ ３１０．６９ １３．９３
验是否存在趋势的统计量为：
Ｕ＝
!
１／２
（ １）
［Ｖａｒ（ !） ］
式中，
!＝ ４Ｐ － １
（ ２）
ｎ（ ｎ－ １）
对于随机序列有 Ｅ（ !） ＝０，则：
Ｖａｒ（ !） ＝ ２（ ２ｎ＋５）
（ ３）
９ｎ（ ｎ－ １）
当 ｎ 增加时， Ｕ 很快收敛于标准正态分布。当原假设为该序列无
趋势时， 一般采用双侧检验， 在给定信度水平 α下， 若 Ｕ ＞Ｕα／２， 则拒绝原假设， 即认为趋势是存在的， 否则接受原假设， 认为趋
表 ２ 年径流序列特性表
１９５９￣２０００
Ｑ（ 万 ｍ３）
Ｃｖ
８６３．４
０．４８
１９５９￣１９７７
Ｑ（ 万 ｍ３）
Ｃｖ
１０３２．５
０．４８
１９７８￣２０００
Ｑ（ 万 ｍ３）
Ｃｖ
７０９．６５
０．３６
３．２ 径流的年内变化 径 流 的 年 内 变 化 主 要 是 受 到 降 雨 季 节 性 变 化 、温 度 、流 域 管
面积随时间的变化。从图 ２ 可以看出， ７０ 年代坝地面积 显 著 增
加， １９７７ 年到达最大， １９７８ 年左右大量淤堤坝被水毁， 面积显著
减少。
２．２ 趋势检验
目前，在时间序列数理统计分析方法中，坎德尔秩次相关检
验 法 ［５］是 分 析 时 间 序 列 平 稳 性 、检 验 趋 势 性 的 有 效 手 段 。
图 １ 岔巴沟流域图 黄河中游大规模水土保持措施是人类活动导致黄河中游 径流情势变化的一个主要因素。本文选取了无定河水土保持重 点治理区岔巴沟为研究区域， 根据岔巴沟流域 １９５９～２０００ 年内 年径流、年降雨序列及各 项 水 保 措 施 的 占 地 面 积 随 时 间 的 变 化 情况来进行研究。岔巴沟流域水土保持治理工作始于 １９５６ 年， 主 要 措 施 主 要 为 梯 田 、造 林 、种 草 和 淤 地 坝 。淤 地 坝 是 工 程 措 施 中的主要部分， 对径流的影响相对比较迅速、明显。它的发展过
１ 引言
黄河是中 国 第 二 大 河 ， 流 经 ９ 省 ， 长 ５ ４６４ｋｍ， 流 域 面 积 ７５．２ 万 ｋｍ２， 从 西 到 东 横 跨 青 藏 高 原 、内 蒙 古 高 原 、黄 土 高 原 和 黄淮海平原四个地貌单元。流域地势西高东低， 西部河源地区由 一系列高山组成， 常年积雪， 冰川地貌发育； 东部主要由黄河冲 积平原组成； 中部地区流经世界上最大的黄土高原， 由于黄土土 质疏松、垂 直 节 理 发 育 ， 水 土 流 失 严 重 ， 导 致 了 严 重 的 生 态 恶 化 和经济损失。２０ 世纪 ６０ 年代以来黄河中游内进行了诸 如 修 建 水 库 、引 水 灌 溉 等 大 量 的 水 事 活 动 和 土 地 开 发 、水 土 保 持 等 一 系 列大规模生产活动， 不同程度地改变了流域下垫面条件， 使得它 的径流情势在过去的几十年里经历了较大的变化， 包括年径流 量 的 减 少 ，及 更 严 重 的 断 流 和 季 节 性 断 流 ［１－ ３］。
第２７卷 第 ５期 ２００７年１０月
水文 ＪＯＵＲＮＡＬ ＯＦ Ｃ水ＨＩＮ文Ａ ＨＹＤＲＯＬＯＧＹ
Ｖｏｌ．２７ Ｎｏ．５ Ｏｃｔ．第， ２２０７０卷７
近 ４０ 年黄河中游径流情势变化分析
王 蕊 １， 夏 军 ２
（ １． 武汉大学水资源与水电工程科学国家重点实验室， 湖北武汉 ４３００７２； ２．中国科学院地理科学与资源研究所， 北京 １００１０１）
对序列 ｘ１， ｘ２， …， ｘｎ， 先确定所有对偶值（ ｘｉ， ｘｊ， ｊ＞ｉ） 中 ｘｉ＜ｘｊ 的 出现个数（ 设为 Ｐ） 。如果按顺序前进的值全部大于前一个值， 则
是 一 种 上 升 趋 势， Ｐ＝（ ｎ－ １） ＋（ ｎ－ ２） ＋…１， 系 为 等 差 级 数 ， 其 总 和
为 Ｐ＝（ ｎ－ １） ｎ／２。如果序列全部倒过来， 则 Ｐ＝０， 即为下降趋势。检
显著干扰， 年径流量减少了 ３０％以上， 且径流年际、年内变化均有所减小。此外， 流域径流情势随淤地坝有起有落的发展过
程出现了一系列相应的变化。
关键词： 径流情势； 黄河中游； 流量历时曲线； 变化分析
中图分类号： Ｐ３３３．１
文献标识码： Ａ
文章编号： １０００－ ０８５２（２００７）０５－ ００７４－ ０４
摘 要： ２０ 世纪 ６０ 年代以来， 黄河中游地区一系列大规模人类生产活动如引水灌溉、水土保持等，使得 黄 河 中 下 游 的 径 流
情势在过去的几十年里经历了较大的变化， 造成流域水土流失、生态恶化和经济损失严重。本文以黄河中游的岔巴沟流域
为例，研究了 １９５９～２０００ 年黄河中游地区人类活动对径流情势的影响。结果表明： 流域径流在 １９７８ 年受到水土保持活动的
１９７８￣２０００
１９５９￣２０００
统计量值
－ ２．４６
－ ０．７８
－ ０．２
－ １．５７
坎德尔法
显著性
＊＊
ｎｓ
ｎｓ
ｎｓ
注： " 取＝０．１， ＊＊ 表示趋势非常显著， ｎｓ 表示趋势不显著。
年降雨量 １９５９￣１９７７
０．４１ ｎｓ
１９７８￣２０００
０．２３ ｎｓ
７６
水文
第２７卷
到了 ０．３６。可见， １９７８ 年以后， 年径流量和径流年际变化都比以 前有所减小。
表 ３ 岔巴沟流域 １９７８ 年前后月径流序列特性表
指标
Ｑ５０ （万 ｍ３） Ｃｖ
Ｑ５（万 ｍ３） Ｑ９５（万 ｍ３）
岔巴沟
１９５６ － １９７７ ４９．２５
１９７８ － ２０００ ３８．８８
１．４３
１．３７பைடு நூலகம்
３１０．６９
１７４．５７２
１３．９３
１６．１１３
但对于 １０％～９０％的中、小来水， 各年代径流量与 ６０ 年代比变化 很小。
收稿日期： ２００６－ １１－ １５ 基金项目： 国家自然科学基金项目（４０６７１０３５） 作者简介： 王蕊（１９８２－ ）， 女， 吉林通化人， 武汉大学硕士研究生， 主要从事水文及水资源研究。
第５期
王蕊等： 近 ４０ 年黄河中游径流情势变化分析
７５
程有起有落［４］， 图 ２ 显示了岔巴沟流域 １９５６～２０００ 年淤地坝占地
ｔ＝ ｘ" １－ ｘ"２
（ ４）
#Ｓ· １ ＋ １ ｎ１ ｎ２
且
２
２
#Ｓ＝ ｎ１ｓ１ ＋ｎ２ｓ２ ｎ１＋ｎ２－ ２
（ ５）
式中： ｘ" １、ｘ" ２ 和 ｓ１ 、ｓ２ 分别代表两个子序列样 本 的 平 均 值 和 均 方 差的无偏估计值。
由于突变点和它的位置是未知的， 滑动 ｔ 检验不是将整个
序列分为两个子序列进行 ｔ 检验， 而是在原序列上滑动突变点，
图 ３ 岔巴沟流域 １９７８ 年前后两阶段月径流过程图
３．４ 流量历时曲线的变化 图 ４ 是 岔 巴 沟 流 域 １９７８ 年 前 后 两 阶 段 的 月 流 量 历 时 曲
线。从图中可以看出， 同样频率下径流大多在减少， 频率在 １５％～９５％之 间 的 来 水 减 少 比 较 稳 定 。在 极 大 或 极 小 来 水 时 ， 径 流的变化幅度较大。表 ３ 对两个阶段月径流序列的一些特性 进 行了比较。１９７８ 年以后 ， 变 差 系 数 Ｃｖ 比 以 前 减 少 了 ４．１９％， Ｑ５０ 减少 ２１．０６％， Ｑ５ 减少 ４３．８１％， Ｑ９５ 增加了 １５．６７％。
等于该流量所对应的时间之间的相关关系， 流量历时曲线的形
状由降雨类型、流域大小和流域的 地 形 特 征 来 决 定 。ＦＤＣ 在 研
究径流情势中应用很广， 它综合地描述了某流域的径流从枯水
到洪水整个阶段的全部特征， 可以较好地反映流域的降雨径流
特性。它也可以研究某频率下径流的变化， 进而来评估由于人
１９８０ｓ
１９９０ｓ
图 ４ １９７８ 年前后流量历时曲线对比图
为了进一步探明径流变化规律， 画出了 ６０￣９０ 年代的 ４ 条 流量历时曲线， 并用曲线表示出各频率下 ７０～９０ 年代相对于 ６０ 年 代 的 所 减 少 的 径 流 量 ， 见 图 ５、图 ６。 图 形 表 明 ： 相 对 于 ６０ 年 代 ， ９０ 年 代 的 曲 线 变 化 最 大 ， ７０ 年 代 和 ８０ 年 代 的 曲 线 较 接 近 ，
由 表 ２ 中 计 算 结 果 可 知 ： １９７８～２０００ 年 的 年 平 均 径 流 量 比 １９５９～１９７８ 年减少了 ３０％以上， 年径 流 变 差 系 数 也 由 ０．４８ 减 少
序列名称 时段
表 １ １９５９￣２０００ 年年径流趋势检验结果表
年径流量
１９５９￣２０００
１９５９￣１９７７
类 活 动 和 气 候 变 化 导 致 的 径 流 情 势 的 改 变 ［８］。
３ 计算结果
３．１ 径流的年际变化 径流量的变化是流域径流情势变化的一个重要方面。采用
坎德尔秩次相关法对年径流序列和年降雨量进行趋势检验， 统 计 量 计 算 结 果 见 表 １。 当 取 显 著 性 水 平 α＝０．０５ 时 ， 临 界 值 为 １．６４５， 则 年 径 流 序 列 有 显 著 的 下 降 趋 势 ， 年 降 雨 量 序 列 无 显 著 趋势。然后用滑动 ｔ 检验法检验年径流序列的突变点， 计算可 得 ： 当 取 信 度 水 平 α＝０．１ 时 ， 年 径 流 序 列 在 １９７８ 发 生 突 变 。 若 １９７８ 是年径流序列的突 变 点，则 其 前 后 的 水 文 序 列 应 无 明 显 趋 势 性 。为 进 一 步 验 证 突 变 点 的 可 靠 性 ，对 突 变 点 前 后 两 个 时 段 的 水文序列进行趋势性检验， 表 １ 中的检验结果进一步说明以 １９７８ 年为年径流序列的突变点是合理的。
表 ４ 列出了各年代的一些径流特性， 数据表明： 与其他年 代相比， ９０ 年代月径流序列的变差系数 Ｃｖ 最大。月平均径流量 先逐渐减少到 ９０ 年代有所增加， 且 ９０ 年 代 的 Ｑ５ 和 Ｑ９５ 恢 复 到 与 ６０ 年代很接近。
表 ４ ６０￣９０ 年代年径流量特性表
年代
１９６０ｓ
１９７０ｓ