Multi-objective real-time integrated solar-wind-thermal power dispatch by using meta-heuristic technique

: The elevated demand for electrical power, expeditious expenditure of fossil fuels, and degradation of the environment because of power generation have renewed attentiveness to renewable energy resources (RER). The rapid augmentation of RER increases the convolutions in leveling the demand and generation of electrical power. In this paper, an elaborated (cid:2009) -constrained simplex method (ACSM) is recommended for multi-objective power dispatch problems. This methodology is devised after synthesizing the non-linear simplex method (SM) with the (cid:2009) -constrained method (ACM) and the evolutionary method (EM). ACSM can transfigure an optimization technique for the constrained problems by reinstating standard juxtapositions with (cid:2009) -level collations. The insertion of mutations and multi-simplexes can explore the periphery of the workable zone. It can also manage the fastness of convergence and therefore, the high precision solution can be obtained. A real-time multi-objective coordinated solar-wind-thermal power scheduling problem is framed. Two conflicting objectives (operating cost and emission) are satisfied. The case studies are carried out for Muppandal (Tamil Nadu), Jaisalmer (Rajasthan), and Okha (Gujarat), India. The annual solar and wind data are analyzed by using Normal Distribution and Weibull Distribution Density Factor, respectively. The presented technique is inspected on numerous archetype functions and systems. The results depict the prevalence of ACSM over particle swarm optimization (PSO), simplex method with mutations (SMM), SM, and EM.


Introduction
Since the industrial uprising, the worldwide energy recoupment has been governed fundamentally by fossil fuels. This has crucial implications for the atmosphere. The expanded employment of RER can assist in the de-carbonization of the energy system in the future. This clean energy can help in dropping the inimical fossil fuel use and energy imports. Hence, it can portray an essential role in garnering the environment green and creating economic evolution.
With the expansion in the share of renewable energy associated with the power grid, the effectual collaboration of the function of different energy sources has emerged as a fresh challenge to the power system scheduling. The integrated operation of the RER-based power generation system may enhance the conflicts between electrical power generation and varying power outputs.
India is a country of geographical diversities. It has a large number of treasures of RER due to its magnificent topographical location. Its huge 'Thar Desert' of Rajasthan has high wind speeds and intense solar radiation. Its vast coasts are a great source of wind, solar, and tidal energies. It has sky-high 'Himalayan Mountain Ranges', which are the origin of thousands of water bodies and forests. Its large plains and plateaus have solar and wind energies, in abundance. To meet the power need and to turn down the role of fossil fuels in the power generation system, a large number of RER-based systems have already been implanted, in the country. Up to December 2019, the total RER-based instated capacity of India was 84 GW, with a set target of inducting 175 GW, by 2022 [1]. Therefore, the deployment of many more RER-based power generating systems is required, so that the rule of fossil fuels can be overthrown and quality power can be delivered to every needy, without much disturbing nature.
In the past, Liaquat et al. [2] have proposed multi-update position criteria for enhancing the investigation characteristics of the traditional firely technique while incorporating the effect of the globally best result on the fluctuation of the fireflies in the exploration zone of the objective function. They have put in the dynamic search space squeezing to compress the fireflies movement inside the definite boundaries to circumvent their oscillatory movement obtained while getting on for the global best solution by finding out the best trajectory for one and all fireflies. Rahimi and co-researchers [3] have elaborated a stochastic thermal and electric load scheduling problem considering the security constraints and also uncertainties of loads, RERs (wind and solar), and market price. They have used a scenario reduction approach to model all uncertain parameters. Naverson et al. [4] have adopted the continuous-time framework to design flexible hydropower sources negotiating with thermal generators (slow-ramping) to minimize the operation cost of the system. They have demonstrated their study through a small-scale case study in which a hydropower plant is connected to a thermal power plant with a manageable high voltage direct current cable.
Narang et al. [5] have applied the predator-prey optimization method for power scheduling of variable/fixed-head hydrothermal system. Predator assists to sustain heterogeneity in the swarm and also avert ill-timed convergence to the localized sub-optimal. They have used the variable elimination technique to control the equality constraint by abolishing variable exactness. Researchers in [6] have employed an adaptive predator-prey optimization to evaluate thermal power scheduling problems in a multi-objective framework. They have maintained the velocity of prey within limits by respectively. Das and co-researchers [23] have evaluated a hydrothermal scheduling problem deploying quasi-reflected symbiotic organisms search. This algorithm has been comprised of symbiotic organisms search to refine the execution of the prescribed technique.
He et al. [24] have used an upgraded combined binary and real number differential evolution technique based upon SHADE to present a model of coordinated power generation scheduling of hydro-thermal-wind system including spinning reserve. They have demonstrated the proposed model with an example and case study. Researchers in [25] have solved a multi-objective economic load dispatch problem with emended salp swarm algorithm. This algorithm includes the solitary and colonial phases of the reproduction cycle of life of salp. They have handled the equality constrained and prescribed functioning zone constraints.
Li et al. [26] have taken a large-scale hydro-wind-solar field in southwestern China to design an optimal power generation scheduling problem. They have maximized the total generated power and the minimum monthly collected output for the entire scheduling interim and minimized the environmental over & short discharge. Panda and Tripathy [27] have employed a new evolutionary hybrid algorithm for environmental optimal power flow problems including wind and thermal power generation systems. They have considered operational cost, emission cost, real power loss, and installation cost of FACTS devices to maintain a stable voltage.
Takahama and Sakai [28] have worked on the -constrained simplex method (ACSM), to solve the constrained optimization problem of the real world. They have instated three modifications in the nonlinear simplex search method to obtain the borderline of the feasible zone, moderate the convergence speed, increase the accuracy, and enhance the overall efficiency of the system. Brar et al. [29] have suggested multi-objective fuzzy satisfying power generation scheduling by using simplex weightage pattern search. They have minimized four contradictory constraints and obtained real and reactive line flows by using generalized Z-bus distribution factors.
In this paper, a futuristic practice described as ACSM is executed to resolve a multi-objective real-time coordinated solar-wind-thermal power scheduling problem. It is a reconditioned unification of the SM introduced by Nelder and Mead. It is an improved conversion technique for constrained optimization. In this method, the non-linear simplex method is perceived as an evolutionary method in which a specific choice, substitution approach, and exceptional variation operator are employed to get high convergence speed, accuracy, and efficiency. It has been invented after hybridizing an established SM with certain other procedures (like-EM, α-constrained method, etc.). To frame this optimization method, three changes in the ordinary SM are executed: (i) α-level comparisons, (ii) the worst point's mutation, and (iii) use of multi-simplexes. In this study, three places from different parts of India are sorted out, where a coordinated solar-wind-thermal power system can operate efficiently. These marked out places are Muppandal (Tamil Nadu), Jaisalmer (Rajasthan), and Okha (Gujarat). A multi-objective coordinated solar-wind-thermal scheduling problem is formulated and optimized for the contemplated sites for two test systems, by using ACSM. To reflect the ascendancy of the suggested operating procedure, the outturns are differentiated with PSO, SMM, SM, and EM.

Optimal problem formulation
In the real world, the coordinated multi-objective optimization problems (CMOP) prerequisite the optimization of many contradictory constraints, concurrently. In this paper, two objectives of thermal and RER systems are discerned. These objectives are total functioning cost and emission (  ,  ,  & ). Therefore, the CMOP is formulated as the minimization of two objectives subject to many equality and inequality constraints. All the objectives are evaluated discretely and then they are solved concomitantly using multi-objective configuration. Two objectives of interest (cost and emission) are of conflicting nature, especially in the case of thermal power generation system. An optimal solution to one can be attained at the cost of the other. Therefore, they are solved simultaneously to achieve the best compromise solution. These objectives can be stated as:

Optimal economy dispatch
The operating cost of a coordinated solar-wind-thermal power generating system depends upon the cost of fossil fuel and the functioning cost including the uncertainty cost of an RER-based power system. This objective can be conceived as the minimization of the total functioning cost of the system. The economy objective of the contemplated system can be examined as [30][31][32][33][34]: where , , and are the number of thermal generators, wind generators, and solar units, respectively.
is the power output of the thermal generator in MW. , and are the cost coefficients of the thermal generator. is the wind power cost of the wind generator.
is the solar power cost of the solar unit.

Environmental objectives
Contrasting with RER-based plants the substantial environmental pollution is originated from the thermal power generation system, which is comprised mostly of nitrogen oxides ( ), sulfur dioxide ( ), and carbon dioxide ( ). In this paper, , , & are specified as the emission determining index and treated as a single objective instead of three objectives. The economy and emission functions can be directly associated through the persistent factor known as emission rate per Mkcal for the defined grade and categorization of fossil fuel. The total thermal emission content is taken as the quadratic functions of thermal power output and can be expressed as [30][31][32][33][34]: where , & are the emission coefficients of the thermal generator and is the emission ( , , & ).

Economic-environmental optimization problem
The goal of the multi-objective coordinated optimization problem for solar-wind-thermal power system is the acquisition of the optimal power dispatch by effectuating the minimization of incongruous objectives, simultaneously. The multi-objective power scheduling problem can be stated as [30][31][32][33][34]: Subject to: i. The equality constraint: The total generated power of the solar-wind-thermal system must be equal to the addition of power demand and transmission losses. Therefore, the load demand equality constraint of the developed problem can be defined as [30][31][32][33]: where is the scheduled power of the wind generator in MW. is the scheduled power of the solar unit in MW.
is the system power demand in MW. is the total system transmission losses in MW. ii. Power generation limits of generating units: The decision variables of thermal, wind, and solar systems ( , , & ) must lie between the power generation limits of the respective generating unit. The lower and upper generation limits enforced on thermal, wind, and solar power generating units are [30][31][32][33]: where and are lower and upper limits of the power output of the thermal generator in MW, respectively. is the rated power output of the wind generator and is the rated power output of the solar unit in MW.

Model of uncertainty of wind generators
Electrical power generation has arisen as the principal implementation of wind energy, globally. This energy renders an accepted contemporary power generation source and a vital participant in the world's energy trade. Wind power generation exceedingly depends upon wind speeds. To procure the precise solution of wind power dispatch prognostication of wind power is decisive. In this paper, the Weibull Distribution Density Factor is used to examine irregular wind data. The wind speed variations are demonstrated by using the Probability Density Function (PDF) and it can be evaluated as [32][33][34][35][36]: where is the annual average wind speed in m/sec. k is the shape factor. c is the scale factor in m/sec. The shape factor is a parameter that displays the span of allocation of wind speeds. It can be obtained as [32,34]: where and are the mean wind speed and the mode wind speed, in m/sec, respectively. The scale factor displays the capability of the wind power of that location. It can be defined as [32,34]: The Gamma Function has frequently applied extension of the factorial functions to the complex numbers and can be observed as [32,34]: Wind energy approximation is decisive to guarantee grid regulation and optimal wind power dispatch. Wind velocity distribution for a specific wind power zone can be designed by using probability distribution functions (PDF). The PDF of wind power can be expressed as [32,34]: where , , & are the cut-in speed, the rated speed, and the cut-out speed of the wind generator in m/sec, respectively. is the operating wind speed in m/sec.
The available wind power at a particular location depends upon the specifications of the wind generator and operating wind speeds during the considered period. The available wind power for the wind generator, at different wind velocities, can be calculated as [32][33][34][35]: Electrical power systems which assimilate RER have to deal with unreliability about the accessibility of load or injected power. This causes the consideration of uncertainty costs in the representation of stochastic economic dispatch. The observation of these costs is vital for the accepted management of RER and the proper issuance of the available energy amount for the power system. The actual cost of wind power is often found more than its anticipated cost. The direct cost function of the wind generator can be evaluated as [32,34]: where is the direct cost coefficient of the wind generator.
When the actual wind power is found less than the planned wind power, the operator has to a pay penalty cost, which is called the overestimation cost. The overestimation cost function of the wind generator is determined as [32,34]: where is the overestimation cost coefficient of the wind generator.
On the other hand, underestimation cost is fine for not utilizing the available wind power for the certain duration. The underestimation cost function of the wind generator can be obtained as [32,34]: where is the underestimation cost coefficient of the wind generator.
The total operating wind power cost of the wind generator is equal to the sum of the direct cost, the overestimation cost, and the underestimation cost, for a specific time. The total operating cost function of the wind generator can be stated as [32][33][34]:

Model of uncertainty of solar units (PV)
Solar energy can be pivotal to the clean energy future. The sun daily radiates far more energy than the power requirements of all the human beings on earth. Solar radiations vary with the topography and climate of a certain area. In this paper, to analyze the irregular solar data, normal distribution is utilized. The Probability Density Function (PDF) of solar irradiance can be calculated as [31,34]: where is the solar irradiance at a given time, is the mean of solar irradiance over the year, and is the standard deviation of solar irradiance, in kWh/m 2 /day.
The available power of the solar unit can be evaluated as [31,33,34]: where is the operating temperature and is the reference temperature of the solar unit, in ºC. is the temperature coefficient in /ºC. is the maximum value of solar radiation incident under standard conditions in kWh/m 2 /day. Solar radiation is a broad expression for the electromagnetic radiation discharged by the sun. These can be seized and converted into useful formations of energy, such as electricity and heat, employing different technologies. The solar radiations incident on an inclined plane is expressed as [31,33,34,37]: where Ø is geographical latitude, A is the angle of the tilt of the solar collector, δ is the sun's declination, and ω is the hour angle, in degrees.
The declination angle of the sun varies seasonally because of the Earth's tilt on its rotation axis and its rotation around the sun. This angle would always be 0° if the Earth were not leaning on its rotation axis. As the Earth is sloped by 23.45° and the angle of declination depends upon this amount [37]. The angle of the sun's declination can be obtained as [31,33,34]: where is the number of the day of the year. Similar to the wind power system, due to the unsure conduct of the sun, the forecasted solar power may not always be equal to the scheduled solar power. The operating cost of the solar power system also depends upon the direct cost and the uncertainty cost of the solar unit. The direct cost function of the solar unit can be determined as [31,34]: where is the direct cost coefficient of the solar unit.
The penalization for deploying another energy resource or for not supplying energy is called overestimation cost as discussed in the wind system. The overestimation cost function of the solar power unit is given as [31,34]: where is the overestimation cost coefficient of the solar unit.
The castigation for not utilizing all the available power or the underestimation cost function of the solar power unit can be obtained as [31,34]: where is the underestimation cost coefficient of the solar unit.
The total operating cost of the solar unit can be evaluated as [31,34]:

Transmission losses
These losses directly rely upon the characteristics of the network and the operation mode. They are mostly caused by energy dissipation in the conductors, appliances used in transmission lines, etc. Regardless of how the electrical power system is modeled, these losses are unpreventable and must be fabricated before depiction can be evaluated. In this paper, the transmission losses of the coordinated solar-wind-thermal power system are obtained by making use of Kron's approximated loss formula through beta-coefficients and these can be expressed as [30][31][32][33][34]: where , , and are the transmission losses due to the thermal system, the wind system, and the solar system respectively, in MW.

Solution Procedure
In this paper, ACSM is applied to optimize the proposed real-time constrained power scheduling problem. In this technique, SM is upgraded by inlaying α-level comparisons rather than ordinary comparisons, mutations of the worst points, and multi-simplexes as a substitute for a single simplex.
The -level comparisons are used to transform algorithms for unconstrained optimization problems into algorithms for constrained optimization problems. By employing these comparisons, the search points are compared based on the pre-defined satisfaction level of their constraints. It means that the points are differentiated on the grounds of their constraint infringement.
The optimal solutions to the constrained problems are frequently found near the frontier of the feasible region. Therefore, to explore the points around the boundary of the feasible zone, the mutations to the worst points and multi-simplexes are adjoined, in the algorithm. It can also control the convergence speed.

Refinements to the nonlinear simplex method
The nonlinear simplex method is upgraded by performing the following alterations to enhance its efficacy [28]: (i) -level collations: To transfigure a constrained optimization problem to an unconstrained optimization problem, -level comparisons can be applied in place of usual comparisons. Consider ( ) is an objective function and ( ) is its satisfaction level. The -level comparisons are the order relation on ( ( ), ( )) set, in which the viability of is more vital than minimizing ( ). -level collations between any two objective functions & having satisfaction levels & respectively and can be stated as follows: The value of lies between 0 and 1. The level collations are analogous to the usual differentiations if the value of is zero. (ii) Incorporation of mutations: In constrained optimization problems usually, the optimal solutions are found very close to the borderline of the feasible province. While solving these problems with the help of the nonlinear simplex method, when simplex reduces some search points encircling the boundary of the feasible zone are ignored, sporadically. Therefore, to evade such situations mutations are included because they are capable of producing optimal results surrounding the frontier of the feasible stretch of the surveyed points. The least desirable point is exchanged by mutations utilizing Eqs (46,47).
Mutations can also sway the consolidation speed of the algorithm. Therefore, the preferable value of the mutation rate should be picked out. High values of mutation rate gravitate towards a large number of calculations and therefore sometimes computation speed may decay. (iii) Inclusion of multi-simplexes: In the nonlinear simplex method, simplex may overlook affine autonomy occasionally and therefore this technique can't search for the optimal solutions. To handle such situations multi-simplexes are included. The affine autonomous simplexes can work for the optimal solutions, even when some simplexes mislay affine sovereignty. In a nonlinear simplex method, for the decision variable of dimension, initially, + 1 points are searched, whereas to formulate multi-simplexes at least + 2 points are generated.
The count of explored points regulates the diversification of the investigation operation and also the simulation speed. If the number is very small but the convergence speed is lofty the surveyed points usually concur to a confined optimum. If the number is very large and the convergence speed is low, the explored points can't arrive at the global optimum.
(2) Randomly generate the initial search points ( (> + 1)) in the extremities of the search zone. Scrutinize the dimension arbitrarily and now either the upper limit or the lower limit is designated to . The remaining variables are created inconstantly betwixt the upper limit and lower limit of each variable.
(3) Determine the feasible solution by employing the designed casual heuristic search approach.
In this paper, two distinct search proceedings are executed auspiciously to attain the feasible solution that delivers the solar-wind-thermal power generations during appeasing the equality constraints (Eq (3)). One strategy is to encounter the accessible power demand constraint of the system over the inspected interim. The second maneuver is based upon the handling of the uncertainty of the RER system using Eq (7) to Eq (26).
and go back to step 4. The contraction factor controls the convergence speed of the computation process. If it is small the investigating process reaches its centroid very soon. If the convergence speed is high the search operation may omit the global optima and coincide with local optima. If the contraction factor is very large the processing speed turns out to be low therefore search may not attain global optima.
The feasible zone can be extended by moderating the value of and the extended feasible zone can be brought down to the primal by increasing the value of to 1, in a constrained SM. The value of regulates the growing speed of and the speed of bringing down the extended feasible zone. The value of reaches 1 moderately if the value of is less. Here, the possibility that the explored points approach to local optimum is less. If the value of is very small the explored points must examine a large area, therefore efficacy reduces. If the extended zone happens to be large at the half iterations, the explored points might reach to feasible zone quickly therefore, it can omit the global optima [28].

Decision-making
The decision-making has an indefinite character and fuzzy targets for the objective functions. The goals comprise categories of alternatives whose limits are not distinctly defined. The fuzzy goals and fuzzy objectives can be defined accurately as fuzzy sets in the zone of substitutes. Here, a fuzzy decision can be observed as the junction of the specified targets and objectives. These fuzzy aims are adjusted by establishing their membership functions, whose values vary from 0-1. The value 0 of the membership function means irreconcilability and value 1 indicates complete complementarity. It can be defined as [30][31][32][33][34]: where is the objective function. and are the maximum and minimum values of the objective function, respectively.
The Fuzzy Cardinal Priority Ranking (the membership function) of the non-dominated (pareto-optimal) solution to a fuzzy set can be stated as [30,33,34]: where K is the number of non-dominated solutions.

Case studies and results
The optimization problems of the electrical power systems (EPS) are very strenuous to solve because the EPSs are very sizeable, composite, structurally extensively distributed, and are impacted by several unpredicted circumstances. Real-time optimization techniques utilize the accessible computations in the optimization structure and are, thus, competent in conducting the appropriate self-optimizing regulation. In this paper, a real-time multi-objective coordinated solar-wind-thermal power scheduling problem is optimized for three different places in India, by using ACSM. The names of these places are: 1. Muppandal (Tamil Nadu) 2. Jaisalmer (Rajasthan) 3. Okha (Gujarat) The same set of power generating units (PGU) and coefficients (cost and emissions) are used for all three places, to access the performance of PGU in the different geographical and environmental conditions. The functions of cost coefficients of the RER power system are also taken as the same for all three sites. This work is executed with the help of the FORTRAN-90 programming language using ACSM and results are collated with PSO, SMM, SM, and EM.

Geographical positions and renewable energy potentials of Muppandal, Jaisalmer, and Okha
India is the seventh-largest country in the world and it lies on the north of the equator between 8°4' north to 37°6' north latitude and 68°7' east to 97°25' east longitude [38]. It possesses large topographical and meteorological variations. Many parts of the country are rich in solar and wind energies. Muppandal, Jaisalmer, and Okha are among such places.
Muppandal is a village in the Kanyakumari district of Tamil Nadu state of India. It is situated at the southernmost point of India. It has India's largest wind farm, with a 1500 MW installed capacity [39,40]. Jaisalmer city of Rajasthan is placed in the northwestern region of India. It is a segment of the 'Great Indian Thar Desert'. India's second-largest wind farm of 1064 MW capacity, is installed here [41]. Okha is a famous town in the Dwarka district of Gujarat state. It is a port at the west-central tip of India. It has a high solar and wind energy prospect. The solar and wind parameters of these three sites are tabulated in Table 1.  Every single locality on Earth acquires solar radiation at the minimum fragment of the year. The quantity of solar radiation that arrives at any one part on the surface of Earth differs according to geographic location, time of day, season, local weather, local landscape, etc. Due to the round shape of Earth, the sun hits the surface at dissimilar angles, varying from 0° to 90° [37]. When the rays of the sun are vertical the surface of Earth achieves all the possible energy. If the sun's rays are more tilted, they travel through the atmosphere for a longer period and therefore, become more diffused and scattered. Figure 2 displays the variation of solar radiation in three considered zones, over the year. It can be seen that all three locations have a high value of solar radiation (except in the monsoon season, from June to August). The maximum drops of radiation can be seen in Okha during this season. Therefore, it has the minimum value of available solar power, in this span of the year.
The mean and mode wind speeds of the year are nearly 11.50 m/sec, 9.0 m/sec & 8.70 m/sec, and 5.50 m/sec, 3.42 m/sec & 2.90 m/sec, respectively for Muppandal, Jaisalmer, and Okha.

Simulation and results
In fact, because of the intermittency of wind and solar radiation, generating extra power results in increasing output fluctuations. Therefore, in this paper, an integrated scheduling model of a territorial RER-based energy system combined with a conventional thermal power generation system is established to encounter the oscillating power requirements of consumers. The considered coordinated solar-wind-thermal power system contains six generating units (two thermal generators, two solar units, and two wind farms). The fuel costs functions and the pollutant emission ( , & ) functions of two thermal generators are given in Table 2  Coefficient of overestimation cost (Rs/kWh) 12.280 The solar system parameters for the 15 th day of June of each year are enlisted in Table 3 and wind system parameters are charted in Table 4. At 1 PM the value of the hour angle is found as −15°. The values of angle of tilt of solar collector and cost coefficients of solar & wind systems are contemplated as the same for the described sectors of India.
The fuel cost and the emissions ( , & ) of the thermal generating system are obtained by using Eqs (1) and (2). The power balance Eq (3) is solved, subject to the equality constraint and power generation limits of thermal generators, wind generators, and solar units, using Eqs (4-6). The wind data is contemplated according to the Weibull distribution density function. The PDF of wind behavior is observed from Eq (7). The shape factor 'k' is evaluated from the mean and the mode wind speeds by using Eq (8), which is found as 2.229, 2.86 & 3.33 for Muppandal, Jaisalmer, and Okha, respectively. The scale factor 'c' is calculated from Eqs (9) and (10), by using the Gamma Function. These values are obtained as 12.981 m/sec, 10.099 m/sec, and 9.695 m/sec for Muppandal, Jaisalmer, and Okha, respectively.   Figure 3 represents the variation of wind speed frequency distributions with a range of wind speeds, at 1 PM on the 15 th day of June of each year, for Muppandal, Jaisalmer, and Okha. It can be seen that all three places have divergent heights and areas of their wind frequency curves because of their distinct wind distributions. The PDF of wind powers is determined from Eqs (11-13) and available wind powers for different locales are observed from Eq (14). The direct cost, overestimation cost, underestimation cost, and total operating wind power cost for all three wind power systems are evaluated from Eqs (15)(16)(17)(18). Depending upon the regional wind distributions for the considered period of specified zones, the available wind powers are observed as 20.869565 MW, 14.3478 MW, and 13.565217 MW for Muppandal, Jaisalmer, and Okha, respectively.
The solar data is examined under the normal distribution of solar irradiance. The PDF of solar radiation is determined with the help of standard deviations and the mean of solar irradiance of three considered places (Figure 4), using Eq (19). The hourly beam solar irradiance incidents on an inclined plane are calculated with the help of Eq (21). The angle of declination of the sun is obtained from Eq (22). The available solar power depends upon solar radiation and the reference temperature of the examined area. Since Jaisalmer has the highest values of both of these variables for the testing interval, therefore it has the highest available solar power (27.39851 MW), which is succeeded by Muppandal (25.26320660 MW) and Okha (24.02259351 MW). The direct cost, overestimation cost, underestimation cost, and total solar power cost are obtained from Eqs (23)(24)(25)(26). Transmission losses of the system are observed from Eqs (27)(28)(29)(30). Fuzzy cardinal priority ranking of non-dominating solutions is employed to obtain the best compromise solution (BCS), with the help of Eqs (44-45). This test system is comprised of power scheduling of solar-wind-thermal of three inspected places of India for 250 MW power demand. Since Muppandal has the maximum available wind power for the given spell, therefore it has the maximum scheduled wind power (40.1869 MW). It is found that Muppandal has the highest direct cost of wind power (97654.1700 Rs/h), which is followed by Jaisalmer 68422.3300 Rs/h and Okha (65499.0300 Rs/h). After optimizing the power scheduling problem, the values of fuel cost, NOx emission, SO2 emission, CO2 emission, total operating cost of wind system, total operating cost of solar system, transmission losses, and simulation time for Muppandal, Jaisalmer, and Okha are obtained as 1765. The solution of power scheduling problem of test system-I by using ACSM is tabulated in Table 5.  Figure 5 represents the comparison of load shared by thermal, solar, and wind powers for Muppandal, Jaisalmer, and Okha. Jaisalmer has the highest PDF of solar radiation for the considered time. Therefore, the solar power generated here is about 21% of the total power generated by the coordinated system. Load shared by solar power in Muppandal and Okha is 20% and 19%, respectively. Muppandal has the lowest value of shape factor (k = 2.229) and highest value of scale factor (c = 12.981 m/sec). The mean wind speed in Muppandal is 11.50 m/s. Therefore, the generated wind power is about 16% of the total power generated by the system, whereas this value is 11% for Jaisalmer and 10% for Okha. The overall load share of RER-based power is found 36% for Muppandal, 33% for Jaisalmer, and 29% for Okha. Therefore, the thermal power generation is the minimum at Muppandal (64%), for the examined period.

Test system II (power demand = 400)
Test system-II comprehends the same set of six generators as used in test system-I but the power demand is increased from 250 MW to 400 MW. Since the parameters of solar and wind systems are not changed, therefore, the available solar and wind powers are also not changed. Now, the extra load is supplied by the thermal generating system. It is therefore the value of fuel cost, emission, emission, and emission also arise.  Table 6.   Table 7. It can be seen that ACSM possesses the highest value of Cardinal Priority Ranking and least values of fuel cost & emissions for both the test systems. It takes the least time to achieve BCS as compared to the other four techniques and also its performance is not affected by using a large number of decision variables.

Comparisons
The box plots in Figure 5

Conclusions
The integration of large-scale RER is of substantial thrust to electrical energy economizing and limiting emissions. The intensive solar and wind power plant constitutes an auspicious alternative source of RER technology. It also acknowledges the amalgamation of thermal power storage for the accumulation of energy for future utilizations, but RER poses multiple obstacles for power systems because of their changeability, uncertainty, and discontinuity. The employment of RER through prudent scheduling and consigning of power can impart operative pliability into the electrical power system.
In this paper, a multi-objective coordinated solar-wind-thermal power scheduling problem is formulated and optimized for two conflicting economic and environmental objectives. ACSM has been successfully employed in the presented non-linear optimization problem and results are contrasted with some other existing popular population-based techniques. ACSM displays significant competence to runoff from local best solutions because the priority is given to the satisfaction level over the value of the objective function, by applying -level comparisons. The addition of mutation of the least wanted point and multi-simplexes enhances the exactness of the technique. Therefore, it lowers the probability of missing out on the points around the boundary, during the reduction of the simplex. It is an effective method for constrained optimization because ACSM remains ineffective even when parameters are reshaping. It is a very quick, lethal, and stable technique for constrained optimization problems.