A new circulating accumulation emission model for assessing dust emission from open pit mine

To reduce the inaccuracy of using the monitoring data outside the pit to evaluate the unorganized emission dust source of open pit mine, the circulating accumulation emission model is established. Based on the model, the monitoring data in the pit can be converted into the dust emission from the pit. The main conclusions include: (1) the circulating accumulation emission model is suitable for the dust diffusion process in open pit mine. The ratio of diffusion \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μ and the ratio of surplus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon$$\end{document}ε were used to simulate the dust diffusion process in open pit mine, containing emission, retention and diffusion. (2) The initial value of the dust in the pit before the team operation has little influence on the final stable value. (3) When the external dust enters the pit, it will accumulate under the action of eddy current. The dust background value in the pit is different from that outside the pit. (4) The dust emission from the pit can be calculated from the monitoring data in the pit based on the circulating accumulation emission model. The model can deal with environmental changes such as the wind direction and speed, without arranging a lot of external monitoring equipment like the traditional external monitoring methods.

Surface mine is an unorganized dust source. At present, the monitoring method of unorganized dust emission intensity is to set up monitoring equipment separately in the upper and lower wind directions. The difference between the two types of monitoring data is the emission intensity of the unorganized dust source. However, due to the large area and varied shape of the open pit and a large number of dust sources in the pit, this value usually can not correctly reflect the dust pollution in the open pit in many cases. For example, when the wind direction or wind force changes, the upwind direction and downwind direction also change. When the shape of the mine pit and wind direction are in different combinations, the height of dust emitted outside the mine pit changes in a large range, etc. For this type of unorganized dust source of open pit, if we want to get the dust emission of the pit through the monitoring outside the pit, we need to arrange a large number of monitoring points in different positions and different elevations outside the pit, which is difficult to complete in practice.
According to the previous research results 1 , there is no significant difference in dust concentration in most spaces away from dust sources in the pit. The consistent law of dust distribution in the pit can be used to study the change rule of dust concentration in the mine pit, establish the dust diffusion model, and calculate the dust emission of the open pit mine.
As the monitoring data represents the accumulated value of dust emission, it may be unstable even though the emission keeps stable, because of dust cumulative process is affected by turbulence, which is easily affected by geometric dimension of open pit mine and weather conditions, such as wind, temperature, humidity et al. 2 . So the monitoring data cannot be directly used for assessing the dust emission intensity of mining activities.
In this study, the circulating accumulation emission model is established to convert the monitoring data in the pit into the dust emission of the pit. This method is less affected by wind direction, mine pit size and other factors. Hence the number of monitoring equipment needed is less. The dust emission assessment of open pit can be completed more efficiently.

Module of diffusion process.
To model the diffusion process, a flow chart is used by decomposing the process, as shown in Fig. 2. This flow chart is the abstraction of the diffusion process shown in Fig. 1 4 . The meaning of them is similar. But based on the abstraction chart, the logical relationship among the parts of the diffusion process is clearer 5 . This relationship is consistent with the Markov birth and death process 6,7 . This model analyzes the law of dust accumulation and diffusion in the mine pit under the influence of eddy current and its calculation model, so it is suitable for small particle size dust and not for large particle size dust. Therefore, it is not applicable to all PMs. It has good applicability to respirable dust, such as PM2.5, PM10, etc.

The dust emission and diffusion model for open pit mine
The basic form of model. To describe the diffusion process, time series analysis is suitable 8    where ε = ratio of surplus in the pit, 0 < ε < 1.

The variable Q(t) is a time series of concentration of PMs. Derive the above equation by the variable t:
Further, recurse the equation by the variable t: Finally, this equation is got by deriving the variable t to 1: A simplified form of model by dividing the operation time. It is unrealistic to obtain the real value of dust emission from mines with a general area of several square kilometers in real-time. So the model should be properly simplified so that it can not only calculate the dust emission but also ensure a certain accuracy. The basic idea of the simplified model in this paper is to segment the operation time according to the changes of influencing factors. In these divided time periods, since the influencing factors of the variable have not fluctuated greatly, it can naturally be considered that the variable has not fluctuated greatly. Since the main influencing factor of dust emission is the intensity and type of mining operation, the operation time can be segmented according to the actual operation intensity of the mine. In the divided operation time period, the electric shovel peels the coal and rock mass at the same position, the drilling rig drills at the same position, and the truck runs on the same road. Therefore, it can be assumed that the dust emission from mining operations in the pit has not changed during this time period. If the mining operation conditions change, the time period shall also be divided accordingly.
In the divided time period, the emission increment is assumed to be constant, which means that: The model is simplified: Arrange the equation: The above equation is the simplified form of model if the increase is constant. As 0 < ε < 1 , the limit of Q(t) is got: The above equation shows that the model is convergent. The variable lim t→∞ Q(t) represents the final steady state of the simplified model. Background value. The background value is the concentration of PMs existing in the surrounding atmosphere. In the national standard, the background value should be monitored and subtracted for researching the pollution emission source 9 . However, due to the turbulence in the pit, not only the dust generated in the pit will accu- www.nature.com/scientificreports/ mulate, but also the dust in the surrounding environment outside the pit will accumulate after entering the pit. The dust emission is calculated by the method provided in this paper, and the detection data in the pit are used, which is the sum of the accumulated value of dust in the pit and the accumulated value of surrounding dust after entering the pit. Since the dust in the surrounding environment will accumulate and increase after entering the pit, its embodiment on the monitoring instrument in the pit is not equal to the concentration detected outside the pit. To eliminate the effect of background value, the monitoring data should be subtracted by the accumulated value, not the value monitored directly out of pit 10 . Although the dust background concentration outside the pit is changing, it often remains stable for a period of time. Whether it fluctuates greatly can be judged by the consistency test. If the test is passed, the average value in this period shall be used instead. If the inspection fails, the monitoring data is segmented. Recheck whether the data in each segment meets the consistency test. Repeat the above operations until all segments of data meet the consistency test. In the divided time period, the dust background concentration q back is assumed to be constant.
By using the simplified model, when the variable Q(0) is 0, the effect of background value in the pit is clear: where q back (t) = the background value out of pit, Q back (t) = the effect of q back (t) in the pit. The limit of Q back (t) is similar to limit of Q(t).
If the fluctuation range of background value can be seen as constant, the limit of Q back (t) would be a better way to calculate the effect of background value in the pit.
If a mutation appears in the background value, the simplified model cannot work, the basic form of dust diffusion model should be used, and the effect of the mutation is calculated based on the Eqs. 3-5. So, the concentration of PMs monitored in the pit should be handled by this equation: where Q monitor (t) = the concentration of PMs monitored in the pit. The variable Q(t) , calculated by Eq. (12), can be dealt with by the diffusion model for understanding the rules and characteristics in an open pit mine after monitoring data obtained.
Ratio of surplus in the pit. The dust emission of open-pit mine is mainly caused by mining, blasting, transportation and dumping 11 . Although there is dust caused by natural wind, it accounts for a small proportion according to the research of previous scholars 12,13 . The premise of no dust emission proposed here is not that the wind speed is very small, but by analyzing the daily mining operation process of open-pit mine, select the time period when most equipment does not operate. During these non-operation periods, the manual operation is mostly stopped, and the dust generated by manual operation is ignored. During this period, the dust in the pit is gradually dissipated. The external wind speed is the same as the operation time period, that is, the diffusion conditions have not changed. Therefore, the dust monitoring data in the non-operation period can be used to calculate the diffusion coefficient, and the diffusion coefficient is also effective in the operation period with constant diffusion conditions. According to the law of mining operation, these non-operation periods occur many times every day. Such as daily maintenance period, daily evacuation period before blasting operation, shift handover period of different operation shifts, etc. The specific operation is to select one or more non-operation time periods according to the mining operation time arrangement, record the monitoring data and calculate the diffusion coefficient, which respectively represents the diffusion coefficient for a period of time. The escape rate calculated for these time periods acts on half of the adjacent time before and after them.
When q(1) = q(2) = · · · = (t − 1) = q(t) = 0 , which represents that there is no emission, the model is shown as below: Arrange the equation: Then, a very important variable in this model is calculated out.
Value of emission. The variable q(t) represents all emissions in the pit on the time t . Under the simplified model, the Eq. (9) can be used for calculating it. Arrange the Eq. (9): where q year = Annual emissions, V = The air volume of open pit mine, T = Annual working hours.

Conclusions
Based on analyzing the accumulation and emission of dust concentration by dust monitoring data in an open pit, the circulating accumulation emission model is established. The advantage of the model constructed in this paper is to deduce the overall dust emission of the open-pit mine by analyzing the dust concentration change of the flow field in the pit.The model can deal with the change of environment, such as the change of wind speed and wind direction, without arranging a lot of external monitoring equipment like the traditional external monitoring methods. Because the traditional unorganized dust source monitoring method is used, a whole circle of dust monitoring equipment needs to be arranged around the mine to meet the monitoring conditions of different wind directions. 1. The model is suitable for the dust diffusion process in open pit mine, containing emission, retention, and diffusion. 2. The initial value of the dust in the pit before the team operation has little influence on the final stable value. 3. The dust background value of the surrounding environment monitored outside the pit can not be directly used. Because the dust of the surrounding environment may accumulate in the pit with turbulence. To eliminate the effect of background value, the monitoring data should be subtracted by the accumulated value, not the value outside the pit monitored directly. 4. Formula 16 can be used to calculate the dust emission in a certain monitoring period. If annual monitoring is maintained, the annual dust emission from an open pit can be obtained.