CN103001220A - 一种基于最优匹配回路功率的配电网拓扑结构辨识方法 - Google Patents
一种基于最优匹配回路功率的配电网拓扑结构辨识方法 Download PDFInfo
- Publication number
- CN103001220A CN103001220A CN2012105714901A CN201210571490A CN103001220A CN 103001220 A CN103001220 A CN 103001220A CN 2012105714901 A CN2012105714901 A CN 2012105714901A CN 201210571490 A CN201210571490 A CN 201210571490A CN 103001220 A CN103001220 A CN 103001220A
- Authority
- CN
- China
- Prior art keywords
- power
- opt
- distribution network
- delta
- power distribution
- 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
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种基于最优匹配回路功率的配电网拓扑结构辨识方法,属于电力***调度自动化与电网仿真技术领域。将电量计费***采集的负荷数据作为负荷量测;并利用配电网中的功率量测计算配电网各独立回路的最优匹配回路功率;然后根据最优匹配回路功率得到若干种可能的配电网拓扑结构;分别对可能的拓扑结构进行状态估计;最后将状态估计结果与量测最匹配的拓扑结构作为配电网正确的拓扑结构。本发明提出的方法充分利用了配电网中的支路功率量测,以及配电网计量***中的负荷数据,且方法实现简单,计算速度快,结果准确。
Description
技术领域
本发明涉及一种基于最优匹配回路功率的配电网拓扑结构辨识方法,属于电力***调度自动化与电网仿真技术领域。
背景技术
配电网状态估计是配电网管理***的最基础和最重要的功能之一,而正确的网络拓扑结构是保证配电网状态估计合理的前提条件。由于配电网的开关、刀闸设备缺乏实时遥信,大多采用人工维护,这就导致得到的配电网拓扑结构不可信。同时,配电网一般只在馈线根节点、重要联络开关存在实时量测,用户变配备计量终端,因此配电网量测冗余度很低,这就使得配电网拓扑结构辨识存在困难。但随着配电网管理***在配电网中逐步应用,研究一种能应用于配电网的拓扑结构辨识方法非常重要。
发明内容
本发明的目的是提出一种基于最优匹配回路功率的配电网拓扑结构辨识方法,在配电网量测数目有限的情况下,对配电网进行有效的拓扑结构辨识,为配电网的高级应用,如状态估计、潮流计算等,提供准确的拓扑结构模型。
本发明提出的基于最优匹配回路功率的配电网拓扑结构辨识方法,包括以下步骤:
(1)计算配电网中最优匹配回路功率,具体过程如下:
(1-1)设定一个回路矩阵B,用于表示配电网中独立回路与支路之间的关系,设配电网中含有n个独立回路,l条支路,则配电网的回路矩阵B为n×l阶,回路矩阵B中的元素的取值为:
式中,i表示配电网中的第i个独立回路,作为回路矩阵B的第i行,i=1…n,j表示电网的第j条支路,j=1…l,作为回路矩阵B的第j列;当支路j处于独立回路i的道路上,且支路方向与回路方向相同时bij=1;当支路j不在独立回路i上时bij=0,当支路j处于独立回路i的道路上,且支路方向与回路方向相反时,bij=-1;
(1-3)若配电网中只有一个独立回路,则配电网的最优匹配回路功率(ΔPopt,ΔQopt)为:
式中,jm为存在功率量测设备的支路,M为独立回路中的存在功率量测设备的支路集合,表示支路jm上的功率量测值的权重;和分别表示支路jm上功率量测的有功功率值和无功功率值,表示与功率量测设备所在支路相对应的回路矩阵中的元素值;
(1-4)若配电网中含有多个独立回路,各个独立回路的最优匹配回路功率 通过下面的步骤求得:
(1-4-1)根据以下的公式,计算存在功率量测设备的所有支路的最优匹配回路功率之和
式中,M′为支路jm所在的独立回路中存在支路功率量测设备的所有支路的集合;
(1-4-3)用最小二乘法求解步骤(1-4-2)中的方程组,得到各独立回路的最优匹配回路功率
(2)用向量P0、Q0分别表示配电网中所有支路的初始有功功率和初始无功功率,用向量ΔPopt、ΔQopt表示配电网中所有独立回路的最优匹配回路功率,将所有支路的初始有功功率和初始无功功率与所有独立回路的最优匹配回路功率按下面的公式进行叠加,得到叠加后的支路残留有功功率和无功功率P、Q如下;
(3)将配电网中各独立回路中的所有支路按照残留功率的大小从小到大依次排列,且将每个独立回路的前三条支路作为该独立回路中任意一条作为断开支路,若该断开支路的量测设备的量测值不为零,则将该支路作为闭合支路,若该断开支路上的量测设备的量测值为零,则将该支路作为断开支路;
(4)对步骤(3)中的断开支路进行排列组合,得到多种配电网的拓扑结构;
(5)当配电网中的任意两个独立回路之间存在公共支路,且两个独立回路中的断开支路同时为公共支路时,从步骤(4)得到的多种拓扑结构中删除由该状态确定的拓扑结构;
(6)建立一个匹配目标函数:
式中,Zreal表示配电网量测设备的实时量测值,包括支路有功功率量测值、支路无功功率量测值、配电网电压量测值和配电网电流量测值,表示对配电网进行状态估计后的配电网量测设备的量测估计值,H配电网量测设备的实时量测值集合;
(7)根据步骤(5)得到的多个配电网拓扑结构,分别用同一时刻的实时量测值对配电网进行状态估计,将状态估计得到的配电网量测设备的量测估计值和同一时刻的量测设备的实时量测值代入步骤(6)的匹配目标函数中,计算得到多个目标函数值J′,以与多个目标函数值中最小的目标函数值相对应的拓扑结构作为配电网的实际拓扑结构。
本发明提出的基于最优匹配回路功率的配电网拓扑结构辨识方法,其优点是:
1、本发明方法充分利用了配电网中的支路功率量测,可以针对量测有限的配电网进行有效的拓扑结构辨识。
2、本发明方法充分利用了计量***的负荷数据,减小了负荷伪量测的误差,进而提高了配电网拓扑结构辨识的准确度。
3、本发明方法实现简单、计算速度快,尤其适用规模不大的配电网。
附图说明
图1是基本配电网的结构示意图。
图2是含有两个独立回路的配电网结构示意图。
具体实施方式
本发明提出的基于最优匹配回路功率的配电网拓扑结构辨识方法,包括以下步骤:
(1)计算配电网中最优匹配回路功率,具体过程如下:
(1-1)设定一个回路矩阵B,用于表示配电网中独立回路与支路之间的关系,设配电网中含有n个独立回路,l条支路,则配电网的回路矩阵B为n×l阶,回路矩阵B中的元素的取值为:
式中,i表示配电网中的第i个独立回路,作为回路矩阵B的第i行,i=1…n,j表示电网的第j条支路,j=1…l,作为回路矩阵B的第j列;当支路j处于独立回路i的道路上,且支路方向与回路方向相同时bij=1;当支路j不在独立回路i上时bij=0,当支路j处于独立回路i的道路上,且支路方向与回路方向相反时,bij=-1;
以图1所示的基本配电网为例,该配电网共有7条支路,和一个独立回路,所有支路的方向为编号较小的节点指向编号较大的节点,独立回路的方向为逆时针方向。则该配电网的回路矩阵为:
(1-3)若配电网中只有一个独立回路,则配电网的最优匹配回路功率(ΔPopt,ΔQopt)为:
式中,jm为存在功率量测设备的支路,M为独立回路中的存在功率量测设备的支路集合,表示支路jm上的功率量测值的权重;和分别表示支路jm上功率量测的有功功率值和无功功率值,表示与功率量测设备所在支路相对应的回路矩阵中的元素值;
(1-4)若配电网中含有多个独立回路,各个独立回路的最优匹配回路功率 通过下面的步骤求得:
(1-4-1)根据以下的公式,计算存在功率量测设备的所有支路的最优匹配回路功率之和
式中,M′为支路jm所在的独立回路中存在支路功率量测设备的所有支路的集合;
(1-4-3)用最小二乘法求解步骤(1-4-2)中的方程组,得到各独立回路的最优匹配回路功率
(2)用向量P0、Q0分别表示配电网中所有支路的初始有功功率和初始无功功率,用向量ΔPopt、ΔQopt表示配电网中所有独立回路的最优匹配回路功率,将所有支路的初始有功功率和初始无功功率与所有独立回路的最优匹配回路功率按下面的公式进行叠加,得到叠加后的支路残留有功功率和无功功率P、Q如下;
(3)将配电网中各独立回路中的所有支路按照残留功率的大小从小到大依次排列,且将每个独立回路的前三条支路作为该独立回路中任意一条作为断开支路,若该断开支路的量测设备的量测值不为零,则将该支路作为闭合支路,若该断开支路上的量测设备的量测值为零,则将该支路作为断开支路;
(4)对步骤(3)中的断开支路进行排列组合,得到多种配电网的拓扑结构;
以图2所示的含有两个独立回路的配电网为例,经过步骤(3)处理,独立回路L1的断开支路为支路8,独立回路L2的断开支路为支路12和支路13,对独立回路L1和L2中的断开支路进行排列组合,得到两种配电网拓扑结构,一种配电网拓扑结构为支路8、支路12断开,另一种配电网拓扑结构为支路8、支路13断开。
(5)当配电网中的任意两个独立回路之间存在公共支路,且两个独立回路中的断开支路同时为公共支路时,则从步骤(4)得到的多种拓扑结构中删除由该状态确定的拓扑结构;
(6)建立一个匹配目标函数:
式中,Zreal表示配电网量测设备的实时量测值,包括支路有功功率量测值、支路无功功率量测值、配电网电压量测值和配电网电流量测值,表示对配电网进行状态估计后的配电网量测设备的量测估计值,H配电网量测设备的实时量测值集合;
(7)根据步骤(5)得到的多个配电网拓扑结构,分别用同一时刻的实时量测值对配电网进行状态估计,将状态估计得到的配电网量测设备的量测估计值和同一时刻的量测设备的实时量测值代入步骤(6)的匹配目标函数中,计算得到多个目标函数值J′,以与多个目标函数值中最小的目标函数值相对应的拓扑结构作为配电网的实际拓扑结构。
Claims (1)
1.基于最优匹配回路功率的配电网拓扑结构辨识方法,其特征在于该方法包括以下步骤:
(1)计算配电网中最优匹配回路功率,具体过程如下:
(1-1)设定一个回路矩阵B,用于表示配电网中独立回路与支路之间的关系,设配电网中含有n个独立回路,l条支路,则配电网的回路矩阵B为n×l阶,回路矩阵B中的元素的取值为:
式中,i表示配电网中的第i个独立回路,作为回路矩阵B的第i行,i=1…n,j表示电网的第j条支路,j=1…l,作为回路矩阵B的第j列;当支路j处于独立回路i的道路上,且支路方向与回路方向相同时bij=1;当支路j不在独立回路i上时bij=0,当支路j处于独立回路i的道路上,且支路方向与回路方向相反时,bij=-1;
(1-2)闭合配电网的全部开关,将配电网中的电量计费***采集的负荷数据作为负荷量测,使用基于回路分析的潮流计算方法,对配电网进行潮流计算,得到配电网中每条支路的初始支路功率
(1-3)若配电网中只有一个独立回路,则配电网的最优匹配回路功率(ΔPopt,ΔQopt)为:
式中,jm为存在功率量测设备的支路,M为独立回路中的存在功率量测设备的支路集合,表示支路jm上的功率量测值的权重;和分别表示支路jm上功率量测的有功功率值和无功功率值,表示与功率量测设备所在支路相对应的回路矩阵中的元素值;
(1-4)若配电网中含有多个独立回路,各个独立回路的最优匹配回路功率 通过下面的步骤求得:
(1-4-1)根据以下的公式,计算存在功率量测设备的所有支路的最优匹配回路功率之和
式中,M′为支路jm所在的独立回路中存在支路功率量测设备的所有支路的集合;
(1-4-3)用最小二乘法求解步骤(1-4-2)中的方程组,得到各独立回路的最优匹配回路功率
(2)用向量P0、Q0分别表示配电网中所有支路的初始有功功率和初始无功功率,用向量ΔPopt、ΔQopt表示配电网中所有独立回路的最优匹配回路功率,将所有支路的初始有功功率和初始无功功率与所有独立回路的最优匹配回路功率按下面的公式进行叠加,得到叠加后的支路残留有功功率和无功功率P、Q如下;
(3)将配电网中各独立回路中的所有支路按照残留功率的大小从小到大依次排列,且将每个独立回路的前三条支路作为该独立回路中任意一条作为断开支路,若该断开支路的量测设备的量测值不为零,则将该支路作为闭合支路,若该断开支路上的量测设备的量测值为零,则将该支路作为断开支路;
(4)对步骤(3)中的断开支路进行排列组合,得到多种配电网的拓扑结构;
(5)当配电网中的任意两个独立回路之间存在公共支路,且两个独立回路中的断开支路同时为公共支路时,则从步骤(4)得到的多种拓扑结构中删除由该状态确定的拓扑结构;
(6)建立一个匹配目标函数:
式中,Zreal表示配电网量测设备的实时量测值,包括支路有功功率量测值、支路无功功率量测值、配电网电压量测值和配电网电流量测值,表示对配电网进行状态估计后的配电网量测设备的量测估计值,H配电网量测设备的实时量测值集合;
(7)根据步骤(5)得到的多个配电网拓扑结构,分别用同一时刻的实时量测值对配电网进行状态估计,将状态估计得到的配电网量测设备的量测估计值和同一时刻的量测设备的实时量测值代入步骤(6)的匹配目标函数中,计算得到多个目标函数值J′,以与多个目标函数值中最小的目标函数值相对应的拓扑结构作为配电网的实际拓扑结构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210571490.1A CN103001220B (zh) | 2012-12-25 | 2012-12-25 | 一种基于最优匹配回路功率的配电网拓扑结构辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210571490.1A CN103001220B (zh) | 2012-12-25 | 2012-12-25 | 一种基于最优匹配回路功率的配电网拓扑结构辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103001220A true CN103001220A (zh) | 2013-03-27 |
CN103001220B CN103001220B (zh) | 2014-08-20 |
Family
ID=47929490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210571490.1A Active CN103001220B (zh) | 2012-12-25 | 2012-12-25 | 一种基于最优匹配回路功率的配电网拓扑结构辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103001220B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014198158A1 (zh) * | 2013-06-14 | 2014-12-18 | 山东理工大学 | 一种配电网控制应用网络拓扑自动识别方法 |
CN104600699A (zh) * | 2015-01-27 | 2015-05-06 | 清华大学 | 一种基于混合整数二次规划模型的配电网结构估计方法 |
CN110389269A (zh) * | 2019-07-25 | 2019-10-29 | 广州市奔流电力科技有限公司 | 基于电流优化匹配的低压台区拓扑关系识别方法及其装置 |
CN112989542A (zh) * | 2021-05-10 | 2021-06-18 | 广东电网有限责任公司中山供电局 | 一种虚拟电厂拓扑优化方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101964525A (zh) * | 2010-06-25 | 2011-02-02 | 清华大学 | 一种支持大规模电流量测的配电网状态估计方法 |
-
2012
- 2012-12-25 CN CN201210571490.1A patent/CN103001220B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101964525A (zh) * | 2010-06-25 | 2011-02-02 | 清华大学 | 一种支持大规模电流量测的配电网状态估计方法 |
Non-Patent Citations (2)
Title |
---|
HAIBIN WANG,SCHULZ, N.N.: "A revised branch current-based distribution system state estimation algorithm and meter placement impact", 《POWER SYSTEMS, IEEE TRANSACTIONS ON》 * |
吴文传等: "配网潮流回路分析法", 《中国电机工程学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014198158A1 (zh) * | 2013-06-14 | 2014-12-18 | 山东理工大学 | 一种配电网控制应用网络拓扑自动识别方法 |
CN104600699A (zh) * | 2015-01-27 | 2015-05-06 | 清华大学 | 一种基于混合整数二次规划模型的配电网结构估计方法 |
CN110389269A (zh) * | 2019-07-25 | 2019-10-29 | 广州市奔流电力科技有限公司 | 基于电流优化匹配的低压台区拓扑关系识别方法及其装置 |
CN110389269B (zh) * | 2019-07-25 | 2021-08-31 | 广州市奔流电力科技有限公司 | 基于电流优化匹配的低压台区拓扑关系识别方法及其装置 |
CN112989542A (zh) * | 2021-05-10 | 2021-06-18 | 广东电网有限责任公司中山供电局 | 一种虚拟电厂拓扑优化方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103001220B (zh) | 2014-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105791280B (zh) | 一种抵御电力***直流状态估计中数据完整性攻击的方法 | |
CN104332997B (zh) | 一种变电站三相线性广义状态估计方法 | |
CN104600699B (zh) | 一种基于混合整数二次规划模型的配电网结构估计方法 | |
CN103413044B (zh) | 一种基于变电站量测信息的电力***局部拓扑估计方法 | |
CN103093097A (zh) | 基于规范切的电力***脆弱断面识别方法 | |
CN103050969B (zh) | 一种基于馈线树的配电网状态估计方法 | |
CN103840452A (zh) | 一种引入pmu量测信息的大电网状态估计方法 | |
CN103001220B (zh) | 一种基于最优匹配回路功率的配电网拓扑结构辨识方法 | |
CN109752629B (zh) | 一种电网量测问题智能诊断方法及*** | |
CN103020726A (zh) | 面向全pmu量测的抗差状态估计方法 | |
CN103593566B (zh) | 混合二次规划形式的电力***综合状态估计方法 | |
CN104218569A (zh) | 一种大规模电网静态安全校核的评价分析方法 | |
CN103199528A (zh) | 广域电力***状态估计协调方法 | |
CN103607292B (zh) | 面向电力通信网络业务的快速分布式监测方法 | |
CN102901868A (zh) | 一种对电能量采集***数据校核的方法 | |
CN104749488A (zh) | 一种基于连续数据窗的直流线路时域故障测距方法 | |
CN107257130A (zh) | 基于区域量测解耦的低压配电网损耗计算方法 | |
CN106372440B (zh) | 一种并行计算的配电网自适应抗差状态估计方法及装置 | |
CA2602980C (en) | Method and system for real time identification of voltage stability via identification of weakest lines and buses contributing to power system collapse | |
Gao et al. | A method for the topology identification of distribution system | |
CN102280877A (zh) | 一种多量测断面的电力***不良支路参数辨识方法 | |
CN109726880B (zh) | 一种输电线路参数合理性评估方法及*** | |
CN105069707A (zh) | 一种基于正态分布统计的母线功率平衡分析方法 | |
CN102142682B (zh) | 一种基于直流潮流模型的支路开断灵敏度计算方法 | |
CN103838962A (zh) | 一种计及pmu量测的分步线性状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |