Coordinated Secondary Frequency Regulation Strategy of Doubly-Fed Induction Generator and Electric Vehicle

: Wind turbines can participate in frequency regulation by controlling active power output, but the indeterminacy and volatility of wind power result in low reliability of frequency support. Therefore, as a kind of energy storage system, an electric vehicle is adopted to coordinate with wind turbines to regulate system frequency considering its large-scale development. First, based on the reasonable division of wind speed regions and operation point selection of pitch angle, the de-loading strategy of doubly-fed induction generator for reserve capacity under continuously varying wind speed is proposed. Then, through the combination of rotor speed and pitch angle control, frequency regulation model of a doubly-fed induction generator in whole wind speed range is established. Finally, taking into account the driving demand of electric vehicle owners, through the real-time allocation of system frequency regulation task based on frequency regulation capacity, the coordinated control strategy of doubly-fed induction generator and electric vehicle cluster for secondary frequency regulation is put forward. The simulation results show that the coordinated frequency regulation strategy based on real-time allocation can suppress frequency deviation e ﬀ ectively, and the regulation e ﬀ ect is better than the situations of wind turbine coordinating with the conventional unit or coordinating with electric vehicle cluster based on ﬁxed allocation ratio. With the increase in wind power penetration in the power system, WTs are required to have the ability to control active power output. In order to compensate for DFIG’s weakness in FR caused by the indeterminacy and volatility of wind power, an electric vehicle is adopted as an auxiliary storage device to assist WTs to participate in system FR. This study establishes the FR model of DFIG under continuously varying WS in whole WS range, then proposes a coordinated control strategy of DFIG and EV cluster for system secondary FR with the real-time allocation of system FR task. The e ﬀ ectiveness of the proposed FR strategy is veriﬁed through various simulation cases based on a two-area interconnected system. the analysis of simulation following


Introduction
With the penetration rate of wind power in micro-grid increasing, the indeterminacy and volatility of wind power generation will directly affect the balance of supply and demand of active power in the micro-grid system, causing frequency fluctuations. Countries around the world have introduced grid codes, which stipulate that wind farms must be able to regulate active power and provide frequency support [1,2].
At present, the most widely used kind of wind turbines (WTs) is a double-fed induction generator (DFIG), which is connected to the grid through power electronic devices, thus is unable to respond to system frequency deviations without supplementary active power control methods, which includes rotor speed control, pitch angle control, and their combination. In Ref. [3], only pitch angle control is applied in the frequency regulation (FR) by DFIG, while the potential of rotor speed control is not exploited. However, the power regulation through rotor speed control is faster than pitch angle control, so its FR performance is better. The combination of two control methods is adopted in the FR strategy presented by Ref. [4], but the de-loading of DFIG is not studied; thus, DFIG can only participate in system to demand power that cannot be effectively fulfilled or the waste of power output. Moreover, in the most relevant researches, the influence of FR reserve capacity of DFIG or EV on the power output limit is not considered. Therefore, it is necessary to study the allocation method of system FR task considering the FR capacity so as to achieve the mutually coupled coordination of DFIG and EV in FR.
For the sake of solving the problems mentioned above, the main contributions of the paper are listed as follows: • Through reasonable dividing of WS regions and operation point selection of pitch angle, a modified de-loading strategy of DFIG for FR reserve capacity is proposed, which realizes the continuous change of rotor speed or pitch angle and guarantees expected FR capacity in every WS zone.

•
Combining rotor speed and pitch angle control, the FR model of DFIG for both up-regulation and down-regulation under continuously varying WS in whole WS range is established.

•
In the premise of meeting the driving demand of EV owners, based on the real-time allocation strategy of system FR task considering the real-time predicted FR capacity of DFIG and EV cluster, the coordinated control strategy of DFIG and EV in system secondary FR is presented.
The remainder of the paper is organized as follows. Section 2 presents the de-loading strategy and FR model of DFIG based on the combination of rotor speed and pitch angle control. Then, Section 3 elaborates the FR control strategy of EV cluster considering the driving demand of vehicle owners. The coordinated FR strategy of EV and DFIG based on the real-time allocation of system FR task is described in Section 4. The analysis of simulation results is shown in Section 5. Finally, conclusions are drawn in Section 6.

Power Output Control Principle of DFIG
By controlling the pitch angle and the converters on both sides of rotor and grid, the mechanical power output of wind turbine can be regulated through various controllers [15,16], as shown in (1): where P m is the active power output of DFIG, which equals to the mechanical power captured by DFIG while the loss of gearbox and transmission system is ignored; ρ air is the air density, whose value is 1.225 kg/m 3 ; R is the blade radius of DFIG; v w is the wind speed; C p is the blade power coefficient of DFIG, whose value is related to the pitch angle β and the tip-speed ratio λ; c 1 − c 7 are relevant parameters of DFIG characteristics [4], whose value are 0.22, 116, 0.4, 5, 12.5, 0.08, and 0.035, respectively; ω M is the angular velocity of blade, whose value is proportional to the rotor speed ω r , as shown in (2): where n is the gearbox transmission ratio of DFIG, whose value is constant for a fixed WT model. From Equations (1) and (2), it can be derived that P m changes with wind speed v w , rotor speed ω r , and pitch angle β, which can be expressed as P m (v w , ω r , β). Therefore, in the situation of fixed WS, the active power output of DFIG can be regulated through controlling the rotor speed and pitch angle. In this way, DFIG will be able to participate in system FR. As can be derived from (1), DFIG blade power coefficient C p changes with β and λ, and Figure 1 shows the change rule of C p . It can be seen that when pitch angle β = 0 • , there is an optimal tip-speed ratio λ opt corresponding to the maximum blade power coefficient C p,max , which can be reached by regulating the rotor speed to be the optimal rotor speed ω opt . Moreover, the DFIG power output can be reduced by increasing the pitch angle to a certain extent.
angle. In this way, DFIG will be able to participate in system FR. As can be derived from (1), DFIG blade power coefficient p C changes with β and λ , and Figure 1 shows the change rule of p C . It can be seen that when pitch angle β =°0 , there is an optimal tip-speed ratio λ opt corresponding to the maximum blade power coefficient ,max p C , which can be reached by regulating the rotor speed to be the optimal rotor speed ω opt . Moreover, the DFIG power output can be reduced by increasing the pitch angle to a certain extent.

Normal Operation of DFIG
For the normal operation of DFIG, in order to improve the power output efficiency of WT, the maximum power point tracking (MPPT) control strategy is adopted to achieve maximum capture of wind energy [17]. The rotor speed of the DFIG normally operates in the range of 0.7-1.2 p.u., and the safety range can reach 0.6-1.3 p.u. under transient dynamic disturbances [5]. In the MPPT state, the rotor speed of DFIG is adjusted dynamically to be the optimal rotor speed ω opt , so as to achieve the maximum power output power reaches the rated power e P , the corresponding optimal speed ω opt will go beyond the normal operation range; thus, the pitch angle will be increased to constrain rotor speed within its limit, keeping DFIG's power output at the rated power, as shown in Figure 2b, where e v is the rated WS corresponding to e P .  For the normal operation of DFIG, in order to improve the power output efficiency of WT, the maximum power point tracking (MPPT) control strategy is adopted to achieve maximum capture of wind energy [17]. The rotor speed of the DFIG normally operates in the range of 0.7-1.2 p.u., and the safety range can reach 0.6-1.3 p.u. under transient dynamic disturbances [5]. In the MPPT state, the rotor speed of DFIG is adjusted dynamically to be the optimal rotor speed ω opt , so as to achieve the maximum power output P MPPT (v w ) under different WSs, as shown in Figure 2a. If the maximum power reaches the rated power P e , the corresponding optimal speed ω opt will go beyond the normal operation range; thus, the pitch angle will be increased to constrain rotor speed within its limit, keeping DFIG's power output at the rated power, as shown in Figure 2b, where v e is the rated WS corresponding to P e .

Division of Wind Speed Region
For participating in down-regulation, the DFIG can reduce the power output directly through controlling rotor speed or pitch angle. However, DFIG cannot participate in up-regulation when operating in MPPT state since the power output has already reached the maximum value. In order to obtain up-regulation reserve capacity, the DFIG needs to operate in a de-loading state. In order to realize the de-loading of DFIG, the rotor speed needs to deviate from the optimal speed. The suboptimal operation state of rotor corresponds to two rotor speeds under the same WS, i.e., the lowspeed and over-speed suboptimal rotor speed. In the actual situation, the over-speed suboptimal trajectory is usually adopted in order to maintain DFIG's stable operation.
However, in consideration of the operation limit of rotor speed, the de-loading level may not be

Division of Wind Speed Region
For participating in down-regulation, the DFIG can reduce the power output directly through controlling rotor speed or pitch angle. However, DFIG cannot participate in up-regulation when operating in MPPT state since the power output has already reached the maximum value. In order to obtain up-regulation reserve capacity, the DFIG needs to operate in a de-loading state. In order to Appl. Sci. 2019, 9, 2815 5 of 24 realize the de-loading of DFIG, the rotor speed needs to deviate from the optimal speed. The suboptimal operation state of rotor corresponds to two rotor speeds under the same WS, i.e., the low-speed and over-speed suboptimal rotor speed. In the actual situation, the over-speed suboptimal trajectory is usually adopted in order to maintain DFIG's stable operation.
However, in consideration of the operation limit of rotor speed, the de-loading level may not be satisfied if the power output regulation is only based on rotor speed control. In this situation, the pitch angle control will be activated to keep DFIG power output at a designated de-loading level. In order to realize the smooth transition between rotor speed control, pitch angle control, and their combination in the de-loading strategy, reasonable division of the WS region according to the operation characteristic of DFIG is important.
The whole WS range is divided into three regions considering the operation characteristic of DFIG, including the start-up region, the FR region, and the cut-out region. In the start-up region, the power output of DFIG is too low for providing frequency support; in the cut-out region, the DFIG will exit the grid as the WS is too high. Therefore, the DFIG will only participate in FR in the FR region. Further, considering the operation limit of DFIG, the FR region should be divided into three regions, that is, the low, medium, and high WS region. The division principle of WS regions is elaborated below.
Firstly, in order to ensure that the up-regulation capacity can be fully obtained, the WS at the lower boundary of low WS region should meet the requirement that the maximum point of DFIG power output can be reached when the rotor speed is at the lower limit. Secondly, when suboptimal rotor speed corresponding to designated de-loading level has reached the upper limit, the WS will be set as the boundary of low and median WS region, and it is where the pitch angle control will start to work. Then, the rated WS is set as the boundary of median and high WS region, where the upper limit of rotor speed equals to the optimal speed; thus, only pitch angle control will work on the de-loading in high WS region. Finally, the cut-out speed of DFIG v cut−out is set as the higher boundary of high WS region. The division principle and corresponding calculated values are shown in Table 1, where d% is the de-loading rate, which is set to be 10% in this paper [5] and can be adjusted according to the actual situation.

WS Region
Division Principle Boundary Value

. Operation Point Selection of Pitch Angle
It is noteworthy that the pitch angle control starts to work at the boundary of low and median WS regions, but the DFIG power output does not increase monotonically with pitch angle here, as shown in Figure 3a, where the current rotor speed keeps at 1.2 p.u. At the boundary WS (9.1 m/s), there are three pitch angle values (β 0 , β 1 and β 2 ) that can keep the blade power coefficient C p to be the value corresponding to (1 − d%)P MPPT (v w ), that is, 0.9C p,max . It can be seen that when λ decreases with the increase in WS in medium WS region, the change rule of C p with pitch angle is gradually changing from non-monotonic to monotonic. In order to keep the power output at (1 − d%)P MPPT (v w ), the pitch angle will be adjusted to keep C p at 0.9C p,max . However, if β 0 is chosen as the initial operation point of pitch angle in medium WS region, the pitch angle will jump near 9.3 m/s due to the change in variation trend of C p , as shown in Figure 3b. decreases with the increase in WS in medium WS region, the change rule of p with pitch angle is gradually changing from non-monotonic to monotonic. In order to keep the power output at     On the one hand, the frequent adjustment of pitch angle is not recommended considering the life-span and maintenance cost of WT [7]; on the other hand, when the time interval of FR order is small, the pitch angle may not be able to achieve the target value as a result of change in rate limit. Therefore, for the sake of avoiding the sharp change in pitch angle, β 2 is chosen as the initial operation point of pitch angle control in medium WS region. Moreover, in low WS region, λ always keep at 7.887 through the adjustment of rotor speed; thus, the profile of p C is the same as 9.1 m/s. In order to avoid the jump of pitch angle at the boundary of low and medium WS regions, the pitch angle should keep operating at β 2 in low WS region, which can achieve the same power output as β 0 , i.e., 0°. On the one hand, the frequent adjustment of pitch angle is not recommended considering the life-span and maintenance cost of WT [7]; on the other hand, when the time interval of FR order is small, the pitch angle may not be able to achieve the target value as a result of change in rate limit. Therefore, for the sake of avoiding the sharp change in pitch angle, β 2 is chosen as the initial operation point of pitch angle control in medium WS region. Moreover, in low WS region, λ always keep at 7.887 through the adjustment of rotor speed; thus, the profile of C p is the same as 9.1 m/s. In order to avoid the jump of pitch angle at the boundary of low and medium WS regions, the pitch angle should keep operating at β 2 in low WS region, which can achieve the same power output as β 0 , i.e., 0 • .

De-Loading Strategy of DFIG
Based on the division of the WS region and operation point selection of pitch angle, the specific de-loading strategy of DFIG in different WS regions are elaborated below.
Low WS region: The DFIG outputs 90% maximum power by controlling the rotor speed to operate at the over-speed suboptimal trajectory, providing reserve capacity of 10% maximum power for up-regulation. The pitch angle remains at β 2 .
Medium WS region: The rotor speed cannot follow the over-speed suboptimal trajectory due to the upper limit of 1.2 p.u.; thus, the pitch angle is increased to keep 90% maximum power output, providing reserve capacity of 10% maximum power for up-regulation. The rotor speed remains at 1.2 p.u.
High WS region: The maximum power point has reached the rated power. The DFIG outputs 90% rated power by controlling the pitch angle, providing reserve capacity of 10% rated power for up-regulation. The rotor speed remains at 1.2 p.u.
In the de-loading strategy, the rotor speed and pitch angle operation strategy for providing FR reserve capacity is summed up in (3) and (4), respectively: where ω FRC and β FRC are the operating rotor speed and a pitch angle of DFIG for providing FR capacity, respectively. In addition, the operation range of the pitch angle is 0-70 • [5].
Therefore, through real-time measured wind speed, the FR reserve capacity of DFIG for up-regulation and down-regulation can be predicted in real-time by (5) and (6), respectively: where P min is the lower limit of DFIG's power output for grid integration, whose value equals the maximum power at cut-in speed. If the power output is lower than P min , DFIG's mechanical power cannot maintain the power balance, thus, is unable to participate in FR. Moreover, R wind up is set to be positive, and R wind down is set to be negative in this paper. The de-loading operation trajectory of DFIG power output is shown in Figure 4.

FR Model of DFIG in Whole WS Range
Based on the combination of rotor speed and pitch angle control, an FR control model is established for DFIG to participate in system frequency up-regulation and down-regulation. When the system frequency fluctuates, the active power output of DFIG is regulated by controlling rotor speed ω r , pitch angle β , or both according to the current WS region. The FR control principle block diagram of DFIG is shown in Figure 5. The active power adjustment through rotor speed control is realized by controlling the active component of rotor current; thus, it can respond to the system frequency change in a relatively fast manner. While the pitch angle control is implemented by controlling the mechanical components of the pitch angle controller, which is much slower [7]. Therefore, the rotor speed control is always given priority in the participation of FR. Based on the operation characteristic and de-loading strategy of DFIG, the FR control methods applied in different WS regions are shown in Table 2. The procedures of DFIG participating in FR in the three WS regions are elaborated below.

FR Model of DFIG in Whole WS Range
Based on the combination of rotor speed and pitch angle control, an FR control model is established for DFIG to participate in system frequency up-regulation and down-regulation. When the system frequency fluctuates, the active power output of DFIG is regulated by controlling rotor speed ω r , pitch angle β, or both according to the current WS region. The FR control principle block diagram of DFIG is shown in Figure 5. The active power adjustment through rotor speed control is realized by controlling the active component of rotor current; thus, it can respond to the system frequency change in a relatively fast manner. While the pitch angle control is implemented by controlling the mechanical components of the pitch angle controller, which is much slower [7]. Therefore, the rotor speed control is always given priority in the participation of FR.
Based on the operation characteristic and de-loading strategy of DFIG, the FR control methods applied in different WS regions are shown in Table 2. The procedures of DFIG participating in FR in the three WS regions are elaborated below. speed ω r , pitch angle β , or both according to the current WS region. The FR control principle block diagram of DFIG is shown in Figure 5. The active power adjustment through rotor speed control is realized by controlling the active component of rotor current; thus, it can respond to the system frequency change in a relatively fast manner. While the pitch angle control is implemented by controlling the mechanical components of the pitch angle controller, which is much slower [7]. Therefore, the rotor speed control is always given priority in the participation of FR. Based on the operation characteristic and de-loading strategy of DFIG, the FR control methods applied in different WS regions are shown in Table 2. The procedures of DFIG participating in FR in the three WS regions are elaborated below.

Low WS region:
Up-regulation: The DFIG reduces ω r so as to increase power output. When ω r reaches ω opt , if the FR task allocated to DFIG by grid control center still cannot be fulfilled, the deceleration should be stopped. Then, β will be reduced to continue increasing power output. If the FR task still cannot be fulfilled even when β reaches 0 • , β stops decreasing and DFIG stays in the MPPT state, providing frequency support in another way.
Down-regulation: The DFIG raises ω r so as to decrease power output. If the FR task cannot be fulfilled even when ω r reaches 1.2 p.u., the acceleration should be stopped. Then, β will be raised to continue decreasing power output. If the FR task still cannot be fulfilled even when β reaches 70 • or the power output reaches P min , β stops increasing, providing frequency support in another way.

Medium WS region:
Up-regulation: The DFIG reduces ω r so as to increase power output. If the FR task still cannot be fulfilled even when ω r reaches ω opt , the deceleration should be stopped. Then, β will be reduced to continue increasing power output. If the FR task still cannot be fulfilled even when β reaches 0 • , β stops decreasing and DFIG stays in the MPPT state, providing frequency support in another way.
Down-regulation: The DFIG raises β so as to decrease power output. If the FR task cannot be fulfilled even when β reaches 70 • or the power output reaches P min , β stops increasing, providing frequency support in another way. ω r always remains at 1.2 p.u.

High WS region:
Up-regulation: The DFIG reduces β so as to increase power output. If the FR task cannot be fulfilled even when the power output reaches P e , β stops decreasing, providing frequency support in another way. ω r always remains at 1.2 p.u.
Down-regulation: The DFIG raises β so as to decrease power output. If the FR task cannot be fulfilled even when β reaches 70 • or the power output reaches P min , β stops increasing, providing frequency support in another way. ω r always remains at 1.2 p.u.

Rotor Speed Protection
In the process of DFIG participating in FR, especially in the low WS region, the excessively large rotor speed drop may cause the DFIG to be unstable or stalled [12]. Therefore, the rotor speed protection is usually introduced to make the DFIG directly exit FR when the rotor speed reaches the lower limit (0.7 p.u.). However, the protection method, with fixed limit value, does not consider the influence of the external WS condition, and DFIG directly exiting FR is also likely to cause a secondary deviation of the system frequency.
Therefore, a rotor speed protection block is added before the DFIG controller in this paper, as shown in Figure 6, where v 1 = 6.7m/s and v 2 = 9.1m/s. The principle of the protection block is that the FR task of DFIG is gradually reduced with the decrease in WS in low WS region so that the DFIG gradually exits FR with the decrease in rotor speed. The input signal of protection trigger is the system FR task, which ensures the synchronous action of rotor speed protection and DFIG's participation in FR. The function of saturation block is set to make the output equal to zero when the input is negative, so the FR task of DFIG will become zero when the WS is below 6.7 m/s, which means that DFIG is totally exiting FR. In this way, the stability of the wind turbine during the FR process is maintained, and the system frequency is prevented from secondary deviation.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 25 In conclusion, based on the de-loading strategy for FR capacity, considering rotor speed control, the secondary FR model for DFIG to participate in both system frequency up-regulation and downregulation in whole WS range is established through the combination of rotor speed and pitch angle control.

Grouping Strategy of EV
As a kind of mobile energy storage device, EV has the advantage of fast response [18]. EVs can participate in system FR through controlling charging and discharging power [19]. The EV aggregator (EVA) collects the declared information of EVs connected to the grid in real time, including the residence time i T in the grid, the state of charge (SOC), and the credibility i ρ . Based on residence time and credibility, EVs connected to the grid are divided into FR service group and energy demand group. In FR service group, according to the initial SOC when EV becomes connected to the grid, EVs are divided into charging group and discharging group. EVA reports the FR capacity of FR service group to control center in real time. According to the system frequency deviation and the FR capacity of DFIG, the grid control center allocates FR task to EVA; then the EVA instructs EVs in the FR service group to charge or discharge to provide frequency support.
Only when the driving demand of vehicle owner is satisfied, EV will participate in FR; otherwise, the EV will be included in the energy demand group, charging to meet the driving demand. The FR service group j Se at time j is shown in (7): where * T is the threshold of residence time in the grid for FR; down ρ is the lower threshold of T is the necessary charging duration for the i-th EV to meet driving demand. After being included in the FR service group, the EV begins to accept the FR instructions from EVA. In order to protect EV battery, frequent conversions between the charging and discharging state, as well as excessive charging/discharging behaviors, should be avoided [20,21]. Therefore, during the FR process, EVs are divided into discharging group DC and charging group C based on the rules shown in Table 3, where  In conclusion, based on the de-loading strategy for FR capacity, considering rotor speed control, the secondary FR model for DFIG to participate in both system frequency up-regulation and down-regulation in whole WS range is established through the combination of rotor speed and pitch angle control.

Grouping Strategy of EV
As a kind of mobile energy storage device, EV has the advantage of fast response [18]. EVs can participate in system FR through controlling charging and discharging power [19]. The EV aggregator (EVA) collects the declared information of EVs connected to the grid in real time, including the residence time T i in the grid, the state of charge (SOC), and the credibility ρ i . Based on residence time and credibility, EVs connected to the grid are divided into FR service group and energy demand group. In FR service group, according to the initial SOC when EV becomes connected to the grid, EVs are divided into charging group and discharging group. EVA reports the FR capacity of FR service group to control center in real time. According to the system frequency deviation and the FR capacity of DFIG, the grid control center allocates FR task to EVA; then the EVA instructs EVs in the FR service group to charge or discharge to provide frequency support.
Only when the driving demand of vehicle owner is satisfied, EV will participate in FR; otherwise, the EV will be included in the energy demand group, charging to meet the driving demand. The FR service group Se j at time j is shown in (7): where T * is the threshold of residence time in the grid for FR; ρ down is the lower threshold of credibility; d i is the planned departure time of the i-th EV; T D i is the necessary charging duration for the i-th EV to meet driving demand.
After being included in the FR service group, the EV begins to accept the FR instructions from EVA. In order to protect EV battery, frequent conversions between the charging and discharging state, as well as excessive charging/discharging behaviors, should be avoided [20,21]. Therefore, during the FR process, EVs are divided into discharging group DC and charging group C based on the rules shown in Table 3, where SOC i,j is the state of charge of the i-th EV at time j; SOC max and SOC min are the upper and lower limits of SOC, respectively. Table 3. Grouping strategy for EVs (electric vehicles) in FR (frequency regulation) service group.

Group Mode Description
Discharging group 1 The EV is included in discharging group. If SOC i, j > SOC min , the EV remains in discharging group and discharges for FR. 2 The EV is included in the charging group. If SOC i, j ≥ SOC max , the EV switches into the discharging group and discharges for FR.
Charging group 3 The EV is included in the charging group. If SOC i,j < SOC max , the EV remains in charging group and charges for FR. 4 The EV is included in discharging group. If SOC i, j ≤ SOC min , the EV switches into charging group and charges for FR.
Therefore, the discharging and charging groups at time j are shown in (8) and (9), respectively: where P i,j is the charging/discharging power for FR of the i-th EV at time j.
The EV quantity in the discharging and charging groups at time j are shown in (10) and (11), respectively:

Real-Time Prediction of EVA's FR Capacity
As a kind of transportation, EV's traffic characteristics cannot be ignored. It is possible for EV owners to violate dispatching instructions and leaving the grid before the declared departure time, resulting in the inaccurate prediction of EV's FR capacity. When the number of EVs violating the FR dispatching plan reaches a certain level, the FR task allocated to the EVA from the grid cannot be effectively accomplished.
Therefore, in order to solve this problem, the concept of credibility is introduced into the prediction of EV's FR capacity [22]. If an EV obeys the declared plan strictly, leaving the grid according to the declared time, the credibility of this EV is regarded as high; thus, its FR capacity will be fully included in the FR capacity of EVA. However, if an EV is used to suddenly leaving the grid during the declared residence time, the credibility of this EV is regarded as low, and its FR capacity will be included only partially in EVA's FR capacity with the proportion based on the value of credibility. The definition of credibility is shown in (12): where T a i is the historical average actual FR participation time of the EV; T d i is historical average declared FR participation time of the EV; ρ up is the upper threshold of credibility. When the credibility is too high/low, the value is set to be 1/0.
Taking the credibility of EVs into account, the FR capacity of EVA at time j can be predicted in real-time by (13)-(15): where R ev up,j is EVA's real-time predicted capacity for up-regulation; R ev down,j is EVA's real-time predicted capacity for down-regulation; m is the limit of EV's charging/discharging power; R ev j is the real-time predicted FR capacity of EVA; ∆ f is the system frequency deviation.

The Coordinated FR Strategy of DFIG and EV
In order to improve the reliability of frequency support provided by DFIGs, considering the frequency stability improvement from load side, EV is used as an auxiliary energy storage system to coordinate with DFIG to participate in system FR.

Real-Time Allocation of System FR Task
In order to realize the coordination of DFIG and EVA in system FR, a real-time allocation strategy of system FR task is proposed. The de-loading operation will influence the economic efficiency of DFIG [4]. Therefore, taking into account the economical operation of the wind farm, the DFIG is set to have priority over EV in up-regulation. Meanwhile, when the EV is connected to the grid, charging is the first demand for EV owner. In order to reduce EV's charging load for driving demand, the EVA is set to have priority over DFIG in down-regulation.
When the system frequency fluctuates, the system FR task P task will be allocated to DFIG and EVA according to the FR capacity and the type of FR, which includes five modes, as elaborated below. P task is set to be positive for up-regulation and negative for down-regulation.
Mode 1: If P task > 0, the system FR task P task will be allocated to DFIG first. If the up-regulation capacity of DFIG is more than the system FR task, i.e., R wind up ≥ P task , all the FR power output will be provided by DFIG, while EVA will not participate in FR. Thus, the FR tasks of DFIG and EVA are shown in (16): P wind task = P task P ev task = 0 where P wind task and P ev task are the FR tasks allocated to DFIG and EVA, respectively. Mode 2: If P task > 0, P task will be allocated to DFIG first. If R wind up < P task , the DFIG cannot accomplish the system FR task alone, so the shortage will be allocated to EVA. Thus, the FR tasks of DFIG and EVA are shown in (17): Mode 3: If P task < 0, P task will be allocated to EVA first. If the up-regulation capacity of EVA is more than the system FR task, i.e., R ev down ≤ P task , all the FR power output will be provided by EVA, while the DFIG will not participate in FR. Thus, the FR tasks of DFIG and EVA are shown in (18): Mode 4: If P task < 0, P task will be allocated to EVA first. If R ev down > P task , the EVA cannot accomplish the system FR task alone, so the shortage will be allocated to DFIG. Thus, the FR tasks of DFIG and EVA are shown in (19): Mode 5: If P task = 0, DFIG and EVA will not participate in FR. Thus, the FR tasks of DFIG and EVA are shown in (20):

The Procedure of Coordinated FR Strategy
The detailed procedure of DFIG and EVA coordinating to participate in system secondary FR are shown below.
Step 1: According to the WS environment, whether the DFIG can be connected to the grid safely or not will be confirmed. If the security of the DFIG connection can be ensured, the current WS region is determined by the WS, and then the DFIG is controlled to operate in the corresponding de-loading state. At the same time, EVs are connected to the grid according to daily driving rules. Under the dispatch of EVA, EVs in the grid are divided into discharging group, charging group, and energy demand group; then, EV's controllable energy for up-regulation and down-regulation can be calculated in real time.
Step 2: Considering the FR capacity and the type of FR, the system FR task P task will be allocated to DFIG and EVA according to the proposed real-time allocation strategy.
Step 3: Following the FR task allocated by the control center, both the DFIG and EVA output power to participate in FR.
DFIG: Based on the allocated FR task, the DFIG starts to output power by following the commands from the rotor speed and pitch angle controller. When P wind task is accomplished, DFIG stops regulating rotor speed and pitch angle, and the actual power output of DFIG is P wind a . If P wind task has not been accomplished when the rotor speed, pitch angle, or output power has reached the operation limit, the DFIG stops regulating its operation state and remains at current power output P wind a . EVA: Based on the allocated FR task, the EVA instructs the EVs inside the cluster to participate in FR. If the FR capacity of EVA is larger than the FR task of EVA, i.e., |R ev | ≥ P task ev , EVs in the FR service group will be sorted by the credibility ρ i of each EV, and then the EV with a lower credibility will be given higher priority to output power for FR until the FR task is accomplished; if |R ev | < P task ev , all the EVs in the FR service group will output power for FR at the same time. The actual power output of EVA is P ev a .
Step 4: If the system FR task can be accomplished by the actual FR power output of DFIG and EVA, i.e., P wind a + P ev a ≥ |P task |, the conventional unit will not participate in FR. If the actual FR power output of DFIG and EVA cannot fulfill the system FR task, i.e., P wind a + P ev a < |P task |, the shortage will be compensated by the power output of a conventional unit. The FR task allocated to the conventional unit is shown in (21): Then, considering the limitation of ramp rate, the conventional unit output power according to the allocated FR task.
In summary, the flow chart of DFIG and EVA coordinating to participate in system FR is shown in Figure 7, where the up-regulation and down-regulation capacity are set to be positive and negative, respectively, for both EVA and DFIG.
Then, considering the limitation of ramp rate, the conventional unit output power according to the allocated FR task.
In summary, the flow chart of DFIG and EVA coordinating to participate in system FR is shown in Figure 7, where the up-regulation and down-regulation capacity are set to be positive and negative, respectively, for both EVA and DFIG.

Simulation Model and Parameters
The simulation cases in this paper are based on a two-area interconnected system, which contains DFIGs, EV cluster, conventional units, and random loads, as shown in Figure 8

Simulation Model and Parameters
The simulation cases in this paper are based on a two-area interconnected system, which contains DFIGs, EV cluster, conventional units, and random loads, as shown in Figure 8. All the simulation results are based on the data of Area 1. In the system model, M and D are the inertia constant and damping coefficient of the generator, respectively; 1/(s + T delay ) represents the delay effect of system control and communication, and T delay is called the system time delay constant. The load frequency control (LFC) blocks in two areas adopt the combination mode of tie-line bias control (TBC), i.e., TBC-TBC mode [23], and the system FR task signal is obtained by proportional-integral (PI) control. The conventional unit block adopts the model of traditional thermal power unit [24].
The parameters of the system model are listed in Table 4 [25]. The parameters of DFIG are listed in Table 5 [26], and the model of DFIGs in two areas are the same. The parameters of EV are listed in Table 6. The driving behaviors and SOC changes of individual EVs are simulated by the Monte Carlo method. By reference to Beijing transport annual report [27], it is assumed that the time of arrival at workplace and time of departure from workplace obey the normal distribution, and the driving distance obeys the lognormal distribution. The system rated frequency f 0 is 50 Hz. In the simulation results, the base value for the calculation of per-unit value is set to be 5 MW.
In the system model, the load fluctuation and WS are simulated through white noise by reference to Ref. [14]. The profiles of load fluctuation and WS in two areas are shown in Figure 9a,b, respectively. In the figure, "p.u." is the abbreviation for per unit, which means the variable is represented in the form of per-unit value. The parameters of the system model are listed in Table 4 [25]. The parameters of DFIG are listed in Table 5 [26], and the model of DFIGs in two areas are the same. The parameters of EV are listed in Table 6. The driving behaviors and SOC changes of individual EVs are simulated by the Monte Carlo method. By reference to Beijing transport annual report [27], it is assumed that the time of arrival at workplace and time of departure from workplace obey the normal distribution, and the driving distance obeys the lognormal distribution. The system rated frequency 0 f is 50 Hz. In the simulation results, the base value for the calculation of per-unit value is set to be 5 MW.      Cut-out speed 25 m/s In the system model, the load fluctuation and WS are simulated through white noise by reference to Ref. [14]. The profiles of load fluctuation and WS in two areas are shown in Figure 9a,b, respectively. In the figure, "p.u." is the abbreviation for per unit, which means the variable is represented in the form of per-unit value.

DFIG Participating in FR
In order to verify the effectiveness of the proposed FR model for DFIG, two scenarios are set up. In Scenario 1, DFIG participates in FR, and the lack of FR demand power is compensated by the conventional unit. In Scenario 2, DFIG does not participate in FR and operates according to the proposed de-loading strategy. Other conditions of the two scenarios are the same.
The power output of DFIG in two scenarios are shown in Figure 10. Following the proposed de-loading operation strategy for FR capacity, DFIG's power output in Scenario 2 varies with the change in WS. It can be seen from the figure that DFIG's power output in Scenario 1 fluctuates around the de-loading operation curve in Scenario 2, which indicates that DFIG plays an active role in system FR through regulating power output.
proposed de-loading strategy. Other conditions of the two scenarios are the same.
The power output of DFIG in two scenarios are shown in Figure 10. Following the proposed deloading operation strategy for FR capacity, DFIG's power output in Scenario 2 varies with the change in WS. It can be seen from the figure that DFIG's power output in Scenario 1 fluctuates around the de-loading operation curve in Scenario 2, which indicates that DFIG plays an active role in system FR through regulating power output. The FR task of DFIG in Scenario 1 is shown in Figure 11, where the blue curve represents the FR task without rotor speed protection, while the red curve represents the FR task after being processed by rotor speed protection block. It can be seen that in the time period of 2:00-5:00 and 10:30-14:30, the WS is in the low-speed region, as shown in Figure 9b; thus, the rotor speed protection will work and reduce DFIG's FR task with the decrease in WS. In the time period of 0:00-2:00 and 5:00-7:00, the WS is lower than 6.7 m/s, so DFIG completely exits FR, and its FR task is set to be zero. The FR capacity of DFIG for up-regulation and down-regulation are shown in Figure 12a,b, respectively. It can be seen that in low and medium WS region, the up-regulation capacity of DFIG is relatively less when the WS is small and rises with the increase in WS. When the WS enters high WS region, the up-regulation capacity remains at 0.1 p.u. At the same time, the down-regulation capacity of DFIG rises as well with the increase in WS, and DFIG's down-regulation capacity is more sufficient than up-regulation capacity due to the operation characteristic of DFIG and the set value of deloading rate. The FR task of DFIG in Scenario 1 is shown in Figure 11, where the blue curve represents the FR task without rotor speed protection, while the red curve represents the FR task after being processed by rotor speed protection block. It can be seen that in the time period of 2:00-5:00 and 10:30-14:30, the WS is in the low-speed region, as shown in Figure 9b; thus, the rotor speed protection will work and reduce DFIG's FR task with the decrease in WS. In the time period of 0:00-2:00 and 5:00-7:00, the WS is lower than 6.7 m/s, so DFIG completely exits FR, and its FR task is set to be zero.
proposed de-loading strategy. Other conditions of the two scenarios are the same.
The power output of DFIG in two scenarios are shown in Figure 10. Following the proposed deloading operation strategy for FR capacity, DFIG's power output in Scenario 2 varies with the change in WS. It can be seen from the figure that DFIG's power output in Scenario 1 fluctuates around the de-loading operation curve in Scenario 2, which indicates that DFIG plays an active role in system FR through regulating power output. The FR task of DFIG in Scenario 1 is shown in Figure 11, where the blue curve represents the FR task without rotor speed protection, while the red curve represents the FR task after being processed by rotor speed protection block. It can be seen that in the time period of 2:00-5:00 and 10:30-14:30, the WS is in the low-speed region, as shown in Figure 9b; thus, the rotor speed protection will work and reduce DFIG's FR task with the decrease in WS. In the time period of 0:00-2:00 and 5:00-7:00, the WS is lower than 6.7 m/s, so DFIG completely exits FR, and its FR task is set to be zero. The FR capacity of DFIG for up-regulation and down-regulation are shown in Figure 12a,b, respectively. It can be seen that in low and medium WS region, the up-regulation capacity of DFIG is relatively less when the WS is small and rises with the increase in WS. When the WS enters high WS region, the up-regulation capacity remains at 0.1 p.u. At the same time, the down-regulation capacity of DFIG rises as well with the increase in WS, and DFIG's down-regulation capacity is more sufficient than up-regulation capacity due to the operation characteristic of DFIG and the set value of deloading rate. The FR capacity of DFIG for up-regulation and down-regulation are shown in Figure 12a,b, respectively. It can be seen that in low and medium WS region, the up-regulation capacity of DFIG is relatively less when the WS is small and rises with the increase in WS. When the WS enters high WS region, the up-regulation capacity remains at 0.1 p.u. At the same time, the down-regulation capacity of DFIG rises as well with the increase in WS, and DFIG's down-regulation capacity is more sufficient than up-regulation capacity due to the operation characteristic of DFIG and the set value of de-loading rate.
The FR output of DFIG is shown in Figure 13. The comparison of FR output and FR task is shown in Figure 13a. It can be seen that the FR task is well accomplished in the down-regulation situation since the FR capacity for down-regulation is sufficient, but the FR task for up-regulation is not completely accomplished. This is because, in the situation of up-regulation, the FR output of DFIG is restrained by the FR capacity, as shown in Figure 13b.
The operation state changes of DFIG are shown in Figure 14, where the rotor speed or pitch angle equaling to zero indicates that the DFIG has exited FR. The rotor speed change in the FR process is shown in Figure 14a. Corresponding to the WS profile, it can be seen that in low WS region, ω FRC in Scenario 2 changes to keep the power output at 0.9P MPPT (v w ). On this basis, ω r in Scenario 1 changes to regulate DFIG's power output for FR. In medium and high WS regions, ω FRC remains at 1.2 p.u., and ω r only decrease to participate in up-regulation except for transient over-limit. Change in β in the FR process is shown in Figure 14b. In low WS region, β FRC in Scenario 2 keeps at β 2 = 2.3 • ; in medium and high WS regions, β FRC will increase to keep power output at 0.9P MPPT (v w ) and 0.9P e , respectively. In Scenario 1, β changes to participate in both up-regulation and down-regulation. The FR output of DFIG is shown in Figure 13. The comparison of FR output and FR task is shown in Figure 13a. It can be seen that the FR task is well accomplished in the down-regulation situation since the FR capacity for down-regulation is sufficient, but the FR task for up-regulation is not completely accomplished. This is because, in the situation of up-regulation, the FR output of DFIG is restrained by the FR capacity, as shown in Figure 13b.  The operation state changes of DFIG are shown in Figure 14, where the rotor speed or pitch angle equaling to zero indicates that the DFIG has exited FR. The rotor speed change in the FR process The FR output of DFIG is shown in Figure 13. The comparison of FR output and FR task is shown in Figure 13a. It can be seen that the FR task is well accomplished in the down-regulation situation since the FR capacity for down-regulation is sufficient, but the FR task for up-regulation is not completely accomplished. This is because, in the situation of up-regulation, the FR output of DFIG is restrained by the FR capacity, as shown in Figure 13b.  The operation state changes of DFIG are shown in Figure 14, where the rotor speed or pitch angle equaling to zero indicates that the DFIG has exited FR. The rotor speed change in the FR process 1.2 p.u., and ω r only decrease to participate in up-regulation except for transient over-limit. Change in β in the FR process is shown in Figure 14b. In low WS region, β FRC in Scenario 2 keeps at β =°2 2 .3 ; in medium and high WS regions, β FRC will increase to keep power output at The comparison of the FR effect in the two scenarios is shown in Figure 15. The root mean square (RMS) value of frequency deviation in Scenario 1 and Scenario 2 is 0.0346 Hz and 0.0386 Hz, respectively. It can be seen that the FR effect of DFIG participating in FR is obviously better than the situation without frequency support by DFIG, especially in medium and high WS regions. In low WS region, the FR effect is relatively worse as the FR output of DFIG is restrained by both the rotor speed protection and limit of FR capacity. In conclusion, the proposed FR model of DFIG can suppress system frequency deviation effectively under continuously varying WS.  The comparison of the FR effect in the two scenarios is shown in Figure 15. The root mean square (RMS) value of frequency deviation in Scenario 1 and Scenario 2 is 0.0346 Hz and 0.0386 Hz, respectively. It can be seen that the FR effect of DFIG participating in FR is obviously better than the situation without frequency support by DFIG, especially in medium and high WS regions. In low WS region, the FR effect is relatively worse as the FR output of DFIG is restrained by both the rotor speed protection and limit of FR capacity. In conclusion, the proposed FR model of DFIG can suppress system frequency deviation effectively under continuously varying WS.
1.2 p.u., and ω r only decrease to participate in up-regulation except for transient over-limit. Change in β in the FR process is shown in Figure 14b. In low WS region, β FRC in Scenario 2 keeps at β =°2 2 .3 ; in medium and high WS regions, β FRC will increase to keep power output at The comparison of the FR effect in the two scenarios is shown in Figure 15. The root mean square (RMS) value of frequency deviation in Scenario 1 and Scenario 2 is 0.0346 Hz and 0.0386 Hz, respectively. It can be seen that the FR effect of DFIG participating in FR is obviously better than the situation without frequency support by DFIG, especially in medium and high WS regions. In low WS region, the FR effect is relatively worse as the FR output of DFIG is restrained by both the rotor speed protection and limit of FR capacity. In conclusion, the proposed FR model of DFIG can suppress system frequency deviation effectively under continuously varying WS.

DFIG and EV Coordinating to Participate in FR
In order to verify the effectiveness of the coordinated FR strategy of DFIG and EV, Scenario 3 is set up. In Scenario 3, DFIG and EV cluster coordinates to participate in system FR, and the lack of FR demand power is compensated by the conventional unit. Other conditions are the same as with Scenario 1 and Scenario 2. The system FR task in Scenario 3 is shown in Figure 16, which will be allocated to EV and DFIG, respectively, according to the proposed real-time allocation strategy.
In order to verify the effectiveness of the coordinated FR strategy of DFIG and EV, Scenario 3 is set up. In Scenario 3, DFIG and EV cluster coordinates to participate in system FR, and the lack of FR demand power is compensated by the conventional unit. Other conditions are the same as with Scenario 1 and Scenario 2. The system FR task in Scenario 3 is shown in Figure 16, which will be allocated to EV and DFIG, respectively, according to the proposed real-time allocation strategy. The FR capacity of EV cluster for up-regulation and down-regulation is shown in Figure 17a,b, respectively. With EVs accessing and leaving the grid according to the driving rules, the FR capacity of EV cluster will vary with time. The FR output of EV cluster for up-regulation and down-regulation is shown in Figure 18a,b, respectively. It can be seen from the figure that as EVA is given priority in down-regulation in the proposed FR task allocation strategy, EVA's down-regulation task is more than up-regulation task. EV cluster will participate in up-regulation only when the up-regulation capacity of EVA cannot fulfill the FR task. In 9:00-10:00 of Figure 18a as well as 8:30-9:00 and 17:30-19:00 of Figure 18b, the FR output of EV cluster is constrained by the FR capacity.   The FR capacity of EV cluster for up-regulation and down-regulation is shown in Figure 17a,b, respectively. With EVs accessing and leaving the grid according to the driving rules, the FR capacity of EV cluster will vary with time. The FR output of EV cluster for up-regulation and down-regulation is shown in Figure 18a,b, respectively. It can be seen from the figure that as EVA is given priority in down-regulation in the proposed FR task allocation strategy, EVA's down-regulation task is more than up-regulation task. EV cluster will participate in up-regulation only when the up-regulation capacity of EVA cannot fulfill the FR task. In 9:00-10:00 of Figure 18a as well as 8:30-9:00 and 17:30-19:00 of Figure 18b, the FR output of EV cluster is constrained by the FR capacity.
The FR output of DFIG is shown in Figure 19. It can be seen from the figure that as DFIG is given priority in up-regulation in the proposed FR task allocation strategy, DFIG's up-regulation task is more than down-regulation task. DFIG will participate in down-regulation only when the down-regulation capacity of EVA cannot fulfill the FR task. The FR task allocated to DFIG is almost accomplished, except for some up-regulation moment that the FR output of DFIG is constrained by the FR capacity.
If the joint FR output of EV cluster and DFIG cannot fulfill the system FR task, the conventional unit will output power to compensate for the lack of FR demand power. The FR task allocated to the conventional unit is shown in Figure 20.

DFIG and EV Coordinating to Participate in FR
In order to verify the effectiveness of the coordinated FR strategy of DFIG and EV, Scenario 3 is set up. In Scenario 3, DFIG and EV cluster coordinates to participate in system FR, and the lack of FR demand power is compensated by the conventional unit. Other conditions are the same as with Scenario 1 and Scenario 2. The system FR task in Scenario 3 is shown in Figure 16, which will be allocated to EV and DFIG, respectively, according to the proposed real-time allocation strategy. The FR capacity of EV cluster for up-regulation and down-regulation is shown in Figure 17a,b, respectively. With EVs accessing and leaving the grid according to the driving rules, the FR capacity of EV cluster will vary with time. The FR output of EV cluster for up-regulation and down-regulation is shown in Figure 18a,b, respectively. It can be seen from the figure that as EVA is given priority in down-regulation in the proposed FR task allocation strategy, EVA's down-regulation task is more than up-regulation task. EV cluster will participate in up-regulation only when the up-regulation capacity of EVA cannot fulfill the FR task. In 9:00-10:00 of Figure 18a as well as 8:30-9:00 and 17:30-19:00 of Figure 18b, the FR output of EV cluster is constrained by the FR capacity.   The FR output of DFIG is shown in Figure 19. It can be seen from the figure that as DFIG is given priority in up-regulation in the proposed FR task allocation strategy, DFIG's up-regulation task is more than down-regulation task. DFIG will participate in down-regulation only when the downregulation capacity of EVA cannot fulfill the FR task. The FR task allocated to DFIG is almost accomplished, except for some up-regulation moment that the FR output of DFIG is constrained by the FR capacity. If the joint FR output of EV cluster and DFIG cannot fulfill the system FR task, the conventional unit will output power to compensate for the lack of FR demand power. The FR task allocated to the conventional unit is shown in Figure 20.   The FR output of DFIG is shown in Figure 19. It can be seen from the figure that as DFIG is given priority in up-regulation in the proposed FR task allocation strategy, DFIG's up-regulation task is more than down-regulation task. DFIG will participate in down-regulation only when the downregulation capacity of EVA cannot fulfill the FR task. The FR task allocated to DFIG is almost accomplished, except for some up-regulation moment that the FR output of DFIG is constrained by the FR capacity. If the joint FR output of EV cluster and DFIG cannot fulfill the system FR task, the conventional unit will output power to compensate for the lack of FR demand power. The FR task allocated to the conventional unit is shown in Figure 20.  In Scenario 1, DFIG coordinates with the conventional unit to participate in FR, while in Scenario 3, DFIG coordinates with EVA to participate in FR. The comparison of the FR effect in Scenario 1 and Scenario 3 is shown in Figure 21. The RMS value of frequency deviation in Scenario 1 and Scenario 3 is 0.0346 Hz and 0.0267 Hz, respectively. It can be seen that the FR effect of DFIG coordinating with EVA is obviously better than the situation when DFIG is coordinating with the conventional unit. This is because EVs can respond to the FR signal more quickly than the conventional unit, so the system FR task can be better accomplished. Therefore, through the proposed coordinated FR strategy, EV cluster is able to assist DFIG to participate in FR effectively.  In Scenario 1, DFIG coordinates with the conventional unit to participate in FR, while in Scenario 3, DFIG coordinates with EVA to participate in FR. The comparison of the FR effect in Scenario 1 and Scenario 3 is shown in Figure 21. The RMS value of frequency deviation in Scenario 1 and Scenario 3 is 0.0346 Hz and 0.0267 Hz, respectively. It can be seen that the FR effect of DFIG coordinating with EVA is obviously better than the situation when DFIG is coordinating with the conventional unit. This is because EVs can respond to the FR signal more quickly than the conventional unit, so the system FR task can be better accomplished. Therefore, through the proposed coordinated FR strategy, EV cluster is able to assist DFIG to participate in FR effectively. In Scenario 1, DFIG coordinates with the conventional unit to participate in FR, while in Scenario 3, DFIG coordinates with EVA to participate in FR. The comparison of the FR effect in Scenario 1 and Scenario 3 is shown in Figure 21. The RMS value of frequency deviation in Scenario 1 and Scenario 3 is 0.0346 Hz and 0.0267 Hz, respectively. It can be seen that the FR effect of DFIG coordinating with EVA is obviously better than the situation when DFIG is coordinating with the conventional unit. This is because EVs can respond to the FR signal more quickly than the conventional unit, so the system FR task can be better accomplished. Therefore, through the proposed coordinated FR strategy, EV cluster is able to assist DFIG to participate in FR effectively.
The system frequency deviation (Hz) Figure 21. Comparison of the FR (frequency regulation) Effect in Scenario 1 and Scenario 3.
The proposed coordinated FR strategy is based on the real-time allocation of system FR task, which considers the FR capacity of DFIG and EV cluster to make the allocation more reasonable. In order to verify the superiority of this method, Scenario 4 is set up. In Scenario 4, the system FR task is allocated to the DFIG and EV cluster at a fixed ratio of 5:5. The FR task accomplishment situation in the two scenarios is shown in Figure 22. If the joint FR output of DFIG and EV cannot fulfill the system FR task, the shortage will be allocated to the conventional unit, and the FR task of the conventional unit is shown in Figure 22a. It can be seen that in Scenario 4, the conventional unit's FR task is obviously more than Scenario 3. This is because in the FR strategy, with fixed allocation ratio, when the FR capacity is sufficient, EV cluster and DFIG could have provided larger FR output, but the FR task is restrained by fixed allocation ratio. Due to the high climbing rate of the conventional unit, the FR task allocated to the conventional unit cannot be well accomplished, so there exists the shortage of the system FR demand power, as shown in Figure 22b. The RMS value of FR demand power shortage in Scenario 3 and Scenario 4 is 0.0325 and 0.0342 p.u., respectively. Therefore, the shortage of the system FR demand power in Scenario 4 is larger, resulting in worse FR effect. The proposed coordinated FR strategy is based on the real-time allocation of system FR task, which considers the FR capacity of DFIG and EV cluster to make the allocation more reasonable. In order to verify the superiority of this method, Scenario 4 is set up. In Scenario 4, the system FR task is allocated to the DFIG and EV cluster at a fixed ratio of 5:5. The FR task accomplishment situation in the two scenarios is shown in Figure 22. If the joint FR output of DFIG and EV cannot fulfill the system FR task, the shortage will be allocated to the conventional unit, and the FR task of the conventional unit is shown in Figure 22a. It can be seen that in Scenario 4, the conventional unit's FR task is obviously more than Scenario 3. This is because in the FR strategy, with fixed allocation ratio, when the FR capacity is sufficient, EV cluster and DFIG could have provided larger FR output, but the FR task is restrained by fixed allocation ratio. Due to the high climbing rate of the conventional unit, the FR task allocated to the conventional unit cannot be well accomplished, so there exists the shortage of the system FR demand power, as shown in Figure 22b. The RMS value of FR demand power shortage in Scenario 3 and Scenario 4 is 0.0325 and 0.0342 p.u., respectively. Therefore, the shortage of the system FR demand power in Scenario 4 is larger, resulting in worse FR effect. The comparison of the FR effect in Scenario 3 and Scenario 4 is shown in Figure 23. The RMS value of system frequency deviation in Scenario 3 and Scenario 4 is 0.2671 Hz and 0.3595 Hz, respectively. It can be seen that the FR effect of proposed real-time allocation strategy is better than The comparison of the FR effect in Scenario 3 and Scenario 4 is shown in Figure 23. The RMS value of system frequency deviation in Scenario 3 and Scenario 4 is 0.2671 Hz and 0.3595 Hz, respectively. It can be seen that the FR effect of proposed real-time allocation strategy is better than the situation of fixed allocation ratio. Therefore, the proposed coordinated FR strategy with the real-time allocation of system FR task is more reasonable and effective than the FR strategy with fixed allocation ratio.
(b) Figure 22. FR (frequency regulation) task accomplishment situation in Scenario 3 and Scenario 4. (a) The FR task of a conventional unit; (b) The shortage of system FR demand power.
The comparison of the FR effect in Scenario 3 and Scenario 4 is shown in Figure 23. The RMS value of system frequency deviation in Scenario 3 and Scenario 4 is 0.2671 Hz and 0.3595 Hz, respectively. It can be seen that the FR effect of proposed real-time allocation strategy is better than the situation of fixed allocation ratio. Therefore, the proposed coordinated FR strategy with the realtime allocation of system FR task is more reasonable and effective than the FR strategy with fixed allocation ratio.

Conclusions
With the increase in wind power penetration in the power system, WTs are required to have the ability to control active power output. In order to compensate for DFIG's weakness in FR caused by the indeterminacy and volatility of wind power, an electric vehicle is adopted as an auxiliary storage device to assist WTs to participate in system FR. This study establishes the FR model of DFIG under continuously varying WS in whole WS range, then proposes a coordinated control strategy of DFIG and EV cluster for system secondary FR with the real-time allocation of system FR task. The effectiveness of the proposed FR strategy is verified through various simulation cases based on a two-

Conclusions
With the increase in wind power penetration in the power system, WTs are required to have the ability to control active power output. In order to compensate for DFIG's weakness in FR caused by the indeterminacy and volatility of wind power, an electric vehicle is adopted as an auxiliary storage device to assist WTs to participate in system FR. This study establishes the FR model of DFIG under continuously varying WS in whole WS range, then proposes a coordinated control strategy of DFIG and EV cluster for system secondary FR with the real-time allocation of system FR task. The effectiveness of the proposed FR strategy is verified through various simulation cases based on a two-area interconnected system. Through the analysis of simulation results, the following conclusion can be obtained:

1.
Considering the smooth transition of different control methods, the de-loading strategy enables DFIG to have both up-regulation and down-regulation capacity in every FR region, and the FR capacity increases with WS in low and medium WS regions while remains at the steady value in high WS region.

2.
Through the FR control model based on the combination of rotor speed and pitch angle control, DFIG can effectively suppress system frequency deviation under continuously varying WS in the FR region.

3.
Through the coordinated FR strategy of DFIG and EV, EV cluster can assist DFIG to participate in system secondary FR effective in the premise of meeting driver's demand, and the FR effect is better than DFIG coordinating with the conventional unit. 4.
The real-time allocation strategy of system FR task, considering the FR capacity of DFIG and EV cluster, is reasonable and ensures good FR effect, which is better than the situation of fixed allocation ratio.