Improving The Response Time of Load – Frequency Control of Gasturbine in a Nigerian Refinery

DOI : 10.17577/IJERTV7IS010075

Download Full-Text PDF Cite this Publication

Text Only Version

Improving The Response Time of Load – Frequency Control of Gasturbine in a Nigerian Refinery

C.A.Ugoh, M. Olubiwe, D.S Inaibo

Department of Electrical and Electronic Engineering, Federal University of Technology,

Owerri, Nigeria.

Abstract – Gas turbine as a fast response unit, suitable for improving transient response of the power system at any time when there is problem usually from 15-20 seconds in some cases has been lost to some unforeseen issues due to energy demand, lack of control measures, improper attention , relatively high constraints in increasing and decreasing the power output during the nominal operation. Refinery which is one of the outfit of Nigeria National Petroleum Corporation (NNPC) is one of the most important sources of revenue to the federal government of Nigeria . Hence efficient energy and stability of energy is very paramount. This paper has presented a control system for improving the response time of speed/load-frequency of a typical gas turbine in a refinery in Nigeria. The objective of the study is to reduce the response time of the Gas turbine when there is abnormal increase in load or transient in the system, to enhance stability. . In order to realize the objective, the dynamic equations of a single shaft heavy duty gas turbine (HDGT) were obtained for load- frequency control operation. The dynamic equations were transformed into their equivalent simulink model in Matlab/Simulink environment. A lead phase compensator was initially developed for the control loop and maintaining the optimized performance specification, a proportional, integral and derivative (PID) controller was developed in Matlab/Simulink environment. The governor speed droop of 4% was used for the simulation. An analogue controller was designed. The essence of the analogue controller is only to demonstrate the fact that the control loop developed in this context used a permanent droop of 4% which gave the speed governor gain of 25 (per unit MW/per unit speed). Hence the optimized compensator performance specifications were then maintain and used to develop a PID structure with the same performance but with improved tracking. The designed PID controller was improved by adding a low pass filter to the derivative component. This greatly improved the transient response time. The simulations results showed that the introduction of the nonlinear PID controller significantly improved the transient response time of the system from 20s to 10.5s. This was achieved by introducing a nonlinear PID controller as part of the control loop, and consequently improved the performance thereby reducing the response time greatly.

Key words: Improving, Response-time, controller, Gas-turbine, frequency

I INTRODUCTION

Gas turbines have become more common as a power generating engine due to their flexibility and low specific investment costs as witnessed over the last decades.

Development of gas turbines for power generation has gone simultaneously with the development of new gas turbines for use in aircraft engines and refineries. This increase in the proportion of gas turbine and combined-cycle units in power generation grids will affect the grid characteristics because of the gas turbine performance [1]. In a power grid system, load frequency control (LFC) plays an essential role. It provides better conditions for power exchange and supply in trading electricity. Delay in such systems could reduce system performance and even cause frequency or other parameters instability within the system [2]. The dynamic behavior of the power systems and also their effects in industrial loads depend on disturbances and in particular on changes in the operating point [3]. The goal of the LFC is to maintain zero steady state errors in a multi area interconnected power system or stand alone units. Furthermore, the power system should fulfill the proposed dispatch conditions [4].Accurate modeling of power system components is required in the study of power system stability. The abundance of natural gas, couple with the relatively short time for planning to commercial operation, has contributed to the increase in use of gas turbine owing to the removal of government control from power market environment. Refineries being one of the important revenue generating parastatals of the government need to have efficient power generation to sustain productivity and to reduce downtime.

Governor and turbine are relevant components of the generating units. The gas turbine engine is a complex assembly of different components such as compressors, turbines, combustion chambers, etc., designed on the basis of thermodynamic laws [5]. Gas turbines usually consist of an axial compressor, a combustion chamber and a turbine operating under Brayton cycle [6]. These three elements make-up the thermal block and are complemented by the air intake system, the exhaust system, auxiliaries and controls.

The gas turbine as a fast response component is designed for improving transient response of the power system, but this function has been lost to some extents due to excessive energy demand. This problem also arises as a result of the relatively high constraints in increasing and decreasing the power output during the nominal operation. This problem most often persists if not checked or controlled. This paper is a contribution to the systematic analysis of load-frequency

control concepts for single shaft heavy duty gas turbines in aNigerian refinery. It is concerned with the improvement of the transient response time of load-frequency control of a gas turbine using a nonlinear proportional, integral and

where a , b , and c are the valve positioner constants. The fuel system actuator transfer function is:

derivative (PID) controller integrating a low pass filter.

II. MODEL DYNAMICS OF GAS TURBINE

A typical gas turbine model has basically three control loops

F (s)

1

T fc s 1

(2)

as shown below. The control loop, which is mainly active when the gas turbine is operating at nominal conditions is the load-frequency control (LFC). The input to this control is the speed/load set-point. The control section is restricted by the minimum fuel limit. The speed control or LFC is also the most important aspect during stability studies. The dynamics of the parts of the gas turbine that are essential in analyzing the transient characteristics of a load-frequency control in gas turbine are considered in this work and are presented below.

where T fc , is the fuel system actuator time constant in seconds.

The block diagram of the fuel system is represented in Fig. 2.

Wmin

a

1

Tfcs 1

+

VCE

Temperatu re

Temperature Control

N

km

bs c

+

Valve

Positioner Fuel System

k f actuator

Fig. 2: Block diagram of the fuel system Dynamics

Maximum Value

Acceleration control

Maximum Value

B. Compressor-Turbine Dynamics

The burning of the fuel in the combustor is presented by the following function:

m

m

e sTCR

(3)

Speed load set point

Acceleration set point

where TCR is the combustion reaction time delay constant in seconds.

The transfer function of the hot computation gas expansion is expressed as follows:

1

Fuel System Dynamics

T (s)

TCD s 1

(4)

Total

where

TCD

is the compressor discharge volume time

Gas turbine Dynamics

Fig. 1: Representation of gas turbine model

contant in seconds.

The mechanical torque produced which drives the electric generator is presented by the following equation [7].

A. Fuel System Dynamics

The fuel system consists of the fuel valve and actuator. The

Tm A Bm f

C1 N

(5)

amount of fuel injected for a gas turbine is determined by a valve positioner controlled by the turbine speed controller.

where

m f

is the fuel nominal flow rate, A and B are the

The valve positioner transfer function can be represented by the following transfer function [7].

a

coefficient of output torque which could be obtained by applying the actual data in Table 1.and N is the per unit rotor speed. The value for C in the torque equation (5) varies

V (s)

bs c

(1)

between 0.5 and 0.67 for heavy duty gas turbine (HDGT) [8]. In this context, C is assigned a value of 0.5.

Mechani

1

+ 1

Gas

N

Fig. 3: Block diagram of compressor-turbine dynamics

III CASE STUDY OF SELECTED GE GAS TURBINE

power would be the same. The calculated values of the parameters using the equations obtained from [8] are: m f =

3.33 kg/s, A = -0.4978, and B = 1.5002.

B Configuration of the Control System

The block diagram of the control system for realizing the objective of this work is shown in Figure. 4. In this work, the entire system is represented as a single input single output (SISO) system. The controller, the fuel system and the compressor-turbine assembly are all represented on the

OF A REFINARY IN NIGERIA

A 38MW single shaft heavy duty gas turbine (HDGT) of model GE whose design specifications are given in Table 1 is considered in this context. It uses natural gas or distillate as its primary operating fuel. In this context, the computations performed did not take into account the pressure loss in the entrance of air filters and also at combustor.

Parameters

Symbol

Unit

Value

Heat Rate

HR

KJ/Kwh

11302

Electrical power

PGT

MW

33.4

Exhaust Temperature

TOE

0C

538

Pressure Ratio

PR

11.8

Fuel Flow

Kg/s

139.1

Compressor Discharge Temperature

TCD

0C

343

Table 1: Design Specifications of HDGT Selected

forward part.

PID

+

r

– y

Compressor- Turbine

Rotor

Hi- Lo

X

+

+

Fuel system

The design specifications in Table 1 are provided for nominal operating conditions. Table 2 provides typical operational data selected for computing the turbine and compressor efficiencies. The turbine and compressor efficiencies are inherently assumed to be constant in power, put near nominal operation [8].

Table 2: Typical Operating Data for Computing Turbine and Compressor Efficiencies

Parameters

Unit

Value

Out Put Power

MW

30.2

Turbine Inlet Temperature

0 C

520

Exhaust Gas Temperature

0 C

200

Ambient temperature

0 C

26

Fuel flow

Kg/ s

0.99

Exhaust Mass Flow

Kg/s

85.9

Lower Heating Value of Fuel (LHV)

KJ/Kg

43309

The values of the efficiencies of turbine and compressor obtained using the parameters of Tables 1 and 2 were used to calculated the m f , A and B of the HDGT considered in

this context. The detail analysis on how to calculate the efficiencies of turbine and compressor can be obtained in [8].

A The mechanical torque

In order to extract the parameters: m f , A and B, the HDGT

considered has been assumed to be operating at nominal speed. The mechanical output torque and the mechanical

Fig. 4: Configuration of the control system

IV CONTROLLER DESIGN

During normal operating conditions, the speed/load- frequency control is active. The input to the controller in this case is the speed error or deviation. An analogue compensator or a proportional, integral and derivative (PID) controller may be required. In this work, the load-frequency control has been isolated and studied separately and it is expected that gain of control should be the inverse of the droop. Hence, an analogue lead-phase compensator is first developed in this context and a PID structure of the controller is developed with the same performance criteria.

  1. Analogue Compensator

    In this work, the controller is designed to maintain a reference speed. A nonlinear control approach has been adopted in this context to select the compensator parameters. This approach uses the Matlab/Simulink optimization toolbox. The gains of the compensator are constrained such that the system performance attains specifications of the transient response conditions. The constrained performance conditions are shown in Table 3

    Table 3: Constrain performance conditions

    Performance condition

    Values

    Rise time (s)

    0.2

    Settling time (s)

    10

    % Overshoot

    5

    % Rise

    90

    The analogue compensator is given by:

    W * Xs W

    Table 6: Model parameters and numerical values

    G(S )

    Ys Z

    (6)

    The obtained optimized selected parameters of the controller obtained by constraining system to the specification are stated in Table 4.

    Table 4: Selected gain parameters of the analogue compensator

    Parameters

    Values

    W

    25.0

    X

    0.2

    Y

    1.0

    Z

    5.0

    Substituting the parameters in Table 4 into Equation (6) yields:

    G(S ) 5s 25

    s 5

    (7)

  2. The Designed PID Controller

    The designed nonlinear PID controller is represented in equation (8). The controller was implemented using the Simulink PID block and the gains were selected by constraining the response of the system to meet the designed specifications. The values of the gains of the PID designed in this context are stated in Table 5.

    Table 5: The gains of the designed PID controller

    Parameters

    Gain

    KP

    609.760

    Ki

    7.683×104

    Kd

    112.805

    N

    10

    Parameter

    Symbol

    Unit

    Values

    Governor speed gain=1/droop

    W

    p.u.

    25

    Governor lead time constant

    X

    sec.

    0.2

    Governor lag time constant

    Y

    sec.

    0.05

    Control mode (1=droop, 0= isochronous)

    Z

    1

    Upper limit of fue demand

    Fmax

    p.u.

    1.5

    Lower limit of fuel demand

    Fmin

    p.u.

    0

    Valve positioned

    a

    1

    Valve positioned

    b

    -0.1

    Valve positioned

    c

    1

    Minimum fuel flow

    Wmin

    0.23

    Fuel control time constant

    TFC

    sec.

    0.4

    Fuel system feedback

    KFB

    0

    Combustion reaction time delay constant

    TCR

    sec.

    0.01

    Compressor discharge volume time constant

    ECD

    sec.

    0.2

    Turbine inertia constant

    Tt

    sec.

    15.64

    Gas turbine torque bock parameter

    A

    -0.4978

    Gas turbine torque bock parameter

    B

    1.5002

    Gas turbine torque bock parameter

    C

    0.5

    V SIMULATION RESULTS

    The simulation of the control system was performed in Matlab/Simulink environment. The parameters used in the simulations are presented in Table 6. The simulations were performed to check the performance of the system based on the transient response of the rotor speed (or frequency) of a single-shaft heavy duty gas turbine (HDGT). This is done as at the time the controller was not in the control loop and when it was integrated into the loop for better comparison..

    U (s) K

    I DNs s s N

    (8)

    The simulation results during the time the controller was not in the control loop and when it was, were performed in the droop mode. This control mode is the operational mode by

    where K k p E(s) 10, I ki E(s) 1.26 105 ,

    which the turbine is allowed to react to load variation by changing its speed. The droop of the gas turbine considered

    D kd E(s) 1.85 , and N

    coefficient. This yield:

    10 , which is the filter

    in this context is 4%. This means that for every change of 1% rated speed, the load of the machine is expected to change by 25% of its rated load. In order to give an obvious

    U (s) 10

    1.26 105

    s

    18.5s

    s 10

    (9)

    demonstration of the effectiveness of the designed nonlinear proportional integral and derivative (PID) controller with a low pass filter, simulations were performed for the designed

    The low pass filter in the PID controller will take care of disturbance that may affect the set point signal of the system and consequently eliminate cycling.

  3. Model Parameters and Numerical Values

The values of the parameters used for the simulation are given in Table 6.

analogue compensator such that the constrained conditions were then maintained and used to design the nonlinear PID controller. The nonlinear PID controller was implemented in two phases. In the first case the system was implemented with no low pass filter, and in the second case with a low pass filter whose coefficient, N =10 was integrated.

A Simulation Graphs

Graph of Rotor Speed against Time

1.4

1.4

Rotor Speed (p.u)

1.2

Graph of Rotor Speed against Time

step response of rotor speed

Rotor Speed (p.u)

1.2

1

0.8

0.6

0.4

0.2

step response of rotor speed

1

0.8

0.6

0.4

0.2

0

0 5 10 15 20 25 30 35 40 45 50

Time (s)

Fig. 8: Step response for speed control system (with PID controller)

0

0 5 10 15 20 25 30 35 40 45 50

Time (s)

Fig. 5: Step response for speed control

Graph of Rotor Speed against Time

1.4

system (with PID integrating a filter)

B Discussion

The governor speed droop of 4% is used for the simulation, as shown in speed-load control block. The transient response time of the speed/load-frequency control has been presented in Fig. 5, 6, 7and 8. The operation of the speed control is

1.4

Rotor Speed (p.u)

1.2

1

0.8

0.6

0.4

0.2

1.2

Rotor Speed (p.u)

1

0.8

0.6

0.4

0.2

0

step response of rotor speed

0 5 10 15 20 25 30 35 40 45 50

Time (s)

Fig. 6: Step response for speed control system (with no controller)

Graph of Rotor Speed against Time

based on the speed error formed when the HDGT rotor speed was compared with the referenced speed of one per- unit (1p.u.). The speed error or deviation was fed to the controller which in turn produced the control signal to the fuel system. Figures. 6 and 7 are the simulation results for improved rotor speed. The main objective of this work is to design a proportional, integral and derivative (PID) controller that will improved the transient response of a HDGT in a refinery in Nigeria. The essence of the analogue controller is only to demonstrate the fact that the control loop developed in this context used a permanent droop of 4% which gave the speed governor gain of 25 (per unit MW/per unit speed). It is just a design approach. Hence the optimized compensator performance specifications were then maintain and used to develop a PID structure with the same performance but with improved tracking. In Figure. 8 the designed PID controller was improved by adding a low pass filter to the derivative component. This greatly improved the transient response time of Figure.7. The simulations results show that the introduction of the nonlinear PID controller significantly improved the transient response time of the system from 20s to 10.5s.

step response of rotor speed

0

0 5 10 15 20 25 30 35 40 45 50

Time (s)

Fig. 7: Step response for speed control system (with compensator)

The essence of the low pass filter is to eliminate any disturbance that may enter the system through the set point signal. In all cases, the load/frequency control loop has a settling time of 10.5 seconds and improved overshoot performance. The result shows that the performance of the optimized control system has been significantly improved compared to the loop without the controller. This signifies improved transient response.

Finally, with the nature of the transient response time shown in Fig.5, it is obvious that the system will sluggishly rise above the required rated load or full load (set point = 1p.u.). This, will to a large extent, affect the stability of the system at nominal operating condition. In order to take care of this problem and improve the performance of the system, a nonlinear PID controller was introduced as part of the control loop, and consequently improved the performance.

  1. CONCLUSION

    This paper has presented a control system for improving the transient response time of speed/load-frequency of a gas turbin in a refinery in Nigeria. In order to realize the objective, the dynamic equations of a single shaft heavy duty gas turbine (HDGT) were obtained for load-frequency control operation. The dynamic equations were transformed into their equivalent simulink model in Matlab/Simulink environment. A lead phase compensator was initially developed for the control loop and maintaining the optimized performance specification, a proportional, integral and derivative (PID) controller was developed in Matlab/Simulink environment. A simple procedure in line with articles that have been published in literature was used for parameter estimation of Rowens model for a practical HDGT. The obtained parameters for a 33MW HDGT were used for simulation tests and the obtained results are presented. The gains of the designed nonlinear PID controller was realized using Matlab/Simulink optimization toolbox rather than the Ziegler-Nichols. The nonlinear PID controller with a low pass filter was able to improve the transient response performance of the load-frequency control loop..

  2. REFERENCES

  1. Saeed Bahrami , Ali Ghaffari and Marcus Thern (2013). Improving the Transient Performance of the Gas Turbine by Steam Injection during Frequency Dips, Energies 2013, 6, 5283-5296; doi:10.3390/en6105283

  2. Al-Hamouz Z. M, Abdel-Magid Y. L (1993). Variable-structure load frequency controllers for multi area power systems. Int. J. Electr.

    Power Energy Syst., 15(5): 23-29

  3. Chang C. S, Weihui Fu (1997). Area load frequency control using fuzzygain scheduling of PI controllers. Electric Power Syst. Res., 42: 145-152.

  4. Behnam Jamshidzadep and Shahram Jamali (2011). An analysis on the application of PID controller in the frequency contribution of gas turbine power plant, International Journal of the Physical Sciences Vol. 6(19), pp. 4519-4527, 16 September, 2011 DOI: 10.5897/IJPS11.371, ISSN 1992 – 1950 ©2011 Academic Journals

  5. Cohen H., Rogers G.F.C., Saravanamuttoo H.I.H. (1996). Gas Turbine Theory, 4th Edition, Longman, London.

  6. Harman R. T. (1981). Gas Turbine Engineering: Applications, Cycles and Characteristics, Macmillan, 1981, pp. 270.

  7. Rowen W. I (1983). Simplified mathematical representations of heavy-duty gas turbines, ASME J. Eng. Power, Vol. 105, pp. 865869

  8. BankTavakoli M. R., Vahidi B., and Gawlik W. (2009). An Educational Guide to Extract the Parameters of Heavy Duty Gas Turbines Model in Dynamic Studies Based on Operational Data, IEEE Trans. Power Syst., Vol. 24, No.3. August 2009.

Leave a Reply