CN111337768A - 变压器油中溶解气体的深度并行故障诊断方法及系统 - Google Patents

变压器油中溶解气体的深度并行故障诊断方法及系统 Download PDF

Info

Publication number
CN111337768A
CN111337768A CN202010134616.3A CN202010134616A CN111337768A CN 111337768 A CN111337768 A CN 111337768A CN 202010134616 A CN202010134616 A CN 202010134616A CN 111337768 A CN111337768 A CN 111337768A
Authority
CN
China
Prior art keywords
diagnosis
data
lstm
diagnostic model
dissolved gas
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
CN202010134616.3A
Other languages
English (en)
Other versions
CN111337768B (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202010134616.3A priority Critical patent/CN111337768B/zh
Publication of CN111337768A publication Critical patent/CN111337768A/zh
Application granted granted Critical
Publication of CN111337768B publication Critical patent/CN111337768B/zh
Priority to US17/160,415 priority patent/US11656298B2/en
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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004Gaseous mixtures, e.g. polluted air
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • 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/62Testing of transformers
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Computational Linguistics (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Combustion & Propulsion (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Power Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

本发明公开了一种变压器油中溶解气体的深度并行故障诊断方法及系统,属于变压器故障诊断领域。采集各变电站油中溶解气体监测信息,并对数据进行归一化处理;用油中溶解气体构建特征参量作为LSTM诊断模型的输入,并进一步将数据进行图像化处理作为CNN诊断模型的输入;分别构建LSTM诊断模型和CNN诊断模型,并按比例用数据集分别对诊断模型进行训练及验证;利用DS证据理论算法对两种深度学习模型的softmax层输出进行深度并行融合。本发明综合考虑不同深度学习方法对复杂非线性关系的提取能力,构建了基于CNN‑LSTM深度并行融合的变压器故障诊断模型。该模型不需要经过复杂的特征设计,便能得到较高的准确率,因而模型具有较高的可移植性。

Description

变压器油中溶解气体的深度并行故障诊断方法及系统
技术领域
本发明属于电力变压器故障诊断领域,更具体地,涉及一种变压器油中溶解气体的深度并行故障诊断方法及系统。
背景技术
随着计算机存储及传感器技术的不断发展,电力变压器的在线监测数据将呈现出爆炸式增长态势。作为电力系统的枢纽性装备,电力变压器的健康状态直接影响着整个电网的安全稳定。因此,综合运用人工智能、大数据等先进技术对大容量电力变压器的在线监测数据进行故障诊断与状态预测,是指导变压器运维检修相关工作的一大研究热点。油中溶解气体分析(Dissolved Gas Analysis,DGA)是一项对变压器绝缘油中溶解气体成分及含量进行分析的在线监测技术。目前基于DGA数据的传统诊断方法已进行了大量的研究。传统的比值方法(例如IEC比率法,Rogers比率法和Doernenburg比率法)简单易行,并已在工程实践中广泛使用。但是,仍然存在一些问题,例如诊断边界过于绝对和缺码漏码等。传统智能方法因学习能力、处理效率、特征提取等方面的不足,在推广应用时仍存在一定的局限性,如神经网络容易陷入局部最优解,K近邻法在高维空间中计算效率低下,支持向量机本质是一个二分类器,处理多分类问题较麻烦等等。传统比值方法与传统智能方法的组合诊断也有相关的研究,但效果仍然不够理想。近年来,由于深度学习对复杂非线性特征具有强大的提取能力,因此将其引入变压器故障诊断领域,其故障诊断准确率相较传统机器学习的诊断模型有了显著提升。如目前的以DGA数据的无编码比值作为输入矩阵的基于ReLU-DBN的新型故障诊断方法,相较传统方法,准确率大大提高。也有基于CNN的变压器内部故障诊断,将CNN用于增强差动保护性能,并通过不同的可靠性指标来验证该方法的有效性。
卷积神经网络(Convolutional Neural Networks,CNN)和循环神经网络(Recurrent Neural Network,RNN)是深度学习中两种应用最为广泛的框架。CNN在视觉辨识、图像处理、故障诊断等领域已经应用得较为成熟。而LSTM作为RNN的改进模型,用以弥补RNN模型在反馈信息过程中出现的长期记忆缺失、梯度消散或爆炸的缺陷,如今,在语音辨识、视频分类和股票预测等各个领域得到广泛应用,也取得了不错的效果。
因此,如何将CNN及LSTM引进变压器DGA故障诊断建模,进而提高故障诊断准确率是目前亟需解决的技术难题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提出了一种变压器油中溶解气体的深度并行故障诊断方法及系统,由此解决如何将CNN及LSTM引进变压器DGA故障诊断建模,进而提高故障诊断准确率的技术问题。
为实现上述目的,按照本发明的一个方面,提供了一种变压器油中溶解气体的深度并行故障诊断方法,包括:
(1)获取各变压器油中溶解气体的若干组监测数据,对各组所述监测数据中的溶解气体含量进行分析,获取对应的故障类型标签,对各组所述监测数据进行归一化处理后与对应的故障类型标签组成目标数据集;
(2)将所述目标数据集划分为第一训练集和第一验证集,由所述第一训练集对LSTM诊断模型进行训练,由所述第一验证集对训练后的LSTM诊断模型进行验证;
(3)对所述目标数据集中的各组数据进行图像化处理,得到图像数据集,将所述图像数据集划分为第二训练集和第二验证集,由所述第二训练集对CNN诊断模型进行训练,由所述第二验证集对训练后的CNN诊断模型进行验证;
(4)分别对训练后的LSTM诊断模型及训练后的CNN诊断模型的softmax层的输出进行深度并行融合,根据最大置信原则,输出最终的诊断结果。
优选地,所述各变压器油中溶解气体的若干组监测数据为:datai={ai,1,ai,2,...,ai,j,...,ai,N,si}i∈[1,K],其中,K为K组溶解气体监测数据,ai,j为第i组溶解气体监测数据中第j(j∈[1,N])种气体参量的含量,si为第i组溶解气体监测数据对应的变压器状态,N为气体种类数。
优选地,归一化后得到的数据为:datai'={bi,1,bi,2,...,bi,j,...,bi,N,si}i∈[1,K],其中,K为K组溶解气体监测数据,bi,j为第i组溶解气体监测数据中第j(j∈[1,N])种气体参量的含量的归一化值,si为第i组溶解气体监测数据对应的变压器状态,N为气体种类数。
优选地,所述LSTM诊断模型包括:输入层、LSTM层、全连接层、softmax层及分类输出层,其中,LSTM层拥有若干个隐藏单元,各隐藏单元的状态激活函数为tanh,门激活函数为sigmoid,softmax层将所述LSTM诊断模型中全连接层的输出作为输入向量,输出所述LSTM诊断模型对故障标签的诊断支持度。
优选地,由
Figure BDA0002396896810000031
确定所述LSTM诊断模型中softmax层输出的故障标签的诊断支持度,其中,x1表示所述LSTM诊断模型中全连接层的输出,θi,i=1,2,...,N为LSTM诊断模型中softmax层的权值矩阵,Softmax为激活函数。
优选地,在步骤(3)中,对所述目标数据集中的各组数据进行图像化处理,得到图像数据集,包括:
对所述目标数据集中的各组数据进行图像化处理,将数据的差异分别用颜色进行图像化呈现得到图像数据集A,用高度进行图像化呈现得到图像数据集B,其中,所述图像数据集A表示为:datai”={ci,si}i∈[1,K],所述图像数据集B表示为:datai”'={di,si}i∈[1,K],其中,K为K组溶解气体监测数据,ci为所述图像数据集A中第i组溶解气体监测数据参量转化的图像,di为所述图像数据集B中第i组溶解气体监测数据参量转化的图像,si为第i组溶解气体监测数据对应的变压器状态。
优选地,所述CNN诊断模型中的softmax层将CNN诊断模型中全连接层的输出作为输入向量,由
Figure BDA0002396896810000041
得到故障标签的诊断支持度,其中,x2表示CNN诊断模型中全连接层的输出,θ'i,i=1,2,...,N为CNN诊断模型中softmax层的权值矩阵,Softmax为激活函数。
优选地,步骤(4)包括:
(4.1)分别获取LSTM诊断模型及CNN诊断模型的softmax层对同一组待诊断监测数据对应的故障标签的诊断支持度:ξk,γ=ξk,1,ξk,2,…,ξk,l,其中,k的数值代表不同的诊断方法,l为故障标签总数,γ=1,2,…,l;
(4.2)以不同的诊断方法k为行,诊断支持度ξk,γ为列,组成支持度矩阵{(Hγk,γ),k=1,2,γ=1,...,l},其中,所述支持度矩阵的每个元素表示:第k个诊断方法对故障标签Hγ的支持度为ξk,γ
(4.3)将所述支持度矩阵中各列看作DS证据理论的识别框架Θ,则Θ={Hγ|γ=1,2,…,l}={H1,H2,…,Hl};
(4.4)将不同诊断方法的诊断支持度信息归入到同一识别框架Θ下,得到基本概率分配;
(4.5)由所述基本概率分配得到复合概率分配,进而由所述复合概率分配得到融合后的不同故障标签的置信度,将最大置信度对应的故障标签作为所述待诊断监测数据的诊断结果。
优选地,由
Figure BDA0002396896810000051
得到基本概率分配,其中,ωk为第k个诊断模型的权重值,mk,γ表示第k个诊断模型对于变压器的故障标签Hγ的基本概率分配,mk表示计算基本概率分配的函数,mk,H表示分配给整个故障标签集H的剩余概率。
优选地,由
Figure BDA0002396896810000052
得到复合概率分配,其中,mγ表示故障标签Hγ的复合概率分配,K为中间变量;由
Figure BDA0002396896810000053
计算出故障标签Hγ的归一化概率分配,最终根据最大置信原则,输出最大概率对应的标签作为最终诊断结果。
按照本发明的另一方面,提供了一种变压器油中溶解气体的深度并行故障诊断系统,包括:
数据处理模块,用于获取各变压器油中溶解气体的若干组监测数据,对各组所述监测数据中的溶解气体含量进行分析,获取对应的故障类型标签,对各组所述监测数据进行归一化处理后与对应的故障类型标签组成目标数据集;
第一训练模块,用于将所述目标数据集划分为第一训练集和第一验证集,由所述第一训练集对LSTM诊断模型进行训练,由所述第一验证集对训练后的LSTM诊断模型进行验证;
第二训练模块,用于对所述目标数据集中的各组数据进行图像化处理,得到图像数据集,将所述图像数据集划分为第二训练集和第二验证集,由所述第二训练集对CNN诊断模型进行训练,由所述第二验证集对训练后的CNN诊断模型进行验证;
诊断模块,用于分别对训练后的LSTM诊断模型及训练后的CNN诊断模型的softmax层的输出进行深度并行融合,根据最大置信原则,输出最终的诊断结果。
按照本发明的另一方面,提供了一种计算机可读存储介质,其上存储有程序指令,所述程序指令被处理器执行时实现如上述任一所述的变压器油中溶解气体的深度并行故障诊断方法。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:本发明充分利用了不同深度学习方法对复杂非线性特征的提取能力,提高了变压器故障诊断的准确率。所采用的特征参量无需经过复杂的特征设计,便能得到较高的诊断准确率,具有较高的可推广性。
附图说明
图1是本发明实施例提供的一种变压器油中溶解气体的深度并行故障诊断方法的流程示意图;
图2是本发明实施例提供的一种LSTM网络结构示意图;
图3是本发明实施例提供的一种CNN单元结构示意图;
图4是本发明实施例提供的一种数据图像化及CNN迁移训练示意图;
图5是本发明实施例提供的一种深度并行融合示意图;
图6是本发明实施例提供的一种系统结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
在本发明实例中,“第一”、“第二”等是用于区别不同的对象,而不必用于描述特定的顺序或先后次序。
本发明提出一种变压器油中溶解气体的深度并行故障诊断方法。首先采集了DGA故障数据集,并对数据进行归一化、图像化处理;接着,利用新数据集分别对LSTM、CNN诊断网络进行训练及验证,得到所需的针对DGA的LSTM、CNN故障诊断网络;最终将两种深度学习网络的输出层去掉,利用DS证据理论多源信息融合对softmax层的输出进行深度并行融合,诊断结果取融合结果最大置信度对应的故障标签。这一深度并行融合的方法,充分利用了不同深度学习方法对复杂非线性特征的提取能力,有效提高了变压器故障诊断的准确率。所采用的特征参量无需经过复杂的特征设计,便能得到较高的诊断准确率,具有较高的可推广性。
如图1所示是本发明实施例提供的一种变压器油中溶解气体的深度并行故障诊断方法的流程示意图,包括以下步骤:
S1:采集各变压器油中溶解气体的监测数据,对变压器油中溶解气体含量进行分析,获取对应故障标签,并对监测数据进行归一化处理;
在本发明实施例中,数据可以取自历年相关文献及电力公司实际测试数据,每一组数据包括氢气(H2)、甲烷(CH4)、乙烷(C2H6)、乙烯(C2H4)、乙炔(C2H2)五种关键状态参量的含量及其对应的变压器状态;
在本发明实施例中,步骤S1中溶解气体的监测数据可以表示为:
datai={ai,1,ai,2,...,ai,j,,...,ai,N,si}i∈[1,K]
其中,K为K组溶解气体监测数据,ai,j为第i组溶解气体监测数据中第j(j∈[1,N])种气体参量的含量,si为第i组溶解气体监测数据对应的变压器状态,N为气体种类数;
作为一种可选的实施方式,步骤S1中归一化的方法可以采取公式(1)进行计算:
Figure BDA0002396896810000081
其中,Xi表示数据样本,E(X)表示均值,D(X)表示方差,zi表示正则化值。
归一化后得到的数据表示为:
datai'={bi,1,bi,2,...,bi,j,...bi,N,si}i∈[1,K]
其中,K为K组溶解气体监测数据,bi,j为归一化后第i组溶解气体监测数据中第j(j∈[1,N])种气体参量的含量的归一化值,si为第i组溶解气体监测数据对应的变压器状态。
在本发明实施例中,依据常用的IEC60599标准,将变压器故障模式分为6种,分别为高温过热、中温过热、低温过热、局部放电、低能放电、高能放电,如表1所示:
表1:故障模式编码
Figure BDA0002396896810000082
S2:用油中溶解气体构建特征参量,按比例将归一化后的新数据集划分为训练集和验证集,作为LSTM诊断模型的输入;
将步骤S2中获取的输入数据集划分为两部分,其中,取η%作为训练集,对LSTM网络进行训练,(1-η%)作为测试集验证网络对变压器DGA故障诊断的准确率,例如,80%作为训练集对Bi-LSTM网络进行训练,20%作为验证集验证网络对LSTM诊断模型诊断效果;
S3:构建LSTM诊断模型,输入训练集和验证集进行网络训练;
在步骤S3中,构建基于LSTM的深度学习框架的具体方法为:
构建五层网络,如图2所示,包括:输入层、LSTM层、全连接层、softmax层、分类输出层,其中LSTM层拥有若干个隐藏单元,其状态激活函数为‘tanh’,门激活函数为‘sigmoid’。网络的超参数设置如表2所示;
表2:LSTM网络超参数设置
名称 取值
优化器 Adam
梯度阈值 0.01
隐藏单元数 150
Minibatchsize 100
Maxepoch 500
softmax层将LSTM诊断模型中全连接层的输出作为输入向量,输出LSTM诊断模型对故障标签的诊断支持度,具体如下:
Figure BDA0002396896810000091
其中,x1表示LSTM诊断模型中全连接层的输出,即softmax层的输入,θi,i=1,2,...,N为LSTM诊断模型中softmax层的权值矩阵,Softmax为激活函数。
S4:用油中溶解气体构建特征参量,对数据进行图像化处理,得到图像数据集,按比例将图像数据集划分为训练集和验证集,作为CNN诊断模型的输入;
作为一种可选的实施方式,步骤S4中对数据进行图像化处理,将数据的差异分别用颜色和高度进行图像化呈现,分别得到图像数据集A和图像数据集B,分别取图像数据集A和图像数据集B中的一个图像,分别表示为:
datai”={ci,si}i∈[1,K]
其中,K为K组溶解气体监测数据,ci为图像数据集A中第i组溶解气体监测数据参量转化的图像,si为第i组溶解气体监测数据对应的变压器状态。
datai”'={di,si}i∈[1,K]
其中,K为K组溶解气体监测数据,di为图像数据集B中第i组溶解气体监测数据参量转化的图像,si为第i组溶解气体监测数据对应的变压器状态。
S5:构建CNN诊断模型,输入训练集和验证集,对CNN诊断模型进行训练;倘若数据量足够大,可构建最基本的CNN诊断模型,并利用大数据对CNN诊断模型进行训练及验证,倘若数据量较小,可利用迁移学习对经典模型如MobileNet-V2进行迁移训练及验证;
在本发明实施例中,CNN诊断模型还可以选用其它模型,本发明实施例不做唯一性限定。
如图3所示是本发明实施例提供的一种CNN单元结构示意图,具体为:
MobileNet-V2是基于MobileNet-V1的基本概念构建,MobileNet-V2提出两大富有创新性的设计思想:Inverted residuals和Linear bottlenecks。
Inverted residuals:MobileNet-V2的Depth-wise convolution之前有一个1x1的“扩张”层,其目的是为了在数据进入深度卷积之前扩展数据中的通道数量,丰富特征数量,从而提高精度。
Linear bottlenecks:MobileNet-V2提出在通道数较少的层后,用线性激活函数代替ReLU激活函数。因为“扩张层”的引入,卷积层输出的大量特征需要进行“压缩”操作以减少计算量,随着通道数目减少,倘若激活函数还选取ReLu,特征将会受到破坏。这是由于ReLu对于负的输入,输出全为零;而本来特征就已经被“压缩”,再经过ReLu的话,又要“损失”一部分特征。
具体地,softmax层将CNN诊断模型中全连接层的输出作为输入向量,输出CNN诊断模型对故障标签的诊断支持度,具体如下:
Figure BDA0002396896810000111
其中,x2表示CNN诊断模型中全连接层的输出,即softmax层的输入,θ'i,i=1,2,...,N为CNN诊断模型中softmax层的权值矩阵,Softmax为激活函数。
作为一种可选的实施方式,如图4所示,步骤S5利用迁移学习对经典模型MobileNet-V2进行迁移训练及验证,首先对MobileNet-V2的最后一层全连接层进行替换,令其输出数量等于诊断标签类别数量,在本实施例为6。接着冻结预训练网络前50%层的参数,最后对基于MobileNet-V2模型的迁移训练超参数进行设置如表3所示,输入训练集对网络未冻结部分进行训练。
表3:CNN诊断模型超参数设置
Figure BDA0002396896810000112
Figure BDA0002396896810000121
S6:利用DS证据理论对两种深度学习模型的softmax层输出进行深度并行融合,经输出层输出最终的诊断结果;
作为一种可选的实施方式,步骤S6中的深度并行融合借助DS证据理论对CNN及LSTM诊断网络的softmax层输出的对故障标签的诊断支持度进行融合。如图5所示,步骤S6具体可以采用以下方式实现:
(1)获取各深度学习诊断网络的softmax层输出:即CNN和LSTM对同一组DGA监测数据对应的故障标签的诊断支持度:ξk,γ=ξk,1,ξk,2,…,ξk,l,(k的数值代表不同方法,l为故障标签总数,即γ=1,2,…,l)
(2)以不同方法k为行,诊断支持度ξk,γ为列,组成支持度矩阵,即:{(Hγk,γ),k=1,2,γ=1,...,l},支持度矩阵的每个元素表示:第k个诊断方法对故障标签Hγ的支持度为ξk,γ
(3)将支持度矩阵中各列看作DS证据理论的识别框架Θ,从而有Θ={Hγ|γ=1,2,…,l}={H1,H2,…,Hl};
(4)将不同诊断方法的诊断支持度信息归入到同一识别框架Θ下,计算基本概率分配m(Basic Probability Assignment,BPA):
mk,γ=mk(Hγ)=ωkξk,γ (2)
Figure BDA0002396896810000122
Figure BDA0002396896810000123
Figure BDA0002396896810000124
其中,ωk为各个诊断方法k(k=1、2)的权重值,一般情况下可以取平均值,即:ωk=1/2;mk,γ表示第k种方法对于评估目标(即变压器的故障状态)Hγ的基本概率分配;mk表示计算基本概率分配的函数;mk,H表示分配给整个故障标签集H,而不是特定变压器故障状态的剩余概率;
Figure BDA0002396896810000131
是用于计算基本概率分配的中间变量。
(5)复合概率分配:
Figure BDA0002396896810000132
Figure BDA0002396896810000133
Figure BDA0002396896810000134
Figure BDA0002396896810000135
其中,mγ表示故障标签的复合概率分配,K、
Figure BDA0002396896810000136
是用于计算复合概率分配的中间变量。
最后归一化处理,得到综合诊断结果:
Figure BDA0002396896810000137
Figure BDA0002396896810000138
即ξγ,γ=1,…,l为融合得到的不同故障标签置信度的归一化结果,
Figure BDA0002396896810000139
表示不确定性分布的归一化值。
最后将最大置信度对应的故障标签输出作为最终的诊断结果。
最终得本例模型LSTM故障诊断准确率为93.6%,CNN故障诊断准确率为92.1%,深度并行诊断准确率为96.9%。
如图6所示是本发明实施例提供的一种变压器油中溶解气体的深度并行故障诊断系统结构示意图,包括:
数据处理模块201,用于获取各变压器油中溶解气体的若干组监测数据,对各组监测数据中的溶解气体含量进行分析,获取对应的故障类型标签,对各组监测数据进行归一化处理后与对应的故障类型标签组成目标数据集;
第一训练模块202,用于将目标数据集划分为第一训练集和第一验证集,由第一训练集对LSTM诊断模型进行训练,由第一验证集对训练后的LSTM诊断模型进行验证;
第二训练模块203,用于对目标数据集中的各组数据进行图像化处理,得到图像数据集,将图像数据集划分为第二训练集和第二验证集,由第二训练集对CNN诊断模型进行训练,由第二验证集对训练后的CNN诊断模型进行验证;
诊断模块204,用于分别对训练后的LSTM诊断模型及训练后的CNN诊断模型的softmax层的输出进行深度并行融合,经输出层输出最终的诊断结果。
其中,各模块的具体实施方式可以参考方法实施例的描述,本发明实施例将不再复述。
在本发明的另一实施例中,还提供了一种计算机可读存储介质,其上存储有程序指令,该程序指令被处理器执行时实现上述方法实施例中的变压器油中溶解气体的深度并行故障诊断方法。
需要指出,根据实施的需要,可将本申请中描述的各个步骤/部件拆分为更多步骤/部件,也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件,以实现本发明的目的。
上述根据本发明的方法可在硬件、固件中实现,或者被实现为可存储在记录介质(诸如CD ROM、RAM、软盘、硬盘或磁光盘)中的软件或计算机代码,或者被实现通过网络下载的原始存储在远程记录介质或非暂时机器可读介质中并将被存储在本地记录介质中的计算机代码,从而在此描述的方法可被存储在使用通用计算机、专用处理器或者可编程或专用硬件(诸如ASIC或FPGA)的记录介质上的这样的软件处理。可以理解,计算机、处理器、微处理器控制器或可编程硬件包括可存储或接收软件或计算机代码的存储组件(例如,RAM、ROM、闪存等),当所述软件或计算机代码被计算机、处理器或硬件访问且执行时,实现在此描述的处理方法。此外,当通用计算机访问用于实现在此示出的处理的代码时,代码的执行将通用计算机转换为用于执行在此示出的处理的专用计算机。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种变压器油中溶解气体的深度并行故障诊断方法,其特征在于,包括:
(1)获取各变压器油中溶解气体的若干组监测数据,对各组所述监测数据中的溶解气体含量进行分析,获取对应的故障类型标签,对各组所述监测数据进行归一化处理后与对应的故障类型标签组成目标数据集;
(2)将所述目标数据集划分为第一训练集和第一验证集,由所述第一训练集对LSTM诊断模型进行训练,由所述第一验证集对训练后的LSTM诊断模型进行验证;
(3)对所述目标数据集中的各组数据进行图像化处理,得到图像数据集,将所述图像数据集划分为第二训练集和第二验证集,由所述第二训练集对CNN诊断模型进行训练,由所述第二验证集对训练后的CNN诊断模型进行验证;
(4)分别对训练后的LSTM诊断模型及训练后的CNN诊断模型的softmax层的输出进行深度并行融合,根据最大置信原则,输出最终的诊断结果。
2.根据权利要求1所述的方法,其特征在于,所述各变压器油中溶解气体的若干组监测数据为:datai={ai,1,ai,2,...,ai,j,...,ai,N,si}i∈[1,K],其中,K为K组溶解气体监测数据,ai,j为第i组溶解气体监测数据中第j(j∈[1,N])种气体参量的含量,si为第i组溶解气体监测数据对应的变压器状态,N为气体种类数;归一化后得到的数据为:datai'={bi,1,bi,2,...,bi,j,...,bi,N,si}i∈[1,K],其中,K为K组溶解气体监测数据,bi,j为第i组溶解气体监测数据中第j(j∈[1,N])种气体参量的含量的归一化值,si为第i组溶解气体监测数据对应的变压器状态,N为气体种类数。
3.根据权利要求2所述的方法,其特征在于,所述LSTM诊断模型包括:输入层、LSTM层、全连接层、softmax层及分类输出层,其中,LSTM层由拥有若干个隐藏单元,各隐藏单元的状态激活函数为tanh,门激活函数为sigmoid,softmax层将所述LSTM诊断模型中全连接层的输出作为输入向量,由
Figure FDA0002396896800000021
得到所述LSTM诊断模型对故障标签的诊断支持度,其中,x1表示所述LSTM诊断模型中全连接层的输出,θi,i=1,2,...,N为LSTM诊断模型中softmax层的权值矩阵,Softmax为激活函数。
4.根据权利要求3所述的方法,其特征在于,在步骤(3)中,对所述目标数据集中的各组数据进行图像化处理,得到图像数据集,包括:
对所述目标数据集中的各组数据进行图像化处理,将数据的差异分别用颜色进行图像化呈现得到图像数据集A,用高度进行图像化呈现得到图图像数据集B,其中,所述图像数据集A表示为:datai”={ci,si}i∈[1,K],所述图像数据集B表示为:datai”'={di,si}i∈[1,K],其中,K为K组溶解气体监测数据,ci为所述图像数据集A中第i组溶解气体监测数据参量转化的图像,di为所述图像数据集B中第i组溶解气体监测数据参量转化的图像,si为第i组溶解气体监测数据对应的变压器状态。
5.根据权利要求4所述的方法,其特征在于,所述CNN诊断模型中的softmax层将CNN诊断模型中全连接层的输出作为输入向量,由
Figure FDA0002396896800000022
得到故障标签的诊断支持度,其中,x2表示CNN诊断模型中全连接层的输出,θ'i,i=1,2,...,N为CNN诊断模型中softmax层的权值矩阵,Softmax为激活函数。
6.根据权利要求5所述的方法,其特征在于,步骤(4)包括:
(4.1)分别获取LSTM诊断模型及CNN诊断模型的softmax层对同一组待诊断监测数据对应的故障标签的诊断支持度:ξk,γ=ξk,1,ξk,2,…,ξk,l,其中,k的数值代表不同的诊断模型,l为故障标签总数,γ=1,2,…,l;
(4.2)以不同的诊断方法k为行,诊断支持度ξk,γ为列,组成支持度矩阵{(Hγk,γ),k=1,2,γ=1,...,l},其中,所述支持度矩阵的每个元素表示:第k个诊断模型对故障标签Hγ的支持度为ξk,γ
(4.3)将所述支持度矩阵中各列看作DS证据理论的识别框架Θ,则Θ={Hγ|γ=1,2,…,l}={H1,H2,…,Hl};
(4.4)将不同诊断模型的诊断支持度信息归入到同一识别框架Θ下,得到基本概率分配;
(4.5)由所述基本概率分配得到复合概率分配,进而由所述复合概率分配得到融合后的不同故障标签的置信度,将最大置信度对应的故障标签作为所述待诊断监测数据的诊断结果。
7.根据权利要求6所述的方法,其特征在于,由
Figure FDA0002396896800000031
得到基本概率分配,其中,ωk为第k个诊断模型的权重值,mk,γ表示第k个诊断模型对于变压器的故障标签Hγ的基本概率分配,mk表示计算基本概率分配的函数,mk,H表示分配给整个故障标签集H的剩余概率。
8.根据权利要求7所述的方法,其特征在于,由{Hγ}:
Figure FDA0002396896800000041
得到复合概率分配,其中,mγ表示故障标签Hγ的复合概率分配,K为中间变量;由
Figure FDA0002396896800000042
计算出故障标签Hγ的归一化概率分配ξγ,最终根据最大置信原则,输出最大概率对应的标签作为最终诊断结果。
9.一种变压器油中溶解气体的深度并行故障诊断系统,其特征在于,包括:
数据处理模块,用于获取各变压器油中溶解气体的若干组监测数据,对各组所述监测数据中的溶解气体含量进行分析,获取对应的故障类型标签,对各组所述监测数据进行归一化处理后与对应的故障类型标签组成目标数据集;
第一训练模块,用于将所述目标数据集划分为第一训练集和第一验证集,由所述第一训练集对LSTM诊断模型进行训练,由所述第一验证集对训练后的LSTM诊断模型进行验证;
第二训练模块,用于对所述目标数据集中的各组数据进行图像化处理,得到图像数据集,将所述图像数据集划分为第二训练集和第二验证集,由所述第二训练集对CNN诊断模型进行训练,由所述第二验证集对训练后的CNN诊断模型进行验证;
诊断模块,用于分别对训练后的LSTM诊断模型及训练后的CNN诊断模型的softmax层的输出进行深度并行融合,根据最大置信原则,输出最终的诊断结果。
10.一种计算机可读存储介质,其上存储有程序指令,其特征在于,所述程序指令被处理器执行时实现如权利要求1至8任一所述的变压器油中溶解气体的深度并行故障诊断方法。
CN202010134616.3A 2020-03-02 2020-03-02 变压器油中溶解气体的深度并行故障诊断方法及系统 Active CN111337768B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010134616.3A CN111337768B (zh) 2020-03-02 2020-03-02 变压器油中溶解气体的深度并行故障诊断方法及系统
US17/160,415 US11656298B2 (en) 2020-03-02 2021-01-28 Deep parallel fault diagnosis method and system for dissolved gas in transformer oil

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010134616.3A CN111337768B (zh) 2020-03-02 2020-03-02 变压器油中溶解气体的深度并行故障诊断方法及系统

Publications (2)

Publication Number Publication Date
CN111337768A true CN111337768A (zh) 2020-06-26
CN111337768B CN111337768B (zh) 2021-01-19

Family

ID=71183831

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010134616.3A Active CN111337768B (zh) 2020-03-02 2020-03-02 变压器油中溶解气体的深度并行故障诊断方法及系统

Country Status (2)

Country Link
US (1) US11656298B2 (zh)
CN (1) CN111337768B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111738436A (zh) * 2020-06-28 2020-10-02 电子科技大学中山学院 一种模型蒸馏方法、装置、电子设备及存储介质
CN111950526A (zh) * 2020-09-01 2020-11-17 国网河北省电力有限公司检修分公司 基于深度学习断路器储能机构故障诊断方法
CN112163619A (zh) * 2020-09-27 2021-01-01 北华大学 一种基于二维张量的变压器故障诊断方法
CN112232396A (zh) * 2020-10-08 2021-01-15 西北工业大学 一种基于lstm与d-s证据理论的舰船编队意图融合识别方法
CN112686093A (zh) * 2020-12-02 2021-04-20 重庆邮电大学 一种基于ds证据理论的融合局部放电类型识别方法
CN112749509A (zh) * 2020-12-30 2021-05-04 西华大学 一种基于lstm神经网络的智能变电站故障诊断方法
CN113761804A (zh) * 2021-09-13 2021-12-07 国网江苏省电力有限公司电力科学研究院 一种变压器状态诊断方法、计算机设备及存储介质
CN113782113A (zh) * 2021-09-17 2021-12-10 黄河水利职业技术学院 一种基于深度残差网络下的变压器油中气体故障识别方法
CN114549883A (zh) * 2022-02-24 2022-05-27 北京百度网讯科技有限公司 图像处理方法、深度学习模型的训练方法、装置和设备
CN115201640A (zh) * 2022-06-21 2022-10-18 国网四川省电力公司经济技术研究院 基于fi-cnn的变压器故障诊断方法

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111753895A (zh) * 2020-06-12 2020-10-09 北京小米松果电子有限公司 数据处理方法、装置及存储介质
CN112115417B (zh) * 2020-08-10 2024-02-27 贵州电网有限责任公司 一种基于概率密度聚类的变压器实时故障概率计算方法
CN113889198B (zh) * 2021-09-24 2024-09-24 国网宁夏电力有限公司电力科学研究院 一种基于油色谱时频域信息和残差注意力网络的变压器故障诊断方法及设备
CN113870052A (zh) * 2021-09-28 2021-12-31 国网福建省电力有限公司 基于多输入lstm-cnn的工作票安全措施识别方法及终端
CN113917370B (zh) * 2021-09-29 2023-12-29 浙江工业大学 一种基于油中溶解气体小样本数据的变压器故障诊断方法
CN113959956A (zh) * 2021-10-21 2022-01-21 河北卫讯电力自动化设备有限公司 一种变压器油中溶解气体双腔室光声光谱监测系统
CN114048815B (zh) * 2021-11-12 2023-11-28 国网江苏省电力有限公司常州供电分公司 基于厂站侧的电网运行信息感知系统、感知方法
CN113963757B (zh) * 2021-11-15 2024-08-06 国网山西省电力公司晋城供电公司 基于气体关系和图神经网络的充油电气设备故障诊断方法
CN114117909A (zh) * 2021-11-19 2022-03-01 河南九域恩湃电力技术有限公司 基于统计分布模型的变压器油溶解气体健康状态预测方法
CN114239703B (zh) * 2021-12-01 2024-10-15 湖南大学 主动配电系统故障诊断方法、系统、设备及存储介质
CN114167180B (zh) * 2021-12-03 2023-10-03 国网山西省电力公司晋城供电公司 一种基于图注意力神经网络的充油电气设备故障诊断方法
CN114325493B (zh) * 2021-12-07 2024-04-02 山东和兑智能科技有限公司 基于多维度关联和综合诊断的变压器状态评价方法
CN114596445B (zh) * 2022-01-24 2024-09-10 南京富岛油气智控科技有限公司 一种提升采油机故障诊断精度方法
CN114755529A (zh) * 2022-04-06 2022-07-15 重庆大学 一种基于深度学习的多特征融合单相接地故障类型辨识方法
CN114895222B (zh) * 2022-04-29 2024-08-20 宁波晨江电力集团有限公司 一种识别变压器各类故障及多重故障的诊断方法
CN114819102A (zh) * 2022-05-19 2022-07-29 东南大学溧阳研究院 基于gru的空调设备故障诊断方法
CN114757309B (zh) * 2022-06-16 2022-09-02 青岛理工大学 多物理场监测数据协同融合的工程灾害预警方法及系统
CN115078618B (zh) * 2022-07-13 2023-11-07 广东电网有限责任公司 一种变压器油色谱故障识别方法及相关装置
CN115407753B (zh) * 2022-08-18 2024-02-09 广东元梦泽技术服务有限公司 一种多变量加权集成学习的工业故障诊断方法
CN115753103A (zh) * 2022-08-30 2023-03-07 山东建筑大学 一种基于标准自学习数据增强的故障诊断方法及系统
CN115169506A (zh) * 2022-09-06 2022-10-11 中铁第四勘察设计院集团有限公司 一种供变电关键设备故障快速诊断方法及系统
CN116298674B (zh) * 2023-02-14 2023-12-26 四川轻化工大学 基于MobileNetV3的同杆双回线路故障选相
CN116310760B (zh) * 2023-04-06 2023-09-22 河南禹宏实业有限公司 基于机器视觉的智能水利监测系统
CN117350328B (zh) * 2023-09-11 2024-06-25 江南大学 一种基于lstm-cvae的乳酸菌发酵过程故障检测与诊断方法
CN117349786B (zh) * 2023-09-11 2024-08-09 国网湖北省电力有限公司宜昌供电公司 基于数据均衡的证据融合变压器故障诊断方法
CN117708696A (zh) * 2023-11-29 2024-03-15 国网湖北省电力有限公司超高压公司 一种基于LGEO-gcForest的油浸式变压器故障诊断方法
CN117554302B (zh) * 2024-01-09 2024-04-02 杉木(深圳)生物科技有限公司 液体检测方法、装置、设备及存储介质
CN117725529B (zh) * 2024-02-18 2024-05-24 南京邮电大学 基于多模态自注意力机制的变压器故障诊断方法
CN118425691B (zh) * 2024-07-05 2024-09-24 国网浙江省电力有限公司宁波市鄞州区供电公司 一种输电线路故障定位方法、装置、设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108037378A (zh) * 2017-10-26 2018-05-15 上海交通大学 基于长短时记忆网络的变压器运行状态预测方法及系统
CN109242205A (zh) * 2018-10-09 2019-01-18 国网上海市电力公司 一种变压器状态预测方法及系统
CN109885948A (zh) * 2019-02-28 2019-06-14 沈阳天眼智云信息科技有限公司 电力变压器故障诊断方法
CN110441500A (zh) * 2019-08-15 2019-11-12 昆明理工大学 一种基于长短期记忆网络的变压器油中溶解气体浓度预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108037378A (zh) * 2017-10-26 2018-05-15 上海交通大学 基于长短时记忆网络的变压器运行状态预测方法及系统
CN109242205A (zh) * 2018-10-09 2019-01-18 国网上海市电力公司 一种变压器状态预测方法及系统
CN109885948A (zh) * 2019-02-28 2019-06-14 沈阳天眼智云信息科技有限公司 电力变压器故障诊断方法
CN110441500A (zh) * 2019-08-15 2019-11-12 昆明理工大学 一种基于长短期记忆网络的变压器油中溶解气体浓度预测方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111738436A (zh) * 2020-06-28 2020-10-02 电子科技大学中山学院 一种模型蒸馏方法、装置、电子设备及存储介质
CN111738436B (zh) * 2020-06-28 2023-07-18 电子科技大学中山学院 一种模型蒸馏方法、装置、电子设备及存储介质
CN111950526A (zh) * 2020-09-01 2020-11-17 国网河北省电力有限公司检修分公司 基于深度学习断路器储能机构故障诊断方法
CN112163619A (zh) * 2020-09-27 2021-01-01 北华大学 一种基于二维张量的变压器故障诊断方法
CN112232396B (zh) * 2020-10-08 2023-01-24 西北工业大学 一种基于lstm与d-s证据理论的舰船编队意图融合识别方法
CN112232396A (zh) * 2020-10-08 2021-01-15 西北工业大学 一种基于lstm与d-s证据理论的舰船编队意图融合识别方法
CN112686093A (zh) * 2020-12-02 2021-04-20 重庆邮电大学 一种基于ds证据理论的融合局部放电类型识别方法
CN112749509A (zh) * 2020-12-30 2021-05-04 西华大学 一种基于lstm神经网络的智能变电站故障诊断方法
CN113761804A (zh) * 2021-09-13 2021-12-07 国网江苏省电力有限公司电力科学研究院 一种变压器状态诊断方法、计算机设备及存储介质
CN113782113A (zh) * 2021-09-17 2021-12-10 黄河水利职业技术学院 一种基于深度残差网络下的变压器油中气体故障识别方法
CN114549883A (zh) * 2022-02-24 2022-05-27 北京百度网讯科技有限公司 图像处理方法、深度学习模型的训练方法、装置和设备
CN114549883B (zh) * 2022-02-24 2023-09-05 北京百度网讯科技有限公司 图像处理方法、深度学习模型的训练方法、装置和设备
CN115201640A (zh) * 2022-06-21 2022-10-18 国网四川省电力公司经济技术研究院 基于fi-cnn的变压器故障诊断方法

Also Published As

Publication number Publication date
US11656298B2 (en) 2023-05-23
US20210278478A1 (en) 2021-09-09
CN111337768B (zh) 2021-01-19

Similar Documents

Publication Publication Date Title
CN111337768B (zh) 变压器油中溶解气体的深度并行故障诊断方法及系统
CN113283027B (zh) 一种基于知识图谱和图神经网络的机械故障诊断方法
CN110929918B (zh) 一种基于CNN和LightGBM的10kV馈线故障预测方法
CN109146246B (zh) 一种基于自动编码器和贝叶斯网络的故障检测方法
CN112101480B (zh) 一种多变量聚类与融合的时间序列组合预测方法
CN110929847A (zh) 一种基于深度卷积神经网络的换流变压器故障诊断方法
Fu et al. Spatiotemporal attention networks for wind power forecasting
CN117273422B (zh) 用于数字化生产的供应链协同管控方法及系统
Herman et al. Deep learning for oil palm fruit ripeness classification with DenseNet
CN112147432A (zh) 基于注意力机制的BiLSTM模块、变压器状态诊断方法和系统
Lin et al. CAM-UNET: class activation MAP guided UNET with feedback refinement for defect segmentation
CN112330050A (zh) 一种基于双层XGBoost考虑多特征的电力系统负荷预测方法
CN116562114A (zh) 一种基于图卷积神经网络的电力变压器故障诊断方法
Xu et al. Global attention mechanism based deep learning for remaining useful life prediction of aero-engine
CN116842337A (zh) 基于LightGBM优选特征与COA-CNN模型的变压器故障诊断方法
CN114169091A (zh) 工程机械部件剩余寿命的预测模型建立方法及预测方法
CN117494037A (zh) 一种基于变权vae与双通道特征融合的变压器故障诊断方法
CN117150383A (zh) 一种ShuffleDarkNet37-SE的新能源汽车动力电池故障分类方法
CN114547976B (zh) 基于金字塔变分自编码器的多采样率数据软测量建模方法
CN117235565A (zh) 一种变压器故障诊断模型构建方法和装置
Zhong et al. Fault diagnosis of Marine diesel engine based on deep belief network
CN117150232A (zh) 一种大模型非时序训练数据质量评价方法
CN112001530A (zh) 一种变压器油色谱在线监测装置预测性维护方法及系统
CN114383845B (zh) 一种基于嵌入式零样本学习模型的轴承复合故障诊断方法
CN113780405A (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