Next Article in Journal
Validation of Estimating Stress from Fracture Injection Tests Using Continuous Wavelet Transform with Experimental Data
Next Article in Special Issue
A Comprehensive Review of the Design and Operation Optimization of Energy Hubs and Their Interaction with the Markets and External Networks
Previous Article in Journal
Intelligent Digital Twin Modelling for Hybrid PV-SOFC Power Generation System
Previous Article in Special Issue
A Comprehensive Tool for Scenario Generation of Solar Irradiance Profiles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Static Voltage Stability Zoning Analysis Based on a Sensitivity Index Reflecting the Influence Degree of Photovoltaic Power Output on Voltage Stability

School of Electric Power Engineering, Nanjing Institute of Technology, Nanjing 211167, China
*
Author to whom correspondence should be addressed.
Energies 2023, 16(6), 2808; https://doi.org/10.3390/en16062808
Submission received: 14 February 2023 / Revised: 9 March 2023 / Accepted: 15 March 2023 / Published: 17 March 2023
(This article belongs to the Special Issue Power System Analysis Control and Operation)

Abstract

:
The large-scale integration of photovoltaic (PV) power can bring a greatly negative influence on the grid-connected system’s voltage stability. To study the static voltage stability (SVS) of PV grid-connected systems, the traditional SVS index, L-index, was re-examined. It was firstly derived and proved that the PV active output Ppv is proportional to the voltage phase angle of the PV station’s POI (Point of Interconnection), based on a simplified two-node system integrated with a PV station operating in PV (active power—voltage) mode or PQ (active power—reactive power) mode with unit power factor. Then a novel voltage stability sensitivity index LPAS-index was proposed that takes the derivative of the L-index with respect to the POI’s voltage phase angle, so as to reflect the influence degree of Ppv on the SVS of each load node. A SVS zoning analysis method for the PV grid-connected system was designed according to the classification results of load nodes based on the proposed LPAS-index, the power grid can be zoned into three kinds of areas that reflect different correlations between the SVS and Ppv: strong correlation, moderate correlation and weak correlation. Since the LPAS-index is less impacted by Ppv, the SVS zoning results are relatively unchanged. On the basis of a classic 14-node system with PV, the practicability of the zoning analysis method was verified. The simulation results show that the PV access point in general falls within the strongly or moderately associated area with Ppv. When most of the load nodes fall within the weakly associated area with Ppv, it is not necessary to consider the impact of Ppv and load power is still the main influencing factor on the SVS. In the multi-PV case, owing to the expansion of areas more affected by Ppv, an excessive Ppv can cause adverse influence on the SVS of the whole power grid; and an effective PV power-shedding measure is proposed to solve this problem. The proposed SVS zoning analysis method can be used for reference by power grid dispatchers.

1. Introduction

In order to alleviate problems such as energy shortage and environmental pollution and achieve the goal of ‘dual carbon’, China has made great efforts to develop photovoltaic (PV) power generation technology in recent years [1,2]. With the maturity of PV grid-connected technology and the continuous reduction of installation costs, PV power generation shows a development trend of large-scale grid-connected operation and becomes the main new-energy in the construction of modern power systems [3]. In China, according to the latest data released by the National Energy Administration, the installed capacity of grid-connected PV power generation reached 392.04 GW by the end of 2022, and centralized PV accounts for about 60 percent of the total installed capacity. The PV power output is greatly affected by natural and meteorological factors and has a strong random fluctuation and intermittence [4,5]. Meanwhile, the integration of large-scale PV intensifies the power electronic characteristics and reduces the power system inertia, and brings more effects and challenges to the static and dynamic stability of the power system [6,7,8], in particular the static-state (steady-state), small-disturbance, transient-state and long-term voltage stability [9,10,11].
Static voltage stability (SVS) refers to the voltage stability when the power system sustains various small disturbances (such as small changes in load power) without considering the dynamic characteristics of each electric element [12]. In recent years, the studies on the SVS of traditional AC power systems have been still a hot spot and various machine learning algorithms are applied to the SVS prediction issue [13,14,15].
In the context of the integration of large-scale PV power stations, what are the main factors affecting the SVS of the power grid? How to determine the fluctuation range of SVS critical point, SVS domain and SVS margin; how to analyze the impact of PV power fluctuation on the weakest SVS area; how to design a novel SVS index (or criterion) that can apply to the PV power fluctuation and how to use certain appropriate control strategies to improve the system’s SVS are all the main concerns.
The research on the SVS issue of the large-scale PV station grid-connected system has been reported in relevant references. In Reference [16], the impacts of multiple factors were studied on the SVS, including PV transmission line length, PV active output, and PV topological structure, etc. In Reference [17], the impacts of several factors such as solar irradiance, PV generation power factor, PV installed capacity and PV transmission line impedance, were further studied.
In Reference [18], the SVS margin of the IEEE 14-node system with PV under the two modes of unity power factor operation and constant voltage operation was calculated, and it is considered that the integration of large-scale PV is conducive to improving the SVS margin. However, it was pointed out that too much PV active output can decrease the SVS margin in Reference [9]. In Reference [19], an assessment for the probabilistic SVS margin was studied, by using the probabilistic model of PV/wind power and the Monte-Carlo simulation method.
Several novel indexes that are suited for the SVS analysis in the system with large-scale PV stations were proposed. In Reference [20], the influence of the PCC (Point of Common Coupling) of the PV station on the SVS was studied by defining a sensitivity index of load node voltage—PCC active power that can reflect the impact of PV penetration rate on the SVS. In Reference [21], an improved NVSI-index based on the traditional IVSI-index was used to measure the SVS of PCC, which directly considers the impact of PC/wind injection power. A short-circuit ratio index was proposed in Reference [22], which is more related to the SVS of a weak sending-end system with PV. In Reference [23], a synthetical application framework was proposed to measure the SVS of a PV grid-connected system, considering the critical eigenvalue by modal analysis, the reactive power margin by QV (reactive power—voltage) analysis and the line-loss by power flow analysis.
The SVS control for the system with a large-scale PV station was also explored. In Reference [24], a SVS fuzzy controller was designed, synthetically adopting the load node voltage and the load-margin index considering PV power fluctuation as the fuzzy input variables. In Reference [25], SSSC (Static Synchronous Series Compensator) was used to improve the SVS of weak nodes considering the integration of high penetration PV/wind.
In the above references, the influence law of PV power size and fluctuation on the SVS of load nodes, load areas and the whole power grid is not studied. In the SVS analysis of large-scale PV grid-connected systems, in addition to the frequently-used small disturbance condition of gradually increasing load power, the impact of PV power output should also be considered when the system load is at a certain level.
From the above concerns, a traditional SVS index, L-index [26], is re-examined in this paper. Firstly, the relationship between the L-index and PV active output Ppv is derived. Then a novel sensitivity index of L-index of each load node relative to Ppv is proposed. According to the numerical level of the proposed sensitivity index value, a zoning method for the SVS analysis is proposed and verified. Compared with previous studies, the proposed zoning analysis method can reveal the influence degree on the SVS of different areas in the power grid by Ppv, and some effective measures can be suggested to improve the SVS according to the zoning results.

2. A Novel Sensitivity Index of L-Index Relative to PV Active Output

2.1. The Traditional L-Index

The traditional L-index is a local SVS index and was first proposed by Kessel [26], which is used to monitor and evaluate the SVS of load nodes in the traditional power system.
For a large-scale PV grid-connected system, the node-voltage equation of the system is given below:
I ˙ L I ˙ G 0 = Y LL Y LG Y LC Y GL Y GG Y GC Y CL Y CG Y CC V ˙ L V ˙ G V ˙ C
In Equation (1):
YLL, YLG, YLC, YGL, YGG, YGC, YCL, YCG and YCC correspond to the sub-matrices of the grid-connected system’s node admittance matrix, respectively.
V ˙ L and I ˙ L are the load nodes’ voltage and current vectors, respectively.
V ˙ G and I ˙ G are the power nodes’ voltage and current vectors, respectively, including slack bus, PV (active power—voltage) bus of synchronous generator, and POI (Point of Interconnection) of PV or PV converter outlet bus.
V ˙ C and I ˙ C are the contact nodes’ voltage and current vectors respectively, and the contact node is the node with neither power supply nor load demand.
By eliminating the contact nodes in Equation (1), we can get:
V ˙ L I ˙ G = H I ˙ L V ˙ G = Z LL F LG K GL Y GG I ˙ L V ˙ G
where ZLL, FLG, KGL and YGG are the block sub-matrices of the H-matrix, thereinto FLG is the load participation factor sub-matrix.
Then the L-index of each load node can be given below:
L j = L ˜ j = 1 i a G F ˜ j i V ˙ i V ˙ j
where i and j are the number of power nodes and load nodes, respectively; V ˙ i and V ˙ j are the voltage phasors of node i and node j respectively; αG is the set of power nodes, F ˜ j i is the load participation factor (complex form); L ˜ j is the complex expression of L-index.
Virtually, the L-index is in the complex form L ˜ j , and its modulus is taken as the practical index in order to measure the SVS of each load node. The range of Lj is between 0 and 1. When the value of Lj trends to 0, the SVS of load node j trends to more stability.
In general, the voltage fluctuations of load nodes and power nodes caused by the PV power fluctuation are very small, so the L-index value is less affected by the PV active output Ppv. The synchronous generator (set as PV bus) is the same.

2.2. The Relationship between PV Active Output and POI’s Voltage Phase Angle

In the SVS analysis, the PV power can be regarded as a negative load power, so that the Thevenin equivalent can be carried out from the POI of the PV station. Then a simplified PV station grid-connected two-node system is formed, as shown in Figure 1.
In Figure 1, Ppv and Qpv are the active output and reactive output of PV, respectively; Vpv and δ are the voltage amplitude and phase angle of POI, respectively. Z is the impedance modulus of the equivalent power grid line and θ is the impedance angle. E is the potential of an equivalent electric source.
In this paper, all electrical quantities are in per-unit value (pu) and the phase angle’s unit is rad.
Then we can get:
P pv j Q pv = V pv δ E 0 V pv δ Z θ * = E V pv Z cos δ + θ V pv 2 Z cos θ + j E V pv Z sin δ + θ V pv 2 Z sin θ
The followings can be obtained by arranging Equation (4):
cos δ + θ = V pv E cos θ Z E V pv P pv
sin δ + θ = V pv E sin θ Z E V pv Q pv
In normal operation of the grid-connected system, δ is a positive or negative value, and its absolute value is generally small (around 0 rad); while θ is generally large (close to π/2 rad), therefore (δ + θ) is greater than 0 rad and less than π rad.
The PV station’s operation mode can generally be divided into PV (active power—voltage) mode and PQ (active power—reactive power) mode [18,27]. PV mode is the constant voltage operation mode and the voltage amplitude of the POI or PV inverter outlet bus is set as a constant value. When Ppv changes, the POI’s voltage phase angle will change. PQ mode is the constant generation power factor operation mode. When Ppv changes, the voltage amplitude and voltage phase angle of POI will change.

2.2.1. PV Mode

In Figure 1, if the PV station operates in PV mode, Vpv is constant. According to Equation (5), with the increase of Ppv, cos(δ + θ) gradually decreases, (δ + θ) gradually increases. It can be known that δ gradually increases since θ is a constant, that is, δ is proportional to Ppv.

2.2.2. PQ Mode

If the PV station operates in PQ mode, the situation is more complicated. Now, we only consider that Ppv takes the unit power factor, that is, Qpv = 0 pu. Considering the effect of equivalent line resistance, then θ is greater than 0 rad and less than π/2 rad.
From Equation (6), we can get:
d sin δ + θ d P pv = sin θ E d V pv d P pv
Since d V pv d P pv decreases monotonically from a positive value to a negative value (over 0) with the increase of Ppv (the derivation process is given in Appendix A), sin(δ + θ) first increases and then decreases monotonically. Whereas (δ + θ) is greater than 0 rad and less than π rad, δ changes from small to large, that is, δ is proportional to Ppv.
By summarizing the above results, it can be concluded that the active output Ppv of the PV station operating in PV mode or PQ mode with unit power factor is proportional to the voltage phase angle δ of POI.

2.3. A Novel Voltage Stability Sensitivity Index LPAS

