一种计及规模储能的含风电电力系统黑启动分区方法
1.本发明涉及新型电力系统恢复控制领域，具体为一种计及规模储能的含风电电力系统黑启动分区方法。
背景技术：
2.为实现“双碳”目标，推动能源转型，中央提出构建以新能源为主体的新型电力系统。以风电、光伏为代表的新能源虽然具有启动速度快、启动功率小等优势，但因其固有的出力间歇性、波动性给电网的安全运行造成了极大的挑战，大停电事故时有发生，比如2019年的委内瑞拉大停电、2021年的美国德州大停电。所以研究大电网大停电后的黑启动有着重要的意义。
3.目前我国“三华”特高压同步电网形成“五纵五横”主网架，电网规模比2010年翻一番以上。在大电网、新型电力系统背景下越来越多的区域已经实现了电网互联，提高了供电的可靠性、更合理的利用系统中各类发电厂提高运行经济性。由于互联大电网一般不只有一个黑启动电源且系统规模巨大所以首先按照大电网中黑启动电源的分布情况对系统进行合理分区，各个分区在时间上并行恢复，提高了整体恢复效率，降低了各分区子系统的停电时间。所以科学合理的分区对系统快速完成网架重构有着重要意义。
技术实现要素：
4.(一)解决的技术问题
5.现有研究主要针对传统电力系统分区，对于大量配置储能的新型电力系统分区研究较为薄弱。而在新型电力系统背景下为了应对新能源出力的波动性需要系统中由足够的灵活性资源作支撑，而在黑启动阶段电力系统骨干网架尚未恢复，系统比较脆弱，灵活性资源匮乏难以接入规模新能源。储能是系统中最重要的灵活性资源，为构建新型电力系统需要在系统中配置大量储能。为此需要考虑各个分区的灵活性支撑能力。
6.所以本发明首先使用两阶段优化模型求解分区方案，第一阶段以待恢复发电机获得启动功率时间最短对于系统内的风电场储能资源以及其他待恢复机组进行分区；第二阶段以分区内部联系紧密为目标对系统中剩余的负荷节点进行分区得到分区方案。之后在第一阶段优化模型中添加“割”约束去掉已有分区方案重新优化得到备选分区方案集合；随后使用innc法求解各个分区方案的ppf曲线对各个分区进行准重构，得到各个分区的评价指标值，使用变异系数法筛选出最优分区方案。
7.(二)技术方案
8.本发明提供一种计及规模储能的含风电电力系统黑启动分区方法，包括以下步骤：
9.s1：待恢复机组分区
10.进一步的，s1所述详细步骤如下：
11.s101：依据电力系统黑启动电源的位置与数量确定机组分区约束。
12.s102：以待恢复发电机得到启动功率时间最短为目标，考虑机组分区约束、机组数量约束对系统内的风电场-特级负荷联合系统、风电场及其储能资源、其余机组节点进行分区。本发明以发电机距黑启动电源的电气距离代表待恢复发电机得到启动功率时间。目标函数f1的表示如下式所示:
[0013][0014]
式中，m为分区总数；n
gh,i
为分区i待恢复发电机数量；n
wh,i
为分区i待恢复普通风电场数量；n
ws,i
为分区i风电场-特级负荷系统数量；im表示第i个分区内的黑启动电源节点。d
i,j,im
为节点j到分区i黑启动电源im的电气距离；d
i,j,iml
，d
i,j,imws
分别为风电场-特级负荷系统j中特级负荷和风储系统到im的电气距离。z
i,j
为节点分区决策变量，z
i,j
＝1表示将节点j分到第i个分区。
[0015]
s2：备选分区方案集合求解
[0016]
进一步的，s2所述详细步骤如下：
[0017]
s201：根据每个分区内的火电机组厂用电和黑启动电源的发电能力确定机组启动约束、源荷匹配约束、区内联通约束、以及联合系统的联合运行约束。
[0018]
s202：以分区节点到分区黑启动电源的电气距离最小作为目标函数为目标，考虑上述约束对系统剩余储能、负荷节点进行分区。
[0019][0020]
式中，m为分区总数；n
l,i
为第i个分区的节点数；im表示第i个分区内的黑启动电源节点。
[0021]
s203：在第一阶段的机组分区模型中重置分区模型，本发明添加如下约束去掉已经得到的分区方案重新优化，如此循环直到无法求解。
[0022][0023]
其中：γ为己确定的黑启动分区方案集合，随着优化进行该集合元素数量逐渐增加；γ为某一具体子系统划分方案；z
γim
为方案γ中变量的具体取值。
[0024]
如此循环直到无法求解。得到的方案组成备选分区方案集合。上述过程添加约束之后得到的分区方案至少有一台待恢复机组、储能的分区与已有方案不同，而负荷的分区根据机组的分区不同相应的改变。
[0025]
将得到的方案组成备选分区方案集合。
[0026]
s3：分区方案准重构结果求解
[0027]
进一步的，s3所述详细步骤如下：
[0028]
s301：对分区的机组、储能建立模型。具体包括：火电机组模型、风电机组模型、储能资源模型。
[0029]
s302：确定潮流约束、暂态频率约束约束。
[0030]
本发明采用直流潮流作为约束。
[0031][0032]
式中b
l
、b0分别为支路电纳矩阵和节点电纳矩阵；pb为节点注入有功向量；p
l
为支路功率矩阵。
[0033]
s303：以最大化火电机组出力、最大化火电机组出力为目标，使用innc法求解分区备选准重构方案集合
[0034]
(1)最大化火电机组出力
[0035]
在最大化火电机组出力的目标引导下，黑启动机组或已经恢复的初期带电系统的启动功率将全部用来恢复火电机组的厂用负荷，机组启动决策力图使整个重构时长t期间火电机组的发电量最大，目标函数可表示为，
[0036][0037]
式中，n
gh
为已恢复的火电机组总数。n
gd
为待启动的火电机组总数；k
j,t
表示t时步第j台火电机组是否启动的0/1决策变量；t
0,i
为已恢复机组的启动时刻，t0为当前时刻。
[0038]
(2)最大化风电机组出力
[0039]
为了最大化风电机组出力，可用的启动功率需要按照风电出力的预测波动程度为其提供最大可能的有功备用支持，不再恢复火电机组辅机设备。t为本时步时长，其目标函数为
[0040][0041]
式中第1项表示该时步决策启动的台风电机组按照期望出力的总发电量，其中变量含义同式(18)。第2项表示在ty内已达到最小技术出力的相关火电机组为风电提供的同步备用电量。在网架重构阶段，厂用电负荷是最重要的负荷，第3项表示发电功率优先恢复厂用电负荷。
[0042]
以上述两个目标为目标函数，以潮流约束、暂态频率约束、生产模型作为约束条件组成双目标优化模型按照innc法的步骤进行求解。
[0043]
s304：使用变异系数法考虑风电恢复量、火电恢复量、重要负荷恢复量选择最优准重构方案，并用其结果作为该分区方案的评价指标值。
[0044]
s4：最优分区方案选择
[0045]
新型电力系统要求系统有足够的灵活性消纳新能源。本发明选择准重构阶段的风电恢复量、火电恢复量两个指标表征系统的灵活性。并选择分区电气距离、系统恢复时间、分区恢复时间平衡性三个指标表征分区方案的合理性，使用变异系数法综合考虑以上五个指标，根据各个方案的各个指标值之差使用变异系数法求得各个指标的权重，之后对各方案的指标值进行加权求和，获得方案的综合评价值。综合评价值最大的方案为最优分区方案。
[0046]
(三)有益效果
[0047]
与现有技术相比，本发明提供了一种计及规模储能的含风电电力系统黑启动分区
方法，具备以下有益效果：
[0048]
1、本发明的分区方法可以根据系统的实际情况合理分区，使得网络中的所有发电机组以及灵活性资源合理分配，提高其利用率，使得分区之后的方案在网架重构阶段中恢复的系统发电能力最大。
[0049]
2、相较于传统的网架重构分区，本发明考虑了储能这种最重要的灵活性资源，使得本发明接入的风电比例大大提高，对于新型电力系统背景下的网架重构分区提出了一种参考方法；
[0050]
3、本发明将黑启动风电场与其最近恢复的火电机组组成扩展黑启动电源使得分区方案以及准重构过程更加符合实际；将配置储能但达不到黑启动的风电场与其最近的特级负荷组成风电场-特级负荷系统，大大减少了特级负荷的停电时间减少了大停电带来的经济损失。
附图说明
[0051]
图1为本发明流程示意图；
[0052]
图2为图1火电厂出力曲线。
具体实施方式
[0053]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0054]
请参阅图2，本发明提供一种计及规模储能的含风电电力系统黑启动分区方法：包括以下步骤：
[0055]
s1：待恢复机组分区
[0056]
第一阶段是以待恢复发电机获得启动功率时间最短为目标，考虑机组归属约束、分区机组数量约束对系统内的风电场-特级负荷联合系统、风电场及其储能资源、其余机组节点进行分区。
[0057]
s101：依据电力系统黑启动电源的位置与数量确定机组分区约束。机组分区约束即每个机组只能归属于一个分区
[0058]z1,j
+z
2,j
+
…
+z
i,j
…
+z
m,j
＝1,j∈{1,2,
…
,n}
[0059]zi,j
为节点分区决策变量，z
i,j
＝1表示将节点j分到第i个分区。
[0060]
机组数量约束即机组不能集中于一个分区，如下式所示：
[0061][0062]
s102：以待恢复发电机得到启动功率时间最短为目标，考虑机组分区约束、机组数量约束对风电场-特级负荷联合系统、风电场及其储能资源、其余机组节点。第一阶段以待恢复发电机得到启动功率时间最短为目标进行求解。本发明以发电机距黑启动电源的电气距离代表其恢复时间。目标函数f1的表示如下式所示:
[0063][0064]
式中，m为分区总数；n
gh,i
为分区i待恢复发电机数量；n
wh,i
为分区i待恢复普通风电场数量；n
ws,i
为分区i风电场-特级负荷系统数量；im表示第i个分区内的黑启动电源节点。d
i,j,im
为节点j到分区i黑启动电源im的电气距离；d
i,j,iml
，d
i,j,imws
分别为风电场-特级负荷系统j中特级负荷和风储系统到im的电气距离。z
i,j
为节点分区决策变量，z
i,j
＝1表示将节点j分到第i个分区。
[0065]
s2：备选分区方案集合求解
[0066]
第二阶段以分区内部联系紧密为目标对系统中剩余储能、负荷节点进行分区得到一个分区方案。本发明以分区节点到分区黑启动电源的电气距离最小作为目标函数表示分区内部联系紧密。考虑各种分区约束，建立分区优化模型对系统剩余负荷节点进行分区得到一个分区方案。
[0067]
之后在第一阶段优化模型重置分区模型并去掉已经生成的方案重新求解。反复求解直到无可行分区方案或者分区方案足够多时停止求解，将得到的所有分区方案组成备选分区方案集合。
[0068]
s201：根据每个分区内的火电机组厂用电和黑启动电源的发电能力确定机组启动约束、源荷匹配约束、区内联通约束、以及联合系统的联合运行约束。
[0069]
(1)机组启动约束。
[0070]
对于分区i，其待恢复机组中至少一台机组的启动功率要小于黑启动机组额定输出功率的70％，其余30％额定功率供给其附近的负荷，以保证其稳定运行，即
[0071][0072]
式中：p
bi
为分区i中的黑启动机组额定功率；p
gij，f
为分区i中第j台待恢复机组所需的启动功率。
[0073]
(2)源荷匹配约束
[0074]
分区时需要保证每个分区内均有足够的可调度负荷支撑发电机组安全运行，即分区内火电机组最小出力与风电机组期望出力之和小于等于可调度负荷总量:
[0075][0076]
式中，n
g,m
为第m个分区内的火电机组数量；n
w,k,m
为第m个分区内的风电场数量；n
d,m
第m个分区内的负荷节点数量；p
min,m,j
为第m个分区内第j台火电机组的最小技术出力；为第m个分区内第k个风电场的期望出力；p
dl,m
为第m个分区内第l个负荷节点所带的负荷；
[0077]
(3)区内连通约束
[0078]
本发明使用网络流理论保证各个分区的连通性，需要满足如下条件：网络中仅存在一个源点汇点需有流量到达；满足基尔霍夫第一定律。为此需要将系统中的支路使用两
条方向相反的有向弧表示，构造成有向图。如果支路有向弧两端节点均在同一个分区内，则其一定属于该分区。
[0079]
(4)联合运行系统约束
[0080]
对于特级负荷联合运行系统，其储能配置不足。为了保证特级负荷不断电，需要保证在2h之内纳入稳定恢复的分区中，所以该联合运行系统与其所在分区的黑启动电源的距离不能太远。如式(21)所示：
[0081]di,j
≤d
max,j
ꢀꢀꢀ
(21)
[0082]di,j
表示分区i中第j个特级负荷联合运行系统到分区i中黑启动电源的距离；d
max,j
为第j个特级负荷联合运行系统到其所在分区黑启动电源的最大距离。
[0083]zim
＝z
(i-j)m
[0084]zim
为第i个特级负荷风电场的分区决策变量；z
(i-j)m
为该风电场供电的特级负荷j的分区决策变量。该约束表示特级负荷联合运行系统的两个子节点处于同一个分区。
[0085]
对于扩展黑启动系统，要求其恢复的火电厂处于该黑启动电源形成的分区内。如下式所示：
[0086]zjm
＝1
[0087]
zjm为风电场黑启动电源m直接供电的火电厂j的分区决策变量。
[0088]
s202：以分区节点到分区黑启动电源的电气距离最小作为目标函数为目标，考虑上述约束对系统剩余储能、负荷节点进行分区。
[0089][0090]
式中，m为分区总数；n
l,i
为第i个分区的节点数；im表示第i个分区内的黑启动电源节点。
[0091]
s203：在第一阶段的机组分区模型中重置分区模型，本发明添加如下约束去掉已经得到的分区方案重新优化，如此循环直到无法求解。
[0092][0093]
其中：γ为己确定的黑启动分区方案集合，随着优化进行该集合元素数量逐渐增加；γ为某一具体子系统划分方案；z
γim
为方案γ中变量的具体取值。
[0094]
如此循环直到无法求解。得到的方案组成备选分区方案集合。上述过程添加约束之后得到的分区方案至少有一台待恢复机组、储能的分区与已有方案不同，而负荷的分区根据机组的分区不同相应的改变。
[0095]
将得到的方案组成备选分区方案集合。
[0096]
s3：分区方案准重构结果求解
[0097]
本发明首先对系统中各种资源分别建立模型，同时考虑潮流约束、暂态频率约束、火电机组启动时限约束。之后以最大化火电机组出力、最大化火电机组出力为目标建立双目标准重构模型，之后使用innc法求解得到各个分区的准重构方案。最后以火电机组发电量、风电场发电量、重要负荷恢复量为指标选择最终准重构方案，得到各个分区方案的准重构结果。
[0098]
s301：对分区的机组、储能建立模型。具体包括：火电机组模型、风电机组模型、储
能资源模型。
[0099]
(1)火电机组
[0100]
同一重构时步内的不同机组可分别处于准备启动或启动、升负荷及达到最小技术出力等不同运行状态。火电机组恢复之后其出力如下图所示：
[0101]
网架重构阶段火电机组在达到最小技术出力之前不会预留备用，在达到最小技术出力之后可以牺牲部分出力作为上调备用恢复风电机组。需要注意，由于火电机组出力恢复速度受到其锅炉温度和压力的较大影响，为了最大限度地消除对重构效率的制约或不利影响，机组启动时刻t
ia
应满足式(13)所示的相应启动时间限制。
[0102]
0《t
ia
《t
ch,i
或t
ia
》t
cc,j
ꢀꢀꢀꢀ
(13)
[0103]
其中，t
ch,i
为自然循环机组的最大临界热启动时间，t
cc,j
为直流机组的最小临界启动时间。
[0104]
(2)风电机组
[0105]
当黑启动风电场供电火电机组出力为0时
[0106]
此时黑启动风电场所带储能只能维持本风电场出力，调度权受限。约束表示如下：
[0107][0108]
m表示风电场黑启动电源节点，此约束保证此时黑启动风电场储能的出力正好平抑本身的波动。其余备用约束如下：
[0109][0110]
其中nw表示的是该时步可恢复区域内的风电场数；nj为第j个风电场所拥有的风机台数。为第j个风电场本时步恢复的风机台数。为第j个风电场的期望出力；为第j个风电场在置信度水平下的最小出力；n
gm
为该时步已经达到最小技术出力可以预留上调备用的火电机组台数；p
gup,i,t
为第k台火电机组预留的上调备用容量。p
i,tfang
为t时步储能i的放电功率；ns为储能的个数。
[0111]
如果分区黑启动电源为风电场，其配置储能在其最先恢复的火电机组出力之前只能平抑本风电场出力波动；当其最先恢复的火电机组开始并网出力之后则其配置储能调度权限放开，为系统灵活性需求服务。具体模型如下：
[0112]
其最先恢复的火电机组没有并网出力，此时预留备用模型如下：
[0113][0114]
除此之外还要保证风电场m供电的火电厂必须恢复，如下式所示：
[0115][0116]
表示火电机组g在t时步的恢复决策变量。
[0117]
如果火电机组已经并网出力则预留备用模型如下：
[0118][0119]
(3)储能资源模型
[0120]
储能系统在t时刻的功率输出可以表示为：
[0121][0122]
其中，表示t时刻bess输出有功功率，表示t时刻bess的放电功率，表示t时刻bess的充电功率。如果为正，储能系统放电；否则，储能系统充电。
[0123]
(3.1)储能的功率约束
[0124][0125][0126]
上述两式分别表示储能的充、放电功率不超过其最大、最小功率限制。
[0127]yich
和y
ifang
分别为储能的充放电状态决策变量。y
ich
为1表示储能i充电、为0表示不充电。y
ifang
同理。
[0128]
由于同一时步储能的充放电状态只有一种故还要满足下面的约束
[0129][0130]
(3.2)储能容量约束
[0131][0132][0133]
ηi为储能i的转化效率。为储能i的最小容量一般取额定容量的10％，为储能i的最大容量一般取额定容量的90％。
[0134]
s302：确定潮流约束、暂态频率约束约束。
[0135]
(1)潮流约束
[0136]
本发明采用直流潮流作为约束。
[0137][0138]
式中b
l
、b0分别为支路电纳矩阵和节点电纳矩阵；pb为节点注入有功向量；p
l
为支路功率矩阵。
[0139]
(2)暂态频率偏差限值约束
[0140][0141]
δf
max
为系统允许的最大暂态频率偏差，通常取0.5hz；ks,i为已经并网的第i台机组的单位调节功率；pstr,t为t时步机组最大增发功率。
[0142]
s303：以最大化火电机组出力、最大化火电机组出力为目标，使用innc法求解分区备选准重构方案集合。
[0143]
(1)最大化火电机组出力
[0144]
在最大化火电机组出力的目标引导下，黑启动机组或已经恢复的初期带电系统的启动功率将全部用来恢复火电机组的厂用负荷，机组启动决策力图使整个重构时长t期间火电机组的发电量最大，目标函数可表示为式(20)。
[0145][0146]
式中，n
gh
为已恢复的火电机组总数。n
gd
为待启动的火电机组总数；k
j,t
表示t时步第j台火电机组是否启动的0/1决策变量；t
0,i
为已恢复机组的启动时刻，t0为当前时刻。
[0147]
(2)最大化风电机组出力
[0148]
为了最大化风电机组出力，可用的启动功率需要按照风电出力的预测波动程度为其提供最大可能的有功备用支持，不再恢复火电机组辅机设备。t为本时步时长，其目标函数如式(23)。
[0149][0150]
式中第1项表示该时步决策启动的台风电机组按照期望出力的总发电量，其中变量含义同式(18)。第2项表示在ty内已达到最小技术出力的相关火电机组为风电提供的同步备用电量。在网架重构阶段，厂用电负荷是最重要的负荷，第3项表示发电功率优先恢复厂用电负荷。
[0151]
以上述两个目标为目标函数，考虑以上约束组成双目标优化模型按照innc法的步骤进行求解。
[0152]
s304：使用变异系数法考虑风电恢复量、火电恢复量、重要负荷恢复量选择最优准重构方案，并用其结果作为该分区方案的评价指标值。
[0153]
s4：最优分区方案选择
[0154]
新型电力系统要求系统有足够的灵活性消纳新能源。本发明选择准重构阶段的风电恢复量、火电恢复量两个指标表征系统的灵活性。并选择分区电气距离、系统恢复时间、分区恢复时间平衡性三个指标表征分区方案的合理性，使用变异系数法综合考虑以上五个指标，根据各个方案的各个指标值之差使用变异系数法求得各个指标的权重，之后对各方案的指标值进行加权求和，获得方案的综合评价值。综合评价值最大的方案为最优分区方案。
[0155]
采用变异系数法对最优分区方案的决策过程简述如下：
[0156]
(1)备选分区方案集合a＝{a1,a2,
…
,ai,
…
,an}共n个方案；
[0157]
(2)利用火电出力恢复量、风电出力恢复量、分区电气距离、系统恢复时间、分区恢复时间平衡性构建评价指标集合e＝{e1,e2,e3,e4,e5}，记待评价方案ai针对评价指标ej的取值为y
ij
；
[0158]
(3)构建决策矩阵[y
ij
]n×5并依据式(30)进行归一化，得到规范化决策矩阵
[0159][0160]
(4)确定各评价指标的权重。首先，根据式(31)求出评价指标ej的变异系数vj。进而，采用式(32)求出该指标的权重wj。
[0161][0162][0163]
(5)确定最优重构方案。以规范化决策矩阵为基础，采用评价指标权重对各方案的指标值进行加权求和，获得方案的综合评价值。综合评价值最大的方案为时步的最优重构方案。
[0164]
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
