CN106505523B - A kind of excitation flow recognition method suitable for Traction networks transformer - Google Patents
A kind of excitation flow recognition method suitable for Traction networks transformer Download PDFInfo
- Publication number
- CN106505523B CN106505523B CN201611014755.2A CN201611014755A CN106505523B CN 106505523 B CN106505523 B CN 106505523B CN 201611014755 A CN201611014755 A CN 201611014755A CN 106505523 B CN106505523 B CN 106505523B
- Authority
- CN
- China
- Prior art keywords
- diff
- change rate
- phase current
- value
- current
- 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
Links
- 230000005284 excitation Effects 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000008859 change Effects 0.000 claims abstract description 41
- 238000005070 sampling Methods 0.000 claims abstract description 22
- 150000001875 compounds Chemical class 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000004804 winding Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 9
- 230000008901 benefit Effects 0.000 description 5
- 230000004069 differentiation Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005415 magnetization Effects 0.000 description 3
- 230000007257 malfunction Effects 0.000 description 3
- 238000000518 rheometry Methods 0.000 description 3
- 238000007716 flux method Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
- H02H7/00—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
- H02H7/04—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions for transformers
- H02H7/045—Differential protection of transformers
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Protection Of Transformers (AREA)
Abstract
The invention discloses a kind of excitation flow recognition methods suitable for Traction networks transformer, comprising: step 1: acquisition tractive transformer winding three-phase electric current, and carries out difference processing to each phase current;Step 2: calculating the three-phase current change rate in power frequency period T;Step 3: any to read T/2 cycle data window, and the difference of three-phase current change rate absolute value corresponding to the point-by-point sampling number for sequentially elapsing calculating three-phase current change rate half period adjacent thereto backward and and value;Step 4: utilizing N/2 difference accumulated value Δ i 'diff‑φ‑1, N/2 and value accumulated value Δ i 'diff‑φ‑2, foundationFurther criterion identification is shoved and internal fault.A kind of method using interval angle and asymmetry compound characteristics identification excitation surge current of the present invention, reduces the influence of DC component, and this method principle is simple, calculation amount is small, and recognition accuracy is higher, is easily applied in engineering practice.
Description
[technical field]
The invention belongs to tractive power supply system tractive transformer technical field of relay protection, and in particular to one kind is suitable for leading
Draw the excitation flow recognition method of net transformer.
[background technique]
Tractive transformer supplies entire traction as part most important in tractive power supply system, reliable, safe operation
The safe operation important in inhibiting of electric system.Excitation surge current is the principal element for causing transformer main protection differential protection malfunction
One of.How rapidly and accurately to identify transformer excitation flow and internal fault, helps to promote Traction networks tranformer protection work
The reliability of work.
In recent years, domestic and foreign scholars propose many methods to transformer identification excitation surge current problem, according to utilizing electrical quantity
It can be divided into: 1) current characteristic method;2) voltage characteristic method;3) Current Voltage compound characteristics method.
1, the excitation flow recognition method of current characteristic: mainly including the method for quadratic harmonics, interruption horn cupping, waveform symmetry method
Deng.The method of quadratic harmonics has the characteristics that principle is simply easily achieved, but is influenced by transient current characteristic and transformer magnetization characteristic
It is larger, and the adjusting of secondary harmonic brake threshold is difficult.It is interrupted horn cupping and the excitation surge current time difference is mainly generated according to tractive transformer
Streaming current has apparent interruption corner characteristics, and difference current is then without apparent interruption corner characteristics when failure.Waveform symmetry for
Threshold value chooses none foundation determined, is only chosen by engineering experience.
2, the excitation flow recognition method of voltage characteristic: mainly there are works, transformer equivalent circuit method including magnetic flux method, difference
Deng.Magnetic flux method can preferably reflect the transient process of transformer, but stronger to the characteristic and parameter dependence of transformer itself.
The voltage and current information that difference has works to take full advantage of transformer identifies the state of transformer, but the adjusting of threshold value
It is more difficult.Although transformer equivalent circuit does not need the magnetization curve of transformer, but to the winding parameter of transformer have compared with
High requirement.
3, the excitation flow recognition method of voltage and current compound characteristics: mainly including magnetic flux characteristic, power differential principle etc..
Since transformer remanent magnetism obtains difficult, to obtain in the case of excitation surge current ψ-idCurve will deviate from magnetization curve, so as to cause magnetic
The erroneous judgement of logical Characterization method, exists simultaneously subregion difficulty and the adjusting of restraint coefficient threshold value is complicated.Power differential method needs are avoided gushing
The charging process in the 1st period of transformer when stream causes to differentiate delay;Copper loss accurately calculates difficulty when excitation surge current, therefore utilizes
Threshold value adjusting is more difficult when principle identification is shoved.
In addition, some scholars combine the intellectual technologies such as fuzzy logic, wavelet transformation, neural network to carry out excitation surge current identification
Method, but intelligent algorithm identification shoves there are biggish uncertainty and realizes that algorithm is complicated, computationally intensive, has been difficult in a short time
Effect is used for Practical Project.
[summary of the invention]
For the deficiency of existing trailer system transformer excitation flow identification, the object of the present invention is to provide one kind to be suitable for
The excitation flow recognition method of Traction networks transformer solves asking for transformer differential protection malfunction existing in the prior art
Topic.
Object of the present invention is to what is be achieved through the following technical solutions:
A kind of excitation flow recognition method suitable for Traction networks transformer, comprising the following steps:
Step 1: acquisition transformer A, B, C three-phase current iφ(k), wherein φ=A, B, C, and it is poor to the progress of each phase current
Divide processing;
Step 2: calculating the three-phase current change rate i ' in a power frequency period Tdiff-φ(k);
Step 3: it is any to read T/2 cycle data window, and point-by-point sequentially elapse backward calculates three-phase current change rate i 'diff-φ
(k) and three-phase current change rate corresponding to the sampling number of adjacent integers times half periodCalculate separately change
Rate absolute value differencesAnd change rate absolute value and value
The sampling number that wherein k is sampling number, N is a power frequency calculating cycle, n=1,2 ...;
Step 4: according to the sampling number of three-phase current change rate half period adjacent thereto three-phase electricity rheology calculated
The difference accumulated value Δ i ' of ratediff-φ-1With with value accumulated value Δ i 'diff-φ-2, shove according to the identification of following inrush current distinguishing criterion
And internal fault,
Wherein, KresTo brake threshold;N is integer value;In one cycle, as long as A, B, C three-phase current have phase satisfaction
Above-mentioned inrush current distinguishing criterion, then be judged to shove, transformer differential protection locking;Otherwise it is judged to internal fault.
Step 1 specifically: to the three-phase current i of transformerA(k)、iB(k)、iC(k) it is acquired, calculates each phase differential electrical
Flow idiff-A(k)=iA(k)-iA(k-1)、idiff-B(k)=iB(k)-iB(k-1)、idiff-C(k)=iC(k)-iC(k-1), wherein
K=2,3,4 ....
Step 2 specifically: calculate each phase difference current idiff-A(k)、idiff-B(k)、idiff-C(k) it is calculated in a power frequency
Current changing rate i ' in cycle Tdiff-φ(k):
Current changing rate:
Wherein, k is sampling number, TsFor the sampling period.
Step 3 specifically: calculate step 2 resulting three-phase current change rate half week adjacent thereto in some sampled point
It is poor that three-phase current change rate absolute value corresponding to the sampling number of phase is made, and N/2 calculated value is added up and obtains Δ i 'diff-φ-1,
That is:
Then to corresponding to the sampling number of three-phase current change rate half period adjacent thereto corresponding to this sampled point
Three-phase current change rate absolute value summation, by N/2 calculated value add up obtain Δ i 'diff-φ-2, i.e.,
KresValue is that 0.2, n takes 1~3.
Compared with the existing technology, the invention has the following advantages that
A kind of excitation flow recognition method using interval angle and asymmetry compound characteristics of the present invention, takes full advantage of and leads
Wave character when drawing transformer excitation flow handles the excitation surge current being exaggerated under weak current operating condition and inside using differential calculation
The interval angle and pinnacle waveform feature difference of fault current are expected to improve the inrush current distinguishing of trailer system transformer differential protection
Energy.Excitation flow recognition method of the invention, the interval angle to shove when in view of transformer excitation flow and internal fault and pinnacle
Wave characteristic difference is realized and is shoved and internal fault accurately identifies.Fault current waveform in internal fault, after differential process
In significant periodic sinusoidal wave characteristic in any half of power frequency calculating cycle, at least exist a curent change extreme point and
Electric current has symmetric characteristics;When transformer during no-load closing generates excitation surge current, due to the interval angle and pinnacle wave characteristic to shove,
Electric current after differential is only in peaked wave region there are curent change extreme point and electric current does not have symmetric characteristics.Therefore, sharp
The reliable differentiation shoved with internal fault may be implemented with this curent change feature, reduce the influence of DC component, and
This method principle is simple, and calculation amount is small, and recognition accuracy is higher, is easily applied in engineering practice.It is identified with other excitation surge currents
Principle is compared, and the present invention has following remarkable advantage:
(1) it using the compound characteristics of the interval angle of transformer excitation flow electric current and pinnacle, improves and utilizes magnitude of current spy
The applicability of sign identification excitation surge current;
(2) whether the electric current after being gone out using differential still has whether interval angle and curent change extreme point have symmetry,
It can rapidly and accurately identify excitation surge current.
[Detailed description of the invention]
Fig. 1 is a kind of process using interval angle and the excitation flow recognition method of asymmetry compound characteristics of the present invention
Figure.
[specific embodiment]
With reference to the accompanying drawing, a specific embodiment of the invention is described in detail, but the present invention is not limited to the implementations
Example.In order to make the public have thorough understanding to the present invention, is preferably applied in following present invention and concrete details is described in detail in example.
The principle of the present invention are as follows: static exciter circuit is equivalent to the fault branch of power transformer interior fault, works as transformer
When core sataration, very big excitation surge current can be generated and flowed into differential relay, differential protection misoperation is caused.Due to excitation
Electric current is very big, if escaping its influence with action current, differential protection sensitivity in power transformer interior fault will be reduced.Mirror
In this, the present invention identifies excitation surge current with internal fault current wave character difference using excitation surge current and then prevents excitation from gushing
Stream causes the malfunction of differential protection.Since in internal fault, the fault current after differential process is in any half of power frequency
In calculating cycle, at least there is a curent change extreme point and electric current has symmetric characteristics;When transformer during no-load closing, warp
Electric current after differential process be there are interruption corner characteristics, and only in peaked wave region there are curent change extreme point and electric current not
With symmetric characteristics.Therefore, the reliable differentiation shoved with internal fault may be implemented using this curent change feature.
A kind of method using interval angle and asymmetry compound characteristics identification excitation surge current of the present invention, specific process is such as
Shown in Fig. 1, follow the steps below to implement:
Step 1: acquisition Three-Phase Transformer electric current iφ(k), wherein φ=A, B, C respectively represent Three-Phase Transformer circuit
A wherein phase, i.e. iA(k)、iB(k)、iC(k), each phase difference current i is calculateddiff-A(k)=iA(k)-iA(k-1)、idiff-B(k)-
iB(k)-iB(k-1)、idiff-C(k)=iC(k)-iC(k-1), wherein k=2,3,4 ...;
Since transformer winding is in inductive, excitation surge current and internal fault current contain part aperiodic component, benefit
Reduce the influence of aperiodic component with difference processing.
Step 2: calculating the three-phase current change rate i ' in power frequency period Tdiff-φ(k), specifically:
Calculate each phase difference current idiff-A(k)、idiff-B(k)、idiff-C(k) electric current in a power frequency calculating cycle T
Change rate i 'diff-φ(k) i.e.
Current changing rate:
Wherein, k is sampling number, TsFor the sampling period;
Step 3: it is shoved and internal fault using criterion identification:
Step 2 is calculated to the sampling number of resulting three-phase current change rate half period adjacent thereto in some sampled point
It is poor that corresponding three-phase current change rate absolute value is made, and N/2 calculated value is added up and obtains Δ i 'diff-φ-1I.e.
Then to corresponding to the sampling number of three-phase current change rate half period adjacent thereto corresponding to this sampled point
Three-phase current change rate absolute value summation, by N/2 calculated value add up obtain Δ i 'diff-φ-2I.e.
Step 4 specifically:
According to the sampling number three-phase electricity rheology calculated of the three-phase current change rate half period adjacent thereto of step 3
The difference accumulated value Δ i ' of ratediff-φ-1With with value accumulated value Δ i 'diff-φ-2, shove according to the identification of following inrush current distinguishing criterion
And internal fault:
Wherein, KresTo brake threshold, considering the measurement error of current transformer and calculating error, KresValue is 0.2;n
Round numbers can use 1~3 according to the reliability value of differentiation.
In one cycle, only there are electric current changes in peaked wave region for the electric current due to excitation surge current after differential process
Change extreme point and electric current does not have symmetric characteristics, and at least there is an electricity in fault current of the internal fault after differential process
Rheology extreme point and extreme point has symmetry feature.As long as therefore three-phase current has a phase to meet above-mentioned inrush current distinguishing and sentences
According to, then it is judged to shove, transformer differential protection locking;Otherwise it is judged to internal fault.
The method of the present invention proposes a kind of utilization using tractive transformer excitation surge current and internal fault current feature difference
Interval angle and the excitation surge current of asymmetry feature and internal fault recognition methods.This method generates excitation using tractive transformer
Difference current has apparent interruption corner characteristics and curent change extreme point not to have symmetry feature when shoving, and differential when failure
Without being significantly interrupted corner characteristics and at least there are the characteristics that an extreme point and extreme point has symmetry differentiation in electric current
Excitation surge current and fault current angle characteristic, this method principle is simple, and calculation amount is small, has engineering application value.
More than, only presently preferred embodiments of the present invention is not limited only to practical range of the invention, all according to the invention patent
The equivalence changes and modification that the content of range is done all should be technology scope of the invention.
Claims (2)
1. a kind of excitation flow recognition method suitable for Traction networks transformer, which comprises the following steps:
Step 1: acquisition transformer A, B, C three-phase current iφ(k), wherein φ=A, B, C, and each phase current is carried out at difference
Reason;
Step 2: calculating the three-phase current change rate i ' in a power frequency period Tdiff-φ(k);
Step 3: it is any to read T/2 cycle data window, and point-by-point sequentially elapse backward calculates three-phase current change rate i 'diff-φ(k) and phase
Three-phase current change rate corresponding to the sampling number of adjacent integral multiple half periodIt is exhausted to calculate separately change rate
To value difference valueAnd change rate absolute value and value
The sampling number that wherein k is sampling number, N is a power frequency calculating cycle, n are positive integer;
Step 4: according to the sampling number of three-phase current change rate half period adjacent thereto three-phase current change rate calculated
Difference accumulated value Δ i 'diff-φ-1With with value accumulated value Δ i 'diff-φ-2, whether the electric current after being gone out using differential, which is still had, is interrupted
Whether angle and curent change extreme point have symmetry, identify excitation surge current;According to following inrush current distinguishing criterion identification shove and
Internal fault,
Wherein, KresTo brake threshold;N is integer value;In one cycle, as long as to have a phase to meet above-mentioned for A, B, C three-phase current
Inrush current distinguishing criterion is then judged to shove, transformer differential protection locking;Otherwise it is judged to internal fault;
Step 2 specifically: calculate each phase difference current idiff-A(k)、idiff-B(k)、idiff-C(k) in a power frequency calculating cycle T
Interior current changing rate i 'diff-φ(k):
Current changing rate:
Wherein, k is sampling number, TsFor the sampling period;
Step 3 specifically: step 2 is calculated into resulting three-phase current change rate half period adjacent thereto in some sampled point
It is poor that three-phase current change rate absolute value corresponding to sampling number is made, and N/2 calculated value is added up and obtains Δ i 'diff-φ-1, it may be assumed that
Then to three corresponding to the sampling number of three-phase current change rate half period adjacent thereto corresponding to this sampled point
The summation of phase current change rate absolute value, N/2 calculated value is added up and obtains Δ i 'diff-φ-2, i.e.,
Step 1 specifically: to the three-phase current i of transformerA(k)、iB(k)、iC(k) it is acquired, calculates each phase difference current
idiff-A(k)=iA(k)-iA(k-1)、idiff-B(k)=iB(k)-iB(k-1)、idiff-C(k)=iC(k)-iC(k-1), wherein k
=2,3,4 ....
2. a kind of excitation flow recognition method suitable for Traction networks transformer according to claim 1, which is characterized in that
KresValue is that 0.2, n takes 1~3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611014755.2A CN106505523B (en) | 2016-11-15 | 2016-11-15 | A kind of excitation flow recognition method suitable for Traction networks transformer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611014755.2A CN106505523B (en) | 2016-11-15 | 2016-11-15 | A kind of excitation flow recognition method suitable for Traction networks transformer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106505523A CN106505523A (en) | 2017-03-15 |
CN106505523B true CN106505523B (en) | 2019-05-03 |
Family
ID=58324923
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611014755.2A Active CN106505523B (en) | 2016-11-15 | 2016-11-15 | A kind of excitation flow recognition method suitable for Traction networks transformer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106505523B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108899871B (en) * | 2018-05-08 | 2019-11-19 | 广东电网有限责任公司 | A kind of recognition methods of excitation surge current, device and computer-readable medium |
CN110161297B (en) * | 2019-06-28 | 2021-05-14 | 沈阳工业大学 | Method for calculating winding current of power transformer under reclosing condition |
CN111103481A (en) * | 2019-12-11 | 2020-05-05 | 国家电网有限公司 | Method for identifying magnetizing inrush current of transformer |
CN111580035B (en) * | 2020-05-27 | 2022-12-09 | 国网江苏省电力有限公司 | Statistical identification method for magnetic saturation of voltage transformer |
CN116298580B (en) * | 2022-12-09 | 2024-04-02 | 国家电网有限公司 | Transformer inrush current identification method and system based on waveform deviation and interruption characteristics |
CN116706835B (en) * | 2022-12-29 | 2024-02-23 | 国家电网有限公司 | Method, device, medium and equipment for identifying inrush current of transformer |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10313531A (en) * | 1997-05-07 | 1998-11-24 | Toshiba Corp | Ratio differential relay |
CN101257208A (en) * | 2007-12-21 | 2008-09-03 | 清华大学 | Method for identifying transformer excitation surge current |
CN103675415B (en) * | 2012-08-31 | 2017-03-01 | 西门子公司 | Excitation surge current detection method, braking magnetizing inrush current method and its device |
CN103683198B (en) * | 2013-12-03 | 2017-01-04 | 昆明理工大学 | The excitation surge current method for quickly identifying of consecutive points distance in a kind of plane based on difference current adjacent order difference composition |
CN105653882A (en) * | 2016-03-09 | 2016-06-08 | 西安工程大学 | Method for identifying magnetizing inrush current by utilizing current waveform characteristics |
-
2016
- 2016-11-15 CN CN201611014755.2A patent/CN106505523B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN106505523A (en) | 2017-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106505523B (en) | A kind of excitation flow recognition method suitable for Traction networks transformer | |
Kang et al. | Compensation of the distortion in the secondary current caused by saturation and remanence in a CT | |
Ajaei et al. | Compensation of the current-transformer saturation effects for digital relays | |
Hajipour et al. | Current-transformer saturation compensation for transformer differential relays | |
Wiszniewski et al. | Correction of current transformer transient performance | |
CN103728538B (en) | A kind of ground fault line selecting method of small current neutral grounding system | |
CN103050941B (en) | Morphological gradient-based identification method for magnetizing inrush current of transformer | |
CN107390010B (en) | Method for rapidly detecting trailing current of current transformer | |
CN102353873A (en) | Line selection method realized by using single criterion for low-current ground system during single-phase ground fault | |
He et al. | A new method to identify inrush current based on error estimation | |
CN110579682A (en) | Fault recording data transient homologous comparison method and device | |
CN105653882A (en) | Method for identifying magnetizing inrush current by utilizing current waveform characteristics | |
CN104237731A (en) | Single-phase earth fault line selection method on basis of EEMD (ensemble empirical mode decomposition) and energy method for resonance earthed power distribution network | |
Bi et al. | Correlation analysis of waveforms in non-saturation zone based method to identify the magnetizing inrush in transformer | |
CN108490310A (en) | Method for small electric current grounding system of distribution network earth fault line selection | |
Li et al. | High-speed directional pilot protection for MVDC distribution systems | |
CN106569163A (en) | Method for detecting trailing amount of transformer | |
CN105866511B (en) | A method of the identification tail currents based on current waveform neighbouring sample point ratio | |
Chen et al. | Harmonic components based protection strategy for inverter-interfaced AC microgrid | |
CN106324328B (en) | A kind of transformer excitation flow recognition method based on morphology cascade erosion operation | |
CN109088402A (en) | Transmission line of electricity time domain Full wave shape protective device and its phase Sudden Changing Rate direction determination process | |
Lee et al. | Phasor estimation algorithm based on the least square technique during CT saturation | |
CN104977488A (en) | Differential current gradient angle approximate entropy-based transformer magnetizing inrush current identification method | |
Wang et al. | Protection of line faults in HVDC grids through convexity detection in backward traveling wave voltages | |
CN109245057B (en) | Time domain full-waveform protection device for power transmission line and interphase mutation direction judgment method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |