A New Short-Circuit Current Calculation and Fault Analysis Method Suitable for Doubly Fed Wind Farm Groups

The transient characteristics of wind farms in groups are quite different; in addition, there is a strong coupling between the wind farms and the grid, and these factors make the fault analysis of the grid with wind farm groups complicated. In order to solve this problem, a mathematical model of the converter is established based on the input-output external characteristics of the converter, and a transient model of a doubly fed wind turbine (DFIG) is presented considering the influence of the low-voltage ride-through control (LVRT) of the converter, and the effect mechanism of the LVRT strategy on the short-circuit current is analyzed. Finally, a short-circuit current calculation model of a doubly fed wind turbine with low-voltage crossing control is established. The interaction mechanism between wind farms during the fault is analyzed, and a short-circuit current calculation method of doubly fed wind farm groups is proposed. RTDS is used to verify the accuracy of the proposed short-circuit current calculation method for doubly fed field groups. On this basis, a method of power grid fault analysis after doubly fed field group access is discussed and analyzed.


Introduction
Due to the characteristics of the regional distribution of wind resources in China, wind power generation is often built in two forms: cluster development and centralized grid connection [1].At present, China is planning to build nine ten-million-kilowatt wind power bases [2].If a wind farm group is connected to the grid, it is necessary to carry out corresponding research to meet the needs of large-scale grid connection of new energy [3].
Some field data show that in some areas where wind farms are clustered, wind farms cause large fault currents [4].Therefore, considering that the accuracy of short-circuit current calculation also has a great influence on the relay protection action, it is necessary to propose a new method for calculating short-circuit current.
At present, the existing literature has focused on the problems caused by the integration of wind farm groups into the grid [5].Some literature [6] analyzes the impact of large-capacity wind farms on the transient stability of the power grid [7].
Most of the existing research considers wind farm groups as equivalent to a wind turbine of equal capacity [8].However, the centralized grid-connected wind farms are distributed according to the characteristics of wind resources, the actual distance between farms is relatively far, and the operating conditions of each wind farm are different, which makes their transient characteristics after fault greatly different [9].In addition, there is a strong coupling relationship between each wind farm and the system [10].Therefore, it is significant to propose a method to calculate the short-circuit current considering the coupling relationship between wind farm groups and the system [11].
In addition, in order to meet the requirement of providing a certain voltage support during a fault [12].The converter of the DFIG cannot simply be blocked using crowbar during the fault period, and the rotor excitation needs to be adjusted through the LVRT Sensors 2023, 23, 8372 2 of 16 control strategy to output reactive current [13].In this case, the output short-circuit current of the DFIG is affected by the converter control characteristics [14].Therefore, it is necessary to consider the low-voltage crossing control strategy and the control characteristics of the converter when calculating the short-circuit current [15].
In this paper, a transient model of doubly fed wind turbines is proposed, which takes into account the influence of the low-voltage crossing control strategy [16].Considering the transient model and the influence of converter control characteristics, a DFIG shortcircuit current calculation model is established [17].The coupling relationship between the short-circuit current of the wind farm and the voltage of the grid node is further given, and the interaction mechanism of the short-circuit current within the wind farm groups and between the wind farm and the grid system is analyzed; On this basis, this paper proposes a calculation method of short-circuit current suitable for wind farm group access [18].Secondly, in order to verify the accuracy of the proposed short-circuit current calculation method, a doubly fed field group experimental platform is established through RTDS.Finally, on the basis of further analysis of a doubly fed wind farm equivalent model, this paper proposes a fault analysis method for a group of doubly fed wind farms connected to the electricity grid.

DFIG Transient Model Considering the Control Characteristics of the Converter
In previous studies, it was believed that after a fault occurred, crowbar was switched on and the rotor converter was blocked [19].According to wind power grid-connected specifications, the converter needs to emit a reactive current to support grid voltage recovery during low-voltage traversal [20].In addition, due to the transient process of a DFIG, there will be a certain effect because of the rotor-side converter atresia.Therefore, the converter input and output characteristics should also be considered in the analysis of the electromagnetic transient process.
It has been shown that power imbalances between converters can be suppressed by using the chopper circuit in the DC bus [21,22].Thus, the DC bus voltage during the fault period is guaranteed to be near the rated value [23].Therefore, it is assumed in this paper that DC voltage u dc can remain stable after failure, DC voltage u dc can be obtained using the rotor-side converter, and the excitation voltage u r can be adjusted through the converter duty ratio.Then, the rotor-side excitation voltage u r can be expressed as follows: where A, n, and K dc are the current inverting coefficient, ratio, and duty cycle, respectively; i dc , X rsc , and ∆u rsc refer to the DC bus current, equivalent arcing reactance of converter, and IGBT voltage drop, respectively.In the transient process, the magnetic saturation phenomenon is not considered and the rotation speed is assumed to be constant, then motor convention is adopted [24].The DFIG space vector model under a synchronous rotation coordinate system is as follows [25,26]: Equation (2) gives the relationship between stator voltage u s and rotor voltage u r after converting the stator side; Equation (3) gives the relationship between the DQ-axis flux ψ s and ψ r and the stator and rotor currents.
Sensors 2023, 23, 8372 gives the relationship between stator voltage us and rotor voltage ur converting the stator side; Equation (3) gives the relationship between the DQ-axis flu and ψr and the stator and rotor currents.
The synthesis of Equations ( 1)-( 3) is shown in Figure 1.Here Eg is the equivalent potential at the grid side, Z1L is the equivalent imped from the system to the fault point, the equivalent impedance from the doubly fed w turbine to the fault point is Z2L, the transition impedance is Zf, the filter resistance reactance of the grid-measured converter are Rg and Lg, and the voltage at the end o DFIG machine is us.

Study on Short-Circuit Current Variation Mechanism Considering Low-Voltage Crossin Control
As shown in Figure 1, the DFIG rotor converter can complete the low-voltage c ing by controlling the rotor excitation voltage ur during the grid voltage drop.There the effect of the excitation voltage ur cannot be ignored.
The rotor voltage equation can be obtained by combining Equations ( 2) and ( eliminate the rotor current: Converted to the rotor rotating coordinate system, Formula (4) is simplified a lows: In the formula, u r r ψ r s and ir are the rotator voltage, the stator magnetic chain, an rotator current in the rotating coordinate system.
Considering that the DC attenuation component of the short-circuit current o DFIG group is closely related to the stator-side flux dΨs/dt, the relationship betwee DC attenuation component and the stator-side flux is better discussed.This article not consider the impact of the stator flux transient process with the rotor-side conv in the inertial time of IGBT.Thus, it can be assumed that Kdc immediately changes t corresponding modulation ratio after the fault occurs, and the rotor-side voltage is qu adjusted to the reference value ur, by the rotor converter.Formula (6) can be obtaine combining Formula (1) and Formula (5) as follows.Here E g is the equivalent potential at the grid side, Z 1L is the equivalent impedance from the system to the fault point, the equivalent impedance from the doubly fed wind turbine to the fault point is Z 2L , the transition impedance is Z f , the filter resistance and reactance of the grid-measured converter are R g and L g , and the voltage at the end of the DFIG machine is us.

Study on Short-Circuit Current Variation Mechanism Considering Low-Voltage Crossing Control
As shown in Figure 1, the DFIG rotor converter can complete the low-voltage crossing by controlling the rotor excitation voltage u r during the grid voltage drop.Therefore, the effect of the excitation voltage u r cannot be ignored.
The rotor voltage equation can be obtained by combining Equations ( 2) and ( 3) to eliminate the rotor current: Converted to the rotor rotating coordinate system, Formula (4) is simplified as follows: In the formula, u r r ψ r s and i r are the rotator voltage, the stator magnetic chain, and the rotator current in the rotating coordinate system.
Considering that the DC attenuation component of the short-circuit current of the DFIG group is closely related to the stator-side flux dΨ s /dt, the relationship between the DC attenuation component and the stator-side flux is better discussed.This article does not consider the impact of the stator flux transient process with the rotor-side converter in the inertial time of IGBT.Thus, it can be assumed that K dc immediately changes to the corresponding modulation ratio after the fault occurs, and the rotor-side voltage is quickly adjusted to the reference value u r , by the rotor converter.Formula (6) can be obtained by combining Formula (1) and Formula (5) as follows.
where, the rotor transient inductance L r is L r = L r − L 2 m L −1 s .According to Equation (6), the following can be obtained: i r (t) = (i r0 − i r_ref )e −t/τ r + i r_ref (7) Sensors 2023, 23, 8372 4 of 16 Considering the parameters of the DFIG, L s = L sσ + L m , L r = L rσ + L m and L m ≥≥ L sσ , L m ≥≥ L rσ , then the stator and rotor currents can be expressed by eliminating Formula (3): Equation ( 8) reflects the opposite approximate magnitude of stator current and rotor current during the fault, so Equation ( 7) can be rewritten as follows: where i s0 and i s_ref are stator side short-circuit currents corresponding to the low-through control strategy at the initial and steady time of grid voltage sag, respectively.Equation (9) reflects that the variation rule of short-circuit current is related to the initial value and steady-state value of stator short-circuit current.

Short-Circuit Current Calculation Method of DFIG Group at the Initial Time of Grid Voltage Sag
The stator flux expression can be obtained from Equation (3): The stator voltage equation can be obtained from Equations (2,10) as follows: In the initial stage of the fault, the stator flux dψ r /dt = 0 can be considered.Then, the initial moment of short-circuit current can be simplified into the following form (Equation ( 12)): According to Equation (12), only ψ r0 is an unknown quantity in the fault current of the doubly fed unit at the initial time.
Considering that the rotor flux of the DFIG group cannot be mutated at the fault instant, before the fault, the output active power and reactive power of the DFIG group are as follows: As the stator flux orientation control strategy is adopted on the rotor side of the DFIG group, ψ sd = ψ s = u s /jω s , ψ s is the amplitude of the stator flux.Formula (3) can be rewritten as follows: (14) Sensors 2023, 23, 8372 5 of 16 According to Equations ( 13) and (14), the relationship between rotor flux and DFIG output active power, reactive power, and pre-fault voltage at the initial time of fault can be expressed as follows: ) Assuming that the stator voltage u s0 before the fault is rated, the active and reactive power depend on the setting of the working conditions before the network failure.The rotor flux ψ r is obtained from Equation (15) and the rotor flux ψ r is brought into Equation (12) to obtain the short circuit current i ' so at the initial moment of grid voltage drop.

Short-Circuit Current Calculation of DFIG Group at Fault Steady-State Time
Due to the low-voltage crossing period, the rotor-side converter needs to emit a reactive current to support the grid voltage recovery when the grid voltage drops [27][28][29].When terminal voltage sag is detected after the fault, the rotor converter control strategy is adjusted to output reactive current in the low-voltage traverse control mode, and the rotor excitation current reference value is adjusted to output reactive power according to the wind power grid connection standard [30][31][32][33].
When the fault reaches the steady state, that is, dψ s /dt = 0, the following can be obtained from Formula (2): According to Formulas (3,16), it can be obtained that: where, i r_ref = i rd_ref + ji rq_ref , i rq_ref and i rd_ref are the rotor active and reactive current reference values, respectively.
where, X = jω s L s .According to Equation ( 18), i s_ref is determined by i rq_ref , i rd_ref , R s , X, E g , Z 1L , Z 2L , Z f , L mb , and ω s , where only i rd_ref and i rd_ref are unknown quantities.
According to China's grid-connected standard, the reference values of the field current after the fault can be determined as i d_ref and i q_ref : where i rN and i rmax are the rated current and upper current of the rotor and; K d is the reactive current gain coefficient determined using the grid connection standard.
Sensors 2023, 23, 8372 6 of 16 According to Equations ( 13) and ( 14), the rotor current at the initial time can be obtained by eliminating the stator current: Formula ( 20) is brought into Formula ( 19) to obtain the following: After the grid voltage drops, the rotor excitation current is adjusted according to Equation (21).Since the active power P 0 cannot be mutated after the fault.The reactive current emission is only related to the voltage sag degree.Considering the requirement of the grid-connected code to issue reactive current to support the voltage recovery of the grid, the rotor-side converter should first determine the reference value of reactive current and then determine the reference value of active current.
Equation (21) shows the relationship between the reference value of reactive current and the drop depth, as shown in Figure 2. The dashed line area represents the control dead zone.
Sensors 2023, 23, x FOR PEER REVIEW 7 of 17 Reactive current reference value and the relationship between the degree of voltage sags.
In Figure 2, usN and Δus are the rated terminal voltage and the amount of terminal voltage drop, respectively; irN and Δird_ref are the rated rotor excitation current and reactive current reference value change s, respectively.
Assuming Xm = jωsLm, according to Equations ( 18) and ( 21), the short-circuit steadystate current is As can be seen from Equation ( 22), the short-circuit current is_ref in the fault steady state is only related to the voltage at the initial moment, the electromotive force (Eg) at the grid side, the line impedance parameter, and power us0, P0, Q0.Combining Equations ( 9), ( 12) and ( 22), we can obtain the short-circuit current value of the doubly fed wind turbine during the whole transient process.

Short-Circuit Current Calculation Method for Doubly Fed Field Group Access
Because of the different capacities and, operating conditions and the long distance between wind farms, their transient characteristics after failure are quite different.In addition, there is a strong coupling relationship between wind farms.Therefore, the research on DFIG wind farm groups cannot be simply based on the calculation method of the shortcircuit current of a single DFIG wind turbine.It is necessary to analyze the fault interaction mechanism between each wind farm in the backcourt group.Furthermore, a short-circuit In Figure 2, u sN and ∆u s are the rated terminal voltage and the amount of terminal voltage drop, respectively; i rN and ∆i rd_ref are the rated rotor excitation current and reactive current reference value change s, respectively.
Assuming X m = jω s L m , according to Equations ( 18) and ( 21), the short-circuit steadystate current is As can be seen from Equation (22), the short-circuit current i s_ref in the fault steady state is only related to the voltage at the initial moment, the electromotive force (E g ) at the grid side, the line impedance parameter, and power u s0 , P 0 , Q 0 .Combining Equations ( 9), ( 12) and ( 22), we can obtain the short-circuit current value of the doubly fed wind turbine during the whole transient process.

Short-Circuit Current Calculation Method for Doubly Fed Field Group Access
Because of the different capacities and, operating conditions and the long distance between wind farms, their transient characteristics after failure are quite different.In addition, there is a strong coupling relationship between wind farms.Therefore, the research on DFIG wind farm groups cannot be simply based on the calculation method of the short-circuit current of a single DFIG wind turbine.It is necessary to analyze the fault interaction mechanism between each wind farm in the backcourt group.Furthermore, a short-circuit current calculation method suitable for a doubly fed field group connected to the power grid is proposed.
Without loss of generality, the actual power grid in a certain area shown in Figure 3 is taken as an example for specific analysis.Since the fans in each wind farm in this example are of the same model, their transient characteristics are essentially consistent during the fault.In this paper, each wind farm is equivalent to a doubly fed wind turbine with equal capacity instead.
grid side, the line impedance parameter, and power us0, P0, Q0.Combining Equations ( 9), ( 12) and ( 22), we can obtain the short-circuit current value of the doubly fed wind turbine during the whole transient process.

Short-Circuit Current Calculation Method for Doubly Fed Field Group Access
Because of the different capacities and, operating conditions and the long distance between wind farms, their transient characteristics after failure are quite different.In addition, there is a strong coupling relationship between wind farms.Therefore, the research on DFIG wind farm groups cannot be simply based on the calculation method of the shortcircuit current of a single DFIG wind turbine.It is necessary to analyze the fault interaction mechanism between each wind farm in the backcourt group.Furthermore, a short-circuit current calculation method suitable for a doubly fed field group connected to the power grid is proposed.
Without loss of generality, the actual power grid in a certain area shown in Figure 3 is taken as an example for specific analysis.Since the fans in each wind farm in this example are of the same model, their transient characteristics are essentially consistent during the fault.In this paper, each wind farm is equivalent to a doubly fed wind turbine with equal capacity instead.Among them, E g is the system equivalent potential; Z g and Z ln are the system equivalent impedance and the impedance of the line.I g and I DFIGi are the short-circuit currents provided by the system and each doubly fed wind farm, respectively.
The rotor flux jω s L m L −1 r ψ r0 in Formula ( 12) at the initial time is equivalent to the electric potential E DFIGi , Combined with the fault equivalent circuit, the initial short-circuit current I (0) of each doubly fed wind farm can be calculated.
where Z is the grid impedance matrix; and When the initial short-circuit current I (0) of each wind farm is injected into the grid, it will affect the voltage of each node in the grid.The short-circuit current of each wind farm at the next moment can be obtained from the attenuation law of the node voltage through Equations ( 9) and (22).At the same time, this will also lead to changes in the voltage of each node of the grid, making the wind farm group and the grid interact with each other.
In order to accurately calculate the short-circuit current of the doubly fed wind farm group after being connected to the grid, in addition to considering the coupling relationship between the field group and the grid, the step size ∆t is selected as 0.1 ms, and the shortcircuit current of each wind farm in each fault period is calculated hourly.
Through the above operation solution, together with the grid impedance matrix, the voltage U (1)  i of each grid node in the first ∆t can be obtained.Since the calculation step is short, it can be considered that the voltage in this time period is constant, and the voltage U (1) i in this period is incorporated into Equation ( 22) to obtain I (1)  s_refi .The short-circuit current I (2)  DFIGi of each doubly fed wind farm at the end of the step is further obtained according to Equation (9).
According to the short-circuit current I (1) DFIGi at the end of the first ∆t period of each wind farm, combined with the grid impedance matrix, the system node voltage U (2)  i in the second ∆t can be obtained.Formula ( 22) is incorporated to obtain the I (2)  s_refi of the step, and further, the short-circuit current I (2)  DFIGi of each doubly fed wind farm at the end of the step is obtained according to Formula (9).

Similarly, the node voltage U (k)
i in the kth calculation step ∆t can be brought into Formulas (22) and ( 9) to obtain the short-circuit current I (k) DFIGi of the doubly fed wind farm at the end of the step.andcalculate each node voltage U (k+1) i in the next ∆t time period.When the calculation time is longer than the fault recovery time, that is, when t t a , the fault analysis ends, and the short-circuit current effective value of each DFIG wind farm is output, as shown in Figure 4 for the above calculation process.

Substitute
at this time period into Formula ( 22) to obtain the step .When the power grid voltage asymmetrical drop occurs, PLL can be used to obtain the positive sequence voltage phase and amplitude, so that the active and reactive current reference values can be obtained according to Equation (21).Therefore, when an asymmetric fault occurs, the expression of the positive sequence current of the DFIG is the same as that of Formula (22), which conforms to the law of change of Formula (9).According to the system circuit structure, a negative sequence equivalent circuit can be established.Further, combining the boundary conditions of asymmetric faults, a composite sequence network circuit is established.According to the above calculation method, the asymmetric short-circuit current of the doubly fed wind farm group can be solved.

Doubly Fed Wind Power Grid Group of Failure Analysis
When the DFIG adopts the low-voltage through-through control mode, the inverter can output reactive current support grid voltage support recovery and no longer latch, providing continuous excitation.During the fault, this can be thought of as DFIG excitation to produce continuous power frequency electric potential.
According to the above analysis, at the fault steady-state moment, a steady shortcircuit current that can be made of type (18) is obtained.Let ωsLmirq_ref + jωsLmird_ref be the steady-state equivalent potential EDFIGi of the doubly fed wind turbine, then the DFIG equivalent circuit can be expressed in the form of EDFIGi and steady-state impedance ZDFIGi in series in the fault steady state.
In case of asymmetric fault, the expression of DFIG rotor current the same type (21).It can be seen that in the case of asymmetric faults, only positive sequence potentials still exist in the DFIG group.When the power grid voltage asymmetrical drop occurs, PLL can be used to obtain the positive sequence voltage phase and amplitude, so that the active and reactive current reference values can be obtained according to Equation (21).Therefore, when an asymmetric fault occurs, the expression of the positive sequence current of the DFIG is the same as that of Formula (22), which conforms to the law of change of Formula (9).According to the system circuit structure, a negative sequence equivalent circuit can be established.Further, combining the boundary conditions of asymmetric faults, a composite sequence network circuit is established.According to the above calculation method, the asymmetric short-circuit current of the doubly fed wind farm group can be solved.

Doubly Fed Wind Power Grid Group of Failure Analysis
When the DFIG adopts the low-voltage through-through control mode, the inverter can output reactive current support grid voltage support recovery and no longer latch, providing continuous excitation.During the fault, this can be thought of as DFIG excitation to produce continuous power frequency electric potential.
According to the above analysis, at the fault steady-state moment, a steady shortcircuit current that can be made of type (18) is obtained.Let ω s L m i rq_ref + jω s L m i rd_ref be the steady-state equivalent potential E DFIGi of the doubly fed wind turbine, then the DFIG equivalent circuit can be expressed in the form of E DFIGi and steady-state impedance Z DFIGi in series in the fault steady state.
In case of asymmetric fault, the expression of DFIG rotor current the same type (21).It can be seen that in the case of asymmetric faults, only positive sequence potentials still exist in the DFIG group.
According to the above analysis, when symmetrical and asymmetrical faults occur in the system, the doubly fed wind turbine only has the positive internal sequence potential of power frequency in the fault steady state.According to Equations ( 18) and ( 21), the steady-state internal potential of the doubly fed wind turbine can be obtained as follows.
According to Equation ( 24), the DFIG cannot be treated as a constant voltage source before and after the fault.It is necessary to use the characteristics of the fault steady DFIG equivalent potential to analyze the interaction mechanism between the doubly fed farm group and power grid.It is necessary to propose power grid fault analysis method for after the doubly fed field group is connected.
Suppose that point A in Figure 3 has a two-phase short-circuit of AB, The positive sequence and negative sequence network diagram of the power grid fault is shown in Figure 5. Where, E g and Z g are the equivalent potential and equivalent impedance of the system, respectively; E DFIGi, Z DFIGi , and Z ln are the equivalent potential and steadystate impedance of each doubly fed wind farm and the equivalent impedance of the line L n , respectively.
Sensors 2023, 23, x FOR PEER REVIEW 10 of 17 of power frequency in the fault steady state.According to Equations ( 18) and ( 21), the steady-state internal potential of the doubly fed wind turbine can be obtained as follows., (0 According to Equation ( 24), the DFIG cannot be treated as a constant voltage source before and after the fault.It is necessary to use the characteristics of the fault steady DFIG equivalent potential to analyze the interaction mechanism between the doubly fed farm group and power grid.It is necessary to propose power grid fault analysis method for after the doubly fed field group is connected.
Suppose that point A in Figure 3 has a two-phase short-circuit of AB, The positive sequence and negative sequence network diagram of the power grid fault is shown in Figure 5. Where, Eg and Zg are the equivalent potential and equivalent impedance of the system, respectively; EDFIGi, ZDFIGi, and Zln are the equivalent potential and steady-state impedance of each doubly fed wind farm and the equivalent impedance of the line Ln, respectively.

Zg+ ZL2+ ZL3+
ZL1+ ZDFIG1+ When the DFIG group is connected to the power grid, the system will become a form of multi-power supply.After the power grid fails, considering the mutual coupling between wind farms and the influence of the control strategy, there is a strong non-linear relationship between the terminal voltage and the equivalent internal potential.The magnitude and phase of the fault current of the wind farm group in the steady state will be changed.The coupling relationship is shown in Figure 6.When the DFIG group is connected to the power grid, the system will become a form of multi-power supply.After the power grid fails, considering the mutual coupling between wind farms and the influence of the control strategy, there is a strong non-linear relationship between the terminal voltage and the equivalent internal potential.The magnitude and phase of the fault current of the wind farm group in the steady state will be changed.The coupling relationship is shown in Figure 6.As shown in Figure 6, the equivalent internal potential EDFIGi and steady-state impedance ZDFIGi of the i-th doubly fed wind farm can be obtained from Equation (24) in the fault steady state.Further, the grid admittance matrix Y can be obtained at steady state.
Among them, where: Yij is the admittance between nodes, Yg is the system equivalent admittance, and YDFIGi represents the equivalent admittance of each DFIG.From Formula (25), the short-circuit current output of each doubly fed wind farm can be obtained as follows: where Ui is the voltage of each node; Ig and IDFIGi are the short-circuit current provided by the system and each doubly fed wind farm, respectively.From Formulas ( 18) and ( 24), in the power grid voltage drop steady-state moment in the case of short-circuit current of DFIGi, the relationship with the voltage of the machine is as follows: As shown in Figure 6, the equivalent internal potential E DFIGi and steady-state impedance Z DFIGi of the i-th doubly fed wind farm can be obtained from Equation (24) in the fault steady state.Further, the grid admittance matrix Y can be obtained at steady state.
Among them, where: Y ij is the admittance between nodes, Y g is the system equivalent admittance, and Y DFIGi represents the equivalent admittance of each DFIG.From Formula (25), the short-circuit current output of each doubly fed wind farm can be obtained as follows: where U i is the voltage of each node; I g and I DFIGi are the short-circuit current provided by the system and each doubly fed wind farm, respectively.
From Formulas ( 18) and ( 24), in the power grid voltage drop steady-state moment in the case of short-circuit current of DFIG i , the relationship with the voltage of the machine is as follows: Therefore, this paper uses an iterative correction method to solve the problem.Since the nodal admittance matrix Y is a positive definite symmetric matrix, the iterative process must converge.Therefore, the output current before the fault can be selected as the initial value, and the voltage value of each node and the short-circuit current of each wind farm can be solved.Similarly, the voltages and branch currents of various nodes in the power grid can be obtained through calculations for other types of faults.Formula ( 28) is substituted into Formula (27) for solving.Because these formulas are high-order nonlinear systems and cannot be solved directly, the iterative correction solution is carried out using a computer.Since the output current before the fault is selected as the initial value, and the voltage value of each node is solved, as well as the short-circuit current of each wind farm.Similarly, computer iteration can be used to obtain node voltage and branch current in other types of faults.

Experimental Verification of Calculation Method for Short-Circuit Current of DFIG Wind Farm Group
In order to verify the correctness of the short-circuit current calculation method of doubly fed field groups, a model of a doubly fed field group connected to the power grid was established in RTDS.
Figure 7 shows the distance point of failure, respectively.The electric nearest and farthest doubly fed DFIG1 and, DFIG4 short-circuit current test results are compared with the model calculation.Therefore, this paper uses an iterative correction method to solve the problem.Since the nodal admittance matrix Y is a positive definite symmetric matrix, the iterative process must converge.Therefore, the output current before the fault can be selected as the initial value, and the voltage value of each node and the short-circuit current of each wind farm can be solved.Similarly, the voltages and branch currents of various nodes in the power grid can be obtained through calculations for other types of faults.Formula ( 28) is substituted into Formula (27) for solving.Because these formulas are high-order nonlinear systems and cannot be solved directly, the iterative correction solution is carried out using a computer.Since the output current before the fault is selected as the initial value, and the voltage value of each node is solved, as well as the short-circuit current of each wind farm.Similarly, computer iteration can be used to obtain node voltage and branch current in other types of faults.

Experimental Verification of Calculation Method for Short-Circuit Current of DFIG Wind Farm Group
In order to verify the correctness of the short-circuit current calculation method of doubly fed field groups, a model of a doubly fed field group connected to the power grid was established in RTDS.
Figure 7 shows the distance point of failure, respectively.The electric nearest and farthest doubly fed DFIG1 and, DFIG4 short-circuit current test results are compared with the model calculation.For further analysis, Figure 9 shows the comparison results of short-circui measurement of the positive sequence loop and negative sequence loop of DF DIFG4.For further analysis, Figure 9 shows the comparison results of short-circuit current measurement of the positive sequence loop and negative sequence loop of DFIG 1 and DIFG 4 .Effective values can be extracted from the instantaneous values of the short-circui current of DFIG1 and DFIG4 measured using the full Fourier algorithm.Figure 7a,b show that the effective values of the short-circuit current of DFIG1 and DFIG4 at the 0.5 s faul instant are 2.25 p.u. and 1.76 p.u., respectively.However, the results calculated in this paper are 2.31 p.u. and 1.72 p.u., and the errors are only 2.7% and 2.4%.In the case of faul steady state, the error is only 1.7% and 1.4%.Therefore, the method proposed in this pape is effective.
In order to make the experiment more accurate, we carried out multiple sets of ex perimental tests for different fault point conditions under symmetrical working condi tions and obtained the results shown in Figure 8, whose errors are all less than 4%.

Error at initial
Error at steady Error Error Error at 20ms Error at 50ms Error Error For further analysis, Figure 9 shows the comparison results of short-circuit curren measurement of the positive sequence loop and negative sequence loop of DFIG1 and DIFG4.It can be easily seen from Figure 9 that after the occurrence of the 0.5 S fault, the erro of the comparison results of the positive sequence short-circuit current of DFIG1 an DFIG4 was only 2% and 2.3%, while the error of negative sequence short-circuit curren was only 2.1% and 1.5%.After entering the fault stable state, the error of the positive se quence short-circuit current was only 1.6%,and 1.7%, and the error of the negative se quence short-circuit current was only 1.1% and, 1.6%.Therefore, good results can still b achieved when asymmetric voltage drop conditions occur against the power grid.
In addition, we carried out multiple sets of experimental tests on different fault point for asymmetric working conditions, and obtained the results shown in Figure 10.The er rors were all less than 4%, which verifies that the calculation method in this paper also ha considerable accuracy for asymmetric working conditions.

Error Error
Error Error

Error steady at initial
Error at steady Figure 10.Calculation error between calculation results and experimental results for two-phas short circuit.

Experimental Verification of Fault Analysis Method for Doubly Fed Wind Farm Group
Figure 3 is taken as an example to verify the short-circuit fault on the connected gri side of the doubly fed field group.Tables 1 and 2 show the positive and negative sequenc short-circuit current comparison of the branches corresponding to DFIG1-4 when point A It can be easily seen from Figure 9 that after the occurrence of the 0.5 S fault, the error of the comparison results of the positive sequence short-circuit current of DFIG1 and DFIG4 was only 2% and 2.3%, while the error of negative sequence short-circuit current was only 2.1% and 1.5%.After entering the fault stable state, the error of the positive sequence short-circuit current was only 1.6%,and 1.7%, and the error of the negative sequence shortcircuit current was only 1.1% and, 1.6%.Therefore, good results can still be achieved when asymmetric voltage drop conditions occur against the power grid.
In addition, we carried out multiple sets of experimental tests on different fault points for asymmetric working conditions, and obtained the results shown in Figure 10.The errors were all less than 4%, which verifies that the calculation method in this paper also has considerable accuracy for asymmetric working conditions.It can be easily seen from Figure 9 that after the occurrence of the 0.5 S fault, the error of the comparison results of the positive sequence short-circuit current of DFIG1 and DFIG4 was only 2% and 2.3%, while the error of negative sequence short-circuit current was only 2.1% and 1.5%.After entering the fault stable state, the error of the positive sequence short-circuit current was only 1.6%,and 1.7%, and the error of the negative sequence short-circuit current was only 1.1% and, 1.6%.Therefore, good results can still be achieved when asymmetric voltage drop conditions occur against the power grid.
In addition, we carried out multiple sets of experimental tests on different fault points for asymmetric working conditions, and obtained the results shown in Figure 10.The errors were all less than 4%, which verifies that the calculation method in this paper also has considerable accuracy for asymmetric working conditions.

Experimental Verification of Fault Analysis Method for Doubly Fed Wind Farm Group
Figure 3 is taken as an example to verify the short-circuit fault on the connected grid side of the doubly fed field group.Tables 1 and 2 show the positive and negative sequence short-circuit current comparison of the branches corresponding to DFIG1-4 when point A is in symmetric drop and asymmetric drop, respectively.It can be seen that the shortcircuit current amplitude error is only 2.5% and the phase angle error is only 10.43% under

Experimental Verification of Fault Analysis Method for Doubly Fed Wind Farm Group
Figure 3 is taken as an example to verify the short-circuit fault on the connected grid side of the doubly fed field group.Tables 1 and 2 show the positive and negative sequence short-circuit current comparison of the branches corresponding to DFIG1-4 when point A is in symmetric drop and asymmetric drop, respectively.It can be seen that the shortcircuit current amplitude error is only 2.5% and the phase angle error is only 10.43% under the symmetrical drop condition.This is much smaller than the calculation error of the traditional calculation method using the equivalent synchronous generator.

Conclusions
Because the traditional method is simple and equivalent to the synchronous generator, the fault analysis method is not applicable to the wind power generation mode of cluster development and centralized grid connection in our country, the influence of the LVRT control strategy on short-circuit current calculation is comprehensively considered in this paper, and the RTDS model is established.The coupling relationship between the wind farms and the voltage of the grid node is further given, finally, the grid fault analysis method suitable for doubly fed wind farm group access is proposed.
(1) The control strategy adopted during the fault period will have a great influence on the short-circuit current characteristics of doubly fed wind turbines.China's new wind power grid-connected standard puts forward requirements for the output reactive support current of doubly fed wind turbines during the fault period, and it is necessary to calculate the short-circuit current of doubly fed wind turbines and its influence.(2) There is a strong coupling relationship between the short-circuit current of each wind farm in the wind farm group and the voltage of the grid node.The coupling relationship must be considered to obtain a more accurate output short-circuit current of each wind farm.

Figure 1 .
Figure 1.DFIG equivalent circuit under the condition of power grid voltage sags.

Figure 1 .
Figure 1.DFIG equivalent circuit under the condition of power grid voltage sags.

Figure 2 .
Figure 2. Reactive current reference value and the relationship between the degree of voltage sags.

Figure 3 .
Figure 3. Diagram of power grid with doubly fed wind farm groups.Figure 3. Diagram of power grid with doubly fed wind farm groups.

Figure 3 .
Figure 3. Diagram of power grid with doubly fed wind farm groups.Figure 3. Diagram of power grid with doubly fed wind farm groups.

Figure 4 .
Figure 4. Flow chart of doubly fed wind farm group short-circuit-current calculation method.

Figure 4 .
Figure 4. Flow chart of doubly fed wind farm group short-circuit-current calculation method.

Figure 5 .
Figure 5. Composite order network diagram of power grid with DFIG field group access.(a) Positive sequence network after doubly fed wind farm group is connected to power grid; (b) negative sequence network after doubly fed wind farm group is connected to power grid.

Figure 5 .
Figure 5. Composite order network diagram of power grid with DFIG field group access.(a) Positive sequence network after doubly fed wind farm group is connected to power grid; (b) negative sequence network after doubly fed wind farm group is connected to power grid.

Figure 6 .
Figure 6.Coupling relation between doubly fed wind farm group and network.

Figure 6 .
Figure 6.Coupling relation between doubly fed wind farm group and network.

Figure 7 .
Figure 7.A point voltage three-phase short-circuit fault comparison of short-circuit current calculation and measured value of each wind farm.(a) Comparison of DFIG 1 short-circuit current calculated results and experimental results.(b) Comparison between calculated and measured DFIG 4 shortcircuit current of doubly fed wind farm.

Figure 8 .
Figure 8. Error between the calculated results of the short-circuit current model and the m results in three-phase short-circuit.

Figure 8 .
Figure 8. Error between the calculated results of the short-circuit current model and the measured results in three-phase short-circuit.
tion and measured value of each wind farm.(a) Comparison of DFIG1 short-circuit current calcu lated results and experimental results.(b) Comparison between calculated and measured DFIG short-circuit current of doubly fed wind farm.

Figure 8 .
Figure 8. Error between the calculated results of the short-circuit current model and the measured results in three-phase short-circuit.

Figure 9 .
Figure 9.The measured short-circuit current is compared with the calculated results when the two phase short circuit of the grid voltage occurs at terminal A. (a) Comparison of DFIG1 experimenta and calculated results for positive and negative sequence short-circuit current; (b) comparison o DFIG4 experimental and calculated results for positive and negative sequence short-circuit curren

Figure 9 .
Figure 9.The measured short-circuit current is compared with the calculated results when the twophase short circuit of the grid voltage occurs at terminal A. (a) Comparison of DFIG 1 experimental and calculated results for positive and negative sequence short-circuit current; (b) comparison of DFIG 4 experimental and calculated results for positive and negative sequence short-circuit current.

Figure 9 .
Figure 9.The measured short-circuit current is compared with the calculated results when the twophase short circuit of the grid voltage occurs at terminal A. (a) Comparison of DFIG1 experimental and calculated results for positive and negative sequence short-circuit current; (b) comparison of DFIG4 experimental and calculated results for positive and negative sequence short-circuit current.

Figure 10 .
Figure 10.Calculation error between calculation results and experimental results for two-phase short circuit.

Figure 10 .
Figure 10.Calculation error between calculation results and experimental results for two-phase short circuit.

Table 1 .
The calculated values of the model were compared with the experimental test values in the case of a three-phase short circuit.

Table 2 .
The model calculated values were compared with the experimental test values when the two phases of AB were short-circuited.