CN106886616B - 一种大规模电磁暂态电网仿真的自动分网方法 - Google Patents
一种大规模电磁暂态电网仿真的自动分网方法 Download PDFInfo
- Publication number
- CN106886616B CN106886616B CN201510930042.XA CN201510930042A CN106886616B CN 106886616 B CN106886616 B CN 106886616B CN 201510930042 A CN201510930042 A CN 201510930042A CN 106886616 B CN106886616 B CN 106886616B
- Authority
- CN
- China
- Prior art keywords
- network
- subnets
- sub
- networks
- division
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 33
- 230000001052 transient effect Effects 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000004891 communication Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 16
- 239000003990 capacitor Substances 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种大规模电磁暂态电网仿真的自动分网方法,包括不指定子网数目的自动分网和指定子网数目的自动分网,本发明在电力系统电磁暂态仿真中应用并行处理技术,可以提高仿真速度,满足大规模电力系统实时和超实时仿真对计算速度的要求。通过网络之间联络的节点作为子网络的边界点,把整个电力系统仿真网络自然的划分为各个子网,根据精度需要确定不同网络的仿真步长。本发明方法在网络划分时不需要解耦元件,网络划分简单清楚,通用性强,分网速度快。可以并行仿真不同速率、不同步长的含VSC精确模型的电磁暂态子网进行实时、超实时仿真。
Description
技术领域
本发明涉及一种电磁暂态的自动分网方法,具体涉及一种大规模电磁暂态电网仿真的自动分网方法。
背景技术
电力系统电磁暂态仿真并行方法的设计与实现既要针对电磁暂态仿真的特点,也要充分考虑硬件平台—PC机群的特点,以有效提高电磁暂态仿真的速度。由于电磁暂态实时仿真对计算时序有严格要求,考虑到PC机群的通信速度较低的特点,亟需一种大规模电磁暂态电网仿真的自动分网方法。有分网方式是手动分网,分网不均匀,计算时间长,且不是全网最优化的,这样计算过程中容易导致有的子网计算快,有的子网计算太慢,而总的计算步长必须容纳计算最慢的子网,产生木桶短板效应,而且随着电网规模的变大,容易产生子网之间数据交换不匹配、计算错误。
发明内容
为解决上述现有技术中的不足,本发明的目的是提供一种大规模电磁暂态电网仿真的自动分网方法,解决大规模电磁暂态电网仿真,在电力系统电磁暂态仿真中应用并行处理技术,可以提高仿真速度,满足大规模电力系统实时和超实时仿真对计算速度的要求。
本发明的目的是采用下述技术方案实现的:
本发明提供一种大规模电磁暂态电网仿真的自动分网方法,其改进之处在于,所述方法包括不指定子网数目的自动分网和指定子网数目的自动分网。
进一步地,所述不指定子网数目的自动分网包括:通过长线将全网分解为电气子站,对于每一个电气子站统计其计算时间和通信时间;根据网络分割目标函数进行优化,采用的自动网络分割的目标函数为:
式中,S是目标函数;na、nb、nc为各子网中的节点数目;ns为各子网间的联络线数目;
优化后将各个电气子站组合成一个个子网,原则如下:
1)按照全局最优化原则,每个子网的计算时间+通信时间最小的原则,即每个子网的计算时间+通信时间<=k1*仿真步长,其中k1为小于1的一个系数;系数保存在配置文件中;
2)对于不满足1)原则的子网,采用节点分网方式分解,将电气子站分解为若干个子网每个子网的执行时间=计算时间+通信时间+主控处理时间,其中节点分网方式的单位通信时间与长线分网方式不同,长线分网是电气量中的输电线作为不同子网之间的边界,节点分网是把母线节点作为不同子网之间的联络边界。
进一步地,所述指定子网数目的自动分网包括:通过长线将全网分解成电气子站,通过电气子站的组合或者再次用节点分网方式分解后,按下式(2)和(3)进行优化形成小于或等于k1*仿真步长的各个子网并将式(2)和式(3)进行组合得出所采用的如式(4)和式(5)所示的目标函数:
Min C=αM2+βL2 (2)
式(2)中,MinC是最优目标函数;M是所有子网中的最大节点数目;L是子网之间联络线的数目;α、β均为权重因子;αM2反映单个子网的计算时间,βL2反映单个子网间的通讯时间;
式(3)中,T是所有子网的总母线数目;ni是第i个子网中的母线数目;Fi是可行性标志,如果第i个子网是连通的,则表示第i个子网是可行的,Fi=1;否则不可行,Fi=0;Fgi、Fci、Fti表示子网中发电机、电容器和变压器与设定的边界值比较所得的逻辑值标志,当子网中的发电机数目大于设定的发电机数目时,Fgi=0,否则Fgi=1;当子网中的电容器数目大于设定的电容器数目时,Fci=0,否则Fci=1;当子网中的电容器数目大于设定的电容器数目时,当子网中的变压器数目大于设定的变压器数目时,Fti=0,否则Fti=1;β是罚因子;反映母线在子网间均匀分配的需求,b2考虑子网间联络线数目最小化,反映对动态元件数目的限制及子网是否连通;
如果形成的子网数目大于指定子网数目,提示用户修改步长或子网数目;如果形成子网数目小于指定子网数目,则将通过组合形成的子网解散,直至形成的子网数目与指定的子网数目相同,通过节点分裂形成的子网不作二次分割或组合。
与最接近的现有技术相比,本发明提供的技术方案具有的优异效果是:
本发明提出一种解决大规模电磁暂态电网仿真的自动分网方法,在电力系统电磁暂态仿真中应用并行处理技术,可以提高仿真速度,满足大规模电力系统实时和超实时仿真对计算速度的要求。通过网络之间联络的节点作为子网络的边界点,把整个电力系统仿真网络自然的划分为各个子网,根据精度需要确定不同网络的仿真步长。本发明方法在网络划分时不需要解耦元件,网络划分简单清楚,通用性强,分网速度快。可以并行仿真不同速率、不同步长的含VSC精确模型的电磁暂态子网进行实时、超实时仿真。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
以下描述和附图充分地示出本发明的具体实施方案,以使本领域的技术人员能够实践它们。其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他的改变。实施例仅代表可能的变化。除非明确要求,否则单独的组件和功能是可选的,并且操作的顺序可以变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发明的实施方案的范围包括权利要求书的整个范围,以及权利要求书的所有可获得的等同物。在本文中,本发明的这些实施方案可以被单独地或总地用术语“发明”来表示,这仅仅是为了方便,并且如果事实上公开了超过一个的发明,不是要自动地限制该应用的范围为任何单个发明或发明构思。
本发明提供一种大规模电磁暂态电网仿真的自动分网方法,包括不指定子网数目的自动分网和指定子网数目的自动分网。
不指定子网数目的自动分网包括:通过长线将全网分解为电气子站,对于每一个电气子站统计其计算时间和通信时间;根据网络分割目标函数进行优化,采用的自动网络分割的目标函数为:
式中,S是目标函数;na、nb、nc为各子网中的节点数目;ns为各子网间的联络线数目;
优化后将各个电气子站组合成一个个子网,原则如下:
1)按照全局最优化原则,每个子网的计算时间+通信时间最小的原则,即每个子网的计算时间+通信时间<=k1*仿真步长,其中k1为小于1的一个系数;系数保存在配置文件中;
2)对于不满足1)原则的子网,采用节点分网方式分解,将电气子站分解为若干个子网每个子网的执行时间=计算时间+通信时间+主控处理时间,其中节点分网方式的单位通信时间与长线分网方式不同,长线分网是电气量中的输电线作为不同子网之间的边界,而节点分网是把母线节点作为不同子网之间的联络边界。
指定子网数目的自动分网包括:通过长线将全网分解成电气子站,通过电气子站的组合或者再次用节点分网方式分解后,根据目标函数进行优化,优化后形成小于等于k1*仿真步长的各个子网,网络分割的目标函数如式(2)和式(3),并将式(2)和式(3)进行组合得出所采用的目标函数,如式(4)和式(5)所示:
Min C=αM2+βL2 (2)
式(2)中,MinC是最优目标函数;M是所有子网中的最大节点数目;L是子网之间联络线的数目;α、β均为权重因子;αM2反映单个子网的计算时间,βL2反映单个子网间的通讯时间;
式(3)中,T是所有子网的总母线数目;ni是第i个子网中的母线数目;Fi是可行性标志,如果第i个子网是连通的,则表示第i个子网是可行的,Fi=1;否则不可行,Fi=0;Fgi、Fci、Fti表示子网中发电机、电容器和变压器与设定的边界值比较所得的逻辑值标志,当子网中的发电机数目大于设定的发电机数目时,Fgi=0,否则Fgi=1;当子网中的电容器数目大于设定的电容器数目时,Fci=0,否则Fci=1;当子网中的电容器数目大于设定的电容器数目时,当子网中的变压器数目大于设定的变压器数目时,Fti=0,否则Fti=1;β是罚因子;反映母线在子网间均匀分配的需求,b2考虑子网间联络线数目最小化,反映对动态元件数目的限制及子网是否连通;
如果形成的子网数目大于指定子网数目,提示用户修改步长或子网数目;如果形成子网数目小于指定子网数目,则将通过组合形成的子网解散,直至形成的子网数目与指定的子网数目相同,通过节点分裂形成的子网不作二次分割或组合。
以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (1)
1.一种大规模电磁暂态电网仿真的自动分网方法,其特征在于,所述方法包括不指定子网数目的自动分网和指定子网数目的自动分网;
所述不指定子网数目的自动分网包括:通过长线将全网分解为电气子站,对于每一个电气子站统计其计算时间和通信时间;根据网络分割目标函数进行优化,采用的自动网络分割的目标函数为:
式中,S是目标函数;na、nb、nc为各子网中的节点数目;ns为各子网间的联络线数目;
优化后将各个电气子站组合成一个个子网,原则如下:
1)按照全局最优化原则,每个子网的计算时间+通信时间最小的原则,即每个子网的计算时间+通信时间<=k1*仿真步长,其中k1为小于1的一个系数;系数保存在配置文件中;
2)对于不满足1)原则的子网,采用节点分网方式分解,将电气子站分解为若干个子网每个子网的执行时间=计算时间+通信时间+主控处理时间,其中节点分网方式的单位通信时间与长线分网方式不同,长线分网是电气量中的输电线作为不同子网之间的边界,节点分网是把母线节点作为不同子网之间的联络边界;
所述指定子网数目的自动分网包括:通过长线将全网分解成电气子站,通过电气子站的组合或者再次用节点分网方式分解后,按下式(2)和(3)进行优化形成小于或等于k1*仿真步长的各个子网并将式(2)和式(3)进行组合得出所采用的如式(4)和式(5)所示的目标函数:
Min C=αM2+βL2 (2)
式(2)中,MinC是最优目标函数;M是所有子网中的最大节点数目;L是子网之间联络线的数目;α、β均为权重因子;αM2反映单个子网的计算时间,βL2反映单个子网间的通讯时间;
式(3)中,T是所有子网的总母线数目;ni是第i个子网中的母线数目;Fi是可行性标志,如果第i个子网是连通的,则表示第i个子网是可行的,Fi=1;否则不可行,Fi=0;Fgi、Fci、Fti表示子网中发电机、电容器和变压器与设定的边界值比较所得的逻辑值标志,当子网中的发电机数目大于设定的发电机数目时,Fgi=0,否则Fgi=1;当子网中的电容器数目大于设定的电容器数目时,Fci=0,否则Fci=1;当子网中的电容器数目大于设定的电容器数目时,当子网中的变压器数目大于设定的变压器数目时,Fti=0,否则Fti=1;β是罚因子;反映母线在子网间均匀分配的需求,b2考虑子网间联络线数目最小化,反映对动态元件数目的限制及子网是否连通;
如果形成的子网数目大于指定子网数目,提示用户修改步长或子网数目;如果形成子网数目小于指定子网数目,则将通过组合形成的子网解散,直至形成的子网数目与指定的子网数目相同,通过节点分裂形成的子网不作二次分割或组合。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510930042.XA CN106886616B (zh) | 2015-12-15 | 2015-12-15 | 一种大规模电磁暂态电网仿真的自动分网方法 |
PCT/CN2016/087821 WO2017101295A1 (zh) | 2015-12-15 | 2016-06-30 | 一种电磁暂态电网仿真的自动分网方法和计算机存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510930042.XA CN106886616B (zh) | 2015-12-15 | 2015-12-15 | 一种大规模电磁暂态电网仿真的自动分网方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106886616A CN106886616A (zh) | 2017-06-23 |
CN106886616B true CN106886616B (zh) | 2020-06-16 |
Family
ID=59055662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510930042.XA Active CN106886616B (zh) | 2015-12-15 | 2015-12-15 | 一种大规模电磁暂态电网仿真的自动分网方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN106886616B (zh) |
WO (1) | WO2017101295A1 (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108241771A (zh) * | 2017-10-27 | 2018-07-03 | 国网江苏省电力公司技能培训中心 | 基于电磁暂态仿真平台的二次设备建模方法 |
CN109190320A (zh) * | 2018-11-02 | 2019-01-11 | 贵州电网有限责任公司 | 适用于adpss双馈风机闭环试验的并行异构仿真方法 |
CN109766573B (zh) * | 2018-11-30 | 2023-04-21 | 中国电力科学研究院有限公司 | 一种新能源场站的电磁暂态实时仿真方法和装置 |
CN109783845B (zh) * | 2018-12-05 | 2023-07-18 | 中国电力科学研究院有限公司 | 边界点分群解耦的机电-电磁混合仿真分网方法及装置 |
CN109657332B (zh) * | 2018-12-14 | 2023-09-01 | 中国电力科学研究院有限公司 | 一种大规模电网电磁暂态自动建模的分网解耦方法及系统 |
CN109766586B (zh) * | 2018-12-14 | 2023-10-10 | 中国电力科学研究院有限公司 | 一种自动生成大规模电网电磁暂态仿真模型的方法及系统 |
CN109948185B (zh) * | 2019-02-01 | 2020-10-16 | 全球能源互联网研究院有限公司 | 一种电力系统的解耦仿真方法 |
CN110196997B (zh) * | 2019-04-30 | 2024-04-26 | 中国电力科学研究院有限公司 | 一种电磁暂态异步并行计算方法及装置 |
CN110380888B (zh) * | 2019-05-29 | 2021-02-23 | 华为技术有限公司 | 一种网络异常检测方法和装置 |
CN110472338B (zh) * | 2019-08-16 | 2021-07-27 | 上海交通大学 | 适用于现场可编程逻辑阵列的改进电磁暂态仿真方法 |
CN110569558B (zh) * | 2019-08-16 | 2021-09-07 | 上海交通大学 | 适用于微电网实时仿真的混合电磁暂态仿真方法 |
CN111709107B (zh) * | 2020-06-16 | 2023-03-24 | 全球能源互联网研究院有限公司 | 一种含大规模电力电子设备的电力网络仿真方法及系统 |
CN112163315B (zh) * | 2020-08-28 | 2023-03-28 | 国网天津市电力公司电力科学研究院 | 一种交直流混合电网机电-电磁暂态仿真方法 |
CN114779878B (zh) * | 2022-04-25 | 2023-03-03 | 南方电网科学研究院有限责任公司 | 一种电力系统仿真步长时间的计算系统及方法 |
CN115021256B (zh) * | 2022-08-08 | 2022-10-21 | 四川大学 | 一种大规模交直流输电系统电磁暂态模型自动生成方法 |
CN115422761B (zh) * | 2022-09-16 | 2023-10-20 | 南方电网科学研究院有限责任公司 | 一种电磁暂态仿真中并行仿真计算时间的计算方法及装置 |
CN117217002B (zh) * | 2023-09-14 | 2024-04-26 | 国网四川省电力公司电力科学研究院 | 一种基于多步长的电磁暂态并行仿真方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101221599A (zh) * | 2007-12-12 | 2008-07-16 | 国电南瑞科技股份有限公司 | 基于计算模型拼接的互联系统分布式潮流计算方法 |
US7920468B2 (en) * | 2002-03-01 | 2011-04-05 | Cisco Technology, Inc. | Method and system for constraint-based traffic flow optimisation system |
CN102508967A (zh) * | 2011-11-04 | 2012-06-20 | 南方电网科学研究院有限责任公司 | 一种大电网的电磁/机电暂态混合仿真平启动方法 |
CN102916450A (zh) * | 2012-09-28 | 2013-02-06 | 南方电网科学研究院有限责任公司 | 一种三相不对称交直流电力系统的混合实时仿真方法 |
CN104866665A (zh) * | 2015-05-19 | 2015-08-26 | 清华大学 | 基于接口等值与交互的含电力电子设备的混合仿真方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06261453A (ja) * | 1993-03-08 | 1994-09-16 | Hitachi Ltd | 電力系統シミュレーション装置 |
CN1321490C (zh) * | 2003-09-28 | 2007-06-13 | 中国电力科学研究院 | 电力系统数字仿真方法 |
CN101719182B (zh) * | 2009-12-11 | 2012-01-18 | 中国电力科学研究院 | 一种交直流电力系统分割并行电磁暂态数字仿真方法 |
CN103605852B (zh) * | 2013-11-25 | 2016-09-21 | 国家电网公司 | 一种用于大规模电网机电暂态实时仿真的并行拓扑方法 |
-
2015
- 2015-12-15 CN CN201510930042.XA patent/CN106886616B/zh active Active
-
2016
- 2016-06-30 WO PCT/CN2016/087821 patent/WO2017101295A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7920468B2 (en) * | 2002-03-01 | 2011-04-05 | Cisco Technology, Inc. | Method and system for constraint-based traffic flow optimisation system |
CN101221599A (zh) * | 2007-12-12 | 2008-07-16 | 国电南瑞科技股份有限公司 | 基于计算模型拼接的互联系统分布式潮流计算方法 |
CN102508967A (zh) * | 2011-11-04 | 2012-06-20 | 南方电网科学研究院有限责任公司 | 一种大电网的电磁/机电暂态混合仿真平启动方法 |
CN102916450A (zh) * | 2012-09-28 | 2013-02-06 | 南方电网科学研究院有限责任公司 | 一种三相不对称交直流电力系统的混合实时仿真方法 |
CN104866665A (zh) * | 2015-05-19 | 2015-08-26 | 清华大学 | 基于接口等值与交互的含电力电子设备的混合仿真方法 |
Non-Patent Citations (1)
Title |
---|
大规模电力系统机电暂态实时仿真算法及软件的研究;李亚楼;《中国优秀博硕士学位论文全文数据库(博士)工程科技II辑》;20040315(第01期);C042-34 * |
Also Published As
Publication number | Publication date |
---|---|
WO2017101295A1 (zh) | 2017-06-22 |
CN106886616A (zh) | 2017-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106886616B (zh) | 一种大规模电磁暂态电网仿真的自动分网方法 | |
CN105023058B (zh) | 一种同时考虑开关动作的配电网智能软开关运行优化方法 | |
CN104794541B (zh) | 一种基于模拟退火和锥优化的配电网运行优化方法 | |
CN102750416B (zh) | 一种含开关特性电路的电磁暂态仿真的拓扑分网方法 | |
WO2015184188A1 (en) | Systems and methods for convex relaxations and linear approximations for optimal power flow in multiphase radial networks | |
CN109494719B (zh) | 一种中低压混合配电网层次化阻抗分析方法 | |
CN105354656A (zh) | 基于分区解耦的配电网状态估计的分布式并行计算方法及系统 | |
CN110908799A (zh) | 一种分布式训练中的通信方法、装置、设备、介质 | |
CN109948185B (zh) | 一种电力系统的解耦仿真方法 | |
CN108599172A (zh) | 一种基于人工神经网络的输配网全局潮流计算方法 | |
CN105391056B (zh) | 一种考虑不平衡通信网络的电力系统分布式经济调度方法 | |
CN103795057A (zh) | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 | |
CN106886617B (zh) | 一种含多vsc的多速率电磁暂态分网方法 | |
CN110120673B (zh) | 基于戴维南等值参数辨识的分布式输配协同无功优化方法及系统 | |
CN108123450A (zh) | 一种孤立直流微网自适应下垂控制方法 | |
Li et al. | A new simplification scheme for crosstalk prediction of complex cable bundles within a cylindrical cavity | |
CN112531715B (zh) | 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法 | |
CN107947180B (zh) | 一种适用于发输配一体化潮流计算的拓扑分析方法 | |
Wang et al. | A gossip-based asynchronous distributed algorithm for economic dispatch problem with transmission losses | |
CN107977791A (zh) | 一种考虑网损的一致性算法的电力系统经济运行调度方法 | |
Chen et al. | A density-based controller placement algorithm for software defined networks | |
CN105656682B (zh) | 一种sdn控制器选址方法 | |
CN110932335B (zh) | 一种基于Petri网的船舶电力系统发电调度管理方法 | |
CN109687484B (zh) | 一种电气化铁路外部电网接入方案优化设计方法 | |
Afzalan et al. | Placement and sizing of DG using PSO&HBMO algorithms in radial distribution networks |
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 |