Disclosure of Invention
The invention aims to provide a method for searching an optimal active splitting section of a power system, which can quickly and accurately determine the active splitting section and simultaneously take the two problems of unbalanced power of an island after splitting and tidal current impact on the system caused by splitting operation into consideration.
In order to achieve the above purpose, the solution of the invention is:
an active splitting optimal section searching method for a power system comprises the following steps:
step 1, establishing a system authorized undirected graph according to the running state information of the power system;
step 2, calculating the electrical connection from each load node to each generator group by using a Floyd algorithm for obtaining the shortest distance in the graph theory;
step 3, dividing the load nodes into general nodes and public nodes according to the difference of the electrical contact degree between each load node and different generator groups;
and 4, solving an optimal splitting section in the common node area by taking the minimum unbalanced power of the split system island as a target.
In the step 1, the ratio X of the reactance value of the line to the transmitted active power is taken into consideration of the energy transfer size and distance ij /P ij And measuring the tightness degree of the electrical connection between the nodes, wherein the smaller the value of the tightness degree is, namely the smaller the reactance is, the larger the transmitted active power is, the tighter the electrical connection between the adjacent nodes is, and further establishing a system weighted undirected graph adjacency matrix.
The details of the step 2 are as follows: for the power system weighted undirected graph, when the edge weight is set as the ratio of the line reactance to the transmission power, the average minimum edge weight summation of the load node to each generator in the generator group is obtained through the Floyd algorithm to measure the electrical connection strength of the load node to the generator group, as shown in the following formula:
in the formula, N m Representing the number of generators in the cluster m, W ij And representing the minimum edge weight accumulated sum between the node i and the generator node j in the cluster m.
In step 3, it is assumed that m is present for both clusters 1 And m 2 Let us orderWherein, the first and the second end of the pipe are connected with each other, are respectively a cluster m 1 、m 2 Strength of electrical connection of; taking a threshold value of epsilon, if gamma&Epsilon, the node i belongs to a public node with similar relation degree with the two groups; gamma ray to&And e, i belongs to a common node, i is classified into a cluster which is closely connected with the common node by comparing the positive and negative of gamma, wherein the threshold value e is specifically set according to statistical data and operation experience, and the common node is about 20 percent of the load node.
The value of epsilon is such that the common node is 20% of the load nodes.
In the step 4, the splitting section is obtained in the common node area by taking the minimum unbalanced power as a target: endowing the generator node with a weight P G The load node weight is-P L And sequentially searching common nodes from a certain cluster and a general node set area contained in the cluster, carrying out objective function check on the common nodes, if unbalanced power in the area is reduced, bringing the common nodes into the area until the common nodes meeting the conditions do not exist, and obtaining a final splitting section.
After adopting the scheme, compared with the prior art, the invention has the beneficial effects that:
the search space of the splitting section is reduced based on the electrical connection, so that the searching speed of the splitting section is improved. The method accurately measures the electrical connection between the nodes according to the ratio of the line reactance to the transmission power, does not cause large tidal current impact because the active power is introduced when the nodes are divided according to the electrical connection, takes the two problems of the unbalanced power of the island after splitting and the tidal current impact caused by the splitting operation to the system into consideration, and can better adapt to the requirement of the power system on the splitting scheme.
Detailed Description
The technical solution and the advantages of the present invention will be described in detail with reference to the accompanying drawings.
As shown in fig. 1, the present invention provides a method for searching an optimal active splitting section of a power system, which comprises the following steps:
step 1, establishing a system authorized undirected graph according to the running state information of the power system;
considering the energy transfer size and distance, and taking the ratio X of the line reactance value to the transmitted active power size ij /P ij And measuring the tightness degree of the electrical connection between the nodes, wherein the smaller the value of the tightness degree is, namely the smaller the reactance is, the larger the transmitted active power is, the tighter the electrical connection between the adjacent nodes is, and further establishing a system weighted undirected graph adjacency matrix.
Constructing an adjacent matrix of a weighted undirected graph of the power system, and if any two nodes i and j in the system are adjacent, forming an edge of the undirected graph by a line connecting the two nodes, wherein an edge weight is set as a ratio of a line reactance to active power transmitted between i and j; if the two nodes are not adjacent, the corresponding value in the adjacent matrix is infinite; if two nodes actually represent the same node, the value in the adjacency matrix takes 0.
Step 2, calculating the electrical connection from each load node to each generator group by using a Floyd algorithm for obtaining the shortest distance in the graph theory;
the Floyd algorithm is a classical method for solving the shortest distance between two nodes in graph theory, and the core idea is to select an arc section with the smallest distance from unmarked points, further perform cyclic comparison and finally determine the shortest distance between the two nodes. For the power system weighted undirected graph, when the edge weight is set as the ratio of the line reactance to the transmission power, the minimum edge weight sum between two nodes obtained by the Floyd algorithm reflects the electrical connection strength between the two nodes to a certain extent. And then, the average minimum edge weight summation of the load node to each generator in the generator group can be obtained to measure the electrical connection strength from the load node to the generator group, as shown in formula (2):
in the formula, N m Representing the number of generators in the cluster m, W ij And representing the minimum edge weight accumulated sum between the node i and the generator node j in the cluster m.
Step 3, dividing the load nodes into general nodes and common nodes according to the difference of the electrical contact degree between each load node and different generator groups;
the load nodes in the system after the instability can be divided into two types, one type is a general node which has strong electrical connection with a certain cluster, and the other type is a common node which has similar electrical connection degree with two or more clusters. The majority of common nodes are divided into clusters closely connected with the common nodes, and the searching and splitting section substantially determines the splitting scheme to reasonably divide the common nodes into the corresponding clusters, thereby effectively reducing the searching space.
Suppose that for two clusters m 1 And m 2 Let us orderTaking a threshold value of epsilon, if gamma&If the node i belongs to a public node with similar relation degree with the two computer groups, the node i is a node with similar relation degree with the two computer groups; gamma ray&And e, i belongs to a common node, i is classified into a cluster which is closely connected with the common node by comparing the positive and negative of gamma, wherein the threshold value e is specifically set according to statistical data and operation experience, and the common node is about 20 percent of the load node.
And 4, solving an optimal splitting section in the common node area by taking the minimum unbalanced power of the split system island as a target.
The objective functions of the splitting section search mainly include two types, one is minimum unbalanced power between the generator and the load in the split region, and the other is minimum tidal current impact, and expressions are respectively shown in formulas (3) and (4):
min|P G -P L | (3)
in the formula, P G 、P L Respectively representing generator power and load power in the split region, P ij Representing the power transmission size between node i and node j, wherein i and j belong to different areas V 1 、V 2 . The minimum unbalanced power ensures less generator tripping and load shedding of the system after the splitting, and the minimum tide impact aims to reduce the impact on the system caused by the splitting operation. In the simplified search space, namely, the common node area, the invention obtains the splitting section by taking the minimum unbalanced power as a target: endowing the generator node with a weight P G The load node weight is-P L And sequentially searching common nodes from a certain cluster and a general node set area contained in the cluster, carrying out objective function check on the common nodes, if unbalanced power in the area is reduced, bringing the common nodes into the area until the common nodes meeting the conditions do not exist, and obtaining a final splitting section.
Taking the new england 39-node standard system as an example for testing, the system structure is shown in fig. 2, and the system is a standard system for analyzing the transient safety of the system, and comprises 10 generators and 46 power lines.
A fault condition: the 0s three-phase short circuit at the line 23-24, the 0.4s fault disappears, the generator under the fault is divided into two groups, wherein the generators 33-36 are the advanced group.
Firstly, an adjacency matrix is calculated according to line reactance and power flow, a system weighted undirected graph is established, the average minimum edge weight cumulative sum from a load node to two clusters is calculated respectively by adopting a Floyd algorithm, and the difference between the edge weight cumulative sum from each load node to the two clusters is obtained, as shown in Table 1.
TABLE 1 difference between the cumulative sums of the minimum edge weights
Taking the threshold e =0.035, the common nodes 2, 3, 10, 11, 12, 13, 14, 25 are available, while the clusters 33-36 comprise the common nodes 15, 16, 17, 18, 19, 20, 21, 22, 23, 24. The splitting of the common nodes into respective clusters with the aim of minimum unbalanced power, the resulting splitting pattern is shown in fig. 3.
In the splitting mode, the minimum unbalanced power of the two islands is 0.329p.u and 0.130p.u respectively, and the transmission power of a line is considered when the electrical connection between the nodes is measured, so that the splitting scheme can effectively avoid the line with larger power flow, and the condition that the system is subjected to larger power flow impact due to the splitting operation is avoided, and the power flow impact under the scheme is 0.805p.u. And because the invention reduces the search space of the column-searching section based on the electrical connection, the solving complexity is greatly reduced, the time has obvious advantages, the solving time in the process is about 11ms, and the invention can better adapt to the requirement of the actual system on the response rate.
In summary, the invention discloses a method for searching the optimal splitting section of an active splitting of a power system, so as to quickly and accurately search the optimal splitting section after the system is unstable. Firstly, measuring the electrical connection from a load node to a generator group according to an electrical connection index between nodes fusing line reactance and transmission power, and dividing the load node into a general node and a common node so as to reduce a search space of a splitting section; and then, quickly determining an optimal splitting section by taking the minimum unbalanced power of the split system island as a target. Compared with the existing splitting section searching method, the method provided by the invention has higher efficiency, can take two problems of unbalanced power of an island after splitting and tidal current impact on a system caused by splitting operation into consideration, and can better adapt to the requirement of a power system on a splitting scheme.
The above embodiments are only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited thereby, and any modifications made on the basis of the technical scheme according to the technical idea of the present invention fall within the protection scope of the present invention.