ORIGINAL RESEARCH article

Front. Energy Res., 07 May 2024

Sec. Process and Energy Systems Engineering

Volume 12 - 2024 | https://doi.org/10.3389/fenrg.2024.1390073

Cluster partition-based two-layer expansion planning of grid–resource–storage for distribution networks

  • SY

    Song Yang *

  • CW

    Chenglong Wang

  • SS

    Shumin Sun

  • YC

    Yan Cheng

  • PY

    Peng Yu

  • State Grid Shandong Electric Power Research Institute, Jinan, China

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

.

CL

is the investment cost parameter of lines.

Dij

is the line length of line

i–j

.

xij

is 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.

Ce

is the electricity price.

Nn,k

is the node numbers of cluster

k

.

is the resistance of the line

i

j

.

iij,t

is the squared value of the current in the line

i

j

at 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

ψj

is the upstream node set of node

j.

is the reactance of line

i

j

.

Pij,t

and

Qij,t

are the active and reactive power through line

i

j

at time

t

, respectively.

is the ES charging and discharging efficiency.

and

are the charging and discharging power of ES at node

j

at time

t

, respectively.

and

are the active and reactive power of loads at node

j

at time

t

, respectively.

is the actual active power generated by the PV at node

j

at time

t

.

is the additional active power required to be generated by the PV at node

j

at time

t

.

φj

is the downstream node set of node

j.

is the actual reactive power generated by the PV inverter at node

j

at time

t

.

is the additional reactive power required to be generated by the PV inverter at node

j

at time

t. vi,t

is the square of the voltage at node

i

at time

t.

is the square of the voltage at node

i

that is constrained by the line variable

xij.

W is a large constant.

vmin

and

vmax

are the upper and lower voltage limits, respectively.

Imax

is 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

F1

is the annual investment and operation costs of PVs in cluster

k

.

F2

is the annual revenue of PVs in cluster

k

.

F3

is the annual investment costs of ESs in cluster

k

.

CPV

is the investment cost parameter of PVs.

COMPV

is the annual fixed maintenance cost parameter of PVs.

Пk

is the node set of cluster

k

.

is the original installed PV capacity of node

j

in cluster

k

.

is the planned PV capacity of node

j

in cluster

k

.

Csell

and

Csub

are 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.

r

is the discounted number of years.

cbat

is the investment cost parameter of ESs.

Pbatt,j

is 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

U

represents the boxed uncertainty set.

,

, and

are the predicted values of active load power, reactive load power, and PV outputs of node

j

at time

t

in cluster

k,

respectively.

,

, and

are the fluctuation deviations of active load power, reactive load power, and PV output of node

j

at time

t

in cluster

k

, respectively.

Na,Nb,

and

Ne

are 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

ГPV

are 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,t

is the state of charge of ES at node

j

in cluster

k

at time

t

.

socmin

and

socmax

are the minimum and maximum charge states of ES, respectively.

PCH,max

and

PDS,max

are the maximum power limits for charging and discharging power of ES, respectively.

is the maximum load demand of node

j

in 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 locationInstalled capacity (kW)Node locationInstalled capacity (kW)Node locationInstalled capacity (kW)
75001650021500
950017500271,500

Installed capacity of photovoltaics (PVs).

FIGURE 1

TABLE 2

LinesDistance/kmLineDistance/kmLineDistance/km
7–1115.6430–2717.6216–2723.64
9–1122.3714–216.5023–215.30

Lines to be planned.

TABLE 3

Start nodeFinal nodeLine impedance (Ω/km)Start nodeFinal nodeLine impedance (Ω/km)
120.192 + 20.533 j4290.515 + 0 j
130.192 + 0 j4300.13 + 0 j
170.192 + 13.200 j5220.13 + 6.482 j
180.61 + 1.820 j7110.299 + 0.823 j
1204.681 + 8.314 j9110.12 + 0.471 j
1327.686 + 14.750 j14210.186 + 0.261 j
261.084 + 3.236 j15160.21 + 0.730 j
292.132 + 3.786 j16270.9 + 0.804 j
2101.703 + 4.690 j16310.731 + 1.298 j
2120.13 + 10.323 j23213.773 + 5.297 j
2130.238 + 16.259 j24291.24 + 2.662 j
2180.238 + 0 j28292.388 + 4.242 j
345.628 + 16.797 j28300.515 + 29.563 j
355.616 + 9.975 j30272.495 + 4.432 j
3140.238 + 10.209 j33294.231 + 3.621 j
3150.476 + 0 j34303.681 + 9.235 j
3174.456 + 7.915 j3537.214 + 3.751 j
3192.7 + 3.790 j4284.396 + 6.172 j
3240.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 algorithm0.69530.86390.68660.7233
Traditional genetic algorithms0.61320.80710.61200.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

CombinationNumber of clusters
10.10.40.40.10.632 10.896 40.711 70.701 67
20.20.30.30.20.677 90.889 10.679 30.708 36
30.30.20.20.30.742 40.833 00.657 10.739 36
40.40.10.10.40.789 10.798 30.625 30.756 15

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

CaseLine connection resultsCluster planning results
ConnectivityDisconnection
17–119–11Cluster 1 (7,8,11,20,32); cluster 2 (2,6,9,10,12,13,18)
23–2114–21Cluster 3 (3,14,17,19,24,26,35); cluster 4 (5,21,22,23,25)
30–2716–27Cluster 5 (15,16,31); cluster 6 (4,27,28,29,30,33,34)
27–119–11Cluster 1 (7,8,11,20,32); cluster 2 (2,6,9,10, 12,13,18)
14–2123–21Cluster 3 (3,14,17,19,21,24,26,35); cluster 4 (5, 22,23,25)
30–2716–27Cluster 5 (15,16,31); cluster 6 (4,27,28,29,30,33,34)
37–119–11Cluster 1 (7,8,11,20,32); cluster 2 (2,6,9,10,12,13,18)
14–2123–21Cluster 3 (3,14,17,19,21,24,26,35); cluster 4 (5,22,23,25)
16–2730–27Cluster 5 (15,16,27,31); cluster 6 (4,28,29,30,33,34)

Planning results of different schemes.

TABLE 7

Comparison itemCase 1Case 2Case 3
Annual investment costs of photovoltaics (PVs)587.54 $553.48 $528.22 $
Annual operation costs of PVs414.93 $232.90 $187.04 $
Annual income of PVs240.43 $247.92 $285.45 $
Annual investment costs of lines622.98 $856.43 $914.28 $
Network losses139.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 itemProposed methodCentralized planning methodDifference rate (%)
Annual investment costs of photovoltaics (PVs)569.86 $541.11 $−5.3
Annual operation costs of PVs293.10 $239.87 $−14.1
Annual income of PVs240.43 $277.97 $13.5
Annual investment costs of lines692.89 $932.30 $25.7
Network losses150.14 $177.58 $15.5
Annual investment costs of energy storages (ESs)1,360.22 $1,689.09 $10.3
Calculation time20.41 s43.62 s53.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), 65836592. 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), 205217. 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), 21872201. 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), 108738108738.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), 7482. 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), 109156109156.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), 821836. 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), 561576. 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), 380394. 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), 18201835. 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), 5564. 10.7500/AEPS20220501001

  • 12

    KoutsoukisN.GeorgilakisP.HatziargyriouN. (2018). Multistage coordinated planning of active distribution networks. IEEE Trans. Power Syst.33 (1), 3244. 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), 109145109145.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), 7683. 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), 109370109370.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), 5158. 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), 136146. 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), 32893303. 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), 3241. 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), 146155. 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), 121128. 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), 11351147. 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), 168176. 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), 10231032. 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), 106455106455.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), 14251433. 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), 3044. 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), 2934. 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), 109103109103.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), 108585108585.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), 4051. 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), 8392. 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), 147155. 10.7500/AEPS20170101002

  • 34

    XuL.YangJ.XiongZ. (2021a). A method of power distribution network classification based on improved grey cluster. Electr. Eng. (13), 8184. 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 24131. 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), 18231829. 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), 109043109043.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), 3544. 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), 24192434. 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), 111119. 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

Reviewed by

Wenliang Yin, The University of Sydney, Australia

Tianbai Deng, Anhui University, China

Updates

Copyright

*Correspondence: Song Yang,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics