^{1}

^{1}

^{1}

The interest on the use of renewable energy resources is increasing, especially towards wind and hydro powers, which should be efficiently converted into electric energy via suitable technology tools. To this aim, self-tuning control techniques represent viable strategies that can be employed for this purpose, due to the features of these nonlinear dynamic processes working over a wide range of operating conditions, driven by stochastic inputs, excitations and disturbances. Some of the considered methods were already verified on wind turbine systems, and important advantages may thus derive from the appropriate implementation of the same control schemes for hydroelectric plants. This represents the key point of the work, which provides some guidelines on the design and the application of these control strategies to these energy conversion systems. In fact, it seems that investigations related with both wind and hydraulic energies present a reduced number of common aspects, thus leading to little exchange and share of possible common points. This consideration is particularly valid with reference to the more established wind area when compared to hydroelectric systems. In this way, this work recalls the models of wind turbine and hydroelectric system, and investigates the application of different control solutions. Another important point of this investigation regards the analysis of the exploited benchmark models, their control objectives, and the development of the control solutions. The working conditions of these energy conversion systems will also be taken into account in order to highlight the reliability and robustness characteristics of the developed control strategies, especially interesting for remote and relatively inaccessible location of many installations.

The trend to reduce the use of fossil fuels, motivated by the need to meet greenhouse gas emission limits, has driven much interest on renewable energy resources, in order also to cover global energy requirements. Wind turbine systems, which now represent a mature technology, have had much more development with respect to other energy conversion systems, e.g. for biomass, solar, and hydropower [

One common aspect regarding the design of the renewable energy conversion system concerns the optimality and the efficiency of its converter. However, as wind and hydraulic resources are free, the key point is represented by the minimisation of the cost per kWh, also considering the lifetime of the deployments. Moreover, when the capital, operational and commissioning/decommissioning costs are fixed, the value of the energy sold (i.e. the energy receipts) has to be maximised. This represents a fundamental economic objective, which should be carefully taken into account also by the design of the control system.

It is worth noting that preliminary works highlighted interactions between the renewable energy conversion system design and the strategies exploited to control them, see for example, [

It is quite established that the mathematical descriptions of both wind turbine processes and hydroelectric plants are represented by nonlinear dynamic processes working over a wide range of operating conditions and excitations. These systems are also required to operate under specific physical constraints, such as displacements, velocities, accelerations, torques and forces. Therefore, these systems can operate effectively with economically attractive and high operational lifetimes only if their working conditions are carefully fulfilled. With these issues in mind, the paper recalls the mathematical description of a wind turbine system and a hydroelectric plant, by using a wind turbine benchmark and a hydroelectric simulator, respectively. The former process was proposed for the purpose of an international competition started in 2009 and described in detail [

After these remarks, by means of the analysis of the proposed modelling and control topics, the work will sketch common and different aspects of wind turbines and hydraulic systems, which will be exploited for the design of the control technique. On one hand, hydroelectric power plants result to be more established and even more common than wind turbine processes, but the control aspects for the latter systems have been examined more in depth in the last decades.

With reference to wind turbine systems, it can be observed that modern installations exploit control techniques and technologies in order to obtain the needed goals and performance achievements. These plants can implement their regulation via “passive” control methods, such as the plants with fixed-pitch, and stall control machines. The blades of these systems are deployed in order to limit their power via the blade stall, when the wind speed exceeds its rated value. These systems do not need any pitch control mechanism, as addressed e.g. in [

It is worth noting that a limited number of works have addressed the development of self-tuning control techniques for hydroelectric plants, as shown e.g. in [

In the light of these considerations, the paper also proposes to analyse those control aspects that might be similar between wind turbine and hydroelectric systems, with the aim of exploiting some solutions, developed in the wind turbine domain, and to apply them within the other concerning hydroplants. This approach could be used to stimulate novel research topics and the development of innovative techniques in a multidisciplinary control community, and the most important achievements will be summarised in this paper. In particular, suitable analytical models of these energy conversion plants should be able to provide the overall dynamic behaviour of the monitored processes, thus leading important impacts on the development of the control techniques. Moreover, the work introduces some kind of common rules for tuning the different controllers, for both wind turbine and hydroelectric plants. Therefore, the paper shows that the parameters of these controllers are obtained by exploiting the same tuning strategies. This represents one of the key contributions of the study.

Note that some previous studies by the same authors addressed several topics presented in this paper. For example, the work [

It is worth highlighting the main contribution of the paper, which aims at providing some guidelines on the design and the application of self-tuning control strategies to two energy conversion systems. Some of these techniques were already verified on wind turbine systems, and important advantages may thus derive from the appropriate implementation of the same control methods for hydroelectric plants. In fact, it seems that investigations related with both wind and hydraulic energies present a reduced number of common aspects, thus leading to little exchange and share of possible common points. This consideration is particularly valid with reference to the more established wind area when compared to hydroelectric systems. In this way, the paper summarises also the most common models used for describing wind turbine and hydroelectric systems. Moreover, it analyses the application of the different control solutions to these energy conversion systems. The aim is thus to exploit common points in the control objectives and the achievable results from the application of different solutions. The working conditions of these energy conversion systems will be also taken into account in order to highlight the reliability and robustness characteristics of the developed control strategies.

Finally, the paper has the following structure. Section 2 provides the brief presentation of the benchmark and simulation models used for describing the accurate behaviour of the dynamic processes. Section 3 discusses the specific requirements of the control systems exploited to control these energy conversion plants. Section 3 summarises the design of the proposed model-based and data-driven control techniques, taking into account the control objectives and the available tools. In Section 4, these self-tuning control strategies are implemented and compared, with respect to the achievable reliability and robustness features. Section 5 ends the paper summarising the main achievements of the paper, and drawing some concluding remarks.

This section recalls the simulators used for describing the dynamic behaviour of the wind turbine and the hydroelectric processes considered in this paper. Moreover, the baseline control schemes developed for the regulation of the wind turbine benchmark are also summarised in Section 2.1. On the other hand, the hydroelectric simulator, together with its reference governor, is recalled in Section 2.2.

Industrial wind turbine installations are normally equipped with large rotors, flexible blades and light load-carrying structures, which work in uncertain environments, often placed in remote and inaccessible places [

In particular regarding wind turbines, this work focuses on a horizontal-axis device, which nowadays represents the most common type of installation for large-scale deployments. Moreover, the three-bladed horizontal axis wind turbine model reported in this work follows the principle that the wind power activates the wind turbine blades, thus producing the rotation of the low speed rotor shaft. In order to increase its rotational speed generally required by the generator, a gear-box with a drive-train is included in the system. A more detailed description of this benchmark that was proposed for different purposes is provided in [

The wind turbine simulator presents 2 controlled outputs, i.e. the generator rotational speed ω g ( t ) and its generated power P g ( t ) . The wind turbine model is controlled by means of two actuated inputs, i.e. the generator torque τ g ( t ) and the blade pitch angle β ( t ) . The latter signal controls the blade actuators, which can be implemented by hydraulic or electric drives. The benchmark considered in this work includes a hydraulic circuit actuating the

wind turbine blades [

Several other measurements are acquired from the wind turbine benchmark: the signal ω r ( t ) represents the rotor speed and τ r ( t ) is the reference torque. Moreover, aerodynamic torque signal τ a e r o ( t ) is computed from the wind speed v ( t ) , which is usually available with limited accuracy. In fact, the wind field is not uniform around the wind turbine rotor plane, especially for large rotor systems. Moreover, anemometers measuring this variable are mounted behind the rotor on the nacelle. Therefore, the wind speed measurement v ( t ) is affected by the interference between the blades and the nacelle, as well as the turbulence around the rotor plane. Furthermore, when these instantaneous wind fields are considered across the rotor plane, the wind variable v ( t ) may change in space and time, and it is especially true in large rotor installations. The alteration of the wind speed measurement v ( t ) with respect to its nominal value around the rotor plane represents an uncertainty in the wind turbine model and a disturbance term in control design. On the other hand, the aerodynamic torque depends on another factor, C p , representing the power coefficient, as shown by Equation (1):

τ a e r o ( t ) = ρ A C p ( β ( t ) , λ ( t ) ) v 3 ( t ) 2 ω r ( t ) (1)

ρ being air density, A the area swept by the turbine blades during their rotation, whilst λ ( t ) represents an important variable, i.e. the tip-speed ratio of the blade, which is given by the relation of Equation (2):

λ ( t ) = ω r ( t ) R v ( t ) (2)

where R is the rotor radius. The nonlinear static function C p ( ⋅ ) represents the power coefficient, which is usually modelled via a two-dimensional map (or look-up table) [

It is worth noting that diagram in

The wind turbine benchmark considered in this work includes a simple two-body linear model of the third order that is exploited to describe the dynamic behaviour of the drive-train. It implements also a simple first-order linear dynamic model of the electric generator and a second order dynamic description of the pitching system, as addressed in more detail in [

continuous-time representation of the wind turbine benchmark is represented by the general model in form of Equation (3):

{ x ˙ ( t ) = f c ( x ( t ) , u ( t ) ) y ( t ) = x ( t ) (3)

with u ( t ) = [ τ r ( t ) β ( t ) ] T and y ( t ) = [ ω g ( t ) P g ( t ) ] T are the manipulated input signals and the controlled output measurements, respectively. f c ( ⋅ ) is described by means of a continuous-time nonlinear function that will be exploited for representing the complete dynamic behaviour of the controlled process. Moreover, since this paper will analyse several data-driven control approaches, this system will be used to acquire a number of N sampled data sequences u ( k ) and y ( k ) , with k = 1 , 2 , ⋯ , N . Furthermore, the variables and parameters of the wind turbine benchmark submodels (see e.g.

As already highlighted by

τ r = K o p t ω r 2 (4)

In this situation, the optimal tracking of the power reference is obtained, as soon as the wind speed v ( t ) increases, and the working condition moves to the control region 3. The control task aims also at tracking the power reference P r , which is achieved by modifying β , while C p is decreasing. The advanced control strategies considered in this work tries to maintain the generator speed ω g at its nominal value ω n o m by changing both β and τ g .

Therefore, the control system operating in region 2 exploits the relations in the form of Equations (5) when implemented as difference equations [

{ β ( k ) = β ( k − 1 ) + k p e ( k ) + ( k i T s − k p ) e ( k − 1 ) e ( k ) = ω g ( k ) − ω n o m (5)

where k = 1 , 2 , ⋯ , N corresponds to the sample indices, and the variable ω n o m is the given reference generator speed, depending on the wind turbine plant. For the case of the wind turbine system considered in this work, P r = 4.8 MW is the rated power, and ω n o m = 162.5 rad / s . The standard PI governor parameters used for the speed control task were settled to k i = 0.5 and k p = 3 , with

Variable | v ( t ) | ω r | ω g | τ g | P g | β |
---|---|---|---|---|---|---|

Std. Dev. Value | 0.5 m/s | 0.025 rad/s | 0.05 rad/s | 90 Nm | 10^{3} W | 0.2 deg |

sampling time T s = 0.01 s [

Concerning the regulation of the second input τ g , a further standard PI governor is implemented in the wind turbine benchmark, similarly to the one of Equation (5), which is described again in its discrete-time formulation of Equations (6):

{ τ r ( k ) = τ r ( k − 1 ) + k p e ( k ) + ( k i T s − k p ) e ( k − 1 ) e ( k ) = P g ( k ) − P r (6)

This standard PI regulator exploited in the benchmark for the power control task has its parameters settled to k i = 0.014 and k p = 447 × 10 − 6 , as proposed in [

Finally, Section 4 will consider the performances of these baseline controllers summarised by the overall laws of Equations (4)-(6) proposed in [

It is well-established that hydroelectric systems transform hydraulic renewable source into useful energies, mostly electric but also mechanical one. However, as for wind turbines, they must operate according to different load situations. In general, hydroelectric plants must operate despite of possible variations in the hydropower flow, and in particular in both planned or nominal conditions and accidental or unplanned situations. Moreover, routine operations such as start-up, shutdown, load rejection and acceptance may induce important hydraulic transients, possible leading to dangerous high pressure and sub-pressure variations and oscillations in the hydraulic system. These situations must be analysed in order to avoid possible mechanical malfunctions and failures. The same simulation codes already exploited for the development of the wind turbine benchmark described in Section 2.1, i.e. Matlab and Simulink, are the tools exploited for modelling, simulating, and analysing the behaviour of hydroelectric plants that exhibit important nonlinear dynamics. Therefore, hydropower plants have to include special control techniques to guarantee stable and safe working conditions. The same self-tuning control methodologies already developed for wind turbine systems, as summarised in Section 3, will be thus considered for the hydroelectric process described in this paper.

With reference to the hydroelectric system, which is recalled in this work for analysis and comparison purposes, consists of a high water head and a long penstock. It includes also upstream and downstream surge tanks, where a Francis hydraulic turbine is included [

The scheme of this hydroelectric simulator including two surge tanks and the Francis hydraulic turbine considered in this work is recalled in

The hydraulic system considered in this paper was modified by the authors in order to include the Francis hydraulic turbine, as presented in [

h q = − T w s − H f (7)

Moreover, Equation (7) represents the transfer function between the flow rate deviation and the water pressure deviation valid for a simple penstock. The variable h represents the water pressure relative deviation, whilst q is the flow rate relative deviation. The term H f represents the hydraulic loss, with s the Laplace operator, and T w the water inertia time expressed by the relation of Equation (8):

T w = L Q r g A H r (8)

Note that the time variable T w of the water inertia described by the relation of Equation (8) is a function of the hydraulic variable, such as the penstock length L, the rated flow rate Q r , the gravity acceleration g, the cross-section area

A, and the rated water pressure H r . The classic plant represented in

h 1 q 1 = − T w 1 s − H f 1 (9)

On the other hand, the downstream tail water tunnel connects the downstream surge tank with the tail water lake. The downstream tail water tunnel outlet is assumed to coincide with the tailwater lake, with constant outlet water pressure deviation. In this way, the transfer function between the flow rate and the water pressure deviations of the downstream tail water tunnel inlet is represented in the form of Equation (10):

h 5 q 5 = − T w 5 s − H f 5 (10)

Usually, the draft tube water inertia is considered within the penstock. Therefore, the transfer function between the flow rate and the water pressure deviations within the penstock are expressed by the relation of Equation (11):

h t = h 2 − h 4 + h 3 (11)

with:

h 3 q 3 = − T w 3 s − H f 3 (12)

The relations describing the surge tanks are formulated from the flow continuity at the two junctions, by neglecting the hydraulic losses at surge tank orifices, and represented via the relations of Equations (13):

{ A 2 H r Q r d h 2 d t = q 2 = q 1 − q 3 A 4 H r Q r d h 4 d t = q 4 = q 3 − q 5 (13)

In this situation, the surge tank filling time has the form of Equation (14):

T s = A H r Q r (14)

The mathematical model and the performance curves of the Francis turbine considered in this work were obtained in order to describe the dynamic behaviour of a realistic hydroelectric process. To this aim, the values of the most important variables the hydraulic system and the Francis hydraulic turbine, which represent the overall hydroelectric process simulator working at rated conditions are summarised in

Note that, with reference to the values summarised in

After these considerations, in the following the procedure for computing the non-dimensional performance curves of the hydraulic turbine considered in this work is briefly recalled. In particular, the non-dimensional water flow rate Q / Q r is expressed as a function of the non-dimensional rotational speed n / n r , and represented by the second order polynomial of Equation (15):

Q Q r = G [ a 1 ( n n r ) 2 + b 1 ( n n r ) + c 1 ] = f 1 ( n , G ) (15)

Moreover, the relation of Equation (15) includes the wicked gate opening, described by the non-dimensional parameter G, varying from 0 to 100%. In particular,

The hydroelectric simulator assumes that the turbine efficiency is constant and equal to its rated value η r , i.e. 0.9, as reported in

On the other hand, the non-dimensional turbine torque M results a function of the water flow rate Q, the water level H and the rotational speed n, as

Variable | Description | Value |
---|---|---|

H r | Reservoir water level | 400 m |

Q r | Water flow rate | 36.13 m^{3}/s |

P r | Hydraulic turbine power | 127.6 MW |

n r | Turbine rated rotational speed | 500 rpm |

η r | Efficiency rated value | 0.90 |

M r | Turbine-rated torque | 2437 kNm |

highlighted by the relation of Equation (16):

M M r = Q Q r H H r n n r = f 2 ( Q , n , G ) (16)

Moreover, the combination of the relations of Equations (15) and (16) highlights that the turbine torque M is a function of the water flow rate Q, the rotational speed n and wicked gate opening G.

Finally, the overall model of the hydroelectric simulator is described by the relations of Equations (17)-(20), which express the non-dimensional variables with respect to their relative deviations:

Q Q r = 1 + q t (17)

H H r = 1 + h t (18)

n n r = 1 + x (19)

G = 1 + y (20)

with q t is the turbine flow rate relative deviation, h t the turbine water pressure relative deviation, x the turbine speed relative deviation, and y the wicket gate servomotor stroke relative deviation. Moreover, the relation of Equation (20) allows only negative values of y.

On the other hand, when the generator unit and its network are considered, and in particular the generator unit is connected only to an isolated load, the load characteristic of the generator unit is described by the dynamic model of Equation (21):

x m t − m g 0 = 1 T a s + e g (21)

with m g 0 being the load torque, T a representing the generator unit mechanical time, whilst the parameter e g is the load self-regulation factor. The variables and parameters of the hydroelectric model were selected according to the work [

With reference to the control strategies for classic hydroelectric plants, standard PID regulators are used to compensate the hydraulic turbine speed. Therefore, the actuated signal u is computed as sum of the proportional, integral, and differential terms of the error x in Equation (19), expressed in the form of Equation (22):

u = x ( K p + K i s + K d s 1 + T n s ) (22)

with K p being the proportional gain, K i the integral gain, and K d the derivative gain. T n is the parameter of the derivative filter time constant. The hydroelectric simulator considered in this work exploits an electric servomotor that is used as a governor.

The servomechanism implemented in the hydroelectric simulator is described as a first-order model, which relates the control signal u with the wicket gate servomotor stroke y according to Equation (23):

y u = 1 T y s + 1 (23)

with T y representing the wicket gate servomotor response time.

This concludes the description of the complete nonlinear simulator of a typical hydroelectric plant consisting of two surge tanks and a Francis hydraulic turbine, as represented in

Finally, it is worth noting that some relations of the hydroelectric system have been linearised, see e.g. Equations (7) and (17). However, this simplified model has been considered for comparison purpose, as the nonlinear parts of the processes under investigation are closer, as highlighted by Equations (1) and (15).

This section recalls the self-tuning control methodologies that will be designed and compared when applied to the considered energy conversion benchmark and simulator.

In general, control systems exploit design algorithms that force a dynamic model to track prescribed references or set-point, such that fixed objectives or behaviour modes are achieved. In this way, the classic control problem is formulated as tracking task, where the system output has to follow the set-point, thus representing the final objective. Tasks expressed in this form are also present in energy conversion applications, for example the speed control of both wind and hydraulic turbines. However, it could be useful to improve the problem descriptions and to give a deeper insight into possible solutions, in order to achieve all potentials of control theory when applied to energy conversion systems.

On one hand, data-driven techniques rely on the availability of the input-output data acquired from the monitored plant. These data are use for the on-line estimation of a suitable model of the dynamic process, which is thus exploited for the identification of the control law to be applied to the controlled system. On the other hand, model-based methods require the mathematical description of the monitored plant, in the form of input-output or state-space representations. These forms are thus employed for the analytical derivation of the mathematical function of the controller, again in form of input-output or state-space relations. Both the derived controllers are obtained in order to reach prescribed performances.

First, with reference to the process output, the desired transient or steady-state responses can be considered, as for the case of self-tuning PID regulators summarised in Section 3.1. On the other hand, if the frequency behaviour is taken into account, the desired closed-loop poles can be fixed as roots of the closed-loop transfer function. This represents the design approach used by the adaptive strategy considered in Section 3.3. Moreover, when robust performances are included, the minimisation of the sensitivity of the closed-loop system with respect to the model-reality mismatch or external disturbances can be considered. This approach is related for example to the fuzzy logic methodology

reported in Section 3.2. Some other strategies provide solutions to this optimisation problem when it is defined at each time step, as for the case of the Model Predictive Control (MPC) with disturbance decoupling considered in Section 3.4. The considered strategy integrates the advantages of the MPC solution with the disturbance compensation feature.

It is worth noting that model-based control designs rely on the mathematical descriptions of the process models, in order to derive the control laws. The design of standard PID regulators and Model Predictive Control (MPC) methods follows a model-based approach, which will be illustrated in Sections 3.1 and 3.4, respectively. However, the need for these high-fidelity mathematical descriptions can require much more effort than the derivation of the controller models. Therefore, dynamic system identification methodologies have been successfully proposed in order to determine the so-called black-box representations, which is also used for the self-tuning PID design. Usually, these descriptions do not present structural relationships to the physical processes. On the other hand, dynamic system identification schemes can be also exploited for deriving adaptive controller prototypes, which are thus able to adapt themselves with respect to unknown conditions or time-varying systems. By means of this “self-tuning mode”, adaptive control strategies relying on linear models of the controlled process are also able to track changes of the plant. Examples of these data-driven approaches are represented by the fuzzy logic and adaptive controllers recalled in Sections 3.2 and 3.3, respectively. On the other hand, the MPC strategy exploits the proposed disturbance compensation method, which is thus able to cope with uncertainty and model-reality mismatch effects.

In general, the mathematical formulation of the control law can be provided as linear or nonlinear dynamic function F in the form of Equation (24):

u ( t ) = F ( y ( t ) ) (24)

with y ( t ) being the monitored output, whilst u ( t ) is the control input. The control techniques proposed for the systems under investigation should lead to the computation of the control law of Equation (24) generating the input u ( t ) that allows to track the given reference or set-point r ( t ) for the controlled output y ( t ) .

Industrial processes commonly exploit closed-loop including standard PID controllers, due to their simple structure and parameter tuning [

u ( t ) = K p e ( t ) + K i ∫ 0 t e ( τ ) d τ + K d d e ( t ) d t (25)

with K p , K i , K d being the PID proportional, integral, and derivative gains, respectively. The most common strategy exploited for the computation of the optimal parameters of the PID governor uses proper Ziegler-Nichols formulas, as described in [

Note finally that the PID block in

Fuzzy Logic Control (FLC) solutions are often exploited when the dynamics of the monitored process are uncertain and can present nonlinear characteristics. The design method proposed in this work exploits the direct identification of rule-based Takagi-Sugeno (TS) fuzzy prototypes. Moreover, the fuzzy model structure, i.e. the number of rules, the antecedents, the consequents and the fuzzy membership functions can be estimated by means of the Adaptive Neuro-Fuzzy Inference System (ANFIS) toolbox implemented in the Simulink environment [

The TS fuzzy prototype relies on a number of rules R i , whose consequents are deterministic functions f i ( ⋅ ) in the form of Equation (26):

R i : IF x is A i THEN u i = f i ( x ) (26)

where the index i = 1 , 2 , ⋯ , K describes the number of rules K, x is the input vector containing the antecedent variables, i.e. the model inputs, whilst u i represents the consequent output. The fuzzy set A i describing the antecedents in the i-th rule is described by its (multivariable) membership function μ A i ( x ) → [ 0,1 ] . The relation f i ( x ) assumes the form of parametric affine model represented by the i-th relation of Equation (27):

u i = a i T x + b i (27)

with the vector a i and the scalar b i being the i-th submodel parameters. The vector x consists of a proper number n of delayed samples of input and output signals acquired from the monitored process. Therefore, the term a i T x is an Auto-Regressive eXogenous (ARX) parametric dynamic model of order n, and b i a bias.

The output u of the TS fuzzy prototype is computed as weighted average of all rule outputs u i in the form of Equation (28):

u = ∑ i = 1 K μ A i ( x ) y i ( x ) ∑ i = 1 K μ A i ( x ) (28)

The estimation scheme implemented by the ANFIS tool follows the classic dynamic system identification experiment. First, the structure of the TS fuzzy prototype is defined by selecting a suitable order n, the shape representing the membership functions μ A i , and the proper number of clusters K. Therefore, the input-output data sequences acquired from the monitored system are exploited by ANFIS for estimating the TS model parameters and its rules R i after the selection of a suitable error criterion. The optimal values of the controller parameters represented by the variables a i and b i of (27) are thus estimated [

The work proposes also a strategy different from ANFIS that can be exploited for the estimation of the parameters of the fuzzy controller. This method relies on the Fuzzy Modelling and Identification (FMID) toolbox designed in the Matlab and Simulink environments as described in [

Note that the overall digital control scheme consisting of the discrete-time fuzzy regulator of Equation (28) and continuous-time nonlinear system of Equation (24) includes also Digital-to-Analog (D/A) and Analog-to-Digital (A/D) converters, as shown in

With reference to

The adaptive control technique proposed in this work relies on the recursive estimation of a 2nd order discrete-time transfer function G ( z ) with time-varying parameters described by Equation (29):

G ( z ) = β 1 z − 1 + β 2 z − 2 1 + α 1 z − 1 + α 2 z − 2 (29)

where α i and β i are identified on-line at each sampling time t k = k T , with k = 1 , 2 , ⋯ , N , for N samples, and T being the sampling interval. z − 1 indicates the unit delay operator. A viable and direct way for deriving the model parameters in Equation (29) that is proposed in this work is based on the Recursive Least-Square Method (RLSM) with directional forgetting factor, which was presented in [

Once the parameters of the model of Equation (29) have been derived, this paper proposes to compute the adaptive controller in the form of Equation (30):

u k = q 0 e k + q 1 e k − 1 + q 2 e k − 2 + ( 1 − γ ) u k − 1 + γ u k − 2 (30)

with e k and u k represent the sampled values of the tracking error e ( t ) and the control signal u k at the time t k , respectively. With reference to the description

of Equation (30), by following a modified Ziegler-Nichols criterion, q 0 , q 1 , q 2 , and γ represent the adaptive controller parameters, which are derived by solving a Diophantine equation. As described in [

{ q 0 = 1 β 1 ( d 1 + 1 − α 1 − γ ) γ = s 1 r 1 β 2 α 2 q 1 = α 2 β 2 − s 1 r 1 ( β 1 β 2 − α 1 α 2 + 1 ) q 2 = s 1 r 1 (31)

where:

{ r 1 = ( b 1 + b 2 ) ( a 1 b 2 b 1 − a 2 b 1 2 − b 2 2 ) s 1 = a 2 ( ( b 1 + b 2 ) ( a 1 b 2 − a 2 b 1 ) + b 2 ( b 1 d 2 − b 2 d 1 − b 2 ) ) (32)

Note that the design technique proposed in this work and represented by the relations of Equations (31) and (32) assumes that the behaviour of the overall closed-loop system can be approximated by a 2^{nd} order transfer function with characteristic polynomial represented by Equation (33):

D ( s ) = s 2 + 2 δ ω s + ω 2 (33)

with δ and ω being the damping factor and natural frequency, respectively. s is the derivative operator. Furthermore, if δ ≤ 1 , the following relations are used [

{ d 1 = − 2 e − δ ω T cos ( ω T 1 − δ 2 ) d 2 = e − 2 δ ω T (34)

This paper suggested this adaptive control technique since both the recursive estimation procedure of Equation (29) and the on-line computation of the control law of Equation (30) are available from the digital Self-Tuning Controller Simulink Library (STCSL) described in [

The on-line control law of Equation (30) is used for the regulation of the continuous-time nonlinear system of Equation (24) by including D/A and A/D converters, as highlighted in the scheme of

Note finally that the adaptive control sketched in

The general structure of the proposed MPC is illustrated in

MPC design itself.

The overall scheme is thus represented aim by the MPC design with disturbance compensation, such that the compensated system has response very similar to the nominal system and the constraints are not violated. The fault compensation problem within the MPC framework is defined as follows. Given a state-space representation of the considered system affected by disturbance or uncertainty has the following form:

{ x ˙ l = A l x l + B l u + B d d + w y l = C l x l + v (35)

and its nominal reference model:

{ x ˙ r = A l x r + B l u r y r = C l x r (36)

the disturbance compensation problem is solved by finding the control input u that minimises the cost function:

J = ∫ t t + N c Δ t ( ‖ x l − x r ‖ Q 2 + ‖ u ˙ ‖ R 2 ) d τ (37)

given the reference input u r .

In Equation (35) the matrices A l , B l , B d and C l are of proper dimensions. The vector y l represents the output measurements, x l is the state of the model with disturbance, whilst x r is the reference state, and y r the reference output, corresponding to the reference input u r of the nominal model. The vectors w and v include the model mismatch and the measurement error, respectively. d represents the equivalent disturbance signal. In Equation (37) t is the current time, Δ t is the control interval, and N c is the length of the control horizon. Q and R are suitable weighting matrices. Note that the model of Equation (35) can be derived by nonlinear model linearisation or identification procedures, as suggested in Sections 3.1 and 3.3, respectively.

This work proposes to solve the problem in two steps: the reconstruction of the disturbance d, i.e. d ^ , provided by the disturbance estimation module, and the MPC tool. Due to the model-reality mismatch and the measurement error in (35), the Kalman filter (38) is used to provide the estimation of the state vector x l , the output y l of the system affected by the estimated disturbance d ^ :

{ x ˙ l = A l x l + B l u − B l d ^ + K f ( y l − C l x l ) y l = C l x l (38)

where K f is the Kalman filter gain. In this way, based on the estimations d ^ and x l , an MPC is designed, which contains the reference model of Equation (36) and the filtered system of Equation (38), with d ^ provided by the Kalman filter. Moreover, the MPC has the objective function:

∫ t t + N c Δ t [ ( x l − x r ) T Q ( x l − x r ) + u ˙ T R u ˙ ] d τ (39)

in which x l and x r are the states of the filtered and the reference models, respectively. The integrated MPC with the Kalman filter solves this general disturbance compensation problem, as long as the estimations of both the state and the disturbance are correct. An illustration of the structure of the fault compensated MPC is shown in

The global estimation and control scheme is a nonlinear MPC problem with the nominal model for the considered energy conversion systems of Equation (35), the disturbance d with its estimator, and the Kalman filter of Equation (38) as prediction model. The local observability of the model of Equation (35) is essential for state estimation, which is easily verified. The implementation of the proposed disturbance compensation strategy has been integrated into the MPC Toolbox of the Simulink environment.

This section presents the simulations achieved in the Matlab and Simulink environments implementing the control techniques and tools recalled in Section 3. The obtained results are evaluated via the percent Normalised Sum of Squared Error ( N S S E % ) performance function in the form of Equation (40):

N S S E % = 100 ∑ k = 1 N ( r k − o k ) 2 ∑ k = 1 N r k 2 (40)

with r k being the sampled reference or set-point r ( t ) , whilst o k is the sampled continuous-time signal representing the generic controlled output y ( t ) of the process. In particular, this signal is represented by the wind turbine generator angular velocity ω g ( t ) in Equation (3), and the hydraulic turbine rotational speed n in Equation (19) for the hydroelectric plant.

Note that the wind turbine benchmark and the hydroelectric plant simulator of Section 2 allow the generation of several input-output data sequences due to different wind speed v ( t ) effects (see e.g. (1)) and hydraulic transient under variable loads m g 0 (see e.g. (21)), respectively. Moreover, in order to obtain comparable working situations, the wind turbine benchmark has been operating from partial to full load conditions, as highlighted in

It is worth highlighting that the simulations considered in this work take into account disturbance and uncertainty effects. In fact, the hydroelectric plant considers a load disturbance, whilst the turbine simulator is driven by wind, which represents the main disturbance source. Moreover, the uncertainty effect has been analysed in Section 4.2.

In detail, with reference to the picture in

Moreover,

On the other hand, the picture in

Finally, the picture of

horizon N c = 2 . The weighting factors have been settled to w y k = 0.1 and w u k = 1 , in order to reduce possible abrupt changes of the control input. Note that, in this case, the MPC technique has led to the best results, since it exploits a disturbance decoupling strategy, whilst its parameters have been iteratively adapted in the Simulink environment in order to optimise the MPC cost function of Equation (37), as addressed in Section 3.4.

The second test case regards the hydroelectric plant simulator, where the hydraulic system with its turbine speed governor generates hydraulic transients due to the load changes. As already recalled in Section 2.2, an effective behaviour of a classic PID governor addressed e.g. in [

Under these assumptions,

In more detail,

x = [ e k , e k − 1 , e k − 2 , e k − 3 , u k − 1 , u k − 2 , u k − 3 ] . Moreover, as recalled in Section 3.2, the data-driven FMID and ANFIS tools implemented in the Simulink toolboxes are able to provide the estimates of the shapes of the membership functions μ A i used in Equation (28).

On the other hand,

Finally, regarding the MPC technique with disturbance decoupling proposed in Section 3.4,

After these considerations, it is worth noting that some of the control techniques recalled in this paper rely on self-tuning and adaptive methodologies, that are based on data-driven algorithms. This means that they do not need for the knowledge of a high-fidelity description of the controlled process, thus providing a viable and direct implementation.

In order to provide a quantitative comparison of the tracking capabilities obtained by the considered control techniques for the wind turbine benchmark, the first row in

In particular, the NSSE% values in the first row of

Simulated system | Wind Turbine | Hydro Plant |
---|---|---|

Working condition | From partial to full load | From start-up to full load |

Standard PID | 11.5% | 6.2% |

Autotuning PID | 7.3% | 4.9% |

Fuzzy Controller | 5.7% | 3.1% |

Adaptive Controller | 4.1% | 1.8% |

MPC Scheme | 2.8% | 0.9% |

technique with disturbance decoupling has led to the best results, as reported in the first row of

On the other hand, the results achieved by the validation of the considered control techniques to the hydroelectric plant simulator are summarised in the second row of

With reference again to

Finally, in order to highlight some further features of the considered, the controlled inputs applied to the wind turbine system are depicted and compared in

By considering these control inputs, with reference to the data-driven methodologies, and in particular to the design of the fuzzy controllers, off-line optimisation strategies allow to reach quite good results. However, control inputs are subjected to faster variations. Other control techniques can take advantage of more complicated and not direct design methodologies, as highlighted by the MPC scheme. In this case, due to the input constraint, its changes are reduced. This feature is attractive for wind turbine systems, where variations of the control inputs must be limited. This represents another important benefit of MPC with disturbance decoupling, which integrates the advantages of the classic MPC scheme with disturbance compensation effects. Therefore, with reference to these two control methods, they can appear rather straightforward, even if further optimisation and estimation strategies have to be applied.

This section analyses the robustness properties of the developed controllers when parameter variations and measurement errors are considered. This further investigation relies on the Monte-Carlo tool, since the control behaviour and the tracking capabilities depend on both the model-reality mismatch effects and the input-output uncertainty levels. Therefore, this analysis has been implemented by describing the parameters of both the wind turbine system and hydroelectric plant models as Gaussian stochastic processes with average values corresponding to the nominal ones summarised in

Moreover,

On the other hand,

Therefore, the average values of NSSE% index have been thus evaluated by means of 1000 Monte-Carlo simulations. They have been reported in

Variable | R | χ |
---|---|---|

Nominal value | 57.5 m | 0.6 |

Variable | ω n | B d t |

Nominal value | 11.11 rad s^{-1} | 775.49 N m s rad^{-1} |

Variable | B r | B g |

Nominal value | 7.11 N m s rad^{-1} | 45.6 N m s rad^{-1} |

Variable | K d t | η d t |

Nominal value | 2.7×10^{9} N m rad^{-1} | 0.97 |

Variable | J g | J r |

Nominal value | 390 kg m^{2} | 55×10^{6} kg m^{2} |

Variable | a | b | c |
---|---|---|---|

Nominal value | −0.08 | 0.14 | 0.94 |

Variable | T c | T s 2 | T s 4 |

Nominal value | 20 s | 476.05 s | 5000 s |

Variable | H f 3 | H f 5 | T a |

Nominal value | 0.0481 m | 0.0047 m | 5.9 s |

Variable | T w 3 | T w 5 | H f 1 |

Nominal value | 0.83 s | 0.1 s | 0.0481 m |

Standard PID | Autotuning PID | Fuzzy PID | Adaptive PID | MPC Scheme |
---|---|---|---|---|

13.8% | 9.2% | 7.6% | 5.3% | 3.9% |

Standard PID | Autotuning PID | Fuzzy PID | Adaptive PID | MPC Scheme |
---|---|---|---|---|

9.1% | 7.4% | 5.6% | 3.5% | 2.2% |

simulator, respectively.

It is worth noting that the results summarised in

The work considered two renewable energy conversion systems, such as a wind turbine benchmark and a hydroelectric plant simulator. The most important modelling aspects and the baseline control strategies were also summarised. In particular, the three-bladed horizontal axis wind turbine benchmark reported in this work consisted of simple models of the gear-box, the drive-train, and the electric generator/converter. On the other hand, the hydroelectric plant simulator included a high water head, a long penstock with upstream and downstream surge tanks, and a Francis hydraulic turbine. Standard PID governors were earlier developed for these processes, which were rather simple and straightforward, but with limited achievable performances. Therefore, the paper proposed different control techniques relying on model-based and data-driven approaches. Their performances were analysed first. Then, the robustness characteristics of these solutions were also verified and validated with respect to parameter variations of the plant models and measurement errors, via the Monte-Carlo tool. The achieved results highlighted that data-driven approaches, such as the fuzzy regulators were able to provide good tracking performances. However, they were easily outperformed by adaptive and model predictive control schemes, representing data-driven and model-based approaches that require optimisation stages, adaptation procedures and disturbance compensation methods. Future investigations will consider the verification and the validation of the considered control techniques when applied to higher fidelity simulators of energy conversion systems.

The authors declare no conflicts of interest regarding the publication of this paper.

Simani, S., Alvisi, S. and Venturini, M. (2019) Self-Tuning Control Techniques for Wind Turbine and Hydroelectric Plant Systems. Journal of Power and Energy Engineering, 7, 27-61. https://doi.org/10.4236/jpee.2019.71003