CN112595928A - 一种适用于单极大地运行的柔直系统接地故障测距的方法 - Google Patents

一种适用于单极大地运行的柔直系统接地故障测距的方法 Download PDF

Info

Publication number
CN112595928A
CN112595928A CN202011448732.9A CN202011448732A CN112595928A CN 112595928 A CN112595928 A CN 112595928A CN 202011448732 A CN202011448732 A CN 202011448732A CN 112595928 A CN112595928 A CN 112595928A
Authority
CN
China
Prior art keywords
mki
distance
fault
transmission line
calculating
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
CN202011448732.9A
Other languages
English (en)
Other versions
CN112595928B (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN202011448732.9A priority Critical patent/CN112595928B/zh
Publication of CN112595928A publication Critical patent/CN112595928A/zh
Application granted granted Critical
Publication of CN112595928B publication Critical patent/CN112595928B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08—Locating faults in cables, transmission lines, or networks
    • G01R31/088—Aspects of digital computing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08—Locating faults in cables, transmission lines, or networks
    • G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/083—Locating faults in cables, transmission lines, or networks according to type of conductors in cables, e.g. underground
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/52—Testing for short-circuits, leakage current or ground faults
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52—Outage or fault management, e.g. fault detection or location

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Locating Faults (AREA)

Abstract

本发明提供一种适用于单极大地运行的柔直系统接地故障测距的方法，包括以下步骤：确定直流输电线路的总长度，利用卡松公式计算线路特征阻抗以及其地模分量和差模分量；假设故障距离为lmki＝i·Δl；计算该假设故障距离lmki下每个采样点对应测距判据Difi(t)的值，然后求出该假设故障距离下所有采样点测距判据值的方差和Si 2并记录；以1为单位长度不断增加i的值并重复步骤(3)和步骤(4)，直到lmki＝i·Δl≥l时结束循环；在所有假设的故障距离中，计算出方差和Si 2最小时所对应的假设故障距离lmki就是实际故障距离。

Description

一种适用于单极大地运行的柔直系统接地故障测距的方法
技术领域
本发明属于电力系统及其自动化继电保护技术领域，涉及一种适用于柔性直流输电系统单极大地运行模式下金属性接地故障测距的方法。
背景技术
与传统高压直流输电方式相比，柔性直流输电方式在更具灵活性的同时兼具经济性好，环保性能优良等特点，在经济方面和技术方面都具有独特的优势。超、特高压直流输电系统在单极闭锁或检修等特殊工况下会运行于单极-大地方式，此时接地极的额定入地电流可达3kA以上。此外，在实际高压柔性直流工程中的输电线路一般采用架空输电线路，虽然采用架空输电线路能够大幅度节省输电成本，但是相对于直流电缆而言其发生各种直流故障的概率也相对较大。直流架空输电线路往往距离很长，跨越地区地形复杂多样且故障率高，准确可靠的故障测距技术对减轻巡线人员的工作负担，提高供电可靠性尤为重要。
发明内容
本发明公开了一种适用于柔性直流输电系统单极大地运行模式下金属性接地故障测距的方法，通过本方法的应用可以实现柔性直流输电系统单极大地运行模式下金属性接地故障测距。技术方案如下：
一种适用于单极大地运行的柔直系统接地故障测距的方法，包括以下步骤：
(1)确定直流输电线路的总长度l，根据直流输电线路几何参数，利用卡松公式计算出线路的特征阻抗zc(t)。
(2)保持与zc(t)相同的采样频率对故障电压um(t)和故障电流im(t)进行采样，数据窗长与Zc(t)长度相同，其采样点数记为M；
(3)假设故障距离为lmki＝i·Δl，其中Δl是设定的单位长度，其值按照波速v乘以采样周期T来设置；i是常整数且i的取值范围是
Figure BDA0002825883910000011
利用卡松公式计算得到长度为lmki的输电线路的传播系数amki(t)以及长度为2lmki的输电线路的传播系数a2mki(t)；
(4)计算该假设故障距离lmki下每个采样点对应测距判据Difi(t)的值，然后求出该假设故障距离下所有采样点测距判据值的方差和Si 2并记录；
测距判据为：
Difi(t)＝|Fwm(t)*a2mki(t)+Bwm(t)|
式中，
Figure BDA0002825883910000012
Fwm(t)和Bwm(t)分别为线路上的正向行波和反向行波，um(t)和im(t)分别表示直流输电线路首端m处的测量电压和测量电流；a2mki(t)表示长度2lmki的直流输电线路对应的传播系数；
所有采样点测距判据值的方差和：
Figure BDA0002825883910000021
其中，
Figure BDA0002825883910000022
(5)以1为单位长度不断增加i的值并重复步骤(3)和步骤(4)，直到lmki＝i·Δl≥l时结束循环，并进行步骤(6)；
(6)在所有假设的故障距离中，计算出方差和Si 2最小时所对应的假设故障距离lmki就是实际故障距离。
与现有技术相比，本发明仅利用单端电气量数据，原理简单易实现，成本低且具有较高的故障测距精度。另外，本发明不基于通信技术，无需安装通信模块且无通信同步问题。
附图说明
图1单极大地运行模式下的双端柔性直流输电系统
图2算法程序框图
具体实施方式
本发明提出一种适用于柔性直流输电系统单极大地运行模式下金属性接地故障测距的方法。本发明基于直流输电线路的频率相关模型，根据行波传播公式推导出柔性直流输电系统单极大地运方式下发生接地故障时，直流输电线路首端测量电压电流与传播系数和特征阻抗之间的关系，通过假设故障距离并判断是否满足关系式的方式实现故障测距，具有较高的测距精度。
本发明基于直流输电线路的频率相关模型，直流输电线路的特征阻抗Zc(f)和传播系数A(f)计算公式如下：
Figure BDA0002825883910000023
Figure BDA0002825883910000024
式中，R(f)和L(f)分别是单位长度线路的电阻和电感，其值会随着频率的变化而发生变化，具有频率相关特性。G(f)是单位长度线路的电导，其值很小，实际计算中可忽略不计。C是单位长度线路的电容，其值不会随着频率的变化而发生改变，因此不具有频率相关特性。当线路长度l和频率f给定时，R(f)，L(f)和C可以通过卡松公式计算得到。由此可以得到线路的特征阻抗的频域值Zc(f)和特定长度下的传播系数的频域值A(f)，再将其变化到时域即可得到特定频率下线路特征阻抗的时域值zc(t)以及特定长度和频率下线路传播系数的时域值a(t)。
技术方案主要包括以下具体步骤：
步骤A：根据行波传播公式推导出柔性直流输电系统单极大地运方式下发生金属性接地故障时，直流输电线路首端测量电压电流与传播系数和特征阻抗之间的关系。
以图1所示的单极大地运行模式下的双端柔性直流输电系统模型为基础进行分析，图1中线路总长为lmn，输电线路首端m侧到故障点k处的距离为lmk，im(t)和um(t)分别是m端的测量电流和测量电压，in(t)是n端电流，uk(t)是故障点k处的电压，imk(t)是线路首端m传向故障点k的电流。
当线路中某一点k处发生接地故障时，由行波的传播公式可以得到输电线路首端m侧的正向电压行波和故障点k处的正向电压行波的关系为：
uk(t)+zc(t)*imk(t)＝[um(t)+zc(t)*im(t)]*amk(t) (3)
输电线路首端m侧的反向电压行波和故障点k处的反向电压行波的关系为：
uk(t)-zc(t)*imk(t)＝[um(t)-zc(t)*im(t)]/*amk(t) (4)
式中amk(t)是故障点k到输电线路首端m侧的传播系数的时域表达式。zc(t)是输电线路特征阻抗的时域值。
若输电线路在k点处发生的是金属性接地故障，那么故障点的电压uk(t)会迅速降为零，此时式(3)和式(4)变为：
zc(t)*imk(t)＝[um(t)+zc(t)*im(t)]*amk(t) (5)
-zc(t)*imk(t)＝[um(t)-zc(t)*im(t)]/*amk(t) (6)
式(5)和式(6)相加可以得到：
-[um(t)-zc(t)*im(t)]/*amk(t)＝[um(t)+zc(t)*im(t)]*amk(t) (7)
因此在发生金属性接地故障时有：
-[um(t)-zc(t)*im(t)]＝[um(t)+zc(t)*im(t)]*a2mk(t) (8)
式中*和/*分别表示卷积与反卷积运算。a2mk(t)是长度为2lmk的输电线路对应的传播系数的时域表达式。
令-[um(t)-zc(t)*im(t)]＝Bwm(t)，[um(t)+zc(t)*im(t)]＝Fwm(t)，其中Bwm(t)和Fwm(t)分别表示线路上的反向行波和正向行波，则公式(8)变为：
Fwm(t)*a2mk(t)+Bwm(t)＝0 (9)
在公式(9)中，反向行波的负值中不包含传播系数项，因此其值受故障距离的影响很小。正向行波与传播系数的卷积中包含传播系数项，因此其值会随着故障距离的改变而发生大幅度的变化。
步骤B：根据发生金属性接地故障后正向行波与反向行波之间的关系，得到求解故障距离的方法。
在求解故障距离时，可通过假设不同的故障距离lmk并利用卡松公式求出zc(t)和此故障距离lmk下的传播系数a(t)，进而求得此故障距离lmk下等式(9)左侧的值。在不同的假设故障距离求出的结果中，公式(9)中左侧值的方差和最小时所对应的假设故障距离就是实际故障距离。
下面介绍具体的求解过程，计算流程图如图2所示。
(1)确定直流输电线路的总长度l，根据直流输电线路几何参数，利用卡松公式计算出线路的特征阻抗zc(t)。
(2)保持与zc(t)相同的采样频率对故障电压um(t)和故障电流im(t)进行采样，数据窗长与Zc(t)长度相同，其采样点数记为M。
(3)假设故障距离为lmki＝i·Δl，其中Δl是设定的单位长度，其值按照波速v乘以采样周期T来设置。i是常整数且i的取值范围是
Figure BDA0002825883910000041
利用卡松公式计算得到长度为lmki的输电线路的传播系数amki(t)以及长度为2lmki的输电线路的传播系数a2mki(t)。
(4)计算该假设故障距离lmki下每个采样点对应测距判据Difi(t)的值，然后求出该假设故障距离下所有采样点测距判据值的方差和Si 2并记录。
测距判据为：
Difi(t)＝|Fwm(t)*a2mki(t)+Bwm(t)|
式中
Figure BDA0002825883910000042
Fwm(t)和Bwm(t)分别为线路上的正向行波和反向行波.
um(t)和im(t)分别表示直流输电线路首端m处的测量电压和测量电流；a2mki(t)表示长度2lmki的直流输电线路对应的传播系数。
所有采样点测距判据值的方差和为：
Figure BDA0002825883910000043
其中，
Figure BDA0002825883910000044
(5)以1为单位长度不断增加i的值并重复步骤(3)和步骤(4)，直到lmki＝i·Δl≥l时结束循环，并进行步骤(6)。
(6)在所有假设的故障距离中，计算出方差和Si 2最小时所对应的假设故障距离lmki就是实际故障距离。

Claims (1)

1.一种适用于单极大地运行的柔直系统接地故障测距的方法，包括以下步骤：
(1)确定直流输电线路的总长度l，根据直流输电线路几何参数，利用卡松公式计算出线路的特征阻抗zc(t)。
(2)保持与zc(t)相同的采样频率对故障电压um(t)和故障电流im(t)进行采样，数据窗长与Zc(t)长度相同，其采样点数记为M；
(3)假设故障距离为lmki＝i·Δl，其中Δl是设定的单位长度，其值按照波速v乘以采样周期T来设置；i是常整数且i的取值范围是
Figure FDA0002825883900000011
利用卡松公式计算得到长度为lmki的输电线路的传播系数amki(t)以及长度为2lmki的输电线路的传播系数a2mki(t)；
(4)计算该假设故障距离lmki下每个采样点对应测距判据Difi(t)的值，然后求出该假设故障距离下所有采样点测距判据值的方差和Si 2并记录；
测距判据为：
Difi(t)＝|Fwm(t)*a2mki(t)+Bwm(t)|
式中，
Figure FDA0002825883900000012
Fwm(t)和Bwm(t)分别为线路上的正向行波和反向行波，um(t)和im(t)分别表示直流输电线路首端m处的测量电压和测量电流；a2mki(t)表示长度2lmki的直流输电线路对应的传播系数；
所有采样点测距判据值的方差和：
Figure FDA0002825883900000013
其中，
Figure FDA0002825883900000014
(5)以1为单位长度不断增加i的值并重复步骤(3)和步骤(4)，直到lmki＝i·Δl≥l时结束循环，并进行步骤(6)；
(6)在所有假设的故障距离中，计算出方差和Si 2最小时所对应的假设故障距离lmki就是实际故障距离。
CN202011448732.9A 2020-12-09 2020-12-09 一种适用于单极大地运行的柔直系统接地故障测距的方法 Active CN112595928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011448732.9A CN112595928B (zh) 2020-12-09 2020-12-09 一种适用于单极大地运行的柔直系统接地故障测距的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011448732.9A CN112595928B (zh) 2020-12-09 2020-12-09 一种适用于单极大地运行的柔直系统接地故障测距的方法

Publications (2)

Publication Number Publication Date
CN112595928A true CN112595928A (zh) 2021-04-02
CN112595928B CN112595928B (zh) 2022-07-05

Family

ID=75192014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011448732.9A Active CN112595928B (zh) 2020-12-09 2020-12-09 一种适用于单极大地运行的柔直系统接地故障测距的方法

Country Status (1)

Country Link
CN (1) CN112595928B (zh)

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1123639A (ja) * 1997-07-03 1999-01-29 Toko Electric Co Ltd 片極間欠弧光地絡試験装置
CN201993433U (zh) * 2010-11-15 2011-09-28 中国南方电网有限责任公司超高压输电公司检修试验中心 高压直流输电接地极线路行波测距系统
CN102288869A (zh) * 2011-05-10 2011-12-21 山东大学 一种输电线路单端行波故障测距方法
CN102508121A (zh) * 2011-11-08 2012-06-20 西安交通大学 多端柔性直流输电系统的直流线路单端故障定位方法
CN102570428A (zh) * 2012-02-29 2012-07-11 山东电力集团公司莱芜供电公司 基于电子式互感器微分输出的故障测距与距离保护方法
CN103293439A (zh) * 2013-05-19 2013-09-11 国家电网公司 基于分布参数测量阻抗幅值特性线路单相接地故障单端测距方法
JP2015137865A (ja) * 2014-01-20 2015-07-30 Ｊｘ日鉱日石エネルギー株式会社 地絡検出装置
CN105388398A (zh) * 2015-12-25 2016-03-09 国网四川省电力公司电力科学研究院 一种特高压直流输电系统接地极引线故障监测方法
CN105548802A (zh) * 2015-12-04 2016-05-04 昆明理工大学 一种基于故障行波沿线分布特性的t接线路三端不同步故障测距方法
CN105738770A (zh) * 2016-04-01 2016-07-06 昆明理工大学 一种基于故障行波沿线分布特性的线缆混合线路单端行波测距方法
CN105891672A (zh) * 2016-04-01 2016-08-24 昆明理工大学 一种基于故障行波沿线分布特性的含upfc线路单端行波测距方法
CN106019088A (zh) * 2016-05-17 2016-10-12 中国南方电网有限责任公司超高压输电公司检修试验中心 一种直流接地极线路测距方法
CN109031049A (zh) * 2018-09-07 2018-12-18 昆明理工大学 一种基于单侧谐波量的高压直流接地极线路故障测距方法
CN109142974A (zh) * 2018-08-31 2019-01-04 天津大学 一种适用于柔性直流馈入的交流线路单端故障测距方法
CN109541392A (zh) * 2018-10-15 2019-03-29 天津大学 一种适用于柔性直流输电系统的单端故障测距方法
CN110514964A (zh) * 2019-09-05 2019-11-29 国网内蒙古东部电力有限公司检修分公司 一种特高压直流输电系统接地极线路故障测距方法

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1123639A (ja) * 1997-07-03 1999-01-29 Toko Electric Co Ltd 片極間欠弧光地絡試験装置
CN201993433U (zh) * 2010-11-15 2011-09-28 中国南方电网有限责任公司超高压输电公司检修试验中心 高压直流输电接地极线路行波测距系统
CN102288869A (zh) * 2011-05-10 2011-12-21 山东大学 一种输电线路单端行波故障测距方法
CN102508121A (zh) * 2011-11-08 2012-06-20 西安交通大学 多端柔性直流输电系统的直流线路单端故障定位方法
CN102570428A (zh) * 2012-02-29 2012-07-11 山东电力集团公司莱芜供电公司 基于电子式互感器微分输出的故障测距与距离保护方法
CN103293439A (zh) * 2013-05-19 2013-09-11 国家电网公司 基于分布参数测量阻抗幅值特性线路单相接地故障单端测距方法
JP2015137865A (ja) * 2014-01-20 2015-07-30 Ｊｘ日鉱日石エネルギー株式会社 地絡検出装置
CN105548802A (zh) * 2015-12-04 2016-05-04 昆明理工大学 一种基于故障行波沿线分布特性的t接线路三端不同步故障测距方法
CN105388398A (zh) * 2015-12-25 2016-03-09 国网四川省电力公司电力科学研究院 一种特高压直流输电系统接地极引线故障监测方法
CN105738770A (zh) * 2016-04-01 2016-07-06 昆明理工大学 一种基于故障行波沿线分布特性的线缆混合线路单端行波测距方法
CN105891672A (zh) * 2016-04-01 2016-08-24 昆明理工大学 一种基于故障行波沿线分布特性的含upfc线路单端行波测距方法
CN106019088A (zh) * 2016-05-17 2016-10-12 中国南方电网有限责任公司超高压输电公司检修试验中心 一种直流接地极线路测距方法
CN109142974A (zh) * 2018-08-31 2019-01-04 天津大学 一种适用于柔性直流馈入的交流线路单端故障测距方法
CN109031049A (zh) * 2018-09-07 2018-12-18 昆明理工大学 一种基于单侧谐波量的高压直流接地极线路故障测距方法
CN109541392A (zh) * 2018-10-15 2019-03-29 天津大学 一种适用于柔性直流输电系统的单端故障测距方法
CN110514964A (zh) * 2019-09-05 2019-11-29 国网内蒙古东部电力有限公司检修分公司 一种特高压直流输电系统接地极线路故障测距方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
冯志翔: "《基于行波原理的配电网故障选线与定位方法》", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *
张帆: "《基于单端暂态行波的接地故障测距与保护研究》", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *

Also Published As

Publication number Publication date
CN112595928B (zh) 2022-07-05

Similar Documents

Publication Publication Date Title
CN103176107B (zh) 高压直流输电线路混合故障测距方法
CN109541392B (zh) 一种适用于柔性直流输电系统的单端故障测距方法
CN111313932B (zh) 一种基于电磁时间反演的传递函数相关性的传输线及网络故障定位方法
CN101666848B (zh) 基于行波的故障定位方法
CN103901324A (zh) 一种基于单端故障信息的配电网混合线路组合式测距方法
CN105182184B (zh) 一种基于沿线方向行波分解的配网线缆混合直配线路故障测距方法
CN111766470B (zh) 高压直流输电线路的故障定位方法、系统及直流输电线路
CN108196166B (zh) 基于故障波形起始点到达时差的双端行波测距方法
Ha et al. A novel principle of single-ended fault location technique for EHV transmission lines
CN109387744B (zh) 基于奇异值分解的配网线路故障点定位方法及装置
CN109444657B (zh) 一种配电网高阻接地故障区段定位方法
CN111948493A (zh) 一种mmc-hvdc直流输电线路故障定位方法
CN109100605B (zh) 利用故障边界条件的高压电缆单相接地故障的单端定位方法
CN103869221A (zh) 基于sv网络采样的双回线路单相接地故障距离测量方法
CN110161375B (zh) 一种基于分布电阻参数的高压直流输电线路计算模型
CN104062550A (zh) 基于二分搜索的双回线路非同名相跨线接地故障测距方法
CN110749835A (zh) 基于卡尔曼滤波器的输电线路故障定位方法
CN110514964B (zh) 一种特高压直流输电系统接地极线路故障测距方法
CN109324262A (zh) 一种基于tt变换和波速优化的输电线路故障测距方法
CN112595928B (zh) 一种适用于单极大地运行的柔直系统接地故障测距的方法
CN112595929B (zh) 一种适用于真双极柔直系统单极接地故障测距的方法
CN110163148A (zh) 一种电动汽车直流充电畸变信号自适应识别方法
CN113945806A (zh) 一种柔性直流输电线路单端故障测距方法及系统
CN113406441A (zh) 基于聚类和迭代算法的柔性直流电网故障测距方法
CN115542087B (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