CN106770015A - A kind of oil property detection method based on the similar differentiation of principal component analysis - Google Patents

A kind of oil property detection method based on the similar differentiation of principal component analysis Download PDF

Info

Publication number
CN106770015A
CN106770015A CN201710017913.8A CN201710017913A CN106770015A CN 106770015 A CN106770015 A CN 106770015A CN 201710017913 A CN201710017913 A CN 201710017913A CN 106770015 A CN106770015 A CN 106770015A
Authority
CN
China
Prior art keywords
sample
principal component
component analysis
similar
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
CN201710017913.8A
Other languages
Chinese (zh)
Other versions
CN106770015B (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.)
NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd
Original Assignee
NANJING RICHISLAND INFORMATION ENGINEERING 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 NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd filed Critical NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd
Priority to CN201710017913.8A priority Critical patent/CN106770015B/en
Publication of CN106770015A publication Critical patent/CN106770015A/en
Application granted granted Critical
Publication of CN106770015B publication Critical patent/CN106770015B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3577Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The present invention proposes a kind of oil property detection method based on the similar differentiation of principal component analysis, after the near infrared spectrum to sample carries out conventional pretreatment, the score matrix of spectroscopic data is obtained by principal component analysis, first three score vector for choosing score matrix constitutes new score matrix, according to calibration set Sample Storehouse and the score vector drawing three-dimensional principal component analysis figure of sample to be tested, length and width at high proportion 3 are drawn centered on sample to be tested in figure:2:1 cube frame so that number of samples in cube frame is 50 ± 5, and using these samples as sample to be tested similar sample, then the spectroscopic data according to similar sample set up partial least square model, finally treat survey sample properties and be predicted.The method energy is quick, Accurate Prediction oil property, is favorably improved Business Economic Benefit.

Description

A kind of oil property detection method based on the similar differentiation of principal component analysis
Technical field
The present invention is a kind of oil property method for quick, specifically a kind of oil based on the similar differentiation of principal component analysis Moral character quality detection method.
Background technology
At present, near-infrared spectral analytical method is widely used in oil property analysis, with traditional laboratory assay Method is compared, and the method has the advantages that fast analyze speed, high precision and expends few.
Modeling method based on fractional sample, is a kind of effective ways that can improve model accuracy, its basic thought It is:Chosen and one group of most like sample of sample to be tested from calibration set sample based on spectrum, it is then (i.e. local by these samples Sample) obtain final predicting the outcome by statistical analysis or the bearing calibration of classics.Modeling strategy based on fractional sample is fitted For the correction of Nonlinear system, while the advantage of Sample Storehouse can be made full use of, it is to avoid traditional factor-analysis approach is because of sample sets Into the disadvantage for waiting variation to need frequent updating model.
Our early stages have pointed out a kind of octane value detection method based on similar differentiation, and the method is pre-processed to spectrum Afterwards, the score matrix of absorbance matrix is calculated using principal component analysis (PCA) method, using accumulation contribution rate 85%-95% Corresponding score vector constitutes new score matrix, and sample to be tested and spectrum number are calculated using Euclidean distance based on new score matrix According to the spectrum intervals of sample in storehouse, and as the criterion of similar sample (i.e. fractional sample) is searched, select light spectrum distance Partial least square model is set up as calibration samples from the similar sample less than threshold value, and sample to be tested is predicted.
The above method typically chooses the first two score vector of score matrix, calculates to be measured equivalent in two-dimensional coordinate figure The spectrum intervals of sample in sample and spectra database.Inevitably, the first principal component and Second principal component, of some samples It is very close, but the 3rd principal component has difference, and this causes the position distribution illusion in two-dimensional principal component analysis figure, in turn results in phase Choose incorrect like sample, the precision for predicting the outcome is not ideal enough.
The content of the invention
In order to solve the above problems, the present invention proposes a kind of oil property detection based on the similar differentiation of principal component analysis Method.
The present invention is comprised the following steps that:
(1) near infrared spectrum of oil product sample to be measured is obtained;
(2) conventional pretreatment is carried out to sample spectrum in oil product spectrum to be measured and calibration set;
(3) pretreated all spectroscopic datas are carried out into principal component analysis;
(4) according to the score matrix after principal component analysis, m score vector before choosing draws principal component analysis figure;
(5) according to principal component analysis figure, p similar sample of sample to be tested is chosen according to certain rule;
(6) model is set up using offset minimum binary according to similar sample;
(7) model by building up is predicted to the property of sample to be tested.
Preprocessing procedures are using baseline correction, wave band interception and vector normalizing.
On the premise of the principal character for being extracted spectrum, this method is wished accurate and chooses most like with sample to be tested facing Nearly sample, therefore this method depicts three-dimensional principal component analysis figure, i.e., with first principal component as transverse axis, Second principal component, is vertical Axle, the 3rd principal component is vertical pivot.
Because first principal component represents the maximum direction of absorbance matrix variation, Second principal component, takes second place, the 3rd principal component Third, therefore this method draws cube frame in principal component analysis figure, the length and width of cube frame are at high proportion 3:2:1, will Positioned at cube inframe sample as sample to be tested similar sample.
This method chooses 50 similar samples, i.e. p=50 ± 5.
Beneficial effect:
Detection method provided by the present invention is based on oil product near infrared spectrum, using principal component analysis combination offset minimum binary Method is modeled, and realizes the quick detection of oil property.Compared with general modeling method, this method can quickly and more Accurate Prediction The octane number of oil property, such as gasoline product, density, saturated vapour pressure, contribute to the vehicle air-conditioning of Petrochemical Enterprises, And then improve the economic benefit of enterprise.
Brief description of the drawings
Fig. 1 is based on the oil property detection method flow chart of the similar differentiation of principal component analysis
Distribution of Fig. 2 octane numbers to be measured in principal component analysis figure
Specific embodiment
The present invention is further illustrated with case study on implementation below in conjunction with the accompanying drawings.
The present invention introduces the oil property detection method based on the similar differentiation of principal component analysis by taking certain 92# product oil as an example. Table 1 is the numbering and its corresponding research octane number (RON) (RON) of certain all sample of 92# product oils.
Certain the 92# product oils sample number of table 1 and corresponding RON
In table 1, the sample of numbering 95#-1~290 is used as calibration set sample, the sample conduct of numbering 95#-291~300 Sample to be tested.Below by taking numbering 95#-291 samples to be tested as an example, the detailed process of present invention prediction octane number is illustrated:
The first step:Sample to numbering 95#-1~291 carries out Pretreated spectra, including baseline correction, spectrogram interception and arrow Amount normalization.
First, baseline correction uses two-point method baseline correction, chooses 6400cm-1And 9200cm-1Two wave numbers o'clock are used as two bases Point, the absorbance after baseline correction is calculated by following formula:
In formula, xiIt is gasoline in the wave number of near infrared spectrum;kxi+ b was 6400cm-1And 9200cm-12 points of straight line Equation, wherein k are the straight slope, and b is the Linear intercept;yiRepresent former spectrogram in wave number xiUnder absorbance;yi *Represent base Spectrogram is in wave number x after line correctioniUnder absorbance.
Secondly, 4000cm is intercepted-1~4800cm-1Spectrogram in wave number section.
Finally, vector normalization is calculated using following formula:
In formula, XijRefer to absorbance of i-th sample under wave number j;Refer to i-th absorbance values of sample;m It is the number of wave number point;Xij *Represent absorbance of i-th sample after vector normalization under wave number j.
Second step:After carrying out above-mentioned pretreatment to the spectral data of the gasoline sample of numbering 95#-1~291, by it is main into Divide analysis to obtain score matrix, and choose preceding 3 score vectors, as shown in table 2.
The numbering of table 2 is the score vector of the spectrum of the gasoline sample of 95#-1~291
Fig. 2 be the gasoline sample to be measured of numbering 92#-291 in the position of PCA distribution maps, the triangle mark seen in Fig. 2, most Whole cube frame has also shown in figure.
3rd step:The sample of cube inframe as the similar sample of 92#-291 will be located at.The numbering of these similar samples And property laboratory values are as shown in table 3.
The similar sample number and RON of the gasoline sample of the numbering 92#-291 of table 3
4th step:The absorbance data of the similar sample in table 3 sets up model using offset minimum binary, then right The RON of 92#-291 samples is predicted, and it is 92.16 to predict the outcome, and octane number laboratory values are 92.2, and predicated error is only 0.04.
In order to prediction effect of the invention is better described, to above-mentioned 92# product oils Sample Storehouse and sample to be tested, with me Early stage application a kind of side that is proposed of patent " octane value detection method based on similar differentiation " (A of CN 104990893) Method is contrasted.Table 4 gives predicting the outcome for two methods and compares.
Predicting the outcome for the two methods of table 4 is compared
Appreciation gist root-mean-square error for final result of the invention, i.e. RMSE.RMSE is calculated by following formula:
In formula, np is the number of gasoline sample to be measured;Refer to i-th octane number predicted value of gasoline to be measured;xiRefer to i-th The octane number laboratory values of gasoline to be measured.RMSE value is smaller, illustrates that the accuracy of prediction is higher, and prediction effect is better.
It is computed, the RMSE of new method is 0.2086, and former method is 0.2232.It can be seen that, although former method is more conventional Method for quick has more excellent performance, and the method for the invention precision on the basis of former method has and further carries Rise.
Although technical scheme disclosed in the A of CN 104990893 is generally to choose sample in two-dimensional coordinate figure, Seem the technical inspiration for giving and sample selection being carried out with three-dimensional coordinate figure.But be emphasized that the A institutes of CN 104990893 Selection the first two score vector carries out sample selection and is drawn in 85%~95% by calculating accumulation contribution rate, and non-straight Selecting take (if exist the first two score vector accumulation contribution rate be unsatisfactory for requirement extend to first three, four ... individual score vectors May);And the scheme of the application is clearly to carry out solid space with first three score vector to divide to choose sample.The two thinking Entirely different, can put together carries out Contrast on effect, but to each other can't generation technology enlightenment.

Claims (7)

1. a kind of oil property detection method based on the similar differentiation of principal component analysis, it is characterised in that have steps of:
(1) near infrared spectrum of oil product sample to be measured is obtained;
(2) conventional pretreatment is carried out to sample spectrum in oil product spectrum to be measured and calibration set;
(3) pretreated all spectroscopic datas are carried out into principal component analysis;
(4) according to the score matrix after principal component analysis, m score vector before choosing draws principal component analysis figure;
(5) according to principal component analysis figure, p similar sample of sample to be tested is chosen according to certain rule;
(6) model is set up using offset minimum binary according to similar sample;
(7) model by building up is predicted to the property of sample to be tested.
2. a kind of oil property detection method based on the similar differentiation of principal component analysis according to claim 1, its feature Be m=3, i.e., with first principal component as transverse axis, Second principal component, is the longitudinal axis, the 3rd principal component is vertical pivot, drawing three-dimensional it is main into Divide analysis chart.
3. a kind of oil property detection method based on the similar differentiation of principal component analysis according to claim 2, its feature It is that certain rule is:Centered on sample to be tested, cube frame is drawn in principal component analysis figure, positioned at cube frame Interior sample as sample to be tested similar sample.
4. a kind of oil property detection method based on the similar differentiation of principal component analysis according to claim 3, its feature The length and width for being cube frame are at high proportion 3:2:1.
5. a kind of oil property detection method based on the similar differentiation of principal component analysis according to claim 1, its feature It is p=50 ± 5.
6. a kind of oil property detection method based on the similar differentiation of principal component analysis according to claim 1, its feature It is that the conventional pretreatment includes baseline correction, spectrogram interception and vector normalization.
7. a kind of oil property detection method based on the similar differentiation of principal component analysis according to claim 6, its feature It is the conventional pretreatment,
First:Baseline correction uses two-point method baseline correction, chooses 6400cm-1And 9200cm-1Two wave numbers o'clock are used as two basic points, base Absorbance after line correction is calculated by following formula:
y i * = y i - ( kx i + b ) - - - ( 1 )
In formula, xiIt is gasoline in the wave number of near infrared spectrum;kxi+ b was 6400cm-1And 9200cm-12 points of straight line side Journey, wherein k are the straight slope, and b is the Linear intercept;yiRepresent former spectrogram in wave number xiUnder absorbance;Represent baseline school Spectrogram is in wave number x after justiUnder absorbance;
Secondly, 4000cm is intercepted-1~4800cm-1Spectrogram in wave number section;
Finally, vector normalization is calculated using following formula:
X i j * = X i j - X ‾ i Σ j = 1 m ( X i j - X ‾ i ) 2 - - - ( 2 )
In formula, XijRefer to absorbance of i-th sample under wave number j;Refer to i-th absorbance values of sample;M is wave number The number of point;Xij *Represent absorbance of i-th sample after vector normalization under wave number j.
CN201710017913.8A 2017-01-10 2017-01-10 Oil product property detection method based on principal component analysis similarity discrimination Active CN106770015B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710017913.8A CN106770015B (en) 2017-01-10 2017-01-10 Oil product property detection method based on principal component analysis similarity discrimination

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710017913.8A CN106770015B (en) 2017-01-10 2017-01-10 Oil product property detection method based on principal component analysis similarity discrimination

Publications (2)

Publication Number Publication Date
CN106770015A true CN106770015A (en) 2017-05-31
CN106770015B CN106770015B (en) 2020-12-15

Family

ID=58948996

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710017913.8A Active CN106770015B (en) 2017-01-10 2017-01-10 Oil product property detection method based on principal component analysis similarity discrimination

Country Status (1)

Country Link
CN (1) CN106770015B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107179293A (en) * 2017-06-23 2017-09-19 南京富岛信息工程有限公司 A kind of assessment method of oil property uncertainty
CN107271400A (en) * 2017-06-23 2017-10-20 南京富岛信息工程有限公司 A kind of method of automatic addition calibration set sample
CN107356535A (en) * 2017-06-12 2017-11-17 湖北久之洋红外系统股份有限公司 A kind of marine oil overflow detection method based on spectral imaging technology
CN107505282A (en) * 2017-08-28 2017-12-22 南京富岛信息工程有限公司 A kind of method for improving oil product near-infrared modeling robustness
CN108226093A (en) * 2018-01-11 2018-06-29 南京富岛信息工程有限公司 A kind of atmospheric and vacuum distillation unit model parameter automatically selects and bearing calibration
CN111474134A (en) * 2020-04-24 2020-07-31 驻马店华中正大有限公司 Method for controlling butyric acid fermentation by using online near infrared
CN113433088A (en) * 2021-06-25 2021-09-24 南京富岛信息工程有限公司 Fine monitoring method for oil mixing section of crude oil long-distance pipeline
CN115060687A (en) * 2022-08-18 2022-09-16 南京富岛信息工程有限公司 Tax administration method for finished oil production enterprise

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140197316A1 (en) * 2012-11-30 2014-07-17 Suncor Energy Inc. Measurement and control of bitumen-containing process streams
CN104990893A (en) * 2015-06-24 2015-10-21 南京富岛信息工程有限公司 Gasoline octane number detecting method based on similar discriminance
CN105891141A (en) * 2016-03-30 2016-08-24 南京富岛信息工程有限公司 Method for rapidly measuring gasoline property data
CN105954223A (en) * 2016-04-28 2016-09-21 南京富岛信息工程有限公司 Method for improving prediction accuracy of gasoline properties

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140197316A1 (en) * 2012-11-30 2014-07-17 Suncor Energy Inc. Measurement and control of bitumen-containing process streams
CN104990893A (en) * 2015-06-24 2015-10-21 南京富岛信息工程有限公司 Gasoline octane number detecting method based on similar discriminance
CN105891141A (en) * 2016-03-30 2016-08-24 南京富岛信息工程有限公司 Method for rapidly measuring gasoline property data
CN105954223A (en) * 2016-04-28 2016-09-21 南京富岛信息工程有限公司 Method for improving prediction accuracy of gasoline properties

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈颂超 等: "基于局部加权回归的土壤全氮含量可见_近红外光谱反演", 《土壤学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107356535A (en) * 2017-06-12 2017-11-17 湖北久之洋红外系统股份有限公司 A kind of marine oil overflow detection method based on spectral imaging technology
CN107179293A (en) * 2017-06-23 2017-09-19 南京富岛信息工程有限公司 A kind of assessment method of oil property uncertainty
CN107271400A (en) * 2017-06-23 2017-10-20 南京富岛信息工程有限公司 A kind of method of automatic addition calibration set sample
CN107505282A (en) * 2017-08-28 2017-12-22 南京富岛信息工程有限公司 A kind of method for improving oil product near-infrared modeling robustness
CN108226093A (en) * 2018-01-11 2018-06-29 南京富岛信息工程有限公司 A kind of atmospheric and vacuum distillation unit model parameter automatically selects and bearing calibration
CN111474134A (en) * 2020-04-24 2020-07-31 驻马店华中正大有限公司 Method for controlling butyric acid fermentation by using online near infrared
CN113433088A (en) * 2021-06-25 2021-09-24 南京富岛信息工程有限公司 Fine monitoring method for oil mixing section of crude oil long-distance pipeline
CN115060687A (en) * 2022-08-18 2022-09-16 南京富岛信息工程有限公司 Tax administration method for finished oil production enterprise
CN115060687B (en) * 2022-08-18 2022-11-08 南京富岛信息工程有限公司 Tax administration method for finished oil production enterprise

Also Published As

Publication number Publication date
CN106770015B (en) 2020-12-15

Similar Documents

Publication Publication Date Title
CN106770015A (en) A kind of oil property detection method based on the similar differentiation of principal component analysis
CN104990894B (en) A kind of gasoline property detection method based on weighting absorbance and similar sample
CN104089911B (en) Spectral model transmission method based on one-variable linear regression
Fortin et al. Chironomid-environment relations in northern North America
CN103364364B (en) Crude oil property rapid detection method based on recombination prediction technology
CN109345007B (en) Advantageous reservoir development area prediction method based on XGboost feature selection
CN101576485A (en) Analytical method of multi-source spectrum fusion water quality
CN103487411A (en) Method for recognizing steel grade by combining random forest algorithm with laser-induced breakdown spectroscopy
CN104990893B (en) A kind of gasoline octane value detection method based on similar differentiation
Palacios et al. Modeling the temperature‐nitrate relationship in the coastal upwelling domain of the California Current
CN104020135A (en) Calibration model establishing method based on near infrared spectrum
Wang et al. Review of the chemometrics application in oil-oil and oil-source rock correlations
CN105424641A (en) Crude oil type near infrared spectrum identification method
CN108573105A (en) The method for building up of soil heavy metal content detection model based on depth confidence network
CN106153871A (en) A kind of OIL SOURCE CORRELATION method
CN102954946B (en) By the method for infrared spectrum measurement sulfur content in crude oil
CN104502302A (en) Terahertz time-domain-waveform multiparameter-combined quantitative analysis method for mixed oil
CN107505282A (en) A kind of method for improving oil product near-infrared modeling robustness
Bender et al. Microstructure alignment of wood density profiles: an approach to equalize radial differences in growth rate
CN107389645B (en) The method that the Fisher model that wavelet transform parses oil product fluorescent characteristic identifies marine oil overflow
CN103049664A (en) Temperature interpolation method based on position classification
CN103927438A (en) Successive projection algorithm based near-infrared wavelength variable selecting method
CN113720952B (en) Method, device, equipment and medium for generating image plate for reservoir interpretation evaluation
Coy et al. Identifying patterns in multicomponent signals by extended cross correlation
Lehtinen et al. Nucleation rate and vapor concentration estimations using a least squares aerosol dynamics method

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