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:
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:
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.