CN115879734A - 一种混合式抽水蓄能电站与风电联合运行调度方法 - Google Patents

一种混合式抽水蓄能电站与风电联合运行调度方法 Download PDF

Info

Publication number
CN115879734A
CN115879734A CN202310008506.6A CN202310008506A CN115879734A CN 115879734 A CN115879734 A CN 115879734A CN 202310008506 A CN202310008506 A CN 202310008506A CN 115879734 A CN115879734 A CN 115879734A
Authority
CN
China
Prior art keywords
power
unit
pumped
storage
constraint
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.)
Pending
Application number
CN202310008506.6A
Other languages
English (en)
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.)
Tsinghua University
State Grid Sichuan Electric Power Co Ltd
Sichuan Energy Internet Research Institute EIRI Tsinghua University
Original Assignee
Tsinghua University
State Grid Sichuan Electric Power Co Ltd
Sichuan Energy Internet Research Institute EIRI 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 Tsinghua University, State Grid Sichuan Electric Power Co Ltd, Sichuan Energy Internet Research Institute EIRI Tsinghua University filed Critical Tsinghua University
Priority to CN202310008506.6A priority Critical patent/CN115879734A/zh
Publication of CN115879734A publication Critical patent/CN115879734A/zh
Pending legal-status Critical Current

Links

Images

Abstract

本发明的实施例提供了一种混合式抽水蓄能电站与风电联合运行调度方法,涉及新能源技术领域。混合式抽水蓄能电站与风电联合运行调度方法包括:S1:建立混合式抽水蓄能电站与风电联合运行调度模型;其中,调度模型以联合体整体收益最大为目标函数,用公式描述为:maxF=F1+F2‑F3‑F4,F为调度期的联合体整体运行效益,单位为元;F1、F2、F3、F4分别为联合体按发电计划的基础收益、实际出力高于发电计划的额外收益、开停机成本和实际出力低于发电计划的惩罚费用;S2:求解混合式抽水蓄能电站与风电联合运行调度模型。该方法能够提高联合体整体发电收益。

Description

一种混合式抽水蓄能电站与风电联合运行调度方法
技术领域
本发明涉及新能源技术领域,具体而言,涉及一种混合式抽水蓄能电站与风电联合运行调度方法。
背景技术
目前,抽水蓄能与风电等新能源的联合运行调度已有较多研究,例如,以系统负荷跟踪为目标的风电与抽水蓄能联合调度,以风电-抽水蓄能联合体输出功率平滑为目标联合调度,以减少电网弃风为目标的风电与抽水蓄能协调运行,以系统运行成本最低为目标的含抽水蓄能的多能源协调调度。然而,针对混合式抽水蓄能电站的运行调度的研究相对较少,已有的研究主要集中在考虑混合式抽水蓄能的梯级水电站群运行调度,但类似常规水电融合改造的混合式抽水蓄能与新能源的运行调度建模研究鲜见报道,与常规水电和抽水蓄能调度建模相比具有新的特点和难点,主要体现在:
(1)传统抽水蓄能电站的上下库一般为独立修建的蓄水水库,无天然径流来水,日常运行时仅简单地考虑了水力约束,而融合改造的混合式抽水蓄能电站利用已有常规电站的水库,相比传统抽水蓄能电站需要考虑更为复杂的水力联系;
(2)混合式抽水蓄能具备鲜明的“量调并重”特点,一方面既要确保电站的发电经济性,另一方面也需要协调促进新能源的消纳,提升联合系统的可调度性,如何发挥混合式抽水蓄能的优势选择合适的建模目标也是关键;
(3)混合式抽水蓄能电站包括只“发”的常规水电机组和“抽-发”兼具的抽水蓄能机组,如何建模体现差异化运行特性并实现运行工况的关联切换,并保证求解效率,也是模型建模和求解中面临的难点问题。
发明内容
本发明的目的包括提供了一种混合式抽水蓄能电站与风电联合运行调度方法,其能够提高联合体整体发电收益。
本发明的实施例可以这样实现:
本发明的实施例提供一种混合式抽水蓄能电站与风电联合运行调度方法,方法包括:
S1:建立混合式抽水蓄能电站与风电联合运行调度模型;
其中,调度模型以联合体整体收益最大为目标函数,用公式描述为:
maxF=F1+F2-F3-F4 (1)
其中,F为调度期的联合体整体运行效益,单位为元;F1、F2、F3、F4分别为联合体按发电计划的基础收益、实际出力高于发电计划的额外收益、开停机成本和实际出力低于发电计划的惩罚费用;
S2:求解混合式抽水蓄能电站与风电联合运行调度模型。
本发明实施例提供的混合式抽水蓄能电站与风电联合运行调度方法的有益效果包括:
以混合式抽水蓄能电站和风电所构成联合体整体的收益最大为目标函数,并通过引入奖惩电价机制,可以有效激励联合体提高风电的预测精度并制定合理的发电计划,从而提高了发电计划的可执行性。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例提供的混合式抽水蓄能电站与风电联合运行调度方法的流程图;
图2为风电机组预测出力及联合体发电计划示意图;
图3为混合式抽水蓄能电站上游水库水位变化过程示意图;
图4为抽水蓄能机组的抽水和发电过程的示意图;
图5为常规水电机组的出力过程的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
需要说明的是,在不冲突的情况下,本发明的实施例中的特征可以相互结合。
为探索水电站融合改造后的混合式抽水蓄能电站的运行模式,本实施例提出了一种混合式抽水蓄能电站与风电联合运行调度模型(以下简称:调度模型)。该调度模型以联合体整体收益最大为目标,以机组为最小调度单元,针对常规水电机组和抽水蓄能机组的差异化运行特性分别精细化建模,并引入状态变量实现运行状态的解耦与关联切换。在模型求解方面,通过线性化方法及建模技巧将原模型转换为MILP模型,最后在JAVA环境中采用CPLEX工具进行求解。以西南某流域电站为参考构建的应用示例验证了本文模型和求解方法的有效性,可为推进常规水电站的融合改造提供借鉴。
请参考图1,本实施例提供了一种混合式抽水蓄能电站与风电联合运行调度方法(以下简称:方法),方法包括以下步骤:
S1:建立混合式抽水蓄能电站与风电联合运行调度模型。
具体的,S1包括:
S11:建立调度模型的目标函数。
混合式抽水蓄能电站与风电联合运行,既可以利用抽水蓄能机组的快速调节能力跟随风电出力波动,降低计划执行偏差考核,又可以同时发挥常规水电机组的调峰和发电作用。因此,调度模型以联合体整体收益最大为目标函数,用公式描述为:
maxF=F1+F2-F3-F4 (1)
其中,F为调度期的联合体整体运行效益,单位为元;F1、F2、F3、F4分别为联合体按发电计划的基础收益、实际出力高于发电计划的额外收益、开停机成本和实际出力低于发电计划的惩罚费用,计算公式分别如下:
Figure BDA0004036776290000041
Figure BDA0004036776290000042
Figure BDA0004036776290000043
Figure BDA0004036776290000044
其中,ct
Figure BDA0004036776290000045
分别为时段t的上网电价、出力高于发电计划的电价、出力低于发电计划的惩罚电价,单位为元/kWh;Pt plan为联合体发电计划,单位为MW;Pt wind为时段t的预测风电出力,单位为MW;/>
Figure BDA0004036776290000046
为时段t机组i的发电功率;/>
Figure BDA0004036776290000051
为时段t机组j的抽水功率,单位为MW;N、Npump分别为总机组台数和抽水蓄能机组的台数;con、coff分别为机组的启成本和停成本,单位为元/台次;/>
Figure BDA0004036776290000052
分别为t时段启的机组数量和停的机组数量,单位为台;t∈{1,2,…,T}为时段编号,T为调度期内总的时段数。
S12:建立目标函数的约束条件。
为精细化调度模型,目标函数的约束条件包括电站约束、常规水电机组约束和抽水蓄能机组约束。
其中,电站约束包括以下约束:
(1)上下库水量平衡约束
Figure BDA00040367762900000511
/>
Figure BDA0004036776290000053
其中:
Figure BDA0004036776290000054
分别表示混合式抽水蓄能电站在时段t的上库库容和下库库容,单位为m3;/>
Figure BDA0004036776290000055
分别为时段t上库的来水流量、下库来水流量、发电流量和抽水流量,单位为m3/s;Δt为计算时段步长,单位为h。
(2)库容大小约束
Figure BDA0004036776290000056
Figure BDA00040367762900000512
Figure BDA00040367762900000513
Figure BDA00040367762900000514
其中:
Figure BDA0004036776290000057
分别为抽水蓄能电站上库的最小库容和最大库容,单位为m3
Figure BDA0004036776290000058
分别为抽水蓄能电站下库的最小库容和最大库容,单位为m3;/>
Figure BDA0004036776290000059
分别为上库的初始库容和下库的初始库容;/>
Figure BDA00040367762900000510
分别为第1时段上库的库容和下库的库容。
(3)水位-库容关系约束
Figure BDA0004036776290000061
Figure BDA00040367762900000611
其中:
Figure BDA0004036776290000062
分别为水库上库的水位-库容关系曲线函数和下库的水位-库容关系曲线函数;/>
Figure BDA0004036776290000063
分别为时段t的上游水库的水位和下游水库的水位,单位为m。
(4)电站机组总数约束
Npump+Ntrad=N (14)
其中:Npump、Ntrad、N分别为抽水蓄能机组台数、常规水电机组台数以及总机组台数,单位为台。
(5)电站发电与抽水互斥约束
Figure BDA0004036776290000064
Figure BDA00040367762900000612
Figure BDA00040367762900000613
/>
其中:
Figure BDA0004036776290000065
分别为机组在时段t的发电功率和抽水功率,单位为MW;/>
Figure BDA0004036776290000066
为时段t常规水电机组i的发电功率;/>
Figure BDA0004036776290000067
为时段t抽水蓄能电机组j的发电功率;任意时刻,电站不能出现同时抽水和发电的情况,发电与抽水状态互斥。
(6)调度期末水位控制约束
Figure BDA0004036776290000068
其中:
Figure BDA0004036776290000069
为调度期末的上游水库的水位;/>
Figure BDA00040367762900000610
为第1时段的上游水库的水位;δ为调度期内允许的水位变幅,目的是保证下一个调度期的调度需求。
(7)互补系统发电计划约束
Figure BDA0004036776290000071
其中:ε是给定的发电计划偏差系数,公式(19)保证了联合体实际出力与发电计划的偏差不会过大。
其中,常规水电机组约束包括以下约束:
(1)出力上下限约束
Figure BDA0004036776290000072
其中:
Figure BDA0004036776290000073
分别为第i台常规水电机组的出力的下限制及其上限制,MW;
Figure BDA0004036776290000079
为第i台常规水电机组在时段t的开停机状态变量,/>
Figure BDA0004036776290000074
1表示开机状态,0表示停机状态。
(2)发电流量约束
Figure BDA0004036776290000075
其中:
Figure BDA0004036776290000076
分别为第i台常规水电机组发电流量、发电流量上限制和发电流量下限制,单位为m3/s。
(3)振动区约束
Figure BDA0004036776290000077
其中:
Figure BDA0004036776290000078
分别为第i台常规水电机组的第k个振动区的出力上限、下限,单位为MW。
(4)出力爬坡约束
Figure BDA0004036776290000081
/>
其中:
Figure BDA0004036776290000082
为第i台常规水电机组的爬坡能力,单位为MW/h。
(5)机组开停机持续时间约束
Figure BDA0004036776290000083
其中:
Figure BDA0004036776290000084
为第i台常规水电机组在时段t的启动操作变量,1表示启动操作;/>
Figure BDA0004036776290000085
为第i台常规水电机组在时段t的停机操作变量,1表示停机操作;/>
Figure BDA0004036776290000086
分别第i台常规水电机组的最小开机、停机持续时段数;/>
Figure BDA0004036776290000087
代表调度期内的最大开机次数。
(6)机组出力波动限制约束
Figure BDA0004036776290000088
其中:te为常规水电机组在一轮出力升降过程中最高点或最低点需持续的最少时段数,te>1。
(7)机组发电水头约束
Figure BDA0004036776290000091
其中:Hi,t
Figure BDA0004036776290000092
分别为第i台机组在时段t的发电水头和水头损失,单位为m;/>
Figure BDA0004036776290000093
为时段t的上游水库的水位;/>
Figure BDA0004036776290000094
为时段t的下游水库的水位。
(8)机组动力特性关系约束
Figure BDA0004036776290000095
/>
其中:
Figure BDA0004036776290000096
为第i台常规水电机组的出力-水头-流量的关系函数。
抽水蓄能机组与常规水电机组相比,在启动灵活性、爬坡速度等性能上更具优势,能够更好匹配风电等随机性能源的快速波动。抽水蓄能机组在发电状态时,一般没有爬坡和最小开停机时间限制,除满足式(20)、式(21)、式(22)、式(24)、式(26)、式(27)的约束之外,在抽水状态下抽水蓄能机组约束还包括以下约束:
(1)抽水功率约束
Figure BDA0004036776290000097
其中:
Figure BDA0004036776290000098
分别为第j台抽水蓄能机组在时刻t的抽水功率、抽水功率下限和抽水功率上限,单位为MW;/>
Figure BDA0004036776290000099
为抽水蓄能机组的抽水状态变量,
Figure BDA0004036776290000101
1表示处于抽水状态,0表示不处于抽水状态。
(2)抽水流量约束
Figure BDA0004036776290000102
其中:
Figure BDA0004036776290000103
分为第j台抽水蓄能机组在时刻t的抽水流量、抽水流量上限和抽水流量下限,单位为m3/s。
(3)抽水功率特性曲线约束
Figure BDA0004036776290000104
其中:fj(·)为第j台抽水蓄能机组的抽水功率-水头-流量的关系函数;
Figure BDA0004036776290000105
为第j台抽水蓄能机组在时段t的抽水功率;Hj,t为第j台抽水蓄能机组在时段t的发电水头,单位为m;/>
Figure BDA0004036776290000106
为第j台抽水蓄能机组在时段t的发电流量。
(4)抽水和发电互斥约束
Figure BDA0004036776290000107
对于任意抽水蓄能机组,同一时刻不能出现同时发电和抽水的情况,发电与抽水状态互斥。
(5)启停次数约束
Figure BDA0004036776290000111
其中:
Figure BDA0004036776290000112
分别表示时段t+1、t运行的抽水蓄能机组台数;/>
Figure BDA0004036776290000113
为抽水蓄能机组总在调度期内的启停次数限制。
(6)机组状态切换约束
Figure BDA0004036776290000114
其中:
Figure BDA0004036776290000115
为第i台抽水蓄能机组在时段t的开停机状态变量,/>
Figure BDA0004036776290000116
1表示开机状态,0表示停机状态。
如此,通过上述约束保证了抽水蓄能机组在抽水和发电状态切换时必须经历停机状态,避免出现上一时段抽水,下一时段发电,或者上一时段发电,下一时段抽水的情况。
S2:求解混合式抽水蓄能电站与风电联合运行调度模型。
经分析,由公式(1)-(33)描述的调度模型是一个高维度、多变量、多约束的混合整数非线性约束规划(mixed integer nonlinear programming,简称:MINLP)问题,非线性约束包括式(3)、(5)、(12)、(13)、(15)、(22)、(25)、(27)、(30)、(31)。考虑到直接求解面临的求解效率低、初始解影响大等问题。
本实施例将原MINLP问题转化为MILP问题进行求解。MILP模型具有求解算法成熟、计算效率高、输出结果稳定等突出优势,已在水电领域得到广泛应用,其核心关键在于非线性约束的处理,建模技巧和线性化方法的好坏直接影响求解效率和精度。下面重点介绍非线性约束的处理方法和建模技巧,即S2具体包括以下步骤:
S21:线性化目标函数。
目标函数中式(3)和式(5)包含的max{·}函数导致了目标函数的非线性,难以直接求解。以式(3)为例,任意时刻t,本实施例通过引入2个0-1变量zt,m和3个连续变量wt,n实现转换,具体地:
Figure BDA0004036776290000121
Figure BDA0004036776290000122
wt,1≤zt,1 (36)
wt,2≤zt,1+zt,2 (37)
wt,3≤zt,2 (38)
Figure BDA0004036776290000123
/>
Figure BDA0004036776290000124
其中:wt,n实现函数的线性化,辅助变量zt,m限定wt,n的取值;bt,n为分位点,即联合体在时刻t的出力,本模型中,bt,1=0,bt,2=Pt plan,bt,3=bigM,bigM为大于Pt plan的实数;f(bt,n)为分位点对应的函数值,即分位点对应出力的高于发电计划的电量,f(bt,1)=f(bt,2)=0,f(bt,3)=bigM-Pt plan。因此,式(3)可以用式(34)-式(40)等价转换,式(5)的转换与此类似。
S22:线性化目标函数的约束条件。
具体的,S22包括:
S221:线性化机组出力波动限制约束。
由式(25)可知,常规水电机组的出力波动限制约束是和时段相关的非线性约束,需要考虑相邻多个时段的出力变化情况。本文通过引入功率上下调节指标变量、以及上下调节总次数变量进行线性化建模,描述如下:
Figure BDA0004036776290000131
αi,ti,t≤1 (42)
Figure BDA0004036776290000132
其中:αi,t∈{0,1}、βi,t∈{0,1}为第i台常规水电机组在时段t的功率向上调节指标变量、向下调节指标变量,αi,t=1表示时段t+1功率向下调节;βi,t=1表示时段t+1功率向上调节;当功率不发生变化时,αi,t=0或βi,t=0;MAXαβ表示功率调整(向上和向下)时段数上限。因此,为避免机组出力频繁出现相邻时段的上调和下调情况,保证机组每次调整之后,均可以稳定出力一定时段,功率上下调节指标还应满足以下公式约束:
Figure BDA0004036776290000141
如此,式(25)的出力波动限制非线性约束可用式(41)-式(44)等价转换。
S222:线性化运行工况关联与互斥约束。
由水电站融合改造的混合式抽水蓄能电站,机组类型包括常规水电机组和抽水蓄能机组。常规水电机组有发电、停机两种状态,抽水蓄能机组有发电、抽水、停机三种状态,如式(15)和(31)可知,电站和机组的运行工况都存在互斥关系,电站的运行状态又和机组的运行状态耦合关联,如何有效建模直接影响求解结果和效率。本实施例提出了以下基于状态变量解耦的建模和线性化方法,具体描述如下:
Figure BDA0004036776290000142
Figure BDA0004036776290000143
Figure BDA0004036776290000144
Figure BDA0004036776290000151
Figure BDA0004036776290000152
Figure BDA0004036776290000153
Figure BDA0004036776290000154
Figure BDA0004036776290000155
其中:
Figure BDA0004036776290000156
为抽水蓄能机组的状态指示变量,
Figure BDA0004036776290000157
表示第j台抽水蓄能机组在时段t处于发电状态,/>
Figure BDA0004036776290000158
表示第j台抽水蓄能机组在时段t处于抽水状态,则式(45)可约束抽水蓄能机组发电与抽水状态的互斥;
Figure BDA0004036776290000159
为电站的状态指示变量,/>
Figure BDA00040367762900001510
表示电站在时段t处于发电状态,/>
Figure BDA00040367762900001511
表示电站在时段t处于抽水状态,则式(46)可约束电站发电与抽水状态的互斥;/>
Figure BDA00040367762900001512
分别表示时段t的发电机组台数、抽水机组台数和停机机组台数,通过式(47)-(52)实现了机组状态与电站状态关联性的解耦和线性化建模。因此,式(15)和(31)的电站和机组状态非线性约束可用式(45)-(52)等价转换。
经过上述模型转换,将原MINLP模型转换为标准的MILP模型,求解流程如下。
Step1基础数据初始化,包括风电场预测出力、系统负荷、机组出力的上下限、抽水蓄能机组抽水功率上下限、水库库容上下限、振动区、爬坡能力、最小开停机持续时间、稳定出力持续时间等基础信息。
Step2约束线性化处理。采用S2中的建模技巧和线性化方法,对模型中涉及到的非线性约束进行线性化转化,构建标准的MILP模型。
Step3模型求解。CPLEX优化软件提供了Java接口jar包,可在Java环境中完成模型的建模编码和接口调用,实现MILP模型的高效求解。
Step4结果输出与分析。输出目标函数最优解、计算时间、机组各时段状态、机组出力、机组抽水功率、抽水和发电流量等结果信息。
实施例
以西南地区某梯级水电站的融合改造为参考构建应用实例,验证所提模型和求解方法的有效性。该水电站目前已有4台常规水电机组,拟融合改造新增2台抽水蓄能机组,形成混合式抽水蓄能电站,并与周边风电场联合调度运行。常规水电机组和抽水蓄能机组的特征参数如表1所示,以15min为1个时段,整个调度期共96个时段。需要说明的是,抽水蓄能机组在抽水工况时,一般运行在最优功率点附近,为提升抽水蓄能机组在抽水工况的调节能力,变速抽水蓄能技术是重要方向,以变速抽水蓄能机组进行建模,设置最小发电出力为额定出力的30%,最小抽水功率为额定抽水功率的70%。风电机组预测出力及联合体发电计划如图2所示。
表1混合式抽水蓄能电站机组特征参数
Figure BDA0004036776290000161
设置惩罚电价的目的是为了有效激励提高风电预测精度和联合体发电计划可执行性。不同电网和市场环境有不同的惩罚机制,本实施例中参考国外日前电力市场的典型经验公式,设置ω=0.4,如公式(53)-(54)所示。
考虑我国日前电力市场尚未全面开展,假设峰平谷上网电价如表2所示。
Figure BDA0004036776290000171
Figure BDA0004036776290000172
表2不同时段的上网电价
Figure BDA0004036776290000173
不考虑常规水电机组的启停成本,抽水蓄能机组启停成本为3000元/(台·次);最小开停机持续时段
Figure BDA0004036776290000174
调度期末上库水位偏差控制在δ=0.1m以内;常规水电机组稳定出力持续时段te=4;发电计划最大允许偏差ε为10%;上库径流采用典型日平均流量,为120m3/s。模型采用JAVA1.7版本编码,CPLEX软件版本为V12.10.0.0,设置收敛精度为1.0e-4,运行于DELL工作站,基本配置为4核8线程、8.0G内存、主频3.6GHz。
本实施例的结果分析如下:
(1)联合体出力与收益
联合体发电总收益为938.95万元,其中按发电计划执行的发电收益为944.83万元,高于发电计划的额外发电收益为0.85万元,开停机成本为2.1万元,低于发电计划的惩罚为4.64万元。通过混合式抽水蓄能电站和风电的互补协调,发电计划执行平均偏差为0.7%,满足所设置的偏差要求。
在凌晨00:00-07:00前风电出力高峰和系统负荷低谷阶段,2台抽水蓄能机组全部启动,并运行于抽水工况;随着08:00早高峰电负荷的增长及风电出力的降低,抽水蓄能机组由抽水工况转换为发电工况,特别在早晚负荷高峰时段,常规水电机组开始投入运行,保证了高峰用电负荷的同时,增加了高峰出力,提高了联合体整体发电收益,充分体现了常规水电机组融合改造后的“量调并重”的优势。
(2)水库水位变化情况
混合式抽水蓄能电站的上游水位如图3所示。可以看出,在夜晚时段,抽水蓄能机组抽水,水位上升,伴随白天抽水蓄能机组和常规水电机组的出力,水位逐渐下降,满足调度期末水位控制要求。
(3)抽水蓄能机组运行结果
抽水蓄能机组的抽水和发电功率如图4所示,图4中出力0KW以下的折线表示抽水功率,出力0KW以上的表示发电功率。在晚上负荷低谷阶段,2台抽水蓄能机组均处于抽水状态,有别于常规抽水蓄能机组以固定功率抽水不同,变速抽水蓄能机组在抽水状态也具备调节能力,能够很好地与波动性风电互补,跟踪负荷运行;白天发电状态下,充分发挥了灵活调节优势,一方面跟随风电的波动,避免了常规水电机组的频繁调节,另一方面,与常规水电机组一起共同支撑发电计划执行。同时也可以看出,抽水蓄能机组的运行工况、开停机持续时段、发电振动区等均满足所设置约束要求。
(4)常规水电机组运行结果
常规水电机组的出力过程如图5所示。从出力过程来看,在夜间00:00-07:00负荷低谷时段,常规水电机组保持停机状态,主要由抽水蓄能机组承担低谷调峰作用;在早高峰10:00和晚高峰19:00前后时段,当抽水蓄能机组出力基本满发还不能满足负荷需求时,常规水电机组启动,有效支撑晚高峰负荷。早上07:30左右,由于抽水蓄能机组在进行抽水到发电工况的转换,且必须经过停机工况过渡(式(33)),4#常规水电机组短时间开机出力,保证了负荷需求。此外,可以直观看出,机组出力有效避开了振动区,同时也满足所设置的开停机持续时间(1h)、出力波动限制约束(1h)、爬坡约束(60MW/时段)等要求,保证了机组的安全稳定运行,验证了文中约束处理的有效性。
本实施例提供的混合式抽水蓄能电站与风电联合运行调度方法的有益效果包括:
(1)混合式抽水蓄能电站包括常规水电机组和抽水蓄能机组,二者的运行特性差异较大。以机组为单元进行精细化建模,根据运行特性的不同分别设置运行约束,并通过引入机组状态变量实现运行工况的解耦与关联切换,实例表明,常规水电机组和抽水蓄能机组的运行均满足所设置约束要求,符合实际调度需求。
(2)混合式抽水蓄能电站具有典型的“量调并重”的优势,在夜晚负荷低谷时段,抽水蓄能机组处于抽水工况,有效支撑了低谷调峰并存储电量;在白天,抽水蓄能机组主要处于发电工况,承担着灵活快速调节的角色,一方面跟随风电出力保证了发电计划的执行,另一方面也避免了常规水电机组频繁调节带来的出力波动。
(3)以混合式抽水蓄能电站和风电所构成联合体整体的收益最大为目标函数,并通过引入奖惩电价机制,可以有效激励联合体提高风电的预测精度并制定合理的发电计划,从而提高了发电计划的可执行性。
以上,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (9)

1.一种混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,所述方法包括:
S1:建立混合式抽水蓄能电站与风电联合运行调度模型;
其中,所述调度模型以联合体整体收益最大为目标函数,用公式描述为:
maxF=F1+F2-F3-F4 (1)
其中,F为调度期的联合体整体运行效益,单位为元;F1、F2、F3、F4分别为联合体按发电计划的基础收益、实际出力高于发电计划的额外收益、开停机成本和实际出力低于发电计划的惩罚费用;
S2:求解混合式抽水蓄能电站与风电联合运行调度模型。
2.根据权利要求1所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,F1、F2、F3、F4计算公式分别如下:
Figure FDA0004036776280000011
Figure FDA0004036776280000012
Figure FDA0004036776280000013
Figure FDA0004036776280000014
其中,ct
Figure FDA0004036776280000015
分别为时段t的上网电价、出力高于发电计划的电价、出力低于发电计划的惩罚电价,单位为元/kWh;Pt plan为联合体发电计划,单位为MW;Pt wind为时段t的预测风电出力,单位为MW;/>
Figure FDA0004036776280000021
为时段t机组i的发电功率;/>
Figure FDA0004036776280000022
为时段t机组j的抽水功率,单位为MW;N、Npump分别为总机组台数和抽水蓄能机组的台数;con、coff分别为机组的启成本和停成本,单位为元/台次;/>
Figure FDA0004036776280000023
分别为t时段启的机组数量和停的机组数量,单位为台;t∈{1,2,…,T}为时段编号,T为调度期内总的时段数。
3.根据权利要求2所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,所述目标函数的约束条件包括电站约束、常规水电机组约束和抽水蓄能机组约束。
4.根据权利要求3所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,所述电站约束包括上下库水量平衡约束、库容大小约束、水位-库容关系约束、电站机组总数约束、电站发电与抽水互斥约束、调度期末水位控制约束和互补系统发电计划约束。
5.根据权利要求4所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,所述上下库水量平衡约束用公式描述为:
Figure FDA0004036776280000024
Figure FDA0004036776280000025
其中:Vt up、Vt down分别表示混合式抽水蓄能电站在时段t的上库库容和下库库容,单位为m3
Figure FDA0004036776280000026
分别为时段t上库的来水流量、下库来水流量、发电流量和抽水流量,单位为m3/s;Δt为计算时段步长,单位为h;
所述库容大小约束用公式描述为:
Figure FDA0004036776280000027
Figure FDA0004036776280000028
Figure FDA0004036776280000029
Figure FDA0004036776280000031
其中:
Figure FDA0004036776280000032
分别为抽水蓄能电站上库的最小库容和最大库容,单位为m3
Figure FDA0004036776280000033
分别为抽水蓄能电站下库的最小库容和最大库容,单位为m3;/>
Figure FDA0004036776280000034
分别为上库的初始库容和下库的初始库容;V1 up、V1 down分别为第1时段上库的库容和下库的库容;
所述水位-库容关系约束用公式描述为:
Figure FDA0004036776280000035
Figure FDA0004036776280000036
其中:
Figure FDA0004036776280000037
分别为水库上库的水位-库容关系曲线函数和下库的水位-库容关系曲线函数;/>
Figure FDA0004036776280000038
分别为时段t的上游水库的水位和下游水库的水位,单位为m;
所述电站机组总数约束用公式描述为:
Npump+Ntrad=N (14)
其中:Npump、Ntrad、N分别为抽水蓄能机组台数、常规水电机组台数以及总机组台数,单位为台;
所述电站发电与抽水互斥约束用公式描述为:
Figure FDA0004036776280000039
Figure FDA00040367762800000310
Figure FDA00040367762800000311
/>
其中:Pt gen、Pt pump分别为机组在时段t的发电功率和抽水功率,单位为MW;
Figure FDA00040367762800000312
为时段t常规水电机组i的发电功率;/>
Figure FDA00040367762800000313
为时段t抽水蓄能机组j的发电功率;
所述调度期末水位控制约束用公式描述为:
Figure FDA00040367762800000314
其中:
Figure FDA00040367762800000315
为调度期末的上游水库的水位;/>
Figure FDA00040367762800000316
为第1时段的上游水库的水位;δ为调度期内允许的水位变幅;
所述互补系统发电计划约束用公式描述为:
Figure FDA0004036776280000041
其中:ε是给定的发电计划偏差系数。
6.根据权利要求3所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,所述常规水电机组约束包括出力上下限约束、发电流量约束、振动区约束、出力爬坡约束、机组开停机持续时间约束、机组出力波动限制约束、机组发电水头约束和机组动力特性关系约束。
7.根据权利要求6所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,所述出力上下限约束用公式描述为:
Figure FDA0004036776280000042
其中:
Figure FDA0004036776280000043
分别为第i台常规水电机组的出力的下限制及其上限制,MW;/>
Figure FDA0004036776280000044
为第i台常规水电机组在时段t的开停机状态变量,/>
Figure FDA0004036776280000045
1表示开机状态,0表示停机状态;Ntrad为常规水电机组台数;Pt gen为机组在时段t的发电功率;
所述发电流量约束用公式描述为:
Figure FDA0004036776280000046
其中:
Figure FDA0004036776280000047
分别为第i台常规水电机组发电流量、发电流量上限制和发电流量下限制,单位为m3/s;
所述振动区约束用公式描述为:
Figure FDA0004036776280000048
其中:
Figure FDA0004036776280000051
分别为第i台常规水电机组的第k个振动区的出力上限、下限,单位为MW;
所述出力爬坡约束用公式描述为:
Figure FDA0004036776280000052
其中:ΔPi gen为第i台常规水电机组的爬坡能力,单位为MW/h;
所述机组开停机持续时间约束用公式描述为:
Figure FDA0004036776280000053
Figure FDA0004036776280000054
Figure FDA0004036776280000055
Figure FDA0004036776280000056
Figure FDA0004036776280000057
其中:
Figure FDA0004036776280000058
为第i台常规水电机组在时段t的启动操作变量,1表示启动操作;/>
Figure FDA0004036776280000059
为第i台常规水电机组在时段t的停机操作变量,1表示停机操作;Ti on、Ti off分别第i台常规水电机组的最小开机、停机持续时段数;/>
Figure FDA00040367762800000510
代表调度期内的最大开机次数;t∈{1,2,…,T}为时段编号,T为调度期内总的时段数;
所述机组出力波动限制约束用公式描述为:
Figure FDA00040367762800000511
其中:te为常规水电机组在一轮出力升降过程中最高点或最低点需持续的最少时段数,te>1;
所述机组发电水头约束用公式描述为:
Figure FDA0004036776280000061
其中:Hi,t
Figure FDA0004036776280000062
分别为第i台机组在时段t的发电水头和水头损失,单位为m;/>
Figure FDA0004036776280000063
为时段t的上游水库的水位;/>
Figure FDA0004036776280000064
为时段t的下游水库的水位;
所述机组动力特性关系约束用公式描述为:
Figure FDA0004036776280000065
其中:fi NHQ(·)为第i台常规水电机组的出力-水头-流量的关系函数。
8.根据权利要求7所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,在发电状态下所述抽水蓄能机组约束包括式(20)、式(21)、式(22)、式(24)、式(26)和式(27),在抽水状态下所述抽水蓄能机组约束包括抽水功率约束、抽水流量约束、抽水功率特性曲线约束、抽水和发电互斥约束、启停次数约束和机组状态切换约束。
9.根据权利要求8所述的混合式抽水蓄能电站与风电联合运行调度方法,其特征在于,所述抽水功率约束用公式描述为:
Figure FDA0004036776280000066
其中:
Figure FDA0004036776280000067
分别为第j台抽水蓄能机组在时刻t的抽水功率、抽水功率下限和抽水功率上限,单位为MW;/>
Figure FDA0004036776280000068
为抽水蓄能机组的抽水状态变量,
Figure FDA0004036776280000069
1表示处于抽水状态,0表示不处于抽水状态;Npump为抽水蓄能机组台数;
所述抽水流量约束用公式描述为:
Figure FDA0004036776280000071
其中:
Figure FDA0004036776280000072
分为第j台抽水蓄能机组在时刻t的抽水流量、抽水流量上限和抽水流量下限,单位为m3/s;
所述抽水功率特性曲线约束用公式描述为:
Figure FDA0004036776280000073
其中:fj(·)为第j台抽水蓄能机组的抽水功率-水头-流量的关系函数;
Figure FDA0004036776280000074
为第j台抽水蓄能机组在时段t的抽水功率;Hj,t为第j台抽水蓄能机组在时段t的发电水头,单位为m;/>
Figure FDA0004036776280000075
为第j台抽水蓄能机组在时段t的发电流量;
所述抽水和发电互斥约束用公式描述为:
Figure FDA0004036776280000076
其中:
Figure FDA0004036776280000077
为时段t抽水蓄能机组j的发电功率;t∈{1,2,…,T}为时段编号,T为调度期内总的时段数;
所述启停次数约束用公式描述为:
Figure FDA0004036776280000078
其中:
Figure FDA0004036776280000079
分别表示时段t+1、t运行的抽水蓄能机组台数;/>
Figure FDA00040367762800000710
为抽水蓄能机组总在调度期内的启停次数限制;
所述机组状态切换约束用公式描述为:
Figure FDA0004036776280000081
其中:
Figure FDA0004036776280000082
为第i台抽水蓄能机组在时段t的开停机状态变量,/>
Figure FDA0004036776280000083
1表示开机状态,0表示停机状态。/>
CN202310008506.6A 2023-01-04 2023-01-04 一种混合式抽水蓄能电站与风电联合运行调度方法 Pending CN115879734A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310008506.6A CN115879734A (zh) 2023-01-04 2023-01-04 一种混合式抽水蓄能电站与风电联合运行调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310008506.6A CN115879734A (zh) 2023-01-04 2023-01-04 一种混合式抽水蓄能电站与风电联合运行调度方法

Publications (1)

Publication Number Publication Date
CN115879734A true CN115879734A (zh) 2023-03-31

Family

ID=85758037

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310008506.6A Pending CN115879734A (zh) 2023-01-04 2023-01-04 一种混合式抽水蓄能电站与风电联合运行调度方法

Country Status (1)

Country Link
CN (1) CN115879734A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116231767A (zh) * 2023-05-11 2023-06-06 南方电网调峰调频发电有限公司储能科研院 一种梯级水电站的多能互补调度方法和系统
CN117114330A (zh) * 2023-08-31 2023-11-24 湖北清江水电开发有限责任公司 基于虚拟电厂的抽水蓄能和梯级水电联合调度方法
CN117791663A (zh) * 2024-02-27 2024-03-29 华北电力大学 一种梯级混合式抽水蓄能电站调控方法、系统及电子设备

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116231767A (zh) * 2023-05-11 2023-06-06 南方电网调峰调频发电有限公司储能科研院 一种梯级水电站的多能互补调度方法和系统
CN117114330A (zh) * 2023-08-31 2023-11-24 湖北清江水电开发有限责任公司 基于虚拟电厂的抽水蓄能和梯级水电联合调度方法
CN117791663A (zh) * 2024-02-27 2024-03-29 华北电力大学 一种梯级混合式抽水蓄能电站调控方法、系统及电子设备

Similar Documents

Publication Publication Date Title
CN115879734A (zh) 一种混合式抽水蓄能电站与风电联合运行调度方法
CN109284878B (zh) 一种考虑风电、核电、抽水蓄能协调的多源优化调度方法
CN107153885B (zh) 考虑火电机组深度调峰的实时发电计划优化方法
CN112103943B (zh) 日前电力现货市场出清的安全校核方法、装置及存储介质
CN112952805A (zh) 一种考虑柔性氢需求的电氢能源系统调度方法
CN107276122B (zh) 适应大规模可再生能源并网的调峰资源调用决策方法
CN111431213A (zh) 激励风电场和抽水蓄能电站联合运行的厂网协调方法及其联合调度方法
CN111882131B (zh) 一种基于双工况最优效率跟踪的变速海水抽蓄电站容量优化方法
Zhang et al. An ultra-short-term scheduling model for cascade hydropower regulated by multilevel dispatch centers suppressing wind power volatility
CN112909933A (zh) 现货市场环境下含抽水蓄能机组的日内滚动优化调度方法
CN113937819A (zh) 一种多能源短期优化调度方法
CN112910013A (zh) 考虑“深度调峰消纳-煤耗”联合约束的机组优化调度方法
Kloess Electric storage technologies for the future power system—An economic assessment
CN110854933A (zh) 利用灵活性资源的月度机组组合优化方法
CN114254937A (zh) 以机组为调度单元的梯级水电及光伏短期互补调度方法和系统
Rizvi et al. Optimal scheduling of virtual power plants utilizing wind power and electric vehicles
CN112688328B (zh) 一种交直流混合微电网时间协调能量优化配置方法
Tian et al. Research on day-ahead optimal dispatch considering the mutual constraint of pumped storage unit
Zhu et al. Dynamic economic dispatching strategy based on multi-time-scale complementarity of various heterogeneous energy
Xu et al. A two-layer multi-energy management system for microgrids with solar, wind, and geothermal renewable energy
CN114465226A (zh) 一种电力系统多级备用获取联合优化模型的建立方法
Hao et al. Research on power trade potential and power balance between Lancang-Mekong countries and southern China under long-term operation simulation
CN114204599A (zh) 一种抑制配电网电压波动的分布式储能布点运行控制方法
CN113629781A (zh) 一种基于empc的多vpp日前日内协调优化方法及装置
CN113077096B (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