CN105929264A - 一种750kV输电线路雷电绕击性能评估方法 - Google Patents
一种750kV输电线路雷电绕击性能评估方法 Download PDFInfo
- Publication number
- CN105929264A CN105929264A CN201610237831.XA CN201610237831A CN105929264A CN 105929264 A CN105929264 A CN 105929264A CN 201610237831 A CN201610237831 A CN 201610237831A CN 105929264 A CN105929264 A CN 105929264A
- Authority
- CN
- China
- Prior art keywords
- transmission line
- electricity
- thunderbolt
- span
- line
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Suspension Of Electric Lines Or Cables (AREA)
Abstract
本发明公开了一种750kV输电线路雷电绕击性能评估方法,包括以下步骤:1)获取进行雷电绕击性能校核的750kV输电线路信息;2)将750kV输电线路的档距划分为若干分段,得各分段处导地线的相对几何位置信息;3)给定750kV输电线路的地理信息;4)将750kV输电线路中的导线工作电压相位划分为若干区段,得到不同相位导线工作电压下各分段线路的雷电绕击跳闸率;5)不同相位导线工作电压下各分段线路的雷电绕击跳闸率确定输电线路的雷电绕击薄弱点,完成750kV输电线路雷电绕击性能评估。本发明能够准确的评估出750kV输电线路雷电绕击性能。
Description
技术领域
本发明属于750kV电力工程仿真计算领域,涉及一种750kV输电线路雷电绕击性能评估方法。
背景技术
ATP-EMTP程序是国际通用的先进的图形化的电磁暂态计算程序。该程序由世界各国参与研发,程序的准确性得到世界公认,并由IEC71-4导则予以推荐和确认,国内外研究机构均采用该程序进行电磁暂态的仿真计算。ATP-EMTP的基本算法:根据元件的不同特性,建立相应的代数方程,常微分方程和偏微分方程,利用梯形积分法将电感、电容、电源等集中参数元件化成电阻性网络,对于传输线等分布性参数利用其上的波过程的特征线方程,经过一定的转换,把分布参数的线段也等效为电阻性网络,则其相应的方程也变为代数方程,进一步形成节点导纳矩阵;然后采用优化结点编号技术和稀疏矩阵算法,以节点电压为未知量,利用矩阵三角分解求解,最后求得各支路的电流、电压和所有消耗的功率、能量。在稳态计算中应将非线性元件线性化,包括利用简单的迭代进行潮流计算。
电气几何模型(EGM)是当今国际上流行的分析绕击率的方法。它以雷击机理的现代知识作为基础。其基本原理建立在下列基本概念和假设基础上。
(1)由雷云向地面发展的先导通道头部到达被击物体的临界击穿距离-击距以前,击中点是不确定的。先到哪个物体的击距内,即向该物体放电。
(2)击距rs是雷电流幅值I的函数。各国采用此函数关系不完全相同。考虑了下列4种函数关系式,取其最严重者,并考虑10%误差。
(a)rs=7.1I0.75(Whitehead)
(b)rs=8I0.65(IEEE工作组)
(c)rs=10I0.65(美国)
(d)rs=6.72I0.8(日本)
(3)先导接近地面时的入射角ψ服从某一给定的概率分布函数。其分布概率
图1为雷击线路的电气几何模型,其中,S为地线,C为导线,α为保护角,对于不同的雷电流幅值Ii,由上述关系式可算出相应的击距rsi。当先导头部进入BiCi弧面,放电将击向地线。当先导头部进入CiDi弧面则击中导线,即发生绕击。所以CiDi弧面称为暴露面。当先导头部进入DiEi平面,则击中大地。
随着雷电流幅值的增大,暴露弧CiDi逐渐缩小。当雷电流增大到Im时暴露弧缩小为零。此时雷或击中地线,或击中大地,不再发生绕击。Im称为最大绕击电流。相应的击距称为最大击距Rsm,一般情况下击距并非所有的绕击都会引起绝缘的闪络,只有当雷电流在导线上引起的电压Uc(包括工作电压)大于绝缘放电电压时,即Uc=I0Zc/2.2+Em≥U50才会闪络,I0为最小危险绕击电流。对应雷电流I及击距rs,在一定的rs下,暴露弧在地面的投影见,参考图2,线路地平面上(水平方向,地面倾斜角为零时)的相应暴露面积为:
考虑不同雷电流的出现概率,即
其中,P(I)为雷电流I的概率分布密度,为危险绕击率,式中W为引雷宽度。当地面倾斜角不为零时,其计算要复杂些,但基本原理是相同的。
近年来,相关的研究院所及高校对750kV线路的雷电绕击性能进行了仿真研究,但其提出的仿真计算方法大多是在一般线路的计算方法上进行杆塔模型的修正,且未形成完整、统一的仿真方法,且在仿真计算中没有考虑多种因素下线路雷电绕击性能影响,然而仿真计算得到结果精度不高,其雷电性能评估结论也未应用至750kV线路的设计中。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供了一种750kV输电线路雷电绕击性能评估方法,该方法能够准确的评估出750kV输电线路雷电绕击性能。
为达到上述目的,本发明所述的750kV输电线路雷电绕击性能评估方法包括以下步骤:
1)获取进行雷电绕击性能校核的750kV输电线路信息,所述750kV输电线路信息包括750kV输电线路的基本信息、750kV输电线路的结构特征及750kV输电线路的绝缘特征;
2)将750kV输电线路的档距划分为若干分段,再根据750kV输电线路的基本信息及750kV输电线路的结构特征得各分段处导地线的相对几何位置信息;
3)根据将750kV输电线路的档距划分为若干分段的结果给定750kV输电线路的地理信息;
4)将750kV输电线路中的导线工作电压相位划分为若干区段,再根据步骤2)得到的各分段处导地线相对几何位置信息以及步骤3)得到的750kV输电线路的地理信息、750kV输电线路的绝缘特征计算得到不同相位导线工作电压下各分段线路的雷电绕击跳闸率,其中,
不同相位导线工作电压下各分段线路的雷电绕击跳闸率为各分段线路在不同导线工作电压下雷电绕击跳闸率的加权平均值;
不同相位导线工作电压下各档距的雷电绕击跳闸率为对应档距内各分段线路雷电绕击跳闸率的加权平均值;
750kV输电线路的雷电绕击跳闸率为750kV输电线路中所有档距雷电绕击跳闸率的加权平均值;
5)根据步骤4)得到的各分段线路在不同导线工作电压下雷电绕击跳闸率的加权平均值、各档距内各分段线路雷电绕击跳闸率的加权平均值以及750kV输电线路中所有档距雷电绕击跳闸率的加权平均值得各档距的雷电绕击跳闸率相对偏差,然后根据各档距的雷电绕击跳闸率相对偏差确定输电线路的雷电绕击薄弱点,完成750kV输电线路雷电绕击性能评估。
输电线路的基本信息包括输电线路的名称、电压等级、总长度及各档距线路内的分段数;
输电线路的结构特征包括输电线路中杆塔塔型、导地线型号、导地线弧垂、档距长度及绝缘子串结构串长。
输电线路的绝缘特征包括绝缘子串闪络电压及最短空气间隙闪络电压。
输电线路地理信息包括各杆塔处的海拔高度、沿线各分段处的海拔高度、垂直线路方向上距杆塔不同距离处的海拔高度及垂直线路方向上距沿线各分段不同距离位置处的海拔高度。
步骤2)中将750kV输电线路的档距划分为10个分段。
步骤4)中导线工作电压相位的变化范围为0°-360°。
本发明具有以下有益效果:
本发明所述的750kV输电线路雷电绕击性能评估方法在操作时,针对750kV输电线路杆塔的塔高较高、导线跨越档距较长的特点,将750kV输电线路的档距划分为若干分段,使750kV输电线路进行细分,提高地线跟杆塔连接时在工频电压作用下的环流电流仿真计算精度,同时将导线工作电压相位划分为若干区段,计算不同相位导线工作电压下各分段线路的雷电绕击跳闸率,然后再根据不同相位导线工作电压下各分段线路的雷电绕击跳闸率确定输电线路的雷电绕击薄弱点,完成750kV输电线路雷电绕击性能评估。需要说明的是,本发明对750kV输电线路及导线工作电压相位进行细分,从而有效的提高仿真的精度,实现750kV输电线路雷电绕击性能准确评估,计算结果能够适用于指导工程的科学设计。
附图说明
图1为雷击线路的电气几何模型;
图2为计算雷绕击线路导线概率的电气几何模型示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
本发明所述的750kV输电线路雷电绕击性能评估方法包括以下步骤:
1)获取进行雷电绕击性能校核的750kV输电线路信息,所述750kV输电线路信息包括750kV输电线路的基本信息、750kV输电线路的结构特征及750kV输电线路的绝缘特征,其中,输电线路的基本信息包括输电线路的名称、电压等级、总长度及各档距线路内的分段数;输电线路的结构特征包括输电线路中杆塔塔型、导地线型号、导地线弧垂、档距长度及绝缘子串结构串长;输电线路的绝缘特征包括绝缘子串闪络电压及最短空气间隙闪络电压;
2)将750kV输电线路的档距划分为若干分段,再根据750kV输电线路的基本信息及750kV输电线路的结构特征得各分段处导地线的相对几何位置信息,其中,将750kV输电线路的档距划分为10个分段;
3)根据将750kV输电线路的档距划分为若干分段的结果给定750kV输电线路的地理信息,其中,输电线路地理信息包括各杆塔处的海拔高度、沿线各分段处的海拔高度、垂直线路方向上距杆塔不同距离处的海拔高度及垂直线路方向上距沿线各分段不同距离位置处的海拔高度;
4)将750kV输电线路中的导线工作电压相位划分为若干区段,其中,导线工作电压相位的变化范围为0°-360°,再根据步骤2)得到的各分段处导地线相对几何位置信息以及步骤3)得到的750kV输电线路的地理信息、750kV输电线路的绝缘特征计算得到不同相位导线工作电压下各分段线路的雷电绕击跳闸率,其中,
不同相位导线工作电压下各分段线路的雷电绕击跳闸率为各分段线路在不同导线工作电压下雷电绕击跳闸率的加权平均值;
不同相位导线工作电压下各档距的雷电绕击跳闸率为对应档距内各分段线路雷电绕击跳闸率的加权平均值;
750kV输电线路的雷电绕击跳闸率为750kV输电线路中所有档距雷电绕击跳闸率的加权平均值;
5)根据步骤4)得到的各分段线路在不同导线工作电压下雷电绕击跳闸率的加权平均值、各档距内各分段线路雷电绕击跳闸率的加权平均值以及750kV输电线路中所有档距雷电绕击跳闸率的加权平均值得各档距的雷电绕击跳闸率相对偏差,然后根据各档距的雷电绕击跳闸率相对偏差确定输电线路的雷电绕击薄弱点,完成750kV输电线路雷电绕击性能评估。
需要说明的是,本发明针对750kV输电线路杆塔的塔高较高、导线跨越档距较长、杆塔引起雷电绕击概率较高、且发生雷击事故后不易维修、更换电力设备的特点。提出了调整工作电压相位角使得雷电绕击线路时所有工频电压幅值涵盖在仿真计算中,更精确的评估输电线路的雷电绕击性能,使计算结果适用于指导工程的科学设计。同时针对750kV输电线路的杆塔高度较高的特点,将杆塔各部分、各段的不同波阻抗进行细分,提高地线跟杆塔连接时在工频电压作用下的环流电流仿真计算精度,最后根据评估的结果找出整条750kV输电线路中防雷薄弱点,并提出不同地区、不同条件下,750kV输电线路应采用的差异化防雷措施,优化工程设计,节省工程投资。
Claims (7)
1.一种750kV输电线路雷电绕击性能评估方法,其特征在于,包括以下步骤:
1)获取进行雷电绕击性能校核的750kV输电线路信息,所述750kV输电线路信息包括750kV输电线路的基本信息、750kV输电线路的结构特征及750kV输电线路的绝缘特征;
2)将750kV输电线路的档距划分为若干分段,再根据750kV输电线路的基本信息及750kV输电线路的结构特征得各分段处导地线的相对几何位置信息;
3)根据将750kV输电线路的档距划分为若干分段的结果给定750kV输电线路的地理信息;
4)将750kV输电线路中的导线工作电压相位划分为若干区段,再根据步骤2)得到的各分段处导地线相对几何位置信息以及步骤3)得到的750kV输电线路的地理信息、750kV输电线路的绝缘特征计算得到不同相位导线工作电压下各分段线路的雷电绕击跳闸率,其中,
不同相位导线工作电压下各分段线路的雷电绕击跳闸率为各分段线路在不同导线工作电压下雷电绕击跳闸率的加权平均值;
不同相位导线工作电压下各档距的雷电绕击跳闸率为对应档距内各分段线路雷电绕击跳闸率的加权平均值;
750kV输电线路的雷电绕击跳闸率为750kV输电线路中所有档距雷电绕击跳闸率的加权平均值;
5)根据步骤4)得到的各分段线路在不同导线工作电压下雷电绕击跳闸率的加权平均值、各档距内各分段线路雷电绕击跳闸率的加权平均值以及750kV输电线路中所有档距雷电绕击跳闸率的加权平均值得各档距的雷电绕击跳闸率相对偏差,然后根据各档距的雷电绕击跳闸率相对偏差确定输电线路的雷电绕击薄弱点,完成750kV输电线路雷电绕击性能评估。
2.根据权利要求1所述的750kV输电线路雷电绕击性能评估方法,其特征在于,输电线路的基本信息包括输电线路的名称、电压等级、总长度及各档距线路内的分段数。
3.根据权利要求1所述的750kV输电线路雷电绕击性能评估方法,其特征在于,输电线路的结构特征包括输电线路中杆塔塔型、导地线型号、导地线弧垂、档距长度及绝缘子串结构串长。
4.根据权利要求1所述的750kV输电线路雷电绕击性能评估方法,其特征在于,输电线路的绝缘特征包括绝缘子串闪络电压及最短空气间隙闪络电压。
5.根据权利要求1所述的750kV输电线路雷电绕击性能评估方法,其特征在于,输电线路地理信息包括各杆塔处的海拔高度、沿线各分段处的海拔高度、垂直线路方向上距杆塔不同距离处的海拔高度及垂直线路方向上距沿线各分段不同距离位置处的海拔高度。
6.根据权利要求1所述的750kV输电线路雷电绕击性能评估方法,其特征在于,步骤2)中将750kV输电线路的档距划分为10个分段。
7.根据权利要求1所述的750kV输电线路雷电绕击性能评估方法,其特征在于,步骤4)中导线工作电压相位的变化范围为0°-360°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610237831.XA CN105929264B (zh) | 2016-04-15 | 2016-04-15 | 一种750kV输电线路雷电绕击性能评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610237831.XA CN105929264B (zh) | 2016-04-15 | 2016-04-15 | 一种750kV输电线路雷电绕击性能评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105929264A true CN105929264A (zh) | 2016-09-07 |
CN105929264B CN105929264B (zh) | 2018-08-03 |
Family
ID=56839305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610237831.XA Active CN105929264B (zh) | 2016-04-15 | 2016-04-15 | 一种750kV输电线路雷电绕击性能评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105929264B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107703381A (zh) * | 2017-09-19 | 2018-02-16 | 国网浙江省电力公司 | 一种邻近输电线路间雷击屏蔽效应的计算方法 |
CN109521333A (zh) * | 2018-04-08 | 2019-03-26 | 国网甘肃省电力公司检修公司 | 计及导线间屏蔽作用的同塔双回线路绕击耐雷性能分析方法 |
CN111814355A (zh) * | 2020-07-24 | 2020-10-23 | 南方电网科学研究院有限责任公司 | 高压输电线路雷电绕击的防护方法、装置及终端设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030069657A1 (en) * | 2001-10-09 | 2003-04-10 | Verdecchio Thomas A. | Portable protective air gap tool and method |
JP2004012292A (ja) * | 2002-06-06 | 2004-01-15 | Tokyo Electric Power Co Inc:The | 送電線の故障点標定方法およびそれを用いた故障点標定システム |
CN101093240A (zh) * | 2006-06-23 | 2007-12-26 | 郭玉章 | 电力输电线路雷击点的确定方法 |
CN101216864A (zh) * | 2008-01-09 | 2008-07-09 | 国网武汉高压研究院 | 大跨越线路的总雷击跳闸率仿真计算方法 |
CN101290336A (zh) * | 2008-06-18 | 2008-10-22 | 昆明理工大学 | 一种交流输电线路雷电绕击与反击的识别方法 |
CN103488815A (zh) * | 2013-08-26 | 2014-01-01 | 国家电网公司 | 一种输电线路雷电绕击风险评估方法 |
CN103646148A (zh) * | 2013-12-20 | 2014-03-19 | 国家电网公司 | 一种特高压输电线路雷电反击性能仿真计算方法 |
CN103810339A (zh) * | 2014-02-14 | 2014-05-21 | 国家电网公司 | 一种特高压输电线路大跨越段雷电反击性能的确定方法 |
-
2016
- 2016-04-15 CN CN201610237831.XA patent/CN105929264B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030069657A1 (en) * | 2001-10-09 | 2003-04-10 | Verdecchio Thomas A. | Portable protective air gap tool and method |
JP2004012292A (ja) * | 2002-06-06 | 2004-01-15 | Tokyo Electric Power Co Inc:The | 送電線の故障点標定方法およびそれを用いた故障点標定システム |
CN101093240A (zh) * | 2006-06-23 | 2007-12-26 | 郭玉章 | 电力输电线路雷击点的确定方法 |
CN101216864A (zh) * | 2008-01-09 | 2008-07-09 | 国网武汉高压研究院 | 大跨越线路的总雷击跳闸率仿真计算方法 |
CN101290336A (zh) * | 2008-06-18 | 2008-10-22 | 昆明理工大学 | 一种交流输电线路雷电绕击与反击的识别方法 |
CN103488815A (zh) * | 2013-08-26 | 2014-01-01 | 国家电网公司 | 一种输电线路雷电绕击风险评估方法 |
CN103646148A (zh) * | 2013-12-20 | 2014-03-19 | 国家电网公司 | 一种特高压输电线路雷电反击性能仿真计算方法 |
CN103810339A (zh) * | 2014-02-14 | 2014-05-21 | 国家电网公司 | 一种特高压输电线路大跨越段雷电反击性能的确定方法 |
Non-Patent Citations (9)
Title |
---|
A. J. ERIKSSON: "An Improved Electrogeometric Model for Transmission Line Shielding Analysis", 《IEEE POWER ENGINEERING REVIEW》 * |
刘勇: "区域输电网雷击闪络风险评估技术研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
司马文霞 等: "复杂地形下输电线路绕击跳闸次数3维计算方法", 《高电压技术》 * |
司马文霞 等: "输电线路雷电绕击评估方法及展望", 《高电压技术》 * |
宗天元: "蒙特卡洛法计算750kV输电线路绕击跳闸率的研究", 《宁夏电力》 * |
曾嵘 等: "输电线路雷电防护技术研究)二):分析方法", 《高电压技术》 * |
柯维新 等: "同塔四回输电线路的雷击跳闸率计算方法研究", 《电力学报》 * |
袁海燕 等: "导线工作电压对特高压输电线路雷电绕击性能的影响", 《中国电机工程学会高压专委会2007年学术年会论文集》 * |
邓红雷 等: "基于层次化、差异化的架空输电线路雷击风险评估", 《电力系统保护与控制》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107703381A (zh) * | 2017-09-19 | 2018-02-16 | 国网浙江省电力公司 | 一种邻近输电线路间雷击屏蔽效应的计算方法 |
CN107703381B (zh) * | 2017-09-19 | 2020-04-21 | 国网浙江省电力公司 | 一种邻近输电线路间雷击屏蔽效应的计算方法 |
CN109521333A (zh) * | 2018-04-08 | 2019-03-26 | 国网甘肃省电力公司检修公司 | 计及导线间屏蔽作用的同塔双回线路绕击耐雷性能分析方法 |
CN111814355A (zh) * | 2020-07-24 | 2020-10-23 | 南方电网科学研究院有限责任公司 | 高压输电线路雷电绕击的防护方法、装置及终端设备 |
CN111814355B (zh) * | 2020-07-24 | 2024-04-26 | 南方电网科学研究院有限责任公司 | 高压输电线路雷电绕击的防护方法、装置及终端设备 |
Also Published As
Publication number | Publication date |
---|---|
CN105929264B (zh) | 2018-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101900779B (zh) | 一种结合地形地貌测量高压输电线路雷电绕击跳闸率的方法 | |
CN102435921B (zh) | 同塔双回输电线路绝缘及耐雷电冲击性能的判定方法 | |
CN102841280B (zh) | 500kV同塔四回路输电线雷击跳闸率仿真方法 | |
CN105242133B (zh) | 一种改进配电线路雷电跳闸率计算方法 | |
CN108693444A (zh) | 一种基于电磁信号时间反演的输电线路故障点定位方法 | |
CN100361123C (zh) | 判定输电线路防雷性能的全线路、多参数综合优化方法 | |
CN103837769B (zh) | 一种输电线路雷害预警方法及其系统 | |
CN109753703A (zh) | 一种配电线路防雷等级评估方法 | |
CN104573245B (zh) | 一种特高压输电线路中地线损耗的计算方法 | |
CN105224797B (zh) | 一种特高压大跨越线路反击跳闸率计算方法 | |
CN103293419A (zh) | 一种接地装置冲击特性的评估方法 | |
CN105868872A (zh) | 一种配电网雷电灾害故障预测方法 | |
CN103488815A (zh) | 一种输电线路雷电绕击风险评估方法 | |
CN112257238B (zh) | 多重雷下的断路器绝缘受损风险评估方法及其装置 | |
CN105929264A (zh) | 一种750kV输电线路雷电绕击性能评估方法 | |
CN111695775A (zh) | 基于偏好比率-熵值法最优组合赋权的输电线路雷击风险评估方法 | |
CN102298105A (zh) | 一种线路雷电绕击性能校核方法 | |
CN102680837B (zh) | 输电线路周边配电线路感应雷跳闸风险评估方法与装置 | |
CN107703381B (zh) | 一种邻近输电线路间雷击屏蔽效应的计算方法 | |
CN110687371B (zh) | 一种用于确定同塔多回线路的雷电绕击性能的方法及系统 | |
CN117291418A (zh) | 一种基于分布式监测数据的线路雷害风险评估方法及装置 | |
CN105279342A (zh) | 地铁接触网导线的分段方法和系统、仿真方法和系统 | |
CN107783017A (zh) | 一种交流输电导线电晕损失计算方法 | |
Filipe et al. | A methodology for estimating transmission lines lightning performance using a statistical approach | |
Trotsenko et al. | Overview of factors affecting the estimation of lightning shielding performance of overhead transmission lines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |