Page 24 - 《中国电力》2026年第3期
P. 24

2026  年 第 59 卷



                  日内调度的其他约束,如功率平衡约束、新                                                T ∑
                                                                                
                                                                                
                                                                                
                                                                                 min   F t (x DN )
                                                                                
              能源出力约束等应与日前调度保持一致。需要注                                             
                                                                                
                                                                                 x DN
                                                                                
                                                                                    t=1                (43)
              意的是,日内调度时对动态碳排放因子进行了修                                                
                                                                                
                                                                                
                                                                                    g(x DN )≤0
                                                                                
                                                                                   
                                                                                
                                                                                
              正,需考虑其修正后新增的约束式(6)。                                                  
                                                                                s.t. 
                                                                                
                                                                                     h(x DN ) = 0
              2.3.2    下层模型
                                                                式中:    x DN 为配电网内部的运行变量;            F t (·)为配电
                  1)目标函数。VPP        的日内调度以最小化总运
                                                                网运行成本函数;          g(·)为配电网中的不等式约束;
              行成本为目标,目标函数为
                                                                h(·)为配电网中的等式约束。
                           ∑ K
                           T 1 ∑
                 C VPP  =min    (C buy,in  +C gt,in  +C PVr,in +    在引入响应变量        r DN 后,式(43)可被改写为
                   in             k,t    k,t   k,t
                                                                             
                           t=1 k=1
                                                                             
                                                                                 T ∑
                                 ′
                                                                              min
                                                                             
                        Wr,in  pcd ,in  cVPP,in                                    F t (x ∗  )+π(r DN )
                       C    +C      +C     )          (42)                             DN
                                                                             
                        k,t    k,t    k,t                                     x DN
                                                                             
                                                                                 t=1                   (44)
                                                                             
              式中:    C t buy,in 为日内调度时刻   t 第  k 个  VPP  与配电                      ∗ DN )≤0
                                                                                g(x
                                                                             
                                                                             
                                                                             
                                                                             
                                                                               
                                                                             s.t.   ∗
                                                                             
              网 交 互 的 购 电 成 本 ;     C gt,in 为 日 内 调 度 时 刻  t 第                  h(x DN ) = 0
                                     t
              k 个  VPP  中燃气轮机的运行成本;           C PVr,in 为日内调
                                                k,t                            r DN = e VPP −e DN = 0   (45)
              度 时 刻  t 第  k 个  VPP  的 弃 光 成 本 ;  C cVPP,in 为 日 内
                                                k,t             式中:   π(r DN )为松弛函数,且     x ∗ DN  = [x DN ,r DN ]; e VPP 、
                                                      ′
              调度时刻     t 第  k 个  VPP  的碳排放成本;     C pcd ,in 为日
                                                   t            e DN 分 别 为 在  VPP  和 配 电 网 中 引 入 的 响 应 变 量 ;
              内调度时刻       t 由于配电网潮流约束导致的第                k 个
                                                                本文中    e DN 为配电网因     VPP  调度导致的有功强制
              VPP  有功强制改变成本;上述各成本详细的计算                                  pcd
                                                                该变量    P j,t  ;  e VPP 为  VPP  中由于配电网潮流约束
              方法与日前调度时相同。                                                              pcd ′
                                                                导致的有功强制改变量            P   。
                  2)约束条件。VPP        日内燃气轮机的出力爬坡                                        j,t
                                                                    将式(44)看作一次性优化问题为
              约束和日内调度时储能设备的偏差约束见文献                      [24],
                                                                                  
                                                                                   min f(x)
              其余约束条件应当与日前调度时保持一致,同时                                               
                                                                                  
                                                                                  
                                                                                  
                                                                                                      (46)
                                                                                     g(x)≤0
              需考虑动态碳排放因子修正后新增的约束。                                                   
                                                                                  
                                                                                  

                                                                                  s.t. 
                                                                                    
                                                                                     h(x) = 0
              3    模型求解                                             根据文献      [30],引入罚函数公式后,响应变
                                                                量 r DN = e VPP −e DN = 0可被转换为不等式形式,即
                  由于上述虚拟电厂-配电网协同调度的日前和                                                
                                                                                      e VPP −e DN ≤0
                                                                                                       (47)
              日内模型均为双层,且            VPP  与配电网两者的模型                     r DN = e VPP −e DN ≡  e DN −e VPP ≤0
                                                                                      
              间存在一定的耦合关系,无法对以上双层模型直                                 式(44)可通过乘数法将与响应变量                  r DN 相关
              接求解。
                                                                的不等式约束从模型的约束条件中去除,即
                  目标级联分析可以将双层模型中上层的目标                                      
                                                                               T ∑
                                                                           
                                                                           
              传播到下层,并通过下层的优化结果反馈调整上                                         min   F t (x ∗  )+µ DN ψ(r DN )
                                                                           
                                                                           
                                                                                     DN
                                                                            x DN
                                                                               t=1
                                                                           
              层决策,实现双向协同优化。这种分层协调方式                                                                     (48)
                                                                           
                                                                              
                                                                           
                                                                              g(x DN )≤0
                                                                           
              天然契合     VPP  与配电网的“指令-响应”交互需                                   
                                                                           
                                                                           s.t. 
                                                                           
                                                                               h(x ∗ DN ) = 0
              求 : 一 方 面 保 留    VPP  对 分 布 式 资 源 的 自 主 调 度
              权,适应新能源出力不确定性等本地约束;另一                             式中:    µ DN 为罚参数;    ψ(·)为罚函数,且必须二次
              方面通过迭代修正确保配电网的全局安全与经济                             可微。需要注意的是,式(44)中                 g(x ∗ DN )≤0变更
              最优性。因此,本文使用目标级联分析法对所提                             为 g(x DN )≤0,是因为通过乘数法已经将响应变量
              双层模型进行转换。                                         r DN 相关的不等式约束已从模型的约束条件中去除。

              3.1    基于目标级联分析法的模型转换                                 VPP  下 层 模 型 的 转 换 方 法 与 上 层 模 型 相 同 。
                  首先以日前调度为例进行模型转换,上层模                           因此,基于目标级联分析法的日前调度双层模型
              型中配电网的运行模型可以表示为                                   被分解为
               20
   19   20   21   22   23   24   25   26   27   28   29