Open Access
Issue
Mechanics & Industry
Volume 22, 2021
Article Number 8
Number of page(s) 11
DOI https://doi.org/10.1051/meca/2021005
Published online 08 March 2021

© A.A. Golneshan and H. Nemati, Hosted by EDP Sciences 2021

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

In spite of the considerable improvement in electricity generation by new and green sources of energy, it is still undeniable that electrical power mainly depends on fossil fuels. On the other hand, it seems that the adverse effect of using fossil fuels on nature and climate besides the depletion of fossil fuel reserves are avoidable. So, it is required to use fossil fuels more efficiently.

Methods used in optimizing power plants may be divided into two distinct categories; methods that may be used during the design and methods which are applied during operation. Several references may be found that proposed the optimization of power plants [110]. Most of these works are based on exergy destruction minimization. The Exergy of a system at a specific thermodynamic state is the maximum available work that a system provides when reaches the equilibrium state relative to the surroundings.

It is known that [8,9] boiler systems as well as the turbines are the main sources of exergy destruction. Kaushik et al. [7] concluded that the largest exergy loss occurs during the combustion process. Hou et al. [11] found that in a steam power plant, most of the exergy destructs in steam turbines and in the boilers. In fact, irreversibilities in the combustion process and heat leak to ambient causes a large amount of exergy destruction in the boiler. The main contribution of the combustion process in exergy destruction was approved by Ahmadi and Toghraie [6]. They also showed that in the second stage, the steam turbine is the most considerable source of destroying exergy.

More and less, the same results may be found as researches conclusions. As it can be observed, these valuable results mostly concern with the first category aforementioned above and must be considered in design steps. Because considering these results lead to major changes in power plants and are usually avoided.

The present study, unlike those studies, is concerned with the second category. Maintenance policy and repair strategy have been always big questions during overhaul. Which equipment does need special care? For which one routine inspection or repairmen is enough? These questions are critical questions that can affect the performance of a power plant. However, those aforementioned studies may not answer these questions simply.

In such conditions, comparing similar plants may be a powerful technique to improve the performance of under operation plants. To show the ability of this method, the Fars power station is selected as the pilot. The station has six similar gas turbine power cycles. The performance of the main components was compared with each other for one day. Data were gathered every 5 min in one day and thoroughly analyzed. Results of the comparison showed the strength and or weakness of each component. Later, it was deduced that which one was working more efficiently and which one needed maintenance. Based on this comparison, some advice has been suggested for further improvements. Owing to this fact that the comparison took place over one-day length; the impact of the ambient temperature was studied in reality.

This study is classified into five parts. After the introduction, in the second part, plant details are specified. Required equations are derived and presented in part 3. Since the evaporative coolers were at the service on that afternoon, the common assumption of dry air and using low heating value are not applicable. Part 4 is devoted to results and discussions. In this part, the comparisons of similar components from different viewpoints are presented. As a final point, a brief conclusion is presented in part 5. Moreover, complementary information is summarized in the appendices.

2 Plant details

The power station is comprised of 123.4 MW six gas. Combustion chambers are feed with natural gas. Ambient air is fed into the compressor through a filter and a small pressure drop occurs in the filtration process. At the ambient temperature around 40 °C, the evaporative cooler installed after the air filter comes into service. Spraying water droplets into the air stream decrease compressor inlet temperature and increase its humidity.

In the compressor, air pressure increases and is delivered to the combustion chamber. In the combustion chamber, the temperature increases in an approximately isobaric process. Outlet gas generates work in the turbine and then goes to the stack.

3 Analysis

Thermodynamic analysis for the power plant is presented in the following. Unless otherwise specified, the variation of kinetic energy as well as the potential energy is neglected. Also, flow is assumed steady-state. For a steady-state process, the mass balance and energy balance for each component can be written as: o m o = i m i (1) i m h + Q = o m h + W (2)

The required data like the pressure and temperature of all components of six gas plants were gathered on July, 27th. Data gathering started from 09:00 till 16:00 for every five minutes. A datasheet of measured data for gas unit number 1 is shown in Appendix A.

Figure 1 shows a schematic diagram of the power plant. In all components, air or flue gas is considered as an ideal gas.

Using the evaporative coolers, the air moisture content is considerably high. So, fuel low heating value beside the dry air assumption is not applicable, and more precise equations are needed in combustion simulation.

thumbnail Fig. 1

A schematic diagram of the power plant.

3.1 Evaporative cooler

Enthalpy of moisture and air mixture is: h 1 = h a i r 1 + h w a t e r 1 18 ω 1 . (3)in which ω1 is the specific humidity of inlet air to the compressor. The mass-weighted average of the specific heat capacity of the mixture is assumed to be only a function of temperature: c p 1 m i x = c p a i r ( T 1 ) + c p w a t e r ( T 1 ) ω 1 1 + ω 1 (4) c p 2 m i x = c p a i r ( T 12 ) + c p w a t e r ( T 12 ) ω 1 1 + ω 1 (5)in which T 12 is the average temperature between T1 and T2. R m i x = R a i r + R w a t e r ω 1 1 + ω 1 (6)

For times that evaporative cooler is shut down, ω = ω0 and otherwise ω = ω1: ω 1 = ω 0 [ h w a t e r ( P w , T 0 ) h f , w a t e r ( T 0 ) ] c p a i r ( T 01 ) [ T 1 T 0 ] [ h g , w a t e r ( T 1 ) h f , w a t e r ( T 0 ) ] (7) hg,water and hf,water are the saturated steam and saturated liquid specific enthalpy, respectively. In Eqution (8), the water partial pressure, Pw, is calculated by [12]: P w = ϕ 0 P s a t , w a t e r ( T 0 ) (8) ϕ 0 = ω 0 P 0 ( 0.622 + ω 0 ) P s a t , w a t e r ( T 0 ) (9)

In which Psat,water(T0) is the water saturated pressure at temperature T0. The inlet water vapor mass flow rate to the gas unit is: m m o i s t , 0 = m d r y õ a i r ω 0 (10)

Also, the vapor mass flow rate inlet to the compressor is: m m o i s t , 1 = m d r y õ a i r ω 1 (11)

Therefore, the total air mass flow rate is: m a i r = m d r y õ a i r + m m o i s t , 1 (12)

3.2 Compressor

The compressor isentropic pressure ratio can be calculated as: ( P 2 s P 1 ) = [ T 2 T 1 ] c p 2 m i x R m i x (13)in which P2s is compressor isentropic outlet pressure.

3.3 Combustion chamber

Based on the fuel analysis, used natural gas components are (Tab. 1):

Analysis of exhausted flue gas from the stack shows the complete combustion. Therefore, the flue gas components are (Tab. 2):

The combustion in the combustion chamber can be simulated thermodynamically as: i n i [ h f + { h ( T 3 ) h ( 298.15 K ) } ] i j n j [ h f + { h ( T 3 ) h ( 298.15 K ) } ] j = Q (14)

In the above equation, hf is the formation enthalpy and i and j, indicates to the ith component of inlet and the jth component of the outlet of the combustion chamber respectively. Q is the released heat for each mole of fuel.

Knowing xi as mass fraction, Ji as mole fraction and Wi as molar mass of ith component and Wfuel as fuel molar mass, the thermodynamic properties of flue gas is: c p 1 f l u = { ( x C o 2 c pC o 2 + x h 2 o c p h 2 o + ( x N 2 + 3 . 76 C W N 2 ) c p N 2 + 7.63 C ω 1 W H 2 o c p h 2 o + ( C - J O 2 ) W O 2 c p O 2 ) } ( 1 + ω 1 ) ( 137.33 C ) + W f u e l (15)

All above used properties are at the outlet temperature of the turbine. C is the available air mole (including excess air) for each mole of the fuel. Gas constant (R) of flue gas mixture is calculated as the weighted function of individual components: R f l u = { R C o2 x C o 2 + R H 2 o x H 2 o + ( x N 2 + 3 . 76 C W N 2 ) R N 2 + 7 . 63 C ω 1 W H 2 o R H 2 o + ( C - J O 2 ) W o 2 R O 2 } ( 1+ ω 1 ) ( 137 .33 C ) +100 W F u e l (16)in the same way, cp 2flu is evaluated at the outlet temperature of the combustion chamber.

Table 1

Fuel content.

Table 2

The flue gas contents.

3.4 Turbine

The isentropic temperature at the turbine outlet is evaluated as: T S , 4 = ( P 1 P 2 ) R f l u c p 1 f l u [ T 3 T 34 ] c p 2 f l u c p 1 f l u T 34 (17)in which T 34 is the average temperature between T4 and T3. Consequently, the isentropic efficiency of the turbine is: η T = c p 2 f l u [ T 3 - T 34 ] + c p 1 f l u [ T 34 - T 4 ] c p 2 f l u [ T 3 - T 34 ] + c p 1 f l u [ T 34 - T S , 4 ] (18)

It is assumed that P2 ≈ P3 and P1 ≈ P4. Moreover: h4=JCo2hCo2+JH2ohH2O+(C-JO2)hO2+(JN2+3.76C)hN2+Cω1137.3318hH2O137.33C(19)

in the above equation, all the thermodynamic properties are at T4 temperature. In the same method, h3 may is evaluated combustion chamber temperature, T3. Temperature functionality of all mentioned thermodynamic variables is presented in Appendix B.

The mass flow rate of the dry air is calculated using known fuel component and mass flow rate of the fuel, m F u e l ,: m d r y õ a i r = m F u e l W F u e l C × 28 .851 × 4 .76 (20)So, the electricity output power is: W E = m d r y õ a i r [ ( h 3 h 4 ) ( h 2 h 1 ) ] η G e n (21)

ηGen is the efficiency of the power generator.

Unfortunately, combustion chamber exit temperature was not accessible. So, it was calculated indirectly by a try and error algorithm. For this reason, for a first guess of T3, C and after that, m d r y õ a i r were calculated at and based on these values, W E was calculated employing equtaion (21). This method was continued using a better guess for T3, till the measured value of W E and the calculated value was the same.

4 Results and discussions

4.1 The evaporative coolers

Ambient and inlet compressors temperatures are shown in Figure 2. According to this figure, the range of the ambient temperature is between around 30 °C at 9:00 and around 40 °C at noon. Interestingly, before starting an evaporative cooler, the air inlet temperature to the compressor of unit 1 is about 3– 4 °C warmer than other units. This shows that the location of the compressor unit 1 causes that air enters the compressor with higher temperature and this can be an object of the unit improvement. At 12:00, evaporative coolers start working. These evaporative coolers cool the air entrance temperature and reduce compressors consumption power. Worthwhile to note that, all evaporative coolers do not work in the same level of performance. Because there are a considerable differences among coolers outlet temperatures. The evaporative cooler of the 1st unit reduces the temperature from 40 to 22 °C, while cooler number 6 or number 3, can decrease temperature from 38 to 30 °C. So, this difference in outlet temperature is subject to question.

The differences in cooler performances shall be investigated through their structure since the air inlet temperature to all the coolers is approximately the same. For this reason, the relative humilities of coolers outlet air are presented in Figure 3. It can be observed that, for cooler number 1, the outlet relative humidity is around 93% while for cooler numbers 2, 4, and 5, this value something between 60% and 70%. For cooler 3 and 6 the condition is even worse and the relative humidity is around 45%.

thumbnail Fig. 2

Comparison of compressors inlet temperature and ambient temperature through one day.

thumbnail Fig. 3

Relative humidity of inlet air to compressors in one day.

4.2 Compressors

Variations of compression ratio for all six compressors against outlet to inlet temperature ratio are presented in Figure 4. Aggregated points are related to times when evaporative coolers are in service and the compressors air inlet temperatures are independent of ambient temperature and are relatively constant. In the same manner, scatter points are related to the morning working condition. At first glance, it clear that decreasing the inlet temperature increases the compressor discharge pressure. Moreover, comparing the pressure ratio for the same temperature ratio in both conditions (with and without evaporative cooler), the pressure ratio when the evaporative cooler is in the service is much more.

However, this variation is not the same for all compressors. Although the inlet conditions of compressors are not the same, some parts of these deviations are due to dissimilarities in operation despite this fact that the compressors are the same.

To find out, which compressor works better and also to separate the contribution of inlet condition; κc may be a good parameter for judgment. κc is introduced as the ratio of isentropic to real compressor outlet pressure. Consequently, it shows how far the compressor outlet pressure is from isentropic outlet pressure: κ c = P 2 s P 2 (22)

It is clear that in an ideal compressor, the κc is equal to unity. Larger κc shows that the compressor is far from its ideal isentropic operation. Variations of κc against temperature ratio are shown in Figure 5 for both conditions, when evaporative coolers are in the service and when they are out of service. Figure 5 shows that, when the evaporative coolers are off, compressors numbers 2, 4, and 6 work efficiently in comparison to others.

Also, compressor number 2 performance is superior to the others in both working conditions. So, it is worth wide to compare other compressors with compressor number 2 during the overhaul.

All in all, it can be claimed that if all the coolers worked like cooler compressor number 2 and all the compressors like compressor number 1, then the compressor discharge pressure in all units would reach more than 9.1 bars; which means a considerable increase in output power. This judgment can be supported by Figure 6. This figure shows generated power outputs against the compressor discharge pressure. Based on Figure 6, an increase of 5% in compressor number 6 discharge pressure results in 7% increase in generated power, or on the other hand, 7% increase in discharge pressure of the compressor number 2, increases output power 15%. These examples show how the output power is sensitive to compressor discharge pressure.

thumbnail Fig. 4

Compression ratio against temperature ratio.

thumbnail Fig. 5

Variations of κc against temperature ratio.

thumbnail Fig. 6

Plants power outputs against compressor discharge pressure.

4.3 Turbine

Turbine is entropic efficiency is defined by equation (18). To have good judgment about the performances of gas turbines, the efficiency of each turbine has been shown in Figure 7 against its compressor temperature ratio, T 3 / T 4 , in both working conditions, when evaporative coolers are working and when they are out of service. This figure shows how far the efficiency of a turbine is from the ideal point. First of all, it shows that in both working conditions, by increasing the temperature ratio, T 3 / T 4 , the turbine efficiency gets closer to its ideal efficiency. Moreover, a meaningful difference may be observed between Turbine 1 with the others. So, it can be recommended that all turbines get compared deeply with turbine 1 during overhaul.

thumbnail Fig. 7

Turbines isentropic efficiencies against temperature ratio.

4.4 Power plants efficiencies

The efficiency of each power plant was calculated. The efficiency is a good index to compare units, overall. The results of this comparison are presented in Figure 8. Since evaporative coolers are not in the service all day, units efficiencies are classified into three groups:

  • When the evaporative coolers are in service (in the afternoon),

  • when the evaporative coolers are off (in the morning)

  • and one for one day (the combination of two above).

Based on Figure 8, except for unit 3, the performances of all other units are approximately the same. However, the overall efficiency of unit 3 is slightly lower. Moreover, since the evaporative cooler in unit one has a higher performance when the evaporative coolers are in service, this unit has a higher efficiency than the others.

thumbnail Fig. 8

Gas units efficiencies.

5 Conclusion

In this study, the required parameters for monitoring a power plant was introduced and calculated for six similar gas turbine. These parameters were compared in one day. As a result, it was cleared that which equipment works well and which one needs more care.

Some concluded results in these comparisons are:

  • All evaporative coolers do not show the same performance. Cooler number 1 works better than the others while cooler number 3 or cooler number 6 is not as well as cooler number 1.

  • Reducing the inlet temperature of the compressor increases the compressor discharge pressure.

  • The compressor number 1 shall be verified thoroughly during overhaul.

  • Performance of the compressor number 2 is better than others. So, it is worthwhile to compare other compressors with compressor number 2 during the overhaul.

  • For a compressor in the same temperature ratio, the pressure ratio when the evaporative cooler is in the service is much more in comparison to the condition in which the evaporative cooler is off.

  • It must be tried to reform all compressors to reach the discharge pressure of compressor number 2 and all coolers work like cooler number 1. Under this condition, compressor discharge pressure would reach 9.1 bars and consequently, a considerable gain in power.

  • Since turbine 1 acts better than the other turbines, it is recommended that all turbines get compared deeply with turbine 6 during overhaul.

  • Unit number 4 has lower efficiency among all units and unit numbers 6 and 1 have the best efficiencies.

  • Due to the high performance of evaporative cooler number 1, the overall efficiency of unit 1, when the cooler is in service, is the highest.

  • Using the evaporative coolers increase the plant efficacies significantly.

The above findings show the benefits of comparing similar cycles. These results simply show the ability of this method in cycle performance improvement. Valuable repair or maintenance guidelines can be in hand if one casts periodically this method.

Nomenclature

cp(T): Specific heat capacity at temperature T (kJ/kg K)

h : Specific enthalpy (kJ/kg)

J : Mole fraction

m : Mass flow rate (kg/s)

Q : Rate of heat transfer (kW)

R : Gas constant (kJ/kmol K)

T : Temperature (K or °C)

W : Molar mass (kg/ kmol)

W : Power (kW)

x : Mass fraction

Greek symbols

ϕ : Relative humidity

ω : Specific humidity

ηT : Turbine isentropic efficiency

ηTh : Unit thermodynamic efficiency

Subscripts

C : Compressor

Gen : Generator

i : Inlet

mix : Mixture

moist : Moisture

o : Outlet

S : Isentropic

sat : Saturated

T : Turbine

Appendix

Appendix A

Sample of measuring data for gas unit number 1

Appendix B Temperature functionality of thermodynamic variables

(T in Kelvin) θ=T/100

for air [10,14] cpair=(1.04841-3.8371θ103+9.4537θ21055.4931θ3107+7.9218θ41010 (kJ/kgK)(B1)

for other gases [12] cpN2=[39.06512.79θ1.5+1072.7θ2820.4θ3]/28.013 (kJ/kgK)(B2) cpO2=[37.432 + 0.020102θ1.5178.57 θ1.5+236.88 θ2]/31.999 (kJ/kgK)(B3) cpCO2=[3.7357 + 30.529θ0.54.1034θ+0.024198θ2]/44.01 (kJ/kgK)(B4) cpH2O=[143.05183.54θ0.25+82.751 θ0.53.6989 θ]/18.015 (kJ/kgK)(B5) cpCH4=[672.87 + 439.74 θ0.2524.875 θ0.75+323.88θ0.5]/16.04 (kJ/kgK)(B6) cpC2H6=[6.895 + 17.26θ0.6402θ2+ 0.00728θ3]/30.07 (kJ/kgK)(B7) cpC3H8=[4.042+30.46θ1.571θ2+0.03171θ3]/44.09 (kJ/kgK)(B8) cpC4H10=[3.954+37.12θ1.833θ2+0.03498θ3]/58.124 (kJ/kgK)(B9)

Regression based on thermodynamic table [12], hair=4.666416+0.96831θ+ 3.961×105θ2+3.2648×108θ3 (kJ/kg)(B10) hN2={40848385.458+13666978θ1.0024978/(4711.4676+θ1.0024978)for θ48670+ 2923.5942θ-12.044984θ2+2.136434θ3for θ<10exp(8.07894354.8921249θ+1.0352864ln(θ))for θ10 (kJ/kmol)(B11) hO2={8683.2986+2915.5857θ11.791232θ2 + 3.6922661θ3for θ<48349+2672θ+43θ2for θ47939.7229+2451.347θ+80.132035θ21.8838384θ3for θ<108879.333+2882.7507θ+31.329828θ20.28252785θ3for θ10 (kJ/kmol)(B12) hCO2={9363.977+2939.3076θ83.011986θ2+50.577468θ3forθ<49525+2646θ+184θ2forθ49584.9286+2538.8571θ+240.08333θ26.41667θ3forθ<10exp(8.65466615.9395707θ+1.0233514ln(θ))forθ10 (kJ/kmol)(B13) hH2O={9903.9842+3250.3251θ+42.806253θ26.2183832θ3forθ<49598+3094θ+42θ2forθ49494.7937+3083.6892θ+31.208332θ2+1.9537038θ30.04166667θ4forθ<1010205.43+2943.4614θ+73.011089θ20.63909098θ3forθ10 (kJ/kmol)(B14)

References

  1. R. Bontempo, M. Manna, Efficiency optimisation of advanced gas turbine recuperative-cycles, Proceedings of the Institution of Mechanical Engineers, Part A: Journal of Power and Energy (2019) [Google Scholar]
  2. A.A. Golneshan, H. Nemati, Energy and exergy comparison of three similar combined power cycles, Mechanics & Industry 20 , 504 (2019) [EDP Sciences] [Google Scholar]
  3. S.M. Bapat, G.D. Gokak, Exergetic evaluation and optimization of Combined Heat and Power (CHP) plant of 20.7 MW capacities under varying load conditions: a case study exergetic evaluation and optimization of Combined Heat and Power (CHP) plant of 20.7 MW capacities under varying load conditions: a case study, European Journal of Electrical Engineering 21 , 207–219 (2019) [Google Scholar]
  4. A.A. Golneshan, H. Nemati, Evaluation of six gas turbine evaporative cooling for Fars-Iran, Annales de Chimie, Science des Materiaux 42, 281 (2018) [Google Scholar]
  5. A.S. Karakurt, Performance analysis of a steam turbine power plant at part load conditions, Journal of Thermal Engineering 3 , 1121–1128 (2017) [Google Scholar]
  6. G.R. Ahmadi, D. Toghraie, Energy and exergy analysis of Montazeri Steam Power Plant in Iran, Renewable and Sustainable Energy Reviews 56 , 454–463 (2016) [Google Scholar]
  7. S. Kaushik, V.S. Reddy, S. Tyagi, Energy and exergy analyses of thermal power plants: a review, Renewable and Sustainable Energy Reviews 15 , 1857–1872 (2011) [Google Scholar]
  8. P. Regulagadda, I. Dincer, G.F. Naterer, Exergy analysis of a thermal power plant with measured boiler and turbine losses, Applied Thermal Engineering 30 , 970–976 (2010) [Google Scholar]
  9. I.H. Aljundi, Energy and exergy analysis of a steam power plant in Jordan, Applied Thermal Engineering 29 , 324–328 (2009) [Google Scholar]
  10. M. Ameri, P. Ahmadi, S. Khanmohammadi, Exergy analysis of a 420 MW combined cycle power plant, International Journal of Energy Research 32 , 175–183 (2008) [Google Scholar]
  11. D. Hou, S. Shao, Y. Zhang, S.L. Liu, Y. Chen, S.S. Zhang, Exergy analysis of a thermal power plant using a modeling approach, Clean Technologies and Environmental Policy 14 , 805–813 (2012) [Google Scholar]
  12. C. Borgnakke, R.E. Sonntag, Fundamentals of thermodynamics, Wiley Global Education, 2016 [Google Scholar]
  13. M.J. Moran, H.N. Shapiro, D.D. Boettner, M.B. Bailey, Fundamentals of Engineering Thermodynamics, Wiley, 2010 [Google Scholar]
  14. H. Kurt, Z. Recebli, E. Gedik, Performance analysis of open cycle gas turbines, International Journal of Energy Research 33 , 285–294 (2009) [Google Scholar]

Cite this article as: A.A. Golneshan, H. Nemati, Comparison of six gas turbine power cycle, a key to improve power plants, Mechanics & Industry 22, 8 (2021)

All Tables

Table 1

Fuel content.

Table 2

The flue gas contents.

Appendix A

Sample of measuring data for gas unit number 1

All Figures

thumbnail Fig. 1

A schematic diagram of the power plant.

In the text
thumbnail Fig. 2

Comparison of compressors inlet temperature and ambient temperature through one day.

In the text
thumbnail Fig. 3

Relative humidity of inlet air to compressors in one day.

In the text
thumbnail Fig. 4

Compression ratio against temperature ratio.

In the text
thumbnail Fig. 5

Variations of κc against temperature ratio.

In the text
thumbnail Fig. 6

Plants power outputs against compressor discharge pressure.

In the text
thumbnail Fig. 7

Turbines isentropic efficiencies against temperature ratio.

In the text
thumbnail Fig. 8

Gas units efficiencies.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.