Next Article in Journal
The Potential of the Biodigester as a Useful Tool in Coffee Farms
Previous Article in Journal
Impact of Tourist Behavior on the Discharge of Sunscreen Contamination in Aquatic Parks, Sinkholes, and Beaches of the Mexican Caribbean
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimization-Based Strategy for Solving Optimal Power Flow Problems in a Power System Integrated with Stochastic Solar and Wind Power Energy

1
Department of Electrical Engineering, Wah Engineering College, University of Wah, Wah Cantt 47040, Pakistan
2
Department of Electrical and Electronics, Global College of Engineering and Technology, Muscat 112, Oman
3
Department of Unmanned Vehicle Engineering, Sejong University, Seoul 05006, Korea
4
Department of Electrical Engineering, University of Lahore (Islamabad Campus), Islamabad 54590, Pakistan
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2021, 11(15), 6883; https://doi.org/10.3390/app11156883
Submission received: 2 July 2021 / Revised: 19 July 2021 / Accepted: 20 July 2021 / Published: 27 July 2021
(This article belongs to the Section Energy Science and Technology)

Abstract

:
In an effort to reduce greenhouse gas emissions, experts are looking to substitute fossil fuel energy with renewable energy for environmentally sustainable and emission free societies. This paper presents the hybridization of particle swarm optimization (PSO) with grey wolf optimization (GWO), namely a hybrid PSO-GWO algorithm for the solution of optimal power flow (OPF) problems integrated with stochastic solar photovoltaics (SPV) and wind turbines (WT) to enhance global search capabilities towards an optimal solution. A solution approach is used in which SPV and WT output powers are estimated using lognormal and Weibull probability distribution functions respectively, after simulation of 8000 Monte Carlo scenarios. The control variables include the forecast real power generation of SPV and WT, real power of thermal generators except slack-bus, and voltages of all voltage generation buses. The total generation cost of the system is considered the main objective function to be optimized, including the penalty and reserve cost for underestimation and overestimation of SPV and WT, respectively. The proposed solution approach for OPF problems is verified on the modified IEEE 30 bus test system. The performance and robustness of the proposed hybrid PSO-GWO algorithm in solving the OPF problem is assessed by comparing the results with five other metaheuristic optimization algorithms for the same test system, under the same control variables and system constraints. Simulation results confirm that the hybrid PSO-GWO algorithm performs well compared to other algorithms and shows that it can be an efficient choice for the solution of OPF problems.

1. Introduction

The rapid increase in energy demand and higher prices of conventional thermal generators led the world to a challenging problem of greenhouse gas emissions and an unstable economy. In this regard, energy production sectors are trying to integrate environmentally friendly and cost-effective renewable energy sources (RES) into the grid [1]. The integration of RES provides economic, social and environmental benefits for a sustainable society. Nowadays solar photovoltaics (SPV) and wind turbines (WT) are the two most integrated RES into the grid [2]. The incorporation of RES into the system is a tough and complicated task due to the sporadic nature of solar photovoltaics (SPV) and WT [3]. Due to the intermittent nature of RES, such as SPV and WT sources, it is hard to transmit power over a large distance between generation and consumers. This hurdle can be overcome by integrating enormous storage systems into the transmission network that will also help to reduces fluctuations in the network [4]. Energy storage devices, such as those that are lithium and sodium based [5,6], can be integrated along with conventional thermal and RES. Lithium-ion batteries have gained a lot of significance due to their high energy density, long life span and low self-discharge property, which can be added in the future expansion of microgrids [7]. OPF solutions provide the best operating values of control parameters, which are generators of real power and voltages, to reduce the total generation cost of the system [8]. Optimizing the generation cost system constraints, which are power flow balance, bus voltages, line capacity, and generator capability, is to be satisfied. OPF solutions must incorporate the uncertainties of these RES, which makes OPF a more complex and nonlinear problem.
For OPF problem formulation and solutions, many optimization algorithms and techniques have been developed and used in the past couple of years. These methods can be classified into classical methods, evolutionary based methods, and advanced metaheuristic algorithm-based methods. The classical method consists of linear programming [9], quadratic programming [10], non-linear programming [11], interior point method [12], and dynamic programming etc. For large scale problems, these classical methods are inadequate due to their complex and extensive calculation and their high dependency on initial guess values.
Evolutionary based algorithms copy the evolutionary idea of generating next generations from the parentages by mating to obtain the optimal solution. One of the most popular evolutionary based algorithms is the Genetic Algorithm, which has been used to solve OPF as mentioned in Dashtdar et al. [13], Moeini-Aghtaie et al. [14], and Bakirtzis et al. [15]. Another algorithm is the Differential Evolution (DE), that falls under this category and has been successfully implemented for constraint handling while optimizing the different objectives of OPF [16].
To overcome the discrepancies of the classical and evolutionary based methods, various metaheuristic algorithms have been effectively utilized for OPF solutions. Examples of various metaheuristic algorithms are particle swarm optimization (PSO), the moth-swarm algorithm (MSA), whale optimization algorithm (WOA), grey wolf optimization (GWO), and elephant herding optimization (EHO) etc. Physics based algorithms that are stimulated from different nature/physics laws and swarm-based algorithms which are inspired from the social manner of a cluster of animals have been included in the group of metaheuristic algorithms. The Gravitational Search Algorithm (GSA) has been suggested by Rashedi et al. in [17], that imitates the concept of Newton’s gravity law. Mirjalili et al., in [18], developed a Multi-Verse Optimization (MVO) algorithm which mimics the idea of cosmology that is comprised of wormholes, white holes and black holes. Physics based algorithms that have been utilized for the solution of OPF problems are GSA, presented by researchers in [19,20]. The Electromagnetism-like Mechanism Method (EMM) presented by the authors in [21] and Pinheiro et al. in [22] used the Doubly Modified logarithmic Barrier Function (M2BF) algorithm. Other swarm-based algorithms that have been utilized for OPF solutions exist such as PSO, presented by Kennedy and Eberhart in [23], Moth-flame Optimization (MFO) [24], Ant Lion Optimization (ALO) [25], GWO [26], and the Grasshopper Optimization Algorithm (GOA) [27]. The authors in [28,29,30] used PSO and its enhanced variants to provide an optimal solution for OPF problems. A novel meta-heuristic moth swarm algorithm (MSA) has been used for OPF solutions, which is inspired by the path finding processes of moths [31]. The performance of the moth swarm algorithm has been improved with arithmetic crossover (MSA-AC) and has been successfully utilized in the optimization of OPF problems for HVDC systems [32]. The results obtained show the better exploration and speedy convergence capability of MSA-AC for a global optimum. A modified sine cosine algorithm (MSCA) has been used for OPF solutions by the researcher in [33]. The effectiveness of modified MSCA has been verified in terms of a decrease in computational time and improvement of finding global OPF solutions with different objective functions. In reference [34], a novel improved artificial bee colony (IABC) optimization algorithm has been used for OPF solutions. The poor exploitation capability of the original ABC technique has been improved by an orthogonal learning approach, which guesses the optimum solution among the solution vectors created on restricted trials rather than a complete search. The dominance of IABC has been verified on IEEE 30 and 118 bus test systems, while optimizing the OPF solution, where IABC achieves a lower cost in comparison to other heuristics techniques for all test cases.
To improve the global searching ability, hybridization of two metaheuristic algorithms has been conducted to obtain the global optimum results; a hybrid of dragon fly and PSO (DA-PSO) have been proposed in reference [35] in which the exploitation and exploration capability of both algorithms are enhanced. The competence of the proposed algorithm (DA-PSO) has been confirmed on IEEE 30 and 57 bus test systems, while taking into account multi-objective and single-objective cases of power losses, cost and emissions. The dominance of the suggested algorithm has been proven by simulation results. However, computation time is increased due to the sequential operation of the hybrid technique. The authors in [36] present a hybrid DE along with a harmony search algorithm (Hybrid DE-HS) for OPF solutions. The usefulness of the applied hybrid DE-HS algorithm has been verified on IEEE 30, 118 and 300 bus test systems with reactive power compensation and transformer tap settings as a control variable. Other hybrid algorithms utilized for OPF solutions are hybrid PSO-GSA [37], imperialist competitive algorithm (ICA) hybridized with sequential quadratic programming by Ben Hmida et al. in [38], hybrid of WOA with Pattern Search Algorithm (PS) presented by authors in [39], hybrid of PSO and Artificial Physics Optimization (HPSOAPO) presented by Teeparthi et al. in [40], hybrid of PSO, and Shuffle Frog Leaping algorithm (HPSO-SFLA) presented by Narimani et al. in [41]. Other techniques that have been used with different multi-objective functions for OPF solutions are fuzzy logic-based techniques. A fuzzy logic-based approach has been used for the contingency constrained OPF problem formulation and solution in a decomposed form [42]. The authors in [43] used fuzzy based membership functions to form fuzzy optimal power problem formulations with an artificial bee colony (ABC) algorithm. In [44], the authors proposed an OPF solution using fuzzy logic with integrated GA and PSO. A hybrid of fuzzy based PSO techniques has also been used in [45], while considering the uncertainties in load demand and wind speed. The uncertainties in load demand and wind speed have been forecasted using fuzzy set notation.
Generally, there is not a single algorithm that can be considered better than others to solve all of the objectives of OPF problems, due to the inconsistent nature of OPF solutions. Thus, the performance of an algorithm or technique depends upon the nature and space of the problem that needs to be considered. Therefore, there must exist a space for new algorithms which can efficiently solve various objectives of OPF solutions. The objective of this paper is to accomplish a renewable integrated energy system with the use of WT and SPV sources to fulfill the energy requirements in the best possible way with a reduction in emissions for a sustainable environment. In this paper, utilization of the HPSO-GWO algorithms for OPF solutions incorporated with uncertain RES is presented. The results achieved with the suggested hybrid PSO-GWO algorithm (here after will be called HPSO-GWO) will be compared with PSO, GWO, success history based adaptive DE combined with superiority of feasible technique (SHADE-SF) [8], and with those results of GOA and Barnacles Mating Optimizer (BMO) reported in the literature [46]. The rest of the sections in this paper are set out as follows: Section 2 provides OPF problem mathematical formulations. Section 3 formulates frequency distributions of stochastic SPV and WT. Section 4 presents details of selected metaheuristic algorithms. Section 5 provides a detail discussion on the simulation results of each case. Section 6 provides a conclusion and future work.

2. Mathematical Formulation of OPF Problem

The standard IEEE 30 bus system consists of twenty-four load buses and six generator buses to which thermal generators are connected. In the proposed modified system, three thermal generators are replaced by two WT and a SPV source to formulate OPF, which is an optimization problem with inequality and equality constraints. The generation at a bus should not violate the line capacity constraints while replacing RES; in this study, the same method of replacing thermal generators with RES has been adopted as in in reference [8,46]. To forecast the output power of SPV and WT sources, lognormal and Weibull probability density functions (PDF) have been used, respectively. The application of the proposed HPSO-GWO algorithm is tested on the modified IEEE 30-bus system (hereafter test system) while considering different objective functions such as total generation cost, power losses, emissions, and cumulative voltage deviation in the system. The control parameters and other details of the test system under study are summarized in Table 1. Along with thermal generators, the adopted test system also includes SPV and WT, which are presented in Figure 1. The uncertain output power of wind and SPV sources must be compensated by the output power of thermal generators or reserve sources. The total generation cost of the entire system represents a combination of costs of all thermal generators, SPV and WT, plus reserve and penalty costs for overestimation and underestimation of wind and SPV sources, respectively. The subsequent sections explain cost function formulations for each type of source.

2.1. Cost Function Formulation for Thermal Generators

For the operation of thermal generators, they require fossil fuels and their output power depends on the consumption of fossil fuels. The fuel cost ($/h) and output power (MW) of the thermal generators have a relationship which can be expressed as below:
C T G = i = 1 n T a i + b i + c i P T G i 2
where C T G represents the total cost of n thermal generators producing P T G (MW) real power, while a i , b i and c i represent fuel cost coefficients for ith thermal generators whose values are mentioned in Table 2.
There is a greater divergence in the cost function of thermal generators with multi valve steam turbines. Therefore, for modeling a precise and more accurate cost function for thermal generators, valve point effects must be considered. Valve point effect is a sinusoidal function and must be included in the basic cost of thermal generators as an absolute value. The equation for the cost function, after counting the valve point effect, now becomes:
C T G = i = 1 n T a i + b i + c i P T G i 2 + | d i × sin { e i × ( P T G i min P T G i ) } |
where e i and d i represent the valve point effect coefficients of steam turbines. All cost and valve point coefficients used in Equations (1) and (2) are taken from reference [8] and are summarized in Table 2 along with emission coefficients.

2.2. Cost Function Formulation for Emission and Carbon Tax

During the operation of thermal generators, harmful emissions are released into the atmosphere due to consumption of fossil fuels. The release of these emissions increases with the escalation of output power from thermal generators. The expression for calculating the emission in ton/h is given below:
E e m i = [ i = 1 n T ( α i + β i P T G i + γ i P T G i 2 ) × 0.01 + ω i e ( μ T i P T G i ) ]
Emission coefficient values for thermal generators are given in Table 2 and are the same as in reference [8]. To control greenhouse gas emission, various countries are now putting immense pressure on power production sectors due the greenhouse effect [48]. As a result, the world has seen a shift towards solar photovoltaics, due to its tremendous potential across the globe and green renewable nature [49]. To reduce harmful emissions and encourage the integration of environmentally friendly and cost-effective RES such as SPV and wind, carbon tax is applied on a per unit quantity of produced emissions. The emission cost ($/h) for thermal generators is formulated as below:
C e m i = C t × E e m i

2.3. Cost Function Formulation for SPV and Wind Turbine

RES involves no running costs, unlike fossil fuel based thermal generation plants. RES per unit price is determined by either reimbursement cost to the initial expenditure or renovation and maintenance cost in case of possession by independent systems operators (ISO). RES may be operated and owned by private sectors; in such cases, ISO must pay a per unit price proportionate to the scheduled output energy of these plants. Direct cost of jth wind energy plants is mathematically modeled as a function of scheduled output power:
C w , j = g w , j × P w s c h , j
where P w s c h , j represents the schedule output power for jth wind turbine and g w , j represents the direct cost coefficient, linked with the same wind turbine.
Similarly, a direct cost assigned to kth SPV sources producing P s s c h , k scheduled power is:
C s , k = h k × P s s c h , k
where h k represents the direct cost coefficient of the kth SPV source.

2.4. Mathematical Modeling of Uncertainties in RES

There might be a situation in which the RES real power is below the forecasted power, which is referred to as over estimation. The system must be capable of having a reserve source to deliver the continuous power supply to consumers. The cost assigned to generating units of the reserve source, supplied to fulfill the deficit of overestimated units, is named the reserve cost.
Reserve cost for the overestimated units can be expressed as:
C R w , j = K R w , j ( P w s c h , j P w a c , j ) = K R w , j 0 P w s c h , j ( P w s c h , j P w a c , j ) f ( p w , j ) d p w , j
where P w s c h , j represents the scheduled power output from jth WT power and P w a c , j represents the actual power output from jth WT. K R w , j represents the reserve cost coefficient of jth WT. f ( p w , j ) represents PDF of WT, where details about probabilities calculation for different wind speeds can be obtained in Biswas et al. [8].
The reserve cost for the kth SPV source may be expressed as below:
C R s , k = K R s , k ( P s s c h , k P s a c , k ) = K R s , k f s ( P s a c , k < P s s c h , k ) [ P s s c h , k E ( P s a c , k < P s s c h , k ) ]
where P s s c h , k represents the scheduled output power from kth SPV source, P s a c , k represents actual power output from kth SPV source and K R s , k represents the reserve cost coefficient of the kth SPV source. f s ( P s a c , k < P s s c h , k ) represents the SPV source deficient amount of power than the scheduled SPV source power P s s c h , k . E ( P s a c , k < P s s c h , k ) represents the probability of SPV power below P s s c h , k .
On other hand, there might be a situation in which the estimated output power from SPV and WT is more than the estimated power. Under such conditions, the output power from RES is underestimated. The excess amount of energy is expected to be useless if not stored or used properly. ISO must pay a penalty cost for the additional amount of energy.
Penalty cost for the additional amount of energy for jth WT may be expressed in the form:
C P w , j = K P w , j ( P w a c , j P w s c h , j ) = K P w , j P w s c h , j P w r , j ( P w , j P w s c h , j ) f w ( p w , j ) d p w , j
Pwr,j represents the rated output power of jth WT, and K P w , j represents penalty cost coefficient of jth WT.
Similarly, ISO should pay a penalty cost for the additional amount of solar energy of kth SPV sources, and may be expressed below as:
C P s , k = K P s , k ( P s a c , k P s s c h , k ) = K P s , k f s ( P s a c , k > P s s c h , k ) [ E ( P s a c , k > P s s c h , k ) P s s c h , k ]
where K P s , k represents the penalty cost coefficient of the kth SPV source, and f s ( P s a c , k > P s s c h , k ) represents the possibility of the left-over amount SPV power. E ( P s a c , k > P s s c h , k ) P s s c h , k represents the probability of P s a c , k above P s s c h , k . The cost function evaluation of SPV and wind turbines are based upon the distribution of the lognormal PDF and Weibull PDF, respectively.

2.5. Objective Function Formulation

For the OPF solution we formulate two main objective functions, one without carbon tax and the other with the carbon tax. Both objective functions incorporate all of the cost functions of all generators.
First objective function:
F o b j 1 = C T G + j = 1 N w [ C w , j + C R w , j + C P w , j ] + k = 1 N s [ C s , k + C R s , k + C P s , k ]
where Ns and Nw represents the total number of SPV sources and WT in the network, respectively.
Second objective function:
F o b j 2 = F o b j 1 + C e m i
While optimizing these OPF objective functions, system inequality, equality, and line capacity constraints must be satisfied.

2.6. OPF Equality Constraints

OPF has two equality constraints, that are real (MW) and reactive power (MVAr) balances in the network. The total real and reactive power produced must compensate for the total power loss and demand in the network.
P G i P D i V i j = 1 N B V j [ G i j cos ( δ i j ) + B i j sin ( δ i j ) ] = 0 i N B
Q G i Q D i V i j = 1 N B V j [ G i j sin ( δ i j ) B i j cos ( δ i j ) ] = 0 i N B
Here, N B represents the total number of buses while B i j and G i j represent the susceptance and transfer conductance between bus i and j, respectively.

2.7. Inequality Constraints in OPF

For finding a solution to the OPF problem, there are several inequality constraints that should be satisfied. OPF inequality constraints are the generator voltages, real and reactive power limits, load bus voltages limits, and line capacity limits.
Generator inequality constraints are expressed as follow:
P T G i min P T G i P T G i max       i = 1 N TG
P W G j min P W G j P W G j max       j = 1 N WG
P S G k min P S G k P S G k max       k = 1 N SG
Q T G i min Q T G i Q T G i max       i = 1 N TG
Q W G j min Q W G j Q W G j max       j = 1 N WG
Q S G k min Q S G k Q S G k max       k = 1 N SG
V G i min V G i V G i max       i = 1 N G
Line and load bus inequality constraints:
V L m min V L m V L m max       m = 1 N l
S ln S ln max       n = 1 N l

2.8. Power Losses and Voltage Deviation

The most significant factor besides the reduction of total generation cost is the reduction of power loss in the network. Due to inherent line resistance, power losses are inevitable in transmission systems. The total network loss can be calculated as:
P l o s s = i = 1 n l j 1 n l G i j V i 2 + V j 2 2 V i V j cos ( δ i j )
The difference between the voltage angles of two buses is represented by δ i j = δ i δ j ; G i j represents transfer conductance and nl shows the total number of lines in the network.
The voltage profile of the system has been expressed in terms of voltage deviation (VD), which is a gauge of the voltage quality within the system. Voltage deviation is computed by the following expression:
V D = p = 1 n l | 1 V l p |

3. Mathematical Formulation for Stochastic SPV and WT Power

The adopted IEEE 30 bus system is changed at bus number 5 and 11 to WT, and at bus number 13 a SPV power source is integrated. Weibull PDF and wind speed distribution possess a similar shape [50]. A Weibull PDF with shape factor ( k ) and scale factor ( c ) can be used to obtain the probability of wind speed (m/s), as follows:
f v ( w ) = ( k c ) ( v c ) k 1 e ( v / c ) k       for   0 < v <
The output power of WT, which is a function of wind speed, can be modeled as follows:
P w ( v ) = { 0 , for   v < v i n   and   v < v o u t P w r , for   v i n v v r P w r , for   v r v v o u t
The values for the shape factor ( k ) and scale factor ( c ) are taken from reference [8] and are mentioned in Table 3. A Weibull fitting curve and wind speed frequency distribution curve are illustrated in Figure 2 and Figure 3, which are obtained after simulation of 8000 Monte Carlo scenarios.
Output power from SPV sources relies on solar irradiance that follows the lognormal PDF [51]. The lognormal PDF, having a mean ( μ ) and standard deviation ( σ ), may be used to obtain the probability for solar irradiance ( G ) by utilizing the following expression:
f G ( S ) = 1 G σ 2 π exp { ( ln x μ ) 2 2 σ 2 }   for   G > 0
The output power of the SPV source, which is a function of solar irradiance ( G ), can be mathematically modeled as follows:
P S P V ( G ) = { P s r ( G 2 G s t d R c ) for   0 < G < R c P s r ( G 2 G s t d R c ) for   G R c
The parameters for the lognormal PDF are taken from reference [8] and are presented in Table 4. After simulation of 8000 Monte Carlo scenarios, solar irradiance frequency distribution and a lognormal fitting curve is obtained and is illustrated in Figure 4, while Figure 5 shows the real power distribution of the SPV source at bus 13.

4. Optimization Algorithm for OPF Solutions

Several optimization algorithms and techniques are currently being developed and have been applied for OPF solutions. The common objective of all of these nature inspired algorithms is trying to achieve the best global optimum with a fast convergence speed. In this paper a HPSO-GWO algorithm is utilized to provide the best OPF solution for the adopted test system.

4.1. Particle Swarm Optimization

The PSO is a probabilistic, population-based algorithm that is stimulated from the social behavior of bird flocking. While searching for food, they either fly together or in a scattered way before reaching the food. During the search for food, there must be a bird that can sense the location of the food precisely, where it can be found. The one having the correct information of the food position communicates the message to the group. During this search process, they change their positions continuously until they have found the food.
This behavior is inspired from the social behavior of birds, to find the best optimum value of the function/problem, where every bird of the flock is named a particle. The position of every particle in the PSO algorithm can be updated with the following equation:
v i k + 1 = v i k + c 1 r 1 ( p i k x i k ) + c 2 r 2 ( g b x i k ) x i k + 1 = x i k + v i k + 1

4.2. Grey Wolf Optimization

Mirjalili et al. [26] in 2014 developed a GWO algorithm, which copies the leadership hierarchy and hunting mechanism of grey wolves. In the GWO algorithm, the hunting mechanism is implemented in three main steps to perform optimization, which are: searching, encircling, and attacking the prey. In the GWO algorithm, the leadership hierarchy is simulated by employing the grey wolves into four groups that are: alpha (α), beta (β), delta (ẟ) and omega (γ) wolves. The alpha (α) wolves, which are a male and female, lead the pack and are liable for making the decisions of sleeping, hunting and many other aspects. The second upper level of wolves in the pack are known as beta (β) wolves, which are the secondary wolves that can help alphas (α) in decision making and other pack actions. The beta (β) wolves reinforce the alphas’ (α) orders and provide feedback to alphas (α) during the pack actions. The third level of wolves in the pack are known as omegas (γ), which perform the character of a scapegoat. The omega (γ) wolves must always comply to all other leading alpha (α) and beta (β) wolves. The omegas (γ) are the last wolves that are permitted to eat, and do not have an important role in the hierarchy; however, if they lose the omegas (γ), the entire pack faces internal struggle and trouble. If a wolf does not belong to an alpha (α), beta (β) or omega (γ), she/he may be named as a delta (ẟ) or subordinate. Delta (ẟ) wolves dominate omegas, but they must comply to alpha (α) and beta (β) wolves in the command. To perform optimization, three main steps are to be executed, which are searching, encircling and hunting the prey.
The encircling movement of each wolf in the group is calculated through the following equation:
x ( t + 1 ) = x p ( t ) a · d d = | c · x p ( t ) x ( t ) |
where c and a are calculated as:
c = 2 · r 2 a = 2 l · r 1
Hunting behavior is mathematically expressed as:
x 1 = x α a 1 · d α d α = | c 1 · x α x | x 2 = x β a 2 · d β d β = | c 2 · x β x | x 3 = x δ a 3 · d δ d δ = | c 3 · x δ x | x 1 + x 2 + x 3 3 a ( · ) = 2 l · r 1 l c ( · ) = 2 · r 2 }
If A represents an arbitrary value between −2 a and 2 a , then an arbitrary value of |A| < 1 means that the wolves are enforced to encircle and attack the prey. While |A| > 1 means the wolves are forced to diverge from the prey. Exploration is the searching ability of GWO, whereas exploitation is the ability to attack the prey.

4.3. HPSO-GWO Optimizer

To find the global optimum solution with an excellent convergence capability is the common aim of all nature inspired metaheuristics algorithms. To find an optimal solution with better convergence performance, an algorithm should be furnished with the exploitation and exploration capabilities. Exploration is the capability of an algorithm to search the complete space of the solution, where exploitation is the capability to converge to the utmost result of the solution. The aim of all nature inspired hybrid algorithms is to establish a balance between these two capabilities to find a global optimal solution.
PSO has the best exploitation capability, but may become stuck at a local optimum, whereas the GWO has a best exploration capability of searching the whole search space, but has a slow convergence and poor local searching ability. A HPSO-GWO algorithm is developed in [52], in which the exploitation capability in PSO is improved with the exploration capability of the GWO to produce a combined algorithm strength.
In the developed HPSO-GWO algorithm, the position of delta, beta and alpha in the search space are updated with the following consecutive equations:
d δ = | c 3 · x δ ω * x | d β = | c 2 · x β ω * x | d α = | c 1 · x α ω * x |
In the above equations, both capabilities of grey wolves are balanced in the search area by an inertia constant. The flow chart of proposed HPSO-GWO algorithm is presented in Figure 6.
In order to combine both variants (GWO and PSO) the equations for velocity and position are modified as follows:
x i k + 1 = x i k + v i k + 1 v i k + 1 = ω * ( v i k + c 1 r 1 ( x 1 x i k ) + c 2 r 2 ( x 2 x i k ) + c 3 r 3 ( x 3 x i k ) )

5. Simulation Results and Discussion

The suggested HPSO-GWO algorithm is utilized for OPF solutions using MATLAB software. The simulation results are executed on a system with Intel (R) Core (TM) i5-4210U CPU @ 1.70 GHz 2.4 GHz with 4.0 GB RAM.
Simulations are carried out for several cases of adopted test systems, in which an SPV and two WT are integrated. Simulation results for each case obtained with utilization of the HPSO-GWO algorithm are presented with necessary description in this section. In the first two cases, variation in generation costs for SPV and WT are examined with variation in respective schedule powers and PDF parameters. In the third case, reduction of total generation cost is achieved without considering carbon tax and valve point effects for thermal generators. In the fourth case, total generation costs with valve point effects for thermal generators are considered. The fifth case optimizes the total generation cost with valve point effects of thermal generators, and a carbon tax of 20$/ton applied on emission. Each optimization algorithm is run five times, and a total of 500 iterations are executed in a whole run of the algorithm. During the complete run of each optimization algorithm for five times, the minimum cost of generation thus found and consequent control variable values are recorded.

5.1. Case 1: Total Cost, Reserve Cost, Direct Cost and Penalty Cost vs. PDF Parameters

In this case, total cost, penalty cost, direct cost and reserve cost of two WTs are plotted against Weibull PDF parameters, keeping the shape parameter ( k = 2) constant while the scale parameter ( c ) is changed from 2 to 16 and respective changes in costs are plotted at a fixed scheduled WT power. A suitable schedule WT power of 25 MW and 20 MW are selected for WT 1 and WT 2, respectively. The schedule power selected is about one-third of the total installed capability, which is reasonable because practical WTs have a capacity between 30% and 45% [53]. The curves for WT 1 and WT 2 are plotted in Figure 7 and Figure 8, respectively. The minimum value of total cost occurs around a value of 8 of the scale parameter ( c ); as the value of c increases, penalty cost increases while the reserve cost becomes constant after a certain value of c .
In a similar way, the SPV penalty cost, reserve cost and direct cost are plotted against lognormal PDF parameters from 2 to 7 at a fixed value of standard deviation ( σ = 0.6). A suitable value of scheduled SPV power is chosen at 20 MW. The minimum value of total cost occurs at a lognormal mean value of 5.5; the penalty and reserve cost have the same values at μ = 5.7. After this, the point penalty cost increases, due to which a sharp increase occurs in total cost. SPV output power is highly dependent on lognormal mean μ i.e., at a low value of μ , a high reserve is necessary and vice versa. The variation in SPV cost vs. lognormal mean ( μ ) is plotted in Figure 9.

5.2. Case 2: Total Cost, Penalty Cost, Reserve Cost and Direct Cost vs. Schedule Power

In this case, the scheduled wind power of two WT connected at bus 5 and 11 are varied from zero to WT rated power; variation of total cost, direct cost, reserve cost and penalty cost are plotted in Figure 8 and Figure 9. The coefficients of direct cost, reserve cost and penalty cost are taken from [8] and provided in Table 3. Simulation results show that reserve cost decreases as scheduled WT power increases, while penalty costs increase as scheduled power increases. There is a linear relationship between direct cost and wind schedule power. The total cost increased as schedule power increased, as it is the sum of direct cost, reserve cost and penalty cost.
Similarly, SPV total cost, direct cost, reserve cost and penalty costs are plotted against SPV scheduled power in Figure 10. The SPV different parameters and cost coefficients are provided in Table 4. The total minimum cost occurs somewhere between 15 MW and 20 MW in the case of solar power, while in the case of WT sources there is a linear relation between cost and wind schedule power, which can be observed from Figure 10, Figure 11 and Figure 12.

5.3. Case 3: Minimization of Total Generation Cost without Valve Point Effect

In this case, total generation costs without considering valve point effects of thermal generators is optimized. Cost coefficients and all other parameters for wind and SPV sources remain the same, as in Case 2. The optimized values obtained for all thermal generators, RES and other computed parameters are given in Table 5. The obtained results from PSO, GWO, HPSO-GWO and SHADE-SF are summarized in Table 5 for comparison. The convergence curve of four selected algorithms is shown in Figure 13. The voltage deviation VD, power losses, emissions and total generation costs are also provided in the table, and the best value obtained is highlighted in bold. The minimum generation cost obtained by HPSO-GWO is 772.736 $/h, and the worst result obtained by SHADE-SF is 773.979 $/h. The difference between the total generation cost obtained by HPSO-GWO and SHADE-SF is 1.243 $/h, which becomes a significant amount of (1.243 × 24 × 365) 10,889 $/Y. The minimum power losses and emissions are achieved by HPSO-GWO, while a minimum voltage deviation is achieved by SHADE-SF. Generator bus voltages, real powers and reactive powers for all algorithms are combined in Figure 14, Figure 15 and Figure 16, respectively. All of the values obtained by these algorithms are within specified limits, as shown clearly. The load bus voltages profile obtained through all algorithms is combined in Figure 17. The five run results of each algorithm are combined in Figure 18, and the minimum value achieved by each algorithm in the respective run of the simulation is indicated by bold.

5.4. Case 4: Minimization of Total Generation Cost with Valve Point Effect

In this case, total generation cost with valve point effect of thermal generators is considered as the objection function for optimization. Cost coefficients and all other parameters for wind and SPV sources remain the same as in Case 2. The obtained results from PSO, GWO, HPSO-GWO and SHADE-SF are summarized in Table 6 for comparison with the results obtained by BMO and GOA [46]. Comparing the results of Table 6 with Table 5, it is clear that penetration from thermal generators decreases as we add up valve point effects in the cost function. The convergence curve of four selected algorithms is shown in Figure 19. The voltage deviation VD, power losses, emissions and total generation costs are also provided in the table, and the best value obtained is highlighted in bold. The minimum generation cost obtained by HPSO-GWO is 780.587 $/h, and the worst result obtained by GOA is 785.711 $/h. The difference between total generation cost obtained by HPSO-GWO and GOA is 5.124 $/h, which become a significant amount of (5.124 × 24 × 365) 44,886.24 $/Y. The minimum power losses and emissions are achieved by PSO, while a minimum voltage deviation is achieved by SHADE-SF. Generator bus voltages, real powers and reactive powers for all algorithms are combined in Figure 20, Figure 21 and Figure 22, respectively. All of the values obtained by these algorithms are within the defined limits, as shown clearly. The load bus voltages obtained by all algorithms are combined in Figure 23. The five run results of each algorithm are combined in Figure 24, and the minimum value achieved by each algorithm in the respective run of each simulation is indicated by bold.

5.5. Case 5: Minimization of Total Generation Cost with Valve Point Effect and Carbon Tax

In this case, total generation cost with valve point effects and carbon tax imposed on emission of thermal generators is considered as the main objective function for optimization. The rate of carbon tax is kept the same, as in [8] and [46], that is 20$/ton. Cost coefficients and all other parameters for wind and SPV sources remained the same, as in Case 2. The obtained results from PSO, GWO, HPSO-GWO and SHADE-SF are summarized in Table 7 for a comparison with the results achieved by BMO and GOA [46]. Due to carbon tax imposed on emissions of thermal generators, their penetration decreases whereas RES penetration increases in the system; this can be noticed by comparing the results of Table 5, Table 6 and Table 7. The convergence curve of four selected algorithms is shown in Figure 25. The voltage deviation VD, power losses, emissions and total generation costs are also provided in the table, and the best value obtained is highlighted in bold. The minimum generation cost obtained by HPSO-GWO is 809.277 $/h, and the worst result obtained by GOA is 822.307 $/h. The difference between total generation cost obtained by HPSO-GWO and GOA is 13.03 $/h, which become a significant amount of (13.03 × 24 × 365) 114,142.8 $/Y. The minimum power losses and emissions are achieved by GWO while a minimum voltage deviation is achieved by SHADE-SF. Generator bus voltages, real powers and reactive powers for all algorithms are combined in Figure 26, Figure 27 and Figure 28, respectively. All of the values obtained by these algorithms are within specified limits, as shown clearly. The load bus voltages obtained by all algorithms are illustrated in Figure 29. The five run results of each algorithm are combined in Figure 30, and the minimum value achieved by each algorithm in the respective run of each simulation is indicated by bold.

6. Conclusions

This paper presents a solution to OPF problems, integrated with RES using a HPSO-GWO algorithm. The stochastic nature of SPV and WT is modeled with Lognormal and Weibull probability density functions, respectively. Important system constraints such as equality constraints, inequality constraints and line capacity constraints have been taken into consideration. The usefulness and robustness of the proposed HPSO-GWO algorithm is tested on adopted test systems integrated with RES. Based upon the simulation results obtained, it has been concluded that:
  • The proposed HPSO-GWO algorithm has been effectively applied to adopted test systems to find the optimal values of the control settings.
  • Minimization of total system generation costs without considering valve point effect, with valve effect, and valve effect with carbon tax imposed on emission) is an objective function which has been achieved using the HPSO-GWO algorithm.
  • The dominance of the HPSO-GWO algorithm in comparison with the PSO, GWO, SHADE-SF, GOA and BMO algorithms has been confirmed.
  • The HPSO-GWO has a speedy and steady convergence characteristic curve in comparison with PSO, GWO and SHADE-SF.
  • Penetration from renewable energy sources further increased, and emissions have been reduced by 0.864 ton/h after imposing a carbon tax on emissions.
In future, the HPSO-GWO algorithm will be effectively utilized in other power system expansions and planning studies, with the possible addition of hydro power generation in the system.

Author Contributions

Conceptualization, M.R., A.H. and M.U.A.; formal analysis, M.R. and M.I.M.; funding acquisition, S.J.H. and M.U.A.; investigation, M.R.; methodology, M.R. and A.Z.; project administration, A.H. and S.J.H.; Software, M.R.; supervision, A.H. and A.Z.; validation, S.J.H.; writing–original draft, M.R. and M.I.M.; writing–review & editing, A.H., S.J.H., M.U.A. and A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the Global college of Engineering and Technology, Muscat for funding this research under an internal research funding grant.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hannan, M.; Tan, S.Y.; Al-Shetwi, A.Q.; Jern, K.P.; Begum, R. Optimized controller for renewable energy sources integration into microgrid: Functions, constraints and suggestions. J. Clean. Prod. 2020, 256, 120419. [Google Scholar] [CrossRef]
  2. Zia, U.; Wang, S.; Radosavljević, J.; Lai, J. A solution to the optimal power flow problem considering WT and PV generation. IEEE Access 2019, 7, 46763–46772. [Google Scholar]
  3. Alsharif, M.H.; Kim, J.; Kim, J.H. Opportunities and Challenges of Solar and Wind Energy in South Korea: A Review. Sustainability 2018, 10, 1822. [Google Scholar] [CrossRef] [Green Version]
  4. Chandy, K.M.; Low, S.H.; Topcu, U.; Xu, H. A simple optimal power flow model with energy storage. In Proceedings of the 49th IEEE Conference on Decision and Control (CDC), Atlanta, GA, USA, 15–17 December 2010; pp. 1051–1057. [Google Scholar]
  5. Minakshi, M.; Sharma, N.; Ralph, D.; Appadoo, D.; Nallathamby, K. Synthesis and Characterization of Li(Co0.5Ni0.5)PO4 Cathode for Li-Ion Aqueous Battery Applications. Electrochem. Solid-State Lett. 2011, 14, A86–A89. [Google Scholar] [CrossRef]
  6. Sundaram, M.M.; Watcharatharapong, T.; Duraisamy, S.; Rao, P.T.; Munichandraiah, N.; Chakraborty, S.; Ahuja, R. Synthesis, and crystal and electronic structure of sodium metal phosphate for use as a hybrid capacitor in non-aqueous electrolyte. Dalton Trans. 2015, 44, 20108–20120. [Google Scholar] [CrossRef] [Green Version]
  7. Divakaran, A.M.; Minakshi, M.; Bahri, P.A.; Paul, S.; Kumari, P.; Divakaran, A.M.; Manjunatha, K.N. Rational design on materials for developing next generation lith-ium-ion secondary battery. Prog. Solid State Chem. 2020, 62, 100298. [Google Scholar] [CrossRef]
  8. Biswas, P.P.; Suganthan, P.; Amaratunga, G.A. Optimal power flow solutions incorporating stochastic wind and solar power. Energy Convers. Manag. 2017, 148, 1194–1207. [Google Scholar] [CrossRef]
  9. Mota-Palomino, R.; Quintana, V.H. Sparse Reactive Power Scheduling by a Penalty Function—Linear Programming Technique. IEEE Trans. Power Syst. 1986, 1, 31–39. [Google Scholar] [CrossRef]
  10. Burchett, R.C.; Happ, H.H.; Vierath, D.R. Quadratically convergent optimal power flow. IEEE Power Eng. Rev. 1984, 4, 34–35. [Google Scholar] [CrossRef]
  11. Habibollahzadeh, H.; Luo, G.-X.; Semlyen, A. Hydrothermal optimal power flow based on a combined linear and nonlinear programming methodology. IEEE Trans. Power Syst. 1989, 4, 530–537. [Google Scholar] [CrossRef]
  12. Jabr, R.A. A Primal-Dual Interior-Point Method to Solve the Optimal Power Flow Dispatching Problem. Optim. Eng. 2003, 4, 309–336. [Google Scholar] [CrossRef]
  13. Dashtdar, M.; Najafi, M.; Esmaeilbeig, M. Calculating the locational marginal price and solving optimal power flow problem based on congestion management using GA-GSF algorithm. Electr. Eng. 2020, 102, 1549–1566. [Google Scholar] [CrossRef]
  14. Moeini-Aghtaie, M.; Abbaspour, A.; Fotuhi-Firuzabad, M.; Hajipour, E. A decomposed solution to multiple-energy carriers optimal power flow. IEEE Trans. Power Syst. 2014, 29, 707–716. [Google Scholar] [CrossRef]
  15. Bakirtzis, A.; Biskas, P.; Zoumas, C.; Petridis, V. Optimal power flow by enhanced genetic algorithm. IEEE Trans. Power Syst. 2002, 17, 229–236. [Google Scholar] [CrossRef]
  16. Biswas, P.P.; Suganthan, P.; Mallipeddi, R.; Amaratunga, G.A. Optimal power flow solutions using differential evolution algorithm integrated with effective constraint handling techniques. Eng. Appl. Artif. Intell. 2018, 68, 81–100. [Google Scholar] [CrossRef]
  17. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A Gravitational Search Algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  18. Mirjalili, S.; Mirjalili, S.M.; Hatamlou, A. Multi-verse optimizer: A nature-inspired algorithm for global optimization. Neural Comput. Appl. 2016, 27, 495–513. [Google Scholar] [CrossRef]
  19. Duman, S.; Güvenç, U.; Sönmez, Y.; Yörükeren, N. Optimal power flow using gravitational search algorithm. Energy Convers. Manag. 2012, 59, 86–95. [Google Scholar] [CrossRef]
  20. Bhowmik, A.R.; Chakraborty, A.K. Solution of optimal power flow using nondominated sorting multi objective gravitational search algorithm. Int. J. Electron. Power Energy Syst. 2014, 62, 323–334. [Google Scholar] [CrossRef]
  21. Bouchekara, H.R.E.-H.; Abido, M.A.; Chaib, A.E. Optimal Power Flow Using an Improved Electromagnetism-like Mechanism Method. Electr. Power Compon. Syst. 2016, 44, 434–449. [Google Scholar] [CrossRef]
  22. Pinheiro, R.B.; Nepomuceno, L.; Balbo, A.R. Solving largescale reactive optimal power flow problems by a primaldual M2BF approach. Opt. Eng. 2020, 21, 485–515. [Google Scholar] [CrossRef]
  23. Eberhart, R.; Kennedy, J. A new optimizer using particle swarm theory. In Proceedings of the Sixth International Symposium on Micro Machine and Human Science (MHS’95), Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar]
  24. Mirjalili, S. Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm. Knowl. Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  25. Mirjalili, S. The Ant Lion Optimizer. Adv. Eng. Softw. 2015, 83, 80–98. [Google Scholar] [CrossRef]
  26. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  27. Saremi, S.; Mirjalili, S.; Lewis, A. Grasshopper Optimisation Algorithm: Theory and application. Adv. Eng. Softw. 2017, 105, 30–47. [Google Scholar] [CrossRef] [Green Version]
  28. Abido, M. Optimal power flow using particle swarm optimization. Int. J. Electr. Power Energy Syst. 2002, 24, 563–571. [Google Scholar] [CrossRef]
  29. Naderi, E.; Pourakbari-Kasmaei, M.; Abdi, H. An efficient particle swarm optimization algorithm to solve optimal power flow problem integrated with FACTS devices. Appl. Soft Comput. 2019, 80, 243–262. [Google Scholar] [CrossRef]
  30. Singh, R.P.; Mukherjee, V.; Ghoshal, S. Particle swarm optimization with an aging leader and challengers algorithm for the solution of optimal power flow problem. Appl. Soft Comput. 2016, 40, 161–177. [Google Scholar] [CrossRef]
  31. Mohamed, A.-A.A.; Mohamed, Y.S.; El-Gaafary, A.A.; Hemeida, A.M. Optimal power flow using moth swarm algorithm. Electr. Power Syst. Res. 2017, 142, 190–206. [Google Scholar] [CrossRef]
  32. Duman, S. A Modified Moth Swarm Algorithm Based on an Arithmetic Crossover for Constrained Optimization and Optimal Power Flow Problems. IEEE Access 2018, 6, 45394–45416. [Google Scholar] [CrossRef]
  33. Attia, A.-F.; El Sehiemy, R.A.; Hasanien, H.M. Optimal power flow solution in power systems using a novel Sine-Cosine algorithm. Int. J. Electr. Power Energy Syst. 2018, 99, 331–343. [Google Scholar] [CrossRef]
  34. Wenlei, B.; Eke, I.; Lee, K.Y. An improved artificial bee colony optimization algorithm based on orthogonal learning for optimal power flow problem. Control Eng. Pract. 2017, 61, 163–172. [Google Scholar]
  35. Khunkitti, S.; Siritaratiwat, A.; Premrudeepreechacharn, S.; Chatthaworn, R.; Watson, N.R. A hybrid DA-PSO optimization algorithm for multiobjective optimal power flow problems. Energies 2018, 11, 2270. [Google Scholar] [CrossRef] [Green Version]
  36. Reddy, S.S. Optimal power flow using hybrid differential evolution and harmony search algorithm. Int. J. Mach. Learn. Cybern. 2018, 10, 1077–1091. [Google Scholar] [CrossRef]
  37. Duman, S.; Li, J.; Wu, L.; Guvenc, U. Optimal power flow with stochastic wind power and FACTS devices: A modified hybrid PSOGSA with chaotic maps approach. Neural Comput. Appl. 2019, 32, 8463–8492. [Google Scholar] [CrossRef]
  38. Hmida, J.B.; Chambers, T.; Lee, J. Solving constrained optimal power flow with renewables using hybrid modified imperialist competitive algorithm and sequential quadratic programming. Electron. Power Syst. Res. 2019, 177, 105989. [Google Scholar] [CrossRef]
  39. Shilaja, C.; Arunprasanth, T. Optimal power flow using moth swarm algorithm with gravitational search algorithm considering wind power. Future Gener. Comput. Syst. 2019, 98, 708–715. [Google Scholar]
  40. Teeparthi, K.; Kumar, D.V. Multi-objective hybrid PSO-APO algorithm based security constrained optimal power flow with wind and thermal generators. Eng. Sci. Technol. Int. J. 2017, 20, 411–426. [Google Scholar] [CrossRef]
  41. Narimani, M.R.; Azizipanah-Abarghooee, R.; Zoghdar-MoghadamShahrekohne, B.; Gholami, K. A novel approach to multiobjective optimal power flow by a new hybrid optimization algorithm considering generator constraints and multi-fuel type. Energy 2013, 49, 119–136. [Google Scholar] [CrossRef]
  42. Ramesh, V.; Li, X. A fuzzy multiobjective approach to contingency constrained OPF. IEEE Trans. Power Syst. 1997, 12, 1348–1354. [Google Scholar] [CrossRef]
  43. He, X.H.; Wang, W.; Wang, Y.N.; Kong, J.; Geng, J.; Fan, S.B. Fuzzy optimal power flow with multi-objective based on artificial bee colony algorithm in power system. Appl. Mech. Mater. 2014, 448, 2473–2477. [Google Scholar]
  44. Kumar, S.; Chaturvedi, D. Optimal power flow solution using fuzzy evolutionary and swarm optimization. Int. J. Electr. Power Energy Syst. 2013, 47, 416–423. [Google Scholar] [CrossRef]
  45. Liang, R.-H.; Tsai, S.-R.; Chen, Y.-T.; Tseng, W.-T. Optimal power flow by a fuzzy based hybrid particle swarm optimization approach. Electr. Power Syst. Res. 2011, 81, 1466–1474. [Google Scholar] [CrossRef]
  46. Sulaiman, M.H.; Mustaffa, Z. Optimal power flow incorporating stochastic wind and solar generation by metaheuristic optimizers. Microsyst. Technol. 2020, 1–15. [Google Scholar] [CrossRef]
  47. Alsac, O.; Stott, B. Optimal Load Flow with Steady-State Security. IEEE Trans. Power Appar. Syst. 1974, PAS-93, 745–751. [Google Scholar] [CrossRef] [Green Version]
  48. Yao, F.; Dong, Z.Y.; Meng, K.; Xu, Z.; Iu, H.H.; Wong, K.P. Quantum-inspired particle swarm optimization for power system opera-tions considering wind power uncertainty and carbon tax in Australia. IEEE Trans. Ind. Inform. 2012, 8, 880–888. [Google Scholar] [CrossRef]
  49. Ahmed, W.; Sheikh, J.A.; Ahmad, S.; Farjana, S.H.; Mahmud, M.A.P. Impact of PV system orientation angle accuracy on greenhouse gases mitigation. Case Stud. Therm. Eng. 2021, 23, 100815. [Google Scholar] [CrossRef]
  50. Roy, R.; Jadhav, H. Optimal power flow solution of power system incorporating stochastic wind power using Gbest guided artificial bee colony algorithm. Int. J. Electr. Power Energy Syst. 2015, 64, 562–578. [Google Scholar] [CrossRef]
  51. Chang, T.P. Investigation on frequency distribution of global radiation using different probability density functions. Int. J. Appl. Sci. Eng. 2010, 8, 99–107. [Google Scholar]
  52. Singh, N.; Singh, S.B. Hybrid algorithm of particle swarm optimization and grey wolf optimizer for improving convergence performance. J. Appl. Math. 2017, 2017, 2030489. [Google Scholar] [CrossRef]
  53. Black, V. Cost and Performance Data for Power Generation Technologies; National Renewable Energy Laboratory: Golden, CO, USA, 2012.
Figure 1. Adapted IEEE 30 bus network with integrated WT and SPV sources.
Figure 1. Adapted IEEE 30 bus network with integrated WT and SPV sources.
Applsci 11 06883 g001
Figure 2. Wind speed distribution at bus #5.
Figure 2. Wind speed distribution at bus #5.
Applsci 11 06883 g002
Figure 3. Wind speed distribution at bus #11.
Figure 3. Wind speed distribution at bus #11.
Applsci 11 06883 g003
Figure 4. Distribution of solar irradiance at bus #13.
Figure 4. Distribution of solar irradiance at bus #13.
Applsci 11 06883 g004
Figure 5. Distribution of SPV real power at bus #13.
Figure 5. Distribution of SPV real power at bus #13.
Applsci 11 06883 g005
Figure 6. Flow chart of the proposed HPSO-GWO algorithm.
Figure 6. Flow chart of the proposed HPSO-GWO algorithm.
Applsci 11 06883 g006
Figure 7. Cost vs. Weibull scale parameter ( c ) at bus #5.
Figure 7. Cost vs. Weibull scale parameter ( c ) at bus #5.
Applsci 11 06883 g007
Figure 8. Cost vs. Weibull scale parameter ( c ) at bus #11.
Figure 8. Cost vs. Weibull scale parameter ( c ) at bus #11.
Applsci 11 06883 g008
Figure 9. Cost vs. lognormal mean ( μ ) at bus #13.
Figure 9. Cost vs. lognormal mean ( μ ) at bus #13.
Applsci 11 06883 g009
Figure 10. Cost vs. wind schedule power (MW) (bus #5).
Figure 10. Cost vs. wind schedule power (MW) (bus #5).
Applsci 11 06883 g010
Figure 11. Cost vs. wind schedule power (MW) (bus#11).
Figure 11. Cost vs. wind schedule power (MW) (bus#11).
Applsci 11 06883 g011
Figure 12. Cost vs. SPV schedule power (MW) (bus#13).
Figure 12. Cost vs. SPV schedule power (MW) (bus#13).
Applsci 11 06883 g012
Figure 13. OPF convergence curves comparison (Case 3).
Figure 13. OPF convergence curves comparison (Case 3).
Applsci 11 06883 g013
Figure 14. Generator bus voltages comparison (Case 3).
Figure 14. Generator bus voltages comparison (Case 3).
Applsci 11 06883 g014
Figure 15. Generators real power comparison (Case 3).
Figure 15. Generators real power comparison (Case 3).
Applsci 11 06883 g015
Figure 16. Generators reactive powers comparison (Case 3).
Figure 16. Generators reactive powers comparison (Case 3).
Applsci 11 06883 g016
Figure 17. Load buses voltage profile comparison (Case 3).
Figure 17. Load buses voltage profile comparison (Case 3).
Applsci 11 06883 g017
Figure 18. Cost comparison of 5 runs (Case 3).
Figure 18. Cost comparison of 5 runs (Case 3).
Applsci 11 06883 g018
Figure 19. OPF convergence curves comparison (Case 4).
Figure 19. OPF convergence curves comparison (Case 4).
Applsci 11 06883 g019
Figure 20. Generator bus voltages comparison (Case 4).
Figure 20. Generator bus voltages comparison (Case 4).
Applsci 11 06883 g020
Figure 21. Generators real power comparison (Case 4).
Figure 21. Generators real power comparison (Case 4).
Applsci 11 06883 g021
Figure 22. Generators reactive powers comparison (Case 4).
Figure 22. Generators reactive powers comparison (Case 4).
Applsci 11 06883 g022
Figure 23. Load buses voltage profile comparison (Case 4).
Figure 23. Load buses voltage profile comparison (Case 4).
Applsci 11 06883 g023
Figure 24. Cost comparison of 5 runs (Case 4).
Figure 24. Cost comparison of 5 runs (Case 4).
Applsci 11 06883 g024
Figure 25. OPF convergence curves comparison (Case 5).
Figure 25. OPF convergence curves comparison (Case 5).
Applsci 11 06883 g025
Figure 26. Generator bus voltages comparison (Case 5).
Figure 26. Generator bus voltages comparison (Case 5).
Applsci 11 06883 g026
Figure 27. Generators real power comparison (Case 5).
Figure 27. Generators real power comparison (Case 5).
Applsci 11 06883 g027
Figure 28. Generators reactive powers comparison (Case 5).
Figure 28. Generators reactive powers comparison (Case 5).
Applsci 11 06883 g028
Figure 29. Load buses voltage profile comparison (Case 5).
Figure 29. Load buses voltage profile comparison (Case 5).
Applsci 11 06883 g029
Figure 30. Cost comparison of 5 runs (Case 5).
Figure 30. Cost comparison of 5 runs (Case 5).
Applsci 11 06883 g030
Table 1. Adopted test system details.
Table 1. Adopted test system details.
ParametersQuantityDescription/Ranges
Branches41[47]
Buses30[47]
Thermal generators3At bus #1 (slack), #2 and #8
SPV plant1At bus #13
WTs2At bus #5 and #11
Load connected-283.4 MW and 126.3 MVAr
Shunt compensators2At bus# 5(0.19 MVAr) and at bus# 24(0.04 MVAr)
Load buses24(0.95–1.05) p.u.
Voltages of voltage generation buses6(0.95–1.10) p.u.
Control variables11Schedule active power for 5 sources (except slack).
Generation buses voltages.
Table 2. Details of thermal generators cost and emission coefficients.
Table 2. Details of thermal generators cost and emission coefficients.
Bus # P T G i min
(MW)
P T G i max
(MW)
abcdeγβαωμT
150200020.00375180.0376.49−5.5544.0916.6670.0002
2208001.750.0175160.0385.638−6.0472.5433.3330.0005
8103503.250.00834120.0453.38−3.555.32620.002
Table 3. WT PDF parameters and cost coefficients.
Table 3. WT PDF parameters and cost coefficients.
WT
No.
Bus
No.
Prw
(MW)
Scale Factor
(c)
Scale Factor
(k)
Weibull MeanDirect Cost CoefficientReserve Cost CoefficientPenalty Cost Coefficient
157592 v = 7.976 m/s g w , 5 = 1.6 K R w , 5 = 3 K P w , 5 = 1.5
21160102 v = 8.862 m/s g w , 11 = 1.75 K R w , 11 = 3 K P w , 11 = 1.5
Table 4. SPV source PDF parameters and cost coefficients.
Table 4. SPV source PDF parameters and cost coefficients.
Bus
No.
Prs
(MW)
Mean
(μ)
Standard Deviation (σ)Lognormal MeanDirect Cost CoefficientReserve Cost CoefficientPenalty Cost Coefficient
135060.6 G = 483 W/m2 h k , 13 = 1.6 K R s , 13 = 3 K P s , 13 = 1.5
Table 5. Simulation results for Case 3.
Table 5. Simulation results for Case 3.
ParametersMin. ValueMax. ValuePSOGWOHPSO-GWOSHADE-SF
Swing generatorPTG1(MW)50140136.056134.457133.046135.247
Control variablesPTG2(MW)208038.81238.31238.41838.691
PWT5(MW)07539.64239.33739.34439.694
PTG8(MW)103510.00010.00810.00010.000
PWT11(MW)06033.69633.27733.25733.673
PSPV13(MW)05031.03733.76835.17432.203
V1(p.u.)0.951.101.11.11.11.073
V2(p.u.)0.951.101.0901.0871.0901.058
V5(p.u.)0.951.101.0691.0681.0681.034
V8(p.u.)0.951.101.0941.0881.0981.045
V11(p.u)0.951.101.11.11.11.099
V13(p.u.)0.951.101.0941.0931.0941.049
Generator reactive power QTG1(MVAr)−20150−11.682−18.768−6.721−2.494
QTG2(MVAr)−206018.86327.8084.77312.444
QWT5(MVAr)−303525.51023.52635.00022.954
QTG8(MVAr)−154040.00040.00040.00040.000
QWT11(MVAr)−253019.14719.22117.86230.000
QSPV13(MVAr)−202521.70621.75922.23514.879
Objective functionsGen. cost ($/h)--773.213773.566772.736773.979
Ploss (MW)--5.8415.8215.7636.107
Emission (t/h)--1.8901.7121.5701.798
V.D (p.u.)--1.0330.9991.0470.452
Table 6. Simulation results for Case 4.
Table 6. Simulation results for Case 4.
ParametersMin. ValueMax. ValuePSOGWOHPSO-GWOSHADE-SFBMO [46]GOA [46]
Swing generatorPTG1(MW)50140134.908135.123135.053134.908134.908134.909
Control variablesPTG2(MW)208029.05727.09027.63027.96626.60228.066
PWT5(MW)07543.87542.76943.29043.40643.81742.766
PTG8(MW)103510.00010.00010.00010.00010.00012.241
PWT11(MW)06037.19236.82237.22436.72736.04736.468
PSPV13(MW)05033.85037.30135.84536.17937.81235.505
V1(p.u.)0.951.101.1001.1001.0991.07041.0811.024
V2(p.u.)0.951.101.0891.0891.0881.05650.9501.097
V5(p.u.)0.951.101.0701.0711.0721.03481.0451.002
V8(p.u.)0.951.101.0861.0821.0991.09451.0491.075
V11(p.u)0.951.101.1001.0971.0991.09961.1001.061
V13(p.u.)0.951.101.0911.0991.0951.05311.0681.064
Generator reactive powers QTG1(MVAr)−20150−10.613−10.334−13.645−5.10018.050−20.000
QTG2(MVAr)−206018.26214.14422.18513.073−20.00060.000
QWT5(MVAr)−303525.05027.89425.07122.55230.344−5.628
QTG8(MVAr)−154040.00040.00040.00040.00040.00040.000
QWT11(MVAr)−253019.49617.90920.02330.00027.86919.217
QSPV13(MVAr)−202520.48723.05419.25916.49819.99825.000
Objective functionsGen. cost ($/h)--781.321780.924780.587781.293781.652785.711
Ploss (MW)--5.4825.5165.5325.785--
Emission (t/h)--1.7621.7861.7781.762--
Carbon tax--------
V.D (p.u.)--1.0211.0660.9780.466--
Table 7. Simulation results for Case 5.
Table 7. Simulation results for Case 5.
ParametersMin. ValueMax. ValuePSOGWOHPSO-GWOSHADE-SFBMO [46]GOA [46]
Swing generatorPTG1(MW)50140123.663122.783123.956123.123123.127134.908
Control variablesPTG2(MW)208033.35033.42433.60531.99431.94725.215
PWT5(MW)07546.06145.09846.22645.47045.40341.755
PTG8(MW)103510.00010.01010.00010.00010.00014.648
PWT11(MW)06038.83837.73939.03338.31038.27034.946
PSPV13(MW)05036.61639.37235.78539.78339.86538.053
V1(p.u.)0.951.101.1001.1001.1001.0701.0781.092
V2(p.u.)0.951.100.9501.0891.0901.0561.0640.958
V5(p.u.)0.951.101.1001.0721.0711.0351.0431.071
V8(p.u.)0.951.101.1001.0971.0881.1001.0471.095
V11(p.u)0.951.101.1001.0981.1001.0981.1001.002
V13(p.u.)0.951.101.1001.0961.0951.0511.0601.052
Generator reactive powersQTG1(MVAr)−2015012.477−14.26612.264−3.215−1.84941.855
QTG2(MVAr)−2060−20.00020.714−20.00010.73812.406−20.000
QWT5(MVAr)−303535.00023.31035.00022.23222.91835.000
QTG8(MVAr)−154040.00040.00040.00040.00035.68640.000
QWT11(MVAr)−253019.78918.08919.91630.00028.506−0.900
QSPV13(MVAr)−202524.81923.34024.97915.95017.09421.178
Objective functionsGen. cost ($/h)--810.071810.154809.277810.404810.798822.307
Ploss (MW)--5.1285.0295.1325.281--
Emission (t/h)--0.8980.8540.9140.871--
Carbon tax--17.9617.0818.2817.42--
V.D (p.u.)--1.0331.0971.0270.462--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Riaz, M.; Hanif, A.; Hussain, S.J.; Memon, M.I.; Ali, M.U.; Zafar, A. An Optimization-Based Strategy for Solving Optimal Power Flow Problems in a Power System Integrated with Stochastic Solar and Wind Power Energy. Appl. Sci. 2021, 11, 6883. https://doi.org/10.3390/app11156883

AMA Style

Riaz M, Hanif A, Hussain SJ, Memon MI, Ali MU, Zafar A. An Optimization-Based Strategy for Solving Optimal Power Flow Problems in a Power System Integrated with Stochastic Solar and Wind Power Energy. Applied Sciences. 2021; 11(15):6883. https://doi.org/10.3390/app11156883

Chicago/Turabian Style

Riaz, Muhammad, Aamir Hanif, Shaik Javeed Hussain, Muhammad Irfan Memon, Muhammad Umair Ali, and Amad Zafar. 2021. "An Optimization-Based Strategy for Solving Optimal Power Flow Problems in a Power System Integrated with Stochastic Solar and Wind Power Energy" Applied Sciences 11, no. 15: 6883. https://doi.org/10.3390/app11156883

APA Style

Riaz, M., Hanif, A., Hussain, S. J., Memon, M. I., Ali, M. U., & Zafar, A. (2021). An Optimization-Based Strategy for Solving Optimal Power Flow Problems in a Power System Integrated with Stochastic Solar and Wind Power Energy. Applied Sciences, 11(15), 6883. https://doi.org/10.3390/app11156883

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