CN115640982B - 基于抽水蓄能优先调节的多能互补系统日前优化调度方法 - Google Patents

基于抽水蓄能优先调节的多能互补系统日前优化调度方法 Download PDF

Info

Publication number
CN115640982B
CN115640982B CN202211442469.1A CN202211442469A CN115640982B CN 115640982 B CN115640982 B CN 115640982B CN 202211442469 A CN202211442469 A CN 202211442469A CN 115640982 B CN115640982 B CN 115640982B
Authority
CN
China
Prior art keywords
constraint
thermal power
complementary system
energy
day
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
CN202211442469.1A
Other languages
English (en)
Other versions
CN115640982A (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202211442469.1A priority Critical patent/CN115640982B/zh
Publication of CN115640982A publication Critical patent/CN115640982A/zh
Application granted granted Critical
Publication of CN115640982B publication Critical patent/CN115640982B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明属于电力系统优化运行技术领域,公开了基于抽水蓄能优先调节的多能互补系统日前优化调度方法,包括建立风‑光‑水‑火‑抽水蓄能多能互补系统的日前优化调度模型,将运行成本最小和碳排放量最低作为目标函数,多能互补系统中各能源调度的顺序依次为:风电与光伏、抽水蓄能、常规水电站、火电厂;基于各能源的基础数据,采用多目标混合粒子群算法对日前优化调度模型进行求解得到Pareto解集,基于隶属度函数选取一个最优折中解作为日前优化调度方案。本发明可降低多能互补系统的运行成本和碳排放量,释放常规水电和火电机组调节压力,能有效地应用于“水风光火储一体化系统”的运行管理中。

Description

基于抽水蓄能优先调节的多能互补系统日前优化调度方法
技术领域
本发明属于电力系统优化运行技术领域,更具体地,涉及一种基于抽水蓄能优先调节的多能互补系统日前优化调度方法。
背景技术
目前,火力发电在中国电力系统中仍占主要地位,若仅依靠电力系统中火电机组的灵活性对可再生能源进行调节,必将造成火电机组频繁启停,严重威胁电网运行的安全性与经济性,且不利于低碳减排工作。由于水电技术成熟、具有优质的调节能力,近年来已有水火联合调度的一些尝试,但水资源不确定性和有限性使水电不足以满足人们日益增长的调节需求。破解上述困境的一种方法是利用抽水蓄能电站灵活的存储特性,使能量能够从整个系统中重新分配。因此含有抽水蓄能的风-光-水-火-抽水蓄能多能互补系统协同优化调度研究备受关注。
归纳现有技术,我们发现风-光-水-火-抽水蓄能多能互补系统协同优化调度研究主要存在以下不足:1)研究对象方面:很多研究关注风-水-火、风-光-水-火等系统,针对风-光-水-火-抽蓄多能互补系统优化运行研究较少,不完全适用于我国能源发展现状和未来趋势;2)调度目标方面:虽然部分研究关注了风-光-水-火-抽蓄多能互补系统协同优化调度问题,但调度目标仍限于经济调度,考虑碳排放的研究较少;3)模型求解方面:目前对常规水电和抽水蓄能同时存在的多能互补调度问题研究较少,少量研究将抽水蓄能考虑为梯级水电站调度过程的一个分支,通过水力联系使两者形成整体;4)子系统特性分析方面:许多研究仅根据目标函数值来评判多能互补系统优劣、变化,对子系统运行、调节强度缺乏评价指标和定量分析。
发明内容
为克服上述现有技术存在的不足,本发明提供一种基于抽水蓄能优先调节的多能互补系统日前优化调度方法,能够降低风-光-水-火-抽水蓄能多能互补系统的运行成本和碳排放量,释放常规水电和火电机组调节压力。
本发明提供一种基于抽水蓄能优先调节的多能互补系统日前优化调度方法,包括以下步骤:
建立风-光-水-火-抽水蓄能多能互补系统的日前优化调度模型;所述日前优化调度模型包括目标函数和约束条件;将多能互补系统的运行成本最小和所述多能互补系统的碳排放量最低作为所述目标函数;所述多能互补系统中各能源调度的顺序依次为:风电与光伏、抽水蓄能、常规水电站、火电厂;
获取各能源的基础数据;
基于所述各能源的基础数据,采用多目标混合粒子群算法对所述日前优化调度模型进行求解,求解的过程中基于所述约束条件进行约束处理,得到Pareto解集;
基于隶属度函数,从所述Pareto解集中选取一个最优折中解,并作为日前优化调度方案。
优选的,所述目标函数包括第一目标函数和第二目标函数;
所述第一目标函数如下:
Figure DEST_PATH_IMAGE001
Figure DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE003
式中,
Figure DEST_PATH_IMAGE004
为第一目标函数,对应多能互补系统的运行成本;N为火电厂的火电机组 的数量;T为时段数目;
Figure DEST_PATH_IMAGE005
为第j台火电机组在i时段的启停状态,
Figure DEST_PATH_IMAGE006
表示机组开启,
Figure DEST_PATH_IMAGE007
表示机组关闭;
Figure DEST_PATH_IMAGE008
为第j台火电机组在i-1时段的启停状态;
Figure DEST_PATH_IMAGE009
为第j台火电机组在 i时段的启动成本;
Figure DEST_PATH_IMAGE010
为第j台火电机组在i时段的出力;
Figure DEST_PATH_IMAGE011
Figure DEST_PATH_IMAGE012
Figure DEST_PATH_IMAGE013
分别为第j台火电机组 的运行成本系数;
Figure DEST_PATH_IMAGE014
表示第j台火电机组的火电热启动成本;
Figure DEST_PATH_IMAGE015
表示第j台火电机组 的火电冷启动成本;
Figure DEST_PATH_IMAGE016
为第j台火电机组的最小停机时长,
Figure DEST_PATH_IMAGE017
为第j台火电机组在i时段 的累计停机时长,
Figure DEST_PATH_IMAGE018
为第j台火电机组的冷启动时长;
所述第二目标函数如下:
Figure DEST_PATH_IMAGE019
式中,
Figure DEST_PATH_IMAGE020
为第二目标函数,对应多能互补系统的碳排放量;
Figure DEST_PATH_IMAGE021
Figure DEST_PATH_IMAGE022
Figure DEST_PATH_IMAGE023
分别为第j 台火电机组的碳排放系数。
优选的,所述日前优化调度模型还包括评价指标,所述评价指标用于衡量所述多能互补系统的运行状况和调节强度。
优选的,所述评价指标包括:火电厂运行小时数、常规水电站调节里程和抽水蓄能电站启停次数。
优选的,所述常规水电站调节里程表示为:
Figure DEST_PATH_IMAGE024
式中,
Figure DEST_PATH_IMAGE025
为常规水电站调节里程,
Figure DEST_PATH_IMAGE026
为水电机组的额定功率,T为时段数目, 为水电机组在i时刻的出力,
Figure DEST_PATH_IMAGE027
为水电机组在i+1时刻的出力。
优选的,所述各能源的基础数据包括风电出力、光伏出力、多能互补系统的负荷、常规水电站的水库入流量。
优选的,所述约束条件包括:所述多能互补系统的系统约束条件、新能源约束条件、水电站约束条件、火电厂约束条件和抽水蓄能电站约束条件;
所述多能互补系统的系统约束条件包括系统功率平衡约束和系统旋转备用约束;所述新能源约束条件包括风电出力约束和光伏出力约束;所述水电站约束条件包括水电站的出力约束、水库库容约束和水库流量约束;所述火电厂约束条件包括火电机组出力约束、爬坡速率约束和最小启停时间约束;所述抽水蓄能电站约束条件包括抽水蓄能电站的发电功率约束、抽水功率约束、工况约束和库容约束。
优选的,对所述日前优化调度模型进行求解时,将所述日前优化调度模型划分为四个子模型进行求解,所述四个子模型分别为风电与光伏子模型、抽水蓄能子模型、常规水电站子模型、火电厂子模型。
优选的,采用多目标混合粒子群算法对所述日前优化调度模型进行求解包括如下步骤:
步骤1、设定风电、光伏全额入网,计算所述多能互补系统除去新能源后的剩余负荷;
步骤2、建立所述多能互补系统的决策变量矩阵并设置参数;所述决策变量矩阵包括抽水蓄能的抽水功率变量和发电功率变量、常规水电站的流量变量、火电厂的机组运行状态变量;所述参数包括种群大小、迭代次数、学习因子、惯性权重;
步骤3、种群初始化;
步骤4、进行约束处理;
步骤5、将各参数超出约束值构建为惩罚函数,并添加至所述目标函数;
步骤6、计算所述目标函数的适应度值;
步骤7、根据所述适应度值确定非支配解,更新粒子速度和位置,并继续执行步骤4至步骤7,直至达到最大迭代次数时输出Pareto解集。
优选的,所述步骤4中进行约束处理包括如下子步骤:
步骤4.1、针对抽水蓄能的库容约束、抽水功率约束、发电功率约束进行处理;
步骤4.2、针对水电站的水库库容约束、出力约束进行处理;
步骤4.3、针对火电厂的爬坡速率约束、最小启停时间约束、出力约束进行处理;
步骤4.4、通过系统功率平衡约束将各个子模型进行耦合。
本发明中提供的一个或多个技术方案,至少具有如下技术效果或优点:
在本发明中,构建了风-光-水-火-抽水蓄能多能互补系统的日前优化调度模型,将多能互补系统的运行成本最小和碳排放量最低作为目标函数,多能互补系统中各能源调度的顺序依次为:风电与光伏、抽水蓄能、常规水电站、火电厂,即本发明提出的是一种在常规水电站和抽水蓄能同时存在的情况下抽水蓄能优先调节的调度策略,本发明在新能源外优先调节抽水蓄能,能够释放常规水电和火电机组调节压力,减少因机组频繁调节而造成的损耗,提高多能互补系统的经济性、降低碳排放量。本发明基于各能源的基础数据,采用多目标混合粒子群算法(MOHPSO)对日前优化调度模型进行求解,求解的过程中基于日前优化调度模型的约束条件进行约束处理得到Pareto解集;最后基于隶属度函数,从Pareto解集中选取一个最优折中解,并作为日前优化调度方案。此外,本发明还提出了将日前优化调度模型分解为4个子模型的求解方法,并引入了常规水电站调节里程等衡量子系统运行情况和调节强度的指标,有助于对子系统的运行、调节强度进行评价和分析。本发明提供的调度方法能有效地应用于“水风光火储一体化系统”和“综合能源基地”的运行管理中。
附图说明
图1为本发明实施例提供的一种基于抽水蓄能优先调节的多能互补系统日前优化调度方法的原理流程图;
图2为风电出力值、光伏出力值、多能互补系统剩余负荷;
图3为3个实验组(方案S1、方案S2、方案S3)的Pareto解集;
图4为方案S1求出的20组解归一化隶属度值;
图5为3个实验组对应的多能互补系统出力计划;
图6为3个实验组对应的常规水电站下泄流量;
图7为3个实验组对应的常规水电站水库库容;
图8为3个实验组对应的目标函数值和各子系统评价指标值;其中,图8中的(a)为系统运行成本,图8中的(b)为碳排放量,图8中的(c)为常规水电调节里程,图8中的(d)为火电机组运行时间。
具体实施方式
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
本实施例提供一种基于抽水蓄能优先调节的多能互补系统日前优化调度方法,参见图1,包括以下步骤:
步骤1:建立风-光-水-火-抽水蓄能多能互补系统的日前优化调度模型。
所述多能互补系统中各能源调度的顺序依次为:风电与光伏、抽水蓄能、常规水电站、火电厂,即本发明建立的是基于抽水蓄能优先调节的风-光-水-火-抽水蓄能多能互补系统日前优化调度模型。
本发明为响应节能发电调度政策并促进风电和太阳能发展,多能互补系统中各能源调度顺序依次为:不可调可再生能源(风电和光伏)、可调可再生能源(抽水蓄能电站与常规水电站)、可调不可再生能源(火电厂)。
所述日前优化调度模型包括目标函数和约束条件。
(1)目标函数。
将多能互补系统的运行成本最小和所述多能互补系统的碳排放量最低作为所述目标函数。
具体的,所述目标函数包括第一目标函数和第二目标函数。
所述第一目标函数如下:
Figure DEST_PATH_IMAGE028
Figure DEST_PATH_IMAGE029
Figure DEST_PATH_IMAGE030
式中,
Figure DEST_PATH_IMAGE031
为第一目标函数,对应多能互补系统的运行成本;N为火电厂的火电机组 的数量;T为时段数目;
Figure DEST_PATH_IMAGE032
为第j台火电机组在i时段的启停状态,
Figure DEST_PATH_IMAGE033
表示机组开启,
Figure DEST_PATH_IMAGE034
表示机组关闭;
Figure DEST_PATH_IMAGE035
为第j台火电机组在i-1时段的启停状态;
Figure DEST_PATH_IMAGE036
为第j台火电机组在i 时段的启动成本;
Figure DEST_PATH_IMAGE037
为第j台火电机组在i时段的出力;
Figure DEST_PATH_IMAGE038
Figure DEST_PATH_IMAGE039
Figure DEST_PATH_IMAGE040
分别为第j台火电机组的 运行成本系数;
Figure 214581DEST_PATH_IMAGE014
表示第j台火电机组的火电热启动成本;
Figure DEST_PATH_IMAGE041
表示第j台火电机组的 火电冷启动成本;
Figure DEST_PATH_IMAGE042
为第j台火电机组的最小停机时长,
Figure DEST_PATH_IMAGE043
为第j台火电机组在i时段的 累计停机时长,
Figure DEST_PATH_IMAGE044
为第j台火电机组的冷启动时长。
由于多能互补系统中可再生能源和常规水电、抽水蓄能均不消耗燃料,且常规水电站和抽蓄的成本主要为建设成本,简化起见,在短期调度时其运行成本可忽略,故多能互补系统运行成本主要为火电机组的燃煤成本和启停成本。
所述第二目标函数如下:
Figure DEST_PATH_IMAGE045
式中,
Figure DEST_PATH_IMAGE046
为第二目标函数,对应多能互补系统的碳排放量;
Figure DEST_PATH_IMAGE047
Figure DEST_PATH_IMAGE048
Figure DEST_PATH_IMAGE049
分别为第j台 火电机组的碳排放系数。
可再生能源、常规水电和抽水蓄能均为清洁能源,多能互补系统的碳排放主要源于火电厂燃烧燃料。为减缓温室效应和实现碳中和目标,本发明考虑多能互补系统的碳排放量最低,即火电厂碳排放量最低。
(2)约束条件。
所述约束条件包括:所述多能互补系统的系统约束条件、新能源约束条件、水电站约束条件、火电厂约束条件和抽水蓄能电站约束条件。
(2.1)多能互补系统的系统约束条件,包括系统功率平衡约束和系统旋转备用约束。
由于电能不具有大量储存的特性,故多能互补系统应实现供电与负荷需求之间的平衡。不考虑系统网损,系统功率平衡方程为:
Figure DEST_PATH_IMAGE050
其中,i为时间尺度(h);
Figure DEST_PATH_IMAGE051
为i时段风电场实际调度出力(MW);
Figure DEST_PATH_IMAGE052
为i时刻光伏 电站实际调度出力(MW);
Figure DEST_PATH_IMAGE053
为i时段水电站出力(MW);
Figure DEST_PATH_IMAGE054
为i时段火电厂出力(MW);
Figure DEST_PATH_IMAGE055
为 i时段抽水蓄能电站出力(MW),
Figure DEST_PATH_IMAGE056
表示放水发电,
Figure DEST_PATH_IMAGE057
表示抽水储能;
Figure DEST_PATH_IMAGE058
为i时段多 能互补系统的负荷值(MW)。
旋转备用可分为上旋转和下旋转,分别表示发电输出和最大/最小可调功率输出之间的差异。
Figure DEST_PATH_IMAGE059
Figure DEST_PATH_IMAGE060
其中,
Figure DEST_PATH_IMAGE061
为旋转备用容量(MW);
Figure DEST_PATH_IMAGE062
Figure DEST_PATH_IMAGE063
分别为火电厂最小、最大出力(MW)。
(2.2)新能源约束条件,包括风电出力约束和光伏出力约束。
新能源(风电和光伏)的出力限制约束如下:
Figure DEST_PATH_IMAGE064
Figure DEST_PATH_IMAGE065
其中,
Figure DEST_PATH_IMAGE066
Figure DEST_PATH_IMAGE067
分别为风电场、光伏电站最大出力(MW)。
(2.3)水电站约束条件,包括水电站的出力约束、水库库容约束和水库流量约束。
即设置水电站出力限制和水量平衡约束。水电站出力与水轮机流量和水库库容呈非线性关系,可将其表示为二元二次多项式:
Figure DEST_PATH_IMAGE068
Figure DEST_PATH_IMAGE069
其中,
Figure DEST_PATH_IMAGE070
为i时段水库库容(104m3);
Figure DEST_PATH_IMAGE071
为i时段水电站下泄流量(104m3);
Figure DEST_PATH_IMAGE072
为水力 发电系数,n取1~6。
水电站出力限制约束如下:
Figure DEST_PATH_IMAGE073
其中,
Figure DEST_PATH_IMAGE074
Figure DEST_PATH_IMAGE075
分别为i时段水电站最大、最小出力。
水电站水量平衡约束如下:
(a)水库库容约束:
Figure DEST_PATH_IMAGE076
Figure DEST_PATH_IMAGE077
Figure DEST_PATH_IMAGE078
Figure DEST_PATH_IMAGE079
其中,
Figure DEST_PATH_IMAGE080
Figure DEST_PATH_IMAGE081
分别为i时段水电站最大、最小库容(104m3),
Figure DEST_PATH_IMAGE082
为i时段水库入流 量(104m3),
Figure DEST_PATH_IMAGE083
为i时段水库溢流量(104m3),
Figure DEST_PATH_IMAGE084
Figure DEST_PATH_IMAGE085
分别为水电站初始、末库容(104m3),
Figure DEST_PATH_IMAGE086
为水电站始末库容比例系数。
(b)水库流量约束:
Figure DEST_PATH_IMAGE087
其中,
Figure DEST_PATH_IMAGE088
Figure DEST_PATH_IMAGE089
分别为水电站最小、最大流量(104m3)。
(2.4)火电厂约束条件,包括火电机组出力约束、爬坡速率约束和最小启停时间约束。
(a)火电机组出力限制约束:
Figure DEST_PATH_IMAGE090
其中,
Figure DEST_PATH_IMAGE091
Figure DEST_PATH_IMAGE092
分别为第j台火电机组出力最小、最大值(MW)。
(b)爬坡速率约束:
Figure DEST_PATH_IMAGE093
其中,
Figure DEST_PATH_IMAGE094
Figure DEST_PATH_IMAGE095
分别为第j台火电机组最大增、减出力速率(MW/h),
Figure DEST_PATH_IMAGE096
为时间间 隔。
(c)最小启停时间约束:
Figure DEST_PATH_IMAGE097
其中,
Figure DEST_PATH_IMAGE098
为第j台机组最小运行时长(h),
Figure DEST_PATH_IMAGE099
为第j台机组至i时刻运行时长 (h),
Figure DEST_PATH_IMAGE100
为第j台机组最小停机时长(h),
Figure DEST_PATH_IMAGE101
为第j台机组至i时刻停机时长(h)。
(2.5)抽水蓄能电站约束条件,包括抽水蓄能电站的发电功率约束、抽水功率约束、工况约束和库容约束。
将上水库库容转化为电能(MWh),下水库库容考虑为无穷大(∞),后文所提抽水蓄能水库均指上水库。
(a)放水发电功率约束(
Figure DEST_PATH_IMAGE102
):
Figure DEST_PATH_IMAGE103
Figure DEST_PATH_IMAGE104
Figure DEST_PATH_IMAGE105
Figure DEST_PATH_IMAGE106
其中,
Figure DEST_PATH_IMAGE107
Figure DEST_PATH_IMAGE108
分别为抽水蓄能电站最大、最小发电功率(MW),
Figure DEST_PATH_IMAGE109
为i时段抽水 蓄能电站上水库库容(MWh),
Figure DEST_PATH_IMAGE110
为时间间隔(h),
Figure DEST_PATH_IMAGE111
为抽水蓄能电站抽水储能效率,
Figure DEST_PATH_IMAGE112
为抽 水蓄能电站放水发电效率,
Figure DEST_PATH_IMAGE113
Figure DEST_PATH_IMAGE114
分别为抽水蓄能电站在第i时段的抽水功率、发电功 率(MW),
Figure DEST_PATH_IMAGE115
为抽水蓄能水库最大储能量(MWh)。
(b)抽水储能功率约束(
Figure DEST_PATH_IMAGE116
):
Figure DEST_PATH_IMAGE117
Figure DEST_PATH_IMAGE118
其中,
Figure DEST_PATH_IMAGE119
Figure DEST_PATH_IMAGE120
为抽水蓄能电站最大、最小抽水功率(MW)。
(c)工况限制约束:
抽水蓄能有放水发电和抽水储能两种工况,两种工况不同时发生。
Figure DEST_PATH_IMAGE121
(d)库容限制约束:
对于日调节抽水蓄能电站,为保证抽水蓄能下一个调度周期的正常运行,抽水蓄能初始库容和末库容应保持在一定范围。
Figure DEST_PATH_IMAGE122
其中,γ为始末库容比例系数;
Figure DEST_PATH_IMAGE123
Figure DEST_PATH_IMAGE124
为抽水蓄能电站初始、末库容(MWh)。
优选的方案中,本发明提出的所述日前优化调度模型还包括评价指标,所述评价指标用于衡量所述多能互补系统的运行状况和调节强度。
所述评价指标包括:火电厂运行小时数、常规水电站调节里程和抽水蓄能电站启停次数。下面对各个评价指标进行说明。
(a)火电厂运行小时数(
Figure DEST_PATH_IMAGE125
);
Figure DEST_PATH_IMAGE126
其中,TTO为火电厂运行小时数(h)。火电厂运行小时数越大意味着火电运行时间越长、机组损耗也越大,该指标越大意味着火电机组参与调节的时间就越多,就越不利。
(b)常规水电站调节里程(
Figure DEST_PATH_IMAGE127
);
Figure DEST_PATH_IMAGE128
其中,
Figure DEST_PATH_IMAGE129
为常规水电站调节里程(MW2),
Figure 687980DEST_PATH_IMAGE129
越大表示机组调节强度越大,对机 组造成的磨损越大;
Figure DEST_PATH_IMAGE130
为水电机组的额定功率,T为时段数目,
Figure DEST_PATH_IMAGE131
为水电机组在i时刻的 出力,
Figure DEST_PATH_IMAGE132
为水电机组在i+1时刻的出力。
(c)抽水蓄能电站启停次数(
Figure DEST_PATH_IMAGE133
)。
Figure DEST_PATH_IMAGE134
Figure DEST_PATH_IMAGE135
Figure DEST_PATH_IMAGE136
Figure DEST_PATH_IMAGE137
其中,
Figure DEST_PATH_IMAGE138
为抽水蓄能电站启停次数,抽水蓄能机组过于频繁启停会对抽水蓄能 机组造成损害,也会增加水头损失,一般日调节抽水蓄能电站设计运行强度为2抽1发,即启 停3次;
Figure DEST_PATH_IMAGE139
Figure DEST_PATH_IMAGE140
分别为抽水蓄能发电、抽水状态,1/0代表开/关。
Figure DEST_PATH_IMAGE141
Figure 169515DEST_PATH_IMAGE140
的i为0表示一 天(24h)中1:00的前一时刻,i为25表示一天中24:00的后一时刻。
步骤2:获取各能源的基础数据;基于所述各能源的基础数据,采用多目标混合粒子群算法对所述日前优化调度模型进行求解,求解的过程中基于所述约束条件进行约束处理,得到Pareto解集。
其中,所述各能源的基础数据包括风电出力、光伏出力、多能互补系统的负荷、常规水电站的水库入流量。
采用多目标混合粒子群算法对所述日前优化调度模型进行求解包括如下步骤:
Step 1:设定风电、光伏全额入网,计算所述多能互补系统除去新能源后的剩余负荷。
本发明提供的抽水蓄能优先调节模式是抽水蓄能优先参与调节净负荷(风光并网后的剩余负荷)、水电与火电后参与调节的运行模式。风电、光伏全额入网后多能互补系统剩余负荷为:
Figure DEST_PATH_IMAGE142
其中,
Figure DEST_PATH_IMAGE143
为风电、光伏全额入网后多能互补系统的剩余负荷,
Figure DEST_PATH_IMAGE144
为多能互补系统 的负荷,
Figure DEST_PATH_IMAGE145
为风电场出力,
Figure DEST_PATH_IMAGE146
为光伏电站出力。
Step 2:建立所述多能互补系统的决策变量矩阵并设置参数;所述决策变量矩阵包括抽水蓄能的抽水功率变量和发电功率变量、常规水电站的流量变量、火电厂的机组运行状态变量(二进制型决策变量);所述参数包括种群大小、迭代次数、学习因子、惯性权重等。
所述决策变量矩阵可表示为:
Figure DEST_PATH_IMAGE147
例如,
Figure DEST_PATH_IMAGE148
为抽水蓄能在时段3的功率(放水发电功率或抽水储能功率),
Figure DEST_PATH_IMAGE149
为水电 站在时段3的下泻流量,
Figure DEST_PATH_IMAGE150
为第1台火电机组在时段3的启停状态。
Step 3:种群初始化。
Step 3.1:在约束范围内采用混沌算法原理生成初始种群;
Step 3.2:设置迭代计数器为n=1。
Step 4:进行约束处理。
Step 4.1:针对抽水蓄能的库容约束、抽水功率约束、发电功率约束进行处理。
设定抽水蓄能优先参与调节,并定义抽水蓄能调节后的剩余负荷为净负荷(分配给常规水电和火电机组的负荷),具体处理步骤为:
Step 4.1.1:根据约束条件中放水发电功率约束求解得到抽水蓄能库容E;
Step 4.1.2:根据约束条件中库容限制约束求解得出E的违反程度ΔE;
Figure DEST_PATH_IMAGE151
Step 4.1.3:若|ΔE|<δ,则无需处理,否则需回溯处理,具体为:
1)ΔE<-δ,计算平均违反程度,并调整抽水储能和放水发电功率:
Figure DEST_PATH_IMAGE152
Figure DEST_PATH_IMAGE153
2)ΔE>δ,调整抽水储能和放水发电功率:
Figure DEST_PATH_IMAGE154
Step 4.1.4:根据约束条件中库容限制约束计算和检查E,返回至Step 4.1.2;
Step 4.1.5:输出净负荷值。
Figure DEST_PATH_IMAGE155
其中,
Figure DEST_PATH_IMAGE156
为除去抽水蓄能电站出力后多能互补系统的剩余负荷,
Figure DEST_PATH_IMAGE157
为抽水蓄能电 站出力。
Step 4.2:针对水电站的水库库容约束、出力约束进行处理。
常规水电站约束处理方式与抽水蓄能电站相同,常规水电调节后的剩余负荷为:
Figure DEST_PATH_IMAGE158
其中,
Figure DEST_PATH_IMAGE159
为除去常规水电站出力后多能互补系统的剩余负荷,
Figure DEST_PATH_IMAGE160
为常规水电站 出力。
Step 4.3:针对火电厂的爬坡速率约束、最小启停时间约束、出力约束进行处理。
Step 4.3.1:根据约束条件中爬坡速率约束和最小启停时间约束确定火电机组状 态
Figure DEST_PATH_IMAGE161
Step 4.3.2:根据约束条件中出力限制约束和爬坡速率约束,与抽水蓄能约束处理方式相同,通过回溯处理迭代求解出负荷分配情况。
Step 4.4:通过系统功率平衡约束将各个子模型进行耦合。
即针对抽水蓄能、常规水电站和火电厂的复杂约束进行处理,并通过多能互补系统功率平衡约束将四个子系统进行耦合。
Step 5:将各参数超出约束值构建为惩罚函数,并添加至所述目标函数。
Step 6:计算所述目标函数的适应度值。
Step 7:根据所述适应度值确定非支配解,更新粒子速度和位置,并继续执行Step4至Step7,直至达到最大迭代次数时输出Pareto解集。
优选方案中,对所述日前优化调度模型进行求解时,可将所述日前优化调度模型划分为四个子模型进行求解,所述四个子模型分别为风电与光伏子模型、抽水蓄能子模型、常规水电站子模型、火电厂子模型。由于本发明所建模型是一个大规模、多维数、多时段、多约束非连续问题,因此根据各能源调度顺序将多能互补模型划分四个子模型能够降低多目标优化问题庞大的维数,简化求解的复杂性。
步骤3:基于隶属度函数,从所述Pareto解集中选取一个最优折中解,并作为日前优化调度方案。
由于多目标优化时,各个目标不能够同时达到最优解,所以使用MOHPSO算法得出的结果是一个Pareto解集,只能从Pareto最优解集中折中选取一个综合结果最好的值作为最优解。一般根据实际情况可以有侧重的选取,本发明参考模糊控制中的隶属度函数的概念,采用归一化隶属度的方法从Pareto最优解集中选取一个最优折中解,使得更有说服力。
Step 3.1:隶属度函数的取值范围为0~1,目标函数值越小则对应的隶属度越大,越符合选择要求,也符合MOHPSO算法的选择规律。本发明采用线性隶属度函数:
Figure DEST_PATH_IMAGE162
其中,
Figure DEST_PATH_IMAGE163
为第j个目标函数值,
Figure DEST_PATH_IMAGE164
,这里由于目标函数最值无法提前确 定,故取计算出的最终的Pareto最优解集中各目标函数的最值作为对应范围。
Step 3.2:计算出每个非支配解的隶属度
Figure DEST_PATH_IMAGE165
以及归一化隶属度值
Figure DEST_PATH_IMAGE166
,计算式如 下:
Figure DEST_PATH_IMAGE167
Figure DEST_PATH_IMAGE168
其中,m为非支配解的个数,z为目标函数的个数。
Step 3.3:计算后选取归一化隶属度值最大的解即为最优解。
综上,本发明提出一种基于抽水蓄能优先调节的多能互补系统日前优化调度方法,可以优化电源结构、缓解可再生能源对电网的冲击、提高多能互补系统的经济性、降低碳排放量、释放了常规水电和火电机组的调节压力、减少了因机组频繁调节而造成的损耗,为多能互补系统调度运行提供支撑。
为验证本发明的效果,下面以某能源基地(风电装机容量500MW,光伏装机容量300MW,水电装机容量300MW,火电总装机容量1662MW,抽水蓄能装机容量150MW)为多能互补模拟系统,定量分析抽水蓄能是否参与调节、抽蓄调节次序对多能互补系统运行的影响,下面设置3个实验组:无抽蓄参加-S1、抽蓄优先调节-S2(即本发明)和常规水电优先调节-S3的运行模式,并进行详细说明:
步骤(1):设定风电、光伏全额入网。
风电出力值、光伏出力值、多能互补系统剩余负荷如图2所示,常规水电站24小时入流量如表1所示。
表1常规水电站24h入流量
Figure DEST_PATH_IMAGE169
步骤(2):设置MOHPSO算法、抽水蓄能电站参数。
MOHPSO算法和抽水蓄能电站的参数分别如表2、表3所示。
表2 MOHPSO参数
Figure DEST_PATH_IMAGE170
表3抽水蓄能电站参数
Figure DEST_PATH_IMAGE171
步骤(3):基于MOHPSO求解并输出Pareto最优解
基于本发明所建模型和求解方法,运用MOHPSO算法求解得出3个实验组的Pareto解集如图3所示,由该图可以直观看出抽水蓄能参与调度时可降低多能互补系统运行成本和碳排放量,抽蓄优先调节的运行模式优于常规水电优先调节的运行模式。
步骤(4):最优解选取及结果分析
(1)最优解选取
以图3中S1的Pareto最优解选取方法求出20组解归一化隶属度值如图4所示,由于归一化隶属度值越大约优,故选择第3组为S1的相对最优解。根据图3确定出3种方案各自的相对最优解,各方案详细情况为:多能互补系统出力计划和运行细节分别如图5和表4所示,三个实验组中常规水电站下泄流量情况图6如所示,三个实验组中常规水电站水库库容情况如图7所示,三个实验组中目标函数值和各子系统评价指标值如图8所示,其中图8中的(a)为系统运行成本,图8中的(b)为碳排放量,图8中的(c)为常规水电调节里程,图8中的(d)为火电机组运行时间。
表4多能互补系统运行细节参数
Figure DEST_PATH_IMAGE172
(2)结果分析。
1)抽水蓄能是否参与调节。
本节针对S1(无抽蓄)和S2(抽蓄优先调节)2种运行模式进行分析。由图5分析可知,不同运行模式时火电机组的出力差别较大,为应对可再生能源的波动性和调峰需求,S1模式先后启动了全部火电机组,并且在21点负荷最高峰时10台机组均在出力;S2模式先后启动了1-7台机组,同样在21点7台机组全部出力,用于调峰需要,这一情况也可以从表4看出。表4包括3个实验组各时段多能互补系统中火电机组启停情况(1/0代表开/关)、启停成本和燃料成本,分析可知,火电机组1-3由于其优越的经济性和较低的碳排放量,在S1、S2中均保持开启状态以提供基荷。机组8-10因其效率低、污染大,非必要时段均保持关闭状态,但S1模式由于调峰需要在19:00-22:00时段开启了这3台低性能机组,而S2模式8-10号机组始终保持关闭状态。
由图5分析可知,抽水蓄能的抽水与发电工况与剩余负荷(除去风光)的波动情况高度相关,其主要在1:00-5:00的负荷低谷时段抽水储能,在18:00-22:00的负荷高峰时段放水发电,经抽水蓄能调节后的净负荷方差为10069、最大峰谷差为345MW,较多能互补系统未接入抽蓄时,分别下降了63.52%、45.15%,这充分体现出抽水蓄能的调峰填谷作用。由图8分析可知,抽水蓄能的加入导致多能互补系统运行成本降低11698$(启停成本和燃料成本分别降低470$和11228$)、碳排放量减少2524kg、常规水电调节里程减少25276MW2、火电机组运行时间减少12h。
2)抽蓄优先调节模式的优势:不同调节次序的对比。
由于常规水电和抽水蓄能均为可调可再生能源,两者调节顺序对多能互补系统的影响以及抽水蓄能的加入对常规水电站的影响有待研究,此节主要针对抽水蓄能不同调节顺序进行对比分析,即针对S2和S3模式进行分析。由图5分析可知,S2与S3模式下抽水蓄能出力情况相似,均为负荷低谷抽水、负荷高峰发电,启停次数均为3次(2抽1发)。结合图6、图7和表1可以看出,S3模式常规水电站下泄流量、水库库容变化与负荷变化和出力情况保持一致趋势,在1:00-6:00的负荷低谷时段由于剩余负荷(除去风光)呈逐渐增加趋势,常规水电站下泄流量也随之变化导致水库库容呈减小趋势;在18:00-22:00的负荷高峰时段常规水电站为响应调峰需求,下泄流量大于入流量,使得水库库容也呈现减小趋势。正是由于常规水电站调节(尤其是调峰)的需要,导致其调节强度较大,后期通过抽水蓄能一定的调节使剩余负荷(除去风光、常规水电、抽水蓄能出力)方差降至9271,最大峰谷差降至340MW。但S2模式由于抽水蓄能的调峰填谷作用使净负荷更平稳,使常规水电站不需要频繁改变出力以应对负荷波动,经抽蓄和常规水电调节后的剩余负荷(除去风光、常规水电、抽水蓄能出力)方差为8831,最大峰谷差为336MW。此时,常规水电站除调峰时段仍需响应外,其它时段出力较为平稳。常规水电站的出力波动情况和调节强度也可从图8的调节里程看出,S2模式调节里程为63643MW2,S3模式则为97689MW2
通过对上述3种运行模式进行对比分析,可知本发明具有如下优点和效果:
(1)较常规水电优先调节模式,本发明所提出的调度模式通过优化电源结构使得多能互补系统经济性提高0.98%,碳排放量降低0.64%,子系统中的常规水电站调节里程和火电机组运行小时数分别降低34.85%和4.72%,而抽水蓄能均为2抽1发;
(2)本发明提出的抽水蓄能优先调节的多能互补运行模式可以优化电源结构、缓解可再生能源对电网的冲击、提高多能互补系统的经济性、降低碳排放量、释放了常规水电和火电机组的调节压力、减少了因机组频繁调节而造成的损耗;
(3)针对抽水蓄能与常规水电站无直接水力联系的多能互补系统优化调度求解问题,本发明通过将多能互补系统调度模型分解为4个子系统,并分部计算的求解方法可以降低多目标优化问题庞大的维数和规模,简化求解的复杂性。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (9)

1.基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,包括以下步骤:
建立风-光-水-火-抽水蓄能多能互补系统的日前优化调度模型;所述日前优化调度模型包括目标函数和约束条件;将多能互补系统的运行成本最小和所述多能互补系统的碳排放量最低作为所述目标函数;所述多能互补系统中各能源调度的顺序依次为:风电与光伏、抽水蓄能、常规水电站、火电厂;
获取各能源的基础数据;
基于所述各能源的基础数据,采用多目标混合粒子群算法对所述日前优化调度模型进行求解,求解的过程中基于所述约束条件进行约束处理,得到Pareto解集;
基于隶属度函数,从所述Pareto解集中选取一个最优折中解,并作为日前优化调度方案;
所述目标函数包括第一目标函数和第二目标函数;
所述第一目标函数如下:
Figure QLYQS_1
Figure QLYQS_2
Figure QLYQS_3
式中,
Figure QLYQS_5
为第一目标函数,对应多能互补系统的运行成本;N为火电厂的火电机组的数量;T为时段数目;/>
Figure QLYQS_7
为第j台火电机组在i时段的启停状态,/>
Figure QLYQS_9
表示机组开启,/>
Figure QLYQS_12
表示机组关闭;/>
Figure QLYQS_14
为第j台火电机组在i-1时段的启停状态;/>
Figure QLYQS_15
为第j台火电机组在i时段的启动成本;/>
Figure QLYQS_17
为第j台火电机组在i时段的出力;/>
Figure QLYQS_4
、/>
Figure QLYQS_6
、/>
Figure QLYQS_8
分别为第j台火电机组的运行成本系数;/>
Figure QLYQS_10
表示第j台火电机组的火电热启动成本;/>
Figure QLYQS_11
表示第j台火电机组的火电冷启动成本;/>
Figure QLYQS_13
为第j台火电机组的最小停机时长,/>
Figure QLYQS_16
为第j台火电机组在i时段的累计停机时长,/>
Figure QLYQS_18
为第j台火电机组的冷启动时长;
所述第二目标函数如下:
Figure QLYQS_19
式中,
Figure QLYQS_20
为第二目标函数,对应多能互补系统的碳排放量;/>
Figure QLYQS_21
、/>
Figure QLYQS_22
、/>
Figure QLYQS_23
分别为第j台火电机组的碳排放系数。
2.根据权利要求1所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,所述日前优化调度模型还包括评价指标,所述评价指标用于衡量所述多能互补系统的运行状况和调节强度。
3.根据权利要求2所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,所述评价指标包括:火电厂运行小时数、常规水电站调节里程和抽水蓄能电站启停次数。
4.根据权利要求3所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,所述常规水电站调节里程表示为:
Figure QLYQS_24
式中,
Figure QLYQS_25
为常规水电站调节里程,/>
Figure QLYQS_26
为水电机组的额定功率,T为时段数目,/>
Figure QLYQS_27
为水电机组在i时刻的出力,/>
Figure QLYQS_28
为水电机组在i+1时刻的出力。
5.根据权利要求1所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,所述各能源的基础数据包括风电出力、光伏出力、多能互补系统的负荷、常规水电站的水库入流量。
6.根据权利要求1所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,所述约束条件包括:所述多能互补系统的系统约束条件、新能源约束条件、水电站约束条件、火电厂约束条件和抽水蓄能电站约束条件;
所述多能互补系统的系统约束条件包括系统功率平衡约束和系统旋转备用约束;所述新能源约束条件包括风电出力约束和光伏出力约束;所述水电站约束条件包括水电站的出力约束、水库库容约束和水库流量约束;所述火电厂约束条件包括火电机组出力约束、爬坡速率约束和最小启停时间约束;所述抽水蓄能电站约束条件包括抽水蓄能电站的发电功率约束、抽水功率约束、工况约束和库容约束。
7.根据权利要求1所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,对所述日前优化调度模型进行求解时,将所述日前优化调度模型划分为四个子模型进行求解,所述四个子模型分别为风电与光伏子模型、抽水蓄能子模型、常规水电站子模型、火电厂子模型。
8.根据权利要求7所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,采用多目标混合粒子群算法对所述日前优化调度模型进行求解包括如下步骤:
步骤1、设定风电、光伏全额入网,计算所述多能互补系统除去新能源后的剩余负荷;
步骤2、建立所述多能互补系统的决策变量矩阵并设置参数;所述决策变量矩阵包括抽水蓄能的抽水功率变量和发电功率变量、常规水电站的流量变量、火电厂的机组运行状态变量;所述参数包括种群大小、迭代次数、学习因子、惯性权重;
步骤3、种群初始化;
步骤4、进行约束处理;
步骤5、将各参数超出约束值构建为惩罚函数,并添加至所述目标函数;
步骤6、计算所述目标函数的适应度值;
步骤7、根据所述适应度值确定非支配解,更新粒子速度和位置,并继续执行步骤 4至步骤7,直至达到最大迭代次数时输出Pareto解集。
9.根据权利要求8所述的基于抽水蓄能优先调节的多能互补系统日前优化调度方法,其特征在于,所述步骤4中进行约束处理包括如下子步骤:
步骤4.1、针对抽水蓄能的库容约束、抽水功率约束、发电功率约束进行处理;
步骤4.2、针对水电站的水库库容约束、出力约束进行处理;
步骤4.3、针对火电厂的爬坡速率约束、最小启停时间约束、出力约束进行处理;
步骤4.4、通过系统功率平衡约束将各个子模型进行耦合。
CN202211442469.1A 2022-11-18 2022-11-18 基于抽水蓄能优先调节的多能互补系统日前优化调度方法 Active CN115640982B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211442469.1A CN115640982B (zh) 2022-11-18 2022-11-18 基于抽水蓄能优先调节的多能互补系统日前优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211442469.1A CN115640982B (zh) 2022-11-18 2022-11-18 基于抽水蓄能优先调节的多能互补系统日前优化调度方法

Publications (2)

Publication Number Publication Date
CN115640982A CN115640982A (zh) 2023-01-24
CN115640982B true CN115640982B (zh) 2023-03-31

Family

ID=84948674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211442469.1A Active CN115640982B (zh) 2022-11-18 2022-11-18 基于抽水蓄能优先调节的多能互补系统日前优化调度方法

Country Status (1)

Country Link
CN (1) CN115640982B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116231767B (zh) * 2023-05-11 2023-07-14 南方电网调峰调频发电有限公司储能科研院 一种梯级水电站的多能互补调度方法和系统
CN116667395B (zh) * 2023-05-31 2023-11-21 河海大学 基于梯级水电改造的水风光蓄互补泵站容量配置方法
CN116780657B (zh) * 2023-08-17 2024-01-19 长江三峡集团实业发展(北京)有限公司 一种水风蓄互补发电系统的调度方法、装置、设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109767105A (zh) * 2018-12-29 2019-05-17 东北电力大学 一种基于风光水火储联合系统的多能互补协调发电调度方法
CN109830168A (zh) * 2019-02-27 2019-05-31 南方电网调峰调频发电有限公司 一种可模拟海浪水位波动的抽水蓄能模型试验下游水箱
CN112803494A (zh) * 2021-02-19 2021-05-14 国网湖南省电力有限公司 一种含风光水火的多目标agc协调优化方法及系统
CN114243692A (zh) * 2021-12-15 2022-03-25 深圳供电局有限公司 一种源-网-荷协调优化调度方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109767105A (zh) * 2018-12-29 2019-05-17 东北电力大学 一种基于风光水火储联合系统的多能互补协调发电调度方法
CN109830168A (zh) * 2019-02-27 2019-05-31 南方电网调峰调频发电有限公司 一种可模拟海浪水位波动的抽水蓄能模型试验下游水箱
CN112803494A (zh) * 2021-02-19 2021-05-14 国网湖南省电力有限公司 一种含风光水火的多目标agc协调优化方法及系统
CN114243692A (zh) * 2021-12-15 2022-03-25 深圳供电局有限公司 一种源-网-荷协调优化调度方法

Also Published As

Publication number Publication date
CN115640982A (zh) 2023-01-24

Similar Documents

Publication Publication Date Title
CN115640982B (zh) 基于抽水蓄能优先调节的多能互补系统日前优化调度方法
CN103151803B (zh) 一种含风电系统机组及备用配置的优化方法
CN113095791B (zh) 一种综合能源系统运行方法及系统
CN107732949B (zh) 一种综合风电全年多季节特性的储能布点定容方法
CN114221338B (zh) 考虑电源灵活性和互补性的多能源电力系统优化调度方法
CN111210079B (zh) 一种分布式能源虚拟电厂运行优化方法及系统
CN110245794B (zh) 考虑灵活性的多能源汇集中心火储容量双层优化方法
CN112270433B (zh) 考虑可再生能源不确定性和用户满意度的微电网优化方法
CN111668878A (zh) 一种可再生微能源网的优化配置方法和系统
CN115102159A (zh) 计及碳排放流的共享储能优化运行方法
CN114362153B (zh) 一种并网型风光储系统多目标容量优化配置方法及系统
CN111049179A (zh) 一种考虑不确定性的新能源发电系统多目标优化调度方法
CN112307603B (zh) 考虑大规模风电接入的混合储能容量优化配置方法及系统
CN114884101B (zh) 一种基于自适应模型控制预测的抽水蓄能调度方法
CN115940284B (zh) 一种考虑分时电价的新能源制氢系统的运行控制策略
CN109840621A (zh) 考虑储能系统影响的并网型微电网日前多目标优化方法
CN115513997A (zh) 一种中长时间尺度下新能源电源实时备用率的确定方法
CN115310714A (zh) 一种园区能源优化调控方法、调控系统及存储介质
CN114725961A (zh) 平抑风电波动的制氢系统容量分层优化配置方法
CN115021336A (zh) 一种面向高比例新能源电力系统惯量不足场景的优化方法
CN114498690A (zh) 支持大规模可再生能源消纳的多元复合储能优化配置方法
CN112054553A (zh) 电热气互联系统协调优化运行方法、系统、介质及设备
CN112332456A (zh) 一种蓄热式电采暖负荷参与电网运行的优化控制方法及装置
CN116307935B (zh) 一种电力系统运行模拟的抽水蓄能价值评价方法及系统
Chen et al. Coordinated Source-Grid-Storage Planning of High Renewable Energy Proportion Power System Based on Refined Operation Simulation

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