Quantification and Mitigation of Unfairness in Active Power Curtailment of Rooftop Photovoltaic Systems Using Sensitivity Based Coordinated Control

: With increasing photovoltaic (PV) penetration in low voltage networks (LVNs), voltage regulation is a challenge. Active power curtailment (APC) is one possible solution for mitigating over voltages resulting from active power injection in LVNs. There is an inherent unfairness in the APC scheme. When generation is high and consumption is low, the voltages at the end of the feeder tend to be the highest. This results in high curtailment of active power output of the inverters located at the end of the feeder and low or even no curtailment for the inverts located closer to the transformer. A secondary voltage controller has been implemented to mitigate this unfairness in APC based voltage support schemes. The focus of this work is to quantify this unfairness and develop methods that enable residential PV owners serviced by the same feeder to participate equally in voltage regulation in the LVN.


Introduction
Crystalline silicon photovoltaic (PV) systems currently dominate the PV market.Even though crystalline silicon PVs have high efficiency compared to other technologies, high production costs have motivated researchers to find a more cost effective alternative.One such solution is the thin film technology, which requires lesser material to produce and is cheaper as a result.They currently however have much lower efficiencies when compared with crystalline silicon technologies [1].Advances in technology, ever increasing demand for green energy and favorable government policies have resulted in rapid increase of PV systems in a number of energy markets around the world [2].A considerable percentage of the PV systems have been installed in low voltage network (LVN).High PV penetration results in a number of engineering challenges like harmonic distortion, light flicker due to voltage fluctuation, reverse power flow, and over voltage [3][4][5].These challenges are a major limiting factor in increasing PV penetration in distribution networks.
Distribution systems have been designed to deliver power from the high voltage side of the electrical grid to end customers connected in the LVN.Current distribution systems at the time of conception were designed to cater for the voltage drop caused by consumer load and inherent impedance of the distribution lines.Injection of power at the LVN results in over voltage which the networks were not designed to cater for.Traditionally, network voltage has been regulated by the on-load tap changer (OLTC) of the medium voltage/high voltage (HV/MV) transformers [6].
Energies 2016, 9, 436 2 of 16 Distributed generation at medium and LVNs reduces the voltage regulation capability of such a scheme [7].At times when there is high generation from PVs and low load consumption, power is injected into the grid and line impedance causes the voltage at the point of common coupling (PCC) to rise [8].This may lead to violation of voltage bands prescribed by standards such as EN 50160 in Europe and ANSI C84.1 in the US [9,10].Distribution system operators (DSO) therefore require additional control over voltage regulation.
Over voltage in distribution networks is associated with a number of problems in power systems.It increases the stress on insulation of the electrical appliances connected to the grid by consumers.This can significantly reduce the life of a product.In worst case scenarios, insulation failure can occur resulting in faults.A number of electrical appliances have built-in over voltage protection that shuts down the appliance in case of a severe over voltage [11,12].Over voltage can potentially trigger protection relays resulting in unintentional islanding of a feeder from the main grid [13].
In literature, a number of schemes have been suggested to improve voltage regulation in the distribution network and thereby making it possible to increase the PV penetration.Turitsyn et al. [14] proposed reactive power compensation (RPC) by PV owners.Injection of reactive power however results in higher switching losses in an inverter.In several countries-e.g., Canada, India, and Bangladesh-RPC is not legal by rooftop PV owners and inverters are only allowed to operate at unity power factor [15][16][17].In these countries, active power curtailment (APC) is the only way of regulating voltage at the PCC.It is possible for feeders in LVNs to have R/X > 1 in such a scenario, voltage is more sensitive to active power injection rather than RPC [8].On the other hand, APC means loss of revenue for PV owners.Based on their work in [18] Long and Ochoa have concluded that installing an OLTC on distribution transformer can improve voltage regulation and increase PV hosting capacity of a LVN up to 20% depending on network topology and PV distribution along the feeders.Aggarwal et al. [19] investigated the possibility of using a distribution static synchronous compensator (D-STATCOM) in LVNs with high PV penetration.Results showed that not only does use of D-STATCOM improves voltage regulation, it also suppresses harmonic distortion [19].A coordinated control strategy for OLTC and D-STATCOM has been presented in [20].The authors have concluded that proposed scheme is capable of mitigating both over voltages and voltage imbalances via active and reactive power control respectively.Use of storage for voltage rise mitigation has been the subject of investigation in a number of papers [21,22].Coordinated control for APC and/or RPC of PVs is another solution to the voltage regulation problem proposed by a number of papers [23][24][25][26][27]. Implementation of these strategies would require communication infrastructure.
The Voltage at the PCC is a function of the active and reactive power being injected at that point and the grid impedance as seen by the inverter at PCC.As the electrical distance from a transformer increases, so does the voltage sensitivity.At times of high generation, inverters connected at PCC with higher voltage sensitivity may violate voltage limits.As a result, the inverters that are violating voltage limits have their output curtailed.Essentially, PV owners located at the end of the feeder are more susceptible to voltage violations and getting their power curtailed resulting in loss of revenue.DSOs commonly calculate a feeders' hosting capacity using the worst case over voltage scenario which might further limit the size of PV that can be installed by the house owner connect near the end of the feeder.Worst case voltage is calculated using a scenario assuming no load and all single phase PV inverter units are connected on the same phase.
A number of different solutions have been presented in literature to mitigate this inherent unfairness in droop based APC schemes for voltage control.In Japan for example, as of 2012, if power from a PV inverter is curtailed for more than 30 days a year (8%), the owner has to be compensated by the DSO [28].In [26], authors have suggested using batteries at times of low consumption and high distributed generator (DG) generation.An added advantage of using battery storage is that it minimizes fluctuations in DG output.Researchers have also suggested using demand side management and increasing self-consumption, thereby reducing power injection during the time power curtailment is an issue [29].unfairness by experimenting with different droop curves for DG owners and achieved excellent results.Perera et al. [30] developed an iterative algorithm that improves fairness while keeping voltage within limits.Zhao et al. [31] have addressed unfair APC by formulating it as part of the objective function in an optimization problem.The proposed algorithm uses load and solar irradiance forecast profiles to calculate active power dispatch values for PVs while minimizing both unfairness in curtailment and total curtailment and satisfying voltage constraints.References [8,32] have exploited the voltage sensitivity matrix to calculate potential voltage rise and reduction in active power injection required to ensure that voltage remains within limits.
The scope of this work is to quantify the unfairness in local APC schemes and investigate the possibility of increasing fairness using coordinated control.Quantification of unfairness is important as it facilitates in evaluating the effectiveness of the proposed solution and comparing different solutions proposed in literature.This work further investigates the consequences of improving fairness e.g., increasing curtailment.The proposed algorithm also aims to reduce excessive curtailment by using an S function.Structure of the paper is as follows: Section 2 explains the principle of APC and presents problem formulation; Section 3 the formulation of the key performance indices (KPIs) used for comparison has been presented; Section 4 details the proposed algorithm; Section 5 details network model and simulation setup; simulation results are presented and discussed in detail in Section 6; finally, in Section 7, conclusions have been drawn on the basis of this work.

