Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

On-line high-accuracy particulate matter monitoring technology using multi-channel scattering signals

Open Access Open Access

Abstract

Particulate matter has adverse effects on the environment and human health, thus emission monitoring of particulate matter is essential for environment and human health protection. Optical methods are popular for on-line particulate matter emission monitoring due to their low cost, high sensitivity and easy maintainability. However, the measurement accuracy is susceptible to the particle size distribution of the particulate matter. To resolve this problem, a new optical method using multi-channel scattering signals and a proof-of-concept prototype sensor are proposed in this paper. According to multi-channel scattering signals, which reflect the change of particle size distribution, the prototype sensor adaptively sets the mass scattering coefficient to improve the mass concentration measurement accuracy. Compared with the state-of-the-art optical technologies, simulation results show that the relative standard deviation was reduced from 242% to 4% by our method. In the tests of our prototype sensor, the maximum relative measurement errors are 10% for di-ethylhexyl-sebacat (DEHS) monodisperse aerosols and 11% for coal smoke. Given that it is low cost, highly sensitive and easy to maintain, the new method has significant potential for precise measurement of particulate matter in mobile or stationary pollution monitoring applications.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Particulate matter (PM) has drawn extensive attention due to its adverse effects on environment and human health. Recent reports from Harvard University and Germany indicate that the loss of life as a result of atmospheric pollutants is approximately 8 million people per year. Thus, methods for reliably monitoring these PM are of critical importance [15]. Large-scale air pollution studies have revealed that the mass concentration of PM was positively related to the mortality [68]. It is essential to monitor the mass concentration of PM emitted from the stationary sources, such as thermal power plants, utility boilers, and fluidized beds for semi-dry desulfurization and so on. The manual gravimetric method is the standard method for aerosol mass concentration measurement. It needs several days for sampling and analysis procedure, which is labor-consuming and expensive [9]. Therefore, several technologies are developed for automatic particulate monitoring. The commercial instruments are mainly designed based on tapered element oscillating microbalances (TEOM) or Beta-attenuation technologies [10,11]. These automatic monitoring instruments are very expensive and need hours to collect the particles for each measurement [12]. Furthermore, these instruments require periodic replacement of the particle collection filters every 2-3 weeks, which significantly shortens the maintenance period [13]. Recently, several new methods are developed for particle analysis. Thermal optical methods provide nearly real time measurements of filter sample data [14]. Light absorption methods, photoacoustic methods, and laser-induced incandescence methods provide fast measurements of soot volume fraction with high resolution and accuracy down to low concentrations [15,16]. However, these methods are not suitable for commercial application due to high cost and poor maintainability.

Compared with the monitoring technologies mentioned above, light scattering methods provide highly-sensitive and rapid measurement of PM mass concentration. Meanwhile, as a non-contact measurement, the light scattering methods do not need particle collection filters, and can effectively extend the maintenance cycle to several months. However, the usage of light scattering methods is limited by the measurement accuracy. Existing optical light scattering methods adopt a linear model to calculate the PM mass concentration from the light signal, but the conversion coefficient is susceptible to the physical properties of the PM, especially the particle size distribution [1719]. Due to the variation of combustion condition and collection efficiency of the dust removing equipment, the particle size distribution of the PM eventually released into the atmosphere will be different. Our simulation results in Section 3.1 show that the expected measurement error of light scattering based technology reaches to 454% owing to the change of particle size distribution. As a compromise, Thermo Fisher corporation develops a hybrid particulate matter continuous emission monitoring system (PM-CEMS) Model 3880i combining light scattering technology for on-line monitoring and TEOM technology for periodical calibration [20]. Comparing with the instruments using single measurement technology, this kind of hybrid instrument is even more expensive and also requires frequent maintenance.

According to Mie theory, the light scattering signals contain the information of mass concentration and particle size distribution. Several methods have been developed to improve the mass concentration measurement accuracy only using optical signals. Meng Jiang et.al established a correction model of PM mass concentration measurement according to the mass media diameter, which is derived from the light extinction and scattering signals [21]. Li Wang et.al demonstrated a spectral light extinction method to measure the concentration and particle size distribution [22]. However, the light extinction method usually needs a long light path to ensure adequate sensitivity. Ring-down is an effective approach to reduce the size of the instruments and improves the measurement sensitivity by folding the optical path. To further simplify the structure and reduce the size of the sensor, the light scattering method is adopted. In our previous work, multi-wavelength technology with empirically selected wavelengths is proposed to measure the aerosol mass concentration and estimate the particle size distribution [23]. Dong Chen proposed a laboratory analysis method using a fixed detector and a rotating detector to obtain mass concentration and particle size respectively [24]. It implies that the light scattering signals at different observing angles also contain the information of mass concentration and particle size distribution.

In this work, a new method is proposed to improve the PM mass concentration measurement accuracy. Particle size distribution tracking technology is developed using multi-channel scattering signals with optimized wavelengths and observing angles. Instead of a preset constant value, the mass scattering coefficient is adaptively set to reduce the measurement errors caused by the change of particle size distribution. A proof-of-concept prototype sensor with 3 optimized light scattering channels is designed to verify the method proposed in this paper, which has a simple optical structure and only a few optical components. The experimental results of the prototype sensor show that our sensor is very sensitive and can precisely measure the mass concentration of DEHS aerosols and coal smoke.

The contribution of this study are as follows:

  • • We proposed a non-linear model to measure the aerosol mass concentration via optical signals, where the mass scattering coefficient is adaptively set with the change of particle size distribution.
  • • We developed a particle size distribution tracking technology using multi-channel scattering signals, where the optical parameters (number of measurement channel, wavelengths and observing angles) are comprehensively optimized to minimize the tracking error.
  • • We designed a proof-of-concept sensor, which is highly-sensitive, low-cost and easy to maintain. The simulation and experimental results show that our sensor can precisely measure the mass concentration of PM with different particle size distributions.

The following paper is organized into four sections: Section 2 describes the sensing method of PM mass concentration with adaptively adjusted mass scattering coefficient using multi-channel scattering signals. Section 3 shows experimental results of simulation, DEHS monodisperse aerosols and coal smoke. Section 4 concludes the whole paper.

2. Methods

2.1 Existing linear model for particulate matter mass concentration measurement

The PM mass concentration CM is proportional to the volume concentration CV, where the mass density ρ is a constant value obtained by laboratory calibration. Thus CM can be described by the particle size distribution f(x):

$${{C}_M} = \rho {\textrm{C}_\textrm{V}} = \rho {{C}_N}\int {f(x)(\pi {x^3}/6)dx}$$
where CN is the PM number concentration, x is the particle diameter. The mass concentration can be obtained by light scattering technologies. As illustrated by Van de Hulst et. al. [2527], the scattering light can be described by Mie scattering theory, where the particles are considered as sphere or optical equivalent sphere. The scattering light intensity P is also related to the particle size distribution:
$$P\; = {K}{{C}_N}\int {f\; (x)q(x,m,\lambda ,\theta )dx}$$
where q(x, m, λ, θ) describes the light intensity scattered by a single particle with diameter x based on Mie theory, m is the refractive index, λ is the wavelength of incident light, θ is the observing angle. The channel coefficient K is used to describe the systematic factors that would affect the measurement of scattering light intensity P, such as the incident light intensity, the volume of sampling cross section, the amplification coefficient of the amplifier, the response sensitivity of the photodiode, and so on.

To measure the PM mass concentration, existing PM-CEMS adopts a linear conversion mode with a preset empirical mass scattering coefficient T:

$${C_M}\; = \rho TP$$

As discussed above, Eq. (3) holds only if q(x, m, λ, θ) = πx3/6 T. However, as shown in Fig. 1(a), the scattering light intensity per volume concentration πx3/6q(x, m, λ, θ) is not constant with the change of the particle size. Thus, the simplified linear model Eq. (3) will cause mass concentration measurement errors. As an example shown in Fig. 1(b), when the scattering light intensities are consistent, the mass concentration of the aerosol sample B with bigger particle size reaches up to 3.56 times that of the aerosol sample A with smaller particle size. Since the mass scattering coefficient is a constant value in the linear model, existing PM-CEMS may fail to measure the mass concentration correctly when the PM size distribution changes. Specifically, the mass concentration of the small particle size aerosol will be overestimated, while that of the big particle size aerosol will be underestimated.

 figure: Fig. 1.

Fig. 1. The mass concentration measurement using the existing linear model will cause errors due to the variation of mass scattering coefficient. (a) The relative volume concentration per unit scattering intensity of a single particle with the change of the particle size. (b) The mass concentration of aerosols with different particle size distribution when the scattering light intensities are consistent.

Download Full Size | PDF

2.2 Non-linear model with adaptively adjusted mass scattering coefficient

As discussed above, the mass scattering coefficient T is related to the particle size distribution f(x) rather than a constant value, thus Eq. (3) is rewritten as:

$${C_M}\; = \rho T(f(x))P$$

As to PM emitted from stationary sources, the particle size distribution can be described by the lognormal distribution model [28,29]:

$${f_{[\mu ,\sigma ]}}(x) = \frac{1}{{\sqrt {2\pi } x\ln \sigma }}\textrm{exp} \left[ { - \frac{{{{({\ln x - \ln \mu } )}^2}}}{{2{{\ln }^2}\sigma }}} \right]$$
where μ denotes the count median diameter and σ describes the geometric standard deviation. Thus, the mass scattering coefficient T is associated with [μ, σ], namely T(μ, σ).

On the other hand, as illustrated by Eq. (2), the scattering light intensity P is also affected by the particle size distribution f[μ, σ](x). With only one measurement channel, existing PM-CEMS can only estimate the concentration via the linear model, while no information of particle size distribution can be obtained. To eliminate the effect of concentration, we use the ratio vector R(μ, σ) to track the change of particle size distribution, where the ith element Ri(μ, σ) is defined as the ratio of the scattering signals Pi with optical parameter set [λi, θi] to the basic scattering signal PB with [λB, θB]:

$${R_i}({\mu ,\sigma } )\;\textrm{ = }\frac{{{K_i}\int {{f_{[\mu ,\sigma ]}}(x )q({x,m,{\lambda_i},{\theta_i}} )dx} }}{{{K_B}\int {{f_{[\mu ,\sigma ]}}(x )q({x,m,{\lambda_B},{\theta_B}} )dx} }}\textrm{ }$$
Ki and KB are constant values once the sensor is fabricated, which can be calibrated by experiments. Since the refractive index m and optical parameters [λi, θi] are pre-known parameters, Ri(μ, σ) is only associated with [μ, σ], thus we established a conversion model G(·) from R(μ, σ) to T(μ, σ) by Mie theory (see example in Section 1 of Supplement 1).
$$T(\mu ,\sigma ) = G({\boldsymbol R}(\mu ,\sigma ))$$

Thus, the non-linear model from scattering signals to PM mass concentration CM is proposed as:

$${\textrm{C}_M} = \rho T(\mu ,\sigma ){P_B} = \rho G({\boldsymbol R}(\mu ,\sigma )){P_B}$$

As discussed above, the flow chart of mass concentration sensing method with multi-channel scattering signals is shown in Fig. 2. By measuring the scattering signals set P = {PB, P1, …, Pi}, our method calculates the ratio vector R = {R1, R2 …, Ri} to track the change of particle size distribution, and then adaptively set the mass scattering coefficient T by Eq. (7). Thus, the mass concentration of aerosol samples with different particle size distributions can be precisely measured according to Eq. (8).

 figure: Fig. 2.

Fig. 2. The flow chart of dynamically adjusted particle mass concentration sensing method using the non-linear model.

Download Full Size | PDF

2.3 Optical parameters optimization for particle size distribution tracking technology

The scattering ratio vector R(μ, σ) needs sufficient information to track the change of particle size distribution. Otherwise, measurement error will be introduced when R(μ, σ) of different particle size distributions are too similar to be differentiated. Specifically, the aerosol samples [μj, σj] and [μk, σk] are considered to be differentiable if

$$abs(1\textrm{ } - {R_i}({\mu _\textrm{j}},{\sigma _j})\textrm{ }/{R_i}({\mu _k},{\sigma _k}))\textrm{ } > {D_t},\forall \; {R_i} \in R$$
where the discrimination threshold Dt is determined by the systematic measurement errors of the scattering light, such as signal noise, stray light, truncation errors of ADC (analog-digital converter) and so on. As illustrated by Eq. (9), once the difference between any two corresponding elements Ri(μj, σj) and Ri(μk, σk) exceeds Dt, aerosol samples [μj, σj] and [μk, σk] can be discriminated by our sensor. Increasing measurement channels of the scattering light can provide more information to track the particle size distribution, but also increase the complexity of the optical system, leading to higher cost and lower stability. Therefore, optical parameter optimization of the scattering ratio vector [λ, $\boldsymbol{\theta}$] = {[λt, θt]|t=1,…,n} is essential to maximize the measurement accuracy with as few measurement channels n as possible. It should be noted that the scattering ratio vector R(μ, σ) can only be obtained when n≥2, while n=1 denotes the linear model of existing PM-CEMS with a constant mass scattering coefficient.

As described by Eq. (9), if the ith PM sample is distinguishable from all of the other PM samples, our non-linear model uses the practical mass scattering coefficient of ith PM sample TGTi as the mass scattering coefficient of the prototype sensor. Otherwise, the PM samples which cannot be distinguished from the ith PM sample make up an undistinguishable group [μS, σS] (including the ith PM sample). The mass scattering coefficients of all undistinguishable PM samples TS are the potential output of G(R(μi, σi)). As a compromise, Tavg is defined as the average values of all mass scattering coefficient of PM samples in [μS, σS] according to Eq. (10).

$$\begin{aligned} {T_{avg}}(\mu ,\sigma ) &= average({T_{GT}}({\mu _S},{\sigma _S}))\\ \textrm{ } &= average(G({\boldsymbol R}({\mu _S},{\sigma _S})))\textrm{ } \end{aligned}$$
where TGT(μS, σS) is the ground-truth mass scattering coefficient of the undistinguishable PM samples in [μS, σS]. For all k kinds of different particle size distributions, a relative standard deviation RSDm is defined to evaluate the mass concentration expected measurement error of the optical system:
$$\begin{aligned} RS{D_m}(\lambda ,{\kern 1pt} {\kern 1pt} \theta ) &= \sqrt {\frac{1}{k}\sum\limits_{S = 1}^k {{{\left( {\frac{{\rho {T_{avg}}({\mu_S},{\sigma_S})P - \rho {T_{GT}}({\mu_S},{\sigma_S})P}}{{\rho {T_{GT}}({\mu_S},{\sigma_S})P}}} \right)}^2}} } \\ &= \sqrt {\frac{1}{k}\sum\limits_{S = 1}^k {{{\left( {\frac{{{T_{avg}}({\mu_S},{\sigma_S}) - {T_{GT}}({\mu_S},{\sigma_S})}}{{{T_{GT}}({\mu_S},{\sigma_S})}}} \right)}^2}} } \end{aligned}$$

In order to minimize the measurement error, the optical parameter set [λ, $\boldsymbol{\theta}$] of the optical system are optimized according to Eq. (12).

$$[\lambda ,{\kern 1pt} {\kern 1pt} \theta ] = {\kern 1pt} {\kern 1pt} argmin{\kern 1pt} {\kern 1pt} RS{D_m}(\lambda ,{\kern 1pt} {\kern 1pt} \theta )$$

3. Typographical results and discussion

3.1 Optimum optical parameters of the particulate matter mass concentration monitoring system

The optical parameter set [λ, θ] is optimized by Eq. (12) according to the properties of PM. As a typical aerosol pollutant, the coal smoke emitted from the thermoelectricity power plants are investigated. With the help of dust removal equipment, particles eventually emitted from the stationary air pollution sources usually are smaller than 10 μm [3032], thus the particle size range is set to 100nm to 10000nm with an increment of 10nm in simulation. Previous studies [21,3335] indicate that the typical range of μ and σ is 200nm - 2500nm and 1.5 - 2.0, respectively. Since the wavelengths of 450nm, 650nm and 950nm are very close, the corresponding refractive indices are considered as the same. The particles are supposed to be uniformly mixed or consist of single substance, thus the refractive index of coal smoke is set to 1.7 + 0.1i. The parameters of the coal smoke are summarized in Table 1.

Tables Icon

Table 1. Parameters of the particulate matter emitted from the stationary sources

On the other hand, the design of PM monitoring system is limited by the availability of the optical components and mechanical structure. The wavelengths of the incident light are selected from 450 nm, 650 nm, and 950 nm in consideration of cost. Due to the limitation of mechanical design and the interference caused by the stray light, the observing angle is restricted from 40° to 140° with an increment of 5°. Meanwhile, the discrimination threshold Dt is selected to be 5% in the simulation, which is fair enough to cover the systematic measurement errors of the scattering light.

According to the ranges of optical parameters discussed above and PM samples listed in Table 1, [λ, $\boldsymbol{\theta}$] is optimized by Eq. (12) with traversal algorithm. As the optimum simulation results shown in Table 2, RSDm is reduced with increasing number of measurement channels.

Tables Icon

Table 2. Expected accuracy of the PM mass concentration by multi-channel scattering signals.

For further illustrating the mechanism of the non-linear conversion model, the adaptive mass scattering coefficient and relative error of each aerosol sample are shown in Fig. 3, where TGT and Tavg are normalized by the maximum value of TGT and the relative error REM is defined as:

$$\begin{aligned} \textrm{R}{\textrm{E}_M}(\mu ,\sigma ) &= ({C_{M\_avg}} - {C_{M\_GT}})/{C_{M\_GT}}\\ &= (\rho {T_{avg}}(\mu ,\sigma )P - \rho {T_{GT}}(\mu ,\sigma )P)/\rho {T_{GT}}(\mu ,\sigma )P\\ &= ({T_{avg}}(\mu ,\sigma ) - {T_{GT}}(\mu ,\sigma ))/{T_{GT}}(\mu ,\sigma ) \end{aligned}$$

 figure: Fig. 3.

Fig. 3. The reconstructed conversion factor T and expected measurement error of mass concentration with different measurement channels in simulation experiment. (a) Normalized TGT verses Tavg with 1 measurement channel. (b) REM with 1 measurement channel. (c) Normalized TGT verses Tavg with 2 measurement channels. (d) REM with 2 measurement channels. (e) Normalized TGT verses Tavg with 3 measurement channels. (f) REM with 3 measurement channels.

Download Full Size | PDF

As illustrated by Eq (13), REM is defined as the relative error of the mass concentration and can be simplified as the function of TGT and Tavg.

As the linear model in existing PM-CEMS (n=1) illustrated by Fig. 3(a) the particle size distribution cannot be tracked. Therefore, the mass scattering coefficient of each aerosol sample is the average mass scattering coefficient of all aerosol samples according to Eq. (10), where the relative measurement error ranges from −85% to 454%. With an additional measurement channel (n = 2) in Fig. 3(b), the scattering ratio vector is calculated to estimate the particle size distribution and adaptively set the mass scattering coefficient according to Eq. (7). The expected measurement error is significantly reduced by non-linear model with relative measurement error from −37% to 71%. It is worth noting that the measurement accuracy changes with the standard deviation σ, since 2 channels cannot provide sufficient information to estimate μ and σ simultaneously. Therefore, 3 measurement channels in Fig. 3(c) further reduce the measurement error, where the relative error ranges from −12% to 19%. As illustrated in Table 2, more scattering light channels (n>3) further reduce the measurement error slightly, but also increase the system complexity. In consideration of measurement accuracy and structural complexity, the mass concentration monitoring device with 3 measurement channels reaches a proper balance between performances and costs, where the optimum parameter sets are [950 nm - 40°, 950 nm 140°, 450 nm - 135°].

According to the simulation optimization results, we design the proof-of-concept prototype sensor as shown in Fig. 4. Compared with a light extinction method, the optical path of light scattering method is very short. Hence, the structure of the prototype sensor is simple and small with a height 70 mm and a diameter 80 mm. A 450 nm laser and a 950 nm laser are adopted as the light sources, where the polarization angle of each laser is 45° and the power is 5 mW, (More details of polarization angle are discussed in Section 2 of Supplement 1). The collimated incident lights travel through the aerosol sample and then are absorbed by the light traps with absorptive filters. The particles are considered to be evenly distributed in the sampling cross section. Thus, the intensities of the scattering lights are affected by the total light intensity and unaffected by the light intensity distributions of the incident lights. Two photodiodes are used to receive the scattering light signals at the specific observing angles. With limited light collecting aperture and small half-power angle of the PD, the view scopes of scattering light receivers are considered to be collimated as well (More design details of the prototype sensor can be found in Section 3 and 4 of Supplement 1). The integral amplify circuits and 16-bit Analog-Digital converters are adopted to provide highly sensitive measurement. In order to keep the sample flow rate constant, isokinetic sampling structures are adopted in the design of inlet pipe and outlet pipe. Our sensor provides real-time, on-line and non-contact measurement of PM and extends the maintenance period without using particle collection filter.

 figure: Fig. 4.

Fig. 4. The overall structure of the proof-of-concept prototype sensor. (a) and (b) are the 3D mechanical design of the prototype sensor. (c) and (d) are the pictures of the prototype sensor. 1. 450 nm laser, 2. PD1, 3. 950 nm laser, 4. 450 nm light trap, 5. 950 nm light trap, 6. PD2, 7. Flue, 8. Inlet, 9. Outlet.

Download Full Size | PDF

3.2 Simulation experiments of particulate matter with random particle size distribution

A simulation experiment using 1000 random-generated particle size distributions are carried out to evaluate the performance of our prototype sensor. The ranges of the simulation parameters are consistent with those listed in Table 1, while the increment strides are not limited. Specifically, the count median diameter μ is random generated from 200nm to 2500nm and the geometric standard deviation σ is random generated from 1.5 to 2.0. The scattering ratio vectors of these simulation sample are calculated according to Mie theory, and thus the mass scattering coefficients are obtained.

As the simulation results shown in Fig. 5, the mass scattering coefficient is adaptively set with the change of particle size distribution. The relative deviation of the mass concentration ranges from –14% to 21% and corresponding relative standard error is 4%. It implies that the parameter intervals in Table 1 are appropriate and the non-linear model is suitable for any aerosol sample within the parameter range.

 figure: Fig. 5.

Fig. 5. The reconstructed conversion factor T and expected mass concentration measurement error of our prototype sensor to 1000 simulation aerosol pollutants with random particle size distributions. (a) Normalized TGT verses Tavg of our prototype sensor. (b) REM of our prototype sensor.

Download Full Size | PDF

3.3 Laboratory experiments of monodisperse DEHS aerosol samples and coal smoke

A laboratory validation measurements were conducted to measure the monodisperse DEHS aerosol samples with different particle sizes. As the experimental platform shown in Fig. 6, the aerosols samples are generated by the monodisperse aerosol generator (MAG, TSI MAG-3475), and then pumped into the dilution chamber for testing. Meanwhile, as the reference instruments, scanning mobility particle sizer (SMPS, TSI SMPS-3938) and the aerodynamic particle sizer (APS, TSI APS-3321) are used to measure the particle size distributions. With the knowledge of particle size given by the reference instruments and the refractive index of DEHS 1.45, the channel coefficients K are calibrated according to the measurement results and simulation results.

 figure: Fig. 6.

Fig. 6. The aerosol testing and particle sizing platform. (a)The block diagram of aerosols testing platform. (b) The picture of aerosols testing platform.

Download Full Size | PDF

As the experimental results shown in Fig. 7(a) and (b), the output of existing linear model is 328% higher than the mass concentration given by the SMPS. Meanwhile, our prototype sensor can significantly improve the mass concentration measurement accuracy by adaptively set the mass scattering coefficient with the change of particle size distribution, where the results are in agreement to the SMPS within 10%.

 figure: Fig. 7.

Fig. 7. Experimental results of single channel system, our prototype sensor and the reference instruments. The mass concentrations (a) and measurement deviation (b) of DEHS aerosol samples. The measurement results of mass concentrations (c) and measurement deviation (d) of coal smoke.

Download Full Size | PDF

As one of the major air pollutants, coal smoke are also tested by our prototype sensor. The experiments are repeated multiple times with different fuel quantities and aging times to obtain the coal smoke with different particle size distributions. As shown in Fig. 7(c) and (d), the deviation of the linear model to coal smoke ranges from 81% to 131%, while that of our prototype sensor is constrained within −9% to 11%. Correspondingly, the sensitivity of our sensor varies from 0.003 mg/m3 to 0.082 mg/m3 with the change of the mass scattering coefficient. Thus, our sensor can provide highly sensitive and precise measurement of coal smoke mass concentration.

4. Conclusion

In this paper, we propose a new optical method to precisely measure the mass concentration of PM. Compared with the state of art PM-CEMS using the linear model, our method adopts a non-linear model to adaptively set the mass scattering coefficient with the change of particle size distribution using the multi-channel scattering technology. A proof-of-concept prototype sensor with optimized parameter sets [950nm - 40°, 950nm - 140°, 450nm - 135°] was designed, which is simple-structure, low-cost and easy to maintain. According to the tests of the monodisperse DEHS aerosol samples and the coal smoke, our prototype sensor has significantly improved the mass concentration measurement accuracy in comparison with existing PM-CEMS using linear model. This new optical method can provide highly-sensitive, on-line and precise measurement of the PM mass concentration for both stationary and mobile aerosol pollution monitoring applications.

Funding

National Natural Science Foundation of China (61873322, 6210011211); Science Technology Foundation of Ministry of Emergency Management of China (2020XFZD13); Fundamental Research Funds for the Central Universities of Huazhong University of Science and Technology (2020kfyXJJS102).

Acknowledgments

We thank Dr. Zheng Dou for the help in the design of mechanical structure and electronic system.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

References

1. K. Vohra, A. Vodonos, J. Schwartz, E. A. Marias, M. P. Sulprizio, and L. J. Mickley, “Global mortality from outdoor fine particle pollution generated by fossil fuel combustion: Results from GEOS-Chem,” Environmental Research. 195, 110754 (2021). [CrossRef]  

2. W. W. Delp and B. C. Singer, “Wildfire Smoke Adjustment Factors for Low Cost and Professional PM2.5 Monitors with Optical Sensors,” Sensors. 20(13), 3683 (2020). [CrossRef]  

3. X. Y. Li, S. Q. Huang, A. Q. Jiao, X. H. Yang, J. F. Yun, Y. X. Wang, X. W. Xue, Y. Y. Chu, F. F. Liu, Y. S. Liu, M. Ren, X. Chen, N. Li, Y. A. Lu, Z. F. Mao, L. Q. Tian, and H. Xiang, “Association between Ambient Fine Particulate Matter and Preterm Birth or Term Low Birth Weight: An Updated Systematic Review and Meta-Analysis,” Environmental Pollution. 227, 596–605 (2017). [CrossRef]  

4. R. Lu, Y. Li, W. Li, Z. Xie, C. Fan, P. Liu, and S. Deng, “Bacterial Community Structure in Atmospheric Particulate Matters of Different Sizes during the Haze Days in Xi’an, China,” Science of the Total Environment. 637-638, 244–252 (2018). [CrossRef]  

5. W. C. V. Wang, T. H. Lin, C. H. Liu, C. W. Su, and S. C. C. Lung, “Fusion of Environmental Sensing on Pm2.5 and Deep Learning on Vehicle Detecting for Acquiring Roadside Pm2.5 Concentration Increments,” Sensors. 20(17), 4679 (2020). [CrossRef]  

6. Y. Guo, J. Teixeira, and N. Ryti, “Ambient Particulate Air Pollution and Daily Mortality in 652 Cities,” New England Journal of Medicine. 381(8), 705–715 (2019). [CrossRef]  

7. H. Y. Qi, C. Wang, W. X. Dong, and L. Meng, “Analyzing on Monitoring Technology Situation of Fine Particulate Matters Emitted from Stationary Source,” Environmental Science and Technology. 29, 69–74 (2016).

8. B. B. Wu, X. X. Bail, W. Liu, C. Y. Zhu, and H. Z. Tian, “Variation Characteristics of Final Size-Segregated PM Emissions from Ultralow Emission Coal-Fired Power Plants in China,” Environmental Pollution. 259, 113886 (2020). [CrossRef]  

9. Q. Zhifu, C. Shou, Z. Zou, W. D. Li, B. Chen, and C. Liu, “Research on test method of low concentration particulate matter in flue gas with low temperature and high humidity,” Environmental Pollution & Control. 07, 94–98 (2017).

10. L. Fu, T. Nunifu, and B. Leung, “A Two-Step Approach for Relating Tapered Element Oscillating Microbalance and Dichotomous Air Sampler PM2.5 Measurements,” Journal of the Air and Waste Management Association. 64(10), 1195–1203 (2014). [CrossRef]  

11. K. Takahashi, H. Minoura, and K. Sakamoto, “Examination of Discrepancies between Beta-Attenuation and Gravimetric Methods for the Monitoring of Particulate Matter,” Atmospheric Environment. 42(21), 5232–5240 (2008). [CrossRef]  

12. K. Gong, T. Shi, P. A. Ramachandran, K. W. Hutchenson, and B. Subramaniam, “Adsorption/Desorption Studies of 224-Trimethylpentane in β-Zeolite and Mesoporous Materials Using a Tapered Element Oscillating Microbalance (TEOM),” Ind. Eng. Chem. Res. 48(21), 9490–9497 (2009). [CrossRef]  

