黄河源区玛曲至唐乃亥区间水文模型适用性研究.pdf

收稿日期：基金项目：国家自然科学基金黄河水科学研究联合基金资助项目（）；青海黄河上游水电开发有限责任公司科技项目（）作者简介：胡一栋（），男，陕西柞水人，高级工程师，硕士，主要从事水电运行管理和水库调度工作 通信作者：史玉品（），女，山东成武人，正高级工程师，硕士，主要从事洪水分析及预报工作：【水文泥沙】黄河源区玛曲至唐乃亥区间水文模型适用性研究胡一栋，史玉品，康耀辉，沈延青，刘玉环（青海黄河上游水电开发有限责任公司，青海 西宁；黄河水利委员会 水文局，河南 郑州）摘 要：黄河源区是黄河流域重要产流区，该区域的洪水预报对龙羊峡、刘家峡等黄河上游梯级水库群的防洪调度，乃至整个黄河流域的水旱灾害防御具有举足轻重的作用。以黄河源区主要产流区玛曲至唐乃亥区间为研究对象，研究了萨克拉门托模型、陆浑模型及垂向混合产流模型在该区域的适用性，结果表明：种模型模拟的洪峰流量、次洪水量的合格率均大于，确定系数均值都在 以上，各模型都能较好地模拟该区域的洪水过程；种模型均不能对所有场次洪水的模拟精度达到满意的效果，因此黄河源区的洪水作业预报应采用多模型组合的集合预报方法。关键词：萨克拉门托模型；陆浑模型；垂向混合产流模型；黄河源区中图分类号：；文献标志码：引用格式：胡一栋，史玉品，康耀辉，等黄河源区玛曲至唐乃亥区间水文模型适用性研究人民黄河，（）：，（，；，）：，：；水文模型分析是水文规律研究、洪水分析及洪水预报的主要技术手段。随着计算机技术、地理信息系统和遥感技术的发展，流域水文模拟经历了由系统模型模拟到概念性模型模拟再到物理模型模拟、由集总式模拟到分布式模拟的发展历程。不同水文模型结构、原理及参数意义不同，适用区域及条件也不同，在湿润半湿润地区蓄满产流模型的应用效果较好，而干旱半干旱地区超渗产流模型的适用性较强。受人类活动等因素影响，流域内植被、下垫面及河道条件不断变化，且每次降水产流过程也有一定差异，选取单一水文模型对多场次洪水进行模拟预报时，有些场次洪水的模拟精度较高，而某些场次洪水模拟效果不太理想。为了在实际洪水预报中减小模拟误差，应采用多模型组合的集合预报方法。黄河源区是黄河流域重要产水区和水源涵养区，蕴藏着丰富的水资源。唐乃亥水文站既是黄河源区出口控制站，也是黄河上游多年调节性水库龙羊峡水库入库站。如果黄河上游汛期发生长历时、大流量洪水过程，在没有精准洪水预报为水库调度提供技术支撑的情况下，那么龙羊峡水库可能出现大流量泄洪、弃水，既影响黄河中下游防汛工作，也会因弃水而造成大的经济损失，因此唐乃亥水文站洪水精确预报对龙羊峡、刘家峡等黄河上游梯级水库群的防洪调度乃至整个黄河流域水旱灾害防御具有举足轻重的作用。第 卷第 期 人 民 黄 河 ，年 月 ，流域概况黄河源区是指黄河上游唐乃亥水文站以上区域，为典型的高原大陆性气候区。本文研究区域为黄河源区重要产流区玛曲至唐乃亥之间的区域，集水面积为 万，流域地形地貌以山地、丘陵和盆地为主，土地覆盖以草地、灌木和裸土为主，河道切割较深，水系发达，较大支流有泽曲、切木曲、大河坝河等。该区域处于东亚季风区边缘，西太平洋副热带高压的强度和位置变化导致东南暖湿气流对黄河源区水分的补充有很大不确定性，致使黄河源区降水在年度和季节尺度上的振荡较为明显。西太平洋副热带高压强大西伸，且停滞少动，使青藏高原上空偏南气流明显增强，加之孟加拉湾暖湿气流不断北上，若遇冷空气自北南下，则往往形成持续性强连阴雨天气。受此影响，黄河源区洪水具有持续时间长、洪量大、缓涨缓落等特性。玛曲至唐乃亥区间流域示意见图。图 玛曲至唐乃亥区间流域示意 研究方法针对黄河源区复杂的下垫面条件及气候特点，选用考虑蓄满和超渗两种产流机制的萨克拉门托模型、陆浑模型及垂向混合产流模型，对黄河源区主要产流区玛曲唐乃亥的洪水过程进行模拟，并对模拟结果进行对比分析，从而优选模拟效果较好的模型，以便为该区域洪水分析与预报、水资源优化配置及上游梯级水库调度服务。萨克拉门托模型是一个集总性的多参数模型，具有超渗和蓄满产流特性。垂向混合产流模型是超渗产流和蓄满产流在垂向上组合的模型。由于这两种模型的应用范围比较广，因此本文对于模型结构及原理不再详细介绍。陆浑模型是黄河水利委员会水文局研发的兼有超渗与蓄满两种产流机制的经验水文模型，该模型最早试用于伊河上游陆浑水库的入库洪水预报，因此命名为陆浑模型，该模型在黄河三门峡至花园口区间应用效果较好。下渗曲线采用霍顿下渗能力曲线，下渗能力面积分配曲线采用抛物线型，其公式形式为 （）（）（）（）（）（）（）（）（）式中：为研究时段径流深，为研究时段降水量，为下渗能力消退系数，为计算时长，为流域平均最大蓄水量，为研究时段流域累计平均蓄水量，为抛物线指数，为下渗能力分配曲线上与 相对应的纵坐标。三种模型的汇流均采用三水源滞后演算法及马斯京根分段连续演算法。模型适用性分析整理分析玛曲至唐乃亥区间 年 个雨量站（年以后）、个水文站的气象、水文资料，分别运用萨克拉门托模型、陆浑模型及垂向混合产流模型模拟研究区域洪水过程，采用人工试错和自动优选两种完全耦合方式率定 种模型的最优参数，采用专家经验法对参数的敏感性进行分析，找出模型需要优选的参数。以洪峰流量、洪量（径流量）及确定系数为评价指标分析 种模型的模拟效果，模拟结果的精度评定根据水文情报预报规范（）进行。参数率定利用收集的气象、水文资料进行萨克拉门托模型、陆浑模型及垂向混合产流模型的参数率定。种模型的流域汇流参数为地表水线性水库汇流系数、壤中流线性水库汇流系数、地下水线性水库汇流系数、河网汇流的滞后时间，马斯京根分段连续演算参数为分段马斯京根法演算参数、马斯京根法流量比重系数、河道分段数。虽然 种模型采用相同的汇流结构，但水源结构的划分不同，计算产流量也不同，因此汇流参数的取值也不同。河网汇流滞后时间 根据经验取值，河道分段数 依据洪水传播时间取值，其余 个参数均为敏感性参数，需要进行率定优选，优选结果见表。人 民 黄 河 年第 期表 种模型的汇流参数优选结果参数萨克拉门托模型陆浑模型垂向混合产流模型 萨克拉门托模型其他参数率定结果见表。萨克拉门托模型参数较多，其中不敏感的参数可以根据经验、物理特征给出参数值，也可以借用黄河流域其他相似地区已率定好的参数值。比如直接进入下层自由水的渗透水比例 一般取值为；下层自由水中不蒸发的比例；对于不闭合的地面水出流，沿河槽损失的流量 ；不闭合地下水出流与河道基流比例 。模型中对径流总量影响敏感的参数有、，对透水地面上的径流影响敏感的参数有、，对不透水地面上的径流影响敏感的参数有、，其中可以估计初始值的参数有、。萨克拉门托模型中有些参数是互相影响的，在率定优选参数时必须综合考虑，将参数值控制在合理的范围内，避免模拟结果出现虚假精度。表 萨克拉门托模型参数率定结果参数参数含义取值不透水地面占全流域的比例变化的不透水地面占全流域的比例河网、湖泊、水库及水生植物面积占比上层张力水最大值 上层自由水最大值 上层自由水日出流系数最大下渗率渗透指数直接进入下层自由水的渗透水比例下层张力水最大值 下层浅层自由水最大值 快速地下水日出流系数下层深层自由水最大值 下层深层自由水日出流系数下层自由水中不蒸发的比例沿河槽损失的流量 （）不闭合地下水出流与河道基流比例 陆浑模型参数有流域平均最大蓄水量，下渗强度消退系数，下渗能力面积分配曲线指数。相对不敏感，没有对其优选，根据黄河源区的气候及下垫面条件取 。、属于敏感参数，经过参数率定优选确定其值分别为、。垂向混合产流模型参数率定结果见表。根据专家经验，垂向混合产流模型中、是不敏感参数，一般根据经验及物理意义进行取值。参数 与流域特性密切相关，取值要保证土壤含水量 的合理性；参数 反映流域面上蓄水容量分布的不均匀性，取值范围一般为；参数 反映流域下渗能力分布的不均匀性，取值范围一般为；深层蒸散发折算系数 对模拟结果的影响与流域的湿润程度有关，在半干旱半湿润地区较敏感，因此对参数 的取值进行优选率定。模型中敏感参数还有、，其中 表示自由水出流的快慢，一般取 ，因此该模型需要优选率定的参数有、。表 垂向混合产流模型参数率定结果参数参数含义取值流域蒸散发折算系数流域平均蓄水容量 流域上层蓄水容量 流域下层蓄水容量 深层蒸散发折算系数流域蓄水容量分布曲线指数表层土自由水蓄水容量 土壤对下渗率影响系数流域下渗率分布曲线指数自由水箱地下水出流系数自由水箱壤中流出流系数流域实际下渗率（）模拟结果分析 种模型对黄河源区 年 场次洪水的模拟效果见图 图。）萨克拉门托模型。由洪水过程模拟结果可知，场次洪水中洪峰流量相对误差超过 的有 场次，合格率为，相对误差最大值为、平均值为。次洪水量相对误差在 以内，合格率为；次洪水量相对误差超过 的共有 场次，占；次洪水量相对误差最大值为，平均值为。确定系数 场次达到，最大值为，最小值为，平均值为，该模型的模拟结果与实测洪水过程拟合较好。综合分析，洪水模拟的合格率为，确定系数均值为，按水文情报预报规范（）精度评定标准，该模型的模拟精度达到乙等。）陆浑模型。由洪水过程模拟结果可知，场次洪水中洪峰流量相对误差超过 的有 场次，合格率为 ，相对误差最大值为 、平均值为。次洪水量有 场次的相对误差超过，合人 民 黄 河 年第 期格率为；次洪水量相对误差超过 的共有 场次，占；次洪水量相对误差最大值为，平均值为。确定系数 场次达到，最大值为，最小值为，平均值为，该模型的模拟结果与实测洪水过程拟合较好。综合分析，洪水模拟的合格率为，确定系数均值为，按水文情报预报规范（）精度评定标准，该模型的模拟精度达到乙等。图 洪峰流量相对误差图 次洪水量相对误差图 确定系数）垂向混合产流模型。由洪水过程模拟结果可知，场次洪水中洪峰流量相对误差超过 的有 场次，合格率为，相对误差最大值为、平均值为。次洪水量相对误差均在 以内，合格率为；场次洪水次洪水量相对误差超过，占；次洪水量相对误差最大值为，平均值为。确定系数 场次洪水达到，最大值为，最小值为，平均值为，模型模拟过程与实测洪水过程拟合较好。综合分析，洪水模拟的合格率为，确定系数均值为，按水文情报预报规范（）精度评定标准，该模型的模拟精度达到乙等。由以上分析可知，萨克拉门托模型、陆浑模型及垂向混合产流模型在玛曲至唐乃亥区间模拟精度都能满足洪水预报精度要求，尤其是对次洪水量的模拟计算，除陆浑模型模拟结果有 场次相对误差超过 外，其他两个模型模拟结果的相对误差均在 以内。陆浑模型的确定系数均值最小，由图 可以看出，陆浑模型的确定系数比其他两种模型的总体偏小。典型洪水模拟结果分析受连阴雨天气影响，黄河源区洪水过程一般是震荡性缓涨缓落，以单峰型和复式峰型居多。分别选择复式峰型的洪水过程与单峰型洪水过程进行洪水模拟结果分析。其中，年 种模型模拟的洪峰流量相对误差均较大，因此单峰型洪水过程选择 年场次洪水。同时选择 年复式峰型洪水过程分析 种模型的模拟结果。年洪水过程模拟结果见图。种模型模拟的唐乃亥水文站洪峰流量误差较大，均偏小 以上。影响模型模拟结果的因素有很多，比如模型结构、参数、降水及流域下垫面条件等，其中降水是影响水文模型模拟结果最主要、最直接的因素。由图 的洪水过程线可以看出，峰前涨水过程与峰后退水过程拟合较好，但洪峰部分偏差较大。该场次洪水玛曲水文站实测洪峰流量为 ，唐乃亥水文站洪峰流量为 ，玛曲至唐乃亥区间加水产生的最大流量为 ，而模型模拟结果误差最大达，最大流量差为 ，占区间来水最大流量的，因此此次模拟结果偏小的主要影响因素是降水。年玛曲至唐乃亥区间 万 的流域上只有玛曲、军功、唐乃亥、东倾沟 个雨量站的降水资料，面积大、历时长、强度小的连阴雨天气，对模型模拟结果的影响相对较小，一旦出现离雨量站较远的局部强降水过程，雨量站就无法监测到真实的降水资料，计算的降水量误差较大，造成洪水预报结果出现系统偏差。图 单峰型洪水过程模拟结果人 民 黄 河 年第 期 年洪水过程模拟结果见图，由图 可以看出，年洪峰流量最接近实测值的为萨克拉门托模型的模拟结果，其模拟的两个洪峰流量的相对误差分别为、；垂向混合产流模型的洪峰流量模拟结果最小，洪峰流量相对误差分别为 和。涨水过程 种模型模拟结果与实测值拟合较好，退水过程垂向混合产流模型与萨克拉门托模型的模拟结果小于实测值；陆浑模型的确定系数为，图 复式峰型洪水过程模拟结果其他两个模型的分别为、，说明陆浑模型模拟的该场次洪水过程更接近实测过程。年汛期模型应用 年汛期黄河源区出现了长历时、大流量洪水过程，唐乃亥站分别在 年 月 日及 月 日出现 次编号洪水，洪峰流量分别为 、，月 日至 月 日流量在 以上历时达（期间仅 月 日、日 日均流量为 ），最大 洪量为 亿，月总径流量 亿、较多年均值偏多。利用已经率定好参数的 种模型，模拟 年汛期唐乃亥站 次编号洪水，号洪水和 号洪水是一个连续过程，洪水过程线不再分开显示，以 号洪水退水至最低点（月 日 时）为节点计算 次编号洪水的洪峰流量、洪量误差与确定系数等，结果见表、图。表 年两次编号洪水模拟结果模型洪水编号洪峰流量（）实测值计算值洪峰流量相对误差次洪水量 亿 实测值计算值次洪水量相对误差确定系数垂向混合产流模型 号洪水 号洪水 萨克拉门托模型 号洪水 号洪水 陆浑模型 号洪水 号洪水 图 年 号、号洪水过程模拟结果 对 年唐乃亥站 次编号洪水模拟结果显示，垂向混合产流模型的模拟结果最优，洪峰流量相对误差小于，次洪水量相对误差小于，确定系数在以上；陆浑模型与萨克拉门托模型的模拟结果相对差一些，尤其是 号洪水，萨克拉门托模型拟合效果相对较差，确定系数为，陆浑模型模拟 号洪水的洪峰流量比实测值偏小。结论以玛曲至唐乃亥区间 年的气象、水文资料为基础，选择萨克拉门托模型、陆浑模型以及垂向混合产流模型模拟唐乃亥站的 场次洪水过程，并定量分析了模型的模拟效果，次洪水量的合格率 种模型均达到，洪峰流量的合格率均达到，确定系数均超过，依据洪水预报精度评定标准，种模型都可用于该区域的洪水作业预报。不同模型在黄河源区的洪水模拟中表现出不同的效果，在实际洪水预报中为了减小模拟误差，应采用多模型组合的集合预报方法。降水是影响洪水预报精度最主要、最直接的因素，但是黄河源区复杂的高山地形及高寒气候条件，使得雨量计的安装与维护极为困难，该区域 万 的流域面积目前只有 个雨量站，离水文预报要求的 多个相差较大。虽然黄河源区的降水以大面积、低强度、长历时的连阴雨为主，但也会发生区域性强降水，如果因雨量站的稀少而没有监测到强降水信息，模型的模拟结果就会出现系统性偏差。参考文献：朱炬明，周买春不同水文模型在双桥流域的应用比较人民黄河，（）：汪妮，解建仓，冯黎，等面向对象技术在流域洪水预报系统中的应用西安理工大学学报，（）：（下转第 页）人 民 黄 河 年第 期 ，：，：，（）：，：，：，（）：，：，：，：，（）：马承新，邓海燕，王维平山东省与其他区域间虚拟水贸易定量研究灌溉排水学报，（）：马承新，王维平，刘品山东省虚拟水贸易问题研究：基于投入产出分析方法 中国农村水利水电，（）：李凤丽，曲士松，王维平，等 年山东省虚拟水贸易变化及典型区生态环境响应灌溉排水学报，（）：【责任编辑 张华兴】（上接第 页）黄志强，李柏宏，凌仙华水文模型组合预报应用研究浙江水利水电专科学校学报，（）：雍斌，张建云，王国庆黄河源区水文预报的关键科学问题水科学进展，（）：马明卫，程玉佳，崔惠娟，等黄河唐乃亥站径流干旱多变量联合特征分析华北水利水电大学学报（自然科学版），（）：赵芳芳气候变化对黄河源区径流量的影响分析北京：北京师范大学，：熊剑锋，张树全萨克拉门托水文模型在天生桥洪水预报的应用云南水力发电，（）：包为民，王从良垂向混合产流模型及应用水文，（）：张芳珠，马卫东，李作安，等降雨径流模型在故县水库入库预报中的应用研究水资源与水工程学报，（）：瞿思敏，包为民，张明，等新安江模型与垂向混合产流模型的比较河海大学学报（自然科学版），（）：张强昌江洪水预报模型研究南昌：南昌大学，：刘金平，乐嘉祥萨克拉门托模型参数初值分析方法研究水科学进展，（）：【责任编辑 吕艳梅】人 民 黄 河 年第 期