CN114636900A - 一种配电网多重故障诊断方法及系统 - Google Patents

一种配电网多重故障诊断方法及系统 Download PDF

Info

Publication number
CN114636900A
CN114636900A CN202210506365.6A CN202210506365A CN114636900A CN 114636900 A CN114636900 A CN 114636900A CN 202210506365 A CN202210506365 A CN 202210506365A CN 114636900 A CN114636900 A CN 114636900A
Authority
CN
China
Prior art keywords
fault
power distribution
distribution network
line
branch
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
CN202210506365.6A
Other languages
English (en)
Other versions
CN114636900B (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.)
Dongguan Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Dongguan Power Supply Bureau of Guangdong Power Grid 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 Dongguan Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Dongguan Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority to CN202210506365.6A priority Critical patent/CN114636900B/zh
Publication of CN114636900A publication Critical patent/CN114636900A/zh
Application granted granted Critical
Publication of CN114636900B publication Critical patent/CN114636900B/zh
Priority to PCT/CN2022/133889 priority patent/WO2023216553A1/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/086Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • 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
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage or fault management, e.g. fault detection or location

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Locating Faults (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种配电网多重故障诊断方法及系统,所述方法包括:S1、利用MATLAB平台对配电网线路进行短路故障分析,得出基于电气量的故障信息决策表;S2、利用Simulink平台对故障信息决策表进行建模仿真;S3、采用小波变换法对神经网络的输出训练数据进行去噪、采集处理,形成相关故障信息决策表作为神经网络的训练样本;S4、采用改进的Artificial Tree智能优化算法优化神经网络的权值和阈值,选择部分数据作为故障数据将训练好的神经网络用来检测故障,不仅能提高电网故障诊断准确性与可靠性,还可在故障信息存在误动时对故障进行检测,从而达到定位电网故障区域的效果,有效增强了电网故障检测的能力。

Description

一种配电网多重故障诊断方法及系统
技术领域
本发明涉及配电网故障诊断技术领域,具体涉及一种配电网多重故障诊断方法及系统。
背景技术
当今时代,电力能源已成为人类生活中不可缺少的重要资源,在电网运行时一旦出现故障将会给社会生产和人民生活带来巨大的损失,因此当电网发生故障时,应采取必要的措施,快速准确地定位故障发生的区域,找出具体的故障线路,提高电网故障后恢复供电的速度和检测精度。
现有的配电网故障检测系统还存在一下缺陷:
(1)实际配电网中的采集信号装置和其他保护装置容易受到人为或者各种电磁设备之间产生的相互影响,使得配电网中存在着很多不确定因素,现有的故障检测方法针对线路动态变化时,无法及时对故障定位检测,在故障信息存在误动信息的情况下,容错能力较差。
(2)目前针对故障检测的拓扑结构规划完成后,当遇到实际电网的结构和相关保护器件发生改变之后,将很难对网络的拓扑结构进行改变,而是需要重新进行设计,在实际操作中会变得非常繁琐。
发明内容
本发明的目的在于提供一种配电网多重故障诊断方法及系统,以解决现有技术中针对线路动态变化时,无法及时对故障定位检测,在故障信息存在误动信息的情况下,容错能力较差,拓扑结构规划繁琐的技术问题。
为解决上述技术问题,本发明具体提供下述技术方案:
一种配电网多重故障诊断方法,包括以下步骤:
步骤S1、利用MATLAB平台对所述配电网线路进行短路故障分析,得出基于电气量的故障信息决策表;
步骤S2、利用Simulink平台对所述基于电气量的故障信息决策表进行建模仿真,获取配电网线路神经网络的输出训练数据;
步骤S3、采用小波变换法对所述神经网络的输出训练数据进行去噪、采集处理,形成相关故障信息决策表作为神经网络的训练样本;
步骤S4、采用改进的Artificial Tree 智能优化算法优化神经网络的权值和阈值,通过迭代达到相关误差精度的要求,选择部分数据作为故障数据并利用训练好的神经网络来检测故障。
作为本发明的一种优选方案,所述步骤S1中,所述MATLAB平台采集配电网线路发生小电流接地短路故障后产生的零序电流、零序功率和零序导纳信号,基于这三种信息形成基于电气量的故障信息决策表。
作为本发明的一种优选方案,依据所述零序电流、零序功率和零序导纳信号用离散采样时间和连续采样时间相混合的采样模式进行故障分析,具体步骤如下:
步骤S201、在每条线路两个终端中的一个终端上采样并存储正负极零序电流;
步骤S202、选择固定长度的时间窗内的采样电流作为线路参考电流,对动态采样电流与线路参考电流进行比较,使用改进Pearson方法实时求解二者的相关系数大小;
步骤S203、设定相关系数阈值,判断相关系数与所述设定阈值的大小,并输出故障所在线路及故障类型。
作为本发明的一种优选方案,所述步骤S1中,所述改进Pearson方法采用线路相关系数的累加值进行故障检测的相关系数分析,所述改进Pearson方法具体步骤如下:
步骤S2021、定义具有固定长度的时间窗内可移动窗口值
Figure 852123DEST_PATH_IMAGE001
Figure 395362DEST_PATH_IMAGE002
,所述
Figure 175100DEST_PATH_IMAGE003
Figure 694943DEST_PATH_IMAGE004
窗口的末端固定时刻分别为
Figure 509315DEST_PATH_IMAGE005
Figure 708609DEST_PATH_IMAGE006
,将所述
Figure 975642DEST_PATH_IMAGE007
中的采样电流与稳态参考电流进行比较,检测直流微电网的故障,所述
Figure 174542DEST_PATH_IMAGE008
Figure 968054DEST_PATH_IMAGE004
的关系式为:
Figure 587517DEST_PATH_IMAGE009
其中
Figure 341846DEST_PATH_IMAGE010
为线路电流的采样频率,
Figure 78858DEST_PATH_IMAGE011
为采样点的数量,
Figure 602243DEST_PATH_IMAGE012
从故障发生时刻
Figure 156721DEST_PATH_IMAGE013
到故障被检测到时刻
Figure 398347DEST_PATH_IMAGE014
之间的时长;
步骤S2022、采用线路相关系数的累加值进行故障检测的相关系数分析,将稳态参考电流曲线和采样电流曲线的第
Figure 673470DEST_PATH_IMAGE015
个对比值可分别表示为
Figure 316941DEST_PATH_IMAGE016
Figure 652108DEST_PATH_IMAGE017
Figure 381029DEST_PATH_IMAGE018
其中
Figure 210576DEST_PATH_IMAGE019
为第
Figure 708554DEST_PATH_IMAGE020
个暂态线路电流,
Figure 106299DEST_PATH_IMAGE021
为第
Figure 56938DEST_PATH_IMAGE020
个采样线路电流;
步骤S2023、引入调整因子
Figure 939443DEST_PATH_IMAGE022
优化Pearson相关系数调整因子
Figure 291927DEST_PATH_IMAGE023
表达式为:
Figure 234475DEST_PATH_IMAGE024
其中,
Figure 672410DEST_PATH_IMAGE025
为数学常数,表示为自然对数的底数,
Figure 217661DEST_PATH_IMAGE026
为两条曲线所对比点的总数;
步骤S2023、依据所述调整因子
Figure 690230DEST_PATH_IMAGE027
优化相关系数
Figure 803680DEST_PATH_IMAGE028
,所述相关系数
Figure 728911DEST_PATH_IMAGE028
的表达式为:
Figure 953218DEST_PATH_IMAGE029
Figure 14715DEST_PATH_IMAGE030
作为本发明的一种优选方案,所述步骤S2中,依据所述故障分析获取故障位置和故障电阻通过Simulink平台构建模型输出训练数据,采用小波变换法处理训练数据。
作为本发明的一种优选方案,所述小波变换法处理训练数据具体步骤为:
步骤S301、对含噪声的训练数据信号
Figure 925165DEST_PATH_IMAGE031
作小波变换,获取一组小波分解系数;
步骤S302、采样软阈值函数对得到的小波分解系数进行阈值处理,得到一组估计小波系数,所述软阈值函数
Figure 337692DEST_PATH_IMAGE032
表达式为:
Figure 365690DEST_PATH_IMAGE033
Figure 547273DEST_PATH_IMAGE034
其中,
Figure 736946DEST_PATH_IMAGE035
为符号函数,
Figure 761403DEST_PATH_IMAGE036
为小波系数,
Figure 327513DEST_PATH_IMAGE037
为阈值,
Figure 629182DEST_PATH_IMAGE038
为噪声强度,
Figure 989756DEST_PATH_IMAGE039
为信号长度,
Figure 376875DEST_PATH_IMAGE040
为幅值因子,
Figure 746676DEST_PATH_IMAGE041
Figure 239933DEST_PATH_IMAGE042
为调节因子;
步骤S303、通过调节
Figure 36988DEST_PATH_IMAGE043
使得在阈值点
Figure 645824DEST_PATH_IMAGE044
处连续且降低对原信号产生的震荡误差,提取去噪后的信号。
作为本发明的一种优选方案,在步骤S4中,将所述震荡误差的均方误差采用Artificial Tree 算法优化获取线路神经网络,所述Artificial Tree 算法优化步骤如下:
步骤S401、以所述震荡误差的均方误差MSE作为目标函数,将收集到的数据进行整理,准备好训练样本和测试样本;
步骤S402、设置Artificial Tree 算法的相关参数,确定分支个数
Figure 943950DEST_PATH_IMAGE045
、空间维数
Figure 954631DEST_PATH_IMAGE046
、迭代次数
Figure 922587DEST_PATH_IMAGE047
,并对分支进行初始化操作,生成初始分支;
步骤S402、计算每个分支所对应的函数值,根据函数值选出最好分支
Figure 18719DEST_PATH_IMAGE048
和最好分支所对应的函数值
Figure 730323DEST_PATH_IMAGE049
步骤S403、判断是否满足误差精度的要求或者是否达到了最大迭代次数,如果都不满足要求则准备进入IAT算法的迭代寻优过程;
步骤S404、如果搜索次数达到了搜索次数的上限,将进行随机操作生成新的分支,将新产生的分支与旧的分支进行比较,如果新的分支优于旧的分支则新的分支将取代旧分支,否则仍以旧的分支为最优分支;
步骤S405、再次判断最优分支是否满足误差精度或者当前迭代次数是否达到了设定的最大迭代次数
Figure 595511DEST_PATH_IMAGE047
作为本发明的一种优选方案,所述Artificial Tree 算法的相关参数依据配电网线路保护装置的动作信息及对应发生故障的区域设定,形成相应的神经网络训练样本。
作为本发明的一种优选方案,所述配电网线路保护装置的动作信息主要采集基于开关量的电网信号,用0、1对电网开关动作进行编号。
一种配电网多重故障诊断方法的系统,包括:电路采样模块、线路保护装置、故障分析模块、仿真模块以及检测输出模块
所述电路采样模块,用于采集和预处理导入到所述MATLAB平台和Simulink平台的各项数据;
所述线路保护装置,用于配电网发生故障时会进行跳闸或者报警等相应的保护动作,根据保护装置不同的保护动作,判断故障发生的大致区间,并且将电网的区间进行编号;
所述故障分析模块,用于分析故障发生位置及故障电阻阻值,判断故障大致区域;
所述仿真模块,利用Simulink平台通过Artificial Tree智能优化算法获取故障检测数据,显示仿真结果。
所述检测输出模块,依据仿真结果选用部分数据检测线路神经网络对电网故障检测的准确率分析适应性报告。
本发明与现有技术相比较具有如下有益效果:
本发明利用MATLAB平台用改进后的小波阈值去噪方法去除信号中的噪声源建立电网故障模型,提高信号采集准确度,运行仿真软件进行仿真,并建立基于开关量的电网故障模型得到电网故障信息决策表,训练神经网络,检测故障发生的区间,对复杂电网采取分而治之的策略,使得线路拓扑结构能够灵活多变,这样不仅能提高电网故障诊断的准确性与可靠性,还可以在故障信息存在误动信息的情况下对故障进行检测,从而达到电网故障区域定位的效果,有效的增强了电网故障检测的能力。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。
图1为本发明实施例提供的配电网多重故障诊断方法流程图;
图2为本发明实施例提供的动态故障检测时间轴示意图;
图3为本发明实施例提供的故障检测系统结构框图。
图中的标号分别表示如下:
1-电路采样模块;2-线路保护装置;3-故障分析模块;4-仿真模块;5-检测输出模块。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1-3所示,本发明提供了一种配电网多重故障诊断方法,包括以下步骤:
步骤S1、利用MATLAB平台对所述配电网线路进行短路故障分析,得出基于电气量的故障信息决策表;
步骤S2、利用Simulink平台对所述基于电气量的故障信息决策表进行建模仿真,获取配电网线路神经网络的输出训练数据;
步骤S3、采用小波变换法对所述神经网络的输出训练数据进行去噪、采集处理,形成相关故障信息决策表作为神经网络的训练样本;
步骤S4、采用改进的Artificial Tree 智能优化算法优化神经网络的权值和阈值,通过迭代达到相关误差精度的要求,选择部分数据作为故障数据并利用训练好的神经网络来检测故障。
本实施例中,根据配电网线路数据建立电网故障模型,收集所需要的神经网络输入训练的数据,在收集信息的同时采用小波变换理论对收集的数据进行去噪处理,去除干扰数据准确性的噪声源,然后将去噪后的数据进行采集处理,形成相关故障信息决策表作为神经网络的训练样本,然后应用Artificial Tree智能优化算法来优化神经网络的权值和阈值,通过迭代达到相关误差精度的要求,选择部分数据作为故障数据将训练好的神经网络用来检测故障。
所述步骤S1中,所述MATLAB平台采集配电网线路发生小电流接地短路故障后产生的零序电流、零序功率和零序导纳信号,基于这三种信息形成基于电气量的故障信息决策表。
本实施例中,在应用神经网络解决电网故障检测问题时,将保护装置的动作信息作为神经网络的输入值,将对应发生故障的区域作为输出值,形成相应的神经网络训练样本,从而用来建立检测模型。
依据所述零序电流、零序功率和零序导纳信号用离散采样时间和连续采样时间相混合的采样模式进行故障分析,具体步骤如下:
步骤S201、在每条线路两个终端中的一个终端上采样并存储正负极零序电流;
步骤S202、选择固定长度的时间窗内的采样电流作为线路参考电流,对动态采样电流与线路参考电流进行比较,使用改进Pearson方法实时求解二者的相关系数大小;
步骤S203、设定相关系数阈值,判断相关系数与所述设定阈值的大小,并输出故障所在线路及故障类型。
本实施例中,如果正负两极中的一个极的电流相关程度低于阈值,且该线路中两个极点的电流和不为零,则认为该线路发生极地短路故障;如果线路正负极的电流相关程度均低于阈值,且该线路中两个极点的电流和近似为零,则认为该线路发生极间短路故障;否则,认为该线路没有发生故障。
本实施例中,电网故障保护由故障检测、故障隔离和故障定位三部分组成,如图2所示:假定在
Figure 360467DEST_PATH_IMAGE050
时刻发生短路故障,故障线路的电流迅速增大。故障检测设备会在之后的某时刻检测到故障,定义该时刻为
Figure 943895DEST_PATH_IMAGE051
,从
Figure 193611DEST_PATH_IMAGE051
时刻开始,故障隔离和故障定位同步进行:一方面故障隔离设备开始动作,并于之后的
Figure 178884DEST_PATH_IMAGE052
时刻成功开断故障线路;另一方面故障定位设备开始预测故障的发生位置,并于之后的
Figure 488643DEST_PATH_IMAGE053
时刻给出预测的故障发生位置,分析故障保护中各部分的时间顺序,可以发现:
一、故障隔离需要故障被检测到后才会动作,由于直流微电网故障隔离需要在几个毫秒内完成,因此系统对故障检测的快速性要求很高;
二、故障定位是在故障被检测到后,与故障隔离同步开始进行,且故障定位与故障隔离两个过程相互解耦,因此系统对故障定位的快速性要求不高;
三、由于故障隔离设备动作后线路电流变化与故障隔离设备的动作情况相关且难以预测,因此只有故障被检测到前的系统状态量是可供故障检测与故障定位的优质数据。
所述步骤S202中,所述改进Pearson方法采用线路相关系数的累加值进行故障检测的相关系数分析,所述改进Pearson方法具体步骤如下:
步骤S2021、定义具有固定长度的时间窗内可移动窗口值
Figure 559367DEST_PATH_IMAGE054
Figure 612774DEST_PATH_IMAGE055
,所述
Figure 203286DEST_PATH_IMAGE056
Figure 683946DEST_PATH_IMAGE002
窗口的末端固定时刻分别为
Figure 976387DEST_PATH_IMAGE057
Figure 833485DEST_PATH_IMAGE058
,将所述
Figure 527771DEST_PATH_IMAGE056
中的采样电流与稳态参考电流进行比较,检测直流微电网的故障,所述
Figure 179332DEST_PATH_IMAGE008
Figure 224649DEST_PATH_IMAGE004
的关系式为:
Figure 744492DEST_PATH_IMAGE059
其中
Figure 293285DEST_PATH_IMAGE060
为线路电流的采样频率,
Figure 381327DEST_PATH_IMAGE061
为采样点的数量,
Figure 648360DEST_PATH_IMAGE062
从故障发生时刻
Figure 847260DEST_PATH_IMAGE063
到故障被检测到时刻
Figure 876658DEST_PATH_IMAGE064
之间的时长;
步骤S2022、采用线路相关系数的累加值进行故障检测的相关系数分析,将稳态参考电流曲线和采样电流曲线的第
Figure 135601DEST_PATH_IMAGE015
个对比值可分别表示为
Figure 889931DEST_PATH_IMAGE065
Figure 626942DEST_PATH_IMAGE066
Figure 274961DEST_PATH_IMAGE067
其中
Figure 704806DEST_PATH_IMAGE068
为第
Figure 680852DEST_PATH_IMAGE015
个暂态线路电流,
Figure 221555DEST_PATH_IMAGE069
为第
Figure 865026DEST_PATH_IMAGE015
个采样线路电流;
步骤S2023、引入调整因子
Figure 200192DEST_PATH_IMAGE070
优化Pearson相关系数调整因子
Figure 555212DEST_PATH_IMAGE027
表达式为:
Figure 634027DEST_PATH_IMAGE071
其中,
Figure 132004DEST_PATH_IMAGE025
为数学常数,表示为自然对数的底数,
Figure 762706DEST_PATH_IMAGE026
为两条曲线所对比点的总数;
步骤S2023、依据所述调整因子
Figure 978923DEST_PATH_IMAGE070
优化相关系数
Figure 861429DEST_PATH_IMAGE072
,所述相关系数
Figure 964645DEST_PATH_IMAGE072
的表达式为:
Figure 907193DEST_PATH_IMAGE073
Figure 345128DEST_PATH_IMAGE074
本实施例中,在步骤S2022中,如果直接对采样电流和稳态计算电流进行对比,相关系数将随电流纹波的变化而变化,可能出现大量的误检测,为防止误检测的发生,采用线路的累加值进行故障检测的相关系数分析,消除由这种现象引起的干扰,导致的故障定位准确性不高的现象。
本实施例中,在故障检测中对采样电流与稳态参考电流进行比较,使用改进Pearson 方法实时求解二者的相关系数大小,当相关系数降低到设定阈值以下时,认为发生短路故障,并输出故障所在线路及故障类型;
本实施例中,在故障定位中采用基于遗传算法的迭代方法来实现所预测故障位置和故障阻抗的生成与更新,计算所预测故障位置和故障阻抗对应的故障电流曲线,并采用改进Pearson方法求解采样曲线与所计算曲线的相关系数大小,当相关系数升高到设定阈值以上时,认为所预测故障位置和故障阻抗足够准确,并输出预测的故障位置和故障阻抗。
所述步骤S2中,依据所述故障分析获取故障位置和故障电阻通过Simulink平台构建模型输出训练数据,采用小波变换法处理训练数据。
所述小波变换法处理训练数据具体步骤为:
步骤S301、对含噪声的训练数据信号
Figure 765745DEST_PATH_IMAGE075
作小波变换,获取一组小波分解系数;
步骤S302、采样软阈值函数对得到的小波分解系数进行阈值处理,得到一组估计小波系数,所述软阈值函数
Figure 238315DEST_PATH_IMAGE076
表达式为:
Figure 712283DEST_PATH_IMAGE077
Figure 637514DEST_PATH_IMAGE078
其中,
Figure 861822DEST_PATH_IMAGE079
为符号函数,
Figure 188898DEST_PATH_IMAGE036
为小波系数,
Figure 473249DEST_PATH_IMAGE044
为阈值,
Figure 10410DEST_PATH_IMAGE038
为噪声强度,
Figure 772829DEST_PATH_IMAGE080
为信号长度,
Figure 954412DEST_PATH_IMAGE081
为幅值因子,
Figure 409664DEST_PATH_IMAGE082
Figure 309487DEST_PATH_IMAGE083
为调节因子;
步骤S303、通过调节
Figure 875597DEST_PATH_IMAGE084
使得在阈值点
Figure 537785DEST_PATH_IMAGE085
处连续且降低对原信号产生的震荡误差,提取去噪后的信号。
本实施例中,噪声源引起的小波系数的值将会小于预先设置的临界阈值,而原始信号所引起的小波系数将会大于预先设置的临界阈值,利用这个特性,将小波系数小于临界阈值的噪声信号去除,将小波系数大于临界阈值的原始信号保留,然后进行小波逆变换对信号进行重构,达到去噪的目的。
本实施例中,通过调节
Figure 163939DEST_PATH_IMAGE086
Figure 551058DEST_PATH_IMAGE087
的值来增强阈值函数使用的灵活性,参数
Figure 655280DEST_PATH_IMAGE088
决定阈值的逼近程度,通过调节
Figure 811455DEST_PATH_IMAGE088
使得在阈值点
Figure 733143DEST_PATH_IMAGE089
处连续且降低对原信号产生的震荡误差,使得提取到的去噪后的原始信号更为准确,可信度更高,由此提高了电网故障检测中对于故障信号去噪的精度,从而提高了电网故障检测的准确性。
在步骤S4中,将所述震荡误差的均方误差采用Artificial Tree 算法优化获取线路神经网络,所述Artificial Tree 算法优化步骤如下:
步骤S401、以所述震荡误差的均方误差MSE作为目标函数,将收集到的数据进行整理,准备好训练样本和测试样本;
步骤S402、设置Artificial Tree 算法的相关参数,确定分支个数
Figure 341979DEST_PATH_IMAGE090
、空间维数
Figure 515471DEST_PATH_IMAGE091
、迭代次数
Figure 526153DEST_PATH_IMAGE092
,并对分支进行初始化操作,生成初始分支;
步骤S402、计算每个分支所对应的函数值,根据函数值选出最好分支
Figure 120207DEST_PATH_IMAGE048
和最好分支所对应的函数值
Figure 340973DEST_PATH_IMAGE093
步骤S403、判断是否满足误差精度的要求或者是否达到了最大迭代次数,如果都不满足要求则准备进入IAT算法的迭代寻优过程;
步骤S404、如果搜索次数达到了搜索次数的上限,将进行随机操作生成新的分支,将新产生的分支与旧的分支进行比较,如果新的分支优于旧的分支则新的分支将取代旧分支,否则仍以旧的分支为最优分支;
步骤S405、再次判断最优分支是否满足误差精度或者当前迭代次数是否达到了设定的最大迭代次数
Figure 786998DEST_PATH_IMAGE047
所述Artificial Tree 算法的相关参数依据配电网线路保护装置的动作信息及对应发生故障的区域设定,形成相应的神经网络训练样本。
所述配电网线路保护装置的动作信息主要采集基于开关量的电网信号,用0、1对电网开关动作进行编号。
一种配电网多重故障诊断方法的系统,包括:电路采样模块1、线路保护装置2、故障分析模块3、仿真模块4以及检测输出模块5
所述电路采样模块1,用于采集和预处理导入到所述MATLAB平台和Simulink平台的各项数据;
所述线路保护装置2,用于配电网发生故障时会进行跳闸或者报警等相应的保护动作,根据保护装置不同的保护动作,判断故障发生的大致区间,并且将电网的区间进行编号;
所述故障分析模块3,用于分析故障发生位置及故障电阻阻值,判断故障大致区域;
所述仿真模块4,利用Simulink平台通过Artificial Tree智能优化算法获取故障检测数据,显示仿真结果。
所述检测输出模块5,依据仿真结果选用部分数据检测线路神经网络对电网故障检测的准确率分析适应性报告。
本发明利用MATLAB平台用改进后的小波阈值去噪方法去除信号中的噪声源建立电网故障模型,提高信号采集准确度,运行仿真软件进行仿真,并建立基于开关量的电网故障模型得到电网故障信息决策表,训练神经网络,检测故障发生的区间,对复杂电网采取分而治之的策略,使得线路拓扑结构能够灵活多变,这样不仅能提高电网故障诊断的准确性与可靠性,还可以在故障信息存在误动信息的情况下对故障进行检测,从而达到电网故障区域定位的效果,有效的增强了电网故障检测的能力。
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。

Claims (10)

1.一种配电网多重故障诊断方法,其特征在于,包括以下步骤:
步骤S1、利用MATLAB平台对配电网线路进行短路故障分析,得出基于电气量的故障信息决策表;
步骤S2、利用Simulink平台对所述基于电气量的故障信息决策表进行建模仿真,获取配电网线路神经网络的输出训练数据;
步骤S3、采用小波变换法对所述神经网络的输出训练数据进行去噪、采集处理,形成相关故障信息决策表作为神经网络的训练样本;
步骤S4、采用改进的Artificial Tree 智能优化算法优化神经网络的权值和阈值,通过迭代达到相关误差精度的要求,选择部分数据作为故障数据并利用训练好的神经网络检测故障。
2.根据权利要求1所述的配电网多重故障诊断方法,其特征在于:所述步骤S1中,所述MATLAB平台采集配电网线路发生小电流接地短路故障后产生的零序电流、零序功率和零序导纳信号,基于这三种信息形成基于电气量的故障信息决策表。
3.根据权利要求2所述的配电网多重故障诊断方法,其特征在于:依据所述零序电流、零序功率和零序导纳信号用离散采样时间和连续采样时间相混合的采样模式进行故障分析,具体步骤如下:
步骤S201、在每条线路两个终端中的一个终端上采样并存储正负极零序电流;
步骤S202、选择固定长度的时间窗内的采样电流作为线路参考电流,对动态采样电流与线路参考电流进行比较,使用改进Pearson方法实时求解二者的相关系数大小;
步骤S203、设定相关系数阈值,判断相关系数与设定阈值的大小,并输出故障所在线路及故障类型。
4.根据权利要求3所述的配电网多重故障诊断方法,其特征在于:所述步骤S202中,所述改进Pearson方法采用线路相关系数的累加值进行故障检测的相关系数分析,所述改进Pearson方法具体步骤如下:
步骤S2021、定义具有固定长度的时间窗内可移动窗口值
Figure 145548DEST_PATH_IMAGE001
Figure 284406DEST_PATH_IMAGE002
,所述
Figure 470495DEST_PATH_IMAGE003
Figure 720211DEST_PATH_IMAGE004
窗口的末端固定时刻分别为
Figure 705484DEST_PATH_IMAGE005
Figure 15243DEST_PATH_IMAGE006
,将所述
Figure 85967DEST_PATH_IMAGE007
中的采样电流与稳态参考电流进行比较,检测直流微电网的故障,所述
Figure 139374DEST_PATH_IMAGE008
Figure 103787DEST_PATH_IMAGE002
的关系式为:
Figure 584447DEST_PATH_IMAGE009
其中
Figure 876888DEST_PATH_IMAGE010
为线路电流的采样频率,
Figure 733986DEST_PATH_IMAGE011
为采样点的数量,
Figure 428272DEST_PATH_IMAGE012
从故障发生时刻
Figure 705932DEST_PATH_IMAGE013
到故障被检测到时刻
Figure 751249DEST_PATH_IMAGE014
之间的时长;
步骤S2022、采用线路相关系数的累加值进行故障检测的相关系数分析,将稳态参考电流曲线和采样电流曲线的第
Figure 146458DEST_PATH_IMAGE015
个对比值可分别表示为
Figure 695251DEST_PATH_IMAGE016
Figure 783293DEST_PATH_IMAGE017
Figure 50326DEST_PATH_IMAGE018
其中
Figure 108281DEST_PATH_IMAGE019
为第
Figure 777159DEST_PATH_IMAGE015
个暂态线路电流,
Figure 36102DEST_PATH_IMAGE020
为第
Figure 524853DEST_PATH_IMAGE015
个采样线路电流;
步骤S2023、引入调整因子
Figure 153542DEST_PATH_IMAGE021
优化Pearson相关系数调整因子
Figure 676928DEST_PATH_IMAGE022
表达式为:
Figure 841193DEST_PATH_IMAGE023
其中,
Figure 82818DEST_PATH_IMAGE024
为数学常数,表示为自然对数的底数,
Figure 623521DEST_PATH_IMAGE025
为两条曲线所对比点的总数;
步骤S2023、依据所述调整因子
Figure 126046DEST_PATH_IMAGE021
优化相关系数
Figure 726792DEST_PATH_IMAGE026
,所述相关系数
Figure 190134DEST_PATH_IMAGE026
的表达式为:
Figure 534528DEST_PATH_IMAGE027
Figure 32505DEST_PATH_IMAGE028
5.根据权利要求4所述的配电网多重故障诊断方法,其特征在于:所述步骤S2中,依据所述故障分析获取故障位置和故障电阻,通过Simulink平台构建模型输出训练数据,采用小波变换法处理训练数据。
6.根据权利要求5所述的配电网多重故障诊断方法,其特征在于:所述小波变换法处理训练数据的具体步骤为:
步骤S301、对含噪声的训练数据信号
Figure 538573DEST_PATH_IMAGE029
作小波变换,获取一组小波分解系数;
步骤S302、采样软阈值函数对得到的小波分解系数进行阈值处理,得到一组估计小波系数,所述软阈值函数
Figure 380889DEST_PATH_IMAGE030
表达式为:
Figure 263395DEST_PATH_IMAGE031
Figure 615879DEST_PATH_IMAGE032
其中,
Figure 292848DEST_PATH_IMAGE033
为符号函数,
Figure 730782DEST_PATH_IMAGE034
为小波系数,
Figure 541612DEST_PATH_IMAGE035
为阈值,
Figure 748603DEST_PATH_IMAGE036
为噪声强度,
Figure 862052DEST_PATH_IMAGE037
为信号长度,
Figure 787283DEST_PATH_IMAGE038
为幅值因子,
Figure 11591DEST_PATH_IMAGE039
Figure 338667DEST_PATH_IMAGE040
为调节因子;
步骤S303、通过调节
Figure 623018DEST_PATH_IMAGE041
使得在阈值点
Figure 396064DEST_PATH_IMAGE035
处连续且降低对原信号产生的震荡误差,提取去噪后的信号。
7.根据权利要求6所述的配电网多重故障诊断方法,其特征在于:在步骤S4中,将所述震荡误差的均方误差采用Artificial Tree 算法优化获取线路神经网络,所述ArtificialTree 算法优化步骤如下:
步骤S401、以所述震荡误差的均方误差MSE作为目标函数,将收集到的数据进行整理,准备好训练样本和测试样本;
步骤S402、设置Artificial Tree 算法的相关参数,确定分支个数
Figure 424063DEST_PATH_IMAGE042
、空间维数
Figure 730279DEST_PATH_IMAGE043
、迭代次数
Figure 185531DEST_PATH_IMAGE044
,并对分支进行初始化操作,生成初始分支;
步骤S402、计算每个分支所对应的函数值,根据函数值选出最好分支
Figure 85354DEST_PATH_IMAGE045
和最好分支所对应的函数值
Figure 651465DEST_PATH_IMAGE046
步骤S403、判断是否满足误差精度的要求或者是否达到了最大迭代次数,如果都不满足要求则准备进入IAT算法的迭代寻优过程;
步骤S404、如果搜索次数达到了搜索次数的上限,将进行随机操作生成新的分支,将新产生的分支与旧的分支进行比较,如果新的分支优于旧的分支则新的分支将取代旧分支,否则仍以旧的分支为最优分支;
步骤S405、再次判断最优分支是否满足误差精度或者当前迭代次数是否达到了设定的最大迭代次数
Figure 313653DEST_PATH_IMAGE047
8.根据权利要求7所述的配电网多重故障诊断方法,其特征在于:所述ArtificialTree 算法的相关参数依据配电网线路保护装置的动作信息及对应发生故障的区域设定,形成相应的神经网络训练样本。
9.根据权利要求8所述的配电网多重故障诊断方法,其特征在于:所述配电网线路保护装置的动作信息包括基于开关量的电网信号,用0、1对电网开关动作进行编号。
10.一种采用权利要求1-9任一项所述的配电网多重故障诊断方法的系统,其特征在于,包括:电路采样模块(1)、线路保护装置(2)、故障分析模块(3)、仿真模块(4)以及检测输出模块(5);
所述电路采样模块(1),用于采集和预处理导入到所述MATLAB平台和Simulink平台的各项数据;
所述线路保护装置(2),用于配电网发生故障时会进行跳闸或者报警保护动作,根据保护装置不同的保护动作,判断故障发生的区间,并且将电网的区间进行编号;
所述故障分析模块(3),用于分析故障发生位置及故障电阻阻值,判断故障区域;
所述仿真模块(4),利用Simulink平台通过Artificial Tree智能优化算法获取故障检测数据,显示仿真结果;
所述检测输出模块(5),依据仿真结果选用部分数据检测线路神经网络对电网故障检测的准确率分析适应性报告。
CN202210506365.6A 2022-05-11 2022-05-11 一种配电网多重故障诊断方法 Active CN114636900B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210506365.6A CN114636900B (zh) 2022-05-11 2022-05-11 一种配电网多重故障诊断方法
PCT/CN2022/133889 WO2023216553A1 (zh) 2022-05-11 2022-11-24 一种配电网多重故障诊断方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210506365.6A CN114636900B (zh) 2022-05-11 2022-05-11 一种配电网多重故障诊断方法

Publications (2)

Publication Number Publication Date
CN114636900A true CN114636900A (zh) 2022-06-17
CN114636900B CN114636900B (zh) 2022-09-30

Family

ID=81953165

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210506365.6A Active CN114636900B (zh) 2022-05-11 2022-05-11 一种配电网多重故障诊断方法

Country Status (2)

Country Link
CN (1) CN114636900B (zh)
WO (1) WO2023216553A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023216553A1 (zh) * 2022-05-11 2023-11-16 广东电网有限责任公司东莞供电局 一种配电网多重故障诊断方法及系统

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117310353B (zh) * 2023-11-30 2024-02-09 淮安苏达电气有限公司 一种变电站一、二次回路通流加压故障测试方法及系统
CN117630579A (zh) * 2023-12-06 2024-03-01 海南电力产业发展有限责任公司 一种基于分布式行波检测的配电网故障精确定位方法
CN117388641B (zh) * 2023-12-11 2024-03-08 武汉邢仪新未来电力科技股份有限公司 一种基于分布式物联网终端的线路故障远程监测系统
CN117390407B (zh) * 2023-12-13 2024-04-05 国网山东省电力公司济南供电公司 变电站设备的故障识别方法、系统、介质和设备
CN117538689A (zh) * 2024-01-10 2024-02-09 高唐县恒诚建筑工程有限公司 一种基于物联网的电路故障检测方法、系统、装置及介质
CN117668622B (zh) * 2024-02-01 2024-05-10 山东能源数智云科技有限公司 设备故障诊断模型的训练方法、故障诊断方法及装置
CN117909909B (zh) * 2024-03-19 2024-06-14 青岛鼎信通讯股份有限公司 一种配电网弧光接地故障识别方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508076A (zh) * 2011-11-09 2012-06-20 东北大学 一种基于多智能体系统和小波分析的故障诊断装置及方法
CN103713237A (zh) * 2013-12-30 2014-04-09 华北电力大学 一种电力系统输电线路短路故障诊断方法
CN104502795A (zh) * 2014-11-26 2015-04-08 国家电网公司 一种适用于微电网的智能故障诊断方法
US20160041216A1 (en) * 2013-03-29 2016-02-11 Beijing Inhand Networks Technology Co., Ltd. Method and system for detecting and locating single-phase ground fault on low current grounded power-distribution network
CN111239549A (zh) * 2020-02-18 2020-06-05 国网信通亿力科技有限责任公司 一种基于离散小波变换的配电故障快速定位方法
CN111537830A (zh) * 2020-03-26 2020-08-14 南方电网科学研究院有限责任公司 一种基于云边架构和小波神经网络的配电网故障诊断方法
CN111537841A (zh) * 2020-06-30 2020-08-14 上海交通大学 一种适用于接地故障类型识别的优化方法及系统
US20200300907A1 (en) * 2016-12-29 2020-09-24 Hefei University Of Technology Analog-circuit fault diagnosis method based on continuous wavelet analysis and elm network
CN111896890A (zh) * 2020-08-06 2020-11-06 山东大学 基于希尔伯特-黄变换的微电网线路故障诊断方法及系统
CN113902250A (zh) * 2021-09-03 2022-01-07 广东电网有限责任公司 一种基于人工神经网络的风险评估分析方法
CN113902050A (zh) * 2021-08-30 2022-01-07 广东电网有限责任公司 一种适用于电网故障预案信息的检测提取方法
WO2022064350A1 (en) * 2020-09-22 2022-03-31 Victoria Link Limited Electricity distribution network fault detection system and method

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05298105A (ja) * 1992-04-22 1993-11-12 Mitsubishi Electric Corp 故障原因推論装置
DE19905814A1 (de) * 1999-02-12 2000-08-17 Deutsche Telekom Ag Verfahren zur Überwachung der Übertragungsqualität eines optischen Übertragungssystems, insbesondere eines optischen Wellenlängenmultiplexnetzes
CN103778479A (zh) * 2014-01-10 2014-05-07 国网上海市电力公司 自适应信息容错保护方法
CN104360227B (zh) * 2014-10-31 2017-05-24 国家电网公司 基于行波法和暂态主频法的变电站电缆出线故障监测方法
CN110457750B (zh) * 2019-07-09 2023-10-24 国家电网有限公司 基于神经网络响应面的有载分接开关弹簧储能不足故障识别方法
CN110533549B (zh) * 2019-07-12 2022-04-29 杭州电子科技大学 一种基于ga优化rs-bp神经网络电网故障定位方法
CN111339872A (zh) * 2020-02-18 2020-06-26 国网信通亿力科技有限责任公司 一种基于分类模型的电网故障分类方法
CN112749509B (zh) * 2020-12-30 2022-06-10 西华大学 一种基于lstm神经网络的智能变电站故障诊断方法
CN114636900B (zh) * 2022-05-11 2022-09-30 广东电网有限责任公司东莞供电局 一种配电网多重故障诊断方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508076A (zh) * 2011-11-09 2012-06-20 东北大学 一种基于多智能体系统和小波分析的故障诊断装置及方法
US20160041216A1 (en) * 2013-03-29 2016-02-11 Beijing Inhand Networks Technology Co., Ltd. Method and system for detecting and locating single-phase ground fault on low current grounded power-distribution network
CN103713237A (zh) * 2013-12-30 2014-04-09 华北电力大学 一种电力系统输电线路短路故障诊断方法
CN104502795A (zh) * 2014-11-26 2015-04-08 国家电网公司 一种适用于微电网的智能故障诊断方法
US20200300907A1 (en) * 2016-12-29 2020-09-24 Hefei University Of Technology Analog-circuit fault diagnosis method based on continuous wavelet analysis and elm network
CN111239549A (zh) * 2020-02-18 2020-06-05 国网信通亿力科技有限责任公司 一种基于离散小波变换的配电故障快速定位方法
CN111537830A (zh) * 2020-03-26 2020-08-14 南方电网科学研究院有限责任公司 一种基于云边架构和小波神经网络的配电网故障诊断方法
CN111537841A (zh) * 2020-06-30 2020-08-14 上海交通大学 一种适用于接地故障类型识别的优化方法及系统
CN111896890A (zh) * 2020-08-06 2020-11-06 山东大学 基于希尔伯特-黄变换的微电网线路故障诊断方法及系统
WO2022064350A1 (en) * 2020-09-22 2022-03-31 Victoria Link Limited Electricity distribution network fault detection system and method
CN113902050A (zh) * 2021-08-30 2022-01-07 广东电网有限责任公司 一种适用于电网故障预案信息的检测提取方法
CN113902250A (zh) * 2021-09-03 2022-01-07 广东电网有限责任公司 一种基于人工神经网络的风险评估分析方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
QIN HAO 等: "High impedance fault line selection method for resonant grounding system based on wavelet packet analysis", 《 2018 CHINA INTERNATIONAL CONFERENCE ON ELECTRICITY DISTRIBUTION (CICED)》 *
ZHILING YIN 等: "Research on Network Topology Analysis Method of Distribution Management Based on GIS", 《2021 INTERNATIONAL CONFERENCE ON INTELLIGENT TRANSPORTATION, BIG DATA & SMART CITY (ICITBS)》 *
孙其东: "小波包能量和LM-BP神经网络的发电机机端故障类型识别", 《电子质量》 *
罗金满 等: "考虑源网荷储协调的主动配电网优化调度方法研究", 《电力系统保护与控制》 *
訾小娜: "配电网单相接地故障定位方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023216553A1 (zh) * 2022-05-11 2023-11-16 广东电网有限责任公司东莞供电局 一种配电网多重故障诊断方法及系统

Also Published As

Publication number Publication date
WO2023216553A1 (zh) 2023-11-16
CN114636900B (zh) 2022-09-30

Similar Documents

Publication Publication Date Title
CN114636900B (zh) 一种配电网多重故障诊断方法
CN109543317B (zh) 一种pemfc剩余使用寿命预测的方法及装置
CN106932686B (zh) 基于最大最小特征值法的电网异常快速检测方法
Asada et al. Identifying multiple interacting bad data in power system state estimation
CN112415414A (zh) 一种锂离子电池剩余使用寿命预测方法
CN110907754B (zh) 一种基于psd-bpa的故障线路严重程度评估方法
Purushothama et al. ANN applications in fault locators
Mirshekali et al. Linear and nonlinear fault location in smart distribution network under line parameter uncertainty
CN113363980B (zh) 一种适用于低压配电网拓扑自动识别方法及设备
CN111426905B (zh) 一种配电网同母线变关系异常诊断方法、装置及系统
CN112180204A (zh) 一种基于电气量信息的电网线路故障诊断方法
CN112098889B (zh) 一种基于神经网络和特征矩阵的单相接地故障定位方法
Wei et al. Fault detection based on the generalized S-transform with a variable factor for resonant grounding distribution networks
Kezunovic Translational knowledge: From collecting data to making decisions in a smart grid
CN110726898A (zh) 一种配电网故障类型识别方法
CN113937764A (zh) 一种低压配电网高频量测数据处理与拓扑识别方法
CN110108971B (zh) 一种架空裸导线经树木接地故障的电弧检测方法
Coban et al. Artificial neural network based fault location on 230 kv transmission line using voltage and current signals
CN106019043B (zh) 一种基于模糊积分多源信息融合理论的电网故障诊断方法
CN111935762B (zh) 一种5g承载网下基于ewt和cnn的配网故障诊断方法和系统
CN106646106B (zh) 基于变点探测技术的电网故障检测方法
KR20150037415A (ko) 배전계통의 상태추정 장치
Nagam et al. Gaussian process regression based fault location in DC microgrid
CN115795360A (zh) 基于人工神经网络的电缆故障检测方法
Bohórquez et al. One-ended fault location method based on machine learning models

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