Non-Dominated Sorting-Based Hybrid Optimization Technique for Multi-Objective Hydrothermal Scheduling

: Short-term hydrothermal scheduling problem plays an important role in maintaining a high degree of economy and reliability in power system operational planning. Since electric power generation from fossil ﬁred plants forms a major part of hydrothermal generation mix, therefore their emission contributions cannot be neglected. Hence, multi-objective short term hydrothermal scheduling is formulated as a bi-objective optimization problem by considering (a) minimizing economical power generation cost, (b) minimizing environmental emission pollution, and (c) simultaneously minimizing both the conﬂicting objective functions. This paper presents a non-dominated sorting disruption-based oppositional gravitational search algorithm (NSDOGSA) to solve multi-objective short-term hydrothermal scheduling (MSHTS) problems and reveals that (i) the short-term hydrothermal scheduling problem is extended to a multi-objective short-term hydrothermal scheduling problem by considering economical production cost (EPC) and environmental pollution (EEP) simultaneously while satisfying various diverse constraints; (ii) by introducing the concept of non-dominated sorting (NS) in gravitational search algorithm (GSA), it can optimize two considered objectives such as EPC and EEP simultaneously and can also obtain a group of conﬂicting solutions in one trial simulation; (iii) in NSDOGSA, the objective function in terms ﬁtness for mass calculation has been represented by its rank instead of its EPC & EEP values by using the NS approach; (iv) an elite external archive set is deﬁned to keep the NS solutions with the idea of spread indicator; (v) the optimal schedule value is extracted by using fuzzy decision approach; (vi) a consistent handling strategy has been adopted to handle effectively the system constraints; (vii) ﬁnally, the NSDOGSA approach is veriﬁed on two test systems with valve point loading effects and transmission loss, and (viii) computational discussion show that the NSDOGSA gives improved optimal results in comparison to other existing methods, which qualiﬁes that the NSDOGSA is an effective and competitive optimization approach for solving complex MSHTS problems.


Introduction
In modern power systems, the efficient dispatch of available power resources for satisfying the increased load has become an important aspect of economic operation of the power system.In this respect, hydro energy sources are considered to be promising alternative energy sources integrated with thermal plants to provide a clean, smooth, environment-friendly energy operation and also meet other requirements such as flood control and irrigation.Usually, the operation of a hydrothermal integrated system is more complex than the scheduling of an all-thermal generation system as the hydroelectric plants are coupled both electrically and hydraulically.Therefore, the optimum scheduling of hydrothermal generation plants is of utmost importance in the current scenario due to its economical aspect linked with the power system operation.
In the process of power system operation and planning, hydrothermal scheduling (HTS) is of a vital role to maintain a high degree of economy and reliability.The storage capacity of the reservoir is the deciding factor in whether long-term or short-term scheduling is desired for the operation planning of a hydrothermal power system.In long-term scheduling, the discharge rates at hydro stations and generation of thermal stations are determined on monthly intervals for an optimization period of several years so that the economical production cost of thermal generation is minimum.Such type of problem is a stochastic one because the load demand and inflows at various reservoirs over the period are not known exactly.On the other hand, the short-term scheduling problem is concerned with an available water allocation over a shorter period such as an hourly period that minimizes the overall thermal generation cost.The inclusion of time delays for water travel time between cascaded reservoirs is usually not necessary for the long-term problem; however, it is mandatory to consider them for short-term hydrothermal scheduling (SHTS) problems.SHTS issues are divided into fixed-head and variable-head problems.The SHTS issue may be divided into two categories based on the head of the reservoir: fixed head and variable head.Fixed head reservoirs have a big capacity, but variable head reservoirs have a restricted quantity of water.Reservoir head is almost constant for high-capacity reservoirs in SHTS problems due to low fluctuation of inflows.Even so, because of the considerable effect of inflows, head variation cannot be overlooked in small scale reservoirs.The hydro generation is assumed to be a function of water discharge quantities in fixed head SHTS formulation.The variable head scheduling issue, on the other hand, is concerned with the nonlinear function of hydropower in terms of reservoir storage and water discharge.Moreover, reservoir volume restrictions like as maximum and minimum limits are generally taken into account in variable head formulation.Variable head plant models also incorporate flow balance equations that influence reservoir dynamics.Therefore, the problem of SHTS scheduling generation with adjustable head hydro plants is relatively a more complicated problem and requires special attention.In addition to the low operating cost of power generation, emission pollution must be lowered for eco-friendly power generation.Both these factors, i.e., cost and emission pollution are conflicting in nature.Therefore, the SHTS problem is reframed to the multi-objective short-term hydrothermal scheduling (MSHTS) subject to satisfy all constraints.
Many investigators have suggested different multi-objective procedures to solve the MSHTS problem over the past decades such as a simulated annealing-based goalattainment approach [1], a fuzzy satisfying method based on evolutionary programming [2], differential evolution [3,4], particle swarm optimization [5], quantum-behaved particle swarm optimization [6,7], a quadratic approximation-based differential evolution with valuable trade-off [8], predator-prey optimization [9,10], self-organizing hierarchical particle swarm optimization technique with time-varying acceleration coefficients [11], hybrid chemical reaction optimization [12], surrogate differential evolution [13], cuckoo bird-inspired meta-heuristic algorithm [14], etc.These entire methods convert the multiobjective optimization problem into a single objective problem by using weights, price penalty factors, or trade-offs.Although these approaches can obtain compromise solutions for the MSHTS problem, some drawbacks of these are (i) there is only one solution produced at each simulation run and the whole optimization process has to be repeated when the weights are changed; (ii) these methods may fail in yielding a uniformly distributed nondominated front when objective functions are non-convex.To overwhelm these difficulties, recently, trade-off approaches have been established such as multi-objective differential evolution algorithm [15][16][17], multi-objective cultural algorithm based on particle swarm optimization [18], hybrid multi-objective cultural algorithm [19], multi-objective artificial bee colony [20], normal boundary intersection method [21], multi-objective quantum-behaved particle swarm optimization [22], parallel multi-objective differential evolution [23], etc.However, these methods are having one other disadvantage premature convergence, more execution time weakly obtained non-dominated solutions, and greater effort is required from decision-makers to find compromise solutions.
Over past few years, the competency of gravitational search algorithm has been proved for solving single-and multi-objective problems .This research article presents a non-dominated sorting disruption-based oppositional gravitational search algorithm (NSDOGSA) for solving single/multi-objective hydrothermal scheduling problems and differs from existing literature based on following points:

•
Initially, an initialization process is escalated over complete search space using the concept of oppositional-based learning.

•
Then, an integrated non-dominated sorting procedure deals with conflicting objective functions and attains a set of non-dominated solutions in a single run.

•
The NS solutions are then stored in a limited length external archive, and a spread indicator metric is incorporated to update the archive.Once more, a disruption operator is used to update each agent's position in order to prevent premature convergence.

•
A fuzzy decision approach is also utilized to select better option from agent set.In addition, an effective constraint handling technique for dealing with load balancing limitations and end reservoir storage volumes in a hydrothermal scheduling problem is described.

•
Finally, NSDOGGSA is carried out on two MSHTS test systems and results are compared with other existing techniques.Thus, the proposed approach attains feasible solutions effectively in terms of compromise solutions with less computational time for solving SHTS/MSHTS problems.
The structure of the paper is arranged as follows: Section 2 presents the problem statement of MSHTS mathematical formulations.Section 3 provides the proposed NSDOGSA methodology for solving MSHTS benchmark problems.Section 4 explains the computational discussion on MSHTS problems trailed by conclusions & future work in Section 5.

Problem Statement
Multi-objective short-term hydrothermal scheduling problem has been formulated to minimize both the power generation cost and emission pollution concurrently, simultaneously apart from fulfilling equality and inequality constraints.The mathematical expressions of objective functions and constraints are mentioned below:

Objective Functions
Traditionally, thermal plants consume fuel for electricity generation, which invites a cost per unit of power generated.Further, the total generation costs for power plants include fuel, labor, and maintenance costs.Unlike capital costs, which are "fixed" (do not vary with the level of output), a plant's total power generation cost depends on how much electricity in that plant produces.In practice, whenever the steam admission valve of the thermal unit starts to open, the economical production cost increases sharply.This phenomenon is called a valve point loading effect which can be represented by the non-smooth generation cost function.The valve point loading effects can be modeled by adding a sinusoidal function to the generation cost function.Consequently, the economical production cost function of thermal power plants with the consideration of valve point loading effects can be described as a quadratic function along with a sinusoidal function, which is represented as follows [2]: where EPC represents cost of economic power generation, Ng indicate thermal plants, T indicates scheduled periods, P t i indicates i th thermal generation at t th a scheduled period, P min i indicates minimum generation limit of i th plant, and a i , b i , c i , d i , e i are i th thermal plant cost coefficients [33].
Moreover, the majority of air pollution is caused by thermal generating units, which are mostly composed of nitrogen oxides NOx, sulphur oxides SO 2 , and carbon dioxide CO 2 , resulting in environmental degradation due to the high concentration of pollutants that they create.The emission pollution functions for a thermal plant can be directly related to the environmental production cost function via the emission pollution rate per megajoule (1 Btu = 1055.06J), that acts as a constant factor for a given type or grade of fuel, yielding quadratic NOx, SO 2 , and CO 2 emission pollution functions in terms of active power generation.Thermal power stations emit the most pollutants into the atmosphere when compared to hydropower plants.As a result, the total environmental emission pollution (EEP) is the mixture of a quadratic and an exponential function, as shown below.[2]: where α i , β i , γ i , η i , δ i indicates coefficients i th thermal unit.

Constraints
The MSHTS is constrained by several hydraulic and thermal factors, such as a load balance constraint, generation boundaries, water discharge boundaries, variable reservoir boundaries, the cascading of the hydraulic system with time delays effect, and starting and ending water storage volumes.All of these constraints have detailed mathematical expressions, which are discussed in [33].

Proposed Methodology for Solving Multi-Objective Short-Term Hydrothermal Scheduling
The gravitational search algorithm (GSA) is receiving more attention and applies to any single objective optimization problem regardless of the linearity, nonlinearity, or non-convexity of its objective function and constraints.However, the GSA has lack of global optimality, as it can often get stuck in local optima.It also tends to be computationally expensive and sensitive to parameter selection.To overcome this inadequacy, the authors have demonstrated different conceptual methods such as oppositional GSA (OGSA), disruption-based GSA (DGSA), and disruption in OGSA (DOGSA) to enhance the performance of GSA convergence for solving single objective SHTS problems [25,26].However, GSA can be used for multi-objective optimization by assigning each agent a different objective function and allowing them to interact with each other.The agents' masses can be used as a measure of the importance of their respective objective functions, and the goal of the optimization would be to find the set of solutions that provides the best trade-off between the conflicting objectives.As a result, attempts were made in the current work towards the creation of a multiobjective approach by considering the fitness value for mass calculation of each particle defined by its nondomination rank rather than the objective function value.The non-dominated sorting procedure yields the non-domination rank.The different steps and mathematical expressions of the proposed technique are mentioned in the subsequent sections.

Population Initialization
In the initialization process of the proposed methodology, the decision variables called as agents are representing the water discharge rates of the hydro plants and power generations of thermal plants at each time interval.Hence, every initial individual is obtained by creating all the decision variables randomly in feasible range with the hydro and thermal power generation outputs, which represents the agents of the population.
Additionally, opposite number may be described as the mirror position of agent from the centre of the search area.Tizhoosh invented opposition-based learning, which is mentioned in [26], opposition-based learning is a novel concept in computer-based techniques, which considers the present agent along with its opposite agent to achieve a better optimal solution.It has been proved in the literature that an opposite candidate solution has a better chance to be closer to the global optimum solution than a random agent solution.Therefore, the random population of agents and opposition-based agents X i & OX i are designated as follows: where X i &OX i indicates i th random and opposite agents, U are i th a random and opposite location at d th dimension, D indicates the overall defined space, and x d L , x d U are lower and upper bounds.In MSHTS problems, over the scheduled time horizon, the water discharge and thermal generation are regarded as separate quantities.q t j = q min j + rand × q max j − q min j P Hence, an initialization structure of MSHTS problem is represented as: , . . ., q T 1 , . . . . . ., q 1 Nh , q 2 Nh , . . .q T Nh , P 1  1 , P A collective population {X, OX} is established by randomly generated agents along with opposite position.After the initialization process, the constraint handling strategies are adopted from [26,29] to satisfy the system balance and end storage volume restrictions for MSHTS issues.

Update External Archive and Calculate Fitness and Mass
The population comprising of the initial and opposite populations {X, OX} is formed and categorized into various non-domination layers using the NS procedure as presented in [32] and the current population is ranked on the fitness value.The agent of a lesser rank is selected and the agents having more crowding distance remain in the current population.This iteration process of the proposed NSDOGSA is carried out for N several agents.The elite external archive is updated using the Pareto dominance concept [32] and the length of the external archive is controlled using the spread indicator metric [27].Further, the NS rank obtained by the non-dominated sorting operation is used to determine the fitness for mass calculation.Figure 1 displays the link between all of the agents' fitness and rank.
and the length of the external archive is controlled usi Further, the NS rank obtained by the non-dominated mine the fitness for mass calculation.Figure 1 display fitness and rank.Fitness value is calculated based on NS rank and are determined given as below [26]: fies the normalized mass, fitness value, worst fitness tion.

Update the Agent's Acceleration, Velocity, and Position
The agent's acceleration in every repetition is com Fitness value is calculated based on NS rank and their corresponding agent's masses are determined given as below [26]: where the normalized mass, fitness value, worst fitness and best of i th agent at k th iteration.

Update the Agent's Acceleration, Velocity, and Position
The agent's acceleration in every repetition is computed as: where, the constant of gravitational is indicates a Euclidean distance between two agents i th and j th at k th genearion, and gbest is the best fitness value.The agent's velocity and the position for the (k + 1) th generation are calculated as below: where rand i indicates a random number between intervals represent the velocity and position [33].

Moving Agents with Disruption Operator
The disruption operator is employed in the oppositional gravitational search algorithm (DOGSA) to avoid the algorithm from becoming locked in local optima and to improve local search [25].The following criterion must be satisfied for disruption of the solutions: The position of the solution satisfying the disruption criterion is updated as follows: where R ij is Euclidean distance between agents i and j, C 0 is constant, C is threshold value, rand −0.25 R ij , 0.25 R ij indicates a random number among the intervals −0.25R ij , 0.25R ij and ρ = 10 −16 is constant value.Thereafter, the NS method in addition to the elite agent's preservation is used for solving the multi-objective problems and fuzzy decision strategy implemented to obtain an optimal compromised result [32,33].

NSDOGSA Procedure for MSHTS Problem
The computational framework of the proposed NSDOGSA method for solving the MSHTS problem is shown in Figure 2 and described as follows: Step 1: Enter the data from the hydrothermal test system and define the NSDOGSA's settings.Create an opposition-based population by randomly initializing the water discharges and thermal power generation outputs for each agent in the population.
Step 2: Constraint handling procedure is performed to satisfy end reservoir volumes and system load balance constraints.
Step 3: Apply the NS strategy and update an external archive group.
Step 4: Update the agent's acceleration, velocity, and position.
Step 5: Moving agent's updating using disruption operator.
Step 6: When prescribed iteration number is reached, stop the operation and store the trade-off solutions in an external archive; otherwise, go to Step 2.
Step 7: Identify the non-dominated solution from the trade-off solutions by using a fuzzy decision strategy.

Discussion and Simulation Studies
NSDOGSA technique as well as other versions have been successfully used to address MSHTS issues on MATLAB 8.5 platform and prepared in PC (core i7-3rd Gen, 2.67 GHz, 8GB RAM).The suggested NSDOGSA approach is demonstrated on two MSHTS problems.The first one consists of 4 hydro and 3 thermal plants without considering transmission losses and another with the same system by including transmission losses.In order to study the feasibility and robust analysis of NSDOGSA, there are two standard test

Discussion and Simulation Studies
NSDOGSA technique as well as other versions have been successfully used to address MSHTS issues on MATLAB 8.5 platform and prepared in PC (core i7-3rd Gen, 2.67 GHz, 8GB RAM).The suggested NSDOGSA approach is demonstrated on two MSHTS problems.The first one consists of 4 hydro and 3 thermal plants without considering transmission losses and another with the same system by including transmission losses.In order to Energies 2023, 16, 2316 9 of 25 study the feasibility and robust analysis of NSDOGSA, there are two standard test systems with different complex properties adopted from [22].The complete hydraulic network and the complete hydrothermal information given in Appendix A. EPC and EEP functions are optimized separately by using DOGSA to enlighten the conflicting relation between these objective functions for MSHTS problems.
It is understood from  7 and 8 shows the optimal hydro water discharges and hydrothermal power outputs attained using DOGSA for separate optimization of EPC and EEP.The convergence characteristics for individual minimization of EPC and EEP with DOGSA approach are revealed in Figure 5.For MSHTS problem, the detailed objective function values of non-dominated solutions obtained with NSDOGSA approach are listed in Table 9.A fuzzy strategy is used to identify viable result from tradeoff solutions and corresponding optimal hydrothermal schedules are presented in Table 10 and Figure 6, which indicates the satisfaction of all kinds of constraints by using the constraints handling method.For MSHTS problem, the detailed objective function values of non-dominated solutions obtained with NSDOGSA approach are listed in Table 9.A fuzzy strategy is used to identify viable result from tradeoff solutions and corresponding optimal hydrothermal schedules are presented in Table 10 and Figure 6, which indicates the satisfaction of all kinds of constraints by using the constraints handling method.

Statistical and Non-Parametric Analysis
Further, the statistical analysis is conducted to show the robustness of the proposed approach for separate minimization of EPC and EEP in terms of the best, average, worst, and standard deviation values, and for MSHTS problem with the highest cardinal priority ranking value which are presented in Tables 11 and 12, for Test system-I and II.From statistical analysis, it is clear that the superlative, middling, and wickedest objective values and standard deviation obtained with DOGSA after 50 runs are better than that obtained with GSA and its variants, which depicts the robustness of the solutions found with DOGSA approach.For MSHTS problem, the highest membership value obtained with the fuzzy decision approach is compared with NSGSA, NSOGSA, NSDGSA, and it is concluded that the NSDOGSA approach attains a greater satisfaction in terms of the quality of the solution.It can also be mentioned that the results obtained with the proposed approach show minimum standard deviation in values, which ensures the precision of the obtained results.To observe the robustness of the proposed algorithm for solving MSHTS problems, a well-known Kruskal-Wallis non-parametric test is conducted for all the test systems.It is a non-parametric method to test the null hypothesis that all the populations have the same median.This test can be used to compare three or more samples and assumed that the populations under comparison need not be normally distributed.Simulation results are obtained from 50 trial runs of the proposed algorithm and compared with the other methods for all the cases mentioned in box plots of Figure 7 for Test system-I and Figure 8 for Test system-II.It is observed that the probability of rejection of the null hypothesis is very low and variation is very small for the proposed method.
Energies 2023, 16, x FOR PEER REVIEW 19 of 24 Test system-I and Figure 8 for Test system-II.It is observed that the probability of rejection of the null hypothesis is very low and variation is very small for the proposed method.

Conclusions
In this paper, the hydrothermal scheduling problem is formulated in three different scenarios, namely: (a) minimizing economical production cost, (b) minimizing environmental emission pollution, and (c) minimizing both the objective functions while considering the water transport delay between hydraulic connected reservoirs, valve point load- Test system-I and Figure 8 for Test system-II.It is observed that the probability of rejection of the null hypothesis is very low and variation is very small for the proposed method.

Conclusions
In this paper, the hydrothermal scheduling problem is formulated in three different scenarios, namely: (a) minimizing economical production cost, (b) minimizing environmental emission pollution, and (c) minimizing both the objective functions while considering the water transport delay between hydraulic connected reservoirs, valve point load-

Conclusions
In this paper, the hydrothermal scheduling problem is formulated in three different scenarios, namely: (a) minimizing economical production cost, (b) minimizing environmental emission pollution, and (c) minimizing both the objective functions while considering the water transport delay between hydraulic connected reservoirs, valve point loading effect of thermal units, and transmission line losses among generation units.The major contribution of the proposed work concludes that (i) an oppositionalbased learning concept during the initiation phase and a disruption operator during the updating phase to achieve diversity and prevent premature convergence; (ii) the concept of a non-dominated sorting procedure has been presented to optimize both objective functions simultaneously to achieve a conflicting solutions in same computational run; (iii) an archive has been formed to store trade-off and a spread metric utilized to update the archive; (iv) a fuzzy strategy is adopted to select the best feasible solution; (v) an effective constraint approach has been formulated to satisfy the load balance and end reservoir storages; (vi) carried out on two MSHTS problems, the first one consisting of four hydro and three thermal plants without considering transmission losses, and another with including transmission losses results have been compared with the other existing approaches, and (vii) statistical analysis has been conducted to show the robustness of the proposed approach for separate minimization of EPC and EEP in terms of the best, average, worst, and standard deviation values, and for MSHTS problems with the highest ranking value, (viii) the simulation results with the proposed NSDOGSA revealed that it is an effective and competitive solution approach for solving complex MSHTS problems.
Further, the areas of the research and the solution techniques presented so far can be further extended to wider regions.In spite of all the above-mentioned contributions, much more intensive work still needs to be carried out, which can be broadly suggested for the following aspects such as (i) the work can be extended to consider the hydro and thermal unit commitment, which involves the optimization of continuous as well as discrete variables during the analysis, (ii) by considering additional complex constraints such as minimum up and down times of generating units, spinning reserve, ramp rates of thermal plants and prohibited zones of hydro plants, (iii) the scheduling problems can be further extended to consider the renewable energy such as tidal energy, wind energy, solar energy, etc., (iv) the present research work can be extended to include pumped storage plants for all scheduling problems, (v) the multi-objective optimization problems can also be handled using other aggregation methods such as ε-constraint method, normal boundary intersection method, lexicographic approach, etc., to generate non-dominated solutions.

Figure 1 .
Figure 1.Relationship between the non-domination rank an

Figure 1 .
Figure 1.Relationship between the non-domination rank and fitness.

Figure 3 .
Figure 3. Performance characteristics of EPC and EEP with DOGSA for MSHTS of Test System-I.

Figure 3 .
Figure 3. Performance characteristics of EPC and EEP with DOGSA for MSHTS of Test System-I.

Figure 4 .
Figure 4. Storage volume of reservoir with NSDOGSA for MSHTS of Test System-I.

Figure 4 .
Figure 4. Storage volume of reservoir with NSDOGSA for MSHTS of Test System-I.

Figure 5 .
Figure 5. Performance of convergence in EPC and EEP with DOGSA for MSHTS of Test system-II.

Figure 6 .
Figure 6.Hydro reservoir storage volumes with NSDOGSA for MSHTS of Test system-II.

Figure 6 .
Figure 6.Hydro reservoir storage volumes with NSDOGSA for MSHTS of Test system-II.

Figure 7 .
Figure 7. Box-type plots of Kruskal-Wallis test for Test System-I.

Figure 8 .
Figure 8. Box-type plots of Kruskal-Wallis test for Test system-II.

Figure 7 .
Figure 7. Box-type plots of Kruskal-Wallis test for Test System-I.

Figure 7 .
Figure 7. Box-type plots of Kruskal-Wallis test for Test System-I.

Figure 8 .
Figure 8. Box-type plots of Kruskal-Wallis test for Test system-II.

Figure 8 .
Figure 8. Box-type plots of Kruskal-Wallis test for Test system-II.

Table 1 .
Performance evaluation of different methods for MSHTS of Test system-I.

Table 2 .
EPC minimization schedule results with DOGSA of Test System-I.

Table 3 .
EEP minimization schedule results with DOGSA of Test system-I.

Table 4 .
Non-dominated solutions for MSHTS of Test System-I achieved using various methods.
Figure 3. Performance characteristics of EPC and EEP with DOGSA for MSHTS of Test System-I.

Table 4 .
Non-dominated solutions for MSHTS of Test System-I achieved using various methods.

Table 5 .
Schedule results for the MSHTS of Test System-I should be compromised with NSDOGSA.
Table 6 that the single EPC optimization: EPC is 42,738.57$and EEP is 20,617.73lb; for single EEP optimization: EPC is enlarged to 49,472.11$ and EEP is declined to 16,939.41lb, and for MSHTS problem compromised with EPC is 43,664.12$and 18,211.05lb.Tables

Table 6 .
Performance evaluation of Test System-II's multiple approaches.

Table 7 .
EPC minimization schedule results with DOGSA for Test system-II.

Table 8 .
EEP minimization schedule results of DOGSA for Test system-II.
Figure 5. Performance of convergence in EPC and EEP with DOGSA for MSHTS of Test system-II.

Table 7 .
EPC minimization schedule results with DOGSA for Test system-II.

Table 8 .
EEP minimization schedule results of DOGSA for Test system-II.

Table 9 .
NS results attained with various techniques for MSHTS of Test system-II.

Table 10 .
Optimal MSHTS results with NSDOGSA for Test system-II.

Table 11 .
Test System-I Statistical analysis with various approaches.

Table A3 .
Reservoir inflows of Test System-I and II.

Table A4 .
MSHTS problem Test System-I and II load demands.

Table A5 .
Thermal generator data of Test System-I and II.