CN111313416A - 一种多源协同的智能配电网故障恢复次序优化决策方法 - Google Patents

一种多源协同的智能配电网故障恢复次序优化决策方法 Download PDF

Info

Publication number
CN111313416A
CN111313416A CN202010243290.8A CN202010243290A CN111313416A CN 111313416 A CN111313416 A CN 111313416A CN 202010243290 A CN202010243290 A CN 202010243290A CN 111313416 A CN111313416 A CN 111313416A
Authority
CN
China
Prior art keywords
node
line
passing
power
distribution network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010243290.8A
Other languages
English (en)
Other versions
CN111313416B (zh
Inventor
许寅
王颖
和敬涵
王小君
马佳骏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Jiaotong University
Original Assignee
Beijing Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Jiaotong University filed Critical Beijing Jiaotong University
Priority to CN202010243290.8A priority Critical patent/CN111313416B/zh
Publication of CN111313416A publication Critical patent/CN111313416A/zh
Application granted granted Critical
Publication of CN111313416B publication Critical patent/CN111313416B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供的一种多源协同的智能配电网故障恢复次序优化决策方法,包括根据配电网的节点分布情况,建立配电网故障恢复的目标函数和相应的优化变量;基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果。本发明提供的方法没通过计及故障恢复过程中系统三相不对称潮流约束以及其他各项稳态运行约束、辐射状拓扑约束和操作过程相关约束,提出更为完善的配电网故障恢复次序优化混合整数线性规划模型,形成了一种智能配电网故障恢复次序优化决策方法,以确保每次恢复操作后均能满足配电网各项稳态约束,确保系统的稳定运行。

Description

一种多源协同的智能配电网故障恢复次序优化决策方法
技术领域
本发明涉及配电网故障恢复技术领域,尤其涉及一种多源协同的智能配电网故障恢复次序优化决策方法。
背景技术
近年来世界范围内自然灾害、恐怖袭击等极端事件频发,对电力设施造成严重损害,并带来了巨大的经济损失。配电网韧性描述配电网在极端事件中采取主动措施保证重要负荷供电,并迅速响应恢复断电负荷的能力。随着分布式发电技术的不断发展,利用配电网内分布式资源,如分布式电源、储能系统和微电网等为重要负荷恢复供电已成为提升配电网韧性的有效方法之一。本发明重点关注利用配电网本地分布式资源恢复重要负荷的最优恢复操作次序的决策方法。
目前已有部分技术对恢复操作次序的确定方法进行了探索。例如在一些现有技术中利用改进粒子群算法和最短路径算法求解通过分布式电源恢复配电网的操作顺序。在另一些现有技术中将问题建模为混合整数线性规划,得到分布式电源与线路的一系列控制操作。上述现有技术的缺点在于,在其恢复过程中涉及同时闭合多个线路开关或同时连接多个负荷,与实际操作有所差别。
发明内容
本发明的实施例提供了一种多源协同的智能配电网故障恢复次序优化决策方法,用于解决申请人发现的现有技术中,未考虑三相不对称潮流约束条件,使实际恢复过程与实际操作有较大差别的技术问题。
为了实现上述目的,本发明采取了如下技术方案。
一种多源协同的智能配电网故障恢复次序优化决策方法,其特征在于,包括:
建立配电网故障恢复的目标函数,该目标函数具有优化变量;
基于目标函数,通过加入三相不对称潮流约束、稳态运行约束和过程操作约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果。
优选地,目标函数为:
Figure BDA0002433268680000021
优化变量为:公式(1)中,N为配电网中所有节点集合;Y为整个故障恢复过程中所有操作的集合;y表示某一次操作,第0次操作表示配电网故障隔离后恢复的初始状态;y的最大值为ymax=nload+nL+ngen+nes,其中nload表示带有开关的负荷数量,nL表示带开关的线路数量,ngen表示分布式电源的数量,nes表示储能系统的数量;ωi为节点i所连接的负荷的权重系数;
Figure BDA0002433268680000022
为整数变量,表示第y次操作后,节点i所连接的负荷是否被恢复;
Figure BDA0002433268680000023
表示负荷被恢复,否则负荷未被恢复。
优选地,基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果的步骤中,三相不对称潮流约束包括:
通过
Figure BDA0002433268680000024
表示流入节点i的功率之和等于流出节点i的总和;
通过
Figure BDA0002433268680000025
表示
Figure BDA0002433268680000026
Figure BDA0002433268680000027
之间的关系;
通过
Figure BDA0002433268680000028
表示线路(i,j)上的两端节点的电压差;
通过
Figure BDA0002433268680000029
表示线路(i,j)是否连接的两种情况,当
Figure BDA00024332686800000210
时,该式为支路首末两端电压模值平方差的关系,当
Figure BDA00024332686800000211
时,支路两端电压无约束;
通过
Figure BDA00024332686800000212
Figure BDA00024332686800000213
Figure BDA00024332686800000214
表明节点i的注入功率;若节点i与分布式电源连接,则节点i的注入功率等于连接到节点i的分布式电源的发电量减去节点i的负荷功率;若节点i与储能连接,则节点i的注入功率等于连接到节点i的储能系统的发电量减去节点i的负荷功率;若节点i不与分布式电源或储能系统连接,则节点i的注入功率等于负的节点i的负荷功率;
上述公式中,G为配电网中所有可控分布式电源的集合;E为配电网中所有储能系统集合;L为配电网中所有支路集合;函数DIAG(·)返回以输入向量为对角元素的对角矩阵;
Figure BDA0002433268680000031
表示第y次操作后从节点h流向i的复功率,为3×3维矩阵;
Figure BDA0002433268680000032
表示
Figure BDA0002433268680000033
的对角线元素,为3×1维向量;δ为线性化转换矩阵;
Figure BDA0002433268680000034
为中间辅助变量;
Figure BDA0002433268680000035
是节点i第y次操作后的电压,为3×1维向量;
Figure BDA0002433268680000036
为i节点第y次操作后的电压向量与其共轭转置相乘得到的3×3维矩阵,
Figure BDA0002433268680000037
是一个对角矩阵,其对角线元素为节点电压各相幅值的平方;zij表示线路(i,j)的3×3维阻抗矩阵;
Figure BDA0002433268680000038
表示第y次操作后节点i的注入功率,为3×1维向量,
Figure BDA0002433268680000039
表示节点i的负荷功率需求,为3×1维向量;
Figure BDA00024332686800000310
表示第y次操作后节点i分布式电源出力,为3×1维向量;
Figure BDA00024332686800000311
表示第y次操作后节点i储能出力,为3×1维向量;
Figure BDA00024332686800000312
表示第y次操作后线路状态,
Figure BDA00024332686800000313
表示线路被连接;否则线路未被连接;M为值很大的正实数。
优选地,基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果的步骤中,稳态运行约束包括运行安全约束和拓扑约束;过程操作约束包括初始状态约束;
运行安全约束包括:
通过
Figure BDA00024332686800000314
表示线路功率约束,具体表示线路(i,j)是否连接的两种情况,当
Figure BDA00024332686800000315
时,未对
Figure BDA00024332686800000316
进行限制,当
Figure BDA00024332686800000317
时,限制
Figure BDA00024332686800000318
使得其值为0;
通过
Figure BDA00024332686800000319
表示节点电压上下限约束,若节点没有通电,则限制其电压为0;
通过
Figure BDA0002433268680000041
Figure BDA0002433268680000042
表示节点电压上下限约束,若节点没有通电,则限制其电压为0;
通过
Figure BDA0002433268680000043
Figure BDA0002433268680000044
表示储能的容量约束;
公式(9)-(14)中,函数diag(·)返回由输入矩阵的对角元素组成的向量;函数Re(·)返回输入矩阵的实部;函数Im(·)返回输入矩阵的虚部;函数ones(·)返回维度是输入数的单位矩阵;
Figure BDA0002433268680000045
Figure BDA0002433268680000046
表示节点电压各相幅值最小值的平方与幅值最大值的平方,为3×1维向量;
Figure BDA0002433268680000047
表示第y次操作后节点i的通电状态,
Figure BDA0002433268680000048
表示节点i在第y次操作后通电;否则节点i没有通电;
Figure BDA0002433268680000049
表示第y次操作后分布式电源i的连接状态,
Figure BDA00024332686800000410
表示分布式电源i接入电网,否则其没有接入电网;
Figure BDA00024332686800000411
表示第y次操作后储能i的连接状态,
Figure BDA00024332686800000412
表示储能i接入电网,否则其没有接入电网;
Figure BDA00024332686800000413
表示分布式电源与储能系统的有功功率最大值和无功功率最大值,为3×1维向量;
拓扑约束包括:
通过
Figure BDA00024332686800000414
表示在图论的树的定义中,每一个可以形成环的线路组合中,至少有一条线路是断开的;
通过
Figure BDA00024332686800000415
表示树的连通性,表示树的必要不充分条件,所有通电节点之和减去所有通电线路之和小于等于1;
通过
Figure BDA00024332686800000416
Figure BDA00024332686800000417
Figure BDA00024332686800000418
表示线路通电与节点通电状态和线路连通状态的关系;
公式(15)-(20)中,K为配电网中全部可用线路连接后形成的环的集合;Ck为第k个环中所包含的线路的集合;
Figure BDA00024332686800000419
为第k个环中所包含的线路数量,
Figure BDA00024332686800000420
表示第y次操作后节点线路(i,j)是否通电,
Figure BDA00024332686800000421
表示线路(i,j)通电,否则表示线路(i,j)未通电;
过程操作约束包括:
通过
Figure BDA0002433268680000051
表示只有节点通电时,该节点的负荷才能够接入配电网;
通过
Figure BDA0002433268680000052
表示只有节点通电时,ES才能启动;
通过
Figure BDA0002433268680000053
表示除根节点之外的DG,只有在所连接节点通电时,才能启动;
通过
Figure BDA0002433268680000054
表示当具有黑启动能力且容量最大的DG接入配电网时,其所连接的节点也同时通电;
通过
Figure BDA0002433268680000055
表示节点若通电,就不能断电;
通过
Figure BDA0002433268680000056
Figure BDA0002433268680000057
表示DG、ES、负荷若接入系统就不能再断开;
通过
Figure BDA0002433268680000058
表示每次只进行一项操作,Lc,Lo分别为2种线路类型的集合,发生故障后,初始状态为闭合的线路集合与初始状态为断开的线路集合;
通过
Figure BDA0002433268680000059
表示对于初始状态为断开的线路,上一个操作后的线路连接状态总是闭合或与下一个操作后的线路连接状态相同;
通过
Figure BDA00024332686800000510
表示对于初始状态为闭合的线路,上一个操作后的线路连接状态总是断开或与下一个操作后的线路连接状态相同;
通过
Figure BDA00024332686800000511
表示对于当线路状态为连接时,线路两端的通电节点状态总是相同的;
通过
Figure BDA00024332686800000512
表示对于初始状态为断开的线路,当线路状态没有改变的时候,无约束,当线路状态改变时,在改变之前必然是两端节点有且只有一端是通电的;
通过
Figure BDA0002433268680000061
表示对于初始状态为闭合的线路,当线路状态没有改变的时候,无约束,当线路状态改变时,在改变之前必然是两端节点都没有通电的;
初始状态约束包括:
通过
Figure BDA0002433268680000062
Figure BDA0002433268680000063
表示将初始线路的状态设置为闭合或断开;
通过
Figure BDA0002433268680000064
表示在所有节点在初始状态下均未通电;
通过
Figure BDA0002433268680000065
表示当容量最大的黑启动DG启动时,其所连接的节点电压为恒定电压。
由上述本发明的实施例提供的技术方案可以看出,本发明提供的一种多源协同的智能配电网故障恢复次序优化决策方法,通过计及故障恢复过程中系统三相不对称潮流约束以及其他各项稳态运行约束、辐射状拓扑约束和操作过程相关约束,提出更为完善的配电网故障恢复次序优化混合整数线性规划模型,形成了一种智能配电网故障恢复次序优化决策方法,以确保每次恢复操作后均能满足配电网各项稳态约束,确保系统的稳定运行。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种多源协同的智能配电网故障恢复次序优化决策方法的处理流程图;
图2为本发明提供的一种多源协同的智能配电网故障恢复次序优化决策方法的测试系统拓扑图;
图3为本发明提供的一种多源协同的智能配电网故障恢复次序优化决策方法的恢复策略求解结果图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的任一单元和全部组合。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
为便于对本发明实施例的理解,下面将结合附图以几个具体实施例为例做进一步的解释说明,且各个实施例并不构成对本发明实施例的限定。
参见图1,本发明提供一种多源协同的智能配电网故障恢复次序优化决策方法,包括如下步骤:
建立配电网故障恢复的目标函数,该目标函数具有优化变量;
基于目标函数,通过加入三相不对称潮流约束、稳态运行约束和过程操作约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果。
进一步的,在本发明提供的实施例中,目标函数为:
Figure BDA0002433268680000081
优化变量为:公式(1)中,N为配电网中所有节点集合;Y为整个故障恢复过程中所有操作的集合;y表示某一次操作,第0次操作表示配电网故障隔离后恢复的初始状态;y的最大值为ymax=nload+nL+ngen+nes,其中nload表示带有开关的负荷数量,nL表示带开关的线路数量,ngen表示分布式电源的数量,nes表示储能系统的数量;ωi为节点i所连接的负荷的权重系数;
Figure BDA0002433268680000082
为整数变量,表示第y次操作后,节点i所连接的负荷是否被恢复;
Figure BDA0002433268680000083
表示负荷被恢复,否则负荷未被恢复。
进一步的,上述的基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果的步骤中,三相不对称潮流约束包括:
通过
Figure BDA0002433268680000084
表示流入节点i的功率之和等于流出节点i的总和;
通过
Figure BDA0002433268680000085
表示
Figure BDA0002433268680000086
Figure BDA0002433268680000087
之间的关系;
通过
Figure BDA0002433268680000088
表示通过大“M”方法线路(i,j)上的两端节点的电压差;
通过
Figure BDA0002433268680000089
表示线路(i,j)是否连接的两种情况,当
Figure BDA00024332686800000810
时,该式为支路首末两端电压模值平方差的关系,当
Figure BDA00024332686800000811
时,支路两端电压无约束;
上述公式(2)-(5)是用于三相不平衡潮流约束的线性近似模型;
通过
Figure BDA00024332686800000812
Figure BDA00024332686800000813
表明节点i的注入功率;若节点i与分布式电源连接,则节点i的注入功率等于连接到节点i的分布式电源的发电量减去节点i的负荷功率;若节点i与储能连接,则节点i的注入功率等于连接到节点i的储能系统的发电量减去节点i的负荷功率;若节点i不与分布式电源或储能系统连接,则节点i的注入功率等于负的节点i的负荷功率;
上述公式中,G为配电网中所有可控分布式电源(DG)的集合;E为配电网中所有储能系统集合;L为配电网中所有支路集合;函数DIAG(·)返回以输入向量为对角元素的对角矩阵;
Figure BDA0002433268680000091
表示第y次操作后从节点h流向i的复功率,为3×3维矩阵;
Figure BDA0002433268680000092
表示
Figure BDA0002433268680000093
的对角线元素,为3×1维向量;δ为线性化转换矩阵;
Figure BDA0002433268680000094
为中间辅助变量;
Figure BDA0002433268680000095
是节点i第y次操作后的电压,为3×1维向量;
Figure BDA0002433268680000096
为i节点第y次操作后的电压向量与其共轭转置相乘得到的3×3维矩阵,
Figure BDA0002433268680000097
实际上是一个对角矩阵,其对角线元素为节点电压各相幅值的平方;zij表示线路(i,j)的3×3维阻抗矩阵;
Figure BDA0002433268680000098
表示第y次操作后节点i的注入功率,为3×1维向量,
Figure BDA0002433268680000099
表示节点i的负荷功率需求,为3×1维向量;
Figure BDA00024332686800000910
表示第y次操作后节点i分布式电源出力,为3×1维向量;
Figure BDA00024332686800000911
表示第y次操作后节点i储能出力,为3×1维向量;
Figure BDA00024332686800000912
表示第y次操作后线路状态,
Figure BDA00024332686800000913
表示线路被连接;否则线路未被连接;M为值很大的正实数。
更进一步的,上述的基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果的步骤中,稳态运行约束包括运行安全约束和拓扑约束;过程操作约束还包括初始状态约束;
运行安全约束包括:
通过
Figure BDA00024332686800000914
表示线路功率约束,同样使用大
“M”方法,表示线路(i,j)是否连接的两种情况,当
Figure BDA00024332686800000915
时,未对
Figure BDA00024332686800000916
进行限制,当
Figure BDA00024332686800000917
时,限制
Figure BDA00024332686800000918
使得其值为0;
通过
Figure BDA00024332686800000919
表示节点电压上下限约束,若节点没有通电,则限制其电压为0;
通过
Figure BDA00024332686800000920
Figure BDA00024332686800000921
Figure BDA00024332686800000922
表示节点电压上下限约束,若节点没有通电,则限制其电压为0;
通过
Figure BDA00024332686800000923
Figure BDA00024332686800000924
Figure BDA00024332686800000925
表示储能的容量约束;
公式(9)-(14)中,函数diag(·)返回由输入矩阵的对角元素组成的向量;函数Re(·)返回输入矩阵的实部;函数Im(·)返回输入矩阵的虚部;函数ones(·)返回维度是输入数的单位矩阵;
Figure BDA00024332686800000926
Figure BDA00024332686800000927
表示节点电压各相幅值最小值的平方与幅值最大值的平方,为3×1维向量;
Figure BDA0002433268680000101
表示第y次操作后节点i的通电状态,
Figure BDA0002433268680000102
表示节点i在第y次操作后通电;否则节点i没有通电;
Figure BDA0002433268680000103
表示第y次操作后分布式电源i的连接状态,
Figure BDA0002433268680000104
表示分布式电源i接入电网,否则其没有接入电网;
Figure BDA0002433268680000105
表示第y次操作后储能i的连接状态,
Figure BDA0002433268680000106
表示储能i接入电网,否则其没有接入电网;
Figure BDA0002433268680000107
表示分布式电源与储能系统的有功功率最大值和无功功率最大值,为3×1维向量;
拓扑约束包括:
在本实施例中,辐射状拓扑可由图论中的树表示;
通过
Figure BDA0002433268680000108
表示在图论的树基本定义,即不存在环,每一个可以形成环的线路组合中,至少有一条线路是断开的;
通过
Figure BDA0002433268680000109
表示树的连通性,表示树的必要不充分条件,所有通电节点之和减去所有通电线路之和小于等于1;
通过
Figure BDA00024332686800001010
Figure BDA00024332686800001011
Figure BDA00024332686800001012
表示线路通电与节点通电状态和线路连通状态的关系;
公式(15)-(20)中,K为配电网中全部可用线路连接后形成的环的集合;Ck为第k个环中所包含的线路的集合;
Figure BDA00024332686800001013
为第k个环中所包含的线路数量,
Figure BDA00024332686800001014
表示第y次操作后节点线路(i,j)是否通电,
Figure BDA00024332686800001015
表示线路(i,j)通电,否则表示线路(i,j)未通电;
过程操作约束包括:
通过
Figure BDA00024332686800001016
表示只有节点通电时,该节点的负荷才能够接入配电网;
通过
Figure BDA00024332686800001017
表示只有节点通电时,ES才能启动;
通过
Figure BDA00024332686800001018
表示除根节点之外的DG,只有在所连接节点通电时,才能启动;
通过
Figure BDA00024332686800001019
表示当具有黑启动能力且容量最大的DG接入配电网时,其所连接的节点也同时通电;
通过
Figure BDA00024332686800001020
表示节点一旦通电,就不能断电;
通过
Figure BDA0002433268680000111
Figure BDA0002433268680000112
表示DG、ES、负荷一旦接入系统就不能再断开;
通过
Figure BDA0002433268680000113
表示每次只进行一项操作,同时,由于配电网发生故障的位置不同与最终恢复所形成的拓扑变化,在本实施例中将线路分为2种类型,式中,Lc,Lo分别为2种线路类型的集合,发生故障后,初始状态为闭合的线路集合与初始状态为断开的线路集合;
通过
Figure BDA0002433268680000114
表示对于初始状态为断开的线路,上一个操作后的线路连接状态总是闭合或与下一个操作后的线路连接状态相同;
通过
Figure BDA0002433268680000115
表示对于初始状态为闭合的线路,上一个操作后的线路连接状态总是断开或与下一个操作后的线路连接状态相同;
通过
Figure BDA0002433268680000116
表示对于当线路状态为连接时,线路两端的通电节点状态总是相同的;
通过
Figure BDA0002433268680000117
表示对于初始状态为断开的线路,当线路状态没有改变的时候,无约束,当线路状态改变时,在改变之前必然是两端节点有且只有一端是通电的;
通过
Figure BDA0002433268680000118
表示对于初始状态为闭合的线路,当线路状态没有改变的时候,无约束,当线路状态改变时,在改变之前必然是两端节点都没有通电的;
初始状态约束包括:
通过
Figure BDA0002433268680000119
Figure BDA00024332686800001110
表示将初始线路的状态设置为闭合或断开;
通过
Figure BDA00024332686800001111
表示在所有节点在初始状态下均未通电;
通过
Figure BDA00024332686800001112
表示当容量最大的黑启动DG启动时,其所连接的节点电压为恒定电压。
通过加入约束条件,建立配电网故障恢复次序优化模型;通过优化软件进行求解目标函数,获得相应的配电网故障恢复次序优化结果。
具体应用如在本发明提供的一个优选实施例中,如图2所示的测试系统中,一共有12个节点,包括三个分布式电源和一个储能系统,负荷分为三个等级,一级重要负荷的权重系数为100,二级重要负荷的权重系数为10,普通负荷的权重系数为0.2。测试场景为极端事件后,电力基础设施被严重毁坏,输电网送电通路故障,整个配电区域无法从变电站获得电能,线路650-632与线路684-652发生故障并已被隔离,线路654-632、632-633、632-671、611-684、671-684因保护动作而断开连接。
根据图2所示测试系统的信息和场景信息,建立智能配电网故障恢复次序优化模型。通过优化软件进行求解。求解结果如图3、表1所示。目标函数的最优值为1750,2个一级负荷和1个二级负荷被恢复,恢复操作步骤如表1。
Figure BDA0002433268680000121
表1恢复操作序列求解结果
综上所述,本发明提供的一种多源协同的智能配电网故障恢复次序优化决策方法,包括根据配电网的节点分布情况,建立配电网故障恢复的目标函数和相应的优化变量;基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果。本发明提供的方法没通过计及故障恢复过程中系统三相不对称潮流约束以及其他各项稳态运行约束、辐射状拓扑约束和操作过程相关约束,提出更为完善的配电网故障恢复次序优化混合整数线性规划模型,形成了一种智能配电网故障恢复次序优化决策方法,以确保每次恢复操作后均能满足配电网各项稳态约束,确保系统的稳定运行。
本领域普通技术人员可以理解:附图只是一个实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置或系统实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的装置及系统实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (4)

1.一种多源协同的智能配电网故障恢复次序优化决策方法,其特征在于,包括:
建立配电网故障恢复的目标函数,该目标函数具有优化变量;
基于目标函数,通过加入三相不对称潮流约束、稳态运行约束和过程操作约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果。
2.根据权利要求1所述的方法,其特征在于,所述目标函数为:
Figure FDA0002433268670000011
所述优化变量为:公式(1)中,N为配电网中所有节点集合;Y为整个故障恢复过程中所有操作的集合;y表示某一次操作,第0次操作表示配电网故障隔离后恢复的初始状态;y的最大值为ymax=nload+nL+ngen+nes,其中nload表示带有开关的负荷数量,nL表示带开关的线路数量,ngen表示分布式电源的数量,nes表示储能系统的数量;ωi为节点i所连接的负荷的权重系数;
Figure FDA0002433268670000012
为整数变量,表示第y次操作后,节点i所连接的负荷是否被恢复;
Figure FDA0002433268670000013
表示负荷被恢复,否则负荷未被恢复。
3.根据权利要求2所述的方法,其特征在于,所述的基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果的步骤中,三相不对称潮流约束包括:
通过
Figure FDA0002433268670000014
表示流入节点i的功率之和等于流出节点i的总和;
通过
Figure FDA0002433268670000015
表示
Figure FDA0002433268670000016
Figure FDA0002433268670000017
之间的关系;
通过
Figure FDA0002433268670000018
表示线路(i,j)上的两端节点的电压差;
通过
Figure FDA0002433268670000021
表示线路(i,j)是否连接的两种情况,当
Figure FDA0002433268670000022
时,该式为支路首末两端电压模值平方差的关系,当
Figure FDA0002433268670000023
时,支路两端电压无约束;
通过
Figure FDA0002433268670000024
Figure FDA0002433268670000025
Figure FDA0002433268670000026
表明节点i的注入功率;若节点i与分布式电源连接,则节点i的注入功率等于连接到节点i的分布式电源的发电量减去节点i的负荷功率;若节点i与储能连接,则节点i的注入功率等于连接到节点i的储能系统的发电量减去节点i的负荷功率;若节点i不与分布式电源或储能系统连接,则节点i的注入功率等于负的节点i的负荷功率;
上述公式中,G为配电网中所有可控分布式电源的集合;E为配电网中所有储能系统集合;L为配电网中所有支路集合;函数DIAG(·)返回以输入向量为对角元素的对角矩阵;
Figure FDA0002433268670000027
表示第y次操作后从节点h流向i的复功率,为3×3维矩阵;
Figure FDA0002433268670000028
表示
Figure FDA0002433268670000029
的对角线元素,为3×1维向量;δ为线性化转换矩阵;
Figure FDA00024332686700000210
为中间辅助变量;
Figure FDA00024332686700000211
是节点i第y次操作后的电压,为3×1维向量;
Figure FDA00024332686700000212
为i节点第y次操作后的电压向量与其共轭转置相乘得到的3×3维矩阵,
Figure FDA00024332686700000213
是一个对角矩阵,其对角线元素为节点电压各相幅值的平方;zij表示线路(i,j)的3×3维阻抗矩阵;
Figure FDA00024332686700000214
表示第y次操作后节点i的注入功率,为3×1维向量,
Figure FDA00024332686700000215
表示节点i的负荷功率需求,为3×1维向量;
Figure FDA00024332686700000216
表示第y次操作后节点i分布式电源出力,为3×1维向量;
Figure FDA00024332686700000217
表示第y次操作后节点i储能出力,为3×1维向量;
Figure FDA00024332686700000218
表示第y次操作后线路状态,
Figure FDA00024332686700000219
表示线路被连接;否则线路未被连接;M为值很大的正实数。
4.根据权利要求2所述的方法,其特征在于,所述的基于目标函数,通过加入三相不对称潮流约束和稳态运行约束,建立配电网故障恢复次序优化模型,获得配电网故障恢复次序优化结果的步骤中,稳态运行约束包括运行安全约束和拓扑约束;过程操作约束包括初始状态约束;
运行安全约束包括:
通过
Figure FDA0002433268670000031
表示线路功率约束,具体表示线路(i,j)是否连接的两种情况,当
Figure FDA0002433268670000032
时,未对
Figure FDA0002433268670000033
进行限制,当
Figure FDA0002433268670000034
时,限制
Figure FDA0002433268670000035
使得其值为0;
通过
Figure FDA0002433268670000036
表示节点电压上下限约束,若节点没有通电,则限制其电压为0;
通过
Figure FDA0002433268670000037
Figure FDA0002433268670000038
表示节点电压上下限约束,若节点没有通电,则限制其电压为0;
通过
Figure FDA0002433268670000039
Figure FDA00024332686700000310
表示储能的容量约束;
公式(9)-(14)中,函数diag(·)返回由输入矩阵的对角元素组成的向量;函数Re(·)返回输入矩阵的实部;函数Im(·)返回输入矩阵的虚部;函数ones(·)返回维度是输入数的单位矩阵;
Figure FDA00024332686700000311
Figure FDA00024332686700000312
表示节点电压各相幅值最小值的平方与幅值最大值的平方,为3×1维向量;
Figure FDA00024332686700000313
表示第y次操作后节点i的通电状态,
Figure FDA00024332686700000314
表示节点i在第y次操作后通电;否则节点i没有通电;
Figure FDA00024332686700000315
表示第y次操作后分布式电源i的连接状态,
Figure FDA00024332686700000316
表示分布式电源i接入电网,否则其没有接入电网;
Figure FDA00024332686700000317
表示第y次操作后储能i的连接状态,
Figure FDA00024332686700000318
表示储能i接入电网,否则其没有接入电网;
Figure FDA00024332686700000319
表示分布式电源与储能系统的有功功率最大值和无功功率最大值,为3×1维向量;
拓扑约束包括:
通过
Figure FDA00024332686700000320
表示在图论的树的定义中,每一个可以形成环的线路组合中,至少有一条线路是断开的;
通过
Figure FDA00024332686700000321
表示树的连通性,表示树的必要不充分条件,所有通电节点之和减去所有通电线路之和小于等于1;
通过
Figure FDA0002433268670000041
Figure FDA0002433268670000042
Figure FDA0002433268670000043
表示线路通电与节点通电状态和线路连通状态的关系;
公式(15)-(20)中,K为配电网中全部可用线路连接后形成的环的集合;Ck为第k个环中所包含的线路的集合;
Figure FDA0002433268670000044
为第k个环中所包含的线路数量,
Figure FDA0002433268670000045
表示第y次操作后节点线路(i,j)是否通电,
Figure FDA0002433268670000046
表示线路(i,j)通电,否则表示线路(i,j)未通电;
过程操作约束包括:
通过
Figure FDA0002433268670000047
表示只有节点通电时,该节点的负荷才能够接入配电网;
通过
Figure FDA0002433268670000048
表示只有节点通电时,ES才能启动;
通过
Figure FDA0002433268670000049
表示除根节点之外的DG,只有在所连接节点通电时,才能启动;
通过
Figure FDA00024332686700000410
表示当具有黑启动能力且容量最大的DG接入配电网时,其所连接的节点也同时通电;
通过
Figure FDA00024332686700000411
表示节点若通电,就不能断电;
通过
Figure FDA00024332686700000412
Figure FDA00024332686700000413
表示DG、ES、负荷若接入系统就不能再断开;
通过
Figure FDA00024332686700000414
表示每次只进行一项操作,Lc,Lo分别为2种线路类型的集合,发生故障后,初始状态为闭合的线路集合与初始状态为断开的线路集合;
通过
Figure FDA00024332686700000415
表示对于初始状态为断开的线路,上一个操作后的线路连接状态总是闭合或与下一个操作后的线路连接状态相同;
通过
Figure FDA00024332686700000416
表示对于初始状态为闭合的线路,上一个操作后的线路连接状态总是断开或与下一个操作后的线路连接状态相同;
通过
Figure FDA0002433268670000051
表示对于当线路状态为连接时,线路两端的通电节点状态总是相同的;
通过
Figure FDA0002433268670000052
表示对于初始状态为断开的线路,当线路状态没有改变的时候,无约束,当线路状态改变时,在改变之前必然是两端节点有且只有一端是通电的;
通过
Figure FDA0002433268670000053
表示对于初始状态为闭合的线路,当线路状态没有改变的时候,无约束,当线路状态改变时,在改变之前必然是两端节点都没有通电的;
初始状态约束包括:
通过
Figure FDA0002433268670000054
Figure FDA0002433268670000055
表示将初始线路的状态设置为闭合或断开;
通过
Figure FDA0002433268670000056
表示在所有节点在初始状态下均未通电;
通过
Figure FDA0002433268670000057
表示当容量最大的黑启动DG启动时,其所连接的节点电压为恒定电压。
CN202010243290.8A 2020-03-31 2020-03-31 一种多源协同的智能配电网故障恢复次序优化决策方法 Active CN111313416B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010243290.8A CN111313416B (zh) 2020-03-31 2020-03-31 一种多源协同的智能配电网故障恢复次序优化决策方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010243290.8A CN111313416B (zh) 2020-03-31 2020-03-31 一种多源协同的智能配电网故障恢复次序优化决策方法

Publications (2)

Publication Number Publication Date
CN111313416A true CN111313416A (zh) 2020-06-19
CN111313416B CN111313416B (zh) 2022-04-15

Family

ID=71148248

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010243290.8A Active CN111313416B (zh) 2020-03-31 2020-03-31 一种多源协同的智能配电网故障恢复次序优化决策方法

Country Status (1)

Country Link
CN (1) CN111313416B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113868589A (zh) * 2021-09-29 2021-12-31 山东大学 考虑应急资源调度的弹性配电网多点故障抢修方法及系统
CN114050573A (zh) * 2021-11-29 2022-02-15 杭州电力设备制造有限公司 一种有源配电网的故障恢复控制方法、装置、设备及介质
CN114156847A (zh) * 2021-12-08 2022-03-08 国网辽宁省电力有限公司朝阳供电公司 一种配电网接地故障混合优化决策处理方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103579990A (zh) * 2013-10-17 2014-02-12 华南理工大学 基于混合决策图的含分布式电源配电系统故障恢复方法
US9088160B2 (en) * 2007-12-12 2015-07-21 Abb Research Ltd. Load restoration for feeder automation in electric power distribution systems

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9088160B2 (en) * 2007-12-12 2015-07-21 Abb Research Ltd. Load restoration for feeder automation in electric power distribution systems
CN103579990A (zh) * 2013-10-17 2014-02-12 华南理工大学 基于混合决策图的含分布式电源配电系统故障恢复方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JIAXU LI 等: "Critical load restoration method for unbalanced distribution systems based on linear approximation", 《2018 IEEE 2ND INTERNATIONAL ELECTRICAL AND ENERGY CONFERENCE》 *
徐玉琴 等: "基于多智能体遗传算法并考虑分布式电源的配电网大面积断电供电恢复算法", 《电工技术学报》 *
陈帅: "大停电背景下配电网孤岛划分与多时段故障恢复的鲁棒优化", 《中国优秀硕士学位论文全文数据库(电子期刊)工程科技Ⅱ辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113868589A (zh) * 2021-09-29 2021-12-31 山东大学 考虑应急资源调度的弹性配电网多点故障抢修方法及系统
CN113868589B (zh) * 2021-09-29 2024-05-17 山东大学 考虑应急资源调度的弹性配电网多点故障抢修方法及系统
CN114050573A (zh) * 2021-11-29 2022-02-15 杭州电力设备制造有限公司 一种有源配电网的故障恢复控制方法、装置、设备及介质
CN114156847A (zh) * 2021-12-08 2022-03-08 国网辽宁省电力有限公司朝阳供电公司 一种配电网接地故障混合优化决策处理方法

Also Published As

Publication number Publication date
CN111313416B (zh) 2022-04-15

Similar Documents

Publication Publication Date Title
CN111313416B (zh) 一种多源协同的智能配电网故障恢复次序优化决策方法
CN112271727B (zh) 含柔性软开关的韧性配电网的故障恢复方法
CN108182498A (zh) 配电网故障恢复性重构方法
CN108683173A (zh) 船舶中压直流配电网络故障工况粒子群重构方法
Biswas et al. Optimal power flow solutions using algorithm success history based adaptive differential evolution with linear population reduction
CN108539781B (zh) 提高恢复过程安全性的扩展黑启动方案二层规划优化方法
CN114050573A (zh) 一种有源配电网的故障恢复控制方法、装置、设备及介质
Janssens et al. Loop flows in a ring AC power system
CN108695865B (zh) 一种多源协同的三相不对称配电网故障恢复策略生成方法
CN106487001B (zh) 一种独立电力系统智能重构方法
Sheng et al. Collaborative Fault Recovery and Network Reconstruction Method for Cyber-physical-systems Based on Double Layer Optimization
CN114784796A (zh) 一种基于多端sop的柔性互联配电系统多阶段恢复方法
Nakka et al. Design and simulation of monarch butterfly optimized SP-UPFC for power quality compensation
CN113595073A (zh) 基于故障信息获取的不平衡配电网负荷恢复方法及系统
Bessani et al. Robustness of reconfigurable complex systems by a multi-agent simulation: Application on power distribution systems
Warudkar et al. Particle swarm optimization based transmission expansion planning
CN110165712A (zh) 一种基于网络流约束推导的骨干网架规划建模方法
Xiao et al. A robust mixed-integer second-order cone programming for service restoration of distribution network
JP5014254B2 (ja) 系統分離判定方法およびそのプログラム
Quang et al. Distribution Network Reconfiguration Using Genetic Algorithm Considering Load Profile and the Penetration of Distributed Generation
Mahdavi et al. Genetic algorithm application in distribution system reconfiguration
Wu et al. Service restoration of active distribution network considering the islanded operation of distributed generation and micro-grid
Tsang et al. ANN controlled battery energy storage system for enhancing power system stability
CN114336749B (zh) 一种配电网优化方法、系统、装置及存储介质
CN109473980B (zh) 一种配电网的重构方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant