CN115078909A - Non-ideal pulse leading edge detection method and device based on energy inflection point - Google Patents

Non-ideal pulse leading edge detection method and device based on energy inflection point Download PDF

Info

Publication number
CN115078909A
CN115078909A CN202210694580.3A CN202210694580A CN115078909A CN 115078909 A CN115078909 A CN 115078909A CN 202210694580 A CN202210694580 A CN 202210694580A CN 115078909 A CN115078909 A CN 115078909A
Authority
CN
China
Prior art keywords
inflection point
leading edge
energy
window
ideal pulse
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
Application number
CN202210694580.3A
Other languages
Chinese (zh)
Other versions
CN115078909B (en
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 Jiangsu Electric Power Co ltd Suqian Power Supply Branch
Suqian Electric Power Design Institute Co ltd
Original Assignee
State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch
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 Jiangsu Electric Power Co ltd Suqian Power Supply Branch filed Critical State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch
Priority to CN202210694580.3A priority Critical patent/CN115078909B/en
Publication of CN115078909A publication Critical patent/CN115078909A/en
Application granted granted Critical
Publication of CN115078909B publication Critical patent/CN115078909B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/083Locating faults in cables, transmission lines, or networks according to type of conductors in cables, e.g. underground
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/11Locating faults in cables, transmission lines, or networks using pulse reflection methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations

Abstract

The invention discloses a non-ideal pulse leading edge detection method and device based on an energy inflection point. For a waveform with a non-ideal pulse leading edge, the method regards the inflection point of the accumulated energy of the waveform as the starting moment of the pulse leading edge, and achieves the purpose of non-ideal pulse leading edge detection through the detection of the inflection point of the accumulated energy. The invention realizes the nonideal pulse leading edge detection method based on the energy inflection point and designs a matching device, and the method has the advantages of strong anti-noise capability and high accuracy, and effectively reduces the nonideal pulse leading edge detection error.

Description

Non-ideal pulse leading edge detection method and device based on energy inflection point
Technical Field
The invention relates to the field of edge detection algorithms, in particular to a non-ideal pulse leading edge detection method and device based on an energy inflection point.
Background
The cable is widely applied to power distribution systems of power grids and power supply systems in industries such as railways, aviation, manufacturing and the like. In order to avoid the influence of external force such as lightning stroke, wind damage and the like, the cable is usually buried in the system or underground, so that the fault of the cable cannot be identified and positioned by routing inspection like an overhead line, but a special method is needed.
The traveling wave method is a widely used cable fault location method, and detects and locates faults by using reflected traveling waves at discontinuous points of a cable, and the arrival time of reflected waves, namely the leading edge of the reflected waves, needs to be accurately measured. However, due to the dispersion effect of the cable itself, the bandwidth limitation of the receiving device, the signal interference and transmission on the propagation path, etc., the reflected wave may be distorted during the propagation process, and the pulse front edge no longer satisfies the characteristics of the ideal step signal, resulting in a deviation of the cable fault location algorithm based on the traveling wave arrival time. Therefore, a new method for detecting the leading edge of the non-ideal pulse with strong noise immunity and high accuracy is needed to reduce the error of the detection of the leading edge of the non-ideal pulse.
Disclosure of Invention
In order to solve the above mentioned drawbacks in the background art, the present invention provides a method and a device for detecting a leading edge of a non-ideal pulse based on an energy inflection point, wherein the method has the advantages of strong anti-noise capability and high accuracy, and effectively reduces the error of the detection of the leading edge of the non-ideal pulse.
The invention is realized by adopting the following technical scheme:
a non-ideal pulse front edge detection method based on energy inflection point regards the inflection point of the accumulated energy of the waveform as the starting time of the pulse front edge of a waveform with a non-ideal pulse front edge, and the purpose of non-ideal pulse front edge detection is achieved through the detection of the inflection point of the accumulated energy.
The method specifically comprises the following steps:
s1, sampling to obtain an original waveform;
s2, detecting the inflection point of the waveform accumulated energy;
and S3, obtaining the non-ideal pulse front edge starting moment according to the accumulated energy inflection point.
In step S2, the detection of the inflection point of the accumulated energy is specifically to perform inflection point detection on the accumulated energy by using a sliding window algorithm, and specifically includes:
s21, setting two adjacent windows to slide along the time direction of the waveform;
s22, calculating the difference between the first window and the second window;
and S23, detecting inflection points according to the differences.
The difference between the first window and the second window is specifically defined as:
d(f a..t ,f t..b )=|E a..t -E t..b | (1)
where d denotes the difference between the first and second windows, the subscript "a.. t" denotes "from time a to time t" and a and t denote the left and right edges of a time window, respectively, f denotes the signal itself, and E denotes the accumulated energy.
Preferably, the method for calculating the accumulated energy E specifically includes:
Figure BDA0003697922290000021
in practical applications, the integrated term can be replaced according to the effect and the calculation amount, including but not limited to f (tau) 2 And f (τ) | and the like.
In step S23, the definition of the inflection point is specifically: the position where the difference between the first window and the second window is the largest is considered as an inflection point.
In step S23, the step of detecting an inflection point according to the difference includes:
given a sliding time window width w, the knee detection function z (t) defining the sliding window is:
z(t)=d(f t-w..t ,f t..t+w ) (3)
the time at which the inflection point is located, i.e. the starting time of the leading edge of the non-ideal pulse, is:
Figure BDA0003697922290000031
in the formula, t dect Is the starting moment of the detected non-ideal pulse front edge.
The invention has the beneficial effects that:
the invention realizes the accurate detection of the leading edge of the non-ideal pulse based on the energy inflection point, improves the anti-noise capability of the detection of the leading edge of the non-ideal pulse, and effectively reduces the error of the detection of the leading edge of the non-ideal pulse.
Drawings
The invention will be further described with reference to the accompanying drawings.
FIG. 1 is a general flow diagram of the energy inflection point-based nonideal pulse front edge detection method and apparatus of the present invention;
FIG. 2 is a block diagram of a hardware device for the method and apparatus for energy inflection point based non-ideal leading edge detection of pulses of the present invention;
FIG. 3 is a diagram illustrating the effect of the non-ideal pulse front edge detection method and apparatus based on energy inflection point.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, a non-ideal pulse leading edge detection method based on energy inflection point includes the following steps:
s1, sampling to obtain an original waveform;
s2, setting a first sliding time window and a second sliding time window, and calculating the difference;
s3, solving the maximum difference position to obtain an accumulated energy inflection point;
and S4, obtaining the starting time of the non-ideal pulse front edge.
As shown in fig. 2, the non-ideal pulse leading edge detection device based on the energy inflection point comprises a voltage divider, a signal acquisition device, a signal transmission device and a signal processing device.
In the embodiment, in order to realize the attenuation of the high-voltage signal and meet the withstand voltage range of the signal acquisition equipment, a high-voltage broadband resistance-capacitance voltage divider is used for dividing the voltage of the original signal; collecting signals by using a high-speed data acquisition card; the signal transmission equipment transmits signals by using a wireless transmission module to realize remote monitoring and uploading of data; and signal processing is carried out on an STM32 board card, and the leading edge of the non-ideal pulse is detected, so that the detection of the non-ideal leading edge of the ns-level high-voltage pulse is realized.
The high-speed data acquisition card can adopt a commercially available high-speed data acquisition card with the sampling frequency of hundred megaHz grade, and uses LabView or MATLAB to drive a board card so as to effectively detect the nonideal pulse leading edge at ns grade.
The signal processing device adopts an MCU (micro controller unit) to reduce the consumption of computing resources.
A nonideal pulse leading edge detection method based on an energy inflection point specifically comprises the following steps: for a waveform with a non-ideal pulse leading edge, the inflection point of the accumulated energy of the waveform is regarded as the starting moment of the pulse leading edge, and the detection of the inflection point of the accumulated energy is used for achieving the purpose of non-ideal pulse leading edge detection.
In this embodiment, the detection of the inflection point of the accumulated energy specifically adopts a sliding window algorithm to perform inflection point detection on the accumulated energy:
two adjacent windows are arranged to slide along the time direction of the waveform, the difference between the first window and the second window is calculated, and an inflection point is detected according to the difference.
The difference between the first window and the second window is specifically defined as:
d(f a..t ,f t..b )=|E a..t -E t..b | (1)
where d denotes the difference between the first and second windows, the subscript "a.. t" denotes "from time a to time t" and a and t denote the left and right edges of a time window, respectively, f denotes the signal itself, and E denotes the accumulated energy.
In this embodiment, the method for calculating the accumulated energy E specifically includes:
Figure BDA0003697922290000051
in practical applications, the integrated term can be replaced according to the effect and the calculation amount, including but not limited to f (tau) 2 And f (τ) | and the like.
The definition of the inflection point is specifically as follows: the position where the difference between the first window and the second window is the largest is considered as an inflection point.
Detecting inflection points according to the difference, specifically:
given a sliding time window width w, which in this embodiment is chosen to be 100 Δ t, where Δ t is the sampling interval, the knee detection function z (t) defining the sliding window is:
z(t)=d(f t-w..t ,f t..t+w ) (3)
the moment of inflection point, i.e. the start moment t of the leading edge of the non-ideal pulse dect Comprises the following steps:
Figure BDA0003697922290000052
the method is used for detecting a waveform with a non-ideal pulse leading edge, the rising rate of the non-ideal pulse leading edge is 25ns, the signal-to-noise ratio of the waveform is 20dB, and the waveform is compared with a normalized threshold value comparison method and a Chichi information criterion method, wherein the normalized threshold value is selected to be 25%, and the implementation effect graph shown in figure 3 is obtained.
The above examples are merely illustrative of the specific procedures of the method of the present invention, and the scope of the claims of the present invention should not be limited thereto, and any modifications made based on the technical ideas of the present invention are within the scope of the present invention.

Claims (9)

1. A nonideal pulse leading edge detection method based on energy inflection points is characterized in that: for a waveform with a non-ideal pulse leading edge, the inflection point of the accumulated energy of the waveform is regarded as the starting moment of the pulse leading edge, and the detection of the inflection point of the accumulated energy is used for achieving the purpose of non-ideal pulse leading edge detection.
2. The method for detecting the non-ideal leading edge of the pulse based on the energy inflection point as claimed in claim 1, which comprises the following steps:
s1, sampling to obtain an original waveform;
s2, detecting the inflection point of the waveform accumulated energy;
and S3, obtaining the non-ideal pulse front edge starting moment according to the accumulated energy inflection point.
3. The method for detecting non-ideal pulse front edge based on energy inflection point as claimed in claim 2, wherein in said step S2, the detection of accumulated energy inflection point is performed by using a sliding window algorithm, specifically:
s21, setting two adjacent windows to slide along the time direction of the waveform;
s22, calculating the difference between the first window and the second window;
and S23, detecting an inflection point according to the difference.
4. The method according to claim 3, wherein the difference between the first window and the second window is defined as:
d(f a..t ,f t..b )=|E a..t -E t..b | (1)
where d denotes the difference between the first and second windows, the subscript "a.. t" denotes "from time a to time t" and a and t denote the left and right edges of a time window, respectively, f denotes the signal itself, and E denotes the accumulated energy.
5. The method for detecting the non-ideal pulse front edge based on the energy inflection point as claimed in claim 4, wherein the method for calculating the accumulated energy E specifically comprises:
Figure FDA0003697922280000011
6. the method for detecting non-ideal pulse front edge based on energy inflection point as claimed in claim 3, wherein in step S23, the inflection point is defined as: the position where the difference between the first window and the second window is the largest is considered as an inflection point.
7. The method for detecting a non-ideal pulse front edge based on energy inflection point as claimed in claim 6, wherein the step of detecting the inflection point according to the difference in step S23 includes:
given a sliding time window width w, the knee detection function z (t) defining the sliding window is:
z(t)=d(f t-w..t ,f t..t+w ) (3)
the time at which the inflection point is located, i.e. the starting time of the leading edge of the non-ideal pulse, is:
Figure FDA0003697922280000021
in the formula, t dect Is the starting moment of the detected non-ideal pulse front edge.
8. A detection device for the non-ideal pulse leading edge detection method based on energy inflection point as claimed in claim 1, which comprises a voltage divider, a signal acquisition device, a signal transmission device and a signal processing device;
the voltage divider is used for dividing the voltage of an original signal, so that the attenuation of a high-voltage signal is realized, and the withstand voltage range of the signal acquisition equipment is met;
the signal acquisition equipment is used for acquiring signals and adopts a high-speed data acquisition card;
the signal transmission equipment is used for realizing remote monitoring and uploading of data;
the signal processing device is used for detecting the leading edge of the non-ideal pulse, so that the non-ideal leading edge detection aiming at the ns-level high-voltage pulse is realized.
9. The apparatus according to claim 8, wherein the high speed data acquisition card is a high speed data acquisition card with a sampling frequency of the order of hundreds of megahertz.
CN202210694580.3A 2022-06-16 2022-06-16 Non-ideal pulse leading edge detection method and device based on energy inflection point Active CN115078909B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210694580.3A CN115078909B (en) 2022-06-16 2022-06-16 Non-ideal pulse leading edge detection method and device based on energy inflection point

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210694580.3A CN115078909B (en) 2022-06-16 2022-06-16 Non-ideal pulse leading edge detection method and device based on energy inflection point

Publications (2)

Publication Number Publication Date
CN115078909A true CN115078909A (en) 2022-09-20
CN115078909B CN115078909B (en) 2024-03-15

Family

ID=83252852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210694580.3A Active CN115078909B (en) 2022-06-16 2022-06-16 Non-ideal pulse leading edge detection method and device based on energy inflection point

Country Status (1)

Country Link
CN (1) CN115078909B (en)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6509740B1 (en) * 1999-04-09 2003-01-21 Jovial Test Equipment, Inc. Methods and apparatus for time domain reflectometry
JP2009216587A (en) * 2008-03-11 2009-09-24 Central Res Inst Of Electric Power Ind Transmission line failure location indicating method, apparatus and program
CN103323032A (en) * 2013-05-17 2013-09-25 厦门科华恒盛股份有限公司 Step signal detection method for smoothing of double sliding windows with dead zone unit
CN105223481A (en) * 2015-10-22 2016-01-06 西安交通大学 Based on the partial-discharge ultrahigh-frequency signal initial time defining method of difference energy function
CN105319487A (en) * 2015-10-30 2016-02-10 国家电网公司 Transformer station partial discharge signal detection and positioning system and method
CN105574284A (en) * 2015-12-29 2016-05-11 山东鲁能软件技术有限公司 Power equipment fault diagnosis method based on tendency characteristic point
CN107102244A (en) * 2017-05-25 2017-08-29 国网浙江省电力公司电力科学研究院 A kind of discharge source localization method of GIS ultrahigh frequency local discharge on-line monitoring device
CN107655533A (en) * 2017-10-31 2018-02-02 合肥工业大学 A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integration
CN109828187A (en) * 2018-11-19 2019-05-31 国网浙江省电力有限公司电力科学研究院 A kind of time difference extraction and analytical method positioning transformer equipment partial discharges fault
CN110658418A (en) * 2019-09-30 2020-01-07 山东信通电子股份有限公司 Cable fault detection method and device
CN110940933A (en) * 2019-12-02 2020-03-31 国网河北省电力有限公司电力科学研究院 Comprehensive calculation method for measuring rising edge starting time of steep pulse
KR102106307B1 (en) * 2019-11-11 2020-05-04 한국건설기술연구원 Energy data processing method and device
CN113639208A (en) * 2021-07-16 2021-11-12 汕头大学 Negative pressure wave signal inflection point positioning method and system based on extreme approximation
CN114386324A (en) * 2021-12-27 2022-04-22 上海电力大学 Ultra-short-term wind power segmented prediction method based on turning period identification

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6509740B1 (en) * 1999-04-09 2003-01-21 Jovial Test Equipment, Inc. Methods and apparatus for time domain reflectometry
JP2009216587A (en) * 2008-03-11 2009-09-24 Central Res Inst Of Electric Power Ind Transmission line failure location indicating method, apparatus and program
CN103323032A (en) * 2013-05-17 2013-09-25 厦门科华恒盛股份有限公司 Step signal detection method for smoothing of double sliding windows with dead zone unit
CN105223481A (en) * 2015-10-22 2016-01-06 西安交通大学 Based on the partial-discharge ultrahigh-frequency signal initial time defining method of difference energy function
CN105319487A (en) * 2015-10-30 2016-02-10 国家电网公司 Transformer station partial discharge signal detection and positioning system and method
CN105574284A (en) * 2015-12-29 2016-05-11 山东鲁能软件技术有限公司 Power equipment fault diagnosis method based on tendency characteristic point
CN107102244A (en) * 2017-05-25 2017-08-29 国网浙江省电力公司电力科学研究院 A kind of discharge source localization method of GIS ultrahigh frequency local discharge on-line monitoring device
CN107655533A (en) * 2017-10-31 2018-02-02 合肥工业大学 A kind of Ultrasonic Wave Flowmeter signal processing method and system based on backward energy integration
CN109828187A (en) * 2018-11-19 2019-05-31 国网浙江省电力有限公司电力科学研究院 A kind of time difference extraction and analytical method positioning transformer equipment partial discharges fault
CN110658418A (en) * 2019-09-30 2020-01-07 山东信通电子股份有限公司 Cable fault detection method and device
KR102106307B1 (en) * 2019-11-11 2020-05-04 한국건설기술연구원 Energy data processing method and device
CN110940933A (en) * 2019-12-02 2020-03-31 国网河北省电力有限公司电力科学研究院 Comprehensive calculation method for measuring rising edge starting time of steep pulse
CN113639208A (en) * 2021-07-16 2021-11-12 汕头大学 Negative pressure wave signal inflection point positioning method and system based on extreme approximation
CN114386324A (en) * 2021-12-27 2022-04-22 上海电力大学 Ultra-short-term wind power segmented prediction method based on turning period identification

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHARLES TRUONG: "Selective review of offline change point detection methods", 《SIGNAL PROCESSING》, no. 167, 19 September 2019 (2019-09-19), pages 1 - 20 *
李天辉: "最小能量法因子对RSO 法检测匝间短路故障的时间差 定位影响分析", 《高压电器》, vol. 55, no. 5, 16 May 2019 (2019-05-16), pages 54 - 61 *

Also Published As

Publication number Publication date
CN115078909B (en) 2024-03-15

Similar Documents

Publication Publication Date Title
CN101509949B (en) Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line
CN103176107B (en) High-voltage direct-current power transmission line hybrid fault ranging method
CN103901324A (en) Method for hybrid circuit combined distance measurement of power distribution network based on single-ended fault information
CN210923925U (en) C interface cable fault detection device of LEU
CN103809173A (en) Detection and tracking integration method for frame constant false-alarm target
CN111386468B (en) Travelling wave based fault localization using unsynchronized measurements for transmission lines
CN104406510A (en) Cable length measurement method based on EE (empirical envelope) method
CN115097253A (en) MMC-HVDC direct current transmission line fault distance measurement method and system
CN107179473B (en) A kind of power transmission line fault locating method
CN109324262A (en) A kind of fault positioning method for transmission line based on TT transformation and velocity of wave optimization
CN113406436B (en) Traveling wave fault location method and system for alternating-current and direct-current transmission line based on 5G communication
CN115047284A (en) Fault distance measuring method and system for high-voltage direct-current transmission line
CN115078909A (en) Non-ideal pulse leading edge detection method and device based on energy inflection point
CN108828398B (en) Single-core cable partial discharge positioning method based on structural similarity algorithm
CN111766470B (en) Fault positioning method and system for high-voltage direct-current transmission line and direct-current transmission line
CN109901016B (en) Double-end ranging method for eliminating wave velocity influence
CN102033229B (en) Tracking method of synthetic aperture radar signals
CN114675145B (en) High-frequency partial discharge double-end monitoring partial discharge source positioning method for high-voltage cable
CN111766477B (en) Method and device for rapidly detecting and identifying line faults of direct-current circuit-breaker-free power distribution network
CN116298657A (en) Multi-branch power distribution network fault positioning method under weak fault traveling wave
Wijaya et al. Review of transmission line fault location using travelling wave method
CN105467270B (en) Single Terminal Traveling Wave Fault Location back wave identification algorithm based on frequency spectrum similarity evaluation
Wang et al. Change point detection for piecewise envelope current signal based on wavelet transform
CN108572301B (en) Three-point current acquisition-based hybrid conductor fault positioning method and system
CN115023620A (en) Traveling wave based fault localization using parameter independence of unsynchronized measurements

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
TA01 Transfer of patent application right

Effective date of registration: 20240208

Address after: 223800 no.2481 Fazhan Avenue, Sucheng District, Suqian City, Jiangsu Province

Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. SUQIAN POWER SUPPLY BRANCH

Country or region after: China

Applicant after: SUQIAN ELECTRIC POWER DESIGN INSTITUTE Co.,Ltd.

Address before: 223800 no.2481 Fazhan Avenue, Sucheng District, Suqian City, Jiangsu Province

Applicant before: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. SUQIAN POWER SUPPLY BRANCH

Country or region before: China

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant