Self-Unit Commitment of Combined-Cycle Units with Real Operational Constraints

Abstract

This paper highlights the importance of accurately modeling the operational constraints of Combined-Cycle Gas Turbines (CCGTs) within a unit-commitment framework. In practice, in Colombia, when given an initial dispatch by the Independent System Operator, CCGT plants are operated according to the results of heuristic simulation codes. Such heuristics often omit technical operating constraints, including hot, warm, or cold startup ramps; the minimum operation hours required for a gas turbine to start a steam turbine; the relationship between the dispatched number of steam and gas turbines; the load distribution among gas turbines; and supplementary fires. Most unit-commitment models in the literature represent standard technical constraints like startup, shutdown, up/down ramps, and in some cases, supplementary fires. However, they typically overlook other real-life CCGT operating constraints, which were considered in this work. These constraints are crucial in integrated energy systems to avoid equipment damage, which can potentially put CCGT plants out of service and ultimately lead to lower operating costs.

1. Introduction

1.1. Motivation

In modern power systems, power plant operators and Independent System Operators (ISOs) rely on mathematical models to simulate various operating conditions for specific scenarios, ensuring power system security and reliability. These models are crucial, especially for Combined-Cycle Gas Turbine (CCGT) plants, which are prevalent worldwide due to their high efficiency and flexibility in integrating renewable energy resources [1,2,3]. Accurate representation of CCGTs’ operational elements is essential for simulating the correct power output available to the ISO for meeting demand and avoiding equipment damage [4,5,6].

The Colombian power system, with its diverse CCGT plants, exemplifies this need. The largest plant, TEBSA, located in northern Colombia with a 5 × 2 configuration (5 gas turbines and 2 steam turbines), plays a vital role in system reliability and security. Its actual power output depends on the thermal states of the gas turbines and their interaction with the steam turbines, which require specific temperature and pressure conditions to operate effectively.

Currently, most Colombian CCGT plants, including TEBSA, lack a mathematical optimization model for detailed operation planning based on dispatch. Moreover, the Colombian ISO lacks a faithful technical representation of these plants. This gap often leads to operational inefficiencies and technical issues, such as equipment failures or deviations from unit-commitment programs. For instance, TEBSA’s operations often assume a hot startup condition or overlook the minimum number of gas units needed for steam turbine operation, leading to preventable equipment issues [7] (see Figure 1). Addressing these challenges requires a comprehensive optimization model that accurately represents the complex operating conditions of CCGTs.

Therefore, this paper proposes a Self-Unit-Commitment (SEUC) model, formulated as a Mixed-Integer Programming (MIP) problem, to address these challenges and to accurately represent the detailed operating conditions of a CCGT plant given a specific dispatch. The model introduces two novel aspects: firstly, it considers the minimum operation hours required for a gas turbine to start a steam turbine, a factor not previously included in existing Unit-Commitment (UC) models to the best of our knowledge; secondly, it accounts for the load distribution among gas turbines. A crucial feature of our model is its treatment of the CCGT plant, representing gas and steam turbines separately, rather than as a single entity. This granular representation enables the ISO to execute a more-precise unit-commitment plan for these plants, ensuring better accuracy in active and reactive power output. Such accuracy is vital for supporting grid constraints and adhering to the operational rules of the CCGT, thereby reducing program deviations (Although reactive power output is not directly considered, the model’s active power outputs for each CCGT unit can be utilized in power system analysis software by the ISO, where the capability curve of each unit at different voltage levels is modeled). Furthermore, this model can be integrated into a broader UC framework for dispatching CCGT plants within a general power system.

1.2. Literature Review

Modeling approaches for Combined-Cycle Gas Turbines (CCGTs) in the literature vary significantly. One method involves representing CCGT components individually, focusing on the separate functionality of gas and steam turbines, as seen in studies like [8,9,10,11,12,13]. Another common approach models the entire CCGT plant using configurations or modes, as detailed in [2,14,15,16,17,18,19,20,21,22]. Our work adopted an individual component representation, but included constraints typical of CCGT configurations, such as power output limits, startup/shutdown ramps, and the number of gas and steam units. This approach is similar to that in [12], but our model uniquely combines both component and mode models.

CCGTs are often modeled using either configurations or individual representations of gas and steam turbines. Works like [2,14,18] propose a Mixed-Integer Programming (MIP) model for Unit Commitment (UC), considering CCGT operation modes and addressing ramping and minimum uptime/downtime constraints. They also consider startup/shutdown costs associated with mode transitions, improving computational performance. However, these models typically overlook operational constraints like the minimum operation hours required for a gas turbine to start a steam turbine and load distribution requirements.

Refs. [8,9] propose componentwise formulations for CCGTs, incorporating them into UC problems. These models address the interdependence of gas and steam turbines, with the steam turbine’s power output dependent on the gas turbine’s performance. Ref. [8] focuses on ramping constraints on a minute-by-minute basis, aligning with the hourly dispatch requirements of the Colombian energy market. However, this model does not consider crucial constraints like minimum startup hours, differentiation between hot and cold startups for steam turbines, and the load distribution among gas turbines, which are essential for maintaining operational integrity and preventing equipment failures.

Our work presents a hybrid model that captures the benefits of both component- and configuration-based approaches. Unlike [23], which proposes a similar hybrid model, our approach is tailored to the Colombian energy market framework, considering the specific operational constraints and bidding policies for CCGT units. We incorporated supplementary fires, an element also modeled in [8,23], which enhance the steam turbine’s power output without increasing the gas turbine’s output. This was achieved using the Heat Recovery Steam Generator (HRSG), with the assumption that steam quality is maintained when the operational constraints are met.

1.3. Contributions

The main contributions of this paper are as follows:

• The introduction of the SEUC model for CCGT hourly dispatch, accommodating individual the minimum uptime/downtime constraints for each gas and steam unit, ramp rates, and operational rules for steam turbines, including the requisite number of gas turbines for optimal steam production.
• The incorporation of steam turbine startup processes in the model, taking into account the thermal state of each unit (hot or cold) and the minimum operation hours required for a gas turbine to initiate a steam turbine.
• The implementation of a load distribution constraint among gas turbines. This feature is crucial for maintaining uniform temperature conditions across the steam turbine rotors, thereby preventing thermal stress and ensuring consistent steam production for the steam turbines.

2. Mathematical Formulation

In this section, we present the novel mathematical formulation of the Self-Unit-Commitment (SEUC) model of a CCGT.

2.1. Objective Function

The self-unit commitment of the CCGT presented in this paper aims at minimizing the total operational cost of the CCGT given by the following components: (i) the cost of the CCGT in each hourly period; (ii) the cost of non-served energy, the energy related to the wasted steam, and the auxiliary consumption; (iii) the startup costs of the gas turbines; (iv) the costs related to the additional fire, as follows:

$$\begin{array}{c}{\displaystyle \sum _{t\in \ll }}{\displaystyle \sum _{c\in NC}}PCC·(g_{tc}+a{f}_{tc})+{\displaystyle \sum _{t\in \ll }}PBC·(g{r}_t+f_t)+\hfill \\ {\displaystyle \sum _{t\in \ll }}{\displaystyle \sum _{c\in NC}}CSC·y_{tc}+{\displaystyle \sum _{t\in \ll }}{\displaystyle \sum _{s\in NS}}CSS·y_{ts}+\hfill \\ {\displaystyle \sum _{t\in \ll }}{\displaystyle \sum _{c\in NC}}{\displaystyle \sum _{cc\ne c\in NC}}DSC·vhd{r}_{t,c,cc}+{\displaystyle \sum _{t\in \ll }}CCC·c{c}_t\hfill \end{array}$$

(1)

2.2. CCGT Operational Constraints

This section describes the constraints related to the energy balance and limits to the operational constraints for each unit and the CCGT.

2.2.1. Energy Balance Constraints

Constraint (2) outlines the load to be met, where $L_t$ represents the initial input load and $g{r}_t$ denotes the energy that the CCGT is unable to serve. The variable $e_t$ is a slack variable introduced to prevent infeasibilities during the startup and shutdown ramps when the CCGT output falls below the minimum output of the configuration. Lastly, $gcc{t}_t$ signifies the total output of the CCGT.

$$\begin{array}{c}\hfill gcc{t}_t+g{r}_t-e_t=L_t\forall t\end{array}$$

(2)

$$\begin{array}{c}\hfill 0\le e_t\le \overline{GCC}\forall t\end{array}$$

(3)

Equations (4)–(6) delineate the operational limits, specifically the minimum and maximum output of the gas turbines $g_{tc}$, steam turbines $g_{ts}$, and the CCGT $gc{c}_t$. The term $gc{c}_t$ specifically refers to the total output when the CCGT is operating in combined-cycle mode.

2.2.2. Gas Turbines’ Operational Limits

Equations (4) and (5) define the maximum and minimum output limits for each gas and steam turbine unit, respectively.

$$\begin{array}{c}\hfill \underline{G_c}·u_{tc}\le g_{tc}\le \overline{G_c}·u_{tc}\forall t,c\end{array}$$

(4)

2.2.3. Steam Turbines’ Operational Limits

The minimum output limit outlined in Equation (5) is associated with Equation (32), which equals 1 to accommodate either a hot or cold startup of a steam turbine. Following this logic, for a hot or cold startup of a steam turbine, the term $(u_{ts}-y_{ts})$ equals zero. This allows the output of the steam turbine to match the output value determined for each respective thermal status startup.

$$\begin{array}{c}\hfill \underline{G_s}·(u_{ts}-y_{ts})\le g_{ts}\le \overline{G_s}·u_{ts}\forall t,s\end{array}$$

(5)

2.2.4. Combined-Cycle Operational and Commitment Limits

$$\begin{array}{c}\hfill \underline{GCC}·c{c}_t\le gc{c}_t\le \overline{GCC}·c{c}_t\forall t\end{array}$$

(6)

$$\begin{array}{c}\hfill c{c}_t-{\displaystyle \sum _{s\in NS}}{u}_{ts}\le 0\forall t\end{array}$$

(7)

$$\begin{array}{c}\hfill MUG·c{c}_t-\sum _{c\in NC}{u}_{tc}\le 0\forall t\end{array}$$

(8)

$$\begin{array}{c}\hfill MUG·u_{ts}\le \sum _{c\in NC}{u}_{tc}\forall t,s\end{array}$$

(9)

$$\begin{array}{c}\hfill \sum _{s\in NS}{u}_{ts}\le \sum _{c\in NC}{u}_{tc}\forall t\end{array}$$

(10)

$$\begin{array}{c}\hfill v{h}_t\le \underline{GCC}·(1-c{c}_t)\forall t\end{array}$$

(11)

Constraints (7) through (10) define the minimum number of gas units required for a coupled operation in a combined cycle with steam turbines. Specifically, the inequality (10) ensures that the number of online steam units does not exceed the number of gas turbine units.

Constraint (7) ensures that only one steam unit will be dispatched if the plant is operating in combined-cycle mode. Lastly, the constraints (8) and (9) ensure that a combined-cycle operation is only possible when the minimum number of gas turbines (MUG) is dispatched.

2.3. Constraints Related to The Combined-Cycle Operation

The following constraints represent the combined-cycle operation associated with the coupled operation between gas turbines and steam turbines.

2.3.1. CCGT Output Constraints

$$\begin{array}{c}\hfill gc{c}_t+v{h}_t=\sum _{c\in NC}{g}_{tc}+\sum _{s\in NS}{g}_{ts}-aux{g}_t\forall t\end{array}$$

(12)

$$\begin{array}{c}\hfill gcc{t}_t=gc{c}_t+dd{e}_t+nd{e}_t\forall t\end{array}$$

(13)

Equations (12) and (13) represent the total generation of the CCGT including the contribution of the gas turbines ${\sum }_{c\in NC}{g}_{tc}$ and the steam turbines ${\sum }_{s\in NS}{g}_{ts}$ and subtracting the auxiliary energy consumption $aux{g}_t$.

On the other hand, Equation (13) describes the total output of the CCGT $gcc{t}_t$. This is equal to the output of the CCGT $gc{c}_t$ when operating in combined-cycle mode, plus the startup and shutdown output ramps ($dd{e}_t$ and $nd{e}_t$, respectively).

2.3.2. Steam–Gas Coupling Operation Constraints

$$\begin{array}{cc}\hfill \sum _{s\in NS}{g}_{ts}+gvs{c}_t=STF·\sum _{c\in NC}{g}_{tc}+\sum _{c\in NC}a{f}_{ts}& \forall t\hfill \end{array}$$

(14)

$$\begin{array}{c}\hfill a{f}_{ts}\le PAF·u_{tc}\forall t,s\end{array}$$

(15)

$$\begin{array}{c}\hfill gvs{c}_{ts}\le \sum _{s\in NS}\overline{G_s}·(u_{ts}+\overline{u_{ts}})\forall t\end{array}$$

(16)

$$\begin{array}{c}\hfill \sum _{s\in NS}{g}_{ts}\ge \overline{u_{ts}}·\sum _{s\in NS}\overline{G_s}\forall t\end{array}$$

(17)

Equation (14) relates the energy produced by steam turbines as a function of the energy produced by gas turbines. The STF (For the sake of simplifying our model, we assumed a constant steam-to-gas output factor, directing our focus towards resolving the operational constraints linked to the steam turbines. Indeed, a more-precise approach accounting for variable steam-to-gas output ratios could be proposed as an area for future research). relates the amount of energy produced for each MW produced by the gas turbine. Also, the energy produced by the steam turbines depend on the status of the additional fire $af$, which depends on the status of the gas turbine, i.e., there are supplementary fires for each gas turbine that can increase the output of the steam turbine unit until the maximum capacity of the unit. Each additional fire is represented by the inequality (14). Also, for a certain level of load, the control valves of the steam turbines are fully open, and the HRSGs are then used to increase or decrease the pressure rating in the turbines.

In the same Equation (14), the steam waste $gvs{c}_t$ is represented, which represents the amount of energy that is no longer delivered when it is not possible to take advantage of all the steam generated by the gas turbines. This variable is defined by Equation (16), and its quantity is limited by the maximum capacity that a steam turbine can deliver, as long as it is online.

2.3.3. CCGT Auxiliary Consumption

$$\begin{array}{cc}\hfill & aux{g}_t=AUXCC·c{c}_t+AUXGT·\sum _{c\in NC}{u}_{tc}+\hfill \\ \hfill & AUXST·\sum _{s\in NS}{u}_{ts}\forall t\hfill \end{array}$$

(18)

Equation (18) represents the maximum and minimum limits of auxiliary consumption $aux{g}_t$, both for steam turbines and for gas turbines. It is important to note at this juncture that auxiliary consumption is considered only when each unit, be it a gas or steam turbine, is in operation. For this reason, each value of auxiliary consumption is multiplied by the on/off binary variable of each unit. The same principle applies to the combined-cycle mode operation $c{c}_t$.

2.4. Minimum Uptime and Downtime Constraints

The following constraints represent the minimum uptime and downtime for each gas and steam unit and CCGT, adopted from [4,5]. Equation (19) describes the startup and shutdown processes for each gas turbine, steam turbine, and the entire CCGT plant.

2.4.1. Gas Turbines Minimum Uptime and Downtime Constraints

Delving into the details of Equation (19), one can observe the following: If the variable $u_{t,c/s/cc}$ transitions from 0 at period $t-1$ to 1 at period t, this indicates a startup of the unit or plant, denoted as $y_{t,c/s/cc}=1$. Conversely, if the unit or plant transitions from $u_{t-1,c/s/cc}=1$ in $t-1$ to $u_{t,c/s/cc}=0$ at period t, the equation signifies a shutdown, represented by $z_{t,c/s/cc}=1$.

$$\begin{array}{c}\hfill y_{t,c/s/cc}-z_{t,c/s/cc}=u_{t,c/s/cc}-u_{t-1,c/s/cc}\forall t,c,s,cc\end{array}$$

(19)

For the first period, variable $u_{t-1}$ has to be replaced in constraint (19) by $OnOff$, considering the initial conditions of the gas and steam turbines.

2.4.2. Startup and Shutdown Mutually Exclusive Variables

$$\begin{array}{c}\hfill y_{t,c/s/cc}+z_{t,c/s/cc}\le 1\forall t,c,s,cc\end{array}$$

(20)

Equation (20) ensures that a startup and shutdown do not occur simultaneously within the same period.

2.4.3. Status to Reach The Minimum Uptime/Downtime

$$\begin{array}{c}\hfill u_{t,c/s/cc}=OnOfft=1;\forall c,s,cc,\\ \hfill Lupmin+Ldownmin>t\end{array}$$

(21)

Equation (21) sets the variable $u_{t,c/s/cc}$ to the value of the initial status at period 1 if the minimum uptime and the minimum downtime have not been met. To clarify, if the unit was “on” during the last period of the previous day and has not yet fulfilled the required number of periods in that state, it must remain online until it completes the minimum uptime. The same logic applies when the unit is in the “off” condition. These conditions are effectively addressed by the constraints (22) and (23), which are detailed in the subsequent sections.

2.4.4. Minimum Uptime and Downtime Condition Constraint

$$\begin{array}{c}\hfill u_{t,c/s/cc}\ge \sum _{c,s\in N}{y}_{i,c/s/cc}\forall t,c,s,cc\\ \hfill i\ge t-UT+CountOn·max(0,2-i)+1,\\ \hfill i\le t\end{array}$$

(22)

Equation (22) sets the variable $u_{t,c/s/cc}$ to be greater than 0 to ensure the unit meets the minimum uptime, UT.

$$\begin{array}{c}\hfill 1-u_{t,c/s/cc}\ge \sum _{i\in \ll }{y}_{i,/c/s/cc}\forall t,c,s,cc,\\ \hfill i\ge t,\\ \hfill i\ge (t-DT+1)\end{array}$$

(23)

Otherwise, Equation (23) sets the left side of the equation $1-u_{t,c/s/cc}$ to be greater than 0 to ensure the unit meets the minimum downtime, DT.

2.5. Ramps’ Constraints

The following equations represent the ramping constraints for the CCGT. Constraints (24)–(30) represent the startup ($dd{e}_t$) and shutdown ($nd{e}_t$) constraints proposed by [24]. As previously noted, the Colombian ISO adheres to the startup and shutdown guidelines outlined in [7]. According to these guidelines, thermal power plants are required to follow a predetermined startup sequence, using hourly energy blocks that progress from 0 to the plant’s minimum power output. Similarly, during shutdown, the plant should adhere to the declared hourly energy blocks, descending from the minimum generation output to 0. It is essential to specify that these startup ramps are determined based on the thermal status of the plant, namely: hot, warm, and cold conditions, which are described as follows.

Two types of overlaps between states arise during the application of startup ramps, especially in periods adjacent to those defining each state. The first overlap concerns the change in thermal state amidst the startup ramp process. For instance, while applying fixed increase blocks for the ‘hot’ condition, the resource might transition to the ‘warm’ state. The second overlap stems from the varied count of fixed increase blocks assigned to each thermal state, leading to a situation where the commencement of fixed blocks for the succeeding state overlaps with the current state. This is evident when, to implement the fixed augmentation blocks for the ‘cold’ state, initiation in the ‘warm’ state becomes necessary.

2.5.1. Hot Startup Ramp Constraints

For periods when the CCGT is deemed to be in a ‘hot’ thermal state ($t+toff\le te1$), the programmed generation conforming to the startup ramp is governed by Equation (24). This constraint integrates two components: the first delineates the ‘hot’ energy blocks, while the second encompasses the ‘warm’ fixed energy blocks that need to be factored in during a ‘hot’ startup process for the CCGT.

$$\begin{array}{c}\hfill dd{e}_t={\displaystyle \sum _{i=1}^{nSe1}}Se1·y_{t+i}+{\displaystyle \sum _{j=1}^{nSe2}}Se2·y_{t+j}\\ \hfill \forall t+toff\le te1\end{array}$$

(24)

$$\begin{array}{c}\hfill y_t=0toff\le te1;t\le nSe1\end{array}$$

(25)

If the CCGT begins in an offline state and the conditions for a ‘hot’ startup ramp remain viable (i.e., $toff\le te1$), the plant can initiate a startup during periods when those fixed energy blocks are achievable. To facilitate this, the constraints (25) ensure the variable $y_t=0$ for periods where the CCGT cannot fulfill all the requisite startup energy blocks.

2.5.2. Warm Startup Ramp Constraints

For periods when the CCGT is deemed to be in a ‘warm’ thermal state ($te1<t+toff\le te2$), the programmed generation conforming to the startup ramp is governed by Equation (26). This constraint integrates three components: the first delineates the ‘hot’ energy blocks that began to be applied in the ‘warm’ state; the second represents the fixed augmentation blocks to make the ‘warm’ startup; the third represent the ‘cold’ fixed energy blocks that need to be factored in during a ‘warm’ startup process for the CCGT.

$$\begin{array}{c}\hfill dd{e}_t={\displaystyle \sum _{i=1}^{nSe1}}Se1·y_{t+i}+{\displaystyle \sum _{j=1}^{nSe2}}Se2·y_{t+j}+\\ \hfill {\displaystyle \sum _{k=1}^{nSe3}}Se3·y_{t+k}\forall te1<t+toff\le te2\end{array}$$

(26)

$$\begin{array}{c}\hfill y_t=0\forall te1<t+toff\le te2;t\le nSe2\end{array}$$

(27)

If the CCGT begins in an offline state and the conditions for a ‘warm’ startup ramp remain viable (i.e., $te1<toff\le te1$), the plant can initiate a startup during periods when those fixed energy blocks are achievable. To facilitate this, the constraints (27) ensure the variable $y_t=0$ for periods where the CCGT cannot fulfill all the requisite startup energy blocks.

2.5.3. Cold Startup Ramp Constraints

For periods when the CCGT is deemed to be in a ‘cold’ thermal state ($t+toff>te2$), the programmed generation conforming to the startup ramp is governed by Equation (28). This constraint integrates two components: the first delineates the ‘warm’ energy blocks that began to be applied in the ‘cold’ state, and the second represents the fixed augmentation blocks to make the ‘cold’ startup.

$$\begin{array}{c}\hfill dd{e}_t={\displaystyle \sum _{j=1}^{nSe2}}Se2·y_{t+j}+\\ \hfill {\displaystyle \sum _{k=1}^{nSe3}}Se3·y_{t+k}\forall t+toff>te2\end{array}$$

(28)

$$\begin{array}{c}\hfill y_t=0\forall toff>te2;t\le nSe3\end{array}$$

(29)

If the CCGT begins in an offline state and the conditions for a ‘cold’ startup ramp remain viable (i.e., $toff>te2$), the plant can initiate a startup during periods when those fixed energy blocks are achievable. To facilitate this, the constraints (29) ensure the variable $y_t=0$ for periods where the CCGT cannot fulfill all the requisite startup energy blocks.

2.5.4. Shutdown Ramp Constraint

Equation (30) characterizes the shutdown process using fixed energy blocks.

$$\begin{array}{c}\hfill nd{e}_t=\sum _{i=1}^{nSr}S{r}_i·zcct-i+1\forall tt\le i\end{array}$$

(30)

2.5.5. Operating Ramp Constraints

Constraint (31) represent the ramping constraint between two consecutive hours.

$$\begin{array}{c}\hfill -RD\le gcc{t}_t-gcc{t}_{t-1}\le RU\forall t\end{array}$$

(31)

Equation (31) stipulates that the variation in the output from one hour to the next is constrained to RU and RD for upward and downward changes, respectively. Here, RU and RD represent the maximum block of energy that the CCGT can increase or decrease from one hour to the next, based on the configuration or mode selected a priori for operation.

2.6. Steam Turbine Startup

Constraint (32) introduces two binary variables that distinguish between a cold (C) and hot (H) start of the steam turbine. Constraint (33) stipulates that a hot start cannot be performed unless the steam turbine has been dispatched within the preceding 9 h.

Constraint (34) asserts that a hot start for the steam turbine is not possible unless at least one of the gas turbines was connected in the previous hour.

$$\begin{array}{c}\hfill y_{t,s}=y_{t,s}^C+y_{t,s}^H\forall t,s\end{array}$$

(32)

$$\begin{array}{c}\hfill y_{t,s}^H\le A{A}_{t,s}+\sum _{j=t-8}^{t-1}{u}_{j,s}\forall t,s\end{array}$$

(33)

$$\begin{array}{c}\hfill y_{t,s}^H\le \sum _c{u}_{t-1,c}\forall t>1,s\end{array}$$

(34)

$$\begin{array}{c}\hfill y_{t,s}^H\le \sum _{c}onof{f}_c^{t_0}\forall t=1,s\end{array}$$

(35)

$$\begin{array}{c}\hfill B{B}_{t,c}+\sum _{j=t-2}^t{u}_{j,c}\ge KMH·b_{t,c}\forall t,c\end{array}$$

(36)

$$\begin{array}{c}\hfill y_{t,s}^C\le \sum _c{b}_{t,c}\forall t,c\end{array}$$

(37)

$$\begin{array}{c}\hfill y_{t,s},y_{t,s}^C,y_{t,s}^H,b_{t,c}\in \{0,1\}\forall t,s,c\end{array}$$

(38)

Constraint (36) introduces an auxiliary binary variable, $b_{t,c}$, which is set to 1 if a gas turbine has been connected during the last $KMH$ hours and 0 otherwise. Constraint (37) models the condition that a cold start (or any startup, for that matter) cannot be performed unless at least one of the gas turbines has been connected for the past 6 h. Finally, the constraint (38) designates the new variables as binary.

The startup of a steam turbine depends on the unit’s temperature in each period, which, in turn, depends on the number of hours the unit has been offline. Once the number of hours that the steam turbine is in a cold state (denoted as $KST$) is defined, it is necessary to determine the number of hours that must have elapsed since the first gas turbine started. This is crucial to ensure the temperature required for the steam turbine startup is reached. We define $KGC$ for a cold startup and $KGH$ for a hot startup.

The following constraints depict how the steam turbine should be started, taking into account the unit’s state (cold or hot startup) as explained above.

Another crucial aspect to consider for the steam turbine startup is the power output in the first time period. For a cold startup, the power output should be at a low load to prevent thermo-mechanical damage. The subsequent constraints represent the power output of the steam turbine unit at startup, taking into account the type of startup, i.e., cold startup ($GTSC$) or hot startup ($GTSH$):

$$\begin{array}{c}\hfill g_{t,s}\le GSTC·y_{t,s}^C+\\ \hfill GSTH·y_{t,s}^H+{\overline{G}}_s·(u_{t,s}-y_{t,s})\forall t,s\end{array}$$

(39)

$$\begin{array}{c}\hfill g_{t,s}\ge GSTC·y_{t,s}^C+GSTH·y_{t,s}^H\forall t,s\end{array}$$

(40)

Figure 2, adopted from [13], summarizes the previous steam turbines startup explained in the set of equations from (32) to (40).

2.7. Load Distribution between Gas Turbines’ Constraints

Load Distribution Constraint

The following set of constraints capture the fact that CCGT operators strive to achieve a similar power output among the gas turbines that operate above their technical minimum. For that purpose, we introduce the following constraints, which establish the absolute value of the difference in power output between two different gas turbines.

$$\begin{array}{c}\hfill gcd{r}_{tc}-gcd{r}_{tcc}\le {\Delta }_{t,c,cc}\forall t,cc\ne c\end{array}$$

(41)

$$\begin{array}{c}\hfill {\Delta }_{t,c,cc}\ge 0\forall t,cc\ne c\end{array}$$

(42)

Note that this difference $gcd{r}_{tc}-gcd{r}_{tcc}$ could be either positive or negative; however, if $gcd{r}_{tc}-gcd{r}_{tcc}$ is negative, then $gcd{r}_{tcc}-gcd{r}_{tc}$ will be positive, yielding a positive lower bound for ${\Delta }_{t,cc,c}$. We leave one of the sums on purpose in order to not count the power deviation twice.

This difference in power output could be penalized in the objective function value; however, we can only take it into account if both gas turbines are actually above the technical minimum (It should not be taken into account when one gas turbine is off, and the other one is above the technical minimum,). Hence, we define binary variable ${\delta }_{t,c,cc}$, which takes a value of 1 if both gas turbines are above the technical minimum:

$$\begin{array}{c}\hfill u_{tc}+u_{tcc}\le 1+{\delta }_{t,c,cc}\forall t,cc\ne c\end{array}$$

(43)

$$\begin{array}{c}\hfill {\Delta }_{t,c,cc}-\overline{GCC}·(1-{\delta }_{t,c,cc})\le vhd{r}_{t,c,cc}\forall t,cc\ne c\end{array}$$

(44)

$$\begin{array}{c}\hfill {\delta }_{t,c,cc}\in \{0,1\}\forall t,cc\ne c\end{array}$$

(45)

$$\begin{array}{c}\hfill \overline{GCC}\ge vhd{r}_{t,c,cc}\ge 0\forall t,cc\ne c\end{array}$$

(46)

However, when the load distribution among the gas turbines is unequal, the steam delivered through the common collector exhibits varying thermal characteristics. This leads to temperature differences between the surface and the center of the steam turbine rotors, which can cause long-term damage.

In practical applications, maintaining an equal load distribution may not always be feasible. For instance, depending on the grid’s condition, the ISO may require more output power from a specific gas turbine unit than from others. Despite this, the ISO strives to adhere to this constraint as closely as possible during operations.

3. Results and Discussion

This section presents the analysis of two distinct numerical case studies:

(I) The first case involves a CCGT that is operational at the end of the previous day and, then, shuts down at the beginning of the analyzed day, only to restart and reach the maximum output capacity by the day’s end.

(II) The second case examines a scenario where the CCGT is initially offline on the analyzed day, but starts up later, eventually achieving maximum output capacity by the end of the day.

The input parameters for both cases are detailed in

Table 1. CCGT parameters.

Variable	Value	Unit
$\overline{GCC}$	800	MW
$\underline{GCC}$	210	MW
$PAF$	15	MW
$AUXCC$	5	MW
$AUXGT$	0.45	MW
$AUXST$	2	MW
$RD/RU$	335	MWh
$PCC$	120	USD/MWh
$PBC$	500	USD/MWh
$CSC$	15,000	USD
$MUG$	2	p.u.
$STF$	0.613	p.u.
$NC$	5	p.u.
$NS$	2	p.u.
$t1$	t $<=$ 16	Hours
$t2$	16 < t $<=$ 30	Hours
$t3$	t > 30	Hours
$KGC$	3	Hours

. Additionally, the startup and shutdown ramping processes for these scenarios are outlined in

Table 2. Startup and shutdown ramps.

Hour	H-Startup	W-Startup	C-Startup	Shutdown
H1	50	50	50	210
H2	100	100	100	100
H3	150	100	100	50
H4	210	150	100	0
H5	0	210	150	0
H6	0	0	210	0

.

For the cases simulated in this work, $NC$ gas turbines and $NS$ steam turbines have the same characteristics. Figure 3 shows the general scheme of the CCGT plant considered for the simulations.

Table 3. Gas turbines.

Variable	Value	Units
$\overline{G}$	100	MW
$\underline{G}$	50	MW
TC	5	MW/min
TD	5	MW/min

and

Table 4. Steam turbines.

Variable	Value	Units
$\overline{G}$	170	MW
$\underline{G}$	80	MW
$GSTH$	80	MW
$GSTC$	30	MW

show the characteristics of the gas and steam turbine units.

3.1. Case I

The initial conditions for the CCGT are outlined in

Table 5. Initial conditions of units—Case I.

Unit	$\mathit{ton}/\mathit{off}$ (Hours)	$\mathit{Gt}0$ (MW)
GT1	8	67
GT2	0	0
GT3	0	0
GT4	0	0
GT5	8	67
ST1	8	83
ST2	0	0

. In Figure 4, the actual CCGT operation (discontinuous red line) is shown to deviate from the initial dispatch L due to the ramping constraints. Specifically, the model initiates a shutdown ramp in period 5, rendering the CCGT offline until period 15, followed by a hot startup ramp. Notably, the model ramps up in period 19, achieving the maximum CCGT capacity from periods 20 to 22, as necessitated by the constraint (31).

Figure 5 displays the unitwise dispatch and auxiliary consumption. The model maintains units GT1, GT5, and ST1 online to meet the initial dispatch requirements. For starting up, it activates units GT3, GT5, and GT11, with ST2 undergoing a cold startup, necessitating at least 3 h of operation from the gas units.

Adherence to a heuristic-model-based dispatch (yellow line) leads to significant deviations from scheduled production, primarily due to the misrepresentation of the CCGT’s operational rules. Realistically, the CCGT cannot accurately follow such heuristic dispatches because of technical constraints like the operational hours’ relationship between gas and steam turbines. As per [25], deviations exceeding five percent between the heuristic and actual generation are penalized. Assuming the model’s results (discontinuous red line) represent actual generation and the penalty cost is $PCC$, as per [25], this would result in a daily penalty of USD 60,957 for Case I. Such a penalty, recurring daily in scenarios like those depicted in Figure 6, represents a substantial financial burden. Furthermore, these deviations necessitate the use of costly reserves by the system operator to balance the grid. Employing a more-realistic model, like the one proposed in this paper, can prevent these costly penalties and reserve activations by enabling more-accurate dispatch decisions.

To achieve maximum capacity during periods 20 to 22, the model strategically dispatches all units, utilizing supplementary fires to ensure the steam turbines deliver the required energy. This highlights the model’s capability to optimize unit-specific operations, a significant advancement over the heuristic simulation model. The heuristic approach, in contrast, represents the CCGT output in aggregate, failing to detail the individual contributions of each gas and steam turbine. This underscores the superiority of the proposed SEUC optimization model, which offers a more-granular and -realistic representation of CCGT operations, crucial for accurate energy management and decision-making.

3.2. Case II

In this scenario, all units begin offline, as indicated in

Table 6. Initial conditions of units—Case II.

Unit	$\mathit{ton}/\mathit{off}$ (Hours)	$\mathit{Gt}0$ (MW)
GT1	8	0
GT2	8	0
GT3	8	0
GT4	8	0
GT5	8	0
ST1	8	0
ST2	8	0

. The results for Case II, shown in Figure 6, revealed that a warm startup is needed from periods 14 to 18 to meet the initial dispatch L. As in Case I, an increased ramp in period 19 is crucial to achieve maximum capacity from periods 20 to 22. Contrarily, the heuristic model opts for a hot startup, neglecting the prior state of the units, which could lead to long-term equipment damage and plant unavailability, thus affecting power system reliability (see Figure 1).

Calculating the penalties for Case II, as we did in Case I, we found daily penalties amounting to USD 66,093 due to discrepancies between the heuristic schedule and the actual CCGT operation. These penalties can significantly impact the power plant’s profitability and system operations, depending on the frequency of such startups annually.

As depicted in Figure 7, the startup ramp involved units GT2, GT3, and GT4, with ST2 starting in cold mode in period 18. To achieve the CCGT’s maximum capacity, all units were utilized at full output, including one supplementary fire, contributing an additional 4.75 MW.

4. Conclusions

In this paper, we proposed a self-unit-commitment model to optimize the dispatch of a CCGT with real operational constraints. This approach seeks to improve upon heuristic models in use currently by CCGT operators in the Colombian electric power system. We highlighted the importance for ISOs and generation companies to employ precise tools, as the model proposed in this paper, when planning operating decisions, as they can avoid economic penalties for CCGT operators, as well as deviations from the scheduled output for the ISO.

Apart from other standard constraints in the literature, we proposed an original formulation for individual gas and steam turbine units that guarantee specific characteristics of the steam. Those characteristics, which are actually in place in large CCGTs, are necessary to minimize the impact of the thermo-mechanical fatigue produced by the energy output changes required by the system operator. Employing a realistic optimization model helps to increase the useful time of the CCGT units and the reliability of the CCGT, minimizing future failures and avoiding penalties due to deviations to the program, once the output of the proposed model in this work can be followed by the CCGT in real-time. We also proposed a novel operating constraint that allows for an even load distribution among individual gas turbines—a constraint that is being imposed in real-life CCGTs.

In future research, we want to extend this work from a self- to a full-unit commitment, considering all power plants of the system. Such a model would help the ISO improve the solution of the dispatch in the Colombian power system, where CCGT plants play an important role.