13. A. Winkel, J. L. Rubio, I. Huis, J. Vonk, and N. Ogink, “Equivalence testing of filter-based, beta-attenuation, TEOM, and light-scattering devices for measurement of PM10 concentration in animal houses,” Journal of Aerosol Science. 80, 11–26 (2015). [CrossRef]  

14. T. Ahmed, V. A. Dutkiewicz, A. Shareef, G. Tuncel, S. Tuncel, and L. Husain, “Measurement of black carbon (BC) by an optical method and a thermal-optical method: Intercomparison for four sites,” Atmos. Environ. 43(40), 6305–6311 (2009). [CrossRef]  

15. J. W. Cremer, P. A. Covert, E. A. Parmentier, and R. Signorell, “Direct measurement of photoacoustic signal sensitivity to aerosol particle size,” J. Phys. Chem. Lett. 8(14), 3398–3403 (2017). [CrossRef]  

16. AVL List GmbH, “Avl Micro Soot Sensor Operating Manual, Rev. 14,” Graz, Austria, 2018.

17. Z. Sui, Y. Zhang, Y. Peng, P. Norris, Y. Cao, and W. Pan, “Fine Particulate Matter Emission and Size Distribution Characteristics in an Ultra-Low Emission Power Plant,” Fuel 185, 863–871 (2016). [CrossRef]  

18. W. F. Tu, J. D. Lin, and Y. L. Wu, “An Analysis of Extinction Coefficients of Particles and Water Moisture in the Stack after Flue Gas Desulfurization at a Coal-Fired Power Plant,” Journal of the Air and Waste Management Association. 61(8), 815–825 (2011). [CrossRef]  

19. K. Xu, H. Zhang, W. F. Shen, Y. Zhang, and J. F. Lyu, “Soot Formation and Distribution in Coal Jet Flames over a Broad Range of Coal Concentration,” Energy Fuels 34(6), 7545–7553 (2020). [CrossRef]  

20. Brown and John, “Using an Optical PM CEMS with Wet FGD for MATS Compliance,” Power 159, 36–42 (2015).

21. M. Jiang, X. W. Liu, J. K. Han, Z. F. Wang, and M. H. Xu, “Combined Light Extinction and Scattering Measurement for Measuring a Low-Particulate-Mass Concentration with a White Cell-Based Optical System,” Energy Fuels 34(6), 7726–7734 (2020). [CrossRef]  

22. L. Wang and X. Jian, “Determination of particle size distribution by light extinction method using improved pattern search algorithm with Tikhonov smoothing functional,” Journal of Modern Optics. 59(21), 1829–1840 (2012). [CrossRef]  

23. X. Xiao, S. Wang, M. Zhu, T. Deng, A. Chen, and J. Zeng, “Novel three-wavelength optical sensor for measuring distributed mass concentration of aerosols from stationary sources,” Opt. Express 29(5), 6407 (2021). [CrossRef]  

24. D. Chen, X. W. Liu, J. K. Han, M. Jiang, Z. F. Wang, and J. X. Qi, “A New Angular Light Scattering Measurement of Particulate Matter Mass Concentration for Homogeneous Spherical Particles,” Sensors 19(10), 2243 (2019). [CrossRef]  

25. H. C. V. D. Hulst and G. W. Kattawar, “Exact Spread Function for a Pulsed Collimated Beam in a Medium with Small-Angle Scattering,” Appl. Opt. 33(24), 5820–5829 (1994). [CrossRef]  

26. R. T. Wang and H. C. V. D. Hulst, “Application of the Exact Solution for Scattering by an Infinite Cylinder to the Estimation of Scattering by a Finite Cylinder,” Appl. Opt. 34(15), 2811 (1995). [CrossRef]  

27. H. C. V. D. Hulst and R. T. Wang, “Glare points,” Appl. Opt. 30(33), 4755–4763 (1991). [CrossRef]  

28. T. Deng, J. Zeng, S. Wang, S. Yan, and A. Chen, “An Optical Fire Detector with Enhanced Response Sensitivities for Black Smoke Based on the Polarized Light Scattering,” Measurement Science and Technology. 30(11), 115203 (2019). [CrossRef]  

29. S. Wang, X. Xiao, T. Deng, A. Chen, and M. Zhu, “A Sauter Mean Diameter Sensor for Fire Smoke Detection,” Sensors and Actuators, B: Chemical. 281, 920–932 (2019). [CrossRef]  

30. Y. F. Li, W. D. Zhao, S. H. Xu, and W. C. Xia, “Changes of Size, Ash and Density of Coal Particles on the Column Axis of a Liquid-Solid Fluidized Bed,” Powder Technology. 245, 251–254 (2013). [CrossRef]  

31. E. Mészáros, T. Barcza, A. Gelencser, J. Hlavay, and K. Polyak, “Size Distributions of Inorganic and Organic Species in the Atmospheric Aerosol in Hungary,” Journal of Aerosol Science. 28(7), 1163–1175 (1997). [CrossRef]  

32. S. K. Sahu, M. Tiwari, R. C. Bhangare, and G. G. Pandit, “Enrichment and Particle Size Dependence of Polonium and Other Naturally Occurring Radionuclides in Coal Ash,” Journal of Environmental Radioactivity. 138, 421–426 (2014). [CrossRef]  

33. D. Chen, X. W. Liu, J. K. Han, M. Jiang, Y. S. Xu, and M. H. Xu, “Measurements of Particulate Matter Concentration by the Light Scattering Method: Optimization of the Detection Angle,” Fuel Processing Technology. 179, 124–134 (2018). [CrossRef]  

34. R. L. Ming, Y. Q. Zheng, L. L. Hua, and T. H. Ping, “The study on complex refractive indices of ash particles and bulk,” Journal of engineering thermophysics. 301(5), 13–19 (1997).

35. Y. S. Xu, X. W. Liu, H. Wang, Y. F. Zhang, J. X. Qi, and M. H. Xu, “Investigation of Simultaneously Reducing the Emission of Ultrafine Particulate Matter and Heavy Metals by Adding Modified Attapulgite during Coal Combustion,” Energy Fuels 33(2), 1518–1526 (2019). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Supplement 1

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. The mass concentration measurement using the existing linear model will cause errors due to the variation of mass scattering coefficient. (a) The relative volume concentration per unit scattering intensity of a single particle with the change of the particle size. (b) The mass concentration of aerosols with different particle size distribution when the scattering light intensities are consistent.
Fig. 2.
Fig. 2. The flow chart of dynamically adjusted particle mass concentration sensing method using the non-linear model.
Fig. 3.
Fig. 3. The reconstructed conversion factor T and expected measurement error of mass concentration with different measurement channels in simulation experiment. (a) Normalized TGT verses Tavg with 1 measurement channel. (b) REM with 1 measurement channel. (c) Normalized TGT verses Tavg with 2 measurement channels. (d) REM with 2 measurement channels. (e) Normalized TGT verses Tavg with 3 measurement channels. (f) REM with 3 measurement channels.
Fig. 4.
Fig. 4. The overall structure of the proof-of-concept prototype sensor. (a) and (b) are the 3D mechanical design of the prototype sensor. (c) and (d) are the pictures of the prototype sensor. 1. 450 nm laser, 2. PD1, 3. 950 nm laser, 4. 450 nm light trap, 5. 950 nm light trap, 6. PD2, 7. Flue, 8. Inlet, 9. Outlet.
Fig. 5.
Fig. 5. The reconstructed conversion factor T and expected mass concentration measurement error of our prototype sensor to 1000 simulation aerosol pollutants with random particle size distributions. (a) Normalized TGT verses Tavg of our prototype sensor. (b) REM of our prototype sensor.
Fig. 6.
Fig. 6. The aerosol testing and particle sizing platform. (a)The block diagram of aerosols testing platform. (b) The picture of aerosols testing platform.
Fig. 7.
Fig. 7. Experimental results of single channel system, our prototype sensor and the reference instruments. The mass concentrations (a) and measurement deviation (b) of DEHS aerosol samples. The measurement results of mass concentrations (c) and measurement deviation (d) of coal smoke.

Tables (2)

Tables Icon

Table 1. Parameters of the particulate matter emitted from the stationary sources

Tables Icon

Table 2. Expected accuracy of the PM mass concentration by multi-channel scattering signals.

Equations (13)

Equations on this page are rendered with MathJax. Learn more.

C M = ρ C V = ρ C N f ( x ) ( π x 3 / 6 ) d x
P = K C N f ( x ) q ( x , m , λ , θ ) d x
C M = ρ T P
C M = ρ T ( f ( x ) ) P
f [ μ , σ ] ( x ) = 1 2 π x ln σ exp [ ( ln x ln μ ) 2 2 ln 2 σ ]
R i ( μ , σ )  =  K i f [ μ , σ ] ( x ) q ( x , m , λ i , θ i ) d x K B f [ μ , σ ] ( x ) q ( x , m , λ B , θ B ) d x  
T ( μ , σ ) = G ( R ( μ , σ ) )
C M = ρ T ( μ , σ ) P B = ρ G ( R ( μ , σ ) ) P B
a b s ( 1   R i ( μ j , σ j )   / R i ( μ k , σ k ) )   > D t , R i R
T a v g ( μ , σ ) = a v e r a g e ( T G T ( μ S , σ S ) )   = a v e r a g e ( G ( R ( μ S , σ S ) ) )  
R S D m ( λ , θ ) = 1 k S = 1 k ( ρ T a v g ( μ S , σ S ) P ρ T G T ( μ S , σ S ) P ρ T G T ( μ S , σ S ) P ) 2 = 1 k S = 1 k ( T a v g ( μ S , σ S ) T G T ( μ S , σ S ) T G T ( μ S , σ S ) ) 2
[ λ , θ ] = a r g m i n R S D m ( λ , θ )
R E M ( μ , σ ) = ( C M _ a v g C M _ G T ) / C M _ G T = ( ρ T a v g ( μ , σ ) P ρ T G T ( μ , σ ) P ) / ρ T G T ( μ , σ ) P = ( T a v g ( μ , σ ) T G T ( μ , σ ) ) / T G T ( μ , σ )
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.