A Multi-Objective Scheduling Optimization Model for a Multi-Energy Complementary System Considering Different Operation Strategies

In order to reduce the amount of abandoned clean energy, the complementary characterization of wind power plants (WPPs), photovoltaic power plants (PVs), hydropower stations (HSs), and thermal power plants (TPPs) combined with energy storage devices (ESDs) is considered, and they are integrated into a multi-energy complementary system (MECS). Firstly, a scenario-generating technique is proposed for uncertainty factors using the Wasserstein method and the improved K-medoids theory. Then, a multi-objective model and solution algorithm are constructed under the objectives of attaining the maximum operation revenue, the minimum abandoned energy cost, and the minimum output fluctuations. Meanwhile, the influence of different ESD operation modes on MECS operation is discussed, specifically, the longest life cycle (LLC) and the optimum economic efficiency (OEE). Thirdly, in order to solve the multi-objective model, a solution algorithm is proposed by using the rough set method to convert the multi-objective model into a single objective model based on the payoff table. Moreover, the complementary features of the MECS are evaluated in terms of the load tracking degree, HS secondary peaking capacity, and units of coal consumption. Finally, the improved IEEE 14-bus system is chosen for the simulation analysis. The results show that (1) the proposed uncertainty simulation method can efficiently generate the most representative scenarios; (2) MECSs can utilize complementary power sources, the OEE mode can better optimize MECS scheduling, and the LLC mode can ensure the ESDs’ life cycles; (3) the scheduling scheme of MECS operation reach the optimal level when the capacity ratio of ESD:WPP–PV iso [0.62, 0.77] in the OEE mode and [1, 1.08] in the LLC mode on a typical summer day, and the ratio is [0.92, 1] in the OEE mode and [1.23, 1.31] in the LLC mode on a typical winter day. Therefore, the proposed model provides effective decision-making support for designing the optimal plan for MECS operation.


Introduction
The rapid development of the Chinese economy has resulted in a sharp increase in energy consumption.However, China's fossil-based energy consumption structure has also brought about a severe energy crisis as well as environmental issues.The optimization of the energy consumption structure has already been shown to be an important method for alleviating the above contradiction [1].Consequently, the importance of wind, solar, and other green renewable energy sources has become increasingly prominent.The Chinese government recently implemented numerous monetary compensation policies to guarantee the sustainable development of renewable energy [2].The cumulative installed capacities of wind power and photovoltaic power had already reached 169 × 10 6 kW and 77.42 × 10 6 kW, respectively in 2016.However, owing to the obvious uncertainty and volatility of wind power and solar power, the retained generation of power should be abandoned to ensure stable system operation [3].In only the first half of 2017, the total abandoned wind and photovoltaic power generation reached 23.5 × 10 9 kW•h and 3.7 × 10 9 kW•h, respectively.Obviously, the solution to the abandoned energy problem has already become an important issue for the promotion of clean energy development and the optimization of the energy consumption structure.
Because to the grid-connected system of uncontrollable power sources brings great challenges to the stable operation of power systems, three available research methods have been adopted by domestic and foreign scholars and workers, namely, improving the forecast accuracy, improving energy storage device (ESD) technology, and utilizing complementary power sources [4].Firstly, improving the forecast accuracy of uncontrollable clean power sources could reduce the output deviation of power generation.Song et al. [5] designed a new forecast model with an advanced intelligent algorithm, a data preprocessing technique, and a non-negative constraint theory.Mojtaba et al. [6] applied the mixture models of Kohonen's self-organizing map and the radial basis functions as the clustering method and the classification method, respectively.Then, ESD could reduce the fluctuations of the power output of intermittent clean energy by utilizing the characteristics of discharging-charging power [7].Kim et al. [8] designed a new distributed energy generation system by integrating renewable power, air energy storage, and natural gas combustion.Atherton et al. [9] analyzed the feasibility of economics and technology when ESDs are applied to a wind farm.Furthermore, Tan et al. [10] discussed the collaborative synthetic effect of ESDs and the demand response in promoting clean energy power generation.The above forecast methods have already been used in practical applications; however, the achievement of high accuracy relies on having an excessively large amount of data, which limits the application space.Meanwhile, ESDs could promote the output of the uncontrollable power source; however, the high investment costs limit the actual application space.Therefore, determining how to utilize complementary but different power sources is key for improving the reliability of the power supply.
Correspondingly, developing an optimal capacity configuration scheme for wind power, photovoltaic power, and conventional power sources could better utilize the complementary characteristics of different power sources [11].Derrouazin et al. [12] developed a multi-input-output fuzzy logic smart controller to achieve the optimal plan for a hybrid power system.Zhang et al. [13] constructed a short-term optimal scheduling model for a hybrid power system with wind power plants WPP and photovoltaic power plants (PV).Hassan et al. [14] put forward a new hybrid power system including WPP, PV, and fuel cell power generation.Mason [15] proposed a capacity configuration model for small-island electricity systems with WPP, PV, and ESDs.Similarly, Ghada et al. [16] proposed a hybrid generation system with WPP, PV, and flywheel energy storage systems.Iman et al. [17] discussed the influence of the synergic operation of WPP, HS, and PV on net power output curve smoothing.Okan et al. [18] put forward a new energy coordinated model for a virtual power plant and discussed the impact of plug-in hybrid vehicles with a high proportion of wind power penetration.Kaabeche et al. [19] constructed an optimization configuration model for a WPP-PV hybrid energy system and analyzed its economy and reliability.These studies have demonstrated the significance of complementary power sources.Aboelsood et al. [20] chose WPP, solar collectors (SC), combined heat, and ESDs to get the optimal operation cost.The above research discussed the optimal capacity planning configuration and coordinated control of hybrid energy systems with different power sources; however, the complementarity of different power sources from the scheduling perspective needs to be discussed further.
According to the above analysis, some research has already investigated the randomness and intermittent nature of uncontrollable power sources with the three proposed available ways.Firstly, a series of advanced intelligent methods were put forward to improve the forecast accuracy and achieve a better practical effect.However, an excessively large amount of data is needed to achieve high accuracy which consumes too much effort.In order to overcome the above problem, the uncertainty could be simulated by a random sampling technique based on the probability distributions of different uncertainty factors.Secondly, some of the literature has discussed the optimization effect of ESDs on multi-energy complementary system (MECS) operation, and different optimization models have been constructed.However, almost all studies have assumed that the ESDs will fully respond to MECS scheduling to achieve the maximum use of clean energy without considering the influences of different operation modes.In fact, ESD could operate in two modes, namely, the longest life cycle (LLC) and the optimum economic efficiency (OEE).The optimal scheduling strategy of MECS operation will show great differences under different operation modes.That means that a discussion about ESDs' operation modes is necessary.Finally, most of the research studied the optimal capacity planning configuration and coordinated control of hybrid energy systems and analyzed the complementarity characteristics of different power sources.However, the studied hybrid systems have been limited to only two or three power sources.In fact, WPP, photovoltaic power plants (PV), hydropower stations (HS), thermal power plants (TPP), and energy storage devices (ESDs) show a great complementarity.Therefore, it is important to construct a new hybrid system with the above power sources and discuss the optimal operation strategy.All the above analyses motivated us to develop a multi-objective scheduling optimization model for MECS operation.The main works of the paper are summarized as follows: • An uncertainty analysis and simulation method are proposed to generate a typical distribution scenario for the uncertainty factors of MECS operation based on the Wasserstein method, and the K-distance and the K-medoids algorithms.The method includes three steps, namely, the discretization of the continuous probability distribution functions, the generation of the initial simulation scenarios, and the selection of the most representative scenarios.

•
A multi-objective scheduling model and solution algorithm are proposed by considering different ESD operation modes under three objective functions, namely, the maximum operation revenue, the minimum abandoned energy cost, and the minimum output fluctuations.Then, the multi-objective model is weighted into the single objective mode by the rough set theory based on the payoff table.

•
A complementary evaluation index system is given to evaluate the optimal degree for the whole MECS operation, including the load tracking degree, the HS secondary peaking capacity, and TPP units of coal consumption.The optimal capacity ratios of ESD:WPP and ESD:PV are calculated by a sensitivity analysis to provide reliable decision-making support.
The rest of the paper is organized as follows: Section 2 presents the basic structure of MECS, and Section 3 constructs an uncertainty analysis and simulation method.Then, a multi-objective optimization model for MECS operation is proposed, and a three-step solution algorithm is constructed based on the payoff table and the rough set theory in Section 4. Section 5 chooses the load tracking degree, HS secondary peaking capacity, and TPP units of coal consumption as the complementarity indexes.Furthermore, the improved IEEE14 bus system is used as a simulation system in Section 6. Section 7 highlights the conclusions of the paper.

MECS Structure Description
WPPs and PVs (or HSs) are significantly complementary.The output of WPPs is low in the daytime, whereas the output of PVs is high.The output of WPPs is high at night-time (or in rainy weather), whereas the output of PVs is low.Similarly, the output of WPPs is lowest in the summer and highest in the winter and spring, and the output of HSs is highest in summer and lowest in winter [21].Therefore, the complementary power generation (CPG) of WPPs and PVs ensures a continuous power supply, whereas the CPG of WPPs and HSs could exploit more wind power and improve the HS peaking capacity.Moreover, TPPs could handle the output fluctuations of WPPs and PVs.The bundling operation of WPPs, PVs, and TPPs could smooth out the output curve.As the short-term operation characteristics of WPPs, PVs, HSs, and TPPs exhibit great differences and a strong complementarity, a sensible joint operation scheduling could result in significant comprehensive benefits.
In this study, a MECS is proposed by integrating WPPs, PVs, TPPs, HSs, and ESDs.For analyzing the complementary of different power sources, the paper assumes all the power sources are independent.Then, the MECS is controlled by an independent operation center (IOC) based on the prediction of weather conditions, which is priority-scheduled to meet terminal customer load demand.The surplus output of the MECS can be transferred to other regions to earn more economic revenue.Moreover, when the available output of MECS cannot meet the load demand, the IOC can purchase power from a generation company (GC).Figure 1 shows the basic structure of the MECS operation.winter [21].Therefore, the complementary power generation (CPG) of WPPs and PVs ensures a continuous power supply, whereas the CPG of WPPs and HSs could exploit more wind power and improve the HS peaking capacity.Moreover, TPPs could handle the output fluctuations of WPPs and PVs.The bundling operation of WPPs, PVs, and TPPs could smooth out the output curve.As the short-term operation characteristics of WPPs, PVs, HSs, and TPPs exhibit great differences and a strong complementarity, a sensible joint operation scheduling could result in significant comprehensive benefits.
In this study, a MECS is proposed by integrating WPPs, PVs, TPPs, HSs, and ESDs.For analyzing the complementary of different power sources, the paper assumes all the power sources are independent.Then, the MECS is controlled by an independent operation center (IOC) based on the prediction of weather conditions, which is priority-scheduled to meet terminal customer load demand.The surplus output of the MECS can be transferred to other regions to earn more economic revenue.Moreover, when the available output of MECS cannot meet the load demand, the IOC can purchase power from a generation company (GC).Figure 1 shows the basic structure of the MECS operation.In the MECS, the WPP and the PV are uncontrolled power sources, whereas the HS, the TPP, and the ESD are controlled power sources.The ESD charges when the load demand is lower than the available power output; otherwise, the ESD discharges (when the load demand is higher than the actual power output).In this study, the maximum benefits, the minimum abandoned energy, and the minimum output fluctuations of the MECS operation are taken as the objective functions for obtaining comprehensive scheduling optimization.

Uncertainty Analysis
In general, the uncertainties of a power system mainly come from two aspects: the power sources and the load demand.For the uncertainties from the load demand side, the load prediction deviation can be described by a normal distribution function.Then, the distribution functions for load demand will also follow a normal distribution function.During the day-ahead scheduling, HS, TPP, and ESD could all provide a reserve service for the load demand.Furthermore, daily weather forecasting technology is mature, so the think power system can better cope with the uncertainty of the load demand.In particular, as smart grid technology becomes more and more advanced, the demand response will be another resource that can be used to overcome the uncertainty from load demand.
The uncertainties from the power source side are mainly caused by uncontrollable clean energy sources, such as wind power and photovoltaic power.The outputs of WPP and PV are not controllable because the natural wind speed and solar irradiance are influenced by the weather In the MECS, the WPP and the PV are uncontrolled power sources, whereas the HS, the TPP, and the ESD are controlled power sources.The ESD charges when the load demand is lower than the available power output; otherwise, the ESD discharges (when the load demand is higher than the actual power output).In this study, the maximum benefits, the minimum abandoned energy, and the minimum output fluctuations of the MECS operation are taken as the objective functions for obtaining comprehensive scheduling optimization.

Uncertainty Analysis
In general, the uncertainties of a power system mainly come from two aspects: the power sources and the load demand.For the uncertainties from the load demand side, the load prediction deviation can be described by a normal distribution function.Then, the distribution functions for load demand will also follow a normal distribution function.During the day-ahead scheduling, HS, TPP, and ESD could all provide a reserve service for the load demand.Furthermore, daily weather forecasting technology is mature, so the think power system can better cope with the uncertainty of the load demand.In particular, as smart grid technology becomes more and more advanced, the demand response will be another resource that can be used to overcome the uncertainty from load demand.
The uncertainties from the power source side are mainly caused by uncontrollable clean energy sources, such as wind power and photovoltaic power.The outputs of WPP and PV are not controllable because the natural wind speed and solar irradiance are influenced by the weather conditions.The output of HS is uncertain from the perspective of the seasonal characteristics, and it is controllable from the perspective of the complementary energy characteristics.If the TPP output is sufficient to offset the insufficient output of HS, the uncertainty of the HS output can be ignored.Therefore, this paper mainly discusses the uncertainties due to WPP and PV.

Uncertainty Model
Based on the above analysis, only the uncertainties from WPP and PV are considered.In order to simulate the uncertainties, suitable probability distribution functions are needed.Once the simulation results of the wind speed and solar irradiance have been attained, the power outputs of WPP and PV can be calculated. (

1) WPP Uncertainty Model
The output of WPP has a high level of uncertainty.The simulation of wind speed is the key to calculating the output of WPP.It has been proven that the Rayleigh distribution function can be used to simulate wind speed as follows: where v t is the wind speed at time t.ϕ t and ϑ t are the shape factor and the scale factor at time t.Then, the WPP output can be calculated by Equation (2), as follows: where g * WPP,t is the available output of the WPP at time t.g R is the rated output of WPP.v t is the real-time wind speed at time t.v in , v out and v rated are the cut-in, cut-out, and rated speeds, respectively.
(2) PV uncertainty model Similarly, existing research has proven that the beta distribution function can be used to describe solar irradiance as follows: where θ t is the solar irradiance at time t.α t and β t are the shape parameters of the beta distribution at time t, which can be calculated by the mean value and the standard deviation of irradiance [15].After simulating the solar irradiance, the PV output can be calculated as follows: where g * PV,t is the available output of PV at time t.η PV and S PV are the efficiency and total area, respectively.θ t is the solar radiation intensity at time t.

Uncertainty Simulation
The PDFs of the WPP output, PV output, and load are continuous probability distribution functions (CPDFs).In general, discrete PDFs are required for MECS scheduling.Therefore, a discretization process for continuous PDFs must be completed before the uncertainty simulation.Thus, the uncertainty simulation comprises three steps: PDF discretization, scenario simulation, and scenario reduction.
(1) PDF discretization The discretization process involves dividing the continuous function curve into a number of segments corresponding to certain intervals.In general, these intervals can be equal or unequal.In the former case, all scenarios can be considered, and, in fact, the probabilities of some scenarios may not reflect the real probabilities.In the latter case, the probability can be captured more accurately.In this study, the Wasserstein method is adopted for the discretization of PDFs.More details on this method may be found in Reference [22].
It is assumed that f c (x) is a CPDF to be discretized using N points.Then, the optimal discrete points z n (n = 1, 2, . . ., N) under the Wasserstein distance indicator are obtained by Then, the discrete probability of z n is calculated as follows: where p n d represents the discrete probability of z n .When n = 0, p 0 d represents the lower bound of the variable x, and its value is represents the upper bound of x, and its value is equal to By combining Equations ( 1) and (3) with Equation ( 6), the optimal discrete points can be calculated.Correspondingly, the discrete PDFs for the outputs of the WPP, PV, and load can be obtained for the scenario simulation.Figure 2 shows the discretization process (6 discrete points).Thus, the uncertainty simulation comprises three steps: PDF discretization, scenario simulation, and scenario reduction.
(1) PDF discretization The discretization process involves dividing the continuous function curve into a number of segments corresponding to certain intervals.In general, these intervals can be equal or unequal.In the former case, all scenarios can be considered, and, in fact, the probabilities of some scenarios may not reflect the real probabilities.In the latter case, the probability can be captured more accurately.In this study, the Wasserstein method is adopted for the discretization of PDFs.More details on this method may be found in Reference [22].
It is assumed that ( ) f x is a CPDF to be discretized using N points.Then, the optimal discrete points ( ) under the Wasserstein distance indicator are obtained by ( ) ( ) Then, the discrete probability of n z is calculated as follows: ( ) where n d p represents the discrete probability of n z .When n = 0, 0 d p represents the lower bound of the variable x, and its value is ( )


. By combining Equations ( 1) and (3) with Equation ( 6), the optimal discrete points can be calculated.Correspondingly, the discrete PDFs for the outputs of the WPP, PV, and load can be obtained for the scenario simulation.Figure 2 shows the discretization process (6 discrete points).(2) Scenario simulation After the PDFs of the WPP output, PV output, and load have been discretized, the simulation scenarios can be obtained by a sampling method, such as Monte Carlo or Latin hypercube sampling.Monte Carlo is highly versatile and accurate.When the distribution has low probability regions, small samples are generally not representative.The Latin hypercube sampling method divides the PDF into several equal intervals to obtain the sample from each interval [12].However, the probabilities of some scenarios may not reflect the real probabilities [23].Based on the above (2) Scenario simulation After the PDFs of the WPP output, PV output, and load have been discretized, the simulation scenarios can be obtained by a sampling method, such as Monte Carlo or Latin hypercube sampling.Monte Carlo is highly versatile and accurate.When the distribution has low probability regions, small samples are generally not representative.The Latin hypercube sampling method divides the PDF into several equal intervals to obtain the sample from each interval [12].However, the probabilities of some scenarios may not reflect the real probabilities [23].Based on the above analysis, in this study, the Wasserstein method is used to divide the continuous probability distribution curves.Then, the scenario tree method is used to describe the different states of power output, namely, the high, mean, and low (z = 1, 2, 3, respectively).In Equation ( 5), the optimal discrete points for the WPP and PV outputs are calculated and the simulation interval of the PDF is accordingly divided.The simulation scenario set is gained by using the interval simulation method.Then, we assume g z t is the expected output in each state at time t.Equation (6) yields the probability p z t of simulation scenarios in each state at time t.Consequently, the scenario set and probability are S = { g z t |t = 1, 2, . . ., T} and p = ∏ t∈T p z t , respectively.More details on the scenario simulation can be found in [24].
(3) Scenario reduction After the continuous distributions of the random variables have been optimally discretized in each period, the scenario set formed by the entire scheduling interval is large, which results in dimension explosion.If all scenarios are applied to make the optimal scheduling scheme, the effort required to solve the proposed mode will be large.The number of simulation scenarios will determine the number of model iteration calculations.If the scenarios are not reduced, some approximate repeated calculations will occur when some scenarios are similar.However, if the number of scenarios is reduced and the most representative scenarios are chosen, then the scenario reduction will decrease the calculation effort.Based on the above analysis, some scenarios with high similarity should be deleted.The main theory of scenario reduction involves deleting similar scenarios by comparing each scenario with the other scenarios.Correspondingly, this paper uses the K-distance method to measure the distances of different scenarios to delete redundant scenarios [25].It is assumed that d u s , u k is the distance between scenarios u s and u k in scenario set S. If scenario set K is used to replace the original scenario set with k scenarios, d u s , u k is calculated in Equation ( 7): Then, we let q k and p s be the probabilities of scenario k and scenario s, respectively.The optimal scenario set K should follow the same principle, namely, the probability information should be closest to S.
In traditional K-medoids methods, scenario u s is compared with each scenario u k to order the group of scenarios closest to u s , which requires a large amount of calculation [26].In this study, the relative distance across which the ordering is performed is constructed through a comparison with the same reference point u 0 u 0 = 0 .The detailed process of sorting and grouping can be found in Reference [27].After the group order has been obtained, u 0 should be updated by where u l is the updated clustering center.Then, the updated ∑ and if the value meets the requirements, the final reduced scenario set is obtained.Otherwise, the updated clustering center needs to be updated again.To improve the calculation speed, the simultaneous backward technique is applied to the scenario reduction process.Figure 3 shows the uncertainty simulation flowchart.
is calculated, and if the value meets the requirements, the final reduced scenario set is obtained.
Otherwise, the updated clustering center needs to be updated again.To improve the calculation speed, the simultaneous backward technique is applied to the scenario reduction process.Figure 3 shows the uncertainty simulation flowchart.

Multi-Objective Scheduling Model
To satisfy different operational requirements and consider different benefits, the maximum operation revenue, the minimum abandoned energy cost, and the minimum output fluctuations were chosen as the MECS operation objectives in this study.

Objective Functions
MECSs can utilize the complementarity of different power sources, make full use of the adjustable power sources' peaking capacities, and promote the grid connection of clean energy sources, such as wind and photovoltaic power.Therefore, the MECS operation revenue comprises the economic revenue, environmental revenue, and peak shaving revenue.
where t is the index for time.k is the index for the scenario.k γ is the probability of scenario k.
, and , G t R denote the operation revenue for the WPP, the PV, the ESD, the HS, and the TPP at time t, respectively.The operation revenue for the WPP, the PV, and the ESD have already been studied in Reference [10].As is the case with the WPP and the PV, the cost of the initial investment, operation, and maintenance of the HS is classified as the sunk cost, whereas the energy consumption cost of the HS can also be ignored.Therefore, the revenue of the HS operation equals the product of the power price and power output.Thus, the revenue of the TPP operation comprises three parts: the power generation cost, the startup-shutdown cost, and the pollutant emission cost.The detailed calculations are as follows: ( )

Multi-Objective Scheduling Model
To satisfy different operational requirements and consider different benefits, the maximum operation revenue, the minimum abandoned energy cost, and the minimum output fluctuations were chosen as the MECS operation objectives in this study.

Objective Functions
MECSs can utilize the complementarity of different power sources, make full use of the adjustable power sources' peaking capacities, and promote the grid connection of clean energy sources, such as wind and photovoltaic power.Therefore, the MECS operation revenue comprises the economic revenue, environmental revenue, and peak shaving revenue.
where t is the index for time.k is the index for the scenario.γ k is the probability of scenario k.R WPP,t , R PV,t , R ESS,t , R HS,t , and R G,t denote the operation revenue for the WPP, the PV, the ESD, the HS, and the TPP at time t, respectively.The operation revenue for the WPP, the PV, and the ESD have already been studied in Reference [10].As is the case with the WPP and the PV, the cost of the initial investment, operation, and maintenance of the HS is classified as the sunk cost, whereas the energy consumption cost of the HS can also be ignored.Therefore, the revenue of the HS operation equals the product of the power price and power output.Thus, the revenue of the TPP operation comprises three parts: the power generation cost, the startup-shutdown cost, and the pollutant emission cost.The detailed calculations are as follows: where π pg TPP,t , π ss TPP,t , and π pe TPP,t denote the power generation cost, startup-shutdown cost, and pollutant emissions cost of the TPP at time t, respectively.a TPP , b TPP , and c TPP are the cost coefficients of the TPP power generation.ρ coal is the coal price.u TPP,t and D TPP,t are the status and cost of the TPP startup-shutdown at time t.N cold TPP and N hot TPP are startup costs of the TPP operation in the cold mode and hot mode.T min TPP and T off CGT,t are the minimum allowable and continuous downtime of TPP at time t, respectively.T cold TPP is the required time for TPP startup in cold-start mode.a TPP , b TPP , and c TPP are the cost coefficients of TPP power generation.
where, F 2 is the abandoned energy cost of WPP, PV, and HS during the scheduling period.For MECS, determining how to use the complementary characteristics of different power sources is the key to reducing the amount of abandoned energy.In fact, if the abandoned energy could be absorbed, more economic revenue could be gained.Correspondingly, this part of the economic revenue can be regarded as the opportunity cost.ρ WPP,t , ρ PV,t , and ρ HS,t are the power prices of WPP, PV, and HS, respectively, at time t.g * WPP , g * PV , and g * HS,t are the available outputs of WPP, PV, and HS at time t, respectively.
where g av is the average value of the MECS power output fluctuations.For the MECS, the uncertainties of WPP and PV will bring a great level of uncertainty in the net power output.Therefore, determining how to measure the power output fluctuations is key to guaranteeing the stable operation of the MECS.In the MECS, ESD could charge the power of WPP and PV in the valley period and discharge their power in the peak period.Therefore, the standard deviation of the total net power output of WPP and PV is taken as the measure indicator.g dis ESD,t and g ch ESD,t are the power levels of ESD discharging-charging at time t, respectively.F 3 is the standard deviation of the power output fluctuations of the MECS.When F 3 is small, the output fluctuations of the WPP and PV are small, which is more favorable for the safe operation of the MECS.

Constraint Conditions
To ensure the safe and stable operation of the MECS, some necessary constraint conditions should be considered, for example, the load supply and demand balance constraints, the TPP operation constraints, the ESD operation constraints, and the system reserve constraints.A detailed description of the MECS operation constraints is as follows.
(1) Supply and demand balance constraint where ϕ WPP , ϕ PV , ϕ HS , and ϕ TPP are the power loss rates of the WPP, PV, HS, and TPP, respectively.g GC,t is the exchange power from the GC at time t.When g GC,t is positive, the MECS will purchase power from the GC to satisfy its power demand.When the g GC,t is negative, the MECS will transfer surplus power to the utility grid.
(2) HS power output constraint where ε is an error term calculated by history data, ρ 0 HS , ρ 1 HS , ρ 2 HS , ρ 3 HS , ρ 4 HS , and ρ 5 HS are the power output coefficients of the HS.Q HS,t and V HS,t are the output capacity and the reservoir capacity of HS at time t.
where Q min HS is the minimum water flow for the HS power output.Q max HS is the maximum water flow for the HS power output.V min HS is the lower bound of the reservoir capacity.V max HS,t is the upper bound of the reservoir capacity.The reservoir capacity of the HS at the initial and end time points of the scheduling period should meet the following constraint: where Q t is the water inflow at time t, Q t is the water outflow at time t.S HS,t is the amount of wastewater at time t, and H is the number of upstream hydropower stations.τ h is the delayed water flow between the HS and the upstream hydropower station h.V T HS is the water storage capacity at time T. V 0 HS is the water storage capacity at the initial time point.Water runoff into the reservoir comprises the water supply and the outflow from the upstream reservoir, as follows: where HS−1 and S τ−t up HS−1 are the power generation flow and wastewater flow, respectively, of the upstream hydropower station.t up is the time lag corresponding to the flow from the upstream hydropower station to the downstream hydropower station. (

3) ESD operation constraints
The ESD can complementarily operate with the WPP and the PV by utilizing its discharge-charge cycle.Two operation modes exist, namely, the longest life cycle (LLC) mode and the optimum economic efficiency (OEE) mode.Thus, the main constraints are as follows: where g ch ESS,t and g dis ESS,t are the upper bounds for ESS charging-discharging at time t, respectively.
where ρ dis ESS,t and ρ ch ESS,t are the power loss rates of the ESD during discharge and charge, respectively, at time t.The main objectives of LLC are to optimize the ESD's discharge-charge cycles, improve the ESD's operation life, and simultaneously increase the system's peak shaving capacity.That is, where g chr,R ESS,t and g dis,R ESS,t are the rated power levels of the ESD charge and discharge, respectively.(4) System reserve constraints: where g max TPP,t and g max HS,t are the maximum outputs of the TPP and the HS at time t, respectively.g min TPP,t and g min HS,t are the minimum outputs of the TPP and HS at time t, respectively.r i (i = 1, 2, 3, 4, 5) represents the up reserve coefficients of the load, WPP, and PV and the down reserve coefficients of the WPP and PV, respectively.
(5) Other power source constraints In the MECS, the other power sources are the TPP, the WPP, and the PV.The operation constraints of the WPP and the PV should satisfy the maximum output constraints, namely, the outputs of the WPP and the PV should not be greater than the system's allowable output, which has already been studied in Reference [21].The constraints of the TPP operation comprise the maximum/minimum output constraints, the upwards/downwards climbing constraints, and start-shutdown time constraints, which have already been studied in Reference [28].

Multi-Objective Model Solution
We proposed a multi-objective optimization model for the MECS operation.There are two methods for solving such a model.The first is to solve it using group intelligence algorithms, for example, genetic algorithms, particle swarm algorithms, and ant colony algorithms.As an intelligence algorithm can globally search for the optimal solution, the programming code is excessively complex, particularly in mathematical models with integers and continuous and discrete variables, thus yielding non-optimal results [10].
The second method is to weight multi-objective functions into an integrated single-objective function.In general, the weighted coefficients can be obtained subjectively or via weight determination methods for different objective functions [28].In a previous study, the weight coefficients were mainly calculated by using the entropy weight method [29].The applicability of this method is limited by the uncertainty factors and imperfect information existing in the MECS operation.To overcome this problem, this paper uses the rough theory to handle the information in the proposed multi-objective model.The calculation process comprises three steps, namely, the construction of the payoff table, the calculation of the weight coefficients, and the conversion of the multi-objective functions.

Payoff Table
The payoff table mainly reflects the relationships between different objective functions [28].The values of the objective functions were calculated when the MECS operated under other optimization objectives.Firstly, the optimal scheduling schemes were calculated as individual optimization functions.Then, the values of the other objective functions F i , i = 1, . . ., i − 1, . . ., i + 1 . . ., I were calculated under the optimized objective functions F i , i = 1, 2, . . ., I. Assume F ij , j = 1, 2, . . ., I is the values of the other objectives under F i , and F * ii is the optimal value of objective F i .Then, the payoff table can be gained.Table 1 is the payoff table of the proposed multi-objective model.1, the minimum value and the maximum value can be gained for different objective functions.Correspondingly, the information needed to calculate the weight can be collected.Since different objectives have different magnitudes, dimension reduction is necessary before calculating the weight coefficients.The fuzzy satisfaction theory can be used to estimate the optimization degrees of different objective functions [29].In general, the objective functions exist in two optimization directions, namely, the maximum and the minimum.Therefore, this paper applies the decreasing gradient half membership functions and increasing half line membership functions to handle the objectives.
When the optimization direction of the objective function is in the minimal direction: When the optimization direction of the objective function is in the maximal direction: where F i min is the minimum value of the objective function F i (•).F i max is the maximum value of the objective function F i (•).f i is the membership function.

Weight Calculation
Rough sets can be used to process incomplete and inaccurate data by learning, induction, and excavation and can be transformed into relatively clear and concise data systems to support decision analyses.The detailed theory about the rough set theory is detailed in the literature [26].When the rough set theory is used to calculate the weight, the weight problem is transformed into a multi-attributed evaluation problem.Correspondingly, the knowledge dependency and attribute importance evaluation methods were constructed for the weight calculations.The detailed process is as follows: (1) Construction of the relationship data model It is assumed that the weight coefficient of f i is 1/I.Then, the value of the integrated single objective F can be calculated.F is regarded as a decision attribute, and D = F is regarded as a decision attribute set.u k is the optimal value of the comprehensive objective, which is an information point of the study object F. Obviously, u j = ( f 1j , f 2j , . . ., f mj ; F j ).U = u 1 , u 2 , . . ., u j is the main domain, namely, the sample set.Therefore, the attributes of u k are f i (u j ) = v ij and F i (u j ) = F j .It is assumed that R V and R D are knowledge bases.
(2) Calculation of the dependence degree of R V on R D , that is, where ρ(•) represents the set base, and p v (D) is the knowledge of all individuals using the classification U/C in the universe U, which can correctly classify the equivalence classes of object sets in the U/D.
(3) Calculation of the dependence degree of where p V−|v i | (D) is the knowledge of all individuals using the classification U/C in the universe U after removing the indicator V i .
(4) Calculation of the weight value of objective i, namely, where σ D (D) calculates the importance of objective function I. λ i is the weight of objective i.

Weighted Single Objective
After getting the weights of the different objective functions, the integrated single objective function can be gained, as shown in Equation (39): where λ i and λ i are the weights of the objective functions with the maximum optimization direction and the minimum optimization direction.The integrated single objective function F is thus obtained in Equation (40).

Complementarity Evaluation Indexes
MECSs can make full use of the complementarity of different power sources and the peaking capacity of TPP and HS to improve the grid connection of clean energy.In this section, the complementarity evaluation indexes of the MECS are constructed.Moreover, complementation scenarios with different power combinations are considered so that the comprehensive operation benefits of MECSs can be comparatively analyzed.In a previous study [26], the load tracking degree, HS secondary peaking capacity, and TPP unit coal consumption were used as evaluation indexes.
(1) The load tracking degree (LTD) To fully utilize the complementarity of different power sources, the load tracking degree (LTD) is used to evaluate whether the complementary ratio of power sources is optimal.LTD shows the degree of dispersion of the revised load curve compared with the average value, which can reflect the smoothing effect of the complementary functions of different power sources on the load demand.When the LTD is closer to 1, the capacity of the bundled power output to track the load change is higher and the load curve is smoother [26].The LTD is defined as where L t represents the load required to absorb the power output of the MECS at time t, L av t represents the average load, G MECS,t represents the total output of the MECS at time t, and G av MECS,t represents the total average output.The load curve may become smoother by optimizing G MECS,t , thus making the output curve trend more consistent with L t .
(2) The HS secondary peaking capacity (SPC) HSs are characterized by low power generation emissions and high regulation, which allow their joint operation with WPPs and PVs.Except for the HS used for the bundling complementary of WPP and PV, the scheduling plan of the other units should be arranged by considering the HSs' adjustment capability.HSs could be used to perform secondary peaking on the revised load curve, which could promote the load supplied by TPP to be as smooth as possible.This could reduce the start-stop number of TPP and improve the overall operating efficiency.The net load curve for the TPP would be gentler by considering the HSs' adjustability and, thus, improving the operational efficiency.To distribute the residual load in each period as evenly as possible, the peaking mean square deviation is used to describe the HS peaking capacity.It is defined by where L2 t is the surplus load after removing the HS output, which is equal to (L t − g HS,t ), and Lav t is the average surplus load after removing the HS output.In general, the load demand curve is gentler when ϕ is smaller, thus reducing the start-shutdown frequency and coal consumption.
(3) TPP unit coal consumption The MECS can fully utilize the regulation characteristics of the hydropower and thermal power units to promote the grid connection of clean energy sources.The large-scale grid connection of clean energy sources could reduce the coal consumption of TPP; correspondingly, the reasonable arrangement of the left load in TPP could reduce coal consumption, which could reduce pollution emission and power generation costs.The power generation unit coal consumption rate γ is used to evaluate the operational benefits of the MECS and is defined by

Basic Data
The improved IEEE 14-bus system with four thermal power units was used as the simulation system [28].A WPP, a PV, an HS, and ESDs were connected with the simulation system to analyze the complementarity of the different power sources.The rated capacities of the WPP, the PV, and the HS were 900 MW, 400 MW, and 600 MW, respectively.Moreover, 400 MW ESDs were considered to analyze the optimal effects in the OEE and LLC modes.The rated discharge-charge power of the ESDs was 100 MW, the maximum discharge-charge power was 150 MW, and the initial energy storage was zero.The total maximum and minimum loads on a typical load day were 18,500 MW and 1400 MW, respectively, in summer and 1800 MW and 1350 MW, respectively, in winter.The loss of the tested system was 5.5%.The operating parameters and pollutant emission parameters of thermal power units can be found in Reference [25].Currently, most thermal power tariffs have already included SO 2 , NO x , soot, and other pollutant costs in China; thus, the carbon cost was primarily considered in this study, and the price of CO 2 emission was set to ¥55/ton [29].Table 2 shows the operation parameters of the TPPs.The grid connection prices of the WPP, PV, HS, and TPPs were ¥610/MW•h, ¥1000/MW•h, ¥320/MW•h, and ¥450/MW•h, respectively [26].The average price of customer power consumption was ¥590/MW•h [8].It was assumed that the ESDs enjoyed real-time pricing.The division of the peak-valley period and peak-valley electricity price can be found in Reference [10].It was also assumed that the coal price was ¥500/tce.Load demand curves for a typical summer day and a typical winter day were chosen as the input data, and the water resources situation on a typical load day was chosen as the HS's input data [21].Figure 4 shows the load demand and the HS's available output on a typical load day.
The parameters of the wind turbine were v in = 3 m/s and v rated = 14 m/s [21].The shape parameter and scale parameter on a typical load day were 3.8 and 9.4, respectively, in summer, and 5.2 and 11.8, respectively, in winter.In accordance with the literature [26], the illumination intensity parameters α and β on a typical load day were 0.5 and 9.12, respectively, in summer, and 0.3 and 8.54, respectively, in winter.Then, the continuous probability distribution curve of wind speed and solar irradiance was obtained.As this study was primarily concerned with the complementarity of different power sources and because the exact method for selecting the number of discrete points was irrelevant here, ten discrete points were used for both curves.The proposed uncertainty simulation method was applied to obtain 50 groups of simulation scenarios for wind speed and solar irradiance.Then, the improved K-medoids method was used to reduce the number of scenarios, and the outputs of ten typical scenarios were calculated.Subsequently, the average outputs of the different scenarios were taken as the input data.Figure 5 shows the simulation scenarios for WPP and PV on a typical summer day. Figure 6 shows the average output of the WPP and the PV on typical load days.
After this basic data had been input, the model was solved by the GAMS software application using the CPLEX 11.0 linear solver from the ILOG_solver.The proposed model involved a nonlinear objective function and constraints belonging to the class of mixed integer nonlinear programming problems.As GAMS is better suited to mixed integer programming problems, the objective function and the constraints were linearized [28].Figure 7 shows the optimization flowchart for the MECS scheduling model.The grid connection prices of the WPP, PV, HS, and TPPs were ¥610/MW•h, ¥1000/MW•h, ¥320/MW•h, and ¥450/MW•h, respectively [26].The average price of customer power consumption was ¥590/MW•h [8].It was assumed that the ESDs enjoyed real-time pricing.The division of the peak-valley period and peak-valley electricity price can be found in Reference [10].It was also assumed that the coal price was ¥500/tce.Load demand curves for a typical summer day and a typical winter day were chosen as the input data, and the water resources situation on a typical load day was chosen as the HS's input data [21].Figure 4 shows the load demand and the HS's available output on a typical load day.
The parameters of the wind turbine were 3m / in v s = and 14 m / rated v s = [21].The shape parameter and scale parameter on a typical load day were 3.8 and 9.4, respectively, in summer, and 5.2 and 11.8, respectively, in winter.In accordance with the literature [26], the illumination intensity parameters α and β on a typical load day were 0.5 and 9.12, respectively, in summer, and 0.3 and 8.54, respectively, in winter.Then, the continuous probability distribution curve of wind speed and solar irradiance was obtained.As this study was primarily concerned with the complementarity of different power sources and because the exact method for selecting the number of discrete points was irrelevant here, ten discrete points were used for both curves.The proposed uncertainty simulation method was applied to obtain 50 groups of simulation scenarios for wind speed and solar irradiance.Then, the improved K-medoids method was used to reduce the number of scenarios, and the outputs of ten typical scenarios were calculated.Subsequently, the average outputs of the different scenarios were taken as the input data.Figure 5 shows the simulation scenarios for WPP and PV on a typical summer day. Figure 6 shows the average output of the WPP and the PV on typical load days.
After this basic data had been input, the model was solved by the GAMS software application using the CPLEX 11.0 linear solver from the ILOG_solver.The proposed model involved a nonlinear objective function and constraints belonging to the class of mixed integer nonlinear programming problems.As GAMS is better suited to mixed integer programming problems, the objective function and the constraints were linearized [28].Figure 7 shows the optimization flowchart for the MECS scheduling model.

Weighting Calculation
It is necessary to obtain suitable weight coefficients for the objective functions before solving the multi-objective model.Firstly, the proposed model was solved with individual objective functions and the values of all objectives were gained.Then, the payoff table for the MECS operation was gained, as shown in Table 3.According to Table 3, taking the scheduling results on a typical summer day as an example, in order to distinguish the different optimal degrees, the value of the objective is set as 1 when the value of the objective after fuzzy processing belongs to (0,1).The values of the objective are set as 0 and 2 when the values of the objective after fuzzy processing equal 0 and 1, respectively.Then, the importance of the indexes can be obtained through a comparison of the results with the index and without the index.In the above evaluation problem, six problems exist (F 1 , F 2 , F 3 in the OEE mode and LLC mode).We choose F 11 , F 12 , F 13 and F 21 , F 22 , F 23 to replace the above six indexes.At this time, the importance of the indexes is high when the change of the sample is obvious if the index is deleted.Then, the problems can be divided as follows: {{F 11 }, {F 12 }, {F 13 }, {F 21 }, {F 22 }, {F 23 }}.
Additionally, the information of the index set is as follows: Then, if F 1 is removed, the equivalence relationship of the indexes in different problems can be divided as follows: {{F 11 , F 21 , F 22 }, {F 12 }, {F 13 }, {F 23 }}.
At this time, the information of the index set is as follows: Correspondingly, the information of the index set after removing F 2 and F 3 can be calculated.The values are 0.67 and 0.72; therefore, the importance level of each index is 0.17, 0.17, and 0.11, respectively.Thus, the weighting of the different indexes equals 0.375, 0.375, and 0.25, respectively.Similarly, the weighting of the different indexes equals 0.5, 0.25, and 0.25, respectively, using the above steps.Then, the comprehensive objective function is obtained by Equation ( 40) and used to optimize MECS scheduling.Figure 8 shows the values of the different objective functions under the comprehensive objective function on typical days.The comparative analysis results for MECS operation were obtained under different operation modes of the ESD.Compared with the LLC mode, if an ESD was in the OEE mode, the outputs of WPP and PV increased by 149.84 MW•h and 54.98 MW•h, respectively, for a typical summer's day, and by 267.38 MW•h and 45.42 MW•h, respectively, for a typical winter's day.Likewise, F1 increased by 8400 ¥, and F2 and F3 decreased by 100 thousand ¥ and 3.24 MW•h, respectively, whereas, F1 increased by 23.8 thousand ¥, and F2 and F3 decreased by 212.9 thousand ¥ and 11.30 MW•h on two typical days.This is due to the fact that the discharge-charge power of the ESD is not limited to the rated power in the OEE mode; thus, the ESD will charge more power in the valley period and similarly, discharge more power in the peak period.

Scheduling Results
To analyze the complementarity of different power sources, MECS scheduling was optimized for a typical day.On a typical summer day, the available outputs of HS and PV are high, whereas the available output of WPP is low.The situation is completely reversed for a typical winter day.Consequently, the MECS operation was optimized under different ESD modes.Figures 9 and 10 show the results under the OEE mode and LLC mode, respectively.In Figure 9, it can be seen that when the ESD operated in the OEE mode, the load demand was mostly satisfied by the TPP and HS on a typical summer day, and by TPP and WPP on a typical winter day.It should be noted that the PV was also a major power source with a power output of 1673 MW•h from 11:00 h to 15:00 h.The WPP was another major power source with a power output of 2638 MW•h from 16:00 h to 21:00 h.Moreover, since the load demand is greater for a typical winter day than that a typical summer day, TPP2 was always scheduled to have a power output of 5000 MW•h.Then, the ESD was scheduled to overcome the intermittency and volatility of the WPP power output without considering the limitations of the rated power.The total power output of the The comparative analysis results for MECS operation were obtained under different operation modes of the ESD.Compared with the LLC mode, if an ESD was in the OEE mode, the outputs of WPP and PV increased by 149.84 MW•h and 54.98 MW•h, respectively, for a typical summer's day, and by 267.38 MW•h and 45.42 MW•h, respectively, for a typical winter's day.Likewise, F 1 increased by 8400 ¥, and F 2 and F 3 decreased by 100 thousand ¥ and 3.24 MW•h, respectively, whereas, F 1 increased by 23.8 thousand ¥, and F 2 and F 3 decreased by 212.9 thousand ¥ and 11.30 MW•h on two typical days.This is due to the fact that the discharge-charge power of the ESD is not limited to the rated power in the OEE mode; thus, the ESD will charge more power in the valley period and similarly, discharge more power in the peak period.

Scheduling Results
To analyze the complementarity of different power sources, MECS scheduling was optimized for a typical day.On a typical summer day, the available outputs of HS and PV are high, whereas the available output of WPP is low.The situation is completely reversed for a typical winter day.Consequently, the MECS operation was optimized under different ESD modes.Figures 9 and 10 show the results under the OEE mode and LLC mode, respectively.optimize MECS scheduling.Figure 8 shows the values of the different objective functions under the comprehensive objective function on typical days.The comparative analysis results for MECS operation were obtained under different operation modes of the ESD.Compared with the LLC mode, if an ESD was in the OEE mode, the outputs of WPP and PV increased by 149.84 MW•h and 54.98 MW•h, respectively, for a typical summer's day, and by 267.38 MW•h and 45.42 MW•h, respectively, for a typical winter's day.Likewise, F1 increased by 8400 ¥, and F2 and F3 decreased by 100 thousand ¥ and 3.24 MW•h, respectively, whereas, F1 increased by 23.8 thousand ¥, and F2 and F3 decreased by 212.9 thousand ¥ and 11.30 MW•h on two typical days.This is due to the fact that the discharge-charge power of the ESD is not limited to the rated power in the OEE mode; thus, the ESD will charge more power in the valley period and similarly, discharge more power in the peak period.

Scheduling Results
To analyze the complementarity of different power sources, MECS scheduling was optimized for a typical day.On a typical summer day, the available outputs of HS and PV are high, whereas the available output of WPP is low.The situation is completely reversed for a typical winter day.Consequently, the MECS operation was optimized under different ESD modes.Figures 9 and 10 show the results under the OEE mode and LLC mode, respectively.In Figure 9, it can be seen that when the ESD operated in the OEE mode, the load demand was mostly satisfied by the TPP and HS on a typical summer day, and by TPP and WPP on a typical winter day.It should be noted that the PV was also a major power source with a power output of 1673 MW•h from 11:00 h to 15:00 h.The WPP was another major power source with a power output of 2638 MW•h from 16:00 h to 21:00 h.Moreover, since the load demand is greater for a typical winter day than that a typical summer day, TPP2 was always scheduled to have a power output of 5000 MW•h.Then, the ESD was scheduled to overcome the intermittency and volatility of the WPP power output without considering the limitations of the rated power.The total power output of the In Figure 9, it can be seen that when the ESD operated in the OEE mode, the load demand was mostly satisfied by the TPP and HS on a typical summer day, and by TPP and WPP on a typical winter day.It should be noted that the PV was also a major power source with a power output of 1673 MW•h from 11:00 h to 15:00 h.The WPP was another major power source with a power output of 2638 MW•h from 16:00 h to 21:00 h.Moreover, since the load demand is greater for a typical winter day than that a typical summer day, TPP2 was always scheduled to have a power output of 5000 MW•h.Then, the ESD was scheduled to overcome the intermittency and volatility of the WPP power output without considering the limitations of the rated power.The total power output of the ESD discharge-charge was ±1853 MW•h on a typical winter day, i.e., ±848 MW•h higher than on a typical summer day.Overall, to satisfy the load demand and system operation constraints, the MECS takes advantage of the thermal-hydro-photovoltaic complementarity (or thermal-hydro-wind complementarity in certain periods) on a typical summer day and the thermal-wind complementarity on a typical winter day.
ESD discharge-charge was ±1853 MW•h on a typical winter day, i.e., ±848 MW•h higher than on a typical summer day.Overall, to satisfy the load demand and system operation constraints, the MECS takes advantage of the thermal-hydro-photovoltaic complementarity (or thermal-hydrowind complementarity in certain periods) on a typical summer day and the thermal-wind complementarity on a typical winter day.When the ESD operated in the LLC mode, the power output was limited by the rated power, i.e., 100 MW, as determined by Equations ( 29)-(31).The total discharge-charge power of the ESD was ±1100 MW•h, i.e., ±205 MW•h lower, on a typical summer day, and ±1200 MW•h, i.e., ±653 MW•h, on a typical winter day than in OEE mode.It should be noted that the available outputs of the WPP and the HS were low in the valley period; therefore, TPP4 was scheduled to satisfy the surplus load demand with a power output of 200 MW•h.The ESD was primarily discharged for economic benefit, whereas TPP2 was not scheduled in the peak period on a typical summer day.The optimal operation strategy should be chosen considering the specific situation.Figure 11 shows the complementarity indexes of MECS operation under different ESD operation modes.First, the influences of different ESD operation modes were analyzed.As the power output of the ESD in the OEE mode is higher, the outputs of the WPP and PV were also higher than in the LLC mode.Thus, when the MECS was scheduled on a typical summer day, the values of LTD and SPC were 5% and 6.78 MW•h higher, and the coal consumption was 5.86 g/kW•h lower in the OEE mode than in the LLC mode.When the MECS was scheduled on a typical winter day, the HS and TPP3 were required more to absorb WPP output more by providing a reserve service.This decreased the load tracking ability; thus, LTD was 8% lower in the OEE mode than in the LLC mode.Table 4 shows the power output results of MECS scheduling on a typical day.
Subsequently, the results of the MECS operation on different typical days were analyzed.As the available WPP output is higher on a typical winter day, more reserve service was required from the ESD, TPP3, and TPP4.The actual curves of load demand and MECS power output exhibited When the ESD operated in the LLC mode, the power output was limited by the rated power, i.e., 100 MW, as determined by Equations ( 29)-(31).The total discharge-charge power of the ESD was ±1100 MW•h, i.e., ±205 MW•h lower, on a typical summer day, and ±1200 MW•h, i.e., ±653 MW•h, on a typical winter day than in OEE mode.It should be noted that the available outputs of the WPP and the HS were low in the valley period; therefore, TPP4 was scheduled to satisfy the surplus load demand with a power output of 200 MW•h.The ESD was primarily discharged for economic benefit, whereas TPP2 was not scheduled in the peak period on a typical summer day.The optimal operation strategy should be chosen considering the specific situation.Figure 11 shows the complementarity indexes of MECS operation under different ESD operation modes.
ESD discharge-charge was ±1853 MW•h on a typical winter day, i.e., ±848 MW•h higher than on a typical summer day.Overall, to satisfy the load demand and system operation constraints, the MECS takes advantage of the thermal-hydro-photovoltaic complementarity (or thermal-hydrowind complementarity in certain periods) on a typical summer day and the thermal-wind complementarity on a typical winter day.When the ESD operated in the LLC mode, the power output was limited by the rated power, i.e., 100 MW, as determined by Equations ( 29)-(31).The total discharge-charge power of the ESD was ±1100 MW•h, i.e., ±205 MW•h lower, on a typical summer day, and ±1200 MW•h, i.e., ±653 MW•h, on a typical winter day than in OEE mode.It should be noted that the available outputs of the WPP and the HS were low in the valley period; therefore, TPP4 was scheduled to satisfy the surplus load demand with a power output of 200 MW•h.The ESD was primarily discharged for economic benefit, whereas TPP2 was not scheduled in the peak period on a typical summer day.The optimal operation strategy should be chosen considering the specific situation.Figure 11 shows the complementarity indexes of MECS operation under different ESD operation modes.First, the influences of different ESD operation modes were analyzed.As the power output of the ESD in the OEE mode is higher, the outputs of the WPP and PV were also higher than in the LLC mode.Thus, when the MECS was scheduled on a typical summer day, the values of LTD and SPC were 5% and 6.78 MW•h higher, and the coal consumption was 5.86 g/kW•h lower in the OEE mode than in the LLC mode.When the MECS was scheduled on a typical winter day, the HS and TPP3 were required more to absorb WPP output more by providing a reserve service.This decreased the load tracking ability; thus, LTD was 8% lower in the OEE mode than in the LLC mode.Table 4 shows the power output results of MECS scheduling on a typical day.
Subsequently, the results of the MECS operation on different typical days were analyzed.As the available WPP output is higher on a typical winter day, more reserve service was required from the ESD, TPP3, and TPP4.The actual curves of load demand and MECS power output exhibited First, the influences of different ESD operation modes were analyzed.As the power output of the ESD in the OEE mode is higher, the outputs of the WPP and PV were also higher than in the LLC mode.Thus, when the MECS was scheduled on a typical summer day, the values of LTD and SPC were 5% and 6.78 MW•h higher, and the coal consumption was 5.86 g/kW•h lower in the OEE mode than in the LLC mode.When the MECS was scheduled on a typical winter day, the HS and TPP3 were required more to absorb WPP output more by providing a reserve service.This decreased the load tracking ability; thus, LTD was 8% lower in the OEE mode than in the LLC mode.Table 4 shows the power output results of MECS scheduling on a typical day.
Subsequently, the results of the MECS operation on different typical days were analyzed.As the available WPP output is higher on a typical winter day, more reserve service was required from the ESD, TPP3, and TPP4.The actual curves of load demand and MECS power output exhibited strong fluctuations with high variance.Thus, on a typical winter day, LTD was 25% lower and 8% lower, SPC was 10.84 MW•h and 7.36 MW•h higher, and UCC was 2.55 g/kW•h and 0.94 g/kW•h higher than on a typical summer day in the OEE mode and LLC mode, respectively.Overall, different ESD operation modes and seasons directly influence the MECS scheduling.Figure 12 shows the power output structure of MECS scheduling under the OEE mode.According to Figure 12, the load demand is mainly supplied by the TPPs, WPP, and HS.The complementarity between WPP and HS was significant.On a typical summer day, the outputs of the HS were 5487 MW•h and 5478 MW•h higher than the outputs of WPP in the OEE/LLC modes, respectively.However, on a typical winter day, the output of HS were 6996 MW•h and 6729 MW•h lower than the outputs of WPP in the OEE/LLC modes, respectively.Moreover, complementarity also appeared in the thermal power units.For the maximum economic benefit, TPP1 and TPP2 were primarily scheduled for high capacity power output.In contrast, TPP3 and TPP4 were primarily scheduled for the low capacity reserve service.According to Figure 12, the load demand is mainly supplied by the TPPs, WPP, and HS.The complementarity between WPP and HS was significant.On a typical summer day, the outputs of the HS were 5487 MW•h and 5478 MW•h higher than the outputs of WPP in the OEE/LLC modes, respectively.However, on a typical winter day, the output of HS were 6996 MW•h and 6729 MW•h lower than the outputs of WPP in the OEE/LLC modes, respectively.Moreover, complementarity also appeared in the thermal power units.For the maximum economic benefit, TPP1 and TPP2 were primarily scheduled for high capacity power output.In contrast, TPP3 and TPP4 were primarily scheduled for the low capacity reserve service.Based on the above analysis, the output structures of different power sources in the MECS operation plan were determined.On both typical summer and winter days, the TPP, the WPP, and the HS were the major power sources.The output proportions of the WPP and HS were 14% lower and 17% higher on a typical summer day than that on a typical winter day.The total output proportions of the PV and the ESD were the same, 9%.Overall, the complementarity of TPPs, WPP, PV, and HS was significant.The TPPs and the HS were the major power sources for a typical summer day, whereas the TPPs and WPP were the major power sources for a typical winter day.Complementarity should be utilized to ensure safe and stable operation.

Results Analysis
The above analysis has already demonstrated the influences of different ESD modes on MECS operation.To prove the optimization effect of the ESD on MECS operation, the MECS was scheduled without the ESD to comparatively analyze the results of different operation scenarios.The comparative analysis index encompassed the objective function values, the pollutant emissions, the complementarity indexes, and the load demand curve.Table 5 shows the comparative results of MECS scheduling under different operation scenarios.Based on the above analysis, the output structures of different power sources in the MECS operation plan were determined.On both typical summer and winter days, the TPP, the WPP, and the HS were the major power sources.The output proportions of the WPP and HS were 14% lower and 17% higher on a typical summer day than that on a typical winter day.The total output proportions of the PV and the ESD were the same, 9%.Overall, the complementarity of TPPs, WPP, PV, and HS was significant.The TPPs and the HS were the major power sources for a typical summer day, whereas the TPPs and WPP were the major power sources for a typical winter day.Complementarity should be utilized to ensure safe and stable operation.

Results Analysis
The above analysis has already demonstrated the influences of different ESD modes on MECS operation.To prove the optimization effect of the ESD on MECS operation, the MECS was scheduled without the ESD to comparatively analyze the results of different operation scenarios.The comparative analysis index encompassed the objective function values, the pollutant emissions, the complementarity indexes, and the load demand curve.Table 5 shows the comparative results of MECS scheduling under different operation scenarios.When the MECS was scheduled without the ESD, the values of all objective functions worsened.For example, in the LLC mode, the value of F1 decreased by 43.1 thousand ¥ on a typical summer day, and by 47.7 thousand ¥ on a typical winter day, compared with the corresponding values in the OEE mode.In the OEE mode, the value of F3 increased by 3.78 MW on a typical summer day and by 1.52 MW on a typical winter day, compared with the corresponding values in the LLC mode.The outputs of WPP and PV decreased without ESD, as the load demand had to be partly satisfied by thermal power units.Therefore, pollution emissions increased significantly.For example, the total CO2 emission increased by 378 tons and by 189 tons on two typical days, compared with the values in the LLC mode.Consequently, the values of the complementarity indexes also worsened.The LTD reached its lowest value, i.e., 30%, and the SPC and UCC reached their highest values, i.e., 26.17 MW•h and 329.9 g/kW•h.Figure 13 shows the scheduling results of the MECS operation without the ESD on a typical day.According to Figure 13, when ESD was not considered, TPP4 was used to provide a reserve service on a typical summer day and part-time service on a typical winter day.TPP3 was also scheduled for full-time load demand satisfaction, and part-time reserve service.For example, TPP3 was scheduled with a power output of 199 MW•h from 3:00 h to 5:00 h as a reserve service.Moreover, the abandoned power outputs of WPP and PV were 1123.8MW•h and 412.35MW•h, respectively, on a typical summer day and 2005.35MW•h and 340.65 MW•h on a typical winter day.Thus, the total power generation cost reached 6148.47 tce and 6442.46 tce on the two typical days, respectively.Obviously, the power output structure without the ESD was worse than that with the ESD. Figure 14 shows the net output curve of the thermal power units under different operational strategies.
It can be seen that as the ESD charges when load demand is low and discharges when load demand is high, the load demand curve is smoother.Combined with Table 5, if the ESD is not considered, the peak-valley ratio of the net load demand (thermal power out) curve reached the highest value of 2.24 on a typical summer day, and 1.37 on a typical winter day.However, the corresponding ratios in the OEE mode were 1.75 and 1.26, and 1.95 and 1.31 in the LLC mode.According to Figure 13, when ESD was not considered, TPP4 was used to provide a reserve service on a typical summer day and part-time service on a typical winter day.TPP3 was also scheduled for full-time load demand satisfaction, and part-time reserve service.For example, TPP3 was scheduled with a power output of 199 MW•h from 3:00 h to 5:00 h as a reserve service.Moreover, the abandoned power outputs of WPP and PV were 1123.8MW•h and 412.35MW•h, respectively, on a typical summer day and 2005.35MW•h and 340.65 MW•h on a typical winter day.Thus, the total power generation cost reached 6148.47 tce and 6442.46 tce on the two typical days, respectively.Obviously, the power output structure without the ESD was worse than that with the ESD. Figure 14 shows the net output curve of the thermal power units under different operational strategies.
It can be seen that as the ESD charges when load demand is low and discharges when load demand is high, the load demand curve is smoother.Combined with Table 5, if the ESD is not considered, the peak-valley ratio of the net load demand (thermal power out) curve reached the highest value of 2.24 on a typical summer day, and 1.37 on a typical winter day.However, the corresponding ratios in the OEE mode were 1.75 and 1.26, and 1.95 and 1.31 in the LLC mode.Overall, to smooth the load demand curve to improve the space of WPP and PV, ESD should be integrated into the MECS.
To analyze the influence of the capacity of the ESD on MECS operation, a sensitivity analysis of the ESD capacity was carried out.The values of the objective functions and the complementarity indexes were calculated for different ESD capacities.Assuming that the maximum and rated discharge-charge power of the ESD are positively correlated with the ESD capacity, the optimal MECS operation can be obtained for different ESD capacities.Figure 15 shows the objective function value of the MECS operation with different ESD capacities.Overall, to smooth the load demand curve to improve the space of WPP and PV, ESD should be integrated into the MECS.
To analyze the influence of the capacity of the ESD on MECS operation, a sensitivity analysis of the ESD capacity was carried out.The values of the objective functions and the complementarity indexes were calculated for different ESD capacities.Assuming that the maximum and rated discharge-charge power of the ESD are positively correlated with the ESD capacity, the optimal MECS operation can be obtained for different ESD capacities.Figure 15 shows the objective function value of the MECS operation with different ESD capacities.It can be seen that when the MECS was scheduled on a typical summer day, the ESD capacity demand for MECS operation was low, as the available WPP output was relatively low.When the ESD capacity reached 1000 MW in the OEE mode, the values of objective functions reached the limit values, namely, F1, F2, and F3 were equal to 1553.8 thousand Yuan, 708.2 thousand Yuan, and 24.32 MW, respectively.The values of the objective functions were nearly unchanged when the capacity was more than 1000 MW, namely, the capacity ratios of the ESD to the WPP and the ESD to the PV reached 0.77.The ratio was 1.08 in the LLC mode.
However, when the ESD was scheduled on a typical winter day, the ESD capacity demand for the MECS operation was relatively high as the available WPP output was high.When the ESD capacity reached 1400 MW in the OEE mode, the values of the objective functions reached the same limit values as on a typical summer day.Thus, the corresponding ratios were 0.92 and 1.23.Overall, the values of the objective functions were better when the ESD capacity was high.However, they reached the limit values when the ESD capacity reached the required capacity. Figure 16 shows the values of the complementarity indexes under different ESD capacities.Overall, to smooth the load demand curve to improve the space of WPP and PV, ESD should be integrated into the MECS.
To analyze the influence of the capacity of the ESD on MECS operation, a sensitivity analysis of the ESD capacity was carried out.The values of the objective functions and the complementarity indexes were calculated for different ESD capacities.Assuming that the maximum and rated discharge-charge power of the ESD are positively correlated with the ESD capacity, the optimal MECS operation can be obtained for different ESD capacities.Figure 15 shows the objective function value of the MECS operation with different ESD capacities.It can be seen that when the MECS was scheduled on a typical summer day, the ESD capacity demand for MECS operation was low, as the available WPP output was relatively low.When the ESD capacity reached 1000 MW in the OEE mode, the values of objective functions reached the limit values, namely, F1, F2, and F3 were equal to 1553.8 thousand Yuan, 708.2 thousand Yuan, and 24.32 MW, respectively.The values of the objective functions were nearly unchanged when the capacity was more than 1000 MW, namely, the capacity ratios of the ESD to the WPP and the ESD to the PV reached 0.77.The ratio was 1.08 in the LLC mode.
However, when the ESD was scheduled on a typical winter day, the ESD capacity demand for the MECS operation was relatively high as the available WPP output was high.When the ESD capacity reached 1400 MW in the OEE mode, the values of the objective functions reached the same limit values as on a typical summer day.Thus, the corresponding ratios were 0.92 and 1.23.Overall, the values of the objective functions were better when the ESD capacity was high.However, they reached the limit values when the ESD capacity reached the required capacity. Figure 16 shows the values of the complementarity indexes under different ESD capacities.It can be seen that when the MECS was scheduled on a typical summer day, the ESD capacity demand for MECS operation was low, as the available WPP output was relatively low.When the ESD capacity reached 1000 MW in the OEE mode, the values of objective functions reached the limit values, namely, F 1 , F 2 , and F 3 were equal to 1553.8 thousand Yuan, 708.2 thousand Yuan, and 24.32 MW, respectively.The values of the objective functions were nearly unchanged when the capacity was more than 1000 MW, namely, the capacity ratios of the ESD to the WPP and the ESD to the PV reached 0.77.The ratio was 1.08 in the LLC mode.
However, when the ESD was scheduled on a typical winter day, the ESD capacity demand for the MECS operation was relatively high as the available WPP output was high.When the ESD capacity reached 1400 MW in the OEE mode, the values of the objective functions reached the same limit values as on a typical summer day.Thus, the corresponding ratios were 0.92 and 1.23.Overall, the values of the objective functions were better when the ESD capacity was high.However, they reached the limit values when the ESD capacity reached the required capacity. Figure 16 shows the values of the complementarity indexes under different ESD capacities.It can be seen that the values of the complementarity indexes improved as the ESD capacity increased.It should be noted that this trend was not always continuous.When the MECS was scheduled on a typical summer day, the change was small once the ESD capacity reached 800 MW in the OEE mode and 1300 MW in the LLC mode, i.e., the capacity ratios of the ESD to the WPP and the ESD to the PV were 0.62 and 1, respectively.However, when the MECS was scheduled on a typical winter day, the change was small once the ESD capacity reached 1300 MW in the OEE mode and 1700 MW in the LLC mode, i.e., the capacity ratios of the ESD to the WPP and the ESD to the PV were 1 and 1.31, respectively.Overall, the change trend of the complementarity indexes was nearly identical to that of the objective functions.When the capacity ratios of the ESD to the WPP and the ESD to the PV belonged to [0.62, 0.77] in the OEE mode and [1, 1.08] in the LLC mode, the optimal results of MECS scheduling were obtained on a typical summer day, whereas for capacity ratios in [0.92, 1] in the OEE mode and [1.23, 1.31] in the LLC mode, the optimal results were obtained on a typical winter day.Therefore, the optimal ESD capacity and operation strategy for MECS operation should be determined according to the actual operating environment.

Conclusions
WPPs, PVs, HSs, TPPs, and ESDs were integrated into a MECS.This paper proposed a multi-objective model for MECS operation considering different operation constraint conditions under three objective functions: the maximum operation revenue (F1), the minimum abandoned energy cost (F2), and the minimum output fluctuation (F3).Then, the rough set theory was applied to handle the information of the payoff table of the objective functions to calculate the weight coefficients.Finally, the improved IEEE 14-bus system was used as the simulation system.The following conclusions were gained: (1) The proposed uncertainty analysis and simulation method generated the most appropriate scenarios for wind power, photovoltaic power, and other uncertainty factors.The Wasserstein method achieved a better discretization by dividing the CPDF into several unequal subintervals.The scenario tree method comprehensively simulated the uncertainty factors by considering different occurrence states, whereas the improved K-medoids method selected the most representative scenarios more efficiently by calculating the relative K-distances of different scenarios and ensuring the updated clustering center.
(2) The MECS fully utilized the complementarity of the different power sources, and the proposed multi-objective model obtained the optimal scheme for MECS scheduling.TPPs, HSs, and PVs were the major power sources used to satisfy the load demand on a typical summer day, whereas TPPs and WPPs were the major power sources used for load demand on a typical winter day.Compared with the results of the single optimization objective, the scheduling results better balanced the revenue and risk when multiple optimization objectives were taken into consideration.
(3) The operation results of MECS scheduling improved as the ESD capacity increased, and they were restricted when the capacity reached the required value.The change trend of the complementarity indexes was nearly identical to that of the objective functions.The ESD optimized the MECS scheduling by utilizing the characteristics of power discharge-charge.The ESD operation mode directly influenced the optimality of the results.In the OEE mode, the discharge-charge power was better utilized than in the LLC mode; moreover, the LLC mode ensured the ESD's It can be seen that the values of the complementarity indexes improved as the ESD capacity increased.It should be noted that this trend was not always continuous.When the MECS was scheduled on a typical summer day, the change was small once the ESD capacity reached 800 MW in the OEE mode and 1300 MW in the LLC mode, i.e., the capacity ratios of the ESD to the WPP and the ESD to the PV were 0.62 and 1, respectively.However, when the MECS was scheduled on a typical winter day, the change was small once the ESD capacity reached 1300 MW in the OEE mode and 1700 MW in the LLC mode, i.e., the capacity ratios of the ESD to the WPP and the ESD to the PV were 1 and 1.31, respectively.Overall, the change trend of the complementarity indexes was nearly identical to that of the objective functions.When the capacity ratios of the ESD to the WPP and the ESD to the PV belonged to [0.62, 0.77] in the OEE mode and [1, 1.08] in the LLC mode, the optimal results of MECS scheduling were obtained on a typical summer day, whereas for capacity ratios in [0.92, 1] in the OEE mode and [1.23, 1.31] in the LLC mode, the optimal results were obtained on a typical winter day.Therefore, the optimal ESD capacity and operation strategy for MECS operation should be determined according to the actual operating environment.

Conclusions
WPPs, PVs, HSs, TPPs, and ESDs were integrated into a MECS.This paper proposed a multi-objective model for MECS operation considering different operation constraint conditions under three objective functions: the maximum operation revenue (F 1 ), the minimum abandoned energy cost (F 2 ), and the minimum output fluctuation (F 3 ).Then, the rough set theory was applied to handle the information of the payoff table of the objective functions to calculate the weight coefficients.Finally, the improved IEEE 14-bus system was used as the simulation system.The following conclusions were gained: (1) The proposed uncertainty analysis and simulation method generated the most appropriate scenarios for wind power, photovoltaic power, and other uncertainty factors.The Wasserstein method achieved a better discretization by dividing the CPDF into several unequal subintervals.The scenario tree method comprehensively simulated the uncertainty factors by considering different occurrence states, whereas the improved K-medoids method selected the most representative scenarios more efficiently by calculating the relative K-distances of different scenarios and ensuring the updated clustering center.
(2) The MECS fully utilized the complementarity of the different power sources, and the proposed multi-objective model obtained the optimal scheme for MECS scheduling.TPPs, HSs, and PVs were the major power sources used to satisfy the load demand on a typical summer day, whereas TPPs and WPPs were the major power sources used for load demand on a typical winter day.Compared with the results of the single optimization objective, the scheduling results better balanced the revenue and risk when multiple optimization objectives were taken into consideration.
(3) The operation results of MECS scheduling improved as the ESD capacity increased, and they were restricted when the capacity reached the required value.The change trend of the complementarity indexes was nearly identical to that of the objective functions.The ESD optimized the MECS scheduling by utilizing the characteristics of power discharge-charge.The ESD operation mode directly influenced the optimality of the results.In the OEE mode, the discharge-charge power was better utilized than in the LLC mode; moreover, the LLC mode ensured the ESD's operation life cycle.

Figure 1 .
Figure 1.The basic structure of the multi-energy complementary system (MECS) operation.

Figure 1 .
Figure 1.The basic structure of the multi-energy complementary system (MECS) operation.

Figure 2 .
Figure 2. The discretization of the continuous probability distribution functions (CPDFs) (six discrete points).

Figure 2 .
Figure 2. The discretization of the continuous probability distribution functions (CPDFs) (six discrete points).

Figure 4 .
Figure 4.The load demand and the available output of the hydropower station (HS) on a typical load day.

Figure 4 .
Figure 4.The load demand and the available output of the hydropower station (HS) on a typical load day.

Figure 5 .
Figure 5.The ten simulation scenarios of WPPs and PVs on a typical load day (summer).

Figure 6 .
Figure 6.The average outputs of WPP and PV on a typical load day (summer and winter).

Figure 7 .
Figure 7.The optimization flowchart for the MECS scheduling model.

Figure 5 .
Figure 5.The ten simulation scenarios of WPPs and PVs on a typical load day (summer).

25 Figure 5 .
Figure 5.The ten simulation scenarios of WPPs and PVs on a typical load day (summer).

Figure 6 .
Figure 6.The average outputs of WPP and PV on a typical load day (summer and winter).

Figure 7 .
Figure 7.The optimization flowchart for the MECS scheduling model.

Figure 6 .
Figure 6.The average outputs of WPP and PV on a typical load day (summer and winter).

25 Figure 5 .
Figure 5.The ten simulation scenarios of WPPs and PVs on a typical load day (summer).

Figure 6 .
Figure 6.The average outputs of WPP and PV on a typical load day (summer and winter).

Figure 7 .
Figure 7.The optimization flowchart for the MECS scheduling model.Figure 7. The optimization flowchart for the MECS scheduling model.

Figure 7 .
Figure 7.The optimization flowchart for the MECS scheduling model.Figure 7. The optimization flowchart for the MECS scheduling model.

Figure 8 .
Figure 8.The different objective values under the comprehensive objective function on typical days.

Figure 9 .
Figure 9.The MECS scheduling results on typical days under the OEE mode: summer (left) and winter (right).

Figure 8 .
Figure 8.The different objective values under the comprehensive objective function on typical days.

Figure 8 .
Figure 8.The different objective values under the comprehensive objective function on typical days.

Figure 9 .
Figure 9.The MECS scheduling results on typical days under the OEE mode: summer (left) and winter (right).

Figure 9 .
Figure 9.The MECS scheduling results on typical days under the OEE mode: summer (left) and winter (right).

Figure 10 .
Figure 10.The MECS scheduling results on a typical day under the LLC mode: summer (left) and winter (right).

Figure 11 .
Figure 11.The complementarity indexes of MECS operation under different energy storage devices (ESD) operation modes.

Figure 10 .
Figure 10.The MECS scheduling results on a typical day under the LLC mode: summer (left) and winter (right).

Figure 10 .
Figure 10.The MECS scheduling results on a typical day under the LLC mode: summer (left) and winter (right).

Figure 11 .
Figure 11.The complementarity indexes of MECS operation under different energy storage devices (ESD) operation modes.

Figure 11 .
Figure 11.The complementarity indexes of MECS operation under different energy storage devices (ESD) operation modes.

Figure 12 .
Figure 12.The power output structure of MECS scheduling under the OEE mode: (a) summer and (b) winter.

Figure 12 .
Figure 12.The power output structure of MECS scheduling under the OEE mode: (a) summer and (b) winter.

Figure 13 .
Figure 13.The scheduling results of the MECS operation without ESD on a typical day: (a) summer and (b) winter.

Figure 13 .
Figure 13.The scheduling results of the MECS operation without ESD on a typical day: (a) summer and (b) winter.

Figure 14 .
Figure 14.The net output curve of the thermal power units under different operational strategies.

Figure 15 .
Figure 15.The objective function values of MECS operation under different ESD capacities.

Figure 14 .
Figure 14.The net output curve of the thermal power units under different operational strategies.

Figure 14 .
Figure 14.The net output curve of the thermal power units under different operational strategies.

Figure 15 .
Figure 15.The objective function values of MECS operation under different ESD capacities.

Figure 15 .
Figure 15.The objective function values of MECS operation under different ESD capacities.

Figure 16 .
Figure 16.The complementarity indexes of MECS operation under different ESD capacities.

Figure 16 .
Figure 16.The complementarity indexes of MECS operation under different ESD capacities.

Table 1 .
The payoff table of the proposed multi-objective model.
IIAccording to Table

Table 2 .
The key parameters of the thermal power plant (TPP) operation.

Table 2 .
The key parameters of the thermal power plant (TPP) operation.

Table 3 .
The payoff table of the MECS operation in Case 1.

Table 4 .
The power output results of MECS scheduling on a typical day.ESD: energy storage device; HS: hydropower station.
Appl.Sci.2018,8,xFORPEER REVIEW 20 of 25 strong fluctuations with high variance.Thus, on a typical winter day, LTD was 25% lower and 8% lower, SPC was 10.84 MW•h and 7.36 MW•h higher, and UCC was 2.55 g/kW•h and 0.94 g/kW•h higher than on a typical summer day in the OEE mode and LLC mode, respectively.Overall, different ESD operation modes and seasons directly influence the MECS scheduling.Figure12shows the power output structure of MECS scheduling under the OEE mode.

Table 4 .
The power output results of MECS scheduling on a typical day.ESD: energy storage device; HS: hydropower station.

Table 5 .
The comparative results of MECS scheduling under different operation scenarios.