日前现货市场安全约束机组组合数学模型
1.目标函数:
日前现货市场竞价交易出清数学模型的目标函数:
m a x ( ∑ i = 1 N ∑ t = 1 T [ C i , t ( P i , t ) + C i , t U ] − ∑ u = 1 M ∑ t = 1 T ( C u , t ( L u , t ) ) + ∑ s = 1 N S ∑ t = 1 T K [ S L s + + S L s − ] ) max(\sum_{i=1}^{N}\sum_{t=1}^{T}[C_{i,t}(P_{i,t})+C_{i,t}^{U}]-\sum_{u=1}^{M}\sum_{t=1}^{T}(C_{u,t}(L_{u,t}))+\sum_{s=1}^{NS}\sum_{t=1}^{T}K[SL_s^+ + SL^-_s]) max(i=1∑Nt=1∑T[Ci,t(Pi,t)+Ci,tU]−u=1∑Mt=1∑T(Cu,t(Lu,t))+s=1∑NSt=1∑TK[SLs++SLs−])
其中:
- T T T 表示所考虑的时段综述,若每15分钟一个时段,考虑96个时段
- N N N 表示机组总台数
- M M M 表示参与日前现货市场报价的市场用户数
- P i , t P_{i,t} Pi,t 表示机组 i i i 在 t t t 时段的出力(输出功率)
- L u , t L_{u,t} Lu,t 表示市场用户 u u u 在 t t t 时段的用电
- C i , t ( P i , t ) C_{i,t}(P_{i,t}) Ci,t(Pi,t) 、 C i , t U C_{i,t}^{U} Ci,tU 分别为机组 i i i 在时段 t t t 的运行费用、启动费用,其中*组运行费用 C i , t ( P i , t ) C_{i,t}(P_{i,t}) Ci,t(Pi,t) *是与机组申报的各段出力区间和对应能量价格有关的多段线性函数;机组启动费用 C i , t U C_{i,t}^{U} Ci,tU 是与机组停机时间有关的函数,表示机组在不同状态(冷态、温态、热态)下的启动费用;
- C u , t ( L u , t ) C_{u,t} (L_{u,t}) Cu,t(Lu,t) 代表市场用户分段报价曲线
- K K K 表示市场出清优化的网络潮流约束松弛罚因子
- S L s + SL_s^+ SLs+ , S L s − SL_s^{-} SLs− 分别为断面 s s s 的正反向潮流松弛变量,
- N S NS NS 为断面总数
2 约束条件
(1)系统负荷平衡约束
对于每个时段 t t t ,负荷平衡约束可以描述为
∑ i = 1 N P i , t + ∑ j = 1 N T T j , t = D t + ∑ u = 1 M L u , t \sum_{i=1}^{N}P_{i,t} + \sum_{j=1}^{NT}T_{j,t} = D_{t} + \sum_{u=1}^{M}L_{u,t} i=1∑NPi,t+j=1∑NTTj,t=Dt+u=1∑MLu,t
- P i , t P_{i,t} Pi,t 表示机组 i i i 在时段 t t t 的出力
- T j , t T_{j,t} Tj,t 表示联络线 j j j 在时段t的计划功率(别的地方送入为正,输出到别的地方为负)
- N T NT NT 为联络线的总数
- D t D_t Dt 为非市场用户在时段 t t t 的用电
- L u , t L_{u,t} Lu,t 表示用户 u u u 在 t t t 时段的用电
(2)系统正备用容量约束
在确保系统平衡的前提下,为了防止系统负荷预测偏差以及各种实际运行事故带来的系统供需不平衡波动,一般系统需要留有一定的容量备用
需要保证每天的总开机容量满足系统的最小备用容量:系统正备用容量约束可以描述为
∑ i = 1 N α i , t P i , t m a x ≥ D t + ∑ u = 1 M L u , t − ∑ j = 1 N T T j , t + R t U + ( 1 − μ t w ) ∑ j = 1 N P j , t w + ( 1 − μ t s ) ∑ k = 1 N P k , t s \sum_{i=1}^N\alpha_{i,t}P_{i,t}^{max} \geq D_{t} + \sum_{u=1}^{M}L_{u,t} - \sum_{j=1}^{NT}T_{j,t} + R_t^U + (1-\mu_t^w)\sum_{j=1}^{N}P_{j,t}^w + (1-\mu_t^s)\sum_{k=1}^{N}P_{k,t}^s i=1∑Nαi,tPi,tmax≥Dt+u=1∑MLu,t−j=1∑NTTj,t+RtU+(1−μtw)j=1∑NPj,tw+(1−μts)k=1∑NPk,ts
- α i , t \alpha_{i,t} αi,t 表示机组 i i i 在时段 t t t 的启停状态,0表示停机,1 表示开机
- P i , t m a x P_{i,t}^{max} Pi,tmax 表示机组 i i i 在时段 t t t 的最大出力
- L u , t L_{u,t} Lu,t 表示用户 u u u 在 t t t 时段的用电
- R t U R_t^U RtU 表示时段 t t t 的系统正备用容量要求
- μ t w \mu_t^w μtw 表示 t t t 时段的风电置信因子
- P j , t w P_{j,t}^w Pj,tw 表示为风电机组 j j j 在时段 t t t 的预测出力
- μ t s \mu_t^s μts 表示 t t t 时段的光伏置信因子
- P k , t w P_{k,t}^w Pk,tw 表示为光伏机组 k k k 在时段 t t t 的预测出力
(3)系统负备用容量约束
∑ i = 1 N α i , t P i , t m i n ≤ D t + ∑ u = 1 M L u , t − ∑ j = 1 N T T j , t − R t D \sum_{i=1}^N\alpha_{i,t}P_{i,t}^{min} \leq D_{t} + \sum_{u=1}^{M}L_{u,t} - \sum_{j=1}^{NT}T_{j,t} - R_t^D i=1∑Nαi,tPi,tmin≤Dt+u=1∑MLu,t−j=1∑NTTj,t−RtD
- P i , t m i n P_{i,t}^{min} Pi,tmin 为机组 i i i 在时段 t t t 的最小出力
- L u , t L_{u,t} Lu,t 用户 u u u 在 t t t 时段的用电
- R t D R_t^D RtD 表示时段 t t t 的系统正备用容量要求
(4)特殊机组状态约束
1,算法和人工判断确定为必开机组的,应处于开机状态:
α
i
,
t
=
1
,
∀
i
∈
I
s
1
\alpha_{i,t} = 1, \quad \forall i \in I_{s1}
αi,t=1,∀i∈Is1
- I s 1 I_{s1} Is1 必开机组集合
2,算法和人工判断确定为必停机组的,应处于关机状态:
α i , t = 0 , ∀ i ∈ I s 2 \alpha_{i,t} = 0, \quad \forall i \in I_{s2} αi,t=0,∀i∈Is2
I s 2 I_{s2} Is2 必停机组集合
(5)机组出力上下限约束
机组的出力应该处于其最大/最小出力范围之内,
α i , t P i , t m i n ≤ P i , t ≤ α i , t P i , t m a x , ∀ i ∈ N , t ∈ T \alpha_{i,t} P_{i,t} ^ {min} \leq P_{i,t} \leq \alpha_{i,t} P_{i,t}^{max}, \qquad \forall i \in N,t\in T αi,tPi,tmin≤Pi,t≤αi,tPi,tmax,∀i∈N,t∈T
对于必开机组,在其必开时段内,要求 α i , t = 1 \alpha_{i,t} = 1 αi,t=1 ,若有最低出力要求,则在上式中 P i , t P_{i,t} Pi,t 取为对应时段的必开最低出力
(6)机组群出力上限约束
机组群的出力应该处于其最大/最小出力方位内
P j , t m i n ≤ ∑ i ∈ j P i , t ≤ P j , t m a x P_{j,t} ^ {min} \leq \sum_{i\in j }P_{i,t} \leq P_{j,t}^{max} Pj,tmin≤i∈j∑Pi,t≤Pj,tmax
- P j , t m a x P_{j,t}^{max} Pj,tmax , P j , t m i n P_{j,t}^{min} Pj,tmin 为机组群 j j j 在时段 t t t 的最大最小出力
(7)机组爬坡约束
爬坡:功率的变化率,模型中简化,一旦开机功率就可以达到最小出力
上爬坡:
P i , t − P i , t − 1 ≤ Δ P i U α i , t − 1 + P i , t m i n ( α i , t − α i , t − 1 ) + P i , t m a x ( 1 − α i , t ) P_{i,t}- P_{i,t-1} \leq \Delta P_i^U \alpha_{i,t-1} + P_{i,t}^{min}(\alpha_{i,t}-\alpha_{i,t-1}) + P_{i,t}^{max}(1-\alpha_{i,t}) Pi,t−Pi,t−1≤ΔPiUαi,t−1+Pi,tmin(αi,t−αi,t−1)+Pi,tmax(1−αi,t)
下爬坡:
P i , t − 1 − P i , t ≤ Δ P i D α i , t − P i , t m i n ( α i , t − α i , t − 1 ) + P i , t m a x ( 1 − α i , t − 1 ) P_{i,t-1}- P_{i,t} \leq \Delta P_i^D \alpha_{i,t} - P_{i,t}^{min}(\alpha_{i,t}-\alpha_{i,t-1}) + P_{i,t}^{max}(1-\alpha_{i,t-1}) Pi,t−1−Pi,t≤ΔPiDαi,t−Pi,tmin(αi,t−αi,t−1)+Pi,tmax(1−αi,t−1)
- Δ P i U \Delta P_i^U ΔPiU 机组 i i i 最大上爬坡速率
- Δ P i D \Delta P_i^D ΔPiD 机组 i i i 最大下爬坡速率
(8)机组最小连续开停时间约束
由于火电机组的物理属性及实际运行需要,要求火电机组满足最小连续开机/停机时间。
T i , t D − ( α i , t − α i , t − 1 ) T D ≥ 0 T_{i,t}^D - (\alpha_{i,t} - \alpha_{i,t-1})T_D \geq 0 Ti,tD−(αi,t−αi,t−1)TD≥0
T i , t U − ( α i , t − 1 − α i , t ) T U ≥ 0 T_{i,t}^U - (\alpha_{i,t-1} - \alpha_{i,t})T_U \geq 0 Ti,tU−(αi,t−1−αi,t)TU≥0
- T U T_U TU , T D T_D TD 为机组最小连续开机时间和最小连续停机时间(是时段的倍数?)
- T i , t U T_{i,t}^U Ti,tU , T i , t D T_{i,t}^D Ti,tD 为机组 i i i 在时段 t t t 已经连续开机的时间和连续停机的时间
T i , t U = ∑ k = t − T U t − 1 α i , k T_{i,t}^U = \sum_{k=t-T_U}^{t-1}\alpha_{i,k} Ti,tU=k=t−TU∑t−1αi,k
T i , t D = ∑ k = t − T D t − 1 ( 1 − α i , k ) T_{i,t}^D = \sum_{k=t-T_D}^{t-1}(1-\alpha_{i,k}) Ti,tD=k=t−TD∑t−1(1−αi,k)
(9)网络约束
支路潮流约束
线路和变压器容量极限、断面的传输极限约束等
X t , l = ∑ i = 1 N G l − i P i , t + ∑ j = 1 N T G 1 − j T j , t − ∑ u = 1 M G l − u L u , t − ∑ k = 1 K G l − k D k , t , ∀ t ∈ T , ∀ l X_{t,l} = \sum_{i=1}^{N}G_{l-i}P_{i,t} + \sum_{j=1}^{NT}G_{1-j}T_{j,t} - \sum_{u=1}^M G_{l-u}L_{u,t}-\sum_{k=1}^K G_{l-k}D_{k,t},\quad \forall t\in T, \forall l Xt,l=i=1∑NGl−iPi,t+j=1∑NTG1−jTj,t−u=1∑MGl−uLu,t−k=1∑KGl−kDk,t,∀t∈T,∀l
- X t , l X_{t,l} Xt,l 为支路 l l l 的潮流
- N N N 机组总数
- N T NT NT 联络线总数
- M M M 市场用户数
- K K K 系统节点数
- G l − i G_{l-i} Gl−i 机组 i i i 所在节点对线路 l l l 的发电机输出功率转移分布因子
- G l − j G_{l-j} Gl−j 联络线 j j j 所在节点对线路 l l l 的发电机输出功率转移分布因子
- G l − u G_{l-u} Gl−u 市场用户 u u u 所在节点对线路 l l l 的发电机输出功率转移分布因子
- G l − k G_{l-k} Gl−k 节点 k k k 所在节点对线路 l l l 的发电机输出功率转移分布因子
- D k , t D_{k,t} Dk,t 为接点 k k k 在时段 t t t 的母线负荷值
X m a x , l − ≤ X t , l − S L s + + S L s − ≤ X m a x , l + X_{max,l}^- \leq X_{t,l} - SL_s^{+} + SL_s^{-} \leq X_{max,l}^+ Xmax,l−≤Xt,l−SLs++SLs−≤Xmax,l+
- X m a x , l + X_{max,l}^+ Xmax,l+ 线路潮流正向极限
- X m a x , l − X_{max,l}^- Xmax,l− 线路潮流反向极限
线路、断面、变压器的传输容量都可以抽象为线路的传输容量约束。即各个时段的线路传输功率不可超过线路正反向容量
在实际应用中,线路的正反向容量一般都留有一定的裕度,即 X m a x , l + X_{max,l}^+ Xmax,l+ , X m a x , l − X_{max,l}^- Xmax,l− 乘以一定比例最为线路的正反向传输极限
断面潮流约束:
X m a x , s − ≤ ∑ l = 1 N L X t , l − S L s + + S L s − ≤ X m a x , s + X_{max,s}^- \leq \sum _{l=1}^{NL}X_{t,l} - SL_s^{+} + SL_s^{-} \leq X_{max,s}^+ Xmax,s−≤l=1∑NLXt,l−SLs++SLs−≤Xmax,s+
- ∑ l = 1 N L X t , l \sum _{l=1}^{NL}X_{t,l} ∑l=1NLXt,l 断面 s s s 的潮流
- N L NL NL 组成断面 s s s 的支路集合
- X m a x , s + X_{max,s}^+ Xmax,s+ 断面潮流正向极限
- X m a x , s − X_{max,s}^- Xmax,s− 断面潮流反向极限
机组出力表达式 ∀ i , t \forall i,t ∀i,t
P i , t = ∑ m = 1 N M P i , t , m P_{i,t} = \sum_{m=1}^{NM}P_{i,t,m} Pi,t=m=1∑NMPi,t,m
P i , m m i n ≤ P i , t , m ≤ P i , m m a x P_{i,m}^{min} \leq P_{i,t,m}\leq P_{i,m}^{max} Pi,mmin≤Pi,t,m≤Pi,mmax
- N M NM NM 机组报价总段数
- P i , t , m P_{i,t,m} Pi,t,m 机组 i i i 在时段 t t t 第 m m m 个出力区间中的中标电力
- P i , m m a x P_{i,m}^{max} Pi,mmax 、 P i , m m i n P_{i,m}^{min} Pi,mmin 机组 i i i 申报的第 m m m 个出力区间上下界
机组运行费用表达式:
C i , t ( P i , t ) = ∑ m = 1 N M C i , m P i , t , m C_{i,t}(P_{i,t}) = \sum_{m=1}^{NM}C_{i,m}P_{i,t,m} Ci,t(Pi,t)=m=1∑NMCi,mPi,t,m
- N M NM NM 为机组报价总段数
- C i , t , m C_{i,t,m} Ci,t,m 为机组 i i i 申报的第 $ m $ 个出力区间对应的能量价格
(10)同电厂一天最多开停一台机约束
同电厂下的 N N N 台可优化开机的自由优化机组一天最多开一台机约束:
∑ i = 1 N ∑ t = 1 96 U p i , t ≤ 1 \sum_{i=1}^{N} \sum_{t=1}^{96} Up_{i,t} \leq 1 i=1∑Nt=1∑96Upi,t≤1
- U p i , t Up_{i,t} Upi,t 为机组 i i i 在时段 t t t 的开机动作状态
同电厂下的 N N N 台可优化停机的自由优化机组一天最多停
∑ i = 1 N ∑ t = 1 96 D n i , t ≤ 1 \sum_{i=1}^N\sum_{t=1}^{96}Dn_{i,t}\leq 1 i=1∑Nt=1∑96Dni,t≤1
- D n i , t Dn_{i,t} Dni,t 为机组 i i i 在时段 t t t 的停机动态状态
(11)同电厂自由优化机组置换开停状态约束
同电厂下的可优化开机的自由优化机组和可优化停机的自由优化机组之间的约束
(同一天最多开停一个机组?)
∑ t = 1 96 ( U p i , t + D n j , t ) ≤ 1 \sum_{t=1}^{96} (Up_{i,t} + Dn_{j,t})\leq 1 t=1∑96(Upi,t+Dnj,t)≤1
- i i i 表示可优化开机的自由优化机组
- j j j 表示可优化停机的自由优化机组
- U p i , t Up_{i,t} Upi,t 为机组 i i i 在时段 t t t 的开机动作状态
- D n j , t Dn_{j,t} Dnj,t 为机组 j j j 在时段 t t t 的停机动作状态
该约束有 N ∗ M N*M N∗M 条, N N N 表示可优化开机的自由优化机组数(集合?), M M M 表示可优化停机的自由优化机组数
(12) 电厂最大最小运行方式
对于每一个时段
t
t
t ,同电厂下开机火电数需要满足该电厂的最大,最小运行机组数
P
p
,
t
m
i
n
≤
∑
i
∈
P
p
α
i
,
t
≤
P
p
,
t
m
a
x
P_{p,t}^{min} \leq \sum_{i\in P_{p}} \alpha _{i,t} \leq P_{p,t}^{max}
Pp,tmin≤i∈Pp∑αi,t≤Pp,tmax
- i i i 表示电厂 p p p 内机组,
- α i , t \alpha _{i,t} αi,t 机组 i i i 在时段 t t t 的启停状态
- P p , t m i n P_{p,t}^{min} Pp,tmin 表示电厂 p p p 在时段 t t t 内最小运行机组数
- P p , t m a x P_{p,t}^{max} Pp,tmax 表示电厂 p p p 在时段 t t t 内最小运行机组数
(13) 机组群最大最小运行方式
对于每个时段
t
t
t ,机组群下开机火电数需要满足该机组群的最大,最小运行机组数
U
G
g
,
t
m
i
n
≤
∑
i
∈
U
G
g
α
i
,
t
≤
U
G
g
,
t
m
a
x
UG_{g,t}^{min} \leq \sum _{i\in UG_g} \alpha _{i,t} \leq UG_{g,t}^{max}
UGg,tmin≤i∈UGg∑αi,t≤UGg,tmax
- i i i 表示机组群 g g g 内机组,
- α i , t \alpha _{i,t} αi,t 机组 i i i 在时段 t t t 的启停状态
- U G g , t m i n UG_{g,t}^{min} UGg,tmin 表示机组群 g g g 在时段 t t t 内最小运行机组数
- U G g , t m a x UG_{g,t}^{max} UGg,tmax 表示机组群 g g g 在时段 t t t 内最小运行机组数
(14)最小开机容量
对于每个时段
t
t
t 系统开机火电机组容量需要满足最小开机容量需求
∑
(
C
i
m
a
x
×
α
i
,
t
)
≥
T
C
t
m
i
n
\sum (C_{i}^{max}\times \alpha _{i,t}) \geq TC_{t}^{min}
∑(Cimax×αi,t)≥TCtmin
- C i m a x C_{i}^{max} Cimax 机组 i i i 的额定容量
- T C t m i n TC_{t}^{min} TCtmin 系统在 t t t 时段的最小火电开机容量需求
(15) 市场用户运行约束
市场用户用电限制约束
L u , t m i n ≤ L u , t ≤ L u , t m a x , ∀ t ∈ T , u ∈ M L_{u,t}^{min} \leq L_{u,t} \leq L_{u,t}^{max}, \quad \forall t\in T, u\in M Lu,tmin≤Lu,t≤Lu,tmax,∀t∈T,u∈M
- L u , t L_{u,t} Lu,t 为第 u u u 个市场用户 t t t 时段的用电
- 市场用户的用电上限 L u , t m a x L_{u,t}^{max} Lu,tmax 为 t t t 时段的报价最后一段分段终止电力
- L u , t m i n L_{u,t}^{min} Lu,tmin 为市场下限为0