Next Article in Journal
Ranking Importance of Topographical Surface and Subsurface Parameters on Paludification in Northern Boreal Forests Using Very High Resolution Remotely Sensed Datasets
Previous Article in Journal
Holistic Approach to Phosphorus Recovery from Urban Wastewater: Enhanced Biological Removal Combined with Precipitation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Salp–Swarm Optimizer and Accurate Forecasting Model for Dynamic Economic Dispatch in Sustainable Power Systems

1
College of Engineering at Wadi Addawaser, Prince Sattam Bin Abdulaziz University, KSA, Al-Kharj 16278, Saudi Arabia
2
Department of Electrical Engineering, Aswan University, Aswan 81542, Egypt
3
Department of Electrical Engineering and Automation, Aalto University, FI-00076 Espoo, Finland
4
Departament d’Enginyeria Informàtica i Matemàtiques, Universitat Rovira i Virgili, 43007 Tarragona, Spain
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(2), 576; https://doi.org/10.3390/su12020576
Submission received: 16 December 2019 / Revised: 5 January 2020 / Accepted: 9 January 2020 / Published: 12 January 2020

Abstract

:
Worldwide, the penetrations of photovoltaic (PV) and energy storage systems are increased in power systems. Due to the intermittent nature of PVs, these sustainable power systems require efficient managing and prediction techniques to ensure economic and secure operations. In this paper, a comprehensive dynamic economic dispatch (DED) framework is proposed that includes fuel-based generators, PV, and energy storage devices in sustainable power systems, considering various profiles of PV (clear and cloudy). The DED model aims at minimizing the total fuel cost of power generation stations while considering various constraints of generation stations, the power system, PV, and energy storage systems. An improved optimization algorithm is proposed to solve the DED optimization problem for a sustainable power system. In particular, a mutation mechanism is combined with a salp–swarm algorithm (SSA) to enhance the exploitation of the search space so that it provides a better population to get the optimal global solution. In addition, we propose a DED handling strategy that involves the use of PV power and load forecasting models based on deep learning techniques. The improved SSA algorithm is validated by ten benchmark problems and applied to the DED optimization problem for a hybrid power system that includes 40 thermal generators and PV and energy storage systems. The experimental results demonstrate the efficiency of the proposed framework with different penetrations of PV.

1. Introduction

Economic dispatch (ED) methods aim to schedule generating units and allocate the demand power among them to determine the best-generating scenarios. ED can benefit power utilities in various ways by systematically minimizing the cost of energy production consistent with the load demand. For this purpose, ED typically increases the usage of the most efficient generators, which can yield lower fuel costs and reduced carbon emissions. A complex mathematical model is solved through multiple computations to satisfy the demand while achieving the minimum generating costs of fuel-based generation stations. These computations are restricted by various constraints of power systems [1,2]. The typical constraints of the ED problem are the capacity of generators, the ramp-rate of generating units, and the power balance.
Conventionally, an approximate quadratic function is employed to mathematically formulate the ED problem in order to reduce the computational complexity. Nevertheless, in practice, the input–output curves of generators have nonlinear characteristic due to the multi-valve steam turbines (so-called valve-point effects). In addition, various faults in the machines prohibit some generators from operating in some zones (i.e., prohibited zones). The existence of the valve-point effects and prohibited zones make the solution space of the ED problem highly nonlinear, and thus increasing the complexity of the optimization process [3].
The static economic dispatch (SED) provides an economical solution to the power of generation stations at a certain load level. Differently, dynamic economic dispatch (DED) optimizes power generation for multiple load levels over a period of time (e.g., one day) [4]. Notably, DED is a more realistic procedure compared with SED because of considering the variation in the demand load and the variation of power generation over the studied period.
After formulating the DED problem mathematically, an optimization method is required to concurrently the fuel cost and carbon emission while handling all constraints. Although the traditional power systems have satisfied the demand loads for a long time, they will not be able to meet the current challenges alone, knowing that the fossil fuel is predicted to run out [5]. The integration of renewable energy sources (RES) into power systems can be attributed to environmental, economic, and social benefits. Driven by these benefits, the future power systems are predicted to be fed by energy sources that are totally renewable. Therefore, modern power systems combine non-conventional resources, such as RES and energy storage systems in order to provide sustainable power systems capable of meeting the significant increase in demand [6,7,8,9,10,11]. For a realistic DED framework and sustainable power systems, there is an essential need to incorporate RES and energy storage systems in the optimization model. The declining trend of the costs of battery energy storage with enhanced efficiency, along with an increasing need to alleviate the intermittent RES generation increased the penetrations of such energy storage systems in the transmission levels. A benefits of the storage systems is to curtail some of the variability challenges related to the RES power, and so smooth the fluctuated generation while providing further charging/discharging control options. Indeed, the DED optimization process can allow the optimal scheduling of the power outputs of generating units and RES according to the load demands over predefined time intervals. Concretely, RES contributes to reducing fuel costs and providing sustainable power and clean energy. Examples of RES are solar energy, wind energy, and tidal energy [12,13,14]. Indeed, solar energy is a necessary form of renewable energy. The earth receives a massive amount of energy from the sun that can be converted into clean electricity directly through photovoltaics (PVs) [15,16] or indirectly through concentrating solar power (CSP) [17]. Although the capital cost of PV farms is high, they are a preferable choice for generating clean electricity [18]. The integration of PV to power systems guarantees efficient power delivery and reduces the amount of CO2 emissions, thus protecting the surrounding environment.
In the last years, several heuristic optimization algorithms have been proposed for solving engineering problems in general due to their high performance and simplicity. These heuristic algorithms are mimicking natural phenomena or the social behaviors of creatures. For instance, particle swarm optimization (PSO) [19], genetic algorithm (GA) [20], ant–lion optimization (ALO) [21], ant colony optimization (ACO) [22], and grey wolf optimization (GWO) [23] are applied to solve the SED problem. The authors of [24] evaluated the performance of moth–flame optimization (MFO), moth swarm algorithm (MSA), GWO, ALO, sine cosine algorithm (SCA), and multi-verse optimization (MVO) with applying mutation operators in solving the SED problem.
Regarding DED, many optimization algorithms have been applied to solve the DED problem, such as symbiotic organisms search (SOS) algorithm, which combines GA, PSO, and SOS in a tri-base population [25]. In [26], the GA algorithm is implemented to optimize the demand side management and the DED as a complementary stage. An accelerated approach is proposed in [27] to solve the DED problem with high computational speed. A hybrid PSO algorithm called BBPSO is presented in [28] for solving the DED problem. In addition, chaotic differential bee colony optimization algorithm (CDBCO) [29], optimality condition decomposition (OCD) technique [30], biogeography-based optimization [31], differential evolution algorithm (DEA) [32,33], and hybrid flower pollination algorithm (HFPA) [34] are implemented to solve the DED problem. Furthermore, hybrid genetic algorithm and bacterial foraging (HGABF) approach [35], weighted probabilistic neural network and biogeography based optimization (WPNN–BBO) [36], multidisciplinary collaborative optimization (MCO) [37], quasi-oppositional group search optimization (QOGSO) [38], and improved real coded genetic algorithm (IRCGA) [39] are also applied to solve the DED problem. In [40], the MILP-IPM approach is applied to solve the DED, and it combines the mixed-integer linear programming (MILP) with the interior point method (IPM).
Indeed, the DED optimization problem is still a challenging issue because of the nonlinear and non-convex characteristics of DED objective functions. The DED optimization process also involves massive calculations to find the optimal values of tens of variables and parameters (high dimensionality) while satisfying the constraints of the power system. The nonlinear, non-convex, and non-differentiable characteristics of the DED, as well as the constraints, may squeeze the solution space. In addition, the intermittent nature of PV increases the difficulty of the DED optimization problem as it adds fluctuations to the power system. These aspects might push the search agents to deviate from the global optima and yield diverse local optima. To cope with these challenges, a robust optimizer is required to provide an accurate solution with low computational complexity.
To cope with the issues above, in this paper, we propose a comprehensive DED framework, deep learning-based forecasting models, and an improved optimizer. Specifically, these major three contributions can be summarized as follows:
  • Comprehensive DED framework: A comprehensive DED framework is formulated that includes fuel-based generators, PV, and storage devices in a sustainable power system, considering clear and cloudy profiles of PV.
  • Improved optimizer: We propose an improved salp–swarm optimizer that helps manage the global exploration of the DED algorithm and reach reasonable DED solutions. Specifically, we apply a mutation operator to the salp swarm optimizer to increase the exploitation of the search space for improved solutions. The proposed algorithm is validated with ten benchmark problems and then used to optimize the DED problem for a sustainable power system with PV within the studied period.
  • Deep learning-based forecasting models: We propose a DED handling strategy that involves the use of PV power and load forecasting models based on deep learning techniques.
The rest of this paper is organized as follows. Section 2 presents the DED framework. Section 3 explains the proposed DED solution algorithm and the deep learning-based forecasting models. Section 4 provides the results. Section 5 concludes the paper and gives some lines of future work.

2. Comprehensive DED Framework

The main objective of the DED problem is to minimize the costs of the fuel consumed by generator units. This optimization problem is heavily restricted by various practical considerations, e.g., generation constraints of units, ramp-rate constraints of generating units, and power mismatches constraints. Unlike the basic SED problems, DED aims to optimally dispatch the output power of all generators in a number of time instants. Furthermore, we have here considered the PV generation units and energy storage systems, complying with the global policy for establishing sustainable power systems. Figure 1 presents a general structure of sustainable power systems in which renewable energy sources and energy storage systems are interconnected, besides the fuel-based generation stations, to feed various loads.
To address the DED problem, most approaches published in the literature can be itemized into (1) meta-heuristic based optimization methods and (2) conventional optimization-based methods. Examples of meta-heuristic based optimization methods are [25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40]. Conventional optimization based methods involve: linear programming algorithm [41], lambda iteration method [42], and the interior algorithm [43]. These later methods are usually computationally efficient; however, they are suitable mainly with convex functions [44].
The optimal dispatch can be achieved through maximizing the total fuel-cost savings (TFS) of the fuel-based generator units, or in other words minimizing the total fuel-cost (TFC). During time intervals t { 1 , 2 , , T } , the corresponding demand power at these occasions must be optimally allocated among the generating units. The TFC is proportional to the number of thermal generators, and TFS can be formulated as follows:
m a x T F S = T F C 1 T F C 2 ,
in which
T F C = t = 1 T i = 1 n a i ( P i t ) 2 + b i P i t + c i + | e i | sin ( f i ( P i m i n P i t ) ) ,
where T F C 1 is the total fuel cost of all generator units to supply P d without PV and energy storage systems, and T F C 2 is the total fuel cost of all generator units to supply ( P d + P p v A v , t P p v c u r t , t ± P s ) with PV and energy storage systems. P d , P p v A v , t , P p v c u r t , t , and ± P s t are, respectively, the total predicted load power, available active power generation of PV based on the predicted environmental conditions (i.e., solar radiation), active power curtailment of PV, and the charging/discharging power of the energy storage systems at time instant t { 1 , 2 , , T } . For generator i t h , a i , b i , and c i are its cost coefficients, and e i , and f i are the coefficients of valve-point effects, n is the number of thermal generators to be scheduled, t represents the current time, T is the total time of the studied period, and P i t is the thermal power generated from the i t h generator at time t. The valve-point effect can be considered as a practical operation constraint of thermal generators. This effect introduces a ripple in the heat rate function, yielding a discontinuous nonlinear fuel cost function that has multiple local minima. In Equation (2), a second order quadratic cost function is formulated with a rectified sinusoidal term for precise modeling of the cost function of generators considering the valve-point effect.
The DED is mathematically formulated as:
DED : max T F S P i , P p v c u r t , t , P s , s . t . system constraints , s . t . PV constraints , s . t . Energy storage constraints .
The equality and inequality constraints of the power system, PV, and energy storage systems are listed below:
  • Equality constraints
    i = 1 n P i t + P p v t ± P s t = P d t + P l o s s t , t { 1 , 2 , , T } ,
    P p v t = P p v A v , t P p v c u r t , t , t { 1 , 2 , , T } ,
    P l o s s t = i = 1 n j = 1 n B i j P j t P i t , t { 1 , 2 , , T } .
  • Inequality constraints
    P i m i n P i t P i m a x , i , t { 1 , 2 , , T } ,
    P i t P i t 1 U R i , i , t { 1 , 2 , , T } ,
    P i t 1 P i t D R i , i , t { 1 , 2 , , T } ,
    P p v c u r t , t F p v c u r t , t P p v A v , t , t { 1 , 2 , , T } ,
    P s m i n P s t P s m a x , t { 1 , 2 , , T } ,
    S O C s m i n S O C s t S O C s m a x , t { 1 , 2 , , T } .
Constraint (4) represents the balance between the active powers of all thermal generators, the PV system, and the energy storage systems with respect to the total load demand and active power losses in the power system at each time instant t. P P V t is the output power from the PV unit at time t expressed by Label (5), P s t and + P s t are the charging and discharging power of the storage device at time t, respectively, P d t is the total demand power at time t, and P l o s s represents the transmission losses computed by Label (6) in which B i j represents B-coefficients.
Constraints (7), (8), and (9) represent the upper and lower operational boundaries, the ramp-up limit, and the ramp-down limit of each i t h thermal generator at each time instant t, respectively. P i m i n and P i m a x are the minimum and maximum limits of the output power from the i t h thermal generator, respectively. U R i and D R i represent the ramp-up and ramp-down boundaries of the i t h thermal generator, respectively.
Constraints (10) can set the maximum allowed curtailed power of PV according to the regulations of utilities. F p v c u r t , t is a factor where its value ranges from 0 to 1. In this work, this latter factor is set to zero to prevent the active power curtailment. Regarding the energy storage device, it has two operational constraints represented by (11) and (12). P s m i n and P s m a x are the minimum and maximum charging rate limits of the storage energy devices, respectively. S O C s m i n and S O C s m a x are the minimum and maximum state of charge limits of the storage energy devices, respectively. Note that the decision variables in the DED optimization problem are: the output power from each i t h thermal generator ( P i t ) , curtailed active PV power ( P P V c u r t , t ) , charging rate of the storage energy devices ( P s t ) , and state of charge of the storage energy devices ( S O C s t ) for all time instants during the dispatching period, which is the next day in this work.
Besides the consideration of valve point effects, the prohibited operating zones of generators can model their persistent physical operation boundaries (e.g., generator faults, excessive vibrations, turbine constraints). Therefore, the feasible zones for the ith generator can be mathematically formulated as follows:
P i m i n P i t P i , 1 L , P i , z 1 U P i t P i , z L , P i , z U P i t P i m a x ,
where z is the number of the prohibited operating zones for the i t h generator. P i , z L and P i , z U represent the lower and upper (MW) power boundaries of z t h prohibited operating zone for the i t h generator, respectively. Figure 2 illustrates the basic generator costs curve, the costs considering valve point effect, and the impacts of the prohibited operating zones on the two curves. Figure 2 implies that, if a generator has z prohibited zones, its operating region will be split into isolated feasible sub-regions, yielding multiple decision spaces for the DED problem. As noticed, by considering both valve-point effects and the prohibited operating zones of generators, the optimization model became more complex and non-convex, and so required robust optimization algorithms to be accurately solved.
The output power of a PV module depends on the ambient temperature and solar irradiation. It can be modeled as follows [45,46]:
T c t = T a t + N T 20 0.8 · R ,
I t = I s c = α i ( T c t T s ) ,
V t = V o c α v ( T c t ) ,
F F = V M P P I M P P V o c I S C ,
P p v A v , t = m c F F V t I t ,
where T c t is the cell temperature and T a t is the ambient temperature at time t, N T refers to the nominal operating cell temperature, R is the solar irradiation, I s c is the short circuit current, α i is the temperature coefficient of current, T s is the temperature of the standard test conditions, V o c is the open-circuit voltage, α v is the temperature coefficient of voltage, F F is the fill factor, m c is the number of cells in the module, and V M P P and I M P P are the voltage and current at the maximum power point, respectively. In this paper, one hour resolution of the data and the dispatch cycle have been considered. However, the proposed model is general, and can be applied to other resolutions.

3. Proposed DED Algorithm

The step-by-step procedure for the proposed method can be summarized as follows:
  • Step 1: Read system data, including cost coefficients of generation stations, power limits of generation stations, ramp rate limits for each generation station, and B-loss coefficients. Read the historical datasets of load and PV solar radiation, and data of energy storage systems.
  • Step 2: Set the population of the improved SSA algorithm, number of agents, and the maximum number of iterations.
  • Step 3: Forecast the PV radiation and load for the period in which the system is required to be optimally dispatched.
  • Step 4: Run the improved SSA algorithm considering the mutation operator and the handling strategy of the various constraints.
  • Step 5: Save and print the calculated results, including the scheduling of the generation stations, and energy storage systems, and the total costs.

3.1. Improved Salp–Swarm Algorithm (ISSA)

Swarm intelligence algorithms are widely used for solving complex optimization problems. Salp swarm algorithm (SSA) is a bio-inspired optimization algorithm where it mimics the swarming behavior of salps that belong to the family Salpidae [47]. Salps have a transparent body and look like jellyfishes, and they have a similar movement behavior in water. Salps almost form swarms called salp chains in deep oceans. The salp chain population consists of leaders and followers; the leader is the first salp in the swarm, where it tracks the food source and guides the remainder of the swarm followers. The position of all salps are saved in a two-dimensional matrix labeled S. The SSA algorithm assumes that there is a food source called F S in the search space as the swarm target. The positions of salps can be represented as:
S = s 1 , 1 s 1 , 2 s 1 , m s 2 , 1 s 2 , 2 s 2 , m s n , 1 s n , 2 s n , m ,
where n is the number of salps, and m is the number of variables. The leader’s position is updated according to the following formula:
S m 1 = F S m + x 1 ( ( u b m l b m ) x 2 + l b m ) x 3 0 , F S m x 1 ( ( u b m l b m ) x 2 + l b m ) x 3 < 0 ,
where, in the m t h dimension, S m 1 is the 1st salp, while F S m represents the food source, u b m and l b m are the upper and lower boundaries, respectively. X2 and X3 dictate if the next position in jth dimension should be towards positive infinity or negative infinity as well as the step size.
In SSA, x 1 is responsible for balancing exploration and exploitation. It can be calculated as follows:
x 1 = 2 exp ( 4 l L ) 2 .
In this equation, l represents the current iteration, and L represents the total number of iterations, x 2 , and x 3 are random numbers uniformly generated in the interval [0,1]. They control the step size and the next position in the m t h dimension.
The position of the followers is updated according to Newton’s law of motion as follows:
S m n = 1 2 a t 2 + v 0 t , n 2 ,
where n 2 , S m n is the n t h salp in the m t h dimension, t is time, v 0 is the initial velocity, and a = v final v 0 where v = x x 0 t . Since the time in the optimization problem represents iteration, the discrepancy between iterations equals 1 , and assuming v 0 = 0 , Equation (22) can be reformulated as follows:
S m n = 1 2 ( S m n + S m n 1 ) , n 2 ,
where S m n is the n t h salp in the m t h dimension.
The SSA optimization method can be trapped into local minima, which can lead to inaccurate results, especially with non-convex complex optimization problems, such as DED. To precisely solve the DED problem, we improve the performance of SSA by adopting a mutation operator. The improved SSA with the mutation operator is able to improve the convergence ability of the basic SSA algorithm by adjusting the searching direction in space in an adaptive way. Specifically, the mutation operator utilized in this paper is expressed as follows:
P m u t i = ( P r 1 i + P r 2 i + P r 3 i ) / 3 + ( η 2 η 1 ) ( P r 1 i P r 2 i ) + ( η 3 η 2 ) ( P r 2 i P r 3 i ) + ( η 1 η 3 ) ( P r 3 i P r 1 i ) ,
η 1 = | f ( P r 1 i ) | η * , η 2 = | f ( P r 2 i ) | η * , η 3 = | f ( P r 3 i ) | η * ,
where P r 1 i , P r 2 i , P r 3 i are the best three positions in the current population (i.e., the positions that have the top fitness values), η * = | f ( P r 1 i ) + f ( P r 2 i ) + f ( P r 3 i ) | , and f ( P ) is the fitness value.

3.2. Satisfying Various Constraints

When employing the improved SSA algorithm employed for solving the DED model, many updated solutions may be infeasible, especially in the early stage of the optimizer. To solve this issue, we utilize a reforming algorithm to enforce the updated solutions to be within their allowed limits, i.e., getting into the feasible search space. Consider S = [ P i P P V C u r t P s ] as an updated set of all candidate solutions, and it is formulated as:
P i = [ P 1 1 , , P n 1 n , , P 1 t , , P n t n , , P 1 T , , P n T n ] ,
P P V c u r t = [ P P V c u r t , 1 , , P P V C u r t , T T ] ,
P s = [ P s 1 , , P s T T ] ,
where P i , P P V c u r t , and P s represent, respectively, the set of power values of thermal generators, PV, and energy storage systems at all time periods t { 1 , 2 , , T } . First, the reforming algorithm enforces all variables with respected to their upper and lower boundaries for the variables at first time instant, i.e., t = 1. For instance, the upper and lower boundaries ( P i L o w , t , P i U P , t ) for the i t h generator at time period t can be formulated as follows:
P i L o w , t = P i min , if t = 1 , max P i min , P i t 1 D R i , otherwise ,
P i U p , t = P i max , if t = 1 , min P i max , P i t 1 D R i , otherwise .
Similar formulations are considered to enforce the control variables of PV and storage energy systems with respect to their upper and lower boundaries for the variables at first time instant ( t = 1 ). Then, for each time intervals t { 2 , , T } , all the feasible search space of the upper and the lower boundaries t + 1 at are explored. If the corresponding variables violate the explored feasible search space boundaries, their values are corrected to these upper/lower boundaries; otherwise, their values are kept the same.

3.3. Forecasting Models

Indeed, PV power forecasting is complex due to the fluctuating profiles of the weather (e.g., solar irradiance and temperature). Recurrent neural networks have been used in many applications, such as the optimal PV planning (i.e., determining the optimal sizes of PV plants with considering their intermittent generation), the optimal scheduling of thermal generators considering the predicted values of the PV generation for reducing operational costs of the grid, and the optimal management of PV plants while mitigating the operational problems, such as voltage violations and reverse power flow. Specifically, a deep long-short term memory (LSTM) recurrent neural network has been employed in the literature to forecast the intermittent PV power [48,49,50] and solar irradiance [51,52,53]. The merit of the LSTM is that it can model the temporal changes in the data due to their recurrent architecture and memory units. Unlike the traditional recurrent neural networks, LSTMs were designed to avoid the long-term dependency problem. LSTM can capture abstract concepts in the time-series forecasting problems.
Here, a DED handling strategy is introduced that involves the use of solar radiation (R) and load forecasting models. To build a solar radiation forecasting model, we propose the use of LSTM. Here, solar radiation forecasting is formulated as a regression problem. To train the solar radiation forecasting models, year-long solar radiation datasets with time resolution 1 h are used. The solar radiation dataset is restructured as follows: input and target. The input is the solar radiation at time steps { t i L B , t i 2 , t i 1 , t i } , and the target is the solar radiation at time step t i + 1 , where L B is the width of the lag window. The load dataset is restructured in the same way.
Figure 3 shows the block diagram of LSTM, where it receives an input sequence in which the activation units are used to trigger the gates. Each LSTM block has a cell with a state c t at time step t. The input gate i t , forget gate f t , and output gate o t are used to manage the reading or updating processes of this cell. The LSTM operation is expressed as follows:
i t = σ ( W x i X t + W h i h t 1 + W c i c t 1 + b i ) ,
f t = σ ( W x f X t + W h f h t 1 + W c f c t 1 + b f ) ,
c t = f t c t 1 + i t tanh ( W x c X t + W h c h t 1 + b c ) ,
o t = σ ( W x o X t + W h o h t 1 + W c o c t + b o ) ,
h t = o t tanh ( c t ) ,
where W are weights of the LSTM model that can be learned in the training stage. In this paper, we adopt our LSTM model proposed in [48] to construct a solar radiation forecasting model and load forecasting model, separately. The two LSTM models have a visible layer with one input, a hidden layer with four LSTM blocks, and an output layer that gives the target (solar radiation or load). The sigmoid activation function is used. The LSTM models are trained with the mean-squared error loss function, ADAM optimizer, and 100 epochs with a batch size of 1.
We use an adaptive moment estimation (ADAM) optimizer, a method for efficient stochastic optimization that only requires first-order gradients with little memory requirements. ADAM computes separate learning rates for each weight from estimates of first and second moments of the gradients as well as an exponentially decaying average of previous gradients. It is memory efficient since it does not keep a history of anything (just the rolling averages). We use an ADAM optimizer because it is computationally efficient (small memory requirements), has straightforward implementation, invariant to diagonally rescale the gradients, and suitable for optimizing the problems that are large in terms of data and for solving problems with dense noisy (or sparse gradients). The hyper-parameters of ADAM need little tuning. ADAM is also suitable for non-stationary objectives, which is the case of solar radiation forecasting.

4. Results and Discussion

To verify the effectiveness of the proposed method, we evaluate it using ten benchmark problems [47] (see Table 1) and the DED problem. In this study, we simulate a sustainable power system containing 40 thermal generators, PV, and energy storage systems. We provide the data of thermal generators in Table 2. We compare the results of the improved SSA (ISSA) with other optimization algorithms: SSA, MFO [54], and MVO [55]. The main inspiration for MFO is the navigation method of moths. In short, moths fly at night by keeping a fixed angle concerning the moon. The MFO mathematically models the navigation behaviour of moths to perform optimization. In turn, the main inspiration for the MVO algorithm relies on some concepts in cosmology: (1) white hole, (2) black hole, and (3) wormhole. The mathematical models of these concepts are implemented to perform optimization. It is a fact that there is no meta-heuristic best suited for treating all optimization models, complying with the No Free Lunch (NFL) theorem [56]. Therefore, we can conclude that each method of the three algorithms (SSA, MFO and MVO) gives a supervisor performance for specified types of functions, not all of them. In this work, we have developed an improved version of SSA to specifically solve the DED model. The parameters of the ISSA optimizer are set as follows: the population size is 30, the dimension is 10, and the maximum iterations number is 1000. All experiments have been carried out in MATLAB 2019a.

4.1. Analyzing the Performance of ISSA with Ten Benchmark Problems

The proposed method is implemented on ten benchmark problems, which are shown in Table 1. These problems have a challenging search space similar to the real complex search space. Indeed, they can verify the performance of algorithms and check their balancing strategy between the exploration and the exploitation phases. In Table 3, we show the average (Ave) and standard deviation (Std) values of the fitness function for 10 runs. The dimension of F1 is 30, the dimension of F2, F3, F4, F5, F6, F7, F8, and F9 is 30, and the dimension of F10 is 100. As shown, ISSA achieves results better than SSA with the benchmark problems (F1, F5, F7, and F9) in terms of Ave and Std, thanks to the mutation operator. Furthermore, ISSA outperforms MFO and MVO for the ten functions. In the case of F10, which is a highly nonlinear with a dimension of 100. With such complex optimization problem, ISSA achieves an Ave and Std values of 6.3844 and 1.3105, respectively. In turn, SSA and MVO achieves higher Ave values. MFO gives the worst Ave value (19.849) compared to the other optimizer. Driven by this superior performance, the proposed ISSA can be applied for the large scale optimization problems with several variables and wide search space, such as the DED problem. Note that the computational time of these methods is within 5 s.

4.2. Analyzing the Performance of an LSTM Forecasting Model

To evaluate the performance of the proposed method, two yearly PV datasets (dataset1 and dataset2) are used in this paper. We divide each dataset into training and testing datasets. A total of 70% of the samples are used to train the forecasting model, while the remaining samples are used for testing the model. We used the root mean square error (RMSE) to evaluate the performance of the forecasting models. RMSE can be defined as follows:
R M S E = 1 N i = 1 N ( D ^ i D i ) 2 .
In this equation, D i and D ^ i are the ith predicted and actual values, respectively, and N is the size of the testing dataset.
The performance of the LSTM forecasting method is compared with four PV forecasting methods: multiple linear regression (MLR), bagged regression trees (BRT), and ANN methods. These methods have been widely used in the literature for forecasting PV generation [58]. Thus, we used them in this paper to demonstrate the effectiveness of the LSTM model. As we can see in Table 4, LSTM gives small RMSE values with dataset1 and dataset2 (82.15 and 136.87). MLR and BRT give the highest RMSE values. With dataset1, the RMSE values of MLR and BRT are 384.90 and 494.46, respectively, while they are 329.11 and 416.212 with dataset2, respectively. Indeed, MLR and BRT were developed for stationary time series forecasting and thus they are not suitable for forecasting solar radiation. In the case of ANN, we have tried different configurations such as using one or two layers while changing the number of neurons from 1 to 50 with a step of 1. We found that the ANN model gives its best results with two layers and seven neurons. The RMSE values of ANN are 377.072 and 348.931 with dataset1 and dataset2, respectively. The forecasting RMSE values of LSTM are lower than MLR, BRT, and ANN methods. Indeed, MLR, BRT, and ANN methods do not have memory units, and thus they cannot handle the fluctuations in solar radiation. However, the architecture of the ANN method is a similar to LSTM, but it does not contain memory units or a recurrent architecture. In turn, LSTM uses the information learned in the previous time steps to estimate the current value, leading to accurate forecasting results.
Figure 4 shows the predicted solar irradiance for one day using the LSTM model for the dataset1 and dataset2. As shown, the LSTM can accurately predict the solar irradiance for the two datasets where the actual and predicted values are highly correlated. Notably, this high accuracy will have pronounced positive impacts on the DED problem since the decision variables will be properly computed considering accurate forecasting results. In this study, we demonstrate the effectiveness of deep learning approaches against traditional machine learning algorithms. Therefore, our recommendation is to utilize deep machine learning to accurately forecast PV power. In our paper, we utilize deep LSTM as an example, but other deep methods can give similar accuracy rates, but we do not investigate other deep learning, which is left as a future work.

4.3. Analyzing the Performance of ISSA with the DED Problem

Here, we simulate a sustainable power system containing 40 thermal generators, PV, and energy storage systems. The studied period is 24 h in which PV and load profiles are predicted using LSTM. The demand power is shown in Figure 5, and the output power of PV is shown in Figure 6 during clear and cloudy days. Below, we present two test cases of the proposed method with different weather conditions (clear and cloudy days). The capacity of the energy storage system (Sodium Sulfur Battery system) is 128 MWh. The parameters of the energy storage system are assumed as follows [59]: P s m i n = 0 , P s m a x = 8 M W , S O C s m i n = 50%, and S O C s m a x = 100%. The maximum allowed curtailed power ( F p v c u r t , t ) is set to zero.
Case 1: the performance of the proposed method with clear weather.Table 5 demonstrates the effectiveness of the proposed algorithm to minimize the total fuel costs with 0, 10, 20, 30, 40, and 50% penetrations of PV concerning the load. As we can see in this table, ISSA gives the lowest fuel cost for the six PV penetrations compared to SSA, MFO, and MVO. With ISSA, the fuel cost at 0% penetration of PV (i.e., without PV) is 2.7687*106 $ while it drops to 2.4228*106 $ at 50% penetration. It is worth noting that higher penetration levels of the PV can lead to reduced fuel costs. In the case of cloudy weather, it is expected that the output power of thermal generation units and energy storage devices will compensate the decrease of the PV power to satisfy the load, and so the fuel cost will be increased. Note that, if some clouds exist on a particular day, they will not affect the accuracy of the LSTM prediction model since the model is trained using a data of one year containing different weather conditions (clear, rainy, overcast, and cloudy).
Indeed, with higher penetrations of PV, the power generated from the thermal generators is less than the ones with lower penetrations. This power is optimally allocated among the thermal generators by ISSA through various generating scenarios over the studied period. Figure 7 shows the generated power of the 40 thermal generators over 24 h (see color bar) at the PV penetration of 50%. The x-axis of this figure represents the 40 generations (1, 2, …, 40) where the corresponding generated power values are shown. Each color on the color map represents the optimal generated power of the 40 thermal generators at each time interval (1, 2, …, 24). Note that each layer has a different color and a width. The color represents the time instance while the width represents the MW generated power of each generator. The width of each layer indicates the boundaries of the generated power of each generator, which can be easily tracked in this figure. For instance, the widths of layers corresponding to generators from 1 to 5 are narrow compared with generators from 15 to 20. This figure can quantify the daily generation of the different generators. For example, we can notice that generator number 20 has a very high production profile compared with generator 10. This can help the system operators to specify the generators that will be expected to be heavily loaded among the other generators according to the forecasting PV and load profiles as well as the solved DED problem.
Figure 8 shows the charging and discharging of the storage device for the different penetration cases. The ISSA algorithm can optimally schedule both the thermal generators and energy storage systems simultaneously to minimize the fuel cost. Figure 9 shows the convergence characteristics of ISSA. As we can see, the ISSA algorithm achieves better convergence characteristics, especially at 0%, 20%, and 30% penetration levels of the PV, where it converges rapidly compared to SSA.
Case 2: the performance of the proposed method with cloudy weather. In the case of cloudy weather, the output power of PV (Figure 6) is less than that of the clear weather as the output PV generation mainly depends on the solar irradiation. Therefore, the required powers from the thermal generators are increased causing higher fuel costs compared to case 1 (clear weather). Figure 10 shows the allocated load demand among the 40 thermal generators. Similar to Figure 7, the color map shows that the generators are within their operating boundaries (represented by layer width). It is clear that the scheduling of the 40 thermal generators with the cloudy weather profile is far different compared to the scheduling of the clear weather. Figure 11 shows the charging and the discharging of the storage device. Importantly, we notice higher storage rates with higher PV penetrations. Table 6 shows a comparison between the resulted fuel cost using ISSA, SSA, MFO, and MVO during cloudy days. Similar to the results of case 1 (clear day), ISSA achieves the lowest fuel cost. Figure 12 shows the convergence characteristics of ISSA, SSA, MFO, and MVO, in which ISSA has a stable and faster convergence. The improvement in performance when applying ISSA (with the mutation operator) is more significant in the case of 20%, 30%, and 50% PV penetrations than the other levels.
In this work, the optimization process is solved once, not repeated for each hour. This process aims to schedule the operations of generators day ahead considering daily forecast PV and load profiles. Therefore, the optimizer will be run offline (i.e., computational speed is not a vital factor), and where the computed optimal scheduling will be applied in the next day. However, the computational time of the optimizer when solving the DED is less than 10 min. These simulations reveal the effectiveness of the proposed algorithm to solve the DED problem with intermittent PV and energy storage systems. This analysis demonstrates that the proposed DED model can guarantee a minimum operational cost (total fuel costs) for sustainable power systems.

5. Conclusions and Future Work

Sustainable power systems provide the solution for satisfying the increased demand while reducing the harmful impacts of fuel-based generation stations. In this paper, a comprehensive DED framework has been proposed considering various profiles of PV (clear and cloudy) and energy storage systems. The DED model aims at minimizing the total fuel cost of power generation stations while considering various constraints of generation stations, PV, and energy storage systems. An improved optimization algorithm has been proposed to solve this optimization model. In particular, a mutation mechanism has been combined with SSA to produce to enhance the exploitation of the search space so that it provides a better population to get the optimal global solution. In addition, we have proposed a DED handling strategy that involves the use of PV power and load forecasting models based on deep learning models. The proposed algorithm has been validated with ten benchmark problems and then applied to the DED problem. A power system of forty thermal generators and PV is simulated considering the storage device with various PV penetrations (0%, 10%, 20%, 30%, 40%, and 50%). Two test cases have been analyzed: clear and cloudy weather conditions. ISSA obtains accurate results compared to three other optimization algorithms (SSA, MFO, and MVO). In future work, the proposed comprehensive DED model will be integrated with other sustainable energy resources (e.g., hydro energy, wind energy).

Author Contributions

Conceptualization, K.M. and M.A.-N.; methodology, K.M., E.M., and M.A.-N.; software, M.A.-N.; validation; formal analysis, K.M., E.M., and M.A.-N.; investigation, K.M., M.A.-N., and Z.M.A.; resources, K.M., M.A.-N., and Z.M.A.; data curation, E.M.; writing—original draft preparation, K.M., E.M., and M.A.-N.; writing—review and editing, K.M., M.A.-N., and Z.M.A.; visualization; supervision, K.M., M.A.-N. and Z.M.A.; project administration, K.M., M.A.-N., and Z.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, S.; Tan, S.; Xu, J.X. Consensus based approach for economic dispatch problem in a smart grid. IEEE Trans. Power Syst. 2013, 28, 4416–4426. [Google Scholar] [CrossRef]
  2. Abujarad, S.Y.; Mustafa, M.; Jamian, J. Recent approaches of unit commitment in the presence of intermittent renewable energy resources: A review. Renew. Sustain. Energy Rev. 2017, 70, 215–223. [Google Scholar] [CrossRef]
  3. Lee, F.N.; Breipohl, A.M. Reserve constrained economic dispatch with prohibited operating zones. IEEE Trans. Power Syst. 1993, 8, 246–254. [Google Scholar] [CrossRef]
  4. Zaman, M.; Elsayed, S.M.; Ray, T.; Sarker, R.A. Evolutionary algorithms for dynamic economic dispatch problems. IEEE Trans. Power Syst. 2015, 31, 1486–1495. [Google Scholar] [CrossRef]
  5. Rahman, M.M.; Mostafiz, S.B.; Paatero, J.V.; Lahdelma, R. Extension of energy crops on surplus agricultural lands: A potentially viable option in developing countries while fossil fuel reserves are diminishing. Renew. Sustain. Energy Rev. 2014, 29, 108–119. [Google Scholar] [CrossRef] [Green Version]
  6. Abas, N.; Kalair, A.; Khan, N. Review of fossil fuels and future energy technologies. Futures 2015, 69, 31–49. [Google Scholar] [CrossRef]
  7. Mohamed, A.A.A.; Ali, S.; Alkhalaf, S.; Senjyu, T.; Hemeida, A.M. Optimal Allocation of Hybrid Renewable Energy System by Multi-Objective Water Cycle Algorithm. Sustainability 2019, 11, 6550. [Google Scholar] [CrossRef] [Green Version]
  8. Alkhalaf, S.; Senjyu, T.; Saleh, A.A.; Hemeida, A.M.; Mohamed, A.A.A. A MODA and MODE Comparison for Optimal Allocation of Distributed Generations with Different Load Levels. Sustainability 2019, 11, 5323. [Google Scholar] [CrossRef] [Green Version]
  9. Ahmed, E.M.; Aly, M.; Elmelegi, A.; Alharbi, A.G.; Ali, Z.M. Multifunctional Distributed MPPT Controller for 3P4W Grid-Connected PV Systems in Distribution Network with Unbalanced Loads. Energies 2019, 12, 4799. [Google Scholar] [CrossRef] [Green Version]
  10. Abdel-Nasser, M.; Mahmoud, K.; Kashef, H. A Novel Smart Grid State Estimation Method Based on Neural Networks. IJIMAI 2018, 5, 92–100. [Google Scholar] [CrossRef] [Green Version]
  11. Mahmoud, K.; Hussein, M.M.; Abdel-Nasser, M.; Lehtonen, M. Optimal Voltage Control in Distribution Systems With Intermittent PV Using Multiobjective Grey-Wolf-Lévy Optimizer. IEEE Syst. J. 2019, 1–11. [Google Scholar] [CrossRef] [Green Version]
  12. Clancy, J.; Gaffney, F.; Deane, J.; Curtis, J.; Gallachóir, B.Ó. Fossil fuel and CO2 emissions savings on a high renewable electricity system–a single year case study for Ireland. Energy Policy 2015, 83, 151–164. [Google Scholar] [CrossRef]
  13. Mohamed, F.; AbdelNasser, M.; Mahmoud, K.; Kamel, S. Accurate economic dispatch solution using hybrid whale-wolf optimization method. In Proceedings of the 2017 Nineteenth International Middle East Power Systems Conference (MEPCON), Cairo, Egypt, 19–21 December 2017; pp. 922–927. [Google Scholar]
  14. Mahmoud, K.; Abdel-Nasser, M. Efficient SPF approach based on regression and correction models for active distribution systems. IET Renew. Power Gener. 2017, 11, 1778–1784. [Google Scholar] [CrossRef]
  15. Novakovic, B.; Nasiri, A. 1-Introduction to electrical energy systems. In Electric Renewable Energy Systems; Rashid, M.H., Ed.; Academic Press: Boston, MA, USA, 2016; pp. 1–20. [Google Scholar] [CrossRef]
  16. Mahmoud, K.; Abdel-Nasser, M. Fast-yet-Accurate Energy Loss Assessment Approach for Analyzing/Sizing PV in Distribution Systems using Machine Learning. IEEE Trans. Sustain. Energy 2018, 1025–1033. [Google Scholar] [CrossRef]
  17. Cirocco, L.R.; Belusko, M.; Bruno, F.; Boland, J.; Pudney, P. Maximising revenue via optimal control of a concentrating solar thermal power plant with limited storage capacity. IET Renew. Power Gener. 2016, 10, 729–734. [Google Scholar] [CrossRef]
  18. Shaahid, S.M. Economic feasibility of decentralized hybrid photovoltaic-diesel technology in saudi arabia. Therm. Sci. 2017, 21, 745–756. [Google Scholar] [CrossRef] [Green Version]
  19. Mahor, A.; Prasad, V.; Rangnekar, S. Economic dispatch using particle swarm optimization: A review. Renew. Sustain. Energy Rev. 2009, 13, 2134–2141. [Google Scholar] [CrossRef]
  20. Rastgoufard, S.; Iqbal, S.; Hoque, M.T.; Charalampidis, D. Genetic algorithm variant based effective solutions for economic dispatch problems. In Proceedings of the 2018 IEEE Texas Power and Energy Conference (TPEC), College Station, TX, USA, 8–9 February 2018; pp. 1–6. [Google Scholar]
  21. Kamboj, V.K.; Bath, S.; Dhillon, J. Solution of non-convex economic load dispatch problem using Grey Wolf Optimizer. Neural Comput. Appl. 2016, 27, 1301–1316. [Google Scholar] [CrossRef]
  22. Pothiya, S.; Ngamroo, I.; Kongprawechnon, W. Ant colony optimisation for economic dispatch problem with non-smooth cost functions. Int. J. Electr. Power Energy Syst. 2010, 32, 478–487. [Google Scholar] [CrossRef]
  23. Jayabarathi, T.; Raghunathan, T.; Adarsh, B.; Suganthan, P.N. Economic dispatch using hybrid grey wolf optimizer. Energy 2016, 111, 630–641. [Google Scholar] [CrossRef]
  24. Mostafa, E.; Abdel-Nasser, M.; Mahmoud, K. Performance evaluation of metaheuristic optimization methods with mutation operators for combined economic and emission dispatch. In Proceedings of the 2017 Nineteenth International Middle East Power Systems Conference (MEPCON), Cairo, Egypt, 19–21 December 2017; pp. 1004–1009. [Google Scholar]
  25. Duman, S. Symbiotic organisms search algorithm for optimal power flow problem based on valve-point effect and prohibited zones. Neural Comput. Appl. 2017, 28, 3571–3585. [Google Scholar]
  26. Mellouk, L.; Boulmalf, M.; Aaroud, A.; Zine-Dine, K.; Benhaddou, D. Genetic Algorithm to Solve Demand Side Management and Economic Dispatch Problem. Procedia Comput. Sci. 2018, 130, 611–618. [Google Scholar] [CrossRef]
  27. Chen, G.; Li, C.; Dong, Z. Parallel and distributed computation for dynamical economic dispatch. IEEE Trans. Smart Grid 2016, 8, 1026–1027. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Gong, D.W.; Geng, N.; Sun, X.Y. Hybrid bare-bones PSO for dynamic economic dispatch with valve-point effects. Appl. Soft Comput. 2014, 18, 248–260. [Google Scholar] [CrossRef]
  29. Lu, P.; Zhou, J.; Zhang, H.; Zhang, R.; Wang, C. Chaotic differential bee colony optimization algorithm for dynamic economic dispatch problem with valve-point effects. Int. J. Electr. Power Energy Syst. 2014, 62, 130–143. [Google Scholar] [CrossRef]
  30. Rabiee, A.; Mohammadi-Ivatloo, B.; Moradi-Dalvand, M. Fast dynamic economic power dispatch problems solution via optimality condition decomposition. IEEE Trans. Power Syst. 2013, 29, 982–983. [Google Scholar] [CrossRef]
  31. Ma, H.; Yang, Z.; You, P.; Fei, M. Multi-objective biogeography-based optimization for dynamic economic emission load dispatch considering plug-in electric vehicles charging. Energy 2017, 135, 101–111. [Google Scholar] [CrossRef]
  32. Jebaraj, L.; Venkatesan, C.; Soubache, I.; Rajan, C.C.A. Application of differential evolution algorithm in static and dynamic economic or emission dispatch problem: A review. Renew. Sustain. Energy Rev. 2017, 77, 1206–1220. [Google Scholar] [CrossRef]
  33. Hongfeng, Z. Dynamic economic dispatch based on improved differential evolution algorithm. Clust. Comput. 2018, 1–8. [Google Scholar] [CrossRef]
  34. Dubey, H.M.; Pandit, M.; Panigrahi, B. Hybrid flower pollination algorithm with time-varying fuzzy selection mechanism for wind integrated multi-objective dynamic economic dispatch. Renew. Energy 2015, 83, 188–202. [Google Scholar] [CrossRef]
  35. Elattar, E.E. A hybrid genetic algorithm and bacterial foraging approach for dynamic economic dispatch problem. Int. J. Electr. Power Energy Syst. 2015, 69, 18–26. [Google Scholar] [CrossRef]
  36. Nanjundappan, D. Hybrid weighted probabilistic neural network and biogeography based optimization for dynamic economic dispatch of integrated multiple-fuel and wind power plants. Int. J. Electr. Power Energy Syst. 2016, 77, 385–394. [Google Scholar]
  37. Xie, M.; Luo, W.; Cheng, P.; Ke, S.; Ji, X.; Liu, M. Multidisciplinary collaborative optimisation-based scenarios decoupling dynamic economic dispatch with wind power. IET Renew. Power Gener. 2018, 12, 727–734. [Google Scholar] [CrossRef]
  38. Basu, M. Quasi-oppositional group search optimization for multi-area dynamic economic dispatch. Int. J. Electr. Power Energy Syst. 2016, 78, 356–367. [Google Scholar] [CrossRef]
  39. Pattanaik, J.K.; Basu, M.; Dash, D.P. Improved real coded genetic algorithm for dynamic economic dispatch. J. Electr. Syst. Inf. Technol. 2018, 5, 349–362. [Google Scholar] [CrossRef]
  40. Pan, S.; Jian, J.; Yang, L. A hybrid MILP and IPM approach for dynamic economic dispatch with valve-point effects. Int. J. Electr. Power Energy Syst. 2018, 97, 290–298. [Google Scholar] [CrossRef] [Green Version]
  41. Somuah, C.; Khunaizi, N. Application of linear programming redispatch technique to dynamic generation allocation. IEEE Trans. Power Syst. 1990, 5, 20–26. [Google Scholar] [CrossRef]
  42. Lin, C.E.; Viviani, G. Hierarchical economic dispatch for piecewise quadratic cost functions. IEEE Trans. Power Appar. Syst. 1984, 1170–1175. [Google Scholar] [CrossRef]
  43. Irisarri, G.; Kimball, L.; Clements, K.; Bagchi, A.; Davis, P. Economic dispatch with network and ramping constraints via interior point methods. IEEE Trans. Power Syst. 1998, 13, 236–242. [Google Scholar] [CrossRef]
  44. Binetti, G.; Davoudi, A.; Naso, D.; Turchiano, B.; Lewis, F.L. A distributed auction-based algorithm for the nonconvex economic dispatch problem. IEEE Trans. Ind. Inform. 2013, 10, 1124–1132. [Google Scholar] [CrossRef]
  45. Atwa, Y.; El-Saadany, E.; Salama, M.; Seethapathy, R. Optimal renewable resources mix for distribution system energy loss minimization. IEEE Trans. Power Syst. 2009, 25, 360–370. [Google Scholar] [CrossRef]
  46. Ali, A.; Raisz, D.; Mahmoud, K.; Lehtonen, M. Optimal placement and sizing of uncertain PVs considering stochastic nature of PEVs. IEEE Trans. Sustain. Energy 2019. [Google Scholar] [CrossRef] [Green Version]
  47. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
  48. Abdel-Nasser, M.; Mahmoud, K. Accurate photovoltaic power forecasting models using deep LSTM-RNN. Neural Comput. Appl. 2019, 31, 1–14. [Google Scholar] [CrossRef]
  49. Kim, S.G.; Jung, J.Y.; Sim, M.K. A Two-Step Approach to Solar Power Generation Prediction Based on Weather Data Using Machine Learning. Sustainability 2019, 11, 1501. [Google Scholar] [CrossRef] [Green Version]
  50. Li, G.; Wang, H.; Zhang, S.; Xin, J.; Liu, H. Recurrent neural networks based photovoltaic power forecasting approach. Energies 2019, 12, 2538. [Google Scholar] [CrossRef] [Green Version]
  51. Yu, Y.; Cao, J.; Zhu, J. An LSTM Short-Term Solar Irradiance Forecasting Under Complicated Weather Conditions. IEEE Access 2019, 7, 145651–145666. [Google Scholar] [CrossRef]
  52. Jeon, B.K.; Lee, K.H.; Kim, E.J. Development of a Prediction Model of Solar Irradiances Using LSTM for Use in Building Predictive Control. J. Korean Sol. Energy Soc. 2019, 39, 41–52. [Google Scholar]
  53. Husein, M.; Chung, I.Y. Day-Ahead Solar Irradiance Forecasting for Microgrids Using a Long Short-Term Memory Recurrent Neural Network: A Deep Learning Approach. Energies 2019, 12, 1856. [Google Scholar] [CrossRef] [Green Version]
  54. Mirjalili, S. Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm. Knowl.-Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  55. 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]
  56. Wolpert, D.H.; Macready, W.G. No Free Lunch Theorems for Optimization. IEEE Trans. Evol. Comput. 1997, 1, 67. [Google Scholar] [CrossRef] [Green Version]
  57. Güvenç, U.; Sönmez, Y.; Duman, S.; Yörükeren, N. Combined economic and emission dispatch solution using gravitational search algorithm. Sci. Iran. 2012, 19, 1754–1762. [Google Scholar] [CrossRef] [Green Version]
  58. Sharifzadeh, M.; Sikinioti-Lock, A.; Shah, N. Machine-learning methods for integrated renewable power generation: A comparative study of artificial neural networks, support vector regression, and Gaussian Process Regression. Renew. Sustain. Energy Rev. 2019, 108, 513–538. [Google Scholar] [CrossRef]
  59. You, S.; Rasmussen, C.N. Generic modelling framework for economic analysis of battery systems. In Proceedings of the IET Renewable Power Generation Conference, Edinburgh, UK, 6–8 September 2011. [Google Scholar]
Figure 1. An overview of the structure of sustainable power systems.
Figure 1. An overview of the structure of sustainable power systems.
Sustainability 12 00576 g001
Figure 2. Structure of sustainable power systems.
Figure 2. Structure of sustainable power systems.
Sustainability 12 00576 g002
Figure 3. The block diagram of LSTM.
Figure 3. The block diagram of LSTM.
Sustainability 12 00576 g003
Figure 4. The predicted solar irradiance using LSTM for one day.
Figure 4. The predicted solar irradiance using LSTM for one day.
Sustainability 12 00576 g004
Figure 5. Demand load during the day.
Figure 5. Demand load during the day.
Sustainability 12 00576 g005
Figure 6. PV output power during clear and cloudy days.
Figure 6. PV output power during clear and cloudy days.
Sustainability 12 00576 g006
Figure 7. Generated power of the 40 thermal generators with case 1 at 50% penetration.
Figure 7. Generated power of the 40 thermal generators with case 1 at 50% penetration.
Sustainability 12 00576 g007
Figure 8. Charging and discharging of storage device with clear weather.
Figure 8. Charging and discharging of storage device with clear weather.
Sustainability 12 00576 g008
Figure 9. Convergence characteristics of the ISSA with case 1.
Figure 9. Convergence characteristics of the ISSA with case 1.
Sustainability 12 00576 g009
Figure 10. Generated power of the 40 thermal generators with case 2 at 50% penetration.
Figure 10. Generated power of the 40 thermal generators with case 2 at 50% penetration.
Sustainability 12 00576 g010
Figure 11. Charging and discharging of storage device with cloudy weather.
Figure 11. Charging and discharging of storage device with cloudy weather.
Sustainability 12 00576 g011
Figure 12. Convergence characteristics of the ISSA with case 2.
Figure 12. Convergence characteristics of the ISSA with case 2.
Sustainability 12 00576 g012
Table 1. The benchmark problems.
Table 1. The benchmark problems.
FunctionDimensionLimits f min
F 1 ( x ) = i = 1 n x i 2 30[−100,100]0
F 2 ( x ) = i = 1 n | x i | + i = 1 n i = 1 n | x i | 10[−10,10]0
F 3 ( x ) = i = 1 n ( j = 1 i x j ) 2 10[−100,100]0
F 4 ( x ) = max | x i | , 1 i n 10[−100,100]0
F 5 ( x ) = i = 1 n 1 [ 100 ( x i + 1 x i 2 ) 2 ] + ( x i 1 ) 2 10[−30,30]0
F 6 ( x ) = i = 1 n ( [ X i + 0.5 ] ) 2 10[−100,100]0
F 7 ( x ) = i = 1 n i x i 4 + r a n d o m 0 , 1 10[−1.28,1.28]0
F 8 ( x ) = i = 1 n x i sin ( | x i | ) 10[−500,500]0
F 9 ( x ) = i = 1 n [ x i 2 10 cos ( 2 π x i ) + 10 ] 10[−5.12,5.12]0
F 10 ( x ) = i = 1 n 20 exp ( 0.2 1 n i = 1 n x i 2 ) exp ( 1 n i = 1 n cos ( 2 π x i ) ) + 20 + e 100[−32,32]0
Table 2. Characteristics of the 40 thermal generators [57].
Table 2. Characteristics of the 40 thermal generators [57].
UnitCost CoefficientsPower Limits
a ($/MW2 h)b ($/MWh)c ($/h)e ($/h)f (rad/MW)Pmin (MW)Pmax (MW)
10.006906.7394.7051000.08436114
20.006906.7394.7051000.08436114
30.020287.07309.5401000.08460120
40.009428.18369.0301500.06380190
50.011405.35148.8901200.0774797
60.011428.05222.3301000.08468140
70.003578.03287.7102000.042110300
80.004926.99391.9802000.042135300
90.005736.60455.7602000.042135300
100.0060512.9722.8202000.042130300
110.0051512.9635.2002000.04294375
120.0056912.8654.6902000.04294375
130.0042112.5913.4003000.035125500
140.007528.841760.4003000.035125500
150.007528.841760.4003000.035125500
160.007528.841760.4003000.035125500
170.003137.97647.8503000.035220500
180.003137.95647.8503000.035220500
190.003137.97647.8503000.035242550
200.003137.97647.8503000.035242550
210.002986.63785.9603000.035254550
220.002986.63785.9603000.035254550
230.002846.66794.5303000.035254550
240.002846.66794.5303000.035254550
250.002777.10801.3203000.035254550
260.002777.10801.3203000.035254550
270.521243.331055.1001200.07710150
280.521243.331055.1001200.07710150
290.521243.331055.1001200.07710150
300.011405.35148.8901200.0774797
310.001606.43222.9201500.06360190
320.001606.43222.9201500.06360190
330.001606.43222.9201500.06360190
340.000108.95107.8702000.04290200
350.000108.62116.5802000.04290200
360.000108.62116.5802000.04290200
370.016105.88307.450800.09825110
380.016105.88307.450800.09825110
390.016105.88307.450800.09825110
400.003137.97647.8303000.035242550
Table 3. Comparing the performance of the ISSA, SSA, MFO, and MVO with ten benchmark problems.
Table 3. Comparing the performance of the ISSA, SSA, MFO, and MVO with ten benchmark problems.
FunctionOptimization Method
ISSASSAMFOMVO
AveStdAveStdAveStdAveStd
F10.000000.000000.000010.000000.000140.000171.19690.1407
F20.000010.000000.000000.000000.000000.000000.03570.0128
F30.000000.000000.00000.00000.000000.000010.11520.0771
F40.000020.000000.000020.000000.60321.36820.09270.035
F54.165003.16217.21122.43678.15027.718488.0909125.6254
F60.000000.000000.00000.00000.000000.000000.01560.0043
F70.001300.00140.00480.003000.004900.002600.00350.0018
F8−2955.8184.4003−2746.9239.6771−3305.9243.2327−3006337.0665
F90.298500.671514.92444.547417.41186.735914.63344.7858
F106.38441.31057.49991.551319.84390.22597.17376.7536
Table 4. Comparison with related methods.
Table 4. Comparison with related methods.
MethodMLRBRTANNLSTM
RMSE of dataset1384.8951494.4633377.07282.15
RMSE of dataset2329.11416.212348.931136.87
Recurrent?×××
Can it remember?×××
Table 5. Total fuel cost for case 1 (clear weather).
Table 5. Total fuel cost for case 1 (clear weather).
PenetrationFuel Cost*106 ($)
ISSASSAMFOMVO
0%2.76872.84513.09042.8945
10%2.70172.81622.91532.8945
20%2.67902.77832.85182.8945
30%2.57142.73852.79492.7446
40%2.51552.61312.76142.7446
50%2.42282.53952.58922.7446
Table 6. Total fuel cost for case 2 (cloudy weather).
Table 6. Total fuel cost for case 2 (cloudy weather).
PenetrationFuel Cost*106 ($)
ISSASSAMFOMVO
0%2.76872.84513.09042.8945
10%2.72612.82982.88242.8824
20%2.68492.78612.79452.7945
30%2.67102.78612.78232.7823
40%2.52442.69052.71712.7171
50%2.48772.64952.66412.6641

Share and Cite

MDPI and ACS Style

Mahmoud, K.; Abdel-Nasser, M.; Mustafa, E.; M. Ali, Z. Improved Salp–Swarm Optimizer and Accurate Forecasting Model for Dynamic Economic Dispatch in Sustainable Power Systems. Sustainability 2020, 12, 576. https://doi.org/10.3390/su12020576

AMA Style

Mahmoud K, Abdel-Nasser M, Mustafa E, M. Ali Z. Improved Salp–Swarm Optimizer and Accurate Forecasting Model for Dynamic Economic Dispatch in Sustainable Power Systems. Sustainability. 2020; 12(2):576. https://doi.org/10.3390/su12020576

Chicago/Turabian Style

Mahmoud, Karar, Mohamed Abdel-Nasser, Eman Mustafa, and Ziad M. Ali. 2020. "Improved Salp–Swarm Optimizer and Accurate Forecasting Model for Dynamic Economic Dispatch in Sustainable Power Systems" Sustainability 12, no. 2: 576. https://doi.org/10.3390/su12020576

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