基于多源遥感数据的黄河数字流域模型改进.pdf

年 月水利学报 第 卷第 期文章编号：（）收稿日期：；网络首发日期：网络首发地址：?基金项目：国家自然科学联合基金项目（，）；清华大学水沙科学与水利水电工程国家重点实验室课题（，）作者简介：薛源（），博士，助理研究员，主要从事遥感水文学的研究。：通信作者：覃超（），博士，助理研究员，主要从事土壤侵蚀、河床演变的研究。：基于多源遥感数据的黄河数字流域模型改进薛源，覃超，吴保生，张?，李丹，傅旭东，（清华大学 水圈科学与水利工程全国重点实验室，北京 ；清华大学 水利部水圈科学重点实验室，北京 ；清华大学 水利水电工程系，北京 ；中国煤炭科工集团有限公司 煤炭科学研究总院有限公司 应急科学研究院，北京 ）摘要：山区河流水系众多，往往因地势险要难以现场量测获得河流基础信息，是典型的缺资料地区。本文结合多源遥感提取的 个河流断面及 个水文站实测断面，建立了黄河中游各级别河流的断面概化模型。融合传统线状 水系及河流表面、概化断面构建河道边界形态。针对黄河数字流域模型，基于不同级别河流的河道边界条件建立了断面判别模块，改进其汇流模型。以黄河一级支流皇甫川流域为例，以 年的资料率定模型，以 年汛期日径流及年内 场典型洪水过程验证模拟效果。结果表明，改进后的模型输入参数均可由遥感观测获取，且日流量模拟的纳什效率系数高于 ，流速及水深的模拟误差分别为 、，全时期洪峰流量的模拟误差平均为 、峰现时间误差约 ，均较改进前明显提升。研究成果可为山区河流水文过程模拟及河流地貌演化等研究提供重要的基础数据和技术参考。关键词：山区河流；多源遥感；河流断面；黄河数字流域模型中图分类号：文献标识码：?研究背景我国地形复杂，山区约占国土面积的 ，发育着众多的山区河流。但山区河流往往因地势险要，无法通过现场量测获得流量等基础资料，是典型的水文缺资料地区，制约了水资源调查、河流水沙通量等重要研究的开展。为缓解这一现状，学者们多通过卫星遥感观测、经验估算或反演山区河流的径流过程 ，但所得数据往往针对单一断面或河段，丰度有限、精度较差，方法也难以推广 。因此，利用水文模型能够模拟流域径流过程的优势，结合多源遥感数据大范围提取的河流关键参数，模拟山区等缺资料地区河流的径流过程是当今的研究热点 。分布式水文模型，能较充分考虑水沙运动的物理机制及下垫面的空间差异，反映了接近真实的水文过程 。其中黄河数字流域模型（，），以多源遥感及地理信息系统等技术模拟全流域的降雨 径流过程 。模型依托 数据，以流域分级理论为依据，将全流域分为 个尺度：坡面 沟道单元、单一子流域、多个子流域的组合和全流域，通过基于坡面 沟道单元的产流、汇流模型，在子流域、流区和全流域水系模拟水沙过程 。近年来，被诸多学者改进并应用，如李铁键、刘家宏等 利用 模拟了黄河流域的径流过程，验证了 模拟不同下垫面条件、多种类型河流水沙过程可靠性；张昂等 在考虑降雨、河宽等参数空间分布的基础上改进了 ，模拟了黄河源区 年的径流过程，进一步提高了缺资料地区径流模拟的精度，文献 也通过考虑植被、淤地坝的影响改进 。由于推广性较好，还被用于模拟和预测珠江流域、闽江流域、南方山区河流、长江上游等区域的径流过程 。尽管 能模拟不同尺度、类型流域较真实的径流过程，但仍有较大提升空间。如 将河道断面简化为“”形，虽然能较好反映低级别河流的断面特征 ，但不能用以模拟高级别河流 ，即使通过模型率定、调参等方式降低因断面形态失真带来的误差 ，仍会使模拟结果不能反映真实的下垫面条件、水流情况，因此对土壤特性、水力要素的模拟偏离实际。可见，虽然 的模拟效果整体较好，但仍需基于河流几何信息，获取更真实的断面形貌沿程变化，改进模型、提升模拟效果。本文以位于黄土高原的黄河一级支流皇甫川流域为研究区域，以多源遥感数据为主要数据源，提取了河流表面信息及各级河流的断面形态，与 河网提取工具得到的线状河流水系组合，得到了流域河道的边界条件。基于各级河流的断面形态，建立了 断面判别模块，改进了其汇流模型，提升了径流过程的模拟效果，提高了洪水过程与水力要素的模拟精度。本研究可为探究基于多源遥感数据获得较完整的河道边界条件、提升水文模拟效果及河流地貌演化等研究提供参考。基于断面判别模块的模型改进 是以 河网为核心，基于物理机理、分布式参数和连续计算时段，能较准确描述黄河流域下垫面环境，较精确模拟黄河径流过程的分布式水文模型 。模拟流域径流模拟时，涉及其产流、汇流模型，本文主要针对汇流模型进行改进（图 ）。图 产流及改进的汇流模型框架 产流、汇流模型及数据输入坡面产流是流域径流过程的起点，为了适应黄河中游超渗产流的特点，产流模型采用土水势概念模拟地表入渗过程的变化，获得较精确的产流过程 。其后，以每一个坡面及与其相邻沟道组成 的最小汇流单元 坡面 沟道单元，并假定坡面 沟道单元产流直接作用于单元出口，按自上游至下游的拓扑关系流入下一级河网中完成水流的汇集与演进 。支流入汇直接由各级径流叠加，水流演进过程通过断面判别模块调用扩散波、运动波法计算。基于 产流、汇流模型的径流模拟需输入流域降雨、地形几何参数、下垫面参数，参数率定是兼顾了模型实际意义的部分下垫面参数 。断面判别模块山区河流多不具备实测断面资料，仅能基于中、低分辨率的 提取有限几何参数，使得改进前的汇流模型对水流运动的模拟主要基于简化后对称的“”型断面，忽略了断面随河流级别由窄“”型至宽“”型的沿程变化。虽然原模型对尺度本就较小的低级别河流仍有较好的模拟效果 ，但为最大限度保留原模型对低级别河流的模拟效果，并使其能考虑高级别河流的河道特征，建立了 断面判别模块，自动考虑各级河流的断面形态。断面判别模块通过综合判别河流级别及断面形态，调用扩散波法 、运动波法 模拟不同级别河流径流（图 ）。与以往 基于“”型断面的扩散波法模拟水流演进不同，本研究建立的断面判别模块则主要基于概化断面的边界条件和适用于一般梯形断面的优化的运动波法 ，自动判别断面形态并调用对应的演进方程。研究区域及数据处理 研究区域概况皇甫川流域位于黄土高原北麓（，），面积为 ，于陕西省榆林市府谷县汇入黄河 。流域地处半干旱高原地区，砒砂岩大面积裸露，沙地面积大，土壤类型以栗钙土、黄绵土和风沙土为主 。流域内多年平均降雨量为 且多以短时强降雨为主，年均蒸发量达 ，年流域内主要分布有草地、农作物、灌木、沙地植被等 。流域内设皇甫、沙圪堵水文站和 个雨量站，分布见图 （），年皇甫川流域 和 分辨率的土地利用叠加数据见图 （）。研究数据及预处理为满足改进后模型需要（表 ），本研究采用了多种类型的数据，包括国产高分 号卫星（）、国产资源 号卫星（）影像，?号卫星影像、系列卫星影像、产品、数据、高分辨率 产品、国产卫星融合降雨数据等多源遥感数据。水文年鉴数据涵盖 年皇甫川流域皇甫站和沙圪堵站实测资料，以及用于概化河流断面的 年黄河中游紧邻于皇甫川、分布于 个流域的 个水文站的实测断面（图 （）。图 皇甫川流域基本情况 平滩流量及枯水流量确定河流断面几何信息提取及形态概化的关键是断面边坡坡度，该坡度可由最低、最高?平滩水位对应的河宽确定。这两类河宽资料往往需通过实测获取 ，十分匮乏。为了能通过遥感影像获取缺资料地区这两类河宽数据，需确定平滩及枯水流量日期。采用皇甫、沙圪堵两站 年的汛期径流及场次洪水资料，并验证资料一致性 。由于皇甫站断面不具备滩槽结构，难以直接确定平滩流量及日期，因此依据缺资料流域平滩流量估算方法及流域干支流关系 ，以沙圪堵站断面的平滩流量、日期代表皇甫川流域。确定水位为 左右的洪水流量为沙圪堵站平滩流量，流量为?，重现期约 年。最枯流量的确定，为保证一致性仍以沙圪堵站为准，为?。卫星遥感数据及预处理（）光学卫星影像获取及预处理：为通过遥感影像提取流域平滩、枯水河宽，选用能在山区更完整提取细小河流的高分辨率光学卫星影像，故联合使用 和 影像 。根据平滩日期，选择了 景 级 影像，成像时间为分别为 年 月 日、月 日，年 月 日、月 日。选择了 景 级 影像，成像时间分别为 年 月 日和 年 月 日。热红外数据在夏季能有效区分水体及其他地物 ，因此采集了同期的 景 产品。在枯水日期下，影像足够覆盖，直接使用 景 级影像，成像时间为 年 月 日、月 日。（）测高卫星数据：基于星载雷达观测的水位将用于提取缺资料地区河流断面形态。通过观测平滩、枯水河宽及其对应的水位信息，即可获取概化河道断面形态 。以?卫星产品为数据源，依据两类流量日期提取对应水位。经查询 年 月 年 月，卫星的第 号轨道位于流域上空，但受限于时空分辨率，选取与平滩流量相近的 年 月 日的 卫星数据，与枯水流量相近的 年 月 日的 卫星数据提取相应水位。土地利用数据及采样通过本团队开发的河宽自动提取算法 （）准确提取河流中心线及河宽，该算法因基于随机森林及神经网络算法，故需依赖准确的机器学习样本及合适的分类指数 。为确定研究区河流周围地物类型，基于土地利用类型数据 分析了 年皇甫川流域地物类型（图 （），结果显示 年流域内土地利用类型变化不大，确定了紧邻河流且易混淆的 类地物，遵循随机原则 完成了 余万学习样本的采集，样本像元总面积约占流域的 ，其中河流、农田耕地、山体林地、草场、城市的样本比例为 。年研究团队在研究区域内进行了 次野外考察（图 （），考察中记录下 个样本点的位置，测量了平滩河宽。因研究区域地物类型变化不大，故这 个实测样本可作为验证样本。数据及河网提取 河网作为 的结构核心以及河道边界条件的重要组成，其精度至关重要。低级别河流多因河宽极细，难以通过卫星产品提取 ，因此其断面概化、河宽提取均依赖高精度 数据辅助。数据分为以下几类：（）高分辨率 数据：团队于 年使用点位精度 、距离精度 、扫描密度 的徕卡三维激光扫描仪（），经过配准、点云处理等流程 生成了岔巴沟流域 数据。于 年、年联合使用机载激光扫仪与实时动态差分定位技术（）生成了皇甫川流域上游 分辨率的 和桥沟流域、油房渠流域 的 数据。（）基础 数据：用于提取黄河中游、皇甫川流域河网的 基础数据为 年公布的 版 ，该数据经过 、人工处理，高程精度优于 ，局部可达 。在纸坊沟流域使用由陕西省测绘局提供的 数据。基于 ，以清华大学开发的河网提取工具 （）提取黄河中游线状河流水系。经验证，能在提取流域水系时最大程度的消除平行错误，配合 使用提取结果更准 ，。使用 时流域最小汇流面积经预实验设定为 。融合降雨数据 需降雨数据驱动 ，以缺资料地区有限的雨量站插值估算流域降水进而模拟径流过程往往使模拟失真 。为解决这一问题，本研究采用了分布式卫星降雨融合数据 （），通过与皇甫川流域 年的实测资料比对，验证了 数据精度及可用性 。河流断面形态的概化以一般梯形等概化河流断面形态已有一定的研究基础 ，因此通过融合有限实测资料和多源遥感数据，建立了研究区各级河流断面概化模型：（）绘制 年皇甫川流域及与其紧邻流域共 个水文站的河流实测断面，按梯形断面概化并验证；（）将实测断面对应的最低河流级别作为断面概化模型中高、低级别河流的分界；（）基于研究区域及紧邻流域的高分辨率 ，提取了 个河流断面，按三角形（底宽为 的梯形）概化低级别河流断面并验证；（）基于多源遥感资料，提取了皇甫川流域的 个梯形断面。基于上述流程，完成了对皇甫川流域各级河流断面的概化。基于实测数据概化断面形态在黄河中游基于 数据共获得 级 河流，依此将断面分级。由高分辨率 所提断面知 级低级别河流断面近似“”型；水文站实测断面对应河流级别为 级，近似呈梯形。统计实测断面，综合比对断面所处级别后，按梯形（三角形为底宽为 的梯形）概化 ，。改进后的 模拟径流过程时，将按河流级别调用概化断面。概化后的断面几何参数见表 。表 水文站实测及高分辨率 概化断面几何参数级别位置左岸坡度 右岸坡度 级别位置左岸坡度 右岸坡度 桥沟流域油房渠流域纸坊沟流域岔巴沟流域纸坊沟流域岔巴沟流域岔巴沟流域皇甫川源区子长站青阳岔站李家河站贾家沟站杏河站枣园站马湖峪站曹坪站马王庙站 高家堡站高家川站申家湾站韩家峁站安塞站赵石窑站横山站殿市站转龙湾站王道恒塔站神木站温家川站延川站延安站甘谷驿站绥德站大宁站 皇甫站沙圪堵站高石崖站新庙站 丁家沟站白家川站河曲站府谷站 注：桥沟、油房渠、纸坊沟、岔巴沟流域及皇甫川源区断面为基于高分辨率 数据提取得到。基于遥感观测提取断面形态现有研究对缺资料河流断面的提取多基于河相关系 或由水槽实验推导出的经验公式 ，研究区域多为大江大河的下游，对山区河流、低级别河流的研究相对较少。皇甫川流域暂无相关研究且仅有 个实测断面，为更准确地概化各级断面，除参考实测资料外，还尝试基于遥感观测的河宽、水位提取河流断面。多源遥感观测的河道断面提取原理遥感观测断面形态，需提取断面关键几何参数，再依据图所示方法获取。由于研究区域内的低级别河流河道极窄，即便通过高分辨率卫星也难以提取 ，因此本方法适用于级别相对较高河流。基于多源遥感观测提取河流断面形态的步骤如下：（）基于 （）河流表面信息提取方法 及 河宽自动提取算法，提取流域高分辨率河流表面信息、河流中心线及平滩、最枯河宽。（）基于 算法 提取经河流表面信息去噪的研究区域?卫星的平滩及枯水水位。（）确定 卫星轨道与步骤（）提取的河流表面信息的交线，交线的端点坐标为 、且认为、处水位相等。读取该处水位 （图 ）。（）按水位递减顺序，从左岸至右岸依次连接端点得图 中最低水位以上断面。图 多源遥感观测的河流断面示意 河宽自动提取基于 及 法提取河宽的流程为：（）基于 法提取平滩、最枯流量下的河流表面信息。（）基于 算法，提取对应河流表面信息的河流中心线及水面边界并校正。（）基于校正后的河流中心线及水面边界，自动提取河宽。平滩河宽及枯水河宽提取结果基于 法，通过预实验确定的流域丰、枯水期能较好区分水体及其他地物的指数，依据这些指数提取流域内水体表面信息（包含河流、湖泊、水库等水体的初步产品）。如丰水期（月）选用了三价铁离子指数（）、归一化水体指数（）、增强水体指数（）、色彩渲染指数 （）种指数、种 热红外数据和 种灰度共生矩阵（）纹理指数（绿光、近红外波段均值，红光、近红外波段熵）提取平滩流量下的水体表面信息。枯水期使用形状指数（）、氧化铁指数（）、显色指数（）、及蓝光、绿光、红光波段的 种纹理指数提取最枯流量下的水体表面信息。水体表面信息提取完成后，需基于由流域 河流水系生产的缓冲区（半径 ）去除湖泊、水库、水塘等非河流水体和其他地物，再通过形态学算法、适当人工编辑后，得到精确的河流表面信息 。基于精确的河流表面信息，算法结合 的子流域划分结果共提取了 个平滩、枯水河宽数据对（图 （）。提取的平滩河宽通过 组野外实测数据及 组经高分辨率 断面算得的平滩河宽检验，枯水期最枯水位河宽则通过水文站及 历史影像检验。平滩河宽提取的、平均偏差 分别为 、，最小河宽对应 的 级河流；枯水河宽提取的、分别为 、，最小河宽对应 的 级河流。提取所得流域枯、平滩河宽之比为 ，实测值为 ，整体结果较准，提取结果的误差分析示意见图 （）。平滩、枯水河宽数据对将用以提取相应水位。平滩及枯水水位的自动提取提取 系列卫星 号轨道与皇甫川流域河流交线，据交线与河流边界交点确定水位提取点，假定每时期、每断面左右交点水位相等，概化河流断面。受限于 卫星的测量幅宽和过境时间等，在 年 月 日、年 月 日仅从皇甫川流域自动提取到 个断面，如图 （），断面关键参数见表 。图 河宽提取结果及检验表 遥感观测的河道断面断面编号 河流级别经度纬度平滩河宽（水位）?枯水河宽（水位）?左岸坡度 右岸坡度 级 （）（）级 （）（）级 （）（）皇甫川流域断面概化结果验证综合实测断面概化结果（表 ）及多源遥感断面概化结果（表 ）结果，将皇甫川流域 级河道的左、右岸坡度调整为 、，级河道调整为 、。概化后的各级断面与皇甫川流域实测断面对比见图 ，验证断面区域见图 。图 皇甫川流域断面概化结果示意以两岸坡度平均偏差定量评价，验证得 级低级别河道概化的坡度平均偏差 ，级河道概化边坡坡度平均偏差 ，故概化断面在一定程度上能较好反映断面真实形态。图 （）（）为全流域概化模型示意。将沙圪堵、皇甫站按本模型概化后（图 （）可知除皇甫站因设在乡镇，左岸受人为活动影响使坡角接近 外，概化断面与实测断面边坡较吻合，可供模型使用。概化模型由皇甫川及与其紧邻流域建立，概化断面的相关参数可供黄河中游上段其他无资料流域参考。因河流表面信息由 水系辅助提取，包含了 拓扑信息在内的所有参数，概化断面与河流表面信息的组合可提供流域河道连续的边界条件。改进的数字流域模型模拟结果 研究区域参数率定基于皇甫川流域河道边界条件，建立 运行所需数据库。模型率定基于 降雨数据、皇甫及沙圪堵站 年汛期逐日流量，主要评价指标为径流模拟纳什效率系数（）及 水文情报预报规范 中的洪峰流量相对误差、峰现时间误差。根据双层参数率定法 ，模型率定结果见表 。表 参数率定结果参数类型参数名称参数作用率定范围率定结果下垫面参数 表层土壤垂向饱和导水率?（?）表、下层土垂向饱和导水率?（?）表层土的水平渗透速率?（?）中层土的水平渗透速率?（?）表层土的田间持水量?（?）表层土的自由含水量?（?）中层土的田间持水量?（?）中层土的自由含水量?（?）表层土初始含水量?（?）以流域出口的皇甫站为例，率定期径流过程模拟结果如图 。改进后模型模拟 年汛期日径流的平均 为 ，从区可得改进后模型模拟日径流过程趋势更平缓、结果更准。模拟结果及验证以率定后参数模拟 年 月 日 日、月 日 日、月 日月 日的径流过程及 年洪水过程，仍以流域出口的皇甫站为例，模拟结果如图 （）（）所示，水力要素模拟较改进前有较大提升。径流模拟 达 ，且模型在区径流过程更接近实测过程（图 （），涨、落水趋势更平稳。进一步分析 年模拟所得断面水深 及平均流速 ，与实测数据验证率定期和模拟期所得水深 的平均误差为 ，平均测小 、为 、为 ；流速 模拟的平均误差为 ，平均测小?、为 、为 ，水力要素模拟结果均较准确。对洪水过程的模拟，以模拟难度较大的多峰洪水过程展示效果。模拟时间步长为匹配水文站洪水要素表的 。模拟洪峰流量相对误差分别为 、，均处在 水文情报预报规范规定的 范围内，模拟效果见图 ，改进后的模型能更真实模拟主洪峰过程中的子过程，且峰现时间与实际接近。模拟 年发生的共计 场洪水时，率定期 场洪水峰现时间较真实值平均提前约 ；模拟期则为 ，均处于 范围内。模拟的峰现时间整体提前可能源于两方面：（）卫星降雨融合产品在皇甫川流域仍存在一定的偏差；（）建立数字流域模型网格时密度仍然稀疏，导致汇流时间缩短。图 率定期夏汛期径流过程结果及对比图 年夏汛期径流过程及水力要素模拟图 典型洪水过程模拟效果及对比 结论与展望本文基于 河流表面信息提取方法及 河宽提取方法，提取了皇甫川流域 对平滩、最枯河宽数据对。结合有限的水文站实测资料及多源遥感数据，建立了皇甫川流域各级河流的断面概化模型，得到了流域各级较完整的河道边界条件。根据其包含的断面几何信息，建立了 断面判别模块，改进了 汇流模型。以 年为率定期率定 参数，模拟了 年皇甫川流域的汛期日径流过程、水力要素及年内的 场洪水过程。改进后的模型模拟效果提升明显，获得了更准确的水力要素。研究成果为水文模型改进提供了重要基础数据和参考，为模拟山区等缺资料地区河流的径流过程及洪水预报提供技术支持。取得的主要结论如下：（）建立了基于多源遥感数据的河道边界条件提取方法。该方法以传统线状 河流水系为基础，通过 、法自动提取流域内的河宽数据，结合 系列卫星提取的水位数据获得了河流断面。融合 个多源遥感提取断面、个水文站实测断面及 个高分辨率 断面，以一般梯形为基础建立了皇甫川流域 级河流的断面形态概化模型。经实测资料验证，最枯河宽提取的 和 分别为 、；断面概化模型的左右岸坡度与实测坡度的相对误差分别为 、。（）建立数字流域模型断面判别模块，该模块可根据河流级别自动识别河道断面形态，并调用对应的水流演进方程综合模拟流域内各级河流的径流过程。对 年日径流过程模拟及验证的 均优于 ；流速模拟的相对误差、分别为 、；水深模拟的相对误差、分别为 、。（）对 年 场典型多峰洪水的模拟结果表明，改进后的模型模拟洪峰流量的相对误差为 ，峰现时间误差约为 。整个率定、模拟期，场洪水的洪峰流量相对误差约 ，峰现时间提前约 。虽然本研究改进了 的汇流模型，提升了径流、洪水过程及水力要素的模拟效果，但未能联合不易受云雨干扰的雷达卫星，暂未实现对河流的全天候观测、提取。对断面的概化尚基于较单一的一般梯形，暂未考虑更复杂的复式断面及最枯水位以下形态，较适用于河流的中上游段，对最枯水位下仍存在较大水深的河道断面适用性有待验证。对 的改进暂引入了较简单的运动波法，还应继续尝试建立更适合体现山区河流洪水及其输沙特征的计算方法。参考文献：龙笛，李雪莹，李兴东，等遥感反演 年青藏高原水储量变化及其驱动机制 水科学进展，（）：李丹，吴保生，陈博伟，等基于卫星遥感的水体信息提取研究进展与展望 清华大学学报（自然科学版），（）：岩腊，龙笛，白亮亮，等基于多源信息的水资源立体监测研究综述 遥感学报，（）：胡春宏，郑春苗，王光谦，等“西南河流源区径流变化和适应性利用”重大研究计划进展综述 水科学进展，（）：姚檀栋，王伟财，安宝晟，等 年青藏高原科学考察研究历程 地理学报，（）：徐宗学，程磊分布式水文模型研究与应用进展 水利学报，（）：王光谦，刘家宏黄河数字流域模型 水利水电技术，（）：李铁键流域泥沙动力学机理与过程模拟 北京：清华大学，刘家宏黄河数字流域模型 北京：清华大学，张昂，李铁键，傅汪，等基于卫星降水融合产品的黄河源区汛期径流模拟 应用基础与工程科学学报，（）：王晨沣，傅旭东，张?，等黄土高原植被作用下黄河数字流域模型坡面侵蚀模块改进 清华大学学报（自然科学版），（）：陈瑾晗，夏军强，李铁键，等数字流域模型在长江上游龙河流域的应用 水力发电学报，（）：罗火钱，侯才水，马广恩 闽江沙溪数字流域模型建立及其应用研究 水利科技与经济，（）：尹小玲基于数字流域模型的珠江补淡压咸水库调度研究 北京：清华大学，王兴奎，王光谦，吴保生，等河流动力学 北京：科学出版社，王光谦，李铁键，贺莉，等黄土丘陵沟壑区沟道的水沙运动模拟 泥沙研究，（）：宋凤军皇甫川流域植被减水减沙效应的模拟研究 杨凌：西北农林科技大学，杨大文，李罛，倪广恒，等分布式水文模型在黄河流域的应用 地理学报，（）：，：，（）：薛源，覃超，吴保生，等基于多源国产高分辨率遥感影像的山区河流信息自动提取 清华大学学报（自然科学版），（）：，（）：覃超，吴保生，汪舸，等青藏高原山区河流广义河相关系与多频率沿程河相关系 水利学报，（）：，（）：李玲玲若干典型地物温度变化特征试验研究 沈阳：东北大学，（?）：，（）：，（）：，：，（）：，：，（）：，（）：，（）：，：田凤云，吴志勇，徐征光，等黄土高原典型流域网格化逐时降水融合产品 的精度分析 水电能源科学，（）：，：，：胡春宏，陈建国，刘大滨，等水沙变异条件下黄河下游河道横断面形态特征研究 水利学报，（）：于红波改进随机森林算法在北方河流测流断面形态预测中的应用 水利技术监督，（）：黄琦基于多源遥感观测的西南缺资料流域河道径流量反演研究 北京：清华大学，（）：，（?）：中国国家标准化管理委员会水文情报预报规范：北京：中国标准出版社，（）：，（，；，；，；，）：，（），：；（责任编辑：于福亮）