一种多组分气体分析方法与流程

1.本发明涉及环境监测技术领域，具体为一种多组分气体分析方法。


背景技术：

2.针对应急监测过程中现有分析方法时效性差、分析效率低、选择性差、抗干扰能力差、检测成本高等问题，采用紫外-可见连续光谱法与现代计量学方法相结合开发的现场自动监测的便携式紫外差分吸收光谱仪，实现vocs、nox等多参数同时测量。
3.紫外-可见连续光谱分析技术中常采用紫外差分吸收方法，从采集的光谱数据建立数学模型计算得出气体浓度。现有常用的紫外差分吸收气体浓度计算方法主要包括以下几类：
①
采用多项式拟合得到吸收光谱慢变部分，扣除后得到差分吸收光谱，然后采用特征谱段的累积差分吸光度进行最小二乘拟合得到浓度反演公式；
②
对差分吸收光谱采用多元线性回归方法进行分解求解浓度；
③
采用傅里叶变换方法提取高频部分反变换得到差分吸收光谱，然后再进行累积吸光度最小二乘拟合或者多元线性回归求解浓度。
4.但前述第
①
种方法由于采用累积差分吸光度计算浓度，不能有效利用光谱形状信息，无检测气体时引入的光谱噪声的计算值也总为正值，因此抬高了检测限；第
②
种方法充分利用了光谱形状信息，能有效避免噪声干扰和气体间交叉干扰，但对光谱的波长漂移非常敏感，一两个像素的光谱漂移能使反演精度下降较多；第
③
种方法同样不能避免前述两种方法所遇到的检测下限高或对波长漂移敏感的问题；此外，多组分气体间的非线性交叉干扰是气体分析的难点问题。


技术实现要素：

5.针对上述情况，为克服现有技术的缺陷，本发明提供一种多组分气体分析方法，有效的解决了上述技术背景中提到的问题。
6.为实现上述目的，本发明提供如下技术方案：本发明包括获取光谱并计算吸光度、提取特征值并进行线性交叉干扰消除、二级非线性交叉干扰消除和特征值转化为浓度值，所述获取光谱并计算吸光度：采集待测气体光谱，利用特征谱线将采集的待测气体光谱和参考光谱的光谱波长对齐；
7.计算待测气体光谱的气体吸光度曲线，计算公式如下：
[0008][0009]
其中，a为待测气体光谱的气体吸光度曲线，为参考光谱，为待测气体光谱；
[0010]
所述提取特征值并进行线性交叉干扰消除：是指采用傅里叶变换-多元线性拟合方法提取所述待测气体的吸光度曲线中各种气体的特征值，并消除气体间线性交叉干扰，包括采用傅里叶变换方法得到所述待测气体中so2、no2、no、nh3、苯、甲苯的吸收频谱，选取所述吸收频谱的特征高频吸收频谱进行最小二乘多元线性回归得到特征值；
[0011]
所述二级非线性交叉干扰消除：包括前述消除气体间线性交叉干扰后仍残留的气
体间干扰，采用多项式拟合得到干扰气体特征值与对待测气体特征贡献值的关系，并从待测气体特征值中扣除干扰气体贡献值，得到待测气体最终的特征值；
[0012]
所述特征值转化为浓度值：指建立待测气体特征值与浓度值之间的关系曲线，用于从特征值计算得到浓度值。
[0013]
根据上述技术方案：所述利用特征谱线将采集的待测气体光谱和参考光谱的光谱波长对齐，包括：
[0014]
选择待测气体光谱中至少一个较尖锐且不被待测物吸收或吸收较少的峰作为特征谱线；
[0015]
通过寻找所述至少一个特征谱线，将采集的待测气体光谱和参考光谱的峰值位置进行对齐，多个特征谱线分别针对不同波长范围的光谱进行对齐，从而分段实现整个光谱的对齐。
[0016]
根据上述技术方案：所述采用傅里叶变换方法得到所述待测气体吸收频谱，包括：
[0017]
选取待测气体so2中吸收谱段为280-320nm的吸光度曲线进行傅里叶变换；
[0018]
选取待测气体no2中吸收谱段为420-450nm的吸光度曲线进行傅里叶变换；
[0019]
选取待测气体no中吸收谱段为200-230nm的吸光度曲线进行傅里叶变换；
[0020]
选取待测气体nh3中吸收谱段为200-230nm的吸光度曲线进行傅里叶变换；
[0021]
选取待测气体苯，甲苯中吸收谱段为260-285nm的吸光度曲线进行傅里叶变换。
[0022]
根据上述技术方案：所述吸收频谱的特征高频吸收频谱进行最小二乘多元线性回归得到特征值，包括：将待测气体及干扰气体的特征高频吸收频谱组建成特征矩阵，利用最小二乘多元线性回归方法得到特征值并消除气体间线性交叉干扰。
[0023]
根据上述技术方案：所述步骤特征值转化为浓度值中当计算出的浓度值小于设定的阈值时，根据吸收光谱谱型匹配重新计算浓度值，包括建立吸收光谱谱型特征矩阵，采用最小二乘线性回归计算气体浓度。
[0024]
有益效果：本发明通过获取待测气体的吸光度曲线；采用傅里叶变换-多元线性回归方法计算所述待测气体的特征值并消除气体间线性交叉干扰，对残留的非线性交叉干扰采用多项式拟合的方法消除二级非线性交叉干扰，实现多种气体的同时计算，消除气体间交叉干扰，有效降低检测限，解决波长漂移对反演精度的影响问题。
[0025]
有益效果：。
附图说明
[0026]
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
[0027]
图1是本发明整流程图。
具体实施方式
[0028]
下面结合附图1对本发明的具体实施方式做进一步详细说明。
[0029]
实施例一，由图1给出，本发明提供一种多组分气体分析方法，包括获取光谱并计算吸光度、提取特征值并进行线性交叉干扰消除、二级非线性交叉干扰消除和特征值转化为浓度值，1)所述获取光谱并计算吸光度：采集待测气体光谱，利用特征谱线将采集的待测
气体光谱和参考光谱的光谱波长对齐；
[0030]
计算待测气体光谱的气体吸光度曲线，计算公式如下：
[0031][0032]
其中，a为待测气体光谱的气体吸光度曲线，为参考光谱，为待测气体光谱；
[0033]
2)所述提取特征值并进行线性交叉干扰消除：是指采用傅里叶变换-多元线性拟合方法提取所述待测气体的吸光度曲线中各种气体的特征值，并消除气体间线性交叉干扰，包括采用傅里叶变换方法得到所述待测气体中so2、no2、no、nh3、苯、甲苯的吸收频谱，选取待测气体so2中吸收谱段为280nm的吸光度曲线进行傅里叶变换；选取待测气体no2中吸收谱段为420nm的吸光度曲线进行傅里叶变换；选取待测气体no中吸收谱段为200nm的吸光度曲线进行傅里叶变换；选取待测气体nh3中吸收谱段为200nm的吸光度曲线进行傅里叶变换；选取待测气体苯，甲苯中吸收谱段为260nm的吸光度曲线进行傅里叶变换，选取所述吸收频谱的特征高频吸收频谱进行最小二乘多元线性回归得到特征值；
[0034]
3)所述二级非线性交叉干扰消除：包括前述消除气体间线性交叉干扰后仍残留的气体间干扰，采用多项式拟合得到干扰气体特征值与对待测气体特征贡献值的关系，并从待测气体特征值中扣除干扰气体贡献值，得到待测气体最终的特征值；
[0035]
4)所述特征值转化为浓度值：指建立待测气体特征值与浓度值之间的关系曲线，用于从特征值计算得到浓度值。
[0036]
所述利用特征谱线将采集的待测气体光谱和参考光谱的光谱波长对齐，包括：
[0037]
选择待测气体光谱中至少一个较尖锐且不被待测物吸收或吸收较少的峰作为特征谱线；
[0038]
通过寻找所述至少一个特征谱线，将采集的待测气体光谱和参考光谱的峰值位置进行对齐，多个特征谱线分别针对不同波长范围的光谱进行对齐，从而分段实现整个光谱的对齐。
[0039]
所述吸收频谱的特征高频吸收频谱进行最小二乘多元线性回归得到特征值，包括：将待测气体及干扰气体的特征高频吸收频谱组建成特征矩阵，利用最小二乘多元线性回归方法得到特征值并消除气体间线性交叉干扰。
[0040]
所述步骤特征值转化为浓度值中当计算出的浓度值小于设定的阈值时，根据吸收光谱谱型匹配重新计算浓度值，包括建立吸收光谱谱型特征矩阵，采用最小二乘线性回归计算气体浓度。
[0041]
实施例二，由图1给出，本发明提供一种多组分气体分析方法，包括获取光谱并计算吸光度、提取特征值并进行线性交叉干扰消除、二级非线性交叉干扰消除和特征值转化为浓度值，1)所述获取光谱并计算吸光度：采集待测气体光谱，利用特征谱线将采集的待测气体光谱和参考光谱的光谱波长对齐；
[0042]
计算待测气体光谱的气体吸光度曲线，计算公式如下：
[0043][0044]
其中，a为待测气体光谱的气体吸光度曲线，为参考光谱，为待测气体光谱；
[0045]
2)所述提取特征值并进行线性交叉干扰消除：是指采用傅里叶变换-多元线性拟
合方法提取所述待测气体的吸光度曲线中各种气体的特征值，并消除气体间线性交叉干扰，包括采用傅里叶变换方法得到所述待测气体中so2、no2、no、nh3、苯、甲苯的吸收频谱，选取待测气体so2中吸收谱段为320nm的吸光度曲线进行傅里叶变换；选取待测气体no2中吸收谱段为450nm的吸光度曲线进行傅里叶变换；选取待测气体no中吸收谱段为230nm的吸光度曲线进行傅里叶变换；选取待测气体nh3中吸收谱段为230nm的吸光度曲线进行傅里叶变换；选取待测气体苯，甲苯中吸收谱段为285nm的吸光度曲线进行傅里叶变换，选取所述吸收频谱的特征高频吸收频谱进行最小二乘多元线性回归得到特征值；
[0046]
3)所述二级非线性交叉干扰消除：包括前述消除气体间线性交叉干扰后仍残留的气体间干扰，采用多项式拟合得到干扰气体特征值与对待测气体特征贡献值的关系，并从待测气体特征值中扣除干扰气体贡献值，得到待测气体最终的特征值；
[0047]
4)所述特征值转化为浓度值：指建立待测气体特征值与浓度值之间的关系曲线，用于从特征值计算得到浓度值。
[0048]
所述利用特征谱线将采集的待测气体光谱和参考光谱的光谱波长对齐，包括：
[0049]
选择待测气体光谱中至少一个较尖锐的峰作为特征谱线；
[0050]
通过寻找所述至少一个特征谱线，将采集的待测气体光谱和参考光谱的峰值位置进行对齐，多个特征谱线分别针对不同波长范围的光谱进行对齐，从而分段实现整个光谱的对齐。
[0051]
所述吸收频谱的特征高频吸收频谱进行最小二乘多元线性回归得到特征值，包括：将待测气体及干扰气体的特征高频吸收频谱组建成特征矩阵，利用最小二乘多元线性回归方法得到特征值并消除气体间线性交叉干扰。
[0052]
所述步骤特征值转化为浓度值中当计算出的浓度值小于设定的阈值时，根据吸收光谱谱型匹配重新计算浓度值，包括建立吸收光谱谱型特征矩阵，采用最小二乘线性回归计算气体浓度。
[0053]
有益效果：本发明通过获取待测气体的吸光度曲线；采用傅里叶变换-多元线性回归方法计算所述待测气体的特征值并消除气体间线性交叉干扰，对残留的非线性交叉干扰采用多项式拟合的方法消除二级非线性交叉干扰，实现多种气体的同时计算，消除气体间交叉干扰，有效降低检测限，解决波长漂移对反演精度的影响问题。
[0054]
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。