CN107482593B - 一种基于能量频段的柔性直流线路电流差动保护方法 - Google Patents

一种基于能量频段的柔性直流线路电流差动保护方法 Download PDF

Info

Publication number
CN107482593B
CN107482593B CN201710557216.1A CN201710557216A CN107482593B CN 107482593 B CN107482593 B CN 107482593B CN 201710557216 A CN201710557216 A CN 201710557216A CN 107482593 B CN107482593 B CN 107482593B
Authority
CN
China
Prior art keywords
current
line
criterion
formula
energy
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
CN201710557216.1A
Other languages
English (en)
Other versions
CN107482593A (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.)
North China Electric Power University
Original Assignee
North China Electric Power 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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201710557216.1A priority Critical patent/CN107482593B/zh
Publication of CN107482593A publication Critical patent/CN107482593A/zh
Application granted granted Critical
Publication of CN107482593B publication Critical patent/CN107482593B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/261Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations
    • H02H7/262Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations involving transmissions of switching or blocking orders
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/261Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations
    • H02H7/263Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations involving transmissions of measured values
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/268Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured for dc systems

Landscapes

  • Emergency Protection Circuit Devices (AREA)

Abstract

本发明公开了一种基于能量频段的柔性直流线路电流差动保护方法,其特征在于,包括以下步骤:步骤一:在直流线路两端的断路器处各设置一个电流传感器;步骤二:分别在两个电流传感器处设置保护装置,保护装置之间通过光纤通信从而获取直流线路对端电流瞬时值,并计算差动电流瞬时值;步骤三:根据故障后电流变化量的幅值构成保护启动判据公式,当启动判据满足,则进入步骤四;步骤四:根据双极故障判据公式,判定系统是否发生双极短路故障,当判据满足,则进入步骤五;步骤五:根据区内故障判据公式,判定是否为本线路区内故障,当判据满足,则保护出口跳闸。

Description

一种基于能量频段的柔性直流线路电流差动保护方法
技术领域
本发明涉及电力系统继电保护技术领域,特别是涉及基于能量频段的柔性直流线路电流差动保护方法。
背景技术
随着电压源换流器(voltage sourced converter,VSC)的发展以及电源和负荷中电流变换中直流环节的增加,柔性直流系统由于其具有控制灵活、便于接纳分布式电源与直流负荷、供电容量大、电能质量好等优点已引起国内外广泛关注。然而缺乏成熟的直流线路保护技术是柔性直流系统面临的主要挑战之一。柔性直流系统线路保护的难点在于既要在短时间内快速识别直流故障以防止设备损坏,又要准确选择故障线路以减少停电范围。
目前柔性直流线路电流差动保护方法包括:
1.单端量保护难以实现清晰的保护边界,不能较好的满足柔性直流线路保护快速性与选择性需求。
2.方向纵联保护主要基于电流突变量判定电流方向,然而故障后VSC快速闭锁,系统拓扑发生根本变化,叠加原理将不再适用。
3.基于电流突变量的方向纵联保护从原理上是不适用于柔性直流系统的。
4.电流差动保护原理简单可靠,然而直流线路电流差动保护易受线路分布电容暂态电流的影响。如果电流差动保护的定值过低,则区外故障时线路分布电容的暂态电流易导致保护误动,保护的可靠性难以满足;如果定值过高,则区内高阻故障容易拒动,保护的灵敏性难以满足。
因此希望有一种柔性直流线路电流差动保护方法可以克服或至少减轻现有技术中直流差动保护易受分布电容暂态电流影响的问题。
发明内容
本发明的目的在于提供一种基于能量频段的柔性直流线路电流差动保护方法,利用区内故障时差流能量主要集中于低频段、区外故障时集中于高频段的特点,识别区内故障。
本发明提供一种基于能量频段的柔性直流线路电流差动保护方法,包括以下步骤:
步骤一:在所述直流线路两端的断路器处各设置一个电流传感器;
步骤二:分别在两个电流传感器处设置保护装置,保护装置之间通过光纤通信从而获取所述直流线路对端电流瞬时值,并计算差动电流瞬时值;
步骤三:根据故障后电流变化量的幅值构成保护启动判据公式,当启动判据满足,则进入步骤四;
步骤四:根据双极故障判据公式,判定系统是否发生双极短路故障,当判据满足,则进入步骤五;
步骤五:根据区内故障判据公式,判定是否为本线路区内故障,当判据满足,则保护出口跳闸。
优选地,所述步骤一中的电流传感器的测量频段要求覆盖0-50kHz全频段,所述电流传感器的安装极性为电流从母线指向所述直流线路为正。
优选地,所述步骤二中的保护装置为继电保护器。
优选地,所述步骤二中的差动电流瞬时值的计算公式为:
icd=iM+iN
式中iM和iN分别是所述直流线路两端的电流瞬时值。
优选地,所述步骤三中的启动判据公式为:
式中为正极或负极电流变化量,通过所述直流线路两端的电流瞬时值减去其1ms之前的值计算获得;In为所述直流线路的电流额定值。
优选地,所述步骤四中区内故障判据公式为:
式中udc为极间电压的瞬时值;Un为极间电压额定值;ip和in分别是本端的正极电流和负极电流的瞬时值。
优选地,所述步骤五中利用差动电流的1.56kHz以下的低频能量与1.56kHz以上的高频能量的比值识别区内故障,区内故障判据公式为:
EL/EH>1
式中EL表示差动电流经小波变换分解重构出的1.56kHz以下的低频信号的能量,EH表示差动电流经小波变换分解重构出的1.56kHz以上的高频信号的能量,并且定义时域信号平方的时间积分为该时域信号的能量。
本发明公开了一种基于能量频段的柔性直流线路电流差动保护方法,本方法将分布电容的影响变为保护的制动量,因此区外故障发生时分布电容电流不会造成保护误动,而且无需定值整定。
附图说明
图1是二进正交小波对信号分解示意图;
图2是基于能量频段的柔性直流线路电流差动保护方法的工作流程图;
图3是基于能量频段的柔性直流线路电流差动保护系统的保护测量点示意图;
图4是区内故障时差动电流波形图;
图5是区外故障时差动电流波形图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。下面结合附图对本发明的实施例进行详细说明。
以两电平VSC作为换流器的±10kV柔性直流配电网系统为实施例,说明本发明。
首先分析电源和换流器产生的故障电流的频段分布特征。
双极短路故障的发展过程可分为三个阶段:电容放电阶段、二极管自由导通阶段和交流侧馈入阶段。
在电容放电阶段,等值电路为二阶振荡电路,故障电流的主要成分为衰减振荡分量。根据电路分析可知,此阶段的故障电流振荡的频率为:
式中x为故障距离,C为直流并联电容器的容值,Lu和Ru分别是线路单位长度的电感值和电阻值。
从上式可知,电流频率与故障距离和并联电容器的容值成反比。由于并联电容器的容值一般较大,因此故障电流的振荡频率一般较低。对于本文的系统参数,0.5km处发生双极短路故障,电流的振荡频率约为875Hz。当故障距离更大时,电流的振荡频率将更低。对于故障距离小于0.5km的故障,电容放电阶段存在的时间极短,对一个时间窗内的故障电流能量分布影响不大。因此,在电容放电阶段,故障电流的能量集中在数百Hz的频段。
在二极管自由导通阶段,等值电路为一阶衰减电路。故障点电流的表达式为:
式中I0为本阶段的初始故障电流。
由于故障电流呈指数衰减,按正常线路参数估算,二极管自由导通阶段故障电流的能量集中在更低的频段。
在交流馈入阶段,交流侧通过二极管整流桥(闭锁后的换流器)向故障点馈入电流,根据二极管整流桥的特性,故障点电流主要包含直流分量和六次谐波分量,其能量主要集中于300Hz(六次谐波的频率)以下。
综上,在其整个故障过程中,故障电流的能量主要集中于1000Hz以下的低频段。
再分析分布电容暂态电流的能量频段分布特征。
在故障暂态阶段,故障行波的传播导致了电压的变化,进而引起了线路分布电容的充放电,产生了高频的分布电容暂态电流。因此分布电容暂态电流的频率与行波的固有频率一致。行波的固有频率与故障距离和边界条件有关,可表示为如下数学关系:
式中v为行波波速,d为故障距离;θS和θF分别是线路始端和故障点处的反射系数角。线路始端的反射系数角θS的取值范围为0到π;无论是金属性故障还是经过渡电阻故障,故障点处的反射系数角θF都为0。
因此行波固有频率的理论最低值为:
考虑到±10kV直流配电网的送电距离一般在40km以内,认为故障距离d不会大于40km。对于行波高频分量,其波速v约为光速的97%~99%。可计算出固有频率主频不会低于1800Hz。因此,线路分布电容暂态电流的能量集中于1800Hz以上的频段。
综上,对于区内双极故障,线路差动电流的能量分布主要由故障电流决定,其能量集中于1000Hz以下的低频段。对于区外故障,线路差动电流只含有分布电容暂态电流,其能量集中于1800Hz以上的高频段。
如图1所示,由于小波分析在时域和频域上同时具有良好的局部化性质,能对不同的频率成分采用逐渐精细的采样步长,其对检测高频和低频信号同时有效,因此本发明选择离散小波变换(discrete wavelet transform,DWT)对差流信号进行分析。
小波变换可按一定的尺度将不同频率的信号划分到不同的频段中,然后对各个子频段进行重构,从而分离出各个频段的时域信号。对于采样频率为50kHz的信号,利用小波变换对信号进行4层分解重构,每层的频带范围如下:第1层(d1):12.5kHz~25kHz;第2层(d2):6.25kHz~12.5kHz;第3层(d3):3.125kHz~6.25kHz;第4层(d4):1.56kHz~3.125kHz;第4层(a4):0kHz~1.56kHz。因此,利用小波变换可实现对差动电流信号的能量频段分析。
如图2和3所示,基于能量频段的柔性直流线路电流差动保护方法的工作流程包括如下步骤:
步骤11:在直流线路两端各装设一个电流传感器,通过光纤通信,保护装置可计算出差动电流瞬时值。电流传感器的安装极性为电流从母线指向线路为正。差动电流的计算公式为:
icd=iM+iN
式中iM和iN分别是线路两端电流的瞬时值。
步骤12:根据电流变化量的幅值构成保护启动判据,公式为:
式中为正极或负极电流变化量,通过电流瞬时值减去其1ms之前的值计算获得;In为线路的电流额定值。
当启动判据满足,则进入下一步骤。
步骤13:根据双极故障判据,判定是否发生双极短路故障,区内故障判据公式为:
式中udc为极间电压的瞬时值;Un为极间电压额定值;ip和in分别是本端的正极电流和负极电流的瞬时值。
当判据满足,则进入下一步骤
步骤14:提取故障时刻后2ms的差动电流数据,进行小波分解重构,计算出低频能量与高频能量。定义1.56kHz以下为低频段,1.56kHz以上为高频段,时域信号平方的时间积分为该信号的能量。
步骤15:根据区内故障判据,判定是否为本线路区内故障,双极故障判据公式为:
EL/EH>1
式中EL表示小波变换分解重构出的低频信号(1.56kHz以下)的能量,EH表示分解重构出的高频信号(1.56kHz以上)的能量。
当判据满足,则保护出口跳闸。
下面给出本发明在PSCAD/EMTDC上搭建的±10kV柔性直流系统的测试结果。
如图4所示,对于区内故障,图4(a)为差动电流的原始时域信号;图4(b)为通过离散小波变换对差流信号进行分解、重构后获得的低频时域信号,其频段为0-1.56kHz;图4(c)为分解、重构后获得的差流高频时域信号,其频段为1.56-25kHz。经计算,故障后2ms内差动电流低频信号的能量为48080.549A2·s,高频信号的能量为124.336A2·s,两者的比值为386.698,远大于1,满足区内故障判据,保护正确动作。
如图5所示,对于区外故障,图5(a)为差动电流的原始时域信号;图5(b)为通过离散小波变换对差流信号进行分解、重构后获得的低频时域信号,其频段为0-1.56kHz;图5(c)为分解、重构后获得的差流高频时域信号,其频段为1.56-25kHz。经计算,故障后2ms内差动电流低频信号的能量为1.421A2·s,高频信号的能量为64.250A2·s,两者的比值为0.022,远小于1,不满足区内故障判据,保护不会误动。
最后需要指出的是:以上实施例仅用以说明本发明的技术方案,而非对其限制。尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种基于能量频段的柔性直流线路电流差动保护方法,其特征在于,包括以下步骤:
步骤一:在所述直流线路两端的断路器处各设置一个电流传感器;
步骤二:分别在两个电流传感器处设置保护装置,保护装置之间通过光纤通信从而获取所述直流线路对端电流瞬时值,并计算差动电流瞬时值;
步骤三:根据故障后电流变化量的幅值构成保护启动判据公式,当启动判据满足,则进入步骤四;
步骤四:根据双极故障判据公式,判定系统是否发生双极短路故障,当判据满足,则进入步骤五;
步骤五:利用差动电流的1.56kHz以下的低频能量与1.56kHz以上的高频能量的比值识别区内故障,根据区内故障判据公式,判定是否为本线路区内故障,当判据满足,则保护出口跳闸。
2.如权利要求1所述的基于能量频段的柔性直流线路电流差动保护方法,其特征在于:所述步骤一中的电流传感器的测量频段要求覆盖0-50kHz全频段,所述电流传感器的安装极性为电流从母线指向所述直流线路为正。
3.如权利要求1所述的基于能量频段的柔性直流线路电流差动保护方法,其特征在于:所述步骤二中的保护装置为继电保护器。
4.如权利要求1所述的基于能量频段的柔性直流线路电流差动保护方法,其特征在于:所述步骤二中的差动电流瞬时值的计算公式为:
icd=iM+iN
式中iM和iN分别是所述直流线路两端的电流瞬时值。
5.如权利要求4所述的基于能量频段的柔性直流线路电流差动保护方法,其特征在于:所述步骤三中的启动判据公式为:
式中为正极或负极电流变化量,通过所述直流线路两端的电流瞬时值减去其1ms之前的值计算获得;In为所述直流线路的电流额定值。
6.如权利要求1所述的基于能量频段的柔性直流线路电流差动保护方法,其特征在于:所述步骤四中区内故障判据公式为:
式中udc为极间电压的瞬时值;Un为极间电压额定值;ip和in分别是本端的正极电流和负极电流的瞬时值。
7.如权利要求1所述的基于能量频段的柔性直流线路电流差动保护方法,其特征在于:所述步骤五中区内故障判据公式为:
EL/EH>1
式中EL表示差动电流经小波变换分解重构出的1.56kHz以下的低频信号的能量,EH表示差动电流经小波变换分解重构出的1.56kHz以上的高频信号的能量,并且定义时域信号平方的时间积分为该时域信号的能量。
CN201710557216.1A 2017-07-10 2017-07-10 一种基于能量频段的柔性直流线路电流差动保护方法 Active CN107482593B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710557216.1A CN107482593B (zh) 2017-07-10 2017-07-10 一种基于能量频段的柔性直流线路电流差动保护方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710557216.1A CN107482593B (zh) 2017-07-10 2017-07-10 一种基于能量频段的柔性直流线路电流差动保护方法

Publications (2)

Publication Number Publication Date
CN107482593A CN107482593A (zh) 2017-12-15
CN107482593B true CN107482593B (zh) 2019-05-14

Family

ID=60596259

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710557216.1A Active CN107482593B (zh) 2017-07-10 2017-07-10 一种基于能量频段的柔性直流线路电流差动保护方法

Country Status (1)

Country Link
CN (1) CN107482593B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108306263B (zh) * 2018-01-15 2019-09-20 清华大学 电力线路行波差动保护方法、装置、设备及介质
CN108258665B (zh) * 2018-01-15 2019-08-09 清华大学 并联电抗补偿线路行波差动保护方法、装置、设备及介质
CN108270207B (zh) * 2018-01-15 2019-10-18 清华大学 串联电容补偿线路行波差动保护方法、装置、设备及介质
CN109638789B (zh) * 2018-10-29 2020-07-07 全球能源互联网研究院有限公司 柔性交直流互联系统的直流双极短路故障保护方法及装置
CN110364997B (zh) * 2019-07-03 2022-09-20 中国电力科学研究院有限公司 一种基于单端电流差分量的直流输电线路后备保护方法及系统
CN110601153B (zh) * 2019-08-19 2020-12-25 西安交通大学 直流配电网单端量的继电保护方法
CN110824294B (zh) * 2019-10-18 2022-02-18 济南大学 一种柔性直流配电网双极短路故障区段识别及测距方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2306607A2 (en) * 2009-08-31 2011-04-06 EMforce B.V. Device and method to protect an electric power distribution network against current faults
CN102005739A (zh) * 2010-08-27 2011-04-06 昆明理工大学 一种识别特高压直流输电线路区内外故障小波能量的方法
CN102570419A (zh) * 2011-12-28 2012-07-11 西安交通大学 一种基于电流量的输电线路纵联保护方法
CN103633629A (zh) * 2013-11-11 2014-03-12 同济大学 基于小波变换和能量谱分析的高压直流输电线路的保护方法
CN104820159A (zh) * 2015-04-30 2015-08-05 国家电网公司 一种柔性直流输电系统直流双极短路故障辨识方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2306607A2 (en) * 2009-08-31 2011-04-06 EMforce B.V. Device and method to protect an electric power distribution network against current faults
CN102005739A (zh) * 2010-08-27 2011-04-06 昆明理工大学 一种识别特高压直流输电线路区内外故障小波能量的方法
CN102570419A (zh) * 2011-12-28 2012-07-11 西安交通大学 一种基于电流量的输电线路纵联保护方法
CN103633629A (zh) * 2013-11-11 2014-03-12 同济大学 基于小波变换和能量谱分析的高压直流输电线路的保护方法
CN104820159A (zh) * 2015-04-30 2015-08-05 国家电网公司 一种柔性直流输电系统直流双极短路故障辨识方法

Also Published As

Publication number Publication date
CN107482593A (zh) 2017-12-15

Similar Documents

Publication Publication Date Title
CN107482593B (zh) 一种基于能量频段的柔性直流线路电流差动保护方法
CN103529332B (zh) 一种基于电压相关度和小波变换暂态能量分布特性的特高压直流输电线路雷击干扰识别方法
CN108054736B (zh) 一种基于电压极波的多端柔性直流电网直流线路保护方法
CN110376490A (zh) 基于同步挤压小波变换的高压直流输电线路故障测距方法
CN107064741A (zh) 一种配电网线路异名相两点相继接地故障选线方法
Zand et al. Fault locating transmission lines with thyristor-controlled series capacitors By fuzzy logic method
CN108896869A (zh) 一种基于s变换的高压直流输电线路雷击干扰识别方法
CN101915888B (zh) ±800kV直流输电线路雷击干扰的可拓融合识别方法
Biswas et al. A dual-time transform assisted intelligent relaying scheme for the STATCOM-compensated transmission line connecting wind farm
CN110231552A (zh) 煤矿电缆绝缘监测选线方法、信号注入电路及相应装置
CN108899879A (zh) 基于零序电流投影系数的小电阻接地系统接地保护方法
CN108258662A (zh) 一种多端柔性直流配电线路暂态保护方法及装置
CN105896485A (zh) 基于新能源并网的特高压长距离直流线路的保护算法
CN106611957A (zh) 电力系统次/超同步振荡广域监测方法和系统
CN102279309A (zh) 一种特高压直流线路保护启动判据的方法
CN109270333B (zh) 一种识别雷击引起的电压暂降的方法
CN108270239A (zh) 一种含分布式电源的配电网电能质量扰动源方向判定方法
Hossam-Eldin et al. Combined traveling wave and fuzzy logic based fault location in multi-terminal HVDC systems
CN104638671A (zh) 基于引入粒子群算法及相位扰动的孤岛检测方法
CN110350493A (zh) 基于线路电流二阶导数的中压柔性直流系统故障检测方法
CN107315127B (zh) 高阻接地故障线路检测方法和系统
Saber A backup protection algorithm for bipolar line-commutated converter HVDC lines
Tajani et al. A novel differential protection scheme for AC microgrids based on discrete wavelet transform
CN102082420B (zh) 输电线路的纵差保护方法
CN109387713A (zh) 一种分布式并网孤岛检测的混合方法

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