Abstract
Adequate and accurate models of thermal power plants (TPP) and hydropower plants (HPP) are of great importance for the optimal economic functioning of electric power systems. The extracted models have a direct impact on economic dispatch calculations. In this research, a modified optimization algorithm called an improved gradient-based optimizer (IGBO) is deployed for the optimal extraction of TPP and HPP input-output parameters. Firstly, the IGBO is tested for use with well-known benchmark functions and shows outstanding performance over the original GBO and other competitive algorithms. For the input-output parameters extraction of TPP and HPP, the sum of the absolute error (SAE) is utilized as a fitness function to be minimized. Secondly, nine models of TPP and HPP are employed for parameter identification using IGBO. Simulation outcomes prove the capability of IGBO to accurately extract input-output parameters of TPP and HPP. Moreover, the convergence characteristics of IGBO are remarkable among investigated optimization algorithms.
1. Introduction
The reliable, secure, and economical operation of electrical power systems is the major concern of power systems planners [1, 2]. One of the planning tools is the economic dispatch (ED) of generation units. Its economic scheduling is proportional to operating cost [3, 4]. Hence, the ED aims to determine the generation share of power plants that will keep a certain fitness function minimized [5, 6]. The fitness function can be, for example, the operation cost and emission cost [7–9]. Moreover, the overall accuracy of the ED is entirely dependent on numerous elements such as the chosen model, input data and parameters employed, and even the technique that is used for the solution of this optimization problem [10, 11]. The input-output curves of thermal power plants (TPP) or nonthermal like hydropower plants (HPP) are influenced by a variety of factors, including ambient operating temperature and generating unit ageing [12]. Because the accuracy of the estimated parameters impacts the final accuracy of the ED problem, it is critical to precisely estimate the numerical values of these coefficients.
1.1. Literature Review
A long time ago, several mathematical techniques were deployed for the parameter identification of input-output curves of TPP and HPP. Among these techniques are the least square error, Marquardt algorithm, Powell regression algorithm, least absolute value, and Kalman filters that have been investigated in [13–17]. Unfortunately, these mathematical procedures necessitate a significant amount of computational work, resulting in a long computation time. Furthermore, in the presence of significant errors, these classical estimators underachieve. Also, the recursive nature of dynamic filters like Kalman filters necessitates a huge amount of data to achieve a steady-state solution.
It is worth mentioning that the input-output curves of power plants are nonconvex and highly nonlinear, and their discontinuity may result from the valve point effect [18]. This makes the use of conventional techniques of estimation like the Newton–Raphson method [19] leaving much to be desired. The metaheuristic optimization techniques surpass the mathematical ones as per their fast convergence, robustness, and avoidance trapping in regional optimal points. Recently, many metaheuristic algorithms are deployed for the parameter extraction of input-output curves of power plants like particle swarm optimization (PSO) [20], least error square (LES) [20], artificial bee colony (ABC) [21], differential evolution (DE) [22], cuckoo search (CS) [23], radial movement optimization (RMO) [24], genetic algorithm [25], teaching learning-based optimization (TLBO) [26], crow search algorithm (CSA) [27], and whale optimization algorithm (WOA) [28].
1.2. Motivation
Despite the succinct literature representation, the No Free Lunch theorem steers us that the extraction of TPP and HPP input-output parameters is probably improved using modern optimization approaches. Hence, in this study, we proposed a gradient-based optimizer (GBO) to extract TPP and HPP input-output parameters. The GBO algorithm is considered a metaheuristic optimization technique inspired by gradient descent and the Newton method [29]. GBO has the benefit of rapid convergence due to gradient search rules. In order to decrease the probability of escaping from local optima, the improved GBO (IGBO) technique is a modified version of GBO proposed in this paper. The improved optimization algorithms show enhanced performance where they are employed, for example, solving optimum reactive power dispatch problem [30], day-ahead scheduling in microgrids [31], and unit commitment [32]. The IGBO technique is an enhancement applied to the GBO using utilization from one of the recent metaheuristic algorithms that are named Marine Predators Algorithm (MPA) [33] to enhance the balance between global search and exploration and increase the strength of the proposed IGBO technique for several high-dimensional optimization problems.
1.3. Contributions and Paper Organization
To the best of the authors’ knowledge, there is no research work on the application of GBO for the parameter extraction of input-output curves of TPP and HPP as a resource of renewable energy. Also, the outstanding performance of its improved version for the application (IGBO) to one of the nonconvex, nonlinear engineering optimization problems motivates the authors for this performance investigation.
The following are the paper’s contributions:(i)Proposing an innovative approach known as IGBO via improving the original GBO.(ii)The execution of the suggested algorithm is certified via fair comparisons among the suggested IGBO and other earlier optimizers, and the conventional GBO is based on employing 23 benchmark functions.(iii)Utilization of suggested IGBO for modelling TPP and HPP.(iv)Fair comparison among the suggested optimizer and other earlier optimizers for modelling TPP and HPP.
The rest of this manuscript can be organized as follows. The mathematical formulation of the input-output curves of both thermal and hydropower plants and the utilized fitness function are illustrated in Sections 2 and 3, respectively. Section 4 presents the idea of GBO and its modified version, while Section 5 demonstrates the simulation outcomes. The main findings from this research are summarized in Section 6.
2. Modelling Thermal and Hydropower Plants
The input-output curves of mathematical models of both TPP and HPP are introduced in this section. As per literature, there are five mathematical models which can express the smooth and nonsmooth incremental fuel cost; corresponds to the active power generated by a thermal unit as can be depicted in Figure 1.

2.1. Smooth Models of Thermal Units’ Fuel Cost
The thermal units’ fuel cost can be smoothly modelled via linear, quadratic, or cubic formulas as follows:
Model 1:
Model 2:
Model 3:where represent the generator cost parameters.
2.2. Nonsmooth Models of Thermal Units’ Fuel Cost
The nonsmoothness of the model of thermal units’ fuel cost due to valve loading can be presented via quadratic or cubic formula as follows:
Model 4:
Model 5:where is the thermal generating unit’s lowest generated active power.
2.3. Thermal Units’ Emission Cost
The relationship between the generated power of the thermal unit and the related emissions that is emitted into the atmosphere can be modelled using one of two models as given in equations (6) and (7):
Model 6:
Model 7:where represent the unit emission cost parameters.
2.4. Input-Output Models of HPP Cost
The HPP as one of the renewable resources completely differs from thermal power generation units. The input-output models of these plants are concerned with the relation between the reservoir actual head and the water discharge rate of the unit . This can be attained with one of the following models:
Model 8:where represent the Glimn‐Kirchmayer model parameters of the unit.Model 9:with represents Hamilton‐Lamont’s model parameters of unit.
3. Formulating the Objective Function
The accurate estimation of the input-output coefficients of TPP or HPP calls for the use of optimization algorithms. The input-output relations can be arranged in the following form:where represents the input-output curve model of a generation plant (thermal/hydro) and can be, for example, the fuel cost measurement vector of length N that is combined with an error vector . Also, denotes the plant parameters, which are specified asfor models 1–5,for models 6-7,for models 8-9.
The IGBO will deploy the fitness function (sum of the absolute error, SAE) in (14) that is required to be minimized to attain the accurate coefficients of the input-output curves of TPP and HPP.
4. The Proposed Optimization Algorithm
This section describes the concept of the conventional GBO algorithm. Then, the process of the proposed improved GBO technique is presented.
4.1. The Original GBO Algorithm
The conventional gradient-based optimizer (GBO) algorithm [29] combines gradient and population-based approaches, and it employs Newton’s method that requires the search direction to observe the search domain with the use of a collection of vectors and two main operators, namely, gradient search rule (GSR) and local escaping operators (LEO).
4.1.1. Initialization Process
The first stage of the GBO algorithm includes the number of populations in the D-dimensional search space. The initial position is randomly formed as follows:where is the nth vector, refer to the bounds of the solution space in each problem, and is a random number which is defined between [0, 1].
4.1.2. Gradient Search Rule (GSR) Process
The GSR is according to the gradient-based approach, and the targets of the GSR method are increasing the convergence rate and exploration development. Thus, the new location can be calculated as follows:
Equation (16) will be edited to contain the population-based search principle that is described bywhere denotes a random number with a normal distribution, are the best and worst solutions attained through the procedure of optimization, is a number in the range [0, 0.1], and denotes the variation in location at each iteration. From these equations, the gradient search rule (GSR) is defined aswhere is the randomly produced parameter, and it can be calculated as follows:where is a sine function for the transference from exploration to exploitation, , are equal to 0.2 and 1.2, respectively, denotes the present number of iterations, and refers to the maximum value of iterations. The modification between the optimal solution and a randomly selected location can be given bywhere is a random vector with N dimensions, ) represent the different integers randomly selected from [1, N], and is the size of step. The new location is updated according to the GSR from the following equation:
The direction of movement (DM) is added for improved exploitation of the neighboring area of that can be calculated as follows:
Consequently, the new location can be calculated after considering the GSR and DM from the following equation:
The GBO used another location to improve the local search using putting the best-so-far solution () instead of the location . The new location () can be calculated as follows:where
Based on the locations and the current location , the new location at the next iteration () is formulated as
4.1.3. Local Escaping Operator (LEO)
The LEO is used to improve the effectiveness of the GBO technique and to escape the local optima for finding the solution to complex problems. The LEO produces a suitable solution () by using numerous solutions that include , solutions and , two random solutions and , and a new randomly generated solution (). The solution is given aswhere is a uniformly distributed random number in the range of [−1,1], denotes a random number from a normal distribution with a mean of 0 and a standard deviation of 1, refers to the probability, and , , and are random values produced as follows:where is a random number in the range of [0, 1] and represents a number in the range of [0, 1]. The above equations are simply clarified as follows:where is a binary parameter with a value of 0 or 1. If parameter < 0.5, the value of = 1; otherwise, = 0. The solution is created as follows:where refers to a randomly produced solution, denotes a randomly selected solution, and denotes a random number in the range of [0, 1]. Equation (30) is adjusted as follows:where refers to a binary parameter with a value of 0 or 1. If < 0.5, the value of = 1; otherwise, = 0.
4.2. Improved Gradient-Based Optimizer (IGBO)
The original GBO algorithm suffers from the local minima when handling nonlinear objective functions. Also, the convergence rate and accuracy are low because of the randomness [34]. To overcome all these shortcomings and develop the solution competence as well as the speed of convergence and to improve the strength of the proposed IGBO technique for several high-dimensional optimization problems, the way is to utilize one of the recent metaheuristic algorithms that are named the Marine Predators Algorithm (MPA). This MPA is presented in [30]. The fittest solution is nominated as the best solution to construct a matrix which is called Elite. The Elite is calculated from the following equation:where is the best solution that is replicated times to construct the matrix. The modification process is divided into two phases:
Phase 1: the mathematical model of this phase is given asPhase 1: the population is moving to exploit the detected best locations. Therefore, populations have divided into two groups: the first group follows the Levy flight, and the second group follows the Brownian distribution during their movement, as described in the following:where is a vector that follows the Levy flight distribution, and it can be calculated as follows:where denotes the Levy flight exponent and u and are calculated as
The standard deviations and are explained aswhere is the Gamma function.
A flowchart of the suggested IGBO algorithm is revealed in Figure 2.

5. The Simulation Results and the Discussion
5.1. The Performance of IGBO
The proposed IGBO technique competence and ability are tested on the numerous benchmark functions, using the statistical measurement including the minimum values, mean values, median values, maximum values, and standard deviation (STD) for optimal solutions attained using the IGBO technique and the other recent algorithms. The data of the 23 benchmark functions are taken from [35]. The results achieved with the IGBO technique are compared with five recent techniques, hunter–prey optimization (HPO) [36], whale optimization algorithm (WOA) [37], supply demand-based optimization (SDO) [38], grey wolf optimizer (GWO) [39], and the original (GBO). All these algorithms were implemented for the maximum number of iterations of 200 and population size of 50 for 20 independent trails using MATLAB R2016a. All computations were executed on a Core i5-4210U CPU@ 2.40 GHz speed and 8 GB RAM. Figure 3 displays the qualitative metrics on F1, F1, F3, F5, F8, F9, F12, F15, and F18: 2D views of the functions, convergence curve, average fitness history, and search history.

Tables 1–3 tabulate the statistical analysis results of the IGBO algorithms and five optimization algorithms when used for three types of benchmark functions (unimodal functions, multimodal functions, and composite functions, resp.). The optimum values were attained using the IGBO, GBO, HPO, WOA, SDO, and GWO techniques displayed in bold. It is obviously clear that the IGBO technique reaches the optimum values for most of these benchmark functions. The convergence rates of all techniques for those functions are illustrated in Figure 4, and the Boxplots for those techniques are shown in Figure 5. From Figures 4 and 5, it is seen that the IGBO technique reached a stable point for all functions and the boxplots of the proposed IGBO technique are very narrow and stable for most functions compared to the other techniques.




Table 4 contains the obtained results using the Wilcoxon signed-rank test. The proposed IGBO algorithm is compared with five well-known algorithms in the literature, including GBO, HPO, WOA, SDO, and GWO, respectively, for the 23 benchmark functions.
5.2. Case Studies
This section contains the generated results of employing GBO and IGBO to identify parameters of nine models of TPP and HPP. Fuel costs of TPP are modelled using models 1 to 5. Emissions of TPP are modelled using models 6 and 7. Relationships of HPP are modelled using models 8 and 9. We use the top and bottom bounds for the parameters of models of TPP and HPP in [28]. The parameters of GBO and IGBO are listed in Table 5.
Fifty independent runs for all tested models are operated. All models are implemented using MATLAB-R2020b on a laptop with Intel Core i7-1065G7 CPU at 1.3 GHz (8 CPUs), ∼1.5 GHz, and 16 GB RAM.
All nine models are compared with other approaches which were employed in the literature in accordance with their results.
5.2.1. Clarification of the Tested Models
(1) Models of Fuel Costs of Thermal Units (Models 1–5). The models of fuel cost of three power systems that own thermal generating units are studied. Models 1–3 are employed for a 3‐unit thermal system (TS3). Models 4 and 5 are employed for 2‐unit (TS2) and 13‐unit thermal systems (TS13), respectively.
TS3 comprises three thermal power generating units (units 1–3) which are coal-, oil-, and gas-fired, respectively. The actual values of generated power (P) and the fuel cost (Factual) are in [40]. The relationships between P and Factual in models 1–3 are linear, quadratic, and cubic, respectively.
TS2 comprises two thermal units (units 1 and 2). The actual values of P and Factual are in [20]. Model 4 is used for two units and includes five parameters, which should be extracted.
TS13 comprises thirteen thermal units (units 1–13). The actual values of P and Factual are in [24]. Six parameters of model 5, which is used for thirteen units, should be extracted.
(2) Models of Emissions of Thermal Units (Models 6 and 7). A power system of five thermal generating units (TS5) is studied. The actual values of P and emissions () are in [26]. The relationships between and in models 6 and 7 are quadratic and quadratic plus an exponential term, respectively. The parameters to be extracted are three and five in models 6 and 7, respectively.
(3) Models of Hydropower Units (Models 8 and 9). A power system of one hydrogenerating unit (HS1) is studied. The actual values of P, effective head (H), and the rate of water discharge (Q) are in the Newton–Raphson methodology [19]. The relationships between P, H, and Q in models 8 and 9 are based on Glimn‐Kirchmayer and Hamilton‐Lamont’s models, respectively. There are six parameters in both models 8 and 9.
5.2.2. Results with Discussion
Tables 6–14 include optimized parameters of models 1–9 using IGBO, respectively. Additionally, attained SAE of models 1–9 using GBO, IGBO, and other approaches are written in Tables 6–14, respectively. Comparisons of IGBO with other employed approaches in the literature, NR [19], LES [20], PSO [20], ABC [21], DE [22], CS [23], RMO [24], TLBO [26], CSA [27], WOA [28], reveal that IGBO has the minimum SAE. The SAE of GBO are the same as those of IGBO in models 1–3, which have linear, quadratic, and cubic equations; that is, they are either linear or low nonlinear, while IGBO produces SAEs with fewer values in models 4–9 which have high nonlinearity. Therefore, we can conclude that the effectiveness of IGBO is more apparent when optimizing high nonlinear problems. Subsequently, the statistical indexes have been computed, that is, minimum (min) and maximum (max) values of SAE and standard deviation (STD). Smaller STD values assert the efficacy of IGBO in modelling TPP and HPP.
Regarding the SAE convergence curves which were obtained using IGBO in Figures 6–14, we can observe that quickness, smoothness, and absence of oscillations are consistent features of the SAE convergence curves until attaining the optimal SAE. The population and maximum iterations differ from one model to another due to the difference in nonlinearity between models. As the nonlinearity increases, the population and maximum iteration increase. In model 7, unit 1 requires population and maximum iterations more than other units since its data are very different from other units.







(a)

(b)


6. Conclusions
An innovative improvement for conventional GBO has been suggested in this study. The IGBO technique is considered an enhancement for the conventional GBO to improve the balance between exploration and exploitation and enhance the algorithm performance. Then, the suggested IGBO was tested for 23 benchmark functions, and it revealed good functioning compared to the conventional GBO and other earlier algorithms. Economic dispatch is one of the powerful electric power planning tools that is highly coupled with the accurate input-output curves of generation units. Extracting the parameters of input-output curves of TPP and HPP using IGBO has been illustrated in this research. With the aid of the sum of the absolute errors (SAE) that has been utilized as a cost function to be minimized, the IGBO has shown excellent performance in identifying the TPP and HPP curves parameters through conducted case studies. Moreover, with the comparison with GBO and previously utilized optimization algorithms, the IGBO convergence curve is rapidly converged to the lowest value.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the approval and support of this research study by the Deanship of Scientific Research at Northern Border University, Arar, Saudi Arabia (Grant no. ENGA-2022-11-1409).