[发明专利]一种风光火储协调优化运行方法

摘要

一种风光火储协调优化运行方法，首先获得火电机组、背压式热电联产机组及抽汽式热电联产机组开机状态序列，其次获得风电机组及光伏机组可用功率序列，然后以电力系统运行成本最小为目标，考虑电力系统实时功率平衡约束、机组功率约束、火电机组爬坡约束、热电联产机组供热约束、抽水蓄能电站约束及电力系统备用约束，建立风光火储协调运行模型，最后，对不同模拟时段内电力系统状态求解风光火储协调运行模型，得到机组功率安排；若模型不存在可行解，则引入弃风功率及弃光功率，对风光火储协调运行模型修改后再进行求解；本发明实现以电力系统运行成本最小的风光火储协调优化运行，能够为电力系统运行和调度提供指导。

主权项

1.一种风光火储协调优化运行方法，其特征在于，包括以下步骤:1)输入原始数据，包括电力系统电负荷数据、热负荷数据、火电机组强迫停运率及平均修复时间、背压式热电联产机组强迫停运率及平均修复时间、抽汽式热电联产机组强迫停运率及平均修复时间、历史风速数据及历史太阳辐照度数据；定义最大模拟周期数Dmax和一个模拟周期内最大模拟时段数Tmax，且设定模拟周期D＝1，模拟时段T＝1；2)基于状态持续时间采样法获得模拟周期D内火电机组开机状态序列、背压式热电联产机组开机状态序列及抽汽式热电联产机组开机状态序列；3)基于逆变换法获得模拟周期D内风电机组可用功率序列及光伏机组可用功率序列；4)获取模拟周期D的模拟时段T内的电力系统状态，包括电负荷数据、热负荷数据、火电机组开机状态、背压式热电联产机组开机状态及抽汽式热电联产机组开机状态、风电机组可用功率及光伏机组可用功率；5)对步骤4)中的电力系统状态求解风光火储协调运行模型，若模型存在可行解，则说明模拟周期D的模拟时段T内电力系统能够完全消纳风电机组可用功率及光伏机组可用功率，进入步骤8)；若模型不存在可行解，说明模拟周期D的模拟时段T内电力系统不能完全消纳风电机组可用功率及光伏机组可用功率，进入步骤6)；风光火储协调运行模型由目标函数和约束条件组成，具体模型为：5.1)目标函数风光火储协调运行模型的目标函数为模拟周期D内电力系统运行成本最小，电力系统运行成本包括火电机组运行成本、背压式热电联产机组运行成本、抽汽式热电联产机组运行成本及抽水蓄能机组运行成本，具体为：火电机组i在模拟时段T内的运行成本为：背压式热电联产机组m在模拟时段T内的运行成本为：抽汽式热电联产机组n在模拟时段T内的运行成本为：抽水蓄能电站运行成本为：式中：N_th为火电机组数，为背压式热电联产机组数，为抽汽式热电联产机组数，N_p为抽水蓄能电站数，为火电机组i在模拟时段T内的功率，为背压式热电联产机组m在模拟时段T内的供电功率，为抽汽式热电联产机组n在模拟时段T内的供电功率，为抽汽式热电联产机组n在模拟时段T内的供热功率，a_i,b_i和c_i为火电机组耗量特性参数，为背压式热电联产机组耗量特性参数，为抽汽式热电联产机组耗量特性参数，为抽水蓄能电站u启机一次的成本，为抽水蓄能电站u停机一次的成本，为抽水蓄能电站u在模拟时段T内的抽水运行状态指示变量，为抽水蓄能电站u在模拟时段T+1内的发电运行状态指示变量，当抽水蓄能电站u在模拟时段T内处于抽水运行状态时，值为1，值为0；当抽水蓄能机组电站u在模拟时段T内处于发电运行状态时，值为1，值为0；5.2)约束条件风光火储协调运行模型的约束条件包括电力系统实时功率平衡约束、机组功率约束、火电机组爬坡约束、热电联产机组供热约束、抽水蓄能电站约束及电力系统备用约束，具体为：5.2.1)电力系统实时功率平衡约束电力系统实时功率平衡约束包括实时电功率平衡约束和实时热功率平衡约束；实时电功率平衡约束为：实时热功率平衡约束为：式中：为抽水蓄能电站u在模拟时段T内的抽水功率或发电功率，为风电机组j在模拟时段T内的可用功率，为光伏机组k在模拟时段T内的可用功率，P_L,T为电力系统在模拟时段T内的电负荷，H_L,T为电力系统在模拟时段T内的热负荷，N_w为风电机组数，N_s为光伏机组数，为背压式热电联产机组m在模拟时段T内的供热功率；5.2.2)机组功率约束机组功率约束包括火电机组功率约束、背压式热电联产机组供电功率约束、抽汽式热电联产机组供电功率约束、抽汽式热电联产机组供热功率约束及抽水蓄能电站功率约束；火电机组功率约束为：背压式热电联产机组供电功率约束为：抽汽式热电联产机组供电功率约束为：抽汽式热电联产机组供热功率约束为：抽水蓄能电站功率约束为：式中：为火电机组i的最小功率，为火电机组i的最大功率，为背压式热电联产机组m的最小供电功率，为背压式热电联产机组m的最大供电功率，为抽汽式热电联产机组n的最小供电功率，为抽汽式热电联产机组n的最大供电功率，H_nmin为抽汽式热电联产机组n的最小供热功率，H_nmax为抽汽式热电联产机组n的最大供热功率，为抽水蓄能电站u的最小抽水功率，为抽水蓄能电站u的最大抽水功率，为抽水蓄能电站u的最小发电功率，为抽水蓄能电站u的最大发电功率；5.2.3)火电机组爬坡约束式中：Rup,i为火电机组i的向上爬坡率，Rdown,i为火电机组i的向下爬坡率，Δt为每个模拟时段步长；5.2.4)热电联产机组供热约束热电联产机组供热约束包括背压式热电联产机组供热约束和抽汽式热电联产机组供热约束；背压式热电联产机组m的供热约束为线性等式约束，具体为：抽汽式热电联产机组n的供热约束为线性不等式约束，具体为：式中：和分别为背压式热电联产机组电/热线性关系参数，由背压式热电联产机组m的运行工况图确定，及分别为抽汽式热电联产机组电/热线性关系参数，由抽汽式热电联产机组n的运行工况图确定；5.2.5)抽水蓄能电站约束抽水蓄能电站约束包括上水库库容约束、抽水蓄能电站能量守恒约束、抽水蓄能电站末状态约束及经济性约束；上水库库容约束为：式中：E_umin为抽水蓄能电站u的上水库最小储能，E_umax为抽水蓄能电站u的上水库最大储能，为抽水蓄能电站u在模拟时段T结束时的上水库储能，计算公式为：式中：E_u,0为抽水蓄能电站u上水库的初始储能，为抽水蓄能电站u在抽水运行状态下的能量转换效率，为抽水蓄能电站u在发电运行状态下的能量转换效率；抽水蓄能电站能量守恒约束为：抽水蓄能电站末状态约束为：经济性约束为：且v≠w5.2.6)系统备用约束式中：R为电力系统备用容量；6)引入风电机组弃风功率及光伏机组弃光功率，对步骤5)中的风光火储协调运行模型进行修改；6.1)修改步骤5.1)中的目标函数，修改后的风光火储协调运行模型的目标函数为：式中：为弃风惩罚系数，为弃光惩罚系数，为风电机组j在模拟时段T内的弃风功率，为光伏机组k在模拟时段T内的弃光功率；6.2)在步骤5.2)约束条件的基础上增加风电机组弃风功率约束与光伏机组弃光功率约束；风电机组j弃风功率约束为：光伏机组k弃光功率约束为：6.3)在步骤5.2)约束条件的基础上增加风电机组可用功率、风电机组实际消纳功率与风电机组弃风功率之间的关系约束及光伏机组可用功率、光伏机组实际消纳功率与光伏机组弃光功率之间的关系约束；风电机组j可用功率、风电机组j实际消纳功率与风电机组j弃风功率之间的关系约束为：光伏机组k可用功率、光伏机组k实际消纳功率与光伏机组k弃光功率之间的关系约束为：其中：为模拟时段T内风电机组j实际消纳功率，为模拟时段T内光伏机组k实际消纳功率；6.4)修改步骤5.2.1)中电力系统实时电功率平衡约束，修改后的电力系统实时电功率平衡约束为：7)对步骤4)中的电力系统状态求解步骤6)中修改后的风光火储协调运行模型，得到火电机组功率、背压式热电联产机组供电功率及供热功率、抽汽式热电联产机组供电功率及供热功率、抽水蓄能机组抽水功率或发电功率、风电机组实际消纳功率及光伏机组实际消纳功率；8)统计模拟周期D的模拟时段T内火电机组功率、背压式热电联产机组供电功率及供热功率、抽汽式热电联产机组供电功率及供热功率、抽水蓄能机组抽水功率或发电功率、风电机组可用功率及实际消纳功率、光伏机组可用功率及实际消纳功率；9)如果模拟周期D的模拟时段T等于最大模拟时段数Tmax，则进入步骤11)，否则，返回步骤4)，进入下一模拟时段；10)如果模拟周期D等于最大模拟周期数Dmax，则进入步骤11)，否则，返回步骤2)，进入下一模拟周期；11)输出各模拟周期内火电机组时序功率、背压式热电联产机组时序供电功率及时序供热功率、抽汽式热电联产机组时序供电功率及时序供热功率、抽水蓄能机组时序抽水功率及时序发电功率、风电机组可用功率及实际消纳功率、光伏机组可用功率及实际消纳功率。