一种基于电压调控的虚拟电厂分区构建模型及构建方法与流程

[0001]
本发明公开了一种基于电压调控的虚拟电厂分区构建模型及构建方法，属于电力系统调压技术领域。


背景技术：

[0002]
近年来，随着经济发展电力负荷特性恶化，负荷增速高于电量增速，加之分布式风电和光伏的发展，电网电源结构改变，配电网运行方式也随之发生变化，部分节点易出现电压越限问题。随着虚拟电厂(vpp)技术的发展，用户侧分布式风电和光伏以及可调负荷资源通过聚合构建虚拟电厂参与电网运行，由于分布式风电和光伏本身配置有无功设备，存在无功出力调节空间，可帮助调节部分节点电压，当节点有功负荷发生变化时同样会影响节点的电压，如何基于实际电网分区情况优化构建分区vpp来实现电压调整尚未得到很好的解决。


技术实现要素：

[0003]
本发明为了优化构建分区vpp实现电压调整，提供了一种基于电压调控的虚拟电厂分区构建模型及构建方法，该模型或方法能够优化构建分区vpp，助力电力系统电压调整。
[0004]
本发明所采取的技术方案为：一种基于电压调控的虚拟电厂分区构建模型，包括上层模型和下层模型，所述下层模型在以最小化调用成本为目标的基础上调用分区内节点vpp
m
的调节容量，基于对各节点vpp
m
容量的优化配置构建分区vpp，并将各分区vpp的配置输送至所述上层模型，所述上层模型在下层模型各分区vpp构建的基础上以整体潮流不均衡最小为目标函数优化各线路潮流改变，并向所述下层模型反馈电力系统中整体潮流变化。
[0005]
进一步的，所述下层模型的构建方法为：
[0006]
s01根据网络实时量测的有功、无功、电压、相角等数据和电阻、电抗等系统拓扑结构参数通过潮流计算得到各分区雅克比矩阵；
[0007]
s02根据各分区雅克比矩阵获得各分区电压灵敏度矩阵；
[0008]
s03构建vpp
m
，以各个节点的有功和无功可调资源构建节点虚拟电厂，即vppm，并根据vpp
m
出力范围确定各节点电压改变量；
[0009]
s04根据各节点电压改变量获得调压后各节点的节点电压，且调压后的各节点电压需满足节点电压要求；
[0010]
s05以区域内各节点vpp
m
的总调用成本最小为目标函数以所述节点电压要求与各节点虚拟电厂vpp
m
的有功和无功调节容量为约束条件构建分区vpp。
[0011]
进一步的，所述上层模型以潮流熵为指标衡量系统中线路潮流分布的均衡度，通过所述潮流熵得到潮流不均衡严重度函数定量衡量系统中线路潮流分布的均衡度。
[0012]
进一步的，步骤s01中雅克比矩阵具体为
[0013][0014]
式中，p为各节点注入的有功功率，q为各节点注入的无功功率，v为各节点的电压幅值，δ为各节点的相角。
[0015]
进一步的，步骤s02中电压灵敏度矩阵为雅克比矩阵的逆矩阵，所述电压灵
[0016]
敏度矩阵为
[0017]
式中，p为各节点注入的有功功率，q为各节点注入的无功功率，v为各节点的电压幅值，δ为各节点的相角。
[0018]
进一步的，步骤s03中节点电压改变量范围为
[0019][0020]
式中，i和m均为节点编号(i＝1,2,k,n；m＝1,2,k,n)；和q
vpp,m
分别为节点vpp
m
的有功调节容量和无功调节容量。
[0021]
进一步的，步骤s04中调压后各节点需满足的节点电压要求为
[0022][0023]
u
i,min
≤u
i
≤u
i,max
[0024]
式中，为调压前节点i的电压，u
i
为调压后节点i的电压，u
i,min
为节点i的电压上限；u
i,max
为节点i的电压下限，δv
i
为调压前和调压后节点电压的变化量。
[0025]
进一步的，步骤s05中的目标函数为：
[0026][0027]
式中，f为区域内各节点vpp
m
的总调用费用；c
load
、c
wind
和c
pv
分别为调用可调负荷、风电机组无功和分布式光伏无功的单位成本；δp
load,m
、q
wind,m
和q
pv,m
分别为节点vpp
m
内部电设备的虚拟有功出力、风电机组和分布式光伏的无功出力。
[0028]
进一步的，所述潮流不均衡严重度函数获取方法为：
[0029]
a计算线路l的负载率μ
l
(t)
[0030][0031]
式中，p
l0
(t)为t时段线路l的潮流值；p
lmax
为线路l的最大功率传输容量；m为线路总条数；
[0032]
b计算vpp参与调压后系统的潮流熵h
[0033][0034]
c计算潮流熵最大值h
max
[0035][0036]
d潮流不均衡严重度函数为
[0037][0038]
式中，h为vpp参与调压后系统的潮流熵；h0为电力系统初始状态下的潮流熵；h
max
为潮流熵最大值。
[0039]
本发明中一种基于电压调控的虚拟电厂分区构建方法，包括构建上层模型和构建下层模型，所述下层模型在以最小化调用成本为目标的基础上调用分区内节点vpp
m
的调节容量，基于对各节点vpp
m
容量的优化配置构建分区vpp，并将各分区vpp的配置输送至所述上层模型，所述上层模型以整体潮流不均衡最小为目标函数优化各线路潮流改变，并向所述下层模型反馈电力系统中整体潮流变化。
[0040]
本发明所产生的有益效果包括：
[0041]
(1)vpp调用的是现有资源，可实现资源的充分利用；
[0042]
(2)采用vpp调压可以基于现有资源完成电压调整，减少无功电源的投资；
[0043]
(3)采用双层优化模型可以在完成调压的同时降低系统潮流的不平衡度。
附图说明
[0044]
图1 vpp参与调压前节点电压；
[0045]
图2 vpp参与调压后节点电压；
[0046]
图3风电机组基本结构示意图；
[0047]
图4基于电压调控的分区vpp构建结构图。
具体实施方式
[0048]
下面结合附图和具体实施例对本发明做进一步详细的解释说明，但应当理解为本发明的保护范围并不受具体实施例的限制。
[0049]
(1)vpp调压原理
[0050]
vpp通过聚合用户侧分布式电源和可调负荷等资源形成，由于其内部资源的多样性，可以通过调整不同设备的有功或无功出力来帮助系统进行调压。
[0051]
系统运行中节点b进行调压前如图1所示，即
[0052][0053]
式中，u
a
(t)为节点a的电压；u
b
(t)为节点b的电压；p
b
(t)和q
b
(t)分别为节点b处的有功功率和无功功率；r和x分别为线路的电阻和电抗。
[0054]
vpp由于内部资源丰富即可发出有功也可发出无功，当其接入节点b对其进行调压时，其接入后如图2所示
[0055]
vpp接入节点b参与调压后：
[0056][0057]
式中，为vpp参与调压后的节点b的电压；p
vpp
(t)和q
vpp
(t)分别为vpp提供的有功功率和无功功率。
[0058]
(2)vpp提供调压能力
[0059]
1)无功出力设备
[0060]
a.分布式风电机组
[0061]
以常见的双馈感应风机为例进行分析，风机基本结构如图3所示，
[0062]
图3中p
m
为风电机组输入机械功率，由捕获风力机获得的风能大小决定；p
s
和q
s
分别为定子的注入有功和无功功率；p
g
和q
g
分别为风电机组注入电网的有功和无功功率，p
r
和q
r
为转子侧变换器注入转子的有功和无功功率；p
c
和q
c
为网侧变换器从电网输入的有功和无功功率。
[0063]
忽略定转子电阻，转子侧最大电流约束下，定子侧无功调节范围为
[0064][0065]
式中，l
s
和l
m
分别为定子电感和激磁电感；i
r
为定转子电流有效值；s为转差率，s＝(ω
1-ω
r
)/ω1；ω1和ω
r
分别为同步旋转角速度和转子旋转角速度；u
s
为定子电压有效值。
[0066]
考虑定子侧最大电流i
s,max
限制，定子侧无功调节范围为
[0067][0068]
综合上述两种情况，分布式风机定子侧无功调节范为
[0069][0070]
网侧换流器无功调节能力主要受限于换流器容量s
c,max
：
[0071][0072]
则网侧换流器无功极限为
[0073][0074]
分布式双馈风电机组无功调节极限为
[0075][0076]
b.分布式光伏
[0077]
光伏逆变器可利用自身无功控制功能为配电系统提供电压支撑，可调无功范围与逆变器容量关系为
[0078][0079]
式中，为逆变器最大无功输出容量；p
pv
为光伏逆变器有功出力；s
inv
为逆变器容量，约为额定有功容量的1.0～1.1倍。
[0080]
2)有功出力资源
[0081]
虚拟电厂主要是通过削减内部用户用电负荷，如空调、热水器等可调性负荷，以此来实现对外的虚拟有功出力，其有功出力可调范围为
[0082][0083]
式中，δp
max
为最大有功可调量；p
b
为未参与vpp时的用电负荷；为参与vpp后最小的用电负荷。
[0084]
首先以各分区内存在有功或无功调节潜力的节点构建节点vpp，综上可得每个节点vpp的整体出力为
[0085]
s
vpp,m
＝δp
load,m
+j(q
wind,m
+q
pv,m
)
[0086]
δp
load,m
≤δp
max,m
[0087]
q
min,m
≤q
wind,m
≤q
max,m
[0088][0089]
式中，s
vpp,m
为节点vpp
m
出力视在功率；δp
load,m
、q
wind,m
和q
pv,m
分别为节点vpp
m
内部用电设备的虚拟有功出力、风电机组和分布式光伏的无功出力。
[0090]
(3)一种基于电压调控的虚拟电厂分区构建模型
--
vpp调压双层模型
[0091]
下层模型
[0092]
下层模型在以最小化调用成本为目标的基础上调用分区内节点vpp的调节容量，基于对各节点vpp容量的优化配置构建分区vpp。
[0093]
1)计算各分区雅克比矩阵
[0094]
电力系统中节点功率变化与母线电压和相角的关系为
[0095][0096][0097]
式中，v、δ为各节点电压幅值与相角；p、q为各节点注入的有功、无功功率。δp和δq分别为用功功率和无功功率改变量；δδ和δv节点电影相角和幅值的改变量。
[0098]
根据式(1)和式(2)得到雅克比矩阵j。
[0099]
2)计算各分区电压灵敏度矩阵
[0100]
各节点的电压灵敏度取决于网络阻抗，即由线路参数和拓扑结构决定。通过改变vpp所在节点注入的有功或无功功率来调节该节点的电压幅值或相角。系统中各节点的电压灵敏度可以由潮流计算中雅克比矩阵j的逆阵获得，则电压灵敏矩阵可表述为
[0101][0102]
式中，v、δ为各节点电压幅值与相角；p、q为各节点注入的有功、无功功率。灵敏度子矩阵分别表示由于vpp输出的有功和无功变化所导致的各节点电压幅值变化。每个灵敏度子矩阵都是n
×
n，其中n表示(分区内)配电网中节点数。
[0103]
3)构建vpp
m
，并根据vpp
m
出力范围确定各节点电压改变量；
[0104]
当具备调节能力的节点m上节点vpp
m
有功功率与无功功率输出发生变化分别为p
vpp,m
和q
vpp,m
时，系统中各节点电压变化量表示为
[0105][0106]
式中，i和m均为节点编号(i＝1,2,k,n；m＝1,2,k,n),v
i
为节点i的电压幅值，p
m
为节点m注入的有功功率；p
vpp,m
和q
vpp,m
分别为节点vpp
m
的有功调节容量和无功调节容量，和通过电压灵敏度矩阵获得。
[0107]
由上式可以看出vpp可对电压灵敏度较高的节点提供较好的调压服务，对电压灵敏度低的节点，vpp的调压效果很小。
[0108]
4)获得调压后各节点电压
[0109][0110]
调压后各节点应满足节点电压要求
[0111]
u
i,min
≤u
i
≤u
i,max
[0112]
式中，和u
i
分别为调压前后节点i的电压；u
i,min
和u
i,max
分别为节点i的电压上下限。
[0113]
5)构建目标函数
[0114]
下层以区域内节点vpp调用成本最小化为目标函数
[0115][0116]
式中，f为区域内调用节点vpp的费用；c
load
、c
wind
和c
pv
分别为调用可调负荷、风电机组无功和分布式光伏无功的单位成本。
[0117]
2)上层模型
[0118]
潮流分布是否均衡对连锁故障的传播有很大的影响，因此是系统潮流状态分析需关注的重要参量之一。在下层模型完成区域节点vpp容量配置的基础上，由于各分区节点vpp接入容量的不同，系统各区域实时功率变化也不同，进而导致系统潮流改变情况不同，上层以整体潮流不均衡最小为目标函数
[0119]
mins(h)
[0120]
式中，s(h)为潮流不均衡度。
[0121]
本专利以潮流熵指标来定量衡量系统中线路潮流分布的均衡度。线路l的负载率为
[0122][0123]
式中，p
l0
(t)为t时段线路l的潮流值；p
lmax
为线路l的最大功率传输容量；m为线路总条数。
[0124]
给定一个常数数列a＝{a1,a2,k,a
k
,k,a
n
}，用l
k
表示负载率μ
i
∈(a
k
,a
+k1
的输电线路条数，把各个负载区间内的线路条数比例化得：
[0125][0126]
由此，系统的潮流熵为
[0127][0128]
当系统中所有支路的负载率都在同一个区间时，h为零，这时系统的潮流分布处于最均衡状态，即线路承载的潮流与其容量成正比；而当系统中所有线路的负载率都不在同一个区间时，h便达到最大值：
[0129][0130]
此时电网的潮流分布最不均衡，某些线路负载率很低，某些线路负载率很高，当有扰动出现时，高负载率的线路容易过载而引发系统故障。
[0131]
则潮流不均衡严重度函数可表示为
[0132][0133]
式中，h为vpp参与调压后系统的潮流熵；h0为系统初始状态下的潮流熵；h
max
为潮流熵最大值。
[0134]
本发明通过分析用户侧分布式风电和光伏以及可调负荷资源的有功和无功调节潜力形成节点vpp，进一步构建调压双层优化模型，下层模型以电压灵敏度为基础以调用成本最小为目标函数构建分区vpp，实现电压的分区调整，最后上层模型在下层各分区vpp构建的基础上以系统整体潮流不均衡最小为目标函数与下层目标函数进行迭代优化，在尽量保证成本最小和潮流平衡的情况下完成分区vpp的构建，实现电压调控。
[0135]
本发明中一种基于电压调控的虚拟电厂分区构建方法，包括构建上层模型和构建下层模型，上层模型的构建与下层模型的构建同基于电压调控的虚拟电厂分区构建模型中所述。
[0136]
已经通过参考少量实施方式描述了本发明。然而，本领域技术人员所公知的，正如附带的专利权利要求所限定的，出了本发明以上公开的其他的实施例等同的落在本发明的范围内。
[0137]
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0138]
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0139]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0140]
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。