一种区域电量传输可行域确定方法及系统

摘要

本发明涉及一种区域电量传输可行域确定方法及系统，该方法包括以下步骤：S1、构建区域电网计及新能源不确定性的机组组合约束模型；S2、基于步骤S1构建的模型，采用全枚举法确定计及新能源不确定性及机组离散运行特性的区域电量传输可行域。该方法及系统有利于自动、准确地获取区域电量传输可行域。

权利要求

1.一种区域电量传输可行域确定方法，其特征在于，包括以下步骤：S1、构建区域电网计及新能源不确定性的机组组合约束模型；

S2、基于步骤S1构建的模型，采用全枚举法确定计及新能源不确定性及机组离散运行特性的区域电量传输可行域；

所述步骤S1中，考虑联络线全天电量，建立区域电网计及新能源不确定性的机组组合约束，包括系统约束、联络线约束、单机约束以及不确定性变量的区间表征，所述机组组合约束即为机组组合约束模型；

所述系统约束包括功率平衡约束和线路潮流约束，所述功率平衡约束为：式中，pG,t、pB,t、pD,t、pR,t分别为t时刻的机组出力、联络线功率、负荷和新能源出力；1G、

1B、1D、1R分别为pG,t、pB,t、pD,t、pR,t维度相匹配的元素均为1的列向量；T为调度时段集合；

所述线路潮流约束为：

式中，S为转移分布因子矩阵；AG、AB、AD、AR分别为pG,t、pB,t、pD,t、pR,t的节点关联矩阵；

和F分别为线路潮流的上限和下限；

所述联络线约束包括：

a.第b条联络线全天电量和各时段联络线功率的关系式中，B为边界联络线集合；

b.联络线功率容量约束

式中， 和PB分别为联络线传输上限和下限；

c.交割曲线约束

qB＝f(pB)           (5)式中，qB为所有时刻所有联络线传输电量的变量集合；pB为所有时刻所有联络线传输功率的变量集合；f表征交割曲线上电力与电量的关系；

所述单机约束包括：

a.系统启停约束

式中：ug,t为表征机组g在时刻t是否处于运行状态的0‑1变量；vg,t为表征机组g在时刻t是否由关闭状态变为开启状态的0‑1变量；G为机组集合；

b.机组容量约束

式中， 和Pg分别为机组g出力上限和下限；

c.机组最小启停时间约束

式中，Lg和lg分别为机组g的最小启动时间和最小关闭时间；

d.机组爬坡约束

式中，Rg为机组g的爬坡速率；Vg为机组g开启或关闭速率；

所述不确定性变量的区间表征为：

式中， 和PR,t分别为新能源出力的预测值上限和下限；

综合所有约束条件，得到如下表达式：式中：uG为所有机组在所有时刻是否处于运行状态的0‑1变量集合；vG为所有机组在所有时刻是否处于启动的0‑1变量集合；pG为所有机组在所有时刻的出力变量集合；pR为所有新能源场站在所有时刻的出力集合； 和PR分别为的预测值上限和下限；式(13a)等价描绘约束(1)‑(11)，其中Bq、Bu、Bv、Bp和BR为系数矩阵，C为常数向量；式(13b)即表征约束；

所述步骤S2的具体实现方法为：

将计及新能源不确定性及机组离散运行特性的区域电量传输可行域表述为：电量若属于可行域R，则当新能源出力在 范围内变化时，均存在一组满足式(13a)式约束的可行解，即：因此，若给定机组启停状态j，记其满足式(13a)的约束为：j

BqqB+BppG+BRpR≤C      (15)记机组启停状态j下的电量交易可行域为：则R和Rj存在如下关系：

R＝∪j∈JRj        (17)式中，J为所有可能的机组启停状态集合；

在机组状态j下，需满足下述约束：

则式所示的可等价表征为：

综上，得出计及新能源不确定性及机组离散运行特性的区域电量传输可行域即为：所有机组启停状态下、所有新能源极端出力场景下的省间电量交易可行空间的并集。

2.一种区域电量传输可行域确定系统，其特征在于，包括存储器、处理器以及存储于存储器上并能够被处理器运行的计算机程序指令，当处理器运行该计算机程序指令时，能够实现如权利要求1所述的方法步骤。

