基于参数移植法的无资料地区洪水预报分析

基于参数移植法的无资料地区洪水预报分析
赵芳; 华东; 田质胜; 唐克银; 张晨
【期刊名称】《《水科学与工程技术》》
【年(卷),期】2019(000)005
【总页数】4页(P29-32)
【关键词】参数移植; 无资料地区; 洪水预报; 新安江模型; SCE-UA算法
【作者】赵芳; 华东; 田质胜; 唐克银; 张晨
【作者单位】山东省水利勘测设计院济南 250014; 水利部信息中心北京100032; 国家水资源监控能力建设项目办公室北京100032
【正文语种】中文
【中图分类】TV122
无资料地区的洪水预报，对于区域防洪和水利工作的开展具有重大意义。

目前针对无资料地区，常用洪水预报方法主要有：①参数移植法，采取临近有资料流域的参数移植到无资料地区；②模型法，建立具有物理机制的水文预报模型推广至无资料地区；③3S技术，获取无资料地区空间地理信息，进行三维建模预测。

这3种方
法均有各自优势，张漫莉［1］采取模型法构建了武江流域新安江洪水预报模型，并成功应用到无资料子流域，预报精度达乙级以上；宋昭义［2］采取参数移植法，利用冷水河临近流域的水文资料构建了新安江模型，最终成功运用在冷水河流域，预报精度为丙级；冯娇娇［3］基于DEM数据和遥感数据，利用GIS提取流域特
征建立辽宁东部短缺资料地区的新安江模型。

尽管模型法和3S技术比参数移植法精度和准确性更高，但建立具有物理机制的水文预报模型往往需要长久测算，而3S技术成本较高且技术性强，因此，目前较多用的仍是参数移植法。

基于此，本文以浑河上游清原下水库流域为研究区，探究了参数移植法的可行性和有效性。

1 研究思路
本文以缺少资料的浑河上游清原下水库流域为研究区，探索参数移植法在无资料地区的适用性和可行性。

根据流域下垫面属性存在的差异，降雨产流3个随机性，
流域形状和雨量站点分布状况，以及现有预报模型的使用效果与经验确定洪水预报方案，即在清原下水库流域的临近子流域中运用新安江模型进行建模，以SCE-UA 算法对新安江模型参数进行优化率定，最终结合空间临近法和流域物理特性相似法将率定后的参数移植到清原下水库所在子流域上，分别计算产流和汇流，最后得到出流过程。

1.1 研究区概况
清原下水库位于辽宁省浑河上游右岸支流树基沟河内，坝址流域面积65.3km2。

水库坝址附近无气象观测站，东侧20km处有清原气象站，1957年设站观测至今，资料系列较长，根清原气象站1981～2013年降水观测资料分析，多年平均降水
量792.9mm。

水库坝址附近有树基沟雨量站，1964～2012年多年平均降水量862.6mm。

清源下水库流域无水文站和雨量站，临近流域有南口前水文站，南口前水文站位于浑河干流左侧的南口前支流，为径流实验站，南口前站流域面积186km2，与清
源下水库所在树基沟河流域面积163.62km2最接近，为本阶段径流、洪水计算的主要依据站。

流域内水系如图1。

图1 研究区流域水系
1.2 数据来源及处理
1.2.1 数据来源
通过当地水文局收集了以下资料：
（1）南口前流域内3处雨量站1968～1988年降水量摘录资料和日降水量资料。

（2）流域内南口前水文站1968～1988年的洪水摘录资料和日均流量资料。

（3）本次未收集到南口前流域内蒸发资料。

（4）南口前流域内历次暴雨洪水调查资料若干。

（5）对南口前流域各个单元的流域特征值包括河长、河流纵比降、流域面积等在大比尺地形图上进行了重新量算。

（6）用泰森多边形法对南口前流域内各雨量站面积权重进行了量算。

1.2.2 子流域划分
根据南口前站地理位置，以南口前站作为南口前流域的预报断面，清原抽水蓄能电站下水库入库断面为研究子流域的预报断面。

在ARCGIS中的水文分析工具对流域DEM底图进行处理。

通过填洼——流向分析——汇流——河网提取——河网连接——河网分级——捕捉倾泻点——划分子流域流程得到南口前流域文件和研究子流域文件。

南口前流域划分、研究子流域划分和两个子流域的位置关系如图2。

图2 研究区流域与南口前流域位置
1.2.3 基于GIS处理的泰森多边形
鉴于南口前流域中有阿尔当站、海阳站和暖泉南沟站3 个雨量站，故创建雨量站SHP 图层，在ARCGIS中通过分析工具——领域分析——创建泰森多边形流程打开创建泰森多边形工具，设置参数，输入要素选择构造的雨量站，并设置环境后泰森多边形图层，如图3。

图3 裁剪后泰森多边形图层
2 参数率定与结果分析
2.1 参数率定
本次模型率定，充分利用流域内暴雨洪水调查资料和南口前站实测洪水资料，为各单元的参数率定提供了保证，使率定的参数更加合理。

率定时，首先选用近期雨量和流量资料控制较好、观测精度较高洪水进行参数率定，为早期洪水参数率定积累经验、奠定基础。

最终本次选用11场洪水进行预报模型参数率定，用SCE-UA算法进行新安江模型参数优化率定，率定结果如表1。

表1 南口前流域新安江模型率定结果时期编号实测值预报值误差/洪峰流量/
（m3/s）峰现时间确定性系数% 预报值实测值误差/h 1# 045.42 026.22 42.28 1970-08-03 12：00：00 1970-08-03 09：00：00 3 0.76 2# 130.00 111.26 14.41 1971-08-03 12：00：00 1971-08-01 10：00：00 50 0.51 3# 032.00 022.47 29.78 1973-07-17 09：00：00 1973-07-17 7：00：00 2 0.66率定期
4# 026.90 027.56 02.44 1973-08-02 00：00：00 1973-08-01 18：00：00 6 0.52 5# 307.00 163.34 46.80 1975-08-22 03：00：00 1975-08-21 21：00：00 6 0.31 6# 031.10 024.66 20.72 1975-07-09 23：00：00 1975-07-10 03：00：00 4 0.84 7# 192.00 145.01 24.47 1975-08-01 0：00：00 1975-08-10 0：00：00 0 0.77 8# 038.00 051.85 36.44 1977-08-03 19：00：00 1977-08-03 16：00：00 3 0.68验证期1# 028.70 030.12 04.94 1975-07-23 00：00：00 1975-07-23 00：00：00 0 0.67 2# 130.00 082.18 36.79 1979-07-17 02：00：00 1979-07-17 01：00：00 1 0.78 3# 192.00 135.05 29.66 1975-07-31 23：00：00 1975-07-31 23：00：00 0 0.85
从表1看出率定期，1#，5#洪水的洪峰预报值偏差较大，均达到35%以上，2#
洪水的峰现时间异常，误差高达50h。

验证期内1#洪水的洪峰拟合效果相对较好，整体的峰现时间都拟合的较为准确。

模型最终确定模型系数如表2。

表2 新安江模型洪水参数参数物理意义率定值SM 自由水蓄水容量 37.60 WM
张力水容量 177.32 WUM 上层张力水容量 18.59 WLM 下层张力水容量 68.04 K 蒸发能力折算系数 0.38 B张力水蓄水容量曲线方次 0.29 C深层蒸散发系数 0.11 IMP 流域不透水面积比例 0.02 KG 自由水蓄水库地下水日出流系数 0.25 KSS 自
由水蓄水库壤中流小时出流系数 0.32 EX 自由水蓄水容量曲线指数 0.16 KKSS 壤
中流小时消退系数 1.16 KKG 地下径流小时消退系数 0.90 CS 河网蓄水量退水系
数 0.96 l河网汇流滞时 2.93
2.2 成果分析
选用研究区4场实测洪水来验证利用参数移植法建立的新安江模型，验证结果如
表3，洪水过程线如图4。

