Static Equivalence Method of Power Grid Based on Genetic Algorithm

: When the electromagnetic transient simulation of a large-scale power grid is carried out, because the simulation is limited by the scale of simulation software, it is necessary to divide the grid into internal and external grids, and the external grid is modeled with equivalent simpliﬁcation. Aiming at resolving the difﬁculties of the traditional grid equivalence method, such as a cumbersome calculation process and harsh calculation conditions, a static grid equivalence method based on the genetic algorithm is proposed in this paper. The method ﬁrst constructs an external grid equivalence network, which includes coupling branches between boundary nodes, and each boundary node is connected to the external grid equivalence power supply through the branches. Then, the external grid equivalence model is used to establish an optimization model for solving the equivalence network parameters based on the information of the internal network, the difference of the power input from the external grid to the boundary nodes before and after equivalence is used as the objective function, and the genetic algorithm is used to realize the iterative optimization of the objective function to obtain a set of optimal post-equivalence network model parameters so that the state variables of the internal power grid before and after the equivalence are matched. The accuracy and effectiveness of the proposed method are veriﬁed through simulation with the CERPI36V7 node system.


Introduction
With the construction of ultra-high-voltage lines and the gradual interconnection of power systems [1], the scale of China's power grid, containing thousands of nodes, continues to expand.Limited by the computing power of computers and the computing scale of simulation software, it will become more and more difficult to directly build the entire power system in simulation software and perform simulation calculations.Therefore, in actual engineering calculations, only the detailed modeling of the research area is needed, while the rest of the grid can be simplified equivalently [2,3].However, in the context of the electricity market, each power company in China belongs to different interests that are both cooperating and competing with each other.Each power operator is often reluctant to disclose complete and accurate information.Due to the inconsistency of information collected by different local power grids, it is difficult to obtain information transmitted from adjacent power grids in real time [4].Under this condition of a lack of external grid information, it is of great significance to study the static equivalent method of power grids based on internal network information only.
At present, the static equivalent method of power grids can be divided into two categories according to the acquisition of external network information.
(1) The static equivalent method of power grids with known external information.The research of this kind of method is relatively mature.Ward equivalence and its improved method [5], Thevenin equivalence and its improved method [6], and the radial equivalence independent method and its improved method [7] can be used to make a more accurate equivalence of the external network, but it cannot solve the problem of unknown or partially known equivalence of the external network information.Greater research is required to determine how to ensure the accuracy of the power grid equivalent based on only one state of the internal power grid when the state of the internal power grid changes greatly.
(2) The static equivalent method of power grids with unknown or partially known external network information.These methods are based on the internal network information, and the equivalent parameters of the external power grid model are estimated by various methods.Reference [8] proposed a single-port network equivalence method based on internal network information.The external power grid is equivalent to a single-port Thevenin equivalent model with a voltage source in series with an impedance, and the equivalent parameters are solved by using the multi-period measurement voltage and current of the boundary nodes.However, this method is only applicable to the equivalence of single-port networks.References [9,10] proposed a multi-port network equivalent method based on internal network information which utilizes the measurement data before and after the internal network branch is disconnected to solve the equivalent parameters.This method requires harsh equivalent conditions and is difficult to apply to the actual system.In order to avoid a series of problems caused by unplanned disconnection of the internal network, reference [11] proposed to simulate the disconnection of the internal network branch by node power perturbation so as to obtain effective supplementary measurement information, which is advantageous in engineering applications.Reference [12] proposed to obtain supplementary measurement information based on historical data instead of internal network branch breaking, which improves the practicability of the method.However, the measurement data at different times lack explicit correlation, and the reliability of the network equivalent parameters obtained from the solution needs to be verified.Reference [13] proposed the Thevenin equivalent impedance calculation method with stratified estimation, which stratifies the samples according to the load-side parameters and transforms the equivalent impedance calculation problem into a statistical problem.This significantly reduces the computational error, but the workload of the equivalent calculation is large when the number of nodes at the grid boundary is large.Although many achievements have been made in the static equivalent method of power grids based on internal network information, there are still many limitations and complicated calculation processes when using the existing methods to solve the equivalent parameters.The optimization algorithm is a global optimization random search algorithm.In principle, it can find the global optimal solution of the equivalent parameters with a large probability, which can provide a new idea for the static equivalent method of power grids.
Therefore, this paper proposes a static equivalent method of power grids based on the genetic algorithm.The external power grid equivalent network of this method includes the coupling branches between the boundary nodes, which can effectively reflect the electrical coupling relationship between the external power grids.This method is based on the measurement data of the internal network at one moment, which avoids the impact on the grid caused by an interruption or disturbance on the power grid.The algorithm uses the difference in the power input from the external grid to the boundary node before and after equipartitioning as the objective function.The genetic algorithm is used to solve the equivalent parameters of the equivalent external power grid model, which solves the problems such as harsh requirements for solving the equivalent parameters and a complicated solving process.The correctness and practicability of the method proposed in this paper are verified by the CERPI36v7 system.
The rest of this paper is organized as follows.In Section 2, we present a brief description of the structure and parameters of the equivalent power grid model.In the Section 3, the optimization model of power grid parameters after equivalence and the process of solving equivalent parameters based on genetic algorithm are introduced.Later, results and discussion are provided in Section 4. Finally, Section 5 concludes the paper.

Static Equivalent Network Model of Power Grids
According to the needs of the study, the whole system is divided into the internal power grid, the boundary node, and the external power grid shown in Figure 1.The internal power grid is the part of the power grid that needs to be studied and analyzed.It is not necessary to model other regions in detail, and this part is referred to as the external power grid.The external power grid is connected to each boundary node in the form of an equivalent branch of voltage source series impedance, and the electrical coupling relationship between the external power grids is reflected by establishing the coupling branch between the boundary nodes [14].The equivalent network is shown in Figure 2.
The rest of this paper is organized as follows.In Section 2, we present a brief description of the structure and parameters of the equivalent power grid model.In the Section 3, the optimization model of power grid parameters after equivalence and the process of solving equivalent parameters based on genetic algorithm are introduced.Later, results and discussion are provided in Section 4. Finally, Section 5 concludes the paper.

Static Equivalent Network Model of Power Grids
According to the needs of the study, the whole system is divided into the internal power grid, the boundary node, and the external power grid shown in Figure 1.The internal power grid is the part of the power grid that needs to be studied and analyzed.It is not necessary to model other regions in detail, and this part is referred to as the external power grid.The external power grid is connected to each boundary node in the form of an equivalent branch of voltage source series impedance, and the electrical coupling relationship between the external power grids is reflected by establishing the coupling branch between the boundary nodes [14].The equivalent network is shown in Figure 2. In Figure 2, n is the number of boundary nodes, Ui∠θi denotes the voltage and phase angle of boundary node i, Pi and Qi are the active and reactive power input from the external grid to the internal grid through boundary node i, respectively, and i = 1, 2, 3, ..., n; Zij = Rij + jXij is the impedance parameter of the coupling branch between the boundary node i and the boundary node j, where Rij and Xij represent the resistance and reactance, respectively; Zi = Ri + jXi is the impedance parameter of the extended branch between the boundary node i and its corresponding equivalent generator; and Ei and αi are the voltage amplitude and phase angle of the equivalent generator, respectively.The rest of this paper is organized as follows.In Section 2, we present a brief description of the structure and parameters of the equivalent power grid model.In the Section 3, the optimization model of power grid parameters after equivalence and the process of solving equivalent parameters based on genetic algorithm are introduced.Later, results and discussion are provided in Section 4. Finally, Section 5 concludes the paper.

Static Equivalent Network Model of Power Grids
According to the needs of the study, the whole system is divided into the internal power grid, the boundary node, and the external power grid shown in Figure 1.The internal power grid is the part of the power grid that needs to be studied and analyzed.It is not necessary to model other regions in detail, and this part is referred to as the external power grid.The external power grid is connected to each boundary node in the form of an equivalent branch of voltage source series impedance, and the electrical coupling relationship between the external power grids is reflected by establishing the coupling branch between the boundary nodes [14].The equivalent network is shown in Figure 2. In Figure 2, n is the number of boundary nodes, Ui∠θi denotes the voltage and phase angle of boundary node i, Pi and Qi are the active and reactive power input from the external grid to the internal grid through boundary node i, respectively, and i = 1, 2, 3, ..., n; Zij = Rij + jXij is the impedance parameter of the coupling branch between the boundary node i and the boundary node j, where Rij and Xij represent the resistance and reactance, respectively; Zi = Ri + jXi is the impedance parameter of the extended branch between the boundary node i and its corresponding equivalent generator; and Ei and αi are the voltage amplitude and phase angle of the equivalent generator, respectively.In Figure 2, n is the number of boundary nodes, U i ∠θ i denotes the voltage and phase angle of boundary node i, P i and Q i are the active and reactive power input from the external grid to the internal grid through boundary node i, respectively, and i = 1, 2, 3, . .., n; Z ij = R ij + jX ij is the impedance parameter of the coupling branch between the boundary node i and the boundary node j, where R ij and X ij represent the resistance and reactance, respectively; Z i = R i + jX i is the impedance parameter of the extended branch between the boundary node i and its corresponding equivalent generator; and E i and α i are the voltage amplitude and phase angle of the equivalent generator, respectively.

Optimization Model for Network Parameter Solving
For the equivalent power grid model shown in Figure 2, the power input from the external power grid to the internal power grid through the boundary node must be consistent before and after the equivalence: where ∆P and ∆Q are the active and reactive power difference between the external grid input to the internal grid before and after the equivalence, respectively; P ce,i and Q ce,i are the active power and reactive power input from the external power grid to the internal power grid before the equivalence; P i and Q i are the active power and reactive power input from the external power grid to the internal power grid after the equivalence, and their calculation equations are shown in Equation ( 2) and Equation (3), respectively [15].
where θ ir is the difference between the voltage phase angle of the boundary node i and the voltage phase angle of the equivalent generator; θ ij is the difference between the voltage of boundary node i and the voltage phase angle of boundary node j; and Y ij = G ij + jB ij is the admittance parameter of the coupling branch between the boundary node i and the boundary node j, where G ij and B ij denote the resistance and reactance, respectively, and Y i = G i + jB i is the admittance parameter of the extended branch between the boundary node i and its corresponding equivalent generator.
In Equations ( 2) and ( 3), the variables to be estimated are E i , θ i , R i , X i , R ij , and X ij .Because most of the impedance parameters are directly input when the grid is constructed based on the simulation software, the network impedance parameters are selected for the power grid equivalence.This can minimize the equivalent error caused by the difference in the number of decimal points when the network equivalent parameters are calculated, and the simulation software inputs the network equivalent parameters.
Based on the measured data of the internal network, the absolute value and the minimum of the difference between the active power and the reactive power of the external power grid input to the internal power grid through the boundary node before and after the equivalence are constructed as the optimization objective, and the voltage of the equivalent generator and the impedance of the equivalent branch are within a certain range as the constraints.The specific expression is as follows: where E min and E max represent the lower and upper limits of the equivalent generator voltage, respectively; R min and R max represent the lower and upper limits of the equivalent branch resistance, respectively; and X min and X max represent the lower and upper limits of the equivalent branch reactance, respectively.

