CN108957097A - Method for measuring resistive current fundamental wave of metal oxide arrester - Google Patents

Method for measuring resistive current fundamental wave of metal oxide arrester Download PDF

Info

Publication number
CN108957097A
CN108957097A CN201810704845.7A CN201810704845A CN108957097A CN 108957097 A CN108957097 A CN 108957097A CN 201810704845 A CN201810704845 A CN 201810704845A CN 108957097 A CN108957097 A CN 108957097A
Authority
CN
China
Prior art keywords
leakage current
fundamental wave
formula
moa
phase angle
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
CN201810704845.7A
Other languages
Chinese (zh)
Other versions
CN108957097B (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.)
Jiangsu University of Technology
Original Assignee
Jiangsu University of Technology
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 Jiangsu University of Technology filed Critical Jiangsu University of Technology
Priority to CN201810704845.7A priority Critical patent/CN108957097B/en
Publication of CN108957097A publication Critical patent/CN108957097A/en
Application granted granted Critical
Publication of CN108957097B publication Critical patent/CN108957097B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Electric Properties And Detecting Electric Faults (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)

Abstract

The invention relates to a method for measuring resistive current fundamental wave of a metal oxide arrester, wherein the parameters of the resistive leakage current fundamental wave are obtained by adopting a projection method, and the specific process is as follows: simultaneously sampling a power grid voltage and an MOA leakage current signal; using quasi-synchronous DFT harmonic phase angle linear correctionMethod for obtaining fundamental wave initial phase angle of power grid voltage by square methodMethod for obtaining fundamental wave amplitude I of MOA leakage current by applying quasi-synchronous DFT harmonic phase angle linear correction method1And initial phase angleObtaining the projection angle of the leakage current fundamental wave on the grid voltage fundamental wave according to the projection methodObtaining MOA resistive leakage current fundamental wave

Description

The measurement method of the resistive current first harmonics of metal oxide arrester
The application be application No. is: 201510257527.7, invention and created name is " a kind of MOA resistive current first harmonics Measurement method ", the applying date are as follows: the divisional application of the application for a patent for invention on May 19th, 2015.
Technical field
The present invention relates to the technical field of MOA testing current in resistance property, specifically a kind of high-precision MOA resistive current first harmonics Measurement method.
Background technique
Metal oxide arrester (hereinafter referred to as MOA) obtains in the power system because of its superior overvoltage protection characteristic It is widely applied, but the aging of MOA resistor disc and warp are heated and impact failure can cause failure, seriously may result in its explosion, Arrester breakdown also results in substation bus bar short circuit, influences system safety operation, it is therefore necessary to carry out to running MOA tight Lattice effectively detect and periodic preventative test.In the detection and test of Zinc-Oxide Arrester, letting out under ac operation voltage Leakage current measurement is an important project, and the size of the resistive current first harmonics ingredient in leakage current can accurate reflection oxygen Change the defects of zinc arrester dampness, the aging of valve block and built-in electrical insulation are damaged.Therefore, to the parameters such as resistive leakage current fundamental wave into The accurate measurement of row can relatively accurately differentiate MOA performance.
Frequency analysis technology is answered in various fields such as electric energy quality monitoring, electronic product production testing, electric appliances monitorings It is the important technical for carrying out power system monitor, quality inspection, monitoring of tools with extensive.Frequency analysis is most widely used at present Technology be discrete Fourier transform (DFT) and Fast Fourier Transform (FFT) (FFT).Quasi-synchronous sampling technique is mutually tied with DFT technique The frequency analysis technology of conjunction can be improved the precision of frequency analysis, formula are as follows:
In formula: k is the number (such as fundamental wave k=1,3 subharmonic k=3) for needing the harmonic wave obtained;Sin and cos are positive respectively String and cosine function;And akAnd bkThe respectively real and imaginary parts of k subharmonic;N is the number of iterations;W is determined by integration method, is adopted When with muiltiple-trapezoid integration method, W=nN;γiFor a weighting coefficient;For the sum of all weighting coefficients;f(i) For the ith sample value of analysis waveform;N is sampling number in the period.
In engineer application, frequency analysis always carry out the sampling of finite point be difficult to stricti jurise synchronize adopt Sample.In this way, when the plesiochronous DFT of application carries out frequency analysis, will exist the leakage of the long range as caused by truncation effect and The short range leakage as caused by fence effect, so that analysis result precision is not high or even insincere.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of current in resistance property of the higher metal oxide arrester of precision The measurement method of fundamental wave obtains high-precision frequency analysis to be efficiently modified the analytical error of plesiochronous DFT frequency analysis technology As a result, to improve the reliability of the MOA testing current in resistance property based on frequency analysis theory.
Realize that the technical solution of the object of the invention is to provide a kind of survey of the resistive current first harmonics of metal oxide arrester Amount method, comprising the following steps:
(1) the W+2 sample point data of synchronized sampling mains voltage signal V and MOA leakage current signal I at equal intervals: { fV (i),fI(i), i=0,1 ..., W+1 };
(2) plesiochronous DFT formula is applied since the sampled point i=0 of the mains voltage signal V:
W+1 data of analysis obtain the mains voltage signal V Fundamental informationWith
Plesiochronous DFT formula is applied from the sampled point i=1 of the mains voltage signal V:
W+1 data of analysis obtain the mains voltage signal V Fundamental informationWith
Using formula:Calculate the frequency drift μ of the mains voltage signal VV
Using formulaCalculate the fundamental wave initial phase angle of the mains voltage signal V;
Using formulaLinearly correct the fundamental wave initial phase angle of the mains voltage signal V.
(3) plesiochronous DFT formula is applied since the sampled point i=0 of the MOA leakage current signal I:
W+1 data of analysis obtain the MOA leakage current letter The fundamental information of number IWith
Plesiochronous DFT formula is applied from the sampled point i=1 of the MOA leakage current signal I:
W+1 data of analysis obtain the MOA leakage current letter The fundamental information of number IWith
Using formula:Calculate the frequency drift μ of the MOA leakage current signal I1
Using formulaCalculate the fundamental wave initial phase angle of the MOA leakage current signal I;
Using formulaLinearly correct the fundamental wave initial phase angle of the MOA leakage current signal I.
Using formulaCalculate the fundamental voltage amplitude of the MOA leakage current signal I;
(4) formula is appliedCalculate projected angle of the MOA leakage current fundamental wave on network voltage fundamental wave;
(5) formula is appliedCalculate MOA resistive leakage current fundamental wave.
The main reason for plesiochronous DFT frequency analysis can effectively inhibit long range and leak, spectrum leakage is signal Short range leakage caused by frequency drift, the present invention, which discloses one kind, can effectively inhibit the humorous phase angle of short range leakage linearly to correct Method, to obtain high-precision humorous phase angle information and dielectric dissipation factor.
N is the sampling number in an ideal period.The synchronized sampling at equal intervals is according to progress frequency analysis The cycle T and frequency f (such as power frequency component frequency f is 50Hz, period 20mS) of ideal signal, sample N point in one cycle, I.e. sample frequency is fs=Nf, and N >=64.
Described W+2 sample point data of sampling is accordingly selected according to selected integration method, according to multiple Change trapezoidal integration method, then W=nN;According to complexification rectangular integration method, then W=n (N-1);According to iterative Simpson product Divide method, then W=n (N-1)/2;Then according to sample frequency fs=Nf obtains sample point data sequence;N is the number of iterations, one As n >=3.
An iteration coefficient gammaiIt is determined by integration method, ideal period sampled point N and the number of iterations n, specific derivation process Referring to document [wear some problem [J] the electrical measurement and instrument in the application of quasi-synchro sampling in elder generation, 1988, (2): 2-7.].
For the sum of all weighting coefficients.
The drift μ of signal frequencyvAnd μ1It is according to sampling number N in neighbouring sample point fundamental wave phase angle difference and ideal period Fixed relationship and obtain, the drift of signal frequency can also be used for amendment fundamental wave and higher hamonic wave frequency f1With higher hamonic wave Frequency fk
The present invention has the effect of positive: (1) resistive current first harmonics measurement method of the invention, humorous for improving The quality of wave analysis obtains fundamental voltage amplitude, the phase angle of high-precision network voltage and MOA leakage current, and then obtains high-precision Resistive current first harmonics information, judge MOA performance provide foundation.
(2) method of the present invention fundamentally solves the problems, such as that the humorous phase angle analysis precision of plesiochronous DFT is low, and nothing Complicated inverting and amendment need to be carried out, algorithm is simple.
(3) relative to plesiochronous DFT, frequency analysis technology of the present invention only needs to increase a sampled point and just solves Plesiochronous DFT analytical error big problem, it is easy to accomplish.
(4) existing instrument and equipment is improved using the present invention, be technically feasible, and do not need to increase any hard Part expense can be such that analysis result can be improved to 10-8Grade.
(5) this method is similarly also applied for carrying out successive ignition rather than the frequency analysis process of an iteration, at this time only Needing an iteration to resolve into successive ignition realization can.As an iteration with successive ignition is substantially, only It is that when calculating, successive ignition carries out decoupled method, and an iteration is that the process of successive ignition is merged into iteration coefficient γi In once calculate complete, so the present disclosure applies equally to successive ignition processes.
Specific embodiment
A kind of measurement method of the resistive current first harmonics of metal oxide arrester of the invention, comprising the following steps:
(1) the W+2 sample point data of synchronized sampling mains voltage signal V and MOA leakage current signal I at equal intervals: { fV (i),fI(i), i=0,1 ..., W+1 };
Described W+2 sample point data of sampling is accordingly selected according to selected integration method, according to multiple Change trapezoidal integration method, then W=nN;According to complexification rectangular integration method, then W=n (N-1);According to iterative Simpson product Divide method, then W=n (N-1)/2;Then according to sample frequency fs=Nf obtains sample point data sequence;N is the number of iterations, one As n >=3.
(2) plesiochronous DFT formula is applied since the sampled point i=0 of the mains voltage signal V:
W+1 data of analysis obtain the mains voltage signal The fundamental information of VWith
Plesiochronous DFT formula is applied from the sampled point i=1 of the mains voltage signal V:
W+1 data of analysis obtain the mains voltage signal V Fundamental informationWith
Using formula:Calculate the frequency drift μ of the mains voltage signal Vv
Using formulaCalculate the fundamental wave initial phase angle of the mains voltage signal V;
Using formulaLinearly correct the fundamental wave initial phase angle of the mains voltage signal V.
Then, plesiochronous DFT formula is applied since the sampled point i=0 of the MOA leakage current signal I:
W+1 data of analysis obtain the MOA leakage current letter The fundamental information of number IWith
Plesiochronous DFT formula is applied from the sampled point i=1 of the MOA leakage current signal I:
W+1 data of analysis obtain the MOA leakage current letter The fundamental information of number IWith
Using formula:Calculate the frequency drift μ of the MOA leakage current signal I1
Using formulaCalculate the fundamental wave initial phase angle of the MOA leakage current signal I;
Using formulaLinearly correct the fundamental wave initial phase angle of the MOA leakage current signal I.
Using formulaCalculate the fundamental voltage amplitude of the MOA leakage current signal I;
Again, using formulaCalculate projected angle of the MOA leakage current fundamental wave on network voltage fundamental wave;
Finally, using formulaCalculate MOA resistive leakage current fundamental wave.
Those skilled in the art it should be appreciated that more than embodiment be intended merely to illustrate the present invention, and It is not intended as limitation of the invention, the present invention can also be changing into more modes, as long as in connotation model of the invention In enclosing, variation, the modification of embodiment described above will all be fallen within the scope of claims of the present invention.

Claims (1)

1. a kind of measurement method of the resistive current first harmonics of metal oxide arrester, it is characterised in that the following steps are included:
(1) the W+2 sample point data of synchronized sampling mains voltage signal V and MOA leakage current signal I at equal intervals: { fV(i), fI(i), i=0,1 ..., W+1 };
(2) plesiochronous DFT formula is applied since the sampled point i=0 of the mains voltage signal V:
W+1 data of analysis obtain the base of the mains voltage signal V Wave informationWith
Plesiochronous DFT formula is applied from the sampled point i=1 of the mains voltage signal V:
W+1 data of analysis obtain the base of the mains voltage signal V Wave informationWith
Using formula:Calculate the frequency drift μ of the mains voltage signal VV
Using formulaCalculate the fundamental wave initial phase angle of the mains voltage signal V;
Using formulaLinearly correct the fundamental wave initial phase angle of the mains voltage signal V;
(3) plesiochronous DFT formula is applied since the sampled point i=0 of the MOA leakage current signal I:
W+1 data of analysis obtain the MOA leakage current signal I Fundamental informationWith
Plesiochronous DFT formula is applied from the sampled point i=1 of the MOA leakage current signal I:
W+1 data of analysis obtain the MOA leakage current signal I Fundamental informationWith
Using formula:Calculate the frequency drift μ of the MOA leakage current signal II
Using formulaCalculate the fundamental wave initial phase angle of the MOA leakage current signal I;
Using formulaLinearly correct the fundamental wave initial phase angle of the MOA leakage current signal I;
Using formulaCalculate the fundamental voltage amplitude of the MOA leakage current signal I;
(4) formula is appliedCalculate projected angle of the MOA leakage current fundamental wave on network voltage fundamental wave;
(5) formula is appliedCalculate MOA resistive leakage current fundamental wave;
The synchronized sampling at equal intervals is according to the cycle T and frequency f of the ideal signal for carrying out frequency analysis, power frequency component frequency Rate f is 50Hz, cycle T 20mS, samples N point in one cycle, i.e. sample frequency is fs=Nf, and N >=64;
Described W+2 sample point data of sampling is to use iterative Simpson integration method, then W=n (N-1)/2);For the sum of all weighting coefficients.
CN201810704845.7A 2015-05-19 2015-05-19 Method for measuring resistive current fundamental wave of metal oxide arrester Active CN108957097B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810704845.7A CN108957097B (en) 2015-05-19 2015-05-19 Method for measuring resistive current fundamental wave of metal oxide arrester

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201810704845.7A CN108957097B (en) 2015-05-19 2015-05-19 Method for measuring resistive current fundamental wave of metal oxide arrester
CN201510257527.7A CN104849530B (en) 2015-05-19 2015-05-19 MOA resistive current fundamental wave measuring method

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201510257527.7A Division CN104849530B (en) 2015-05-19 2015-05-19 MOA resistive current fundamental wave measuring method

Publications (2)

Publication Number Publication Date
CN108957097A true CN108957097A (en) 2018-12-07
CN108957097B CN108957097B (en) 2020-10-09

Family

ID=53849323

Family Applications (3)

Application Number Title Priority Date Filing Date
CN201510257527.7A Active CN104849530B (en) 2015-05-19 2015-05-19 MOA resistive current fundamental wave measuring method
CN201810704845.7A Active CN108957097B (en) 2015-05-19 2015-05-19 Method for measuring resistive current fundamental wave of metal oxide arrester
CN201810704840.4A Active CN109490603B (en) 2015-05-19 2015-05-19 Method for measuring resistive current fundamental wave

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201510257527.7A Active CN104849530B (en) 2015-05-19 2015-05-19 MOA resistive current fundamental wave measuring method

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201810704840.4A Active CN109490603B (en) 2015-05-19 2015-05-19 Method for measuring resistive current fundamental wave

Country Status (1)

Country Link
CN (3) CN104849530B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110763930A (en) * 2019-09-28 2020-02-07 沈阳工程学院 Arrester resistive current on-line monitoring system based on Blackman bimodal interpolation harmonic analysis method
CN111323665A (en) * 2020-03-18 2020-06-23 合肥瀚度电力科技有限公司 Lightning arrester monitoring device, method and system based on GPS time service ratio correction
CN114994573A (en) * 2022-06-02 2022-09-02 国网安徽省电力有限公司马鞍山供电公司 Metal oxide lightning arrester fault diagnosis method based on artificial intelligence

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105044524B (en) * 2015-08-10 2018-12-04 许继集团有限公司 A kind of arrester monitoring method and system suitable for intelligent substation

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1235675A (en) * 1996-09-05 1999-11-17 西门子公司 Determination of fundamental and harmonic oscillations of measured electrical quantity
US6195328B1 (en) * 1998-04-15 2001-02-27 The United States Of America As Represented By The Secretary Of The Air Force Block adjustment of synchronizing signal for phase-coded signal tracking
CN102393488A (en) * 2011-08-24 2012-03-28 江苏技术师范学院 Harmonic analysis method
CN102435816A (en) * 2011-10-21 2012-05-02 江苏技术师范学院 Working method of MOA resistive current tester
CN102901856A (en) * 2012-09-17 2013-01-30 吉林省电力有限公司长春供电公司 Cable line arrester resistive current detection method based on phase search
CN103592512A (en) * 2011-10-21 2014-02-19 蒋春花 Harmonic analysis method of power quality harmonic analyzer
CN104062528A (en) * 2014-07-04 2014-09-24 武汉大学 Signal harmonic analysis method and system based on Hanning product window

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5832413A (en) * 1995-12-18 1998-11-03 Abb Power T&D Company Inc. Generator protection system and method for phasor estimation and frequency tracking during frequency ramping
US5832414A (en) * 1995-12-18 1998-11-03 Abb Power T&D Company Inc. Generator protection system and method of compensating for errors in phasor estimation due to oscillations in discrete Fourier transform
CN101430349B (en) * 2008-12-09 2012-04-18 东北电力科学研究院有限公司 Novel method for measuring resistive current of gapless metal oxide arrester
CN102435815B (en) * 2011-10-21 2014-02-19 江苏理工学院 Operating method of resistive current on-line monitoring system of metal oxide arrester (MOA)
CN103439566B (en) * 2011-10-21 2014-12-24 常州顺创电气科技有限公司 Operating method of MOA resistive current tester with relatively high precision
CN102721850A (en) * 2012-06-08 2012-10-10 上海市电力公司 Measuring method for fundamental wave component of resistive current of zinc oxide arrester
CN103207307B (en) * 2013-03-13 2015-07-08 福建省电力有限公司 Test method of leakage current component of metal oxide arrester
CN103267896B (en) * 2013-05-07 2016-04-20 江苏理工学院 Method for measuring initial phase angle of periodic signal
CN103257273B (en) * 2013-05-07 2016-05-04 江苏理工学院 Method for measuring phase difference of signals with same frequency period
CN203241464U (en) * 2013-05-07 2013-10-16 常州顺创电气科技有限公司 Current measurement device for lightning arresters and capacitive devices

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1235675A (en) * 1996-09-05 1999-11-17 西门子公司 Determination of fundamental and harmonic oscillations of measured electrical quantity
US6195328B1 (en) * 1998-04-15 2001-02-27 The United States Of America As Represented By The Secretary Of The Air Force Block adjustment of synchronizing signal for phase-coded signal tracking
CN102393488A (en) * 2011-08-24 2012-03-28 江苏技术师范学院 Harmonic analysis method
CN102435816A (en) * 2011-10-21 2012-05-02 江苏技术师范学院 Working method of MOA resistive current tester
CN103592512A (en) * 2011-10-21 2014-02-19 蒋春花 Harmonic analysis method of power quality harmonic analyzer
CN102901856A (en) * 2012-09-17 2013-01-30 吉林省电力有限公司长春供电公司 Cable line arrester resistive current detection method based on phase search
CN104062528A (en) * 2014-07-04 2014-09-24 武汉大学 Signal harmonic analysis method and system based on Hanning product window

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XIANZHONG DAI 等: "QUASI-SYNCHRONOUS SAMPLING ALGORITHM AND ITS APPLICATIONS-3.HIGH ACCURATE MEASUREMENT OF FREQUENCY, FREQUENCY DEVIATION AND PHASE ANGLE DIFFERENCE IN POWER SYSTEMS", 《PROCEEDINGS OF THE 1993 IEEE INSTRUMENTATION AND MEASUREMENT TECHNOLOGY CONFERENCE》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110763930A (en) * 2019-09-28 2020-02-07 沈阳工程学院 Arrester resistive current on-line monitoring system based on Blackman bimodal interpolation harmonic analysis method
CN111323665A (en) * 2020-03-18 2020-06-23 合肥瀚度电力科技有限公司 Lightning arrester monitoring device, method and system based on GPS time service ratio correction
CN111323665B (en) * 2020-03-18 2022-06-28 合肥瀚度电力科技有限公司 Lightning arrester monitoring device, method and system based on GPS time service ratio correction
CN114994573A (en) * 2022-06-02 2022-09-02 国网安徽省电力有限公司马鞍山供电公司 Metal oxide lightning arrester fault diagnosis method based on artificial intelligence

Also Published As

Publication number Publication date
CN104849530B (en) 2018-06-08
CN108957097B (en) 2020-10-09
CN109490603A (en) 2019-03-19
CN104849530A (en) 2015-08-19
CN109490603B (en) 2020-12-04

Similar Documents

Publication Publication Date Title
CN103439567B (en) A kind of method of work of MOA current in resistance property on-line monitoring system
CN106093593B (en) Condenser type electric power equipment dielectric loss online test method and system
Yao et al. Fast S-transform for time-varying voltage flicker analysis
CN104849530B (en) MOA resistive current fundamental wave measuring method
CN103439566B (en) Operating method of MOA resistive current tester with relatively high precision
CN102508026B (en) Harmonic wave analysis method for electric energy quality harmonic wave analyzer
CN113366327A (en) Arrester resistive leakage current detection device and detection method thereof
Zhao et al. Investigation on digital algorithm for on-line monitoring and diagnostics of metal oxide surge arrester based on an accurate model
CN104849569B (en) Dielectric loss measuring method
CN114002475B (en) Online monitoring method for resistive current of lightning arrester
Phung Computer-based partial discharge detection and characterisation
CN103969507B (en) Power quality harmonic analysis method
CN102435816B (en) Working method of MOA (Metal Oxide Arrester) resistive current tester
Zaro et al. Characterization of short-duration voltage events
CN103983852B (en) Harmonic analysis method of power quality harmonic analyzer
CN104849551B (en) Harmonic phase angle analysis method
CN103604989A (en) Harmonic wave analysis method of power-quality harmonic wave analyzer
CN103592512A (en) Harmonic analysis method of power quality harmonic analyzer
Xingang et al. Supraharmonics measurement algorithm based on CS-SAMP
CN110083891A (en) A kind of calculation method of the plesiochronous linear correction factor of DFT amplitude
KR20160110772A (en) Measuring algorithm for electric power element
Zorrozua et al. Estimation of power system harmonics and interharmonics in the presence of aperiodic components
Dasari Phase Locked Loop Based Signal Processing Approach for the Health Monitoring of Power Systems from their RF Emissions
Yi-Xiong et al. Research on the simple interpolated FFT algorithm for harmonic power energy measurement
CN105866622B (en) Two ends of electric transmission line fault distance-finding method is realized using zero-sequence component

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