CN112083286B - 一种风电场集电线路单相接地故障选线方法 - Google Patents

一种风电场集电线路单相接地故障选线方法 Download PDF

Info

Publication number
CN112083286B
CN112083286B CN202010926138.XA CN202010926138A CN112083286B CN 112083286 B CN112083286 B CN 112083286B CN 202010926138 A CN202010926138 A CN 202010926138A CN 112083286 B CN112083286 B CN 112083286B
Authority
CN
China
Prior art keywords
current
zero sequence
matrix
sequence current
correlation coefficient
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
CN202010926138.XA
Other languages
English (en)
Other versions
CN112083286A (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.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
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 Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202010926138.XA priority Critical patent/CN112083286B/zh
Publication of CN112083286A publication Critical patent/CN112083286A/zh
Application granted granted Critical
Publication of CN112083286B publication Critical patent/CN112083286B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/085Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/52Testing for short-circuits, leakage current or ground faults
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/58Testing of lines, cables or conductors

Abstract

本发明涉及一种风电场集电线路单相接地故障选线方法,属于电力系统继电保护技术领域。若判断系统发生单相接地故障,对各集电线路的零序电流进行采样,选取暂态零序电流进行2层小波分解。提取低频带的零序电流小波系数,并计算各集电线路之间的Hausdorff双向距离,形成各集电线路的Hausdorff距离矩阵。对Hausdorff距离矩阵每列进行求和,形成各集电线路零序电流综合关联系数矩阵。选取零序电流综合关联系数最大的线路,计算该集电线路零序电流综合关联系数商,形成该集电线路零序电流综合关联系数商矩阵。根据该矩阵的互关联系数是否大于阈值,判断该集电线路是否发生单相接地故障。理论分析、仿真及实际录播数据验证本发明效果良好。

Description

一种风电场集电线路单相接地故障选线方法
技术领域
本发明涉及一种风电场集电线路单相接地故障选线方法,属于电力系统继电保护技术领域。
背景技术
从本质上讲,风电场集电系统就是一个有源配电系统。早期风电场集电线路系统多沿用小电流接地方式,但区别于配电网保护的侧重点,近年来,越来越多的风电场将小电流接地方式更换为小电阻接地方式。相较于小电阻接地方式,小电流接地系统发生单相接地故障时,接地电流更加微弱,多数情况下电气量临界于甚至未达到保护整定值,仅依靠继电保护无法保证故障可靠切除。因此,许多学者研究了多种单相接地故障选线及定位方法,对继电保护形成辅助的故障辨识方式,保证故障以可靠切除。
而对于小电阻接地系统,低阻故障时接地电流高达上百或上千安培,基于零序分量的大多数继电保护都有足够的选择性和灵敏性,正确动作,无需进行选线。然而,常规风电场集电线路的保护整定值较高且为定值,在系统发生高阻故障时,系统保护越级动作、或不可靠动作。因此,利用集中保护装置,对集电线路信息集中比较分析,甄选出故障线路,配合常规保护对其进行切除,是一个行之有效的方法。
发明内容
本发明要解决的技术问题是针对风电场小电阻接地系统集电线路高阻故障时零序保护不可靠的问题,提出一种风电场集电线路单相接地故障选线方法,用以解决上述问题。
本发明的技术方案是:一种风电场集电线路单相接地故障选线方法,利用零序电压突变量判断风电场小电阻接地集电系统是否异常运行,若判断系统发生单相接地故障,对各集电线路的零序电流进行采样,选取暂态零序电流进行2层小波分解,提取低频带的零序电流小波系数,并计算各集电线路之间的Hausdorff双向距离,形成各线路的Hausdorff距离矩阵,对Hausdorff距离矩阵每列进行求和,形成各集电线路零序电流综合关联系数矩阵,选取零序电流综合关联系数最大的线路,计算该集电线路零序电流综合关联系数商,形成该集电线路零序电流综合关联系数商矩阵,若该矩阵的互关联系数均大于阈值,则该条集电线路发生单相接地故障,若该矩阵的互关联系数均小于阈值,则汇集母线发生单相接地故障。
具体步骤为:
Step1:提取风电场集电系统中性点零序电压,定义突变量启动判据为:
Figure GDA0003210236530000011
式(1)中,k表示当前采样点,K表示一定时窗内的采样点数,k≥K,Δu(k)=u0(k)-u0(k-1), u0(k)为零序电压采样值,当零序电压突变量超过设定的阈值时,判断风电集电系统发生了单相接地故障。
Step2:提取各集电线路出口保护处的故障后5ms(T/4工频周期)的零序电流数据,对各集电线路零序电流进行小波2层分解与重构,提取低频带零序电流小波系数,根据式(2) 求取第i和j条集电线路低频带零序电流小波系数的Hausdorff距离,形成各线路的Hausdorff距离矩阵H。
H(i,j)=max[h(i0i,i0j),h(i0j,i0i)] (2)
式(2)中,h(i0i,i0j) 、 h(i0j,i0i) 分别为第i到第j条集电线路、第j到第i条集电线路零序电流小波系数的Hausdorff单向距离。
Figure GDA0003210236530000022
Step3:通过式(4)求取各集电线路零序电流波形综合关联系数Hj,并根据式(5)形成零序电流波形综合关联系数矩阵HΣ
Figure GDA0003210236530000023
HΣ=[H1 … Hj … Hk … Hm],j=1,…,j,…,k,…m (5)
Step4:假设第k条集电线路对应的元素Hk为矩阵HΣ的最大值,根据式(6)求取各集电线路综合关联系数商Skj,并根据式(7)综合关联系数商矩阵Sk
Figure GDA0003210236530000024
Sk=[Sk1 … Skj … 1 … Skm],j=1,…,j,…,k,…m (7)
Step5:设定阈值Sth,若矩阵Sk互关联系数均大于Sth,则第k条集电线路故障;若矩阵Sk互关联系数均小于Sth,则为汇集母线故障。经大量仿真分析,将阈值设为1.5即可对集电线路故障或汇集线路故障有效区分。
本发明的有益效果是:
1、本发明利用零序电压突变量作为故障选线启动,极大提高了选线装置的灵敏性。
2、本发明采用故障后T/4时窗的零序电流,此时电流互感器采集到的零序电流波形还未发生畸变,避免了波形畸变带来的影响。
3、本发明使用的选线算法相较相关系数分析法,能更好地反映各集电线路零序电流波形的相似度,性能较相关系数法占优。
附图说明
图1是本发明实施例中风电场集电系统仿真模型图;
图2是本发明实施例1中暂态零序电流小波系数波形图;
图3是本发明实施例2中暂态零序电流小波系数波形图。
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步说明。
一种风电场集电线路单相接地故障选线方法,利用零序电压突变量判断风电场小电阻接地集电系统是否异常运行,若判断系统发生单相接地故障,对各集电线路的零序电流进行采样,选取暂态零序电流进行2层小波分解,提取低频带的零序电流小波系数,并计算各集电线路之间的Hausdorff双向距离,形成各集电线路的Hausdorff距离矩阵,对Hausdorff距离矩阵每列进行求和,形成各集电线路零序电流综合关联系数矩阵,选取零序电流综合关联系数最大的线路,计算该集电线路零序电流综合关联系数商,形成该集电线路零序电流综合关联系数商矩阵,若该矩阵的互关联系数均大于阈值,则该集电线路发生单相接地故障,若该矩阵的互关联系数均小于阈值,则汇集母线发生单相接地故障。
具体步骤为:
Step1:提取风电场集电系统中性点零序电压,定义突变量启动判据为:
Figure GDA0003210236530000031
式(1)中,k表示当前采样点,K表示一定时窗内的采样点数,k≥K,Δu(k)=u0(k)-u0(k-1), u0(k)为零序电压采样值,当零序电压突变量超过设定的阈值时,判断风电集电系统发生了单相接地故障。
Step2:提取各集电线路出口保护处的故障后5ms(T/4工频周期)的零序电流数据,对各集电线路零序电流进行小波2层分解与重构,提取低频带零序电流小波系数,根据式(2) 求取第i和j条集电线路低频带零序电流小波系数的Hausdorff距离,形成各集电线路的Hausdorff距离矩阵H。
H(i,j)=max[h(i0i,i0j),h(i0j,i0i)] (2)
式(2)中,h(i0i,i0j) 、 h(i0j,i0i) 分别为第i到第j条集电线路、第j到第i条集电线路零序电流小波系数的Hausdorff单向距离。
Figure GDA0003210236530000042
Step3:通过式(4)求取各集电线路零序电流波形综合关联系数Hj,并根据式(5)形成零序电流波形综合关联系数矩阵HΣ
Figure GDA0003210236530000043
HΣ=[H1 … Hj … Hk … Hm],j=1,…,j,…,k,…m (5)
Step4:假设第k条线路对应的元素Hk为矩阵HΣ的最大值,根据式(6)求取各集电线路综合关联系数商Skj,并根据式(7)综合关联系数商矩阵Sk
Figure GDA0003210236530000044
Sk=[Sk1 … Skj … 1 … Skm],j=1,…,j,…,k,…m (7)
Step5:设定阈值Sth,若矩阵Sk互关联系数均大于Sth,则第k条集电线路故障;若矩阵Sk互关联系数均小于Sth,则为汇集母线故障。
实施例1:如图1所示,构建典型风电场集电系统仿真模型,模型具体参数如下:
(1)35kV汇集母线上有4条集电线路L1~L4,其中L1~L3每条集电线路有3台双馈风电机组,L4含有分支,共有8台双馈风电机组;共接入17台双馈风机。双馈风电机组经由0.69kV/38.5kV的Y/△型箱变升压连接到集电线路。每台风机容量为2MW,风电机组总容量为17*2MW=34MW。
(2)集电线路L1~L4主干线路总长分别为:24km、17km、14km、23km,其中L1~L4上游为电缆线路,分别长20km、13km、10km、11km,下游风机之间为架空线路,每两台风电机组间距为2km。集电线路L4分支1、分支2长度分别为8km、4km。
(3)35kV母线经由38.5kV/220kV主变升压至220kV,主变容量为100MVA,Ynd11连接。主变低压侧通过Z型变压器经50Ω电阻接地,主变高压侧经送出线路与220kV大电网系统相连。
设置距离汇集母线13km,集电线路L4的分支1(距分支1出口2km)上发生A相接地故障,过渡电阻为1000Ω,故障时刻为0.23s,持续时长为0.05s。L1~L4集电线路出口0.01km分别有W1、W2、W3、W4个检测点,采样率为10kHz。本例的风电场集电线路单相接地故障选线方法的具体步骤如下:
按照步骤(1),采集中性点零序电压和各集电线路零序电流,计算零序电压突变量E(k)=203897V2,零序电压突变量大于设定的阈值,判定风电场集电系统发生单相接地故障。
确定集电系统发生单相接地故障后,按照步骤(2),选取故障后5ms各馈线的零序电流进行db4小波2层分解,暂态零序电流0~2.5kHz频带的小波系数波形如图2所示,采用暂态零序电流0~2.5kHz频带的小波系数计算各线路之间的Hausdorff双向距离,集电线路L1~L4的Hausdorff距离矩阵H及综合关联系数矩阵HΣ
Figure GDA0003210236530000051
HΣ=[7.4406 7.4406 7.7562 20.5521]
由上述矩阵可看出,第4条集电线路的零序电流综合关联系数最大,计算各集电线路的综合关联系数商S4j,并形成综合关联系数商矩阵S4
S4=[2.7622 2.7622 2.64981]
由上述矩阵可知,S4的互关联系数均大于Sth,根据步骤(4)可知,可判断集电线路L4为故障集电线路;L1~L3为健全集电线路。
实施例2:2018年2月9日,云南省曲靖市支锅山风电场35kV母线上某集电线路发生单相接地故障。该风电场中性点接地方式为小电阻接地,35kV汇集母线I段上共有5条集电线路。本例的风电场集电线路单相接地故障选线方法的具体步骤如下:
按照步骤(1),采集中性点零序电压和各集电线路零序电流,计算零序电压突变量E(k)=14608V2,零序电压突变量大于设定的阈值,判定风电场集电系统发生单相接地故障。
确定集电系统发生单相接地故障后,按照步骤(2),选取故障后5ms各集电线路的零序电流进行db4小波2层分解,暂态零序电流低频带的小波系数波形如图3所示,采用暂态零序电流低频带的小波系数计算各线路之间的Hausdorff双向距离,L1~L5的Hausdorff距离矩阵H及综合关联系数矩阵HΣ
Figure GDA0003210236530000061
HΣ=[3.2046 3.4103 3.2755 10.2535 4.2325]
由上述矩阵可看出,第4条集电线路的零序电流综合关联系数最大,计算各集电线路的综合关联系数商S4j,并形成综合关联系数商矩阵S4
S4=[3.1996 3.0066 3.1304 1 2.4225]
由上述矩阵可知,S4互关联系数均大于Sth,根据步骤(4)可知,可判断集电线路L4为故障集电线路;L1~L3、L5为健集电全线路。与实际故障集电线路相符。
以上结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (1)

1.一种风电场集电线路单相接地故障选线方法,其特征在于:利用零序电压突变量判断风电场小电阻接地集电系统是否异常运行,若判断系统发生单相接地故障,对各集电线路的零序电流进行采样,选取暂态零序电流进行2层小波分解,提取低频带的零序电流小波系数,并计算各集电线路之间的Hausdorff双向距离,形成各集电线路的Hausdorff距离矩阵,对Hausdorff距离矩阵每列进行求和,形成各条集电线路的零序电流综合关联系数矩阵,选取零序电流综合关联系数最大的集电线路,计算该条集电线路的零序电流综合关联系数商,形成该条集电线路的零序电流综合关联系数商矩阵,若该矩阵的互关联系数均大于阈值,则该条集电线路发生单相接地故障,若该矩阵的互关联系数均小于阈值,则汇集母线发生单相接地故障;
具体步骤为:
Step1:提取风电场集电系统中性点零序电压,定义突变量启动判据为:
Figure FDA0003210236520000011
式(1)中,k表示当前采样点,K表示一定时窗内的采样点数,k≥K,Δu(n)=u0(n)-u0(n-1),u0(n)为零序电压采样值,当零序电压突变量超过设定的阈值时,判断风电集电系统发生了单相接地故障;
Step2:提取各集电线路出口保护处的故障后5ms的零序电流数据,对各集电线路零序电流进行小波2层分解与重构,提取低频带零序电流小波系数,根据式(2)求取第i和j条集电线路低频带零序电流小波系数的Hausdorff距离,形成各集电 线路的Hausdorff距离矩阵H;
H(i,j)=max[h(i0i,i0j),h(i0j,i0i)] (2)
式(2)中,h(i0i,i0j) 、 h(i0j,i0i)分别为第i到第j条集电线路、第j到第i条集电线路零序电流小波系数的Hausdorff单向距离;
Figure FDA0003210236520000012
Step3:通过式(4)求取各集电线路零序电流波形综合关联系数Hj,并根据式(5)形成零序电流波形综合关联系数矩阵HΣ
Figure FDA0003210236520000021
HΣ=[H1 … Hj … Hk … Hm],j=1,…,j,…,k,…m (5)
Step4:假设第k条集电 线路对应的元素Hk为矩阵HΣ的最大值,根据式(6)求取各集电线路综合关联系数商Skj,并根据式(7)计算综合关联系数商矩阵Sk
Figure FDA0003210236520000022
Sk=[Sk1 … Skj … 1 … Skm],j=1,…,j,…,k,…m (7)
Step5:设定阈值Sth,若矩阵Sk互关联系数均大于Sth,则第k条集电线路故障;若矩阵Sk互关联系数均小于Sth,则为汇集母线故障。
CN202010926138.XA 2020-09-07 2020-09-07 一种风电场集电线路单相接地故障选线方法 Active CN112083286B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010926138.XA CN112083286B (zh) 2020-09-07 2020-09-07 一种风电场集电线路单相接地故障选线方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010926138.XA CN112083286B (zh) 2020-09-07 2020-09-07 一种风电场集电线路单相接地故障选线方法

Publications (2)

Publication Number Publication Date
CN112083286A CN112083286A (zh) 2020-12-15
CN112083286B true CN112083286B (zh) 2021-09-28

Family

ID=73733211

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010926138.XA Active CN112083286B (zh) 2020-09-07 2020-09-07 一种风电场集电线路单相接地故障选线方法

Country Status (1)

Country Link
CN (1) CN112083286B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113189446B (zh) * 2021-04-28 2022-09-06 国网山东省电力公司泰安供电公司 一种小电流接地故障选线方法
CN113804978B (zh) * 2021-09-01 2024-01-23 青岛鼎信通讯股份有限公司 一种基于故障指示器的雷击波形识别方法
CN117375066A (zh) * 2023-09-05 2024-01-09 长江勘测规划设计研究有限责任公司 一种新能源电站集电线路单回多分支优化系统及其方法
CN117686843B (zh) * 2024-01-31 2024-04-26 昆明理工大学 基于波形高次累加重构的配电网故障选线方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101162838A (zh) * 2007-11-29 2008-04-16 昆明理工大学 一种利用小波包分解和相关分析的小电流接地系统故障选线方法
CN101242097A (zh) * 2008-03-11 2008-08-13 昆明理工大学 利用零模电流测后模拟的小电流接地系统故障选线方法
WO2014039317A2 (en) * 2012-09-06 2014-03-13 Leonard Flom Iris identification system and method
CN105529688A (zh) * 2016-02-25 2016-04-27 三峡大学 基于Hausdorff距离算法的变压器励磁涌流和故障差流识别方法
CN108572303A (zh) * 2018-04-25 2018-09-25 西南交通大学 一种基于随机矩阵和Hausdorff距离的配电网故障选线方法
CN109541395A (zh) * 2018-11-23 2019-03-29 三峡大学 一种特征频带下基于Hausdorff的配网区段定位方法
CN110988604A (zh) * 2019-12-30 2020-04-10 武汉理工大学 一种配电网单相接地故障选相方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101162838A (zh) * 2007-11-29 2008-04-16 昆明理工大学 一种利用小波包分解和相关分析的小电流接地系统故障选线方法
CN101242097A (zh) * 2008-03-11 2008-08-13 昆明理工大学 利用零模电流测后模拟的小电流接地系统故障选线方法
WO2014039317A2 (en) * 2012-09-06 2014-03-13 Leonard Flom Iris identification system and method
CN105529688A (zh) * 2016-02-25 2016-04-27 三峡大学 基于Hausdorff距离算法的变压器励磁涌流和故障差流识别方法
CN108572303A (zh) * 2018-04-25 2018-09-25 西南交通大学 一种基于随机矩阵和Hausdorff距离的配电网故障选线方法
CN109541395A (zh) * 2018-11-23 2019-03-29 三峡大学 一种特征频带下基于Hausdorff的配网区段定位方法
CN110988604A (zh) * 2019-12-30 2020-04-10 武汉理工大学 一种配电网单相接地故障选相方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于相似度比较的环状柔直配电网线路单极接地故障保护;戴志辉 等;《电力系统自动化》;20191210;第43卷(第23期);第107-115页 *

Also Published As

Publication number Publication date
CN112083286A (zh) 2020-12-15

Similar Documents

Publication Publication Date Title
CN112083286B (zh) 一种风电场集电线路单相接地故障选线方法
CN112147462A (zh) 一种基于深度学习的输电线路故障辨识方法
Devi et al. Detection of transmission line faults using discrete wavelet transform
CN109901013B (zh) 一种利用电流、电压突变量极性来判别配电网故障方向的方法
CN105119255A (zh) 一种基于故障态的光伏微电网故障隔离方法
CN103163430A (zh) 一种利用复小波和ann结合的谐振接地系统故障选线方法
CN107765139A (zh) 一种高准确率的谐振接地系统单相接地故障选线方法
Beniwal et al. Detection of anomalies in the quality of electricity supply
CN112816831B (zh) 风电场集电线单相接地故障定位方法
CN104808109A (zh) 基于录波数据的高压输电线路故障识别方法和系统
CN106230378A (zh) 一种光伏电站组串故障的诊断方法
Saleem et al. Faults detection and classification of HVDC transmission lines of using discrete wavelet transform
CN112083270B (zh) 一种基于相关系数的风电场集电线路单相接地故障选线方法
CN106932641A (zh) 一种基于改进hht变换的电网故障诊断方法
CN112505489A (zh) 风电场集电线故障定位方法
CN113447758B (zh) 一种风电场多分支集电线路单相接地故障测距方法
CN113777440A (zh) 一种光伏电站集电线路单相接地故障的定位方法
Suman et al. Detection of transmission line faults in the presence of solar PV generation using discrete wavelet
CN107462810A (zh) 一种适用于有源配电网的故障区段定位方法
CN116973683A (zh) 一种风电场集电线路单相接地故障区段定位方法及系统
CN112557950A (zh) 基于矩阵相似性的配电网谐振接地系统故障选线方法
CN112083284B (zh) 一种风电场集电线路故障区段与故障分支识别方法
CN112098773A (zh) 一种含分布式电源配电网相间短路故障定位方法
CN112649756A (zh) 风电场集电线单相接地故障测距方法、系统、介质及设备
Kulshrestha et al. A discrete wavelet transform and rule based decision tree based technique for identification of fault in utility grid network with wind energy

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