Next Article in Journal
A Novel Contact Stiffness Model for Grinding Joint Surface Based on the Generalized Ubiquitiformal Sierpinski Carpet Theory
Next Article in Special Issue
The Regulation of Superconducting Magnetic Energy Storages with a Neural-Tuned Fractional Order PID Controller Based on Brain Emotional Learning
Previous Article in Journal
A Novel Stochastic Model for Human Norovirus Dynamics: Vaccination Impact with Lévy Noise
Previous Article in Special Issue
A High-Performance Fractional Order Controller Based on Chaotic Manta-Ray Foraging and Artificial Ecosystem-Based Optimization Algorithms Applied to Dual Active Bridge Converter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimizing Economic Dispatch with Renewable Energy and Natural Gas Using Fractional-Order Fish Migration Algorithm

by
Abdallah Aldosary
Computer Engineering Department, College of Engineering, Prince Sattam bin Abdulaziz University, Wadi Addwasir 11991, Ar Riyadh, Saudi Arabia
Fractal Fract. 2024, 8(6), 350; https://doi.org/10.3390/fractalfract8060350
Submission received: 7 May 2024 / Revised: 1 June 2024 / Accepted: 4 June 2024 / Published: 12 June 2024

Abstract

:
This work presents a model for solving the Economic-Environmental Dispatch (EED) challenge, which addresses the integration of thermal, renewable energy schemes, and natural gas (NG) units, that consider both toxin emission and fuel costs as its primary objectives. Three cases are examined using the IEEE 30-bus system, where thermal units (TUs) are replaced with NGs to minimize toxin emissions and fuel costs. The system constraints include equality and inequality conditions. A detailed modeling of NGs is performed, which also incorporates the pressure pipelines and the flow velocity of gas as procedure limitations. To obtain Pareto optimal solutions for fuel costs and emissions, three optimization algorithms, namely Fractional-Order Fish Migration Optimization (FOFMO), Coati Optimization Algorithm (COA), and Non-Dominated Sorting Genetic Algorithm (NSGA-II) are employed. Three cases are investigated to validate the effectiveness of the proposed model when applied to the IEEE 30-bus system with the integration of renewable energy sources (RESs) and natural gas units. The results from Case III, where NGs are installed in place of two thermal units (TUs), demonstrate that the economic dispatching approach presented in this study significantly reduces emission levels to 0.4232 t/h and achieves a lower fuel cost of 796.478 USD/MWh. Furthermore, the findings indicate that FOFMO outperforms COA and NSGA-II in effectively addressing the EED problem.

1. Introduction

As a reaction to diverse environmental, social, economic, technical, and political concerns, there has been a notable growth in the capacity of renewable energy systems (RESs) within contemporary power grids [1,2]. Consequently, the environmental economic dispatch (EED) problem has emerged as a crucial optimization challenge in power system operation and planning. This is particularly significant due to the growing apprehension about global climate change and environmental contamination associated with traditional fossil fuel-based electricity production [3,4]. In the realm of power generation, the traditional EED issue revolves around determining the optimal output levels of generating units to meet the load demand while minimizing operational expenses. Simultaneously, the problem must adhere to system constraints and address environmental concerns to minimize pollution from thermal power plants. Mathematically, the EED issue is non-convex [5,6], involving inconsistent objectives and non-linear restrictions derived from power flow limitations, grid compliance requirements, valve point impacts, and prohibited operating zones for units. Efficient algorithms are necessary to find the optimal trade-off between environmental and economic considerations. Moreover, the integration of RESs into electricity systems further adds complexity to the problem [7,8].
In the last decades, there has been a growing trend of integrating cost and emissions considerations in the development and management of electric systems, leading to the emergence of the EED challenge [9]. The objective of EED is to determine the optimal allocation of power generation among different generating units that are both cost-effective and environmentally responsible. However, EED poses significant difficulties as the fuel cost function of generators often exhibits nonlinearity and non-convexity. To tackle this challenge, various optimization methods have been proposed in the literature. These approaches can be grouped into three classifications:
  • Conventional techniques that address convex problems using established approaches.
  • Unconventional methods that focus on practical and non-convex problems, offering alternative solutions.
  • Hybrid techniques overcome the limitations of individual nonconventional algorithms by combining their strengths to solve complex problems. These hybrid methods have demonstrated success in obtaining globally optimal solutions for EED problems under various constraints.
Iterative or conventional techniques are commonly used as conventional approaches in addressing the EED problem. These techniques include the linear programming approach [10], gradient iterative approach [11], quadratic programming approach [12], Lagrange’s method (LM) [13,14], hybrid mixed-integer linear programming and the interior point technique [15], and Newton–Raphson approach [16]. In a study [17], the EED issue has been reformulated into a quadratic programming problem, which can be further expressed as a semidefinite programming formulation. Through the utilization of convex iteration and branch-and-bound techniques, this problem can be solved iteratively. However, the efficiency of these techniques heavily relies on the initial conditions, often requiring a large number of iterations to achieve satisfactory results, and they may struggle to converge to a global solution. Consequently, these methods are quite intricate, involving computational complexities and complex mathematical expressions. As a result, unconventional methods have emerged as effective alternatives for addressing the optimal EED problem.
In contrast, researchers have proposed various metaheuristic methods as alternative approaches to tackle this challenge. The advantages of metaheuristic algorithms have been well-established in addressing complex optimization problems [18]. In a study [19], the authors introduced a varied version of the sine cosine optimization (SCO) algorithm to solve the optimal EED problem. They made improvements in two key aspects. Firstly, an enhanced elite leadership technique was incorporated throughout the particle position update phase, resulting in improved search capabilities of the algorithm. Secondly, a combination of crossover and optimal choices was utilized to block the algorithm from being trapped in local optima. Additionally, a dimension-by-dimension variation technique was implemented to enhance the optimization accuracy of SCO and increase population diversity. The presented algorithm was tested on two systems: the IEEE 30-bus and the 10-unit test schemes [19]. Another approach, inspired by human behavior during search and rescue operations, was employed in addressing the EED problem using the search and rescue optimizer [20]. Furthermore, a squirrel search optimization method with function was projected to solve the EED issue [21]. The weighted sum method was utilized to convert the multi-objective function into a single-objective function, demonstrating the efficiency of the algorithm in resolving the multi-objective power system optimization issue through experimental tests. Additionally, in a different study [22], the Monte Carlo method was applied to a power dispatching strategy for the baseline case. This strategy prioritized wind and solar PV generators to assist hydropower plants in meeting demand while considering the loss of load probability (LOLP) and expected demand not supplied (EDNS).
The genetic algorithm (GA) has been projected as a solution to the optimal EED problem [23]. However, dependent on the size of the scheme under analysis, the computational time required for execution may lead to suboptimal solutions. In a study [24], the Moth Swarm Algorithm was tested on two EED test systems comprising a combination of thermal and PV plants across a 24 h period, while considering spinning reserve allocation. To address the coupled EED problem, both GA and particle swarm optimization (PSO) were proposed [25]. These methods were implemented on a self-sufficient power plant in Pakistan, considering varying load requirements. The results obtained confirmed that PSO outperformed GA. In another study [26], the authors introduced a multi-objective learning backtracking search approach for the EED issue. This method integrated leader selection and guidance as learning methodologies to enhance the uniformity and diversity of the Pareto frontier. Furthermore, for the coupled EED problem, a multi-objective membrane search approach was proposed [27]. The evaluation of this algorithm across multiple test cases demonstrated its high spatial detection capability and ability to generate improved solutions for the EED problem.
In order to address the EED issue, a varied version of the shuffle frog-jump approach (SFLA) was proposed [28]. This modified version incorporates a local search method based on the inertia equation, which enhances the performance of the original SFLA. Furthermore, the integration of crossover and mutation operators from the GA improves the global search process. In a separate study [29], the Shark Optimization Algorithm (SOA) was developed to allocate dispersed generations in distribution systems. The objective was to minimize losses while simultaneously improving bus voltage values and stability. However, this analysis did not consider emissions reduction. To tackle the EED problem, a hybrid approach combining a modified version of the artificial bee colony technique and a Non-Dominated Sorting Genetic Algorithm-II (NSGA-II) was suggested [30]. This hybrid method was employed in three schemes with 10, 20, and 40 producing units, demonstrating its effectiveness. Another proposed hybrid strategy, presented in [31], combines the particle diversity of PSO with the quick convergence of differential evolution (DE) techniques. The crossover operators from GA were also incorporated into this hybrid method. A parameter adaptive control technique was manipulated to modernize the crossover probability and improve optimization findings. Additionally, the issue of multi-objective optimization was addressed by introducing a penalty factor.
In order to predict future load conditions, a functional link artificial neural network (ANN) model based on single-layer black widow optimization was presented by the authors in [32]. Additionally, a nondominated sorting multi-objective instructional learning-based optimizer was presented to address the EED issue in a grid with a PV source, taking into account the anticipated future load. The effectiveness of the recommended method was assessed under two distinct circumstances with varying sun radiation. In a study [33], an ANN controller was utilized to install a photovoltaic distributed generator in an IEEE 33-bus distribution feeder. However, it should be noted that the ANN-based technology provided was only suitable for a single unit. Furthermore, various optimization approaches have been suggested to tackle the EED issue. These methods include the flower pollination algorithm [34], pigeon-inspired optimizer [35], improved bare-bone multi-objective PSO [36], lightning flash algorithm [37], chameleon swarm algorithm [38], niching penalized chimp algorithm [39], multi-verse optimization algorithm [40], bat algorithm [41], teaching–learning-based algorithm [42], hybrid crow search–JAYA algorithm [43], simulated annealing algorithm [44], improved marine predators algorithm [45], osprey optimization algorithm [9], quantum-behaved bat algorithm [46], hybrid differential evaluation and crow search approach [47], adaptive Hooke and Jeeves algorithm [48], mantis search algorithm [49], and hybrid heap-based and jellyfish search approach [50].
In this research, the conventional IEEE 30-bus system was modified to incorporate a restricted set of thermal units (TUs), allowing for the integration of solar PV, natural gas (NG), and wind generation units. The study extensively examines the stochastic characteristics of RESs such as wind and solar. Probability density functions (PDFs) including Weibull and lognormal distributions are utilized to capture their uncertainty. To account for the variability and intermittency of these RESs, the proposed cost model includes penalty costs for underestimation and reserve costs for overestimation. Subsequently, a restricted optimization issue is developed to minimize the EED issue while satisfying various constraints. These constraints encompass power flow equality and inequality limits, NG limitations, and prohibited operating zones (POZs) bounds. To solve this EED issue, three strategies are employed: Fractional-Order Fish Migration Optimization (FOFMO) [51], Coati Optimization Algorithm (COA) [52], and NSGA-II [53]. It is worth noting that most multi-objective optimization techniques are typically tailored for unconstrained optimization problems. The key contributions of this study can be summed up as follows:
  • Formulation of the EED issue considers a combination of thermal, solar, wind, and NG units.
  • Detailed stochastic study of RESs using appropriate PDFs.
  • Inclusion of all system constraints, including security, power flow conditions (equality and inequality), POZs, and NG constraints, are formulated in the EED issue.
  • Utilization of optimization techniques, specifically FOFMO, COA, and NSGA-II are employed.
  • Comparative analysis of the solutions attained from the utilized optimization techniques.
The rest of the paper is organized as follows:
  • Section 2 introduces the system configuration being studied.
  • Section 3 includes a mathematical analysis of the EED problem, taking into account the incorporation of RESs and NG. The interpretation of the EED issue, system constraints, and the optimization methods employed in this work are also presented in this section.
  • Section 4 depicts and elucidates the obtained illustrations and results.
  • Finally, in Section 5, conclusions drawn from the study and suggestions for future research are provided.

2. System Investigated

The primary objective of power grids is to optimize generator output to meet demand requirements while minimizing fuel costs, with emissions not being a factor. Presently, every nation is implementing various strategies and innovative approaches to reduce pollutants and safeguard the atmosphere in alignment with national and global environmental conservation standards. Furthermore, the minimization of fuel costs for generation units is a critical factor in addressing the EED problem and meeting load requirements. Simultaneously, reducing emissions (referred to as the environmental objective) holds great importance in mitigating climate change and environmental pollution. Traditionally, the EED problem solely focused on the economic objective concern. However, due to the urgency of combating climate change, greater attention has been given to adjusting greenhouse gas productions. Hence, it is highly significant to address the EED challenge, which involves balancing both economic and environmental goals. To realize this, the integration of clean and cost-effective RESs into electrical networks becomes a crucial enabler for minimizing environmental emissions.
This study examines a modified version of the IEEE standard 30-bus scheme, which incorporates both TUs and RESs. The system model, illustrated in Figure 1, includes three TUs located at buses 1, 2, and 8, as well as three distinct RESs, namely wind and two PV systems situated at buses 5, 11, and 13, respectively. The specifications of the scheme under investigation can be found in Table 1.
To achieve the simultaneous reduction of emissions and fuel cost, three distinct cases are examined. In each case, three RESs comprising wind and two PV units are implemented at buses 5, 11, and 13, respectively. Three optimization methods of FOFMO, COA, and NSGA-II are applied to the problem under investigation. Figure 2 illustrates the methodology for the presented techniques. The breakdown of these cases is as follows:
  • Case I: Utilizing three TUs and three RESs [55].
  • Case II: Changing the fuel of TU at bus 1 into NG.
  • Case III: Changing the fuel of TUs at buses 2 and 8 into NGs.

3. Formulation of the Optimization Problem

3.1. Objective Functions

The EED issue to be addressed is one in which the first objective function reflects the economic elements ( C t o t ) and the second represents the environmental implications ( E t o t ) , as described in Equation (1).
Minimize E E D = m i n ( C t o t , E t o t )

3.1.1. Fuel Costs

The overall cost of produced electricity is the sum of the costs of TUs, RESs, and NGs (if NGs are included), as given in Equation (2).
C t o t = C t o t P T U + C t o t P R E S s + α × C t o t ( P N G )
where C t o t P T U denotes the total TUs cost, C t o t P R E S denotes the total RESs cost, and C t o t P N G denotes the total NGs cost, where α =0 with no NGs involved and α = 1 if NGs are involved.

Fuel Cost Study of TUs

The TU cost should take into account the quantity of steam to the turbine’s blades, as well as any rapid changes in the valve’s condition. Steam in these facilities is regulated by valves that direct the turbine via an independent set of nozzles. To attain maximum efficiency, this set of nozzles has to be operated at full power [56]. For obligatory production, these valves open progressively. Equation (3) represents the cost function of TUs [57].
C t o t P T U = i = 1 N T U a T U i + b T U i P T U i + c T U i P T U i 2 + d T U i sin e T U i × P T U i m i n P T U i
where a T U i , b T U i , and c T U i denote the cost factors of the ith thermal generator unit ( P T U i ) . d T U i and e T U i represent the valve point impacts. P T U i m i n denotes the minimum power of P T U i during the generator operation. These factors are described in Table 2.

Fuel Cost Study of RESs

The cost of the RESs, in USD/MWh, represents the total costs of WPs C t o t P W P and PVs C t o t P P V which can be expressed as in Equation (4):
C t o t P R E S = C t o t P W P + C t o t P P V
However, each RES possesses a distinct cost function. In order to address the intermittent nature of RESs, two approaches can be employed. First, standby power units (SPUs) can be installed to compensate for the deficit in generated power compared to the scheduled amount. Second, energy storage (ES) can be fitted to store any surplus power generated beyond the planned capacity.
To establish cost expressions, casual data for wind speed and solar irradiance can be fitted using the Weibull and lognormal distributions respectively. These distributions are utilized in the subsequent subsections to formulate the cost expressions.
A.
Cost estimation of wind plants
The cost calculation for WPs incorporates various factors, including the initial investment costs and the expenses associated with standby and energy storage (ES) units. The direct cost, denoted as C d W P , accounts for the initial investment, operational, and maintenance costs. Equation (5) defines this cost relationship as follows:
C d W P = K d W P P W P s c h
where K d W P and P W s c h denote the direct cost factor and the scheduled power of the WPs, respectively. When the actual power provided by the turbines drops below the scheduled power, the system uses reserve units to satisfy demand. The reserve capacity C r W P cost is expressed in Equation (6).
C r W P = K r W P P W P s c h P W P a c t = K r W P 0 P W a c t P W P s c h p W P f w p W P d p w
The cost factor for the standby elements is denoted as K r W P , while P W P a c t represents the actual power delivered by the wind turbines. Similarly, when the actual power, P W P a c t , exceeds the P W P s c h , the cost related to the storage units can be expressed by Equation (7) as follows:
C s W P = K s W P P W P a c t P W P s c h = K s W P P W P s c h P W P r p W P P W P s c h f w p W P d p w
where the variable P W P r , f w p W P signifies the rated wind power and the PDF of the wind speed. Typically, the PDF of the wind power is modeled using the widely recognized Weibull distribution, which effectively captures the arbitrary occurrence of different wind velocity (v) levels [58,59]. In Figure 3a, the Weibull PDF of the wind velocity statistics is depicted, generated through 8000 Monte Carlo situations that consider the scale and shape factors (represented as α and β). In this case, α is set to 9, and β is set to 2. The probability f v v of a specific (v) can be described using Equation (8) [58]:
f v v = β α v α β 1 e v α β
The cost factors associated with wind power are presented in Appendix A. The equation describing the power output of wind turbines, which is contingent upon the wind speed, is provided as Equation (9) [59]:
p W P = 0 v o u t v v i n P W P r v v i n v r v i n v i n v v r P W P r v r v v o u t
where v i n ,   v r , and v o u t indicate the cut-in speed, nominal speed, and cutout speed of the WPs, respectively. The wind power prospect can be stated as in Equation (10):
f w p W P = β v r v i n α β × P W r v i n + p W P P W P r v r v i n β 1 × e x p v i n + p W P P W P r v r v i n α β
Lastly, the total cost of the WP is stated in Equation (11):
C t o t W P = C d W P + C r W P + C s W P
B.
Cost Estimation of the PV Plant
Similarly, the overall cost equation for the PV plant can be formulated following the same approach employed to estimate the cost equations for the WP.
The direct cost, denoted as C d P V , signifies the initial investment, operational, and maintenance expenses associated with PVs. Equation (12) defines this cost relationship as follows:
C d P V = K d P V P P V s c h
In this equation, K d P V and P P V s c h correspond to the direct cost factor and the scheduled power of the PV, respectively. When the actual power ( P P V a c t ) drops below the scheduled power ( P P V s c h ) , the system uses reserve units to satisfy demand [60]. The reserve capacity C r P V cost is expressed in Equation (13).
C r P V = K r P V P P V s c h P P V a c t = K r P V P P V s c h p P V × f P V p P V
K r P V represents the cost factor of the standby elements. If P P V a c t exceeds P P V s c h , the cost related with the storage units can be expressed by Equation (14) as follows
C s P V = K s P V P P V a c t P P V s c h = K s P V p P V P P V s c h × f P V p P V
The power provided by the backup and storage elements is proportional to the solar irradiance PDF G . The PDF of G is fitted using a lognormal model [60,61,62], as shown in Figure 3b. Equation (15) gives the probability of G with lognormal fit values of μ = 5.61 and σ = 0.61 ; therefore:
f P V G = 1 G σ 2 π e x p l n G μ 2 2 σ 2 ,           G > 0
The accessible power from the PV scheme can be revealed as follows:
p P V G = P P V r G 2 G s t d ,         0 < G < R c P P V r G G s t d ,                   G R c
where the standard solar irradiation G s t d is equivalent to 1000 W/m2. Through the operation irradiance R c is adjust as 120 W/m2. P P V r denotes the nominal output power of the PV units. The cost factors of PV can be explained in Appendix A.
Ultimately, the total cost of the PV is stated as follows:
C t o t P V = C d P V + C r P V + C s P V

Fuel Cost Study of NGs

In this study, NGs are used rather than TUs to reduce fuel costs and pollution emissions as much as feasible. The overall cost of NGs is comprised of many spending components: original cost, operation and maintenance cost, and fuel cost. Therefore, the overall cost of the NGs C t o t P N G is given in Equation (18) [63,64], thus:
  C t o t P N G = i = 1 N N G g N G i × P N G i + i = 1 N N G a N G i + b N G i P N G i + c N G i P N G i 2
where a N G i , b N G i , and c N G i signify the cost factors of the ith NG ( P N G i ) , and g N G i indicates the factor of the initial and operating costs of NGs. These factors are detailed in Table 3.
To summarize, C t o t of the complete scheme can be written as in Equation (19).
C t o t = i = 1 N T U a T U i + b T U i P T U i + c T U i P T U i 2 + d T U i × sin e T U i × P T U i m i n P T U i   + C d W P + C r W P + C s W P + C d P V + C r P V + C s P V + i = 1 N N G g N G i P N G i + i = 1 N N G a N G i + b N G i P N G i + c N G i P N G i 2

3.1.2. Emission Levels

Equation (20) describes the total emissions E t o t resulting from thermal and natural gas units, as the RESs do not emit any emissions. Only the emissions of TUs and NGs units are accounted for in this presentation.
E t o t = i = 1 N T U E T U i + i = 1 N N G E N G i
where E T U i and E N G i denote the overall emissions of TU and NG, respectively.

Emission Study of TUs

Equation (21) provides a method to calculate the total emission E t o t ( T U ) in tons per hour (t/h) of environmentally harmful gases such as SOx and NOx, which are directly linked to the generated output power.
E t o t ( T U ) = i = 1 N T U φ T U i + ( ψ T U i   × P T U i ) + ( ω T U i × P T U i 2 ) + τ T U i × e ξ T U i P T U i
where φ T U i , ψ T U i , ω T U i , τ T U i , and ξ T U i denote the factors of emission levels correlated with the ith T U . These factors can be identified in Table 2.

Emission Study of NGs

Equation (22) [66] describes the total emissions resulting from NGs, which is equivalent to the total emissions resulting from the TUs.
E t o t ( N G s ) = i = 1 N N G φ N G i + ( ψ N G i × P N G i ) + ( ω N G i × P N G i 2 )
where φ N G i , ψ N G i , and ω N G i denote the factors of emission levels correlated with the ith N G . These factors can be identified in Table 3.

3.2. Limitations

When solving the objective function for any configuration, it is crucial to take into account the following key restrictions.

3.2.1. Power Balance Constraints

The power balance constraints pertain to both active and reactive powers which can be described as follows:
P G U = P L i + P L o s s i
Q G U = Q L i + Q L o s s i

3.2.2. Active and Reactive Powers Limits

The operational restrictions of all utilized elements of the active and reactive powers, involving thermal, wind, solar, and NGs, are expressed by Equations (25) to (32) respectively.
P T U i m i n P T U i P T U i m a x           i N T U
P W P m i n P W P P W P m a x
P P V m i n P P V P P V m a x
P N G m i n P N G P N G m a x
Q T U i m i n Q T U i Q T U i m a x           i N T U
Q W P m i n Q W P Q W P m a x
Q P V m i n Q P V Q P V m a x
Q N G m i n Q N G Q N G m a x

3.2.3. Prohibited Operating Zones Boundaries

In order to account for physical limitations experienced by TUs, such as stillness in shaft bearings or failures in pumps and boilers, Power Operating Zones (POZs) can be permitted within specific operating districts. These POZs result in a non-continuous operation of TUs. Equation (33) provides a description of the POZs.
P T U i m i n P O Z , j P O Z T U i j P T U i m a x P O Z , j
where P T U i m i n P O Z , j and P T U i m a x P O Z , j denote the min and max restrictions of the jth POZ of the ith TUs, respectively.

3.2.4. Security Limitations

Equations (34) and (35) depict the voltage security constrictions for the alternators’ bus and the load bus voltage, respectively. Additionally, Equation (36) outlines the limitations on the capacity of the branches.
V G i m i n V G i V G i m a x           i N G
V L i m i n V L i V L i m a x           i N L
S L i S L i m a x                           i n l
where V G i denotes the voltage of the ith generator bus, while V L i represents the voltage of the ith load bus. N G , N L , and n l signify the count of alternator, load, and branches buses in the utility, respectively.
Furthermore, two additional parameters, network power loss ( P l o s s ) and voltage deviation (VD), are considered as network constraints. They can be computed using Equations (37) and (38), respectively.
P l o s s = q = 1 n l G q ( i j ) V i 2 + V j 2 2 V i V j cos δ i j
V D = p = 1 N L V L p 1
where G q ( i j ) and δ i j indicate the transconductance of branch q connected to bus i and bus j and the phase difference between δ i and δ j , respectively.

3.2.5. Natural Gas Boundaries

To be more precise, it is necessary to examine the impact of the required volume of natural gas ( V N G ) on the remaining loads within the natural gas supply network [67]. The specific volume of natural gas needed to generate a particular amount of electrical power can be determined using Equation (39).
V N G i = 0.278 η N G i H H V p N G i
The efficiency of the gas turbine is denoted as η N G i , and HHV signifies the high heat value of NG. The characteristics of the NG, including its distance, diameter, material kind, pressure, and mean flow rate ( q ~ x y , t ), play a crucial role prior to bringing the NG to a gas turbine. This relationship can be expressed mathematically through Equation (40).
q ~ x y , t = ± C N G T O p o η p p x 2 p y 2 D p 5 S N G Z T a v L p f
q ~ x y , t = q x y , t i n + q x y , t o u t 2
where C N G signifies the constant factor, T o is the typical temperature, p o indicates the absolute pressure under atmospheric environment, and p x and p y signify the absolute upstream pressure and absolute downstream pressure, respectively. D P signifies the inner diameter of the pipe in (mm), S N G is the specific gravity of NG, and Z is the factor that captures the contrast between the actual state of the gas and its ideal state. T a v represents the average temperature. L P refers to the span of the pipeline in meters, while f represents the hydraulic friction factor.
When the flow speed surpasses 20 m/s, the movement of dust atoms can have detrimental effects on cooking purposes and may even cause erosion on the internal surface of the pipeline. Therefore, the flow velocity ( U N G ) in (m/s) can be determined using Equation (42).
U N G = 353 × q N G × p o D p 2 p p 2 3730   f   L p q ~ x y , t 2 D p 5
The pipeline’s flow velocity and pressure are interrelated and act as constraints within the system during the transition from thermal to gas generation units. Both the U N G and pipeline pressure ( p p ) need to be taken into account [66,67].
Equations (43) and (44) define the equality restrictions for the line pack of pipeline x-y at hour t ( L x y , t ). Equation (45) specifies that the initial and final estimates of the line pack are equal. Additionally, Equations (46) to (49) describe the inequality restrictions pertaining to the flow speed, pressure, flow rate, and air compressor limitations, respectively.
L x y ,     t = G x y p x ,   t + p y ,   t 2
L x y ,     t = L x y ,     t 1 + q x y ,   t q y ,   t
t q x y , t i n = t q x y , t o u t
U N G 20   m / sec
p p i p e m i n p p i p e p p i p e m a x
q p i p e m i n q p i p e q p i p e m a x
p x ,   t Γ C   p y ,     t

3.3. Optimization Techniques

3.3.1. Fish Migration Optimization (FMO) Algorithm

In 2010, the FMO approach was introduced as a swarm intelligence technique [68,69]. It emulates the growth, migration procedures, and predation approach observed in fish biology. A notable distinction between FMO and other meta-heuristic approaches is that its optimization formularies were derived from insights provided by biologists. In comparison to PSO, FMO demonstrates superior accuracy and reasonable time consumption. However, when it comes to low-dimensional complex functions, FMO’s optimization performance is not particularly outstanding. This is primarily because the integer order speed update in FMO tends to overlook the optimal solution, resulting in suboptimal optimization outcomes. On the other hand, the fractional order speed update in FMO employs fractions to adjust the step size and incorporates historical speed information, enabling the algorithm to yield more precise results.
Each species in nature possesses unique survival strategies, methods of predation, and reproductive mechanisms. Furthermore, they must adapt to challenging environments, as they constantly face the risk of being captured by natural predators. Not every fish can successfully reach adulthood. Biologists have observed that fish swim in water for various reasons. Drawing inspiration from the grayling fish, the FMO algorithm is introduced, leveraging its characteristics. The life cycle of the grayling fish is divided into five stages, which serve as the basis for the procedure.
  • Stage 0+: Newborn and young (age ranging from 0 to 1 year).
  • Stage 1+: Juvenile (age ranging from 1 to 2 years).
  • Stage 2+: Sub-adult (age ranging from 2 to 3 years).
  • Stage 3+: Adult (age ranging from 3 to 4 years).
  • Stage 4+: Adult (age ranging from 4 to 5 years).
Each stage of the fish’s life cycle is associated with distinct fecundity rates. The life cycle graph of the grayling, depicted in Figure 4, illustrates this. Within the graph, F2, F3, and F4 represent the fecundity rates corresponding to stages 2+, 3+, and 4+ respectively.
The FMO algorithm is specifically developed to optimize processes through two main components: the swim and the migration processes. The swim procedure simulates the behavior of graylings as they navigate and grow in the water while seeking food sources. During this procedure, energy consumption aligns with the fish’s change, as defined by Equation (50).
E r ,   d = r a n d . E
E r ,   d represents the energy expended in the d dimension, while rand denotes a randomly generated number. E stands as a constant that states the maximum energy consumption in a single dimension. For this particular procedure, we assign a value of 2 to E. The operative correlation among the remoteness traveled and energy consumption can be expressed as:
d i s o f f s e t ,   d = E r ,   d . U s , d a + b . ( U s , d ) x
The moving distance is represented by disoffset; d, while U s , d denotes the swimming velocity in the d dimension. Constants a, b, and x can be involved, where a represents the standard metabolic rate, b is a scaling constant, and x is the velocity exponent. The values for a, b, and x are documented as 2.25, 36.2, and 2.23, respectively. To calculate the new position, the value of every dimension is updated using Equation (52), where offset represents the moving space of the particle.
p n e w = p o l d + p o f f s e t
If the fitness value of p n e w is superior to that of p o l d , p o l d will be revised using p n e w according to Equation (52). At the same time, the speed will be adjusted using Equation (53).
U s = 2 . U s
As the fish reach maturity, a portion of them return to their birthplace for the purpose of reproducing offspring, which is known as the migration process. This migration process occurs exclusively in stages 2+, 3+, and 4+, as fish in stages 0+ and 1+ are unable to reproduce. The fecundity rates for these three stages are 5%, 10%, and 100% respectively. When fish discover a new potential location, the coordinates are renewed using Equation (54).
p n e w = ( d m a x d m i n ) . r a n d + d m i n
d m a x and d m i n represent the upper and lower limits, respectively.
To determine the fitness value of another candidate, the velocity will be adjusted as follows:
U = π . U s ,                           F ( P ) < F ( P b e s t ) U s ,                                                     o t h e r w i s e
where U s represents the initial velocity.

Fractional Calculus (FC)

FC is an extension of the traditional perception of calculus, encompassing fractional derivatives and fractional integrals, much like classical calculus [70,71]. However, the key distinction lies in the order of derivatives and integrals: while classical calculus deals exclusively with integer orders, FC allows for fractional orders. This fundamental difference enables FC to precisely capture the memory and inherited traits of diverse constituents and their evolutionary procedures. Since the inception of FC, its associated philosophy has found successful applications across numerous fields.
The FC, originating from classical calculus with a history spanning over 300 years, had limited advancements until modern times caused by the emergence of employed disciplines like fluid mechanics, cybernetics, and biology. However, its practical significance gradually became apparent, leading to increased scholarly interest in FC. FC is a branch of mathematical analysis that offers several advantages. Firstly, it provides a mathematical perspective on the inevitability of historical development. Additionally, FC is known for its memorability. Moreover, compared to nonlinear models, the concise expression of the FC model allows for a better fit with the real world. There are three definitions of fractional derivatives (FDs), each associated with an order α.
The first definition is the Riemann–Liouville FD, and it is described as follows [70]:
R L   D α , t α f ( x ) = d n f d t n ,                                                                                                                                                       α = n , 1 Γ ( n α )   d n d x n α x x t n α 1 f n t ,                       n 1 < α < n
The second definition is the Caputo FD, and it is described as follows [70]:
C   D α , t α f ( x ) = d n f d t n ,                                                                                                                                                 α = n , 1 Γ ( n α )   α x x t n α 1 f n t ,                       n 1 < α < n
The third definition is the Grünwald–Letnikov FD, and it is described as follows [70]:
G L   D α f ( x ) = lim h 0 1 h α k = 0 + 1 k Γ ( α + 1 ) f ( x k h ) Γ ( k + 1 ) Γ ( α k + 1 )
To approximate the discrete time situation using Grünwald–Letnikov FD, the following expression can be used (referred to as Equation (59)).
G L   D α f t = 1 T α k = 0 r 1 k Γ ( α + 1 ) f ( t k T ) Γ ( k + 1 ) Γ ( α k + 1 )
The time increment is denoted as T, and the truncation order as r in the approximation of the discrete time situation using the Grünwald–Letnikov FD (referring to Equation (59)). The Gamma function, denoted as Γ (x), is present in all three FD definitions mentioned above. The three FC definitions mentioned above exhibit the following properties:
(a)
Linearity D α α f x + b g ( x ) = α D α f x + b D α g ( x )
(b)
The index theory D α + β f x = D α D β f x
(c)
Generalized Leibniz rule D α f x . g ( x ) = i = 0 D i f x . D α i g x
where a and b are constants.

Proposed Fractional-Order Fish Migration Optimization Algorithm

This section describes in detail a novel algorithm known as Fractional-Order FMO (FOFMO). The FOFMO algorithm combines the FMO method with the idea of the Grünwald–Letnikov FD. There are two major differences between the FMO technique and the FOFMO technique. First, the FOFMO employs a new technique for fractional-order velocity. Second, the FOFMO produces new offspring placements using the global optimal particles.
A.
Fractional-Order Velocity
The velocity update equations in the FMO algorithm are given by (53) and (55). While the FMO procedure excels in penetrating for global optimum solutions, it exhibits a limited exploitation ability as it requires considerable time for exploration. To enhance the exploitation capability, the FC is employed for velocity updates. Assuming a time interval of 1, the following expression is obtained:
d o f f s e t = E r . U s t a + b . ( U s t ) x
U s t = P t P t 1
Based on Equation (60), with a = 2.25 ; b = 36.2 ; x = 2.23 , it is evident that the denominator is significantly larger than the numerator in the overall fraction. As the iteration progresses, the fish’s speed gradually decreases, leading to a potential algorithm stagnation.
From (58), let α = 1 , we get
G L   D α f t + 1 = f t + 1 f t
Equation (61) denotes the derivative of order 1 in the discrete case. For Γ (x), we get
Γ α + 1 = α Γ α , Γ α + 1 = α ( α 1 ) Γ α 1 , Γ α + 1 = α ( α 1 ) ( α 2 ) Γ α 2 , Γ α + 1 = α ( α 1 ) ( α 2 ) ( α 3 ) Γ α 3
In order to generalize, let T = 1; r = 4 [72,73,74], Equation (62) can be attained.
G L   D α f t = 1 T α k = 0 r 1 k Γ ( α + 1 ) f t + 1 k T Γ ( k + 1 ) Γ ( α k + 1 ) = f t + 1 α f t 1 0.5 α 1 α f t 1 1 6 α 1 α 2 α f t 2 1 24 α 1 α 2 α 3 α f t 3
In the proposed algorithm, with a fish population size of ps and dimension Dim, the position matrix P is defined as P(ps; Dim) = [p1; p2; · · ·; pps]T, where pi = [pi;1; pi;2; · · ·; pi;Dim] represents the position of particle i. Similarly, Ppre = [Ppre1; Ppre2; Ppre3; Ppre4] denotes the past positions of the particles utilized for fractional-order speed calculation. Precisely, p 1 p r e h ( p s ;   D i m ) = p 1 p r e h , p 2 p r e h , . , p p s p r e h   T , p i p r e h = p i , 1 p r e h , p i , 2 p r e h , . , p i , D i m p r e h   , where h = 1, 2, 3, 4; i = 1, 2,…, ps. Thus, the velocity of the particle is updated by (64), where
U s , d = P i , d α p i , d p r e 1 1 2   1 α p i , d p r e 2 1 6   1 α 2 α p i , d p r e 3 1 24   1 α 2 α 3 α p i , d p r e 4
p i , d n e w = E r ,   d . U s , d a + b . ( U s , d ) x , p i , d p r e 4 = p i , d p r e 3 = p i , d p r e 2 = p i , d p r e 1 = p i , d n e w
B.
New Positions of the Offspring
Furthermore, graylings migrate back to replicate new offspring once they reach maturity as a group. Graylings need to breed in a location that enhances their chances of survival.
Hence, the positions of the new offspring should be near the global finest particle. Equation (54) is substituted with Equation (65).
p i n e w = p g b e s t + r a n d .   ( p i o l d p g b e s t )
This section provides a detailed description of the innovative algorithm utilizing fractional-order speed and new positions. The flowchart of the FOFMO can be found in Figure 5. Algorithm 1 presents the pseudo code of the FOFMO procedure, with Ratefecundity representing the fecundity rates at specific stages, Emin as the minimum energy threshold for elimination, xg,i indicating the growth status, and xeng,i representing the energy of particle i.
Algorithm 1 FOFMO
  • The dimension Dim, population size ps, a, b, x,
  • Maximum iteration Maxitn,
  • Initialize the searching space V, itn = 1, position matrix Pitn
  • Past position matrix Ppre,itn, the energy X e n g i t n
  • Compute fitness function values f ( P i i t n ) , set grow status X g i t n = 0 +
  • While itn < MaxitnΙΙ!stopCriterion do
  • for t = 1 to ps do
  •       for d = 1 to Dim do
  •       Calculate consuming energy by Equation (50)
  •       Calculate fractional-order position P i , d i t n , update the history by (64)
  •       Calculate the fitness value of new position P i , d i t n ,
  • if  f ( P i i t n ) <   f ( P g b e s t i t n ) then
  •              f ( P g b e s t i t n ) = f ( P i i t n )
  •              P g b e s t i t n = P i i t n
  •             Increased energy x e n g , i = x e n g , i + r a n d . E m a x
  •       End if
  •             Consuming energy by x e n g , i = x e n g , i f ( P i i t n ) / i f ( P i i t n )
  •          If   x e n g , i < E m i n , then
  •              x g , i = 5, the grayling died
  •       End if
  •        If   x g , i = 0 + x g , i = 1 + then
  •              x g , i = x g , i + 1 +
  •       Else if  x g , i = 2 + x g , i = 3 + then
  •       If rand()<Ratefecundity then
  •             Immigrate and produce offspring by Equation (65)
  •              x g , i = 0
  •       Else
  •              x g , i = x g , i + 1 +
  •       End if
  •       Else if  x g , i = 4 + then
  •                 Immigrate and produce offspring by Equation (65)
  •                  x g , i = 0
  •       End if
  • End for
  •                 Itn = itn + 1
  • End while
  • Output:
  • The global optimal solution: p g b e s t ,   f ( p g b e s t )

3.3.2. Coati Optimization Algorithm (COA) and Non-Dominated Sorting Genetic Algorithm (NSGA-II) Techniques

The COA and NSGA-II are optimization algorithms employed within the study to deal with the EED problem [52,53].
The COA is an unconventional optimization approach that specializes in realistic and non-convex issues [52,75]. It offers an alternative approach to dealing with complex optimization-demanding situations. In the context of the EED problem, COA is utilized as one of the optimization algorithms to acquire Pareto premiere solutions for fuel costs and emissions. It contributes to finding the superior trade-off between monetary and environmental concerns.
On the other hand, NSGA-II is a hybrid optimization technique that combines the strengths of more than one nonconventional algorithm to resolve complex troubles [53]. It is a well-known set of rules for multi-goal optimization. In the study, NSGA-II is employed as one of the optimization algorithms to attain Pareto optimal solutions for the EED issue. It helps to concurrently optimize a couple of goals, which include minimizing toxin emissions and fuel costs, even considering system constraints.
Both COA and NSGA-II play significant roles in the proposed model for optimizing economic dispatch with RESs and NGs. They contribute to locating efficient alternatives that address the mixing of thermal, renewable schemes, and natural gas units while considering environmental and financial objectives.

4. Simulation Results and Comparative Analysis

Three cases are evaluated in order to minimize both emissions and fuel costs. In general, each case includes three RESs: wind and two PVs at buses 5, 11, and 13. This can be described as follows:
  • Case I: involves three TUs and three RESs [55].
  • Case II: Changing the fuel of TU at bus 1 to NG.
  • Case III: Changing the fuel of TUs at buses 2 and 8 to NGs.

4.1. First Case

In this reported literature case, a comparison was made between the Pareto fronts (PFs) generated by the SMODE and MOEA/D techniques, as documented in [55] and presented in Table 4. The study reveals that SMODE exhibits superior diversity compared to MOEA/D, particularly regarding the cost objective. Notably, SMODE achieves fewer emission levels, specifically 0.4722 t/h, while MOEA/D attains a slightly lower fuel cost of 919.041 USD/MWh.
The variables attained utilizing the proposed FOFMO method are compared with four other considered optimizers, as shown in Table 4. Notably, the proposed FOFMO achieves a minimum cost of 911.854 USD/MWh, a minimum power loss of 5.1254 MW, and a minimum VD of 0.412 pu. While the COA technique achieves a minimum emissions level of 0.4248 t/h.
Furthermore, the FOFMO demonstrates fast solution convergence, with a computational time of approximately 120.798 s, outperforming the other algorithms, followed by the COA algorithm with a computational time of around 133.063 s. (See Figure 6).

4.2. Second Case

In this case, it is replaced the largest TU at bus 1 to NG incorporating stochastic RESs. Figure 7 exhibits the PFs of FOFMO, COA, and NSGA-II methods. Additionally, Table 5 shows the statistical simulation findings for solutions from three optimization strategies used in this event using stochastic RESs.
The variables obtained using the proposed FOFMO method are compared with two other considered optimizers, as illustrated in Table 5. Notably, the proposed FOFMO achieves a minimum cost of 798.182 USD/MWh, a minimum power loss of 4.501 MW, and minimum VD of 0.567 pu. Furthermore, the FOFMO demonstrates fast solution convergence, with a computational time of approximately 128.248 s, outperforming the other algorithms, followed by the COA algorithm with a computational time of around 145.735 s.
In terms of emission minimization, COA obtains lower emission levels of 0.4523 t/h, whilst FOFMO reaches emission levels of 0.5222 t/h. As a result, COA has a higher variety of PF than FOFMO in the direction of the cost function. In comparison to overall fuel expenses, FOFMO obtains a lower fuel cost value of 798.182 USD/MWh, whereas COA achieves a fuel cost of 804.584 USD/MWh.

4.3. Third Case

In this case, the largest thermal unit at bus 1 is retained, while the TUs at buses 2 and 8 are replaced with NGs incorporating stochastic RESs. Figure 8 showcases the PFs achieved through the utilization of the FOFMO, COA, and NSGA-II techniques. Furthermore, Table 6 presents the statistical simulation findings of solutions obtained from these three optimization methods employed in this case, considering the incorporation of stochastic RESs.
The variables obtained using the proposed FOFMO method are compared with two other considered optimizers, as illustrated in Table 6. Notably, the proposed FOFMO achieves a minimum cost of 796.478 USD/MWh, a minimum power loss of 3.214 MW, and a minimum VD of 0.509 pu. Furthermore, the FOFMO demonstrates fast solution convergence, with a computational time of approximately 136.116 s, outperforming the other algorithms, followed by the COA algorithm with a computational time of around 151.245 s. In terms of convergence, FOFMO slightly outperforms COA and NSGA-II, as shown in Figure 8.
Similarly, when comparing the minimization of emission levels, NSGA-II attains lower emission levels at a value of 0.4181 t/h, while FOFMO and COA obtain emission levels valued at 0.5591 t/h and 0.4232 t/h, respectively. Hence, NSGA-II reveals better diversity in the PFs compared to FOFMO in terms of the cost function. In contrast, regarding overall fuel costs, FOFMO attains a lower fuel cost valued at 796.478 USD/MWh, while COA and NSGA-II obtain a fuel cost of 807.782 USD/MWh and 816.451 USD/MWh, respectively.

4.4. Comparative Study

A comparative study is presented in two cases, Cases II and III, where TUs are replaced with NGs. These cases offer a framework for comparison, encompassing two perspectives: comparing the algorithms within each case and comparing the cases themselves.
The results indicate that the FOFMO technique exhibits high emission levels, making it environmentally unfavorable. However, its cost objective may encourage investment. In contrast, COA and NSGA-II achieve the lowest emission levels but have a higher cost objective in Case II and Case III, respectively.
The complexity of these findings arises from the multidisciplinary nature of the objectives. To resolve this complexity, a holistic approach is suggested, where the two cases are compared to determine which one aligns with governmental regulations and provides benefits to the community.
The exclusion of Case I from the analysis is justified due to its significantly higher emission levels and fuel costs compared to the proposed cases (Cases II and III), as illustrated in Table 7.
Upon examining the outcomes of Case I, it becomes apparent that FOFMO yields the highest emission levels of 0.5473 t/h. Conversely, in Case II, COA achieves the lowest emission levels of 0.4523 t/h, while FOFMO attains the lowest fuel costs of 798.182 USD/MWh. In Case III, COA achieves the lowest emission levels of 0.4232 t/h, while FOFMO attains the lowest fuel costs of 796.478 USD/MWh.
The analysis of the three cases reveals that Case III stands out for achieving the lowest fuel costs and emission levels. This outcome is attributed to the combination of the COA technique, which minimizes emissions, and the FOFMO technique, which reduces fuel costs. As a result, Case III successfully meets the system requirements by striking a balance between environmental and economic considerations.
The question remains as to which optimization technique in Case III provides the best compromise solution. Armed with the data provided, the power system operator can make an informed decision regarding the optimal optimization technique.
In summary, Case III offers the most favorable solution by retaining the largest thermal unit at bus 1 and replacing thermal units at buses 2 and 8 with NG units, incorporating stochastic renewable energy sources. This approach not only minimizes emissions but also reduces fuel costs, making it the most effective solution among the three cases.
Figure 9 and Figure 10 illustrate the load bus voltage profiles for two proposed cases, integrating RESs and NGs, with the worst VD values using three optimization methods. In Case II, the worst VD values obtained using FOFMO and COA are 0.56 pu and 0.79 pu, respectively. Similarly, in Case III, the worst VD values are 0.5 pu and 0.64 pu, respectively, using FOFMO and COA. The variation in VD values is attributed to the higher diversity of COA, which results in smaller emission levels or larger cost objectives.

5. Conclusions

This study presents an economic-environmental dispatch model that aims to determine the optimal solutions for an integrated IEEE 30-bus system comprising conventional thermal units, natural gas units, and renewable energy sources such as wind and two photovoltaic systems. The model takes into account both the emission and total cost functions to evaluate and identify the optimal operating points that simultaneously minimize emissions and total costs.
To achieve the simultaneous minimization of emissions and fuel or generation costs, three distinct scenarios were tested. Each scenario consisted of three RESs, specifically wind and two PV units located at buses 5, 11, and 13, respectively. The first scenario involved three TUs and three RESs. In Case II, the fuel of the TU at bus 1 was replaced with NG. Case III replaced the fuels of TUs at buses 2 and 8 with NGs. Through the obtained results, it was determined that the third scenario provided the optimal compromise solution, achieving the lowest emissions and total costs.
To obtain Pareto optimal solutions simultaneously, three optimization techniques were employed: FOFMO, COA, and NSGA-II. All system constraints, including equality, inequality, and natural gas limits, were successfully met. A comparative analysis was conducted between the recent optimization methods to determine the optimal values for the EED pollutant emissions and fuel costs. The results demonstrated that FOFMO outperformed COA and NSGA-II.
Further investigation can be conducted on the EED formulation by exploring alternative optimization techniques such as reinforcement learning (RL), and others. It would be valuable to combine these techniques with suitable constraint-handling methods. Additionally, addressing the dynamic EED issue, which considers variations in load demands over time, generator ramping rates, and uncertainties associated with Renewable Energy Sources (RESs) and system limitations, presents an important area for future research.

Funding

This study is supported via funding from the Prince Sattam bin Abdulaziz University (project number PSAU/2024/R/1445).

Data Availability Statement

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

Acknowledgments

The author acknowledges the Prince Sattam bin Abdulaziz University for their technical and financial support.

Conflicts of Interest

The author declares no conflicts of interest.

Abbreviations

ANNArtificial Neural Network
COACoati Optimization Algorithm
DEDifferential Evolution
EDNSExpected Demand Not Supplied
EEDEnvironmental Economic Dispatch
ESEnergy Storage
FCFractional Calculus
FDsFractional Derivatives
FOFMOFractional-Order Fish Migration Optimization
GAGenetic Algorithm
HHVHigh Heat Value of Natural Gas
LMLagrange’s Method
LOLPLoss Of Load Probability
NGNatural Gas
NSGA-IINon-Dominated Sorting Genetic Algorithm-II
PDFs Probability Density Functions
POZs Prohibited Operating Zones
PSOParticle Swarm Optimization
PVPhotovoltaic
RESRenewable Energy Sources
SCOSine Cosine Optimization
SFLAShuffle Frog-Jump Approach
SOAShark Optimization Algorithm
TUsThermal Units
VDVoltage Deviation
WTWind Turbine
Nomenclature
α Scale factor of the wind turbine
β Shape factor of the wind turbine
C d P V Direct cost of the photovoltaic system
C d W P Direct cost of the wind turbine
C N G Constant-coefficient of natural gas
C r P V Reserve capacity cost of the photovoltaic system
C r W P Reserve capacity cost of the wind turbine
C S P V Storage units cost of the photovoltaic system
C S W P Storage units cost of the wind turbine
C t o t Total cost of the fuel or generation
C t o t P V Total cost of the photovoltaic generation unit
C t o t W P Total cost of the wind turbine generation unit
C t o t P N G Total cost of the natural gas unit
C t o t P T U Total cost of the thermal unit
C t o t P R E S Total cost of renewable energy sources
δ i j Phase difference between the buses i and j
d P l i n e Internal diameter of the pipe in millimeters
η N G Efficiency of the gas turbine
E t o t Total emission
f v v Probability of wind speed
f Friction factor
g N G i Initial and operation costs coefficient for ith natural gas units
GSolar irradiance
G s t d Standard solar irradiance
G q ( i j ) Transconductance of branch q connected to bus i and bus j
K d W P Direct cost parameter of the wind turbine
K r W P Reserve capacity cost parameter of the wind turbine
K S W P Storage unit cost parameter of the wind turbine
L P l i n e Length of pipeline in meters
N G Number of generator buses
N L Number of load buses
n l Number of branches in the network
P x Absolute upstream (inlet) pressure
P y Absolute downstream (outlet) pressure
P b Absolute pressure
P l o s s Network power loss
P p Pipeline pressure
P P V a c t Actual power of the photovoltaic system
P P V r Rated power of the photovoltaic system
P P V s c h Scheduled power of the photovoltaic system
P T U i m i n Minimum power of the ith thermal power generator unit
P W a c t Actual power of the wind turbine
P W P r Rated power of the wind turbine
P W P s c h Scheduled power of the wind turbine
S L i Ith branches’ capacity limit
S N G Specific gravity of natural gas
T N G Average temperature of the flowing gas in kelvin
T s Standard temperature in kelvin
U N G Flow velocity of the natural gas in m/sec
v Wind speed
V G i Voltage of the ith on generator bus
v i n Cut-in speed of the wind turbine
V L i Voltage of the ith load bus
V N G Volume on the remind loads of the natural gas
v o u t Cut-out speed of the wind turbine
v r Rated speed of the wind turbine
Z Average compressibility factor of natural gas

Appendix A

Table A1. Cost factors (direct, reserve, and standby) for stochastic RESs.
Table A1. Cost factors (direct, reserve, and standby) for stochastic RESs.
Wind (Bus 5)Solar (Bus 11 and Bus 13)
Direct cost factors (USD/MW) K d W T     = 1.7 K d P V   = 1.6
Reserve cost factors (USD/MW) K r W T   = 3 K r P V   = 3
Penalty cost factors (USD/MW) K s W T   = 1.4 K s P V   = 1.4

References

  1. Khalid, M. Smart Grids and Renewable Energy Systems: Perspectives and Grid Integration Challenges. Energy Strategy Rev. 2024, 51, 101299. [Google Scholar] [CrossRef]
  2. Trojovský, P.; Trojovská, E.; Akbari, E. Economical-Environmental-Technical Optimal Power Flow Solutions Using a Novel Self-Adaptive Wild Geese Algorithm with Stochastic Wind and Solar Power. Sci. Rep. 2024, 14, 4135. [Google Scholar] [CrossRef] [PubMed]
  3. Kabeyi, M.J.B.; Olanrewaju, O.A. Environmental Impacts of Power Plants and Energy Conversion Systems. In Green Energy and Technology; Chen, L., Ed.; Springer Nature: Cham, Switzerland, 2024; Volume Part F2329, pp. 445–469. ISBN 978-3-031-49787-2. [Google Scholar]
  4. Singh, M.K.; Malek, J.; Sharma, H.K.; Kumar, R. Converting the Threats of Fossil Fuel-Based Energy Generation into Opportunities for Renewable Energy Development in India. Renew. Energy 2024, 224, 120153. [Google Scholar] [CrossRef]
  5. Marzbani, F.; Abdelfatah, A. Economic Dispatch Optimization Strategies and Problem Formulation: A Comprehensive Review. Energies 2024, 17, 550. [Google Scholar] [CrossRef]
  6. Mundotiya, P.; Tiwari, H.P.; Alotaibi, M.A.; Singh, S.; Gupta, N.K.; Malik, H.; Marquez, F.P.G. A Novel Hybrid Approach for Hydrothermal Scheduling Using Mathematical and Metaheuristic Search Methods. IEEE Access 2024, 12, 22420–22432. [Google Scholar] [CrossRef]
  7. Manoharan, P.; Chandrasekaran, K.; Chandran, R.; Ravichandran, S.; Mohammad, S.; Jangir, P. An Effective Strategy for Unit Commitment of Microgrid Power Systems Integrated with Renewable Energy Sources Including Effects of Battery Degradation and Uncertainties. Environ. Sci. Pollut. Res. 2024, 31, 11037–11080. [Google Scholar] [CrossRef] [PubMed]
  8. Keshuov, S.; Moldybayeva, N. Selection of Optimal Structure of an Energy-Supply System for Objects in an Agro-Industrial Complex Based on Renewable-Energy Sources. Renew. Energy Focus 2024, 48, 100533. [Google Scholar] [CrossRef]
  9. Ismaeel, A.A.K.; Houssein, E.H.; Khafaga, D.S.; Abdullah Aldakheel, E.; AbdElrazek, A.S.; Said, M. Performance of Osprey Optimization Algorithm for Solving Economic Load Dispatch Problem. Mathematics 2023, 11, 4107. [Google Scholar] [CrossRef]
  10. Pan, S.; Jian, J.; Chen, H.; Yang, L. A Full Mixed-Integer Linear Programming Formulation for Economic Dispatch with Valve-Point Effects, Transmission Loss and Prohibited Operating Zones. Electr. Power Syst. Res. 2020, 180, 106061. [Google Scholar] [CrossRef]
  11. Chang, X.; Xu, Y.; Sun, H.; Khan, I. A Distributed Robust Optimization Approach for the Economic Dispatch of Flexible Resources. Int. J. Electr. Power Energy Syst. 2021, 124, 106360. [Google Scholar] [CrossRef]
  12. Sadeghian, O.; Moradzadeh, A.; Mohammadi-Ivatloo, B.; Abapour, M.; Marquez, F.P.G. Generation Units Maintenance in Combined Heat and Power Integrated Systems Using the Mixed Integer Quadratic Programming Approach. Energies 2020, 13, 2840. [Google Scholar] [CrossRef]
  13. Hassan, M.H.; Kamel, S.; Shaikh, M.S.; Alquthami, T.; Hussien, A.G. Supply-Demand Optimizer for Economic Emission Dispatch Incorporating Price Penalty Factor and Variable Load Demand Levels. IET Gener. Transm. Distrib. 2023, 17, 3211–3231. [Google Scholar] [CrossRef]
  14. Ahmad, A.; Kashif, S.A.R.; Ashraf, A.; Gulzar, M.M.; Alqahtani, M.; Khalid, M. Coordinated Economic Operation of Hydrothermal Units with HVDC Link Based on Lagrange Multipliers. Mathematics 2023, 11, 1610. [Google Scholar] [CrossRef]
  15. 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]
  16. Yang, Q.; Chen, G.; Wang, T. ADMM-Based Distributed Algorithm for Economic Dispatch in Power Systems with Both Packet Drops and Communication Delays. IEEE/CAA J. Autom. Sin. 2020, 7, 842–852. [Google Scholar] [CrossRef]
  17. Alawode, K.O.; Jubril, A.M.; Kehinde, L.O.; Ogunbona, P.O. Semidefinite Programming Solution of Economic Dispatch Problem with Non-Smooth, Non-Convex Cost Functions. Electr. Power Syst. Res. 2018, 164, 178–187. [Google Scholar] [CrossRef]
  18. Tabassum, M.F.; Akram, S.; Hassan, S.; Karim, R.; Naik, P.A.; Farman, M.; Yavuz, M.; Naik, M.U.D.; Ahmad, H. Differential Gradient Evolution plus Algorithm for Constraint Optimization Problems: A Hybrid Approach. Int. J. Optim. Control Theor. Appl. 2021, 11, 158–177. [Google Scholar] [CrossRef]
  19. Liu, J.; Zhao, F.; Li, Y.; Zhou, H. A New Global Sine Cosine Algorithm for Solving Economic Emission Dispatch Problem. Inf. Sci. 2023, 648, 119569. [Google Scholar] [CrossRef]
  20. Said, M.; Houssein, E.H.; Deb, S.; Ghoniem, R.M.; Elsayed, A.G. Economic Load Dispatch Problem Based on Search and Rescue Optimization Algorithm. IEEE Access 2022, 10, 47109–47123. [Google Scholar] [CrossRef]
  21. Sakthivel, V.P.; Suman, M.; Sathya, P.D. Combined Economic and Emission Power Dispatch Problems through Multi-Objective Squirrel Search Algorithm. Appl. Soft Comput. 2021, 100, 106950. [Google Scholar] [CrossRef]
  22. Rateele, T.N.; Thamae, L.Z. An Optimization Approach for the Economic Dispatch Incorporating Renewable Energy Resources into Lesotho Power Sources Portfolio. Heliyon 2023, 9, e14748. [Google Scholar] [CrossRef] [PubMed]
  23. Basu, M. Combined Heat and Power Economic Emission Dispatch Using Nondominated Sorting Genetic Algorithm-II. Int. J. Electr. Power Energy Syst. 2013, 53, 135–141. [Google Scholar] [CrossRef]
  24. Ajayi, O.; Heymann, R. Day-Ahead Combined Economic and Emission Dispatch with Spinning Reserve Consideration Using Moth Swarm Algorithm for a Data Centre Load. Heliyon 2021, 7, e08054. [Google Scholar] [CrossRef] [PubMed]
  25. Hussain, S.; Al-Hitmi, M.; Khaliq, S.; Hussain, A.; Saqib, M.A. Implementation and Comparison of Particle Swarm Optimization and Genetic Algorithm Techniques in Combined Economic Emission Dispatch of an Independent Power Plant. Energies 2019, 12, 2037. [Google Scholar] [CrossRef]
  26. Xu, X.; Hu, Z.; Su, Q.; Xiong, Z.; Liu, M. Multi-Objective Learning Backtracking Search Algorithm for Economic Emission Dispatch Problem. Soft Comput. 2021, 25, 2433–2452. [Google Scholar] [CrossRef]
  27. Lai, W.; Zheng, X.; Song, Q.; Hu, F.; Tao, Q.; Chen, H. Multi-Objective Membrane Search Algorithm: A New Solution for Economic Emission Dispatch. Appl. Energy 2022, 326, 119969. [Google Scholar] [CrossRef]
  28. Elattar, E.E. Environmental Economic Dispatch with Heat Optimization in the Presence of Renewable Energy Based on Modified Shuffle Frog Leaping Algorithm. Energy 2019, 171, 256–269. [Google Scholar] [CrossRef]
  29. Ali, E.S.; Abd Elazim, S.M.; Hakmi, S.H.; Mosaad, M.I. Optimal Allocation and Size of Renewable Energy Sources as Distributed Generations Using Shark Optimization Algorithm in Radial Distribution Systems. Energies 2023, 16, 3983. [Google Scholar] [CrossRef]
  30. Sutar, M.; Jadhav, H.T. A Modified Artificial Bee Colony Algorithm Based on a Non-Dominated Sorting Genetic Approach for Combined Economic-Emission Load Dispatch Problem. Appl. Soft Comput. 2023, 144, 110433. [Google Scholar] [CrossRef]
  31. Xin-gang, Z.; Ji, L.; Jin, M.; Ying, Z. An Improved Quantum Particle Swarm Optimization Algorithm for Environmental Economic Dispatch. Expert Syst. Appl. 2020, 152, 113370. [Google Scholar] [CrossRef]
  32. Mishra, S.K.; Gupta, V.K.; Kumar, R.; Swain, S.K.; Mohanta, D.K. Multi-Objective Optimization of Economic Emission Load Dispatch Incorporating Load Forecasting and Solar Photovoltaic Sources for Carbon Neutrality. Electr. Power Syst. Res. 2023, 223, 109700. [Google Scholar] [CrossRef]
  33. Zdiri, M.A.; Ben Salem, F.; Abdallah, H.H.; Alaas, Z.; Dhouib, B. Load Flow Analysis and the Impact of a Solar PV Generator in a Radial Distribution Network. Eng. Technol. Appl. Sci. Res. 2023, 13, 10078–10085. [Google Scholar] [CrossRef]
  34. Abdelaziz, A.Y.; Ali, E.S.; Abd Elazim, S.M. Combined Economic and Emission Dispatch Solution Using Flower Pollination Algorithm. Int. J. Electr. Power Energy Syst. 2016, 80, 264–274. [Google Scholar] [CrossRef]
  35. Ramalingam, R.; Karunanidy, D.; Alshamrani, S.S.; Rashid, M.; Mathumohan, S.; Dumka, A. Oppositional Pigeon-Inspired Optimizer for Solving the Non-Convex Economic Load Dispatch Problem in Power Systems. Mathematics 2022, 10, 3315. [Google Scholar] [CrossRef]
  36. Xiong, G.; Shuai, M.; Hu, X. Combined Heat and Power Economic Emission Dispatch Using Improved Bare-Bone Multi-Objective Particle Swarm Optimization. Energy 2022, 244, 123108. [Google Scholar] [CrossRef]
  37. Kheshti, M.; Kang, X.; Li, J.; Regulski, P.; Terzija, V. Lightning Flash Algorithm for Solving Non-Convex Combined Emission Economic Dispatch with Generator Constraints. IET Gener. Transm. Distrib. 2018, 12, 104–116. [Google Scholar] [CrossRef]
  38. Said, M.; El-Rifaie, A.M.; Tolba, M.A.; Houssein, E.H.; Deb, S. An Efficient Chameleon Swarm Algorithm for Economic Load Dispatch Problem. Mathematics 2021, 9, 2770. [Google Scholar] [CrossRef]
  39. Zhang, L.; Khishe, M.; Mohammadi, M.; Mohammed, A.H. Environmental Economic Dispatch Optimization Using Niching Penalized Chimp Algorithm. Energy 2022, 261, 125259. [Google Scholar] [CrossRef]
  40. Sundaram, A. Multiobjective Multi-Verse Optimization Algorithm to Solve Combined Economic, Heat and Power Emission Dispatch Problems. Appl. Soft Comput. J. 2020, 91, 106195. [Google Scholar] [CrossRef]
  41. Liang, H.; Liu, Y.; Li, F.; Shen, Y. A Multiobjective Hybrid Bat Algorithm for Combined Economic/Emission Dispatch. Int. J. Electr. Power Energy Syst. 2018, 101, 103–115. [Google Scholar] [CrossRef]
  42. Sarhan, S.; Shaheen, A.; El-Sehiemy, R.; Gafar, M. A Multi-Objective Teaching–Learning Studying-Based Algorithm for Large-Scale Dispatching of Combined Electrical Power and Heat Energies. Mathematics 2022, 10, 2278. [Google Scholar] [CrossRef]
  43. Basak, S.; Bhattacharyya, B.; Dey, B. Combined Economic Emission Dispatch on Dynamic Systems Using Hybrid CSA-JAYA Algorithm. Int. J. Syst. Assur. Eng. Manag. 2022, 13, 2269–2290. [Google Scholar] [CrossRef]
  44. Ziane, I.; Benhamida, F.; Graa, A. Simulated Annealing Algorithm for Combined Economic and Emission Power Dispatch Using Max/Max Price Penalty Factor. Neural Comput. Appl. 2017, 28, 197–205. [Google Scholar] [CrossRef]
  45. Shaheen, A.M.; Elsayed, A.M.; Ginidi, A.R.; EL-Sehiemy, R.A.; Alharthi, M.M.; Ghoneim, S.S.M. A Novel Improved Marine Predators Algorithm for Combined Heat and Power Economic Dispatch Problem. Alex. Eng. J. 2022, 61, 1834–1851. [Google Scholar] [CrossRef]
  46. Mahdi, F.P.; Vasant, P.; Abdullah-Al-Wadud, M.; Kallimani, V.; Watada, J. Quantum-Behaved Bat Algorithm for Many-Objective Combined Economic Emission Dispatch Problem Using Cubic Criterion Function. Neural Comput. Appl. 2019, 31, 5857–5869. [Google Scholar] [CrossRef]
  47. Bhargava, G.; Yadav, N.K. Solving Combined Economic Emission Dispatch Model via Hybrid Differential Evaluation and Crow Search Algorithm. Evol. Intell. 2022, 15, 1161–1169. [Google Scholar] [CrossRef]
  48. Farhan Tabassum, M.; Saeed, M.; Ahmad Chaudhry, N.; Ali, J.; Farman, M.; Akram, S. Evolutionary Simplex Adaptive Hooke-Jeeves Algorithm for Economic Load Dispatch Problem Considering Valve Point Loading Effects. Ain Shams Eng. J. 2021, 12, 1001–1015. [Google Scholar] [CrossRef]
  49. Moustafa, G.; Alnami, H.; Hakmi, S.H.; Shaheen, A.M.; Ginidi, A.R.; Elshahed, M.A.; Mansour, H.S.E. A Novel Mantis Search Algorithm for Economic Dispatch in Combined Heat and Power Systems. IEEE Access 2024, 12, 2674–2689. [Google Scholar] [CrossRef]
  50. Ginidi, A.; Elsayed, A.; Shaheen, A.; Elattar, E.; El-Sehiemy, R. An Innovative Hybrid Heap-Based and Jellyfish Search Algorithm for Combined Heat and Power Economic Dispatch in Electrical Grids. Mathematics 2021, 9, 2053. [Google Scholar] [CrossRef]
  51. Guo, B.; Zhuang, Z.; Pan, J.S.; Chu, S.C. Optimal Design and Simulation for PID Controller Using Fractional-Order Fish Migration Optimization Algorithm. IEEE Access 2021, 9, 8808–8819. [Google Scholar] [CrossRef]
  52. Dehghani, M.; Montazeri, Z.; Trojovská, E.; Trojovský, P. Coati Optimization Algorithm: A New Bio-Inspired Metaheuristic Algorithm for Solving Optimization Problems. Knowl. Based Syst. 2023, 259, 110011. [Google Scholar] [CrossRef]
  53. Araujo, E.X.S.; Cerbantes, M.C.; Mantovani, J.R.S. Optimal Power Flow with Renewable Generation: A Modified NSGA-II-Based Probabilistic Solution Approach. J. Control Autom. Electr. Syst. 2020, 31, 979–989. [Google Scholar] [CrossRef]
  54. Ilyas, M.A.; Abbas, G.; Alquthami, T.; Awais, M.; Rasheed, M.B. Multi-Objective Optimal Power Flow with Integration of Renewable Energy Sources Using Fuzzy Membership Function. IEEE Access 2020, 8, 143185–143200. [Google Scholar] [CrossRef]
  55. Omar, A.I.; Ali, Z.M.; Al-Gabalawy, M.; Abdel Aleem, S.H.E.; Al-Dhaifallah, M. Multi-Objective Environmental Economic Dispatch of an Electricity System Considering Integrated Natural Gas Units and Variable Renewable Energy Sources. Mathematics 2020, 8, 1100. [Google Scholar] [CrossRef]
  56. Yanto, D.T.P.; Akhmadeev, R.; Hamad, H.S.; Alawadi, A.H.R.; Abdullayev, A.B.; Romero-Parra, R.M.; Fooladi, H. Development and Investigation of a Pollutants Emission Reduction Process from a Coal-Gasifcation Power Plant Integrated with Fuel Cell and Solar Energy. Int. J. Low-Carbon Technol. 2023, 18, 1120–1133. [Google Scholar] [CrossRef]
  57. Li, J.; Li, X.; Yan, P.; Zhou, G.; Liu, J.; Yu, D. Thermodynamics, Flexibility and Techno-Economics Assessment of a Novel Integration of Coal-Fired Combined Heating and Power Generation Unit and Compressed Air Energy Storage. Appl. Energy 2023, 339, 120924. [Google Scholar] [CrossRef]
  58. Liu, J.; Song, D.; Li, Q.; Yang, J.; Hu, Y.; Fang, F.; Hoon Joo, Y. Life Cycle Cost Modelling and Economic Analysis of Wind Power: A State of Art Review. Energy Convers. Manag. 2023, 277, 116628. [Google Scholar] [CrossRef]
  59. Aldin, N.A.N.; Abdellatif, W.S.E.; Salem Elbarbary, Z.M.; Omar, A.I.; Mahmoud, M.M. Robust Speed Controller for PMSG Wind System Based on Harris Hawks Optimization via Wind Speed Estimation: A Real Case Study. IEEE Access 2023, 11, 5929–5943. [Google Scholar] [CrossRef]
  60. Mousa, A.G.; Abdel Aleem, S.H.E.; Ibrahim, A.M. Mathematical analysis of maximum power points and currents based maximum power point tracking in solar photovoltaic system: A solar powered water pump application. Int. Rev. Electr. Eng 2016, 11, 97. [Google Scholar] [CrossRef]
  61. Omar, A.I.; Mohsen, M.; Abd-Allah, M.A.; Salem Elbarbary, Z.M.; Said, A. Induced Overvoltage Caused by Indirect Lightning Strikes in Large Photovoltaic Power Plants and Effective Attenuation Techniques. IEEE Access 2022, 10, 112934–112947. [Google Scholar] [CrossRef]
  62. Ali, Z.M.; Aleem, S.H.E.A.; Omar, A.I.; Mahmoud, B.S. Economical-Environmental-Technical Operation of Power Networks with High Penetration of Renewable Energy Systems Using Multi-Objective Coronavirus Herd Immunity Algorithm. Mathematics 2022, 10, 1201. [Google Scholar] [CrossRef]
  63. Zhou, B.; Ai, X.; Fang, J.; Li, K.; Yao, W.; Chen, Z.; Wen, J. Function-Space Optimization to Coordinate Multi-Energy Storage across the Integrated Electricity and Natural Gas System. Int. J. Electr. Power Energy Syst. 2023, 151, 109181. [Google Scholar] [CrossRef]
  64. Khatibi, M.; Rabiee, A.; Bagheri, A. Integrated Electricity and Gas Systems Planning: New Opportunities, and a Detailed Assessment of Relevant Issues. Sustainability 2023, 15, 6602. [Google Scholar] [CrossRef]
  65. Algabalawy, M.A.; Abdelaziz, A.Y.; Mekhamer, S.F.; Abdel Aleem, S.H.E. Considerations on Optimal Design of Hybrid Power Generation Systems Using Whale and Sine Cosine Optimization Algorithms. J. Electr. Syst. Inf. Technol. 2018, 5, 312–325. [Google Scholar] [CrossRef]
  66. Li, G.; Zhang, R.; Jiang, T.; Chen, H.; Bai, L.; Li, X. Security-Constrained Bi-Level Economic Dispatch Model for Integrated Natural Gas and Electricity Systems Considering Wind Power and Power-to-Gas Process. Appl. Energy 2017, 194, 696–704. [Google Scholar] [CrossRef]
  67. Avalos, R.; Fitzgerald, T.; Rucker, R.R. Measuring the Effects of Natural Gas Pipeline Constraints on Regional Pricing and Market Integration. Energy Econ. 2016, 60, 217–231. [Google Scholar] [CrossRef]
  68. Pan, J.S.; Tsai, P.W.; Liao, Y. Bin Fish Migration Optimization Based on the Fishy Biology. In Proceedings of the 2010 4th International Conference on Genetic and Evolutionary Computing, ICGEC 2010, Shenzhen, China, 13–15 December 2010; pp. 783–786. [Google Scholar]
  69. Zhou, J.; Ali, M.A.; Sharma, K.; Alizadeh, A.; Almojil, S.F.; Almohana, A.I.; Fahmi Alali, A.; Almoalimi, K.T.; Ibrahim, B.F. Improved Fish Migration Optimization Method to Identify PEMFC Parameters. Int. J. Hydrogen Energy 2023, 48, 20028–20040. [Google Scholar] [CrossRef]
  70. Huang, G.; Qin, H.; Chen, Q.; Shi, Z.; Jiang, S.; Huang, C. Research on Application of Fractional Calculus Operator in Image Underlying Processing. Fractal Fract. 2024, 8, 37. [Google Scholar] [CrossRef]
  71. Vieira, L.C.; Costa, R.S.; Valério, D. An Overview of Mathematical Modelling in Cancer Research: Fractional Calculus as Modelling Tool. Fractal Fract. 2023, 7, 595. [Google Scholar] [CrossRef]
  72. Mousavi, Y.; Alfi, A. Fractional Calculus-Based Firefly Algorithm Applied to Parameter Estimation of Chaotic Systems. Chaos Solitons Fractals 2018, 114, 202–215. [Google Scholar] [CrossRef]
  73. Mousavi, Y.; Alfi, A.; Kucukdemiral, I.B. Enhanced Fractional Chaotic Whale Optimization Algorithm for Parameter Identification of Isolated Wind-Diesel Power Systems. IEEE Access 2020, 8, 140862–140875. [Google Scholar] [CrossRef]
  74. Chou, F.I.; Huang, T.H.; Yang, P.Y.; Lin, C.H.; Lin, T.C.; Ho, W.H.; Chou, J.H. Controllability of Fractional-Order Particle Swarm Optimizer and Its Application in the Classification of Heart Disease. Appl. Sci. 2021, 11, 11517. [Google Scholar] [CrossRef]
  75. Ćalasan, M.; Aleem, S.H.A.; Zobaa, A.F. A new approach for parameters estimation of double and triple diode models of photovoltaic cells based on iterative Lambert W function. Sol. Energy 2021, 218, 392–412. [Google Scholar] [CrossRef]
Figure 1. Modified IEEE 30-bus system investigated.
Figure 1. Modified IEEE 30-bus system investigated.
Fractalfract 08 00350 g001
Figure 2. The methodology for the presented techniques.
Figure 2. The methodology for the presented techniques.
Fractalfract 08 00350 g002
Figure 3. PDFs: (a) Weibull distribution of wind velocity and (b) lognormal distribution of solar irradiation.
Figure 3. PDFs: (a) Weibull distribution of wind velocity and (b) lognormal distribution of solar irradiation.
Fractalfract 08 00350 g003aFractalfract 08 00350 g003b
Figure 4. Life cycle graph of the grayling.
Figure 4. Life cycle graph of the grayling.
Fractalfract 08 00350 g004
Figure 5. The flowchart of the presented FOFMO technique.
Figure 5. The flowchart of the presented FOFMO technique.
Fractalfract 08 00350 g005
Figure 6. Results of PFs utilizing FOFMO, COA, and NSGA-II approaches for Case I.
Figure 6. Results of PFs utilizing FOFMO, COA, and NSGA-II approaches for Case I.
Fractalfract 08 00350 g006
Figure 7. Results of PFs utilizing FOFMO, COA, and NSGA-II approaches for Case II.
Figure 7. Results of PFs utilizing FOFMO, COA, and NSGA-II approaches for Case II.
Fractalfract 08 00350 g007
Figure 8. Results of PFs utilizing FOFMO, COA, and NSGA-II techniques for Case III.
Figure 8. Results of PFs utilizing FOFMO, COA, and NSGA-II techniques for Case III.
Fractalfract 08 00350 g008
Figure 9. Load bus voltages for worst VD values utilizing the three optimizations in Case II.
Figure 9. Load bus voltages for worst VD values utilizing the three optimizations in Case II.
Fractalfract 08 00350 g009
Figure 10. Load bus voltages for worst VD values utilizing the three optimizations in Case III.
Figure 10. Load bus voltages for worst VD values utilizing the three optimizations in Case III.
Fractalfract 08 00350 g010
Table 1. The model descriptions [54].
Table 1. The model descriptions [54].
ItemQuantitySpecifications
Units63 TUs and 3 RESs
TUs/NGs3Bus 1 (slack), bus 2, and bus 8
Wind turbine (WT)25Bus 5, 75 MW
PV2Bus 11 and 13, 50 MW
Active loads-283.4 MW
Reactive loads-126.2 MVAr
Number of load buses24-
Table 2. Emission and cost factors of the TUs [57].
Table 2. Emission and cost factors of the TUs [57].
Emission Factors
GeneratorBus φ T U ψ T U ω T U τ T U ξ T U
(t/h)(t/pu. MWh)(t/pu. MW2h)(t/h)(pu. MW−1)
TU110.04091−0.055540.06490.00026.667
TU220.02543−0.060470.056380.00053.333
TU380.05326−0.03550.03380.0022
Cost Factors
GeneratorBus a T U b T U c T U d T U e T U
(USD/h)(USD/MWh)(USD/MW2h)(USD/h)(MW−1)
TU113020.00375180.037
TU22251.750.0175160.038
TU38203.250.00834120.045
Table 3. Emission and cost factors of the NGs [65].
Table 3. Emission and cost factors of the NGs [65].
Emission Factors
GeneratorBus φ N G ψ N G ω N G
(t/h)(t/pu. MWh)(t/pu. MW2h)
NG110.02091−0.075540.04490
NG220.02543−0.050470.03638
NG380.03326−0.055500.01380
Cost Factors
GeneratorBus a N G b N G c N G
(USD/h)(USD/MWh)(USD/MW2h)
NG11141.060.00175
NG22151.050.0105
NG38171.250.02434
Table 4. Complete statistical findings of the optimization approaches for Case I.
Table 4. Complete statistical findings of the optimization approaches for Case I.
VariablesUnitsMin.Max.MOEA/DSMODEFOFMOCOANSGA-II
PTU1MW50140117.117111.912122.054111.247120.045
PTU22080656542.57442.84943.490
PTU3103518.40423.55619.0217.08916.025
PWP07555.44854.05956.44857.05958.448
Ppv (Bus 11)05017.65018.43720.65019.43721.650
Ppv (Bus 13)05015.32715.75617.32718.75618.327
Ploss 5.54305.31495.12545.48615.5258
Q1MVAr−501402.1292.789−1.035−0.944−4.775
Q2−206021.42034.50521.40214.35114.801
Q5−157037.72836.17025.87035.00135.002
Q8−306027.10320.37740.00140.00140.002
Q11−203024.91223.42021.00520.190719.183
Q13−202520.32915.63022.27821.33425.002
V1pu0.961.101.0771.0771.0771.0781.079
V20.961.101.06491.0671.0701.0721.074
V50.961.101.04451.03631.0641.0651.067
V80.961.101.04031.03631.04281.0431.045
V110.961.101.08791.07791.0421.0411.045
V130.961.101.06031.04331.0271.0271.031
VD 0.45310.4220.4120.46840.4861
Total cost (USD/MWh) 919.041927.050911.854912.554914.247
Emission (t/h) 0.62220.47220.54730.42480.5487
The boldface value denotes the optimal outcome obtained through the proposed algorithms.
Table 5. Complete statistical findings of the optimization approaches for Case II.
Table 5. Complete statistical findings of the optimization approaches for Case II.
Variables Min.Max.FOFMOCOANSGA-II
PNG1MW5014072.84164.87875.541
PTU2208080.00274.22142.574
PTU3103535.00134.98210.578
PWP07533.74431.35132.744
Ppv (Bus 11)05039.73739.60539.247
Ppv (Bus 13)05028.43137.83432.481
Ploss 4.5014.8925.682
Q1MVAr−5014027.312−29.199−11.36590
Q2−206016.33267.74421.86743
Q5−157030.88164.07225.95090
Q8−306051.18328.01240.00000
Q11−20302.45211.82120.85680
Q13−20253.457−8.32122.28300
V1pu0.961.101.0911.0621.07869
V20.961.101.0851.1021.07306
V50.961.101.0551.0841.06652
V80.961.101.0541.0331.04903
V110.961.101.0111.0321.04444
V130.961.101.0120.9711.03090
VD 0.5670.7940.94580
Wgencost 102.394124.889126.9591
Sgencost (Bus 11) 147.248150.245114.9636
Sgencost (Bus 13) 97.385143.72139.9283
mass flow (m3) 63.55271.36578.345
P1_cost 9.53410.70510.975
Total cost (USD/h) 798.182804.584812.4878
Emission (t/h) 0.52220.45230.6823
Fuelvlvcost 340.821327.402353.139
The boldface value denotes the optimal outcome obtained through the proposed algorithms.
Table 6. Complete statistical findings of the optimization approach for Case III.
Table 6. Complete statistical findings of the optimization approach for Case III.
VariablesUnitsMinMaxFOFMOCOANSGA-II
PTU1MW50140121.48252.811112.811
PNG2208080.00280.00180.001
PNG3104535.00543.19841.178
Pw07569.56452.29662.246
Ppv (Bus 11)05023.93732.95942.159
Ppv (Bus 13)05018.55628.58834.588
Ploss 3.2144.0445.374
Q1MVAr−5014030.127−5.72950.127
Q2−2060−5.8676.156−9.847
Q5−154017.85736.58921.357
Q8−305048.55849.17943.518
Q11−20255.52317.4648.573
Q13−202515.68726.71221.648
V1pu0.961.101.0581.0061.01
V20.961.101.0591.0071.004
V50.961.101.0080.9911.07
V80.961.101.0131.0141.019
V110.961.101.0271.0691.078
V130.961.101.0871.0581.096
VD 0.5090.6440.601
Wgencost 99.534193.20399.534
Sgencost (Bus 11) 32.90251.60632.902
Sgencost (Bus 13) 27.89239.78127.892
mass flow_02 172.832172.837172.832
mass flow_08 75.62193.31775.621
Total cost (USD/MWh) 796.478807.782816.451
Emission (t/h) 0.55910.42320.4181
Fuelvlvcost 336.871147.932325.171
The boldface value denotes the optimal outcome obtained through the proposed algorithms.
Table 7. Comparative study between optimization methods for three cases.
Table 7. Comparative study between optimization methods for three cases.
ComparisonCase ICase IICase III
ApproachesFOFMOCOAFOFMOCOAFOFMOCOA
No. of iterations300300300300300300
No. of population100100100100100100
Computation time (s)120.798133.063128.248145.735136.116151.245
Total cost (USD/MWh)911.854912.554798.182804.584796.478807.782
Emission (t/h)0.54730.42480.52220.45230.55910.4232
The boldface value denotes the optimal outcome obtained through the proposed algorithms.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aldosary, A. Optimizing Economic Dispatch with Renewable Energy and Natural Gas Using Fractional-Order Fish Migration Algorithm. Fractal Fract. 2024, 8, 350. https://doi.org/10.3390/fractalfract8060350

AMA Style

Aldosary A. Optimizing Economic Dispatch with Renewable Energy and Natural Gas Using Fractional-Order Fish Migration Algorithm. Fractal and Fractional. 2024; 8(6):350. https://doi.org/10.3390/fractalfract8060350

Chicago/Turabian Style

Aldosary, Abdallah. 2024. "Optimizing Economic Dispatch with Renewable Energy and Natural Gas Using Fractional-Order Fish Migration Algorithm" Fractal and Fractional 8, no. 6: 350. https://doi.org/10.3390/fractalfract8060350

APA Style

Aldosary, A. (2024). Optimizing Economic Dispatch with Renewable Energy and Natural Gas Using Fractional-Order Fish Migration Algorithm. Fractal and Fractional, 8(6), 350. https://doi.org/10.3390/fractalfract8060350

Article Metrics

Back to TopTop