现货市场环境下含抽水蓄能机组的日内滚动优化调度方法

1.本发明涉及电力系统发电调度领域，具体涉及一种利用在现货市场环境下研究含抽水蓄能机组的日内滚动优化调度方法。


背景技术：

2.在现有电网调控体系下，抽水蓄能电站的调度运行由电网调控中心“按需调度”，抽水蓄能电站时刻处于热备用状态，运行管理人员也需要24小时处于备战状态，造成人员劳动强度和设备运行成本增大，不利于抽水蓄能电站的整体效益提升。同时，受新能源出力预测和负荷预测水平影响，日前市场与实时市场出清结果偏差增大，需要更加完备的日内市场参与，减轻实时市场的调整负担。关于抽水蓄能机组参与含火电/可再生能源电力系统的协调优化调度策略的研究工作大部分以理论研究为主，因此，在电力市场环境下满足实际应用要求的抽水蓄能机组的日内滚动优化运行方法和策略还需深入研究，实现发电计划跟随日前制定计划。


技术实现要素：

3.本发明所要解决的技术问题是提供一种现货市场环境下研究含抽水蓄能机组的日内优化调度方法，解决现有的抽水蓄能电站获得的收益难以覆盖投资和运行成本，电网调用积极性不高，同时新能源出力预测、负荷预测水平受限，日前制定的发电计划在系统实际运行时跟随效果并不理想，引入市场机制，基于所设计的日内市场交易流程，在日前电量出清计划的基础上。利用更准确的新能源出力和负荷超短期预测(15min)信息，构建了含抽水蓄能机组的日内滚动优化调度模型，并对相关市场主体所申报的电量价格信息线性化处理，制定机组一个交易周期内(4h)的计划运行曲线。设计了滚动偏差限制，保障日前计划的有效跟随，减少实时调度的调节成本。在满足用电需求的前提下，保障电网安全稳定运行，促进可再生能源发展和消纳。
4.为解决上述技术问题，本发明所采用的技术方案是：
5.现货市场环境下含抽水蓄能机组的日内滚动优化调度方法，其特征在于，包括以下步骤：
6.步骤一、设计日内现货交易流程；
7.步骤二、获取电力系统的基础数据；
8.步骤三、建立含抽水蓄能机组的日内滚动优化模型；
9.步骤四、对各市场主体申报的电量价格分段函数作线性化处理；
10.步骤五、对日内滚动优化模型进行求解，得到该交易周期内的机组计划运行曲线。
11.上述的步骤一中，日内现货市场交易流程包括交易准备、市场出清、交易信息发布及执行三个阶段。
12.上述的步骤二中，基础数据包括火电机组以15min为一点的日前机组96点出清计划以及运行状态信息、火电机组信息、新能源场站信息、抽水蓄能电站信息、区域负荷预测
信息、日内市场的市场主体申报的电量价格信息。
13.上述步骤三中，建立含抽水蓄能机组的日内滚动优化模型的目标函数如下：
[0014][0015]
其中，p
i,t
为火电机组i在时段t的日内计划出力值，为火电机组i在时段t的增、减出力，f
i,t
为新能源站i在时段t的日内出清出力值，为新能源站i在时段t的增、减出力，为抽水蓄能机组i在发电状态时在时段t的增、减出力，p
g,i,t
为抽水蓄能机组i在发电状态时在时段t的日内出清计划出力值，利用机组日前制定计划和日内计划的差值可以表示该机组在某一时刻的增减出力，c
i,t
(p
i,t
)、c
i,t
(f
i,t
)、c
i,t
(p
g,i,t
)是与火电机组i、新能源站i、抽水蓄能机组i申报的各段出力区间和对应的能量价格有关的多段线性函数，为抽水蓄能机组发电的状态变量，其中1为发电状态，0为停机状态，n
g
为系统火电机组数，n
m
为系统新能源站数量，n
h
为抽水蓄能机组数，t1为一个日内交易周期16点的出清时段。
[0016]
上述步骤三中，所建模型的约束条件包括：系统功率平衡约束、机组出力约束、火电机组爬坡约束、系统旋转备用约束、抽水蓄能机组发电功率约束、水库容量约束、水库容量平衡约束、滚动计划值偏差约束。
[0017]
1)上述系统功率平衡约束的表达式为：
[0018][0019]
式中p
l,i,t
为抽水蓄能机组i在时段t的抽水功率，为抽水蓄能机组抽水的状态变量；1为抽水状态，0为停机状态，l
in,t
、l
off,t
分别为该交易时段t的某区域负荷超短期预测值及外购电计划值。
[0020]
2)上述机组出力约束的表达式为：
[0021]
i
i,t
*p
imin
≤p
i,t
≤i
i,t
*p
imax
[0022]
f
i,t
≤q
i,t
[0023]
式中p
imax
、p
imin
为火电机组i的出力上下限，q
i,t
为新能源站i在时段t的最大预测出力；i
i,t
为火电机组i在时段t的启停状态变量；0为停机状态，1为运行状态；由日前计划决定。
[0024]
3)上述火电机组爬坡约束的表达式为：
[0025]
p
i,t
‑
p
i,t
‑1≤
△
p
iu
[0026]
p
i,t
‑1‑
p
i,t
≤
△
p
id
[0027]
△
p
iu
、
△
p
id
为火电机组i的最大上、下爬坡速率。
[0028]
4)上述系统旋转备用约束的表达式为：
[0029][0030][0031]
式中l
u
％、l
d
％为负荷高于、低于预期百分比，f
u
％、f
d
％为新能源出力高于、低于预期百分比。
[0032]
5)上述抽水蓄能机组，在日内滚动优化调度阶段，主要以抽水蓄能机组来平衡新能源出力预测误差，因此需要对抽水蓄能电站的出力及库容等加以约束：
[0033]
5.1)抽水蓄能机组为满功率抽水，仅考虑发电时机组出力约束，约束的表达式为：
[0034][0035]
式中为抽水蓄能机组i发电时的最大出力；
[0036]
5.2)水库容量约束的表达式为：
[0037]
e
min
≤e
t
≤e
max
[0038]
式中e
t
为抽水蓄能电站在时段t的水库容量，利用电量来衡量，e
max
、e
min
为抽水蓄能电站上水库的最大、最小库容；
[0039]
5.3)水库容量平衡约束的表达式为：
[0040][0041]
式中，η
l
、η
g
为抽水蓄能机组的抽水、发电效率；
[0042]
5.4)抽水蓄能机组不能同时发电和抽水，抽发互斥约束，约束的表达式为：
[0043][0044]
6)滚动调度的修正值与原计划值具有关联性，需将每台机组修正后的出力值与日前计划出力偏差控制在一定范围内，保证日内调度和日前调度之间的协调性，滚动计划偏差约束的表达式为：
[0045][0046]
为机组i日前计划出力值。p
ie
为机组i的额定功率，ξ为约束乘子，调整其大小可以改变日内调度最大允许的修正偏差。
[0047]
上述的步骤四中，考虑到常规火电机组的发电成本具有二次函数形式的表达式，电力现货市场允许发电机组提交多段报价，即出力－价格阶梯曲线(分段函数)，且各段报价单调非递减；所建立的模型主要利用0
‑
1整数线性规划方法求解，考虑到采用线性规划来解决问题得到的优化结果比较理想，在出清机组计划运行曲线时需要对电量价格分段函数作线性化处理，以一台机组为例，其电量价格具体线性化过程如下：
[0048][0049]
其中，为机组第i段出力区间对应的起点、终点，c
i
为机组第i段出力区间对应的电量价格，按容量大小依次递增，c为机组在某时刻出力对应线性化后的电量价格信息，考虑n＝3段，即将机组出力范围分为3个区间段进行申报，引入0
‑
1整数变量k
i
对电量价格线性化。
[0050]
优选地，上述的含抽水蓄能机组的日内滚动优化模型的求解，采用0
‑
1混合整数线性规划方法。
[0051]
本发明提供的一种现货市场环境下研究含抽水蓄能机组的日内优化调度方法，具有以下有益效果：
[0052]
1)本发明考虑了可再生能源出力预测的精度会随着预测时间缩短而提高的特点，利用日前电能量交易出清结果为计算依据，提出了滚动计划修订模型，通过设计滚动偏差值约束，较好的对发电计划跟随问题进行处理，较好消纳新能源；
[0053]
2)本发明考虑了现阶段抽水蓄能电站经营模式导致成本回收困难背景，将其引入现货市场，提出的含抽水蓄能机组的日内滚动优化模型中考虑了抽水蓄能运行效益，抽水蓄能在按需调度过程中保证了新能源最大消纳量；
[0054]
3)本发明是在现货市场背景下组织展开的，涉及相关市场主体申报电量价格信息过程，引入0
‑
1变量后较好的处理了发用电方分段报价的问题。
附图说明
[0055]
下面结合附图和实施例对本发明作进一步说明：
[0056]
图1为现货市场环境下含抽水蓄能机组的日内滚动优化调度的流程图；
[0057]
图2为日内市场交易组织流程图；
[0058]
图3为日内滚动时序图；
[0059]
图4为反映新能源消纳情况示意图；
[0060]
图5为滚动修订后的机组运行曲线示意图；
[0061]
图6为抽水蓄能每台机组运行状态和出力情况示意图。
具体实施方式
[0062]
以下结合附图和实施例详细说明本发明技术方案。
[0063]
现货市场环境下含抽水蓄能机组的日内滚动优化调度方法，其特征在于，包括以下步骤：
[0064]
步骤一、设计日内现货交易流程；
[0065]
1)交易准备阶段主要工作：日前交易出清闭市后，所有直调发电机组、场(站)需在
日内市场向市场运营机构提供其运行参数，各市场主体基于日前现货市场交易出清结果，针对日内电网负荷波动、临时安全调整等需要，申报运行日上调、下调容量及价格信息。市场运营机构综合有关信息，提前4h向相关市场主体发布日内现货市场的运行参数及边界条件信息；
[0066]
2)市场出清阶段主要工作：以4h为一个日内交易周期，每15min为一个交易出清时段，每个交易周期含有16个交易出清时段，日内市场每天开放6次，以前两个交易周期为例，图3为具体滚动时序图，t为交易执行时刻，基于日前市场优化出清结果、实时更新的超短期负荷预测和新能源出力预测信息，综合考虑相关网络安全约束，进行日内优化出清，在(t
‑
140)min前，出清该日内交易时段的机组运行计划曲线(16点)，在(t
‑
110)min前，市场运营机构对该时段的机组出力曲线进行交流潮流安全校核，若不满足该约束，则在模型中添加相应约束条件，重新进行上述过程直到满足，得到该日内交易时段的分时发电出力曲线；
[0067]
3)信息发布及执行阶段主要工作：(t
‑
80)min前，市场技术支持系统向相关市场主体发布日内交易出清结果。各市场主体审核后返回确认信息，对出清结果有异议的市场主体可向市场运营机构提出解释申请；
[0068]
步骤二、获取电力系统的基础数据；
[0069]
1)火电机组日前96点运行计划每台火电机组运行状态i
i,t
；
[0070]
2)新能源场(站)信息：新能源场(站)数量n
m
，新能源场(站)编号i，新能源超短期功率预测信息q
i,t
,新能源出力高于、低于预期百分比f
u
％、f
d
％；
[0071]
3)某区域负荷超短期预测值l
in,t
、外购电计划值l
off,t
,负荷高于、低于预期百分比l
u
％、l
d
％；
[0072]
4)电力系统中火电机组信息：火电机组总台数n
g
，火电机组编号i，每台火电机组的最大有功出力p
imax
、最小有功出力p
imin
以及额定功率p
ie
，上爬坡速率
△
p
iu
、下爬坡速率
△
p
id
；
[0073]
5)抽水蓄能电站信息：抽水蓄能机组数n
h
，抽水蓄能机组编号i，抽水蓄能机组最大发电功率抽水功率抽水蓄能电站上水库最大库容e
max
、最小库容e
min
、抽水效率η
l
、发电效率η
g
；
[0074]
6)市场主体申报的电量价格信息：火电机组c
i,t
(p
i,t
)、新能源场站c
i,t
(f
i,t
)、抽水蓄能机组c
i,t
(p
g,i,t
)；
[0075]
步骤三、建立含抽水蓄能机组的日内滚动优化模型,该模型由目标函数和约束条件构成；
[0076]
步骤四、对各市场主体申报的电量价格分段函数作线性化处理；
[0077]
步骤五、通过线性规划求解模型，该模型为0
‑
1混合整数线型规划模型，调用成熟的线型规划软件包进行求解，如glpk,得到该交易周期内的机组计划运行曲线。
[0078]
上述步骤三中，确定市场环境下含抽水蓄能机组的日内滚动优化模型目标函数中，考虑系统整体效益最大，保证了抽水蓄能电站的一部分收益,同时，以日前计划出清的机组运行结果信息为计算依据，并在约束条件中添加了滚动计划修正值偏差约束，在使新能源被系统最大化消纳的情况下可保证修正后的机组出力计划较好的跟随日前计划，减少部分调节成本,最后，对于日内市场价格信息，也做了线性化处理，较好的处理了市场主体
分段报价的问题；
[0079]
综上，目标函数综合考虑系统运行整体效益最大，目标函数如下：
[0080][0081]
其中，p
i,t
为火电机组i在时段t的日内计划出力值，为火电机组i在时段t的增、减出力。f
i,t
为新能源站i在时段t的日内计划出力值，为新能源站i在时段t的增、减出力，为抽水蓄能机组i在发电状态时在时段t的增、减出力，p
g,i,t
为抽水蓄能机组i在发电状态时在时段t的日内计划出力值，利用机组日前制定计划和日内计划的差值可以表示该机组在某一时刻的增减出力，c
i,t
(p
i,t
)、c
i,t
(f
i,t
)、c
i,t
(p
g,i,t
)是与火电机组i、新能源站i、抽水蓄能机组i申报的各段出力区间和对应的能量价格有关的多段线性函数，为抽水蓄能机组发电的状态变量，其中1为发电状态，0为停机状态，n
g
为系统火电机组数，n
m
为系统新能源站数量，n
h
为抽水蓄能机组数，t1为一个日内交易周期16点的出清时段。
[0082]
上述步骤三中，所建模型的约束条件包括：系统功率平衡约束、机组出力约束、火电机组爬坡约束、系统旋转备用约束、抽水蓄能机组发电功率约束、水库容量约束、水库容量平衡约束、滚动计划值偏差约束。
[0083]
1)上述系统功率平衡约束的表达式为：
[0084][0085]
式中p
l,i,t
为抽水蓄能机组i在抽水状态时的抽水功率，为抽水蓄能机组抽水的状态变量；1为抽水状态，0为停机状态，l
in,t
、l
off,t
分别为该交易时段t的某区域负荷超短期预测值及外购电计划值。
[0086]
2)上述机组出力约束的表达式为：
[0087]
i
i,t
*p
imin
≤p
i,t
≤i
i,t
*p
imax
[0088]
f
i,t
≤q
i,t
[0089]
式中p
imax
、p
imin
为火电机组i的出力上下限，q
i,t
为新能源站i在时段t的最大预测出力；i
i,t
为火电机组i在时段t运行的状态变量；0为停机状态，1为运行状态；由日前计划决定。
[0090]
3)上述火电机组爬坡约束的表达式为：
[0091]
p
i,t
‑
p
i,t
‑1≤
△
p
iu
[0092]
p
i,t
‑1‑
p
i,t
≤
△
p
id
[0093]
△
p
iu
、
△
p
id
为火电机组i的最大上、下爬坡速率。
[0094]
4)上述系统旋转备用约束的表达式为：
[0095][0096][0097]
式中l
u
％、l
d
％为负荷高于、低于预期百分比，f
u
％、f
d
％为新能源出力高于、低于预期百分比。
[0098]
5)上述抽水蓄能机组，在日内滚动优化调度阶段，主要以抽水蓄能机组来平衡新能源出力预测误差，因此需要对抽水蓄能电站的出力及库容等加以约束：
[0099]
5.1)抽蓄机组为满功率抽水，仅考虑发电状态机组出力约束，约束的表达式为：
[0100][0101]
式中为抽水蓄能机组i发电时的最大出力；
[0102]
5.2)水库容量约束的表达式为：
[0103]
e
min
≤e
t
≤e
max
[0104]
式中e
t
为抽水蓄能电站t时刻上水库的容量，利用电量来衡量，e
max
、e
min
为抽水蓄能电站上水库的最大、最小库容；
[0105]
5.3)水库容量平衡约束的表达式为：
[0106][0107]
式中，η
l
、η
g
为抽水蓄能机组的抽水、发电效率；
[0108]
5.4)抽水蓄能机组不能同时发电和抽水，抽发互斥约束，约束的表达式为：
[0109][0110]
6)滚动调度的修正值与原计划值具有关联性，需将每台机组修正后的出力值与日前计划出力偏差控制在一定范围内，保证日内调度和日前调度之间的协调性，滚动计划偏差约束的表达式为：
[0111][0112]
式中为机组i日前计划出力值。p
ie
为机组i的额定功率，ξ为约束乘子，调整其大小可以改变日内调度最大允许的修正偏差。
[0113]
上述的步骤四中，电量价格分段函数线性化处理如下所示：
[0114][0115]
式中，为机组申报的第i段出力区间对应的起点、终点，c
i
为机组第i段出
力区间对应的电量价格，按容量大小依次递增，c为机组在某时刻出力对应线性化后的电量价格信息，考虑n＝3段，引入0
‑
1整数变量k
i
对电量价格线性化。
[0116]
优选地，上述的含抽水蓄能机组的日内滚动优化模型的求解，采用0
‑
1混合整数线性规划方法。
[0117]
采用某省级电网数据为计算依据，该电网潮流母线节点数1925个，线路851条，变压器934台，发电机561台(含风光电厂等值发电机)；等值风电场224个，光伏电站95个，直调风电装机容量为18446.34mw，直调光伏装机容量为7240.1mw。该地区建有一座抽蓄电站，安装有4台300mw抽水蓄能机组，总装机容量为1200mw；最大储能量为满发7.4h，满功率抽水需要9.15h。滚动周期设为8:15～12:00，即时刻33～48，从d5000平台自动获取优化算法所需要的各类数据后对所提模型及优化方法进行验证，所建优化模型为0
‑
1混合整数线性规划模型，共有57466个变量，其中整数变量9200个，58193个约束条件，采用自由软件项目gnu的线性规划软件包glpk进行求解，因为问题规模较大，且多为0
‑
1变量，采用了fischetti
‑
monaci邻近搜索启发式算法，该算法能够快速改进混合整数规划问题的可行解，使得一些需要很长时间才能得到最优解的问题快速获得次优解，大大提高了求解速率。
[0118]
采用本发明的现货市场环境下含抽水蓄能机组的日内滚优化模型进行求解，所得新能源消纳情况、机组修正后的计划运行曲线和抽水蓄能机组运行状态结果如图4、图5和图6所示。
[0119]
从图4可以看出滚动周期内的所有时段新能源超短期预测出力和实际出力一致，一个周期内并无弃风、光现象，说明了所建模型在大部分时间下可对可再生能源消纳情况较好。
[0120]
风电场、光伏电站、抽水蓄能电站运行曲线以及滚动修正后的火电机组计划运行曲线如图5所示，可以看出在9:00～12:00,负荷降低且新能源总出力基本保持稳定，系统为了吸收更多的新能源，抽水蓄能电站机组基本处于抽水状态运行，火电机组出力较稳定，降低了电机组调峰压力，实现各种能源的优化配置。
[0121]
抽水蓄能电站4台机组运行状态和出力情况如图6所示，上水库储水量的变化是通过电量来衡量的，由4台机组运行状态共同决定，从图中可以看出在8:15
‑
8:45蓄能蓄电站处于发电状态，在其他时段由于新能源出力缓慢增加，抽水蓄能电站安排相应机组抽水运行，随着新能源出力情况的变化，抽水蓄能机组按需调度来保证新能源被尽数消纳。