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 PDF

Info

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
Application number
CN201711086022.4A
Other languages
Chinese (zh)
Other versions
CN108180152B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201711086022.4A priority Critical patent/CN108180152B/en
Publication of CN108180152A publication Critical patent/CN108180152A/en
Application granted granted Critical
Publication of CN108180152B publication Critical patent/CN108180152B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04DNON-POSITIVE-DISPLACEMENT PUMPS
    • F04D17/00Radial-flow pumps, e.g. centrifugal pumps; Helico-centrifugal pumps
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04DNON-POSITIVE-DISPLACEMENT PUMPS
    • F04D27/00Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
    • F04D27/001Testing 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

A kind of wind turbine Weak fault detection method based on vibration signal cyclo-stationary
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.
CN201711086022.4A 2017-11-07 2017-11-07 A kind of blower Weak fault detection method based on vibration signal cyclo-stationary Active CN108180152B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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