Network Parameter Solving Method Based on Genetic Algorithm
When the number of boundary nodes is n, the network equivalent parameters to be solved for variables E i , θ i , R i , and X i each have n, respectively, and the network equivalent parameters to be solved for variables R ij and X ij each have n(n − 1), respectively, so the total number of variables to be solved is n 2 + 3n.The traditional network parameter solving method based on internal network information essentially functions by obtaining the boundary measurement information of multiple time periods.There are 2n measurement equations available for each time, and there are 2mn measurement equations for m time periods.To make Equation (4) solvable, the number of measurement equations must be greater than the number of variables to be estimated: With an increase in the number of boundary nodes, the number of measurement equations required to solve the equivalent network parameters also increases linearly.However, when the boundary node measurements of several time periods in proximity do not differ much, the iterative convergence is poor or non-convergent [16], and the value of the boundary nodes obtained by the internal power grid branch breaking or by a disturbance will affect the safe and stable operation of the power grid.In order to overcome this problem, this paper uses the optimization algorithm to solve the equivalent parameters of the equivalent network model.This algorithm only needs the measurement data of the internal power grid in a time period.A flowchart for solving model parameter essence based on the optimization algorithm is shown in Figure 3.
parameters to be solved for variables Rij and Xij each have n(n − 1), respectively, so the total number of variables to be solved is n 2 + 3n.The traditional network parameter solving method based on internal network information essentially functions by obtaining the boundary measurement information of multiple time periods.There are 2n measurement equations available for each time, and there are 2mn measurement equations for m time periods.To make Equation (4) solvable, the number of measurement equations must be greater than the number of variables to be estimated: With an increase in the number of boundary nodes, the number of measurement equations required to solve the equivalent network parameters also increases linearly.However, when the boundary node measurements of several time periods in proximity do not differ much, the iterative convergence is poor or non-convergent [16], and the value of the boundary nodes obtained by the internal power grid branch breaking or by a disturbance will affect the safe and stable operation of the power grid.In order to overcome this problem, this paper uses the optimization algorithm to solve the equivalent parameters of the equivalent network model.This algorithm only needs the measurement data of the internal power grid in a time period.A flowchart for solving model parameter essence based on the optimization algorithm is shown in Figure 3.

Actual Value of Power Flow
Optimization Algorithm

Adjust Model Parameters
Objective Function J The genetic algorithm is an optimization algorithm inspired by evolutionary theory and genetics mechanisms [17,18].It has the advantages of robustness, simple structure, good global convergence, etc.In this paper, the genetic algorithm is used to solve the parameters of the equivalent network model.The main idea is to find the parameter group that makes the objective function obtain the optimal solution through multiple iterations.The calculation process is shown in Figure 4.The main calculation steps are as follows: 1. Input the measured data.Because it is assumed that the power flow data of the internal power grid is unchanged before and after the equivalence, the measured data include the voltage and phase angle of the boundary node and the power flow information of the internal power grid; 2. Set the population size, the maximum number of iterations, and the range of parameters to be optimized and initialize the population.Unsuitable selection of an initial value may lead to problems such as convergence difficulties or non-convergence of subsequent current calculations based on this equivalent network.Since the operation mode of the power grid must be between the maximum and minimum operation modes, the maximum and minimum values of the generator voltage and branch impedance values of the internal power grid under these two operation modes are referred to, and then the upper and lower limits of the voltage of the equivalent generator and the impedance value of the equivalent branch are set; 3. Calculate the fitness value of the population.According to the objective function Equation ( 4), the fitness function is defined as follows: The genetic algorithm is an optimization algorithm inspired by evolutionary theory and genetics mechanisms [17,18].It has the advantages of robustness, simple structure, good global convergence, etc.In this paper, the genetic algorithm is used to solve the parameters of the equivalent network model.The main idea is to find the parameter group that makes the objective function obtain the optimal solution through multiple iterations.The calculation process is shown in Figure 4.The main calculation steps are as follows: 1.
Input the measured data.Because it is assumed that the power flow data of the internal power grid is unchanged before and after the equivalence, the measured data include the voltage and phase angle of the boundary node and the power flow information of the internal power grid; 2.
Set the population size, the maximum number of iterations, and the range of parameters to be optimized and initialize the population.Unsuitable selection of an initial value may lead to problems such as convergence difficulties or non-convergence of subsequent current calculations based on this equivalent network.Since the operation mode of the power grid must be between the maximum and minimum operation modes, the maximum and minimum values of the generator voltage and branch impedance values of the internal power grid under these two operation modes are referred to, and then the upper and lower limits of the voltage of the equivalent generator and the impedance value of the equivalent branch are set; 3.
Calculate the fitness value of the population.According to the objective function Equation (4), the fitness function is defined as follows: where f is the fitness value (the larger the f is, the higher the equivalent accuracy is), and n is the number of boundary nodes.

4.
Roulette selection, crossover, and mutation are performed on the population to obtain the next generation of the population; 5.
According to the number of boundary nodes and the structure of the grid model after the equivalence, considering the estimation accuracy of the equivalence parameter and the actual engineering experience to determine the iteration termination value, detect whether the fitness value of the current iteration number is smaller than the set value.If it meets, the iteration will be ended, and the optimal solution of the iterative process should be outputted.If it does not meet, then it will return to continue the iteration.
where f is the fitness value (the larger the f is, the higher the equivalent accuracy is), and n is the number of boundary nodes.
4. Roulette selection, crossover, and mutation are performed on the population to obtain the next generation of the population; 5.According to the number of boundary nodes and the structure of the grid model after the equivalence, considering the estimation accuracy of the equivalence parameter and the actual engineering experience to determine the iteration termination value, detect whether the fitness value of the current iteration number is smaller than the set value.If it meets, the iteration will be ended, and the optimal solution of the iterative process should be outputted.If it does not meet, then it will return to continue the iteration.

Input internal network measurement data
Set the population size, the maximum number of iterations, the objective function, and initialize the population

Example Analysis
This section addresses the CERPI36v7 node system as an arithmetic example to verify the effectiveness of the model and method proposed in this paper.The CERPI36v7 node system is a common power system simulation model.It consists of 36 nodes, including 8 generator nodes and 8 load nodes, and a total of 37 branches, representing a typical power system.In the simulation, firstly, the equivalent parameters of the external network are calculated based on the genetic algorithm, and then the accuracy of this paper's method is verified by comparing it with the accurate power flow results of the power grid before the equivalence.Finally, by comparing the opening and closing of the grid branches inside the grid before and after the equalization with the results of the current calculation for changing the load of the grid inside the grid, we verify that the external grid equalization model obtained by this method can be adjusted accordingly under the influence of the random measurement error.At present, most of the practical calculations are based on the direct hanging equivalent machine method for power grid equivalence [19][20][21].In this paper, the direct hook-up equalizer method is compared with the method proposed

Example Analysis
This section addresses the CERPI36v7 node system as an arithmetic example to verify the effectiveness of the model and method proposed in this paper.The CERPI36v7 node system is a common power system simulation model.It consists of 36 nodes, including 8 generator nodes and 8 load nodes, and a total of 37 branches, representing a typical power system.In the simulation, firstly, the equivalent parameters of the external network are calculated based on the genetic algorithm, and then the accuracy of this paper's method is verified by comparing it with the accurate power flow results of the power grid before the equivalence.Finally, by comparing the opening and closing of the grid branches inside the grid before and after the equalization with the results of the current calculation for changing the load of the grid inside the grid, we verify that the external grid equalization model obtained by this method can be adjusted accordingly under the influence of the random measurement error.At present, most of the practical calculations are based on the direct hanging equivalent machine method for power grid equivalence [19][20][21].In this paper, the direct hook-up equalizer method is compared with the method proposed herein in order to comparatively verify the equalization accuracy of the method proposed in this paper.
herein in order to comparatively verify the equalization accuracy of the method proposed in this paper.

Simulation Objects and Conditions
The simulation example system is the CERPI36v7 node system.As shown in Figure 5, four boundary nodes are set up in the CERPI36v7 node system, namely, boundary nodes 9, 14, 16, and 19.Nodes 6,7,8,11,12,13,25,26,27,28,29,30,31,33,34,51, and 52 belong to external power grid nodes, and other nodes belong to internal power grid nodes.The purpose of the static equivalence of a power system is mostly to analyze the static security of the power system, which is used to judge whether the system will overload or reach its voltage limit after the expected accident.The analysis results are related to the absolute error of power flow calculation.Therefore, this paper uses the average relative error of power flow and the maximum error of power flow to measure the equivalent accuracy of the equivalent method.
The average relative error of ξave is as follows: ,eq ,eq ,eq ,eq ave ----= ave + ave + ave + ave where Pl, Pl,eq, Ql, and Ql,eq represent the active power and reactive power values transmitted by the internal power grid lines before and after the equivalence, respectively; and Ul, Ul,eq, θl, and θl,eq represent the voltage amplitude and phase angle of internal power grid nodes before and after equivalence, respectively.Ave represents the mean value of the expression.
The maximum relative error of ξmax is as follows: ,eq ,eq ,eq ,eq max - where max denotes the maximum value of the expression.The purpose of the static equivalence of a power system is mostly to analyze the static security of the power system, which is used to judge whether the system will overload or reach its voltage limit after the expected accident.The analysis results are related to the absolute error of power flow calculation.Therefore, this paper uses the average relative error of power flow and the maximum error of power flow to measure the equivalent accuracy of the equivalent method.
The average relative error of ξ ave is as follows: where P l , P l,eq , Q l , and Q l,eq represent the active power and reactive power values transmitted by the internal power grid lines before and after the equivalence, respectively; and U l , U l,eq , θ l , and θ l,eq represent the voltage amplitude and phase angle of internal power grid nodes before and after equivalence, respectively.Ave represents the mean value of the expression.The maximum relative error of ξ max is as follows: where max denotes the maximum value of the expression.

Results of Network Parameter Solving Based on Genetic Algorithm
Firstly, the network parameters of the CERPI36v7 node system are solved based on the genetic algorithm.The results of the external network parameter solving iterations for the CERPI36v7 node system are shown in Figure 6.The direct hanging equivalent machine method uses the equivalent machine to characterize the power input from the external power grid to the internal power grid, and the power of the equivalent machine is the sum of the power input from the external power grid to the internal power grid through the boundary nodes, as shown in Table 1.

Results of Network Parameter Solving Based on Genetic Algorithm
Firstly, the network parameters of the CERPI36v7 node system are solved based on the genetic algorithm.The results of the external network parameter solving iterations for the CERPI36v7 node system are shown in Figure 6.The direct hanging equivalent machine method uses the equivalent machine to characterize the power input from the external power grid to the internal power grid, and the power of the equivalent machine is the sum of the power input from the external power grid to the internal power grid through the boundary nodes, as shown in Table 1.The obtained parameters of the external grid network are substituted into the equivalent grid model, and the power flow is calculated.The average relative error of the power flow and the maximum relative error of the power flow are shown in the following Table 2.It can be seen from Table 2 that the average relative error and the maximum relative error of the power flow between the method proposed in this paper and the method of directly hanging the equivalent machine are kept within 1%, which meets the  The obtained parameters of the external grid network are substituted into the equivalent grid model, and the power flow is calculated.The average relative error of the power flow and the maximum relative error of the power flow are shown in the following Table 2.It can be seen from Table 2 that the average relative error and the maximum relative error of the power flow between the method proposed in this paper and the method of directly hanging the equivalent machine are kept within 1%, which meets the requirements of the equivalent accuracy.The equating process of this paper's method is to equate a state of the grid into the form of an equipotential machine plus an equipotential branch, while the process of the direct hanging equivalent machine method is to directly replace the power input from the external power grid to the boundary node with an equal power equivalent machine and substitute it into Equation (6).Therefore, the fitness value calculated by the direct hanging equivalent method is infinity.When the internal power grid state does not change, the power flow error of the direct hanging equivalent method is relatively low.

Load Power Growth
Table 3 provides the average relative error of the tidal current, the maximum relative error of the tidal current of the CERPI36v7 node system for this paper's method, and the direct hooked equalizer method for the internal grid node 21 load power growth.It can be seen from the table that the average relative error of the proposed model and method is kept within 1.1%, and the maximum relative error of the power flow is kept within 4%, which meets the requirements of the equivalent accuracy.However, the average relative error of the maximum power flow of the direct hanging equivalent method exceeds 5%, the maximum relative error of the power flow exceeds 16%, and the power flow error of the internal network may be greatly increased due to an error of the equivalent network.Although the direct hanging equivalent method is equivalent based on the power flow results in the initial state, the equivalent results can accurately reflect the external power grid without any change in the power grid.However, when the internal state of the power grid changes, the average relative error and the maximum relative error of the power flow of the equivalent results of the direct hanging equivalent method begin to increase significantly, which cannot accurately reflect the external power grid when the power grid changes.Compared with the direct hanging equivalent machine method, the grid equivalent model established by the proposed method can still maintain high accuracy when the grid state changes and is more suitable for the actual situation of power flow state changes.

Line Breaking
Table 4 shows the equivalent model error of the CERPI36v7 node system breaking branches 9-23, branches 22-23, branches 21-22, and branches 16-18 under four kinds of internal power grid state changes.It can be seen from the table that the method proposed in this paper achieves a lower equivalent error than the direct hanging equivalent method, whether it is the relative error of the power flow or the maximum error of the power flow.Although the equivalent error is large when the internal power grid topology changes, (e.g., the average relative error of the power flow reaches 96.62% when line 22-23 is disconnected), and although the maximum relative error of the power flow exceeds 100%, the proposed method in this paper can be carried out in real time, and the interval of each equivalence is short.When the topology of the internal power grid changes, the next equivalence will begin soon, which can meet the requirements of power flow calculation accuracy.In addition, whether it is the method proposed in this paper or the direct hanging equivalent machine method, the equivalent error caused by line interruption is large.This is because when the line is disconnected, the power flow near the boundary node will change greatly, and the greater the power transmitted on the line before the interruption, the greater the fluctuation in the power flow redistribution.
From the above simulation results, it can be concluded that the method proposed in this paper is superior to the direct hanging equivalent machine method used in the actual power grid.When using the direct hanging equivalent machine method, one inputs the constant equivalent power into the internal power grid through the boundary node.When the state of the internal power grid changes, it is bound to produce a relatively large equivalent error.The model used in this paper reflects the electrical coupling relationship between the external power grids by establishing the coupling branches between the boundary nodes.The external power grid has a certain supporting effect on the power and voltage of the internal power grid.When the state of the internal power grid changes, the equivalent machine of the external power grid can make corresponding adjustments to reduce the equivalent error.

Conclusions
In this paper, a static equivalent method of power grids based on the genetic algorithm is proposed.Based on the measurement information of the internal power grid in a certain period of time, the method constructs the grid equivalence network and the optimization model for solving the network parameters.The difference in the power input from the external power grid to the boundary node before and after the equivalence is utilized as the objective function.The equivalent parameters of the external network are solved by the genetic algorithm, which avoids the influence of an unplanned branch interruption of the internal network on the power grid and effectively improves the adaptive range of solving the equivalent parameters of the external grid.Simulation of the CERPI36v7 node system verified the accuracy and effectiveness of the method proposed in this paper.The method proposed in this paper is simple to calculate and its potential application is wide; further, it is expected to be applied in practical engineering.

Figure 3 .
Figure 3. Flowchart for solving model parameter essence based on optimization algorithm.

Figure 3 .
Figure 3. Flowchart for solving model parameter essence based on optimization algorithm.

Figure 4 .
Figure 4.The flowchart for solving model parameters based on genetic algorithm.

Figure 4 .
Figure 4.The flowchart for solving model parameters based on genetic algorithm.

Figure 6 .
Figure 6.Iterative results of equivalent parameters of external network of CERPI36v7 node system.

Figure 6 .
Figure 6.Iterative results of equivalent parameters of external network of CERPI36v7 node system.

Table 1 .
The external network equivalent parameters calculated by the direct hanging equivalent machine method.

Table 2 .
Equivalent error of power flow results of different equivalent methods.

Table 1 .
The external network equivalent parameters calculated by the direct hanging equivalent machine method.

Table 2 .
Equivalent error of power flow results of different equivalent methods.

Table 3 .
The equivalent error of the internal grid node 21 when the load power increases.

Table 4 .
Equivalent error of internal power grid line breaking.