MPC Based Coordinated Active and Reactive Power Control Strategy of DFIG Wind Farm with Distributed ESSs

: The ESS is considered as an effective tool for enhancing the ﬂexibility and controllability of a wind farm, and the optimal control scheme of a wind farm with distributed ESSs is vital to the stable operation of wind power generation. In this paper, a coordinated active and reactive power control strategy based on model predictive control (MPC) is proposed for doubly fed induction generator (DFIG)-based wind farm (WF) with distributed energy storage systems (ESSs). The proposed control scheme coordinates the active and reactive power output among DFIG wind turbines (WTs), grid-side converters (GSCs), and distributed ESSs inside the WF, and the aim is to decrease fatigue loads of WTs, make the WT terminal voltage inside the extent practicable, and take the WF economic operation into consideration. Moreover, the best reactive power references of DFIG stator and GSC are produced independently based on their dynamics. At last, the control scheme generates optimal power references for all ESS to make the SOC of each ESS converge to their average state. With the distributed ESSs, the WF controller regulates the WTs inside WF more ﬂexibly. A WF composed of 10 DFIG WTs was utilized to verify the control performance of the proposed coordinated active and reactive power control strategy.


Introduction
In order to deal with the environmental issue and depletion of fossil fuels, wind power has been developing rapidly over the world [1][2][3]. The increasing installed capacity of wind farms (WFs) brings challenges to power system operation and control [4]. WFs are required to operate as conventional power plants to preserve safety of the power system with high wind power penetration rate. To solve these challenges, the WFs have to use reasonable control strategy to meet grid operating demands from the transmission system operators (TSOs) [5]. The DFIG-based wind turbine (WT) has been broadly applied for in advanced WFs benefit from excellent controllability and little converter capacity. By means of power electronic converters, the reactive power of DFIG-based WF could be separated adjusted and the power system reactive power support could be offered [6].
In the conventional active power dispatch strategy of WFs, the power reference is proportionally distributed to individual WTs on the basis of either WTs available power or actual output power. Without considering the operating condition of each WT, the fatigue load increases which results in reduction of WT lifetime [7,8]. For minimizing fatigue loads, in [9], a MPC-based active power control strategy for WFs was proposed to decrease the WTs fatigue loads.
Due to the energy storage system (ESS) has good performance of adaptable charging/discharging, it can be used as an useful tool to improve the adaptability and controllability of WF. For a WF with a centralized ESS [10,11], the requirements for electronic devices of centralized ESS will be rapidly increasing with WF size increasing. If the centralized ESS fails, it will result in security operation problems for the wind farm. For a WF with distributed ESSs, each WT inside the WF is equipped with a small capacity ESS on the DC side of the WT converter. One ESS failure will not have a big impact on the WF operation. WT with distributed ESS has motivated a great number of studies [12][13][14]. In [13], a constant power control strategy was presented for a DFIG WT WF, at the same time, ESS is installed on every DFIG. The scheme enabled the WF to effectively regulate active power of the grid. Reference [14] has proposed a model and control framework to attain coordinated run for grid connected DFIG-based WT and ESS.
For the purpose of the DFIG-based WF safe operation, voltage performance also should be addressed. Every DFIG-based WT inside DFIG-based WFs is prepared with power electronic converter. The DFIG-based WF can utilize the capabilities of the electronic converters to regulate the bus voltages inside the WF. A few control modes involve voltage, reactive power at the point of connection (POC) and power factor have been designated in a lot of grid codes [15]. Voltage control mode usually exhibit excellent execution of transmission systems [16,17]. In [18], A optimized voltage control strategy of offshore WFs cluster was proposed, the pilot bus voltage is controlled under steady-state operation by assigning reactive power references to every WT according to their accessible reactive power margin. An ideal reactive power dispatch strategy was presented in [19] to keep WF electrical losses at a minimum, containing transmission cables losses, WT transformers, and the wind energy generation systems. In [20], a hierarchical automatic voltage control was presented to keep every WT terminal voltage inside the attainable extend.
As mentioned above, there are some currently studies researches have made incredible contributions to the WFs voltage/reactive power control. But most of references described the dynamic model of WT as a first-order lag function, and did not consider the real WT generator and converter dynamic model [21][22][23]. Therefore, an MPC-based active and reactive power integrated control strategy for DFIG WF with distributed ESSs is proposed in this paper, which has two stages. The control objectives of the first stage are controlling the DFIGs inside the WF to follow the active power reference issued by TSO, decrease fatigue loads through minimizing the change of WT thrust force and shaft torque, meanwhile generate a ESSs active power reference. Under the second stage, the controller offer optimized reactive power references of DFIG stator and grid-side converter to adjust the voltage of each bus in the reasonable range, and reduce the WF network losses. Moreover, controller generates the active power references for ESSs to keep its SOCs close to a common state and consider the WF economic operation.
The most important contribution of this paper is the MPC based coordinated active and reactive power control strategy for distributed large-scale WFs with ESSs. The ideal reactive power references of DFIG stator and GSC are generated separately to keep bus voltage deviation and network power losses at a minimum. With the distributed ESSs, the WF controller regulates WTs more flexibly. The WTs and ESSs active and reactive power output could be best distributed to promote the optimal operation of WF. Compared with conventional reactive power control for DFIGs WF, reactive power references are generated separately on the basis of stator and GSC dynamics to achieve a better performance.
The paper includes the following parts. Section 2 introduces a brief overview of the coordinated active and reactive power control structure. Section 3 introduces the prediction model for RSC, GSC, ESS and WT in detail. Section 4 presents the MPC based mathematical formulation of coordinated power control scheme is described. Section V introduces and discusses the simulation results, and gives conclusions. This article takes DFIG-based WF with distributed ESS as the research object, different from the traditional control method that control WT active power independently to reduce fatigue load or control WT reactive power independently to stabilize WT terminal voltage, the proposed control method, the proposed control method can coordinate active and reactive power control at the same time, reduce WT fatigue load and stabilize WT voltage at the same time, and reduce the WF network losses. meanwhile, tracking the dispatch command and maintaining the state-of-charge (SOC) of ESSs within a specified range.

Configuration of DFIG Wind Farm
A representative topology of a WF is shown in Figure 1. It is linked to a outside AC grid by transmission cable. The HV/MV on-load tap changer (OLTC) is located between point of connection (POC) in HV side and collection point in MV side. Several feeders are connected to the MV collection point, and each feeder is connected to series WTs. reduce the WF network losses. meanwhile, tracking the dispatch command and maintaining the state-of-charge (SOC) of ESSs within a specified range.

Configuration of DFIG Wind Farm
A representative topology of a WF is shown in Figure 1. It is linked to a outside AC grid by transmission cable. The HV/MV on-load tap changer (OLTC) is located between point of connection (POC) in HV side and collection point in MV side. Several feeders are connected to the MV collection point, and each feeder is connected to series WTs.

Concept of Coordinated Power Control Scheme
The structure of the presented active and reactive power coordinated control strategy is shown in Figure 2. The WF total active power reference is set by TSO and send to the WF controller. The MPC-based WF controller implements a two-stage optimization scheme.
Configuration of a wind farm.

Concept of Coordinated Power Control Scheme
The structure of the presented active and reactive power coordinated control strategy is shown in Figure 2. The WF total active power reference is set by TSO and send to the WF controller. The MPC-based WF controller implements a two-stage optimization scheme.

Configuration of DFIG Wind Farm
A representative topology of a WF is shown in Figure 1. It is linked to a outside AC grid by transmission cable. The HV/MV on-load tap changer (OLTC) is located between point of connection (POC) in HV side and collection point in MV side. Several feeders are connected to the MV collection point, and each feeder is connected to series WTs.

Concept of Coordinated Power Control Scheme
The structure of the presented active and reactive power coordinated control strategy is shown in Figure 2. The WF total active power reference is set by TSO and send to the WF controller. The MPC-based WF controller implements a two-stage optimization scheme. , , Figure 2. Structure of control scheme. The first stage controls all DFIGs inside the WF to follow the active power reference issued by TSO, and decreasing WTs fatigue loads, brought by the changes of WTs shaft torque and thrust force. And the total active power reference value of ESSs and the best active power reference of every DFIG WT are obtained by the controller.
In the second stage control, controller make the ESSs follow its power reference from the first stage control, while keep its SOCs to an average state. Meanwhile controller regulates the reactive power output of DFIG stator and GSC to minimize bus voltage deviation and decrease network power losses of WF, and consider the influence of active power outputs of WTs and ESSs. Figure 3 shows the structure of the DFIG WT equipped with an ESS. DFIG stator joins to the AC grid directly. DFIG rotor is linked to the AC network by a back-to-back PWM converter, including a rotor side converter (RSC), a GSC and an ESS.ESS contains an energy storage unit (ESU) and a DC/DC converter. The ESU is linked to a DC bus by a DC/DC converter to adjust charge/discharge power of ESU to track the active power reference from the WF controller. The first stage controls all DFIGs inside the WF to follow the active power referen issued by TSO, and decreasing WTs fatigue loads, brought by the changes of WTs sh torque and thrust force. And the total active power reference value of ESSs and the b active power reference of every DFIG WT are obtained by the controller.

WF Model
In the second stage control, controller make the ESSs follow its power reference fro the first stage control, while keep its SOCs to an average state. Meanwhile controller re ulates the reactive power output of DFIG stator and GSC to minimize bus voltage dev tion and decrease network power losses of WF, and consider the influence of active pow outputs of WTs and ESSs. Figure 3 shows the structure of the DFIG WT equipped with an ESS. DFIG stator joi to the AC grid directly. DFIG rotor is linked to the AC network by a back-to-back PW converter, including a rotor side converter (RSC), a GSC and an ESS.ESS contains an e ergy storage unit (ESU) and a DC/DC converter. The ESU is linked to a DC bus by DC/DC converter to adjust charge/discharge power of ESU to track the active power r erence from the WF controller.

RSC Model
The RSC is mainly used to control active and reactive power output of DFIG. R could decouple active and reactive power by adjusting rotor current in a synchronou rotating reference frame oriented by the stator flux. Figure 4 shows the RSC control stru ture.
DFIG active power output is: stator active power output is the first item in function (1), and GSC active power outp is the second item. DFIG reactive power is DFIG stator reactive power QS. QS is controlled through a justing DFIG rotor d-axis current idr, its calculation formula can be derived and express as:

RSC Model
The RSC is mainly used to control active and reactive power output of DFIG. RSC could decouple active and reactive power by adjusting rotor current in a synchronously rotating reference frame oriented by the stator flux. Figure 4 shows the RSC control structure. where ψS is the stator flux, Lm is the mutual inductance, ωS is the supply angular speed, LS is DFIG stator inductance. Uqr is q-axis rotor voltage, iqr is rotor q-axis current.    DFIG active power output is: stator active power output is the first item in function (1), and GSC active power output is the second item. DFIG reactive power is DFIG stator reactive power Q S . Q S is controlled through adjusting DFIG rotor d-axis current i dr , its calculation formula can be derived and expressed as: where ψ S is the stator flux, L m is the mutual inductance, ω S is the supply angular speed, L S is DFIG stator inductance. U qr is q-axis rotor voltage, i qr is rotor q-axis current. Normally, The time constant of power outer loop is much smaller than that of the current inner loop The first-order lag function could describe the dynamic behavior of the current loop. Hence, the model of RSC active power and reactive power could be simplified re which is shown in Figures 5 and 6 spectively. Where k P r_P , k Q r_P , k P r_i , and k Q r_i are the proportional gain and integral gain of PI controllers, respectively. T ir is the time constant of current loop. T fr denote filter time constant. Superscript '*'indicates reference value. where ψS is the stator flux, Lm is the mutual inductance, ωS is the supply angular speed, LS is DFIG stator inductance. Uqr is q-axis rotor voltage, iqr is rotor q-axis current.       Define ∆i qr , ∆i dr , ∆P g , ∆Q S , ∆P ref g , ∆P int , ∆Q ref S , and ∆Q S int are the incremental values between its current value and the initial value at operation point, respectively. Where P int is the error integral of P ref g and P g , Q S int is the integral of error between Q ref S and Q S . Then, RSC active power incremental state space model is as follows, The state space matrix form is, where The RSC reactive power incremental state space model is derived as follows, The state space matrix form is, where

GSC Model
The DC bus voltage is adjusted by GSC, which offering reactive power support for WF AC grid at the same time. GSC can decouple DC voltage and reactive power control by adjusting current in a synchronously rotating reference frame oriented by the grid voltage [24]. Figure 7 shows the GSC control structure. GSC reactive power in the grid-voltage-oriented synchronously rotating reference frame can be modeled as, where Um is the grid phase voltage amplitude. The same as the RSC, the GSC model is shown in Figure 8, k  GSC reactive power in the grid-voltage-oriented synchronously rotating reference frame can be modeled as, where U m is the grid phase voltage amplitude. The same as the RSC, the GSC model is shown in Figure 8, k Q g_P , k Q g_i is the proportional gain and integral of PI controller, T ig denote current loop time constant, T fg denote GSC filter time constant. Superscript '*'indicates reference value.
where Um is the grid phase voltage amplitude. The same as the RSC, the GSC model is shown in Figure 8, k Q g_P , k Q g_i is the proportional gain and integral of PI controller, Tig denote current loop time constant, Tfg denote GSC filter time constant. Superscript '*'indicates reference value.  Define the ∆i q , ∆Q C , ∆Q ref C , and ∆Q C int are the incremental values between its current value and the initial value at operation point, respectively, and Q C int is the integral of error between Q ref C and Q C . Superscript '*'indicates reference value. GSC incremental state space model is derived as follows, The state space matrix form can be get, where

ESS Model
The calculation formula of ESS stored energy C ESS is, where C ESS,0 denotes the initial energy, and P ESS denotes the charge/discharge power of ESS.
To enhance the DC side current response speed, a double closed loop control structure diagram of DC/DC converter, as shown in Figure 9. Figure 10 shows the ESS simplified block diagram. Where k P d_P and k P d_i denotes the proportional gain and integral of the PI controller. U ESS denotes ESU voltage. T id denotes the current loop time constant, T fd denotes ESS filter time constant. Superscript '*' indicates reference value.
where CESS,0 denotes the initial energy, and PESS denotes the charge/discharge power of ESS.
To enhance the DC side current response speed, a double closed loop control structure diagram of DC/DC converter, as shown in Figure 9. Figure 10 shows the ESS simplified block diagram. Where k P d_P and k P d_i denotes the proportional gain and integral of the PI controller. UESS denotes ESU voltage. Tid denotes the current loop time constant, Tfd denotes ESS filter time constant. Superscript '*'indicates reference value.   where CESS,0 denotes the initial energy, and PESS denotes the charge/discharge power of ESS.
To enhance the DC side current response speed, a double closed loop control structure diagram of DC/DC converter, as shown in Figure 9. Figure 10 shows the ESS simplified block diagram. Where k P d_P and k P d_i denotes the proportional gain and integral of the PI controller. UESS denotes ESU voltage. Tid denotes the current loop time constant, Tfd denotes ESS filter time constant. Superscript '*'indicates reference value.
The state space expression of the ESS is, where

WT Model
The mechanical system model of the WT can be obtained according to the analysis in [25].
Let K RSC = L m ψ S ω S T fr L S + 3U qr 2T fr L S , and combine the RSC active power model and the WT model, the incremental state-space matrix of DFIG active power dynamic model is, where ∆x W = ∆P g , ∆P int , ∆i qr , ∆β ref , ∆β, ∆ω g , ∆ω f T , where θ denotes pitch angle, β = K C θ, K C = K 0 + 2K 1 θ, where K 0 and K 1 are constants, T θ denotes time constant of pitch servo function, T a denotes aerodynamic torque, and T g denotes generator torque, K p and K i denote the proportional gain and integral gain of pitch controller, ω f is filtered speed at which the generator speed ω g is filtered by the low-pass filter, T f is filter time constant, η g denotes gear box ratio, J r denotes rotor mass, J g denotes generator mass, J t = J r + η 2 g J g is the equivalent mass, µ g is the generator efficiency; T S denotes shaft torque, F t denotes thrust force, K βT , K ωT , K βF and K ωF are the coefficients obtained the Taylor expansion of T S and F t at the operating point, and subscript 0 is used to denote the initial value in control period.

Whole System Model
The continuous state-space matrix of a WF consist of N T WTs can be get as follows, where The discrete state space matrix is,

MPC Based Coordinated Active and Reactive Power Control for Wind Farm with ESSs
The DFIG-based WF requires the WF controller consider the coordinated control of the all the energy sources inside the WF. All WTs, RSCs, GSCs, and VSCs within the WF are optimally coordinated to reduce fatigue loads, keep the each bus voltage in a reasonable range and decrease network power losses. In order to decouple active and reactive power optimal control and keep fatigue loads at a minimum, the control strategy is divided into two steps. every DFIG active power reference and a total ESSs active power reference are produced by the controller in the first stage, which aims to follow the dispatch command and reduce the fatigue loads by minimizing variations in TS and Ft of WTs. The controller coordinates control the reactive power of DFIG stator and GSC inside this WF in the second stage, aims to keep the all buses voltage of WF in a feasible range, and reduces WF network power losses. Moreover, the WF controller regulates ESSs charge/discharge power to follow the ESSs active power reference obtained in the first step and keeps every ESS SOC value converging to average.

Objective Function
The WF controller reduces fatigue loads of WTs. In addition, in order to enable the ESS to operate stably for a long time, the CESS of ESS is controlled to keep at the intermediate level. Then the objective function can beget as follows, where N P denotes MPC predictive steps, and C total ESS,mid denotes intermediate level of total ESS capacity. λ C , λ T , λ F are the weighting coefficients for the total ESSs energy management, variations of T S and F t , separately.

Constraints
The active power constraint of the WF is, The constraint of a WT is as follows, The constraint of total ESSs is as follows,

Objective Function
The first objective is voltage regulation, which aims to keep the each MV bus voltage and WT terminal bus voltage in a reasonable range [26]. Therefore, the objective function about voltage regulation is, Keep the total network power losses at a minimum is the second objective, which is, where P pre LOSS is the predictive value of network power losses, where sensitivity coefficients regarding power injections in (34), (35), In addition, (37) can be obtained by using the analytical method in [26].
The third objective is to ensure the SOC of each ESS converging to an average value, which is, According to (33), (36), (38), the objective function of the second step is, where λ V , λ L , λ S are the weighting coefficients for Obj V , Obj L , Obj S , respectively.

Constraints
For a distributed ESS, its active power reference constraint is as follow, where Q min S,i /Q min C,i and Q max S,i /Q max C,i are the minimum and maximum Var capacity of the DFIG stator and GSC, separately.
The entire charge/discharge active power output should be set equal to the total ESS active power reference obtained by the first stage as follow,

Test System
A WF composed of 10 × 5 MW DFIG WTs with distributed ESSs is used to verify presented control strategy in the section. The WF modeling that considers the effects of turbulence and wakes is created by SimWindFarm, a toolbox for dynamic WF modeling, simulation and control. The control period T C , prediction horizon T P and the prediction sampling period ∆T P are set as 1 s, 5 s, and 0.1 s, separately. For the purpose of verify the effectiveness of the presented control scheme, the simulation results are contrast with the one without ESSs and reactive power control.

Control Performance
In the MATLAB/Simulink, the simulation time is set as 550 s. The total WF available active power and TSO dispatch demand are shown in Figure 11. During t = 0-280 s, the WF available power fluctuates between 38-50 MW. At t = 280 s, the total available active power begin increasing and reaches the maximum available power 50 MW at 300 s. The available active power gradually reduces after t = 480 s.  Figure  12, thanks to the control effect of ESSs on WF active power adjustment, Ft variation with ESSs is smaller than that without ESSs. In Figure 13, during t = 0-340 s, the fluctuations of TS with ESSs is much less than that without ESSs. After t = 340 s, wind farm is operating in the maximum power point tracking (MPPT) mode, the changes in TS with ESSs is also less than that without ESSs. The fluctuations of thrust force are also less than that without ESSs. The control strategy proposed in this paper can effectively decrease the fatigue loads of DFIG WT by minimizing the changes of Ft and TS, then improve the security of the WF. Figure 11. WF total available and reference active power. Figures 12 and 13 show the thrust force and shaft torque changes of WT4. In Figure 12, thanks to the control effect of ESSs on WF active power adjustment, F t variation with ESSs is smaller than that without ESSs. In Figure 13, during t = 0-340 s, the fluctuations of T S with ESSs is much less than that without ESSs. After t = 340 s, wind farm is operating in the maximum power point tracking (MPPT) mode, the changes in T S with ESSs is also less than that without ESSs. The fluctuations of thrust force are also less than that without ESSs. The control strategy proposed in this paper can effectively decrease the fatigue loads of DFIG WT by minimizing the changes of F t and T S , then improve the security of the WF.

The First Stage Control Performance
12, thanks to the control effect of ESSs on WF active power adjustment, Ft variation with ESSs is smaller than that without ESSs. In Figure 13, during t = 0-340 s, the fluctuations of TS with ESSs is much less than that without ESSs. After t = 340 s, wind farm is operating in the maximum power point tracking (MPPT) mode, the changes in TS with ESSs is also less than that without ESSs. The fluctuations of thrust force are also less than that without ESSs. The control strategy proposed in this paper can effectively decrease the fatigue loads of DFIG WT by minimizing the changes of Ft and TS, then improve the security of the WF.    Figures 12 and 13 show the thrust force and shaft torque changes of WT4. In Figure  12, thanks to the control effect of ESSs on WF active power adjustment, Ft variation with ESSs is smaller than that without ESSs. In Figure 13, during t = 0-340 s, the fluctuations of TS with ESSs is much less than that without ESSs. After t = 340 s, wind farm is operating in the maximum power point tracking (MPPT) mode, the changes in TS with ESSs is also less than that without ESSs. The fluctuations of thrust force are also less than that without ESSs. The control strategy proposed in this paper can effectively decrease the fatigue loads of DFIG WT by minimizing the changes of Ft and TS, then improve the security of the WF.

Simulation Analysis of the Second Stage
The WT 4 terminal voltage is shown in Figure 14. With the reactive power control, the voltage could be kept within the reasonable range effectively. The voltage deviation from nominal value with the reactive power control is significantly smaller than that without this control, which guarantee the WF secure operation.

Simulation Analysis of the Second Stage
The WT 4 terminal voltage is shown in Figure 14. With the reactive power control, the voltage could be kept within the reasonable range effectively. The voltage deviation from nominal value with the reactive power control is significantly smaller than that without this control, which guarantee the WF secure operation.   Figures 15 and 16 show the active power and reactive power output of WT 4, respectively. During t = 0-300 s, the WT 4 stator and GSC absorb reactive power separately to regulate the WF voltages. After t = 300 s, due to the increasing active power output, the reactive power absorbed by stator is decreasing. At this time, the GSC supply reactive power. Due to GSC reactive power output limits, the WT 4 terminal voltage is increased, but it is still within the feasible operating range. During t = 0-300 s, the WT 4 stator and GSC absorb reactive power separately to regulate the WF voltages. After t = 300 s, due to the increasing active power output, the reactive power absorbed by stator is decreasing. At this time, the GSC supply reactive power. Due to GSC reactive power output limits, the WT 4 terminal voltage is increased, but it is still within the feasible operating range.  The WF total network power losses is shown in Figure 17. The network power losses with reactive power control is smaller than that without reactive power control.   The WF total network power losses is shown in Figure 17. The network power losses with reactive power control is smaller than that without reactive power control. The WF total network power losses is shown in Figure 17. The network power losses with reactive power control is smaller than that without reactive power control. Figure 17. Network power losses of wind farm. Figure 18 shows the SOC of each ESS. All the ESSs SOC are kept to the average value. From t = 0-350 s, the SOC fluctuates around 0.5, which means the proposed scheme manage the ESSs to keep the SOCs at the medium level. From t = 350-500 s, the WF output rated active power. The ESSs also needs output a spot of active power to decrease the WTs   Figure 18 shows the SOC of each ESS. All the ESSs SOC are kept to the average value. From t = 0-350 s, the SOC fluctuates around 0.5, which means the proposed scheme manage the ESSs to keep the SOCs at the medium level. From t = 350-500 s, the WF output rated active power. The ESSs also needs output a spot of active power to decrease the WTs fatigue loads. The SOCs are decreased slightly, but they are kept closing to the medium level. After 500 s, the available power gradually decreases, and the dispatch command remain unchanged. The available power is less than the dispatch command. Therefore, all ESSs discharge power to track the dispatch command, causing SOCs decrease gradually.  Figure 19 shows the WF total active power output that equipped with distributed ESSs. With the ESSs, the WF total active power output could accurately follow the active power dispatch command from the TSO. Figure 19. Active power output of wind farm.

Conclusions
In this paper, an MPC based coordinated active and reactive power control strategy is proposed for DFIG wind farm with distributed ESSs, which coordinates the active and reactive power output among DFIG WTs, GSC and distributed ESSs inside wind farm to track the dispatch command the TSO set, reduce the fatigue loads of the WTs, and manage the distributed ESSs. Moreover, the control scheme can provide the best reactive power references for DFIG stator and GSC separately to regulate all the bus voltages and network power losses while track the dispatch command from TSO. Case studies show the control strategy could effectively decrease the fatigue loads, regulate the voltage in the feasible range, decrease the WF network power losses and ensure the SOC of each ESS converging to an average value. The presented coordinated active and reactive power control strategy is suitable for real-time control of large-scale wind farm.   Figure 19 shows the WF total active power output that equipped with distributed ESSs. With the ESSs, the WF total active power output could accurately follow the active power dispatch command from the TSO.  Figure 19 shows the WF total active power output that equipped with distributed ESSs. With the ESSs, the WF total active power output could accurately follow the active power dispatch command from the TSO. Figure 19. Active power output of wind farm.

Conclusions
In this paper, an MPC based coordinated active and reactive power control strategy is proposed for DFIG wind farm with distributed ESSs, which coordinates the active and reactive power output among DFIG WTs, GSC and distributed ESSs inside wind farm to track the dispatch command the TSO set, reduce the fatigue loads of the WTs, and manage the distributed ESSs. Moreover, the control scheme can provide the best reactive power references for DFIG stator and GSC separately to regulate all the bus voltages and network power losses while track the dispatch command from TSO. Case studies show the control strategy could effectively decrease the fatigue loads, regulate the voltage in the feasible range, decrease the WF network power losses and ensure the SOC of each ESS converging to an average value. The presented coordinated active and reactive power control strategy is suitable for real-time control of large-scale wind farm.

Conclusions
In this paper, an MPC based coordinated active and reactive power control strategy is proposed for DFIG wind farm with distributed ESSs, which coordinates the active and reactive power output among DFIG WTs, GSC and distributed ESSs inside wind farm to track the dispatch command the TSO set, reduce the fatigue loads of the WTs, and manage the distributed ESSs. Moreover, the control scheme can provide the best reactive power references for DFIG stator and GSC separately to regulate all the bus voltages and network power losses while track the dispatch command from TSO. Case studies show the control strategy could effectively decrease the fatigue loads, regulate the voltage in the feasible range, decrease the WF network power losses and ensure the SOC of each ESS converging to an average value. The presented coordinated active and reactive power control strategy is suitable for real-time control of large-scale wind farm. Data Availability Statement: The study did not report any data.

Conflicts of Interest:
The authors declare no conflict of interest.