Next Article in Journal
Tooth Discoloration after Regenerative Endodontic Procedures with Calcium Silicate-Based Cements—An Ex Vivo Study
Next Article in Special Issue
A Spring Search Algorithm Applied to Engineering Optimization Problems
Previous Article in Journal
Climatic Indices over the Mediterranean Sea: A Review
Previous Article in Special Issue
Robust Parking Block Segmentation from a Surveillance Camera Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New “Doctor and Patient” Optimization Algorithm: An Application to Energy Commitment Problem

by
Mohammad Dehghani
1,
Mohammad Mardaneh
1,
Josep M. Guerrero
2,
Om Parkash Malik
3,
Ricardo A. Ramirez-Mendoza
4,*,
José Matas
5,
Juan C. Vasquez
2 and
Lizeth Parra-Arroyo
4
1
Department of Electrical and Electronics Engineering, Shiraz University of Technology, Shiraz 71557-13876, Iran
2
CROM Center for Research on Microgrids, Department of Energy Technology, Aalborg University, 9220 Aalborg, Denmark
3
Department of Electrical and Computer Engineering, University of Calgary, Calgary, AB T2N 1N4, Canada
4
School of Engineering and Sciences, Tecnologico de Monterrey, 64849 Monterrey, Mexico
5
Electric Engineering Department, Polytechnic University of Catalonia (EEBE-UPC), 08019 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(17), 5791; https://doi.org/10.3390/app10175791
Submission received: 24 July 2020 / Revised: 10 August 2020 / Accepted: 12 August 2020 / Published: 21 August 2020

Abstract

:
Regular assessments of events taking place around the globe can be a conduit for the development of new ideas, contributing to the research world. In this study, the authors present a new optimization algorithm named doctor and patient optimization (DPO). DPO is designed by simulating the process of treating patients by a physician. The treatment process has three phases, including vaccination, drug administration, and surgery. The efficiency of the proposed algorithm in solving optimization problems compared to eight other optimization algorithms on a benchmark standard test function with 23 objective functions is been evaluated. The results obtained from this comparison indicate the superiority and quality of DPO in solving optimization problems in various sciences. The proposed algorithm is successfully applied to solve the energy commitment problem for a power system supplied by a multiple energy carriers system.

1. Introduction

1.1. Motivation

Energy commitment (EC), the concept of choosing an adequate energy carrier operation, poses an important challenge in energy studies. Primary energy carriers are those that are extracted directly from natural resources, such as coal, oil, and natural gas, while secondary energy carriers are derived from primary energy [1]. In order to keep with the network’s energy demand, energy carriers are optimized considering the technical and economical constraints [2]. In fact, EC is a constrained optimization problem that can be solved using optimization algorithms [3].
Optimization algorithms perform well in solving a variety of problems. In order to achieve the appropriate pattern of utilization of energy carriers, the EC problem was assessed using suitable optimization tools.

1.2. Contribution

This paper proposes a new optimization algorithm named Doctor and Patient Optimizer (DPO) that obtains the optimal solution to an EC problem in power systems. The study aimed to achieve the following:
  • Design and present a novel optimization algorithm named “Doctor and Patient” Optimization.
  • Evaluate the proposed DPO algorithm on a set of benchmark test functions with 23 objective functions.
  • Compare the efficiency of the DPO to eight other optimization algorithms.
  • Study the EC issue on a standard energy grid with twenty-six power plants in different sectors of energy consumption (commercial, transportation, industrial, agriculture, residential, and public).
  • Apply DPO to EC problem solving.
  • Investigate the export and import of energy carriers in the EC problem.
  • Investigate oil refining in the EC problem.
  • Determine the appropriate pattern of energy carrier use to supply energy demand.

1.3. Paper Structure

The rest of paper is organized as follows. Section 2 reviews the studies conducted by the researchers. Section 3 introduces doctor and patient optimization, followed by the formulation of the energy commitment problem in Section 4. The benchmarking of DPO on twenty-three test functions and simulation of applying the proposed method on the EC problem is presented in Section 5, and, finally, conclusions are given in Section 6.

2. Background

Several research papers are published using different classical optimization algorithms to handle the optimization problem. The classical methods, such as the Lagrangian approach [4] Dynamic Programming (DP) [5] and Quadratic Programming (QP) [6], fail to optimize problems globally, which has led to the development of multiple new alternatives. Many heuristic and meta-heuristic optimization algorithms inspired by nature were developed in the search for alternatives.
New optimizing techniques inspired by major activities of living beings offer a wide range of problem-solving possibilities. Some are based on life style, movement patterns, or activities, like hunting, searching for food, etc. This has resulted in the development of many methods, such as in Reference [7], where the strategy for grey wolf optimization (GWO) was formulated based on the hunting of grey wolfs. Lion optimization algorithms (LOA) [8] were proposed based on the simulation of the lion life style; ant colony optimization (ACO) [9] was proposed based on movement pattern of ants; and donkey theorem optimization (DTO) [10] was presented based on behavior of donkeys searching for food. In general, optimization algorithms can be divided into four categories as physics-based, swarm-based, evolutionary-based, and game-based algorithms.
Physics-based algorithms are developed based on phenomena and laws of physics [11]. The Spring search algorithms (SSA) [12] is a physics-based algorithm which simulates Hooke’s law. The Water cycle algorithms (WCA) [13] is proposed based on the natural event of the water follow cycle from rivers and streams into the sea. Gravitational search algorithms (GSA) [14] are based on gravitational force modeling between bodies. Some of the other algorithms that fall into this category are: simulated annealing (SA) [15], curved space optimization (CSO) [16], galaxy-based search algorithm (GbSA) [17], artificial chemical reaction optimization algorithms (ACROA) [18], central force optimization (CFO) [19], and small world optimization algorithms (SWOA) [20].
Swarm-based algorithms have been suggested based on collectives of living things. Particle swarm optimization (PSO) [21], derived from the bird group’s social behavior during migration, is a common swarm-based algorithm. Another optimization process is the grasshopper optimization algorithm (GOA) [22] which simulates the grasshopper behavior. Marine predators algorithms (MPA) [23] are based on the biological interaction between predator and prey in the ocean. Some of the other algorithms that fall into this category are: grey wolf optimization (GWO) [7], lion optimization algorithm (LOA) [8], ant colony optimization (ACO) [9], donkey theorem optimization (DTO) [10], cuckoo search (CS) [24], artificial bee colony (ABC) [25], ant lion optimizer (ALO) [26], whale optimization algorithm (WOA) [27], and bat inspired algorithm (BA) [28].
Evolutionary-based algorithms use biologically based processes, such as mutation, reproduction, selection, and recombination. Genetic algorithm (GA) [29] is the most famous type of algorithm in this category and is based on the theory of Darwinian evolution. Some other algorithms in this category are: evolution strategy (ES) [30], differential evolution (DE) [31], biogeography-based optimizer (BBO) [32], and genetic programming (GP) [33].
Game-based algorithms have introduced new optimization techniques by simulating rules of different games. The dice game optimizer (DGO) [34] is a game-based algorithm that has been proposed based on the rules governing the game of dice and the impact the players have on each other. Another algorithm in this category is the orientation search algorithm (OSA) [35] that has been inspired by the game of orientation in which players move in the direction of a referee. Shell game optimization (SGO) [36] is a game-based algorithm proposed which is based on a simulation of the rules of the shell game.
Energy commitment (EC) sets the best template for using energy carriers because the technical limitations are dealt with first and the economic challenges after. Adjusting energy carriers to the highest demand would be unnecessary and costly. Indeed, energy carriers should be used optimally, as the proper management of energy resources can save considerable money. First, the energy demand must be determined in the EC issue. Similar to the unit commitment (UC) problem, this energy demand could span 24 h. In the UC issue, the demand for electricity must be fulfilled with the appropriate unit combination for every hour of the study.
UC involves adjusting thermal generators in order to meet the projected demand and minimize the cost of system operation [37]. UC is accountable in the selection of the units which can be set to operate economically [38]. UC also contributes to the power calculation of each unit based on total demand [39]. In power systems, it is important to create a table of optimum generating units with minimum fuel and transaction costs corresponding to the load requirements [40]. In order to solve the UC problem, both intelligent and classical techniques have been proposed [41]. a mixed-integer linear programming (MILP) model to figure out the transmission-constrained direct current (DC)-based unit commitment (UC) problem using the generalized generation distribution factors (GGDF) for modeling the transmission network constraints is proposed in Reference [42]. Intelligent techniques are an important choice in the engineering field due to their ability to optimize multi-range local optimal points [43]. The memetic binary differential evolution algorithm (MDPE) has been proposed to solve a profit-based UC problem [44]. An uncertain UC problem study is suggested in the presence of energy storage systems using list-based genetic algorithm-priority [45]. Quantum binary particle swarm optimization (QBPSO) algorithms are proposed to reduce operation cost in the UC problem [46]. Other algorithms, such as the whale optimization algorithm (WOA) [47], gray wolf algorithm (GWO) [48], shuffled frog-leaping algorithm [49], improved genetic algorithm [50], and simulated annealing [51], have also been suggested to find the solution of UC problem. The various studies in operation of power systems, such as energy reservation review [52], energy storage systems [53], and the impact of renewable energy sources [54], are analyzed by researchers.

3. Doctor and Patient Optimization (DPO)

In this section, the Doctor and Patient Optimization (DPO) algorithm is introduced to solve optimization problems. DPO are designed using simulation of patients’ treatment steps. The proposed algorithm has three phases, including: (a) vaccination, (b) drug administration, and (c) surgery. This process is such that population is vaccinated first to prevent infection. In the second phase, appropriate medication is prescribed to treat patients. Finally, in the third phase, surgery is performed on patients with a serious condition.

3.1. Mathematical Modeling

The population in DPO are patients who need to be treated by a doctor. This population of patients is specified in Equation (1).
P = [ P 1 P i P N         |         p 1 1 p 1 m p i 1 p i d p i m p N 1 p N m ] .
where P is the patients population, P i is the ith patient, p i d is the dth feature of ith patient, N is the number of patient (population), and m is the number of variables.
This population is treated and updated in three phases. The required information in this process is calculated by Equations (2)–(5).
d o s a g e i = 2 F i n F b e s t n ,
F i n = f i t i f w o r s t j = 1 N ( f i t j f w o r s t ) ,
f w o r s t = max ( f i t )   &   P w o r s t = P ( l o c a t i o n ( f w o r s t ) ) ,
f b e s t = min ( f i t )   &   P b e s t = P ( l o c a t i o n ( f b e s t ) ) .
Here, d o s a g e i is the dosage of vaccine or drug for ith patient, F i n is the normalized fitness of ith patient, F b e s t n is the normalized fitness of best patient, f w o r s t is the fitness function of worst patient,   f b e s t is the fitness function of best patient,   P w o r s t is the position of worst patient, and P b e s t is the position of best patient.

3.1.1. Phase A: Vaccination

An important step in the community health process is vaccination. This phase is simulated by Equations (6)–(7).
V i d = r a n d × ( d o s a g e i × p i d p w o r s t d ) ,
V i d = r a n d × ( d o s a g e i × p i d p w o r s t d ) .
Here, V i d is the dth dimension of vaccine for ith patient, r a n d is a random number in the interval [ 0 1 ] , and p w o r s t d is the dth dimension of worst patient.

3.1.2. Phase B: Drug Administration

In this phase of the patient treatment process, the doctor prescribes each patient pharmaceuticals according to the patient’s condition. Drug administration is simulated by Equations (8)–(9).
d i d = r a n d × ( p b e s t d   d o s a g e i × p i d ) ,
P i = { P i + d i ,     f i t ( P i + d i ) f i t i P i ,     e l s e .
Here, d i d is the dth dimension of a drug for the ith patient, and p b e s t d is the dth dimension of best patient.

3.1.3. Phase C: Surgery

Vaccination and medication are not enough for patients with serious conditions. In such cases, the patient’s condition will improve with surgery. This phase of treatment is modeled by Equation (10).
P i = { 0.6 × P i + 0.4 × P b e s t ,     F b e s t n F i n 0.9 F b e s t n P i ,     e l s e .

3.2. Implementation of DPO

After designing the proposed DPO algorithm, it can be used to solve optimization problems. Implementation of DPO is expressed in Algorithm 1.
Algorithm 1. The pseudo code of DPO
Start DPO
1System tuning and parameters determination.
2Formation of the initial population of patients: P.
3For iteration = 1: iteration max
4Fitness function evaluation.
5Updating f w o r s t and P w o r s t based (4).
6Updating f b e s t and P b e s t based (5).
7Updating F i n based (3).
8For i = 1:N
9Updating d o s a g e i based (2).
10Updating P i based phase a.
11Updating P i based phase b.
12Updating P i based phase c.
13End for i
14Saving f b e s t and   P b e s t .
15End for iteration
16Return best solution.
End DPO

4. Energy Commitment (EC) Problem

The EC analysis should be performed in a suitable area, such as the energy grid, which includes the public, commercial, residential, industrial, agricultural, and transportation sectors.
In the energy grid, the energy demand is determined as the sum of the demand in the various grid subdivisions using Equation (11).
E C f = E C 1 + E C 2 + + E C N =   i = 1 N E C i ,
where E C f is the total energy demand, E C i is the energy demand of the i-th sector of grid, and N is the number of different sectors of the energy grid.
In various sectors, the energy consumption is expressed in Equation (12):
E 1 =   [ E C 1   E C 2 E C i   E C N ] T .
Here, E 1 is the energy demand matrix in the various energy sectors.
Final energy consumption based on different energy carriers is determined by Equation (13):
E 2 =   T 1 , 2   ×   E 1 ,
where E 2 is the final energy consumption based on different energy carriers, and T 1 , 2 is the transform matrix of different energy sectors to different energy carriers.
Energy loss is modeled using Equation (14).
E 3 =   T 2 , 3   ×   E 2 .
Here, E 3 is the final energy consumption based on different energy carriers considering losses, and T 2 , 3 is the efficiency matrix.
Input fuels to generation unit in order to electrical energy demand supply are calculated by Equations (15)–(16).
E u =   T u   ×   E e ,
E e 1 =   T u , f   ×   E u ,
where E u is the value of generation of different units, T u is the separation matrix of electricity generated by different units that is specified by UC solving, E e is the total electrical energy demand, E e 1 is the input fuel to different units, and T u , f is the unit efficiency matrix.
The input of energy carriers to the units are calculated by Equation (17).
E e 2 =   T f , c   ×   E e 1 ,
where E e 2 is the value of energy carriers for electricity generation, and T f , c is the conversion matrix of input fuel to energy carriers.
In this stage after conversion of electrical energy demand to source energy carriers, final energy consumption is calculated using Equation (18).
E 4 =   E 3 + E e 2 E e .
E 4 is the final energy consumption after converting electrical energy demand to an input from energy carriers to units.
At this stage, Equation (19) is used to simulate the process of refining crude oil.
E p 1 =   T p   ×   E p .
Here, E p 1 represents the energy carriers produced by refining the oil, T p is the separation matrix of products produced from the refining process, and E p is the maximum capacity of refineries.
Final energy consumption considering the refining crude oil process is determined using Equation (20).
E 5 =   E 4 + E p E p 1 .
Here, E 5 is the final energy consumption after refining crude oil. Actually E 5 determines energy carriers to supply energy demand.
Finally, the import and export of energy carriers is determined using Equation (21).
E 6 =   E 5 P ,
where P is the domestic production of energy carriers, and E 6   is import and/or export of energy carriers. In E 6 , a negative sign denotes an export, while a positive sign means the import of energy carriers.

5. Simulation Study and Discussion

5.1. Case Study A: Benchmark Test Functions

In this section, the performance of DPO is evaluated on a standard set of benchmark test functions which have been used by the researchers in various earlier studies [55,56]. These benchmark functions includes twenty-three test functions that are categorized into Unimodal [57,58], Multimodal [58,59], and Fixed-dimension Multimodal [58] functions. The description of these test functions is found in Appendix A and in Table A1, Table A2 and Table A3.

5.1.1. Experimental Setup

The performance of the DPO is compared with the following eight optimization algorithms: Genetic Algorithm (GA) [60], Particle Swarm Optimization (PSO) [61], Gravitational Search Algorithm (GSA) [14], Teaching Learning Based Optimization (TLBO) [62], Grey Wolf Optimizer (GWO) [7], Grasshopper Optimization Algorithm (GOA) [22], Whale Optimization Algorithm (WOA) [27], and Marine Predators Algorithm (MPA) [23].
The proposed algorithm is implemented 30 times for each benchmark test function to obtain the average (avg), standard deviation (std), best, and worst values. In each run, the number of maximum iterations performed is fixed at 1000 for all the twenty-three benchmark test functions. The population size (N) is fixed at 50. The algorithm is implemented in MATLAB R2017b version using a 64-bit Core i7 processor with 3.20 GHz and 16 GB main memory.

5.1.2. Benchmarking Results of Unimodal Test Function

This group of functions is used to evaluate the exploitation ability of algorithms. The results of the implementation of the DPO and other mentioned algorithms on these test functions are presented in Table 1. DPO is clearly superior to all other compared algorithms in all F1 to F7 test functions.

5.1.3. Benchmarking Results of Multimodal Test Function

In this type of test functions, the number of local solutions are increased exponentially with the increasing dimensions of functions. As a result, it is very difficult to achieve the optimal response in this type of test functions. Table 2 shows the results of implementing and comparing the proposed algorithm and other eight optimization algorithms on this group of test functions, including F8 to F13.

5.1.4. Benchmarking Results of Fixed-Dimension Multimodal Test Function

The characteristic of this group of objective functions is the low number of local responses and dimensions. The results of the evaluation and optimization of these objective functions are given in Table 3. The ability of DPO to access the optimal answer is evident compared to other algorithms.

5.2. Case Study B: EC Problem

In this section, after implementing the DPO on benchmark test function and showing its strong ability in solving optimization problems, the proposed optimization algorithm is applied to the EC problem to determine the appropriate pattern of use of energy carriers.
The EC is implemented on an energy network with 26 power plants for a 24-h study period. The energy network included residential, commercial, public, industrial, transportation, and agriculture sectors and is supplied by various energy carriers. The energy demand in this network is shown for different sections in Table 4. The profile of this energy demand is displayed intuitively in Figure 1. All the other information surrounding the energy network is supplied in Appendix B and in Table A8, Table A9 and Table A10. The MBOE (millions of barrels of oil equivalent) unit is applied as the energy unit in this paper.

5.2.1. Objective Function and Constraints

In the present study, the objective function for solving the EC problem is considered to reduce the cost of supplying energy demand. This objective function for 24-h study period is expressed by Equation (22). Additionally, to optimize the EC’s objective function, the constraints related to the start-up cost of power plants and their authorized production range, specified in Equations (23)–(25), must be considered.
F o b j e c t i v e = m i n { t = 1 T [ i = 1 N c c a r r i e r i t × p r i c e i + i = 1 N g S C i t + i = 1 N g C i u i t ]   } ,
S C i t = { S C i ,     u i t > u i t 1 0 ,     e l s e ,
P g i m i n P g i P g i m a x ,
i = 1 N g P g i t = l o a d t .

5.2.2. DPO Implementation to EC Problem

The purpose of implementing the EC problem is to supply energy demand by determining the most appropriate use of energy carriers, considering technical and economical constraints. In the study of each hour of the 24-h period the first step, after the required energy conversions, was to determine possible combinations of power plants based on the required electrical energy demand. Therefore, all possible combinations of power plants are determined for each hour of the study period. The second step involved determining a suitable pattern for energy carrier use for the entire study period, as well as the optimal combination of power plants for each hour, based on the objective function and using the proposed optimization algorithm. This convenient pattern of energy carrier usage is actually the main output of EC problem.
The EC problem is coded in MATLAB and executed on a system with a quad-core 3.3 GHz processor and 8 GB of RAM. The pseudo code of EC problem solution using DPO is specified in Algorithm 2.
Algorithm 2. DPO implementation to EC problem
START
1:Problem information.
2:Inputs data: E 1 s t u d y   p e r i o d   , T 1 , 2
3:For Hour = 1: Study period (24 h)
4: E 1 = E 1 s t u d y   p e r i o d   ( Hour , : ) .
5: E 2 calculation based (13).
6: E 3 calculation based (14).
7: E e =   E 3 ( e d , 1 ) and e d = row number of electrical demand in E 3 .
8:END Hour
9:Determine possible combinations of power plants for electrical demand supplying.
10:DPO
11:Initial population formation based on possible combinations of units.
12:ITERATION = 1:T
13:For i = 1:Npopulatio
14:Combination = population (i,:).
15:IF this combination is possible.
16:UC Problem solving.
17:input energy to power plants calculation.
18:END UC solving.
19: E 4 calculation based (15) to (18).
20:Refinery simulation based (19).
21: E 5 calculation based (20).
22: E 6 calculation based (21).
23:Fitness calculation based (22).
24:Else if the combination is impossible.
25:Fitness = 1 × 10.
26:END if
27:END FOR
28:Updating f w o r s t and P w o r s t based (4).
29:Updating f b e s t and P b e s t based (5).
30:Updating based (3).
31:FOR i = 1:N
32:Updating d o s a g e i based (2).
33:Updating P i based phase a. (6) and (7).
34:Updating P i based phase b. (8) and (9).
35:Updating P i based phase c. (10).
36:END FOR
 
37:END ITERATION
 
38:EC outputs (for every hour and whole period of study).
39:Determining the pattern of energy carriers using.
40:Determining the UC output (power plant production).
41:Convergence curve.
42:Cost of energy supply.
43:Import and export of energy carriers.
END

5.2.3. Results and Discussion

The proposed DPO algorithm was implemented on the power system in order to achieve optimal results in an economical manner for the introduced energy commitment problem. The purpose of this operation is to reduce operating costs in order to supply energy demand. The important output of the energy commitment problem, the determination of the amount of different energy carriers for each hour of the study period, is specified in Table 5. The convergence curve (as an important indicator in the evaluation of optimization algorithms) of the implementation of the DPO on the EC problem is drawn in Figure 2. This curve shows the precise behavior of the algorithm while reaching the appropriate response, indicating the exploitation, exploration and power of the proposed algorithm. Another important output of the EC is to determine the appropriate pattern for the on and off state of power plant units for each hour of period of study to supply the electrical demand which is specified in Table 6. Additionally, the hourly production rate of the units during the study period, the output of the UC problem, is presented in Table 7. Finally, the import and export values of the energy carriers based on domestic production are specified in Table 8. According to this table, Petroleum (399,217), Fuel oil (29,054.3), Gas oil (478.824), Kerosene (297.468), and Coke gas (92.2559) are in the export section and liquid gas (2451.409), Gasoline (20,384.46), plane fuel (2934.952), natural gas (12,502.36), and coal (906.6509) are in the import section.

5.2.4. Comparison DPO and Other Algorithms on EC Problem

In order to evaluate the performance of the proposed algorithm in solving the EC, the other eight algorithms mentioned in this paper have been implemented on the EC problem. The results of this simulation are presented in Table 9. This table specifies the value of the objective function for each of the optimization algorithms. The proposed DPO algorithm is the best optimizer among the compared algorithms with the value of the objective function equal to 2.1153 × 107 Dollar. WOA with the value of the objective function 2.1739 × 107 Dollar, MPA with the value of the objective function 2.2365 × 107 Dollar, GWO with the value of the objective function 2.4257 × 107 Dollar, GOA with the value of the objective function 2.7592 × 107 Dollar, TLBO with the value of the objective function 3.2648 × 107 Dollar, GSA with the value of the objective function 6.7624 × 107 Dollar, PSO with The value of the target function is 5.2158 × 108 Dollar, and the GA with the value of the target function of 8.5146 × 108 Dollar are ranked second to ninth, respectively. Based on the results, the proposed algorithm has a high ability to solve the EC problem and is much more competitive than the other eight algorithms.

6. Conclusions

A new doctor and patient optimization (DPO) Algorithm was introduced based on a simulation of the patient treatment process. This treatment process has three phases including vaccination, drug administration, and surgery. To evaluate the effectiveness and performance of the DPO, two case studies were considered. In case study A, the performance and effectiveness of the proposed DPO algorithm was evaluated on a benchmark standard test function with twenty-three objective functions and compared to eight other algorithms. These results show the exploitation and exploration capacity of the proposed algorithm in solving optimization problems. In case study B, the proposed DPO algorithm was implemented on the energy commitment (EC) problem in a power system with twenty-six power plants and various energy sectors, including residential, commercial, public, industrial, transportation, and agriculture sectors. The purpose of the EC was to determine the appropriate pattern of use of energy carriers to supply energy demand and minimize operation costs considering the technical constraints. The DPO with high exploitation and exploration capacity was well implemented on the EC problem, and its results were determined including the appropriate pattern of use of energy carriers, proper composition, and production of power plants, as well as the amount of import and export of energy carriers.
In future works, the authors propose several study ideas, such as solving the EC problem using other optimization algorithms and techniques, creating a binary variant of the DPO which has an important potential contribution, and applying DPO to overcome many-objective real-life optimization problems, as well as multi-objective problems.

Author Contributions

Conceptualization, M.D., M.M., J.C.V. and J.M.G.; methodology, M.D.; software, M.D.; validation, R.A.R.-M., J.M.G., L.P.-A., and O.P.M.; formal analysis, J.M., R.A.R.-M., L.P.A., and O.P.M.; investigation, M.D., J.C.V. and O.P.M.; resources, R.A.R.-M.; data curation, R.A.R.-M.; writing—original draft preparation, M.D. and M.M.; writing—review and editing, J.M., J.C.V., O.P.M., L.P.-A., and R.A.R.-M.; visualization, M.D.; supervision, M.M.; project administration, M.D. and M.M.; funding acquisition, R.A.R.-M. All authors have read and agreed to the published version of the manuscript.

Funding

The current project was funded by Tecnologico de Monterrey and FEMSA Foundation (grant CAMPUSCITY project).

Conflicts of Interest

The authors declare no conflict of interest. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Table A1. Unit information.
Table A1. Unit information.
Row.Power PlantCapacity of Unit (MW)EfficiencyConstant CostPriorityMUT (Hour)MDT (Hour)Cold StartInitial ConditionsHot Start (Dollar)Cold Start (Dollar)
MinMax
1Thermal1004000.36831218−54108001500
2Thermal1004000.34531028−54107751500
3Combined Cycle1403500.45517738−54107251200
4Thermal68.951970.31726045−4287501300
5Gas68.951970.326055−4287001100
6Combined Cycle68.951970.4726065−428650950
7Thermal54.251550.3514375−328600850
8Gas54.251550.2514385−328550900
9Combined Cycle54.251550.514395−32−8500700
10Thermal54.251550.358143105−32−8450800
11Thermal251000.32218114−21−8200400
12Gas251000.27218124−21−8600900
13Combined Cycle251000.25218134−21−8250500
14Gas15.2760.381143−21−8400600
15Combined Cycle15.2760.381153−21−8250400
16Thermal15.2760.2981163−21−8400600
17Thermal15.2760.2981173−21−8300500
18Thermal4200.29118181−10−4300450
19Combined Cycle4200.291118191−10−4200350
20Gas4200.275118201−10−4200400
21Gas4200.27118211−10−1150300
22Thermal2.4120.2624221−10−350200
23Thermal2.4120.2524231−10−2100250
24Combined Cycle2.4120.2324241−10−1150300
25Combined Cycle2.4120.2224251−10−2100200
26Gas2.4120.224261−10−3150250
Table A2. T1,2 matrix.
Table A2. T1,2 matrix.
Residential, Commercial and PublicIndustrialTransportationAgricultureOtherNon-Energy
Petroleum000000
Liquid gas0.0510.0130.01000
Fuel oil0.0230.2120.014000
Gas oil0.0550.0870.3630.68900
Kerosene0.1410.00200.01800
Gasoline0.0020.0020.5730.00300
Plane fuel000.031000
Other products000000.402
Natural gas0.5640.5210.007000.497
Coke gas00.0210000
Coal0.000300000.101
Non-Commercial fuels0.06400000
Electricity(power)0.1020.1420.00040.2910
Table A3. Matrix Tp.
Table A3. Matrix Tp.
Petroleum0
liquid Gas0.032
Fuel Oil0.293
Gas Oil0.293
Kerosene0.099
Gasoline0.157
plane Fuel0
Other Products0.058
Natural Gas0
Coke Gas0
Coal0
Non-Commercial Fuels0
Electricity(power)0
Table A4. Conversion matrix input energy to fuel power plants.
Table A4. Conversion matrix input energy to fuel power plants.
Power PlantThermal UnitCombined Cycle UnitGas Unit
Fuel Oil0.25400
Gas Oil0.0030.0820.166
Natural Gas0.7430.9180.834
Table A5. Domestic supplies of energy carriers.
Table A5. Domestic supplies of energy carriers.
RowEnergy CarrierEnergy (Boe)
1Petroleum25,747.64405
2liquid Gas0
3Fuel Oil0
4Gas Oil0
5Kerosene0
6Gasoline0
7Plane Fuel0
8Other Products0
9Natural Gas9861.294929
10Coke Gas65.15249127
11Coal97.72873691
12Non-Commercial Fuels394.0174472
13Electricity(power)0
Table A6. Heating value [63] and energy rates [64].
Table A6. Heating value [63] and energy rates [64].
Energy CarrierHeating ValueEnergy Rates
Petroleum38.5 M J L i t 48 dollar/boe
Liquid Gas46.15 M J K g 374 dollar/tone
Fuel Oil42.18 M J K g 180 dollar/tone
Gas Oil43.38 M J K g 350 dollar/tone
Kerosene43.32 M J K g 500 dollar/tone
Gasoline44.75 M J K g 450 dollar/tone
Plane Fuel45.03 M J K g 555 dollar/tone
Natural Gas39 M J m 3 237 dollar/1e3m3
Coke Gas16.9 M J K g 157 dollar/tone
Coal26.75 M J K g 61 dollar/tone
Table A7. Matrix T23.
Table A7. Matrix T23.
Petroleum1000000000000
Liquid Gas0100000000000
Fuel Oil0010000000000
Gas Oil0001000000000
Kerosene0000100000000
Gasoline0000010000000
plane Fuel0000001000000
Other Products0000000100000
Natural Gas000000001.16010000
Coke Gas0000000001000
Coal0000000000100
Non-Commercial Fuels0000000000010
Electricity(power)0000000000001.3158

Appendix B

Table A8. Unimodal test functions.
Table A8. Unimodal test functions.
[ 100 , 100 ] m F 1 ( x ) = i = 1 m x i 2
[ 10 , 10 ] m F 2 ( x ) = i = 1 m | x i | + i = 1 m | x i |
[ 100 , 100 ] m F 3 ( x ) = i = 1 m ( j = 1 i x i ) 2
[ 100 , 100 ] m F 4 ( x ) = max { | x i | , 1 i m }
[ 100 , 100 ] m F 5 ( x ) = i = 1 m 1 [ 100 ( x i + 1 x i 2 ) 2 + ( x i 1 ) 2 ) ]
[ 100 , 100 ] m F 6 ( x ) = i = 1 m ( [ x i + 0.5 ] ) 2
[ 1.28 , 1.28 ] m F 7 ( x ) = i = 1 m i x i 4 + r a n d o m ( 0 , 1 )
Table A9. Multimodal test functions.
Table A9. Multimodal test functions.
[ 500 , 500 ] m F 8 ( x ) = i = 1 m x i   sin ( | x i | )
[ 5.12 , 5.12 ] m F 9 ( x ) = i = 1 m [   x i 2 10 cos ( 2 π x i ) + 10 ]
[ 32 , 32 ] m F 10 ( x ) = 20 exp ( 0.2 1 m i = 1 m x i 2 ) exp ( 1 m i = 1 m cos ( 2 π x i ) ) + 20 + e
[ 600 , 600 ] m F 11 ( x ) = 1 4000 i = 1 m x i 2 i = 1 m cos ( x i i ) + 1
[ 50 , 50 ] m F 12 ( x ) = π m   { 10 sin ( π y 1 ) + i = 1 m ( y i 1 ) 2 [ 1 + 10 sin 2 ( π y i + 1 ) ] + ( y n 1 ) 2 } + i = 1 m u ( x i , 10 , 100 , 4 )
u ( x i , a , i , n ) = { k ( x i a ) n                               x i > a 0                                       a < x i < a k ( x i a ) n                       x i < a
[ 50 , 50 ] m F 13 ( x ) = 0.1 {   sin 2 ( 3 π x 1 ) + i = 1 m ( x i 1 ) 2 [ 1 + sin 2 ( 3 π x i + 1 ) ] + ( x n 1 ) 2 [ 1 + sin 2 ( 2 π x m ) ] } + i = 1 m u ( x i , 5 , 100 , 4 )
Table A10. Multimodal test functions with fixed dimension.
Table A10. Multimodal test functions with fixed dimension.
[ 65.53 , 65.53 ] 2 . F 14 ( x ) = ( 1 500 + j = 1 25 1 j + i = 1 2 ( x i a i j ) 6 ) 1
[ 5 , 5 ] 4 F 15 ( x ) = i = 1 11 [ a i x 1 ( b i 2 + b i x 2 ) b i 2 + b i x 3 + x 4 ] 2
[ 5 , 5 ] 2 F 16 ( x ) = 4 x 1 2 2.1 x 1 4 + 1 3 x 1 6 + x 1 x 2 4 x 2 2 + 4 x 2 4
[−5,10]   ×   [0,15] F 17 ( x ) = ( x 2 5.1 4 π 2 x 1 2 + 5 π x 1 6 ) 2 + 10 ( 1 1 8 π ) cos x 1 + 10
[ 5 , 5 ] 2 F 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 ) ]
[ 0 , 1 ] 3 F 19 ( x ) = i = 1 4 c i exp ( j = 1 3 a i j ( x j P i j ) 2 )
[ 0 , 1 ] 6 F 20 ( x ) = i = 1 4 c i exp ( j = 1 6 a i j ( x j P i j ) 2 )
[ 0 , 10 ] 4 F 21 ( x ) = i = 1 5 [ ( X a i ) ( X a i ) T + 6 c i ] 1
[ 0 , 10 ] 4 F 22 ( x ) = i = 1 7 [ ( X a i ) ( X a i ) T + 6 c i ] 1
[ 0 , 10 ] 4 F 23 ( x ) = i = 1 10 [ ( X a i ) ( X a i ) T + 6 c i ] 1

References

  1. IEA. Energy Statistics Manual; OECD Publishing: Paris, France, 2004. [Google Scholar]
  2. Dehghani, M.; Montazeri, Z.; Ehsanifar, A.; Seifi, A.R.; Ebadi, J.M.; Grechko, O. Planning of energy carriers based on final energy consumption using dynamic programming and particle swarm optimization. Elec. Eng. Electromech. 2018, 5, 62–71. [Google Scholar] [CrossRef] [Green Version]
  3. Dehghani, M.; Montazeri, Z.; Malik, O. Energy commitment: A planning of energy carrier based on energy consumption. Электрoтехника и Электрoмеханика 2019, 6. [Google Scholar] [CrossRef]
  4. Smith, D.K.; Bertsekas, D.P. Nonlinear Programming. J. Oper. Res. Soc. 1997, 48, 334. [Google Scholar] [CrossRef]
  5. Bellman, R. Dynamic programming. Science 1966, 153, 34–37. [Google Scholar] [CrossRef] [PubMed]
  6. Frank, M.; Wolfe, P. An algorithm for quadratic programming. Nav. Res. Logist. Q. 1956, 3, 95–110. [Google Scholar] [CrossRef]
  7. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2013, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  8. Yazdani, M.; Jolai, F. Lion Optimization Algorithm (LOA): A nature-inspired metaheuristic algorithm. J. Comput. Des. Eng. 2015, 3, 24–36. [Google Scholar] [CrossRef] [Green Version]
  9. Dorigo, M.; Stützle, T. Ant Colony Optimization: Overview and Recent Advances. Handb. Metaheuristics. Int. Ser. Oper. Res. Manag. Sci. 2019, 272, 311–351. [Google Scholar]
  10. Dehghani, M.; Mardaneh, M.; Malik, O.P.; NouraeiPour, S.M. DTO: Donkey Theorem Optimization. In Proceedings of the 2019 27th Iranian Conference on Electrical Engineering (ICEE), Yazd, Iran, 30 April–2 May 2019; Institute of Electrical and Electronics Engineers: New Jersey, NJ, USA, 2019; pp. 1855–1859. [Google Scholar]
  11. Barry, J.; Thron, C. A Computational Physics-Based Algorithm for Target Coverage Problems. In Advances in Nature-Inspired Computing and Applications; Shandilya, S., Nagar, A., Eds.; Springer: Cham, Germany, 2019; pp. 269–290. [Google Scholar]
  12. Dehghani, M.; Montazeri, Z.; Dehghani, A.; Seifi, A. Spring search algorithm: A new meta-heuristic optimization algorithm inspired by Hooke’s law. In Proceedings of the 4th International Conference on Knowledge-Based Engineering and Innovation (KBEI), Berlin, Germany, 21–23 December 2017; pp. 0210–0214. [Google Scholar]
  13. Eskandar, H.; Sadollah, A.; Bahreininejad, A.; Hamdi, M. Water cycle algorithm–A novel metaheuristic optimization method for solving constrained engineering optimization problems. Comput. Struct. 2012, 110, 151–166. [Google Scholar] [CrossRef]
  14. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A Gravitational Search Algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  15. Kirkpatrick, S.; Gelatt, J.C.D.; Vecchi, M.P. Optimization by Simulated Annealing. World Sci. Lect. Notes Phys. 1986, 220, 339–348. [Google Scholar] [CrossRef]
  16. Moghaddam, F.F.; Moghaddam, R.F.; Cheriet, M. Curved space optimization: A random search based on general relativity theory. arXiv 2012, arXiv:1208.2214. [Google Scholar]
  17. Shah-Hosseini, H. Principal components analysis by the galaxy-based search algorithm: A novel metaheuristic for continuous optimisation. Int. J. Comput. Sci. Eng. 2011, 6, 132. [Google Scholar] [CrossRef]
  18. Alatas, B. ACROA: Artificial Chemical Reaction Optimization Algorithm for global optimization. Expert Syst. Appl. 2011, 38, 13170–13180. [Google Scholar] [CrossRef]
  19. Formato, R.A. Central Force Optimization: A New Nature Inspired Computational Framework for Multidimensional Search and Optimization. In Nature Inspired Cooperative Strategies for Optimization (NICSO 2007); Krasnogon, N., Nicosia, V., Pavone, M., Pelta, D.A., Eds.; Springer Science and Business Media LLC: Berlin, Germany, 2008; Volume 129, pp. 221–238. [Google Scholar]
  20. Du, H.; Wu, X.; Zhuang, J. Small-World Optimization Algorithm for Function Optimization. In Proceedings of the Computer Vision, Xi’an, China, 24–28 September 2006; Springer Science and Business Media LLC: Berlin, Germany, 2006; pp. 264–273. [Google Scholar]
  21. Bansal, J.C. Particle Swarm Optimization. In Evolutionary and Swarm Intelligence Algorithms; Springer Science and Business Media LLC: Berlin, Germany, 2018; pp. 11–23. [Google Scholar]
  22. Saremi, S.; Mirjalili, S.; Lewis, A. Grasshopper Optimisation Algorithm: Theory and application. Adv. Eng. Softw. 2017, 105, 30–47. [Google Scholar] [CrossRef]
  23. Faramarzi, A.; Heidarinejad, M.; Mirjalili, S.; Gandomi, A.H. Marine Predators Algorithm: A nature-inspired metaheuristic. Expert Syst. Appl. 2020, 152, 113377. [Google Scholar] [CrossRef]
  24. Gandomi, A.H.; Yang, X.-S.; Alavi, A.H. Cuckoo search algorithm: A metaheuristic approach to solve structural optimization problems. Eng. Comput. 2011, 29, 17–35. [Google Scholar] [CrossRef]
  25. Karaboga, D.; Basturk, B. Artificial Bee Colony (ABC) Optimization Algorithm for Solving Constrained Optimization Problems. Comput. Vis. 2007, 4529, 789–798. [Google Scholar] [CrossRef]
  26. Mirjalili, S. The Ant Lion Optimizer. Adv. Eng. Softw. 2015, 83, 80–98. [Google Scholar] [CrossRef]
  27. Mirjalili, S.; Lewis, A. The Whale Optimization Algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  28. Yang, X.-S. A New Metaheuristic Bat-Inspired Algorithm. In Nature Inspired Cooperative Strategies for Optimization (NICSO 2010); Krasnogon, N., Nicosia, V., Pavone, M., Pelta, D.A., Eds.; Springer Science and Business Media LLC: Berlin, Germany, 2010; Volume 284, pp. 65–74. [Google Scholar]
  29. Castillo, O.; Aguilar, L.T. Genetic Algorithms. In Type-2 Fuzzy Logic. in Control. of Nonsmooth Systems; Springer: Cham, Switzerland, 2019; pp. 23–39. [Google Scholar]
  30. Beyer, H.-G.; Schwefel, H.-P. Evolution strategies—A comprehensive introduction. Nat. Comput. 2002, 1, 3–52. [Google Scholar] [CrossRef]
  31. Storn, R.; Price, K. Differential Evolution-A Simple and Efficient Adaptive Scheme for Global Optimization Over Continuous Spaces; Berkeley: Berkeley, CA, USA, 1995. [Google Scholar]
  32. Mirjalili, S. BiogeographyB-Based Optimisation. In Evolutionary Algorithms and Neural Networks; Krasnogon, N., Nicosia, V., Pavone, M., Pelta, D.A., Eds.; Springer Science and Business Media LLC: Berlin, Germany, 2018; pp. 57–72. [Google Scholar]
  33. Koza, J.R. Genetic programming: A Paradigm for Genetically Breeding Populations of Computer Programs to Solve Problems; Stanford University, Department of Computer Science: Stanford, CA, USA, 1990. [Google Scholar]
  34. Dehghani, M.; Montazeri, Z.; Malik, O.P. DGO: Dice Game Optimizer. Gazi Univ. J. Sci. 2019, 32, 871–882. [Google Scholar] [CrossRef] [Green Version]
  35. Dehghani, M.; Montazeri, Z.; Malik, O.P.; Ehsanifar, A.; Dehghani, A. 0OSA: Orientation Search Algorithm. Int. J. Ind. Elect. Control. Optim. 2019, 2, 99–112. [Google Scholar]
  36. Dehghani, M.; Shiraz University of Technology; Montazeri, Z.; Malik, O.; Givi, H.; Guerrero, J. University of Calgary; University of Shahreza; Aalborg University Shell Game Optimization: A Novel Game-Based Algorithm. Int. J. Intell. Eng. Syst. 2020, 13, 10. [Google Scholar] [CrossRef]
  37. Wood, A.J.; Wollenberg, B.F. Power Generation, Operation, and Control; John Wiley & Sons: New Jersey, NJ, USA, 2012. [Google Scholar]
  38. Abdou, I.; Tkiouat, M. Unit Commitment Problem in Electrical Power System: A Literature Review. Int. J. Electr. Comput. Eng. 2018, 8, 1357–1372. [Google Scholar] [CrossRef]
  39. Gögler, P.; Dorfner, M.; Hamacher, T. Hybrid Robust/Stochastic Unit Commitment With Iterative Partitions of the Continuous Uncertainty Set. Front. Energy Res. 2018, 6, 71. [Google Scholar] [CrossRef]
  40. Tiwari, S.; Dwivedi, B.; Dave, M. A two stage solution methodology for deterministic unit commitment problem. In Proceedings of the 2016 IEEE Uttar Pradesh Section International Conference on Electrical, Computer and Electronics Engineering (UPCON), Varanasi, India, 9–11 December 2016; pp. 317–322. [Google Scholar]
  41. Del Nozal, A.R.; Tapia, A.; Alvarado-Barrios, L.; Reina, D.G. Application of Genetic Algorithms for Unit Commitment and Economic Dispatch Problems in Microgrids. In Nature Inspired Computing for Data Science; Krasnogon, N., Nicosia, V., Pavone, M., Pelta, D.A., Eds.; Springer Science and Business Media LLC: Berlin, Germany, 2019; pp. 139–167. [Google Scholar]
  42. Gutiérrez-Alcaraz, G.; Hinojosa, V. Using Generalized Generation Distribution Factors in a MILP Model to Solve the Transmission-Constrained Unit Commitment Problem. Energies 2018, 11, 2232. [Google Scholar] [CrossRef] [Green Version]
  43. Hussein, B.M.; Jaber, A.S. Unit commitment based on modified firefly algorithm. Meas. Control. 2020, 53, 320–327. [Google Scholar] [CrossRef] [Green Version]
  44. Dhaliwal, J.S.; Dhillon, J.S. Profit based unit commitment using memetic binary differential evolution algorithm. Appl. Soft Comput. 2019, 81, 105502. [Google Scholar] [CrossRef]
  45. Nikzad, H.R.; Abdi, H. A robust unit commitment based on GA-PL strategy by applying TOAT and considering emission costs and energy storage systems. Electr. Power Syst. Res. 2020, 180, 106154. [Google Scholar] [CrossRef]
  46. Hussain, A.N.; Ismail, A.A. Operation cost reduction in unit commitment problem using improved quantum binary PSO algorithm. Int. J. Electr. Comput. Eng. 2020, 10, 1149–1155. [Google Scholar] [CrossRef]
  47. Strikanth, R.K.; Panwar, L.; Panigrahi, B.K.; Kumar, R. Binary whale optimization algorithm: A new metaheuristic approach for profit-based unit commitment problems in competitive electricity markets. Eng. Optim. 2018, 51, 369–389. [Google Scholar] [CrossRef]
  48. Panwar, L.K.; Reddy, S.; Verma, A.; Panigrahi, B.K.; Kumar, R. Binary Grey Wolf Optimizer for large scale unit commitment problem. Swarm Evol. Comput. 2018, 38, 251–266. [Google Scholar] [CrossRef]
  49. Ebrahimi, J.; Hosseinian, S.H.; Gharehpetian, G.B. Unit Commitment Problem Solution Using Shuffled Frog Leaping Algorithm. IEEE Trans. Power Syst. 2010, 26, 573–581. [Google Scholar] [CrossRef]
  50. Jo, K.-H.; Kim, M.-K. Improved Genetic Algorithm-Based Unit Commitment Considering Uncertainty Integration Method. Energies 2018, 11, 1387. [Google Scholar] [CrossRef] [Green Version]
  51. Simopoulos, D.; Kavatza, S.; Vournas, C. Reliability Constrained Unit Commitment Using Simulated Annealing. Ieee Trans. Power Syst. 2006, 21, 1699–1706. [Google Scholar] [CrossRef]
  52. Carrión, M.; Zarate-Minano, R.; Domínguez, R. A Practical Formulation for Ex-Ante Scheduling of Energy and Reserve in Renewable-Dominated Power Systems: Case Study of the Iberian Peninsula. Energies 2018, 11, 1939. [Google Scholar] [CrossRef] [Green Version]
  53. Li, J.; Niu, D.; Wu, M.; Wang, Y.; Li, F.; Dong, H. Research on Battery Energy Storage as Backup Power in the Operation Optimization of a Regional Integrated Energy System. Energies 2018, 11, 2990. [Google Scholar] [CrossRef] [Green Version]
  54. Dominković, D.F.; Stark, G.; Hodge, B.-M.; Pedersen, A.S. Integrated Energy Planning with a High Share of Variable Renewable Energy Sources for a Caribbean Island. Energies 2018, 11, 2193. [Google Scholar] [CrossRef] [Green Version]
  55. Kaur, S.; Awasthi, L.K.; Sangal, A.; Dhiman, G. Tunicate Swarm Algorithm: A new bio-inspired based metaheuristic paradigm for global optimization. Eng. Appl. Artif. Intell. 2020, 90, 103541. [Google Scholar] [CrossRef]
  56. Dhiman, G.; Garg, M.; Nagar, A.K.; Kumar, V.; Dehghani, M. A Novel Algorithm for Global Optimization: Rat Swarm Optimizer. J. Ambient Int. Human. Comp. 2020, 11, 1868–5145. [Google Scholar]
  57. Digalakis, J.; Margaritis, K.G. On benchmarking functions for genetic algorithms. Int. J. Comput. Math. 2001, 77, 481–506. [Google Scholar] [CrossRef]
  58. Wang, G.; Gandomi, A.H.; Yang, X.-S.; Alavi, A.H. A novel improved accelerated particle swarm optimization algorithm for global numerical optimization. Eng. Comput. 2014, 31, 1198–1220. [Google Scholar] [CrossRef]
  59. Yang, X.-S. Firefly algorithm, stochastic test functions and design optimisation. Int. J. Bio-Inspired Comput. 2010, 2, 78. [Google Scholar] [CrossRef]
  60. Mirjalili, S. Genetic Algorithm. In Evolutionary Algorithms and Neural Networks; Krasnogon, N., Nicosia, V., Pavone, M., Pelta, D.A., Eds.; Springer: Berlin, Germany, 2019; pp. 43–55. [Google Scholar]
  61. Mirjalili, S. Optimisation. In Evolutionary Algorithms and Neural Networks; Krasnogon, N., Nicosia, V., Pavone, M., Pelta, D.A., Eds.; Springer: Berlin, Germany, 2019; pp. 15–31. [Google Scholar]
  62. Rao, R.V.; Savsani, V.J.; Vakharia, D. Teaching–learning-based optimization: A novel method for constrained mechanical design optimization problems. Comput. Des. 2011, 43, 303–315. [Google Scholar] [CrossRef]
  63. Birol, F. International Energy Agency. Global Energy Review Report; IEA: Paris, France, 2004; pp. 1–50. [Google Scholar]
  64. U.S. Energy Information Administration (EIA). Available online: http://www.eia.gov (accessed on 12 August 2020).
Figure 1. Energy demand profile.
Figure 1. Energy demand profile.
Applsci 10 05791 g001
Figure 2. Convergence curve of doctor and patient optimization (DPO) on energy commitment (EC) solving.
Figure 2. Convergence curve of doctor and patient optimization (DPO) on energy commitment (EC) solving.
Applsci 10 05791 g002
Table 1. Optimization results on unimodal test functions.
Table 1. Optimization results on unimodal test functions.
DPOMPAWOAGWOGOATLBOGSAPSOGA
03.27 × 10−211.41 × 10−306.59 × 10−282.81 × 10−13.55 × 10−21.16 × 10−164.98 × 10−91.95 × 10−12AveF1
04.61 × 10−214.91 × 10−306.34 × 10−51.11 × 10−11.06 × 10−16.10 × 10−171.40 × 10−82.01 × 10−11std
5.20 × 10−1851.57 × 10−121.06 × 10−217.18 × 10−173.96 × 10−13.23 × 10−51.70 × 10−17.29 × 10−46.53 × 10−18AveF2
01.42 × 10−122.39 × 10−212.90 × 10−21.41 × 10−18.57 × 10−59.29 × 10−11.84 × 10−35.10 × 10−17std
1.13 × 10−1188.64 × 10−25.39 × 10−73.29 × 10−64.31 × 104.91 × 1034.16 × 1021.4 × 107.70 × 10−10AveF3
5.32 × 10−1181.444 × 10−12.93 × 10−67.918.973.89 × 1031.56 × 1027.137.36 × 10−9std
1.48 × 10−1522.60 × 10−87.25 × 10−28.73 × 10−18.80 × 10−11.87 × 101.126.00 × 10−19.17 × 10AveF4
09.25 × 10−93.97 × 10−11.19 × 10−12.50 × 10−18.219.89 × 10−11.72 × 10−15.67 × 10std
25.106144.6049 × 102.79 × 108.91 × 1021.18 × 1027.37 × 1023.85 × 104.93 × 105.57 × 102AveF5
1.43 × 10−144.22 × 10−17.63 × 10−12.97 × 1021.43 × 1021.98 × 1033.47 × 103.89 × 104.16 × 10std
03.98 × 10−13.118.18 × 10−173.15 × 10−14.881.08 × 10−169.23 × 10−93.15 × 10−1AveF6
01.91 × 10−15.32 × 10−11.70 × 10−189.98 × 10−29.75 × 10−14.00 × 10−171.78 × 10−89.98 × 10−2std
4.15 × 10−51.80 × 10−31.42 × 10−35.37 × 10−32.02 × 10−23.88 × 10−27.68 × 10−16.92 × 10−26.79 × 10−4AveF7
1.82 × 10−201.00 × 10−31.14 × 10−31.89 × 10−17.43 × 10−35.79 × 10−22.772.87 × 10−23.29 × 10−3std
Table 2. Optimization results on multimodal test functions.
Table 2. Optimization results on multimodal test functions.
DPOMPAWOAGWOGOATLBOGSAPSOGA
−8548.93−8.36 × 102−5.10 × 102−6.12 × 10−6.92 × 102−3.81 × 102−2.75 × 102−5.01 × 102−5.11 × 102AveF8
8.13 × 10−138.11 × 1026.95 × 1023.94 × 109.19 × 102.83 × 105.72 × 104.28 × 104.37 × 10std
0003.10 × 10−11.01 × 1022.23 × 103.35 × 101.20 × 10−11.23 × 10AveF9
0003.91 × 101.89 × 103.25 × 101.19 × 104.01 × 104.11 × 10std
4.44 × 10−159.69 × 10−127.401.06 × 10−131.151.55 × 108.25 × 10−95.20 × 10−115.31 × 10−11AveF10
7.06 × 10−316.13 × 10−129.894.34 × 10−27.87 × 10−18.111.90 × 10−91.08 × 10−101.11 × 10−10std
002.89 × 10−42.49 × 10−35.74 × 10−13.01 × 10−18.193.24 × 10−63.31 × 10−6AveF11
001.58 × 10−31.34 × 10−41.12 × 10−12.89 × 10−13.704.11 × 10−54.23 × 10−5std
1.35 × 10−38.50 × 10−33.39 × 10−11.34 × 10−21.275.21 × 102.65 × 10−18.93 × 10−89.16 × 10−8AveF12
9.31 × 10−185.20 × 10−32.14 × 10−16.23 × 10−21.022.47 × 1023.14 × 10−14.77 × 10−74.88 × 10−7std
7.44 × 10−19.90 × 10−11.896.54 × 10−16.60 × 10−22.81 × 1025.738.26 × 10−19.39 × 10−1AveF13
6.95 × 10−161.93 × 10−12.66 × 10−14.47 × 10−34.33 × 10−28.63 × 1028.954.39 × 10−24.49 × 10−2std
Table 3. Optimization results on multimodal test functions with low dimension.
Table 3. Optimization results on multimodal test functions with low dimension.
DPOMPAWOAGWOGOATLBOGSAPSOGA
9.98 × 10−19.98 × 10−12.11 × 101.26 × 109.98 × 106.79 × 103.61 × 102.77 × 104.39 × 10AveF14
1.02 × 10−152.47 × 10−132.49 × 106.86 × 10−19.14 × 10−11.12 × 102.96 × 102.32 × 104.41 × 10−2std
3.11 × 10−48.21 × 10−33.66 × 10−31.01 × 10−27.15 × 10−25.15 × 10−26.84 × 10−29.09 × 10−37.36 × 10−2AveF15
2.42 × 10−194.09 × 10−157.60 × 10−23.75 × 10−31.26 × 10−13.45 × 10−37.37 × 10−22.38 × 10−32.39 × 10−3std
−1.03 × 10−1.02 × 10−1.02 × 10−1.02 × 10−1.02 × 10−1.01 × 10−1.02 × 10−1.02 × 10−1.02 × 10AveF16
3.97 × 10−164.46 × 10−167.02 × 10−93.23 × 10−54.74 × 10−83.64 × 10−80.00 × 100.00 × 104.19 × 10−7std
3.98 × 10−13.98 × 10−13.98 × 10−13.98 × 10−13.98 × 10−13.98 × 10−13.98 × 10−13.98 × 10−13.98 × 10−1AveF17
9.93 × 10−179.12 × 10−157.00 × 10−57.61 × 10−41.15 × 10−79.45 × 10−151.13 × 10−169.03 × 10−163.71 × 10−17std
3.00 × 103.00 × 103.00 × 103.00 × 103.00 × 103.00 × 103.00 × 103.00 × 103.00 × 10AveF18
8.94 × 10−161.95 × 10−157.16 × 10−62.25 × 10−51.48 × 101.94 × 10−103.24 × 10−26.59 × 10−56.33 × 10−7std
−3.86 × 10−3.86 × 10−3.84 × 10−3.75 × 10−3.77 × 10−3.73 × 10−3.86 × 10−3.80 × 10−3.81 × 10AveF19
2.68 × 10−152.42 × 10−71.57 × 10−32.55 × 10−33.53 × 10−79.69 × 10−44.15 × 10−13.37 × 10−154.37 × 10−10std
−3.32 × 10−3.32 × 10−2.98 × 10−2.84 × 10−3.23 × 10−2.17 × 10−1.47 × 10−3.32 × 10−2.39 × 10AveF20
1.29 × 10−151.14 × 10−113.76 × 10−13.71 × 10−15.37 × 10−21.64 × 10−15.32 × 10−12.66 × 10−14.37 × 10−1std
−10.15 × 10−8.11 × 10−7.05 × 10−2.28 × 10−7.38 × 10−7.33 × 10−4.57 × 10−7.54 × 10−5.19 × 10AveF21
4.57 × 10−152.53 × 10−113.62 × 101.80 × 102.91 × 101.29 × 101.30 × 102.77 × 102.34 × 10std
−1.04 × 10−1.00 × 10−8.18 × 10−3.99 × 10−8.50 × 10−1.00 × 10−6.58 × 10−8.55 × 10−2.97 × 10AveF22
2.78 × 10−152.81 × 10−113.82 × 101.99 × 103.02 × 102.89 × 10−42.64 × 103.08 × 101.37 × 10−2std
−10.53 × 10−10.41 × 10−9.34 × 10−4.49 × 10−8.41 × 10−2.46 × 10−9.37 × 10−9.19 × 10−3.10 × 10AveF23
2.98 × 10−153.89 × 10−112.41 × 10−41.96 × 103.13 × 101.19 × 102.75 × 102.52 × 102.37 × 10std
Table 4. Final energy consumption (barrels of oil equivalent (BOE)).
Table 4. Final energy consumption (barrels of oil equivalent (BOE)).
Hour 1 2 3 4 5 6 7 8
Residential, Commercial, and Public4609.3734690.7154582.2594609.3734744.9435016.0825422.7926588.692
Industrial2169.2522207.5332156.4922169.2522233.0542360.6572552.0623100.755
Transportation2931.1422982.8682913.92931.1423017.3523189.7723448.4024189.808
Agriculture384.9789391.7726382.7143384.9789396.3018418.9476452.9163550.2934
Other28.8157929.324328.6462828.8157929.6633131.3583533.9009241.18962
Non-Energy983.19461000.545977.4111983.19461012.1121069.9471156.71405.39
Hour 9 10 11 12 13 14 15 16
Residential, Commercial, and Public6886.9467049.6297239.4277022.5157022.5156914.067103.8577185.199
Industrial3241.1183317.683407.0023304.923304.923253.8793343.2013381.482
Transportation4379.4714482.9234603.6174465.6814465.6814396.7134517.4074569.133
Agriculture575.2038588.7913604.6433586.5267586.5267577.4683593.3204600.1142
Other43.0541744.071245.2577343.901743.901743.2236844.4102144.91872
Non-Energy1469.0081503.7091544.1941497.9261497.9261474.7921515.2761532.627
Hour 17 18 19 20 21 22 23 24
Residential, Commercial, and Public6914.066859.8326778.496914.067049.6296724.2625965.0714988.968
Industrial3253.8793228.3583190.0773253.8793317.683164.5562807.2682347.897
Transportation4396.7134362.2294310.5034396.7134482.9234276.0183793.2423172.53
Agriculture577.4683572.9392566.1454577.4683588.7913561.6163498.208416.683
Other43.2236842.8846742.3761643.2236844.071242.0371537.2910231.18885
Non-Energy1474.7921463.2251445.8741474.7921503.7091434.3071272.3691064.164
Table 5. The need of energy carriers (BOE).
Table 5. The need of energy carriers (BOE).
Hour 1 2 3 4 5 6 7 8
Liquid Gas292.5897297.7531290.8686292.5897301.1953318.4065344.2232418.2312
Fuel Oil1107.2231133.4251098.5651114.3231132.1741206.2421284.841520.263
Gas Oil1931.7911963.2361921.2921930.4191997.8482110.832310.592837.866
Kerosene661.1897672.8578657.3004661.1897680.6365719.53777.8703945.1124
Gasoline1694.2561724.1551684.291694.2561744.0871843.751993.2432421.79
Plane Fuel90.8653992.468990.3308990.8653993.537998.88293106.9005129.8841
Natural Gas7559.8587691.9547515.8447564.337822.2838285.0979021.85111,060.32
Coke Gas45.554346.358245.2863345.554346.8941349.573853.5932965.11585
Coal100.6855102.4623100.0932100.6855103.6468109.5695118.4535143.921
Hour 9 10 1112 13 14 15 16
Liquid Gas437.1635447.4902459.538445.7691445.7691438.8846450.9324456.0957
Fuel Oil1603.1731653.5361693.2151651.4061649.9661612.9751662.1161685.183
Gas Oil2960.8473026.3163117.6533013.4923013.892971.3413054.7123091.246
Kerosene987.89521011.2311038.4571007.3421007.342991.78461019.011030.678
Gasoline2531.4182591.2162660.9792581.2492581.2492541.3852611.1482641.047
Plane Fuel135.7636138.9706142.7121138.4361138.4361136.2981140.0396141.6431
Natural Gas11,593.9311,881.3112249.311,829.511,829.9411,641.0911982.2712134.19
Coke Gas68.0634869.6712871.5470569.4033169.4033168.3314570.2072171.01111
Coal150.4359153.9895158.1354153.3973153.3973151.0282155.1741156.9509
Hour 17 18 19 20 21 22 23 24
Liquid Gas438.8846435.4424430.279438.8846447.4902426.8368378.6455316.6853
Fuel Oil1614.7721593.3711563.7591619.0751654.9761550.4611418.2041209.212
Gas Oil2969.8522950.3522918.8682967.7273025.9182894.2682530.5522097.053
Kerosene991.7846984.0059972.3378991.78461011.231964.5591855.6573715.6406
Gasoline2541.3852521.4522491.5542541.3852591.2162471.6212192.5671833.783
Plane Fuel136.2981135.2291133.6256136.2981138.9706132.5566117.590598.34843
Natural Gas11,638.7511,546.7711,404.711,636.1411,880.8711,306.049885.8698211.221
Coke Gas68.3314567.7955266.9916268.3314569.6712866.4556858.9526249.30583
Coal151.0282149.8437148.0669151.0282153.9895146.8823130.2988108.9772
Table 6. Appropriate combination of units and total cost for energy supply.
Table 6. Appropriate combination of units and total cost for energy supply.
Hour123456789101112Cost ($)
Combination101010111317192223222622
Hour1314151617181920212223242.1153 × 107
Combination212322242023221822171512
Table 7. Unit commitment (UC) result (MW).
Table 7. Unit commitment (UC) result (MW).
Hourunit 1Unit 2unit 3unit 4unit 5unit 6unit 7unit 8unit 9unit 10unit 11unit 12unit 13
140040035019719719765.6652654.2554.2554.25000
2400400350197197197100.419654.2554.2554.25000
3400400350197197196.830554.2554.2554.2554.25000
4400400350197197183.415354.2554.2554.2554.252500
5400400350197197191.339254.2554.2554.2554.25252525
6400400350197197197103.637254.2554.2554.25252525
740040035019719719715515567.9091354.25252525
8400400350197197197155155155155100100100
9400400350197197197155155155155100100100
10400400350197197197155155155155100100100
11400400350197197197155155155155100100100
12400400350197197197155155155155100100100
13400400350197197197155155155155100100100
14400400350197197197155155155155100100100
15400400350197197197155155155155100100100
16400400350197197197155155155155100100100
17400400350197197197155155155155100100100
18400400350197197197155155155155100100100
19400400350197197197155155155155100100100
20400400350197197197155155155155100100100
21400400350197197197155155155155100100100
22400400350197197197155155155155100100100
2340040035019719719715515515515510032.2550525
2440040035019719719715577.102454.2554.2525250
Hourunit 14unit 15unit 16unit 17unit 18unit 19unit 20unit 21unit 22unit 23unit 24unit 25unit 26
10000000000000
20000000000000
30000000000000
40000000000000
50000000000000
615.215.215.215.2000000000
715.215.215.215.2440000000
87629.3053515.215.244442.40000
976767632.738144442.42.4000
10767676762019.046874400000
1176767676202020201212129.7404422.4
1276767676205.062077442.40000
1376767676207.4620774400000
1476767644.3228944442.42.4000
15767676762020207.8164642.40000
1676767676202020201210.570852.400
1776767653.12289444000000
1876767621.153344442.42.4000
19767649.5989215.244442.40000
2076767661.12289400000000
21767676762016.64687442.40000
22767644.8293215.2000000000
2315.215.200000000000
240000000000000
Table 8. Import and export of carriers (BOE).
Table 8. Import and export of carriers (BOE).
HourImportExport
Petroleum0−399,217
Liquid Gas2451.4090
Fuel Oil0−29,054.3
Gas Oil0−478.824
Kerosene0−297.468
Gasoline20,384.460
Plane Fuel29,34.9520
Natural Gas12,502.360
Coke Gas0−92.2559
Coal906.65090
Table 9. Results for DPO and other algorithms in EC problem.
Table 9. Results for DPO and other algorithms in EC problem.
AlgorithmAvg (Dollar)Std (Dollar)Rank
GA8.5146   × 10 8 2.6145   × 10 6 9
PSO5.2158   × 10 8 1.2485   × 10 6 8
GSA6.7624   × 10 7 5.2176   × 10 4 7
TLBO3.2648   × 10 7 7.5423   × 10 3 6
GOA2.7592   × 10 7 8.6427   × 10 2 5
GWO2.4257   × 10 7 6.5654   × 10 2 4
WOA2.1739   × 10 7 2.7865   × 10 2 2
MPA2.2365   × 10 7 1.4552   × 10 2 3
DPO2.1153   × 10 7 7.51421

Share and Cite

MDPI and ACS Style

Dehghani, M.; Mardaneh, M.; Guerrero, J.M.; Malik, O.P.; Ramirez-Mendoza, R.A.; Matas, J.; Vasquez, J.C.; Parra-Arroyo, L. A New “Doctor and Patient” Optimization Algorithm: An Application to Energy Commitment Problem. Appl. Sci. 2020, 10, 5791. https://doi.org/10.3390/app10175791

AMA Style

Dehghani M, Mardaneh M, Guerrero JM, Malik OP, Ramirez-Mendoza RA, Matas J, Vasquez JC, Parra-Arroyo L. A New “Doctor and Patient” Optimization Algorithm: An Application to Energy Commitment Problem. Applied Sciences. 2020; 10(17):5791. https://doi.org/10.3390/app10175791

Chicago/Turabian Style

Dehghani, Mohammad, Mohammad Mardaneh, Josep M. Guerrero, Om Parkash Malik, Ricardo A. Ramirez-Mendoza, José Matas, Juan C. Vasquez, and Lizeth Parra-Arroyo. 2020. "A New “Doctor and Patient” Optimization Algorithm: An Application to Energy Commitment Problem" Applied Sciences 10, no. 17: 5791. https://doi.org/10.3390/app10175791

APA Style

Dehghani, M., Mardaneh, M., Guerrero, J. M., Malik, O. P., Ramirez-Mendoza, R. A., Matas, J., Vasquez, J. C., & Parra-Arroyo, L. (2020). A New “Doctor and Patient” Optimization Algorithm: An Application to Energy Commitment Problem. Applied Sciences, 10(17), 5791. https://doi.org/10.3390/app10175791

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