基于图像分析的掘进工作面粉尘颗粒检测方法

龚晓燕 冯浩 付浩然 陈龙 常虎强 刘壮壮 贺子纶 裴晓泽 薛河

文章編号：1671?251X（2024）04?0055?09

DOI：10.13272/j.issn.1671-251x.2023100074

摘要：基于光散射原理测定粉尘质量浓度只能定时定点手动检测，实时性差，且只能检测出粉尘质量浓度，并不能给出粒径分布范围。目前基于图像分析的粉尘颗粒检测研究主要是针对粉尘质量浓度或粒径分布进行单方面研究，并不能实现粉尘质量浓度和粒径分布范围的同时检测。针对上述问题，提出了一种基于图像分析的掘进工作面粉尘颗粒检测方法，探究图像特征与粉尘质量浓度、粒径分布间的关系。通过粉尘样本收集及图像采集装置，采集粉尘颗粒图像并获取采集图像时的粉尘质量浓度。编写粉尘样本图像处理算法，提取图像的灰度特征、纹理特征、几何特征相关参数。对提取的图像特征与实测粉尘质量浓度进行相关性分析，选取相关性较大的图像特征作为参数建立回归数学模型。提取粉尘颗粒对象像素点个数，结合转换系数，基于几何当量等效面积径计算粉尘粒径大小及分布范围。实验结果表明：实测粉尘质量浓度与建立的图像特征多元非线性回归模型数学模型计算值间的平均相对误差为12.37%，标准实测粒径与几何当量等效面积径得到的粒径分布间的最大相对误差为8.63%，平均相对误差为6.37%，验证了基于图像特征的粉尘质量浓度回归数学模型和基于几何当量等效面积径分布数学模型的准确性。

关键词：掘进工作面；粉尘质量浓度；粉尘粒径分布范围；图像分析；几何当量等效面积径；皮尔逊相关系数

中图分类号：TD714  文献标志码：A

A method for detecting dust particles in excavation working face based on image analysis

GONG Xiaoyan1， FENG Hao1， FU Haoran1， CHEN Long1， CHANG Huqiang2， LIU Zhuangzhuang1，HE Zilun1， PEI Xiaoze1， XUE He1

（1. College of Mechanical Engineering， Xi'an University of Science and Technology， Xi'an 710054， China；

2. Shaanxi Coal Group ShenmuNingtiaota Mining Co.， Ltd.， Shenmu 719300， China）

Abstract： Based on the principle of light scattering， measuring dust concentration can only be done manually at fixed times and locations， with poor real-time performance. It can only detect dust mass concentration and cannot provide a range of particle size distribution. At present， research on dust particle detection based on image analysis mainly focuses on unilateral research on dust mass concentration or particle size distribution. It cannot achieve simultaneous detection of dust mass concentration and particle size distribution range. In order to solve the above problems， a method for detecting dust particles in excavation working face based on image analysis is proposed. It explores the relationship between image features and dust mass concentration and particle size distribution. By using a dust sample collection and image acquisition device， dust particle images are collectedand the dust mass concentration at the time of image acquisition is obtained. An image processing algorithm for dust samples， is developed to extract parameters related to grayscale features， texture features， and geometric features of the image. The correlation analysis between the extracted image features and the measured dust mass concentration is performed， and the image features with high correlation is selected as parameters to establish a regression mathematical model. The method extracts the number of pixels in the dust particle object. Combining with the conversion coefficient， the method calculates the size and distribution range of the dust particle based on the geometric equivalent area diameter. The experimental results show that the average relative error between the measured dust mass concentration and the calculated values of the established image feature multiple nonlinear regression model mathematical model is 12.37%. The maximum relative error between the standard measured particle size and the geometric equivalent area size obtained from the particle size distribution is 8.63%， and the average relative error is 6.37%. This verifies the accuracy of the image feature based dust mass concentration regression mathematical model and the geometric equivalent area diameter distribution mathematical model.

Key words： excavation working face; dust mass concentration; dust particle size distribution range; image analysis; geometric equivalent equivalent area diameter; Pearson correlation coefficient

0引言

随着近年来煤矿井下机械化水平的不断提高，掘进工作面工作效率提高的同时产尘量也在不断增大，井下粉尘问题越来越严峻[1-4]。粉尘质量浓度过高不仅存在安全隐患，也会对工作人员身体健康造成伤害。因此粉尘颗粒的检测对于实际生产有着重要作用。近年来，国内外学者对粉尘颗粒的测定展开了大量的研究[5-7]。根据粉尘检测的原理不同，主要可分为取样法（滤膜称重法[8]、微量振荡天平法[9]和β射线法[10]）和非取样法（光散射法[11]、电荷感应法[12]和超声波衰减法[13]）。以上检测方法中，光散射法因具有精度高、可靠性高的优点而得到广泛应用，但其主要为手动定点操作，程序繁琐、实时性差，不能连续检测，且只能检测出粉尘质量浓度，并不能给出粒径分布范围。

随着计算机与机器视觉技术的进步，基于图像分析的粉尘颗粒测定技术得到快速发展。陈锋等[14]探索了粉尘图像的平均灰度值与实测粉尘质量浓度之间的关系。刘丹丹等[15]提出了基于线激光辅助 CCD 图像的面维度粉尘质量浓度分布检测方法，通过颗粒群光散射理论和 Beer?Lambert 衰减定律求解粉尘质量浓度分布。叶方平等[16]研究了粉尘质量浓度与图像透光率之间的映射关系，得到了基于图像透光率特征值计算的粉尘质量浓度测量算法。刘伟华等[17]对获取到的粉尘图像进行分析，得出粉尘相关几何特征，克服了传统煤尘定时定点测量的缺陷。A. Konoshonkin 等[18]分别研究了2个不同波长的全光散矩阵，得到了光散消光系数与粉尘粒子空间取向数的关系，表明后向散射光矩陣的元素服从幂函数。 D. Baaga 等[19]通过测量粉尘平均质量浓度、瞬

浓度，研究了 PM10、PM4和 PM2.5粉尘颗粒的分布特征，建立了基于幂函数的功能模型，并预测同等条件下不同类型的粉尘质量浓度分布。M. C. Baddock 等[20]通过对能见度和相关矿物浓度的详细案例进行研究，指出能见度与粉尘质量浓度之间存在数学关系。目前对于通过图像分析检测粉尘颗粒的研究还处于探索阶段，且相关的研究主要是基于粉尘质量浓度或粒径分布进行单方面研究，并不能实现粉尘质量浓度和粒径分布范围的同时检测。

针对上述问题，本文提出了一种基于图像分析的掘进工作面粉尘颗粒检测方法。通过设计粉尘样本收集装置和图像采集装置，获取粉尘样本图像并编写图像处理算法。提取粉尘样本图像特征进行分析，探究图像特征与粉尘质量浓度、粒径分布间的关系。建立相关图像特征与粉尘质量浓度及粉尘颗粒粒径分布的数学模型。

1粉尘检测系统

1.1系统结构原理

粉尘检测系统主要由样本收集、图像采集和图像处理3个部分构成。通过粉尘样本收集装置和图像采集装置获取粉尘样本图像。通过编写图像处理算法相关程序，对传输到数据分析模块中的样本图像进行处理，提取粉尘样本图像的相关特征，建立图像特征与粉尘质量浓度及粉尘颗粒粒径分布的数学模

1.2粉尘样本收集装置

粉尘样本收集装置由电磁开关、抽气泵、排气管、进气管等构成，如图1所示。通过控制抽气泵开启，将环境中的粉尘通过进气管传送至密闭结构空间中的样本收集平台（样本收集平台是由载玻片和盖玻片组成的采样空间，电磁铁控制载玻片和盖玻片的开闭）。上电时，遮光板前移，载玻片和盖玻片相互接触，断电时，遮光板后移，载玻片和盖玻片分离。采样完毕后，开启抽气泵，将密闭空间中的采样含尘气体经排气管排出密闭空间，完成1次更新粉尘的过程。

1.3粉尘样本图像采集装置

在粉尘样本图像采集装置中，主要是通过工业相机及镜头相关硬件进行图像采集，因此需参考粉尘粒径范围、最小分辨率等物性参数选取合适的图像传感器、镜头及打光方式。

综合本文检测需求及粉尘颗粒的相关物性参数，本文选择的工业相机图像传感器为 CMOS，靶面尺寸为1/2.5英寸，分辨率为2592×1944，像元尺寸为2.2μm ×2.2μm。选用 C 型接口，镜头放大倍数为10，且满足相机的分辨率要求。选取 LED 为光源镜头，采用环形正面打光方式。

在图像处理中，通常使用像素长度来测量图像中物体的尺寸，但要得到真实的物理尺寸需进行相机标定。经标定，本文所构建的图像采集系统中转换系数为0.00235 mm/pixel。

1.4图像处理

由于获取到的原始图像噪声分布明显，粉尘特征分布并不显著，所以需对原始图像进行处理，如图2所示。首先，采用非局部均值（Non?Local Means， NLM）滤波算法对获得的粉尘颗粒图像进行去噪，同时最大程度保持图像的细节特征。然后，采用多尺度 Retinex算法对图像进行增强，提高图像中粉尘颗粒和背景的对比度及各参量的准确度，以解决获取的图像存在亮度差异、不能有效显示因光照不均匀而淹没粉尘颗粒细节的问题。最后，运用分水岭算法对图像进行分割，将图像中的粉尘与背景分开，分割出图像中的实测粉尘颗粒对象。

由图2可看出，与原始灰度图像相比，经过NLM 算法去噪后的图像去噪效果明显，纹理保持得非常理想，且粉尘边缘更加光滑自然。经过 Retinex 增强算法处理后的粉尘图像明显改善了光照不均匀的问题，为后续图像的正确分割提供了可靠的保障。经过分水岭算法处理后，粉尘与背景部分分离，可准确分割出实测粉尘颗粒对象。

2粉尘图像的特征提取

从粉尘质量浓度及粒径分布检测的实际研究内容出发，对经过处理的粉尘样本图像进行特征提取。

2.1粉尘图像的灰度特征提取

由于粉尘颗粒样本与背景存在灰度差异，全局灰度可在一定程度上反映粉尘质量浓度，即全局灰度较低的样本可能含有较多粉尘颗粒。所以，本文选取了粉尘样本图像的灰度均值作为灰度特征，分析灰度均值与粉尘质量浓度之间的关系。

灰度直方图描述的是图像中各灰度级像素个数，反映图像中各灰度等级出现的频率。粉尘样本图像经过预处理后的灰度直方图如图3所示。

灰度均值 m 反映了一幅图像中的平均灰度值，其计算方法为

式中：k 为图像中所有可能的灰度值；L 为灰度等级； zk 为各个灰度级；p（zk）为各个灰度级出现的概率。

粉尘样本经过预处理后的平均灰度特征曲线如图4所示，可看出100张图像的灰度均值的范围主要集中在94～200。

2.2图像的纹理特征提取

因本文研究的粉尘样本颗粒图像中蕴含着丰富的细节信息，通过计算每个像素与周围像素的灰度值出现频率来提取图像的纹理信息，并通过计算粉尘样本颗粒图像纹理特征的灰度共生矩阵来获取其纹理特征。

1）纹理对比度表示整幅图像中像素和其相邻像素之间的亮度对比值，反映图像的清晰度和纹理沟纹深浅的程度。纹理对比度值越大表明图像越清晰，纹理沟纹越深；纹理对比度越小则表明纹理沟纹浅，图像视觉效果模糊。

Con =（i-j）2 F（i; j） （2）

式中：Con 为纹理对比度；i; j分别为图像的行像素点和列像素点；F（i; j）为灰度共生矩阵。

粉尘样本经过预处理后的纹理对比度特征曲线如图5所示，可看出100张图像的纹理对比度主要处于6～8。

2）纹理同质性表示灰度共生矩阵中元素相对于灰度共生矩阵对角线分布的紧密度，反映图像纹理区域变化。纹理同质性的值大则图像纹理各区域间变化很小，局部均匀；纹理同质性的值小则图像纹理局部变化大，局部不均匀。

IDM =  （3）

式中 IDM 为纹理同质性。

粉尘样本经过预处理后的纹理同质性特征曲线如图6所示，可看出100张图像的纹理同质性的值主要集中在0.3～0.5。

3）纹理相关性表示整幅图像中像素与其相邻像素之间的互相关度。纹理相关值大，则灰度分布均匀；相反，则灰度分布不均。

式中：Cor 为纹理相关性；mr 为沿归一化的各行计算的灰度均值；mc 为沿各列计算的灰度均值；δr 和δc 分别为沿各行和各列计算的标准差。

粉尘样本经过预处理后的纹理相关性特征曲线如图7所示，可看出100张图像的纹理相关性主要集中在0.5～1。

4）角二阶矩描述纹理特征的均匀与粗细程度，角二阶矩越大，图像的纹理分布越均一。

式中 ASM 为角二阶矩。

粉尘样本经过预处理后的角二阶矩特征曲线如图8所示，可看出100张图像的角二阶矩均处于0～0.7，绝大多数小于0.1。

2.3图像的几何特征提取

图像的几何特征指图像中目标区域的周长、面积、长轴、短轴等相关特征。为了更好地探究粉尘图像的特征与粉尘质量浓度之间的关系，提取通过图像分割得到的粉尘像素数量与整体图像像素数量的比值作为研究粉尘质量浓度的几何特征。

C = A/S   （6）

式中：C为粉尘像素数量与整体图像像素数量的比值；A为粉尘颗粒对象所占像素点的个数；S 为整体图像像素数量。

粉尘样本经过预处理后的粉尘像素数量与整体图像像素数量的比值几何特征曲线如图9所示，可看出100张图像中粉尘像素数量与整体图像像素数量的比值主要集中在0～30%。

2.4图像获取时粉尘实测浓度

采集图像的同时，采用标准粉尘质量浓度检测仪器实测掘进工作面的粉尘质量浓度，得到100张图像获取时的实测粉尘质量浓度，記录每个图像所对应的实测粉尘质量浓度的变换曲线，如图10所示。可看出最大实测粉尘质量浓度为296.5 mg/m3，最小为76.5 mg/m3。

3.1基于皮尔逊相关系数的相关性分析

皮尔逊相关系数常用来反映两变量相关关系密切程度，本文选用协方差矩阵方法研究实测粉尘质量浓度与图像的灰度、纹理、几何特征之间的相关性。相关性分析结果见表2。

由表2可看出，灰度均值、粉尘像素数量与整体图像像素数量的比值与实测粉尘质量浓度的相关性大，且灰度均值与实测粉尘质量浓度呈负相关，与粉尘像素数量与整体图像像素数量的比值呈正相关。纹理相关性和纹理对比度与实测粉尘质量浓度呈正相关，且相关性较大。纹理特征中的角二阶矩、纹理同质性与实测粉尘质量浓度的相关性较小，且与实测粉尘质量浓度呈负相关。故本文选取灰度均值 m、粉尘像素数量与整体图像像素数量的比值 b、纹理相关性 n、纹理对比度 r 四个特征参数，分析其与实测粉尘质量浓度 Q 的内在规律。

3.2粉尘质量浓度与相关参数多元回归模型的建立及求解

将实测粉尘质量浓度 Q 作为最终参数，选取灰度均值 m、粉尘像素数量与整体图像像素数量的比值 b、纹理相关性 n、纹理对比度 r 作为自变量，并将以上图像特征的范围作为约束条件建立数学模型。

结合优化模型与样本数据，建立图像特征参数中的灰度均值 m、粉尘像素数量与整体图像像素数量的比值 b、纹理相关性 n、纹理对比度 r 与实测粉尘质量浓度 Q 的多元线性回归模型和多元非线性回归模型。

多元线性回归模型为

多元非线性回归模型为

表3为线性回归模型和非线性回归模型方差分析结果。 F 值为显著性差异水平，P 值为检验水平。 F 值越大，说明变量对于模型的影响就越显著，P 值的显著水平小于0.05时，认为该变量对于模型有显著水平。R 为输入变量对输出变量的解释程度，值越大，说明拟合程度越好。由表3可看出，多元线性和多元非线性回归模型的 P 值均小于0.05。但多元非线性回归模型的 F 值和 R 均優于多元线性回归模型。因此本文选取多元非线性回归模型作为粉尘质量浓度检测的数学模型。

3.3粉尘的粒径分布

颗粒的粒径为颗粒所占空间大小的尺度，需选取合适的粒径定义方式。本文选用几何当量等效面积径来表征粉尘的粒径分布。选用直方图来显示不同区间范围内粒径的分布数量。结合实际图像中的粉尘颗粒数目，选取10μm作为粉尘粒径区间长度，计算不同区间范围内粒径大小的占比情况。采用累计分布来说明粒径区间小于或等于某一特定粒径的整体概率。

选取经过预处理后的1张图像来获取粉尘对象面积。将得到的对象像素个数与实际的尺寸进行转化得到实际的面积尺寸，根据计算公式得到粉尘等效面积径。

式中：DV 为几何当量等效面积径；V为颗粒的面积。

粉尘各粒径区间范围内的分布个数及区间占比的累计分布如图11所示。可看出粒径在60μm 以下的粉尘样本占整体分布数量的80%以上。其中粒径在[30，40）μm 范围占比最多，为14.26%。在[90， 100）μm 范围内占比最少，为1.79%。同时，粒径在100μm 以上的粉尘样本占整体分布数量的7.38%。

4实验验证

4.1实验平台搭建

为验证粉尘图像特征参数与实测质量浓度及粉尘粒径分布模型的准确性，在前期搭建的煤矿局部通风设备智能调节测试平台基础上，加入粉尘样本和图像采集装置，如图12所示。为减小实验的随机误差，采用多次测量求平均值的方法。

4.2粉尘质量浓度测试与验证

通过粉尘图像获取相关特征参数，结合建立的数学模型，得到粉尘质量浓度与粉尘粒径分布，同时采用标准仪器检测实验平台巷道内实际的粉尘质量浓度与粒径分布，粉尘实测质量浓度与模型计算浓度对比见表4。

由表4可看出，经过100次实验所测的实测粉尘质量浓度与多元非线性回归模型数学模型计算值的平均相对误差为12.37%，处于合理区间范围内，验证了多元非线性粉尘质量浓度数学模型的有效性和准确性。

4.3粒径分布测试与验证

经过5组粒径分布对照实验（图13），可看出标准实测数据与用几何当量等效面积径得到的粒径分布间的最大相对误差为8.63%，平均相对误差为6.37%，处于合理区间范围内。验证了粒径分布数学模型的有效性和准确性。

5结论

1）通过皮尔逊相关系数对实测粉尘质量浓度与图像特征进行相关性分析，验证了实测粉尘质量浓度与灰度均值、粉尘像素数量与整体图像像素数量的比值、纹理相关性、纹理对比度之间存在强相关性。

2）建立实测粉尘质量浓度与图像特征间的多元非线性回归模型，通过计算的粉尘质量浓度与实测粉尘质量浓度之间对比分析，得到两者间的平均相对误差为12.37%。

3）根据几何当量等效面积径理论，计算出粉尘粒径分布结果，并与实测粒径分布范围进行对比，区间最大相对误差为8.63%，平均相对误差为6.37%。

参考文献（References）：

[1]李德文，赵政，郭胜均，等.“十三五”煤矿粉尘职业危害防治技术及发展方向[J].矿业安全与环保，2022，49（4）：51-58.

LI Dewen，ZHAO Zheng，GUO Shengjun，et al.'13th Five-Year Plan' coal mine dust occupational hazard prevention and control technology and development direction[J]. Mining  Safety & Environmental Protection，2022，49（4）：51-58.

[2]王双明，申艳军，宋世杰，等.“双碳”目标下煤炭能源地位变化与绿色低碳开发[J].煤炭学报，2023，48（7）：2599-2612.

WANG Shuangming， SHEN Yanjun， SONG Shijie， et al. Change of coal energy status and green and low- carbon development under the 'dual carbon' goal[J]. Journal of China Coal Society，2023，48（7）：2599-2612.

[3]程卫民，周刚，陈连军，等.我国煤矿粉尘防治理论与技术20年研究进展及展望[J].煤炭科学技术，2020，48（2）：1-20.

CHENG Weimin，ZHOU Gang，CHEN Lianjun，et al. Research progress and prospect of dust control theory and technology in China's coal mines in the past 20 years[J]. Coal Science and Technology，2020，48（2）：1-20.

[4]袁亮.煤矿粉尘防控与职业安全健康科学构想[J].煤炭学报，2020，45（1）：1-7.

YUAN Liang. Scientific conception of coal mine dust control and occupational safety[J]. Journal of China Coal Society，2020，45（1）：1-7.

[5]陈继民，陈鹤天.激光在粉尘检测领域的进展与应用[J].应用激光，2018，38（3）：496-501.

CHEN Jimin，CHEN Hetian. The development and application of laser in dust detection[J]. Applied Laser，2018，38（3）：496-501.

[6]王智超，吳占松，杨英霞，等.粉尘质量浓度测试方法的实验研究[J].清华大学学报（自然科学版），2013，53（3）：366-370.

WANG Zhichao，WU Zhansong，YANG Yingxia，et al. Experimental study of test method for dust mass concentration[J]. Journal of Tsinghua University （Science and Technology），2013，53（3）：366-370.

[7] HALTERMAN A， SORSAN A， PETERS T M. Comparison of respirable mass concentrations measured by a personal dust monitor and a personal dataRAM to gravimetric  measurements [J]. Annals  of  Work Exposures and Health，2017，62（1）：62-71.

[8]王凯，张霞，苏瑾.地铁站中 PM（2.5）和 PM（10）光散射法与滤膜称重法比对[J].实用预防医学，2021，28（12）：1539-1542.

WANG Kai，ZHANG Xia，SU Jin. Comparison of PM（2.5） and PM（10） mass concentration in subway measured by light scattering method and gravimetric method[J]. Practical Preventive Medicine，2021，28（12）：1539-1542.

[9]简亚婷，高军，肖卫民，等.微量振荡天平法确定激光粉尘仪大气颗粒物转换系数[J].井冈山大学学报（自然科学版），2014，35（3）：1-6.

JIAN Yating，GAO Jun，XIAO Weimin，et al. Using tapered element oscillating microbalance method to determine the conversion coefficient of atmospheric particulate matter tested by laser dust monitor[J]. Journal of Jinggangshan University（Natural Science），2014，35（3）：1-6.

[10]靳燕，王晓东.基于β射线和光散射法的颗粒物测量系统[J].山西电子技术，2022（4）：27-28，35.

JIN Yan，WANG Xiaodong. An atmospheric particulate matter measurement system based on β-ray and light scattering method[J]. Shanxi Electronic Technology，2022（4）：27-28，35.

[11]陈建阁，李德文，许江，等.基于光散射法无动力粉尘质量浓度检测技术[J].煤炭学报，2023，48（增刊1）：149-158.

CHEN Jiange，LI Dewen，XU Jiang，et al. Detection technology of unpowered dust concentration based on light scattering method[J]. Journal of China Coal Society，2023，48（S1）：149-158.

[12]颜鸽来，刘国庆.基于主动荷电的电荷感应粉尘浓度检测技术[J].电子设计工程，2023，31（10）：111-116.

YAN Gelai， LIU Guoqing. Charge-induced dust concentration detection technology based on active charging[J]. Electronic Design Engineering，2023，31（10）：111-116.

[13]丁喜波，王珊，邹一风.基于超声波声衰减的甲烷浓度测量方法[J].应用基础与工程科学学报，2018，26（3）：672-679.

DING Xibo，WANG Shan，ZOU Yifeng. The method of methane concentration measurement based on ultrasonic sound attenuation[J]. Journal of Basic Science and Engineering，2018，26（3）：672-679.

[14]陈锋，夏凤毅，罗爱爱.基于图像法的粉尘浓度检测[C].《环境工程》2019年全国学术年会，北京，2019：179-183.

CHEN Feng， XIA Fengyi， LUO Aiai. Dust concentration detection based on image method[C].2019 National Academic Annual Conference by Environmental Engineering，Beijing，2019：179-183.

[15]劉丹丹，侯展博，胡丙阳，等.基于 CCD 图像的粉尘浓度面维度分布测量方法[J].黑龙江科技大学学报，2023，33（1）：99-102.

LIU Dandan， HOU Zhanbo， HU Bingyang， et al. Surface dimension distribution measurement method of dust concentration based on CCD image[J]. Journal of Heilongjiang University of Science and Technology，2023，33（1）：99-102.

[16]叶方平，方朝阳，徐显金，等.基于图像透光率的粉尘浓度测量算法研究[J].应用光学，2022，43（3）：496-502.

YE Fangping，FANG Chaoyang，XU Xianjin，et al. Dust concentration measurement algorithm based on image transmittance[J]. Journal of Applied Optics，2022，43（3）：496-502.

[17]刘伟华，袁伟.基于图像分析的煤尘在线监测系统[J].电子测量技术，2010，33（1）：125-128，146.

LIU Weihua，YUAN Wei. Investigation of on-line monitoring instrument used in coal dusts based on image analysis[J]. Electronic Measurement Technology，2010，33（1）：125-128，146.

[18] KONOSHONKIN A，KUSTOVA N，BOROVOI A， et al. Using the physical optics approximation for estimating the light scattering properties of large dust particles forlidar applications[C]. The 29th Interational Laser Radar Conference，2020. DOI：10.1051/J.EPJCONF/202023708025.

[19] BAAGA D，KALITA M，DOBRZANIECKI P，et al. Analysis and forecasting of PM2.5，PM4，and PM10 dust concentrations，based on in situ tests in hard coal mines[J]. Energies，2021，14. DOI：10.3390/en 14175527.

[20] BADDOCK M C，STRONG C L，KEYS J F，et al. A visibility and total suspended dust relationship[J]. Atmospheric Environment，2014，89：329-336.