Computational Intelligence and Neuroscience

Computational Intelligence and Neuroscience / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 6192980 | 12 pages |

Optimal Operation of the Hybrid Electricity Generation System Using Multiverse Optimization Algorithm

Academic Editor: Antonino Laudani
Received12 Sep 2018
Revised31 Jan 2019
Accepted11 Feb 2019
Published11 Mar 2019


The ongoing load-shedding and energy crises due to mismanagement of energy produced by different sources in Pakistan and increasing dependency on those sources which produce energy using expensive fuels have contributed to rise in load shedding and price of energy per kilo watt hour. In this paper, we have presented the linear programming model of 95 energy production systems in Pakistan. An improved multiverse optimizer is implemented to generate a dataset of 100000 different solutions, which are suggesting to fulfill the overall demand of energy in the country ranging from 9587 MW to 27208 MW. We found that, if some of the power-generating systems are down due to some technical problems, still we can get our demand by following another solution from the dataset, which is partially utilizing the particular faulty power system. According to different case studies, taken in the present study, based on the reports about the electricity short falls been published in news from time to time, we have presented our solutions, respectively, for each case. It is interesting to note that it is easy to reduce the load shedding in the country, by following the solutions presented in our dataset. Graphical analysis is presented to further elaborate our findings. By comparing our results with state-of-the-art algorithms, it is interesting to note that an improved multiverse optimizer is better in getting solutions with lower power generation costs.

1. Introduction

Proper management and finding out optimal solutions to utilize the available power production sources is an important optimization problem in Pakistan, which is the main theme of this paper [111]. According to the reports in [12, 13], the demand of the power is increasing exponentially and which is causing an increase in load shedding each year [14]. The cause of this shortage is mismanagement of power production sources.

According to a report published by National Electric Power Regulatory Authority (NEPRA) in 2016 [15], the total capacity of 27240 MW is installed in the country, whereas the demand rises to 22000 MW in summer [16]. These production units can produce enough energy if all are managed optimally. In [17], these production units were normally producing up to 13240 MW with 4760 MW of shortfall. In nomenclature, we have indicated all the 95 installed capacities of electricity production with majority of them depending on thermal sources [15]. In [14], it is reported that the gap between demand and production is expected to increase in the coming years. It is further reported in the literature [14, 18] that urban areas are facing a load shading in summer season as of 10–12 hours, and in rural areas, it is even worse as 16–18 hours a day. Urbanization and rapid increase in population have caused an increase in the industrial zones [19], and it is worth noting that the demand of electricity was 16000 (MW) during summer of 2012 while the short fall was 8500 (MW) [20].

In this study, we have explored a scientific decision-making approach in terms of mathematical programming for providing a strategy to utilize the available production sources of Pakistan to meet the demand of energy in different seasons. However, in [21], a general mechanism was presented in terms of 5 major energy production types as to . They have fixed the per unit price to estimate the overall mathematical model. We have extended the model and included 95 production sources, which generates electricity using different means (oil, gas, wind, hydropower, and coal). We have implemented per unit price according to the actual price of fuel used to run these sources. The production capacities of each source are given in Table 1 [22]. Our model will help to reduce the shortfall by utilizing the power production sources properly.

Power unitLower limitUpper limit


In the last few years, evolutionary algorithms (EAs) are implemented to solve different real-world optimization problems. EAs simulate social behavior or natural procedures in order to handle optimization problems with best solution. Among the popular EAs are the Bat algorithm [23, 24], grasshopper optimization algorithm (GOA) [25], and firefly algorithm (FA) [26, 27]. The efficiency of metaheuristics is better as compared to classical optimization techniques in solving optimization problems with iterations and random search behavior. The main idea behind all EAs is the survival of the fittest, which in return increases the fitness of individuals in population. EAs are implemented with different frameworks, which includes population structure and search equations [28]. EAs with different frameworks generate random populations. Then, all solutions are evaluated according to a given objective function. The search space is explored and exploited in two phases. EAs get stuck in local optima due to their random search. Several papers are published in the literature to overcome this problem in EAs [29, 30].

