CN111222719A - Power grid region division method and device - Google Patents

Power grid region division method and device Download PDF

Info

Publication number
CN111222719A
CN111222719A CN202010123708.1A CN202010123708A CN111222719A CN 111222719 A CN111222719 A CN 111222719A CN 202010123708 A CN202010123708 A CN 202010123708A CN 111222719 A CN111222719 A CN 111222719A
Authority
CN
China
Prior art keywords
node
electrical distance
voltage
power grid
cluster center
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.)
Withdrawn
Application number
CN202010123708.1A
Other languages
Chinese (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.)
Electric Power Research Institute of Yunnan Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Yunnan Power Grid Co Ltd
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 Electric Power Research Institute of Yunnan Power Grid Co Ltd filed Critical Electric Power Research Institute of Yunnan Power Grid Co Ltd
Priority to CN202010123708.1A priority Critical patent/CN111222719A/en
Publication of CN111222719A publication Critical patent/CN111222719A/en
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • 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"
    • G06Q10/043—Optimisation of two dimensional placement, e.g. cutting of clothes or wood
    • 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • 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

Landscapes

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

Abstract

The application discloses a power grid region division method and a device, wherein the method comprises the following steps: acquiring a topological structure of a power grid; determining nodes in the power grid according to the topological structure; wherein, the node includes: a voltage pivot and nodes other than the voltage pivot; constructing an admittance matrix according to the nodes; determining a first electrical distance between each node and a voltage center point according to the admittance matrix and the node parameters; dividing the power grid into a plurality of areas according to the first electrical distance; determining the cluster center of each region according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center; and calculating a second electrical distance between the node and the cluster center, and adjusting the node in each area according to the second electrical distance and the minimum variance in each area to finish the area division of the power grid. Therefore, the electric characteristics among the nodes are utilized, the region division is carried out on the power grid, and the application range is wider.

Description

