CN117556641A - 基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法 - Google Patents

基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法 Download PDF

Info

Publication number
CN117556641A
CN117556641A CN202410040594.2A CN202410040594A CN117556641A CN 117556641 A CN117556641 A CN 117556641A CN 202410040594 A CN202410040594 A CN 202410040594A CN 117556641 A CN117556641 A CN 117556641A
Authority
CN
China
Prior art keywords
boost
full
simulation
bridge converter
electromagnetic transient
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
CN202410040594.2A
Other languages
English (en)
Other versions
CN117556641B (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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN202410040594.2A priority Critical patent/CN117556641B/zh
Publication of CN117556641A publication Critical patent/CN117556641A/zh
Application granted granted Critical
Publication of CN117556641B publication Critical patent/CN117556641B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid distribution networks

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)
  • Dc-Dc Converters (AREA)

Abstract

本发明涉及电力系统技术领域,公开了基于Boost‑全桥变流器解耦的恒定等效电导电磁暂态建模方法,包括:采集新能源并网系统的参数,所述的参数包括太阳能光伏电池阵列、风力发电机、储能装置、电力电子设备的电阻、电感和电容数量以及电阻值、电容值、电感值等数据;将新能源并网系统分成新能源部分和Boost‑全桥变流器部分;对新能源部分进行电磁暂态仿真;对Boost‑全桥变流器部分进行等效电导的电磁暂态建模,保持导纳矩阵恒定,确保整个系统仿真的快速性以及稳定性。通过本发明所提供的技术方案,可以实现提高仿真效率,同时保持足够的精度。

Description

基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模 方法
技术领域
本发明涉及电力系统技术领域,具体是基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法。
背景技术
随着可再生能源的大规模开发及特高压交直流输电技术的迅速发展,新一代电力系统呈现出电源种类丰富、电能转换形式复杂等特征。在发电侧,风电、光伏等大规模可再生能源通过变流器集中并入传统交流电网;在输电侧,西部能源基地通过特高压直流(UHVDC)及特高压交流(UHVAC)混合组网的方式对三华地区送电;在配电侧,大量分布式电源通过变流器接口组成微电网系统,逐步形成功率双向流动的主动配电网。目前,电力电子变流器作为连接交直流电网的能量变换装置已在电力系统中获得广泛应用,发电、输电和配电各环节 “电力电子化”的趋势日益显著。
电力电子技术在能源转换和控制领域的广泛应用促使了对电力电子设备的建模和仿真方法不断发展。Boost-全桥变流器作为电力电子系统中的核心组件,对于实现能量转换、电压控制和频率调节等任务至关重要。在这一领域,电磁暂态建模和仿真技术是解决电力电子系统设计、分析和优化问题的关键。
在电力电子领域,研究者们正在寻求新的建模和仿真方法,以解决传统方法所面临的挑战。传统的方法在模拟系统时可能存在复杂性和不稳定性问题。等效电导电磁暂态建模方法已成为一个备受关注的研究领域。该方法旨在降低计算复杂性,提高仿真效率,同时保持足够的精度。本发明旨在解决这些问题,提供一种新的方法,基于Boost-全桥变流器解耦,以改善新能源并网系统的电磁暂态建模仿真的情况。
发明内容
本发明的目的在于克服现有技术的不足,提供基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法,包括如下步骤:
步骤一,采集新能源并网系统的参数,所述的参数包括太阳能光伏电池阵列、风力发电机、储能装置、电力电子设备的电阻、电感和电容数量以及电阻值、电容值、电感值的数据;
步骤二,将新能源并网系统分成新能源部分和Boost-全桥变流器部分;
步骤三,分别对新能源部分和Boost-全桥变流器部分进行仿真,包括对新能源部分进行电磁暂态仿真和对Boost-全桥变流器部分进行等效电导的电磁暂态建模并进行仿真;根据新能源部分和Boost-全桥变流器部分的仿真结果,得到离散化模型的输出结果。
进一步的,所述的对Boost-全桥变流器部分进行等效电导的电磁暂态建模,包括对Boost-全桥变流器部分中的电容与电感分别进行动态建模,将Boost-全桥变流器部分中的含有电容部分与含有电感部分相互解耦,得到包含电感的离散化模型与包含电容的离散化模型 ,在仿真时间T内,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值;仿真时间结束时,将得到的桥臂电流以及电容电压作为Boost-全桥变流器的输出结果。
进一步的,所述的在仿真时间T内,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值,包括:
设当前时刻为0,当前时刻的电感电流已经根据初值给出;更新时刻的电容电压;仿真时刻再推进半个步长至/>;再计算出/>时刻电感电流/>;每个时步的计算流程都与之相同,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值。其中,/>为仿真步长。
进一步的,还包括在电磁暂态仿真过程中,用等效电阻来表示开关管,等效电阻即为二值电阻,开关关断为大电阻,导通则为小电阻。
本发明的有益效果是:本发明通过对Boost-全桥变流器进行电磁暂态建模,旨在降低计算复杂性,提高仿真效率,同时保持足够的精度。
附图说明
图1为基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法的流程示意图;
图2为光伏发电系统拓扑图;
图3为光伏发电系统解耦等效电路示意图;
图4为Boost-全桥变流器的恒导纳电磁暂态建模仿真方法流程图;
图5为Boost变换器等效电路
图6为全桥变流器等效电路示意图;
图7为对新能源部分进行电磁暂态仿真的流程示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
为了使本发明的目的,技术方案及优点更加清楚明白,结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。需要说明的是,术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。
而且,术语“包括”,“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程,方法,物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程,方法,物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程,方法,物品或者设备中还存在另外的相同要素。
以下结合实施例对本发明的特征和性能作进一步的详细描述。
如图1所示,基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法,包括如下步骤:
步骤一,采集新能源并网系统的参数,所述的参数包括太阳能光伏电池阵列、风力发电机、储能装置、电力电子设备的电阻、电感和电容数量以及电阻值、电容值、电感值的数据;
步骤二,将新能源并网系统分成新能源部分和Boost-全桥变流器部分;
步骤三,分别对新能源部分和Boost-全桥变流器部分进行仿真,包括对新能源部分进行电磁暂态仿真,如图7所示,和对Boost-全桥变流器部分进行等效电导的电磁暂态建模并进行仿真;根据新能源部分和Boost-全桥变流器部分的仿真结果,得到离散化模型的输出结果。
所述的对Boost-全桥变流器部分进行等效电导的电磁暂态建模,包括对Boost-全桥变流器部分中的电容与电感分别进行动态建模,将Boost-全桥变流器部分中的含有电容部分与含有电感部分相互解耦,得到包含电感的离散化模型与包含电容的离散化模型 ,在仿真时间T内,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值;仿真时间结束时,将得到的桥臂电流以及电容电压作为Boost-全桥变流器的输出结果。
所述的在仿真时间T内,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值,包括:
设当前时刻为0,当前时刻的电感电流已经根据初值给出;更新时刻的电容电压;仿真时刻再推进半个步长至/>;再计算出/>时刻电感电流/>;每个时步的计算流程都与之相同,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值。其中,/>为仿真步长。
进一步的,还包括在电磁暂态仿真过程中,用等效电阻来表示开关管,等效电阻即为二值电阻,开关关断为大电阻,导通则为小电阻。
具体的, 本发明提供基于Boost-全桥变流器解耦的新能源并网系统恒定等效电导电磁暂态建模仿真方法。顾名思义就是将新能源系统各部分相互解耦,再依次进行仿真计算。如图2为整个新能源系统的拓扑图(以光伏发电单元为例)。
将图2所示的光伏发电系统电力电子部分进行相互的解耦,如图3所示,其中黑色虚线框表示电感部分的解耦等效模型,蓝色实线框表示电容部分解耦等效模型。
本发明基于Boost-全桥变流器恒定等效电导电磁暂态建模仿真方法将系统分解开,电源侧与电力电子部分相互解耦,分模块进行计算。
图4为本发明提供的基于Boost-全桥变流器解耦的新能源并网系统恒定等效电导电磁暂态建模仿真方法中流程图。如图所示,本发明所提供的电磁暂态建模方法包括:
步骤1、对上述Boost-全桥变流器模型中的电容与电感分别进行动态建模,得到包含电感的离散化模型与包含电容的离散化模型;
步骤2、在仿真时间内,按照交替计算方式计算上述电感与电容离散化模型中的电流与电压值;
按照交替计算方式计算上述电感与电容离散化模型,具体包括:
设当前时刻为,当前时刻的电感电流已经在上一时步求得;更新/>时刻的电容电压/>;仿真时刻再推进半个步长至/>;再计算出/>时刻电感电流/>。其中,/>为仿真步长。
通过电感与电容错开半个步长即的时间,电感电流的动态方程与电容电压的动态方程相互解耦。
步骤3、仿真时间结束,将所求得的桥臂电流以及电容电压/>作为整个离散化模型的输出结果。
具体的,在电磁暂态仿真过程中,可以用等效电阻来表示开关管,该等效电阻即为二值电阻,开关关断为大电阻,导通则为小电阻。
图5所示为将Boost子模块中的IGBT与二极管简化成电阻后的电路图。
在左边端口对列写电压方程可得:
其中L 0为图5中的桥臂电感,R 1R 2为IGBT/二极管开关组的等效电阻,为左侧1、2两端口间的电压差,/>为电容两端电压,/>为流入电感的电流
在原式中设
可以将上式改写为
以下对该式进行至/>积分
利用梯形积分法将上式离散化为:
式中,表示积分步长;V c (t)为t时刻电容两端电压;/>为/>时刻电容两端电压;
梯形积分法与中矩形积分法有相似的精度,在上式中将部分近似为/>
式中,i arm (t)为t时刻流入电感的电流;为/>时刻流入电感的电流。
将上式转化为诺顿等效电路的形式
其中
V arm (t)为t时刻1、2两端口间的电压差。
对右端口的电容列写电流方程
其中C表示图5中端口电容值;i o为流出端口3的电流;
;/>
以下对该式进行至/>积分,并利用梯形积分法将上式离散化为
由于上式化简为:
梯形积分法与中矩形积分法有相似的精度,在上式中将部分近似为
将上式转化为诺顿等效电路的形式
其中
图6为将上图全桥变流器结构图中的IGBT与二极管简化成电阻后的电路图
列写电压方程
式中,L a 、L b 、L c 表示图6中的3个桥臂电感;R 1R 2R 3R 4R 5R 6为对应的六个IGBT/二极管开关组的等效电阻;i a 、i b 、i c 表示流入桥臂的电流;V c2 表示图6中右侧电容两端电压;V a 、V b 、V c 表示三个桥臂的端口a、b、c的端口对地电压;i d 表示流入端口d的电流。
在原式中设
电感方程化简为
按照半隐式延迟解耦方法上式化为
为例对电感方程化简为:
其中将化为/>
电容方程化简为(其中将化为/>):
其中,设
将表达式转化为诺顿等效电路的形式
在进行新能源并网系统电磁暂态系统计算时,新能源系统包含大量的非线性电阻以及电力电子变换器,采用传统的详细模型建模,步骤繁琐,耗时较长。本发明将Boost-全桥变流器简化为电阻和电流源并联的诺顿等效电路,实现系统的解耦,新能源部分与电力电子部分相互分开并行计算,大大提升计算的效率;同时,该方法避免了大量的电导矩阵的LU分解,保证系统的节点导纳矩阵不变,极大提升了仿真速度。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。

Claims (4)

1.基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法,其特征在于,包括如下步骤:
步骤一,采集新能源并网系统的参数,所述的参数包括太阳能光伏电池阵列、风力发电机、储能装置、电力电子设备的电阻、电感和电容数量以及电阻值、电容值、电感值的数据;
步骤二,将新能源并网系统分成新能源部分和Boost-全桥变流器部分;
步骤三,分别对新能源部分和Boost-全桥变流器部分进行仿真,包括对新能源部分进行电磁暂态仿真和对Boost-全桥变流器部分进行等效电导的电磁暂态建模并进行仿真;根据新能源部分和Boost-全桥变流器部分的仿真结果,得到离散化模型的输出结果。
2.根据权利要求1所述的基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法,其特征在于,所述的对Boost-全桥变流器部分进行等效电导的电磁暂态建模,包括对Boost-全桥变流器部分中的电容与电感分别进行动态建模,将Boost-全桥变流器部分中的含有电容部分与含有电感部分相互解耦,得到包含电感的离散化模型与包含电容的离散化模型 ,在仿真时间T内,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值;仿真时间结束时,将得到的桥臂电流以及电容电压作为Boost-全桥变流器的输出结果。
3.根据权利要求2所述的基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法,其特征在于,所述的在仿真时间T内,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值,包括:
设当前时刻为0,当前时刻的电感电流已经根据初值给出;更新时刻的电容电压/>;仿真时刻再推进半个步长至/>;再计算出/>时刻电感电流/>;每个时步的计算流程都与之相同,按照交替计算方式分别得到包含电感的离散化模型与包含电容的离散化模型中的电流与电压值其中,/>为仿真步长。
4.根据权利要求2所述的基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法,其特征在于,还包括在电磁暂态仿真过程中,用等效电阻来表示开关管,等效电阻即为二值电阻,开关关断为大电阻,导通则为小电阻。
CN202410040594.2A 2024-01-11 2024-01-11 基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法 Active CN117556641B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410040594.2A CN117556641B (zh) 2024-01-11 2024-01-11 基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410040594.2A CN117556641B (zh) 2024-01-11 2024-01-11 基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法

Publications (2)

Publication Number Publication Date
CN117556641A true CN117556641A (zh) 2024-02-13
CN117556641B CN117556641B (zh) 2024-03-26

Family

ID=89823657

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410040594.2A Active CN117556641B (zh) 2024-01-11 2024-01-11 基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法

Country Status (1)

Country Link
CN (1) CN117556641B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106911132A (zh) * 2017-04-27 2017-06-30 武汉理工大学 基于mmc的分布式潮流控制器的电磁暂态数学模型、控制系统及建模方法
CN110472265A (zh) * 2019-06-24 2019-11-19 华北电力大学 一种isop型级联型电力电子变压器的电磁暂态等效建模方法
CN110489806A (zh) * 2019-07-23 2019-11-22 清华大学 包含多电压源型变流器电网的电磁暂态建模和计算方法
US20210050722A1 (en) * 2019-08-16 2021-02-18 Shanghai Jiao Tong University Hybrid electromagnetic transient simulation method for microgrid real-time simulation
CN112597731A (zh) * 2021-03-04 2021-04-02 清华四川能源互联网研究院 变流器电磁暂态仿真方法、装置、电子设备及存储介质
CN113190993A (zh) * 2021-04-27 2021-07-30 中国科学院电工研究所 高压大功率储能变流器电磁暂态快速仿真建模方法及系统
US20220269837A1 (en) * 2020-06-24 2022-08-25 Southeast University Converter parameterized constant admittance modeling method based on cross initialization
US20220366115A1 (en) * 2021-05-17 2022-11-17 North China Electric Power University General decoupling method and system for electromagnetic transient simulation of voltage source converter
CN117217151A (zh) * 2023-09-19 2023-12-12 华北电力大学 一种超级电容储能型mmc的多工况电磁暂态等效建模方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106911132A (zh) * 2017-04-27 2017-06-30 武汉理工大学 基于mmc的分布式潮流控制器的电磁暂态数学模型、控制系统及建模方法
CN110472265A (zh) * 2019-06-24 2019-11-19 华北电力大学 一种isop型级联型电力电子变压器的电磁暂态等效建模方法
CN110489806A (zh) * 2019-07-23 2019-11-22 清华大学 包含多电压源型变流器电网的电磁暂态建模和计算方法
US20210050722A1 (en) * 2019-08-16 2021-02-18 Shanghai Jiao Tong University Hybrid electromagnetic transient simulation method for microgrid real-time simulation
US20220269837A1 (en) * 2020-06-24 2022-08-25 Southeast University Converter parameterized constant admittance modeling method based on cross initialization
CN112597731A (zh) * 2021-03-04 2021-04-02 清华四川能源互联网研究院 变流器电磁暂态仿真方法、装置、电子设备及存储介质
CN113190993A (zh) * 2021-04-27 2021-07-30 中国科学院电工研究所 高压大功率储能变流器电磁暂态快速仿真建模方法及系统
US20220366115A1 (en) * 2021-05-17 2022-11-17 North China Electric Power University General decoupling method and system for electromagnetic transient simulation of voltage source converter
CN117217151A (zh) * 2023-09-19 2023-12-12 华北电力大学 一种超级电容储能型mmc的多工况电磁暂态等效建模方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐义良;赵成勇;赵禹辰;石璐;许建中;: "双端口子模块MMC电磁暂态通用等效建模方法", 中国电机工程学报, no. 20, 11 July 2018 (2018-07-11), pages 163 - 174 *

Also Published As

Publication number Publication date
CN117556641B (zh) 2024-03-26

Similar Documents

Publication Publication Date Title
CN106374767B (zh) 一种考虑二次回路的模块化多电平换流器子模块仿真模型
CN106649927B (zh) 一种基于fpga的电力电子元件实时仿真组合建模方法
CN103116665A (zh) 一种mmc拓扑变换器高效电磁暂态仿真方法
Liu et al. Feasible power-flow solution analysis of DC microgrids under droop control
CN103605850A (zh) 一种带子模块闭锁功能的mmc等效建模方法
CN106655195B (zh) 一种主动式配电网高频谐波潮流的计算方法
CN108629499A (zh) 一种基于二阶锥理论的配电网光伏电站选址定容方法
CN117556641B (zh) 基于Boost-全桥变流器解耦的恒定等效电导电磁暂态建模方法
CN110895644A (zh) 一种大规模变流器群并网的仿真系统及变流器的等效电路
CN105896613B (zh) 一种考虑通讯时滞的微电网分布式有限时间控制方法
CN110134004B (zh) 一种基于电力弹簧电路结构的pi控制参数整定方法
CN107221931A (zh) 基于有源电力滤波器的z源逆变器并网控制方法
CN111262466A (zh) 基于奇异摄动的模块化多电平变流器的建模控制方法
CN105743120A (zh) 单级式光伏z源逆变器变论域模糊控制mppt系统及控制方法
CN106844911A (zh) 一种大规模风电汇集系统静态电压稳定分析方法
CN113224788A (zh) 含有风电接入的交直流混联电网静态电压稳定评估方法
CN108649596B (zh) 一种适于负荷建模的电池储能系统动态模型
CN105739594B (zh) 单级式光伏z源逆变器模糊控制mppt方法
CN114021513B (zh) 单相变流器的非线性虚拟振荡器能量成型控制方法
CN112417667B (zh) 一种基于mmc高效电磁暂态桥臂等效模型的仿真方法
CN113890097B (zh) 一种变流器并网稳定性分析方法、装置、系统及存储介质
CN116760039B (zh) 节点类型扩展的交直流潮流计算方法、系统、设备及介质
CN115598994A (zh) 一种永磁直驱式风力发电系统的多速率实时仿真方法
CN112966463A (zh) 一种变流器群分散式并网系统高效仿真方法
CN107508286A (zh) 一种psasp软件潮流作业批量处理优化方法

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