CN109270333A - 一种识别雷击引起的电压暂降的方法 - Google Patents
一种识别雷击引起的电压暂降的方法 Download PDFInfo
- Publication number
- CN109270333A CN109270333A CN201811296162.9A CN201811296162A CN109270333A CN 109270333 A CN109270333 A CN 109270333A CN 201811296162 A CN201811296162 A CN 201811296162A CN 109270333 A CN109270333 A CN 109270333A
- Authority
- CN
- China
- Prior art keywords
- voltage dip
- lightning
- voltage
- struck
- caused
- 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
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
-
- 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
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
本发明公开一种识别雷击引起的电压暂降的方法,包括步骤:S1:从变电站的电能质量监测装置获取母线电压信号;S2:计算有效值电压,判断电压暂降的起始点和结束点;S3:计算电压暂降幅值Vsag,若Vsag大于阈值且电压暂降结束点有高频分量,则为其它扰动原因;若Vsag小于阈值或电压暂降结束点没有高频分量,则进入步骤S4;S4:对母线电压信号进行小波变换得到细节分量并计算信号能量;S5:若有高频段的信号能量,则为雷击引起的故障;若没有高频段的信号能量,则为普通短路故障;S6:判断结束。本发明利用现有的电能质量监测装置采集的电压信号就能进行识别,无需增设行波检测装置或其他检测装置,节约了大量成本,更好利用电能质量监测装置采集的数据。
Description
技术领域
本发明涉及输电线路故障识别领域,更具体地,涉及一种识别雷击引起的电压暂降的方法。
背景技术
电压暂降是电压有效值在短时间内减小,是由短路、过负荷和大型电机启动等引起的电压扰动事件。其中,系统故障是导致电压暂降的最主要原因,而雷击故障在电网故障中所占比例很大,是造成母线电压暂降的重要原因。不同扰动源引起的电压暂降的持续时间和暂降幅值不同,对用户和电网的影响也不一样,对应的补偿措施也不一样。因此,准确的电压暂降源识别对于改善电能质量具有重要意义。
雷击引起的故障和普通短路故障引起的电压暂降有所差异,造成的影响也不相同。大量的学者对此展开研究,提出了许多识别雷击和普通短路故障的方法。但不论是基于电流行波信号,还是电压行波信号,都有一定的局限性且需装设专门的行波检测装置,成本较高技术复杂。
现有技术主要基于S变换、小波变换、HHT变换来区分雷电冲击与普通故障,其中S变换计算过程较为复杂,速度慢,比较适合仿真计算;小波变换和HHT变换速度较快,但现有技术又多是基于电流波形对雷击引起的故障进行识别,与如今实际电网中的电能质量监测系统主要采集的是电压数据相冲突。因此需要一种基于电能质量监测装置采集的电压数据采用小波变换或HHT变换的雷击识别新方法。
发明内容
本发明为克服上述现有技术所述的至少一种缺陷,提供一种识别雷击引起的电压暂降的方法。
本发明旨在至少在一定程度上解决上述技术问题。
本发明的首要目的是利用现有的电能质量监测装置采集的电压信号就能进行识别,无需增设行波检测装置或其他检测装置,节约了大量成本,也更好的利用了电能质量监测装置采集的数据。
为解决上述技术问题,本发明的技术方案如下:
一种识别雷击引起的电压暂降的方法,包括以下步骤:
S1:从变电站的电能质量监测装置获取母线电压信号;
S2:计算有效值电压,判断电压暂降的起始点和结束点;
S3:计算电压暂降幅值Vsag,若Vsag大于阈值且电压暂降结束点有高频分量,则为其它扰动原因,进入步骤S6;若Vsag小于阈值或电压暂降结束点没有高频分量,则进入步骤S4;
S4:对电压暂降段的母线电压信号进行小波变换得到各频段的细节分量并计算各频段的信号能量;
S5:若有高频段的信号能量,则为雷击引起的故障;若没有高频段的信号能量,则为普通短路故障;
S6:判断结束;
上述方法利用雷击与短路故障造成的电压信号经小波变换后能量分布的巨大差异分辨出雷击与短路故障,更加直观,简单高效,同时对高频分量出现时刻进行检测与小波变换相结合的方法进行雷击识别,与其它扰动原因分开,同时小波变换计算速度快,简单高效,提高了雷击识别效率。
优选地,步骤S2中判断电压暂降的起始点和结束点的具体依据为:
a)有效值电压降至0.9p.u.时的一刻为电压暂降的起始点;
b)有效值电压升至0.9p.u.时的一刻为电压暂降的结束点。
其中标幺值(p.u.)是电力系统分析和工程计算中常用的数值标记方法,表示各物理量及参数的相对值。
优选地,步骤S3中所述阈值为0.7p.u.。
优选地,步骤S4的小波变换为db6小波6层变换,得到细节分量D1、D2、D3、D4、D5、D6,其频率大小为D6<D5<D4<D3<D2<D1。
优选地,步骤S4计算各频段的信号能量的具体计算方法为:
定义单一尺度的小波能量为该尺度下重构信号的平方和:
式中,Ej为Dj细节分量的信号能量,Cj(k)为小波系数;
则信号总能量为
各频带能量所占比重得到向量组P=[P1 P2 P3 P4 P5 P6]。
优选地,步骤S5中若有高频段的信号能量,则为雷击引起的故障;若没有高频段的信号能量,则为普通短路故障的具体步骤为:
若P1>0.1且P2>0.1,则为雷击引起的故障;若P1≤0.1或P2≤0.1,则为普通短路故障;
由于雷击引起的故障的能量在低频带(P6)和高频带(P1、P2)都存在,而普通故障的能量集中在低频带(P6)。其中雷击时P1、P2的最小值都在0.1以上,而普通故障时P1、P2的最大值仍为0,基于这种差异,可以分辨雷击和普通故障引起的电压暂降。
优选地,步骤S3中其它扰动原因为变压器激磁。
与现有技术相比,本发明技术方案的有益效果是:
本发明利用现有的电能质量监测装置采集的电压信号就能进行识别,无需增设行波检测装置或其他检测装置,节约了大量成本,也更好的利用了电能质量监测装置采集的数据。利用小波变换后能量分布的巨大差异就能够准确的识别出雷击与短路故障,更加直观,简单高效。采用高频分量出现时刻的检测与小波变换相结合的方法进行雷击识别,其中小波变换计算速度快,简单高效,提高了雷击识别效率。
附图说明
图1为本发明提供的一种识别雷击引起的电压暂降的方法流程图。
图2为雷击与普通短路故障母线电压暂降波形对比图。
图3为雷击引起的电压暂降波形经小波变换后得到的D1细节分量图。
图4为变压器激磁引起的电压暂降波形经小波变换后得到的D1细节分量图形。
图5为雷击引起的能量分布图。
图6为普通短路的电压暂降的能量分布图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
一种识别雷击引起的电压暂降的方法,如图1,包括以下步骤:
S1:从变电站的电能质量监测装置获取母线电压信号;
S2:计算有效值电压,判断电压暂降的起始点和结束点;
S3:计算电压暂降幅值Vsag,若Vsag大于阈值且电压暂降结束点有高频分量,则为其它扰动原因,进入步骤S6;若Vsag小于阈值或电压暂降结束点没有高频分量,则进入步骤S4;
S4:对电压暂降段的母线电压信号进行小波变换得到各频段的细节分量并计算各频段的信号能量;
S5:若有高频段的信号能量,则为雷击引起的故障;若没有高频段的信号能量,则为普通短路故障;
S6:判断结束;
步骤S2中判断电压暂降的起始点和结束点的具体依据为:
a)有效值电压降至0.9p.u.时的一刻为电压暂降的起始点;
b)有效值电压升至0.9p.u.时的一刻为电压暂降的结束点。
其中标幺值(p.u.)是电力系统分析和工程计算中常用的数值标记方法,表示各物理量及参数的相对值;
步骤S3中所述阈值为0.7p.u.;
步骤S4的小波变换为db6小波6层变换,得到细节分量D1、D2、D3、D4、D5、D6,其频率大小为D6<D5<D4<D3<D2<D1;
步骤S4计算各频段的信号能量的具体计算方法为:
定义单一尺度的小波能量为该尺度下重构信号的平方和:
式中,Ej为Dj细节分量的信号能量,Cj(k)为小波系数;
则信号总能量为
各频带能量所占比重得到向量组P=[P1 P2 P3 P4 P5 P6];
步骤S5中若有高频段的信号能量,则为雷击引起的故障;若没有高频段的信号能量,则为普通短路故障的具体步骤为:
若P1>0.1且P2>0.1,则为雷击引起的故障;若P1≤0.1或P2≤0.1,则为普通短路故障;
步骤S3中其它扰动原因为变压器激磁。
在具体实施过程中,如图2,雷电直击输电线路后,可能会造成电压暂降问题。雷电行波在系统中传播,系统中各节点电压由于行波传播与折反射会上升波动,在雷击造成系统故障后,工频短路电流必将引起电压暂降在系统中传播,因此节点电压在上升后再下降。而线路发生普通短路故障时,母线电压波形不会上升,而是直接下降发生电压暂降。雷击是造成故障的原因,不影响短路电流的大小,因此两种原因造成的同等故障,故障点位置、短路电流大小、类型都一致,因此电压暂降幅值一致;
如图3和图4,雷击引起的电压暂降波形经小波变换后得到的D1细节分量只在暂降发生时刻出现了高频分量,由变压器激磁引起的电压暂降波形经小波变换后得到的D1细节分量图形,在整个电压暂降期间都存在着高频分量,因此判断暂降幅值和高频分量出现的时刻可以区分出故障和其他扰动原因;
如图5、图6和表一可以得出雷击和普通故障导致的,其母线电压的能量分布存在很大的差异。雷击引起的电压暂降的能量在低频带(P6)和高频带(P1、P2)都存在,而普通故障的能量集中在低频带(P6)。其中雷击时P1、P2的最小值都在0.1以上,而普通故障时P1、P2的最大值仍为0,基于这种差异,可以分辨雷击和普通故障引起的电压暂降。
表一母线电压能量占比
P1 | P2 | P3 | P4 | P5 | P6 | ||
雷击 | 平均值 | 0.27 | 0.25 | 0.15 | 0.09 | 0.02 | 0.18 |
最小值 | 0.11 | 0.10 | 0.02 | 0.01 | 0.01 | 0.16 | |
普通故障 | 平均值 | 0 | 0 | 0 | 0 | 0.02 | 0.97 |
最大值 | 0 | 0 | 0 | 0 | 0.35 | 0.99 |
相同或相似的标号对应相同或相似的部件;
附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制;
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (7)
1.一种识别雷击引起的电压暂降的方法,其特征在于,包括以下步骤:
S1:从变电站的电能质量监测装置获取母线电压信号;
S2:计算有效值电压,判断电压暂降的起始点和结束点;
S3:计算电压暂降幅值Vsag,若Vsag大于阈值且电压暂降结束点有高频分量,则为其它扰动原因,进入步骤S6;若Vsag小于阈值或电压暂降结束点没有高频分量,则进入步骤S4;
S4:对电压暂降段的母线电压信号进行小波变换得到各频段的细节分量并计算各频段的信号能量;
S5:若有高频段的信号能量,则为雷击引起的故障;若没有高频段的信号能量,则为普通短路故障;
S6:判断结束。
2.根据权利要求1所述的识别雷击引起的电压暂降的方法,其特征在于,步骤S2中判断电压暂降的起始点和结束点的具体依据为:
a)有效值电压降至0.9p.u.时的一刻为电压暂降的起始点;
b)有效值电压升至0.9p.u.时的一刻为电压暂降的结束点。
3.根据权利要求1所述的一种识别雷击引起的电压暂降的方法,其特征在于,步骤S3中所述阈值为0.7p.u.。
4.根据权利要求1所述的识别雷击引起的电压暂降的方法,其特征在于,步骤S4的小波变换为db6小波6层变换,得到细节分量D1、D2、D3、D4、D5、D6,其频率大小为D6<D5<D4<D3<D2<D1。
5.根据权利要求4所述的识别雷击引起的电压暂降的方法,其特征在于,步骤S4计算各频段的信号能量的具体计算方法为:
定义单一尺度的小波能量为该尺度下重构信号的平方和:
式中,Ej为Dj细节分量的信号能量,Cj(k)为小波系数;
则信号总能量为
各频带能量所占比重得到向量组P=[P1 P2 P3 P4 P5 P6]。
6.根据权利要求5所述的识别雷击引起的电压暂降的方法,其特征在于,步骤S5中若有高频段的信号能量,则为雷击引起的故障;若没有高频段的信号能量,则为普通短路故障的具体步骤为:
若P1>0.1且P2>0.1,则为雷击引起的故障;若P1≤0.1或P2≤0.1,则为普通短路故障。
7.根据权利要求1所述的识别雷击引起的电压暂降的方法,其特征在于,步骤S3中其它扰动原因为变压器激磁。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811296162.9A CN109270333B (zh) | 2018-11-01 | 2018-11-01 | 一种识别雷击引起的电压暂降的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811296162.9A CN109270333B (zh) | 2018-11-01 | 2018-11-01 | 一种识别雷击引起的电压暂降的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109270333A true CN109270333A (zh) | 2019-01-25 |
CN109270333B CN109270333B (zh) | 2019-10-22 |
Family
ID=65191040
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811296162.9A Active CN109270333B (zh) | 2018-11-01 | 2018-11-01 | 一种识别雷击引起的电压暂降的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109270333B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111007357A (zh) * | 2019-12-18 | 2020-04-14 | 广东电网有限责任公司 | 一种考虑负荷电压暂降耐受能力曲线的暂降域识别方法 |
CN111639554A (zh) * | 2020-05-15 | 2020-09-08 | 神华包神铁路集团有限责任公司 | 牵引供电系统雷击识别方法、装置和设备 |
CN113804979A (zh) * | 2021-09-03 | 2021-12-17 | 青岛鼎信通讯股份有限公司 | 一种基于馈线终端的雷击波形识别方法 |
CN113946959A (zh) * | 2021-10-18 | 2022-01-18 | 国网河南省电力公司电力科学研究院 | 电压暂降数据段提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102135560A (zh) * | 2011-02-23 | 2011-07-27 | 山东大学 | 雷电波侵入变电站的扰动识别方法 |
CN103424600A (zh) * | 2013-08-20 | 2013-12-04 | 昆明理工大学 | 一种基于Hilbert-Huang变换和小波包能量谱的电压暂降源识别方法 |
CN103529332A (zh) * | 2013-10-28 | 2014-01-22 | 昆明理工大学 | 一种基于电压相关度和小波变换暂态能量分布特性的特高压直流输电线路雷击干扰识别方法 |
US20150094975A1 (en) * | 2013-10-01 | 2015-04-02 | King Fahd University Of Petroleum And Minerals | Wavelet transform system and method for voltage events detection and classification |
CN106443335A (zh) * | 2016-09-20 | 2017-02-22 | 广西电网有限责任公司钦州供电局 | 雷击故障识别方法及系统 |
-
2018
- 2018-11-01 CN CN201811296162.9A patent/CN109270333B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102135560A (zh) * | 2011-02-23 | 2011-07-27 | 山东大学 | 雷电波侵入变电站的扰动识别方法 |
CN103424600A (zh) * | 2013-08-20 | 2013-12-04 | 昆明理工大学 | 一种基于Hilbert-Huang变换和小波包能量谱的电压暂降源识别方法 |
US20150094975A1 (en) * | 2013-10-01 | 2015-04-02 | King Fahd University Of Petroleum And Minerals | Wavelet transform system and method for voltage events detection and classification |
CN103529332A (zh) * | 2013-10-28 | 2014-01-22 | 昆明理工大学 | 一种基于电压相关度和小波变换暂态能量分布特性的特高压直流输电线路雷击干扰识别方法 |
CN106443335A (zh) * | 2016-09-20 | 2017-02-22 | 广西电网有限责任公司钦州供电局 | 雷击故障识别方法及系统 |
Non-Patent Citations (5)
Title |
---|
GANG ZHENG等: "CLASSIFICATION OF VOLTAGE SAG BASED ON WAVELET TRANSFORM AND WAVELET NETWORK", 《PROCEEDINGS OF THE THIRD INTERNATIONAL CONFERENCE ON MACHINE LEARNING AND CYBERNETICS》 * |
M. MANJULA等: "Wavelet Transform for Classification of Voltage Sag Causes using Probabilistic Neural Network", 《INTERNATIONAL JOURNAL OF ELECTRICAL ENGINEERING》 * |
何梓田: "深圳电网电压暂降分析", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
无: "讲解!江苏电网电压暂降在线监测系统", 《北极星电力网》 * |
王克星等: "基于小波变换的配电网电压暂降的干扰源辨识", 《中国电机工程学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111007357A (zh) * | 2019-12-18 | 2020-04-14 | 广东电网有限责任公司 | 一种考虑负荷电压暂降耐受能力曲线的暂降域识别方法 |
CN111007357B (zh) * | 2019-12-18 | 2021-12-24 | 广东电网有限责任公司 | 一种考虑负荷电压暂降耐受能力曲线的暂降域识别方法 |
CN111639554A (zh) * | 2020-05-15 | 2020-09-08 | 神华包神铁路集团有限责任公司 | 牵引供电系统雷击识别方法、装置和设备 |
CN111639554B (zh) * | 2020-05-15 | 2023-04-21 | 神华包神铁路集团有限责任公司 | 牵引供电系统雷击识别方法、装置和设备 |
CN113804979A (zh) * | 2021-09-03 | 2021-12-17 | 青岛鼎信通讯股份有限公司 | 一种基于馈线终端的雷击波形识别方法 |
CN113804979B (zh) * | 2021-09-03 | 2023-06-16 | 青岛鼎信通讯股份有限公司 | 一种基于馈线终端的雷击波形识别方法 |
CN113946959A (zh) * | 2021-10-18 | 2022-01-18 | 国网河南省电力公司电力科学研究院 | 电压暂降数据段提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109270333B (zh) | 2019-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109270333B (zh) | 一种识别雷击引起的电压暂降的方法 | |
Wang et al. | Location of single phase to ground faults in distribution networks based on synchronous transients energy analysis | |
Li et al. | A new fault detection and fault location method for multi-terminal high voltage direct current of offshore wind farm | |
Luo et al. | A directional protection scheme for HVDC transmission lines based on reactive energy | |
Nanayakkara et al. | Traveling-wave-based line fault location in star-connected multiterminal HVDC systems | |
CN106093700B (zh) | 一种基于电压行波原理的故障录波装置及测距方法 | |
Galvez et al. | Fault location in active distribution networks containing distributed energy resources (DERs) | |
CN107482593B (zh) | 一种基于能量频段的柔性直流线路电流差动保护方法 | |
CN104237731B (zh) | 基于eemd与能量法的谐振接地配电网单相接地故障选线方法 | |
CN204422697U (zh) | 配电网故障定位装置 | |
CN109478778A (zh) | 用于检测三相配电网络中的故障的方法和装置 | |
CN102928704A (zh) | 一种变电站接地网腐蚀故障点智能诊断方法 | |
CN108020754A (zh) | 基于波形重构的单端行波故障测距方法 | |
WO2019097312A1 (en) | Parameter free traveling wave based fault location for power transmission lines | |
CN104316833A (zh) | 一种中压配电网故障行波选线方法 | |
CN109490706A (zh) | 一种多分支输电线路故障定位方法 | |
Hossam-Eldin et al. | Combined traveling wave and fuzzy logic based fault location in multi-terminal HVDC systems | |
CN109188193A (zh) | 一种基于特征频带收敛性闵氏距离的配电网故障选线方法 | |
CN105071447A (zh) | 一种交直流混联多馈入直流受端电网强度获得方法 | |
CN102540015B (zh) | 小电流接地系统单相接地测距方法和测距装置 | |
CN103424627A (zh) | 双端测量平行电网线路零序阻抗的方法 | |
CN109387713A (zh) | 一种分布式并网孤岛检测的混合方法 | |
Banafer et al. | Traveling wave-based primary protection and fault localization scheme for MTDC grid considering IEC 61869-9 measurement standard | |
CN103235198B (zh) | 一种同塔双回输电线路的导线相序排列的确定方法 | |
CN115136439A (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 |