Analysis of Distributed Generation Accommodation in Flexible Distribution Networks

This work proposes a method to analyze distributed generation (DG) accommodation in a flexible distribution network (FDN). Firstly, the DG-load matching degree is proposed to quantitatively describe the power balance degree of DG and load in a distribution network. Secondly, the accommodation ratio of DG is proposed and divided into a DG-load accommodation ratio and DG-network-load accommodation ratio, to distinguish whether the index takes the network operational constraints into account. We derive the DG-load accommodation ratio directly from the matching degree and propose the simulation model of sequential production to solve the DG-network-load accommodation ratio. Finally, the cases of FDN under the scenarios of different matching degrees are studied and compared with those of the traditional rigid distribution network. The results show that the improvement of the accommodation ratio by upgrading the rigid distribution network to an FDN is conditional, which is related to not only the matching degree of the whole network but also that of each local network. The DG-network-load accommodation ratio will tend to the DG-load accommodation ratio if proper planning or optimization measures are taken. We find that the capacity of branches adjacent to the DG bus mainly limit the DG accommodation in the FDN, and it is recommended to relocate the DGs and enlarge the capacity of those branches.


INTRODUCTION
The installed capacity of distributed generations (DGs) is rapidly increasing in distribution networks in recent years. Due to the intermittency and uncertainty of DGs, it is difficult to accommodate a large number of DGs in the traditional distribution networks (Stram, 2016).
The development of power electronics technology provides a new idea to solve this problem (Huang et al., 2011;Rueda-Medina and Padilha-Feltrin, 2013). The flexible power electronic equipment, such as the soft open point (SOP) (Zhu et al., 2018), is with great power flow control capability. SOPs refer to the flexible switches in distribution networks usually using back-to-back voltage source converters. SOPs are installed previously at normally open points (NOPs) of the distribution network (Cao et al., 2016). An SOP is able to provide the dynamic and continuous active/reactive power flow controllability and limit the short current (Escalera et al., 2020). The application of the SOP will promote the flexibility and controllability of the distribution system. The SOP improves the capability of the distribution network to host high-penetration DGs (Shafik et al., 2019). The SOP relives the negative effect of intermittent DGs on the operation of distribution networks (Guo et al., 2020). In addition, through controlling the SOPs, the network losses can be reduced and the voltage levels can be well maintained (Flottemesch and Rother, 2004) without involving open/close switch operations. The SOPs can also avoid the power outage when the system is in fault and improve the system reliability (Barragan et al., 2012). Some pilot projects have demonstrated that the SOP is beneficial to a distribution network in many aspects, such as optimal power flow, load balancing, voltage regulation, power supply restoration, and accommodation of DG, which indicates its potential in the future distribution network (Western Power Distribution, 2016;SP Energy Networks, 2016).
The article (Xiao et al., 2018b) proposes the system-level concept of the flexible distribution network (FDN). An FDN is defined as the distribution network with flexible power flow controllability using multiterminal SOPs (including twoterminal SOPs). The article (Xiao et al., 2017a) proposes the total supply capability (TSC) of the FDN. The articles (Xiao et al., 2017b) and (Xiao et al., 2020) adopt the region method to observe the FDN and propose the model and observation approach of the dispatchable region of the FDN. However, the DG accommodation in the FDN has not been studied yet in the existing researches. This work analyzes the DG accommodation in the FDN, which will provide fundamental results for the planning scheme of the FDN to accommodate more DGs.
As for the assessment of the DG accommodation capability of distribution networks, the classical indices include the maximum admissible wind/solar power, wind/solar curtailment rate (Liu and Chu 2018;Ye et al., 2019), etc. However, these indices cannot reflect the essence of DG accommodation capability, which is determined by the location and power balance of DG and load, as well as network operational constraints. This paper proposes the concept of DG-load matching degree, to quantitatively describe the relative value of the power energy output from DG and load consumption in a distribution network. The DG accommodation ratio is proposed and divided into two subindices, including the DG-load accommodation ratio and DG-network-load accommodation ratio. DG-load accommodation ratio is derived directly from matching degree and DG-network-load accommodation ratio is solved by simulation model of sequential production. The proposed DG accommodation analysis method is simple and effective. It directly reflects the effect of DG accommodation capability caused by location and power balance between DG and load, as well as network operational constraints. The accommodation ratio of the traditional rigid distribution network and FDN are calculated and analyzed on case studies. The results discover the bottleneck of DG accommodation in the FDN, and the measures are proposed to remove these bottlenecks.

