Abstract
In order to realize the optimal planning of grid–resource–storage for distribution networks (DNs) with high penetrated distributed photovoltaics (PVs), a cluster partition-based two-layer expansion planning for DNs is proposed. First, a comprehensive cluster partition index-based cluster partition method is proposed, which involves the indexes such as electrical distance, power balance of the cluster, and cluster size. Second, a cluster partition-based two-layer expansion planning model is proposed. In the upper layer, a line planning model for clusters is established to carry out the planning of cluster connection lines. In the lower layer, a robust source-storage planning model is established with the uncertainty of PVs and loads, and then, the optimal location and capacity of PVs and energy storages (ESs) can be obtained. In addition, the uncertainty regulation parameter is utilized to control the range of uncertainty sets, which can reduce the conservatism of the optimization. Finally, the proposed method is carried out in a real DN in China, which can effectively improve the economy of DN planning.
1 Introduction
With the rapid growth of energy demand, photovoltaics (PVs) are developing rapidly in China. The large amount of distributed PVs has significantly changed the power flow of the distribution network (DN) (), which poses new challenges to DN planning and operation (). How to carry out optimal DN planning is the key to realize the economic operation of the DNs with large-scale distributed PVs.
Currently, the main idea for planning the DNs with distributed PVs is to build a centralized planning model () by strengthening or extending the lines of the DNs (), which takes into account the investment of PVs () and the operating costs (). The centralized planning model is suitable for the DN planning when the proportion of distributed PVs is low. However, when large-scale distributed PV is connected to the DNs, the dimensionality of the variables in the centralized planning model increases significantly (), and the planning model becomes too complex to be solved (). To solve the challenges of centralized optimization, a cluster partition-based planning method provides a new way for DN planning. The cluster partition-based planning method can not only decompose the centralized optimization problem into simple sub-problems of cluster optimization but also maximize the degree of power matching between PVs and the load within the cluster during the planning process, which can greatly increase the PV consumption ().
Cluster partition-based DN planning mainly includes two aspects of cluster partition and cluster planning. In terms of cluster partition, existing research mainly establishes cluster partition indexes based on the grid structure () and the power balance in the clusters (). Cluster partition is optimized by particle swarm algorithms (), clustering algorithms (), and community detection algorithms (). By improving the community algorithm, the division of reactive and active clusters considering the power balance and node coupling degree is realized by . The gray clustering method based on the improved whitening weight function is used to partition the distribution network by , and the index weight is obtained by comprehensively applying the analytic hierarchy process and the entropy weight method. The modular index based on the electrical distance and the active power balance index are used as comprehensive division indexes by , and the distributed photovoltaic generation in the distribution network is divided into clusters by using genetic algorithms. Based on the theory of the modularity function model in complex networks, a voltage coordination control method of partitioning the aggregated domain of reactive voltage sensitivity weights and active network loss-voltage sensitivity weights of power systems is proposed by . The nodes with a strong coupling relationship are merged to determine the initial number of partitioning by , and then, the final partitioning result is determined according to the affiliation between each load node and each reactive power source. Combining the K-means clustering algorithm and optimized PSO algorithm for voltage regulation within the cluster ensures that the voltage crossing problem is solved by . A cluster partition index system considering the structural and functional properties is proposed by , and the modular index that takes into account the characteristics of electric and heat networks is used on the structural property to describe the connection strength between different network nodes. However, the existing cluster partition indexes only concern the active power balance in each cluster, ignoring the impact of reactive power. Meanwhile, the existing cluster partition indexes ignore the influence of cluster size on the planning results, which can easily lead to large differences in the cluster size, even leading to isolated nodes (). In addition, the existing cluster partition methods have insufficient computational accuracy, and for a complex cluster partitioning index, the optimization results tend to fall into local optimum solutions.
The current active distribution network (ADN) planning strategy usually includes the reinforcement or expansion of distribution networks and DG integration under the active management of DG outputs (). A two-level robust optimal feeder routing model for the planning of radial distribution networks is proposed by , where power demand is uncertain. The robust model is solved by implementing the column and constraint generation strategy. A method based on calculating the probability of electric vehicles (EVs) entering each parking lot is proposed by for the long-term planning of EV parking lots. An integrated power and gas systems of IPGS considering cascading effects for enhancing resilience is proposed by , and the two-phase framework containing phases of “demand reachability evaluation” and “integrated planning” is proposed. A framework for the optimal planning of battery swapping stations (BSSs) in centralized charging mode is proposed by , and in this mode, the batteries are charged at a central charging station. Possible equipment measures are classified into several categories by , formulating the “low-voltage system configuration determination problem;” in addition, a solution algorithm based on the practical priorities of classified measures is proposed. The resilience-oriented distribution network planning problem utilizing a novel three-stage hybrid framework is proposed by , and the decision-making on the line hardening and DG placement is carried out in the first stage. In the second stage, emergency and normal operation optimization is conducted. A collaborative stochastic expansion planning model of a cyber–physical system with resilience constraints is proposed by , and the model can reduce the coupling risk and enhance the resilience under extreme scenarios. An appropriate probabilistic wind power capacity expansion planning method for a bundled wind–thermal generation system with retrofitted coal-fired units is reformulated as a mixed-integer second-order cone programming problem by . In terms of cluster planning, the existing research mainly considers deterministic scenarios as the research background (), ignoring the impact of source-load uncertainty (). Aiming at solving the problems of resource waste caused by the large-scale access of distributed generators to distribution networks and improving the economy of energy storage systems, a cluster energy-storage control strategy for prompting the distributed generation accommodation and improving the economy of energy storage systems is proposed by . In considering the optimization of load distribution among units and introducing consumption costs, a grid evaluation index system including the coordination index of the power transmission and distribution network is constructed by . A novel cluster-based distributed generation planning approach is proposed by , and the distribution network is divided into several partitions considering the system network structure and the load characteristics, thus conducting a hierarchical and partitioned network structure. A planning model of renewable energy access is established by based on cluster partition considering the investments and power generation interests of power producers and the power match degree within clusters. With the increasing proportion of distributed PVs, the source-load uncertainty increases the difficulty in modeling the uncertainty of DNs () and increases the dimensionality of variables in the planning model (). How to establish a cluster planning model based on the source-load uncertainty () and simplify the traditional centralized planning model need to be further researched.
Based on the above analysis, this paper proposes a cluster partition-based two-layer expansion planning model of grid–resource–storage for DNs. The main contributions of this paper are summarized as follows:
(1) To deal with poor power balance and unbalanced cluster size in existing cluster partition, a comprehensive cluster partition index is proposed, which includes the modularity index, power balance index, and nodal size index. In addition, based on the comprehensive cluster partition index, an improved genetic algorithm is proposed to partition the DN into some clusters.
(2) To deal with complex models in centralized planning methods, a cluster partition-based two-layer expansion planning model is established for the DNs. In the upper layer, a line planning model is established to carry out the planning of cluster connection lines. In the lower layer, the PV and ES planning model within a cluster is established, which can realize the optimal planning of PVs and ESs in each cluster.
(3) To reduce the conservatism of the traditional robust optimization, a box uncertainty set is utilized to characterize the uncertainty of loads and PVs, and an uncertainty regulation parameter is used to control the range of uncertainty sets, which can reduce the conservatism of the optimization and simplify the calculation process.
The remainder of this paper is organized as follows: a comprehensive cluster partition index-based cluster partition method is proposed in Section 2; a cluster partition-based two-layer expansion planning method is proposed in Section 3; in Section 4, the case study is analyzed; and the conclusion is given in Section 5.
2 Comprehensive cluster partition index-based cluster partition method
2.1 Comprehensive cluster partition index
As the existing cluster partition index is not comprehensive, based on the DN structure and cluster function, a comprehensive cluster partition index is proposed to complete the cluster partition in this paper. The proposed comprehensive cluster partition index includes the modularity index, power balance index, and nodal size index.
2.1.1 Modularity index
The coupling degree between nodes can be measured by a modularity index based on voltage sensitivity, which is expressed as
where ρm is the modularity index. is the edge weight between node i and node j. S is the node set of DNs. is the sum of the edge weight of all networks. is the sum of the edge weights that are connected to node i. is the judgment function of the cluster; when node i and node j are in the same cluster, then = 1; otherwise, = 0. is determined by the electrical distance, and the electrical distance can indicate the electrical coupling degree between two nodes in the network. For a DN with N nodes, the electrical distance based on the reactive voltage sensitivity matrix is expressed as
where is the element in row i and column j of the reactive voltage sensitivity matrix, which represents the sensitivity of reactive power generation at node i to the voltage at node j. is the electrical distance between the two nodes based on the reactive voltage sensitivity matrix; the larger the value of , the smaller the electrical distance between the two nodes. Similarly, the electrical distance based on the active voltage sensitivity matrix can be obtained. The node voltages are affected by active and reactive power variations, and then, the electrical distance based on the sensitivity matrix is expressed as follows:
The relationship between the edge weight and electrical distance is that the larger the edge weight, the smaller the electrical distance. Then, the mathematical expression between the edge weight and electrical distance can be obtained as follows:
where max(Lij) is the maximum value of the elements in the electrical distance matrix.
2.1.2 Power balance index
In order to avoid large-scale power transfer between clusters, the PV output and load demand within a cluster should be as equal as possible. In order to evaluate the ability of clusters to hold the distributed PVs, this paper proposes the power balance index. The active power balance index φP is given as follows:where Nk is the number of clusters. is the net power of cluster k at time t. T is the period of optimization. Similarly, the reactive power balance index φQ is established as follows:
where Qsup is the maximum value of reactive power supply within the cluster and Qneed is the value of reactive power demand within the cluster.
2.1.3 Nodal size index
A reasonable cluster size will directly impact the complexity of the subsequent cluster planning, as well as avoid the differences in the complexity of optimization among different clusters caused by unbalanced cluster size. In addition, a reasonable cluster size can avoid the isolated nodes in the planning process. In order to balance the size of each cluster, a nodal size index is established as follows:where is the cluster that node i belongs to. is the affiliation degree of node x to cluster y. is the sum of the edge number in the cluster that node i belongs to. are the clusters that do not contain node i. is the sum of the edge number in clusters that do not contain cluster .
The indexes shown in Eq 1, (7), (8), and (9) are combined into a comprehensive cluster partition index for the DNs, which is expressed as follows:where ω1, ω2, ω3, and ω4 are the weights of each index. In the process of cluster partition, different weights can be set for each index depending on different needs.
2.2 Improved genetic algorithm-based cluster partition method
To carry out the cluster partition, a hybrid genetic-simulated annealing (HGSA) algorithm is utilized. The HGSA algorithm uses the annealing selection as the individual replacement strategy, while the global information obtained by the genetic algorithm can be completely used, and the premature convergence of the genetic algorithm can be avoided. Then, the global convergence of the algorithm is enhanced. The cluster partition is implemented as follows:
Step 1: Initial optimization parameters are set: population size n, initial temperature T0, termination temperature Tend, temperature cooling factor r, maximum number of genetic generations M, and objective function for the cluster partitioning index .
Step 2: The number of temperature updates is set equal to 0. Considering that the cluster partition is carried out based on the original network connectivity, this paper uses the unweighted adjacency matrix to represent the connectivity of the network, where indicates that nodes i and j are connected and indicates that nodes i and j are not connected. The clustering partition is randomly modified when = 1. The result of the modification is 0 or 1, and 0 means that the two nodes are disconnected from each other, and then, the two nodes belong to different clusters. In the iteration of the algorithm, in order to ensure that the result after crossover still satisfies the meaning of the adjacency matrix, only the upper triangle of the adjacency matrix is chosen to carry out the crossover. After the crossover step, the upper triangle of the matrix is symmetrically transferred to the lower triangle to form the newborn individuals, thus generating the initial population Pl(k).
Step 3: The individuals that satisfy the cluster constraints are screened. Considering that the reverse power flow only occurs within the cluster, the net power within each cluster needs to be greater than 0 at each moment among the cluster partition. For the individuals who cannot satisfy the constraint, the population eliminates these individuals.
Step 4: Eq. (12) is chosen as an indicator to calculate the fitness of individuals. Considering that the value of fitness intuitively reflects the superiority or inferiority of cluster partition, the individuals with greater fitness are replicated to the offspring to form the new population .
Step 5: The crossover and mutation are carried out for using the traditional genetic algorithm to obtain the population . The mutated populations are limited to accepting bad solutions by the simulated annealing algorithm to form new populations .
Step 6: is set. If the genetic algebra accumulates to the maximum number M, then step 7 is repeated; otherwise, step 3 is repeated.
Step 7: The temperature is updated, i.e., Tl=rT0, k = 0, . If the convergence condition Tl< Tend, then the computation to output the optimal solution is terminated. Otherwise, the temperature reduction operation is performed, i.e., Tl+1= rTl, and step 3 is repeated.
3 Cluster partition-based two-layer expansion planning method
In this paper, a cluster partition-based two-layer expansion planning model is proposed, which involves the uncertainty of PVs and loads. In the upper layer, a line planning model is established with the objective of minimizing the line investment and network loss costs. In the lower layer, a source-storage planning model is proposed for PVs and ESs with the objective function of minimizing source-storage investment and operation costs within a cluster. Meanwhile, the box-type uncertainty set is utilized to characterize the uncertainty of PVs and loads in the lower layer, and an uncertainty parameter is used to control the range of uncertainty sets, which can reduce the conservatism of the optimization.
3.1 Upper-layer line planning model
In the upper layer, the objective function of the line planning model is established as follows:
where
is the node set of the branch network.
is the child node set of node
i.
CLis the investment cost parameter of lines.
Dijis the line length of line
i–j.
xijis the 0–1 variable, where
xij= 1 means that the line is installed and
xij= 0 means that the line is not installed.
is the bank rate.
is the payback period.
Ceis the electricity price.
Nn,kis the node numbers of cluster
k.
is the resistance of the line
i–
j.
iij,tis the squared value of the current in the line
i–
jat time
t.The constraints include a power flow constraint and a penetration rate constraint of PVs in the clusters.
(1) Power flow constraint
The power flow containing line variables is constrained by second-order conic relaxation (SOCR) (
) as follows:
where
ψjis the upstream node set of node
j.is the reactance of line
i–
j.
Pij,tand
Qij,tare the active and reactive power through line
i–
jat time
t, respectively.
is the ES charging and discharging efficiency.
and
are the charging and discharging power of ES at node
jat time
t, respectively.
and
are the active and reactive power of loads at node
jat time
t, respectively.
is the actual active power generated by the PV at node
jat time
t.
is the additional active power required to be generated by the PV at node
jat time
t.
φjis the downstream node set of node
j.is the actual reactive power generated by the PV inverter at node
jat time
t.
is the additional reactive power required to be generated by the PV inverter at node
jat time
t. vi,tis the square of the voltage at node
iat time
t.is the square of the voltage at node
ithat is constrained by the line variable
xij.W is a large constant.
vminand
vmaxare the upper and lower voltage limits, respectively.
Imaxis the maximum current limit.
(2) Penetration rate constraint of PVs
Many isolated nodes exist in the DN that need to be connected to the system, and the penetration rate of PVs should be constrained when carrying out line planning to access these isolated nodes, which is defined as follows:
where is the planed capacity of the PV of node j in the cluster and is the load apparent power of node j in the cluster.
3.2 Lower-layer PV and ES planning model
In the lower layer, the PV and ES planning model is established to identify the “worst scenario " of the uncertain variables, and based on the worst scenario, the proposed model minimizes the investment and operation costs of source storage within the cluster. Therefore, the objective function for the PV and ES planning model within the cluster is established as
where
F1is the annual investment and operation costs of PVs in cluster
k.
F2is the annual revenue of PVs in cluster
k.
F3is the annual investment costs of ESs in cluster
k.
CPVis the investment cost parameter of PVs.
COMPVis the annual fixed maintenance cost parameter of PVs.
Пkis the node set of cluster
k.
is the original installed PV capacity of node
jin cluster
k.
is the planned PV capacity of node
jin cluster
k.
Cselland
Csubare the feed-in tariff and subsidized tariff of the PVs, respectively.
is the output limit of PVs per megawatt for a given sunshine condition at time
t.is the discount rate.
ris the discounted number of years.
cbatis the investment cost parameter of ESs.
Pbatt,jis the allocated capacity of ES at node
j.The constraints include uncertainty constraints for PVs and loads, PV and ES capacity constraints, and power flow constraints.
(1) Uncertainty constraints for PVs and loads
The box uncertainty set is utilized to characterize the uncertainty range of active load power, reactive load power, and PV outputs. Meanwhile, an uncertainty parameter is utilized that can be set to adjust the conservativeness of the optimal solution. The larger the uncertainty parameter, the more conservative the solution. The specific formula is as follows:
where
Urepresents the boxed uncertainty set.
,
, and
are the predicted values of active load power, reactive load power, and PV outputs of node
jat time
tin cluster
k,respectively.
,
, and
are the fluctuation deviations of active load power, reactive load power, and PV output of node
jat time
tin cluster
k, respectively.
Na,Nb,and
Neare the total number of active load nodes, the total number of reactive load nodes, and the total number of PV nodes in cluster
k, respectively.
ГPL,
ГQL, and
ГPVare the uncertain adjustment parameters for active load power, reactive load power, and PV output S, respectively, which are integers from 0 to
Na,0 to
Nb, and 0 to
Ne,respectively. The decision-maker can choose a variety of uncertainty regulation parameters to adjust the scheme flexibly; the larger the value of each uncertainty regulation parameter, the more conservative the resulting planning scheme.
(2) PV and ES constraints
where
is the simulated operating time of the ESs in cluster
k.
is the annual PV generation time.
socj,tis the state of charge of ES at node
jin cluster
kat time
t.
socminand
socmaxare the minimum and maximum charge states of ES, respectively.
PCH,maxand
PDS,maxare the maximum power limits for charging and discharging power of ES, respectively.
is the maximum load demand of node
jin the cluster.
(3) Power flow constraints
Eqs (14), (15), and (17) are referred to for the power flow constraints.
3.3 Iterative solving process
In this paper, an iterative solution method is utilized to solve the proposed planning models. For the upper layer, the planning model can be expressed by a specific form as follows:
where x is the optimization vectors. c is the coefficient matrices corresponding to the objective functions. P, k, B, C, Rx, S* u, G, and H are the coefficient matrices corresponding to the variables under the constraints. α, D, O, V, and W are the constant column vectors. l is the current number of iterations. k is the maximum iteration. yl is the variable at the lth iteration. u* l is the value of the uncertain variable u in the “worst” scenario after the lth iteration.
The specific form of the lower-layer model is
where y is the optimization vectors. is the feasible domain of y under the given x and u. x* is the optimization vectors obtained from the upper layer. d is the coefficient matrices corresponding to the objective functions. Su and R* u are the coefficient matrices corresponding to the variables under the constraints. Given a set of u, the inner min. problem becomes a second-order cone programming problem, which can be transformed into a “max” form by the duality theory:where[γ,ν,π,μ,μ1 ] are the dual variables in the lower-layer model.
After the above transformation, the proposed model can be solved by the iterative solution method as follows:
Step 1: Given a set of u values as the initial worst-case scenario, a lower bound is set on the operating cost LB =, an upper bound UB = , and the number of iterations l = 1.
Step 2: The upper model is solved based on the worst-case scenario ul to obtain the optimal solution x* l and α* l. The value of α* l is used as the new lower bound LB = max(LB, α* l).
Step 3: The lower layer is optimized based on the optimization results of the upper layer, and the optimized results fl (x* l) and the worst-case scenario x* l are obtained. The upper bound is updated as UB = min(UB, fl (x* l)).
Step 4: If UB-LB < ε, where ε is a threshold of convergence, then, the optimal solutions can be obtained, and the iteration is stopped. Otherwise, the variable yl+1 and the following constraints are added:
Let l = l + 1, and step 2 is repeated until the algorithm converges.
4 Case study
4.1 Case study system
In order to verify the effectiveness of the proposed method, an actual 35 kV/10 kV DN in China is utilized for analysis. The total loads in this DN are 7.3 MVA, and the total capacity of PVs in this DN is 4 MW. The installed PV capacity is shown in Table 1. The topology of this system is shown in Figure 1, the lines to be planned are shown in Table 2, and the line parameters of the network are shown in Table 3.
TABLE 1
| Node location | Installed capacity (kW) | Node location | Installed capacity (kW) | Node location | Installed capacity (kW) |
|---|---|---|---|---|---|
| 7 | 500 | 16 | 500 | 21 | 500 |
| 9 | 500 | 17 | 500 | 27 | 1,500 |
Installed capacity of photovoltaics (PVs).
FIGURE 1
TABLE 2
| Lines | Distance/km | Line | Distance/km | Line | Distance/km |
|---|---|---|---|---|---|
| 7–11 | 15.64 | 30–27 | 17.62 | 16–27 | 23.64 |
| 9–11 | 22.37 | 14–21 | 6.50 | 23–21 | 5.30 |
Lines to be planned.
TABLE 3
| Start node | Final node | Line impedance (Ω/km) | Start node | Final node | Line impedance (Ω/km) |
|---|---|---|---|---|---|
| 1 | 2 | 0.192 + 20.533 j | 4 | 29 | 0.515 + 0 j |
| 1 | 3 | 0.192 + 0 j | 4 | 30 | 0.13 + 0 j |
| 1 | 7 | 0.192 + 13.200 j | 5 | 22 | 0.13 + 6.482 j |
| 1 | 8 | 0.61 + 1.820 j | 7 | 11 | 0.299 + 0.823 j |
| 1 | 20 | 4.681 + 8.314 j | 9 | 11 | 0.12 + 0.471 j |
| 1 | 32 | 7.686 + 14.750 j | 14 | 21 | 0.186 + 0.261 j |
| 2 | 6 | 1.084 + 3.236 j | 15 | 16 | 0.21 + 0.730 j |
| 2 | 9 | 2.132 + 3.786 j | 16 | 27 | 0.9 + 0.804 j |
| 2 | 10 | 1.703 + 4.690 j | 16 | 31 | 0.731 + 1.298 j |
| 2 | 12 | 0.13 + 10.323 j | 23 | 21 | 3.773 + 5.297 j |
| 2 | 13 | 0.238 + 16.259 j | 24 | 29 | 1.24 + 2.662 j |
| 2 | 18 | 0.238 + 0 j | 28 | 29 | 2.388 + 4.242 j |
| 3 | 4 | 5.628 + 16.797 j | 28 | 30 | 0.515 + 29.563 j |
| 3 | 5 | 5.616 + 9.975 j | 30 | 27 | 2.495 + 4.432 j |
| 3 | 14 | 0.238 + 10.209 j | 33 | 29 | 4.231 + 3.621 j |
| 3 | 15 | 0.476 + 0 j | 34 | 30 | 3.681 + 9.235 j |
| 3 | 17 | 4.456 + 7.915 j | 35 | 3 | 7.214 + 3.751 j |
| 3 | 19 | 2.7 + 3.790 j | 4 | 28 | 4.396 + 6.172 j |
| 3 | 24 | 0.64 + 1.290 j |
Line parameters of the network.
The simulated genetic annealing algorithm used in this paper sets the population size n = 40, the initial temperature T0 = 100°C, the termination temperature Tend = 1°C, the temperature cooling factor r = 0.8, the maximum number of genetic generations N = 500, the crossover probability pc = 0.4, and the variation probability pm = 0.2. This paper uses MATLAB version 2019 to program the algorithm.
4.2 Analysis of cluster partition
In order to illustrate the superiority of the proposed cluster partition index in this paper, the modularity function used by is selected to compare with the proposed cluster partition index. The weights of the proposed cluster partition index in this paper are taken as ω1 = ω2 = ω3 = ω4 = 0.25. The cluster partition results under the proposed cluster partition index are shown in Figure 2, and the cluster partitioning results under the modularity function are shown in Figure 3.
FIGURE 2
FIGURE 3
Figure 2 and Figure 3 show that the size difference among sub-networks attained by the proposed cluster partition index is smaller than that attained by the modularity function, and isolated nodes forming separate clusters (cluster 3 and cluster 6) occur in modularity function. In addition, there are PVs in each cluster attained by the proposed method, while no PVs exist in cluster 1 attained by the modularity function. The differences between the two methods prove that the cluster partition index can result in a more reasonable cluster partition for the subsequent cluster control and operation.
In order to illustrate the superiority of the proposed cluster partitioning algorithm, the traditional genetic algorithm is selected to be compared. The comparison of cluster partition results obtained by the two algorithms is shown in Table 4. Table 4 shows that the optimization results of each cluster partitioning index of the proposed algorithm are greater than those of the traditional genetic algorithm, which indicates that the proposed cluster partitioning algorithm has better optimization performance.
TABLE 4
| Result | ||||
|---|---|---|---|---|
| Proposed algorithm | 0.6953 | 0.8639 | 0.6866 | 0.7233 |
| Traditional genetic algorithms | 0.6132 | 0.8071 | 0.6120 | 0.6982 |
Comparison of cluster partition results obtained by different algorithms.
Different combinations of partition index weights are given in Table 5. As shown in Table 5, when increasing the weight of the modularity index, the nodes within the cluster are better connected, but the cluster power complementarity decreases. When increasing the weight of the active and reactive power balance index, the cluster power complementarity characteristics improve, but the nodes within the cluster are significantly less connected. When increasing the weight of the nodal size index, the number of clusters changes. When increasing the weight of the nodal size index, the number of clusters changes, and the cluster partition is more focused on the change in cluster size.
TABLE 5
| Combination | Number of clusters | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| 1 | 0.1 | 0.4 | 0.4 | 0.1 | 0.632 1 | 0.896 4 | 0.711 7 | 0.701 6 | 7 |
| 2 | 0.2 | 0.3 | 0.3 | 0.2 | 0.677 9 | 0.889 1 | 0.679 3 | 0.708 3 | 6 |
| 3 | 0.3 | 0.2 | 0.2 | 0.3 | 0.742 4 | 0.833 0 | 0.657 1 | 0.739 3 | 6 |
| 4 | 0.4 | 0.1 | 0.1 | 0.4 | 0.789 1 | 0.798 3 | 0.625 3 | 0.756 1 | 5 |
Influence of different weights on cluster partition.
4.3 Analysis of the proposed planning strategy
In order to verify the flexibility of the proposed two-layer expansion planning model, three sets of uncertain regulation parameters are selected and compared in terms of network connection, planned capacity of PVs and ESs, and total planning costs. Three sets of uncertain regulation parameters are set as follows:
Case 1:ГPL = ГQL = ГPV = 0 is set, and in this case, the active load power, reactive load power, and PV output are all equal to the predicted values.
Case 2:ГPL = ГQL = 16 and ГPV = 4 are set, and in this case, the active and reactive load power of 16 nodes is taken to the minimum value of the prediction interval, while the PV outputs of 4 nodes are taken to the maximum value of the prediction interval.
Case 3:ГPL = ГQL = 27 and ГPV = 6 are set. In this case, the active and reactive load power of 27 nodes is taken to the minimum value of the prediction interval, while the PV outputs of 4 nodes are taken to the maximum value of the prediction interval, which is the worst scenario.The comparisons of the three cases in planned PV capacity and ES capacity are shown in Figure 4 and Figure 4, and the planning results and total planning costs are shown in Table 6 and Table 7, respectively.Figure 4 shows that with the increase in uncertainty regulation parameters, the planned PV capacity decreases. The reason is that the scenario becomes increasingly severe, and to maintain the safe operation of the cluster, the corresponding PVs are reduced.As shown in Table 6, with the increase in uncertainty regulation parameters, the connection of nodes 21 and 27 is different under the three cases. The reason is that as the operating scenario becomes more severe, and the load power increases, the planning results tend to favor a reduction in the power supply range.As shown in Figure 5, as the conservatism of the planning increases, the ES installed capacity also increases. The reason is that the active power and reactive power of the load in cases 2 and 3 are smaller than those of case 1, and the PV outputs are larger; hence, more ESs are needed to mitigate the uncertainty of PVs and loads.In addition, Table 7 shows that as the uncertainty regulation parameter increases, although the annual investment and operation costs of PVs correspondingly decrease, the annual income of PVs decreases, and the annual investment costs of networks, network losses, and the annual investment costs of ESs correspondingly increase. The reason is that the more uncertainty of DNs considered in the planning process, the more conservative the resulting solution becomes, and the corresponding total cost increases.In order to further illustrate the superiority of the proposed planning method, the centralized planning method is selected to compare with the proposed planning method. The uncertain regulation parameters in this comparison are set as ГPL = ГQL = 12 and ГPV = 3, and in this case, the active and reactive load power of 12 nodes is taken to the minimum value of the prediction interval, while the PV outputs of 3 nodes are taken to the maximum value of the prediction interval. The network connection of the two methods is shown in Figure 6 and Figure 7. The comparison of the calculation time and total planning cost of the two methods is shown in Table 6.Figure 6 and Figure 7 show that the connection of node 11 and node 27 under the proposed method is different from the centralized planning scheme. The length of the power supply lines for node 11 and node 27 under the proposed method has been reduced by 6.73 km and 6.02 km, respectively; the power supply range is greatly shortened, which can effectively reduce the line investment costs and network loss costs, as well as improve the economic efficiency of the planning scheme.The calculation method for the difference rate given in Table 8 is to subtract the results of the proposed method from those of the centralized planning method and divide them by the results of the centralized planning method. As shown in Table 8, in the proposed method, the annual PV investment costs and PV operation costs, as well as the annual income of PVs, increased by 5.3% and 14.1%, respectively, compared to the centralized planning method, and the annual investment costs of lines, network losses, and annual investment costs of ESs decreased by 13.5%, 25.7%, and 15.5%, respectively. These different rates indicate that the proposed method can improve the overall economic efficiency of the planning, as well as reduce the costs of planning. In addition, the calculation time of the proposed method is 20.41 s, but that of the centralized planning method is 43.62 s. Owing to the fact that the centralized planning method requires a whole optimization process, the huge variables complex the optimization. However, the proposed method can partition the complex network into smaller sub-networks, and the complex problem is divided into several relatively simple sub-problems that can be solved rapidly and easily by a parallel computing way. Therefore, the computing time can be greatly shortened. The optimization speed of the proposed method is improved by 53.2% compared to centralized planning methods, which indicates that proposed method is more suitable for DN planning with large-scale PVs.
FIGURE 4
TABLE 6
| Case | Line connection results | Cluster planning results | |
|---|---|---|---|
| Connectivity | Disconnection | ||
| 1 | 7–11 | 9–11 | Cluster 1 (7,8,11,20,32); cluster 2 (2,6,9,10,12,13,18) |
| 23–21 | 14–21 | Cluster 3 (3,14,17,19,24,26,35); cluster 4 (5,21,22,23,25) | |
| 30–27 | 16–27 | Cluster 5 (15,16,31); cluster 6 (4,27,28,29,30,33,34) | |
| 2 | 7–11 | 9–11 | Cluster 1 (7,8,11,20,32); cluster 2 (2,6,9,10, 12,13,18) |
| 14–21 | 23–21 | Cluster 3 (3,14,17,19,21,24,26,35); cluster 4 (5, 22,23,25) | |
| 30–27 | 16–27 | Cluster 5 (15,16,31); cluster 6 (4,27,28,29,30,33,34) | |
| 3 | 7–11 | 9–11 | Cluster 1 (7,8,11,20,32); cluster 2 (2,6,9,10,12,13,18) |
| 14–21 | 23–21 | Cluster 3 (3,14,17,19,21,24,26,35); cluster 4 (5,22,23,25) | |
| 16–27 | 30–27 | Cluster 5 (15,16,27,31); cluster 6 (4,28,29,30,33,34) | |
Planning results of different schemes.
TABLE 7
| Comparison item | Case 1 | Case 2 | Case 3 |
|---|---|---|---|
| Annual investment costs of photovoltaics (PVs) | 587.54 $ | 553.48 $ | 528.22 $ |
| Annual operation costs of PVs | 414.93 $ | 232.90 $ | 187.04 $ |
| Annual income of PVs | 240.43 $ | 247.92 $ | 285.45 $ |
| Annual investment costs of lines | 622.98 $ | 856.43 $ | 914.28 $ |
| Network losses | 139.82 $ | 154.68 $ | 185.82 $ |
| Annual investment costs of energy storages (ESs) | 1,062.34 $ | 1,579.61 $ | 1,726.73 $ |
Planning costs in different cases.
FIGURE 5
FIGURE 6
FIGURE 7
TABLE 8
| Comparison item | Proposed method | Centralized planning method | Difference rate (%) |
|---|---|---|---|
| Annual investment costs of photovoltaics (PVs) | 569.86 $ | 541.11 $ | −5.3 |
| Annual operation costs of PVs | 293.10 $ | 239.87 $ | −14.1 |
| Annual income of PVs | 240.43 $ | 277.97 $ | 13.5 |
| Annual investment costs of lines | 692.89 $ | 932.30 $ | 25.7 |
| Network losses | 150.14 $ | 177.58 $ | 15.5 |
| Annual investment costs of energy storages (ESs) | 1,360.22 $ | 1,689.09 $ | 10.3 |
| Calculation time | 20.41 s | 43.62 s | 53.2 |
Comparison between the proposed method and centralized planning method.
5 Conclusion
In this paper, a cluster partition-based two-layer expansion planning for DNs is proposed, and an actual 35 kV/10 kV DN in China is utilized for analysis. The results show that
1) To deal with poor power balance and unbalanced cluster size in existing cluster partition, a comprehensive cluster partition index is proposed in this paper, which includes the modularity index, power balance index, and nodal size index, and in order to avoid the cluster partition falling into local optimum, an improved genetic algorithm is utilized to carry out the network partition. The differences in cluster size among the clusters attained by the proposed cluster partition are smaller, and every cluster has PVs.
2) To deal with complex models in centralized planning methods, a cluster partition-based two-layer expansion planning model is established for the DNs, which decomposes the complex centralized planning model into cluster planning. The proposed method can improve the overall economic efficiency of the planning, as well as reduce the costs of planning. Meanwhile, the optimization speed is also improved.
3) To reduce the conservatism of traditional robust optimization, a box uncertainty set is utilized to characterize the uncertainty of loads and PVs, and an uncertainty regulation parameter is used to control the range of uncertainty sets, which can reduce the conservatism of the optimization, as well as simplify the calculation process.
In the actual operation, the PV outputs are uncontrollable and may exceed the worst scenario set in the planning stage. Then, the PV outputs within the cluster may not be fully absorbed. In future research, a cluster partition-based scheduling method for DNs should be further researched to enhance the source-load complementarity of the DNs.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
SY: writing–original draft and writing–review and editing. SS: writing–original draft. YC: writing–review and editing. PY: writing–review and editing. CW: writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the project supported by the Science and Technology Project of State Grid Shandong Electric Power Company “Research and application of access and control technology targeting the distributed photovoltaics developed by the entire county (city, district)” (520626220014).
Conflict of interest
The authors declare that this study received funding from State Grid Shandong Electric Power Company. The funder had the following involvement in the study: The funder was involved in the study design.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
BiR.LiuX.DingM.et al (2019). A method of renewable energy power generation cluster classification with the goal of improving the consumption capacity. Chin. J. Electr. Eng.39 (22), 6583–6592. 10.13334/j.0258-8013.pcsee.182512
2
CaiZ.ShuangH.ShanJ. (2022). Optimal scheduling strategy for electric bus cluster power exchange considering operation cost. Power Syst. Autom.46 (17), 205–217. 10.7500/AEPS20210916007
3
DingM.FangH.BiR.et al (2019). Distributed photovoltaic and energy storage siting and capacity planning for distribution networks based on cluster delineation. Chin. J. Electr. Eng.39 (8), 2187–2201. 10.13334/j.0258-8013.pcsee.180757
4
FaramarziD.RastegarH.RiahyG. H.Doagou-MojarradH. (2023). A three-stage hybrid stochastic/IGDT framework for resilience-oriented distribution network planning. Int. J. Electr. Power Energy Syst.146 (5), 108738–108738.15. 10.1016/j.ijepes.2022.108738
5
GeJ.LiyY.PangD.et al (2024). Cluster division voltage control strategy of photovoltaic distribution network with high permeability. High. Volt. Eng.50 (1), 74–82. 10.13336/j.1003-6520.hve.20230816
6
Haji-AghajaniE.HasanzadehS.Heydarian-ForushaniE. (2023). A novel framework for planning of EV parking lots in distribution networks with high PV penetration. Electr. Power Syst. Res.217 (4), 109156–109156.12. 10.1016/j.epsr.2023.109156
7
HuD.DingM.BiR.et al (2020). Analysis of the impact of photovoltaic and wind power complementarity on the planning of high penetration renewable energy cluster access. Chin. J. Electr. Eng.40 (3), 821–836. 10.13334/j.0258-8013.pcsee.182308
8
HuD.SunL.DingM. (2023). Integrated planning of an active distribution network and DG integration in clusters considering a novel formulation for reliability assessment. CSEE J. Power Energy Syst.9 (2), 561–576. 10.17775/CSEEJPES.2020.00150
9
JiY.ChenX.HeP.LiuX.WuX.ZhaoC. (2023). A novel voltage/var sensitivity calculation method to partition the distribution network containing renewable energy. Recent Adv. Electr. Electron. Eng.16 (4), 380–394. 10.2174/2352096516666221130150549
10
JiangY.RenZ.LiQ.et al (2022). New energy consumption strategy in distribution networks considering coordinated scheduling of multi-flexibility resources. J. Electr. Eng. Technol.37 (7), 1820–1835. 10.19595/j.cnki.1000-6753.tces.211464
11
KongX.LiuC.ChenS.et al (2022). A multi-time-node response potential assessment method for adjustable resource clusters considering dynamic processes. Power Syst. Autom.46 (18), 55–64. 10.7500/AEPS20220501001
12
KoutsoukisN.GeorgilakisP.HatziargyriouN. (2018). Multistage coordinated planning of active distribution networks. IEEE Trans. Power Syst.33 (1), 32–44. 10.1109/tpwrs.2017.2699696
13
LeiC.WangQ.ZhouG.BuS.LinT.et al (2023). Probabilistic wind power expansion planning of bundled wind-thermal generation system with retrofitted coal-fired plants using load transfer optimization. Int. J. Electr. Power Energy Syst.151 (9), 109145–109145.18. 10.1016/j.ijepes.2023.109145
14
LiC.DongZ.LiJ.et al (2021). Optimal control strategy for distributed energy storage clusters to enhance new energy consumption capacity of distribution networks. Power Syst. Autom.45 (23), 76–83. 10.7500/AEPS20210224004
15
LiP.ShenJ.WuZ.YinM.DongY.HanJ. (2023). Optimal real-time Voltage/Var control for distribution network: droop-control based multi-agent deep reinforcement learning. Int. J. Electr. Power Energy Syst.153 (11), 109370–109370.11. 10.1016/j.ijepes.2023.109370
16
LiS.YangX. (2022). Siting and capacity determination of photovoltaic energy storage based on bi-directional dynamic reconfiguration and cluster delineation. Power Syst. Prot. Control50 (3), 51–58. 10.19783/j.cnki.pspc.210385
17
LiY.LuN.LiuX.et al (2022c). Output evaluation method of distributed photovoltaic cluster considering renewable energy accommodation and power loss of network. Electr. Power Constr.43 (10), 136–146. 10.12204/j.issn.1000-7229.2022.10.013
18
LiY.YaoT.QiaoX.et al (2022b). Optimal allocation of distributed PV and energy storage based on joint timing scenarios and source-grid-load synergy. J. Electr. Eng. Technol.37 (13), 3289–3303. 10.19595/j.cnki.1000-6753.tces.210712
19
LiZ.WangW.HanS.et al (2022a). Study on voltage adaptation of distributed photovoltaic access to distribution network considering reactive power support. Power Syst. Prot. Control50 (11), 32–41. 10.19783/j.cnki.pspc.211781
20
LiuR.ShengW.MaX.et al (2021). Research on the planning of distributed photovoltaic access to distribution networks based on multiple swarm genetic algorithms. J. Sol. Energy42 (6), 146–155. 10.19912/j.0254-0096.tynxb.2019-0370
21
LiuZ.ZhangT.WangY. (2022). Model-based predictive control for multi-scenario variable time-scale optimal scheduling of active distribution network. Power Autom. Equip.42 (4), 121–128. 10.16081/j.epae.202201001
22
MukherjeeB.SossanF. (2023). Optimal Planning of single-port and multi-port charging stations for electric vehicles in medium voltage distribution networks. IEEE Trans. Smart Grid14 (2), 1135–1147. 10.1109/tsg.2022.3204150
23
PanM.LiuN.LeiJ. (2021). A dynamic classification method for distributed energy clusters containing cogeneration units. Power Syst. Autom.45 (1), 168–176. 10.7500/AEPS20200217013
24
SasakiY.FukubaS.YokotaH.YorinoN.FukubaS.et al (2023). Comprehensive solution method considering load distribution for determining low-voltage network configurations. IEEE Trans. Electr. Electron. Eng.18 (6), 1023–1032. 10.1002/tee.23808
25
ShakerM.FarzinH.MashhourE. (2023). Joint planning of electric vehicle battery swapping stations and distribution grid with centralized charging. J. Energy Storage58 (2), 106455–106455.11. 10.1016/j.est.2022.106455
26
ShenX.ShahidehpourM.ZhuS.HanY.ZhengJ. (2018). Multi-stage planning of active distribution networks considering the co-optimization of operation strategies. IEEE Trans. Smart Grid9 (2), 1425–1433. 10.1109/tsg.2016.2591586
27
SuS.LeiJ.YanY.PanS.YangY.et al (2023). Voltage regulation strategy of distribution network with decentralized wind power based on cluster partition. Recent Adv. Electr. Electron. Eng.16 (1), 30–44. 10.2174/2352096515666220902125455
28
WangL.ZhangF.KouL.et al (2021). Scaled distributed PV power cluster classification based on Fast Unfolding clustering algorithm. J. Sol. Energy42 (10), 29–34. 10.19912/j.0254-0096.tynxb.2018-0896
29
WangY.YangY.XuQ. (2023b). Integrated planning of natural gas and electricity distribution systems for enhancing resilience. Int. J. Electr. Power Energy Syst.151 (9), 109103–109103.11. 10.1016/j.ijepes.2023.109103
30
WangZ.TanW.LiH.GeJ.WangW. (2023a). A voltage coordination control strategy based on the reactive power-active network loss partitioned aggregation domain. Int. J. Electr. Power Energy Syst.144 (1), 108585–108585.10. 10.1016/j.ijepes.2022.108585
31
WuH.SunL.XiangS.et al (2022a). Active distribution network expansion planning with high-dimensional time-series correlation of renewable energy and load. Power Syst. Autom.46 (16), 40–51. 10.7500/AEPS20220104004
32
WuL.XinJ.WangC. (2022b). Installed capacity forecasting of distributed photovoltaic taking into account user herding psychology. Power Syst. Autom.46 (14), 83–92. 10.7500/AEPS20210902006
33
XiaoC.ZhaoB.ZhouJ.et al (2017). Voltage control of high proportional distributed PV clusters in distribution networks based on network partitioning. Power Syst. Autom.41 (21), 147–155. 10.7500/AEPS20170101002
34
XuL.YangJ.XiongZ. (2021a). A method of power distribution network classification based on improved grey cluster. Electr. Eng. (13), 81–84. 10.19768/j.cnki.dgjs.2021.13.021
35
XuX.ZhengX.WangS.et al (2021b). Coordinated planning method for transmission and distribution networks based on improved genetic annealing algorithm. Power Syst. Prot. Control49 (15), 1 24–131. 10.19783/j.cnki.pspc.201236
36
YangJ.HaoJ.BoZ. (2017). Hierarchical control strategy for reactive voltage of wind farm clusters based on adjacency empirical particle swarm algorithm. Power Grid Technol.41 (6), 1823–1829. 10.13335/j.1000-3673.pst.2016.2237
37
ZdraveskiV.VuleticJ.AngelovJ.TodorovskiM. (2023). Radial distribution network planning under uncertainty by implementing robust optimization. Int. J. Electr. Power Energy Syst.149 (7), 109043–109043.12. 10.1016/j.ijepes.2023.109043
38
ZhangT.XieM.WangP.et al (2021). System dynamics simulation of the shared value of distributed photovoltaic and its impact on the distribution network. Power Syst. Autom.45 (18), 35–44. 10.7500/AEPS20200901004
39
ZhangY.LiC.WanH.ShiQ.LiuW.XuY. (2023). Collaborative stochastic expansion planning of cyber-physical system considering extreme scenarios. IET Generation, Transm. Distribution17 (10), 2419–2434. 10.1049/gtd2.12819
40
ZhuJ.GuW.ZhangH.et al (2018). An optimization method for siting and capacity determination of distributed power supplies considering dynamic reconfiguration of networks. Power Syst. Autom.42 (5), 111–119. 10.7500/APES20170605016
Summary
Keywords
cluster partition, distribution network, expansion planning, distributed photovoltaics, energy storages
Citation
Yang S, Wang C, Sun S, Cheng Y and Yu P (2024) Cluster partition-based two-layer expansion planning of grid–resource–storage for distribution networks. Front. Energy Res. 12:1390073. doi: 10.3389/fenrg.2024.1390073
Received
22 February 2024
Accepted
08 April 2024
Published
07 May 2024
Volume
12 - 2024
Edited by
Minglei Bao, Zhejiang University, China
Updates
Copyright
© 2024 Yang, Wang, Sun, Cheng and Yu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Song Yang, yangsong19971@163.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.