基于ERHIs法的雅鲁藏布江流域水文情势变化研究.pdf

第 卷第 期 年 月中国水利水电科学研究院学报（中英文），收稿日期：；网络首发时间：网络首发地址：?基金项目：华能集团总部科技项目（）；流域水循环模拟与调控国家重点实验室重点项目（）；长江勘测规划设计研究有限责任公司自主创新项目（）作者简介：杨靖（），博士生，高级工程师，主要从事水利水电工程管理研究。：文章编号：（）基于 法的雅鲁藏布江流域水文情势变化研究杨靖，蔺子琪，赵庆绪，谢静，廖卫红（华能西藏水电安全工程技术研究中心，四川 成都 ；四川大学 水利水电学院，四川 成都 ；华中科技大学 土木与水利工程学院，湖北 武汉 ；长江勘测规划设计研究有限责任公司，湖北 武汉 ；水利部水网工程与调度重点实验室，湖北 武汉 ；中国水利水电科学研究院 流域水循环模拟与调控国家重点实验室，北京 ）摘要：为探究变化环境下雅鲁藏布江流域环境流量组分与水文情势变化情况，采用改进后生态最相关水文指标 （）、（）法定量分析了 年间不同阶段不同环境组分下的雅鲁藏布江流域内水文情势的变化情况，并结合 种突变检验方法对指标进行验证。结果表明，通过 插补数据是可行的，参数 优化 后的 系 数最 大为 ，最小为 ，均在合理范围内，在缺资料地区运用较好。针对雅鲁藏布江 个断面强变异指标改变度计算结果均表明雅江中下游地区的水文改变度为高度改变。在指标选取方面，生态水文改变过程的指标体系（，）参 数 之 间 的 相 关 性 较 高，指 标 间 的 冗 余 问 题 十 分 突 出。指标具有一定的合理性，不仅有效缓解了 的冗余性问题，还有利于抓住最关键的生态水文变量，对生态水文研究、流域水资源管理和生态保护都具有重要的参考价值和借鉴意义，为短缺水文资料地区的生态评估提供参考。关键词：水文改变度；水文情势；雅鲁藏布江；生态最相关水文指标；模型中图分类号：文献标识码：?研究背景雅鲁藏布江是世界上海拔最高、西藏最大的一条河流，也是我国西南地区非常重要的跨界河流之一。雅鲁藏布江流域径流主要由降水、冰雪融水和地下水补给组成 。由于流域面积广阔，流域内的降水量、蒸发量以及地形、冰雪覆盖面积、植被等条件差异很大，流域内的径流补给类型复杂多样 。精确描述径流特征及其变化的基础是获得大量高精度流量数据，特别是在水文站点缺乏的地区。分布式水文模型为流域径流数据的延展和加密提供了重要的支撑 ，基于此，可量化评估人类活动、下垫面和气候变化等变化条件对水文特性影响程度，从而可为区域水资源合理开发利用提供重要的依据 。针对气候变化和水资源开发利用而导致的水文改变度的分析一直为生态水文研究领域的热点问题 ，等 在 年提出了一套能够评估水生态水文改变过程的指标体系（，），目前被广泛应用于生态水文领域 。但 等 学者分析结果表明，众多水文变量之间存在十分突出的冗余性问题，而 指标与大多数水文变量之间存在较高的信息重叠现象。相对于早期研究中只考虑单一水文变量，尽管 体系较完善，但是水文变量仍然较多，变量之间也存在着自相关性和严重的冗余性问题。法目前被认为可以相对最全面、最系统的来评估水文情势变化特征。此外，雅鲁藏布江地区的水文情势研究较多，但具体对于环境流量组分的变化分析较少 ，本文采用改进后的 法 降低生态水文变量中的冗余性，保留原始数据中的关键信息，分析雅鲁藏布江下游、尼洋河出口断面受上游水库影响的水文改变度，聚焦于对流量组分分级以及缺资料地区的水文情势评估。以期为分析雅鲁藏布江地区水资源管理和生态治理提供科学参考。研究概况 研究区概况雅鲁藏布江（、）发源于西藏自治区西南部喜马拉雅山中段北麓的杰马央宗冰川，是世界上海拔最高、西藏最大的一条河流。雅鲁藏布江从源头 入海口全长约 ，其中，中国境内长约 。入海口处总流域面积 万，在中国境内流域面积约 万，雅鲁藏布江径流总体以降水补给为主，流域内多年平均降水量 ，且年际变化较小，年内分配不均。流域多年平均净流量约 亿，年径流深与年降水量的分布趋势具有较好的一致性 。测站分布位置如图 所示。图 研究区地理位置和观测站点示意图 研究数据奴下 水 文 站 位 于 雅 鲁 藏 布 江 流 域 干 流，控 制 流 域 面 积 占 其 国 内 总 流 域 面 积 的 。尼洋河出口断面的更张测站为雅鲁藏布江的上游控制站，是青藏高原重要的水文站点之一 。本文根据雅鲁藏布江流域水文 站 点 分 布 及 各 站 基本 资 料 条 件，选 用 奴 下、更张 站作为径流计算的依据站，进行分析评价采用的是 年的径流日值数据。然 而，由 于资料获取难度 问 题，本 文 分 析 的 几 个 断 面 的 数 据 均 存 在 不 同 程 度 的 资 料 缺 失，如 奴 下 站 点 缺 年、年的实测资料。因 此，本 文 采 用 分 布 式 水 文 模 型 进 行 主要站点逐日流量数据的插补还原。研究方法论文基于 等 ，提出的 指标进行分析，通过重新构建 指标体系将降水时期划分为枯水期（月和 月）和丰水期（月 月）两个时段，将 指标划分为不同时段内 指标，进一步优化 体系计算生态最相关水文指标 ，从而达到减少水文指标的同时尽可能多地保留原数据关键信息 。然后基于 进一步优选参数，采用 种突变检验诊断指标变化规律，最后应用水文改变度描述雅鲁藏布江径流变化情况。水文改变指标（）提出的 体系，该方法定量分析了水文变异程度，目前该指标体系被普遍用来代表水文情势变化特征，该指标体系包含 个指标，目前被认为可以相对最全面、最系统的来评估水文情势变化特征。该体系共分为 组，详细指标划分可参照文献 。由于研究时段内雅鲁藏布江未出现断流现象，故本研究不考虑“断流天数”这一变量 。突变诊断采用贝叶斯推断 、曼 肯德尔（）法 、勒帕热（）法 、（）分割算法 以及有序聚类法 对改进后的指标突变点进行识别，用以排除单一突变诊断方法的局限性以及数据漂移问题。（）贝叶斯推断贝叶斯推断 是利用贝叶斯定理来修正特定假设的可能性，该方法可以充分利用先验信息和样本信息，得到统计结论。计算公式为：（，）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）法 法 作为常见的非参数统计检验方法，无需样本遵从一定的分布，适用于类型变量和顺序变量，计算简便，可以明确突变开始的时间并指出突变区域。计算公式如下：，（），当 ，当 ，（）（）（槡），（）其中，是第 时刻数值个数的累计数，（），（）分别是 的均值和方差，统计量 为标准正态分布。按时间序列 逆序，再重复上述过程，同时使 ，其中 ，。（）勒帕热法勒帕热（）法 是一种具有两个样本不受分布影响的非参数判别法，计算公式如下：（）（）（）（）（）（）（）为统计量，式（）是勒帕热统计量。注意到当样本量足够大时，渐近具有自由度为 的分布。其余参数含义同上。（）分割算法 分割算法 主要针对非平稳时间序列进行突变检测，其核心思想是通过显著性检验找到突变点的左右两部分序列的差异最大的位置。计算公式如下：珚珔（）（）左 右?左右()?（）珚、珔分别为突变点左右两侧的均值，这里 和 分别为左、右两边子序列的标准差，左和 右分别为左、右子序列的样本量个数，为参数。其余参数含义同上。（）有序聚类法有序聚类分析法 进行突变点的识别，该方法核心思想是在分类时不打乱次序的情况下，通过计算突变点前后点位的离差平方和寻求最优分割点。计算公式为：（）（珚）（珔）则当 （）时的 为最优分割点，推断为突变点。其余参数含义同上。水文改变度基于 指标体系，等 在 年提出了用于判断 指标改变程度的变化范围法（，），并定义了 阈值概念，即在 的基础上，将各指标发生几率的 和 作为 的上下阈值，用来作为评价各指标受影响程度的标准 。为了量化 受干扰后的变化程度，采用水文改变度来评估 ，其定义如下：；()式中：为第 个 指标的水文变化程度；为第 个 指标变化在之后（之前）仍在 阈值范围内的实际观测年数；为第 个 突变（前）后预计落在 阈值内的年数，可以表示为 ，其中 为突变前（后）第 个 落在 阈值内的比例（本文取 ），是受突变（前）和后影响的时间序列中记录的总年数。为流域整体的水文改变度，规定水文变化程度在 之间为无变化或低程度变化，为中度变化，之间是高度变化。模型验证与指标确定 测站验证依据流域和数据资料情况，论文对奴下、更张断面的 到 年水文日模拟精度进行了参数率定，将优化得到的参数对流域各主要断面 年的逐日水文过程进行了模拟，各站 系数评价见表 ，说明论文采用的 模型用于日径流资料的延展和加密是可行的。依据流域和数据资料情况，本研究对奴下、更张断面挑选洪水场次进行了参数率定和验证，并对精度进行了评价，针对流域特点，主要针对洪峰误差精度对方案合格程度进行了评定，其中水库入库数据主要通过出库和蓄量反推，本身就存在误差，因此对其洪峰许可误差做了适当放大。选取 年间的 组径流过程进行率定，径流过程拟合效果较好，洪峰流量误差和峰现时间误差基本符合要求。对 年间的 组径流过程样本进行评定，洪水场次合格率为 。模型采用 全局参数优化算法进行参数率定 。表 模型模拟径流 系数表 控制站名称所在河流紧邻下游控制站 系数控制站名称所在河流紧邻下游控制站 系数控制站名称所在河流紧邻下游控制站 系数江孜年楚河日喀则 羊八井堆龙曲羊村 唐加拉萨河拉萨 旁多拉萨河唐加 拉萨拉萨河羊八井 奴各沙雅鲁藏布江羊村 拉孜雅鲁藏布江无 泽当（二）雅砻河羊村 更张尼洋河奴下 日喀则年楚河奴各沙 羊村雅鲁藏布江奴下 奴下雅鲁藏布江德兴 巴河桥巴河更张 工布江达尼洋河巴河桥 更张（二）尼洋河奴下 的冗余性分析 体系包含 个指标，由于研究时段内雅鲁藏布江未出现断流现象，故本研究不考虑“断流天数”这一变量 。论文对所选择的 个 指标进行冗余性分析，主要针对指标间的相关系数，认为相关系数小于 为冗余度较小。从图 可以看出，部分水文变量之前展现了较强的相关性，月流量，月流量，月流量，年最小 日、日、日、日流量之间的相关性较强，高流量年内发生次数，高流量年内平均历时，上升率，下降率之间相关性较强。从图 可以看出，年最小 日、日、日流量，年最大 日、日、日、日、日流量之间也具有较强的相关性，其相关系数甚至超过了 。结果表明，雅鲁藏布江的站点断面的 指标变量相互之间显然存有高度自关性，即指标冗余性问题非常明显。图 奴下断面 指标相关关系热力图 图 更张断面 指标相关关系热力图 指标体系重建鉴于 变量之间存在的冗余性问题比较突出，本文通过重新构建 指标体系将降水时期划分为枯水期（月和 月）和丰水期（月 月）两个时段，将 指标划分为不同时段内 指标，进一步优化 体系计算 ，在缓解 变量冗余性时，减少水文指标，在保证全面描述流域水文情势变化时尽可能多地保留原数据关键信息，通过涵盖丰枯时段以及指标所在组别信息来确定。流域的径流量以及径流特征受到季节变化的影响，由于不同季节降水量的不同，可将降水时期分为丰水期和枯水期。由于一年内枯水期历史长短与流域自然地理及气象条件相关，根据丰水期和枯水期的特点将雅鲁藏布江奴下断面 月至 月划分为丰水期，将 月、月划分为枯水期（表 ）。表 重建后 指标 不同时期组别重分配水文改变指标枯水期第一组：月流量 月流量和 月流量第二组：年极端流量及历时最小 、日流量、基流指数第三组：年极端流量出现时间年最小流量出现时间第四组：低流量发生频率和历时低流量年内发生次数、低流量年内平均历时第五组：流量改变率和频率上升率、下降率、流量逆转次数丰水期第一组：月流量 月流量第二组：年极端流量及历时最大 、日流量第三组：年极端流量出现时间年最大流量出现时间第四组：高流量发生频率和历时高流量年内发生次数、高流量年内平均历时第五组：流量改变率和频率上升率、下降率、流量逆转次数为使筛选出的 指标具有较高的可靠性，本文通过选取样本年并挑选各样本年所对应的相似年，计算相似年各 指标的贡献率，选取贡献率排在前六位且重复出现次数最多的 指标为 指标。相似年通常指在表征同一流域有两年或多个年份具有相同的水文指标或径流特征或其他流域特征时，都可以称这两年或多个年份之间在某一特征前提下存在相似性，即这两年或多个年互为相似年 。相似年在各个领域的应用十分广泛，在天文、气象、水利等各行业的预报都有应用。本文从奴下、更张断面 年数据中分别随机选取十年作为相应断面的样本年，将能够反应径流特征的水文统计值作为选取样本年所对应相似年的基准，运用二维欧氏距离筛选方法结合统计值对选取好的样本年进行相似年份的选取 。为对应不同时期 指标体系，将 年流量数据按枯水期（月和 月）和丰水期（月 月）分开，分别计算 年数据的枯水期统计值和丰水期统计值。统计值包括：均值、方差、极差、中位数、分位数、分位数、变异系数、偏度和峰度。图 指标确定步骤 指标确定通过计算各样本年所对应十个相似年枯水期和丰水期 指标的变异系数，将十个相似年中各年枯水期和丰水期变异系数前六个指标挑选出来并记录各指标在十年中出现的重叠次数，挑选出枯水期和丰水期出现次数排在前六的指标，将这十二个指标再次按照出现次数排序，查看前六个指标在一级指标中的覆盖度是否合理。通过计算各样本年所对应十个相似年枯水期和丰水期 指标的变异系数，将十个相似年中各年枯水期和丰水期变异系数前六个指标挑选出来并记录各指标在十年中出现的重叠次数，挑选出枯水期和丰水期出现次数排在前六的指标，将这十二个指标再次按照出现次数排序，查看前六个指标在一级指标中的覆盖度是否合理。经合理性确认，将前六个指标确定为 。在合理性判定的基础上计算各指标的相关性，经计算后发现只有两个指标存在相关性，其他指标之间的相关系数都小于 为不相关，将这些指标确定为筛选指标。具体步骤见图 。指标的变异系数反应该指标在所有 指标中的鉴别能力。变异系数越大，表明该指标的信息分辨能力就越强，从而保证了筛选出的指标对流域影响最大 。初步筛选出的 指标见表 、表 。表 更张断面 指标 指标变异系数累计次数年最大流量出现时间低脉冲流量平均历时高脉冲流量平均历时年最小流量出现时间上升率 月表 奴下断面 指标 指标变异系数累计次数月最小 日 最小 日年最小流量出现时间高脉冲流量平均历时上升率 指标的合理性分析为验证 个断面选出的 指标相互独立且具有较高的代表性，现对选出的 指标进行合理性分析，本文中主要通过 指标的覆盖性分析和冗余性分析来进行合理性分析。在减少 变量冗余性的同时，以保证 指标可以均匀覆盖不同时期以及不同组别，有效保留原始信息。（）覆盖性分析覆盖性分析通过对各 指标在组别中的覆盖范围以及各指标在不同时期的覆盖范围进行指标的合理性分析，研究断面指标覆盖范围分析结果见表 、表 。表 更张断面覆盖性分析 指标出现时期第三组年最大流量出现时间丰水期第四组低脉冲流量平均历时枯水期第四组高脉冲流量平均历时枯水期第三组年最小流量出现时间枯水期第五组上升率枯水期第一组 月丰水期表 奴下断面覆盖性分析 指标出现时期第一组月枯水期第二组最小 日枯水期第二组最小 日丰水期第三组年最小流量出现时间丰水期第四组高脉冲流量平均历时丰水期第五组上升率枯水期由表可以得出，更张断面六个 指标涉及到了第一组、第三组至第五组且达到丰枯水期全覆盖，奴下断面的六个 指标涉及到了第一组至第五组指标且达到丰枯水期全覆盖。可以看出雅鲁藏布江下游的 个断面挑选出的 指标在组别中的覆盖率以及各指标在不同时期的覆盖率都很高，得出挑选出的 指标具有一定的合理性。（）冗余性分析为验证 个断面选出的 指标相互独立且具有较高的代表性，现对选出的 指标进行相关性分析，确定 指标各指标间的相关性，图 展示了各指标间的相关系数。结果表明：大部分指标之间存在较小的相关性，即信息并不存在冗余性问题，则可最终确定为 指标。图 更张断面和奴下断面 相关关系热力图 计算结果分析 初步诊断初步诊断是对已经挑选出的 指标进行指标变异强度的确定，将强变异指标确定为最终用于识别变异年份的指标。变异系统的初步诊断采用 系数法 分别对雅鲁藏布江的更张水文站、奴下水文站的 指标进行变异检验。经检验得出尼洋河站的 月指标 系数为 属于强变异，奴下站的 月指标 系数为 属于强变异。计算结果见表 、表 。表 更张站 指标变异强度 指标年最大流量出现时间年最小流量出现时间低脉冲流量平均历时高脉冲流量平均历时上升率 月 系数 表 奴下站 指标变异强度 指标最小 日月年最小流量出现时间高脉冲流量平均历时上升率最小 日 系数 详细诊断根据初步诊断结果针对各断面挑选出的强变异指标进行详细诊断。本文利用水文变异诊断系统的跳跃诊断法和气候突变检测方法包括 方法、有序聚类法、贝叶斯推断方法、方法和 （）分割法等方法进行突变年份的确定。更张断面有 种方法诊断 月指标在 年发生显著突变，有 种方法诊断其在 年发生显著突变。奴下断面有 种方法诊断 月指标在 年发生显著突变。具体诊断结果见表 。表 雅鲁藏布江各断面强变异指标的突变诊断结果 诊断方法更张断面奴下断面 方法 贝叶斯推断方法 、有序聚类法 勒帕热（）法 、分割算法 、结果表明，不同突变方法对突变点的识别有所差异，由于 方法需要人为确定子序列长度，可能会由于突变点的漂移造成误差，只用一种突变方法难以全面的评判指标突变的结果。综上所述，运用多种突变检验方法分析更张、奴下站 年径流数据得出以下结论：对于更张断面：月指标突变时段为 年。对于奴下断面：月指标突变时段综合序列为 年。造成奴下断面指标综合突变时间较长的原因可能为上游支流拉萨河直孔及旁多水库对下游的调蓄作用，年冲久水库投入运行后，改变了下游河段的天然径流过程。因而，奴下断面通过 种突变方法中体现出这种突变状态，其中 系数指标显示年最小流量出现时间、年最小 日流量指标均有中等程度的改变。强变异指标变化程度本文通过判断变异前后时间序列的长短来确定稳定时段，通过以上对奴下和更张断面水文情势突变点的研究，认为时间序列较长一段的水文情势更为稳定 。将更张断面 年作为稳定时段，将奴下断面 年作为稳定时段。以 等 定义的阈值概表 雅鲁藏布江 个断面强变异指标改变度计算结果 断面强变异指标改变度?改变程度更张 月 高度奴下月 高度念即将稳定时段强变异指标发生几率的 作为 的上阈值，发生几率的 作为 的下阈值作为标准。论文分别计算各断面强变异指标包括更张断面 月指标、奴下断面 月指标的水文改变度（表 ）。结果显示，两个断面的强变异指标都发生了高度改变。综上得出，更张断面和奴下断面整体水文改变度为高度改变。结论（）通过 插补雅江流域站点数据是可行的，其中，参数优化后的 系数最大为 ，最小为 ，均在合理范围内，在缺资料地区运用较好。此外，方法有效的通过减少变量简化水文响应推求过程，以获得水文情势变化情况，提高水资源调配效率，对于水资源管理和生态系统保护都具有重要的参考价值和借鉴意义。（）不仅有效缓解了 的冗余性问题，还有利于抓住最关键的生态水文变量。从 指标相关性热图可以看出，参数之间的相关性较高，说明指标间的冗余问题十分突出。指标突变结果分析表明，运用欧式距离结合相似年筛选 指标与多种突变结果的横向对比以及 自身的纵向对比，都表明所选取的生态最相关指标是合理的。（）气候变化和人类活动对环境组分流量和水文状况的构成都有影响，而大规模水利建设是其中的主导因素。针对更张断面和奴下断面的水文改变度分析可以得出，在受到多布水电站建设的影响，尼洋河下游水文改变度差异显著，这和张冬艳、刘海平等人的研究结果相一致 。（）上游水电站建设运行对水文情势有不同程度的影响，水文情势变化总体属高度改变，非汛期所受影响程度比汛期更大。干旱季节的径流变化较大，除受到上游水电站影响外，气温的增温效应导致冰川融雪的补给量增加也在影响雅鲁藏布江的年内径流分配，但这种变化还需相关的资料和分析去验证。参考文献：李浩，牛乾坤，王宣宣，等 年雅鲁藏布江流域径流演变规律分析 水土保持学报，（）：（，（）：（）张仪辉，刘昌明，梁康，等雅鲁藏布江流域降水时空变化特征 地理学报，（）：（，（）：（）班春广，徐宗学，左德鹏，等高寒流域水文模拟与径流水源解析 以雅鲁藏布江帕隆藏布上游流域为例 北京师范大学学报（自然科学版），（）：（，：（），（）：（）关铜垒，刘佳嘉，周祖昊，等考虑山区 平原区的分布式水文模型计算单元划分方法研究 水文，（）：（，（）：（），（）：葛金金，张汶海，渠晓东，等湄公河下游生态水文情势变化研究 人民长江，（）：（，（）：（）康泽璇，王芳，刘扬，等基于 法的大通河上中游水文节律变化 南水北调与水利科技（中英文），（）：（，（）：（）李光录，樊立娟基于 法的长江源区生态水文情势变化 长江科学院院报，（）：，（，?，（）：，（），（）：，：，（）：，（）：郭文献，焦旭洋，周昊彤，等嘉陵江水文情势变化及其鱼类影响研究 长江流域资源与环境，（）：（，（）：（）邹杰平，陈明，刘友存，等赣江上游典型流域水文情势分析 安徽农业大学学报，（）：（，（）：（）王鸿翔，查胡飞，卓志宇，等基于 法四水流域水文情势变化评估 中国水利水电科学研究院学报，（）：（，（）：（），（）：常俊超基于?法的黄河中游环境流量组分与水文情势变化研究 水电能源科学，（）：（?，（）：（）程俊翔，徐力刚，姜加虎生态最相关水文指标的优选及其在洞庭湖环境流量估算中的应用 湖泊科学，（）：（，（）：（）刘天仇雅鲁藏布江水文特征 地理学报，（）：（，（）：（）王浩，王芳，严登华，等中华水塔面临的挑战与应对战略 中国水利水电科学研究院学报（中英文），（）：（，（）：（）雷晓辉，廖卫红，蒋云钟，等分布式水文模型 模型 水利信息化，（）：（，（）：（）雷晓辉，廖卫红，蒋云钟，等分布式水文模型 （）：理论方法 水利学报，（）：（，：，（）：（）雷晓辉，蒋云钟，王浩，等分布式水文模型 （）：应用实例 水利学报，（）：，（，（）：，（），（）：，？，（）：，：，（）：?李道扬，何幼桦混合非参数回归的贝叶斯推断 上海大学学报（自然科学版），（）：（，（）：（），（?）：董弟文，阿布都热合曼哈力克，王大伟，等近 年和田河源流区径流特征及对气候变化的响应 中国水利水电科学研究院学报，（）：（，（）：（）李海彬，张小峰，胡春宏，等基于 分割算法的河川年输沙量突变分析 水利学报，（）：（，（）：（）李保敏，张俊芝，吴旭，等基于 法和有序聚类法的径流变化特征研究 海河水利，（）：，（，（）：，（）张婷，魏凤英华南地区汛期极端降水的概率分布特征 气象学报，（）：（，（）：（）刘克琳，王银堂，胡四一，等 最优分割法在汛期分期中的应用 水利水电科技进展，（）：，（，（）：，（）曾金凤，刘祖文，刘友存，等基于 法的东江源区生态水文情势变化分析 水土保持通报，（）：（，（）：（）李舒，齐青松基于 法的人类活动对窟野河流域水文情势影响评估 水电能源科学，（）：，（，（）：，（）涂玉律，李英海，郭家力，等基于 法的隔河岩水库下游水文情势分析 人民珠江，（）：（，（）：（）刘贵花，朱婧蠧，熊梦雅，等基于变动范围法（）的信江水文改变及生态流量研究 水文，（）：（，（）：（），：，（）：，?（），（）：，（）：高祥涛数据挖掘在水文相似年查找中的应用研究 计算机工程与应用，（）：（，（）：（）杜瑞莉，闫伟杰，冯慧敏综合优势相似法在短期气候预测中的应用 河南气象，（）：（，（）：（）石宝峰，迟国泰基于信息含量最大的绿色产业评价指标筛选模型及应用 系统工程理论与实践，（）：（，（）：（）谢平，陈广才，雷红富基于 系数的水文变异分析方法 应用基础与工程科学学报，（）：（，（）：（）张东艳，吴运卿，李妮基于 检验的尼洋河流域水文变量演变趋势分析 中国农村水利水电，（）：（，（）：（）刘海平，郭传波，叶少文，等人类干扰下西藏尼洋河生态系统的健康评估?年中国水产学会学术年会，南昌，（，：?，（），（，；，；，；，；，；，）：，（）（），（），：；（责任编辑：祁伟）