# Research on Optimization of Mining Substitution in Wangzhuang Mine Based on System Dynamics

Authors: Shuheng Zhong; Yuchen Guo; Yanjiao Li; Ming Yang
Journal: Advances in Civil Engineering
 (2021)

---

## Abstract

The unbalanced mining replacement is a major problem restricting efficient mining. The optimization of mining replacement process can achieve perfect coordination of mining face and tunneling in time and ensure the efficient production of mine to the greatest extent. Based on this, this paper takes Wangzhuang Mine as the research background, applies system dynamics to mining replacement research, builds the Wangzhuang Mine mining-driving system simulation model, and uses this model to dynamically simulate the working face 9102 and its replaced heading face 7106. The research has found that there is an imbalance in the replacement between the two working faces, and the replacement process can be optimized by advancing the driving period or increasing the driving team to work in parallel at the appropriate time point, so as to meet the normal demand of mining replacement and provide guidance for the mining plan arrangement of other similar working faces.


---

## Body

## 1. Introduction

With the development of big data technology and mechanization and intelligence, the advancing speed of the working face has become faster and faster. Although the advancing speed of heading face has been improved, it still lags behind the working face, which causes the imbalance of mining replacement and affects the production efficiency of the mine. Therefore, in order to improve the coal output and to ensure the high efficiency of mine production, it is necessary to optimize the mining process [1–3]. Aiming at the problem that the excavation lags behind the mining, the main factors restricting the speed of roadway driving were summarized and the future development direction was proposed [4–6]. Some scholars have put forward prospects and suggestions for the improvement of driving speed from the aspects of tunneling technology and equipment [7–9]. Some scholars have established an intelligent decision system for mining replacement with the help of computers [10–12].References [4–6] only put forward the problem of comprehensive roadway excavation in our country, indicate the future development direction of tunneling, or just point out the overall process cycle optimization. References [7–9] explain the current situation of rapid tunneling technology, discuss key technical issues, or quickly optimize tunneling technology through theoretical calculation. References [10–12] abandon the traditional manual planning and gradually use the cross-platform and versatile software in the computer to prepare the mining succession planning system. The above research is committed to explore how to improve the speed of excavation to cooperate with mining face completing the excavation replacement both from the perspective of technology and intelligent decision and alleviate the imbalance problem of excavation replacement to a certain extent. However, it is rarely considered from the perspective of mining replacement optimization. The optimization of mining and driving replacement is to arrange the driving tasks reasonably in terms of time to meet the replacement demand and avoid the waste of resources, which is the prerequisite for the realization of intelligent mining. Based on this, this paper applies system dynamics to the study of optimization of mining and driving replacement to improve mine production efficiency better.


## 2. Model Building

System dynamics is a discipline which could study the dynamic development process of complex systems and uses computer simulation technology to simulate the dynamic development process of the system based on feedback control theory [13–15]. The modeling process mainly includes two aspects, model building (system boundary determination, causality analysis, and structural flow diagram drawing) and simulation (parameter setting, simulation analysis, countermeasures, and suggestions) [16, 17].


### 2.1. System Boundary Determination

There are two mining levels of +630 and +540 in Wangzhuang Mine. The mining coal seam is No. 3 coal seam. Working plane which has been replaced recently includes 9102 mining face, 7106 ventilating roadway, and 7106 transportation roadway. The 9102 mining face is located in the 91 mining area. In the working face, there are 3 faults and 2 collapse columns (fault F304H = 5.0 m ∠65°, fault F317 H = 8.5 m ∠40°, fault F318 H = 1.4 m ∠40°, collapse column H51, and collapse column X14). The 7106 working face passes through five anomaly areas. They are the comprehensive response of coal seam thinning, crack development, concealed geological structure, or coal structure abnormality. From the perspective of current situation of the mine production, mining replacement is unbalanced. Through the analysis, the factors affecting the driving speed of heading face 7106 and the advancing speed of working face 9102 are summarized into the following six aspects, as shown in Table 1. In the table, “√” means that it will affect the system, and “×” means that it will not affect the system.Table 1
Key influencing factors of mining system.
The key factorsDetailed indexDriving systemCoal mining systemGeological structure levelDistance of geological structural zone (m)√√Mine pressure level of underground roadwayAdvance support distance of roadway (m)×√Two-sided displacement (mm)×√Working environment levelWater inflow per unit time at mining face (m3/h)√√Roof crushing degree√√Roof falling capacity (mm)√√Roof subsidence rate (mm/min)√√Return air flow gas concentration at the mining face (%)√√Regular circulation levelNormal circulation rate (%)√√Mechanical and electrical equipment failure levelElectromechanical equipment accident rate (%)√√Power on levelOperating rate (%)√√


### 2.2. Causality Analysis

Causality analysis can effectively analyze the relationship between various factors in the system [18, 19]. The causality diagram of the excavating system is based on the coal mining system and the driving system, illustrating the influence of various factors on the rapid driving of the roadway, the advancing speed of the mining face, and the interaction between the two systems. The factors that affect the excavating system mainly include geological structure, operating environment, mechanical equipment, regular circulation level, etc. The influence of each factor on the system is shown in Figure 1. The causal relationship between each factor is represented by a causal chain composed of causal arrows with positive (+) and negative (−) signs. Taking “Airflow Gas Concentration of Heading Face” as an example, the Airflow Gas Concentration of Driving Face ↑⟶ (−) Gas Level of Driving Face ↓⟶ (−) Operating Environment Level of Driving Face ↓⟶ (−) Tunneling footage ↓⟶ (−) Driving distance of Driving Face↓.Figure 1
Causality diagram of mining-driving system.


### 2.3. System Structure Flow Diagram

Based on the analysis of the causality of the excavating system, the structure flow diagram of the excavating system is drawn with Vensim software, as shown in Figure2. The structure flow diagram contains two rate variables: the driving speed and the advancing speed of the working face, and two state variables (cumulative quantities): the driving length and the advancing distance of the working face.Figure 2
Structure flow diagram of mining-driving system.


## 2.1. System Boundary Determination

There are two mining levels of +630 and +540 in Wangzhuang Mine. The mining coal seam is No. 3 coal seam. Working plane which has been replaced recently includes 9102 mining face, 7106 ventilating roadway, and 7106 transportation roadway. The 9102 mining face is located in the 91 mining area. In the working face, there are 3 faults and 2 collapse columns (fault F304H = 5.0 m ∠65°, fault F317 H = 8.5 m ∠40°, fault F318 H = 1.4 m ∠40°, collapse column H51, and collapse column X14). The 7106 working face passes through five anomaly areas. They are the comprehensive response of coal seam thinning, crack development, concealed geological structure, or coal structure abnormality. From the perspective of current situation of the mine production, mining replacement is unbalanced. Through the analysis, the factors affecting the driving speed of heading face 7106 and the advancing speed of working face 9102 are summarized into the following six aspects, as shown in Table 1. In the table, “√” means that it will affect the system, and “×” means that it will not affect the system.Table 1
Key influencing factors of mining system.
The key factorsDetailed indexDriving systemCoal mining systemGeological structure levelDistance of geological structural zone (m)√√Mine pressure level of underground roadwayAdvance support distance of roadway (m)×√Two-sided displacement (mm)×√Working environment levelWater inflow per unit time at mining face (m3/h)√√Roof crushing degree√√Roof falling capacity (mm)√√Roof subsidence rate (mm/min)√√Return air flow gas concentration at the mining face (%)√√Regular circulation levelNormal circulation rate (%)√√Mechanical and electrical equipment failure levelElectromechanical equipment accident rate (%)√√Power on levelOperating rate (%)√√


## 2.2. Causality Analysis

Causality analysis can effectively analyze the relationship between various factors in the system [18, 19]. The causality diagram of the excavating system is based on the coal mining system and the driving system, illustrating the influence of various factors on the rapid driving of the roadway, the advancing speed of the mining face, and the interaction between the two systems. The factors that affect the excavating system mainly include geological structure, operating environment, mechanical equipment, regular circulation level, etc. The influence of each factor on the system is shown in Figure 1. The causal relationship between each factor is represented by a causal chain composed of causal arrows with positive (+) and negative (−) signs. Taking “Airflow Gas Concentration of Heading Face” as an example, the Airflow Gas Concentration of Driving Face ↑⟶ (−) Gas Level of Driving Face ↓⟶ (−) Operating Environment Level of Driving Face ↓⟶ (−) Tunneling footage ↓⟶ (−) Driving distance of Driving Face↓.Figure 1
Causality diagram of mining-driving system.


## 2.3. System Structure Flow Diagram

Based on the analysis of the causality of the excavating system, the structure flow diagram of the excavating system is drawn with Vensim software, as shown in Figure2. The structure flow diagram contains two rate variables: the driving speed and the advancing speed of the working face, and two state variables (cumulative quantities): the driving length and the advancing distance of the working face.Figure 2
Structure flow diagram of mining-driving system.


## 3. Simulation

The above established model is applied to the 9102 working face and 7106 heading face of Wangzhuang Mine. The 9102 working face has a total length of 3,500 meters, the actual mining distance is 2450 meters, and the mining time is 24 months. The total length of the 7106 ventilation roadway and transport roadway lane is 2500 meters, and the driving time is 27 months. Therefore, the system simulation time is set as 27 months and simulation step length is 1 month.


### 3.1. Parameter Setting

#### 3.1.1. Weight Determination

The weight values of each index are determined by the group AHP method except the two indexes of advanced roadway support level and the moving level of two sides of roadway, which were given by experts. The weight results are shown in Table2.Table 2
Simulation model weight assignment.
IndexWjWeight valueInfluence of geological structure level of heading face on heading speedW10.3176Influence of working environment level of heading face on heading speedW20.3175Influence of normal level of heading face on heading speedW30.123Influence of starting level of heading face on heading speedW40.123Influence of mechanical equipment fault level coal mining on heading face speed in heading faceW50.1189Mine pressure level of heading faceW60.3008Gas level of heading faceW70.4901Mine water level of heading faceW80.2091Roof breakage degree of heading roadwayW90.4012Roof subsidence level of heading faceW100.2908Approaching level of two sides of tunneling roadwayW110.308Influence of geological structure level on advancing speed of coal mining faceW120.2986Influence of working environment level of coal mining face on advancing speed of working faceW130.1962Influence of the normal level of the working face on the advancing speed of the working faceW140.0616Influence of mining face level on working face advancing speedW150.0792Influence of mechanical equipment failure level of coal mining on working face advancing speedW160.0835Influence of roadway end management level on working face advancing speedW170.099Influence of driving speed on advancing speed of working faceW180.1819Mine pressure level of coal mining faceW190.3325Gas level in coal mining faceW200.5278Mine water level of coal mining faceW210.1397Roof breakage degree of coal mining faceW220.3062Roof subsidence level of coal mining faceW230.3012Roof deformation level of coal mining faceW240.3926Roadway advanced support levelW250.35The approaching level of two sides of roadwayW260.65


#### 3.1.2. Basic Parameter Setting

The basic parameters mainly include the distance from the tunneling face to the geological structure zone, the failure rate of the tunneling face mechanical equipment, the gas concentration of the return air flow from the heading face, the distance from the mining face to the geological structure zone, the normal circulation rate of the mining face, etc.


#### 3.1.3. Assignment of Auxiliary Variables

The level of each influencing factor changes as the value of the factor changes. This paper uses table function to express the relationship between the level and the actual value of the influencing factors. The level of influencing factors was quantified with 40∼100 scores. The influence degree is represented by the promotion or weakening effect on the numerical value of the affected factor, whose value may be greater than or less than 1. Considering the research results of many scholars, the paper defines the range of influence to be between 0.5 and 1.5 [20].


### 3.2. Simulation Analysis

#### 3.2.1. Model Checking

Vensim software has its own “check model” function, which can check the consistency of model units. The simulation model in this paper runs well, and the consistency of the model units is verified by the “check model” function. After running the model, it can be obtained that the simulated output value of driving length is 2,441.765 meters, and the error rate is −2.33% compared with the actual value of 2500 meters. Compared with the actual value of 2450 meters, the simulated output value of the advance distance of the coal mining face is 2526.87 meters, and the error rate is +3.14%. The absolute error value of both is less than 5%, indicating that the model has a good effect and is qualified.


#### 3.2.2. Analysis of Simulation Results

(1)
The relationship between the level of each influencing factor and the driving speed is shown in Figures3(a) and 3(b). It can be seen from the figure that all influencing factors have an impact on the driving speed. The level of geological structure reached its lowest point in the 5th and 19th months, and the level of other influencing factors reached their lowest in the 6th and 20th months. Driving speed also reached its lowest point in the 6th and 20th months correspondingly. The curve trend of the driving speed was roughly the same as the horizontal curve trend of each influencing factor.(2)
The relationship between the level of each influencing factor and the advancing speed of the working face is shown in Figures4(a) and 4(b). It can be seen from the figure that various factors have an impact on the advancing speed of the working face. The level of geological structure remained at a relatively low level from 4th to 9th and 16th to 18th months, while other factors also fluctuated at a low level during this period. The driving speed reached its lowest point in the 6th and 17th months, and the trend of the tunneling speed curve of the working face was roughly consistent with that of the horizontal curve of the influencing factors.(3)
The simulation results of the mining-driving system are shown in Figures5(a)–5(d). From Figure 5(a), it can be seen that the advancing speed of the 9102 working face has been significantly slower since April. Also, for nearly 5 months thereafter, its speed remained at a low level, rose in the 11th month, and then significantly decreased in the 17th month. But the duration is short. It returned to above the normal advancing speed in the 20th month. The main reason is that the 9102 working face was getting closer and closer to the F304 fault in the fourth month, and the span of F304 is 490 meters. Therefore, under the influence of the fault, the driving speed remained at a low state for a long time. In the 17th month, the working face encountered collapse columns H51 and X14 one after another. The two collapse columns were adjacent and the span was relatively small. Therefore, the impact time on the speed was relatively short. Correspondingly, the working face advance distance curve slowed down as its speed decreased. It can be seen from Figure 5(b) that the driving speed of the heading face decreased significantly in the 4th month, reached the lowest point in the 6th month, then decreased significantly in the 16th month, and reached the lowest point in the 20th month. The lowest speed was lower than that in the 6th month. The main reason is that the heading face is getting closer to anomaly area 3 in 4 and 16 months, and the return air roadway traverses the anomaly area, so it has a greater influence on it. As a result, the driving speed is lower and the duration is longer in the 20th month, thus affecting the normal replacement between roadway driving. The driving footage is the cumulative amount of the driving speed. According to the driving footage curve, it can also be seen that the trend of the driving footage curve becomes slow in the month with a lower driving speed. It can be seen from Figure 5(c) that under normal circumstances, the driving speed is lower than the advancing speed of the working face, and it has a certain impact on the advancing speed of the working face. The trends are not exactly the same because the advancing speed of working face is also affected by other factors. Figure 5(d) shows that the fully mechanized mining team completed the entire construction period in the 24th month, while the tunneling team completed the entire construction period in the 27th month, which is 3 months later than the mining face, and there is an imbalance problem between mining and driving.Figure 3
Relationship between driving speed and influencing factors.
(a)(b)Figure 4
Relationship between the level of each influencing factor and the advancing speed of coal mining face.
(a)(b)Figure 5
Simulation results of mining-driving system: (a) relationship between advancing speed and advancing distance of working face; (b) relationship between driving speed and driving footage; (c) relationship between driving footage and advancing distance of working face; (d) relationship between driving speed and advancing speed of working face.
(a)(b)(c)(d)


### 3.3. The Optimization of Mining and Driving Replacement

It can be seen from the simulation results that there is an imbalance problem in the replacement of mining and driving. Due to the relatively complicated process of the tunneling face itself, the normal driving speed is relatively slow. Secondly, the tunneling face encountered anomaly area 3 in the 4th and 16th months, which affected the normal driving speed. For the reason that most of the factors affecting the driving speed are the changes of geological conditions, resulting in the change of other factors, and most of them belong to the mine hydrogeological factors, it is difficult to improve the driving speed by changing these factors. In order to ensure the balance of mining replacement, the following two aspects can be considered.(1)
Advance tunneling period: from Figure5(d), it can be seen that the coal mining team completes the construction period 3 months earlier than the tunneling team. If the tunneling team starts working three months earlier, the construction period will be completed at the same time as the mining team to ensure the normal replacement of mining.(2)
Arranging more tunneling teams at the right time: when the resources are sufficient, two driving teams can be reasonably arranged to work. The tunneling working face is composed of transport lane and return air lane, and the two roadways do not affect each other when working, which can meet the requirements of two tunneling teams working at the same time. Through further simulation, the paper found that in the 11th month, a new driving team was added to start work in the transport lane, and the two teams worked in parallel from 11 to 13 months, and after the 13th month (return air lane driving task is completed), the remaining roadway driving tasks are completed independently by the second team. The specific simulation results are shown in Figure6.Figure 6
Simulation results of mining plan arrangement: (a) relationship between original driving speed and current driving speed; (b) relationship between current driving speed and driving footage; (c) relationship between original driving footage and current driving footage; (d) relationship between working face advancing distance and current driving footage.
(a)(b)(c)(d)It can be seen from Figure6(a) that in the 11th month, due to the simultaneous operation of the two teams, the driving speed has increased significantly. After the 13th month, since only one team is operating, the driving speed has dropped to the original level. It can be seen from Figure 6(b) that at the current driving speed, the driving team’s driving footage reached 2,441.765 meters in the 24th month, which is still nearly 60 meters away from completing the task. However, after the completion of the construction period, the working face will take time to move. Therefore, it has basically completed the driving task. Figure 6(c) shows that the improved driving plan completes the task 3 months earlier than the original plan. Figure 6(d) shows that the replacement of mining and driving has basically worked properly.


## 3.1. Parameter Setting

### 3.1.1. Weight Determination

The weight values of each index are determined by the group AHP method except the two indexes of advanced roadway support level and the moving level of two sides of roadway, which were given by experts. The weight results are shown in Table2.Table 2
Simulation model weight assignment.
IndexWjWeight valueInfluence of geological structure level of heading face on heading speedW10.3176Influence of working environment level of heading face on heading speedW20.3175Influence of normal level of heading face on heading speedW30.123Influence of starting level of heading face on heading speedW40.123Influence of mechanical equipment fault level coal mining on heading face speed in heading faceW50.1189Mine pressure level of heading faceW60.3008Gas level of heading faceW70.4901Mine water level of heading faceW80.2091Roof breakage degree of heading roadwayW90.4012Roof subsidence level of heading faceW100.2908Approaching level of two sides of tunneling roadwayW110.308Influence of geological structure level on advancing speed of coal mining faceW120.2986Influence of working environment level of coal mining face on advancing speed of working faceW130.1962Influence of the normal level of the working face on the advancing speed of the working faceW140.0616Influence of mining face level on working face advancing speedW150.0792Influence of mechanical equipment failure level of coal mining on working face advancing speedW160.0835Influence of roadway end management level on working face advancing speedW170.099Influence of driving speed on advancing speed of working faceW180.1819Mine pressure level of coal mining faceW190.3325Gas level in coal mining faceW200.5278Mine water level of coal mining faceW210.1397Roof breakage degree of coal mining faceW220.3062Roof subsidence level of coal mining faceW230.3012Roof deformation level of coal mining faceW240.3926Roadway advanced support levelW250.35The approaching level of two sides of roadwayW260.65


### 3.1.2. Basic Parameter Setting

The basic parameters mainly include the distance from the tunneling face to the geological structure zone, the failure rate of the tunneling face mechanical equipment, the gas concentration of the return air flow from the heading face, the distance from the mining face to the geological structure zone, the normal circulation rate of the mining face, etc.


### 3.1.3. Assignment of Auxiliary Variables

The level of each influencing factor changes as the value of the factor changes. This paper uses table function to express the relationship between the level and the actual value of the influencing factors. The level of influencing factors was quantified with 40∼100 scores. The influence degree is represented by the promotion or weakening effect on the numerical value of the affected factor, whose value may be greater than or less than 1. Considering the research results of many scholars, the paper defines the range of influence to be between 0.5 and 1.5 [20].


## 3.1.1. Weight Determination

The weight values of each index are determined by the group AHP method except the two indexes of advanced roadway support level and the moving level of two sides of roadway, which were given by experts. The weight results are shown in Table2.Table 2
Simulation model weight assignment.
IndexWjWeight valueInfluence of geological structure level of heading face on heading speedW10.3176Influence of working environment level of heading face on heading speedW20.3175Influence of normal level of heading face on heading speedW30.123Influence of starting level of heading face on heading speedW40.123Influence of mechanical equipment fault level coal mining on heading face speed in heading faceW50.1189Mine pressure level of heading faceW60.3008Gas level of heading faceW70.4901Mine water level of heading faceW80.2091Roof breakage degree of heading roadwayW90.4012Roof subsidence level of heading faceW100.2908Approaching level of two sides of tunneling roadwayW110.308Influence of geological structure level on advancing speed of coal mining faceW120.2986Influence of working environment level of coal mining face on advancing speed of working faceW130.1962Influence of the normal level of the working face on the advancing speed of the working faceW140.0616Influence of mining face level on working face advancing speedW150.0792Influence of mechanical equipment failure level of coal mining on working face advancing speedW160.0835Influence of roadway end management level on working face advancing speedW170.099Influence of driving speed on advancing speed of working faceW180.1819Mine pressure level of coal mining faceW190.3325Gas level in coal mining faceW200.5278Mine water level of coal mining faceW210.1397Roof breakage degree of coal mining faceW220.3062Roof subsidence level of coal mining faceW230.3012Roof deformation level of coal mining faceW240.3926Roadway advanced support levelW250.35The approaching level of two sides of roadwayW260.65


## 3.1.2. Basic Parameter Setting

The basic parameters mainly include the distance from the tunneling face to the geological structure zone, the failure rate of the tunneling face mechanical equipment, the gas concentration of the return air flow from the heading face, the distance from the mining face to the geological structure zone, the normal circulation rate of the mining face, etc.


## 3.1.3. Assignment of Auxiliary Variables

The level of each influencing factor changes as the value of the factor changes. This paper uses table function to express the relationship between the level and the actual value of the influencing factors. The level of influencing factors was quantified with 40∼100 scores. The influence degree is represented by the promotion or weakening effect on the numerical value of the affected factor, whose value may be greater than or less than 1. Considering the research results of many scholars, the paper defines the range of influence to be between 0.5 and 1.5 [20].


## 3.2. Simulation Analysis

### 3.2.1. Model Checking

Vensim software has its own “check model” function, which can check the consistency of model units. The simulation model in this paper runs well, and the consistency of the model units is verified by the “check model” function. After running the model, it can be obtained that the simulated output value of driving length is 2,441.765 meters, and the error rate is −2.33% compared with the actual value of 2500 meters. Compared with the actual value of 2450 meters, the simulated output value of the advance distance of the coal mining face is 2526.87 meters, and the error rate is +3.14%. The absolute error value of both is less than 5%, indicating that the model has a good effect and is qualified.


### 3.2.2. Analysis of Simulation Results

(1)
The relationship between the level of each influencing factor and the driving speed is shown in Figures3(a) and 3(b). It can be seen from the figure that all influencing factors have an impact on the driving speed. The level of geological structure reached its lowest point in the 5th and 19th months, and the level of other influencing factors reached their lowest in the 6th and 20th months. Driving speed also reached its lowest point in the 6th and 20th months correspondingly. The curve trend of the driving speed was roughly the same as the horizontal curve trend of each influencing factor.(2)
The relationship between the level of each influencing factor and the advancing speed of the working face is shown in Figures4(a) and 4(b). It can be seen from the figure that various factors have an impact on the advancing speed of the working face. The level of geological structure remained at a relatively low level from 4th to 9th and 16th to 18th months, while other factors also fluctuated at a low level during this period. The driving speed reached its lowest point in the 6th and 17th months, and the trend of the tunneling speed curve of the working face was roughly consistent with that of the horizontal curve of the influencing factors.(3)
The simulation results of the mining-driving system are shown in Figures5(a)–5(d). From Figure 5(a), it can be seen that the advancing speed of the 9102 working face has been significantly slower since April. Also, for nearly 5 months thereafter, its speed remained at a low level, rose in the 11th month, and then significantly decreased in the 17th month. But the duration is short. It returned to above the normal advancing speed in the 20th month. The main reason is that the 9102 working face was getting closer and closer to the F304 fault in the fourth month, and the span of F304 is 490 meters. Therefore, under the influence of the fault, the driving speed remained at a low state for a long time. In the 17th month, the working face encountered collapse columns H51 and X14 one after another. The two collapse columns were adjacent and the span was relatively small. Therefore, the impact time on the speed was relatively short. Correspondingly, the working face advance distance curve slowed down as its speed decreased. It can be seen from Figure 5(b) that the driving speed of the heading face decreased significantly in the 4th month, reached the lowest point in the 6th month, then decreased significantly in the 16th month, and reached the lowest point in the 20th month. The lowest speed was lower than that in the 6th month. The main reason is that the heading face is getting closer to anomaly area 3 in 4 and 16 months, and the return air roadway traverses the anomaly area, so it has a greater influence on it. As a result, the driving speed is lower and the duration is longer in the 20th month, thus affecting the normal replacement between roadway driving. The driving footage is the cumulative amount of the driving speed. According to the driving footage curve, it can also be seen that the trend of the driving footage curve becomes slow in the month with a lower driving speed. It can be seen from Figure 5(c) that under normal circumstances, the driving speed is lower than the advancing speed of the working face, and it has a certain impact on the advancing speed of the working face. The trends are not exactly the same because the advancing speed of working face is also affected by other factors. Figure 5(d) shows that the fully mechanized mining team completed the entire construction period in the 24th month, while the tunneling team completed the entire construction period in the 27th month, which is 3 months later than the mining face, and there is an imbalance problem between mining and driving.Figure 3
Relationship between driving speed and influencing factors.
(a)(b)Figure 4
Relationship between the level of each influencing factor and the advancing speed of coal mining face.
(a)(b)Figure 5
Simulation results of mining-driving system: (a) relationship between advancing speed and advancing distance of working face; (b) relationship between driving speed and driving footage; (c) relationship between driving footage and advancing distance of working face; (d) relationship between driving speed and advancing speed of working face.
(a)(b)(c)(d)


## 3.2.1. Model Checking

Vensim software has its own “check model” function, which can check the consistency of model units. The simulation model in this paper runs well, and the consistency of the model units is verified by the “check model” function. After running the model, it can be obtained that the simulated output value of driving length is 2,441.765 meters, and the error rate is −2.33% compared with the actual value of 2500 meters. Compared with the actual value of 2450 meters, the simulated output value of the advance distance of the coal mining face is 2526.87 meters, and the error rate is +3.14%. The absolute error value of both is less than 5%, indicating that the model has a good effect and is qualified.


## 3.2.2. Analysis of Simulation Results

(1)
The relationship between the level of each influencing factor and the driving speed is shown in Figures3(a) and 3(b). It can be seen from the figure that all influencing factors have an impact on the driving speed. The level of geological structure reached its lowest point in the 5th and 19th months, and the level of other influencing factors reached their lowest in the 6th and 20th months. Driving speed also reached its lowest point in the 6th and 20th months correspondingly. The curve trend of the driving speed was roughly the same as the horizontal curve trend of each influencing factor.(2)
The relationship between the level of each influencing factor and the advancing speed of the working face is shown in Figures4(a) and 4(b). It can be seen from the figure that various factors have an impact on the advancing speed of the working face. The level of geological structure remained at a relatively low level from 4th to 9th and 16th to 18th months, while other factors also fluctuated at a low level during this period. The driving speed reached its lowest point in the 6th and 17th months, and the trend of the tunneling speed curve of the working face was roughly consistent with that of the horizontal curve of the influencing factors.(3)
The simulation results of the mining-driving system are shown in Figures5(a)–5(d). From Figure 5(a), it can be seen that the advancing speed of the 9102 working face has been significantly slower since April. Also, for nearly 5 months thereafter, its speed remained at a low level, rose in the 11th month, and then significantly decreased in the 17th month. But the duration is short. It returned to above the normal advancing speed in the 20th month. The main reason is that the 9102 working face was getting closer and closer to the F304 fault in the fourth month, and the span of F304 is 490 meters. Therefore, under the influence of the fault, the driving speed remained at a low state for a long time. In the 17th month, the working face encountered collapse columns H51 and X14 one after another. The two collapse columns were adjacent and the span was relatively small. Therefore, the impact time on the speed was relatively short. Correspondingly, the working face advance distance curve slowed down as its speed decreased. It can be seen from Figure 5(b) that the driving speed of the heading face decreased significantly in the 4th month, reached the lowest point in the 6th month, then decreased significantly in the 16th month, and reached the lowest point in the 20th month. The lowest speed was lower than that in the 6th month. The main reason is that the heading face is getting closer to anomaly area 3 in 4 and 16 months, and the return air roadway traverses the anomaly area, so it has a greater influence on it. As a result, the driving speed is lower and the duration is longer in the 20th month, thus affecting the normal replacement between roadway driving. The driving footage is the cumulative amount of the driving speed. According to the driving footage curve, it can also be seen that the trend of the driving footage curve becomes slow in the month with a lower driving speed. It can be seen from Figure 5(c) that under normal circumstances, the driving speed is lower than the advancing speed of the working face, and it has a certain impact on the advancing speed of the working face. The trends are not exactly the same because the advancing speed of working face is also affected by other factors. Figure 5(d) shows that the fully mechanized mining team completed the entire construction period in the 24th month, while the tunneling team completed the entire construction period in the 27th month, which is 3 months later than the mining face, and there is an imbalance problem between mining and driving.Figure 3
Relationship between driving speed and influencing factors.
(a)(b)Figure 4
Relationship between the level of each influencing factor and the advancing speed of coal mining face.
(a)(b)Figure 5
Simulation results of mining-driving system: (a) relationship between advancing speed and advancing distance of working face; (b) relationship between driving speed and driving footage; (c) relationship between driving footage and advancing distance of working face; (d) relationship between driving speed and advancing speed of working face.
(a)(b)(c)(d)


## 3.3. The Optimization of Mining and Driving Replacement

It can be seen from the simulation results that there is an imbalance problem in the replacement of mining and driving. Due to the relatively complicated process of the tunneling face itself, the normal driving speed is relatively slow. Secondly, the tunneling face encountered anomaly area 3 in the 4th and 16th months, which affected the normal driving speed. For the reason that most of the factors affecting the driving speed are the changes of geological conditions, resulting in the change of other factors, and most of them belong to the mine hydrogeological factors, it is difficult to improve the driving speed by changing these factors. In order to ensure the balance of mining replacement, the following two aspects can be considered.(1)
Advance tunneling period: from Figure5(d), it can be seen that the coal mining team completes the construction period 3 months earlier than the tunneling team. If the tunneling team starts working three months earlier, the construction period will be completed at the same time as the mining team to ensure the normal replacement of mining.(2)
Arranging more tunneling teams at the right time: when the resources are sufficient, two driving teams can be reasonably arranged to work. The tunneling working face is composed of transport lane and return air lane, and the two roadways do not affect each other when working, which can meet the requirements of two tunneling teams working at the same time. Through further simulation, the paper found that in the 11th month, a new driving team was added to start work in the transport lane, and the two teams worked in parallel from 11 to 13 months, and after the 13th month (return air lane driving task is completed), the remaining roadway driving tasks are completed independently by the second team. The specific simulation results are shown in Figure6.Figure 6
Simulation results of mining plan arrangement: (a) relationship between original driving speed and current driving speed; (b) relationship between current driving speed and driving footage; (c) relationship between original driving footage and current driving footage; (d) relationship between working face advancing distance and current driving footage.
(a)(b)(c)(d)It can be seen from Figure6(a) that in the 11th month, due to the simultaneous operation of the two teams, the driving speed has increased significantly. After the 13th month, since only one team is operating, the driving speed has dropped to the original level. It can be seen from Figure 6(b) that at the current driving speed, the driving team’s driving footage reached 2,441.765 meters in the 24th month, which is still nearly 60 meters away from completing the task. However, after the completion of the construction period, the working face will take time to move. Therefore, it has basically completed the driving task. Figure 6(c) shows that the improved driving plan completes the task 3 months earlier than the original plan. Figure 6(d) shows that the replacement of mining and driving has basically worked properly.


## 4. Conclusions

Based on the system dynamics modeling, the dynamic simulation structure model of mining-driving system in Wangzhuang Mine was established with the help of Vensim software, and the dynamic simulation of 9102 mining face and 7106 heading face was carried out. The results show that there exists the phenomenon of replacement imbalance between mining and driving. Advancing the driving schedule or adding a new driving team can meet the demand of normal replacement. It can provide guidance for making driving and mining planning in other similar working faces. Optimizing the driving process reasonably helps to realize the perfect coordination of mining face and roadway driving in time and avoid the emergence of mining face replacement tension and resource waste phenomenon, which ensures the mine efficiency and production safety to the utmost extent.



---

*Source: 9636178-2021-09-10.xml*