CN109193799B - 一种基于图论的配电网多种量测量的最优配置方法 - Google Patents

一种基于图论的配电网多种量测量的最优配置方法 Download PDF

Info

Publication number
CN109193799B
CN109193799B CN201811043944.1A CN201811043944A CN109193799B CN 109193799 B CN109193799 B CN 109193799B CN 201811043944 A CN201811043944 A CN 201811043944A CN 109193799 B CN109193799 B CN 109193799B
Authority
CN
China
Prior art keywords
node
distribution network
power distribution
measurement
nodes
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.)
Active
Application number
CN201811043944.1A
Other languages
English (en)
Other versions
CN109193799A (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.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power 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 State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, North China Electric Power University filed Critical State Grid Corp of China SGCC
Priority to CN201811043944.1A priority Critical patent/CN109193799B/zh
Publication of CN109193799A publication Critical patent/CN109193799A/zh
Application granted granted Critical
Publication of CN109193799B publication Critical patent/CN109193799B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明属于电力系统调度自动化技术领域,尤其涉及一种基于图论的配电网多种量测量的最优配置方法,包括:步骤1:建立配电网两节点之间有功功率和电压相角之间的量测模型;步骤2:根据配电网节点之间的连接关系形成配电网的拓扑关系图;步骤3:获得拓扑关系图上有连接的节点之间的支路有功功率量测成本以及每个节点的注入有功功率量测成本;步骤4:采用改进的Kruskal算法进行加边加点操作,逐步形成广义最小生成树作为该配电网多种量测量的最优配置方案。本发明在配置过程中具有良好的经济性,高效的计算效率以及完备的观测性,具有良好的工程应用前景,在满足整个系统可观测性下,得到了需要的最小经济投入,非常适宜于实际工程应用。

Description

一种基于图论的配电网多种量测量的最优配置方法
技术领域
本发明属于电力系统调度自动化技术领域,尤其涉及一种基于图论的配电网多种量测量的最优配置方法。
背景技术
电力系统状态估计是能量管理系统的基础和核心。现在几乎每一个大型调度中心都安装了状态估计器,状态估计已成为电网安全运行的基石。自1970年国外学者首次提出状态估计以来,人们对状态估计的研究和应用已经有40多年的历史了,这期间涌现出了各种各样的状态估计方法。在传统的配电网量测量的的配置中,人们往往是根据以往的经验进行配置,这样会带来一定的量测量重复以及不必要的成本。
发明内容
针对上述问题,本发明提出了一种基于图论的配电网多种量测量的最优配置方法,包括:
步骤1:建立配电网两节点之间有功功率和电压相角之间的量测模型;
步骤2:根据配电网节点之间的连接关系形成配电网的拓扑关系图;
步骤3:获得拓扑关系图上有连接的节点之间的支路有功功率量测成本以及每个节点的注入有功功率量测成本;
步骤4:采用改进的Kruskal算法进行加边加点操作,逐步形成广义最小生成树作为该配电网多种量测量的最优配置方案。
所述量测模型为:
Figure BDA0001792799400000011
其中Pij为节点i与节点j之间的有功交换功率,θi与θj分别为节点i和节点j的相角,xij为节点i与节点j之间的支路电抗。
所述两节点之间构成一个量测岛,若获取了某一个节点的相角和节点之间传递的有功功率,则该量测岛可观;若与一个节点连接的支路数大于等于引入该节点后系统孤立岛数时,则整个系统可观。
所述步骤4包括:若当前有功功率量测成本最小为支路有功功率则保证支路上的两个节点不在同一颗树上,然后合并两节点;若当前有功功率量测成本最小为节点注入功率,则保证当前节点的支路数大于等于引入节点后当前量测孤岛的数目,然后合并两个量测孤岛;从而形成该配电网的广义最小生成树,以及得到最小量测成本和配电网多种量测最优配置方案。
本发明的有益效果:本发明在配置过程中具有良好的经济性,高效的计算效率以及完备的观测性,具有良好的工程应用前景,在满足整个系统的可观测性下,得到了需要的最小的经济投入,非常适宜于实际工程应用。
附图说明
图1为实施例中配电网节点拓扑关系图。
图2为实施例中节点1、2之间的支路双向有功功率量测成本图。
图3为实施例中节点1、3之间的支路双向有功功率量测成本图。
图4为实施例中节点2、3之间的支路双向有功功率量测成本图
图5为实施例中节点3、4之间的支路双向有功功率量测成本图
图6为实施例中节点3、5之间的支路双向有功功率量测成本图
图7为实施例中节点4、6之间的支路双向有功功率量测成本图
图8为实施例中节点5、6之间的支路双向有功功率量测成本图
图9为实施例中节点1的注入有功功率量测成本图
图10为实施例中节点2的注入有功功率量测成本图
图11为实施例中节点3的注入有功功率量测成本图
图12为实施例中节点4的注入有功功率量测成本图
图13为实施例中节点5的注入有功功率量测成本图
图14为实施例中节点6的注入有功功率量测成本图
图15为实施例中配电网节点可观测图
图16为实施例中Kruskal算法第一次加边
图17为实施例中Kruskal算法第二次加边
图18为实施例中Kruskal算法第三次加边
图19为实施例中Kruskal算法第四次加边
图20为实施例中Kruskal算法第五次加边
图21为实施例中配电网多种量测量的最优配置方案
图22为本发明的配电网多种量测量的最优配置方案流程图
具体实施方式
下面结合附图,对实施例作详细说明。
如图22所示,本发明配电网量测量的最优配置方案(Optimal configurationscheme for distribution network measurement)包括下列步骤:
步骤A:提供一种通过有功功率量测电压相角,无功功率量测电压幅值的量测模型;具体的,本发明的提出通过有功功率量测电压相角,无功功率量测电压幅值的量测模型;如下所示:
Figure BDA0001792799400000031
其中Pij为节点i与节点j之间的有功交换功率,θi与θj分别为节点i和节点j的相角,xij为节点i与节点j之间的支路电抗。根据电力系统基础知识,电压幅值与无功功率的关系显然易见。
步骤B:利用改进最小生成树算法,对所述配电网多种量测量的最优配置方案求解。
注意到配电网多种量测量的最优配置方案是一个基于图论的问题,适宜用一种改进过的最小生成树算法进行求解。为使本领域技术人员更好地理解本发明,首先最优配置方案求解的推导过程如下:
如图1所示为配电网节点拓扑关系图。获得拓扑关系图上面有连接的节点之间的支路有功功率量测成本以及每个节点的注入有功功率量测成本,如图2~图14所示
因为有功功率与无功功率之间的关系是仅仅数学上正余弦的关系,因此对于支路功率,只用考虑有功功率的双向量测成本,其中Pij代表量测从节点i(i=1,2,3,4,5,6)流向节点j(j=1,2,3,4,5,6)的有功功率量测成本,Pji代表量测从节点j(j=1,2,3,4,5,6)流向节点i(i=1,2,3,4,5,6)的有功功率量测成本,ij两个节点之间通过有功功率和无功功率构成了一个对于向角和电压幅值的可观测岛;对于节点注入功率,只用考虑有功功率的注入量测成本,其中Pi即代表量测从节点i(i=1,2,3,4,5,6)注入系统的有功功率量测成本,如果在构建最小生成树时,引入节点i以后,节点i的支路数大于等于当前系统的孤岛数,则i节点通过有功功率和无功功率与之前的系统构成了一个对于电压相角和幅值的可观测岛。
如图15所示为配电网节点可观测图。由图可见,不同的节点之前可以通过支路功率或者节点注入功率构成一个可观测图,保证一个全局可观的状态,但不同的节点之间还存在一定的冗余观测。
采用一种改进的Kruskal算法(加边算法),逐步形成一种广义最小生成树,每一步都确保最新系统的可观测性,该广义最小生成树既是该配电网多种量测量的最优配置方案
在配电网多种量测量的最优配置方案的求解推导过程之后,将最小生成树的求解步骤归纳如下:
步骤B41:选择当前成本最小的量测量(节点1,2之间的支路有功功率,P21=1);如果是支路有功功率并保证支路上的两个节点(1,2)不在同一颗树上,然后合并节点(1,2),如图16所示;
步骤B42:选择当前成本最小的量测量(节点1的注入有功功率,P1=2.5);如果是节点注入功率,并保证当前节点的支路数(当前支路数为2)大于等于引入节点后当前量测孤岛的数目(量测孤岛1(1,2节点),量测孤岛2(节点3)),然后合并两个量测孤岛,如图17所示:
步骤B43:选择当前成本最小的量测量(节点3,4之间的支路有功功率,P34=7);如果是支路有功功率并保证支路上的两个节点(3,4)不在同一颗树上,然后合并节点(3,4),如图18所示:
步骤B44:选择当前成本最小的量测量(节点3的注入有功功率,P3=6.5);如果是节点注入功率,并保证当前节点的支路数(当前支路数为4)大于等于引入节点后当前量测孤岛的数目(量测孤岛1(1,2,3,4节点),量测孤岛2(节点5)),然后合并两个量测孤岛,如图19所示:
步骤B45:选择当前成本最小的量测量(节点4的注入有功功率,P4=9.5);如果是节点注入功率,并保证当前节点的支路数(当前支路数为2)大于等于引入节点后当前量测孤岛的数目(量测孤岛1(1,2,3,4,5节点),量测孤岛2(节点6)),然后合并两个量测孤岛,如图20所示:
步骤B46:形成该配电网的广义最小生成树,以及得到最小量测成本(P21+P1+P3+P43+P4=1+2.5+6.5+7+9.5=26.5)和配电网多种量测最优配置方案,如图21所示。
此实施例仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (1)

1.一种基于图论的配电网多种量测量的最优配置方法,其特征在于,包括:
步骤1:建立配电网两节点之间有功功率和电压相角之间的量测模型;
步骤2:根据配电网节点之间的连接关系形成配电网的拓扑关系图;
步骤3:获得拓扑关系图上有连接的节点之间的支路有功功率量测成本以及每个节点的注入有功功率量测成本;
步骤4:采用改进的Kruskal算法进行加边加点操作,逐步形成广义最小生成树作为该配电网多种量测量的最优配置方案;
所述量测模型为:
Figure FDA0003242971430000011
其中Pij为节点i与节点j之间的有功交换功率,θi与θj分别为节点i和节点j的相角,xij为节点i与节点j之间的支路电抗;
所述两节点之间构成一个量测岛,若获取了某一个节点的相角和节点之间传递的有功功率,则该量测岛可观;若与一个节点连接的支路数大于等于引入该节点后系统孤立岛数时,则整个系统可观;
所述步骤4包括:若当前有功功率量测成本最小为支路有功功率则保证支路上的两个节点不在同一颗树上,然后合并两节点;若当前有功功率量测成本最小为节点注入功率,则保证当前节点的支路数大于等于引入节点后当前量测孤岛的数目,然后合并两个量测孤岛;从而形成该配电网的广义最小生成树,以及得到最小量测成本和配电网多种量测最优配置方案。
CN201811043944.1A 2018-09-07 2018-09-07 一种基于图论的配电网多种量测量的最优配置方法 Active CN109193799B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811043944.1A CN109193799B (zh) 2018-09-07 2018-09-07 一种基于图论的配电网多种量测量的最优配置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811043944.1A CN109193799B (zh) 2018-09-07 2018-09-07 一种基于图论的配电网多种量测量的最优配置方法

Publications (2)

Publication Number Publication Date
CN109193799A CN109193799A (zh) 2019-01-11
CN109193799B true CN109193799B (zh) 2022-01-04

Family

ID=64915290

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811043944.1A Active CN109193799B (zh) 2018-09-07 2018-09-07 一种基于图论的配电网多种量测量的最优配置方法

Country Status (1)

Country Link
CN (1) CN109193799B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105514970A (zh) * 2015-10-22 2016-04-20 国家电网公司 一种基于图论的upfc监控子站布点算法
CN107017631A (zh) * 2017-05-22 2017-08-04 国网福建省电力有限公司 一种基于线性电路的三相配电网可观测性分析方法
CN107453351A (zh) * 2017-07-12 2017-12-08 河海大学 一种基于节点注入功率的配电网运行拓扑辨识方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102982379B (zh) * 2012-11-16 2016-08-03 清华大学 基于集合论估计模型的电力系统状态估计方法及系统
CN104332997B (zh) * 2014-11-18 2016-03-02 国电南瑞科技股份有限公司 一种变电站三相线性广义状态估计方法
CN104600699B (zh) * 2015-01-27 2016-08-17 清华大学 一种基于混合整数二次规划模型的配电网结构估计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105514970A (zh) * 2015-10-22 2016-04-20 国家电网公司 一种基于图论的upfc监控子站布点算法
CN107017631A (zh) * 2017-05-22 2017-08-04 国网福建省电力有限公司 一种基于线性电路的三相配电网可观测性分析方法
CN107453351A (zh) * 2017-07-12 2017-12-08 河海大学 一种基于节点注入功率的配电网运行拓扑辨识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
最小生成树算法在配电网重构中的应用;李超 等;《四川电力技术》;20141231;第37卷(第6期);第1-4页 *
电力系统中PMU最优配置的研究;蔡田田 等;《电网技术》;20060731;第30卷(第13期);第32-37页 *

Also Published As

Publication number Publication date
CN109193799A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
Schenato et al. Bayesian linear state estimation using smart meters and pmus measurements in distribution grids
CN101924364B (zh) 变电站-调度中心两级分布式电网的非线性状态估计方法
CN102403720B (zh) 一种基于暂态电压安全裕度的超实时重合时序整定方法
CN102521452B (zh) 一种大电网合环计算系统
CN103632313A (zh) 一种基于pmu数据的电网动态可观方法
CN110797874B (zh) 含电力电子变压器的交直流混合配电网状态估计方法
CN108933445A (zh) 基于回路电流法的直流电网短路电流计算方法
CN109918762B (zh) 一种交直流电网混合仿真对比模型构建方法
CN103020726A (zh) 面向全pmu量测的抗差状态估计方法
CN102420427A (zh) 一种考虑外网等值的区域电网电压稳定裕度计算方法
CN103077480A (zh) 一种电力系统的安全校核方法
CN103954885A (zh) 基于分布参数的双回线故障单端定位系统及定位方法
CN107392498A (zh) 一种电网拓扑分析系统
CN103018633B (zh) 基于不对称阻抗模型短路电流故障特征的含dg配电网故障区间判定方法
CN104113061A (zh) 一种含分布式电源的配电网三相潮流计算方法
CN109858061B (zh) 面向电压功率灵敏度估计的配电网络等值与化简方法
CN103595039B (zh) 电网系统分布式合环风险分析方法
CN107046285A (zh) 一种基于混合量测的配电网状态评估方法
CN109193799B (zh) 一种基于图论的配电网多种量测量的最优配置方法
CN104573384B (zh) 一种适用于配电网在线分析的改进直流潮流算法
CN104007364A (zh) 一种基于电压相位关系的单端行波测距方法
CN112234598B (zh) 一种电磁暂态仿真初始化方法
CN107196307B (zh) 一种变压器故障后电网有功潮流快速估算的方法
CN105514978A (zh) 一种混合整数二次规划形式的抗差状态估计方法
CN104505866A (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