风电与储能系统互补下的火电机组组合优化决策方法与流程

本发明涉及电力系统运行与控制技术领域，尤其涉及一种风电与储能系统互补下的火电机组组合优化决策方法。

背景技术：

机组组合是电力系统运行与控制中关键而核心的问题，其实质是制定日或周前的满足发电、负荷及备用的机组启停计划。

近年来，迫于能源需求和环境保护的双重压力，风电等可再生能源发电得到了迅猛发展。其在发挥节能减排效益的同时，由于其存在固有的强不确定性，使电力系统安全经济运行面临挑战。为适应这一挑战，储能等各类互补设施不断涌现，其目的是促进可再生能源发电的消纳，如何实现这一最佳消纳，成为机组组合决策面临的新问题。

针对此问题，“利用风蓄联合削峰的电力系统经济调度”和“含风储混合系统的多目标机组组合优化模型及求解”建立了含风储系统多目标机组组合模型，分析了引入储能系统对火电机组运行成本、能源利用效率等因素的影响。为增强机组组合决策结果应对不确定性的有效性，“robustunitcommitmentwithwindpowerandpumpedstoragehydro”针对含风电和抽水蓄能系统提出了鲁棒机组组合决策方法，能够有效计及风电极限误差的影响，实现抽水蓄能机组对风电波动的抑制。“风电-火电-抽水蓄能联合优化机组组合模型”采用置信区间法对风电不确定性进行建模，考虑了相邻时段风电极限场景间系统调节能力的牵制关系，可发挥抽水蓄能机组削峰填谷作用，促进风电的消纳。依据随机规划理论，“储能系统对含风电的机组组合问题影响分析”针对集中式储能建立了含风储系统随机机组组合模型，储能系统的峰谷调节作用不仅能够降低系统运行成本，同时能够从一定程度上减轻风电不确定性的影响，减少机组组合无解的情况。“计及风险备用的大规模风储联合系统广域协调调度”针对配置在各个风电场中的分散式储能，提出了计及风险备用约束的随机机组组合模型，可实现不同风电场储能系统间的协调。

从目前的研究状况看，储能等补偿设施的加入无疑提升了机组组合应对不确定性的能力。然而，由其性质可知，储能等补偿设施是依赖电网而存在的，其可控性源自对电网中电能的时空平移，在时变过程中尽可能发挥其位置、数量的有效改变，才能够起到应对不确定性的作用。若在研究周期内将其按预期进行锁定，转而当预期出现偏差时，完全由火电机组承担，导致机组必须同时应对负荷与可再生能源两类预期外的不确定性，造成机组提供更多的备用，或者出现切负荷、弃风光发电现象。

技术实现要素：

本发明的目的就是为了解决上述问题，提供一种风电与储能系统互补下的火电机组组合优化决策方法，建立了机组组合两层优化决策模型，通过系统分解协调的原理，下层问题在上层问题给定的优化信息引导下，通过上、下层问题的交替迭代对该模型予以求解，决策储能系统充/放电功率、调控范围及机组启停方案，从而更有效提升风电的接纳能力。

为了实现上述目的，本发明采用如下技术方案：

风电与储能系统互补下的火电机组组合优化决策方法，包括上层和下层两层问题优化的过程：

上层问题以设定周期内火电机组运行成本最小为目标函数，建立上层问题约束条件；

下层问题以储能系统对电网中电能时空平移和提供备用所得收益最大为目标函数，建立下层问题约束条件；

利用上层问题和下层问题交替迭代的方式求解。

上层问题的目标函数具体为：其中，n为机组数量；ai、bi、ci分别为机组i成本特性系数；pi,t为机组i在t时段输出功率；ui,t为机组i在t时段启停状态，0表示停运，1表示运行；si,t表示机组i在t时段启动成本。

上层问题约束条件包括：功率平衡约束条件，机组输出功率上下限约束条件，机组爬坡速率约束条件，机组最小开停机时间约束条件以及备用约束条件。

下层问题的目标函数具体为：其中λt为t时段边际发电成本，μt、γt分别表示储能系统提供上调备用和下调备用的价格，分别为t时段储能系统充、放电功率计划值，分别为t时段储能应对不确定性的调控范围。

取值根据下层优化问题决策结果确定，在上层问题中为给定常数；

取值也是根据下层优化问题决策结果确定，在上层问题中为给定常数。

下层问题约束条件包括：储能系统充/放电功率约束条件，储能系统充/放电功率调控范围约束条件，储能系统存储能量约束条件以及储能系统存储能量调控范围约束条件。

λt反映t时段发电与负荷有功平衡的紧张情况，在负荷高峰时段λt较大，负荷低谷时λt较小；

μt、γt数值越大，储能系统预留的调控范围越大，储能系统提供备用的收益越大，同时上层问题中机组所需提供的备用容量越小。

在求解过程中：根据上层问题决策结果，计算t时段边际发电成本、储能系统提供上调备用和下调备用的影子价格，对下层问题中储能系统充/放电功率及调控范围进行优化决策。

在求解过程中：依据储能系统与风电互补的调控策略，对风电波动进行补偿，提高其不确定性的可信度。

在求解过程中：若机组运行成本在迭代过程中不再下降，或达到给定的迭代步数或计算时间，则认为迭代收敛；否则，返回到上层问题重新进行优化计算。

本发明的有益效果：

(1)本发明模型中储能系统不仅可参与期望的功率平衡约束，还能够通过与风电的互补起到备用的作用，有直接消除不确定性的能力；

(2)在提升接纳不确定性能力的基础上，避免火电机组频繁的启停；

(3)在决策中考虑了电力系统频率的调节效应，使决策免除保守性，更符合实际。

附图说明

图1为应对风电不确定性的基本原理；

图2储能系统存储电量；

图3为各时段储能系统存储电量计划值；

图4为储能系统存储电量实际变化情况；

图5为为机组组合总运行成本随储能系统容量的变化趋势。

具体实施方式

下面结合附图与实施例对本发明作进一步说明。

1储能系统的互补作用

本发明重点考虑风电，风电功率预测误差的分布与预测时间尺度有关，机组组合决策中预测时间尺度一般至少为24h，此时风电功率可认为服从正态分布。由预测技术获得其概率分布规律后，考虑风电波动需满足的置信度要求，可获得调度决策需满足的风电功率预测值、最大值以及最小值。

风电功率不确定性强，若仅依靠常规火电机组应对风电波动，则造成其备用增加，频繁启动等现象，不仅代价高，也会造成弃风电或切负荷等情况出现。

储能系统运行灵活，在时变过程中应尽可能发挥其位置、数量的有效改变，以起到减轻或消除被动量不确定性的作用。为此，本发明提出了储能与不确定性互补的机制：储能的决策不仅参与期望功率平衡约束，还将起到备用的作用，以减少常规火电机组备用容量，提高其运行效率。图1以储能应对风电不确定性为例进行说明。

图1给出了风电功率概率分布密度曲线，对应t时段储能应对风电不确定性的调控范围。

由图1可见，阴影部分所对应的风电不确定性可通过储能在运行过程中调整充/放电功率予以平抑，而由常规火电机组集中应对风电预测值邻域内有限的不确定性。若储能在满足基本用电需求的同时，仍存在较高冗余的调节空间，那么在理想情况下，储能可完全消除风电的不确定性。

对应上述思路，在不考虑波动对应概率数值前提下，储能与风电互补过程中的调控策略可表示为，

式(1)中：pess,t为调度给定的t时段储能系统充电功率计划值，负值表示放电；为t时段储能系统补偿风电功率波动后的实际充电功率；为随机波动量，表示t时段可能的风电功率；分别为t时段风电功率预测在一定置信范围内的最大值和最小值；为便于表示，对应图1，式(1)中at等于bt等于

2机组组合的两层优化决策模型

2.1上层问题

上层问题以研究周期内火电机组运行成本最小为目标，其目标函数可表示为，

式(2)中，n为机组数量；ai、bi、ci分别为机组i成本特性系数；pi,t为机组i在t时段输出功率；ui,t为机组i在t时段启停状态，0表示停运，1表示运行；si,t表示机组i在t时段启动成本。

约束条件：

(1)功率平衡约束

式(3)中，pw,t、pd,t为t时段风电、常规负荷预测值；分别为t时段储能系统充、放电功率计划值，其取值根据下层优化问题决策结果确定，在上层问题中为给定常数。

(2)机组输出功率上下限约束

式(4)中，pi^max、pi^min分别为机组i允许的最大、最小输出功率。

(3)机组爬坡速率约束

pi,t+1-pi,t≤ri^up·ui,t+pi^max(1-ui,t)(5)

pi,t-pi,t+1≤ri^dn·ui,t+pi^max(1-ui,t+1)(6)

式(5)、式(6)中，分别为机组i输出功率最大的向上、向下调节速率。

(4)机组最小开停机时间约束

式(7)、式(8)，分别为机组i已开机、停机的时段数；ti,on、ti,off分别为机组i最小开、停机时间。

(5)备用约束

式(9)—式(13)中，分别为机组i在t时段为应对风电和负荷不确定性所达到的最大功率值和最小功率值；分别为t时段负荷功率在一定置信范围内的最大值、最小值；βt为t时段系统频率调节系数，以mw/hz为单位；ki为机组i频率调节系数；为s情景下t时段电网频率相对额定频率的偏差，s＝1表示负荷位于波动范围上限而风电位于波动范围下限的严峻情景，s＝2表示负荷位于波动范围下限而风电位于波动范围上限的严峻情景；δf^max为正常运行时电网频率允许的最大偏差；δi为应对风电波动，agc机组i输出功率基点最大调整量，对于非agc机组，δi等于0；分别为t时段储能应对不确定性的调控范围，其取值根据下层优化问题决策结果确定，在上层问题中为给定常数。

传统机组组合决策中应对负荷与风电波动主要对应agc机组二次调节能力和调节容量的配置。实际伴随负荷与风电波动，电网频率允许偏离额定值运行，系统中运行的所有机组以及负荷均可通过参与频率调节应对不确定性。进一步考虑一次调节作用，相比视电网频率处于额定值不变，仅考虑agc机组调节的方式，可增强系统应对不确定性的能力，避免调度无解而实际可行的情况。t时段系统频率调节系数可表示为，

式(14)中，dt为t时段负荷频率调节系数。

2.2下层问题

下层问题目标函数可表示为，

式(15)中，λt、μt、γt为上层问题给定的协调信息。

λt取为t时段边际发电成本，可根据上层问题决策结果求得。由于λt可反映t时段发电与负荷有功平衡的紧张情况，一般在负荷高峰时段λt较大，负荷低谷时λt较小。将λt作为协调信息，引导储能系统充/放电功率的优化决策，使之在负荷高峰(λt较大)时放电，在负荷低谷(λt较小)时充电，按此不仅能够促进上层问题有功平衡的实现，而且能够使储能系统在对电网中电能时空平移过程中实现自身收益最大化。类似，μt、γt分别表示储能系统提供上调备用和下调备用价格，该数值越大，储能系统预留的调控范围越大，储能系统提供备用的收益越大，同时上层问题中机组所需提供的备用容量越小，可有效避免机组频繁启停。

按上述，由于本发明研究的是储能系统配置给定下的短期运行优化问题，因此下层问题不考虑储能系统配置的投资成本，而是以储能系统对电网中电能时空平移和提供备用所得收益最大为目标，在这一目标实现过程中，借助上层问题给定的协调信息的引导，可促进上层问题有功平衡的实现以及减轻火电机组的备用负担，提高机组组合决策的经济性。

约束条件主要包括：

(1)储能系统充/放电功率约束

式(16)—式(18)中，分别为储能系统允许的最大放电、充电功率；分别为储能系统处于放电或充电的标志，式(18)保证储能系统不能够同时处于充电或放电状态。

(2)储能系统充/放电功率调控范围约束

当储能系统放电或充电功率等于其最大允许值时，储能系统将丧失上调或下调的调控能力。本发明通过式(19)、式(20)对储能系统调控范围与其计划充/放电功率的牵制关系进行了描述。

当t时段储能系统处于充电状态时，应有在应对风电不确定性过程中储能系统可采取的调控手段包括：风电低于预测值或负荷高于预测值时减少充电功率，为避免储能系统在短时内反复充/放电影响电池寿命，充电功率最多减少至0，调控范围受式(19)的限制；风电高于预测值或负荷低于预测值时增加充电功率，充电功率最大为调控范围受式(20)的限制。

类似，当t时段储能系统处于放电状态时，应有风电低于预测值或负荷高于预测值时储能系统可增加放电功率，放电功率最大为调控范围受式(19)的限制；风电高于预测值或负荷低于预测值时应减少放电功率，为避免短时内反复充/放电，放电功率最多减少至0，调控范围受式(20)的限制。

(3)储能系统存储能量约束

e^min≤eess,t≤e^max(22)

式(21)、式(22)中，eess,t为t时段储能系统存储能量值；η^ch、η^d分别为储能系统充电、放电效率；e^max、e^min为储能系统存储能量允许的最大值、最小值。

(4)储能系统存储能量调控范围约束

由于储能系统存储能量处于最大值或最小值时，将丧失部分调控能力，为此应使存储能量预留一定的调节空间，需满足如下约束，

式(23)、式(24)为储能系统能够在时变过程中发挥应对不确定性的作用预留了存储能量的调节空间。

2.3求解流程

本发明模型采用上、下层问题交替迭代的方式求解，具体求解流程如下：

(1)对上层问题中机组启停计划、备用配置进行优化决策，储能系统充/放电功率、风电补偿后的功率等下层问题互补结果在上层问题中为给定常数。

(2)根据上层问题决策结果，计算λt、μt、γt。对下层问题中储能系统充/放电功率及调控范围进行优化决策。

(3)依据储能系统与风电互补的调控策略，对风电波动进行补偿，提高其不确定性的可信度。

(4)若机组运行成本在迭代过程中不再下降，或达到给定的迭代步数或计算时间，则认为迭代收敛。否则，返回到上层问题重新进行优化计算。

3算例分析

3.1算例概述

采用10机组系统对本发明模型进行仿真验证。火电机组容量、爬坡速率等参数见”acomputationallyefficientmixed-integerlinearformulationforthethermalunitcommitmentproblem”，风电数据和常规负荷数据见表1，储能系统参数见表2。将机组二次耗量特性函数进行分段线性化(三分段，分段大小相等)。系统频率为50hz，正常运行时频率允许变化范围为50±0.1hz。

表1负荷和风电功率数据

表2储能系统参数

为说明本发明模型有效性，考虑以下三种情况：

(1)情况1：机组组合决策中不考虑储能系统情况；

(2)情况2：机组组合决策中考虑储能系统，但储能系统不参与系统备用，通过协调优化确定其充/放电功率后直接下发执行；

(3)情况3：本发明方法。

3.2情况1决策结果

在该种情况下，根据计算结果，机组1、6、7在所有时段均投入运行，其它7台机组启停计划如表3所示。

表3情况1对应的机组启停计划

总运行成本为472414美元，共涉及13次机组启动过程。在研究周期内agc机组(机组6—10)启停频繁，在时段12-13需要5台agc机组全部投入运行，而在时段16-19仅需要2台agc机组。可见，为满足风电和常规负荷不确定性，仅依靠火电的机组组合方法需要在研究周期内进行较为频繁的机组启停，运行成本较高。

3.3情况2决策结果

在该种情况下，不考虑储能系统与风电互补应对不确定性的作用，风电和常规负荷不确定性完全由火电机组承担。根据计算结果，机组1、6、7在所有时段均投入运行，机组3在所有时段均停运，其它6台机组启停计划如表4所示。

表4情况2对应的机组启停计划

总运行成本为464784美元，共涉及11次机组启动过程。相较不考虑储能系统参与调度决策的情况1，运行成本和机组启停次数有所减少。但由于备用仍完全由火电机组承担，在运行过程中agc机组(机组6—10)启停仍较为频繁。

3.4情况3决策结果

采用本发明模型，经过3次交替迭代后收敛。根据本发明模型计算结果，机组1在所有时段均投入运行，机组3、9、10在所有时段均停运，其它6台机组的启停计划如表5所示。

表5情况3对应的机组启停计划

总运行成本为420671美元，共涉及9次机组启停过程，相较前两种情况均有较为明显的减少。同时，在研究周期内需要投入运行的agc机组数量明显减少，在部分时段(如时段13—19)可完全通过风储互补的方式，并结合系统频率调节效应实现不确定性的消纳，减少了agc机组调节负担。

根据本发明模型计算结果，各时段储能系统存储电量计划值如图3所示。

图3中，虚线分别为储能系统存储电量允许的最大值、最小值。由图3可知，储能系统在净负荷较小的时段2～5充电，在净负荷较高的时段9～11、时段20～21放电，可起到峰谷调节的作用。同时，储能系统存储电量基本维持在允许范围的中间值附近，从而保证储能系统留有调节余量，能够通过调整充电或放电功率应对风电不确定性。

为说明本发明模型决策结果的有效性，通过生产模拟的方式随机生成一个风电和负荷波动的情景，该情景具体数据见附录表6，以此为例，对储能系统在互补过程中的调控策略以及电网频率变化情况进行仿真计算。随着储能系统备用容量的释放，各时段存储电量实际变化情况如图4所示。

表6负荷和风电功率随机模拟数据

由图4可知，按照本发明模型决策结果，储能系统在实际运行中可顺利发挥应对不确定性的作用，未出现调节能力受限的情况。同时，仅需要在时段2、时段13、时段15调用电网频率调节资源，在频率允许范围内偏离额定值运行。对比表4、表5可知，在上述时段利用电网频率调节效应，可充分发挥非agc机组及负荷本身的调节作用，减少了agc机组调节容量配置。

3.5应对不确定性能力对比

为进一步说明本发明模型应对不确定性的能力，将风电功率相对预测值的偏差以一定步长逐渐增加，计算结果如表7所示。

表7应对不确定性能力对比

由表7中情况1可知，随着风电不确定性的增强，仅依靠常规火电机组难以实现风电的消纳。对比情况1、情况2可知，储能系统的引入从一定程度上增加了系统应对不确定性的能力，但情况2对应运行成本明显高于本发明方法，并且随着风电不确定性的增强，出现机组组合无可行解的情况。而在表中所示的不确定性范围内，本发明方法均可得到可行解，增强了系统应对不确定性的能力，有利于可再生能源发电的消纳。

3.6储能系统容量对机组组合决策的影响

为分析本发明模型下储能系统容量对机组组合决策的影响，对表2所示的除充放电效率参数外的其他储能系统参数按等比例放缩，然后依据本发明模型进行计算，所得的机组组合总运行成本随储能系统容量变化的趋势如图5所示。

由图5可知，机组组合总运行成本随着储能系统容量的增加呈现单调递减的趋势，且这一变化过程是非线性的，原因在于该问题的非凸性，当储能容量达到某些临界值时，如60mwh，会引起机组启停的变化，从而造成总运行成本波动性的变化；当然，如果储能系统容量位于机组启停相对稳定的区间，如330～390mwh，由于总运行成本的降低来源于机组输出功率的优化，未涉及机组启停调整，使这一变化的趋势则相对平缓。最后，当储能系统容量达到某一限值(本发明为570mwh)，由于系统应对不确定性的能力相当充裕甚至出现冗余，再增加储能系统容量也无法引起电网运行经济性的进一步提升。可见，储能系统容量配置需要与电网不确定运行环境相适应，既要避免因容量不足引起的电网运行经济性降低，又要避免因配置容量过大出现冗余造成投资浪费。

针对如何在火电机组组合决策中利用储能系统应对可再生能源发电不确定性的问题，本发明提出储能系统应对风电不确定性的互补机制，建立了火电机组组合两层优化决策的模型。上层问题以机组组合成本最小为目标，下层问题以储能系统与风电、负荷间的互补程度为满足对象，其中计及了agc机组和非agc机组的特性，以及系统频率调节效应的作用。通过系统分解协调的原理，下层问题在上层问题给定的优化信息引导下，通过上、下层问题的交替迭代对该模型予以求解，决策储能系统充/放电功率、调控范围及机组启停方案，从而更有效提升风电的接纳能力。

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。