Optimization of Power Plant Operation Scheduling in Tahuna 1 Isolated System to Minimize the Levelized Cost of Electricity


On this article

Gilang Cahyo Nugroho<sup>1.2</sup> and Poetro Lebdo Sambegoro<sup>1</sup>

<sup>1</sup>Faculty of Mechanical and Aerospace Engineering, Institut Teknologi Bandung <sup>2</sup>PT Perusahaan Listrik Negara (Persero), Bandung, Indonesia <sup>1</sup>Faculty of Mechanical and Aerospace Engineering, Institut Teknologi Bandung, Bandung, Indonesia gilang.cahyon@pln.co.id, poetrols@itb.ac.id

Abstract: The Tahuna 1 System which is located on Sangihe Island, North Sulawesi, Indonesia consists of 4 power plants, namely the Tahuna Diesel Power Plant, Tamako Diesel Power Plant, Sangihe Solar Power Plant, and Ulung Peliang Micro Hydro Power Plant are synchronized through the 20 kV medium voltage network with a total net power capacity of 15,358 kW. The cumulative value of the Levelized Cost of Electricity (LCOE) of the Tahuna 1 System in 2023 reached Rp4,210.86/kWh, because there is no normal operation scheduling for the generating units per hour as a reference in operating the system. Until now, the dispatcher of the Tahuna Command Center must coordinate with the power plant operator to operate the generating units without any priority order of operation, but only the availability and readiness of each unit. In this research, a scheme for hourly generating unit operation scheduling was designed on the Tahuna 1 System to minimize the LCOE by considering the system load profile, generating unit capacity, and spinning reserve with a modification of the Merit Order and Priority Listing method to increase the selectivity of determining the priority list of generating unit operations based on the cheapest LCOE ranking (Rp/kWh) according to the range of generating zones created from the intersection of the cost function equation between all diesel generating units. The design of the hourly operation scheduling scheme of generating units was carried out on 12 units of Tahuna Diesel Power Plant, and 6 units of Tamako Diesel Power Plant as independent variables. Based on the results of the operation scheduling scheme with consideration of two spinning reserve units, there is a decrease in the LCOE value of the Tahuna 1 System from the initial Rp4,210.86/kWh to Rp2,913.66/kWh or 30.81% more economical compared to the current conditions. This scheduling scheme is also a better choice to implement when viewed from the operational reliability of the 20 kV distribution system.

Keywords: isolated system; power plant operation scheduling; levelized cost of electricity

1. Introduction

Tahuna (3.6111059, 125.4751878) is the capital of Sangihe Islands Regency located on Sangihe Island, North Sulawesi, Indonesia. There are twelve sub-districts on Sangihe Island with a total area of 692.67 km² [1]. Sangihe Island has a population of 130,458 people as of 2023 [1], with a population growth ratio of 0.29% from 2020 to 2023 [1]. The electricity system on this island is categorized as an isolated system because it is not yet connected to the Minahasa System which is the grid of North Sulawesi and Gorontalo. The Tahuna 1 System consists of 4 power plants, namely the Tahuna Diesel Power Plant, Tamako Diesel Power Plant, Sangihe Solar Power Plant, and Ulung Peliang Micro Hydro Power Plant which are synchronized through the 20 kV medium voltage network with a total power capacity of 15,358 kW [2]. During 2023, the peak load of the Tahuna 1 System was 5,710 kW which occurred at 07:00 PM [3] because many household consumers use electricity for lighting simultaneously.

The cumulative LCOE value of the Tahuna 1 System in 2023 reached Rp4,210.86/kWh [4] or 81.65% higher than the average LCOE of the Northern Sulawesi region, which was Rp2,318/kWh [4]. This is because there is no normal operation scheduling for the generating units per hour as a reference in operating the system. Until now, the dispatcher of the Tahuna Command Center must coordinate with the power plant operator to operate the generating units without any priority order of operation, but only the availability and readiness of each unit.

Received: August 9th, 2025. Accepted: December 16th, 2025

DOI: 10.15676/ijeei.2025.17.4.9

Research on the scheduling of thermal generating unit operations has been conducted throughout the world. In previous research, Ananda R et al. [5] analyzed the operation scheduling of 4 geothermal generating units with maximum power capacities varying from 32 MW to 55.5 MW to meet the system's power needs for 24 hours with a system load of 159.59 MW to 163.79 MW using the Forward Dynamic Programming method, the results obtained were a reduction in operating costs from Rp 4,885,715 to Rp3,036,492 or 30% more economical compared to current conditions. Another research by Raj, Vinod and Saurabh Chanana [6] discussed the scheduling of 10 thermal generating unit operations with varying maximum power capacities from 55 MW to 455 MW to meet the system load from 700 MW to 1,500 MW using the Priority Listing and Lagrange Relaxation methods, the results obtained there is a decrease in the value of total generation cost of $1,480 in the daily period using the Lagrange Relaxation method compared to the Priority Listing method which was initially from $566,760 to $565,280. Meanwhile Hanafi, Ikhsan and Rinaldy Dalimi [7] conducted research related to the scheduling of operations of 6 thermal generating units using the Merit Order method with variations in maximum power capacity from 32 MW to 650 MW to supply the system with load fluctuations from 1,317 MW to 1,774 MW, a decrease in the value of generation costs in the daily period was obtained from the initial $1,988,410 to $1,696,770.

From several research that have been described above, it was found that each design of the generating unit operation schedule aims to minimize the generation cost which consists of objective function equations and constraint functions based on existing parameters. In this research, a scheme for hourly operation scheduling of generating units in the Tahuna 1 System will be designed to minimize LCOE by considering the system load profile, generating unit capacity, and spinning reserve with a modification of the Merit Order and Priority Listing method to increase the selectivity of determining the priority list of generating unit operations based on the cheapest LCOE ranking (Rp/kWh) according to the range of generating zones created from the intersection of the cost function equation between all diesel generating units. This method is also used because it can be formulated using a formula with an appropriate algorithm.

1. Methodology

To optimize the operating pattern of all generating units in the Tahuna 1 System, a scheduling scheme is required by determining the on and off status and the loading value of generating units by considering the system load profile, generating unit capacity, and spinning reserve based on performance data from each generating unit. To assist the calculation and data processing process starting from the regression of Specific Fuel Consumption (SFC) performance data and fuel consumption of each generating unit, the formulation of generating unit cost functions, generation costs ramping starting from minimum to maximum capacity of each diesel generating unit as an independent variable, until the final hourly operation scheduling scheme in this research, Microsoft Excel software is used by compiling several formulas in each process. The output results of this research will be in the form of a file with the format (*.xlxs) from Microsoft Excel so that it can be used by all stakeholders easily without any additional plugins to be installed, starting from dispatcher of the Tahuna Command Center to operators at each power plant in determining the priority order of unit operations in real time according to the estimated hourly system load. The process of working on this research is explained in the research flow diagram as shown in Figure 1 as follows.

1

Figure 1. Research flow diagram

A. System Load Profile

One important aspect in designing a generating unit operation scheduling scheme is to know the average trend of the system load per hour within a certain period. This aims to ensure that the designed scheme can fulfill the system's hourly power needs for 24 hours. In this research, only data for the year 2023 was used to represent the hourly system load profile in the Tahuna 1 System. This is because in the previous year there were no additional generating units in the Tahuna Diesel Power Plant (Unit #10, Unit #11, and Unit #12) and Tamako Diesel Power Plant (Unit #5, Unit #6, and Unit #7), so that optimal scheduling could not be made. In Figure 2 below is the average system load profile per hour of the Tahuna 1 System period of 2023 [3].

5

Figure 2. Average system load profile per hour of the Tahuna 1 System 2023 [3]

From Figure 2, it is explained that the average load on the Tahuna 1 System for the period of 2023 is 4,578 kW with the peak load value of 5,710 kW at 07:00 PM [3]. The peak load of the system occurs at that hour because many household consumers use electricity for lighting simultaneously. It is also known from Figure 2 that the order of the largest composition of the Tahuna 1 System load is borne by the Tahuna Diesel Power Plant, then continued by the Tamako Diesel Power Plant, Sangihe Solar Power Plant, and Ulung Peliang Micro Hydro Power Plant.

Based on Figure 2, it can also be seen that the operating hours of the Sangihe Solar Power Plant only occur at 06:00 AM-05.00 PM along with the availability of solar irradiation. Based on direct measurement data at the location, during 2023, the average daily irradiation at the Sangihe Solar Power Plant reached 4.49 kWh/m2 /day with the highest solar irradiation occurring at 11:00 AM with a value of 0.604 kW/m2 [8]. Since the Sangihe Solar Power Plant and Ulung Peliang Micro Hydro Power Plant are renewable energy-based power plants which cause the system load to be unpredictable (intermittent), the average load of these power plants in this research is used as a control variable.

B. Existing Generating Unit

The Tahuna 1 System consists of the Tahuna Diesel Power Plant, Tamako Diesel Power Plant, Sangihe Solar Power Plant, and Ulung Peliang Micro Hydro Power Plant which are synchronized through the 20 kV medium voltage with a total installed power of 21,101 kW. However, due to age factors, the total net power that can be generated by all these generating units is 15,358 kW as described in Table 1 below [2].

Table 1. Generating units installed on Tahuna 1 System [2]

Power Capacity (kW)
NoGenerating UnitsInstalledMinimumNetto
#1Unit #1 Diesel Generator (DG)
Tahuna
765145425
#2Unit #4 DG
Tahuna
50035150
#3Unit #5 DG Tahuna50050300
#4Unit #6 DG Tahuna1,300168850
#5Unit #7 DG Tahuna1,3001961,000
#6Unit #8 DG Tahuna1,120377850
#7Unit #9 DG Tahuna1,120450850
#8Unit #10 DG Tahuna
(Ex UIW MMU)
1,292207904
#9Unit #11 DG Tahuna
(Ex UIW MMU)
1,292378904
#10Unit #12 DG Tahuna (Ex UIW MMU)1,292202904
#11Unit #1 DG Tahuna (KIP 1)1,6002251,100
#12Unit #1 DG Tahuna (KIP 2)3,7202032,250
#13Unit #1 DG Tamako700172630
#14Unit #2 DG Tamako50023450
#15Unit #4 DG Tamako1,3002551,170
#16Unit #5 DG Tamako
(Genset Mobile)
20020180
#17Unit #6 DG Tamako
(Genset Mobile)
20040180
#18Unit #7 DG Tamako
(Genset Mobile)
1002090

Table 1. Continued

#19Unit #1 Micro Hydro (MH) Ulung Peliang
(Upel)
1,0001201,000
#20Unit #1 Photovoltaic (PV) Sangihe1,300141,170
Total (kW)21,101-15,358

From Table 1, it is known that the Tahuna Diesel Power Plant consists of 12 generating units, the Tamako Diesel Power Plant consists of 6 generating units, and the Ulung Peliang Micro Hydro Power Plant and Sangihe Solar Power Plant each consist of one generating unit.

In this research, the design of the hourly generating unit operation scheduling scheme will only be carried out on the Tahuna Diesel Power Plant and Tamako Diesel Power Plant as independent variables. While the Sangihe Solar Power Plant and Ulung Peliang Micro Hydro Power Plant are used as control variables, because the loading of these generating units on the Tahuna 1 System depends on natural conditions so that it cannot be predicted (intermittent). Figure 3 below is the Single Line Diagram (SLD) of the Tahuna 1 System [9].

5

Figure 3. Single line diagram 20 kV of the Tahuna 1 System [9]

C. Economic Dispatch Modeling

1) Generation Cost Function

In designing a generating unit scheduling scheme, there are two main topics, namely determining the generating unit to be operated and the total generation costs of the operating unit [10]. In the process of designing the operating scheduling scheme for generating units, the generation cost function of each thermal generating unit is used in the form of a secondorder polynomial as in (1) below [11].

\[F_i(P_{i,t}) = a_i + b_i P_{i,t} + c_i P_{i,t}^2 \tag{1}\]

where:

Fi(Pi,t) = cost function of generating unit-i

ai, bi, ci = coefficient of cost function of generating unit-i

Pi,t = generated power of generating unit-i at hour-t (kW)

Based on the results of the regression of SFC performance data and fuel consumption of each generating unit, by referring to (1), the cost function for each installed diesel generating units of the Tahuna 1 System can be obtained as in Table 2.

Table 2. Generating unit cost function in Tahuna 1 System

NoDiesel Generating UnitsGeneration Cost Function
#1Unit #1 DG Tahuna2
F1
= 110,601.086 + 3,501.193P1
+ 0.1197P1
#2Unit #4 DG Tahuna2
F2
= 39,539.790 + 4,289.739P2 - 1.1966P2
#3Unit #5 DG Tahuna2
F3
= 62,254.468 + 3,641.193P3 - 1.1966P3
#4Unit #6 DG Tahuna2
F4
= 287,095.439 + 3,427.005P4
-
0.7179P4
#5Unit #7 DG Tahuna2
F5
= 299,001.410 + 3,482.048P5 - 0.8376P5
#6Unit #8 DG Tahuna2
F6
= 240,680.101 + 3,332.475P6 - 0.3590P6
#7Unit #9 DG Tahuna2
F7
= 78,134.281 + 3,502.390P7 - 0.3590P7
#8Unit #10 DG Tahuna
(Ex UIW MMU)
2
F8
= 425,743.164 + 4,087.517P8 - 1.1966P8
#9Unit #11 DG Tahuna
(Ex PLN UIW MMU)
2
F9
= 411,527.794 + 4,105.466P9 - 1.1966P9
#10Unit #12 DG Tahuna
(Ex PLN
UIW MMU)
= 186,953.659 + 3,707.005P10 - 0.5983P102
F10
#11Unit #1 DG Tahuna (KIP 1)= 228,104.045 + 3,303.757P11 - 0.7179P112
F11
#12Unit #1 DG Tahuna (KIP 2)= 310,117.639 + 3,432.988P12 - 0.3590P122
F12
#13Unit #1 DG Tamako0.0957P132
F13
= 6,358.626 + 2,788.031P13
-
#14Unit #2 DG Tamako= 289,991.163 + 4,439.312P14 - 2.3932P142
F14
#15Unit #4 DG Tamako0.8376P152
F15
= 201,288.688 + 3,605.296P15
-
#16Unit #5 DG Tamako
(Genset Mobile)
= 92,091.190 + 2,924.442P16 - 1.1966P162
F16
#17Unit #6 DG Tamako
(Genset Mobile)
= 35,273.982 + 3,954.697P17 - 5.9829P172
F17
#18Unit #7 DG Tamako
(Genset Mobile)
= 33,700.479 + 3,509.569P18 - 5.9829P182
F18

Based on Table 2, the Fi variable represents the total generation cost of each diesel generating unit (Rp), while the Pi variable represents the power generated of each diesel generating unit (kW). Apart from the Tahuna 1 System, the cost function of the diesel generating unit can be varied according to data on other isolated systems in Indonesia to be further calculated using the same research method. So, it can be calculated further to minimize generation costs according to the needs of each system as a continuous improvement.

2) Objective Function and Constraint Function Modeling

Optimization of the generating unit operation scheduling scheme aims to minimize the LCOE. Based on the generation cost function for each generating unit from Table 2, the objective function can be formulated through (2).

\[\min F(P_{i,t}, U_{i,t}) = \sum_{i=1}^{N} \sum_{t=1}^{T} [F_i(P_{i,t}) U_{i,t}]\] (2)

where:

min F(Pi,t) = objective function

Fi = cost function of generating unit-i

Pi,t = generated power of generating unit-i at hour-t (kW)

Ui,t = on/off status of generating unit-i at hour-t

There are several constraints that must be met in designing the diesel generating unit operation scheduling scheme in the Tahuna 1 System which are divided into system constraints and diesel generating unit constraints. When viewed from a system perspective, the constraints used in this research are shown in (3), (4), and (5) below.

\[P_{load,t} = P_{PV,t} + P_{MH,t} + \sum_{i} P_{DG,t}\] (3)

where:

Pload,t = system load at hour-t (kW)

PPV,t = output power of solar power plant at hour-t (kW)

PMH,t = output power of micro hydro power plant at hour-t (kW) PDG,t = output power of diesel power plant at hour-t (kW)

Equation (3) is used as a constraint of the total load of the Tahuna 1 System at hour-t for 24 hours. The system load at hour-t is the accumulation of the output power of the Tahuna Diesel Power Plant, Tamako Diesel Power Plant, Sangihe Solar Power Plant, and Ulung Peliang Micro Hydro Power Plant. Meanwhile, other system constraints are explained in (4) and (5) as follows.

\[P_{reserve\ DG,t} \ge P_{PV,t} + P_{MH,t} \tag{4}\]

\[\sum_{i=1}^{N} (P_{i \max DG} - P_{i,t DG}) U_{i,t DG} \ge P_{PV,t} + P_{MH,t}\] (5)

where:

Preserve DG,t = spinning reserve of diesel generating unit at hour-t (kW) Pi max DG = maximum operating power of diesel generating unit (kW) Pi,t DG = generated power of diesel generating unit at hour-t (kW) Ui,t DG = on/off status of diesel generating unit-i at hour-t

Equation (4) is the constraint of the total spinning reserves of the diesel generating unit which must be more than or equal to the accumulated power generated by the solar and micro hydro power plant at hour-t. Spinning reserves of diesel generating unit are needed as a form of mitigation if the power generated from solar and micro hydro power plant decreases or is lost due to the dependence of renewable energy generators on natural conditions. Meanwhile, (5) explained that spinning reserves of diesel generating unit are the accumulation of the difference in its maximum capacity to the power generated at hour-t. Another constraint that must be met in this research is the net power capacity of the diesel generating unit as shown in Table 2, which is modeled into the (6).

\[P_{i\,min} \le P_{i,t} \le P_{i\,max} \tag{6}\]

where:

Pi min = minimum output power of generating unit (kW) Pi,t = generated power of generating unit at hour-t (kW) Pi max = maximum output power of generating unit (kW)

3) Priority Order of Generating Unit Operations

Dynamic changes in the priority list of diesel generating unit operations along with the increase in power generated due to the characteristics of non-linear generation costs per kWh. In this research, the unit Rp/kWh is used in the generation cost ramping of diesel generating unit to make it easier to determine the generation cost ranking when producing the same power. The results of the roots of the equation of the generation cost function per kWh of diesel generating units that intersect are used as the limits of the generation zone range, so that variations in operation combinations can be obtained based on range of the generation zone. In Figure 4 below is the intersection of 18 graphs of the generation cost function per kWh of the diesel generating unit.

3

Figure 4. Intersection of the generation cost function per kWh

Figure 4 explains that the graphs of the generation cost function per kWh of each diesel generating unit intersect each other as the power generated increases. From the two results of the roots of the generation cost function equation per kWh obtained, the root with a positive value is selected to be used as the boundary of the generation zone range. In each generation zone, there are variations in the combination of the priority order of diesel generating unit operations according to the amount of power generated, along with calculations of generation costs per kWh. The list of priority order of operations in each generation zone is then used as a reference in preparing the schedule of the generating unit operation schedule per hour in the Tahuna 1 System to obtain a more optimal LCOE value compared to the current condition.

2. Result and Analysis

A. Scheduling Scheme with Consideration of Spinning Reserves on One Operating Unit

After obtaining a list of priority order of the diesel generating unit operations in the Tahuna 1 System according to their respective generation zone boundaries, an hourly operation scheduling scheme can be prepared to fulfill the system's power needs. As a comparison of research results, the scheduling of diesel generating unit operations will be divided into two schemes, namely consideration of spinning reserve on one and two operating units. As explained in (3) which is one of the constraint functions, the total load of the Tahuna 1 System at hour-t is the accumulation of the output power of the solar power plant, micro hydro power plant, and diesel power plant. The hourly loading value of the diesel generating unit can be obtained which will be optimized through the operation scheduling scheme as in Table 3 below [3].

Table 3. Load profile per hour of the Tahuna 1 System 2023 [3]

Load (kW)
TimeTotalSolar + Micro HydroDiesel
System(Control Variables)(Independent Variables)
01:004,2131474,066
02:004,0811373,944
03:004,0161343,882
04:003,9941293,865
05:004,0831323,951
06:004,1501523,998
07:004,0862893,797
08:004,2064333,773
09:004,3685273,841
10:004,4685823,886
11:004,5795973,982
12:004,5735743,999
13:004,6205594,061
14:004,6425184,124
15:004,5944474,147
16:004,7143154,399
17:004,6861784,508
18:005,4471855,262
19:005,7102185,492
20:005,5022065,296
21:005,2491955,054
22:004,9091754,734
23:004,6051664,439
00:004,3791584,221

Based on Table 3, the hourly loading value on the diesel generating unit is obtained which will be used to schedule operations as independent variables. Based on the ramping data of the cheapest to most expensive LCOE values of the generating unit and the generation zone range, as well as the constraint function of the generating unit's net power capacity from Table 1, the operation scheduling scheme for the diesel generating unit of the Tahuna 1 System can be designed with consideration of the spinning reserves on one operating unit in the Table 4 below.

Table 4. Operation scheduling scheme with one spinning reserve unit in Tahuna 1 System

Operating Generating Units
TimeIIIIIIIVVVIVIIVIII
01:00#11#13#12#17#16OFFOFFOFF
02:00#11#13#12#18OFFOFFOFFOFF
03:00#11#13#12#18OFFOFFOFFOFF
04:00#11#13#12OFFOFFOFFOFFOFF
05:00#11#13#12#18OFFOFFOFFOFF
06:00#11#13#12#17OFFOFFOFFOFF
07:00#11#13#12#18OFFOFFOFFOFF
08:00#11#13#12#17#16OFFOFFOFF
09:00#11#13#12#17#16#18OFFOFF

Table 4. Continued

Operating Generating Units
TimeIIIIIIIVVVIVIIVIII
10:00#11#13#12#17#16#18#14OFF
11:00#11#13#12#17#16#18#3OFF
12:00#11#13#12#17#16#18#3OFF
13:00#11#13#12#17#16#18#3OFF
14:00#11#13#12#17#16#18#3OFF
15:00#11#13#12#17#16#18#3OFF
16:00#11#13#12#17#16#18#3OFF
17:00#11#13#12#17#16#18#3OFF
18:00#11#13#12#15#17#18#14OFF
19:00#11#13#12#15#17#16#18#3
20:00#11#13#12#15#17#16OFFOFF
21:00#11#13#12#15#18OFFOFFOFF
22:00#11#13#12#5OFFOFFOFFOFF
23:00#11#13#12#17#16#18#3OFF
00:00#11#13#12#17#16#18OFFOFF

Based on Table 4, it is known that Unit #11 and Unit #13 are always the priority generating units to be operated in the hourly scheduling scheme with consideration of the spinning reserve on one operating unit. This is because both units have the lowest LCOE value sequentially within their generation zone range. In this scheme, Unit #11 functions as the spinning reserve because it is always the priority of operation and fulfills the capacity constraint function as (6). Meanwhile, the loading value of each generating unit operating per hour is shown in Table 5.

Table 5. Loading value of the operating generating units with one spinning reserve

Loading Value of Operating Generating Units (kW)
TimeTotal
IIIIIIIVVVIVIIVIII(kW)
01:009536302,25018053OFFOFFOFF4,066
02:009746302,25090OFFOFFOFFOFF3,944
03:009666302,25036OFFOFFOFFOFF3,882
04:009856302,250OFFOFFOFFOFFOFF3,865
05:009816302,25090OFFOFFOFFOFF3,951
06:009486302,250170OFFOFFOFFOFF3,998
07:008276302,25090OFFOFFOFFOFF3,797
08:006676302,25018046OFFOFFOFF3,773
09:005736302,25018018028OFFOFF3,841
10:005186302,2501801809038OFF3,886
11:005036302,25018018090149OFF3,982
12:005266302,25018018090143OFF3,999
13:005416302,25018018090190OFF4,061
14:005826302,25018018090212OFF4,124
15:006536302,25018018090164OFF4,147
16:007856302,25018018090284OFF4,399
17:009226302,25018018090256OFF4,508

Table 5. Continued

TimeLoadingValue of Operating Generating Units (kW)Total
1 IIIICIIIIIIIVVVIVIIVIII(kW)
18:009156302,25011701809027OFF5,262
19:008826302,2501170180180901105,492
20:008946302,2501170180172OFFOFF5,296
21:009146302,250117090OFFOFFOFF5,054
22:009256302,250929OFFOFFOFFOFF4,734
23:009346302,25018018090175OFF4,439
00:009426302,25018018039OFFOFF4,221

From Table 5, it is known that the operation scheduling scheme with consideration of the spinning reserve on Unit #11 can be met. This can be proven by the total scheduled hourly load (kW) in accordance with the loading value of the diesel generating unit of Tahuna 1 System as an independent variable as shown in Table 3 before. After determining the loading value per hour, the generation cost per hour for each diesel generating unit can be calculated in Table 6.

Table 6. Generation cost of the operating generating units with one spinning reserve

TimeTotal CostLCOE
TimeIIIIIIIVVVIVIIVIII(Rp)(Rp/kWh)
01:002,724,5381,724,8256,217,035553,273243,725OFFOFFOFF11,463,3962,819.33
02:002,764,8641,724,8256,217,035301,100OFFOFFOFFOFF11,007,8232,791.03
03:002,749,5761,724,8256,217,035152,291OFFOFFOFFOFF10,843,7272,793.34
04:002,785,7341,724,8256,217,035OFFOFFOFFOFFOFF10,727,5932,775.57
05:002,778,1651,724,8256,217,035301,100OFFOFFOFFOFF11,021,1252,789.45
06:002,714,8431,724,8256,217,035534,667OFFOFFOFFOFF11,191,3692,799.24
07:002,469,2861,724,8256,217,035301,100OFFOFFOFFOFF10,712,2462,821.24
08:002,112,3031,724,8256,217,035553,273224,084OFFOFFOFF10,831,5192,870.80
09:001,885,4341,724,8256,217,035553,273579,721127,278OFFOFF11,087,5662,886.64
10:001,746,8081,724,8256,217,035553,273579,721301,100455,229OFF11,577,9922,979.41
11:001,708,2471,724,8256,217,035553,273579,721301,100578,227OFF11,662,4282,928.79
12:001,767,2411,724,8256,217,035553,273579,721301,100558,476OFF11,701,6722,926.15
13:001,805,3071,724,8256,217,035553,273579,721301,100710,885OFF11,892,1462,928.38
14:001,907,7051,724,8256,217,035553,273579,721301,100780,408OFF12,064,0672,925.33
15:002,079,3181,724,8256,217,035553,273579,721301,100627,227OFF12,082,5002,913.55
16:002,379,1361,724,8256,217,035553,273579,721301,100999,842OFF12,754,9332,899.51
17:002,663,8521,724,8256,217,035553,273579,721301,100915,981OFF12,955,7882,873.95
18:002,649,9581,724,8256,217,0353,272,886553,273301,100408,108OFF15,127,1852,874.80
19:002,583,5091,724,8256,217,0353,272,886553,273579,721301,100448,30715,680,6562,855.18
20:002,607,8531,724,8256,217,0353,272,886553,273559,696OFFOFF14,935,5672,820.16
21:002,647,9671,724,8256,217,0353,272,886301,100OFFOFFOFF14,163,8132,802.50
22:002,669,7851,724,8256,217,0352,810,935OFFOFFOFFOFF13,422,5802,835.36
23:002,687,5071,724,8256,217,035553,273579,721301,100662,818OFF12,726,2802,866.92
00:002,703,162OFF11,939,4902,828.59
293,573,4612,857.97
7Total Daily Energy Producction (kWh))102,7212,031.71

Table 6 is the result of the calculation of the generation cost of the diesel generating units of Tahuna 1 System per hour considering the spinning reserve on one operating unit. In one day, the total generation cost required from all diesel generating units in the Tahuna 1 System is Rp293,573,461 to produce 102,721 kWh of electrical energy. Based on these two calculation results, the daily LCOE value of the diesel generating unit can also be obtained, namely Rp2,857.97/kWh.

B. Scheduling Scheme with Consideration of Spinning Reserves on Two Operating Units

After the scheduling scheme for the diesel generating units of Tahuna 1 System has been prepared with consideration of the spinning reserves on one operating unit, the next step is to prepare a scheduling scheme for the operation of the diesel generating units according to the hourly system loading requirements with consideration of the spinning reserves on two operating units as in Table 7 below.

Table 7. Operation scheduling scheme with two spinning reserve units in Tahuna 1 System

Operating Generating Units
TimeIIIIIIIVVVIVIIVIII
01:00#11#13#12#17#16#18OFFOFF
02:00#11#13#12#17#16OFFOFFOFF
03:00#11#13#12#17OFFOFFOFFOFF
04:00#11#13#12#18#17OFFOFFOFF
05:00#11#13#12#17#16OFFOFFOFF
06:00#11#13#12#17#16OFFOFFOFF
07:00#11#13#12#17#16#18OFFOFF
08:00#11#13#12#17#16#18#3OFF
09:00#11#13#12#5OFFOFFOFFOFF
10:00#11#13#12#15OFFOFFOFFOFF
11:00#11#13#12#15#18OFFOFFOFF
12:00#11#13#12#15OFFOFFOFFOFF
13:00#11#13#12#15#18OFFOFFOFF
14:00#11#13#12#15OFFOFFOFFOFF
15:00#11#13#12#15OFFOFFOFFOFF
16:00#11#13#12#15OFFOFFOFFOFF
17:00#11#13#12#17#16#18#3#2
18:00#11#13#12#15#17#16#18#14
19:00#11#13#12#15#17#16#18#3
20:00#11#13#12#15#17#16#18#3
21:00#11#13#12#15#17#18#14OFF
22:00#11#13#12#15OFFOFFOFFOFF
23:00#11#13#12#17#16#18#3#2
00:00#11#13#12#17#16#18#3OFF

Based on Table 7, it is known that in this scheduling scheme Unit #11 and Unit #12 are used as spinning reserves because they are always the priority of operations for 24 hours and fulfill the constraint function of the power capacity as (6). Meanwhile, Unit #13 cannot be operated as a spinning reserve unit because this diesel generating unit does not fulfill the constraint function. The loading value for each diesel generating unit operating per hour is shown in Table 8 as follows.

Table 8. Loading value of the operating generating units with two spinning reserves

Loading Value of the operating generating units with two spinning reserve Loading Value of Operating Generating Units (kW) Total
TimeLoadingTotal
TimeIIIIIIIVVVIVIIVIII(kW)
01:009536302,10318018020OFFOFF4,066
02:009636302,11318058OFFOFFOFF3,944
03:009666302,116170OFFOFFOFFOFF3,882
04:009716302,1219053OFFOFFOFF3,865
05:009686302,11818055OFFOFFOFF3,951
06:009486302,098180142OFFOFFOFF3,998
07:008116301,96118018035OFFOFF3,797
08:006676301,81718018090209OFF3,773
09:005736301,723915OFFOFFOFFOFF3,841
10:005186301,6681,070OFFOFFOFFOFF3,886
11:005036301,6531,17026OFFOFFOFF3,982
12:005266301,6761,167OFFOFFOFFOFF3,999
13:005416301,6911,17029OFFOFFOFF4,061
14:005876301,7371,170OFFOFFOFFOFF4,124
15:006536301,8031,061OFFOFFOFFOFF4,147
16:007856301,9351,049OFFOFFOFFOFF4,399
17:009226302,072180180903001344,508
18:009156302,0651,17018018090325,262
19:008966302,0461,170180180903005,492
20:008946302,0441,170180180901085,296
21:009056302,0551,1701809024OFF5,054
22:009256302,0751,104OFFOFFOFFOFF4,734
23:009346302,08418018090300414,439
00:009426302,09218018090107OFF4,221

It is known from Table 8, the hourly scheduling scheme with consideration of spinning reserves on two operating units, namely Unit #11 and Unit #12, can be fulfilled. This can be proven by the total scheduled hourly load (kW) in accordance with the loading value of the diesel generating unit of Tahuna 1 System as an independent variable. The generation cost per hour for each diesel generating unit is calculated in Table 9.

Table 9. Generation cost of the operating generating units with two spinning reserves

Tuble 7. Generation cost of the operating generating aims with two spinning reserves
TimeGeneration Costof OperatingGeneratingg Units (Rp)Total CostLCOE
TimeIIIIIIIVVVIVIIVIII(Rp)(Rp/kWh)
01:002,724,5381,724,8255,942,090553,273579,721101,499OFFOFF11,625,9462,859.31
02:002,743,8201,724,8255,961,285553,273257,684OFFOFFOFF11,240,8872,850.12
03:002,749,5761,724,8255,967,030534,667OFFOFFOFFOFF10,976,0972,827.43
04:002,759,1421,724,8255,976,590301,100228,067OFFOFFOFF10,989,7232,843.40
05:002,753,4071,724,8255,970,856553,273249,316OFFOFFOFF11,251,6772,847.80
06:002,714,8431,724,8255,932,465553,273483,234OFFOFFOFF11,408,6412,853.59
07:002,435,2421,724,8255,661,765553,273579,721149,206OFFOFF11,104,0332,924.42
08:002,112,3031,724,8255,362,708553,273579,721301,100770,996OFF11,404,9273,022.77
09:001,885,4341,724,8255,159,4592,783,810OFFOFFOFFOFF11,553,5283,007.95
TD 11 1\(\sim\)~1
hle.9Continuıed
TimeGeeneration CosTotal CostLCOE
TillicIIIIIIIVVVIVIIVIII(Rp)(Rp/kWh)
10:001,746,8081,724,8255,037,5963,099,980OFFOFFOFFOFF11,609,2082,987.44
11:001,708,2471,724,8255,003,9833,272,886120,905OFFOFFOFF11,830,8452,971.08
12:001,767,2411,724,8255,055,4563,267,942OFFOFFOFFOFF11,815,4642,954.60
13:001,805,3071,724,8255,088,8213,272,886130,446OFFOFFOFF12,022,2852,960.42
14:001,920,0271,724,8255,190,1333,272,886OFFOFFOFFOFF12,107,8702,935.95
15:002,079,3181,724,8255,332,8393,083,597OFFOFFOFFOFF12,220,5782,946.85
16:002,379,1361,724,8255,608,8703,061,541OFFOFFOFFOFF12,774,3722,903.93
17:002,663,8521,724,8255,882,127553,273579,721301,1001,046,920592,87913,344,6992,960.23
18:002,649,9581,724,8255,868,4923,272,886553,273579,721301,100429,59915,379,8542,922.82
19:002,611,8911,724,8255,831,3043,272,886553,273579,721301,1001,046,92015,921,9202,899.11
20:002,607,8531,724,8255,827,3753,272,886553,273579,721301,100441,54615,308,5802,890.59
21:002,629,9871,724,8255,848,9523,272,886553,273301,100395,156OFF14,726,1792,913.77
22:002,669,7851,724,8255,887,9603,160,647OFFOFFOFFOFF13,443,2182,839.72
23:002,687,5071,724,8255,905,420553,273579,721301,1001,046,920213,40813,012,1752,931.33
00:002,703,1621,724,8255,920,892553,273579,721301,100438,162OFF12,221,1362,895.32
299,293,8422,913.66
Total DailyEnergy Produuction (kWh)102,7212,913.00

Table 9 is the result of the calculation of the hourly generation costs of the generating units of Tahuna 1 System with consideration of the spinning reserves on the two operating units. The total daily generation cost required from all diesel generating units in the Tahuna 1 System is Rp299,293,842 to produce 102,721 kWh of electrical energy. Based on the calculation results, the daily LCOE value can also be obtained, namely Rp2,913.66/kWh.

C. Comparison of Levelized Cost of Electricity

The diesel generating unit operation scheduling scheme per hour has been designed in subchapter 3A and 3B in accordance with all the constraint functions that have been set. Next, a comparison graph can be made of the average daily LCOE value between the existing conditions and two new schemes as shown in Figure 5 below.

6

Figure 5. Comparison of average daily LCOE value of the Tahuna 1 System between schemes

Based on Figure 5, the LCOE value of the Tahuna 1 System experienced a significant decrease after designing two types of hourly diesel generating unit operation scheduling schemes. In the first new scheduling scheme, namely by considering one spinning reserve unit, there was a decrease in the LCOE value from the initial Rp4,210.86/kWh to Rp2,857.97/kWh or 32.13% more efficient. Meanwhile, in the diesel generating unit operation scheduling scheme by considering two spinning reserve units, the calculated LCOE value was Rp2,913.66/kWh or 30.81% more economical when compared to the average cumulative LCOE value of the Tahuna 1 System for the period 2023.

Since the average LCOE value of the Tahuna 1 System in existing condition is the result of calculations based on cumulative data of 2023, by assuming the total daily electrical energy production between all the diesel generating unit operational scheduling schemes is the same, namely 102,721 kWh, the daily generation cost can then be calculated as in Figure 6 below.

3

Figure 6. Comparison of daily generation costs of the Tahuna 1 System between schemes

Figure 6 explains that there was a decrease in the total daily generation costs in the Tahuna 1 System after designing two types of hourly diesel generating unit operation scheduling schemes. In the first new scheduling scheme, which considers one spinning reserve unit, there was a decrease in the total daily generation costs from the initial Rp432,543,816 to Rp293,573,461 in producing 102,721 kWh of electrical energy. Meanwhile, from the calculation results of the diesel generating unit operation scheduling scheme considering two spinning reserve units, there was a decrease in the total daily generation costs to Rp299,293,842.

In the diesel generating unit operation scheduling scheme with the consideration of one spinning reserve unit, only Unit #11 is used as a spinning reserve because this unit is always the main priority in the operation scheduling scheme. Based on performance data, during 2023 there were a total of 7 blackouts in the Tahuna 1 System caused by a decrease in the Sangihe Solar Power Palant load due to a sudden decrease in solar irradiation accompanied by a malfunction of the spinning reserve unit in responding to system frequency fluctuations [12]. Therefore, in the scheduling scheme with the consideration of two spinning reserve units, Unit #11 and Unit #12 will be used as spinning reserves because both engine units are the priority of operation for 24 hours and fulfill the constraint function of the capable power capacity.

With the difference in LCOE value of two new scheduling schemes are relatively small, namely Rp55.69/kWh, the scheduling scheme with consideration of two spinning reserve units is a better choice to implement when viewed from the aspect of the reliability of the operation of the 20 kV distribution system. This is because there is still a backup diesel generating unit (Unit #12) as a form of mitigation if the main spinning reserve (Unit #11) fails to respond to the solar and/or micro hydro power plant which experiences a decrease or loss of load.

D. Research Results Product

The results of this research are in the form of a file with the format (*.xlxs) from Microsoft Excel so that it can be used by all stakeholders easily without any additional plugins to be installed, starting from dispatcher of the Tahuna Command Center to operators at each power plant in determining the priority order of unit operations in real time according to the estimated hourly system load. Figure 7 below shows a display of Microsoft Excel which can be used in the operation of the Tahuna 1 System as the result of this research.

Load (kW)Unit Commitment
TimeSystemSpin
Reserve
UC DGIIIIIIIVVVIVIIVIII
01:0042131474066#11#13#12#17#16#18OFFOFF
02:0040811373944#11#13#12#17#16OFFOFFOFF
03:0040161343882#11#13#12#17OFFOFFOFFOFF
04:0039941293865#11#13#12#18#17OFFOFFOFF
05:0040831323951#11#13#12#17#16OFFOFFOFF
06:0041501523998#11#13#12#17#16OFFOFFOFF
07:0040862893797#11#13#12#17#16#18OFFOFF
08:0042064333773#11#13#12#17#16#18#3OFF
09:0043685273841#11#13#12#5OFFOFFOFFOFF
10:0044685823886#11#13#12#15OFFOFFOFFOFF
11:0045795973982#11#13#12#15#18OFFOFFOFF
12:0045735743999#11#13#12#15OFFOFFOFFOFF
13:0046205594061#11#13#12#15#18OFFOFFOFF
14:0046425184124#11#13#12#15OFFOFFOFFOFF
15:0045944474147#11#13#12#15OFFOFFOFFOFF
16:0047143154399#11#13#12#15OFFOFFOFFOFF
17:0046861784508#11#13#12#17#16#18#3#2
18:0054471855262#11#13#12#15#17#16#18#14
19:0057102185492#11#13#12#15#17#16#18#3
20:0055022065296#11#13#12#15#17#16#18#3
21:0052491955054#11#13#12#15#17#18#14OFF
22:0049091754734#11#13#12#15OFFOFFOFFOFF
23:0046051664439#11#13#12#17#16#18#3#2
00:0043791584221#11#13#12#17#16#18#3OFF

Figure 7. Display of research results product (*.xlxs)

The red box in Figure 7 represents the hourly load estimation form for the Tahuna 1 System, which can be filled in dynamically and in real time by the dispatcher of the Tahuna Command Center and power plant operator as needed. Furthermore, a formula implemented in Microsoft Excel will automatically determine the priority order for generating unit operations to achieve the lowest generation cost without compromising system reliability by providing spinning reserve units.

3. Conclusions

Based on the results of the design, simulation, and analysis that have been carried out to reduce the LCOE value in the Tahuna 1 System, the following conclusions can be obtained:

1. Two hourly diesel generating unit operation scheduling schemes for the Tahuna 1 System have been designed on 12 generating units of Tahuna Diesel Power Plant and 6 generating units of Tamako Diesel Power Plant by considering the system load profile, generating unit capacity, and spinning reserve as a constraint function.

  • 2. There is a decrease in the total daily generation costs of the Tahuna 1 System on the scheduling scheme with one spinning reserve unit from Rp432,543,816 to Rp293,573,461 to produce 102,721 kWh of electrical energy, and the decrease in the LCOE value from Rp4,210.86/kWh to Rp2,857.97/kWh or 32.13% more economical compared to current conditions. Meanwhile, the scheduling scheme with two spinning reserve units shows that there is a decrease in the total daily generation costs from Rp432,543,816 to Rp299,293,842, and the decrease in the LCOE value from Rp4,210.86/kWh to Rp2,913.66/kWh or 30.81% more efficient compared to current conditions.
  • 3. With the difference in LCOE value of two new scheduling schemes are relatively small, namely Rp55.69/kWh, the scheduling scheme with consideration of two spinning reserve units is a better choice to implement from the aspect of the reliability of the operation of the 20 kV distribution system, because there is still a backup generating unit (Unit #12) which functions as mitigation if the main spinning reserve (Unit #11) fails to respond to the solar and/or micro hydro power plant when decrease or loss of load.
  • 4. Apart from the Tahuna 1 System, the cost function of the diesel generating unit can be varied according to data on other isolated systems in Indonesia to be further calculated using the same research method. So, it can be calculated further to minimize generation costs according to the needs of each system as a continuous improvement.

4. References

  • [1] Badan Pusat Statistik, Kepulauan Sangihe Regency in Figures 2024. Tahuna: BPS-Statistics Kepulauan Sangihe Regency, 2024.
  • [2] PLN UID Suluttenggo, "Laporan Pengusahaan Pembangkitan PLN UID Suluttenggo 2023," PT PLN (Persero), Manado, Indonesia, 2024.
  • [3] PLN UP3 Tahuna, "Profil Beban Sistem Tahuna 1 Periode Tahun 2023," PT PLN (Persero), Tahuna, Indonesia, 2024.
  • [4] PLN UID Suluttenggo, "Biaya Pokok Pembangkitan Sistem Sangihe Tahun 2023," PT PLN (Persero), Manado, Indonesia, 2024.
  • [5] R. D. Ananda, N. Sartika, and L. Kamelia, "Implementation of Forward Dynamic Programming in Solving Thermal Generation Scheduling," in 2023 10th International Conference on Information Technology, Computer, and Electrical Engineering, ICITACEE 2023, Institute of Electrical and Electronics Engineers Inc., 2023, pp. 293– 297. doi: 10.1109/ICITACEE58587.2023.10276771.
  • [6] V. M. Raj and S. Chanana, "Analysis of Unit Commitment Problem Through Lagrange Relaxation and Priority Listing Method". IEEE, 2014.
  • [7] Hanafi, I. F. & Dalimi, R., "Economic Load Dispatch Optimation of Thermal Power Plant Based on Merit Order and Bat Algorithm", in ICIRD 2019: 2nd IEEE International Conference on Innovative Research and Development: 28 June - 30 June 2019, Universitas Indonesia, Depok, Indonesia. IEEE, 2019.
  • [8] PLN UPDK Minahasa, "Iradiasi Matahari PLTS Sangihe Periode Tahun 2023," PT PLN (Persero), Manado, North Sulawesi, Indonesia, 2024.
  • [9] PLN UP3 Tahuna, "SLD SISTEM SANGIHE UPDATE FEBRUARI 2024 TANPA GARDU," PT PLN (Persero), Tahuna, Indonesia, 2024.
  • [10] Lestari, T. D., Wibowo, R. S., & Aryani, N. K., "Unit Commitment Menggunakan Metode Mixed Integer Linear Programming". JURNAL TEKNIK POMITS Vol. 124, No. 1, 1-7, 2021.
  • [11] Wood, A. J., Wollenberg, B. F., & Sheble G. B., Power Generation Operation and Control, 3rd edition, 2014.
  • [12] PLN UP3 Tahuna, "Rekap Blackout Sistem Tahuna 1 Periode Tahun 2023," PT PLN (Persero), Tahuna, Indonesia, 2024.

Gilang Cahyo Nugroho earned his Bachelor's Degree in Electrical Engineering from Universitas Diponegoro, Master of Science in Mechanical Engineering from Institut Teknologi Bandung. He is currently working at PT PLN (Persero) UP2D Suluttenggo as a Team Leader of Electromechanical Maintenance.

Poetro Lebdo Sambegoro earned his Bachelor's Degree in Mechanical Engineering from Institut Teknologi Bandung, Master of Science and Doctor of Philosophy in Mechanical Engineering from Massachusetts Institute of Technology. He is an Assistant Professor of Heat Transfer at the Faculty of Mechanical and Aerospace Engineering, Institut Teknologi Bandung.