Page 97 - 2023年第54卷第8期
P. 97

站、光伏电站的装机成本,元?kW;C                   、C 分别为风能电站、光伏电站的运行成本,元?(MWh);
                                                 om,w  on,s
                                                                                            为风电、光伏多
               w   s                                w    s
              I、I为风电和光伏的装机容量,MW;P 、P 为风电和光伏发电量,kWh;β w                                  、β s
              年平均出力系数;d为折现率;T为一年的运行小时数。
                                r
                  在各经济技术参数恒定的情况下,联立式( 1)、(6)—(7)后得,互补系统净效益为一多元函数如
              式( 8),自变量为风电装机 I和光伏装机 I。为求解最大净效益值所对应的装机容量组合即最优装机
                                        w
                                                     s
                    
                         
              容量( I,I),对两自变量分别求偏导可求解最优风光装机容量函数式如式(9)。
                         s
                    w
                                 2           2
                                        BTI + [(B - cB - C
                                                                                           )a T - C ]I
                                                                                 s
                                                                                   2
                                                                      in,w
                               w
                         1
                                                                           w
                    ws
                  NR =- ka β w B TI- ka β s s  s  w  1  w  om,w )a β w T - C ]I+ [(B- cB- C om,s β s  in,s  s
                                                                                     s
                                 w
                                     2
                                                                                                        (8)
                                                      (B - cB - C   )T - C
                                                   a β w  w  1  w  om,w  in,w
                                                
                                               I =
                                                w
                                                            1
                                                                  w
                                                          2ka β w B T
                                                                                                        (9)
                                                     (B- cB- C     )T - C
                                                   a β s  s  2  s  om,s  in,s
                                                
                                               I =
                                                s
                                                               BT
                                                            2
                                                          2ka β s s
                        Y     1
                                                      
                      ∑
              式中:a =                为多年折现系数;I、I 分别为最优风电装机和最优光伏装机,MW。
                                                          
                                                      w
                                                          s
                       y =1(1+d) y - 1
                                r
                  可以看出,最优装机函数形式十分简洁,在各规划参数已知的情况下,只需获取该地区风光多年
              平均出力系数和求解弃电率函数,即可计算与原有水电站相配套的风光装机容量。如计算结果为负
              值,表明从经济角度出发该地区不适合风光互补电站的建设。
              2.3 经济技术参数敏感性 在装机规划中,经济技术参数主要包括风光上网电价 B 、B,风光装机
                                                                                                s
                                                                                           w
              成本 C 、C ,风光运行维护成本 C                om,w 、C 和多年折现系数 a。在实际情况下,各参数时刻受到
                                                        on,s
                    in,w
                          in,s
              政策、经济、技术等因素影响而并非恒定不变,同时式(9)表明各参数的变化会改变系统的最优装机
              容量,同时对互补系统的经济效益产生显著影响。因此,有必要进行敏感性分析获得对最优装机和系
              统净效益敏感和不敏感的参数,以应对未来时期内各参数变动带来的影响。
                  为准确量化各参数在原有情况变动下对最优装机容量的影响,对式( 9)在最优装机容量值处泰勒
              展开得式( 10)—(13),代入已知其他参数可直接解析相关变量的敏感性。联立式(8)(9)可得到系统
              净效益与经济技术参数的关系如式( 14)。
                                             C   T + C         C    T + C
                                          a β w  om,w  in,w  a β w  om,w  in,w
                                       
                                                                             2
                                    Δ I =          2   Δ B -          3   Δ B + o( Δ B )
                                                                             w
                                                                                    w
                                       w
                                                          w
                                            2ka β w B T       2ka β w B T
                                              1    w            1     w
                                                                                                       (10)
                                            C   T + C         C   T + C
                                          a β s om,s  in,s  a β s om,s  in,s
                                       
                                                                          2
                                    Δ I =         2   Δ B-          3   Δ B+ o( Δ B)
                                                         s
                                                                          s
                                                                                 s
                                       s
                                                                  BT
                                                 BT
                                                               2
                                              2
                                           2ka β s s         2ka β s s
                                               1
                                       
                                    Δ I =-          Δ C in,w
                                       w
                                             1    w
                                           2 ka β w B T
                                                                                                       (11)
                                               1
                                       
                                    Δ I =-          Δ C in,s
                                       s
                                                BT
                                           2 ka β s s
                                             2
                                             1
                                       
                                    Δ I =-      Δ C om,w
                                       w
                                           2kB
                                             1 w
                                                                                                       (12)
                                             1
                                       
                                    Δ I =-      Δ C om,s
                                       s
                                           2 kB
                                             2 s
                                             C
                                       
                                    Δ I =     in,w  Δ a
                                       w
                                                  2
                                               B aT
                                            1   w
                                          2kβ w
                                                                                                       (13)
                                             C in,s
                                       
                                    Δ I =        2  Δ a
                                       s
                                              BaT
                                            2
                                          2kβ s s
                                          (B - cB - C            2    (B- cB- C     )T - C ] 2
                                                             in,w
                                     [a β w  w  1 w  om,w )T - C ] [a β s  s  2  s  om,s  in,s
                            maxNR =                              +                                     (14)
                                  ws
                                                                                 BT
                                               4ka β w B T                 4ka β s s
                                                 1    w                      2
                                                                                                —  9 8 1 —
   92   93   94   95   96   97   98   99   100   101   102