CN112198465B - 一种变压器的剩余磁通的检测方法、介质及系统 - Google Patents

一种变压器的剩余磁通的检测方法、介质及系统 Download PDF

Info

Publication number
CN112198465B
CN112198465B CN202010791070.9A CN202010791070A CN112198465B CN 112198465 B CN112198465 B CN 112198465B CN 202010791070 A CN202010791070 A CN 202010791070A CN 112198465 B CN112198465 B CN 112198465B
Authority
CN
China
Prior art keywords
transformer
magnetic
residual
maximum
dynamic
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
CN202010791070.9A
Other languages
English (en)
Other versions
CN112198465A (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.)
Anhui Huidian Technology Co ltd
Wuhan NARI Ltd
Xian Jiaotong University
State Grid Ningxia Electric Power Co Ltd
Electric Power Research Institute of State Grid Ningxia Electric Power Co Ltd
Original Assignee
Anhui Huidian Technology Co ltd
Wuhan NARI Ltd
Xian Jiaotong University
State Grid Ningxia Electric Power Co Ltd
Electric Power Research Institute of State Grid Ningxia Electric Power 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 Anhui Huidian Technology Co ltd, Wuhan NARI Ltd, Xian Jiaotong University, State Grid Ningxia Electric Power Co Ltd, Electric Power Research Institute of State Grid Ningxia Electric Power Co Ltd filed Critical Anhui Huidian Technology Co ltd
Priority to CN202010791070.9A priority Critical patent/CN112198465B/zh
Publication of CN112198465A publication Critical patent/CN112198465A/zh
Application granted granted Critical
Publication of CN112198465B publication Critical patent/CN112198465B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/02Measuring direction or magnitude of magnetic fields or magnetic flux

Abstract

本发明公开一种变压器的剩余磁通的检测方法、介质及系统。该检测方法包括:获取变压器铁芯的饱和磁滞回线的饱和点处的最大磁感应强度以及最大磁感应强度对应的阈值磁场强度;计算变压器铁芯无剩磁时的动态磁阻最大值;通过变压器绕组端子外接的可调交流电压源施加初始相角为0°的一个周期正弦电压激励,并测量一个周期内的电压值和电流值;根据测量得到的电压值和电流值,以及,饱和磁滞回线中最大磁感应强度对应的阈值磁场强度,计算得到剩磁条件下的变压器铁芯的动态磁阻;根据变压器铁芯的动态磁阻、无剩磁时的动态磁阻最大值和最大磁感应强度,计算得到变压器的剩余磁通。本发明准确检测变压器铁芯的剩余磁通。

Description

一种变压器的剩余磁通的检测方法、介质及系统
技术领域
本发明涉及变压器技术领域,尤其涉及一种变压器的剩余磁通的检测方法、介质及系统。
背景技术
变压器空载投运时,由于铁芯中可能存在的剩余磁通与合闸预期磁通不相等,会导致励磁涌流出现,其幅值可能达到6-8倍额定电流,对电网运行产生危害,引起继电保护装置误动作。利用选相合闸技术能够有效抑制励磁涌流的产生。目前常用的选相合闸技术,不考虑剩余磁通而直接在电压幅值处选相合闸,对励磁涌流有一定的抑制效果。但是,由于合闸预击穿以及剩余磁通的影响,其效果往往不够理想。故选相合闸技术需要对铁芯剩余磁通进行有效检测。由于铁芯磁路封闭,无法对其进行直接测量,所以目前已有方法均为间接测量。一般的间接测量方法,往往是对局部磁滞回线进行测量后,利用外加激励分析电压、电流变化率等在磁滞回线中定位剩余磁通大小,但是,激励过后剩余磁通大小本身也会发生变化,即测量所得到的剩余磁通与最终受影响的剩余磁通不相等。并且现有方案对局部磁滞回线的分析往往直接使用饱和磁滞回线来代替,然而铁芯局部磁滞回线与饱和磁滞回线的B-H关系并不相同,所以计算结果存在误差。
发明内容
本发明实施例提供了一种变压器的剩余磁通的检测方法、介质及系统,以解决现有技术检测得到的剩余磁通存在误差的问题。
第一方面,提供了一种变压器的剩余磁通的检测方法,包括:获取变压器铁芯的饱和磁滞回线的饱和点处的最大磁感应强度以及最大磁感应强度对应的阈值磁场强度;计算变压器铁芯无剩磁时的动态磁阻最大值;通过变压器绕组端子外接的可调交流电压源施加初始相角为0°的一个周期正弦电压激励,并测量一个周期内的电压值和电流值;根据测量得到的电压值和电流值,以及,饱和磁滞回线中最大磁感应强度对应的阈值磁场强度,计算得到剩磁条件下的变压器铁芯的动态磁阻;根据变压器铁芯的动态磁阻、无剩磁时的动态磁阻最大值和最大磁感应强度,计算得到变压器的剩余磁通。
第二方面,提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序指令;所述计算机程序指令被处理器执行时实现如第一方面实施例所述的变压器的剩余磁通的检测方法。
第三方面,提供了一种变压器的剩余磁通的检测系统,包括:如第二方面实施例所述的计算机可读存储介质。
本发明实施例,可准确检测变压器铁芯的剩余磁通,并不引起变压器的剩余磁通变化,特别适用于采用选相分合闸技术的空载变压器投切领域。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对本发明实施例的描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例的变压器的剩余磁通的检测方法的流程图;
图2是本发明实施例的可调交流电压源激励一个周期局部磁滞回线变化图与饱和磁滞回线的对比示意图;
图3是本发明实施例的采用可调交流源的剩余磁通检测装置的原理图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获取的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种变压器的剩余磁通的检测方法。该检测方法针对选相投切空载变压器投切时,剩余磁通与合闸预期磁通不相符时导致的励磁涌流现象。该方法利用可调交流电源进行激励,利用变压器动态磁阻与磁通的对应关系,对变压器断电后进行一个周期正弦电压激励,测量其电流、电压,并计算出对应的铁芯磁场强度与磁感应强度关系,准确检测变压器铁芯剩余磁通数值,并不引起变压器剩余磁通变化,特别适用于采用选相分合闸技术的空载变压器投切领域。
如图1所示,该检测方法具体包括如下的步骤:
步骤S1:获取变压器铁芯的饱和磁滞回线的饱和点处的最大磁感应强度以及最大磁感应强度对应的阈值磁场强度。
一般的,变压器铁芯在出厂时,会标有最大磁感应强度以及最大磁感应强度对应的阈值磁场强度等参数,因此,可直接读取这些数据。
步骤S2:计算变压器铁芯无剩磁时的动态磁阻最大值。
具体的,该步骤包括如下的过程:
(1)在变压器铁芯的剩余磁通为零时,通过变压器绕组端子外接的可调交流电压源施加正弦电压激励。
一般的,刚出厂的变压器铁芯是没有剩余磁通的,可在这个时候施加电压激励。如图3所示,通过可调交流电压源101对铁芯绕组进行交流电压源激励。应当理解的是,该电压激励对初始相角没有要求。
(2)测量得到初始时刻的电压值和电流值。
如图3所示,利用电压互感器102、电流互感器103、采集卡104获得绕组端电压与绕组中的电流。
(3)采用
Figure BDA0002623760420000041
计算得到动态磁阻最大值。
其中,Rm表示无剩磁时的动态磁阻最大值,Ut0表示初始时刻的电压值,it0表示初始时刻的电流值,H0表示饱和磁滞回线中最大磁感应强度Bs对应的阈值磁场强度,N表示变压器励磁线圈的匝数,S表示铁芯有效截面积,l表示变压器铁芯的有效磁路长度。H0,N,S,l都可直接通过变压器铁芯出厂时标注的参数直接得到。
该计算公式由下述几个公式推导得到:
Figure BDA0002623760420000042
Figure BDA0002623760420000043
Figure BDA0002623760420000044
其中,公式(2)可根据磁通变化量
Figure BDA0002623760420000045
得到。
通过上述三个公式可以联立推导得到
Figure BDA0002623760420000046
该公式中,r(B)表示变压器铁芯的动态磁阻。U表示电压值,i表示电流值,H0表示饱和磁滞回线中最大磁感应强度Bs对应的阈值磁场强度,N表示变压器励磁线圈的匝数,S表示铁芯有效截面积,l表示变压器铁芯的有效磁路长度。
此外,r(B)和Rm的关系可通过下式表示:
Figure BDA0002623760420000047
一般的,r(B)取值与铁芯材料和结构有关,对于不同的变压器铁芯,需要专门进行测量。通过上式可知,当刚开始给绕组通电时,磁感应强度B为0,此时变压器铁芯动态磁阻r(B)最大,即此时r(B)=Rm。因此,动态磁阻最大值可采用电压激励的初始时刻的电压值和电流值通过上述的公式计算得到。如图3所示,利用数字信号处理模块105进行相应的计算以及后续的数据计算。
步骤S3:通过变压器绕组端子外接的可调交流电压源施加初始相角为0°的一个周期正弦电压激励,并测量一个周期内的电压值和电流值。
该电压激励和测量也可通过前述的图3所示的装置进行施加和测量。
如图2所示,在变压器断电后,对于含有剩余磁通的铁芯进行一个周期的交流源激励,电压从零开始增大,忽略铜损铁损。由于电压等级远低于变压器额定电压,所以将电压源激励下的局部磁滞回线近似为其磁化曲线,即一个周期之内,磁滞回线变化路径为从R点出发沿下曲线(即上升沿曲线)到达P点,再沿上曲线(即下降沿曲线)回到R点,等效为R点出发的直线到达P点再沿原路返回R点。磁感应强度与磁场强度在交流源激励正半周期内增加,在负半周期内减小,最终回到R点。
此外,由于磁通由电压积分可得,所以由
Figure BDA0002623760420000051
可知,对剩磁进行从0°开始的一个周期正弦电压激励测量并不会改变初始剩余磁通。
综上,本步骤施加初始相角为0°的一个周期正弦电压激励,在一次测量周期内,可以不改变剩余磁通大小,完成对剩余磁通的测量,在后续步骤中通过采集的电压值和电流值进行相应的计算最终得到剩余磁通。
步骤S4:根据测量得到的电压值和电流值,以及,饱和磁滞回线中最大磁感应强度对应的阈值磁场强度,计算得到剩磁条件下的变压器铁芯的动态磁阻。
具体的,计算得到剩磁条件下的变压器铁芯动态磁阻的计算式为(该公式的推导原理如步骤S2中所示,在此不再赘述):
Figure BDA0002623760420000061
其中,r(B)表示变压器铁芯的动态磁阻,U(t)表示t时刻测量得到的电压值,i(t)表示t时刻测量得到的电流值,H0表示饱和磁滞回线中最大磁感应强度Bs对应的阈值磁场强度,N表示变压器励磁线圈的匝数,S表示变压器铁芯有效截面积,l表示变压器铁芯的有效磁路长度。应当理解的是,测量得到的电压值和电流值可以为正弦电压激励半个周期时对应的电压值和电流值。
步骤S5:根据变压器铁芯的动态磁阻、无剩磁时的动态磁阻最大值和最大磁感应强度,计算得到变压器的剩余磁通。
具体的,计算得到变压器的剩余磁通的计算式为:
Figure BDA0002623760420000062
其中,BR表示变压器的剩余磁通,BS表示最大磁感应强度,r(B)表示变压器铁芯的动态磁阻,Rm表示无剩磁时的动态磁阻最大值。
如前所述,r(B)和Rm的关系可通过下式表示:
Figure BDA0002623760420000063
该式中的B表示磁感应强度。对于动态磁阻r(B)可以是磁感应强度B的函数。此外,r(B)也可以是常数,这是因为本发明实施例的激励所采用的电压等级远低于变压器额定电压,所以上式中的B取值为变压器的剩余磁通BR,则r(B)为常数。因此,
Figure BDA0002623760420000071
可记为
Figure BDA0002623760420000072
从而可利用动态磁阻与磁感应强度的关系计算变压器的剩余磁通。如前所述,由于电压源激励为一个完整周期,所以磁通变化量为0,测量结束剩余磁通仍然维持初始值。
对10kV冲击变压器进行测量,并对其进行了合闸涌流测试,利用合闸涌流推导计算得出的剩磁值作为参考,对本发明实施例的通过电压激励并计算得到剩余磁通的方法进行了验证,实验结果如表1所示。
表1合闸涌流验证结果
合闸涌流(A) 涌流计算剩磁(p.u.) 测量剩磁(p.u.)
620 -0.73 -0.68
520 -0.64 -0.63
300 -0.39 -0.42
-140 0.69 0.67
-260 -0.20 -0.23
-430 0.28 0.30
198 0.73 0.74
-280 0.68 0.66
340 -0.59 -0.56
460 0.38 0.40
320 -0.52 -0.53
500 -0.69 -0.66
通过表1的结果,可以看出本发明实施例的方法具有较高的准确性,因此,具有实用性。
本发明实施例还公开了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序指令;所述计算机程序指令被处理器执行时实现如上述实施例所述的变压器的剩余磁通的检测方法。
本发明实施例还公开了一种变压器的剩余磁通的检测系统,包括:如上述实施例所述的计算机可读存储介质。
综上,本发明实施例,可准确检测变压器铁芯的剩余磁通,并不引起变压器的剩余磁通变化,特别适用于采用选相分合闸技术的空载变压器投切领域。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (5)

1.一种变压器的剩余磁通的检测方法,其特征在于,包括:
获取变压器铁芯的饱和磁滞回线的饱和点处的最大磁感应强度以及最大磁感应强度对应的阈值磁场强度;
计算变压器铁芯无剩磁时的动态磁阻最大值;
通过变压器绕组端子外接的可调交流电压源施加初始相角为0°的一个周期正弦电压激励,并测量一个周期内的电压值和电流值;
根据测量得到的电压值和电流值,以及,饱和磁滞回线中最大磁感应强度对应的阈值磁场强度,计算得到剩磁条件下的变压器铁芯的动态磁阻;
根据变压器铁芯的动态磁阻、无剩磁时的动态磁阻最大值和最大磁感应强度,计算得到变压器的剩余磁通,所述计算得到变压器的剩余磁通的计算式为:
Figure FDA0003656196570000011
其中,BR表示变压器的剩余磁通,BS表示最大磁感应强度,r(B)表示变压器铁芯的动态磁阻,Rm表示无剩磁时的动态磁阻最大值,所述计算得到剩磁条件下的变压器铁芯的动态磁阻的计算式为:
Figure FDA0003656196570000012
其中,r(B)表示变压器铁芯的动态磁阻,U(t)表示t时刻测量得到的电压值,i(t)表示t时刻测量得到的电流值,H0表示饱和磁滞回线中最大磁感应强度Bs对应的阈值磁场强度,N表示变压器励磁线圈的匝数,S表示变压器铁芯有效截面积,l表示变压器铁芯的有效磁路长度。
2.根据权利要求1所述的变压器的剩余磁通的检测方法,其特征在于,所述计算变压器铁芯无剩磁时的动态磁阻最大值的步骤,包括:
在变压器铁芯的剩余磁通为零时,通过变压器绕组端子外接的可调交流电压源施加正弦电压激励;
测量得到初始时刻的电压值和电流值;
采用
Figure FDA0003656196570000021
计算得到动态磁阻最大值;
其中,Rm表示无剩磁时的动态磁阻最大值,Ut0表示初始时刻的电压值,it0表示初始时刻的电流值,H0表示饱和磁滞回线中最大磁感应强度Bs对应的阈值磁场强度,N表示变压器励磁线圈的匝数,S表示铁芯有效截面积,l表示变压器铁芯的有效磁路长度。
3.根据权利要求1所述的变压器的剩余磁通的检测方法,其特征在于:测量得到的电压值和电流值为正弦电压激励半个周期时对应的电压值和电流值。
4.一种计算机可读存储介质,其特征在于:所述计算机可读存储介质上存储有计算机程序指令;所述计算机程序指令被处理器执行时实现如权利要求1~3中任一项所述的变压器的剩余磁通的检测方法。
5.一种变压器的剩余磁通的检测系统,其特征在于,包括:如权利要求4所述的计算机可读存储介质。
CN202010791070.9A 2020-08-07 2020-08-07 一种变压器的剩余磁通的检测方法、介质及系统 Active CN112198465B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010791070.9A CN112198465B (zh) 2020-08-07 2020-08-07 一种变压器的剩余磁通的检测方法、介质及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010791070.9A CN112198465B (zh) 2020-08-07 2020-08-07 一种变压器的剩余磁通的检测方法、介质及系统

Publications (2)

Publication Number Publication Date
CN112198465A CN112198465A (zh) 2021-01-08
CN112198465B true CN112198465B (zh) 2022-08-09

Family

ID=74005530

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010791070.9A Active CN112198465B (zh) 2020-08-07 2020-08-07 一种变压器的剩余磁通的检测方法、介质及系统

Country Status (1)

Country Link
CN (1) CN112198465B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05322575A (ja) * 1992-05-18 1993-12-07 Casio Comput Co Ltd 電子式方位計
JP2009055750A (ja) * 2007-08-29 2009-03-12 Nidec Servo Corp Pmクローポール型モータとその製作方法
CN102403083A (zh) * 2010-09-07 2012-04-04 中国石油天然气集团公司 一种管道端口磁滞缓冲消磁方法及消磁器
CN102565727A (zh) * 2012-02-20 2012-07-11 江苏多维科技有限公司 用于测量磁场的磁电阻传感器
EP2843433A1 (en) * 2012-11-16 2015-03-04 Korea Basic Science Institute Spin torque-type magnetic sensor
CN107677975A (zh) * 2017-10-10 2018-02-09 广东电网有限责任公司电力科学研究院 一种变压器饱和状态下的励磁曲线的测量方法
CN109444776A (zh) * 2018-11-05 2019-03-08 国网江苏省电力有限公司电力科学研究院 三相三柱变压器铁芯剩磁测算方法、系统及存储介质
EP3597899A1 (en) * 2013-07-29 2020-01-22 Hitachi Automotive Systems, Ltd. Drive device for fuel injection device, and fuel injection system

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FI20030333A (fi) * 2002-03-15 2003-09-16 Adaptamat Tech Oy MSM-anturi
US20060225973A1 (en) * 2005-03-30 2006-10-12 Dimig Steven J Residual magnetic devices and methods

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05322575A (ja) * 1992-05-18 1993-12-07 Casio Comput Co Ltd 電子式方位計
JP2009055750A (ja) * 2007-08-29 2009-03-12 Nidec Servo Corp Pmクローポール型モータとその製作方法
CN102403083A (zh) * 2010-09-07 2012-04-04 中国石油天然气集团公司 一种管道端口磁滞缓冲消磁方法及消磁器
CN102565727A (zh) * 2012-02-20 2012-07-11 江苏多维科技有限公司 用于测量磁场的磁电阻传感器
EP2843433A1 (en) * 2012-11-16 2015-03-04 Korea Basic Science Institute Spin torque-type magnetic sensor
EP3597899A1 (en) * 2013-07-29 2020-01-22 Hitachi Automotive Systems, Ltd. Drive device for fuel injection device, and fuel injection system
CN107677975A (zh) * 2017-10-10 2018-02-09 广东电网有限责任公司电力科学研究院 一种变压器饱和状态下的励磁曲线的测量方法
CN109444776A (zh) * 2018-11-05 2019-03-08 国网江苏省电力有限公司电力科学研究院 三相三柱变压器铁芯剩磁测算方法、系统及存储介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Ke Wang 等.Research on Residual Flux Characteristics of Transformer with Single-Phase Four-Limb Core Under Different DC Excitation Current.《2020 IEEE International Conference on High Voltage Engineering and Application (ICHVE)》.2020, *
杨向东 等.螺旋埋弧焊管剩磁机理及剩磁综合控制措施.《焊管》.2015,第38卷(第03期), *
荣毅 等.真空断路器截流对变压器铁心剩磁的影响.《宁夏电力》.2019,(第03期), *

Also Published As

Publication number Publication date
CN112198465A (zh) 2021-01-08

Similar Documents

Publication Publication Date Title
CN100477437C (zh) 剩余电流设备
Locci et al. A digital compensation method for improving current transformer accuracy
Locci et al. Hysteresis and eddy currents compensation in current transformers
CN105425175A (zh) 一种基于极性反转的直流电压源的铁芯剩磁消除和测量方法
De Leon et al. Accurate measurement of the air-core inductance of iron-core transformers with a non-ideal low-power rectifier
Shepard et al. An overview of Rogowski coil current sensing technology
Emanuel et al. Current harmonics measurement by means of current transformers
Bhowmick et al. Online detection of an interturn winding fault in single-phase distribution transformers using a terminal measurement-based modeling technique
CN103913710A (zh) 一种变压器铁芯剩余磁感应强度的测量方法
Vukosavić et al. High-precision sensing of DC bias in AC grids
CN105182261B (zh) 线圈内部铁芯磁场强度非接触式测量方法
Wath et al. Error compensation method for current transformer
Vulin et al. Determining the residual magnetic flux value of a single-phase transformer using a minor hysteresis loop
Huo et al. Residual flux measurement of the single-phase transformer based on transient current method
US8934210B1 (en) Demagnetization using a determined estimated magnetic state
CN111965569B (zh) 在线运行变压器的直流偏磁水平测量方法、装置及系统
Bessolitsyn et al. Experimental study of current error of up to 50 hz current-measuring transformer
CN112198465B (zh) 一种变压器的剩余磁通的检测方法、介质及系统
CN104614688B (zh) 用于直流偏磁动态磁通测量的c型传感器及其检测方法
CN112345966A (zh) 一种通过剩余电流互感器检测直流漏电的方法及其装置
Cataliotti et al. Characterization of clamp-on current transformers under nonsinusoidal conditions
Petrun et al. Evaluation of iron core quality for resistance spot welding transformers using current controlled supply
Wang et al. Study of current transformer calibrating system based on equivalent model
US4827214A (en) Method for measuring losses of a magnetic core
Rajendran et al. Saturation analysis on current transformer

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