CN112906189B - 一种计算线路母线至刀闸引线对地电容的方法 - Google Patents

一种计算线路母线至刀闸引线对地电容的方法 Download PDF

Info

Publication number
CN112906189B
CN112906189B CN202110066494.3A CN202110066494A CN112906189B CN 112906189 B CN112906189 B CN 112906189B CN 202110066494 A CN202110066494 A CN 202110066494A CN 112906189 B CN112906189 B CN 112906189B
Authority
CN
China
Prior art keywords
capacitance
switch
ground
voltage
bus
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
CN202110066494.3A
Other languages
English (en)
Other versions
CN112906189A (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.)
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hebei 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 State Grid Corp of China SGCC, State Grid Hebei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202110066494.3A priority Critical patent/CN112906189B/zh
Publication of CN112906189A publication Critical patent/CN112906189A/zh
Application granted granted Critical
Publication of CN112906189B publication Critical patent/CN112906189B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G01R27/26Measuring inductance or capacitance; Measuring quality factor, e.g. by using the resonance method; Measuring loss factor; Measuring dielectric constants ; Measuring impedance or related variables
    • G01R27/2605Measuring capacitance
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明提供了一种计算线路母线侧至刀闸引线对地电容的方法,本技术方案包括线路母线对地电容计算、母线开关引线对地电容计算、停送电过程电压实测数据和结论输出;本发明考虑了母线侧至刀闸的引线电容,为搭建同塔双回输电线路的感应电模型、感应电压仿真,提供了更准确的数据,与实际测量值偏差更小,为感应电压的计算提供一个理论基础,从而应对感应电压问题,给检修人员提供一个安全的保障。

Description

一种计算线路母线至刀闸引线对地电容的方法
技术领域
本发明涉及电力***领域,具体地说,涉及一种计算线路母线至刀闸引线对地电容的方法。
背景技术
同塔多回架空输电线路的出现,虽然减少了线路的建设改造成本,以及缓解了输电走廊用地的问题,但随之带来了线路之间空间间距的减小,由于线路间距离近,因此容易在线路上感应出大的感应电压和感应电流。特别是在其中一回线路进行停电检修时,由于感应电压、感应电流的存在,对检修人员的安全是一个严重的威胁,线路间的耦合参数主要是根据线路之间的间距、线路本身的参数所决定的,因此要计算出线路之间的感应电压,还需对互感、电容等线路之间参数的计算原理进行了解。
目前耦合参数计算方法多基于简单的短线集中参数模型(忽略了对地电容),且未考虑线路两端采样的不同步问题,并且正向计算线路中刀闸开关的引线对地电容比较困难,没有能够直接计算引线对地电容的方法。
发明内容
本发明的目的在于建立一种计算线路母线至刀闸引线对地电容的方法,为感应电压的计算提供一个理论基础,从而应对感应电压问题,给检修人员提供一个安全的保障。
为了实现上述技术方案,本发明提供了一种计算母线至刀闸引线对地电容的方法,所述计算方法包括如下步骤:
步骤1):线路母线对地电容计算;向量测量单元(PMU)和数据采集与监控(SCADA)量测相结合,SCADA数据采样间隔普遍为1-10秒,线路两端量SCADA测量可能存在不同步时间较长,甚至量测量分属不同运行工况的情形,从而影响线路辨识精度。PMU数据间隔一般为10ms或20ms,其两端量测量基本在在同一时间断面,且具有提供相角测量的优势,因而采用PMU数据的在线量测得到了越来越广泛的应用。在双回并行输电线路的正、负序集中参数电路模型中,不同回路间不存在耦合,即各回线路的正、负序集中参数可以进行独立计算。每回线中正负零序间参数不存在序间耦合,互相独立;两回线路间存在正序耦合参数、负序耦合参数及零序耦合参数,但正负零序也都是相互独立,不存在序间耦合。主要是利用去除相角信息的PMU数据通过两端获取的PMU数据中的电压、电流幅值和功率信息推导相关辨识方程,获取了单时刻的非线性辨识方程,并利用多时刻的量测数据进行最小二乘辨识;
步骤2):母线开关引线对地电容计算;通过送电或停电过程每一个母线开关状态变化后实测的SCADA数据录波电压数据,从线路中只有一个开关处于热备用状态开始计算,该状态的电压实测值与单相相电压的比值等于该开关自身电容与开关自身对地电容、母线自身对地电容、母线上所带电压互感器电容、开关引线对地电容之和,只有开关引线对地电容未知,由此计算出开关引线对地电容的值,以此类推,每有一个开关转热备用,即可求得该开关引线对地电容;
步骤3):停送电过程电压实测数据;导出停送电过程中测量的电压数据,同时通过反推算法搭建模型并仿真电压,将实测电压与仿真电压进行对比,计算出误差;
步骤4):结论输出;根据计算得出母线开关引线对地电容,其有助于输电线路模型的搭建。
优选的,所述步骤1具体分为以下3个部分:
1)釆用正序分量描述,针对输电线路单回线路,建立集中参数π型等值模型;
2)考虑任一k时刻,基于输电线路两端PMU量测的量,去除其中的相角信息,只采用电压和电流幅值量测和有功及无功量测进行参数辨识;此时,假设N端电压相角为0,M端电压相角为a,则有:
Figure BDA0002904217690000031
Figure BDA0002904217690000032
其中,UM、UN分别为M、N端电压的幅值,UNR、UNI分别为N端电压相量实部和虚部,UMR、UMI分别为M端电压相量实部和虚部;
3)对于采用单时刻的辨识,通过计算可以得到k时刻电阻、电抗及对地导纳的辨识值:
Figure BDA0002904217690000041
Figure BDA0002904217690000042
YC.k=yc.k
通过电抗的辨识值从而的得到电容矩阵,其中对角线的值即为输电线路的对地电容。
优选的,所述步骤2具体分为以下3个部分:
1)线路母线及各开关均冷备用时,设UB为运行母线电压,Ub为停运母线感应电压,CBb为两母线间互电容,Cb为停运母线自身对地电容;其中UB为已知值;Ub为已知实测值;CBb可由所述步骤1的线路母线对地电容计算得到,约为0;Cb可由所述步骤1的线路母线对地电容计算得到;
2)线路母线间有一开关转热备用时,则计算中需加入一个开关断口电容和一个开关引线对地电容;设C5011为加入的开关断口电容,C引5011为加入的开关引线对地电容;其中C5011为已知值,C引5011为待求值;依据计算公式
Figure BDA0002904217690000043
即可求得C引5011的值;
3)再将一开关转热备用时,与上一步骤相比,则需要在计算中再加入一个开关断口电容和一个开关引线对地电容;设C5042为再次加入的开关断口电容,C引5042为再次加入的开关引线对地电容;其中C5042为已知值,C引5042为待求值;同上一步骤的方法可求得C引5042的值。
优选的,所述步骤4中由于第一个开关状态改变时电压过小,第一个开关状态的电压值不计入结果。
同现有技术相比,本发明的有益效果体现在:
(1)本发明的技术方案考虑了母线侧至刀闸的引线电容,为搭建同塔双回输电线路的感应电模型、感应电压仿真,提供了更准确的数据,与实际测量值偏差更小,为感应电压的计算提供一个理论基础,从而应对感应电压问题,给检修人员提供一个安全的保障.
(2)本发明的技术方案根据计算得出母线开关引线对地电容,计算结果搭建模型与实测值的误差较小,对输电线路的各项参数计算有很大帮助,有助于进行输电线路模型搭建及一系列后续研究。
附图说明
图1为本发明的一种计算线路母线侧至刀闸引线对地电容的方法流程图;
图2为线路母线静电感应电压模型示意图;
图3为本发明中所述步骤1的输电线路单回线路的集中参数π型等值模型示意图;
图4为本发明中所述步骤2的母线及各开关均冷备用时的模型示意图;
图5为本发明中所述步骤2的母线间有一开关转热备用时的模型示意图;
图6为本发明中所述步骤2的再将一开关转热备用时的模型示意图。
具体实施方式
为了能够进一步了解本发明的结构、特征及其他目的,现结合所附较佳实施例附以附图详细说明如下,本附图所说明的实施例仅用于说明本发明的技术方案,并非限定本发明。
实施例一:
如图1所示,图1为本发明的一种计算线路母线侧至刀闸引线对地电容的方法流程图;所述计算方法包括如下步骤:
步骤1):线路母线对地电容计算;向量测量单元(PMU)和数据采集与监控(SCADA)量测相结合;SCADA数据采样间隔为1-10秒,PMU数据间隔为10ms或20ms;利用去除相角信息的PMU数据,通过两端获取的PMU数据中的电压、电流幅值和功率信息推导相关辨识方程,获取单时刻的非线性辨识方程,并利用多时刻的量测数据进行最小二乘辨识;
步骤2):母线开关引线对地电容计算;通过送电或停电过程每一个母线开关状态变化后实测的SCADA数据录波电压数据,从线路中只有一个开关处于热备用状态开始计算,该状态的电压实测值与单相相电压的比值等于该开关自身电容与开关自身对地电容、母线自身对地电容、母线上所带电压互感器电容、开关引线对地电容之和,只有开关引线对地电容未知,由此计算出开关引线对地电容的值,以此类推,每有一个开关转热备用,即可求得该开关引线对地电容;
步骤3):停送电过程电压实测数据;导出停送电过程中测量的电压数据,同时通过反推算法搭建模型并仿真电压,将实测电压与仿真电压进行对比,计算出误差;
步骤4):结论输出;根据计算得出母线开关引线对地电容,其有助于输电线路模型的搭建。
实施例二:
请参见图2,图2为线路母线静电感应电压模型示意图;本发明的技术原理为:
在图2中,a为正常运行母线,b为停运母线,Usb为b母线感应电压,Ua为a母线电压,Cab为两条母线间互电容(耦合电容),Cb为b母线自身对地电容;线路母线静电感应电压的数学模型为
Figure BDA0002904217690000071
本发明计算的感应电压分为简单两线交流架空线路感应电压和停运检修回路(首末双端均不接地)感应电压。其中,简单两线交流架空线路感应电压包括电磁感应电压和静电感应电压,电磁感应电压由两回线间互电感引起,静电感应电压由两回线间互电容引起;停运检修回路(首末双端均不接地)感应电压只包括静电感应电压,停运检修线路中,首末双端未接地,不构成回路,感应电流为零,因此感应电压只为静电感应电压。
实施例三:
如图3所示,图3为本发明中所述步骤1的输电线路单回线路的集中参数π型等值模型示意图;所述步骤1具体分为以下3个步骤:
1)釆用正序分量描述,针对输电线路单回线路,建立集中参数π型等值模型,如图3所示;
2)考虑任一k时刻,基于输电线路两端PMU量测的量,去除其中的相角信息,只采用电压和电流幅值量测和有功及无功量测进行参数辨识;此时,假设N端电压相角为0,M端电压相角为a,则有:
Figure BDA0002904217690000081
Figure BDA0002904217690000082
其中,UM、UN分别为M、N端电压的幅值,UNR、UNI分别为N端电压相量实部和虚部,UMR、UMI分别为M端电压相量实部和虚部;
3)对于采用单时刻的辨识,通过计算可以得到k时刻电阻、电抗及对地导纳的辨识值:
Figure BDA0002904217690000083
Figure BDA0002904217690000084
YC.k=yc.k
通过电抗的辨识值从而的得到电容矩阵,其中对角线的值即为输电线路的对地电容。
实施例四:
请参见图4、图5和图6,图4为本发明中所述步骤2的母线及各开关均冷备用时的模型示意图;图5为本发明中所述步骤2的母线间有一开关转热备用时的模型示意图;图6为本发明中所述步骤2的再将一开关转热备用时的模型示意图;在所述步骤2母线开关引线对地电容的计算中,以母线静电感应模型为基础,得到各串开关引线对地电容计算方法;所述步骤2具体分为以下3个部分:
1)当线路母线及各开关均冷备用时,线路模型如图4所示;设UB为运行母线电压,Ub为停运母线感应电压,CBb为两母线间互电容,Cb为停运母线自身对地电容;其中UB为已知值;Ub为已知实测值;CBb可由所述步骤1的线路母线对地电容计算得到,约为0;Cb可由所述步骤1的线路母线对地电容计算得到;
2)当线路母线间有一开关转热备用时,线路模型如图5所示;计算中需加入一个开关断口电容和一个开关引线对地电容;设C5011为加入的开关断口电容,C引5011为加入的开关引线对地电容;其中C5011为已知值,C引5011为待求值;依据计算公式
Figure BDA0002904217690000091
即可求得C引5011的值;
3)再将一开关转热备用时,线路模型如图6所示;与上一步骤相比,需要在计算中再加入一个开关断口电容和一个开关引线对地电容;设C5042为再次加入的开关断口电容,C引5042为再次加入的开关引线对地电容;其中C5042为已知值,C引5042为待求值;同上一步骤的方法可求得C引5042的值。
实施例五:
在本发明中所述步骤3的停送电过程电压实测数据;导出停送电过程中测量的电压数据,同时通过反推算法搭建模型并仿真电压,将实测电压与仿真电压进行对比,计算出误差;以石北站母线送电过程为例,通过以上反推算法搭建模型与实测电压对比如下:
Figure BDA0002904217690000092
Figure BDA0002904217690000101
Figure BDA0002904217690000102
Figure BDA0002904217690000111
Figure BDA0002904217690000112
Figure BDA0002904217690000121
实施例六:
在本发明中所述步骤4的结论输出;由于第一个开关状态改变时电压过小,此时感应电压不计入结果;其余开关状态改变时根据计算得出母线开关引线对地电容;由于本发明的计算结果与实测值的误差相对较小,对输电线路的各项参数计算有很大的帮助,有助于进行输电线路模型搭建及一系列后续的研究。
与传统技术方案相比,本发明的技术效果为:
本发明的技术方案考虑了母线侧至刀闸的引线电容,为搭建同塔双回输电线路的感应电模型、感应电压仿真,提供了更准确的数据,与实际测量值偏差更小,为感应电压的计算提供一个理论基础,从而应对感应电压问题,给检修人员提供一个安全的保障。
需要声明的是,上述发明内容及具体实施方式意在证明本发明所提供技术方案的实际应用,不应解释为对本发明保护范围的限定。本领域技术人员在本发明的精神和原理内,当可作各种修改、等同替换或改进。本发明的保护范围以所附权利要求书为准。

Claims (4)

1.一种计算线路母线侧至刀闸引线对地电容的方法,其特征在于,所述计算方法包括如下步骤:
步骤1):线路母线对地电容计算;相量测量单元PMU和数据采集与监控SCADA量测相结合;SCADA数据采样间隔为1-10秒,PMU数据间隔为10ms或20ms;利用去除相角信息的PMU数据,通过两端获取的PMU数据中的电压、电流幅值和功率信息推导相关辨识方程,获取单时刻的非线性辨识方程,并利用多时刻的量测数据进行最小二乘辨识;
步骤2):母线开关引线对地电容计算;通过送电或停电过程每一个母线开关状态变化后实测的SCADA数据录波电压数据,从线路中只有一个开关处于热备用状态开始计算,该状态的电压实测值与单相相电压的比值等于该开关自身电容与开关自身对地电容、母线自身对地电容、母线上所带电压互感器电容、开关引线对地电容之和,只有开关引线对地电容未知,由此计算出开关引线对地电容的值,以此类推,每有一个开关转热备用,即可求得该开关引线对地电容;
步骤3):停送电过程电压实测数据;导出停送电过程中测量的电压数据,同时通过反推算法搭建模型并仿真电压,将实测电压与仿真电压进行对比,计算出误差;
步骤4):结论输出;根据计算得出母线开关引线对地电容,其有助于输电线路模型的搭建。
2.如权利要求1所述的计算方法,其特征在于,所述步骤1具体分为以下3个部分:
1)釆用正序分量描述,针对输电线路单回线路,建立集中参数π型等值模型;
2)考虑任一k时刻,基于输电线路两端PMU量测的量,去除其中的相角信息,只采用电压和电流幅值量测和有功及无功量测进行参数辨识;此时,假设N端电压相角为0,M端电压相角为α,则有:
Figure FDA0003617522440000021
Figure FDA0003617522440000022
其中,UM、UN分别为M、N端电压的幅值,UNR、UNI分别为N端电压相量实部和虚部,UMR、UMI分别为M端电压相量实部和虚部;
3)对于采用单时刻的辨识,通过计算可以得到k时刻电阻、电抗及对地导纳的辨识值:
Figure FDA0003617522440000023
Figure FDA0003617522440000024
YC.k=yc.k
通过电抗的辨识值从而的得到电容矩阵,其中对角线的值即为输电线路的对地电容。
3.如权利要求1所述的计算方法,其特征在于,所述步骤2具体分为以下3个部分:
1)线路母线及各开关均冷备用时,设UB为运行母线电压,Ub为停运母线感应电压,CBb为两母线间互电容,Cb为停运母线自身对地电容;其中UB为已知值;Ub为已知实测值;CBb可由所述步骤1的线路母线对地电容计算得到,约为0;Cb可由所述步骤1的线路母线对地电容计算得到;
2)线路母线间有一开关转热备用时,则计算中需加入一个开关断口电容和一个开关引线对地电容;设C5011为加入的开关断口电容,C引5011为加入的开关引线对地电容;其中C5011为已知值,C引5011为待求值;依据计算公式
Figure FDA0003617522440000031
即可求得C引5011的值;
3)再将一开关转热备用时,与上一步骤相比,则需要在计算中再加入一个开关断口电容和一个开关引线对地电容;设C5042为再次加入的开关断口电容,C引5042为再次加入的开关引线对地电容;其中C5042为已知值,C引5042为待求值;同上一步骤的方法可求得C引5042的值。
4.如权利要求1所述的计算方法,其特征在于,所述步骤4中由于第一个开关状态改变时电压过小,第一个开关状态的电压值不计入结果。
CN202110066494.3A 2021-01-19 2021-01-19 一种计算线路母线至刀闸引线对地电容的方法 Active CN112906189B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110066494.3A CN112906189B (zh) 2021-01-19 2021-01-19 一种计算线路母线至刀闸引线对地电容的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110066494.3A CN112906189B (zh) 2021-01-19 2021-01-19 一种计算线路母线至刀闸引线对地电容的方法

Publications (2)

Publication Number Publication Date
CN112906189A CN112906189A (zh) 2021-06-04
CN112906189B true CN112906189B (zh) 2022-06-24

Family

ID=76115091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110066494.3A Active CN112906189B (zh) 2021-01-19 2021-01-19 一种计算线路母线至刀闸引线对地电容的方法

Country Status (1)

Country Link
CN (1) CN112906189B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5623768B2 (ja) * 2010-03-29 2014-11-12 中国電力株式会社 対地静電容量測定装置および対地静電容量測定方法
CN108181500B (zh) * 2017-12-28 2020-02-14 中国南方电网有限责任公司超高压输电公司南宁局 一种输变电设备合闸无压定值确定方法
CN108680795B (zh) * 2018-05-15 2020-08-04 石家庄科林电气股份有限公司 基于一阶电路全响应的直流***对地电容计算方法
CN110501619B (zh) * 2019-08-30 2021-11-12 国网山西省电力公司电力科学研究院 一种高频响应分压测量装置

Also Published As

Publication number Publication date
CN112906189A (zh) 2021-06-04

Similar Documents

Publication Publication Date Title
Liao et al. Online optimal transmission line parameter estimation for relaying applications
CN110850236A (zh) 一种基于参数估计的配电网故障定位方法
Cordova et al. Fault location identification in smart distribution networks with distributed generation
CN105388396B (zh) 一种用序有功增量电流方向追溯电压暂降源的方法
CN107271841B (zh) 基于正序电流差的同塔双回t型输电线路故障测距方法
CN102901891A (zh) 一种电力***变压器保护装置的柔性仿真测试平台
Li et al. Fault location algorithm in distribution networks considering distributed capacitive current
Yarong et al. Time difference method to calculate the optimal differential point of half-wavelength AC transmission line differential protection
CN112906189B (zh) 一种计算线路母线至刀闸引线对地电容的方法
Zhang et al. A segmented network method based faulted line selection strategy for single-phase earth fault in small current grounding distribution network
Li et al. Augmented state estimation method for fault location based on on-line parameter identification of PMU measurement data
Xie et al. Method of T-connection line parameter estimation in distribution network based on PMU and AMI
Liao et al. Identification of Fault Line Selection and Section for Single-Phase Ground Fault in Small Current Grounding System
Quan et al. The real-time digital simulation system for distribution network
CN113933648A (zh) 配电网继电保护向量检查和动作逻辑校验测试装置及方法
Liu et al. Fault detection and location of microgrid based on distributed decision
Presadă et al. An algorithm for improving the power system state estimation using PMU measurements
Fan et al. Fault interval judgment of urban distribution grid based on edge computing of distribution Internet of Things
Yong et al. Recognition of Single-Phase-to-Ground Fault Section of Distribution Line Based Phase Current Variation
Su et al. Investigation of the use of electromagnetic transient models for transient stability simulation
Li et al. A wide‐area relaying protection system based on multi‐point measurement information
Lin et al. An Accurate Fault Location Method for Distribution Lines Based on Coupling Parameter Identification Using Two Terminals Data
Chen et al. The development of anti-ground fault suppression test device for smart meters
Shang et al. A new phase mode transformation matrix
Zelenskii et al. Detection and identification of bad power measurements in distribution networks using state estimation by synchronized voltages and currents

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