Sample Entropy Based Net Load Tracing Dispatch of New Energy Power System

: The high-proportion of renewable energies is gradually becoming one of the main power supply sources and bringing strong uncertainties to the power grid. In this paper, a sample entropy (SampEn) based net load tracing dispatch strategy with a speciﬁc thermal generating mode is proposed. In this strategy, renewable energies are fully and preferentially consumed by electric loads, turned to net loads, to maximize the utilization of renewable energies. SampEn theory is utilized to evaluate the complexity of net load time series, based on which, the traditional power generators trace the complexity of the net load ﬂexibly. According to the SampEn, a speciﬁc generating model of thermal generators is determined and the cooperation between thermal generators and pumped storage is realized, aiming at reducing the ramp power of thermal generators and increasing the throughput of pumped storage. The experiment simulation is developed on the 10-unit test system. Results show that the ramping power of the thermal generators are reduced 43% and 13% in the two cases together with the throughput of pumped storage is increased 44% and 27% on the premise that the economy of the system is maintained and renewable energies are fully consumed. Therefore, the efﬁciency and reasonability of the proposed dispatch strategy are conﬁrmed.


Introduction
The modern power system is transforming to the third-generation power system, or called new energy power system, which mainly depends on those clean and renewable non-fossil energies [1], such as wind power [2], solar energy [3], and so forth. In the new energy power system, the intermittent energies such as wind power and photovoltaic share the responsibility of power load demand together with the traditional power sources. The thermal power output in the traditional day-ahead power dispatch with high uncertainties will ramp up and down frequently. Thus, the inlet valve of steam turbine will be adjusted accordingly, which will lead to the difficulties on actual operation and the aggravation of mechanical wear [4,5]. In the long run, the power generation efficiency will be decreased and the life of the generators will be shortened. Therefore, the reasonable determination of unit generating mode, which can not only respond to these strong uncertainties, but reduce the power output fluctuation and increase the utilization efficiency of non-fossil energies as well, plays an important role in the new power system dispatch.
The renewable energies with high proportion have gradually been as important as traditional power sources and the proportion is expected to reach up to 60% by 2050 [6]. Moreover, the distributed generators, electric vehicles, energy storages and so forth on the load side increase the uncertainties of the system as well. Hence, researchers introduced the concept of net load [7,8]. The net load means of the system in the premise that the economy of the system is ensured. In order to make full use of the renewable energy, wind power and photovoltaic power are consumed by electric loads first, turned to net loads. Focusing on the day-ahead power dispatch, the SampEn theory is used to evaluate the complexity level of net loads. On these bases, the traditional dispatch is transformed to a specific generating mode according to the SampEn, in which the thermal ramping mode and time duration are determined. The high operation efficiency and less ramping power of thermal generators are achieved. Moreover, the high throughput of pumped storage is realized. A simulation model based on a 10-unit power system is developed and tested in MATLAB (R2018a, MathWorks, Natick, MA, USA). This paper is organized as follows. Section 2 analyzed the characteristics of the net loads and introduced the SampEn calculation of net loads. In Section 3, a net load tracing dispatch strategy based on SampEn is studied in detail. The case study and discussion are conducted in Section 4. Finally, the conclusions are stated in Section 5.

Net Loads Description
The variation of load data and wind power in a specific province from 3 March to 24 April 2017 are shown in Figures 1 and 2, respectively. Based on those data, the electric loads are 12.5% scale-down to match the wind data for simulating a power grid operation condition with high-percentage wind power. After the scale-down, the trend of curve remains unchanged and the differences between wind power and electric loads are defined as net loads shown in Figure 3. As shown in Figure 3, the net loads are highly stochastic. No obvious regulation or trend can be found in these figures. If the power dispatch strategies remain the traditional ones, such highly stochastic net loads will result in large fluctuations on the power supply side, causing thermal generation power outputs to change frequently and decrease the efficiency and stability of the system. This paper focuses on the dispatch on the net loads with strong uncertainty and large fluctuation. The time series is used to describe the net loads and the SampEn is utilized to evaluate the complexity level of the net loads.
Energies 2018, 11, x FOR PEER REVIEW 3 of 23 stability of the system in the premise that the economy of the system is ensured. In order to make full use of the renewable energy, wind power and photovoltaic power are consumed by electric loads first, turned to net loads. Focusing on the day-ahead power dispatch, the SampEn theory is used to evaluate the complexity level of net loads. On these bases, the traditional dispatch is transformed to a specific generating mode according to the SampEn, in which the thermal ramping mode and time duration are determined. The high operation efficiency and less ramping power of thermal generators are achieved. Moreover, the high throughput of pumped storage is realized. A simulation model based on a 10-unit power system is developed and tested in MATLAB (R2018a, MathWorks, Natick, MA, USA). This paper is organized as follows. Section 2 analyzed the characteristics of the net loads and introduced the SampEn calculation of net loads. In Section 3, a net load tracing dispatch strategy based on SampEn is studied in detail. The case study and discussion are conducted in Section 4. Finally, the conclusions are stated in Section 5.

Net Loads Description
The variation of load data and wind power in a specific province from 3rd March to 24th April, 2017 are shown in Figures 1 and 2, respectively. Based on those data, the electric loads are 12.5% scale-down to match the wind data for simulating a power grid operation condition with high-percentage wind power. After the scale-down, the trend of curve remains unchanged and the differences between wind power and electric loads are defined as net loads shown in Figure 3. As shown in Figure 3, the net loads are highly stochastic. No obvious regulation or trend can be found in these figures. If the power dispatch strategies remain the traditional ones, such highly stochastic net loads will result in large fluctuations on the power supply side, causing thermal generation power outputs to change frequently and decrease the efficiency and stability of the system.   This paper focuses on the dispatch on the net loads with strong uncertainty and large fluctuation. The time series is used to describe the net loads and the SampEn is utilized to evaluate the complexity level of the net loads.

SampEn of Net Loads in New Energy Power System
According to the theory of SampEn, the value of SampEn is positive correlation to the complexity level of sample series. That is, the time series with larger SampEn value are more complicated than the one with lower SampEn value. The calculation of SampEn does not include the comparison of its own data series. Therefore, the calculation is independent to the length of the series. The comparison of two time series will still be the same in any of their subsequence dimension and with similar tolerance level. In addition, the SampEn calculation can remain unaffected with over one-third data loss [43].
In the theory of SampEn, the complexity means the generation rate of new data mode with time. In a time series with larger SampEn, new data modes will be generated with higher frequently and vice versa [44]. By defining the sign change of the slope as the new mode, the generation rate of new mode is reflected through the polarity changing frequency of the slope. The change of slope sign indicates the direction of the net load variation has been altered. The changing frequency of the slope sign reflects not only the changing speed of value but also the changing speed of changing direction. On the contrary, if the slope signs are not changed in a time period, the net loads will be changed monotonously or limited changes in direction can be found. In this paper, the slope of the  This paper focuses on the dispatch on the net loads with strong uncertainty and large fluctuation. The time series is used to describe the net loads and the SampEn is utilized to evaluate the complexity level of the net loads.

SampEn of Net Loads in New Energy Power System
According to the theory of SampEn, the value of SampEn is positive correlation to the complexity level of sample series. That is, the time series with larger SampEn value are more complicated than the one with lower SampEn value. The calculation of SampEn does not include the comparison of its own data series. Therefore, the calculation is independent to the length of the series. The comparison of two time series will still be the same in any of their subsequence dimension and with similar tolerance level. In addition, the SampEn calculation can remain unaffected with over one-third data loss [43].
In the theory of SampEn, the complexity means the generation rate of new data mode with time. In a time series with larger SampEn, new data modes will be generated with higher frequently and vice versa [44]. By defining the sign change of the slope as the new mode, the generation rate of new mode is reflected through the polarity changing frequency of the slope. The change of slope sign indicates the direction of the net load variation has been altered. The changing frequency of the slope sign reflects not only the changing speed of value but also the changing speed of changing direction. On the contrary, if the slope signs are not changed in a time period, the net loads will be changed monotonously or limited changes in direction can be found. In this paper, the slope of the

SampEn of Net Loads in New Energy Power System
According to the theory of SampEn, the value of SampEn is positive correlation to the complexity level of sample series. That is, the time series with larger SampEn value are more complicated than the one with lower SampEn value. The calculation of SampEn does not include the comparison of its own data series. Therefore, the calculation is independent to the length of the series. The comparison of two time series will still be the same in any of their subsequence dimension and with similar tolerance level. In addition, the SampEn calculation can remain unaffected with over one-third data loss [43].
In the theory of SampEn, the complexity means the generation rate of new data mode with time. In a time series with larger SampEn, new data modes will be generated with higher frequently and vice versa [44]. By defining the sign change of the slope as the new mode, the generation rate of new mode is reflected through the polarity changing frequency of the slope. The change of slope sign indicates the direction of the net load variation has been altered. The changing frequency of the slope sign reflects not only the changing speed of value but also the changing speed of changing direction. On the contrary, if the slope signs are not changed in a time period, the net loads will be changed monotonously or limited changes in direction can be found. In this paper, the slope of the adjacent points, the changing rate of the slope sign and the valley-to-peak of the time series are calculated and analyzed. Through the rolling calculation of the SampEn and its auxiliary value of the time series, the corresponding scheduling period is divided and the strategy of the dispatch are confirmed. Based on the description above, the computational procedure of the SampEn NetLoad of the net load time series is shown as follows: 1. The whole net load time series is given: {NL l } = {nl(1), nl(2), . . . , nl(N)}, where the N is the total amount of states; 2. The slope of the adjacent points in net load {NL l } is shown in Equation (1).
where, l = 1, 2, . . . , N − 1; 3. The slope signs are calculated as Equation (2) shows and the changing markers are determined as Equation (3) shows. symbol are divided in the same time period, and the piecewise time points t g , g ∈ N are achieved. Thus, the subsequence of net load is expressed as {nl(g), nl(g + 1), . . . , nl(g + T NL )}, where T NL is the length of the net load subsequence. 6. The time series P NetLoad i = {p NL (1), p NL (2), . . . , p NL (T NL )} is assumed as the subsequence of the net load at time t g . A subsequence with dimension of m NL is generated as Equation (6) shows where, i = 1, 2, . . . , T NL − m NL + 1 and m NL is usually equal to 2 or 3. 7. The D m NL NLij (P NL (i), P NL (j)) is defined as the distance of the subsequence P NL (i) and P NL (j), where j = 1, 2, . . . , T NL − m NL + 1, j = i. The distance means the maximum value of the difference between the corresponding elements in these two subsequences. To each i the distance between P NL (i) and P NL (j) is shown in Equation (7).
Energies 2019, 12,193 6 of 23 Thus, the distance of the subsequence P NL (i) and P NL (j) is calculated as Equation (8) shows.
8. The tolerance level r NL is a threshold defining when two states are similar, where r NL > 0. In this paper, the r NL was set 0.2 × STD NL , STD NL is the standard deviation of NL. To each i, the amount of D m NL NLij (P NL (i), P NL (j)) < r NL are calculated. B m NL NLi (r NL ) is defined as Equation (9): where T NL − m NL is the total distance of the sample. 9. The actual average over all the vectors B m NL NLi (r NL ) of is computed as Equation (10) shows.
10. The dimension of subsequence is increased to m + 1. Repeat steps 6 to 10 and achieve B m NL +1 Avg (r NL ). Theoretically, the sample entropy is defined as In practice, the length of time series T NL is finite. Therefore, the estimation of SampEn is denoted as Equation (12).

SampEn Application of Net Loads
The net load data of a province on 3 March and 4 April are shown in Table 1. The net load curve and the adjacent point slope of 11 April and 3 March are shown in Figures 4 and 5, respectively. According to the calculation process of SampEn in Section 2.2, the subsequence series are divided into different time length and calculated correspondingly. As Table 1 shows, the ratio of slope sign changing frequency to valley-to-peak is proportional to the SampEn, which also can be verified by Figures 4 and 5. The data in time interval of 1 to 18 and 42 to 69 in Figure 4 show that the slope sign of adjacent points change frequently while the ratio of slope sign changing amount to valley-to-peak stays large. Thus, the curve fluctuates irregularly in a small range and the corresponding SampEn is large. The above result analysis is also suitable for the data in time interval 46 to 96 in Figure 5. In addition, the data in time interval of 19 to 41 and 70 to 96 in Figure 4 show that the slope sign of adjacent points changing amount are less and the ratio of slope sign changing amount to valley-to-peak is small. Thus, the curve appears to monotonic variation in a certain time period and only a single or small numbers of sharply data changes and changing tendency can be found. The complexity level of the subsequence is low and the corresponding SampEn is small. The above result analysis is also suitable for the data in time intervals 1 to 45 in Figure 5.

Generating Mode of Thermal Generators
The complexity level of the net loads can be reflected by SampEn. The net load with larger SampEn will be at a higher complexity level. In a certain time periods, the data of net loads change frequently and there is no obviously data increase or decrease. Thus, the power outputs of thermal

Generating Mode of Thermal Generators
The complexity level of the net loads can be reflected by SampEn. The net load with larger SampEn will be at a higher complexity level. In a certain time periods, the data of net loads change frequently and there is no obviously data increase or decrease. Thus, the power outputs of thermal generators prefer to remain constant and the pumped storage is as the main supply power source. On the contrary, the net load with lower SampEn will be at a lower complexity level. In certain time periods, there are obviously increase and decrease of the net loads. In addition, the variations will last periods of time. Although the pumped storage is available to respond to the increase and decrease of the net loads, the power outputs cannot last for a long time. Thus, the thermal generators and the pumped storage are in cooperation as the main supply power sources.
The generating mode of thermal generators based on SampEn is determined as follows: The SampEn SampEn NetLoad t g t l of net load series P NLt g t l in time periods t g to t l is calculated as Section 2.2 shows. The generating mode setting is depended on the value of SampEn NetLoad If SampEn NetLoad t g t l is large, the complexity level in this time period is high and the fluctuation of thermal generation are minimized. If SampEn NetLoad t g t l is small, the complexity level in this time period is low and the thermal generators are operated normally to response to the net loads. The generating mode is settled as Equation (13) shows.
where P pi,t is the real power output of unit pi at time t. P min pi and P max pi are the minimum and maximum power outputs of the ith thermal generator, respectively.

Objective Functions
The thermal generation cost F thermal is considered as the objective function in this paper, which is calculated by where T is total number of time periods; ng is the number of dispatchable units; a pi , b pi , and c pi are the coefficients of fuel cost functions for units pi.

Power Balance Equations
The power balance equations at time t are formulated as where ng ∑ pi=1 P pi,t is the total generator power output during the tth time period; P Lt is the total load demand at time t. P Hgeni,t is the generating power of pumped storage Hgeni at time t and P H pumpi,t is the pumping power of pumped storage H pumpi at time t. P NLt is the total net load demands at time t. P pre wpi,t is the predicted wind power of wind farm wpi at time t and P pre solari,t is the predicted photovoltaic power of photovoltaic station solari at time t. np, nw, and ns are the total number of the pumped storage, wind farms, and photovoltaic power stations, respectively.

2.
The constraints of pumped storage The pump and generating power constraints of pumped storage are shown in Equations (17) and (18).
where P min Hgeni and P max Hgeni are the minimum and maximum generating power of pumped storage Hgeni, respectively. P min H pumpi and P max H pumpi are the minimum and maximum pumping power of pumped storage H pumpi, respectively. The constraints of the reservoir energy conversion are shown as follows: where H Pumpi,t+1 and H Pumpi,t are the storage condition of the reservoir at time t + 1 and time t. ∆t is the time interval. γ Pump is the pumping efficiency and γ Gen is the generating efficiency. H min Pumpi and H max Pumpi the minimum and maximum power storage of reservoir. Moreover, the pumping and generating processes of the pumped storage cannot be carried out at the same time, which is as Equation (21) shows. P Hgeni,t ·P H pumpi,t = 0 3.

Constraints of Thermal Generators
The power output constraints of thermal generators are shown in Equation (22).
The generating ramp rate limits are formulated by where UR pi and DR pi are the ramp-up and ramp-down rate limits of the ith thermal generator, respectively. The spinning reserve chance constraints are formulated by Equations (24) and (25). The up and down reserve should cover the errors between the actual wind power output and the predicted wind power output. In order to avoid the waste of reserve sources, the reserve needs not to cover the whole error range and only to satisfy a certain probability. In addition, the upper limits of up and down reserve are shown in Equations (26) and (27).
Pr  (27) where ρ is the confidence coefficient and P wpi,t is the actual wind power of wind farm wpi at time t. P URi,t and P DRi,t are the actual up and down reserve of thermal generator pi.

Stochastic Variables
The probability distribution of wind power is estimated based on the Beta distribution [45] and shown in Figure 6, which is described by where p wp is a normalized wind power output, B(α, β) is the Beta distribution function and α, β are the distribution shape parameters. The calculated formulas are where P wp is the actual wind power output; P min wp and P max wp are the minimum and maximum outputs of wind turbines, respectively. P Pr  Figure 6. Beta distribution curve.

Power Dispatch Strategy Process Based On SampEn
The power dispatch strategy is mainly divided in five parts: 1. Renewable energy is connected to the power grid and consumed by electric loads firstly. Thus, the net load time series is generated.
2. The characteristics and numeric features of the net loads containing the slope of the adjacent points shown in Equation (1), the total number of slope sign changes, ratio of sign changing amount to net load amount shown in Equation (4) and the ratio of proportion to the valley-to-peak of the The fractile is used to solve the above chance constraints [46]. where Y is the random variable and Y ρ is the fractile of ρ. In Figure 6, the Y ρ is the upper fractile and the Y 1−ρ is the lower fractile. f (Y) is the density function of Y. Equations (24) and (25) can be transformed into Equations (34) and (35), corresponding to the type of Equation (33). When the ρ is determined, the Y ρ and Y 1−ρ can be received by calculating the inverse function through MATLAB. Thus, when the function (36) is satisfied, the constraint (34) is ensured. When the function (37) is satisfied, the constraint (35) is ensured.

Power Dispatch Strategy Process Based On SampEn
The power dispatch strategy is mainly divided in five parts: 1. Renewable energy is connected to the power grid and consumed by electric loads firstly. Thus, the net load time series is generated.
2. The characteristics and numeric features of the net loads containing the slope of the adjacent points shown in Equation (1), the total number of slope sign changes, ratio of sign changing amount to net load amount shown in Equation (4) and the ratio of proportion to the valley-to-peak of the net loads shown in Equation (5) are analyzed.
3. The characteristics of net load mentioned in Step 2 are rolling calculated. According to the results of SampEn NetLoad t g t l , the closer results are divided into one subsequence. Thus, the net load time series are divided into a few certain subsequences.
4. The SampEn of the subsequences are calculated and the time frame is determined according to the point-in-time of the subsequences. Moreover, the generating mode of thermal generators is confirmed according to the SampEn as Equation (13) shows.
5. The power dispatch strategy based on SampEn is conducted. A prime-dual interior point method is used to solve the optimization problem.
The flow chart of the power dispatch strategy is shown in Figure 7.

Case Study
In order to demonstrate the effectiveness of the proposed power dispatch strategy, a 10- Coe  are the up and down ramping cost coefficients, respectively. The load demand changing curve is shown in Figure 8, and the maximum load is 2220 MW. The installed capacity of pumped storage is 700 MW and the total

Case Study
In order to demonstrate the effectiveness of the proposed power dispatch strategy, a 10-unit test system [47] is developed in this work. The unit data and cost coefficients are modified from [47,48] and listed in Table 2, where P max and P min are the maximum and minimum power limit of the units, respectively; a, b, and c are the coefficients of fuel cost functions for units; UR and DR are the ramp-up and ramp-down rate limits of the units, respectively; Coe URi and Coe DRi are the up and down reserve cost coefficients, respectively; Coe ramp−upi and Coe ramp−downi are the up and down ramping cost coefficients, respectively. The load demand changing curve is shown in Figure 8, and Energies 2019, 12,193 13 of 23 the maximum load is 2220 MW. The installed capacity of pumped storage is 700 MW and the total available operational thermal power is 2358 MW. There are two wind farms with 500 MW installed capacity, respectively and the wind power penetration is 45%. The predicted wind power curves and the limit curves are shown in Figure 9. The installed capacity of photovoltaic station is 600 MW and the penetration is 27% as shown in Figure 10. The percentage of renewable energy in this model is 72% and the net load curve is shown in Figure 11. available operational thermal power is 2358 MW. There are two wind farms with 500 MW installed capacity, respectively and the wind power penetration is 45%. The predicted wind power curves and the limit curves are shown in Figure 9. The installed capacity of photovoltaic station is 600 MW and the penetration is 27% as shown in Figure 10. The percentage of renewable energy in this model is 72% and the net load curve is shown in Figure 11.    available operational thermal power is 2358 MW. There are two wind farms with 500 MW installed capacity, respectively and the wind power penetration is 45%. The predicted wind power curves and the limit curves are shown in Figure 9. The installed capacity of photovoltaic station is 600 MW and the penetration is 27% as shown in Figure 10. The percentage of renewable energy in this model is 72% and the net load curve is shown in Figure 11.

SampEn Calculation
The characteristics and numeric features of the net loads are analyzed and the SampEn are calculated. The results are shown in Table 3.

Result Comparison and Analysis of Cases
The simulation in this paper is conducted in four cases.

SampEn Calculation
The characteristics and numeric features of the net loads are analyzed and the SampEn are calculated. The results are shown in Table 3.

Result Comparison and Analysis of Cases
The simulation in this paper is conducted in four cases.

SampEn Calculation
The characteristics and numeric features of the net loads are analyzed and the SampEn are calculated. The results are shown in Table 3.

Result Comparison and Analysis of Cases
The simulation in this paper is conducted in four cases.
• Case 1: power dispatch without SampEn and the wind power reserve confidence degree is 0.9. • Case 2: power dispatch based on SampEn at wind power reserve confidence degree of 0.9. • Case 3: power dispatch without SampEn and the wind power reserve confidence degree is 0.95. • Case 4: power dispatch based on SampEn at wind power reserve confidence degree of 0.95.

Results in Case 1
The power outputs of thermal generators in Case 1 are shown in Figure 12. The power outputs of pumped storage are shown in Figure 13, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power. The power outputs of thermal generators in Case 1 are shown in Figure 12. The power outputs of pumped storage are shown in Figure 13, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power.

Results in Case 2
The power outputs of thermal generators based on the SampEn in Case 2 are shown in Figure  14. The power outputs of pumped storage are shown in Figure 15, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power. The power output comparisons of pumped storage in Case 1 and Case 2 are shown in Figure 16. The power outputs of thermal generators in Case 1 are shown in Figure 12. The power outputs of pumped storage are shown in Figure 13, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power.

Results in Case 2
The power outputs of thermal generators based on the SampEn in Case 2 are shown in Figure  14. The power outputs of pumped storage are shown in Figure 15, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power. The power output comparisons of pumped storage in Case 1 and Case 2 are shown in Figure 16.

Results in Case 2
The power outputs of thermal generators based on the SampEn in Case 2 are shown in Figure 14. The power outputs of pumped storage are shown in Figure 15, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power. The power output comparisons of pumped storage in Case 1 and Case 2 are shown in Figure 16.     The power outputs of thermal generators in Case 3 are shown in Figure 17. The power outputs of pumped storage are shown in Figure 18, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating  The power outputs of thermal generators in Case 3 are shown in Figure 17. The power outputs of pumped storage are shown in Figure 18, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power.

Results in Case 3
The power outputs of thermal generators in Case 3 are shown in Figure 17. The power outputs of pumped storage are shown in Figure 18, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power.

Results in Case 4
The power outputs of thermal generators based on the SampEn in Case 4 are shown in Figure  19. The power outputs of pumped storage are shown in Figure 20, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power. The power output comparisons of pumped storage in Case 3 and Case 4 are shown in Figure 21.

Results in Case 4
The power outputs of thermal generators based on the SampEn in Case 4 are shown in Figure 19. The power outputs of pumped storage are shown in Figure 20, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power. The power output comparisons of pumped storage in Case 3 and Case 4 are shown in Figure 21.
The power outputs of thermal generators based on the SampEn in Case 4 are shown in Figure  19. The power outputs of pumped storage are shown in Figure 20, where the positive power means the operation state is as loads absorbing power from grid and the negative power means the operation state is generating power. The power output comparisons of pumped storage in Case 3 and Case 4 are shown in Figure 21.   The operation costs, total up and down ramping power of thermal generators in Case 1 and Case 2 are shown in Table 4 together with the throughput power of pumped storage.   The operation costs, total up and down ramping power of thermal generators in Case 1 and Case 2 are shown in Table 4 together with the throughput power of pumped storage.

Result Comparison of Case 1 and Case 2
The operation costs, total up and down ramping power of thermal generators in Case 1 and Case 2 are shown in Table 4 together with the throughput power of pumped storage. The operation costs, throughput power of pumped storage, total up and down ramping power of thermal generators in Case 3 and Case 4 are shown in Table 5.

Power Outputs Analysis of Thermal Generators
The power output of thermal generators in Case 1-Case 4 are shown in Figures 12, 14, 17 and 19, respectively. The power outputs of Unit 8 and Unit 9 are similar thus the power output curve of Unit 9 is covered by the one of Unit 8. The optimization in Case 1 and Case 3 aim at minimizing the operational costs. Apart from the traditional operation constraints, the power output plan also pursues the maximum economy and responses to the uncertainty of renewable power energies. Through the net load tracing and SampEn calculation, the power outputs are as Figures 14 and 19 shows. For the time periods with high complexity level net loads, the power outputs of thermal generators are limited according to SampEn. For example, as Table 3 shows, in the second time period (duration is 5 h) and the fourth time period (duration is 7 h), the SampEn are large. The thermal generators are limited and forced to operate continuously and stability. Therefore, the operation is kept stable and the generator efficiency is increased together with the ramping power in these time periods reduced.
In the first, third, and fifth time periods, the SampEn is small, and there will be obviously power increase or decrease and last for hours. Thus, thermal generators need to respond to these power changes. Through the adjustment based on SampEn, the generating mode of thermal generators is settled as "3-5-3-7-6" in this simulation. Moreover, the time scales of the generating mode are corresponding to the results of SampEn.

Power Output Analysis of Pumped Storage
As shown in Figures 13 and 18, at the initial period of time, in order to respond to the load increases in the second time period, the pumped storage in Case 1 and Case 3 should be activated and absorb power from power grid. The absorption of the pumped storage would store the redundant power from the grid for the net load peak time, during which the most amount of energies are required. Thus, the overall system efficiency is increased. In Case 2 and Case 4 as shown in Figures 15 and 20, at the initial period of time, the pumped storage absorbs power from power grid. However, the power outputs of thermal generators are kept stable or changes slightly. Thus, according to the Equations (6) and (7), the power is absorbed from the renewable energy in the net loads. As Figures 16 and 21 show, during the time periods with high net load complexity level, the throughput of pumped storage is increased in Case 2 and Case 4. The power output adjustment of pumped storage response to the fluctuation of net loads instead of thermal generators. In addition, according to the generating mode of "3-5-3-7-6", during the "5" and "7" time periods, the thermal power is kept stable and the changes of pumped storage take full response to the fluctuation of net loads. Thus, this strategy increases the interaction between pumped storage and renewable energies and stabilizes the thermal power, which further improves the stability of the system. Table 4 shows that, the total operation cost of thermal generators in Case 2 is slightly less than that in Case 1. The result in Case 2 based on SampEn is only 0.04% less than in Case 1, but the operation status of thermal generators has been optimized. In Case 2, the total up and down ramping power have been both optimized up to around 43% compared to the results in Case 1. Moreover, the throughput of pumped storage in Case 2 is 44% more than that in Case 1. Table 5 shows that, the total operation cost of thermal generators in Case 4 is less than that in Case 3. The result in Case 4 based on SampEn is only 0.25% less than in Case 3, but the total up and down ramping power have been both optimized up to around 13% compared to the results in Case 3. Moreover, the throughput of pumped storage in Case 4 is 28% more than that in Case 3.

Power Output Analysis of the Cases with and Without SampEn
The overall results show that through the net load tracing and SampEn calculation, the original operation status with large fluctuation is transformed into the generating mode of "3-5-3-7-6". Thus, the continuity and stability operation of thermal generators are realized and the operation efficiency of the system is improved. The power dispatch strategy based on SampEn can reduce the ramping power of thermal generators while ensuring both the economy and safety of the system.

Conclusions
This paper has presented a net load tracing dispatch strategy based on SampEn focusing on the high-percentage renewable energy connected to the power grid. A specific generating mode of thermal generators is determined by the net load tracing and SampEn. The conclusions based on the simulation results are as follows: 1. In the strategy, the renewable energies have priority to connect to power grid and the net loads are the dispatch objects. The renewable energies consumed by traditional loads totally and the curtailment of the wind and photovoltaic power are prevented. Thus, the utilization efficiency of renewable energies is maximized.
2. According to the net load tracing, the SampEn reflecting the complexity of net loads is calculated depending on the analysis of the characteristics and numeric features of net loads. Moreover, based on SampEn, the generating mode of thermal generators, the main power sources response to the fluctuation of renewable energies and their generating and ramping characteristics in different time periods are determined.
3. The interaction between pumped storage and renewable energies is increased in the dispatch based on SampEn. The pumped storage is responsible for the fluctuation of renewable energies instead of thermal generators in the certain time periods, during which the net loads are with large SampEn. Thus, the throughput of pumped storage is increased 44% in Case 2 and 27% in Case 4.
4. In this proposed strategy, the power outputs of thermal generators are smoothed and the generating time periods are optimized. The duration of stable operation of the thermal generators is increased and the ramping power during the time periods, in which the net loads are with large SampEn, are reduced 43% in Case 2 and 13% in Case 4. Overall, the stability of the system is improved with the economy of the system is ensured.