Page 55 - 2023年第54卷第4期
P. 55
图 3 基于反向拟合算法的实时校正技术示意
RMSE(k)和 MAE(k)按下式计算:
m
1 )
t - k
RMSE(k) = ∑ (Z- Z ) 2 (8)
槡
t
t
m t =1
1 m )
MAE(k) = ∑ Z- Z t - k (9)
t
t
m t =1
3.5 研究方案设置 表 2列出了本文考虑的 14种方案。其中,“IM” 代表预报调度集成法,作为水
位预报的传统方法,“ TM” 则代表考虑时变参数的预报调度集成法。“IM” 或 “TM” 后一位数字,
“0” 代表未将不同预见期的预报水位纳入目标函数;“1” 代表将不同预见期的预报水位纳入目标函
数。小数点后一位,“1” 和 “2” 分别表示采用最新已知的误差和采用迭代估计的误差进行基于反向
拟合算法的实时校正,“3” 表示采用最新已知误差进行基于联合率定的实时校正。
表 2 水库水位实时预报方案设置
实时校正技术
考虑预见期的
反向拟合 反向拟合 联合率定
预报技术
(采用最新已知的误差) (采用迭代估计的误差) (采用最新已知的误差)
IM
IM0.1 √
IM0.2 √
非时变参数方案 IM1 √
IM1.1 √ √
IM1.2 √ √
IM1.3 √ √
TM
TM0.1 √
TM0.2 √
时变参数方案 TM1 √
TM1.1 √ √
TM1.2 √ √
TM1.3 √ √
4 结果与讨论
4.1 时变参数函数式的确定 构建预报调度集成模型,输入水布垭水库流域的气象驱动数据,采用
Morris方法进行敏感性分析,得出新安江模型的敏感参数为蒸散发折算系数(KE)、表征流域蓄水容量
不均匀性的指数( B)和壤中流消退系数(CI),在时变参数估计过程中非敏感参数保持不变。采用 SSC -
DP以月尺度识别水文模型敏感参数 KE、B和 CI的变化过程,设置了 SSC - DP不同的参数方案,水位
模拟结果以 NSE评价,并汇总于表 3。可知,采用时变参数方法能显著降低水位模拟的误差。
— 4 3 1 —