本发明公开了一种基于约束聚合的虚拟电厂出力可行域辨识方法及装置，涉及电力系统调度运行技术领域，其中，该方法包括：获取虚拟电厂内部配电系统运行基础数据；根据虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型；根据虚拟电厂静态安全约束模型构建虚拟电厂出力可行域模型；根据虚拟电厂内部配电系统运行基础数据和可行域顶点枚举算法对虚拟电厂出力可行域模型进行计算得到虚拟电厂出力可行域。该方法可以刻画安全运行允许的虚拟电厂出力可行域，使得虚拟电厂可以参与电力系统优化调度运行和电力现货市场出清，确保调度和出清结果不会违反虚拟电厂内部的安全运行约束，进而可提升系统的安全性和经济性。
1.一种基于约束聚合的虚拟电厂出力可行域辨识方法，其特征在于，包括以下步骤：S1，获取虚拟电厂内部配电系统运行基础数据；
S2，根据所述虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型；
S3，根据所述虚拟电厂静态安全约束模型构建虚拟电厂出力可行域模型；
所述虚拟电厂静态安全约束模型，包括：其中，式(1)和式(2)分别为节点i和节点j之间的配电线路的有功和无功潮流方程约束， 分别为配电线路的有功和无功潮流，bij、gij分别为配电线路的串联电纳、串联电导，vi、vj分别为节点i和节点j电压幅值，θi、θj分别为节点i和节点j电压相角， 是节点编号集合， 表示与节点i通过配电线路直接相连的节点编号的集合；
式(3)和式(4)分别是节点有功和无功平衡方程约束，Pi、 分别是节点i处分布式发电d
的有功出力和无功出力，Pi 、 分别是节点i的有功负荷和无功负荷，P0、Q0分别是虚拟电厂的有功和无功出力；
式(5)表示节点i和节点j之间的配电线路、配电变压器的传输容量约束，式中 分别是节点i和节点j之间的配电线路和配电变压器的视在功率传输上限；
式(6)表示配网节点电压上下限约束，式中 分别为节点i允许的电压幅值上限和下限；
式(7)为分布式发电的有功和无功出力上下限约束，式中 是节点i的分布式发电的有功出力上下限， 是节点i的分布式发电的无功出力的上下限；
S31，根据所述虚拟电厂静态安全约束模型确定关于(v,θ,P ,Q ,P0,Q0)的虚拟电厂静态安全约束可行域记Ω：
Ω＝{(v,θ,P ,Q ,P0,Q0)|s.t.(1)～(7)} (8)g g g
其中，v、θ、P、Q分别是以vi、θi、Pi、 为元素的向量；
S32，消去虚拟电厂内部状态变量v、θ和内部决策变量P 、Q ，将虚拟电厂静态安全约束聚合为关于虚拟电厂有功和无功出力的约束，得到虚拟电厂出力可行域模型：其中，Φ为所述虚拟电厂出力可行域；
S4，根据所述虚拟电厂内部配电系统运行基础数据和可行域顶点枚举算法对所述虚拟电厂出力可行域模型进行计算得到虚拟电厂出力可行域；
S42，根据所述虚拟电厂内部配电系统运行基础数据计算初始顶点集合和初始近似多边形；
所述虚拟电厂出力可行域Φ的初始顶点集包含4个可行域顶点，通过求解如下优化问题计算得到：
其中， 是优化问题目标函数系数，对 依次取值(1,1)、(1,‑1)、(‑1,
1)、(‑1,‑1)，通过非线性优化算法或调用非线性优化求解器对式(10)进行求解，对应的优化问题最优解记为 vm为第m个初始顶点；将初始顶点构成集合V0＝{v1,v2,v3,v4}，将可行域顶点集合初始化为V＝V0；将初始顶点构成的多边形各个边界的方程的集T
S44，通过向外平移F 中的边界搜索所述虚拟电厂出力可行域Φ的边界点，对于F中的第k个边界，边界的方程为 边界两端的顶点分别为 顶点搜索的计算模型如下：
通过非线性优化算法或调用非线性优化求解器求解顶点搜索的计算模型，求解得到顶点最优解记为 顶点最优值记为S45，将顶点最优解加进可行域顶点集合V；
S46，根据设定的所述虚拟电厂出力可行域的误差限值和步骤S44中的顶点搜索的改进量判断是否需要增加待平移边界；
S461，步骤S44中顶点搜索的改进量为：(k) (k)
S462，在Δgk>ε时，分别计算vk与 确定的直线的方程f1 、 将f1 、 加N T
入新增边界集合F，返回步骤S44，使用F中的第k+1个边界再进行顶点搜索；
S463，在Δgk≤ε时，返回步骤S44，使用F中的第k+1个边界进行顶点搜索；
S47，对F中的所有边界方程进行顶点搜索后，判断F是否为空集；
S471，如果F不为空集，将待搜索边界集F 更新为F ，即F←F ，并返回步骤S43对更新后T
S472，如果F为空集，则结束，根据顶点集合V中的顶点确定多边形，根据确定的多边形得到所述虚拟电厂出力可行域。
2.根据权利要求1所述的基于约束聚合的虚拟电厂出力可行域辨识方法，其特征在于，所述虚拟电厂内部配电系统运行基础数据，包括：配网连接拓扑、配电线路电导和电纳参数、变压器的电导和电纳参数、配电线路传输容量极限、安全运行允许的节点电压上下限、配网内部分布式电源有功和无功出力上下限、配网内部负荷节点的有功和无功负荷预测。
3.一种基于约束聚合的虚拟电厂出力可行域辨识装置，其特征在于，包括：获取模块，用于获取虚拟电厂内部配电系统运行基础数据；
第一构建模块，用于根据所述虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型；
第二构建模块，用于根据所述虚拟电厂静态安全约束模型构建虚拟电厂出力可行域模型；所述虚拟电厂静态安全约束模型，包括：其中，式(1)和式(2)分别为节点i和节点j之间的配电线路的有功和无功潮流方程约束， 分别为配电线路的有功和无功潮流，bij、gij分别为配电线路的串联电纳、串联电导，vi、vj分别为节点i和节点j电压幅值，θi、θj分别为节点i和节点j电压相角， 是节点编号集合， 表示与节点i通过配电线路直接相连的节点编号的集合；
式(3)和式(4)分别是节点有功和无功平衡方程约束，Pi、 分别是节点i处分布式发电d
的有功出力和无功出力，Pi 、 分别是节点i的有功负荷和无功负荷，P0、Q0分别是虚拟电厂的有功和无功出力；
式(5)表示节点i和节点j之间的配电线路、配电变压器的传输容量约束，式中 分别是节点i和节点j之间的配电线路和配电变压器的视在功率传输上限；
式(6)表示配网节点电压上下限约束，式中 分别为节点i允许的电压幅值上限和下限；
式(7)为分布式发电的有功和无功出力上下限约束，式中 是节点i的分布式发电的有功出力上下限， 是节点i的分布式发电的无功出力的上下限；
根据所述虚拟电厂静态安全约束模型确定关于(v,θ,P ,Q ,P0,Q0)的虚拟电厂静态安全约束可行域记Ω：
Ω＝{(v,θ,P ,Q ,P0,Q0)|s.t.(1)～(7)} (8)g g g
其中，v、θ、P、Q分别是以vi、θi、Pi、 为元素的向量；
消去虚拟电厂内部状态变量v、θ和内部决策变量P 、Q ，将虚拟电厂静态安全约束聚合为关于虚拟电厂有功和无功出力的约束，得到虚拟电厂出力可行域模型：其中，Φ为所述虚拟电厂出力可行域；
计算模块，用于根据所述虚拟电厂内部配电系统运行基础数据和可行域顶点枚举算法对所述虚拟电厂出力可行域模型进行计算得到虚拟电厂出力可行域；
设定单元，用于设定所述虚拟电厂出力可行域的误差限值ε；
初始顶点与多边形计算单元，用于根据所述虚拟电厂内部配电系统运行基础数据计算初始顶点集合和初始近似多边形；
所述初始顶点与多边形计算单元具体用于，所述虚拟电厂出力可行域Φ的初始顶点集包含4个可行域顶点，通过求解如下优化问题计算得到：其中， 是优化问题目标函数系数，对 依次取值(1,1)、(1,‑1)、(‑1,
1)、(‑1,‑1)，通过非线性优化算法或调用非线性优化求解器对式(10)进行求解，对应的优化问题最优解记为 vm为第m个初始顶点；将初始顶点构成集合V0＝{v1,v2,v3,v4}，将可行域顶点集合初始化为V＝V0；将初始顶点构成的多边形各个边界的方程的集T
新增边界初始化单元，用于初始化新增边界集合F为空集，即T
顶点搜索单元，用于通过向外平移F中的边界搜索所述虚拟电厂出力可行域Φ的边界T
点，对于F 中的第k个边界，边界的方程为 边界两端的顶点分别为 顶点搜索的计算模型如下：通过非线性优化算法或调用非线性优化求解器求解顶点搜索的计算模型，求解得到顶点最优解记为 顶点最优值记为加入单元，用于将顶点最优解加进可行域顶点集合V；
第一判断单元，根据设定的所述虚拟电厂出力可行域的误差限值和所述顶点搜索单元中的顶点搜索的改进量判断是否需要增加待平移边界；
所述第一判断单元，具体用于，顶点搜索的改进量为：(k) (k)
在Δgk>ε时，分别计算vk与 确定的直线的方程f1 、 将f1 、 加入新N T
增边界集合F，返回所述顶点搜索单元，使用F中的第k+1个边界再进行顶点搜索；
在Δgk≤ε时，返回所述顶点搜索单元，使用F中的第k+1个边界进行顶点搜索；
第二判断单元，对F中的所有边界方程进行顶点搜索后，判断F是否为空集；
如果F不为空集，将待搜索边界集F更新为F，即F←F ，并返回所述新增边界初始化单T
如果F为空集，则结束，根据顶点集合V中的顶点确定多边形，根据确定的多边形得到所述虚拟电厂出力可行域。
4.根据权利要求3所述的基于约束聚合的虚拟电厂出力可行域辨识装置，其特征在于，所述虚拟电厂内部配电系统运行基础数据，包括：配网连接拓扑、配电线路的电导和电纳参数、变压器的电导和电纳参数、配电线路传输容量极限、安全运行允许的节点电压上下限、配网内部分布式电源有功和无功出力上下限、配网内部负荷节点的有功和无功负荷预测。
基于约束聚合的虚拟电厂出力可行域辨识方法及装置
技术领域
[0001] 本发明涉及电力系统调度运行技术领域，特别涉及一种基于约束聚合的虚拟电厂出力可行域辨识方法及装置。
背景技术
[0002] 随着分布式风电、光伏、小型燃气轮机、需求响应资源、用户侧储能、电动汽车等各类分布式资源的快速发展和大量接入，传统的输电网与配电网解耦的调度运行模式已经无
法适应，这制约了分布式资源的有效利用。同时，当前的电力批发市场仅允许容量较大的发
电机组参与，配网内的分布式资源难以直接参与电力市场竞争和资源优化配置。此外，配网
内分布式资源运行方式不合理，将可能造成配电网电压越限、潮流阻塞等安全问题。在这一
背景下，如何高效聚合配网内的分布式资源参与输电网优化运行、参与电力系统大范围资
源优化配置、同时保证配网安全可靠运行，成为我国电力系统调度运行以及分布式资源发
展迫切需要解决的现实问题。
[0003] 虚拟电厂是指将含有分布式发电设施的配电系统等值为一个发电厂，其可以根据电网调度机构的调度指令或电力现货市场的出清结果调节配电系统在并网关口处与输电
网交换的有功和无功功率，从而使得配网内的分布式资源可以参与输电网优化运行以及现
货市场资源优化配置。将虚拟电厂视作传统发电厂参与输电系统和电力现货市场优化，需
要知道虚拟电厂可以提供的有功和无功出力可行范围。已有文献研究了实现虚拟电厂实现
分布式发电与配电网协同优化的定价策略和方法，但其没有充分考虑虚拟电厂内的各类安
全运行约束。已有文献提出了考虑网络约束的虚拟电厂参与电力市场优化模型，但这一模
型需要假设其它各发电设备的成本和物理参数已知，这无法适应实际应用需求。已有文献
研究了利用虚拟电厂提供无功支撑和无功辅助服务的机制，但其没有考虑虚拟电厂在不同
的有功出力下允许的无功出力范围变化。
[0004] 现有文献和技术尚无法计算和刻画虚拟电厂出力可行域、无法在虚拟电厂参与输电网和现货市场的优化过程中简洁、有效地考虑虚拟电厂的安全运行约束。
发明内容
[0005] 本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
[0006] 为此，本发明的一个目的在于提出一种基于约束聚合的虚拟电厂出力可行域辨识方法，该方法将虚拟电厂内各类分布式资源的运行约束以及配电网安全运行约束聚合为关
于虚拟电厂有功和无功出力的显式约束条件，可以刻画安全运行允许的虚拟电厂出力可行
域，确保调度和出清结果不会违反虚拟电厂内部的安全运行约束，提升系统的安全性和经
济性。
[0007] 本发明的另一个目的在于提出一种基于约束聚合的虚拟电厂出力可行域辨识方法装置。
[0008] 为达到上述目的，本发明一方面实施例提出了一种基于约束聚合的虚拟电厂出力可行域辨识方法，包括：
[0009] S1，获取虚拟电厂内部配电系统运行基础数据；
[0010] S2，根据所述虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型；
[0011] S3，根据所述虚拟电厂静态安全约束模型构建虚拟电厂出力可行域模型；
[0012] S4，根据所述虚拟电厂内部配电系统运行基础数据和可行域顶点枚举算法对所述虚拟电厂出力可行域模型进行计算得到虚拟电厂出力可行域。
[0013] 本发明实施例的基于约束聚合的虚拟电厂出力可行域辨识方法，通过将含有分布式电源和需求响应资源的配电网等值为一个虚拟电厂参与输电网的优化调度运行；基于静
态安全约束模型，通过约束聚合构建虚拟电厂出力可行域模型；通过一种顶点枚举算法，计
算虚拟电厂出力可行域的边界顶点，从而得到描述虚拟电厂出力可行域的线性不等式约
束。将虚拟电厂内各类分布式资源的运行约束以及配电网安全运行约束聚合为关于虚拟电
厂有功和无功出力的显式约束条件，可以刻画安全运行允许的虚拟电厂出力可行域，使得
虚拟电厂可以像同步发电机一样参与电力系统优化调度运行和电力现货市场出清，确保调
度和出清结果不会违反虚拟电厂内部的安全运行约束，进而可提升系统的安全性和经济
性。
[0014] 另外，根据本发明上述实施例的基于约束聚合的虚拟电厂出力可行域辨识方法还可以具有以下附加的技术特征：
[0015] 进一步地，在本发明的一个实施例中，所述虚拟电厂内部配电系统运行基础数据，包括：
[0016] 配网连接拓扑、配电线路的电导和电纳参数、变压器的电导和电纳参数、配电线路传输容量极限、安全运行允许的节点电压上下限、配网内部分布式电源有功和无功出力上
下限、配网内部负荷节点的有功和无功负荷预测。
[0017] 进一步地，在本发明的一个实施例中，所述虚拟电厂静态安全约束模型，包括：
[0018]
[0019]
[0020]
[0021]
[0022]
[0023]
[0024]
[0025] 其中，式(1)和式(2)分别为节点i和节点j之间的配电线路的有功和无功潮流方程约束， 分别为配电线路的有功和无功潮流，bij、gij分别为配电电路的串联电纳、串
联电导，vi、vj分别为节点i和节点j电压幅值，θi、θj分别为节点i和节点j电压相角， 是节
点编号集合， 表示与节点i通过配电线路直接相连的节点编号的集合；
[0026] 式(3)和式(4)分别是节点有功和无功平衡方程约束，Pig、 分别是节点i处分布d
式发电的有功出力和无功出力，Pi、 分别是节点i的有功负荷和无功负荷，P0、Q0分别是虚
拟电厂的有功和无功出力；
[0027] 式(5)表示节点i和节点j之间的配电线路、配电变压器的传输容量约束，式中分别是节点i和节点j之间的配电线路和配电变压器的视在功率传输上限，N是对传输容
量约束进行分段线性化的分段数；
[0028] 式(6)表示配网节点电压上下限约束，式中 分别为节点i允许的电压幅值上限和下限；
[0029] 式(7)为分布式发电的有功和无功出力上下限约束，式中 是节点i的分布式发电的有功出力上下限， 是节点i的分布式发电的无功出力的上下限。
[0030] 进一步地，在本发明的一个实施例中，所述S3进一步包括：
[0031] S31，根据所述虚拟电厂静态安全约束模型确定关于(v,θ,Pg,Qg,P0,Q0)的虚拟电厂静态安全约束可行域记Ω：
[0032] Ω＝{(v,θ,Pg,Qg,P0,Q0)|s.t.(1)～(7)} (8)g g g
[0033] 其中，v、θ、P、Q分别是以vi、θi、Pi、 为元素的向量；
[0034] S32，消去虚拟电厂内部状态变量v、θ和内部决策变量Pg、Qg，将虚拟电厂静态安全约束聚合为关于虚拟电厂有功和无功出力的约束，得到虚拟电厂出力可行域模型：
[0035]
[0036] 其中，Φ为所述虚拟电厂出力可行域。
[0037] 进一步地，在本发明的一个实施例中，所述S4进一步包括：
[0038] S41，设定所述虚拟电厂出力可行域的误差限值ε；
[0039] S42，根据所述虚拟电厂内部配电系统运行基础数据计算初始顶点集合和初始近似多边形；
[0040] 所述虚拟电厂出力可行域Φ的初始顶点集包含4个可行域顶点，通过求解如下优化问题计算得到：
[0041]
[0042] 其中， 是优化问题目标函数系数，对 依次取值(1,1)、(1,‑1)、 (‑1,1)、(‑1,‑1)，通过非线性优化算法或调用非线性优化求解器对式(10)进行求解，
对应的优化问题最优解记为 vm为第m个初始顶点；将初始顶点构成集合 V0
＝{v1,v2,v3,v4}，将可行域顶点集合初始化为V＝V0；将初始顶点构成的多边形各个边界的
T
方程的集合记为F；对V0中顶点计算平均值，即
N
[0043] S43，初始化新增边界集合F为空集，即
[0044] S44，通过向外平移FT中的边界搜索所述虚拟电厂出力可行域Φ的边界点，对于FT(k)
中的第k个边界，边界的方程为f : 边界两端的顶点分别为
顶点搜索的计算模型如下：
[0045]
[0046] 通过非线性优化算法或调用非线性优化求解器求解顶点搜索的计算模型，求解得到顶点最优解记为 顶点最优值记为
[0047] S45，将顶点最优解加进可行域顶点集合V；
[0048] S46，根据设定的所述虚拟电厂出力可行域的误差限值和步骤S44中的顶点搜索的改进量判断是否需要增加待平移边界；
[0049] S461，步骤S44中顶点搜索的改进量为：
[0050](k) (k) (k)
[0051] S462，在Δgk>ε时，分别计算vk与 确定的直线的方程f1 、f2 ，将f1 、 (k) N T
f2 加入新增边界集合F，返回步骤S44，使用F中的第k+1个边界再进行顶点搜索；
[0052] S463，在Δgk≤ε时，返回步骤S44，使用FT中的第k+1个边界进行顶点搜索；
[0053] S47，对FT中的所有边界方程进行顶点搜索后，判断FN是否为空集；
[0054] S471，如果FN不为空集，将待搜索边界集FT更新为FN，即FT←FN，并返回步骤 S43对T
更新后的F中的边界进行顶点搜索；
[0055] S472，如果FN为空集，则结束，根据顶点集合V中的顶点确定多边形，根据确定的多边形得到所述虚拟电厂出力可行域。
[0056] 为达到上述目的，本发明另一方面实施例提出了一种基于约束聚合的虚拟电厂出力可行域辨识装置，包括：
[0057] 获取模块，用于获取虚拟电厂内部配电系统运行基础数据；
[0058] 第一构建模块，用于根据所述虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型；
[0059] 第二构建模块，用于根据所述虚拟电厂静态安全约束模型构建虚拟电厂出力可行域模型；
[0060] 计算模块，用于根据所述虚拟电厂内部配电系统运行基础数据和可行域顶点枚举算法对所述虚拟电厂出力可行域模型进行计算得到虚拟电厂出力可行域。
[0061] 本发明实施例的基于约束聚合的虚拟电厂出力可行域辨识装置，通过将含有分布式电源和需求响应资源的配电网等值为一个虚拟电厂参与输电网的优化调度运行；基于静
态安全约束模型，通过约束聚合构建虚拟电厂出力可行域模型；通过一种顶点枚举算法，计
算虚拟电厂出力可行域的边界顶点，从而得到描述虚拟电厂出力可行域的线性不等式约
束。将虚拟电厂内各类分布式资源的运行约束以及配电网安全运行约束聚合为关于虚拟电
厂有功和无功出力的显式约束条件，可以刻画安全运行允许的虚拟电厂出力可行域，使得
虚拟电厂可以像同步发电机一样参与电力系统优化调度运行和电力现货市场出清，确保调
度和出清结果不会违反虚拟电厂内部的安全运行约束，进而可提升系统的安全性和经济
性。
[0062] 另外，根据本发明上述实施例的基于约束聚合的虚拟电厂出力可行域辨识装置还可以具有以下附加的技术特征：
[0063] 进一步，在本发明的一个实施例中，所述虚拟电厂内部配电系统运行基础数据，包括：
[0064] 配网连接拓扑、配电线路的电导和电纳参数、变压器的电导和电纳参数、配电线路传输容量极限、安全运行允许的节点电压上下限、配网内部分布式电源有功和无功出力上
下限、配网内部负荷节点的有功和无功负荷预测。
[0065] 进一步，在本发明的一个实施例中，所述虚拟电厂静态安全约束模型，包括：
[0066]
[0067]
[0068]
[0069]
[0070]
[0071]
[0072]
[0073] 其中，式(1)和式(2)分别为节点i和节点j之间的配电线路的有功和无功潮流方程约束， 分别为配电线路的有功和无功潮流，bij、gij分别为配电电路的串联电纳、串
联电导，vi、vj分别为节点i和节点j电压幅值，θi、θj分别为节点i和节点j电压相角， 是节
点编号集合， 表示与节点i通过配电线路直接相连的节点编号的集合；
g
[0074] 式(3)和式(4)分别是节点有功和无功平衡方程约束，Pi 、 分别是节点i处分布d d
式发电的有功出力和无功出力，Pi、Qi分别是节点i的有功负荷和无功负荷，P0、Q0分别是虚
拟电厂的有功和无功出力；
[0075] 式(5)表示节点i和节点j之间的配电线路、配电变压器的传输容量约束，式中分别是节点i和节点j之间的配电线路和配电变压器的视在功率传输上限，N是对传输容
量约束进行分段线性化的分段数；
[0076] 式(6)表示配网节点电压上下限约束，式中 分别为节点i允许的电压幅值上限和下限；
[0077] 式(7)为分布式发电的有功和无功出力上下限约束，式中 是节点i的分布式发电的有功出力上下限， 是节点i的分布式发电的无功出力的上下限。
[0078] 进一步，在本发明的一个实施例中，所述第二构建模块进一步用于：
[0079] 根据所述虚拟电厂静态安全约束模型确定关于(v,θ,Pg,Qg,P0,Q0)的虚拟电厂静态安全约束可行域记Ω：
[0080] Ω＝{(v,θ,Pg,Qg,P0,Q0)|s.t.(1)～(7)} (8)g g g
[0081] 其中，v、θ、P、Q分别是以vi、θi、Pi、 为元素的向量；
[0082] 消去虚拟电厂内部状态变量v、θ和内部决策变量Pg、Qg，将虚拟电厂静态安全约束聚合为关于虚拟电厂有功和无功出力的约束，得到虚拟电厂出力可行域模型：
[0083]
[0084] 其中，Φ为所述虚拟电厂出力可行域。
[0085] 进一步，在本发明的一个实施例中，所述计算模块，包括：
[0086] 设定单元，用于设定所述虚拟电厂出力可行域的误差限值ε；
[0087] 初始顶点与多边形计算单元，用于根据所述虚拟电厂内部配电系统运行基础数据计算初始顶点集合和初始近似多边形；
[0088] 所述初始顶点与多边形计算单元具体用于，所述虚拟电厂出力可行域Φ的初始顶点集包含4个可行域顶点，通过求解如下优化问题计算得到：
[0089]
[0090] 其中， 是优化问题目标函数系数，对 依次取值(1,1)、(1,‑1)、 (‑1,1)、(‑1,‑1)，通过非线性优化算法或调用非线性优化求解器对式(10)进行求解，
对应的优化问题最优解记为 vm为第m个初始顶点；将初始顶点构成集合 V0
＝{v1,v2,v3,v4}，将可行域顶点集合初始化为V＝V0；将初始顶点构成的多边形各个边界的
T
方程的集合记为F；对V0中顶点计算平均值，即
[0091] 新增边界初始化单元，用于初始化新增边界集合FN为空集，即
[0092] 顶点搜索单元，用于通过向外平移FT中的边界搜索所述虚拟电厂出力可行域Φ的T (k)
边界点，对于F中的第k个边界，边界的方程为f : 边界两端的顶点
分别为 顶点搜索的计算模型如下：
[0093]
[0094] 通过非线性优化算法或调用非线性优化求解器求解顶点搜索的计算模型，求解得到顶点最优解记为 顶点最优值记为
[0095] 加入单元，用于将顶点最优解加进可行域顶点集合V；
[0096] 第一判断单元，根据设定的所述虚拟电厂出力可行域的误差限值和所述顶点搜索单元中的顶点搜索的改进量判断是否需要增加待平移边界；
[0097] 所述第一判断单元，具体用于，顶点搜索的改进量为：
[0098](k) (k) (k) (k)
[0099] 在Δgk>ε时，分别计算vk与 确定的直线的方程f1 、f2 ，将f1 、f2 加N T
入新增边界集合F，返回所述顶点搜索单元，使用F中的第k+1个边界再进行顶点搜索；
[0100] 在Δgk≤ε时，返回所述顶点搜索单元，使用FT中的第k+1个边界进行顶点搜索；
[0101] 第二判断单元，对FT中的所有边界方程进行顶点搜索后，判断FN是否为空集；
[0102] 如果FN不为空集，将待搜索边界集FT更新为FN，即FT←FN，并返回所述新增边界初始T
化单元对更新后的F中的边界进行顶点搜索；
[0103] 如果FN为空集，则结束，根据顶点集合V中的顶点确定多边形，根据确定的多边形得到所述虚拟电厂出力可行域。
[0104] 本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
[0105] 本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
[0106] 图1为根据本发明一个实施例的基于约束聚合的虚拟电厂出力可行域辨识方法流程图；
[0107] 图2为根据本发明一个实施例的通过顶点枚举算法计算虚拟电厂出力可行域的具体流程图；
[0108] 图3为根据本发明又一个实施例的基于约束聚合的虚拟电厂出力可行域辨识方法流程图；
[0109] 图4为根据本发明一个实施例的基于约束聚合的虚拟电厂出力可行域辨识装置结构示意图。
具体实施方式
[0110] 下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附
图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
[0111] 下面参照附图描述根据本发明实施例提出的基于约束聚合的虚拟电厂出力可行域辨识方法及装置。
[0112] 首先将参照附图描述根据本发明实施例提出的基于约束聚合的虚拟电厂出力可行域辨识方法。
[0113] 图1为根据本发明一个实施例的基于约束聚合的虚拟电厂出力可行域辨识方法流程图。
[0114] 如图1所示，该基于约束聚合的虚拟电厂出力可行域辨识方法包括以下步骤：
[0115] 步骤S1，获取虚拟电厂内部配电系统运行基础数据。
[0116] 具体地，虚拟电厂是指将含有分布式光伏、分布式风电、小型燃气轮机、需求响应资源等分布式发电设施的配电系统等值为一个发电厂，可以根据电网调度机构的调度指令
或电力现货市场的出清结果调节配电系统在并网关口处与输电网交换的有功和无功功率。
[0117] 获取的虚拟电厂内部配电系统运行基础数据包括：虚拟电厂内部配电系统运行基础数据包括配网连接拓扑、配电线路的电导和电纳参数、变压器的电导和电纳参数、配电线
路传输容量极限、安全运行允许的节点电压上下限、配网内部分布式电源有功和无功出力
上下限、配网内部负荷节点的有功和无功负荷预测。
[0118] 步骤S2，根据虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型。
[0119] 在获取虚拟电厂内部配电系统运行基础数据之后，根据虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型。
[0120] 虚拟电厂静态安全约束模型，表达式如下：
[0121]
[0122]
[0123]
[0124]
[0125]
[0126]
[0127]
[0128] 其中，式(1)和式(2)分别为节点i和节点j之间的配电线路ij的有功和无功潮流方程约束， 分别为配电线路的有功和无功潮流，bij、gij分别为配电电路的串联电纳、
串联电导，vi、vj分别为节点i和节点j电压幅值，θi、θj分别为节点i和节点j电压相角， 是
节点编号集合， 表示与节点i通过配电线路直接相连的节点编号的集合；
[0129] 式(3)和式(4)分别是节点有功和无功平衡方程约束，Pig、 分别是节点i处分布d
式发电的有功出力和无功出力，Pi、 分别是节点i的有功负荷和无功负荷，P0、Q0分别是虚
拟电厂的有功和无功出力；
[0130] 式(5)表示配电线路ij、配电变压器的传输容量约束，式中 分别是配电线路 ij和配电变压器的视在功率传输上限，N是对传输容量约束进行分段线性化的分段数；
[0131] 式(6)表示配网节点电压上下限约束，式中 分别为节点i允许的电压幅值上限和下限；
[0132] 式(7)为分布式发电的有功和无功出力上下限约束，式中 是节点i的分布式发电的有功出力上下限， 是节点i的分布式发电的无功出力的上下限。
[0133] 步骤S3，根据虚拟电厂静态安全约束模型构建虚拟电厂出力可行域模型。
[0134] 具体地，根据上面得到的虚拟电厂静态安全约束模型，及约束聚合理论构建电厂出力可行域模型，包括：
[0135] S31，根据虚拟电厂静态安全约束模型(式(1)‑式(7))确定关于(v,θ,Pg,Qg,P0,Q0) 的虚拟电厂静态安全约束可行域记Ω：
[0136] Ω＝{(v,θ,Pg,Qg,P0,Q0)|s.t.(1)～(7)} (8)
[0137] 其中，v、θ、Pg、Qg分别是以vi、θi、Pig、 为元素的向量；
[0138] S32，消去虚拟电厂内部状态变量v、θ和内部决策变量Pg、Qg，将虚拟电厂静态安全约束聚合为关于虚拟电厂有功和无功出力的约束，得到虚拟电厂出力可行域模型：
[0139]
[0140] 其中，Φ为虚拟电厂出力可行域，其是二维空间中的一个有界区域；该模型将虚拟电厂内各类分布式资源的运行约束以及配网安全运行约束聚合为关于虚拟电厂有功和无
功出力的可行域，其可以保证任何满足约束Φ的虚拟电厂出力(P0,Q0)均可以被虚拟电厂内
部的配电系统执行且不会违反静态安全约束。
[0141] 步骤S4，根据虚拟电厂内部配电系统运行基础数据和可行域顶点枚举算法对虚拟电厂出力可行域模型进行计算得到虚拟电厂出力可行域。
[0142] 采用一种可行域顶点枚举算法计算步骤S3中构建的虚拟电厂出力可行域：该可行域顶点枚举算法通过向外平移已有近似多边形的边、搜索Φ的新顶点，并使用新搜索得到
的顶点更新近似多边形，由此不断逼近真实的虚拟电厂出力可行域；如图2所示，使用该顶
点枚举算法计算虚拟电厂出力可行域的具体流程叙述如下：
[0143] S41，算法的输入包括虚拟电厂内部配电系统运行基础数据、预先设定的虚拟电厂出力可行域的误差限值ε(典型取值为0.01～0.05，可根据计算的精度需求设定)。
[0144] S42，根据虚拟电厂内部配电系统运行基础数据计算初始顶点集合和初始近似多边形。
[0145] 虚拟电厂出力可行域Φ的初始顶点集包含4个可行域顶点，通过求解如下优化问题计算得到：
[0146]
[0147] 其中， 是优化问题目标函数系数，对 依次取值(1,1)、(1,‑1)、 (‑1,1)、(‑1,‑1)，通过非线性优化算法或调用非线性优化求解器对式(10)进行求解，对应
的优化问题最优解记为 vm为第m个初始顶点；将初始顶点构成集合 V0＝
{v1,v2,v3,v4}，将可行域顶点集合初始化为V＝V0；将初始顶点构成的多边形各个边界的方
T
程的集合记为F；对V0中顶点计算平均值，即
[0148] S43，初始化新增边界集合FN为空集，即
[0149] S44，顶点搜索：通过向外平移FT中的边界搜索虚拟电厂出力可行域Φ的边界点，T (k)
对于F 中的第k个边界，边界的方程为f : 边界两端的顶点分别为
顶点搜索的计算模型如下：
[0150]
[0151] 通过非线性优化算法(如原对偶内点法)或调用非线性优化求解器(如IPOPT)求解顶点搜索的计算模型，求解得到顶点最优解记为 顶点最优值记为
[0152] S45，更新顶点集：将顶点最优解加进可行域顶点集合V。
[0153] S46，根据设定的虚拟电厂出力可行域的误差限值和步骤S44中的顶点搜索的改进量判断是否需要增加待平移边界。
[0154] S461，步骤S44中顶点搜索的改进量为：
[0155](k) (k) (k)
[0156] S462，在Δgk>ε时，分别计算vk与 确定的直线的方程f1 、f2 ，将f1 、 (k) N T
f2 加入新增边界集合F，返回步骤S44，使用F中的第k+1个边界再进行顶点搜索；
[0157] S463，在Δgk≤ε时，返回步骤S44，使用FT中的第k+1个边界进行顶点搜索；
[0158] S47，判断算法是否终止：对FT中的所有边界方程进行顶点搜索后，通过判断FN是否为空集判断算法是否终止；
[0159] S471，如果FN不为空集，即 将待搜索边界集FT更新为FN，即FT←FN，并返回T
步骤S43对更新后的F中的边界进行顶点搜索；
N
[0160] S472，如果F 为空集，即 则算法终止，根据顶点集合V中的顶点确定多边形，根据确定的多边形得到虚拟电厂出力可行域。
[0161] 顶点集合V中的顶点确定的多边形即为虚拟电厂出力可行域的计算结果，该多边形的各边界不等式即为表征虚拟电厂出力可行域的线性不等式约束。
[0162] 如图3所示，为本发明实施例的整体流程图，通过获取虚拟电厂内部配电系统运行基础数据，构建基于交流潮流方程的虚拟电厂静态安全约束模型，再基于约束聚合理论构
建虚拟电厂出力可行域模型，采用一种可行域顶点枚举算法计算虚拟电厂出力可行域。
[0163] 通过上述方法，填补了在虚拟电厂参与输电网和现货市场的优化过程中无法有效考虑虚拟电厂出力可行域的技术空白。将虚拟电厂内各类分布式资源的运行约束以及配网
安全运行约束聚合为关于虚拟电厂有功和无功出力的显式约束条件，刻画安全运行允许的
虚拟电厂出力可行域，从而使得虚拟电厂可以像同步发电机一样参与电力系统优化调度运
行和电力现货市场出清，可保证对虚拟电厂的调度和出清结果不会违反虚拟电厂安全运行
约束，有利于提升电力系统运行的安全性、经济性。
[0164] 其次参照附图描述根据本发明实施例提出的基于约束聚合的虚拟电厂出力可行域辨识装置。
[0165] 图4为根据本发明一个实施例的基于约束聚合的虚拟电厂出力可行域辨识装置结构示意图。
[0166] 如图4所示，该基于约束聚合的虚拟电厂出力可行域辨识装置包括：获取模块100、第一构建模块200、第二构建模块300和计算模块400。
[0167] 获取模块100，用于获取虚拟电厂内部配电系统运行基础数据。
[0168] 第一构建模块200，用于根据虚拟电厂内部配电系统运行基础数据构建基于交流潮流方程的虚拟电厂静态安全约束模型。
[0169] 第二构建模块300，用于根据虚拟电厂静态安全约束模型构建虚拟电厂出力可行域模型。
[0170] 计算模块400，用于根据虚拟电厂内部配电系统运行基础数据和可行域顶点枚举算法对虚拟电厂出力可行域模型进行计算得到虚拟电厂出力可行域。
[0171] 该装置使得虚拟电厂可以像同步发电机一样参与电力系统优化调度运行和电力现货市场出清，可保证对虚拟电厂的调度和出清结果不会违反虚拟电厂安全运行约束，有
利于提升电力系统运行的安全性、经济性。
[0172] 进一步，在本发明的一个实施例中，虚拟电厂内部配电系统运行基础数据，包括：
[0173] 配网连接拓扑、配电线路的电导和电纳参数、变压器的电导和电纳参数、配电线路传输容量极限、安全运行允许的节点电压上下限、配网内部分布式电源有功和无功出力上
下限、配网内部负荷节点的有功和无功负荷预测。
[0174] 进一步，在本发明的一个实施例中，虚拟电厂静态安全约束模型，包括：
[0175]
[0176]
[0177]
[0178]
[0179]
[0180]
[0181]
[0182] 其中，式(1)和式(2)分别为节点i和节点j之间的配电线路ij的有功和无功潮流方程约束， 分别为配电线路的有功和无功潮流，bij、gij分别为配电电路的串联电纳、
串联电导，vi、vj分别为节点i和节点j电压幅值，θi、θj分别为节点i和节点j电压相角， 是
节点编号集合， 表示与节点i通过配电线路直接相连的节点编号的集合；
[0183] 式(3)和式(4)分别是节点有功和无功平衡方程约束， 分别是节点i处分布d
式发电的有功出力和无功出力，Pi、 分别是节点i的有功负荷和无功负荷，P0、Q0分别是虚
拟电厂的有功和无功出力；
[0184] 式(5)表示配电线路ij、配电变压器的传输容量约束，式中 分别是配电线路 ij和配电变压器的视在功率传输上限，N是对传输容量约束进行分段线性化的分段数；
[0185] 式(6)表示配网节点电压上下限约束，式中 分别为节点i允许的电压幅值上限和下限；
[0186] 式(7)为分布式发电的有功和无功出力上下限约束，式中 是节点i的分布式发电的有功出力上下限， 是节点i的分布式发电的无功出力的上下限。
[0187] 进一步，在本发明的一个实施例中，第二构建模块进一步用于：
[0188] 根据虚拟电厂静态安全约束模型确定关于(v,θ,Pg,Qg,P0,Q0)的虚拟电厂静态安全约束可行域记Ω：
[0189] Ω＝{(v,θ,Pg,Qg,P0,Q0)|s.t.(1)～(7)} (8)g g g
[0190] 其中，v、θ、P、Q分别是以vi、θi、Pi、 为元素的向量；
[0191] 消去虚拟电厂内部状态变量v、θ和内部决策变量Pg、Qg，将虚拟电厂静态安全约束聚合为关于虚拟电厂有功和无功出力的约束，得到虚拟电厂出力可行域模型：
[0192]
[0193] 其中，Φ为虚拟电厂出力可行域。
[0194] 进一步，在本发明的一个实施例中，计算模块，包括：
[0195] 设定单元，用于设定虚拟电厂出力可行域的误差限值ε；
[0196] 初始顶点与多边形计算单元，用于根据虚拟电厂内部配电系统运行基础数据计算初始顶点集合和初始近似多边形；
[0197] 初始顶点与多边形计算单元具体用于，虚拟电厂出力可行域Φ的初始顶点集包含4个可行域顶点，通过求解如下优化问题计算得到：
[0198]
[0199] 其中， 是优化问题目标函数系数，对 依次取值(1,1)、(1,‑1)、 (‑1,1)、(‑1,‑1)，通过非线性优化算法或调用非线性优化求解器对式(10)进行求解，
对应的优化问题最优解记为 vm为第m个初始顶点；将初始顶点构成集合 V0
＝{v1,v2,v3,v4}，将可行域顶点集合初始化为V＝V0；将初始顶点构成的多边形各个边界的
T
方程的集合记为F；对V0中顶点计算平均值，即
[0200] 新增边界初始化单元，用于初始化新增边界集合FN为空集，即
[0201] 顶点搜索单元，用于通过向外平移FT中的边界搜索虚拟电厂出力可行域Φ的边界T (k)
点，对于F中的第k个边界，边界的方程为f : 边界两端的顶点分别
为 顶点搜索的计算模型如下：
[0202]
[0203] 通过非线性优化算法或调用非线性优化求解器求解顶点搜索的计算模型，求解得到顶点最优解记为 顶点最优值记为
[0204] 加入单元，用于将顶点最优解加进可行域顶点集合V；
[0205] 第一判断单元，根据设定的虚拟电厂出力可行域的误差限值和顶点搜索单元中的顶点搜索的改进量判断是否需要增加待平移边界；
[0206] 第一判断单元，具体用于，顶点搜索的改进量为：
[0207](k) (k) (k) (k)
[0208] 在Δgk>ε时，分别计算vk与 确定的直线的方程f1 、f2 ，将f1 、f2 加N T
入新增边界集合F，返回顶点搜索单元，使用F中的第k+1个边界再进行顶点搜索；
[0209] 在Δgk≤ε时，返回顶点搜索单元，使用FT中的第k+1个边界进行顶点搜索；
[0210] 第二判断单元，对FT中的所有边界方程进行顶点搜索后，判断FN是否为空集；
[0211] 如果FN不为空集，将待搜索边界集FT更新为FN，即FT←FN，并返回新增边界初始化单T
元对更新后的F中的边界进行顶点搜索；
[0212] 如果FN为空集，则结束，根据顶点集合V中的顶点确定多边形，根据确定的多边形得到虚拟电厂出力可行域。
[0213] 需要说明的是，前述对基于约束聚合的虚拟电厂出力可行域辨识方法实施例的解释说明也适用于该实施例的装置，此处不再赘述。
[0214] 根据本发明实施例提出的基于约束聚合的虚拟电厂出力可行域辨识装置，通过将含有分布式电源和需求响应资源的配电网等值为一个虚拟电厂参与输电网的优化调度运
行；基于静态安全约束模型，通过约束聚合构建虚拟电厂出力可行域模型；通过一种顶点枚
举算法，计算虚拟电厂出力可行域的边界顶点，从而得到描述虚拟电厂出力可行域的线性
不等式约束。将虚拟电厂内各类分布式资源的运行约束以及配电网安全运行约束聚合为关
于虚拟电厂有功和无功出力的显式约束条件，可以刻画安全运行允许的虚拟电厂出力可行
域，使得虚拟电厂可以像同步发电机一样参与电力系统优化调度运行和电力现货市场出
清，确保调度和出清结果不会违反虚拟电厂内部的安全运行约束，进而可提升系统的安全
性和经济性。
[0215] 此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者
隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三
个等，除非另有明确具体的限定。
[0216] 在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特
点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不
必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任
一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技
术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结
合和组合。
[0217] 尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述
实施例进行变化、修改、替换和变型。
