Collaborative Mining Sequence Optimization for Multiple Stopes under Intensive Mining

*e optimization of a mining sequence not only reduces stress concentration in surrounding rock but also prevents underground debris flows, significantly improving safety. Firstly, the 870–898 m level of the eastern mining area in the Tiaoshuihe phosphate mine was divided into 25 ore blocks, and six different mining sequences were designed for this area. *en, it was calculated that five ore blocks must be processed simultaneously to reach the annual production output. *e distances between the five simultaneously mined ore blocks will inevitably affect the efficiency of the equipment for any scheme. So, a collaborative model considering both the area stability and production capacity was established by combining the distance between the centers of the five ore blocks as an index. Differences in stability, deformation, and plastic zone size between the schemes are compared. *e calculation results show that a mining scheme with a convex stepped shape produces the best results. *ese results provide a general method for entropy-based mining sequence optimization and an optimal solution for the Tiaoshuihe phosphate mine.


Introduction
Filling mining is a popular contemporary technique that is increasingly advocated by many governments. Ground pressure in these types of mines will increase considerably as mining continues to deeper depths, so the corresponding requirements for ground pressure control will continue to expand [1,2]. If measures are only taken after ground pressure triggers surrounding rock damage, production costs will increase and the production efficiency will be reduced. To better compensate for ground pressure, active control must be used instead of passive control. Implementing a reasonable mining sequence is an active ground pressure control method [3].
In 1963, Goodman and Brown put forward the concept of rock engineering excavation steps, which is described as the mining sequence now [4]. In 1972, Zouanalyzed the influence of construction steps on the construction of underground tunnels and proposed the concept of simulated excavation steps for the first time but did not combine finite element calculations and simulations [5]. Liu successfully optimized the mining sequence of the Dongguashan deep well by considering the production capacity, stress, displacement, and plastic zone of the surrounding rock [6]. Kang simulated different mining sequences in the middle level of the Xinli mining area of Sanshandao Gold Mine, obtaining characteristic values of the stress and displacement of rock under different mining steps. is work suggested that it is safe and efficient to mine one panel at the interval of two panels [7]. Based on the time-varying evolution of the mechanical parameters of a rock mass, Hu designed five mining sequence schemes in the deep section of the Gaofeng mine. Using FLAC3D, the spatiotemporal evolutions of rock mass mechanical parameters were simulated under different mining sequences. e results show that case 2 represents the best mining sequence [8].
e Tiaoshuihe phosphate mine is a gently inclined thin ore body with an average thickness of 3.22 m. is mine contains a mining area of approximately 23.42 km 2 and an estimated reserve of 200 million t [9,10]. To realize production, it is necessary to plan intensive mining, adjust the mining areas of simultaneous production in one location, and then maximize the utilization rate of personnel and equipment [11]. At present, excavating is mainly concentrated in the 870-898 m level of the eastern mining area. Optimization of the mining sequence in this area is necessary by considering regional stability and productivity objectives [12].

Geological Characteristics of Mining Area
Tiaoshuihe phosphate mining area is located in the north margin of Huangling fault area. e strata are inclined to NNE monoclinic structure. e dip angle is generally 4-8°.
e attitude of the stratum becomes steep due to the influence of fault structure. ere are 6 faults in the mining area, and all faults are normal. Among them, F1 and F4 are NNW trending faults. F3, F5, and F6 are NW trending faults, and F2 is the NNE trending fault. F1 is the eastern natural boundary of the mining area. F5 is the southern natural boundary of the mining area. e F4 fault divides the mining area into two natural blocks, the east and the west. It causes great damage to the continuity of the ore body. e model of ore body is shown in Figure 1.
870-898 m level is the main research area. 870-898 m level belongs to the main industrial phosphate bed, which is stable and widely developed. It is distributed at the bottom of Huji member and occurs in layers. e roof is thin to middle micritic dolomite in the lower submember of Huji member, intercalated with thin layer of gray (weathering color is grayish yellow) micritic mudstone. And the floor is thick bedded siliceous block in the upper submember of Zhangcunping member. e average thickness of 870-898 m level is 4.61 m. e grade of P 2 O 5 is 21.22-34.02%, and the average is 27.21%. Figure 2, the ore body was meshed by using a DIMINE-CAD-MIDAS/GTS coupling operation and then imported into the numerical simulation software FLAC3D for calculation. e fractures in the 870-898 m level of the eastern mining area are relatively developed, and the mining sequence is related to the overall stability and mining safety of the area [9]. According to geological occurrence conditions, a three-dimensional ore body model was established, and 25 numbered mining areas were identified. e width of each mining area was 60 m, and the length was 100 m. According to Saint-Venant's principle, an area 3 times greater than the ore body was selected to determine the model boundary [13,14].

Scheme Design.
Six schemes of mining sequences were designed for optimization. To ensure the comparability of the six schemes, all schemes completely mined and filled the 25 mining areas. Each scheme includes 50 calculation steps to simulate the process of mining and filling. en, the displacement, stress, and plastic zone results were analyzed. e six mining sequence schemes are shown in Figure 3. e numbers in the table indicate the mining sequence: scheme 1 mines along the short axis of the mining area, from east to west; scheme 2 mines along the long axis of the mining area, from south to north; scheme 3 mines from southwest to northeast; scheme 4 mines from northeast to southwest; scheme 5 mines northward a convex stepped shape; and scheme 6 mines in a "V" shape.
Numerical simulation of the mining area requires a mathematical model of the rock mechanical properties, which is usually the defined stress-strain relationship [15]. According to a geological investigation at the Tiaoshuihe phosphate mine and laboratory mechanical experimentation, the initial stress is considered to reflect only the selfweight of the rock mass. e Mohr-Coulomb constitutive model is commonly used, as follows [16]: where c is the cohesion in MPa, σ 1 and σ 3 are the principal stresses in MPa, and φ is the internal friction angle,°.

Results
Six command flows were compiled according to the mining sequence schemes described above [17]. Combined with the rock mass mechanical parameters in Table 1, the aforementioned three-dimensional ore body grid model was imported into the finite element numerical simulation software, FLAC3D. e maximum displacement, maximum stress, and plastic area volume were obtained after each mining step of each scheme.

Calculation Results of Six
Schemes. e displacement, tensile stress, compressive stress, shear stress, and plastic zone results calculated for six schemes are shown below. e overall trends of the six scheme results were similar, and the displacement results are shown in Figure 4.
Based on the displacement change curve of each step in schemes 1-6 ( Figure 4), the maximum displacement is 10.3 cm and occurs in step 19 of scheme 6. In this step, the displacement increases significantly, and the following displacement remains stable. e maximum displacement of the roof does not occur at the current goaf but instead occurs in the roof of the large filling body. erefore, the displacement of 10.3 cm in the roof of the large filling body will affect the whole project. e displacement of the roof in the 19th step is between 6-7 cm, which is also dangerous. During the actual engineering excavation, the area will be mined block by block. e displacement is expected to be less than 6 cm, so it is necessary to use bolts and anchor mesh for support after excavation. e plastic tensile failure zone curve under six schemes can be seen in Figure 5. e plastic tensile failure zone changes greatly after each step of mining. e maximum value is found to be 380182 m 3 , which arises after step 25 of scheme 4. e interface between the goaf and pillar exhibits high tensile stress, and the plastic tensile failure zone occurs mainly at the backfill. e maximum tensile plastic zone gradually expands with each step of mining for the first five steps and then expands more slowly after subsequent steps.

Analysis of the Maximum Values of Calculation Results
Using Different Schemes. Changes in the displacement, stress, and plastic zone after each step of mining in each scheme are shown below. Table 2 highlights the maximum displacements, stresses, and the final plastic zone total volumes under different schemes. In Table 2, D represents the displacement, TS represents the tensile stress, CS represents the compressive stress, SS represents the shear stress, SF represents plastic shear failure zone volume, TF represents tensile failure zone, TSF represents the total volume of the plastic shear failure zone, and TTF represents the total volume of the plastic tensile failure zone.
Among the six schemes, the displacement of scheme 6 is the largest, while that of scheme 5 is the smallest throughout, with relatively stable displacement in every step. e tensile stress of scheme 6 is also the greatest, while those of schemes 5 and 2 are greater than those other schemes. Scheme 6 has the greatest compressive and shear stresses, followed by scheme 1. e plastic shear failure zone of scheme 2 is the largest, followed by scheme 4. Scheme 2 produces the largest plastic zone of tensile failure, followed by scheme 3. e total volume of the plastic shear failure zone in scheme 1 is the largest, while that of scheme 6 is the smallest. Scheme 2 produces the largest total volume of the plastic tensile failure zone, while scheme 5 produces the smallest. e shear stress evolution trend is similar to the compressive stress. e maximum shear stress is 19.98 MPa after step 13, produced at the joint of the long axis side roof, pillar, and filling body.

Calculation of the Production Capacity of the Tiaoshuihe
Phosphate Mine. e current production plan at the Tiaoshuihe phosphate mine is 1 million t/a, and extracted ore is transported and exported by automobiles. Due to the high altitude of Zhangcunping Town, Yichang City, snow Advances in Civil Engineering covering in winter shuts down transportation and production. erefore, producing more than 4000 t/d is required to achieve the annual production target. e ore grade stratification of the mine is relatively apparent. e upper layer is high-grade ore, and the lower layer is low-grade ore. is ore is mined in layers to reduce the mixing rate. Each operation can explore only half of the working face. e footage of one blasting is 2.83 m, and the density of the phosphate rock is 3.021 g/cm 3 . e production of one blasting in a single stope is 89.77 t.
Because the mining method of excavating alternate strips is adopted, each ore block can be divided into eight strips. Each ore block can have at most four strips being mined at the same time, so the single production capacity of each ore block is 359.08 t. A two-shift system is used at the Tiaoshuihe phosphate mine, where each shift completes one working procedure of drilling, blasting, and ore drawing and transportation. e daily production capacity of each ore block is 718.16 t. erefore, according to the production target of 4000 t/d, at least 5.6 ore blocks need to be mined at the same time to meet the production capacity requirements. However, a small part of the 870-898 m level in the eastern mining area has a good ore grade and does not need to be mined in layers. erefore, to achieve the annual production, only five ore blocks need to be mined at the same time in this area.
Among the six different mining sequences, there are considerable differences in the distances between the five ore block centers. e discreteness of these five blocks has a significant influence on the cooperative use of equipment and operation cost in the mining process. e center of each group of five ore blocks mined at the same time is taken as the base point, and the sum of the distances from the center of each of the five ore blocks to the base point is called the center distance. en, the center distance is used to measure the discreteness of five simultaneously mined ore blocks. e   higher the discreteness, the lower the coordination degree of a given scheme. Figure 6 shows the center distance calculation for the first group of ore blocks in scheme 1. e distances from the five ore blocks to the base point are 0 m, 60 m, 60 m, 120 m, and 120 m. erefore, the center distance of the first group in scheme 1 is 360 m. Because the center distances of the five groups in scheme 1 are all 360 m, the sum of the center distances of scheme 1 is 1800 m.    Figure 6: Calculating the center distance of the first step in scheme 1.     Advances in Civil Engineering e center distance for the other schemes is calculated similarly. e sums of the five groups of center distances for each scheme are shown in Table 3.

Calculation of Coordination
Degree. According to previous calculation indicators, a collaborative model of mining sequence and production capacity target was established [18,19], as shown in Figure 7. e collaborative entropy method was used to calculate the coordination degree of various schemes. e value of entropy and the coordination degree of the system are inversely related. So, the larger the coordination entropy, the lower the coordination degree of the system. e smaller the coordination entropy is, the greater the coordination degree of the system is. According to equation (2), the basic elements in the collaboration model were evaluated and scored, and the calculation results are shown in Table 4 [20]: where θ is the score of the basic element and A i is the value of the basic element. Entropy value calculations of the basic elements follow equation (3) and are shown in Table 5.
where H(s) is the value of coordination entropy. Finally, according to equation (4), the coordination degrees of the basic elements of each scheme were calculated [17]: where C i is the coordination degree of the basic elements. e calculated values of the coordination degrees of the basic elements are shown in Table 6. e coordination degree of each scheme can be obtained by summing the coordination degrees of the basic elements. e calculation results are shown in Table 7.
From the results of the coordination entropy method, the highest coordination degree, 0.646, is found in scheme 5. e second highest is calculated from scheme 1, with a value of 0.644. erefore, the use of a convex stepped shape in scheme 5 is the best mining scheme.

Conclusions
An accurate model of the 870-898 m level in the eastern mining area of the Tiaoshuihe phosphate mine was built, and the model grid was divided for mining sequence optimization. Command flows were written according to the six mining sequences. e results of the six schemes were calculated by using the Mohr-Coulomb constitutive model. e annual production capacity of the mine was also calculated.
A collaborative model of the mining sequence and productivity target was established. Via the coordination entropy method, the final collaborative degree of the six schemes was determined. It was concluded that scheme 5 is the best, while scheme 1 is the second-best. erefore, it was recommended that the 870-898 m level of the eastern mining area is mined using a convex stepped shape. e theoretical research result of this study is that the idea of coordination is introduced into the process of mining sequence optimization and filling optimization. e basic parameters are brought into the coordination model by coordination entropy, and the coordination model includes all kinds of factors that affect the stability in the mining process. Economic indicators can also be taken into account in the future.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.