Page 111 - 2024年第22卷第1期

分的变化分析较少 ［１６］ ，本文采用改进后的 ＩＨＡ法 ［１７］ 降低生态水文变量中的冗余性，保留原始数据中
的关键信息，分析雅鲁藏布江下游、尼洋河出口断面受上游水库影响的水文改变度，聚焦于对流量组
分分级以及缺资料地区的水文情势评估。以期为分析雅鲁藏布江地区水资源管理和生态治理提供科学
参考。

２ 研究概况

２．１ 研究区概况 雅鲁藏布江（２８°Ｎ—３１°Ｎ、８２°Ｅ—９７°Ｅ）发源于西藏自治区西南部喜马拉雅山中段
北麓的杰马央宗冰川，是世界上海拔最高、西藏最大的一条河流。雅鲁藏布江从源头 ～入海口全长约
２
３１６２ｋｍ，其中，中国境内长约 ２０５７ｋｍ。入海口处总流域面积 １７５万ｋｍ ，在中国境内流域面积约 ３３．
２
２０ 万ｋｍ ，雅鲁藏布江径流总体以降水补给为主，流域内多年平均降水量 ３００～５００ｍｍ，且年际变化
３
较小，年内分配不均。流域多年平均净流量约 １６５５亿ｍ ，年径流深与年降水量的分布趋势具有较好
的一致性 ［１８］ 。测站分布位置如图 １所示。

图 １ 研究区地理位置和观测站点示意图
Ｆｉｇ．１ Ｓｃｈｅｍａｔｉｃｄｉａｇｒａｍｏｆｇｅｏｇｒａｐｈｉｃａｌｌｏｃａｔｉｏｎａｎｄｏｂｓｅｒｖａｔｉｏｎｓｔａｔｉｏｎｓｉｎｔｈｅｓｔｕｄｙａｒｅａ

２．２ 研究数据 奴 下 水 文 站 位 于 雅 鲁 藏 布 江 流 域 干 流 ， 控 制 流 域 面 积 占 其 国 内 总 流 域 面 积 的
５７．６％。 尼 洋 河 出 口 断 面 的 更 张 测 站 为 雅 鲁 藏 布 江 的 上 游 控 制 站 ， 是 青 藏 高 原 重 要 的 水 文 站 点
之 一 ［１９］ 。 本 文 根 据 雅 鲁 藏 布 江 流 域 水 文 站 点 分 布 及 各 站 基 本 资 料 条 件 ， 选 用 奴 下 、 更 张 站 作
为 径 流 计 算 的 依 据 站 ， 进 行 分 析 评 价 采 用 的 是 １９６５—２０１８年 的 径 流 日 值 数 据 。 然 而 ， 由 于 资
料 获 取 难 度 问 题 ， 本 文 分 析 的 几 个 断 面 的 数 据 均 存 在 不 同 程 度 的 资 料 缺 失 ， 如 奴 下 站 点 缺
１９５８—１９５９年 、１９６９年 的 实 测 资 料 。 因 此 ， 本 文 采 用 分 布 式 水 文 模 型 ＥａｓｙＤＨＭ ［２０ － ２２］ 进 行 主
要 站 点 逐 日 流 量 数 据 的 插 补 还 原 。
２．３ 研究方法 论文基于 Ｒｉｃｈｔｅｒ等 ［９，２３ － ２４］ 提出的 ＩＨＡ指标进行分析，通过重新构建 ＩＨＡ指标体系将
降水时期划分为枯水期（ １—５月和 １１—１２月）和丰水期（６月—１０月）两个时段，将 ＩＨＡ指标划分为不
同时段内 ＩＨＡ指标，进一步优化 ＩＨＡ体系计算生态最相关水文指标 ＥＲＨＩｓ，从而达到减少水文指标的
同时尽可能多地保留原数据关键信息 ［１７］ 。然后基于 ＥＲＨＩｓ进一步优选参数，采用 ６种突变检验诊断
指标变化规律，最后应用水文改变度描述雅鲁藏布江径流变化情况。
２．３．１ 水文改变指标（ＩＨＡ） Ｒｉｃｈｔｅｒ提出的 ＩＨＡ体系，该方法定量分析了水文变异程度，目前该指标
体系被普遍用来代表水文情势变化特征，该指标体系包含 ３３个指标，目前被认为可以相对最全面、
最系统的来评估水文情势变化特征。该体系共分为 ５组，详细指标划分可参照文献 ［２３］ 。由于研究时段
内雅鲁藏布江未出现断流现象，故本研究不考虑 “断流天数” 这一变量 ［２５］ 。
２．３．２ 突变诊断 采用贝叶斯推断 ［２６］ 、曼－ 肯德尔（Ｍａｎｎ － Ｋｅｎｄａｌｌ）法 ［２７ － ２８］ 、勒帕热（Ｌｅ － Ｐａｇｅ）法 ［２８］ 、
ＢＧ（Ｂｅｒｎａｏｌａ － Ｇａｌｖａｎ）分割算法 ［２９］ 以及有序聚类法 ［３０］ 对改进后的指标突变点进行识别，用以排除单一
突变诊断方法的局限性以及数据漂移问题。

— １ ０ ９ —