## Servicios Personalizados

## Revista

## Articulo

## Indicadores

- Citado por SciELO

## Links relacionados

- Similares en SciELO

## Compartir

## Latin American applied research

##
*versión impresa* ISSN 0327-0793

### Lat. Am. appl. res. v.33 n.4 Bahía Blanca oct./dic. 2003

**Identification of partially known models of the Susqueda hidroelectric power plant**

**O. Quiroga ^{1}, J. Riera^{2} and C. Batlle^{3}**

^{1} *Dep de Ing. Ind. (FIQ), Univ. Nacional del Litoral, 3000 Santa Fe, Argentina oscar_quiroga@ieee.org*

^{2}

*IRII, Univ. Politècnica de Catalunya - CSIC, 08034 Barcelona, Spain*

riera@iri.upc.es

riera@iri.upc.es

^{3}

*Dep. de Matem. Aplic., Univ. Polit. de Catalunya, 08034, Barcelona, Spain*

carles@mat.upc.es

carles@mat.upc.es

* Abstract* ¾

**This paper presents the identification of a hydroelectric power turbine dynamics. Knowledge of power plant behaviour is fundamental to obtain reliable and efficient operation of power systems. Starting from models already proposed, some modifications are suggested in order to adjust real plant response, recorded from different conditions and situations, to model behaviour.**

* Keywords* ¾

**Identification. Hydroelectric Power Plant. Models.**

**I. INTRODUCTION**

The AGC (Automatic Generation Control) performance of power systems is strongly influenced by the dynamic characteristics of its power plants. It is therefore of certain importance to have accurate models of the plants that contribute significantly to the AGC. There are in the literature some well-established structures for those models, but it is necessary, in each particular case, to identify the model parameters. In this paper the identification of a hydroelectric power turbine dynamic model has been performed for the Susqueda power plant, which is a hydroelectric plant belonging to the Endesa Group in Spain. Models taken as an initial step are described, although some modifications have been carried out in order to adjust the response of the real plant to model behaviour. The identification approach makes use a group of models with complex hydraulic dynamics, since all of them are nonlinear, and the models consider surge tank effects.

Dynamic behaviour of the power plant has been recorded in different conditions and situations, and consists of the gate opening and the measured electric power, which may be considered as the measured mechanical power generated by the turbine.

The identification has used registers that correspond to normal work conditions and they have been chosen for having the most complete frequency spectrum, hence they guarantee appropriate identification results.

This paper is organised as follows: Section II describes the physical characteristics of the power station. Section III presents the dynamic equations considering a general nonlinear model with surge tank effects, and proposes some adjustments of the model. Section IV shows the adjustment of the equations. Section V presents the results of the simulation using alternative models. Section VI describes a comparative study of the behaviour of the models. Finally, Section VII summarises the conclusions of this paper.

**II. CHARACTERISTICS OF THE HYDROELECTRIC POWER PLANT**

The Susqueda power station is situated next to Susquedas reservoir in the province of Girona (Spain), which is supplied by the river Ter. The total installed power is 86 MW with three units (2 x 37 MW + 1 x 12 MW) and with an annual production of 180 GWh.

This work deals with one of the 37 MW units of the power station. Figure 1 presents a diagram of the main elements of the Susqueda power plant No 2 and shows the heads and flows that intervene in a hydroelectric plant model that considers surge tank effects.

Fig. 1. Graphic of the hydroelectric power plant.

The hydroelectric power plant No 2 of Susqueda has the following characteristics (Tables 1 and 2):

Table 1: Plant Characteristics.

Table 2: Characteristics of the conduits.

**III. GENERAL NONLINEAR EQUATIONS**

The initial approach employed in this work is to consider the most general dynamic equations and tune the resulting model to match the reality. The general dynamic equations are taken from different models from (IEEE Working Group, 1992; Kundur, 1994; Quiroga and Riera, 1999; Quiroga, 2000). The general equations are:

Equation of continuity

. (1)

Dynamics of the tunnel

, (2) |

. (3)

Dynamics of the surge tank

. (4) |

Dynamics of the penstock

, (5)

(6)

. (7)

Mechanical power

. (8)

The variables and parameters and their meaning are described in Table 3.

Table 3: List of variables and parameters.

**A. Formulas**

The parameters for the power plant are deduced from values of Tables 1 and 2 using the following formulas:

Elastic time in the penstock:

.

Water starting time in the penstock:

.

Hydraulic surge impedance of the penstock:

.

Water starting time in the tunnel:

.

Storage constant of surge tank:

.

Relationship between the normalised flow and the normalised water velocity in the tunnel or penstock:

.

Table 3 also shows the values of the main parameters of the plant of Susqueda.

**B. Block Diagrams**

A block diagram for the general nonlinear model with surge tank effects is shown in Fig. 2.

Fig. 2: Functional diagram for the general nonlinear model with surge tank effects.

To make the model equations useful, a key point is converting Eqn. (6) to a more treatable function. Figure 3 depicts the blocks diagram employed to calculate the exact hyperbolic tangent function.

Fig. 3: Blocks diagram to calculate the expression *z _{p}*×

*tanh*(

*T*×

_{ep}*s*).

The hyperbolic tangent can also be approximated by the alternative diagrams shown in Fig. 4.

Fig. 4: Blocks diagrams used to calculate the approximations n=0, n=1 and n=2 of the hyperbolic tangent.

**IV. ADJUSTMENT OF THE EQUATIONS**

In order to match model behaviour to real plant response, two kind of adjustments must be performed. The first is needed to adjust the static gain of the model to the Susqueda plant. The second one is related with the pressure oscillations due to the surge tank effects.

**A. Output power adjustment**

To adjust the output power () to real values it is necessary to multiply the value given by Eqn. (8), by a nonlinear function that represents the efficiency of the turbine. This function depends on the gate opening and its shape is similar to the efficiency curve of a Francis hydraulic turbine.

. (9)

Table 4 gives the values of deduced from steady state values of experimental registers considering that *P _{mech}* =

*P*. This last hypothesis holds because the efficiency of the generator is very high and can be considered almost constant for the whole range of operation. The nonlinear function is not considered in the general models given by (IEEE Working Group, 1992; Quiroga and Riera, 1999).

_{elec} Table 4: Values of the function for different gate positions deduced from experimental tests.

Moreover, the nonlinear function is calculated by means of the method of the least squares. The 1^{st}, 2^{nd}, 3^{rd} and 5^{th} order polynomials can be easily found using for example the function *polyfit* from the MATLAB mathematical package. This function finds the coefficients of a polynomial of degree m that fits the values given by the first and third column of Table 4 in the least square sense. Table 5 depicts these polynomials.

Table 5: Polynomials approximating .

Fig. 5: Plot of the nonlinear function calculated by the quadratic and the cubic polynomials.

Figure 5 shows the quadratic and the cubic polynomials that represent versus the steady state values taken from the third column of Table 4.

**B. Surge tank natural period**

The period (*T*) of the pressure waves due to the surge tank is a function of three physical parameters of the hydroelectric plant: the cross section area of the surge tank (*A _{s}*), the cross section area of the tunnel (

*A*) and the length of the tunnel (

_{c}*L*).

_{c}. (10) |

In order to adjust the oscillation period of the model with the period of the real plant it is necessary to slightly adjust the model parameters. Since the surge tank in Susqueda is cylindrical and very large, there exists the possibility that the diameter of the surge tank could vary a small percentage along its 100m of height.

Therefore, its diameter f * _{s}* may be used to adjust the surge tank natural period (

*T*). The modified value of the diameter is only five per cent less than the original diameter. The physical parameters f

*or*

_{c}*L*

_{c}(diameter and length of the tunnel) may also be used to adjust

*T*; however, the modified values for these parameters are greater than the five per cent obtained for f

*. Apart from this, the surge tank has not got an orifice; for this reason the value of the loss coefficient*

_{s}*f*is equal to zero.

_{0} **V. IDENTIFICATION RESULTS**

Table 6 presents a group of four nonlinear models for the power plant. These models differ in the way the block "In-Out" is implemented.

Table 6: Description of the Models A.

Figure 6 shows the measured electric power labelled as *P _{mech}*(register), and the mechanical power calculated by means of the Models A. Figures 7 and 9 show the gate opening, the measured electric power

*P*(register) and the mechanical power

_{mech}*P*(Models A).

_{mech}

Fig. 6: Identification of Susqueda using the Models A.

In Figs. 6, 7, 8, 9, 10 and 11 *P _{mech}* (Models A) represent the same curves for Models A1, A2, A3 and A4 since the differences are indistinguishable, however, more precise details are given in Section VI.

Fig. 7: Identification of Susqueda using the Models A.

Fig. 8: Detail of Fig. 7.

Fig. 9: Identification of Susqueda using the Models A.

Figures 8 and 10 show the detail of Figs. 7 and 9, respectively.

Fig. 10: Detail of Fig. 9.

In Figs. 6, 7, 8, 9, 10 and 11 the degree m=5 of the polynomial is used.

**VI. COMPARATIVE STUDY**

This section proposes a comparative study of the responses of the models, although the four models give good identification results of the Susqueda hydropower plant. A quadratic error index is defined by means of the following expression:

,

where *Ax* represents the models A1, A2, A3 or A4.

The register chosen to calculate the quadratic error is depicted in Fig. 11, where the hydraulic system reaches the steady state value after a variation in the gate opening. Table 7 describes the quadratic error of the models for polynomials with degree 2 and 5, approximating the nonlinear function .

The comparative study between the real records and the response of the models shows that the difference in the value of the quadratic error is between 0.15 to 0.39 for the Models A, giving quite a good representation of the real plant. Models A2, A3 and A4 have a quadratic error value that is a sixty per cent greater than the value of the quadratic error of the model A1, while the differences among the models A2, A3 and A4 are around 0.1 per cent.

Fig. 11: Identification of Susqueda using the Models A.

Table 7: Quadratic error found using the models A, where the nonlinear function is approximated by two polynomials with degree 2 and 5.

Model A1 has the least quadratic error between 0.1530 (for m=5) and 0.1554 (for m=2) and thus the best approximation to the real power plant, although the models A2, A3 and A4 are also good approximations.

**VII. CONCLUSIONS**

This paper has presented an identification process of the hydroelectric power station of Susqueda using the models A. These models are based on models taken from (IEEE Working Group, 1992 and Quiroga and Riera, 1999; Quiroga, 2000), and have been adjusted in order to match the power plant behaviour.

The identification process follows two adjustment procedures. The first one tunes the static gain of the Susqueda power plant by taking into account a nonlinear function in the calculation of the mechanical power. The shape of this function is similar to the efficiency curve of a Francis turbine.

The second one adjusts the oscillation time period (*T*), caused by the surge tank, by a five per cent modification of the measured value of the surge tank diameter.

The comparative study shows on one hand that the models A1, A2 and A3 need more complex expressions to calculate the hyperbolic tangent. On the other hand, the model A4 is the easiest to simulate since it considers non-elastic water columns and the hyperbolic tangent function is represented by means of a simple derivative function. Furthermore, the model A4 may be expressed as a nonlinear system in the state space and thus may be used in the nonlinear controller design.

**ACKNOWLEDGMENT**

This work has been partially supported by the Spanish CICYT DPI2000-1509-C03-02.

**REFERENCES**

1. IEEE Working Group, "Hydraulic Turbine and Turbine Control Models for System Dynamic Studies", *IEEE Transactions on Power Systems*, **7**, No.1, 167-179 (1992). [ Links ]

2. Kundur, P., *Power System Stability and Control*, McGraw-Hill, New York (1994). [ Links ]

3. Quiroga, O. D. and J. Riera, "Control Models for Hydroelectric Power Groups ", *SAAEI99*, Madrid, Spain, 645-648 (1999). (In Spanish). [ Links ]

4. Quiroga, O. D., "Modelling and Nonlinear Control of Voltage Frequency of Hydroelectric Power Plants", *Thesis Dissertation, Universitat Politècnica de Catalunya*, Barcelona, Spain (2000). [ Links ]