CN115296345A - 一种燃气发电机组启停-出力-备用联合优化方法及装置 - Google Patents

一种燃气发电机组启停-出力-备用联合优化方法及装置 Download PDF

Info

Publication number
CN115296345A
CN115296345A CN202210648680.2A CN202210648680A CN115296345A CN 115296345 A CN115296345 A CN 115296345A CN 202210648680 A CN202210648680 A CN 202210648680A CN 115296345 A CN115296345 A CN 115296345A
Authority
CN
China
Prior art keywords
output
day
unit
constraint
new energy
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210648680.2A
Other languages
English (en)
Other versions
CN115296345B (zh
Inventor
郑晓东
陈明媛
周保荣
祁乐
张旻钰
宋吉峰
杨有慧
禤培正
程兰芬
王嘉阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CSG Electric Power Research Institute
Guangxi Power Grid Co Ltd
Original Assignee
CSG Electric Power Research Institute
Guangxi Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CSG Electric Power Research Institute, Guangxi Power Grid Co Ltd filed Critical CSG Electric Power Research Institute
Priority to CN202210648680.2A priority Critical patent/CN115296345B/zh
Publication of CN115296345A publication Critical patent/CN115296345A/zh
Application granted granted Critical
Publication of CN115296345B publication Critical patent/CN115296345B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

本发明涉及电力系统运行控制技术领域,公开了一种燃气发电机组启停‑出力‑备用联合优化方法及装置。本发明通过获取电力系统的有功调度参数作为优化的基础数据,通过获取新能源出力预测误差的历史数据,构造计划日的新能源出力预测误差凸包络场景集,进而构建计及多阶段调度非预期性的燃气发电机日前有功计划模型,构建考虑该凸包络场景集的时序解耦日内有功调度模型,进而根据日前有功计划模型和日内有功调度模型得到燃气发电机组启停、出力和备用联合优化模型,最后根据对该联合优化模型的求解结果编制燃气发电机组启停、出力和备用计划。本发明提出的方法及装置能够适应概率化电力平衡模式,实现燃气发电机组日前计划的科学和精准编制。

Description

一种燃气发电机组启停-出力-备用联合优化方法及装置
技术领域
本发明涉及电力系统运行控制技术领域,尤其涉及一种燃气发电机组启 停-出力-备用联合优化方法及装置。
背景技术
大力发展风电和光伏等新能源发电是构建新型电力系统的必要举措。然 而,风电和光伏等新能源的发电具有间歇性、波动性及随机性等特征。为可靠 地保证电力系统的实时有功平衡,电力系统需要配套一定数量的灵活调节资 源。燃气发电机组作为灵活调节资源的一种,具有启停灵活、爬坡速度快、调 峰能力强等优点,其相比燃煤发电机组更具清洁性,且可以通过富余新能源制 氢或制甲烷等方式减少发电的碳排放。因此,燃气发电机组将在新型电力系统 中扮演重要的角色。
构建新型电力系统,除了需要在规划阶段配置足够的灵活调节资源,还需 要在运行阶段的电力平衡过程中,充分考虑新能源发电的强不确定性及其给 灵活调节资源调控引入的复杂性,并相应地改进有功计划和调度的模型或方 法。在日前发电和备用计划编制过程中,考虑新能源发电的多个场景进行概率 化的电力平衡,可以得到更加经济和可靠的调度方案,促进新能源的高效消纳。 概率化电力平衡的难点之一是非预期性的建模。非预期性的含义是,在一个多 时段的随机决策过程中,以当前时段为参考点,具有同样历史随机量实现值的 场景必须对应同一个决策方案;非预期性也可以理解为,在一个多时段的随机 决策问题中,决策方案的制定不能够以掌握整个时域过程的随机量实现值为 假设和前提。
燃气发电机组的启停、出力和备用联合优化是一个典型的多时段决策问 题。具体而言,在十五分钟时间尺度上,燃气发电机组的日内启停需要考虑最 小开停机时间约束,燃气发电机组提供灵活爬坡备用的能力还受到启停和相 邻时段出力计划的影响。目前,燃气发电机组的日前发电和备用计划编制均以 确定性的预测场景为基础和假设,并没有考虑随机场景下燃气发电机组调控 的非预期性。因此,在概率化电力平衡模式下,如何在燃气发电机组启停、出 力和备用的联合优化过程中考虑非预期性,是实现日前发电和备用计划科学、 精确编制必须解决的技术难题。
发明内容
本发明提供了一种燃气发电机组启停-出力-备用联合优化方法及装置,解 决了目前燃气发电机组的日前发电和备用计划编制未考虑到随机场景下燃气 发电机组调控的非预期性,计划编制的精确性有待提高的技术问题。
本发明第一方面提供一种燃气发电机组启停-出力-备用联合优化方法,包 括:
获取电力系统的有功调度参数以及新能源关于出力预测值和出力真实值 的历史数据;
引入燃气发电机组关于启停、出力和备用的决策变量,基于所述决策变量 和所述有功调度参数,构建包含机组最小开停机时间精确约束和灵活爬坡备 用约束在内的日前有功计划约束,根据所述日前有功计划约束建立燃气发电 机组日前有功计划模型;
根据所述历史数据确定计划日的新能源出力预测误差集合,将所述新能 源出力预测误差集合的高维凸包作为凸包络场景集,根据所述有功调度参数 构建考虑所述凸包络场景集的时序解耦日内有功调度模型;
根据新能源出力削减约束、功率平衡约束和线路潮流约束建立系统侧模 型,联合所述燃气发电机组日前有功计划模型、所述时序解耦日内有功调度模 型和所述系统侧模型,以计划成本加新能源出力预测场景下的运行成本之和 最小为优化目标,构建燃气发电机组启停、出力和备用的联合优化模型;
采用场景生成法求解所述联合优化模型,根据得到的求解结果生成燃气 发电机组的启停、出力和备用计划。
根据本发明第一方面的一种能够实现的方式,所述根据所述历史数据确 定计划日的新能源出力预测误差集合,包括:
在新能源的所有历史出力预测值中确定与计划日新能源出力预测值相似 度最大的数据子集;
将所述数据子集对应的历史出力真实值作为计划日可能的出力真实值, 计算计划日可能的出力真实值与对应的计划日出力预测值的差值;
将所述差值作为新能源出力预测误差,构建计划日的新能源出力预测误 差集合。
根据本发明第一方面的一种能够实现的方式,所述引入燃气发电机组关 于启停、出力和备用的决策变量,基于所述决策变量和所述有功调度参数,构 建包含机组最小开停机时间精确约束和灵活爬坡备用约束在内的日前有功计 划约束,包括:
引入二进制变量
Figure BDA0003687047600000031
Figure BDA0003687047600000032
Figure BDA0003687047600000033
Figure BDA0003687047600000034
分别表示机组g在时段t开机状 态变量实际可取得的最小值和最大值;引入日前的连续变量
Figure BDA0003687047600000035
Figure BDA0003687047600000036
Figure BDA0003687047600000037
Figure BDA0003687047600000038
分别表示机组g在日内时段t出力的基准上限和基准下限;引入日前的连 续变量
Figure BDA0003687047600000039
Figure BDA00036870476000000310
Figure BDA00036870476000000311
Figure BDA00036870476000000312
分别表示机组g在时段t提供的向上灵活爬坡备用 容量和向下灵活爬坡备用容量;
构建机组最小开机时间精确约束:
Figure BDA00036870476000000313
式中,
Figure BDA00036870476000000314
表示机组g在时段t+τ开机状态变量实际可取得的最小值,
Figure BDA00036870476000000315
表示机组g在时段t-1开机状态变量实际可取得的最小值,
Figure BDA00036870476000000316
表示燃 气发电机组的指标集合,T表示计划时段的指标集合,
Figure BDA00036870476000000317
表示机组g的最小 开机时间;
构建机组最小停机时间精确约束:
Figure BDA00036870476000000318
式中,
Figure BDA00036870476000000319
表示机组g在时段t+τ开机状态变量实际可取得的最大值,
Figure BDA00036870476000000320
表示机组g在时段t-1开机状态变量实际可取得的最大值,
Figure BDA00036870476000000321
表示机组 g的最小停机时间;
构建向上灵活爬坡备用约束:
Figure BDA00036870476000000322
式中,
Figure BDA00036870476000000323
表示机组g在日内时段t-1出力的基准下限;
构建向下灵活爬坡备用约束:
Figure BDA0003687047600000041
式中,
Figure BDA0003687047600000042
示机组g在日内时段t-1出力的基准上限。
根据本发明第一方面的一种能够实现的方式,所述引入燃气发电机组关 于启停、出力和备用的决策变量,基于所述决策变量和所述有功调度参数,构 建包含机组最小开停机时间精确约束和灵活爬坡备用约束在内的日前有功计 划约束,还包括:
引入日前的二进制变量yg,t,当yg,t取1时表示机组g在日内的时段t允许停 机,取0时表示机组g在日内的时段t不允许停机或者因为基准计划已经是停 机状态而无法进行关停;引入日前的连续变量
Figure BDA0003687047600000043
Figure BDA0003687047600000044
表示机组g在日内时 段t出力的实际下限;
构建机组开机状态变量最小值的表达式:
Figure BDA0003687047600000045
式中,xg,t表示机组g在时段t开停机状态的基准计划;
构建机组开机状态变量最大值的表达式:
Figure BDA0003687047600000046
构建机组出力上下限关联关系约束:
Figure BDA0003687047600000047
式中,P g 表示机组g的最小技术出力,
Figure BDA0003687047600000048
表示机组g的装机容量;
构建机组出力上爬坡约束:
Figure BDA0003687047600000049
其中,
Figure BDA00036870476000000410
表示机组g在日内时段t-1出力的实际下限,
Figure BDA00036870476000000411
表示机组g 的最大上爬坡能力;
构建机组出力下爬坡约束:
Figure BDA00036870476000000412
式中,
Figure BDA00036870476000000413
表示机组g的最大下爬坡能力;
构建机组出力实际下限的表达式:
Figure BDA0003687047600000051
构建旋转备用约束:
Figure BDA0003687047600000052
式中,
Figure BDA0003687047600000053
表示机组g在时段t提供的旋转备用容量。
根据本发明第一方面的一种能够实现的方式,所述根据所述有功调度参 数构建考虑所述凸包络场景集的时序解耦日内有功调度模型,包括:
引入日内实时运行的二进制变量zg,t,当zg,t取1时表示机组g在时段t进行 关停,取0时表示机组g在时段t保持开机运行状态;引入新能源出力预测误差 随机向量ξ;引入连续变量pg,t(ξ),pg,t(ξ)表示在场景ξ发生的情况下,机组g在 时段t的实际出力;
构建机组停机约束:
Figure BDA0003687047600000054
构建机组出力上下限约束:
Figure BDA0003687047600000055
根据本发明第一方面的一种能够实现的方式,所述根据新能源出力削减 约束、功率平衡约束和线路潮流约束建立系统侧模型,包括:
构建新能源出力削减约束为:
Figure BDA0003687047600000056
式中,
Figure BDA0003687047600000057
表示在新能源出力场景ξ下系统需要削减的新能源功率, Ww,t表示新能源场站w在时段t的预测出力,ξw,t表示新能源场站w在时段t的实 际出力预测误差,
Figure BDA0003687047600000058
表示新能源场站的指标集合;
构建功率平衡约束为:
Figure BDA0003687047600000061
式中,Dd,t表示负荷节点d在时段t的负荷需求,
Figure BDA0003687047600000062
表示负荷节点的指标集 合;
构建线路潮流约束为:
Figure BDA0003687047600000063
式中,Fl表示线路l的潮流限值,
Figure BDA0003687047600000064
表示机组到线路l的功率转移因子,
Figure BDA0003687047600000065
表示负荷到线路l的功率转移因子,
Figure BDA0003687047600000066
表示新能源场站到线路l的功率转移 因子,
Figure BDA0003687047600000067
表示输电线路的指标集合。
根据本发明第一方面的一种能够实现的方式,所述采用场景生成法求解 所述联合优化模型,包括:
S701,考虑预测场景,形成一个松弛模型,调用最优化求解器获得一个待 校验的日前计划方案;
S702,给定日前计划方案,在新能源出力预测值的基础上叠加预测误差场 景集中的误差值,形成多个新能源出力场景,并调用最优化求解器逐个验证日 内实时运行是否存在可行方案;
S703,如果存在不可行方案,则将产生不可行方案的新能源出力场景扩充 到所述松弛模型中,更新得到新的待校验日前计划方案;
S704,重复步骤S702-S703,直至日前计划方案通过校验。
本发明第二方面提供一种燃气发电机组启停-出力-备用联合优化装置,包 括:
获取模块,用于获取电力系统的有功调度参数以及新能源关于出力预测 值和出力真实值的历史数据;
第一模型构建模块,用于引入燃气发电机组关于启停、出力和备用的决策 变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停机时间 精确约束和灵活爬坡备用约束在内的日前有功计划约束,根据所述日前有功 计划约束建立燃气发电机组日前有功计划模型;
第二模型构建模块,用于根据所述历史数据确定计划日的新能源出力预 测误差集合,将所述新能源出力预测误差集合的高维凸包作为凸包络场景集, 根据所述有功调度参数构建考虑所述凸包络场景集的时序解耦日内有功调度 模型;
模型联合模块,用于根据新能源出力削减约束、功率平衡约束和线路潮流 约束建立系统侧模型,联合所述燃气发电机组日前有功计划模型、所述时序解 耦日内有功调度模型和所述系统侧模型,以计划成本加新能源出力预测场景 下的运行成本之和最小为优化目标,构建燃气发电机组启停、出力和备用的联 合优化模型;
优化求解模块,用于采用场景生成法求解所述联合优化模型,根据得到的 求解结果生成燃气发电机组的启停、出力和备用计划。
根据本发明第二方面的一种能够实现的方式,所述第二模型构建模块包 括:
确定单元,用于在新能源的所有历史出力预测值中确定与计划日新能源 出力预测值相似度最大的数据子集;
计算单元,用于将所述数据子集对应的历史出力真实值作为计划日可能 的出力真实值,计算计划日可能的出力真实值与对应的计划日出力预测值的 差值;
构建单元,用于将所述差值作为新能源出力预测误差,构建计划日的新能 源出力预测误差集合。
根据本发明第二方面的一种能够实现的方式,所述第一模型构建模块包 括:
第一变量引入单元,用于引入二进制变量
Figure BDA0003687047600000071
Figure BDA0003687047600000072
Figure BDA0003687047600000073
Figure BDA0003687047600000074
分别表 示机组g在时段t开机状态变量实际可取得的最小值和最大值;引入日前的连 续变量
Figure BDA0003687047600000075
Figure BDA0003687047600000076
Figure BDA0003687047600000077
Figure BDA0003687047600000078
分别表示机组g在日内时段t出力的基准上限和基 准下限;引入日前的连续变量
Figure BDA0003687047600000079
Figure BDA00036870476000000710
Figure BDA00036870476000000711
Figure BDA00036870476000000712
分别表示机组g在时段t提供 的向上灵活爬坡备用容量和向下灵活爬坡备用容量;
第一构建单元,用于构建机组最小开机时间精确约束:
Figure BDA0003687047600000081
式中,
Figure BDA0003687047600000082
表示机组g在时段t+τ开机状态变量实际可取得的最小值,
Figure BDA0003687047600000083
表示机组g在时段t-1开机状态变量实际可取得的最小值,
Figure BDA0003687047600000084
表示燃 气发电机组的指标集合,T表示计划时段的指标集合,
Figure BDA0003687047600000085
表示机组g的最小 开机时间;
第二构建单元,用于构建机组最小停机时间精确约束:
Figure BDA0003687047600000086
式中,
Figure BDA0003687047600000087
表示机组g在时段t+τ开机状态变量实际可取得的最大值,
Figure BDA0003687047600000088
表示机组g在时段t-1开机状态变量实际可取得的最大值,
Figure BDA0003687047600000089
表示机组 g的最小停机时间;
第三构建单元,用于构建向上灵活爬坡备用约束:
Figure BDA00036870476000000810
式中,
Figure BDA00036870476000000811
表示机组g在日内时段t-1出力的基准下限;
第四构建单元,用于构建向下灵活爬坡备用约束:
Figure BDA00036870476000000812
式中,
Figure BDA00036870476000000813
示机组g在日内时段t-1出力的基准上限。
根据本发明第二方面的一种能够实现的方式,所述第一模型构建模块还 包括:
第二变量引入单元,用于引入日前的二进制变量yg,t,当yg,t取1时表示机 组g在日内的时段t允许停机,取0时表示机组g在日内的时段t不允许停机或者 因为基准计划已经是停机状态而无法进行关停;引入日前的连续变量
Figure BDA00036870476000000814
Figure BDA00036870476000000815
表示机组g在日内时段t出力的实际下限;
第五构建单元,用于构建机组开机状态变量最小值的表达式:
Figure BDA00036870476000000816
式中,xg,t表示机组g在时段t开停机状态的基准计划;
第六构建单元,用于构建机组开机状态变量最大值的表达式:
Figure BDA0003687047600000091
第七构建单元,用于构建机组出力上下限关联关系约束:
Figure BDA0003687047600000092
式中,P g 表示机组g的最小技术出力,
Figure BDA0003687047600000093
表示机组g的装机容量;
第八构建单元,用于构建机组出力上爬坡约束:
Figure BDA0003687047600000094
其中,
Figure BDA0003687047600000095
表示机组g在日内时段t-1出力的实际下限,
Figure BDA0003687047600000096
表示机组g 的最大上爬坡能力;
第九构建单元,用于构建机组出力下爬坡约束:
Figure BDA0003687047600000097
式中,
Figure BDA0003687047600000098
表示机组g的最大下爬坡能力;
第十构建单元,用于构建机组出力实际下限的表达式:
Figure BDA0003687047600000099
第十一构建单元,用于构建旋转备用约束:
Figure BDA00036870476000000910
式中,
Figure BDA00036870476000000911
表示机组g在时段t提供的旋转备用容量。
根据本发明第二方面的一种能够实现的方式,所述第二模型构建模块包 括:
第三变量引入单元,用于引入日内实时运行的二进制变量zg,t,当zg,t取1 时表示机组g在时段t进行关停,取0时表示机组g在时段t保持开机运行状态; 引入新能源出力预测误差随机向量ξ;引入连续变量pg,t(ξ),pg,t(ξ)表示在场 景ξ发生的情况下,机组g在时段t的实际出力;
第十二构建单元,用于构建机组停机约束:
Figure BDA00036870476000000912
第十三构建单元,用于构建机组出力上下限约束:
Figure BDA0003687047600000101
根据本发明第二方面的一种能够实现的方式,所述模型联合模块包括:
第十四构建单元,用于构建新能源出力削减约束为:
Figure BDA0003687047600000102
式中,
Figure BDA0003687047600000103
表示在新能源出力场景ξ下系统需要削减的新能源功率, Ww,t表示新能源场站w在时段t的预测出力,ξw,t表示新能源场站w在时段t的实 际出力预测误差,
Figure BDA00036870476000001011
表示新能源场站的指标集合;
第十五构建单元,用于构建功率平衡约束为:
Figure BDA0003687047600000104
式中,Dd,t表示负荷节点d在时段t的负荷需求,
Figure BDA0003687047600000105
表示负荷节点的指标集 合;
第十六构建单元,用于构建线路潮流约束为:
Figure BDA0003687047600000106
式中,Fl表示线路l的潮流限值,
Figure BDA0003687047600000107
表示机组到线路l的功率转移因子,
Figure BDA0003687047600000108
表示负荷到线路l的功率转移因子,
Figure BDA0003687047600000109
表示新能源场站到线路l的功率转移 因子,
Figure BDA00036870476000001010
表示输电线路的指标集合。
根据本发明第二方面的一种能够实现的方式,所述优化求解模块具体用 于:
S701,考虑预测场景,形成一个松弛模型,调用最优化求解器获得一个待 校验的日前计划方案;
S702,给定日前计划方案,在新能源出力预测值的基础上叠加预测误差场 景集中的误差值,形成多个新能源出力场景,并调用最优化求解器逐个验证日 内实时运行是否存在可行方案;
S703,如果存在不可行方案,则将产生不可行方案的新能源出力场景扩充 到所述松弛模型中,更新得到新的待校验日前计划方案;
S704,重复步骤S702-S703,直至日前计划方案通过校验。
本发明第三方面提供了一种燃气发电机组启停-出力-备用联合优化装置, 包括:
存储器,用于存储指令;其中,所述指令用于实现如上任意一项能够实现 的方式所述的燃气发电机组启停-出力-备用联合优化方法;
处理器,用于执行所述存储器中的指令。
本发明第四方面一种计算机可读存储介质,所述计算机可读存储介质上 存储有计算机程序,所述计算机程序被处理器执行时实现如上任意一项能够 实现的方式所述的燃气发电机组启停-出力-备用联合优化方法。
从以上技术方案可以看出,本发明具有以下优点:
1)本发明采用历史数据构建新能源出力预测误差的凸包络场景集,作为 概率化电力平衡的基础,有助于高效应对新能源的预测误差;
2)在日前计划模型中引入了燃气发电机组是否允许日内停机的决策变量 和最小开停机时间的精确约束,实现日内调度模型的时序解耦,进而有效计及 了调控的非预期性,保证概率化电力平衡模式下燃气发电机组模型的精确性;
3)利用机组出力上下限基准决策变量构造向上和向下灵活爬坡备用容量, 为概率化电力平衡模式下燃气发电机组调节贡献的评估提供了基础;
4)采用场景生成法求解日前-日内两阶段的燃气发电机组启停、出力和备 用联合优化模型,有效提高了计算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实 施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描 述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不 付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1为本发明一个可选实施例提供的一种燃气发电机组启停-出力-备用 联合优化方法的流程图;
图2为本发明一个可选实施例提供的凸包络场景集的示意图;
图3为本发明一个可选实施例提供的燃气发电机组启停-出力-备用联合 优化装置的原理框图。
附图标记:
1-获取模块;2-第一模型构建模块;3-第二模型构建模块;4-模型联合模 块;5-优化求解模块。
具体实施方式
本发明实施例提供了一种燃气发电机组启停-出力-备用联合优化方法及 装置,用于解决目前燃气发电机组的日前发电和备用计划编制未考虑到随机 场景下燃气发电机组调控的非预期性,计划编制的精确性有待提高的技术问 题。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结 合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描 述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施 例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提 下所获得的所有其它实施例,都属于本发明保护的范围。
本发明提供了一种燃气发电机组启停-出力-备用联合优化方法。
请参阅图1,图1示出了本发明实施例提供的一种燃气发电机组启停-出 力-备用联合优化方法的流程图。
本发明实施例提供的一种燃气发电机组启停-出力-备用联合优化方法,包 括步骤S1-S5。
步骤S1,获取电力系统的有功调度参数以及新能源关于出力预测值和出 力真实值的历史数据。
其中,新能源关于出力预测值和出力真实值的历史数据,包括各风电、光 伏场站所有的历史出力预测值和真实值时序数据。
步骤S2,引入燃气发电机组关于启停、出力和备用的决策变量,基于所 述决策变量和所述有功调度参数,构建包含机组最小开停机时间精确约束和 灵活爬坡备用约束在内的日前有功计划约束,根据所述日前有功计划约束建 立燃气发电机组日前有功计划模型。
在一种能够实现的方式中,所述引入燃气发电机组关于启停、出力和备用 的决策变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停 机时间精确约束和灵活爬坡备用约束在内的日前有功计划约束,包括:
引入二进制变量
Figure BDA0003687047600000131
Figure BDA0003687047600000132
Figure BDA0003687047600000133
Figure BDA0003687047600000134
分别表示机组g在时段t开机状 态变量实际可取得的最小值和最大值;引入日前的连续变量
Figure BDA0003687047600000135
Figure BDA0003687047600000136
Figure BDA0003687047600000137
Figure BDA0003687047600000138
分别表示机组g在日内时段t出力的基准上限和基准下限;引入日前的连 续变量
Figure BDA0003687047600000139
Figure BDA00036870476000001310
Figure BDA00036870476000001311
Figure BDA00036870476000001312
分别表示机组g在时段t提供的向上灵活爬坡备用 容量和向下灵活爬坡备用容量;
构建机组最小开机时间精确约束:
Figure BDA00036870476000001313
式中,
Figure BDA00036870476000001314
表示机组g在时段t+τ开机状态变量实际可取得的最小值,
Figure BDA00036870476000001315
表示机组g在时段t-1开机状态变量实际可取得的最小值,
Figure BDA00036870476000001316
表示燃 气发电机组的指标集合,T表示计划时段的指标集合,
Figure BDA00036870476000001317
表示机组g的最小 开机时间;
构建机组最小停机时间精确约束:
Figure BDA00036870476000001318
式中,
Figure BDA00036870476000001319
表示机组g在时段t+τ开机状态变量实际可取得的最大值,
Figure BDA00036870476000001320
表示机组g在时段t-1开机状态变量实际可取得的最大值,
Figure BDA00036870476000001321
表示机组 g的最小停机时间;
构建向上灵活爬坡备用约束:
Figure BDA00036870476000001322
式中,
Figure BDA00036870476000001323
表示机组g在日内时段t-1出力的基准下限;
构建向下灵活爬坡备用约束:
Figure BDA00036870476000001324
式中,
Figure BDA0003687047600000141
示机组g在日内时段t-1出力的基准上限。
该向上灵活爬坡备用约束的表达式表明,当
Figure BDA0003687047600000142
为负数时,
Figure BDA0003687047600000143
取 0;当
Figure BDA0003687047600000144
为正数时,
Figure BDA0003687047600000145
Figure BDA0003687047600000146
该向下灵活爬坡备用约束的表达式表明,当
Figure BDA0003687047600000147
为负数时,
Figure BDA0003687047600000148
取 0;当
Figure BDA0003687047600000149
为正数时,
Figure BDA00036870476000001410
Figure BDA00036870476000001411
本发明实施例中,构建机组最小开机时间精确约束和机组最小停机时间 精确约束,是本发明考虑概率化电力平衡模式下燃气发电机组非预期性的关 键技术。本发明进一步创建了向上灵活爬坡备用约束和向下灵活爬坡备用约 束,该灵活爬坡备用约束是已有的机组组合模型中未有的约束,构建活爬坡备 用是本发明考虑概率化电力平衡模式下燃气发电机组灵活性,以及量化评估 燃气发电机组调节贡献的关键技术。通过构建包含机组最小开停机时间精确 约束和灵活爬坡备用约束在内的日前有功计划约束,有利于提高概率化电力 平衡模式下燃气发电机组模型的精确性。
在一种能够实现的方式中,所述引入燃气发电机组关于启停、出力和备用 的决策变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停 机时间精确约束和灵活爬坡备用约束在内的日前有功计划约束,还包括:
引入日前的二进制变量yg,t,当yg,t取1时表示机组g在日内的时段t允许停 机,取0时表示机组g在日内的时段t不允许停机或者因为基准计划已经是停 机状态而无法进行关停;引入日前的连续变量
Figure BDA00036870476000001412
Figure BDA00036870476000001413
表示机组g在日内时 段t出力的实际下限;
构建机组开机状态变量最小值的表达式:
Figure BDA00036870476000001414
式中,xg,t表示机组g在时段t开停机状态的基准计划;
构建机组开机状态变量最大值的表达式:
Figure BDA00036870476000001415
构建机组出力上下限关联关系约束:
Figure BDA0003687047600000151
式中,P g 表示机组g的最小技术出力,
Figure BDA0003687047600000152
表示机组g的装机容量;
构建机组出力上爬坡约束:
Figure BDA0003687047600000153
其中,
Figure BDA0003687047600000154
表示机组g在日内时段t-1出力的实际下限,
Figure BDA0003687047600000155
表示机组g 的最大上爬坡能力;
构建机组出力下爬坡约束:
Figure BDA0003687047600000156
式中,
Figure BDA0003687047600000157
表示机组g的最大下爬坡能力;
构建机组出力实际下限的表达式:
Figure BDA0003687047600000158
构建旋转备用约束:
Figure BDA0003687047600000159
式中,
Figure BDA00036870476000001510
表示机组g在时段t提供的旋转备用容量。
本实施例中,机组开机状态变量最小值的表达式建立了机组g在时段t的 开机状态变量实际可取得的最小值与xg,t、yg,t之间的关联关系,即当机组g在 时段t的机组状态是开机,且允许在日内停机时,开机状态变量
Figure BDA00036870476000001511
的值是0。 其中,变量xg,t是常规机组组合模型中已有的决策变量。
本实施例中,机组开机状态变量最大值的表达式建立了机组g在时段t的 开机状态变量实际可取得的最大值与开停机状态的基准计划变量xg,t的关系; 机组出力上下限关联关系约束建立了机组g在日内时段t出力上下限和机组开 机状态基准计划之间的关系。
本实施例中,机组出力实际下限的表达式建立了机组g在日内时段t的出 力实际下限变量与基准下限变量、是否允许日内停机变量之间的关联关系,即 当允许机组g在日内时段t停机时,出力的实际下限为0,否则出力的实际下限 与基准下限一致。
本发明实施例,在构建机组最小开停机时间精确约束和灵活爬坡备用约 束的基础上,进一步构建了其他日前有功计划约束,完善了燃气发电机日前有 功计划模型,使得该日前有功计划模型能够在计及非预期性的基础上更加贴 合实际应用需求。
步骤S3,根据所述历史数据确定计划日的新能源出力预测误差集合,将 所述新能源出力预测误差集合的高维凸包作为凸包络场景集,根据所述有功 调度参数构建考虑所述凸包络场景集的时序解耦日内有功调度模型。
在一种能够实现的方式中,所述根据所述历史数据确定计划日的新能源 出力预测误差集合,包括:
在新能源的所有历史出力预测值中确定与计划日新能源出力预测值相似 度最大的数据子集;
将所述数据子集对应的历史出力真实值作为计划日可能的出力真实值, 计算计划日可能的出力真实值与对应的计划日出力预测值的差值;
将所述差值作为新能源出力预测误差,构建计划日的新能源出力预测误 差集合。
其中,可以以曼哈顿距离为度量,确定历史出力预测值与计划日新能源出 力预测值的相似度。
需要说明的是,还可以以其他合适的距离参数例如欧式距离等为度量,来 确定历史出力预测值与计划日新能源出力预测值的相似度。本发明实施例不 限定于此。
本发明实施例,在所有历史出力预测值中找到与计划日新能源出力预测 值较为接近的一个数据子集,进而将该数据子集对应的历史真实值作为计划 日可能的出力真实值,通过求差值的方式,得到计划日的新能源出力预测误差 集合,最后将计划日的新能源出力预测误差集合的高维凸包作为所述的凸包 络场景集。作为示例,构建的凸包络场景集如图2所示。
本发明实施例中,计划日的新能源出力预测误差集合的构建方法简单便 捷,且贴合新能源的实际出力情况。
在一种能够实现的方式中,所述根据所述有功调度参数构建考虑所述凸 包络场景集的时序解耦日内有功调度模型,包括:
引入日内实时运行的二进制变量zg,t,当zg,t取1时表示机组g在时段t进行 关停,取0时表示机组g在时段t保持开机运行状态;引入新能源出力预测误差 随机向量ξ;引入连续变量pg,t(ξ),pg,t(ξ)表示在场景ξ发生的情况下,机组g在 时段t的实际出力;
构建机组停机约束:
Figure BDA0003687047600000171
构建机组出力上下限约束:
Figure BDA0003687047600000172
本发明实施例中,该机组停机约束的不等式表明,仅当机组g的时段t在日 前计划中被设定为允许停机时,日内实时运行才可以停机。以上两条不等式用 于描述机组出力pg,t(ξ)与日前计划变量
Figure BDA0003687047600000173
及日内调控变量zg,t之间的关 联关系。当zg,t取1,机组g在时段t进行关停时,pg,t(ξ)必须为0;当zg,t取0, 机组g在时段t保持开机,pg,t(ξ)可以在
Figure BDA0003687047600000174
Figure BDA0003687047600000175
的范围内任意调控。
步骤S4,根据新能源出力削减约束、功率平衡约束和线路潮流约束建立 系统侧模型,联合所述燃气发电机组日前有功计划模型、所述时序解耦日内有 功调度模型和所述系统侧模型,以计划成本加新能源出力预测场景下的运行 成本之和最小为优化目标,构建燃气发电机组启停、出力和备用的联合优化模 型。
在一种能够实现的方式中,所述根据新能源出力削减约束、功率平衡约束 和线路潮流约束建立系统侧模型,包括:
构建新能源出力削减约束为:
Figure BDA0003687047600000176
式中,
Figure BDA0003687047600000177
表示在新能源出力场景ξ下系统需要削减的新能源功率, Ww,t表示新能源场站w在时段t的预测出力,ξw,t表示新能源场站w在时段t的实 际出力预测误差,
Figure BDA0003687047600000178
表示新能源场站的指标集合;
构建功率平衡约束为:
Figure BDA0003687047600000181
式中,Dd,t表示负荷节点d在时段t的负荷需求,
Figure BDA0003687047600000182
表示负荷节点的指标集 合;
构建线路潮流约束为:
Figure BDA0003687047600000183
式中,Fl表示线路l的潮流限值,
Figure BDA0003687047600000184
表示机组到线路l的功率转移因子,
Figure BDA0003687047600000185
表示负荷到线路l的功率转移因子,
Figure BDA0003687047600000186
表示新能源场站到线路l的功率转移 因子,
Figure BDA0003687047600000187
表示输电线路的指标集合。
步骤S5,采用场景生成法求解所述联合优化模型,根据得到的求解结果 生成燃气发电机组的启停、出力和备用计划。
在一种能够实现的方式中,所述采用场景生成法求解所述联合优化模型, 包括:
S701,考虑预测场景,形成一个松弛模型,调用最优化求解器获得一个待 校验的日前计划方案;
S702,给定日前计划方案,在新能源出力预测值的基础上叠加预测误差场 景集中的误差值,形成多个新能源出力场景,并调用最优化求解器逐个验证日 内实时运行是否存在可行方案;
S703,如果存在不可行方案,则将产生不可行方案的新能源出力场景扩充 到所述松弛模型中,更新得到新的待校验日前计划方案;
S704,重复步骤S702-S703,直至日前计划方案通过校验。
作为优选,可采用gurobi求解器作为该最优化求解器。
作为示例,采用场景生成法求解所述联合优化模型,得到的求解结果如表 1所示。
表1:
Figure BDA0003687047600000191
本发明实施例中,采用场景生成法求解日前-日内两阶段的燃气发电机组 启停、出力和备用联合优化模型,能够有效提高计算效率。
本发明还提供了一种燃气发电机组启停-出力-备用联合优化装置,该装置 可以用于实施本发明上述任一项实施例所述的燃气发电机组启停-出力-备用 联合优化方法。
请参阅图3,图3示出了本发明实施例提供的一种燃气发电机组启停-出 力-备用联合优化装置的原理框图。
本发明实施例提供了一种燃气发电机组启停-出力-备用联合优化装置,包 括:
获取模块1,用于获取电力系统的有功调度参数以及新能源关于出力预测 值和出力真实值的历史数据;
第一模型构建模块2,用于引入燃气发电机组关于启停、出力和备用的决 策变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停机时 间精确约束和灵活爬坡备用约束在内的日前有功计划约束,根据所述日前有 功计划约束建立燃气发电机组日前有功计划模型;
第二模型构建模块3,用于根据所述历史数据确定计划日的新能源出力预 测误差集合,将所述新能源出力预测误差集合的高维凸包作为凸包络场景集, 根据所述有功调度参数构建考虑所述凸包络场景集的时序解耦日内有功调度 模型;
模型联合模块4,用于根据新能源出力削减约束、功率平衡约束和线路潮 流约束建立系统侧模型,联合所述燃气发电机组日前有功计划模型、所述时序 解耦日内有功调度模型和所述系统侧模型,以计划成本加新能源出力预测场 景下的运行成本之和最小为优化目标,构建燃气发电机组启停、出力和备用的 联合优化模型;
优化求解模块5,用于采用场景生成法求解所述联合优化模型,根据得到 的求解结果生成燃气发电机组的启停、出力和备用计划。
在一种能够实现的方式中,所述第二模型构建模块3包括:
确定单元,用于在新能源的所有历史出力预测值中确定与计划日新能源 出力预测值相似度最大的数据子集;
计算单元,用于将所述数据子集对应的历史出力真实值作为计划日可能 的出力真实值,计算计划日可能的出力真实值与对应的计划日出力预测值的 差值;
构建单元,用于将所述差值作为新能源出力预测误差,构建计划日的新能 源出力预测误差集合。
在一种能够实现的方式中,所述第一模型构建模块2包括:
第一变量引入单元,用于引入二进制变量
Figure BDA0003687047600000201
Figure BDA0003687047600000202
Figure BDA0003687047600000203
Figure BDA0003687047600000204
分别表 示机组g在时段t开机状态变量实际可取得的最小值和最大值;引入日前的连 续变量
Figure BDA0003687047600000205
Figure BDA0003687047600000206
Figure BDA0003687047600000207
Figure BDA0003687047600000208
分别表示机组g在日内时段t出力的基准上限和基 准下限;引入日前的连续变量
Figure BDA0003687047600000209
Figure BDA00036870476000002010
Figure BDA00036870476000002011
Figure BDA00036870476000002012
分别表示机组g在时段t提供 的向上灵活爬坡备用容量和向下灵活爬坡备用容量;
第一构建单元,用于构建机组最小开机时间精确约束:
Figure BDA0003687047600000211
式中,
Figure BDA0003687047600000212
表示机组g在时段t+τ开机状态变量实际可取得的最小值,
Figure BDA0003687047600000213
表示机组g在时段t-1开机状态变量实际可取得的最小值,
Figure BDA0003687047600000214
表示燃 气发电机组的指标集合,T表示计划时段的指标集合,
Figure BDA0003687047600000215
表示机组g的最小 开机时间;
第二构建单元,用于构建机组最小停机时间精确约束:
Figure BDA0003687047600000216
式中,
Figure BDA0003687047600000217
表示机组g在时段t+τ开机状态变量实际可取得的最大值,
Figure BDA0003687047600000218
表示机组g在时段t-1开机状态变量实际可取得的最大值,
Figure BDA0003687047600000219
表示机组 g的最小停机时间;
第三构建单元,用于构建向上灵活爬坡备用约束:
Figure BDA00036870476000002110
式中,
Figure BDA00036870476000002111
表示机组g在日内时段t-1出力的基准下限;
第四构建单元,用于构建向下灵活爬坡备用约束:
Figure BDA00036870476000002112
式中,
Figure BDA00036870476000002113
示机组g在日内时段t-1出力的基准上限。
在一种能够实现的方式中,所述第一模型构建模块2还包括:
第二变量引入单元,用于引入日前的二进制变量yg,t,当yg,t取1时表示机 组g在日内的时段t允许停机,取0时表示机组g在日内的时段t不允许停机或者 因为基准计划已经是停机状态而无法进行关停;引入日前的连续变量
Figure BDA00036870476000002114
Figure BDA00036870476000002115
表示机组g在日内时段t出力的实际下限;
第五构建单元,用于构建机组开机状态变量最小值的表达式:
Figure BDA00036870476000002116
式中,xg,t表示机组g在时段t开停机状态的基准计划;
第六构建单元,用于构建机组开机状态变量最大值的表达式:
Figure BDA0003687047600000221
第七构建单元,用于构建机组出力上下限关联关系约束:
Figure BDA0003687047600000222
式中,P g 表示机组g的最小技术出力,
Figure BDA0003687047600000223
表示机组g的装机容量;
第八构建单元,用于构建机组出力上爬坡约束:
Figure BDA0003687047600000224
其中,
Figure BDA0003687047600000225
表示机组g在日内时段t-1出力的实际下限,
Figure BDA0003687047600000226
表示机组g 的最大上爬坡能力;
第九构建单元,用于构建机组出力下爬坡约束:
Figure BDA0003687047600000227
式中,
Figure BDA0003687047600000228
表示机组g的最大下爬坡能力;
第十构建单元,用于构建机组出力实际下限的表达式:
Figure BDA0003687047600000229
第十一构建单元,用于构建旋转备用约束:
Figure BDA00036870476000002210
式中,
Figure BDA00036870476000002211
表示机组g在时段t提供的旋转备用容量。
在一种能够实现的方式中,所述第二模型构建模块3包括:
第三变量引入单元,用于引入日内实时运行的二进制变量zg,t,当zg,t取1 时表示机组g在时段t进行关停,取0时表示机组g在时段t保持开机运行状态; 引入新能源出力预测误差随机向量ξ;引入连续变量pg,t(ξ),pg,t(ξ)表示在场 景ξ发生的情况下,机组g在时段t的实际出力;
第十二构建单元,用于构建机组停机约束:
Figure BDA00036870476000002212
第十三构建单元,用于构建机组出力上下限约束:
Figure BDA0003687047600000231
在一种能够实现的方式中,所述模型联合模块4包括:
第十四构建单元,用于构建新能源出力削减约束为:
Figure BDA0003687047600000232
式中,
Figure BDA0003687047600000233
表示在新能源出力场景ξ下系统需要削减的新能源功率, Ww,t表示新能源场站w在时段t的预测出力,ξw,t表示新能源场站w在时段t的实 际出力预测误差,
Figure BDA00036870476000002311
表示新能源场站的指标集合;
第十五构建单元,用于构建功率平衡约束为:
Figure BDA0003687047600000234
式中,Dd,t表示负荷节点d在时段t的负荷需求,
Figure BDA0003687047600000235
表示负荷节点的指标集 合;
第十六构建单元,用于构建线路潮流约束为:
Figure BDA0003687047600000236
式中,Fl表示线路l的潮流限值,
Figure BDA0003687047600000237
表示机组到线路l的功率转移因子,
Figure BDA0003687047600000238
表示负荷到线路l的功率转移因子,
Figure BDA0003687047600000239
表示新能源场站到线路l的功率转移 因子,
Figure BDA00036870476000002310
表示输电线路的指标集合。
在一种能够实现的方式中,所述优化求解模块5具体用于:
S701,考虑预测场景,形成一个松弛模型,调用最优化求解器获得一个待 校验的日前计划方案;
S702,给定日前计划方案,在新能源出力预测值的基础上叠加预测误差场 景集中的误差值,形成多个新能源出力场景,并调用最优化求解器逐个验证日 内实时运行是否存在可行方案;
S703,如果存在不可行方案,则将产生不可行方案的新能源出力场景扩充 到所述松弛模型中,更新得到新的待校验日前计划方案;
S704,重复步骤S702-S703,直至日前计划方案通过校验。
本发明还提供了一种燃气发电机组启停-出力-备用联合优化装置,包括:
存储器,用于存储指令;其中,所述指令用于实现如上任意一项实施例所 述的燃气发电机组启停-出力-备用联合优化方法;
处理器,用于执行所述存储器中的指令。
本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质上 存储有计算机程序,所述计算机程序被处理器执行时实现如上任意一项实施 例所述的燃气发电机组启停-出力-备用联合优化方法。
本发明上述实施例,能够适应概率化电力平衡模式,实现燃气发电机组日 前计划的科学和精准编制,解决目前燃气发电机组的日前发电和备用计划编 制未考虑到随机场景下燃气发电机组调控的非预期性,计划编制的精确性有 待提高的技术问题,具体的优点包括:
1)本发明采用历史数据构建新能源出力预测误差的凸包络场景集,作为 概率化电力平衡的基础,有助于高效应对新能源的预测误差;
2)在日前计划模型中引入了燃气发电机组是否允许日内停机的决策变量 和最小开停机时间的精确约束,实现日内调度模型的时序解耦,进而有效计及 了调控的非预期性,保证概率化电力平衡模式下燃气发电机组模型的精确性;
3)利用机组出力上下限基准决策变量构造向上和向下灵活爬坡备用容量, 为概率化电力平衡模式下燃气发电机组调节贡献的评估提供了基础;
4)采用场景生成法求解日前-日内两阶段的燃气发电机组启停、出力和备 用联合优化模型,有效提高了计算效率。
所属领域的技术人员可以清楚地了解到,为描述地方便和简洁,上述描述 的系统、装置和模块的具体工作过程,可以参考前述方法实施例中的对应过程, 上述描述的系统、装置和模块的具体有益效果,可以参考前述方法实施例中的 对应有益效果,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统、装置和方 法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性 的,例如,所述模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另 外的划分方式,例如多个模块或组件可以结合或者可以集成到另一个系统,或 一些特征可以忽略,或不执行。
所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作 为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或 者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者 全部模块来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能模块可以集成在一个处理模块中, 也可以是各个模块单独物理存在,也可以两个或两个以上模块集成在一个模 块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块 的形式实现。
所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售 或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发 明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全 部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存 储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务 器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而 前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、 随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可 以存储程序代码的介质。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽 管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当 理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部 分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质 脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种燃气发电机组启停-出力-备用联合优化方法,其特征在于,包括:
获取电力系统的有功调度参数以及新能源关于出力预测值和出力真实值的历史数据;
引入燃气发电机组关于启停、出力和备用的决策变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停机时间精确约束和灵活爬坡备用约束在内的日前有功计划约束,根据所述日前有功计划约束建立燃气发电机组日前有功计划模型;
根据所述历史数据确定计划日的新能源出力预测误差集合,将所述新能源出力预测误差集合的高维凸包作为凸包络场景集,根据所述有功调度参数构建考虑所述凸包络场景集的时序解耦日内有功调度模型;
根据新能源出力削减约束、功率平衡约束和线路潮流约束建立系统侧模型,联合所述燃气发电机组日前有功计划模型、所述时序解耦日内有功调度模型和所述系统侧模型,以计划成本加新能源出力预测场景下的运行成本之和最小为优化目标,构建燃气发电机组启停、出力和备用的联合优化模型;
采用场景生成法求解所述联合优化模型,根据得到的求解结果生成燃气发电机组的启停、出力和备用计划。
2.根据权利要求1所述的燃气发电机组启停-出力-备用联合优化方法,其特征在于,所述根据所述历史数据确定计划日的新能源出力预测误差集合,包括:
在新能源的所有历史出力预测值中确定与计划日新能源出力预测值相似度最大的数据子集;
将所述数据子集对应的历史出力真实值作为计划日可能的出力真实值,计算计划日可能的出力真实值与对应的计划日出力预测值的差值;
将所述差值作为新能源出力预测误差,构建计划日的新能源出力预测误差集合。
3.根据权利要求1所述的燃气发电机组启停-出力-备用联合优化方法,其特征在于,所述引入燃气发电机组关于启停、出力和备用的决策变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停机时间精确约束和灵活爬坡备用约束在内的日前有功计划约束,包括:
引入二进制变量
Figure FDA0003687047590000021
Figure FDA0003687047590000022
Figure FDA0003687047590000023
Figure FDA0003687047590000024
分别表示机组g在时段t开机状态变量实际可取得的最小值和最大值;引入日前的连续变量
Figure FDA0003687047590000025
Figure FDA0003687047590000026
Figure FDA0003687047590000027
Figure FDA0003687047590000028
分别表示机组g在日内时段t出力的基准上限和基准下限;引入日前的连续变量
Figure FDA0003687047590000029
Figure FDA00036870475900000210
Figure FDA00036870475900000211
Figure FDA00036870475900000212
分别表示机组g在时段t提供的向上灵活爬坡备用容量和向下灵活爬坡备用容量;
构建机组最小开机时间精确约束:
Figure FDA00036870475900000213
式中,
Figure FDA00036870475900000214
表示机组g在时段t+τ开机状态变量实际可取得的最小值,
Figure FDA00036870475900000215
表示机组g在时段t-1开机状态变量实际可取得的最小值,
Figure FDA00036870475900000216
表示燃气发电机组的指标集合,T表示计划时段的指标集合,
Figure FDA00036870475900000217
表示机组g的最小开机时间;
构建机组最小停机时间精确约束:
Figure FDA00036870475900000218
式中,
Figure FDA00036870475900000219
表示机组g在时段t+τ开机状态变量实际可取得的最大值,
Figure FDA00036870475900000220
表示机组g在时段t-1开机状态变量实际可取得的最大值,
Figure FDA00036870475900000221
表示机组g的最小停机时间;
构建向上灵活爬坡备用约束:
Figure FDA00036870475900000222
式中,
Figure FDA00036870475900000223
表示机组g在日内时段t-1出力的基准下限;
构建向下灵活爬坡备用约束:
Figure FDA00036870475900000224
式中,
Figure FDA00036870475900000225
示机组g在日内时段t-1出力的基准上限。
4.根据权利要求3所述的燃气发电机组启停-出力-备用联合优化方法,其特征在于,所述引入燃气发电机组关于启停、出力和备用的决策变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停机时间精确约束和灵活爬坡备用约束在内的日前有功计划约束,还包括:
引入日前的二进制变量yg,t,当yg,t取1时表示机组g在日内的时段t允许停机,取0时表示机组g在日内的时段t不允许停机或者因为基准计划已经是停机状态而无法进行关停;引入日前的连续变量
Figure FDA0003687047590000031
Figure FDA0003687047590000032
表示机组g在日内时段t出力的实际下限;
构建机组开机状态变量最小值的表达式:
Figure FDA0003687047590000033
式中,xg,t表示机组g在时段t开停机状态的基准计划;
构建机组开机状态变量最大值的表达式:
Figure FDA0003687047590000034
构建机组出力上下限关联关系约束:
Figure FDA0003687047590000035
式中,P g 表示机组g的最小技术出力,
Figure FDA0003687047590000036
表示机组g的装机容量;
构建机组出力上爬坡约束:
Figure FDA0003687047590000037
其中,
Figure FDA0003687047590000038
表示机组g在日内时段t-1出力的实际下限,
Figure FDA0003687047590000039
表示机组g的最大上爬坡能力;
构建机组出力下爬坡约束:
Figure FDA00036870475900000310
式中,
Figure FDA00036870475900000311
表示机组g的最大下爬坡能力;
构建机组出力实际下限的表达式:
Figure FDA00036870475900000312
构建旋转备用约束:
Figure FDA00036870475900000313
式中,
Figure FDA00036870475900000314
表示机组g在时段t提供的旋转备用容量。
5.根据权利要求4所述的燃气发电机组启停-出力-备用联合优化方法,其特征在于,所述根据所述有功调度参数构建考虑所述凸包络场景集的时序解耦日内有功调度模型,包括:
引入日内实时运行的二进制变量zg,t,当zg,t取1时表示机组g在时段t进行关停,取0时表示机组g在时段t保持开机运行状态;引入新能源出力预测误差随机向量ξ;引入连续变量pg,t(ξ),pg,t(ξ)表示在场景ξ发生的情况下,机组g在时段t的实际出力;
构建机组停机约束:
Figure FDA0003687047590000041
构建机组出力上下限约束:
Figure FDA0003687047590000042
6.根据权利要求5所述的燃气发电机组启停-出力-备用联合优化方法,其特征在于,所述根据新能源出力削减约束、功率平衡约束和线路潮流约束建立系统侧模型,包括:
构建新能源出力削减约束为:
Figure FDA0003687047590000043
式中,
Figure FDA0003687047590000044
表示在新能源出力场景ξ下系统需要削减的新能源功率,Ww,t表示新能源场站w在时段t的预测出力,ξw,t表示新能源场站w在时段t的实际出力预测误差,
Figure FDA0003687047590000045
表示新能源场站的指标集合;
构建功率平衡约束为:
Figure FDA0003687047590000046
式中,Dd,t表示负荷节点d在时段t的负荷需求,
Figure FDA0003687047590000047
表示负荷节点的指标集合;
构建线路潮流约束为:
Figure FDA0003687047590000051
式中,Fl表示线路l的潮流限值,
Figure FDA0003687047590000052
表示机组到线路l的功率转移因子,
Figure FDA0003687047590000053
表示负荷到线路l的功率转移因子,
Figure FDA0003687047590000054
表示新能源场站到线路l的功率转移因子,
Figure FDA0003687047590000055
表示输电线路的指标集合。
7.根据权利要求1所述的燃气发电机组启停-出力-备用联合优化方法,其特征在于,所述采用场景生成法求解所述联合优化模型,包括:
S701,考虑预测场景,形成一个松弛模型,调用最优化求解器获得一个待校验的日前计划方案;
S702,给定日前计划方案,在新能源出力预测值的基础上叠加预测误差场景集中的误差值,形成多个新能源出力场景,并调用最优化求解器逐个验证日内实时运行是否存在可行方案;
S703,如果存在不可行方案,则将产生不可行方案的新能源出力场景扩充到所述松弛模型中,更新得到新的待校验日前计划方案;
S704,重复步骤S702-S703,直至日前计划方案通过校验。
8.一种燃气发电机组启停-出力-备用联合优化装置,其特征在于,包括:
获取模块,用于获取电力系统的有功调度参数以及新能源关于出力预测值和出力真实值的历史数据;
第一模型构建模块,用于引入燃气发电机组关于启停、出力和备用的决策变量,基于所述决策变量和所述有功调度参数,构建包含机组最小开停机时间精确约束和灵活爬坡备用约束在内的日前有功计划约束,根据所述日前有功计划约束建立燃气发电机组日前有功计划模型;
第二模型构建模块,用于根据所述历史数据确定计划日的新能源出力预测误差集合,将所述新能源出力预测误差集合的高维凸包作为凸包络场景集,根据所述有功调度参数构建考虑所述凸包络场景集的时序解耦日内有功调度模型;
模型联合模块,用于根据新能源出力削减约束、功率平衡约束和线路潮流约束建立系统侧模型,联合所述燃气发电机组日前有功计划模型、所述时序解耦日内有功调度模型和所述系统侧模型,以计划成本加新能源出力预测场景下的运行成本之和最小为优化目标,构建燃气发电机组启停、出力和备用的联合优化模型;
优化求解模块,用于采用场景生成法求解所述联合优化模型,根据得到的求解结果生成燃气发电机组的启停、出力和备用计划。
9.一种燃气发电机组启停-出力-备用联合优化装置,其特征在于,包括:
存储器,用于存储指令;其中,所述指令用于实现如权利要求1-7任意一项所述的燃气发电机组启停-出力-备用联合优化方法;
处理器,用于执行所述存储器中的指令。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1-7任意一项所述的燃气发电机组启停-出力-备用联合优化方法。
CN202210648680.2A 2022-06-09 2022-06-09 一种燃气发电机组启停-出力-备用联合优化方法及装置 Active CN115296345B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210648680.2A CN115296345B (zh) 2022-06-09 2022-06-09 一种燃气发电机组启停-出力-备用联合优化方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210648680.2A CN115296345B (zh) 2022-06-09 2022-06-09 一种燃气发电机组启停-出力-备用联合优化方法及装置

Publications (2)

Publication Number Publication Date
CN115296345A true CN115296345A (zh) 2022-11-04
CN115296345B CN115296345B (zh) 2023-08-22

Family

ID=83820400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210648680.2A Active CN115296345B (zh) 2022-06-09 2022-06-09 一种燃气发电机组启停-出力-备用联合优化方法及装置

Country Status (1)

Country Link
CN (1) CN115296345B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116526468A (zh) * 2023-05-09 2023-08-01 国网湖北省电力有限公司经济技术研究院 考虑多类型备用辅助服务的高比例可再生能源电力系统优化调度方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006180673A (ja) * 2004-12-24 2006-07-06 Mitsubishi Electric Corp 発電機の運用計画作成装置
CN109242366A (zh) * 2018-11-06 2019-01-18 国电南瑞科技股份有限公司 一种电-气互联综合能源系统的多时段潮流优化方法
CN111525625A (zh) * 2020-04-01 2020-08-11 国网天津市电力公司 一种考虑燃气机组的日前计划制定方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006180673A (ja) * 2004-12-24 2006-07-06 Mitsubishi Electric Corp 発電機の運用計画作成装置
CN109242366A (zh) * 2018-11-06 2019-01-18 国电南瑞科技股份有限公司 一种电-气互联综合能源系统的多时段潮流优化方法
CN111525625A (zh) * 2020-04-01 2020-08-11 国网天津市电力公司 一种考虑燃气机组的日前计划制定方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张旭、王洪涛: "高比例可再生能源电力系统的输配协同优化调度方法", 《电力系统自动化》, vol. 43, no. 3, pages 67 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116526468A (zh) * 2023-05-09 2023-08-01 国网湖北省电力有限公司经济技术研究院 考虑多类型备用辅助服务的高比例可再生能源电力系统优化调度方法
CN116526468B (zh) * 2023-05-09 2024-04-26 国网湖北省电力有限公司经济技术研究院 考虑多类型备用辅助服务的高比例可再生能源电力系统优化调度方法

Also Published As

Publication number Publication date
CN115296345B (zh) 2023-08-22

Similar Documents

Publication Publication Date Title
Li et al. A coordinated dispatch method with pumped-storage and battery-storage for compensating the variation of wind power
Lima et al. On the computational studies of deterministic global optimization of head dependent short-term hydro scheduling
CN111555281B (zh) 一种电力系统灵活性资源配置仿真方法及装置
CN110854932B (zh) 一种交直流配电网多时间尺度优化调度方法及系统
Tejada-Arango et al. Representation of storage operations in network-constrained optimization models for medium-and long-term operation
JP4187620B2 (ja) 発電機の起動停止計画作成方法とその装置及び起動停止計画作成装置の処理プログラムを記録する記録媒体。
CN111092455A (zh) 一种储能系统与已恢复机组联合运行的负荷恢复优化方法
CN110854931A (zh) 一种抽水蓄能机组日前发电计划编制方法、系统以及设备
CN112510690B (zh) 考虑风火储联合和需求响应奖惩的优化调度方法及系统
CN115296345A (zh) 一种燃气发电机组启停-出力-备用联合优化方法及装置
CN114784882A (zh) 机组组合优化处理方法及装置
CN113364043A (zh) 一种基于条件风险价值的微电网群优化方法
CN109787217B (zh) 基于风电多状态模型及机会成本修正的备用出清方法
Mahlke A scenario tree-based decomposition for solving multistage stochastic programs: with application in energy production
CN115907140A (zh) 电力现货出清方案优化方法、装置、计算机设备和介质
CN115549137A (zh) 分布式电网调控系统及调控方法
CN114094609A (zh) 电网储能系统的优化配置方法
Fu et al. Study on wind power consumption capacity of power grid considering risk and unit commitment
CN113489066A (zh) 计及供需不确定性的含储能电网区间供电可靠性评估方法
CN112054553A (zh) 电热气互联系统协调优化运行方法、系统、介质及设备
Yang et al. A Multistage Semianticipativity Model for UC via Distributionally Robust Optimization Combined With Principal Component Analysis
CN118054480A (zh) 一种综合能源电力系统优化调度方法、装置、设备及介质
CN116128193A (zh) 一种基于区块链技术的园区微网调度方法及调度系统
Collins Modelling the interactions between power systems and energy systems
Yang et al. A Fully Adaptive DRO Multistage Framework Based on MDR for Generation Scheduling under Uncertainty

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant