CN109861289B - 一种基于含储能单元的电力系统可用传输容量的控制方法 - Google Patents

一种基于含储能单元的电力系统可用传输容量的控制方法 Download PDF

Info

Publication number
CN109861289B
CN109861289B CN201910190005.8A CN201910190005A CN109861289B CN 109861289 B CN109861289 B CN 109861289B CN 201910190005 A CN201910190005 A CN 201910190005A CN 109861289 B CN109861289 B CN 109861289B
Authority
CN
China
Prior art keywords
power
power system
energy storage
storage unit
generating unit
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.)
Active
Application number
CN201910190005.8A
Other languages
English (en)
Other versions
CN109861289A (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.)
Jiangsu Qingke Electric Co ltd
Tsinghua University
Original Assignee
Jiangsu Qingke Electric Co ltd
Tsinghua University
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 Jiangsu Qingke Electric Co ltd, Tsinghua University filed Critical Jiangsu Qingke Electric Co ltd
Priority to CN201910190005.8A priority Critical patent/CN109861289B/zh
Publication of CN109861289A publication Critical patent/CN109861289A/zh
Application granted granted Critical
Publication of CN109861289B publication Critical patent/CN109861289B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明涉及一种基于含储能单元的电力系统可用传输容量的控制方法,属于电力系统规划与运行控制技术领域。本发明方法首先获取电力系统中储能单元、发电机组、输电线路的数据,建立两阶段的电力系统可用传输容量计算模型,进而对模型各阶段依次迭代求解。第一阶段建立电力系统日前市场模型,求解得到系统中各发电机组的启停计划和电能外送计划;第二阶段建立电力市场实时市场模型,迭代求解每小时发电机组出力与储能单元充放电功率,最终得到电力系统中发电机组、储能单元的调度运行方案。本发明通过迭代计算确定发电机组与储能单元的运行方案,使储能单元的调节作用得到充分利用,进而精准控制与提升电力系统可用传输容量。

Description

一种基于含储能单元的电力系统可用传输容量的控制方法
技术领域
本发明涉及一种基于含储能单元的电力系统可用传输容量的控制方法,属于电力系统规划与运行控制技术领域。
背景技术
电力系统的可用传输容量(Available Transfer Capability,ATC),是指在既有输电计划的基础上,输电线路所剩余的电能传输能力,是电力市场中保障电力自由可靠交易的一个重要因素。随着化石燃料短缺和环境污染等问题的日益严重,可再生能源发电机组正大量接入电力系统。由于风电、太阳能等可再生能源与负荷中心的地域分布不同,电能跨区交易的需求不断增加,对电力系统的可用传输容量提出了更高的要求,例如我国已建成多条高压直流输电线路以满足西南地区水电外送、西北地区风电外送等输电需求。
传统的电力系统可用传输容量的计算方法,通常只考虑固定负荷与发电机组的影响。随着储能与分布式电源等技术的高速发展,储能单元等需求侧灵活响应资源的调节能力不断提升,采用传统方法计算电力系统可用传输容量将出现较大偏差,影响电网规划,增加输电成本。因此,充分利用电力系统中储能单元的调节能力,实现发电机组与储能单元的协同运行,才能精准控制和提升电力系统的可用传输容量。
发明内容
本发明的目的是提出一种基于含储能单元的电力系统可用传输容量的控制方法,通过建立两阶段的电力系统数学模型和迭代计算方法,确定发电机组与储能单元的协同运行方案,精准控制和提升电力系统可用传输容量,用于优化发电机组与储能单元的运行。
本发明提出的基于含储能单元的电力系统可用传输容量的控制方法,包括以下步骤:
(1)从电力系统控制中心获取电力系统中储能单元、发电机组及输电线路的参数,储能单元参数分别包括:最大充电功率、最大放电功率、容量、充电效率、放电效率、最大充电状态和最小充电状态,发电机组参数分别包括:所有发电机组的装机容量、火电机组的煤耗成本、火电机组的最小运行时间和火电机组的最小冷却时间,输电线路参数包括:电力系统的输电线路传输容量,电能外送线路的最大传输功率,电能外送线路的最小传输容量和节点负荷;
(2)建立一个电力系统日前市场模型,并求解,得到各发电机组的启停计划,具体过程如下:
(2-1)建立电力系统日前市场模型的目标函数,表达式如下:
Figure BDA0001994105460000021
其中,
Figure BDA0001994105460000022
表示火电机组运行与启停的总成本;
Figure BDA0001994105460000023
表示电能外送的总收益,Φ表示集合,上标S、T、G、I分别表示电力系统运行场景、控制时段、火电机组和电力系统外送输电线路,ΦT,Φc和ΦI分别从电力系统控制中心获取,γs为电力系统运行场景s的权重,ΦS和γs从电力系统控制中心获取,
Figure BDA0001994105460000024
为第g个火电机组在第s个电力系统运行场景下第t控制时段的输出功率,Yg,t为第g个火电机组在第t控制时段的启动状态值,若火电机组在第t-1控制时段为关闭,在第t控制时段为运行,则Yg,t为1,若火电机组在第t-1控制时段为运行,在第t控制时段为关闭,则Yg,t为0,若火电机组在第t-1控制时段为运行,在第t控制时段仍为运行,则Yg,t为0,若火电机组在第t-1控制时段为关闭,在第t控制时段仍为关闭,则Yg,t为0;Zg,t为第g个火电机组在第t控制时段的关机状态值,若火电机组在第t-1控制时段为运行,在第t控制时段为关闭,则Zg,t为1,若火电机组在第t-1控制时段为关闭,在第t控制时段为运行,则Zg,t为0,若火电机组在第t-1控制时段为运行,在第t控制时段仍为运行,则Zg,t为0,若火电机组在第t-1控制时段为关闭,在第t控制时段仍为关闭,则Zg,t为0;
Figure BDA0001994105460000025
为第g个火电机组在第s个电力系统运行场景下第t控制时段的运行成本,
Figure BDA0001994105460000026
为第g个火电机组在第t控制时段的启动成本,
Figure BDA0001994105460000027
为第g个火电机组在第t控制时段的关机成本,
Figure BDA0001994105460000028
为电力系统第i条电能外送线路在第t控制时段的传输功率,
Figure BDA0001994105460000029
为第i条电能外送线路在第t控制时段电能外送的单位收益;
(2-2)建立电力系统日前市场模型的约束条件,具体如下:
(2-2-1)电力系统火电机组的功率约束如下:
Figure BDA00019941054600000210
其中,
Figure BDA0001994105460000031
Figure BDA0001994105460000032
分别为设定的第g个火电机组的最小功率限制和最大功率限制,Ug,t为第g个火电机组在第t控制时段的运行状态值,若火电机组在第t控制时段运行,则Ug,t为1,若火电机组在第t控制时段关机,则Ug,t为0;
(2-2-2)电力系统火电机组工作状态的逻辑约束,即在第t控制时段火电机组g只能处于运行或关机状态中的一种,表达式如下:
Figure BDA0001994105460000033
Figure BDA0001994105460000034
(2-2-3)电力系统火电机组的启动和关机时间约束为:
Figure BDA0001994105460000035
Figure BDA0001994105460000036
其中,
Figure BDA0001994105460000037
表示第g个火电机组连续开机的最少时间,
Figure BDA0001994105460000038
表示第g个火电机组保持关机状态最少时间;
(2-2-4)电力系统中可再生能源发电机组的功率约束,表达式为:
Figure BDA0001994105460000039
其中,上标R表示可再生能源发电机组,
Figure BDA00019941054600000310
为第r个可再生能源发电机组第s个电力系统运行场景下在第t控制时段的最大输出功率的预测值,
Figure BDA00019941054600000311
为第r个可再生能源发电机组第s个电力系统运行场景下在第t控制时段的输出功率;
(2-2-5)电力系统可再生能源的消纳率约束,表达式为:
Figure BDA00019941054600000312
其中,RA为设定的可再生能源消纳率的约束下限,ΦR为可再生能源的发电机组集合;
(2-2-6)电力系统有功功率平衡约束为:
Figure BDA00019941054600000313
其中,上标D表示电力系统节点负荷,上标CHA表示电力系统储能单元充电,上标DIS表示储能单元放电,上标ES表示储能单元,ΦD和ΦES分别为节点负荷集合和储能单元集合,
Figure BDA0001994105460000041
为电力系统中节点b处第s个场景下第t控制时段的负荷,
Figure BDA0001994105460000042
表示第j个储能单元在第s个电力系统运行场景下第t个控制时段的充电功率,
Figure BDA0001994105460000043
表示第j个储能单元在第s个电力系统运行场景下第t个控制时段的放电功率;
(2-2-7)电力系统直流潮流约束,,表达式为:
Figure BDA0001994105460000044
其中,L表示电力系统中的输电线路,
Figure BDA0001994105460000045
为电力系统中第l条输电线路在第s个电力系统运行场景下第t控制时段的传输功率,Fl-b表示第l条输电线路对第b个节点的发电转移分布因子,Fl-b从电力系统控制中心获取,
Figure BDA0001994105460000046
Figure BDA0001994105460000047
分别为电力系统的第b节点处火电机组、可再生能源发电机组、功率外送线路和储能单元的集合;
(2-2-8)电力系统线路功率约束;
Figure BDA0001994105460000048
其中,
Figure BDA0001994105460000049
为电力系统第l条输电线路的传输容量;
(2-2-9)电力系统电能外送线路功率约束;
Figure BDA00019941054600000410
其中,
Figure BDA00019941054600000411
为电力系统第i条电能外送线路的最大传输功率,
Figure BDA00019941054600000412
为电力系统第i条电能外送线路的最小传输功率;
(2-2-10)电力系统储能单元的运行约束为;
其中,储能单元的充放电功率需满足其最大功率约束为:
Figure BDA00019941054600000413
Figure BDA00019941054600000414
储能单元运行状态的逻辑约束,表达式为:
Figure BDA00019941054600000415
Figure BDA00019941054600000416
储能单元的电量与充放电功率之间的关系约束为:
Figure BDA0001994105460000051
储能单元的电量需满足其最大、最小充电状态的约束为:
Figure BDA0001994105460000052
储能单元每日开始与结束时的电量需保持相等,以保证储能单元的长期稳定运行,表达式为:
Figure BDA0001994105460000053
其中,
Figure BDA0001994105460000054
为第j个储能单元在第s个电力系统运行场景下第t控制时段的充电/放电功率,
Figure BDA0001994105460000055
为第j个储能单元在第s个电力系统运行场景下第t控制时段充电/放电状态值,若储能单元处于充电状态,则
Figure BDA0001994105460000056
为1,
Figure BDA0001994105460000057
为0,若储能单元处于放电状态,则
Figure BDA0001994105460000058
为0,
Figure BDA0001994105460000059
为1,若储能单元处于既不充电又不放电状态,则
Figure BDA00019941054600000510
为0,
Figure BDA00019941054600000511
为0;
Figure BDA00019941054600000512
为第j个储能单元的最大输出功率,Ej,s,t为第j个储能单元在第s个电力系统运行场景下第t控制时段所存储的能量,ηj为第j个储能单元的充电效率或放电效率,Ej,max为第j个储能单元的电能容量,SOCj,min/max为第j个储能单元的最小/最大充电状态;
(2-3)利用混合整数规划求解方法,求解由上述目标函数和约束条件组成的电力系统日前市场模型,得到电力系统火电机组的启停计划和电能外送计划,即电力系统火电机组的启动状态值
Figure BDA00019941054600000513
火电机组的关机状态值
Figure BDA00019941054600000514
火电机组的运行状态值
Figure BDA00019941054600000515
和电能外送功率
Figure BDA00019941054600000516
(3)根据步骤(2-3)的电力系统火电机组的启动状态值
Figure BDA00019941054600000517
火电机组的关机状态值
Figure BDA00019941054600000518
火电机组的运行状态值
Figure BDA00019941054600000519
和电能外送功率
Figure BDA00019941054600000520
建立电力系统实时市场模型,并求解,得到发电机组与储能单元各小时的运行功率,包括以下步骤:
(3-1)建立电力系统实时市场模型的目标函数,表达式如下:
Figure BDA00019941054600000521
其中,上标-RT表示实时市场,
Figure BDA00019941054600000522
表示电力系统实时市场模型中第i条电能外送线路第t控制时段的输电功率,
Figure BDA00019941054600000523
表示日前市场模型计算得到的第i条电能外送线路第t控制时段传输功率,
Figure BDA0001994105460000061
表示电力系统中第i条电能外送线路在第t控制时段的可用传输容量;
(3-2)构建电力系统实时市场模型的约束条件,具体如下:
(3-2-1)电力系统火电机组的功率约束,表达式为:
Figure BDA0001994105460000062
其中,
Figure BDA0001994105460000063
为电力系统日前市场模型求解得到的第g个火电机组在第t控制时段运行的状态值,若火电机组在第t时段运行,
Figure BDA0001994105460000064
为1,若火电机组在第t时段关闭,
Figure BDA0001994105460000065
为0,
Figure BDA0001994105460000066
为电力系统实时市场模型中第g个火电机组第t控制时段的功率;
(3-2-2)电力系统可再生能源发电机组的功率约束,表达式为:
Figure BDA0001994105460000067
其中,
Figure BDA0001994105460000068
为电力系统实时市场模型中第r个可再生能源发电机组第t控制时段的最大输出功率,
Figure BDA0001994105460000069
为电力系统实时市场模型中第r个可再生能源发电机组在第t控制时段的输出功率;
(3-2-3)电力系统可再生能源的消纳率约束:
Figure BDA00019941054600000610
(3-2-4)电力系统功率平衡约束:
Figure BDA00019941054600000611
其中,
Figure BDA00019941054600000612
为电力系统实时市场模型中节点b处在第t控制时段的节点负荷,
Figure BDA00019941054600000613
为电力系统实时场景下第j个储能单元在第t控制时段的充电/放电功率;
(3-2-5)电力系统直流潮流约束,表达式为:
Figure BDA00019941054600000614
其中,
Figure BDA00019941054600000615
表示电力系统实时市场模型中第l条输电线路在第t控制时段的传输功率;
(3-2-6)电力系统线路功率约束,表达式为:
Figure BDA00019941054600000616
(3-2-7)储能单元的运行约束:
其中,储能单元的充放电功率需满足的最大功率约束为:
Figure BDA0001994105460000071
Figure BDA0001994105460000072
储能单元运行状态的逻辑约束,表达式为:
Figure BDA0001994105460000073
Figure BDA0001994105460000074
储能单元的电量与充电、放电功率之间的关系约束为:
Figure BDA0001994105460000075
储能单元的电量需满足最大充电状态和最小充电状态的约束为:
Figure BDA0001994105460000076
储能单元每日开始与结束时的电量需保持相等,以保证储能单元的长期稳定运行,表达式为:
Figure BDA0001994105460000077
其中,
Figure BDA0001994105460000078
为电力系统实时市场模型中第j个储能单元第t控制时段充电/放电的状态值,若储能单元处于充电状态,则
Figure BDA0001994105460000079
为1,
Figure BDA00019941054600000710
为0,若储能单元处于放电状态,则
Figure BDA00019941054600000711
为0,
Figure BDA00019941054600000712
为1,储能单元处于既不充电也不放电的状态,则
Figure BDA00019941054600000713
为0,
Figure BDA00019941054600000714
为0;
Figure BDA00019941054600000715
为电力系统实时市场模型中第j个储能单元在第t控制时段所存储的能量。
(3-3)采用混合整数规划的求解方法,求解由步骤(3-1)的目标函数和步骤(3-2)的约束条件组成的电力系统实时市场模型,确定火电机组在每个控制时段的输出功率
Figure BDA00019941054600000716
储能单元每个控制时段的充电功率
Figure BDA00019941054600000717
和储能单元每个控制时段的放电功率
Figure BDA00019941054600000718
实现基于含储能单元的电力系统可用传输容量的控制。
本发明提出的基于含储能单元的电力系统可用传输容量的控制方法,其特点和优点是:
本发明方法首先获取电力系统中储能单元、发电机组、输电线路的数据,进而建立两阶段的电力系统数学模型,通过迭代计算方法确定发电机组与储能单元的协同运行方案,即第一阶段在日前市场中求解得到系统中各发电机组的启停计划,第二阶段在实时市场中迭代求解每小时发电机组出力与储能单元充放电功率,最终得到系统中各发电机组、储能单元的调度运行方案。本发明综合考虑电力系统运行成本控制与可用传输容量评估,充分利用储能单元“削峰填谷”的调节作用,在一定程度上减轻电网规划压力,延缓新建输电线路,实现电网资产的高效利用。本发明充分利用已有储能单元的调节作用,确定发电机组与储能单元的协同运行方案,使电网规划更加科学合理,降低电网规划与系统运行的成本。
附图说明
图1是本发明的一种基于含储能单元的电力系统可用传输容量的控制方法的流程框图。
具体实施方式
本发明提出的基于含储能单元的电力系统可用传输容量的控制方法,其流程框图如图1所示,特征在于该方法包括以下步骤:
(1)从电力系统控制中心获取电力系统中储能单元、发电机组及输电线路的参数,储能单元参数分别包括:最大充电功率、最大放电功率、容量、充电效率、放电效率、最大充电状态和最小充电状态,发电机组参数分别包括:所有发电机组的装机容量、火电机组的煤耗成本、火电机组的最小运行时间和火电机组的最小冷却时间,输电线路参数包括:电力系统的输电线路传输容量,电能外送线路的最大传输功率,电能外送线路的最小传输容量和节点负荷;
本发明中建立的电力系统可用传输容量评估模型分为两阶段,第一阶段是日前市场模型,求解系统中各发电机组的启停计划;第二阶段是实时市场模型,求解每小时发电机组出力与储能单元充放电功率。
(2)建立一个电力系统日前市场模型,并求解,得到各发电机组的启停计划,具体过程如下:
(2-1)建立电力系统日前市场模型的目标函数,表达式如下:
Figure BDA0001994105460000081
该目标函数表示使电力系统发电总成本最小化,其中发电总成本为火电机组的运行及启停总成本减去电能外送的总收益;
其中,
Figure BDA0001994105460000091
表示火电机组运行与启停的总成本;
Figure BDA0001994105460000092
表示电能外送的总收益,Φ表示集合,上标S、T、G、I分别表示电力系统运行场景、控制时段、火电机组和电力系统外送输电线路,ΦT,ΦG和ΦI分别从电力系统控制中心获取,γs为电力系统运行场景s的权重,本发明采用已有的聚类算法,对电力系统的日负荷及风电、光伏、水电等可再生能源出力曲线进行聚类,得到ΦS个聚类中心和每个电力系统运行场景的权重γs,ΦS和γs从电力系统控制中心获取,
Figure BDA0001994105460000093
为第g个火电机组在第s个电力系统运行场景下第t控制时段的输出功率,Yg,t为第g个火电机组在第t控制时段的启动状态值,若火电机组在第t-1控制时段为关闭,在第t控制时段为运行,则Yg,t为1,若火电机组在第t-1控制时段为运行,在第t控制时段为关闭,则Yg,t为0,若火电机组在第t-1控制时段为运行,在第t控制时段仍为运行,则Yg,t为0,若火电机组在第t-1控制时段为关闭,在第t控制时段仍为关闭,则Yg,t为0;Zg,t为第g个火电机组在第t控制时段的关机状态值,若火电机组在第t-1控制时段为运行,在第t控制时段为关闭,则Zg,t为1,若火电机组在第t-1控制时段为关闭,在第t控制时段为运行,则Zg,t为0,若火电机组在第t-1控制时段为运行,在第t控制时段仍为运行,则Zg,t为0,若火电机组在第t-1控制时段为关闭,在第t控制时段仍为关闭,则Zg,t为0;
Figure BDA00019941054600000911
为第g个火电机组在第s个电力系统运行场景下第t控制时段的运行成本,
Figure BDA0001994105460000095
为第g个火电机组在第t控制时段的启动成本,
Figure BDA0001994105460000096
为第g个火电机组在第t控制时段的关机成本,
Figure BDA0001994105460000097
为电力系统第i条电能外送线路在第t控制时段的传输功率,
Figure BDA00019941054600000912
为第i条电能外送线路在第t控制时段电能外送的单位收益;
(2-2)建立电力系统日前市场模型的约束条件,具体如下:
(2-2-1)电力系统火电机组的功率约束如下:
Figure BDA0001994105460000098
其中,
Figure BDA0001994105460000099
Figure BDA00019941054600000910
分别为设定的第g个火电机组的最小功率限制和最大功率限制,Ug,t为第g个火电机组在第t控制时段的运行状态值,若火电机组在第t控制时段运行,则Ug,t为1,若火电机组在第t控制时段关机,则Ug,t为0;
(2-2-2)电力系统火电机组工作状态的逻辑约束,即在第t控制时段火电机组g只能处于运行或关机状态中的一种,表达式如下:
Figure BDA0001994105460000101
Figure BDA0001994105460000102
(2-2-3)电力系统火电机组的启动和关机时间约束为:
Figure BDA0001994105460000103
Figure BDA0001994105460000104
其中,
Figure BDA0001994105460000105
表示第g个火电机组连续开机的最少时间,
Figure BDA0001994105460000106
表示第g个火电机组保持关机状态最少时间;
(2-2-4)电力系统中可再生能源发电机组的功率约束,该约束表示第r个可再生能源发电机组输出功率受到电力系统控制中心给出的预测值的约束,表达式为:
Figure BDA0001994105460000107
其中,上标R表示可再生能源发电机组,
Figure BDA0001994105460000108
为第r个可再生能源发电机组第s个电力系统运行场景下在第t控制时段的最大输出功率的预测值,
Figure BDA0001994105460000109
为第r个可再生能源发电机组第s个电力系统运行场景下在第t控制时段的输出功率;
(2-2-5)电力系统可再生能源的消纳率约束,该约束表示可再生能源的消纳率需高于所设定的约束下限,表达式为:
Figure BDA00019941054600001010
其中,RA为设定的可再生能源消纳率的约束下限,ΦR为可再生能源的发电机组集合;
(2-2-6)电力系统有功功率平衡约束为:
Figure BDA00019941054600001011
其中,上标D表示电力系统节点负荷,上标CHA表示电力系统储能单元充电,上标DIS表示储能单元放电,上标ES表示储能单元,ΦD和ΦES分别为节点负荷集合和储能单元集合,
Figure BDA00019941054600001012
为电力系统中节点b处第s个场景下第t控制时段的负荷,
Figure BDA00019941054600001013
表示第j个储能单元在第s个电力系统运行场景下第t个控制时段的充电功率,
Figure BDA0001994105460000111
表示第j个储能单元在第s个电力系统运行场景下第t个控制时段的放电功率;
(2-2-7)电力系统直流潮流约束,该约束以电力系统直流潮流的形式表示输电线路满足的电路潮流约束,表达式为:
Figure BDA0001994105460000112
其中,L表示电力系统中的输电线路,
Figure BDA0001994105460000113
为电力系统中第l条输电线路在第s个电力系统运行场景下第t控制时段的传输功率,Fl-b表示第l条输电线路对第b个节点的发电转移分布因子(generation shift distribution factor,GSDF),Fl-b从电力系统控制中心获取,
Figure BDA0001994105460000114
Figure BDA0001994105460000115
Figure BDA0001994105460000116
分别为电力系统的第b节点处火电机组、可再生能源发电机组、功率外送线路和储能单元的集合;
(2-2-8)电力系统线路功率约束;
Figure BDA0001994105460000117
其中,
Figure BDA0001994105460000118
为电力系统第l条输电线路的传输容量;该约束表示输电线路输电功率受其传输容量的限制;
(2-2-9)电力系统电能外送线路功率约束;
Figure BDA0001994105460000119
其中,
Figure BDA00019941054600001110
为电力系统第i条电能外送线路的最大传输功率,
Figure BDA00019941054600001111
为电力系统第i条电能外送线路的最小传输功率;该约束表示电能外送线路的传输功率受所设定的传输功率上下限的限制;
(2-2-10)电力系统储能单元的运行约束为;
其中,储能单元的充放电功率需满足其最大功率约束为:
Figure BDA00019941054600001112
Figure BDA00019941054600001113
储能单元运行状态的逻辑约束,即在同一电力系统运行场景同一控制时段下储能单元只能处于充电
Figure BDA00019941054600001114
放电
Figure BDA00019941054600001115
或既不充电也不放电
Figure BDA00019941054600001116
状态中的一种,表达式为:
Figure BDA0001994105460000121
Figure BDA0001994105460000122
储能单元的电量与充放电功率之间的关系约束为:
Figure BDA0001994105460000123
储能单元的电量需满足其最大、最小充电状态的约束为:
Figure BDA0001994105460000124
储能单元每日开始与结束时的电量需保持相等,以保证储能单元的长期稳定运行,表达式为:
Figure BDA0001994105460000125
其中,
Figure BDA0001994105460000126
为第j个储能单元在第s个电力系统运行场景下第t控制时段的充电/放电功率,
Figure BDA0001994105460000127
为第j个储能单元在第s个电力系统运行场景下第t控制时段充电/放电状态值,若储能单元处于充电状态,则
Figure BDA0001994105460000128
为1,
Figure BDA0001994105460000129
为0,若储能单元处于放电状态,则
Figure BDA00019941054600001210
为0,
Figure BDA00019941054600001211
为1,若储能单元处于既不充电又不放电状态,则
Figure BDA00019941054600001212
为0,
Figure BDA00019941054600001213
为0;
Figure BDA00019941054600001214
为第j个储能单元的最大输出功率,Ej,s,t为第j个储能单元在第s个电力系统运行场景下第t控制时段所存储的能量,ηj为第j个储能单元的充电效率或放电效率,Ej,max为第j个储能单元的电能容量,SOCj,min/max为第j个储能单元的最小/最大充电状态,例如SOCj,min=5%,SOCj,max=95%,则表示第j个储能单元的最低可放电至剩余电量为容量的5%,最高可充电至电量为容量的95%。
(2-3)利用混合整数规划求解方法,求解由上述目标函数和约束条件组成的电力系统日前市场模型,得到电力系统火电机组的启停计划和电能外送计划,即电力系统火电机组的启动状态值
Figure BDA00019941054600001215
火电机组的关机状态值
Figure BDA00019941054600001216
火电机组的运行状态值
Figure BDA00019941054600001217
和电能外送功率
Figure BDA00019941054600001218
(3)根据步骤(2-3)的电力系统火电机组的启动状态值
Figure BDA00019941054600001219
火电机组的关机状态值
Figure BDA00019941054600001220
火电机组的运行状态值
Figure BDA00019941054600001221
和电能外送功率
Figure BDA00019941054600001222
建立电力系统实时市场模型,并求解,得到发电机组与储能单元各小时的运行功率,包括以下步骤:
(3-1)建立电力系统实时市场模型的目标函数,表达式如下:
Figure BDA0001994105460000131
该目标函数表示电力系统功率外送线路的可用传输容量总和最大化;
其中,上标-RT表示实时市场,
Figure BDA0001994105460000132
表示电力系统实时市场模型中第i条电能外送线路第t控制时段的输电功率,
Figure BDA0001994105460000133
表示日前市场模型计算得到的第i条电能外送线路第t控制时段传输功率,
Figure BDA0001994105460000134
表示电力系统中第i条电能外送线路在第t控制时段的可用传输容量;
(3-2)构建电力系统实时市场模型的约束条件,具体如下:
(3-2-1)电力系统火电机组的功率约束,该约束表示火电机组的实时运行功率需满足功率上下限约束,表达式为:
Figure BDA0001994105460000135
其中,
Figure BDA0001994105460000136
为电力系统日前市场模型求解得到的第g个火电机组在第t控制时段运行的状态值,若火电机组在第t时段运行,
Figure BDA0001994105460000137
为1,若火电机组在第t时段关闭,
Figure BDA0001994105460000138
为0,
Figure BDA0001994105460000139
为电力系统实时市场模型中第g个火电机组第t控制时段的功率;
(3-2-2)电力系统可再生能源发电机组的功率约束,该约束表示可再生能源发电机组输出功率受到其最大输出功率的约束,表达式为:
Figure BDA00019941054600001310
其中,
Figure BDA00019941054600001311
为电力系统实时市场模型中第r个可再生能源发电机组第t控制时段的最大输出功率,
Figure BDA00019941054600001312
为电力系统实时市场模型中第r个可再生能源发电机组在第t控制时段的输出功率;
(3-2-3)电力系统可再生能源的消纳率约束:
Figure BDA00019941054600001313
该约束表示可再生能源的消纳率需高于所设定的约束下限;
(3-2-4)电力系统功率平衡约束:
Figure BDA00019941054600001314
其中,
Figure BDA0001994105460000141
为电力系统实时市场模型中节点b处在第t控制时段的节点负荷,
Figure BDA0001994105460000142
为电力系统实时场景下第j个储能单元在第t控制时段的充电/放电功率;
(3-2-5)电力系统直流潮流约束,该约束以直流潮流的形式表示电力系统输电线路满足的电路潮流约束,表达式为:
Figure BDA0001994105460000143
其中,
Figure BDA0001994105460000144
表示电力系统实时市场模型中第l条输电线路在第t控制时段的传输功率;
(3-2-6)电力系统线路功率约束,该约束表示输电线路输电容量受传输功率上限的限制,表达式为:
Figure BDA0001994105460000145
(3-2-7)储能单元的运行约束:
其中,储能单元的充放电功率需满足的最大功率约束为:
Figure BDA0001994105460000146
Figure BDA0001994105460000147
储能单元运行状态的逻辑约束,即在同一电力系统运行场景同一控制时段下储能单元只能处于充电
Figure BDA0001994105460000148
放电
Figure BDA0001994105460000149
或既不充电也不放电
Figure BDA00019941054600001410
状态中的一种,表达式为:
Figure BDA00019941054600001411
Figure BDA00019941054600001412
储能单元的电量与充电、放电功率之间的关系约束为:
Figure BDA00019941054600001413
储能单元的电量需满足最大充电状态和最小充电状态的约束为:
Figure BDA00019941054600001414
储能单元每日开始与结束时的电量需保持相等,以保证储能单元的长期稳定运行,表达式为:
Figure BDA00019941054600001415
其中,
Figure BDA0001994105460000151
为电力系统实时市场模型中第j个储能单元第t控制时段充电/放电的状态值,若储能单元处于充电状态,则
Figure BDA0001994105460000152
为1,
Figure BDA0001994105460000153
为0,若储能单元处于放电状态,则
Figure BDA0001994105460000154
为0,
Figure BDA0001994105460000155
为1,储能单元处于既不充电也不放电的状态,则
Figure BDA0001994105460000156
为0,
Figure BDA0001994105460000157
为0;
Figure BDA0001994105460000158
为电力系统实时市场模型中第j个储能单元在第t控制时段所存储的能量。
(3-3)采用混合整数规划的求解方法,求解由步骤(3-1)的目标函数和步骤(3-2)的约束条件组成的电力系统实时市场模型,确定火电机组在每个控制时段的输出功率
Figure BDA0001994105460000159
储能单元在每个控制时段的充电功率
Figure BDA00019941054600001510
和储能单元在每个控制时段的放电功率
Figure BDA00019941054600001511
实现基于含储能单元的电力系统可用传输容量的控制。
根据日前市场与实时市场的两阶段模型的优化结果,确定外送线路的输电计划,确定发电机组与储能单元的运行方案,计算得到电力系统可用传输容量结果。至此,本发明所提方法实施完毕。以上实施步骤仅用以说明而非限制本发明的技术方案。不脱离本发明精神和范围的任何修改或局部替换,均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种基于含储能单元的电力系统可用传输容量的控制方法,其特征在于该方法包括以下步骤:
(1)从电力系统控制中心获取电力系统中储能单元、发电机组及输电线路的参数,储能单元参数分别包括:最大充电功率、最大放电功率、容量、充电效率、放电效率、最大充电状态和最小充电状态,发电机组参数分别包括:所有发电机组的装机容量、火电机组的煤耗成本、火电机组的最小运行时间和火电机组的最小冷却时间,输电线路参数包括:电力系统的输电线路传输容量,电能外送线路的最大传输功率,电能外送线路的最小传输容量和节点负荷;
(2)建立一个电力系统日前市场模型,并求解,得到各发电机组的启停计划,具体过程如下:
(2-1)建立电力系统日前市场模型的目标函数,表达式如下:
Figure FDA0002477137530000011
其中,
Figure FDA0002477137530000012
表示火电机组运行与启停的总成本;
Figure FDA0002477137530000013
表示电能外送的总收益,Φ表示集合,上标S、T、G、I分别表示电力系统运行场景、控制时段、火电机组和电力系统外送输电线路,ΦTG和ΦI分别从电力系统控制中心获取,γs为电力系统运行场景s的权重,Φs和γs从电力系统控制中心获取,
Figure FDA0002477137530000014
为第g个火电机组在第s个电力系统运行场景下第t控制时段的输出功率,Yg,t为第g个火电机组在第t控制时段的启动状态值,若火电机组在第t-1控制时段为关闭,在第t控制时段为运行,则Yg,t为1,若火电机组在第t-1控制时段为运行,在第t控制时段为关闭,则Yg,t为0,若火电机组在第t-1控制时段为运行,在第t控制时段仍为运行,则Yg,t为0,若火电机组在第t-1控制时段为关闭,在第t控制时段仍为关闭,则Yg,t为0;Zg,t为第g个火电机组在第t控制时段的关机状态值,若火电机组在第t-1控制时段为运行,在第t控制时段为关闭,则Zg,t为1,若火电机组在第t-1控制时段为关闭,在第t控制时段为运行,则Zg,t为0,若火电机组在第t-1控制时段为运行,在第t控制时段仍为运行,则Zg,t为0,若火电机组在第t-1控制时段为关闭,在第t控制时段仍为关闭,则Zg,为0;
Figure FDA0002477137530000021
为第g个火电机组在第s个电力系统运行场景下第t控制时段的运行成本,
Figure FDA0002477137530000022
为第g个火电机组在第t控制时段的启动成本,
Figure FDA0002477137530000023
为第g个火电机组在第t控制时段的关机成本,
Figure FDA0002477137530000024
为电力系统第i条电能外送线路在第t控制时段的传输功率,
Figure FDA0002477137530000025
为第i条电能外送线路在第t控制时段电能外送的单位收益;
(2-2)建立电力系统日前市场模型的约束条件,具体如下:
(2-2-1)电力系统火电机组的功率约束如下:
Figure FDA0002477137530000026
其中,
Figure FDA0002477137530000027
Figure FDA0002477137530000028
分别为设定的第g个火电机组的最小功率限制和最大功率限制,Ug,t为第g个火电机组在第t控制时段的运行状态值,若火电机组在第t控制时段运行,则Ug,t为1,若火电机组在第t控制时段关机,则Ug,t为0;
(2-2-2)电力系统火电机组工作状态的逻辑约束,即在第t控制时段火电机组g只能处于运行或关机状态中的一种,表达式如下:
Figure FDA0002477137530000029
Figure FDA00024771375300000210
(2-2-3)电力系统火电机组的启动和关机时间约束为:
Figure FDA00024771375300000211
Figure FDA00024771375300000212
其中,
Figure FDA00024771375300000213
表示第g个火电机组连续开机的最少时间,
Figure FDA00024771375300000214
表示第g个火电机组保持关机状态最少时间;
(2-2-4)电力系统中可再生能源发电机组的功率约束,表达式为:
Figure FDA00024771375300000215
其中,上标R表示可再生能源发电机组,
Figure FDA00024771375300000216
为第r个可再生能源发电机组第s个电力系统运行场景下在第t控制时段的最大输出功率的预测值,
Figure FDA00024771375300000217
为第r个可再生能源发电机组第s个电力系统运行场景下在第t控制时段的输出功率;
(2-2-5)电力系统可再生能源的消纳率约束,表达式为:
Figure FDA0002477137530000031
其中,RA为设定的可再生能源消纳率的约束下限,ΦR为可再生能源的发电机组集合;
(2-2-6)电力系统有功功率平衡约束为:
Figure FDA0002477137530000032
其中,上标D表示电力系统节点负荷,上标CHA表示电力系统储能单元充电,上标DIS表示储能单元放电,上标ES表示储能单元,ΦD和ΦES分别为节点负荷集合和储能单元集合,
Figure FDA0002477137530000033
为电力系统中节点b处第s个场景下第t控制时段的负荷,
Figure FDA0002477137530000034
表示第j个储能单元在第s个电力系统运行场景下第t个控制时段的充电功率,
Figure FDA0002477137530000035
表示第j个储能单元在第s个电力系统运行场景下第t个控制时段的放电功率;
(2-2-7)电力系统直流潮流约束,表达式为:
Figure FDA0002477137530000036
其中,L表示电力系统中的输电线路,
Figure FDA0002477137530000037
为电力系统中第l条输电线路在第s个电力系统运行场景下第t控制时段的传输功率,Fl-b表示第l条输电线路对第b个节点的发电转移分布因子,Fl-b从电力系统控制中心获取,
Figure FDA0002477137530000038
Figure FDA0002477137530000039
分别为电力系统的第b节点处火电机组、可再生能源发电机组、功率外送线路和储能单元的集合;
(2-2-8)电力系统线路功率约束;
Figure FDA00024771375300000310
其中,
Figure FDA00024771375300000311
为电力系统第l条输电线路的最大传输容量;
(2-2-9)电力系统电能外送线路功率约束;
Figure FDA00024771375300000312
其中,
Figure FDA00024771375300000313
为电力系统第i条电能外送线路的最大传输功率,
Figure FDA00024771375300000314
为电力系统第i条电能外送线路的最小传输功率;
(2-2-10)电力系统储能单元的运行约束为;
其中,储能单元的充放电功率需满足其最大功率约束为:
Figure FDA0002477137530000041
Figure FDA0002477137530000042
储能单元运行状态的逻辑约束,表达式为:
Figure FDA0002477137530000043
Figure FDA0002477137530000044
储能单元的电量与充放电功率之间的关系约束为:
Figure FDA0002477137530000045
储能单元的电量需满足其最大、最小充电状态的约束为:
Figure FDA0002477137530000046
储能单元每日开始与结束时的电量需保持相等,以保证储能单元的长期稳定运行,表达式为:
Figure FDA0002477137530000047
其中,
Figure FDA0002477137530000048
为第j个储能单元在第s个电力系统运行场景下第t控制时段的充电/放电功率,
Figure FDA0002477137530000049
为第j个储能单元在第s个电力系统运行场景下第t控制时段充电/放电状态值,若储能单元处于充电状态,则
Figure FDA00024771375300000410
为1,
Figure FDA00024771375300000411
为0,若储能单元处于放电状态,则
Figure FDA00024771375300000412
为0,
Figure FDA00024771375300000413
为1,若储能单元处于既不充电又不放电状态,则
Figure FDA00024771375300000414
为0,
Figure FDA00024771375300000415
为0;
Figure FDA00024771375300000416
为第j个储能单元的最大输出功率,Ej,s,t为第j个储能单元在第s个电力系统运行场景下第t控制时段所存储的能量,ηj为第j个储能单元的充电效率或放电效率,Ej,max为第j个储能单元的最大电能容量,SOCj,min/max为第j个储能单元的最小/最大充电状态;
(2-3)利用混合整数规划求解方法,求解由上述目标函数和约束条件组成的电力系统日前市场模型,得到电力系统火电机组的启停计划和电能外送计划,即电力系统火电机组的启动状态值
Figure FDA00024771375300000417
火电机组的关机状态值
Figure FDA00024771375300000418
火电机组的运行状态值
Figure FDA00024771375300000419
和电能外送功率
Figure FDA00024771375300000420
(3)根据步骤(2-3)的电力系统火电机组的启动状态值
Figure FDA00024771375300000421
火电机组的关机状态值
Figure FDA0002477137530000051
火电机组的运行状态值
Figure FDA0002477137530000052
和电能外送功率
Figure FDA0002477137530000053
建立电力系统实时市场模型,并求解,得到发电机组与储能单元各小时的运行功率,包括以下步骤:
(3-1)建立电力系统实时市场模型的目标函数,表达式如下:
Figure FDA0002477137530000054
其中,上标-RT表示实时市场,
Figure FDA0002477137530000055
表示电力系统实时市场模型中第i条电能外送线路第t控制时段的输电功率,
Figure FDA0002477137530000056
表示日前市场模型计算得到的第i条电能外送线路第t控制时段传输功率,
Figure FDA0002477137530000057
表示电力系统中第i条电能外送线路在第t控制时段的可用传输容量;
(3-2)构建电力系统实时市场模型的约束条件,具体如下:
(3-2-1)电力系统火电机组的功率约束,表达式为:
Figure FDA0002477137530000058
其中,
Figure FDA0002477137530000059
为电力系统日前市场模型求解得到的第g个火电机组在第t控制时段运行的状态值,若火电机组在第t时段运行,
Figure FDA00024771375300000510
为1,若火电机组在第t时段关闭,
Figure FDA00024771375300000511
为0,
Figure FDA00024771375300000512
为电力系统实时市场模型中第g个火电机组第t控制时段的功率;
(3-2-2)电力系统可再生能源发电机组的功率约束,表达式为:
Figure FDA00024771375300000513
其中,
Figure FDA00024771375300000514
为电力系统实时市场模型中第r个可再生能源发电机组第t控制时段的最大输出功率,
Figure FDA00024771375300000515
为电力系统实时市场模型中第r个可再生能源发电机组在第t控制时段的输出功率;
(3-2-3)电力系统可再生能源的消纳率约束:
Figure FDA00024771375300000516
(3-2-4)电力系统功率平衡约束:
Figure FDA00024771375300000517
其中,
Figure FDA00024771375300000518
为电力系统实时场景下第j个储能单元在第t控制时段的充电/放电功率;
(3-2-5)电力系统直流潮流约束,表达式为:
Figure FDA0002477137530000061
其中,
Figure FDA0002477137530000062
表示电力系统实时市场模型中第l条输电线路在第t控制时段的传输功率;
(3-2-6)电力系统线路功率约束,表达式为:
Figure FDA0002477137530000063
(3-2-7)储能单元的运行约束:
其中,储能单元的充放电功率需满足的最大功率约束为:
Figure FDA0002477137530000064
Figure FDA0002477137530000065
储能单元运行状态的逻辑约束,表达式为:
Figure FDA0002477137530000066
Figure FDA0002477137530000067
储能单元的电量与充电、放电功率之间的关系约束为:
Figure FDA0002477137530000068
储能单元的电量需满足最大充电状态和最小充电状态的约束为:
Figure FDA0002477137530000069
储能单元每日开始与结束时的电量需保持相等,以保证储能单元的长期稳定运行,表达式为:
Figure FDA00024771375300000610
其中,
Figure FDA00024771375300000611
为电力系统实时市场模型中第j个储能单元第t控制时段充电/放电的状态值,若储能单元处于充电状态,则
Figure FDA00024771375300000612
为1,
Figure FDA00024771375300000613
为0,若储能单元处于放电状态,则
Figure FDA00024771375300000614
为0,
Figure FDA00024771375300000615
为1,储能单元处于既不充电也不放电的状态,则
Figure FDA00024771375300000616
为0,
Figure FDA00024771375300000617
为0;
Figure FDA00024771375300000618
为电力系统实时市场模型中第j个储能单元在第t控制时段所存储的能量;
(3-3)采用混合整数规划的求解方法,求解由步骤(3-1)的目标函数和步骤(3-2)的约束条件组成的电力系统实时市场模型,确定火电机组在每个控制时段的输出功率
Figure FDA0002477137530000071
储能单元在每个控制时段的充电功率
Figure FDA0002477137530000072
和储能单元在每个控制时段的放电功率
Figure FDA0002477137530000073
实现基于含储能单元的电力系统可用传输容量的控制。
CN201910190005.8A 2019-03-13 2019-03-13 一种基于含储能单元的电力系统可用传输容量的控制方法 Active CN109861289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910190005.8A CN109861289B (zh) 2019-03-13 2019-03-13 一种基于含储能单元的电力系统可用传输容量的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910190005.8A CN109861289B (zh) 2019-03-13 2019-03-13 一种基于含储能单元的电力系统可用传输容量的控制方法

Publications (2)

Publication Number Publication Date
CN109861289A CN109861289A (zh) 2019-06-07
CN109861289B true CN109861289B (zh) 2020-07-10

Family

ID=66900644

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910190005.8A Active CN109861289B (zh) 2019-03-13 2019-03-13 一种基于含储能单元的电力系统可用传输容量的控制方法

Country Status (1)

Country Link
CN (1) CN109861289B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110346649B (zh) * 2019-07-09 2021-07-09 国网河北省电力有限公司衡水供电分公司 一种储能缓存器容量的确定方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170338659A1 (en) * 2014-10-28 2017-11-23 Sinewatts, Inc. Systems and methods for dispatching maximum available capacity for photovoltaic power plants
CN104779631B (zh) * 2014-12-31 2017-07-18 国家电网公司 基于风电预测功率的储能跟踪风电出力计划方法及其系统
US9871396B2 (en) * 2016-06-21 2018-01-16 General Electric Company System and method for controlling charge of an energy storage device from a renewable energy source
CN106683003B (zh) * 2016-12-23 2022-09-30 中国电力科学研究院 一种电力可用传输容量确定方法及装置
CN106779498A (zh) * 2017-01-24 2017-05-31 国网福建省电力有限公司 评估含风光储的大用户参与日前市场电力交易风险的方法
CN109193623B (zh) * 2018-09-10 2021-11-26 广东工业大学 一种结合电动汽车的日前市场和实时市场联合优化方法

Also Published As

Publication number Publication date
CN109861289A (zh) 2019-06-07

Similar Documents

Publication Publication Date Title
Sun et al. An optimal combined operation scheme for pumped storage and hybrid wind-photovoltaic complementary power generation system
Roslan et al. Scheduling controller for microgrids energy management system using optimization algorithm in achieving cost saving and emission reduction
CN109193812B (zh) 一种园区光储荷微电网经济调度实现方法
WO2022100091A1 (zh) 一种广义源储系统调度的集中控制方法
Kazempour et al. Electric energy storage systems in a market-based economy: Comparison of emerging and traditional technologies
Chakraborty et al. PV-microgrid operational cost minimization by neural forecasting and heuristic optimization
CN109217290B (zh) 计及电动汽车充放电的微网能量优化管理方法
CN111525627B (zh) 含抽蓄与新能源发电的柔性直流输电系统日前调度方法
CN107425534B (zh) 一种基于优化蓄电池充放电策略的微电网调度方法
CN101662244B (zh) 功耗控制方法和设备
CN103217900A (zh) 基于实时电价的中压微网混沌pso最优潮流实现方法
CN110633854A (zh) 一种考虑储能电池多重分段服务全生命周期优化规划方法
CN111934360A (zh) 基于模型预测控制的虚拟电厂-储能系统能量协同优化调控方法
CN109861289B (zh) 一种基于含储能单元的电力系统可用传输容量的控制方法
CN110571868B (zh) 微电网的优化配置方法
CN112311017A (zh) 一种虚拟电厂与主网最优协同调度方法
CN112307603A (zh) 考虑大规模风电接入的混合储能容量优化配置方法及系统
El-Hendawi et al. Enhanced MG with optimum operational cost of pumping water distribution systems
CN110826210A (zh) 基于功率互联的多区域楼宇虚拟电厂建模及优化协调方法
CN115995832A (zh) 一种计及柔性负荷的户用氢能系统优化运行方法及装置
Pasban-Gajan et al. Optimal scheduling of renewable-based energy hubs considering time-of-use pricing scheme
CN115882460A (zh) 一种考虑需求侧管理的两阶段新能源微电网优化调度方法
CN114465226A (zh) 一种电力系统多级备用获取联合优化模型的建立方法
Hongli et al. Day-ahead optimal dispatch of regional power grid based on electric vehicle participation in peak shaving pricing strategy
CN116979570B (zh) 微电网混合型储能方法、系统及装置

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