改进垂向混合产流模型在营口洪水模拟中的应用

1 研究流域概况

近些年来，随着极端气候天气频发，营口地区的暴雨洪水频率也逐步增多，为有效应对洪水，需要对洪水进行精准预报，提高区域防洪决策的科学性。营口地区常用的洪水预报方案为降雨径流经验关系，这种方式需要采用较长的实测数据系列，且需要对降雨径流经验关系进行多次修订［1-2］。当前，水文模型的多样化已逐步为区域洪水的精准预报提供了较好的模拟平台［3］，这其中改进的垂向混合产流模型由于垂直方向可综合考虑超渗和蓄满产流机制，在北方地区适用性较好［4-9］，但在营口地区还未得到相关应用，本文结合改进的垂向混合产流模型，以营口地区典型水文站熊岳站15场实测洪水数据，分析该模型的适用性。研究成果对于营口地区洪水预报方案制定具有重要参考价值。

本文以熊岳站为研究实例，熊岳站位于熊岳河上，熊岳河位于盖州市熊岳镇。上游分南北两支，南支（主支）发源于杨运乡老帽山，北支发源于老平顶，河长44.2km，流域面积352km2，河道平均比降0.315%，熊岳水文站断面以上面积307km2，河长南支29.2km、北支28.2km。该流域植被良好，包气带较薄，地下水埋藏较深，熊岳镇东南的熊岳河谷有丰富的地下温泉。该流域处于中纬度地区，属暖温带半湿润季风性气候，四季分明。春季少雨多风，夏季高温多雨，秋季天高气爽，冬季寒冷干燥。本流域暴雨多数为西风带系统造成。本流域暴雨中心往往出现在上游八道河一带，八道河24 h最大降雨量289.9mm （1981年7月27日），本流域年最大降雨量1088.9mm（1964年），最小降雨量355.4mm（1944年）多年平均降雨量为702mm。

2 改进垂向混合产流模型原理

垂向混合产流中蓄满、超渗的流域面积比例是随前期土壤含水量和下渗水量的变化而改变，如式（1）：

“孩子，你身上带着一股不祥的气息，这股气息正在朝着你的内心侵蚀，让你变成邪恶的魔鬼！”天葬师的语气骤然变得激烈，一字一字击打在青辰的胸膛，“这些天，你究竟遇到了什么！”

式中 α为蓄满产流的面积比例系数； FA为实际下渗量（mm）；Wmm 为流域最大蓄水量（mm）； α为相应于初始土壤平均含水量为W时的纵坐标值；B为流域蓄水容量分布曲线指数。采用具有流域分布特征的格林—安普特下渗曲线，计算如式（2）～式（4）：

结合熊岳站1964～2017年15场典型洪水，对比改进垂向混合产流模型在营口地区洪水模拟精度，结果如表2，并选取6场典型洪水分析改进前后垂向混合产流模型洪水过程的模拟结果，如图2。

在产流计算的基础上，对其水源进行划分，进入地面以下的水流RR首先补充自由水蓄水量S（mm），然后按照出流系数对壤中流和地下径流进行划分，划分方程如式（3）～式（5）：

式中 FM为流域平均下渗能力 （mm）；FC为稳定下渗率（mm/s）；WM为流域平均蓄水容量（mm）；W为流域实际土壤含水量（mm）；KF为土壤缺水量对下渗率影响的灵敏系数；BF为反映下渗能力空间分布特征的参数；PE为扣除雨间蒸发的降雨（mm）；RS为地面径流（mm）。

地面径流：RS=PE-FA

传统垂向混合产流在水源划分时地面以下径流主要划分为壤中流和地下径流两种水源，而由于在干旱半干旱地区，其土壤蓄水容量一般较大，且由于降雨和下垫面时空分布的不均匀性，土壤缺水量较小的区域遭遇一定量级的降水后其较易形成饱和地面径流。为此运用三水源新安江结构对垂向混合产流模型的水源结构进行改进，划分为饱和地表、壤中及地下径流3种水源，如图1。

式中 St为t时段达到自由水的蓄水量（mm）；St-1为t-1时段达到自由水的蓄水量（mm）；RIt及RGt分别为t时段壤中流及地下径流；KI，KG分别为两种水源的出流系数。

结合熊岳站1964～2017年实测径流资料，对改进模型参数进行率定，参数率定结果如表1。

图1 新安江三水源划分结构

3 模型应用

社区大学是社会参与力的培植及养成场所，所有的教学与活动内容紧扣着社区及地方发展的需求。而根据教师的需求及兴趣，可以组建不同的专业学习社群，教师可以随时开展专业经验交流及分享讨论。社区大学的课程开发和设计，教师都要付出长时间的精力和心血才能完成，靠教师个体往往很难达成目标，可通过各专业学习社群的经验交流分享，提升其社区专业知识技能，共同打造社区大学的特色课程。社区大学的教师若能积极参与专业学习社群，形成学习型组织，彼此分享教学经验和教学成长故事，将有助于社区大学教师的教学专业发展［5］。

表1 改进模型参数率定成果

参数符号 参数意义 参数值Kc 蒸发折算系数 000.85 UM 上层张力水容量/mm 030 LM 下层张力水容量/mm 070 C 深层蒸散发折算系数 000.22 F0 初始下渗率 022.3 FC 稳定下渗率 021.2 WM 流域平均张力水容量/mm 100 BF 流域下渗曲线分布方次 000.27 B 张力水蓄水容量曲线方次 000.43 IM 不透水面积占全流域面积的比例 000.05 SM 表层自由水蓄水容量/mm 050 EX 表层自由水蓄水容量曲线方次 000.74 KG 表层自由水蓄水库对地下水的日出流系数 000.38 KI 表层自由水蓄水库对壤中流的日出流系数 000.19 CI 壤中流消退系数 000.61 GG 地下水消退系数 000.94 L 河网滞时/h 003.54 XE 马斯京根法演算参数 000.53 KE 马斯京根法演算参数/h 001.72 N 马斯京根划分河段数目 014 IM 不透水面积所占比例 000.35

实际下渗量：

表2 改进前后垂向混合产流模型模拟精度对比

序号模型 原模型 改进模型000 0.353 0.672 -22.3 14.5 25.3 14.8 -6 -2 001 0.451 0.771 -18.5 15.6 28.5 15.2 -4 -2 002 0.352 0.706 -19.2 -17.5 -29.6 -16.5 -3 2 003 0.289 0.713 -21.5 -18.9 -31.5 13.8 -2 0 004 0.461 0.725 -28.5 10.5 22.3 12.9 -5 -3 005 0.551 0.743 22.9 11.6 24.5 14.3 2 -4 006 0.452 0.752 -20.2 -14.5 25.6 -15.2 -1 -6确定性系数 径流深相对误差/% 洪峰误差/% 峰现时间误差/h原模型 改进模型 原模型 改进模型 原模型 改进

续表2

序号模型 原模型 改进模型007 0.383 0.775 -21.5 -13.5 27.5 -16.5 -6 2 008 0.396 0.718 22.7 -15.3 22.3 -18.3 -5 -1 009 0.413 0.737 23.5 10.5 -21.5 -14.5 -6 1 010 0.405 0.743 24.3 -11.5 22.9 16.2 -5 0 011 0.514 0.756 -26.5 -12.9 -24.5 15.3 -4 -2 012 0.554 0.743 -37.5 13.4 -26.3 14.5 -3 0 013 0.525 0.731 22.5 -15.2 -28.5 15.9 -2 -1 014 0.483 0.774 21.5 14.3 -24.7 12.8 2 2 015 0.515 0.752 -24.5 15.7 26.5 13.7 -2 -5确定性系数 径流深相对误差/% 洪峰误差/% 峰现时间误差/h原模型 改进模型 原模型 改进模型 原模型 改进平均值或合格率 0.443 0.738 23.6 14.1 25.8 15.0 40% 80%

图2 垂向混合产流模型典型洪水模拟

从表2可看出，改进后垂向混合产流相比于改进前，在熊岳站洪水模拟各项评价指标中都有明显改善，其中确定性系数是表示洪水过程模拟精度的高低。改进的垂向混合产流模型的15场洪水确定性系数均值为0.738，相比于原模型提高0.295，在洪峰和洪量误差上，改进模型也较原模型有较为明显的改善，原模型洪峰和洪量误差均在20%以上，而改进模型在20%以内，相比于原模型，改进模型洪峰和洪量均值分别提高9.5%和10.8%。从峰现时间合格率上，改进的垂向混合产流模型也较原模型有明显改善，峰现时间合格率提高40%。这主要是改进的垂向混合模型综合考虑饱和地面径流，较原模型在水源结构上有所改善，更适合于北方地区有降雨径流的变化特性。从图2中选取的6场典型洪水可看出，各场次洪水实测值和模拟值的吻合度都较高，改进的垂向混合产流模型适合于营口地区的洪水模拟。

会议要求，整县提升合作社发展质量要以贯彻落实农民专业合作社法为主线，坚持规范与创新并举，强化指导扶持服务，树立合作社高质量发展的县域样板。一要增强单体合作社自身能力，强化规范建设，拓展业务范围，丰富出资方式，增强服务功能。二要培育发展农民合作社联合社，推进区域性联合、行业性联合、产业链联合、多主体融合。三要提升县域指导服务能力，健全规章制度，构建扶持政策体系，加强示范引领，探索准入和退出机制，搭建合作社服务平台。

4 结语

（1） 改进的垂向混合产流模型由于综合考虑饱和地面径流，改善传统模型的水源划分结果，在熊岳站模拟精度都较原模型有较大程度改善，更适用于营口地区的洪水模拟。

（2） 由于本文选取站点较少，在以后研究中还需引入更多的站点，从而充分论证改进垂向混合产流模型在营口地区洪水模拟的适用性。

针对临床中患有食管癌的老年患者治疗方式选择来讲,给予患者管状胃治疗的效果较好,患者不良反应发生几率较低,且生活质量会得到显著提高,该方案可作为此类患者的首选治疗方案进行推广。

参考文献：

［1］杨大卓.营口地区地表水可利用量估算［A］.2011（第五届）水业高级技术论坛论文集［C］.

［2］郭晓亮.感潮河段潮位与洪水预报及上游水库控制模式研究［D］.大连：大连理工大学，2009.

［3］聂大鹏，秦秀梅，薛丽.营口地区水资源现状分析［J］.东北水利水电，2013，31（3）：42-44.

［4］杨伟，赵天宇，杜彦臻，毕庶刚，林洪孝，王刚.基于改进的垂向混合产流模型的水资源评价模型研究及应用 ［J］.中国农村水利水电，2018（8）：95-99.

［5］李大洋，梁忠民，侯博，李彬权，王军.垂向混合模型在湫水河洪水预报预警中的应用［J］.人民黄河，2018，40（6）：24-28.

［6］韩丽娜.改进的垂向混合产流模型在辽西干旱半干旱区域洪水模拟中的应用［J］.吉林水利，2017（5）：44-48.

［7］毛北平.垂向混合产流模型在无资料地区山洪灾害临界雨量计算中的应用［J］.应用基础与工程科学学报，2016，24（4）：720-730.

［8］闻绍珂，蔡涛.改进的垂向混合产流模型在大凌河流域的运用［J］.东北水利水电，2015，33（4）：45-48，72.

［9］王庆平，沈国华，王红艳.垂向混合产流模型在不同地区的应用与改进［J］.节水灌溉，2012（5）：11-15.