The simulated PV hosting capacity scenarios.
Abstract
Maximizing the share of renewable resources in the electric energy supply is a major challenge in the design of the future energy system. Regarding the low voltage (LV) level, the main focus is on the integration of distributed photovoltaic (PV) generation. Nowadays, the lack of monitoring and visibility, combined with the uncoordinated integration of distributed generation, often leads system operators to an impasse. As a matter of fact, the numerous dispersed PV units cause distinct power quality and costefficiency problems that restrain the further integration of PV units. The PV hosting capacity is a tool for addressing such power system performance and profitability issues so that the different stakeholders can discuss on a common ground. Photovoltaic hosting capacity of a feeder is the maximum amount of PV generation that can be connected to it without resulting in unacceptable power quality. This chapter demonstrates the usefulness of smart metering (SM) data in determining the maximum PV hosting capacity of an LV distribution feeder. Basically, the chapter introduces a probabilistic tool that estimates PV hosting capacity by using customerspecific energy flow data, recorded by SM devices. The probabilistic evaluation and the use of historical SM data yield a reliable estimation that considers the volatile character of distributed generation and loads as well as technical constraints of the network (voltage magnitude, phase unbalance, congestion risk). As a case study, an existing LV feeder in Belgium is analysed. The feeder is located in an area with high PV penetration and large deployment of SM devices.
Keywords
 distributed generation
 low voltage network
 probabilistic simulation
 smart meters
 photovoltaic
1. Introduction
Maximizing the share of renewable resources in the electric energy supply is a major challenge in the design of the future energy system. Concerning the low voltage (LV) distribution system, this objective aims at increasing the selfsufficiency of LV feeders, based on local resources, while responding to the climate change. In such feeders, distributed photovoltaic (PV) generation is the mostly met distributed energy resource (DER).
So far, the biggest share of distributed PV units came with no previous planning or reinforcement of the network while monitoring data in the residential or commercial sector were absent almost everywhere in Europe. Given the lack of controllability in common LV networks, the uncoordinated integration of PV units often leads to distinct power quality issues. As a result, the connection of new PV units and therefore the increase of renewable energy share slows down. Adding to this fact the growing volatility of electricity consumption in the distribution network, the adoption of a streamlined planning approach for the future energy system becomes urgent.
In this evolving framework, distribution system operators (DSOs) are called to safeguard a stable and secure power supply in all possible demand conditions while fostering the massive integration of DER generation. In costefficiency terms, this fact highlights the necessity of leaving behind deterministic worst case planning approach. This traditionally applied approach focuses on the least favourable network operation states, which are very rare. Naturally, it leads to very restrictive decisions in terms of PV hosting capacity or to costly network reinforcements.
Given the current uncertainty of DSO costs and revenues, new planning tools are required for considering the constant variability of the energy network [1]. This argument becomes even more solid in view of the upcoming integration of electric vehicles and the development of flexibility services. As a matter of fact, both are seen as basic components of the future energy model. The large deployment of smart metering (SM) devices in the residential and commercial sector will drastically enlarge the potential of costeffective planning approach. Indeed, customerspecific data will result in a better insight of the power distribution system.
Considering the above facts and the probabilistic character of the EN 50160 technical standard [2], [3] (which addresses the LV network), this chapter presents a feeder and customerspecific probabilistic method that estimates the DER hosting capacity of an LV feeder. Practically it introduces a probabilistic tool that uses customerspecific energy flow data recorded by SM devices (that are installed in the studied feeder), over a period of 2 to 4 years depending on the customer. The deployment and probabilistic elaboration of longterm energy measurements (SM readings) yield a reliable estimation that considers the volatile character of distributed generation and loads as well as several operational metrics.
Section 2 of this chapter outlines some of the existing scientific contributions that address this subject and presents the drivers for developing the proposed methodology and algorithm. Section 3 presents the overall structure of the developed algorithm and Section 4 thoroughly describes the primordial contribution of customerspecific SM readings in this development. Section 5 explains the computation process of the maximum acceptable PV hosting capacity.
As a case study, a real LV feeder in Belgium is analysed. The feeder is located in an area with high PV penetration and large deployment of SM devices. When the probabilistic character of EN 50160 standard’s voltage limits is considered, the estimated PV hosting capacity is proved to be much higher than the one obtained with a deterministic approach, based on worst case energy flow profiles. At the same, the use of the SM readings verifies the computation of the technical constraints in the feeder that cannot be considered with a probabilistic approach (violation of the maximum current capacity of the lines).
2. Current Framework
In many regions worldwide, DER integration is hampered due to slow or over rigid hosting capacity review processes. As a result, customers who want to invest and play an active role in managing their energy usage are increasingly unable, in expediency and costefficiency terms, to do so. In this context, a streamlined approach together with the expansion of allowable DER integration approvals seem to be a necessity [4].
However, the expansion of allowable approvals depends heavily on DER admissible penetration levels, which are determined by local DSOs. In order to increase penetration levels while facilitating the application review process, DSOs should incorporate automated DER hosting capacity analyses. An example process flow for incorporating such analysis into the DER integration review process is outlined in Figure 1.
Currently, many energy utilities are adapting their DER hosting capacity review so as to remove or update restrictive maximum allowable limits [5]. To this end, the Electric Power Research Institute (EPRI) presents a set of models that could be used by DSOs or electric utilities [6], [7]. These feederbased methodologies are very solid computation examples that take account of all steady state operational criteria.
Focusing on PV hosting capacity, the EPRI’s report presents stochastic analysis as a highly appropriate tool for determining feeder hosting capacity for distributed PV units. The stochastic deployment concerns the position and size of future PV units while the steady state estimation of the feeder is done with deterministic approach. Indeed, the analysed state estimation scenarios are based on four worst case load/ PV generation profiles.
In the same vein, a set of studies addressing the European framework and the EN 50160 standard highlight the efficiency of stochastic and probabilistic analysis in determining hosting capacity or otherwise the impact of PV generation in LV feeders [5]–[12].
Meanwhile, the European Photovoltaic Industry Association (EPIA) and the technical standard EN 50160 suggest that distribution networks should be designed on a probabilistic basis. For example, EN 50160 standard deals with the voltage characteristics of LV feeders in probabilistic terms. It gives recommendations that, for a percentage of measurements (e.g. 95%) over a given time, the voltage value must be within specified limits.
Most of the existing methodologies deploy the stochastic analysis regarding the size and position of PV units and not the load/generation profiles of customers. However, the ongoing integration of SM devices in LV networks enlarges the potential of using feederspecific or even customerspecific data for modelling energy flows. According to [9], performing longterm measurements in the LV network is highly valuable and strongly recommended, not only for estimating the maximum PV hosting capacity, but also for voltage coordination of the network in general. Measurements of the voltage magnitude for a large time of customers may be timeconsuming and expensive; however, many countries have already been installing energy meters that also allow the recording of voltage, current, active power and other metrics.
Considering these facts, the EPRI’s report [7] estimates PV hosting capacity using feederspecific data to create either absolute worst case scenarios (maximum recorded generationminimum recorded load) or load/PV timeofday coincident worst case scenarios. Therefore, feederspecific data are indeed used; however, the steady state estimation of the feeder is still done with a deterministic approach. Consequently, this approach does not consider the fact that the timeofday in which worst case values apply for a specific customer does not necessarily coincide with the one of other customers connected to the same feeder. Nevertheless, the operational criteria of the feeder are determined both by the individual user’s demand and by the simultaneous demands of other network users. Since the demands of every user and the degree of coincidence between them constantly vary, so does the operation of the feeder [3].
The above argument demonstrates that although customerspecific SM data are primordial for creating reliable network models, there is another challenge that needs to be addressed. The latter lies in the fact that customers follow volumewise (kWh) or capacitywise (kW) an almost stable daily pattern. However, this pattern does not necessarily remain the same on the time axis. In longterm decision making, profiles should be based on the recorded ones considering all possible deviations. Those deviations could be inserted either as random statistical errors or by making random possible combinations of the recorded values or even by combining both approaches.
Consequently, reliable models that use customerspecific real SM readings and take into account load/PV time and customervariability are necessary for applying a less conservative and more costeffective hosting capacity review. Probabilistic and particularly Monte Carlo approach are very suitable to address this modelling challenge.
3. The PV Hosting Capacity Computation Tool
Hosting capacity is defined as the maximum amount of PV that can be accommodated in the feeder without impacting system operation (reliability, power quality, etc.) under existing control and infrastructure configurations [7]. This chapter presents a tool that uses probabilistic state estimation, 15min customerspecific SM energy flow readings and feederspecific technical parameters to estimate the PV hosting capacity of a given LV feeder.
The proposed methodology aims to address the central block of Figure 1 (“
Although the EN 50160 standard sets the same voltage limits in all European countries (except from cases where stricter limits are locally imposed), the maximum line capacity heavily depends on the respective DSOs. In certain countries, line sections are chosen based on a longterm strategy that aims at minimising voltage and congestion risk even if loads and generation increase importantly in the future. However, such approach leads to higher initial investment which is not necessarily costeffective. In other cases, line sections are chosen based on actual conditions or short term future scenarios so that customised solutions are applied as soon as problems arise.
Apart from steady state constraint management, there are other considerations that could be accounted for, such as transformer aging factor, line losses, etc. Such criteria are usually considered in an overall costbenefit analysis (CBA); however, at present they are not addressed by the EN 50160 standard. Depending on the country and the applied DSO tariff methodology (“costplus”, “revenue cap”, etc.), DSOs are incentivised to reduce certain operation costs that can or cannot be integrated in their tariffs. Thus, the impact of such criteria on decision making, varies in function of the distribution utility. For this reason, this chapter computes PV hosting capacity focusing on commonly adopted EN 50160 standard criteria and line capacity issues. Line losses in the feeder, during PV injection hours, are also addressed however their rise is not imposed as a constraint to the further increase of admissible hosting capacity.
3.1. Overview of the simulation tool
As previously said, this chapter presents a probabilistic algorithm that determines the PV hosting capacity of an LV feeder by elaborating feederspecific SM measurements. The SM measurements are the necessary input for performing a reliable steady state analysis of various possible energy flow scenarios in the studied feeder. The flowchart in Figure 2 presents the structure of the simulation algorithm, which is entirely developed in MATLAB®. The energy exchange scenarios are generated by the Monte Carlo algorithm sampling from the historic SM data of the feeder [13], [14]. The power flow analysis is performed with the threephase algorithm that is presented in [14] and outlined in Appendix A. Both balanced and unbalanced situations can be considered in this study.
3.2. Feeder model
The feeder model is constructed based on the technical parameters of the lines, the position of the customers, the installed PV power per node, the voltage at the MV/LV transformer secondary output and the respective set points and bandwidths in case voltage control algorithms are integrated. The feeder model also assigns the load/PV generation SM datasets to the respective customers. This necessary information is available to the DSO.
Regarding the PV hosting capacity computation, the possible future locations of the PV units have to be specified in the feeder model. This analysis is not based on stochastic random distribution of PV units along the feeder. A set of scenarios regarding the positions of future PV nodes is specified and each one of them is studied separately so as to focus on its specific impact on the feeder, which is possible thanks to the customerspecific SM datasets.
The technical constraints that must be respected for the current situation and for future scenarios are the ones specified in local, regional or national directives. However, these operational constraints can be determined in a more restrictive manner, depending on the case. In the EU framework, the steady state constraints are set by the EN 50160 standard. Regarding voltage magnitude and unbalance, 95percentile limits are suggested. For considering this EU standard, the simulation tool verifies that the following criteria apply for the whole system (in current and future installed PV power scenarios):
where
The thermal limits of the cables are also considered in the computation. The current carrying capacities of the lines should not exceed the DSO requirements or the recommended values in technical standards such as [15]. The load flow analysis of each system state is performed with the threephase algorithm that is explained in the Appendix [14].
4. The use of SM measurements
4.1. Customer profiles and feeder state modelling based on historic SM datasets
The load/ PV profiles of existing customers are created by using their respective SM recorded datasets. Practically, each dataset consists of values of 15min PV energy injection to the grid (
where
The two 15min resolution datasets of
The TDP of energy consumption
As a next step, the customerspecific TDPs are statistically transformed to 2 x 96 Cumulative Distribution Functions (CDFs) of Probability, one for PV generation and one for energy consumption (=2) multiplied by the number of quartershourly time steps in a day (=96), representing each one of the studied months. This transformation is made by applying the basic statistical formula, for each 15min dataset:
Therefore, the “typical day” of each customer (for the respective month) can be illustrated by two diagrams like the one presented in Figure 4, one for PV injection and one for net energy consumption. Exactly the same methodology is applied to build TDPs for the r.m.s. voltage at the secondary output of the MV/LV transformer, using 15min data, also recorded at the specific MV/LV substation.
The created statistical distributions (CDFs) of the time varying parameters (PV injection, net energy consumption and voltage at the MV/LV transformer) are used for defining multiple possible network states corresponding to each 15min time step. Practically, a MC algorithm is applied for randomly sampling the values of the variable parameters at each node of the studied LV feeder as explained in [13], [14] . The combination of the nodal sampled values defines each network state that will be afterwards analysed by the power flow algorithm (Appendix). Subsequently, for each state, a feeder with
where
The power flow analysis of the feeder requires considering each system state as instantaneous, and therefore the sampled energy values have to be transformed into instantaneous power values (
If
The probabilistic deployment of this simulation tool relies on the principle that load/PV generation profiles of customers are highly timevarying. The generation of the system states is therefore based on a very large number of random combinations of customers’ energy flow values. This timevariability induces another variability that concerns the time coincidence of the load profiles of various customers. Both arguments are very important when assessing the impact of PV generation on a LV network. Indeed, the consideration of this variability, both in the time axis and regarding customers coincidence, makes more realistic the simulation of the network operation. Such an approach can lead to less restrictive and more costeffective decisions that do not rely on rare extreme cases but on the most frequent ones.
4.2. Generation profiles of future PV nodes
A key component in accurately assessing the impact of future PV units is reliably representing their generation profiles. Based on the findings of several studies, geographically close customers are entirely correlated as far as their PV generation profiles are concerned [16]. For this reason, this study considers that the generation profiles of future PV customers will be very similar, along the time axis, to the ones of the existing PV units.
As previously explained, the load/PV generation profiles of customers with SM devices are made of 96 Cumulative Distribution Functions (CDFs) of probability built with the 15min recorded datasets. Concerning PV generation, such CDFs are apparently not available for the future PV units. For this reason, the available SM datasets are used in this case to create a reference CDF, based on the 15min generation SM datasets of the existing PV owners [17]. This reference CDF is used to simulate the timevariability of PV generation at the future PV nodes.
In reality, customers that are connected to the same LV feeder can have different PV units’ sizes. Assuming an equivalent statistical distribution of their PV power profiles due to geographical proximity, the principle is to create a standardized reference CDF for PV generation in the specific feeder, based on the measurements of the available SM devices. Initially, the CDF for the 15min PV energy generation
where
Once this is done, the 15min CDFs of every user are aggregated, as graphically outlined in Figure 5, in order to create one reference CDF that can represent all PV owners in the specific feeder. For creating the CDF of each particular future PV owner, this reference CDF should be normalised in function of his annual PV generation. For existing PV owners, such information is usually available to the DSO even if the customer is not monitored by an SM device. In case of future PV nodes, such information is apparently not available since no PV unit is connected. Consequently, the reference CDF is normalised with the annual PV generation of an existing PV unit (in the feeder or in proximity) multiplied by a reference factor
5. PV hosting capacity computation
Practically, the algorithm starts with the probabilistic analysis of the current situation (existing PV units), by simulating a large number
The probabilities
where
Once relation (6) is applied, the new installed PV power
where
Once the generation profiles have been set up for the future PV nodes, the algorithm repeats steps 2 and 3 for analysing the current configuration
The described iterations stop as soon as the operational constraints are for the first time exceeded at least at one of the nodes. The PV size of some units could probably increase even more, given that the operational constraints at their PCC are not violated. However, this study treats the LV feeder as a whole since the violation of limits at one node is always affected by the energy flow at all nodes. The
The aggregated PV hosting capacity of the feeder is computed by adding
where
6. CASE STUDY: An LV feeder in Belgium
6.1. Description of the simulation
This section describes the application of the previously described analysis tool for computing the PV hosting capacity of an LV feeder in Flobecq. Flobecq is a municipal area in Belgium with high penetration of distributed PV generation (≈25% of Flobecq LV network customers) and large deployment of SM devices. Thanks to an official research fellowship between the local DSO and the authors’ affiliation, the technical parameters of the feeder and SM datasets of the respective customers have been communicated strictly for research purposes. The datasets that were used in this case study cover a total period of one year (2013).
The topology of the simulated threephase feeder is presented in Figure 6. Currently, four PV units are installed in the feeder which supplies a total of 16 residential customers. These PV units are located at nodes 4,5,12 and 14, by means of singlephase inverters, and their installed PV power is respectively 5kVA, 10kVA, 2.63kVA and 5kVA. A spatial correlation study had already been performed for the specific feeder and the generation profiles of the customers were proved to be entirely correlated [16]. This consideration is taken into account in this analysis, including the future PV nodes. Practically, this means that for every simulated system state, the randomly sampled probability value for defining PV generation is common for all PV units.
Concerning operational constraints, the ones of EN 50160 standard have been considered in the simulation. Therefore, compliance with the group of equations (1) has been verified for each system state, as far as voltage magnitude and unbalance are concerned. The maximum current capacity of the lines has been determined based on table [15]. The PV size increase step is defined equal to 1kVA and the power factor of all PV inverters is considered equal to 1, unless reactive power control is considered in the simulation.
A set of different scenarios have been simulated regarding the position and phase connection of future PV units as well as the action of voltage control schemes. The analysed scenarios are listed in Table 1.
Concerning scenarios AD, only the onoff control scheme is considered, which is currently implemented by most DSOs in Europe. This control scheme enables a total cutoff of the PV unit (in most cases during 3 minutes) as soon as the voltage limit has been locally exceeded for a period longer than 10 minutes. This analysis considers each simulated state as instantaneous. Therefore, each violation of the 95percentile limit of EN 50160 standard is counted in the probabilities even though in reality it might had lasted less than 10 minutes. This means that the computed maximum PV hosting capacity is possibly slightly lower than the one that the feeder can really support.
No  Description 

A  12 new PV units at nodes 2, 3, 6, 7, 8, 10, 11, 13, 15, 17, 18, and 19. The PV units at nodes 8, 11, 17, 18, and 19 are connected to phase A, the PV unit at node 3 is connected to phase B and the PV units at nodes 2, 6, 7, 10, 13, and 15 are connected to phase C. 
B  Similarly to scenario A but all new PV units connected to phase B, except from PV unit at node 15 that is connected to phase A. 
C  3 new PV units connected to nodes 13, 18 and 19 (end of the line). All three new units connected to phase B. 
D  1 new PV unit connected to node 13 (phase B). 
E  Similarly to scenario A but considering 100percentile and not 95percentile operational limits. Practically the PV hosting capacity is not increased as soon as voltage and VUF limits are exceeded at least once in the feeder. 
F  Similarly to scenario A but considering the action of threephase damping control integrated in the new PV inverters. In this case, the new PV units need to be connected by means of threephase PV inverters. 
G  Similarly to scenario A but considering the action of reactive power control of (CEI 2012) 
The control scheme applied in scenario F is the threephase damping control scheme which behaves resistively towards the negative and zerosequence voltage component, without modifying the injected power, so as to eliminate phase voltage unbalance [18]. This control scheme requires a threephase PV inverter and it is very promising in terms of voltage magnitude and unbalance mitigation. It is actually implemented in a EU pilot program (FP7 INCREASE Project). The third control scheme is reactive power control in the way it is implemented in the Italian distribution system [20] concerning new PV units in the LV network.
6.2. Comparing with a deterministic approach
One of the main purposes of this study is to investigate, up to which extent, a probabilistic method based on customerspecific SM readings leads to a less restrictive computation of PV hosting capacity, compared to a deterministic approach. For this purpose, a deterministic approach has been implemented simulating worst case energy flow profiles. The load profiles of all customers and the PV generation profiles of existing PV units have been also based on SM recorded data. The deterministic steady state analysis has been conducted for scenarios AD, F, G. Scenario E is not mentioned because, although SM readings are used, only 100percentile limits are considered which means that probabilities are not accounted for in the computation of hosting capacity. Thus, this scenario is practically a deterministic scenario.
The following load/ PV generation profiles have been considered in the deterministic approach:
Maximum PV power per node (installed PV power) – Minimum recorded load per node; absolute values, irrespective of time coincidence among customers
Maximum PV power recorded in the feeder – Coincident PV generation/load values for the other nodes.
Minimum recorded load in the feeder during PV injection hours – Coincident PV generation/load values for the other nodes.
6.3. Results and discussion
The probabilistic hosting capacity review results are illustrated in Figure 7 and analytically listed in Table 2. The aggregated maximum admissible PV hosting capacity in the feeder is presented for each individual scenario, considering separately the EN 50160 standard’s voltage limits and the maximum current capacity of the lines. This separate presentation has been chosen because voltage limits are treated with a probabilistic approach in the EN 50160 standard while congestion risk is treated by each DSO with a different approach. Most of them apply a deterministic approach that considers an upper (100percentile) current limit. The violation which forbade further increase of the PV hosting capacity is also presented and quantified for each scenario. The aggregated PV hosting capacity obtained with deterministic analysis is presented in Figure 7 and Table 3 for all treated scenarios (§5.2).





Aggregated PV hosting capacity (kVA)  Violation  Aggregated PV hosting capacity (kVA)  Violation  
A  154.63kVA (11kVA/ per new PV node + existing PV units) 
 and 19 (phase (B)) resulted 5.7% and 6.4% respectively (>5%, which is the value accepted by the EN 50160 standard) 
70.63kVA (4kVA/ per new PV node + existing PV units) 

B  144.63kVA (10kVA/ per new PV node + existing PV units) 
 13,14,15 (phase (C), resulted 5.4%, 6.16% and 6.18% respectively (>5%, which is the value accepted by the EN 50160 standard) 
58.63kVA (3kVA/ per new PV node + existing PV units) 
(50% higher than the maximum current capacity of the lines) 
C  82.63kVA (20kVA/ per new PV node + existing PV units) 
 19 (phase (B)) resulted 6.3% (>5%, which is the value accepted by the EN 50160 standard) 
43.63kVA (7kVA/ per new PV node + existing PV units) 
(10.5% higher than the maximum current capacity of the lines) 
D  65.63kVA (43kVA/ per new PV node + existing PV units) 
 19 (phase (B)) resulted 5.15% (>5%, which is the value accepted by the EN 50160 standard) 
37.63kVA (15kVA/ per new PV node + existing PV units) 
(6.2% higher than the maximum current capacity of the lines) 
E  94.63kVA (6kVA/ per new PV node + existing PV units) 
 13,14, (phase (C)), resulted 0.0001% in both cases (>0%, which is the condition in scenario E) 
70.63kVA (4kVA/ per new PV node + existing PV units) 

F  202.63kVA (15kVA/ per new PV node + existing PV units) 
 219 (at all three phases) resulted from 5.5% to 28% (>5%, which is the value accepted by the EN 50160 standard)  
70.63kVA (4kVA/ per new PV node + existing PV units) 

G  154.63kVA (11kVA/ per new PV node + existing PV units) 
 19 (phase (B)) resulted 5.17% (>5%, which is the value accepted by the EN 50160 standard) 
70.63kVA (4kVA/ per new PV node + existing PV units) 


(considering EN 50160 standard voltage limits) 







A  70.63kVA (4kVA/ per new PV node + existing PV units) 
82.63kVA (5kVA/ per new PV node + existing PV units) 
82.63kVA (5kVA/ per new PV node + existing PV units) 
Overvoltage at all new PV nodes 
B  58.63kVA (3kVA/ new PV node + existing PV units) 
58.63kVA (3kVA/ new PV node + existing PV units) 
58.63kVA (3kVA/ new PV node + existing PV units) 
Overvoltage at all new PV nodes 
C  43.63kVA (7kVA/ new PV node + existing PV units) 
43.63kVA (7kVA/ new PV node + existing PV units) 
43.63kVA (7kVA/ new PV node + existing PV units) 
Overvoltage at all new PV nodes 
D  37.63kVA (15kVA/ new PV node + existing PV units) 
37.63kVA (15kVA/ new PV node + existing PV units) 
37.63kVA (15kVA/ new PV node + existing PV units) 
Overvoltage at all new PV nodes 
Firstly, one should note that the results considering the maximum current capacity coincide for the probabilistic and the deterministic computations since this metric is not addressed with probabilistic terms. In case the maximum value is exceeded at one point of the feeder during just one of the simulated states, the PV hosting capacity is not further increased in the respective simulation. Also, the result of scenario E (applying 100percentile limits) is close to the ones of the deterministic scenarios A.I, A.II and A.III which analyse the same topology as scenario E but with a deterministic approach. Based on these remarks, one can reasonably assume that the probabilistic computation covers (samples and analyses) almost the whole range of possible system states, including the ones recorded in reality (the combination of coincidently recorded values) which are treated in the deterministic scenarios A.II and A.III.
However, accounting only for voltage violation, the restrictive condition of scenario E based on which voltage limits must never be exceeded (in none of the simulated states), results in a quite lower admissible PV hosting capacity compared to scenario A (same topology as scenario E). Basically, in scenario E, PV hosting capacity could not further increase because the computed
Figure 9 clearly demonstrates how the probabilistic consideration of overvoltage risk affects the computation. This figure shows the evolution of the CDF of phase voltage (B) at node 19 while the total installed PV power increases in scenario A. Based on the probabilistic analysis of the feeder, when total installed PV power increases by 144kVA (12kVA per new PV unit), phase voltage (B) at node 19 respects the defined limits in 94.6% of the simulated states. However, EN 50160 defines that the limits should be respected in at least 95% of cases. Thus, the maximum PV power that can be added to the feeder, considering this configuration, is 132kVA (11kVA per new PV unit).
The above arguments should be considered in a costbenefit analysis (CBA) that compares costs for exceeding operational limits to losses due to an eventual penalty for low DER integration or loss of potential revenue for customers and energy utilities. Such considerations may allow a much more costeffective PV integration strategy which also respects the applied standard’s criteria. At the same time, the identification of critical points regarding congestion risk should also be considered. Both arguments highlight the usefulness of considering SM historic datasets in similar studies, so that critical points and probabilities are carefully mapped and quantified.
For highlighting the costeffectiveness of deploying longterm measurements in the LV network and analysing it with a probabilistic approach, a more detailed computation of line losses in the feeder was performed for scenario A. Assuming that the computed maximum admissible PV power is installed (=154,63kVA if one considers only the voltage limits), the study focuses on the total energy losses along the lines of the feeder during hours of high PV injection in a typical day. The worst case approach considers only one system state which will more likely take place during hours with the highest PV injection. Based on the available historic data for the feeder, this period is between 12:00AM and 18:30PM on a typical July day. The sum of energy losses has been computed along the feeder for the considered period, for each simulated day. Figure 10 illustrates the statistical distribution (CDF) of the computed line losses, obtained with the probabilistic approach.
The probabilistic approach and the consideration of the SM measurements demonstrated that total energy losses in the feeder vary significantly, depending on the system state. Consequently, in 95% of the simulated days, total energy losses during high PV injection hours (12:00AM to 18:30PM) do not exceed 35kWh in a day. In the deterministic approach which assumes the worst case scenario taking place all along the high PV injection period, the respective energy losses result equal to 148kWh. This important difference is due to the fact that the probabilistic approach considers the extremely low frequency of worst case scenarios to take place simultaneously for all feeder users. Considering such probabilities, the DSO could manage a less conservative and more costeffective longterm strategy.
Undoubtedly, the computed PV hosting capacity values depend on the load profiles of the customers that are located in the feeder. However, the results clearly indicate in relative terms, that smaller distributed PV units have a much smoother impact than the bigger ones concentrated in one small area of the feeder. This fact is demonstrated by the comparison of scenario A to scenarios C and D. Moreover, as previously mentioned, in several cases the maximum admissible installed power per PV unit connected to the LV network is equal to 10kVA. In such cases, scenarios C and D might not be appropriate based on the probabilistic simulation results. As a matter of fact, the admissible total installed power would have to limit to 32.63kVA (for scenario D) although the network would be able to support 37.63kVA. The difference between the PV hosting capacity computed with the probabilistic and the deterministic approach for these cases (considering only voltage limits) is not as big as for scenarios A and B. Indeed, in scenarios A and B, the volatile character and the extremely rare coincidence of worst case values for 12 units cannot be reliably represented by a deterministic model.
Regarding the distribution of units among phases, the comparison of scenarios A and B shows that the existing phase unbalance affected the computation. Indeed, the violated parameter in this case is voltage magnitude of phase (C) although all new PV units are connected to phase (B). Therefore, the unfair distribution of new PV units among phases did not directly affect
Scenario F demonstrated that the connection of new PV units by means of threephase inverters integrating threephase damping control can increase the aggregated hosting capacity by 36%, if the probabilistic character of EN 50160 standard’s voltage limits is considered. Thanks to the resistive behaviour of this control scheme towards the zero and negativesequence voltage component, the deviation of voltage magnitude and unbalance becomes much smoother compared to the currently applied onoff control. Thus, the risk of exceeding the defined limits is reduced and a bigger share of PV generation can be integrated.
Based on the results of scenario G, reactive power control does not result in higher PV hosting capacity compared to scenario A (onoff control). Voltage profile in the feeder is however improved compared to scenario A. As a matter of fact, voltage limits are not violated in scenario G whereas the maximum current capacity limit is exceeded for the same amount of PV integration compared to scenario A.
In the first two cases (scenarios A and B), comparing the probabilistic simulation results to the respective ones of the deterministic approach, an important difference in the aggregated admissible hosting capacity is observed. At this point, it is important to mention that the violated parameter in the deterministic approaches is mainly the voltage magnitude and secondly the maximum current capacity of the lines. The deterministic approach led to 74146% lower aggregated PV hosting capacity (compared to the one computed with the probabilistic approach) due to a violation that according to the probabilistic elaboration of the historic SM dataset took place for much less than 5% of the simulated system states. Indeed, based on Figure 8, the addition of 12 new PV units of 4kVA each (result of deterministic scenario A.I) generated an overvoltage risk that is lower than 1%.
In such cases, the probabilistic analysis demonstrates that the simulated worst case scenarios are extremely unlikely to happen in the studied feeder. However, such worst case approach is currently implemented by most DSOs when performing hosting capacity reviews. As a result, decisions for connecting new PV units in certain networks are very often extremely restrictive, with a big impact on the costefficiency of the network.
When it comes to scenarios C and D, one can note that if less but bigger size PV units are connected, the results of the probabilistic and the deterministic approach do not differ significantly. This result proves that in case of many distributed PV units, an approach that considers all worst case customers’ profiles coinciding in time is mostly extreme. Deterministic approach cannot accurately simulate the volatile character of PV generation and the random loading parameters of residential and small commercial customers.
A general remark would concern the design strategy of distribution feeders like the studied one. The studied feeder currently hosts 22.63kVA of distributed PV generation and supplies 19 residential customers. The analysis of the current conditions (based on the historic SM datasets) demonstrated that both voltage violation risk and congestion risk are very low. Moreover, the above probabilistic loadflow analysis demonstrated that congestion and voltage problems will only appear if 48kVA and 132kVA respectively of distributed PV generation (scenario A) are further integrated. This remark highlights the costefficiency of designing distribution networks based on the most frequent system states or on wellstudied future scenarios. This approach can lead to customised solutions and help to avoid overdimensioning and costly initial investments for the DSO.
Based on the above analysis, certain renewable integration scenarios could increase to an important extent the selfsufficiency of feeders like the studied one. As a result, their dependency on big conventional power plants, connected at the transmission level, could be efficiently reduced. However, big conventional plants are important for maintaining grid stability. In a high DER integration scenario, without large and reactive storage facilities and/or flexibility services, the amount of RES should be carefully reviewed. To this end, costs induced by the use of grid services, including insurance against periods when it is not possible to consume own generated electricity, should be considered and reflected in the bill of generator owners [1]. Reliable feasibility studies and comprehensive CBAs are necessary for evaluating various strategies in the decision making process.
6.4. The role of customerspecific SM data in PV hosting capacity reviews in LV networks
The above analysis is based on the use of customerspecific SM energy flow readings. Various maximum PV hosting capacity scenarios have been analysed by applying a probabilistic steady state analysis of the feeder on a 15min time scale, sampling from the available SM data. In this way, the real probability of worst case scenarios has been accounted for and as a result, a probabilistic view of several technical metrics has been enabled (voltage and current magnitudes, voltage phase unbalance, line losses). Even if DSO longterm planning is based on deterministic approach, the use of SM datasets can validate the considered worst case scenarios. Besides, the wide deployment of SM devices can offer other possibilities such as better coordination and control of technical parameters of the LV network as well as better visibility and adaptability to actual load and generation profiles of LV customers. The longterm planning of LV networks can become more customised to local conditions and therefore more costeffective.
7. Conclusions
This chapter addresses the problem of determining the maximum PV hosting capacity that can be accommodated in a LV distribution feeder, while respecting local technical standards. To this purpose, a probabilistic simulation tool that uses as input customerspecific SM energy flow data and feederspecific parameters is presented. A PV hosting capacity review for a municipal area in Belgium is used as a case study for evaluating the usefulness and reliability of the proposed tool. The study outcome demonstrates that it is to the interest of the DSO and of the grid users to deploy probabilistic analysis that considers the timevariability of load/PV generation, both in the time axis and between different customers’ profiles. This variability of network state can be taken into account thanks to the deployment of longterm SM measurements in the studied network. Consequently, the further deployment of SM devices is strongly recommended for achieving a more costeffective longterm planning and coordination of the LV network.
Acknowledgments
The authors of this chapter acknowledge the support of ORES, the DSO who manages the electricity and natural gas distribution grids in 193 communes in Wallonia (Belgium), in terms of funding and SM data supply, both indispensable elements for conducting this research work.
References
 1.
EDSO, “European Distribution System Operators for Smart Grids Adapting distribution network tariffs to a decentralised energy future,” 2015.  2.
EN50160, “Voltage characteristics of electricity supplied by public electricity networks,” 2012.  3.
M. B.P. Antoni Klajn, “Application Note Standard EN50160 Voltage Characteristics of Electricity Supplied by Public Electricity Networks,” March, 2013.  4.
Solar City Grid Engineering, “Integrated Distribution Planning: A holistic approach to meeting grid needs and expanding customer choice by unlocking the benefits of distributed energy resources,” 2015.  5.
J. C. Hernandez, F. J. RuizRodriguez, and F. Jurado, “Technical impact of photovoltaicdistributed generation on radial distribution systems: Stochastic simulations for a feeder in Spain,” Int. J. Electr. Power Energy Syst. , vol. 50, no. 1, pp. 25–32, 2013.  6.
F. J. RuizRodriguez, J. C. Hernández, and F. Jurado, “Probabilistic load flow for photovoltaic distributed generation using the CornishFisher expansion,” Electr. Power Syst. Res. , vol. 89, pp. 129–138, 2012.  7.
P. Zhang, W. Li, S. Li, Y. Wang, and W. Xiao, “Reliability assessment of photovoltaic power systems: Review of current status and future perspectives,” Applied Energy , vol. 104. pp. 822–833, 2013.  8.
R. Billinton and Bagen, “Generating capacity adequacy evaluation of small standalone power systems containing solar energy,” Reliab. Eng. Syst. Saf. , vol. 91, no. 4, pp. 438–443, 2006.  9.
Bollen M.H.J. and F. Hassan, Integration of Distributed Generation in the Power System , IEEE Press. Wiley, 2011.  10.
R. Billinton and R. Karki, “Reliability/cost implications of utilizing photovoltaics in small isolated power systems,” Reliab. Eng. Syst. Saf. , vol. 79, no. 1, pp. 11–16, 2003.  11.
European Photovoltaic Industry Association, “Connecting the Sun. Solar Photovoltaics on the Road to LargeScale Grid Integration,” 2012.  12.
European Photovoltaic Industry Association, “Global Market Outlook for Photovoltaics,” 2014.  13.
F. Vallee, V. Klonari, T. Lisiecki, O. Durieux, F. Moiny, and J. Lobry, “Development of a probabilistic tool using Monte Carlo simulation and smart meters measurements for the long term analysis of low voltage distribution grids with photovoltaic generation,” Int. J. Electr. Power Energy Syst. , vol. 53, pp. 468–477, 2013.  14.
V.Klonari, JF.Toubeau, Z. De Greve, O.Durieux, J.Lobry, F.Vallee, “Probabilistic Simulation Framework of Balanced and Unbalanced Low Voltage Networks,” Int. J. Electr. Power Energy Syst . (Article in Press), 2016.  15.
IEC, “IEC 60364552 Table A5210.” .  16.
F. Vallée, F. Moutier, V. Klonari, J. Toubeau, and J. Lobry, “Clustering of Photovoltaic Generation for the Consideration of Time Changing Geographical Correlation in Probabilistic Analysis of Low Voltage Distribution Systems,” in 5th Solar Integration Workshop , 2015.  17.
S. Lefebvre, “Contribution à l’utilisation d’un outil d’analyse technicoéconomique de réseaux de distribution Basse Tension en l’absence de compteurs électriques intelligents,”, UMONS, 2015.  18.
B. Meersman, B. Renders, L. Degroote, T. Vandoorn, and L. Vandevelde, “Threephase inverterconnected DGunits and voltage unbalance,” Electr. Power Syst. Res. , vol. 81, no. 4, pp. 899–906, 2011.  19.
FP7 INCREASE Project, “http://www.projectincrease.eu/.” .  20.
C.E.I. 0–21, “Reference technical rules for the connection of active and passive users to the LV electrical utilities,” 2012.  21.
M. AbdelAkher, K. M. Nor, and A. H. AbdulRashid, “Development of unbalanced threephase distribution power flow analysis using sequence and phase components,” in 2008 12th International Middle East Power System Conference, MEPCON 2008 , 2008, pp. 406–411.  22.
M. Abdel Akher, K. M. Nor, and a. H. a Rashid, “Improved threephase powerflow methods using sequence components,” IEEE Trans. Power Syst. , vol. 20, no. 3, pp. 1389–1397, 2005.  23.
I. T. Papaioannou, A. Purvins, and E. Tzimas, “Demand shifting analysis at high penetration of distributed generation in low voltage grids,” Int. J. Electr. Power Energy Syst. , vol. 44, no. 1, pp. 540–546, 2013.