CN114089220B - 直流充电桩整流器功率管开路故障诊断及故障点定位方法 - Google Patents
直流充电桩整流器功率管开路故障诊断及故障点定位方法 Download PDFInfo
- Publication number
- CN114089220B CN114089220B CN202111407262.6A CN202111407262A CN114089220B CN 114089220 B CN114089220 B CN 114089220B CN 202111407262 A CN202111407262 A CN 202111407262A CN 114089220 B CN114089220 B CN 114089220B
- Authority
- CN
- China
- Prior art keywords
- fault
- phase current
- power tube
- open
- rectifier
- 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
Links
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/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
- G01R31/54—Testing for continuity
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Rectifiers (AREA)
Abstract
本申请提供一种直流充电桩整流器功率管开路故障诊断及故障点定位方法,所述方法包括如下步骤:S1:采集直流充电桩整流器交流侧相电流数据;S2:确定相电流数据窗口的瞬时频率的时间序列;S3:确定整流器桥臂故障特征量Flj,若整流器桥臂故障特征量Flj为1,存在故障,并进入下一步,若整流器桥臂故障特征量Flj为0,则无故障,并结束;S4:确定发生开路故障的相电流均值iave‑l;S5:确定开路故障功率管定位特征量Rl,若开路故障功率管定位特征量Rl=1,则表示功率管上管故障,若开路故障功率管定位特征量Rl=0,则表示功率管双管故障,若开路故障功率管定位特征量Rl=‑1,则表示功率管下管故障。本申请实现对整流器故障桥臂中开路故障功率管的定位。
Description
技术领域
本发明涉及充电桩故障诊断技术领域,尤其涉及一种直流充电桩整流器功率管开路故障诊断及故障点定位方法。
背景技术
电动汽车充电桩已成为重要的基础设施,直流充电桩是未来建设和发展重要方向,由于其变流器功率器件的脆弱性和控制的复杂性导致故障比例较高,是充电桩的薄弱环节。其中,功率器件发生开路故障时,具备故障特征不明显和不易发现特点。功率器件开路故障若不及时发现和检修,将对电网产生谐波影响,并诱发正常功率器件加速损坏,最终导致充电桩突发故障停机。然而,目前还没有有效易行的故障诊断技术。
因此,有必要研究直流充电桩整流器功率管开路故障诊断方法。为直流充电桩的智能运维提供技术支撑,对降低了直流充电桩运维费用和提高直流充电桩可靠性有重要研究意义。
发明内容
有鉴于此,本发明提供一种直流充电桩整流器功率管开路故障诊断及故障点定位方法,其特征在于:所述方法包括如下步骤:
S1:采集直流充电桩整流器交流侧相电流数据;
S2:基于交流侧相电流数据确定相电流数据窗口的瞬时频率的时间序列;
S3:基于相电流数据窗口的瞬时频率的时间序列确定整流器桥臂故障特征量Flj,若整流器桥臂故障特征量Flj为1,存在故障,并进入下一步,若整流器桥臂故障特征量Flj为0,则无故障,并结束;
S4:确定发生开路故障的相电流均值iave-l;
S5:基于故障相电流均值iave-l确定开路故障功率管定位特征量Rl,若开路故障功率管定位特征量Rl=1,则表示功率管上管故障,若开路故障功率管定位特征量Rl=0,则表示功率管双管故障,若开路故障功率管定位特征量Rl=-1,则表示功率管下管故障。
进一步,步骤S2包括如下步骤:
S21:将步骤S1采集的相电流数据进行希尔伯特变换获得相电流信号瞬时频率f(t),将相电流信号瞬时频率f(t)按采集时间先后顺序排列获得相电流瞬时频率的时间序列L,若时间序列L中包括n个时刻数据点,即L={i(1),i(2),i(3),…,i(n)},其中,L表示相电流时序数据,i(1)表示第1个时刻数据点,i(2)表示第2个时刻数据点,i(3)表示第3个时刻数据点,i(n)表示第n个时刻数据点,n表示相电流数据点的变量;
S22:将包括有m个时序数据点确定为一个数据窗口,确定{i(1),i(2),i(3),…,i(m)}为L的第一个数据窗口I0,接下来,每次将窗口按时间顺序在L上滑动k个数据,直至滑动s次(s≤((n-m)/k)+1),依次获得窗口I1、I2、…Is,其中,s表示第s个窗口,n表示相电流数据点的变量,m表示一个窗口的数据点的数量,k表示每次滑动的数据点数量,其中,n>>m>>k>0,且都为正整数;
S23:基于希尔伯特变换求取每个数据窗口的瞬时频率,进而确定各窗口瞬时频率时间序列{f0,f1,f2,…,fs};
S24:重复步骤S21-S23,获得l相相电流数据窗口的瞬时频率{fl0,fl1,fl2,…,fls},其中,l相相电流包括a相电流、b相电流和c相电流。
进一步,步骤S21中所述相电流信号瞬时频率f(t)采用如下方法确定:
其中,f(t)表示相电流信号瞬时频率,φ(t)为信号的相位信息。
进一步,所述步骤S3包括对相电流数据窗口的瞬时频率的时间序列进行归一化处理和确定整流器桥臂故障特征量Flj;
所述归一化处理采用如下方法确定:
其中,μlj表示归一化处理后的相电流数据窗口的瞬时频率,flj相电流数据窗口的瞬时频率,fN为额定频率;
所述整流器桥臂故障特征量Flj采用如下方法确定:
其中,k1为判断阈值I,μlj表示归一化处理后的相电流数据窗口的瞬时频率。
进一步,所述开路故障相电流均值iave-l采用如下方法确定:
其中,iave-l表示发生开路故障的相电流均值,q表示开路故障桥臂发生故障时刻前后的整数倍电流周期内的q个数据点,h表示数据点变量,il(h)表示相电流。
进一步,所述开路故障功率管定位特征量Rl采用如下方法确定:
其中,Rl表示开路故障功率管定位特征量,iave-l表示发生开路故障的相电流均值,k2为判断阈值Ⅱ。
本发明的有益技术效果:本发明提出的基于相电流瞬时频率的直流充电桩整流器功率管开路故障诊断方法可应用于直流充电桩整流器功率管的故障诊断中,仅需采用交流侧相电流数据,就可以实现直流充电桩整流器功率管开路故障诊断,可为直流充电桩的智能运维提供技术支撑,对降低直流充电桩运维费用和提高直流充电桩可靠性有重要应用价值。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1直流充电桩整流器功率管开路故障诊断流程图。
图2直流充电桩整流器电路结构图。
图3直流充电桩整流器功率管VT5开路故障诊断结果。
具体实施方式
以下结合说明书附图对本发明做出进一步的说明:
本发明提供一种直流充电桩整流器功率管开路故障诊断及故障点定位方法,其特征在于:所述方法包括如下步骤:如图1所示,
S1:采集直流充电桩整流器交流侧相电流数据;
S2:基于交流侧相电流数据确定相电流数据窗口的瞬时频率的时间序列;
S3:基于相电流数据窗口的瞬时频率的时间序列确定整流器桥臂故障特征量Flj,若整流器桥臂故障特征量Flj为1,存在故障,并进入下一步,若整流器桥臂故障特征量Flj为0,则无故障,并结束;
S4:确定发生开路故障的相电流均值iave-l;
S5:基于故障相电流均值iave-l确定开路故障功率管定位特征量Rl,若开路故障功率管定位特征量Rl=1,则表示功率管上管故障,若开路故障功率管定位特征量Rl=0,则表示功率管双管故障,若开路故障功率管定位特征量Rl=-1,则表示功率管下管故障。
上述技术方案本发明提出了基于相电流瞬时频率的直流充电桩整流器功率管开路故障诊断方法,采集直流充电桩整流器交流侧相电流信号数据,希尔伯特变换和滑动窗口均值的瞬时频率方法,获取相电流瞬时频率;其次,将瞬时频率归一化,通过构造的桥臂开路故障识别特征量,获得整流器开路故障桥臂;最后,利用相电流均值,通过故障桥臂功率管定位特征量,实现对整流器故障桥臂中开路故障功率管的定位。
步骤S2包括如下步骤:
S21:将步骤S1采集的相电流数据进行希尔伯特变换获得相电流信号瞬时频率f(t),将相电流信号瞬时频率f(t)按采集时间先后顺序排列获得相电流瞬时频率的时间序列L,若时间序列L中包括n个时刻数据点,即L={i(1),i(2),i(3),…,i(n)},其中,L表示相电流时序数据,i(1)表示第1个时刻数据点,i(2)表示第2个时刻数据点,i(3)表示第3个时刻数据点,i(n)表示第n个时刻数据点,n表示相电流数据点的变量;
S22:将包括有m个时序数据点确定为一个数据窗口,确定{i(1),i(2),i(3),…,i(m)}为L的第一个数据窗口I0,接下来,每次将窗口按时间顺序在L上滑动k个数据,直至滑动s次(s≤((n-m)/k)+1),依次获得窗口I1、I2、…Is,其中,s表示第s个窗口,n表示相电流数据点的变量,m表示一个窗口的数据点的数量,k表示每次滑动的数据点数量,其中,n>>m>>k>0,且都为正整数;
S23:基于希尔伯特变换求取每个数据窗口的瞬时频率,进而确定各窗口瞬时频率时间序列{f0,f1,f2,…,fs};
S24:重复步骤S21-S23,获得l相相电流数据窗口的瞬时频率{fl0,fl1,fl2,…,fls},其中,l相相电流包括a相电流、b相电流和c相电流。
步骤S21中所述相电流信号瞬时频率f(t)采用如下方法确定:
其中,f(t)表示相电流信号瞬时频率,φ(t)为信号的相位信息。
所述步骤S3包括对相电流数据窗口的瞬时频率的时间序列进行归一化处理和确定整流器桥臂故障特征量Flj;
所述归一化处理采用如下方法确定:
其中,μlj表示归一化处理后的相电流数据窗口的瞬时频率,flj相电流数据窗口的瞬时频率,fN为额定频率;
所述整流器桥臂故障特征量Flj采用如下方法确定:
其中,k1为判断阈值I,μlj表示归一化处理后的相电流数据窗口的瞬时频率。判断阈值I通过经验设定,在本实施例中k1为0.9。
在本实施例中,所述开路故障相电流均值iave-l采用如下方法确定:
其中,iave-l表示发生开路故障的相电流均值,q表示开路故障桥臂发生故障时刻前后的整数倍电流周期内的q个数据点,h表示数据点变量,il(h)表示相电流。
在本实施例中,所述开路故障功率管定位特征量Rl采用如下方法确定:
其中,Rl表示开路故障功率管定位特征量,iave-l表示发生开路故障的相电流均值,k2为判断阈值Ⅱ。判断阈值Ⅱ通过经验设定,在本实施例中k2取0.1。在本实施例中根据经验设定的判断阈值I和判断阈值Ⅱ均为本领域技术人员通过大量实验后确定的经验数值。
以实际直流充电桩整流器为例,其电路图如图2所示。采集整流器交流侧某相电流信号数据L,利用希尔伯特变换和滑动窗口的瞬时频率计算方法,得到通过滑动窗口的瞬时频率,归一化瞬时频率,通过故障桥臂开路故障识别特征量Flj和相电流均值iave-l的开路故障功率管定位特征量Rl,实现直流充电桩功率整流器功率管开路故障识别和定位。接下来,以图2中C相桥臂功率管VT5发生开路故障为例,说明本发明的有效性。
瞬时频率计算部分。采集直流充电桩整流器交流侧相电流信号数据,采用希尔伯特变换和滑动窗口的瞬时频率计算方法,获得相电流瞬时频率。例如,获得通过滑动窗口的瞬时频率序列。过程为:采集整流器交流侧某相相电流信号数据L,利用希尔伯特变换和滑动窗口的瞬时频率计算方法,得到通过滑动窗口的瞬时频率。
开路故障识别部分。对相电流瞬时频率归一化处理,通过桥臂开路故障识别特征量,获得整流器开路故障桥臂。例如,当图2中C相桥臂功率管VT5发生开路故障时,利用所采集直流充电桩整流器交流侧相电流信号数据,采用希尔伯特变换和滑动窗口的瞬时频率计算方法,获得相电流瞬时频率,再通过式(3)计算得到开路故障识别特征量Flj,见图3所示,在时间为1秒时,由1秒时的滑动窗口数据,计算得到Fa=0、Fb=0、Fc=1,由此便可以识别出整流器中含有功率管VT2和VT5的C相桥臂在1秒时出现了开路故障。
开路故障定位部分。若存在整流器开路故障桥臂,计算该相相电流均值,通过故障桥臂功率管定位特征量Rc,实现对整流器故障桥臂中开路故障功率管的定位。例如,由S2已经识别出整流器中C相桥臂出现了开路故障,根据式(4)和式(5)先后计算获得故障桥臂电流均值和开路故障桥臂功率管定位特征量Rc=1,从时间1秒时开始,C相上桥臂功率管VT5出现了开路故障,有效了实现了开路故障的定位。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种直流充电桩整流器功率管开路故障诊断及故障点定位方法,其特征在于:所述方法包括如下步骤:
S1:采集直流充电桩整流器交流侧相电流数据;
S2:基于交流侧相电流数据确定相电流数据窗口的瞬时频率的时间序列;
S3:基于相电流数据窗口的瞬时频率的时间序列确定整流器桥臂故障特征量Flj,若整流器桥臂故障特征量Flj为1,存在故障,并进入下一步,若整流器桥臂故障特征量Flj为0,则无故障,并结束;
S4:确定发生开路故障的相电流均值iave-l;
S5:基于故障相电流均值iave-l确定开路故障功率管定位特征量Rl,若开路故障功率管定位特征量Rl=1,则表示功率管上管故障,若开路故障功率管定位特征量Rl=0,则表示功率管双管故障,若开路故障功率管定位特征量Rl=-1,则表示功率管下管故障;
步骤S2包括如下步骤:
S21:将步骤S1采集的相电流数据进行希尔伯特变换获得相电流信号瞬时频率f(t),将相电流信号瞬时频率f(t)按采集时间先后顺序排列获得相电流瞬时频率的时间序列L,若时间序列L中包括n个时刻数据点,即L={i(1),i(2),i(3),…,i(n)},其中,L表示相电流时序数据,i(1)表示第1个时刻数据点,i(2)表示第2个时刻数据点,i(3)表示第3个时刻数据点,i(n)表示第n个时刻数据点,n表示相电流数据点的变量;
S22:将包括有m个时序数据点确定为一个数据窗口,确定{i(1),i(2),i(3),…,i(m)}为L的第一个数据窗口I0,接下来,每次将窗口按时间顺序在L上滑动k个数据,直至滑动s次(s≤((n-m)/k)+1),依次获得窗口I1、I2、…Is,其中,s表示第s个窗口,n表示相电流数据点的变量,m表示一个窗口的数据点的数量,k表示每次滑动的数据点数量,其中,n>>m>>k>0,且都为正整数;
S23:基于希尔伯特变换求取每个数据窗口的瞬时频率,进而确定各窗口瞬时频率时间序列{f0,f1,f2,…,fs};
S24:重复步骤S21-S23,获得l相相电流数据窗口的瞬时频率{fl0,fl1,fl2,…,fls},其中,l相相电流包括a相电流、b相电流和c相电流;
步骤S21中所述相电流信号瞬时频率f(t)采用如下方法确定:
其中,f(t)表示相电流信号瞬时频率,φ(t)为信号的相位信息;
所述步骤S3包括对相电流数据窗口的瞬时频率的时间序列进行归一化处理和确定整流器桥臂故障特征量Flj;
所述归一化处理采用如下方法确定:
其中,μlj表示归一化处理后的相电流数据窗口的瞬时频率,flj相电流数据窗口的瞬时频率,fN为额定频率;
所述整流器桥臂故障特征量Flj采用如下方法确定:
其中,k1为判断阈值I,μlj表示归一化处理后的相电流数据窗口的瞬时频率。
2.根据权利要求1所述直流充电桩整流器功率管开路故障诊断及故障点定位方法,其特征在于:所述开路故障相电流均值iave-l采用如下方法确定:
其中,iave-l表示发生开路故障的相电流均值,q表示开路故障桥臂发生故障时刻前后的整数倍电流周期内的q个数据点,h表示数据点变量,il(h)表示相电流。
3.根据权利要求2所述直流充电桩整流器功率管开路故障诊断及故障点定位方法,其特征在于:所述开路故障功率管定位特征量Rl采用如下方法确定:
其中,Rl表示开路故障功率管定位特征量,iave-l表示发生开路故障的相电流均值,k2为判断阈值Ⅱ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111407262.6A CN114089220B (zh) | 2021-11-24 | 2021-11-24 | 直流充电桩整流器功率管开路故障诊断及故障点定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111407262.6A CN114089220B (zh) | 2021-11-24 | 2021-11-24 | 直流充电桩整流器功率管开路故障诊断及故障点定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114089220A CN114089220A (zh) | 2022-02-25 |
CN114089220B true CN114089220B (zh) | 2023-08-25 |
Family
ID=80304170
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111407262.6A Active CN114089220B (zh) | 2021-11-24 | 2021-11-24 | 直流充电桩整流器功率管开路故障诊断及故障点定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114089220B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106019056A (zh) * | 2016-06-23 | 2016-10-12 | 华中科技大学 | 一种三相整流器功率管开路故障实时检测方法 |
CN106680647A (zh) * | 2016-12-06 | 2017-05-17 | 西南交通大学 | 一种动车组逆变器igbt单管断路故障诊断方法 |
CN107102247A (zh) * | 2017-04-18 | 2017-08-29 | 西南交通大学 | 一种crh3型动车组牵引逆变器igbt开路故障诊断方法 |
CN110244169A (zh) * | 2019-07-17 | 2019-09-17 | 西南交通大学 | 一种双重化脉冲整流器功率管多管开路故障快速诊断方法 |
CN111077471A (zh) * | 2019-12-30 | 2020-04-28 | 合肥工业大学 | 基于瞬时频率的npc三电平逆变器开路故障诊断方法 |
CN111562517A (zh) * | 2020-06-10 | 2020-08-21 | 合肥工业大学 | 一种npc三电平逆变器开关管开路故障诊断方法 |
CN112083353A (zh) * | 2020-07-22 | 2020-12-15 | 国网上海市电力公司 | 一种基于开关模态特征对变流器开路故障进行检测的方法及系统 |
CN114879087A (zh) * | 2022-05-12 | 2022-08-09 | 河北工业大学 | 一种用于vienna整流器开路故障综合辨识的方法 |
CN115774218A (zh) * | 2022-11-16 | 2023-03-10 | 重庆理工大学 | 基于Spearman的直流充电桩整流器开路故障诊断方法 |
-
2021
- 2021-11-24 CN CN202111407262.6A patent/CN114089220B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106019056A (zh) * | 2016-06-23 | 2016-10-12 | 华中科技大学 | 一种三相整流器功率管开路故障实时检测方法 |
CN106680647A (zh) * | 2016-12-06 | 2017-05-17 | 西南交通大学 | 一种动车组逆变器igbt单管断路故障诊断方法 |
CN107102247A (zh) * | 2017-04-18 | 2017-08-29 | 西南交通大学 | 一种crh3型动车组牵引逆变器igbt开路故障诊断方法 |
CN110244169A (zh) * | 2019-07-17 | 2019-09-17 | 西南交通大学 | 一种双重化脉冲整流器功率管多管开路故障快速诊断方法 |
CN111077471A (zh) * | 2019-12-30 | 2020-04-28 | 合肥工业大学 | 基于瞬时频率的npc三电平逆变器开路故障诊断方法 |
CN111562517A (zh) * | 2020-06-10 | 2020-08-21 | 合肥工业大学 | 一种npc三电平逆变器开关管开路故障诊断方法 |
CN112083353A (zh) * | 2020-07-22 | 2020-12-15 | 国网上海市电力公司 | 一种基于开关模态特征对变流器开路故障进行检测的方法及系统 |
CN114879087A (zh) * | 2022-05-12 | 2022-08-09 | 河北工业大学 | 一种用于vienna整流器开路故障综合辨识的方法 |
CN115774218A (zh) * | 2022-11-16 | 2023-03-10 | 重庆理工大学 | 基于Spearman的直流充电桩整流器开路故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114089220A (zh) | 2022-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021184820A1 (zh) | 一种两电平三相电压源逆变器桥臂开路故障诊断方法 | |
CN108196154B (zh) | 航空三级式同步电机旋转整流器故障检测及故障定位方法 | |
CN106841949B (zh) | 三相异步交流电机定子绝缘在线监测方法和装置 | |
CN103675544B (zh) | 基于优化算法的电力系统故障信号检测与波形识别方法 | |
CN106383289B (zh) | 开关磁阻电机功率变换器相关性分析故障诊断方法 | |
CN112162218B (zh) | 一种直驱永磁风电机组变流器多功率管开路故障诊断方法 | |
CN101344549B (zh) | 一种三相电路缺相检测的方法及其装置 | |
CN102393507A (zh) | 一种电机参数检测方法及电机参数检测装置 | |
CN102937672A (zh) | 一种用于并网设备的三相交流电相序检测方法及其检测装置 | |
CN109738778B (zh) | 逆变器开路诊断方法、装置、终端设备及计算机可读介质 | |
CN114089220B (zh) | 直流充电桩整流器功率管开路故障诊断及故障点定位方法 | |
CN114236374A (zh) | 一种整流器开路故障的实时诊断方法 | |
CN112345885A (zh) | 一种高压直流输电线路故障检测方法 | |
CN109298361B (zh) | 基于暂态录波型故障指示器的故障分析预处理方法 | |
Veerendra et al. | A novel fault‐detection methodology of proposed reduced switch MLI fed induction motor drive using discrete wavelet transforms | |
CN112198458B (zh) | 一种三相电压源逆变器开路故障实时检测方法和系统 | |
CN115774218A (zh) | 基于Spearman的直流充电桩整流器开路故障诊断方法 | |
CN114609435A (zh) | 一种电压暂降检测和分类辨识方法 | |
CN108919120A (zh) | 一种能源互联网功率变换器故障诊断方法 | |
CN109245057B (zh) | 输电线路时域全波形保护装置及相间突变量方向判断方法 | |
CN112731103A (zh) | 一种双级矩阵变换器的故障诊断方法 | |
CN117092554B (zh) | 逆变器耦合故障分析方法、装置、电子设备及存储介质 | |
CN113514776A (zh) | 一种异步电机的匝间短路故障诊断方法、装置、设备 | |
Kherif et al. | Fault Diagnosis Method for Cascaded H-bridge Multilevel Inverters under Switching Device Failure | |
Jiajia et al. | Research on fault diagnosis of DC charging pile power device based on wavelet packet and elman neural network |
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 |