CN103149513B - 一种变压器局部放电超声源重建的定位方法及装置 - Google Patents

一种变压器局部放电超声源重建的定位方法及装置 Download PDF

Info

Publication number
CN103149513B
CN103149513B CN201310060568.8A CN201310060568A CN103149513B CN 103149513 B CN103149513 B CN 103149513B CN 201310060568 A CN201310060568 A CN 201310060568A CN 103149513 B CN103149513 B CN 103149513B
Authority
CN
China
Prior art keywords
transformer
source
sound
ultrasonic
discharge
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
CN201310060568.8A
Other languages
English (en)
Other versions
CN103149513A (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.)
Institute of Electrical Engineering of CAS
Original Assignee
Institute of Electrical Engineering of CAS
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 Institute of Electrical Engineering of CAS filed Critical Institute of Electrical Engineering of CAS
Priority to CN201310060568.8A priority Critical patent/CN103149513B/zh
Publication of CN103149513A publication Critical patent/CN103149513A/zh
Application granted granted Critical
Publication of CN103149513B publication Critical patent/CN103149513B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Relating To Insulation (AREA)
  • Locating Faults (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

一种变压器局部放电超声源重建的定位方法及装置,变压器局部放电的脉冲电流激发洛伦兹力超声信号,在变压器内部安装由多个超声传感器组成的传感器阵列,记录超声传感器接收到的超声信号,根据变压器内声波的传播模型,采用时间反转的逆源分析方法进行超声源的反演,进行局部放电点的定位,同时获取放电量信息,并估计放电量大小。

Description

一种变压器局部放电超声源重建的定位方法及装置
技术领域
本发明涉及一种变压器局部放电检测的方法,特别涉及一种变压器局部放电的超声检测方法及其装置。
背景技术
目前的大型电力变压器多为油浸式电力变压器,其绝缘结构主要由油、纸、纸板和其它固体绝缘等构成的固体-油绝缘结构。虽然在设计上具有足够的电气强度和优良的机械性能,但是制造过程中的偶然因素会造成一些先天性局部缺陷,如气泡、裂缝、悬浮导电质点和电极毛刺等。正是这些缺陷会造成绝缘体内部或表面出现某些区域电场强度高于平均电场强度,当这些区域的击穿场强低于平均击穿场强时,将会首先发生放电、而其它区域仍保持绝缘特性,从而形成局部放电。局部放电是变压器长期运行中绝缘劣化的一个重要征兆。电气设备的绝缘结构如果长时间发生持续的局部放电,绝缘介电性能可能会严重受损,如果局部放电故障一直未被发现和处理,最终可能导致电气设备发生故障,引起局部甚至全部地区的灾难性停电事故。
解决变压器故障问题最为有效的方法就是变压器局部放电在线检测,局部放电测量是监测变压器绝缘状况的重要诊断工具。研究发现,在绝缘材料内部局部放电过程中,通常伴随着很多物理现象,如电磁现象(电脉冲、电磁辐射等)和非电磁现象(气体、压力变化、光、热等)。利用放电过程产生的电信息或非电信息,可以对局部放电源进行定位,提供局部放电的凭据和类型。根据测量信息的不同,可以将局部放电检测方法分为电测法和非电测法两大类。
超声波法是一种非常重要的非电测类局部放电检测方法,主要是通过超声波传感器接收电气设备内部局部放电产生的超声波。超声检测法已经被广泛应用于地震监测、无损检测和全球定位系统等领域。在电力方面,超声波法是目前变压器绝缘检测领域的重要发展方向,它具有抗电气干扰能力强、定位方便等特点。关于超声波法局部放电检测的研究,
按超声探头位置,超声波法可分为内置式和外置式两类。内置式传感器接收到的波形幅值比外置式大很多,利于检测小放电量工况。由于变压器内部绝缘结构复杂,导致局部放电声信号的传播情况也无法准确预知,等值声速难以精确计算。此外,由于传感器接收的复合声波前沿不能准确确定,加之定位算法不够完善,使得局部放电信号提取与故障点定位精度受到影响。目前,信号前沿确定方法、定位算法是超声定位研究的热点和发展方向。国内外有关学者对于如何有效解决上述问题进行了深入研究,并提出了一系列新的理论和方法。
综合国内外对超声波法的研究,可分为基于声波时差和基于声场的两类:
基于声波时差的方法居多,如三角定位法、双曲面定位法、V型曲线法、顺序定位法和模式识别法等方法。三角定位法以局部放电电信号作为触发基准,通过三个(或以上)传感器测量超声波信号与电信号的时间差,以此作为超声信号的传播时间,并建立球面方程求解局部放电点。双曲面定位法,通过四个(或以上)超声传感器,以最先接收到的超声脉冲为时间基准,测量其他超声信号与标准信号的时延,由此建立双曲面方程组求解局部放电点。在上述方法中,通常假定声速为常数,实际上由于变压器结构复杂,各种声介质对声波及声速的影响都不一样,因此定位误差很大。尽管研究人员在基于最小二乘原理方法上进行了改进,一定程度上降低了传播路径对声速的影响,提高了定位精度,但改进方法仍然简化认为各路超声波信号声速相等,没有从根本上解决问题。在定位算法上的改进,还包括混沌遗传算法和改进的粒子群优化算法等。总体而言,这类方法的共性之处都是利用声波的时延性,由于变压器内部结构复杂,各种声介质对声速影响不同,超声波传播时发生多次折、反射,再加上现场各种噪声的影响,算法不尽完善等,实际的时延估计不能准确的得到用于定位的时间信息,给定位造成一定的误差,甚至难以成功定位。
基于声场类的定位方法还不多见。目前比较有特色的工作是声压数据库查表法。该方法的主要原理是:直接将变压器的局部放电点设定为超声波发射源,并进行声场分析计算,求出绕组上任何位置的点声源在超声波传感器阵列上激发的声压。同时,建立一个声压数据库,将超声波传感器阵列检测到的超声波信号声压和数据库里的数据进行比对,得出放电点的空间几何位置,从而实现局部放电的精确定位。这种方法的优势在于可以进行多点定位,不足之处在于,为构建数据库,需要对多源、任何位置、各种放电量等参数进行交叉离散。为实现高精度的定位,每个参数的采样率多,整体而言需要海量的样本量,构建过程所需时间难以承受,甚至是不可能的。
发明内容
为了克服上述现有方法的不足,本发明提出一种变压器局部放电超声源重建的定位方法。
本发明利用声源重建的方法进行超声源定位。变压器局部放电产生脉冲电流,脉冲电流产生脉冲磁场,脉冲电流和脉冲磁场二者相互作用,产生洛伦兹力,激发超声信号,变压器的放电点为超声源点。超声信号在变压器内部从放电点传播开来,由安装在变压器内部的多个超声传感器组成的传感器阵列接收超声信号,根据变压器内声波的传播模型,采用时间反转的逆源分析,反演重建声源,进行局部放电点的定位,同时获取放电量信息,估计放电量大小。
本发明变压器局部放电超声源重建的定位方法包含两个步骤:
步骤1、布置在变压器内部的超声传感器阵列获取超声传感器接收到的声场信息数据;
步骤2、根据变压器内部的结构和介质特性,建立变压器内部的声场的数学物理模型,利用时间反转的逆源方法,由超声传感器阵列接收到的超声信号,反演重建声源,从而估计超声源的位置和超声源强度,实现变压器局部放电的定位和放电量估计。
本发明利用时间反转的声源重建方法,是利用声场的互易原理实现声源的定位。超声源、变压器声传播模型、超声传感器阵列是定位的三要素。变压器局部放电时,根据法拉第电磁感定律,局部放电过程中产生的脉冲电流产生脉冲磁场,二者相互作用产生洛伦兹力,激发超声。超声从变压器的放电点,即超声源发射出来,经过变压器内部,传输至超声传感器接收阵列,由超声传感器接收阵列记录声压信号。超声传感器接收阵列接收到的声压信号重新发射出去,形成逆发射,再经过变压器内部的重新传播之后,在原超声源位置形成聚焦,据此可以实现声源的定位。聚焦形成的声源的强度大小反映了变压器的放电点放电量大小,根据聚焦形成声源的强度,再通过校准,可估算放电量。
超声传感器接收阵列记录到的声压的逆发射过程是利用数值计算方法实现。在变压器内的介质中,声压的传播特性可以利用时域的声波方程来描述, 为声压,t为时间,为变压器局部放电的脉冲电流,为脉冲电流产生的脉冲磁场,c为变压器内介质中的声速,为脉冲电流与脉冲磁场产生的洛伦兹力,作为局部放电激发的超声源。
首先,计算机记录了分布于变压器内四周的超声传感器接收到的声压信号。
然后,利用时间反转法由各个超声传感器接收到的超声信号计算反向声场,方法如下:
根据变压器内的空间结构、变压器内各种介质,如油、绕组、铁芯、外壳的材料特性,通过理论和实验确定求解声波方程所需要的传播速度、边界条件,即变压器的声传播模型。反演声场可以通过正向声场求解:
式中,T0为传感器信号的记录时间,为正向声场的时间反转, 为反演场的格林函数,为源点位置,为声场检测位置,c为变压器内介质中的声速。利用闭合面上的正向声压信号求解闭合面内部的反演声场,在本发明中,超声源的重建是利用分布在变压器内四周的超声传感器测得的声压信号逆源重建变压器内声场。
由记录的声压信号计算反演声场,实现声场的聚集和声源的重建。根据聚焦声场的位置和强度确定放电点位置和估算放电量大小。
应用本发明超声源重建的变压器局部放电的定位方法的装置包含超声传感器阵列、滤波放大器、数据采集器和计算机。超声传感器接收阵列由多个内置于变压器中的超声传感器组成,每个超声传感器与一个滤波放大器相连,超声传感器接收到的超声信号经过带通滤波和放大,滤除不相关的干扰信号,提取有用的通带内的超声信号,并放大微弱的超声信号。多个滤波放大器与数据采集器相连。数据采集器与计算机相连。数据采集器将信号模数转换并传输给计算机,计算机利用数值计算方法实现超声源的重建。
本发明利用声场逆源分析方法考虑超声在变压器条件下的传播过程,解决了声波时差法由于声射线多次反射和折射引起的定位精度低等问题,同时克服了声压数据库查表法数据量大、定位效率低的不足,并可以同时确定放电位置和放电量信息。
附图说明
图1本发明变压器超声源定位系统框图;
图2本发明时间反转法示意图;
图3本发明时间反转法重建声源框图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步说明。
本发明利用声源重建的方法进行超声源定位。变压器局部放电的脉冲电流激发超声信号,变压器的放电点为超声源点。超声信号在变压器内部从放电点传播开来。由安装在变压器内部的多个超声传感器组成的传感器阵列记录超声传感器接收到的超声信号,根据声压波动方程和变压器声传播模型,利用时间反转法由各个超声传感器接收到的超声信号计算反向声场。再将反向声场叠加得到聚焦声场,实现声源的重建。根据聚焦声场的位置和强度确定放电点位置和估算放电量大小。
利用本发明的变压器局部放电超声源定位装置如图1所示。所述的装置包含变压器、超声传感器阵列、滤波放大器、数据采集器和计算机。变压器内部主要由绕组和变压器油组成。变压器内部四周安装有由多个超声传感器组成的超声接收阵列。所述的超声传感器的信号通过电缆连接到滤波放大器的输入端,滤波放大器的输出连接到数据采集器,数据采集器和计算机连接,将所采集的数据传输到计算机。变压器内部局放放电时,激发超声信号,超声信号在变压器内部传播,到达变压器四周的超声传感器,超声传感器将接收到的声压信号转换为电压信号,并经过滤波放大器放大滤波之后,输入到数据采集器,经数据采集器进行模拟-数字转换,转化的数字信号传输到计算机存储和处理。
所述的时间反转法如图2所示。变压器局部放电时,放电点为超声源点,局部放电激发的超声信号从放电点开始传播,经过变压器的内部,到达变压器四周的超声传感器阵列,超声传感器阵列接收并记录到所述的超声信号。这是超声信号传播的正过程,如图2a所示。本发明利用时间反转法逆发射重建超声源,如图2b所示,将各个传感器阵列记录的超声信号再重新发射回去,根据声场的互易定理,重新发射形成的声波会在原声源位置聚焦,实现声场重建,据此可以确定超声源位置,从而确定放电点位置,并根据聚焦的声源强度的大小,估算放电量大小。
所述的利用时间反转法的逆发射实现声场重建是由计算机采用数值计算的方法实现的。如图3所示,变压器放电点激发超声,超声信号在变压器内部传播之后到达超声传感器阵列,接收到的信号经过放大滤波之后经数据采集器转换成数字信号存储到计算机中,计算机利用数值计算的方法实现逆发射,如图3的虚线框所示,首先,记录超声传感器阵列中各个传感器接收到超声信号P1,P2,...Pn。然后,根据声压波动方程和变压器声传播模型,利用时间反转法由各个超声传感器接收到的信号P1,P2,...Pn计算反向声场P’1,P’2,...P’n。再将反向声场叠加得到聚焦声场,实现声源的重建。根据聚焦声场的位置和强度确定放电点位置和估计放电量大小。

Claims (3)

1.一种变压器局部放电超声源重建的定位方法,其特征在于所述的方法为:变压器局部放电产生脉冲电流,脉冲电流产生脉冲磁场,脉冲电流和脉冲磁场的相互作用,产生洛伦兹力,激发超声信号;变压器的放电点为超声源点;超声信号在变压器内部从放电点传播开来,由安装在变压器内部的多个超声传感器组成的传感器阵列接收所述的超声信号;根据变压器内声波的传播模型,采用时间反转的逆源分析,由超声传感器阵列接收到的超声信号反演重建声源,进行局部放电点的定位,同时获取放电量信息,并估计放电量大小。
2.根据权利要求1所述的变压器局部放电超声源重建的定位方法,其特征在于所述的变压器局部放电定位方法包含两个步骤:
步骤1、由布置在变压器内部的超声传感器阵列获取超声传感器接收到的声场信息数据;
步骤2、根据变压器内部的结构和介质特性,建立变压器内部的声场的数学物理模型,利用时间反转的逆源方法,由超声传感器阵列接收到的超声信号,反演重建声源,从而估计超声源的位置和超声源强度,实现变压器局部放电的定位和放电量估计。
3.根据权利要求2所述的变压器局部放电超声源重建的定位方法,其特征在于所述的步骤2中,利用时间反转的逆源分析,由超声传感器阵列接收到的超声信号反演重建声源的方法为:利用声场的互易原理实现声源的检测与定位;超声源、变压器声传播模型、超声传感器阵列是定位的三要素;局部放电时,超声从声源发射出来,经过变压器内的传输,到达超声传感器接收阵列,从而记录到声压信号;超声传感器接收阵列接收到的声压信号重新发射出去,形成逆发射,经过变压器内部的重新传播之后,在原来声源的位置形成聚焦,实现声源的定位;洛伦兹力激发的超声携带放电点信息;利用超声传感器接收到的超声信号逆发射,聚焦形成的声源的强度大小反映放电量大小;根据聚焦形成声源的强度,再通过校准,实现放电量的估计;超声传感器接收阵列记录到的声压的逆发射过程是在计算机中利用数值计算方法实现的:首先,记录超声传感器阵列中各个传感器接收到超声信号P1,P2,...Pn;然后,根据声压波动方程和变压器声传播模型,利用时间反转法由各个超声传感器接收到的信号P1,P2,...Pn计算反向声场P’1,P’2,...P’n;再将反向声场叠加得到聚焦声场,实现声源的重建;根据聚焦声场的位置和强度确定放电点位置和估计放电量大小。
CN201310060568.8A 2013-02-26 2013-02-26 一种变压器局部放电超声源重建的定位方法及装置 Active CN103149513B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310060568.8A CN103149513B (zh) 2013-02-26 2013-02-26 一种变压器局部放电超声源重建的定位方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310060568.8A CN103149513B (zh) 2013-02-26 2013-02-26 一种变压器局部放电超声源重建的定位方法及装置

Publications (2)

Publication Number Publication Date
CN103149513A CN103149513A (zh) 2013-06-12
CN103149513B true CN103149513B (zh) 2015-06-24

Family

ID=48547704

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310060568.8A Active CN103149513B (zh) 2013-02-26 2013-02-26 一种变压器局部放电超声源重建的定位方法及装置

Country Status (1)

Country Link
CN (1) CN103149513B (zh)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969560B (zh) * 2014-05-28 2016-10-26 北达正视(北京)科技有限公司 一种用于局部放电故障检测的可视化超声波检测系统
CN103969561B (zh) * 2014-05-28 2017-01-25 上海艾飞能源科技有限公司 一种用于局部放电故障检测的可视化超声波检测方法
CN105334480B (zh) * 2014-08-11 2018-09-18 国家电网公司 四阵元的传感器阵列空间位置的确定方法和装置
CN104950233B (zh) * 2015-06-18 2017-09-01 国家电网公司 基于场源逆问题反演的开关柜局部放电源定位与重构方法
CN105548824B (zh) * 2015-12-02 2018-05-22 华北电力大学(保定) 一种电气设备局放源定位装置及定位方法
CN105719003A (zh) * 2016-01-13 2016-06-29 中国南方电网有限责任公司超高压输电公司 一种基于量子遗传算法的换流变压器局放超声定位方法
FI127711B (en) * 2016-05-13 2018-12-31 Altum Tech Oy METHOD FOR CLEANING THE APPLIANCE
JP6815764B2 (ja) * 2016-06-28 2021-01-20 株式会社日立製作所 溶接監視システム
CN107390100B (zh) * 2017-08-03 2019-07-19 四川大学 一种基于时间反演的电力电缆局部放电定位方法
CN109948791B (zh) * 2017-12-21 2022-03-29 河北科技大学 利用遗传算法优化bp神经网络的方法及其在定位中的应用
CN108152687B (zh) * 2017-12-25 2019-12-31 云南电网有限责任公司昆明供电局 一种利用超声反演的电力变压器局部放电定位方法
CN108387817B (zh) * 2018-01-23 2020-10-27 中国科学院电工研究所 一种基于注入电流式多场耦合的接地网故障诊断方法及装置
CN109188359B (zh) * 2018-09-13 2020-07-21 中国气象科学研究院 闪电定位方法及装置
CN108983032B (zh) * 2018-09-17 2021-02-12 华北电力大学(保定) 一种基于时间反转的电力电缆开路故障定位方法
CN109901031B (zh) * 2019-02-27 2021-06-11 西安电子科技大学 用于局部放电信号的信噪分离方法、信息数据处理终端
CN110703150B (zh) * 2019-10-10 2021-08-03 重庆大学 基于声振耦合信号的变压器内部故障检测与定位方法
CN111157845B (zh) * 2020-01-07 2021-10-26 国网辽宁省电力有限公司朝阳供电公司 一种适用于配电网的基于时域反演技术的故障定位方法
CN111413649B (zh) * 2020-04-23 2021-07-06 中国科学技术大学 基于近场宽带波束形成的大型电抗器故障检测方法及系统
CN112083297B (zh) * 2020-08-20 2021-08-13 西安交通大学 一种基于声压的抗干扰的变压器局部放电定位方法及装置
CN112466627B (zh) * 2020-11-11 2021-12-24 西安西电变压器有限责任公司 一种变压器油箱的设计方法、系统及变压器油箱
CN113051829B (zh) * 2021-03-31 2023-06-02 西南大学 一种应用空间分析理论改进的变压器Duval Pentagon1故障诊断方法
CN113238125A (zh) * 2021-04-30 2021-08-10 中国矿业大学 一种磁矢量分布场重构与电缆局部放电定位方法
CN114325275A (zh) * 2022-02-07 2022-04-12 国网宁夏电力有限公司电力科学研究院 利用信道模型实现路径反演的站域局放定位方法及系统
CN118584399A (zh) * 2024-08-02 2024-09-03 红光电气集团有限公司 一种10kv变压器台区成套设备局部放电位置的定位方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1888924A (zh) * 2006-07-27 2007-01-03 华北电力大学 变压器局部放电超宽带传感器阵列定位系统及其方法
CN101169390A (zh) * 2007-10-12 2008-04-30 南京航空航天大学 一种工程结构损伤主动监测中Lamb波时间反转聚焦方法
CN102621223A (zh) * 2011-01-31 2012-08-01 中国科学院声学研究所 一种基于声波时间反转法的超声扫描与检测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1888924A (zh) * 2006-07-27 2007-01-03 华北电力大学 变压器局部放电超宽带传感器阵列定位系统及其方法
CN101169390A (zh) * 2007-10-12 2008-04-30 南京航空航天大学 一种工程结构损伤主动监测中Lamb波时间反转聚焦方法
CN102621223A (zh) * 2011-01-31 2012-08-01 中国科学院声学研究所 一种基于声波时间反转法的超声扫描与检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
L.E.Lundgaard.Partial Discharge-Part XIV:Acoustic Partial Discharge Detection-Practical Application.《IEEE Electrical Insulation Magazine》.1992,第8卷(第5期),第34-43页. *
基于超声技术的变压器内部局部放电定位研究;刘蓉;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20100715(第07期);第C042-57页 *

Also Published As

Publication number Publication date
CN103149513A (zh) 2013-06-12

Similar Documents

Publication Publication Date Title
CN103149513B (zh) 一种变压器局部放电超声源重建的定位方法及装置
Mondal et al. Partial discharge localization in a power transformer: Methods, trends, and future research
Zheng et al. Location of PDs inside transformer windings using UHF methods
CN103344934B (zh) 电力变压器局部放电超声传感器的检测校验方法及系统
CN203249989U (zh) 变压器局部放电在线监测系统
CN110703150B (zh) 基于声振耦合信号的变压器内部故障检测与定位方法
CN102288625A (zh) 一种渠道混凝土衬砌无损检测系统
CN103424675A (zh) 超高频天线阵列局部放电检测系统
CN101666850A (zh) 一种电力电缆绝缘缺陷的声电检测装置
CN103076587B (zh) 传感器接收性能的标定方法、系统及超声波脉冲发生器
CN105301114A (zh) 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法
CN104614648A (zh) 一种电声联合直流局部放电检测装置
CN103558525A (zh) 一种开关柜局部放电检测方法
CN108152687B (zh) 一种利用超声反演的电力变压器局部放电定位方法
CN103412053B (zh) 一种基于双声发射传感阵列和波束形成的无需波速的声发射源定位方法
CN107192930A (zh) 一种变压器局部放电定位方法
CN104034801B (zh) 基于合成时反的结构损伤迭代聚焦成像监测方法
CN108037410A (zh) 一种基于可控响应功率的局部放电超声波定位方法与装置
CN105486924B (zh) 基于磁声电效应的非接触式导体电导率测量方法
US20200378825A1 (en) Transducer transfer impedance calibration device based on spatial frequency domai smoothing technology
CN111458566A (zh) 一种储能材料电导率的非接触式检测方法及系统
CN118483529A (zh) 基于声-电联合的变压器局部放电源空间定位方法及系统
CN112816835B (zh) 基于电声联合检测信号传播时延补偿的局部放电定位方法
CN114839269A (zh) Gis固体绝缘件内部缺陷应力无损检测方法及装置
CN105974281A (zh) 一种用于变压器局部放电定位的检测系统和检测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant