上图所示为微电网结构,由可控分布式电源、可再生分布式电源、储能及负荷集成而成。
1.可控式分布电源
微电网中的可控分布式电源主要由微型燃气轮机、燃料电池等,其发电成本可以用线性函数表示。其成本函数可以表示为:
C
G
(
t
)
=
[
a
P
G
(
t
)
+
b
]
Δ
t
C_G(t) = [aP_G(t)+b]\Delta t
CG(t)=[aPG(t)+b]Δt
- C G ( t ) C_G(t) CG(t): 表示微型燃气轮机在时段 t t t的发电成本
- a , b a,b a,b:是成本系数
- P G ( t ) P_G(t) PG(t):时段 t t t微型燃气轮机的输出功率
- Δ t \Delta t Δt : 调度步长
由于微型燃气轮机的功率响应速度较快,因此不考虑其爬坡率约束,仅考虑输出功率约束:
P
G
min
≤
P
G
(
t
)
≤
P
G
max
P_G^{\min} \leq P_G(t) \leq P_G^{\max}
PGmin≤PG(t)≤PGmax
- P G min P_G^{\min} PGmin:微型燃气轮机的最小输出功率
- P G max P_G^{\max} PGmax:微型燃气轮机的最大输出功率
2.储能
储能的运行成本主要考虑其一次投资成本和运维成本,在投资回收期内
t
t
t时段的平均充放电成本
C
S
(
t
)
C_S(t)
CS(t)可以表示为:
C
S
(
t
)
=
K
S
[
P
S
d
i
s
(
t
)
/
η
+
P
S
c
h
(
t
)
η
]
Δ
t
C_S(t) = K_S[P_S^{dis}(t)/\eta + P_S^{ch}(t)\eta]\Delta t
CS(t)=KS[PSdis(t)/η+PSch(t)η]Δt
- K S K_S KS:折算后的单位充放电成本
- P S c h ( t ) P_S^{ch}(t) PSch(t):表示在 t t t时段的充电功率
- P S d i s ( t ) P_S^{dis}(t) PSdis(t):表示在 t t t时段的放电功率
- η \eta η:储能单元的充放电效率
储能单元在运行过程中需满足的约束包括:
- 储能的充电和放电功率约束,
P
S
m
a
x
P_S^{max}
PSmax储能允许的最大充放电功率,
U
S
(
t
)
U_S(t)
US(t)表示储能的充放电状态,1:表示放电,0:表示充电。
0 ≤ P S d i s ( t ) ≤ U S ( t ) P S m a x 0 ≤ P S c h ( t ) ≤ [ 1 − U S ( t ) ] P S m a x 0\leq P_S^{dis}(t)\leq U_S(t)P_S^{max}\\ 0\leq P_S^{ch}(t)\leq [1-U_S(t)]P_S^{max} 0≤PSdis(t)≤US(t)PSmax0≤PSch(t)≤[1−US(t)]PSmax - 储能在调度的始末时刻容量相等,有利于储能的循环调度,
N
T
N_T
NT为调度周期,可以取值为24h
η ∑ t = 1 N T [ P S c h ( t ) Δ t ] − 1 η ∑ t = 1 N T [ P S d i s ( t ) Δ t ] = 0 \eta \sum_{t=1}^{N_T}[P_S^{ch}(t)\Delta t] - \frac{1}{\eta} \sum_{t=1}^{N_T}[P_S^{dis}(t)\Delta t] = 0 ηt=1∑NT[PSch(t)Δt]−η1t=1∑NT[PSdis(t)Δt]=0 - 储能各时段的剩余容量约束,
E
S
(
0
)
E_S(0)
ES(0)为储能在调度初始时刻的容量,
E
S
min
,
E
S
max
E_S{\min},E_S{\max}
ESmin,ESmax为储能在调度过程中允许的最小/最大剩余容量,该约束的主要目的在于防止储能过充/过放,延长储能的使用寿命。
E S min ≤ E S ( 0 ) + η ∑ t ′ = 1 t [ P S c h ( t ′ ) Δ t ] − 1 η ∑ t ′ = 1 t [ P S d i s ( t ′ ) Δ t ] ≤ E S m a x E_S{\min} \leq E_S(0) + \eta \sum_{t'=1}^{t} [P_S^{ch}(t')\Delta t]-\frac{1}{\eta} \sum_{t'=1}^{t}[P_S^{dis}(t')\Delta t] \leq E_S^{max} ESmin≤ES(0)+ηt′=1∑t[PSch(t′)Δt]−η1t′=1∑t[PSdis(t′)Δt]≤ESmax
3.需求响应负荷
考虑微电网内包含可平移负荷的情况,其在提供需求响应服务过程中的用电特性可用以下约束表示:
∑
t
=
1
N
T
P
D
R
(
t
)
Δ
t
=
D
D
R
D
D
R
m
i
n
(
t
)
≤
P
D
R
(
t
)
Δ
t
≤
D
D
R
m
a
x
(
t
)
\sum_{t=1}^{N_T} P_{DR}(t)\Delta t = D_{DR}\\ D_{DR}^{min}(t)\leq P_{DR}(t)\Delta t \leq D_{DR}^{max}(t)
t=1∑NTPDR(t)Δt=DDRDDRmin(t)≤PDR(t)Δt≤DDRmax(t)
- P D R ( t ) P_{DR}(t) PDR(t)为 t t t时段微电网对需求响应负荷的实际调度功率
- D D R D_{DR} DDR为需求响应负荷在调度周期内的总用电需求
- D D R m i n ( t ) , D D R m a x ( t ) D_{DR}^{min}(t),D_{DR}^{max}(t) DDRmin(t),DDRmax(t)为需求响应负荷在 t t t时段的最小/最大用电需求,与用户对适应度的要求有关。
用电计划的改变会影响用户舒适,微电网需要给予适度的补偿,
t
t
t时段所付出的调度成本
C
D
R
(
t
)
C_{DR}(t)
CDR(t)可表示为
C
D
R
(
t
)
=
K
D
R
∣
P
D
R
(
t
)
−
P
D
R
∗
(
t
)
∣
Δ
t
C_{DR}(t) = K_{DR}|P_{DR}(t) - P_{DR}^*(t)|\Delta t
CDR(t)=KDR∣PDR(t)−PDR∗(t)∣Δt
- $ K_{DR}$:需求响应负荷的单位调度成本
- $ P_{DR}^*(t) : : :t$时段需求响应负荷的期望用电功率
将上式线性化:
P
D
R
1
(
t
)
=
max
{
0
,
P
D
R
(
t
)
−
P
D
R
∗
(
t
)
}
P
D
R
2
(
t
)
=
max
{
0
,
P
D
R
∗
(
t
)
−
P
D
R
(
t
)
}
P_{DR1}(t) = \max\{0,P_{DR}(t) - P_{DR}^*(t)\}\\ P_{DR2}(t) = \max\{0,P_{DR}^*(t) - P_{DR}(t) \}
PDR1(t)=max{0,PDR(t)−PDR∗(t)}PDR2(t)=max{0,PDR∗(t)−PDR(t)}
则上面绝对值式子可以线性化为:
C
D
R
(
t
)
=
K
D
R
(
P
D
R
1
(
t
)
+
P
D
R
2
(
t
)
)
Δ
t
P
D
R
(
t
)
−
P
D
R
∗
(
t
)
=
P
D
R
1
(
t
)
−
P
D
R
2
(
t
)
C_{DR}(t) = K_{DR}(P_{DR1}(t) + P_{DR2}(t))\Delta t \\ P_{DR}(t) - P_{DR}^*(t) = P_{DR1}(t) - P_{DR2}(t)
CDR(t)=KDR(PDR1(t)+PDR2(t))ΔtPDR(t)−PDR∗(t)=PDR1(t)−PDR2(t)
4. 配电网交互功率
当微电网内各发电单元无法满足负荷需求时,需要向配电网购电,反之,微电网可将富余的电能出售给配电网获取收益。微电网和配电网之间的交互功率需要满足平衡约束:
P
M
b
u
y
(
t
)
−
P
M
s
e
l
l
(
t
)
=
P
S
c
h
(
t
)
+
P
D
R
(
t
)
+
P
L
(
t
)
−
P
G
(
t
)
−
P
S
d
i
s
(
t
)
−
P
P
V
(
t
)
P_M^{buy}(t) - P_M^{sell}(t) = P_S^{ch}(t) + P_{DR}(t) + P_L(t) - P_G(t) - P_S^{dis}(t) - P_{PV}(t)
PMbuy(t)−PMsell(t)=PSch(t)+PDR(t)+PL(t)−PG(t)−PSdis(t)−PPV(t)
- P M b u y ( t ) P_M^{buy}(t) PMbuy(t): t t t时刻微电网向电网购买的功率
- P M s e l l ( t ) P_M^{sell}(t) PMsell(t): t t t时刻微电网向电网出售的功率
- P L ( t ) P_L(t) PL(t): t t t时刻微电网的常规负荷
- P P V ( t ) P_{PV}(t) PPV(t): t t t时刻微电网的光伏输出功率
还需满足:
0
≤
P
M
b
u
y
(
t
)
≤
U
M
(
t
)
P
M
m
a
x
0
≤
P
M
s
e
l
l
(
t
)
≤
[
1
−
U
M
(
t
)
]
M
m
a
x
0 \leq P_M^{buy}(t) \leq U_M(t)P_M^{max}\\ 0 \leq P_M^{sell}(t) \leq [1-U_M(t)]_M^{max}
0≤PMbuy(t)≤UM(t)PMmax0≤PMsell(t)≤[1−UM(t)]Mmax
- P M max P_M^{\max} PMmax:微电网和配电网交换功率的最大值,取值需要考虑配电网和微电网连接处变压器的容量及具体政策等因素
- U M ( t ) U_M(t) UM(t):1表示微电网向配电网购电,0表示微电网向配电网售电
在
t
t
t时刻,微电网和配电网的交互成本
C
M
(
t
)
C_M(t)
CM(t)可表示为:
C
M
(
t
)
=
λ
(
t
)
[
P
M
b
u
y
(
t
)
−
P
M
s
e
l
l
(
t
)
]
Δ
t
C_M(t) = \lambda(t)[P_M^{buy}(t) - P_M^{sell}(t)]\Delta t
CM(t)=λ(t)[PMbuy(t)−PMsell(t)]Δt
- λ ( t ) \lambda(t) λ(t):配电网日勤交易电价