CN110108654A - 一种混合气体成分的非线性分析方法 - Google Patents

一种混合气体成分的非线性分析方法 Download PDF

Info

Publication number
CN110108654A
CN110108654A CN201910401120.5A CN201910401120A CN110108654A CN 110108654 A CN110108654 A CN 110108654A CN 201910401120 A CN201910401120 A CN 201910401120A CN 110108654 A CN110108654 A CN 110108654A
Authority
CN
China
Prior art keywords
gas
mixed gas
absorbance
concentration
wavelength
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
CN201910401120.5A
Other languages
English (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 Institute of Technology
Original Assignee
Nanjing Institute of Technology
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 Institute of Technology filed Critical Nanjing Institute of Technology
Priority to CN201910401120.5A priority Critical patent/CN110108654A/zh
Publication of CN110108654A publication Critical patent/CN110108654A/zh
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
    • 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/314Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths
    • 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
    • G01N21/314Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths
    • G01N2021/3148Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths using three or more wavelengths

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

本发明公开了一种混合气体成分的非线性分析方法,包括如下步骤:S1预估混合气体的组成成分,获取其中每种气体成分的光谱情况,建立二次拟合公式;S2利用MATLAB进行二次方程拟合,得到混合气体浓度与吸光度的非线性关系;S3测量混合气体在不同波长上的吸光度;S4将混合气体在不同波长上的吸光度带入到S2中得到的气体浓度与吸光度的非线性关系,计算得到混合气体中各气体的浓度。本发明的方法有较高的精度,并且对于不包含的气体本方法最终的计算结果自动为零,可以实现对混合气体的盲检测。

Description

一种混合气体成分的非线性分析方法
技术领域
本发明属于混合气体分析技术领域,具体涉及一种混合气体成分的非线性分析方法。
背景技术
近几年,气体检测技术快速发展,在矿井、大气环境监测、医疗、食品、航空航天方面都有广泛的应用。目前对于单一气体的浓度检测已经非常成熟,但是在实际环境中往往多种气体是共存的,故气体的种类和浓度很难从单一传感器的响应来判断。
当前,对多组分气体的定量分析常采用气相色谱(gas chromatogram,GC)和气相-质谱联用(gas chromatogram-mass gas chromatogram,GC-MS)方法、电化学探测法和吸收光谱法。电化学方法需要对待测气体进行采样,不适用于现场监测;GC和GC-MS法需要专门的进样装置,其过程复杂、时间较长,无法满足实时自动化连续监测的需要。光谱分析技术与化学计量学方法的结合极大地提高了定量分析的灵敏度、准确性和可靠性,适用于现场快速检测和实时在线分析,但是,多组分气体的带状吸收谱造成多成分谱的混叠,需用化学计量学方法建立多元回归模型进行分析和预测。常用的校正方法包括多元线性回归(multivariate linear regression,MLR),主成分回归(principle componentregression,PCR),偏最小二乘法(partial least square,PLS),人工神经网络。最小二乘算法、主成分分析算法和主成分回归算法等都是相关性的分析方法,不适用于非线性的、不相关的变量分析,因此这些算法在处理多组分气体时,很难达到预期的实验效果,分析结果精度不高。
发明内容
本发明的针对现有技术中的不足,提出利用已知的单成分气体光谱情况,建立浓度与吸光度的非线性关系,再利用朗伯定律中光谱的线性叠加原理,把混合气体的吸收光谱,看成是单种气体吸收光谱的线性组合,从而解出混合气体中各组份的浓度。
为实现上述目的,本发明采用以下技术方案:
一种混合气体成分的非线性分析方法,包括如下步骤:
S1预估混合气体的组成成分,获取其中每种气体成分的光谱情况,建立二次拟合公式;
S2利用MATLAB进行二次方程拟合,得到混合气体浓度与吸光度的非线性关系;
S3测量混合气体在不同波长上的吸光度;
S4将混合气体在不同波长上的吸光度带入到S2中得到的气体浓度与吸光度的非线性关系,计算得到混合气体中各气体的浓度。
为优化上述技术方案,采取的具体措施还包括:
步骤S1中每种气体成分的光谱情况包括每种气体的浓度、波长、吸光度。
步骤S1包括测量每种气体在不同浓度对应不同波长情况下的吸光度,并构建相应的二次拟合公式;若有n种预估气体,则测量n种波长情况下,每一种波长对应的每种气体几种不同浓度的吸光度。
步骤S2中利用MATLAB进行二次方程拟合,得到二次拟合公式的系数,进而得到混合气体浓度与吸光度的非线性关系。
在某一波长上,混合气体浓度与吸光度的非线性关系为
其中,A(λ)为混合气体在波长λ上的吸光度,Ai(λ)为混合气体中的第i种气体在波长λ上的吸光度,fi表示第i种气体浓度与吸光度之间的函数关系,ci代表混合气体中第i种气体的浓度。
最终得到的混合气体浓度与吸光度的非线性关系为
其中,fin,ci)为混合气体中第i种气体在波长λn时,浓度ci与吸光度A(λn)的拟合公式,只有一个未知数ci;A(λn)通过测量获得,而方程中的未知数为n种气体的浓度,共计n个方程,通过求解方程组获得每种气体的浓度。混合气体中不包含的预估气体成分浓度计算结果为0。
本发明的有益效果是:本发明对于混合气体只需知道可能含有的成分,利用已知的单成分气体光谱情况,建立浓度与吸光度的非线性关系,再利用朗伯定律中光谱的线性叠加原理,把混合气体的吸收光谱,看成是单种气体吸收光谱的线性组合,从而解出混合气体中各组份的浓度;对于不包含的气体本方法最终的计算结果自动为零,可以实现对混合气体的盲检测,该方法有较高的精度。
附图说明
图1是本发明的非线性分析方法的流程图。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
本发明提供的混合气体成分的非线性分析方法,包括如下步骤:
S1预估混合气体的组成成分,获取其中每种气体成分的光谱情况,建立二次拟合公式。
在实际测量中,往往无法知道混合气体是否含有某种气体。但是对于一定的气体(如汽车尾气、燃煤电厂烟气等),其可能的存在的气体成分是已知的。首先预估待分析的混合气体中可能含有的气体成分,获取其中每种气体成分的光谱情况,每种气体成分的光谱情况包括每种气体的浓度、波长、吸光度等参数;
在实际测量中吸光度-浓度曲线可能很难获取解析式,对精确估计造成难度。处理时可以采用高阶拟合的方法,拟合的阶数和运算量成正比,为了提高时效性同时保持精确度,采用二阶拟合。
测量每种气体在不同浓度对应不同波长情况下的吸光度,并构建相应的二次拟合公式;若预估气体有n种,则测量n种波长情况下,每一种波长对应的每种气体几种不同浓度的吸光度。
S2利用MATLAB进行二次方程拟合,得到二次拟合公式的系数,进而得到混合气体浓度与吸光度的非线性关系;
对于混合气体,由朗伯-比尔定律可知,混合气体在某一波长上的吸光度可以由其组成成分的吸光度线性叠加而来,同时,每个气体成分的吸光度与成分的浓度相关,可以表示为浓度的函数,如式1所示,在某一波长上,混合气体浓度与吸光度的非线性关系为
其中,A(λ)为混合气体在波长λ上的吸光度,Ai(λ)为混合气体中的第i种气体在波长λ上的吸光度,fi表示第i种气体浓度与吸光度之间的函数关系,ci代表混合气体中第i种气体的浓度。
由于不同气体组份,可能出现光谱的重叠,导致在测量混合气体时,无法仅利用单一波长测量某种气体的浓度。通常,采用利用已知单一气体浓度,测量多个相应的吸光度,构成回归方程,常用的方法包括多元线性分析、主成分分析、最小二乘法、偏最小二乘回归分析等。这些方法都是基础朗伯-比尔定律,认为气体吸光度与浓度之间是线性关系。但是在实际测量中,有很多的因素导致了浓度和吸光度之间的关系偏离了比尔朗伯定律出现了非线性关系,如气体浓度升高时,气体分子间的平均距离减小了,影响到了相邻质点间电荷分布情况,使分子的吸光能力发生改变,通常吸光能力减弱了。同时光谱仪响应的非线性也会影响到测量结果。
针对这种情况,朗伯-比尔定律中气体的吸光度与气体的浓度成正比关系已经不能准确描述该过程。
故对于多种成分的混合气体,选择每种气体成分的特征波长,并利用线性叠加的方法构成下列方程组:
最终得到的混合气体浓度与吸光度的非线性关系为
其中,fin,ci)为混合气体中第i种气体在波长λn时,浓度ci与吸光度A(λn)的拟合公式,只有一个未知数ci;A(λn)通过测量获得,而方程中的未知数为n个气体的浓度,共计n个方程,通过求解方程组就可以获得每种气体的浓度。
S3测量混合气体在不同波长上的吸光度;
S4将混合气体在不同波长上的吸光度带入到S2中得到的气体浓度与吸光度的非线性关系式(3)中,计算得到混合气体中各气体的浓度ci;对于混合气体中不包含的预估气体成分,浓度计算结果为0。
具体实施例如下:
采用SO2、NO2作为实验对象,这两种气体在190-290nm波段都有吸收,利用光谱吸收进行测量时,存在相互的干扰,需要进行混合气体的成分分析。
分别测量SO2、NO2在浓度为100ppm、300ppm、500ppm时,波长分别为273.33nm和231.33nm的吸光度,然后建立二次拟合公式:A(λ)=a0+a1c+a2c2,式中,A(λ)表示某一波长下的吸光度,a0、a1、a2为拟合系数,c为气体浓度。
表1气体吸光度
表1表示,对于SO2气体,在波长为273.33nm、浓度浓度为100ppm时,吸光度为0.017154,以此类推。
利用MATLAB进行二次方程拟合,其系数如下表
表2拟合系数
利用表2的拟合系数,建立气体浓度与吸光度的公式,代入公式(3)得:
在气室中充入8组不同浓度的SO2、NO2混合气体(具体浓度如表3所示),利用光栅分光测量在273.33nm和231.33nm波长上的吸光度,如表4所示。
编号 SO<sub>2</sub> NO<sub>2</sub> 编号 SO<sub>2</sub> NO<sub>2</sub>
1 90ppm 510ppm 5 330ppm 270ppm
2 150ppm 450ppm 6 390ppm 210ppm
3 210ppm 390ppm 7 450ppm 150ppm
4 270ppm 330ppm 8 510ppm 90ppm
表3混合气体各成分的浓度
表4混合气体的吸光度
将表4的结果代入公式4,利用MATLAB求解,去除不合理的情况(如结果为复数、负数的情况),结果如表5所示。
表5混合气体测量结果
通过表5可以看出,8组混合气体中,计算得到的SO2与NO2浓度分别与表3中充入的气体浓度基本相同,验证了本方法有效且具有较高的精确度,操作简单。
本发明的又一实施例如下:
在实际测量中,往往无法知道混合气体是否含有某种气体。但是对于一定的气体(如汽车尾气、燃煤电厂烟气等),其可能的存在的气体成分是已知的。在使用前可以将可能的气体都建立相应的拟合公式,,然后测量混合气体的吸光度,计算得到各种气体的浓度。
实验采用SO2、NO2、NO、NH3四种气体作为混合气体的可能成分,选择其中两种气体混合,作为测试气体,并利用光谱仪测量得到在273.33nm、231.33nm、225.88nm、208.23nm的吸光度,最后通过计算得到相应的气体浓度,结果如下表所示。
表6种类未知混合气体的测量结果(单位:ppm)
从表6中可以看到,当混合气体中不含有某种成分时,计算结果接近为零,而其他气体的误差小于3%。表明本文的方法可以用于特定气体源混合气体的盲检测,如汽车尾气、燃煤电厂烟气等,其组成的成分是已知的。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。

Claims (7)

1.一种混合气体成分的非线性分析方法,其特征在于,包括如下步骤:
S1预估混合气体的组成成分,获取其中每种气体成分的光谱情况,建立二次拟合公式;
S2利用MATLAB进行二次方程拟合,得到混合气体浓度与吸光度的非线性关系;
S3测量混合气体在不同波长上的吸光度;
S4将混合气体在不同波长上的吸光度带入到S2中得到的气体浓度与吸光度的非线性关系,计算得到混合气体中各气体的浓度。
2.根据权利要求1所述的混合气体成分的非线性分析方法,其特征在于,S1中每种气体成分的光谱情况包括每种气体的浓度、波长、吸光度。
3.根据权利要求2所述的混合气体成分的非线性分析方法,其特征在于,S1包括测量每种气体在不同浓度对应不同波长情况下的吸光度,并构建相应的二次拟合公式;若有n种预估气体,则测量n种波长情况下,每一种波长对应的每种气体几种不同浓度的吸光度。
4.根据权利要求2所述的混合气体成分的非线性分析方法,其特征在于,S2中利用MATLAB进行二次方程拟合,得到二次拟合公式的系数,进而得到混合气体浓度与吸光度的非线性关系。
5.根据权利要求2所述的混合气体成分的非线性分析方法,其特征在于,在某一波长上,混合气体浓度与吸光度的非线性关系为
其中,A(λ)为混合气体在波长λ上的吸光度,Ai(λ)为混合气体中的第i种气体在波长λ上的吸光度,fi表示第i种气体浓度与吸光度之间的函数关系,ci代表混合气体中第i种气体的浓度。
6.根据权利要求5所述的混合气体成分的非线性分析方法,其特征在于,最终得到的混合气体浓度与吸光度的非线性关系为
其中,fin,ci)为混合气体中第i种气体在波长λn时,浓度ci与吸光度A(λn)的拟合公式,只有一个未知数ci;A(λn)通过测量获得,而方程中的未知数为n个气体的浓度,共计n个方程,通过求解方程组获得每种气体的浓度。
7.根据权利要求1所述的混合气体成分的非线性分析方法,其特征在于,混合气体中不包含的预估气体成分浓度计算结果为0。
CN201910401120.5A 2019-05-14 2019-05-14 一种混合气体成分的非线性分析方法 Pending CN110108654A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910401120.5A CN110108654A (zh) 2019-05-14 2019-05-14 一种混合气体成分的非线性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910401120.5A CN110108654A (zh) 2019-05-14 2019-05-14 一种混合气体成分的非线性分析方法

Publications (1)

Publication Number Publication Date
CN110108654A true CN110108654A (zh) 2019-08-09

Family

ID=67490077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910401120.5A Pending CN110108654A (zh) 2019-05-14 2019-05-14 一种混合气体成分的非线性分析方法

Country Status (1)

Country Link
CN (1) CN110108654A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112378873A (zh) * 2020-10-29 2021-02-19 湖北锐意自控系统有限公司 紫外气体分析方法及紫外气体分析仪
CN117147475A (zh) * 2023-10-30 2023-12-01 杭州泽天春来科技有限公司 气体分析仪的多目标气体分析方法、系统及可读介质
CN117347301A (zh) * 2023-09-28 2024-01-05 湖南五凌电力科技有限公司 混合气体检测装置、方法、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105938090A (zh) * 2016-05-20 2016-09-14 璁告椽 一种多光谱检测混合液体的混合比例的方法及其设备
CN108680521A (zh) * 2018-05-18 2018-10-19 黎明职业大学 一种成分分析方法及颜色监测方法
CN109060702A (zh) * 2018-07-27 2018-12-21 中国石油集团长城钻探工程有限公司录井公司 红外光谱非线性定量分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105938090A (zh) * 2016-05-20 2016-09-14 璁告椽 一种多光谱检测混合液体的混合比例的方法及其设备
CN108680521A (zh) * 2018-05-18 2018-10-19 黎明职业大学 一种成分分析方法及颜色监测方法
CN109060702A (zh) * 2018-07-27 2018-12-21 中国石油集团长城钻探工程有限公司录井公司 红外光谱非线性定量分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
唐有祺: "《当代化学前沿》", 30 June 1997, 中国致公出版社 *
苏少林: "《仪器分析》", 31 August 2007, 中国环境科学出版社 *
陈智勇 等: "二次曲线拟合在比色分析中的应用", 《中国生物制品学杂志》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112378873A (zh) * 2020-10-29 2021-02-19 湖北锐意自控系统有限公司 紫外气体分析方法及紫外气体分析仪
CN112378873B (zh) * 2020-10-29 2021-11-16 湖北锐意自控系统有限公司 紫外气体分析方法及紫外气体分析仪
CN117347301A (zh) * 2023-09-28 2024-01-05 湖南五凌电力科技有限公司 混合气体检测装置、方法、电子设备及存储介质
CN117147475A (zh) * 2023-10-30 2023-12-01 杭州泽天春来科技有限公司 气体分析仪的多目标气体分析方法、系统及可读介质
CN117147475B (zh) * 2023-10-30 2024-01-30 杭州泽天春来科技有限公司 气体分析仪的多目标气体分析方法、系统及可读介质

Similar Documents

Publication Publication Date Title
CN104568836B (zh) 基于多种光谱技术融合的低浓度、多组分气体检测方法
CN110108654A (zh) 一种混合气体成分的非线性分析方法
EP1861691B1 (en) Method to reduce background noise in a spectrum
US7251037B2 (en) Method to reduce background noise in a spectrum
CN104819958B (zh) 傅里叶变换红外光谱气体检测中自动消除水汽干扰的方法与装置
CN102183468B (zh) 多组分气体分析的干扰修正与浓度反演方法
Stec et al. Quantification of fire gases by FTIR: experimental characterisation of calibration systems
CN101825567A (zh) 一种近红外光谱和拉曼光谱波长的筛选方法
CN103592257B (zh) 一种用近红外光谱快速测定木器涂料稀释剂中甲苯、乙苯和二甲苯含量的方法
CN102435556A (zh) 用于复杂非均相混合物体系的准确光谱定量分析方法
CN105203497B (zh) 一种由近红外光谱预测脱硫胺液中硫化氢含量的方法
CN105136682B (zh) 一种燃油中芳烃含量快速、准确测定的新方法
CN102495042A (zh) 一种粉末混合物的拉曼光谱准确定量分析方法
EP3492908B1 (en) Augmented raman analysis of a gas mixture
US3727048A (en) Chemical tracer method of and structure for determination of instantaneous and total fluid flow mass and volume
CN109115700B (zh) 炼锌废水金属离子检测的显色剂用量优化方法
Tang et al. On-line multi-component alkane mixture quantitative analysis using Fourier transform infrared spectrometer
CN101907563A (zh) 基于紫外发光二极管的二氧化硫分析仪及分析方法
CN105241823A (zh) 基于稀疏表示的火电厂烟气光谱定量分析方法
CN104777125B (zh) 火场气体产物及气体灭火剂定量分析方法及分析系统
Abbaspour et al. Simultaneous determination of antimony and bismuth by beta-correction spectrophotometry and an artificial neural network algorithm
CN115575337B (zh) 高精度大气co2浓度观测数据标较方法、系统及设备
CN107356556A (zh) 一种近红外光谱定量分析的双集成建模方法
CN103776789B (zh) 一种基于紫外可见光谱的未知污染物预警方法
JP2003014627A (ja) 多成分分析計の検量線作成方法

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190809