CN108647407B - 一种转炉炼钢烟气分析定碳方法 - Google Patents

一种转炉炼钢烟气分析定碳方法 Download PDF

Info

Publication number
CN108647407B
CN108647407B CN201810372680.8A CN201810372680A CN108647407B CN 108647407 B CN108647407 B CN 108647407B CN 201810372680 A CN201810372680 A CN 201810372680A CN 108647407 B CN108647407 B CN 108647407B
Authority
CN
China
Prior art keywords
converter
flue gas
carbon
smelting
molten
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
CN201810372680.8A
Other languages
English (en)
Other versions
CN108647407A (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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201810372680.8A priority Critical patent/CN108647407B/zh
Publication of CN108647407A publication Critical patent/CN108647407A/zh
Application granted granted Critical
Publication of CN108647407B publication Critical patent/CN108647407B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21CPROCESSING OF PIG-IRON, e.g. REFINING, MANUFACTURE OF WROUGHT-IRON OR STEEL; TREATMENT IN MOLTEN STATE OF FERROUS ALLOYS
    • C21C5/00Manufacture of carbon-steel, e.g. plain mild steel, medium carbon steel or cast steel or stainless steel
    • C21C5/28Manufacture of steel in the converter

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Manufacturing & Machinery (AREA)
  • Computer Hardware Design (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Carbon Steel Or Casting Steel Manufacturing (AREA)

Abstract

本发明提供一种转炉炼钢烟气分析定碳方法,属于钢铁冶金领域。本方法包括步骤为:首先,依据质量守恒原理,结合转炉冶炼过程中检测设备获取的烟气CO、CO2百分含量和烟气流量计算出熔池脱碳速率;其次,将熔池混匀度的概念引入到终点碳曲线拟合指数模型;然后,根据冶金熔渣分子理论结合生产钢种特性,获取冶炼终渣成分和冶炼终点钢水温度,计算出冶炼熔池极限碳含量W[C]0;最终,得到转炉炼钢烟气分析定碳模型。本方法中的模型引入熔池混匀度概念,充分考虑了氧枪枪位、顶吹氧气流量以及底吹气体流量操作工艺参数对转炉熔池脱碳过程的影响;并根据实际冶炼钢种利用熔渣活度计算得出熔池极限碳含量w[C]0,准确性更高。

Description

一种转炉炼钢烟气分析定碳方法
技术领域
本发明涉及钢铁冶金领域,特别是指一种转炉炼钢烟气分析定碳方法。
背景技术
氧气转炉炼钢是目前世界上最主要的炼钢方法,我国转炉产钢量占总产钢量的90%以上,脱碳反应贯穿转炉炼钢过程的始终。转炉吹炼是否到达冶炼终点的一个关键指标就是转炉炼钢熔池钢水碳含量是否满足出钢要求;在这方面,目前各钢厂大中型转炉广泛采用副枪技术来获得冶炼终点的钢水碳含量。然而,采用副枪不能连续获得转炉熔池内如脱碳速率等关键冶炼信息,并且副枪技术运行成本较高。随着计算机和各类检测技术的不断进步,转炉烟气分析技术由于不受转炉炉口尺寸大小限制的优势,能连续监测炉内炼钢反应进程,且对转炉节约成本、缩短冶炼周期以及冶炼终点碳-温预测等具有重要指导意义,因此近年来重新得到人们重视。
目前,大多数转炉烟气分析控制模型采用过程碳积分模型或终点碳曲线拟合模型。由于国内废钢管理过于粗放,废钢原料初始成分没有准确的信息,加上分析检测和称量误差等引起的铁水碳含量的误差,使得基于碳质量守恒的积分模型中初始碳含量的计算误差较大(有时可达±0.50%),而中低碳钢终点碳含量的控制精度要求为±0.02%。因此,只有消除初始碳含量的误差才能达到高精度控制。转炉冶炼终点碳曲线拟合模型重点关注吹炼后期,当前常见的终点碳曲线拟合模型有三次方模型和指数模型,这两种转炉冶炼终点碳曲线拟合模型虽然避开了初始碳含量难以精准确定的问题,但是并未考虑实际吹炼过程氧枪枪位、顶吹氧气流量以及底吹气体流量等操作工艺参数对转炉熔池脱碳速率的影响。
发明内容
本发明要解决的技术问题是提供一种转炉炼钢烟气分析定碳方法。该方法是在指数定碳模型的基础上引入了熔池混匀度的概念,完善了对冶炼过程熔池脱碳行为的表征,提高了计算的精度。所谓熔池混匀度,是指某一时刻转炉顶吹射流、底吹流股和熔池CO气泡对熔池的搅拌混匀强度,用来表征顶吹射流、底吹流股和熔池CO气泡对熔池脱碳反应的影响。通过水模型实验验证表明:一定范围内,枪位、顶吹气量、底吹气量和熔池CO气泡与熔池搅拌混匀强度成线性关系。该方法充分考虑了枪位、顶吹流量及底吹流量等操作工艺参数对熔池搅拌的影响,明显提高了在转炉冶炼末期碳含量预测的精确度。避免了基于碳质量守恒的积分模型对初始碳含量难以精准确定的问题,同时又改进了终点碳曲线拟合模型未考虑吹炼过程氧枪枪位、顶吹氧气流量以及底吹气体流量等操作工艺参数对转炉熔池脱碳速率的影响。
该方法包括以下步骤:
S1:依据质量守恒原理,结合转炉冶炼过程中检测设备获取的烟气CO、CO2百分含量和烟气流量计算出熔池脱碳速率;
S2:将熔池混匀度的概念引入到终点碳曲线拟合指数模型;
S3:根据冶金熔渣分子理论结合生产钢种特性,获取冶炼终渣成分和冶炼终点钢水温度,计算出冶炼熔池极限碳含量W[C]0
S4:得到转炉炼钢烟气分析定碳模型,得到转炉冶炼终点碳含量,实现对转炉冶炼终点钢水含碳量的预测。
其中:
S1中熔池脱碳速率是根据转炉冶炼过程碳平衡,结合烟气中CO、CO2百分含量和烟气流量计算得出,具体公式为:
Figure BDA0001638962510000021
式中,
Figure BDA0001638962510000022
表示脱碳速率,%·s-1;Qgas表示烟气流量,Nm3·s-1
Figure BDA0001638962510000023
Figure BDA0001638962510000024
分别表示烟气中CO、CO2百分含量,%;Wsteel表示钢液总质量,kg。
S4中转炉炼钢烟气分析定碳方法数学模型计算公式如下:
Figure BDA0001638962510000025
式中,w[C]0表示转炉熔池极限碳含量,%;h表示氧枪枪位,m;α表示常系数;hmin表示转炉冶炼过程氧枪的最低操作枪位,m;Qtop表示顶吹氧气流量,Nm3·s-1;Qtopmax表示转炉冶炼过程最大顶吹氧流量,Nm3·s-1;Qbottom表示底吹气体流量,Nm3·s-1;Qbottommax表示转炉冶炼过程最大底吹气体流量,Nm3·s-1
Figure BDA0001638962510000031
表示熔池脱碳速率;αk1和k2分为待定参数,在转炉冶炼后期,可通过非线性拟合的方法利用脱碳速率结合枪位、顶吹流量和底吹流量来求解。
S3中冶炼熔池极限碳含量W[C]0表示转炉冶炼过程钢水所能达到的最低碳含量,根据熔渣分子理论求解W[C]0如下所示:
Figure BDA0001638962510000032
式中,aFeO表示熔渣中FeO的活度,T表示冶炼终点出钢的平均温度。
本发明的上述技术方案的有益效果如下:
本发明引入熔池混匀度概念,充分考虑了氧枪枪位、顶吹氧气流量以及底吹气体流量操作工艺参数对转炉熔池脱碳过程的影响;且根据实际冶炼钢种利用熔渣活度计算得出熔池极限碳含量w[C]0,准确性更高。
附图说明
图1为本发明一种转炉炼钢烟气分析定碳方法的数学模型计算流程图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明提供一种转炉炼钢烟气分析定碳方法。
该方法包括以下步骤:
S1:依据质量守恒原理,结合转炉冶炼过程中检测设备获取的烟气CO、CO2百分含量和烟气流量计算出熔池脱碳速率;
S2:将熔池混匀度的概念引入到终点碳曲线拟合指数模型;
S3:根据冶金熔渣分子理论结合生产钢种特性,获取冶炼终渣成分和冶炼终点钢水温度,计算出冶炼熔池极限碳含量W[C]0
S4:得到转炉炼钢烟气分析定碳模型,得到转炉冶炼终点碳含量,实现对转炉冶炼终点钢水含碳量的预测。
如图1所示,在实际应用中,该方法的定碳模型的生成包括以下步骤:
(1)采集转炉炼钢历史生产数据,包含钢种、原材料等基础数据以及氧枪枪位、氧气流量等工艺参数过程数据,然后,对生产数据进行预处理,去除数据表中由于生产事故或人为等原因产生的异常数据;
(2)根据冶炼末期枪位、顶吹流量、底吹流量,烟气中CO、CO2百分含量和烟气流量,历史生产数据中终渣成分、终点出钢的平均温度,建立定碳模型;
(3)按照建立的定碳模型,重新开始吹炼,并对生产数据进行采集;
(4)计算判断终点碳含量是否合格,如不合格,将生产数据迭代至S1,以对定碳模型进行校准。
上述S1中熔池脱碳速率是根据转炉冶炼过程碳平衡,结合烟气中CO、CO2百分含量和烟气流量计算得出,具体公式为:
Figure BDA0001638962510000041
式中,
Figure BDA0001638962510000042
表示脱碳速率,%·s-1;Qgas表示烟气流量,Nm3·s-1
Figure BDA0001638962510000043
Figure BDA0001638962510000044
分别表示烟气中CO、CO2百分含量,%;Wsteel表示钢液总质量,kg。
在转炉吹炼后期,熔池的脱碳速率逐渐变慢,钢液碳含量降低变缓逐渐接近稳定状态,当碳含量降低到一定程度时,脱碳速率与碳含量w[C]之间有如下式:
Figure BDA0001638962510000045
整理得:
Figure BDA0001638962510000046
式中,w[C]0表示熔池极限碳含量,%;k1表示脱碳速率系数,k2表示氧气脱碳利用系数。
将熔池混匀度概念引入冶炼终点碳曲线拟合指数模型,通过水模型实验验证表明:在一定范围内,枪位、顶吹气量、底吹气量和熔池CO气泡对熔池搅拌混匀强度成线性关系。其表达式如下:
Figure BDA0001638962510000051
式中:η表示熔池混匀度,取值范围内0~100%;h表示氧枪枪位,m;α表示常系数;hmin表示转炉冶炼过程氧枪的最低操作枪位,m;Qtop表示顶吹氧气流量,Nm3·s-1;Qtopmax表示转炉冶炼过程最大顶吹氧流量,Nm3·s-1;Qbottom表示底吹气体流量,Nm3·s-1;Qbottommax表示转炉冶炼过程最大底吹气体流量,Nm3·s-1;εco表示转炉熔池内CO气体所占的百分含量,%。需要说明的是,转炉吹炼冶炼末期脱碳反应速率相对偏小,因此忽略该时期熔池内CO气泡对熔池搅拌脱碳速率的影响。
在式(3)中将k1替换为k1η,整理可得:
Figure BDA0001638962510000052
式中,αk1和k2分为待定参数。在转炉冶炼后期,可通过非线性拟合的方法利用脱碳速率结合枪位、顶吹流量和底吹流量来求解。
上述S4中转炉炼钢烟气分析定碳方法数学模型计算公式即为式(5)所示。
熔池极限碳含量w[C]0表示转炉冶炼过程钢水所能达到的最低碳含量,根据熔渣分子理论求解w[C]0如下所示:
转炉吹炼末期C与Fe元素的选择性氧化反应如式(6)所示:
[C]+(FeO)=COg+[Fe] (6)
Figure BDA0001638962510000053
在转炉吹炼末期,认为反应接近平衡状态,fC≈fO≈1,aFe=1,Pco/Pθ=1,由此,式(7)可化简为:
Figure BDA0001638962510000054
将生产钢种的冶炼终渣成分和冶炼终点钢水温度代入公式(8)中求出w[C]0,进而代入转炉炼钢烟气分析定碳模型,可得到转炉冶炼终点碳含量,实现对转炉冶炼终点钢水含碳量的预测。
本发明的具体实施例如下:
以某钢厂210吨转炉生产SPHC钢种时为实施载体,当转炉吹炼接近终点时,通过检测设备获取氧枪枪位、顶吹氧气流量、底吹气体流量、烟气CO和CO2百分含量以及烟气流量,根据冶金熔渣分子理论和历史数据,计算得到炉渣中FeO的活度为0.241。出钢平均温度为1686℃时,C、Fe元素选择性氧化的熔池极限碳含量w[C]0为0.033%。表1为SPHC钢的终渣成分,表2为本发明实施后的实验结果。从实验结果来看,该模型提供的转炉冶炼末期碳含量预报误差在±0.02%之间时有较高的命中率。
表1 SPHC钢的终渣成分,%
渣组元 CaO MgO MnO SiO<sub>2</sub> P<sub>2</sub>O<sub>5</sub> FeO
平均值 46.67 11.14 4.89 17.65 2.57 12.82
表2本发明实施后的试验结果
Figure BDA0001638962510000061
Figure BDA0001638962510000071
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (2)

1.一种转炉炼钢烟气分析定碳方法,其特征在于:包括以下步骤:
S1:依据质量守恒原理,结合转炉冶炼过程中检测设备获取的烟气CO、CO2百分含量和烟气流量计算出熔池脱碳速率;
S2:将熔池混匀度的概念引入到终点碳曲线拟合指数模型;
S3:根据冶金熔渣分子理论结合生产钢种特性,获取冶炼终渣成分和冶炼终点钢水温度,计算出冶炼熔池极限碳含量W[C]0
S4:得到转炉炼钢烟气分析定碳模型,得到转炉冶炼终点碳含量,实现对转炉冶炼终点钢水含碳量的预测;
所述S1中熔池脱碳速率根据转炉冶炼过程碳平衡,结合烟气中CO、CO2百分含量和烟气流量计算得出,具体公式为:
Figure FDA0002554199030000011
式中,
Figure FDA0002554199030000012
表示脱碳速率,%·s-1;Qgas表示烟气流量,Nm3·s-1
Figure FDA0002554199030000013
分别表示烟气中CO、CO2百分含量,%;Wsteel表示钢液总质量,kg;
所述S2中将熔池混匀度的概念引入到终点碳曲线拟合指数模型的具体方式为:
熔池混匀度的公式如下:
Figure FDA0002554199030000014
终点碳曲线拟合指数模型如下:
Figure FDA0002554199030000015
将熔池混匀度η对脱碳速率系数K1进行修正,得到ηK1代替碳曲线拟合模型中的K1,整理得到所述S4中转炉炼钢烟气分析定碳模型计算公式如下:
Figure FDA0002554199030000016
式中,w[C]0表示转炉熔池极限碳含量,%;h表示氧枪枪位,m;α表示常系数;hmin表示转炉冶炼过程氧枪的最低操作枪位,m;Qtop表示顶吹氧气流量,Nm3·s-1;Qtopmax表示转炉冶炼过程最大顶吹氧流量,Nm3·s-1;Qbottom表示底吹气体流量,Nm3·s-1;Qbottommax表示转炉冶炼过程最大底吹气体流量,Nm3·s-1;表示熔池脱碳速率;αk1和k2为待定参数,在转炉冶炼后期,通过非线性拟合
Figure FDA0002554199030000021
的方法利用脱碳速率结合枪位、顶吹流量和底吹流量来求解。
2.根据权利要求1所述的转炉炼钢烟气分析定碳方法,其特征在于:所述S3中冶炼熔池极限碳含量W[C]0表示转炉冶炼过程钢水所能达到的最低碳含量,根据熔渣分子理论求解W[C]0如下所示:
Figure FDA0002554199030000022
式中,aFeO表示熔渣中FeO的活度,T表示冶炼终点出钢的平均温度。
CN201810372680.8A 2018-04-24 2018-04-24 一种转炉炼钢烟气分析定碳方法 Active CN108647407B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810372680.8A CN108647407B (zh) 2018-04-24 2018-04-24 一种转炉炼钢烟气分析定碳方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810372680.8A CN108647407B (zh) 2018-04-24 2018-04-24 一种转炉炼钢烟气分析定碳方法

Publications (2)

Publication Number Publication Date
CN108647407A CN108647407A (zh) 2018-10-12
CN108647407B true CN108647407B (zh) 2020-08-25

Family

ID=63747025

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810372680.8A Active CN108647407B (zh) 2018-04-24 2018-04-24 一种转炉炼钢烟气分析定碳方法

Country Status (1)

Country Link
CN (1) CN108647407B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111079537B (zh) * 2019-11-18 2023-09-26 中冶赛迪技术研究中心有限公司 转炉冶炼工况的识别方法、系统、机器可读介质及设备
CN110991089B (zh) * 2019-12-26 2020-07-31 北京科技大学 一种转炉炼钢冶炼后期碳含量预报方法
CN113076505B (zh) * 2020-01-05 2024-06-04 上海梅山钢铁股份有限公司 一种转炉钢水脱碳速率计算方法
CN113512620B (zh) * 2020-04-10 2024-04-26 北京凯德恒源科技发展有限公司 气体分析+副枪的转炉冶炼全过程终点碳动态控制方法
CN111518980B (zh) * 2020-04-23 2021-11-30 中冶南方工程技术有限公司 一种转炉终点碳含量预报模型的修正方法和系统
CN112680566B (zh) * 2020-12-18 2022-10-21 北京首钢自动化信息技术有限公司 一种精炼炉脱碳终点检测方法及检测系统
WO2022198594A1 (zh) * 2021-03-25 2022-09-29 北京凯德恒源科技发展有限公司 基于气体分析的转炉冶炼全过程终点碳动态控制方法
CN113462846B (zh) * 2021-06-25 2022-11-11 马鞍山钢铁股份有限公司 一种副枪定碳故障应急处理方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101713010A (zh) * 2009-11-17 2010-05-26 邯钢集团邯宝钢铁有限公司 一种利用数学模型指导回收转炉煤气的方法
CN102344986A (zh) * 2011-11-11 2012-02-08 田陆 转炉炼钢终点控制的方法、装置和系统
CN102399933A (zh) * 2010-09-07 2012-04-04 鞍钢股份有限公司 一种转炉吹炼低碳钢氧枪自动控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6440355B1 (en) * 2000-09-06 2002-08-27 Bethlehem Steel Corporation Apparatus for measuring bath level in a basic oxygen furnace to determine lance height adjustment

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101713010A (zh) * 2009-11-17 2010-05-26 邯钢集团邯宝钢铁有限公司 一种利用数学模型指导回收转炉煤气的方法
CN102399933A (zh) * 2010-09-07 2012-04-04 鞍钢股份有限公司 一种转炉吹炼低碳钢氧枪自动控制方法
CN102344986A (zh) * 2011-11-11 2012-02-08 田陆 转炉炼钢终点控制的方法、装置和系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
在150t转炉上用质谱仪进行钢水连续定碳;胡志刚 等;《钢铁研究学报》;20030815;第15卷(第4期);第62-63页第2.1-2.2节 *
基于烟气成分分析的电弧炉炼钢脱碳模型;何春来 等;《北京科技大学学报》;20101123;第32卷(第12期);第1539页第2.2节 *
转炉冶炼中高碳钢过程及终点控制研究;汪宙;《中国博士学位论文全文数据库 工程科技I辑》;20160515(第5期);第55-70页第四章第4.1.3-4.3.2节及图4-7 *

Also Published As

Publication number Publication date
CN108647407A (zh) 2018-10-12

Similar Documents

Publication Publication Date Title
CN108647407B (zh) 一种转炉炼钢烟气分析定碳方法
CN109447346B (zh) 基于灰色预测与神经网络组合模型转炉氧耗量预测方法
CN102344986B (zh) 转炉炼钢终点控制的方法、装置和系统
WO2018012257A1 (ja) 溶鋼中りん濃度推定方法及び転炉吹錬制御装置
CN104630410A (zh) 一种基于数据解析的转炉炼钢质量实时动态预测方法
CN113919206A (zh) 一种炼钢合金化窄成分控制方法
US20150337404A1 (en) Method and device for predicting, controlling and/or regulating steelworks processes
CN106834593A (zh) 一种以参考炉次法确定rh精炼炉脱碳工艺数据的方法
CN101308154A (zh) 一种转炉炼钢钢水连续定碳方法
CN101592650B (zh) 电炉炼钢炉内钢水碳含量的连续测定方法
CN115522012B (zh) 一种大型转炉控制tsc阶段吹氧量的方法
CN113192568A (zh) 一种精炼炉脱硫终点预报方法和系统
CN114854929A (zh) 一种转炉喷吹co2钢液成分与温度实时动态预测的方法
CN101403567A (zh) 基于svm的电弧炉终点温度预报系统
CN113961865A (zh) 一种大型转炉精确控制tsc阶段调温剂加入量的方法
JP2012136767A (ja) 転炉りん濃度推定方法
CN115927784B (zh) 一种基于co2动态混合喷吹的转炉炼钢终点控制方法
CN111518980B (zh) 一种转炉终点碳含量预报模型的修正方法和系统
CN103276136A (zh) 一种基于副枪系统的转炉炼钢钢水定磷方法
Mao-qiang et al. Dynamic, data-driven prediction model of molten steel temperature in the second blowing stage of a converter
CN108588324B (zh) 转炉炼钢通过CO2高强度输入控制渣中(FeO)和粉尘产生的方法
CN113736948B (zh) 一种dc04钢转炉冶炼终点不等样的脱磷控制方法
CN101592651B (zh) 电炉炼钢炉内钢水碳含量的测定方法
CN104846150A (zh) 低铝模铸钢冶炼方法
CN113076505B (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