CONCEPTS OF MATCHING DEGREE AND ACCOMMODATION RATIO
The primary constraint of DG accommodation capability is to keep the real-time power balance of production and consumption. For an ideal distribution system without network operational constraints (i.e., thermal capacity constraints and voltage profile constraint), the DG accommodation capability is determined only by the load power. Therefore, the core idea of DG accommodation is to utilize load. The concept of DG-load matching degree is proposed in this section. Furthermore, the DG-load accommodation ratio and DG-network-load accommodation ratio are proposed to assess the influence of load on the DG accommodation.

Matching Degree
It is said the DG and load are "matched" in a distribution network during a period T only if the following condition is satisfied: (1) Otherwise, it is said the DG and load are "unmatched" if In Eqs. 1, 2, T is the period of observation, and the unit is hour. P ge DG (t) is the sum of maximum available output power of all the DGs at the moment t. |P L (t)| is the sum of consumption power of all the load at the moment t. P ge DG (t) is usually influenced by the weather conditions, such as light intensity and wind speed, as well as by the installed capacity. Equations 1, 2 show that if the DG and load are unmatched, definitely, there is a moment t, the output power of DGs has to be limited.
If the DG and load are matched in a distribution network, that is, Eq. 1 is satisfied, the matching degree, represented by E Ω , is defined as where, T 0 P ge DG (t) dt means the sum of maximum available output power energy of all the DGs during period T.
The DG and load will be better matched if E Ω is larger.
If the DG and load are unmatched in a distribution network, that is, Eq. 1 is not satisfied, the matching degree is defined as , |P L (t)|)dt means the sum of the actual power energy curtailment of all the DGs during period T. According to Eqs. 1-4, it can be seen that "E Ω ≥ 0" is equivalent to "matched," while "E Ω < 0" is equivalent to "unmatched".
The DG output power is characterized by intermittency and uncertainty due to the uncertain weather conditions. Therefore, it is necessary to take uncertain weather conditions into account in the DG-load matching degree, which makes it to be calculated by an integral expression. Generally, the calculation speed can be improved by simplifying the integral. For example, the DG output curve can be simplified by assuming the DG output power is constant during a certain time. The result is relatively conservative, but the calculation speed is improved. In practical application, the speed and accuracy of calculation can be reasonably balanced according to the requirement of planning issues.
In distribution network planning, the DG-load matching degree can be applied to the primary capacity planning for the DG integrated in distribution networks. On the basis of the typical load curve of the distribution network, DG installed capacity can be obtained with the goal of maximum DG-load matching degree. The result is the maximum DG installed capacity such that the output power of the DG can be totally accommodated. When DGs are connected to the distribution network, the DG locations can be optimized by taking them as variables of the DG planning model.

Accommodation Ratio
The accommodation ratio is used in this work to assess the DG accommodation capability of a distribution network. The accommodation ratio, represented by λ DG , is defined as the proportion of the actual output power energy to the maximum available output power energy of all the DGs in a distribution network during period T, as is shown in the following equation: where W ge DG is the sum of the maximum available output power energy of all the DGs during period T. W DG is the sum of the actual output power energy of all the DGs during period T, which is also called accommodation power energy of the DG. |P DG (t)| is the actual output power of all the DGs at the moment t. λ DG ∈[0,1]. The DG accommodation capability of a distribution network will be better if λ DG is larger.
Load power and network operational constraints are the key influences of DG accommodation capability. To analyze the influences separately, the accommodation ratio is divided into two subindices, including the DG-load accommodation ratio and DG-network-load accommodation ratio.
The DG-load accommodation ratio, represented by λ GL , is equal to the accommodation ratio without any network operational constraints. The definition of the DG-load accommodation ratio is consistent with the matching degree, which reflects the overall DG-load power balance. The concept of DG-load accommodation ratio is similar to the "substation capacity-load ratio" (Xiao et al., 2018a), which is a classical index for distribution network planning, and also neglects the network operational constraints.
The DG-network-load accommodation ratio, represented by λ GNL , is equal to the accommodation ratio considering network operational constraints. The definition of the DGnetwork-load accommodation ratio reflects the effects of network operational constraints on the DG accommodation capability, which is effective supplementary to the DG-load accommodation ratio.
The DG-load accommodation ratio (λ GL ) and DG-networkload accommodation ratio (λ GNL ) are depicted in Figure 1. The area between the blue dot line (λ DG 100%) and red curve (λ GL ) is the power energy curtailment of the DG due to the load power limit. The area between the red curve (λ GL ) and green curve (λ GNL ) is the power energy curtailment of the DG due to the network operational constraints limit.
The DG-network-load accommodation ratio curve will tend to the DG-load accommodation ratio curve if proper planning or optimization measures are taken, such as regulating voltage and expanding feeder capacity. Particularly, λ GNL λ GL if the network operational constraints are completely eliminated.
The proposed accommodation ratio is similar to the existing DG penetration (Anderson et al., 2009). The DG penetration can be divided into power penetration, capacity penetration, and energy penetration. The power penetration and capacity penetration are the power ratios of DG to load (Anderson et al., 2009); the energy penetration (Moghaddam et al., 2018), which is most similar to the concept of accommodation ratio, is the energy ratio of DG to load. The accommodation ratio and energy penetration both utilize the power energy of the DG and load to analyze the DG accommodation capability. The difference is that the accommodation ratio, which reflects the utilization of DG, is the proportion of the actual output power energy to the maximum available output power energy of the DG, while the energy penetration, which reflects the supply capability of the DG in the distribution network, is the proportion of the maximum available output power energy of the DG to load consumption. Compared with energy penetration, the accommodation ratio can separately reflect the effects of load power and network operational constraints on the DG accommodation capability directly.

