Next Article in Journal
A Schelling Extended Model in Networks—Characterization of Ghettos in Washington D.C.
Next Article in Special Issue
Geometric Properties of Some Generalized Mathieu Power Series inside the Unit Disk
Previous Article in Journal
Bayesian and Non-Bayesian Analysis of Exponentiated Exponential Stress–Strength Model Based on Generalized Progressive Hybrid Censoring Process
Previous Article in Special Issue
A Fuzzy Random Boundary Value Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Amended Whale Optimization Algorithm for Optimal Bidding in Day Ahead Electricity Market

by
Kavita Jain
1,
Akash Saxena
1,
Ahmad M. Alshamrani
2,
Adel Fahad Alrasheedi
2,
Khalid Abdulaziz Alnowibet
2 and
Ali Wagdy Mohamed
3,*
1
Department of Electrical Engineering, Swami Keshvanand Institute of Technology, Management and Gramothan, Jaipur 302017, Rajasthan, India
2
Statistics and Operations Research Department, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
3
Operations Research Department, Faculty of Graduate Studies for Statistical Research, Cairo University, Giza 12613, Egypt
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(9), 456; https://doi.org/10.3390/axioms11090456
Submission received: 21 July 2022 / Revised: 23 August 2022 / Accepted: 26 August 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Mathematical Analysis and Applications III)

Abstract

:
Successful privatization in other sectors leads to a restructuring in the power sector. The same practice has been adopted in the electrical industry with a deregulated electricity market (EM). This enables competition among generating companies (Genco’s) for maximizing their profit and it plays a central role. With this aim, each Genco gives a higher bid that may result in a risk of losing the opportunity to get selected at auction. The big challenge in front of a Genco is to acquire an optimal bid and this process is known as the Optimal Bidding Strategy (OBS) of a Genco. In this manuscript, a new variant of whale optimization (WOA) termed the Amended Whale Optimization Algorithm (AWOA) is proposed, to attain the OBS of thermal Genco in an EM. Once the effectiveness of new AWOA is proved on 23 benchmark functions, it is applied to five Genco strategic bidding problems in a spot market with uniform price. The results obtained from the proposed AWOA are compared with other competitive algorithms. The results reflect that AWOA outperforms in terms of the profit and convergence rate. Simulations also indicate that the proposed AWOA can successfully be used for an OBS in the EM.

1. Introduction

After achieving successful results of privatization in several sectors, i.e., telecommunication, toll plaza, airlines, and many more, the reformation of the power industry was also started. The reformation of the power industry is termed the restructuring or deregulation of the electricity market (EM). The reason for deregulation in the EM is to restrict the monopolies of government or government authorities and provide a competitive platform for suppliers and buyers [1]. A competitive platform in the EM forces the generators to evaluate the cost in such a manner that they are in a risk-free zone. To reduce the risk of loosening the game with uncertainties of monopolistic market structure leads the EM to innovate a new structure of market termed oligopolistic market structure [2].
Due to certain limitations, i.e., large investment size, transmission constraints, transmission losses, etc., there are a limited number of buyers and sellers in an oligopolistic market. The aim of both buyer and seller in an oligopolistic market is to maximize their profit. All applicants submit their bids (for quantity (MW) and price ($/MW)) in a sealed envelope to the system operator (SO). The SO will finalize the market clearing price (MCP) after receiving the bids from all applicants (supplier and consumer). MCP is the effective price of the market that is found by the crossing point of the supply and load curve. In pay as bid market, market price is set through uniform price clearing mechanism; greater is the bid revenue greater the profit. Thus, every generating company (Genco) bids higher to attain higher profit but they jeopardy of losing the competition with this higher bid. Thus, for a Genco to catch an optimal bid in an oligopolistic market is a complex problem and is known as the strategic bidding problem of a Genco. Many researchers were published their work on strategic bidding problems due to their stochastic rather than deterministic nature. In these research articles, stochastic optimization approaches were used to solve this stochastic problem.
The three ways as given in the literature to solve strategic bidding problems in the EM are the game theory-based, dynamic based approaches and stochastic based approaches. Game theory approaches assumes that rival GENCO’s cost functions and complete bid information are public. This is practically not true. Additionally, multiple Nash equilibriums required for large number of players. Dynamic optimization techniques such as Lagrange Relaxation, Dynamic Programming, etc. These techniques fail for realistic non-differentiable, multi constraint and multi-objective problems and require nonlinear simplification, if adopted. Stochastic based approach gives accurate results, fast convergence, global optimum solution and reliable solution tools in an EM [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22].
In [3], the authors consider a methodology called the fuzzy adaptive particle swarm optimization for a thermal generator in a uniform price spot market, taking into account a precise model of nonlinear operating cost function and unit commitment minimum up/down limitations. The normal PDF is used to model the bidding behavior of other competing Genco’s. In [4], researchers offer particle swarm optimization (PSO) algorithms for determining market price and volumes in a competitive power market in this work. To locate solutions, the first approach combines a traditional PSO algorithm. The second approach combines the PSO strategy with a decomposition technique. This new decomposition-based PSO outperforms the traditional PSO significantly. In this research [5], a new agent-based simulation model based on the Ant Colony Optimization (ACO) algorithm is developed to compare three different wholesale electricity markets clearing strategies, namely uniform, pay-as-bid, and extended Vickery rules. In this study [6], a unique computational intelligence technique for solving the Nash optimization issue is presented. This novel process is based on the PSO algorithm, which employs the SA method to prevent particles from becoming caught in local minima or maxima and improve particle velocity functions. Other computer intelligence techniques such as PSO, Genetic Algorithm (GA), and a mathematical method (GAMS/DICOPT) are compared to the results of this operation. The IEEE 39-bus test system is used to demonstrate and validate the suggested technique’s outcomes.
In order to optimize its own profit as a market participant, the article [7] provides a new approach for bidding strategy in a day-ahead market from the perspective of a generating business (GENCO). The fuzzy adaptive gravitational search algorithm (FAGSA) is used in a unique stochastic optimization approach to tackle the optimal bidding strategy problem in a pool based power market [8]. In this research [9], a unique algorithm based on the Shuffled Frog Leaping Algorithm is used to address the optimal bidding strategy problem (SFLA). It is a memetic meta-heuristic that does a heuristic search to find a global optimal solution. It combines the advantages of the Memetic Algorithm (MA) based on genetics with the Particle Swarm Optimization (PSO) based on social behavior. As a result, it has a more precise search, which prevents premature convergence and operator selection. As a result, the suggested method overcomes the limitations of the Genetic Algorithm (GA) and the PSO method in terms of operator selection and premature convergence. In this study [10], a new strategy for developing optimal double-sided bidding strategies in security-constrained power sector is described, with pollution emission as a secondary goal. Both Generation Companies (Gencos) and Distribution Companies (DisCos) in the suggested algorithm seek to maximize their profit by implementing optimal strategies, despite the fact that they have imperfect knowledge about the rivals and the market mechanism of payment is locational marginal pricing. The optimal bidding strategy is developed using a hybrid technique based on information gap decision theory (IGDT) and modified particle swarm optimization (MPSO) in this work [11].
In order to handle the profit maximizing process in a continuously changing market, a novel form of Grey Wolf Optimizer (GWO) called the Intelligent Grey Wolf Optimizer (IGWO) is developed by the authors [12]. The Krill Herd algorithm (KHA) is used to develop an optimal bidding strategy in this article [13]. Supplier and buyer bidding coefficients are carefully chosen. The proposed KHA’s code was written in MATLAB. It was put through its paces on an IEEE 30 bus power system. The Invasive Weed Optimization technique was used to solve the Optimal Bidding Strategy problem in this paper [14]. The utilities compete in order to maximize their profits. The proposed technique was written in MATLAB and tested using the IEEE 30 bus standard. Authors [15] present an alternate methodology for determining revenue-maximizing strategic bids when the opponents’ bidding strategy is uncertain. To achieve an optimal bidding strategy in the electricity market, a hybrid architecture combining metaheuristic and supervised learning is proposed in this study [16]. The Salp Swarm Algorithm (SSA) is combined with a neural network in the suggested architecture (NN). The suggested architecture is compared to the results of SSA and Opposition-based SSA on the IEEE-14 bus system, IEEE-30 bus system, and 75-bus Indian Practical System (OSSA). To ensure increased effectiveness, a selected learning approach for strategic bidding is presented in the paper [17]. The suggested system uses an ensemble technique, in which many machine learning algorithms are used to predict the price and give a bidding recommendation. The most appropriate ones will be chosen to dominate the bidding approach as the clearing iteration advances.
In this work [18], author presents a model of neural network using Harris Hawk optimizer (HHO) to solve the problem of optimal bidding in the EM. The issue that arises when a group of small prosumers participate in the energy market is discussed in this study [19]. The aggregator takes advantage of the appliances’ flexibility to lower market net costs. There are two optimization techniques suggested. In order to reduce the cost of acquiring energy, how can a time-shiftable load, which may itself be made up of a number of smaller time-shiftable subloads, submit its demand bids to the day-ahead and real-time markets? This is the topic that this study [20] aims to address. In this study [21], authors tackle the issue of competitive bidding for a big price-maker regulatory resource in performance-based regulation markets. In order to assist an aggregator of prosumers in defining bids for the day-ahead energy and secondary reserve markets, this study [22] offers a two-stage stochastic optimization model.
In this research work, a new variant of Whale Optimization Algorithm (WOA) [18], named AWOA, is proposed for solving the optimal bidding problem in a day-ahead EM. WOA is a recently developed meta-heuristic algorithm [23] and as seen from the past research papers that, WOA performs very well on real life applications [24,25,26,27,28]. The WOA is a revolutionary nature-inspired meta-heuristic optimization algorithm that replicates the social behavior of humpback whales. The bubble-net hunting methodology motivated the algorithm. Twenty-nine mathematical optimization problems and six structural design challenges are used to evaluate WOA in [23]. To estimate short-term wind power, a hybrid forecasting model based on Complementary Ensemble Empirical Mode Decomposition (CEEMD) and Whale Optimization Algorithm (WOA)-Kernel Extreme Learning Machine (KELM) is designed to deal with the intermittent and fluctuating characteristics of wind power time series signals [24]. CEEMD first reduced the non-stationary wind power time series into a number of generally stationary components. In this study [25], authors suggest a hybrid model that is an evolution of the hGADE algorithm for addressing the Unit Commitment Scheduling Problem, a mixed-integer optimization problem. For the computation of the overall operation cost of power system operation, the Whale Optimization Algorithm was used. In [26], simulations are run on a test smart grid with loads varying in two service zones, one for residential consumers and the other for commercial users. By comparing the findings with spider monkey optimization and biogeography-based optimization, WOA demonstrates its usefulness. Simulation results show that the proposed demand side management solutions save money while lowering the smart grid’s peak load demand. The electric power system is examined in two phases in this paper [27]: decentralized and centralized ways to reduce operational costs. The Tuned Whale Optimization Algorithm (TWOA), a new artificial intelligence technique, is used to solve these phases. The IEEE 48-bus power system is used to accomplish these concepts. The IEEE 48-bus system is made up of two zones connected by transmission lines. Variations in TWOA’s regulating factors are also discussed.
The whale optimization technique for loss minimization employing FACTS devices in the transmission system is reported in this article [28]. This investigation will use a thyristor controlled series compensator (TCSC). In this research, WOA is used to determine the appropriate FACTS device size for power system loss minimization. To verify the effectiveness of the suggested technique, an IEEE 30-bus RTS was employed as the test system. Opposition-based leaning [29,30,31,32] and CM operator [33,34,35,36] are fused with WOA and experimented with over 23 benchmark functions (unimodal, mutimodal, and fixed-dimension multimodal).
The benefits of OEL and Cauchy mutation that are listed below have encouraged authors to use these in WOA in light of this research review. These qualities are listed below:
  • The curse of dimensionality problem can be solved with the help of the OEL paradigm. Due to the issue of formulating strategic bidding, the huge search space and stochastic nature of the variables make the curse of dimensionality inevitable (rival bids). The job of locating a global optimum in dynamic simulations might thus be challenging. A unique solution to this issue is provided by OEL, which also offers a way out of the neighborhood minima trap. By generating opposing points in the search space, OEL improve any algorithm’s exploration capabilities;
  • The characteristics of probable candidates that can address the strategic bidding dilemma should assist them in avoiding premature convergence. Due to the inclusion of stochastic variables throughout the simulation process, the issue of premature convergence in the strategic biding problem is significant. The introduction of OEL improves convergence speed while also guarding against premature convergence of the solver;
  • By boosting the exploratory power of whales with Cauchy distribution, Cauchy operator aids in preventing the stagnation in local optimums. Thus, the greedy selection maintains a healthy balance between the current and prior placements of whales while the Cauchy operator aids in enhancing the capacity of whales in terms of exploration and exploitation.
The WOA modified version of WOA is then applied for the bidding problem in dynamic EM. To solve this problem there are two critical parameters such as convergence of the algorithm and clarification superiority. Convergence problems are intently related to the fitness value and computation time which may be stimulating revenue and workout market circumstances. To train the rival’s performance four probability distributions are used namely: Normal, Lognormal, Gamma, and Weibull PDF [37] that is built from past market data analysis. Framing a bidding method with incomplete information about rival behavior is a huge assignment for the design engineer. Monte Carlo (MC) simulations [38] are considered powerful gear as they can be hired as sampling, optimization, and assessment gear. For the optimization of strategic bidding trouble, those MC simulations are employed wherein the goal feature is deterministic and randomness is brought artificially to greater emerald search [39]. The above-discussed application and unique feature of MC simulations inspire authors to employ MC simulations in strategic bidding problems. Slow convergence and being stuck in local optima are issues with WOA. The amended WOA method is a novel, nature-inspired heuristic technique that is proposed in this research as a means of overcoming these shortcomings when solving the strategic bidding dilemma. This strategy serves as an alternative to other current, recent algorithms. After having a bird’s eye view of the literature about the strategic bidding problem following research objectives are outlined for this manuscript:
(a)
To test the Amended Whale Optimization Algorithm (AWOA) on benchmark functions and resolve issue of bidding in the EM by confirming enhanced exploration and an optimum exploitation of search space through Oppositional Enabled Learning (OEL) and Cauchy Mutation (CM) Operator;
(b)
To search the performance of this novel variant i.e., AWOA with parent WOA (Whale Optimization Algorithm) and some recently developed algorithms are applied on benchmark functions;
(c)
To achieve different statistical assessments with Wilcoxon rank sum, box plot analysis and convergence investigation for the testing the effectiveness of the developed model;
(d)
To construct rival bidding prices using (normal, lognormal, gamma and Weibull) PDF, interpret them using the MC approach, and design a bidding strategy for the day-ahead market by entrancingly taking into account all inter-temporal limits;
(e)
To represent a fair evaluation between the outcomes acquired through optimization procedure based on profit, MCP (Market Clearing Price) calculation and solution quality.
The rest of the article is prepared as shown in Figure 1.

2. Problem Statement

To decipher the issue of strategic bidding in this article, we assumed m+1 Generating Companies (Genco’s). Genco-k is the generator whose profit has to be maximized by finding the optimal bids with m competitors in the energy market. In the EM, all the m+1 Genco’s and consumers submit their bids in a sealed envelope (contains quantity (MW) and price ($\MW)) to the system operator (SO). After the last date of submission of bids from supplier and consumer, the SO arranges the supplier’s bids in increasing order, and consumer’s bids in decreasing order where X-axis represents quantity (MW) and Y-axis shows the price($\MW). The point where both curves intersect each other is the called the equilibrium point and after drawing a horizontal line from that point to Y-axis, market clearing price (MCP) is received.
A Genco can bid for multiple (maximum l) blocks. It is preferable to bid on multiple blocks rather than just one at a time to ensure financial security. The rival’s bidding block size is supposed to be acknowledged from the past data and their bidding prices are deliberate by Probability Distribution Factor (PDF) through probability statistical analysis of historical bidding data. In this paper, we generate the input data using four different types of PDFs [32]. Thus, we create four cases to solve the problem.
  • Case I. The Normal distribution;
  • Case II. The Lognormal distribution;
  • Case III. The Gamma distribution;
  • Case IV. The Weibull distribution.
Before a unit is committed/not-committed, there is already a least-defined time subsequent of which it can be not-committed/committed again. Inter-temporal operating limits of Genco-K, such as minimum/up and minimum/downtimes have been regarded in the effort. Considering non-differentiable, non-convex cost function, nonlinear (exponential) start/up cost function, and constant shut/down cost, the operating cost function for the lth block of Genco-K is expressed as:
c l ( t ) = c l ( t ) p r + c l s u { u l ( t ) ( 1 u l ( t 1 ) ) } + c l s d { ( 1 u l ( t ) ) u l ( t 1 ) }
where
c l ( t ) p r = c 0 ( Q l ( t ) ) 2 + c 1 ( Q l ( t ) ) + c 2 + | c 4 sin ( c 4 Q min Q l ( t ) ) ) |
c l s u = h c s u + c c s d ( 1 exp ( T o f f T ) )
Due to the consecutive opening of a large number of valves to obtain ever-increasing output of the unit, input-output characteristics for large thermal generators are not always smooth [40]. A rippling effect on the unit curve is common as each steam admission valve in a turbine begins to open. [41] Approximated the rippling effect of valve point loading as a periodic rectified sinusoidal function. [42] Established the effect of valve point loading on economic dispatched output of units, confirming the need of applying a precise production cost function in strategic bidding. Equation (2) signifies this sinusoidal nonlinear characteristic, in which c0, c1, and c2 are cost coefficients and c3 and c4 are the coefficients of the valve point loading impact. An exponential characteristic is considered in Equation (3), to signify the association between the start-up value and the shut-down time. Although, any present start/up cost characteristic can be used. The value received in begin-up and begin/down is united inside the running cost characteristic so that the actual advantage is taken into consideration and the Genco’s are committed/not-committed accordingly. The OBS for Genco-k can be achieved with profit maximization in the terms of output power dispatched ( Q l ( t ) ) and MCP ( M C P ( t ) ). The product of Q l ( t ) and M C P ( t ) is defined as revenue acquired. The increasing profit of lth blocks of the Genco-k overtime period “T” is uttered as:
max P L ( t ) f ( M C P ( t ) , Q l ( t ) ) = t = 1 T l = 1 L ( M C P ( t ) Q l ( t ) c l ( t ) )
Subject to constraints
  • Generation limits
Q min u l ( t ) Q l ( t ) Q max u l ( t ) , t T .
2.
Minimum uptime
( 1 u l ( t + 1 ) ) M l u t H l ( t ) o n ,   i f     u l ( t ) = 1 .
3.
Minimum downtime
u l ( t + 1 ) M l d t H l ( t ) o f f ,   i f     u l ( t ) = 0 .
4.
Limitations on the bid price
c l ( t ) p l ( t ) p , ¯     t T .
The optimization limits described in (4)–(8) may be explained to acquire the finest block bid price of the lth block of Genco-K at hour t, signified as pl(t). In (4), pl(t) and P l m   do now not clearly emerge but those are indirectly involved within the procedure of decisive MCP. Equations (5)–(7) are the operating constraints, while (8) looks like it represents the bid price limit, i.e., p ¯ . With the usage of the PDF’s defined above competitor’s bidding prices may be received from past bidding records. Formation of the Optimal Bidding Strategy (OBS) for Genco-k, with objective function (4) and constraints (5)–(8), be remodeled a stochastic optimization concern, to be resolved by MC founded WOA. The data taken in Table 1 is used to generate the competitor’s bidding data for all the cases explained above. We generate 1000 competitor’s bidding samples in 3 blocks as shown in Figure 2.
In following section, research methodology is presented.

3. Research Methodology

The formulated problem in Section 2 has been solved with the Monte Carlo (MC) approach; which is a method to attain a probabilistic estimate of a mathematical hassle through the usage of statistical samples. This approach runs stochastic simulation with random numbers and accordingly computes the equation to determine the result. This MC simulation is integrated with, a modified version of the WOA, to attain the OBS for a Genco.

3.1. Monte Carlo (MC) Approach

MC methodology proceeds as follows:
  • For all the competitors participating in the EM, create a huge number of random test trials of block bid values considering the probability distribution equation and constraints of the PDFs;
  • With those illustrations of block bid expenses of all of the competitors, decide the huge variety of trial outcomes;
  • The average of all of the trial results gives the anticipated cost.
Details of the procedure:
  • Stipulate the number of M simulations allowed, M;
  • Set simulation counter m = 0;
  • Create random values of bid prices for each lth block of m rivals using Normal, Gamma, Weibull and Lognormal distribution function;
  • Formulate WOA to pursue the finest bid price for every lth block of Genco-k and record the optimal value;
  • Keep posted m = m + 1;
  • If m < M then is present at (step 3), else go to see (step 7);
  • Determine the predictive assessment of optimum bid price i.e., the average of (m = 1, 2… M). This final price is known as the finest bidding price (pl(t)) for the lth block of Genco-k at hour t.

3.2. Whale Optimization Algorithm (WOA)

WOA is a novel metaheuristic algorithm developed by Mirjalili and Lewis in the year 2016. The inspiration of WOA is the unique hunting performance of humpback whales. Whales are taken as smart mammals with emotion. Humpback whales are known for hunting zooplankton or small fish that live close to the sea’s surface. The bubbles net feeding methodology of hunting employed by such humpbacks. Whales do this by swimming around the target and blowing bubbles in a circular or a nine-shaped pattern, as depicted in Figure 3.
The mathematical modeling to perform optimization is given as follows:

3.2.1. Enclosing Prey

Humpback whale searches the residence of the target and encircles them. WOA adopts that the present top applicant result is the goal target or is near to the best. The other pursuit representatives attempt to appraise their locations to the best exploration applicant. The performance modeled is as:
D = C . X ( t ) X ( t )
X ( t + 1 ) = X ( t ) A . D  
where current iteration is i, coefficient vectors are A and C , position vector of the optimal result obtained is X ( i ) and position vector is X ( i ) . || is the absolute value and . is the element by element growth.
It is worth noting here that, if there is a better solution, X* should be modified in each iteration. A and C vectors are designed as follows:
A = 2 a . r a
C = 2 r
where in both exploration and exploitation stages a is linearly reduced from 2 to 0 over the sequence of iterations and r is a random vector in [0, 1].

3.2.2. Bubble-Net Hunting Routine (Exploitation Phase)

In the bubble net hunting routine two attitudes are cast-off:
  • Shrinking encircling prey is attained by declining the value of a Equation (8) thus the range A is also reduced. A is a random value with the interval of [−α, α] where the value of α declined from 2 to 0 over the time of iterations. Setting random values for A in [−1, 1], the updated position of the search agent is well-defined wherever amid the novel and present best agent;
  • The spiral position update methodology initially calculates the distance between the whale and the target location and then created a spiral equation to mimic the helix-shaped drive of humpbacks.
X ( t + 1 ) = D . e b t . cos ( 2 π l ) + X ( t )
D = | X ( t ) X ( t ) |
where b is constant and l is a random number in [−1, 1].
The mathematical model is as follows:
X ( t + 1 ) = X r a n d A . D       i f     P < 0.5
X ( t + 1 ) = D . e b t . cos ( 2 π l ) + X ( t )     i f     P 0.5
where P is the random number of [0, 1].

3.2.3. Prey’s Searching (Exploration Phase)

Humpback whale arbitrarily searches according to the location of each other. Therefore, we use the random values between −1 to 1. This mechanism | A | > 1 highlights exploration and allows the WOA to achieve a global search. The calculated model is as trails:
D = | C . X r a n d X |
X ( t + 1 ) = X r a n d A . D
Here, X r a n d   is a random position vector selected from the existing population.
From the set of random solutions, the WOA was modified. During iteration, the search agents evaluate their position using the modeling described above. The WOA is a universal optimizer. WOA can easily switch between exploration and exploitation due to the adaptive asymmetry of the searching vector. Furthermore, the WOA only has two parameters that can be changed. WOA’s high exploration ability is due to the whales’ position update system Equation (18). The use of Equations (16) and (17) emphasizes high exploitation and convergence Equation (13). The WOA algorithm is capable of delivering high local optima evasion and convergence speed throughout the iteration sequence, as exposed by the following equations.

3.2.4. Flow Chart of WOA

The flowchart of mother WOA is shown in Figure 4.

3.3. Opposition Enabled Learning (OEL)

Hamid R. Tizoosh 2005 introduced the OEL approach [31]. Rahnamayan et al. introduced the concept of opposition-based learning with the metaheuristic approaches to solving the optimization problems [32,34]. This concept is further explored by Wang et al. in [33]. Many research works are conducted using OEL.
According to OBL, assume that Z is the solution for a given problem; then, the opposite of Z will be the other candidate solution. In this case, the chances to obtain the optimal solution will be increased.
Opposite number: let Z Є (x, y) be a real number; then, the opposite of x is given by:
Z = x + y Z
In Equation (19), Z′ denotes the opposite solution, X denotes the current finest solution, and ‘x’ and ‘y’ are two constants. Consider Z = (Z1, Z2, …. Zd) is a point in the d-dimensional space, where (Z1, Z2, …. Zd) ∈ R and Zd ∈ (xd, yd), d = (1, 2, 3, …. D). The opposite number can be described as Z′ = Z1′, Z2′, …. Zd′. In the d-dimension search space, the above equation can be rewritten as:
Z d = x d + y d Z d
The OEL can be defined as: Assume Z = (Z1, Z2, …. Zd) is a point in the d-dimensional space (i.e., a candidate solution) and f (Z) is a fitness function that is used to evaluate the fitness of candidate solutions. According to the above definitions Z′ = Z1′, Z2′, …. Zd′ is the opposite of Z = (Z1, Z2, …. Zd). If f (Z) is better, then update Z; otherwise, Z′. It is also mentioned that both the Z and Z′ are simultaneously computed and keep the best one. The variables xd and yd denote the minimum and maximum values of the dth dimension.

3.4. Cauchy Mutation (CM) Ooperator in WOA

The CM operator is used to prevent the WOA algorithm from falling into the local optima, especially in the exploration phase. Many researchers have introduced the concept of CM and effectively merged it with metaheuristic algorithms such as PSO [34,35], differential evolutionary (DE) [36,37], and KHA [38]. The idea behind the inclusion of a CM operator with heuristic approaches is to prevent the local optima and to maximize the population diversity. To achieve the same, the best position of the whale is mutated. The Cauchy 1-D density function is expressed by:
f ( x ) = 1 π τ τ 2 + x 2 , < x <
where τ > 0 is a scale parameter [32]. The Cauchy distribution function is:
F τ ( x ) = 1 2 + 1 π arctan ( x τ )
The operator of CM is cast-off in AWOA is explained as follows:
W ( d ) = ( g = 1 t p o p V [ g ] [ d ] ) / t p o p
where V [ g ] [ d ] signifies the velocity vector of the gth cat in the dth dimension, d = 1, 2, … tpop, and W(d) is a weight vector in the range of [−Wmax, Wmax], and Wmax is set to 1 in this.
g b e s t ( d ) = g b e s t ( d ) + W ( d ) C ( X min X max )
where C is the Cauchy distribution function and the scale parameter for it is t = 1, and C ( X m i n X m a x ) is a random number within [ X m i n X m a x ] , that is a defined domain of a test function.

3.5. Amended Whale Optimization Algorithm (AWOA)

This section states the structure of our offered algorithm. To make the WOA algorithm more effective and efficient for strategic bidding problems, a few modifications are inculcated into the conventional WOA algorithm. The opposition enabled a method of learning to encourage whales to migrate towards the finest solutions, while the operator of the CM is capable of introducing population diversity. An efficient balance between exploration and exploitation is also given by this operator. A detailed description of these modifications is described below.

The Opposition Enabled WOA

The OEL concept is integrated with the classical WOA in the opposition-enabled WOA to enhance its search ability. The crucial steps are given below for the OWOA methods.
(Step 1) Firstly, by using random distribution, the initial half population is produced. The other half population (opop) is primed as per the initial half population (pop) in terms of OEL as given in Section 3.3;
(Step 2) In the next step, after initialization, the position of the whale is updated for the initial half of the population (pop), as shown in Section 3.1. As expressed in Section 3.3, according to the initial half of the population (pop), the position of the whale is updated for the remaining half of the population (opop);
(Step 3) Afterward the apprising of the solutions in the population, the 2 sub-populations are collected in one population. This procedure should take the population size unaffected in all of the optimization methods. Moreover, in terms of fitness size, we can sort the population and discover the best one. The technique is then iterated.
The foremost steps of the AWOA procedure are as follows:
  • Step 1: Initialization
Random distribution is used to produce the initial half of the population (pop), which includes tpop/2 individuals. In terms of OWOA, the remaining tpop/2 individuals (opop) are initiated as the initial half population (pop) as illustrated in Section 3.3. In this study, maximum population size is an even number.
  • Step 2: Assessment
Every single individual in the population is ranked as per their location.
  • Step 3: The WOA procedure
The three movements in the WOA approach, as labeled in Section 3.2, modify the positions of the tpop/2 individuals in pop. The following is a description of the key phase in the WOA process:
For d = 1: tpop (all whales in pop) do
Perform the following search scheming.
1. Prey’s encircling
2. Bubble net hunting technique
3. Prey’s searching
After 3, modify the whale’s position in the search space.
  • Step 4: The OWOA procedure
The position of single in the associated tpop, for the last
tpop/2 person in opop, is changed by the OEL rules as indicated in Section 3.3. The following is a description of the primary phase in the OEL process.
For i = tpop/2 + 1: tpop (all krill in opop) do
Compute Z t p o p + d as per Z d   Equation (20)
If f ( Z t p o p + d ) > f ( Z d )
Update Z t p o p / 2 + d by using, Z t p o p + d
End if
End for d
  • Step 5: Amalgamation
After that, all of the entities’ populations (pop and opop) are adjusted, and the pop and opop are merged into a single population.
  • Step 6: The CM operator
For all the entities in the population (pop and opop), the CM operator is performed as exposed in Section 3.4.
  • Step 7: Searching the finest solution
Find out the finest solution ever found and assess the average concert of the Population.
  • Step 8: Halt or not
If the discontinuing criteria (finding the optimal result) are fulfilled, the OWOA algorithm halts and the outcome is the finest solution, otherwise, return to step 2.
The flowchart of the AWOA procedure is shown in Figure 5.

4. Benchmarking of AWOA on Standard Benchmark Functions

To prove the efficacy of developed AWOA, benchmarking of AWOA is described in this section.
The suggested method is benchmarked on 23 common benchmark functions in order to first and foremost demonstrate the superiority of the proposed variation. Figure 6 displays these benchmark functions. These benchmark function’s global optima and solution search boundaries were already defined. By comparing the proposed variant statistically to the parent algorithm, it is possible to assess the exploration and exploitation phenomenon of the AWOA.
Table 2, Table 3 and Table 4 include a list of these benchmark functions. The suggested variant is implemented using MATLAB 2017 [32], which runs at 2.00 GHz on an i5 processor with 8 GB of RAM. The number of iterations and population size are kept constant for all metaheuristic algorithms to allow for a fair comparison of the suggested variant (i.e., a maximum of 500 iterations and 50 search agents). The three groups into which these common benchmark functions are divided are as follows:
(a)
Unimodal functions (G1-G7) (Table 2):-A function g(x) is a unimodal function if for some value m, it is monotonically increasing for x ≤ m and monotonically decreasing for x ≥ m. In that case, the maximum value of g(x) is f(m) and there are no other local maxima;
(b)
Multimodal functions (G8-G13) (Table 3):-A function is said to be multimodal function if it has two or more than two local minima or maxima;
(c)
Fixed dimensions multimodal functions (G14-G23) (Table 4):-A function is said to be multimodal function if it has two or more than two local minima or maxima with fixed dimension.
On each benchmark function, the suggested variation is simulated 30 times. Table 5, Table 6 and Table 7 present the statistical findings (mean and standard deviation). The AWOA algorithm and WOA are compared to validate the results.

4.1. Decisive Evaluation of AWOA on Unimodal Test Functions (G1 to G7)

The benchmarking of the created variant’s exploitation potential is appropriate for unimodal functions. The first seven functions, denoted by G1 through G7, are unimodal functions. The simulation results shown in Table 5 indicate that AWOA can deliver outcomes that are extremely competitive. The study of the findings led to the following conclusions.
The results of the suggested variation on the unimodal benchmark functions are shown in Table 5. Two statistical parameters (mean and standard deviation) from the 30 independent runs were calculated as the basis for the analysis. It can be seen from the results of Table 5 that the mean values for functions G1, G2, G5 and G7 are the best. Therefore, it can be said that OEL has significantly improved the exploitation virtue of the WOA. Additionally, the created variant’s standard deviation values are competitive.
A pairwise non parametric Wilcoxon rank sum test [34] with a 5% threshold of significance was performed to guarantee the results’ statistical significance. It has been noted that the resulting p values for all other opponents are less than 0.05 for functions G1, G2, G5, and G7. This shows that the proposed variant performs well and is statistically distinct from others [35].

4.2. Decisive Evaluation of AWOA on Multi-Modal Test Functions (Exploration Behavior of Functions G8 to G13)

Table 6 shows the result of the multimodal test functions. Multimodal features are taken into consideration to assess the investigation functionality of the systems as these features have numerous local optima’s, the wide variety of which will increase as the scale of trouble will increases. These issues are remarkable to discover the exploration behavior of the algorithm.
The CM operator is employed in AWOA to improve the original WOA’s exploration behavior. It can be seen that the proposed method is capable of improving the original WOA’s exploratory behavior. After that, the proposed AWOA’s convergence is accelerated by the OEL search agent. Exploration and exploitation are effectively blended when these two tactics are combined. The average and standard deviation in the Tables shows that this strategy improves the exploration behavior of the original WOA, as the results for functions G8, G12, and G13 show significant improvement, except for function G11, which shows no improvement, and the results for functions G9 and G10 being the same as the original WOA. The Wilcoxon rank sum test’s p-values signify the validity of innovative AWOA the algorithm is shown in Table 6.
It can be found that the offered algorithm is successful in improving the explorative behavior of the authentic WOA. After that, the opposition-based search agent boosts up the convergence of the offered AWOA. The blend of these techniques offers an influential amalgamation of exploration and exploitation. The average and standard deviation exposed within the Tables, in reality, depict that this strategy is refining the exploration behavior of the novel WOA as the results for functions G8, G12 and G13 signify momentous development compared with function G11 that shows no progress, while the results of features G9 and G10 are equal to those of authentic WOA. The p-values of the Wilcoxon rank sum test in Table 3 show the significance of the proposed AWOA set of rules.

4.3. Decisive Evaluation of AWOA on Fixed Modal Test Functions (Exploration Behavior of Functions G14 to G23)

Table 7 displays the findings of 500 iterations of fixed dimension test functions. These functions are also multimodal, but their dimensions are fixed, whereas the dimensions of multimodal functions can vary depending on the designer’s needs. As a result, their exploration behavior differs from the fixed dimension multimodal capabilities G14 to G23 in some ways. For both AWOA and standard WOA, the outcomes of function G16 are the same. The investigation of function G14 is unchanged, but the remainder of the AWOA functions performs better than the standard WOA. The Wilcoxon rank sum test p-values in Table 7 also show that the suggested AWOA method outperforms the original.

4.4. Evaluation of Algorithm

Three forms of evaluation are performed to verify the efficiency of the algorithm:
  • Convergence Evaluation
  • Statistical Evaluation
  • Data Distribution Evaluation

4.4.1. Convergence Evaluation

To evaluate the convergence efficiency of the proposed variant, some functions have been elected from unimodal, multimodal, and fixed dimensions. The convergence graph in Figure 7 shows that proposed AWOA converges at an expressively improved rate than the conventional WOA algorithm. This evaluation illustrates that the AWOA algorithm is capable of achieving global optima.
These conclusions are drawn from this analysis:
  • From the convergence curves it is observed that the proposed OEL and CM mechanism helps WOA to escape from local minima trap. The outcome of this mechanism emerges as high profit for all the cases;
  • It is observed from the figure that WOA has poor convergence properties for this particular problem as the profit yield by the algorithm is minimal as compared with other algorithms.

4.4.2. Statistical Evaluation with Wilcoxon Rank Sum Test

Wilcoxon rank sum test is used to examine the null hypothesis in WOA and AWOA. In this test, we assume two independent samples of optimization run and compare the distinctness of those samples with the proposed AWOA. Based on the p-values achieved, we define the statistical difference (95% confidence interval) as:
(a) Not significant when p-value > 0.10;
(b) Slightly significant when p-value ≤ 0.10;
(c) Null when p-value = 0.05;
(d) Notable when p-value ≤ 0.05;
(e) Highly notable when p-value ≤ 0.01.
The results of this test have already been depicted in Table 5, Table 6 and Table 7. These results indicate that the proposed AWOA outperforms the majority of the functions and if this test is repeated, then also the proposed AWOA will show the same competitive performance as suggested by Wilcoxon rank sum p-values. Hence, with the results it can be concluded that modification proposed in the algorithm are meaningful and yields statistically diverse results when these tests will be performed again.

4.4.3. Data Distribution Evaluation with Boxplot

Boxplot is a tool that focuses on visualization used for the persistence of relating the distribution for each independent run of an objective function value data taken from the conventional WOA and the proposed AWOA.
These conclusions are drawn from this analysis:
  • In Figure 8, four box plots are drawn for functions G3, G8, G13, and G22. This analysis validates the fitness value distribution for unimodal, multimodal, and fixed multimodal dimensions functions. When opposed to ordinary WOA, the proposed AWOA’s interquartile range and median are lower overall specified benchmark functions;
  • This implies that the output of AWOA fall in a comparatively narrow range as the conventional WOA. The significant enhancement attained by AWOA is due to opposition-based theory and the search capability of dynamic CM operators.

5. Application of AWOA on Strategic Bidding Problem

The inclusion of OEL has significantly improved the exploration and exploitation capabilities of WOA, as is clear from the findings presented in the preceding section. This section now explores how the developed variant can be used to the strategic bidding challenge. We used the proposed version on a case study involving power systems to demonstrate the effectiveness of the proposed variant. This test case simulates a strategic bidding situation where Generating company-k (Genco-k) competes in an auction with four competitor for the IEEE-14 bus test system.
To test the efficacy of the proposed AWOA, a numerical example is presented based on the problem formulation in Section 2. In rapidly evolving environments, this problem is formulated and bidding strategies for a day-ahead market are built for multi-hourly trading. Figure 9 displays a daily load curve for 24 h.
The parameters of all three blocks of Genco-k are given in Table 8 [3].
  • Case I: Normal PDF
The bidding performance of competitors is represented in this test case as illustrated in Figure 2. The opponent bid limit size, mean and standard deviations for all blocks for a normal distribution are provided in [40]. In this case, we used a normal probability distribution. The problem of attaining OBS is solved by AWOA, WOA, GWO, SCA, MFO, GSA and ALO. After successful testing, the results are shown in Figure 10, in this we observed that the outcomes of AWOA are competing and yield more profit for the Genco-k for a multi exchanging hour in a day-ahead market. Optimal block bid price and MCP of Genco-k using a normal distribution with AWOA algorithm shown in Figure 11. In this figure bidding prices of Genco-k in 3 blocks are given with MCP.
The profit curve and MCP obtained by AWOA is shown in Figure 11. From the figure, it is concluded that AWOA outpaces over the competitors as the profit computed by this algorithm is suggestively greater. A steep surge in the profit is perceived in the 10th hour and 12th hour as the profit of the Genco-k becomes $15,478 and $18,667 respectively. Fall in MCP results in fall in profit this phenomenon can be observed in the results of block 2, where the profit reduces drastically at 13th hour. This fall in MCP is observed from (36.87 $/MWh) to (25 $/MWh). At 14th hour again the profit goes more due to the high MCP $32.13. Cumulative profit calculated through AWOA is $180,616.
Table 9 shows the LD calculated using the WOA method for each generator taking part in an auction with a standard MCP. The N-D status in the table denotes units that were not deployed because of a high bid offer. Table 9 also displays the conducted results for the Genco-k from algorithm AWOA using standard PDF.
i.
The third block of Genco-k is not dispatched during the hours with a negative benefit (from 1 to 8 h) due to its high production cost and low system demand;
ii.
Due to the third block’s prolonged shutdown, cold startup costs are included in its production costs when it is committed at nine hours (8 h);
iii.
At the end of 12th h, 3rd block is again non-dispatched due to low system demand, and minimum down time constraint is active (4 h);
iv.
Third block is again dispatched at 15th h, and hot start-up cost is accounted in the making cost, because it has been shut-down for a short time (2 h);
v.
Third block is again non-dispatched from 20 to 24 h due to low system demand;
vi.
Optimal bid price of 3rd Block is shown zero during 1–8 h, 13–14 h, 18 h and 20–24 h, when it is non-dispatched.
  • Case II: Lognormal PDF
In this case study, we used a lognormal probability distribution. After successful testing, the results are shown in Figure 12. The optimal block bid price and MCP of Genco-k using a lognormal distribution with AWOA algorithm are shown in Figure 13.
The profit curve derived by AWOA is shown in Figure 12. This number makes it obvious that AWOA outperforms the rest of the opposition since the profit determined by this algorithm is much higher. The profit of the Genco-k rises dramatically in the 10th and 15th hours, reaching $13,489 and $13,515, respectively. For block 2, the profit decreases sharply in the 13th hour. The total profit as determined by AWOA is $185,362.5.
Table 10 presents the LD obtained from AWOA for all the Genco’s competing in a market under a uniform MCP. In Table 10 N-D signifies the non-dispatched units due to offering of high bids by the Genco. The amount of unit dispatched attained from AWOA using lognormal PDF for the Genco-k are also shown in Table 10.
i.
The third block of Genco-k is non-dispatched in the hours of negative profit (from 1 to 7 h) because of its great production cost and small system load;
ii.
Because the third block has been shut-down for a while, the cost of a cold start-up is included in its production costs when it is dispatched at 8 h (7 h);
iii.
The third Block is once more not dispatched at the end of 12th hours due to low system demand, and the minimum downtime constraint kicks in at 13th hours;
iv.
Due to a short period of shutdown, the third block’s hot start-up cost is included in this hour’s output costs when it is dispatched at 14th h (1h);
v.
The third block is again non-dispatched from 20th h to 24th h due to decrement in the system load;
vi.
Optimal bid price of the third Block is shown as zero during 1–10 h, 13–14 h, 16–24 h, when it is non-dispatched.
  • Case III: Gamma PDF
In this case study, we used a gamma probability distribution for constructing the rival behavior. After successful testing, the results are shown in Figure 14. Optimal block bid price and MCP of Genco-k using gamma distribution with AWOA algorithm are shown in Figure 15.
The profit curve obtained by AWOA is shown in Figure 15. This number makes it obvious that AWOA outperforms the rest of the opposition since the profit determined by this algorithm is much higher. The profit of the Genco-k rises sharply to $15,478 in the tenth h and $18,667 in the 12th h, respectively. For block 2, the profit decreases significantly at the 13th hour as a result of the MCP dropping from 28.51 to 24.03 dollars per megawatt hour. Due to the high MCP $29.39, the profit increases once more at the 14th hour. The total profit as determined by AWOA is $180,616.
Table 11 displays the LD determined by the WOA algorithm for each generator taking part in an auction with a standardized MCP. The N-D status in Table 11 denotes units that were not dispatched because of a high bid offer. Table 11 also displays the transmitted results for the Genco-k from algorithm AWOA utilizing gamma PDF.
i.
The third block of Genco-k is non-dispatched in the hours of negative benefit (from 1 to 9 h) because of its great production cost and small system demand;
ii.
When the third block is dispatched at 10th hours, the cost of a cold start-up is taken into consideration because it has been idle for a while (9 h);
iii.
Due to low system demand at the end of the 12th hour, the third block is once more not dispatched, and the minimum downtime constraint is in effect (3 h);
iv.
The third block is re-dispatched at the 13th hour, and as it was briefly shutdown; the hot start-up cost is included in the production cost of this hour;
v.
Due to less system load, the third block was once again not dispatched from 20 to 24 h;
vi.
The optimal bid price of the third block is shown as zero during 1–9 h, 13–14 h and 18–24 h, when it is non-dispatched.
  • Case IV: Weibull PDF
In this case study, we used the Weibull probability distribution. After successful testing, the results are shown in Figure 16. The optimal block bid price and the MCP of Genco-k using Weibull distribution with AWOA algorithm are shown in Figure 17.
The profit curve obtained by AWOA is shown in Figure 17. This number makes it obvious that AWOA outperforms the rest of the opposition since the profit determined by this algorithm is much higher. The profit of the Genco-k rises sharply to $15,478 in the 10th h and $18,667 in the twelfth hour, respectively. For block 2, the profit drops significantly at the 13th h as a result of the MCP falling from 34 to 21.42 dollars per megawatt hour. Due to the high MCP $31.43 at the fourteenth hour, the profit increases once more. The total profit as determined by AWOA is $180,616.
Table 12 displays the LD determined by the WOA algorithm for each generator taking part in an auction with a standardized MCP. The N-D status in Table 12 denotes units that were not dispatched because of a high bid offer. Table 12 also displays the transmitted results for the Genco-k from algorithm AWOA utilizing Weibull PDF.
i.
Block 3 of Genco-k is not supplied during the hours of adverse benefit due to its high manufacturing cost and less system needs (from 1 to 8 h);
ii.
Because the block 3 has been shut-down for a while, the cost of a cold start-up is included in its production costs when it is dispatched at nine hours (8 h);
iii.
At the conclusion of 12th h, the third block is once more not delivered due to decrement in system demand, and the minimal downtime constraint is in effect (4 h);
iv.
The third block is re-dispatched at 15th h, and because it was shut-down for a brief period of time, the hot start-up cost is included in the cost of production for this hr;
v.
The third block is again non-dispatched from the 20th to 24th h due to less system load;
vi.
Optimal bid price of the third block is exposed zero during 1–8 h, 13–14 h and 20–24 h when it is non-dispatched.
The supremacy of the AWOA method is established through the evaluation of simulation results with former methods. Results of Figure 18 evidence that the finest block bid price resolute by the AWOA delivers high profits than that gained by the former methods such as WOA, GWO, SCA, ALO, MFO and GSA using four different PDFs i.e., Lognormal, Normal, Gamma, and Weibull. Thus, it approves that the AWOA is well skilled in describing the results near-global OBS. After analyzing the dispatched status of the units as suggested by AWOA, it is observed that the AWOA is able to have more dispatched units as compared to other optimization algorithms. Hence, the results of AWOA has been showcased.

6. Conclusions and Future Scope

The EM has been reformed with every passing year and, with the use of recent optimization algorithms, many nonlinear optimization problems of EM have been addressed very efficiently. These algorithms can discover the most efficient strategies for difficult problems, particularly in the EM. In this work, the OBS problem for a Genco in a uniform price spot market is considered. The problem of strategic bidding is solved with the help of a stochastic optimization method by incorporating MC simulation. A new variant of WOA, termed AWOA, has been proposed with the amalgamation of CM operator and OEL techniques. The AWOA is initially benchmarked on standard mathematical functions then it is implemented on OBS for the day-ahead EM by using four different cases. The superiority of the algorithm is examined over recent metaheuristic optimizers such as WOA, GWO, SCA, ALO, MFO and GSA algorithms. It has been revealed that the changes advised in WOA make AWOA a more appropriate technique for all of the cases in a day-ahead EM. Simulation values advise the sturdiness and practicality of the AWOA. The AWOA can be also a good method for Genco in the open power market. The future focus of this research lies in the implementation of the proposed AWOA for formulating an adaptive system with various cases and comparison with the bi-level programming model.

Author Contributions

Conceptualization, A.S. and K.J.; methodology, K.J.; software, K.J. and A.S.; validation, A.F.A., A.M.A. and K.A.A.; formal analysis, A.F.A. and A.M.A.; investigation, A.F.A., A.M.A. and K.A.A.; resources, A.F.A., A.M.A. and K.A.A.; data curation, A.F.A., A.M.A. and K.A.A.; writing—original draft preparation, A.S. and K.J.; writing—review and editing, A.S. and K.J.; visualization, A.F.A., A.M.A. and K.A.A.; supervision, A.S. and A.W.M.; project administration, K.A.A.; funding acquisition, A.F.A., A.M.A. and K.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

The research is funded by Researchers Supporting Program at King Saud University, (RSP-2021/323).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors present their appreciation to King Saud University for funding the publication of this research through Researchers Supporting Program (RSP-2021/323), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

The symbolizations used in this manuscript are specified under:
QmaxCapacity of the lth block of Genco-k [MW].
QminLeast output of lth block of Genco-k [MW].
M l u t Least uptime of lth block of Genco-k [hr].
M l d t Least downtime of lth block of Genco-k [hr].
h c s u Cost for cold start/up [$], measured while the Genco has been shut down for a large time.
c c s u Cost for cold start/up [$], measured while the Genco has been shut down for a large time.
T c Cooling time constant [hr].
H l ( t ) ON The amount of hours the lth block of Genco-k has been unceasingly ON at the last of hour t [hr].
H l ( t ) OFF The amount of hours the lth block of Genco-k has been unceasingly OFF at the last of hour t [hr].
u l ( t ) Binary variable, that is 1, if the lth block is dispatched at hour t; otherwise, 0.
T s u OFF The amount of hours the Genco has been OFF, at the time of startup [hr].

References

  1. David, A.K.; Wen, F. Strategic bidding in competitive electricity markets: A literature survey. In Proceedings of the 2000 Power Engineering Society Summary Meeting, Seattle, WA, USA, 16–20 July 2000; Volume 4, pp. 2168–2173. [Google Scholar]
  2. Kirschen, D.S.; Strbac, G. Fundamentals of Power System Economics; John Wiley & Sons Ltd.: London, UK, 2004. [Google Scholar]
  3. Bajpai, P.; Singh, S.N. Fuzzy Adaptive particle swarm optimization for bidding strategy in uniform price spot market. IEEE Trans. Power Syst. 2007, 22, 2152–2160. [Google Scholar]
  4. Yucekaya, A.D.; Valenzuela, J.; Dozier, G. Strategic bidding in electricity markets using particle swarm optimization. Electr. Power Syst. Res. 2009, 79, 335–345. [Google Scholar]
  5. Azadeh, A.; Skandari, M.; Maleki-Shoja, B. An integrated ant colony optimization approach to compare strategies of clearing market in electricity markets: Agent-based simulation. Energy Policy 2010, 38, 6307–6319. [Google Scholar]
  6. Soleymani, S. Bidding strategy of generation companies using PSO combined with SA method in the pay as bid markets. Int. J. Electr. Power Energy Syst. 2011, 33, 1272–1278. [Google Scholar] [CrossRef]
  7. Azadeh, A.; Ghaderi, S.F.; Nokhandan, B.P.; Sheikhalishahi, M. A new genetic algorithm approach for optimizing bidding strategy viewpoint of profit maximization of a generation company. Expert Syst. Appl. 2012, 39, 1565–1574. [Google Scholar]
  8. Kumar, J.V.; Kumar, D.V.; Edukondalu, K. Strategic bidding using fuzzy adaptive gravitational search algorithm in a pool based electricity market. Appl. Soft Comput. 2013, 13, 2445–2455. [Google Scholar]
  9. Kumar, J.V.; Kumar, D.V. Generation bidding strategy in a pool-based electricity market using shuffled frog leaping algorithm. Appl. Soft Comput. 2014, 21, 407–414. [Google Scholar] [CrossRef]
  10. Shivaie, M.; Ameli, M.T. An environmental/techno-economic approach for bidding strategy in security-constrained electricity markets by a bi-level harmony search algorithm. Renew. Energy 2015, 83, 881–896. [Google Scholar] [CrossRef]
  11. Nojavan, S.; Zare, K.; Ashpazi, M.A. A hybrid approach based on IGDT–MPSO method for optimal bidding strategy of price-taker generation station in day-ahead electricity market. Int. J. Electr. Power Energy Syst. 2015, 69, 335–343. [Google Scholar]
  12. Saxena, A.; Soni, B.P.; Kumar, R.; Gupta, V. Intelligent Grey Wolf Optimizer—Development and application for strategic bidding in uniform price spot energy market. Appl. Soft Comput. 2018, 69, 1–13. [Google Scholar] [CrossRef]
  13. Karri, C.; Rajababu, D.; Raghuram, K. Optimal Bidding Strategy in Deregulated Power Market Using Krill Herd Algorithm. In Applications of Artificial Intelligence Techniques in Engineering; Springer: Singapore, 2018; pp. 43–51. [Google Scholar]
  14. Sudhakar, A.V.V.; Karri, C.; Laxmi, A.J. Optimal Bidding Strategy in Deregulated Power Market Using Invasive Weed Optimization. In Applications of Artificial Intelligence Techniques in Engineering; Springer: Singapore, 2018; pp. 421–429. [Google Scholar]
  15. Fanzeres, B.; Ahmed, S.; Street, A. Robust strategic bidding in auction-based markets. Eur. J. Oper. Res. 2019, 272, 1158–1172. [Google Scholar]
  16. Jain, K.; Saxena, A. Evolutionary Neural Network based hybrid architecture for strategic bidding in electricity market. In Proceedings of the 2021 IEEE 2nd International Conference on Smart Technologies for Power, Energy and Control (STPEC), Bilaspur, India, 19–22 December 2021. [Google Scholar]
  17. Yang, Y.; Ji, T.; Jing, Z. Selective learning for strategic bidding in uniform pricing electricity spot market. CSEE J. Power Energy Syst. 2021, 7, 1334–1344. [Google Scholar]
  18. Jain, K.; Jasser, M.B.; Hamzah, M.; Saxena, A.; Mohamed, A.W. Harris Hawk Optimization-Based Deep Neural Networks Architecture for Optimal Bidding in the Electricity Market. Mathematics 2022, 10, 2094. [Google Scholar] [CrossRef]
  19. Iria, J.; Soares, F.; Matos, M. Optimal supply and demand bidding strategy for an aggregator of small prosumers. Appl. Energy 2018, 213, 658–669. [Google Scholar]
  20. Mohsenian-Rad, H. Optimal demand bidding for time-shiftable loads. IEEE Trans. Power Syst. 2014, 30, 939–951. [Google Scholar]
  21. Sadeghi-Mobarakeh, A.; Mohsenian-Rad, H. Optimal bidding in performance-based regulation markets: An MPEC analysis with system dynamics. IEEE Trans. Power Syst. 2016, 32, 1282–1292. [Google Scholar]
  22. Iria, J.; Soares, F.; Matos, M. Optimal bidding strategy for an aggregator of prosumers in energy and secondary reserve markets. Appl. Energy 2019, 238, 1361–1372. [Google Scholar]
  23. Mirjalili, S.; Lewis, A. The Whale Optimization Algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar]
  24. Ding, Y.; Chen, Z.; Zhang, H.; Wang, X.; Guo, Y. A short-term wind power prediction model based on CEEMD and WOA-KELM. Renew. Energy 2022, 189, 188–198. [Google Scholar] [CrossRef]
  25. Singh, A.; Khamparia, A. A hybrid whale optimization-differential evolution and genetic algorithm based approach to solve unit commitment scheduling problem: WODEGA. Sustain. Comput. Inform. Syst. 2020, 28, 100442. [Google Scholar]
  26. Sharma, A.K.; Saxena, A. A demand side management control strategy using Whale optimization algorithm. SN Appl. Sci. 2019, 1, 870. [Google Scholar]
  27. Hassan, N.M.; Swief, R.A.; Kamh, M.Z.; Hasanien, H.M.; Abdelaziz, A.Y. Centralized/decentralized optimal load flow based on tuned whale optimization algorithm. In Proceedings of the 2020 International Conference on Innovative Trends in Communication and Computer Engineering (ITCE), Aswan, Egypt, 8–9 February 2020. [Google Scholar]
  28. Shahbudin, I.S.; Musirin, I.; Suliman, S.; Harun, A.F.; Mustaffa, S.A.S.; Suyono, H.; Ghani, N.A. FACTS device installation in transmission system using whale optimization algorithm. Bull. Electr. Eng. Inform. 2019, 8, 30–38. [Google Scholar]
  29. Tizhoosh, H.R. Opposition-based learning: A new scheme for Machine Intelligence. In Proceedings of the Computational Intelligence for Modeling, Control and Automation, and International Conference on Intelligent Agents, Web Technologies and Internet Commerce, Vienna, Austria, 28–30 November 2005; Volume 1, pp. 695–701. [Google Scholar]
  30. Rahnamayan, S.; Tizhoosh, H.R.; Salama, M.M.A. Opposition-based differential evolution algorithms. In Proceedings of the IEEE Congress on Evolutionary Computation, Vancouver, BC, Canada, 16–21 July 2006; pp. 2010–2017. [Google Scholar]
  31. Xu, Q.; Wang, L.; Wang, N.; Hei, X.; Zhao, L. A review of opposition-based learning from 2005 to 2012. Eng. Appl. Artif. Intell. 2014, 29, 1–12. [Google Scholar]
  32. Wang, H.; Li, H.; Liu, Y.; Li, C.; Zeng, S. Opposition—Based particle swarm algorithm with Cauchy mutation. In Proceedings of the 2007 IEEE Congress on Evolutionary Computation (CEC2007), Singapore, 25–28 September 2007; pp. 4750–4756. [Google Scholar]
  33. Wu, Q. Hybrid forecasting model based on support vector machine and particle swarm optimization with adaptive and Cauchy mutation. Expert Syst. Appl. 2011, 38, 9070–9075. [Google Scholar]
  34. Qin, H.; Zhou, J.; Lu, Y.; Wang, Y.; Zhang, Y. Multi-objective differential evolution with adaptive Cauchy mutation for short-term multi-objective optimal hydro- thermal scheduling. Energy Convers. Manag. 2010, 51, 788–794. [Google Scholar]
  35. Ali, M.; Pant, M. Improving the performance of differential evolution algorithm using Cauchy mutation. Soft Comput. 2010, 15, 991–1007. [Google Scholar]
  36. Wang, G.-G.; Deb, S.; Gandomi, A.H.; Alavi, A.H. Opposition-based krill herd algorithm with Cauchy mutation and position clamping. Neurocomputing 2016, 177, 147–157. [Google Scholar]
  37. Bialek, W.; Callan, C.G.; Strong, S.P. Field theories for learning probability distributions. Phys. Rev. Lett. 1996, 77, 4693–4697. [Google Scholar]
  38. Boyle, P.P. Options: A Monte Carlo approach. J. Financ. Econ. 1977, 4, 323–338. [Google Scholar]
  39. Kroese, D.P.; Brereton, T.; Taimre, T.; Botev, Z.I. Why the monte carlo method is so important today. Wiley Interdiscip. Rev. Comput. Stat. 2014, 6, 386–392. [Google Scholar]
  40. Wood, A.J.; Wollenberg, B.F. Power Generation Operation and Control; Wiley: New York, NJ, USA, 1996. [Google Scholar]
  41. Walters, D.C.; Sheble, G.B. Genetic algorithm solution of economic dispatch with valve point loading. IEEE Trans. Power Syst. 1993, 8, 1325–1332. [Google Scholar] [CrossRef]
  42. Ross, D.W.; Kim, S. Dynamic economic dispatch of generation. IEEE Trans. Power Appar. Syst. 1980, PAS-99, 2060–2068. [Google Scholar] [CrossRef]
Figure 1. Paper Contribution Chart.
Figure 1. Paper Contribution Chart.
Axioms 11 00456 g001
Figure 2. Cost Modal of Different Distributions.
Figure 2. Cost Modal of Different Distributions.
Axioms 11 00456 g002
Figure 3. Bubble Net Feeding Humpback Whale.
Figure 3. Bubble Net Feeding Humpback Whale.
Axioms 11 00456 g003
Figure 4. Flow Chart of WOA.
Figure 4. Flow Chart of WOA.
Axioms 11 00456 g004
Figure 5. Flowchart of AWOA.
Figure 5. Flowchart of AWOA.
Axioms 11 00456 g005
Figure 6. Shape Curve of Benchmark Function.
Figure 6. Shape Curve of Benchmark Function.
Axioms 11 00456 g006
Figure 7. Convergence Behavior of AWOA and WOA.
Figure 7. Convergence Behavior of AWOA and WOA.
Axioms 11 00456 g007
Figure 8. Boxplots for Distribution of Objective Function Values.
Figure 8. Boxplots for Distribution of Objective Function Values.
Axioms 11 00456 g008
Figure 9. Power Load Forecast by the System Operator.
Figure 9. Power Load Forecast by the System Operator.
Axioms 11 00456 g009
Figure 10. Cumulative Profit Calculated of the Algorithms using Normal PDF.
Figure 10. Cumulative Profit Calculated of the Algorithms using Normal PDF.
Axioms 11 00456 g010
Figure 11. AWOA (Normal Distribution) Block Bid Price and MCP.
Figure 11. AWOA (Normal Distribution) Block Bid Price and MCP.
Axioms 11 00456 g011
Figure 12. Cumulative Profit Calculated of the Algorithms using Lognormal PDF.
Figure 12. Cumulative Profit Calculated of the Algorithms using Lognormal PDF.
Axioms 11 00456 g012
Figure 13. AWOA (Lognormal Distribution) Block Bid Price and MCP.
Figure 13. AWOA (Lognormal Distribution) Block Bid Price and MCP.
Axioms 11 00456 g013
Figure 14. Cumulative Profit Calculated of the Algorithms using Gamma PDF.
Figure 14. Cumulative Profit Calculated of the Algorithms using Gamma PDF.
Axioms 11 00456 g014
Figure 15. AWOA (Gamma Distribution Block Bid Price and MCP.
Figure 15. AWOA (Gamma Distribution Block Bid Price and MCP.
Axioms 11 00456 g015
Figure 16. Cumulative Profit Calculated of the Algorithms using Weibull PDF.
Figure 16. Cumulative Profit Calculated of the Algorithms using Weibull PDF.
Axioms 11 00456 g016
Figure 17. AWOA (Weibull Distribution) Block Bid Price and MCP.
Figure 17. AWOA (Weibull Distribution) Block Bid Price and MCP.
Axioms 11 00456 g017
Figure 18. Total Profit of 24 Hours Using Different Approaches.
Figure 18. Total Profit of 24 Hours Using Different Approaches.
Axioms 11 00456 g018
Table 1. Data of Competitor’s Bidding Parameters.
Table 1. Data of Competitor’s Bidding Parameters.
BLOCK 1BLOCK 3BLOCK 3
QIµinσinQIµinσinQIµinσin
RIVAL 1200102.5300203400303
RIVAL 2300153400302500504
RIVAL 3250102300152.5300202.5
RIVAL 4300204350255450405
Table 2. Unimodal Benchmark Functions.
Table 2. Unimodal Benchmark Functions.
FunctionDimRangeMin. Value
G 1 ( x ) = i = 1 n x i 2 30[−100, 100]0
G 2 ( x ) = i = 1 n | x i | + i = 1 n | x i | 30[−10, 10]0
G 3 ( x ) = i = 1 n ( j 1 i x j ) 2 30[−100, 100]0
G 4 ( x ) = m A x i { | x i | , 1 i n } 30[−100, 100]0
G 5 ( x ) = i = 1 n 1 [ 100 ( x i + 1 x i 2 ) 2 + ( x i 1 ) 2 ] 30[−30, 30]0
G 6 ( x ) = i = 1 n 1 ( x i + 0.5 ) 2 30[−100, 100]0
G 7 ( x ) = i = 1 n 1 i x i 4 + r a n d o m [ 0 , 1 ] 30[−1.28, 1.28]0
Table 3. Multimodal Benchmark Functions.
Table 3. Multimodal Benchmark Functions.
FunctionDimRangeMin. Value
G 8 ( x ) = i = 1 n x i sin ( | x i | ) 30[−500, 500]−418.9829 × 5
G 9 ( x ) = i = 1 n [ x i 2 10 c o s ( 2 π x i + 10 ) ] 30[−5.12, 5.12]0
G 10 ( x ) = 20 e x p ( 0.2 1 n i = 1 n x i 2 ) exp ( 1 n i = 1 n cos ( 2 π x i ) ) + 20 + e 30[−32, 32]0
G 11 ( x ) = 1 4000 i = 1 n x i 2 i = 1 n cos ( x i i ) + 1 30[−600, 600]0
G 12 ( x ) = π n { 10 sin ( π y i ) + i = 1 n 1 ( y i 1 ) 2 [ 1 + 10 s i n 2 ( π y i + 1 ) + ( y n 1 ) 2 ] } + i = 1 n u ( x i , 10 , 100 , 4 )
y i = 1 + x i + 1 4
u ( X i , A , k , m ) = { k ( X i A ) m x > a 0 a < x i < a k ( x i a ) m x i a }
30[−50, 50]0
G 13 ( x ) = 0.1 { s i n 2 ( 3 π x 1 ) + i = 1 n ( x i 1 ) 2 [ 1 + sin 2 ( 3 π x i + 1 ) ] + ( x n 1 ) 2 [ 1 + sin 2 ( 2 π x n ) ] } + i = 1 n u ( x i , 5 , 100 , 4 ) 30[−50, 50]0
Table 4. Fixed-Dimension Multimodal Benchmark Functions.
Table 4. Fixed-Dimension Multimodal Benchmark Functions.
FunctionDimRangeMin. Value
G 14 ( x ) = i = 1 n sin ( x i ) . ( sin ( i x i 2 π ) ) 2 m , m = 10 2[−65, 65]1
G 15 ( x ) = [ e i = 1 n ( x i β ) 2 m 2 e i = 1 n x i 2 ] . i = 1 n cos 2 x i , m = 5 4[−5, 5]0.00030
G 16 ( x ) = { [ i = 1 n sin 2 ( x i ) ] e x p ( i = 1 n x i 2 ) } . e x p [ i = 1 n sin 2 | x i | ] 2[−5, 5]−1.0316
G 17 ( x ) = ( x 2 5.1 4 π 2 x 1 2 + 5 π x 1 6 ) 2 + 10 ( 1 1 8 π ) c o s x 1 + 10 2[−5, 5]0.398
G 18 ( x ) = [ 1 + ( x 1 + x 2 + 1 ) 2 ( 19 14 x 1 + 3 x 1 2 14 x 2 + 6 x 1 x 2 + 3 x 2 2 ] [ 30 + ( 2 x 1 3 x 2 ) 2 ( 18 32 x 1 + 12 x 1 2 + 48 x 2 36 x 1 x 2 + 27 x 2 2 ) ] 2[−2, 2]3
G 19 ( x ) = i = 1 4 c i exp ( j = 1 3 a i j ( x j p i j ) 2 ) 3[1, 3]−3.86
G 20 ( x ) = i = 1 4 c i exp ( j = 1 6 a i j ( x j p i j ) 2 ) 6[0, 1]−3.32
G 21 ( x ) = i = 1 5 [ ( X A i ) ( X A i ) T + c i ] 1 4[0, 10]−10.1532
G 22 ( x ) = i = 1 7 [ ( X A i ) ( X A i ) T + c i ] 1 4[0, 10]−10.4028
G 23 ( x ) = i = 1 10 [ ( X A i ) ( X A i ) T + c i ] 1 4[0, 10]−10.5363
Table 5. Average, Standard Deviation, and p-Value of AWOA and other comparable algorithms for Uni-modal Functions.
Table 5. Average, Standard Deviation, and p-Value of AWOA and other comparable algorithms for Uni-modal Functions.
FunctionsStatistical
Parameters
Algorithms
AWOAWOASCAGWOALOMFOGSA
G1Mean1.91 × 10−1032.35 × 10−721.29 × 10−123.82 × 10−332.09 × 10−032.92 × 10−033.96 × 10−17
St.dev3.53 × 10−1044.46 × 10−736.55 × 10−128.63 × 10−331.90 × 10−033.63 × 10−131.33 × 10−17
p-valueNA3.02 × 10−117.22 × 10−016.20 × 10−018.97 × 10−011.88 × 10−031.48 × 10−07
G2Mean2.45 × 10−666.89 × 10−501.16 × 10−096.30 × 10−206.17 × 10+011.673.22 × 10−08
St.dev4.47 × 10−671.26 × 10−502.94 × 10−094.51 × 10−205.11 × 10+014.616.53 × 10−09
p-valueNA3.02 × 10−117.99 × 10−016.50 × 10−012.36 × 10−017.89 × 10−026.75 × 10−01
G3Mean9.16 × 10+045.82 × 10+044.17 × 10−042.60 × 10−085.02 × 10+031.676.00 × 10+02
St.dev1.47 × 10+041.06 × 10+041.72 × 10−037.06 × 10−081.49 × 10+031.73 × 10+031.80 × 10+02
p-valueNA2.44 × 10−096.64 × 10−027.64 × 10−027.84 × 10−011.23 × 10−026.80 × 10−08
G4Mean8.83 × 10+019.30 × 10+012.72 × 10−042.07 × 10−082.02 × 10+013.653.25
St.dev2.90 × 10+012.84 × 10+018.69 × 10−041.78 × 10−084.425.961.38
p-value8.42 × 10−01NA4.69 × 10−024.09 × 10−013.69 × 10−024.78 × 10−016.80 × 10−08
G5Mean2.87 × 10+012.88 × 10+013.13 × 10+012.65 × 10+012.07 × 10+023.27 × 10+033.51 × 10+01
St.dev4.985.01 × 10−015.51 × 10−015.35 × 10−012.31 × 10+021.64 × 10+042.26 × 10+01
p-valueNA5.57 × 10−103.60 × 10−021.20 × 10−064.60 × 10−026.90 × 10−021.80 × 10−06
G6Mean2.78 × 10−018.11 × 10−013.85 × 10−015.12 × 10−011.13 × 10−031.19 × 10−132.33 × 10−01
St.dev7.02 × 10−022.33 × 10−011.14 × 10−013.26 × 10−015.35 × 10−041.60 × 10−134.30 × 10−01
p-valueNA3.19 × 10−097.56 × 10−031.64 × 10−011.26 × 10−017.90 × 10−026.80 × 10−08
G7Mean1.08 × 10−041.47 × 10−021.59 × 10−031.34 × 10−032.93 × 10−017.58 × 10−032.67 × 10−02
St.dev1.13 × 10−033.19 × 10−031.39 × 10−038.20 × 10−031.43 × 10−014.66 × 10−031.26 × 10−02
p-valueNA1.15 × 10−074.48 × 10−023.23 × 10−011.63 × 10−031.70 × 10−026.66 × 10−08
Table 6. Average, Standard Deviation, and p-Value of AWOA and other comparable algorithms for Multimodal Functions.
Table 6. Average, Standard Deviation, and p-Value of AWOA and other comparable algorithms for Multimodal Functions.
FunctionsStatistical
Parameters
Algorithms
AWOAWOASCAGWOALOMFOGSA
G8Mean−9.89 × 10+03−6.57 × 10+03−2.18 × 10+03−6.03 × 10+03−5.49 × 10+03−3.31 × 10+03−2.83 × 10+03
St.dev6.01 × 10+021.83 × 10+031.44 × 10+027.40 × 10+028.41 × 10+013.59 × 10+023.24 × 10+02
p-valueNA3.02 × 10−111.79 × 10−033.37 × 10−011.48 × 10−038.96 × 10−023.79 × 10−01
G9Mean0.000.007.41 × 10−012.007.29 × 10+012.31 × 10+011.81 × 10+01
St.dev0.000.003.743.031.75 × 10+011.21 × 10+013.71
p-value1.62 × 10−01NA2.40 × 10−031.35 × 10−054.57 × 10−011.24 × 10−028.29 × 10−08
G10Mean7.99 × 10−157.99 × 10−156.44 × 10−044.39 × 10−144.386.71 × 10−024.66 × 10−09
St.dev2.03 × 10−152.09 × 10−153.53 × 10−034.70 × 10−153.213.68 × 10−016.96 × 10−10
p-value7.85 × 10−01NA6.80 × 10−027.35 × 10−014.79 × 10−027.80 × 10−026.80 × 10−08
G11Mean2.14 × 10−012.76 × 10−013.18 × 10−015.42 × 10−016.05 × 10−022.58 × 10−011.74 × 10+01
St.dev3.01 × 10−024.12 × 10−021.92 × 10−013.95 × 10−023.28 × 10−023.25 × 10−023.50
p-value4.55 × 10−01NA1.60 × 10−045.40 × 10−011.46 × 10−011.25 × 10−021.60 × 10−05
G12Mean1.94 × 10−026.65 × 10−028.71 × 10−022.73 × 10−021.92 × 10+011.87 × 10−016.04 × 10−01
St.dev3.73 × 10−031.46 × 10−024.64 × 10−021.01 × 10−021.01 × 10+014.59 × 10−015.79 × 10−01
p-valueNA4.31 × 10−081.12 × 10−031.20 × 10−031.48 × 10−027.89 × 10−031.16 × 10−04
G13Mean4.45 × 10−021.422.73 × 10−013.82 × 10−012.39 × 10+015.14 × 10−021.98
St.dev1.04 × 10−013.59 × 10−011.57 × 10−011.69 × 10−011.32 × 10+011.64 × 10−012.21
p-valueNA3.02 × 10−118.97 × 10−032.40 × 10−023.65 × 10−011.46 × 10−011.56 × 10−04
Table 7. Average, Standard Deviation, and p-Value of AWOA and other comparable algorithms for Multi-Modal Functions with Fixed Dimension.
Table 7. Average, Standard Deviation, and p-Value of AWOA and other comparable algorithms for Multi-Modal Functions with Fixed Dimension.
FunctionsStatistical
Parameters
Algorithms
AWOAWOASCAGWOALOMFOGSA
G14Mean1.08 × 10+015.931.462.761.402.315.33
St.dev2.031.208.53 × 10−013.296.94 × 10−012.013.94
p-value1.001.103.58 × 10−014.80 × 10−03N/A1.24 × 10−010.425
G15Mean2.25 × 10−048.27 × 10−031.02 × 10−033.02 × 10−032.79 × 10−039.78 × 10−043.75 × 10−03
St.dev2.57 × 10−041.46 × 10−033.73 × 10−046.92 × 10−032.89 × 10−033.48 × 10−042.05 × 10−03
p-value6.86 × 10−016.76 × 10−011.80 × 10−036.80 × 10−031.65 × 10−03N/A2.62 × 10−01
G16Mean−1.03−1.03−1.03−1.03−1.03−1.03−1.03
St.dev1.94 × 10−087.62 × 10−102.64 × 10−059.00 × 10−092.79 × 10−126.78 × 10−165.22 × 10−16
p-value6.67 × 10−01NA4.88 × 10−036.71 × 10−029.65 × 10−031.80 × 10−021.35 × 10−03
G17Mean3.98 × 10−013.98 × 10−013.98 × 10−013.98 × 10−013.98 × 10−013.98 × 10−013.98 × 10−01
St.dev2.41 × 10−063.01 × 10−041.63 × 10−031.50 × 10−061.34 × 10−130.000.00
p-valueNA7.00 × 10−011.97 × 10−023.94 × 10−021.48 × 10−033.21 × 10−021.23 × 10−07
G18Mean3.003.003.003.003.003.003.00
St.dev7.74 × 10−059.55 × 10−053.28 × 10−051.38 × 10−056.86 × 10−132.54 × 10−151.25 × 10−02
p-valueNA9.37 × 10−015.70 × 10−026.80 × 10−041.47 × 10−024.14 × 10−021.79 × 10−04
G19Mean−3.66−3.69−3.86−3.86−3.86−3.86−3.86
St.dev3.73 × 10−023.88 × 10−022.81 × 10−031.74 × 10−031.78 × 10−132.71 × 10−153.86 × 10−02
p-value7.00 × 10−01NA3.59 × 10−014.68 × 10−054.19 × 10−017.96 × 10−028.59 × 10−02
G20Mean−2.71−2.84−2.91−3.32−3.27−3.22−3.32
St.dev1.08 × 10−021.17 × 10−013.07 × 10−011.19 × 10−016.16 × 10−029.65 × 10−021.19 × 10−01
p-value9.37 × 10−017.37 × 10−011.46 × 10−017.90 × 10−05N/A4.57 × 10−027.56 × 10−01
G21Mean−8.94−2.63−3.25−6.34−6.37−6.13−6.34
St.dev2.21 × 10−012.521.753.662.723.263.66
p-valueNA2.86 × 10−022.37 × 10−021.77 × 10−062.31 × 10−011.80 × 10−029.03 × 10−01
G22Mean4.25 × 10−01−3.84−3.70−9.97−5.10−7.42 −9.97
St.dev1.213.201.801.673.011.21 × 10+011.67
p-valueNA2.86 × 10−027.42 × 10−021.92 × 10−057.90 × 10−012.39 × 10−013.94 × 10−01
Table 8. Data Of Genco-K Power Blocks.
Table 8. Data Of Genco-K Power Blocks.
C0 (MW2h)C1
($/MWh)
C2 ($/h)C3 ($/h)C4
(rad./MW)
Qmax (MW)Qmin (MW)MUT (h)MDT
(h)
h ($)δ ($)τ (h) c i d
($)
Block10.004827.97781500.0632005011100015001100
Block20.0019415.853102000.04240010011150025001200
Block30.00156232.925613000.031560010011200040008400
Table 9. Dispatch Obtained by Different Algorithms Using Normal Distribution for Genco-k.
Table 9. Dispatch Obtained by Different Algorithms Using Normal Distribution for Genco-k.
LOAD DISPATCH OF NORMAL DISTRIBUTION OF 5 GENCO AN-D 3 BLOCK
HOURLOADRival 1Rival 2Rival 3Rival 4GENCO-R
B1B2B3B1B2B3B1B2B3B1B2B3B1B2B3
200300400300400500250300300300350450200400600
11500200N-DN-D300N-DN-D250N-DN-D150N-DN-D200400N-D
21500200N-DN-D150N-DN-D250300N-DN-DN-DN-D200400N-D
31500200N-DN-D300N-DN-D250N-DN-D150N-DN-D200400N-D
41500200150N-D300N-DN-D250N-DN-DN-DN-DN-D200400N-D
5200020050N-D300N-DN-D250300300N-DN-DN-D200400N-D
62000200300N-D300N-DN-D250N-D300300N-DN-D200400N-D
72000200N-DN-D300N-DN-D25030030050N-DN-D200400N-D
82500200300N-D300N-DN-D250300250300N-DN-D200400N-D
93000200300N-D300400N-D250300N-D300350N-D200400N-D
103500200300400300400N-D250300300300N-DN-D200400150
113500200300200300400N-D250300300300350N-D200400N-D
123500200300N-D300400N-D250300300300350N-D200400200
132500200300N-D300N-DN-D250300N-D300350N-D200400400
143000200300400300N-DN-D250300N-D300350N-D200400N-D
153500200300400300N-DN-D250300300300350N-D200400200
163500200300350300250N-D250300300300350N-D200400N-D
173500200300400300200N-D250300300300350N-D200400N-D
183000200300N-D300100N-D250300300300350N-D200400N-D
193000200300N-D300400N-D250300300300N-DN-D20040050
202500200300N-D300N-DN-D250300300250N-DN-D200400N-D
212000200N-DN-DN-DN-DN-D250300N-D300350N-D200400N-D
222000200N-DN-D300N-DN-D250300N-DN-D350N-D200400N-D
231500200N-DN-D300N-DN-D250N-DN-D300N-DN-D200250N-D
241500200N-DN-D300N-DN-D250250N-DN-DN-DN-D200300N-D
Table 10. Dispatch Obtained by Different Algorithms by using Lognormal PDF for Genco-K.
Table 10. Dispatch Obtained by Different Algorithms by using Lognormal PDF for Genco-K.
LOAD DISPATCH OF LOGNORMAL DISTRIBUTION OF 5 GENCO AN-D 3 BLOCK
HOURLOADRival 1Rival 2Rival 3Rival 4GENCO-R
B1B2B3B1B2B3B1B2B3B1B2B3B1B2B3
200300400300400500250300300300350450200400600
11500200100N-D300N-DN-DN-D300N-DN-DN-DN-D200400N-D
21500200N-DN-D300N-DN-D100300N-DN-DN-DN-D200400N-D
31500200100N-D300N-DN-DN-D300N-DN-DN-DN-D200400N-D
41500200N-DN-D300N-DN-DN-D300N-D100N-DN-D200400N-D
52000200300N-D300N-DN-D250300N-DN-D50N-D200400N-D
62000200N-DN-D300N-DN-D250300N-D30050N-D200400N-D
72000200N-DN-D300N-DN-D250300N-DN-D350N-D200400N-D
82500200300N-D300N-DN-D250300N-D300250N-D200400N-D
9300020030040030050N-D250300300300N-DN-D200400N-D
103500200300200300400N-D250300300300350N-D200400N-D
113500200300400300400N-D250300300300N-DN-D200400150
123500200300400300N-DN-D250300300300350N-D200400200
132500200300N-D300N-DN-D250300N-D300250N-D200400N-D
143000200300400300N-DN-D25030030030050N-D200400N-D
153500200300400300N-DN-D250300300300350N-D200400200
163500200300400300200N-D250300300300350N-D200400N-D
173500200300200300400N-D250300300300350N-D200400N-D
183000200300400300N-DN-D250300N-D300350N-D200400N-D
193000200300400300N-DN-D250300N-D300350N-D200400N-D
202500200300N-DN-DN-DN-D250300300300250N-D200400N-D
21200020050N-D300N-DN-D250300N-D300N-DN-D200400N-D
222000200N-DN-D300N-DN-D25030030050N-DN-D200400N-D
231500200N-DN-D300N-DN-D100300N-DN-DN-DN-D200400N-D
241500200N-DN-D300N-DN-D100300N-DN-DN-DN-D200400N-D
Table 11. Dispatch Obtained by Different Algorithms by using Gamma PDF for Genco-K.
Table 11. Dispatch Obtained by Different Algorithms by using Gamma PDF for Genco-K.
LOAD DISPATCH OF GAMMA DISTRIBUTION OF 5 GENCO AN-D 3 BLOCK
HOURLOADRival 1Rival 2Rival 3Rival 4GENCO-R
B1B2B3B1B2B3B1B2B3B1B2B3B1B2B3
200300400300400500250300300300350450200400600
11500200N-DN-D150N-DN-D250300N-DN-DN-DN-D200400N-D
21500200N-DN-DN-DN-DN-D250300N-D150N-DN-D200400N-D
31500200N-DN-D300N-DN-D250300N-DN-DN-DN-D200250N-D
41500200N-DN-D300N-DN-D250300N-DN-DN-DN-D200250N-D
52000200300N-D300N-DN-D25030050N-DN-DN-D200400N-D
62000200N-DN-D300N-DN-D25030050300N-DN-D200400N-D
72000200300N-D300N-DN-D250300N-D50N-DN-D200400N-D
82500200150N-D300N-DN-D250300300300100N-D200400N-D
93000200300100300N-DN-D250300300300350N-D200400N-D
103500200300400300200N-D250300300300350N-D200400N-D
113500200300400300200N-D250300300300350N-D200400N-D
123500200300400300200N-D250300300300350N-D200400N-D
132500200300N-D300N-DN-D250300250300N-DN-D200400N-D
143000200300N-D300100N-D250300300300350N-D200400N-D
153500200300200300400N-D250300300300350N-D200400N-D
163500200300400300200N-D250300300300350N-D200400N-D
173500200300400300200N-D250300300300350N-D200400N-D
183000200300N-D300100N-D250300300300350N-D200400N-D
193000200300100300N-DN-D250300300300350N-D200400N-D
202500200300N-D300N-DN-D250300300300N-DN-D200350N-D
21200020050N-D300N-DN-D250300N-D300N-DN-D200400N-D
222000200N-DN-D300N-DN-D250300N-DN-D350N-D200400N-D
231500200N-DN-D300N-DN-D250300N-DN-DN-DN-D200250N-D
241500200N-DN-D300N-DN-D250300N-DN-DN-DN-D200250N-D
Table 12. Dispatch Obtained by Different Algorithms by using Weibull PDF for Genco-K.
Table 12. Dispatch Obtained by Different Algorithms by using Weibull PDF for Genco-K.
LOAD DISPATCH OF WEIBULL DISTRIBUTION OF 5 GENCO AND 3 BLOCK
HOURLOADRival 1Rival 2Rival 3Rival 4GENCO-R
B1B2B3B1B2B3B1B2B3B1B2B3B1B2B3
200300400300400500250300300300350450200400600
11500200NDND300NDND250300NDNDNDND200250ND
21500200NDND300NDND250300NDNDNDND200250ND
31500200NDND300NDND250300NDNDNDND200250ND
41500200NDNDNDNDND250300300NDNDND200250ND
52000200ND300300NDND250300ND50NDND200400ND
62000200NDND300NDND25030050300NDND200400ND
72000200300ND300NDND25030050NDNDND200400ND
82500200300ND300NDND250300250300NDND200400ND
93000200300400300NDND250300ND300350ND200400ND
103500200300400300200ND250300300300350ND200400ND
113500200300200300400ND250300300300350ND200400ND
123500200300200300400ND250300300300350ND200400ND
132500200300ND300NDND250300250300NDND200400ND
143000200300100300NDND250300300300350ND200400ND
153500200300400300200ND250300300300350ND200400ND
163500200300ND300400ND250300300300350ND200400200
173500200300200300400ND250300300300350ND200400ND
183000200300ND300NDND250300300300350ND200400100
193000200300100300NDND250300300300350ND200400ND
202500200300ND300NDND250300250300NDND200400ND
212000200300ND300NDND250300ND300NDND200150ND
22200020050ND300NDND250300ND300NDND200400ND
231500200NDND300NDND250300NDNDNDND200250ND
241500200NDND300NDND250300NDNDNDND200250ND
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jain, K.; Saxena, A.; Alshamrani, A.M.; Alrasheedi, A.F.; Alnowibet, K.A.; Mohamed, A.W. An Amended Whale Optimization Algorithm for Optimal Bidding in Day Ahead Electricity Market. Axioms 2022, 11, 456. https://doi.org/10.3390/axioms11090456

AMA Style

Jain K, Saxena A, Alshamrani AM, Alrasheedi AF, Alnowibet KA, Mohamed AW. An Amended Whale Optimization Algorithm for Optimal Bidding in Day Ahead Electricity Market. Axioms. 2022; 11(9):456. https://doi.org/10.3390/axioms11090456

Chicago/Turabian Style

Jain, Kavita, Akash Saxena, Ahmad M. Alshamrani, Adel Fahad Alrasheedi, Khalid Abdulaziz Alnowibet, and Ali Wagdy Mohamed. 2022. "An Amended Whale Optimization Algorithm for Optimal Bidding in Day Ahead Electricity Market" Axioms 11, no. 9: 456. https://doi.org/10.3390/axioms11090456

APA Style

Jain, K., Saxena, A., Alshamrani, A. M., Alrasheedi, A. F., Alnowibet, K. A., & Mohamed, A. W. (2022). An Amended Whale Optimization Algorithm for Optimal Bidding in Day Ahead Electricity Market. Axioms, 11(9), 456. https://doi.org/10.3390/axioms11090456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop