CN105547944B - A kind of inhalable particles object detecting method based on stacking diffraction - Google Patents
A kind of inhalable particles object detecting method based on stacking diffraction Download PDFInfo
- Publication number
- CN105547944B CN105547944B CN201610012762.2A CN201610012762A CN105547944B CN 105547944 B CN105547944 B CN 105547944B CN 201610012762 A CN201610012762 A CN 201610012762A CN 105547944 B CN105547944 B CN 105547944B
- Authority
- CN
- China
- Prior art keywords
- particles
- laser beam
- inhalable
- air
- particulate matter
- 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.)
- Expired - Fee Related
Links
- 239000002245 particle Substances 0.000 title claims 7
- 238000000034 method Methods 0.000 title claims 3
- 238000001514 detection method Methods 0.000 claims 7
- 238000005516 engineering process Methods 0.000 claims 3
- 230000005540 biological transmission Effects 0.000 claims 2
- 239000013618 particulate matter Substances 0.000 claims 2
- 230000003287 optical effect Effects 0.000 claims 1
- 238000005086 pumping Methods 0.000 claims 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
- G01N15/0606—Investigating concentration of particle suspensions by collecting particles on a support
- G01N15/0612—Optical scan of the deposits
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth 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
本发明公开了一种基于层叠衍射的可吸入颗粒物检测方法,属于可吸入颗粒物检测领域。其方法是利用激光器和光学透镜产生直径为1~10mm,波长为200~1000nm的激光束并将其照向待测空气;激光经过空气中可吸入颗粒物时产生衍射,在激光器对侧形成干涉图样;由于可吸入颗粒物在重力作用下存在着沉降速率,利用帧率不低于600fps的高速电子感光元件可以接收一系列有关这些颗粒物的干涉图样;对于PM2.5、PM10等一般的可吸入颗粒物,这些干涉图样相互间有所叠加,故可利用层叠衍射算法将这些颗粒物全息形貌计算出来;利用实时层叠衍射算法,在显示设备上可以及时显示可吸入颗粒物图像,并经处理器统计得出不同大小颗粒相应的浓度值。该方法无需PM2.5、PM10等切割器,能够对可吸入颗粒物进行准实时三维成像,有助于在线研究可吸入颗粒物的行为。
The invention discloses a method for detecting inhalable particulate matter based on stacked diffraction, and belongs to the field of inhalable particulate matter detection. The method is to use a laser and an optical lens to generate a laser beam with a diameter of 1-10mm and a wavelength of 200-1000nm and shine it on the air to be measured; when the laser passes through inhalable particles in the air, it will diffract and form an interference pattern on the opposite side of the laser ;Due to the sedimentation rate of inhalable particles under the action of gravity, a series of interference patterns related to these particles can be received by using a high-speed electronic photosensitive element with a frame rate of not less than 600fps; for general inhalable particles such as PM2.5 and PM10, These interference patterns are superimposed on each other, so the holographic shape of these particles can be calculated by using the stacked diffraction algorithm; using the real-time stacked diffraction algorithm, the image of the inhalable particles can be displayed on the display device in time, and the difference can be obtained through the statistics of the processor. Concentration values corresponding to large and small particles. This method does not require cutters such as PM2.5 and PM10, and can perform quasi-real-time three-dimensional imaging of inhalable particulate matter, which is helpful for online research on the behavior of inhalable particulate matter.
Description
技术领域technical field
本发明属于可吸入颗粒物检测领域,具体涉及利用层叠衍射成像技术对空气中可吸入颗粒物进行光强度、相位信息再现,实现在线全息成像,并对不同粒径可吸入颗粒物浓度进行计算。The invention belongs to the field of inhalable particulate matter detection, and specifically relates to reproducing light intensity and phase information of inhalable particulate matter in the air by using stacked diffraction imaging technology, realizing online holographic imaging, and calculating the concentration of inhalable particulate matter with different particle diameters.
背景技术Background technique
雾霾已成为严重影响人们生产生活的一种天气现象,可吸入颗粒物是雾霾的核心污染源,其可对人体健康造成严重危害,尤其是直径小于2.5μm的颗粒(PM2.5),由于可以直接进入人体肺部,危害更大。为保障人民群众正常生活,开发能够对可吸入颗粒物,尤其是PM2.5的实时、准确的监测设备显得十分紧迫。当前对可吸入颗粒物的监测方法主要有以下几种。Smog has become a weather phenomenon that seriously affects people's production and life. Inhalable particulate matter is the core pollution source of smog, which can cause serious harm to human health, especially particles with a diameter of less than 2.5 μm (PM2.5). Directly entering the lungs of the human body is even more harmful. In order to ensure the normal life of the people, it is very urgent to develop real-time and accurate monitoring equipment for inhalable particulate matter, especially PM2.5. The current monitoring methods for inhalable particulate matter mainly include the following.
(1) 滤膜称重法,该方法让空气经不同粒径颗粒切割器后,通过事先放置的滤纸,PM2.5、PM10等不同粒径的可吸入颗粒物将收集到相应的滤纸上;经过一定时间后,取出滤纸,通过比较滤纸吸附前后的质量的变化来推测可吸入颗粒物浓度。(1) Filter membrane weighing method, this method allows the air to pass through the pre-placed filter paper after passing through the particle cutters of different particle sizes, and the inhalable particles of different particle sizes such as PM2.5 and PM10 will be collected on the corresponding filter paper; After a certain period of time, the filter paper was taken out, and the concentration of inhalable particulate matter was estimated by comparing the mass change of the filter paper before and after adsorption.
(2) 射线吸收法,该方法采用和滤膜称重法相同的技术获得沉积有颗粒物的滤纸,用β射线照射;β射线穿过滤纸和颗粒物的衰减程度与颗粒物质量成正比,据此根据射线吸收情况计算出颗粒物质量并计算其浓度(农永光 等,价值工程,2013(35),pp.304-305)。(2) Ray absorption method, which adopts the same technology as the filter membrane weighing method to obtain filter paper with deposited particles, and irradiates it with β-rays; The mass of particulate matter is calculated from the ray absorption and its concentration is calculated (Nong Yongguang et al., Value Engineering, 2013 (35), pp.304-305).
(3) 微量振荡天平法,该方法采用锥形空心玻璃管,粗头固定,细头装有滤芯。空气从粗头进,通过相应粒径的切割器,该粒径可吸入颗粒物被截留在滤芯上。在电场作用下,细头以一定频率振荡,振荡频率与细头重量平方根成反比,根据振荡频率的变化,可以算出收集到颗粒物的质量并能根据进气量推算出其浓度(农永光 等,价值工程,2013(35),pp.304-305)。(3) The micro-oscillating balance method, which uses a tapered hollow glass tube, the thick end is fixed, and the thin end is equipped with a filter element. The air enters from the thick head and passes through the cutter of the corresponding particle size, and the inhalable particles of this particle size are trapped on the filter element. Under the action of the electric field, the fine head oscillates at a certain frequency, and the oscillation frequency is inversely proportional to the square root of the fine head weight. According to the change of the oscillation frequency, the mass of the collected particles can be calculated and its concentration can be calculated according to the intake air volume (Nong Yongguang et al., Value Engineering, 2013(35), pp.304-305).
(4) 激光散射法, 该方法基于Mie散射理论,利用激光做光源,采用光电器件接收光强度,根据散射光通量和粒径关系来确定空气中颗粒体密度,并推算出颗粒浓度(葛宝臻等,天津大学学报(自然科学与工程技术版),2013(1),pp.22-28)。(4) Laser scattering method, which is based on Mie scattering theory, uses laser light as the light source, uses photoelectric devices to receive light intensity, and determines the particle density in the air according to the relationship between scattered light flux and particle size, and calculates the particle concentration (Ge Baozhen et al., Journal of Tianjin University (Natural Science and Engineering Technology Edition), 2013 (1), pp.22-28).
上述方法均存在一定的不足:There are certain deficiencies in the above methods:
滤膜称重法、射线吸收法方法采样时间较长,不适合实时性测量的要求;微量天平法受湿度、温度的干扰较大;激光散射法虽然能够做到对可吸入颗粒物的实时测量,但是由于激光散射与颗粒物折射率有关,该方法所得结果准确性欠佳。The filter weighing method and the ray absorption method have a long sampling time, which is not suitable for real-time measurement requirements; the microbalance method is greatly disturbed by humidity and temperature; although the laser scattering method can achieve real-time measurement of inhalable particles, However, due to the relationship between laser light scattering and the refractive index of particles, the accuracy of the results obtained by this method is not good.
上述方法均无法对可吸入颗粒物进行在线成像,对PM2.5、PM10等不同粒径的可吸入颗粒物测量也必须依赖于切割器;这一方面使得难以直观的看到可吸入颗粒的运动行为,另一方面也增加了装置设计的复杂度。None of the above methods can perform online imaging of inhalable particles, and the measurement of inhalable particles of different particle sizes such as PM2.5 and PM10 must also rely on a cutter; this aspect makes it difficult to intuitively see the movement behavior of inhalable particles. On the other hand, it also increases the complexity of device design.
除上述方法外,现在人们也探索使用数字全息术来测量PM2.5等可吸入颗粒物;该方法基于全息技术、CCD、CMOS等光电器件以及计算机技术,使可吸入颗粒物实时成像成为可能。目前用以观测可吸入颗粒物的数字全息术是同轴数字全息术,该方法利用颗粒物边缘的散射光作为参考光来记录可吸入颗粒物的相位信息(王萍,光学仪器,2015(4),pp.287-292)。这种方法虽然可以对颗粒物进行全息成像,但是由于物光和参考光方向相同,全息图频谱信息交叠严重,所得像的分辨率较低。In addition to the above methods, people are now also exploring the use of digital holography to measure PM2.5 and other inhalable particulate matter; this method is based on holographic technology, CCD, CMOS and other photoelectric devices and computer technology, making real-time imaging of inhalable particulate matter possible. The digital holography currently used to observe respirable particles is coaxial digital holography, which uses the scattered light at the edge of the particles as a reference light to record the phase information of respirable particles (Wang Ping, Optical Instrument, 2015 (4), pp .287-292). Although this method can perform holographic imaging of particles, because the object light and the reference light have the same direction, the spectral information of the hologram overlaps seriously, and the resolution of the obtained image is low.
近些年兴起的基于层叠衍射(Ptychography)的显微技术也可以用来记录并显示物体全息图像。该方法首先记录物体一系列衍射图,每个衍射图至少同另外一个衍射图有交叠区域;之后通过迭代计算来得到物体信息,该信息满足所有衍射图的要求,从而再现出物体全息像。该方法可以得到比传统数字全息术分辨率更高的图像,甚至可以超过瑞利判据所确定的光学极限(Marrison et al, Sci. Rep, 2013(3), No. 2369)。一般的层叠衍射显微术通过移动光源和样品的相互位置来得到一系列符合要求的衍射图,通常这借助于微机械装置来实现。而对于空气中PM2.5、PM10等可吸入颗粒物而言,其在空气中并非不动,它们存在着一定的沉降速率,实验结果表明,对于粒径在2μm左右的颗粒其在200μm/s的量级 (冷吟 等,可吸入颗粒物在驻波声场中运动的可视化研究,工程热物理学报,2007(增刊1),pp.209-212);根据斯托克斯沉降公式,沉降速率v和颗粒半径r的平方成正比,那么粒径为10μm的颗粒其沉降速率为5mm/s的量级。对这些颗粒进行衍射并采用拍摄速率不低于600fps的高速光电传感装置对其衍射图进行记录,那么相邻两帧衍射图一般会有交叠区域,并且拍摄速率越高,交叠区域越大,因此这种方法满足层叠衍射算法所需的条件,故可以利用层叠衍射算法对可吸入颗粒物进行再现成像。该过程无需将颗粒取下,因此此成像为在线成像;结合在光场中的定位技术,还可以实现对颗粒物的实时在线监测,通过计算机处理,颗粒将其形貌和行为实时显示出来。Ptychography-based microscopic techniques emerging in recent years can also be used to record and display holographic images of objects. This method first records a series of diffraction patterns of the object, and each diffraction pattern has at least an overlapping area with another diffraction pattern; then, the object information is obtained through iterative calculations, and the information meets the requirements of all diffraction patterns, thereby recreating the holographic image of the object. This method can obtain images with higher resolution than traditional digital holography, even exceeding the optical limit determined by the Rayleigh criterion (Marrison et al, Sci. Rep, 2013(3), No. 2369). The general stacked diffraction microscopy obtains a series of diffraction patterns that meet the requirements by moving the mutual position of the light source and the sample, usually by means of a micromechanical device. For inhalable particles such as PM2.5 and PM10 in the air, they are not motionless in the air, and they have a certain settling rate. Magnitude (Leng Yin et al., Visualization Study on the Movement of Inhalable Particles in the Standing Wave Sound Field, Journal of Engineering Thermophysics, 2007 (Supplement 1), pp.209-212); according to the Stokes settlement formula , the settling velocity v is proportional to the square of the particle radius r , then the settling velocity of a particle with a particle size of 10 μm is on the order of 5 mm/s. Diffraction these particles and use a high-speed photoelectric sensing device with a shooting rate of not less than 600fps to record their diffraction patterns, then two adjacent frames of diffraction patterns generally have overlapping areas, and the higher the shooting rate, the more overlapping areas Therefore, this method satisfies the conditions required by the stacked diffraction algorithm, so the stacked diffraction algorithm can be used to reproduce the imaging of inhalable particles. This process does not need to remove the particles, so this imaging is online imaging; combined with the positioning technology in the light field, real-time online monitoring of particles can also be realized. Through computer processing, the particles can display their shape and behavior in real time.
发明内容Contents of the invention
本发明的目的在于提出一种基于层叠衍射术的可吸入颗粒物检测方法,利用该方法能够对可吸入颗粒物进行全息成像、实时三维再现,并对不同粒径可吸入颗粒物浓度进行计算。该方法无需PM2.5、PM10等切割器,仪器构造简便;且能够对可吸入颗粒物进行准实时三维动态显示,不仅精度高,而且可以在线研究可吸入颗粒物的行为。The purpose of the present invention is to propose a detection method for inhalable particulate matter based on stacked diffraction, which can be used for holographic imaging, real-time three-dimensional reproduction of inhalable particulate matter, and calculation of the concentration of inhalable particulate matter with different particle sizes. This method does not require PM2.5, PM10 and other cutters, and the instrument structure is simple; and it can perform quasi-real-time three-dimensional dynamic display of inhalable particulate matter, which not only has high precision, but also can study the behavior of inhalable particulate matter online.
本发明技术方案如下文所示。The technical solution of the present invention is as follows.
一种基于层叠衍射显微术的可吸入颗粒物检测方法,包括以下步骤:(1)如附图1所示,利用激光器1和光学透镜2产生球面发散激光束并将其照向待测空气,该激光束一般为发散光斑,照射空气时光斑横截面积在1mm2到100mm2的范围内,激光波长为200~1000nm;待测空气一般在含输入、泵出装置以及气压计31的透明气室3中,由输入装置输入该气室;泵出装置可用来排出气室中原有空气,防止残留空气带来测量误差;测量时输入输出装置保持封闭,使气室中空气静止;气室被激光束照射区域的两个面平行,且与激光束中心轴方向垂直,沿该方向气室内空间厚度在0.5mm到1mm之间;此气室距感光元件探头的距离不小于2cm,且和距激光束球心距离之比不小于10。A method for detecting inhalable particulate matter based on stacked diffraction microscopy, comprising the following steps: (1) As shown in Figure 1, a laser 1 and an optical lens 2 are used to generate a spherically diverging laser beam and irradiate it to the air to be tested, The laser beam is generally a divergent spot, the cross-sectional area of the spot is in the range of 1 mm 2 to 100 mm 2 when irradiating air, and the laser wavelength is 200-1000 nm; In chamber 3, the air chamber is input by the input device; the pumping device can be used to discharge the original air in the air chamber to prevent the residual air from causing measurement errors; the input and output devices are kept closed during measurement, so that the air in the air chamber is still; the air chamber is The two surfaces of the laser beam irradiated area are parallel and perpendicular to the central axis of the laser beam. The thickness of the space in the air chamber along this direction is between 0.5mm and 1mm; the distance between the air chamber and the photosensitive element probe is not less than 2cm, and the distance between the air chamber and the The ratio of the distance between the centers of the laser beams is not less than 10.
(2)利用高速CCD等电子感光元件接收通过待测空气的激光束光斑,并将其传输到含层叠衍射再现系统的处理器,传输装置采用有线传播或者蓝牙、Wifi等无线传播技术。(2) Use high-speed CCD and other electronic photosensitive elements to receive the laser beam spot passing through the air to be measured, and transmit it to the processor containing the laminated diffraction reproduction system. The transmission device adopts wired transmission or wireless transmission technology such as Bluetooth and Wifi.
(3)对激光束光斑中所含可吸入颗粒物信息扣除背景信息后进行层叠衍射再现处理;其主要计算过程为:(3) After deducting the background information of the inhalable particulate matter information contained in the laser beam spot, perform lamination diffraction reproduction processing; the main calculation process is:
第一步,根据颗粒在光斑中衍射图样确定所探测颗粒衍射斑在整个光斑中的相对位置,并将不同颗粒的衍射斑分开,根据一系列衍射图样,确定其位置的变化;The first step is to determine the relative position of the diffraction spot of the detected particle in the entire spot according to the diffraction pattern of the particle in the light spot, and separate the diffraction spots of different particles, and determine the change of its position according to a series of diffraction patterns;
第二步,针对第一个颗粒,根据其衍射图的位置确定层叠衍射计算所需入射光波函数的初始值,其中R j 表示第j个衍射图中颗粒相对于其在第1个衍射图中位置的偏移;由于入射波为发散球面激光束,波长、距离气室位置已知,因此该入射波波函数已知,只需知道颗粒物在光斑中位置,此处的入射波函数便能求得,由于气室玻璃的影响,真实波函数与该波函数略有区别,可在随后迭代计算中进行修正;In the second step, for the first particle, determine the incident light wave function required for the stacked diffraction calculation according to the position of its diffraction pattern The initial value of , where R j represents the offset of the particle in the j- th diffraction pattern relative to its position in the first diffraction pattern; since the incident wave is a diverging spherical laser beam, the wavelength and distance from the gas cell are known, so the The incident wave function is known, and the incident wave function here can be obtained only by knowing the position of the particle in the spot. Due to the influence of the gas chamber glass, the real wave function is slightly different from this wave function, which can be calculated iteratively later make corrections in
第三步,层叠衍射计算所需颗粒形状相位分布(即物函数)初始值设定为圆形,透光率一般设为1,相位设为0;The third step is to calculate the phase distribution of the required particle shape (i.e., the matter function ) is initially set to a circle, the light transmittance is generally set to 1, and the phase is set to 0;
第四步,根据入射波波函数和颗粒函数的初始值,得到颗粒入射函数初始值,其中角标0表示初始值;The fourth step is to obtain the initial value of the particle incident function according to the initial value of the incident wave function and the particle function , where the subscript 0 represents the initial value;
第五步,得到更新的颗粒入射函数值;由于颗粒尺寸相对于其到感光元件的距离足够小,故衍射图样可满足夫琅禾费衍射条件,此时衍射图样强度,F{}表示傅里叶变换,操作过程中采用快速傅里叶变换;,由夫琅禾费衍射公式可以计算出此时在感光元件接收端的波函数为;将U 0的相位部分保留,强度部分用I的平方根代替,生成新的接收端波函数,即;对进行傅里叶逆变换,得到新的颗粒入射函数初始值;The fifth step is to obtain the updated particle incident function value; since the particle size is small enough relative to the distance from the photosensitive element, the diffraction pattern can satisfy the Fraunhofer diffraction condition, and the intensity of the diffraction pattern at this time , F {} represents Fourier transform, fast Fourier transform is used in the operation process;, according to the Fraunhofer diffraction formula, the wave function at the receiving end of the photosensitive element at this time can be calculated as ;Retain the phase part of U 0 , replace the intensity part with the square root of I , and generate a new wave function at the receiving end ,Right now ;right Perform inverse Fourier transform to obtain a new initial value of the particle incident function ;
第六步,更新物函数O和入射波函数P,更新所需公式分别为:The sixth step is to update the object function O and the incident wave function P , and the formulas required for the update are:
第七步,比较是否达到精度要求,设定参量RMS 如下The seventh step is to compare whether the accuracy requirements are met, and set the parameter RMS as follows
如果RMS小于设定的误差值Err,则和P new分别为物函数和入射光函数,反之,则将O0和P0用Onew和Pnew代替,重复第四至七步,直至获得符合条件的物函数和入射光函数;Err一般取值为10-3至10-5,其取值越小,精度越高,迭代时间越长;If RMS is less than the set error value Err , then and P new are the object function and incident light function respectively, otherwise, replace O 0 and P 0 with O new and P new , and repeat steps 4 to 7 until the qualified object function and incident light function are obtained; Err is generally The value is 10 -3 to 10 -5 , the smaller the value, the higher the accuracy and the longer the iteration time;
第八步,针对其他感兴趣的颗粒,重复第二至七步,直至将所有感兴趣的颗粒形貌像均获取为止;其中各感兴趣颗粒可以采用并行计算的方法,对各颗粒同时计算。In the eighth step, for other interested particles, repeat the second to seventh steps until all the shape images of the interested particles are obtained; the parallel calculation method can be used for each interested particle to calculate for each particle at the same time.
(4)将所得颗粒的形状相位分布(即物函数)显示出来,即得到其三维图像;根据其三维图像,统计不同粒径可吸入颗粒物浓度,并将再现图像、颗粒物浓度传输到显示设备中,传输装置采用有线传播或蓝牙、Wifi等无线传播技术。结合光斑所照射空气区域的体积可计算出不同粒径颗粒物体浓度,然后根据国标WS/T206-2001将其转化为相应的质量浓度。(4) Display the shape and phase distribution (i.e., the physical function) of the obtained particles to obtain its three-dimensional image; according to the three-dimensional image, count the concentration of inhalable particles with different particle sizes, and transmit the reproduced image and particle concentration to the display device , the transmission device adopts wired transmission or wireless transmission technology such as Bluetooth and Wifi. Combined with the volume of the air area irradiated by the light spot, the concentration of particles with different particle sizes can be calculated, and then converted into the corresponding mass concentration according to the national standard WS/T206-2001.
(5)显示设备将再现全息图像、不同粒径可吸入颗粒物浓度等信息显示在屏幕。(5) The display device will display information such as the reproduced holographic image and the concentration of inhalable particles with different particle sizes on the screen.
本发明的技术效果为:(1)本发明通过对可吸入颗粒物可以实现在线三维成像,分辨率高;采用波长为525nm的激光,细处理(即Err取值为10-5)时分辨率可达200nm;快速处理(即Err取值为10-3)分辨率为1μm。The technical effects of the present invention are as follows: (1) The present invention can realize online three-dimensional imaging of inhalable particulate matter with high resolution; using a laser with a wavelength of 525nm, the resolution can reach Up to 200nm; fast processing (that is, the value of Err is 10 -3 ) and a resolution of 1μm.
(2)本发明对可吸入颗粒物形貌再现过程的处理时间不超过5min,采用并行运算,可同时处理60~100个颗粒,此时颗粒数对处理时间无明显影响,因此本实验对可吸入颗粒物的准实时监测;结合光场对颗粒衍射图斑的定位,可以显示颗粒运动轨迹,这可为研究其行为提供帮助。(2) The processing time of the present invention for reappearing the appearance of inhalable particulate matter is no more than 5 minutes, and 60 to 100 particles can be processed at the same time by using parallel computing. Quasi-real-time monitoring of particles; combined with the positioning of the particle diffraction pattern by the light field, the particle trajectory can be displayed, which can help to study its behavior.
(3)本发明可以对不同粒径的可吸入颗粒物浓度同时进行统计计算,无需使用切割器。(3) The present invention can simultaneously perform statistical calculations on the concentration of inhalable particulate matter with different particle sizes without using a cutter.
附图说明Description of drawings
图1 一种基于层叠衍射成像技术的可吸入颗粒物检测方法所用装置示意图, 1为激光器,2为透镜组,22为滤光小孔,3为气室(含吸收和泵出装置的阀门),31为气压表,4为高速CCD,5为便携式电子计算机(含处理器和显示设备)。Figure 1. Schematic diagram of the device used in a detection method for inhalable particulate matter based on stacked diffraction imaging technology, 1 is a laser, 2 is a lens group, 22 is a small filter hole, 3 is an air chamber (including the valve of the absorption and pumping device), 31 is a barometer, 4 is a high-speed CCD, and 5 is a portable electronic computer (including a processor and a display device).
图2 人工模拟环境中可吸入颗粒物监测的衍射图样,四幅图按1-4次序分别为相邻四个时间点图样,图中标尺代表20μm。Figure 2 Diffraction patterns of inhalable particulate matter monitoring in artificially simulated environments. The four pictures are the pictures of four adjacent time points in the order of 1-4, and the scale bar in the figure represents 20 μm.
图3 人工模拟环境中可吸入颗粒再现图,子图为最左侧颗粒(子图正下方颗粒)再现图的放大图,主图标尺代表20μm,子图标尺代表1μm。Figure 3 Reproduction of inhalable particles in the artificial simulated environment, the sub-figure is an enlarged view of the reproduction of the leftmost particle (the particle directly below the sub-figure), the scale of the main graph represents 20 μm, and the scale of the sub-graph represents 1 μm.
具体实施方式Detailed ways
下面通过实例对本发明做进一步说明。需要注意的是,公布实施例的目的在于帮助进一步理解本发明,但是本领域的技术人员可以理解:在不脱离本发明及所附权利要求精神和范围内,各种替换和修改都是可能的。因此,本发明不应局限于实施例所公开的内容,本发明要求保护的范围以权利要求书界定的范围为准。The present invention will be further described below by example. It should be noted that the purpose of the published embodiments is to help further understand the present invention, but those skilled in the art can understand that various replacements and modifications are possible without departing from the spirit and scope of the present invention and the appended claims . Therefore, the present invention should not be limited to the content disclosed in the embodiments, and the protection scope of the present invention is subject to the scope defined in the claims.
实施例:基于层叠衍射算法对模拟环境中可吸入颗粒物的实时检测Example: Real-time detection of inhalable particulate matter in a simulated environment based on the stacked diffraction algorithm
搭建光路图并连接、调试相应设备Build the optical path diagram and connect and debug the corresponding equipment
依照附图1所示光路图搭建检测可吸入颗粒物的光路并连接相应的设备,所用气室在激光照射区域内部厚度为1mm,接收设备4为拍摄速率为1000fps的高速CCD。激光经过透镜组2成球面波,球心位置在光学小孔22处。According to the optical path diagram shown in Figure 1, the optical path for detecting inhalable particulate matter was built and connected to the corresponding equipment. The thickness of the gas chamber used in the laser irradiation area was 1mm, and the receiving device 4 was a high-speed CCD with a shooting rate of 1000fps. The laser light passes through the lens group 2 to form a spherical wave, and the center of the sphere is at the optical small hole 22 .
获取背景信息Get background information
泵出气室3内空气;待气压表示数稳定在-0.1Mpa时,表明气室内空气抽取完毕;将波长为645nm的激光通过扩束后沿气室厚度方向照向气室,气室前表面(即靠近激光器的表面)光斑直径约为1cm,通过CCD接收光斑,将此光斑存储作为背景信息。Pump out the air in the air chamber 3; when the air pressure indicator is stable at -0.1Mpa, it indicates that the air in the air chamber has been extracted; the laser with a wavelength of 645nm is irradiated to the air chamber along the thickness direction of the air chamber after beam expansion, and the front surface of the air chamber ( That is, the diameter of the spot near the laser is about 1cm, and the spot is received by the CCD, which is stored as background information.
环境平衡environmental balance
将气室3输入装置和输出装置打开,待气压表31示数不变时,表示气室内空气即为环境中气体,封闭输入和输出装置,并静置30s;环境采用人工模拟环境,即在密闭环境中点燃香烟产生颗粒。Open the input device and output device of the air chamber 3, and when the indication of the barometer 31 remains unchanged, it means that the air in the air chamber is the gas in the environment, close the input and output devices, and let it stand for 30 seconds; Lighting a cigarette in a closed environment produces particles.
获得所探测空气衍射图Obtaining the Diffraction Pattern of the Detected Air
将激光从获取背景信息的照射位置垂直射入气室,并通过高速CCD接收衍射光斑,每两个相邻光斑接收时间不超过1/240s,所得空气的典型衍射图如图2所示,其表示相邻4个时刻的衍射图样。The laser is vertically injected into the air chamber from the irradiation position where the background information is obtained, and the diffraction spot is received by a high-speed CCD. The receiving time of each two adjacent spots does not exceed 1/240s. The typical diffraction pattern of the obtained air is shown in Figure 2. Represents the diffraction patterns at four adjacent moments.
获得所探测空气的再现像Obtain a reconstructed image of the detected air
根据空气的衍射图进行层叠衍射计算。通过典型的衍射图2可以看出,4个时刻的灰尘衍射斑有重叠,因此可以通过层叠衍射算法进行计算。计算方法按照技术方案中步骤(3)中的方法进行,其第六步中的α、β取值为1。经过迭代所得到图2对应的可吸入颗粒物像如图3所示,图3的子图为最左侧颗粒的放大像,由此可以看出该方法能够对PM2.5等颗粒进行三维成像。Stacked diffraction calculations are performed from the diffraction pattern of air. It can be seen from the typical diffraction pattern 2 that the dust diffraction spots at four moments overlap, so it can be calculated by the stacked diffraction algorithm. The calculation method is carried out according to the method in step (3) of the technical solution, and the values of α and β in the sixth step are 1. The inhalable particle image corresponding to Figure 2 obtained after iteration is shown in Figure 3, and the sub-image of Figure 3 is the enlarged image of the leftmost particle, which shows that this method can perform three-dimensional imaging of PM2.5 and other particles.
获得可吸入颗粒物浓度Obtaining the concentration of respirable particulate matter
处理器自动对将样本再现像按其大小进行分类统计,扣除背景信息后,统计出其个数并计算其体积大小,然后和所照射空气体积比较,并换算为标准状况下空气体积,计算出体浓度;根据国标WS/T206-2001换算出相应的质量浓度和污染指数。The processor automatically classifies and counts the reproduced image of the sample according to its size. After deducting the background information, counts its number and calculates its volume, and then compares it with the irradiated air volume, and converts it to the air volume under standard conditions, and calculates body concentration; according to the national standard WS/T206-2001, the corresponding mass concentration and pollution index are converted.
信息显示the message states
将颗粒物再现像和颗粒物浓度、污染指数等信息显示在便携式计算机5的屏幕上;不断刷新屏幕,可以看到颗粒物动态运动过程。The reappearance image of the particles, the concentration of the particles, the pollution index and other information are displayed on the screen of the portable computer 5; the screen is constantly refreshed, and the dynamic movement process of the particles can be seen.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610012762.2A CN105547944B (en) | 2016-01-11 | 2016-01-11 | A kind of inhalable particles object detecting method based on stacking diffraction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610012762.2A CN105547944B (en) | 2016-01-11 | 2016-01-11 | A kind of inhalable particles object detecting method based on stacking diffraction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105547944A CN105547944A (en) | 2016-05-04 |
CN105547944B true CN105547944B (en) | 2018-04-10 |
Family
ID=55827286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610012762.2A Expired - Fee Related CN105547944B (en) | 2016-01-11 | 2016-01-11 | A kind of inhalable particles object detecting method based on stacking diffraction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105547944B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107328694A (en) * | 2017-08-03 | 2017-11-07 | 京东方科技集团股份有限公司 | A kind of particle detection and its detection method |
GB2569088B (en) * | 2017-08-16 | 2021-05-19 | Corey Kaizen Reaux Savonte | Matter-based photon quantum encryption device |
CN108254295B (en) * | 2018-01-15 | 2020-07-24 | 南京大学 | Method and device for positioning and representing spherical particles |
CN109269952A (en) * | 2018-10-16 | 2019-01-25 | 东南大学 | A kind of measuring device and method of filamentary material three dimensional particles concentration information in air flow bed rectangle riser |
CN110057729A (en) * | 2019-03-12 | 2019-07-26 | 天津大学 | A kind of human body particulate matter passively settles concentration determination method |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3873206A (en) * | 1973-10-03 | 1975-03-25 | Leeds & Northrup Co | Method for determining a specific characteristic of fluid suspended particles |
SU1689801A1 (en) * | 1989-07-19 | 1991-11-07 | Киевский Институт Инженеров Гражданской Авиации Им.60-Летия Ссср | Device for determination of aerosol particle size |
EP1591806A1 (en) * | 2004-04-28 | 2005-11-02 | Canon Kabushiki Kaisha | Diffractive optical element |
CN101842751A (en) * | 2007-10-30 | 2010-09-22 | 纽约大学 | Tracking and characterizing particles with holographic video microscopy |
CN103531427A (en) * | 2012-06-29 | 2014-01-22 | Fei公司 | Method for preparing and imaging a lamella in a particle-optical apparatus |
US20150036038A1 (en) * | 2013-07-31 | 2015-02-05 | California Institute Of Technology | Aperture scanning fourier ptychographic imaging |
CN104897538A (en) * | 2015-06-16 | 2015-09-09 | 天津大学 | Construction site raising dust detecting device based on digital holography |
CN105259668A (en) * | 2015-10-12 | 2016-01-20 | 中国科学院大学 | Black support based lamination imaging technology |
-
2016
- 2016-01-11 CN CN201610012762.2A patent/CN105547944B/en not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3873206A (en) * | 1973-10-03 | 1975-03-25 | Leeds & Northrup Co | Method for determining a specific characteristic of fluid suspended particles |
SU1689801A1 (en) * | 1989-07-19 | 1991-11-07 | Киевский Институт Инженеров Гражданской Авиации Им.60-Летия Ссср | Device for determination of aerosol particle size |
EP1591806A1 (en) * | 2004-04-28 | 2005-11-02 | Canon Kabushiki Kaisha | Diffractive optical element |
CN101842751A (en) * | 2007-10-30 | 2010-09-22 | 纽约大学 | Tracking and characterizing particles with holographic video microscopy |
CN103531427A (en) * | 2012-06-29 | 2014-01-22 | Fei公司 | Method for preparing and imaging a lamella in a particle-optical apparatus |
US20150036038A1 (en) * | 2013-07-31 | 2015-02-05 | California Institute Of Technology | Aperture scanning fourier ptychographic imaging |
CN104897538A (en) * | 2015-06-16 | 2015-09-09 | 天津大学 | Construction site raising dust detecting device based on digital holography |
CN105259668A (en) * | 2015-10-12 | 2016-01-20 | 中国科学院大学 | Black support based lamination imaging technology |
Non-Patent Citations (4)
Title |
---|
Direct extraction of the mean particle size from a digital hologram;Loïc Deni et.al;《APPLIED OPTICS》;20060210;第45卷(第5期);第944-952页 * |
一种高分辨率3 维图像的自适应降噪算法;向志聪 等;《激光技术》;20150930;第39卷(第5期);第697-701页 * |
同轴数字全息术在大气可吸入颗粒物PM2.5测量中的应用;王萍;《光学仪器》;20150831;第37卷(第4期);第287-292页 * |
实用化数字全息记录和再现系统的设计与实现;郭中甲;《中国优秀硕士学位论文全文数据库基础科学辑》;20130415(第4期);第1-57页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105547944A (en) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105547944B (en) | A kind of inhalable particles object detecting method based on stacking diffraction | |
Northcross et al. | A low-cost particle counter as a realtime fine-particle mass monitor | |
CN104316443B (en) | A PM 2.5 Concentration Monitoring Method Based on CCD Backscattering | |
CN105334144B (en) | A kind of monodisperse aerosol grain size and apparatus for measuring concentration based on light scattering | |
CN105606499B (en) | Suspended particulate matter mass concentration real-time detection device, and measuring method | |
CN105891063A (en) | Multi-angle dynamic light scattering particle size distribution measuring device and method | |
CN108519311B (en) | Real-time smoke particle density detection device and method | |
CN106124410A (en) | Single photoacoustic cell measures the new method of aerosol multi-wavelength absorptance simultaneously | |
CN106979909A (en) | Particulate matter classification and concentration detection sensor and method based on luminous flux ratio | |
CN111595801A (en) | Remote sensing identification and estimation method for brown carbon component of whole-layer atmospheric aerosol | |
CN112033913A (en) | Apparatus and method for measuring water content of nano-single-particles based on surface plasmon resonance imaging | |
CN112782121B (en) | A multi-angle optical particle counting and refractive index online measurement device and method | |
CN111947783B (en) | Polarization transmission characteristic test and simulation verification method for non-spherical multiple scattering medium | |
CN110441324B (en) | Detection device and detection method for aerosol particles based on Sagnac interferometer | |
CN113776997B (en) | Atmospheric particulate matter Mueller matrix measurement method with environment interference suppression function | |
CN104458512A (en) | Method for measuring fractal dimension of particle group | |
CN107894395A (en) | A kind of Aerosol Extinction Coefficients measuring method based on Research on Cavity Ring Down Spectroscopy | |
CN100344959C (en) | Device and method for checking quality of fibre-optical | |
Yin et al. | Airborne particulate matter distribution in urban green space is size-dependent | |
CN113405961A (en) | Dust concentration detection method based on multi-angle light scattering | |
Zhao et al. | Optimal static light scattering detection angle for particulate matter size and concentration measurement | |
Yao et al. | Measurements of the Optical Scattering Properties of Single Suspended Particles and Implications for Atmospheric Studies: A Review | |
Hickey | Methods of aerosol particle size characterization | |
CN208537351U (en) | A kind of smoke particle density real-time detection apparatus | |
CN106442248A (en) | Measurement method of emission concentrations of particulate matters in gases and system of measurement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180410 Termination date: 20200111 |
|
CF01 | Termination of patent right due to non-payment of annual fee |