The multiverse optimization algorithm (MVO) is a nature-inspired optimization technique given in [31]. The key idea of MVO is inspired from the theory of multiverses. Since its invention, it is applied to solve several real-world problems. In [31], it is shown that MVO obtained very competitive results compared with other metaheuristic optimization algorithms. Also, Faris et al. [32] employed the MVO for training the multilayer perceptions in neural network. The efficiency of their techniques was evaluated on different medical datasets from UCI repository. The outcome of their experiments is then compared to different metaheuristics, namely, Particle Swarm Optimization (PSO), Genetic Algorithm (GA), Differential Evolution (DE) Algorithm, and Cuckoo Search (CS). It was observed that MVO was improved in terms of convergence and avoidance of local optima. MVO still lacks behind in accuracy of results and convergence speed. Fewer studies, like Levy-flight-based MVO [33] and a quantum version of MVO [34], are recently proposed to further enhance the capabilities of MVO. Experimental results show that they have improved the solutions to some extent.

Currently, as can be seen in the former review, there are limited works on improving the performance of MVO. This work presents a new method of initialization with the MVO algorithm called improved multiverse optimization algorithm (IMVO). The experimental results on mathematical model of economic dispatch problem of electricity generation system in Pakistan and its three case studies show that IMVO is very competitive over FA, Bat algorithm, and GOA.

Furthermore, we have presented solutions to five case studies for proper management and production by all the 95 energy producers in the country. Our model can be implemented to reduce the load shedding and distribute the burden of production on respective power production units.

This work is organized as follows: Section 2 describes the concept of multiverse optimizer and improvements in it. The mathematical model for the problem is described in Section 3. Experimental settings, results, and discussions based on numerical simulations are given in Section 4. Conclusion of our findings is presented in Section 5.

2. An Improved Multiverse Optimizer

The MVO algorithm simulates the theory of multiverse, where three verses play main part in the whole algorithmic procedure: these are white holes, black holes, and wormholes. This is a population-based algorithm, where the population is searched in two phases: exploration in which the search space is visited very well and exploitation performs the search in neighborhood of a solution extensively. White hole and black hole, in the theory of multiverses, are used as exploration agents, and wormhole is acting as exploitation agent in this algorithm. In Figure 1, black, yellow, green, and pink are the objects which are moving through the white/black hole, and white points represent objects which move through the wormhole. There are 5 assumptions that are applied in MVO to the population/universe [31]:(a)The probability of having the white holes is directly proportional to the inflation rate.(b)Individuals (universes) having higher inflation rate send objects through white holes with higher probability.(c)Individuals (universes) having higher inflation rate send objects through white holes with higher probability.(d)Individuals (universes) having lower rate of inflation send objects through black holes with higher probability.(e)The objects present in all universes are moved randomly towards the current best individual (universe) through wormholes. This random process occurs without taking into account the inflation rate. A sketch of this algorithm is shown in Figure 1.

The initialization phase of MVO can be written as follows:where represents the matrix of candidate solutions, is the dimensions of the problem, and is the number of universes (population of solutions) which can be represented as follows:where is the parameter of universe, denotes the universe, NI () is the scaled inflation rate of the universe, , and is the parameter of universe selected by a roulette wheel selection mechanism. The objects are changed among the universes without perturbation. To maintain the balance between exploration and exploitation during the searching process in MVO, each universe is randomly treated to have wormholes to transport its objects through space. In order to ensure the exploitation around the current best solutions, particular wormhole tunnels are always established between a universe and the best universe formed so far, which can be represented mathematically as follows:where is the parameter of best universe formed so far, (travelling distance rate) and (wormhole existence probability) are coefficients, and are the lower bound and upper bound of variable, respectively, is the parameter of universe and r2, r3, and r4 are random numbers between 0 and 1. and are treated as adaptive formula as follows:where is the minimum (in this paper, it is set to 0.2), is the maximum (in this paper, it is set to 1), is the current iteration, is the maximum iterations, and is the exploitation accuracy over the iterations (in this paper, it is set to 6). Higher the value of , the sooner and more accurate exploitation/local search is performed [31]. It is worth to highlight that the MVO algorithm depends on number of iterations, number of universes, roulette wheel mechanism, and universe sorting mechanism. The quick sort algorithm is used to sort universe at each iteration, and roulette wheel selection is run for each variable in every universe in all iterations. Detailed description of MVO can be seen in [31].

3. Our Proposed Initialization Strategy

According to the trend in stochastic techniques, random initialization in each generation is the frequently used method to initialize algorithms without knowing any priory information about the problem in hand. The idea of random initialization is often pushing the algorithm to explore the given domain very well, while the exploitation around best solutions is not carried out by this type of initialization. We have used an initialization strategy in which IMVO is balanced in terms of exploration and exploitation. We have designed a two-stage initialization strategy: in the first stage, the algorithm is initialized randomly for a certain number of scattered solutions in the whole search domain. In stage two, after the algorithm completes a certain number of iterations and collects best points in the search domain, we initialize the IMVO with population of these best solutions to concentrate the algorithm on exploitation around the best solutions in the space. Hence, in current studies, we have assigned the function evaluations, fixed to solve a problem, in two parts. The first part initializes with random populations and collects the best solutions of each iteration by using the following equation:

In the second phase, the algorithm is directed to search in the neighborhood of the best solutions obtained in the first part and exploits the search space around the population of these best solutions.

4. Modeling of Electricity Production Cost as a Single Objective Linear Programming Problem

In this paper, we have presented an optimal mixture of energy production system as a linear programming (LP) model consisting of different decision variables, representing different electricity production sources of Pakistan, see Nomenclature [21]. Each variable is bounded and generates electricity through different fuel types like hydro, thermal, nuclear, wind, and solar energies. The suggested LP model includes different production sources as and represent different costs per for each source. The objective function represents the total cost () incurred to meet the demand of energy in different situations. Mathematically,where represents the cost per and are the different sources given in nomenclature. The objective of this study is to produce a data set of solutions so that we can select the required optimal solution fulfilling the given demand with lowest cost per as compared to other solutions in the data set. For more details, about the data set generated, the reader may refer to [35]. The problem is to search for best solutions among several candidate solutions and to meet the required power demand and satisfy constraints on each variable. These constraints are given in Equation (7) and Table 1 as follows.


Price in PKR390.831986252.289853201.3324084111.5048004
Variable (MW)15017.401515006.439815002.2666715067.46055


Price in PKR327.227212236.951036233.3069629155.6167319
Variable (MW)19076.394619017.23619023.0274919044.11301

Total demand constraint:

5. Experimental Settings, Results, and Discussion

We have implemented IMVO which is coded in Matlab R2016 [31]. The main features of our computer included 8 GB RAM with Intel® Core™ i7-7500U CPU @ 2.7 GHz 2.9 GHz with a 64-bit operating system Windows 10. For a fair comparison, the population size was taken as 100 and the number of iterations was 100. We have repeated our simulations 10 times. To elaborate the efficiency of IMVO, we have considered three special cases. The experimental outcome is compared with state-of-the-art algorithms: firefly algorithm, bat algorithm, and grasshopper optimization algorithm, as in Tables 2, 3, and 4. Also, we have picked 5 solutions from overall 100000 solutions and presented these solutions as different case studies with their respective costs per MWh (Table 5), which are based on the power demand in summers as provided in [20, 3638]. The power demands and decision variables, which represent the 95 production sources in Pakistan and cost incurred along with the total power produced, are given in Table 5. The individual variations of 20 best variables required to produce a particular solution in all case studies are furnished in Table 6, and a weighted comparison of the 95 energy sources is given in Figure 2. Looking at these graphs, the reader can deduce the importance of each production source of electricity in national grid of the country. One can deduce which solution to choose if a particular source is down due to some fault or which sources are very impotent and acting as a backbone in generating the required electricity. We have depicted five case studies graphically as shown in Figure 3. These graphs are collective representation of 95 production sources acting to produce a certain demand. It is interesting to note that the production of electricity from available sources is manageable, if proper techniques, like IMVO algorithm, are implemented to solve the problem. It is worth to note that we have pointed out the top 20 production sources of electricity in Pakistan according to the five case studies listed as in Table 5. These sources play vital part in the energy production as shown in Table 6. It is obvious that certain sources like are producing higher demand which are overburdened and could be relaxed by installing or expanding alternate sources. In Table 7, prices depending on different fuel types are listed for a megawatt of electricity. The convergence plots for these three case studies are depicted in Figures 46. It is interesting to note that, in all cases with demands of 15000, 19000, and 23000 megawatts, IMVO produced better solutions with lowest costs in millions (Table 24).


Price in PKR350.831986263.667565276.6948704187.660591
Variable (MW)23061.047623022.725423044.3165623064.1461

Case Study 1Case Study 2Case Study 3Case Study 4Case Study 5

Price in PKR1.12E + 081.51E + 081.84E + 081.97E + 082.16E + 08
Variables (megawatts)15003.161900323000.0825000.327208.36
TD ()2379.1371566.352938.16434783476.16
MD ()818.6771000905.7829763.4957983.48
GBHPP ()516.533714501023.3521347.3061421.054
WD ()130.7975210.4076177.7411106.4834181.5488
CB ()38.7471818475.42442174.4151163.8171
DKD ()42.56789116.3349116.020460.36495127.4962
AKHPP ()25.029090.415286119.003433.8710147.99439
KKHP ()6.0488615.9254441.720647219.27703
JagHP ()14.179293012.856058.99806621.38508
JabHP ()8.37593510.5708318.069.4384970.704657
RG ()6.8660513.3305270.3801664.9875559.745063
DHP ()1.42818713.5305612.65834202.044366
GHD ()6.052250.2807792.2192179.4011915.5921
NHP ()6.64702311.404537.5955076.01087313.86329
SHP ()11.5540411.20188.39161811.1670813.093
KHP ()1.8697161.1371321.5347352.6928041.023365
RKHP ()0.905093110.6516260.854247
CHPP ()0.8641390.8158020.8158020.9659260.830172
GTPP ()1305.3311655165516551645.54
MTPP ()640.0394546.27135013501350
JTPP ()149.5429850850850849.0554
FGTPP ()148.9675181.2005188.0609103.1949217.891
MGTPP ()4.38599974.8375155.4596145.8149168.1718
MGTPP ()127.05665.27735166.1534.1201110.9528
KGTPP ()112.8329116.05930.988320.567566149.5035
LTPP ()40.371490.88369224.3518120.5104117.4427
FSPP ()39.1961354.272585956.6365955.09061
SGTPP ()0.91407621.8678514.6641722.2310831.08445
PGTPP ()27.3005914.459970.5899288.48400923.07876
PTPP ()12.70353176.69568916.8505311.95603
KGTPS ()91.79519228.2708247127.6751246.8065
KPC ()91.128611000.07650316.213921.324039
GTPS ()21.3985930.868090.60831741.915131.135833
TPS ()1019.436852.5889921.79661132.2521260
CCPP 1 ()363.8285375.9421485.8476538.1132552.4153
CCPP 2 ()67.16314359.9164551.1955798.5884899.1415
AES LL ()12.48976362330.1038341.4665346.4409
AES PG ()87.79272289.7815365244.1401364.9033
AEL ()3.047575292912.0848519.5939
AP ()146.5024194.7207214.9626225225
AGL ()153.9627107.4664152.7011114.3378113.5217
CHMC PL ()75.13341265.198273.3111255.6046322.7874
DHA C L ()68.734550.4760996.54502862.9675770.30826
EPC ()23.4047798.65637.10004141.3346152.0636
EPQL ()130.7633214.0389222.3562226.5226.2002
FKPCL ()95.2459157151.264327.8556155.6898
FTSM ()56.16674104.4668100.824976.57812107.0448
FPCDL ()10.1152164.3581127.2432106.039177.7063
GHLLP ()101.70256.26072109.8266140.2346157.6599
GEPtvL ()162.3922188144.4615179.5249131.9891
GEL ()8.19373689.19682201.584.89962201.5
GAEL ()100.559292.3557257.18626110.482582.02294
HCPC ()111.335858.26757110.977243.47083111.5484
HPGC ()169.73130.807693192.3602125.0309215.355
HUBCO NPP ()147.8875225212.279209.5134223.5541
HUBCOHPP ()148.063412921167.51512921292
Iptvl ()21.53738158.6019141.6202164.2585164.9125
JPG ()84.1313712012065.22738119.424
KEL ()43.5068847.52342103.833230.58684130.8207
KAPCL ()811.44511223.2181185.7116381638
LPL ()153.00460.5240294.67205196.4468182.2104
LPTL ()124.189668.191750.95543138.8238198.6912
LEPCL ()255.6012589.7144578.5459660658.6634
NCP ()158.253820095.85135200151.3611
NPL ()32.48763200168.4331168.8323199.1059
OPCL ()205.28650.89022132.758665.89102220.0356
REPGCL ()30.20849100.207552.99079164163.4477
RP ()327.1605276.0325302.1523412411.8322
SPC ()72.0817575.0288391.21424109.1339112.6283
SPPQ ()101.0073225124.0716211.05221.962
SESL ()208.102326.559866.00281215.5041200.3566
SSEL ()220.540114.1361321.1493310.6587343.2993
SNPCL ()39.357344.145360.39879120.1177551.89066
SE ()77.2816611.0750666.3794423.8481616.50386
SEPC ()12.0147371.1123840.133847.4356183.28035
SPGL ()71.38235.8575610.8696988.37954113.884
TEL ()110.227125.38428113.567626.502894.01525
UchPL ()735.278128.16507990990989.7866
KANUPP ()50.415192.10906137105.7683134.7101
CHASNUPP-1 ()282.938686.3647275.8089325268.7866
CHASNUPP-2 ()200.2746243.2484149.4219325320.8773
CHASNUPP-3 ()191.8579340319.874301.6324214.0955
NPP ()18.16426173.693374.8168336.2965420.7756
QezSP ()49.71705106.903370.03342147.3463126.7687
YEL ()15.586822.33451146.2307424.4809112.86986
MWPCL ()41.0143421.3561627.9464236.4524245.21504
TGL ()9.46635725.200024947.7904647.90783
GAWPL ()15.3591918.9172245.1451615.2492627.94244
MWEL ()14.9682721.997755247.8787146.00231
FFCEL ()30.2686142.68928505015.58795
ZEP ()20.8876814.3734644.6704530.4588140.17728
TWEL ()15.2846812.16743029.5247928.65781
HCDPL ()9.87625520.0036247.2520531.903141.52962
FWEL 1 ()43.6920129.8009545.9802829.6855622.61687
FWEPL 2 ()6.404915506.87209321.6533440.28912

No. 1No. 2No. 3No. 4No. 5No. 6No. 7No. 8No. 9No. 10No. 11No. 12No. 13No. 14No. 15No. 16No. 17No. 18No. 19No. 20

Fuel typePrice in PKR per megawatt

Furnace oil11050
High-speed diesel17960
Regasified liquefied natural gas11270

6. Conclusions

In this paper, an improved multiverse optimization algorithm is implemented for solving a linear programming model for proper management of electricity production sources in Pakistan. To address the proper management of power production sources, several bound constraints, such as constraint for total demand, and bounds on all variables are considered. The objective function is linear with 95 decision variables and different coefficients. The IMVO algorithm is easy to implement with few parameters. We have found solutions for five different case studies and analyzed all solutions graphically. A link to the data set of 100000 different solutions is provided in this paper. In our results, we have obtained better solutions to meet the power demand in summer. It is recommended that current resources can produce the required power demand but some of the resources are required to be expanded in capacity along with installation of new sources. In future, we intend to include other factors, like power loss, in our model.


TD:Tarbela Dam ()
MD:Mangla Dam ()
GBHPP:Ghazi Barota Hydropower Project ()
WD:Warsak Dam ()
CB:Chashma Barrage ()
DKD:Duber Khwar Dam ()
AKPH:Allai Khwar Hydropower Plant ()
KKPL:Khan Khwar Hydropower Plant ()
JHP:Jagran Hydropower Plant ()
JabHP:Jabban Hydropower Plant ()
RB:Rasal Barrage ()
DHP:Dargai Power Plant ()
GZD:Gomal Zam Dam ()
NHP:Nandipur Hydropower Plant ()
SHP:Shadiwal Hydropower Plant ()
KHP:Khuram Hydropower Plant ()
RKHP:Renela Khuard Hydropower Plant ()
CHPP:Chatral Hydropower Plant ()
GTPP:Guduu Thermal Power Plant ()
MTPP:Muzaffargarh Thermal Power Plant ()
JTPP:Jamshoro Thermal Power Plant ()
FGTPP:Faisalabad Gas Turbine Power Plant ()
MGTPP:Multan Gas Turbine Power Plant ()
KGTPP:Kotri Gas Turbine Power Plant ()
LTPP:Larkana Thermal Power Plant ()
FSPP:Faisalabad Steam Power Plant ()
SGTPP:Shahdra Gas Turbine Power Plant ()
PGTPP:Panjgur Gas Turbine Power Plant ()
QTPP:Quetta Thermal Power Plant ()
PTPP:Pasni Thermal Power Plant ()
KPC:Korangi Power Complex ()
KGTPS:Korang Gas Turbine Power Plant Station ()
GTPS:Gas Turbine Power Station ()
KPS:Thermal Power Plant ()
CCPP1:Combined Cycle Power Plant 1 ()
CCPP2:Combined Cycle Power Plant 2 ()
AESLL:AES Lalpir Limited ()
AESPG:AES Pak Gen ()
AEL:Altern Energy Limited ()
AP:Altes Power ()
AGL:Attock Gen Limited ()
CMECP:CMEC Power Pvt. ()
DHACL:DHA Cogen Limited ()
EPC:Eastern Power Company ()
EPQL:Engro Powergen Qaidpur Limited ()
FKPP:Fauji Kabirwala Power Plant ()
FTSM:First-Star Modaraba ()
FPCDL:Foundation power company Daharki Limited ()
GHLLP:Grang Holding Limited Power Plant ()
GEPtvl:Green Electric Pvt. Limited ()
GEL:Gujranwala Energy Limited ()
DAEL:Gul Ahmad Energy Limited ()
HCPC:Habibullah Coastal Power Company ()
HPGC:Halmore Power Generation Company ()
HUBCONPP:HUBCO Narowal Power Plant ()
HUBCOHPP:HUBCO Hub Power Plant ()
IPtvL:Intergen PTV Limited ()
JPG:Japan Power Plant ()
KEL:Kohinor Energy Limited ()
KAPCL:Kot Addu Power Company Limited ()
LPL:Liberty Power Limited ()
LPTL:Liberty Power Tech. Limited ()
LEPCL:Lucky Electric Power Company Limited ()
NCP:Nishat Chunian Power ()
NPL:Nishat Power Limited ()
OPCPrvL:Orient Power Company (Pvt.) Limited ()
REPGC:Radian Energy Power Generation Company ()
RPK:Rousch Power, Khanewal ()
SPC:Saba Power Company ()
SPPQ:Saif Power Plant, Qadirabad ()
SECL:Sapphire Electric Company Limited ()
SSEL:Siddiqsons Energy Limited ()
SNPCL:Sindh Nooriabad Power Company (Pvt.) Limited ()
SE:Sitara Energy ()
SEPC:Southern Electric Power Company Limited ()
SPGL:Star Power Generation Limited ()
TEL:Tapal Energy Limited ()
UchPL:Uch (Uch-I and Uch-II) Power Limited ()
NPP:Nandipur Power Project ()
QeaSP:Quaid-e-Azam Solar Park ()
YEL:Yunus Energy Limited ()
MWPCL:Metro Wind Power Co Limited ()
TGL:Tenaga Generai Limited ()
GAWPL:Gul Ahmed Wind Power Limited ()
MWEL:Master Wind Energy Limited ()
FFCEL:FFC Energy Limited, Jhimpir ()
ZEPJ:Zorlu Enerji Pakistan, Jhimpir ()
TWEL:Tapal Wind Energy Limited ()
HCDPL:HydroChina Dawood Power Limited ()
FEW-1L:Foundation Wind Energy-I Limited ()
FEW-2PT:Foundation Wind Energy-II Private Limited ().

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

All authors contributed equally to the writing of this paper. And all authors read and approved the final manuscript.


This paper was supported by the Department of Mathematics and the Office of Research, Innovation, and Commercialization of Abdul Wali Khan University Mardan, Pakistan.


  1. A. Kumar and G. Shankar, “Quasi-oppositional harmony search algorithm based optimal dynamic load frequency control of a hybrid tidal-diesel power generation system,” IET Generation, Transmission & Distribution, vol. 12, no. 5, pp. 1099–1108, 2018. View at: Publisher Site | Google Scholar
  2. M. Mehdinejad, B. Mohammadi-Ivatloo, and R. Dadashzadeh-Bonab, “Energy production cost minimization in a combined heat and power generation systems using cuckoo optimization algorithm,” Energy Efficiency, vol. 10, no. 1, pp. 81–96, 2017. View at: Publisher Site | Google Scholar
  3. K. Y. Lee, C. Y. Chung, B. J. Huang et al., “A novel algorithm for single-axis maximum power generation sun trackers,” Energy Conversion and Management, vol. 149, pp. 543–552, 2017. View at: Publisher Site | Google Scholar
  4. L. S. M. Guedes, P. de Mendonca Maia, A. C. Lisboa, D. A. G. Vieira, and R. R. Saldanha, “A unit commitment algorithm and a compact MILP model for short-term hydro-power generation scheduling,” IEEE Transactions on Power Systems, vol. 32, no. 5, pp. 3381–3390, 2017. View at: Publisher Site | Google Scholar
  5. F. Zaman, S. M. Elsayed, T. Ray, and R. A. Sarker, “Evolutionary algorithms for power generation planning with uncertain renewable energy,” Energy, vol. 112, pp. 408–419, 2016. View at: Publisher Site | Google Scholar
  6. M. E. Nazari and M. M. Ardehali, “Optimal coordination of renewable wind and pumped storage with thermal power generation for maximizing economic profit with considerations for environmental emission based on newly developed heuristic optimization algorithm,” Journal of Renewable and Sustainable Energy, vol. 8, no. 6, Article ID 065905, 2016. View at: Publisher Site | Google Scholar
  7. S. Hemalatha and S. Chandramohan, “Cost optimization of power generation systems using bat algorithm for remote health facilities,” Journal of Medical Imaging and Health Informatics, vol. 6, no. 7, pp. 1646–1651, 2016. View at: Publisher Site | Google Scholar
  8. E. Assareh and M. Biglari, “A novel approach to capture the maximum power generation from wind turbines using hybrid MLP neural network and bees algorithm (HNNBA),” IETE Journal of Research, vol. 62, no. 3, pp. 368–378, 2016. View at: Publisher Site | Google Scholar
  9. N. Nikmehr and S. Najafi-Ravadanegh, “Optimal operation of distributed generations in micro-grids under uncertainties in load and renewable power generation using heuristic algorithm,” IET Renewable Power Generation, vol. 9, no. 8, pp. 982–990, 2015. View at: Publisher Site | Google Scholar
  10. W. Sun, J. Wang, and H. Chang, “Forecasting annual power generation using a harmony search algorithm-based joint parameters optimization combination model,” Energies, vol. 5, no. 10, pp. 3948–3971, 2012. View at: Publisher Site | Google Scholar
  11. J. Jasper, “Cost optimization of power generation using a differential evolution algorithm enhanced with neighbourhood search operation,” International Review of Electrical Engineering (IREE), vol. 7, no. 5, pp. 5854–5865, 2012. View at: Google Scholar
  12. G. Ali, “Factors affecting public investment in manufacturing sector of Pakistan,” European Journal of Economic Studies, vol. 13, no. 3, pp. 122–130, 2015. View at: Publisher Site | Google Scholar
  13. D. H. Husaini and H. H. Lean, “Does electricity drive the development of manufacturing sector in Malaysia?” Frontiers in Energy Research, vol. 3, p. 18, 2015. View at: Publisher Site | Google Scholar
  14. M. M. Rafique and S. Rehman, “National energy scenario of Pakistan-current status, future alternatives, and institutional infrastructure: an overview,” Renewable and Sustainable Energy Reviews, vol. 69, pp. 156–167, 2017. View at: Publisher Site | Google Scholar
  15. NEPRA, State of Industry Report, 2016,
  16. Electricity demand surges to 22,000 MW with rise in mercury, in Dawn, 2017,
  17. J. L. C. Meza, M. B. Yildirim, and A. S. M. Masud, “A multiobjective evolutionary programming algorithm and its applications to power generation expansion planning,” IEEE Transactions on Systems, Man, and Cybernetics-Part A: Systems and Humans, vol. 39, no. 5, pp. 1086–1096, 2009. View at: Publisher Site | Google Scholar
  18. M. Asif, “Sustainable energy options for Pakistan,” Renewable and Sustainable Energy Reviews, vol. 13, no. 4, pp. 903–909, 2009. View at: Publisher Site | Google Scholar
  19. S. Wasti, Economic Survey of Pakistan 2014-15, Government of Pakistan, Islamabad, Pakistan, 2015.
  20. M. Kugelman, Pakistan’s Energy Crisis, National Bureau of Asian Research, Washington, DC, USA, 2013.
  21. H. Zameer and Y. Wang, “Energy production system optimization: evidence from Pakistan,” Renewable and Sustainable Energy Reviews, vol. 82, pp. 886–893, 2018. View at: Publisher Site | Google Scholar
  22. B. Zafar, “Fuel Price Adjustment,” View at: Google Scholar
  23. X. S. Yang, “A new metaheuristic bat-inspired algorithm,” in Nature Inspired Cooperative Strategies for Optimization (NICSO 2010), pp. 65–74, Springer, Berlin, Germany, 2010. View at: Google Scholar
  24. B. Zhu, W. Zhu, Z. Liu, Q. Duan, and L. Cao, “A novel quantum-behaved bat algorithm with mean best position directed for numerical optimization,” Computational intelligence and neuroscience, vol. 2016, Article ID 6097484, 17 pages, 2016. View at: Publisher Site | Google Scholar
  25. S. Saremi, S. Mirjalili, and A. Lewis, “Grasshopper optimisation algorithm: theory and application,” Advances in Engineering Software, vol. 105, pp. 30–47, 2017. View at: Publisher Site | Google Scholar
  26. X. S. Yang, “Firefly algorithm, stochastic test functions and design optimisation,” 2010, View at: Google Scholar
  27. C. F. Chao and M. H. Horng, “The construction of support vector machine classifier using the firefly algorithm,” Computational Intelligence and Neuroscience, vol. 2015, Article ID 212719, 8 pages, 2015. View at: Publisher Site | Google Scholar
  28. F. Özkaynak, “A novel method to improve the performance of chaos based evolutionary algorithms,” Optik, vol. 126, no. 24, pp. 5434–5438, 2015. View at: Publisher Site | Google Scholar
  29. M. Sulaiman, A. Ahmad, A. Khan, and S. Muhammad, “Hybridized symbiotic organism search algorithm for the optimal operation of directional overcurrent relays,” Complexity, vol. 2018, Article ID 4605769, 11 pages, 2018. View at: Publisher Site | Google Scholar
  30. M. Sulaiman, S. Muhammad, and A. Khan, “Improved solutions for the optimal coordination of docrs using firefly algorithm,” Complexity, vol. 2018, Article ID 7039790, 15 pages, 2018. View at: Publisher Site | Google Scholar
  31. S. Mirjalili, Waseem, S. M. Mirjalili, and A. Hatamlou, “Multi-verse optimizer: a nature-inspired algorithm for global optimization,” Neural Computing and Applications, vol. 27, no. 2, pp. 495–513, 2016. View at: Publisher Site | Google Scholar
  32. H. Faris, I. Aljarah, and S. Mirjalili, “Training feedforward neural networks using multi-verse optimizer for binary classification problems,” Applied Intelligence, vol. 45, no. 2, pp. 322–332, 2016. View at: Publisher Site | Google Scholar
  33. C. Hu, Z. Li, T. Zhou, A. Zhu, and C. Xu, “A multi-verse optimizer with levy flights for numerical optimization and its application in test scheduling for network-on-chip,” PLoS One, vol. 11, no. 12, Article ID e0167341, 2016. View at: Publisher Site | Google Scholar
  34. G. I. Sayed, A. Darwish, and A. E. Hassanien, “Quantum multiverse optimization algorithm for optimization problems,” Neural Computing and Applications, vol. 28, pp. 1–18, 2017. View at: Google Scholar
  35. M. Sulaiman, “Optimal operation of the hybrid electricity generation system using multi-verse optimization algorithm,” Mendeley Data, vol. 1, 2019. View at: Google Scholar
  36. K. Kiani, Power cuts return as shortfall touches 7,000 MW, in Daily Dawn, 2017,
  37. K. Kiani, Record power generation of 18,900 MW achieved, in Daily Dawn, 2017,
  38. T. N. S. S. Reporter, Electricity shortfall widens to 6,000 MW, in Dawn, 2017,

Copyright © 2019 Muhammad Sulaiman et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

888 Views | 336 Downloads | 5 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.