Method and apparatus of frequency regulation of power system involving renewable energy power generation, device, and storage medium
A method and an apparatus of frequency regulation of a power system involving renewable energy power generation, a computer device, and a non-transitory computer readable storage medium are provided. The method includes: constructing a system frequency dynamic model according to parameters associated with power generator sets in the power system, where the power generator sets comprise a renewable energy power generator set and a conventional energy power generator set; calculating secure operation indexes of the power system according to the system frequency dynamic model of the power system; and obtaining system comprehensive cost indexes of the power system, constructing a reserve allocation model of the power generator sets according to the system comprehensive cost indexes and the secure operation indexes of the power system, and regulating a system frequency of the power system according to the reserve allocation model.
1 . A method of frequency regulation of a power system involving renewable energy power generation, comprising:
constructing a system frequency dynamic model according to parameters associated with power generator sets in the power system, comprising:
obtaining a largest imbalanced power, parameters of the conventional energy power generator set, and parameters of the renewable energy power generator set;
obtaining an equivalent inertia time constant and an equivalent damping coefficient in a preset time period according to the parameters of the conventional energy power generator set and the parameters of the renewable energy power generator set; and
constructing the system frequency dynamic model according to the equivalent inertia time constant, the equivalent damping coefficient, the largest imbalanced power, the parameters of the conventional energy power generator set, and the parameters of the renewable energy power generator set; the power generator sets comprising a renewable energy power generator set and a conventional energy power generator set;
calculating secure operation indexes of the power system according to the system frequency dynamic model of the power system, comprising:
calculating an absolute value of a maximum Rate-of-Change-of-Frequency (RoCoF) of the power system according to a post-fault instantaneous power change amount of the renewable energy power generator set in the preset time period, the equivalent inertia time constant and the largest imbalanced power in the preset time period;
calculating a steady-state power deviation of the conventional energy power generator set, a steady-state power deviation of the renewable energy power generator set, and an absolute value of a steady-state frequency deviation of the power system, according to the equivalent damping coefficient, the largest imbalanced power, the parameters of the conventional energy power generator set, and the parameters of the renewable energy power generator set in the preset time period; and
calculating a maximum frequency deviation of the power system, according to the parameters of the conventional energy power generator set, the parameters of the renewable energy power generator set, and the largest imbalanced power combining with a preset piecewise linear function; and
obtaining system comprehensive cost indexes of the power system, constructing a reserve allocation model of the power generator sets according to the system comprehensive cost indexes and the secure operation indexes of the power system, and regulating a system frequency of the power system according to the reserve allocation model.
2 . The method of claim 1 , wherein the calculating the maximum frequency deviation of the power system, according to the parameters of the conventional energy power generator set, the parameters of the renewable energy power generator set, and the largest imbalanced power combining with the preset piecewise linear function, comprises:
determining a space division of a definition domain of the preset piecewise linear function, generating data samples of the preset piecewise linear function, determining parameters values of the preset piecewise linear function based on the space division of the definition domain of the preset piecewise linear function and the data samples, and constructing the preset piecewise linear function;
constructing linear constraint conditions of the preset piecewise linear function; and
calculating the maximum frequency deviation of the power system according to the preset piecewise linear function and the linear constraint conditions.
3 . The method of claim 1 , wherein the obtaining the system comprehensive cost indexes of the power system, constructing the reserve allocation model of the power generator sets according to the system comprehensive cost indexes and the secure operation indexes of the power system, and regulating the system frequency of the power system according to the reserve allocation model, comprises:
constructing the system comprehensive cost indexes based on the parameters of the conventional energy power generator set and the parameters of the renewable energy power generator set, and constructing an optimization objective function based on the system comprehensive cost indexes;
constructing constraint conditions of the secure operation indexes of the power system;
constructing the reserve allocation model of the power generator sets according to the optimization objective function and the constraint conditions; and
calculating an optimal solution of the reserve allocation model of the power generator sets, and adjusting a reserve capacity of the renewable energy power generator set and a reserve capacity of the conventional energy power generator set based on the optimal solution to regulate the system frequency of the power system.
4 . The method of claim 3 , wherein the constructing the constraint conditions of the secure operation indexes of the power system, comprises:
constructing combination constraint conditions and operation constraint conditions of the conventional energy power generator set, and operation constraint conditions of the renewable energy power generator set;
constructing a power balance constraint condition of the power system and a constraint condition of the reserve capacity of the power system after a tertiary frequency regulation;
constructing a line power flow constraint condition of the power system in a normal operation condition and a line power flow constraint condition of the power system after the primary frequency regulation; and
constructing constraint conditions of a secondary frequency regulation of the power system, and constraint conditions of a frequency security of the power system in a dynamic of the primary frequency regulation.
5 . The method of claim 1 , wherein the system frequency dynamic model after an accident in the time period k comprises:
2
H
sys
(
k
)
Δ
f
(
t
)
=
-
D
sys
(
k
)
Δ
f
(
t
)
+
∑
i
∈
N
G
Δ
P
i
g
e
n
(
t
)
+
∑
j
∈
N
W
Δ
P
j
w
i
n
d
(
t
)
-
P
l
o
s
s
(
k
)
H
sys
(
k
)
=
∑
i
∈
N
G
v
i
,
k
g
e
n
H
i
g
e
n
D
sys
(
k
)
=
∑
i
∈
N
G
v
i
,
k
g
e
n
d
i
g
e
n
{
τ
i
gen
Δ
P
i
gen
(
t
)
=
-
Δ
P
i
gen
(
t
)
-
α
i
gen
Δ
f
(
t
)
,
if
v
i
,
k
gen
=
1
,
∀
i
∈
N
G
Δ
P
i
gen
(
t
)
=
0
,
if
v
i
,
k
gen
=
0
,
∀
i
∈
N
G
❘
"\[LeftBracketingBar]"
Δ
P
i
g
e
n
(
t
)
❘
"\[RightBracketingBar]"
≤
PR
i
,
k
g
e
n
,
∀
i
∈
N
G
{
Δ
P
j
w
i
n
d
(
t
)
=
-
k
j
inertia
Δ
f
(
t
)
-
k
j
d
r
o
o
p
Δ
f
(
t
)
,
if
v
j
,
k
wind
=
1
,
∀
j
∈
N
W
Δ
P
j
w
i
n
d
(
t
)
=
0
,
if
v
j
,
k
wind
=
0
,
∀
j
∈
N
W
|
Δ
P
j
w
i
n
d
(
t
)
|
≤
P
R
j
,
k
w
i
n
d
,
∀
j
∈
N
W
wherein the parameters of the conventional energy power generator set comprises:
H
i
g
e
n
representing an inertia time constant,
d
i
g
e
n
represent a damping coefficient,
τ
i
g
e
n
representing a time constant,
α
i
gen
representing a speed governor coefficient, N G representing a group of conventional energy power generator sets, and n G representing a number of conventional energy power generator sets in the group;
wherein the parameters of the renewable energy power generator set comprise:
k
j
intertia
representing a virtual inertia time constant of the renewable energy power j,
k
j
droop
representing a droop control coefficient of the renewable energy power generator set j, N W representing a group of renewable energy power generator sets, n W representing a number of renewable energy power generator sets in the group;
H
sys
(
k
)
represents an equivalent inertia time constant,
D
sys
(
k
)
represents an equivalent damping coefficient, Δf(t) represents a frequency deviation of a frequency of a center of inertia of the power system,
Δ
P
i
gen
(
t
)
represents a power regulation amount of the conventional energy power generator set i,
Δ
P
j
wind
(
t
)
represents a power regulation amount of the renewable energy power generator set j,
P
loss
(
k
)
represent a largest imbalanced power, a Boolean variable
v
i
,
k
gen
represents whether the conventional energy power generator set i participates in the primary frequency regulation or not in the time period k, a Boolean variable
v
j
,
k
wind
represents whether the renewable energy power generator set j participates in the primary frequency regulation or not in the time period k,
PR
i
,
k
gen
represents a primary frequency-regulation reserve capacity of the conventional energy power generator set i in the time period k, and
PR
j
,
k
wind
represents a primary frequency-regulation reserve capacity of the renewable energy power generator set j in the time period k.
6 . The method of claim 1 , wherein the calculating the absolute value of the maximum RoCoF of the power system
RoCoF
max
(
k
)
comprise:
-
2
H
sys
(
k
)
RoCoF
max
(
k
)
=
-
P
loss
(
k
)
+
∑
j
Δ
P
wini
_
j
(
k
)
Δ
P
wini
_
j
(
k
)
=
v
j
,
k
wind
min
{
k
j
intertia
RoCoF
max
(
k
)
,
PR
j
,
k
wind
}
,
∀
j
wherein
Δ
P
wini
_
j
(
k
)
represents a post-fault instantaneous power change amount of the renewable energy power generator set in the preset time period k, an instantaneous power support is achieved by a virtual inertia control,
P
l
o
s
s
(
k
)
represents the largest imbalanced power,
H
s
y
s
(
k
)
represents the equivalent inertia time constant, a Boolean variable
v
j
,
k
w
i
n
d
represents whether the renewable energy power generator set j participates in the primary frequency regulation or not in the time period k,
k
j
inertia
represents a virtual inertia time constant of the renewable energy power generator set j, and
P
R
j
,
k
w
i
n
d
represent a primary frequency-regulation reserve capacity of the renewable energy power generator set j in the time period k.
7 . The method of claim 1 , wherein equations for calculating the steady-state power deviation of the conventional energy power generator set, the steady-state power deviation of the renewable energy power generator set, and the absolute value of the steady-state frequency deviation of the power system comprise:
D
sys
(
k
)
Δ
f
s
s
(
k
)
=
-
P
l
o
s
s
(
k
)
+
∑
i
Δ
P
gss
_
i
(
k
)
+
∑
j
Δ
P
wss
_
j
(
k
)
Δ
P
gss
_
i
(
k
)
=
v
i
,
k
g
e
n
min
{
-
α
i
g
e
n
Δ
f
s
s
(
k
)
,
P
R
i
,
k
g
e
n
}
,
∀
i
Δ
P
wss
_
j
(
k
)
=
v
j
,
k
w
i
n
d
min
{
-
k
j
d
r
o
o
p
Δ
f
s
s
(
k
)
,
P
R
j
,
k
w
i
n
d
}
,
∀
j
wherein
Δ
P
gss
_
i
(
k
)
represents the steady-state power deviation of the conventional energy power generator set,
Δ
P
wss
_
j
(
k
)
represents the steady-state power deviation of the renewable energy power generator set,
Δ
f
s
s
(
k
)
represents the absolute value of the steady-state frequency deviation of the power system,
D
s
y
s
(
k
)
represents an equivalent damping coefficient of the preset time period k,
P
l
o
s
s
(
k
)
represents the largest imbalanced power,
α
i
g
e
n
represents a speed governor coefficient,
k
j
droop
represents a droop control coefficient of the renewable energy power generator set j,
P
R
i
,
k
g
e
n
represents a primary frequency-regulation reserve capacity of the conventional energy power generator set in the time period k, and
P
R
j
,
k
w
i
n
d
represents a primary frequency-regulation reserve capacity of the renewable energy power generator set j in the time period k, a Boolean variable
v
i
,
k
gen
represents whether the conventional energy power generator set i participates in the primary frequency regulation or not in the time period k, and a Boolean variable
v
j
,
k
w
i
n
d
represents whether the renewable energy power generator set j participates in the primary frequency regulation or not in the time period k.
8 . The method of claim 3 , wherein the optimization objective function comprises:
min
∑
k
(
∑
i
(
C
i
fixed
u
i
,
k
+
C
i
S
U
z
i
,
k
S
U
+
C
i
s
D
z
i
,
k
S
D
+
C
i
i
n
c
r
P
i
,
k
g
e
n
)
+
∑
j
i
C
j
p
e
n
(
P
j
,
k
m
p
p
t
-
P
j
,
k
w
i
n
d
-
P
R
j
,
k
w
i
n
d
-
S
R
j
,
k
w
i
n
d
-
T
R
j
,
k
w
i
n
d
)
)
over
{
u
i
,
k
,
z
i
,
k
S
U
,
z
i
,
k
S
D
,
v
i
,
k
gen
,
v
j
,
k
wind
,
P
i
,
k
gen
,
P
j
,
k
wind
,
z
i
,
k
gen
,
z
j
,
k
wind
PR
i
,
k
gen
,
PR
j
,
k
wind
,
SR
i
,
k
gen
,
SR
j
,
k
wind
,
TR
i
,
k
gen
,
TR
j
,
k
wind
}
wherein the parameters of the conventional energy power generator set comprise:
C
i
f
i
x
e
d
representing a fixed cost coefficient of power generation,
C
i
S
U
representing a set start-up cost coefficient,
C
i
S
D
representing a set shutdown cost coefficient,
C
i
incr
representing a variable cost coefficient of power generation, u i,k representing an on/off state of the conventional energy power generator set i in the time period k, and a Boolean variable
v
i
,
k
g
e
n
representing whether the conventional energy power generator set i participates in the primary frequency regulation or not in the time period k; and
wherein the parameters of the renewable energy power generator set comprise:
C
j
p
e
n
representing a wind curtailment penalty coefficient of the renewable energy power generator set j,
P
j
,
k
mppt
representing a predicted value of a maximum power point tracking of the renewable energy power generator set j in the time period k, a Boolean variable
v
j
,
k
w
i
n
d
renewable energy power generator set j participates in the primary frequency regulation or not in the time period k;
the system comprehensive cost indexes comprise: a decision variable u i,k representing the on/off state of the conventional energy power generator set i in the time k,
z
i
,
k
S
U
and
z
i
,
k
S
D
representing startup and shutdown actions of the conventional energy power generator set i in the time period k, respectively,
P
i
,
k
g
e
n
representing a planned output of the conventional energy power generator set i in the time period k,
P
j
,
k
w
i
n
d
representing an actual output of the renewable energy power generator set j in the time period k,
PR
i
,
k
gen
,
SR
i
,
k
gen
and
TR
i
,
k
gen
representing a primary frequency-regulation reserve capacity, a secondary frequency-regulation reserve capacity, and a tertiary frequency-regulation reserve capacity of the conventional energy power generator set i in the time period k, respectively,
PR
j
,
k
wind
,
SR
j
,
k
wind
and
TR
j
,
k
wind
representing a primary frequency-regulation reserve capacity, a secondary frequency-regulation reserve capacity, and a tertiary frequency-regulation reserve capacity of the renewable energy power generator set j in the time period k, respectively,
z
i
,
k
gen
representing a post-fault secondary frequency-regulation reserve deployment of the conventional energy power generator set i in the time period k, and
z
i
,
k
wind
representing a post-fault secondary frequency-regulation reserve deployment of the renewable energy power generator set j in the time period k.
9 . The method of claim 4 , wherein
the power balance constraint condition of the power system is:
∑
i
P
i
,
k
gen
+
∑
j
P
j
,
k
wind
=
∑
d
P
d
,
k
load
∀
k
,
wherein
P
i
,
k
gen
represents a planned output of the conventional energy power generator set i in the time period k,
P
j
,
k
wind
represents an actual output of the renewable energy power generator set j in the time period k,
P
d
,
k
load
represents a predicted value of a load d in the time period k; and
the constraint condition of the reserve capacity of the power system after a tertiary frequency regulation is:
∑
i
TR
i
,
k
gen
+
∑
j
TR
j
,
k
wind
=
5
%
∑
d
P
d
,
k
load
∀
k
,
wherein
TR
i
,
k
gen
represents a tertiary frequency-regulation reserve capacity of the conventional energy power generator set i in the time period k,
TR
j
,
k
wind
represents a tertiary frequency-regulation reserve capacity of the renewable energy power generator set j in the time period k.
10 . A computer device comprising a memory and a processor, wherein a computer program is stored in the memory, and the processor, when executing the computer program, performs steps of the method of claim 1 .
11 . A non-transitory computer readable storage medium, having a computer program stored thereon, wherein the computer program, when executed by a processor, causes the processor to perform steps of the method of claim 1 .
12 . A method of frequency regulation of a power system involving renewable energy power generation, comprising:
constructing a system frequency dynamic model according to parameters associated with power generator sets in the power system, the power generator sets comprising a renewable energy power generator set and a conventional energy power generator set;
calculating secure operation indexes of the power system according to the system frequency dynamic model of the power system; and
obtaining system comprehensive cost indexes of the power system, constructing a reserve allocation model of the power generator sets according to the system comprehensive cost indexes and the secure operation indexes of the power system, and regulating a system frequency of the power system according to the reserve allocation model, comprise:
constructing the system comprehensive cost indexes based on the parameters of the conventional energy power generator set and the parameters of the renewable energy power generator set, and constructing an optimization objective function based on the system comprehensive cost indexes;
constructing constraint conditions of the secure operation indexes of the power system;
constructing the reserve allocation model of the power generator sets according to the optimization objective function and the constraint conditions; and
calculating an optimal solution of the reserve allocation model of the power generator sets, and adjusting a reserve capacity of the renewable energy power generator set and a reserve capacity of the conventional energy power generator set based on the optimal solution to regulate the system frequency of the power system.
13 . The method of claim 12 , wherein the constructing the constraint conditions of the secure operation indexes of the power system, comprises:
constructing combination constraint conditions and operation constraint conditions of the conventional energy power generator set, and operation constraint conditions of the renewable energy power generator set;
constructing a power balance constraint condition of the power system and a constraint condition of the reserve capacity of the power system after a tertiary frequency regulation;
constructing a line power flow constraint condition of the power system in a normal operation condition and a line power flow constraint condition of the power system after the primary frequency regulation; and
constructing constraint conditions of a secondary frequency regulation of the power system, and constraint conditions of a frequency security of the power system in a dynamic of the primary frequency regulation.
14 . The method of claim 12 , wherein the optimization objective function comprises:
min
∑
k
(
∑
i
(
C
i
fixed
u
i
,
k
+
C
i
SU
z
i
,
k
SU
+
C
i
SD
z
i
,
k
SD
+
C
i
,
k
incr
P
i
,
k
gen
)
+
∑
j
C
j
pen
(
P
j
,
k
mrrp
-
P
j
,
k
wind
-
PR
j
,
k
wind
-
SR
j
,
k
wind
-
TR
j
,
k
wind
)
over
{
u
i
,
k
,
z
i
,
k
SU
,
z
i
,
k
SD
,
v
i
,
k
gen
,
v
j
,
k
wind
,
P
i
,
k
gen
,
P
j
,
k
wind
,
z
i
,
k
gen
,
z
j
,
k
wind
PR
i
,
k
gen
,
PR
j
,
k
wind
,
SR
i
,
k
gen
,
SR
j
,
k
wind
,
TR
i
,
k
gen
,
TR
j
,
k
wind
}
wherein the parameters of the conventional energy power generator set comprise:
C
i
fixed
representing a fixed cost coefficient of power generation,
C
i
SU
representing a set start-up cost coefficient,
C
i
SD
representing a set shutdown cost coefficient,
C
i
incr
representing a variable cost coefficient of power generation, u i,k representing an on/off state of the conventional energy power generator set i in the time period k, and a Boolean variable
v
i
,
k
gen
representing whether the conventional energy power generator set i participates in the primary frequency regulation or not in the time period k; and
wherein the parameters of the renewable energy power generator set comprise:
C
j
p
e
n
representing a wind curtailment penalty coefficient of the renewable energy power generator set j,
P
j
,
k
mppt
representing a predicted value of a maximum power point tracking of the renewable energy power generator set j in the time period k, a Boolean variable
v
j
,
k
w
i
n
d
representing whether the renewable energy power generator set j participates in the primary frequency regulation or not in the time period k;
the system comprehensive cost indexes comprise: a decision variable u i,k representing the on/off state of the conventional energy power generator set i in the time period k,
z
i
,
k
S
U
and
z
i
,
k
S
D
representing startup and shutdown actions of the conventional energy power generator set i in the time period k, respectively,
P
i
,
k
g
e
n
representing a planned output of the conventional energy power generator set i in the time period k,
P
j
,
k
w
i
n
d
representing an actual output of the renewable energy power generator set j in the time period k,
P
R
i
,
k
g
e
n
,
S
R
i
,
k
g
e
n
and
TR
i
,
k
g
e
n
representing a primary frequency-regulation reserve capacity, a secondary frequency-regulation reserve capacity, and a tertiary frequency-regulation reserve capacity of the conventional energy power generator set i in the time period k, respectively,
P
R
j
,
k
w
i
n
d
,
S
R
j
,
k
w
i
n
d
and
TR
j
k
w
i
n
d
representing a primary frequency-regulation reserve capacity, a secondary frequency-regulation reserve capacity, and a tertiary frequency-regulation reserve capacity of the renewable energy power generator set j in the time period k, respectively,
z
i
,
k
gen
representing a post-fault secondary frequency-regulation reserve deployment of the conventional energy power generator set i in the time period k, and
z
j
,
k
w
i
n
d
representing a post-fault secondary frequency-regulation reserve deployment of the renewable energy power generator set j in the time period k.
15 . The method of claim 13 , wherein
the power balance constraint condition of the power system is:
∑
i
P
i
,
k
g
e
n
+
∑
j
P
j
,
k
w
i
n
d
=
∑
d
P
d
,
k
l
o
a
d
∀
k
,
wherein
p
i
,
k
g
e
n
represents a planned output of the conventional energy power generator set i in the time period k,
P
j
,
k
w
i
n
d
represent an actual output of the renewable energy power generator set j in the time period k,
P
d
,
k
l
o
a
d
represent a predicted value of a load d in the time period k; and
the constraint condition of the reserve capacity of the power system after a tertiary frequency regulation is:
∑
i
T
R
i
,
k
g
e
n
+
∑
j
T
R
j
,
k
w
i
n
d
=
5
%
∑
d
P
d
,
k
l
o
a
d
∀
k
,
wherein
T
R
i
,
k
g
e
n
represents a tertiary frequency-regulation reserve capacity of the conventional energy power generator set i in the time period k,
T
R
j
,
k
w
i
n
d
represent a tertiary frequency-regulation reserve capacity of the renewable energy power generator set j in the time period k.