1. Introduction
In the last decades, the installed capacity of distributed energy resources (DER) increased rapidly [
1]. Distribution networks, initially planned to supply customers, are now connecting a large number of DER which leads to voltage rise and loading constraints [
2]. In [
3], PV was stated to be in the two most installed sources of electricity in the European Union in 2013.
In [
4], it is reported that solar PV is already covering more than 7% of the electricity demand in three countries (Italy, Greece and Germany). Further, in 2014, the power generation capacities for Wind and solar PV in the EU 28 increased by 11,791 MW and 6574 MW, respectively. Globally, about 40 GW solar PV have been installed in 2014. Moreover, the total amount of installed Solar Power in the world reached 178 GW and covers more than 1% of the world electricity demand.
However, to reduce CO
2 emission even more renewable energy resources are required [
5,
6,
7,
8]. Acknowledging the challenges for distribution grids, efforts are needed to limit investments in network reinforcements required to integrate the expected amount of DER.
The detailed investigation of the hosting capacity of the entire supply of a DSO may be a computationally burden. In [
9] the LV-grid data of the largest Dutch DSO is classified. This DSO operates about 88,000 LV-feeders and authors report, that the individual analysis of all LV-feeders is enormously time consuming. In [
10] for example, only 1.35% of the grid data of a DSO was analyzed and three main feeder topologies were presented. This issue can be solved in several ways. First, a common approach in the literature is to cluster networks and focus on representative networks of each cluster. Results are then scaled up for the entire supply area. Several reference networks were defined in the past. For example, in [
11], representative networks were defined for Bavaria based on real network data and this set of representative networks were updated in [
12]. In [
13] it was shown that clustering of distribution networks by means of electrical/topological indices can be reduced to several relevant parameters and a criterion to define a reasonable number of clusters is given. In [
14], a Matlab GUI is presented, which generates statistically correct distribution grid models. Authors point out that such generators are crucial for new applications development such as monitoring and control.
Third, calculations can be distributed to several machines or cores. In [
15], the issue that network calculation tools mostly run on single core architectures running calculations sequentially is addressed. For this reason a very promising high performance computation framework called GridPACK (Battelle Memorial Institute, Columbus, OH, USA) was developed. For a 156-bus system and 1024 contingencies, the overall solution time could be reduced from 2 h to less than a minute.
The second approach to reduce computation time is to restrict the number of nodes in a power flow calculation. According to [
16], low voltage networks can be reduced to a two-bus equivalent model for pre-screening of large supply areas to identify networks where more likely voltage violations may occur. Moreover, “continuous feeders” are used for studies on specific inverter functionalities to control the voltage [
17]. In [
18], benchmark feeders are proposed using a similar approach.
Generic networks have received attention for benchmarking in the frame of specific DSO-revenue regulations. Thereby, tools and methods known from graph theory such as the minimum spanning tree may be suitable for the benchmark of distribution grids. In [
19], the Prize-Collecting Steiner Tree algorithm is utilized to connect customers (e.g., optic fiber, district heating sector) cost-efficiently to the network. Moreover, in [
20,
21], the Facility Location Problem is applied to LV-networks to assist network planners and operators to optimally place secondary substations. In [
22], the “betweenness” is used as a measure to identify vulnerable lines according to either topological point of view or by considering the actual power flows.
In [
23], Kendall’s correlation was used to reduce redundancy between input parameters. Among 13 input parameters, the total cable length of feeders, the share of cables and aerial lines, connection points per length, the yearly consumption and the installed capacity of PV-systems were identified as most important parameters. In [
24], distribution grids of western Australia are classified and prototypical feeders are selected at both MV and LV-level. On LV-level, 26 parameters were investigated, such as secondary transformer rating, feeder length, energy consumption (for residential, commercial and industrial respectively), number of customers and many more. The identified most significant seven parameters are similar to previously mentioned studies: underground cable length, load capacity of residential customers, number of residential customers, overhead-line length and others. The cluster analysis on LV feeder level resulted in an optimal cluster number of eight and a description of the LV feeder subsets is given.
This paper presents findings from an environment that is capable of analyzing a high number of networks and performing various in-depth studies. During the project IGREENGrid ((IntegratinG Renewables in the EuropEaN Electricity Grid) [
25], the complete distribution grid data of two Austrian DSOs with more than 14,000 LV-networks were provided from their GIS systems and imported into PowerFactory (DIgSILENT, Gomaringen, Baden-Württemberg, Germany) [
26]. Thereby, the provided data by the DSOs had a different format. In the case of DSO A, the network data stored in a GIS database were exported in the DIgSILENT DGS format. In the case of DSO B, the data were provided from the GIS database as Neplan export (*.cde, *.edt, and *.ndt files). In
Figure 1, the import process, which was automated with the DIgSILENT Programming Language (DPL) is depicted. In total about 14,600 networks with approximately 54,100 feeders were imported. For about 41,000 feeders, load flow and short circuit calculations were performed. After a plausibility check, more than 36,000 feeders were assessed.
A high number of feeder parameters such as the number of load, number of nodes, average distance to neighboring nodes, average number of nodes, equivalent sum impedance, ampacity of lines, feeder length, etc. were investigated. The hosting capability (HC) of LV-feeders was investigated considering reactive power control strategies and several DER-scenarios for a high number of voltage limits. Based on these results, feeders can be classified according to the limiting constraint (voltage or loading). Further, a sensitivity analysis of the limiting constraint on the admissible voltage rise is performed. Additionally, parameters of feeders are investigated and probability tests for selected parameters are performed. This analysis is supplemented by a linear (Pearson) and a rank (Spearman) correlation analysis of feeder parameters and admissible DER penetration levels according to the considered scenarios.
Based on these results, the impact of the DER-scenario on the constraint reason is evaluated. Moreover, suggestions to either increase the admissible voltage band or apply a reactive power control strategy to integrate a higher number of DER are given. Furthermore, for the grid data of one investigated DSO whether the minimal spanning tree is a suitable tool to benchmark DSOs is analyzed.
The paper is organized as follows: First the Materials and Methods are described in
Section 2. Results on feeder classification and high-level indicators are presented in
Section 3. In
Section 4 a comprehensive study on the DER-scenarios and reactive power control strategies is given. In
Section 5, the correlation between feeder parameters and the HC is investigated. The paper is concluded in
Section 6.
2. Materials and Methods
In this section, the hosting capability (HC) is introduced, followed by scenario definitions regarding the DER-distribution, reactive power control strategy and admissible voltage rise. After that, the two most important constraints limiting the HC are discussed. Next, the assessed feeder parameters and the algorithms to perform the simulation are presented. Finally, two benchmarking measures are presented.
In this work, the amount of DER that can be hosted in LV-feeders is investigated. For that, the hosting capability is defined:
“The hosting capability of a LV-feeder is the amount of DER that can be installed for a given DER-scenario while meeting voltage and loading prerequisites in the feeder.” The hosting capacity evaluation considers as system boundary the LV feeders (not considering the distribution transformer). Moreover, loads are not explicitly taken into account (see the justification below) and further factors that might limit the hosting capacity such as the impact on protection schemes are not considered here (considered as secondary factors).
Load assumptions have in general a high impact on the achievable HC. In [
27], it is reported that “The worst cases of voltage rise occur when the PV-DG is at or close to its maximum output, and the difference between the load and the generation is at its highest” and further “Using deterministic models with average values for the loads underestimates the voltage impact of increasing PV-DG because the real loads on a feeder vary”. Generally, the base load can be well monitored at secondary transformers supplying a particular area with a high number of customers. These customers are usually supplied by more than one feeder. However, given the small number of loads and generators connected to a single LV feeder, the existence of a certain base load level may not be given permanently at particular nodes in the feeder. In [
28], a stochastic analysis was performed on the smart meter data of 1077 customers (households and PV generators) measured over a period of one year. Authors evaluated the peak and average load depending on the number of households in the respected area that needs to be considered in network planning. For a feeder with five households a significantly higher peak load (15 kW) and lower base load (0.5 kW) compared to 100 households (200 kW, 22 kW respectively) has to be considered. These values are measured at the beginning of the feeder. Therefore, reasonable base load assumptions could not be defined at specific nodes (e.g., end node). The simultaneity factor for the balanced PV infeed as well as the orientation and tilt is assumed to be identical for all generators (worst case). Under these conditions, the hosting capability of a feeder heavily depends on the location of the generation along the feeder as well as on the allowed voltage rise and loading limit of lines.
A high penetration level of DER in distribution grids may lead to protection issues that need to be avoided. Authors in [
29] conclude that the fault devices need to handle the additional fault current by DER and that for fuse-fuse coordination a higher fault current by DER is beneficial, while for the fuse-recloser coordination DER the fault injection should be minimal. In [
30], it is stated that with the increasing DER penetration, the radial structure of distribution feeders will be changed to a meshed configuration. Further, a methodology is proposed to determine the size of DG sources while the existing protection scheme can be maintained. Authors in [
31] point out that a large-scale penetration of DER may lead to protection blinding and sympathetic tripping and other failures in meshed MV-grids. However, to which extent these issues can be translated to radial LV-feeders, requires further investigations.
According to the hosting capability definition, the remaining constraint that could be violated is that the total hosting capability of all feeders may exceed the rating of the secondary transformer supplying these feeders. However, whether this issue arises depends on several circumstances:
Could the full hosting capability potential be utilized (e.g., available DER roof-top potential)?
Which time frame can be expected for the full deployment of DER?
What is the actual minimal load at the transformer and expected growth of consumption?
Are network topology changes foreseen (additional feeders, feeders supplied by other secondary substations)?
Hence, there are remaining reserves regarding the full deployment of the hosting capability of all LV-feeders supplied by a secondary substation.
Three DER scenarios were considered (
Figure 2 left):
“Uniform”: Generation resources placed uniformly along the feeder (at all connection boxes) and then scaled-up to reach the hosting capability (loading limit or voltage limit is reached).
“Weighted”: Generation resources distributed according to the summed annual consumption at the connection boxes) and then scaled-up to reach the hosting capability. This scenario is relevant when assuming that in the future most of DER deployment will be close to the loads and driven by self-consumption.
“End of feeder” (eof): Generation resources connected at the “end node” which is defined as the node with the highest voltage for the uniform scenario (this scenario can be considered as “worst-case”).
Beside the DER scenarios, different reactive power control concepts have been considered (
Figure 2 middle): a fixed PF = 0.9 and Q(U)-control. For the Q(U)-control, a minimal PF of 0.9 is allowed. Furthermore, voltage regulated distribution transformer VRDT (transformer with On-Load-Tap-Changer (OLTC)) are considered in this study since they lead to an increase of the hosting capacity by decoupling medium and low voltage networks. This decoupling allows using a larger share of the voltage band for the LV network than the share currently allocated (without VRDT). Hence, a voltage rise of up to 8% of the nominal voltage (
Figure 2 right) may be admissible.
The considered Q(U) control is depicted in
Figure 3 and consists of a dead band, droop and a saturation area for over and under voltage. Thereby the set points of the Q(U)-control are configured according to the admissible voltage rise. For the dead band, 20% of the voltage band is allocated, while for the saturation area 30% is reserved. The remaining voltage band is dedicated to the droop control.
The concept of hosting capacity introduced in [
32] is restricted to the two most relevant limitations in distribution networks: the maximal admissible voltage rise and the maximal admissible loading. An illustration of this is shown in
Figure 4. If the hosting capacity of a voltage constrained feeder is increased by means of reactive power control, the loading reserve is reduced (Feeder 1). However, if the loading reserve is rather small, as in the case of Feeder 2, the full potential of reactive power control strategies cannot be deployed since the feeder would be overloaded by the additional reactive power flow.
The problem formulation for the optimization of the HC of a feeder is given in Equation (1), where
f is a feeder of Network
N,
ui the voltage at bus
i,
lj the loading of line
j of the feeder. The admissible loading
Lmax of lines is defined as 100% of the respective nominal current. The voltage limit
umax is determined by the investigated scenario (1.01 to 1.08 p.u.).
Based on the assumption in
Figure 2 and the problem formulation in Equation (1), the HC is obtained with an adapted Newton–Raphson (NR) algorithm to reach one of the two constraints depicted in
Figure 4 (voltage or loading). For that, two NR-optimizations are performed to calculate the highest HC without violating the admissible voltage rise and the loading of network elements (
Figure 5). In each step, the lower estimated HC is selected for the next iteration.
By means of active power curtailment (P(U)-control), the HC could be increased even further. However, the objective is to maximize the HC of feeders according to the DER-distribution. A P(U)-control would curtail the injected power unequally, since the power would in tendency be more curtailed at the most distant installation where the highest voltage rise occurs. This would cause varying injected powers for the DER-scenario uniform. Hence, considering a P(U)-control in the optimization of the HC would lead to a contradiction that cannot be resolved while maximizing the injected power for all generators in the feeder in a fair way. The amount of curtailed energy however depends on several other factors such as location, connection point, orientation and tilt of the installation, load behavior, etc.
Parameters used for classification in [
11,
33] are considered as well as several other parameters are investigated, which might be suitable to estimate the HC without power flow calculations.
Table 1 gives an overview on investigated parameters. Parameters used in literature to classify networks are underlined. A key design element of the platform is that the data analysis runs on feeder level. This allows incomplete data leading only to the exclusion of affected feeders, while remaining feeders in the same network are processed. Some parameters are suitable to check the plausibility of the feeder data. Parameters listed to validate the classification are available for validation purposes. However, these parameters are listed by benefits, e.g., a classification to estimate the hosting capability has the highest benefits. The equivalent load location is the node in the feeder, where the total load has to be placed to cause a similar voltage drop compared to the initial situation. The aggregated obtained hosting capability of feeders may exceed the nominal rating of the secondary transformer. Therefore, the transformer data are also assessed in the static analysis.
The formulas for not evident parameters are given in Equations (2)–(7), where N is the number of Nodes in the feeder.
In
Figure 6, the implemented algorithm is depicted. Each network is assessed for given input parameters (Voltage Limits, control strategies and DER distribution). After a plausibility test, feeders are defined and the secondary substation of the networks is replaced by a LV slack node. Omitting the transformer allows to study only feeder relevant parameters and avoid mischaracterization of identical feeders caused by different secondary transformers types. First, for each feeder, parameters that are invariant of any assumptions are calculated (static parameters such as the feeder length,
ZK and
ZΣ). Parameters such as the
ZK or
Inom can be used to check the plausibility of network data (incomplete or corrupt). Next, the HC for each scenario according to
Figure 2 is calculated and parameters based on these assumptions are obtained (scenario results). Thereby,
i defines the control strategy (uncontrolled, cosϕ = 0.9, VoltVAr) and
j the DER-scenario (uniform, eof, and weighted). For each scenario, the NR-algorithm is utilized to obtain the maximal HC under the given voltage (
umax) constraint without overloading any network element.
The algorithm is designed to export one file per network and the combination of voltage limit, control strategy, distribution scenario as well as one file per parameter group (static feeder data, parameters requiring a power flow or short-circuit calculation, etc.). In the current configuration, there are 85 files per network available for analysis (
Figure 7). Hence, this allows running or repeating parts of the study in a modular way. Therefore, two ways of parallelizing the computation to reduce the running time are possible: distributing the grid data or running specific scenarios for the same networks on several machines.
The running time for all presented scenarios in
Figure 6, for 14,600 networks of two DSOs depends on the number of available cores and machines. For each feeder in the networks a NR-optimization with several iterations is performed to calculate the HC according to the scenarios. Running on a single core on a single machine (worst case) e.g., on a standard laptop with an Intel
® Core™ i7-2640M @2.8Ghz processor (Intel, Santa Clara, CA, USA), about 40–45 networks can be processed per hour (networks have on average 4.4 feeders and feeders have on average 22 nodes). For each scenario and feeder, several contingencies are calculated with the presented NR-optimization to obtain the maximal HC.
Some regulatory concepts require the comparison of different DSOs across the country to benchmark the capital expenditure (CAPEX) and operational expenditure (OPEX). For example, the weighted transformed connection density (trfNAD) is used as one component to compare DSO in the benchmarking process [
34] in Austria. Equation (8) shows a simplified version of this parameter, where A is the supplied area and N is the number of connection points in the considered area.
4. Results on the Hosting Capability of Feeders
Reactive power control strategies aiming to reduce the voltage rise introduced by DER infeed are causing an additional loading of lines. Before deploying such smart grid solutions, a sufficient reserve in terms of the line rating should be remaining. In
Figure 16, the share of voltage (blue) and current (red) constrained feeders on the “U/I plane” is depicted for the two investigated DSOs. Each point on the main part of the figures (upper right) represents a feeder which is colored according to the constraint for an admissible voltage rise of 3%.
In the case of DSO A,
Figure 16a shows that about 90% of the feeders are voltage-constrained, confirming that meeting the voltage limits is dominantly limiting the hosting capability for the great majority of feeders Additionally, the plot on the left shows that most of the voltage-constrained feeders (blue) are “far from the corner” (voltage and current constraint) which means that the maximal loading of most of the voltage-constrained feeders is rather far from the 100% loading limit. For 50% of the voltage-constrained feeders, the maximal loading is below 37%, which means that there is a large reserve in respect to overloading lines. Accordingly, there is a priori a large deployment potential of smart grids solutions aiming at controlling the voltage (without observing the loading) without taking the risk or running into overloading network elements.
In the case of DSO B (
Figure 16b), the share of loading constrained feeders is significantly higher compared to DSO A. The share of loading constrained feeders is 21% and voltage constrained feeders 77%. Further, fewer feeders have a high loading reserve compared to DSO A. The loading reserve for about 60% of the feeders is less than 50%. Therefore, the share of feeders where the full potential of reactive power control strategies cannot be reached is higher.
One can expect that the DER scenario and the allowed voltage limit have a high impact on the hosting capability of feeders. Furthermore, the constraint reason (voltage or loading) may vary with these assumptions as well.
Figure 17 demonstrates that the constraint limiting the hosting capability is rather independent for lower range of voltage limits from the investigated DER scenarios. At
umax = 1.01 p.u., almost 100% of the feeders are voltage constrained. For higher voltage limits, the share of voltage constrained feeders is steadily decreasing.
In
Figure 17a, less than 40% of the feeders remain voltage constrained at
umax = 1.08 p.u. for the DER-scenario uniform. Up to an admissible voltage rise of 4%, the constraint reason is rather independent from the considered DER-scenario. At an admissible voltage limit of 8%, the difference of the constraint reason for the DER scenarios uniform and eof is less than 10%.
With an activated VoltVAr-control (according to [
38,
39]) for the DER-scenario uniform, a little more feeders remain voltage constrained compared to a PF = 0.9 (
Figure 17b). An application of this approach is presented in [
35]. With reactive power control the share of loading constrained feeders at
umax = 1.08 p.u. reaches 80%, which is significantly higher compared to
Figure 17a. In conclusion, the DER-scenario eof is suitable to estimate if a feeder will be either voltage or loading constrained for a given admissible voltage rise under the DER-scenario uniform. In particular, for PF = 1 an estimation up to an admissible voltage rise of 4% is accurate. Considering PF = 0.9 (
Figure 17b), for the full range of voltage limits accurate results regarding the constraint reason can be expected.
In
Figure 18, the HC results for two DER-scenarios (uniform and end of feeder) with and without reactive power control strategies (PF = 1, PF = 0.9 and VoltVAr) are depicted.
In
Figure 18a, the distribution of the HC of all feeders for different voltage levels are depicted for the DER-scenario uniform and PF = 1 for both DSOs. Up to 1.07 p.u., a steadily increase of the HC can be observed. At 1.07 p.u. and 1.08 p.u., the 95th percentile remains stable, but the median and the 5th percentile are still increasing by about 10 kW and 5 kW, respectively per voltage step. Concerning the allowed voltage rise, allowing a voltage rise of 1.02 p.u. instead of 1.01 p.u. leads to a significant HC increase. In that case, the median value can be doubled from about 25 kW to more than 50 kW. If the admissible voltage rise can be increased to 4%, more than 50% of the feeders have a HC higher than 100 kW. For higher admissible voltage limits a higher share of feeders become loading constrained. Hence, a voltage increase has no impact on the achievable HC anymore, since a cable in the feeder may be fully loaded. Therefore, the maximal HC is reached.
Next, the HC study results for the DER-scenario uniform and a PF = 0.9 control are depicted in
Figure 18b. In this case, the median is always higher compared to the uncontrolled case. The 95th percentile is lower except for 1.01 p.u., 1.02 p.u. and 1.04 p.u., but the 5th percentile is always higher except for the voltage limits of 1.01 p.u. and 1.02 p.u. For the voltage limit 1.05 p.u. and 1.06 p.u., as well as for 1.07 p.u. and 1.08 p.u., similar HC distributions are observed. Compared to the VoltVAr control (DER-scenario uniform), a lower 95th percentile is reached, except for 1.01 p.u. The median HC of both control schemes are comparable for most of the voltage limits. However, regarding the median HC, the VoltVAr control outperforms the PF = 0.9 for higher voltage limits. At an admissible voltage limit of 1.04 p.u., the median is at a comparable level of 120 kW.
Assuming that only one DER is installed in the feeder and is located at the feeder end node, the HC is significantly lower compared to the DER-scenario uniform (
Figure 18c). For example, assuming an admissible voltage level of 1.03 p.u., the median value is less than half of the median value for the respective DER-scenario uniform. Further the 95th percentile above voltage limits of 1.04 p.u. remains almost at the same level. However, the 5th percentile and the median are increasing. In addition, in this scenario, a voltage limit of 1.01 p.u. cannot be recommended at all. For the DER-scenario eof, for none of the voltage limits a median value of 100 kW is reached. The calculated HC for this DER-scenario is the minimal HC of the feeder that can be deployed at any node without violating voltage constraints. If the main path from the substation to the end node contains the cable with the lowest ampacity in the feeder, the loading constraint will not be violated as well. Therefore this scenario could be used as benchmark level for the upstream grid to answer the question which amount of DER could be deployed in all supplied feeders without network reinforcements or reactive power control strategies.
In
Figure 18d the results for the DER-scenario uniform and VoltVAr-control are shown. Similar to the uncontrolled DER-scenario uniform, a steadily increase of the HC is monitored such that the median HC of the VoltVAr scenario is higher for all voltage limits. However, the HC value of the 99th percentile is not reached (loading constrained feeders). Similar to the DER-scenario uniform, a voltage limit of 1.01 p.u. cannot be recommended. For the voltage limit 1.04 p.u., the median value even reaches 120 kW.
Based on the median HC findings for the DER-scenario uniform in
Figure 18, the following suggestions can be given. In most of the cases, extending the voltage band, if possible, by 0.01 p.u., leads to a higher HC value than applying either a PF = 0.9 or a VoltVAr control. At an admissible voltage limit of 1.05 p.u., a PF = 0.9 or a VoltVAr control leads to a similar HC value compared to the HC at an admissible voltage limit of 1.06 p.u. (PF = 1). To increase the median HC at 1.06 p.u. (with PF = 1), increasing the voltage band to 1.07 p.u., or applying a VoltVAr-control (at 1.06 p.u.) can be suggested. At an admissible voltage limit of 1.07 p.u., a VoltVAr-control leads to a higher median HC compared to a voltage band extension to 1.08 p.u. The 5th percentile of the HC with reactive power (PF = 0.9 or VoltVAr—
Figure 18b,d) at 1.04 p.u. equals the 5th percentile of the HC at 1.05 p.u. in
Figure 18a. Hence, below 1.05 p.u., an extension of the voltage band leads to a higher HC compared to applying a reactive power control. Above 1.05 p.u., the 5th percentile of HC with reactive power control is equal or higher compared to the uncontrolled case with an increased admissible voltage rise. The HC values of the 95th percentiles of the scenarios with a reactive power control are lower compared to the uncontrolled case (PF = 1). Hence, in a first step, it is suggested to increase the admissible voltage rise, if possible, before applying a reactive power control. Due to the increasing share of loading constrained feeders at higher admissible voltage limits, higher HC values can be reached with a VoltVAr-control compared to a PF = 0.9 control.
5. Correlation Results of the Hosting Capability and Particular Parameters of Feeders
In this section, the results of the correlation between feeder parameters and the introduced hosting capability of feeders for a high number of scenarios are presented.
For all parameters summarized in
Table 1 the correlation with the hosting capability scenarios was calculated. This approach allows identifying parameters that may be used to estimate the hosting capability without load flow calculations, and are relevant for feeder taxonomy studies.
The correlation was calculated for the DER scenarios uniform, eof and additionally for the DER-scenario uniform with activated VoltVAr-control. For each scenario, the hosting capability for a specified allocated voltage band the correlation was performed (
Table 3). Further, for the DER scenarios uniform and eof a power factor of 1 and 0.9 was applied (
Table 4). Thereby, the combination VoltVAr and HC @ PF = 0.9 is not possible: VoltVAr-results are replicated to the lower left part of the figure.
In
Figure 19, the Pearson correlation of all available parameters for two DER-scenarios and the VoltVAr control and voltage limits with the calculated hosting capability is depicted according to the description in the previous tables (
Table 3 and
Table 4). Thereby, the absolute correlation coefficient between zero and one is depicted. In this case, a value of one equals perfect positive or negative correlation and a value of zero indicates uncorrelated parameters. Each column shows the correlation coefficient for all DER-scenarios and voltage limits with one particular parameter with and without a fixed cosϕ-control. This plot demonstrates that, none of the parameters has a high correlation with the hosting capability. The highest correlation coefficient can be found for the scenario PF = 0.9, an admissible voltage limit of
umax = 1.08 p.u. and the DER scenario eof. In this particular case, the correlation coefficient of the related hosting capability and the parameter Inom
Min (cable in the feeder with the lowest ampacity) is 0.699. At a voltage level of 1.08 p.u., a higher share of the feeders is loading constrained, which explains the relatively high correlation coefficient compared to the other parameters.
In addition to the Pearson correlation, the Spearman correlation (rank correlation) was also performed for the same input data (
Figure 20). For the DER-scenario eof and an admissible voltage limit of 1.01 p.u., perfect correlation is observed for dvdP,
ZK and
RK. This perfect correlation for this particular scenario is expected, since at an admissible voltage limit of 1.01 p.u. for the DER-scenario eof, almost all feeders are voltage constrained. Then, the hosting capability is determined by the resistance at the end node. Consequently, with higher admissible voltage limits, the share of loading constrained feeders is increasing and the correlation coefficient decreases. Further investigations on feeder taxonomy are ongoing.
Apparently, according to
Table 2, the most correlated parameters
RK,
ZK and dvdP for voltage constrained feeders also follow a Weibull distribution for both DSOs. With each of these parameters, the voltage constrained HC can be estimated for feeders.
6. Conclusions
The presented results show that the developed simulation platform is a powerful tool to analyze the entire LV-grid data of DSOs on feeder level. Due to the modular script design, the calculation can be highly parallelized, keeping the execution time of scripts even for a high number of networks in an acceptable time frame. The platform is capable to calculate the maximal hosting capability per feeder for a high number of scenarios. Moreover, the presented simulation platform was used to identify which constraint limits the hosting capability of the feeders, and whether options other than network reinforcement to increase the hosting capacity are reasonable [
35].
To run power flow studies across high-, medium- and low-voltage grids, another option would be to reduce the number of nodes. Therefore, simplified LV-network models are needed to reduce the computational burden. In a future work, simplified LV-feeder models will be defined and validated with the assessed networks. Apparently, out of all investigated parameters, the InomMin (Pearson correlation and higher admissible voltage limits) and the end node resistance (Spearman correlation and lower admissible voltage limits) are the parameters most correlated to the introduced HC, which need to be considered.
The hosting capability results were presented for the DER-scenarios uniform and eof for eight admissible voltage limits ranging from 1% to 8% of the nominal voltage. Moreover, the HC was evaluated with and without reactive power control strategies. In particular, the VoltVAr-control and a fixed cosϕ-control were investigated additionally to the uncontrolled case. Due to the significant HC increase from 1.01 p.u. to 1.02 p.u., an admissible voltage rise limit of 1% cannot be recommended at all. The HC results showed, that increasing the admissible voltage rise for in-feed is only reasonable if a high share of feeders remain voltage constrained. If an admissible voltage rise of 4% is acceptable, for the DER-scenario, uniform more than 50% of the feeders have a HC higher than 100 kW. With a VoltVAr-control, even a median value of 120 kW can be reached. In the worst case scenario, DER-scenario eof, significantly lower median values are achieved. Further, utilizing a reactive power control for loading constrained feeders, leads to a significant HC reduction.
The results on benchmarking showed that different indicators (e.g., minimum spanning tree (without costs) or weighted transformed connection density) could be used to benchmark networks. However, investigations put some bias in evidence when comparing very small and very large networks. These very high-level indicators should therefore be carefully used and further investigations with more DSOs would allow drawing more general conclusions.
The correlation of feeder parameters with the calculated hosting capability values showed that for most of the scenarios no high correlation coefficient could be found. None of the calculated hosting capability values correlate highly with statistical parameters used in literature to classify networks. Further, only a few parameters follow a probability distribution for both investigated DSOs. The described probability plots for selected parameters and feeder statistics however may be used to generate synthetic realistic feeders. Furthermore, the hosting capability results matched with solar resource maps could enable cost-efficient DER deployment plans without investments in grid reinforcement on LV feeder level.