Power grid region division method and device
Technical Field
The application relates to the technical field of power grid management, in particular to a power grid region division method and device.
Background
Along with the continuous rising of power demand and the rising of the permeability of distributed energy at the user side, the scale of the power system is continuously enlarged, and all levels of power grids in the power system are mutually connected, so that the structure of the power grid in the power system is more complex. The complicated power grid structure increases the control difficulty of the power system, namely, the operation analysis, the dispatching management and the load management difficulty of the power system are gradually increased. In order to better realize the control of the power system, regional management of a power grid is introduced. Therefore, regional management of the power grid is of great significance.
In grid-based regional management, a node having a strong coupling property is generally used as the same region, and a weak coupling node is divided into different regions, so as to improve control of a power system. At present, a part of power grid structures adopt power grid areas divided according to regions, but the power grid areas divided according to the regions only consider the natural characteristics of the power grid structures and ignore the electrical characteristics, so that the application of the area division has limitation.
Disclosure of Invention
The embodiment of the application discloses a power grid region division method and device, and aims to solve the problem that in the prior art, only the natural characteristic of a power grid structure is considered in the power grid region division according to regions, the electrical characteristic is ignored, and the application of the region division is limited.
In a first aspect, an embodiment of the present application provides a power grid area division method, where the method includes:
acquiring a topological structure of a power grid;
determining nodes in the power grid according to the topological structure; wherein the node comprises: a voltage pivot and nodes other than the voltage pivot;
constructing an admittance matrix according to the nodes;
determining a first electrical distance between each node and a voltage center point according to the admittance matrix and the node parameters;
dividing the power grid into a plurality of areas according to the first electrical distance;
determining the cluster center of each region according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center;
and calculating a second electrical distance between the node and the cluster center, and adjusting the node in each area according to the second electrical distance and the minimum variance in each area to finish the area division of the power grid.
With reference to the first aspect, in one implementation, determining a first electrical distance between each node and a voltage pivot point according to the admittance matrix and the node parameters includes:
according to the imaginary part of the admittance matrix and the node voltage, calculating the sensitivity of reactive power to voltage by adopting the following formula;
S＝[-B”]-1/U；
wherein, S is a voltage-to-reactive sensitivity matrix, B' is an admittance matrix element imaginary matrix, and U is voltage;
calculating a first electrical distance from the sensitivity using the following equation:
Figure BDA0002393779430000021
wherein σijIs the electrical distance, σ, from node i to node jjiIs the electrical distance from node j to node i, SijSensitivity of reactive power to voltage, S, for nodes i to jjiThe reactive versus voltage sensitivity of nodes j to i.
With reference to the first aspect, in an implementation manner, the most upstream node is a tidal current power originating node in a topology.
With reference to the first aspect, in an implementation manner, the cluster center of each region is determined according to the initial cluster center and the clustering index, and is calculated by using the following formula:
Figure BDA0002393779430000022
wherein m is a node,
Figure BDA0002393779430000023
is the cluster center of the nth region, when n is 2,
Figure BDA0002393779430000024
is the initial cluster center, when > 2,
Figure BDA0002393779430000025
is the cluster center of the (n-1) th area iterated according to the initial cluster center, and theta% is a clustering index, and the value of the theta% is set according to the precision required by power grid or load management.
With reference to the first aspect, in an implementation manner, adjusting nodes in each region according to the second electrical distance and a minimum variance in each region to complete power grid region division includes:
obtaining a minimum variance according to the sum of squares of errors of the electrical distance and the cluster center in each region;
obtaining an objective function according to the second electrical distance and the minimum variance, performing region division on the power grid according to the objective function, and calculating the power grid region division objective function by adopting the following formula:
Figure BDA0002393779430000026
F(σ,n)＝min∑|J(σ,n)|2；
wherein σi0Representing a first electrical distance, J (σ, n) being the sum of the squares of the errors in each sub-region; f (σ, n) is the objective function.
In a second aspect, an embodiment of the present application provides a power grid area dividing apparatus, where the apparatus includes:
the acquisition module is used for acquiring a topological structure of the power grid;
the node determining module is used for determining nodes in the power grid according to the topological structure; wherein the node comprises: a voltage pivot and nodes other than the voltage pivot;
the admittance matrix construction module is used for constructing an admittance matrix according to the nodes;
the first electrical distance determining module is used for determining a first electrical distance between each node and a voltage pivot point according to the admittance matrix and the node parameters;
the dividing module is used for dividing the power grid into a plurality of areas according to the first electrical distance;
the cluster center determining module is used for determining the cluster center of each area according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center;
and the region division module is used for calculating a second electrical distance between the node and the cluster center, and adjusting the node in each region according to the second electrical distance and the minimum variance in each region to finish the region division of the power grid.
With reference to the second aspect, in one implementation, the first electrical distance determination module includes:
the first calculation unit is used for calculating the sensitivity of reactive power to voltage by adopting the following formula according to the imaginary part of the admittance matrix and the node voltage;
S＝[-B”]-1/U；
wherein, S is a voltage-to-reactive sensitivity matrix, B' is an admittance matrix element imaginary matrix, and U is voltage;
a second calculating unit, configured to calculate a first electrical distance according to the sensitivity, using the following formula:
Figure BDA0002393779430000031
wherein σijIs the electrical distance, σ, from node i to node jjiIs the electrical distance from node j to node i, SijSensitivity of reactive power to voltage, S, for nodes i to jjiThe reactive versus voltage sensitivity of nodes j to i.
With reference to the second aspect, in one implementation manner, the cluster center determining module includes: the third calculation unit is used for calculating by adopting the following formula according to the initial cluster center and the clustering index to obtain the cluster center of each area;
Figure BDA0002393779430000032
wherein m is a node,
Figure BDA0002393779430000033
is the cluster center of the nth region, when n is 2,
Figure BDA0002393779430000034
is the initial cluster center, when > 2,
Figure BDA0002393779430000035
is the cluster center of the (n-1) th area iterated according to the initial cluster center, and theta% is a clustering index, and the value of the theta% is set according to the precision required by power grid or load management.
With reference to the second aspect, in one implementation manner, the region dividing module includes:
the minimum variance acquiring unit is used for acquiring minimum variance according to the sum of the square of the error between the electrical distance in each area and the center of the cluster;
the target function obtaining unit is used for obtaining a target function according to the second electrical distance and the minimum variance;
the region division unit is used for performing region division on the power grid according to the target function;
the power grid region division target function is calculated by adopting the following formula:
Figure BDA0002393779430000036
F(σ,n)＝min∑|J(σ,n)|2；
wherein σi0Representing a first electrical distance, J (σ, n) being the sum of the squares of the errors in each sub-region; f (σ, n) is the objective function.
The embodiment of the application discloses a power grid region division method and a device, wherein the method comprises the following steps: firstly, acquiring a topological structure of a power grid; then determining nodes in the power grid according to the topological structure; wherein the node comprises: a voltage pivot and nodes other than the voltage pivot; constructing an admittance matrix according to the nodes; determining a first electrical distance between each node and a voltage center point according to the admittance matrix and the node parameters; dividing the power grid into a plurality of areas according to the first electrical distance; determining the cluster center of each region according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center; and calculating a second electrical distance between the node and the cluster center, and adjusting the node in each area according to the second electrical distance and the minimum variance in each area to finish the area division of the power grid.
According to the scheme, the power grid is subjected to preliminary regional division by determining the first electrical distance between each node and the voltage pivot point to obtain a plurality of sub-regions, and then further adjustment is performed on the sub-regions according to a hierarchical clustering method and the error square sum in the cluster, so that the regional division of the power grid is more accurate.
Drawings
In order to more clearly explain the technical solution of the present application, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious to those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a power grid region dividing method provided in an embodiment of the present application;
fig. 2 is a schematic structural diagram of a power grid area dividing apparatus provided in an embodiment of the present application.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, the present application is described in further detail with reference to the accompanying drawings and the detailed description.
The embodiment of the application discloses a power grid region division method and device, and aims to solve the problem that in the prior art, only the natural characteristic of a power grid structure is considered in the power grid region division according to regions, the electrical characteristic is ignored, and the application of the region division is limited.
Referring to fig. 1, there is shown a method of grid area partitioning, the method comprising the steps of:
and step S1, acquiring the topological structure of the power grid.
In this step, the topological structure of the power grid to be partitioned is obtained first, so as to partition the power grid region according to the topological structure, and as to how to obtain and the specific form of the topological structure, the existing method and structure may be adopted, and the present application is not limited specifically.
Step S2, determining nodes in the power grid according to the topological structure; wherein the node comprises: a voltage pivot point and nodes other than the voltage pivot point.
The nodes in the power grid are generally classified into PQ nodes and PV nodes. The load node is mostly a PQ node, and the voltage pivot point is a PV node.
The voltage pivot point aims at regulating and controlling the first (last) end voltage of the key line, so that the subsequent node voltage is improved, and the voltage pivot point is selected as a reference node, so that the voltage pivot point has a control advantage which is difficult to compare with.
In this step, the voltage pivot and the nodes other than the voltage pivot are determined according to the topological relation in the topological structure in step S1.
And step S3, constructing an admittance matrix according to the nodes.
In this step, the specific method for constructing the admittance matrix may adopt an existing method, and the present application is not particularly limited.
And step S4, determining a first electric distance between each node and the voltage pivot point according to the admittance matrix and the node parameters.
Wherein the node parameters may include node voltage and power.
The connection impedance between two nodes on a power system is defined as the electrical distance, the physical meaning of which is to indicate whether the electrical connection between the two nodes is tight, i.e. whether the electrical connection between the two nodes is tightThat is, the electrical distance is an important index reflecting the degree of electrical connection between nodes, for example, the electrical distance between node i and node j may be represented by the voltage change amplitude Δ U of node iiThe sensitivity of the power injection variable quantity at the node j and the size of the electrical distance represent the strength of electrical connection between the nodes.
In this embodiment, the node power flow formula matrix form is as follows:
(ΔP/U＝-B'UΔσ)t×t(1)；
(ΔQ/U＝-B”ΔU)t×t(2)；
the delta P is the active variable quantity injected into the load node, and the delta Q is the reactive variable quantity injected into the load node; delta sigma is the voltage phase angle variation of the load node, and delta U is the voltage amplitude variation of the load node; u is an initial voltage; b ', B' are respectively the imaginary parts of the admittance matrix of the corresponding elements; t is the matrix order.
And step S5, dividing the power grid into a plurality of areas according to the first electrical distance.
In this embodiment, the first electrical distance obtained in step S4 is used to perform preliminary region division on the power grid according to the value, so as to obtain a plurality of sub-regions.
Step S6, determining the cluster center of each region according to the initial cluster center and the clustering index; wherein the initial cluster center is an electrical distance between the most upstream node and the voltage center point.
Optionally, the most upstream node is a tidal current power originating node in the topology.
In this step, the electrical distance between the most upstream node and the voltage center is first used as the initial cluster center of the hierarchical clustering method, and then the cluster centers of the regions divided in step S5 are determined according to the initial cluster centers and the clustering indexes, wherein the clustering indexes are set according to the accuracy required by the power grid or load management.
And step S7, calculating a second electrical distance between the node and the cluster center, and adjusting the nodes in each area according to the second electrical distance and the minimum variance in each area to finish the division of the power grid area.
In this step, first, a second electrical distance between each node in step S2 and the cluster center of each region obtained in step S6 is calculated, the second electrical distance is compared with the minimum variance in each region, if the second electrical distance is smaller than the minimum variance, the node is merged into the region, otherwise, no adjustment is made.
And S6-S7 are implemented by adopting a hierarchical clustering method and the error square sum in the cluster, and the region primarily divided in the step S5 is further adjusted, so that the power grid region division is more accurate and the application range is wider.
The embodiment of the application discloses a power grid region division method and a device, wherein the method comprises the following steps: firstly, acquiring a topological structure of a power grid; then determining nodes in the power grid according to the topological structure; wherein the node comprises: a voltage pivot and nodes other than the voltage pivot; constructing an admittance matrix according to the nodes; determining a first electrical distance between each node and a voltage center point according to the admittance matrix and the node parameters; dividing the power grid into a plurality of areas according to the first electrical distance; determining the cluster center of each region according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center; and calculating a second electrical distance between the node and the cluster center, and adjusting the node in each area according to the second electrical distance and the minimum variance in each area to finish the area division of the power grid.
According to the scheme, the power grid is subjected to preliminary regional division by determining the first electrical distance between each node and the voltage pivot point to obtain a plurality of sub-regions, and then further adjustment is performed on the sub-regions according to a hierarchical clustering method and the error square sum in the cluster, so that the regional division of the power grid is more accurate.
Optionally, determining a first electrical distance between each node and the voltage pivot point according to the admittance matrix and the node parameters may be implemented by:
according to the imaginary part of the admittance matrix and the node voltage, calculating the sensitivity of reactive power to voltage by adopting the following formula;
S＝[-B”]-1/U (3)；
wherein, S is a voltage-to-reactive sensitivity matrix, B' is an admittance matrix element imaginary matrix, and U is voltage;
calculating a first electrical distance from the sensitivity using the following equation:
Figure BDA0002393779430000061
wherein σijIs the electrical distance, σ, from node i to node jjiIs the electrical distance from node j to node i, SijSensitivity of reactive power to voltage, S, for nodes i to jjiThe reactive versus voltage sensitivity of nodes j to i.
Optionally, the cluster center of each region is determined according to the initial cluster center and the clustering index, and the following formula is adopted for calculation:
Figure BDA0002393779430000062
wherein m is a node,
Figure BDA0002393779430000063
is the cluster center of the nth region, when n is 2,
Figure BDA0002393779430000064
is the initial cluster center, when > 2,
Figure BDA0002393779430000065
is the cluster center of the (n-1) th area iterated according to the initial cluster center, and theta% is a clustering index, and the value of the theta% is set according to the precision required by power grid or load management.
In this embodiment, the set formed by the first electrical distances is sorted, then the cluster center of each region is calculated, assuming that the number of the voltage pivot point is 0, the electrical distance σ from each load node i to the pivot point can be solved according to the formula (4)i0And obtaining a node electrical distance set a as follows:
a＝[σi0]T(6)；
sorting the set a according to the value size to obtain a new set β as follows:
β＝[bm]T(7)；
wherein the content of the first and second substances,
Figure BDA0002393779430000066
Figure BDA0002393779430000072
and (3) determining the cluster centers of the regions one by one according to a formula (7), and stopping calculation until all the regions are confirmed, and then calculating by adopting a formula (5) according to the sorted set beta to obtain the cluster centers of the regions.
Optionally, the nodes in each region are adjusted according to the second electrical distance and the minimum variance in each region to complete the power grid region division, and the method may be implemented by:
obtaining a minimum variance according to the sum of squares of errors of the electrical distance and the cluster center in each region;
obtaining an objective function according to the second electrical distance and the minimum variance, performing region division on the power grid according to the objective function, and calculating the power grid region division objective function by adopting the following formula:
Figure BDA0002393779430000071
F(σ,n)＝min∑|J(σ,n)|2(9)；
wherein σi0Representing a first electrical distance, J (σ, n) being the sum of the squares of the errors in each sub-region; f (σ, n) is the objective function.
In the embodiment, the nodes except the cluster center are allocated to the most similar region according to the error sum of squares formula in the cluster, and the optimal solution is obtained through cyclic selection and is iterated step by step, so that the division precision is further improved, and the application range is wider.
Referring to fig. 2, there is shown a grid area partitioning apparatus, the apparatus comprising:
the acquisition module 10 is used for acquiring a topological structure of a power grid;
the node determining module 20 is configured to determine a node in the power grid according to the topology; wherein the node comprises: a voltage pivot and nodes other than the voltage pivot;
an admittance matrix construction module 30, configured to construct an admittance matrix according to the nodes;
a first electrical distance determining module 40, configured to determine a first electrical distance between each node and the voltage pivot point according to the admittance matrix and the node parameter;
a dividing module 50, configured to divide the power grid into a plurality of regions according to the first electrical distance;
a cluster center determining module 60, configured to determine a cluster center of each region according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center;
and the region division module 70 is configured to calculate a second electrical distance between the node and the cluster center, and adjust the node in each region according to the second electrical distance and the minimum variance in each region to complete power grid region division.
Optionally, the first electrical distance determination module comprises:
the first calculation unit is used for calculating the sensitivity of reactive power to voltage by adopting the following formula according to the imaginary part of the admittance matrix and the node voltage;
S＝[-B”]-1/U (10)；
wherein, S is a voltage-to-reactive sensitivity matrix, B' is an admittance matrix element imaginary matrix, and U is voltage;
a second calculating unit, configured to calculate a first electrical distance according to the sensitivity, using the following formula:
Figure BDA0002393779430000081
wherein σijIs the electrical distance, σ, from node i to node jjiIs the electrical distance from node j to node i, SijSensitivity of reactive power to voltage, S, for nodes i to jjiThe reactive versus voltage sensitivity of nodes j to i.
Optionally, the cluster center determining module includes: the third calculation unit is used for calculating by adopting the following formula according to the initial cluster center and the clustering index to obtain the cluster center of each area;
Figure BDA0002393779430000082
wherein m is a node,
Figure BDA0002393779430000083
is the cluster center of the nth region, when n is 2,
Figure BDA0002393779430000084
is the initial cluster center, when > 2,
Figure BDA0002393779430000085
is the cluster center of the (n-1) th area iterated according to the initial cluster center, and theta% is a clustering index, and the value of the theta% is set according to the precision required by power grid or load management.
Optionally, the area dividing module includes:
the minimum variance acquiring unit is used for acquiring minimum variance according to the sum of the square of the error between the electrical distance in each area and the center of the cluster;
the target function obtaining unit is used for obtaining a target function according to the second electrical distance and the minimum variance;
the region division unit is used for performing region division on the power grid according to the target function;
the power grid region division target function is calculated by adopting the following formula:
Figure BDA0002393779430000086
F(σ,n)＝min∑|J(σ,n)|2(14)；
wherein σi0Representing a first electrical distance, J (σ, n) being the sum of the squares of the errors in each sub-region; f (σ, n) is the objective function.
The method and the device for dividing the power grid area can realize the area division management of the power system, and are mainly used in the technical fields of power system analysis, scheduling management, power quality supervision, load management and the like.
The same and similar parts in the various embodiments in this specification may be referred to each other. In particular, as for the apparatus embodiment, since it is substantially similar to the method embodiment, the description is simple, and the relevant points can be referred to the description in the method embodiment.
The present application has been described in detail with reference to specific embodiments and illustrative examples, but the description is not intended to limit the application. Those skilled in the art will appreciate that various equivalent substitutions, modifications or improvements may be made to the presently disclosed embodiments and implementations thereof without departing from the spirit and scope of the present disclosure, and these fall within the scope of the present disclosure. The protection scope of this application is subject to the appended claims.

Claims (9)

1. A power grid area division method is characterized by comprising the following steps:
acquiring a topological structure of a power grid;
determining nodes in the power grid according to the topological structure; wherein the node comprises: a voltage pivot and nodes other than the voltage pivot;
constructing an admittance matrix according to the nodes;
determining a first electrical distance between each node and a voltage center point according to the admittance matrix and the node parameters;
dividing the power grid into a plurality of areas according to the first electrical distance;
determining the cluster center of each region according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center;
and calculating a second electrical distance between the node and the cluster center, and adjusting the node in each area according to the second electrical distance and the minimum variance in each area to finish the area division of the power grid.
2. The method of claim 1, wherein determining a first electrical distance between each node and a voltage center point based on the admittance matrix and node parameters comprises:
according to the imaginary part of the admittance matrix and the node voltage, calculating the sensitivity of reactive power to voltage by adopting the following formula;
S＝[-B”]-1/U；
wherein, S is a voltage-to-reactive sensitivity matrix, B' is an admittance matrix element imaginary matrix, and U is voltage;
calculating a first electrical distance from the sensitivity using the following equation:
Figure FDA0002393779420000011
wherein σijIs the electrical distance, σ, from node i to node jjiIs the electrical distance from node j to node i, SijSensitivity of reactive power to voltage, S, for nodes i to jjiThe reactive versus voltage sensitivity of nodes j to i.
3. Method according to claim 1 or 2, wherein the most upstream node is a tidal power originating node in the topology.
4. The method of claim 1, wherein the determining the cluster center of each region according to the initial cluster center and the clustering index is calculated by using the following formula:
Figure FDA0002393779420000012
wherein m is a node,
Figure FDA0002393779420000013
is the cluster center of the nth region, when n is 2,
Figure FDA0002393779420000014
is the initial cluster center, when > 2,
Figure FDA0002393779420000015
is the cluster center of the (n-1) th area iterated according to the initial cluster center, and theta% is a clustering index, and the value of the theta% is set according to the precision required by power grid or load management.
5. The method of claim 4, wherein adjusting nodes in each zone according to the second electrical distance and the minimum variance in each zone to complete grid zone division comprises:
obtaining a minimum variance according to the sum of squares of errors of the electrical distance and the cluster center in each region;
obtaining an objective function according to the second electrical distance and the minimum variance, performing region division on the power grid according to the objective function, and calculating the power grid region division objective function by adopting the following formula:
Figure FDA0002393779420000021
F(σ,n)＝min∑|J(σ,n)|2；
wherein σi0Representing a first electrical distance, J (σ, n) being the sum of the squares of the errors in each sub-region; f (σ, n) is the objective function.
6. A grid area partitioning apparatus, comprising:
the acquisition module is used for acquiring a topological structure of the power grid;
the node determining module is used for determining nodes in the power grid according to the topological structure; wherein the node comprises: a voltage pivot and nodes other than the voltage pivot;
the admittance matrix construction module is used for constructing an admittance matrix according to the nodes;
the first electrical distance determining module is used for determining a first electrical distance between each node and a voltage pivot point according to the admittance matrix and the node parameters;
the dividing module is used for dividing the power grid into a plurality of areas according to the first electrical distance;
the cluster center determining module is used for determining the cluster center of each area according to the initial cluster center and the clustering index; wherein the initial cluster center is the electrical distance between the most upstream node and the voltage center;
and the region division module is used for calculating a second electrical distance between the node and the cluster center, and adjusting the node in each region according to the second electrical distance and the minimum variance in each region to finish the region division of the power grid.
7. The apparatus of claim 6, wherein the first electrical distance determination module comprises:
the first calculation unit is used for calculating the sensitivity of reactive power to voltage by adopting the following formula according to the imaginary part of the admittance matrix and the node voltage;
S＝[-B”]-1/U；
wherein, S is a voltage-to-reactive sensitivity matrix, B' is an admittance matrix element imaginary matrix, and U is voltage;
a second calculating unit, configured to calculate a first electrical distance according to the sensitivity, using the following formula:
Figure FDA0002393779420000022
wherein σijIs the electrical distance, σ, from node i to node jjiIs the electrical distance from node j to node i, SijSensitivity of reactive power to voltage, S, for nodes i to jjiThe reactive versus voltage sensitivity of nodes j to i.
8. The apparatus of claim 6, wherein the cluster center determining module comprises: the third calculation unit is used for calculating by adopting the following formula according to the initial cluster center and the clustering index to obtain the cluster center of each area;
Figure FDA0002393779420000023
wherein m is a node,
Figure FDA0002393779420000031
is the cluster center of the nth region, when n is 2,
Figure FDA0002393779420000032
is the initial cluster center, when > 2,
Figure FDA0002393779420000033
is the cluster center of the (n-1) th area iterated according to the initial cluster center, and theta% is a clustering index, and the value of the theta% is set according to the precision required by power grid or load management.
9. The apparatus of claim 8, wherein the region partitioning module comprises:
the minimum variance acquiring unit is used for acquiring minimum variance according to the sum of the square of the error between the electrical distance in each area and the center of the cluster;
the target function obtaining unit is used for obtaining a target function according to the second electrical distance and the minimum variance;
the region division unit is used for performing region division on the power grid according to the target function;
the power grid region division target function is calculated by adopting the following formula:
Figure FDA0002393779420000034
F(σ,n)＝min∑|J(σ,n)|2；
wherein σi0Representing a first electrical distance, J (σ, n) being the sum of the squares of the errors in each sub-region; f (σ, n) is the objective function.
CN202010123708.1A 2020-02-27 2020-02-27 Power grid region division method and device Withdrawn CN111222719A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010123708.1A CN111222719A (en) 2020-02-27 2020-02-27 Power grid region division method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010123708.1A CN111222719A (en) 2020-02-27 2020-02-27 Power grid region division method and device

Publications (1)

Publication Number Publication Date
CN111222719A true CN111222719A (en) 2020-06-02

Family

ID=70828487

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010123708.1A Withdrawn CN111222719A (en) 2020-02-27 2020-02-27 Power grid region division method and device

Country Status (1)

Country Link
CN (1) CN111222719A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111985576A (en) * 2020-09-02 2020-11-24 南宁师范大学 Shop address selection method based on decision tree
CN113078653A (en) * 2021-04-25 2021-07-06 广东电网有限责任公司湛江供电局 Power distribution network partitioning method based on electrical distance and K-means clustering algorithm

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111985576A (en) * 2020-09-02 2020-11-24 南宁师范大学 Shop address selection method based on decision tree
CN111985576B (en) * 2020-09-02 2023-11-03 南宁师范大学 Shop site selection method based on decision tree
CN113078653A (en) * 2021-04-25 2021-07-06 广东电网有限责任公司湛江供电局 Power distribution network partitioning method based on electrical distance and K-means clustering algorithm

Similar Documents

Publication Publication Date Title
CN109272176B (en) Calculation method is predicted to platform area line loss per unit using K-means clustering algorithm
CN111222719A (en) Power grid region division method and device
CN107529644B (en) Linear approximation method for static voltage stability domain boundary of power system
CN106548418A (en) Power system small interference stability appraisal procedure
CN107069716A (en) A kind of robust real-time scheduling method counted and inject transfer distribution factor evaluated error
CN110210690B (en) Optimal configuration method for micro synchronous phasor measurement unit of power distribution system
CN105447082A (en) Distributed clustering method for mass load curves
CN110601174A (en) Load modeling and online correction method based on deep learning
CN111969604A (en) Dynamic calculation method and device for line loss of transformer area based on measured data
CN107025609A (en) Electrical Power System Dynamic method for estimating state based on singular value decomposition CDKF
CN114386700A (en) Power transmission system planning method and device considering wind-light correlation and storage medium
CN109921426A (en) Wind-electricity integration system probability load flow calculation method based on CV-KDE
CN113468792A (en) Parameter correction method and device of electromagnetic transient simulation model and electronic equipment
CN110880757B (en) Optimal partition-based distributed interval state estimation method for active power distribution network
CN111244944A (en) AC/DC power distribution network reliability assessment method based on influence increment
CN105552885A (en) Method and system for improving power distribution network state estimation observability
CN109698505B (en) Regulation and control quantitative mapping calculation method for large power grid static voltage stability online prevention and control
CN107276093B (en) The Probabilistic Load calculation method cut down based on scene
CN113078653A (en) Power distribution network partitioning method based on electrical distance and K-means clustering algorithm
CN110571788A (en) static voltage stability domain boundary coefficient calculation method based on dynamic equivalent circuit
CN104993503A (en) Island microgrid frequency control method
CN101995845A (en) Field programmable gate array (FPGA)-based auto-leveling control system and method
CN110289646B (en) Intelligent soft switch local control strategy optimization method based on meta-model
CN110675276A (en) Method and system for inversion droop control of direct current power transmission system
CN113690930A (en) NSGA-III algorithm-based medium and long term locating and sizing method for distributed photovoltaic power supply

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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20200602