CN105572067A - Flue gas concentration measuring method based on spectrum analysis - Google Patents

Flue gas concentration measuring method based on spectrum analysis Download PDF

Info

Publication number
CN105572067A
CN105572067A CN201510931551.4A CN201510931551A CN105572067A CN 105572067 A CN105572067 A CN 105572067A CN 201510931551 A CN201510931551 A CN 201510931551A CN 105572067 A CN105572067 A CN 105572067A
Authority
CN
China
Prior art keywords
concentration value
optical thickness
ultraviolet band
concentration
actual
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
Application number
CN201510931551.4A
Other languages
Chinese (zh)
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.)
Chongqing Chuanyi Automation Co Ltd
Original Assignee
Chongqing Chuanyi Automation 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 Chongqing Chuanyi Automation Co Ltd filed Critical Chongqing Chuanyi Automation Co Ltd
Priority to CN201510931551.4A priority Critical patent/CN105572067A/en
Publication of CN105572067A publication Critical patent/CN105572067A/en
Pending legal-status Critical Current

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/33Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using ultraviolet 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
    • 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
    • G01N2021/3129Determining multicomponents by multiwavelength light

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention discloses a flue gas concentration measuring method based on spectrum analysis. According to the flue gas concentration measuring method disclosed by the invention, by carrying out band division treatment on flue gas, an instrument SO2 concentration value can be inverted in a first ultraviolet band, an actual standard SO2 concentration value is utilized for calibrating the instrument SO2 concentration value so as to obtain an accurate actual SO2 concentration value, the differential optical thickness of SO2 in a second ultraviolet band can be calculated by bringing in the actual SO2 concentration value, the differential optical thickness of NO can be calculated, and an accurate actual NO concentration value can be obtained by combining the differential optical thickness of the NO in the second ultraviolet band, so that mutual interference between the SO2 and the NO is eliminated, the measuring accuracy of SO2 and NO concentration in short-optical length and low-concentration flue gas is increased, the optical length is prevented from being increased or a high-accuracy spectrograph is prevented from being used, and the equipment cost is reduced.

Description

Based on the flue gas concentration measuring method of spectral analysis
Technical field
The present invention relates to flue gas concentration measuring field, be specifically related to a kind of flue gas concentration measuring method based on spectral analysis.
Background technology
Monitoring the harmful gas in flue gas is an importance of efforts at environmental protection, current flue gas monitoring technology is mainly divided into modern chemistry measuring technique and spectral measurement methods, difference absorption spectrum technology (DOAS, DifferentialOpticalAbsorptionSpectroscopy) as the representative of spectral measurement methods, because its measuring principle is simple, fast response time, the advantages such as the real time on-line monitoring of gas can be realized, it is widely used in the on-line monitoring of flue gas concentration, such as measuring the gas concentration of SO2 and NO in flue gas; Based on the gas concentration measuring method of DOAS, the gasmetry of long light path high concentration has achieved very large achievement, but for the gas concentration measurement of low concentration, short light path, because the signal to noise ratio (S/N ratio) of detection signal is low, the error of its measurement result is large, in the prior art, improve the measuring accuracy of light concentration gas, by increasing light path or using high-precision spectrometer, but can bring again the increase of cost like this.
Therefore, for overcoming the above problems, needing a kind of flue gas concentration measuring method based on spectral analysis, under ensureing short light path while measuring equipment cost can be reduced, to the gas of low concentration, there is higher measuring accuracy.
Summary of the invention
In view of this, the object of the invention is to overcome defect of the prior art, the flue gas concentration measuring method based on spectral analysis is provided, under ensureing short light path while measuring equipment cost can be reduced, to the gas of low concentration, there is higher measuring accuracy.
Flue gas concentration measuring method based on spectral analysis of the present invention, comprises the steps:
A. in the first ultraviolet band, obtain SO2 differential absorption cross-section, obtain instrument SO2 concentration value by DOAS algorithm and least square method, with actual standard SO2 concentration value calibrating instrument SO2 concentration value, set up the linear relationship between actual SO2 concentration value and instrument SO2 concentration value;
B. in the second ultraviolet band, obtain the differential optical thickness that SO2 and NO mixed gas is total, calculate the differential optical thickness of SO2 in the second ultraviolet band according to the actual SO2 concentration value in step a, deduct the differential optical thickness obtaining NO in the second ultraviolet band with total differential optical thickness.
C. in described second ultraviolet band, obtain the differential optical thickness of NO differential absorption cross-section NO in integrating step b, by DOAS algorithm and least square method computing equipment NO concentration value, utilize actual standard NO concentration value calibrating instrument NO concentration value, set up the linear relationship between actual NO concentration value and instrument NO concentration value;
Wherein, the light SO2 in described first ultraviolet band can absorb and NO can not absorb, and in described second ultraviolet band, light SO2 and NO all can absorb.
Further, described first ultraviolet band is positioned at wavelength at 285nm-310nm, and described second ultraviolet band is positioned at wavelength at 200-230nm.
Further, wherein in step a and step c, actual SO2 concentration value and actual NO concentration value are tried to achieve according to each self-corresponding instrument SO2 concentration value and instrument NO concentration value after all adopting demarcate at 2 respectively.
Further, the method wherein obtaining the differential absorption cross-section of SO2 and NO in step a and c is: transfer corresponding absorption cross section data from HITRAN database, and low order matching is carried out to absorption cross section data, obtain the low frequency part data of absorption cross section data, absorption cross section data are deducted low frequency part data and obtain corresponding differential absorption cross-section data.
Further, wherein in step b, obtaining total differential optical thickness approach is: in the first ultraviolet band, according to zero point spectrum and absorption spectrum obtain total optical thickness, low order matching is carried out to total optical thickness and obtains low order match value, utilize total optical thickness to deduct low order match value and obtain total differential optical thickness; In like manner, in the second ultraviolet band, obtain the differential optical thickness of SO2.
The invention has the beneficial effects as follows: a kind of flue gas concentration measuring method based on spectral analysis disclosed by the invention, by flue gas carried out subrane process, inverting instrument SO2 concentration value in the first ultraviolet band, and utilize actual standard SO2 concentration value calibrating instrument SO2 concentration value and obtain accurate actual SO2 concentration value, and then actual SO2 concentration value brought in the second ultraviolet band, calculate SO2 differential optical thickness in the second ultraviolet band, calculate NO differential optical thickness, NO differential optical thickness in conjunction with the second ultraviolet band obtains accurate actual NO concentration value, the mutual interference of SO2 and NO got rid of, improve short light path, the measuring accuracy of SO2 and NO concentration in low-concentration flue gas, avoid increasing light path or using high-precision spectrometer, reduce equipment cost.
Accompanying drawing explanation
Below in conjunction with drawings and Examples, the invention will be further described:
Fig. 1 is the calculation flow chart of step a in the present invention;
Fig. 2 is the calculation flow chart of step b in the present invention;
Fig. 3 is the calculation flow chart of step c in the present invention;
The measured value of the SO2 concentration that Fig. 4 is the normal concentration that uses the algorithm in the present invention to calculate is 49ppm;
The measured value of the NO concentration that Fig. 5 is the normal concentration that uses the algorithm in the present invention to calculate is 25ppm.
Embodiment
Fig. 1 is structural representation of the present invention, Fig. 1 is the calculation flow chart of step a in the present invention, Fig. 2 is the calculation flow chart of step b in the present invention, Fig. 3 is the calculation flow chart of step c in the present invention, the measured value of the SO2 concentration that Fig. 4 is the normal concentration that uses the algorithm in the present invention to calculate is 49ppm, the measured value of the NO concentration that Fig. 5 is the normal concentration that uses the algorithm in the present invention to calculate is 25ppm, the ordinate of Fig. 4 and Fig. 5 represents concentration value and unit is ppm, and horizontal ordinate represents pendulous frequency; As shown in the figure, the flue gas concentration measuring method based on spectral analysis in the present embodiment, comprises the steps:
A. in the first ultraviolet band, obtain SO2 differential absorption cross-section, instrument SO2 concentration value is obtained by DOAS algorithm and least square method, with actual standard SO2 concentration value calibrating instrument SO2 concentration value, set up the linear relationship between actual SO2 concentration value and instrument SO2 concentration value, wherein obtaining instrument SO2 concentration value by DOAS algorithm and least square method is prior art, does not repeat them here;
B. in the second ultraviolet band, obtain the differential optical thickness that SO2 and NO mixed gas is total, calculate the differential optical thickness of SO2 in the second ultraviolet band according to the actual SO2 concentration value in step a, deduct the differential optical thickness obtaining NO in the second ultraviolet band with total differential optical thickness; Optical thickness is composed (or zero point spectrum, spectrum when not having the spectrum of gas absorption or generally refer to logical nitrogen) by light and to be obtained divided by absorption spectrum (having spectrum when absorbing gas) that (wherein optical thickness equals ln (I 0/ I), I 0for initial beam intensity, the light intensity after I smoke absorption); Differential optical thickness by deducting SO2 in the second ultraviolet band eliminates SO2 to the measurements interference of NO, improve NO measuring accuracy, and be the inverse operation of step a algorithm by actual SO2 concentration value calculating differential optical thickness of SO2 in the second ultraviolet band, do not repeat them here.
C. in described second ultraviolet band, obtain the differential optical thickness of NO differential absorption cross-section NO in integrating step b, by DOAS algorithm and least square method computing equipment NO concentration value, utilize actual standard NO concentration value calibrating instrument NO concentration value, set up the linear relationship between actual NO concentration value and instrument NO concentration value;
Wherein, the light SO2 in described first ultraviolet band can absorb and NO can not absorb, and in described second ultraviolet band, light SO2 and NO all can absorb, by flue gas carried out subrane process, inverting instrument SO2 concentration value in the first ultraviolet band, and utilize actual standard SO2 concentration value calibrating instrument SO2 concentration value and obtain accurate actual SO2 concentration value, and then actual SO2 concentration value brought in the second ultraviolet band, calculate SO2 differential optical thickness in the second ultraviolet band, calculate NO differential optical thickness, NO differential optical thickness in conjunction with the second ultraviolet band obtains accurate actual NO concentration value, the mutual interference of SO2 and NO got rid of, improve short light path, the measuring accuracy of SO2 and NO concentration in low-concentration flue gas, avoid increasing light path or using high-precision spectrometer, reduce equipment cost.
In the present embodiment, described first ultraviolet band is positioned at wavelength at 285nm-310nm, and described second ultraviolet band is positioned at wavelength at 200-230nm; The present embodiment first ultraviolet band is preferably 285nm-310nm, and the second ultraviolet band is preferably 200-230nm, makes measurement result more accurate.
In the present embodiment, wherein in step a and step c, actual SO2 concentration value and actual NO concentration value are tried to achieve according to each self-corresponding instrument SO2 concentration value and instrument NO concentration value after all adopting demarcate at 2 respectively; Set up the straight-line equation of instrument concentration value and actual concentrations value by 2 demarcation, make computation process simple and convenient, and fast response time.
In the present embodiment, the method wherein obtaining the differential absorption cross-section of SO2 and NO in step a and c is: transfer corresponding absorption cross section data from HITRAN database, and low order matching is carried out to absorption cross section data, obtain the low frequency part data of absorption cross section data, absorption cross section data are deducted low frequency part data and obtain corresponding differential absorption cross-section data; Described HITRAN database is existing database, and data are complete and accuracy is high, and low order matching can adopt the process of existing low order fitting software, do not repeat them here.
In the present embodiment, wherein in step b, obtaining total differential optical thickness approach is: in the first ultraviolet band, according to zero point spectrum and absorption spectrum obtain total optical thickness, low order matching is carried out to total optical thickness and obtains low order match value, utilize total optical thickness to deduct low order match value and obtain total differential optical thickness; In like manner, in the second ultraviolet band, obtain the differential optical thickness of SO2; Low order matching can adopt the process of existing low order fitting software, does not repeat them here.
What finally illustrate is, above embodiment is only in order to illustrate technical scheme of the present invention and unrestricted, although with reference to preferred embodiment to invention has been detailed description, those of ordinary skill in the art is to be understood that, can modify to technical scheme of the present invention or equivalent replacement, and not departing from aim and the scope of technical solution of the present invention, it all should be encompassed in the middle of right of the present invention.

Claims (5)

1., based on a flue gas concentration measuring method for spectral analysis, it is characterized in that: comprise the steps:
A. in the first ultraviolet band, obtain SO2 differential absorption cross-section, obtain instrument SO2 concentration value by DOAS algorithm and least square method, with actual standard SO2 concentration value calibrating instrument SO2 concentration value, set up the linear relationship between actual SO2 concentration value and instrument SO2 concentration value;
B. in the second ultraviolet band, obtain the differential optical thickness that SO2 and NO mixed gas is total, calculate the differential optical thickness of SO2 in the second ultraviolet band according to the actual SO2 concentration value in step a, deduct the differential optical thickness obtaining NO in the second ultraviolet band with total differential optical thickness.
C. in described second ultraviolet band, obtain the differential optical thickness of NO differential absorption cross-section NO in integrating step b, by DOAS algorithm and least square method computing equipment NO concentration value, utilize actual standard NO concentration value calibrating instrument NO concentration value, set up the linear relationship between actual NO concentration value and instrument NO concentration value;
Wherein, the light SO2 in described first ultraviolet band can absorb and NO can not absorb, and in described second ultraviolet band, light SO2 and NO all can absorb.
2. the flue gas concentration measuring method based on spectral analysis according to claim 1, it is characterized in that: described first ultraviolet band is positioned at wavelength at 285nm-310nm, described second ultraviolet band is positioned at wavelength at 200-230nm.
3. the flue gas concentration measuring method based on spectral analysis according to claim 2, it is characterized in that: wherein in step a and step c, actual SO2 concentration value and actual NO concentration value are tried to achieve according to each self-corresponding instrument SO2 concentration value and instrument NO concentration value after all adopting demarcate at 2 respectively.
4. the flue gas concentration measuring method based on spectral analysis according to claim 1, it is characterized in that: the method wherein obtaining the differential absorption cross-section of SO2 and NO in step a and c is: transfer corresponding absorption cross section data from HITRAN database, and low order matching is carried out to absorption cross section data, obtain the low frequency part data of absorption cross section data, absorption cross section data are deducted low frequency part data and obtain corresponding differential absorption cross-section data.
5. the flue gas concentration measuring method based on spectral analysis according to claim 4, it is characterized in that: wherein in step b, obtaining total differential optical thickness approach is: in the first ultraviolet band, according to zero point spectrum and absorption spectrum obtain total optical thickness, low order matching is carried out to total optical thickness and obtains low order match value, utilize total optical thickness to deduct low order match value and obtain total differential optical thickness; In like manner, in the second ultraviolet band, obtain the differential optical thickness of SO2.
CN201510931551.4A 2015-12-14 2015-12-14 Flue gas concentration measuring method based on spectrum analysis Pending CN105572067A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510931551.4A CN105572067A (en) 2015-12-14 2015-12-14 Flue gas concentration measuring method based on spectrum analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510931551.4A CN105572067A (en) 2015-12-14 2015-12-14 Flue gas concentration measuring method based on spectrum analysis

Publications (1)

Publication Number Publication Date
CN105572067A true CN105572067A (en) 2016-05-11

Family

ID=55882446

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510931551.4A Pending CN105572067A (en) 2015-12-14 2015-12-14 Flue gas concentration measuring method based on spectrum analysis

Country Status (1)

Country Link
CN (1) CN105572067A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107271392A (en) * 2017-06-12 2017-10-20 重庆大学 One kind is based on gas differential absorption cross-section pattern measurement low concentration SO2Method
CN107796777A (en) * 2017-10-20 2018-03-13 西安思坦科技有限公司 A kind of data processing method of low concentration ultraviolet difference gas analyzer
CN112378873A (en) * 2020-10-29 2021-02-19 湖北锐意自控系统有限公司 Ultraviolet gas analysis method and ultraviolet gas analyzer
CN114018853A (en) * 2021-11-24 2022-02-08 青岛崂应海纳光电环保集团有限公司 Photometer gas chamber and gas analysis module

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101290289A (en) * 2008-05-30 2008-10-22 天津市蓝宇科工贸有限公司 Ultraviolet differential flue gas concentration measuring systems calibration method and enforcement device
CN201307092Y (en) * 2008-12-10 2009-09-09 杨秀芹 A ultraviolet absorption method portable stack gas, auto car tail gas measuring apparatus
JP4750580B2 (en) * 2006-02-27 2011-08-17 荏原実業株式会社 Method and apparatus for measuring ozone concentration
CN102788764A (en) * 2012-08-21 2012-11-21 南京埃森环境技术有限公司 Ultraviolet analyzer and detection method for low concentration smoke
CN104568836A (en) * 2015-01-26 2015-04-29 南京国电环保科技有限公司 Low-concentration and multi-component gas detection method based on integration of multiple spectrum technologies

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4750580B2 (en) * 2006-02-27 2011-08-17 荏原実業株式会社 Method and apparatus for measuring ozone concentration
CN101290289A (en) * 2008-05-30 2008-10-22 天津市蓝宇科工贸有限公司 Ultraviolet differential flue gas concentration measuring systems calibration method and enforcement device
CN201307092Y (en) * 2008-12-10 2009-09-09 杨秀芹 A ultraviolet absorption method portable stack gas, auto car tail gas measuring apparatus
CN102788764A (en) * 2012-08-21 2012-11-21 南京埃森环境技术有限公司 Ultraviolet analyzer and detection method for low concentration smoke
CN104568836A (en) * 2015-01-26 2015-04-29 南京国电环保科技有限公司 Low-concentration and multi-component gas detection method based on integration of multiple spectrum technologies

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107271392A (en) * 2017-06-12 2017-10-20 重庆大学 One kind is based on gas differential absorption cross-section pattern measurement low concentration SO2Method
CN107271392B (en) * 2017-06-12 2019-04-26 重庆大学 One kind being based on gas differential absorption cross-section pattern measurement low concentration SO2Method
CN107796777A (en) * 2017-10-20 2018-03-13 西安思坦科技有限公司 A kind of data processing method of low concentration ultraviolet difference gas analyzer
CN112378873A (en) * 2020-10-29 2021-02-19 湖北锐意自控系统有限公司 Ultraviolet gas analysis method and ultraviolet gas analyzer
CN112378873B (en) * 2020-10-29 2021-11-16 湖北锐意自控系统有限公司 Ultraviolet gas analysis method and ultraviolet gas analyzer
CN114018853A (en) * 2021-11-24 2022-02-08 青岛崂应海纳光电环保集团有限公司 Photometer gas chamber and gas analysis module
CN114018853B (en) * 2021-11-24 2024-03-26 青岛崂应海纳光电环保集团有限公司 Photometer air chamber and gas analysis module

Similar Documents

Publication Publication Date Title
CN104568836B (en) Low-concentration and multi-component gas detection method based on integration of multiple spectrum technologies
CN102004097B (en) Coal quality on-line detecting method based on dominating factor and combined with partial least squares method
CN102539377B (en) Intermediate infrared absorption spectra based method for multi-component mixed gas qualitative and quantitative analysis
CN102680020B (en) Gas parameter online measurement method based on wavelength modulation spectroscopy
CN105572067A (en) Flue gas concentration measuring method based on spectrum analysis
CN108181266B (en) TD L AS gas concentration detection method
JP6981817B2 (en) Spectroscopic analyzer and spectroscopic analysis method
US10274422B2 (en) Gas analysis apparatus and gas analysis method
CN102183468A (en) Interference correction and concentration inversion method of multi-component gas analysis
CN104198416A (en) Real-time compensation method of measurement errors caused by spectrograph wavelength drift
CN103940767B (en) Based on the gas concentration inversion method of multiple manifold study
CN104330378A (en) Method for correcting wave number drift of Fourier transform infrared spectrometer
CN105628649A (en) Method for extracting absorption characteristic peak in gas field monitoring
CN105548057A (en) Flue gas analysis and measurement method implemented through ultraviolet spectrum
CN103499391A (en) Spectrum measuring system
CN103616338A (en) Reconstruction method for atmosphere trace gas spatial distribution by differential optical absorption spectroscopy tomoscan
CN103558182B (en) A kind of method for laser gas in-line analyzer determination gas concentration
CN107478593A (en) A kind of low concentration of NO and SO2The concentration detection method of mixed gas
CN102495014A (en) Method for correcting spectral shift in differential optical absorption spectroscopy (DOAS) measurement
CN104048922A (en) Method for measuring polarization degree and polarization angle of fluorescence spectrum
CN104596962B (en) Cross calibrating method and system for different gases based on laser gas sensor
CN112666104A (en) DOAS-based gas concentration inversion method
CN104458630A (en) Data processing method and system for ultraviolet differential gas analyzer
CN104267018A (en) Method for processing gas concentration signal in Raman gas analyzer
US7751051B2 (en) Method for cross interference correction for correlation spectroscopy

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: 20160511

RJ01 Rejection of invention patent application after publication