CN104833878B - Power oscillation of power system determination methods based on synchronous phasor measuring device - Google Patents

Power oscillation of power system determination methods based on synchronous phasor measuring device Download PDF

Info

Publication number
CN104833878B
CN104833878B CN201510218681.3A CN201510218681A CN104833878B CN 104833878 B CN104833878 B CN 104833878B CN 201510218681 A CN201510218681 A CN 201510218681A CN 104833878 B CN104833878 B CN 104833878B
Authority
CN
China
Prior art keywords
power
mtd
oscillation
mtr
measuring device
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
Application number
CN201510218681.3A
Other languages
Chinese (zh)
Other versions
CN104833878A (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 Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
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 Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, Nari Technology Co Ltd, NARI Nanjing Control System Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510218681.3A priority Critical patent/CN104833878B/en
Publication of CN104833878A publication Critical patent/CN104833878A/en
Application granted granted Critical
Publication of CN104833878B publication Critical patent/CN104833878B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses a kind of power oscillation of power system determination methods based on synchronous phasor measuring device, including by analog signalses after anti-aliasing analog filtering, sampled value is stamped into accurate absolute time mark in analog-to-digital conversion;Active power, order components degree of asymmetry and voltage and power factor product are calculated using synchronization three-phase voltage current;Fast Fourier Transform (FFT) is carried out to the result of active power, obtains the frequency spectrum of active power;If spectrum energy of the frequency spectrum of active power between 0.1 2Hz is more than corresponding preset value, and voltage and power factor product are in corresponding preset range, order components degree of asymmetry is again smaller than corresponding preset value simultaneously, then showing power system, there occurs oscillation of power, outputting alarm.The present invention realizes the accurate judgement to oscillation of power, avoids short trouble being mistaken for system oscillation, or occurs to fail to judge during oscillation of power, available for raising synchronous phasor measuring device oscillation of power alarm performance.

Description

Power oscillation of power system determination methods based on synchronous phasor measuring device
Technical field
The present invention relates to a kind of power oscillation of power system determination methods based on synchronous phasor measuring device, belong to electric power System detectio technical field.
Background technology
Power oscillation of power system have become influence modern large-scale power system safe and stable operation major issue it One.In recent years, there occurs multiple oscillation of power phenomenon for China's power system.The main feature of these vibrations is as follows:1) oscillation frequency Rate is low, between 0.2~2Hz, covers main low-frequency oscillation frequency range;2) amplitude is larger, and peak-to-peak value is tens to hundreds of Between megawatt;3) duration is grown, and at tens seconds between dozens of minutes, most of duration of oscillation is in a few minutes.
Synchronous phasor measurement unit (phasor measurement unit, PMU) is used as WAMS (wide Area measurement system, WAMS) data source, can timely find power network oscillation of power have it is important Meaning.
Waveform during power oscillation of power system, usually exponential damping waveform, includes multi-frequency component.At present The upper existing oscillation of power determination methods of PMU are primarily present problems with:
1) common oscillation of power determination methods are led to overpowering fluctuation size and realized, load fluctuation, short trouble and Wrong report is easily caused during circuit unit switching;
2) algorithm is difficult to quantify to the degree of oscillation of power at present, in order to avoid wrong report is then likely to cause oscillation amplitude not Failed to report when big, vice versa.
The content of the invention
It is an object of the invention to overcome deficiency of the prior art, there is provided a kind of electricity based on synchronous phasor measuring device Force system oscillation of power determination methods, solve power oscillation of power system determination methods in the prior art and easily occur reporting by mistake, fail to report Technical problem.
In order to solve the above technical problems, the technical solution adopted in the present invention is:Electricity based on synchronous phasor measuring device Force system oscillation of power determination methods, comprise the following steps:
Step 1:By the analog signalses of synchronous phasor measuring device input after anti-aliasing analog filtering, mould is carried out Number conversion, stamps accurate absolute time mark, and be stored in data buffer zone in analog-to-digital conversion by sampled value;
Step 2:The sampled value of data buffer zone according to caused by step 1, filled using synchronization synchronous phasor measurement The three-phase voltage current put calculates active-power P, order components degree of asymmetry and voltage and power factor product U cos φ;
Step 3:The result of the active-power P calculated step 2 carries out Fast Fourier Transform (FFT), obtains active power Frequency spectrum;
Step 4:Active-power P, order components degree of asymmetry, voltage and the power factor product calculated according to step 2 Ucos φ, the frequency spectrum of the active power obtained with reference to step 3 judge whether power system occurs oscillation of power:If wattful power Spectrum energy of the frequency spectrum of rate between 0.1-2Hz is more than corresponding preset value, and voltage and power factor product Ucos φ are in In corresponding preset range, while order components degree of asymmetry (| I0|+|I2|)/m|I1| again smaller than corresponding preset value, then show electric power There occurs oscillation of power, outputting alarm for system;Otherwise, show that oscillation of power does not occur for power system.
Sampled value is stamped into accurate absolute time calibration method in step 1 is:
By GPS/ big dipper clocks pair when signal access synchronous phasor measuring device;
When analog-digital converter samples, sampled value is stamped to the clock signal of GPS/ big dipper clocks output, and is stored in data Buffering area.
Active-power P is calculated using formula (1) in step 2:
P (t)=ua(t)ia(t)+ub(t)ib(t)+uc(t)ic(t); (1)
In formula:P (t) represents the active power of t, ua(t)、ub(t)、uc(t) t A phase, B phase, C phase is represented respectively Instantaneous voltage, ia(t)、ib(t)、ic(t) transient current of t A phase, B phase, C phase is represented respectively.
Fast Fourier Transform (FFT) is carried out using the change of Cooley-Tu Ji fast Fouriers to the result of active-power P in step 3 Scaling method.
The computational methods of voltage and power factor product Ucos φ described in step 2 are as follows:Voltage and electricity are calculated first The amplitude phase of stream, the phase difference of voltage and electric current is then calculated, be i.e. φ, powerfactorcosφ is tried to achieve, finally by power factor Cos φ are multiplied to obtain Ucos φ with voltage magnitude.
Order components degree of asymmetry described in step 2 are calculated using formula (2):
(|I0|+|I2|)/m|I1| (2)
In formula, m takes 0.5-0.7 for that can adjust definite value;Calculation formula it is as follows:
Compared with prior art, the beneficial effect that is reached of the present invention is:Multiplied using power spectrum, voltage and power factor Product, order components degree of asymmetry are judged oscillation of power, avoid causing when load fluctuation, short trouble and circuit unit switching Wrong report, has also carried out a certain degree of quantization to oscillation of power, it is therefore prevented that oscillation of power is failed to report, and is surveyed available for synchronized phasor is improved Measure rating of set oscillating alarm performance.
Brief description of the drawings
Fig. 1 is the flow chart of the inventive method.
Fig. 2 is the system construction drawing of simulation calculation.
Fig. 3 is that the voltage and current waveform that sampling obtains is synchronized to the system architecture shown in Fig. 2.
Fig. 4 is the active power oscillogram that simulation calculation acquisition is carried out to the system architecture shown in Fig. 2.
Fig. 5 is the active power spectrogram that simulation calculation acquisition is carried out to the system architecture shown in Fig. 2.
Fig. 6 is that the perunit value positive sequence voltage and power factor of the system architecture progress simulation calculation acquisition shown in Fig. 2 are multiplied Product oscillogram.
Fig. 7 is the order components degree of asymmetry figure that simulation calculation acquisition is carried out to the system architecture shown in Fig. 2.
Embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following examples are only used for clearly illustrating the present invention Technical scheme, and can not be limited the scope of the invention with this.
As shown in figure 1, being the flow chart of the inventive method, Fig. 2 is the system construction drawing of simulation calculation.Based on synchronized phasor The power oscillation of power system determination methods of measurement apparatus, comprise the following steps:
Step 1:By the analog signalses of synchronous phasor measuring device input after anti-aliasing analog filtering, mould is carried out Number conversion, by GPS/ big dipper clocks pair when signal access synchronous phasor measuring device, when analog-digital converter samples, will sample Value stamps the clock signal of GPS/ big dipper clocks output, and is stored in data buffer zone.As shown in figure 3, it is to the system shown in Fig. 2 Structure synchronizes the voltage and current waveform that sampling obtains.
Step 2:The sampled value of data buffer zone according to caused by step 1, filled using synchronization synchronous phasor measurement The three-phase voltage current put calculates active-power P, order components degree of asymmetry and voltage and power factor product Ucos φ.
As shown in figure 4, be the active power oscillogram that simulation calculation acquisition is carried out to the system architecture shown in Fig. 2, it is active Power P is calculated using formula (1):
P (t)=ua(t)ia(t)+ub(t)ib(t)+uc(t)ic(t); (1)
In formula:P (t) represents the active power of t, ua(t)、ub(t)、uc(t) t A phase, B phase, C phase is represented respectively Instantaneous voltage, ia(t)、ib(t)、ic(t) transient current of t A phase, B phase, C phase is represented respectively.
As shown in fig. 7, the order components degree of asymmetry figure of simulation calculation acquisition is carried out to the system architecture shown in Fig. 2, sequence point Amount degree of asymmetry is calculated using formula (2):
(|I0|+|I2|)/m|I1| (2)
In formula, m takes 0.5-0.7, preferably 0.66 for that can adjust definite value;Calculation formula it is as follows:
As shown in fig. 6, to shown in Fig. 2 system architecture carry out simulation calculation acquisition perunit value positive sequence voltage and power because Number product waveform figure, voltage and power factor product Ucos φ computational methods are as follows:The amplitude of voltage and current is calculated first Phase, then calculate the phase difference of voltage and electric current, be i.e. φ, try to achieve powerfactorcosφ, finally by powerfactorcosφ with Voltage magnitude is multiplied to obtain Ucos φ.
Step 3:The result of the active-power P calculated step 2 carries out Cooley-Tu Ji Fast Fourier Transform (FFT)s, obtains The frequency spectrum of active power, as shown in figure 5, carrying out the active power spectrogram of simulation calculation acquisition to the system architecture shown in Fig. 2. Specific calculating process is as follows:
Fourier transformation is shown below:
In formula:xkBeing inputted for Fourier transformation, N is data window length,For twiddle factor, ynFor Fourier transformation N-th of component of output, using n times unit root WNTo representWNProperty be:
1) periodically, WNWith cycle N, i.e.,
2) symmetry:
If 3) p is N approximate number,Serial length N=2 to be transformedr, r is positive integer.According to
The symmetry of unit root above, seeks seriesWhen, summation section can be divided into two parts:
FoddAnd F (k)even(k) it is two respectively about serialOdd number and the serial N/2 point transformation of even number.By This formula can only calculate ykPreceding N/2 point, for rear N/2 point, pay attention to FoddAnd F (k)even(k) all it is letter of the cycle for N/2 Number, by the symmetry of unit root, then there is following transformation for mula:
So, a N point transformation has just resolved into two N/2 point transformation, can continue decomposition after this manner and go down, and obtains final Transformation results.
Step 4:Active-power P, order components degree of asymmetry, voltage and the power factor product calculated according to step 2 Ucos φ, the frequency spectrum of the active power obtained with reference to step 3 judge whether power system occurs oscillation of power:If wattful power Spectrum energy of the frequency spectrum of rate between 0.1-2Hz is more than corresponding preset value, and voltage and power factor product Ucos φ are in In corresponding preset range, while order components degree of asymmetry (| I0|+|I2|)/m|I1| again smaller than corresponding preset value, then show electric power There occurs oscillation of power, outputting alarm for system;Otherwise, show that oscillation of power does not occur for power system.
It the above is only the preferred embodiment of the present invention, it is noted that come for those skilled in the art Say, without departing from the technical principles of the invention, some improvement and deformation can also be made, these are improved and deformation also should It is considered as protection scope of the present invention.

