室内空气颗粒物CCD-激光监测系统及方法与流程

本发明涉及一种室内空气颗粒物垂直剖面和动态变化的测定装置，具体涉及利用ccd相机收集垂直激光光束照射下室内空气颗粒物的散射信号，在线动态获取颗粒物浓度的垂直廓线及动态变化过程。

背景技术：

空气中颗粒物由极小的颗粒和液滴组成，根据粒径可以分为pm2.5和pm10，其意义分别指粒径小于2.5微米和10微米的颗粒物^[1]。其中pm2.5由于粒径极小，可以直接渗透进入细支气管和肺部，对人体健康造成巨大的危害^[2]。因此，pm2.5的长期暴露将增加人类的发病率和死亡率^[3；4]。平均而言，人群在室内停留的时间(22小时以上)远远高于在室外停留的时间，因此，室内空气中颗粒物对人群暴露的重要性远远高于室外。

目前室内空气检测的规范做法是在1.5m高度设置测定仪器或者采样设备。这是基于成年人的平均呼吸高度(口、鼻的位置)大致在1.5m。然而，这显然是针对站姿而言的。事实上，坐姿的呼吸高度通常低于站姿，且很长的睡眠时间内，呼吸高度更低。此外，针对成年人的设计显然不适合儿童或者婴儿。可见，室内空气颗粒物(pm)的垂直分布非常重要。只有在了解pm垂直分布特征的条件下，才可能对包括未成年人在内人群对pm的暴露做出可靠的评估。由于传统方法的固定高度测定，且没有可以获得室内pm浓度垂直分布的有效手段，目前这方面数据为空白。从而导致室内暴露贡献的高不确定性。

室外颗粒物垂直剖面测定技术的传统方法是后向散射激光雷达，传统后向散射激光雷达在接收端附近有一段盲区，同时过渡区较长^[5]，尽管过渡区可以进行几何因子校正^[6]，但是会带入误差，想获得高分辨能力的数据可以通过探空实验，但其高成本又令人望而却步。

近年来有学者对传统激光雷达进行改进，研究出侧向散射激光雷达。侧向散射激光雷达可以很大程度克服传统激光雷达的缺陷。但是其也存在明显的缺点：

(1)在用ccd相机采集激光照片的时候，由于云层以及阳光的影响，照片中拍摄到的激光最高点与实际最高点有一定的偏差，给反演计算的过程带来巨大的误差；

(2)室外侧向激光雷达无法在白天测定，只能在夜晚使用；

(3)室外温度变化将会影响滤光片的透过率曲线，进而给计算带来巨大的误差，甚至可能直接导致错误的计算结果；

(4)室外侧向激光雷达的校正可以通过探空实验或者后向散射激光雷达的结果校正，但是二者的成本令人望而生畏，目前市面上侧向散射激光雷达的价格是30万，而后向散射激光雷达则达到了100万；

(5)室外激光雷达的测量范围在千米级别，但其相应的缺点是在距离校正的过程会引入很大的误差。其次，侧向激光雷达仅仅是名义上的侧向，由于高度的原因，其根本不可能达到真正的侧面拍摄的效果。

因为成本、激光强度以及需要后向散射雷达校正的原因将室外激光雷达直接在室内使用是无法实现的，因此需要重新设计了一套室内空气颗粒物ccd-激光监测系统。

技术实现要素：

本发明的目的在于提供一种ccd-激光监测系统，利用激光器发射垂直激光束，使用ccd相机从侧面记录空气中颗粒物散射信号，在进行背景校正后获得颗粒物浓度垂直分布的廓线信息，再利用不同高度的pm测定装置校正后获得室内空气pm(特定粒径)浓度的垂直分布数据。

本发明首先提供了一种室内空气颗粒物ccd-激光监测系统，包括一个ccd相机，至少一个激光器，以及分布在不同高度的一系列pm在线监测仪，其中，所述一系列pm在线监测仪包括至少3个pm在线监测仪，它们垂直分布在室内空间的不同高度处；所述激光器放置在该系列pm在线监测仪前的地面上，垂直向上发射激光；所述ccd相机用于拍摄所述激光器发射的激光的影像，通过激光影像的灰度值得到室内空气颗粒物浓度的垂直分布情况。

在本发明的室内空气颗粒物ccd-激光监测系统中，垂直分布的所述一系列pm在线监测仪的个数优选为3-10个，更优选为5-9个。所述激光器可以是一个，也可以是多个，以获取平行数据。这多个激光器一字排开放置在pm在线监测仪前，多个激光器所在直线与pm在线监测仪在水平方向上的距离优选为100mm～300mm之间。激光器发射的激光波长因所测定的室内空气颗粒物的种类而定，当所述室内空气颗粒物ccd-激光监测系统测定的是pm2.5浓度时，所述激光器发射波长为532nm的激光，所述pm在线监测仪为pm2.5在线监测仪。

所述ccd相机的镜头面对所述激光器发射的激光，设ccd相机到激光束的垂直距离为d，ccd相机位于距离地面高度为h的位置，根据室内空间的高度设定d和h的数值，保证ccd相机能将整个激光束拍到即可，通常h设定为室内空间高度的一半左右。

参见图1，下面以监测pm2.5的浓度为例，说明本发明的室内空气颗粒物ccd-激光监测系统的实施原理：

为了探测pm2.5浓度与激光能量的关系，气溶胶颗粒的直径与发射激光的波长比大于0.1时，可以使用mie散射理论来解释，因此采用532nm的激光来探测pm2.5的浓度满足mie散射理论的条件。

参考室外ccd-激光气溶胶探测系统的激光雷达方程为^[7]：

e(z,θ)＝kelspτzτrβ(θ,z)dz/r²(式1)

其中，e(z,θ)为高度z处散射角为θ时ccd接收到的侧向散射光能量，el为曝光时间内激光发射器发射的总的激光能量，sp为每个像元的接收面积，τz是激光从发射点到高度z的透过率，τr为激光从高度z处被散射后经过距离r到达ccd处的透过率，β(θ,z)是在散射角θ和高度z上的散射函数；这里dz是散射角宽度dθ对应的长度，k为系统常数。

由几何关系可以知道：

tanθ＝tan(π-θ)＝-tanα＝-d/(z-h)

z＝-d/tanθ+h(式2)

上式中的h表示ccd相机离地面的高度，α表示激光上某点与ccd相机的连线与激光的夹角，对上式两边微分，可得

dz＝-d(-sec²θdθ/tan²θ)＝ddθ/sin²θ＝ddθ/(d²/r²)

dz＝r²dθ/d(式3)

式3中，r为激光顶部到相机的距离，d为ccd相机到出射激光束的垂直距离。

因为在室内中ccd相机距离激光的距离比较小，因此散射光被ccd接收时在通过路径上的衰减可以忽略不记，所以把透射率简化为1，上述式1简化为：

e＝kelspβ(θ)dθ(式4)

单位质量颗粒物侧向角散射系数可以表示为^[8]：

式中，n(α)是α粒径气溶胶颗粒谱分布，βaer(θ,α)为α粒径的气溶胶颗粒角散射系数，α1、α2分别为气溶胶最小和最大粒径，βatm(θ)为分子角散射系数，natm为单位提及的分子数(可近似为常数)。实验中将实际测量水平区域10cm内颗粒物等效为均匀混合的气溶胶粒子，上式近似为：

β(θ)＝βaer(θ)naer+βatm(θ)natm(式6)

单个气溶胶粒子的角散射系数为：

式7中，为气溶胶颗粒所有方向的角散射截面平均值，是一常数；p(θ)为气溶胶散射相函数。

设pm2.5的浓度为m，大气成分组成比例在实验测量期间未发生变化，单位质量pm2.5颗粒的比例系数为c，则气溶胶单位体积的颗粒数为^[22]：

式8中，ρ为气溶胶颗粒密度，所以e(θ)可以表示为：

因此可以发现，ccd接收的光能量与pm2.5浓度呈线性关系，ccd接收的光能量与接收到的电压成正比，也与灰度值成正比，所以ccd像元灰度值与pm2.5浓度呈线性相关。

基于上述原理，本发明提供了一种测定室内空气颗粒物浓度的垂直分布信息的方法，包括如下步骤：

1)在室内设置上述室内空气颗粒物ccd-激光监测系统，在开启激光器之前，通过ccd相机获取无激光背景影像；

2)开启激光器，垂直向上发射激光，用ccd相机获取激光影像，并利用步骤1)获取的无激光背景影像数据进行校正；

3)从步骤2)校正后的影像中提取激光光路中不同高度处的灰度值；

4)在步骤2)ccd相机拍摄时，采用pm在线监测仪在不同高度同步获得实测空气颗粒物浓度数据；

5)根据不同高度和不同时间获取的影像灰度值和步骤4)实测的空气颗粒物浓度数据建立校验模型，并运用二者之间的回归关系获得室内空气颗粒物浓度的垂直分布信息。

上述步骤2)中，所述激光器向上发射波长为532nm的激光；步骤4)中采用pm2.5在线监测仪在不同高度同步获得实测pm2.5浓度数据；在步骤5)得到室内空气pm2.5浓度的垂直分布信息。

优选的，上述步骤3)运用matlab软件从影像中提取激光束所经过的区域的灰度值，取同一高度光束区域灰度值的平均值，得到该高度处的灰度值。

进一步的，在上述测定过程中可以采用多个激光器同时获得平行数据。

本发明的技术优势主要体现在：

目前对于室内空气检测的传统方法是固定高度放置采样器监测，市面上没有对室内空气颗粒物垂直廓线测定的仪器。本发明提供的室内空气颗粒物ccd-激光监测系统弥补了室内空气垂直廓线测定的空白，能够提供颗粒物垂直分布廓线信息，对于未成年人颗粒物暴露测定具有重要意义。

附图说明

图1为本发明的室内空气颗粒物ccd-激光监测系统的结构示意图，其中：1-ccd相机，2-激光器，3-pm在线监测仪。

图2为实施例中ccd相机拍到的激光灰度和pm2.5在线监测仪测定浓度的校准曲线。

图3为实施例中获得的室内空气pm2.5垂直分布廓线。

具体实施方式

下面结合附图，通过实施例进一步描述本发明，但不以任何方式限制本发明的范围。

如图1所示，所述室内空气颗粒物ccd-激光监测系统包括一个ccd相机1，至少一个激光器2，以及分布在不同高度的一系列pm在线监测仪3，其中：所述一系列pm在线监测仪3包括3-10个，优选为5-9个，垂直分布在室内空间的不同高度处；所述激光器2放置在该系列pm在线监测仪3前的地面上，垂直向上发射激光；所述ccd相机镜头面对激光束，到激光束的垂直距离为d，ccd相机位于距离底面高度为h的位置；室内空间的高度记为z，根据房间高度设定h和d的数值，保证ccd相机能将整个激光束拍到即可，通常h设定为房间高度的一半左右即可。

在本发明的一个实施例中，如图1所示，在垂直方向上分布了9个pm2.5在线监测仪3，从下到上分别按第1号～第9号编号，分别位于0、165mm、373mm、631mm、874mm、1149mm、1394mm、1831mm、2373mm、2788mm高度处。设置了3个波长为532nm的激光器2，并排放置在该系列pm2.5在线监测仪3前面，ccd相机与pm2.5在线检测仪在水平方向上的距离为4642mm，相邻激光器2之间的间隔为100mm，激光器2到天花板的距离z＝2893mm。ccd相机1位于高度为h＝1580mm的地方，ccd相机1的镜头到激光束的垂直距离为d＝4205mm。

测定室内空气pm2.5浓度的垂直分布数据的步骤如下：

(1)如图1所示，将ccd相机1(1024*748像素)置于距离垂直激光束4205mm处，在开启激光器之前获取背景影像；

(2)开启激光器2，垂直向上发射波长为532nm的激光，用ccd相机1获取影像，并利用步骤(1)获取的无激光背景数据进行校正。曝光时间低于2s。拍摄时长可根据需求设置。短时间拍摄可以计算均值获得静态数据，长时期拍摄则可以得到动态变化信息。

(3)从ccd相机1获得影像中提取激光光路中不同高度点的灰度值。具体运用matlab提取光束所经过的位置的灰度值，取同一高度光束区域灰度值的平均值，得到不同高度的灰度值。

(4)拍摄期间采用pm2.5在线监测仪3在不同高度同步获得实测pm2.5浓度数据，根据不同高度和不同时间获取的影像灰度值和实测pm2.5浓度数据，建立校验模型,如图2所示，可以得到y＝15.444x+22.159，其中x代表激光影像中的灰度值，y代表pm2.5浓度，从而将灰度值转换为pm2.5浓度值。并最终运用pm2.5在线监测仪的数据和影像灰度之间的回归关系获得室内空气中pm2.5垂直分布廓线的完整信息，如图3所示。测定过程可以采用多个激光器同时获得平行数据。

上述步骤(4)中，在建立校验模型时，首先取对应pm2.5在线检测仪高度的激光影像灰度值(因为pm2.5在线检测仪本身具有一定体积，所以取激光光路上同一高度±30mm范围内的灰度平均值作为其对应灰度值)，然后将灰度值归一化处理，再将得到的九个点的pm2.5浓度值和归一化后的灰度值作校准曲线；通过校准曲线将激光通过路径上的灰度值全部转化为对应的pm2.5浓度值，自此可以得到pm2.5垂直分布浓度廓线。

需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域的技术人员可以理解：在不脱离本发明及所附权利要求的精神和范围内，各种替换和修改都是可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权利要求书界定的范围为准。

参考文献：

[1]liuc,hsupc,leehw,etal.transparentairfilterforhigh-efficiencypm2.5capture[j].naturecommunications,2015,6:6205.

[2]poperc,dockerydw.healtheffectsoffineparticulateairpollution:linesthatconnect[j].airrepair,2006,56(6):707-708.

[3]wus,dengf,weih,etal.associationofcardiopulmonaryhealtheffectswithsource-appointedambientfineparticulateinbeijing,china:acombinedanalysisfromthehealthyvolunteernaturalrelocation(hvnr)study[j].environmentalscience&technology,2014,48(6):3438-48

[4]brookrd,rajagopalans,popeca,etal.particulatematterairpollutionandcardiovasculardisease[j].circulation,2010,121(21):2331-2378.

[5]孟祥谦,胡顺星,王珍珠,等.ccd激光雷达探测边界层气溶胶消光系数垂直分布[j].光学学报,2013,33(8):14-20.

[6]汪少林,曹开法,胡顺星,等.对激光雷达几何因子的分析与测量[j].激光技术,2008,32(2):147-150.

[7]孟祥谦,胡顺星,王英俭,等.基于电荷耦合器件探测气溶胶散射相函数与大气能见度的研究[j].光学学报,2012,32(9):1-6.

[8]胡淼,谢家亮,吴端法,等.基于ccd后向散射激光信号的pm2.5测量研究[j].光学学报,2015,35(2):8-14.