Simulation evaluation model of high voltage ride through capability, simulation evaluation method based on the same and storage medium

Abstract

A simulation evaluation model of a high voltage ride through capability includes a wind turbine system aerodynamic model, a torque control model, a converter model, and a high voltage fault generating device model connected in sequence; the wind turbine system aerodynamic model is configured to calculate an airflow input power; the torque control model is configured to calculate a rotor electromagnetic torque according to the airflow input power; the high voltage fault generating device model is configured to simulate a high voltage fault and output a predetermined voltage on a low voltage side of a transformer; and the converter model is configured to calculate a stator reactive current, an active power and a reactive power of the wind turbine system during the high voltage fault according to the airflow input power, the rotor electromagnetic torque and the predetermined voltage on the low voltage side of the transformer.

Claims

1. A computer-implemented simulation evaluation model of a high voltage ride through capability, comprising: a wind turbine system aerodynamic model, a torque control model, a converter model, and a high voltage fault generating device model connected in sequence; wherein the wind turbine system aerodynamic model is configured to calculate an airflow input power according to an obtained first input initialization evaluation parameter; the torque control model is configured to calculate a rotor electromagnetic torque according to the airflow input power; the high voltage fault generating device model is configured to simulate a high voltage fault according to an obtained second input initialization evaluation parameter, and output a predetermined voltage on a low voltage side of a transformer; and the converter model is configured to calculate a stator reactive current, an active power and a reactive power of a wind turbine system during the simulated high voltage fault according to the calculated airflow input power, the calculated rotor electromagnetic torque and the outputted predetermined voltage on the low voltage side of the transformer, the converter model is further configured to obtain the stator reactive current by inputting a −q-axis stator flux linkage, a predetermined stator self-inductance, a predetermined stator mutual inductance, and a predetermined q-axis rotor current to a predetermined stator reactive current calculation formula for calculation; the converter model is further configured to obtain the active power by inputting a stator rotational velocity, the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer into a predetermined active power calculation formula for calculation; and the converter model is further configured to obtain the reactive power by inputting a d-axis rotor current, a d-axis stator voltage, a stator mutual inductance, a q-axis rotor current, the stator rotational velocity, the airflow input power, and the rotor electromagnetic torque to a predetermined reactive power calculation formula for calculation; and wherein the simulation evaluation model determines whether a wind turbine system has high voltage ride through capability and generates an output that indicates whether the wind turbine system has high voltage ride through capability.

2. The computer-implemented simulation evaluation model of the high voltage ride through capability of claim 1, wherein, the wind turbine system aerodynamic model is further configured to obtain a wind turbine system blade length, a wind velocity, an air density, a wind turbine system blade angular velocity, and a wind energy availability; wherein, the obtained first input initialization evaluation parameter comprises: the wind turbine system blade length, the wind velocity, the air density, the wind turbine system blade angular velocity, and the wind energy availability; a tip speed ratio is obtained according to a ratio of a blade tip linear velocity to the wind velocity; and the airflow input power is obtained by inputting the wind turbine system blade length, the wind velocity, the air density, the wind energy availability, a blade pitch angle, and the tip speed ratio to a predetermined airflow input power calculation formula for calculation.

3. The computer implemented simulation evaluation model of the high voltage ride through capability of claim 2, wherein the torque control model is further configured to input a rotor mechanical inertia, a rotor motion damping coefficient, and a rotor mechanical angular velocity into a rotor electromagnetic torque calculation formula preset in the torque control model for calculation.

4. The computer implemented simulation evaluation model of the high voltage ride through capability of claim 2, wherein the high voltage fault generating device model is further configured to obtain a series reactor impedance value, a grid impedance, and a wind turbine system capacity, wherein, the obtained second input initialization evaluation parameter comprises: the series reactor impedance value, the grid impedance, and the wind turbine system capacity; and the predetermined voltage on the low voltage side of the transformer is obtained by inputting the series reactor impedance value, the grid impedance and the wind turbine system capacity to a predetermined calculation formula of a predetermined voltage on the low voltage side of the transformer configured to simulate the high voltage fault.

5. The computer-implemented simulation evaluation model of the high voltage ride through capability of claim 1, wherein, the airflow input power P.sub.m is calculated according to the following formula: $P_m=\frac{1}{2}\pi r_b^2{v}_w^3\rho c_p\left(\lambda ,\beta \right),$ wherein P.sub.m is the airflow input power, r is the ratio of circumference to diameter, r.sub.b is the wind turbine system blade length, v.sub.w is the wind velocity, ρ is the air density, λ is the tip speed ratio, that is, the ratio of the blade tip linear velocity to the wind velocity, $\lambda =\frac{w_b{r}_b}{v_w},$ w.sub.b is the wind turbine system blade angular velocity, β is the blade pitch angle, and c.sub.p is the wind energy availability.

6. The computer implemented simulation evaluation model of the high voltage ride through capability of claim 5, wherein the torque control model is further configured to input a rotor mechanical inertia, a rotor motion damping coefficient, and a rotor mechanical angular velocity into a rotor electromagnetic torque calculation formula preset in the torque control model for calculation.

7. The computer implemented simulation evaluation model of the high voltage ride through capability of claim 5, wherein the high voltage fault generating device model is further configured to obtain a series reactor impedance value, a grid impedance, and a wind turbine system capacity, wherein, the obtained second input initialization evaluation parameter comprises: the series reactor impedance value, the grid impedance, and the wind turbine system capacity; and the predetermined voltage on the low voltage side of the transformer is obtained by inputting the series reactor impedance value, the grid impedance and the wind turbine system capacity to a predetermined calculation formula of a predetermined voltage on the low voltage side of the transformer configured to simulate the high voltage fault.

8. The computer implemented simulation evaluation model of the high voltage ride through capability of claim 5, wherein the high voltage fault generating device model is further configured to obtain a series reactor impedance value, a grid impedance, and a wind turbine system capacity, wherein, the obtained second input initialization evaluation parameter comprises: the series reactor impedance value, the grid impedance, and the wind turbine system capacity; and the predetermined voltage on the low voltage side of the transformer is obtained by inputting the series reactor impedance value, the grid impedance and the wind turbine system capacity to a predetermined calculation formula of a predetermined voltage on the low voltage side of the transformer configured to simulate the high voltage fault.

9. The computer-implemented simulation evaluation model of the high voltage ride through capability of claim 1, wherein the torque control model is further configured to input a rotor mechanical inertia, a rotor motion damping coefficient, and a rotor mechanical angular velocity into a rotor electromagnetic torque calculation formula preset in the torque control model for calculation.

10. The computer-implemented simulation evaluation model of the high voltage ride through capability of claim 9, wherein, the rotor electromagnetic torque T.sub.e is calculated according to the following formula: $T_e=\frac{P_m}{w_{\mathrm{rm}}}-\left(J_r\frac{{\mathrm{dw}}_{\mathrm{rm}}}{\mathrm{dt}}+k_r{w}_{\mathrm{rm}}\right),$ wherein, T.sub.e is the rotor electromagnetic torque, P.sub.m is the airflow input power, J.sub.r is the rotor mechanical inertia, k.sub.r is the rotor motion damping coefficient, and w.sub.rm is the rotor mechanical angular velocity.

11. The computer-implemented simulation evaluation model of the high voltage ride through capability of claim 1, wherein the high voltage fault generating device model is further configured to obtain a series reactor impedance value, a grid impedance, and a wind turbine system capacity, wherein, the obtained second input initialization evaluation parameter comprises: the series reactor impedance value, the grid impedance, and the wind turbine system capacity; and the predetermined voltage on the low voltage side of the transformer is obtained by inputting the series reactor impedance value, the grid impedance and the wind turbine system capacity to a predetermined calculation formula of a predetermined voltage on the low voltage side of the transformer configured to simulate the high voltage fault.

12. The computer-implemented simulation evaluation model of the high voltage ride through capability of claim 11, wherein, the predetermined voltage u.sub.s on the low voltage side of the transformer is calculated as follows: $X_1=\frac{u_s^2}{3*S_{\mathrm{WT}}}-X_{\mathrm{grid}},$ wherein, u.sub.s is the predetermined voltage on the low voltage side of the transformer, i.e., a stator voltage of a converter, X.sub.1 is the series reactor impedance value, X.sub.grid is the grid impedance, and S.sub.WT is the wind turbine system capacity.

13. The computer-implemented simulation evaluation model of the high voltage ride through capability of claim 1, wherein, a calculation formula of the stator reactive current i.sub.qs is as follows: $i_{\mathrm{qs}}=\frac{{\lambda }_{\mathrm{qs}}-L_M{i}_{\mathrm{qr}}}{L_S},$ wherein, i.sub.qs is a q-axis stator current, i.e. the stator reactive current, λ.sub.qs is the q-axis stator flux linkage, L.sub.S is the predetermined stator self inductance, L.sub.M is the stator mutual inductance, and i.sub.qr is the q-axis rotor current; calculation formulas of the active power P.sub.S and the reactive power Q.sub.S are respectively as follows: $\{\begin{array}{c}{P}_S=\frac{T_e{\omega }_s}{P_m}\\ Q_S=-\frac{T_e}{P_m{i}_{\mathrm{dr}}}\left({\omega }_s{i}_{\mathrm{qr}}+\frac{v_{\mathrm{ds}}}{L_M}\right)\end{array},$ wherein P.sub.S is the active power, Q.sub.S is the reactive power, P.sub.m is the airflow input power, T.sub.e is the rotor electromagnetic torque, and w.sub.s is the stator rotational velocity, i.sub.dr and i.sub.qr are the d-axis rotor current and the q-axis rotor current, respectively, vas is the d-axis stator voltage, and L.sub.M is the stator mutual inductance.

14. A simulation evaluation method of the high voltage ride through capability based on the computer implemented simulation evaluation model of the high voltage ride through capability of claim 1, comprising the following steps: calculating the airflow input power by the wind turbine system aerodynamic model according to the obtained first input initialization evaluation parameter; calculating the rotor electromagnetic torque by the torque control model according to the airflow input power obtained; simulating the high voltage fault and outputting the predetermined voltage on the low voltage side of the transformer by the high voltage fault generating device model according to the obtained second input initialization evaluation parameter; calculating the stator reactive current, the active power and the reactive power of the wind turbine system during the high voltage fault by the converter model according to the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer; and determining whether the wind turbine system has the high voltage ride through capability according to the predetermined voltage on the low voltage side of the transformer, the stator reactive current, the active power, and the reactive power, and generating an output that indicates whether the wind turbine system has high voltage ride through capability.

15. The simulation evaluation method of the high voltage ride through capability based on the simulation evaluation model of the high voltage ride through capability of claim 14, wherein, the step of calculating the airflow input power by the wind turbine system aerodynamic model according to the obtained first input initialization evaluation parameter, comprises: obtaining a wind turbine system blade length, a wind velocity, a air density, a wind turbine system blade angular velocity, and a wind energy availability by the wind turbine system aerodynamic model, wherein, the initialization evaluation parameter comprises: the wind turbine system blade length, the wind velocity, the air density, the wind turbine system blade angular velocity, and the wind energy availability; obtaining a tip speed ratio according to a ratio of a blade tip linear velocity to the wind velocity; and obtaining the airflow input power by inputting the wind turbine system blade length, the wind velocity, the air density, the wind energy availability, a blade pitch angle, and the tip speed ratio to a predetermined airflow input power calculation formula for calculation by the wind turbine system aerodynamic model.

16. The simulation evaluation method of the high voltage ride through capability based on the computer implemented simulation evaluation model of the high voltage ride through capability of claim 14, wherein, the step of calculating the rotor electromagnetic torque by the torque control model according to the airflow input power obtained, comprises: inputting the torque control model inputting a rotor mechanical inertia, a rotor motion damping coefficient, and a rotor mechanical angular velocity to a rotor electromagnetic torque calculation formula preset in the torque control model for calculation.

17. The simulation evaluation method of the high voltage ride through capability based on the computer implemented simulation evaluation model of the high voltage ride through capability of claim 14, wherein, the step of simulating the high voltage fault and outputting the predetermined voltage on the low voltage side of the transformer by the high voltage fault generating device model according to the obtained second input initialization evaluation parameter, comprises: obtaining a series reactor impedance value, a grid impedance, and a wind turbine system capacity by the high voltage fault generating device model, wherein, the initialization evaluation parameter comprises: the series reactor impedance value, the grid impedance, and the wind turbine system capacity; and obtaining the predetermined voltage on the low voltage side of the transformer by inputting the series reactor impedance value, the grid impedance and the wind turbine system capacity to a predetermined calculation formula of a predetermined voltage on the low voltage side of the transformer configured to simulate the high voltage fault.

18. The simulation evaluation method of the high voltage ride through capability based on the computer implemented simulation evaluation model of the high voltage ride through capability of claim 14, wherein, the step of calculating the stator reactive current, the active power and the reactive power of the wind turbine system during the high voltage fault by the converter model according to the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer, comprises: obtaining the stator reactive current by inputting a q-axis stator flux linkage, a predetermined stator self inductance, a predetermined stator mutual inductance, and a predetermined q-axis rotor current to a predetermined stator reactive current calculation formula for calculation in the converter model; obtaining the active power by inputting a stator rotational velocity, the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer into a predetermined active power calculation formula for calculation by the converter model; and obtaining the reactive power by inputting a d-axis rotor current, a d-axis stator voltage, a stator mutual inductance, the q-axis rotor current, the stator rotational velocity, the airflow input power, and the rotor electromagnetic torque to a predetermined reactive power calculation formula for calculation by the converter model.

19. A readable storage medium, comprising executable instructions stored in the readable storage medium, wherein, when the instructions are executed by one or more processors, the simulation evaluation method of the high voltage ride through capability based on the computer implemented simulation evaluation model of the high voltage ride through capability of claim 14 is realized.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The above description is merely the summary of the technical solutions in the present disclosure, and in order to specify the technical means of the present disclosure, the present disclosure will be further described in detail below with reference to the accompanying drawings and specific embodiments.

(2) FIG. 1 shows a curve of a high voltage fault in a test procedure (discussion draft) on wind turbine system fault ride through capability;

(3) FIG. 2 is a structural schematic diagram of a simulation evaluation model of a high voltage ride through capability in the present disclosure;

(4) FIG. 3 shows a simulation result curve of a per unit value of a stator voltage of a converter when an evaluation of a three-phase symmetrical voltage fault is performed by a simulation evaluation method of the present disclosure;

(5) FIG. 4 shows a simulation result curve of a per unit value of a stator reactive current when an evaluation of a three-phase symmetrical voltage fault is performed by a simulation evaluation method of the present disclosure;

(6) FIG. 5 shows a simulation result curve of a per unit value of a reactive power when an evaluation of a three-phase symmetrical voltage fault is performed by a simulation evaluation method of the present disclosure;

(7) FIG. 6 shows a simulation result curve of a per unit value of an active power when an evaluation of a three-phase symmetrical voltage fault is performed by a simulation evaluation method of the present disclosure;

(8) FIG. 7 shows a simulation result curve of a per unit value of a stator voltage of a converter when an evaluation of a two-phase asymmetric voltage fault is performed by a simulation evaluation method of the present disclosure;

(9) FIG. 8 shows a simulation result curve of a per unit value of a stator reactive current when an evaluation of a two-phase asymmetric voltage fault is performed by a simulation evaluation method of the present disclosure;

(10) FIG. 9 shows a simulation result curve of a per unit value of a reactive power when an evaluation of a two-phase asymmetric voltage fault is performed by a simulation evaluation method of the present disclosure; and

(11) FIG. 10 shows a simulation result curve of a per unit value of an active power when an evaluation of a two-phase asymmetric voltage fault is performed by a simulation evaluation method of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(12) Referring to FIG. 2, the present disclosure provides a simulation evaluation model of a high voltage ride through capability, including a wind turbine system aerodynamic model, a torque control model, a converter model, and a high voltage fault generating device model. The simulation evaluation model of the high voltage ride through capability can run on a terminal device, wherein the terminal device may be a hardware device such as a personal computer, a mobile device, etc.

(13) Further, the wind turbine system aerodynamic model may be a software module configured to calculate an airflow input power according to an input initialization evaluation parameter.

(14) In this embodiment, the wind turbine system aerodynamic model can obtain the initialization evaluation parameter first, wherein, the initialization evaluation parameter includes: a wind turbine system blade length, a wind velocity, an air density, a wind turbine system blade angular velocity, and a wind energy availability. Thus, the wind turbine system aerodynamic model first obtains a tip speed ratio according to a ratio of a blade tip linear velocity to the wind velocity in the obtained wind turbine system blade length, the wind velocity, the air density, the wind turbine system blade angular velocity, and the wind energy availability. Subsequently, the wind turbine system aerodynamic model calls the predetermined airflow input power calculation formula and calculates to obtain the airflow input power according to the wind turbine system blade length, the wind velocity, the air density, the wind energy availability, the blade pitch angle, and the tip speed ratio.

(15) Optionally, the air flow input power P.sub.m is calculated according to the following formula:

(16) $\begin{array}{cc}{P}_m=\frac{1}{2}\pi r_b^2{v}_w^3\rho c_p\left(\lambda ,\beta \right)& \left(1\right)\end{array}$

(17) where, P.sub.m is the airflow input power, π is the ratio of circumference to diameter, r.sub.b is the wind turbine system blade length, v.sub.w is the wind velocity, ρ is the air density, λ is the tip speed ratio, that is, a ratio of the blade tip linear velocity to the wind velocity,

(18) $\lambda =\frac{w_b{r}_b}{v_w},$
w.sub.b is the wind turbine system blade angular velocity, β is the blade pitch angle, and c.sub.p is the wind energy availability.

(19) The empirical values of the wind energy availability c.sub.p are shown in Table 1 below:

(20) TABLE-US-00001 TABLE 1 Value table of wind energy availability cp Wind velocity (m/s) Availability 5 0.480443 6 0.480442 7 0.480443 8 0.471194 9 0.369719 10 0.269514

(21) The torque control model is configured to calculate the rotor electromagnetic torque according to the airflow input power output by the wind turbine system aerodynamic model. That is, a wind turbine system impeller, a transmission shaft and a generator rotor are simulated as an inertial body, and a transmission process of a shaft torque is simulated by one first order inertia element.

(22) In this embodiment, the torque control model obtains the rotor electromagnetic torque by the following method.

(23) The torque control model can call a predetermined rotor electromagnetic torque calculation formula, and the rotor electromagnetic torque can be obtained by inputting the obtained rotor mechanical inertia, the rotor motion damping coefficient and the rotor mechanical angular velocity to the torque control model.

(24) Optionally, the rotor electromagnetic torque T.sub.e is calculated according to the following formula in the torque control model:

(25) $\begin{array}{cc}{T}_m=\frac{P_m}{w_{\mathrm{rm}}}& \left(2\right)\\ J_r\frac{{\mathrm{dw}}_{\mathrm{rm}}}{\mathrm{dt}}+k_r{w}_{\mathrm{rm}}=T_m-T_e& \left(3\right)\end{array}$

(26) where, T.sub.e is the rotor electromagnetic torque, T.sub.m is a generator torque, J.sub.r is the rotor mechanical inertia, k.sub.r is the rotor motion damping coefficient, and w.sub.rm is the rotor mechanical angular velocity.

(27) $\begin{array}{cc}{T}_e=\frac{P_m}{w_{\mathrm{rm}}}-\left(J_r\frac{{\mathrm{dw}}_{\mathrm{rm}}}{\mathrm{dt}}+k_r{w}_{\mathrm{rm}}\right)& \left(4\right)\end{array}$

(28) The high voltage fault generating device model is configured to simulate a high voltage fault according to an input initialization evaluation parameter, and calculate a predetermined voltage on a low voltage side of a transformer.

(29) In this embodiment, the high voltage fault generating device model calculates the predetermined voltage on the low voltage side of the transformer by the following method.

(30) First, the initialization evaluation parameter further includes: a series reactor impedance value, a grid impedance, and a wind turbine system capacity. Thus, the high voltage fault generating device model can further obtain the input series reactor impedance value, the input grid impedance and the input wind turbine system capacity. Subsequently, the high voltage fault generating device model can call a transformer low-voltage-side predetermined-voltage calculation formula for calculation configured to simulate the high voltage fault, and the predetermined voltage on the low voltage side of the transformer is obtained based on the series reactor impedance value, the grid impedance and the wind turbine system capacity.

(31) Optionally, the predetermined voltage u.sub.s on the low voltage side of the transformer is calculated as follows:

(32) 0$\begin{array}{cc}{X}_1=\frac{u_s^2}{3*S_{\mathrm{WT}}}-X_{\mathrm{grid}}& \left(5\right)\end{array}$
where, u.sub.s is the predetermined voltage on the low voltage side of the transformer, i.e., a stator voltage of a converter, X.sub.1 is the series reactor impedance value, X.sub.grid is the grid impedance, and S.sub.WT is the wind turbine system capacity. Different predetermined voltages u.sub.s on the low voltage side of the transformer can be generated by adjusting X.sub.1 according to different test requirements.

(33) The converter model can respectively calculate the stator reactive current, the active power and the reactive power of the wind turbine system during the high voltage fault according to the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer.

(34) In this embodiment, the converter model calculates the stator reactive current, the active power and the reactive power of the wind turbine system by the following method.

(35) First, the converter model obtains a q-axis stator flux linkage that can be input by a user. The stator reactive current is obtained by inputting the obtained q-axis stator flux linkage, a predetermined stator self inductance, a predetermined stator mutual inductance and a predetermined q-axis rotor current to a predetermined stator reactive current calculation formula for calculation in the converter model.

(36) Second, the converter model can also obtain a stator rotational velocity, wherein the stator rotational velocity can be input by the user as well. The active power is obtained by inputting the obtained stator rotational velocity, the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer into a predetermined active power calculation formula for calculation by the converter model.

(37) Furthermore, the converter model can also obtain a d-axis rotor current, a d-axis stator voltage, and a stator mutual inductance, wherein the d-axis rotor current, the d-axis stator voltage and the stator mutual inductance can be input by the user. Subsequently, the reactive power is obtained by inputting the d-axis rotor current, the d-axis stator voltage, the stator mutual inductance, the q-axis rotor current, the stator rotational velocity, the airflow input power, and the rotor electromagnetic torque to a predetermined reactive power calculation formula for calculation by the converter model.

(38) Optionally, the converter model is configured to output the stator reactive current i.sub.qs, the active power P.sub.s and the reactive power Q.sub.s of the wind turbine system during the high voltage fault according to the airflow input power P.sub.m, the rotor electromagnetic torque T.sub.e and the predetermined voltage u.sub.s on the low voltage side of the transformer.

(39) The rotor electromagnetic torque T.sub.e can be expressed as follows:
T.sub.e=1.5P.sub.m(i.sub.qsλ.sub.ds−i.sub.dsλ.sub.qs)  (6)
where, i.sub.ds and i.sub.qs are the d-axis stator current and the q-axis stator current respectively, i.e. the stator reactive current and the stator active current. λ.sub.ds and λ.sub.qs are the d-axis stator flux linkage and the q-axis stator flux linkage respectively. The calculation formulas of λ.sub.ds and λ.sub.qs are as follows:

(40) $\begin{array}{cc}\{\begin{array}{c}{\lambda }_{\mathrm{ds}}=L_S{i}_{\mathrm{ds}}+L_M{i}_{\mathrm{dr}}\\ {\lambda }_{\mathrm{qs}}=L_S{i}_{\mathrm{qs}}+L_M{i}_{\mathrm{qr}}\end{array}& \left(7\right)\end{array}$

(41) where, L.sub.S and L.sub.M are the stator self inductance and stator mutual inductance respectively; i.sub.dr and i.sub.qr are the d-axis rotor current and the q-axis rotor current respectively, i.e. the rotor reactive current and the rotor active current. The calculation formulas for i.sub.ds and i.sub.qs are as follows:

(42) $\begin{array}{cc}\{\begin{array}{c}{i}_{\mathrm{ds}}=\frac{{\lambda }_{\mathrm{ds}}-L_M{i}_{\mathrm{dr}}}{L_S}\\ i_{\mathrm{qs}}=\frac{{\lambda }_{\mathrm{qs}}-L_M{i}_{\mathrm{qr}}}{L_S}\end{array}& \left(8\right)\end{array}$

(43) The following formula is obtained by putting the formula (8) is into the formula (6):

(44) $\begin{array}{cc}{T}_e=\frac{1.5P_m{L}_M}{L_S}\left(-i_{\mathrm{qr}}{\lambda }_{\mathrm{ds}}+i_{\mathrm{dr}}{\lambda }_{\mathrm{qs}}\right)& \left(9\right)\end{array}$

(45) As can be seen from the above formula, the rotor electromagnetic torque T.sub.e can be expressed as a function of the rotor current and the stator flux linkage.

(46) Further, in a steady state of a generator, a stator voltage vector {right arrow over (u.sub.S)} can be expressed as follows:
{right arrow over (u.sub.S)}=R.sub.S{right arrow over (I.sub.S)}+jω.sub.S{right arrow over (λ)}.sub.Sλ.sub.qs  (10)
where, {right arrow over (u.sub.S)}⋅ is the stator voltage vector, Rs is a stator resistance, {right arrow over (I.sub.S)}⋅ is a stator current vector, j is a vector imaginary part, ω.sub.s is the stator rotational velocity, and {right arrow over (λ)}.sub.s is a stator flux linkage vector.

(47) After conversion by d-axis coordinate and q-axis coordinates, it can be expressed as follows:
(v.sub.ds+jv.sub.qs)=R.sub.S(i.sub.ds+ji.sub.qs)+jω.sub.s(λ.sub.ds+j)  (11)

(48) where, v.sub.dc is the d-axis stator voltage, and v.sub.qc is the q-axis stator voltage.

(49) As can be seen from the above formula, the d-axis stator flux linkage and the q-axis stator flux linkage can be respectively expressed as follows:

(50) $\begin{array}{cc}\{\begin{array}{c}{\lambda }_{\mathrm{ds}}=\frac{v_{\mathrm{qs}}-R_S{i}_{\mathrm{qs}}}{{\omega }_s}\\ {\lambda }_{\mathrm{qs}}=-\frac{v_{\mathrm{ds}}-R_S{i}_{\mathrm{ds}}}{{\omega }_s}\end{array}& \left(12\right)\end{array}$

(51) The following formula is obtained by putting the formula (12) into the formula (9):

(52) $\begin{array}{cc}{T}_e=\frac{3P_m{L}_M}{2{\omega }_s{L}_S}\left(-i_{\mathrm{qr}}{v}_{\mathrm{qs}}+R_S{i}_{\mathrm{qs}}{i}_{\mathrm{qr}}+R_S{i}_{\mathrm{ds}}{i}_{\mathrm{dr}}-i_{\mathrm{dr}}{v}_{\mathrm{ds}}\right)& \left(13\right)\end{array}$

(53) In an orientation control of the stator voltage, let v.sub.qs=0, the formula can be simplified, and in a doubly fed induction generator, R.sub.s is small, and the equation is finally simplified as follows:

(54) $\begin{array}{cc}{T}_e=-\frac{3P_m{L}_M}{2{\omega }_s{L}_S}{i}_{\mathrm{dr}}{v}_{\mathrm{ds}}& \left(14\right)\end{array}$

(55) As can be seen from the above formula, the rotor electromagnetic torque T.sub.e can be expressed as a function of the d-axis rotor current i.sub.dr and the stator voltage v.sub.ds.

(56) The formulas for calculating the active power and the reactive power are as follows:

(57) $\begin{array}{cc}\{\begin{array}{c}{P}_S=1.5\left(v_{\mathrm{ds}}{i}_{\mathrm{ds}}+v_{\mathrm{qs}}{i}_{\mathrm{qs}}\right)\\ Q_S=1.5\left(v_{\mathrm{qs}}{i}_{\mathrm{ds}}-v_{\mathrm{ds}}{i}_{\mathrm{qs}}\right)\end{array}& \left(15\right)\end{array}$

(58) where, P.sub.S is the active power, Q.sub.S is the reactive power, and the stator voltage orientation control is performed, let v.sub.qs=0, the above formula can be simplified as follows:

(59) $\begin{array}{cc}\{\begin{array}{c}{P}_S=1.5v_{\mathrm{ds}}{i}_{\mathrm{ds}}\\ Q_S=-1.5v_{\mathrm{ds}}{i}_{\mathrm{qs}}\end{array}& \left(16\right)\end{array}$

(60) the following formula is obtained by putting the formula (8) into the formula (16):

(61) $\begin{array}{cc}\{\begin{array}{c}{P}_S=1.5v_{\mathrm{ds}}\left(\frac{{\lambda }_{\mathrm{ds}}-L_M{i}_{\mathrm{dr}}}{L_S}\right)\\ Q_S=-1.5v_{\mathrm{ds}}\left(\frac{{\lambda }_{\mathrm{qs}}-L_M{i}_{\mathrm{qr}}}{L_S}\right)\end{array}& \left(17\right)\end{array}$

(62) the following formula is obtained according to the above formula:

(63) 0$\begin{array}{cc}\{\begin{array}{c}{i}_{\mathrm{dr}}=-\frac{2L_S}{3v_{\mathrm{ds}}{L}_M}{P}_S+\frac{1}{L_M}{\lambda }_{\mathrm{ds}}\\ i_{\mathrm{qr}}=\frac{2L_S}{3v_{\mathrm{ds}}{L}_M}{Q}_S+\frac{1}{L_M}{\lambda }_{\mathrm{qs}}\end{array}& \left(18\right)\end{array}$

(64) the following formula is obtained by putting the formula (12) into the formula (18):

(65) $\begin{array}{cc}\{\begin{array}{c}{i}_{\mathrm{dr}}=-\frac{2L_S}{3v_{\mathrm{ds}}{L}_M}{P}_S+\frac{v_{\mathrm{qs}}-R_S{i}_{\mathrm{qs}}}{{\omega }_s{L}_M}\\ i_{\mathrm{qr}}=\frac{2L_S}{3v_{\mathrm{ds}}{L}_M}{Q}_S-\frac{v_{\mathrm{ds}}-R_S{i}_{\mathrm{ds}}}{{\omega }_s{L}_M}\end{array}& \left(19\right)\end{array}$

(66) when v.sub.qs=0, and the stator resistance R.sub.S is ignored, the following formula can be obtained:

(67) $\begin{array}{cc}\{\begin{array}{c}{i}_{\mathrm{dr}}=-\frac{2L_S}{3v_{\mathrm{ds}}{L}_M}{P}_S\\ i_{\mathrm{qr}}=\frac{2L_S}{3v_{\mathrm{ds}}{L}_M}{Q}_S-\frac{v_{\mathrm{ds}}}{{\omega }_s{L}_M}\end{array}& \left(20\right)\end{array}$

(68) as can be concluded that after the stator voltage is determined, the active power P.sub.S and the reactive power Q.sub.S can be calculated according to the d-axis rotor current and the q-axis rotor current.

(69) The following formula is obtained by putting the formula (14) into the formula (20):

(70) $\begin{array}{cc}\{\begin{array}{c}{P}_S=\frac{T_e{\omega }_s}{P_m}\\ Q_S=-\frac{T_e}{P_m{i}_{\mathrm{dr}}}\left({\omega }_s{i}_{\mathrm{qr}}+\frac{v_{\mathrm{ds}}}{L_M}\right)\end{array}& \left(21\right)\end{array}$

(71) The present disclosure provides a simulation evaluation method of a high voltage ride through capability based on the above simulation evaluation model of the high voltage ride through capability, including the following steps:

(72) S1, the wind turbine system aerodynamic model obtains the airflow input power according to the input initialization evaluation parameter;

(73) S2, the rotor electromagnetic torque is obtained by the torque control model according to the airflow input power obtained;

(74) S3, the high voltage fault generating device model simulates the high voltage fault and outputs the predetermined voltage on the low voltage side of the transformer according to the input initialization evaluation parameter;

(75) S4, the converter model calculates the stator reactive current, the active power and the reactive power of the wind turbine system during the high voltage fault according to the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer; and

(76) S5, whether the wind turbine system has the high voltage ride through capability is determined according to the predetermined voltage on the low voltage side of the transformer, the stator reactive current, the active power, and the reactive power.

(77) Further, as an optional method, the step S1 specifically includes:

(78) S11, the wind turbine system aerodynamic model obtains an input wind turbine system blade length, an input wind velocity, an input air density, an input wind turbine system blade angular velocity, and an input wind energy availability, wherein, the initialization evaluation parameter includes: the wind turbine system blade length, the wind velocity, the air density, the wind turbine system blade angular velocity, and the wind energy availability;

(79) S12, a tip speed ratio is obtained according to a ratio of a blade tip linear velocity to the wind velocity; and

(80) S13, the wind turbine system aerodynamic mode inputs the wind turbine system blade length, the wind velocity, the air density, the wind energy availability, the blade pitch angle, and the tip speed ratio to a predetermined airflow input power calculation formula for calculation to obtain the airflow input power.

(81) Further, as an optional method, the step S2 specifically includes: the obtained rotor mechanical inertia, the rotor motion damping coefficient, and the rotor mechanical angular velocity are input by the torque control module to a rotor electromagnetic torque calculation formula preset in the torque control model for calculation to obtain the rotor electromagnetic torque.

(82) Further, as an optional method, the step S3 specifically includes:

(83) S31, the high voltage fault generating device model obtains an input series reactor impedance value, an input grid impedance, and an input wind turbine system capacity, wherein, the initialization evaluation parameter includes: the series reactor impedance value, the grid impedance, and the wind turbine system capacity; and

(84) S32, the predetermined voltage on the low voltage side of the transformer is obtained by inputting the series reactor impedance value, the grid impedance and the wind turbine system capacity into a predetermined calculation formula of a predetermined voltage on a low voltage side of a transformer configured to simulate the high voltage fault for calculation by the high voltage fault generating device model.

(85) Further, as an optional method, the step S4 specifically includes:

(86) S41, the stator reactive current is obtained by inputting an obtained q-axis stator flux linkage, a predetermined stator self inductance, a predetermined stator mutual inductance and a predetermined q-axis rotor current to a predetermined stator reactive current calculation formula for calculation in the converter model.

(87) S42, the active power is obtained by inputting the obtained stator rotational velocity, the airflow input power, the rotor electromagnetic torque, and the predetermined voltage on the low voltage side of the transformer into a predetermined active power calculation formula for calculation by the converter model; and

(88) S43, the reactive power is obtained by inputting the obtained d-axis rotor current, d-axis stator voltage, the stator mutual inductance, the q-axis rotor current, the stator rotational velocity, the airflow input power, and the rotor electromagnetic torque to a predetermined reactive power calculation formula for calculation by the converter model.

(89) It should be noted that, as those skilled in the art can clearly understand, for the convenience and brevity of the description, the specific implementation process of the above described method can refer to the corresponding processes in the above embodiment about the systems, devices, and units, which is not described in detail herein again.

(90) Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the embodiments of the present application may be implemented with entire hardware, entire software or a combination of software and hardware. Moreover, the embodiments of the present application may be a computer program product form implemented in a computer storage medium (including but not limited to disk storage, compact disc read-only memory (CD-ROM), optical storage, and others) including one or more computer program codes.

(91) The simulation evaluation method of the high voltage ride through capability in this embodiment was performed by a wind turbine system with a rated power of 2000 kW and an impeller diameter of 96 m. The results are shown as follows:

(92) 1. In the simulation evaluation of the three-phase symmetrical voltage fault, the stator voltage of the converter in the wind turbine system was increased to 1.3 PU of the rated voltage, and the duration of the fault was 500 ms:

(93) FIG. 3 shows a simulation result curve of a per unit value of the stator voltage of the converter during the high voltage ride through. Since the wind turbine system output an inductive reactive power, the high voltage occurred at the high voltage side of the grid connection point of the wind power plant, and an action of the line impedance voltage division, the stator voltage of the converter actually merely rose to 1.22 PU of the rated voltage.

(94) FIG. 4 shows a simulation result curve of a per unit value of the stator reactive current during the high voltage ride through. The wind turbine system output an inductive reactive current. After an adjustment for 30 ms, the stator reactive current was stable at 65% of the rated current.

(95) FIG. 5 shows a simulation result curve of a per unit value of the reactive power during the high voltage ride through. The wind turbine system output the inductive reactive power. After an adjustment for 30 ms, the reactive power was stable at 75% of the rated current.

(96) FIG. 6 shows a simulation result curve of a per unit value of the active power during the high voltage ride through. The active power of the wind turbine system rapidly increased to 1.2 times of the rated power when the voltage suddenly changes, and then gradually decreased. When the grid voltage returns to normal, the active power dropped to 90% of the rated power, and was then adjusted to the rated power after 3 cycles.

(97) 2. In the simulation evaluation of the two-phase asymmetric voltage fault, the stator voltage of the converter in the wind turbine system was increased to 1.3 PU of the rated voltage, and the duration of the fault was 500 ms:

(98) FIG. 7 shows a simulation result curve of a per unit value of the stator voltage of the converter during the high voltage ride through. Since the wind turbine system output the inductive reactive power, and action of the line impedance voltage division, the stator voltage of the converter actually merely rose to 1.1 PU of the rated voltage.

(99) FIG. 8 shows a simulation result curve of a per unit value of the stator reactive current during the high voltage ride through. The wind turbine system output the inductive reactive current. After an adjustment for 30 ms, the stator reactive current was stable at 35% of the rated current.

(100) FIG. 9 shows a simulation result curve of a per unit value of the reactive power during the high voltage ride through. The wind turbine system output the inductive reactive power. After an adjustment for 30 ms, the reactive power was stable at 35% of the rated power.

(101) FIG. 10 shows a simulation result curve of a per unit value of the active power during the high voltage ride through. The active power of the wind turbine system rapidly increased to 1.2 times of the rated power when the voltage suddenly changes, and then is stable at 0.4 time of the rated power. When the grid voltage returns to normal, the active power dropped to 10% of the rated power, and was then adjusted to the rated power after 3 cycles.

(102) The simulation evaluation model of the high voltage ride through capability, and the simulation evaluation method based thereon of the present disclosure can simulate a high voltage fault generated in high-level debugging process for direct current of power grid in practice, or cascading high voltage faults generated due to a plurality of reactive power generating devices in a wind power plant when failing to timely exit after the low voltage ride through, which is caused by cascading low voltage faults, and can evaluate the operating state in the fault process of the wind turbine system, determine whether the wind turbine system has the high voltage ride through capability by the simulation method, and can replace the high voltage ride through mobile vehicle test equipment on the spot, which reduces the test cost.

(103) The above descriptions are only preferred embodiments of the present disclosure, and are not intended to limit the present disclosure in any way. Simple modifications, equivalent changes or modifications by using the technical content disclosed above may be made by those skilled in the art, which all fall within the scope of the present disclosure.

INDUSTRIAL APPLICABILITY

(104) The simulation evaluation model of the high voltage ride through capability, and the simulation evaluation method based thereon of the present disclosure can simulate a high voltage fault generated in high-level debugging process for direct current of power grid in practice, or cascading high voltage faults generated due to a plurality of reactive power generating devices in a wind power plant when failing to timely exit after the low voltage ride through, which is caused by cascading low voltage faults, and can evaluate the operating state in the fault process of the wind turbine system, determine whether the wind turbine system has the high voltage ride through capability by the simulation method, and can replace the high voltage ride through mobile vehicle test equipment on the spot, which reduces the test cost.