Numerical Simulation of Mining-Induced Stress Evolution and Fault Slip Behavior in Deep Mining

The ground pressure distributes signiﬁcant variation in underground mining near fault. Fault reactivation is an important factor to induce the rock burst. Therefore, characterizing geological settings in mining areas by the geological information can improve the accuracy of simulation. To investigate the characteristic of mining stress evolution and reactivation of the F16 reverse fault during the retreat Mining-Induced s in Yima coalﬁeld, a three-dimensional digital elevation model based on GIS platform was applied. The 3D geological model includes three working faces, and F16 fault was constructed by AutoCAD software. Then, the 3D geological model was imported into the FLAC 3D code to simulate the potential of mining-induced fault reactivation. The simulation results illustrate that the footwall of F16 fault is a high stress concentration area. Aﬀected by F16 fault and the huge thick gravel rock in the roof, the coal seam near the fault accumulates a large amount of elastic strain energy, which increases the potential of rock burst hazards in the process of mining.


Introduction
With the continuous increase of the depth and intensity of coal mining, the mining engineering layout tends to be complicated and the number of deep mines keeps increasing in China [1]. Rock burst is one of the most hazardous phenomena that accumulates huge elastic energy in surrounding rock and releases instantaneously in the process of mining, which poses a great threat to the safety of personnel and equipment in coal mine. Faults are widely distributed in coal measure strata. In the process of coal seam mining, hidden geological structures such as faults are often encountered. Overburden movement and mining-induced stress distribution show significant variation when the working face in mining is near fault. e considerable field practice shows that deep mining near fault is able to easily induce rock burst and mine earthquakes, and the frequency and degree of damage are significantly increased [2][3][4][5]. e conflict between coal mining and safety production is prominent in deep mining near fault. How to effectively prevent rock burst is an urgent problem to be solved in deep mining. us, studying the law of fault activation induced by mining disturbance has great significance on disaster prevention.
Although a lot of research has been done on the mechanism of rock burst, it is still unable to accurately predict the specific time and location of rock burst event. e most effective measure is to optimize the design of mining method and mining sequence. At the early stage of coal seam mining, the evaluation of rock burst hazard in the mining area is carried out to define the potential rock burst area and the corresponding risk degree, which can provide a basis for the monitoring and prevention of rock burst in the mining process. At present, the commonly used methods of evaluation of rock burst hazard include the comprehensive index method, which combines the actual geological and mining technical conditions of the mining area and the influence weights of various factors, to evaluate the rock burst risk. According to the stress state and the properties of coal and rock mass, there is a possibility of index diagnosis method of rock burst, combined with coal and rock mass, geological structure, overburden structure, and other factors of coupling evaluation method [6][7][8][9][10][11]. ese methods can assess the rock burst hazard of mining area from a macroscopic perspective, but they do not consider the influence of mining disturbance and cannot realize the dynamic assessment of the rock burst risk accompanied by the mining changes.
By means of numerical simulation, the evolution behavior of mining-induced stress under different mining conditions is analyzed, and rock burst area can be divided. Numerical simulation is widely used because it is easy to operate and can predict the rock burst hazard before the implementation of mining projects. However, most of the existing numerical simulation studies on coal seam mining have been properly simplified, instead of using actual stratigraphic data to construct three-dimensional models, which makes the simulation results unable to accurately distribute the ground pressure of coal mining affected by the geological factors.
Numerical simulation technology based on 3D fine geological modeling uses actual stratigraphic data to construct 3D models, which is widely used in oil and gas storage, water conservancy, slope, and open pit mining engineering [12,13]. However, it is seldom used in the analysis of ground pressure distribution in underground coal mining. In this study, the geological mining conditions of No. 21 mining area of Qianqiu mine in Yima coalfield are taken as the engineering background, and the numerical simulation research based on fine modeling is carried out to analyze the mining-induced stress evolution and fault slip behavior during mining approaching to F16 fault, and the dynamic rock burst risk regionalization was further discussed. It can provide a theoretical basis for the deployment design of deep mining near fault. F16 fault is in a high horizontal tectonic stress environment, and at the same time affected by the mining disturbance of each mine, the superposition of two effects lead to fault activation frequently. Fault activation can be caused by tectonic stress nappe. But the activation cycle is long, at least for decades. It is concluded that the main factor of fault activation is coal mining. Statistics show that among the 107 cases of rock bursts, there are only 8 cases that fall in the area, which is 1000 m or farther away from the F16 fault with a mining area more than 50 km 2 . For the mining areas, the range of 6 km 2 is 500 to 1000 m away from F16 fault. Rock bursts occurred 48 times. For the mining area of 3 km 2 , where the distance toward the fault is less than 500 m, rock bursts occurred 51 times [14,15]. Consequently, the distance of mining area between the F16 fault and the risk of rock bursts is highly correlated in an inverse correlation.  Figure 2.

Working
e complex geology and mining conditions in the west wing of No. 21 mining area lead to the frequent occurrence of rock bursts and serious damage in the process of coal seam mining, such as the "6.5" accident of 21201 working face and the "11.3" accident of 21221 working face, resulting in heavy casualties and huge economic losses. e cross section of the roadway shrank instantaneously, the partial cross section became less than 1 m 2 , the coal seam had a dull sound before rock burst, and the gas concentration exceeded the limit after rock burst [16].

Construction of the 3D Numerical Model
Numerical calculation has obvious characteristics of quantification, repeatability, and economy, so it is widely used in ground pressure distribution analysis. e calculated results can make assessment of rock burst risk and also provide basic information or basis for other methods. At the same time, it can better reflect the actual situation of the project if it is supplemented by on-site monitoring. e numerical analysis method based on refined modeling includes the following aspects: determination of engineering geological conditions and in situ stress state, construction of refined numerical model reflecting actual stratum data, selection of mechanical model and surrounding rock mechanical parameters, excavation and support according to mining design, analysis of simulation results after numerical calculation, combined with on-site mine pressure behavior characteristics and various types of rock mass monitoring information, and assessment of rock burst risk in mining areas.

Construction of the Numerical Model.
Based on Auto-CAD platform, a 3D geological model was constructed by using spatial sequence profile and DEM hybrid method [17]. e flow chart of generation of 3D geological model is shown in Figure 3. Firstly, the borehole information, geological profile, and floor contour data are extracted and summarized. e relationship and regularity between the existing data are analyzed. Secondly, the initial data obtained are processed according to certain mathematical methods in order to increase the amount of data. e roof strata, coal seam, and floor are successively constructed to form a complete stratum model. Finally, through a certain optimization combination, the spatial geometry of the orebody and geological structure can be reflected realistically, and the interconnection between the surfaces can form a three-dimensional model. e construction process of 3D geological model is shown in Figure 4. ANSYS software is used to generate 3D point line model, and then the point and line model is transformed into three-dimensional model, and the mesh is divided through the ANSYS-FLAC3D transform program to import the three-dimensional entity model into FLAC3D software. As shown in Figure 5, the final numerical model (length × width × height) is configured at the size of 1150 × 1198 × 979 m.
ere are 264,580 nodes and 622,410 cells. e +X direction is the mining direction of the working face, +Y is the dip direction of the working face, +Z direction is vertical to the direction of the model. Working faces are arranged along the dip direction successively, among which the No. 21221 working face is the closest to the F16 fault. e upper bound of the model Z is free, and the upper boundary is the surface. e vertical model is fixed at the bottom, and the horizontal displacement of the side boundary in X-and Y-directions is limited.

In Situ Stress Inversion.
e in situ stress in the depth of 700 meters in Qianqiu mine is measured by the method of hydraulic fracturing. e in situ stress distribution is shown in Table 1 [18].
From the table, it can be seen that the stress field type is generally σ H > σ V > σ h . Tectonic stress is the main stress in the stress field of mining area. Using in situ stress data, we can calculate the gradient at different location and gradient value. In general, the gradient constraint is applied to the boundary of the model, and then the stress field of the model is generated by the model self-weight stress. e directions of maximum principal stress and the minimum principal stress on the X-Y-plane are coincident with the X-and Yaxes, and the Z-axis coincides with the vertical stress. e existence of the shear stress is neglected.       ere is a basic assumption of the stress loading method in the paper, namely, σ z ′ � σ z � σ 2 which indicates that the intermediate principal stress in the original space coordinate system is the same as the intermediate principal stress in the space coordinate system. e basic assumption in this paper coincides to the in situ stress in Qianqiu mine measured by hydraulic fracture. Inversion of in situ stress loading is adopted to improve the accuracy of simulation results. e stress loading method needs to measure the space stress conversion in the north and east direction for the X-axis and the Y-axis of the coordinate system; namely, the three-dimensional model of the X-axis and the Y-axis direction is established in this paper. By means of the spatial stress conversion formula in elasticity, the measured spatial in situ stress is converted to the coordinate system where the geological model is located, and the in situ stress in this coordinate system has shear stress: y′x + σ y l 2 y′y + σ z l 2 y′z + 2 τ xy l y ′ x l y′ + τ yz l yy′ l y ′ z + τ zx l y′z l y′x , σ z ′ � σ x l 2 z′x + σ y l 2 z′y + σ z l 2 z′z + 2 τ xy l z′x l z′y + τ yz l z′y l z′z + τ zx l z′z l z′x , τ xy � σ x l x′x l y′x + τ xy l y′x l x′y + l x′x l y′y + σ y l xy l yy + τ yz l y′y l xz′ + l x′y l yz′ + σ z l xz l yz + τ zx l yz′ l x′x + l xz′ l y′x , τ zx � σ x l z′x l x′x + τ xy l x′x l zy′ + l z′x l xy′ + σ y l z′y l x′y + τ yz l xy′ l z′z + l zy′ l x′z + σ z l z′z l x′z + τ zx l x′z l z′x + l z′z l x′x , τ y′z′ � σ x l y′x l z′x + τ xy l z′x l y′y + l yx′ l z′y + σ y l y′y l z′y + τ yz l z′y l y′z + l y′y l z′z + σ z l y′z l z′z + τ zx l z′z l y′x + l y′z l z′x ,

Numerical Calculation Scheme.
Due to the complexity of the model, each rock layer is not a horizontal plane, but a curved surface with a certain spatial shape, which also increases the difficulty of fault contact surface. e method of this paper is to merge the hanging wall and footwall strata of the fault, respectively, and then separate the strata after adding interface. In the numerical model, Coulomb failure criterion is used to evaluate shear failure on the contact surface [19]. Mechanical properties of stratus are shown in Table 2. Properties of contact surface are defined as follows: normal stiffness is 2 GPa/m, shear stiffness is 5 GPa/m, internal friction angle is 20 degrees, and cohesion is 0.5 MPa. e mining sequence of the working face is, firstly, mining 21181 working face, then mining 21201 working face, and, finally, mining 21221 working face. Simulation of excavation step distance is set to 40 m. e mining process    Table 3. In the numerical model, monitoring lines are both arranged along the coal seam strike and inclination to monitor the stress distribution characteristics of original rock before mining and the dynamic variation of mininginduced stress in the process of mining.

Original Stress Distribution.
In the simulation, original stress equilibrium is reached, and then the working face is mined. e maximum and minimum principal stress distribution of the model is shown in Figure 6. e principal stress of each stratum is compressive stress, and the stress level is gradually increased from the upper to the lower of the model. e fault has a certain influence on the original stress distribution. e stress value changes and breaks on the fault surface and return to normal value far away from the fault. Figure 7 shows the original rock stress distribution of coal seam. It can be seen that the stress distribution of coal seam is uneven, and the closer to the fault along the coal seam dip direction, the greater the stress value. In order to analyze the stress distribution, the stress value is extracted along the monitoring line of coal seam inclination and drawn into a curve, as shown in Figure 8.
It can be seen from the original rock stress curve that the coal seam stress distribution is discontinuous under the influence of coal seam dip angle and fault. e vertical stress is 18-22 MPa, the maximum horizontal stress is 20-28 MPa, and the minimum horizontal stress is 12-16 MPa. e stress value near the fault is higher in the deep coal seam. e horizontal stress of coal seam is higher than the vertical stress. e distribution of stress field represents the regulation that σ H > σ V > σ h . e stress field in the mining area is mainly the tectonic stress field and the results of the inversion stress are similar to in situ stress measured in field. e following calculation and analysis are based on the results above.

Mining-Induced Stress Evolution.
e vertical stress distribution of coal seam in the process of mining is shown in Figure 9. e vertical stress of 21181 working face, 21201 working face, and 21221 working face is only listed in the limited space when the face is mined to different positions.
Stress concentration is formed in the front and side of the coal wall in the process of advancing the working face, and the goaf is the stress reduction area.
e results are consistent with the stress shell theory [20]. e magnitude of the mining-induced stress, which is formed when the mining face is extracted, gradually increased during the mining process. e measuring point records the change of abutment stress in the process of advancing the working face and matches abutment pressure curves of each working face in the mining process, as illustrated in Figure 10.
It can be seen that the initial vertical stress of 21181 working face is 18.2 MPa, and the peak stress reaches 24.2 MPa after being affected by mining. e stress increase is 32.9%, and the stress concentration coefficient is 1.32. When mining 21201 working face, the original rock stress rises to 23.1 MPa, the peak stress is 32.3 MPa, the stress increase is 39.8%, and the stress concentration coefficient is 1.39. Further approaching fault mining, when mining 21221 working face, the original rock stress rises to 28.7 MPa, the peak stress is 39.2 MPa under the influence of mining, the stress increase is 36.6%, and the stress concentration coefficient is 1.36. When the original rock stress was before the work of 21221, the stress of coal body has exceeded the uniaxial compressive strength of coal seam. It has been tested in laboratory that the strength of No. 2 coal seam in Qianqiu mine is between 15 and 25 MPa, and the average strength is 22.47 MPa, which is close to the stress condition of rock burst. Under the dynamic disturbance, rock burst is easy to occur. e horizontal stress and vertical stress change rules of coal seam during mining were recorded along the coal seam dip line, as shown in Figure 11.
It can be seen that the peak value of horizontal stress in the dipping direction of 21181 working face is 22.5 MPa and the peak value of vertical stress is 27.8 MPa. When mining 21201 working face, the peak value of horizontal stress in the inclined direction of working face is 33.4 MPa, the peak value of vertical stress is 38.9 MPa, the increase rate of horizontal stress is 48.4%, and the increase rate of vertical stress is 39.9%. e peak value of horizontal stress is 23.5 MPa and the peak value of vertical stress is 32.2 MPa at       Figure 12 presents the principal stress and displacement vector when mining 21201 working face. e footwall of F16 fault is a high stress concentration area. e rock mass in the footwall of the fault twists to the goaf with the lateral coal seam of 21201 working face as the support. Affected by F16 fault and the huge thick gravel rock in the roof, the coal seam near the fault accumulates a large amount of elastic strain energy. It can be judged that the coal seam in 21201 working face is in the stress limit state and the risk of rock burst is high.    Advances in Materials Science and Engineering

Fault Slip Behavior under
Mining. Make a profile along the coal seam inclination, and the variation shear stress in surrounding rock during the mining process of the working face is shown in Figure 13.
It can be seen from that during the process of mining near the fault, the change of fault stress at the roof position is mainly caused by the increase of shear stress. When mining the 21181 working face, the shear stress on both hanging wall
Monitoring points are arranged on the fault to obtain the spatial and temporal distribution of normal stress and shear stress in different layer during the process of the mining. Mining sequence of working faces is shown in Table 3, x0 is the original stress equilibrium process, x1-x15 is the mining 21181 working face, y1-y23 is mining 21201 working face, and z1-z22 is mining 21221 working face. During the mining process, the normal stress and shear stress variation of the monitoring points are shown in Figure 14.
e initial state of normal and shear stress of each measuring point increases from the overburden to the floor which conforms to the distribution law of in situ stress. e normal stress is greater than the shear stress. e change time of each layer on fault is different under mining; stress variation of monitoring points is shown in Figure 15. Slip mechanism of the fault can be divided into two types. One type is the normal stress which is continuously reduced, and the shear stress increases rapidly.
is situation mainly occurs in rock strata at higher position, such as A2-A4 points. Another type is both normal stress and shear stress increase, but the increase range of shear stress is larger than that of normal stress. is situation mainly occurs in coal seam area, such as A5 point. e stress difference on the fault   plane is related to the stress evolution distribution of surrounding rock after coal mining. e farther away from the coal seam, the earlier the stress change starts. e change of shear stress is mainly increased, while the change of normal stress is mainly decreased, and the change of shear stress is earlier and greater than that of normal stress; the response of shear stress to mining is more sensitive [21][22][23][24]. e closer the distance from the mining seam is, the higher the normal stress and shear stress are. e reason is that the movement and evolution of overburden caused by mining disturbance is highly related [25][26][27][28][29]. e shear displacement of each point is shown in Figure 16.
With working face close to fault, fault shear displacement is growing rapidly. e closer the fault is to the coal seam, the greater the slip is. e farther away from each other, the more the damage and the slip of rock strata gradually weaken. e fault slip increases continuously, and the risk of fault slip and instability is greater. erefore, it can be concluded that the risk of rock burst induced by fault reactivation is the greatest when mining 21221 working face.

Conclusion
Based on the geological mining conditions of Qianqiu mine as engineering background, the refined numerical simulation research was carried out to analyze the surrounding rock stress evolution in the mining process near F16 fault. e following conclusions were drawn.
(1) According to the geological borehole and F16 fault profile information of Qianqiu mine, a three-dimensional refined numerical model reflecting the actual stratigraphic conditions of No. 21 mining area in Qianqiu mine was constructed. (2) Based on the in situ stress data of No. 21 mining area, the stress is applied to the coordinate system of the numerical model through the spatial stress transformation formula. Numerical simulation of the initial balance calculation of the working face is before mining. e vertical stress is 18-22 MPa, the maximum horizontal stress is 20-28 MPa, and the minimum horizontal stress is 12-16 MPa. e distribution of stress field represents a regulation that σ H > σ V > σ h . e stress field in the mining area is mainly the tectonic stress field, and the results of the inversion stress are similar to in situ stress measured in field, which verifies the rationality of the stress loading mode. (3) It shows that there is a critical distance when mining near fault area. When evaluating the influence of the fault on the working face, it should be comprehensively considered the fault occurrence. Under the influence of fault dip angle and mining, the stress state of each measuring point on the F16 fault plane is different, and the normal stress on the upper fault plane tends to decrease, while the shear stress on the upper fault plane tends to increase. e closer the working face is to the fault, the more obvious the disturbance effect of mining on the fault is. e footwall of F16 fault is a high stress concentration area. Affected by F16 fault and the huge thick gravel rock in the roof, the coal seam near the fault accumulates a large amount of elastic strain energy, which increases the potential of rock burst hazards in the process of mining.
Data Availability e data used for conducting classifications are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.