Next Article in Journal
A Study of Learning Issues in Feedforward Neural Networks
Next Article in Special Issue
Prediction of Uniaxial Compressive Strength in Rocks Based on Extreme Learning Machine Improved with Metaheuristic Algorithm
Previous Article in Journal
Lie Symmetry Classification and Qualitative Analysis for the Fourth-Order Schrödinger Equation
Previous Article in Special Issue
Research on b Value Estimation Based on Apparent Amplitude-Frequency Distribution in Rock Acoustic Emission Tests
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pollutant Migration Pattern during Open-Pit Rock Blasting Based on Digital Image Analysis Technology

1
School of Resources and Safety Engineering, Central South University, Changsha 410083, China
2
State Key Laboratory of Coal Mine Safety Technology, China Coal Technology & Engineering Group Shenyang Research Institute, Shenfu Demonstration Zone, Shenyang 113122, China
*
Authors to whom correspondence should be addressed.
Mathematics 2022, 10(17), 3205; https://doi.org/10.3390/math10173205
Submission received: 21 July 2022 / Revised: 25 August 2022 / Accepted: 1 September 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Mathematical Problems in Rock Mechanics and Rock Engineering)

Abstract

:
Previous studies have revealed that toxic gases and dust (smoke dust) are the most common pollutants generated by the blasting operations in open-pit mines, which might lead to a threat to the environment’s condition, health and safety, and properties protection around the blasting site. In order to deal with the problems, a pollution evaluation system was established based on the fractal dimension theory (Dbox(P)) and grayscale average algorithm (Ga) in digital image-processing technology to recognize and analyze the distributions of the smoke-dust cloud, and subsequently determine the pollution degrees. The computation processes of Dbox(P) and Ga indicate three fitted correlations between the parameters and diffusion time of smoke dust. Then, a pollution index (Pi) is put forward to integrate the global and local features of Dbox(P) and Ga, and develop a hazard classification mechanism for the blasting pollutants. Results obviously denote three diffusion stages of the pollutants, mainly including generation stage, cloud-formation stage, and diffusion stage. In addition, it has been validated that the proposed system can also be utilized in single-point areas within a whole digital image. Besides, there are variation trends of the thresholds T1 and T2 in binarization with the diffusion of pollutants. With this identification and evaluation system, the pollution condition of smoke dust can be obviously determined and analyzed.

1. Introduction

With the evolution in industrial fields, the mining industry is confronted with a huge energy demand, and is gradually emerging as the critical support for technological development and progress [1,2,3]. To adapt to community challenges, ore-deposit extractions are constantly expanding, especially for open-pit mines [4,5]. Blasting is a conventional excavation measure for the ores crushing on the ground, and also an option for constructing roads, tunnels and other constructions in mine projects [6]. However, during blasting operations, there are undesirable consequences engendered, such as seismic waves, explosive flying scatter (fly rocks), air shockwave, noise, toxic gases, blasting dust, etc. [6,7].
Within these hazards, blasting dust, the PM2.5 and PM10, and some toxic gases, including CO, NOx, SO2, H2S, and so on, become prominent components of environmental pollutants in blasting sites [8]. Most of these toxic gases are released in a range of concentrations due to environmental and technical factors [9]. Inhaling a certain amount of these toxic gases during their diffusion period can be fatal [10]. Furthermore, prior studies have denoted that it can reach thousands of milligrams per cubic meter for the local blasting dust concentration in surface mining [11]. The combination of blasting dust with toxic gases, which is collectively referred to as smoke dust, has attracted extensive attention owing to its salient hazardousness on air quality in mining sites and detrimental effects on human health [12,13]. It has been demonstrated that a fairly high percentage of miners suffer from respiratory diseases or reduction in respiratory-system immunity due to long-term exposures to high-concentration dust environments, especially pneumoconiosis [14,15]. In addition, the accumulation of dust reduces the humidity of the atmosphere, which then increases the ambient temperature [16]. Meanwhile, it is noteworthy that smoke dust reduces visibility in the mining site, resulting in high working risks and a significant reduction in production efficiency [17]. Large amounts of dust pollution can depose on the ground, causing vegetation degradation and animal habitat destruction [18]. With the characteristic of high concentration and wide range of diffusion it might spread into the water area for dust clusters, thereby contaminating the downstream areas [8]. Studies over the past decades have provided important information on the generation and diffusion mechanisms of smoke dust in open-pit mines. Numerical simulations, for instance, Fluent software, LS-DYNA etc., were conducted to figure out the smoke dust movement, diffusion rate, and the time consuming to float out of the pit, validating the effects of vortex action and the vertical height in the pit on coal-mine dust, and optimizing blasting parameters at the same time [19,20]. Inspired by the approaches, Tang et al. [21] established a circulating accumulation emission model to analyze dust diffusion, including dust emission, retention, and diffusion, respectively. Besides, pertinent artificial intelligence systems, such as particle swarm optimization (PSO) and the long short-term memory (LSTM), were employed to predict and control dust concentration [22,23,24]. In terms of the evaluation system of efficiency of dust reduction with multiple techniques, there are also studies [25,26] in virtue of some tools (e.g., computational fluid dynamics model (CFD): ANSYS/Fluent 10.0 software, ANSYS, Inc., Pittsburgh, Pennsylvania, USA) mainly aiming at monitoring the dust diffusion, and then guaranteeing the mitigation of the smoke-dust pollution.
Furthermore, there are increasing concerns that, exposed to high-concentration smoke dust in open-pit mines, measurement equipment has been a crucial factor for immediately monitoring the pollution condition. Attalla et al. [9] adopted non-dispersion infrared (NDIR) and differential optical absorption spectroscopy (DOAS) to examine the diffusion of fugitive NOx and other pollutant gases. In addition, dust sample sensors, which are commonly based on ultrasound, optical, and electrochemical sensing elements, have also emerged as conventional methods for assessing contaminant gases and particulate fractions of dust emission [2]. With the widespread application of unmanned aerial vehicles (UAVs), researches were performed to load a series of sensors in the rotary-winged platforms, such as CO2, NO2, SO2, PM10, etc. [27]. When in the application process, it is essential for the sampling collection location and source scales to be considered owing to the differences in fixed-wing and multi-rotor UAVs, and thus the long response time, limited measurement distance, and the accuracy of sensing elements hindered further use of UAVs.
Recent developments in the field of machine vision have led to a renewed interest in various industries due to their merits in non-contact measurements, full-visualization monitoring, and affordable cost [28,29,30,31]. In spite of the fact that the related studies on these vision-based models are currently at an early stage, digital image technology, mainly based on grayscale features and morphological features, has been widely applied in pollution analysis and control [32,33,34,35]. Dust-pollution measurement models established on the basis of grayscale features attempted to develop a relationship between the dust concentration and grayscale features extracted from dust images. Yue et al. [30] presented a photo-based method for analyzing air pollution using grayscale data. Thereafter, a logarithmic relation was acquired between the grayscale values and the powder volume by Grasa et al. [34], providing a calibration curve without any fitting constants. On the other hand, morphological feature models are inherently associated with the connections between the morphological features (such as the particle shape, size distribution, and surface structure) and dust concentration. R. Davies [36] and M. Taylor [37] were convinced that parameters of concentration and morphology were important in particle detection, and the latter presented mathematical models to characterize the particles. M. Paphael et al. [38] estimated the mean particle size utilizing optical parameters, on-line turbidity, and solid concentrations. It is also notable for the fractal dimension, another aspect of morphological features, to describe the space occupancy and the irregularities of granulated materials [39,40,41,42]. To figure out the effects of fractal dimension on the porous texture of coal dust images, Mahamud et al. [41] conducted experiments to measure particles exposed to different temperatures. Zhang et al. [42] apply fractal dimension to record the particle size of crushed specimens and obtain a legible regularity of size distribution in different stress conditions. Nevertheless, traditional grayscale or morphological models only considere the local features of pollutants but ignore the global features of their distributions, resulting in a one-sided evaluation of the pollution environment.
To further present the main advantages and limitations of the technology and equipment of previous studies on the pollutant migration that happened during rock-blasting operations, an overall comparation between them is summarized as the Table 1:
Overall, all the evidence reviewed here seems to require a comprehensive evaluation system and distribution analysis for smoke dust pollution caused by blasting operations in open-pit mines.
In this paper, a vision-based model was established based on the fractal dimension theory (Dbox(P)) and a grayscale average algorithm (Ga) for smoke dust evaluation in open-pit mines. Framed pictures were first collected from the video recorded by a UAV, and then, through a series of image processing operations, Dbox(P) and Ga were separately calculated, the results of which were fitted by logarithmic and linear functions, respectively. For the sake of considering both global and local characteristics of pollutants, and fully adopting simplified expressions to characterize the various generation stages of smoke dust in blasting sites, a pollution index (Pi) was proposed to integrate the Dbox(P) and Ga, developing a hazard classification mechanism corresponding to four-level regions. With this evaluation system, the pollution condition of smoke dust can be obviously developed and analyzed, and eventually adequate preparations in response time and emergency preparedness are capable of eliminating detrimental effects on health protection and property safety, as well as improving the environmental quality around the blasting sites.

2. Image Processing Method Based on Grayscale Average and Fractal Dimension

Digital image processing has been successfully applied in the engineering field depending on the difference in the morphological characteristics or the color features of the research object against the homologous background [44,45,46]. Based on the technical principles, and the available devices, it is possible to achieve the monitoring and evaluation of the smoke-dust. The specific procedures of the algorithm in this research can be illustrated as a flowchart in Figure 1:
As shown in Figure 1, the proposed method can be divided into two parts, including the modules of image-processing operations and properties analysis, which primarily aim to acquire the feature image of the pollutants, and relevant calculation results about their distribution properties. According to the flowchart, pollution conditions in blasting sites can be analyzed and evaluated through a series of digital images. The specific application procedures are briefly presented below:
Step one: collect pollutants images of smoke dust at the blasting site, and then the software reads the images to form a three-dimensional digital matrix Cm×n×3, which includes pixel values c i j k and homologous locations ( i , j ); m, n represent the image size in height and width, and the third dimension of the matrix represents three color channels: R, G and B, respectively.
Step two: based on the image matrix Cm×n×3, extract the main regions where pollutants distribute via Equations (1) and (2) below in accordance with the color characteristics of the smoke dust [47]. As is usual, the smoke dust pollutants in non-coal mines appear a brighter color, and thus, compared to several dark backgrounds, the pixel values of pollutants tend to be higher. Through the screening method, a series of pixels representing superfluous regions will be removed, whereby the images can be prepared for subsequent image processing operations.
C s u m = c i j 1 + c i j 2 + c i j 3
C s u m > a
Step three: fulfill further image-processing operations with pertinent programs. Image graying is enforced using the empirical formula shown in Figure 1, the matrix result of which is denoted as Gm×n, and by virtue of image graying, the color dimensions are converted into two dimensions, viz., the image size. Then, binarization was served to recognize the regions of interest, in which a double-threshold method (T1, T2) was put forward to improve the performance of the feature extractions, and in this way, pixels within the gray image that are, respectively, lower or higher than T1 and T2 are assigned as 0 (black in color), while the intermediate areas between the thresholds are assigned as 1 (white in color). Noise reduction is then adopted to eliminate excrescent noise points caused by the forming transversion of pollutant images [48].
Step four: calculate the Dbox(P) and Ga of the smoke-dust images to quantitively analyze the distribution changes of on-site pollutants. In the process, the two-dimensional box counting method was developed to acquire the fractal property of the pollutants [40]. Within the fractal calculation processes, the binary images are traversed by a slew of windows, the side lengths of which are decided by the original picture size, and are incremental pixel by pixel until the width and height reach half of the original length. Then, for each window, keep a track of the number of square boxes N(ε) intersecting the pollutant parts as well as the lengths of the boxes ε. Afterwards, instead of the mathematical definition of fractal dimension theory as Equation (3), a regression method is used to analyze two groups of the box data through logarithmic processing, which can be expressed as Equation (4). Obviously, the slope of the regressed equation Dbox(P) is the fractal dimension of the smoke dust.
D b o x ( P ) = l i m 0 l o g   N ( ) l o g   1 /
l o g   N ( ) = D b o x ( P ) l o g   1 /   + C   ( C   is   a   regressed   constant )
Meanwhile, Ga of the binary images is calculated based on Equation (5), which can be defined as the division operation between the pixels denoting pollutants and site background, respectively. When in the calculation process, a grayscale distribution curve is plotted to develop the pollutants distribution along the directions of horizontal strike.
G a = i = 1 h j = 1 w b i , j h × w
where h and w represent the height and width of smoke-dust images.
Step five: through the computation of the Dbox(P) and Ga of the pollutant, a pollution index Pi is provided to achieve a comprehensive consideration of Dbox(P) and Ga, and the integration expression can be depicted as Equation (6):
P i = D b o x ( P ) + G a
Then, according to the change tendency of Pi, it can be more accurate to evaluate the pollution condition from a quantitative aspect, and with a series of pollutant images collected from a whole blasting project, the index Pi is capable of dividing into four regions of pollution levels, Level I~IV, viz. Therefore, the pollutant condition of smoke dust can be dramatically characterized no matter the period of monitoring or throughout the whole working process by this research.

3. Algorithm Application of Dbox(P) and Ga with a Fugitive Dust Picture

For the sake of further demonstrating the proposed evaluation method of pollutant distribution in a blasting site, an image sample with fugitive dust shown as Figure 2 was applied to validate the effectiveness of the image-processing algorithm, and measurement of Dbox(P) with Ga. Firstly, the original color image of fugitive dust as Figure 2a was read, where in the image dust particles emerge distinctly with homologous pixels when compared to the black background. With the application of the image-processing method, the region of the floating dust was arguably extracted to produce a binary image, which can be visualized as Figure 2b. Based on the binary image, the calculation processes of Dbox(P), Ga are also displayed in Figure 2c,d, respectively. As for the grayscale distribution, an ergodic window of 438 × 10 pixels in size was employed to cover the binary image along the direction of the horizontal, and with the traversal calculations, the grayscale distribution, cumulative distribution curve, and grayscale average were obtained, the results of which are displayed in Figure 2e. According to the grayscale parameters, Ga of the whole image is 0.2490, and the dust objects primarily mainly distribute at 220~550 pixels of the horizontal strike. Similarly, the calculation procedures of Dbox(P) are devised as Figure 2d using the two-dimensional box counting method. With a series of measuring boxes covering the binary image, the computation procedures are revealed as the fitted equation in Figure 2f, and the Dbox(P) of the dust image is 1.977, the R2 of which is 0.9998, embodying a precise fitting effect when applying it to research objects of fine grains.
In order to verify the accuracy of calculation results in Dbox(P) and Ga with the proposed binary dust image handled from the original image, the grayscale average distribution of the gray image, remaining quite complete image information with a 256-level grayscale form, was devised to analyze the change tendency of dust pixels as a comparison to the binary image, and quantitatively appraise the effectiveness on object segmentation within the application of image processing. The analysis results can be shown as Figure 3.
As is mentioned with the grayscale parameters, dust particles are primarily distributed at 220~550 pixels of the horizontal strike in accordance with the binary image. It can be perceived that the developed trend of the grayscale average from gray image tends to be consistent, performing a fairly smooth waveform shape. In addition, it is almost coincident with the high-distribution region of dust particles in any binary image or gray image. While there is a hysteresis effect for the grayscale parameter of gray image at the boundaries between research objects and background, producing two hysteresis intervals in the rising and falling stages of the grayscale average, the area of the hysteresis intervals accounts for about 5.58% of the total enclosed area by the grayscale distribution line of gray image. It is imperative to note the significance turning points, data turning point G1 (90, 25.16) of gray image embodying a hysteresis effect contrasting to B1 (210, 7.37 × 10−5) of binary image for the first hysteresis interval, whereas G2 (560, 51.10) illustrates an almost the same horizontal position with B2 (560, 3.08 × 10−5). A possible explanation for this might be the low pixel values of edge dust particles in the image (i.e., grayscale value). As a whole, in spite of the existence of hysteresis intervals, the application image of fugitive dust presents a validity in the object recognition of dust clouds and other fine-grained substances from a quantitative aspect.

4. Results on the Dbox(P) with Ga and the Pollution Evaluation System of Blasting Operations

A field experiment for analyzing smoke dust distribution, and evaluating the pollution condition based on the proposed method, was conducted in a deep open-pit mine with multi-step terrain production in Xinjiang, China. Shallow hole blasting was adopted for the ores crushing along the hillside wall, and the whole blasting operations with partial diffusion were recorded by a UAV platform flying directly above the blasting region with a fixed height. The platform can be predetermined in an autonomous direction, and illustrates high-temporal sampled pictures in real time.
A total of 95 images were extracted from the video recorded by the UAV with an interval of 1 s, which are 640 pixels in width and 368 pixels in height, integrally revealing the processes of blasting preparation, exploding and diffusion of smoke dust, respectively, which are shown in Figure 4.
From the images in Figure 4, it can be found that there is an obvious contrast between smoke dust and the field background throughout the periods of blasting operation. At the initial stage, debris and dirt deposits on the ground were lifted up due to the shock wave of blasting, and with the gravity effect, the smoke dust was then dispersed suffering from the wind floating and other meteorological conditions. At this stage, most distribution regions of the pollutants perform as partially bright in color, while regions representing the ground surface and the bench tend to be dim. For a clearer description, a three-dimensional elevation grayscale image was put forward to show the grayscale distributions and the differences in elevation, as depicted in Figure 5. The elevation range in Figure 5e extracted from the region presents an exceedingly grayscale difference between the pollutants and ground surface. By this method, feature extraction of a complicated environment between the pollutants and background can be fulfilled through a series of image-processing algorithms proposed in this research.
To ulteriorly assess the effectiveness of the proposed method in blasting operations, the 95 images from the blasting work were processed by the image-processing algorithm, and then served to determine the Dbox(P) and Ga. Within the pollutant images, the computation processes begin from the 16th image because the 1st~15th images belong to the prior stage of blasting, and smoke dust is barely distributed in these images. Partial image-processing results of the 15th, 30th, 45th, 60th, 75th and 90th images are shown as Figure 6:
Following the image-processing results it can be discerned that at the first stage of blasting, a large amount of smoke and dust are engendered to form smoke-dust clouds, mainly due to the explosive at each charge point, such as the 30th picture. At this stage, dust particles start to float due to the shockwave caused by the explosion, which is much larger than the gravity of the dust particles, friction between dust particles and rock particles, and adhesion of dust against other items. Then, with the completion of the blasting operation, the smoke-dust clouds gradually develop, when in this process the smoke dust continues to move upward as a consequence of the inertia effect, and the motion track of a large amount of smoke dust is superimposed, subsequently generating a mushroom or mushroom-like cloud. In addition, from a qualitative perspective for the pollutant distribution, the amount of smoke dust appears to have an upward trend with the explosion process, and then enters into a volatility phase, (60)~(75), for instance. Eventually, as the smoke dust continues to spread, which mainly interfered with meteorological conditions in this period, especially for the disturbed wind flow, the distribution is confronted with a certain degree of reduction.
To further describe the amount of pollutant, Dbox(P) and Ga for the 15th~95th smoke dust images were calculated, and the results with corresponding calculating time are demonstrated in Table 2.
According to the results of Ga and Dbox(P), a relatively legible data trend is presented, which can be summarized as increasing to a peaking value, maintaining a stable fluctuation within a certain range, and then decreasing with time, respectively. Moreover, it is found that there is a discrepancy in calculating time between Ga and Dbox(P), the reasons for which may be explained by the varying levels of program complexity.
For the sake of further analyzing the Dbox(P) and Ga of the collected frame pictures from the video footage of field blasting, computations of all images by the proposed method were carried out and the specific results with related fitting curves can be denoted in Figure 7:
It can be easily discovered from Figure 7 that there are characteristic distributions for 80 image samples (16th~95th) on the basis of the results, which accord with a pattern consistent with the qualitative description, mainly including three stages: increasing, stable fluctuations, and decreasing. The fitted equations for the Ga and Dbox(P) are listed in Table 3.
As illustrated in the fitted equations, for both of the Dbox(P) and Ga, stage one (G1 and D1) and stage three (G3 and D3) can be appropriately represented by logarithmic fits. Stage two (G2 and D2) in Figure 7 presents varying data fluctuation within a certain range, and therefore, a linear fit was employed to demonstrate the relationship between the pollutant parameters and dust images at different times.
In addition, from the fit coefficients R2 of all fitted equations, it can be found that no matter Dbox(P) and Ga, there are differences in the logarithmic fit of stage one and stage three. For example, for the grayscale average, R2 of stage one and three are 0.9165 and 0.8188, and as regards the fractal dimension, R2 is 0.5707 and 0.4978, respectively. It is distinct that there is a numerical reduction from stage one to stage three, the possible reasons for which are mainly because of the different causes for the performances of the pollutants. Smoke dust in stage one is most likely caused by a dynamite explosion in the field, whereas smoke dust in stage three is more likely caused by a diffusion effect caused by meteorological conditions, particularly wind floating.
It is also important to denote that, with reference to the fit coefficient, the fitted equations of Ga (the maximum is 0.9165) perform better when compared to the Dbox(P) (the maximum is 0.5707), and this phenomenon elucidates that smoke dust in images reflects more local features than global features because when in the process, Ga is oriented to one calculation for a whole digital image, while Dbox(P) is confronted with multiple identical calculation procedures by a series of detecting boxes covering the image.
Based on the proposals of Dbox(P) and Ga, and while combining the benefits of local and global features, a pollution condition index (Pi) is developed to assess quantitative changes and spatial-temporal distributions in smoke-dust images. By this way, an integration of image parameters consisting of Dbox(P) and Ga can be additionally completed, and in sites, the mathematical expression of Pi is much more widely applied in engineering computation to guarantee acquisitions of pertinent information related to the engineering efficiency and practitioner safety. Pi can be presented as Equation (6) above. On the basis of the equation, Pi of the blasting images can be drawn as Figure 8:
The results of the correlational analysis are presented in Figure 8, and the relevant fitted equations with homologous goodness of fits (R2) are also shown. Compared to any of the parameters in Figure 7 and the calculated results of fitted equations in Table 3, data distributions of Pi similarly illustrate the legible diffusion stages of smoke dust: generation stage, cloud-formation stage and diffusion stage. In the graph, the time intervals of stages one~three are 16~44 s, 45~84 s and 85~95 s, respectively. The single most striking observation to emerge from the data comparison is that the data distribution tendency of Pi is similar to the results of Dbox(P) and Ga. Moreover, the R2 of the fitted curves are 0.8365, 0.0022 and 0.7397, embodying a fair logarithmic relationship for stage one and stage three. At the same time, it also appears to have a more stable data fluctuation for stage two in contrast to the single parameter. In brief, the proposal of Pi illustrates an ideal performance throughout blasting processes, and achieves one mathematical expression corresponding to three stages of pollution distribution. To observe the performance of the proposed method of smoke dust evaluation in this research, comparative results are summarized in Table 4.
Three previous studies on estimation of PM2.5, PM10 and NO2 were presented based on the digital image processing, and RMSE along with R squared values are the metrics for evaluation models. From Table 4, although the three studies targeted three specific pollutants, the RMSE values of the proposed method (stage one and stage three) indicate a fair regression effectiveness when compared to the other these models. Similarly, the R squared values also reveal that the Dbox(P)Ga based method is reliable for the analysis about pollutants migrations in blasting operations. Therefore, the RMSE and R squared values of the proposed method were confirmed by the comparisons.
Based on the comprehensive pollution index of blasting engineering, the proposed model can be further utilized in pollution hazard classifications. In this research, an attempt was made to provide a method to evaluate the pollution condition on site in accordance with the fitting curves shown in Figure 8, which were divided into four pollution levels, including the serious, the high, the middle, and the mild, as depicted in Figure 9. Detailed division procedures are also introduced.
(a)
At the start of the blasting (from the 16th image), according to the fitting result of stage one, the Pi is 1.81, and at stage two of the fluctuation period, the logarithmic fitting results of stage one and three were extended to intersect with it to generate two intersection points, the coordinates of which are (39.71, 2.30) and (85.56, 2.31). Thus, the peak value of the model is dictated as the intersection point between stage two and three. In this way, the Pi interval is [1.81, 2.31].
(b)
After determining the Pi interval, four equal parts were then calculated and determined, and the Pi values along with their corresponding occurrence times are Pi1 (20.14 s, 1.935), Pi2 (25.49 s, 2.06), Pi3 (32.19 s, 2.185) and Pi4 (85.56 s, 2.31), respectively.
(c)
Based on the Pi values, the distribution graph was divided into four regions. As a supplementary, the region below 1.81 of the Pi value was also assigned to the first level, while that above 2.31 was assigned to the last level. Therefore, the whole graph of pollution distribution is divided into four levels, and a color pyramid is used to describe it. The detailed division criteria of level I~IV can be presented as Table 5.
By virtue of the critical factor of Pi, the pollution evaluation system can be established for analyzing the distribution conditions of smoke dust, and provide an approach to develop a comprehensive program in pollution prevention and control. Different from previous methods in monitoring smoke dust pollution, the proposed system in this research adopts a series of digital images, and then fulfills pollution classifications through computations of Pi in the software. Subsequently, the final pollution levels at the target time point can be determined.

5. Discussion

5.1. Pollution Evaluation of Single-Point Area with Image Slices

The aforementioned research on Dbox(P), Ga and Pi were conducted by virtue of certain smoke dust images, in which the computation scope was oriented to the full range of one blasting image. However, sometimes it is necessary to monitor and evaluate the pollutants at a single point within a provided image to figure out whether the proposed method is capable of enforcing the task. Figure 10 demonstrates the application procedures for two single-point areas in the 30th image chosen from the same video.
As shown in Figure 10 b, two pictures were extracted from the binary image, which are the bench and diffusion area, and the selected square regions are 50 × 50 and 70 × 70 in size, respectively. Dbox(P) and Ga of the pictures were then calculated, and the fitted results of fractal dimension are illustrated in Figure 10c,d. Distinctly, Ga of area one and two are 0.174 and 0.481, and Dbox(P) are 1.5454 and 1.8398. Thus, the pollution indexes (Pi values) can be calculated, the results of which are 1.7194 and 2.3208. Depending on the proposed pollution-level evaluation system in the previous section, regions one and two are supposed to be classified as level I and level IV in smoke dust pollution. Consequently, in addition to the computation ability for the full range of an image, the proposed system can be also exceedingly applied in the calculation of a single-point pollution area, thereby satisfying relevant monitoring tasks. At the same time, the results of these two regions also demonstrate that it is feasible for the system to accomplish analysis on the pollution condition of smoke dust through digital images recorded at blasting sites.

5.2. Variation Trends on Binarization Thresholds within Image Processing Operations

In the conversion process from image graying to binarization, double threshold binarization was adopted to extract the pollutant features in an image, the thresholds of which are denoted as T1 and T2. Nonetheless, it is not always consistent for the threshold values with the blasting program in progress. Owing to the pollutants with variable quantities, and the complex background in blasting site, T1 and T2 are also required to adapt to the optimal image conversion in binarization operation. The thresholds used in this research are denoted in Table 6.
It is apparent from this table that there are different numbers of time nodes for T1 and T2 when choosing threshold values in the process of binarization, where T1 consists of five nodes while T2 possesses eight nodes. This result may be explained by the fact that there are more changes happening at the border of the highlight pixels, while the dark regions in the image mostly representing the background present fewer changes. In other words, the generation process of smoke dust is inherently more associated with the bright pixels in the pollutant images.
As for T1, at the first time-interval of 1–15 s, the optimal threshold value tends to be chosen as less than 10, which is an important parameter for the segmentation between the pollutants and background. Thereafter, T1 is preferably assigned as 100 in 16–35 s. Obviously, there is a pretty large increase from the first time-interval to the second, the reason for which might be the mixture of items produced in the initial stage, which mainly consist of mud blocks, gravel, dust particles, and so on, illustrated as grayish-black in morphology. For the sake of capturing the features of items while trying to remove the background part of the picture, it is essential for the threshold border T1 to be moved up to 100.
After the first time-interval, the blasting operation began at 16 s, and then there is a steady decline for T1 value until 49 s with a range of 20. A possible explanation for this might be that the generation process of pollutants has a relatively stable impact on image pixels in dark regions. The finding provides some supports for the conceptual premise of the generation stage aforementioned in the presentation of image parameters due to the three time nodes being from 16 s to 48 s, where the time interval defined before the generation stage is 16–44 s. Moreover, at the time interval of 49–95 s, it appears appropriate for the T1 to be less than 40 so as to preserve dark-area features of pollutants, and the implication of this is the possibility that there is also a minor effect on T1 for the pollutant distribution within the cloud-formation stage and diffusion stage. Potentially, this phenomenon indicates that in these two stages, there are few variations in the dark features of smoke dust.
Turning now to the T2, it is obvious for it to have a data trend of decreasing and increasing in succession (>240–140–>170), the reasons for which could be attributed to the generation of a large amount of smoke dust, inducing the bright features to tend towards unified in a blasting image, and simultaneously masking the features of other bright areas in sites so that the grayscale values of bright pixels experience a decrease. Then, with the fading of time, the diffusion of smoke dust is confronted with a gradual reduction, and the grayscale distribution of the original image progressively reappears again.
At the pre-blast stage of 1–16 s, T2 was decided as less than 240 to guarantee the removal of pixels not related to pollutants, which indicates that there are more distributions of bright pixels in the background part of the original image. On the other hand, different from the steady declines in T1, T2 emerges at decremental or incremental intervals throughout the whole time-node. Within the time nodes of 17–22 s, 23–31 s and 32–37 s, the decreasing of T2 demonstrates that with the start of blasting, smoke dust gradually covers the bright areas in the background, contributing to a more centralized grayscale distribution. Additionally, the decremental intervals of T2 depict that the generation of smoke dust slows down at a later stage. The finding further supports the logarithmic relationship discussed in the generation stage mentioned above. As for the time nodes from 38–76 s, the selections of T2 also present stable fluctuations similar to the cloud-formation stage, and in this stage, the steady decline may be explained by the inertial motion of smoke dust, resulting in different development tracks. Then, at 77–95 s, T2 are chosen at more than 170, which are in line with the previous description of the diffusion stage, when in this process the smoke dust is mainly subject to meteorological conditions, suffering from a sluggish spread.

6. Limitations and Future Studies

The insights gained from this study may be of assistance to evaluate the smoke-dust pollution caused by blasting operations in virtue of digital image technology, and has gone some way towards enhancing the understanding of the relationship between the Dbox(P) wtih Ga, and smoke-dust distribution in digital images. Nonetheless, there are limitations in this study which could have affected the parameters computation, such as the image quality, mine topography, illumination condition, shooting distance with angle and so on. Therefore, further research could also be conducted to determine the effectiveness of more possible factors existing in blasting sites on the application of the evaluation system.

7. Conclusions

A pollution evaluation system based on digital image-processing technology with the grayscale average algorithm and the fractal dimension theory was proposed to recognize the features of smoke dust, and determine the pollution conditions generated by blasting operations in open-pit mines. By using this method, the dust distribution, grayscale distribution, fitted curves of fractal dimension, grayscale average and pollution index can be characterized. These approaches can quantitatively analyze the smoke dust characteristics by means of a computer. Based on the present study, the following conclusions can be drawn:
(a)
An image example of fugitive dust is employed to illustrate and validate the effectiveness of the algorithm in dust particles, and in spite of the existence of the hysteresis intervals, results of fractal dimension (0.9998 in R2 of the fitted equation) and distributions with grayscale parameters show that the proposed approach successfully identifies the dust objects.
(b)
A total of 95 pictures of pollutants were extracted for the parameter calculations. Within the image-processing operations, the grayscale difference between the smoke dust and background was revealed, which is the essential distinction for feature extraction of smoke dust. The Dbox(P) and Ga results, as well as the related fitted curves, were then obtained, and it is clear that three development stages are required for dust diffusion to emerge. In addition, compared to Dbox(P), Ga demonstrated a better fitted correlation for the distribution of pollutants, reflecting more local features.
(c)
Based on the Dbox(P) and Ga, Pi as well as the fitted results were proposed, simultaneously integrating the global and local features of pollutant images. The Pi denoted three diffusion stages of smoke dust: generation stage, cloud-formation stage, and diffusion stage corresponding to 16~44 s, 45~84 s and 85~95 s, respectively. Then, an evaluation system associated with four levels of pollution conditions was also obtained, of which the Pi values are <1.935, 1.935~2.06, 2.06~2.185 and >2.185. By this way, the final pollution levels of smoke dust can be determined through the provided images.
(d)
The comparative results of RMSE and R squared values show that the proposed method (stage one and stage three) presents a fair performance when compared to the other models. Meanwhile, it also reveals that the vision-based method is reliable for the analysis about pollutants migrations in blasting operations. Therefore, the RMSE and R squared values of the proposed method were confirmed by the comparisons. Therefore, the above improvement significantly evaluates the on-site smoke-dust pollution with a high accuracy, and achieves a comprehensive bivariate index to avoid possible measurement errors due to univariate index.

Author Contributions

Conceptualization, S.W.; Methodology, J.Y. and J.L.; Supervision, S.W.; Data curation, J.Y.; Formal analysis, J.Y.; Funding acquisition, F.T. and S.W.; Validation, F.T. and S.W.; Writing—original draft, J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (No. 52174099), the Natural Science Foundation of Hunan Province (No. 2021JJ30842), the Natural Science Foundation of Liaoning Province (No. 2021-KF-23-01), and the Fundamental Research Funds for the Central Universities of Central South University (No. 2022ZZTS0510).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to acknowledge the editors and reviewers for their invaluable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, T.; Liu, S. The impacts of coal dust on miners’ health: A review. Environ. Res. 2020, 190, 109849. [Google Scholar] [CrossRef] [PubMed]
  2. Miguel, A.; Felipe, G.; Andrew, F.; Ashray, D. Towards the development of a low cost airborne sensing system to monitor dust particles after blasting at open-pit mine sites. Sensors 2015, 15, 19667–19687. [Google Scholar]
  3. Xie, C.; Nguyen, H.; Bui, X.; Nguyen, V.; Zhou, J. Predicting roof displacement of roadways in underground coal mines using adaptive neuro-fuzzy inference system optimized by various physics-based optimization algorithms. J. Rock Mech. Geotech. Eng. 2021, 13, 1452–1465. [Google Scholar] [CrossRef]
  4. Yuan, M.; Ouyang, J.; Zheng, S. Research on ecological effect assessment method of ecological restoration of open-pit coal mines in alpine regions. Int. J. Environ. Res. Public Health 2022, 19, 7682. [Google Scholar] [CrossRef] [PubMed]
  5. Blom, M.; Pearce, A.R.; Stuckey, P.J. Short-term planning for open pit mines: A review. Int. J. Min. Reclam. Environ. 2019, 33, 318–339. [Google Scholar] [CrossRef]
  6. Vasović, D.; Kostić, S.; Ravilić, M.; Trajković, S. Environmental impact of blasting at Drenovac limestone quarry (Serbia). Environ. Earth Sci. 2014, 72, 3915–3928. [Google Scholar] [CrossRef]
  7. Lu, W.; Luo, Y.; Chen, M.; Shu, D. An introduction to Chinese safety regulations for blasting vibration. Environ. Earth Sci. 2012, 67, 1951–1959. [Google Scholar] [CrossRef]
  8. Abdollahisharif, J.; Bakhtavar, E.; Nourizadeh, H. Green biocompatible approach to reduce the toxic gases and dust caused by the blasting in surface mining. Environ. Earth Sci. 2016, 75, 191. [Google Scholar] [CrossRef]
  9. Attalla, M.I.; Day, S.J.; Lange, T.; Lilley, W.; Morgan, S. NOx emissions from blasting operations in open-cut coal mining. Atmos. Environ. 2008, 42, 7874–7883. [Google Scholar] [CrossRef]
  10. Zvyagintseva, A.V.; Sazonova, S.A.; Kulneva, V.V. Analysis of sources of dust and poisonal gases in the atmosphere formed as a result of explosions at quarries of the mining and integrated works. IOP Conf. Ser. Mater. Sci. Eng. 2020, 962, 42045. [Google Scholar] [CrossRef]
  11. Sa, Z.; Li, F.; Qin, B.; Pan, X. Numerical simulation study of dust concentration distribution regularity in cavern stope. Saf. Sci. 2012, 50, 857–860. [Google Scholar] [CrossRef]
  12. Chang, P.; Chen, Y.; Xu, G.; Huang, J.; Ghosh, A.; Liu, W.V. Numerical study of coal dust behaviors and experimental investigation on coal dust suppression efficiency of surfactant solution by using wind tunnel tests. Energy Sources Part A Recovery Util. Environ. Eff. 2021, 43, 2173–2188. [Google Scholar] [CrossRef]
  13. Entwistle, J.A.; Hursthouse, A.S.; Reis, P.A.M.; Stewart, A.G. Metalliferous mine dust: Human health impacts and the potential determinants of disease in mining communities. Curr. Pollut. Rep. 2019, 5, 67–83. [Google Scholar] [CrossRef]
  14. Cui, K.; Shen, F.; Han, B.; Yuan, J.; Suo, X.; Qin, T.; Liu, H.; Chen, J. Comparison of the cumulative incidence rates of coal workers’ pneumoconiosis between 1970 and 2013 among four state-owned colliery groups in China. Int. J. Environ. Res. Public Health 2015, 12, 7444–7456. [Google Scholar] [CrossRef]
  15. Cheng, W.; Yu, H.; Zhou, G.; Nie, W. The diffusion and pollution mechanisms of airborne dusts in fully-mechanized excavation face at mesoscopic scale based on CFD-DEM. Process Saf. Environ. 2016, 104, 240–253. [Google Scholar] [CrossRef]
  16. Kayet, N.; Pathak, K.; Chakrabarty, A.; Kumar, S.; Chowdary, V.M.; Singh, C.P.; Sahoo, S.; Basumatary, S. Assessment of foliar dust using Hyperion and Landsat satellite imagery for mine environmental monitoring in an open cast iron ore mining areas. J. Clean. Prod. 2019, 218, 993–1006. [Google Scholar] [CrossRef]
  17. Ajrash, M.J.; Zanganeh, J.; Moghtaderi, B. Impact of suspended coal dusts on methane deflagration properties in a large-scale straight duct. J. Hazard. Mater. 2017, 338, 334–342. [Google Scholar] [CrossRef]
  18. Conesa, H.M.; Faz, Á.; Arnaldos, R. Heavy metal accumulation and tolerance in plants from mine tailings of the semiarid Cartagena–La Unión mining district (SE Spain). Sci. Total Environ. 2006, 366, 1–11. [Google Scholar] [CrossRef]
  19. Wanjun, T.; Qingxiang, C. Dust distribution in open-pit mines based on monitoring data and fluent simulation. Environ. Monit. Assess 2018, 190, 11. [Google Scholar] [CrossRef]
  20. Jia, Z.; Song, Z.; Fan, J.; Jiang, J.; Guo, S. Numerical simulation study on dust suppression mechanism of burning rock blasting in open-pit mine. Front. Earth Sci. 2022, 10. [Google Scholar] [CrossRef]
  21. Tang, W.; Li, F. A new circulating accumulation emission model for assessing dust emission from open pit mine. Sci. Rep. 2021, 11, 24243. [Google Scholar] [CrossRef] [PubMed]
  22. Bui, X.; Lee, C.W.; Nguyen, H.; Bui, H.; Long, N.Q.; Le, Q.; Nguyen, V.; Nguyen, N.; Moayedi, H. Estimating PM10 Concentration from Drilling Operations in Open-Pit Mines Using an Assembly of SVR and PSO. Appl. Sci. 2019, 9, 2806. [Google Scholar] [CrossRef]
  23. Li, L.; Zhang, R.; Sun, J.; He, Q.; Kong, L.; Liu, X. Monitoring and prediction of dust concentration in an open-pit mine using a deep-learning algorithm. J. Environ. Health Sci. Eng. 2021, 19, 401–414. [Google Scholar] [CrossRef] [PubMed]
  24. Qi, C.; Zhou, W.; Lu, X.; Luo, H.; Pham, B.T.; Yaseen, Z.M. Particulate matter concentration from open-cut coal mines: A hybrid machine learning estimation. Environ. Pollut. 2020, 263, 114517. [Google Scholar] [CrossRef] [PubMed]
  25. Torno, S.; Toraño, J.; Menéndez, M.; Gent, M. CFD simulation of blasting dust for the design of physical barriers. Environ. Earth Sci. 2011, 64, 73–83. [Google Scholar] [CrossRef]
  26. Wang, Z.; Zhou, W.; Jiskani, I.M.; Ding, X.; Liu, Z.; Qiao, Y.; Luan, B. Dust reduction method based on water infusion blasting in open-pit mines: A step toward green mining. Energy Sources Part A Recovery Util. Environ. Eff. 2021, 1–15. [Google Scholar] [CrossRef]
  27. Watai, T.; Machida, T.; Ishizaki, N.; Inoue, G. A lightweight observation system for atmospheric carbon dioxide concentration using a small unmanned aerial vehicle. J. Atmos. Ocean. Technol. 2006, 23, 700–710. [Google Scholar] [CrossRef]
  28. Lu, G.; Yan, Y.; Riley, G.; Bheemul, H.C. Concurrent measurement of temperature and soot concentration of pulverized coal flames. IEEE Trans. Instrum. Meas. 2002, 51, 990–995. [Google Scholar]
  29. Shirmohammadi, S.; Ferrero, A. Camera as the Instrument: The Rising Trend of Vision Based Measurement. IEEE Trans. Instrum. Meas. 2014, 17, 41–47. [Google Scholar] [CrossRef]
  30. Yue, G.; Gu, K.; Qiao, J. Effective and Efficient Photo-Based PM2.5 Concentration Estimation. IEEE Trans. Instrum. Meas. 2019, 68, 3962–3971. [Google Scholar] [CrossRef]
  31. Joshi, M.R.; Nkenyereye, L.; Joshi, G.P.; Islam, S.M.R.; Abdullah-Al-Wadud, M.; Shrestha, S. Auto-Colorization of Historical Images Using Deep Convolutional Neural Networks. Mathematics 2020, 8, 2258. [Google Scholar] [CrossRef]
  32. Liu, C.; Tsow, F.; Zou, Y.; Tao, N. Particle pollution estimation based on image analysis. PLoS ONE 2016, 11, e145955. [Google Scholar]
  33. Zhang, T.; Dick, R.P. Estimation of multiple atmospheric pollutants through image analysis. In Proceedings of the IEEE International Conference on Image Processing (ICIP), Taipei, Taiwan, 22–25 September 2019; pp. 2060–2064. [Google Scholar]
  34. Grasa, G.; Abanades, J.C. A calibration procedure to obtain solid concentrations from digital images of bulk powders. Powder Technol. 2001, 114, 125–128. [Google Scholar] [CrossRef]
  35. Li, Y.; Huang, J.; Luo, J. Using user generated online photos to estimate and monitor air pollution in major cities. In Proceedings of the ICIMCS ‘15: Proceedings of the 7th International Conference on Internet Multimedia Computing and Service, Zhangjiajie, China, 19–21 August 2015; pp. 1–5. [Google Scholar]
  36. Davies, R. Summary of the particle characterization session. Powder Technol. 1996, 88, 191–196. [Google Scholar] [CrossRef]
  37. Taylor, M.A. Quantitative measures for shape and size of particles. Powder Technol. 2002, 124, 94–100. [Google Scholar] [CrossRef]
  38. Raphael, M.; Rohani, S. On-line estimation of solids concentrations and mean particle size using a turbidimetry method. Powder Technol. 1996, 89, 157–163. [Google Scholar] [CrossRef]
  39. Mikula, W.I.F.R. Fractal dimensions of coal particles. J. Colloid Interf. Sci. 1987, 120, 263–271. [Google Scholar]
  40. Pi, Z.; Zhou, Z.; Li, X.; Wang, S. Digital image processing method for characterization of fractures, fragments, and particles of soil/rock-like materials. Mathematics 2021, 9, 815. [Google Scholar] [CrossRef]
  41. Mahamud, M.; López, Ó.; Pis, J.J.; Pajares, J.A. Textural characterization of coals using fractal analysis. Fuel Process. Technol. 2003, 81, 127–142. [Google Scholar] [CrossRef]
  42. Zhang, B.; Liu, W.; Liu, X. Scale-dependent nature of the surface fractal dimension for bi-and multi-disperse porous solids by mercury porosimetry. Appl. Surf. Sci. 2006, 253, 1349–1355. [Google Scholar] [CrossRef]
  43. Zhang, J.; Li, M.; Liu, Z.; Zhou, N. Fractal characteristics of crushed particles of coal gangue under compaction. Powder Technol. 2017, 305, 12–18. [Google Scholar] [CrossRef]
  44. Albatayneh, O.; Forslöf, L.; Ksaibati, K. Developing and validating an image processing algorithm for evaluating gravel road dust. Int. J. Pavement Res. Technol. 2019, 12, 288–296. [Google Scholar] [CrossRef]
  45. Fu, Y.; Wang, N. Measurement of dust concentration based on VBAI. J. Phys. Conf. Ser. 2013, 418, 12079. [Google Scholar] [CrossRef]
  46. Li, G.; Wu, J.; Luo, Z.; Chen, X. Vision-based measurement of dust concentration by image transmission. IEEE Trans. Instrum. Meas. 2019, 68, 3942–3949. [Google Scholar] [CrossRef]
  47. Wang, S.; Li, X.; Wang, S. Separation and fracturing in overlying strata disturbed by longwall mining in a mineral deposit seam. Eng. Geol. 2017, 226, 257–266. [Google Scholar] [CrossRef]
  48. Gonzalez-Vidal, J.J.; Perez-Pueyo, R.; Soneira, M.J. Automatic morphology-based cubic p-spline fitting methodology for smoothing and baseline-removal of Raman spectra. J. Raman Spectrosc. 2017, 48, 878–883. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the image processing algorithms for extracting features of pollutants, and establishing an evaluation mechanism for hazards classification of smoke dust during blasting operations.
Figure 1. Flowchart of the image processing algorithms for extracting features of pollutants, and establishing an evaluation mechanism for hazards classification of smoke dust during blasting operations.
Mathematics 10 03205 g001
Figure 2. Description of the proposed method for the Dbox(P) and Ga of the fugitive dust picture, including (a) original image of fugitive dust, (b) binary image through a series of image-processing operations, (c,d) calculation process of grayscale properties and Dbox(P) and (e,f) results of grayscale distribution, Dbox(P) and Ga of the fugitive dust.
Figure 2. Description of the proposed method for the Dbox(P) and Ga of the fugitive dust picture, including (a) original image of fugitive dust, (b) binary image through a series of image-processing operations, (c,d) calculation process of grayscale properties and Dbox(P) and (e,f) results of grayscale distribution, Dbox(P) and Ga of the fugitive dust.
Mathematics 10 03205 g002
Figure 3. Result comparation on grayscale average between the binary image and grayscale image, including (a) original image of fugitive dust, (b) grayscale image, (c) grayscale average of binary and gray image for algorithm validation.
Figure 3. Result comparation on grayscale average between the binary image and grayscale image, including (a) original image of fugitive dust, (b) grayscale image, (c) grayscale average of binary and gray image for algorithm validation.
Mathematics 10 03205 g003
Figure 4. Frame pictures of pre-blasting, exploding stage and diffusion stages extracted from the UAV video data.
Figure 4. Frame pictures of pre-blasting, exploding stage and diffusion stages extracted from the UAV video data.
Mathematics 10 03205 g004
Figure 5. Grayscale elevation map of smoke dust with homologous calculation processes, including (a) original image of smoke dust, (b) gray image, (c) image fragment of gray image, (d) digital matrix of the image fragment, and (e) grayscale elevation map based on the image fragment.
Figure 5. Grayscale elevation map of smoke dust with homologous calculation processes, including (a) original image of smoke dust, (b) gray image, (c) image fragment of gray image, (d) digital matrix of the image fragment, and (e) grayscale elevation map based on the image fragment.
Mathematics 10 03205 g005
Figure 6. Partial visualization results of the image-processing operations, including original images and binary images, corresponding to the 15th, 30th, 45th, 60th, 75th and 90th images.
Figure 6. Partial visualization results of the image-processing operations, including original images and binary images, corresponding to the 15th, 30th, 45th, 60th, 75th and 90th images.
Mathematics 10 03205 g006
Figure 7. Calculating results of Ga and Dbox(P) with corresponding fitting curves. G1, G2, G3 and D1, D2, D3 represent fitting curves of grayscale average and fractal dimension, respectively.
Figure 7. Calculating results of Ga and Dbox(P) with corresponding fitting curves. G1, G2, G3 and D1, D2, D3 represent fitting curves of grayscale average and fractal dimension, respectively.
Mathematics 10 03205 g007
Figure 8. Calculating results of pollution index with fitting curves of pollutants distribution.
Figure 8. Calculating results of pollution index with fitting curves of pollutants distribution.
Mathematics 10 03205 g008
Figure 9. Four-level divisions of smoke dust pollution within the fitting plot, including a pollution level-pyramid graph representing the pollution condition.
Figure 9. Four-level divisions of smoke dust pollution within the fitting plot, including a pollution level-pyramid graph representing the pollution condition.
Mathematics 10 03205 g009
Figure 10. Calculation results of Dbox(P) and Ga of two single-point pollution areas in the 30th image, including (a) original image selected at 30 s from the recorded video, (b) binary image through a series of image processing operations as described previously, and (c,d) Dbox(P) and Ga of the pollutants in the selected areas.
Figure 10. Calculation results of Dbox(P) and Ga of two single-point pollution areas in the 30th image, including (a) original image selected at 30 s from the recorded video, (b) binary image through a series of image processing operations as described previously, and (c,d) Dbox(P) and Ga of the pollutants in the selected areas.
Mathematics 10 03205 g010
Table 1. Comparation of previous studies on analysis of smoke dust.
Table 1. Comparation of previous studies on analysis of smoke dust.
MethodsMain AdvantagesMain Limitations
Numerical simulation methods [19,20,21]Short research cycle, safety and low investmentLack of authenticity
Machine learning methods [22,23,24]Prediction effectivenessHigh demand for dataset
Measurement equipment for smoke dust [9,28]Visual and clear data with a reference of existing standardlimited measurement distance
Grayscale feature models [32,35]High flexibility; non-contact measurementsFew field applications for smoke dust
Morphological feature models [36,37,38,39,40,41,42,43]
Table 2. Calculation results of fractal dimension and grayscale average for partial pollutant images.
Table 2. Calculation results of fractal dimension and grayscale average for partial pollutant images.
Pollutant ImageParametersCalculating Time of Ga (s)Calculating Time of Dbox(P) (s)
GaDbox(P)
15th000.34911.5598
30th0.33171.72610.35111.1595
45th0.46551.79570.35211.1858
60th0.50231.83760.35021.1795
75th0.47771.82990.36711.1735
90th0.36851.79130.35411.205
Table 3. Fitting results of Ga and Dbox(P) with related R2 for three stages.
Table 3. Fitting results of Ga and Dbox(P) with related R2 for three stages.
Group NumberFitting Equation of Ga and Dbox(P)R2
G1 (Stage one)G1 = 0.382 ln(x) − 0.92390.9165
G2 (Stage two)G2 = −0.0006 x + 0.52120.0583
G3 (Stage three)G3 = −1.94 ln(x) + 9.10070.8188
D1 (Stage one)D1 = 0.1508 ln(x) + 1.25960.5707
D2 (Stage two)D2 = 0.0008 x + 1.7710.2101
D3 (Stage three)D3 = −0.69 ln(x) + 4.9090.4978
Table 4. Comparation of calculation methods.
Table 4. Comparation of calculation methods.
MethodRMSE (PM2.5/PM10/NO2)R Squared (PM2.5/PM10/NO2)
Dbox(P) − Ga based method in this study (Pi1/3)2.85 × 10³/5.99 × 10−40.8365/0.7397
Liu et al. [32]38.28/62.51/25.880.70/0.462/0.328
Zhang et al. [33]27.53/56.67/24.540.881/0.525/0.393
Li et al. [35]50.67/65.19/29.880.563/0.349/0.070
Table 5. Divisions criteria of blasting pollution-levels with the homologous pollution regions.
Table 5. Divisions criteria of blasting pollution-levels with the homologous pollution regions.
Pi<1.9351.935~2.062.06~2.185>2.185
Pollution-LevelLevel ILevel IILevel IIILevel IV
Table 6. Threshold values of T1 and T2 with the time nodes of image acquiring.
Table 6. Threshold values of T1 and T2 with the time nodes of image acquiring.
ThresholdsThreshold Values of T1 and T2 in Binarization
T11–15 s16–35 s36–44 s45–48 s49–95 s
<101008060<40
T21–16 s17–22 s23–31 s32–37 s38–56 s57–63 s64–76 s77–95 s
>240180160150140150160>170
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yin, J.; Lu, J.; Tian, F.; Wang, S. Pollutant Migration Pattern during Open-Pit Rock Blasting Based on Digital Image Analysis Technology. Mathematics 2022, 10, 3205. https://doi.org/10.3390/math10173205

AMA Style

Yin J, Lu J, Tian F, Wang S. Pollutant Migration Pattern during Open-Pit Rock Blasting Based on Digital Image Analysis Technology. Mathematics. 2022; 10(17):3205. https://doi.org/10.3390/math10173205

Chicago/Turabian Style

Yin, Jiangjiang, Jianyou Lu, Fuchao Tian, and Shaofeng Wang. 2022. "Pollutant Migration Pattern during Open-Pit Rock Blasting Based on Digital Image Analysis Technology" Mathematics 10, no. 17: 3205. https://doi.org/10.3390/math10173205

APA Style

Yin, J., Lu, J., Tian, F., & Wang, S. (2022). Pollutant Migration Pattern during Open-Pit Rock Blasting Based on Digital Image Analysis Technology. Mathematics, 10(17), 3205. https://doi.org/10.3390/math10173205

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop