CN107681690B - 光伏电站暂态等值电势在线辨识方法 - Google Patents

光伏电站暂态等值电势在线辨识方法 Download PDF

Info

Publication number
CN107681690B
CN107681690B CN201710916798.8A CN201710916798A CN107681690B CN 107681690 B CN107681690 B CN 107681690B CN 201710916798 A CN201710916798 A CN 201710916798A CN 107681690 B CN107681690 B CN 107681690B
Authority
CN
China
Prior art keywords
chopper circuit
trigger
photovoltaic power
voltage
power station
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
Application number
CN201710916798.8A
Other languages
English (en)
Other versions
CN107681690A (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.)
Shenyang University of Technology
State Grid Qinghai Electric Power Co Ltd
Electric Power Research Institute of State Grid Qinghai Electric Power Co Ltd
Original Assignee
Shenyang University of Technology
State Grid Qinghai Electric Power Co Ltd
Electric Power Research Institute of State Grid Qinghai Electric Power 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 Shenyang University of Technology, State Grid Qinghai Electric Power Co Ltd, Electric Power Research Institute of State Grid Qinghai Electric Power Co Ltd filed Critical Shenyang University of Technology
Priority to CN201710916798.8A priority Critical patent/CN107681690B/zh
Publication of CN107681690A publication Critical patent/CN107681690A/zh
Application granted granted Critical
Publication of CN107681690B publication Critical patent/CN107681690B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H02J3/383
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01Arrangements for reducing harmonics or ripples
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Inverter Devices (AREA)
  • Control Of Electrical Variables (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明提出一种光伏电站暂态等值电势在线辨识方法,将光伏发电站直流输出端与电网之间依次设置DC/DC斩波电路和逆变器,将逆变器输出端连接电压互感器,通过SPWM触发器触发DC/DC斩波电路;建立DC/DC斩波电路不同触发时刻电压值与触发时刻对应触发角的关系式,确定光伏电站SPWM的触发时间;触发DC/DC斩波电路,确定DC/DC斩波电路斩波电路的输出电压;计算逆变器输出的谐波电压;根据DC/DC斩波电路的输出电压、逆变器输出的谐波电压、测得的电压互感器输出电压以及光照强度变化系数计算光伏电站暂态等值电势;该方法为光伏电站暂态等值电势在线辨识提供了可靠依据。

Description

光伏电站暂态等值电势在线辨识方法
技术领域
本发明属于电网技术领域,具体涉及一种光伏电站暂态等值电势在线辨识方法。
背景技术
由于光伏电站极易受到外部环境影响,并且输出电势高低会与电网之间产生复杂交互影响,给光伏电站及接入电网的安全稳定运行带来严峻挑战,需要对光伏电站电势进行在线监测,计算出光伏电站的等值电势,为电力系统的稳定运行提供可靠保障。
光伏电站并网时对电的质量要求较高,由于传统的逆变器采用SPWM技术,难以对光伏电站暂态电势进行在线辨识,将正弦波等效成为一系列等时间间隔的脉冲波形,通过控制IGBT导通时间,形成需要的正弦波,并且由于光伏阵列的输出电压不稳定,利用DC/DC-DC/AC的谐波抑制方法可以有效的减少此种谐波分量,并且可以输出稳定的直流电压,使得电能的质量得以大幅度提升,此种方法为光伏电站暂态等值电势在线辨识提供了可靠依据。
发明内容
针对现有技术的不足,本发明提出一种光伏电站暂态等值电势在线辨识方法。
一种光伏电站暂态等值电势在线辨识方法,包括以下步骤:
步骤1:将光伏发电站直流输出端与电网之间依次设置DC/DC斩波电路和逆变器,将逆变器输出端连接电压互感器,通过SPWM触发器触发DC/DC斩波电路;
步骤2:获取光伏电站光伏阵列输出电压Ug,建立DC/DC斩波电路不同触发时刻电压值与触发时刻对应触发角的关系式,得到各触发时刻对应触发角,从而确定光伏电站SPWM的触发时间;
所述建立DC/DC斩波电路不同触发时刻电压值与触发时刻对应触发角的关系式如下所示:
DC/DC斩波电路第一触发时刻电压u1与第一触发时刻对应触发角α1的关系式如下所示:
Figure GDA0002593946330000011
其中,ut=Ug,α3为第三触发时刻对应触发角,α5为第五触发时刻对应触发角,α7为第七触发时刻对应触发角,α9为第九触发时刻对应触发角;
DC/DC斩波电路第三触发时刻电压u3与第三触发时刻对应触发角α3的关系式如下所示:
Figure GDA0002593946330000021
DC/DC斩波电路第五触发时刻电压u5与第五触发时刻对应触发角α5的关系式如下所示:
Figure GDA0002593946330000026
DC/DC斩波电路第七触发时刻电压u7与第七触发时刻对应触发角α7的关系式如下所示:
Figure GDA0002593946330000022
DC/DC斩波电路第九触发时刻电压u9与第九触发时刻对应触发角α9的关系式如下所示:
Figure GDA0002593946330000023
步骤3:在步骤2中确定的光伏电站SPWM的触发时间触发DC/DC斩波电路,根据DC/DC斩波电路的变压器一次侧线圈匝数、DC/DC斩波电路的变压器二次侧线圈匝数、以及各触发时刻对应触发角确定DC/DC斩波电路的输出电压Ud0
所述根据DC/DC斩波电路的变压器一次侧线圈匝数、DC/DC斩波电路的变压器二次侧线圈匝数、以及各触发时刻对应触发角确定DC/DC斩波电路的输出电压Ud0的计算公式如下所示:
Figure GDA0002593946330000024
其中,n1为DC/DC斩波电路的变压器一次侧线圈匝数,n2为DC/DC斩波电路的变压器二次侧线圈匝数。
步骤4:根据逆变器的频率调制比N、逆变器的幅值调制比M和光伏电站光伏阵列输出电压Ug计算逆变器输出的谐波电压Uxh1
所述根据逆变器的频率调制比N、逆变器的幅值调制比M和光伏电站光伏阵列输出电压Ug计算逆变器输出的谐波电压Uxh1的计算公式如下所示:
Figure GDA0002593946330000025
其中,m为谐波次数,
Figure GDA0002593946330000031
fc为载波频率,fm为调制波频率,Mm为载波幅值,Mc为调制波幅值,
Figure GDA0002593946330000032
为二次谐波初相位。
步骤5:根据DC/DC斩波电路的输出电压Ud0和逆变器输出的谐波电压Uxh1计算出光伏电站经过逆变器后的电势E1
所述根据DC/DC斩波电路的输出电压Ud0和逆变器输出的谐波电压Uxh1计算出光伏电站经过逆变器后的电势E1的计算公式如下所示:
Figure GDA0002593946330000033
步骤6:判断当前光照强度变化率是否超过设定的光照强度变化阈值,若是,则设定光照强度变化系数l=1,否则,设定光照强度变化系数l=0;
所述设定的光照强度变化阈值为100w/h。
步骤7:根据DC/DC斩波电路的输出电压Ud0、逆变器输出的谐波电压Uxh1、测得的电压互感器输出电压Uh、以及光照强度变化系数l计算光伏电站暂态等值电势E;
所述光伏电站暂态等值电势E的计算公式如下所示:
Figure GDA0002593946330000034
其中,n=6k,k=1,2,...,∞,ω为角频率。
本发明的有益效果:
本发明提出一种光伏电站暂态等值电势在线辨识方法,本发明将正弦波等效成为一系列等时间间隔的脉冲波形,通过控制IGBT导通时间,形成需要的正弦波,并且由于光伏阵列的输出电压不稳定,利用DC/DC-DC/AC的谐波抑制方法可以有效的减少此种谐波分量,并且可以输出稳定的直流电压,使得电能的质量得以大幅度提升,此种方法为光伏电站暂态等值电势在线辨识提供了可靠依据。
附图说明
图1为本实施方式中光伏电站暂态等值电势在线辨识方法的流程框图;
图2为本发明实施方式中光伏电站暂态等值电势在线辨识方法的装置结构框图。
具体实施方式
下面结合附图对本发明具体实施方式加以详细的说明。
本实施方式中,某光伏电站采用的是36V的电池板,其串并联数为10和100;故可算出其光伏电站光伏阵列输出电压Ug为360V,逆变器的幅值调制比M为0.8,逆变器的频率调制比N为40,变压器的变比为1:25,测得的电压互感器输出电压Uh=372V。
一种光伏电站暂态等值电势在线辨识方法,如图1所示,包括以下步骤:
步骤1:将光伏发电站直流输出端与电网之间依次设置DC/DC斩波电路和逆变器,将逆变器输出端连接电压互感器,通过SPWM触发器触发DC/DC斩波电路,如图2所示。
步骤2:获取光伏电站光伏阵列输出电压Ug,建立DC/DC斩波电路不同触发时刻电压值与触发时刻对应触发角的关系式,得到各触发时刻对应触发角,从而确定光伏电站SPWM的触发时间。
本实施方式中,建立DC/DC斩波电路不同触发时刻电压值与触发时刻对应触发角的关系式如下所示:
DC/DC斩波电路第一触发时刻电压u1与第一触发时刻对应触发角α1的关系式如式(1)所示:
Figure GDA0002593946330000041
其中,ut=Ug,α3为第三触发时刻对应触发角,α5为第五触发时刻对应触发角,α7为第七触发时刻对应触发角,α9为第九触发时刻对应触发角;
DC/DC斩波电路第三触发时刻电压u3与第三触发时刻对应触发角α3的关系式如式(2)所示:
Figure GDA0002593946330000042
DC/DC斩波电路第五触发时刻电压u5与第五触发时刻对应触发角α5的关系式如式(3)所示:
Figure GDA0002593946330000043
DC/DC斩波电路第七触发时刻电压u7与第七触发时刻对应触发角α7的关系式如式(4)所示:
Figure GDA0002593946330000051
DC/DC斩波电路第九触发时刻电压u9与第九触发时刻对应触发角α9的关系式如式(5)所示:
Figure GDA0002593946330000052
令DC/DC斩波电路第一触发时刻电压u1等于光伏电站光伏阵列输出电压Ug,令DC/DC斩波电路第三触发时刻电压u3、第五触发时刻电压u5、第七触发时刻电压u7、第九触发时刻电压u9等于0,得到公式如(6)-(10)所示:
Figure GDA0002593946330000053
Figure GDA0002593946330000054
Figure GDA0002593946330000055
Figure GDA0002593946330000056
Figure GDA0002593946330000057
通过上述公式(6)-(10)可得,
Figure GDA0002593946330000058
步骤3:在步骤2中确定的光伏电站SPWM的触发时间触发DC/DC斩波电路,根据DC/DC斩波电路的变压器一次侧线圈匝数、DC/DC斩波电路的变压器二次侧线圈匝数、以及各触发时刻对应触发角确定DC/DC斩波电路的输出电压Ud0
本实施方式中,DC/DC斩波电路的输出电压Ud0的计算公式如式(11)所示:
Figure GDA0002593946330000059
其中,n1为DC/DC斩波电路的变压器一次侧线圈匝数,n2为DC/DC斩波电路的变压器二次侧线圈匝数。
步骤4:根据逆变器的频率调制比N、逆变器的幅值调制比M和光伏电站光伏阵列输出电压Ug计算逆变器输出的谐波电压Uxh1
本实施方式中,逆变器输出的谐波电压Uxh1的计算公式如式(12)所示:
Figure GDA0002593946330000061
其中,m为谐波次数,
Figure GDA0002593946330000062
fc为载波频率,fm为调制波频率,Mm为载波幅值,Mc为调制波幅值,
Figure GDA0002593946330000063
为二次谐波初相位。
步骤5:根据DC/DC斩波电路的输出电压Ud0和逆变器输出的谐波电压Uxh1计算出光伏电站经过逆变器后的电势E1
本实施方式中,光伏电站经过逆变器后的电势E1的计算公式如式(13)所示:
Figure GDA0002593946330000064
步骤6:判断当前光照强度变化率是否超过设定的光照强度变化阈值,若是,则设定光照强度变化系数l=1,否则,设定光照强度变化系数l=0。
本实施方式中,设定的光照强度变化阈值为100w/h。
步骤7:根据DC/DC斩波电路的输出电压Ud0、逆变器输出的谐波电压Uxh1、测得的电压互感器输出电压Uh、以及光照强度变化系数l计算光伏电站暂态等值电势E。
本实施方式中,光伏电站暂态等值电势E的计算公式如式(14)所示:
Figure GDA0002593946330000065
其中,n=6k,k=1,2,...,∞,ω为角频率。

Claims (6)

1.一种光伏电站暂态等值电势在线辨识方法,其特征在于,包括以下步骤:
步骤1:将光伏发电站直流输出端与电网之间依次设置DC/DC斩波电路和逆变器,将逆变器输出端连接电压互感器,通过SPWM触发器触发DC/DC斩波电路;
步骤2:获取光伏电站光伏阵列输出电压Ug,建立DC/DC斩波电路不同触发时刻电压值与触发时刻对应触发角的关系式,得到第一、三、五、七、九触发时刻对应触发角,从而确定光伏电站SPWM的触发时间;
步骤3:在步骤2中确定的光伏电站SPWM的触发时间触发DC/DC斩波电路,根据DC/DC斩波电路的变压器一次侧线圈匝数、DC/DC斩波电路的变压器二次侧线圈匝数、以及各触发时刻对应触发角确定DC/DC斩波电路的输出电压Ud0
步骤4:根据逆变器的频率调制比N、逆变器的幅值调制比M和光伏电站光伏阵列输出电压Ug计算逆变器输出的谐波电压Uxh1
步骤5:根据DC/DC斩波电路的输出电压Ud0和逆变器输出的谐波电压Uxh1计算出光伏电站经过逆变器后的电势E1
步骤6:判断当前光照强度变化率是否超过设定的光照强度变化阈值,若是,则设定光照强度变化系数l=1,否则,设定光照强度变化系数l=0;
步骤7:根据DC/DC斩波电路的输出电压Ud0、逆变器输出的谐波电压Uxh1、测得的电压互感器输出电压Uh、以及光照强度变化系数l计算光伏电站暂态等值电势E;
所述光伏电站暂态等值电势E的计算公式如下所示:
Figure FDA0002593946320000011
其中,n=6k,k=1,2,...,∞,ω为角频率。
2.根据权利要求1所述的光伏电站暂态等值电势在线辨识方法,其特征在于,所述建立DC/DC斩波电路不同触发时刻电压值与触发时刻对应触发角的关系式如下所示:
DC/DC斩波电路第一触发时刻电压u1与第一触发时刻对应触发角α1的关系式如下所示:
Figure FDA0002593946320000012
其中,ut=Ug,α3为第三触发时刻对应触发角,α5为第五触发时刻对应触发角,α7为第七触发时刻对应触发角,α9为第九触发时刻对应触发角;
DC/DC斩波电路第三触发时刻电压u3与第三触发时刻对应触发角α3的关系式如下所示:
Figure FDA0002593946320000021
DC/DC斩波电路第五触发时刻电压u5与第五触发时刻对应触发角α5的关系式如下所示:
Figure FDA0002593946320000022
DC/DC斩波电路第七触发时刻电压u7与第七触发时刻对应触发角α7的关系式如下所示:
Figure FDA0002593946320000023
DC/DC斩波电路第九触发时刻电压u9与第九触发时刻对应触发角α9的关系式如下所示:
Figure FDA0002593946320000024
3.根据权利要求1所述的光伏电站暂态等值电势在线辨识方法,其特征在于,所述根据DC/DC斩波电路的变压器一次侧线圈匝数、DC/DC斩波电路的变压器二次侧线圈匝数、以及各触发时刻对应触发角确定DC/DC斩波电路的输出电压Ud0的计算公式如下所示:
Figure FDA0002593946320000025
其中,n1为DC/DC斩波电路的变压器一次侧线圈匝数,n2为DC/DC斩波电路的变压器二次侧线圈匝数。
4.根据权利要求1所述的光伏电站暂态等值电势在线辨识方法,其特征在于,所述根据逆变器的频率调制比N、逆变器的幅值调制比M和光伏电站光伏阵列输出电压Ug计算逆变器输出的谐波电压Uxh1的计算公式如下所示:
Figure FDA0002593946320000026
其中,m为谐波次数,
Figure FDA0002593946320000027
fc为载波频率,fm为调制波频率,Mm为载波幅值,Mc为调制波幅值,
Figure FDA0002593946320000031
为二次谐波初相位。
5.根据权利要求1所述的光伏电站暂态等值电势在线辨识方法,其特征在于,所述根据DC/DC斩波电路的输出电压Ud0和逆变器输出的谐波电压Uxh1计算出光伏电站经过逆变器后的电势E1的计算公式如下所示:
Figure FDA0002593946320000032
6.根据权利要求1所述的光伏电站暂态等值电势在线辨识方法,其特征在于,所述设定的光照强度变化阈值为100w/h。
CN201710916798.8A 2017-09-30 2017-09-30 光伏电站暂态等值电势在线辨识方法 Active CN107681690B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710916798.8A CN107681690B (zh) 2017-09-30 2017-09-30 光伏电站暂态等值电势在线辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710916798.8A CN107681690B (zh) 2017-09-30 2017-09-30 光伏电站暂态等值电势在线辨识方法

Publications (2)

Publication Number Publication Date
CN107681690A CN107681690A (zh) 2018-02-09
CN107681690B true CN107681690B (zh) 2020-11-13

Family

ID=61139603

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710916798.8A Active CN107681690B (zh) 2017-09-30 2017-09-30 光伏电站暂态等值电势在线辨识方法

Country Status (1)

Country Link
CN (1) CN107681690B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108711858B (zh) * 2018-07-19 2023-08-08 沈阳工业大学 基于跟随光伏电站出力波动的逆变器谐波抑制方法
CN109217375B (zh) * 2018-11-22 2020-05-22 国网陕西省电力公司电力科学研究院 一种基于adpss风力发电站的暂态等值电势辨识方法
CN109167395B (zh) * 2018-11-22 2021-10-15 国网宁夏电力有限公司电力科学研究院 基于adpss的光伏发电系统的暂态等值电势辨识方法
CN111193267A (zh) * 2019-12-03 2020-05-22 国网内蒙古东部电力有限公司经济技术研究院 一种基于云计算的特高压交直流电压稳定性控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001275258A (ja) * 2000-03-24 2001-10-05 Tokico Ltd 太陽電池システム
CN103762871A (zh) * 2013-12-17 2014-04-30 江苏博纬新能源科技有限公司 一种自举驱动单极性spwm调制非隔离并网逆变电路

Also Published As

Publication number Publication date
CN107681690A (zh) 2018-02-09

Similar Documents

Publication Publication Date Title
CN107681690B (zh) 光伏电站暂态等值电势在线辨识方法
Lorenzen et al. Control of a grid connected converter during weak grid conditions
CN109831093B (zh) 一种用于断路器测试中实现功率因数补偿的电流发生装置
TWI444807B (zh) 換流器之類比控制器
CN104124701B (zh) 一种光伏并网逆变器直流分量抑制方法
CN102611143B (zh) 一种三相并网逆变器并网电流的控制方法
CN107153156B (zh) 一种利用并联电容电流检测低压串联直流电弧的方法
CN106950512B (zh) 一种储能变流器并离网特性一体化检测系统及方法
CN106712527B (zh) 一种基于计算法的隔离式光伏阵列模拟器系统
CN103208939A (zh) 基于副边参考电流重构的光伏微逆变器、控制系统和方法
AU2015286473A1 (en) Power conversion device and three-phase AC power supply device
CN103746591A (zh) 一种h6单相非隔离光伏并网逆变器及其调制方法
CN104034984A (zh) 一种柔性直流输电运行试验中工程阀组件的短路试验方法
CN105141209A (zh) 减少隔离ups旁路导通的变压器励磁电流的方法及装置
CN109239570B (zh) 一种二极管正向电流浪涌实验电路
CN106849728A (zh) 带续流开关的箝位型三相非隔离光伏逆变器的控制方法
CN108462213B (zh) 基于守恒功率理论的多功能并网逆变器控制方法及系统
CN104753058B (zh) 一种直流母线电压谐波抑制计算方法
CN107465191B (zh) 光伏电站dc/dc-dc/ac谐波控制方法
Hussain et al. Grid integration of large capacity solar PV plant using multipulse VSC with robust PLL based control
CN106877716A (zh) 一种带续流开关的箝位型三相非隔离光伏逆变器
CN109842316A (zh) Dc-ac变换器的双前馈控制方法
CN105186888A (zh) 一种大功率电磁勘探恒流发射装置及其调制方法
CN104052321A (zh) 一种微逆变器过零点电流畸变抑制的混合控制方法及其系统
Kuzin et al. Power quality improvement in single-phase power supplies for resistance welding

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