洪水位趋势模型（资料同化应用） · Flood Level Trend Model — Data Assimilation

01

模型定义

基于水位涨率、上游来水与降雨输入的水位短期趋势外推，用于预警窗口判断。

02

公式

常用：状态空间或 ARIMA + 上游传播时间滞后项 H(t+Δ) = f(H(t), dH/dt, Q_up(t−τ), P(t−τ))

公式中各量的定义与单位见"输入 / 输出"。

03

输入与输出

项	内容
输入	资料同化场景：当前水位与涨率、上游水位/流量（滞后 τ）、降雨与观测（站点/雷达/卫星）融合更新初始场，常用 Kalman 滤波或变分同化（3D-Var/4D-Var）。
输出	未来 Δ 时段的水位区间估计

04

参考依据

WMO Guide to Hydrological Practices；相关洪水预报文献

05

成熟度与边界

成熟度

Supported

领域

水文 · Hydrology

适用边界SIYB 立场：趋势外推只用于"预警窗口"，不用于发布确定性水位预报。资料同化需刻画观测误差与背景误差协方差，同化窗口内观测质量控制与异常剔除决定分析场质量。

06

机理与过程链

水文模型沿"产流—汇流"两段展开：产流段决定降雨中有多少变成净雨，常用蓄满产流（土壤缺水量满足后产流）或超渗产流（雨强超过下渗能力即产流），SCS-CN 是工程上最常用的简化形式；汇流段把净雨按坡面与河道演进为出口断面流量，线性水库、单位线与马斯京根法是常用工具。

模型精度取决于两件事：参数与下垫面的匹配（CN 值、糙率、坡度、植被），以及降雨输入的空间代表性。设计复合用参数化模型，实时预警用滚动同化的水文模型，两者对数据频率要求不同：前者年尺度，后者 5~15 min。

07

核心假设与成立条件

假设不成立时，模型偏差通常最先体现在趋势而非绝对值上。

假设	工程含义与成立条件
下垫面稳定	土地利用、土壤类型与河道形态在研究期内不变；城市化或大规模整治后参数须重标
降雨代表性	雨量站网密度满足面雨量误差 <10%（平原约 1 站/50~100 km²），暴雨中心未被漏测
参数时不变	CN、糙率按季节或前期湿度分级，单一参数无法覆盖全年
线性汇流近似	中等洪水下汇流关系近似线性；特大洪水与溃决情景需非线性或二维演进模型

08

误差来源与失效模式

出现下列表现时应先复核假设与输入，而不是调整输出。

失效模式	表现与判据
强降雨空间不均	单站暴雨被面平均稀释，洪峰低估 20%~50%，表现为"雨大洪小"
融雪与雨雪混合	纯降雨模型在春汛期系统性低估，需耦合能量平衡融雪模块
人类调度干扰	水库调度、闸门启闭使实测流量与天然汇流不符，须用调度记录剥离
产流机制转换	同一流域不同场次产流机制不同，单一 CN 值误差可达 30%

09

输入数据条件与精度要求

本模型最小可用数据集：资料同化场景：当前水位与涨率、上游水位/流量（滞后 τ）、降雨与观测（站点/雷达/卫星）融合更新初始场，常用 Kalman 滤波或变分同化（3D-Var/4D-Var）。

数据项	采样与精度要求
降雨	5 min 累计，精度 0.2 mm，站网密度满足面雨量误差 <10%
水位/流量	15 min，水位精度 ±1 cm，流量按率定曲线或走航 ADCP
蒸发/土壤水	日尺度，用于产流参数与前期湿度分级
下垫面参数	CN/糙率按土地利用与土壤类型分区，3~5 年重标一次

输入

资料同化场景：当前水位与涨率、上游水位/流量（滞后 τ）、降雨与观测（站点/雷达/卫星）融合更新初始场，常用 Kalman 滤波或变分同化（3D-Var/4D-Var）。

输出

未来 Δ 时段的水位区间估计

10

部署、标定与运行要点

围绕洪水位趋势模型（资料同化应用）的运行前检查清单。

环节	要点
采样频率	降雨 5 min 采样、面雨量 15 min 更新；水位与流量 15 min 同步
率定策略	3~5 年连续序列，分率定期（70%）与验证期（30%）
在线同化	以实测水位滚动修正状态变量，误差超阈值触发模型重算
预警联动	预报流量超警戒 80% 时提升会商频次，同步核对上游雨量与闸控状态

11

验证指标与参考阈值

阈值为工程常用参考区间，具体项目按合同与规范取值。

指标	参考阈值
NSE	>0.75 良好、0.65~0.75 可用、<0.5 不可用于预警
KGE	>0.7
洪峰相对误差	<10%（大洪水）、<20%（中小洪水）
峰现时间误差	<3 h（中等流域）；水量平衡误差 <5%

12

常见问答

为什么洪峰总是偏低？
多为面雨量代表性不足或 CN 偏小；先核雨量站密度与暴雨中心，再按前期湿度分级重标 CN。

NSE 高就说明模型可用吗？
NSE 对洪峰与低流量敏感度不同，需同时看 KGE、洪峰误差与峰现时间；枯季 NSE 常虚高。

实时预警能直接用设计模型吗？
不能。设计模型参数偏保守，预警需滚动同化与实时雨量驱动，并保留人工会商环节。

14

与相邻方法的取舍

同一目标量常有多种建模路线，取舍取决于数据条件与精度要求。

方法	适用	局限
概念性模型（SCS-CN/线性水库）	参数少、可解释、工程通用	极端事件与机制转换时误差大
分布式物理模型	空间过程完整、可做情景分析	数据与算力要求高，参数不确定性大
数据驱动/机器学习	拟合历史洪水能力强	外推与极值能力弱，需物理约束