1. 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 the imbalance between electric loads and intermittent power energies. In this way, the source–load bilateral uncertainty problem is transformed to how to respond to the uncertainty of net loads.
There are many studies focusing on the power balancing problems related to high-proportion renewable energies connected to the power grid. The multi-source coordination of power systems is an effective means. The flexible resources such as hydropower plant, gas turbine power plant, energy storage, and so forth are developed and utilized to respond to the uncertainty and fluctuation of renewable energies cooperatively [9,10]. In [10], a robust unit commitment model considering the cooperation of wind power and pumped storage is proposed, in which the wind power fluctuation is restrained by the pumped storage. Reference [11] utilized the second life battery energy storage to assist thermal generators to respond to the fluctuation of renewable energies. Facing the high-percentage renewable energies, the ramp capability of generators is an important factor when to response to the load demands [12,13,14]. It is argued in [15] that the insufficient ramp resources in the system may lead to the insufficient flexibility to meet variations in the net loads. Reference [16] estimated the potential variability of net loads by the statistical approaches on the ramping characteristics to determine the required ramping capability in the future. The increasing ramping capability can improve the flexibility on supply side. However, the above studies are based on the traditional power dispatch. With the increasing of renewable energies percentage, the amplitude and frequency of generating and ramping power adjustments are still increased, which will cause security risks of system operation. Power dispatch contains different time-scales. The horizon optimization of day-ahead power dispatch is 24 h [17,18]. The advance dispatch operates in the framework of a short time, such as 30 min [17]. The shorter time-scale can be state estimation, which contains load power measurement, estimation, and dispatching [19,20,21,22,23]. In order to deal with the problems of renewable energies connected to the power grid and the ramping constraints of generators, there are studies focusing on the dispatch from the perspective of time dimension. Reference [24] proposes a model considering operational flexibility at sub-hourly time-scales. Reference [25] provides a multi-time scale coordination power dispatch divided into day-ahead, intra-day, and real-time scale. Reference [26] studies a unit commitment strategy aiming at improving the ramping capability of generators. The strategy provides multiple time-scale ramping constraints to satisfy the demand of load prediction, in which the ramping capability can cover the whole varies range of the uncertain renewable energies. However, sufficient ramping reserve will lead to wasting of resources.
Besides, the description method of net load uncertainties affects the determination of dispatch strategy. The uncertainty characteristics of supply side and load side can be described and analyzed respectively [27,28] further to stabilize the uncertainty of net loads. The stochastic programming model [29] or fuzzy program model [27,30] of the net loads are built and solved by fitting optimization methods, such as the sample-based stochastic decomposition approach [31] and the Benders decomposition method [32]. In [33], the Fourier transform is applied to calculate the convolution of load and wind power and compute the probability distribution of joint variables. These studies, however, are still focused on the assumption and simulation of the distribution characteristics of the power and loads. There is not much attention on numeric characteristics analysis. Reference [34] analyzed the duration curves, volatility, and hourly ramps of net loads and indicated that the characteristic of net loads plays an important role in the power grids with high penetration of renewables. The uncertainty and disorder of net load can be analyzed and measured by the concept of entropy [35,36,37]. The approximate entropy can be utilized to analyze the characteristics of data and evaluate the complexity level of the time series, and it is widely used in many fields [38]. However, the consistency of approximate entropy is poor and during the solution process, there are self-data comparisons and deviations. These weaknesses promote the utilization of sample entropy (SampEn) developed by Richman [39,40], which has better applicability and can evaluate the complexity level more accurately [41,42].
In this paper, a sample entropy based net load tracing dispatch strategy with a specific thermal generating mode is provided to increase the operation efficiency of thermal generators and the 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.
2. Characteristics of Net Loads and SampEn Calculation 2.1. 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 Figure 1 and Figure 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. 2.2. 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 theSampEnNetLoadof the net load time series is shown as follows:
1. The whole net load time series is given:{NLl}={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{NLl}is shown in Equation (1).
Slope(l)=nl(l+1)−nl(l)
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.
symbollSlope=sign(Slope(l))
{ifsymbollSlope=symboll+1Slope,translSlope=0ifsymbollSlope≠symboll+1Slope,translSlope=1
4. The ratio of sign changing amount to net load amount is rolling calculated during the timetl~tl+k. The equation is shown as
Ratiol,l+kSlope=countk=1~N−l+1(translSlope, transl+kSlope)k
5. The ratio ofRatiol,l+kSlopeto the valley-to-peak of the net loadsSupl,l+kSampEnis rolling calculated during the timetl~tl+k. The calculated value is as the auxiliary value of the net load time series. The equation is shown in (5).
Supl,l+kSampEn=Ratiol,l+kSlopemaxk=1~N−l+1(nl(l),nl(l+k))−mink=1~N−l+1(nl(l),nl(l+k))
The piecewise estimation of the net load series is conducted based theSupl,l+kSampEn. The similarSupl,l+kSampEnare divided in the same time period, and the piecewise time pointstg,g∈Nare achieved. Thus, the subsequence of net load is expressed as{nl(g),nl(g+1),…,nl(g+TNL)}, whereTNLis the length of the net load subsequence.
6. The time series{PiNetLoad}={pNL(1),pNL(2),…,pNL(TNL)}is assumed as the subsequence of the net load at timetg. A subsequence with dimension ofmNLis generated as Equation (6) shows
PNL(i)=[pNL(i),pNL(i+1),…,pNL(i+mNL−1)]
where,i=1,2,…,TNL−mNL+1andmNLis usually equal to 2 or 3.
7. TheDNLij mNL(PNL(i),PNL(j))is defined as the distance of the subsequencePNL(i)andPNL(j), wherej=1,2,…,TNL−mNL+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 betweenPNL(i)andPNL(j)is shown in Equation (7).
dijmNL(PNL(i),PNL(j))=[|u(i)−u(j)||u(i+1)−u(j+1)|⋮|u(i+mNL−1)−u(j+mNL−1)|]
Thus, the distance of the subsequencePNL(i)andPNL(j)is calculated as Equation (8) shows.
DNLij mNL(PNL(i),PNL(j))=max(dimNL(PNL(i),PNL(j)))
8. The tolerance levelrNLis a threshold defining when two states are similar, whererNL>0. In this paper, therNLwas set0.2×STDNL,STDNLis the standard deviation ofNL. To each i, the amount ofDNLij mNL(PNL(i),PNL(j))<rNLare calculated.BNLimNL(rNL)is defined as Equation (9):
BNLimNL(rNL)=count(DNLij mNL(PNL(i),PNL(j))<rNL)TNL−mNL
whereTNL−mNLis the total distance of the sample.
9. The actual average over all the vectorsBNLimNL(rNL)of is computed as Equation (10) shows.
BAvgmNL(rNL)=∑i=1TNL−mNL+1BNLimNL(rNL)TNL−mNL+1
10. The dimension of subsequence is increased to m + 1. Repeat steps 6 to 10 and achieveBAvgmNL+1(rNL). Theoretically, the sample entropy is defined as
SampEnNetLoad(rNL,mNL)=limTNL→∞{−ln(BAvgmNL+1(rNL)BAvgmNL(rNL))}
In practice, the length of time seriesTNLis finite. Therefore, the estimation of SampEn is denoted as Equation (12).
SampEnNetLoad(rNL,mNL,TNL)=−ln(BAvgmNL+1(rNL)/BAvgmNL(rNL))
2.3. 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 Figure 4 and Figure 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 Figure 4 and Figure 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.
3. New Energy Net Load Tracing Dispatch Strategy Based on SampEn 3.1. 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 SampEnSampEntg tlNetLoadof net load seriesPNLtg tlin time periodstgtotl is calculated as Section 2.2 shows. The generating mode setting is depended on the value ofSampEntg tlNetLoad. IfSampEntg tlNetLoadis large, the complexity level in this time period is high and the fluctuation of thermal generation are minimized. IfSampEntg tlNetLoadis 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.
{SampEntg tlNetLoad is big:|max[Ppi,tg:Ppi,tl]−min[Ppi,tg:Ppi,tl]|≤SampEntg tlNetLoadSampEntg tlNetLoad is small:Ppimin≤Ppi,t≤Ppimax(tg, tl∈t)
wherePpi,tis the real power output of unit pi at time t.PpiminandPpimaxare the minimum and maximum power outputs of the ith thermal generator, respectively.
3.2. Power Dispatch Model Based on SampEn
3.2.1. Objective Functions
The thermal generation costFthermalis considered as the objective function in this paper, which is calculated by
Fthermal=∑t=1T∑ping(api×Ppi,t2+bpi×Ppi,t+cpi)
where T is total number of time periods; ng is the number of dispatchable units;api,bpi, andcpiare the coefficients of fuel cost functions for units pi.
3.2.2. Power Balance Equations and Constraint Functions
1. Power Balance Equations
The power balance equations at time t are formulated as
∑pi=1ngPpi,t+∑Hgeni=1npPHgeni,t-PNLt−∑Hpumpi=1npPHpumpi,t=0
PNLt=PLt−∑wpi=1nwPwpi,tpre−∑solari=1nsPsolari,tpre
where∑pi=1ngPpi,tis the total generator power output during the tth time period;PLtis the total load demand at time t.PHgeni,tis the generating power of pumped storageHgeniat time t andPHpumpi,tis the pumping power of pumped storageHpumpiat time t.PNLtis the total net load demands at time t.Pwpi,tpreis the predicted wind power of wind farm wpi at time t andPsolari,tpreis the predicted photovoltaic power of photovoltaic station solari at time t.np,nw, andnsare 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).
PHgenimin≤PHgeni,t≤PHgenimax
PHpumpimin≤PHpumpi,t≤PHpumpimax
wherePHgeniminandPHgenimaxare the minimum and maximum generating power of pumped storageHgeni, respectively.PHpumpiminandPHpumpimaxare the minimum and maximum pumping power of pumped storageHpumpi, respectively. The constraints of the reservoir energy conversion are shown as follows:
HPumpi,t+1=HPumpi,t+Δt×(PHpumpi,t×γPump−PHgeni,tγGen)
HPumpimin≤HPumpi,t≤HPumpimax
whereHPumpi,t+1andHPumpi,tare the storage condition of the reservoir at time t + 1 and time t.Δtis the time interval.γPumpis the pumping efficiency andγGenis the generating efficiency.HPumpiminandHPumpimaxthe 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.
PHgeni,t·PHpumpi,t=0
3. Constraints of Thermal Generators
The power output constraints of thermal generators are shown in Equation (22).
Ppimin≤Ppi,t≤Ppimax
The generating ramp rate limits are formulated by
{Ppi,t−Ppi,t−1≤URpiPpi,t−1−Ppi,t≤DRpi
whereURpiandDRpiare 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{∑pi=1ngPURi,t≥∑wpi=1nw(Pwpi,tpre−Pwpi,t)}≥ρ
Pr{∑pi=1ngPDRi,t≥∑i=1nw(Pwpi,t−Pwpi,tpre)}≥ρ
PURi,t≤min(Ppimax−Ppi,t,URpi)
PDRi,t≤min(Ppi,t−Ppimin,DRpi)
whereρis the confidence coefficient andPwpi,tis the actual wind power of wind farm wpi at time t.PURi,tandPDRi,tare the actual up and down reserve of thermal generator pi.
3.2.3. 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
fw(pwp)=pwpα−1 (1−pwp)β−1B(α,β)
wherepwpis a normalized wind power output,B(α,β)is the Beta distribution function andα,βare the distribution shape parameters. The calculated formulas are
pwp=Pwp−PwpminPwpmax−Pwpmin,pwp∈[0,1]
B(α,β)=∫01pwpα−1×(1−pwp)β−1dpwp
E(pwppre)=E(PwpprePwpmax)=αα+β
D(pwppre)=D(PwpprePwpmax)=αβ(α+β)2(α+β+1)
wherePwpis the actual wind power output;PwpminandPwpmaxare the minimum and maximum outputs of wind turbines, respectively.Pwppreis the predicted wind power;E(pwppre)is the expectation ofPwppre; andD(pwppre)is the variance ofPwppre.
The fractile is used to solve the above chance constraints [46].
Pr(Y>Yρ)=ρρ∈(0,1)
where Y is the random variable andYρis the fractile ofρ . In Figure 6, theYρis the upper fractile and theY1−ρ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, theYρandY1−ρ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.
Pr{∑wpi=1nwPwpi,t≥−∑pi=1ngPURi,t+∑wpi=1nwPwpi,tpre}≥ρ
Pr{∑wpi=1nwPwpi,t≤∑pi=1ngPDRi,t+∑wpi=1nwPwpi,tpre}≥ρ
−∑pi=1ngPURi,t+∑wpi=1nwPwpi,tpre≤∑wpi=1nwYρi
∑pi=1ngPDRi,t+∑wpi=1nwPwpi,tpre≥∑wpi=1nwY1−ρi
3.3. 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 ofSampEntg tlNetLoad, 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.
4. 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, wherePmaxandPminare 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;CoeURiandCoeDRiare the up and down reserve cost coefficients, respectively;Coeramp−upiandCoeramp−downi 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 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.
4.1. 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.
4.2. 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.
4.2.1. 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.
4.2.2. 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.
4.2.3. 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.
4.2.4. 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.
4.2.5. 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.
4.2.6. Result Comparison of Case 3 and Case 4
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.
4.3. Discussion
4.3.1. Power Outputs Analysis of Thermal Generators
The power output of thermal generators in Case 1–Case 4 are shown in Figure 12, Figure 14, Figure 17 and Figure 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 Figure 14 and Figure 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.
4.3.2. Power Output Analysis of Pumped Storage
As shown in Figure 13 and Figure 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 Figure 15 and Figure 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 Figure 16 and Figure 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.
4.3.3. Power Output Analysis of the Cases with and Without SampEn
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.
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. 5. 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.
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
[Image omitted; see PDF]
| Date | Time | Slope Sign Changing Amount | Percentage of Slope Sign Changing | Valley-to-Peak | Slope Sign Changing Amount/Valley-to-Peak | SampEn |
|---|---|---|---|---|---|---|
| 0411 | 1–18 | 11 | 0.61 | 731.61 | 0.0150 | 0.42 |
| 19–41 | 3 | 0.13 | 457.95 | 0.0066 | 0.18 | |
| 42–69 | 15 | 0.54 | 247.65 | 0.0606 | 0.74 | |
| 70–96 | 1 | 0.04 | 795.91 | 0.0013 | 0.08 | |
| 0303 | 1–45 | 6 | 0.13 | 2582.03 | 0.0023 | 0.04 |
| 46–96 | 18 | 0.35 | 340.37 | 0.0529 | 0.72 |
| Parameters | Unit1 | Unit2 | Unit3 | Unit4 | Unit5 | Unit6 | Unit7 | Unit8 | Unit9 | Unit10 |
|---|---|---|---|---|---|---|---|---|---|---|
| Pmax(MW) | 470 | 460 | 340 | 300 | 243 | 160 | 130 | 120 | 80 | 55 |
| Pmin (MW) | 150 | 135 | 73 | 60 | 73 | 57 | 20 | 47 | 20 | 20 |
| a (10−3 $/MW2h) | 0.43 | 0.63 | 0.39 | 0.70 | 0.79 | 0.56 | 2.11 | 4.80 | 109.08 | 9.51 |
| B ($/MWh) | 21.60 | 21.05 | 20.81 | 23.90 | 21.62 | 17.87 | 16.51 | 23.23 | 19.58 | 22.54 |
| C ($/h) | 958.20 | 1313.6 | 604.97 | 471.60 | 480.29 | 601.75 | 502.70 | 639.40 | 455.60 | 692.40 |
| UR | 120 | 120 | 120 | 100 | 100 | 100 | 50 | 50 | 50 | 50 |
| DR | 120 | 120 | 120 | 100 | 100 | 100 | 50 | 50 | 50 | 50 |
| CoeURi ($/MWh) | 14.7 | 15.5 | 15.2 | 17.8 | 19.3 | 19.8 | 18.7 | 21.7 | 23.4 | 25.2 |
| CoeDRi ($/MWh) | 15.2 | 14.8 | 15.1 | 18.6 | 21.2 | 19.5 | 19 | 22 | 23.1 | 25.6 |
| Coeramp−upi ($/MWh) | 3.13 | 3.08 | 3.75 | 4.17 | 5.88 | 9.71 | 9.09 | 13.7 | 16.67 | 28.57 |
| Coeramp−downi ($/MWh) | 3.13 | 3.08 | 3.75 | 4.17 | 5.88 | 9.71 | 9.09 | 13.7 | 16.67 | 28.57 |
| Time Periods (h) | SampEn | SampEn Proportion |
|---|---|---|
| 1–4 | 0.04 | 5% |
| 4–8 | 0.14 | 16% |
| 8–12 | 0.08 | 2% |
| 12–18 | 0.61 | 71% |
| 18–24 | 0.06 | 6% |
| Parameters | Case 1 | Case 2 | Percentage Optimization of Case 2 Compared to Case 1 |
|---|---|---|---|
| Operation cost (105 $) | 7.1744 | 7.1712 | 0.04% |
| Up ramping power (MW) | 1636.27 | 918.59 | 43.86% |
| Down ramping power (MW) | 1526.61 | 869.37 | 43.05% |
| Throughput of pumped storage (MW) | 1636.03 | 2358.29 | 44.15% |
| Parameters | Case 3 | Case 4 | Percentage Optimization of Case 4 Compared to Case 3 |
|---|---|---|---|
| Operation cost (105 $) | 7.2082 | 7.1903 | 0.25% |
| Up ramping power (MW) | 1932.88 | 1681.89 | 12.99% |
| Down ramping power (MW) | 1860.82 | 1623.47 | 12.76% |
| Throughput of pumped storage (MW) | 1955.61 | 2497.82 | 27.73% |
Author Contributions
Data curation, F.P.; Formal analysis, S.H. and C.D.; Methodology, S.H.; Project administration, W.Z.; Supervision, H.S.; Validation, Z.G. and W.Z.; Writing-Original draft, S.H.; Writing-Review and editing, H.S.
Funding
This research received no external funding.
Acknowledgments
This work is supported by The Science and Technology Project of State Grid (2017YF-20, 2017YF-27).
Conflicts of Interest
The authors declare no conflict of interest.
1. Zhou, X.X.; Chen, S.Y.; Lu, Z.X.; Huang, Y.H.; Ma, S.C.; Zhao, Q. Technology features of the new generation power system in China. Proc. CSEE 2018, 38, 1893-1904.
2. Krishnan, S.H.; Ezhilarasi, D.; Uma, G.; Umapathy, M. Pyroelectric-based solar and wind energy harvesting system. IEEE Trans. Sustain. Energy 2014, 5, 73-81.
3. Hsiao, C.C.; Jhang, J.W.; Siao, A.S. Study on pyroelectric harvesters integrating solar radiation with wind power. Energies 2015, 8, 7465-7477.
4. Xu, F.; Li, L.L.; Chen, Z.X.; Tu, M.F.; Ding, Q. Generation scheduling model and application with fluctuation reduction of unit output. Autom. Electr. Power Syst. 2012, 36, 45-50+102.
5. Wood, A.J.; Wollenberag, B.F. Power Generation Operation and Control; Wiley: New York, NY, USA, 1996.
6. Bai, J.H.; Xin, S.X.; Liu, J.; Zheng, K. Roadmap of realizing the high penetration renewable energy in China. Proc. CSEE 2015, 35, 3699-3705.
7. Huber, M.; Dimkova, D.; Hamacher, T. Integration of wind and solar power in Europe: Assessment of flexibility requirements. Energy 2014, 69, 236-246. [Green Version]
8. Min, C.G.; Kim, M.K. Net load carrying capability of generating units in power systems. Energies 2017, 10, 1221.
9. Zhou, B.R.; Geng, G.C.; Jiang, Q.Y. Hydro-thermal-wind coordination in day-ahead unit commitment. IEEE Trans. Power Syst. 2016, 31, 4626-4637.
10. Jiang, R.; Wang, J.; Guan, Y. Robust unit commitment with wind power and pumped storage hydro. IEEE Trans. Power Syst. 2012, 27, 800-810.
11. Hu, S.B.; Sun, H.; Peng, F.X.; Zhou, W.; Cao, W.P.; Su, A.L.; Chen, X.D.; Sun, M.Z. Optimization strategy for economic power dispatch utilizing retired EV batteries as flexible loads. Energies 2018, 11, 1657.
12. Ortega-Vazquez, M.A.; Kirschen, D.S. Estimating the spinning reserve requirements in systems with significant wind power generation penetration. IEEE Trans. Power Syst. 2009, 24, 114-124.
13. Papavasiliou, A.; Oren, S.S.; O'Neill, R.P. Reserve requirements for wind power integration: A scenario-based stochastic programming framework. IEEE Trans. Power Syst. 2011, 26, 2197-2206.
14. Matos, M.A.; Bessa, R.J. Setting the operating reserve using probabilistic wind power forecasts. IEEE Trans. Power Syst. 2011, 26, 594-603.
15. Lannoye, E.; Flynn, D.; O'Malley, M. Assessment of power system flexibility: A high-level approach. In Proceedings of the IEEE Power & Energy Society General Meeting, San Diego, CA, USA, 22-26 July 2012; pp. 1-8.
16. Wu, Y.K.; Chang, L.T.; Su, P.E.; Hsieh, T.Y.; Jan, B.S. Assessment of potential variability of net load following the integration of 3 GW wind power in Taiwan. Energy Procedia 2016, 100, 117-121.
17. Marannino, P.; Granelli, G.P.; Montagna, M.; Silvestri, A. Different time-scale approaches to the real power dispatch of thermal units. IEEE Trans. Power Syst. 1990, 5, 169-176.
18. Generation Dispatching Plan of Electricity Market on Generation Side. Available online: http://tech.bjx.com.cn/html/20071227/59347.shtml (accessed on 25 December 2018).
19. Cataliotti, A.; Cervellera, C.; Cosentino, V.; Di Cara, D.; Gaggero, M.; Maccio, D.; Marsala, G.; Ragusa, A.; Tine, G. An Improved Load Flow Method for MV Networks Based on LV Load Measurements and Estimations. IEEE Trans. Instrum. Meas. 2018, 68, 430-438.
20. Cataliotti, A.; Cosentino, V.; Di Cara, D.; Tinè, G. LV Measurement Device Placement for Load Flow Analysis in MV Smart Grids. IEEE Trans. Instrum. Meas. 2016, 65, 999-1006.
21. Xygkis, T.C.; Korres, G.N. Optimized measurement allocation for power distribution systems using mixed integer sdp. IEEE Trans. Instrum. Meas. 2017, 66, 2967-2976.
22. Xygkis, T.C.; Korres, G.N.; Manousakis, N.M. Fisher information based meter placement in distribution grids via the d-optimal experimental design. IEEE Trans. Smart Grid 2018, 9, 1452-1461.
23. Pegoraro, P.A.; Angioni, A.; Pau, M.; Monti, A.; Muscas, C.; Ponci, F.; Sulis, S. Bayesian approach for distribution system state estimation with non-gaussian uncertainty models. IEEE Trans. Instrum. Meas. 2017, 66, 2957-2966.
24. Zhang, L.; Yang, J.; Jian, X.H.; Zhang, F.; Han, X.S. Unit Commitment with Energy Storage Considering Operation Flexibility at Sub-hourly Time-scales. Autom. Electr. Power Syst. 2018, 42, 48-56.
25. Zhai, J.Y.; Ren, J.W.; Zhou, M.; Li, Z. Multi-Time Scale Fuzzy Chance Constrained Dynamic Economic Dispatch Model for Power System with Wind Power. Power Syst. Technol. 2016, 40, 1094-1099.
26. Yong, T.; Yao, J.; Yang, S.; Yang, Z. Ramp enhanced unit commitment for energy scheduling with high penetration of renewable generation. In Proceedings of the 2015 IEEE Power & Energy Society General Meeting, Denver, CO, USA, 26-30 July 2015.
27. Zhao, D.M.; Yin, J.F. Fuzzy random chance constrained preemptive goal programming scheduling model considering source-side and load-side uncertainty. Trans. China Electrotech. Soc. 2018, 33, 1076-1085.
28. Li, P.; Wand, X.Y.; Han, P.F. Risk Analysis of Microgrid Optimal Operation Scheduling under Double Uncertainty Environment. Proc. CSEE 2017, 37, 4296-4303.
29. Guo, P.; Wen, J.; Zhu, D.D.; Wang, W.Z.; Liu, W.Y. The Coordination Control Strategy for Large-Scale Wind Power Consumption Based on Source-Load Interactive. Trans. China Electrotech. Soci. 2017, 32, 1-9.
30. Zhang, X.H.; Zhao, C.M.; Liang, J.X.; Li, K.; Zhong, J.Q. Robust Fuzzy Scheduling of Power Systems Considering Bilateral Uncertainties of Generation and Demand Side. Autom. Electr. Power Syst. 2018, 42, 67-75.
31. Gangammanavar, H.; Sen, S.; Zavala, V.M. Stochastic optimization of sub-hourly economic dispatch with wind energy. IEEE Trans. Power Syst. 2016, 31, 949-959.
32. Wang, J.D.; Wang, J.H.; Liu, C.; Ruiz, J.P. Stochastic unit commitment with sub-hourly dispatch constraints. Appl. Energy 2013, 105, 418-422.
33. Zhang, X.Y.; Yang, J.Q.; Zhang, X.J. Dynamic economic dispatch incorporating multiple wind farms based on FFT simplified chance constrained programming. J. Zhengjiang Univ. (Eng. Sci.) 2017, 51, 976-983.
34. Shaker, H.; Zareipour, H.; Wood, D. Impacts of large-scale wind and solar power integration on California's net electrical load. Renew. Sustain. Energy Rev. 2016, 58, 761-774.
35. Yin, Y.; Shang, P.J.; Andrew, C.A.; Peng, C.K. Multiscale joint permutation entropy for complex time series. Physica A 2019, 515, 388-402.
36. Cuadras, A.; Romero, R.; Ovejas, V.J. Entropy characterisation of overstressed capacitors for lifetime prediction. J. Power Sources 2016, 336, 272-278.
37. Hsiao, C.C.; Liang, B.H. The generated entropy monitored by pyroelectric sensors. Sensors 2018, 18, 3320.
38. Zhang, Y.Y. The Theory of Approximate Entropy and its Application. Chin. J. Med. Phys. 2009, 26, 1543-1546.
39. Chou, C.M. Complexity analysis of rainfall and runoff time series based on sample entropy in different temporal scales. Stoch. Environ. Res. Risk Assess. 2014, 28, 1401-1408.
40. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 278, 2039-2049.
41. Zhang, Y.C.; Liu, K.P.; Qin, L.; Fang, N.C. Wind Power Multi-Step Interval Prediction Based on Ensemble Empirical Mode Decomposition-Sample Entropy and Optimized Extreme Learning Machine. Power Syst. Technol. 2016, 40, 2045-2051.
42. Al-Angari, H.M.; Sahakian, A.V. Use of sample entropy approach to study heart rate variability in obstructive sleep apnea syndrome. IEEE Trans. Bio-Med. Eng. 2007, 54, 1900-1904.
43. Xu, W. The Complexity Algorithms Research of Chaotic Sequences Based on Entropy Theory. Ph.D. Thesis, Heilongjiang University, Harbin, China, 2017.
44. Zhang, X.Q.; Liang, J.; Zhang, X.; Zhang, F.; Zhang, L.; Xu, B. Combined Model for Ultra Short-term Wind Power Prediction Based on Sample Entropy and Extreme Learning Machine. Proc. CSEE 2013, 33, 33-40.
45. Cheng, W.S.; Zhang, H.F. A dynamic economic dispatch model incorporating wind power based on chance constrained programming. Energies 2015, 8, 233-256.
46. Shen, Z.; Xie, S.Q.; Pan, C.Y. Probability and Statistics; High Education Press: Beijing, China, June 2008; pp. 139-143, 161-163.
47. Attaviriyanupap, P.; Kita, H.; Tanaka, E.; Hasegawa, J. A hybrid EP and SQP for dynamic economic dispatch with nonsmooth fuel cost function. IEEE Trans. Power Syst. 2002, 17, 411-416.
48. Zhang, X.H.; Zhao, J.Q.; Chen, X.Y. Multi-objective Unit Commitment Fuzzy Modeling and Optimization for Energy-saving and Emission Reduction. Proc. CSEE 2010, 30, 71-76.
Shubo Hu, Feixiang Peng, Zhengnan Gao, Changqiang Ding, Hui Sun* and Wei Zhou
School of Electrical Engineering, Dalian University of Technology, Dalian 116024, China
*Author to whom correspondence should be addressed.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2019. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
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 specific 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 flexibly. According to the SampEn, a specific 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 efficiency and reasonability of the proposed dispatch strategy are confirmed.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
; Peng, Feixiang
; Gao, Zhengnan; Ding, Changqiang; Sun, Hui
; Zhou, Wei