CN115549083A - 考虑配电网故障线路修复顺序的故障线路修复方法 - Google Patents

考虑配电网故障线路修复顺序的故障线路修复方法 Download PDF

Info

Publication number
CN115549083A
CN115549083A CN202211324621.6A CN202211324621A CN115549083A CN 115549083 A CN115549083 A CN 115549083A CN 202211324621 A CN202211324621 A CN 202211324621A CN 115549083 A CN115549083 A CN 115549083A
Authority
CN
China
Prior art keywords
stage
line
state
node
power
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.)
Pending
Application number
CN202211324621.6A
Other languages
English (en)
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.)
Yibin Power Supply Co Of Sichuan Electric Power Corp
Xihua University
Original Assignee
Yibin Power Supply Co Of Sichuan Electric Power Corp
Xihua 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 Yibin Power Supply Co Of Sichuan Electric Power Corp, Xihua University filed Critical Yibin Power Supply Co Of Sichuan Electric Power Corp
Priority to CN202211324621.6A priority Critical patent/CN115549083A/zh
Publication of CN115549083A publication Critical patent/CN115549083A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/20—Administration of product repair or maintenance
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
    • H02J3/00125—Transmission line or load transient problems, e.g. overvoltage, resonance or self-excitation of inductive loads
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Power Engineering (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了考虑配电网故障线路修复顺序的故障线路修复方法，涉及配电网故障线路修复，包括以下步骤：S1、将故障线路修复分为T个阶段，并将修复之前配电网的拓扑结构和失负荷量设置为初始状态；S11、利用图论求解每个阶段的拓扑图，得到满足约束条件的每个阶段的状态集合；S2、计算相邻阶段状态集合中任意两个状态之间的失负荷量和发电成本，所述失负荷量和发电成本作为决策量，得到决策量集合；S3、根据决策量集合，计算决策序列；S4、遍历所述故障线路修复的所有阶段，得到最终的最优决策序列；不仅减少了故障线路修复过程中的失负荷量，还减少了故障线路修复过程中的成本，提高了配电网在极端天气下的供电能力。

Description

考虑配电网故障线路修复顺序的故障线路修复方法
技术领域
本发明涉及配电网故障线路修复，具体涉及考虑配电网故障线路修复顺序的故障线路修复方法。
背景技术
随着全球气候变化，极端天气的发生日益频繁，造成的电力网大停电事故频频发生，带来了巨大的经济损失和社会影响。例如严重的冰冻雨雪凝冻灾害天气，给电力系统安全稳定运行和电力供应带来极大的影响，部分地区中低压电网遭受毁灭性破坏，极端灾害下中低压配电线路倒杆、断杆和断线等现象十分严重，因此，为了减小极端天气导致的电网大面积停电的影响，以及由此带来的停电损失，提高电网对极端天气的抵抗能力显得尤为重要。
虽然配电网故障有灾前的预防措施，但在小概率的极端天气下，配电网仅靠事前的防控措施无法完全抵御灾害；目前，有利用分布式电源和微网，运用网络重构和孤岛划分的方法，通过调整线路开关的状态来达到恢复负荷供电的目的，但未涉及到灾后配电网的故障抢修问题；或者，有考虑故障抢修的供电恢复策略，但仅仅是对网络中联络开关和分段开关的开合顺序做出判断，并没有给出具体故障点的恢复顺序。
发明内容
本发明所要解决的技术问题是提高配电网的供电能力，目的在于提供考虑配电网故障线路修复顺序的故障线路修复方法。
本发明通过下述技术方案实现：
第一方面提供考虑配电网故障线路修复顺序的故障线路修复方法，包括以下步骤：
S1、将故障线路修复分为T个阶段；
S11、利用图论求解每个阶段的拓扑图，得到每个阶段的状态集合；
S2、计算相邻阶段状态集合中任意两个状态之间的失负荷量和发电成本，上述失负荷量和发电成本作为决策量，得到决策量集合；
S3、根据决策量集合，计算决策序列；
S4、遍历上述故障线路修复的所有阶段，得到最终的最优决策序列。
将上述故障线路修复分为T个阶段，在各阶段利用图论求解拓扑图，确定每个阶段可能选择的线路修复的网络拓扑结构解集，上述运用图论求解每个阶段的拓扑图，舍去不符合情况的拓扑，减少运算量；再利用动态规划算法选择出失负荷量和发电成本最小的最优决策序列，使得故障线路修复顺序达到最优，不仅减少了故障线路修复过程中的失负荷量，还减少了故障线路修复过程中的成本，提高了配电网在极端天气下的供电能力。
进一步的，上述状态集合中的状态为每个阶段故障修复线路修复完成之后的配电网网络结构。
进一步的，计算相邻阶段状态集合中任意两个状态之间的失负荷量的公式如下：
Figure BDA0003911917220000021
其中，
Figure BDA0003911917220000022
表示第t-1阶段的第j个状态到第t阶段的第q个状态的失负荷量；
Figure BDA0003911917220000023
表示第t阶段的第q个状态是否选择修复线路li，选择为1，不选择为0；ΩF,t-1表示第t-1阶段的故障线路集合；ΩS,t-1表示第t-1阶段的所有失负荷节点集合；Ω0,t表示在第t阶段接入主网的节点集合；ωi表示节点i的负荷权重；
Figure BDA0003911917220000024
表示节点i的负荷量；
Figure BDA0003911917220000025
表示节点i在第t阶段的第q个状态下是否接入主网，接入为1，未接入为0；
Figure BDA0003911917220000026
表示节点i在第t阶段的第q个状态下是否与分布式发电机组连接，连接为1，未连接为0；
Figure BDA0003911917220000027
表示节点i所带负荷量在第t阶段是否恢复供电，恢复为1，未恢复为0。
进一步的，计算相邻阶段状态集合中任意两个状态之间的发电成本的公式如下：
Figure BDA0003911917220000028
其中，
Figure BDA0003911917220000029
表示第t-1阶段的第j个状态到第t阶段的第q个状态的发电成本；
Figure BDA00039119172200000210
表示第t阶段的第q个状态是否选择修复线路li，选择为1，不选择为0；wM表示分布式发电的发电成本；ΩM表示分布式发电所在节点的集合；
Figure BDA00039119172200000211
表示节点i在第t阶段的输出功率。
进一步的，上述失负荷量作为主决策量，发电成本作为辅助决策量，首先考虑配电网故障线路修复的可靠性，再考虑成本问题，减少用户的损失。
进一步的，根据决策量集合，计算决策序列的公式如下：
Figure BDA00039119172200000212
其中，Wt j,q表示第0阶段的初始状态到第t阶段的第q个状态的最优决策序列；
Figure BDA00039119172200000213
表示第0阶段的初始状态到第t-1阶段的第j个状态的最优决策序列；
Figure BDA0003911917220000031
表示第t-1阶段的第j个状态到第t阶段的第q个状态的失负荷量；
Figure BDA0003911917220000032
表示第t-1阶段的第j个状态到第t阶段的第q个状态的发电成本。
利用动态规划算法的求解顺序，从初始阶段的初始状态递推到第T阶段的最终状态，得到最终的最优决策序列，具有全面性，最大程度上减少了失负荷量和发电成本。
进一步的，遍历上述故障修复线路的所有阶段，得到最终的最优决策序列，包括以下步骤：
判断用于计算决策序列的第t阶段是否为第T阶段；
若不是，则返回S11继续执行；
若是，则输出该决策序列。
进一步的，在执行S1之前，基于双层规划模型训练，得到故障线路修复模型。
进一步的，所述故障线路修复模型的上层以故障线路修复的失负荷量最小为目标，下层以发电成本最小为目标。
以上述失负荷量和发电成本最小为目标，得到最优的故障线路修复顺序，既考虑了可靠问题，也考虑了经费问题。
进一步的，上述故障线路修复模型的约束条件包括：
(1)分布式发电约束，拓扑图中各节点的分布式发电功率不高于对应节点的上限分布式发电功率，不低于对应节点的下限分布式发电功率；
(2)配电网线性化潮流约束，通过潮流计算拓扑图中各节点的电压幅值和相角，其中电压幅值的误差不大于5％，相角的误差不大于0.5度；
(3)节点电压约束，拓扑图中各节点的节点电压不高于对应节点的上限节点电压，不低于对应节点的下限节点电压；
(4)线路容量约束，拓扑图中各线路的有功功率和无功功率的平方和不高于对应线路的功率最大值；
(5)抢修资源约束，在故障线路修复中，每个阶段只能修复h个故障点，任一线路只能被修复一次；
(6)拓扑辐射状约束，
Figure BDA0003911917220000033
Figure BDA0003911917220000034
rank(Ek)＝N-1
其中，Sl表示线路l的状态，1表示线路正常，0表示线路故障；N表示节点数量；ΩG表示拓扑图中包含的线路集合；ΩF表示故障线路的集合；Ek表示拓扑图中各节点和各线路的关联矩阵；rank(Ek)＝N-1表示拓扑辐射状约束。
本发明与现有技术相比，具有如下的优点和有益效果：
将上述故障线路修复分为T个阶段，在各阶段利用图论求解拓扑图，确定每个阶段可能选择的线路修复的网络拓扑结构解集，上述运用图论求解每个阶段的拓扑图，舍去不符合情况的拓扑，减少运算量；
再计算相邻阶段状态集合中任意两个状态之间的失负荷量和发电成本，上述失负荷量和发电成本作为决策量，得到决策量集合；根据决策量集合，计算决策序列；利用动态规划算法选择出失负荷量和发电成本最小的最优决策序列，使得故障线路修复顺序达到最优，不仅减少了故障线路修复过程中的失负荷量，还减少了故障线路修复过程中的成本，利用该故障线路修复方法可以快速修复故障线路，提高了配电网在极端天气下的供电能力。
附图说明
为了更清楚地说明本发明示例性实施方式的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。在附图中：
图1为实施例1提供的主流程图；
图2为实施例3提供的配电网络拓扑结构示意图；
图3为实施例4提供的各阶段状态之间的关系示意图；
图4为实施例5提供的拓扑图中故障点位置示意图；
图5为实施例5提供的相邻阶段任意两个状态之间的失负荷量和DG发电成本的部分结果示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1
本实施例1提供考虑配电网故障线路修复顺序的故障线路修复方法，如图1所示，包括以下步骤：
S1、将故障线路修复分为T个阶段，并将修复之前配电网的拓扑结构和失负荷量设置为初始状态；
S11、利用图论求解每个阶段的拓扑图，得到满足约束条件的每个阶段的状态集合；
S2、计算相邻阶段状态集合中任意两个状态之间的失负荷量和发电成本，上述失负荷量和发电成本作为决策量，得到决策量集合；
S3、根据决策量集合，计算决策序列；
S4、遍历上述故障线路修复的所有阶段，得到最终的最优决策序列。
将上述故障线路修复分为T个阶段，在各阶段利用图论求解拓扑图，确定每个阶段可能选择的线路修复的网络拓扑结构解集，上述运用图论求解每个阶段的拓扑图，舍去不符合情况的拓扑，减少运算量；再利用动态规划算法选择出失负荷量和发电成本最小的最优决策序列，使得故障线路修复顺序达到最优，不仅减少了故障线路修复过程中的失负荷量，还减少了故障线路修复过程中的成本，提高了配电网在极端天气下的供电能力。
具体的实施例，上述状态集合中的状态为每个阶段故障修复线路修复完成之后的配电网网络结构。
具体的实施例，计算相邻阶段状态集合中任意两个状态之间的失负荷量的公式如下：
Figure BDA0003911917220000051
其中，
Figure BDA0003911917220000052
表示第t-1阶段的第j个状态到第t阶段的第q个状态的失负荷量；
Figure BDA0003911917220000053
表示第t阶段的第q个状态是否选择修复线路li，选择为1，不选择为0；ΩF,t-1表示第t-1阶段的故障线路集合；ΩS,t-1表示第t-1阶段的所有失负荷节点集合；Ω0,t表示在第t阶段接入主网的节点集合；ωi表示节点i的负荷权重；
Figure BDA0003911917220000054
表示节点i的负荷量；
Figure BDA0003911917220000055
表示节点i在第t阶段的第q个状态下是否接入主网，接入为1，未接入为0；
Figure BDA0003911917220000056
表示节点i在第t阶段的第q个状态下是否与分布式发电线连接，连接为1，未连接为0；
Figure BDA0003911917220000057
表示节点i所带负荷量在第t阶段是否恢复供电，恢复为1，未恢复为0。
具体的实施例，计算相邻阶段状态集合中任意两个状态之间的发电成本的公式如下：
Figure BDA0003911917220000058
其中，
Figure BDA0003911917220000061
表示第t-1阶段的第j个状态到第t阶段的第q个状态的发电成本；
Figure BDA0003911917220000062
表示第t阶段的第q个状态是否选择修复线路li，选择为1，不选择为0；wM表示分布式发电的发电成本；ΩM表示分布式发电所在节点的集合；
Figure BDA0003911917220000063
表示节点i在第t阶段的输出功率。
具体的实施例，上述失负荷量作为主决策量，发电成本作为辅助决策量，首先考虑配电网故障线路修复的可靠性，再考虑成本问题，减少用户的损失。
具体的实施例，根据决策量集合，计算决策序列的公式如下：
Figure BDA0003911917220000064
其中，Wt j,q表示第0阶段的初始状态到第t阶段的第q个状态的最优决策序列；
Figure BDA0003911917220000065
表示第0阶段的初始状态到第t-1阶段的第j个状态的最优决策序列；
Figure BDA0003911917220000066
表示第t-1阶段的第j个状态到第t阶段的第q个状态的失负荷量；
Figure BDA0003911917220000067
表示第t-1阶段的第j个状态到第t阶段的第q个状态的发电成本。
利用动态规划算法的求解顺序，从初始阶段的初始状态递推到第T阶段的最终状态，得到最终的最优决策序列，具有全面性，最大程度上减少了失负荷量和发电成本。
具体的实施例，遍历上述故障修复线路的所有阶段，得到最终的最优决策序列，包括以下步骤：
判断用于计算决策序列的第t阶段是否为第T阶段；
若不是，则返回S11继续执行；
若是，则输出该决策序列。
具体的实施例，在执行S1之前，基于双层规划模型训练，得到故障线路修复模型。
实施例2
上述故障线路修复的阶段划分是根据故障点数量和运行人员数量定义的。
具体的实施例，上述状态集合为第t阶段的拓扑图解集{Mt,1,Mt,2,…,Mt,i,…,Mt,n}，其中Mt,i表示第t阶段的第i个状态。
具体的实施例，当状态Mt-1,j向第t阶段转移的决策ut-1(j)以及状态Mt-1,j确定之后，第t阶段的状态就被确定了，即状态转移方程如下：
Mt,q＝Tt-1(Mt-1,j,ut-1(j))
其中，Mt,q表示第t阶段的第q个状态，Tt-1表示故障线路修复的第t-1阶段，Mt-1,j表示第t-1阶段的第j个状态；上述决策ut-1(j)包括选择修复的线路和需要操作的联络开关。
具体的实施例，上述决策序列从初始状态到任意一个阶段的任意一个状态的所有决策构成的序列为S1,t＝{u1,u2,...,ut}。
具体的实施例，上述决策量用于判断一个决策优劣的评价指标，状态Mt-1,j和Mt,q之间的决策为ut-1(j)＝q，在状态Mt,q下，用状态Mt-1,j和Mt,q之间的失负荷量作为主决策量，微型燃气轮机和光伏发电成本作为辅助决策量；从初始阶段的初始状态递推到第T阶段的最终状态，得到最优决策序列。
实施例3
本实施例3提供故障线路修复模型，上述的故障线路修复方法用于求解该故障线路修复模型。
具体的实施例，上述故障线路修复模型是基于双层规划模型训练得到的；
上述故障线路修复模型的上层以故障线路修复的失负荷量最小为目标，下层以发电成本最小为目标。
目标函数如下：
Figure BDA0003911917220000071
其中，obj表示目标函数；f表示失负荷量；CD表示发电成本；T表示故障线路修复的总阶段数；
Figure BDA0003911917220000072
表示第t阶段微型燃气轮机和光伏的发电成本；
Figure BDA0003911917220000073
表示第t阶段的失负荷量；Ωu表示上层的故障线路修复顺序集合；Ωd表示下层的故障线路修复顺序集合。
以上述失负荷量和发电成本最小为目标，得到最优的故障线路修复顺序，既考虑了可靠问题，也考虑了经费问题。
具体的实施例，上述故障线路修复模型的约束条件包括：
(1)分布式发电约束，拓扑图中各节点的分布式发电功率不高于对应节点的上限分布式发电功率，不低于对应节点的下限分布式发电功率，公式如下：
Figure BDA0003911917220000081
其中，
Figure BDA0003911917220000082
表示节点i的分布式发电功率；
Figure BDA0003911917220000083
表示节点i的下限分布式发电功率；
Figure BDA0003911917220000084
表示节点i的上限分布式发电功率。
(2)配电网线性化潮流约束，通过潮流计算拓扑图中各节点的电压幅值和相角，其中电压幅值的误差不大于5％，相角的误差不大于0.5度，公式如下：
Figure BDA0003911917220000085
Figure BDA0003911917220000086
Figure BDA0003911917220000087
Figure BDA0003911917220000088
Figure BDA0003911917220000089
Figure BDA00039119172200000810
其中，Pij,t和Qij,t分别表示第t阶段线路lij的有功功率和无功功率；δi,t表示第t阶段节点i的电压相角；Gij和Bij分别表示线路lij导纳的实部和虚部；Vi,t和Vj,t分别表示第t阶段节点i和节点j的电压；
Figure BDA00039119172200000811
和
Figure BDA00039119172200000812
分别表示第t阶段节点i输出的有功功率和无功功率；
Figure BDA00039119172200000813
表示节点i在第t阶段分布式发电的输出功率；
Figure BDA00039119172200000814
和
Figure BDA00039119172200000815
分别表示节点i在第t阶段的有功负荷和无功负荷。
(3)节点电压约束，拓扑图中各节点的节点电压不高于对应节点的上限节点电压，不低于对应节点的下限节点电压，公式如下：
Vmin≤Vi,t≤Vmax
其中，Vi,t表示节点i的节点电压；Vmin表示节点i的下限节点电压；Vmax表示节点i的上限节点电压。
(4)线路容量约束，拓扑图中各线路的有功功率和无功功率的平方和不高于对应线路的功率最大值，公式如下：
Figure BDA0003911917220000091
其中，Sij,tmax表示线路lij的最大功率。
(5)抢修资源约束，公式如下：
Figure BDA0003911917220000092
表示在整个故障线路修复过程中，任一线路li只能被选择修复一次；
Figure BDA0003911917220000093
表示在任意一个第t阶段只能修复h个故障点。
(6)拓扑辐射状约束，根据图论原理，将Ωp定义为图中所有节点的集合，ΩL定义为图中所有线路的集合，在有向图Gk(Ωp,ΩL)中,把没有环路的图定义为树，利用深度优先搜索算法对所有的节点和线路进行编号，求解出有向图Gk中节点和线路的关联矩阵Ek，Ek中的元素包括1、-1和0，Ek＝[e1,e2,...,eB]。如果Ek是n个节点的有向图，则矩阵Ek的秩为n-1是保证有向图为辐射状的充要条件，若一个图能被分成多个有向图，则这些有向图的集合为{G1,G2,...,GNG}，NG表示有向图的数量。
在故障抢修过程中，每个阶段的故障线路修复完成之后，需要对配电网中的联络开关进行开合操作，让尽可能多的失负荷节点与主网或者DG相连，以保证恢复的供电量最大。
如图2所示，在发生故障之后有多种联络开关的操作方式可以形成新的网络拓扑集合，维持网络辐射状的充要条件如下：
Figure BDA0003911917220000094
Figure BDA0003911917220000095
rank(Ek)＝N-1
其中，Sl表示线路l的状态，1表示线路正常，0表示线路故障；N表示节点数量；ΩG表示拓扑图中包含的线路集合；ΩF表示故障线路的集合；Ek表示拓扑图中各节点和各线路的关联矩阵；rank(Ek)表示拓扑辐射状约束；
Figure BDA0003911917220000101
表示生成的拓扑图是连通的。
实施例4
另一种配电网故障修复方法包括以下步骤：
步骤A：将故障线路修复之前配电网的拓扑图和失负荷量作为初始状态；
步骤B：利用图论求解拓扑图，得到每个阶段满足故障线路修复模型的约束条件的状态集合，如图3所示，为每个阶段状态之间的关系示意图；
步骤C：计算出第t-1阶段到第t阶段的任意两个状态之间的失负荷量和发电成本，得到决策量集合；
步骤D：根据公式Mt,q＝Tt-1(Mt-1,j,ut-1(j))，得到从初始状态到第t-1阶段的决策序列；
步骤E：判断第t阶段是否是最终阶段(即第T阶段)，若不是，则返回步骤B；若是，则对最终阶段的决策量集合中的决策量进行选择，得到最终的最优决策序列。
实施例5
以包括3个联络开关(TS)、7个光伏(PV)和4个微型燃气轮机(MT)的改进PG&69节点配电网系统进行案例分析；各系统基准电压为12.66KV，基准功率为10MVA，PV和MT的额定功率都为200kW，PV的出力成本为0元，MT的出力成本为0.4元/kW。灾后网络中有11条故障线路，各个故障点的位置如图4所示。
假设本文的故障恢复分为4个阶段，前三个阶段每阶段恢复3条故障线路，第四阶段恢复剩下的2条线路，网络中各节点的负荷等级和负荷权重如表1所示。
表1
Figure BDA0003911917220000102
运用组合的方法可以列举出每个阶段所有线路修复方式和对应的联络开关的操作方式，再根据图论的方法，结合约束条件筛选得到每个阶段符合实际情况的网络拓扑状态集合，确定拓扑状态之后求出相邻阶段任意两个状态之间的失负荷量和DG发电成本，部分结果如图5所示，可知，同一阶段的两个状态之间可能会计算得到相同的失负荷量，可以选择具有较低的发电成本的线路进行修复；比如，第3阶段的5,7,9-3,8,11和1,3,4-2,9,10，这两种方式的失负荷量都是71.7kW，可以根据MT的发电成本，选择发电成本相对来说更少为84.15元的修复方式，计算得到了所有决策序列和决策量之后，再从最终的结果往前推，得到本算例的最优恢复顺序如表2所示。
表2
阶段 修复线路 TS操作 失负荷量/KW MT成本/元
1 1,2,6 TS1 586 365.7
2 3,4,9 TS1、TS2 0 0
3 5,8,10 / 0 0
4 7,11 / 0 0
通过比对表2和图5可以看出，按照最优顺序抢修线路，经过1个阶段的抢修之后，所有的失负荷点全部恢复供电，而图5中反映出的其他方式在第三阶段抢修完成之后还存在失负荷情况。表2中，在第3阶段抢修完之后MT发电成本为0，说明此时所有负荷都由主网或者PV供电，相对于图5中的其他方法来说具有更好的经济性。从表2可见，第2阶段恢复之后不再进行任何联络开关的操作，减少了开关的开断成本，降低了修复操作的复杂性。
为了更进一步体现出本发明方法的合理性，设置了以下3种恢复策略进行对比说明。
策略1：按照单阶段恢复最优求解恢复顺序；
策略2：不考虑联络线在抢修过程中的作用。
策略3：建立的目标函数只考虑失负荷最少，不考虑发电成本。
3个策略的结果如表3所示，策略1在第一阶段选择恢复线路1,2,3，失负荷量会比本文的方法少，但是策略1的第2阶段线路修复之后仍然存在失负荷点，而且总的失负荷量比用本文的方法多，MT的发电成本也有所增加。因此策略1在可靠性和经济性两个方面都不如本文的方法好。相对于策略1和本文的方法，策略2各阶段的失负荷量和发电成本都大大增加，说明有效利用联络线调整网络的拓扑结构，有利于负荷恢复供电。策略3的整个阶段失负荷量和本文方法的失负荷一样，但是策略3中MT出力成本高于本文所提方法。因此用本专利的方法不仅失负荷最小也能提高经济性，更加符合实际情况。
表3
Figure BDA0003911917220000111
Figure BDA0003911917220000121
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，包括以下步骤：
S1、将故障线路修复分为T个阶段；
S11、利用图论求解每个阶段的拓扑图，得到每个阶段的状态集合；
S2、计算相邻阶段状态集合中任意两个状态之间的失负荷量和发电成本，所述失负荷量和发电成本作为决策量，得到决策量集合；
S3、根据决策量集合，计算决策序列；
S4、遍历所述故障线路修复的所有阶段，得到最终的最优决策序列。
2.根据权利要求1所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，所述状态集合中的状态为每个阶段故障修复线路修复完成之后的配电网网络结构。
3.根据权利要求1所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，计算相邻阶段状态集合中任意两个状态之间的失负荷量的公式如下：
Figure FDA0003911917210000011
其中，
Figure FDA0003911917210000012
表示第t-1阶段的第j个状态到第t阶段的第q个状态的失负荷量；
Figure FDA0003911917210000013
表示第t阶段的第q个状态是否选择修复线路li，选择为1，不选择为0；ΩF,t-1表示第t-1阶段的故障线路集合；ΩS,t-1表示第t-1阶段的所有失负荷节点集合；Ω0,t表示在第t阶段接入主网的节点集合；ωi表示节点i的负荷权重；
Figure FDA0003911917210000014
表示节点i的负荷量；
Figure FDA0003911917210000015
表示节点i在第t阶段的第q个状态下是否接入主网，接入为1，未接入为0；
Figure FDA0003911917210000016
表示节点i在第t阶段的第q个状态下是否与分布式发电机组连接，连接为1，未连接为0；
Figure FDA0003911917210000017
表示节点i所带负荷量在第t阶段是否恢复供电，恢复为1，未恢复为0。
4.根据权利要求3所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，计算相邻阶段状态集合中任意两个状态之间的发电成本的公式如下：
Figure FDA0003911917210000018
其中，
Figure FDA0003911917210000019
表示第t-1阶段的第j个状态到第t阶段的第q个状态的发电成本；
Figure FDA00039119172100000110
表示第t阶段的第q个状态是否选择修复线路li，选择为1，不选择为0；wM表示分布式发电的发电成本；ΩM表示分布式发电所在节点的集合；
Figure FDA0003911917210000021
表示节点i在第t阶段的输出功率。
5.根据权利要求1所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，所述失负荷量作为主决策量，发电成本作为辅助决策量。
6.根据权利要求4所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，根据决策量集合，计算决策序列的公式如下：
Figure FDA0003911917210000022
其中，Wt j,q表示第0阶段的初始状态到第t阶段的第q个状态的最优决策序列；
Figure FDA0003911917210000023
表示第0阶段的初始状态到第t-1阶段的第j个状态的最优决策序列；
Figure FDA0003911917210000024
表示第t-1阶段的第j个状态到第t阶段的第q个状态的失负荷量；
Figure FDA0003911917210000025
表示第t-1阶段的第j个状态到第t阶段的第q个状态的发电成本。
7.根据权利要求1所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，遍历所述故障修复线路的所有阶段，得到最终的最优决策序列，包括以下步骤：
判断用于计算决策序列的第t阶段是否为第T阶段；
若不是，则返回S11继续执行；
若是，则输出该决策序列。
8.根据权利要求1所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，在执行S1之前，基于双层规划模型训练，得到故障线路修复模型。
9.根据权利要求8所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，所述故障线路修复模型的上层以故障线路修复的失负荷量最小为目标，下层以发电成本最小为目标。
10.根据权利要求9所述的考虑配电网故障线路修复顺序的故障线路修复方法，其特征在于，所述故障线路修复模型的约束条件包括：
(1)分布式发电约束，拓扑图中各节点的分布式发电功率不高于对应节点的上限分布式发电功率，不低于对应节点的下限分布式发电功率；
(2)配电网线性化潮流约束，通过潮流计算拓扑图中各节点的电压幅值和相角，其中电压幅值的误差不大于5％，相角的误差不大于0.5度；
(3)节点电压约束，拓扑图中各节点的节点电压不高于对应节点的上限节点电压，不低于对应节点的下限节点电压；
(4)线路容量约束，拓扑图中各线路的有功功率和无功功率的平方和不高于对应线路的功率最大值；
(5)抢修资源约束，在故障线路修复中，每个阶段修复h个故障点，任一线路只能被修复一次；
(6)拓扑辐射状约束，
Figure FDA0003911917210000031
Figure FDA0003911917210000032
rank(Ek)＝N-1
其中，Sl表示线路l的状态，1表示线路正常，0表示线路故障；N表示节点数量；ΩG表示拓扑图中包含的线路集合；ΩF表示故障线路的集合；Ek表示拓扑图中各节点和各线路的关联矩阵；rank(Ek)＝N-1表示拓扑辐射状约束。
CN202211324621.6A 2022-10-27 2022-10-27 考虑配电网故障线路修复顺序的故障线路修复方法 Pending CN115549083A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211324621.6A CN115549083A (zh) 2022-10-27 2022-10-27 考虑配电网故障线路修复顺序的故障线路修复方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211324621.6A CN115549083A (zh) 2022-10-27 2022-10-27 考虑配电网故障线路修复顺序的故障线路修复方法

Publications (1)

Publication Number Publication Date
CN115549083A true CN115549083A (zh) 2022-12-30

Family

ID=84719602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211324621.6A Pending CN115549083A (zh) 2022-10-27 2022-10-27 考虑配电网故障线路修复顺序的故障线路修复方法

Country Status (1)

Country Link
CN (1) CN115549083A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115935711A (zh) * 2023-01-31 2023-04-07 西华大学 一种基于图论的多能互补主动配网系统可靠性评估方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115935711A (zh) * 2023-01-31 2023-04-07 西华大学 一种基于图论的多能互补主动配网系统可靠性评估方法
CN115935711B (zh) * 2023-01-31 2023-10-20 西华大学 一种基于图论的多能互补主动配网系统可靠性评估方法

Similar Documents

Publication Publication Date Title
CN112217202B (zh) 考虑灵活性投资的分布式新能源、储能与配电网规划方法
CN111478358B (zh) 一种计及分布式电源不确定性的配电网鲁棒恢复决策方法
CN113489033A (zh) 面向韧性提升的配电网分布式电源储能控制方法及系统
CN101685968B (zh) 配电网可靠性评估的故障扩散方法
CN104820864B (zh) 含分布式电源的智能配电网全景故障恢复博弈方法
Wang et al. PTDF-based automatic restoration path selection
CN112909942B (zh) 基于混合求解策略的主动配电网故障恢复方法及系统
CN113239512B (zh) 一种考虑韧性的交直流配电网规划方案筛选方法及系统
CN110676839A (zh) 一种交直流混合配电网可靠性评估方法及系统
CN116882139A (zh) 一种考虑网络重构与故障抢修滚动优化的主动配电网元件修复方法及系统
CN114050573B (zh) 一种有源配电网的故障恢复控制方法、装置、设备及介质
Wang et al. Resilience enhancement strategy using microgrids in distribution network
CN115549083A (zh) 考虑配电网故障线路修复顺序的故障线路修复方法
CN114865625A (zh) 一种包含微网的配电网故障恢复方法
CN112769124B (zh) 一种基于潮流转移和追踪的电力系统快速运行风险评估方法
Jiang et al. Power system emergency control to improve short-term voltage stability using deep reinforcement learning algorithm
CN107179688B (zh) 考虑蒙特卡罗状态抽样截尾的电力系统可靠性分析方法
CN111126841A (zh) 一种基于随机潮流的配电网负荷转供能力评价方法
CN116169673A (zh) 配电网的弹性评估方法、系统、电子设备及可读存储介质
CN115800397A (zh) 一种配电台区拓扑优化方法、系统、设备及介质
CN114676913A (zh) 一种结合道路网与配电网的电力应急资源优化方法及系统
CN111951126B (zh) 一种含柔性多状态开关的有源配电网多阶段供电恢复方法
CN111106622B (zh) 一种基于rmpc的主动配电网供电恢复方法
CN112766532A (zh) 一种基于改进混合整数差分进化算法的dg规划方法
CN115396294B (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