CN102298671B - 一种实现电网故障重演的仿真方法 - Google Patents

一种实现电网故障重演的仿真方法 Download PDF

Info

Publication number
CN102298671B
CN102298671B CN201110179311.5A CN201110179311A CN102298671B CN 102298671 B CN102298671 B CN 102298671B CN 201110179311 A CN201110179311 A CN 201110179311A CN 102298671 B CN102298671 B CN 102298671B
Authority
CN
China
Prior art keywords
fault
data
line
model
online
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
CN201110179311.5A
Other languages
English (en)
Other versions
CN102298671A (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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hebei Electric Power Co Ltd
Hebei Electric Power Construction Adjustment Test Institute
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hebei Electric Power Co Ltd
Hebei Electric Power Construction Adjustment Test Institute
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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hebei Electric Power Co Ltd, Hebei Electric Power Construction Adjustment Test Institute filed Critical State Grid Corp of China SGCC
Priority to CN201110179311.5A priority Critical patent/CN102298671B/zh
Priority to PCT/CN2011/083182 priority patent/WO2013000241A1/zh
Publication of CN102298671A publication Critical patent/CN102298671A/zh
Application granted granted Critical
Publication of CN102298671B publication Critical patent/CN102298671B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • 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
    • Y04S40/00Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
    • Y04S40/20Information technology specific aspects, e.g. CAD, simulation, modelling, system security

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种实现电网故障重演的仿真方法,本发明的技术要点是:首先通过在线运行状态潮流数据与电网离线数据的整合,实现了电网故障时的网络拓扑结构及电源点的出力情况,最大程度保证了仿真数据的准确性。再次,通过建立电网故障地点的电磁暂态局部模型,可以准确模拟系统故障的类型、故障持续时间及故障点。最后,通过仿真数据与电网实际录波数据的对比,可以验证利用本发明提出的方法所得到仿真数据的精度。本发明的优点在于:为电网分析人员提供了一种重现电网异常运行方式的方法,可以分析各种保护装置误动因素,找到主要因素,从而提高了电网分析的准确性。

Description

一种实现电网故障重演的仿真方法
技术领域
本发明涉及一种实现电网故障重演的仿真方法,适用于开展电网运行方式分析以及电网继电保护装置动作特性研究技术领域。
背景技术
随着电力系统的发展,电网运行特性日趋复杂。同杆并架双回线路互感、网架结构的变化及电磁环网等因素对继电保护装置的动作正确性产生了较大的影响。电网事故发生后,第一时间分析继电保护装置的动作正确性是技术人员首先需要开展的工作。但是,系统的运行方式是瞬息变化的,而保护误动作的原因是多样的。事故后的故障分析是从多种原因中找到最可能的答案,得到的结论往往是“也许、可能”等不确定性的结论,不利于潜在危险因素的排查及反事故措施的制定推广。
客观需求决定了实现电网故障重演的必要性。实现电网故障重演,可以精确模拟系统故障过程中的各种电气设备电气量的电磁暂态的变化,重现与实际系统故障完全一致的仿真环境。
发明内容
本发明所要解决的问题是提供一种实现电网实际故障重演的仿真方法。
本发明解决其技术问题所采取的技术方案为:一种实现电网故障重演的仿真方法,其特征在于其具体步骤如下:
(1)通过观察电网故障录波图获取电网故障的实际类型及故障发生的时间,通过查勘现场及观察线路保护装置的故障测距记录,确定电网故障的发生地点;
(2)利用电力设备的实际参数及拓扑结构建立电网故障发生地点的局部电磁暂态模型;
(3)利用电网在线运行状态的潮流断面数据方法,根据电网故障发生前最后一个数据断面建立全网机电暂态模型;
(4)利用仿真工具进行全网机电暂态与局部电磁暂态的混合仿真,并在所述局部电磁暂态模型中设置与实际电网故障地点及故障类型一致的故障模型;在全网机电暂态与局部电磁暂态的混合仿真中,所建立的局部电磁暂态模型与全网机电暂态模型的子网应对应一致;
(5)将仿真得到的故障数据与变电站内录波器的实际录波文件进行波形对比,当波形一致时,则验证所建立的所述故障模型是准确的;当差异较大时,则将步骤(4)中建立的所述故障模型的过渡电阻及故障点位置进行微调,最终保证波形一致。
本发明的工作原理如下:
电力系统的运行方式是瞬息变化的。电网分析人员为了捕捉电网的某些特定运行方式(如电网故障状态),必须采取一些特定的办法及特殊的工具。此外,特定运行方式的重现与实际运行方式相比必须保证一定的精度,这样才能作为电网分析人员的有效分析工具。
本发明提供了一种只需较小工作量即可准确实现电网实际故障重演的仿真方法。首先通过在线运行状态潮流数据与电网离线数据的整合,实现了电网故障时的网络拓扑结构及电源点的出力情况,最大程度保证了仿真数据的准确性。再次,通过建立电网故障地点的电磁暂态局部模型,可以准确模拟系统故障的类型、故障持续时间及故障点。最后,通过仿真数据与电网实际录波数据的对比,可以验证利用本发明提出的方法所得到仿真数据的精度。
本发明的有益效果是为电网分析人员提供了一种重现电网异常运行方式的方法,能够分析各种保护装置误动因素,找到主要因素,从而提高了电网分析的准确性。
附图说明:
图1为500kV石北-清苑线路发生C相接地故障的故障录波图。其中,图1-1为石北站故障录波图,图1-2为清苑站故障录波图;
图2为利用故障发生时的电网实时数据断面建立的全网机电暂态模型;
图3为故障发生地点的局部电磁暂态模型;
图4为仿真数据与实际故障数据的对比。其中,图4-1为石北侧C相短路电流对比图,图4-2为清苑侧C相短路电流对比图。
在图3中,1为故障线路一端的机电暂态接口, 2为故障线路另一端的机电暂态接口;在图4中实线为仿真电流波形图,虚线为实际电流波形图。
具体实施方式
下面结合附图和实施例对本发明再做进一步详细的说明。
实施例1:
本实施例是对500KV石北-清苑线路发生C相接地故障的重演仿真方法。
本实施例的方法步骤如下:
(1)通过观察电网故障录波图获取电网故障的实际类型及故障发生的时间,通过查勘现场及观察线路保护装置的故障测距记录,确定电网故障的发生地点;
参见图1,图1为500kV石北-清苑线路发生C相接地故障的故障录波图。通过分析,可从录波图中获知故障发生的时间及C相故障以及电气量的变化(波形)。
(2)利用电力设备(如发电机、变压器、负荷、线路)的实际参数及拓扑结构建立电网故障发生地点的局部电磁暂态模型;
参见图3,图3为手动建立的电网故障点局部电磁暂态模型。图中,BUS1、BUS2分别对应了图2中同名的母线,BUS3、BUS4、BUS6、BUS7是为了添加断路器元件而增设的母线,BUS5是为了添加故障元件Fault1而增设的母线。Breaker1-Breaker4为输电线路Line1、Line2两侧的断路器,在仿真模型中,可以通过断路器元件采集流经故障线路的短路电流。Line1、Line2分别对应了图2中同名的输电线路,每条输电线路均采用π型等效电路。为了增设故障元件Fault1,将Line1分割为两条线路Line1-1及Line1-2,但Line1-1及Line1-2的线路长度之和与Line1相等。Fault1为故障元件,可以通过故障元件设置故障类型(单相接地故障、相间故障、相间接地故障等)。同时,也可以设置过渡电阻、故障发生时间及持续时间。Fault1所连接的母线BUS5左右两侧线路Line1-1及Line1-2的长度代表了电网故障发生的地点。
故障发生在线路Line1侧,通过采集Breaker1、Breaker3的故障相电流并与实际录波数据进行对比,可以验证故障重演仿真模型的仿真精度。
(3)利用公开专利201010553716.6所实现的电网在线运行状态的潮流断面数据方法,根据电网故障发生前最后一个数据断面建立全网机电暂态模型;
所述建立全网机电暂态模型的具体步骤如下:
a、自动或手动从EMS 调度自动化系统中获取状态估计后的省级电网模型、在线运行数据断面,导入到电网数字仿真系统的基础数据库中;
b、基于在线电网模型和运行数据进行拓扑分析和潮流计算,验证在线潮流合理性;
c、进行在线和离线数据整合;
d、自动生成或者按用户指定的间隔生成在线运行方式数据;
所述在线和离线数据整合是通过维护和设置,定时合并离线/ 在线电网数据,形成一套反映当前/ 历史电网工况并能进行暂态稳定分析的计算数据,为在线电网数字仿真提供数据源;
所述在线和离线数据整合的方法如下:
①以离线数据为基础,用在线运行数据进行修正;使整合后的数据,保留离线数据原有的网络结构、元件名称;
② 建立在线/ 离线元件名称映射表,通过所述映射表实现在线/ 离线数据的共用;
参见图2,图2为利用电网故障时刻前的在线潮流数据断面建立的全网机电暂态模型,用来模拟系统的实际运行背景。在机电暂态模型上划分两个子网,将需要建立电磁暂态模型的部分划为一个子网EMT,包括BUS1、BUS2及两条母线的联络线line1、line2。将电网其它设备划为另一个子网ST。图中,BUS1、BUS2分别代表了石北站及清苑站的500kV母线。
(4)利用仿真工具进行全网机电暂态与局部电磁暂态的混合仿真,并在所述局部电磁暂态模型中设置与实际电网故障地点及故障类型一致的故障模型;在全网机电暂态与局部电磁暂态的混合仿真中,所建立的局部电磁暂态模型与全网机电暂态模型的子网应对应一致;
(5)将仿真得到的故障数据与变电站内录波器的实际录波文件进行波形对比,当波形一致时,则验证所建立的所述故障模型是准确的;当差异较大时,则将步骤(4)中建立的所述故障模型的过渡电阻及故障点位置进行微调,最终保证波形一致。
参见图4,图4为仿真电流与电网实际故障电流的对比图。图4-1为图2及图3中的BUS1侧的短路电流波形对比图,图4-2为图2及图3中BUS2侧的短路电流波形对比图。其中,虚线标识的是故障录波图记录的故障电流波形,实线标识的是本实施例获取的仿真电流波形。从图中可以看出,故障电流的拟和程度非常高。这样可以验证,利用本发明提出的方法可以实现精确的电网故障重演。

Claims (1)

1.一种实现电网故障重演的仿真方法,其特征在于其具体步骤如下:
(1)通过观察电网故障录波图获取电网故障的实际类型及故障发生的时间,通过查勘现场及观察线路保护装置的故障测距记录,确定电网故障的发生地点;
(2)利用电力设备的实际参数及拓扑结构建立电网故障发生地点的局部电磁暂态模型;
(3)利用电网在线运行状态的潮流断面数据方法,根据电网故障发生前最后一个数据断面建立全网机电暂态模型;
所述建立全网机电暂态模型的具体步骤如下:
a、 自动或手动从EMS 调度自动化系统中获取状态估计后的省级电网模型、在线运行数据断面,导入到电网数字仿真系统的基础数据库中;
b、基于在线电网模型和运行数据进行拓扑分析和潮流计算,验证在线潮流合理性;
c、进行在线和离线数据整合;
d、 自动生成或者按用户指定的间隔生成在线运行方式数据;
所述在线和离线数据整合是通过维护和设置,定时合并离线/ 在线电网数据,形成一套反映当前/ 历史电网工况并能进行暂态稳定分析的计算数据,为在线电网数字仿真提供数据源;
所述在线和离线数据整合的方法如下:
①以离线数据为基础,用在线运行数据进行修正;使整合后的数据,保留离线数据原有的网络结构、元件名称;
② 建立在线/ 离线元件名称映射表,通过所述映射表实现在线/ 离线数据的共用;
(4)利用仿真工具进行全网机电暂态与局部电磁暂态的混合仿真,并在所述局部电磁暂态模型中设置与实际电网故障地点及故障类型一致的故障模型;在全网机电暂态与局部电磁暂态的混合仿真中,所建立的局部电磁暂态模型与全网机电暂态模型的子网应对应一致;
(5)将仿真得到的故障数据与变电站内录波器的实际录波文件进行波形对比,当波形一致时,则验证所建立的所述故障模型是准确的;当差异较大时,则将步骤(4)中建立的所述故障模型的过渡电阻及故障点位置进行微调,最终保证波形一致。
CN201110179311.5A 2011-06-29 2011-06-29 一种实现电网故障重演的仿真方法 Active CN102298671B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201110179311.5A CN102298671B (zh) 2011-06-29 2011-06-29 一种实现电网故障重演的仿真方法
PCT/CN2011/083182 WO2013000241A1 (zh) 2011-06-29 2011-11-29 一种实现电网故障重演的仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110179311.5A CN102298671B (zh) 2011-06-29 2011-06-29 一种实现电网故障重演的仿真方法

Publications (2)

Publication Number Publication Date
CN102298671A CN102298671A (zh) 2011-12-28
CN102298671B true CN102298671B (zh) 2015-05-13

Family

ID=45359083

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110179311.5A Active CN102298671B (zh) 2011-06-29 2011-06-29 一种实现电网故障重演的仿真方法

Country Status (2)

Country Link
CN (1) CN102298671B (zh)
WO (1) WO2013000241A1 (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103094930B (zh) * 2013-01-18 2014-12-03 南方电网科学研究院有限责任公司 应用在交直流大电网暂态分析中模拟直流响应的方法
CN103825273B (zh) * 2014-03-21 2016-05-25 北京四方继保自动化股份有限公司 一种以馈线联络组为单元的配电网事故全过程反演方法
CN104036430A (zh) * 2014-05-08 2014-09-10 中国石油天然气股份有限公司 炼化企业的电网运行分析方法及系统
CN104049221B (zh) * 2014-07-08 2016-08-24 哈尔滨工业大学 基于滑动窗口和统计信息的电源电压故障诊断方法
CN104269846B (zh) * 2014-09-26 2016-08-17 华北电力大学 电网故障分区方法
CN104698343B (zh) * 2015-03-26 2016-06-08 广东电网有限责任公司电力调度控制中心 基于历史录波数据的电网故障判断方法和系统
CN105445585B (zh) * 2015-11-27 2018-04-13 广州供电局有限公司 电网一次回路的故障诊断方法和系统
CN109116188A (zh) * 2018-09-05 2019-01-01 广东电网有限责任公司 一种电力系统故障推演方法、装置、设备及可读存储介质
CN110348066A (zh) * 2019-06-17 2019-10-18 国网西藏电力有限公司电力科学研究院 用于新能源大规模接入弱电网的继电保护装置校核方法
CN111189440B (zh) * 2019-12-31 2021-09-07 中国电建集团华东勘测设计研究院有限公司 一种基于空间信息模型与实时图像比对的定位导航方法
CN111157940B (zh) * 2020-01-09 2023-05-12 南方电网电力科技股份有限公司 一种柔性直流电子式互感器仿真模拟系统
CN111611663B (zh) * 2020-04-14 2023-05-05 深圳供电局有限公司 一种电力系统仿真故障生成方法及装置
CN112363409B (zh) * 2020-11-10 2022-02-11 中国核动力研究设计院 一种核电厂安全级仪控仿真系统的工况回溯与重演系统
CN113111512B (zh) * 2021-04-13 2022-03-04 中国电力科学研究院有限公司 一种用于大规模电网全电磁暂态仿真模型稳态建立的方法及系统
CN113391156B (zh) * 2021-06-23 2022-11-01 四川华能宝兴河水电有限责任公司 一种电力系统故障再现方法
CN113887005B (zh) * 2021-09-30 2024-02-06 广东电网有限责任公司广州供电局 一种交直流电力系统仿真建模方法、装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101404411A (zh) * 2008-10-30 2009-04-08 中国电力科学研究院 一种负荷建模中的故障拟合方法
CN101846716B (zh) * 2009-03-27 2011-12-28 北海银河科技继保电气有限公司 一种快速捕捉电网故障启动点的方法
CN101539963B (zh) * 2009-04-30 2011-02-02 华北电力大学 一种机电暂态向电磁暂态的模型转换方案及实现方法
CN101640418B (zh) * 2009-09-03 2011-08-24 国网电力科学研究院 电力系统故障的暂态稳定关键输电断面识别方法
CN101957872A (zh) * 2010-07-09 2011-01-26 南方电网技术研究中心 一种交直流电力系统的混合实时仿真方法
CN101957398B (zh) * 2010-09-16 2012-11-28 河北省电力研究院 一种基于机电与电磁暂态混合仿真技术检测计算电网一次时间常数的方法
CN102025151A (zh) * 2010-11-22 2011-04-20 河北省电力研究院 获取电网在线运行状态的潮流断面数据的方法

Also Published As

Publication number Publication date
WO2013000241A1 (zh) 2013-01-03
CN102298671A (zh) 2011-12-28

Similar Documents

Publication Publication Date Title
CN102298671B (zh) 一种实现电网故障重演的仿真方法
CN101968525B (zh) 一种仿真计算与实时监测相结合的配电网故障定位方法
CN103105550B (zh) 一种智能化继电保护装置检测方法及系统
CN104007978B (zh) 一种智能变电站二次设备逻辑仿真系统
Choi et al. Feasibility study: Autonomous state estimation in distribution systems
CN103529418A (zh) 电力系统同步相量测量装置动态测试方法及装置
CN109557423A (zh) 一种配电网故障诊断系统、方法及应用
CN102981079B (zh) 直挂式动态无功发生装置响应波形检测方法
CN105116246A (zh) 一种电网运行保护实时在线分析方法
CN112510654B (zh) 一种母线故障诊断方法、系统、可读介质及设备
CN104571082A (zh) 用于馈线自动化测试的数据再现系统及方法
CN105576660A (zh) 地区配电网在线仿真分析装置及系统
Meliopoulos et al. Delivering accurate and timely data to all
Kincic et al. Bridging the gap between operation and planning models in WECC
Ren et al. Modeling and simulation tools for teaching protective relaying design and application for the smart grid
Andreev et al. Setting of relay protection of electric power systems using its mathematical models
CN203037759U (zh) 动态无功补偿装置响应波形检测装置
Almas et al. Methodologies for power protection relay testing: From conventional to real-time hardware-in-the-loop (HIL) simulation approaches
Saad et al. HIL simulation to assess interaction risks of HVDC systems for upcoming grid development
de Melo et al. Power Quality Monitoring using Synchronized Phasor Measurements: An approach based on hardware-in-the-loop simulations
Kisser et al. A Test-bed System for Validation of Ancillary Services of Wind Power Plants under Realistic Conditions
Zhang et al. Practical issues of implementation of linear state estimator in WECC
CN105119367A (zh) 一种基于动态拓扑的数字化备自投仿真方法
Dufour et al. Testing 750 node distribution grids and devices
e Sousa et al. Viability assessment for centralised protection and control system architectures in medium voltage (MV) substations

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: HEBEI ELECTRIC POWER CORPORATION ELECTRIC POWER RE

Free format text: FORMER OWNER: HEBEI ELECTRIC POWER RESEARCH INSTITUTE

Effective date: 20121220

Owner name: STATE GRID CORPORATION OF CHINA HEBEI ELECTRIC POW

Effective date: 20121220

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20121220

Address after: 050021 No. 238 South Sports street, Yuhua Road, Hebei, Shijiazhuang

Applicant after: Electric Power Research Institute of Hebei Electric Power Corporation

Applicant after: State Grid Corporation of China

Applicant after: Hebei Electric Power Construction & Adjustment Research Institute

Address before: 050021 No. 238 South Sports street, Hebei, Shijiazhuang

Applicant before: Hebei Electric Power Research Institute

C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Fan Hui

Inventor after: Tang Baofeng

Inventor after: Wang Xiaowei

Inventor after: Yang Xiao

Inventor after: Wang Zhilu

Inventor before: Tang Baofeng

Inventor before: Fan Hui

Inventor before: Wang Xiaowei

Inventor before: Yang Xiao

Inventor before: Wang Zhilu

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: TANG BAOFENG FAN HUI WANG XIAOWEI YANG XIAO WANG ZHILU TO: FAN HUI TANG BAOFENG WANG XIAOWEI YANG XIAO WANG ZHILU

C14 Grant of patent or utility model
GR01 Patent grant