Abstract

Distributed rooftop photovoltaic (PV) generators prospered distributed generation (DG) in recent years. Certain randomness of rooftop PV connection may lead to significant PV power imbalance across three phases, especially in low-voltage distribution systems. Due to interphase line coupling, traditional Var compensation methods which typically have competent voltage regulation performance may become less effective in such PV imbalance scenarios. In this paper, the limitation of traditional Var compensation methods in voltage regulation with unbalanced PV power integration is demonstrated and comprehensively analyzed. After describing the voltage regulation challenge, based on the voltage sensitivity analysis, it is revealed that PV power unbalanced level together with equivalent mutual impedance among phase conductors has a significant impact on the effectiveness of traditional Var compensation methods on voltage regulation. On this basis, to improve the performance of voltage regulation methods, some suggestions are proposed for both current system operation and future distribution system planning. Numerical studies demonstrate the effectiveness of the proposed suggestions. Future rooftop PV integration in LV systems can benefit from this research.

1. Introduction

In recent years, various distributed generation and storage systems including photovoltaic, wind power, electric vehicles, etc., have been developed vigorously [16]. As one of the most attractive options of distributed generation (DG) [710], more and more single-phase rooftop photovoltaic (PV) generators are integrated into low-voltage distribution systems. These small-size PV systems installed on rooftops of customers’ houses as distributed generators can support household appliances and feed excessive PV power back to the grid. However, as PV penetration increases, consequent reverse power flow might induce significant voltage rise at the end of distribution feeders [11, 12].

To avoid the potential overvoltage problem, the maximum PV penetration of a distribution system should be carefully assessed. In early research, single-phase equivalent systems are used to estimate possible voltage problems with high PV penetration levels [13, 14]. While in [15], a stochastic method is designed to imitate the random connection of single-phase PV generators and to estimate the corresponding three-phase voltage problem in low-voltage systems.

In order to accommodate more rooftop PV generators, Var compensation devices are required to actively participate in voltage regulation [1619]. Optimal voltage/Var control methods of inverters have been discussed by many previous publications [2025]. However, distribution network optimization requires full or partial observability of the entire network through communication systems and information exchange systems, which are not available for most low-voltage distribution systems nowadays. Furthermore, since optimal methods take time, its respond speed might not be fast enough to follow the variation in sunlight intensity. To simplify the optimal control methods and shorten computing time, some multilevel dispatch is proposed. In [26], zonal voltage control combining day-ahead dispatch and real-time control for distribution networks with high proportion of PV power is researched. In [27], a bilevel voltage/Var optimization to coordinate smart inverters with voltage control devices is carried out. However, a PV system may lose 70% of its generation in a short period due to fast-moving cloud coverage [28], which makes centralized optimal control methods difficult to keep up with fast fluctuating PV power. Therefore, optimal methods to coordinate all distributed PV inverters in distribution network are still in the academic discussion stage.

On the other hand, Var compensation methods that only rely on local measurements (e.g., PV generation and local line-to-ground voltage) also have already been widely implemented in centralized PV plants for the point of common coupling (PCC) voltage regulation [2931]. The main advantage of these methods is their rapid response speed, which makes them suitable to counter fast voltage fluctuations caused by PV power variations and with more attractive options for distributed rooftop PV generators. In this field, previous research mainly focuses on the design of Q(P) and Q(V) curves in order to improve their voltage regulation performance [32, 33].

However, since single-phase rooftop PV generators are integrated into low-voltage distribution systems randomly, PV power penetration across three phases tends to be unbalanced. In such situations, widely implemented locally dependent Var compensation methods that can successfully control the PCC voltage of large-scale PV plants may become less effective in voltage regulation with distributed rooftop PV generators. In this paper, the limitation of locally dependent Var compensation methods in voltage regulation with unbalanced PV power integration is comprehensively analyzed.

The remainder of this paper is organized as follows: Section 2 introduces background of proposed problem. In Section 3, the limitation of locally dependent Var compensation methods in voltage regulation with unbalanced PV power integration is revealed. Section 4 analyzed this voltage regulation problem. Some suggestions are proposed for both current system operation and future distribution system planning in Section 5. Numerical studies are presented in Section 6. Finally, Section 7 concludes the paper.

2. Background Introduction

2.1. Rooftop PV Generators and Reverse Power Flow

Rooftop PV generators are becoming more and more popular in recent days, not only due to its clean and renewable characteristics, but these PV systems can sell excessive power back to the utility after providing power supply for customers’ household appliances.

PV panels generate active power only in the daytime, with a peak value during noon as shown in Figure 1(a). While, for typical residential load profiles such as in Figure 1(b), peak load usually occurs in the morning and evening, with valley load during noon and night. The valley load during noon can be less than 30% of its peak value. Therefore, PV power would exceed load demand during noon, and reverse power flow would be caused (by negative load) in distribution feeders as shown in Figure 1(c).

2.2. PV System Var Generation Control

Schematic diagram of a PV system with corresponding control loops is displayed in Figure 2 [34]. The main control task for a PV system is to regulate its DC-link voltage following a maximum power-point tracking (MPPT) scheme which is designed to maximize the electricity power harvested from PV panels. At the same time, PV inverters may also be required to generate a certain amount of reactive power for system voltage regulation. Specifically, the control tasks as discussed above can be realized following the three coordinated parts.

2.2.1. Phase-Locked loop (PLL)

A PLL component is adopted to synchronize PWM and control schemes to the PCC voltage. In this way, the AC signals are converted into dq-frame correspondence signals, and the controllers can deal with their DC equivalent values instead of the original sinusoidal signals.

2.2.2. DC-Link Voltage-Control Loop

The error between the square of the DC voltage and its corresponding reference value is dealt with by using a feed-forward compensator, which neutralizes the instability and nonlinearity of PV panels and improves the stability of the PV system. is usually obtained from MPPT schemes, with a certain variation range in order to ensure safe operation of inverters. The output of the compensation is enhanced by the feed-forward signal, and the current command is generated.

2.2.3. Current-Control Loop

The current commands are transmitted to a current-control scheme, which is established in a dq-frame drive to track and dq-frame drive to track . It is worth noting that the active power output is controlled by , while is responsible for adjusting the reactive power generation of a PV inverter. Its value of Var generation can be determined by specific Var compensation strategies, which will be introduced in the next section. In addition, the current-control scheme also protects the PV inverter from overload and external faults, as long as the value of is properly limited by saturation blocks.

2.3. Locally Dependent Var Compensation Methods

Reverse power flow caused by a large amount of PV power integration will significantly increase system voltage at the end of distribution feeders, which might induce overvoltage issues. In order to mitigate the voltage rise, PV inverters are required to provide Var compensation.

Due to the effectiveness and easy implementation, Var compensation methods that only rely on local measurements have already been widely applied on voltage regulation for large-scale PV plants so far. Two typical locally dependent Var compensation methods are shown in Figure 3. On the one hand, in Figure 3(a), PV inverters are required to operate with a constant power factor (e.g., 0.9 leading). Namely, the reactive power generation of a PV inverter should always be in proportion to its active power output. In application, local PV active power generation at time instant is measured in real time, and the reactive power generation can be obtained according to and the constant power factor . On the other hand, Figure 3(b) is a power factor droop curve, through which the power factor of a PV inverter varies according to its local voltage. Specifically, if the local voltage is too high, PV inverters will operate with a leading power factor to absorb reactive power from the grid. Conversely, PV inverters will operate with a lagging power factor to inject reactive power into the grid. In real life application, both local voltage and PV generation will be measured in real time. No matter which method is used, the power factor of grid-connected PV inverters is generally restricted within the range from 0.9 lagging to 0.9 leading by the utility [35].

3. Problem Description

3.1. Voltage Regulation in Traditional Distribution Systems

Residential customers are almost directly connected to 415 V low-voltage systems which are fed by 11 kV/415 V distribution transformers. Since the 11 kV/415 V transformer cannot participate in system voltage regulation due to its fixed tap position, voltage fluctuations in the 11 kV side have a significant impact on 415 V systems. A month-long voltage profile recorded from the secondary side of an 11 kV/415 V distribution transformer is displayed in Figure 4, which mainly fluctuates within the range between 1.0 pu and 1.05 pu. Therefore, 415 V systems were designed for safe operation with its all possible load levels as long as the secondary side voltage of the 11 kV/415 V transformer fluctuates within the normal range.

3.2. Voltage Regulation Problem

In traditional distribution systems, residential loads are approximately balanced across three phases. However, rooftop PV generators are usually installed randomly. Therefore, PV power integration tends to be unbalanced especially in low-voltage systems with limited customers. A possible scenario (PV penetration is 25%, 30%, and 45% in Phases A to C, respectively) is assumed in Case 1. Detailed load and PV capacity are shown in Table 1. In such a situation, two locally dependent Var compensation methods shown in Figure 3 are implemented on PV inverters, respectively, in order to test their voltage regulation performance with unbalanced PV power integration. As shown in Table 2, in a low load but high PV generation scenario, voltage at Phase B is much higher than that of the other two phases and neither the constant power factor method nor the power factor droop curve method can eliminate this overvoltage problem (>1.06 [16]) on Phase B.

Nevertheless, there exists an interesting phenomenon that with the PV penetration of all three phases increased to 50% as in Case 2, the overvoltage problem on Phase B could be successfully eliminated by either of those two locally dependent Var compensation methods. Detailed case data and voltage regulation performance are, respectively, listed in Tables 1 and 2.

Generally, since a large amount of PV power integration might cause overvoltage problem, the maximum allowable PV penetration level should be estimated by utilities. However, based on voltage regulation performance in Case 1 and Case 2, if a higher but balanced PV penetration is regarded as the most severe case, the potential overvoltage issue might be under estimated. An overvoltage problem might occur before the estimated PV integration level with a lower but unbalanced PV penetration.

Besides, it is worth mentioning that the phase with highest PV penetration does not necessarily have the highest voltage. As in Table 1, phase C has the highest PV penetration in Case 1 (45%), but its voltage might be the lowest among three phases (shown in Table 2). In some researches, the utility is suggested to connect new installed PV generators to the phase with lowest voltage if all PV inverters operate with a unit power factor. However, with locally dependent Var compensation methods applied on PV inverters, this connection dispatch might lead to even unbalanced PV power integration across three phases and worse overvoltage issues.

4. Problem Analysis

The limitation of locally dependent Var compensation methods in voltage regulation with unbalanced PV power integration is revealed above. In this section, the reason for this voltage regulation problem is analyzed.

4.1. Power Flow Equations

For distribution systems with 3-phase 4-wire overhead lines, the voltage drop along a feeder can be expressed as [36]where represents the phase or neutral-to-ground voltage at Bus 1; represents the phase or neutral-to-ground voltage at Bus 2; () represents the primitive self or mutual impedance of conductors; and represents the conductor current. Since neutral points are grounded at each bus (), equation (1) can be simplified by Kron reduction aswhere equivalent self and mutual impedance can be obtained from the primitive impedance in the impedance matrix of (1) as

Due to the existence of equivalent mutual impedance (, , and ), active and reactive power flow in one phase will cause voltage variations not only on its own phase conductor but also on conductors of the other two phases.

4.2. Voltage Sensitivity Analysis

To investigate the voltage regulation performance with unbalanced PV power integration, the three-phase voltage sensitivity with respect to single-phase active power injection is demonstrated in this section. Without loss of generality, a certain amount of active power is assumed to be injected into Phase C, which results in line current . At the same time, both Phases A and B remain as open circuit. Therefore, with and set to be zero in (2), the downstream three-phase voltage at Bus 2 can be expressed as

The voltage-current relationship given in (4) can be displayed by using a vector diagram shown in Figure 5(a). In this vector diagram, there are something worth pointing out: (1) balanced three-phase voltage source at Bus 1 (, , and ) is assumed; (2) with a negative value represents reverse power flow; (3) since only active power is injected into Phase C of Bus 2, and have a 180-degree phase angle difference ( and will have a 0-degree phase angle difference if only active power is absorbed from the grid); (4) has the same phase angle with , while has a 90-degree phase angle difference with .

As shown in Figure 5(a), there are different voltage responses on different phases to the same reverse-active power flow in Phase C. With active power injected into Phase C, voltage magnitude at both Phases B and C is increased. However, a voltage drop can be observed at Phase A. In this situation, with locally dependent Var compensation methods, PV inverters installed on Phase C will absorb to compensate local voltage rise.

To figure out the impact of reactive power absorption on phase voltage, similar analysis is given. In Figure 5(b), the three-phase voltage variations induced only by single-phase reactive power absorption from Phase C is displayed. In this case, and have a 90-degree phase angle difference due to pure reactive power absorption. As shown in this figure, reactive power absorption from Phase C can substantially reduce the voltage of its own phase. However, it will lift the voltage of both Phase A and Phase B. To summarize, active power injected into Phase C would lift the voltage of Phase B and reactive power absorption by PV inverts on Phase C would induce a further voltage rise on Phase B.

Three-phase voltage variations caused by single-phase active power injection or reactive power absorption are summarized as in Table 3. Now let us revisit the voltage regulation problem proposed in Section 3. In Case 1 of Table 1, much more PV inverters are installed on Phase C compared with that on Phase B. PV inverters operating with leading power factors on Phase C contribute a significant voltage rise at Phase B, which will make the voltage regulation ability of limited PV inverters on Phase B saturated. Although active power injection into Phase A can mitigate such a voltage rise at Phase B, as shown in Table 3, this effect can be feeble due to less PV power integration in Phase A in this case. As a result, overvoltage issue occurs on Phase B when PV inverters generate reactive power only relying on local measurements. While, for Case 2 of Table 1, each phase has higher but the same PV penetration (50%). Consequently, the interaction among phases can be compensated by each other. Hence, voltage rise caused by PV power injection can be sufficiently mitigated by locally dependent Var compensation methods.

To summarize, due to equivalent mutual impedance among phase conductors, active and reactive power flows in one phase can have significant impact on voltage rise or drop in other two phases. Furthermore, unbalanced PV power integration also makes locally dependent Var compensation methods less effective, which significantly challenges the effectiveness of voltage regulation methods. In order to improve the effectiveness, some suggestions are proposed.

5. Suggestions

5.1. Definition of PV Power Imbalance Index

Most previous researches highlight that high PV penetration may cause voltage regulation issues in distribution systems. However, based on Section 4, besides the PV penetration level, the distribution of PV generators across three phases is also an important factor that influences potential overvoltage problems. Therefore, an imbalance index of PV power integration is designed to describe the distribution of PV generators across three phases in this section.

Since the end of a distribution feeder is most vulnerable to severe overvoltage problems caused by excessive PV power injection, all PV power installation capacities at different buses are first converted to the end of a feeder aswhere represents the equivalent PV installation capacity of Phase after conversion; represents the actual PV installation capacity of Phase Bus ; and represent the distance from slack bus to Bus and Bus , respectively. Assume Bus is the farthest bus of the system and all line segments have the same impedance matrix. Sincewhere represents the equivalent PV current caused by equivalent PV power on phase and represents the voltage at phase . Assume , , and are , , and , respectively. By applying the symmetrical component theory, the imbalance index of PV power integration can be defined as the imbalance of equivalent PV current:whereWith such a definition, this imbalance index will be zero if PV power integration is totally balanced across three phases. On the contrary, if all PV generators are connected to one phase, the imbalance index will be one. In other cases, the possible range of the defined PV power imbalance index is (0, 1). The PV penetration level together with the PV power imbalance index can approximately reflect the integration of single-phase rooftop PV generators in low-voltage distribution systems. It is worth mentioning that the value of PV power imbalance index does not absolutely indicate the effectiveness of the reactive power compensation method in voltage regulation. The performance of voltage regulation is also affected by other factors such as mutual impedance.

5.2. Comparison of Locally Dependent Methods

Although the limitation in voltage regulation is inevitable to all locally dependent Var compensation methods, different methods have different voltage regulation performance. For example, when PV penetration is unbalanced, the power factor droop curve scheme results in lower overvoltage in comparison to that of the constant power factor scheme as in Case 1 of Table 2. Actually, when photovoltaic imbalance occurs, the power factor droop curve scheme does have a better performance in the voltage regulation than the constant power factor method.

According to the three-phase voltage sensitivity summarized in Table 3, properly reducing the reactive power absorption from Phase A and Phase C is an alternative option to mitigate the overvoltage in Phase B when the reactive power generation in Phase B becomes saturated due to the power factor constraint. Compared with the constant power factor (0.9 leading) scheme, the power factor droop curve scheme allows PV inverters to operate with a power factor higher than 0.9 before the voltage violates its limit. As a result, less reactive power will be absorbed from Phases A and C, which will lead to a lower voltage at Phase B according to Table 3. Therefore, the power factor droop curve scheme has better voltage regulation performance compared with that of the constant power factor scheme when PV penetration is unbalanced.

In order to further mitigate the overvoltage issue of Phase B, the dead band (distance between and ) of a power factor droop curve can be widened as shown in Figure 6. The corresponding voltage regulation performance is compared in Table 4 with different parameters of power factor droop curves and the same unbalanced PV power integration (Case 1 of Table 1).

As the dead band of a power factor droop curve becomes wider, less reactive power will be absorbed with the same local voltage before voltage violation. Therefore, the reactive power absorption from Phases A and C can be further reduced, which will lead to a voltage rise at both Phases A and C as in Table 4. At the same time, overvoltage at Phase B can be mitigated. To avoid the oscillation with a steep power factor droop curve, a first-order inertia element can be added in the voltage/Var controller [36].

5.3. Suggestions for Future Distribution System Planning

As revealed in Section 4, unbalanced PV power integration together with the equivalent mutual impedance among phase conductors make locally dependent Var compensation methods less effective in overvoltage mitigation. Therefore, if the mutual impedance can be properly reduced, those locally dependent methods may successfully regulate system voltage with unbalanced PV power integration.

According to (3), equivalent mutual impedance is dependent on primitive self-impedance of neutral and primitive mutual impedance , which can be expressed as in (9) and (10), respectively [37]. In these two equations, represents the resistance of neutral conductor; represents the system frequency; represents the earth resistivity; represents the geometric mean radius of conductor ; represents the distance between conductor and conductor . A proper design of the line configuration on a utility pole can reduce the equivalent mutual impedance . For example, if the distance between two phase conductor scans be properly increased and a neutral conductor with larger GMR value is selected in the future distribution system planning, smaller and can be obtained. Consequently, the equivalent mutual impedance can be correspondingly reduced according to (3). Therefore, such distribution systems can be more bearable to unbalanced PV power integration with locally dependent Var compensation methods:

6. Case Studies

A 415 V low-voltage distribution system with 64 customers shown in Figure 7 and the recorded three-phase load data shown in Figure 1(b) are used in all simulations in this section. Each bus in Figure 6 represents a utility pole which services 2 to 7 residential customers. The number in triangles represents the amount of individual houses connected to the corresponding utility pole. Voltage profiles of Bus 9 are demonstrated as results of following time-series simulations since Bus 9 is most vulnerable to voltage regulation problems.

6.1. Before PV Power Integration

415 V distribution systems were designed to operate with all possible load levels and upstream (11 kV side) voltage fluctuations. Assuming that the secondary side voltage of an 11 kV/415 V transformer varies between 1.0 pu and 1.05 pu, the corresponding three-phase voltage profiles of Bus 9 in one day with highest and lowest upstream voltages are shown in Figures 8(a) and 8(b), respectively. As shown in this figure, the 415 V system can operate within its allowable voltage range before the PV power integration. In order to focus on overvoltage issues, the voltage at the secondary side of the 415 V transformer is assumed to be 1.05 pu for all following simulations.

6.2. PV Power Integration with Different Imbalance Indices

As analyzed in this paper, besides PV power penetration, PV power imbalance index also has a significant impact on distribution system voltage regulation. If PV installation capacity has an approximately balanced distribution across three phases (24 kW, 25 kW, and 26 kW in Phase s A to C, with a total capacity of 75 kW), the imbalance index of PV power is 0.023 according to the definition in Section 5.1. In this situation, locally dependent Var compensation methods are very effective in mitigating voltage rise induced by PV power injection. For simplification, all houses are assumed to have the same PV installation capacity for the same phase, and all rooftop PV generators have the same normalized PV power profile as shown in Figure 1(a).

Figures 9(a) and 9(b) demonstrate the corresponding voltage profiles when the constant power factor (0.9 leading) scheme and the power factor droop curve scheme (, , , and are , , , and , respectively) are applied on all PV inverters, respectively. As in this figure, no overvoltage issue occurs with either locally dependent Var compensation method.

However, since rooftop PV generators are randomly distributed across three phases, it is inevitable to face unbalanced PV power integration in distribution systems. With the same total installation capacity (75 kW), if PV installation capacity in each phase is 17 kW, 24 kW, and 34 kW, respectively, the limitation of locally dependent Var compensation methods arises. In this situation, the PV power imbalance index is 0.2.

Corresponding time-series simulation results with the constant power factor scheme and the power factor droop curve scheme are shown in Figures 10(a) and 10(b), respectively. Both locally dependent Var compensation methods will make the voltage at Phase B much higher than other two phases, and overvoltage issue can be observed during noon when PV generation is high and load demand level is low. These results are consistent to the analysis in Sections 3 and 4.

6.3. Power Factor Droop Curve with a Wide Dead Band

Overvoltage problem can be mitigated in a certain extent if locally dependent Var compensation follows a power factor droop curve with a wide dead band. Figures 11(a) and 11(b) demonstrates the voltage regulation performance of power factor droop curves with different parameters. Time-series simulation results verify that power factor droop curve with a wider dead band can have a better voltage regulation performance in the situation with unbalanced PV power integration across three phases.

7. Conclusion

With random connection of rooftop PV generators, PV penetration tends to be unbalanced across three phases especially in low-voltage distribution systems. In such situations, locally dependent Var compensation methods may become less effective in overvoltage mitigation.

The voltage regulation problem is analyzed in this paper. The analysis result indicates that the equivalent mutual impedance among phase conductors together with unbalanced PV power integration make locally dependent Var compensation methods less effective.

On this basis, some suggestions are proposed for utilities in both current system operation and future planning. Specifically, (1) utilities are suggested to use both PV penetration and the PV imbalance index to describe the integration of rooftop PV generators; (2) the power factor droop curve with a wider dead band is suggested to be applied on rooftop PV generators due to its better voltage regulation performance with unbalanced PV integration; (3) for future distribution system planning, utilities are suggested to design future distribution feeders with less equivalent mutual impedance among phase conductors in order to accommodate more randomly connected rooftop PV generators.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.