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 —
   50   51   52   53   54   55   56   57   58   59   60