CALCULATION OF THE DG-LOAD ACCOMMODATION RATIO
Since the network operational constraints are not considered in the DG-load accommodation ratio, it can be directly deduced from the matching degree. Firstly, a distribution network is usually not fully connective. However, it can be divided into several fully connective local networks. The division boundary is the normally open tie switches. As shown in Figure 2, the distribution network is divided into four local networks Ω 1 , . . ., Ω 4 .
For a fully connective local network Ω, based on Eqs. 3-5, the DG-load accommodation ratio (λ GL ) can be formulated as where C is a constant. The detailed deduction is shown as follows: (1) When E Ω ≥ 0, that is, the DG and load are matched, (2) When E Ω < 0, that is, the DG and load are unmatched, where C is a constant, It can be seen from Eq. 6 that (1) when E Ω ≥ 0, that is, the DG and load are matched, the output power of the DG can be totally accommodated, and the DGload accommodation ratio is 100%.
(2) when E Ω < 0, that is, the DG and load are unmatched, the DG-load accommodation ratio is determined by the matching degree and sum of the maximum available output power energy of DGs.
The FDN is fully connective because all the mechanical normally open tie switches are replaced with SOPs. Therefore, the DG-load accommodation ratio of the FDN can be directly calculated with Eq. 6.
However, for a rigid distribution network, which is not fully connective, the DG-load accommodation ratio should be calculated as follows. Assuming that a rigid distribution network is composed of N fully connective local networks Ω 1 ,. . ., Ω i ,. . ., Ω N , (1) the DG-load accommodation ratio of each local network is calculated according to Eq. 6, and λ GL,i , i 1, . . . , N is obtained.
(2) λ GL,i , i 1, . . . , N is applied to Eq. 5, and W DG,i , i 1, . . . , N is obtained, where W DG,i is the sum of the actual output power energy of DGs in local network i during period T. Then, W DG of the rigid distribution network can be calculated as (3) W DG is applied to Eq. 5, and the DG-load accommodation ratio of the rigid distribution network is obtained. (1) DC power flow (Purchala et al., 2005) is used. First, the network loss can be included in the power flow of the feeder outlets (Xiao et al., 2011) because the feeders are usually short in length and the network loss ratio is small in the urban power grid. Second, the voltage constraints can be neglected because the system and the DGs are all capable of regulating the bus voltage; thus, the voltage can be kept within the security limits (Shi et al., 2016).
(2) The power cannot flow reversely from 10 kV feeders to an upper-level substation because the upper-level power grid needs great changes for reverse power flow, such as relay protection settings, which is hardly realized in the short term (Fernandez et al., 2020).

Simulation Model of Sequential Production of the FDN
(1) Objective function The objective is to maximize the sum of the actual output power energy of all the DGs during period T, where P DG,k is the actual output power of the DG installed on node k at the moment t. G is the set of all the DG nodes.
(2) Network operational constraints The network operational constraints include power flow constraints and security constraints. The power which flows out of bus, such as load consumption power, is noted as positive, while the power which injects into of bus, such as DG output power, is noted as negative. The power which flows from a substation to a feeder terminal is noted as positive. The power which flows from a feeder to the connected SOP is noted as positive. Due to the assumptions in Assumptions, the constraints are simplified as follows.
The branch flow P Bi,j can be expressed as the algebraic sum of the net power of downstream nodes and the power injected into the downstream SOPs. The power flow equations are formulated as follows: Fi ∈ Ω(SOPm) where B i,j is the branch j of feeder F i . P Bi,j is the power flow of B i,j . ± indicates the direction of power. Ω(B i,j ) represents the set of downstream nodes and SOPs of B i,j . L and G represent the sets of all the load and DG nodes, respectively. |P k | is the net power of node k. P L,k is the load power of node k. ΔP i,m is the power flow between F i and SOP m . Ω(SOP m ) represents the set of feeders connected to SOP m . Equation 11 is the power flow calculation of the branch. Equation 12 is the node power equation. Equation 13 is the equilibrium of active power from each terminal of the SOP. Due to the assumptions in Assumptions, the voltage constraints are neglected. The security constraints of the FDN with DG installed are mainly thermal capacity constraints, including the feeder capacity constraints, node power constraints, DG output constraints, SOP capacity constraints, and the reverse power flow constraints.
where C B i,j represents the capacity of B i,j . B is the set of all the branches. P max k is the maximum permitted power of node k. P ge DG,k is the maximum available output power of the DG installed on node k. C SOP m is the capacity of the SOP m terminal. S is the set of all the SOPs. P

Algorithm
The model in Simulation Model of Sequential Production of the FDN is linear, which can be solved by the linear programming software after identical deformation (Xiao et al., 2011). The Frontiers in Energy Research | www.frontiersin.org April 2021 | Volume 9 | Article 634770 detailed deformation is shown in Supplementary Material 1. In this paper, LINGO is used to solve the model to obtain W DG . W DG is applied to Eq. 5, and the DG-network-load accommodation ratio is obtained.

CASE STUDY
In this section, cases of FDN with different matching degrees are studied. For three scenarios with different DG planning, the DG-load accommodation ratio and DG-network-load during 1 h are calculated, and their trends during one day are analyzed. Meanwhile, the DG-load accommodation ratio and DG-network-load accommodation ratio of a rigid distribution network are also analyzed as a comparison. The bottlenecks of the FDN to improve DG accommodation capability are analyzed, and the improvement measurements are proposed.

Case Grid
The case grid of the FDN is shown in Figure 3. For comparison, the case grid of a rigid distribution network (RDN) with the same topology is used, and the FDN is divided into two local networks corresponding to the RDN. The total daily power energy of the load is 597.

Matching Degree Calculation
The daily curve of the maximum available output power of DGs and actual consumption power of the load is shown in Figure 4.

Daily Trend of the Accommodation Ratio
According to Figure 4, for the case grids of both the FDN and rigid distribution network (RDN), the DG-load accommodation ratio and DG-network-load accommodation ratio during each hour are calculated using Eqs. 6, 9. Then, the DG-load  accommodation ratio (GLAR, represented by λ GL ) and DG-network-load accommodation ratio (GNLAR, represented by λ GNL ) of all 24 h are used to form the daily accommodation ratio curve, as is shown in Figure 5.
It can be seen from Figure 5 that (1) the DG-load accommodation ratio of the FDN is higher than that of the rigid distribution network during 66.7% of a day's period. The reason is that the FDN can adjust the power flow flexibly to increase DG accommodation.
(2) the DG-network-load accommodation ratio of the FDN is higher than that of the rigid distribution network during 58.3% of a day's period. The reason is that the FDN increases the DG accommodation by adjusting the power flow, while the increase extent of DG accommodation is limited by network operational constraints. For the PV accommodation of local network 1, it is focused in midday that the DG-network-load accommodation ratio of the FDN is higher than that of the rigid distribution network. For the WT accommodation of local network 2, it is focused from dusk to the next morning such that the DGnetwork-load accommodation ratio of the FDN is higher than that of the rigid distribution network.
(3) the DG-load accommodation ratio of the FDN can reach 100% anytime. The reason is that the FDN is capable to fully accommodate the DGs by adjusting the power flow flexibly, when the load and DG of the whole network are matched. (4) the DG-network-load accommodation ratio of the FDN is lower than the DG-network-load accommodation ratio during 37.5% of a day's period. The reason is that network operational constraints limit the FDN adjusting power flow.