说明书

一种区域电量传输可行域确定方法及系统

技术领域

[0001] 本发明属于电力技术领域，具体涉及一种区域电量传输可行域确定方法及系统。

背景技术

[0002] 随着电力负荷的快速增长和大规模新能源并网，单一区域电网已难以实现源‑荷平衡。为实现电力资源的广域优化配置，互联电网通过联络线实现功率传输，提升电力资源利用效率。在我国新能源大规模并网的背景下，随着复奉、锦苏等特高压跨区域直流工程的相继投入，我国跨区电力资源交换不断增强。跨区电力资源交换需以保证区域内电网运行安全为基本前提，但现有方法难以满足精准描绘区域电量传输可行域的工业要求：1)考虑以机组组合约束为表征的离散运行特性；2)计及以区间数形式表征的新能源不确定性。

发明内容

[0003] 本发明的目的在于提供一种区域电量传输可行域确定方法及系统，该方法及系统有利于自动、准确地获取区域电量传输可行域。

[0004] 为实现上述目的，本发明采用的技术方案是：一种区域电量传输可行域确定方法，包括以下步骤：

[0005] S1、构建区域电网计及新能源不确定性的机组组合约束模型；

[0006] S2、基于步骤S1构建的模型，采用全枚举法确定计及新能源不确定性及机组离散运行特性的区域电量传输可行域。

[0007] 进一步地，所述步骤S1中，考虑联络线全天电量，建立某区域电网计及新能源不确定性的机组组合约束，包括系统约束、联络线约束、单机约束以及不确定性变量的区间表征，所述机组组合约束即为机组组合约束模型。

[0008] 进一步地，所述系统约束包括功率平衡约束和线路潮流约束，所述功率平衡约束为：

[0009]

[0010] 式中，pG,t、pB,t、pD,t、pR,t分别为t时刻的机组出力、联络线功率、负荷和新能源出力；1G、1B、1D、1R分别为pG,t、pB,t、pD,t、pR,t维度相匹配的元素均为1的列向量；T为调度时段集合；

[0011] 所述线路潮流约束为：

[0012]

[0013] 式中，S为转移分布因子矩阵；AG、AB、AD、AR分别为pG,t、pB,t、pD,t、pR,t的节点关联矩阵； 和F分别为线路潮流的上限和下限。

[0014] 进一步地，所述联络线约束包括：

[0015] a.第b条联络线全天电量和各时段联络线功率的关系

[0016]

[0017] 式中，B为边界联络线集合；

[0018] b.联络线功率容量约束

[0019]

[0020] 式中， 和PB分别为联络线传输上限和下限；

[0021] c.交割曲线约束

[0022] qB＝f(pB)   (5)

[0023] 式中，qB为所有时刻所有联络线传输电量的变量集合；pB为所有时刻所有联络线传输功率的变量集合；f表征交割曲线上电力与电量的关系。

[0024] 进一步地，所述单机约束包括：

[0025] a.系统启停约束

[0026]

[0027] 式中：ug,t为表征机组g在时刻t是否处于运行状态的0‑1变量；vg,t为表征机组g在时刻t是否由关闭状态变为开启状态的0‑1变量；G为机组集合；

[0028] b.机组容量约束

[0029]

[0030] 式中， 和Pg分别为机组g出力上限和下限；

[0031] c.机组最小启停时间约束

[0032]

[0033]

[0034] 式中，Lg和lg分别为机组g的最小启动时间和最小关闭时间；

[0035] d.机组爬坡约束

[0036]

[0037]

[0038] 式中，Rg为机组g的爬坡速率；Vg为机组g开启或关闭速率。

[0039] 进一步地，所述不确定性变量的区间表征为：

[0040]

[0041] 式中， 和PR,t分别为新能源出力的预测值上限和下限。

[0042] 进一步地，综合所有约束条件，得到如下表达式：

[0043]

[0044] 式中：uG为所有机组在所有时刻是否处于运行状态的0‑1变量集合；vG为所有机组在所有时刻是否处于启动的0‑1变量集合；pG为所有机组在所有时刻的出力变量集合；pR为所有新能源场站在所有时刻的出力集合； 和PR分别为的预测值上限和下限；式(13a)等价描绘约束(1)‑(11)，其中Bq、Bu、Bv、Bp和BR为系数矩阵，C为常数向量；式(13b)即表征约束。

[0045] 进一步地，所述步骤S2的具体实现方法为：

[0046] 将计及新能源不确定性及机组离散运行特性的区域电量传输可行域表述为：电量若属于可行域R，则当新能源出力在 范围内变化时，均存在一组满足式(13a)式约束的可行解，即：

[0047]

[0048] 因此，若给定机组启停状态j，记其满足式(13a)的约束为：

[0049] BqqB+BppG+BRpR≤Cj  (15)

[0050] 记机组启停状态j下的电量交易可行域为：

[0051]

[0052] 则R和Rj存在如下关系：

[0053] R＝∪j∈JRj  (17)

[0054] 式中，J为所有可能的机组启停状态集合；

[0055] 在机组状态j下，需满足下述约束：

[0056]

[0057] 则式所示的可等价表征为：

[0058]

[0059] 综上，得出计及新能源不确定性及机组离散运行特性的区域电量传输可行域即为：所有机组启停状态下、所有新能源极端出力场景下的省间电量交易可行空间的并集。

[0060] 本发明还提供了一种区域电量传输可行域确定系统，包括存储器、处理器以及存储于存储器上并能够被处理器运行的计算机程序指令，当处理器运行该计算机程序指令时，能够实现如权利要求1‑8所述的方法步骤。

[0061] 相较于现有技术，本发明具有以下有益效果：基于考虑新能源不确定性的机组组合约束，提出基于枚举新能源极端出力场景和机组启停状态的区域电量传输可行域确定方法，从而可以自动、快速地获得计及新能源不确定性及机组离散运行特性的区域电量传输可行域，提高了获取的区域电量传输可行域的精度。

附图说明

[0062] 图1是本发明实施例的方法实现原理图。

[0063] 图2是本发明实施例中区域电量传输可行域计算流程图。

[0064] 图3是本发明实施例中获得的区域电量传输可行域。

具体实施方式

[0065] 下面结合附图及实施例对本发明做进一步说明。

[0066] 应该指出，以下详细说明都是示例性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。

[0067] 需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。

[0068] 如图1所示，本实施例提供了一种区域电量传输可行域确定方法，包括以下步骤：

[0069] S1、构建区域电网计及新能源不确定性的机组组合约束模型。

[0070] 所述步骤S1中，考虑联络线全天电量，建立某区域电网计及新能源不确定性的机组组合约束，包括系统约束、联络线约束、单机约束以及不确定性变量的区间表征，所述机组组合约束即为机组组合约束模型。

[0071] 所述系统约束包括功率平衡约束和线路潮流约束，所述功率平衡约束为：

[0072]

[0073] 式中，pG,t、pB,t、pD,t、pR,t分别为t时刻的机组出力、联络线功率、负荷和新能源出力；1G、1B、1D、1R分别为pG,t、pB,t、pD,t、pR,t维度相匹配的元素均为1的列向量；T为调度时段集合；

[0074] 所述线路潮流约束为：

[0075]

[0076] 式中，S为转移分布因子矩阵；AG、AB、AD、AR分别为pG,t、pB,t、pD,t、pR,t的节点关联矩阵； 和F分别为线路潮流的上限和下限。

[0077] 所述联络线约束包括：

[0078] a.第b条联络线全天电量和各时段联络线功率的关系

[0079]

[0080] 式中，B为边界联络线集合；

[0081] b.联络线功率容量约束

[0082]

[0083] 式中， 和PB分别为联络线传输上限和下限；

[0084] c.交割曲线约束

[0085] qB＝f(pB)  (5)

[0086] 式中，qB为所有时刻所有联络线传输电量的变量集合；pB为所有时刻所有联络线传输功率的变量集合；f表征交割曲线上电力与电量的关系。

[0087] 所述单机约束包括：

[0088] a.系统启停约束

[0089]

[0090] 式中：ug,t为表征机组g在时刻t是否处于运行状态的0‑1变量；vg,t为表征机组g在时刻t是否由关闭状态变为开启状态的0‑1变量；G为机组集合；

[0091] b.机组容量约束

[0092]

[0093] 式中， 和Pg分别为机组g出力上限和下限；

[0094] c.机组最小启停时间约束

[0095]

[0096]

[0097] 式中，Lg和lg分别为机组g的最小启动时间和最小关闭时间；

[0098] d.机组爬坡约束

[0099]

[0100]

[0101] 式中，Rg为机组g的爬坡速率；Vg为机组g开启或关闭速率。

[0102] 所述不确定性变量的区间表征为：

[0103]

[0104] 式中， 和PR,t分别为新能源出力的预测值上限和下限。

[0105] 综合所有约束条件，得到如下表达式：

[0106]

[0107] 式中：uG为所有机组在所有时刻是否处于运行状态的0‑1变量集合；vG为所有机组在所有时刻是否处于启动的0‑1变量集合；pG为所有机组在所有时刻的出力变量集合；pR为所有新能源场站在所有时刻的出力集合； 和PR分别为的预测值上限和下限；式(13a)等价描绘约束(1)‑(11)，其中Bq、Bu、Bv、Bp和BR为系数矩阵，C为常数向量；式(13b)即表征约束。

[0108] S2、基于步骤S1构建的模型，采用全枚举法确定计及新能源不确定性及机组离散运行特性的区域电量传输可行域。

[0109] 所述步骤S2的具体实现方法为：

[0110] 将计及新能源不确定性及机组离散运行特性的区域电量传输可行域表述为：电量若属于可行域R，则当新能源出力在 范围内变化时，均存在一组满足式(13a)式约束的可行解，即：

[0111]

[0112] 因此，若给定机组启停状态j，记其满足式(13a)的约束为：

[0113] BqqB+BppG+BRpR≤Cj  (15)

[0114] 记机组启停状态j下的电量交易可行域为：

[0115]

[0116] 则R和Rj存在如下关系：

[0117] R＝∪j∈JRj  (17)

[0118] 式中，J为所有可能的机组启停状态集合；

[0119] 在机组状态j下，需满足下述约束：

[0120]

[0121] 则式所示的可等价表征为：

[0122]

[0123] 综上，得出计及新能源不确定性及机组离散运行特性的区域电量传输可行域即为：所有机组启停状态下、所有新能源极端出力场景下的省间电量交易可行空间的并集。

[0124] 本实施例还提供了用于实现上述方法的区域电量传输可行域确定系统，其特征在于，包括存储器、处理器以及存储于存储器上并能够被处理器运行的计算机程序指令，当处理器运行该计算机程序指令时，能够实现如权利要求1‑8所述的方法步骤。

[0125] 以IEEE 9节点系统为例，验证本发明方法的正确性。

[0126] 在节点5和节点9分别连接一条联络线，在节点1接入新能源机组1台，并以2时段时间尺度为例。

[0127] 在2时段IEEE 9节点测试系统中，共有两台火电机组和1台新能源机组。根据推论可知：计及新能源不确定性及机组离散运行特性的区域电量传输可行域为，所有机组启停状态下、所有新能源极端出力场景下的区域电量传输可行域的并集。因此，共需枚举机组状2 2 1*2
态2*2＝16个，新能源极端场景2 ＝2个。

[0128] 1)可行域的确定

[0129] 在16种机组状态中，仅有2种机组启停状态下的电量交易可行域为非空集。采用图2所示方法，可得计及新能源不确定性及机组离散运行特性的区域电量传输可行域如图3所示。

[0130] 2)可行域的有效性证明

[0131] 为验证图3所得可行域有效性，本专利采用蒙特卡洛法进行验证。随机抽取图3区域内的10000个样本，并对新能源出力在其预测区间内同时随机采样。对于每个样本，总存在一组可行解使得约束(13a)可行的，即验证了本专利所提推论的正确性。

[0132] 本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD‑ROM、光学存储器等)上实施的计算机程序产品的形式。

[0133] 本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

[0134] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

[0135] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

[0136] 以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。