Abstract
In the process of open-pit mining, the system parameters determine the economic benefit and production efficiency of the mine. Conventional optimization involves building a system model for the process parameters. However, complex large-scale systems such as open-pit mining are difficult to model, resulting in a failure to obtain effective solutions. This paper describes a system simulation method for the process parameters involved in open-pit mining. The nature and interaction of each component of the system are analyzed in detail, and the logical flow of each layer of the system is determined. Taking the basic operational linkages of the equipment as the system drivers, we obtained the operational flow of dragline information. The barycentric circular projection method is used to simplify the control logic of the system, and a system storage state model is constructed to identify dynamic changes in the system and obtain the operation parameters of the dragline. A discrete event system is used for quantitative modeling, and the event step method is employed to advance the simulation process and obtain decision information. Finally, simulations are performed using various system parameters. The simulation results show that the maximum efficiency is achieved when the dragline height is ∼13 m, giving a capacity of 4276.52 m3/h. Error analysis indicates that the modeling error is minimized using a simulation correction coefficient of α = 0.94.
1. Introduction
Open-pit mining is the main means of extracting coal resources in the United States, Australia, and many other developed countries. It has been the focus of research on coal mining by researchers worldwide. In recent years, researchers around the world have made a lot of research achievements in the optimization of open-pit mining technology, open-pit mining production planning, and so on. In terms of process optimization of open-pit mining, Ma et al. [1] and Rui [2] proposed 3D visualization mining process optimization technology mainly supported by 3D laser scanning technology, which improved the informatization and intelligence degree of open-pit mining. Song et al. [3] made a systematic optimization and classification study on the green mining process of open-pit coal mine by constructing a reasonable matching relationship model of all links of the process system. Zheng [4] proposed a comprehensive process system optimization method under rigid process constraints in view of the complex geological conditions of open-pit coal mines. In terms of the optimization of open-pit process parameters, Zhao et al. [5] and Ma et al. [6] established the functional relationship between the height of cast blasting steps, the thickness of raw coal, and the production capacity of the dragline and optimized and established the process parameters of cast blasting and dragline stripping technology. Xiao et al. [7, 8] proposed the calculation method of the deformation time and operation cycle of the dragline and established the relationship model between the production efficiency of the dragline and the height and width of the operation platform. Liu [9] proposed three feasible dumping operation modes of dragline and analyzed and evaluated each operation mode according to operation steps, operation parameters, and applicable conditions. For the global optimization of complex systems, significant breakthroughs have been made in the military and medical fields. Pervaiz et al. [10], Bangyal et al. [11, 12], and others have adopted particle swarm optimization, low-difference sequence algorithms, improved bat algorithms (IBA), and other advanced algorithms to improve the detection accuracy of medical diseases and solve the global optimization problems of complex systems. Such methods can be mapped to the optimization of open-pit mine systems.
In terms of mine production planning, Sepulveda et al. [13], Blom et al. [14], and other scholars have developed models that maximize the expected profit of open-pit mines using metaheuristics and simulation tools. For production scheduling, Rezakhah and Newman [15] established the graph structure of the maximum flow in multiple periods and used a genetic algorithm to control the flow within each arc to generate an effective production plan. Open-pit mining is a complex system, and the basic units of the system structure change in time and space according to the mode of the mining action. The characteristics of open-pit stopes are clearly affected by dynamic changes in geological conditions, and open-pit processes involve various modes of cooperation and coupled parameters. Therefore, open-pit mining systems are very difficult to optimize [16, 17].
In view of the above problems, a series of operational optimization methods have been introduced, such as planning theory, graph theory, and queuing theory [18]. The aim of these methods is to approximate various engineering problems so as to construct model forms that can be readily solved, such as linear algebra, linear differential, and linear difference equations. However, there is always a contradiction between the accuracy of the model and the effectiveness of the method, and approximations may lead to solutions that are far from the real solution. Additionally, the real system structure, parameters, and characteristics are mostly nonlinear, described by variable coefficients and variable-structure partial differential equations, for which there is a lack of general and accurate analytical methods. When the system cannot be solved by establishing a mathematical model, it is necessary to describe the evolution of the system and perform numerical simulations, so as to obtain a generalized model that is closer to the actual situation. Quantitative analysis can then be performed through simulation experiments to obtain more accurate solutions [19, 20].
Taking the dragline stripping technology used in Heidaigou open-pit mine, Inner Mongolia, China, as an example, we simulated the efficiency under different flat plate heights and obtained the optimized precise overcast stripping system. The results of this study can be used to improve the operational efficiency and economic benefit of dragline stripping in open-pit mines, and lay the foundation for further application and promotion of dragline stripping technology. This study also provides a new idea for the optimization of open-pit mining processes.
2. Engineering Background
In the Heidaigou open-pit mine, Inner Mongolia, China, the rock (∼40 m thick) above the roof of the coal seam is stripped by blast casting, overcast stripping, and single-bucket truck comprehensive mining. During the dragline operation in the overcast stripping system, a bulldozer and a single-bucket truck drive to the top of the blasted stockpile and discharge the materials from all parts of the working face to the dump in turn (as shown in Figure 1). The designed production capacity of the dragline is 25 million m3/a, but the actual operation volume is less than 17 million m3/a. The overcast stripping system needs to be equipped with an additional single-bucket-truck-assisted stripping volume of nearly 13 million m3/a. The overcast stripping cost of the dragline is about 5 yuan/m3 lower than that of the single-bucket-truck process. Therefore, if the operational efficiency of the dragline equipment could be improved by optimizing the parameters of the dragline working surface, the operational capacity of the equipment would be improved, and the overall operating cost of the system would be greatly reduced.

In Figure 1, the loess area indicates that the loess in the stope of the open-pit mine is stripped by the single-bucket truck. The upper rock area shows the open-pit stope rock being stripped by a single-bucket truck, and the lower rock area shows the stope rock of the open-pit mine being stripped by a single-bucket truck assisted by a dragline. Blast casting refers to the material pile formed after rock step blasting. The dump area of the dragline creates the mound. H1 represents the step height of the single-bucket truck operation area, H2 represents the height of the cast blasting step, and H3 represents the height of the coal seam step.
There are many parameters in the overcast stripping system, among which the height of the flat plate is the most critical. As the flatbed height increases, the overall operating efficiency of the overcast stripping system will be affected as shown in Figure 2.

As can be seen from Figure 2, a change in the standing flat plate height of the dragline shovel will result in both efficiency improvements and efficiency reductions, and the overall efficiency change of the system is difficult to predict. However, the operational logic of the overcast stripping system is complex, and there are a large number of random and fuzzy factors. Thus, it is difficult to use conventional optimization methods to determine the optimal working face parameters.
3. System Simulation Method
3.1. Simulation Optimization Solution Ideas and Steps
In view of the complex internal logical structure of the dragline stripping technology, the interrelation between each process and each level, and the limitations of conventional system modeling, the optimal solution cannot be obtained. Thus, a system optimization method is proposed. The core idea and main steps of the proposed method are as follows:
Step 1. According to the random and fuzzy production characteristics of the dragline stripping technology, a reasonable system simulation model and method are selected.
Step 2. Detailed analysis of the dragline stripping technology elements is carried out to identify all levels, links, and relationships within the system. This allows the system simulation modeling elements to be identified and the events driving the extraction system to be determined.
Step 3. A digital model of the production process of the dragline stripping technology is constructed according to the production factors and their mutual relations with the dragline stripping technology. However, considering the complex internal structure of the dragline stripping technology, the system simulation model needs to be simplified to improve the efficiency of system simulation.
Step 4. In the process of system simulation, all internal elements change dynamically. To accurately simulate the changes in all elements, it is necessary to further build a system state model, quantify the system states under different conditions, and effectively control the system states based on engineering practice.
Step 5. System simulations are implemented, and the system operation process is repeated with various parameter combinations to obtain decision-making information.
3.2. System Simulation Model Types and System Simulation Methods
3.2.1. System Simulation Models
Simulation models can be either continuous or discrete. In the former, the system state changes continuously with time, whereas in the latter, changes to the system state are driven by random events. The ultimate goal of the simulation is to investigate the system efficiency under different working face parameters, that is, to optimize the design scheme so that the same tasks can be completed in the shortest time. Therefore, the time consumption of each basic operation in the system can be regarded as the lowest random event in the simulation process, and this is the construction objective of the quantization model of discrete events. The operation logic of the system constitutes a discrete event system, which is driven by random events and should be described by a discrete event system simulation model [21].
3.2.2. Clock Advance Method
Discrete events are simulated in a process-oriented way and executed in an event-oriented way. Therefore, a timing variable needs to be set to record the occurrence time of a certain low-level random event [22]. In the simulation of a clock advancing during overcast stripping, the time step method calculates the probability of new events occurring at each moment, while the event step method only calculates the probability of current events taking a certain amount of time. Thus, the event step method should be used to advance the analog clock to simplify the calculation scheme.
3.3. System Production Process Analysis and Deconstruction
In the overcast stripping production system, the material migrates from the work area to the pumping area under the external action of the dragline. The overcast stripping working face, the mound, and the dragline equipment constitute the main elements of the system. From the stope structure, the working face can be divided into many parts, and the shape of each part is extremely irregular (as shown in Figure 3), leading to large fluctuations in operational efficiency of different items of equipment. In terms of the equipment operation process, different machine positions may have the dragline at the advanced ditch, discharging coal ditch material to the mound. In addition, there is a strict sequence regarding digging and discharging [23]. Therefore, at every moment of system operation, the geometry of the working face and the mound, the position of the dragline, and the position of material excavation affect the system state. The positions of the dragline and material excavation are specified by operation rules, and can therefore be described by system control logic. The geometry of the working face and mound are subject to passive changes in the system and can only be determined through quantitative calculations using the system state model.

(a)

(b)
Analysis of the system characteristics indicates that there are three levels of logical flow in the process of overcast stripping:(1)Place the dragline in a specific machine position.(2)When the dragline is in a certain machine position, excavate materials in different areas in turn (as shown in Figure 4(a)).(3)In one operation cycle, repeat the dragline maneuvering actions in a certain order (as shown in Figure 4(b)).

(a)

(b)
3.4. Digital Model of System Production Process
In essence, overcast stripping is the operation of moving a large amount of earth and stone. The changes in the material storage space of the working face and pumping area define the whole system process. Suppose that, at any time, the operating state of the system is described by a group of state variables, namely,where X is the state vector and xi is the state variable, i = 1, 2, …, n, during the operation of the system. Different state vectors X form a point in the storage state space Rn. To realize the orderly state transition of the system in the process of simulation, it is necessary to exert more control on the system. Let the system control vector be
According to the state vector and control vector, the triplet [24] of the system state space can be constructed, where is the initial state vector and is the target state vector, describing the storage space state of materials between the working face and the dump before and after the system simulation. C(U) refers to the state transfer operation generated by u and represents the control logic of the system. This expresses the procedural and control knowledge of the system transfer and is determined by the basic operation process.
3.4.1. Barycentric Circular Projection Model of Material Migration
Practical engineering problems are often complex. To enable effective system analysis and synthesis, it is necessary to simplify the problem. In this case, simplification enables the system control logic (i.e., digging position, abandoning position, and overcast stripping track of the dragline) and the system state space (i.e., working face and pumping area of material storage space) to be modeled.
During digging operations, the dragline is usually pushed into the soil and rock for 5–20 m to complete the material loading process. Therefore, the original distribution of the material in each bucket is a strip on the line segment from the operational center point of the dragline along its operational radius. The material on the line segment is generally discarded in the same position; that is, the operational rotation angle of the material on each line segment is the same. Therefore, we establish m center points passing through the dragline operation position that are perpendicular to the horizontal plane . Let Wi represent the overcast stripping by cutting the surface joining Si to Si+1; this is the basic unit of the cutting face and is projected onto a circle with radius R. Similarly, the mound is cut and its spatial capacity is projected onto the circle. When the rays are dense enough, the rotation angles of all materials in the same block can be regarded as the same, thus establishing the gravity circle projection model of the backfill operation (as shown in Figure 5).

The barycentric circular projection model can be described by a directed graph system considering the single-loop transportation of the overcast stripping system (each rotation of the body is no more than 180°). Taking Wi as the starting point and as the end point in Figure 5, D = (V, A) is a simple graph without loops:
In Figure 5, represents the shovel dig point; represents the abandonment point of the dragline shovel; θ1 represents the total angle of the excavator area; θ2 represents the total angle from the excavation area to the dumping area of the dragline shovel; θ3 represents the total angle of the dumping area of the dragline shovel; α represents the angle between the dig sites; and β represents the angle between discard points.
We denote the material volume and internal row capacity as W, and the material capacity of each node as the single transport volume W/m of the equipment. Then, the arc capacity is W/m, and the degree of the node . The angular interval of material nodes in the mining area is α, and the central angle of the fan-shaped area formed on the center of the gravity circle is θ1. The angular interval of the nodes in the abandoned region is β, and the central angle of the region is θ3. The central angle of excavation and discarding in the middle area is θ2, and so the total rotation angle Ψ of single-ring operation is
According to (4), when the size and location of the material in the working face and the space of the dump are determined, the total rotation angle Ψ of the material is a constant value.
3.4.2. Optimization of System Control Logic
The relationship between the time consuming and the rotation angle of the material transportation link (full load rotation, no load return) is nonlinear with respect to time. According to a large number of statistical results, is a concave function, so it is necessary to discuss the influence of different material transfer sequences and paths on the total time . Let , , . Assume that , , and . According to the mean value theorem, there is some , such that
As is a strictly concave function, decreases monotonically on [a, b], i.e., :
With being a fixed value, when increases, decreases. Thus, in general, we have
That is, “gap discarding” (i.e., the distribution of a single rotation angle tends to the maximum and minimum values of the defined domain) gives the lowest total transport time. Considering practical problems, the order of excavation and discarding at each node can be obtained from
Therefore, the materials at point Wi in Figure 5 will be discarded at point .
3.5. System State Model
3.5.1. 3D Simulation Model of Basic Operation Unit
In the simulations, the variations in the material storage space of the working face and pumping area must be quantitatively calculated to identify the dynamic changes between the two and build the system state model. However, the overcast stripping of the basic unit in the initial and final simulation states is not regular, and the conventional planar graphic method cannot calculate the material space and volume. Therefore, it is necessary to build a 3D simulation model of the basic unit of the overcast stripping system (as shown in Figure 6) and provide a quantitative basis for the construction of the system state model. The model building process is outlined below:(1)The size of the upper top surface of the advance ditch unit can be directly obtained by a planar drawing. The slope angles of high step and the step angles of the dragline operation are 65° and 38°, respectively. The arc surface formed by bucket cutting is simplified into an inclined plane with a natural repose angle of 38°, so as to construct the advance ditch unit.(2)The geometric shape of the lower part of the step element body is relatively complex. In this study, Boolean calculations are applied to the established units of the advance ditch, coal seam, and mound, and a 3D model of the step element body is formed after cutting.(3)The unit body of the mound is formed by the sequential accumulation of materials at multiple operation sites. Iterative modeling is applied until the volume and position of the mound are stable and have a certain periodicity, giving the real form of the actual mound.(4)According to the center of gravity circular projection model, block cutting and projection are applied to pull the dragline operation rotation center to the center of the circle in the constructed 3D model. The included angle of each section is 5°. Thus, the spatial positions of the volume and center of gravity of the unit block of the working face after cutting are obtained (as shown in Figure 6(b)). The corresponding abandonment position of each cut block is then determined using the optimized excavation–abandonment rule.

(a)

(b)
Note that, in the process of 3D model construction, key issues such as the size of the basic unit, upper and lower layer division, Boolean operation principle, and cutting position are closely related to the equipment operation characteristics and operation mode. For detailed modeling methods and principles, please refer to [25–28].
3.5.2. System State Model
The system state model is constructed according to the barycentric circular projection model:where and are state variables, represents the current remaining material amount of block after cutting the basic unit of the working face, and represents the current material storage of the corresponding abandoned position of block .
Obviously, the initial state vector of the system isand the system target state vector is
The basic unit of the system state change is the material in a single bucket. Thus, we set an operation cycle to dump a bucket of material from block of the working face to abandonment point . The control vector of the system iswhereby the system state changes to
The system state vector is transformed from to , and the number of transitions is
The average time required for a state change (the average time of the underlying driver event loop) is
The operating efficiency of the equipment iswhere t is the system operation time, which is recorded by the timing variable.
3.6. System Simulation Model Construction and Implementation
For discrete events, the average index of statistical data does not always reflect the regularity of actual production conditions. Therefore, it is necessary to build a probability distribution model for the lowest level of random events and then reproduce the temporal regularity of random events through a random number generator to restore the actual operation conditions as realistically as possible. On this basis, the timing variable plus the corresponding event time are used to determine the time of the next event in the simulation process.
Analysis of the system characteristics indicates that the time required for digging and bucket adjustment is related to the digging position, while the time required for bucket rotation is related to the digging position and the rotation angle. Therefore, for the digging and bucket adjustment links, the probability distribution function is , whereas for the rotation links, the probability distribution function is . Once the system simulation model has been constructed, the system simulation enters the implementation stage.
Figure 7 shows a flowchart of the simulation calculations, where is the rotation angle of equipment operation; represents the different excavation positions (inside of advance ditch, outside of advance ditch, upper layer inside of step, lower layer inside of step, upper layer outside of step, lower layer outside of step, coal ditch); and represents different operation procedures (digging, full load rotation, discharge, no-load return, bucket adjustment).

Obviously, the dragline stripping technology is relatively complex, making it very difficult to solve the model directly. Analysis of this complex process and the characteristics and relationships of each link has led to the establishment of a quantitative model of the process. In this model, each link behaves according to experimental or quantitative analysis, and specific and detailed information is required to obtain suitable decisions.
4. System Simulation Results
Long-term data collection indicates that the temporal distribution of each facet of dragline mining is as shown in Figures 8–12. The probability distribution functions of the operation time were established on this basis [29, 30], and the results are presented in Table 1 for the example of a full load rotation link. During the simulations, the different functions in Table 1 are called according to the system logic.





We simulated the operation of the dragline in one cycle under working flat plate heights of 11–15 m and a full bucket coefficient of 0.75. The calculation results are presented in Tables 2 and 3.
The simulation results show the following:(1)When the step height rises from 11 m to 15 m, the operation rotation angle decreases at first and then increases. However, the rotation angle has little influence on the average time of the rotation link.(2)The adjustment time is the main factor affecting the average time of an operation cycle. This is because, as the step height changes, the initial state of the material distribution in the working face changes. Additionally, with a 15 m bench height, the materials in the coal groove reach a total of 29969 m3, nearly double that with an 11 m bench height. The average time required to adjust the bucket in these two areas is 7–8 s, whereas the average adjustment time in other areas is only ∼4 s.(3)When the dragline height of the scraper is around 13 m, the operating efficiency of the equipment reaches a maximum of 4276.52 m3/h. According to field data from Heidaigou open-pit mine, the production capacity of the scraper reaches a maximum of 4026 m3/h when the flatbed height is 13.3 m [31]. The simulation results for the operating efficiency are therefore consistent with the actual situation.
5. Discussion and Analysis of Results
For the overcast stripping system, the flatbed height is an extremely important parameter in the dragline operation. Therefore, this study simulated the one-cycle operation of the dragline under working flat plate heights of 11–15 m. The simulation results in Tables 2 and 3, combined with the production characteristics of the overcast stripping system, can be analyzed as follows:(1)When the number of bulldozers and the scale of the single-bucket-truck auxiliary dragline operation are reduced, the layered advance speed on the cast blasting stockpile accelerates. Simultaneously, the step advance speed of the dragline operation decreases, the coal seam advance speed is reduced, the overall stripping cost of the system is reduced, and the raw coal output is limited.(2)As the width of the operation steps decreases, the basic unit shape of the dragline operation changes, and the scale of the secondary dragline operation decreases. Thus, the amount of material in the coal ditch decreases, and the total system output decreases.(3)When the center line of the dragline operation is far away from the pumping area, the average rotation angle of equipment operation increases, and the cycle time of one operation increases.(4)When the dragline cycle advance distance decreases, the cycle advance time increases, and the time required to move equipment and move the shovel increases.
Error analysis of the simulation results is presented in Table 4. The conclusions are as follows:(1)The simulation system is an abstract and idealized description of the actual system. However, there are a large number of uncontrollable factors in the actual system, and it is difficult to operate efficiently according to the standardized system logic process. Therefore, the simulation results are more efficient than the actual situation.(2)During the actual operation of the overcast stripping system, the duration of each system state X is discretely distributed, and its control vector U is also discretely distributed. For example, in the process of equipment operation, the full bucket rate of each operation cycle should be a discrete event. However, in this study, the constant value of the full bucket coefficient was set to 0.75 based on experience. This will result in some errors.(3)To improve the authenticity of the simulation process and the accuracy of the simulation solutions, the error sources can be considered comprehensively, and the simulation results can be corrected. The error can be greatly reduced when the correction coefficient is set to 0.94.
6. Conclusion
This paper has described an optimization method for the parameters of an open-pit mining system based on simulation technology. Quantitative analysis was performed by repeating the evolution of the system to obtain more accurate solutions.
The core ideas and simulation steps were established by focusing on dragline stripping technology. The components of the overcast stripping system were deconstructed, and the process flow of the system was determined. The events driving the system were also identified. A simplified simulation model and control logic were then constructed, and the dynamic laws of the simulation system were obtained. Finally, simulations of the system operation process were repeated to obtain decision-making information.
The Heidaigou open-pit mine overcast stripping system was selected as a specific case study. The logical process of the system was decomposed into three levels, and the basic operation links of the dragline were taken as the driving events of system. A barycentric circular projection method was proposed to simplify the control logic of the system. To accurately determine the dynamics of each element, a system storage state model was built to quantify the system state under different conditions. Finally, the discrete event system was modeled, whereby the simulation process was advanced by the event step method.
For the overcast stripping system, the flatbed height is a key parameter in the operation of the dragline. Therefore, this study simulated the operation of the dragline in one cycle under working flat plate heights of 11–15 m. The simulation results show that the operating efficiency of the equipment is maximized when the working flat plate height of the dragline is around 13 m. Error analysis of the simulation results found that a correction coefficient of α = 0.94 significantly reduces the modeling error.
Data Availability
The basic data involved in this study were obtained by the author and the production site of Heidaigou open-pit coal mine for a long time.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
The work was supported by the basic scientific research service fees of central universities (3142019007).