CN105811450A - 一种优化的中长期机组组合方法 - Google Patents

一种优化的中长期机组组合方法 Download PDF

Info

Publication number
CN105811450A
CN105811450A CN201410849635.9A CN201410849635A CN105811450A CN 105811450 A CN105811450 A CN 105811450A CN 201410849635 A CN201410849635 A CN 201410849635A CN 105811450 A CN105811450 A CN 105811450A
Authority
CN
China
Prior art keywords
unit
term
optimization
rate
long
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410849635.9A
Other languages
English (en)
Other versions
CN105811450B (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.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
North China Grid Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
North China Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, North China Grid Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201410849635.9A priority Critical patent/CN105811450B/zh
Publication of CN105811450A publication Critical patent/CN105811450A/zh
Application granted granted Critical
Publication of CN105811450B publication Critical patent/CN105811450B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种优化的中长期机组组合方法,通过负荷率和三公进度联合优化;所述方法包括:建立负荷率虚拟偏差成本,然后按照传统的中长期机组组合优化法,以三公进度均衡为目标,进行优化计算,得到机组的负荷率偏差成本;把得到机组负荷率偏差成本进行适当压缩,作为约束纳入到优化计算,重新优化计算得到第一套优化方案;调整压缩系数,重新进行优化计算,可以得到多套优化方案;决策者可以从中选择一项作为最优方案。所述方法解决了传统的中长期机组组合方法中仅对三公进度进行优化控制,对负荷率却不能精细优化的问题。

Description

一种优化的中长期机组组合方法
技术领域:
本发明涉及优化的中长期机组组合方法,更具体涉及一种负荷率和三公进度联合优化的中长期机组组合方法。
背景技术:
我国传统的中长期机组组合通常以三公进度均衡为目标进行优化计算,由于优化计算中需要统筹兼顾电力、电量平衡、电网安全约束等,考虑实际生产中机组的运行特性,模型已经很复杂,对机组的负荷率难以再进行精细化考虑。但是机组的负荷率水平往往对机组的污染排放、发电成本有很大的影响,机组的负荷率水平偏高,机组的煤耗成本低、污染排放也会降低,但是实际生产中一般不会安排机组满发,因为如果机组满发,则会降低系统备用,所以如何在兼顾系统电力、电量平衡、系统备用、网络安全等约束的前提下,实现负荷率和三公进度联合优化,是一个亟待解决的问题。本发明的方法提出了负荷率和三公进度联合优化中长期机组组合方法,是电网公司安排中长期机组组合方式的辅助工具,为电网降低污染排放、提高调度管理水平提供有效的技术支持。
发明内容:
本发明的目的是提供一种优化的中长期机组组合方法,所述方法解决了传统的中长期机组组合方法中仅对三公进度进行优化控制,却不能对负荷率精细优化的问题。
为实现上述目的,本发明采用以下技术方案:一种优化的中长期机组组合方法,通过负荷率和三公进度联合优化;所述方法包括以下步骤:
(1)建立负荷率虚拟偏差成本;
(2)通过中长期机组组合方式进行优化得到机组的负荷率偏差成本;
(3)采用约束法进行多目标优化;
(4)联合优化。
本发明提供的一种优化的中长期机组组合方法,所述步骤(1)的偏差成本为机组的负荷率与理想负荷率之间存在偏差产生并通过下式建立:
R(i,t)-R'(i,t)=Δ+ R(i,t)-Δ- R(i,t)(1)
RatioBia=(Δ+ R(i,t)+Δ- R(i,t))*Inccost(i,t)(2)
其中,R′(i,t)表示机组i在t时段的理想负荷率;R(i,t)表示机组i在t时段的负荷率;Δ+ R(i,p)、Δ- R(i,p)分别表示机组i的t时段的负荷率正偏差和负偏差;RatioBia表示机组i在时段t的负荷率偏差;Inccost(i,t)为负荷率偏差微增成本。
本发明提供的一种优化的中长期机组组合方法,所述步骤(2)中的中长期机组组合方式包括以下步骤:
三公进度均衡目标;
电力系统约束。
本发明提供的另一优选的一种优化的中长期机组组合方法,所述三公进度均衡目标通过下式确定:
min f = Σ i = 1 N | Σ t = 1 T E ( i , t ) - E g ( i ) | E g ( i ) - - - ( 3 )
其中,E(i,t)表示机组i在t时段的优化电量;Eg(i)表示机组i的三公计划电量;f表示所有机组的三公进度偏差。
本发明提供的再一优选的一种优化的中长期机组组合方法,所述电力系统约束包括电力平衡约束、电量平衡约束、电力电量关联约束、电力系统备用约束、机组出力上下限约束、机组最小运行时间和机组最小停运时间约束、机组爬坡约束和电网安全约束。
本发明提供的又一优选的一种优化的中长期机组组合方法,所述电力平衡约束通过下式(4)确定:
Σ i I p i ( t ) = p d ( t ) , t = 1,2 , . . . . . T - - - ( 4 )
其中:pi(t)表示机组i在t时段的物理出力;pd(t)表示t时段系统的电力负荷;
所述电量平衡约束通过下式(5)确定:
Σ i I E ( i , t ) = E d ( t ) , t = 1,2 , . . . . . T - - - ( 5 )
其中:Ed(t)表示t时段的预估电量负荷;
所述电力电量关联约束通过下式(6)确定:
E(i,t)=pi(t)*R(i,t)*PrdMin(t)/60(6)
其中:R(i,t)表示机组i在t时段的负荷率;PrdMin(t)表示t时段的优化分钟数;
所述系统备用约束通过下式(7)和(8)确定:
Σ i I ( p max ( i , t ) * u i ( t ) - p i ( t ) ) ≥ R up ( t ) - - - ( 7 )
Σ i = 1 I ( p i ( t ) - p min ( i , t ) * u i ( t ) ) ≥ R down ( t ) - - - ( 8 )
其中,Pmax(i,t)和Pmin(i,t)分别为机组i在t时刻的最大出力和最小出力;ui(t)为机组i在t时刻的运行状态,0表示停机,1表示运行;Rup(t)和Rdown(t)分别为系统在t时刻的上旋、下旋备用需求;
所述机组出力上下限约束通过下式(9)确定:
p i ‾ ( t ) ≤ p i ( t ) ≤ p i ‾ u i ( t ) - - - ( 9 )
其中,p i 分别为机组i输出功率的上下限;
所述机组最小运行时间和机组最小停运时间约束通过下式(10)确定:
u i ( t ) = 1,1 ≤ x i ( t ) ≤ τ i ‾ u i ( t ) = - 1 , - τ ‾ i ≤ x i ( t ) ≤ - 1 - - - ( 10 )
其中,xi(t)为机组i在t时刻已经连续运行的时间或停运的时间;τ i 分别为机组i最小运行时间和停运时间
所述机组爬坡约束通过下式(11)确定:
i≤pi(t)-pi(t-1)≤Δi(11)
其中,pi(t-1)为机组i在t-1时刻的有功功率,Δi为机组i每时刻可加减负荷的最大值;-Δi为机组i每时刻可加减负荷的最小值;
所述电网安全约束通过下式(12)确定:
P ij ( t ) ≤ P ij ‾ - - - ( 12 )
其中,Pij(t)为支路ij在t时刻的潮流功率,为支路ij在t时刻的潮流功率上限。
本发明提供的又一优选的一种优化的中长期机组组合方法,所述步骤(3)的优化过程为把机组负荷率偏差成本RatioBia进行压缩,形成约束纳入优化计算。
本发明提供的又一优选的一种优化的中长期机组组合方法,所述压缩过程通过下式确定:
RatioBia≤B*factor(13)
0≤factor≤1(14)
其中:B为以三公进度均衡为目标计算得到的机组负荷率偏差值;
factor表示在机组负荷率偏差成本的压缩值。
本发明提供的又一优选的一种优化的中长期机组组合方法,所述步骤(4)的优化过程为为把机组负荷率偏差成本约束纳入到所述步骤(2)中长期机组组合进行重新优化计算,得到三公进度和负荷率的联合优化结果。
本发明提供的又一优选的一种优化的中长期机组组合方法,通过调整压缩系数factor,得到多套不同的联合优化方案,调度决策者从中选择最佳方案。
和最接近的现有技术比,本发明提供技术方案具有以下优异效果
1、本发明的方法,通过引入负荷率实现了中长期电力电量联合优化,满足中长期机组组合生产需求;
2、本发明的方法,确保机组三公进度均衡的同时,还实现了负荷率优化,降低了污染排放和机组的运行成本。;
3、本发明的方法避免了负荷率区间的强制性约束,优化方法容易求解;
4、本发明的方法可以通过调整压缩参数得到多套方案进行比对分析;
5、本发明的方法是电网公司安排中长期机组组合方式的辅助工具,为电网降低污染排放、提高调度管理水平提供有效的技术支持。
附图说明
图1为本发明的机组负荷率偏差与偏差成本关系图;
图2为本发明的方法流程图。
具体实施方式
下面结合实施例对发明作进一步的详细说明。
实施例1:
如图1-2所示,本例的发明提供一种负荷率和三公进度联合优化中长期机组组合方法;步骤包括首先建立负荷率虚拟偏差成本,然后按照中长期机组组合方式进行优化计算,得到机组的负荷率偏差成本;把得到机组负荷率偏差成本进行适当压缩,作为约束纳入到优化计算,得到第一套优化方案;调整压缩系数,重新进行优化计算,可以得到多套优化方案;决策者可以从中选择一项作为最优方案。
1、负荷率偏差虚拟成本定义
如图1所示,针对机组的负荷率引入虚拟偏差成本,机组的负荷率与理想负荷率(或理想负荷率区间)之间存在偏差,都会产生偏差成本。
R(i,t)-R'(i,t)=Δ+ R(i,t)-Δ- R(i,t)(1)
RatioBia=(Δ+ R(i,t)+Δ- R(i,t))*Inccost(i,t)(2)
R'(i,t)分别表示机组i在t时段的理想负荷率。
R(i,t)表示机组i在t时段的负荷率。
Δ+ R(i,p)、Δ- R(i,p)分别表示机组i的t时段的负荷率正偏差和负偏差;
RatioBia表示机组i在时段t的负荷率偏差;
2、中长期机组组合
1)三公进度均衡目标
min f = Σ i = 1 N | Σ t = 1 T E ( i , t ) - E g ( i ) | E g ( i ) - - - ( 3 )
E(i,t)表示机组i在t时段的优化电量;Eg(i)表示机组i的三公计划电量;f表示所有机组的三公进度偏差;
2)系统约束
电力平衡约束:
Σ i I p i ( t ) = p d ( t ) , t = 1,2 , . . . . . T - - - ( 4 )
其中:pi(t)表示机组i在t时段的物理出力;pd(t)表示t时段系统的电力负荷;
电量平衡约束:
Σ i I E ( i , t ) = E d ( t ) , t = 1,2 , . . . . . T - - - ( 5 )
其中:Ed(t)表示t时段的预估电量负荷;
电力电量关联约束:
E(i,t)=pi(t)*R(i,t)*PrdMin(t)/60(6)
其中:PrdMin(t)表示t时段的优化分钟数;
系统备用约束:
Σ i I ( p max ( i , t ) * u i ( t ) - p i ( t ) ) ≥ R up ( t ) - - - ( 7 )
Σ i = 1 I ( p i ( t ) - p min ( i , t ) * u i ( t ) ) ≥ R down ( t ) - - - ( 8 )
其中,Pmax(i,t)和Pmin(i,t)分别为机组i在t时刻的最大出力和最小出力;Rup(t)和Rdw(t)分别为系统在t时刻的上旋、下旋备用需求;
机组出力上下限约束:
p i ‾ ( t ) ≤ p i ( t ) ≤ p i ‾ u i ( t ) - - - ( 9 )
其中,p i 分别为机组i输出功率的上下限;
机组最小运行时间和机组最小停运时间约束:
u i ( t ) = 1,1 ≤ x i ( t ) ≤ τ i ‾ u i ( t ) = - 1 , - τ ‾ i ≤ x i ( t ) ≤ - 1 - - - ( 10 )
其中,τ i 分别为机组i最小运行时间和停运时间
机组爬坡约束:
i≤pi(t)-pi(t-1)≤Δi(11)
其中,pi(t-1)为机组i在t-1时刻的有功功率,Δi为机组i每时刻可加减负荷的最大值
电网安全约束:
P ij ( t ) ≤ P ij ‾ - - - ( 12 )
其中,Pij(t)为支路ij在t时刻的潮流功率,为支路ij在t时刻的潮流功率上限;
按照三公进度均衡目标,考虑电力、电量平衡、系统备用、电网安全、机组出力上下限、爬坡率、机组最小运行时间和机组最小停运时间等约束,应用COMPLEX商用算法包,采用混合整数规划法进行优化计算,得到优化结果,计算机组负荷率偏差成本RatioBia的最大值B。
3、采用约束法进行多目标优化
把机组负荷率偏差成本RatioBia进行适当的压缩,形成约束纳入优化计算。这里避免将负荷率按照负荷率区间进行强制性约束,这种强制性约束难以设定合理的负荷率区间,并且容易无解。
RatioBia≤B*factor(13)
0≤factor≤1(14)
其中:factor表示在机组负荷率偏差成本的压缩值。
4、联合优化计算
把机组负荷率偏差成本约束纳入到传统的中长期机组组合重新进行优化计算,得到三公进度和负荷率的联合优化结果。可以调整压缩系数factor,得到多套不同的联合优化方案,调度决策者可以从中选择最佳方案。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员尽管参照上述实施例应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。

Claims (10)

1.一种优化的中长期机组组合方法,其特征在于:通过负荷率和三公进度联合优化;所述方法包括以下步骤:
(1)建立负荷率虚拟偏差成本;
(2)通过中长期机组组合方式进行优化得到机组的负荷率偏差成本;
(3)采用约束法进行多目标优化;
(4)联合优化。
2.如权利要求1所述的一种优化的中长期机组组合方法,其特征在于:所述步骤(1)的偏差成本为机组的负荷率与理想负荷率之间存在偏差产生并通过下式建立:
R(i,t)-R'(i,t)=Δ+ R(i,t)-Δ- R(i,t)(1)
RatioBia=(Δ+ R(i,t)+Δ- R(i,t))*Inccost(i,t)(2)
其中,R′(i,t)表示机组i在t时段的理想负荷率;R(i,t)表示机组i在t时段的负荷率;Δ+ R(i,p)、Δ- R(i,p)分别表示机组i的t时段的负荷率正偏差和负偏差;RatioBia表示机组i在时段t的负荷率偏差;Inccost(i,t)为负荷率偏差微增成本。
3.如权利要求1所述的一种优化的中长期机组组合方法,其特征在于:所述步骤(2)中的中长期机组组合方式包括以下步骤:
三公进度均衡目标;
电力系统约束。
4.如权利要求3所述的一种优化的中长期机组组合方法,其特征在于:所述三公进度均衡目标通过下式确定:
min f = Σ i = 1 N | Σ t = 1 T E ( i , t ) - E g ( i ) | E g ( i ) - - - ( 3 )
其中,E(i,t)表示机组i在t时段的优化电量;Eg(i)表示机组i的三公计划电量;f表示所有机组的三公进度偏差。
5.如权利要求3所述的一种优化的中长期机组组合方法,其特征在于:所述电力系统约束包括电力平衡约束、电量平衡约束、电力电量关联约束、电力系统备用约束、机组出力上下限约束、机组最小运行时间和机组最小停运时间约束、机组爬坡约束和电网安全约束。
6.如权利要求5所述的一种优化的中长期机组组合方法,其特征在于:所述电力平衡约束通过下式(4)确定:
Σ i I p i ( t ) = p d ( t ) , t = 1,2 , . . . . . T - - - ( 4 )
其中:pi(t)表示机组i在t时段的物理出力;pd(t)表示t时段系统的电力负荷;
所述电量平衡约束通过下式(5)确定:
Σ i I E ( i , t ) = E d ( t ) , t = 1,2 , . . . . . T - - - ( 5 )
其中:Ed(t)表示t时段的预估电量负荷;
所述电力电量关联约束通过下式(6)确定:
E(i,t)=pi(t)*R(i,t)*PrdMin(t)/60(6)
其中:R(i,t)表示机组i在t时段的负荷率;PrdMin(t)表示t时段的优化分钟数;
所述系统备用约束通过下式(7)和(8)确定:
Σ i I ( p max ( i , t ) * u i ( t ) - p i ( t ) ) ≥ R up ( t ) - - - ( 7 )
Σ i I ( p i ( t ) - p min ( i , t ) * u i ( t ) ) ≥ R down ( t ) - - - ( 8 )
其中,Pmax(i,t)和Pmin(i,t)分别为机组i在t时刻的最大出力和最小出力;ui(t)为机组i在t时刻的运行状态,0表示停机,1表示运行;Rup(t)和Rdown(t)分别为系统在t时刻的上旋、下旋备用需求;
所述机组出力上下限约束通过下式(9)确定:
p i ‾ u i ( t ) ≤ p i ( t ) ≤ p i ‾ u i ( t ) - - - ( 9 )
其中,p i 分别为机组i输出功率的上下限;
所述机组最小运行时间和机组最小停运时间约束通过下式(10)确定:
u i ( t ) = 1,1 ≤ x i ( t ) ≤ τ i ‾ u i ( t ) = - 1 , - τ ‾ i ≤ x i ( t ) ≤ - 1 - - - ( 10 )
其中,xi(t)为机组i在t时刻已经连续运行的时间或停运的时间;τ i 分别为机组i最小运行时间和停运时间
所述机组爬坡约束通过下式(11)确定:
i≤pi(t)-pi(t-1)≤Δi(11)
其中,pi(t-1)为机组i在t-1时刻的有功功率,Δi为机组i每时刻可加减负荷的最大值;-Δi为机组i每时刻可加减负荷的最小值;
所述电网安全约束通过下式(12)确定:
P ij ( t ) ≤ P ij ‾ - - - ( 12 )
其中,Pij(t)为支路ij在t时刻的潮流功率,为支路ij在t时刻的潮流功率上限。
7.如权利要求2所述的一种优化的中长期机组组合方法,其特征在于:所述步骤(3)的优化过程为把机组负荷率偏差成本RatioBia进行压缩,形成约束纳入优化计算。
8.如权利要求7所述的一种优化的中长期机组组合方法,其特征在于:所述压缩过程通过下式确定:
RatioBia≤B*factor(13)
0≤factor≤1(14)
其中:B为以三公进度均衡为目标计算得到的机组负荷率偏差值;factor表示在机组负荷率偏差成本的压缩值。
9.如权利要求8所述的一种优化的中长期机组组合方法,其特征在于:所述步骤(4)的优化过程为为把机组负荷率偏差成本约束纳入到所述步骤(2)中长期机组组合进行重新优化计算,得到三公进度和负荷率的联合优化结果。
10.如权利要求9所述的一种优化的中长期机组组合方法,其特征在于:通过调整压缩系数factor,得到多套不同的联合优化方案,调度决策者从中选择最佳方案。
CN201410849635.9A 2014-12-29 2014-12-29 一种优化的中长期机组组合方法 Active CN105811450B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410849635.9A CN105811450B (zh) 2014-12-29 2014-12-29 一种优化的中长期机组组合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410849635.9A CN105811450B (zh) 2014-12-29 2014-12-29 一种优化的中长期机组组合方法

Publications (2)

Publication Number Publication Date
CN105811450A true CN105811450A (zh) 2016-07-27
CN105811450B CN105811450B (zh) 2019-01-04

Family

ID=56421504

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410849635.9A Active CN105811450B (zh) 2014-12-29 2014-12-29 一种优化的中长期机组组合方法

Country Status (1)

Country Link
CN (1) CN105811450B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002019209A2 (en) * 2000-08-25 2002-03-07 Sca Technologies, Llc Production and distribution supply chain optimization software
CN102097866A (zh) * 2011-03-28 2011-06-15 国电南瑞科技股份有限公司 中长期机组组合优化方法
CN103337039A (zh) * 2013-04-28 2013-10-02 国家电网公司 一种多目标多约束的中长期购电决策方法
CN103345712A (zh) * 2013-06-19 2013-10-09 国家电网公司 一种中长期跨区跨省交易优化与模式评估方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002019209A2 (en) * 2000-08-25 2002-03-07 Sca Technologies, Llc Production and distribution supply chain optimization software
CN102097866A (zh) * 2011-03-28 2011-06-15 国电南瑞科技股份有限公司 中长期机组组合优化方法
CN103337039A (zh) * 2013-04-28 2013-10-02 国家电网公司 一种多目标多约束的中长期购电决策方法
CN103345712A (zh) * 2013-06-19 2013-10-09 国家电网公司 一种中长期跨区跨省交易优化与模式评估方法

Also Published As

Publication number Publication date
CN105811450B (zh) 2019-01-04

Similar Documents

Publication Publication Date Title
CN104134994B (zh) 一种利用agc进行稳定断面潮流越限的校正控制方法
CN104578176A (zh) 一种计及直流互动的发电计划编制方法
CN106458618A (zh) 用于电网的间歇式可再生能量的负载平衡的系统和方法
CN109118024B (zh) 一种考虑火电机组多阶段状态转移的电力系统多资源调度方法
CN104505864B (zh) 用于消纳分布式发电的需求响应控制策略仿真系统和方法
CN105305419A (zh) 一种含有压缩空气储能的独立微网容量优化配置方法
CN103066626A (zh) 一种多源协调日前发电计划方法
CN106786799A (zh) 一种直流联络线功率阶梯化发电计划优化方法
CN106505633A (zh) 一种风光接入容量确定方法及装置
CN111030110A (zh) 一种考虑电转气消纳风电的电力-天然气耦合系统鲁棒协同调度方法
CN108879657B (zh) 基于风电容量可信度的电力电量平衡优化方法
CN106340893B (zh) 一种基于虚拟发电厂的分布式光伏发电功率控制方法
CN108429280B (zh) 一种无源电网广域虚拟频率控制方法及系统
CN107769266A (zh) 一种多时间尺度发电与备用联合优化方法
US11262713B2 (en) Method for calculating control parameters of heating supply power of heating network
CN104035475B (zh) 兼顾电网安全和弃光最小的光伏电站有功自动控制方法
CN104239960A (zh) 考虑抽水蓄能机组的发电计划优化方法
CN105098839A (zh) 一种基于风电不确定出力的风电并网协调优化方法
CN117254503A (zh) 基于风电消纳的源-网-储-荷双层协同低碳调度方法、设备及介质
CN106877382B (zh) 一种基于直流功率跟随风电功率波动的发电机组组合方法
CN105811450A (zh) 一种优化的中长期机组组合方法
CN112242710A (zh) 一种基于场景分析的新能源跨区域消纳方法和系统
CN112491091B (zh) 一种大规模区域电网日前安全约束机组组合快速求解方法
CN104753086B (zh) 一种满足多种场景鲁棒运行的最小旋转备用调度方法
CN108288854A (zh) 一种源网荷协调控制方法及系统

Legal Events

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