Control Strategy of a Hybrid Energy Storage System to Smooth Photovoltaic Power Fluctuations Considering Photovoltaic Output Power Curtailment

The power fluctuations of grid-connected photovoltaic (PV) systems have negative impacts on the power quality and stability of the utility grid. In this study, the combinations of a battery/supercapacitor hybrid energy storage system (HESS) and the PV power curtailment are used to smooth PV power fluctuations. A PV power curtailment algorithm is developed to limit PV power when power fluctuation exceeds the power capacity of the HESS. A multi-objective optimization model is established to dispatch the HESS power, considering energy losses and the state of charge (SOC) of the supercapacitor. To prevent the SOCs of the HESS from approaching their lower limits, a SOC correction strategy is proposed to correct the SOCs of the HESS. Moreover, this paper also investigates the performances (such as the smoothing effects, losses and lifetime of energy storage, and system net profits) of two different smoothing strategies, including the method of using the HESS and the proposed strategy. Finally, numerous simulations are carried out based on data obtained from a 750 kWp PV plant. Simulation results indicate that the proposed method is more economical and can effectively smooth power fluctuations compared with the method of using the HESS.


Introduction
Photovoltaic (PV) plays a vital role in the field of renewable energy, and a large number of PV plants have been built in various countries, such as China, US, and Germany [1].One of the main characteristics of PV plants is the significant variability of their output power because of the changes in solar radiation.For example, statistics show that the maximum power fluctuation of a PV plant can reach 60%/s of its installed capacity [2].As the penetration of PV increases, if the PV power is directly injected into the utility grid without any control, such high-frequency power fluctuations will cause voltage fluctuations, thereby reducing the reliability and power quality of the power system [3,4].Therefore, the power fluctuations of PV should be limited to a certain range, for example, the Puerto Rico Electric Power Authority (PREPA) requires that the maximum power fluctuation of a PV plant should be less than 10%/min of its installed capacity [5].
Many previous works have investigated numerous methods that can be used to smooth the PV power fluctuations.It has been shown that connecting several PV plants that are far apart can reduce the overall volatility of PV power [6], and it is a simple method to smooth power fluctuations.The controllable resistors can also be used to smooth power fluctuations, and the excess power fluctuations are absorbed by controlling the resistance values [7].The power variation of different PV array configurations has found to be different, and PV power fluctuations can be decreased by increasing the maximum array dimension [8].In addition, PV power curtailment [3] (or PV ramp-rate control [9,10]) is a convenient approach to smooth PV power fluctuations.The main idea of this approach is that the upward power fluctuation can be reduced by adjusting the actual operating point of the PV converter below its maximum power point (MPP) [3,9,10].However, this smoothing method cannot smooth downward power fluctuations caused by the fall of solar radiation.Moreover, the energy losses caused by limiting PV power need further study.
Recently, there has been an increasing interest in using battery energy storage (BES) or a battery/supercapacitor hybrid energy storage system (HESS) to smooth PV power fluctuations at the point of common coupling (PCC) [5,[9][10][11][12].For example, a step control strategy for BES was proposed [5] to smooth the short-term power fluctuations of PV.A state of charge (SOC)-based smoothing strategy was developed in [12], which fully discussed how to dispatch the output power of multiple independent energy storage units.However, the disadvantages of batteries are a lower power rate and low lifetime.Thus, when BES is used to smooth PV power fluctuations, the battery capacity needs to be large, resulting in high costs.Fortunately, previous studies [3,13] have found that combining BES and the PV power curtailment is an economical solution for smoothing PV power fluctuations compared with the method of using BES because this approach can significantly reduce the capacities of batteries and prolong the service life of batteries.However, the main weakness of this approach is that limiting PV power will decrease the revenues of selling PV energy [3,7,13].
Moreover, the HESS has better comprehensive performances (e.g., power rate, lifetime, and efficiency) than that of BES [14].Therefore, a power-sharing strategy of the HESS was developed to smooth the power fluctuations of PV [15], in which power spikes are absorbed by supercapacitor energy storage (SCES), and BES only needs to provide stable energy, hence the battery capacity can be reduced and its service lifetime can also be prolonged [15].However, the number of simulation cases of this study [15] is relatively small, and these cases may not represent the worst PV power fluctuations.
However, what is not yet clear is how the HESS benefits from the PV power curtailment when the HESS is used to smooth PV power fluctuations.There are two primary aims of this study: (1) To investigate the coordinated control strategy for the HESS and PV to smooth PV power fluctuations, and (2) to analyze the performances of the proposed smoothing strategy and the smoothing method of using the HESS.
The main contributions of this work are as follows: (1) A control strategy for the HESS and PV is proposed to smooth power fluctuations of PV.When the power demand, which is used to smooth the upward power fluctuation, exceeds the maximum allowable power of the HESS, the excess power demand will be absorbed by reducing the PV power.Otherwise, only the HESS is used to smooth both the upward and downward power fluctuations.(2) A multi-objective optimization model is developed to dispatch the power demand of the HESS, which consists of two objectives: (1) To minimize the overall losses of the HESS; and (2) to keep the SOC of SCES fluctuating around 50%, considering the HESS needs to smooth the sudden changes in PV power.A simple prediction model is introduced to determine the weights of two sub-objective functions.Also, an SOC correction strategy is proposed to correct the SOCs of the HESS, when power fluctuations do not exceed the allowable limits.(3) According to 100 active power curves of a 750kW PV plant, numerous simulations are carried out to analyze the performances of the proposed smoothing strategies and the smoothing method of using the HESS [15].Simulation results indicate that the performances of the proposed method are much better than the smoothing method of using the HESS.
This study contributes to the field of large-scale grid-connected PV applications.It is beyond the scope of this study to examine the allocation method of the HESS capacities.This paper is structured as follows: Section 2 presents a general description of the power fluctuations of grid-connected PV plants.The details of the coordinated control strategy for the HESS and PV to smooth PV power fluctuations are discussed in Section 3, including the control strategies of the HESS and the PV power curtailment algorithm.Case studies and simulation results are presented in Section 4. The comprehensive performances of two different smoothing strategies are analyzed and discussed in Section 5. Finally, Section 6 presents the conclusions of this work.

System Structures of a PV Plant
The system structures of a grid-connected PV plant with an HESS are shown in Figure 1.P PV denotes the active power of PV, and C pv denotes the installed capacity of a PV plant.The energy and power capacities of BES are denoted by E BES and P BES , respectively.The energy and power capacities of SCES are denoted by E SCES and P SCES , respectively.The SOCs of the battery and the supercapacitor are denoted by S 1 and S 2 , respectively.The total power of the HESS is denoted by P HESS , which is the sum of the BES power, P 1 , and the SCES power, P 2 .The power injected into the utility grid is denoted by P Grid , and P Grid = P PV + P HESS .In addition, the energy management system (EMS) monitors the real-time operating status of the system and controls the power of PV, BES, and SCES.
Sustainability 2019, 11 FOR PEER REVIEW 3 This study contributes to the field of large-scale grid-connected PV applications.It is beyond the scope of this study to examine the allocation method of the HESS capacities.
This paper is structured as follows: Section 2 presents a general description of the power fluctuations of grid-connected PV plants.The details of the coordinated control strategy for the HESS and PV to smooth PV power fluctuations are discussed in Section 3, including the control strategies of the HESS and the PV power curtailment algorithm.Case studies and simulation results are presented in Section 4. The comprehensive performances of two different smoothing strategies are analyzed and discussed in Section 5. Finally, Section 6 presents the conclusions of this work.

System Structures of a PV Plant
The system structures of a grid-connected PV plant with an HESS are shown in Figure 1.PPV denotes the active power of PV, and Cpv denotes the installed capacity of a PV plant.The energy and power capacities of BES are denoted by EBES and PBES, respectively.The energy and power capacities of SCES are denoted by ESCES and PSCES, respectively.The SOCs of the battery and the supercapacitor are denoted by S1 and S2, respectively.The total power of the HESS is denoted by PHESS, which is the sum of the BES power, P1, and the SCES power, P2.The power injected into the utility grid is denoted by PGrid, and PGrid = PPV + PHESS.In addition, the energy management system (EMS) monitors the real-time operating status of the system and controls the power of PV, BES, and SCES.

Definition of the PV Power Fluctuations
In this study, we mainly discuss the power fluctuation of PV in one minute.As shown in Figure .2, the power fluctuation can be defined as the difference between the maximum and minimum power values measured at the PCC in 1 min.

Definition of the PV Power Fluctuations
In this study, we mainly discuss the power fluctuation of PV in one minute.As shown in Figure 2, the power fluctuation can be defined as the difference between the maximum and minimum power values measured at the PCC in 1 min.
Sustainability 2019, 11 FOR PEER REVIEW 3 This study contributes to the field of large-scale grid-connected PV applications.It is beyond the scope of this study to examine the allocation method of the HESS capacities.
This paper is structured as follows: Section 2 presents a general description of the power fluctuations of grid-connected PV plants.The details of the coordinated control strategy for the HESS and PV to smooth PV power fluctuations are discussed in Section 3, including the control strategies of the HESS and the PV power curtailment algorithm.Case studies and simulation results are presented in Section 4. The comprehensive performances of two different smoothing strategies are analyzed and discussed in Section 5. Finally, Section 6 presents the conclusions of this work.

System Structures of a PV Plant
The system structures of a grid-connected PV plant with an HESS are shown in Figure 1.PPV denotes the active power of PV, and Cpv denotes the installed capacity of a PV plant.The energy and power capacities of BES are denoted by EBES and PBES, respectively.The energy and power capacities of SCES are denoted by ESCES and PSCES, respectively.The SOCs of the battery and the supercapacitor are denoted by S1 and S2, respectively.The total power of the HESS is denoted by PHESS, which is the sum of the BES power, P1, and the SCES power, P2.The power injected into the utility grid is denoted by PGrid, and PGrid = PPV + PHESS.In addition, the energy management system (EMS) monitors the real-time operating status of the system and controls the power of PV, BES, and SCES.

Definition of the PV Power Fluctuations
In this study, we mainly discuss the power fluctuation of PV in one minute.As shown in Figure .2, the power fluctuation can be defined as the difference between the maximum and minimum power values measured at the PCC in 1 min.At time t, P Grid,before (t) denotes the power value at the PCC if no correction action is provided, and P Grid,after (t) denotes the power value at the PCC after providing the correction action.The power fluctuations of PV before and after the suppression are denoted by ∆P Grid,before (t) and ∆P Grid,after (t), respectively.According to the historical data of P Grid , ∆P Grid,before (t) can be calculated as follows: ∆P Grid,before (t) = P max,before (t) − P min,before (t), i ≥ j P min,before (t) − P max,before (t), i < j s.t.
where P max,before (t) and P min,before (t) represent the maximum and minimum values of the set A, respectively; i and j represent the indices of the maximum and minimum values of the set A, respectively; function [x, y] = f max (•) returns the largest element of the input argument, where x and y denote the corresponding value and index of the largest element; function [x, y] = f min (•) returns the value and index of the smallest element of the input argument; A is a set of power values measured at the PCC before the correction action is provided, and the time scale of the set A is 1 min; t 0 is the start time; ∆t is the sampling interval and ∆t = 5s.
In the same way, ∆P Grid,after (t) can be calculated by (2) after the smoothing method is applied.
∆P Grid,after (t) = P max,after (t) − P min,after (t), m ≥ n P min,after (t) where B is a set of power values measured at the PCC after providing the correction action, and the time scale of the set B is 1 min; P max,after (t) and P min,after (t) represent the maximum and minimum values of the set B, respectively; m and n represent the indices of the maximum and minimum values of the set B, respectively.

Requirements of the Power Fluctuations of a Grid-Connected PV Plant
This paper assumes that the maximum fluctuation of power injected into the grid from PV should be less than 10 %C pv /min, because this requirement is adopted by the PREPA to restrict the power fluctuations of PV [5].Moreover, it is also adopted by the State Grid Corporation of China to restrict the fluctuations of wind power [16].

Coordinated Control Strategy for the HESS and PV to Smooth PV Power Fluctuations
First, Section 3.1 describes the detailed flowchart of the proposed smoothing strategy.Then, the method used to control the PV power is discussed in Section 3.2.Finally, we explain the control strategies of the HESS in Section 3.3.

Flowchart of the Proposed Smoothing Strategy
This paper proposes a coordinated control strategy for the HESS and PV to smooth PV power fluctuations.The main ideas of this strategy consist of three aspects: (1) Typically, when the PV power fluctuations do not exceed the maximum allowable power of the HESS, the proposed strategy only utilizes the HESS to smooth these power fluctuations while keeping the PV converter operating in the MPPT mode.Namely, the downward power fluctuations are smoothed by discharging (the HESS power is positive, P HESS > 0), in contrast, the upward power fluctuations are smoothed by charging (the HESS power is negative, P HESS < 0).
(2) However, if the upward power fluctuation exceeds the upper power limit of the HESS, the HESS will be charged using the maximum allowable power.At the same time, the excess power fluctuation will be absorbed by reducing the power generated by the PV converter.
(3) In addition, if the downward power fluctuation exceeds the lower power limit of the HESS, in this case, the HESS will be discharged using the maximum allowable power.Despite this, the power fluctuation of PV may not meet the requirements.However, this situation can be avoided by properly allocating the capacities of the HESS and controlling the SOCs of the HESS, which will be discussed later.
In summary, Figure 3 shows the flowchart of the proposed smoothing strategy, which consists of five key steps: Sustainability 2019, 11 FOR PEER REVIEW 5 First, section 3.1 describes the detailed flowchart of the proposed smoothing strategy.Then, the method used to control the PV power is discussed in section 3.2.Finally, we explain the control strategies of the HESS in section 3.3.

Flowchart of the Proposed Smoothing Strategy
This paper proposes a coordinated control strategy for the HESS and PV to smooth PV power fluctuations.The main ideas of this strategy consist of three aspects: 1) Typically, when the PV power fluctuations do not exceed the maximum allowable power of the HESS, the proposed strategy only utilizes the HESS to smooth these power fluctuations while keeping the PV converter operating in the MPPT mode.Namely, the downward power fluctuations are smoothed by discharging (the HESS power is positive, PHESS > 0), in contrast, the upward power fluctuations are smoothed by charging (the HESS power is negative, PHESS < 0).
2) However, if the upward power fluctuation exceeds the upper power limit of the HESS, the HESS will be charged using the maximum allowable power.At the same time, the excess power fluctuation will be absorbed by reducing the power generated by the PV converter.
3) In addition, if the downward power fluctuation exceeds the lower power limit of the HESS, in this case, the HESS will be discharged using the maximum allowable power.Despite this, the power fluctuation of PV may not meet the requirements.However, this situation can be avoided by properly allocating the capacities of the HESS and controlling the SOCs of the HESS, which will be discussed later.
In summary, Figure 3 shows the flowchart of the proposed smoothing strategy, which consists of five key steps:  The first step, the EMS calculates the real-time power fluctuation of PV, ∆P Grid,before (t), by (1), according to the real-time data, P Grid,before (t), and the historical data of P Grid .
The second step, according to the requirements of PV power fluctuations, the power demand of the HESS P * HESS (t), which is used to smooth the current power fluctuation, can be calculated by (4).
Considering the constraints of the HESS, the maximum allowable charging and discharging power of HESS, which is denoted by P min HESS (t) and P max HESS (t), can be calculated by ( 5) and (6).
where the functions of max{x 1 , x 2 } and min{x 1 , x 2 } return the maximum and minimum values between x 1 and x 2 , respectively; S 1,max , and S 1,min are the upper and lower limits of the SOC of the battery, respectively; S 2,max , and S 2,min are the upper and lower limits of the SOC of the supercapacitor, respectively; in this work, S 1,min = 20%, S 1,max = 90%, S 2,min = 5%, and S 2,max = 100%.The third step, according to the power demand of the HESS P * HESS (t) and the maximum power capacity of the HESS, the EMS needs to determine whether to limit the PV power and calculate the actual power command of the HESS.Therefore, a PV power curtailment algorithm is developed to control the PV power, and the details of this algorithm are discussed in the following section.
In addition, after implementing the PV power curtailment algorithm, the power demand of HESS, P * HESS (t), needs to be modified, and the modified power demand is denoted by P * HESS,M (t), which is calculated and explained in Section 3.2.
The fourth step, the HESS can respond immediately after receiving the real power command.There are two possible situations: (1) If the modified power demand, P * HESS,M (t), is not zero, the P * HESS,M (t) will be allocated to BES and SCES using the multi-objective optimization model, which is discussed in Section 3.3.1.
(2) On the other hand, if P * HESS,M (t) is zero, there is no need for the HESS to output power to smooth the current power fluctuation.Therefore, an SOC correction strategy of the HESS is developed to correct the SOCs of BES and SCES, because their SOCs may approach the lower limits due to the PV power curtailment and operating losses of the HESS.The details of this strategy are discussed in Section 3.3.5.
Finally, the EMS sends the power commands to the PV converter, the BES converter, and the SCES converter.Through the close cooperation between the PV converter and the HESS, the power fluctuation at time t will meet the requirements.After that, the EMS will start the next control loop and repeat the above-described steps.

