CN113013876A - 一种电网转动惯量计算方法 - Google Patents
一种电网转动惯量计算方法 Download PDFInfo
- Publication number
- CN113013876A CN113013876A CN202110232229.8A CN202110232229A CN113013876A CN 113013876 A CN113013876 A CN 113013876A CN 202110232229 A CN202110232229 A CN 202110232229A CN 113013876 A CN113013876 A CN 113013876A
- Authority
- CN
- China
- Prior art keywords
- power
- rotational inertia
- inertia
- power grid
- calculation method
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- 230000003068 static effect Effects 0.000 claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 238000010248 power generation Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 238000012795 verification Methods 0.000 claims description 5
- 230000006735 deficit Effects 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000011158 quantitative evaluation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- HQABUPZFAYXKJW-UHFFFAOYSA-N butan-1-amine Chemical compound CCCCN HQABUPZFAYXKJW-UHFFFAOYSA-N 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000003796 beauty Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000037213 diet Effects 0.000 description 1
- 235000005911 diet Nutrition 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明属于电网运行状态量化评估技术领域，尤其涉及一种电网转动惯量计算方法，具体是一种能源互联网转动惯量计算方法。包括转动惯量经典算法；考虑频率静态特性后的改进公式；考虑缺失功率对应的转动惯量后的最终公式；电网实际中的转动惯量计算方法验证。本发明算法配合电力系统仿真软件计算得到的电网转动惯量可有效分析评估任何运行模式下的电网频率稳定特性，具有计算准确性高、计算速度快、误差小的显著特点，可广泛推广应用。
Description
技术领域
本发明属于电网运行状态量化评估技术领域，尤其涉及一种电网转动惯量计算方法，具体是一种能源互联网转动惯量计算方法。
背景技术
近年来，国外发生了几起由系统转动惯量较小的电网设备故障诱发的大停电事故，发人深省。2016年9月28日，澳大利亚南澳州电网发生了世界上第一次由于转动惯量偏低电网的风电大规模脱网诱发的历时50小时的全州大停电事故。事故前新能源发电占比高达48.36％，严重拉低南澳电网系统转动惯量，新能源脱网后电网频率迅速下跌直至停电，停电影响170万人。2019年8月9日英国发生的大面积停电事故约有100万人受到停电影响。官方报告中指出，初步判断事故的原因可能是风电机组涉网性能偏低、系统转动惯量过低和一次调频能力不足。
如文献1：方勇杰.英国“8.9”停电事故对频率稳定控制技术的启示[J].电力系统自动化，2019，43(24):1-5.
文献2：孙华东，许涛，郭强，等.英国“8.9”大停电事故分析及对中国电网的启示[J].中国电机工程学报，2019,39(21):6183-6191。
文献3：Technical report on the events of 9August 2019[R].(2019年8月9日事故技术报告)。
文献4：Appendices to the technical report on the events of 9August2019[R]。(2019年8月9日事故技术报告附录)。
我国新能源总装机规模居世界第一位。常规机组发电占比减少造成了同步交流电网的转动惯量相对持续减小，在理论上削弱了我国电力系统承受有功冲击、频率波动和抗扰动的能力。目前，我国尚未出台统一的转动惯量计算规则。
如文献5：陈国平，李明节，许涛，刘明松.关于新能源发展的技术瓶颈研究[J].中国电机工程学报,2017,37(1):20-26.提出高占比新能源替代转动惯量导致电网抗扰动能力下降，但并未给出转动惯量具体算法。
如以下文献6-13：
文献6：李明节.大规模特高压交直流混联电网特性分析与运行控制[J].电网技术，2016，40(4)：985-991.
文献7：刘吉臻.大规模新能源电力安全高效利用基础问题[J].中国电机工程学报，2013，33(16)：1-8.
文献8：薛禹胜，雷兴，薛峰，等.关于风电不确定性对电力系统影响的评述[J].中国电机工程学报，2014，34(29)：5029-5040.
文献9：丁明，王伟胜，王秀丽，等.大规模光伏发电对电力系统影响综述[J].中国电机工程学报，2014，34(1)：1-14.
文献10：刘纯，马烁，董存，等.欧洲3·20日食对含大规模光伏发电的电网运行影响及启示[J].电网技术，2015，39(7)：1765-1772.
文献11：曾鸣，杨雍琦，李源非，等.能源互联网背景下新能源电力系统运营模式及关键技术初探[J].中国电机工程学报，2016，36(3)：681-691.
文献12：白建华，辛颂旭，刘俊，等.中国实现高比例可再生能源发展路径研究[J].中国电机工程学报，2015，35(14)：3699-3705.
文献13：刘文颖，文晶，谢昶，等.考虑风电消纳的电力系统源荷协调多目标优化方法[J].中国电机工程学报，2015，35(5)：1079-1088.
上述文献6-13提出了转动惯量较小的新能源并网给电网带来的安全稳定问题，但没有针对局域电网的转动惯量进行梳理和计算，更没有电网实际运行参数支持转动惯量理论计算值的验证。随着传统火电机组发电空间被进一步挤压，多源协同下的电力系统转动惯量特性不得而知，迫切需要开展量化摸底。
发明内容
针对上述现有技术中存在的不足之处，本发明提供了一种电网转动惯量计算方法。其目的是为了实现在计划层面能够对转动惯量进行摸底研究，在运行层面能够为实时控制策略提供参考的一种计算方法的发明目的。
本发明为实现上述目的所采用的技术方案是：
一种电网转动惯量计算方法，包括以下步骤：
步骤1：转动惯量经典算法；
步骤2：考虑频率静态特性后的改进公式；
步骤3：考虑缺失功率对应的转动惯量后的最终公式；
步骤4：电网实际中的转动惯量计算方法验证。
进一步的，所述转动惯量经典算法包括：
系统转动惯量是指整个电网或者局部电网的惯性常数，定义为H，单位是MW.s有名值或s标幺值；系统转动惯量通过人为制造功率缺额实验求得，即考虑某局域电网最严重情况，切除一台该区域已并网的具有最大转动惯量的发电机，通过求得系统频率变化计算转动惯量，简称仿真计算法；
基于仿真计算法的电力系统转动惯性常数H推导过程如下：
上述式(1)中，df/dt是RoCoF，Rate of change of frequency频率变化率；H是Synchronous inertia英美国家的系统转动惯性常数；Pg-Pc即是电力系统发电有功功率Pg和消耗的有功功率Pc之差Contingency size故障损失功率；f0表示系统故障前初始频率；如果默认f0为50Hz，得：
进一步的，所述考虑频率静态特性后的改进公式包括：
将式(2)折算，同时考虑功率缺失后的频率静态特性，得：
上述式(3)中，
是频率变化率的最大值；H是系统转动惯性常数；ΔP是故障损失功率，其中Pg-Pc是故障功率缺额，Pl是全网负荷功率，Kf是负荷频率静态特性常数；f是
最大值时对应的频率值。
进一步的，所述考虑缺失功率对应的转动惯量后的最终公式包括：
若f0为50Hz，并考虑脱网的单台发电机对应的转动惯量，则仿真计算法得到的电力系统转动惯量H为：
式中H’为损失的功率缺额对应的转动惯量，单位MW.s，由于公式(3)用的是有名值，因此公式(4)中H单位为MW.s。
进一步的，所述损失的功率缺额中指脱网的单台发电机的损失的功率缺额。
进一步的，所述电网实际中的转动惯量计算方法验证包括：
将发电机组脱网之后短时间内
最大值和脱网的发电机转动惯量的数据代入公式(4)，得到电力系统转动惯量H。
进一步的，所述所述电网实际中的转动惯量计算方法验证包括：
某电网1000MW发电机组脱网之后短时间内
最大值0.08Hz/s，脱网的1000MW发电机转动惯量约为15000MW.s，将以上数据代入公式(4)，得到电力系统转动惯量H为327500MW.s。
一种计算机存储介质，所述计算机存储介质上存有计算机程序，所述计算机程序被处理器执行时实现所述的一种电网转动惯量计算方法的步骤。
本发明具有以下有益效果及优点：
本发明算法配合电力系统仿真软件计算得到的电网转动惯量可有效分析评估任何运行模式下的电网频率稳定特性，具有计算准确性高、计算速度快、误差小等优势，可广泛推广应用。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为本发明计算方法流程图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面将结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
下面参照图1描述本发明一些实施例的技术方案。
实施例1
本发明是一种电网转动惯量计算方法，如图1所示，图1为本发明计算方法流程图。
本发明计算方法首先通过已有通用的转动惯量算法，得出50Hz下的系统转动惯性常数表达式；之后考虑功率缺失后的频率静态特性，推导出新的表达式；最后将仿真中脱网的发电机转动惯量加入总体转动惯量，得出一种基于电力系统仿真电网转动惯量算法。可在计划层面对转动惯量进行摸底研究，在运行层面为实时控制策略提供参考。
本发明一种电网转动惯量计算方法，具体包括以下步骤：
步骤1：转动惯量经典算法。
系统转动惯量是指整个电网或者局部电网的惯性常数，将其定义为H，单位是MW.s(有名值)或s(标幺值)。系统转动惯量可以通过人为制造功率缺额实验求得，即考虑某局域电网最严重情况，切除一台该区域已并网的具有最大转动惯量的发电机，通过求得系统频率变化计算转动惯量，简称仿真计算法。
以下是基于仿真计算法的电力系统转动惯性常数H推导过程：
上述式(1)中，df/dt即是RoCoF，Rate of change of frequency频率变化率；H即是Synchronous inertia英美国家的系统转动惯性常数；Pg-Pc即是电力系统发电有功功率Pg和消耗的有功功率Pc之差Contingency size故障损失功率；f0表示系统故障前初始频率；如果默认f0为50Hz，得：
步骤2：考虑频率静态特性后的改进公式。
将式(2)折算，同时还要考虑功率缺失后的频率静态特性，得：
上述式(3)中，
是频率变化率的最大值；H是系统转动惯性常数；ΔP是故障损失功率，其中Pg-Pc是故障功率缺额，Pl是全网负荷功率，Kf是负荷频率静态特性常数；f是
最大值时对应的频率值。
步骤3：考虑缺失功率对应的转动惯量后的最终公式。
若f0为50Hz，并考虑脱网的单台发电机对应的转动惯量，则仿真计算法得到的电力系统转动惯量H为：
式中H’为损失的功率缺额对应的转动惯量，单位MW.s，由于公式(3)用的是有名值，因此公式(4)中电力系统转动惯量H单位为MW.s。
其中，所述损失的功率缺额中指脱网的单台发电机的损失的功率缺额。
步骤4：电网实际中的转动惯量计算方法验证。
某电网1000MW发电机组发生脱网事件，脱网之后短时间内
最大值达到0.08Hz/s，脱网的1000MW发电机转动惯量约为15000MW.s，将以上数据代入公式(4)，得到电力系统转动惯量H为327500MW.s。
本发明中使用的技术术语如下：
电网转动惯量：转动惯量是刚体绕轴旋转时惯性的量度，电网转动惯量是特定运行方式下电力系统所有在线发电机组和负荷转动惯量之和。
实施例2
基于同一发明构思，本发明实施例还提供了一种计算机存储介质，所述计算机存储介质上存有计算机程序，所述计算机程序被处理器执行时实现实施例1所述的一种电网转动惯量计算方法的步骤。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。
Claims (8)
1.一种电网转动惯量计算方法，其特征是：包括以下步骤：
步骤1：转动惯量经典算法；
步骤2：考虑频率静态特性后的改进公式；
步骤3：考虑缺失功率对应的转动惯量后的最终公式；
步骤4：电网实际中的转动惯量计算方法验证。
2.根据权利要求1所述的一种电网转动惯量计算方法，其特征是：所述转动惯量经典算法包括：
系统转动惯量是指整个电网或者局部电网的惯性常数，定义为H，单位是MW.s有名值或s标幺值；系统转动惯量通过人为制造功率缺额实验求得，即考虑某局域电网最严重情况，切除一台该区域已并网的具有最大转动惯量的发电机，通过求得系统频率变化计算转动惯量，简称仿真计算法；
基于仿真计算法的电力系统转动惯性常数H推导过程如下：
上述式(1)中，df/dt是RoCoF，Rate of change of frequency频率变化率；H是Synchronous inertia英美国家的系统转动惯性常数；Pg-Pc即是电力系统发电有功功率Pg和消耗的有功功率Pc之差Contingency size故障损失功率；f0表示系统故障前初始频率；如果默认f0为50Hz，得：
3.根据权利要求1所述的一种电网转动惯量计算方法，其特征是：所述考虑频率静态特性后的改进公式包括：
将式(2)折算，同时考虑功率缺失后的频率静态特性，得：
上述式(3)中，
是频率变化率的最大值；H是系统转动惯性常数；ΔP是故障损失功率，其中Pg-Pc是故障功率缺额，Pl是全网负荷功率，Kf是负荷频率静态特性常数；f是
最大值时对应的频率值。
4.根据权利要求1所述的一种电网转动惯量计算方法，其特征是：所述考虑缺失功率对应的转动惯量后的最终公式包括：
若f0为50Hz，并考虑脱网的单台发电机对应的转动惯量，则仿真计算法得到的电力系统转动惯量H为：
式中H’为损失的功率缺额对应的转动惯量，单位MW.s，由于公式(3)用的是有名值，因此公式(4)中H单位为MW.s。
5.根据权利要求1所述的一种电网转动惯量计算方法，其特征是：所述损失的功率缺额中指脱网的单台发电机的损失的功率缺额。
6.根据权利要求1所述的一种电网转动惯量计算方法，其特征是：所述电网实际中的转动惯量计算方法验证包括：
将发电机组脱网之后短时间内
最大值和脱网的发电机转动惯量的数据代入公式(4)，得到电力系统转动惯量H。
7.根据权利要求1所述的一种电网转动惯量计算方法，其特征是：所述所述电网实际中的转动惯量计算方法验证包括：
某电网1000MW发电机组脱网之后短时间内
最大值0.08Hz/s，脱网的1000MW发电机转动惯量约为15000MW.s，将以上数据代入公式(4)，得到电力系统转动惯量H为327500MW.s。
8.一种计算机存储介质，其特征是：所述计算机存储介质上存有计算机程序，所述计算机程序被处理器执行时实现权利要求1-7所述的一种电网转动惯量计算方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110232229.8A CN113013876A (zh) | 2021-03-02 | 2021-03-02 | 一种电网转动惯量计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110232229.8A CN113013876A (zh) | 2021-03-02 | 2021-03-02 | 一种电网转动惯量计算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113013876A true CN113013876A (zh) | 2021-06-22 |
Family
ID=76402899
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110232229.8A Pending CN113013876A (zh) | 2021-03-02 | 2021-03-02 | 一种电网转动惯量计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113013876A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115630875A (zh) * | 2022-12-05 | 2023-01-20 | 中国电力科学研究院有限公司 | 一种同步电网的转动惯量评估方法及系统 |
-
2021
- 2021-03-02 CN CN202110232229.8A patent/CN113013876A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115630875A (zh) * | 2022-12-05 | 2023-01-20 | 中国电力科学研究院有限公司 | 一种同步电网的转动惯量评估方法及系统 |
| CN115630875B (zh) * | 2022-12-05 | 2023-03-07 | 中国电力科学研究院有限公司 | 一种同步电网的转动惯量评估方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Altin et al. | Optimization of short-term overproduction response of variable speed wind turbines | |
| CN104078999B (zh) | 一种双馈风电机组超速脱网临界时间的计算方法 | |
| CN106611965B (zh) | 预防大规模风电频繁穿越的风电场协调控制方法及系统 | |
| CN102915471A (zh) | 一种风电扰动在线安全预警的方法 | |
| CN112491043B (zh) | 一种新能源富集电网电源规划方法及系统 | |
| CN112003282B (zh) | 一种基于电网调峰能力预测新能源装机规模的方法 | |
| CN111371123B (zh) | 多类型电源汇集直流外送系统切机策略优化方法和装置 | |
| Li et al. | Transient voltage control of sending-end wind farm using a synchronous condenser under commutation failure of HVDC transmission system | |
| CN111027179B (zh) | 一种计及辅助调频服务的双馈风电场的等值建模方法 | |
| Yang et al. | Day-ahead optimal dispatch model for coupled system considering ladder-type ramping rate and flexible spinning reserve of thermal power units | |
| Weng et al. | Effects of wind power variability and intermittency on power flow | |
| CN109950929B (zh) | 一种应对大规模新能源机组低电压穿越的紧急控制方法及其系统 | |
| CN113013876A (zh) | 一种电网转动惯量计算方法 | |
| CN111130135A (zh) | 一种适用于高比例新能源接入下的电力系统惯量计算方法 | |
| Wang et al. | Development of wind-energy modeling technology and standards | |
| CN111146779A (zh) | 一种电力设备相继故障的大电网柔性安全控制方法及系统 | |
| Panwar et al. | Significance of dynamic and transient analysis in the design and operation of hybrid energy systems | |
| Hao et al. | Grid Frequency Trajectory Fitting under Transient Power Impact Caused by Highly Proportional Wind Power Fault Ride | |
| Dai et al. | Review of Frequency Response Analysis and Evaluation Methods for New Power System | |
| CN112198452B (zh) | 一种适用于工程实用化的新能源短路电流表达式的构建方法 | |
| CN114221339B (zh) | 一种高风电占比系统的暂态稳定性判断方法 | |
| Jain et al. | Variable renewable generation and grid operation | |
| Sixuan et al. | Capacity Optimization Analysis of Dispersed Wind Farm | |
| Li et al. | A disturbance source location method on the low frequency oscillation with time-varying steady-state points | |
| Sun et al. | Impacts of wind power intermittency on power system transient stability |
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 |