Analysis of the Daily Accommodation Ratio
According to Figure 4, for the case grids of both the FDN and rigid distribution network (RDN), the DG-load accommodation ratio (GLAR) and DG-network-load accommodation ratio (GNLAR) during a whole day are calculated using Eqs. 6, 9, which is shown in Table 2. It can be seen from Table 2 that (1) since matching degree of scenario 1 is positive, the daily DGload accommodation ratio of the FDN is 100%. (2) due to the network operational constraints, the daily DGnetwork-load accommodation ratio of local network 1 and local network 2 of the FDN case is reduced by 4.6% and 1.2%, respectively. The daily DG-network-load accommodation ratio of the whole network of the FDN is reduced by 2.9%. (3) since matching degrees of local network 1 and local network 2 are all negative, the rigid distribution network case cannot fully accommodate all DGs, and the daily DG-load accommodation ratio of the rigid distribution network is lower than 100%. (4) due to the network operational constraints, the daily DGnetwork-load accommodation ratio of local network 1 and local network 2 of the rigid distribution network case is reduced by 4.8% and 0.7%, respectively. The daily DGnetwork-load accommodation ratio of the whole network of the rigid distribution network is reduced by 2.7%.  (5) compared with the rigid distribution network, the FDN can improve the daily DG-load accommodation ratio by 9.1% and daily DG-network-load accommodation ratio by 8.9%.

Matching Degree Calculation
The daily curve of the maximum available output power of DGs and actual consumption power of the load is shown in Figure 6. According to Figure 6, the matching degree (MD) is calculated using Eqs. 3, 4. For the whole network, MD −16.8%. For local network 1, MD −20.3%. For local network 2, MD −48.2%.

Daily Trend of the Accommodation Ratio
According to Figure 6, for the case grids of both the FDN and rigid distribution network (RDN), the DG-load accommodation ratio and DG-network-load accommodation ratio during each hour are calculated using Eqs. 6, 9. Then, the DG-load accommodation ratio (GLAR, represented by λ GL ) and DG-network-load accommodation ratio (GNLAR, represented by λ GNL ) of all 24 h are used to form the daily accommodation ratio curve, as is shown in Figure 7.
It can be seen from Figure 7 that (1) the DG-load accommodation ratio of the FDN is higher than that of the rigid distribution network during 83.3% of a day's period. The reason is that the FDN can adjust the power flow flexibly to increase DG accommodation.
(2) the DG-network-load accommodation ratio of the FDN is higher than that of the rigid distribution network during 83.3% of a day's period. The reason is that the FDN increases the DG accommodation by adjusting the power flow, while the increase extent of DG accommodation is limited by network operational constraints. For the PV accommodation of local network 1, it is focused that the DG-network-load accommodation ratio of the FDN is just a little higher than that of the rigid distribution network in daytime. For the WT accommodation of local network 2, it is focused from dusk to the next morning that the DGnetwork-load accommodation ratio of the FDN is much higher than that of the rigid distribution network.
(3) the DG-load accommodation ratio of the FDN cannot reach 100% during 70.8% of a day's period. The reason is that the FDN cannot fully accommodate the DGs by adjusting the power flow, when the load and DG of the whole network are unmatched. (4) the DG-network-load accommodation ratio of the FDN is lower than the DG-load accommodation ratio during 37.5% of a day's period. The reason is that network operational constraints limit the FDN adjusting power flow.

Analysis of the Daily Accommodation Ratio
According to Figure 6, for the case grids of both the FDN and rigid distribution network (RDN), the DG-load accommodation ratio (GLAR) and DG-network-load accommodation ratio (GNLAR) during a whole day are calculated using Eqs. 6, 9, which is shown in Table 3. It can be seen from Table 3 that (1) since the matching degree of scenario 2 is negative, the daily DG-load accommodation ratio of the FDN is 81.9% (<100%). (2) due to the network operational constraints, the daily DGnetwork-load accommodation ratio of local network 1 and local network 2 of the FDN case is reduced by 2.2% and 3.5%, respectively. The daily DG-network-load accommodation ratio of the whole network of the FDN is reduced by 3.0%. (3) since matching degrees of local network 1 and local network 2 are all negative, the rigid distribution network case cannot fully accommodate all DGs, and the daily DG-load accommodation ratio of the rigid distribution network is lower than 100%. (4) due to the network operational constraints, the daily DGnetwork-load accommodation ratio of local network 1 and local network 2 of the rigid distribution network case is reduced by 1.3% and 0.4%, respectively. The daily DGnetwork-load accommodation ratio of the whole network of the rigid distribution network is reduced by 0.8%. (5) compared with the rigid distribution network, the FDN can improve the daily DG-load accommodation ratio by 16.3% and daily DG-network-load accommodation ratio by 14.1%.

Matching Degree Calculation
The daily curve of the maximum available output power of DGs and actual consumption power of the load is shown in Figure 8.

Daily Trend of the Accommodation Ratio
According to Figure 8, for the case grids of both the FDN and rigid distribution network (RDN), the DG-load accommodation ratio and DG-network-load accommodation ratio during each hour are calculated using Eqs. 6, 9. Then, the DG-load accommodation ratio (GLAR, represented by λ GL ) and DGnetwork-load accommodation ratio (GNLAR, represented by λ GNL ) of all 24 h are used to form the daily accommodation ratio curve, as is shown in Figure 9.
It can be seen from Figure 9 that (1) the DG-load accommodation ratio of the FDN is equal to that of the rigid distribution network anytime. The reason is that the rigid distribution network is capable of accommodating the DGs; thus, the FDN cannot increase DG accommodation by adjusting the power flow.
(2) the DG-network-load accommodation ratio of the FDN is higher than that of the rigid distribution network during 16.7% of a day's period. The reason is that DG accommodation in the rigid distribution network is limited by network operational constraints, and it is increased in the FDN by adjusting the power flow. For the PV accommodation of local network 1, it is focused in midday that the DG-network-load accommodation ratio of the FDN is higher than that of the rigid distribution network. For the WT accommodation of local network 2, it is focused from dusk to the next morning that the DGnetwork-load accommodation ratio of the FDN is higher than that of the rigid distribution network. (3) the DG-load accommodation ratio of the FDN can reach 100% anytime. The reason is that the FDN is capable of fully accommodating the DGs, when the load and DG of the whole network are matched. (4) the DG-network-load accommodation ratio of the FDN is lower than the DG-load accommodation ratio during 20.8% of a day's period. The reason is that network operational constraints limit the FDN adjusting power flow.

Analysis of the Daily Accommodation Ratio
According to Figure 8, for the case grids of both the FDN and rigid distribution network (RDN), the DG-load accommodation ratio (GLAR) and DG-network-load accommodation ratio (GNLAR) during a whole day are calculated using Eqs. 6, 9, which is shown in Table 4.
It can be seen from Table 4 that (1) since the matching degree of scenario 3 is positive, the daily DG-load accommodation ratio of the FDN is 100%. (2) due to the network operational constraints, the daily DGnetwork-load accommodation ratio of local network 1 and local network 2 of the FDN case is reduced by 2.0% and 1.0%, respectively. The daily DG-network-load accommodation ratio of the whole network of the FDN is reduced by 1.4%.
(3) since matching degrees of local network 1 and local network 2 are all positive, the daily DG-load accommodation ratio of the rigid distribution network is 100%. (4) due to the network operational constraints, the daily DGnetwork-load accommodation ratio of local network 1 and local network 2 of the rigid distribution network case is reduced by 7.2% and 3.2%, respectively. The daily DGnetwork-load accommodation ratio of the whole network of the rigid distribution network is reduced by 4.8%. (5) compared with the rigid distribution network, the FDN cannot improve the daily DG-load accommodation ratio under scenario 3.

Overall Analysis for 3 Scenarios
In Figure 10, the results of the FDN and rigid distribution network (RDN) in three scenarios are summarized, including DG-load matching degrees (MDs, represented by E Ω ), DG-load accommodation ratios (GLARs, represented by λ GL ), and DGnetwork-load accommodation ratios (GNLARs, represented by λ GNL ).
According to the conclusion of Scenarios 1-3 and Figure 10, it can be seen that   (1) the DG-load accommodation ratio of a distribution network is determined by the matching degree of each local network. Particularly, since the FDN is overall connective, the DGload accommodation ratio of the FDN is only determined by the overall matching degree.
(2) Only when at least one local network of a distribution network is DG-load unmatched can the DG-load accommodation ratio of the FDN be larger than that of the rigid distribution network. The increase extent of the DGload accommodation ratio is determined by the matching degrees of the whole network and local networks together. Particularly, if the load and DG of every local network are matched, that is, the DG-load accommodation ratio of every local network reaches its maximum value 100%, it is evident that the DG-load accommodation ratio cannot be further improved by any measures, including the FDN upgrading. In this case, the rigid distribution network is capable of accommodating the DGs due to the perfect DG-load coordinated planning, although it is unusual in a real rigid distribution network. (3) the increase extents of the hourly DG-load accommodation ratio and hourly DG-network-load accommodation ratio by the FDN fluctuate during one day. It can be seen that the improvements of the hourly DG-load accommodation ratio and hourly DG-network-load accommodation ratio are more obvious in midday and night because the output power of PVs and WTs is sequentially complementary and PVs and WTs are just located in different local networks in the cases. (4) the DG-network-load accommodation ratio is always lower than the DG-load accommodation ratio in a real FDN. This is because the network operational constraints, such as thermal capacity constraints, limit the DG accommodation.
According to the analysis of all three scenarios, we find that it is the branches adjacent to the DG bus that mainly limit the DG accommodation.

Measures to Improve the DG-Network-Load Accommodation Ratio of the FDN
To remove the bottlenecks of DG accommodation in the FDN, the following measures are proposed: (1) If the locating of DGs is under discussion, it is recommended that the DGs should be integrated to the areas with enough load to ensure the output power of DGs can be accepted instead of exporting.
(2) If the locating of DGs is finished or measure 1 is not enough to remove the bottlenecks of DG accommodation, it is recommended that the capacity of branches adjacent to the DG bus should be properly enlarged according to the maximum output power of DGs. Besides, the capacity of SOP terminals should not be smaller than that of adjacent feeders.
For the cases in this work, the specific measures are as follows: Measure 1: change the location of DGs, enlarge the capacity of feeders, and increase the capacity of the SOP terminal to 10.5 MVA, as is shown in Table 5.
Measure 2: based on measure 1, adjust the incremental capacity. The capacity of branch B 1,5 is increase by 0.9 MVA, B 1,7 0.4 MVA, B 1,8 4.8 MVA, B 2,4 0.7 MVA, and B 2,7 8.8 MVA. The capacity of the rest of the branches is increased as shown in Table 5. The details are shown in Supplementary Table S3.
The cases with overall higher capacity of DGs are also studied to verify the conclusion mentioned above, as is shown in Supplementary Material 4.

CONCLUSIONS
The FDN is a looped-operational distribution network using advanced power electronics technology, such as the SOP. This article studies the intermittent DGs accommodation in the FDN (1) The matching degree is proposed to describe the relation of the load and location and size of the DG.
(2) The proposed index, accommodation ratio, is further divided into the DG-network-load accommodation ratio and DG-load accommodation ratio, to distinguish whether the index takes the network operational constraints into account. (3) The calculation of the DG-load accommodation ratio is proposed based on the matching degree; the calculation of the DG-network-load accommodation ratio is proposed based on the simulation model of sequential production. (4) The accommodation ratio of the FDN and rigid distribution network is studied on cases with different matching degrees, and it is concluded that: (a) the DG-load accommodation ratio is determined by the matching degree only. The DG-load accommodation ratio of the FDN is determined by the overall matching degree of the whole network. The DG-load accommodation ratio of the rigid distribution network is determined by the overall matching degree of the whole network, as well as the matching degrees of local networks. (b) when at least one local network of a distribution network is DG-load unmatched, the DG-load accommodation ratio of the FDN will be larger than that of the rigid distribution network because the SOP can redistribute the power flow among different local networks. (c) the DG-network-load accommodation ratio is determined by the matching degree and network operational constraints. The DG-network-load accommodation ratio will tend to the DG-load accommodation ratio if proper planning or optimization measures are taken, such as regulating voltage and expanding feeder capacity. (d) the branches adjacent to the DG bus mainly limit the DG accommodation in the FDN. It is recommended to relocate the DGs and enlarge the capacity of feeders.
Compared with traditional distribution network analysis, the proposed DG accommodation analysis method can reflect the essence of DG accommodation capability clearly, which is determined by the location and power balance of the DG and load, as well as network operational constraints. Under the guidance of DG-network-load coordination and planning, the advantages of the FDN can be maximized. Future research will consider the reverse power flow in the 35 kV and 110 kV network and the influences of voltage constraints and energy storage.

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.