Problem Formulation
Power flow through a line can be calculated using the Equations ( 1) and (2).Where, R + jX is the impedance of the line, u 1 and u 2 are source and line end voltage, respectively, and Φ is the phase angle between the two voltages.
In LVNs R/X ratio is typically high.High R/X ratio results in limited impact of reactive power on voltage regulation making APC a more effective option for voltage control.The phase angle Φ is generally small and sinΦ « Φ and cosΦ « 1 are reasonable approximations.Using these approximations, the change in voltage at the PCC can be formulated as Equation (3) [24].For PV inverter units operating at unity power factor Q = 0 this expression can be reduced to Equation (4).
A local droop based voltage controller has been implemented in PowerFactory (DIgSILENT, Gomaringen, Germany) [33].Below a critical voltage u cri inverter output follows maximum power point tracking (mppt).The active power is curtailed linearly once the voltage at the bus exceeds the critical voltage.If the voltage at PCC exceeds the upper voltage limit u thr the active power output of the inverter is reduced to the minimum permissible value P min .

Voltage Sensitivity Estimation
In this paper, additional logic has been implemented to estimate the voltage sensitivity of the PCC to active power (du/dp).When the Calc signal is set high, a pulse generator generates pulses for a predefined time interval.If P cal is greater than a critical value P cri , the inverter set point is switched to zero temporarily.The change in voltage is then recorded and du/dp is estimated using finite difference.The process is repeated to and the results are filtered to minimize estimation errors.For this work, P cri has been set to 70% of the rated power.P cri is intentionally set high so that a significant change in the voltage can be recorded and the impact of other inverters and the load present in the system can be neglected.The main thought process behind estimation of voltage sensitivity is to provide a truly plug-and-play solution for voltage control.Figure 1 presents a graphical overview of the PV inverter.system can be neglected.The main thought process behind estimation of voltage sensitivity is to provide a truly plug-and-play solution for voltage control.Figure 1 presents a graphical overview of the PV inverter.

Formulation of Key Performance Indices
To measure the performance of the coordinated control scheme and for the purpose of comparison with currently existing schemes, two KPIs have been formulated, namely, total energy curtailed (TEC) and curtailment unfairness index (CUI).

Total Curtailed Energy
TEC is an important KPI for comparison as any proposed coordination based control scheme should not significantly increase loss of green energy as it results in loss of revenue for the PV system owner and loss of green energy for the DSO.
where N is the number of PV systems installed on a particular feeder., , are the maximum active power output, actual active power output, and the curtailed energy of the ith inverter.TEC is therefore TEC from all inverters present in the feeder performing APC.

Unfairness in Active Power Curtailment
Fairness in APC can be examined from two perspectives.The first perspective being the loss of revenue for all PV system owners should be equal irrespective of the nominal rating of the PV system.CUI1 index formulated in Equation (9) has been used in this work to quantify the fairness from this perspective and compare different control strategies.
where is the average energy curtailed during a 24-h period from all PV systems connected to a feeder.The calculated index CUI1 is the standard deviation of curtailed energies, hence, a measure of spread of TEC during the day from the installed PV systems.If equal energy is curtailed from each PV system, CUI1 will be zero implying fair curtailment.If however the installed PV systems differ in rating a PV system with a smaller rating would incur higher percentage reduction in revenue.

Formulation of Key Performance Indices
To measure the performance of the coordinated control scheme and for the purpose of comparison with currently existing schemes, two KPIs have been formulated, namely, total energy curtailed (TEC) and curtailment unfairness index (CUI).

Total Curtailed Energy
TEC is an important KPI for comparison as any proposed coordination based control scheme should not significantly increase loss of green energy as it results in loss of revenue for the PV system owner and loss of green energy for the DSO.
TEC " where N is the number of PV systems installed on a particular feeder.P mppt i , P inv i , E cur i are the maximum active power output, actual active power output, and the curtailed energy of the ith inverter.TEC is therefore TEC from all inverters present in the feeder performing APC.

Unfairness in Active Power Curtailment
Fairness in APC can be examined from two perspectives.The first perspective being the loss of revenue for all PV system owners should be equal irrespective of the nominal rating of the PV system.
Energies 2016, 9, 436 5 of 16 CUI1 index formulated in Equation ( 9) has been used in this work to quantify the fairness from this perspective and compare different control strategies.
where E cur is the average energy curtailed during a 24-h period from all PV systems connected to a feeder.The calculated index CUI1 is the standard deviation of curtailed energies, hence, a measure of spread of TEC during the day from the installed PV systems.If equal energy is curtailed from each PV system, CUI1 will be zero implying fair curtailment.If however the installed PV systems differ in rating a PV system with a smaller rating would incur higher percentage reduction in revenue.Another way to accommodate the unfairness is to ensure that percentage reduction in revenue for every PV system owner is the same.Curtailed power can be normalized by dividing instantaneous curtailed power P cur i by the rated output of the inverter P rated i .Integrating normalized curtailed power gives normalized curtailed energy and as P rated i is a constant it can be calculated using Equation (10).
where P rated i is the rated power for the ith inverter and e cur is average of the normalized curtailed energies of all PV systems present in a feeder.The second index formulated in this work (CUI2) is a measure of spread of percentage of the energy curtailed from each PV rather than TEC from the PVs.Low values of CUI1 and CUI2 indicate small spread hence higher fairness and vice versa.

Proposed Algorithm
In this paper, a coordinated control scheme has been proposed in which a feeder in a low voltage susceptible to over voltages has a coordinating controller for voltage regulation.This controller ensures that every PV system owner participates equally in voltage regulation of a feeder.This work assumes a bidirectional ideal communication channel between the coordinating controller and smart inverters.The controller has two modes of operation.In normal operation, when maximum voltage at every point in feeder is less than u lb , no energy is curtailed.In distressed mode, the controller generates set points for every inverter.This section details the working of the proposed controller.

Calculation of Sensitivity Matrix
The impact of variation of active power injected by the inverters can be measured quantitatively using the sensitivity matrix [34].In this work, inverters have been modeled to operate at unity power factor, consequently the impact of reactive power on voltage can be neglected.
Voltage sensitivity of a node to active power injection at any node in a radial distribution network can be calculated using Equation (13) [34,35].
Energies 2016, 9, 436 6 of 16 The expressions suggest that longer the feeder length more susceptible it is to over voltages.The resistance matrix can be expanded and rewritten as a function of line length and resistance per kilometer.
For i = j, L ij is the minimum electrical distance from the transformer to the ith node.In case i ‰ j, L ij is the maximum overlap of the paths formed from the transformer to the ith and jth node.r ij is the average resistance per kilometer of the branches belonging to L ij .For i > j or j > i, if overlapping path length is much greater than non-overlapping path lengths the off-diagonal elements of the sensitivity matrix can be approximated as the diagonal elements of the matrix.
As each inverter is capable of calculating the voltage sensitivity of the connected node to active power (diagonal elements of the matrix), rest of the matrix can be populated.

Active Power Reduction Calculation
At initial run, the coordinating controller requests and receives rated powers, estimated sensitivities and node ID from every inverter connected to the feeder.Rated output power of the PV systems are stored in the matrix RP kˆl , where k is the number of unique nodes in a feeder with a PV connection and l is the maximum number of PVs connected to a single node.The sensitivities arranged in an ascending order form the diagonal part of the sensitivity matrix S kˆk .Non-diagonal elements of the sensitivity matrix are populated using Equation (15).Voltage at the secondary side of the MV/LV (low voltage) transformer u trlv is calculated using the current tap position TP, percentage voltage change per tap u tap , and nominal voltage u nom .u trlv " u nom ´utap ˆTP (16) Maximum voltage rise at all nodes connected to PV systems is calculated using the sensitivity matrix S ij and the rated powers P rated i .
The voltages measured at PCC of each inverter are sent to the coordinating controller and continuously monitored.Critical voltage for each inverter is set using Equations ( 18) and (19), where u max is the maximum of the measured voltages.
u max " max `u1,2,...,k ˘(18) If i is the index of the node with the maximum measured voltage, reduction in voltage required (∆u) to ensure no over voltage occurs is calculated using Equation (20).
Energies 2016, 9, 436 Reduction required in active power injection for each inverter is calculated using either Equation ( 21) or (22), Equation ( 23) depending on which fairness perspective is currently under investigation.

Equal Loss of Revenue for Each Photovoltaic Owner
∆P max " ∆u where m i is number of non-zero elements in the ith row of RP matrix (number of PV connected to a node).It should be noted that ∆P max is a single value hence the same for every inverter irrespective of its power rating.

Equal Percentage Reduction in Revenue
z " ∆u where z is the percentage reduction in active power for each inverter.In this case, ∆P max is k ˆl matrix.
In the final step, S function is used to reduce excessive curtailment. where: where A is a constant that controls the rate at which ∆P ij increases from zero to ∆P max ij .This ensures that maximum possible curtailment occurs only when necessary.Figure 2 presents graphic overview of the proposed scheme.
where z is the percentage reduction in active power for each inverter.In this case, ∆P max is k × l matrix.
In the final step, S function is used to reduce excessive curtailment. where: where A is a constant that controls the rate at which Δ increases from zero to Δ .This ensures that maximum possible curtailment occurs only when necessary.Figure 2 presents graphic overview of the proposed scheme.

Network Model
Unlike industrial or commercial feeders, active power demand in residential feeders typically peaks at night which does not coincide with peak generation, which is usually mid-day.This is why

Network Model
Unlike industrial or commercial feeders, active power demand in residential feeders typically peaks at night which does not coincide with peak generation, which is usually mid-day.This is why residential feeders are more prone to over voltages.
The test case chosen for this work is the Cigre low voltage distribution network [36].It consists of a three feeders network connected to medium voltage network through a 400 KVA transformer.The cumulative residential load connected to a node is represented by a lumped load connected to the feeder shown in Figure 2. In this study, each of the five lumped loads connected to the residential feeder have a PV system installed.The PV systems have been randomly sized between 25 kVA and 50 kVA.PV inverters provide voltage support using APC.Table 1 lists the rated power of the inverters.For the purpose of simulation, one hour average profiles have been used for the loads and one minute average profiles have been used for the PV systems.Normalized load and PV profiles used for the study are shown in Figure 3.For the purpose of simulation, one hour average profiles have been used for the loads and one minute average profiles have been used for the PV systems.Normalized load and PV profiles used for the study are shown in Figure 3.

PowerFactory Python Interfacing
The distribution network and the local voltage controller detailed in Section 2.1 have been implemented in PowerFactory [33], a popular power network simulation tool.The secondary voltage controller, detailed in Section 4, has been implemented in Python [37] because of PowerFactory's limitations while dealing with matrices.During RMS simulation it is possible for controllers modelled in DSL (DigSilent programming language) blocks to call functions from an external C ++ library.A C ++ library has been used to implement sockets to communicate with Python.In this work, socket based communication approach has been preferred as it can easily be extended to incorporate a network simulator such as OMNET ++ to study the impact of communication on voltage control for future work.Figure 4 is a graphical illustration of coupling scheme used.

PowerFactory Python Interfacing
The distribution network and the local voltage controller detailed in Section 2 have been implemented in PowerFactory [33], a popular power network simulation tool.The secondary voltage controller, detailed in Section 4, has been implemented in Python [37] because of PowerFactory's limitations while dealing with matrices.During RMS simulation it is possible for controllers modelled in DSL (DigSilent programming language) blocks to call functions from an external C ++ library.A C ++ library has been used to implement sockets to communicate with Python.In this work, socket based communication approach has been preferred as it can easily be extended to incorporate a network simulator such as OMNET ++ to study the impact of communication on voltage control for future work.Figure 4 is a graphical illustration of coupling scheme used.
Energies 2016, 9, 436 9 of 16 limitations while dealing with matrices.During RMS simulation it is possible for controllers modelled in DSL (DigSilent programming language) blocks to call functions from an external C ++ library.A C ++ library has been used to implement sockets to communicate with Python.In this work, socket based communication approach has been preferred as it can easily be extended to incorporate a network simulator such as OMNET ++ to study the impact of communication on voltage control for future work.Figure 4 is a graphical illustration of coupling scheme used.

Results and Discussion
In the first step, local voltage sensitivities (diagonal values of the sensitivity matrix ) to active power are calculated by the PV inverter using the logic presented in Figure 1.Once Calc signal is set high, inverter momentarily disconnects every 20 min 25 times.Change in output power and the voltage at the PCC is used to estimate the local sensitivity.The coordinating controller ensures no two inverters calculate the sensitivity at the same time.Figure 5 shows inverter output, voltage, and du/dp estimation results for inverter PV-Res03.

Results and Discussion
In the first step, local voltage sensitivities (diagonal values of the sensitivity matrix S ij ) to active power are calculated by the PV inverter using the logic presented in Figure 1.Once Calc signal is set high, inverter momentarily disconnects every 20 min 25 times.Change in output power and the voltage at the PCC is used to estimate the local sensitivity.The coordinating controller ensures no two inverters calculate the sensitivity at the same time.Figure 5 shows inverter output, voltage, and du/dp estimation results for inverter PV-Res03.The base case results presented in Figure 7 are for PV systems without local and coordinating voltage controller.For a scenario with high PV generation and low consumption, the voltages measured at the residential feeder exceed the permissible voltage band and hence require voltage controllers for better regulation.The base case results presented in Figure 7 are for PV systems without local and coordinating voltage controller.For a scenario with high PV generation and low consumption, the voltages measured at the residential feeder exceed the permissible voltage band and hence require voltage controllers for better regulation.The base case results presented in Figure 7 are for PV systems without local and coordinating voltage controller.For a scenario with high PV generation and low consumption, the voltages measured at the residential feeder exceed the permissible voltage band and hence require voltage controllers for better regulation.In this paper, two design approaches have been implemented for voltage regulation.Each approach is discussed in detail in the subsequent subsections.In this paper, two design approaches have been implemented for voltage regulation.Each approach is discussed in detail in the subsequent subsections.

Local Voltage Regulation
In this approach, each inverter is capable of voltage regulation using APC.This approach requires only local information hence no communication is required.For the experiment, it has been assumed that maximum permissible reduction in active power output of the inverter (P min ) is 50% of the rated power.The critical (u cri ) and threshold voltages (u thr ) have been set at 1.02 p.u. and 1.05 p.u. respectively for every inverter on the residential feeder.Voltage profiles for the inverters are presented in Figure 8a.

Local Voltage Regulation
In this approach, each inverter is capable of voltage regulation using APC.This approach requires only local information hence no communication is required.For the experiment, it has been assumed that maximum permissible reduction in active power output of the inverter (Pmin) is 50% of the rated power.The critical (u cri ) and threshold voltages (u thr ) have been set at 1.02 p.u. and 1.05 p.u. respectively for every inverter on the residential feeder.Voltage profiles for the inverters are presented in Figure 8a.

Coordinated Voltage Regulation
In this work two coordinating controller models have been implemented to minimize the two perspectives of unfairness presented in Section 3.2.Coordinating control schemes implemented in this paper require information from every PV inverter installed on the feeder, which requires a communication network.In this paper, an ideal communication channel has been assumed.For simulation purposes, the upper u ub and lower u lb limits of the voltage regulation band for the coordinating controller have been set at 1.02 and 1.05, respectively.

Equal Loss of Revenue for Each Photovoltaic Owner
The first controller (CC1) generates set points for critical voltage critical u cri in p.u. and the reduction in active power ∆P in kW. Figure 9a shows the plot of the critical voltages calculated for

Coordinated Voltage Regulation
In this work two coordinating controller models have been implemented to minimize the two perspectives of unfairness presented in Section 3.2.Coordinating control schemes implemented in this paper require information from every PV inverter installed on the feeder, which requires a communication network.In this paper, an ideal communication channel has been assumed.
For simulation purposes, the upper u ub and lower u lb limits of the voltage regulation band for the coordinating controller have been set at 1.02 and 1.05, respectively.

Equal Loss of Revenue for Each Photovoltaic Owner
The first controller (CC1) generates set points for critical voltage critical u cri in p.u. and the reduction in active power ∆P in kW. Figure 9a shows the plot of the critical voltages calculated for each PV inverter.u cri is calculated only when u max is greater than u lb , otherwise, u cri is set to 1.05 which means no APC.Figure 9b is the plot for maximum expected over voltage (∆u).It is import to note that even though u max is greater than u lb between the time 2 a.m. and 7 a.m., maximum expected voltage rise is ∆u which is less than zero for the duration, hence no power is curtailed during the above mentioned period.

Coordinated Voltage Regulation
In this work two coordinating controller models have been implemented to minimize the two perspectives of unfairness presented in Section 3.2.Coordinating control schemes implemented in this paper require information from every PV inverter installed on the feeder, which requires a communication network.In this paper, an ideal communication channel has been assumed.For simulation purposes, the upper u ub and lower u lb limits of the voltage regulation band for the coordinating controller have been set at 1.02 and 1.05, respectively.

Equal Loss of Revenue for Each Photovoltaic Owner
The first controller (CC1) generates set points for critical voltage critical u cri in p.u. and the reduction in active power ∆P in kW. Figure 9a shows the plot of the critical voltages calculated for each PV inverter.u cri is calculated only when u max is greater than u lb , otherwise, u cri is set to 1.05 which means no APC.Figure 9b is the plot for maximum expected over voltage (∆u).It is import to note that even though u max is greater than u lb between the time 2 a.m. and 7 a.m., maximum expected voltage rise is ∆u which is less than zero for the duration, hence no power is curtailed during the above mentioned period.The power curtailment signal (∆P) for each PV inverter is the function of the maximum permissible curtailment (∆P max ), maximum voltage in the network (u max ) and the upper and lower The power curtailment signal (∆P) for each PV inverter is the function of the maximum permissible (∆P max ), maximum voltage in the network (u max ) and the upper and lower voltage limits.Figure 10 shows the impact of the S function used in Equations ( 24) and ( 25) on APC.As u max approaches u ub , ∆P approaches ∆P max .This improves efficiency by minimizing unnecessary APC.It is important to note that variable A in Equation ( 25) can be used to adjust the slope of the S function.In this study, the value of A is 5 for all the experiments.The value of ∆P calculated by CC1 is the same for all PV inverters which ensures equal curtailment for every PV system.The use of S function reduces curtailment set point for PV installed at Res01 and reduces E cur from 221 kWh (theoretical maximum) to 148 kWh a reduction of 32.6% (the grey shaded region in Figure 10).voltage limits.Figure 10 shows the impact of the S function used in Equations ( 24) and ( 25) on APC.As u max approaches u ub , ∆P approaches ∆P max .This improves efficiency by minimizing unnecessary APC.It is important to note that variable A in Equation ( 25) can be used to adjust the slope of the S function.In this study, the value of A is 5 for all the experiments.The value of ∆P calculated by CC1 is the same for all PV inverters which ensures equal curtailment for every PV system.The use of S function reduces curtailment set point for PV installed at Res01 and reduces E cur from 221 kWh (theoretical maximum) to 148 kWh a reduction of 32.6% (the grey shaded region in Figure 10).

Equal Percentage Reduction in Revenue
The second controller (CC2) like the controller implemented for equal revenue loss also generates set points for critical voltage critical u cri in p.u. and the reduction in active power Δ in kW.In this case however, ∆P is matrix and the set point for each inverter is a function of its rated power.Figure 12a shows critical voltages calculated for each PV inverter.Figure 12b shows maximum reduction in active power calculated for each inverter to ensure voltage remains within bounds.

Equal Percentage Reduction in Revenue
The second controller (CC2) like the controller implemented for equal revenue loss also generates set points for critical voltage critical u cri in p.u. and the reduction in active power ∆P in kW.In this case however, ∆P is matrix and the set point for each inverter is a function of its rated power.Figure 12a shows critical voltages calculated for each PV inverter.Figure 12b shows maximum reduction in active power calculated for each inverter to ensure voltage remains within bounds.

Equal Percentage Reduction in Revenue
The second controller (CC2) like the controller implemented for equal revenue loss also generates set points for critical voltage critical u cri in p.u. and the reduction in active power Δ in kW.In this case however, ∆P is matrix and the set point for each inverter is a function of its rated power.Figure 12a shows critical voltages calculated for each PV inverter.Figure 12b shows maximum reduction in active power calculated for each inverter to ensure voltage remains within bounds.For PV system installed at Res05 the TEC during the day is 193.10 kWh, 34.2% less than maximum possible curtailment.Figure 13 shows reduction in curtailment signal using S function.Figure 14 shows feeder voltage profiles and power curtailed from each PV inverter.Comparing Figure 8b, Figure 11b, and Figure 14b, impact of different control schemes on curtailed power is visible.
Energies 2016, 9, 436 12 of 15 For PV system installed at Res05 the TEC during the day is 193.10 kWh, 34.2% less than maximum possible curtailment.Figure 13 shows reduction in curtailment signal using S function.Figure 14 shows feeder voltage profiles and power curtailed from each PV inverter.Comparing Figures 8b, 11b, and 14b, impact of different control schemes on curtailed power is visible.Figure 15 shows the energy curtailed from each PV owner for the three control schemes simulated for this work.Voltage regulation schemes using local APC are unfair to PV system owner connected near the end of the feeder.The energy curtailed from PV at Res01 is only 10.64% of the TEC from all inverters.PV at Res05 on the other hand gets 245.7 kWh curtailed which is over 40% of the total curtailed energy.It, therefore, contributes over four times more that PV at Res01 in the voltage regulation process.
In the second scenario (CC1), energy curtailed from each PV ranges between 18.2% and 22.4%.In this scenario, every PV owner contributes approximately equally in the voltage regulation process.One drawback of this scheme however, is that those PV owners with smaller ratings incur a bigger percentage decrease in revenue generated by the PV.At times of peak curtailment PV-Res05 rated at 50 kW gets 20% of the output curtailed, but PV-Res03 rated at 27 kW loses 38% of the revenue generated by the PV at peak output.Figure 15 shows the energy curtailed from each PV owner for the three control schemes simulated for this work.Voltage regulation schemes using local APC are unfair to PV system owner connected near the end of the feeder.The energy curtailed from PV at Res01 is only 10.64% of the TEC from all inverters.PV at Res05 on the other hand gets 245.7 kWh curtailed which is over 40% of the total curtailed energy.It, therefore, contributes over four times more that PV at Res01 in the voltage regulation process.Figure 15 shows the energy curtailed from each PV owner for the three control schemes simulated for this work.Voltage regulation schemes using local APC are unfair to PV system owner connected near the end of the feeder.The energy curtailed from PV at Res01 is only 10.64% of the TEC from all inverters.PV at Res05 on the other hand gets 245.7 kWh curtailed which is over 40% of the total curtailed energy.It, therefore, contributes over four times more that PV at Res01 in the voltage regulation process.
In the second scenario (CC1), energy curtailed from each PV ranges between 18.2% and 22.4%.In this scenario, every PV owner contributes approximately equally in the voltage regulation process.One drawback of this scheme however, is that those PV owners with smaller ratings incur a bigger percentage decrease in revenue generated by the PV.At times of peak curtailment PV-Res05 rated at 50 kW gets 20% of the output curtailed, but PV-Res03 rated at 27 kW loses 38% of the revenue generated by the PV at peak output.In the second scenario (CC1), energy curtailed from each PV ranges between 18.2% and 22.4%.In this scenario, every PV owner contributes approximately equally in the voltage regulation process.One drawback of this scheme however, is that those PV owners with smaller ratings incur a bigger percentage decrease in revenue generated by the PV.At times of peak curtailment PV-Res05 rated at 50 kW gets 20% of the output curtailed, but PV-Res03 rated at 27 kW loses 38% of the revenue generated by the PV at peak output.
In the final scenario (CC2), energy curtailed from each PV is also a function of the rating of the PV system.Energy curtailed from each PV ranges between 13.1% and 26.5% of total curtailed energy.In this study case, All PV system owners do not participate equally in voltage regulation.However, every PV owner on the feeder loses the same percentage of revenue.
To evaluate the performance of control schemes implemented in this paper, KPI formulated in Section 3 have been calculated and listed in Table 2. Loss of green energy increases for control schemes CC1 and CC2 by 18.2% and 14.3%, respectively.CUI1 reduces by 85% for CC1 and by 57% for CC2.CC1 has negligible impact on the second unfairness index while CC2 reduces CUI2 to almost zero meaning APC for every PV owner is the same.

Conclusions
With increasing PV penetration in LVNs, voltage regulation is a challenge.APC is one possible solution for over voltages.PV inverters installed at feeder end are most susceptible to over voltages and are most likely to get their output power curtailed.Local APC schemes are therefore inherently unfair.In this paper, a truly plug-and-play coordinated control scheme has been implemented with two variations.Two perspectives of unfairness have been presented in the paper and two indices have been formulated to quantify the unfairness in APC schemes.The effectiveness of the control schemes has been tested by implementing network in PowerFactory and the controllers in Python.Both variants of the coordinating controllers were able to increase fairness in APC schemes significantly.As proof of concept, the proposed algorithms have been implemented on the Cigre low voltage distribution network.Results show that it is possible to improve fairness in APC but this comes at the cost of increased curtailment for PV owners.In certain scenarios, a significant increase in curtailment may make the proposed methods economically infeasible.Future work will be the implementation of the proposed algorithms on a number of LV networks to better understand the effectiveness and the limitations of the proposed schemes.

Figure 4 .
Figure 4. PowerFactory and Python coupling via sockets using external dynamic link library.

Figure 4 .
Figure 4. PowerFactory and Python coupling via sockets using external dynamic link library.

Figure 4 .
Figure 4. PowerFactory and Python coupling via sockets using external dynamic link library.

Figure 6
Figure 6 presents a comparison of the du/dp calculated at PV connected nodes by PowerFactory and estimated by PV inverters.The average error between the actual and the estimated value of voltage sensitivity to active power is 8.47%.

Figure 6 .
Figure 6.(a) A comparison of actual and estimated du/dp values; and (b) estimation error.

Figure 6 15 Figure 5 .
Figure 6 presents a comparison of the du/dp calculated at PV connected nodes by PowerFactory and estimated by PV inverters.The average error between the actual and the estimated value of voltage sensitivity to active power is 8.47%.

Figure 6
Figure 6 presents a comparison of the du/dp calculated at PV connected nodes by PowerFactory and estimated by PV inverters.The average error between the actual and the estimated value of voltage sensitivity to active power is 8.47%.

Figure 6 .
Figure 6.(a) A comparison of actual and estimated du/dp values; and (b) estimation error.

Figure 6 .
Figure 6.(a) A comparison of actual and estimated du/dp values; and (b) estimation error.

Figure 6 .
Figure 6.(a) A comparison of actual and estimated du/dp values; and (b) estimation error.

Figure 7 .
Figure 7. Voltage profiles at PCC of each inverter with no voltage control function.

Figure 7 .
Figure 7. Voltage profiles at PCC of each inverter with no voltage control function.

Figure 8 .
Figure 8.(a) Voltage profiles at PCC of each inverter with local voltage controller; and (b) curtailed power from each inverter.

Figure
Figure 8b shows the TEC from each PV inverter.The figure shows the inherent unfairness in local APC schemes.PV inverters connected close to the transformer get significantly less energy curtailed in comparison to the ones located near the end of the feeder.

Figure 8 .
Figure 8.(a) Voltage profiles at PCC of each inverter with local voltage controller; and (b) curtailed power from each inverter.

Figure
Figure 8b shows the TEC from each PV inverter.The figure shows the inherent unfairness in local APC schemes.PV inverters connected close to the transformer get significantly less energy curtailed in comparison to the ones located near the end of the feeder.

Figure 9 .
Figure 9. (a) Critical voltages calculated for each PV; and (b) maximum expected over voltage.

Figure 9 .
Figure 9. (a) Critical voltages calculated for each PV; and (b) maximum expected over voltage.

Figure 10 .
Figure 10.Impact of using S function on the curtailment signal.

Figure 11a presents the
Figure 11a presents the voltage profiles at PCC and curtailed energy from each PV inverter during the day.The improvement in fairness in APC is evident from Figure 11b.Every PV inverter participates equally in voltage regulation irrespective of their individual ratings.

Figure 10 .
Figure 10.Impact of using S function on the curtailment signal.

Figure
Figure 11a presents the voltage profiles at PCC and curtailed energy from each PV inverter during the day.The improvement in fairness in APC is evident from Figure 11b.Every PV inverter participates equally in voltage regulation irrespective of their individual ratings.

Figure 10 .
Figure 10.Impact of using S function on the curtailment signal.

Figure
Figure 11a presents the voltage profiles at PCC and curtailed energy from each PV inverter during the day.The improvement in fairness in APC is evident from Figure 11b.Every PV inverter participates equally in voltage regulation irrespective of their individual ratings.

Figure 11 .
Figure 11.(a) Voltage profiles at PCC of each inverter with CC1; and (b) curtailed energy from each inverter.

Figure 12 .
Figure 12.(a) Critical voltages calculated for each PV and (b) maximum curtailment for each inverter.

Figure 11 .
Figure 11.(a) Voltage profiles at PCC of each inverter with CC1; and (b) curtailed energy from each inverter.

Figure 10 .
Figure 10.Impact of using S function on the curtailment signal.

Figure
Figure 11a presents the voltage profiles at PCC and curtailed energy from each PV inverter during the day.The improvement in fairness in APC is evident from Figure 11b.Every PV inverter participates equally in voltage regulation irrespective of their individual ratings.

Figure 11 .
Figure 11.(a) Voltage profiles at PCC of each inverter with CC1; and (b) curtailed energy from each inverter.

Figure 12 .
Figure 12.(a) Critical voltages calculated for each PV and (b) maximum curtailment for each inverter.Figure 12. (a) Critical voltages calculated for each PV and (b) maximum curtailment for each inverter.

Figure 12 .
Figure 12.(a) Critical voltages calculated for each PV and (b) maximum curtailment for each inverter.Figure 12. (a) Critical voltages calculated for each PV and (b) maximum curtailment for each inverter.

Figure 13 .
Figure 13.Impact of S function the on curtailment sent to the Inverter at Res05.Figure 13.Impact of S function the on curtailment sent to the Inverter at Res05.

Figure 13 .
Figure 13.Impact of S function the on curtailment sent to the Inverter at Res05.Figure 13.Impact of S function the on curtailment sent to the Inverter at Res05.

Figure 13 .
Figure 13.Impact of S function the on curtailment sent to the Inverter at Res05.

Figure 14 .
Figure 14.(a) Voltage profiles at PCC of each inverter with CC2 and (b) curtailed energy from each inverter.

Figure 15 .
Figure 15.Energy curtailed from each inverter for the control schemes implemented in the paper.

Figure 14 .
Figure 14.(a) Voltage profiles at PCC of each inverter with CC2 and (b) curtailed energy from each inverter.

Figure 13 .
Figure 13.Impact of S function the on curtailment sent to the Inverter at Res05.

Figure 14 .
Figure 14.(a) Voltage profiles at PCC of each inverter with CC2 and (b) curtailed energy from each inverter.

Figure 15 .
Figure 15.Energy curtailed from each inverter for the control schemes implemented in the paper.Figure 15.Energy curtailed from each inverter for the control schemes implemented in the paper.

Figure 15 .
Figure 15.Energy curtailed from each inverter for the control schemes implemented in the paper.Figure 15.Energy curtailed from each inverter for the control schemes implemented in the paper.
Tonkoski et al. [8] developed a coordinated algorithm to mitigate

Table 2 .
Calculated key performance indices (KPI) for the implemented control schemes.CUI: curtailment unfairness index.