Next Article in Journal
Large Eddy Simulation of Vertical Axis Wind Turbine Wakes
Previous Article in Journal
Impact of Hot-Water Extraction on Acetone-Water Oxygen Delignification of Paulownia Spp. and Lignin Recovery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling and Internal Fuzzy Model Power Control of a Francis Water Turbine

1
Department of Electrical Engineering, Sava Hydroelectric Power Plants Ljubljana d.o.o., Gorenjska Street 46, 1215 Medvode, Slovenia
2
Department of System, Control and Cybernetics, Faculty of Electrical Engineering, Tržaška Street 25, SI-1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Energies 2014, 7(2), 874-889; https://doi.org/10.3390/en7020874
Submission received: 17 November 2013 / Revised: 6 February 2014 / Accepted: 8 February 2014 / Published: 19 February 2014

Abstract

: This paper presents dynamic modelling of a Francis turbine with a surge tank and the control of a hydro power plant (HPP). Non-linear and linear models include technical parameters and show high similarity to measurement data. Turbine power control with an internal model control (IMC) is proposed, based on a turbine fuzzy model. Considering appropriate control responses in the entire area of turbine power, the model parameters of the process are determined from a fuzzy model, which are further included in the internal model controller. The results are compared to a proportional-integral (PI) controller tuned with an integral absolute error (IAE) objective function, and show an improved response of internal model control.

1. Introduction

In HPPs, a digital turbine governor is an indispensable part of the control system. The stability of frequency, active power control, water flow control, turbine start-up procedure and emergency shut-down implemented with algorithms are major functions of digital turbine governor. Connected to Supervisory Control and Data Acquisition (SCADA), the unit controller, the excitation system and the digital voltage controller, it enables the operator to change the operating states of generator, i.e., the rotational speed, active power, reactive power and voltage.

With reference to achieving suitable control results and exploring the dynamic responses of a hydro power plant unit, it is necessary to obtain a mathematical hydraulic model of the hydro turbine and other parts of water system. Different types of hydraulic non-linear models were proposed in [1]: a non-linear model assuming a non-elastic water column, a model including travelling waves, and a model including a surge tank effect. In [27], the authors applied a non-linear or linearised model with a non-elastic water column to study the transient behaviour of a Francis turbine. Chen et al. provided a nonlinear dynamical model with surge tank presented with state space equations, analysed with bifurcation diagrams, phase orbits and spectrograms [8], however a comparison to the real system was neglected. In [9] authors developed interesting non-linear and linearised hydro turbine models with transient stability analysis in the free and open source Power System Analysis Toolbox (PSAT) software, although a conventional controller without a sophisticated tuning method was used.

Many recommendations of modelling, design and testing control systems for hydraulic turbines are described in international standards [10,11] and a few modern methods in [12,13], nevertheless in most cases a proportional-integral-derivative (PID) or proportional-integral (PI) controller is used for speed, active power, water level, flow and wicket gate opening control. However, the parameter determining process is usually taken on site acceptance tests; some recommendation limits are proposed in [11].

In [14] a particle swarm method was applied as a strategy in an improved gravitational search algorithm for identification of an optimal water turbine model and controller parameters. Fang et al. introduced the particle swarm optimization method of controller tuning in water turbine frequency control [15]. In [16] an internal model control (IMC) tuned PID controller in comparison with Ziegler-Nichols and singular frequency based tuning was studied for frequency control of hydropower system with a water hammer effect. The design of a PID controller based on sensitivity margin specifications was introduced in [17], where the author proposed the controller with respect to sensitivity peak specification. As a result in these cases, authors obtained optimal tuned PID parameters for water turbine control that could not necessarily result in optimal transient behaviour due to simplifications of the models and nonlinearity of real processes at different operating points.

In this paper, the focus is on mathematical modelling and the fuzzy power control approach of a Francis turbine of Aggregate 1 at hydro power plant (HPP) Moste. In the second section, the paper includes a comparison of non-linear and linearised models, with and without the surge tank of the Francis Turbine 1 in HPP Moste. In the third section, a fuzzy model of the Francis turbine has been applied based on a first order linearised model. Section 4 presents a fuzzy IMC controller with tuning parameter T to ensure a suitable dynamic response in working points and middle sections. At the end of Section 4, the comparison to the PI controller has been applied, tuned by the integral absolute error (IAE) objective function.

2. Non-Linear and Linearised Model of Turbine Dynamics

2.1. Description of the Hydro Power Plant

Slovenian hydro power production consists of three HPP chains placed on the Sava, Soca and Drava Rivers. As the first HPP on the Sava River chain with a unique accumulative operation type, HPP Moste plays an important role in the Slovenian hydro power plant generation portfolio.

HPP Moste is part of the HPP Moste and HPP Zavrsnica water system depicted in Figure 1. Besides 2,940,000 m3 of useful volume of water storage, a −6.25 m water level deviation is allowed. After the water intake object, there is a surge tank connected by a tunnel. The pressure penstock connects the surge tank with a powerhouse, comprised of two Francis turbines (HPP Moste) and one Francis turbine (HPP Zavrsnica) with a joint tail-water conduit.

2.2. Turbine Model with Water Conduit

Turbine dynamics with the penstock are derived from second Newton's law [1,18] described with Equation (1):

d q ¯ d t = ( h ¯ 0 h ¯ h ¯ f h ¯ s p ) g A l
where q̄ is turbine flow rate in m3/s, with the penstock area section A in m2 and length l in m. In Equation (1), h̄0 is the static head of the water column; h̄ is the head at the turbine admission; h̄f is defined as friction head loss; h̄sp is the tail water head loss in m units and g is the gravity acceleration.

With reference to water inertia time [10]:

T w = Q r g H r i = 1 n L i A i

Equation (1) becomes expressed in per unit:

d q d t = 1 h h f h s p T w

In Equation (2), the water inertia time Tw includes the characteristics of all water passages with area sections Ai; corresponding lengths Li with rated discharge Qr; rated head Hr and gravity g. The change of water flow in the conduit is expressed in Equation (3), where h is the head at turbine; hf is the friction head loss in the penstock and hsp is the tail water head loss per unit. Water flow through the turbine [5] is given by:

q = G h
where G is the wicket gate opening per unit of the Francis turbine and h is the head at the turbine admission per unit.

Turbine power Pm in p.u. is expressed by:

P m = A t h ( q q n l ) D G Δ w
where Δw is the speed deviation; D is the damping factor and qnl is defined as no-load turbine flow. Proportionality factor At as described in [1] is given by Equation (6):
A t = P m , r P g , r h r ( q r q n l )
where Pm,r is the maximum turbine rated power; Pg,r is the generator rated power; hr is the rated head and qr is the rated turbine flow in p.u.

With purpose of analysing dynamic responses, linearisation in the surrounding area of operating point R was made with a Taylor series procedure (considering first two elements), e.g., linearisation of the non-linear head at turbine admission h:

h = q 2 G 2 = h ¯ + ( h q ) R ( q q ¯ ) + ( h G ) R ( G )
with calculation of partial derivation and substitution of:
Δ h = h h ¯ ; Δ q = q q ¯ ; Δ G = G

We obtain linearised Equation (9):

Δ h = 2 h ¯ q ¯ Δ q 2 h ¯ Δ G

The linearisation process for all relations between Equations (3) and (5) was made with the real parameters of the Francis turbine 1 on HPP Moste. According to the synchronised generator and the fixed network frequency Δw = 0 assumption, only the first part of Equation (5) will be concerned in further analysis. The transfer function between turbine power Pm and wicket gate opening G is proposed in Equation (10):

Δ P m Δ G = b 0 s + b 1 a 0 s + a 1
with calculated numerator coefficients:
b 0 = h ¯ ( q ¯ q n l ) A t T w
b 1 = h ¯ h ¯ A t h ¯ A t ( q ¯ q n l ) ( 2 q ¯ Q r 2 ( f p 1 + f p 4 ) + 2 q ¯ a 0 Q r 2 + a 1 Q r H b )
and denominator coefficients in Equations (13) and (14):
a 0 = 1 2 q ¯ h ¯ T w
a 1 = 1 + 1 2 q ¯ h ¯ ( 2 q ¯ Q r 2 ( f p 1 + f p 4 ) + 2 q ¯ a 0 Q r 2 + a 1 Q r H b )

In connection to negative numerator coefficient b0 in Equation (11), the non-minimal phase type of dynamic system is observed. In Figure 2, the matching of the nonlinear and calculated linear models is presented with simulation results at operating point Pm = 0.3422, G = 0.4416, q = 0.431, h = 0.9524, hf = 0.0076 and hsp = 0.015, all p.u. with real ramp reference of the wicket gate opening, in comparison to real system measurement.

2.3. Surge Tank Modelling

In high-head HPPs, the surge tank is often included in the water system for providing the limits of water pressure fluctuation and water speed decreases in water tunnels and penstocks. There are some other important functions of surge tanks that influence the water system [18,19]: hydraulic separation of water tunnel and pressure penstock, decreasing of water hammer effect evaluation, improved results of turbine control, and acceleration of water level stabilisation.

With reference to model improvement, the surge tank with Equations (15) and (16) were included in the previously described Francis turbine model with conduit:

Q s = Q w + Q T
where Qs is the water flow in the tunnel; Qw is the surge tank water flow (positive or negative) and QT is the turbine water flow. Water level of surge tank [1,18] is described in Equation (16):
d z d t = Q s Q T C S
with Cs expressed in Equation (17), where Aw is surge tank section area:
C S = A w h base q base

The surge tank implemented on HPP Moste has two extended horizontal chambers with enlarged upper and bottom surface areas, A1_kom and A3_kom. Therefore, three different areas Aw (A1_kom, A2_kom and A3_kom) were added to the model, depending on absolute surge tank water level calculated with relation:

A w = f ( z abs ) ; z abs = H abs H b , max ( z H b , max )
where Habs is the nominal level of HPP water storage; Hb,max is the maximum available head and zabs is the absolute level of surge tank.

The friction loss of head hs and hfr with friction factors fp2 and fp3 in the water tunnel are presented in Equation (19):

h s = Q s | Q s | f p 2 h f r = f p 3 Q s 2
and the loss of head due to narrowed intake of the surge tank:
h w = Q w | Q w | f 0
where f0 is the intake loss coefficient. The net head of surge tank is expressed in p.u. with Equation (21):
H vod = z h w
which is also presented in differential Equation (22), describing flow in the water tunnel:
d Q s d t = H 0 h s h f r H vod T w 1
and Equation (3) is modified in Equation (23):
d Q T d t = H vod h s h r h s p T w 2
where H0 is static head of water column; Tw1 is the water inertia time of the water passages before the surge tank intake and Tw2 is the water inertia time between the surge tank and the Francis turbine. A dynamic system with differential Equations (16), (22) and (23) was linearised with the same Taylor series procedure as in the model without the surge tank. The third order transfer function of the HPP Moste turbine model with a conduit and surge tank was derived for three different turbine power working points: Pm = 0.3408 p.u., Pm = 0.6657 p.u. and Pm = 0.9684 p.u., with the transfer function shown in Equation (24):
Δ P m Δ G = s 3 d 0 + s 2 d 1 + s d 2 + d 3 s 3 c 0 + s 2 c 1 + s c 2 + c 3

Coefficients of the numerator, denominator and all other parameters in simulation are related to the real system (Francis turbine 1 HPP Moste with penstock and surge tank) and are described in Appendixes A.1. and A.2.

With the real ramp reference (slower opening rate) obtained from the HPP measurement (from 5 to 6 MW), the transient response is practically without undershot (Figure 3). With 1260 s of the simulation time, approximately ±1% of fluctuating power is presented after the reference was changed (ΔG = 0.0903).

Surge tank level Hvod and water flow Qw were also observed (Figure 4) during the transient response in the nonlinear model (ramp reference). The maximum difference of the water level in the surge tank was 0.0096 p.u. (0.68 m).

3. Fuzzy Model

In connection with the high similarity of transient responses (models with/without surge tank) depicted in Figure 3, a first order model (without surge tank) is used for fuzzy modelling. A transfer function derived from Equation (10) with time constants T1, T2 and gain K is presented in Equation (25):

G s ( s ) = K ( 1 T 1 s ) 1 + T 2 ( s )

Due to the nonlinearity of the process, a fuzzy model was realised to determine the parameters in the middle sections of the working points. Membership functions, describing degrees of membership for input parameter Pm and output parameters K, T1 and T2 are presented in Figure 5.

Membership functions with triangular and trapezoidal shapes were determined based on the knowledge acquired from operating experiences and intuitiveness. Fuzzy rules added to the fuzzy model are presented with rule list interpretation:

IF turbine power Pm = small, THEN:
K = medium
T1 = small
T2 = small
IF turbine power Pm = medium, THEN:
K = high
T1 = medium
T2 = medium
IF turbine power Pm = high, THEN:
K = small
T1 = high
T2 = high

Determination of output fuzzy set was made with an aggregation method (maximum) and with a defuzzification process, where the centroid method (also called centre of area, centre of gravity) was used to obtain crisp values of parameters [20]. Turbine model parameters in working points and middle sections obtained with fuzzy model are listed in Table 1.

4. Fuzzy IMC Control and Comparison to PI Control

In connection to IMC calculation, the fuzzy model of the Francis turbine proposed in Section 3 is used. With IMC, a model is parallel-added to the process (Figure 6) and the difference between the output of the process and model is used for the feedback signal.

In an ideal case, when disturbance d = 0 and assuming the model is identical to the process, the feedback is inactive. Conversion to classic control form is made with Equation (26):

G R ( s ) = G IMC ( S ) 1 G IMC ( s ) G M ( s )

In the first step, the minimal phase and all-pass part of the process transfer function are expressed with Equation (27):

G s ( s ) = K ( 1 T 1 s ) 1 + T 2 ( s ) = K ( 1 + T 1 ( s ) ) 1 + T 2 ( s ) 1 T 1 ( s ) 1 + T 1 ( s )

The ideal IMC controller is an inversed minimal phase part that results in mirror-transformed zeros over the ordinate axis:

G IMC ( opt , ideal ) ( s ) = 1 + T 2 ( s ) K ( 1 + T 1 s )

Due to easier realisation and increased robustness, a PT1 (first order) filter with tuning parameter T is added to the ideal IMC controller:

G IMC ( opt , real ) ( s ) = 1 + T 2 ( s ) K ( 1 + T 1 s ) ( 1 + T s )

Considering Equation (26) and Figure 6, conversion to a classic controller GR(s) is expressed with Equation (30):

G R ( s ) = 1 + T 2 s K T T 1 s + ( K T + 2 K T 1 ) 1 s

Referring to different parameters in working points from the fuzzy model, the IMC controller is presented with a matrix of transfer functions MR(s), with different tuning parameter T:

M R ( s ) = [ G R 1 ( s ) ; P m = 0.34081 ; T = 0.3 G R 2 ( s ) ; P m = 0.6657 ; T = 0.2 G R 3 ( s ) ; P m = 0.9684 ; T = 0.08 ]

The results of IMC with a nonlinear model and surge tank are presented in Figure 7, with 0.05 p.u. step response of turbine power and output of controller signal (wicket gate opening) in the first (G = 0.4416, Pm = 0.3408) and second working point (G = 0.6367, Pm = 0.6657). In the third working point (G = 0.9, Pm = 0.9684), a −0.05 p.u. step response was realised. The results shows a small undershoot (first and second working point) and a small overshoot at the third working point.

The fuzzy IMC controller algorithm was additionally formed with a fuzzy system to determine parameter T depending on turbine power Pm; while other parameters K, T1 and T2 are observed from a turbine fuzzy model. Results after the defuzzification process (centroid method) of tuning parameter T are presented in Table 2.

If a hydraulic model of turbine and water passages are known, a PID or PI controller could be used for turbine power control, tuned with IAE expressed with Equation (32) or integral time absolute error (ITAE) objective function (33) with a simulation procedure [11].

| P P r | d t = min
t | P P r | d t = min

In Equations (32) and (33), Pr is the reference power and P is the actual power output. In connection to the IMC result validation, a PI controller with the same model was proposed. The tuning process was made with the IAE objective function, expressed with Equation (32), and tested in the same working points as the fuzzy IMC controller. A similar fuzzy system (with input—Pm, output Kp and KI) was realised with calculated PI parameters, listed in Table 2.

Figure 8 presents turbine wicket gate opening G, turbine flow Q, head at turbine admission h; friction head loss hf; and tail water hsp for both controllers in Pm = 0.484 middle section. The influence of the PI controller proportional part fast response results in a 6% higher variation of head h at turbine admission, which consequently causes a higher under-control effect of turbine power, compared to fuzzy IMC control.

The fuzzy IMC controller in middle section Pm = 0.484 provides 0.200 p.u. undershoot Mp in comparison to 0.915 p.u. with PI controller (Figure 9). Comparison results for the four middle sections of working points are presented in Table 2 (bolded values).

The smallest under-control power effect with 0.200 p.u. undershoot is obtained with IMC in Pm = 0.484 (Table 3); however, the highest under-control effect 1.088 p.u. is observed at Pm = 0.7859 with PI controller. Turbine head drop Δh has the smallest value 0.020 p.u. with the IMC controller in Pm = 0.881; and the highest value 0.086 p.u. with the PI controller in Pm = 0.484 middle section.

In comparison to the PI controller, fuzzy IMC presents better results, with lower undershoots Mp in all four middle sections and with approximately the same settle time on the step reference of turbine power Pm.

Fuzzy IMC also provides smaller variations of head Δh at turbine admission, which ensures fewer oscillations of the HPP water system.

5. Conclusions

In this paper, the dynamic modelling of a Francis turbine 1 (HPP Moste) with and without a surge tank with linearisation and IMC power control was proposed. With the actual HPP technical data (water tunnels and conduits), the presented models show a good match to the real data measurement.

With linearisation, a first-order transfer function with non-minimal phase was used for fuzzy modelling. Concerning appropriately defined input and output membership functions, a fuzzy rule list and the defuzzification process (centroid method), the crisp parameters of the Francis turbine model were calculated.

Based on the fuzzy model, a fuzzy IMC controller was proposed with calculation of tuning parameter T, depending on turbine power. For IMC validation, the PI controller with IAE objective function and fuzzy inference system (similar to IMC) was proposed. A comparison of the fuzzy IMC and PI control with maximum undershoot Mp and head variations Δh criterion in the middle sections of working points shows that IMC has better results due to the smaller undershoot on step response, smaller head variations Δh and easier tuning process (one tuning parameter T) in case of acquired model.

The applied methodology with the extension of working points with different static heads can be used as algorithm in digital turbine governor in order to ensure better results in water turbine power control.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix

A1. Numerator and Denominator Coefficients of Transfer Function Presented in Section 2.3. (24)

d 0 = h ( Q ¯ A T q n l A T ) T w 1 T w 2 H b C s
d 1 = h ( Q ¯ A T q n l A T ) ( C s T w 1 ( 2 Q ¯ Q r 2 f p 1 + 2 Q ¯ a 0 Q r 2 + a 1 Q r ) + 2 T w 2 f p 3 Q s ¯ H b C s ) + h ¯ h ¯ A t H b C s T w 1
d 2 = h ¯ ( Q ¯ A T q n l A T ) ( H b ( T w 1 + T w 2 ) + 2 f p 3 Q s ¯ C s ( 2 Q ¯ Q r 2 f p 1 + 2 Q ¯ a 0 Q r 2 + a 1 Q r ) + 2 h ¯ h ¯ A t f p 3 H b C s Q s ¯
d 3 = h ( Q ¯ A T q n l A T ) ( 2 Q ¯ Q r 2 f p 1 + 2 Q ¯ a 0 Q r 2 + a 1 Q r + 2 f p 3 Q s ¯ H b ) + h ¯ h ¯ A t H b
c 0 = 1 2 Q ¯ h ¯ T w 1 T w 2 H b C s
c 1 = 1 2 Q ¯ h ¯ ( C s T w 1 ( 2 Q ¯ Q r 2 f p 1 + 2 Q ¯ a 0 Q r 2 + a 1 Q r ) + 2 T w 2 f p 3 Q s ¯ H b C s ) H b C s T w 1
c 2 = 1 2 Q ¯ h ¯ ( H b ( T w 1 + T w 2 ) + 2 f p 3 Q s ¯ C s ( 2 Q ¯ Q r 2 f p 1 + 2 Q ¯ a 0 Q r 2 + a 1 Q r ) + 2 f p 3 Q s ¯ H b C s
c 3 = 1 2 Q ¯ h ¯ ( 2 Q ¯ Q r 2 f p 1 + 2 Q ¯ a 0 Q r 2 + a 1 Q r + 2 f p 3 Q s ¯ H b ) + H b

A2. Model Parameters

Table A1. Model parameters.
Table A1. Model parameters.
ParameterSymbolValue
Generator rated powerPg, r9 MVA
Maximum turbine rated powerPm, r7.5 MW
Rated turbine powerPr6.692 MW
Maximum turbine flowQmax14 m3/s
Rated turbine flowQr13 m3/s
No-load flowqnl0.19 p.u.
Rated turbine head p.u.hr0.824 p.u.
Static head p.u.h00.9740 p.u.
Static headhb68.62 m
Gravity constantg9.80665 m/s2
Water inertia timeTw12.2361 s
Water inertia timeTw20.7267 s
Water inertia timeTw32.9628 s
Surge tank intake lossf00.0005
Friction loss penstockfp10.0136
Friction loss water tunnelfp20.0136
Additional friction loss of the water tunnelfp30.01
Additional friction loss of the penstockfp40.003
Surge tank parameterCs225.5150 s
Turbine dampingD0.5 p.u./p.u.
Surge tank upper chamber areaA1_kom316.6 m2
Surge tank main chamber areaA2_kom44.2 m2
Surge tank bottom chamber areaA3_kom276.4 m2

References

  1. De Mello, F.P.; Koessler, R.J.; Agee, J.; Anderson, P.M.; Doudna, J.H.; Fish, J.H., III; Hamm, P.A.L.; Kundur, P.; Lee, D.C.; Rogers, G.J.; et al. Hydraulic turbine and turbine control models for system dynamic studies. IEEE Trans. Power Syst. 1992, 7, 167–179. [Google Scholar]
  2. De Jaeger, E.; Janssens, N.; Malfliet, B.; van de Meulebroeke, F. Hydro turbine model for system dynamic studies. IEEE Trans. Power Syst. 1994, 9, 1709–1715. [Google Scholar]
  3. Ng, T.B.; Walker, G.J.; Sargison, J.E. Modelling of Transient Behaviour in a Francis Turbine Power Plant. Proceedings of the 15th Australasian Fluid Mechanics Conference, Sydney, Australia, 13–17 December 2004.
  4. Fraile-Ardanuy, J.; Wilhelmi, J.R.; Fraile-Mora, J.; Pérez, J.I.; Sarasua, I. A Dynamic Model of Adjustable Speed Hydro Plants. Proceedings of the 9th Congreso Hispano Luso de Ingeniería Eléctrica, Marbella, Spain, 30 June–2 July 2005.
  5. Hannett, L.N.; Feltes, J.W.; Fardanesh, B. Field tests to validate hydro turbine-governor model structure and parameters. IEEE Trans. Power Syst. 1994, 9, 1744–1751. [Google Scholar]
  6. Němec, Z.; Němec, V. Water turbine power control feasible conception. Eng. Mech. 2009, 16, 403–411. [Google Scholar]
  7. Win, T.N. Hydraulic Governor Control System for Francis Turbine. Proceedings of at the GMSARN International Conference on Sustainable Development: Issues and Prospects for the GMS, Kunming, China, 12–14 November 2008.
  8. Chen, D.; Ding, C.; Ma, X.; Yuan, P.; Ba, D. Nonlinear dynamical analysis of hydro-turbine governing system with a surge tank. Appl. Math. Model. 2013, 37, 7611–7623. [Google Scholar]
  9. Li, W.; Vanfretti, L.; Chompoobutrgool, Y. Development and implementation of hydro turbine and governor models in a free and open source software package. Simul. Model. Pract. Theory 2012, 24, 84–102. [Google Scholar]
  10. International Electrotechnical Committee (IEC). IEC 60308: Hydraulic Turbines—Testing of Control Systems, 2nd ed.; IEC: Geneva, Switzerland, 2005. [Google Scholar]
  11. International Electrotechnical Committee (IEC). IEC 61362: Guide to Specification of Hydraulic Turbine Control Systems, 1st ed.; IEC: Geneva, Switzerland, 1998. [Google Scholar]
  12. Mahmoud, M.; Dutton, K.; Denman, M. Design and simulation of a nonlinear fuzzy controller for a hydropower plant. Electr. Power Syst. Res. 2005, 73, 87–99. [Google Scholar]
  13. Tan, W. Unified tuning of PID load frequency controller for power systems via IMC. IEEE Trans. Power Syst. 2010, 25, 341–350. [Google Scholar]
  14. Chen, Z.; Yuan, X.; Tian, H.; Ji, B. Improved gravitational search algorithm for parameter identification of water turbine regulation system. Energy Convers. Manag. 2014, 78, 306–315. [Google Scholar]
  15. Fang, H.; Chen, L.; Shen, Z. Application of an improved PSO algorithm to optimal tuning of PID gains for water turbine governor. Energy Convers. Manag. 2011, 52, 1763–1770. [Google Scholar]
  16. Naik, K.A.; Srikanth, P.; Negi, P. IMC tuned PID governor controller for hydro power plant with water hammer effect. Procedia Technol. 2012, 4, 845–853. [Google Scholar]
  17. Hušek, P. PID controller design for hydraulic turbine based on sensitivity margin specifications. Int. J. Electr. Power Syst. 2014, 55, 460–466. [Google Scholar]
  18. Giesecke, J.; Mosonyi, E. Wasserkraftanlagen: Planung, Bau und Betrieb, 2nd ed.; Springer: Berlin, Germany, 1998; pp. 293–305. [Google Scholar]
  19. Mosonyi, E. Water Power Development: High-Head Power Plants, 3rd ed.; Akadémiai Kiadó: Budapest, Hungary, 1991; Volume 2/B, pp. 975–985. [Google Scholar]
  20. Ross, T.J. Fuzzy Logic with Engineering Applications, 2nd ed.; John Wiley & Sons: Chichester, UK, 2004; p. 101. [Google Scholar]
Figure 1. Hydro power plant (HPP) Moste and HPP Zavrsnica water system.
Figure 1. Hydro power plant (HPP) Moste and HPP Zavrsnica water system.
Energies 07 00874f1 1024
Figure 2. Turbine power of the non-linear model (red solid line), linearised model (blue dashed line) and the real system measurement (green dashed line).
Figure 2. Turbine power of the non-linear model (red solid line), linearised model (blue dashed line) and the real system measurement (green dashed line).
Energies 07 00874f2 1024
Figure 3. Dynamic response on the ramp reference of realised models at turbine power working point 0.6657 p.u.
Figure 3. Dynamic response on the ramp reference of realised models at turbine power working point 0.6657 p.u.
Energies 07 00874f3 1024
Figure 4. Simulated water level and water flow of the surge tank.
Figure 4. Simulated water level and water flow of the surge tank.
Energies 07 00874f4 1024
Figure 5. Membership functions of Francis turbine model.
Figure 5. Membership functions of Francis turbine model.
Energies 07 00874f5 1024
Figure 6. Internal model control.
Figure 6. Internal model control.
Energies 07 00874f6 1024
Figure 7. IMC control of nonlinear turbine model with conduit and surge tank for 3 different working points.
Figure 7. IMC control of nonlinear turbine model with conduit and surge tank for 3 different working points.
Energies 07 00874f7 1024
Figure 8. Wicket gate opening, turbine flow, head at turbine admission, friction head loss, and tail water in p.u. with fuzzy internal model control (IMC) and proportional-integral (PI) control in the middle section at turbine power 0.484 p.u.
Figure 8. Wicket gate opening, turbine flow, head at turbine admission, friction head loss, and tail water in p.u. with fuzzy internal model control (IMC) and proportional-integral (PI) control in the middle section at turbine power 0.484 p.u.
Energies 07 00874f8 1024
Figure 9. Step response of the turbine power with fuzzy IMC and PI controller in the middle sections of working points.
Figure 9. Step response of the turbine power with fuzzy IMC and PI controller in the middle sections of working points.
Energies 07 00874f9 1024
Table 1. Francis turbine model parameters in working points and middle sections, calculated with fuzzy model.
Table 1. Francis turbine model parameters in working points and middle sections, calculated with fuzzy model.
Turbine powerKT1T2
Pm = 0.340811.35830.76100.6614
Pm = 0.4841.35650.89220.7150
Pm = 0.5421.39351.28630.9139
Pm = 0.66571.44701.39180.9506
Pm = 0.78591.35371.66941.0515
Pm = 0.8811.29292.09281.2593
Pm = 0.96841.26712.30731.3348
Table 2. Tuning parameter T for the internal model control (IMC) and proportional-integral (PI) parameters (Kp and KI).
Table 2. Tuning parameter T for the internal model control (IMC) and proportional-integral (PI) parameters (Kp and KI).
Turbine powerT (IMC)KP (PI)KI (PI)
Pm = 0.340810.30000.29830.3154
Pm = 0.4840.27860.29120.3016
Pm = 0.5420.21780.27090.2594
Pm = 0.66570.20000.27640.2437
Pm = 0.78590.16710.25200.2318
Pm = 0.8810.10640.23290.1932
Pm = 0.96840.08000.22250.1759
Table 3. Turbine maximum undershoot Mp and head drop Δh at turbine admission with fuzzy IMC and PI control.
Table 3. Turbine maximum undershoot Mp and head drop Δh at turbine admission with fuzzy IMC and PI control.
Comparison IMC/PI controlPm = 0.484Pm = 0.542Pm = 0.7859Pm = 0.881
Mp IMC (p.u.)0.2000.2020.3400.410
Mp PI (p.u.)0.9150.9901.0881.076
Δh IMC (p.u.)0.0270.0220.0220.020
Δh PI (p.u.)0.0860.0830.0610.050

Share and Cite

MDPI and ACS Style

Nagode, K.; Škrjanc, I. Modelling and Internal Fuzzy Model Power Control of a Francis Water Turbine. Energies 2014, 7, 874-889. https://doi.org/10.3390/en7020874

AMA Style

Nagode K, Škrjanc I. Modelling and Internal Fuzzy Model Power Control of a Francis Water Turbine. Energies. 2014; 7(2):874-889. https://doi.org/10.3390/en7020874

Chicago/Turabian Style

Nagode, Klemen, and Igor Škrjanc. 2014. "Modelling and Internal Fuzzy Model Power Control of a Francis Water Turbine" Energies 7, no. 2: 874-889. https://doi.org/10.3390/en7020874

APA Style

Nagode, K., & Škrjanc, I. (2014). Modelling and Internal Fuzzy Model Power Control of a Francis Water Turbine. Energies, 7(2), 874-889. https://doi.org/10.3390/en7020874

Article Metrics

Back to TopTop