Next Article in Journal
Chemical and Biochemical Processes of Watershed Ecosystems and Their Impacts on Water Quality
Next Article in Special Issue
Exploring the Regulation Reliability of a Pumped Storage Power Plant in a Wind–Solar Hybrid Power Generation System
Previous Article in Journal
Reanalysis of Soil Moisture Used for Rainfall Thresholds for Rainfall-Induced Landslides: The Italian Case Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low Frequency Oscillations in a Hydroelectric Generating System to the Variability of Wind and Solar Power

1
Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas, Ministry of Education, Northwest A&F University, Yangling 712100, China
2
Institute of Water Resources and Hydropower Research, Northwest A&F University, Yangling 712100, China
3
Powerchina Northwest Engineering Corporation Limited, Xi’an 710065, China
*
Authors to whom correspondence should be addressed.
Water 2021, 13(14), 1978; https://doi.org/10.3390/w13141978
Submission received: 18 June 2021 / Revised: 11 July 2021 / Accepted: 11 July 2021 / Published: 19 July 2021

Abstract

:
The penetration of multiple integrated renewable energies to the power grid are relevant for decision making in energy policy, environment and business. Such an electricity penetration is affected by the intermittent and volatile characteristics of integrated energies, mostly significantly related to the safe and stable electricity production and supply in real world. Here, this paper focuses on the low frequency oscillation analysis of the hydropower generation response to the wind and solar variability. To enable this analysis, a hybrid model of hydropower system integrating with the wind and solar power system is presented. The Nyquist and root-locus stability methods are used to investigate the sensitivity performance of the hydropower governor to the fluctuation of the integrated renewable energies. Additionally, to quantify the risk of the hybrid system, the low frequency oscillation response of hydropower system to wind/solar/hydropower quota and transmission line distance ratio is extensively investigated in this study. The results show that under the case of the wind, solar and hydropower ratio is 40:1:150, the optimal values for maximally reducing hydropower low frequency oscillation are finally determined as kp = 0.8, ki = 0.25 and kd = 0.5. Regarding a certain wind/solar/hydropower quota, it is a promising strategy to increase the solar-load transmission line in order to achieve the safe and stable operation of the hybrid system and a relatively excellent dynamic regulation capacity of the hydropower governor. The model, methods and results implemented in this study are exploited to markedly improve new knowledge applications, policy management, low carbon emissions and investment competitiveness of future energy systems.

1. Introduction

Renewable energy-based electricity generation is on the rise to supplement and replace conventional thermal-based power production. The worldwide newly added renewable electricity capacity accounted for about 50% started in 2015—of which wind and solar power shared almost 70% [1]. This stimulates the integration of multiple renewable energy generation, for example, Denmark, Texas, China, and the United States and California are the typical leaders to push the development of this matter [2]. In conclusion, the penetration of renewable energy integration to the power grid, especially for the wind and solar energies, until recently thought to be a long-distant future scenario, is taking place right now [3,4,5].
However, the volatility and fluctuation of wind speed and solar irradiation in different hours and days present some challenges to the electricity operation and management, which raises a new problem on how to balance the production and consumption in the electricity industry [6,7,8]. In other words, the unpredictable wind and solar power injection increases the uncertainty and unbalance of the power grid; meanwhile, this uncertainty significantly leads to the increase of the reserve and imbalance operation costs and the fast growth of global power price fluctuation [9,10,11,12]. One of the most promising solutions to this challenge is the application of energy storage systems (ESSs) [13,14]. The ESS plays a damper and regulator role in grid-connection of renewable energies, and its storage status is generally summarized in three aspects [15]: (i) decoupling the timing of generation and consumption, (ii) ensuring the peak shaving and switching balance from one power energy to another one, and (iii) regulating electricity end-use quality. To underline these goals, so far the hydropower system, which stores the surplus water energy in peak hours and consumes this energy at peak-load hours, is recognized as an ideal ESS for enhancing the economic and environmental profits, the transmission congestion mitigation, and the reduction of wind/solar curtailment in renewable energy integration [16,17].
The participation of a hydropower system is a proper supplement to tackle the uncertainty problems of wind and solar generation. The explanation for this is the capabilities of a hydropower system in the excess energy storage in noticeable amounts and the fast response to the use of operating reserves through the regulation between upstream and downstream reservoirs [18,19]. Hydropower plays an important role in improving energy utilization efficiency and reducing output fluctuation in a hydro–wind–photovoltaic system. In [20], the proposed three indices have the ability to show the stability and complementarity characteristics of the hybrid system with the time scale varying from second to hour. Wang et al. [21] refined the model of a hydro–wind–photovoltaic system and investigates the complementary mechanism underlying multi-energy systems. Hence, the issue of understanding on how wind and solar power injection affects the operational stability of hydropower system, as one of the practical challenges of the electricity network, is the key focus of this work [22,23,24]. Other merits of the hydropower system in the hybrid renewable energy generation are concluded as follows:
  • Fast responses to frequency and voltage variation under the appropriate controller. Energy planers and regulation entities cheer [25,26].
  • Balancing power price volatility and reducing capital and operating expenses. Storage companies and energy investors cheer [27,28].
  • Hydropower greatly promotes the large-scale applications of renewable energies, considerably reducing greenhouse gas (GHG) emissions and improving public health benefits. Governments and environmental organizations cheer [29,30].
Considering the significant contributions of hydropower system to the electrical generation of renewable energy integration, some of the research papers in this regard are from an optimized scheduling viewpoint which involves in both wind and solar joint energies. In [31,32,33], the optimized algorithms aimed at providing the certain operational observations, hydropower generation schedule and scheduling cost assessment for the wind/solar/hydropower hybrid system. However, while accounting for hydropower electricity, these studies do not focus on the impact of the frequent optimal dispatch on the fatigue property of hydropower components. In [34,35,36,37], the synergetic scheduling model and control strategy of wind/solar/hydropower units are the main focuses, but not involving in the detailed hydraulic–electric–mechanical internal characteristics. Moreover, some papers are from the economic, policy and environment viewpoint. In [38,39,40,41,42], the market value and environmental benefits of the hybrid system is successfully quantified from the perspective of data statistics and the corresponding policy support is also recommended. Nevertheless, a quantitative analysis of the coordinating influences between the hydropower system and wind/solar systems is lacking. This paper aims to fill these gaps.
Based on the above considerations, this paper is devoted to the stability interaction between the hydropower system and the intermittent wind/solar power. The purpose of this study is to analyze the internal stability mechanism in the hybrid system and to quantify the influence of wind/solar power uncertainty on the hydropower low frequency oscillation. Innovation aspects of this paper are:
  • A promising hybrid system of the hydropower generation integrating with wind farm and solar photovoltaic system is established using MATLAB/Simulink, in order to enable the stability analysis. This contributes to the current international pool of the integration modelling knowledge.
  • The sensitivity of hydropower low frequency oscillations to its governor regulation capacity is quantified under the volatility influence of wind and solar energies. The main adopted methods include Nyquist response and root-locus analysis.
  • To understand the stability conditions of the hybrid system, the influence of different wind/solar/hydropower quotas (i.e., W: S: H) and the various transmission line distance ratios on the low frequency oscillation mode of hydropower system are also quantified. The assessment indicators in this part include the three-phase parallel RLC load, the grounding transformer, the three-phase PI section line and the wind-farm transmission line, and the assessment criteria are oscillation frequency and damping ratio.
The remainder of this paper is organized as follows. Section 2 presents a detailed model of the integrated wind/solar/hydropower power system. The low frequency oscillation sensitivity of hydropower system to PID governor regulation capacity is then investigated in Section 3. In Section 4, results for the low frequency oscillation response to renewable-quota and transmission line distance ratio are provided and discussed. Finally, conclusions are drawn in Section 5.

2. Mathematical Model of the Hybrid Power System

In this section, a hybrid system comprising the hydropower system, wind farm and solar photovoltaic system is established using MATLAB/Simulink [43] (see the schematic diagram in Figure 1). The hydropower system is composed of the penstock, governor, hydro-turbine, synchronous generator and excitation sector. The wind farm mainly employs a doubly fed induction generator (DFIG) and wind turbine. Additionally, the solar photovoltaic system focuses on the blocks of the photovoltaic array, photovoltaic inverter and maximum power point tracing (MPPT) controller. The detailed modeling process is illustrated as follows.

2.1. Hydropower System

2.1.1. Penstock

The penstock includes two sectors, i.e., the water diversion pipe and steel pressure pipe, and its dynamic behavior is expressed by elastic water hammer equations, as follows:
Q t g A H L + f Q 2 2 A D = 0 a 2 Q L g A H t = 0
where Q, H, L, f, D, A and a denote the flow of a pipe section, head of a pipe section, distance between the studied pipe section and calculated origin, water flow friction coefficient, pipe diameter, pipe area and water-hammer velocity.
If q = Q / Q r , h = H / H r and l = L / L r , Equation (1) can be rewritten as:
T W q t h l + f 1 ( q 0 2 + 2 q 0 q ) = 0 4 T W T r 2 q l h t = 0
where TW is the water flow inertia time constant, T W = L r Q r g A H r . Tr is the water hammer pressure time constant, T r = 2 L r a . The variable f1 is represented by f 1 = f Q r 2 L r 2 D A 2 g H r , wherein parameters Lr, Qr and Hr are the total pipe length, design flow and design head.
Equation (2) experiences the Laplace transform and then is expressed as:
h ( 1 , S ) = c h ( T r 2 s + f ) h ( 0 , S ) + 2 h W s h ( T r 2 S + f ) q ( 0 , S ) q ( 1 , S ) = 1 2 h W s h ( T r 2 s + f ) h ( 0 , S ) + c h ( T r 2 S + f ) q ( 0 , S )
where the hyperbolic function is s h x = e x e x 2 , c h x = e x + e x 2 .
If the upstream pipe head h ( 1 , S ) and the downstream pipe flow q ( 0 , S ) are measured, thus the upstream pipe flow and the downstream pipe head for the water diversion pipe can be obtained as:
h ( 0 , S ) = 1 c h ( T r 2 S + f ) h ( 1 , S ) 2 h W t h ( T r 2 S + f ) q ( 0 , S ) q ( 1 , S ) = 1 c h ( T r 2 S + f ) q ( 0 , S ) + 1 2 h W t h ( T r 2 S + f ) h ( 1 , S )
Similarly, the Laplace transform expression of the steel pressure pipe is obtained as:
h q = Z 01 q tanh T 01 s
where hq, Z01 and T01 are transient head of hydraulic turbine, standardized value of hydraulic surge impedance of pipeline and elastic time, respectively.
Ignoring the higher order term tanh T 01 s in Equation (5), thus the transfer function from the hydro-turbine flow to hydro-turbine head can be rewritten as:
h q s = Z 01 π 2 T 01 s + T 01 3 s 3 π 2 + 4 T 01 2 s 2 q s

2.1.2. Governor

In this work, the governor adopts the proportional–integral–differential (PID) controller, where the proportional loop aims to shorten response time, the integral loop decreases the steady error, and the differential loop improves the overshoot response. Generally, the output signal of a universal PID controller is expressed as:
Y P I D = Y P + Y I + Y D
and
Y P = K P Δ F Y I = K I S Δ I = K I S Δ F + b p ( Y c Y P I D ) Y D = K D S 1 + T 1 v S Δ F
where KP, KI and KD are the proportional adjustment coefficient, integral adjustment coefficient and differential adjustment coefficient, respectively. △F, Yc, bp and T1v are the relative deviation of the input frequency, relative predefined guide vane opening, permanent interpolation coefficient and differential loop time constant, respectively. If the permanent interpolation coefficient (bp) equals to zero, thus Equation (8) is rewritten as:
G ( S ) = Y P I D ( s ) Δ F ( s ) = K P + K I 1 S + K D S
Based on the above considerations, the block of the PID controller is shown in Figure 2.
One of the important sector in hydro-electric subsystem is the hydraulic cylinder. The hydraulic cylinder sector that transforms the output signal of PID controller to the hydraulic signal can be modeled as Figure 3.

2.1.3. Hydro-Turbine

The hydro-turbine model uses the recommendation expression from IEEE Group, which is expressed as:
p m = A t h q q n l D t y Δ ω
where At, q, h, fp, pm, qnl, △ω, Dt, and y denote the hydro-turbine gain, relative deviation of hydro-turbine flow, relative deviation of hydro-turbine head, relative head loss coefficient in penstock, mechanical output power, relative no-load flow, relative deviation of hydro-turbine rotational speed, damping factor, and relative deviation of guide vane opening, respectively. Thus, the block diagram of the hydro-turbine is performed in Figure 4.

2.1.4. Synchronous Generator

The three-phase synchronous generator is achieved in the DQ rotor framework, and its stator winding uses Y type connection form to link with the internal neutral point. The three-phase synchronous generator works at the generator/electric motor mode, and the electrical part of the machine is represented by a sixth-order state-space model and the mechanical part is the same as in the simplified synchronous machine block [44]. The model takes into account the dynamics of the stator, field, and damper windings [45]. The voltage equation and flux equation of the six-order synchronous generator are:
V d = R s i d + d d t φ d ω R φ q V q = R s i q + d d t φ q ω R φ d V f d = R f d i f d + d d t φ f d V k d = R k d i k d + d d t φ k d V k q 1 = R k q 1 i k q 1 + d d t φ k q 1 V k q 2 = R k q 2 i k q 2 + d d t φ k q 2
and
φ d = L d i d + L m d i f d + i k d φ q = L q i q + L m q i k q φ f d = L f d i f d + L m d i d + i k d φ k d = L k d i k d + L m d i d + i f d φ k q 1 = L k q 1 i k q 1 + L m q i q φ k q 2 = L k q 2 i k q 2 + L m q i q
where V, I, R and φ are the voltage, current, resistance and flux linkage of winding, respectively; ω R is rotor rotation angular velocity; Ld, Lq, Lfd, Lkd, Lkq1 and Lkq2 are self-inductance of winding; and Lmd and Lmq are mutual inductance of winding.

2.1.5. Excitation Sector

The excitation system provides the excitation power to regulate the generator voltage and control the reactive power distribution of the generator in parallel operation states. It has a great influence on the dynamic behaviors of the generator and has the ability to improve the stability of the power system. The model of the excitation system is shown in Equation (13). The block diagram of the excitation system is presented in Figure 5.
The   excitation   system   model : U ˙ R = 1 T A K A U r e f + U P S S U t U F U R U ˙ F = 1 T F K F E ˙ f U F E ˙ f = 1 T L U R S E + K L E f
Thus, in light of the above subsystem block modeling, the global hydropower system (run-of-river hydropower plant) is obtained in Figure 6.

2.2. Wind Farm

2.2.1. Wind Turbine Model

The blades in the wind turbine are used to capture the wind, and the output power of the wind turbine is closely related to the wind energy efficiency (Cp) and the swept area (Ar) of the blade. The operational phenomenon of wind turbines can be described by the aerodynamics theory, and based on this theory, the power captured by blades of the wind turbine can be expressed as:
P W = A r ρ a i r C p v ω 3
Where PW, ρ, Cp, Ar, and vω are the mechanical power output of the wind turbine, air density, wind energy efficiency, swept area of the blade, and wind velocity, respectively. Here, the wind energy efficiency Cp can expressed as the ratio between the power captured by the wind turbine (Pw) and the total power through the wind turbine (Pa), i.e.,
C p = P w P a
As discussed in [36], Equation (15) can be expressed in the form of Equation (16), which is easier to perform the calculation and is written as follows:
C p = 0.22 ( 116 λ a 0.4 β 5 ) e 12.5 / λ a λ a = ( 1 λ + 0.08 β 0.035 β 3 + 1 ) 1
where λa is the tip velocity ratio and β is the pitch angle of the blade. Using all these equations, the wind characteristic curves can be obtained and described as shown in Figure 7.

2.2.2. Mechanical Drive Shaft Model

The mechanical drive shaft is an important component of the wind farm (WF), which coupled the wind turbine with the DFIGURE. In this work, a two-mass drive shaft model, as presented in Figure 8, is used, as it describes the energy transfer between medium/high and low speed shafts while comparing with the one-mass model. The dynamic characteristics of a two-mass mechanical drive shaft model can be expressed in terms of the following equations [37]:
ω ˙ t ω ˙ g T ˙ l s = a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ω t ω g T l s + b 11 b 21 b 31 T a + b 12 b 22 b 32 T m
with
b 11 = 1 J r , b 12 = 0 b 21 = 0 , b 22 = 1 J g b 31 = B l s J r , b 32 = B l s n g J g
and
a 11 = B r J r , a 12 = 0 , a 13 = 1 J r a 21 = 0 , a 22 = B g J g , a 23 = 1 n g J g a 31 = ( K l s B l s B r J r ) , a 32 = 1 n g ( B l s B r J g K l s ) , a 33 = B l s ( J r + n g 2 J g n g 2 J g J r )
where Tls, Tt, Tm, ωt, ωg, Jr, Jg, Br, Bg, Bls and ng are the low speed shaft torque, aerodynamic torque, generator electromagnetic torque, rotor speed, generator speed, rotor inertia, generator inertia, rotor external damping, generator external damping, low speed shaft damping, and ratio constant, respectively.
The combined block diagram representation of the wind turbine and mechanical drive shaft model is shown in Figure 9. As shown in Figure 9, the input of the wind turbine are the blade angle and wind speed, and the output is the aerodynamic torque. The input of the drive train are the aerodynamic torque and generator speed, and the output are the wind turbine speed and torque transmitted through the shaft. Then, the torque is transformed to the torque at the input of the generator.

2.2.3. DFIG Model in dq Frame

The DFIG exhibits a flexible working mode due to its AC voltage excitation system which easily achieves the variable-speed constant frequency (VSCF) generation. The Park transformation can be used to represent the dynamical model of the DFIG in a synchronously rotating dq frame and in this frame, the DFIG model can be expressed as [38]:
v s d = R s i s d + d Φ s d d t ω s Φ s q v s q = R s i s q + d Φ s q d t ω s Φ s d v r d = R r i r d + d Φ r d d t ω r Φ r q v r q = R r i r q + d Φ r q d t ω r Φ r d
where vsd,q, vrd,q, isd,q, ird,q, Φsd,q, Φrd,q, Rs, Rr, ωs, and ωr are the stator voltage in dq frame, rotor voltage in dq frame, stator current in dq frame, rotor current in dq frame, stator flux in dq frame, rotor flux in dq frame, stator resistance, rotor resistance, stator speed, and rotor speed, respectively.
The d- and q-axis are magnetically decoupled, and thus the flux-linkage equations can be written as:
Φ s d = L s i s d + M i r d Φ s q = L s i s q + M i r q Φ r d = L r i r d + M i s d Φ r q = L r i r q + M i s q
where Ls, Lr, and M are the stator inductance, rotor inductance, and mutual inductance, respectively.
Additionally, the electromagnetic torque in dq frame can be expressed as:
T g = 3 2 n p [ Φ s d i s q Φ s q i s d ]
where np is the pole pairs of the DFIG.

2.2.4. PWM Converter Model

The PWM converter can be described through a state-average model, which represents the product of the voltage and duty ratio. Based on the principle of the state-space average model, the three-phase PWM converter in a DFIG-based WF is expressed as:
L f d i a d t = R f i a + e a ( d a v d c + v N 0 ) L f d i b d t = R f i b + e b ( d b v d c + v N 0 ) L f d i c d t = R f i c + e c ( d c v d c + v N 0 )
where daT and daT are the conduction times of switch S1 and switch S2, respectively; with T and da as the switching cycle of the PWM switches and the duty ratio of phase da; iabc represents the currents flowing through three phases; eabc represents the voltages at the grid connection point for three phases; Rf and Lf are the filter resistance and inductance, respectively; vdc is the DC-link voltage at the input of the PWM converter; and vN0 is the neutral voltage at the grid connection point. The subscripts a, b, and c with different variables represent the phase a, phase b, and phase c, respectively.
By applying Kirchhoff’s current law at the input side of the PWM converter where the DC-link capacitor is connected, it can be written as:
i d c = d a i a + d b i b + d c i c
where idc is the Dc current at the input of the PWM converter.
The block diagram representation of the state-average PWM converter model is shown in Figure 10, whose outputs are actually average voltage control signals used for controlling the PWM switches in the converter.

2.3. Solar Photovoltaic System

The solar photovoltaic system involves in three components, i.e., the photovoltaic array, photovoltaic inverter and maximum power point tracing (MPPT) controller. Herein, the photovoltaic array contains several parallel modules, and the current-voltage property of each module is expressed by five operating parameters, i.e., current source (IL, which can utilize solar energy to generate current), diode parameters (I0 and NI), series resistance (rs) and shunt resistance (rsh). These parameters are closely to the variation of temperature and irradiation. Besides, the photovoltaic inverter is an important element to ensure system stability during the grid-connected transient process, which can convert the DC to sinusoidal AC in order to successfully link with power grid. Moreover, maximum power point tracing (MPPT) controller is applied in this study to maximize the solar electricity production. The block diagram if the proposed solar photovoltaic system is performed in Figure 11.
Based on the above proposed wind, solar and hydropower subsystems, a hybrid power system is finally obtained by MATLAB/Simulink [43], as shown in Figure 12.

3. Low Frequency Oscillation Response to Hybrid Regulation

The purpose of this section is to investigate the sensitivity of hydropower low frequency oscillations to regulation capacity of the grid-connected hybrid system. The hydropower governor, as an effective coordinator, has a powerful ability to control and adjust the power fluctuation as the intermittent wind and solar powers inject to the power grid. Based on its merit, the PID governor parameters, involving in the proportional adjustment coefficient (kp), the integral adjustment coefficient (ki) and the differential adjustment coefficient (kd), are selected to analyze their regulation performance on low frequency oscillation of hydropower generation based on the Nyquist response and root-locus analysis.

3.1. Nyquist Response to PID Regulation

3.1.1. Nyquist Profile

The stability characteristic of Nyquist response is described as
R(s) = 1/K+ G(s)H(s)
where R(s), K, G(s) and H(s) represent equivalent function, adjustable gains, feedback path system function and forward path system function, respectively. According to the Nyquist encirclement property, the clockwise number around the point of −1/K equals to the difference between the closed-loop poles in the right half plane and the poles from G(s)H(s) plane. In other words, the stable condition for the closed-loop system is that there is no closed-loop pole in the right half plane, while in this situation the open-loop system can maintain stability. Thus, the Nyquist stability criterion is concluded as: To keep the system stability, the clockwise number around the point of −1/K equals to the negative pole numbers in the G(s)H(s) right half plane. If there is no poles in the G(s)H(s) right half plane, the clockwise number around the point of −1/K should be zero in order to keep the system stability.

3.1.2. Influences of Governor Parameters on Nyquist and Step Responses

The studied domains of the proportional adjustment coefficient (kp), the integral adjustment coefficient (ki) and the differential adjustment coefficient (kd) are determined. In this case study, the wind, solar and hydropower ratio keeps 40 vs. 1 vs. 150 to benefits the identification of leading variate. The Nyquist and step response results are shown in Figure 13, Figure 14 and Figure 15.
From Figure 13a, there are similar behaviors for Nyquist results as the proportional adjustment coefficient (kp) increases from 0.8 to 2.4. The Nyquist trajectory is a transversely zygomorphic closed curve starting at the origin, and it surrounds twice in the clockwise direction and then returns back to the origin. The curvature of Nyquist trajectory increases after it first go through the real-axis, but the curvature gradually decreases after the second real-axis traverse. All Nyquist trajectories intersect at the origin or at a point near the origin, whereas the peak response is continuously enhanced with the decrease of the proportional adjustment coefficient (kp). Moreover, it is clearly observed in Figure 13b, there are some similarities for the step response of the hydropower system under a different proportional adjustment coefficient (kp). For example, the step response within kp = 2.4 first increases to the maximum overshoot (i.e., 0.014) and then gradually stabilized at the value of 0.002. However, notr that there is the smaller proportional adjustment coefficient, the smaller maximum overshoot. This reveals that the hydropower block has an excellent regulation capacity to the low frequency oscillation of the hybrid system within kp = 0.8 compared with the situation within kp = 2.4
The Nyquist response in Figure 14a experiences two clockwise rotations to return to the origin, meaning that the Nyquist trajectory is a transversely zygomorphic closed curve for the variation range of the integral adjustment coefficient ki = [0.25, 1.25]. These curves only intersect at the origin, although the curvature of the first circle is greater than that of the second circle for a specific value of ki. Additionally, regarding all different values of ki, their second circles move right along with the real-axis and the corresponding peak response improved from ki = 0.25 to ki = 1.25. As performed in Figure 14b, all step responses have a decaying trend under different ki conditions, differently, the peak value and final value of the step response increase gradually with the increase of integral adjustment coefficient (ki). For instance, the peak value and final value for ki = 1.25 are 0.0115 and 0.0081 in comparison with the relatively smaller values of 0.0098 and 0.001 for ki = 0.25. However, fortunately, all overshoots for various parameter settings of ki are less than the value of 10−3, thus the hydropower system can keep its stability during the selected variation domain of ki.
Figure 15 shows that there is a similar change rule for both Nyquist and step responses for the change of the differential adjustment coefficient between kd = 0.5 and kd = 1.5. For a certain value of kd, its Nyquist trajectory starts and also finally converges at the origin by turning clockwise twice; meanwhile, the curvature of the first circle is greater than that of the second circle. The peak response of various Nyquist trajectories increase as the differential adjustment coefficient rises from kd = 0.5 to kd = 1.5. Additionally, the step responses of the hydropower system perform decaying trends under various kd conditions, and these response curves finally converge to a certain value, around 0.0017. The peak values for different values of kd are extremely close to the mean value of 0.0098, although the peak value increases slightly with the change of the differential adjustment coefficient from kd = 0.5 to kd = 1.5. This implies that the studied value domain of the differential adjustment coefficient between (kd) has a lower sensitivity to the low frequency oscillation of hydropower system.

3.1.3. Root-Locus Profile

To further study the pole and zero point characteristics of the low frequency oscillation stability of hydropower system for the solar/wind/hydropower hybrid system, the root-locus analysis is extensively conducted in this section. Herein, it is emphasized on that the zero point in the root-locus trajectory is the pole from the open-loop system while the pole in the root-locus trajectory equals to the closed-loop pole in this work. The detailed root-locus mechanism can be described as follows:
If the step response of the system is expressed as:
C t = 1 C s = A 0 + i = 1 n A i e s i t
where C(s), A0 and Ai denote the Laplace transformation of the step response, and residues corresponding to zero and pole, respectively. Hence, the general stability rule influenced by the pole and zero point in the root-locus trajectory is summarized two aspects: (1) The stability condition needs all closed-loop poles be located in the left half plane, meaning that all root-locus trajectories are in the left half plane, and (2) the rapidness of system response depends on the distance between the closed-loop pole (si = σ + ) and the imaginary-axis. If such a distance is shorter, there exists a rapider decaying speed of the esit and thus results in an excellent rapidness of system response.

3.1.4. Influences of Governor Parameters on Root-Locus Response

To clearly perform the PID parameters regulation on the root-locus response, the control variate method is adopted. That is, two of the PID parameters remain unchanged, but only one of them changes gradually. Besides, the study is conducted based on the wind, solar and hydropower quota 40:1:150. The root-locus results are shown in Figure 16, Figure 17 and Figure 18.
Both zero points and poles (Figure 16a) lie in the real-axis under the condition of kp = 1.6, ki = 0.75 and kd = 0.5, and the zero point coordinates are (−0.0925, 0) and (−1.69, 0) as well as the pole coordinates are (−1, 0), (−1.05, 0) and (−15, 0). There are three root-locus responses: (i) The first root-locus trajectory begins at the pole (−1, 0) and ends at the zero point (−0.0925, 0). (ii) The second root-locus trajectory starts from the pole (−1.05, 0) and is extended towards the zero point (−1.69, 0), and (iii) the third root-locus trajectory stems from the pole (−15, 0) and goes to the negative infinity. All root-locus trajectories are located in the left half plane, and thus the hydropower system is stable in this situation. With the increase of the differential adjustment coefficient from kd = 0.8 to kd = 2.4, the pole coordinate (Figure 16b) remains unchanged, but the zero point is gradually closer to the imaginary-axis. This indicates that there is an excellent rapidness response of the hydropower system to the disturbance of low frequency oscillations within a smaller kp.
From Figure 17a, it is found that the all poles and zero points are in the real-axis under the condition of kp = 1.6, ki = 0.75 and kd = 1. The zero point coordinates include (−0.255, 0) and (−1.84, 0), and the pole coordinates involve in (−1, 0), (−1.1, 0) and (−15, 0). There are three root-locus trajectories: the first trajectory is between the pole (−1, 0) and zero point (−0.255, 0), the second trajectory is from the pole (−1.1, 0) to the zero point (−1.84, 0), and the last trajectory starts from the pole (−15, 0) and links to the negative infinity. The hydropower system operates safely and stably since these root-locus trajectories are in the left halt plane. Figure 17b illustrates the regulation capacity of the integral adjustment coefficient from ki = 0.25 to ki = 1.25. With the decrease of parameter ki, the pole coordinate keeps constant while the zero point moves closely to the imaginary-axis. This indicates that the hydropower system has a better rapidness as the parameter ki increases to 1.25.
Similarly, the poles and zero points are in the real-axis for kp = 1.6, ki = 0.75 and kd = 1 in Figure 18a. The zero points coordinates are (−0.255, 0) and (−1.84, 0), while the pole coordinates are (−1, 0), (−1.1, 0) and (−15, 0). The first root-locus trajectory is from the pole (−1, 0) to the zero point (−0.255, 0), the second root-locus trajectory is from the pole (−1.1, 0) to the zero point (−1.84, 0), and the last root-locus trajectory experiences from the pole (−15, 0) to negative infinity. All the three trajectories lie in the left half plane, resulting in a steady operational condition of the hydropower system responding to the low frequency oscillations. Moreover, as shown in Figure 18b, the zero point has an opposite trend with the pole with the increase of differential adjustment coefficient (kd), where the zero point moves closely to the imaginary-axis, but the pole is far away from the imaginary-axis. This reveals that there is a rapid response of the hydropower system to the low frequency oscillation within kd = 1 in comparison to the condition of kd = 1.6.

4. Low Frequency Oscillation Response to Renewable-Quota

This section focuses on investigating the influence of different wind/solar/hydropower quota (i.e., W:S:H) on the low frequency oscillation mode and system stability. Simultaneously, the line distance ratio between the solar-load transmission and the wind-load transmission (i.e., SL:WL) is also considered. Based on this, the wind/solar/hydropower quota include 20:1:150, 30:1:150 and 40:1:150 (the combination is ratio of installed capacity), and the line distance are determined as 1:1, 2:1 and 3:1. These renewable quotas are selected in this paper to study the evolution law between the increase of installed capacity of wind power and the low-frequency oscillation of power grid under different renewable quotas. Similarly, the transmission line length ratios are selected in this paper to study the influence of the solar-load line on the low frequency oscillation in the hydropower-based hybrid renewable energy system. The low frequency oscillation mainly be represented by the three-phase parallel RLC load, the grounding transformer, the three-phase PI section line and the wind-farm transmission line. The results are demonstrated in Table 1.
With respect to the line distance ratio between the solar-load transmission and the wind-load transmission (i.e., SL:WL) is 1:1, the evaluation results (Table 1) show that the oscillation frequency of the three-phase parallel RLC load and the wind-farm transmission line roughly keeps unchanged as the wind/solar/hydropower quota (i.e., W:S:H) changes from 20:1:150, 30:1:150 to 40:1:150. Conversely, the oscillation frequency of the grounding transformer and the three-phase PI section line continuously increases in this situation, for instance, the oscillation frequency of grounding transformer is 2.01 vs. 2.63 vs. 3.17, and the oscillation frequency of three-phase PI section line is 7.26 × 102 vs. 7.48 × 102 vs. 7.71 × 102. Besides, the damping ratios for these four stability indicators are constant. This suggests that the appropriate increase of the wind farm installed capacity is beneficial for the electricity production without system stability defects. As regards the condition of SL:WL = 2:1, the oscillation frequency of the three-phase parallel RLC load and the wind-farm transmission line is unchanged while the oscillation frequency of the grounding transformer and the three-phase PI section line obviously increases with the capacity-increasing improvement of wind farm. The oscillation damping ratio of the three-phase parallel RLC load and the grounding transformer has no variation in contrast to the decline damping ratio of the three-phase PI section line and the wind-farm transmission line. It is indicated that the wind farm installed capacity is unable to design too large, otherwise the stability problem is likely to occur under the condition of SL:WL = 2:1. Regarding the condition of SL:WL = 3:1, there is no variation for the oscillation frequency of the three-phase parallel RLC load, the three-phase PI section line and the wind-farm transmission line, but the oscillation frequency of the grounding transformer rises. Differently, the oscillation damping ratio of the three-phase parallel RLC load, the grounding transformer and the wind-farm transmission line totally remains unchanged, whereas the oscillation damping ratio of the grounding transformer first declines and then rises. Thus, it is possible to increase the installed capacity of the wind farm since the hydropower generation has a low sensitivity to the low frequency oscillations in this case.
In addition, as the line distance ratio between the solar-load transmission and the wind-load transmission increases from 1:1 to 3:1, the oscillation frequency and damping ratio of the three-phase parallel RLC load and the grounding transformer are generally unchanged; however, such indicators of the three-phase PI section line has an increase trend and then gradually declines. Additionally, the oscillation frequency of the wind-farm transmission line decreases, although its damping ratio decreases first and then maintains an increase trend. In other words, when the renewable-quota is determined, the system stability and its dynamic behaviors are able to be enhanced by increasing the length of the solar-load transmission line.

5. Conclusions

Renewable energies, especially for wind and solar energies, are major electricity generation strategies to address the constraints of resources and the environment for energy development. However, the intermittency, uncontrollability and randomness demerits of wind/solar generations pose large challenges for the stable integration into the power grid. Based on this, a hybrid system of the hydropower integrating with wind and solar sectors is presented in order to investigate the low frequency oscillation stability of hydropower systems responding to the fluctuation of wind and solar power injections. The main adopted approaches are Nyquist response and root-locus analysis to evaluate the regulation capacity of governor parameters (i.e., the proportional adjustment coefficient kp, the integral adjustment coefficient ki, and the differential adjustment coefficient kd) to the system stability profile. Furthermore, the influence of different renewable-quota and transmission line distance ratios on the low frequency oscillation mode of the hydropower system is quantified. The major conclusions can be drawn as follows.
  • Under the case where the wind, solar and hydropower ratio is 40:1:150, it is interesting that the smaller the governor parameters (kp, ki, and kd), the smaller the Nyquist overshoot and step fluctuation. Herein, the studied domains for kp, ki, and kd are [0.8, 2.4], [0.25, 1.25] and [0.5, 1.5], and thus the optimal values for maximally reducing hydropower low frequency oscillation are finally determined as kp = 0.8, ki = 0.25 and kd = 0.5.
  • The wind/solar/hydropower hybrid system keeps global stability in the studied governor parameter domains since the Nyquist and root-locus low frequency oscillation responses meet the relevant stability criteria, i.e., the clockwise number around the point of −1/K equals to the negative pole numbers in the G(s)H(s) right half plane, as well as all root-locus trajectories are in the left half plane. Despite this merit, the overshoot problem is expected to arouse great attention and discussion to reduce the fatigue damage of hydropower components.
  • Aiming at different wind/solar/hydropower quotas (i.e., 20: 1: 150, 30: 1: 150, and 40: 1: 150), the four quantified indicators (i.e., the three-phase parallel RLC load, the grounding transformer, the three-phase PI section line, and the wind-farm transmission line) show that it is beneficial to increase the wind farm installed capacity to maximize the electricity production without system stability defects under the solar-load and wind-load line ratios for 1:1 and 3:1 excepting for 2:1. This is contributed by the smaller quantified values of oscillation frequency and damping ratio.
  • Regarding a certain wind/solar/hydropower quota, it is a promising strategy to increase the solar-load transmission line in order to achieve the safe and stable operation of the hybrid system and a relatively excellent dynamic regulation capacity of the hydropower governor.
Hydropower systems, such as the vital energy storage system and power balance controller, play a key role towards the successful achievement of the world’s long-term clean energy goals. For this, the results provide important guidance for realizing the maximization of economic and environmental benefits in a hybrid power grid. The proposed solutions evaluate the safety situation of the hybrid system in quantitative terms, which is cheered by the energy stakeholders of governments, energy planers, regulation entities, storage companies and environmental organizations.

Author Contributions

Conceptualization, B.X. and L.L.; methodology, B.X.; software, L.L.; validation, Z.Z., S.X. and H.L.; formal analysis, J.Z.; data curation, W.J.; writing—original draft preparation, D.C.; writing—review and editing, W.J.; visualization, D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Chinese Universities Scientific Fund, grant number 2452020210; scientific research foundation of the Natural Science Foundation of Shaanxi Province of China, grant number 2019JLP-24; Shaanxi Science and Technology Innovation Team, and Water Conservancy Science and Technology Program of Shaanxi Province, grant number 2018slkj-9.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hirth, L. The benefits of flexibility: The value of wind energy with hydropower. Appl. Energy 2016, 181, 210–223. [Google Scholar] [CrossRef]
  2. Lu, X.; McElroy, M.B.; Peng, W.; Liu, S.; Nielsen, C.P.; Wang, H. Challenges faced by China compared with the US in developing wind power. Nat. Energy 2016, 1, 16061. [Google Scholar] [CrossRef]
  3. Hart, E.K.; Stoutenburg, E.D.; Jacobson, M.Z. The Potential of Intermittent Renewables to Meet Electric Power Demand: Current Methods and Emerging Analytical Techniques. Proc. IEEE 2012, 100, 322–334. [Google Scholar] [CrossRef]
  4. Schaber, K.; Steinke, F.; Muehlich, P.; Hamacher, T. Parametric study of variable renewable energy integration in Europe: Advantages and costs of transmission grid extensions. Energy Policy 2012, 42, 498–508. [Google Scholar] [CrossRef] [Green Version]
  5. Kaur, A.; Nonnenmacher, L.; Coimbra, C.F.M. Net load forecasting for high renewable energy penetration grids. Energy 2016, 114, 1073–1084. [Google Scholar] [CrossRef]
  6. Talaat, M.; Farahat, M.A.; Elkholy, M.H. Renewable power integration: Experimental and simulation study to investigate the ability of integrating wave, solar and wind energies. Energy 2019, 170, 668–682. [Google Scholar] [CrossRef]
  7. Zhang, H.; Cao, Y.; Zhang, Y.; Terzija, V. Quantitative synergy assessment of regional wind-solar energy resources based on MERRA reanalysis data. Appl. Energy 2018, 216, 172–182. [Google Scholar] [CrossRef]
  8. Cao, Y.; Zhang, Y.; Zhang, H.; Shi, X.; Terzija, V. Probabilistic Optimal PV Capacity Planning for Wind Farm Expansion Based on NASA Data. IEEE Trans. Sustain. Energy 2017, 8, 1291–1300. [Google Scholar] [CrossRef]
  9. Zini, G.; Rosa, A.D. Hydrogen systems for large-scale photovoltaic plants: Simulation with forecast and real production data. Int. J. Hydrog. Energy 2014, 39, 107–118. [Google Scholar] [CrossRef]
  10. Wang, D.; Qiu, J.; Reedman, L.; Meng, K.; Lai, L.L. Two-stage energy management for networked microgrids with high renewable penetration. Appl. Energy 2018, 226, 39–48. [Google Scholar] [CrossRef]
  11. Vachirasricirikul, S.; Ngamroo, I. Robust LFC in a Smart Grid with Wind Power Penetration by Coordinated V2G Control and Frequency Controller. IEEE Trans. Smart Grid 2014, 5, 371–380. [Google Scholar] [CrossRef]
  12. Zhao, B.; Wang, X.; Lin, D.; Calvin, M.M.; Morgan, J.C.; Qin, R.; Wang, C. Energy Management of Multiple Microgrids Based on a System of Systems Architecture. IEEE Trans. Power Syst. 2018, 33, 6410–6421. [Google Scholar] [CrossRef]
  13. Etherden, N.; Bollen, M.H.J. Dimensioning of Energy Storage for Increased Integration of Wind Power. IEEE Trans. Sustain. Energy 2013, 4, 546–553. [Google Scholar] [CrossRef]
  14. Hill, C.A.; Such, M.C.; Chen, D.; Gonzalez, J.; Grady, W.M. Battery Energy Storage for Enabling Integration of Distributed Solar Power Generation. IEEE Trans. Smart Grid 2012, 3, 850–857. [Google Scholar] [CrossRef]
  15. Notton, G.; Mistrushi, D.; Stoyanov, L.; Berber, P. Operation of a photovoltaic-wind plant with a hydro pumping-storage for electricity peak-shaving in an island context. Sol. Energy 2017, 157, 20–34. [Google Scholar] [CrossRef]
  16. Corsini, A.; Rispoli, F.; Gamberale, M.; Tortora, E. Assessment of H-2- and H2O-based renewable energy-buffering systems in minor islands. Renew. Energy 2009, 34, 279–288. [Google Scholar] [CrossRef]
  17. Sternberg, A.; Bardow, A. Power-to-What?—Environmental assessment of energy storage systems. Energy Environ. Sci. 2015, 8, 389–400. [Google Scholar] [CrossRef]
  18. Fitzgerald, N.; Arantegui, R.L.; McKeogh, E.; Leahy, P. A GIS-based model to calculate the potential for transforming conventional hydropower schemes and non-hydro reservoirs to pumped hydropower schemes. Energy 2012, 41, 483–490. [Google Scholar] [CrossRef]
  19. Su, Y.; Kern, J.D.; Characklis, G.W. The impact of wind power growth and hydrological uncertainty on financial losses from oversupply events in hydropower-dominated systems. Appl. Energy 15 2017, 194, 172–183. [Google Scholar] [CrossRef]
  20. Xiong, H.; Xu, B.; Kheav, K.; Luo, X.; Zhang, X.; Patelli, E.; Guo, P.; Chen, D. Multiscale power fluctuation evaluation of a hydro-wind-photovoltaic system. Renew. Energy 2021, 175, 153–166. [Google Scholar] [CrossRef]
  21. Wang, X.; Mei, Y.; Kong, Y.; Lin, Y.; Wang, H. Improved multi-objective model and analysis of the coordinated operation of a hydro-wind-photovoltaic system. Energy 2017, 134, 813–839. [Google Scholar] [CrossRef]
  22. Xu, B.; Chen, D.; Venkateshkumar, M.; Xiao, Y.; Yue, Y.; Xing, Y.; Li, P. Modeling a pumped storage hydropower integrated to a hybrid power system with solar-wind power and its stability analysis. Appl. Energy 2019, 248, 446–462. [Google Scholar] [CrossRef]
  23. Martinez-Lucas, G.; Ignacio Sarasua, J.; Angel Sanchez-Fernandez, J.; Roman Wilhelmi, J. Frequency control support, of a wind-solar isolated system by a hydropower plant with long tail-race tunnel. Renew. Energy 2016, 90, 362–376. [Google Scholar] [CrossRef] [Green Version]
  24. Goyal, R.; Gandhi, B.K. Review of hydrodynamics instabilities in Francis turbine during off-design and transient operations. Renew. Energy 2018, 116, 697–709. [Google Scholar] [CrossRef]
  25. Zhang, S.; Zhang, Y. Characteristic Analysis and Calculation of Frequencies of Voltages in Out-of-Step Oscillation Power System and a Frequency-Based Out-of-Step Protection. IEEE Trans. Power Syst. 2019, 34, 205–214. [Google Scholar] [CrossRef]
  26. Chen, C.-I.; Chang, G.W. An Efficient Prony-Based Solution Procedure for Tracking of Power System Voltage Variations. IEEE Trans. Ind. Electron. 2013, 60, 2681–2688. [Google Scholar] [CrossRef]
  27. Rintamaki, T.; Siddiqui, A.S.; Salo, A. Does renewable energy generation decrease the volatility of electricity prices? An analysis of Denmark and Germany. Energy Econ. 2017, 62, 270–282. [Google Scholar]
  28. Roozbehani, M.; Dahleh, M.A.; Mitter, S.K. Volatility of Power Grids under Real-Time Pricing. IEEE Trans. Power Syst. 2012, 27, 1926–1940. [Google Scholar] [CrossRef] [Green Version]
  29. El-Fadel, M.; Chedid, R.; Zeinati, M.; Hmaidan, W. Mitigating energy-related GHG emissions through renewable energy. Renew. Energy 2003, 28, 1257–1276. [Google Scholar] [CrossRef]
  30. Tripathi, L.; Mishra, A.K.; Dubey, A.K.; Tripathi, C.B.; Baredar, P. Renewable energy: An overview on its contribution in current energy scenario of India. Renew. Sustain. Energy Rev. 2016, 60, 226–233. [Google Scholar] [CrossRef]
  31. Perez-Diaz, J.I.; Jimenez, J. Contribution of a pumped-storage hydropower plant to reduce the scheduling costs of an isolated power system with high wind power penetration. Energy 2016, 109, 92–104. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, B.; Guo, W.; Yang, J. Analytical solutions for determining extreme water levels in surge tank of hydropower station under combined operating conditions. Commun. Nonlinear Sci. Numer. Simul. 2017, 47, 394–406. [Google Scholar] [CrossRef]
  33. Lai, X.; Li, C.; Guo Wg Xu, Y.; Li, Y. Stability and dynamic characteristics of the nonlinear coupling system of hydropower station and power gird. Commun. Nonlinear Sci. Numer. Simul. 2019, 79, 104919. [Google Scholar] [CrossRef]
  34. Su, C.; Cheng, C.; Wang, P.; Shen, J.; Wu, X. Optimization model for long-distance integrated transmission of wind farms and pumped-storage hydropower plants. Appl. Energy 2019, 242, 285–293. [Google Scholar] [CrossRef]
  35. Helseth, A.; Gjelsvik, A.; Mo, B.; Linnet, U. A model for optimal scheduling of hydro thermal systems including pumped-storage and wind power. IET Gener. Transm. Distrib. 2013, 7, 1426–1434. [Google Scholar] [CrossRef]
  36. Wang, X.; Chen, L.; Chen, Q.; Mei, Y.; Wang, H. Model and Analysis of Integrating Wind and PV Power in Remote and Core Areas with Small Hydropower and Pumped Hydropower Storage. Energies 2018, 11, 3459. [Google Scholar] [CrossRef] [Green Version]
  37. Xu, B.B.; Chen, D.Y.; Zhang, H.; Wang, F.F.; Zhang, X.G.; Wu, Y.H. Hamiltonian model and dynamic analyses for a hydro-turbine governing system with fractional item and time-lag. Commun. Nonlinear Sci. Numer. Simul. 2017, 47, 35–37. [Google Scholar] [CrossRef] [Green Version]
  38. Arena, C.; Genco, M.; Lombardo, A.; Meli, I.; Mazzola, M.R. A Cost-Benefit Based, Parametric Procedure to Screen Existing Irrigation and Municipal Supply Reservoirs for Wind Energy Storage. Water 2018, 10, 1813. [Google Scholar] [CrossRef] [Green Version]
  39. Ouyang, W.; Hao, F.H.; Zhao, C.; Lin, C. Vegetation response to 30 years hydropower cascade exploitation in upper stream of Yellow River. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 1928–1941. [Google Scholar] [CrossRef]
  40. Karhinen, S.; Huuki, H. Private and social benefits of a pumped hydro energy storage with increasing amount of wind power. Energy Econ. 2019, 81, 942–959. [Google Scholar] [CrossRef]
  41. Fahimi, A.; Upham, P. The renewable energy sector in Afghanistan: Policy and potential. Wiley Interdiscip. Rev. Energy Environ. 2018, 7, e280. [Google Scholar] [CrossRef]
  42. Jaime, R.H.; Mahmut, R.; Tang, H.; Erdal, K. Nonlinear control of fixed-wing UAVs in presence of stochastic winds. Commun. Nonlinear Sci. Numer. Simul. 2016, 33, 57–69. [Google Scholar]
  43. Emaminaeini, A.; Franklin, G.F.; Powell, J.D.; MathWorks, Inc. Feedback Control of Dynamic Systems/MATLAB & Simulink Student Version 2012a. J. De Gynécologie Obs. Et Biol. De La Reprod. 2012, 17, 67–74. [Google Scholar]
  44. Krause, P.C.; Wasynczuk, O.; Sudhoff, S.D. Analysis of Electric Machinery, Analysis of Electric Machinery and Drive; Wiley-IEEE Press: Piscataway, NJ, USA, 1986. [Google Scholar]
  45. Kamwa, I.; Pilote, M. Experience with computer-aided graphical analysis of sudden-short-circuit oscillograms of large synchronous machines. IEEE Trans. Energy Conver. 1995, 10, 407–414. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of wind/solar/hydropower integrated system. Images of wind turbine, PV generator, and electric load from the ref. [20].
Figure 1. Schematic diagram of wind/solar/hydropower integrated system. Images of wind turbine, PV generator, and electric load from the ref. [20].
Water 13 01978 g001
Figure 2. Block diagram of the PID controller in hydropower system. Parameters kp, ki and kd are the proportional adjustment coefficient, integral adjustment coefficient and differential adjustment coefficient, respectively.
Figure 2. Block diagram of the PID controller in hydropower system. Parameters kp, ki and kd are the proportional adjustment coefficient, integral adjustment coefficient and differential adjustment coefficient, respectively.
Water 13 01978 g002
Figure 3. Block diagram of the hydraulic cylinder sector.
Figure 3. Block diagram of the hydraulic cylinder sector.
Water 13 01978 g003
Figure 4. Block diagram of the hydro-turbine. Parameters At, q, h, fp, pm, qnl, △ω, Dt, y hfc and hq denote the hydro-turbine gain, relative deviation of hydro-turbine flow, relative deviation of hydro-turbine head, relative head loss coefficient in penstock, mechanical output power, relative no-load flow, relative deviation of hydro-turbine rotational speed, damping factor, relative deviation of guide vane opening, relative head friction loss and hydro-turbine head change caused by flow, respectively.
Figure 4. Block diagram of the hydro-turbine. Parameters At, q, h, fp, pm, qnl, △ω, Dt, y hfc and hq denote the hydro-turbine gain, relative deviation of hydro-turbine flow, relative deviation of hydro-turbine head, relative head loss coefficient in penstock, mechanical output power, relative no-load flow, relative deviation of hydro-turbine rotational speed, damping factor, relative deviation of guide vane opening, relative head friction loss and hydro-turbine head change caused by flow, respectively.
Water 13 01978 g004
Figure 5. The block diagram of the excitation system.
Figure 5. The block diagram of the excitation system.
Water 13 01978 g005
Figure 6. Block of an integrated hydropower system.
Figure 6. Block of an integrated hydropower system.
Water 13 01978 g006
Figure 7. Wind turbine characteristic curves.
Figure 7. Wind turbine characteristic curves.
Water 13 01978 g007
Figure 8. Diagram of a two-mass drive shaft model of a mechanical drive shaft model. In this figure, the parameters Tt, Ths, Tls, Tm, ωt, ωg, ωls, Jr, Jg, Br, Bg, Bls, Kls, and ng represent the aerodynamic torque, high speed shaft torque, low speed shaft torque, generator electromagnetic torque, rotor speed, generator speed, low speed shaft speed, rotor inertia, generator inertia, rotor external damping, generator external damping, low speed shaft damping, low speed shaft stiffness, and ratio constant, respectively.
Figure 8. Diagram of a two-mass drive shaft model of a mechanical drive shaft model. In this figure, the parameters Tt, Ths, Tls, Tm, ωt, ωg, ωls, Jr, Jg, Br, Bg, Bls, Kls, and ng represent the aerodynamic torque, high speed shaft torque, low speed shaft torque, generator electromagnetic torque, rotor speed, generator speed, low speed shaft speed, rotor inertia, generator inertia, rotor external damping, generator external damping, low speed shaft damping, low speed shaft stiffness, and ratio constant, respectively.
Water 13 01978 g008
Figure 9. Block diagram of the wind turbine and mechanical drive shaft model. The symbols in this figure, W_wt, Pitch_deg, WS, Tt, T_wt, ωg, T_hs and Tm denote the wind turbine speed, blade pitch angle, wind speed, torque of the mechanical drive shaft, wind turbine torque, generator speed, torque transmitted through the shaft, and torque at the input of the generator, respectively.
Figure 9. Block diagram of the wind turbine and mechanical drive shaft model. The symbols in this figure, W_wt, Pitch_deg, WS, Tt, T_wt, ωg, T_hs and Tm denote the wind turbine speed, blade pitch angle, wind speed, torque of the mechanical drive shaft, wind turbine torque, generator speed, torque transmitted through the shaft, and torque at the input of the generator, respectively.
Water 13 01978 g009
Figure 10. Block diagram of a state-average PWM converter model. In this figure, the parameters Uctrl_grid_conv, Uctrl_rotor_conv, Uavg_grid_conv, Uavg_rotor_conv, Iabc_grid_conv_pu, Iabc_rotor_pu, Vab_gc, Vbc_gc, Vab_rc, Vbc_rc, Vdc, and Asyncmac_sig represent the grid-side converter voltage control signal, rotor-side converter voltage control signal, grid-side average converter voltage signal, rotor-side average converter voltage signal, grid converter current signal, rotor converter current signal, grid-side voltage between phases a and b, grid-side voltage between phases b and c, rotor-side voltage between phases a and b, rotor-side voltage between phases b and c, DC-link voltage, and the asynchronous machine signal, respectively.
Figure 10. Block diagram of a state-average PWM converter model. In this figure, the parameters Uctrl_grid_conv, Uctrl_rotor_conv, Uavg_grid_conv, Uavg_rotor_conv, Iabc_grid_conv_pu, Iabc_rotor_pu, Vab_gc, Vbc_gc, Vab_rc, Vbc_rc, Vdc, and Asyncmac_sig represent the grid-side converter voltage control signal, rotor-side converter voltage control signal, grid-side average converter voltage signal, rotor-side average converter voltage signal, grid converter current signal, rotor converter current signal, grid-side voltage between phases a and b, grid-side voltage between phases b and c, rotor-side voltage between phases a and b, rotor-side voltage between phases b and c, DC-link voltage, and the asynchronous machine signal, respectively.
Water 13 01978 g010
Figure 11. Block diagram of the solar photovoltaic system applied in this study.
Figure 11. Block diagram of the solar photovoltaic system applied in this study.
Water 13 01978 g011
Figure 12. Block diagram of the wind/solar/hydropower integrated power system.
Figure 12. Block diagram of the wind/solar/hydropower integrated power system.
Water 13 01978 g012
Figure 13. Nyquist and step responses to the variation of the proportional adjustment coefficient, kp = [0.8, 2.4]. (a) Nyquist response, and (b) step response.
Figure 13. Nyquist and step responses to the variation of the proportional adjustment coefficient, kp = [0.8, 2.4]. (a) Nyquist response, and (b) step response.
Water 13 01978 g013
Figure 14. Nyquist and step responses to the variation of the integral adjustment coefficient, ki = [0.25, 1.25]. (a) Nyquist response, and (b) step response.
Figure 14. Nyquist and step responses to the variation of the integral adjustment coefficient, ki = [0.25, 1.25]. (a) Nyquist response, and (b) step response.
Water 13 01978 g014
Figure 15. Nyquist and step responses to the variation of the differential adjustment coefficient, kd = [0.5, 1.5]. (a) Nyquist response, and (b) step response.
Figure 15. Nyquist and step responses to the variation of the differential adjustment coefficient, kd = [0.5, 1.5]. (a) Nyquist response, and (b) step response.
Water 13 01978 g015
Figure 16. Root-locus response to the variation of the proportional adjustment coefficient, kp = [0.8, 2.4]. (a) Root-locus response at kp = 1.6, ki = 0.75 and kd = 0.5, and (b) root-locus response with the variation of kp = 0.8, 1.6 and 2.4.
Figure 16. Root-locus response to the variation of the proportional adjustment coefficient, kp = [0.8, 2.4]. (a) Root-locus response at kp = 1.6, ki = 0.75 and kd = 0.5, and (b) root-locus response with the variation of kp = 0.8, 1.6 and 2.4.
Water 13 01978 g016
Figure 17. Root-locus response to the variation of the integral adjustment coefficient, ki = [0.25, 1.25]. (a) Root-locus response at kp = 1.6, ki = 0.75 and kd = 1, and (b) root-locus response with the variation of ki = 0.25, 0.75 and 1.25.
Figure 17. Root-locus response to the variation of the integral adjustment coefficient, ki = [0.25, 1.25]. (a) Root-locus response at kp = 1.6, ki = 0.75 and kd = 1, and (b) root-locus response with the variation of ki = 0.25, 0.75 and 1.25.
Water 13 01978 g017
Figure 18. Root-locus response to the variation of the differential adjustment coefficient, kd = [0.5, 1.5]. (a) Root-locus response at kp = 1.6, ki = 0.75 and kd = 1, and (b) Root-locus response with the variation of kd = 0.5, 1 and 1.5.
Figure 18. Root-locus response to the variation of the differential adjustment coefficient, kd = [0.5, 1.5]. (a) Root-locus response at kp = 1.6, ki = 0.75 and kd = 1, and (b) Root-locus response with the variation of kd = 0.5, 1 and 1.5.
Water 13 01978 g018
Table 1. Low frequency oscillation response of hydropower system to wind/solar/hydropower quota and line distance ratio.
Table 1. Low frequency oscillation response of hydropower system to wind/solar/hydropower quota and line distance ratio.
Line Length
Ratio (S:W)
W:S:HEigenvalueFrequency/HZDamping Ratio/%Remark
1:120:1:150−1.14 × 10−1 + j00.0182100Three-Phase Parallel RLC Load
−1.26 × 101 + j02.01100Grounding Transformer
−9.18 × 102 − j4.47 × 1037.26 × 10220.1Three-Phase PI Section Line
−3.15 × 103 − j9.43 × 1031.58 × 10331.7Windformsystem/30 km line
30:1:150−1.13 × 10−1 + j00.018100Three-Phase Parallel RLC Load
−1.65 × 101 + j02.63100Grounding Transformer
−9.47 × 102 − j4.61× 1037.48 × 10220.1Three-Phase PI Section Line
−3.15× 103 − j9.43× 1031.58 × 10331.7Windformsystem/30 km line
40:1:150−1.13 × 10−1 + j00.018100Three-Phase Parallel RLC Load
−1.99 × 101 + j03.17100Grounding Transformer
−9.77 × 102 − j4.74 × 1037.71 × 10220.2Three-Phase PI Section Line
−3.15 × 102 − j9.43 × 1031.58 × 10331.7Windformsystem/30 km line
2:120:1:150−1.14 × 10−1 + j00.0182100Three-Phase Parallel RLC Load
−1.26 × 101 − j02.01100Grounding Transformer
−4.62 × 103 − j07.36 × 102100Three-Phase PI Section Line
−3.23 × 102 − j7.11 × 1031.13 × 1034.54Windformsystem/30 km line
30:1:150−1.13 × 10−1 + j00.018100Three-Phase Parallel RLC Load
−1.65 × 101 + j2.37 × 10−102.63100Grounding Transformer
−9.67 × 102 − j4.58 × 1037.45 × 10220.7Three-Phase PI Section Line
−3.22 × 102 − j7.11 × 1031.13 × 1034.53Windformsystem/30 km line
40:1:150−1.13 × 10−1 + j00.018100Three-Phase Parallel RLC Load
−1.99 × 101 + j03.17100Grounding Transformer
−9.98 × 102 − j4.71 × 1037.66 × 10220.7Three-Phase PI Section Line
−3.22 × 102 − j7.11 × 1031.13 × 1034.52Windformsystem/30 km line
3:120:1:150−1.14 × 10−1 + j00.0182100Three-Phase Parallel RLC Load
−1.26 × 101 + j02.01100
−8.82 × 102 − j4.39 × 1037.13 × 10219.7Three-Phase PI Section Line
−5.85 × 102 + j627 × 1031.00 × 1039.28Windformsystem/30 km line
30:1:150−1.13 × 10−1 + j00.018100Three-Phase Parallel RLC Load
−1.65 × 101 + j02.63100Grounding Transformer
−1.72 × 102 − j4.48 × 1037.13 × 1023.85Three-Phase PI Section Line
−5.85 × 102 + j6.27 × 1031.00 × 1039.28Windformsystem/30 km line
40:1:150−1.13 × 10−1 + j00.018100Three-Phase Parallel RLC Load
−1.99 × 101 − j03.17100Grounding Transformer
−1.76 × 102 − j4.47 × 1037.12 × 1023.95Three-Phase PI Section Line
−5.85 × 102 + j6.27 × 1031.00 × 1039.28Windformsystem/30 km line
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, B.; Lei, L.; Zhao, Z.; Jiang, W.; Xiao, S.; Li, H.; Zhang, J.; Chen, D. Low Frequency Oscillations in a Hydroelectric Generating System to the Variability of Wind and Solar Power. Water 2021, 13, 1978. https://doi.org/10.3390/w13141978

AMA Style

Xu B, Lei L, Zhao Z, Jiang W, Xiao S, Li H, Zhang J, Chen D. Low Frequency Oscillations in a Hydroelectric Generating System to the Variability of Wind and Solar Power. Water. 2021; 13(14):1978. https://doi.org/10.3390/w13141978

Chicago/Turabian Style

Xu, Beibei, Liuwei Lei, Ziwen Zhao, Wei Jiang, Shu Xiao, Huanhuan Li, Junzhi Zhang, and Diyi Chen. 2021. "Low Frequency Oscillations in a Hydroelectric Generating System to the Variability of Wind and Solar Power" Water 13, no. 14: 1978. https://doi.org/10.3390/w13141978

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