Method and system for tracking the oscillation path of a direct drive wind farm

Abstract

A method and system for tracking the oscillation path of a direct drive wind farm involve: collecting the port current, the current and state variables of the current inner loop, the dynamic angle and state component of the phase-locked loop and the phase and amplitude of the fault point voltage of each wind turbine generator in a dominant oscillation mode period which is used as acquisition period, and then obtaining corresponding instantaneous value of each oscillation component of the direct drive wind turbine generator; obtaining energy flow of coupling energy transfer path of direct drive wind farm in each acquisition period based on the instantaneous value; tracking the key energy transfer path of the direct drive wind farm according to the energy flow, so as to realize the stability control of the direct drive wind farm.

Claims

1. A method for tracking oscillation path of a direct drive wind farm, comprising: when the direct drive wind farm failing, collecting a port current, current and state variables of a current inner loop, dynamic angle and state component of a phase-locked loop (PLL) of each direct drive wind turbine generator, and phase and amplitude of a voltage of a fault point in an acquisition period which is used as a dominant oscillation mode period, and then obtaining a corresponding instantaneous value of each oscillation component of said each direct drive wind turbine generator; obtaining an energy flow of a coupling energy transfer path of the direct drive wind farm in each acquisition period based on the corresponding instantaneous value of each oscillation component collected in said each acquisition period; tracking a key energy transfer path of the direct drive wind farm according to the energy flow of the coupling energy transfer path of the direct drive wind farm in said each acquisition period, so as to realize stability control of the direct drive wind farm.

2. The method according to claim 1, wherein the energy flow of the coupling energy transfer path between subsystems of the direct drive wind farm in said each acquisition period is obtained by the following way: establishing a dynamic mathematical model of each subsystem of said each direct drive wind turbine generator in said each acquisition period based on the corresponding instantaneous value of each oscillation component collected in said each acquisition period, wherein the subsystems include current loop d, q axis subsystems, grid-connected line d, q axis subsystems, DC voltage outer loop subsystem and phase-locked loop subsystem; obtaining a dynamic coupling energy derivative of each subsystem of said each direct drive wind turbine generator according to the dynamic mathematical model of each subsystem of said each direct drive wind turbine generator in said each acquisition period, and then obtaining a dynamic coupling energy derivative of said each direct drive wind turbine generator in said each acquisition period, a value corresponding to the dynamic coupling energy derivative of said each direct drive wind turbine generator in said each acquisition period is used as an energy flow of a dynamic coupling energy transfer path of said each direct drive wind turbine generator in said each acquisition period; obtaining energy flows of a generator internal coupling path, a generator-generator coupling path and a generator-grid coupling path by dividing the energy flow of the dynamic coupling energy transfer path of said each direct drive wind turbine generator in said each acquisition period, then obtaining an energy flow of each coupling energy transfer path of the direct drive wind farm in said each acquisition period by dividing the energy flows of the generator internal coupling path, the generator-generator coupling path, and the generator-grid coupling path.

3. The method according to claim 2, wherein the tracking the key energy transfer path of the direct drive wind farm is traced in the following ways: based on a set time interval, obtaining energy flow variation of each coupling energy transfer path of the direct drive wind farm in the set time interval by subtracting energy flow of a last acquisition period of said each coupling energy transfer path of the direct drive wind farm in the set time interval from energy flow of a first acquisition period; based on the energy flow variation of each coupling energy transfer path of the direct drive wind farm in the set time interval, obtaining respectively the coupling energy transfer paths with a largest energy flow variation in the energy flows of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path, and taking the coupling energy transfer paths obtained as the key energy transfer paths of the direct drive wind farm in the set time interval.

4. The method according to claim 2, wherein the energy flow of each coupling energy transfer path obtained by dividing the energy flow of the generator internal coupling path are respectively: $\{\begin{array}{c}{V}_{\mathrm{hi}1}=i_{drefi}{u}_{\mathrm{kdi}}^c\cos {}_i\\ V_{\mathrm{hi}2}=-i_{drefi}{u}_{\mathrm{kqi}}^c\sin {}_i\\ \begin{array}{c}{V}_{\mathrm{hi}3}=-\left(L_{\mathrm{xi}}+2L_T\right)\left(\frac{d{i}_{\mathrm{fdi}}^c\cos {}_i}{dt}{i}_{\mathrm{fdi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{dt}{i}_{\mathrm{fdi}}^c\sin {}_i\right)+\\ \left[R_{\mathrm{xi}}{i}_{\mathrm{fdi}}^c\sin {}_i+{}_0\left(L_{\mathrm{xi}}+L_T\right)i_{\mathrm{fdi}}^c\cos {}_i+\left(L_{\mathrm{xi}}+L_T\right)\frac{d{i}_{\mathrm{fdi}}^c\sin {}_i}{dt}\right]X_i\end{array}\\ \begin{array}{c}{V}_{\mathrm{hi}4}=-\left(L_{\mathrm{xi}}+2L_T\right)\left(\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{dt}{i}_{\mathrm{fqi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{dt}{i}_{\mathrm{fqi}}^c\sin {}_i\right)+\\ \left[R_{\mathrm{xi}}{i}_{\mathrm{fqi}}^c\cos {}_i-{}_0\left(L_{\mathrm{xi}}+L_T\right)i_{\mathrm{fqi}}^c\sin {}_i+\left(L_{\mathrm{xi}}+L_T\right)\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{dt}\right]X_i\\ \begin{array}{c}{V}_{\mathrm{hi}5}=-\left(L_{\mathrm{xi}}+2L_T\right)(-\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\sin {}_i-\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\cos {}_i\\ +\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\sin {}_i-\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\cos {}_i)\end{array}\end{array}\end{array}$ Where, V.sub.hi1 represents an energy flow of an energy transfer path of DC voltage outer loop, current loop d axis, and the PLL (phase locked loop) of an i-th wind turbine generator; V.sub.hi2 represents an energy flow of an energy transfer path of the DC voltage outer loop, current loop q axis, and the PLL of the i-th wind turbine generator; V.sub.hi3 represents an energy flow of an energy transfer path of the current loop d axis and the PLL of the i-th wind turbine generator; V.sub.hi4 represents an energy flow of an energy transfer path of the current loop q axis of the i-th wind turbine generator and the PLL of the i-th wind turbine generator, and V.sub.hi5 represents an energy flow of an energy transfer path of each of the current loop d axis and the current loop q axis of the i-th wind turbine generator and the PLL of the i-th wind turbine generator; i.sub.drefi is a d-axis current reference value of the current inner loop of the i-th wind turbine generator, $u_{\mathrm{kdi}}^c$ and $u_{\mathrm{kqi}}^c$ are respectively d and q axis components of the state variable of the current inner loop of the i-th wind turbine generator under a d-q coordinate system of a control system, .sub.i is a dynamic angle of the PLL of the i-th wind turbine generator, R.sub.xi and L.sub.xi are respectively a resistance and an inductance of a grid-connected line of the i-th wind turbine generator, $i_{\mathrm{fdi}}^c$ and $i_{\mathrm{fqi}}^c$ are d and q axis components of an output current of the i-th wind turbine generator under the d-q coordinate system of the control system, L.sub.T is an inductance of a transformer at a grid-connected bus bar outlet, .sub.0 is a grid synchronous angular frequency, and X.sub.i is a state variable of the PLL of the i-th wind turbine generator.

5. The method according to claim 4, wherein the energy flow of each coupling energy transfer path obtained by dividing the energy flow of the generator-generator coupling path are as follows: $\{\begin{array}{c}{V}_{\mathrm{ij}1}={}_0{L}_T{i}_{\mathrm{fdj}}^c\cos {}_j{X}_i+L_T\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\sin {}_j}{\mathrm{dt}}{X}_i\\ V_{\mathrm{ij}2}=-{}_0{L}_T{i}_{\mathrm{fqj}}^c\sin {}_j+L_T\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{X}_i\\ V_{\mathrm{ij}3}={}_0{L}_T{i}_{\mathrm{fdi}}^c\cos {}_i{X}_j+L_T\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{\mathrm{dt}}{X}_j\\ V_{\mathrm{ij}4}=-{}_0{L}_T{i}_{\mathrm{fqi}}^c\sin {}_i{X}_j+L_T\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{X}_j\\ \begin{array}{c}{V}_{\mathrm{ij}5}=-2L_T(\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\sin {}_i+\\ \frac{{\mathrm{di}}_{\mathrm{fdi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\cos {}_j+\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\sin {}_j)\\ \begin{array}{c}{V}_{\mathrm{ij}6}=-2L_T(-\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\sin {}_i-\\ \frac{{\mathrm{di}}_{\mathrm{fdi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\sin {}_j+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\cos {}_j)\\ \begin{array}{c}{V}_{\mathrm{ij}7}=-2L_T(-\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\sin {}_i+\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\cos {}_i-\\ \frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\cos {}_j+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\sin {}_j)\\ \begin{array}{c}{V}_{\mathrm{ij}8}=-2L_T(\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\sin {}_i+\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\cos {}_i+\\ \frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\sin {}_j+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\cos {}_j)\end{array}\end{array}\end{array}\end{array}\end{array}$ where V.sub.ij1 represents an energy flow of an energy transfer path of the phase locked loop of the i-th wind turbine generator, the current loop d axis of the i-th wind turbine generator, and a phase locked loop of a j-th wind turbine generator; V.sub.ij2 represents an energy flow of an energy transfer path of the phase locked loop of the i-th wind turbine generator, a current loop q axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij3 represents an energy flow of an energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij4 represents an energy flow of an energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij5 represents an energy flow of an energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, a current loop d axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij6 represents an energy flow of an energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, the current loop q axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij7 represents an energy flow of an energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, the current loop d axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij8 represents an energy flow of an energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, the current loop q axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; .sub.j is a dynamic angle of the phase locked loop of the j-th wind turbine generator, and $i_{\mathrm{fdj}}^c$ and $i_{\mathrm{fqj}}^c$ are respectively d and q axis components of an output current of the j-th wind turbine generator under the d-q coordinate system of the control system; X.sub.j is state variable of the phase locked loop of the j-th wind turbine generator.

6. The method according to claim 5, wherein the energy flow of each coupling energy transfer path obtained by dividing the energy flow of the generator-grid coupling path are respectively: $\{\begin{array}{c}{V}_{\mathrm{gi}1}=-2u_{\mathrm{md}}{i}_{fdi}^c\cos {}_i\\ V_{\mathrm{gi}2}=-2u_{\mathrm{mq}}{i}_{fdi}^c\sin {}_i\\ V_{gi3}=2u_{md}{i}_{\mathrm{fqi}}^c\sin {}_i\\ V_{gi4}=-2u_{mq}{i}_{fqi}^c\cos {}_i\\ V_{\mathrm{gi}5}=u_{mq}{X}_i\end{array}$ where V.sub.gi1 represents an energy flow of an energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and a d axis of a power grid; V.sub.gi2 represents an energy flow of an energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and a q axis of the power grid; V.sub.gi3 represents an energy flow of an energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the d axis of the power grid; V.sub.gi4 represents an energy flow of an energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the q axis of the power grid; V.sub.gi5 represents an energy flow of an energy transfer path of the phase locked loop of the i-th wind turbine generator and the q axis of the power grid; u.sub.md and u.sub.mq represent d and q axis components of the voltage of the fault point respectively.

7. The method according to claim 1, wherein the dominant oscillation mode period is $\frac{2}{{}_c};$ where .sub.c represents an oscillation frequency of the dominant oscillation mode of a system.

8. The method according to claim 6, wherein the corresponding instantaneous value of each oscillation component of said each direct drive wind turbine generator is expressed as: $\{\begin{array}{c}{i}_{fdi}^c=A_{di}{e}^t\sin \left({}_{c}t+{}_{di}\right)\\ i_{fqi}^c=A_{qi}{e}^t\sin \left({}_{c}t+{}_{qi}\right)\end{array}$ $i_{drefi}=A_{drefi}{e}^t\sin \left({}_{c}t+{}_{drefi}\right)$ $u_{\mathrm{kdi}}^c=A_{\mathrm{ukdi}}{e}^t\sin \left({}_{c}t+{}_{\mathrm{ukdi}}\right)$ ${}_i=A_i{e}^t\sin \left({}_{c}t+{}_i\right)$ $X_i=A_{xi}{e}^t\sin \left({}_{c}t+{}_{xi}\right)$ $\{\begin{array}{c}{u}_{md}=A_{md}{e}^t\sin \left({}_{c}t+{}_{md}\right)\\ u_{mq}=A_{mq}{e}^t\sin \left({}_{c}t+{}_{mq}\right)\end{array}$ where A.sub.di and A.sub.qi respectively represent amplitude of d and q axis components of an output current of the i-th wind turbine generator, is a system damping coefficient, .sub.c is an oscillation frequency under the d-q coordinate system of the control system, .sub.di and .sub.qi respectively represent phase of d and q axis components of the output current of the i-th wind turbine generator, and A.sub.drefi represents amplitude of the d-axis current reference value of the current inner loop of the i-th wind turbine generator; .sub.drefi is phase of the d-axis current reference value of the current inner loop of the i-th wind turbine generator, A.sub.ukdi is amplitude of a d-axis component of a state variable of the current inner loop of the i-th wind turbine generator, .sub.ukdi is phase of the d-axis component of the state variable of the current inner loop of the i-th wind turbine generator, A.sub.i is amplitude of the dynamic angle of the phase locked loop of the i-th wind turbine generator, and .sub.i is phase of the dynamic angle of the phase locked loop of the i-th wind turbine generator; A.sub.xi is amplitude of state variable of the phase locked loop of the i-th wind turbine generator, .sub.xi is phase of the state variable of the phase locked loop of the i-th wind turbine generator, A.sub.md and A.sub.mq are respectively amplitude of the d and q axis components of the voltage of the fault point, and .sub.md and .sub.mq are phase of the d and q axis components of the voltage of the fault point respectively.

9. A system for tracking oscillation path of a direct drive wind farm, comprising: a data acquisition module, configured to collect port current, current and state variable of current inner loop, dynamic angle and state component of a phase-locked loop of each direct drive wind turbine generator, and phase and amplitude of a voltage of a fault point in an acquisition period which is used as a dominant oscillation mode period when the direct drive wind farm fails, and then to obtain a corresponding instantaneous value of each oscillation component of said each direct drive wind turbine generator; a coupling energy transfer path acquisition module, configured to obtain an energy flow of a coupling energy transfer path of the direct drive wind farm in each acquisition period based on the corresponding instantaneous value of each oscillation component collected in said each acquisition period; a key energy transfer path tracking module, configured to track a key energy transfer path of the direct drive wind farm according to the energy flow of the coupling energy transfer path of the direct drive wind farm in said each acquisition period, so as to realize stability control of the direct drive wind farm.

10. The system according to claim 9, wherein the coupling energy transfer path acquisition module obtains the energy flow of the coupling energy transfer path between subsystems of the direct drive wind farm in said each acquisition period through the following ways: based on the corresponding instantaneous value of each oscillation component collected in said each acquisition period, a dynamic mathematical model of each subsystem of said each direct drive wind turbine generator in said each acquisition period is established, wherein the subsystems include current loop d, q axis subsystems, grid-connected line d, q axis subsystems, DC voltage outer loop subsystem and phase-locked loop subsystem; according to the dynamic mathematical model of each subsystem of said each direct drive wind turbine generator in said each acquisition period, a dynamic coupling energy derivative of each subsystem of said each direct drive wind turbine generator is obtained, and the dynamic coupling energy derivative of said each direct drive wind turbine generator in said each acquisition period is obtained; a value corresponding to the dynamic coupling energy derivative of said each direct drive wind turbine generator in said each acquisition period is used as an energy flow of a dynamic coupling energy transfer path of said each direct drive wind turbine generator in said each acquisition period; the energy flow of the dynamic coupling energy transfer path of said each direct drive wind turbine generator in said each acquisition period is divided to obtain an energy flow of a generator internal coupling path, a generator-generator coupling path and a generator-grid coupling path, then the energy flow of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path is divided respectively to obtain an energy flow of each coupling energy transfer path of the direct drive wind farm in said each acquisition period.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The appended drawings are used only for the purpose of showing specific embodiments and are not considered to be a limitation of the present application. Throughout the appended drawings, the same reference symbols represent the same components.

(2) FIG. 1 is a flow diagram of the direct drive wind farm oscillation path tracking method provided in example 1 of the application;

(3) FIG. 2 is a structural diagram of the grid-connected system of the direct drive wind farm provided in example 1 of the application;

(4) FIG. 3 is a coupling energy transfer path diagram of each subsystem in the direct drive wind turbine generator provided in example 1 of the application;

(5) FIG. 4 is a coupling energy transfer path diagram of each subsystem between the direct drive wind turbine generators provided in example 1 of the application;

(6) FIG. 5 is a coupling energy transfer path diagram between a direct drive wind turbine generator and a power grid subsystem for example 1 of the application;

(7) FIG. 6A is a graph of the energy flow V.sub.hi1 of the generator internal coupling energy transfer path for example 3 of the application;

(8) FIG. 6B is a graph of the energy flow V.sub.hi2 of the generator internal coupling energy transfer path provided for example 3 of the application;

(9) FIG. 6C is a diagram of the energy flow V.sub.hi3 of the generator internal coupling energy transfer path provided for example 3 of the application;

(10) FIG. 6D is a diagram of the energy flow V.sub.hi4 of the generator internal coupling energy transfer path provided for example 3 of the application;

(11) FIG. 6E is a graph of the energy flow V.sub.hi5 of the generator internal coupling energy transfer path for example 3 of the application;

(12) FIG. 7A is a graph of the energy flow V.sub.ij1 of the generator-generator coupling energy transfer path provided for example 3 of the application;

(13) FIG. 7B is a graph of the energy flow V.sub.ij2 of the generator-generator coupling energy transfer path provided for example 3 of the application;

(14) FIG. 7C is a graph of the energy flow V.sub.ij3 of the generator-generator coupling energy transfer path provided for example 3 of the application;

(15) FIG. 7D is a graph of the energy flow V.sub.ij4 of the generator-generator coupling energy transfer path provided for example 3 of the application;

(16) FIG. 7E is a graph of the energy flow V.sub.ij5 of the generator-generator coupling energy transfer path provided for example 3 of the application;

(17) FIG. 7F is a graph of the energy flow V.sub.ij6 of the generator-generator coupling energy transfer path provided for example 3 of the application;

(18) FIG. 7G is a graph of the energy flow V.sub.ij7 in the generator-generator coupling energy transfer path provided for example 3 of the application;

(19) FIG. 7H is a graph of the energy flow V.sub.ij8 in the generator-generator coupling energy transfer path provided for example 3 of the application;

(20) FIG. 8A is a graph of the energy flow V.sub.gi1 of the generator-grid coupling energy transfer path for example 3 of the application;

(21) FIG. 8B is a graph of the energy flow V.sub.gi2 of the generator-grid coupling energy transfer path for example 3 of the application;

(22) FIG. 8C is a graph of the energy flow V.sub.gi3 of the generator-grid coupling energy transfer path for example 3 of the application;

(23) FIG. 8D is a graph of the energy flow V.sub.gi4 of the generator-grid coupling energy transfer path for example 3 of the application;

(24) FIG. 8E is a graph of the energy flow V.sub.gi5 of the generator-grid coupling energy transfer path for example 3 of the application.

DETAILED DESCRIPTION

(25) The proposed method is described in detail below in conjunction with the diagrams in the appendix. It should be emphasized that the following description is only exemplary and is not intended to limit the scope of this invention and its application.

Example 1

(26) A specific embodiment of the application discloses a method for tracking the oscillation path of a direct drive wind farm, as shown in FIG. 1, including the following steps:

(27) S1: When a direct drive wind farm fails, the dominant oscillation mode period is taken as the acquisition period to collect the port current, current and state variables in the current inner loop, dynamic angle and state component of phase locked loop, and the phase and amplitude of the fault point voltage of each wind turbine generator, and then the corresponding instantaneous value of each oscillation component of the direct drive wind turbine generator is obtained.

(28) When implemented, the instantaneous value of each oscillation component of the direct drive wind turbine generator is expressed as:

(29) $$\begin{gathered} \{\begin{array}{c}{i}_{\mathrm{fdi}}^c=A_{di}{e}^t\sin \left({}_{c}t+{}_{di}\right)\\ i_{\mathrm{fqi}}^c=A_{qi}{e}^t\sin \left({}_{c}t+{}_{qi}\right)\end{array} \\ {i}_{\mathrm{drefi}}=A_{\mathrm{drefi}}{e}^t\sin \left({}_{c}t+{}_{\mathrm{drefi}}\right) \\ {u}_{\mathrm{kdi}}^c=A_{\mathrm{ukdi}}{e}^t\sin \left({}_{c}t+{}_{\mathrm{ukdi}}\right) \\ {}_i=A_i{e}^t\sin \left({}_{c}t+{}_i\right) \\ {X}_i=A_{xi}{e}^t\sin \left({}_{c}t+{}_{xi}\right)\{\begin{array}{c}{u}_{md}=A_{md}{e}^t\sin \left({}_{c}t+{}_{md}\right)\\ u_{mq}=A_{mq}{e}^t\sin \left({}_{c}t+{}_{mq}\right)\end{array} \end{gathered}$$ where,

(30) $i_{\mathrm{fdi}}^c$
and

(31) $i_{\mathrm{fqi}}^c$
are respectively the d and q axis components of the output current of the i-th wind turbine generator under the d-q coordinate system of the control system; A.sub.di and A.sub.qi are respectively the amplitudes of the d and q axis components of the output current of the i-th wind turbine generator; is the system damping coefficient; .sub.c is the oscillation frequency under the d-q coordinate system of the control system; .sub.di and .sub.qi are respectively the phases of the d and q components of the output current of the i-th wind turbine generator; i.sub.drefi is the d-axis current reference value of the current inner loop of the i-th wind turbine generator, A.sub.drefi is the amplitude of the d-axis current reference value of the current inner loop of the i-th wind turbine generator, .sub.drefi is the phase of the d axis current reference value of the current inner loop of the i-th wind turbine generator, and

(32) $u_{\mathrm{kdi}}^c$
is the d axis component of the state variable of the current inner loop of the i-th wind turbine generator under the d-q coordinate system of the control system; A.sub.ukdi is the amplitude of the d-axis component of the state variable of the current internal loop of the i-th wind turbine generator, .sub.ukdi is the phase of the d-axis component of the state variable of the current internal loop of the i-th wind turbine generator, .sub.i is the dynamic angle of the phase locked loop of the i-th wind turbine generator, A.sub.i is the amplitude of the dynamic angle of the phase locked loop of the i-th wind turbine generator, .sub.i is the phase of the dynamic angle of the phase locked loop of the i-th wind turbine generator, and X.sub.i is the state variable of the phase locked loop of the i-th wind turbine generator. A.sub.xi is the amplitude of the state variable of the phase locked loop of the i-th wind turbine generator, .sub.xi is the phase of the state variable of the phase locked loop of the i-th wind turbine generator, u.sub.md and u.sub.mq are respectively the d and q axis components of the fault point voltage, A.sub.md and A.sub.mq are respectively the amplitudes of the d and q axis components of the fault point voltage, .sub.md and .sub.mq are respectively the phases of the d and q axis components of the fault point voltage. It should be noted that the parameter superscript c represents the variable under the d-q coordinate system of the control system, the absence of superscript c represents the variable under the d-q coordinate system of the power grid, and the subscript i represents the i-th wind turbine generator.

(33) Specifically, data are collected by PMU devices installed on the ports of each wind turbine generator in a direct drive wind farm.

(34) More specifically, the dominant oscillation mode period is

(35) $\frac{2}{{}_c};$
where, .sub.c represents the oscillation frequency of the dominant oscillation mode of the system.

(36) S2: Based on the instantaneous value of each oscillation component collected in each acquisition period, the energy flow of the coupling energy transfer path of the direct drive wind farm in each acquisition period is obtained.

(37) When implemented, the energy flow of the coupling energy transfer path between subsystems of the direct drive wind farm in each acquisition period is obtained by the following ways:

(38) S21: Based on the instantaneous value of each oscillation component collected in each acquisition period, the dynamic mathematical model of each subsystem of each direct drive wind turbine generator in the direct drive wind farm in each acquisition period is established; the subsystems include current loop d and q axis subsystems, grid-connected line d and q axis subsystems, DC voltage outer loop subsystem and phase-locked loop subsystem.

(39) Specifically, as shown in FIG. 2, each direct drive wind turbine is divided into 6 subsystems according to each control link, which are current loop d, q axes subsystems, grid-connected line d, q axes subsystems, DC voltage outer loop subsystem and phase-locked loop subsystem. The dynamic mathematical models of each subsystem are as follows:

(40) The dynamic mathematical models of the current loop d and q axis subsystems are respectively:

(41) $$\begin{gathered} \{\begin{array}{c}\frac{1}{k_{\mathrm{Ii}}}\frac{d{u}_{\mathrm{kdi}}}{dt}=-i_{\mathrm{fdi}}+i_{\mathrm{drefi}}\\ L_f\frac{d{i}_{\mathrm{fdi}}}{dt}=-\left(k_{Pi}+R_f\right)i_{\mathrm{fdi}}+u_{\mathrm{kdi}}+k_{Pi}{i}_{\mathrm{drefi}}+{}_0{L}_f{i}_{\mathrm{fqi}}-u_{gdi}\end{array} \\ \{\begin{array}{c}\frac{1}{k_{\mathrm{Ii}}}\frac{d{u}_{\mathrm{kqi}}}{dt}=-i_{\mathrm{fqi}}+i_{\mathrm{qrefi}}\\ L_f\frac{d{i}_{\mathrm{fqi}}}{dt}=-\left(k_{Pi}+R_f\right)i_{\mathrm{fqi}}+u_{\mathrm{kqi}}+k_{Pi}{i}_{\mathrm{qrefi}}-{}_0{L}_f{i}_{\mathrm{fdi}}-u_{\mathrm{gqi}}\end{array} \end{gathered}$$ where, k.sub.Pi and k.sub.Ii are the proportion and integral coefficient of the inner current loop of the i-th wind turbine generator; R.sub.f and L.sub.f are the resistance and inductance of the filtering link respectively; i.sub.fdi and i.sub.fqi are the d and q axis components of the output current of the i-th wind turbine generator respectively; u.sub.kdi and u.sub.kqi are the d and q axis components of the state variable of the inner current loop of the i-th wind turbine generator respectively. i.sub.drefi and i.sub.qrefi are respectively the d and q axis current reference values of the current inner loop of the i-th wind turbine generator, .sub.0 is the synchronous angular frequency of the grid, and u.sub.gdi and u.sub.gqi are respectively the d and q axis components of the PLL sampling point voltage of the i-th wind turbine generator. The sampling point is a bus bar located at the wind turbine port.

(42) The dynamic mathematical models of grid-connected line d and q axis subsystems are as follows:

(43) $L_{\mathrm{xi}}\frac{d{i}_{\mathrm{fdi}}}{dt}=-R_{\mathrm{xi}}{i}_{\mathrm{fdi}}+{}_0{L}_{\mathrm{xi}}{i}_{\mathrm{fqi}}+u_{\mathrm{fdi}}-u_{sd}{L}_{\mathrm{xi}}\frac{d{i}_{\mathrm{fqi}}}{dt}=-R_{\mathrm{xi}}{i}_{\mathrm{fqi}}-{}_0{L}_{\mathrm{xi}}{i}_{\mathrm{fdi}}+u_{\mathrm{fqi}}-u_{sq}$ where, R.sub.xi and L.sub.xi are respectively the resistance and inductance of the grid-connected line of the i-th wind turbine generator, and u.sub.sd and u.sub.sq are respectively the d and q axis components of the system PCC point voltage.

(44) The dynamic mathematical model of the DC voltage outer loop subsystem is as follows:

(45) $C_f{u}_{dci}\frac{d{u}_{dci}}{dt}=-\left(u_{\mathrm{fdi}}{i}_{\mathrm{fdi}}+u_{\mathrm{fqi}}{i}_{\mathrm{fqi}}\right)+P_m$ where, C.sub.f is the DC capacitance, u.sub.dci is the DC capacitance voltage of the i-th wind turbine generator, u.sub.fdi and u.sub.fqi are the d and q axis components of the terminal voltage of the i-th wind turbine generator respectively, and P.sub.m is the mechanical power transmitted by the converter at the side of the generator.

(46) The dynamic mathematical model of the phase-locked loop subsystem is:

(47) $\{\begin{array}{c}\frac{1}{k_{Ii}}\frac{{\mathrm{dX}}_i}{dt}=u_{\mathrm{gqi}}\\ \frac{d{}_{\mathrm{pll}i}}{dt}=X_i+{}_0+k_{Pi}{u}_{\mathrm{gqi}}\end{array}$ where, k.sub.Pi and k.sub.Ii are the proportion and integral coefficient of the phase locked loop of the i-th wind turbine generator respectively, and .sub.plli is the output phase locked angle of the i-th wind turbine generator.

(48) S22: According to the dynamic mathematical model of each subsystem of each direct drive wind turbine generator in the direct drive wind farm in each acquisition period, the dynamic coupling energy derivative of each subsystem of each direct drive wind turbine generator is obtained, and the dynamic coupling energy derivative of each generator in each acquisition period is obtained, and the corresponding value is used as the energy flow of the dynamic coupling energy transfer path of each generator in each acquisition period.

(49) Specifically, the dynamic coupling energy derivatives of each subsystem of each direct drive wind turbine generator are:

(50) The dynamic coupling energy derivatives

(51) 0$W_{\mathrm{zfdi}}^{\&}$
and

(52) $W_{\mathrm{zfqi}}^{\&}$
of the current loop d and q axis subsystems are respectively expressed as:

(53) $W_{\mathrm{zfdi}}^{\&}=i_{\mathrm{drefi}}{u}_{\mathrm{kdi}}+\left(k_{\mathrm{Pi}}{i}_{\mathrm{drefi}}+{}_0{L}_f{i}_{\mathrm{fqi}}-u_{\mathrm{gdi}}\right)i_{\mathrm{fdi}}{W}_{\mathrm{zfqi}}^{\&}=i_{\mathrm{qrefi}}{u}_{\mathrm{kqi}}+\left(k_{\mathrm{Pi}}{i}_{\mathrm{qrefi}}+{}_0{L}_f{i}_{\mathrm{fdi}}-u_{\mathrm{gqi}}\right)i_{\mathrm{fqi}}$

(54) The dynamic coupling energy derivatives {dot over (W)}.sub.zldi and {dot over (W)}.sub.zlqi of the grid-connected line d and q axis subsystems are respectively expressed as:

(55) ${\stackrel{}{W}}_{zldi}=\left(u_{gdi}-u_{sd}+{}_0{L}_{xi}{i}_{\mathrm{fqi}}\right)i_{\mathrm{fdi}}{\stackrel{}{W}}_{zlqi}=\left(u_{\mathrm{gqi}}-u_{sq}-{}_0{L}_{xi}{i}_{\mathrm{fdi}}\right)i_{\mathrm{fqi}}$

(56) The dynamic coupling energy derivative {dot over (W)}.sub.zdci of the DC voltage outer loop subsystem is expressed as:

(57) $W_{zdci}=-k_{Pi}{i}_{\mathrm{drefi}}{i}_{\mathrm{fdi}}+k_{Pi}{i}_{\mathrm{fdi}}{i}_{\mathrm{fdi}}-i_{\mathrm{fdi}}{u}_{\mathrm{kdi}}-u_{gdi}{i}_{\mathrm{fdi}}-k_{Pi}{i}_{\mathrm{qrefi}}{i}_{\mathrm{fqi}}+k_{Pi}{i}_{\mathrm{fqi}}{i}_{\mathrm{fqi}}-i_{\mathrm{fqi}}{u}_{kqi}-u_{gqi}{i}_{\mathrm{fqi}}+P_m$

(58) The dynamic coupling energy derivative

(59) $W_{zi}^{\&}$
of the phase-locked loop subsystem is expressed as:
W.sup.&.sub.zi=u.sub.gqiX.sub.i

(60) Based on this, the dynamic coupling energy derivative

(61) $W_{\mathrm{zi}}^{\&}$
of the i-th wind turbine generator is obtained, that is, the energy flow of the dynamic coupling energy transfer path of the i-th wind turbine generator, expressed as:

(62) $W_{zi}^{\&}=W_{\mathrm{zfdi}}^{\&}+W_{\mathrm{zfqi}}^{\&}+W_{zldi}^{\&}+W_{zlqi}^{\&}+W_{zdci}^{\&}+W_{zi}^{\&}=\left(i_{\mathrm{drefi}}-i_{\mathrm{fdi}}\right)u_{kdi}+\left(i_{\mathrm{qrefi}}-i_{\mathrm{fqi}}\right)u_{kqi}+k_{Pi}\left(i_{\mathrm{fdi}}^2+i_{\mathrm{fqi}}^2\right)-\left(u_{gdi}+u_{sd}\right)i_{\mathrm{fdi}}-\left(u_{gqi}+u_{sq}\right)i_{\mathrm{fqi}}+u_{gqi}{X}_i+P_m$ therein,

(63) $u_{\mathrm{gdi}}=R_{\mathrm{xi}}{i}_{\mathrm{fdi}}+L_{\mathrm{xi}}\frac{d{i}_{\mathrm{fdi}}}{dt}-{}_0{L}_{\mathrm{xi}}{i}_{\mathrm{fqi}}+L_T\frac{d\left(i_{\mathrm{fdi}}+i_{\mathrm{fdj}}\right)}{dt}-{}_0{L}_T\left(i_{\mathrm{fqi}}+i_{\mathrm{fqj}}\right)+u_{\mathrm{md}}$$u_{\mathrm{gqi}}=R_{\mathrm{xi}}{i}_{\mathrm{fqi}}+L_{\mathrm{xi}}\frac{d{i}_{\mathrm{fqi}}}{dt}+{}_0{L}_{\mathrm{xi}}{i}_{\mathrm{fdi}}+L_T\frac{d\left(i_{\mathrm{fqi}}+i_{\mathrm{fqj}}\right)}{dt}+{}_0{L}_T\left(i_{\mathrm{fdi}}+i_{fdj}\right)+u_{mq}$

(64) Where, L.sub.T represents the inductance of the transformer at the grid-connected bus bar outlet.

(65) It should be noted that the object of the transfer path is determined according to the control link corresponding to the variable contained in the derivative expression, so that the transfer path is represented by the derivative expression, and the specific value of the derivative obtained is the energy flow of the transfer path.

(66) It should be noted that the dynamic coupling energy derivative basis of each subsystem is obtained based on the following analysis:

(67) The dynamic mathematical model of the 6 subsystems can be uniformly expressed as follows:

(68) $\{\begin{array}{c}A\frac{dx}{dt}=-Cy+H_x\\ B\frac{dy}{dt}=-Dy+Ex+H_y\end{array}$

(69) In the formula, A, B, C, D, E, H, and H.sub.y represent coefficients respectively, and x and y are variables respectively.

(70) The first integration method was used to construct the energy function, and the two formulas above were cross-multiplied and integrated to obtain the dynamic energy model corresponding to each subsystem as follows:

(71) 0$\frac{1}{2}AE{x}^2+\frac{1}{2}BC{y}^2+CD{y}^{2}dt-C{H}_{y}ydt-E{H}_{x}x\mathrm{dt}=\mathrm{Const}$ where Const represents the integral constant, AEx.sup.2/2 and BCy.sup.2/2 correspond to the dynamic storage energy injected into the corresponding control links of x and y in the subsystem after the fault, denoted as W.sub.s; CDy.sup.2dt corresponds to the dynamic dissipated energy injected into the corresponding control link of y in the subsystem after the fault, and is denoted as W.sub.d. CH.sub.yydt and EH.sub.xxdt correspond to the dynamic coupling energy generated between the corresponding control links of x and y and other control links injected into the subsystem after the fault, and are denoted as W.sub.z. It is known that they satisfy

(72) $W_s^{\&}=W_d^{\&}+W_z^{\&}$

(73) The response process of wind farm after disturbance is essentially a process of energy storage, dissipation and coupling among all links, which follows the law of energy conservation. According to Lyapunov's second law, the state variables of a disturbed system will deviate from the equilibrium point and obtain a certain dynamic storage energy. If the dynamic storage energy gradually decays to the minimum value over time, the system can gradually recover stability in this state. On the contrary, if the dynamic storage energy of the system continues to increase, the system is unstable. The dynamic dissipative energy is always negative, which is conducive to the stability of the system, and the dynamic coupling energy is the direct reason to determine the stability of the system. When the energy derivative is positive, the stored energy continues to increase, and the system stability decreases, and the larger the value, the worse the system stability. Therefore, the variation of the dynamic coupling energy derivative directly reflects the stability of the system.

(74) The expression of the dynamic coupling energy derivative obtained from the above equation is as follow:

(75) $W_z^{\&}={\mathrm{EH}}_{x}x+{\mathrm{CH}}_{y}y$

(76) Therefore, the dynamic coupling energy derivative can be obtained according to the dynamic mathematical model of each subsystem of each direct drive wind turbine generator.

(77) S23: The energy flow of the dynamic coupling energy transfer path of each generator in each acquisition period is divided to obtain the energy flow of the generator internal coupling path the generator-generator coupling path and the generator-grid coupling path. Then the energy flow of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path is divided respectively to obtain the energy flow of each coupling energy transfer path of the direct drive wind farm in each acquisition period.

(78) Specifically, the generator internal coupling path of each wind turbine generator contains only the relevant variables of the generator. The generator-generator coupling path of each wind turbine generator contains the variables of the wind turbine generator and any other wind turbine generator. The generator-grid coupling path of each wind turbine generator includes the variables of the generator and the power grid.

(79) Therefore, based on the expression of the dynamic coupling energy transfer path

(80) $W_{\mathrm{zi}}^{\&}$
of each generator, the energy flow of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path of each turbine generator is obtained, which is expressed as follows:

(81) The energy flow of the generator internal coupling path is:

(82) $V_{\mathrm{hi}}=i_{dref}\left(u_{\mathrm{kdi}}^c\cos {}_i-u_{\mathrm{kqi}}^c\sin {}_i\right)-i_{\mathrm{fdi}}^c{u}_{\mathrm{kdi}}^c-i_{\mathrm{fqi}}^c{u}_{\mathrm{kqi}}^c+\left(k_{\mathrm{Pi}}-R_{\mathrm{xi}}\right)\left(i_{\mathrm{fdi}}^{c2}+i_{\mathrm{fqi}}^{c2}\right)-\left(L_{\mathrm{xi}}+2L_T\right)\left(\begin{array}{c}\frac{d\left(i_{\mathrm{fdi}}^c\cos {}_i-i_{\mathrm{fqi}}^c\sin {}_i\right)}{dt}\left(i_{\mathrm{fdi}}^c\cos {}_i-i_{\mathrm{fqi}}^c\sin {}_i\right)+\\ \frac{d\left(i_{\mathrm{fqi}}^c\cos {}_i+i_{\mathrm{fdi}}^c\sin {}_i\right)}{dt}\left(i_{\mathrm{fqi}}^c\cos {}_i+i_{\mathrm{fdi}}^c\sin {}_i\right)\end{array}\right)+\left[\begin{array}{c}{R}_{\mathrm{xi}}\left(i_{\mathrm{fqi}}^c\cos {}_i+i_{\mathrm{fdi}}^c\sin {}_i\right)+\left(L_{\mathrm{xi}}+L_T\right)\frac{d\left(i_{\mathrm{fqi}}^c\cos {}_i+i_{\mathrm{fdi}}^c\sin {}_i\right)}{\mathrm{dt}}+\\ {}_0\left(L_{\mathrm{xi}}+L_T\right)\left(i_{\mathrm{fdi}}^c\cos {}_i-i_{\mathrm{fqi}}^c\sin {}_i\right)\end{array}\right]X_i$ where, V.sub.hi represents the energy flow of the coupling path in the generator.

(83) The energy flow of the generator-generator coupling path is:

(84) $V_{\mathrm{ti}-j}=-2L_T\left[\begin{array}{c}\left(\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\cos {}_j}{dt}-\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\sin {}_j}{dt}\right)\left(i_{\mathrm{fdi}}^c\cos {}_i-i_{\mathrm{fqi}}^c\sin {}_i\right)+\\ \left(\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{dt}+\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\sin {}_j}{dt}\right)\left(i_{\mathrm{fqi}}^c\cos {}_i+i_{\mathrm{fdi}}^c\sin {}_i\right)\\ +{}_0\left(i_{\mathrm{fdj}}^c\cos {}_j-i_{\mathrm{fdj}}^c\sin {}_j\right)\left(i_{\mathrm{fqi}}^c\cos {}_i+i_{\mathrm{fdi}}^c\sin {}_i\right)\\ -{}_0\left(i_{\mathrm{fdj}}^c\cos {}_j+i_{\mathrm{fdj}}^c\sin {}_j\right)\left(i_{\mathrm{fdi}}^c\cos {}_i-i_{\mathrm{fqi}}^c\sin {}_i\right)\end{array}\right]+\left[+L_T\frac{d\left(i_{\mathrm{fqj}}^c\cos {}_j+i_{fdj}^c\sin {}_j\right)}{dt}+{}_{0}L\left(i_{fdj}^c\cos {}_j-i_{\mathrm{fqj}}^c\sin {}_j\right)\right]X_i$ where, V.sub.tij represents the energy flow of the generator-generator coupling path.

(85) The energy flow of the generator-grid coupling path is:

(86) $V_{\mathrm{gi}}=-2u_{\mathrm{md}}\left(i_{\mathrm{fdi}}^c\cos {}_i-i_{\mathrm{fqi}}^c\sin {}_i\right)-2u_{mq}\left(i_{\mathrm{fqi}}^c\cos {}_i+i_{\mathrm{fdi}}^c\sin {}_i\right)+u_{mq}{X}_i$ where V.sub.gi represents the energy flow of the generator-grid coupling path.

(87) In the same way, the energy flow of the related energy coupling transfer path of the j-th wind turbine generator in the direct drive wind farm can be obtained, ij.

(88) Specifically, the energy flow of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path is divided, and the energy flow of the coupling energy transfer path of the direct drive wind farm in each acquisition period is as follows:

(89) As shown in FIG. 3, the energy flows of each coupling energy transfer path obtained by dividing the generator internal coupling path are as follows:

(90) $\{\begin{array}{c}{V}_{\mathrm{hi}1}=i_{drefi}{u}_{\mathrm{kdi}}^c\cos {}_i\\ V_{\mathrm{hi}2}=-i_{drefi}{u}_{\mathrm{kqi}}^c\sin {}_i\\ \begin{array}{c}{V}_{\mathrm{hi}3}=-\left(L_{\mathrm{xi}}+2L_T\right)\left(\frac{d{i}_{\mathrm{fdi}}^c\cos {}_i}{dt}{i}_{\mathrm{fdi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{dt}{i}_{\mathrm{fdi}}^c\sin {}_i\right)+\\ \left[R_{\mathrm{xi}}{i}_{\mathrm{fdi}}^c\sin {}_i+{}_0\left(L_{\mathrm{xi}}+L_T\right)i_{\mathrm{fdi}}^c\cos {}_i+\left(L_{\mathrm{xi}}+L_T\right)\frac{d{i}_{\mathrm{fdi}}^c\sin {}_i}{dt}\right]X_i\end{array}\\ \begin{array}{c}{V}_{\mathrm{hi}4}=-\left(L_{\mathrm{xi}}+2L_T\right)\left(\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{dt}{i}_{\mathrm{fqi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{dt}{i}_{\mathrm{fqi}}^c\sin {}_i\right)+\\ \left[R_{\mathrm{xi}}{i}_{\mathrm{fqi}}^c\cos {}_i-{}_0\left(L_{\mathrm{xi}}+L_T\right)i_{\mathrm{fqi}}^c\sin {}_i+\left(L_{\mathrm{xi}}+L_T\right)\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{dt}\right]X_i\\ V_{\mathrm{hi}5}=-\left(L_{\mathrm{xi}}+2L_T\right)\left(\begin{array}{c}-\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\sin {}_i-\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\cos {}_i\\ +\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\sin {}_i-\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\cos {}_i\end{array}\right)\end{array}\end{array}$ where V.sub.hi1 represents the energy flow of the energy transfer path of the DC voltage outer loop, current loop d axis, and PLL of the i-th wind turbine generator; V.sub.hi2 represents the energy flow of the energy transfer path of the DC voltage outer loop, current loop q axis, and PLL of the i-th wind turbine generator; V.sub.hi3 represents the energy flow of the energy transfer path of the current loop d axis and PLL of the i-th wind turbine generator; V.sub.hi4 represents the energy flow of the energy transfer path of the current loop q axis of the i-th wind turbine generator and the phase locked loop of the i-th wind turbine generator, and V.sub.hi5 represents the energy flow of the energy transfer path of the current loop d axis, q axis of the i-th wind turbine generator, and the phase locked loop of the i-th wind turbine generator.

(91) As shown in FIG. 4, the energy flows of each coupling energy transfer path obtained by dividing the generator-generator coupling paths are as follows:

(92) $\{\begin{array}{c}{V}_{\mathrm{ij}1}={}_0{L}_T{i}_{\mathrm{fdj}}^c\cos {}_j{X}_i+L_T\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\sin {}_j}{\mathrm{dt}}{X}_i\\ V_{\mathrm{ij}2}=-{}_0{L}_T{i}_{\mathrm{fqj}}^c\sin {}_j+L_T\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{X}_i\\ V_{\mathrm{ij}3}={}_0{L}_T{i}_{\mathrm{fdi}}^c\cos {}_i{X}_j+L_T\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{\mathrm{dt}}{X}_j\\ V_{\mathrm{ij}4}=-{}_0{L}_T{i}_{\mathrm{fqi}}^c\sin {}_i{X}_j+L_T\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{X}_j\\ \begin{array}{c}{V}_{\mathrm{ij}5}=-2L_T(\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\sin {}_i+\\ \frac{{\mathrm{di}}_{\mathrm{fdi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\cos {}_j+\frac{{\mathrm{di}}_{\mathrm{fdi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\sin {}_j)\\ \begin{array}{c}{V}_{\mathrm{ij}6}=-2L_T(-\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\cos {}_i+\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fdi}}^c\sin {}_i-\\ \frac{{\mathrm{di}}_{\mathrm{fdi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\sin {}_j+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\cos {}_j)\\ \begin{array}{c}{V}_{\mathrm{ij}7}=-2L_T(-\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\sin {}_i+\frac{{\mathrm{di}}_{\mathrm{fdj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\cos {}_i-\\ \frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\cos {}_j+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fdj}}^c\sin {}_j)\\ \begin{array}{c}{V}_{\mathrm{ij}8}=-2L_T(\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\sin {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\sin {}_i+\frac{{\mathrm{di}}_{\mathrm{fqj}}^c\cos {}_j}{\mathrm{dt}}{i}_{\mathrm{fqi}}^c\cos {}_i+\\ \frac{{\mathrm{di}}_{\mathrm{fqi}}^c\sin {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\sin {}_j+\frac{{\mathrm{di}}_{\mathrm{fqi}}^c\cos {}_i}{\mathrm{dt}}{i}_{\mathrm{fqj}}^c\cos {}_j)\end{array}\end{array}\end{array}\end{array}\end{array}$ where V.sub.ij1 represents the energy flow of the energy transfer path of the phase locked loop of the i-th wind turbine generator, the current loop d axis of the i-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij2 represents the energy flow of the energy transfer path of the phase locked loop of the i-th wind turbine generator, the current loop q axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij3 represents the energy flow of the energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij4 represents the energy flow of the energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij5 represents the energy flow of the energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, the current loop d axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij6 represents the energy flow of the energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, the current loop q axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij7 represents the energy flow of the energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, the current loop d axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator; V.sub.ij8 represents the energy flow of the energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, the current loop q axis of the j-th wind turbine generator, and the phase locked loop of the j-th wind turbine generator.

(93) As shown in FIG. 5, the energy flows of each coupling energy transfer path obtained by dividing the generator-grid coupling paths are as follows:

(94) $\{\begin{array}{c}{V}_{\mathrm{gi}1}=-2u_{\mathrm{md}}{i}_{fdi}^c\cos {}_i\\ V_{\mathrm{gi}2}=-2u_{\mathrm{mq}}{i}_{fdi}^c\sin {}_i\\ V_{gi3}=2u_{md}{i}_{\mathrm{fqi}}^c\sin {}_i\\ V_{gi4}=-2u_{mq}{i}_{fqi}^c\cos {}_i\\ V_{\mathrm{gi}5}=u_{mq}{X}_i\end{array}$ where V.sub.gi1 represents the energy flow of the energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the d axis of the grid; V.sub.gi2 represents the energy flow of the energy transfer path of the current loop d axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the q axis of the grid; V.sub.gi3 represents the energy flow of the energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the d axis of the grid; V.sub.gi4 represents the energy flow of the energy transfer path of the current loop q axis of the i-th wind turbine generator, the phase locked loop of the i-th wind turbine generator, and the q axis of the grid; V.sub.gi5 represents the energy flow of the energy transfer path of the phase locked loop of the i-th wind turbine generator and the q axis of the grid.

(95) S3: According to the energy flow of the coupling energy transfer path of the direct drive wind farm in each acquisition period, the key energy transfer path of the direct drive wind farm in each acquisition period is tracked, and the stability control of the direct drive wind farm is realized.

(96) When implemented, key energy transfer paths of direct drive wind farms are tracked by:

(97) Based on the set time interval, the energy flow of the last acquisition period of each coupling energy transfer path of the direct drive wind farm in this time interval is subtracted from the energy flow of the first acquisition period, and the energy flow variation of each coupling energy transfer path of the direct drive wind farm in this time interval is obtained.

(98) Based on the energy flow variation of each coupling energy transfer path of a direct drive wind farm in the set time interval, the coupling energy transfer paths with the largest energy flow variation in the energy flows of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path are obtained respectively, and the obtained coupling energy transfer paths are taken as the key energy transfer paths of the direct drive wind farm in the time interval.

(99) Preferably, the set time interval is two adjacent acquisition periods, and the key energy transfer path of the direct drive wind farm for each acquisition period can be obtained.

(100) It can be understood that the coupling energy transfer path is the corresponding derivative, and the specific value of the derivative represents the increase or decrease of the energy flow of the coupling energy transfer path. By comparing the specific value of each coupling energy transfer path, that is, the increase or decrease of the energy flow, the key energy transfer path affecting the stability of the direct drive wind farm is tracked, and the key control link affecting the stability level of the system is selected out.

(101) Further, by tracking the key energy transfer path, the adverse control link in the direct drive wind farm is defined, and the oscillation is suppressed by adding the energy compensation branch in the corresponding control link, so as to achieve the stability control of the direct drive wind farm.

(102) Compared with the prior art, this example provides a method for tracking the oscillation path of a direct drive wind farm. By taking the dominant oscillation mode period as the acquisition period when a direct drive wind farm fails, the corresponding instantaneous value of each oscillation component is obtained by collecting data, and the energy flow of the coupling energy transfer path between subsystems of the direct drive wind farm in each acquisition period is obtained. Then, the key energy transfer path of the direct drive wind farm in each acquisition period is tracked to achieve accurate and effective tracking of the key energy transfer path of the grid-connected system of the direct drive wind farm, and guiding suggestions are given for the stable and safe operation of the power system. And by dividing the whole wind turbine generator into multiple subsystems, the key control links that affect the stability level of the system are selected out. The coupling energy transfer paths obtained by dynamic coupling energy derivatives reduce the calculation amount, improve the calculation accuracy, and ensure the reliable operation of the system.

Example 2

(103) A specific embodiment of the application discloses a direct drive wind farm oscillation path tracking system, which comprises:

(104) The data acquisition module is used to collect the port current, the current and state variable of the current inner loop, the dynamic angle and state component of the phase-locked loop, and the phase and amplitude of the fault point voltage of each wind turbine generator when a direct drive wind farm fails, and then to obtain the corresponding instantaneous value of each oscillation component of the direct drive wind turbine generator.

(105) The coupling energy transfer path acquisition module is used to obtain the energy flow of the coupling energy transfer path of the direct drive wind farm in each acquisition period based on the instantaneous value of each oscillation component collected in each acquisition period.

(106) The key energy transfer path tracking module is used to track the key energy transfer path of the direct drive wind farm in each acquisition period according to the energy flow of the coupling energy transfer path of the direct drive wind farm in each acquisition period, and then to realize the stability control of the direct drive wind farm.

(107) When implemented, the coupling energy transfer path acquisition module obtains the energy flow of the coupling energy transfer path between subsystems of the direct drive wind farm in each acquisition period through the following ways:

(108) Based on the instantaneous value of each oscillation component collected in each acquisition period, the dynamic mathematical model of each subsystem of each direct drive wind turbine generator in each acquisition period is established. The subsystems include current loop d, q axis subsystems, grid-connected line d, q axis subsystems, DC voltage outer loop subsystem and phase-locked loop subsystem.

(109) According to the dynamic mathematical model of each subsystem of each direct drive wind turbine generator in each acquisition period, the dynamic coupling energy derivative of each subsystem of each direct drive wind turbine generator is obtained, and the dynamic coupling energy derivative of each generator in each acquisition period is obtained. The corresponding value is used as the energy flow of the dynamic coupling energy transfer path of each generator in each acquisition period.

(110) The energy flow of the dynamic coupling energy transfer path of each generator in each acquisition period is divided to obtain the energy flow of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path. Then the energy flow of the generator internal coupling path, the generator-generator coupling path and the generator-grid coupling path is divided respectively to obtain the energy flow of each coupling energy transfer path of the direct drive wind farm in each acquisition period.

(111) The specific implementation process of the embodiment of the application can be referred to the above method embodiment, and the embodiment will not be repeated here.

(112) Since the principle of this embodiment is the same as that of the above method embodiments, the system also has the corresponding technical effects of the above method embodiments.

Example 3

(113) In order to verify the correctness and feasibility of the direct drive wind farm oscillation path tracking method and system provided by examples 1 and 2 of the application, the scheme in the above examples is tested and verified. The grid-connected system of two direct drive wind farms is built from the simulation level, and the impedance of the AC system is changed to obtain the oscillating divergence scene after the system is disturbed. The system simulation structure diagram is shown in FIG. 2. The simulation model is built on the RT-LAB platform. The parameters of grid-connected system of direct drive wind farm are shown in Table 1.

(114) TABLE-US-00001 TABLE 1 Parameters of grid-connected system of direct drive wind farm parameter symbol value rated power P.sub.n/MW 1.5 current loop control parameters k.sub.P1 0.495e3 k.sub.I1 0.121 PLL control parameters K.sub.P 0.67 K.sub.I 38.2 grid-connected inductance L.sub.x1, L.sub.x2/H 0.3e3, 0.2e3 DC voltage outer loop control k.sub.Pu1 0.1 parameters k.sub.Iu1 10 DC voltage U.sub.dc/kV 1.2

(115) The changes of energy flows V.sub.hi1, V.sub.hi2, V.sub.hi3, V.sub.hi4 and V.sub.hi5 of each coupling energy transfer path within a direct drive wind turbine generator are shown in FIGS. 6A-6E, respectively. During the oscillation process, the d and q axis oscillation currents i.sub.fd1 and i.sub.fq1 of PMSG1 gradually increase, and the energy flow of each transfer path also changes. In FIGS. 6A-6E, with the development of oscillation, the energy flow spirals from locus A to locus B. Among them, the value of V.sub.hi1 gradually increases from 12.22 to 14.24, the value of V.sub.hi2 gradually decreases from 0.012 to 0.016, the value of V.sub.hi3 gradually increases from 482.15 to 570.06, and the value of V.sub.hi4 gradually decreases from 312.56 to 329.76. The value of V.sub.hi5 gradually decreases from 0.57 to 0.47. It can be seen that the energy flow V.sub.hi3 of the transfer path has the largest variation, so the coupling energy transfer path under the action of current loop d axis and PLL is the key energy path of the generator-internal coupling path, and the existing stability evaluation method verifies that it is the key factor affecting the stability of the system, which proves the correctness of the determined key energy path.

(116) The changes of dynamic coupling energy flows V.sub.ij1, V.sub.ij2, V.sub.ij3, V.sub.ij4, V.sub.ij5, V.sub.ij6, V.sub.ij7 and V.sub.ij8 in each transfer path between PMSGi and PMSGj generators are shown in FIGS. 7A-7H respectively. As the d and q axis oscillation current i.sub.fd1 and i.sub.fq1 of PMSG1 increase gradually during the oscillation process, the dynamic coupling energy flow of each transfer path also changes. In FIGS. 7A-7H, with the development of oscillation, the energy flow spirals from locus A to locus B. Among them, the value of V.sub.ij1 gradually increases from 104.11 to 129.58, the value of V.sub.ij2 gradually increases from 1.49 to 1.85, the value of V.sub.ij3 gradually increases from 108.01 to 133.96, and the value of V.sub.ij4 gradually increases from 1.67 to 1.92. The value of V.sub.ij5 gradually decreases from 87.12 to 101.18, the value of V.sub.ij6 gradually decreases from 0.035 to 0.02, the value of V.sub.ij7 gradually increases from 0.30 to 0.31, and the value of V.sub.ij8 gradually decreases from 168.06 to 177.73. Compared with the above, it can be seen that the energy flow V.sub.ij3 of the transfer path has the largest variation. Therefore, for the generator-generator coupling, the coupling transfer path between the current loop d axis of one wind turbine generator and the phase-locked loop of another wind turbine generator is the key energy transfer path of the generator-generator coupling path, and the existing stability evaluation method has verified that it is the key factor affecting the stability of the system. The correctness of the determined key energy path is proved.

(117) The changes of dynamic coupling energy flows V.sub.gi1, V.sub.gi2, V.sub.gi3, V.sub.gi4 and V.sub.gi5 in each transfer path of direct drive wind turbine generators are shown in FIGS. 8A-8E, respectively. As the d and q axis oscillation current i.sub.fd1 and i.sub.fq1 of PMSG1 increase gradually during the oscillation process, the dynamic coupling energy flow of each transfer path also changes. In FIGS. 8A-8E, with the development of oscillation, the energy flow spirals from locus A to locus B. Among them, the value of V.sub.gi1 gradually decreases from 34.33 to 40.38, that of V.sub.gi2 gradually decreases from 0.28 to 0.30, that of V.sub.gi3 gradually increases from 3.59 to 4.15, and that of V.sub.gi4 gradually increases from 10.35 to 11.21. The value of V.sub.gi5 gradually increased from 41.16 to 48.53. It can be seen that the energy flow V.sub.gi5 of the transfer path has the largest variation. Therefore, for the coupling of generator and grid, the coupling transfer path between PLL and the q axis of the power grid is the key energy transfer path of the coupling path of generator and grid, and it is verified by the existing stability evaluation method that it is the key factor affecting the stability of the system, which proves the correctness of the determined key energy path.

(118) It is understood by those skilled in the art that all or part of the process to implement the above exemplary methods may be accomplished by instructing the relevant hardware by means of a computer program that may be stored in a computer readable storage medium. The computer readable storage medium is a disk, optical disc, read-only storage memory or random storage memory, etc.

(119) The above is only the best concrete embodiment of the invention, but the scope of protection of the invention is not limited to this, and any change or replacement that can be easily thought of by a technical person familiar with the technical field within the technical scope disclosed by the invention shall be covered by the scope of protection of the invention.