PV Power Curtailment Algorithm
Reducing the reference power value generated by the MPPT algorithm (e.g., the hill climbing method [17]) can easily restrict the upward power fluctuation of PV.Therefore, the PV power curtailment algorithm proposed in this paper is shown in Figure 4.
Reducing the reference power value generated by the MPPT algorithm (e.g., the hill climbing method [17]) can easily restrict the upward power fluctuation of PV.Therefore, the PV power curtailment algorithm proposed in this paper is shown in Figure.The main control processes of this algorithm are as follows: When the power demand of the HESS used to smooth the upward power fluctuation is greater than the maximum allowable charging power of the HESS, i.e., P * HESS (t) < P min HESS (t), the excess power demand is absorbed by reducing the PV power.The comparator controls the output status of the mode selector, if P * HESS (t) < P min HESS (t), then the comparator (port C) outputs signal "1", thus the mode selector switches from mode II to mode I. Consequently, the actual power of the PV converter is reduced by |P * HESS (t) -P min HESS | compared to the original power, P * PV (t), which is calculated by the MPPT algorithm, such as the hill climbing method [17].Otherwise, there is no need to reduce the PV power, and the actual power of PV is equal to P * PV (t).This section has discussed how to control the PV power, according to the real-time state of the HESS and the power demand of the HESS, P * HESS (t).After the implementation of the PV power curtailment algorithm, the P * HESS (t) needs to be modified due to the constraints of the HESS.Therefore, the modified power demand of HESS is denoted by P * HESS,M (t), which can be calculated as follows: 1) If P * HESS (t) < P min HESS (t), the original power demand will be reduced to P min HESS (t) by limiting the power generated by PV.
2) If P * HESS (t) > P max HESS (t), the original power demand will be reduced to P max HESS (t).
3) If P * HESS (t) does not exceed the allowable limits of the HESS at time t, the power demand will not be changed.Having determined the total power that BES and SCES need to output, this paper will now move on to discuss the control strategy of the HESS.

Control Strategy of the HESS
The control strategy of the HESS consists of two keys aims: 1) To optimize the power between BES and SCES; and 2) to correct the SOCs of the battery and the supercapacitor.The specific explanations are as follows: First, if P * HESS,M (t) is not equal to zero, it indicates that the HESS is required to absorb the excess power fluctuation by charging or discharging.Therefore, this paper establishes a multi-objective optimization model to optimize the power of the HESS, considering the overall losses of the HESS and the SOC optimization of SCES.The main control processes of this algorithm are as follows: When the power demand of the HESS used to smooth the upward power fluctuation is greater than the maximum allowable charging power of the HESS, i.e., P * HESS (t) < P min HESS (t), the excess power demand is absorbed by reducing the PV power.The comparator controls the output status of the mode selector, if P * HESS (t) < P min HESS (t), then the comparator (port C) outputs signal "1", thus the mode selector switches from mode II to mode I. Consequently, the actual power of the PV converter is reduced by |P * HESS (t) -P min HESS | compared to the original power, P * PV (t), which is calculated by the MPPT algorithm, such as the hill climbing method [17].Otherwise, there is no need to reduce the PV power, and the actual power of PV is equal to P * PV (t).This section has discussed how to control the PV power, according to the real-time state of the HESS and the power demand of the HESS, P * HESS (t).After the implementation of the PV power curtailment algorithm, the P * HESS (t) needs to be modified due to the constraints of the HESS.Therefore, the modified power demand of HESS is denoted by P * HESS,M (t), which can be calculated as follows: (1) If P * HESS (t) < P min HESS (t), the original power demand will be reduced to P min HESS (t) by limiting the power generated by PV.
(2) If P * HESS (t) > P max HESS (t), the original power demand will be reduced to P max HESS (t).(3) If P * HESS (t) does not exceed the allowable limits of the HESS at time t, the power demand will not be changed.
Having determined the total power that BES and SCES need to output, this paper will now move on to discuss the control strategy of the HESS.

Control Strategy of the HESS
The control strategy of the HESS consists of two keys aims: (1) To optimize the power between BES and SCES; and (2) to correct the SOCs of the battery and the supercapacitor.The specific explanations are as follows: First, if P * HESS,M (t) is not equal to zero, it indicates that the HESS is required to absorb the excess power fluctuation by charging or discharging.Therefore, this paper establishes a multi-objective optimization model to optimize the power of the HESS, considering the overall losses of the HESS and the SOC optimization of SCES.
On the other hand, the HESS does not need to smooth PV power fluctuations when P * HESS,M (t) is equal to zero.Therefore, this paper proposes an SOC correction strategy of HESS to actively control the SOCs of BES and SCES.Because of the PV power curtailment, the total energy absorbed by the HESS may be less than the total energy released.As a result of this, the SOCs of the HESS will gradually decrease and even exceed their lower limits.That has bad effects on smoothing the following power fluctuations.Therefore, a method is needed to improve the SOCs of BES and SCES, once their SOCs are below certain limits.
In the section that follows, first, we carefully describe the objective functions, constraints, and solution algorithms of the multi-objective optimization model of the HESS, which is used to optimize the power between BES and SCES.Then, we explain the details of the SOC correction strategy of the HESS.

Multi-Objective Optimization Model of the HESS
In the problem of dispatching the power of the HESS, this paper mainly considers the losses of the HESS and the SOC optimization of SCES.
First, the operating losses of two energy storage systems are different, because the comprehensive efficiencies of BES and SCES are different.The efficiency of BES (Li-ion battery) is about 85%-90% [14], and the efficiency of SCES is about 90%-95% [14].Therefore, to minimize the overall losses of systems, the efficiencies of the HESS should be considered.
On the other hand, the energy capacity of SCES is usually configured to be smaller, and the power capacity of SCES is usually configured to be larger because of the high costs and high power rates of the supercapacitor [14].Moreover, the lifetime of the supercapacitor is much longer than that of the battery, for example, the number of cycles of the supercapacitor is about 10 to 100 times that of the Li-ion battery [14].Therefore, to fully exploit the characteristics of the supercapacitor, the SCES is encouraged to absorb more components of the PV power fluctuations.However, the SOC of SCES can easily reach the upper or lower limits due to the smaller energy capacity of SCES.Therefore, it is expected that the SOC of SCES can fluctuate around 50%, ensuring that the SCES has sufficient remaining capacity to cope with sudden changes in PV power.
Consequently, the multi-objective optimization model of the HESS is described in (8).The optimization variable is x, which consists of the BES power, P 1 , and the SCES power, P 2 .f 1 (x) and f 2 (x) are sub-objective functions, where f 1 (x) aims to minimize the overall losses of the HESS, and f 2 (x) aims to optimize the SCES SOC.h(x) and g(x) are the equality and inequality constraints of the model, respectively.min Before solving this model, the mathematical models of f 1 (x), f 2 (x), h(x), and g(x) should be presented.

Objective Functions and Constraints
(a) Sub-objective function, f 1 (x).The losses of the HESS consist of operating losses of BES and SCES.To simplify the analysis, this paper uses the charging and discharging power to evaluate the losses of the HESS.Therefore, the sub-objective function, f 1 (x), can be described as follows: where η 1 and η 2 are the efficiencies of BES and SCES, respectively, in this work, η 1 = 90% and η 2 = 95% [18,19].
To minimize the deviation between the SOC of SCES and the reference value of SOC (50%), the sub-objective function, f 2 (x), can be described as follows: where ∆S 2,max represents the maximum difference between the SCES SOC and its reference value (50%).∆S 2,max can be determined according to the upper and lower limits of the SCES SOC, and ∆S 2,max = max{S 2,max -50%, 50% -S 2,min }.In this work, S 2,max = 100%, and S 2,min = 5%, hence ∆S 2,max is 50%.To solve the above-established optimization model, the multi-objective optimization model is transformed into a single-objective optimization model, as shown in (11).Since the value of the sub-objective function, f 2 (x), is [0, 1], no standardization is required.
Objective function: where λ is the weight, and 0 ≤ λ ≤ 1; f 1,min and f 1,max are the minimum and maximum values of the sub-objective function, f 1 , respectively.In addition, the value of weight, λ, is discussed in the following section.Constraints: −P BES ≤ P 1 (t) ≤ P BES −P SCES ≤ P 2 (t) ≤ P SCES ( 13) The constraints consist of four parts: (1) Constraint (12) states that the sum of P 1 (t) and P 2 (t) should equal to P * HESS,M (t) determined by (7).
(2) Constraint (13) requires that the power of BES and SCES cannot exceed their maximum power capacities.
(3) The SOC constraints of the HESS are described in (14).( 4) Constraint (15) indicates that the directions of the power flow of BES and SCES cannot be reversed.For example, BES cannot be charged while SCES is discharging.

Calculation Method of the Weight, λ
In this optimization problem, there is a trade-off between reducing the losses of the HESS and improving the ability of the HESS to cope with the sudden changes in PV power.It is more preferred to use SCES to respond to the power demand of the HESS when the weight, λ, tends to zero, achieving minimum system losses.However, the variation of the SCES SOC will increase, and the mean SOC of SCES will be farther from the reference SOC (50%).Hence, the sudden power fluctuations of PV may not be smoothed completely.On the other hand, the SOC deviation of SCES will decrease when the weight, λ, tends to 1. Hence, the HESS can respond quickly to the sudden power fluctuations of PV.However, the power and frequency of BES will increase, resulting in increased system losses and a reduced service life of BES.
Therefore, this paper dynamically adjusts the value of the weight, λ, using the short-term prediction model, which is used to predict the HESS's energy demand for smoothing the downward power fluctuation, as shown in Figure 5.Note that the energy demand of smoothing the upward power fluctuation is not considered, because the upward power fluctuation can be smoothed by limiting the PV power.An explanation of this control is as follows: actual applications.However, once the HESS's energy demand of smoothing the downward power fluctuation, which is denoted by Eneed(t) and calculated by (16), is greater than 50% of the remaining energy in SCES, then the weight, λ, will be set to 0, otherwise the weight, λ, will always be equal to 0.7.
A short-term prediction mode of Eneed(t) is shown as Figure 5.When PV power suddenly drops at time t, Eneed(t) (i.e., green shaded areas in Figure 5) can be predicted as follows:  ) where Tdelay is the total time that the HESS will be discharged continuously; P max Last (t) represents the maximum power of PV in the past minute; function ceil(x) rounds x to the nearest integer greater than or equal to x; n = Tdelay/Δt.
In addition, the prediction model described in ( 16) is a rather simple prediction strategy.One reason why we did not choose the more advanced tools, such as model predictive control (MPC), is as follows: The advanced tools require a lot of historical data and computation time, which may reduce the real-time performance of the system because this study sets the sampling time to 5 s.Although the prediction method proposed in this paper is simple, the calculation speed is very fast, and the simulation results also proved the effectiveness of the method.Of course, if there is a higher If the PV power is suddenly declined, the HESS needs to be continuously discharged for a while (see Figure 5).However, if the weight, λ, is always fixed, the energy stored in SCES will be quickly emptied, thereby losing the ability to compensate.Actually, in this case, the BES should participate in smoothing that downward power fluctuation as early as possible, rather than starting to discharge when the remaining energy in SCES approaches zero.
Overall, the weight, λ, is usually set to 0.7, which is obtained by the analytic hierarchy process (AHP) method [20].Of course, λ can be flexibly adjusted within the range of 0 to 1 according to the actual applications.However, once the HESS's energy demand of smoothing the downward power fluctuation, which is denoted by E need (t) and calculated by (16), is greater than 50% of the remaining energy in SCES, then the weight, λ, will be set to 0, otherwise the weight, λ, will always be equal to 0.7.
A short-term prediction mode of E need (t) is shown as Figure 5.When PV power suddenly drops at time t, E need (t) (i.e., green shaded areas in Figure 5) can be predicted as follows: where T delay is the total time that the HESS will be discharged continuously; P max Last (t) represents the maximum power of PV in the past minute; function ceil(x) rounds x to the nearest integer greater than or equal to x; n = T delay /∆t.
In addition, the prediction model described in ( 16) is a rather simple prediction strategy.One reason why we did not choose the more advanced tools, such as model predictive control (MPC), is as follows: The advanced tools require a lot of historical data and computation time, which may reduce the real-time performance of the system because this study sets the sampling time to 5 s.Although the prediction method proposed in this paper is simple, the calculation speed is very fast, and the simulation results also proved the effectiveness of the method.Of course, if there is a higher requirement for prediction accuracy, the proposed method may not work.Thus, more advanced prediction tools, such as MPC, can be used.

Solution Method
This paper uses the IBM ILOG CPLEX solver developed by the International Business Machines Corporation (IBM) to solve the multi-objective optimization model [21].The optimization model established in this paper is standard nonlinear programming.Thus, many algorithms can easily solve this problem, such as the artificial intelligence algorithms: Particle Swarm Optimization (PSO) and Genetic Algorithm (GA)- [22].Considering the speed of calculation and the convenience of programming, this paper uses the CPLEX solver to optimize the established model, because it can be called in MATLAB and has good convergence.

SOC Correction Strategy of the HESS
The main idea of the SOC correction strategy of the HESS is that SCES and BES will be charged when the power demands of the HESS used to smooth PV power fluctuations are equal to zero (i.e., P * HESS,M (t) = 0).The specific rules for correcting the SOCs of SCES and BES are as follows: First, the general principle is that SCES will be charged first until its SOC is greater than 30%.Obviously, the charging process is discontinuous, because the HESS SOCs are corrected only when PV power fluctuations do not exceed the limits and the HESS does not need to provide a correction action.Second, if the SOC of SCES is greater than 30%, indicating that SCES does not need to be corrected, then BES will be charged if its SOC is less than 30%.Otherwise, there is no need to correct the BES SOC.Finally, the charging power needs to meet the constraint described in (17), because the current power fluctuation of PV may exceed the limits due to the additional charging power.Also, the charging power needs to meet the constraints described in (13).
where P charge (t) is the maximum charging power that can be used to charge SCES or BES.

Case Studies
In this section, first, three typical scenarios that represent the general cases of PV power fluctuations are established.Then, taking these typical scenarios as an example, the proposed smoothing strategy is verified by numerous simulations.Finally, we analyze the impacts of the PV power curtailment algorithm and the SOC correction strategy on smoothing PV power fluctuations.

Typical Scenarios of PV Power Fluctuations
Based on the database that recorded the active power data of a 750 kW PV plant for one year, these daily active power curves of PV can be divided into three main categories using the k-means clustering algorithm [23]: (1) Scenario I; (2) scenario II; and (3) scenario III, as shown in Figure 6.Due to the limited size of the figure, only three daily power curves were plotted in each scenario.
Machines Corporation (IBM) to solve the multi-objective optimization model [21].The optimization model established in this paper is standard nonlinear programming.Thus, many algorithms can easily solve this problem, such as the artificial intelligence algorithms: Particle Swarm Optimization (PSO) and Genetic Algorithm (GA)- [22].Considering the speed of calculation and the convenience of programming, this paper uses the CPLEX solver to optimize the established model, because it can be called in MATLAB and has good convergence.

SOC Correction Strategy of the HESS
The main idea of the SOC correction strategy of the HESS is that SCES and BES will be charged when the power demands of the HESS used to smooth PV power fluctuations are equal to zero (i.e., P * HESS,M (t) = 0).The specific rules for correcting the SOCs of SCES and BES are as follows: First, the general principle is that SCES will be charged first until its SOC is greater than 30%.Obviously, the charging process is discontinuous, because the HESS SOCs are corrected only when PV power fluctuations do not exceed the limits and the HESS does not need to provide a correction action.Second, if the SOC of SCES is greater than 30%, indicating that SCES does not need to be corrected, then BES will be charged if its SOC is less than 30%.Otherwise, there is no need to correct the BES SOC.Finally, the charging power needs to meet the constraint described in (17), because the current power fluctuation of PV may exceed the limits due to the additional charging power.Also, the charging power needs to meet the constraints described in (13).
( ) where Pcharge(t) is the maximum charging power that can be used to charge SCES or BES.

Case Studies
In this section, first, three typical scenarios that represent the general cases of PV power fluctuations are established.Then, taking these typical scenarios as an example, the proposed smoothing strategy is verified by numerous simulations.Finally, we analyze the impacts of the PV power curtailment algorithm and the SOC correction strategy on smoothing PV power fluctuations.

Typical Scenarios of PV Power Fluctuations
Based on the database that recorded the active power data of a 750 kW PV plant for one year, these daily active power curves of PV can be divided into three main categories using the k-means clustering algorithm [23]: 1) Scenario I; 2) scenario II; and 3) scenario III, as shown in  Among these scenarios, the power fluctuations in scenario I are minimal because it represents the power generated by PV on sunny or rainy days.Scenario III has the worst power fluctuations because the power is generated by PV on cloudy days.Also, the statistical results indicate that scenario II is the most common cases of PV power fluctuations.Therefore, the probability of scenario II is also the highest, accounting for 56% (204 days) of the year.Moreover, the probabilities of scenario I and scenario III are 15% (55 days) and 29% (106 days), respectively.

Simulation and Analysis
Before starting simulations, it is necessary to set simulation parameters, such as the energy and power capacities of BES and SCES, and the initial SOCs of BES and SCES.However, this paper mainly focuses on how to dispatch the power of the HESS after the capacities of the HESS have been determined, and the optimal allocation of the HESS will be discussed carefully in the following research.Therefore, according to previous studies [15,16,24], the simulation parameters were configured as shown in Table 1.The worst scenario of PV power fluctuations (i.e., scenario III) was chosen to verify the correctness and effectiveness of the proposed smoothing strategy.Figure 7 shows the active power injected into the utility grid from the PV plant (i.e., P Grid ), the overall power of the HESS (i.e., P HESS ), and the power generated by the PV converter (i.e., P PV ).Moreover, the fluctuations of the power injected into the grid before and after the suppression can be calculated by ( 1) and ( 2), respectively, as shown in Figure 8.
As can be seen from Figure 7a,b, and Figure 8, after the suppression, the power fluctuations of the grid-connected PV plant are significantly reduced.The results, as shown in Figure 8, indicate that the PV power fluctuations (solid red line) fully meet the requirements described in (3), although the original power fluctuations (solid blue line) are quite serious.The simulation results demonstrate that the proposed smoothing strategy can effectively smooth the power fluctuations of grid-connected PV plants.
It can be seen from Figure 7c,d that the HESS can accurately track the changes in PV power and provide the correct compensation power.Moreover, the PV power curtailment algorithm can effectively reduce the power generated by the PV converter, as shown in Figure 7e,f, when the HESS is unable to fully absorb the upward power fluctuation of PV.Simulation results demonstrate that the proposed method can correctly control the power of the HESS and the PV converter and ensure that the power injected into the utility grid from PV plants can meet the requirements.
Figure 9 shows the power and SOC profiles of the HESS.As can be seen from Figure 9a,c, the power provided by SCES is significantly larger than that of BES because the objective of reducing the losses of the HESS is considered, while allocating the total power demand, P * HESS,M .It can be seen from Figure 9d that the SOC of SCES fluctuates around the reference value of 50%.Although the amplitude of the fluctuation is large, it is still within the allowable range because we also consider the SCES SOC optimization, while optimizing the power demand of the HESS.In summary, these results indicate that the established multi-objective optimization model can correctly dispatch the power of the HESS.
In addition, the SOCs of BES and SCES need to be adjusted to the initial states, i.e., S 1 = 90% and S 2 = 50%, during the night (00:00-5:00), because the HESS needs to smooth the power fluctuations of the next day.Although all the simulations are conducted on data obtained from a 750 kW PV plant, the idea of the proposed smoothing method can also be used in various PV plants.However, it is necessary to properly allocate the energy and power capacities of the HESS, according to the capacity of PV plants and the requirements of power fluctuations, and these studies are ongoing.Figure 9 shows the power and SOC profiles of the HESS.As can be seen from Figure 9(a) and Figure 9(c), the power provided by SCES is significantly larger than that of BES because the objective of reducing the losses of the HESS is considered, while allocating the total power demand, Figure 9 shows the power and SOC profiles of the HESS.As can be seen from Figure 9(a) and Figure 9(c), the power provided by SCES is significantly larger than that of BES because the objective of reducing the losses of the HESS is considered, while allocating the total power demand,   In addition, the SOCs of BES and SCES need to be adjusted to the initial states, i.e., S1 = 90% and S2 = 50%, during the night (00:00-5:00), because the HESS needs to smooth the power fluctuations of the next day.Although all the simulations are conducted on data obtained from a 750 kW PV plant, the idea of the proposed smoothing method can also be used in various PV plants.However, it is necessary to properly allocate the energy and power capacities of the HESS, according to the capacity of PV plants and the requirements of power fluctuations, and these studies are ongoing.

Impacts of the PV Power Curtailment Algorithm on Smoothing Power Fluctuations.
To verify the effectiveness of the PV power curtailment algorithm, another simulation was done based on the same case studied in section 4.2.1.This simulation did not use the PV power curtailment algorithm.However, only the HESS was used to smooth power fluctuations.As a result, the power fluctuation profiles are shown in Figure 10.As can be seen from Figure 10, some fluctuations highlighted by yellow circles exceed the upper limit.In contrast, the power fluctuations shown in Figure 8 are completely within the allowable limits.Therefore, the PV power curtailment algorithm helps the HESS to achieve better smoothing results.Of course, we can also achieve the same results by increasing the capacities of the HESS, but it will significantly increase the initial investment costs of the HESS.

Impacts of the PV Power Curtailment Algorithm on Smoothing Power Fluctuations
To verify the effectiveness of the PV power curtailment algorithm, another simulation was done based on the same case studied in Section 4.2.1.This simulation did not use the PV power curtailment algorithm.However, only the HESS was used to smooth power fluctuations.As a result, the power fluctuation profiles are shown in Figure 10.In addition, the SOCs of BES and SCES need to be adjusted to the initial states, i.e., S1 = 90% and S2 = 50%, during the night (00:00-5:00), because the HESS needs to smooth the power fluctuations of the next day.Although all the simulations are conducted on data obtained from a 750 kW PV plant, the idea of the proposed smoothing method can also be used in various PV plants.However, it is necessary to properly allocate the energy and power capacities of the HESS, according to the capacity of PV plants and the requirements of power fluctuations, and these studies are ongoing.

Impacts of the PV Power Curtailment Algorithm on Smoothing Power Fluctuations.
To verify the effectiveness of the PV power curtailment algorithm, another simulation was done based on the same case studied in section 4.2.1.This simulation did not use the PV power curtailment algorithm.However, only the HESS was used to smooth power fluctuations.As a result, the power fluctuation profiles are shown in Figure 10.As can be seen from Figure 10, some fluctuations highlighted by yellow circles exceed the upper limit.In contrast, the power fluctuations shown in Figure 8 are completely within the allowable limits.Therefore, the PV power curtailment algorithm helps the HESS to achieve better smoothing results.Of course, we can also achieve the same results by increasing the capacities of the HESS, but it will significantly increase the initial investment costs of the HESS.As can be seen from Figure 10, some fluctuations highlighted by yellow circles exceed the upper limit.In contrast, the power fluctuations shown in Figure 8 are completely within the allowable limits.Therefore, the PV power curtailment algorithm helps the HESS to achieve better smoothing results.Of course, we can also achieve the same results by increasing the capacities of the HESS, but it will significantly increase the initial investment costs of the HESS.

Impacts of the SOC Correction Strategy of the HESS on Smoothing Power Fluctuations
This section also used the case discussed in Section 4.2.1, but this simulation did not implement the SOC correction strategy of the HESS.As a result, the power fluctuation profiles are shown in Figure 11, and the SOC profiles are shown in Figure 12  This section also used the case discussed in section 4.2.1, but this simulation did not implement the SOC correction strategy of the HESS.As a result, the power fluctuation profiles are shown in Figure 11, and the SOC profiles are shown in Figure 12.As can be seen from Figure 11, some power fluctuations, which are highlighted by yellow circles, exceed the lower limit.Because, at those moments, the SOCs of SCES reach the lower limit, the HESS can no longer be discharged to smooth those power fluctuations.However, if the SOC correction strategy is applied, SCES and BES will be charged if their SOC is below 30%.As shown in Figure 12(b), the SOC of SCES (solid red line) is corrected several times, such as around 9:40.Overall, it can be seen from Figure 12 that the corrected SOC profiles of BES and SCES (solid red lines) are higher sometimes compared to the uncorrected SOC profiles (solid blue lines).
In summary, these results indicate that the proposed SOC correction strategy can effectively correct the SOCs of the HESS, and it is advantageous for the HESS to smooth PV power fluctuations.

Comparison of Two Smoothing Strategies
To evaluate the advantages and disadvantages of the proposed strategy, this paper compares it with a strategy developed in [15].The reference strategy also uses the HESS to smooth PV power fluctuations, but it does not control the power of PV, but keeps the PV converter operating in the MPPT mode.
In the simulations discussed in section 4.2, the energy and power capacities of the HESS were set to 17 kWh/400 kW (see Table 1).Therefore, the simulation parameters of the reference strategy are also set to be the same as the proposed strategy, as shown in Table 2.This section also used the case discussed in section 4.2.1, but this simulation did not implement the SOC correction strategy of the HESS.As a result, the power fluctuation profiles are shown in Figure 11, and the SOC profiles are shown in Figure 12.As can be seen from Figure 11, some power fluctuations, which are highlighted by yellow circles, exceed the lower limit.Because, at those moments, the SOCs of SCES reach the lower limit, the HESS can no longer be discharged to smooth those power fluctuations.However, if the SOC correction strategy is applied, SCES and BES will be charged if their SOC is below 30%.As shown in Figure 12(b), the SOC of SCES (solid red line) is corrected several times, such as around 9:40.Overall, it can be seen from Figure 12 that the corrected SOC profiles of BES and SCES (solid red lines) are higher sometimes compared to the uncorrected SOC profiles (solid blue lines).
In summary, these results indicate that the proposed SOC correction strategy can effectively correct the SOCs of the HESS, and it is advantageous for the HESS to smooth PV power fluctuations.

Comparison of Two Smoothing Strategies
To evaluate the advantages and disadvantages of the proposed strategy, this paper compares it with a strategy developed in [15].The reference strategy also uses the HESS to smooth PV power fluctuations, but it does not control the power of PV, but keeps the PV converter operating in the MPPT mode.
In the simulations discussed in section 4.2, the energy and power capacities of the HESS were set to 17 kWh/400 kW (see Table 1).Therefore, the simulation parameters of the reference strategy are also set to be the same as the proposed strategy, as shown in Table 2.As can be seen from Figure 11, some power fluctuations, which are highlighted by yellow circles, exceed the lower limit.Because, at those moments, the SOCs of SCES reach the lower limit, the HESS can no longer be discharged to smooth those power fluctuations.However, if the SOC correction strategy is applied, SCES and BES will be charged if their SOC is below 30%.As shown in Figure 12b, the SOC of SCES (solid red line) is corrected several times, such as around 9:40.Overall, it can be seen from Figure 12 that the corrected SOC profiles of BES and SCES (solid red lines) are higher sometimes compared to the uncorrected SOC profiles (solid blue lines).
In summary, these results indicate that the proposed SOC correction strategy can effectively correct the SOCs of the HESS, and it is advantageous for the HESS to smooth PV power fluctuations.

Comparison of Two Smoothing Strategies
To evaluate the advantages and disadvantages of the proposed strategy, this paper compares it with a strategy developed in [15].The reference strategy also uses the HESS to smooth PV power fluctuations, but it does not control the power of PV, but keeps the PV converter operating in the MPPT mode.
In the simulations discussed in Section 4.2, the energy and power capacities of the HESS were set to 17 kWh/400 kW (see Table 1).Therefore, the simulation parameters of the reference strategy are also set to be the same as the proposed strategy, as shown in Table 2.

Performances of Smoothing Power Fluctuations
To evaluate the performances of the above two strategies in various scenarios, 100 PV power curves were selected from Figure 6, including 15 curves selected from scenario I, 56 curves selected from scenario II, and 29 curves selected from scenario III.
After the suppression, the excess power fluctuation is denoted by P over , as shown in (18).Therefore, the daily excess energy, which is denoted by E over , can be calculated by (19).After the suppression, the smaller the value of E over , the better the performance of the smoothing power fluctuations.

Performances of Smoothing Power Fluctuations
To evaluate the performances of the above two strategies in various scenarios, 100 PV power curves were selected from Figure 6, including 15 curves selected from scenario I, 56 curves selected from scenario II, and 29 curves selected from scenario III.
After the suppression, the excess power fluctuation is denoted by Pover, as shown in (18).Therefore, the daily excess energy, which is denoted by Eover, can be calculated by (19).After the suppression, the smaller the value of Eover, the better the performance of the smoothing power fluctuations.
Grid,after limit over Grid,after limit Grid,after limit 0 ( ) ( )= ( ) Therefore, when two smoothing strategies are applied to these 100 cases, the profiles of Eover are shown in Figure 13.The results demonstrate that the performance of the proposed smoothing strategy is better because the daily excess energy (Eover) is significantly less than that of the reference strategy.In fact, when Eover is equal to zero, it means that the power fluctuations after the suppression are completely satisfactory.However, in some cases, Eover is not equal to zero.One explanation of these results is that the HESS capacities were determined based on the case discussed in section 4.2, and for other cases, the capacities may be a little smaller.Note that the capacity allocation of the HESS is beyond the scope of this paper.
In conclusion, under the same conditions (such as capacities of the HESS, requirements of power fluctuations, and the installed capacity of the PV plant), combining the HESS and the PV power curtailment has better performances of smoothing PV power fluctuations than the method of using the HESS.

Lifetime Aging of Batteries and Supercapacitors.
To analyze the lifetime aging of batteries and supercapacitors, the 100% depth of the discharge (DOD) equivalent model developed in [25] was used to evaluate the lifetime aging of batteries, as shown in (20).Moreover, the energy throughput model defined by (21) was used to evaluate the lifetime aging of the supercapacitor.

( )
bat,eq 1, 1, -1 2 0.5 The results demonstrate that the performance of the proposed smoothing strategy is better because the daily excess energy (E over ) is significantly less than that of the reference strategy.In fact, when E over is equal to zero, it means that the power fluctuations after the suppression are completely satisfactory.However, in some cases, E over is not equal to zero.One explanation of these results is that the HESS capacities were determined based on the case discussed in Section 4.2, and for other cases, the capacities may be a little smaller.Note that the capacity allocation of the HESS is beyond the scope of this paper.
In conclusion, under the same conditions (such as capacities of the HESS, requirements of power fluctuations, and the installed capacity of the PV plant), combining the HESS and the PV power curtailment has better performances of smoothing PV power fluctuations than the method of using the HESS.

Lifetime Aging of Batteries and Supercapacitors.
To analyze the lifetime aging of batteries and supercapacitors, the 100% depth of the discharge (DOD) equivalent model developed in [25] was used to evaluate the lifetime aging of batteries, as shown in (20).Moreover, the energy throughput model defined by (21) was used to evaluate the lifetime aging of the supercapacitor.(20) where n bat,eq is the 100% DOD equivalent cycles of a battery; m represents the number of extreme points of the battery's SOC curve, i.e., S 1 (t), t∈T; S 1,i is the SOC value at the extreme point i; k p is a constant, 0.8 ≤ k p ≤ 2.1, in this work, k p = 0.976 [25,26].
n cap,eq = ∆t 24:00 ∑ t=00:00 where n cap,eq is the equivalent cycles of the supercapacitor.Therefore, when two smoothing strategies are applied to these 100 cases, the profiles of n bat,eq , and n cap,eq are shown in Figure 14.
Sustainability 2019, 11 FOR PEER REVIEW 17 where nbat,eq is the 100% DOD equivalent cycles of a battery; m represents the number of extreme points of the battery's SOC curve, i.e., S1(t), t∈T; S1,i is the SOC value at the extreme point i; kp is a constant, 0.8 ≤ kp ≤ 2.1, in this work, kp = 0.976 [25,26].
( ) 24:00 cap,eq 2 SCES 00:00 where ncap,eq is the equivalent cycles of the supercapacitor.Therefore, when two smoothing strategies are applied to these 100 cases, the profiles of nbat,eq, and ncap,eq are shown in Figure .14.As can be seen from Figure 14, it is apparent that the battery cycles of the proposed smoothing strategy are smaller than that of the reference strategy.Therefore, these results indicate that the PV power curtailment can reduce the number of cycles of the battery while using the HESS to smooth PV power fluctuations.
In conclusion, compared with using the HESS to smooth PV power fluctuations (reference strategy), the proposed smoothing strategy can significantly reduce the number of cycles of the battery, so that the actual service life of BES can be prolonged.

Operating Losses of the HESS.
Figure 15 shows the profiles of the operating losses of the HESS when two smoothing strategies are applied to these 100 cases.The simulation results indicate that the proposed smoothing strategy has lower losses because the number of cycles of the battery is reduced.
In conclusion, the operating losses of the HESS can be reduced when using the HESS and the PV power curtailment to smooth PV power fluctuations.

Power Generation of the PV Plant.
To analyze the energy losses caused by limiting PV power, the daily power generation of these 100 cases was calculated, as shown in Figure 16.Before the suppression, the daily power generation of the PV plant is described in Figure 16  As can be seen from Figure 14, it is apparent that the battery cycles of the proposed smoothing strategy are smaller than that of the reference strategy.Therefore, these results indicate that the PV power curtailment can reduce the number of cycles of the battery while using the HESS to smooth PV power fluctuations.
In conclusion, compared with using the HESS to smooth PV power fluctuations (reference strategy), the proposed smoothing strategy can significantly reduce the number of cycles of the battery, so that the actual service life of BES can be prolonged.In conclusion, the operating losses of the HESS can be reduced when using the HESS and the PV power curtailment to smooth PV power fluctuations.

Power Generation of the PV Plant
To analyze the energy losses caused by limiting PV power, the daily power generation of these 100 cases was calculated, as shown in Figure 16.Before the suppression, the daily power generation of the PV plant is described in Figure 16a; after smoothing by the two strategies, the values of the daily power generation of the PV plant were converted into the per unit (p.u.) values for comparison, as shown in Figure 16b.In addition, based on the power generation data obtained from these 100 cases, the annual power generation of this PV plant was estimated.If no correction action is provided, the annual power generation of this PV plant is about 325,863 kWh.The annual power generation of the reference strategy is reduced by 0.32% (324,840.81kWh) because of the operating losses of the HESS, and the proposed strategy is reduced by 0.44% (324,424.20kWh) because of the HESS losses and the PV power curtailment.
In conclusion, for a 750 kW PV plant, when using the HESS and the PV power curtailment to smooth PV power fluctuations, the annual power generation of the PV plant will be reduced by about 0.44%, which is slightly greater than the method of using the HESS.

Lifetime Net Profits of the PV Plant
The lifetime net profit of the PV plant is an important indicator to evaluate the economics of the different smoothing strategies.Therefore, this paper established a lifetime net profit model, which is shown as follows: where RNP(y) represents the cumulative net profits of the PV plant in the y th year (y ≤ 25); i is the index of the year (i ≤ y); Rpv, Com, Cre, and Cfine represent the revenue of selling PV electricity, the operation and maintenance costs of the energy storage system, the replacement costs of the energy storage system, and the fine caused by the power fluctuations that exceed the allowable limits, respectively; Cinv is the initial investment costs of the PV plant and energy storage systems; r is the discount rate, r = 2% [27].Rpv, Com, Cre, Cfine, and Cinv are calculated as follows: In addition, based on the power generation data obtained from these 100 cases, the annual power generation of this PV plant was estimated.If no correction action is provided, the annual power generation of this PV plant is about 325,863 kWh.The annual power generation of the reference strategy is reduced by 0.32% (324,840.81kWh) because of the operating losses of the HESS, and the proposed strategy is reduced by 0.44% (324,424.20kWh) because of the HESS losses and the PV power curtailment.
In conclusion, for a 750 kW PV plant, when using the HESS and the PV power curtailment to smooth PV power fluctuations, the annual power generation of the PV plant will be reduced by about 0.44%, which is slightly greater than the method of using the HESS.

Lifetime Net Profits of the PV Plant
The lifetime net profit of the PV plant is an important indicator to evaluate the economics of the different smoothing strategies.Therefore, this paper established a lifetime net profit model, which is shown as follows: where R NP (y) represents the cumulative net profits of the PV plant in the y th year (y ≤ 25); i is the index of the year (i ≤ y); R pv , C om , C re , and C fine represent the revenue of selling PV electricity, the operation and maintenance costs of the energy storage system, the replacement costs of the energy storage system, and the fine caused by the power fluctuations that exceed the allowable limits, respectively; C inv is the initial investment costs of the PV plant and energy storage systems; r is the discount rate, r = 2% [27].R pv , C om , C re , C fine , and C inv are calculated as follows: R pv = m 1 E PG (23) where E PG denotes the annual power generation of the PV plant.Referring to the report released by the US Energy Information Administration, the feed-in price of solar power, m 1 , is set to US$0.13/kWh [28].
where a is a conversion factor for calculating the operation and maintenance costs of the system, in this paper, a = 0.5% [27].
where E over,d denotes the excess energy calculated by (19) on the d th day; m 6 is the unit price of the excess energy, this paper sets m 6 = 10m 1 (i.e., US$1.3/kWh).Moreover, the other parameters of the lifetime net profit model are shown in Table 3. Overall, Figure 17 shows the profiles of the lifetime net profits of the PV plant when the two smoothing strategies are applied.As can be seen from Figure 17, the simulation results demonstrate that the proposed smoothing strategy is more economical than the reference strategy.= m R E (23) where EPG denotes the annual power generation of the PV plant.Referring to the report released by the US Energy Information Administration, the feed-in price of solar power, m1, is set to US$0.13/kWh [28].where m2, m3, m4, and m5 are the unit price of batteries, supercapacitors, energy storage converters, and PV systems, respectively.
where a is a conversion factor for calculating the operation and maintenance costs of the system, in this paper, a = 0.5% [27].where λ1 and λ2 are 0-1 variables, when the battery's (supercapacitor's) cumulative cycles in the i th year reach the upper limit, then λ1 = 1 (λ2 = 1), otherwise, λ1 = 0 (λ2 = 0).E m (27) where Eover,d denotes the excess energy calculated by (19) on the d th day; m6 is the unit price of the excess energy, this paper sets m6 = 10m1 (i.e., US$1.3/kWh).Moreover, the other parameters of the lifetime net profit model are shown in Table 3. Overall, Figure 17 shows the profiles of the lifetime net profits of the PV plant when the two smoothing strategies are applied.As can be seen from Figure 17, the simulation results demonstrate that the proposed smoothing strategy is more economical than the reference strategy.
In conclusion, combining the HESS and the PV power curtailment is a more economical method of smoothing PV power fluctuations compared to the method of using the HESS, although it will result in more energy losses due to the PV power curtailment.In conclusion, combining the HESS and the PV power curtailment is a more economical method of smoothing PV power fluctuations compared to the method of using the HESS, although it will result in more energy losses due to the PV power curtailment.

Conclusions
In this study, a coordinated control strategy for the HESS and PV to smooth the power fluctuations of PV was developed.The main conclusions are as follows: (1) A PV power curtailment algorithm was developed to automatically limit the power generated by the PV converter when the HESS could not fully absorb upward power fluctuations.A multi-objective optimization model of the HESS was established to minimize the total losses of the HESS and to optimize the SOC of the supercapacitor while dispatching the power demand of the HESS.To prevent the SOC of the supercapacitor from reaching the lower limit of its constraint during the period of smoothing downward power fluctuations, a simple short-term prediction model was presented to dynamically adjust the weights of the sub-objective functions.Also, the HESS SOCs may approach their lower limits due to the curtailment of PV power and energy losses.Therefore, an SOC correction strategy of the HESS was proposed to correct the SOCs of the battery and the supercapacitor, when PV power fluctuations did not exceed the allowable range.Simulation results proved the effectiveness and correctness of the proposed strategies.
(2) Compared to the method of using the HESS to smooth PV power fluctuations, the advantages of the proposed smoothing method are as follows: • The proposed method is more economical and can increase the net profits of the PV and HESS plant.

•
The proposed method can reduce not only the cycles of batteries, but also the energy losses of the HESS.

•
Under the same conditions (such as the capacities of the HESS, requirements of the power fluctuations, and the installed capacity of the PV plant), the proposed method has better performances of smoothing PV power fluctuations than the method of using the HESS.
However, the disadvantage of the proposed smoothing strategy is that limiting the PV power will reduce the power generation of PV plants.For example, the annual power generation of a 750 kWp PV plant will be reduced by 0.44% (about 1439 kWh).
It is recommended that further research be undertaken in the following aspects: (1) An optimal allocation method of the HESS is needed to determine the energy and power capacities of BES and SCES; and (2) it would be interesting to examine the use of flywheel energy storage systems (FESS) [32] as an alternative to the supercapacitor in the HESS, because the advantages of the FESS are a rapid response, low hazardous waste, high power, and long-life compared to batteries.

Figure 1 .
Figure 1.System structures of a grid-connected PV `plant with an HESS.

Figure 2 .
Figure 2. PV power fluctuations in one minute.

Figure 1 .
Figure 1.System structures of a grid-connected PV 'plant with an HESS.

Figure 1 .
Figure 1.System structures of a grid-connected PV `plant with an HESS.

Figure 2 .
Figure 2. PV power fluctuations in one minute.Figure 2. PV power fluctuations in one minute.

Figure 2 .
Figure 2. PV power fluctuations in one minute.Figure 2. PV power fluctuations in one minute.

Figure 3 .
Figure 3. Coordinated control strategy for the HESS and PV to smooth PV power fluctuations.Figure 3. Coordinated control strategy for the HESS and PV to smooth PV power fluctuations.

Figure 3 .
Figure 3. Coordinated control strategy for the HESS and PV to smooth PV power fluctuations.Figure 3. Coordinated control strategy for the HESS and PV to smooth PV power fluctuations. 4.

Figure 4 .
Figure 4. Schematic diagram of the PV power curtailment algorithm.

Figure 4 .
Figure 4. Schematic diagram of the PV power curtailment algorithm.

Figure 5 .
Figure 5.A short-term prediction model is used to predict the energy demand of the HESS.

Figure 5 .
Figure 5.A short-term prediction model is used to predict the energy demand of the HESS.
Figure. 6. Due to the limited size of the figure, only three daily power curves were plotted in each scenario.

Figure 6 .
Figure 6.Typical scenarios of PV power fluctuations: (a) Scenario I; (b) scenario II; and (c) scenario III.

Figure 7 .Figure 8 .
Figure 7. Active power profiles of the PV plant and the HESS: (a) the power injected into the grid before and after the suppression; (b) the curves from 12:00 to 14:00 in (a) are enlarged; (c) the profile of the overall power of the HESS; (d) the curve from 12:00 to 14:00 in (c) is enlarged; (e) the profiles of the power generated by the PV converter; and (f) the curves from 12:00 to 14:00 in (e) are enlarged.

Figure 7 .Figure 7 .Figure 8 .
Figure 7. Active power profiles of the PV plant and the HESS: (a) the power injected into the grid before and after the suppression; (b) the curves from 12:00 to 14:00 in (a) are enlarged; (c) the profile of the overall power of the HESS; (d) the curve from 12:00 to 14:00 in (c) is enlarged; (e) the profiles of the power generated by the PV converter; and (f) the curves from 12:00 to 14:00 in (e) are enlarged.

Figure 8 .
Figure 8. Power fluctuation profiles before and after the suppression.
. It can be seen from Figure9(d) that the SOC of SCES fluctuates around the reference value of 50%.Although the amplitude of the fluctuation is large, it is still within the allowable range because we also consider the SCES SOC optimization, while optimizing the power demand of the HESS.In summary, these results indicate that the established multi-objective optimization model can correctly dispatch the power of the HESS.

Figure 9 .
Figure 9. Power and SOC profiles of the HESS: (a) power of BES; (b) SOC of BES; (c) power of SCES; and (d) SOC of SCES.

Figure 10 .
Figure 10.Power fluctuation profiles before and after the suppression when the PV power curtailment algorithm was not applied.

Figure 9 .
Figure 9. Power and SOC profiles of the HESS: (a) power of BES; (b) SOC of BES; (c) power of SCES; and (d) SOC of SCES.

Figure 9 .
Figure 9. Power and SOC profiles of the HESS: (a) power of BES; (b) SOC of BES; (c) power of SCES; and (d) SOC of SCES.

Figure 10 .
Figure 10.Power fluctuation profiles before and after the suppression when the PV power curtailment algorithm was not applied.

Figure 10 .
Figure 10.Power fluctuation profiles before and after the suppression when the PV power curtailment algorithm was not applied. .

4. 2 . 3 .
Impacts of the SOC Correction Strategy of the HESS on Smoothing Power Fluctuations.

Figure 11 .
Figure 11.Power fluctuation profiles before and after the suppression when the SOC correction strategy of the HESS was not applied.

Figure 12 .
Figure 12.SOC profiles of BES and SCES when the SOC correction strategy of the HESS was not applied: (a) SOC curves of BES; and (b) SOC curves of SCES.

Figure 11 .
Figure 11.Power fluctuation profiles before and after the suppression when the SOC correction strategy of the HESS was not applied.

Figure 11 .
Figure 11.Power fluctuation profiles before and after the suppression when the SOC correction strategy of the HESS was not applied.

Figure 12 .
Figure 12.SOC profiles of BES and SCES when the SOC correction strategy of the HESS was not applied: (a) SOC curves of BES; and (b) SOC curves of SCES.

Figure 12 .
Figure 12.SOC profiles of BES and SCES when the SOC correction strategy of the HESS was not applied: (a) SOC curves of BES; and (b) SOC curves of SCES.

Figure 13 .
Figure 13.Daily excess energy (Eover) profiles of the two smoothing strategies in 100 cases.

Figure 13 .
Figure 13.Daily excess energy (E over ) profiles of the two smoothing strategies in 100 cases.

Figure 14 .
Figure 14.The profiles of nbat,eq and ncap,eq when two smoothing strategies are applied to these 100 cases: (a) the number of cycles of the battery, nbat,eq; and (b) the number of cycles of the supercapacitor, ncap,eq.
(a); after smoothing by the two strategies, the values of the daily power generation of the PV plant were converted into the per unit (p.u.) values for comparison, as shown in Figure 16(b).

Figure 14 .
Figure 14.The profiles of n bat,eq and n cap,eq when two smoothing strategies are applied to these 100 cases: (a) the number of cycles of the battery, n bat,eq ; and (b) the number of cycles of the supercapacitor, n cap,eq .

Figure 15 Figure 15 .
Figure 15 shows the profiles of the operating losses of the HESS when two smoothing strategies are applied to these 100 cases.The simulation results indicate that the proposed smoothing strategy has lower losses because the number of cycles of the battery is reduced.Sustainability 2019, 11 FOR PEER REVIEW 18

Figure 15 .
Figure 15.Profiles of the operating losses of the HESS when two smoothing strategies are applied to these 100 cases.

Figure 15 .
Figure 15.Profiles of the operating losses of the HESS when two smoothing strategies are applied to these 100 cases.

Figure 16 .
Figure 16.Profiles of the power generation of the PV plant when two smoothing strategies are applied to these 100 cases: (a) the daily power generation of the PV plant before the suppression; and (b) the daily power generation of the PV plant (p.u.) under different smoothing strategies.

Figure 16 .
Figure 16.Profiles of the power generation of the PV plant when two smoothing strategies are applied to these 100 cases: (a) the daily power generation of the PV plant before the suppression; and (b) the daily power generation of the PV plant (p.u.) under different smoothing strategies.

Figure 17 .
Figure 17.Profiles of the lifetime net profits of the PV plant when different smoothing strategies are applied.

Table 2 .
Energy and power capacities of the HESS of the two smoothing strategies.

Table 2 .
Energy and power capacities of the HESS of the two smoothing strategies.

Table 2 .
Energy and power capacities of the HESS of the two smoothing strategies.

Table 3 .
Parameters of the lifetime net profit model.

Table 3 .
Parameters of the lifetime net profit model.