Claims (6)

1. the power oscillation of power system determination methods based on synchronous phasor measuring device, it is characterised in that comprise the following steps:
Step 1:By the analog signalses of synchronous phasor measuring device input after anti-aliasing analog filtering, carry out modulus and turn Change, sampled value is stamped into accurate absolute time mark in analog-to-digital conversion, and be stored in data buffer zone;
Step 2:The sampled value of data buffer zone according to caused by step 1, utilizes synchronization synchronous phasor measuring device Three-phase voltage current calculates active-power P, order components degree of asymmetry and voltage and power factor product U cos φ;
Step 3:The result of the active-power P calculated step 2 carries out Fast Fourier Transform (FFT), obtains the frequency of active power Spectrum;
Step 4:Active-power P, order components degree of asymmetry, voltage and the power factor product U cos calculated according to step 2 φ, the frequency spectrum of the active power obtained with reference to step 3 judge whether power system occurs oscillation of power:If active power Spectrum energy of the frequency spectrum between 0.1-2Hz is more than corresponding preset value, and voltage and power factor product U cos φ are in correspondingly In preset range, while order components degree of asymmetry (| I0|+|I2|)/m|I1| again smaller than corresponding preset value, then show power system There occurs oscillation of power, outputting alarm;Otherwise, show that oscillation of power does not occur for power system.
2. the power oscillation of power system determination methods according to claim 1 based on synchronous phasor measuring device, it is special Sign is that sampled value is stamped into accurate absolute time calibration method in step 1 is:
By GPS/ big dipper clocks pair when signal access synchronous phasor measuring device;
When analog-digital converter samples, sampled value is stamped to the clock signal of GPS/ big dipper clocks output, and is stored in data buffering Area.
3. the power oscillation of power system determination methods according to claim 1 based on synchronous phasor measuring device, it is special Sign is that active-power P is calculated using formula (1) in step 2:
P (t)=ua(t)ia(t)+ub(t)ib(t)+uc(t)ic(t); (1)
In formula:P (t) represents the active power of t, ua(t)、ub(t)、uc(t) wink of t A phase, B phase, C phase is represented respectively When voltage, ia(t)、ib(t)、ic(t) transient current of t A phase, B phase, C phase is represented respectively.
4. the power oscillation of power system determination methods according to claim 1 based on synchronous phasor measuring device, it is special Sign is that the result progress Fast Fourier Transform (FFT) in step 3 to active-power P uses Cooley-Tu Ji Fast Fourier Transform (FFT)s Algorithm.
5. the power oscillation of power system determination methods according to claim 1 based on synchronous phasor measuring device, it is special Sign is that the computational methods of voltage and power factor product U cos φ described in step 2 are as follows:Voltage and electricity are calculated first The amplitude phase of stream, the phase difference of voltage and electric current is then calculated, be i.e. φ, powerfactorcosφ is tried to achieve, finally by power factor Cos φ are multiplied to obtain U cos φ with voltage magnitude.
6. the power oscillation of power system determination methods according to claim 1 based on synchronous phasor measuring device, it is special Sign is that order components degree of asymmetry described in step 2 are calculated using formula (2):
(|I0|+|I2|)/m|I1| (2)
In formula, m takes 0.5-0.7 for that can adjust definite value;Calculation formula it is as follows:
<mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <mo>&amp;CenterDot;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mo>/</mo> <mn>3</mn> </mrow> </msup> </mtd> <mtd> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mo>/</mo> <mn>3</mn> </mrow> </msup> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mo>/</mo> <mn>3</mn> </mrow> </msup> </mtd> <mtd> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mo>/</mo> <mn>3</mn> </mrow> </msup> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>b</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>c</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow> 1
CN201510218681.3A 2015-04-30 2015-04-30 Power oscillation of power system determination methods based on synchronous phasor measuring device Active CN104833878B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510218681.3A CN104833878B (en) 2015-04-30 2015-04-30 Power oscillation of power system determination methods based on synchronous phasor measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510218681.3A CN104833878B (en) 2015-04-30 2015-04-30 Power oscillation of power system determination methods based on synchronous phasor measuring device

Publications (2)

Publication Number Publication Date
CN104833878A CN104833878A (en) 2015-08-12
CN104833878B true CN104833878B (en) 2017-12-01

Family

ID=53811882

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510218681.3A Active CN104833878B (en) 2015-04-30 2015-04-30 Power oscillation of power system determination methods based on synchronous phasor measuring device

Country Status (1)

Country Link
CN (1) CN104833878B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106199183B (en) * 2016-08-16 2018-09-28 国电南瑞科技股份有限公司 A kind of PMU and method realizing sub-synchronous oscillation on-line identification alarm

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101944776A (en) * 2010-09-28 2011-01-12 电子科技大学 Low-frequency oscillation monitoring system of electrical network
CN103248071A (en) * 2013-05-15 2013-08-14 国家电网公司 Generator set low frequency oscillation monitoring method based on PMU measuring data
CN103336909A (en) * 2013-07-23 2013-10-02 国家电网公司 Low frequency oscillation identification method for accessing of wind power to power grid
CN104217112A (en) * 2014-09-02 2014-12-17 南京工程学院 Multi-type signal-based power system low-frequency oscillation analysis method
CN104392141A (en) * 2014-12-05 2015-03-04 广东电网有限责任公司电力科学研究院 Method and device for locating low-frequency oscillation disturbance source of electric power system

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201001007D0 (en) * 2010-01-22 2010-03-10 Ulive Entpr Ltd Apparatus and method for measuring a phasor of an electrical power system signal
US9043040B2 (en) * 2012-01-31 2015-05-26 General Electric Company Power system stabilization

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101944776A (en) * 2010-09-28 2011-01-12 电子科技大学 Low-frequency oscillation monitoring system of electrical network
CN103248071A (en) * 2013-05-15 2013-08-14 国家电网公司 Generator set low frequency oscillation monitoring method based on PMU measuring data
CN103336909A (en) * 2013-07-23 2013-10-02 国家电网公司 Low frequency oscillation identification method for accessing of wind power to power grid
CN104217112A (en) * 2014-09-02 2014-12-17 南京工程学院 Multi-type signal-based power system low-frequency oscillation analysis method
CN104392141A (en) * 2014-12-05 2015-03-04 广东电网有限责任公司电力科学研究院 Method and device for locating low-frequency oscillation disturbance source of electric power system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Experimental studies of three dimensional fuzzy logic power system stabilizer on damping of low-frequency global mode of oscillation;Takashi Hiyama等;《Fuzzy Sets and Systems》;19991231;103-111 *
故障系统低频振荡特征值分析方法;马静等;《电力自动化设备》;20140430;第34卷(第4期);13-18 *

Also Published As

Publication number Publication date
CN104833878A (en) 2015-08-12

Similar Documents

Publication Publication Date Title
CN105223418B (en) The measurement method and measuring device of subsynchronous and supersynchronous harmonic phasor
US20140125314A1 (en) Apparatus and method for alternating current physical signals measurement and data acquisition
CN106841778A (en) The processing method of the subsynchronous and supersynchronous harmonic parameters realized based on PMU
US9488681B2 (en) Convolution integral for synchronized phasor
CN103033789A (en) Static and dynamic test system of synchronous phasor measurement unit pressure measuring unit (PMU)
CN106932642B (en) Electric Power Harmonic Analysis method
CN102253282B (en) Method for obtaining continuous frequency spectrum interpolation power harmonic parameter of Nuttall window function
CN103926462B (en) Rapid harmonic wave analyzing method and device of power system
CN104198977B (en) Accuracy detection method based on average power error for analog input combining unit
JP2013044752A (en) Phase identification system and method
CN107064630B (en) A kind of power system frequency measurement method and device
CN102298099A (en) Power distribution and consumption monitoring and measuring apparatus possessing time calibration function
CN104502707A (en) Synchronized phasor measurement method for electrical power system based on cubic spline interpolation
CN106645919A (en) Power grid full spectrum power oscillation phasor synchronization measurement method based on three-phase instantaneous power
CN104181391A (en) Harmonic detection method of digital power meter
CN109490630A (en) A kind of dynamic phasor measurement method based on pencil of matrix
KR101034259B1 (en) Apparatus for Synchronized Phasor Measurement Unit Using The Re-sampling Technique
CN103605038A (en) Protective error detecting and verifying system of electronic current transformer
CN101504442B (en) Real-time detection method for double-SRF dual-feedback aerogenerator rotor current positive and negative components
CN104833878B (en) Power oscillation of power system determination methods based on synchronous phasor measuring device
Rodrigues et al. A Phasor Measurement Unit based on discrete fourier transform using digital signal processor
CN102169138A (en) Method for processing correction of phase difference of power grid test or counting device
CN101949974A (en) Field tester of alternating current electric energy meter
CN202886585U (en) On-line test device for voltage/current transformer
CN102749488A (en) Power grid harmonic wave real-time on-line monitor and method for detecting harmonic wave using same

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant