CN112700094B - 一种基于mpc及loddlc的综合能源系统多时间尺度优化调度方法 - Google Patents
一种基于mpc及loddlc的综合能源系统多时间尺度优化调度方法 Download PDFInfo
- Publication number
- CN112700094B CN112700094B CN202011528062.1A CN202011528062A CN112700094B CN 112700094 B CN112700094 B CN 112700094B CN 202011528062 A CN202011528062 A CN 202011528062A CN 112700094 B CN112700094 B CN 112700094B
- Authority
- CN
- China
- Prior art keywords
- time
- gas
- load
- scheduling
- period
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/38—Arrangements for feeding a single network from two or more generators or sources in parallel; Arrangements for feeding already energised networks from additional generators or sources in parallel
- H02J3/46—Controlling the sharing of generated power between the generators, sources or networks
- H02J3/466—Scheduling or selectively controlling the operation of the generators or sources, e.g. connecting or disconnecting generators to meet a demand
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/38—Arrangements for feeding a single network from two or more generators or sources in parallel; Arrangements for feeding already energised networks from additional generators or sources in parallel
- H02J3/46—Controlling the sharing of generated power between the generators, sources or networks
- H02J3/48—Controlling the sharing of active power
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2101/00—Supply or distribution of decentralised, dispersed or local electric power generation
- H02J2101/20—Dispersed power generation using renewable energy sources
- H02J2101/22—Solar energy
- H02J2101/24—Photovoltaics
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2101/00—Supply or distribution of decentralised, dispersed or local electric power generation
- H02J2101/20—Dispersed power generation using renewable energy sources
- H02J2101/28—Wind energy
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2101/00—Supply or distribution of decentralised, dispersed or local electric power generation
- H02J2101/40—Hybrid power plants, i.e. a plurality of different generation technologies being operated at one power plant
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2103/00—Details of circuit arrangements for mains or AC distribution networks
- H02J2103/30—Simulating, planning, modelling, reliability check or computer assisted design [CAD] of electric power networks
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Power Engineering (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Biodiversity & Conservation Biology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法,包括,遵循日前调度结果中的DLC控制时段及容量需求,采用等效热参数模型对DLC控制时段的反弹负荷进行细节层次仿真,并优化选择空调机组的控制序列;利用多级Douglas‑Peucker法压缩冗余负荷数据信息,修正反弹负荷曲线,建立计及细节层次DLC空调仿真日内调度模型;缩小时间尺度,进一步细化日内反弹负荷曲线,以日内优化调度结果为参考曲线,建立基于MPC的实时修正调度模型;利用LINGO软件结合MATLAB编程,对日前‑日内‑实时调度模型进行求解,得到多时间尺度下综合能源系统优化调度策略。本发明降低模型描述反弹负荷误差对系统调度计划的影响,贴合实际系统运行状态,体现多时间尺度调度的优越性。
Description
技术领域
本发明涉及综合能源优化调度的技术领域,尤其涉及一种基于MPC(ModelPredictive Control,模型预测控制)及LODDLC(Level of detail direct load control,细节层次直接负荷控制)的综合能源系统多时间尺度优化调度方法。
背景技术
已有的综合能源系统优化调度方法仅考虑了风、光等可再生能源出力的随机性对优化调度结果的影响,目前考虑直接负荷控制参与调度的研究多采用单一的简化直接负荷控制模型,针对直接负荷控制多时间尺度精细仿真的研究较少。传统阶段性反弹等简化模型难以准确表征其受控后的反弹负荷特性,不同时间尺度下采用不同细节层次的直接负荷控制仿真模型参与调度,有利于逐级修正反弹负荷曲线,降低简化直接负荷控制模型引起的反弹负荷误差对系统调度计划的影响,提高调度计划的可执行度。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于上述现有综合能源系统调度存在的问题,提出了本发明。
因此,本发明提供了一种基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法,能够解决多时间尺度协调运行的问题。
为解决上述技术问题,本发明提供如下技术方案:包括,遵循日前调度结果中的DLC控制时段及容量需求,采用等效热参数模型对所述DLC控制时段的反弹负荷进行细节层次仿真,并优化选择空调机组的控制序列;利用多级Douglas-Peucker法压缩冗余负荷数据信息,修正反弹负荷曲线,建立计及细节层次DLC空调仿真日内调度模型;缩小时间尺度,进一步细化日内反弹负荷曲线,以日内优化调度结果为参考曲线,建立基于MPC的实时修正调度模型;利用LINGO软件结合MATLAB编程,对日前-日内-实时调度模型进行求解,得到多时间尺度下综合能源系统优化调度策略。
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:进行所述细节层次仿真之前还包括,获取综合能源系统的网架结构参数、能源类型及单位成本参数、用户各类型的功率需求和可再生能源预测出力信息;基于日前预测的风电、光伏、负荷、电价信息数据,采用三阶段反弹负荷模型简化模拟空调参与直接负荷控制反弹特性;根据所述综合能源系统的网架结构,建立考虑简化直接负荷控制的综合能源系统的电、热、气、冷功率平衡方程,同时建立计及三阶段反弹简化DLC日前调度模型。
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:所述网架结构参数包括,电网、可再生发电机组、微燃机、天然气网、燃气锅炉、电锅炉、电转气设备P2G、电制冷机、吸收式制冷机、余热锅炉和储能设备;所述可再生发电机组包括,风电机组和光伏电池;所述的储能设备包括,储电设备、储气设备、储热设备和储冷设备。
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:采用所述三阶段反弹负荷模型来初步模拟空调负荷反弹特性,包括,
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:根据所述综合能源系统的架构信息,列写综合能源系统的电、热、气、冷功率平衡方程,包括,
分别为电锅炉、电转气装置、电制冷机耗电功率;为经余热回收后的微燃机制热功率,分别为燃气锅炉和电锅炉产热功率,为吸收式制冷机输入热功率;为电转气设备产气量, 分别为微燃机、燃气锅炉消耗的天然气量,为t时段储气容量; 分别为电制冷机、吸收式制冷机制冷功率,分别为t时段电负荷需求、热负荷需求、气负荷需求、冷负荷需求。
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:还包括,基于综合能源系统中各微源的功能数学模型,在满足系统各项基本约束条件下,考虑联络线交互成本,购气成本,各微源运行维护成本和直接负荷控制补偿成本为系统经济运行目标,计及三阶段反弹简化DLC日前调度模型包含的目标函数包括,
minCd-a=min(CGRID+CGAS+CMC+CDLC)
其中,T为调度时段总数,CGRID、CGAS、CMC、CDLC分别为RIES与主网交互成本,购买天然气成本,微源运行维护成本,直接负荷控制补偿成本,Psgrid(t)和Pbgrid(t)为t时段向主网售电和购电功率,Cbgrid(t)和Csgrid(t)分别为t时段向主网购电和售电电价;为t时段产气量,Cgas为单位天然气购买价格; CEboil、CGboil、Cp2g、CEC、CAC分别为微燃机、光伏、风机、电锅炉、燃气锅炉、电转气装置、电制冷机、吸收式制冷机单位运行维护成本; 为对应设备工作的输入功率;x代表能源类型,以e,h,c,g表示电热冷气四种能源,为各类储能设备单位运行维护成本,分别表示t时段各类储能设备充放功率;CCP(t)为t时段受控负荷中断的单位补偿价格。
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:包括,对所述DLC控制时段的反弹负荷进行分钟级仿真,并优化选择空调机组的所述控制序列,采用所述等效热参数模型表示空调热交换过程,定义空调负荷的初始状态均匀分布,忽略空调的瞬间启动,电网电压波动和空调压缩机堵转的影响,遵循日前调度计划直接负荷控制的时段及削减容量,对空调群运行特性进行模拟仿真;空调在较小温度区间内运行时,其室内温度变化轨迹可以近似用直线模拟,进而得到线性化的等效热参数模型,即状态队列模型,如下,
其中,Ti t、Ti t+1表示第i台空调在t、t+1时刻的室内温度;室内温度设定值区间为[Tmin,Tmax],记空调运行时间为τon,停机时间为τoff,k表示空调的运行状态,1表示空调启动,0表示空调停机;Δt为仿真时间步长。
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:包括,建立有限时域微燃机、联络线、储能系统的有功出力预测模型,包括,
其中,W0(t)为微燃机、联络线、储能系统有功初始值,由实际测量值校正得到,ΔUs(t+j|t)为t+j时刻控制有功增量,Ws(t+i|t)为t时刻预测t+i时刻有功出力值,M表示预测时域长度。
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:还包括,实时修正所述调度模型以所述日内优化调度结果为参考曲线,基于模型预测控制理论,优化周期为15min,分辨率为5min,进一步细化日内反弹负荷曲线,校正各微源出力,为保证各微源尽量跟踪日内调度计划,保证整体优化结果的经济性,其目标函数为所有场景下优化周期内各微源出力相对日内调度计划参考值差值及实时调整增量的期望之和最小,如下,
作为本发明所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的一种优选方案,其中:超前的MPC控制无法保证可再生能源出力与预测值相同,从而导致下发的可控微源出力与实际有功出力之间存在一定偏差,增加反馈校正环节约束,实现能量的闭环控制,包括,
W0(t+1)=Wtrue(t+1)+σ(t+1)
其中,Wtrue(t+1)为t+1时刻各微源的实际运行功率;σ(t+1)为t+1时刻运行采样误差。
本发明的有益效果:本发明方法基于模型预测控制的计及细节层次直接负荷控制的多时间尺度协调优化调度策略,日前-日内-实时三时间尺度调度相互协调配合,随时间尺度的细化,逐步减小预测误差,不同时间尺度下采用不同细节层次的直接负荷控制仿真模型参与调度,逐级修正反弹负荷曲线,与传统采用单一简化DLC模型相比,降低模型描述反弹负荷误差对系统调度计划的影响,进一步贴合实际系统运行状态,体现多时间尺度调度方式的优越性。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为本发明第一个实施例所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的流程示意图;
图2为本发明第一个实施例所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的综合能源系统结构示意图;
图3为本发明第一个实施例所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的Douglas-Peucker数据抽稀示意图;
图4为本发明第二个实施例所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的风机光伏出力及各类负荷日前预测示意图;
图5为本发明第二个实施例所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的细节层次DLC与三阶段反弹DLC反弹负荷对比曲线示意图。
图6为本发明第二个实施例所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的日内调度联络线计划和实时调度基于模型预测控制滚动优化结果对比示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
本发明结合示意图进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的剖面图会不依一般比例作局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。此外,在实际制作中应包含长度、宽度及深度的三维空间尺寸。
同时在本发明的描述中,需要说明的是,术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一、第二或第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
本发明中除非另有明确的规定和限定,术语“安装、相连、连接”应做广义理解,例如:可以是固定连接、可拆卸连接或一体式连接;同样可以是机械连接、电连接或直接连接,也可以通过中间媒介间接相连,也可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
参照图1、图2和图3,为本发明的第一个实施例,提供了一种基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法,包括:
S1获取综合能源系统的网架结构参数、能源类型及单位成本参数、用户各类型的功率需求和可再生能源预测出力信息。参照图2,其中需要说明的是:
网架结构参数包括,电网、可再生发电机组、微燃机、天然气网、燃气锅炉、电锅炉、电转气设备P2G、电制冷机、吸收式制冷机、余热锅炉和储能设备;
可再生发电机组包括,风电机组和光伏电池;
的储能设备包括,储电设备、储气设备、储热设备和储冷设备。
S2基于日前预测的风电、光伏、负荷、电价信息数据,采用三阶段反弹负荷模型简化模拟空调参与直接负荷控制反弹特性。本步骤需要说明的是:
S3根据综合能源系统的网架结构,建立考虑简化直接负荷控制的综合能源系统的电、热、气、冷功率平衡方程,同时建立计及三阶段反弹简化DLC日前调度模型。其中需要说明的是:
根据综合能源系统的架构信息,列写综合能源系统的电、热、气、冷功率平衡方程,包括,
其中,和分别为t时段光伏、风机和微燃机发电功率,分别为电锅炉、电转气装置、电制冷机耗电功率;为经余热回收后的微燃机制热功率,分别为燃气锅炉和电锅炉产热功率,为吸收式制冷机输入热功率;为电转气设备产气量, 分别为微燃机、燃气锅炉消耗的天然气量,为t时段储气容量; 分别为电制冷机、吸收式制冷机制冷功率,分别为t时段电负荷需求、热负荷需求、气负荷需求、冷负荷需求;
基于综合能源系统中各微源的功能数学模型,在满足系统各项基本约束条件下,考虑联络线交互成本,购气成本,各微源运行维护成本和直接负荷控制补偿成本为系统经济运行目标,计及三阶段反弹简化DLC日前调度模型包含的目标函数包括,
minCd-a=min(CGRID+CGAS+CMC+CDLC)
其中,T为调度时段总数,CGRID、CGAS、CMC、CDLC分别为RIES与主网交互成本,购买天然气成本,微源运行维护成本,直接负荷控制补偿成本,Psgrid(t)和Pbgrid(t)为t时段向主网售电和购电功率,Cbgrid(t)和Csgrid(t)分别为t时段向主网购电和售电电价;为t时段产气量,Cgas为单位天然气购买价格; CEboil、CGboil、Cp2g、CEC、CAC分别为微燃机、光伏、风机、电锅炉、燃气锅炉、电转气装置、电制冷机、吸收式制冷机单位运行维护成本; 为对应设备工作的输入功率;x代表能源类型,以e,h,c,g表示电热冷气四种能源,为各类储能设备单位运行维护成本,分别表示t时段各类储能设备充放功率;CCP(t)为t时段受控负荷中断的单位补偿价格。
S4遵循日前调度结果中的DLC控制时段及容量需求,采用等效热参数模型对DLC控制时段的反弹负荷进行细节层次仿真,并优化选择空调机组的控制序列。本步骤需要说明的是:
对DLC控制时段的反弹负荷进行分钟级仿真,并优化选择空调机组的控制序列,采用等效热参数模型表示空调热交换过程,定义空调负荷的初始状态均匀分布,忽略空调的瞬间启动,电网电压波动和空调压缩机堵转的影响,遵循日前调度计划直接负荷控制的时段及削减容量,对空调群运行特性进行模拟仿真;
空调在较小温度区间内运行时,其室内温度变化轨迹可以近似用直线模拟,进而得到线性化的等效热参数模型,即状态队列模型,如下,
其中,Ti t、Ti t+1表示第i台空调在t、t+1时刻的室内温度;室内温度设定值区间为[Tmin,Tmax],记空调运行时间为τon,停机时间为τoff,k表示空调的运行状态,1表示空调启动,0表示空调停机;Δt为仿真时间步长。
S5利用多级Douglas-Peucker法压缩冗余负荷数据信息,修正反弹负荷曲线,建立计及细节层次DLC空调仿真日内调度模型。参照图3,本步骤还需要说明的是:
基于多级Douglas-Peucker算法数据压缩冗余负荷数据信息,修正反弹负荷曲线:
常规Douglas-Peucker法的基本思路如图3所示,用直线连接曲线的首末两端,找出曲线中所有点到该直线的最大距离值dmax,比较dmax与精度限差D的大小:
(1)如果dmax<D,则保留曲线首末端点;
(2)如果dmax≥D,则保留dmax对应的点,并以该点为界限,将原曲线分为左右两部分,对左右两部分分别重复上述过程;
改进的多级Douglas-Peucker数据抽稀法分为整体抽稀与局部抽稀,将时序仿真结果适度压缩,以便进一步修正负荷曲线返回至RIES日内滚动调度模型中,修正调度结果;
第一级数据抽稀:选取精度限差D对目标曲线整体抽稀,以保证整体曲线不失真;
第二级数据抽稀:选取第j时段精度限差Dj对各时段曲线局部抽稀,以达到局部平整化效果。
计及细节层次DLC空调仿真日内调度每4小时预测日内风、光及各类负荷数据,优化周期为4小时,分辨率为15min,其模型包含的目标函数参见步骤2中日前调度模型,除了满足日前调度中各微源及储能系统的约束外,日内调度各状态变量应满足日前调度结果。
S6缩小时间尺度,进一步细化日内反弹负荷曲线,以日内优化调度结果为参考曲线,建立基于MPC的实时修正调度模型。其中还需要说明的是:
建立有限时域微燃机、联络线、储能系统的有功出力预测模型,包括,
其中,W0(t)为微燃机、联络线、储能系统有功初始值,由实际测量值校正得到,ΔUs(t+j|t)为t+j时刻控制有功增量,Ws(t+i|t)为t时刻预测t+i时刻有功出力值,M表示预测时域长度;
实时修正调度模型以日内优化调度结果为参考曲线,基于模型预测控制理论,优化周期为15min,分辨率为5min,进一步细化日内反弹负荷曲线,校正各微源出力,为保证各微源尽量跟踪日内调度计划,保证整体优化结果的经济性,其目标函数为所有场景下优化周期内各微源出力相对日内调度计划参考值差值及实时调整增量的期望之和最小,如下,
超前的MPC控制无法保证可再生能源出力与预测值相同,从而导致下发的可控微源出力与实际有功出力之间存在一定偏差,增加反馈校正环节约束,实现能量的闭环控制,包括,
W0(t+1)=Wtrue(t+1)+σ(t+1)
其中,Wtrue(t+1)为t+1时刻各微源的实际运行功率;σ(t+1)为t+1时刻运行采样误差。
S7利用LINGO软件结合MATLAB编程,对日前-日内-实时调度模型进行求解,得到多时间尺度下综合能源系统优化调度策略。
不难理解的是,直接负荷控制(direct load control,DLC)是一种激励型需求响应,一般针对居民用户或小型商业用户中具有热储存能力的温控负荷(thermostaticallycontrolled loads,TCL)实施,如空调、电热水器,受控负荷在中断时段从系统中断开,在中断时段结束后重新与系统连接,形成反弹负荷。
而本发明涉及考虑风光出力和负荷的不确定性及温控负荷反弹特性的精细建模程度对区域综合能源系统优化调度的影响,日前调度采用三阶段反弹简化直接负荷控制模型,以系统运行经济性最优为目标求解机组运行计划,日内调度基于细节层次直接负荷控制仿真修正日前反弹负荷曲线,求解各微源的出力计划,实时调度进一步细化日内仿真反弹负荷曲线,结合模型预测控制理论,通过反馈校正和滚动优化调整微源出力,实现多时间尺度的协调运行。
实施例2
参照图4和图5,为本发明的第二个实施例,该实施例不同于第一个实施例的是,提供了一种基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法的验证,包括:
参照图2,以区域综合能源系统为例,系统各机组参数如表1所示,可再生能源出力和各类负荷的日前预测数据参照图4,,日内预测数据以及实时数据根据日前预测数据及各自的概率分布生成,设定光伏、风机和各类负荷需求均遵循正态分布,日前三阶段反弹负荷模型的经验反弹系数α,β,γ分别为0.6,0.2,0.1,空调负荷基本参数如表2所示,设置区域内受控空调共100台,分为4组,1、2组空调设定温度区间[25℃,27℃],3、4组空调设定温度区间[23℃,25℃],设置各台空调负荷的初始状态为均匀分布,受控负荷接受控制的时段为14:00-17:00,定义其对应时段外界环境温度为定值32℃。
表1:区域综合能源系统微源运行参数表。
表2:空调负荷的运行参数表。
| 含义 | 参数 | 数值 |
| 室外温度 | T<sub>0</sub> | 32℃ |
| 空调的能效比 | η | 2.7 |
| 单台空调功率 | P | 2.5kW |
| 等效热阻 | R | 0.18kWh/℃ |
| 等效热容 | C | 5.56℃/kW |
| 仿真时间间隔 | △t | 1min |
| 空调台数 | N | 100 |
| 1、2组温度限值 | [Tmin1,Tmax1] | [25℃,27℃] |
| 3、4组温度限值 | [Tmin2,Tmax2] | [23℃,25℃] |
在日前调度中,采用了三阶段经验系数反弹模型来简化模拟温控负荷的反弹特性,求解出受控时段及削减容量,参照图5,59-62时段为净负荷高峰阶段,风电光伏出力减少且电网电价处于峰时段,系统功率缺额严重,DLC参与系统调度降低系统峰荷,减少功率缺额,在60-65时段出现反弹负荷,反弹负荷波动较小。
在日内调度中,遵循日前调度结果下发的DLC削减时段及容量需求,对空调群负荷进行等效热参数建模仿真,从图5中可以看出,计及细节DLC空调仿真的反弹负荷与简化三阶段反弹负荷有明显差异,反弹负荷反复波动剧烈,造成负荷的二次高峰,究其原因是削减负荷时控制部分空调停机,导致空调机组的温度变化进而影响整个空调群原有的状态队列正常运行造成的,单独采用简化模型描述反弹负荷造成的误差降低调度计划的执行度,增加系统应对突然负荷变化的压力,通过计及细节层次DLC仿真,逐步修正反弹负荷数据,更能贴合实际系统运行状态,提高调度计划的准确性。
实时修正调度模型进一步细化日内反弹负荷曲线,以日内优化调度结果为参考曲线,基于模型预测控制理论,校正各微源出力。日内调度联络线计划和实时调度基于模型预测控制滚动优化的结果对比如图6所示:
从图6中可以看出,基于MPC的滚动优化平缓了日内调度计划曲线。在实时调整阶段,随时间尺度缩小,进一步细化了日内反弹负荷曲线。在日内调度计划的指导下,计及源荷的预测误差实时滚动调整,利用反馈校正构成闭环优化,在保证整体经济的同时,减小了联络线交互功率的波动性。通过对比有无MPC滚动优化联络线功率的波动情况,可以看出基于模型预测控制的滚动优化既有效降低了预测误差对实际系统经济运行的影响。
表3:基于不同DLC模型下调度结果表。
表3中对比了采用简化DLC反弹模型和考虑细节层次DLC模型的日内调度结果及不同模型下计及反弹负荷的负荷曲线特征指标,从结果中可以看出,考虑细节层次DLC模型后增加了调度成本,但更贴合实际运行结果,对比分析关键时段的负荷波动率、负荷峰谷差、反弹负荷峰值等指标,实际中空调群的反弹负荷呈现反复波动特性,反弹负荷波动性更强,峰谷差更大,计及细节层次直接负荷控制模型可较为准确地表征反弹负荷特性,有利于系统制定出更符合实际运行的调度计划,增强应对实际负荷的能力,提高调度计划的可行性。
本发明方法基于模型预测控制的计及细节层次直接负荷控制的多时间尺度协调优化调度策略,日前-日内-实时三时间尺度调度相互协调配合,随时间尺度的细化,逐步减小预测误差,不同时间尺度下采用不同细节层次的直接负荷控制仿真模型参与调度,逐级修正反弹负荷曲线,与传统采用单一简化DLC模型相比,降低模型描述反弹负荷误差对系统调度计划的影响,进一步贴合实际系统运行状态,体现多时间尺度调度方式的优越性。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法,其特征在于:包括,
获取综合能源系统的网架结构参数、能源类型及单位成本参数、用户各类型的功率需求和可再生能源预测出力信息;
基于日前预测的风电、光伏、负荷、电价信息数据,采用三阶段反弹负荷模型简化模拟空调参与直接负荷控制反弹特性;
根据所述综合能源系统的网架结构,建立考虑简化直接负荷控制的综合能源系统的电、热、气、冷功率平衡方程,同时建立计及三阶段反弹简化DLC日前调度模型;
遵循日前调度结果中的DLC控制时段及容量需求,采用等效热参数模型对所述DLC控制时段的反弹负荷进行细节层次仿真,并优化选择空调机组的控制序列;
对所述DLC控制时段的反弹负荷进行分钟级仿真,并优化选择空调机组的所述控制序列,采用所述等效热参数模型表示空调热交换过程,定义空调负荷的初始状态均匀分布,忽略空调的瞬间启动,电网电压波动和空调压缩机堵转的影响,遵循日前调度计划直接负荷控制的时段及削减容量,对空调群运行特性进行模拟仿真;
空调在较小温度区间内运行时,其室内温度变化轨迹用直线模拟,进而得到线性化的等效热参数模型,即状态队列模型,如下,
其中,Ti t、Ti t+1表示第i台空调在t、t+1时刻的室内温度;室内温度设定值区间为[Tmin,Tmax],记空调运行时间为τon,停机时间为τoff,k表示空调的运行状态,1表示空调启动,0表示空调停机;Δt为仿真时间步长;
利用多级Douglas-Peucker法压缩冗余负荷数据信息,修正反弹负荷曲线,建立计及细节层次DLC空调仿真日内调度模型;
缩小时间尺度,进一步细化日内反弹负荷曲线,以日内优化调度结果为参考曲线,建立基于MPC的实时修正调度模型;
实时修正所述调度模型以所述日内优化调度结果为参考曲线,基于模型预测控制理论,优化周期为15min,分辨率为5min,进一步细化日内反弹负荷曲线,校正各微源出力,为保证各微源尽量跟踪日内调度计划,保证整体优化结果的经济性,其目标函数为所有场景下优化周期内各微源出力相对日内调度计划参考值差值及实时调整增量的期望之和最小,如下,
利用LINGO软件结合MATLAB编程,对日前-日内-实时调度模型进行求解,得到多时间尺度下综合能源系统优化调度策略。
2.根据权利要求1所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法,其特征在于:所述网架结构参数包括,电网、可再生发电机组、微燃机、天然气网、燃气锅炉、电锅炉、电转气设备P2G、电制冷机、吸收式制冷机、余热锅炉和储能设备;
所述可再生发电机组包括,风电机组和光伏电池;
所述储能设备包括,储电设备、储气设备、储热设备和储冷设备。
5.根据权利要求4所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法,其特征在于:还包括,
基于综合能源系统中各微源的功能数学模型,在满足系统各项基本约束条件下,考虑联络线交互成本,购气成本,各微源运行维护成本和直接负荷控制补偿成本为系统经济运行目标,计及三阶段反弹简化DLC日前调度模型包含的目标函数包括,
min Cd-a=min(CGRID+CGAS+CMC+CDLC)
其中,T为调度时段总数,CGRID、CGAS、CMC、CDLC分别为RIES与主网交互成本,购买天然气成本,微源运行维护成本,直接负荷控制补偿成本,Psgrid(t)和Pbgrid(t)为t时段向主网售电和购电功率,Cbgrid(t)和Csgrid(t)分别为t时段向主网购电和售电电价;为t时段产气量,Cgas为单位天然气购买价格; CEboil、CGboil、Cp2g、CEC、CAC分别为微燃机、光伏、风机、电锅炉、燃气锅炉、电转气装置、电制冷机、吸收式制冷机单位运行维护成本; 为对应设备工作的输入功率;x代表能源类型,以e,h,c,g表示电热冷气四种能源,为各类储能设备单位运行维护成本,分别表示t时段各类储能设备充放功率;CCP(t)为t时段受控负荷中断的单位补偿价格。
7.根据权利要求6所述的基于MPC及LODDLC的综合能源系统多时间尺度优化调度方法,其特征在于:超前的MPC控制无法保证可再生能源出力与预测值相同,从而导致下发的可控微源出力与实际有功出力之间存在一定偏差,增加反馈校正环节约束,实现能量的闭环控制,包括,
W0(t+1)=Wtrue(t+1)+σ(t+1)
其中,Wtrue(t+1)为t+1时刻各微源的实际运行功率;σ(t+1)为t+1时刻运行采样误差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011528062.1A CN112700094B (zh) | 2020-12-22 | 2020-12-22 | 一种基于mpc及loddlc的综合能源系统多时间尺度优化调度方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011528062.1A CN112700094B (zh) | 2020-12-22 | 2020-12-22 | 一种基于mpc及loddlc的综合能源系统多时间尺度优化调度方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112700094A CN112700094A (zh) | 2021-04-23 |
| CN112700094B true CN112700094B (zh) | 2022-06-17 |
Family
ID=75510552
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011528062.1A Active CN112700094B (zh) | 2020-12-22 | 2020-12-22 | 一种基于mpc及loddlc的综合能源系统多时间尺度优化调度方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112700094B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113739296B (zh) * | 2021-09-08 | 2022-09-06 | 山东佐耀科技有限公司 | 基于模型预测控制的空气源热泵负荷水温控制方法及系统 |
| CN114066031B (zh) * | 2021-11-08 | 2025-03-07 | 国网山东综合能源服务有限公司 | 一种综合能源系统的日内优化调度方法及系统 |
| CN114202102B (zh) * | 2021-11-14 | 2024-08-27 | 国网辽宁省电力有限公司锦州供电公司 | 一种多时间尺度空调负荷预测方法 |
| CN114371619B (zh) * | 2021-12-17 | 2024-04-23 | 上海电力大学 | 一种mgt-cchp变工况动态能效优化控制方法 |
| CN114298440A (zh) * | 2022-01-04 | 2022-04-08 | 西安交通大学 | 建筑含氢多能源系统供需协同运行优化方法及控制装置 |
| CN114707711B (zh) * | 2022-03-23 | 2022-09-16 | 特斯联科技集团有限公司 | 园区制冷机组多时间尺度最优调度方法及系统 |
| CN115526430A (zh) * | 2022-10-31 | 2022-12-27 | 广东电网有限责任公司广州供电局 | 多距离聚类与信息聚合的负荷区间预测方法、系统及介质 |
| CN116047885B (zh) * | 2022-11-01 | 2025-10-28 | 大连理工大学 | 一种基于递阶预测控制模型的螺杆式压缩机组运行优化方法 |
| CN116345564B (zh) * | 2023-02-17 | 2025-09-19 | 国网山东省电力公司德州供电公司 | 综合能源系统多时间尺度分布式协同优化调度方法及系统 |
| CN115864542B (zh) * | 2023-02-24 | 2023-05-05 | 南方电网数字电网研究院有限公司 | 电力机组调度模型的优化方法、装置、设备和存储介质 |
| CN116599060B (zh) * | 2023-07-18 | 2023-11-21 | 山东大学 | 一种主动配电网一体化调度方法、系统、终端设备及介质 |
| CN119154408A (zh) * | 2024-11-14 | 2024-12-17 | 国网(苏州)城市能源研究院有限责任公司 | 一种多源微网优化调度方法及系统 |
| CN120447358B (zh) * | 2025-04-07 | 2025-11-25 | 国能神华九江发电有限责任公司 | 基于智能控制算法的火电机组转态调峰负荷巡航控制方法 |
| CN120414724B (zh) * | 2025-04-25 | 2026-01-27 | 四川蜀道城乡投资集团有限责任公司 | 一种基于智能建筑绿色能源系统的电力管理优化方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109004686A (zh) * | 2018-08-29 | 2018-12-14 | 三峡大学 | 一种考虑冰蓄冷空调多模式的冷热电联供型微网系统 |
| CN110148064A (zh) * | 2019-03-08 | 2019-08-20 | 国网浙江省电力有限公司台州供电公司 | 基于激励型综合需求响应的实时新能源辅助消纳方法 |
| CN110649594A (zh) * | 2019-08-28 | 2020-01-03 | 广州供电局有限公司 | 一种基于多能协同的工业园区综合需求响应调度方法 |
| CN110807588A (zh) * | 2019-10-31 | 2020-02-18 | 国网河南省电力公司经济技术研究院 | 一种多能源耦合综合能源系统的优化调度方法 |
| CN111738529A (zh) * | 2020-07-28 | 2020-10-02 | 南方电网科学研究院有限责任公司 | 基于强化学习的综合能源系统需求响应方法、系统及设备 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110197312A (zh) * | 2019-06-28 | 2019-09-03 | 东南大学 | 一种基于多时间尺度的用户级综合能源系统优化调度方法 |
| CN110417006B (zh) * | 2019-07-24 | 2023-01-06 | 三峡大学 | 一种综合能源系统多时间尺度能量调度方法 |
| CN112072712B (zh) * | 2020-08-31 | 2022-07-12 | 合肥工业大学 | 综合能源系统多时间尺度优化调度方法、系统及存储介质 |
-
2020
- 2020-12-22 CN CN202011528062.1A patent/CN112700094B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109004686A (zh) * | 2018-08-29 | 2018-12-14 | 三峡大学 | 一种考虑冰蓄冷空调多模式的冷热电联供型微网系统 |
| CN110148064A (zh) * | 2019-03-08 | 2019-08-20 | 国网浙江省电力有限公司台州供电公司 | 基于激励型综合需求响应的实时新能源辅助消纳方法 |
| CN110649594A (zh) * | 2019-08-28 | 2020-01-03 | 广州供电局有限公司 | 一种基于多能协同的工业园区综合需求响应调度方法 |
| CN110807588A (zh) * | 2019-10-31 | 2020-02-18 | 国网河南省电力公司经济技术研究院 | 一种多能源耦合综合能源系统的优化调度方法 |
| CN111738529A (zh) * | 2020-07-28 | 2020-10-02 | 南方电网科学研究院有限责任公司 | 基于强化学习的综合能源系统需求响应方法、系统及设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112700094A (zh) | 2021-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112700094B (zh) | 一种基于mpc及loddlc的综合能源系统多时间尺度优化调度方法 | |
| Swaminathan et al. | Sizing and dispatch of an islanded microgrid with energy flexible buildings | |
| CN111950807B (zh) | 计及不确定性与需求响应的综合能源系统优化运行方法 | |
| CN110807588B (zh) | 一种多能源耦合综合能源系统的优化调度方法 | |
| CN109636056B (zh) | 一种基于多智能体技术的多能源微网去中心化优化调度方法 | |
| CN102509175B (zh) | 分布式供电系统可靠性优化方法 | |
| CN109634119B (zh) | 一种基于日内滚动优化的能源互联网优化控制方法 | |
| CN108229025A (zh) | 一种冷热电联供型多微网主动配电系统经济优化调度方法 | |
| CN108133285B (zh) | 一种接入大规模可再生能源的混合能源系统实时调度方法 | |
| CN117592621B (zh) | 一种虚拟电厂集群两阶段调度优化方法 | |
| CN113779783A (zh) | 计及多重不确定性的区域综合能源系统规划与运行联合优化方法 | |
| CN108875992A (zh) | 一种考虑需求响应的虚拟电厂日前优化调度方法 | |
| CN110474370B (zh) | 一种空调可控负荷、光伏储能系统的协同控制系统及方法 | |
| CN109861302B (zh) | 一种基于主从博弈的能源互联网日前优化控制方法 | |
| Wang et al. | Two-stage optimal scheduling of air conditioning resources with high photovoltaic penetrations | |
| CN116468215A (zh) | 计及源荷不确定性的综合能源系统调度方法及装置 | |
| CN117875644A (zh) | 一种综合能源系统多时间尺度优化方法及系统 | |
| CN112836882B (zh) | 考虑设备负载率变化的区域综合能源系统运行优化方法 | |
| CN120638348A (zh) | 一种分布式储能聚合缓解配电网拥塞的调度方法及系统 | |
| CN111523697A (zh) | 一种综合能源服务成本分摊与定价计算方法 | |
| CN107834574A (zh) | 一种分布式能源系统与电网交换功率的控制方法 | |
| Hu et al. | A hierarchical transactive energy management framework for optimizing the reserve profile in district energy systems | |
| Xie et al. | Collaborative scheduling of source-network-load-storage considering thermal inertia of the integrated electricity and district heating systems | |
| CN114399172A (zh) | 多场景适应多能设备组合的园区电热实时调度方法及系统 | |
| CN114021843B (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 |

























































































































