Phase Synchronization Stability of Non-Homogeneous Low-Voltage Distribution Networks with Large-Scale Distributed Generations

: The ideal distributed network composed of distributed generations (DGs) has unweighted and undirected interactions which omit the impact of the power grid structure and actual demand. Apparently, the coupling relationship between DGs, which is determined by line impedance, node voltage, and droop coefﬁcient, is generally non-homogeneous. Motivated by this, this paper investigates the phase synchronization of an islanded network with large-scale DGs in a non-homogeneous condition. Furthermore, we explicitly deduce the critical coupling strength formula for different weighting cases via the synchronization condition. On this basis, three cases of Gaussian distribution, power-law distribution, and frequency-weighted distribution are analyzed. A synthetical analysis is also presented, which helps to identify the order parameter. Finally, this paper employs the numerical simulation methods to test the effectiveness of the critical coupling strength formula and the superiority over the power-law


Introduction
Renewable energy sources have a broad prospect in power systems these days-more and more distributed generations (DGs) are integrated into the modern electricity grid. The distributed generation has the advantages of being environmentally-friendly and sustainable. In addition, two weaknesses include low output power per unit and energy output in remote areas that can not be ignored. For one thing, the low output power per unit of solar power and wind power results in a low capacity for one DG unit. This is why the distribution network requires a large number of installed DGs to provide effective output power. For another thing, the remote distribution leads to power systems usually working in the islanded operation. Different from the normal generators, there is no rotor in the DGs such as batteries and photovoltaic units. That is to say, the frequency stability is easy to be broken in the islanded microgrid composed of DG units without the rotating mass [1][2][3]. Hence, we will address the issue of phase and frequency stability in the islanded system.
Previous studies have based their criteria for selection on the local control of a few DGs in the ideal model [4][5][6]. Reference [7] shows the change of various parameters of generators in a large scale photovoltaic system and presents two different operations to research the voltage stability in an IEEE 30-bus model with a Large-Scale photovoltaic module. Indeed, the large-scale photovoltaic module can be considered as one high-capacity DG in the case of IEEE 30-bus, and no method of improving stability was given in [7]. Conventionally, when the number of composed DG increases, the active low-voltage islanded network will face the issue of high-order, multi-dimensional, and strong coupling. In brief, with a sharp increase in the number of DGs, the system will encounter the problem named "Curse of Dimensionality". Therefore, it is difficult for previous solutions that consider a few DGs to address relevant problems in current research effectively.
We note that recently there has been some advanced complex network theory that provides the stability analysis of large-scale network systems. It has wide applications in cyber science [8], medical science [9], social science [10], and so on. Nevertheless, there are a few related studies on the frequency stability of an active distribution network composed of DG units.
In addition, previous research mainly assumes that the networks are unweighted and undirected; this means that the coupling relationship between nodes has no direction, and the strength of all couplings is consistent. However, this assumption mismatches the characteristics of the real network of many cases. For instance, in an active distribution network, the coupling relationship between DGs which takes droop control is determined by the node voltage, line impedance, and droop coefficient, where the line impedance between DGs is generally determined by the distance between DGs. Apparently, the assumption that all DGs in the active distribution network is equidistant is undoubtedly harsh and impractical. On the contrary, the assumption that the coupling relationship between DGs is non-homogeneous is more realistic. Fortunately, complex network theory has made a difference in weighted networks. Reference [11] investigates the system of Kuramoto oscillators whose weight is related with a degree in the asymmetrical situation and comes to one conclusion that the critical coupling strength has a positive correlation with the degree exponent. Reference [12] indicates that the degree sequence contributes to the ability to synchronize in the mean field. Reference [13] shows that the degree distribution of the fixed system has a negative degree correlation between its nodes, which improves its synchronizability. Besides the weighted methods with the degree, some other papers pay attention to the weighted methods with frequency. Reference [14] introduces a weighting method which cares about the relationship between frequency and link to achieve the phase synchronization. Linear and nonlinear weighting procedures are applied in mean and heterogeneous systems to verify the method. Frequencies of symmetrical oscillators are studied in the bipolar model built by [15], while the relationship between frequency and coupling strength is analyzed. Reference [16] comes to the conclusion that asymmetric weights that depend on loads make contributions to the stability of the whole system. Overall, weighted methods that are based on the inherent qualities of the inverter nodes [17,18], i.e., node degrees [19], node loads [20], and network's global properties in the connection links have been discussed widely to enhance the network synchronization in recent studies.
However, the general impact of weighted coupling on network synchronization is not studied comprehensively, especially the weighted method. These weighted procedures above, unfortunately, are designed only concerned with the inherent attributes of the network, i.e., the size of the nodes degree or the edges between vertices, which limits their applications. Thus, these weighted procedures may have a significant difference from real systems. Therefore, it is important to describe such real weighted distribution methods in weighted networks.
Thereby, it is of great importance to study different weighted distribution methods of a distribution network, which has not been discussed yet. This paper employs some weighted distribution methods that are found in many realistic networks such as the link weight and investigates the phase synchronization performance in a large-scale islanded network to pave the way for the solution of the similar weighted distribution issues in the future. The innovations of this paper are as follows.
(1) We derive non-homogeneous models for simulating different weighted distribution methods in a large-scale islanded distribution network. (2) We analytically study the effects of weighted distribution methods on the phase synchronization stability by means of explicit synchronization conditions. We illustrate our insightful results and the utility of our approach through the critical coupling strength formula.
(3) Different from other studies focusing on the weighted distribution methods which only associate with the inherent properties of the network, this paper adopts several realistic distribution styles as weighted distribution methods to research the synchronization stability in the large-scale distribution network.
The rest of this paper is organized as follows. In Section 2, we introduce correlation complex network theory. Section 3 focuses on the modeling of weighted distribution and derivation of the critical coupling strength formula. The numerical results and simulation analyses in Section 4 are used to compare the network stability of different weighting methods. Section 5 concludes this paper.

Correlation Theory
The structure of an active distribution network can be achieved by mirroring algebraic graph theory. Furthermore, the synchronization condition of an active distribution network will be analyzed by order parameters that are introduced thoroughly in the following. Then, the critical coupling strength of network stability can be obtained with the help of synchronization condition.

Graph Theory
The microgrid can be modeled as a graph G, where V= {θ 1 ,θ 2 , · · · θ N } is the set of edges while E ∈ V × V is the set of edges. Generally, the connected graph can be called an undirected graph when (θ i , θ j ) and (θ j , θ i ) represent the same edge. The value of weight a ij defines whether the graph is weighted (a ij = 1) or not (a ij = 1) and hence we have the associated adjacency matrix A = [a ij ]. Moreover, the Laplace matrix L is a symmetric zero-sum matrix of undirected graphs, where l ii = ∑ i =j a ij , l ij = −a ij . The matrix L satisfies: Whether the graph G is a directed graph or an undirected graph, its associated matrix is based on a directed label, which means edge e = {i, j} ∈ E is represented by an ordered pair (i, j). In the undirected graph, we assign an arbitrary direction to it. For the ordered pair (i, j), we define b je = 1 if j is the sink node of edge e, b ie = −1 if i is the source node of edge e, and b ke = 0 otherwise in the incidence matrix B = [b ie ] ∈ R n×m . In addition, the edge (i, j) satisfy (B T X) e = x j − x i .

Synchronization Conditions
Specifically, the phase and frequency of all DG units tend to stabilize at a fixed value when t → ∞. Therefore, we focus on the synchronization process of phase and frequency to research the system stability intuitively.
If the coupling strength K between the DG nodes in a complex network is larger than the critical value K cr , the DG nodes will divide into two groups, the phase of one cluster is fixed at a stable specific frequency while the other oscillators have a drift.
The order parameter is a visual signal to show the change of phase, whose definition is: The size of r shows the percentage of stable nodes. The value is: where r = 1 means that a system enters a stable state absolutely while r = 0 means that a system becomes unstable. Criterion 1 [21]: For a complex network graph G with an incidence matrix B, each pair of interconnected oscillators {i, j} ∈ ε could enter frequency synchronization and phase convergence θ i − θ j ≤ γ < π/2 when satisfying where L † is the inverse of Laplace matrix L, ω= (ω 1 , · · · ω n ) T is the vector of the natural frequency of each oscillator in the network and · ε,∞ evaluates the worst-case dissimilarity of this weighted projection. This paper will use this synchronization condition to obtain the critical coupling strength of network synchronization under weighted conditions.

Critical Coupling Strength
In the previous subsection, the synchronization conditions recently mentioned in [22] have been introduced. In different occasions, these conditions have been proven to be applicable to oscillator models. This paper takes the Kuramoto oscillator to express the nonlinear ordinary differential equations: The model can be expressed as a matrix: W := diag(w) is the diagonal matrix of weights. In a complex network, which determines the stability of the system is not the absolute value of the phase but the phase difference. The network entry synchronizes only when two phases of coupled nodes ϕ 1 and ϕ 2 are locked at a certain ratio, which means |nϕ 1 − mϕ 2 | <constant. Therefore, the phase difference vector is ϕ(t) := B T θ, and the Equation (6) can be expressed by the phase difference: In the above equation, the synchronous solution can be expressed asφ=0 orφ=ϕ * for a fixed point ϕ * ∈ R n .
In an analysis of [23], Fazlyab found that all fixed points of Equation (7) can be expressed as: where F ∈ R m×m−n+1 is a matrix whose columns span the null space of B, i.e., BF = 0, y ∈ R m−n+1 , and r is an arbitrary real number. When F = 0 and the fixed point ϕ * is unique, the particular solution W r−1 B T (BW r B T ) † ω is independent of r. In the special case of r = 0, we have that This section has proposed the theory about synchronization. The next section will analyze the situation of different weighting methods in star topology through algebraic graph theory, and obtain the critical coupling strength by using the synchronization condition.

Weighted Distribution
The weighted distribution is not only helpful to describe the features of real networks but are also crucial for controlling network synchronization. Three different weighting methods will be discussed to research the phase synchronization of a large-scale islanded network in this section. This issue takes an ideal star topology as the research object. Furthermore, a critical coupling strength formula for the weighted star network would be proposed first in this paper.

Star Topology
The purpose of this investigation is to explore the synchronization stability of large-scale grid-tied DG in islanded operation. Considering that most of the distribution networks are star topology [24,25], this paper chooses star topology as the object, as shown in Figure 1. The star network is shaped by one primary node and many branch nodes connected with it. It is characterized by high reliability, easy operation, and easy maintenance when the branch lines fail. However, this structure will consume a lot of metal and switchgear. Therefore, it is suitable for important construction sites with large capacity equipment or special places with concentrated load and wet, corrosive environments [21,26,27]. Due to the practical application of the power grid, the traditional research object of multi-machine (three to five machines) can no longer meet the actual needs. This is why this paper delivers an explicit analysis of the network stability with large-scale grid-tied DG. The number of DG N in a Large-scale model generally refers to N ≥ 200, where this paper takes N as 201. There are no conventional generators with rotors used in this model, and all DGs are photovoltaic units. The photovoltaic community electrical diagram with star topology is depicted in Figure 2. We assume that the value of the line impedance between load nodes and the master node is consistent, and the inverter of each load node has the same impedance, so the impedance value between load node and the master node is equal. Therefore, Z = Z line + Z inv and the admittance value is In Figure 2, each load node is only connected to the master node L1, and the impedance value between nodes Z ij = Z + Z inv are equal. Therefore, the admittance matrix with star topology is

The General Formula for Weighted Distribution
The DG units with droop control can be described as a Kuramoto oscillator [28][29][30]: where m i p iθ =ω i , ω i is the rated frequency of inverter in the ith DG node. In addition, the value of ω i is chosen from the probability density function g(ω). When the standard frequency is 50 Hz, we take g(ω) as the symmetric normal distribution of 50, therefore g(ω)∼ (50, 1).
In the above equation, the weight between DG i and DG j nodes is obtained by the droop coefficient d i , the voltage of node E i , E j , and the admittance Y ij . At this point, the coupling strength For the sake of simplicity, we divide the weight between the DG nodes into two parts. One is the coupling strength k, and the other is the weight index m i that conforms to some kind of weighted distribution. Therefore, the weight of the directed edge {i, j} is determined by the weight index m i while the weight of the edge {j, i} is determined by the weight index m j , so the weight distribution of one network can be decided by the distribution of weight index m i . In addition, we can synchronously adjust the size of all m i values to change the size of the weight. Thus, we have the weight diagonal matrix W = k × diag(m), where m 1 is the weight when a load node is the sink node and the master node is the source node, m 2 · · · m n is the weight when the master node is the sink node and a load node is the source node. Thus, the weight distribution of the network can be simulated by adjusting the value of the weight index: In a star topology whose N is 201, the incidence matrix is

Weighted Distribution Model
Different weighting methods have different performances in different topology, but the topology of an actual power grid is fixed and difficult to adjust. Thus, the weighting methods, which are associated with the degree distribution and node strength, are not considered in this paper. Only the weighting method associated with the distribution of coupling strength is considered in this paper.
In practice, the error occurs in the installation process often obey Gaussian distribution, the weight distribution of complex network is mostly in line with the power-law distribution, and the frequency-weighted distribution can reflect some peculiarities of the actual network. Thus, this paper discusses the weighting methods of Gaussian distribution, power-law distribution, and frequency-weighted network, and the aforementioned methods can be simulated by adjusting the weight index.

Gaussian Distribution
At first, this paper considers the case that the weight distribution method is the same as the rated frequency of an inverter distribution method so that the distribution function g(m) of m i is the same as the distribution function g(ω) of ω i ; both are Gaussian distribution. Unlike the rated frequency of the inverter being a symmetric normal distribution of 50, the weight distribution is a symmetric normal distribution of the y-axis. Its Kuramoto oscillator model is: a ij = 1 when DG i has coupling relationship with DG j , otherwise a ij = 0.

Power-Law Distribution
Power-law distribution has extensive application in complex network theory. Barabasi and Albert pointed out that the degree distribution of complex networks is the power-law distribution in [31]. Then, Chunguang Li and Guanrong Chen verified that the weight distribution is a power-law distribution in a general complex network. Furthermore, according to the research of Barrat et al. in [32], the node strength distribution in the weighted network obeys the power-law distribution as well. In this paper, only the weighted distribution of weights is discussed. Therefore, it is worth mentioning that the assumed distribution function g(m) of m i satisfies the power-law distribution, which means g(m) ∼ m −β , anḋ

Frequency-Weighted Network
The frequency-weighted network reflects the feature of several natural and social networks. In fact, a power grid could be expressed as a Kuramoto oscillator model, and the weight of the coupling strengths are obtained by their natural frequency. In this case, the weighted coupling coefficients are related to DGs' own inverter rated frequency rather than the topology of a network. We added the inverter rated frequency to the dynamic model as follows:

Critical Coupling Strength Formula
Reference [23] pointed out that any fixed point in Equation (7) has local exponential stability when ϕ * ∞ ≤ π/2. For the sake of simplicity, we decompose the matrix A = k × A , where A is the adjacency matrix entirely determined by weight indexes, and k is the coupling strength. The in-degree matrix D is a diagonal matrix, whose values in each row are the row sum of corresponding rows of A matrix. Certainly, the D matrix can be decomposed into D = k × D , where D is the in-degree matrix entirely determined by weight indexes. At the same time, the Laplacian matrix L = D − A can be decomposed into L = k × L , where L is the Laplacian matrix entirely determined by weight indexes. Furthermore, Equation (9) can be decomposed into: There is sin(B T θ * ) ∞ = 1 when the network enters a critical stable state, and then the critical coupling values of star topology in different weighting methods can be transformed into the following formula: The calculation of the critical coupling value in different weighting methods is described in detail in Section 4.

Simulation
In this section, three weighted methods from above are simulated, and the proposed formula of critical coupling (18) is verified.

Case 1
It is inevitable that deviations will generate in the process of installing the grid-tied DG inverters in an active distribution network, and the distribution of deviations often conforms to the Gaussian distribution.
In this case, the assumed actual parameters are reported in Table 1: Let us now turn to the case of the undirected situation, the weight distribution of edges in this situation conforms to the Gaussian distribution X ∼ N(0, 1), and source point and sink point of each edge are arbitrary. When N takes 201, and the number of edges takes 200 in this case; meanwhile, the distribution is depicted in Figure 3. Furthermore, Figure 4 shows the probability density function of this Gaussian distribution, and it is obvious that the distribution law conforms to the (0,1) distribution.    Thus, the ordinary differential equation of oscillator θ i iṡ where X is the weight of the corresponding edge. Turn now to verifying the critical coupling strength formula of the weighted situation proposed in Section 3. In this case, Furthermore, the related parameters in this case are depicted in Table 1, and the weight 100×0.153+10×1.15 = 180. The phase and frequency stability are shown in Figures 5 and 6, and it can be seen that the phases of DG nodes change at the same rate and the frequency of DG nodes fixed in one constant value.  In this paper, order parameters are used to reflect the synchronization of DGs in the network. Figure 7 shows the relationship between the coupling strength K and the order parameter r when N has different values, which reflect the change of network synchronization with the coupling strength. It is obvious that the network stability is the best when the value of N is 51, followed by 101. When N takes 201, the network has the worst stability. It can be seen that the stability of the network decreases with the increase of N, which is not consistent with the situation in which the stability does not change with the value of N when the star topology is an undirected network.

Weighted Directed Network
Thus far, this paper has focused on an ideal situation that the weighted distribution of the undirected edge in the network is Gaussian distribution. Let us now turn to consider the weighted directed case. In a power grid, each DG can become not only the power transmitter but also the power receiver. Therefore, it is of great practical value to consider the different weights of edges when they take different directions between coupling DG nodes. That is, the edge weight is the weight index of the master node when the master node is the source point, while the edge weight is the weight index of the load node when the load node is the source point: Since the inverter rated frequency of the DG unit also conforms to the Gaussian distribution, this paper assumes that the weight distribution of edges is the same type as the distribution of inverter rated frequency. Namely, the ordinary differential equation with the oscillator θ i iṡ This paper takes the weight index of the load node as the Gaussian distribution and the weight index of the master node as the mean value of the weight index of each load node to realize the power balance. Thus, there is Frequency deviation will lead to power imbalance and even instability of the whole system. The phase and frequency stabilization processes in the directed situation are shown in Figures 8 and 9 when k = 180. Obviously, the coupling strength k in this case is much greater than the critical coupling strength required for network stability. Let us now turn to the critical coupling strength in the case of Gaussian distribution. When N = 201, the value of L matrix can be substituted for formula (18) to obtain the critical coupling strength K C = 6.7; when N is 51 and 101, the critical coupling strength K C = 3.7 and 4.46 can be obtained from the calculation formula, respectively. Figure 10 shows that, when N takes 201, the network enters synchronization rapidly before K < 6.7, while the network is nearly stable after K > 6.7. This result is basically consistent with the result obtained by formula . Similarly, when the value of N is 51 and 101, the network becomes stable after K C = 3.7 and 4.46, respectively. Thus far, two conclusions have been shown. The first is that the validity of the critical coupling strength formula (18) in the directed weighted network is verified, and the second is that the network stability decreases with the increasing number of nodes in the island star topology directed weighted network.

Case 2
In particular, previous studies show that most of the nodes have a degree that is a less than average value while a few nodes have large connections in the real network. The study of Chunguang Li and Guanrong Chen [17] confirmed that the distribution of connection weights would follow the power-law distribution in general actual complex networks. Thus, it is important to research the stability of a weighting case that conforms to the power-law distribution.
In this case, the assumed actual parameters are reported in Table 2: In this case, we assume that the weight distribution of 200 edges conforms to the power-law distribution, as shown in Figure 11. Then, Figure 12 shows the probability density function diagram. It can be seen that, with the increase of abscissa, the probability of weight distribution is lower. Furthermore, the probability of randomly selecting a node with greater weight is negatively correlated with the number of the same weight. The power-law index in Figure 11 is 2. The distribution of nodes depicted as a line whose slope is negative in the double logarithmic diagram, as demonstrated in Figure 12. This distribution has an asymptote of the black dotted line, where the exponent determines the rate of decay.  The weighted distribution in this case follows g(X) ∼ X −β , which means that the ordinary differential equation of oscillator θ i iṡ where X is the weight of the corresponding edge.
Like the weighted undirected network in case 1, the weighted distribution in this case is reflected by the value of the edge weight, so the Laplace matrix is The same as the coupling strength value of case 1, the phase and frequency stabilization process diagram are shown in Figures 13 and 14. In addition, the influence of different power exponents on network stability is shown in Figure 15. Figure 15 shows the network synchronization process when power-law index β grades -5, -3, -2, 2, 3, 5. It is obvious that the critical coupling strength of the system when the power-law index is positive does not differ greatly, while the positive power-law index values have significantly better stability than the negative values. This paper comes to a conclusion that the positive and negative value of power-law index rather than the size of the Power-law index has a key influence on the stability of the network when ignoring the deviations generated in the calculation process.   The simulation analysis of the different number of nodes N is shown in Figure 16. When the index is 2, the value of N is negatively correlated with network stability. The conclusion also reflects the analysis of the relationship between the node number and the network stability in the previous case.

Weighted Directed Network
Turn now to the directed weighted case. Similar to case 1, the method of adjusting the weight index is used to complete the weight distribution conforms to the power-law distribution g(m) ∼ m −β , and the ordinary differential equation is given in the formula (25): Similar to the weighted directed network in case 1, the weight index of the master node in this case is the mean value of the weight index of each load node considering the problem of power balance. Therefore, the Laplace matrix in the network with 201 nodes is: The phase and frequency stabilization process in the directed situation are shown in Figures 17  and 18 when the value of power-law index beta is 2, and coupling strength k is 180. When the value of the power-law index respectively takes -5, -3, -2, 2, 3, 5, the critical coupling strength K C of network is respectively 1.4976, 1.4353, 1.9639, 1.4976, 1.4353, and 0.7771 by formula (18). These critical coupling strengths are nearly consistent with the critical value depicted in Figure 19. In addition, the weighted directed network has similar performance to the weighted undirected network, which means that the positive power-law index value has better network stability performance than the negative one, while the size of the power-law index value has almost no effect. Apparently, in the case of power-law distribution, the stability of the network is obviously better than that in the case of Gaussian distribution. When the Power-law index is 2, the network can be stable only on the condition that the coupling strength is greater than the critical value K C = 0.5475. Furthermore, the relationship between the number of DG nodes and network stability is shown in Figure 20. The negative correlation between DG number and critical coupling strength is still evident.

Case 3
The following case will discuss the Frequency-weighted situation. The natural frequency is inverter rated frequency in the Kuramoto oscillator model. Thus, case 3 chooses the absolute ratio between the inverter rated frequency of one DG node and its in-degree as the weight; there is the ordinary differential equation: The weighted distribution in this case is different from the previous two cases, whose weight index is not a random value which obeys a certain distribution. The Laplace matrix in this case is In case 3, the phase and frequency stabilization process diagram are shown in Figures 21 and 22, which are similar to the above two cases. According to formula (18), the critical coupling strength of this case is K C = 12.89 when N = 201, and the network synchronization progress is shown in Figure 23. It can be seen that the weighted method in this case is not as stable as the Gaussian distribution and the power-law distribution.

Stability Analysis
Through comparing the critical coupling strength of three cases as illustrated in Figure 24, we find that the power-law distribution has the best stability in different weighted distribution methods. Moreover, the second largest eigenvalue λ 2 or the eigenvalue proportion λ N /λ 2 in adjacency matrix A could describe the synchronization ability of the system. The network coherence is proportional to λ 2 and inversely proportional to λ N /λ 2 . Table 3 shows that the λ 2 of power-law distribution has the smallest value, which means it has the best synchronization ability, followed by the Gaussian distribution and frequency-weighted distribution. Thus, the results of Table 3 and Figure 24 confirm each other. In addition, the contrastive analysis of weighted directed network and weighted undirected network show that the weighted directed network which conforms to the Gaussian distribution or power-law distribution is easier to get into synchronization than the weighted undirected network. Furthermore, by studying the network with different numbers of DG, this paper finds that the stability of the weighted network decreases with an increase of N, which is different from the average field. For the weighting method of power-law distribution, the size of the power-law index β is not significant to the stability of the network, but the sign of the power-law index β is not negligible to the stability of the network.

Further Analysis of Weighted Network in the Actually Distributed Network
Actually, the topology structure of an actual power grid network often isn't a standard star topology. The standard star topology is an ideal assumption. In order to verify the stability analysis from the above subsection, we take the dynamic IEEE test systems [33] and the Iceland 189-node distributed grid [34] as a research object. The synchronizing processes of Gaussian distribution, power-law distribution, and frequency-weighted network are shown in Figures 25 and 26.

Dynamic IEEE Test Systems
Three different weighted distribution methods have been applied to the dynamic IEEE models, including 9-bus, 14-bus, 57-bus, and 118-bus modified test systems. The diagram of coupling strength and order parameter is shown in Figure 25, where the labels G, P, and F represent Gaussian distribution, Power-law distribution and frequency-weighted distribution, respectively, and the network with the Power-law distribution method always enter the stable state at first. The value of the critical coupling strength of three weighted distribution methods in the dynamic IEEE test systems is shown in Table 4, where each value is calculated by Equation (18). It shows that the network with the power-law distribution method always has the least critical coupling strengths, which corresponds to the stability performance in Figure 25.   Figure 26 shows the results of the stability performance of three distribution methods in an Iceland 189-nodes distribution grid. It can be seen that the power-law distribution undoubtedly has the best stability performance. Second is the frequency-weighted network. Moreover, the Gaussian distribution needs a long process to enter synchronization. These results further support the conclusions of the above subsection.

Conclusions
Above all, this paper analyzes the stability of large-scale grid-tied DG in a non-homogeneous active power network by using complex networks theory.
Then, three different weighting distribution methods, including Gaussian distribution, power-law distribution, and frequency-weighted distribution, are analyzed. In addition, the critical coupling strength formula for the weighted situation k c = B T (L) † ω ∞ has been proposed and verified.
The critical coupling strength obtained from this formula conforms to the reality when ignoring tiny deviation. Among the three different weighting methods, the power-law distribution has the best stability. Moreover, the difference of power-law index value has little effect on the network stability, but the sign of power-law index has quite an effect on the system stability. The second is the Gaussian distribution, which is obviously less stable than the power-law distribution. Finally, the frequency-weighted network has the worst synchronization performance, which is significantly inferior to the other two. In addition, further analysis of the IEEE dynamic test systems and the Iceland 189-node distributed grid has proved the better stability of power-law distribution in this paper.