表3 洪水预报验证结果洪峰流量/（m3/s）峰现时间洪号实测值预报值误差/%
实测值预报值误差/h确定性系数1#0722 089.70 082.35 08.20 07-22 04：00：00 07-22 03：00：00 1 0.92 2#0728 061.00 070.09 14.90 07-28 17：00：00 07-28 17：00：00 0 0.89 3#0805 136.00 105.97 22.08 08-05 13：00：00
08-05 12：00：00 1 0.80 4#0814 079.50 065.41 17.72 08-14 20：00：00
08-15 03：00：00 7 0.93
图4 验证期洪水过程线
从表3可看出，验证期所选场次洪水，洪峰流量合格率75%，峰现时间合格率100%，确定性系数为乙等，验证期洪水场次的模拟效果相对较好，从图4可看出，模型对于单峰型洪水拟合效果相对较好。

尽管本次预报方案已有初步成果，但限于资料条件，还不能完全满足规范要求，方案存在如下问题：①根据规范规定，一个方案编制分析的洪水场次不低于20次，由于资料有限，新安江模型参数率定，勉强选用了11场洪水，不满足规范要求。

②方案编制过程中，由于雨洪配套资料太少（有些年份雨量站和水文站数据不完整，有缺失），给资料分析带来了很大困难，直接影响了参数率定。

3 结语
（1）利用参数移植法，构建了浑河上游南口前流域新安江模型，并将模型参数移植到清原下水库流域，最终结果验证该方法可行，由于资料有限，本次预报精度不高，还需后续收集更多的暴雨洪水资料，进一步改进模型参数。

（2）无资料地区洪水预报是水文预报工作中的难点，后续还需投入更多精力，结合新方法去深入探究。

【相关文献】
［1］张漫莉.基于参数分解法的无资料流域洪水预报研究［J］.水电能源科学，2013（11）.
［2］宋昭义.新安江模型在无资料流域冷水河的应用研究［J］.环境科学与技术，2016（S2）. ［3］冯娇娇.辽宁东部短缺资料地区洪水预报方法研究［D］.大连：大连理工大学，2018.
［4］吴碧琼，黎小东，敖天其，等.基于参数移植法的SWAT模型在缺资料地区非点源污染模拟中的应用［J］.工程科学与技术，2015（S2）：9-16.
［5］程艳，敖天其，黎小东，等.基于参数移植法的SWAT模型模拟嘉陵江无资料地区径流［J］.农业工程学报，2016，32（13）：81-86.
［6］徐长江，杨无双，汪青静，等.无资料地区水文模型参数估算方法比较研究［J］.中国农村水
利水电，2018，427（5）：109-113.
［7］张露，张佳宾，梁国华，等.基于API模型与新安江模型的察尔森水库洪水预报［J］.南水北
调与水利科技，2015，13（6）：1056-1059.
［8］陈建，王建平，谢小燕，等.考虑人类活动影响的改进新安江模型水文预报［J］.水电能源科学，2014（10）：22-25.
［9］刘佩瑶，郝振纯，王国庆，等.新安江模型和改进BP神经网络模型在闽江水文预报中的应用［J］.水资源与水工程学报，2017，28（1）：40-44.
［10］郭俊，周建中，邹强，等.新安江模型参数多目标优化研究［J］.水文，2013，33（1）：1-7.
［11］Huijuna X U，Chen Y，Zeng B，et al.Application of SCE-UA Algorithm to Parameter Optimization of Liuxihe Model［J］.Tropical Geography，2012，32（1）：32-
37.
［12］Dong J，Zhijia L I，Dai J.Application of SCE-UA algorithm to optimization of
Xin'anjiang model parameters［J］.Journal of Hohai University，2012，40（5）：485-490. ［13］Song X M，Zhan C S，Xia J.Integration of a statistical emulator approach with the SCE-UA method for parameter optimization of a hydrological model ［J］.Chinese Science Bulletin，2012，57（26）：3397-3403.
［14］Kang T，Lee S.Modification of the SCE-UA to Include Constraints by Embedding an Adaptive Penalty Function and Application： Application Approach ［J］.Water Resources Management，2014，28（8）：2145-2159.
［15］Gupta V K.Optimal Use of the SCE -UA Global Optimization Method for Calibrating Watershed Models［J］.Journal of Hydrology，1994，158（3-4）：265-284.。