CN113591258A - 基于节点关联矩阵的配电网络pmu配置优化方法 - Google Patents
基于节点关联矩阵的配电网络pmu配置优化方法 Download PDFInfo
- Publication number
- CN113591258A CN113591258A CN202110871520.XA CN202110871520A CN113591258A CN 113591258 A CN113591258 A CN 113591258A CN 202110871520 A CN202110871520 A CN 202110871520A CN 113591258 A CN113591258 A CN 113591258A
- Authority
- CN
- China
- Prior art keywords
- node
- distribution network
- pmu
- matrix
- 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.)
- Pending
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 186
- 239000011159 matrix material Substances 0.000 title claims abstract description 112
- 238000005457 optimization Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000006870 function Effects 0.000 claims description 56
- 230000002776 aggregation Effects 0.000 claims description 35
- 238000004220 aggregation Methods 0.000 claims description 26
- 238000004590 computer program Methods 0.000 claims description 12
- 238000005054 agglomeration Methods 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 8
- 238000009434 installation Methods 0.000 claims description 7
- 230000002194 synthesizing effect Effects 0.000 claims description 5
- 238000004891 communication Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 3
- 238000010248 power generation Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 5
- 230000008602 contraction Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011900 installation process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Business, Economics & Management (AREA)
- Evolutionary Computation (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Human Resources & Organizations (AREA)
- Primary Health Care (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Marketing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于节点关联矩阵的配电网络PMU配置优化方法，包括：获取配电网络的拓扑数据，构建配电网络系统模型；获取配电网络中的电源节点，构建电源节点矩阵；根据配电网络系统模型中的邻接矩阵和电源节点矩阵，构建配电网络系统的节点关联矩阵；构建配电网络PMU配置的多目标优化模型，求解所述配电网络PMU配置的多目标优化模型，得到配电网络PMU配置的最优方案。本方法可以用尽可能少的PMU来对整个网络的所有节点进行观测，且配置PMU的费用较少；该方法突出了配电网络的拓扑辨识，同时反映了配电网络中的多种分布式新能源发电站、每个用电负荷地理位置与连接情况，为分布式新能源发电及智能电网调度运行提供必要的技术支撑。
Description
技术领域
本发明涉及电力优化领域，具体涉及一种基于节点关联矩阵的配电网络PMU配置优化方法。
背景技术
随着用电用户的增多，以及供电可靠性要求的提高，配电网络互联规模在不断扩大，配电系统的复杂程度也在持续增加。由于接入配电网络的分布式电源以及用户负荷的不确定性，对配电网络的安全稳定运行、供电可靠性带来了极大的挑战，例如：潮流双向化、电压波动加剧等。为了很好地解决这些问题，就必须要求配电网络具有良好的可观测性和可控制性。
当前通信技术的迅猛发展，5G数据通信已经实现商业化，其可靠性高、延时性地的优点为PMU在电网中的广泛使用打下了坚实的基础。同步向量量测单元(PhasorMeasurement Unit，PMU)，以全球定位系统(Global Position System，GPS)为时间基准，可提供高精度、带时标的电压、电流、以及频率信号，因其相量特性、时钟同步特性、以及数据上传的实时特性，而成为电网动态过程监测的基础手段。
当前，PMU已广泛应用于输电网中，但并未广泛应用于配电网络中，这是因为在配电网络中存在数量众多的节点与支路，配电网络结构复杂，且配电网络的规模以及复杂程度正在持续上升。所以想要通过PMU来对整个配电网络进行实时的监测与控制，就必须要在配电网络中大规模地配置PMU，但是考虑到PMU的造价高、以及大规模地配置PMU在短时间内无法完成，所以在配电网络中大规模地配置PMU并不现实，而且配电网络中某个节点配置PMU后，与之相连接的节点都是可观测的，这些节点都不用配置PMU，因此在配电网络中所有节点处都配置PMU也是没有必要的。只需要在配电网络中最关键的节点处安装PMU就能对整个网络的节点起到实时监测和控制，并且还能保持较低的PMU配置费用。
配电网络PMU配置优化的过程中最主要解决好以下三个方面的问题：(1)使配电网络系统中处于PMU实时监测和控制的节点数最多，最终目标是使网络中所有节点都处于PMU的实时监测和控制；(2)使配电网络中安装的PMU数量最少，即用数量最少的PMU来实时监测和控制最多的节点，保证配电网络PMU配置的费用最少；(3)配电网络要保证较高的供电可靠性，所以在安装PMU的同时要使安装的时间短，又或者对配电网络中的分布式能源、联络开关进行调度，保证在PMU安装过程中供电的可靠性。
由于配电网络结构的日趋复杂，以及分布式能源的大规模接入，配电网络的供电可靠性、以及安全稳定的运行受到严峻的考验；所以要在配电网络中配置PMU，对节点进行实时的监测和控制，保证配电网络安全稳定地运行。在配电网络配置PMU的过程中，在配电网络中的哪些位置配置PMU，并在这过程中重点考虑配电网络配置PMU的费用最少、以及配电网络配置PMU后可实时监测和控制的节点数最多这两个方面，成为配电网络调度面临的重要问题。
发明内容
本发明的目的是以配置PMU的费用最少、配置PMU后可实时监测和控制节点最多为目标，提供一种基于节点关联矩阵的配电网络PMU配置优化方法。
为了实现上述任务，本发明采用以下技术方案：
一种基于节点关联矩阵的配电网络PMU配置优化方法，包括：
获取配电网络的拓扑数据，构建配电网络系统模型，包括配电网络中节点、支路关系的关联矩阵及其邻接矩阵；获取配电网络中的电源节点，构建电源节点矩阵；
根据配电网络系统模型中的邻接矩阵和电源节点矩阵，构建配电网络系统的节点关联矩阵；
构建配电网络PMU配置的多目标优化模型，包括：
以配电网络中配置PMU费用最少、可观测节点的数量最大为优化目标分别构建对应的目标函数，并对两种情况下构建的目标函数进行综合得到总的目标函数；
以节点的关联度值表征节点的关键度，构建节点关联度公式；以节点的凝聚度表征节点是否位于配电网络拓扑的中心，构建节点凝聚度公式；将节点关联度公式、节点凝聚度公式通过不同的权重系数结合起来，构建节点的关键度函数；
采用群体智能优化算法求解所述配电网络PMU配置的多目标优化模型，得到配电网络PMU配置的最优方案。
进一步地，所述根据配电网络系统模型中的邻接矩阵和电源节点矩阵，构建配电网络系统的节点关联矩阵，包括：
B＝C+T
其中，C表示邻接矩阵，T表示电源节点矩阵；
所述邻接矩阵C中的元素cij的取值：若节点i与节点j相连则取值为1，否则为0，对角线元素规定为1；所述电源节点矩阵T中元素tii的取值：若节点是电源节点，则取值为1，否则为0；
所述节点关联矩阵B中的元素bij的取值：若节点i与节点j相连则取值为1，否则为0；对于对角线元素，若节点是电源节点，则取值为2，否则取值为1。
进一步地，以配电网络中配置PMU费用最少构建目标函数，具体为：
其中，n表示配电网络系统中节点的总数；ci为第i个节点的权重系数，表示在第i个节点处安装PMU的费用；ki表示第i个节点配置PMU的状态变量，其值为0时，表示第i个节点处没有配置PMU，其值为1时，表示第i个节点处有配置PMU。
进一步地，以可观测节点的数量最大为优化目标构建目标函数，具体为：
其中，gi表示第i个节点的可观测性，如下式所示：
G＝B·K
其中，K＝[k1，k2，…，kn]T为n维列向量，ki表示第i个节点配置PMU的状态变量，其值为0时，表示第i个节点处没有配置PMU，其值为1时，表示第i个节点处有配置PMU；B为n×n阶的节点关联矩阵；G＝[g1，g2，…，gn]T为n维列向量，gi表示第i个节点可观测性的状态变量，其值为0时，表示第i个节点不可以观测，其值为1时，表示第i个节点可以观测。
进一步地，所述对两种情况下构建的目标函数进行综合得到总的目标函数，表示为：
max F＝F2-F1。
进一步地，所述构建节点关联度公式，具体为：
其中，D1(i)表示节点i的关联度；bij为节点关联矩阵中的元素，表示节点i与节点j的连接情况。
进一步地，所述构建节点凝聚度公式，具体包括：
将配电网络表示为图G(V,E)，其中G为无向网络连通图，其中有n个节点，m条支路，V表示节点集合，E表示支路的集合；首先定义图G的凝聚度，为：
其中，
表示网络中所有节点对最小距离的平均值；以Gi表示为节点i收缩后的图，则定义节点i的凝聚度公式，为：
其中，D2(Gi)表示节点i收缩后配电网络的凝聚度。
进一步地，所述节点i收缩是指将网络中与节点i相连接的k个节点都与节点i融合，用一个新节点代替这k+1个节点，原先与k个节点相连接的节点都与这个新节点相连接。
进一步地，所述将节点关联度公式、节点凝聚度公式结合起来，构建节点的关键度函数，表示为：
D(i)＝W1×D1(i)+W2×D2(i)
其中，D(i)表示节点i的关键度，节点关键度的值越大，就证明该节点在网络中处于关键、核心的位置，应该在此节点处配置PMU；W1、W2分别为权重系数，由熵值法来确定。
进一步地，所述采用群体智能优化算法求解所述配电网络PMU配置的多目标优化模型，得到配电网络PMU配置的最优方案，包括：
根据配电网络的拓扑数据，依次确定关联矩阵、邻接矩阵、电源节点矩阵以及节点关联矩阵后，依据节点关联矩阵，将电源节点作为配置PMU优先级最高的节点；然后根据关联度公式、凝聚度公式，求取配电网络中所有节点的关联度和凝聚度，形成数据矩阵，通过熵值法求取关联度、凝聚度的权重系数，再通过关键度函数计算每个节点的关键度；
将节点配置PMU的优先级按照节点的关键度从大到小进行排列，对于关键度相同的节点，利用群体智能优化算法进行排列组合，这样就会得到所有节点的多组配置PMU优先级不同的组合；
将所述组合的PMU安装位置以及各PMU的安装成本导入所述总的目标函数中，利用群体智能优化算法求解出这些组合的目标函数值，则求解出的目标函数值最大的组合即为配电网络PMU配置的最优方案。
一种基于节点关联矩阵的配电网络PMU配置优化装置，包括：
配电网络系统模型构建模块，用于获取配电网络的拓扑数据，构建配电网络系统模型，包括配电网络中节点、支路关系的关联矩阵及其邻接矩阵；获取配电网络中的电源节点，构建电源节点矩阵；
节点关联矩阵构建模块，用于根据配电网络系统模型中的邻接矩阵和电源节点矩阵，构建配电网络系统的节点关联矩阵；
多目标优化模型构建模块，用于构建配电网络PMU配置的多目标优化模型，包括总目标函数构建单元以及约束条件构建单元，其中：
总目标函数构建单元用于以配电网络中配置PMU费用最少、可观测节点的数量最大为优化目标分别构建对应的目标函数，并对两种情况下构建的目标函数进行综合得到总的目标函数；
约束条件构建单元用于以节点的关联度值表征节点的关键度，构建节点关联度公式；以节点的凝聚度表征节点是否位于配电网络拓扑的中心，构建节点凝聚度公式；将节点关联度公式、节点凝聚度公式通过不同的权重系数结合起来，构建节点的关键度函数；
优化问题求解模块，用于采用群体智能优化算法求解所述配电网络PMU配置的多目标优化模型，得到配电网络PMU配置的最优方案。
一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，处理器执行计算机程序时实现前述基于节点关联矩阵的配电网络PMU配置优化方法的步骤。
一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现前述基于节点关联矩阵的配电网络PMU配置优化方法的步骤。
与现有技术相比，本发明具有以下技术特点：
本发明方法可以用尽可能少的PMU来对整个网络的所有节点进行观测，并且保证配置PMU的费用较少；该方法突出了配电网络的拓扑辨识，同时反映了配电网络中的多种分布式新能源发电站、每个用电负荷地理位置与连接情况，以及网络中所有节点的关键度，为分布式新能源发电及智能电网调度运行提供必要的技术支撑。
附图说明
图1为本发明方法的流程示意图。
具体实施方式
在配电网络中配置PMU是智能化配电网络的重要途经之一，PMU具有高频率采集配电网络中数据的能力，数据的延迟性低，所得数据抗干扰能力强，能为配电网络的实时监测与控制带来很好的帮助。通过在配电网络上每个节点处都配置PMU能很好地实现网络完全可观；但由于PMU的配置费用高昂，以及大规模配置PMU在短期内难以实现等原因，使得在每个节点处都配置PMU的方案难以实现；另外，配电网络中某个节点配置PMU后，与之相连接的节点都是可观测的，这些节点都不用配置PMU，所以在配电网络中所有节点处都配置PMU也是没有必要的。
本发明的思路是，建立考虑了配电网络PMU配置的经济性、以及配电网络的全局可观测性的多目标配电网络PMU优化配置模型，通过寻找配电网络中最关键、最核心的节点，在这些节点处配置PMU，就能保证观测节点数量最多且PMU的配置费用最少，在这个过程中寻找最优的配电网络配置PMU的方案。为此，本发明引入节点关联度、节点凝聚度、节点关键度的定义，以在配电网络中配置PMU所用的费用最少、以及配电网络中可观测节点数最多为优化目标，提出基于节点关联矩阵的配电网络PMU配置优化方法，利用节点关联矩阵计算网络中所有节点的关键度，在配电网络中最关键的节点处配置PMU，通过这样的方式，保证了配置PMU费用低的同时还保证了处于可观测节点的数量多，这种PMU配置方案就为最优的PMU配置方案。
参见图1，本发明提供的一种基于节点关联矩阵的配电网络PMU配置优化方法，包括以下步骤：
S1配电网络系统模型的构建
获取配电网络的拓扑数据，构建配电网络系统模型，包括配电网络中节点、支路关系的关联矩阵A及其邻接矩阵C；其中，所述的节点包括配电网络中的变压器、负荷以及分布式电源等；所述的支路包括配电网络中的馈电线路。
C＝A×AT
其中，参数n表示节点数量，矩阵A中的元素aij的取值：若节点i与支路j相连，则aij＝1；若节点i与支路j不相连，则aij＝0；矩阵C中的元素cij的取值：若节点i与节点j相连，则cij＝1；若节点i与节点j不相连，则cij＝0；而对于对角线元素，规定为1，即cii＝cjj＝1。
然后通过IEEE的数据获取配电网络中电源节点，形成电源节点矩阵T：
其中，矩阵T中的元素tii的取值：若节点是电源节点，则tii＝1，若节点不是电源节点，则tii＝0。
S2节点关联矩阵的构建
根据配电网络系统模型中的邻接矩阵C和电源节点矩阵T，构建配电网络系统的节点关联矩阵：
B＝C+T
其中，节点关联矩阵B表示节点间的连接关系，矩阵B中的元素bij的取值：若节点i与节点j相连，则bij＝1；若节点i与节点j不相连，则bij＝0；而对于对角线元素：若节点是电源节点，则bii＝2，若节点不是电源节点，则bii＝1。
S3配电网络PMU配置的多目标优化模型的构建
S3.1总目标函数的构建
以配电网络中配置PMU费用最少、可观测节点的数量最大为优化目标分别构建对应的目标函数，并对两种情况下构建的目标函数进行综合得到总的目标函数。
(1)以配电网络中配置PMU的费用最少为优化目标：
在配电网络系统中节点、线路、线路分支众多，考虑到经济性、以及获得数据的处理量，当前在整个配电网络系统中的每个节点都配置PMU并不现实，需要对PMU的配置进行优化，所以本发明的其中一个目标是使配电网络系统中所配置PMU的费用最少，其目标函数为：
其中，n表示配电网络系统中节点的总数；ci为第i个节点的权重系数，表示在第i个节点处安装PMU的成本，根据实际安装PMU的复杂程度和时间等进行确定；ki表示第i个节点配置PMU的状态变量，其值为0时，表示第i个节点处没有配置PMU，其值为1时，表示第i个节点处有配置PMU。
(2)以配电网络中可观测节点的数量最大为优化目标：
以配电网络中可观测节点的数量最大为优化目标：配电网络通过配置PMU获得较好的可观性，在配电网络系统发生故障而导致大面积停电后，可以迅速地找到故障区域并将其隔离，再找到最近的联络开关使非故障停电区域的停电负荷尽可能的恢复供电，提高供电可靠性，其目标函数为：
其中，n表示配电网络系统中节点的总数；gi表示第i个节点的可观测性，如下式所示：
G＝B·K
其中，K＝[k1，k2，…，kn]T为n维列向量，ki表示第i个节点配置PMU的状态变量，其值为0时，表示第i个节点处没有配置PMU，其值为1时，表示第i个节点处有配置PMU；B为n×n阶的节点关联矩阵；G＝[g1，g2，…，gn]T为n维列向量，gi表示第i个节点可观测性的状态变量，其值为0时，表示第i个节点不可以观测，其值为1时，表示第i个节点可以观测。
(3)综合两个目标函数：
将以配电网络中配置PMU的费用最少为优化目标的目标函数、以配电网络中可观测节点的数量最大为优化目标的目标函数相结合起来，可以得到总的目标函数为：
max F＝F2-F1
S3.2约束条件的构建
以节点的关联度值表征节点的关键度，构建节点关联度公式；以节点的凝聚度表征节点是否位于配电网络拓扑的中心，构建节点凝聚度公式；将节点关联度公式、节点凝聚度公式通过不同的权重系数结合起来，构建节点的关键度函数。
(1)节点的关联度约束条件：
节点的关联度值表示在整个网络中与该节点相连节点的数目。在本发明中，节点的关联度值越大，则认为该节点的关键度越大，配置优先级越高。节点关联度公式为：
其中，D1(i)表示节点i的关联度；bij为节点关联矩阵中的元素，表示节点i与节点j的连接情况。
(2)节点的凝聚度约束条件：
由于配电网络呈放射状的结构特点，导致很多节点的关联度相同，因此节点的关联度值并不能全面衡量节点的关键程度，所以本发明还引入了节点的凝聚度作为判别网络中节点关键程度的另一指标。将配电网络表示为图G(V,E),其中G为无向网络连通图，其中有n个节点，m条支路，V表示节点集合，E表示支路的集合。首先定义图G的凝聚度，为：
其中，
表示网络中所有节点对最小距离的平均值，可以根据网络的节点关联矩阵，运用Dijkstra算法求出网络中所有节点对的最小距离dij，然后就可以求取网络中所有节点对最小距离的平均值。以Gi表示为节点i收缩后的图；其中节点i收缩是指将网络中与节点i相连接的k个节点都与节点i融合，即用一个新节点代替这k+1个节点，原先与k个节点相连接的节点都与这个新节点相连接。
则定义节点i的凝聚度公式，为：
其中，D2(Gi)表示节点i收缩后配电网络的凝聚度；如果节点i处于核心位置，即许多节点对之间的最短路径都要通过这个节点，那么将节点i收缩以后网络的
就会大大减小，从而获得较大的网络凝聚度。
节点的凝聚度表示节点是否处于配电网络拓扑的中心，凝聚度越大，则节点处于网络中的核心位置。
(3)综合两个约束条件：
将节点的关联度公式、节点的凝聚度公式相结合起来，可以得到节点的关键度函数为：
D(i)＝W1×D1(i)+W2×D2(i)
其中，D(i)表示节点i的关键度，节点关键度的值越大，就证明该节点在网络中处于关键、核心的位置，应该在此节点处配置PMU；W1、W2分别为权重系数，由熵值法来确定。
(4)熵值法确定权重系数：
在信息论中，熵是对不确定性的一种度量。信息量越大，不确定性就越小，熵也就越小；信息量越小，不确定性越大，熵也越大。根据熵的特性，我们可以通过计算熵值来判断一个事件的随机性及无序程度，也可以用熵值来判断某个指标的离散程度，指标的离散程度越大，该指标对综合评价的影响越大。
1)数据矩阵：
其中，Xij表示第i个方案的第j项指标的数值，如：Xn1表示第n个节点的节点关联度，Xn2表示第n个节点的节点凝聚度。
2)计算第j项指标下的第i个方案的特征比重：
3)计算第j项指标的熵值：
4)计算第j项指标的差异系数：
hj＝1-ej
对于第j项指标，指标值Xij的差异越大，对方案评价的作用越大，熵值就越大。
5)计算第j项指标的熵权值：
S4优化问题的求解
采用群体智能优化算法，如粒子群算法、遗传算法、模拟退火算法等，可以求解出基于节点关联矩阵的配电网络PMU配置的多目标优化模型。
实施例：
首先导入配电网络的拓扑数据，获取关联矩阵A、电源节点矩阵T，根据步骤S1的公式求取邻接矩阵C，然后根据步骤S2的公式，求取节点关联矩阵B，矩阵B对角线元素值为2表示配电网络中的电源节点，将电源节点认为是配置PMU优先级最高的节点，然后根据步骤S3.2中的公式，求取配电网络中所有节点的关联度和凝聚度，形成数据矩阵S，然后将数据矩阵S导入熵值法中求取权重系数W1、W2，再根据步骤S3.2中的公式计算配电网络中节点的关键度，这些节点配置PMU的优先级按节点的关键度从大到小进行排列，节点关键度越大的节点，配置优先级越高；由于配电网络中节点数量多、以及配电网络呈辐射状，所以有些节点的关键度相同，对于这些节点的配置PMU的优先级，利用群体智能优化算法进行排列组合，这样就会得到所有节点的多组配置PMU优先级不同的组合；当网络中的某个节点配置PMU后，与该节点相连接的节点都处于可观测状态，这些节点不需要再额外装设PMU，将节点关联矩阵中表示这些节点的行向量变为零元素向量，当整个节点关联矩阵变为零元素矩阵时则停止；
将所述组合的PMU安装位置以及各PMU的安装成本导入到S3.1的总的目标函数中，利用群体智能优化算法求解出这些组合的目标函数值，则求解出的目标函数值最大的组合即为配电网络PMU配置的最优方案。
根据本发明的另一方面，提供了一种基于节点关联矩阵的配电网络PMU配置优化装置，包括：
配电网络系统模型构建模块，用于获取配电网络的拓扑数据，构建配电网络系统模型，包括配电网络中节点、支路关系的关联矩阵及其邻接矩阵；获取配电网络中的电源节点，构建电源节点矩阵；
节点关联矩阵构建模块，用于根据配电网络系统模型中的邻接矩阵和电源节点矩阵，构建配电网络系统的节点关联矩阵；
多目标优化模型构建模块，用于构建配电网络PMU配置的多目标优化模型，包括总目标函数构建单元以及约束条件构建单元，其中：
总目标函数构建单元用于以配电网络中配置PMU费用最少、可观测节点的数量最大为优化目标分别构建对应的目标函数，并对两种情况下构建的目标函数进行综合得到总的目标函数；
约束条件构建单元用于以节点的关联度值表征节点的关键度，构建节点关联度公式；以节点的凝聚度表征节点是否位于配电网络拓扑的中心，构建节点凝聚度公式；将节点关联度公式、节点凝聚度公式通过不同的权重系数结合起来，构建节点的关键度函数；
优化问题求解模块，用于采用群体智能优化算法求解所述配电网络PMU配置的多目标优化模型，得到配电网络PMU配置的最优方案。
需要说明的是，上述各个模块的具体功能和相关解释参见前述方法实施例中对应的步骤S1至S4，在此不赘述。
本申请实施例进一步提供一种终端设备，该终端设备可以为计算机、服务器；包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述基于节点关联矩阵的配电网络PMU配置优化方法的步骤，例如，前述的S1至S4。
计算机程序也可以被分割成一个或多个模块/单元，一个或者多个模块/单元被存储在存储器中，并由处理器执行，以完成本申请。一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述计算机程序在终端设备中的执行过程，例如，计算机程序可以被分割为配电网络系统模型构建模块、节点关联矩阵构建模块、多目标优化模型构建模块以及优化问题求解模块，各模块的功能参见前述装置中的描述，不再赘述。
本申请的实施提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现上述基于节点关联矩阵的配电网络PMU配置优化方法的步骤，例如，前述的S1至S4。
以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。
Claims (10)
1.一种基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，包括：
获取配电网络的拓扑数据，构建配电网络系统模型，包括配电网络中节点、支路关系的关联矩阵及其邻接矩阵；获取配电网络中的电源节点，构建电源节点矩阵；
根据配电网络系统模型中的邻接矩阵和电源节点矩阵，构建配电网络系统的节点关联矩阵；
构建配电网络PMU配置的多目标优化模型，包括：
以配电网络中配置PMU费用最少、可观测节点的数量最大为优化目标分别构建对应的目标函数，并对两种情况下构建的目标函数进行综合得到总的目标函数；
以节点的关联度值表征节点的关键度，构建节点关联度公式；以节点的凝聚度表征节点是否位于配电网络拓扑的中心，构建节点凝聚度公式；将节点关联度公式、节点凝聚度公式通过不同的权重系数结合起来，构建节点的关键度函数；
采用群体智能优化算法求解所述配电网络PMU配置的多目标优化模型，得到配电网络PMU配置的最优方案。
2.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，所述根据配电网络系统模型中的邻接矩阵和电源节点矩阵，构建配电网络系统的节点关联矩阵，包括：
B＝C+T
其中，C表示邻接矩阵，T表示电源节点矩阵；
所述邻接矩阵C中的元素cij的取值：若节点i与节点j相连则取值为1，否则为0，对角线元素规定为1；所述电源节点矩阵T中元素tii的取值：若节点是电源节点，则取值为1，否则为0；
所述节点关联矩阵B中的元素bij的取值：若节点i与节点j相连则取值为1，否则为0；对于对角线元素，若节点是电源节点，则取值为2，否则取值为1。
3.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，以配电网络中配置PMU费用最少构建目标函数，具体为：
其中，n表示配电网络系统中节点的总数；ci为第i个节点的权重系数，表示在第i个节点处安装PMU的费用；ki表示第i个节点配置PMU的状态变量，其值为0时，表示第i个节点处没有配置PMU，其值为1时，表示第i个节点处有配置PMU。
4.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，以可观测节点的数量最大为优化目标构建目标函数，具体为：
其中，gi表示第i个节点的可观测性，如下式所示：
G＝B·K
其中，K＝[k1，k2，…，kn]T为n维列向量，ki表示第i个节点配置PMU的状态变量，其值为0时，表示第i个节点处没有配置PMU，其值为1时，表示第i个节点处有配置PMU；B为n×n阶的节点关联矩阵；G＝[g1，g2，…，gn]T为n维列向量，gi表示第i个节点可观测性的状态变量，其值为0时，表示第i个节点不可以观测，其值为1时，表示第i个节点可以观测。
5.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，所述对两种情况下构建的目标函数进行综合得到总的目标函数，表示为：
max F＝F2-F1。
6.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，所述构建节点关联度公式，具体为：
其中，D1(i)表示节点i的关联度；bij为节点关联矩阵中的元素，表示节点i与节点j的连接情况。
7.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，所述构建节点凝聚度公式，具体包括：
将配电网络表示为图G(V,E)，其中G为无向网络连通图，其中有n个节点，m条支路，V表示节点集合，E表示支路的集合；首先定义图G的凝聚度，为：
其中，
表示网络中所有节点对最小距离的平均值；以Gi表示为节点i收缩后的图，则定义节点i的凝聚度公式，为：
其中，D2(Gi)表示节点i收缩后配电网络的凝聚度。
8.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，所述将节点关联度公式、节点凝聚度公式结合起来，构建节点的关键度函数，表示为：
D(i)＝W1×D1(i)+W2×D2(i)
其中，D(i)表示节点i的关键度，节点关键度的值越大，就证明该节点在网络中处于关键、核心的位置，应该在此节点处配置PMU；W1、W2分别为权重系数，由熵值法来确定。
9.根据权利要求1所述的基于节点关联矩阵的配电网络PMU配置优化方法，其特征在于，所述采用群体智能优化算法求解所述配电网络PMU配置的多目标优化模型，得到配电网络PMU配置的最优方案，包括：
根据配电网络的拓扑数据，依次确定关联矩阵、邻接矩阵、电源节点矩阵以及节点关联矩阵后，依据节点关联矩阵，将电源节点作为配置PMU优先级最高的节点；然后根据关联度公式、凝聚度公式，求取配电网络中所有节点的关联度和凝聚度，形成数据矩阵，通过熵值法求取关联度、凝聚度的权重系数，再通过关键度函数计算每个节点的关键度；
将节点配置PMU的优先级按照节点的关键度从大到小进行排列，对于关键度相同的节点，利用群体智能优化算法进行排列组合，这样就会得到所有节点的多组配置PMU优先级不同的组合；
将所述组合的PMU安装位置以及各PMU的安装成本导入所述总的目标函数中，利用群体智能优化算法求解出这些组合的目标函数值，则求解出的目标函数值最大的组合即为配电网络PMU配置的最优方案。
10.一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，处理器执行计算机程序时实现根据权利要求1至9中任一权利要求所述的基于节点关联矩阵的配电网络PMU配置优化方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110871520.XA CN113591258A (zh) | 2021-07-30 | 2021-07-30 | 基于节点关联矩阵的配电网络pmu配置优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110871520.XA CN113591258A (zh) | 2021-07-30 | 2021-07-30 | 基于节点关联矩阵的配电网络pmu配置优化方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113591258A true CN113591258A (zh) | 2021-11-02 |
Family
ID=78252660
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110871520.XA Pending CN113591258A (zh) | 2021-07-30 | 2021-07-30 | 基于节点关联矩阵的配电网络pmu配置优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113591258A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115906353A (zh) * | 2022-11-17 | 2023-04-04 | 国网上海市电力公司 | 一种基于节点评估的配电网pmu优化配置方法 |
-
2021
- 2021-07-30 CN CN202110871520.XA patent/CN113591258A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115906353A (zh) * | 2022-11-17 | 2023-04-04 | 国网上海市电力公司 | 一种基于节点评估的配电网pmu优化配置方法 |
| CN115906353B (zh) * | 2022-11-17 | 2023-08-08 | 国网上海市电力公司 | 一种基于节点评估的配电网pmu优化配置方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hua et al. | Optimal energy management strategies for energy Internet via deep reinforcement learning approach | |
| Müller et al. | Genetic algorithm‐based phasor measurement unit placement method considering observability and security criteria | |
| Xu et al. | Multiagent-based reinforcement learning for optimal reactive power dispatch | |
| Vinothkumar et al. | Distributed generation planning: A new approach based on goal programming | |
| CN104037776A (zh) | 随机惯性因子粒子群优化算法的电网无功容量配置方法 | |
| US20150084432A1 (en) | Power System Management Device and Method | |
| Theodorakatos | A nonlinear well-determined model for power system observability using Interior-Point methods | |
| Nawaz et al. | Distributed stochastic security constrained unit commitment for coordinated operation of transmission and distribution system | |
| CN114915546B (zh) | 一种停电用户定位方法、装置、设备、介质 | |
| CN110299762A (zh) | 基于pmu准实时数据的主动配电网抗差估计方法 | |
| AU2022231757B2 (en) | Method for multi-adaptive optimal μPMU placement in micro-energy network | |
| CN113591258A (zh) | 基于节点关联矩阵的配电网络pmu配置优化方法 | |
| Wang et al. | Parameter identification for a power distribution network based on MCMC algorithm | |
| Khuntia et al. | Risk‐based security assessment of transmission line overloading considering spatio‐temporal dependence of load and wind power using vine copula | |
| Pappala et al. | Predictive optimal control of wind farm reactive sources | |
| Zhang et al. | Multi–objective Cluster Partition Method for Distribution Network Considering Uncertainties of Distributed Generations and Loads | |
| Kumar et al. | An artificial bee colony algorithm based distribution system state estimation including renewable energy sources | |
| Paensuwan et al. | Improved risk-based TTC evaluation with system case partitioning | |
| CN115513945A (zh) | 一种基于物联网的分布式电源管理方法 | |
| Laouid et al. | Optimal PMUs Placement to ensure Power System Observability under Various Contingencies. | |
| Iwata et al. | Multi-population differential evolutionary particle swarm optimization for distribution state estimation using correntropy in electric power systems | |
| Xue et al. | Typical transmission section searching method considering geographical attributes for large power grids | |
| Hou et al. | A novel algorithm for multi-node load forecasting based on big data of distribution network | |
| CN113449908A (zh) | 考虑需求侧资源的配电网规划方法及系统 | |
| Lakhera et al. | Approximating electrical distribution networks via mixed-integer nonlinear programming |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication |