Page 82 - 2025年第56卷第1期

按照河道上下游拓扑关系将 ３４ 个计算单元进一步分成滹沱河与滏阳河 ２ 个计算单元组实施并行
计算（图 ３）。

图 ３ 子牙河模型计算单元拓扑关系图

３．３ 模型校准 水文模型结果采用洪峰相对误差（Ｒｅｌａｔｉｖｅ Ｅｒｒｏｒ ｏｆ Ｆｌｏｏｄ Ｐｅａｋ， ＲＥ ）、 纳什系数（Ｎａｓｈ－
Ｐ
Ｓｕｔｃｌｉｆｆｅ Ｃｏｅｆｆｉｃｉｅｎｔ ｏｆ Ｅｆｆｉｃｉｅｎｃｙ， ＮＳＥ）进行精度评估。
（１）洪峰相对误差（ＲＥ ）：
Ｐ
Ｑ －Ｑ
ＲＥ ＝ ｐ，ｓｉｍ ｐ，ｏｂｓ ×１００％ （７）
Ｐ
Ｑ
ｐ，ｏｂｓ
式中： Ｑ 为模拟洪峰， ｍ ∕ｓ； Ｑ 为实测洪峰， ｍ ∕ｓ。
３
３
ｐ，ｓｉｍ ｐ，ｏｂｓ
（２）纳什系数（ＮＳＥ）：
Ｔ ｔ ｔ
∑ （Ｑ －Ｑ ） ２
ｔ ＝ １ ｏｂｓ ｓｉｍ
ＮＳＥ ＝ １－ （８）
Ｔ ｔ
∑ （Ｑ － ２
Ｑ ）
ｔ ＝ １ ｏｂｓ ｏｂｓ
ｔ ３ ｔ ３ ３
式中： Ｑ 为 ｔ 时刻的观测流量， ｍ ∕ｓ； Ｑ 为 ｔ 时刻的模拟流量， ｍ ∕ｓ； Ｑ 为观测流量的平均值， ｍ ∕ｓ。
ｏｂｓ ｓｉｍ ｏｂｓ
水动力模型基于遥感监测获得的淹没范围， 通过调整河道、 淹没区地表糙率及部分位置微地形对
模型进行校准。
４ 结果分析
４．１ 产流机制分析 选择 ２０２３ 年 ７ 月 １５ 日 ８∶ ００ 至 ２０２３ 年 ８ 月 ８ 日 ８∶ ００ 降雨数据， 进行子牙河逐小

时分布式水文模拟， 计算总耗时 ２ ｍｉｎ。 模型重要参数取值范围见表 １。 本次模拟洪峰平均误差 ７％，
纳什系数平均值 ０．７。 提取子牙河上、 中、 下游 １０ 个重要站点模拟结果进行进一步分析（表 ２， 图 ４）。
表 １ 子牙河 “２３·７” 大洪水复盘模拟模型重要参数取值范围
参数符号 参数名称 单位 参数取值范围
土壤∕水力传导系数 －７ －６
ｋ ｓ ｍ∕ｓ １×１０ ～１×１０
优先流面积比 ０．２～０．８
ｐｅｒ ｐ
壤中流出流线性系数 ０．２～０．８
ｋ ｉ１
优先流出流线性系数
ｋ ｐ１ ０．２～０．８
壤中流出流非线性系数
ｋ ｉ２ ０．２～０．８
优先流出流非线性系数
ｋ ｐ２ ０．２～０．６
— ７ ７ —