CN111506990A - 一种基于1维电弧模型的环保型灭弧气体快速筛选方法 - Google Patents
一种基于1维电弧模型的环保型灭弧气体快速筛选方法 Download PDFInfo
- Publication number
- CN111506990A CN111506990A CN202010258358.XA CN202010258358A CN111506990A CN 111506990 A CN111506990 A CN 111506990A CN 202010258358 A CN202010258358 A CN 202010258358A CN 111506990 A CN111506990 A CN 111506990A
- Authority
- CN
- China
- Prior art keywords
- arc
- gas
- dimensional
- medium
- stage
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012216 screening Methods 0.000 title claims abstract description 21
- 238000011084 recovery Methods 0.000 claims abstract description 98
- 238000011156 evaluation Methods 0.000 claims abstract description 14
- 230000005855 radiation Effects 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 19
- 230000005684 electric field Effects 0.000 claims description 12
- 238000012423 maintenance Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 54
- 238000002474 experimental method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010891 electric arc Methods 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010791 quenching Methods 0.000 description 1
- 230000000171 quenching effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于1维电弧模型的环保型灭弧气体快速筛选方法,包括如下步骤:步骤1,计算候选的环保型气体的电弧等离子体参数,具体包括热力学参数、输运系数和辐射系数;步骤2,建立1维电弧模型,基于步骤1的电弧等离子体参数求解所述1维电弧模型;步骤3,在步骤2求解结果基础上将电弧熄灭过程划分为热恢复阶段、弧前介质恢复阶段和弧后介质恢复阶段三个阶段,定义并计算每个阶段对应的评价参数分别为热恢复率RT、弧前介质恢复率Rpre和弧后介质恢复率Rpost;步骤4,计算步骤3获得的三个阶段速率参数的调和平均数,以此作为气体灭弧性能的综合评价指标,该指标值越大表明气体灭弧性能越好,从而完成对候选环保型灭弧气体的快速筛选。
Description
技术领域
本发明属于电气工程和等离子体技术领域,特别涉及一种气体灭弧性能评价方法。
背景技术
电弧是断路器开断过程中出现的一种放电现象。良好的灭弧气体应具有高的热导率,以快速地从电弧中去除能量,并且应具有高的介电强度,以避免在出现过渡恢复电压时电弧重燃。现有的气体灭弧性能评价方法可以分为实验方法和数值方法。实验方法包括短路故障切换试验、电感耦合热等离子体技术与功率半导体开关技术等,数值方法包括基于2维和3维磁流体动力学模型的方法。实验的方法成本较高,而数值方法由于需要进行2维或3维的磁流体动力学仿真,计算复杂且非常耗时。考虑到工程中寻找或筛选潜在的灭弧气体尤其是环保型SF6替代气体,需要测试大量的候选气体或混合气体,还要不断地改变开断电流、气体压力等条件,而常规的实验方法和数值计算方法均不能满足廉价、快速评价气体灭弧性能的要求,尤其是无法有效应用于寻找潜在的环保型SF6替代气体。因此快速评价气体灭弧性能对初步寻找潜在的灭弧介质具有重要意义。
发明内容
本发明的目的,在于提供一种基于1维电弧模型的环保型灭弧气体快速筛选方法,能够快速从候选气体中筛选出灭弧性能优良的环保型灭弧气体。
为了达成上述目的,本发明的解决方案是:
一种基于1维电弧模型的环保型灭弧气体快速筛选方法,包括如下步骤:
步骤1,计算电弧等离子体参数,具体包括热力学参数、输运系数和辐射系数;
步骤2,建立1维电弧模型,基于步骤1的电弧等离子体参数求解所述1维电弧模型;
步骤3,在步骤2求解结果基础上将电弧熄灭过程划分为热恢复阶段、弧前介质恢复阶段和弧后介质恢复阶段三个阶段,并计算每个阶段的对应评价参数分别为热恢复率RT、弧前介质恢复率Rpre和弧后介质恢复率Rpost;
步骤4,计算步骤3获得的三阶段速率参数的调和平均数,以该调和平均数(平均速率)作为气体灭弧性能的综合评价指标,该指标值越大表明气体灭弧性能越好,从而完成对候选环保型灭弧气体的快速筛选。
上述步骤2的具体内容是:
步骤21,首先进行如下基本假设:电弧等离子体处于局部热力学平衡;电弧等离子体具有圆柱对称性;电弧等离子体中的压强不随时空变化;忽略电弧的轴向流动和轴向速度;
步骤22,基于步骤21的基本假设,建立1维电弧模型,包括:
其中,r是径向距离,ρ是质量密度,vr是径向速度分量,CP是定压比热,T是电弧温度,σ是电导率,k是热导率,i是电流,g是电弧电导,Erad是辐射产生的能量损失;
步骤24,基于有限体积法求解步骤22中的方程。
上述步骤24的具体过程是:
步骤241,首先设置初始电流I=200A,径向长度R=2.5mm;并确定方程收敛条件,设置为两个迭代步之间的温度T和径向速度vr的均方根值分别小于10-3和10-8;
步骤242,定义平均电弧温度,计算公式如下:
其中,R为径向长度。
上述步骤3中,热恢复阶段的划分依据是平均电弧温度Tavg大于4000K的阶段。
上述步骤3中,弧前介质恢复阶段的划分依据是平均温度Tavg小于4000K且电弧电导g大于10-7S·m的阶段。
上述步骤3中,弧后介质恢复阶段的划分依据弧前介质恢复阶段结束到平均电场强度Eavg升高50%的阶段。
上述平均电场强度Eavg的计算方法是:首先求解玻尔兹曼方程得到临界电场强度Ecr(r),然后根据下式计算平均临界电场强度Eavg:
其中,R为径向长度。
上述步骤3中,热恢复率RT的计算公式如下:
上述步骤3中,弧前介质恢复率Rpre的计算公式如下:
上述步骤3中,弧后介质恢复率Rpost的计算公式如下:
上述步骤4的具体内容是:
步骤41,所述步骤4中,三阶段平均速率Ravg的计算公式如下:
其中,RT、Rpre、Rpost分别为候选气体的热恢复率、弧前介质恢复率、弧后介质恢复率。
步骤42,以三阶段平均速率Ravg作为气体灭弧性能的综合评价指标,该指标值越大表明气体灭弧性能越好。根据综合评价指标,从候选气体中找到灭弧性能较好的气体,从而完成对候选环保型灭弧气体的快速筛选。
采用上述方案后,本发明利用1维电弧模型和玻尔兹曼方程分析分别得到电弧热恢复阶段、弧前介质恢复阶段和弧后介质恢复阶段的恢复速率,并以三阶段平均速率作为气体灭弧性能的综合评价指标,克服了现有方法需要昂贵实验或复杂耗时计算的缺点,能够快速从候选气体中筛选出灭弧性能优良的环保型灭弧气体。
附图说明
图1是本发明的流程图;
图2是CO2热恢复阶段和弧前介质恢复阶段示意图;
图3是CO2弧后介质恢复阶段示意图。
具体实施方式
以下将结合附图,对本发明的技术方案及有益效果进行详细说明。
本实施例中,将SF6,CO2,CF4和空气四种气体作为候选气体,分别比较这四种气体的灭弧性能。选择这四种气体作为研究对象,是因为这四种气体性质稳定,可以作为灭弧气体,而且有工程应用,所以将这四种气体作为研究对象,不仅具有代表性,而且比较结果对于工程上灭弧气体的选择有一定指导意义。
请参见图1,其示出了一种基于1维电弧模型的气体灭弧性能评价方法的流程图,该方法包括以下步骤:
(1)计算电弧等离子体热力学参数,包括质量密度、比焓、定压比热;计算电弧等离子体输运系数,包括热导率、电导率;计算电弧等离子体辐射系数;
(2)将电弧等离子体参数输入到1维电弧模型中进行求解;
(21)1维电弧模型基于如下假设而建立:电弧等离子体处于局部热力学平衡(LTE)状态;电弧等离子体具有圆柱对称性;电弧等离子体中的压强不随时空变化;忽略电弧的轴向流动和轴向速度。
(22)在步骤(21)的基本假设基础上建立质量守恒方程、能量守恒方程和欧姆定律方程:
其中,r是径向距离,ρ是质量密度,vr是径向速度分量,CP是定压比热,T是电弧温度,σ是电导率,k是热导率,I是电流,g是电弧电导,Erad是辐射产生的能量损失。
(24)基于有限体积法求解步骤(22)中的方程;
设置初始电流I=200A,径向长度R=2.5mm。
(241)确定方程收敛条件,方程收敛条件设置为两个迭代步之间的温度T和径向速度vr的均方根值分别小于10-3和10-8;
(242)计算平均电弧温度,计算公式如下:
(3)在1维电弧模型求解结果基础上将电弧熄灭过程划分为三个阶段,包括热恢复阶段、弧前介质恢复阶段和弧后介质恢复阶段;
(31)用平均电弧温度Tavg表征热恢复阶段,划分依据为从平均电弧温度Tavg大于4000K的阶段;
(32)用电弧电导g和平均温度Tavg表征弧前介质恢复阶段,划分依据为平均温度Tavg小于4000K并且电弧电导g大于10-7S·m的阶段;
图2以CO2为例给出了气体热恢复阶段和弧前介质恢复阶段划分示意图,其中Stage1是热恢复阶段,Stage2是弧前介质恢复阶段。因此,根据计算可以获得SF6、CO2、CF4和空气的热恢复阶段持续时间分别为190μs、340μs、246μs和447μs,这四种气体的弧前介质恢复持续时间分别为192μs、231μs、456μs和665μs。
在弧后介质恢复阶段,由于电弧等离子体不再处于局部热力学平衡状态,此时1维电弧模型不能够完整描述电弧等离子体的状态,本实施例通过求解玻尔兹曼方程解决该问题。
(33)用平均电场强度Eavg表征弧后介质恢复阶段;
(331)求解玻尔兹曼方程得到临界电场强度Ecr(r);
(332)计算平均临界电场强度Eavg,计算公式如下:
(333)弧后介质恢复阶段为弧前介质恢复阶段结束到平均电场强度Eavg升高50%的阶段;
图3以CO2为例给出了气体弧后介质恢复阶段划分示意图,其中Stage3是弧后介质恢复阶段。因此,根据计算可以获得SF6、CO2、CF4和空气的弧后介质恢复阶段持续时间分别为588μs、2829μs、1148μs和768μs。
(4)计算得到每个阶段的对应评价参数分别为热恢复率RT、弧前介质恢复率Rpre和弧后介质恢复率Rpost;
(41)计算表征电弧热恢复阶段的参数热恢复率RT,计算公式如下:
(42)计算表征电弧弧前介质恢复阶段的弧前介质恢复率Rpre,计算公式如下:
(43)计算表征电弧弧后介质恢复阶段的弧后介质恢复率Rpost,计算公式如下:
(5)根据上一步获得的三阶段速率参数计算平均恢复速率,以平均速率作为气体灭弧性能的综合评价指标,该指标值越大表明气体灭弧性能越好,从而完成对候选环保型灭弧气体的快速筛选。平均恢复速率的计算公式如下:
其中,RT、Rpre、Rpost分别为候选气体的热恢复率、弧前介质恢复率、弧后介质恢复率。
表1四种气体灭弧特征参数
从表1中可知,四种气体热恢复能力排序为:SF6>CF4>CO2>空气;弧前介质恢复能力排序为:SF6>CO2>CF4>空气;弧后介质恢复能力排序为:SF6>空气>CF4>CO2;综合灭弧能力排序为:SF6>CF4>空气>CO2。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (10)
1.一种基于1维电弧模型的环保型灭弧气体快速筛选方法,其特征在于包括如下步骤:
步骤1,计算候选气体的电弧等离子体参数,具体包括热力学参数、输运系数和辐射系数;
步骤2,建立1维电弧模型,基于步骤1的电弧等离子体参数求解所述1维电弧模型;
步骤3,在步骤2求解结果基础上将电弧熄灭过程划分为热恢复阶段、弧前介质恢复阶段和弧后介质恢复阶段三个阶段,并计算每个阶段的对应评价参数,分别为热恢复率RT、弧前介质恢复率Rpre和弧后介质恢复率Rpost;
步骤4,计算步骤3获得的三个阶段评价参数的调和平均数,以调和平均数作为气体灭弧性能的综合评价指标,该指标值越大表明气体灭弧性能越好,从而完成对候选环保型灭弧气体的快速筛选。
2.如权利要求1所述的基于1维电弧模型的环保型灭弧气体快速筛选方法,其特征在于:所述步骤2的具体内容是:
步骤21,首先进行如下基本假设:电弧等离子体处于局部热力学平衡;电弧等离子体具有圆柱对称性;电弧等离子体中的压强不随时空变化;忽略电弧的轴向流动和轴向速度;
步骤22,基于步骤21的基本假设,建立1维电弧模型,包括:
其中,r是径向距离,ρ是质量密度,vr是径向速度分量,CP是定压比热,T是电弧温度,σ是电导率,k是热导率,i是电流,g是电弧电导,Erad是辐射产生的能量损失;
步骤24,基于有限体积法求解步骤22中的方程。
4.如权利要求1所述的基于1维电弧模型的环保型灭弧气体快速筛选方法,其特征在于:所述步骤3中,热恢复阶段的划分依据是平均电弧温度Tavg大于4000K的阶段。
5.如权利要求1所述的基于1维电弧模型的环保型灭弧气体快速筛选方法,其特征在于:所述步骤3中,弧前介质恢复阶段的划分依据是平均温度Tavg小于4000K且电弧电导g大于10-7S·m的阶段,弧后介质恢复阶段的划分依据弧前介质恢复阶段结束到平均电场强度Eavg升高50%的阶段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010258358.XA CN111506990B (zh) | 2020-04-03 | 2020-04-03 | 一种基于1维电弧模型的环保型灭弧气体快速筛选方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010258358.XA CN111506990B (zh) | 2020-04-03 | 2020-04-03 | 一种基于1维电弧模型的环保型灭弧气体快速筛选方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111506990A true CN111506990A (zh) | 2020-08-07 |
CN111506990B CN111506990B (zh) | 2023-07-14 |
Family
ID=71863953
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010258358.XA Active CN111506990B (zh) | 2020-04-03 | 2020-04-03 | 一种基于1维电弧模型的环保型灭弧气体快速筛选方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111506990B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114487009A (zh) * | 2022-02-09 | 2022-05-13 | 广东电网有限责任公司 | 一种气体灭弧性能的判断方法及*** |
CN116150995A (zh) * | 2023-02-21 | 2023-05-23 | 东南大学 | 一种开关电弧模型的快速仿真方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253293A (zh) * | 2011-05-30 | 2011-11-23 | 南京航空航天大学 | 一种直流电弧故障检测方法及装置 |
CN105223427A (zh) * | 2014-06-17 | 2016-01-06 | 西门子公司 | 故障电弧的检测方法和检测装置 |
CN108549021A (zh) * | 2018-06-21 | 2018-09-18 | 沈阳工业大学 | 一种用于环保型气体灭弧性能研究的试验装置及方法 |
CN109255170A (zh) * | 2018-08-28 | 2019-01-22 | 电子科技大学 | 一种基于磁流体动力学仿真的电弧电压梯度建模方法 |
CN110188450A (zh) * | 2019-05-27 | 2019-08-30 | 三峡大学 | 基于电弧能量热等效的开关柜内部短路燃弧压力升计算方法 |
-
2020
- 2020-04-03 CN CN202010258358.XA patent/CN111506990B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253293A (zh) * | 2011-05-30 | 2011-11-23 | 南京航空航天大学 | 一种直流电弧故障检测方法及装置 |
CN105223427A (zh) * | 2014-06-17 | 2016-01-06 | 西门子公司 | 故障电弧的检测方法和检测装置 |
CN108549021A (zh) * | 2018-06-21 | 2018-09-18 | 沈阳工业大学 | 一种用于环保型气体灭弧性能研究的试验装置及方法 |
CN109255170A (zh) * | 2018-08-28 | 2019-01-22 | 电子科技大学 | 一种基于磁流体动力学仿真的电弧电压梯度建模方法 |
CN110188450A (zh) * | 2019-05-27 | 2019-08-30 | 三峡大学 | 基于电弧能量热等效的开关柜内部短路燃弧压力升计算方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114487009A (zh) * | 2022-02-09 | 2022-05-13 | 广东电网有限责任公司 | 一种气体灭弧性能的判断方法及*** |
CN116150995A (zh) * | 2023-02-21 | 2023-05-23 | 东南大学 | 一种开关电弧模型的快速仿真方法 |
CN116150995B (zh) * | 2023-02-21 | 2023-07-25 | 东南大学 | 一种开关电弧模型的快速仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111506990B (zh) | 2023-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lindmayer et al. | The process of arc splitting between metal plates in low voltage arc chutes | |
CN111506990B (zh) | 一种基于1维电弧模型的环保型灭弧气体快速筛选方法 | |
Wang et al. | Electromagnetic-thermal-flow field coupling simulation of 12-kV medium-voltage switchgear | |
Wang et al. | Investigation of SF 6 arc characteristics under shock condition in a supersonic nozzle with hollow contact | |
Zhu et al. | Study of the switching arc characteristics of a 500 kV HVDC self-excited oscillatory metallic return transfer breaker | |
Ye et al. | CFD analysis of the thermal interruption process of gases with low environmental impact in high voltage circuit breakers | |
Choi et al. | Arc gas-flow simulation algorithm considering the effects of nozzle ablation in a self-Blast GCB | |
Bo et al. | Simulation of arc dwelling behavior during bridge-type contacts opening process for high-power relay | |
CN109031110B (zh) | 一种谐波条件下的断路器开断能力评估方法 | |
CN116564676A (zh) | 一种替代sf6的环保绝缘混合气体及其配制方法和应用 | |
Abid et al. | Effect of filling pressure on post-arc gap recovery of N 2 | |
Pu et al. | Study of the Arc Interruption Performance of SF6 Alternative Gases in Load Switch | |
CN115329688A (zh) | 低频电流条件下的大电流真空电弧参数的获取方法及*** | |
Chen et al. | Analysis of arc characteristics of C 4 F 7 N/CO 2 mixed gas ring main unit | |
Wang et al. | Multi physical field simulation of medium voltage switchgear and optimal design | |
Wang et al. | Study on the dielectric recovery performance of high voltage SF6 circuit breaker | |
Liu et al. | Research on Arcing Characteristics of Sulfur Hexafluoride Alternative Gases in Rotating Arc Ring Main Unit | |
Zhong et al. | Simulation and Comparison of Arc Characteristics in C 4 F 7 N, CO 2 and SF 6 | |
Cui et al. | Three-dimensional Simulation of Hot Gas in Cooling Chamber and Tank of High-Voltage Circuit Breaker for alternative gas | |
Park et al. | Evaluation on Short Line Fault Breaking Performance of SF $ _ {6} $ Gas Circuit Breaker Considering Effects of Ablated Nozzle Vapor | |
Zhang et al. | Study of the Switching Arc Characteristics in a HVDC High Speed Switch | |
Nan et al. | Modelling of arcing phenomena during opening contacts in novel circuit breaker | |
Li et al. | Research about influence of contact burning loss on dielectric recovery characteristic for SF 6 circuit breaker | |
Tanaka et al. | Hybrid Thermofluid Modeling with LTE and non-LTE Assumption for Decaying Molecular Gas Arcs | |
Wu et al. | Experimental and theoretical study of decay and post-arc phases of a SF6 transfer arc in DC hybrid breaking |
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 |