Abstract

With the improvement of the permeability of wind and photovoltaic (PV) energy, it has become one of the key problems to maintain the small-signal stability of the power system. Therefore, this paper analyzes the small-signal stability in a power system integrated with wind and solar energy. First, a mathematical model for small-signal stability analysis of power systems including the wind farm and PV station is established. And the characteristic roots of the New England power system integrated with wind energy and PV energy are obtained to study their small-signal stability. In addition, the validity of the theory is verified by the voltage drop of different nodes, which proves that power system integrated with wind-solar renewable energy participating in the frequency regulation can restore the system to the rated frequency in the shortest time and, at the same time, can enhance the robustness of each unit.

1. Introduction

Recently, with the exhaustion of fossil energy and the deterioration of the natural environment, renewable energy has attracted wide attention [1]. Wind energy and solar energy are the most widely used intermittent clean energy, and they are highly complementary in terms of resource and time distribution [2]. If wind and solar energy are integrated to form a wind-solar complementary energy system and participate in the frequency regulation of the power system, the utilization efficiency of intermittent energy can be improved to a certain extent and the global energy shortage can be alleviated [3, 4].

However, the random fluctuation of the output of wind and solar energy causes huge regulatory peak pressure to the power balance of the power system [5, 6]; on the other hand, the power system is disturbed by small-signal all the time during operation [7, 8]. An unstable system is difficult to operate properly in practice [9, 10]. Thus, the analysis of small-signal stability of power system becomes one of the important tasks of power system [11, 12]. Literature [13] establishes a small-signal model of PV generation connected to a weak AC grid. The stability of PV power generation under different power grid strength and control parameters is studied by means of eigenvalue analysis. Literature [14] studies the influence of a large number of wind power generation on small-signal stability and corresponding control strategies to alleviate this negative influence. In [15], the Lyapunov stability criterion is used to analyze the stability research method of the integrated hybrid system. Stability research can be carried out for different renewable energy sources, such as the wind power generation system, photovoltaic system, and micro hydropower system. However, the above analysis regards wind-solar and other renewable energy sources as a perturbation of the power system and does not consider their participation in the frequency regulation of the power system. Therefore, it is not effective in analyzing the stability of the power system in which wind-solar renewable energy participate in the frequency regulation.

Thus, this paper studies the integrated energy system including wind power and PV system with the method of eigenvalue analysis and studies the oscillation modes of the power system when wind and solar power are connected separately, and when the wind farm is connected first and then the PV system is connected. The simulation model of the system is established and the New England power system is used to verify the correctness of the small-signal stability analysis.

The remaining of this paper is organized as follows: Section 2 develops the system modelling. In Section 3, small-signal stability analysis is described. Comprehensive case studies are undertaken in Section 4. The different systems are discussed in Section 5 and Section 6 summarizes the main contributions of the paper.

2. System Modelling

2.1. Multimachine Power System Modelling

The third-order model of the ith generator in a multimachine power system can be expressed by the following formula:where subscript i denotes the variables of the ith machine; is the relative rotor angle; is the generator rotor speed; is the system speed; and are the voltage and transient voltage on the q-axis; is the constant mechanical power input; is the electric power output; is the generator terminal voltage; and are the d-axis and q-axis generator terminal voltages; and are the d-axis synchronous and transient impedances; is the q-axis synchronous impedance; is the rotor inertia; is the d-axis transient short-circuit time constant; and are the d-axis and q-axis generator currents; is the equivalent admittance between the ith and jth nodes; is the susceptance between i and j nodes; is the conductance between i and j nodes; and ufdi and Efdi are the excitation voltage and the initial excitation voltage, respectively.

2.2. System Modelling of DFIG Based Wind Turbine

DFIG is connected to the power system through the voltage source converter, as shown in Figure 1 [16].

The aerodynamic mathematical model of the wind turbine can be described as [16]withwhere ρ is the air density, R denotes the radius of the wind turbine, and means the wind speed. CP (λ, β) is a function of tip-speed-ratio λ and blade pitch angle β representing the power coefficient. A specific wind speed corresponds to a wind turbine rotational speed to obtain CPmax, namely, the maximum power coefficient and therefore tracks the maximum mechanical (wind) power. ωm denotes the wind turbine rotational speed [17].

The 4th-order mathematical model of DFIG can be described aswhere represents the electrical base speed, denotes the synchronous angle speed, and means the rotor angle speed; and denote the equivalent d-axis and q-axis (dq-) internal voltages; and are the dq-stator currents; and represent the dq-stator terminal voltages; and and are the dq-rotor voltages. Lm means mutual inductance.

The pitch angle control system is designed to improve wind energy conversion efficiency and make wind turbine output stable. Its model can be described as follows:where is the reference value of pitch angle; is the inertia time constant of the pitch control system.

The grid-side converter which is directly connected with the power system has the main function of maintaining constant capacitive voltage under the control of the DC regulating system and the function of adjusting the power factor. The DC sides of both converters are supported by a common capacitor. The power equation of the converter can be described as [18]where is the active power of the AC terminal of the machine side converter and is the active power of the AC terminal of the grid-side converter. is the active power of the capacitor tie line; and are the d-q axis components of rotor current respectively; and are the d-q axis components of the system side converter current, respectively; and are the d-q axis components of the system side converter voltage respectively; and are the current and voltage of the DC link in the converter; C is the capacity of the capacitor.

Equation (6) can be rewritten as

2.3. Modelling of PV System

The control structure of the PV system is shown in Figure 2 [19].

According to Kirchhoff’s law, the U-I equation of PV cell can be described aswhere is the short-circuit current; is the open-circuit voltage; is the voltage at maximum power; is the current at maximum power; is the illumination intensity under standard environment, which is 1 kW/m2. is the temperature in the standard environment, which is 25°C. , , , and are, respectively, the correction values of and under different environments. and are temperature compensation coefficients; and is the compensation coefficient of PV irradiation.

In addition, DC/DC converter mainly plays the role of Boost and power transformation, it can be described as

DC link is the intermediate link connecting DC side and AC side, namely, the DC bus capacitance model. According to the capacitance energy and voltage relationship, the DC link model can be described aswhere is the DC side input power of the DC link; is the output power of the DC link inverter side; C is the capacitance value of the DC link capacitance; is the voltage value of the dc link; the is the amount of energy stored on a capacitor.

The wind-solar complementary energy system has three operating states: first, the wind turbine generated independently; second, the PV array independent power generation state; and third, wind-solar complementary power generation. Wind speed, solar radiation, load power consumption, and charging and discharging capacity of the energy storage device all determine the operation state of the wind-solar complementary energy system. Due to the randomness of these factors, the stability of the power system is bound to be affected to some extent. Therefore, it is necessary to analyze the stability of the small-signal of the power system integrated with renewable energy. The control structure diagram of the wind-solar energy system is giving in Figure 3.

3. Small-Signal Stability Analysis

The Lyapunov linearization method is related to the local stability of nonlinear systems. The basic idea is to obtain the local stability of nonlinear systems near their equilibrium operation points from the linear approximation stability property of nonlinear systems [2022].

For the dynamic characteristic differential-algebraic equation of the power system, linearization at the steady-state operating point (x0, y0) can be obtained as follows [2325]:where represents the state variable that describes the dynamic characteristics of the power system in the system of differential equations and represents the operating parameters of the system in algebraic equations. are, respectively, their partial derivatives at steady-state operating point (x0, y0).

Omitting operation parameter y, the following equation can be obtained:with

Matrix A is usually called the state matrix of the system. The stability of the analyzed system at the steady-state operating point (x, y) can be judged by obtaining the eigenvalue of matrix A [2628]:(a)When the real part of all eigenvalues of A is negative, it means that the actual power system can maintain stability when the equilibrium point encounters a small-signal.(b)When at least one real part of all eigenvalues of A is positive, it means that the actual power system will lose stability when it encounters a small-signal at the equilibrium point.(c)When all eigenvalues of A have no positive eigenvalue of the real part, but at least one eigenvalue of the real part is zero, then the linearized system is in A critical stable state, but it cannot be used to judge whether the actual power system is stable at the equilibrium point.(d)A real characteristic root corresponds to a nonoscillating mode. The modes represented by negative real characteristic roots are attenuated, and the greater the absolute value, the faster the corresponding modes decay.(e)Complex characteristic roots always appear as conjugate pairs and can be described as

Complex eigenvalues are always composed of conjugate pairs, which can be described as the negative real part represents the damping oscillation mode [2931]. The positive real part represents the increased oscillation, and the real part of the eigenvalue represents the damping of the system oscillation, while the imaginary part represents the frequency of the system oscillation [32]. The frequency of oscillation can be expressed as [33]

The damping ratio is defined as

It represents the attenuation characteristic of the oscillation amplitude.

4. Case Studies

The proposed methodology is tested on the New England power system, as shown in Figure 4. It consists of 39 buses and 10 generators, and the New York grid connected to the New England power system is represented by the first generator. In addition, detailed system parameters are shown in literature [16]. The proposed methodology has been developed in MATLAB 2017 b environment. In order to analyze the damping characteristics of interconnected systems when wind farm and PV system are connected to the power system, the eigenvalue analysis is carried out for the following four working conditions:(a)Initial system(b)Only wind farms are connected on bus #1 and output 5 MW(c)Only the PV system is connected on bus #2 and output 5 MW(d)Wind farms and the PV system are both connected on bus #1 and bus #2, respectively, and the output is 2.5 MW and 2.5 MW, respectively

Table 1 shows the partial eigenvalues of the system in four cases. It can be seen that when wind farm and PV system are connected separately, their characteristic roots are all far away from the imaginary axis. In particular, after the addition of wind and solar energy, the characteristic root distribution was well improved, which indicates that wind power and PV system independent access system both can significantly improve the stability and, at the same time, are complementary to each other. And Root loci distribution of different conditions is given in Figure 5.

4.1. Bus #3 Voltage Drop

In order to further verify the validity of the characteristic roots above, voltage drop 0.8 p.u. occurred at bus #3 when t = 5 s and recovered after 0.1 s. The corresponding system response is shown in Figure 6. It can be seen that the system that does not involve wind-solar renewable energy in frequency modulation has the worst recovery ability after small-signal, while, with the connection of wind and solar energy, the recovery ability of the system after small-signal is improved. In particular, the system combination of wind and solar energy has the best recovery from small-signal and the ability to adjust the system frequency to near the rated frequency in the shortest amount of time.

4.2. Bus #15 Voltage Drop

In order to verify the recovery ability of the system after receiving small-signal at different positions, voltage drop 0.8 p.u. occurred at bus #15 when t = 5 s and recovered after 0.1 s. And the voltage was restored after 0.1 s. The corresponding system response is shown in Figure 7. It can be found that, with the system combination of wind and PV, the rotor angle difference regulation capacity of generator G1 is significantly improved, its oscillation amplitude is significantly reduced, and it is restored to the rated value in the shortest time. In addition, it has the best regulation ability for active power and reactive power and will adjust the system to the steady-state in the shortest period, so that the system subject to small-signal has the strongest frequency regulation ability.

5. Discussion

5.1. Bus #3 Voltage Drop

In order to further study the positive effect of the energy storage system on the PV station, based on the above case, this paper considers that the energy storage system is configured in the PV station connected to bus #2. In addition, voltage drop 0.8 p.u. occurred at bus #3 when 5 s and recovered after 0.1 s, and its system response is shown in Figure 8. It can be seen that after the PV system is configured with the energy storage system, the stability of small-signal is better. Compared with the initial PV system, it can restore system frequency in a relatively short time.

5.2. Bus #15 Voltage Drop

Consider the following fault: voltage drop 0.8 p.u. occurred at bus #15 with 5 s and recovered after 0.1 s. To verify the stability of the system in case of failure at different locations, its system response is shown in Figure 9. It can be found that the frequency regulation ability of the PV station equipped with the energy storage system is greatly improved, which can well suppress the frequency fluctuation of the power system subjected to small-signal disturbance. In addition, it can help the synchronous generator to recover to the stable state in a short time.

5.3. Comparative Analysis

The integral of absolute error (IAE) of each index in different fault locations is given in Table 2, in which and denotes the reference of variable , respectively. In particular, of the system of combination of wind and PV is merely 61.90%, 74.79%, 78.49%, 82.84%, and 88.27% of that without wind and PV, only wind, only PV, PV station with energy storage system, and PV followed by wind, respectively, acquired in #3 bus voltage drop (bold colour indicates the best results in Table 2).

6. Conclusions

More and more large- and medium-sized renewable energy power stations have been built and connected to the power system, and they account for an increasing proportion of the power system. It affects the stability and damping characteristics of the traditional power system. In this paper, the influence of wind power and photovoltaic energy on the stability of the power system is studied, and the main conclusions are as follows:(a)Based on the calculation of characteristic roots, it is proved that power system integrated with wind and solar energy participating in frequency regulation has better stability.(b)Based on the New England power system, the damping characteristics of the system can be effectively improved and the system can be more stable after the wind-solar renewable energy is incorporated into the power system.(c)Based on the New England power system test, it is verified that the photovoltaic power station can improve its stability to a certain extent after installing the energy storage system. Particularly, acquired by PV with energy storage system is merely 88.30% and 95.40% of that without wind or PV and only PV, respectively, on the case of bus #3 voltage drop.

Data Availability

The data that support the findings of this study are available upon request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the support of the Research and Development Start-Up Foundation of Shantou University (NTF19028) and Research on Regional Power Grid Frequency Regulation with Renewable Energy Participation (YNKJXM20191240).