Page 63 - 2023年第54卷第11期
P. 63

式中 DOY 为 t时刻对应的日期在年内的日序数(dayofyear,DOY),不考虑闰年的情况,DOY 取值为
                       t                                                                            t
              0~365d。
                  水库作为调节径流时空分布不均的工程措施,总是存在蓄水期和放水期。考虑到 cos函数与 sin函
              数具有严格的周期性且仅有一个波峰和一个波谷,可以与水库每年的蓄水期、放水期相对应,因此,
              本文提出了一个以 cos函数为基础的通用 TCC函数:
                                                   1 + λ 1 - λ              2 π
                                                                [
                                         f(DOY) =     +    ·cos(DOY- k)·       ]                        (2)
                                               t                      t
                                                    2    2                  365
              式中:λ为 TCC函数的最小值,范围为 0~1;k为水库目标库容达到最大值的日序数,范围为 0~365。
              如图 1(a)所示,TCC函数的特点是当 DOY= k时,f(·)达到最大值 1;若 k<183,则当 DOY= k + 183
                                                                                                   t
                                                     t
              时,TCC函数达到最小值 λ ,若 k ≥183,则当 DOY= k - 182 时,TCC函数达到最小值 λ 。
                                                             t
              2.1.2 水库出流计算 TCCRM模型的基本结构如图 1(b)所示,将水库出流量 Q 分为得到发挥经济
                                                                                        out
                                                               spill
              效益的有效出流量 Q         yield 以及未发挥经济效益的弃水 Q 两部分。
                                                               out
                                 out
                                                     Q   = Q yield + Q spill                            (3)
                                                      out,t  out ,t  out ,t
                                                                                                      [21]
                  Q yield 主要受水库功能的影响,是水库为了满足灌溉、发电等功能而下泄的流量。参考 Doll 在
                    out ,t
              2003年提出的 ROSM(ReservoirsimulationoutflowModelproposedbyDoll,DRM),基于水库水的滞留时
              间( LagTime,LT)和一个根据水库实际蓄水量计算得到的水库库容校正系数进行 Q                                 yield 的计算。
                                                                                          out,t
                                                        V- V
                                                       {  t  d ·ζ t  V≥V d
                                                                    t
                                                  yield
                                                 Q   = LT  t                                            (4)
                                                  out ,t
                                                           0       V<V
                                                                    t  d
                                                                                      为水库在 t时刻的库容
                                                      t
                       t
              式中:LT 为水库水在 t时刻的滞留时间;V 为水库在 t时刻的实际蓄水量;ζ t
              校正系数。
                  假设 V  targ,t 从两方面对 Q  yield 产生影响:(1)V   targ,t 较高时水库处于倾向于蓄水,Q           yield 较低,反之水库
                                                                                          out ,t
                                        out ,t
              倾向于放水,Q       yield 较高;(2)V 应尽可能向 V      targ,t 靠近,当 V<V targ,t 时,减小 Q yield 以蓄水,反之则增大
                                          t
                            out ,t
                                                                                     out ,t
              Q yield 以放水。
                out ,t
                                                                                   保持不变的条件下,LT 越
                  首先考虑第一点,LT 表征了水库的放水速度,根据式(4),在 V 与 ζ t                                                t
                                                                             t
                                     t
              大,Q   yield 越小,反之则 Q  yield 越大。因此,在 DRM的基础上做出改进,假设 LT 在年内并非保持不变,而
                    out,t           out,t                                           t
              是随着 V      的变化而变化:当水库的 V              较大时,水库在该季节倾向于蓄水,水库水体更新速度变慢,
                      targ,t                     targ,t
              LT 增加;而当水库的 V        targ,t 较小时,水库在该季节倾向于放水,水库水体更新速度加快,LT 下降。如式
                                                                                                t
                t
              ( 1)所示,TCC函数 f(·)指示了 V           的波动变化过程,V           与之成正比。因此,基于 f(·)来计算 LT:
                                            targ,t              targ,t                                 t
                                                  LT= LT+ LT·f(DOY)                                     (5)
                                                    t    0   1       t
              式中 LT、LT 为模型参数,需要率定,参数范围为 0~365d。
                          1
                     0
                  为了使 V 应尽可能向 V            靠近,在 DRM的基础上对 ζ t          进行重构,使其根据 V 与 V              的大小关
                          t             targ,t                                            t    targ,t
              系对 Q  yield 进行调整,当 V<V    targ,t 时,减小 Q yield 以蓄水,当 V>V  targ,t 时则增大 Q yield 以放水。DRM所采用
                                     t
                                                                    t
                                                      out ,t
                    out ,t
                                                                                     out ,t
                                DRM                          DRM        DRM
              库容校正系数记为 ζ t         ,水库的实际蓄水量越大,ζ t               越大,ζ t    的计算如式(6)所示。根据 V            targ,t 重
                                     TCCRM  TCCRM
              构的库容校正系数记为 ζ t             ,ζ t   的计算如式(7)所示。
                                                          (       )                                     (6)
                                                             t
                                                      DRM   V- V d  1.5
                                                         =
                                                     ζ t   V - V d
                                                            max
                                                           (       )                                    (7)
                                                              t
                                                     TCCRM   V- V d  α
                                                         =
                                                    ζ t     V targ,t - V d
              式中:V 为水库的最大库容;α为模型参数,需要率定,参数范围为 0~5。
                      max
                  最后计算弃水 Q       spill ,产生 Q spill 的原因主要有两点,由于水库调度能力有限产生的弃水以及为了维
                                 out,t
                                            out,t
              持下游河道生态环境而下泄的生态流量。Q                    spill 在一定程度上与水库入流 Q 存在正相关关系:在汛期,
                                                                                 in,t
                                                     out ,t
              Q 越大,因水库防洪库容与下泄能力有限而产生的弃水量越大;而在非汛期,虽然水库在此期间产
                in,t
                                                                                              —   1 2 5 —
                                                                                                   3
   58   59   60   61   62   63   64   65   66   67   68