CN104392141A - Method and device for locating low-frequency oscillation disturbance source of electric power system - Google Patents
Method and device for locating low-frequency oscillation disturbance source of electric power system Download PDFInfo
- Publication number
- CN104392141A CN104392141A CN201410740740.9A CN201410740740A CN104392141A CN 104392141 A CN104392141 A CN 104392141A CN 201410740740 A CN201410740740 A CN 201410740740A CN 104392141 A CN104392141 A CN 104392141A
- Authority
- CN
- China
- Prior art keywords
- mode function
- intrinsic mode
- function imf
- generator
- 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.)
- Pending
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
Disclosed are a method and a device for locating a low-frequency oscillation disturbance source of an electric power system. The method includes: A, acquiring real-time operation data of electric generator power angles; B, subjecting the real-time operation data of the electric generator power angles to empirical mode decomposition so as to acquire a dominant oscillation mode IMF (intrinsic mode function); C, subjecting the dominant oscillation mode IMF to Hilbert transformation so as to acquire instantaneous amplitude values of the electric generator power angles; D, determining a disturbance source position according to the sudden change sequence of the instantaneous amplitude values of the electric generator power angles. By the method and the device for locating the low-frequency oscillation disturbance source of the electric power system, sudden change features of the electric generator power angle amplitude values in the initial oscillation phase can be acquired accurately, and thus, a generator set where the disturbance source is located can be recognized quickly.
Description
Technical field
The present invention relates to field of power, particularly relate to the fault detection and diagnosis technology of electric system.
Background technology
The low-frequency oscillation of electric system refers to that low-frequency oscillation refers to electromechanical oscillations lasting between electric generating set rotor kinetic energy and electrical network potential energy, mechanically show as the long-term swing of complementary group of planes rotor relative angle, on electrically, then show as the oscillation of power between this complementary group of planes stator.Low-frequency oscillation can cause interconnection overcurrent tripping or system and system or the step-out between unit and system and off-the-line, and serious threat, to the stability of electric system, should reduce as far as possible or eliminate low-frequency oscillation in the control procedure of therefore electric system.
Force low-frequency oscillation to have clear and definite disturbing source for the electric system caused by durations disturbance, once described disturbing source is eliminated, low-frequency oscillation can be lowered very soon and even eliminate.Relevant electric system is at present eliminated and is forced the method for low-frequency oscillation or device analysis method mainly to concentrate on the type analyzed and can cause the disturbing source of determining to force low-frequency oscillation, and is carried out the method for low-frequency oscillation disturbance source locating by energy function.Existing low-frequency oscillation disturbance source locating method positioning principle is complicated, and computing velocity is slow, is difficult to meet the requirement that low-frequency oscillation of electric power system disturbing source identifies online fast.
Summary of the invention
Given this, the object of the invention is to the position fast, accurately determining low-frequency oscillation of electric power system disturbing source.
In order to realize this object, the technical scheme that the present invention takes is as follows.
A kind of low-frequency oscillation of electric power system disturbance source locating method, described method comprises step:
A, acquisition generator's power and angle real-time running data;
B, empirical mode decomposition (Empirical ModeDecomposition, EMD) is carried out to described generator's power and angle real-time running data obtain control oscillation modes intrinsic mode function (Intrinsic Mode Function, IMF);
C, Hilbert transform is carried out to control oscillation modes intrinsic mode function IMF, obtain the instantaneous amplitude of generator's power and angle;
D, the sequencing determination disturbing source position of undergoing mutation with the instantaneous amplitude of generator's power and angle.
Wherein, generator's power and angle real-time running data is obtained by WAMS WAMS.
The method obtaining control oscillation modes intrinsic mode function IMF is:
B1, with generator's power and angle real-time running data for calculating data x (t), all maximum points and minimum point that utilize 3 SPL to connect to calculate data x (t), form upper and lower envelope v
1(t) and v
2(t), definition
if
meeting intrinsic mode function IMF condition is then first intrinsic mode function IMF component C with it
1, otherwise right
continue screening process, until after k screening
meet intrinsic mode function IMF condition, with
be first intrinsic mode function IMF component C
1;
B2, data and (n-1)th intrinsic mode function IMF component C will be calculated
n-1difference as new calculating data, and repeated execution of steps B1, obtains the n-th intrinsic mode function IMF component C
n, until new calculating data are monotonic quantity, or stop screening when new calculating data are less than predetermined value;
B3, take mode of oscillation intrinsic mode function IMF as the leading factor with the intrinsic mode function IMF component of dominant oscillating mode.
Especially, in step bl is determined., described intrinsic mode function IMF condition is: the average of the extreme point of intrinsic mode function IMF in data segment and zero number is equal or difference 1, intrinsic mode function IMF is order at any time envelope is up and down zero.
On the other hand, in step bl is determined., the described method continuing screening is:
Jth is screened for-1 time
deduct
the average envelope of envelope up and down, form jth time screening
wherein j>=2.
A kind of low-frequency oscillation of electric power system disturbance source locating device, comprising:
Measuring unit, for obtaining generator's power and angle real-time running data;
Control oscillation modes intrinsic mode function IMF acquiring unit, obtains control oscillation modes intrinsic mode function IMF for carrying out empirical mode decomposition to described real-time running data;
Instantaneous amplitude analytic unit, for carrying out Hilbert transform to control oscillation modes intrinsic mode function IMF, obtains the instantaneous amplitude of generator's power and angle;
Disturbance source locating unit, for the sequencing determination disturbing source position of undergoing mutation with the instantaneous amplitude of generator's power and angle.
By adopting low-frequency oscillation of electric power system disturbance source locating method and apparatus of the present invention, empirical mode decomposition EMD method is adopted to obtain the control oscillation modes curve of generator's power and angle.EMD method has the ability of filtering appts high fdrequency component and DC component, and noise robustness is good, and execution speed is very fast, and can analyze the non-stationary signal that parameter of oscillation changes in time.Therefore can obtain the catastrophe characteristics of generator's power and angle amplitude at vibration initial period exactly by carrying out Hilbert transform to the intrinsic mode function IMF proposed, identify the genset at disturbing source place fast.The present invention not only principle is simple, and practical, is convenient to the ONLINE RECOGNITION realizing low-frequency oscillation disturbing source.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of embodiment of the present invention low-frequency oscillation of electric power system disturbance source locating method.
Fig. 2 is the circuit topological structure figure of an embodiment of the present invention application example.
Fig. 3 is the result schematic diagram of embodiment of the present invention low-frequency oscillation of electric power system disturbance source locating method.
Embodiment
Below in conjunction with accompanying drawing, the present invention is elaborated.
The example embodiment that following discloses are detailed.But concrete structure disclosed herein and function detail are only the objects for describing example embodiment.
But should be appreciated that, the present invention is not limited to disclosed concrete example embodiment, but covers all modifications, equivalent and the alternative that fall within the scope of the disclosure.In the description to whole accompanying drawing, identical Reference numeral represents identical element.
Should be appreciated that, term "and/or" as used in this comprises one or morely relevant lists any of item and all combinations simultaneously.Should be appreciated that in addition, when parts or unit are called as " connection " or " coupling " to another parts or unit, it can be directly connected or coupled to miscellaneous part or unit, or also can there is intermediate member or unit.In addition, other words being used for describing relation between parts or unit should be understood according to identical mode (such as, " between " to " directly ", " adjacent " to " direct neighbor " etc.).
In order to realize object of the present invention, the process flow diagram of low-frequency oscillation of electric power system disturbance source locating method of the present invention as shown in Figure 1, specifically comprises the following steps:
A, acquisition generator's power and angle real-time running data x (t);
B, described generator's power and angle real-time running data x (t) carried out to empirical mode decomposition EMD and obtain control oscillation modes intrinsic mode function IMF;
C, Hilbert transform is carried out to control oscillation modes intrinsic mode function IMF, obtain the instantaneous amplitude of generator's power and angle;
D, the sequencing determination disturbing source position of undergoing mutation with the instantaneous amplitude of generator's power and angle.
Generator's power and angle real-time running data can pass through sensor measurement gained, it will be appreciated by those skilled in the art that its usual measuring principle.In one more specifically embodiment, obtain generator's power and angle real-time running data by WAMS WAMS, but those skilled in the art are to be understood that this can not be formed limits application of the present invention.
By above method, can be located rapidly the place of low-frequency oscillation of electric power system disturbing source by the feature of intrinsic mode function IMF, the protection for electric system has obvious meaning.
In the present invention's specific embodiment, the method for the acquisition control oscillation modes intrinsic mode function IMF of step B is:
B1, with generator's power and angle real-time running data for calculating data x (t), all maximum points and minimum point that utilize 3 SPL to connect to calculate data x (t), form upper and lower envelope v
1(t) and v
2(t), definition
if
meeting intrinsic mode function IMF condition is then first intrinsic mode function IMF component C with it
1, otherwise right
continue screening process, until after k screening
meet intrinsic mode function IMF condition, with
be first intrinsic mode function IMF component C
1;
B2, data and (n-1)th intrinsic mode function IMF component C will be calculated
n-1difference as new calculating data, repeated execution of steps B1, obtains the n-th intrinsic mode function IMF component C
n, until new calculating data are monotonic quantity, or stop screening when new calculating data are less than predetermined value, described predetermined value requires to determine according to the realistic accuracy of system, and those skilled in the art all know in calculating and specifically choose mode.
Now calculate data
be the calculating data and the n-th intrinsic mode function IMF component C that calculate for n-th time
ndifference.
B3, take mode of oscillation intrinsic mode function IMF as the leading factor with the intrinsic mode function IMF component of dominant oscillating mode.
In one embodiment, described intrinsic mode function IMF condition is:
1, the extreme point of intrinsic mode function IMF in data segment is equal with zero number or differ from 1;
2, the average of envelope up and down that intrinsic mode function IMF orders at any time is zero.
In one embodiment, as shown in Figure 1, the described method continuing screening is: described by jth-1 screening
deduct
the average envelope of envelope up and down, form jth time screening
wherein j>=2.
In one embodiment of the present invention, describedly carry out Hilbert transform to control oscillation modes intrinsic mode function IMF, the method obtaining the instantaneous amplitude of generator's power and angle is:
Suppose that control oscillation modes intrinsic mode function IMF is X (t), Hilbert carried out to it and is transformed to:
And then the analytic signal that can obtain X (t) is:
Z(t)=X(t)+iY(t)=a(t)e
jφ(t)
In formula, a (t), φ (t) are respectively instantaneous amplitude and instantaneous phase function.
Wherein a (t) computing method are
Due to the starting of oscillation stage in low-frequency oscillation, the change of the merit angle of the genset at disturbing source place will prior to the genset at non-disturbing source place, therefore the instantaneous amplitude at the merit angle of the genset at disturbing source place is first undergone mutation.Therefore in each genset merit angle instantaneous amplitude change curve obtained according to step C, the priority moment of amplitude sudden change can differentiate the position at disturbing source place.
Fig. 2 is the schematic diagram of the circuit topological structure of an embodiment of the present invention application example.This application example is 3 machine 9 node systems, is the disturbing source of 1.0Hz by superposing a frequency on the machine torque of generator G2, thus initiating system generation low-frequency oscillation.
First, obtain three genset merit angle real-time running data x (t) separately, these data are obtained by WAMS WAMS;
Next, the control oscillation modes intrinsic mode function IMF that empirical mode decomposition EMD obtains each generator is carried out to three genset merit angle real-time running data x (t) separately, then the time dependent curve of merit angle instantaneous amplitude that Hilbert converts acquisition three genset is carried out to the control oscillation modes intrinsic mode function IMF of each genset, as shown in Figure 3.
There is the instantaneous amplitude sudden change of merit angle at 0.26s in genset G2 as can be seen from Figure 3, and genset G1, at 0.43s, the instantaneous amplitude sudden change of merit angle occurs, and genset G3, at 0.37s, the instantaneous amplitude sudden change of merit angle occurs, and they finally all reach stable oscillation stationary vibration amplitude.Therefore, obviously can find out that genset G2 occurs obviously early than genset G1 and genset G3 the moment of merit angle instantaneous amplitude sudden change.Therefore my Russian genset G2 place of disturbing source of low-frequency oscillation can be judged.Therefore can obtain the catastrophe characteristics of generator's power and angle instantaneous amplitude at vibration initial period exactly by the present invention, identify the genset at disturbing source place fast.
On the other hand, in order to realize low-frequency oscillation of electric power system disturbance source locating method of the present invention, further comprises a kind of low-frequency oscillation of electric power system disturbance source locating device in embodiment of the present invention, described device comprises:
Measuring unit, for obtaining generator's power and angle real-time running data;
Control oscillation modes intrinsic mode function IMF acquiring unit, obtains control oscillation modes intrinsic mode function IMF for carrying out empirical mode decomposition to described generator's power and angle real-time running data;
Instantaneous amplitude analytic unit, for carrying out Hilbert transform to control oscillation modes intrinsic mode function IMF, obtains the instantaneous amplitude of generator's power and angle;
Disturbance source locating unit, for the sequencing determination disturbing source position of undergoing mutation with the instantaneous amplitude of generator's power and angle.
It should be noted that; above-mentioned embodiment is only the present invention's preferably embodiment; can not limiting the scope of the invention be understood as, not depart under concept thereof of the present invention, all protection scope of the present invention is belonged to modification to any subtle change that the present invention does.
Claims (6)
1. a low-frequency oscillation of electric power system disturbance source locating method, described method comprises step:
A, acquisition generator's power and angle real-time running data;
B, described generator's power and angle real-time running data carried out to empirical mode decomposition and obtain control oscillation modes intrinsic mode function IMF;
C, Hilbert transform is carried out to control oscillation modes intrinsic mode function IMF, obtain the instantaneous amplitude of generator's power and angle;
D, the sequencing determination disturbing source position of undergoing mutation with the instantaneous amplitude of generator's power and angle.
2. the low-frequency oscillation of electric power system disturbance source locating method described in claim 1, is characterized in that, obtains generator's power and angle real-time running data by WAMS WAMS.
3. the low-frequency oscillation of electric power system disturbance source locating method described in claim 1, is characterized in that, the method obtaining control oscillation modes intrinsic mode function IMF is:
B1, with generator's power and angle real-time running data for calculating data x (t), all maximum points and minimum point that utilize 3 SPL to connect to calculate data x (t), form upper and lower envelope v
1(t) and v
2(t), definition
if
meeting intrinsic mode function IMF condition is then first intrinsic mode function IMF component C with it
1, otherwise right
continue screening process, until after k screening
meet intrinsic mode function IMF condition, with
be first intrinsic mode function IMF component C
1;
B2, data and (n-1)th intrinsic mode function IMF component C will be calculated
n-1difference as new calculating data, and repeated execution of steps B1, obtains the n-th intrinsic mode function IMF component C
n, until new calculating data are monotonic quantity, or stop screening when new calculating data are less than predetermined value;
B3, take mode of oscillation intrinsic mode function IMF as the leading factor with the intrinsic mode function IMF component of dominant oscillating mode.
4. the low-frequency oscillation of electric power system disturbance source locating method described in claim 3, it is characterized in that, in step B1, described intrinsic mode function IMF condition is: the average of the extreme point of intrinsic mode function IMF in data segment and zero number is equal or difference 1, intrinsic mode function IMF is order at any time envelope is up and down zero.
5. the low-frequency oscillation of electric power system disturbance source locating method described in claim 3, is characterized in that, in step B1, the described method continuing screening is:
Jth is screened for-1 time
deduct
the average envelope of envelope up and down, form jth time screening
wherein j>=2.
6. a low-frequency oscillation of electric power system disturbance source locating device, comprising:
Measuring unit, for obtaining generator's power and angle real-time running data;
Control oscillation modes intrinsic mode function IMF acquiring unit, obtains control oscillation modes intrinsic mode function IMF for carrying out empirical mode decomposition to described real-time running data;
Instantaneous amplitude analytic unit, for carrying out Hilbert transform to control oscillation modes intrinsic mode function IMF, obtains the instantaneous amplitude of generator's power and angle;
Disturbance source locating unit, for the sequencing determination disturbing source position of undergoing mutation with the instantaneous amplitude of generator's power and angle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410740740.9A CN104392141A (en) | 2014-12-05 | 2014-12-05 | Method and device for locating low-frequency oscillation disturbance source of electric power system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410740740.9A CN104392141A (en) | 2014-12-05 | 2014-12-05 | Method and device for locating low-frequency oscillation disturbance source of electric power system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104392141A true CN104392141A (en) | 2015-03-04 |
Family
ID=52610043
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410740740.9A Pending CN104392141A (en) | 2014-12-05 | 2014-12-05 | Method and device for locating low-frequency oscillation disturbance source of electric power system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104392141A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104833878A (en) * | 2015-04-30 | 2015-08-12 | 国电南瑞科技股份有限公司 | Electric power system power oscillation judgment method based on synchronized phasor measurement unit |
CN106569086A (en) * | 2015-10-08 | 2017-04-19 | 国家电网公司 | Method for on-site positioning of low-frequency oscillation disturbance source based on oscillation phase difference |
CN110492497A (en) * | 2019-08-16 | 2019-11-22 | 国网河南省电力公司电力科学研究院 | A kind of generating set power oscillation error comprehensive diagnosis method and system |
CN112072650A (en) * | 2020-09-02 | 2020-12-11 | 深圳大学 | Power system power oscillation path identification method, system and storage medium |
CN112103939A (en) * | 2019-06-18 | 2020-12-18 | 上海电机学院 | Method and device for analyzing stability of power system |
CN115577226A (en) * | 2022-12-08 | 2023-01-06 | 南方电网数字电网研究院有限公司 | Forced oscillation source positioning method, device, equipment and storage medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030033094A1 (en) * | 2001-02-14 | 2003-02-13 | Huang Norden E. | Empirical mode decomposition for analyzing acoustical signals |
CN101048677A (en) * | 2004-10-27 | 2007-10-03 | 讯宝科技公司 | Inertial drive scanning arrangement and method |
US20090099798A1 (en) * | 2007-10-09 | 2009-04-16 | Yanfeng Gong | Real-Time Power System Oscillation Detection Using Modal Analysis |
CN103928943A (en) * | 2014-04-21 | 2014-07-16 | 国家电网公司 | Research method for influence of doubly-fed wind turbine generator set on interconnected electric power system low frequency oscillation |
-
2014
- 2014-12-05 CN CN201410740740.9A patent/CN104392141A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030033094A1 (en) * | 2001-02-14 | 2003-02-13 | Huang Norden E. | Empirical mode decomposition for analyzing acoustical signals |
CN101048677A (en) * | 2004-10-27 | 2007-10-03 | 讯宝科技公司 | Inertial drive scanning arrangement and method |
US20090099798A1 (en) * | 2007-10-09 | 2009-04-16 | Yanfeng Gong | Real-Time Power System Oscillation Detection Using Modal Analysis |
CN103928943A (en) * | 2014-04-21 | 2014-07-16 | 国家电网公司 | Research method for influence of doubly-fed wind turbine generator set on interconnected electric power system low frequency oscillation |
Non-Patent Citations (2)
Title |
---|
潘秀才: "《基于HHT变换的电力系统低频振荡分析》", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
谭伟等: "《基于希尔伯特-黄变换(HHT)的机组同调研究》", 《南方电网技术》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104833878A (en) * | 2015-04-30 | 2015-08-12 | 国电南瑞科技股份有限公司 | Electric power system power oscillation judgment method based on synchronized phasor measurement unit |
CN104833878B (en) * | 2015-04-30 | 2017-12-01 | 国电南瑞科技股份有限公司 | Power oscillation of power system determination methods based on synchronous phasor measuring device |
CN106569086A (en) * | 2015-10-08 | 2017-04-19 | 国家电网公司 | Method for on-site positioning of low-frequency oscillation disturbance source based on oscillation phase difference |
CN112103939A (en) * | 2019-06-18 | 2020-12-18 | 上海电机学院 | Method and device for analyzing stability of power system |
CN112103939B (en) * | 2019-06-18 | 2024-01-26 | 上海电机学院 | Power system stability analysis method and device |
CN110492497A (en) * | 2019-08-16 | 2019-11-22 | 国网河南省电力公司电力科学研究院 | A kind of generating set power oscillation error comprehensive diagnosis method and system |
CN110492497B (en) * | 2019-08-16 | 2022-04-22 | 国网河南省电力公司电力科学研究院 | Power oscillation comprehensive diagnosis method and system for generator set |
CN112072650A (en) * | 2020-09-02 | 2020-12-11 | 深圳大学 | Power system power oscillation path identification method, system and storage medium |
CN115577226A (en) * | 2022-12-08 | 2023-01-06 | 南方电网数字电网研究院有限公司 | Forced oscillation source positioning method, device, equipment and storage medium |
CN115577226B (en) * | 2022-12-08 | 2023-03-03 | 南方电网数字电网研究院有限公司 | Forced oscillation source positioning method, device, equipment and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104392141A (en) | Method and device for locating low-frequency oscillation disturbance source of electric power system | |
Delgado-Arredondo et al. | Methodology for fault detection in induction motors via sound and vibration signals | |
Riera-Guasp et al. | Diagnosis of induction motor faults via Gabor analysis of the current in transient regime | |
Zanardelli et al. | Identification of intermittent electrical and mechanical faults in permanent-magnet AC drives based on time–frequency analysis | |
EP2790028B1 (en) | Broken rotor bar detection based on current signature analysis of an electric machine | |
Wang et al. | A new method with Hilbert transform and slip-SVD-based noise-suppression algorithm for noisy power quality monitoring | |
CN105547698A (en) | Fault diagnosis method and apparatus for rolling bearing | |
AU2012314960B2 (en) | A method of determining stationary signals for the diagnostics of an electromechanical system | |
KR20150080484A (en) | A Method for the Diagnostics of Electromechanical System Based on Impedance Analysis | |
Delgado-Arredondo et al. | Comparative study of time-frequency decomposition techniques for fault detection in induction motors using vibration analysis during startup transient | |
Climente-Alarcon et al. | Diagnosis of induction motors under varying speed operation by principal slot harmonic tracking | |
CN102353500B (en) | Extraction method of unbalanced signal for dynamic balance measurement | |
CN104462695A (en) | Weak signal detection method based on double-coupling Duffing vibrators and scale varying | |
Xie et al. | Fast-varying AM–FM components extraction based on an adaptive STFT | |
Chen et al. | Nonstationary signal denoising using an envelope-tracking filter | |
CN103913222A (en) | Duffing oscillator weak signal time domain detection method based on phase-locked loop | |
CN104931793A (en) | Substation grounding grid grounding resistance acquisition method | |
CN110907864B (en) | Fault detection method, device and equipment for motor stator winding and storage medium | |
CN105510066A (en) | Adaptive-noise-reduction-algorithm-based same-class rotary machinery system fault diagnosis method | |
CN106772193B (en) | Measuring method using current transformer frequency characteristic measuring device | |
JP6062588B1 (en) | A method for determining the characteristic features of a three-phase induction motor. | |
CN107861033B (en) | Calibration method and system for positioning error of oscillation wave partial discharge detection system | |
Vedreño-Santos et al. | Diagnosis of faults in induction generators under fluctuating load conditions through the instantaneous frequency of the fault components | |
CN108007548B (en) | Method for diagnosing equipment fault through frequency sweep | |
US10094743B2 (en) | Order analysis system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150304 |
|
RJ01 | Rejection of invention patent application after publication |