Direct-drive wind farm parameter tuning method and system considering the interaction between generators
A direct-drive wind farm parameter tuning method and system considering the interaction between generators are provided. The method includes: collecting initial oscillation current of each direct-drive wind turbine port in the direct-drive wind farm; according to collected data, taking period of dominant oscillation mode as the iteration period to calculate the overall farm-grid interaction dynamic energy and parameter tuning index of the direct-drive wind farm in current iteration period; obtaining stability level of the direct-drive wind farm based on the overall farm-grid interaction dynamic energy in the current iteration period; when the system is unstable, the optimization model of key control parameters is established with the minimum value of the parameter tuning index in the current iteration period as the objective function and the range of each key control parameters as the constraint condition to achieve the key control parameters tuning of the direct-drive wind farm.
1 . A direct-drive wind farm parameter tuning method considering the interaction between generators, comprises the following steps:
collecting initial oscillation current of each direct-drive wind turbine port in the direct-drive wind farm by the phase measurement unit (PMU);
according to the collected data, taking a period of dominant oscillation mode as an iteration period to calculate an overall farm-grid interaction dynamic energy and parameter tuning index of the direct-drive wind farm in current iteration period;
obtaining stability level of the direct-drive wind farm based on the overall farm-grid interaction dynamic energy in the current iteration period; when the system is unstable, an optimization model of key control parameters is established with a minimum value of the parameter tuning index in the current iteration period as an objective function and ranges of each key control parameters as constraint conditions, wherein the key control parameters comprise: a current loop proportional gain, a phase-locked loop proportional gain, a grid-connected line inductance, an active current reference value;
determining optimal key control parameter values based on the optimization model of the key control parameters, and tuning the key control parameters of the direct-drive wind farm according to the optimal key control parameter values.
2 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 1 , wherein, based on the collected data, the overall farm-grid interaction dynamic energy and the parameter tuning index of the direct-drive wind farm in the current iteration period are calculated by taking the period of the dominant oscillation mode as the iteration period, comprising:
taking the collected data as initial iteration period data, generator-grid interaction dynamic energy, inter-generator coupling interaction dynamic energy and inter-generator induction interaction dynamic energy of each Permanent Magnetic Synchronous Generator (PMSG) in each iteration period are calculated, where the period of dominant oscillation mode is the iteration period;
based on the generator-grid interaction dynamic energy, the inter-generator coupling interaction dynamic energy and the inter-generator induction interaction dynamic energy of each PMSG in each iteration period, obtaining the overall farm-grid interaction dynamic energy of the direct-drive wind farm in the current iteration period;
based on the generator-grid interaction dynamic energy, the inter-generator coupling interaction dynamic energy of each PMSG in each iteration period, obtaining the parameter tuning index of the direct-drive wind farm in the current iteration period.
3 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 1 , the key control parameter optimization model is expressed as:
{
min
α
S
n
=
η
F
(
α
)
s
.
t
.
P
(
K
)
=
0
K
P_min
j
≤
K
P
j
≤
K
P_max
j
K
P
θ
_min
j
≤
K
P
θ
j
≤
K
P
θ
_max
j
L
xj_min
≤
L
xj
≤
L
xj_max
i
gdj_min
*
≤
i
gdj
*
≤
i
gdj_max
*
α
=
[
K
P
j
,
K
P
θ
j
,
L
xj
,
i
gdj
*
]
,
where S n is the objective function of the current iteration period, nr is the parameter tuning index of the current iteration period, P(K)=0 means that the power flow meets the static security and stability conditions;
K
P
_
min
j
and
K
P
_
max
j
are the lower and upper limit of the current loop proportional gain
K
P
j
of the jth PMSG respectively;
K
P
θ
_
min
j
and
K
P
θ
_
max
j
are the lower and upper limit of the PLL proportion gain
K
P
θ
j
of the jth PMSG respectively; L xj_min and L xj_max are the lower and upper limit of the integration distance of the jth PMSG respectively;
i
gdj_min
*
and
i
gdj_max
*
are the lower and upper limit of the reference value of active current of the jth PMSG respectively; α is set of decision variables.
4 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 3 , wherein, the parameter tuning index ng of the direct-drive wind farm in the current iteration period is expressed as:
η
f
=
Δ
W
Fs
n
+
Δ
W
Fcp
n
.
where
Δ
W
Fs
n
is the overall generator-grid interaction dynamic energy of the direct-drive wind farm after n iteration periods,
Δ
W
Fcp
n
is the overall inter-generator coupling interaction dynamic energy of the direct-drive wind farm after n iteration periods, and n is the number of iteration periods from the initial iteration period to the current iteration period;
the overall farm-grid interaction dynamic energy
Δ
W
F
n
of the direct-drive wind farm after n iteration periods is expressed as:
Δ
W
F
n
=
Δ
W
Fs
n
+
Δ
W
Fcp
n
+
Δ
W
Fin
n
where
Δ
W
Fin
n
is the overall inter-generator induction interaction dynamic energy of the direct-drive wind farm after n iteration periods.
5 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 4 , wherein, the stability level of the direct-drive wind farm is obtained based on the overall farm-grid interaction dynamic energy
Δ
W
F
n
after n iteration periods, comprising:
when
Δ
W
F
n
<
0
,
the direct-drive wind farm has a positive damping effect on oscillation, and the system is in a stable state; the smaller the value, the higher the stability level;
when
Δ
W
F
n
=
0
,
the direct-drive wind farm has no damping effect on the oscillation, and the system is in a critical stable state;
when
Δ
W
F
n
>
0
,
the direct-drive wind farm has a negative damping effect on the oscillation, and the system oscillation divergence will be completely unstable.
6 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 4 , wherein,
the overall generator-grid interaction dynamic energy
Δ
W
Fs
n
of the direct-drive wind farm after n iteration periods is expressed as:
Δ
W
Fs
n
=
∑
k
=
1
n
∑
j
=
1
m
∑
i
=
1
m
i
≠
j
Δ
W
js
(
k
)
where
Δ
W
js
(
k
)
represents the generator-grid interaction dynamic energy of the jth PMSG in the kth iteration period, and m represents total number of PMSGs;
the overall inter-generator coupling interaction dynamic energy
Δ
W
F
c
p
n
of the direct-drive wind farm after n iteration periods is expressed as:
Δ
W
Fcp
n
=
∑
k
=
1
n
∑
j
=
1
m
∑
i
=
1
m
i
≠
j
Δ
W
ij_cp
(
k
)
where
Δ
W
ij_cp
(
k
)
represents the inter-generator coupling interaction dynamic energy between the jth PMSG and the ith PMSG in the kth iteration period;
the overall inter-generator induction interaction dynamic energy of the direct-drive wind farm after n iteration periods
Δ
W
Fin
n
is expressed as:
Δ
W
Fin
n
=
∑
k
=
1
n
∑
j
=
1
m
∑
i
=
1
m
i
≠
j
Δ
W
ij_in
(
k
)
where
Δ
W
ij_in
(
k
)
represents the inter-generator induction interaction dynamic energy between the jth PMSG and the ith PMSG in the kth iteration period.
7 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 6 , wherein, the generator-grid interaction dynamic energy of the jth PMSG in the kth iteration period
Δ
W
js
(
k
)
is expressed as:
Δ
W
js
(
k
)
=
0.5
ω
(
U
Bd
0
j
+
I
gd
0
j
R
s
)
2
+
[
ω
I
gd
0
j
(
L
s
+
L
sj
)
]
2
∑
r
=
1
4
∑
t
=
1
4
(
-
1
)
r
C
rt
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
jd
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
os
(
φ
gdr
jk
-
φ
gdt
jk
-
φ
d
1
t
j
)
+
❘
"\[LeftBracketingBar]"
M
jq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
j
-
φ
gqt
j
-
φ
q
1
t
j
)
]
Δ
t
-
0.5
ωω
0
I
gd
0
j
(
L
s
+
L
xj
)
{
∑
r
=
1
,
3
∑
t
=
1
4
C
rt
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
[
-
❘
"\[LeftBracketingBar]"
M
jd
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
jk
-
φ
gdt
jk
-
φ
d
2
t
j
)
-
❘
"\[LeftBracketingBar]"
M
jq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
jk
-
φ
gqt
jk
-
φ
q
2
t
j
)
]
+
∑
r
=
2
,
4
∑
t
=
1
4
C
rt
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
jd
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
jk
-
φ
gdt
jk
-
φ
g
2
t
j
)
+
❘
"\[LeftBracketingBar]"
M
jq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
jk
-
φ
gqt
jk
-
φ
q
2
t
j
)
]
}
Δ
t
+
0.5
ω
R
s
2
+
ω
2
(
L
s
+
L
sj
)
2
{
∑
r
=
1
,
3
∑
t
=
1
4
(
-
1
)
t
+
1
C
rt
[
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
jk
-
φ
gqt
jk
-
φ
1
t
j
)
-
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
jk
-
φ
gdt
jk
-
φ
1
t
j
)
]
+
∑
r
=
2
,
4
∑
t
=
1
4
(
-
1
)
t
C
rt
[
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
sin
(
φ
gdr
jk
-
φ
gqt
jk
-
φ
1
t
j
)
-
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
jk
-
φ
gdt
jk
-
φ
1
t
j
)
]
}
Δ
t
where
φ
d
1
t
j
=
{
σ
r
-
σ
t
+
φ
c
j
+
φ
mdj
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
c
j
+
φ
mdj
j
+
90
°
t
=
2
,
4
φ
q
1
t
j
=
{
σ
r
-
σ
t
+
φ
c
j
+
φ
mqj
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
c
j
+
φ
mqj
j
+
90
°
t
=
2
,
4
;
φ
d
2
t
j
=
{
σ
r
-
σ
t
+
φ
mdj
+
90
°
t
=
1
,
3
σ
r
-
σ
t
+
φ
mdj
j
t
=
2
,
4
φ
q
2
t
j
=
{
σ
r
-
σ
t
+
φ
mqj
j
+
90
°
t
=
1
,
3
σ
r
-
σ
t
+
φ
mqj
j
t
=
2
,
4
;
φ
1
t
j
=
{
σ
r
-
σ
t
-
arc
tan
[
R
s
/
(
ω
L
s
+
ω
L
xj
)
]
t
=
1
,
3
σ
r
-
σ
t
+
arc
tan
[
(
ω
L
s
+
ω
L
xj
)
/
R
s
]
t
=
2
,
4
φ
c
j
=
{
arc
tan
[
(
U
Bd
0
j
+
I
gd
0
j
R
s
)
/
ω
I
gd
0
j
(
L
s
+
L
xj
)
]
r
=
1
,
3
-
arc
tan
[
ω
I
gd
0
j
(
L
s
+
L
xj
)
/
(
U
Bd
0
j
+
I
gd
0
j
R
s
)
]
r
=
2
,
4
;
σ
r
=
{
φ
sj
r
=
1
,
2
φ
si
r
=
3
,
4
σ
t
=
{
φ
sj
t
=
1
,
2
φ
si
t
=
3
,
4
;
C
rt
=
{
I
sj
I
sj
r
,
t
∈
[
1
,
2
]
I
si
I
si
r
,
t
∈
[
3
,
4
]
I
sj
I
si
others
;
{
M
jd
j
=
M
pllj
(
ω
0
L
s
+
ω
0
L
xj
)
M
jq
j
=
M
pllj
(
R
s
+
sL
s
+
sL
xj
)
M
pllj
=
sK
P
θ
j
+
K
I
θ
j
s
2
;
(
G
d
1
jK
G
d
2
jK
G
d
3
jK
G
d
4
jK
G
q
1
jK
G
q
2
jK
G
q
3
jK
G
q
4
jK
G
d
1
iK
G
d
2
iK
G
d
3
iK
G
d
4
iK
G
q
1
iK
G
q
2
iK
G
q
3
iK
G
q
4
iK
)
=
(
G
d
1
j
G
d
2
j
G
d
3
j
G
d
4
j
G
q
1
j
G
q
2
j
G
q
3
j
G
q
4
j
G
d
1
i
G
d
2
i
G
d
3
i
G
d
4
i
G
q
1
i
G
q
2
i
G
q
3
i
G
q
4
i
)
k
;
{
G
d
1
j
=
F
d
j
[
U
Bd
0
j
+
I
gd
0
j
(
R
s
+
sL
s
+
sL
xj
)
]
G
d
2
j
=
F
d
j
I
gd
0
j
ω
0
(
-
L
s
-
L
xj
)
G
d
3
j
=
F
d
j
I
gd
0
j
(
R
s
+
sL
s
)
G
d
4
j
=
F
d
j
(
-
I
gd
0
j
ω
0
L
s
)
;
{
G
q
1
j
=
F
q
j
(
R
s
+
sL
s
+
sL
xj
)
G
q
2
j
=
F
q
j
ω
0
(
-
L
s
-
L
xj
)
G
q
3
j
=
F
q
j
(
R
s
+
sL
s
)
G
q
4
j
=
F
q
j
(
-
ω
0
L
s
)
;
F
d
j
=
(
K
P
j
+
K
I
j
s
)
(
K
Pu
j
+
K
Iu
j
s
)
K
P
j
+
K
I
j
s
+
R
1
j
+
sL
1
j
1
sC
dc
U
dc
0
;
F
q
j
=
-
(
K
P
j
+
K
I
j
S
)
ω
0
C
f
K
P
j
+
K
I
j
s
+
R
1
j
+
sL
1
j
;
where ω is sub/super synchronous oscillation frequency in dq axis, ω 0 is angular frequency of grid, U Bd0j and I gd0j are steady-state values of d-axis voltage and current at port of the jth PMSG respectively, R s and L s are resistance and inductance of AC line, L xj is inductance of grid-connecting line,
φ
g
d
r
jk
,
φ
g
d
t
jk
,
φ
g
q
r
jk
and
φ
g
q
t
jk
are phases of
G
d
r
jK
,
G
d
t
jK
,
G
q
r
jK
and
G
q
t
jK
respectively, φ sj and φ si are initial phase of initial sub/super synchronous current of the jth and ith PMSG respectively, I sj and I si are amplitudes of the initial sub/super synchronous current of the jth and ith PMSG respectively, s is Laplace operator
φ
mdj
j
and
φ
mqj
j
are phases of
M
jd
j
and
M
j
q
j
respectively,
K
P
j
and
K
I
j
are current loop proportion and integral gains of the jth PMSG respectively,
K
P
θ
j
and
K
I
θ
j
are proportional and integral gains of phase-locked loop of the jth PMSG respectively,
K
P
u
j
and
K
I
u
j
are proportional and integral gains of the voltage outer loop of the jth PMSG respectively, R 1j and L 1j are equivalent resistance and inductance of filter respectively, C dc is DC capacitance, U dc0 is steady-state value of DC bus voltage, C f is filter capacitor.
8 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 6 , wherein, the inter-generator coupling interaction dynamic energy between the jth PMSG and the ith PMSG in the kth iteration period
Δ
W
ij_cp
(
k
)
is expressed as:
Δ
W
ij
_
cp
(
k
)
=
0.5
ω
(
U
Bd
0
j
+
I
gd
0
j
R
s
)
2
+
(
ω
I
gd
0
j
(
L
s
+
L
xj
)
)
2
∑
r
=
1
4
∑
t
=
1
4
(
-
1
)
r
C
rt
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
id
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
jk
-
φ
gdt
ik
-
φ
d
3
t
j
)
+
❘
"\[LeftBracketingBar]"
M
iq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
jk
-
φ
gqt
ik
-
φ
q
3
t
j
)
]
Δ
t
-
0.5
ωω
0
I
gd
0
j
(
L
s
+
L
xj
)
{
∑
r
=
1
,
3
∑
t
=
1
4
C
rt
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
[
-
❘
"\[LeftBracketingBar]"
M
id
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
jk
-
φ
gdt
ik
-
φ
d
4
t
j
)
-
❘
"\[LeftBracketingBar]"
M
iq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
jk
-
φ
gqt
ik
-
φ
q
4
t
)
]
+
∑
r
=
2
,
4
∑
t
=
1
4
C
rt
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
id
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
iK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
jk
-
φ
gdt
ik
-
φ
d
4
t
j
)
+
❘
"\[LeftBracketingBar]"
M
iq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
iK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
jk
-
φ
gqt
ik
-
φ
q
4
t
j
)
]
}
Δ
t
+
0.5
ω
(
I
gd
0
j
R
s
)
2
+
(
ω
I
gd
0
j
L
s
)
2
∑
r
=
1
4
∑
t
=
1
4
(
-
1
)
r
C
rt
❘
"\[LeftBracketingBar]"
G
dr
iK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
id
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
ik
-
φ
gdt
jk
-
φ
d
5
t
j
)
+
❘
"\[LeftBracketingBar]"
M
iq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
ik
-
φ
gqt
jk
-
φ
q
5
t
j
)
]
Δ
t
-
0.5
ωω
0
I
gd
0
j
L
s
{
∑
r
=
1
,
3
∑
t
=
1
4
C
rt
❘
"\[LeftBracketingBar]"
G
qr
iK
❘
"\[RightBracketingBar]"
[
-
❘
"\[LeftBracketingBar]"
M
jd
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
ik
-
φ
gdt
jk
-
φ
d
2
t
j
)
-
❘
"\[LeftBracketingBar]"
M
jq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
ik
-
φ
gqt
jk
-
φ
q
2
t
j
)
]
+
∑
r
=
2
,
4
∑
t
=
1
4
C
rt
❘
"\[LeftBracketingBar]"
G
qr
iK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
jk
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
jK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
ik
-
φ
gdt
jk
-
φ
d
2
t
j
)
+
❘
"\[LeftBracketingBar]"
M
jq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
jK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
ik
-
φ
gqt
jk
-
φ
q
2
t
j
)
]
}
Δ
t
+
0.5
ω
R
s
2
+
(
ω
L
s
)
2
{
∑
r
=
1
,
3
∑
t
=
1
4
(
-
1
)
t
+
1
C
rt
[
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
jk
-
φ
gqt
ik
-
φ
2
t
j
)
-
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
jk
-
φ
gdt
ik
-
φ
2
t
j
)
]
+
∑
r
=
2
,
4
∑
t
=
1
4
(
-
1
)
t
C
rt
[
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
iK
❘
"\[RightBracketingBar]"
sin
(
φ
gdr
jk
-
φ
gqt
ik
-
φ
2
t
j
)
-
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
iK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
jk
-
φ
gdt
ik
-
φ
2
t
j
)
]
}
Δ
t
+
0.5
ωω
0
L
s
{
∑
r
=
1
,
3
∑
t
=
1
4
C
rt
[
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
iK
❘
"\[RightBracketingBar]"
sin
(
φ
gdr
jk
-
φ
gdt
ik
-
φ
d
6
t
j
)
-
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
iK
❘
"\[RightBracketingBar]"
sin
(
φ
gqr
jk
-
φ
gqt
ik
-
φ
g
6
t
j
)
]
+
∑
r
=
2
,
4
∑
t
=
1
4
C
rt
[
❘
"\[LeftBracketingBar]"
G
dr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
dt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gdr
jk
-
φ
gdt
ik
-
φ
d
6
t
)
-
❘
"\[LeftBracketingBar]"
G
qr
jK
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
qt
iK
❘
"\[RightBracketingBar]"
cos
(
φ
gqr
jk
-
φ
gqt
ik
-
φ
d
6
t
j
)
]
}
Δ
t
where
φ
d
3
t
j
=
{
σ
r
-
σ
t
+
φ
c
j
+
φ
mdi
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
c
j
+
φ
mdi
j
+
90
°
t
=
2
,
4
φ
q
3
t
j
=
{
σ
r
-
σ
t
+
φ
c
j
+
φ
mqi
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
c
j
+
φ
mqi
j
+
90
°
t
=
2
,
4
;
φ
d
4
t
j
=
{
σ
r
-
σ
t
+
φ
mdi
j
+
90
°
t
=
1
,
3
σ
r
-
σ
t
+
φ
mdi
j
t
=
2
,
4
φ
q
4
t
j
=
{
σ
r
-
σ
t
+
φ
mqi
j
+
90
°
t
=
1
,
3
σ
r
-
σ
t
+
φ
mqi
j
t
=
2
,
4
;
φ
d
5
t
j
=
{
σ
r
-
σ
t
+
φ
p
j
+
φ
mdj
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
p
j
+
φ
mdj
j
+
90
°
t
=
2
,
4
φ
d
5
t
j
=
{
σ
r
-
σ
t
+
φ
p
j
+
φ
mqj
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
p
j
+
φ
mqj
j
+
90
°
t
=
2
,
4
;
φ
q
6
t
j
=
{
σ
r
-
σ
t
t
=
1
,
3
σ
r
-
σ
t
+
90
°
t
=
2
,
4
φ
2
t
j
=
{
σ
r
-
σ
t
-
arc
tan
[
R
s
/
(
ω
L
s
)
]
t
=
1
,
3
σ
r
-
σ
t
+
arc
tan
[
(
ω
L
s
)
/
R
s
]
t
=
2
,
4
;
φ
p
j
=
{
arc
tan
[
R
s
/
(
ω
L
s
)
]
r
=
1
,
3
-
arc
tan
[
(
ω
L
s
)
/
R
s
]
r
=
2
,
4
;
{
M
id
j
=
M
pllj
ω
0
L
s
M
iq
j
=
M
pllj
(
R
s
+
sL
s
)
.
9 . The direct-drive wind farm parameter tuning method considering the interaction between generators according to claim 6 , wherein, the inter-generator induction interaction dynamic energy between the jth PMSG and the ith PMSG in the kth iteration period
Δ
W
ij_in
(
k
)
is expressed as:
Δ
W
ij_in
(
k
)
=
0.5
ω
I
gd
0
j
R
s
2
+
(
ω
L
s
)
2
∑
r
=
1
4
∑
t
=
1
4
(
-
1
)
r
C
r
t
❘
"\[LeftBracketingBar]"
G
d
r
iK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
id
j
G
dt
iK
❘
"\[LeftBracketingBar]"
cos
(
φ
gdr
ik
-
φ
gdt
ik
-
φ
d
7
t
j
)
+
❘
"\[LeftBracketingBar]"
M
iq
j
G
q
t
iK
❘
"\[LeftBracketingBar]"
cos
(
φ
gdr
ik
-
φ
gqt
ik
-
φ
q
7
t
j
)
]
Δ
t
-
0.5
ω
ω
0
I
gd
0
j
L
s
{
∑
r
=
1
,
3
∑
t
=
1
4
C
r
t
❘
"\[LeftBracketingBar]"
G
q
r
iK
❘
"\[RightBracketingBar]"
[
-
❘
"\[LeftBracketingBar]"
M
id
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
d
t
iK
❘
"\[LeftBracketingBar]"
cos
(
φ
g
q
r
ik
-
φ
g
d
t
ik
-
φ
d
2
t
j
)
-
❘
"\[LeftBracketingBar]"
M
iq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
q
t
iK
❘
"\[RightBracketingBar]"
cos
(
φ
g
q
r
ik
-
φ
g
q
t
ik
-
φ
q
2
t
j
)
]
+
∑
r
=
2
,
4
∑
t
=
1
4
C
rt
❘
"\[LeftBracketingBar]"
G
q
r
iK
❘
"\[RightBracketingBar]"
[
❘
"\[LeftBracketingBar]"
M
id
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
d
t
iK
❘
"\[LeftBracketingBar]"
sin
(
φ
g
q
r
ik
-
φ
g
d
t
ik
-
φ
d
2
t
j
)
+
❘
"\[LeftBracketingBar]"
M
iq
j
❘
"\[RightBracketingBar]"
❘
"\[LeftBracketingBar]"
G
q
t
iK
❘
"\[RightBracketingBar]"
sin
(
φ
g
q
r
ik
-
φ
g
q
t
ik
-
φ
q
2
t
j
)
]
}
Δ
t
where
φ
d
7
t
j
=
{
σ
r
-
σ
t
+
φ
p
j
+
φ
mdi
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
p
j
+
φ
mdi
j
+
90
°
t
=
2
,
4
φ
q
7
t
j
=
{
σ
r
-
σ
t
+
φ
p
j
+
φ
mqi
j
t
=
1
,
3
σ
r
-
σ
t
+
φ
p
j
+
φ
mqi
j
+
90
°
t
=
2
,
4
.
10 . A direct-drive wind farm parameter tuning system considering the interaction between generators, comprises a data acquisition module, a system stability evaluation module and a parameter optimization module;
the data acquisition module is used to collect the initial oscillation current of each PMSG port in the direct-driven wind farm;
the system stability evaluation module is used to, according to the collected data, take a period of the dominant oscillation mode as an iteration period to calculate an overall farm-grid interaction dynamic energy and parameter tuning index of the direct-drive wind farm in the current iteration period; obtaining stability level of the direct-drive wind farm based on the overall farm-grid interaction dynamic energy in the current iteration period;
the parameter optimization module is used to, according to the result of system stability assessment, when the system is unstable, establish an optimization model of key control parameters with the minimum value of the parameter tuning index in the current iteration period as an objective function and the ranges of each key control parameters as the constraint conditions, wherein the key control parameters comprise: a current loop proportional gain, a phase-locked loop proportional gain, a grid-connected line inductance, an active current reference value;
wherein the parameter optimization module is further used to, determine optimal key control parameter values based on the optimization model of the key control parameters, and tune the key control parameters of the direct-drive wind farm according to the optimal key control parameter values.