CN117347297B - 大气nh3分析系统、方法及可读存储介质 - Google Patents
大气nh3分析系统、方法及可读存储介质 Download PDFInfo
- Publication number
- CN117347297B CN117347297B CN202311653398.4A CN202311653398A CN117347297B CN 117347297 B CN117347297 B CN 117347297B CN 202311653398 A CN202311653398 A CN 202311653398A CN 117347297 B CN117347297 B CN 117347297B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- absorbance
- atmospheric
- wavelength
- standard
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 80
- 238000002835 absorbance Methods 0.000 claims abstract description 40
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000003595 spectral effect Effects 0.000 claims abstract description 14
- 230000002452 interceptive effect Effects 0.000 claims abstract description 12
- 239000013618 particulate matter Substances 0.000 claims abstract description 8
- 238000012937 correction Methods 0.000 claims description 20
- 238000010521 absorption reaction Methods 0.000 claims description 17
- 230000008033 biological extinction Effects 0.000 claims description 5
- 238000012886 linear function Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 4
- 230000031700 light absorption Effects 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 20
- 238000005259 measurement Methods 0.000 description 11
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 5
- 238000001179 sorption measurement Methods 0.000 description 5
- MWUXSHHQAYIFBG-UHFFFAOYSA-N Nitric oxide Chemical compound O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 4
- 239000013307 optical fiber Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- PAWQVTBBRAZDMG-UHFFFAOYSA-N 2-(3-bromo-2-fluorophenyl)acetic acid Chemical compound OC(=O)CC1=CC=CC(Br)=C1F PAWQVTBBRAZDMG-UHFFFAOYSA-N 0.000 description 2
- 229910021529 ammonia Inorganic materials 0.000 description 2
- BFNBIHQBYMNNAN-UHFFFAOYSA-N ammonium sulfate Chemical compound N.N.OS(O)(=O)=O BFNBIHQBYMNNAN-UHFFFAOYSA-N 0.000 description 2
- 229910052921 ammonium sulfate Inorganic materials 0.000 description 2
- 235000011130 ammonium sulphate Nutrition 0.000 description 2
- 230000003197 catalytic effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000002253 acid Substances 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001658 differential optical absorption spectrophotometry Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
- Y02A50/20—Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种大气NH3分析系统、方法及可读存储介质,其中系统包括采集模块,用于采集NH3标准吸光度光谱、N2光谱和大气NH3吸收光谱;计算模块,用于基于不同波长对应的N2光谱强度和大气NH3吸收光谱强度计算得到不同波长对应的吸光度;构建模块,用于根据不同波长对应的吸光度、NH3标准吸光度、颗粒物散射强度和干扰气体吸光度构建NH3浓度分析模型;所述计算模块还用于通过最小二乘法获取NH3浓度分析模型的最优解,得到NH3标准吸光度对应的比例系数,并根据比例系数与标准NH3浓度的乘积计算得到大气NH3浓度。本发明可有效去除大气颗粒物以及干扰气体对大气NH3分析的影响,有效提升大气NH3分析的精度。
Description
技术领域
本发明属于分析检测技术领域,具体涉及一种大气NH3分析系统、方法及可读存储介质。
背景技术
氨气污染,是PM2.5指数被持续推高的重要原因。氨气(NH3),是一种无色透明气体,有刺激性臭味,极易溶于水,氨气与酸反应生成的硫酸铵、硝酸铵,这两种铵盐有利于大气中雾霾颗粒的形成。研究表明,重污染天气中硫酸铵、硝酸铵的质量总和约占PM2.5的40-60%,越严重的污染天气,比例越高。因此,NH3引起了国内外研究者的广泛关注,并对其进行了大量的测量分析研究。
准确测量大气中NH3的浓度需要克服以下三点:(1)NH3有较强的吸附性,容易附着在测量分析仪的壁面上,有时会与吸附面材料发生反应;(2)不同地点、时间大气中NH3的浓度差异变化很大;(3)NH3存在形态多样化,气态、液态和颗粒物中均会存在,有时在常温常压条件下就可相互转变,造成测量分析的复杂程度大。目前NH3的测量分析方法主要分为两类:一类是间接法,即通过催化氧化的方式,将样品气中的NH3转换成一氧化氮,通过测量一氧化氮,从而得到NH3在样品气中的浓度;另一类是直接法,当前直接可测量NH3的方法比较多,主要为光谱法;而光谱法的分析中又分为无采样和有采样两类。间接法受催化效率的影响,同时采样过程中容易吸附。光谱测量方法相比间接法,测量大气NH3受影响要小,尤其是无采样不受催化效率及吸附影响。开放DOAS技术是利用空气中气体分子的窄带吸收特性来进行气体成分识别,根据分子的窄带吸收特征光谱,进行反演计算,从而测定所需测量的NH3的浓度。选用开放的光路设计,可以大大提高测量的光程,这样可以得到更加低的检测限。由于开放的光路设计,大气颗粒物散射或吸收都可能测量光强发生变化,而这种变化很多时候被认为是气体吸收,造成气体浓度反演计算出现偏差。
发明内容
基于现有技术中存在的上述缺点和不足,本发明的目的之一是至少解决现有技术中存在的上述问题之一或多个,换言之,本发明的目的之一是提供满足前述需求之一或多个的一种大气NH3分析系统、方法及可读存储介质。
为了达到上述发明目的,本发明采用以下技术方案:
一种大气NH3分析系统,包括:
采集模块,用于采集NH3标准吸光度光谱、N2光谱和大气NH3吸收光谱;
计算模块,用于基于不同波长对应的N2光谱强度和大气NH3吸收光谱强度计算得到不同波长对应的吸光度;
构建模块,用于根据不同波长对应的吸光度、NH3标准吸光度、颗粒物散射强度和干扰气体吸光度构建NH3浓度分析模型;
所述计算模块还用于通过最小二乘法获取NH3浓度分析模型的最优解,得到NH3标准吸光度对应的比例系数,并根据比例系数与标准NH3浓度的乘积计算得到大气NH3浓度。
优选地,不同波长λ对应的吸光度为:
;
其中,为波长λ对应的N2光谱强度,/>为波长λ对应的大气NH3吸收光谱强度。
优选地,所述NH3浓度分析模型为:
;
其中,k为比例系数,为波长λ对应的NH3标准吸光度,为波长λ对应的颗粒物散射强度,β为散射强度权重,μ 550为550nm的消光系数,L为光程,k 1、k 2和d为干扰气体的二阶线性函数的拟合系数。
进一步地,所述的大气NH3分析系统,还包括:
选取模块,用于基于NH3标准吸光度光谱和大气NH3吸光度光谱的最大吸收峰,分别对应选取包含最大吸收峰的连续N个峰和谷对应的波长和/>,i∈[1,N],N为大于2的整数;其中,大气NH3吸光度光谱的表达式为:
;
所述计算模块用于对波长和/>进行重合度计算,得到波长重合度;
判断模块,用于判断波长重合度是否大于预设重合度;若是,则满足要求;若否,则光谱出现漂移。
进一步地,大气NH3分析系统,还包括:
校正模块,用于在光谱出现漂移时,对NH3标准吸光度光谱进行漂移校正;
更新模块,用于对利用校正后的NH3标准吸光度光谱对NH3浓度分析模型进行更新。
进一步地,对NH3标准吸光度光谱进行漂移校正,包括:
计算波长横向偏移量;
若波长横向偏移量大于零,则NH3标准吸光度光谱的横坐标波长向右移动,得到校正后的NH3标准吸光度光谱;
若波长横向偏移量小于零,则NH3标准吸光度光谱的横坐标波长向左移动,得到校正后的NH3标准吸光度光谱。
优选地,所述预设重合度不小于0.99。
本发明还提供一种大气NH3分析方法,应用于如上方案所述的大气NH3分析系统,所述大气NH3分析方法包括以下步骤:
S1、采集NH3标准吸光度光谱、N2光谱和大气NH3吸收光谱,并基于不同波长对应的N2光谱强度和大气NH3吸收光谱强度计算得到不同波长对应的吸光度;
S2、根据不同波长对应的吸光度、NH3标准吸光度、颗粒物散射强度和干扰气体吸光度构建NH3浓度分析模型;
S3、通过最小二乘法获取NH3浓度分析模型的最优解,得到NH3标准吸光度对应的比例系数;根据比例系数与标准NH3浓度的乘积计算得到大气NH3浓度。
进一步地,大气NH3分析方法,还包括以下步骤:
S41、基于NH3标准吸光度光谱和大气NH3吸光度光谱的最大吸收峰,分别对应选取包含最大吸收峰的连续N个峰和谷对应的波长和/>,i∈[1,N],N为大于2的整数;其中,大气NH3吸光度光谱的表达式为:
;
S42、对波长和/>进行重合度计算,得到波长重合度;
S43、判断波长重合度是否大于预设重合度;若是,则满足要求;若否,则光谱出现漂移,对NH3标准吸光度光谱进行漂移校正,利用校正后的NH3标准吸光度光谱对NH3浓度分析模型进行更新,并转至步骤S3。
本发明还提供一种可读存储介质,可读存储介质中存储有指令,当指令在计算机上运行时,使得计算机执行如上任一方案所述的大气NH3分析方法。
本发明与现有技术相比,有益效果是:
(1)本发明可以有效去除大气颗粒物以及干扰气体对大气NH3分析的影响,有效提升大气NH3分析的精度;
(2)本发明在光谱出现漂移时进行漂移校正,保证大气NH3分析的正常运行;
(3)本发明在漂移校正时无需通入标准NH3以重新采集NH3标准吸光度光谱,直接根据波长横向偏移量校正之后的NH3标准吸光度光谱作为新的NH3标准吸光度光谱即可。
附图说明
图1是现有的NH3分析仪的结构示意图;
图2是本发明实施例1的大气NH3分析系统的模块组成图;
图3是本发明实施例1的大气NH3分析方法的流程图
图4是本发明实施例1的光谱漂移的示例图;
图5是本发明实施例1的光谱漂移未校正的测量结果对比图;
图6是本发明实施例1的光谱漂移经过校正后的测量结果对比图。
具体实施方式
为了更清楚地说明本发明实施例,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
实施例1:
本实施例的大气NH3分析系统基于现有的NH3分析仪,如图1所示,NH3分析仪为开放测量光路设计,包括收发端和反射端,收发端包括光源、Y型光纤、遮光片、抛物面镜、光纤固定架、气体室、光谱仪和工控机,反射端包括平面反射镜,上述各器件的布设具体可参考现有技术,在此不赘述。
如图2所示,本实施例的大气NH3分析系统包括采集模块、计算模块、构建模块、选取模块、判断模块、校正模块和更新模块。
本实施例的采集模块用于采集NH3标准吸光度光谱、N2光谱和大气NH3吸收光谱。其中,NH3标准吸光度光谱通过通入标准NH3浓度的NH3进行测量并转换得到,在仪器安装之后标定时测量一次即可。另外,NH3分析仪具有零气管路,空气经过泵抽,先后经过颗粒物过滤器及NH3吸附管,进气气体室,通过电磁阀控制光路,即可获得N2光谱。本实施例的大气NH3吸收光谱根据待测大气环境下开放测量光路即可采集得到。
本实施例的计算模块用于基于不同波长对应的N2光谱强度和大气NH3吸收光谱强度计算得到不同波长对应的吸光度。
具体地,不同波长λ对应的吸光度为:
;
其中,为波长λ对应的N2光谱强度,/>为波长λ对应的大气NH3吸收光谱强度。
本实施例的构建模块用于根据不同波长对应的吸光度、NH3标准吸光度、颗粒物散射强度和干扰气体吸光度构建NH3浓度分析模型。
由于开放测量光路中测量光谱能量大小容易受到颗粒物及其它气体等因素的干扰,NH3浓度分析模型需要去除上述因素影响,故本实施例的大气NH3浓度分析模型具体设计为:
;
其中,k为比例系数,为波长λ对应的NH3标准吸光度,即NH3标准吸光度光谱中横坐标波长λ对应的纵坐标吸光度;/>表示波长λ对应的颗粒物散射强度,β是散射强度权重,/>表示波长λ对应的干扰气体的吸光度。
上述颗粒物散射强度根据波盖耳定律可表示为:
;
其中,为波长λ对应的消光系数,L为开放测量光路的光程。
具体地,消光系数与波长的关系可表示为:
;
其中,μ 550为550nm的消光系数。
上述采用二阶线性函数表示为:
;
其中,k 1、k 2和d为干扰气体的二阶线性函数的拟合系数。
综上,本实施例的NH3浓度分析模型为:
。
上述计算模块还用于通过最小二乘法获取NH3浓度分析模型的最优解,得到NH3标准吸光度对应的比例系数k,并根据比例系数k与标准NH3浓度的乘积计算得到大气NH3浓度C为:/>。
本实施例的选取模块用于基于NH3标准吸光度光谱和大气NH3吸光度光谱的最大吸收峰,分别对应选取包含最大吸收峰的连续N个峰和谷对应的波长和/>,i∈[1,N],N为大于2的整数。
作为示例,本实施例选取包含最大吸收峰的连续10个峰和谷对应的波长和/>,i∈[1,10],即NH3标准吸光度光谱对应的10个波长分别为/>,大气NH3吸光度光谱对应的10个波长分别为/>;
由于通过最小二乘法只能确定NH3浓度分析模型的最优解,并不能确定是唯一解,最终NH3浓度分析模型是约等于的关系。因此,大气NH3吸光度光谱的表达式为:
。
上述计算模块还用于对波长和/>进行重合度计算,得到波长重合度;具体地,重合度计算公式为:
;
其中,为/>与/>的协方差,/>为/>的方差,/>为的方差。
本实施例的判断模块用于判断波长重合度是否大于预设重合度,作为示例,预设重合度取值为0.99;若是,则满足要求;若否,则光谱出现漂移,之后通过校正模块进行校正。
本实施例的校正模块用于在光谱出现漂移时,对NH3标准吸光度光谱进行漂移校正。具体地,对NH3标准吸光度光谱进行漂移校正,包括:
计算波长横向偏移量;
若波长横向偏移量大于零,则NH3标准吸光度光谱的横坐标波长向右移动,得到校正后的NH3标准吸光度光谱;
若波长横向偏移量小于零,则NH3标准吸光度光谱的横坐标波长向左移动,得到校正后的NH3标准吸光度光谱。
本实施例的更新模块用于对利用校正后的NH3标准吸光度光谱对NH3浓度分析模型进行更新。
基于上述大气NH3分析系统,如图3所示,本实施例对应的大气NH3分析方法包括以下步骤:
(1)采集NH3标准吸光度光谱、N2光谱和大气NH3吸收光谱;
(2)基于不同波长对应的N2光谱强度和大气NH3吸收光谱强度计算得到不同波长对应的吸光度;
(3)根据不同波长对应的吸光度、NH3标准吸光度、颗粒物散射强度和干扰气体吸光度构建NH3浓度分析模型;
(4)通过最小二乘法获取NH3浓度分析模型的最优解,得到NH3标准吸光度对应的比例系数;根据比例系数与标准NH3浓度的乘积计算得到大气NH3浓度;
(5)基于NH3标准吸光度光谱和大气NH3吸光度光谱的最大吸收峰,作为示例,分别对应选取包含最大吸收峰的连续10个峰和谷对应的波长和/>,i∈[1,10],即NH3标准吸光度光谱对应的10个波长分别为/>,大气NH3吸光度光谱对应的10个波长分别为/>;
(6)对波长和/>进行重合度计算,得到波长重合度;
(7)判断波长重合度是否大于预设重合度0.99;若是,则满足要求,无漂移;若否,则光谱出现漂移,需要对NH3标准吸光度光谱进行漂移校正,并利用校正后的NH3标准吸光度光谱对NH3浓度分析模型进行更新,并转至上述步骤(4)。
其中,本实施例的大气NH3分析方法各步骤的详细过程具体可参考上述大气NH3分析系统中的详细描述,在此不赘述。
以下通过试验数据来验证上述大气NH3分析系统及方法的效果:
如图4所示,标准NH3浓度为1.01ppm的NH3标准吸光度光谱的偏移示例,即NH3标准吸光度光谱的横向出现向右漂移;
如图5所示,本实施例(即本项目)在未进行漂移校正时,测量结果与OPSIS AR500分析仪的测量结果存在较大的偏差。经过漂移校正之后,如图6所示,测量结果与OPSISAR500分析仪的测量结果吻合,有效保证了大气NH3分析精度。
本实施例的可读存储介质,可读存储介质中存储有指令,当指令在计算机上运行时,使得计算机执行上述大气NH3分析方法,实现智能分析。
实施例2:
本实施例的本实施例的大气NH3分析系统与实施例1的不同之处在于:
省略选取模块、判断模块、校正模块和更新模块,即简化系统构架,满足不同应用的需求;
其他模块可以参考实施例1。
相应地,本实施例的大气NH3分析方法,省略步骤(5)及其之后的步骤,满足不同应用的需求;
其他步骤可以参考实施例1。
本实施例的可读存储介质,可读存储介质中存储有指令,当指令在计算机上运行时,使得计算机执行本实施例的大气NH3分析方法,实现智能分析。
实施例3:
本实施例的大气NH3分析系统与实施例1的不同之处在于:
选取模块在选取连续的峰和谷的数量可根据实际应用精度需求确定,例如3个、5个、6个、8个等,满足不同应用的需求;
其他模块可以参考实施例1。
相应地,本实施例的大气NH3分析方法,选取连续的峰和谷的数量进行适应性变化,满足不同应用的需求;
其他步骤可以参考实施例1。
本实施例的可读存储介质,可读存储介质中存储有指令,当指令在计算机上运行时,使得计算机执行本实施例的大气NH3分析方法,实现智能分析。
以上所述仅是对本发明的优选实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。
Claims (8)
1.一种大气NH3分析系统,其特征在于,包括:
采集模块,用于采集NH3标准吸光度光谱、N2光谱和大气NH3吸收光谱;
计算模块,用于基于不同波长对应的N2光谱强度和大气NH3吸收光谱强度计算得到不同波长对应的吸光度;不同波长λ对应的吸光度为:
其中,为波长λ对应的N2光谱强度,Sλ为波长λ对应的大气NH3吸收光谱强度;
构建模块,用于根据不同波长对应的吸光度、NH3标准吸光度、颗粒物散射强度和干扰气体吸光度构建NH3浓度分析模型;所述NH3浓度分析模型为:
其中,k为比例系数,为波长λ对应的NH3标准吸光度,exp(-(550/λ)1.3*μ550*L)为波长λ对应的颗粒物散射强度,β为散射强度权重,μ550为550nm的消光系数,L为光程,k1、k2和d为干扰气体的二阶线性函数的拟合系数;
所述计算模块还用于通过最小二乘法获取NH3浓度分析模型的最优解,得到NH3标准吸光度对应的比例系数,并根据比例系数与标准NH3浓度的乘积计算得到大气NH3浓度。
2.根据权利要求1所述的大气NH3分析系统,其特征在于,还包括:
选取模块,用于基于NH3标准吸光度光谱和大气NH3吸光度光谱的最大吸收峰,分别对应选取包含最大吸收峰的连续N个峰和谷对应的波长和λi,i∈[1,N],N为大于2的整数;其中,大气NH3吸光度光谱Aλ的表达式为:Aλ=αλ-[β*exp(-(550/λ)1.3*μ550*L)+k1*λ2+k2*λ+d];
所述计算模块用于对波长和λi进行重合度计算,得到波长重合度;
判断模块,用于判断波长重合度是否大于预设重合度;若是,则满足要求;若否,则光谱出现漂移。
3.根据权利要求2所述的大气NH3分析系统,其特征在于,还包括:
校正模块,用于在光谱出现漂移时,对NH3标准吸光度光谱进行漂移校正;
更新模块,用于对利用校正后的NH3标准吸光度光谱对NH3浓度分析模型进行更新。
4.根据权利要求3所述的大气NH3分析系统,其特征在于,对NH3标准吸光度光谱进行漂移校正,包括:
计算波长横向偏移量
若波长横向偏移量大于零,则NH3标准吸光度光谱的横坐标波长向右移动△λ,得到校正后的NH3标准吸光度光谱;
若波长横向偏移量小于零,则NH3标准吸光度光谱的横坐标波长向左移动△λ,得到校正后的NH3标准吸光度光谱。
5.根据权利要求2所述的大气NH3分析系统,其特征在于,所述预设重合度不小于0.99。
6.一种大气NH3分析方法,应用于如权利要求1所述的大气NH3分析系统,其特征在于,所述大气NH3分析方法包括以下步骤:
S1、采集NH3标准吸光度光谱、N2光谱和大气NH3吸收光谱,并基于不同波长对应的N2光谱强度和大气NH3吸收光谱强度计算得到不同波长对应的吸光度;
S2、根据不同波长对应的吸光度、NH3标准吸光度、颗粒物散射强度和干扰气体吸光度构建NH3浓度分析模型;
S3、通过最小二乘法获取NH3浓度分析模型的最优解,得到NH3标准吸光度对应的比例系数;根据比例系数与标准NH3浓度的乘积计算得到大气NH3浓度。
7.根据权利要求6所述的大气NH3分析方法,其特征在于,还包括以下步骤:
S41、基于NH3标准吸光度光谱和大气NH3吸光度光谱的最大吸收峰,分别对应选取包含最大吸收峰的连续N个峰和谷对应的波长和λi,i∈[1,N],N为大于2的整数;其中,大气NH3吸光度光谱Aλ的表达式为:
Aλ=αλ-[β*exp(-(550/λ)1.3*μ550*L)+k1*λ2+k2*λ+d];
S42、对波长和λi进行重合度计算,得到波长重合度;
S43、判断波长重合度是否大于预设重合度;若是,则满足要求;若否,则光谱出现漂移,对NH3标准吸光度光谱进行漂移校正,利用校正后的NH3标准吸光度光谱对NH3浓度分析模型进行更新,并转至步骤S3。
8.一种可读存储介质,可读存储介质中存储有指令,其特征在于,当指令在计算机上运行时,使得计算机执行如权利要求6或7所述的大气NH3分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311653398.4A CN117347297B (zh) | 2023-12-05 | 2023-12-05 | 大气nh3分析系统、方法及可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311653398.4A CN117347297B (zh) | 2023-12-05 | 2023-12-05 | 大气nh3分析系统、方法及可读存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117347297A CN117347297A (zh) | 2024-01-05 |
CN117347297B true CN117347297B (zh) | 2024-04-02 |
Family
ID=89365333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311653398.4A Active CN117347297B (zh) | 2023-12-05 | 2023-12-05 | 大气nh3分析系统、方法及可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117347297B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003014625A (ja) * | 2001-06-28 | 2003-01-15 | Ishikawajima Harima Heavy Ind Co Ltd | So3,nh3同時連続濃度計 |
CN101290289A (zh) * | 2008-05-30 | 2008-10-22 | 天津市蓝宇科工贸有限公司 | 紫外差分烟气浓度测量系统校准方法及实施装置 |
CN101363796A (zh) * | 2007-12-29 | 2009-02-11 | 北京航空航天大学 | 一种二氧化碳气体分析仪及其分析方法 |
CN206074435U (zh) * | 2016-09-26 | 2017-04-05 | 三峡大学 | 一种基于水质cod兼容氨氮浓度实时监测装置 |
CN108489907A (zh) * | 2018-04-10 | 2018-09-04 | 中国科学院合肥物质科学研究院 | 一种氨气检测装置及方法 |
CN114235742A (zh) * | 2021-12-17 | 2022-03-25 | 中国石油大学(华东) | 基于呼吸气体大类标志物复合光谱检测系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10041926B2 (en) * | 2016-07-14 | 2018-08-07 | Chevron U.S.A. Inc. | Method for predicting total petroleum hydrocarbon concentration in soils |
-
2023
- 2023-12-05 CN CN202311653398.4A patent/CN117347297B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003014625A (ja) * | 2001-06-28 | 2003-01-15 | Ishikawajima Harima Heavy Ind Co Ltd | So3,nh3同時連続濃度計 |
CN101363796A (zh) * | 2007-12-29 | 2009-02-11 | 北京航空航天大学 | 一种二氧化碳气体分析仪及其分析方法 |
CN101290289A (zh) * | 2008-05-30 | 2008-10-22 | 天津市蓝宇科工贸有限公司 | 紫外差分烟气浓度测量系统校准方法及实施装置 |
CN206074435U (zh) * | 2016-09-26 | 2017-04-05 | 三峡大学 | 一种基于水质cod兼容氨氮浓度实时监测装置 |
CN108489907A (zh) * | 2018-04-10 | 2018-09-04 | 中国科学院合肥物质科学研究院 | 一种氨气检测装置及方法 |
CN114235742A (zh) * | 2021-12-17 | 2022-03-25 | 中国石油大学(华东) | 基于呼吸气体大类标志物复合光谱检测系统及方法 |
Non-Patent Citations (2)
Title |
---|
F. Aousgi et al..On the Pressure and Temperature Dependence of the Absorption Coefficient of NH3.《International Journal of Spectroscopy》.2011,第2011卷1-7. * |
谢秋菊等.猪舍氨气与二氧化碳浓度变化时序预测模型优化.《农业机械学报》.2023,第54卷(第7期),381-390. * |
Also Published As
Publication number | Publication date |
---|---|
CN117347297A (zh) | 2024-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Notholt et al. | Formation of HNO2 on aerosol surfaces during foggy periods in the presence of NO and NO2 | |
Blaser et al. | Process analytical chemistry | |
Levy et al. | Measurements of nitrous acid (HONO) using ion drift-chemical ionization mass spectrometry during the 2009 SHARP field campaign | |
Liu et al. | In situ monitoring of atmospheric nitrous acid based on multi-pumping flow system and liquid waveguide capillary cell | |
Clemitshaw et al. | A calibrated peroxy radical chemical amplifier for ground‐based tropospheric measurements | |
CN117347297B (zh) | 大气nh3分析系统、方法及可读存储介质 | |
CN201607402U (zh) | 煤质含硫量紫外吸收光谱测量装置 | |
CN110455735B (zh) | 一种一硫化碳吸收截面测量装置及方法 | |
Cleland et al. | Measurement of rotational temperature and dissociation in N2O glow discharges using in situ Fourier transform infrared spectroscopy | |
Wu et al. | Measurement techniques for sulfur trioxide concentration in coal-fired flue gas: a review | |
CN216646260U (zh) | 一种检验红外分析用气室的设备 | |
Huntzicker et al. | Investigation of an ambient interference in the measurement of ozone by ultraviolet absorption photometry | |
Ai et al. | Quantitative measurement model for sulfate in the temperature range of 298.15–343.15 K based on vacuum ultraviolet absorption spectroscopy | |
CN116183537A (zh) | 基于差分消元的抗干扰ndir混合气体检测方法和系统 | |
CN114184566B (zh) | 一种基于紫外吸收光谱法适用于不同温度的硫酸根浓度测量模型及其验证方法 | |
Thatipamala et al. | Spectrophotometric method for high biomass concentration measurements | |
EP4033242A1 (en) | A device and a method for complete carbonaceous aerosol analysis in real time | |
Crawley | Application of non-dispersive infrared (NDIR) spectroscopy to the measurement of atmospheric trace gases | |
CN108896519A (zh) | 双光谱烟气汞分析装置及相应的方法 | |
CN109211816B (zh) | 烟气中含Hg浓度的抗干扰检测方法和检测装置 | |
Krebbers et al. | Ultra-broadband spectroscopy using a 2–11.5 μm IDFG-based supercontinuum source | |
CN112824875A (zh) | 一种紫外差分二氧化氮气体的检测方法 | |
Faguy et al. | Real-time polarization modulation in situ infrared spectroscopy applied to the study of atmospheric corrosion | |
CN117388204A (zh) | 一氧化氮气体分析系统、方法及计算机可读存储介质 | |
Ming-Quan et al. | Increasing the accuracy of atmospheric ammonia concentrations calculated from Open-Path Fourier transform infrared spectra using partial least squares model by scanning and removing interference spectral data |
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 |