一种用于复合气体组分的快速检测方法与流程

1.本发明涉及光学传感器领域，特别涉及一种快速高精度复合气体组份、浓度及悬浮颗粒的颗粒度、浓度的检测方法。
技术背景
2.随着社会的发展、科技的进步，人们对于空气质量及安全性的关注度逐年提高。一方面，随着污染的加剧，各种有毒、有害气体，以及颗粒物都使得生活环境中的空气影响到人们的身体健康；另一方面，对于一些特殊生产场景，如矿井、化工厂、建筑场地等，往往充斥着各类有毒气体及粉尘，严重威胁着人们的身体健康，更严重者危及人身安全。因此，对于能够快速、实时、高精度的复合气体组份浓度、颗粒度的检测方法的需求也越来越高。
3.红外吸收谱检测是一种常用的快速、实时、高精度的检测方法，但是这种方法并不能满足目前实际需求，首先红外吸收谱检测方法无法检测气体中悬浮颗粒物的颗粒度及浓度，对于目前空气悬浮颗粒污染严重的城市，以pm2.5为例，虽然中国城市的pm2.5水平在2018年下降了12％之后，在2019年平均下降了9％。但是，98％的城市依旧超过了世界卫生组织(who)的准则值。可见气体中悬浮颗粒物的颗粒度、浓度检测是一项必要检测项。其次，在空气中悬浮物浓度高、或者高粉尘污染的应用场景中，红外吸收谱检测方法的精度将大幅降低。对于后者，目前常采用三种解决方案：第一种，双波长差分法，这种方法使用与测试光波长相近，但被测气体不吸收的另一波长的参考光，认为散射系数相近而认为散射项近似相等，进行差分计算，散射项将直接消去，这种方法的对散射项的近似处理，在红外吸收谱检测方法所采用的长光程测试条件下，带来的测量误差是无法忽略的；第二种，忽略悬浮颗粒物散射对测量精度的影响，即使是在中度空气悬浮物污染的生活环境中，这种解决方案都会使得气体组份浓度检测的精度大大降低，并不能真正解决具体问题；第三种，是在检测气室前设置漉网，这种解决方案虽然能够减少悬浮颗粒物散射对测量精度的影响，但需要定期更换滤网，增加了检测设备运维成本。
4.目前，针对这些问题的解决方案目前主要有文章“张曦雯,李立京,梁生,张春熹.新型多波长红外同步粉尘、气体浓度传感器[j].红外技术,2009,31(01):35
‑
38.”与专利“us2017/0097301a1”及中所提的，不忽略lambert
‑
beer定律
[0005]
i
out
(λ)＝i
in
(λ)e
‑
k(λ)l
ꢀꢀ
(1)
[0006]
衰减系数k(λ)的散射项，将衰减系数还原为吸收项与散射项的求和项。但是这种方法仍存在严重不足。
[0007]
首先，这种方法将n种波长下n组lambert
‑
beer定律的入射光强i
in
(λ)与透射光强i
out
(λ)联立，解出n个衰减系数。为了使得线性方程组有解，需要满足对于任意波长λ，悬浮颗粒物或气体组份的衰减系数k(λ)仅为吸收项或散射项，另一项为零项。但对于红外光谱吸收法这种为了提高测试精度往往采用很长的测试光路的测试方法，吸收项与散射项是不可能同时为零的。这种这使得模型得到大幅简化的方法，会大大降低计算结果的准确度。同时，这种方法在确定悬浮颗粒物颗粒度与浓度时，通过光吸收法仅对入射光强i
in
(λ)与透射
光强i
out
(λ)的测量而得到悬浮颗粒物颗粒度与对应浓度，则需预先确定不同粒径d悬浮颗粒物的衰减系数k
d
(λ)，显然，在实际快速测量中是不可行的。
[0008]
其次，正如上述文章与专利中提及的，使用红外光谱法针对多气体组份的测量，往往使用多种波长的红外光源进行测试。依靠传统的测试方法若想保证各个波长信号间相互无干扰，建立相互独立、无耦合的lambert
‑
beer定律的线性方程组，在共用测试气室的情况下，则需要设置与红光光源数量相同的数量的带有窄带滤波片的探测器。使得这种测量方法对于装置的要求较高，提高了测试成本。在非共用气室的情况下，有着与共用气室同样的复杂装置要求外，还引入了不同气室间的误差，降低了测试精度。
[0009]
再者，对于不同的气体组份，其摩尔分子吸收系数往往差别较大，以二氧化碳、甲烷、甲醛为例，其摩尔分子吸收系数量级分别为10
‑
18
，10
‑
19
，10
‑
20
，在同样光路条件下进行红外光谱法测量，采用相同分辨率的光强探测方法，其摩尔分子吸收系数量级越小其浓度测量误差越大，这将大大增加小摩尔分子吸收系数的气体组份的浓度测量误差。


技术实现要素：

[0010]
本发明为了克服现有技术的不足，提供了一种快速高精度复合气体组份、浓度及悬浮颗粒的颗粒度、浓度的检测方法。
[0011]
具体而言，本发明提供一种用于复合气体组分的快速检测方法，其特征在于，包括以下步骤：
[0012]
(a)以脉冲序列的形式，将多个波长的测量光序列依次射入测试气体气室，所述测量光序列包含多种用于散射光谱法测量的散射光谱测量光和多种用于红外光谱法测量的红外光谱测量光，每种测量光对应一种物质的指纹波长，所述物质包括悬浮颗粒物和各种气体组分；
[0013]
(b)在预定散射角度测量每种波长下的散射光强度，根据散射光谱测量光的测量光强度和散射光强度，计算所述气室内物质对于至少一种波长的散射系数作为参考散射系数；
[0014]
(c)对(b)中的参考散射系数进行拟合，确定各个散射光谱测量光的波长下，任意两种波长的散射系数关系，构建同质散射系数方程；
[0015]
(d)测量红外光谱测量光的输出光强度，基于各波长的入射光强和出射光强构建红外光谱测量方程组；
[0016]
(e)求解由同质散射系数方程以及红外光谱测量方程组构成的联立方程组，确定各个气体组份浓度。
[0017]
在一种优选实现方式中，所述方法还包括基于各个波长(步骤a中的测量光序列)的散射光强度和入射光强度，利用光散射法计算悬浮颗粒物颗粒度及相应浓度。
[0018]
在另一种优选实现方式中，所述步骤(a)中，脉冲序列重复频率为r，脉冲宽度为τ，且对于每个波长，满足单个脉冲在气室内的耗散时间
[0019]
在另一种优选实现方式中，所述步骤(b)包括选取散射系数较大的组分对应(比如大于预定值或从大到小选取)的波长作为计算参考散射系数的测试光，所述参考散射系数
通过前向散射光强测量得到，其计算公式为其中，i
s||
为前向散射光的光强，i0为入射光的光强，l为散射系数测量时的光程。
[0020]
在另一种优选实现方式中，所述红外光谱测量方程组的表达式为其中，i
in
为输入光强，i
out
为输出光强，气体浓度为c1、c2、c3···
c
n
，摩尔分子吸收系数分别为，摩尔分子吸收系数分别为测试波长为测试波长λ1、λ2、λ3···
λ
n
，l为红外光谱测量时的光程。
[0021]
在另一种优选实现方式中，气体组分的浓度计算公式为，
[0022]
在另一种优选实现方式中，在所述步骤(d)中，红外光谱法所使用的入射光波长与步骤(a)中的入射光波长不同，所述步骤(e)包括对于红外光谱测量方程组中的每个方程，以与该方程所对应的波长最为接近的参考散射系数及同质散射系数方程与其联立进行求解，得到相应的散射项。
[0023]
在另一种优选实现方式中，所述步骤(e)中，为红外光谱测量方程组进行不同气体组份浓度求解时，采用摩尔分子吸收系数自大至小排列的每两两相邻或相近的两种气体组份进行阶梯差分运算。
[0024]
在另一种优选实现方式中，所述步骤(e)中，对于摩尔分子吸收系数低于预定值的气体组份浓度，由已测得的高摩尔分子吸收系数的气体组份浓度差分运算获得，其运算公式为
[0025]
其中，
[0026]
c为待测组分的百分比浓度；a为待测组分对红外光线的摩尔分子吸收系数。
[0027]
在另一种优选实现方式中，所述步骤(a)和步骤(d)中的每一种入射光可以为脉冲光或者一段连续光。
[0028]
需要说明的是，虽然下文实施例中，本发明以pm2.5和pm10，co、co2以及甲醛等作为示例进行的描述，但是本领域技术人员应该理解，本发明可以用于对各种其他有害或无害气体进行检测，包括但不限于氨、臭氧、二氧化氮、二氧化硫等，并且本发明还可以对其他粒径的颗粒物进行检测，比如，直径为0.5微米、1微米、5微米、8微米、15微米等等。
[0029]
本发明不仅可以应用于空气污染的检测，还可以用于对特定气体的浓度测定，比如对空气中含氧量的测定，对麻醉气体浓度的测定。本发明还可以应用于燃气、石油、化工、冶金等存在易燃、易爆、毒性气体的危险场所中对可燃气体浓度的测定。本发明还可以用于
对粉尘等物质浓度的测定。
[0030]
发明原理
[0031]
下面对本发明的原理进行详细说明。本发明首先以重复频率为r，脉冲宽度为τ的脉冲序列的方式，向测试气室中发射若干波长的光脉冲进行测试。此次向气室内发射的光既有适用于光散射法测量的测量光又有适用于红外光谱法测量的测量光。单个波长为λ的脉冲在气室内的耗散时间t满足
[0032]
对于用于光散射法测量的测量光，在气室侧部或者前部进行散射光采集。优选地，各个波长中至少一个波长对应pm2.5和pm10的指纹波长。
[0033]
对于红外光谱测量，也需要考虑散射问题。根据lambert
‑
beer定律，将输出光中衰减系数k(λ)写作分立的吸收项与散射项，则有
[0034]
i
out
(λ)＝i
in
(λ)e
‑
[k(λ(+γ(λ)]l
ꢀꢀ
(2)
[0035]
其中k(λ)为吸收系数，γ(λ)为散射系数，l为测试光程。红外吸收谱检测中，特定空气组份吸收相应的检测波长为λ的检测光，对波长为λ的检测光的光强带来的改变，其lambert
‑
beer定律可以写作
[0036][0037]
其中c为气体组份浓度，a
m
为检测波长下该气体组份的摩尔分子吸收系数，l为测试光程。根据式(3)，确定了检测波长为λ下相应的散射项的大小，即可消除悬浮颗粒物对波长为λ的检测光的散射对测量结果带来的影响，得到准确的空气组份浓度。在这里，散射项主要由空气中的悬浮颗粒物产生。对于悬浮颗粒物，其直径远大于红外吸收谱检测波长，散射适用mie散射模型。在任意直角坐标系中，假设光强为i0的入射光沿z轴正向射入，电矢量沿x轴方向，在距散射颗粒r位置处，散射角为θ，与z轴组成的平面为散射面，为入射光振动面与散射面之前的夹角，其散射光强为i
s
。其mie散射公式为
[0038][0039]
同时，为保证散射的单次不相干，使光散射法测量的测试光程尽量短，比如，短于一定值(该值本领域技术人员可以根据经验或者实验进行设定)。极短的测试光程，气体对于测试光的吸收可以忽略不计。经散射后，前向散射光的光强i
s||
可以与无吸收项的lambert
‑
beer定律i
out
(λ)＝i
in
(λ)e
‑
γ(λ)l
相对应。则可以得到，经过长度为l的散射介质后，对于不同波长下不同散射系数间关系为：
[0040][0041]
即为同质散射系数方程，λ
a
、λ
b
为两种测量波长。根据此公式，相同散射介质条件下，两两一组，将用于红外光谱法测量的不同波长测试光分别代入该公式，得到多个同质散射系数方程，构成同质散射系数方程组。
[0042]
因此，本发明为确保各个波长信号间相互无干扰，以重复频率为r，脉冲宽度为τ的脉冲序列的方式，向测试气室中发射若干波长的光脉冲，进行测试。这样使得单个波长为λ的脉冲拥有更高的峰值功率，在探测过程中拥有更高的对比度。同时只要确保单个波长为λ
的脉冲在气室内的耗散时间t满足
[0043][0044]
即可确保多种测试光脉冲，在经过气室后进行测量时相互独立。
[0045]
首先，使用光散射法进行悬浮颗粒物颗粒度及相应浓度的测量。通过测量前向散射光强，用光散射法获得悬浮颗粒物颗粒度和浓度信息。
[0046]
优选地，同时测量特定散射角下所选取的多波长的侧向散射光强i
s
⊥
(λ)，基于侧向散射光强i
s
⊥
(λ)和相应入射光强，使用光散射法测量得到悬浮颗粒物颗粒度及相应浓度。
[0047]
悬浮颗粒物颗粒度及相应浓度的测量可以采用现有方法进行，比如采用题为“基于mie散射的可吸入颗粒物浓度检测系统中的关键技术”，张棚，太原科技大学的硕士学位论文中记载的方式进行即可。
[0048]
接下来讨论脉冲序列中用于红外光谱检测法测量的测量光，即，检测这些测量光的输出光强度。在红外光谱检测中，对于需要测量的气体组份1、气体组份2、气体组份3
···
气体组份m，有对应测试波长λ1、λ2、λ3···
λ
n
，对应摩尔分子吸收系数分别为，对应摩尔分子吸收系数分别为对应所需测得浓度为c1、c2、c3···
c
n
，且测试波长仅与对应气体组份1、气体组份2、气体组份3
···
气体组份m存在较强吸收，其余组份的吸收可以忽略。检测光程为l(需要说明的是，红外光谱测量时，光程可以与粒子浓度测量时的相等也可以不相等)，得到红外光谱测量方程组
[0049][0050]
从多种波长的测试光中选取用于光谱测量的测试光，若波长为λ，出射光强为i0(λ)，测试光程长度为l，并测量前向散射光强i
s||
(λ)，则得到参考散射系数γ(λ)，其值为
[0051][0052]
进一步地，联立(5)(7)式即可推导出，在相同悬浮颗粒物状态下，对应红外吸收光谱检测中其他检测波长λ下的散射系数，选择与测试波长相近波长的参考系数通过同质散射系数方程进行散射项求解，确定(8)式中每一项对应散射项的大小，即可消除悬浮颗粒物对波长为λ的检测光的散射对测量结果带来的影响。
[0053]
在确定了散射项大小的红外光谱测量方程求解过程中，为了提高气体组份浓度的精度，尤其对于摩尔分子吸收系数小的气体。将各个气体组分按照摩尔分子吸收系数自大至小排列，摩尔分子吸收系数大于设定值(本领域技术人员可以根据具体测量精度进行设定)的气体组分，由与其对应的确定了散射项的红外光谱测量方程直接得到其气体组分浓度，其余的两两相邻或相近的气体组份间进行阶梯差分得到准确的气体组份浓度的结果，其运算公式如下
[0054][0055]
这样，就可以确定每个目标气体组分的浓度。
[0056]
需要指出的是，可以引入待测气体组分不敏感(无吸收)波长的参考光对整个系统进行标定和优化。
[0057]
与现有技术相比，本发明具有的有益效果为：
[0058]
1、本发明提供了一种快速高精度复合气体组份、浓度及悬浮颗粒的颗粒度、浓度的检测方法，本发明方法能够消除散射对红外光谱吸收法测量带来的误差，使得测量结果更加准确。尤其在具有悬浮颗粒物污染的环境中优点更为突出。
[0059]
2、本发明以脉冲序列的方式进行多波长测量，能够实现红外光谱法检测的多波长测量的快速、互不干扰的进行。使得测量过程更加快速，同时降低了对测量设备的要求。
附图说明
[0060]
图1为本发明方法的示意性流程图。
具体实施方式
[0061]
为使发明的目的、技术方案和优点更加清楚，下面将结合发明中的附图，对发明中的技术方案进行清楚地描述，显然，所描述的实施例是发明一部分实施例，而不是全部的实施例。基于发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于发明保护的范围。
[0062]
实施例
[0063]
如图1所示，本实施例的快速高精度复合气体组份、浓度及悬浮颗粒的颗粒度、浓度的检测方法，可以包括以下步骤：
[0064]
步骤a：以重复频率为r，脉冲宽度为τ的脉冲序列的方式，向测试气室中发射多个波长的光脉冲，进行测试。单个波长为λ的脉冲在气室内的耗散时间t满
[0065]
步骤b：使用光散射法进行悬浮颗粒物颗粒度及相应浓度的测量，任意从多个波长的测试光中选取适用于散射光谱法测量的测试光，确定选定波长的散射系数。若测试光波长为λ，出射光强为i0(λ)，测试光程长度为l，并测量前向散射光强i
s||
(λ)，得到多个参考散射系数(每个波长有一个参考散射系数)
[0066][0067]
将所获得的参考散射系数代入式(5)。
[0068]
如发明原理中所述，公式(5)为同时得到的特定粒径的目标颗粒物的同质散射系数方程
[0069][0070]
根据此公式，相同散射介质条件下，两两一组，将不同波长测试光分别代入该公
式，得到多个同质散射系数方程，构成同质散射系数方程组。
[0071]
同时测量特定散射角下所选取的多波长的侧向散射光强i
s
⊥
(λ)，按照上述公式(4)使用光散射法分析得到悬浮颗粒物颗粒度及相应浓度。
[0072]
步骤c：按照上述公式(8)使用红外光谱检测法进行测量，这些测量光可以直接在脉冲序列内与用于光谱法测量的光脉冲顺序发射接收，也可以单独发射。在红外光谱检测中，对于需要测量的气体组份1、气体组份2、气体组份3
···
气体组份m，采用对应测试波长λ1、λ2、λ3···
λ
n
(测试波长中至少包含对应每一种组分的吸收峰的光波长)，对应摩尔分子吸收系数分别为对应所需测得浓度为c1、c2、c3···
c
n
，且测试波长λ1、λ2、λ3···
λ
n
仅与对应气体组份1、气体组份2、气体组份3
···
气体组份n存在较强吸收，其余组份的吸收可以忽略。检测光程为l，得到红外光谱测量方程组
[0073][0074]
步骤d：将(5)(7)(8)联立(将公式5和7带入公式8)，在求解过程中，由于公式(7)中可以获得多个参考散射系数，在选取时，选择参考散射系数与测试波长相近的波长的参考散射系数，确定散射项大小。
[0075]
对于各个气体组分，摩尔分子吸收系数较大的气体由前述确定了散射项的红外光谱测量方程直接得到，得到了摩尔分子吸收系数较大的气体组分浓度之后，带回到红外光谱测量方程组，对摩尔分子吸收系数较小的，由确定了散射项大小的红外光谱测量方程组采用阶梯差分运算得到准确的气体组份浓度结果，其运算公式如下
[0076][0077]
虽然上面结合本发明的优选实施例对本发明的原理进行了详细的描述，本领域技术人员应该理解，上述实施例仅仅是对本发明的示意性实现方式的解释，并非对本发明包含范围的限定。实施例中的细节并不构成对本发明范围的限制，在不背离本发明的精神和范围的情况下，任何基于本发明技术方案的等效变换、简单替换等显而易见的改变，均落在本发明保护范围之内。