Since the PV active output, Ppv is directly proportional to the POI’s voltage phase angle δ, we can take the derivative of L-index with respect to δ to reflect the influence degree of Ppv on the SVS of each load node.
Set F ˜ j i = F j i α j i , V i . = V i δ i and V j . = V j δ j , then Equation (3) can be expressed below:
L j = 1 1 V j i a G F j i V i [ cos ( α j i + δ i δ j ) + j sin ( α j i + δ i δ j ) ]
The sensitivity of the L-index of load node j relative to the voltage phase angle of power node i can be obtained by taking the partial derivative of the complex expression of L-index ( L ˜ j ) with respect to the voltage phase angle δi and taking the modulus value. According to the verification, the partial derivative of the modulus expression of L-index (Lj) with respect to δi is the same as it, so we can obtain:
L j δ i = L ˜ j δ i = F j i V i V j sin ( α j i + δ i δ j ) j cos ( α j i + δ i δ j ) = F j i V i V j
Equation (9) represents the coupling degree between the L-index of each load node and the voltage phase angle of each power node (including synchronous generators and PV stations). The influence of voltage phase angle change Δδi of each power node on the L-index of each load node can be expressed below:
Δ L 1 Δ L 2 Δ L m = L 1 δ 1 L 1 δ 2 L 1 δ n L 2 δ 1 L 2 δ 2 L 2 δ n L m δ 1 L m δ 2 L m δ n Δ δ 1 Δ δ 2 Δ δ n
where m is the maximum number of load nodes and n is the maximum number of power nodes.
For a PV station in the system, set the number of POI as k and its voltage amplitude as Vpvk. Then, a voltage stability sensitivity index (named LPAS-index), namely the derivative of L-index with respect to POI’s voltage phase angle can be defined below:
L P A S j = F j k V pv k V j
where j is the number of load nodes.
Equation (11) reflects the sensitivity of the L-index of each load node relative to POI’s voltage phase angle of each PV station, that is, the sensitivity of the L-index relative to Ppv. By calculating the LPAS-index, the relationship between the SVS of each load node and Ppv can be explored.
In a known power grid, if the power grid’s structure does not change, the load participation factor Fjk will remain unchanged. For a PV station operating in PV mode, since Vpvk is constant, the LPAS-index value is mainly affected by the load node voltage Vj. For a PV station operating in PQ mode with a unit power factor, the LPAS-index value is affected by the POI’s voltage Vpvk and the load node voltage Vj at the same time. It is noteworthy that the voltage fluctuations of load nodes and POI caused by the PV power fluctuation are very small normally, so the LPAS-index value is less affected by Ppv. The same is true for the synchronous generator (set as PV bus).

3. Static Voltage Stability Zoning Analysis Method

In the system with a large-scale PV station, if the PV station operates in PV mode or PQ mode with a unit power factor, for each load node, the sensitivity of the L-index relative to Ppv can be obtained by calculating the corresponding LPAS-index. According to the numerical level of the LPAS-index value, the whole system can be zoned into several areas which can reflect the correlations between the SVS and Ppv. On the other hand, the weakest SVS area and the weakest node can be determined by calculating the L-index value. Then according to the zoning results, the SVS analysis and control aiming at the weakest area and the whole power grid can be conducted.
Figure 2 is the flow chart of the SVS zoning analysis method based on the LPAS-index.
The main steps of the SVS zoning analysis method are as follows:
(1)
Get the sets of power nodes, load nodes and contact nodes by assessing the node types, and calculate the admittance matrix of the grid-connected system, find the sub-matrix FLG of load participation factor and obtain the participation factor of each load node relative to the POI of each PV station.
(2)
According to Equations (3) and (11), the L-index and LPAS-index of each load node can be calculated. By ranking the L-index values of all load nodes, determine the weakest SVS area and the weakest node.
(3)
Load nodes are classified according to the numerical level of LPAS-index value, which can be generally classified into three types of nodes whose SVS is greatly, moderately and less affected by Ppv. According to the classification results for load nodes, the PV grid-connected system can be zoned into three kinds of SVS areas that are strongly, moderately and weakly associated with Ppv.
The numerical level of LPAS-index value can be classified according to the standard: greater than 50% (corresponding to the strong association area), from 20% to 50% (corresponding to the moderate association area), and less than 20% (corresponding to the weak association area). However, the classification standard is not unchangeable and should be determined according to the actual situation of the PV grid-connected system.
(4)
According to the classification results for load nodes, the PV grid-connected system can be zoned into three kinds of SVS areas that are strongly, moderately and weakly associated with Ppv. As the LPAS-index value of the load node is less affected by Ppv and other generators, the zoning results are relatively unchanged.
(5)
Analyze the SVS of the whole power grid and the weakest SVS area on the basis of the zoning results, and find out the rules that are affected by Ppv. Based on the analysis results, some useful control strategies can be proposed to improve the SVS of the weakest area and the whole system.
If most load nodes fall within the weakly associated area with Ppv, the impact of Ppv on the SVS can be neglected, and the change of load power is the main influencing factor on the SVS of the whole power grid. Similar to the traditional power grid, reactive power compensation or other effective measures can be used to improve the SVS.
If more load nodes fall within the strongly or moderately associated area with Ppv, the impact of Ppv on the SVS has to be valued. Load-shedding or other effective measures can be taken when necessary, and we will use a new PV power-shedding measure.

4. Simulation Verification for the Zoning Analysis Method

In order to verify the rationality and practicability of the above zoning analysis method, we take the IEEE 14-node system as the test example, as shown in Figure 3. The 14-node system is a sub-transmission system [18], and its voltage classes are 69 kV (including Node 1~Node 5), 13.8 kV (including Node 6, Node 7, and Node 9~Node 14), and 18 kV (Node 8).
See Figure 3, Node 1 is the slack bus (swing bus), Node 2 is a PV bus. SC1~SC3 are synchronous condensers. Set the base power as 100 MVA and the initial power of the total load as 2.59 + j0.814 pu. Set the voltage amplitude of Node 2 as 1.045 pu and the initial active output of generator Gen 2 (PG2) as 1 pu.
See the dotted line in Figure 3, three plans about the centralized PV station integrated into the 14-node system will be adopted:
  • Plan A: PV station A is integrated into Node 5 by a PV transmission line, and the corresponding POI is POI_A;
  • Plan B: PV station B is integrated into Node 14 by a PV transmission line, and the corresponding POI is POI_B;
  • Plan C: PV station B is integrated into Node 14, the generator Gen 2 is replaced with PV station C (operating in PV mode), and the corresponding POI is Node 2. This is a multi-PV case.
PpvA, PpvB and PpvC represent the active output of PV stations A, B and C, respectively, and the maximum active output of them are 1 pu, 0.3 pu and 1 pu, respectively.

4.1. Verification for the Relationship between PV Active Output and Voltage Phase Angle of POI

Section 2.2 proved that the PV active output Ppv is proportional to the POI’s voltage phase angle δ, now we take Plan A as an example to verify it. When PV station A operates in PV mode, the PV inverter outlet bus voltage amplitude is set as 1.05 pu. When PV station A operates in PQ mode with a unit power factor, its reactive power is set as 0 pu. The PpvAδA (δA is the voltage phase angle of POI_A.) curve is shown in Figure 4, we can see that PpvA is proportional to δA under the two operation modes.

4.2. Index Calculation

4.2.1. Plan A and Plan B

Firstly, Plan A and Plan B are investigated. Now take PV stations A and B operating in PV mode as an example. Set PpvA = 1 pu, PpvB = 0.3 pu (The capacity of PV station integrated into 13.8 kV node can’t be too large). Set the current load multiple of the whole system lm = 1 pu. Calculate the L-index values of all load nodes when the PV station adopts Plan A and Plan B, respectively, and does not integrate into the 14-node system, the results are listed in Table 1.
It can be seen from Table 1 that when the PV station is not integrated into the system, Node 14 is the weakest node and Node 5 is the strongest node. Node 9, Node 10 and Node 14 compose the weakest SVS area.
When Plan A is carried out, the weakest SVS area is still composed of Node 9, Node 10 and Node 14, Node 14 and Node 5 are still the weakest node and the strongest node.
When Plan B is adopted, the weakest SVS area is composed of Node 9 and Node 10. The weakest node becomes Node 9, this is because the SVS of Node 14 is observably enhanced after integrating with PV station B. The strongest node is still Node 5.
It can be seen that the change of Ppv has a minor impact on the L-index from Figure 5 (Node 5 in Plan A and Node 14 in Plan B are taken as examples). According to the computing results, the variable amplitude of the L-index is less than 0.1%. Therefore, when Ppv changes, the determined weakest SVS area can remain consistent.
See Table 1, when a PV station is integrated into the system, compared with the system without a PV station, the L-index values of all load nodes are relatively reduced, indicating that the SVS of load nodes is improved to a certain extent. Moreover, when the PV station is integrated into different load nodes, the weakest SVS area and the weakest node can be changed.
As shown in Table 2, the LPAS-index values of all load nodes, respectively, relative to POI_A and POI_B are sorted from large to small under Plan A and Plan B.
Figure 6 shows the Ppv—LPAS-index curves (Node 5 in Plan A and Node 14 in Plan B are taken for examples). To the same load node, the change of Ppv has a minor impact on the LPAS-index value. Other active power supplies are similar such as a generator. In a general way, the variable amplitude of the LPAS-index is less than 1% when Ppv changes. Therefore, aiming at the changes of Ppv, the classification of load nodes and the SVS zoning results can remain unchanged.
In summary, the maximum active output of the PV station can be directly selected to calculate the LPAS-index and L-index of load nodes, so as to carry out the classification of load nodes and the determination of the weakest SVS area.

4.2.2. Plan C

In Plan C, PV station B has remained, and Gen 2 is replaced with PV station C. PV station C operates in PV mode, and the voltage of Node 2 remains as 1.045 pu. When PpvB = 0.3 pu, PpvC = 1 pu, the values of L-index and LPAS-index relative to Node 2 (namely POI of PV station C) and POI_B are listed in Table 3. Node 9 is the weakest node and Node 12 is the strongest node. Node 9 and Node 10 compose the weakest SVS area.
Similarly, the active output of PV station C (PpvC) has a minor impact on the values of the L-index and LPAS-index. For example, in Plan C, the LPAS-index values of Node 5 relative to Node 2 and Node POI_B are 0.3859 and 0.0193, respectively, when PpvC = 0.3 pu and PpvB = 0.3 pu, and are approximately equal to the corresponding LPAS-index values when PpvC = 1 pu and PpvB = 0.3 pu, as shown in Table 3.

4.3. SVS Zoning and Analysis

According to the data in Table 2 and Table 3, the load nodes can be classified and the 14-node system with PV can be zoned.

4.3.1. Plan A

See Table 2, in Plan A, only the LPAS-index value of Node 5 is greater than 20% and less than 50%, illustrating that the SVS of Node 5 is moderately affected by PpvA. The LPAS-index values of other load nodes (including Node 4, Node 9, Node 10, Node 11, Node 12, Node 13 and Node 14) are all less than 20%, illustrating that the SVS of these nodes are all less affected by PpvA.
Furthermore, the LPAS-index value of Node 5, namely the access point of PV station A, is not too big, indicating that the PV access point is not necessarily greatly affected by Ppv.
Figure 7 shows the zoning chart of Plan A. The system can be zoned as a moderately associated area (including Node 5) and a weakly associated area (including other load nodes) with PpvA. Owing to the LPAS-index and L-index are less affected by Ppv and PG2, the zoning results of Plan A are relatively unchanged.
See Figure 7, the weakest SVS area (composed of Node 9, Node 10 and Node 14) falls within the area that is weakly associated with PpvA, so this area’s SVS is less affected by PpvA. In fact, since the overwhelming majority of load nodes (except Node 5) are falls with the weakly associated area with PpvA, the SVS of the whole power grid is less affected by PpvA, and more attention should focus on the impact of load power change on the SVS.
As the impact of Ppv on the L-index is minimal, we select the system load-margin index ILM to measure the SVS and verify the above conclusion. The load-margin index ILM can be calculated by Equation (12):
I LM = 1 1 λ MAX
where λMAX is the maximum load margin parameter that corresponds to the critical point of the SVS.
The value range of ILM is 0~1. 0 means voltage collapse, and 1 means absolute voltage stability.
Now set PG2 = 0.4 pu. The values of ILM are given in Table 4. When the system load multiple lm is 1 pu, the system operates in a normal state (namely a low load level). When PpvA changes from 0 to 1 pu, the ILM values keep almost unchanged (Adjoining variation amplitude is less than 1%).
When the load multiple lm is 1.5 pu, the system operates in a heavy load (namely a high load level). See Table 4, with the change of PpvA, the ILM values still keep almost unchanged. However, compared to the normal operation state (lm = 1 pu), the ILM values and the SVS margin significantly reduce, indicating that the SVS of the whole system (including the weakest area) is mainly affected by the load power. The load-shedding method can be used to improve the SVS [24].

4.3.2. Plan B

See Table 2, in Plan B, the LPAS-index value of Node 14 (the access point of PV station B) is greater than 50%, so the SVS of Node 14 is greatly affected by PpvB. However, the LPAS-index values of other load nodes (including Node 4, Node 5, Node 9, Node 10, Node 11, Node 12 and Node 13) are all less than 20%, so the SVS of these nodes is less affected by PpvB.
Figure 8 shows the zoning chart of Plan B. Node 14 can be set as a strongly associated area with PpvB, and other load nodes can be set as a weakly associated area with PpvB. The weakest SVS area (composed of Node 9 and Node 10) falls within the area that is weakly associated with PpvB.
See Table 5, with the change of PpvB, the variation amplitudes of ILM value are still very small under normal state and heavy load, and the SVS of the whole power grid is mainly affected by the load power. It is similar to the case in Plan A, more attention should focus on the impact of load power change on the SVS.

4.3.3. Plan C

See Table 3, in Plan C, the LPAS-index values of Node 4 and Node 5 relative to PpvC is greater than 20% and less than 50%, so the two nodes can compose a moderately associated area with PpvC. The LPAS-index values of other load nodes relative to PpvC are all less than 20%, so they can compose a weakly associated area with PpvC. However, the LPAS-index value of Node 14 relative to PpvB is greater than 50%, so we can neglect the impact of PpvC and zone Node 14 as a strongly associated area with PpvB. Ultimately, Node 9~Node 13 can be determined to compose a weakly associated area with PpvB and PpvC. The zoning chart of Plan C is shown in Figure 9.
The weakest SVS area (composed of Node 9 and Node 10) falls within the weakly associated area with PpvB and PpvC, so the impact of Ppv on the weakest area can be neglected.
See Figure 9, what is different from Plan A and Plan B is that the areas more affected by Ppv are significantly enlarged in Plan C. Especially, Node 4 and Node 5 (69 kV voltage class) fall within a moderately associated area with PpvC, so the impact of Ppv on the SVS should be further investigated.
Now we set PpvB = 0.3 pu, and PpvC changes from 0 to 1 pu.
The ILM values are listed in Table 6. It can be seen that when the system operates in a normal state (lm = 1 pu), with the gradual increase of PpvC from 0 to 0.5 pu, the ILM value gradually increases too, indicating that the SVS is gradual enhanced with the increase of PpvC. However, with the gradual increase of PpvC from 0.6 pu to 1 pu, the ILM value gradually decreases and the SVS is gradually weakened. When PpvC reaches the maximum (1 pu), the impact of Ppv on the ILM value is very obvious (decreases to 0.4597), and the SVS becomes as weak as the status that the system operates in heavy load.
See Table 6, when the system operates under heavy load (lm = 1.5 pu), the ILM value gradually increases with the increase of PpvC from 0 to 0.6 pu, and gradually decreases with the increase of PpvC from 0.7 pu to 1 pu, but the variation amplitude is very small.
Now we set PpvC = 1 pu, and PpvB changes from 0 to 0.3 pu. The ILM values are listed in Table 7. It can be seen that excessive PV power can evidently weaken the SVS when the system operates in a normal state (lm = 1 pu). When the system operates under heavy load (lm = 1.5 pu), the increase of PpvB only can slightly improve the SVS and has no negative impact on the SVS.
From the above analyses, it is concluded that in a multi-PV case, excessive PV active power can evidently weaken the SVS when the system operates in a normal state, and only slightly impact the SVS when the system operates under heavy load.
A PV power-shedding method can be used in time to maintain the SVS when the system operates at a low load level and an excessive PV power. For example, when lm = 1 pu, PpvB = 0.3 pu and PpvC = 0.95 pu, the value of ILM is 0.4958 by calculation. Then the PV active power control system can execute a series of power-shedding operations and the shedding quantity is 0.1 pu each time, PpvC will be reduced to 0.55 pu, and the value of ILM will recover to 0.6180. The system’s SVS is improved effectively after power-shedding.
Moreover, reducing the rated installation capacity of PV station is also a valid solution to this problem.

5. Conclusions

Firstly, this paper derives and proves that the POI’s voltage phase angle of the PV station is proportional to the PV active output Ppv, based on a simplified two-node system with PV. From this, a novel sensitivity index LPAS-index of the traditional L-index relative to the voltage phase angle of POI is proposed, and the LPAS-index can be used to reflect the influence degree of Ppv on the SVS of load nodes. A PV grid-connected system can be zoned into different areas: strongly, moderately and weakly associated with Ppv for the SVS according to the numerical level of LPAS-index value. Based on the zoning results, the SVS analysis method is applied, by taking the classic 14-node system integrated with a centralized PV station as an example. The following conclusions can be obtained:
(1)
The LPAS-index value of one load node is less affected by Ppv or PG (active output of generator), so the SVS zoning results are relatively unchanged. However, by changing the location of the PV access point or the numerical classification standard of the LPAS-index value, the SVS zoning results can change.
(2)
The access point of the PV station is not always greatly affected by Ppv, so it can fall within the strongly associated area with Ppv, and also fall within the moderately associated area with Ppv.
(3)
If most of the load nodes except for the PV access point fall within the weakly associated with Ppv, the impact of Ppv on the SVS can be neglected and more attention should be focused on the impact of load power.
(4)
In the multi-PV case, more load nodes may fall within the areas more affected by Ppv. If no excessive PV power flows into the power grid, the increase of Ppv can improve the SVS to a certain extent. However, excessive PV active power can evidently weaken the SVS when the system operates at a low load level, and PV power-shedding can make the system maintain the SVS. On the other hand, excessive PV active power has a minor impact on the SVS when the system operates under heavy load.
(5)
In addition, when the system load is constant, the change in Ppv has a minor impact on the L-index. It is necessary to design a novel SVS index for load nodes that can better adapt to the fluctuation of Ppv in the follow-up work.

Author Contributions

Conceptualization, S.L.; methodology, S.L., Y.L. and Y.G.; software, S.L., Y.L. and Y.G.; validation, S.L. and Y.L.; formal analysis, S.L.; investigation, S.L. and Y.L.; resources, S.L.; data curation, S.L.; writing—original draft preparation, S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Scientific Research Foundation of Nanjing Institute of Technology (ZKJ202102), and the College Student Challenge Cup Support Project of Nanjing Institute of Technology (TZ20190008).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Vectors, matrices and sets
V ˙ L , V ˙ G , V ˙ C Voltage vectors of load nodes, power nodes and contact nodes
I ˙ L , I ˙ G , I ˙ C Current vectors of load nodes, power nodes and contact nodes
YLL, YLG, YLC,
YGL, YGG, YGC,
YCL, YCG, YCC
Sub-matrices of power system’s node admittance matrix
HH-matrix generated from node admittance matrix by a partial inversion
ZLL, FLG, KGLSub-matrices of H-matrix
αGSet of power nodes
Parameters and variables
i, jNumber of power nodes and load nodes
Lj, L ˜ j Modulus form and complex form of L-index of load node j
V ˙ i , V ˙ j Voltage phasors of power node i and load node j [pu, pu]
Vi, δiVoltage amplitude and phase angle of power node i [pu, rad]
Vj, δjVoltage amplitude and phase angle of load node j [pu, rad]
Fji, F ˜ j i Modulus form and complex form of load participation factor of load node j relative to power node i
Ppv, QpvPV active power output and reactive power output [pu, pu]
Vpv, δVoltage amplitude and phase angle of PV station’s POI in a Thevenin equivalent two-node system [pu, rad]
Z, θImpedance modulus and angle of equivalent line [pu, rad]
EPotential of equivalent electric source [pu]
m, nMaximum number of load nodes and maximum number of power nodes
kNumber of PV stations
VpvkVoltage amplitude of POI of PV station k [pu]
FjkLoad participation factor of load node j relative to POI of PV station k
LPASjLPAS-index of load node j
lmCurrent load multiple of the whole power grid [pu]
ILMLoad-margin index
λMAXMaximum load margin parameter [pu]
PGActive power of synchronous generator [pu]
Abbreviations
SVSStatic voltage stability
PVPhotovoltaic, photovoltaic station
POIPoint of Interconnection of PV station
PCCPoint of Common Coupling of PV station
PV busPower bus (node) with constant active power output and voltage amplitude
PV modeOperation mode of PV station with constant active power output and voltage amplitude
PQ modeOperation mode of PV station with constant active power output and reactive power output

Appendix A

The derivation process for monotonicity of d V pv d P pv is as follows.
After eliminating δ from Equations (5) and (6), we can get:
V pv 4 2 Z cos θ P pv + E 2 V pv 2 + Z 2 P pv 2 = 0
Solving Equation (A1) and removing the unreasonable solution, we get:
V pv = 2 Z cos θ P pv + E 2 + Δ 2
where, Δ = 4 Z 2 sin 2 θ P pv 2 + 4 Z E 2 cos θ P pv + E 4 .
Δ should be greater than or equal to 0, when Δ = 0, the system is at the critical point of voltage collapse.
Now set a = 4 Z 2 sin 2 θ , b = 4 Z E 2 cos θ , and set Δ > 0, then we can obtain:
d V pv 2 d P pv = 2 Z cos θ + 2 a P pv + b 2 Δ
d d P pv ( d V pv 2 d P pv ) = 4 a E 4 b 2 4 Δ Δ < 0
It can be seen from Equations (A3) and (A4) that d V pv 2 d P pv decreases monotonically from positive to negative (over 0) with the increase of Ppv.
Since d V pv 2 d P pv = 2 V pv d V pv d P pv , and Vpv > 0, we can obtain the conclusion that d V pv d P pv and d V pv 2 d P pv have the same monotonicity with the increase of Ppv.

References

  1. Cheng, L.; Zang, H.; Wei, Z.; Sun, G. Ultra-short-term forecasting of regional photovoltaic power generation considering multispectral satellite remote sensing data. Proc. CSEE 2022, 42, 7451–7464. [Google Scholar]
  2. Yao, S.; Zhang, C.; Liu, G.; Ma, J.; Wang, Y. Electromagnetic transient decoupling of photovoltaic power generation units and fast simulation method. Autom. Electr. Power Syst. 2022, 46, 170–178. [Google Scholar]
  3. Ding, M.; Xu, Z.; Wang, W.; Wang, X.; Song, Y.; Chen, D. A review on China’s large-scale PV integration: Progress, challenges and recommendations. Renew. Sust. Energ. Rev. 2016, 53, 639–652. [Google Scholar] [CrossRef]
  4. Wang, F.; Li, J.; Zhen, Z.; Wang, C.; Ren, H.; Ma, H.; Zhang, W.; Huang, L. Cloud feature extraction and fluctuation pattern recognition based ultrashort-term regional PV power forecasting. IEEE Trans. Ind. Appl. 2022, 58, 6752–6767. [Google Scholar] [CrossRef]
  5. Martins, C.C.; Sperandio, M.; Welfer, D. Cloud variance impact on the voltage and power flow of distribution networks with photovoltaic generators. IEEE Lat. Am. Trans. 2021, 19, 217–225. [Google Scholar] [CrossRef]
  6. Tamimi, B.; Canizares, C.; Bhattacharya, K. System stability impact of large-scale and distributed solar photovoltaic generation: The case of Ontario, Canada. IEEE Trans. Sustain. Energy 2013, 4, 680–688. [Google Scholar] [CrossRef]
  7. Eftekharnejad, S.; Vittal, V.; Heydt, G.T.; Keel, B.; Loehr, J. Impact of increased penetration of photovoltaic generation on power systems. IEEE Trans. Power Syst. 2013, 28, 893–901. [Google Scholar] [CrossRef]
  8. Liang, Z.; Qin, F.; Cui, F. Impact analysis of annular solar eclipse on June 21, 2020 in China on photovoltaic power generation and power grid operation. Autom. Electr. Power Syst. 2021, 45, 1–7. [Google Scholar]
  9. Li, S.; Wei, Z.; Sun, G.; Gao, P.; Xiao, J. Voltage stability bifurcation of large-scale grid-connected PV system. Electr. Power Autom. Equip. 2016, 36, 17–23. [Google Scholar]
  10. Lammert, G.; Premm, D.; Ospina, L.D.P.; Boemer, J.C.; Braun, M.; Cutsem, T.V. Control of photovoltaic systems for enhanced short-term voltage stability and recovery. IEEE Trans. Energy Conver. 2019, 34, 243–254. [Google Scholar] [CrossRef] [Green Version]
  11. Munkhchuluun, E.; Meegahapola, L. Impact of the solar photovoltaic (PV) generation on long-term voltage stability of a power network. In Proceedings of the 2017 IEEE Innovative Smart Grid Technologies—Asia (ISGT-Asia), Auckland, New Zealand, 4–7 December 2017; pp. 1–6. [Google Scholar]
  12. Kundur, P. Voltage stability. In Power System Stability and Control; McGraw-Hill: New York, NY, USA, 1994; pp. 990–1006. [Google Scholar]
  13. Tang, Y.; Dong, S.; Zhu, C.; Wu, J.; Song, Y. Online prediction method of static voltage stability margin based on Tri-Training-LASSO-BP network. Proc. CSEE 2020, 40, 3824–3834. [Google Scholar]
  14. Li, J.; Liu, D.; An, J.; Li, Z.; Yang, H.; Zhao, G.; Yang, S.; Zheng, H. Static voltage stability margin assessment based on reinforcement learning theory. Proc. CSEE 2020, 40, 5136–5147. [Google Scholar]
  15. Liu, Y.; Wang, T.; Qiu, G.; Wei, W.; Zhou, B.; Liu, T.; Liu, J.; Mei, S. Surrogate modeling method and its analytical algorithm for static voltage stability control embedded with input convex neural network. Electr. Power Autom. Equip. 2023, 43, 151–159. [Google Scholar]
  16. Wang, D.; Yuan, X.; Zhao, M.; Qian, Y. Impact of large-scale photovoltaic generation integration structure on static voltage stability in China’s Qinghai province network. J. Eng. 2017, 2017, 671–675. [Google Scholar] [CrossRef]
  17. Du, X.; Zhou, L.; Guo, K.; Yang, M.; Liu, Q. Static voltage stability analysis of large-scale photovoltaic plants. Power Syst. Technol. 2015, 39, 3427–3434. [Google Scholar]
  18. Kabir, S.; Nadarajah, M.; Bansal, R. Impact of large scale photovoltaic system on static voltage stability in sub-transmission network. In Proceedings of the 2013 IEEE ECCE Asia Downunder, Melbourne, VIC, Australia, 3–6 June 2013; pp. 468–473. [Google Scholar]
  19. Rawat, M.S.; Vadhera, S. Probabilistic steady state voltage stability assessment method for correlated wind energy and solar photovoltaic integrated power systems. Energy Technol. 2020, 9, 2000732. [Google Scholar] [CrossRef]
  20. Liu, Y.; Yao, L.; Liao, S.; Xu, J.; Dun, Y.; Cheng, F.; Cui, W.; Li, F. Study on the impact of photovoltaic penetration on power system static voltage stability. Proc. CSEE 2022, 42, 5484–5496. [Google Scholar]
  21. Li, S.; Ma, Y.; Zhang, C.; Yu, J. Static voltage stability prediction for PCCs of new energy power stations in the grid-connected system. J. Phys. Conf. Ser. 2022, 2369, 012037. [Google Scholar] [CrossRef]
  22. Xiao, F.; Han, M.; Tang, X.; Zhang, X. Voltage stability of weak sending-end system with large-scale grid-connected photovoltaic power plants. Electr. Power 2020, 53, 31–39. [Google Scholar]
  23. Rahman, S.; Saha, S.; Islam, S.N.; Arif, M.T.; Mosadeghy, M.; Haque, M.E.; Oo, A.M.T. Analysis of power grid voltage stability with high penetration of solar PV systems. IEEE Trans. Ind. Appl. 2021, 57, 2245–2257. [Google Scholar] [CrossRef]
  24. Li, S.; Wei, Z.; Ma, Y. Fuzzy load-shedding strategy considering photovoltaic output fluctuation characteristics and static voltage stability. Energies 2018, 11, 779. [Google Scholar] [CrossRef] [Green Version]
  25. Zheng, L.; Ma, D.; Zhou, X.; Dai, J. Static voltage stability characteristics of high permeability new energy grid considering SSSC. Power DSM 2022, 24, 50–56. [Google Scholar]
  26. Kessel, P.; Glavitsch, H. Estimating the voltage stability of a power system. IEEE Trans. Power Deliv. 1986, 1, 346–354. [Google Scholar] [CrossRef]
  27. Li, S.; Wei, Z.; Sun, G.; Wang, Z. Stability research of transient voltage for multi-machine power systems integrated large-scale PV power plant. Acta Energ. Sol. Sin. 2018, 39, 3356–3362. [Google Scholar]
Figure 1. The simplified PV grid-connected two-node system.
Figure 1. The simplified PV grid-connected two-node system.
Energies 16 02808 g001
Figure 2. Flow chart of the SVS zoning analysis method.
Figure 2. Flow chart of the SVS zoning analysis method.
Energies 16 02808 g002
Figure 3. The classic IEEE 14-node system with PV station.
Figure 3. The classic IEEE 14-node system with PV station.
Energies 16 02808 g003
Figure 4. PpvAδA (active output of PV station A—voltage phase angle of POI_A) curve.
Figure 4. PpvAδA (active output of PV station A—voltage phase angle of POI_A) curve.
Energies 16 02808 g004
Figure 5. PV active output Ppv—L-index curves (Load multiple lm = 1 pu).
Figure 5. PV active output Ppv—L-index curves (Load multiple lm = 1 pu).
Energies 16 02808 g005
Figure 6. PV active output Ppv—LPAS-index curves (Load multiple lm = 1 pu).
Figure 6. PV active output Ppv—LPAS-index curves (Load multiple lm = 1 pu).
Energies 16 02808 g006
Figure 7. Static voltage stability zoning chart of Plan A.
Figure 7. Static voltage stability zoning chart of Plan A.
Energies 16 02808 g007
Figure 8. Static voltage stability zoning chart of Plan B.
Figure 8. Static voltage stability zoning chart of Plan B.
Energies 16 02808 g008
Figure 9. Static voltage stability zoning chart of Plan C.
Figure 9. Static voltage stability zoning chart of Plan C.
Energies 16 02808 g009
Table 1. L-index values of load nodes under Plan A and Plan B (Load multiple lm = 1 pu).
Table 1. L-index values of load nodes under Plan A and Plan B (Load multiple lm = 1 pu).
Node NumberL-Index Value
Plan A (PpvA = 1 pu)Plan B (PpvB = 0.3 pu)Without PV
140.07600.02970.0784
90.06460.05320.0681
100.06190.05260.0649
110.03510.03030.0366
130.03160.02090.0321
40.02930.02830.0307
120.02380.01850.0241
50.01970.01940.0209
Table 2. LPAS-index values of load nodes under Plan A and Plan B (Load multiple lm = 1 pu).
Table 2. LPAS-index values of load nodes under Plan A and Plan B (Load multiple lm = 1 pu).
Plan A (PpvA = 1 pu)Plan B (PpvB = 0.3 pu)
Node NumberLPAS-IndexNode NumberLPAS-Index
523.17%1461.23%
414.18%918.71%
95.95%1015.48%
104.92%1314.05%
143.78%117.85%
112.49%126.97%
130.85%43.15%
120.42%51.93%
Table 3. The values of L-index and LPAS-index of load nodes under Plan C (Load multiple lm = 1 pu).
Table 3. The values of L-index and LPAS-index of load nodes under Plan C (Load multiple lm = 1 pu).
Plan C (PpvB = 0.3 pu and PpvC = 1 pu)
Node NumberL-IndexLPAS-Index (to Node 2)LPAS-Index (to POI_B)
50.019438.58%1.93%
40.028338.18%3.15%
90.053214.33%18.71%
100.052611.86%15.48%
110.03036.01%7.85%
140.02973.99%61.23%
130.02090.92%14.05%
120.01850.45%6.97%
Table 4. Values of system load-margin index (Plan A, active output of Gen 2 PG2 = 0.4 pu).
Table 4. Values of system load-margin index (Plan A, active output of Gen 2 PG2 = 0.4 pu).
PpvA/puLoad Multiple lm = 1 puLoad Multiple lm = 1.5 pu
λMAX/puILMλMAX/puILM
02.58460.6131 1.73460.4235
0.12.59570.6147 1.74440.4267
0.22.60680.6164 1.75360.4297
0.32.61630.6178 1.76250.4326
0.42.62510.6191 1.77080.4353
0.52.63300.62031.77870.4378
0.62.64000.62121.78610.4401
0.72.64600.62211.79210.4420
0.82.65080.62281.79760.4437
0.92.65420.62321.80320.4454
12.65600.62351.80850.4471
Table 5. Values of system load-margin index (Plan B, active output of Gen 2 PG2 = 0.4 pu).
Table 5. Values of system load-margin index (Plan B, active output of Gen 2 PG2 = 0.4 pu).
PpvB/puLoad Multiple lm = 1 puLoad Multiple lm = 1.5 pu
λMAX/puILMλMAX/puILM
02.52720.60431.69610.4104
0.12.56950.6108 1.72390.4199
0.22.60890.6167 1.75040.4287
0.32.64590.62211.77510.4367
Table 6. Values of system load-margin index (Plan C, active output of PV station B PpvB = 0.3 pu).
Table 6. Values of system load-margin index (Plan C, active output of PV station B PpvB = 0.3 pu).
PpvC/puLoad Multiple lm = 1 puLoad Multiple lm = 1.5 pu
λMAX/puILMλMAX/puILM
02.58460.6131 1.73780.4246
0.12.61380.6174 1.75050.4287
0.22.63430.6204 1.76080.4321
0.32.64510.6219 1.76920.4348
0.42.64590.6221 1.77510.4367
0.52.63210.6201 1.77870.4378
0.62.59600.6148 1.77910.4379
0.72.51710.6027 1.77610.4370
0.82.34010.5727 1.76850.4345
0.92.13110.5308 1.75570.4304
11.85070.4597 1.73430.4234
Table 7. Values of system load-margin index (Plan C, active output of PV station C PpvC = 1 pu).
Table 7. Values of system load-margin index (Plan C, active output of PV station C PpvC = 1 pu).
PpvB/puLoad Multiple lm = 1 puLoad Multiple lm = 1.5 pu
λMAX/puILMλMAX/puILM
02.10280.5244 1.67600.4033
0.12.03030.5075 1.69780.4110
0.21.95830.4894 1.71730.4177
0.31.85070.4597 1.73430.4234
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, S.; Lu, Y.; Ge, Y. Static Voltage Stability Zoning Analysis Based on a Sensitivity Index Reflecting the Influence Degree of Photovoltaic Power Output on Voltage Stability. Energies 2023, 16, 2808. https://doi.org/10.3390/en16062808

AMA Style

Li S, Lu Y, Ge Y. Static Voltage Stability Zoning Analysis Based on a Sensitivity Index Reflecting the Influence Degree of Photovoltaic Power Output on Voltage Stability. Energies. 2023; 16(6):2808. https://doi.org/10.3390/en16062808

Chicago/Turabian Style

Li, Sheng, Yuting Lu, and Yulin Ge. 2023. "Static Voltage Stability Zoning Analysis Based on a Sensitivity Index Reflecting the Influence Degree of Photovoltaic Power Output on Voltage Stability" Energies 16, no. 6: 2808. https://doi.org/10.3390/en16062808

APA Style

Li, S., Lu, Y., & Ge, Y. (2023). Static Voltage Stability Zoning Analysis Based on a Sensitivity Index Reflecting the Influence Degree of Photovoltaic Power Output on Voltage Stability. Energies, 16(6), 2808. https://doi.org/10.3390/en16062808

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop