Next Article in Journal
Reducing Energy Consumption and Increasing the Performances of AC Motor Drives Using Fuzzy PI Speed Controllers
Next Article in Special Issue
Influence of Flow Divider on Overall Efficiency of a Hydrostatic Drivetrain of a Skid-Steer All-Wheel Drive Multiple-Axle Vehicle
Previous Article in Journal
Methodology for Estimating the Spatial and Temporal Power Demand of Private Electric Vehicles for an Entire Urban Region Using Open Data
Previous Article in Special Issue
Influence of the Reshaped Elbow on the Unsteady Pressure Field in a Simplified Geometry of the Draft Tube
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Frequency Containment Control of Hydropower Plants Using Different Adaptive Methods

by
Doğan Gezer
1,*,
Yiğit Taşcıoğlu
2 and
Kutay Çelebioğlu
3
1
TUBITAK MRC Energy Institute, METU Campus, Ankara 06531, Turkey
2
Department of Mechanical Engineering, TED University, Ankara 06420, Turkey
3
TOBB University of Economics and Technology, ETU Hydro, Ankara 06560, Turkey
*
Author to whom correspondence should be addressed.
Energies 2021, 14(8), 2082; https://doi.org/10.3390/en14082082
Submission received: 13 February 2021 / Revised: 5 April 2021 / Accepted: 7 April 2021 / Published: 8 April 2021

Abstract

:
With the growth in the share of variable renewable energy sources, fluctuations in the power generation caused by these types of power plants can diminish the stability and flexibility of the grid. These two can be enhanced by applying frequency containment using hydropower plants as an operational reserve. The frequency containment in hydropower plants is automatically controlled by speed governors within seconds. Disturbances such as fluctuations in the net head and aging may diminish the performance of the controllers of the speed governors. In this study, model reference adaptive control approaches based on the Massachusetts Institute of Technology (MIT) rule and Lyapunov method were exploited in order to improve the performance of the speed governor for frequency containment control. The active power control with frequency control was enhanced by the aforementioned adaptive control methods. A mathematical model of a hydropower plant with a surge tank and medium penstock was constructed and validated through site measurements of a plant. It was shown that, as they are applicable in real life, both methods perform significantly better compared to conventional proportional-integrator control. Even in first five deviations, the performance of the conventional controller improved by 58.8% using the MIT rule and by 65.9% using the Lyapunov method. When the two adaptive control approaches were compared with each other, the MIT rule outputted better results than the Lyapunov method when the disturbance frequency was higher; however, the latter was more functional for rare disturbances.

1. Introduction

As the penetration of variable renewable energy sources (V-RES) such as wind and photovoltaics in the grid increases to higher proportions, a higher amount of operational reserve will be needed to provide grid stability and flexibility. Grid stability involves the availability of electricity at all times and a power quality covering the voltage and frequency stability. In a similar manner, grid flexibility means the capacity for load and generation fluctuations, especially resulting from the high share of V-RESs. The operation of hydropower plants (HPPs) with storage capacity (e.g., reservoir HPPs) is a key solution for maintaining the power quality and grid flexibility of an interconnected network, relying on their very short start-up and action times.
Considering the frequency stability, there are four steps of action that can be taken against the frequency deviation [1]. The most rapid one is frequency containment control (FCC), formerly known as primary frequency control, which is used for limiting the deviations in the grid frequency within a few seconds. These deviations result from the imbalance between demand (e.g., load fluctuations) and supply (e.g., intermittent generation of V-RESs. The units contributing to the frequency containment control automatically start this action. In HPPs, the speed governor of the unit is responsible for the FCC. Speed governors adjust the turbine speed for grid synchronization. Following the synchronization process, the speed governors regulate the active power output for network operation. As long as the FCC function is on, they are expected to increase/decrease the active power output of the unit depending on the sign of the frequency deviation, i.e., an increase in the frequency, meaning an excess of supply, will result in the automatic action of speed governors by decreasing the active power. As there exist other generating units contributing to the frequency response, each unit takes a portion of the required load according to its permanent speed droop value. An FCC action is traditionally implemented at the generation side in the existing control strategy of the speed governors. Recent studies on FCC have focused on the high penetration of V-RESs. Flywheel systems for inertia control were proposed in [2] for active power support in the case of frequency dips. A fuzzy-logic-based algorithm controls the bidirectional active power flow in a wind farm, according to the system frequency and the rate-of-change of system frequency. The possibility of FCC in smaller-alternative-current (AC) subnetworks that are connected to a high-voltage direct current (HVDC) network was investigated [3]. Different controller schemes were analyzed to benefit from the abundant inertial resources of rotors in wind power plants [4]. FCC studies for HPPs have concentrated on performance analyses of HPP models with/without surge tanks in [5,6] using frequency response, step response, and setpoint change tests.
The state of the art for hydraulic turbine control is digital governors, mostly with Proportional—Integral (PI) controllers and some with modern control techniques [7]. These controllers are tuned under ideal conditions with newly installed equipment. On the classical control side, PI control stability issues [8,9] and techniques for tuning PI parameters [10,11,12,13] have been studied by the majority. On the modern control side, various studies have covered nonlinear, optimal, and robust control using state variable representation. A self-tuning power system stabilizer, which uses the minimization of the quadratic performance index, was presented in [14]. A genetic algorithm was used in [15] to optimize the speed governor performance of an HPP. Speed governor designs with neural networks were described in [16] with an intelligent PI tuning method, in [17] for an isolated power system, and in [18] by the implementation of neuro-controllers. The robust control of hydraulic turbines through frequency-response was presented in [19]. In [20], a robust H∞ control was investigated for a very long headrace tunnel. Ref. [21] studied the multi-frequency dynamic performance of an HPP with a surge tank through the coupling effect among penstock, a water turbine, and the electrical grid. In [22], the chaotic motion of the speed governor was investigated while nonlinear turbine characteristics were present.
The performance of speed governors can be disrupted by the changes in the reservoir and tailwater levels, as well as aging effects, e.g., wear and tear in the flow-regulating mechanism. Model reference adaptive control (MRAC) can improve the performance of speed governors in the presence of these disturbances, as MRAC with the Massachusetts Institute of Technology (MIT) rule and Lyapunov method follows the desired behavior described by a reference model.
MRAC for FCC is only investigated in [23] for micro-grids. However, there is a gap in the literature for MRAC usage for FCC in interconnected networks.
In this study, the usage of MRAC for FCC in an interconnected grid was investigated. First, a dynamic HPP model and a reference model describing the ideal behavior for FCC were formed. Then, adaptive controllers, designed with the MIT rule and the Lyapunov method, were added into the main active power control loop of the speed governor to rapidly adapt the controller to the changing plant conditions. Finally, the simulation results of both controllers were compared.
The rest of the article is organized as follows: the HPP modeling is detailed in Section 2. The proposed controllers including the MIT rule and the Lyapunov method, along with the reference model, are presented in Section 3. The simulations and results are presented in Section 4, and Section 5 concludes the paper.

2. Hydropower Models

The control parameters of speed governors are traditionally determined by using linear models. This approach is convenient when the changes in the loading and speed of the unit are within 10% and 1%, respectively.
On the other hand, significant changes in the active power output and turbine speed of the unit occur during isolated (i.e., without being interconnected with other generating units) or islanded (i.e., interconnection with only a small number of units) operation, as well as during loading, unloading, load rejection, and turbine speed adjustment before network synchronization. It is necessary to simulate the speed governor controller and related systems with nonlinear models for changes larger than 25% in loading and 8% in turbine speed [24]. Nonlinear HPP models include the compressibility of water. An instance of this effect could be the propagating wave and hydraulic pressure oscillations in penstock during load rejections [25].
In this study, a nonlinear model was used to test all possible scenarios and a linear model was used to conduct a stability analysis of the controller. Figure 1 shows the general layout of an HPP including the feedbacks to the speed governor.

2.1. Nonlinear HPP Model

The modeling equations of the depicted components in the waterway, servomotor and guide vanes, turbine, and generator are given in the following subsections [24,25,26,27,28,29,30,31,32]. Many of the quantities below are given in the per-unit system, which is the expression of the quantities as a ratio of the actual value over the base value, and it is dimensionless. The base value is the rated value of the quantity during normal operation.

2.1.1. Waterway (Tunnel, Surge Tank, and Penstock)

The relationships between the water head in the turbine and the reservoir level, including the energy tunnel, surge tank, penstock dynamics, and friction loss, are given as follows:

Tunnel

The rate of change of water velocity through the tunnel is calculated as:
d U w t d t = H r H s H f t t w t
where Uwt is the water velocity through the tunnel in pu, Hr is the reservoir water level in pu, Hs is the water level in the surge tank in pu, Hft is the head loss due to friction in the tunnel in pu, and twt is the water starting time of the tunnel in s, which is calculated as:
t w t = Q 0 H 0 g l t A t
where Q0 is the nominal flow rate in m3/s, H0 is the nominal net head in m, g is the gravitational acceleration in m2/s, and lt and At are the length and cross-sectional area of tunnel portions in m and m2, respectively.

Surge Tank

The water level in the surge tank is found as:
H s = 1 C s · ( U w t   U w p ) · d t
where Cs is the storage capacity of the surge tank in s, and Uwt and Uwp are the water velocities, in pu, through the tunnel and penstock, respectively. The storage capacity is calculated as follows:
C s =   A s   H 0     Q 0
where As is the average cross-sectional area of the surge tank in m2.

Penstock

The rate of change of water velocity through the penstock is defined as:
d U w p d t = H s H t u r H f p   H w t w p
where Htur is the water head in the turbine in pu, Hfp is the head loss due to friction in the penstock in pu, Hw is the head change due to the propagating wave effect in pu, and twp is the water starting time of the penstock in s, which is calculated as follows:
t w p = Q 0 H 0 g l p A p
where Q0 is nominal flow rate in m3/s, H0 is the nominal net head in m, g is the gravitational acceleration in m2/s, and lp and Ap are the length and cross-sectional area of penstock portions in m and m2, respectively.
The water starting time or water time constant is the time needed to accelerate the water in the penstock from zero discharge to the rated discharge under the rated head. This parameter is used in both linear and nonlinear models. Its value depends on the length and cross-sectional area of the penstock, rated head, and discharge.

2.1.2. Turbine and Servo Cylinders

The relationships between the mechanical power output of the turbine and the guide vane opening, including the servo cylinder and the water head in front of the turbine, are given below:
The mechanical power output of the turbine is defined as follows:
P m = K t u r   ·   H t u r · ( U t u r U n l ) =   T m   ·   ω
where Pm is the mechanical power output of the turbine in pu, Ktur is the turbine gain, Htur is the water head in front of the turbine in pu, Utur and Unl are the water velocities, in pu, through the turbine at load and no-load, respectively, Tm is the mechanical torque in pu, and ω is the turbine speed in pu.
Utur depends on Htur and the guide vane opening G, and it is equal to Uwp:
U t u r = U w p = G H t u r
The turbine gain, Ktur, is defined as:
K t u r =   1   G f l   G n l
where Gfl and Gnl are, respectively, the full-load and no-load of the guide vane openings.
d G d t = 1 t s e r v o ( G s e t G ) e s t d e l _ s e r
where G is the guide vane opening in pu, tservo is the time constant of the servomotor in s, Gset is the guide vane opening setpoint, in pu, sent from the speed governor controller, and t d e l _ s e r is the time delay for servo opening in s.

2.1.3. Generator and Excitation System

The equations relating the turbine speed with the electrical torque, including field voltage, are given as follows:
2 H   d ω d t =   T m   T e
where Te is the electromagnetic torque in pu, and H is the inertia constant, in s, which depends on the flywheel effect [29] as:
H = 1 2 ( π 60 ) 2 G D 2   n 2 S
where GD2 is the flywheel effect of the generator in tone.m2, n is the turbine speed in rpm, and S is the nominal apparent power of the generator in kVA.
The generator in the steady-state is electrically modeled by considering the simplifications explained in [30]:
P e =   | E q |   | V s | X d +   X e sin δ
where Eq is the generator internal voltage in pu, Vs is the voltage of the electrical network bus in pu, δ is the angle between E q δ ° and V s 0 ° , Xd is the direct axis reactance in pu, Xe is the bus impedance in pu, and:
E q =   n   M f   V f r f
where Vf is the main field voltage in pu, Mf is the mutual inductance between the field and any armature phase in pu, and rf is the resistance of the field circuit in pu.
The interactions between the subsystems of the generic model are depicted in Figure 2. The implementation of the generic nonlinear model is given in Figure 3.

2.2. Linear HPP Model

Linear models are used for modeling small changes in the turbine speed and the active power of the HPP. Stability issues for the controller are investigated through the linear model. Linear model equations are obtained by the linearization of the aforementioned nonlinear equations around an operating point. The linearization of relationships between the mechanical power output of the turbine, the head in front of the turbine, the velocity of the water entering the turbine, and the guide vane opening rooting from Equations (7) and (8) result in a linear equation linking the mechanical power output of the turbine and the guide vane opening in Equation (15).
Δ P m   Δ G = 1 t w p s 1 + 1 2 t w p s
where t w p is the water starting time of the penstock in s.
When the natural frequency effect is taken into account, the dynamic behavior due to the presence of penstock and the surge tank is obtained as in Equation (16) [33].
Δ P m   Δ G = s 2 ω s 2 t w p s + ω s 2 s 2 + 1 2 ω s 2 t w p s + ω s 2 ω s = g A p l p A s }
where t w p is the water starting time of the penstock in s and ωs is the natural frequency of the mass oscillation due to the surge tank in rad/s. ωs could be calculated by using lp and Ap, which are the length and cross-sectional area of the penstock in m and m2, respectively, and As, which is the cross-sectional area of the surge tank in m2.

2.3. System Description

In this study, Seyhan I HPP was modeled with the characteristic parameters listed in Table 1. Seyhan I HPP is an aged hydropower plant that started operation in 1957. It has three vertical Francis turbines with 22.5 MVA of installed capacity each.
A mathematical model of Seyhan I HPP is based on the following assumptions and modeling principles:
  • The modeled HPP has a short-medium penstock and a surge tank, but it does not have an energy tunnel. Although there are no common or standardized lengths to classify the penstocks, if any fluctuations in the active power output are observed due to wave effects within the penstock, the penstock can be considered as medium-long and the wave effect should be taken into account while modeling the waterway. Hence, wave propagation is not considered;
  • The electrical dynamics of the generator has a very short time constant compared to that of hydrodynamics [34]. Electromagnetic interactions in the generator occur much faster than the speed governor control actions do; hence, electromagnetic generator dynamics are not considered. On the other hand, during the adjustment of the turbine speed before network synchronization, it is critical to pay regard to the mechanical starting time (also known as mechanical inertia time) for tuning speed governor parameters. Furthermore, the speed rise following a load rejection is limited by the rotational inertia of the unit. The turbine rotational inertia is approximately 5% of the generator rotational inertia [27], meaning that the rotational inertia of the unit dominantly depends on the generator characteristic rather than the hydraulic turbine. Therefore, the rotational inertia of the generator is taken into account;
  • The modelling of the electrical grid and load is necessary for simulations of the speed governor in isolated and islanded modes of operation, where small load changes result in significant frequency deviations. An interconnected mode of operation, as preferred in this study, requires the network frequency to be kept constant [35].

2.4. Model Validation

As the validation methodology, the mathematical model was constructed considering the aforementioned assumptions and modeling principles. Then, site tests were conducted for unit start-up, shut-down, and loading. The guide vane opening, the turbine speed, and the active power measurements were recorded during the site tests. Following the site tests, the measured guide vane opening values were given as workspace inputs to the mathematical model in MATLAB/Simulink. The simulations conducted led to the simulated turbine speed and the simulated active power. These simulation results were compared with the actual site measurements in Figure 4a–c during start-up and shut-down, and the unit was disconnected from the grid. The simulated turbine speed profiles were in close agreement with the site measurements. In Figure 4b, the unit was already synchronized with the network; hence, the guide vane opening affected the active power output. Despite the acceptable lead in the simulated response, loading profiles were in correspondence, and they attained the same active power output.
In Figure 5, different models are compared with the measured active power. Each model was given the measured guide vane opening values as input. The first linear model was constructed by using the relationship in Equation (15). The second was established using Equation (16). The third was the nonlinear model generated by Equations (1)–(14). Although the nonlinear model had the most correspondence, in order to use it in stability issues, the model by Equation (16) was the most convenient one.

3. Proposed Adaptive Methods

Speed governors take advantage of active power control in order to follow the active power setpoint given by the operator. Figure 6 shows that the input of the main PI controller (with proportional gain, Kp, and the integral time, Ti) has two branches with different setpoints and feedbacks. One is the active power control branch, where the error is the difference between the given active power setpoint, Pset, and the measured active power, P. The other is the FCC branch, through which the frequency deviation (fset-f) exceeding the speed dead-band is added to the active power control branch after being divided by the permanent speed droop, R. Figure 6 shows the inside of the speed governor in Figure 1.
Adaptive control is used for the automatic adjustment of controllers in real-time to maintain the desired performance. It is an efficient technique for dealing with deviations and uncertainties in slow-varying plant parameters [36,37].
MRAC is a direct adaptive control approach, in which the desired performance of the plant controller is obtained through a dynamic model called the reference model. This controller has an ordinary feedback loop taken from the process and another loop for adjusting the controller parameters. As displayed in Figure 7, the controller—which outputs a processed setpoint through comparing the given setpoint and the feedback signal, the plant in which the process occurs, the reference model by which the desired performance of the plant is modeled, and the adaptation mechanism that compares the outputs of the reference model and the plant—and MRAC are the main parts in the adaptive structure. The mechanism for adjusting the parameters can be developed by the gradient method (e.g., MIT rule) or stability theory (e.g., Lyapunov method).

3.1. Model Reference Adaptive Controller Using MIT Rule

The gradient design method is based on the MIT rule, which is used for the computation of approximate sensitivity functions. The adaptive law is designed such that the controller’s parameters are tuned in a direction, which decreases the quadratic performance function. For this purpose, the partial derivative of the cost function is used. The output of the adaptation mechanism, θ , and its relationship with the learning rate, Γ , the output of the reference model, Y m , and the error, e , are written in Equation (17) [38,39].
P l a n t   T F :   G p = y p u   M o d e l   T F : G m = y m u   C o n t r o l   L a w : u = θ u E r r o r : e =   y p   y m C o s t   F u n c t i o n :   J ( t , θ ) = | e ( t , θ ) | 2 M I T   R u l e :   d d t θ = γ Θ J = γ [ 2 e Θ e ] θ = Γ [ Y m ( t ) e ( t , θ ) d t ] }
Figure 8 shows the block diagram of the proposed controller including MRAC with the MIT rule. In addition to the controller scheme shown in Figure 6, the reference model and adaptation mechanism are embedded. The frequency containment error is also given as an input to the reference model. A transport delay is added to represent the latency in the first movement of the guide vanes. The limiter in the reference model confines the reference model output in the predetermined active power reserve of the power plant. The difference between the ideal behavior through the reference model and the actual primary control action is given as an input to the adaptation mechanism. The output of the PI controller is corrected by the output θ of the adaptation mechanism. In most speed governor systems, the output of the controller is limited by the load limiter. The controller output is sent to the guide vane mechanism actuator.

3.2. Model Reference Adaptive Controller Using Lyapunov Method

The stability method relies on Lyapunov stability theory. This theory assures the stability and convergence of the error between the plant output and the model output. The procedure for developing the MRAC with the Lyapunov method starts with the derivation of the differential equation for the error. Then, a Lyapunov function is found that makes the error go to zero. In order to drive the error to zero, a Lyapunov function is used with a symmetric positive definite matrix P. The derivative of the Lyapunov function is made negative by solving the Lyapunov function with a positive definite matrix Q. The output of the adaptation mechanism, θ , and its relationship with the learning rate, Γ , the output of the reference model, Y m , and the error, e , are given in Equation (18) [40,41,42].
P l a n t   T F :   G p = y p u     M o d e l   T F : G m = y m u   C o n t r o l   L a w : u = θ u E r r o r : e =   y p   y m L y a p u n o v   F u n c t i o n :   V ( x ) = x T P x   D e r i v a t i v e :   V ˙ = 1 2 e T ( t ) Q e ( t ) θ = Γ [ u ( t ) e ( t , θ ) d t ] }
Figure 9 represents the block diagram for the model reference adaptive controller based on Lyapunov stability. The practical difference of this method compared to the MIT rule is using u ( t ) instead of Y m ( t ) as an input to the multiplication block within the adaptation mechanism.

3.3. Reference Model for Frequency Containment Control

Frequency containment control is the fastest control action taken against frequency deviations in the grid. All units, contributing to FCC, automatically give active power support by increasing/decreasing their active power output depending on the sign of the frequency deviation. An FCC action is expected to start immediately following the frequency deviation. All FCC reserves are to be given in 30 s and last for 15 min. In particular, hydropower plants are permitted to start this action with a 4 s delay at maximum [43,44].
As shown in Figure 10, the amount of active power contribution of each unit is determined according to the predefined permanent speed droop value. The responses of all interconnected units contributing to the FCC are coordinated with this value. The amount of active power support, ∆P, of the unit depends on the amount of frequency deviation, ∆f, and the permanent speed droop, R, i.e., a 200 mHz decrease in frequency of a 50 Hz network requires a 0.1 pu increase in active power output with a permanent speed droop of 0.04. The relationship between the frequency deviation, permanent droop, and change in active power output is given in Equation (19).
Δ f = R   Δ P
For the adaptive controller, a reference model is needed to define the desired behavior of the plant. A first-order transfer function, as given in Equation (20), is an appropriate reference model for the frequency containment action. A time constant, t P F C , of 6 s is used in the simulations.
P m o d e l f d e v i a t i o n   =   1 t P F C s + 1   1 R
The defined reference model with its permitted limits is depicted in Figure 11. Within 90 s after a negative frequency deviation, (e.g., f = 49.8 Hz), the actual control action should not exceed the ideal action by 2% of the nominal generator power output [45].

4. Results

4.1. Stability of the Power Control Loop

The parameters of the main PI controller shown in Figure 6 are Kp = 0.25 and Ti = 11 s. Figure 12 shows the root locus of the closed-loop system, with no active MRAC, for the different values of the proportional gain. The low-frequency pair of poles crosses to the right-hand side of the complex plane at Kp = 0.388. As the output of the MRAC effectively amplifies Kp, the closed-loop system is stable up to θ = 1.55.

4.2. Simulation Scenarios

4.2.1. Change in Net Head

The reservoir level of an HPP could change from 0.9 pu, during an arid season, to 1.1 pu, following a rainy or snowy season. The presence of other units in the power plant, which are in operation, also raises the tailwater level or decreases the inlet pressure for a common penstock. The contribution of the MRAC for FCC is more perceivable when the net head drops to 0.9 pu. The speed governor opens the guide vanes more than that in the nominal head to reach the same active power output. This introduces latency in the FCC behavior. The Hr parameter is adjusted to reflect this scenario on the mathematical model.

4.2.2. Degradation in Turbine Efficiency

The turbine efficiency could dramatically drop over the years due to aging. It is necessary to increase the guide vane opening to reach the same value of active power when the efficiency drops. This scenario could be simulated by increasing Gfl, which represents the guide vane opening at the nominal active power output. The increase in this parameter affects the turbine gain, Ktur.

4.2.3. Deceleration in Guide Vane Driving

The elements in the guide vane driving system are the servo cylinder, hydraulic power unit, regulation ring, and the guide vanes. Any aging, clogging, and wear in this system could slow down the movement of the guide vanes. By changing the values of ts and t d e l _ s e r parameters, this scenario will be reflected in the simulation.

4.3. Simulation Results

Figure 13 shows the simulation results for the net head of 0.9 pu for the MRAC. The simulation conditions are as follows: The reference model matches the actual FCC action for the nominal head (i.e., net head of 1.0 pu). The learning rate, Γ, is 1 as long as the reference model corresponds to the actual FCC action. The permanent speed droop value, R, is 0.04. The FCC setpoint value changes between 0 and 0.1 as a result of the 0.004 pu frequency deviation, meaning a 49.8 Hz network frequency. Figure 13a displays the results for MRAC with the MIT rule; similarly, Figure 13b depicts the simulations for MRAC with the Lyapunov method. A pseudo-frequency deviation signal is given to the validated mathematical model. It is seen that when MRAC is not active, the action for FCC does not meet the requirements. It takes nearly 60 s to provide the necessary FCC reserve; however, it is expected to be given in 30 s. However, the FCC action satisfies the requirements as soon as MRAC is activated even when the Γ value is low with a very low number of occurrences (e.g., 0.8 after five cycles). As the number of the occurrences is increased to 50 and the Γ value is raised to 5, the FCC with the MIT rule action is nearly the same as the reference model. When focused on FCC with the Lyapunov method, it is seen that FCC with the Lyapunov method also satisfies the FCC requirements.
The change in θ for both the MIT rule and Lyapunov method is depicted in Figure 14. One should note that θ converged to the same value and the convergence was faster with greater Γ values until a critical value. During the simulations, it was seen that when a critical Γ value was exceeded, θ displayed unstable behavior. The critical Γ value was found to be 5.6 for MRAC with the MIT rule; similarly, the critical Γ value was found to be 7.1 for MRAC with the Lyapunov method. For the MRAC with the MIT rule, the θ value converged to 1.15 independent of the value of the learning rate, Γ. For the MRAC with the Lyapunov method, the θ value converged to 1.45. Considering these saturation values, both MRAC controllers were within the stability margin, as explained in the previous section. However, a greater θ value means a greater change in guide vane opening. As the θ value was greater in the Lyapunov method, it is critical to limit the guide vane opening for secure operation.
The insufficient response, i.e., when the MRAC was off, was set as the baseline performance for both the MIT rule and the Lyapunov method. The results of the simulations with different Γ values and the number of cycles for both methods (i.e., recursive occurrences of frequency deviation) are displayed on Figure 15 and listed in Table 2. The sum of the squared errors (e2) to the response of the reference model was normalized to 1 pu. After 5 cycles with Γ = 0.8, the improvement was 58.8% for the MIT rule and was 65.9% for the Lyapunov method. This shows that for the same number of occurrences and the learning rate, the Lyapunov method had a better improvement in performance. Similarly, after 50 cycles with Γ = 5, the improvement was 97.5% for the MIT rule and was 96.5% for the Lyapunov method. Considering the top performance, the MIT rule performed better than the Lyapunov method.

5. Conclusions

Hydropower generation plays a crucial role in modern electric power systems to retard fluctuations caused by V-RESs. The FCC contribution of HPPs helps in decreasing the imbalances between supply and demand for the whole grid. In this study, model reference adaptive controllers with the MIT rule and Lyapunov method were designed for the speed governor controller of an HPP to improve the FCC performance of the unit.
The very first conclusion is that the performance of the conventional PI controller can be diminished by external disturbances, as shown in the simulations. Therefore, it is useful to add precautionary loops into the controller in order to provide the tuned performance. Secondly, adaptive controllers are very helpful in real life because it is not necessary to adjust adaptive parameters each time the disturbance occurs; however, the very first tuning of the adaptive controllers’ parameters requires considerable experience or a high level of modeling. As the third conclusion, the MRAC with both the MIT rule and Lyapunov method in the active power control loop makes considerable improvements in FCC against slow-varying disturbances such as variations in the net head. Both adaptive methods are also convenient to implement on the off-the-shelf controllers such as programmable logical controllers (PLCs).
The final conclusion is on the comparison of the adaptive approaches. It was seen that the Lyapunov method has a better performance for a small number of occurrences (1–5 times); however, disregarding the time scales, the top performance of the MIT rule was greater than the Lyapunov method. For the same number of occurrences, the θ value of the Lyapunov method reaches higher than the MIT rule, meaning that the guide vanes will be opened more in the Lyapunov method. This can bring a practical precaution of adding a limiter for this potential risk.

Author Contributions

Conceptualization, D.G., Y.T. and K.Ç.; methodology, D.G., Y.T. and K.Ç.; software, D.G.; validation, D.G.; formal analysis, D.G. and Y.T.; investigation, D.G.; resources, D.G.; data curation, D.G. and Y.T.; writing—original draft preparation, D.G., Y.T. and K.Ç.; writing—review and editing, D.G., Y.T. and K.Ç.; visualization, D.G. and Y.T.; supervision, Y.T. and K.Ç.; project administration, D.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

The authors gratefully acknowledge the contributions of TUBITAK MRC through supporting this research study by the project “Seyhan I HPP Rehabilitation Project” numbered 5132804 and conducted at TUBITAK, Scientific and Technological Research Council of Turkey, Marmara Research Center (MRC) Energy Institute.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Network of Transmission System Operators for Electricity, ENTSOE. Available online: https://www.entsoe.eu/fileadmin/user_upload/_library/news/ENTSO-_Stakeholder_Workshop_LFC___RNC_on_12_July/120606_Operational_Reserves_Report_public.pdf (accessed on 8 November 2018).
  2. Jauch, C.; Hippel, S. Hydraulic–pneumatic flywheel system in a wind turbine rotor for inertia control. IET Renew. Power Gener. 2016, 10, 33–41. [Google Scholar] [CrossRef]
  3. Bianchi, F.D.; Domnguez-Garca, J.L. Coordinated frequency control using MT-HVDC grids with wind power plants. IEEE Trans. Sustain. Energy 2016, 7, 213–220. [Google Scholar] [CrossRef]
  4. Wang-Hansen, M.; Josefsson, R.; Mehmedovic, H. Frequency Controlling Wind Power Modeling of Control Strategies. IEEE Trans. Sustain. Energy 2013, 4, 954–959. [Google Scholar] [CrossRef]
  5. Yang, W.; Yang, J.; Guo, W.; Norrlund, P. Response Time for Primary Frequency Control of Hydroelectric Generating Unit. Int. J. Electr. Power Energy Syst. 2016, 74, 16–24. [Google Scholar] [CrossRef] [Green Version]
  6. Saarinen, L.; Norrlund, P.; Lundin, U. Field Measurements and System Identification of Three Frequency Controlling Hydropower Plants. IEEE Trans. Energy Convers. 2015, 30, 1061–1068. [Google Scholar] [CrossRef]
  7. Fasol, K.H. A short history of hydropower control. IEEE Control Syst. 2002, 22, 68–76. [Google Scholar] [CrossRef]
  8. Hagihara, S.; Yokota, H.; Goda, K.; Isobe, K. Stability of a Hydraulic Turbine Generating Unit Controlled by P.I.D. Governor. IEEE Trans. Power Appar. Syst. 1979, PAS-98, 2294–2298. [Google Scholar] [CrossRef]
  9. Phi, D.T.; Bourque, E.J.; Thorne, D.H.; Hill, E.F. Analysis and Application of the Stability Limits of a Hydro-Generating Unit. IEEE Trans. Power Appar. Syst. 1981, PAS-100, 3203–3212. [Google Scholar] [CrossRef]
  10. Wozniak, L. A graphical approach to hydrogenerator governor tuning. IEEE Trans. Energy Convers. 1990, 5, 417–421. [Google Scholar] [CrossRef]
  11. Kamwa, I.; Lefebvre, D.; Loud, L. Small signal analysis of hydro-turbine governors in large interconnected power plants. In Proceedings of the 2002 IEEE Power Engineering Society Winter Meeting. Conference Proceedings (Cat. No.02CH37309), New York, NY, USA, 27–31 January 2002; Volume 2, pp. 1178–1183. [Google Scholar]
  12. Eilts, L.E.; Schleif, F.R. Governing features and performance of the first 600-MW hydrogenerating unit at grand coulee. IEEE Trans. Power Appar. Syst. 1977, 96, 457–466. [Google Scholar] [CrossRef]
  13. Filbert, T.L.; Wozniak, L. Speed loop cancellation governor for hydrogenerators. II. Application. IEEE Trans. Energy Convers. 1988, 3, 91–94. [Google Scholar] [CrossRef]
  14. Lim, C.M. A self-tuning stabiliser for excitation or governor control of power systems. IEEE Trans. Energy Convers. 1989, 4, 152–159. [Google Scholar] [CrossRef]
  15. Lansberry, J.E.; Wozniak, L.; Goldberg, D.E. Optimal hydrogenerator governor tuning with a genetic algorithm. IEEE Trans. Energy Convers. 1992, 7, 623–630. [Google Scholar] [CrossRef]
  16. Yamamoto, T.; Kaneda, M.; Oki, T.; Watanabe, E.; Tanaka, K. Intelligent tuning PID controllers. In Proceedings of the 1995 IEEE International Conference on Systems, Man and Cybernetics. Intelligent Systems for the 21st Century, Vancouver, BC, Canada, 22–25 October 1995; Volume 3, pp. 2610–2615. [Google Scholar]
  17. Padrón, S.; Hernández, M.; Falcón, A. Reducing Under-Frequency Load Shedding in Isolated Power Systems Using Neural Networks. Gran Canaria: A Case Study. IEEE Trans. Power Syst. 2016, 31, 63–71. [Google Scholar] [CrossRef]
  18. Venayagamoorthy, G.K.; Harley, R.G.; Wunsch, D.C. Implementation of adaptive critic-based neurocontrollers for turbogenerators in a multimachine power system. IEEE Trans. Neural Netw. 2003, 14, 1047–1064. [Google Scholar] [CrossRef] [PubMed]
  19. Natarajan, K. Robust PID Controller Design for Hydroturbines. IEEE Trans. Energy Convers. 2005, 20, 661–667. [Google Scholar] [CrossRef]
  20. Fangle, Q.; Guo, W. Robust H∞ control for hydro-turbine governing system of hydropower plant with super long headrace tunnel. Electr. Power Energy Syst. 2021, 124, 106336. [Google Scholar]
  21. Liu, Y.; Guo, W. Multi-frequency dynamic performance of hydropower plant under coupling effect of power grid and turbine regulating system with surge tank. Renew. Energy 2021, 171, 557–581. [Google Scholar] [CrossRef]
  22. Xu, X.; Guo, W. Chaotic behavior of turbine regulating system for hydropower station under effect of nonlinear turbine characteristics. Sustain. Energy Technol. Assess. 2021, 44, 101088. [Google Scholar]
  23. Mahdi, M.M.; Mhawi Thajeel, E.; Ahmad, A.Z. Load Frequency Control for Hybrid Micro-grid Using MRAC with ANN Under-sudden Load Changes. In Proceedings of the 2018 Third Scientific Conference of Electrical Engineering (SCEE), Baghdad, Iraq, 19–20 December 2018; pp. 220–225. [Google Scholar]
  24. Kishor, N.; Saini, R.P.; Singh, S.P. A review on Hydropower Plant Models and Control. Renew. Sustain. Energy Rev. 2007, 11, 776–796. [Google Scholar] [CrossRef]
  25. Ramey, D.G.; Skooglund, J.W. Detailed Hydrogovernor Representation for System Stability Studies. IEEE Trans. Power Appar. Syst. 1970, PAS-89, 106–112. [Google Scholar] [CrossRef]
  26. Mover, Working Group Prime, and Energy Supply. Hydraulic turbine and turbine control models for system dynamic studies. IEEE Trans. Power Syst. 1992, 7, 167–179. [Google Scholar] [CrossRef]
  27. IEEE. Guide for the Application of Turbine Governing Systems for Hydroelectric Generating Units. In IEEE Std 1207-2004; IEEE: New York, NY, USA, 2004; pp. 1–140. [Google Scholar] [CrossRef]
  28. Murat, D.; Kosalay, I.; Gezer, D.; Sahin, C. Validation of hydroelectric power plant model for speed governor development studies. In Proceedings of the 2015 International Conference on Renewable Energy Research and Applications (ICRERA), Palermo, Italy, 22–25 November 2015; pp. 278–282. [Google Scholar]
  29. Piriz, H.D.; Cannatella, A.R.; Guerra, E.; Porcari, D.A. Inertia of Hydrogenerators, Influence on the Dimensioning, Cost, Efficiency and Performance of the Units. In Proceedings of the CIGRE Paris Session, Paris, France, 27–31 August 2012; pp. A1–A102. [Google Scholar]
  30. Fischer, N.; Benmouyal, G.; Samineni, S. Tutorial on the Impact of the Synchronous Generator Model on Protection Studies. In Proceedings of the 35th Annual Western Protective Relay Conference, Spokane, WA, USA, 20–22 October 2008. [Google Scholar]
  31. Kundur, P. Power System Stability and Control; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]
  32. Report, IEEE Committee. Computer representation of excitation systems. IEEE Trans. Power Appar. Syst. 1968, PAS-87, 1460–1464. [Google Scholar] [CrossRef]
  33. Vournas, C.D.; Papaioannou, G. Modelling and stability of a hydro plant with two surge tanks. IEEE Trans. Energy Convers. 1995, 10, 368–375. [Google Scholar] [CrossRef]
  34. Munoz-Hernandez, G.A. Modelling and Controlling Hydropower Plants; Springer: New York, NY, USA, 2013. [Google Scholar]
  35. Cebeci, M.E. The Effects of Hydro Power Plants’ Governor Settings on the Stability of Turkish Power System Frequency. Master’s Thesis, Middle East Technical University, Ankara, Turkey, 2009. [Google Scholar]
  36. Landau, I.D.; Lozano, R.; M’Saad, M.; Karimi, A. Adaptive Control: Algorithms, Analysis and Applications; Springer: New York, NY, USA, 2011. [Google Scholar]
  37. Lavretsky, E. Adaptive control: Introduction, overview, and applications. In Proceedings of the NASA Adaptive Control Workshop, NASA Marshall Space Center, Huntsville, AL, USA, 17–20 March 2009. [Google Scholar]
  38. Ioannou, P.A. Robust Adaptive Control; PTR Prentice-Hall: New York, NY, USA, 1996. [Google Scholar]
  39. Sevcik, K. Model Reference Adaptive Control MRAC Tutorial. Available online: http://www.pages.drexel.edu/~kws23/tutorials/MRAC/MRAC.html (accessed on 15 November 2016).
  40. Astrom, K.J.; Wittenmark, B. Adaptive Control, 2nd ed.; Addison-Wesley Longman Publishing Co., Inc.: Boston, MA, USA, 1994. [Google Scholar]
  41. Perdukova, D.; Fedor, P.; Fedak, V.; Padmanaban, S. Lyapunov Based Reference Model of Tension Control in a Continuous Strip Processing Line with Multi-Motor Drive. Electronics 2019, 8, 60. [Google Scholar] [CrossRef] [Green Version]
  42. Oltean, S.E.; Dulau, M.; Duka, A.V. Model Reference Adaptive Control Design for Slow Processes. A Case Study on Level Process Control. Procedia Technol. 2016, 22, 629–636. [Google Scholar] [CrossRef] [Green Version]
  43. NERC. North American Electric Reliability Corporation, Balancing and Frequency Control; NERC: Princeton, NJ, USA, 2011. [Google Scholar]
  44. UCTE. Union for the Coordination of the Transmission of Electricity, UCTE Operation Handbook; UCTE: Brussels, Belgium, 2004. [Google Scholar]
  45. Energy Market Regulatory Authority, Regulation Concerning Electricity Network. Available online: https://www.epdk.gov.tr/Detay/Icerik/3-6730/elektrik--sebeke-.html (accessed on 15 February 2021).
Figure 1. General layout of a hydropower plant (HPP).
Figure 1. General layout of a hydropower plant (HPP).
Energies 14 02082 g001
Figure 2. The interactions between subsystems in an HPP.
Figure 2. The interactions between subsystems in an HPP.
Energies 14 02082 g002
Figure 3. The schematic of the mathematical model implementation.
Figure 3. The schematic of the mathematical model implementation.
Energies 14 02082 g003
Figure 4. Model validation of Seyhan I HPP with site measurements during (a) start-up, (b) loading, and (c) shutdown.
Figure 4. Model validation of Seyhan I HPP with site measurements during (a) start-up, (b) loading, and (c) shutdown.
Energies 14 02082 g004
Figure 5. Comparison of various models with site measurements.
Figure 5. Comparison of various models with site measurements.
Energies 14 02082 g005
Figure 6. Active power control with frequency containment.
Figure 6. Active power control with frequency containment.
Energies 14 02082 g006
Figure 7. General schematic for model reference adaptive control (MRAC).
Figure 7. General schematic for model reference adaptive control (MRAC).
Energies 14 02082 g007
Figure 8. Block diagram of the controller with MIT rule.
Figure 8. Block diagram of the controller with MIT rule.
Energies 14 02082 g008
Figure 9. Block diagram of the controller with Lyapunov function.
Figure 9. Block diagram of the controller with Lyapunov function.
Energies 14 02082 g009
Figure 10. Change in active power output for different permanent speed droop values.
Figure 10. Change in active power output for different permanent speed droop values.
Energies 14 02082 g010
Figure 11. The step response of the reference model.
Figure 11. The step response of the reference model.
Energies 14 02082 g011
Figure 12. Root locus of the power control loop.
Figure 12. Root locus of the power control loop.
Energies 14 02082 g012
Figure 13. The FCC simulations for 0.9 pu net head: (a) MRAC with MIT rule; (b) MRAC with Lyapunov method.
Figure 13. The FCC simulations for 0.9 pu net head: (a) MRAC with MIT rule; (b) MRAC with Lyapunov method.
Energies 14 02082 g013
Figure 14. Change in θ for (a) MRAC with the MIT rule and (b) MRAC with the Lyapunov method.
Figure 14. Change in θ for (a) MRAC with the MIT rule and (b) MRAC with the Lyapunov method.
Energies 14 02082 g014
Figure 15. Performance index as a function of Γ and # of cycles: (a) MRAC with MIT rule; (b) MRAC with Lyapunov method.
Figure 15. Performance index as a function of Γ and # of cycles: (a) MRAC with MIT rule; (b) MRAC with Lyapunov method.
Energies 14 02082 g015
Table 1. Characteristic parameters of Seyhan I HPP.
Table 1. Characteristic parameters of Seyhan I HPP.
ParameterValueUnit
Nominal net head32m
Nominal flow rate77m3/s
Surge tank storage capacity124.8s
Surge tank cross-sectional area380m2
Penstock length81.7m
Penstock cross-sectional area21.2m2
Penstock friction factor0.05-
Guide vane opening at no-load11%
Guide vane opening at full-load85%
Servomotor time constant4s
Guide vane delay time2s
Flywheel effect of the generator3320tone.m2
Turbine nominal speed125rpm
Apparent power of the generator22.5MVA
Inertia constant of the generator3.14s
Table 2. Results for the contribution of MRAC controllers.
Table 2. Results for the contribution of MRAC controllers.
CaseNumber of Cycles e 2 Improvement wrt Baseline
TypeΓ
MIT rule--1 pubaseline
MIT rule0.850.412 pu58.8%
MIT rule0.8100.272 pu72.8%
MIT rule0.8500.047 pu95.3%
MIT rule250.237 pu76.3%
MIT rule2100.118 pu88.2%
MIT rule2500.026 pu97.4%
MIT rule550.112 pu88.8%
MIT rule5100.041 pu95.9%
MIT rule5500.025 pu97.5%
Lyapunov0.850.341 pu65.9%
Lyapunov0.8100.237 pu76.3%
Lyapunov0.8500.039 pu96.1%
Lyapunov250.189 pu81.1%
Lyapunov2100.111 pu88.9%
Lyapunov2500.031 pu96.9%
Lyapunov550.139 pu86.1%
Lyapunov5100.079 pu92.1%
Lyapunov5500.035 pu96.5%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gezer, D.; Taşcıoğlu, Y.; Çelebioğlu, K. Frequency Containment Control of Hydropower Plants Using Different Adaptive Methods. Energies 2021, 14, 2082. https://doi.org/10.3390/en14082082

AMA Style

Gezer D, Taşcıoğlu Y, Çelebioğlu K. Frequency Containment Control of Hydropower Plants Using Different Adaptive Methods. Energies. 2021; 14(8):2082. https://doi.org/10.3390/en14082082

Chicago/Turabian Style

Gezer, Doğan, Yiğit Taşcıoğlu, and Kutay Çelebioğlu. 2021. "Frequency Containment Control of Hydropower Plants Using Different Adaptive Methods" Energies 14, no. 8: 2082. https://doi.org/10.3390/en14082082

APA Style

Gezer, D., Taşcıoğlu, Y., & Çelebioğlu, K. (2021). Frequency Containment Control of Hydropower Plants Using Different Adaptive Methods. Energies, 14(8), 2082. https://doi.org/10.3390/en14082082

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