Dispatching of Wind/battery Energy Storage Hybrid Systems Using Inner Point Method-based Model Predictive Control

The application of large scale energy storage makes wind farms more dispatchable, which lowers operating risks to the grid from interconnected large scale wind farms. In order to make full use of the flexibility and controllability of energy storage to improve the schedulability of wind farms, this paper presents a rolling and dispatching control strategy with a battery energy storage system (BESS) based on model predictive control (MPC). The proposed control scheme firstly plans expected output, i.e., dispatching order, of a wind/battery energy storage hybrid system based on the predicted output of the wind farm, then calculates the order in the predictive horizon with the receding horizon optimization and the limitations of energy storage such as state of charge and depth of charge/discharge to maintain the combination of active output of the wind farm and the BESS to track dispatching order at the extreme. The paper shows and analyses the effectiveness of the proposed strategy with different sizes of capacity of the BESS based on the actual output of a certain actual wind farm in the northeast of China. The results show that the proposed strategy that controls the BESS could improve the schedulability of the wind farm and maintain smooth output of wind/battery energy storage hybrid system while tracking the dispatching orders. When the capacity of the BESS is 20% or the rated capacity of the wind farm, the mean dispatching error is only 0.153% of the rated capacity of the wind farm.


Introduction
Wind power technology is the most mature and widespread applied one among new energy power generation technologies.Exploration of wind energy on a large scale could substantialy alleviate the problems of energy shortages and environmental pollution across the globe.Up to the end of 2015, the global installed wind capacity had reached 432.4 GW, of which the figure in China was 148 GW, representing the largest share around the world [1].Random output of wind speed lowers the schedulability of active power in a wind farm, which in turn limits the availability of wind energy [2].With increasing permeability of wind power in the state grid, lower schedulability would increase operational risks, especially wind power curtailment [3,4].
Proven large scale battery energy storage system (BESS) technology and its application in the grid provides a new platform to increase the utilization rate of new energy sources [5][6][7][8][9], to improve the schedulability in a wind farm and to guarantee the availability to the grid of high wind power that can operate safely and steadily [10,11].In recent years, domestic and foreign scholars have studied control strategies for BESS that have the ability of charging and discharging power rapidly to smooth fluctuating wind power widely and deeply [12][13][14][15][16][17][18][19][20].So far the most widely applied control algorithm is low pass filtering (LPF) that obtains the wind power target after the actual wind power is sent to the first order filter and then treats the difference between the target and actual wind power as a compensating instruction sent to the BESS.In [16,17], the authors proposed a moving average control algorithm which control principle is similar to that of LPF.In [18], the authors decomposed the wind power with the wavelet transform (WT), and then hybrid energy storage was applied to smooth the corresponding fluctuating component at different frequencies, which alleviated the pressure of fluctuating wind power after the wind plant was integrated into the grid.In [19], the authors proposed a dynamic programming (DP) algorithm with the consideration of state of charge (SOC) of the BESS to smooth fluctuating wind power.In [20], the authors obtained a smoothing target which minimized the action of the BESS, and then Model Predictive Control (MPC) was applied to make the total output smooth.
To sum up, the common objective of the above control algorithms is to smooth the output of wind farms and reduce the deviation between the target and actual wind power to the greatest extent, but no dispatch order is introduced, which limits the contribution of these methods to the schedulability of the wind plant.MPC is wildly applied with good results in the fields of renewable energy generation and control.In [21], current control of a Hybrid Fuel Cell Power system is realized.In [22], frequency/voltage control of a wind-diesel power system based on MPC is presented.In [23], MPC is applied in dispatching active wind power of an off-grid wind turbine/photovoltaic/hydrogen/battery hybrid system.In [24], active power control for a wind/electrolyzer/fuel cell hybrid system based on non-linear MPC is proposed.This paper considers the grid requirements and proposes the MPC for the BESS to improve the schedulability of the wind farm, and then compares the effects of two methods of solving for MPC, the Inner Point (IP) method (IP-MPC) and the Quadratic Programming (QP) method (QP-MPC).Firstly, the wind power dispatch order in the coming dispatching period (30 min) is obtained according to the predictive wind power, and this is regarded as the operational target for the wind/battery energy storage hybrid system, then the control command for the BESS in the dispatching period is acquired through MPC, which can ensure the total output of the wind/battery energy storage hybrid system could be injected into the grid smoothly and safely.
The remainder of the paper is organized as follows.Firstly, the dynamic model of the wind/battery energy storage hybrid system is introduced in Section 2. Second, the dispatching strategy for a wind/battery energy storage hybrid system based on MPC is proposed in Section 3. Then the presented dispatching strategy is investigated by application to a certain wind farm in the northeast of China in Section 4.

System Configuration
The overall structure of the wind/battery energy storage hybrid system is illustrated in Figure 1.The system is consist of a wind farm, a BESS, a converter and a power management unit (PMU), and the system is coupled to the main grid via transmission line.The BESS is connected to the grid at the point of common coupling (PCC) and charges/discharges power via the converter to smooth the wind power injected into the grid.
This paper aims at improving the schedulability of the wind farm, which means to enable the output of the wind/battery storage hybrid system to approach the desired dispatch curve P D as much as possible.In order to achieve the goal, a wind power forecasting model is applied to predict future p-step wind power at each sampling point, and the dispatch order is based on these predicted values.The reasonable action of the BESS (P B ) is calculated by the control strategy in the PMU to compensate for the difference between the dispatch order P D and the actual wind power P W .Under this strategy the power injected to the grid (total output of the wind/battery energy storage hybrid system) P G would follow the dispatch order from the transmission system operator (TSO), and the effects mainly depends on the capacity of the battery and the control algorithm in the PMU.

System Modeling
According to Figure 1 and the analysis in Section 2.1, the model of wind/battery energy storage hybrid system can be formulated as [25]: ( ) ( ) where PG(k) is the total output of the wind/battery energy storage hybrid system; PB(k) is the power control signal for the BESS.When PB(k) > 0, it denotes the BESS discharges power; when PB(k) < 0, it denotes the BESS charges power.PW(k) denotes the actual wind power, EB(k) denotes the remaining energy in the BESS.State of Charge (SOC) of a BESS is the ratio between available energy and its rated capacity, expressed as a percentage.ΔTB is the transfer coefficient from MW to MWh, for instance when the sampling time is 5 min, ΔTB = 5 min/60 min = 1/12.Equation ( 1) can be written in the form of state space as follows: with, state vector: control input: u(k) = PB(k); system output: y(k) = PG(k); combined with the system: The controllability matrix Q = [B1, AB1] is full rank, so the system is controllable.
With consideration of the technology and economy, in actual operation the BESS is limited by the SOC and charging/discharging power:

System Modeling
According to Figure 1 and the analysis in Section 2.1, the model of wind/battery energy storage hybrid system can be formulated as [25]: where P G (k) is the total output of the wind/battery energy storage hybrid system; P B (k) is the power control signal for the BESS.When P B (k) > 0, it denotes the BESS discharges power; when P B (k) < 0, it denotes the BESS charges power.P W (k) denotes the actual wind power, E B (k) denotes the remaining energy in the BESS.State of Charge (SOC) of a BESS is the ratio between available energy and its rated capacity, expressed as a percentage.∆T B is the transfer coefficient from MW to MWh, for instance when the sampling time is 5 min, ∆T B = 5 min/60 min = 1/12.Equation ( 1) can be written in the form of state space as follows: with, state vector: control input: u(k) = P B (k); system output: y(k) = P G (k); combined with the system: Energies 2016, 9, 629 4 of 16 The controllability matrix Q = [B 1 , AB 1 ] is full rank, so the system is controllable.With consideration of the technology and economy, in actual operation the BESS is limited by the SOC and charging/discharging power: The lifetime of a BESS is influenced by the depth of discharge (DoD), and in order to prolong the lifespan of a BESS and reduce costs, there would be a strict constraint on SOC of a BESS at any time k: where, x 2 (k) denotes the remaining energy at time k; C B denotes rate capacity of a BESS; S min denotes lower limit of a BESS, S max denotes upper limit of a BESS.

(B) Charge/Discharge Power Constraint
To ensure the BESS can operate safely, there would be strict constraints on the value of charging/discharging power in the actual BESS system: P c " min ppS max ´S pk ´1qq C B {∆T B , P max_ch q (6) where, P c (k) and P d (k) are the available charging power and discharging power of the BESS, respectively, P max_ch and P max_f are the maximum charging power and discharging power separately.

Wind Power Predicting and Dispatch Curve
This paper aims at improving the schedulability for the wind farm using BESS, i.e., the total active output of the wind/battery energy storage hybrid system (P G ) could trace the dispatch curve.Therefore, a reasonable dispatch order is significant to achieve the goal.
The paper focuses on calculating action order of the BESS, not the method on predicting wind power, therefore the method called ARIMA is adapted to predict wind power [26].The model of ARMA could expressed as: where, ϕ 1 , ϕ 2 , . . ., ϕ p and θ 1 , θ 2 , . . ., θ q are all constants.a t is white noise time series, therefore time series {z t } comply with ARMA with the order of (p, q), described as ARMA(p, q).Wind power is a non-stationary time series that needs to be tackled with difference processing.After being processed d times, {∇ d z t } becomes a stationary ARMA series, therefore: It is ARMA with the order of (p, d, q), described as ARIMA(p, d, q), F is the backward difference operator.The key for predicting wind power with ARMA is to identify the parameters of ϕ 1 , ϕ 2 , ϕ p and θ 1 , θ 2 , θ q .In the paper, least squares estimation is used to identify the parameters which involves minimizing the residual sum of squares.
Energies 2016, 9, 629 With the consideration of computing time and computing accuracy, in the paper the data in the last 12 h with 144 wind power data is sampled to forecast future wind power in 0.5 h with six wind power data.
The dispatch order is made based on wind power forecasting values, as a consequence the precision of prediction in a sense could impact the control effects.The rules that all the generation units submit the dispatch power that the plants have to send into the grid to the TSO in each dispatching interval [9] must be satisfied.Therefore, the dispatch power is designed to meet a constant value in each dispatch interval and is obtained by averaging the predictive wind power in each dispatch interval.In the k th dispatching interval i.e., kT d < t ď (k + 1)T d , the dispatch power is defined as follows: where, T d is the dispatching interval and it is assumed to be 30 min in this paper.The control interval of the BESS is 5 min.The sampling time interval of actual wind power is 5 min and the interval of predictive wind power is also 5 min.The method of acquiring the dispatch order is illustrated in Figure 2.
Energies 2016, 9, 629 5 of 16 constant value in each dispatch interval and is obtained by averaging the predictive wind power in each dispatch interval.In the k th dispatching interval i.e., kTd < t ≤ (k + 1)Td, the dispatch power is defined as follows: where, Td is the dispatching interval and it is assumed to be 30 min in this paper.The control interval of the BESS is 5 min.The sampling time interval of actual wind power is 5 min and the interval of predictive wind power is also 5 min.The method of acquiring the dispatch order is illustrated in Figure 2. At the beginning of each dispatching interval, the dispatch order is obtained by averaging the predicted wind power and is selected as the target for MPC in this paper to enable the total output of the system to track correctly.

Design for MPC Contoller for Dispatching
Model Predictive Control (MPC) is a method to solve online optimal control problems.It can find out present and a series of future control values according to the state x(k) and present and a series of future constraints.The control strategy for MPC proposed in the paper is illustrated in Figure 3 and it is comprised of a predictive model and rolling optimization.At the beginning of each dispatching interval, the dispatch order is obtained by averaging the predicted wind power and is selected as the target for MPC in this paper to enable the total output of the system to track correctly.

Design for MPC Contoller for Dispatching
Model Predictive Control (MPC) is a method to solve online optimal control problems.It can find out present and a series of future control values according to the state x(k) and present and a series of future constraints.The control strategy for MPC proposed in the paper is illustrated in Figure 3 and it is comprised of a predictive model and rolling optimization.

Design for MPC Contoller for Dispatching
Model Predictive Control (MPC) is a method to solve online optimal control problems.It can find out present and a series of future control values according to the state x(k) and present and a series of future constraints.The control strategy for MPC proposed in the paper is illustrated in Figure 3 and it is comprised of a predictive model and rolling optimization.The object of the wind/battery energy storage hybrid system is to enable the total output to track the dispatch order extremely closely with the lowest cost of BESS action, i.e., the charging/discharging power of the BESS.In order to achieve this goal to improve the lifespan of the BESS, the function of control performance of the system could be described as the following: where, (k+i|k) is the predicted value based on the information at time k, m p is the predictive horizon, m c is the control horizon, α P p0, 1q is the balance coefficient of the tracing performance and the control action.In order to simplify the calculations, in this paper the predictive horizon equals the control horizon, i.e., m p = m c = p.Together with (2)-( 7), the MPC formulation for wind/battery storage hybrid systems used to dispatch short-term wind power can be described as (13): where, y ref is reference trajectory, which in this paper is the dispatch order.
Based on the principle of the MPC, the predicted value y(k+i|k) in the predictive horizon could be obtained according to the current system state x i (k+1) at time k and formulation (2), and it could be described as (15): where, y(k+j|k) denotes the predicted system response at k+j being calculated at time k, j = 1, . . ., p; u(k+j|k) denotes the calculated control value at time k, j = 0, 1, . . ., p ´1; I denotes identity matrix.
Energies 2016, 9, 629 7 of 16 Equation ( 15) could be described in matrix form: Ypkq " Gxpkq `KUpkq `LDpkq (16) where: Ypkq " ry pk `1|kq , y pk `2|kq , ¨¨¨, y pk `p|kqs T (17) Upkq " ru pk|kq , u pk `1|kq , ¨¨¨, u pk `p ´1|kqs T (18) Dpkq " rd pk|kq , d pk `1|kq , ¨¨¨, d pk `p ´1|kqs T G " rCA, CA 2 , ¨¨¨, CA p s T ( 20) Furthermore, in control horizon, indicator function of quadratic performance for controller of the combined system could be described as: Equation ( 13) and the constraints ( 14) formulate a quadratic programming problem, so rolling optimization is adapted to obtain the optimal control sequence, only the first of which is sent into the system to work.At time k + 1, the system information updates and the optimal procedure is repeated as above.

Rolling Optimization Using Inner Point Method
In the design for MPC, control sequence (U) is the useful one and the items that have no relationship with U would be neglected, so the control objective would be written as follows: It is standard form of quadratic programming: (25) Energies 2016, 9, 629 8 of 16 The inner point method is adapted in the paper to solve the control value, the basic line of thought is based on the condition of Karush-Kuhn-Tucker (KKT): ´Aineq u `bineq ´s " 0 s i z i " 0, i " 1, 2, . . .m s ě 0; z ě 0 (28) Equation ( 28) is rewritten in the iteration form: where, σ is the diminish iteration variable converging to 0, σ P (0, 1], m is the number of inequality constraints.Therefore, the step of inner point to solve the quadratic programming is as follows: Step 1: Select the initial condition (x(0), s(0), z(0)) satisfying s(0) > 0, z(0) > 0.
Step 2: Obtain the ∆u and ∆z satisfying the following equation: ff and ∆s " ´s `Z´1 pβe ´T∆Zq .
Step 4: Judge the convergence.If it converges, the procedure would stop and the optimal control values are obtained, otherwise, the value are treated as initial value and sent into step 2 to solve iteratively.
The optimal control sequence is solved by IP for the quadratic programming problem, and only the first of which is sent into the BESS.The next time, the procedure is repeated to realize rolling optimization.

Simulation Results
In this paper, actual wind power data from a wind farm in northeast China whose capacity is 45 MW and the sampling time interval is 5 min was adopted.The sodium-sulfur (NaS) battery was selected due to its fast response time.It only takes 10 s to go from 0 to rated power [27], which means it could operate at rated power in 5 min.The technology for the BESS is max-charging power P max_ch and max-discharging power P max_f are both 3 MW; initial SOC is 0.5; the limit of SOC is 0.2-0.8.
Data for MPC: the predictive horizon and control horizon should satisfy the relation: m c ď m p .In this paper wind power is forecasted 30 min in the future according to the wind data during the last 12 h, as a consequence, in order to ensure sound control effects, the predictive horizon and control horizon are both set to be 30 min, i.e., m c = m p = 6, (30 min), α is 0.8.In [12], the authors proposed that when the SOC limit is between 0.3 and 0.7, the capacity of the BESS to dispatch the wind power would be 20%-30% of the rated capacity of the wind farm.As for this paper, the SOC would operate between 0.2-0.8,therefore, the capacity of the BESS would be lower than that in [12].In this paper, capacities of 10% and 20% would be selected to test the proposed strategy.Based on the predicted wind power using the ARIM model developed in [26] and the method in Section 2.1, the dispatch curve in a day was acquired and illustrated in Figure 4.

Capacity of BESS in 4.5 MWh
The dispatch strategy basing on MPC in Section 3 is applied to dispatch a wind farm whose rated capacity is 45 MW with a BESS of 45 MWh rate capacity, to maintain the total output from the wind/battery storage hybrid system to effectively track the dispatch order shown in Figure 4.Meanwhile, the effects of the proposed strategy are compared with those of a traditional solution for MPC, the QP method, which is illustrated in Figure 5.
In Figure 5 for both methods, the effects of the tracing are mostly nice.The IP method has a better performance on the whole even though there would be some worse tracing performance between 15 and 17 h.The statistical analysis on the absolute deviation between diapatched power and dispatch order in Figure 5 is illustrated in Table 1 and the statistical result of relative deviation is shown in Figure 6.

Capacity of BESS in 4.5 MWh
The dispatch strategy basing on MPC in Section 3 is applied to dispatch a wind farm whose rated capacity is 45 MW with a BESS of 45 MWh rate capacity, to maintain the total output from the wind/battery storage hybrid system to effectively track the dispatch order shown in Figure 4.Meanwhile, the effects of the proposed strategy are compared with those of a traditional solution for MPC, the QP method, which is illustrated in Figure 5.

Capacity of BESS in 4.5 MWh
The dispatch strategy basing on MPC in Section 3 is applied to dispatch a wind farm whose rated capacity is 45 MW with a BESS of 45 MWh rate capacity, to maintain the total output from the wind/battery storage hybrid system to effectively track the dispatch order shown in Figure 4.Meanwhile, the effects of the proposed strategy are compared with those of a traditional solution for MPC, the QP method, which is illustrated in Figure 5.
In Figure 5 for both methods, the effects of the tracing are mostly nice.The IP method has a better performance on the whole even though there would be some worse tracing performance between 15 and 17 h.The statistical analysis on the absolute deviation between diapatched power and dispatch order in Figure 5 is illustrated in Table 1 and the statistical result of relative deviation is shown in Figure 6.In Figure 5 for both methods, the effects of the tracing are mostly nice.The IP method has a better performance on the whole even though there would be some worse tracing performance between 15 and 17 h.
The statistical analysis on the absolute deviation between diapatched power and dispatch order in Figure 5 is illustrated in Table 1 and the statistical result of relative deviation is shown in Figure 6.As can be shown in Table 1, for the IN-MPC strategy proposed in this paper, the maximum deviation is 1.478 MW, however, for the QP-MPC, it is 1.537 MW, meanwhile, the mean value of the proposed one is 0.2047 MW which is smaller than the 0.4139 MW solved by the conventional one.Therefore, the dispatching results based on IP-MPC are better due to its algorithm being better than QP-MPC.The results can also be verified in Figure 6, where as can be seen, the frequency of smaller errors in IP-MPC would be less than that in QP-MPC.

Method
Figures 7 and 8 explain the reason for the "bad points".In Figure 6 both curves are strictly limited between 0.2 and 0.8, during 15-17 h, SOC for IP-MPC almost reaches the upper limits, therefore the BESS has no space for BESS to absorb energy.Furthermore, in Figure 8, the action of BESS is little and less than that in QP-MPC.The operating range of the BESS is limited within [−3, 3], most of which are in half the limited power region of [−1.5, 1.5].As a consequence, capacity is an important factor in the performance.In the next part, larger capacity would be selected.As can be shown in Table 1, for the IN-MPC strategy proposed in this paper, the maximum deviation is 1.478 MW, however, for the QP-MPC, it is 1.537 MW, meanwhile, the mean value of the proposed one is 0.2047 MW which is smaller than the 0.4139 MW solved by the conventional one.Therefore, the dispatching results based on IP-MPC are better due to its algorithm being better than QP-MPC.The results can also be verified in Figure 6, where as can be seen, the frequency of smaller errors in IP-MPC would be less than that in QP-MPC.
Figures 7 and 8 explain the reason for the "bad points".In Figure 6 both curves are strictly limited between 0.2 and 0.8, during 15-17 h, SOC for IP-MPC almost reaches the upper limits, therefore the BESS has no space for BESS to absorb energy.Furthermore, in Figure 8, the action of BESS is little and less than that in QP-MPC.The operating range of the BESS is limited within [´3, 3], most of which are in half the limited power region of [´1.5, 1.5].As a consequence, capacity is an important factor in the performance.In the next part, larger capacity would be selected.

Capacity of BESS in 9 MWh
In this section, the capacity of the BESS is doubled, which is 20% rated capacity of the wind farm, i.e., 9 MWh.The proposed strategy is applied to dispatch and control active power from the wind/battery storage hybrid system for one day (24 h), and the results are shown in Figure 9.In Figure 9, both have a better performance than that in Figure 5 because of the larger capacity of the BESS.The total output under the control of IP could trace the dispatch order perfectly, and the performance of IP-MPC is better than that of QP-MPC.

Capacity of BESS in 9 MWh
In this section, the capacity of the BESS is doubled, which is 20% rated capacity of the wind farm, i.e., 9 MWh.The proposed strategy is applied to dispatch and control active power from the wind/battery storage hybrid system for one day (24 h), and the results are shown in Figure 9.In Figure 9, both have a better performance than that in Figure 5 because of the larger capacity of the BESS.The total output under the control of IP could trace the dispatch order perfectly, and the performance of IP-MPC is better than that of QP-MPC.

Capacity of BESS in 9 MWh
In this section, the capacity of the BESS is doubled, which is 20% rated capacity of the wind farm, i.e., 9 MWh.The proposed strategy is applied to dispatch and control active power from the wind/battery storage hybrid system for one day (24 h), and the results are shown in Figure 9.In Figure 9, both have a better performance than that in Figure 5 because of the larger capacity of the BESS.The total output under the control of IP could trace the dispatch order perfectly, and the performance of IP-MPC is better than that of QP-MPC.
The statistical analysis on both absolute deviation and relative deviation from the dispatching order is shown in Table 2 and Figure 10.The comparison between Tables 1 and 2 shows that due to the larger capacity of the BESS, the maximum deviation is only MW with the proposed strategy, and the tracking performance is better than that with 4.5 MWh.The statistical analysis on both absolute deviation and relative deviation from the dispatching order is shown in Table 2 and Figure 10.The comparison between Tables 1 and 2 shows that due to the larger capacity of the BESS, the maximum deviation is only MW with the proposed strategy, and the tracking performance is better than that with 4.5 MWh.As for errors in Figure 12, IP-MPC is less error prone than QP-MPC.Figure 11 explains the effects in Figure 9, where neither has reached the limiting bounds of 0.2 and 0.8, which means the BESS has sufficient energy to charge and discharge power to maintain an excellent performance.In Figure 12, the action of the BESS operates in the limited range, and most time, both run between −1.5 and 1.5.The statistical analysis on both absolute deviation and relative deviation from the dispatching order is shown in Table 2 and Figure 10.The comparison between Tables 1 and 2 shows that due to the larger capacity of the BESS, the maximum deviation is only MW with the proposed strategy, and the tracking performance is better than that with 4.5 MWh.As for errors in Figure 12, IP-MPC is less error prone than QP-MPC.Figure 11 explains the effects in Figure 9, where neither has reached the limiting bounds of 0.2 and 0.8, which means the BESS has sufficient energy to charge and discharge power to maintain an excellent performance.In Figure 12, the action of the BESS operates in the limited range, and most time, both run between −1.5 and 1.5.As for errors in Figure 12, IP-MPC is less error prone than QP-MPC.Figure 11 explains the effects in Figure 9, where neither has reached the limiting bounds of 0.2 and 0.8, which means the BESS has sufficient energy to charge and discharge power to maintain an excellent performance.In Figure 12, the action of the BESS operates in the limited range, and most time, both run between ´1.5 and 1.5.From both simulations, when the capacity of the BESS is small, IP-MPC has a better performance in tracking dispatch orders8 than QP-MPC.When the capacity is bigger, the tracking results of both are almost the same.Due to the high cost of the BESS with big capacity, IP-MPC has an edge in the market.

Conclusions
This paper developed a MPC-based dispatching strategy for wind/battery energy storage hybrid systems.Firstly, the dispatch curve in the dispatching interval is obtained using the predicted wind power, and then the MPC controller calculates the control value for the BESS.The simulation results show that: (a) the proposed strategy based on MPC could enable the active output of a wind/battery energy storage hybrid system to trace the dispatch order, which improves the schedulablity of the wind farm; (b) the IP-MPC has a better performance in tracking dispatch orders than QP-MPC under the same BESS conditions, such as capacity, limited SOC and charging/discharging power limits.This paper focuses on the strategy for BESS based on MPC to  From both simulations, when the capacity of the BESS is small, IP-MPC has a better performance in tracking dispatch orders8 than QP-MPC.When the capacity is bigger, the tracking results of both are almost the same.Due to the high cost of the BESS with big capacity, IP-MPC has an edge in the market.

Conclusions
This paper developed a MPC-based dispatching strategy for wind/battery energy storage hybrid systems.Firstly, the dispatch curve in the dispatching interval is obtained using the predicted wind power, and then the MPC controller calculates the control value for the BESS.The simulation results show that: (a) the proposed strategy based on MPC could enable the active output of a wind/battery energy storage hybrid system to trace the dispatch order, which improves the schedulablity of the wind farm; (b) the IP-MPC has a better performance in tracking dispatch orders than QP-MPC under the same BESS conditions, such as capacity, limited SOC and charging/discharging power limits.This paper focuses on the strategy for BESS based on MPC to From both simulations, when the capacity of the BESS is small, IP-MPC has a better performance in tracking dispatch orders8 than QP-MPC.When the capacity is bigger, the tracking results of both are almost the same.Due to the high cost of the BESS with big capacity, IP-MPC has an edge in the market.

Conclusions
This paper developed a MPC-based dispatching strategy for wind/battery energy storage hybrid systems.Firstly, the dispatch curve in the dispatching interval is obtained using the predicted wind power, and then the MPC controller calculates the control value for the BESS.The simulation results show that: (a) the proposed strategy based on MPC could enable the active output of a wind/battery energy storage hybrid system to trace the dispatch order, which improves the schedulablity of the wind farm; (b) the IP-MPC has a better performance in tracking dispatch orders than QP-MPC under the same BESS conditions, such as capacity, limited SOC and charging/discharging power limits.This paper focuses on the strategy for BESS based on MPC to smooth and dispatch wind power.Our future work will investigate the errors in wind power prediction and the stochastic nature of wind power.

Figure 1 .
Figure 1.Structure chart of wind/battery energy storage hybrid system.

Figure 1 .
Figure 1.Structure chart of wind/battery energy storage hybrid system.

Figure 3 .
Figure 3. Diagram of smoothing and dispatching wind power via MPC.Figure 3. Diagram of smoothing and dispatching wind power via MPC.

Figure 3 .
Figure 3. Diagram of smoothing and dispatching wind power via MPC.Figure 3. Diagram of smoothing and dispatching wind power via MPC.