CN108180152A - A kind of wind turbine Weak fault detection method based on vibration signal cyclo-stationary - Google Patents
A kind of wind turbine Weak fault detection method based on vibration signal cyclo-stationary Download PDFInfo
- Publication number
- CN108180152A CN108180152A CN201711086022.4A CN201711086022A CN108180152A CN 108180152 A CN108180152 A CN 108180152A CN 201711086022 A CN201711086022 A CN 201711086022A CN 108180152 A CN108180152 A CN 108180152A
- Authority
- CN
- China
- Prior art keywords
- wind turbine
- signal
- stationary
- data
- detection method
- 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
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D17/00—Radial-flow pumps, e.g. centrifugal pumps; Helico-centrifugal pumps
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
- F04D27/001—Testing thereof; Determination or simulation of flow characteristics; Stall or surge detection, e.g. condition monitoring
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
A kind of wind turbine Weak fault detection method based on vibration signal cyclo-stationary, includes the following steps:Step 1, acceleration signal is acquired;Step 2, collected signal is imported into program to be run, with the correlative character function check based on cyclostationary characteristic, obtains circulating density spectrum;Step 3, circulating density is composed and normalized, obtain cycle Correlated Spectroscopy;Step 4, it is strengthened envelope spectrum according to cycle Correlated Spectroscopy, preserves the ordinate data for strengthening envelope spectrum, as handle obtained data in real time;Step 5, using relevance function, calculating handles obtained data and the correlation of database data in real time;Step 6, according to relative coefficient failure resolution standard, judge whether this kind of failure be normal or belongs to which kind of failure, and will determine that result adds in database.It can detect and judge the fault type of wind turbine in real time using the present invention, detection is more accurate, has powerful practicability.
Description
Technical field
The invention belongs to signal processing and field of fault detection more particularly to one kind to be based on vibration signal cyclo-stationary
Wind turbine Weak fault detection method.
Background technology
Cyclo-stationary signal processing is a kind of emerging technology of mechanical signal processing risen recently.Cyclo-stationary signal is
The signal of hiding cycle information is included in signal.Cyclo-stationary signal is one kind of non-stationary signal, is examined compared to tradition
Survey mode, the signal generated closer to actual signal, especially rotating machinery.
The common rotating machinery fault detection method of field of signal processing mainly has Fourier transform, Short-time Fourier at present
Transformation, wavelet transformation, Second Generation Wavelet Transformation and multi-wavelet transformation etc., it may be said that be all based on the signature waveform base of inner product principle
Function signal decomposes, it is intended to flexibly with going preferably to handle signal with the basic function that signature waveform matches, extraction failure spy
Sign, so as to fulfill fault diagnosis.
But following shortcoming and defect exists in the prior art:Fourier transform, short time discrete Fourier transform, wavelet transformation,
The method of the fault detects such as Second Generation Wavelet Transformation and multi-wavelet transformation is built upon assume that signal is stationary signal on the basis of,
And in reality it is often non-stationary signal, so as to which these detection methods have unreasonable place, do not conform to reality.Meanwhile these
Traditional detection method is due to theoretic limitation, it is difficult to detect some important features of rotating machinery, such as blade passing frequency
BPF, blade have significant limitation than frequency BRF etc..
Invention content
The present invention provides a kind of wind turbine Weak fault detection methods based on vibration signal cyclo-stationary, can be real-time
Ground detects and judges the fault type of wind turbine, and detection is more accurate, has powerful practicability.
A kind of wind turbine Weak fault detection method based on vibration signal cyclo-stationary, includes the following steps:
Step 1, using the vibration acceleration signal of acceleration transducer acquisition wind turbine;
Step 2, collected acceleration signal is imported into program to be run, with the phase based on cyclostationary characteristic
The characteristic function detection of closing property, obtains circulating density spectrum;
Step 3, circulating density is composed and normalized, obtain cycle Correlated Spectroscopy;
Step 4, it is strengthened envelope spectrum according to cycle Correlated Spectroscopy, preserves the ordinate data for strengthening envelope spectrum, as in fact
When the data that handle;
Step 5, using relevance function, calculating handles obtained data and the correlation of database data in real time;
Step 6, according to relative coefficient failure resolution standard, judge whether this kind of failure be normal or belongs to which kind of failure,
And it will determine that result adds in database.
In step 2, the correlative character function check process based on cyclostationary characteristic is:
Relevant parameter is set in a program, calculates the parameter of cyclo-stationary:
Wherein:T, T is the time;△ f go to zero, and the sequence that T tends to be just infinite cannot exchange;f1、f2Represent two calculated
Frequency;xΔf(t;f1) represent filtering;Represent xΔf(t,f2) conjugate complex number;J represents imaginary unit.
In step 4, the process for the envelope spectrum that strengthened is:The three-dimensional matrice of cycle Correlated Spectroscopy will be obtained, withTable
Show, then keeping cycle frequency α constant, frequency f being integrated, strengthened envelope spectrum after normalizing.
In step 5, the relevance function is
Wherein, A is to handle obtained data in real time, and B is database data;Cov (A, B) is contrast signal and fitted signal
Covariance, calculation formula is as follows:
Cov (A, B)=E (A × B)-E (A) * E (B)
Wherein, E (A × B) is the mathematic expectaion of the product of A and B;E (A) is the mathematic expectaion of A;E (B) is the mathematic expectaion of B.
In step 5, the process for obtaining database data is:Under the conditions of wind turbine normal work and various typical faults, with
Acceleration transducer acquires the acceleration signal of fan vibration, and inputs in line program to be shipped, normalized by what is be calculated
The ordinate output for strengthening envelope spectrum is data file, is stored in Database Folder.
The present invention solves the problems, such as to be assumed to be stationary signal to a certain extent, since wind turbine is as rotating machinery, production
Raw signal is largely cyclo-stationary signal, and the result of detection is closer to practical, more reliable.Meanwhile cyclo-stationary
Handling result overcome traditional detection method, blade passing frequency and blade ratio frequency detecting less than or unconspicuous difficulty, energy
The more features of wind turbine is showed.
In detection process, the time of program operation is short, can achieve the purpose that monitor in real time, according to relative coefficient and one
Fixed experience judges which kind of failure wind turbine belongs to, and by establishing the file of normal condition and fault condition, needs to compare to enrich
Compared with database so that later detection judgement more accurate, even entirely without experience is the detection of wind turbine Weak fault
Provide a kind of good mode.
Description of the drawings
Fig. 1 is the flow chart of the wind turbine Weak fault detection method the present invention is based on vibration signal cyclo-stationary;
Fig. 2 is the reinforcement envelope spectrogram of ideal signal simulation result;
Fig. 3 is the Fast Fourier Transform figure of ideal signal simulation result;
Fig. 4 is the reinforcement envelope diagram after normal wind turbine data processing;
Fig. 5 is the Fast Fourier Transform figure after normal wind turbine data processing;
Fig. 6 is the reinforcement envelope diagram after the wind turbine data processing that a foundation bolt loosens;
Fig. 7 is the Fast Fourier Transform figure after the wind turbine data processing that a foundation bolt loosens.
Specific embodiment
In order to more specifically describe the present invention, below in conjunction with the accompanying drawings and specific embodiment is to technical scheme of the present invention
It is described in detail.
As shown in Figure 1, the wind turbine Weak fault detection method based on vibration signal cyclo-stationary, includes the following steps:
S01, with the vibration acceleration signal of acceleration transducer acquisition wind turbine.
S02 sets corresponding parameter in a program, and collected signal is imported into program, calculates cyclo-stationary
Parameter:
Wherein:T, T is the time;△ f go to zero, and the sequence that T tends to be just infinite cannot exchange;f1、f2Represent two calculated
Frequency;xΔf(t,f1) represent filtering;Represent xΔf(t,f2) conjugate complex number;J represents imaginary unit.
Fan vibration signal can be reduced to:
Wherein v (t) is random signal;AiFor real number, the mould of cosine signal is represented;αiFor the cycle frequency for needing to detect.
For random signal v (t), work as corrv(f1,f2) in f1=f2When corr values be not zero.
Enable α=f1-f2, whenOrOr ± αiWhen, above formula corrx(f1,f2) value be not zero, that is, examine
The cycle frequency for needing to detect is measured, by corr functions, we can detect existing cyclo-stationary information in signal.
For emulation signal x=[1+cos (60 π t)+cos (100 π t)+cos (130 π t)] × N (0,1) due to negative value portion
The information divided can not be shown on the image, then at α=15,20,30,35,50,60,65,80,95,100,115,130, cycle
There is larger peak value, but the cyclical information that image at this time can't be indicated a need for intuitively in density spectra.
S03, obtained amplitude differ greatly, and are 1 by setting maximum value, obtained circulating density is composed and is normalized, is obtained
It is intuitive to show cycle frequency information to cycle Correlated Spectroscopy image;
S04 obtains the three-dimensional matrice of cycle Correlated Spectroscopy, withIt represents, then keeps cycle frequency α constant, f is accumulated
Point, strengthened envelope spectrum s (α) after normalizing:
S (α) is proportional to
Obtained reinforcement envelope spectrum ordinate matrix is saved as data file by S05;
The data file handled in real time is done relative coefficient calculating by S06 with corresponding normal and fault data file:
Wherein, A is the data in the data file handled in real time, and B is that the correspondence preserved in database is normal or failure
Data in file;Cov (A, B) is the covariance of contrast signal and fitted signal, and calculation formula is as follows:
Cov (A, B)=E (A × B)-E (A) * E (B)
Wherein, E (A × B) is the mathematic expectaion of the product of A and B;E (A) is the mathematic expectaion of A;E (B) is the mathematic expectaion of B.
S07 is detected and is judged with the related coefficient breakdown judge standard and experience tentatively established, real-time condition whether failure
And fault type, and the result that will determine that is saved in data, and database is helped to improve the mark of relative coefficient breakdown judge
It is accurate.
In order to which specific manifestation this method is in the advantage and feature of wind turbine Weak fault process field, to ideal signal x=[1+
Cos (60 π t)+cos (100 π t)+cos (130 π t)] × N (0,1) emulated, wherein, N (0,1) represents width in the time domain
Degree obedience mean value is zero, the random signal that variance is 1, and also the vibration acceleration data tested on wind turbine are carried out
Processing and comparison, while have also been made respective Fast Fourier Transform.
The simulation result of ideal signal strengthens envelope spectrum image and Fast Fourier Transform image, as shown in Figures 2 and 3;
The handling result of wind turbine data that the normal wind turbine data of acquisition and a foundation bolt loosen, the reinforcement envelope spectrum finally obtained
Image and Fast Fourier Transform image.Reinforcement envelope spectrum image and Fast Fourier Transform figure after normal wind turbine data processing
As shown in Figure 4 and Figure 5, reinforcement envelope spectrum image and fast Flourier after the wind turbine data processing that a foundation bolt loosens
Changing image is as shown in Figure 6 and Figure 7.
It is not difficult to find out, the frequency and Amplitude Ration that simulation result medium wave peak occurs all meet theory deduction, simultaneously because
The error amount estimated and occurred is also visibly homogeneous in envelope spectrum is strengthened, and demonstrates the feasible part of our programs;Normal wind turbine
Real data processing, obtained image meet actual treatment as a result, detected basic frequency 10Hz and in wind turbine
Due to the frequency multiplication of fundamental frequency that blade rotation actually generates, 20Hz, 30Hz, 40Hz etc..And traditional Fast Fourier Transform, although
Basic frequency 10Hz is equally detected, but certain frequencys multiplication, such as 20Hz, 30Hz, detection result are bad.Strengthen envelope stave
More information is showed and has been easy to differentiate and read, while the datum line of fluctuation very little for strengthening occurring in envelope spectrum is also demonstrated by this
The dependable with function that method handles real data.Meanwhile the processing knot of the Weak fault loosened to a foundation bolt
Fruit also detected basic frequency 10Hz and in wind turbine due to the frequency multiplication of fundamental frequency that actually generates of blade rotation, 20Hz,
30Hz, 40Hz etc. have different in frequency distribution and amplitude size with reinforcement envelope spectrum under normal circumstances, calculate related
Property coefficient is about 0.82, it was demonstrated that the method is used for the practicability and reliability of detection blower fan Weak fault.
Technical scheme of the present invention and advantageous effect is described in detail in above-described specific embodiment, Ying Li
Solution is the foregoing is merely presently most preferred embodiment of the invention, is not intended to restrict the invention, all principle models in the present invention
Interior done any modification, supplementary, and equivalent replacement etc. are enclosed, should all be included in the protection scope of the present invention.
Claims (5)
1. a kind of wind turbine Weak fault detection method based on vibration signal cyclo-stationary, includes the following steps:
Step 1, using the vibration acceleration signal of acceleration transducer acquisition wind turbine;
Step 2, collected acceleration signal is imported into program to be run, with the correlation based on cyclostationary characteristic
Characteristic function detects, and obtains circulating density spectrum;
Step 3, circulating density program calculated composes normalization, obtains cycle Correlated Spectroscopy;
Step 4, it is strengthened envelope spectrum according to cycle Correlated Spectroscopy, preserves the ordinate data for strengthening envelope spectrum, as place in real time
Manage obtained data;
Step 5, using relevance function, calculating handles obtained data and the correlation of database data in real time;
Step 6, according to relative coefficient failure resolution standard, judge whether this kind of failure be normal or belongs to which kind of failure, and will
Judging result adds in database.
2. the wind turbine Weak fault detection method according to claim 1 based on vibration signal cyclo-stationary, feature
It is, in step 2, the correlative character function check method based on cyclostationary characteristic is:
Wherein:T, T is the time;△ f go to zero, and the sequence that T tends to be just infinite cannot exchange;f1、f2Represent two frequencies calculated
Rate;xΔf(t,f1) represent filtering;Represent xΔf(t,f2) conjugate complex number;J represents imaginary unit.
3. the wind turbine Weak fault detection method according to claim 1 based on vibration signal cyclo-stationary, feature
It is, in step 4, the process for the envelope spectrum that strengthened is:The three-dimensional matrice of cycle Correlated Spectroscopy will be obtained, withIt represents, so
It keeps cycle frequency α constant afterwards, frequency f is integrated, strengthened envelope spectrum after normalizing.
4. the wind turbine Weak fault detection method according to claim 1 based on vibration signal cyclo-stationary, feature
It is, in step 5, the relevance function is
Wherein, A is to handle obtained data in real time, and B is database data;Cov (A, B) is the association of contrast signal and fitted signal
Variance, calculation formula are as follows:
Cov (A, B)=E (A × B)-E (A) * E (B)
Wherein, E (A × B) is the mathematic expectaion of the product of A and B;E (A) is the mathematic expectaion of A;E (B) is the mathematic expectaion of B.
5. the wind turbine Weak fault detection method according to claim 1 based on vibration signal cyclo-stationary, feature
It is, in step 5, the process for obtaining database data is:Under the conditions of wind turbine normal work and various typical faults, to accelerate
The acceleration signal of sensor acquisition fan vibration is spent, and is inputted in line program to be shipped, the normalized reinforcement that will be calculated
The ordinate output of envelope spectrum is data file, is stored in Database Folder.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711086022.4A CN108180152B (en) | 2017-11-07 | 2017-11-07 | A kind of blower Weak fault detection method based on vibration signal cyclo-stationary |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711086022.4A CN108180152B (en) | 2017-11-07 | 2017-11-07 | A kind of blower Weak fault detection method based on vibration signal cyclo-stationary |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108180152A true CN108180152A (en) | 2018-06-19 |
CN108180152B CN108180152B (en) | 2019-05-14 |
Family
ID=62544933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711086022.4A Active CN108180152B (en) | 2017-11-07 | 2017-11-07 | A kind of blower Weak fault detection method based on vibration signal cyclo-stationary |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108180152B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085763A (en) * | 2018-06-11 | 2018-12-25 | 浙江上风高科专风实业有限公司 | Extraction of the cyclo-stationary method based on complicated amplitude modulation model to fan vibration feature |
CN109323757A (en) * | 2018-10-29 | 2019-02-12 | 浙江大学 | A method of estimation bubble population is to propeller sound source characteristics frequency inhibiting effect |
CN109341780A (en) * | 2018-11-29 | 2019-02-15 | 浙江省环境保护科学设计研究院 | A kind of more means low cost fan trouble monitoring methods |
CN110320018A (en) * | 2019-07-12 | 2019-10-11 | 北京交通大学 | A kind of combined failure of rotating machinery diagnostic method based on second-order cyclostationary characteristic |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1494476A (en) * | 1973-11-13 | 1977-12-07 | Imp Metal Ind Kynoch Ltd | Fourier analysers |
CN1514209A (en) * | 2003-08-01 | 2004-07-21 | 重庆大学 | Rotary machine failure intelligent diagnosis method and device |
CN104268335A (en) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | Vibration fatigue life predication method and system for micro-packaging assembly |
CN105758604A (en) * | 2015-11-18 | 2016-07-13 | 上海航天精密机械研究所 | Multiaxis vibration test condition cutting method based on acceleration equivalence |
CN106092565A (en) * | 2016-06-08 | 2016-11-09 | 成都阜特科技股份有限公司 | A kind of Vibration Analysis method and system thereof |
CN106482827A (en) * | 2016-11-11 | 2017-03-08 | 北京航空航天大学 | Electronic product based on crosspower spectrum function Modal Parameter Identification vibrates DLP method |
-
2017
- 2017-11-07 CN CN201711086022.4A patent/CN108180152B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1494476A (en) * | 1973-11-13 | 1977-12-07 | Imp Metal Ind Kynoch Ltd | Fourier analysers |
CN1514209A (en) * | 2003-08-01 | 2004-07-21 | 重庆大学 | Rotary machine failure intelligent diagnosis method and device |
CN104268335A (en) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | Vibration fatigue life predication method and system for micro-packaging assembly |
CN105758604A (en) * | 2015-11-18 | 2016-07-13 | 上海航天精密机械研究所 | Multiaxis vibration test condition cutting method based on acceleration equivalence |
CN106092565A (en) * | 2016-06-08 | 2016-11-09 | 成都阜特科技股份有限公司 | A kind of Vibration Analysis method and system thereof |
CN106482827A (en) * | 2016-11-11 | 2017-03-08 | 北京航空航天大学 | Electronic product based on crosspower spectrum function Modal Parameter Identification vibrates DLP method |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085763A (en) * | 2018-06-11 | 2018-12-25 | 浙江上风高科专风实业有限公司 | Extraction of the cyclo-stationary method based on complicated amplitude modulation model to fan vibration feature |
CN109085763B (en) * | 2018-06-11 | 2020-02-21 | 浙江上风高科专风实业有限公司 | Extraction of fan vibration characteristics by using circulation stabilization method based on complex amplitude modulation model |
CN109323757A (en) * | 2018-10-29 | 2019-02-12 | 浙江大学 | A method of estimation bubble population is to propeller sound source characteristics frequency inhibiting effect |
CN109341780A (en) * | 2018-11-29 | 2019-02-15 | 浙江省环境保护科学设计研究院 | A kind of more means low cost fan trouble monitoring methods |
CN110320018A (en) * | 2019-07-12 | 2019-10-11 | 北京交通大学 | A kind of combined failure of rotating machinery diagnostic method based on second-order cyclostationary characteristic |
CN110320018B (en) * | 2019-07-12 | 2020-08-11 | 北京交通大学 | Rotary machine composite fault diagnosis method based on second-order cyclostationarity |
Also Published As
Publication number | Publication date |
---|---|
CN108180152B (en) | 2019-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108180152B (en) | A kind of blower Weak fault detection method based on vibration signal cyclo-stationary | |
Chen et al. | Detection of rub-impact fault for rotor-stator systems: A novel method based on adaptive chirp mode decomposition | |
Li et al. | Application of bandwidth EMD and adaptive multiscale morphology analysis for incipient fault diagnosis of rolling bearings | |
Borghesani et al. | CS2 analysis in presence of non-Gaussian background noise–Effect on traditional estimators and resilience of log-envelope indicators | |
Cui et al. | Quantitative trend fault diagnosis of a rolling bearing based on Sparsogram and Lempel-Ziv | |
CN103575523B (en) | The rotary machinery fault diagnosis method of kurtosis-envelope spectrum analysis is composed based on FastICA- | |
Wu et al. | Characterization of gear faults in variable rotating speed using Hilbert-Huang Transform and instantaneous dimensionless frequency normalization | |
Urbanek et al. | Application of averaged instantaneous power spectrum for diagnostics of machinery operating under non-stationary operational conditions | |
Yongbo et al. | Review of local mean decomposition and its application in fault diagnosis of rotating machinery | |
CN107631877A (en) | A kind of rolling bearing fault collaborative diagnosis method for casing vibration signal | |
CN108072517B (en) | Method for detecting weak fault signal of rotating machine | |
CN109375060B (en) | Method for calculating fault waveform similarity of power distribution network | |
Xue et al. | Automatic diagnosis method for structural fault of rotating machinery based on distinctive frequency components and support vector machines under varied operating conditions | |
Cong et al. | Research of singular value decomposition based on slip matrix for rolling bearing fault diagnosis | |
CN101603854A (en) | The rotating machinery non-stationery vibration signal instantaneous frequency estimation algorithm in start and stop period | |
Li et al. | Incipient fault information determination for rolling element bearing based on synchronous averaging reassigned wavelet scalogram | |
Wu et al. | A carrier wave extraction method for cavitation characterization based on time synchronous average and time-frequency analysis | |
Han et al. | Multistate fault diagnosis strategy for bearings based on an improved convolutional sparse coding with priori periodic filter group | |
JP5780870B2 (en) | Rotating equipment soundness diagnosis apparatus, method and program | |
Zhao et al. | Vibration health monitoring of rolling bearings under variable speed conditions by novel demodulation technique | |
CN103471708B (en) | Rotating machine fault diagnosis method based on nonlinear ICA (Independent Component Analysis) of improved particle swarm | |
Choudhury et al. | A methodology to handle spectral smearing in gearboxes using adaptive mode decomposition and dynamic time warping | |
Abolghasemi et al. | Recursive singular spectrum analysis for induction machines unbalanced rotor fault diagnosis | |
CN107559228B (en) | Method based on bispectral data detection and diagnosis fan trouble | |
Hosseini et al. | Comparison of nonlinearity measures based on time series analysis for nonlinearity detection |
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 |