Onsite Tests and Numerical Simulation of Broken Rock Zones in Surrounding Rocks of Seepage Roadways under Blasting

To study the influence of blasting vibration on the broken rock zone around a seepage roadway and provide guidance for design of the roadway support, the broken rock zones around rock of seepage roadways under production blasting vibration are determined by onsite tests in a mining area, Daye iron mine. During the testing process, it is found that blasting vibration causes internal cracks of surrounding rocks to initiate and develop, the fracture density increases, the acoustic wave velocity of rock mass decreases, and the broken rock zones expand. At the same time, through onsite observation, it is found that blasting vibration results in crack development and formation of a water pathway to lead to surface water into the ground. The mechanical response around rock of the seepage roadway under blasting vibration is simulated by the two-dimensional realistic fracture progress analysis calculation software (RFPA). It is found that internal cracks of roadway surrounding rock initiate, propagate, and join up gradually, and the fracture range is expanding under the seepage water pressure, ground stress, and cyclic loads, and the broken rock zones also expand. The results from numerical simulation are consistent with the results of onsite tests. It is also found that the tensile stress appears around some cracks, leading to part of the cracks more likely to generate shear failure under the seepage water pressure during simulation.


Introduction
During deep mining, the surrounding rock has always high ground stress and high hydraulic pressure [1,2]. As the depth increases, the pressure of the confined water becomes larger. Mining and excavation activities can easily lead to the expansion of fractures or activation of faults, thus forming relatively concentrated seepage channels. When the roof, floor, or wall of a roadway suddenly is destroyed, high-pressure karst water will enter the roadway and cause serious water inrush disaster. Due to the disturbance of blasting operation, the damage of rock mass is aggravated, and the water inrush disaster has obvious characteristics of suddenness and unpredictability [3]. Therefore, it is significant to guide the mine support and prevention of mine water inrush by monitoring the broken zone of the surrounding rock in the roadways under blasting load [4][5][6].
Many previous researchers have provided qualitative descriptions about the broken zone of the surrounding rock. Jiang et al. proposed the theoretical method to predict the development around the rock plastic zone in soft rock tunnels and discussed the effect of the mechanical properties of soft rock on loosening pressure [7]. Han et al. measured the broken rock zone of the surrounding rock in Zhongyi tunnel [8]. Wang et al. investigated the excavation damaged zone around roadways during longwall mining in a coal mine and found that dynamic pressure and faults contribute to the expansion of excavation damaged zones [9]. Yang et al. proposed an analytical time-dependent viscoelastic model for the evolution of fractured zone in the surrounding rock mass based on sliding crack model and equivalent crack method and simulated excavation damage zone in the surrounding rock mass of circular tunnel under lining and hydrostatic condition by using RFPA 2D [10]. The broken rock zones are affected by these factors, such as the ground stress, geological conditions, hydrostatic condition, and excavation methods. However, there are few papers about the combined effect of the ground stress, hydrostatic conditions, and blasting excavation.
For a proper investigation of the width of broken rock zones in the surrounding rock, in situ tests and computer simulations are needed. The tests of broken rock zones are conducted in the seepage roadway of one mining area in Daye iron mine, after blasting is carried out.

Mining Area Condition
The highest elevation in the testing mining area is +250 m, and the sublevel caving method is used in the mining area. The ore body is the concealed one, and its hanging wall rock is metamorphic diorite. Part of the contacting zone between ore and rock is with frequent geological activity, rock crushing, which belongs to the unstable area type IV and V [11,12]. When constructing the roadway, the shotcrete-bolt support, partial the bolt-shotcrete-mesh support, is usually used, that can maintain the stability of the roadway tunnel. But during a segmented ore body is caving, the surrounding rock of the roadway is affected by the long-time and repeated production blasting vibration, seepage water pressure and ground pressure in the surrounding rock of the roadway will change, the roadway support will         2 Geofluids gradually crack and deform, and some surrounding rocks will collapse, resulting in smashing rock drilling equipment, underground transportation will be blocked, especially part of the ore body cannot be recovered due to long roadway collapsing. In order to guide the roadway support and keep the roadway stable before blasting, it is necessary to test the broken zone of the surrounding rock in different periods in order to detect the gradual damage development of the surrounding rock.

Onsite Determination of Surrounding Rock Broken Zone
At present, there are many methods of measuring the surrounding rock broken zone. These methods can mainly be divided into two categories: electromagnetic wave detection and acoustic detection methods [13][14][15][16][17]. The acoustic detection method is to use the sound wave instrument and other equipment to excite sound (ultrasonic) waves through the rock into the receiving device, by analyzing the change of the acoustic parameters to determine the broken rock range. Because the equipment is simple and the price is relatively cheap, the sound wave method is used to detect the broken zone of the surrounding rock.
3.1. Test System. The RSM-SY5 intelligent acoustic test system is adopted through the refraction wave method in testing. The system is mainly composed of RSM-SY5 acoustic instrument, one emissive probe, two receiving probes (one transmitting double receiving form, as shown in Figure 1), a laptop, a pump, a pipe, wires, and a bracket. As shown in Figure 2, single-hole refraction test steps are as follows: (1) Sent the transceiver test probe to the hole test position with steel pipe and fix it with bracket (2) Inflate the balloon below the hole probe by manual airbag, make it close to the hole wall, tighten the screws, and make the air the airtight gasbag to the bottom hole section form a closed space (3) Pump water into the closed space, until it is full and the probe is completely submerged (4) Keep the probe completely submerged state, and stimulate emission probe to test (5) Select the appropriate stable waveform as a result of the test, record in the disk, and complete a test (6) Loosen the screw, discharged air out of the airtight bag, and retreat the rod back to the next measuring point position according to the test step distance, repeat the above process on to the next test 3.2. Measuring Point Placement. The measuring points are arranged in the 48 th extracting drift at the -50 m level of the mining area, as shown in Figure 3. During the testing process, there is normal production of blasting and falling ore in the surrounding area nearby the measuring points, and the testing results can fully reflect the comprehensive effect of adjacent blasting. Three testing sections are arranged in the 48 th extracting drift, from caving ore blasting pile back to the roadway direction; every testing section consists of 5 testing holes in a row, adjacent testing sections are spaced 5 rows of blasting holes apart, and the row space of the blasting holes is 1.8 m. The total length of the measuring area is about 28 m. The testing holes are arranged in the roof of the 48 th extracting drift, and every hole is about 3 m depth and 50 mm in diameter. The layout is shown in Figure 3.

Test Results and Data
Processing. The relationship between the velocity and hole depth is obtained from the data, as shown in Figure 4. It shows the relationship between the velocity and hole depth of No.1 hole of the 3 rd row. When the 48 th extracting drift formed, the surrounding broken rock zone is 1.28 m, and the rock mass acoustic wave velocity is 4400~5400 m/s. After the No.1 row of blasting holes in the 48 th extracting drift blasted, the surrounding broken rock zone is 1.42 m, and the rock mass acoustic wave velocity is 4000~4800 m/s. After several production blasting, the broken  Figure 7: Distribution of the seepage field in the surrounding rock of the seepage tunnel.
3 Geofluids rock zone is 1.50 m, and the rock mass acoustic wave velocity is 1900~2500 m/s.
In the blasting vibration caused by rock crack initiation and propagation, the test wave velocity is reduced, and the repeated production blasting vibration makes crack continue to expand, its scope is increased, resulting in surrounding broken rock zones increased.

Numerical Simulation of Seepage Damage
The upper part of the rock mass is in the fissure zone nearby ground, and the roadway is in a wet state affected by the surface runoff. Due to the effect of blasting excavation, the surrounding rock of seepage roadway has more internal microcracks. The surrounding rock away from the roadway has good integrity, its internal fissure has less fracture compared with the upper zone and the surrounding broken rock zone, and it becomes a less penetrated layer, as shown in Figure 5

Geofluids
Considering the scene, there is only water dropping from the roadway walls before blasting mining; after a period of production blasting, there are two sites obviously changing with the weather: larger water flow when it rains and less when it is sunny. It is estimated that there is a crack through the earth surface. Using the RFPA2D-Flow system, the calculation model is established [18][19][20][21] as shown in Figure 6. According to the measured data and modeling experience [20][21][22][23], the mechanical parameters of the model are shown in Table 1.
The actual stress condition of the research part of the mining area is as follows: the vertical stress from the above is 7.47 MPa, the horizontal pressure is from 16 MPa to 17.1 MPa from the top down and the periodical falling ore blasting vibration from the upper and horizontal direction. Therefore, the initial condition of the calculation model is set as follows: the vertical stress on the above surface is 7.47 MPa and the horizontal pressure is from 16 MPa to 17.1 MPa compliance with liner change. Upper boundary subjected to blasting vibration is assumed as the loading and unloading process, the loading process from the smallest stress to maximum stress is divided into four step, and the unloading process is divided into four steps too; the minimum is 0 MPa, and the maximum is 20 MPa; the inertia acceleration is 9.8 m/s, direction downward; the under surface constraints by displacement. The penetration pressure from the above is about 3 MPa.

The Seepage Field Evolution in Surrounding
Rock of the Roadway. The seepage roadway elevation is -50 m, and the surface vertical lever is +250 m high, so the osmotic pressure head is 300 m. When the production is not on, the water is infiltrated into the rock mass through the main crack and the branch cracks. So, the upper central has a large crack connected to the surface, and water through the crack penetrates into the rock. Because there is 300 m deep from the surface, the penetration pressure is about 3 MPa.
The internal fissure in the surrounding rock nearby the roadway is developed, and rock permeability is good and water in the fissure seeps down quickly into the roadway, so   Shear stress (Pa) Figure 12: Test diagram of the roadway surrounding broken rock zone. 5 Geofluids the water pressure in the surrounding rock nearby the roadway is almost zero.
Through simulation, it is clear to see the activation, extension progress that the main cracks in the upper part of the surrounding rock of the seepage roadway are in under the action of seepage pressure, ground stress, and cyclic loading. Cracks of the upper and lower surrounding rocks subjected to cyclic loading gradually expand and run through at a certain moment, as shown in Figure 7. The surface water pours into the roadway [24], which is consistent with the fact.

Stress Field Evolution in Surrounding Rock of Seepage
Roadway. After the roadway excavation, surrounding rock is in equilibrium under the ground stress and seepage pressure, the maximum principal stress has a stress concentration phenomenon, and the minimum principal stress under the influence of seepage water pressure exists with local tension stress in some defects, as shown in Figure 8. The broken rock zones can be interpreted from the stress distribution. There is a stress-reducing area nearby the surface of the roadway, as shown in Figure 8(a), and the area within the outer boundary of the stress reduction zone is regarded as a broken rock zone.
The cyclic loading is loaded to simulate the vibration load produced by the blasting production. Under cyclic loading, the stress distribution in the first cycle is shown in Figures 9-11. At the starting point of cyclic loading, the maximum principal stress concentrates in the local area, as shown in Figure 8(a). As the stress increases to the peak of cyclic load, the maximum principal stress increases, the stress concentration becomes more obvious, and the stress at the concentration area also increases, as shown in Figure 9(a). Then, the cycle load gradually drops to the wave valley, the maximum principal stress decreases, and damages appear in the stress concentration area and part of stress release, as shown in Figure 9(b). After the cycle of load up to the end of the cycle, the maximum principal stress increases, the left surrounding rock of the roadway is damaged, and the stress concentration area transfers to the internal area of surrounding rock of the roadway.
During the course of the first cycle, the minimum principal stress decreases with cycle load increasing, and the number of the areas in which the tensile stresses exist has a tendency to increase, as shown in Figures 8(b) and 10(a). Under the effect of seepage water pressure, tensile stress appears around some cracks, which leads to some cracks more easily generating shear failure. The minimum principal stress reaches the minimum value at the wave valley of cycle loads, and the number of the areas in which the tensile stresses exist decreases, as shown in Figure 10(b).
The evolution of shear stress is similar to that of the maximum principal stress, as shown in Figure 11. It can be seen that the surrounding rock of roadways is damaged under repeated cyclic loading, and stress concentration area gradually transferred to the deep surrounding rock.

Extension Evolution of Broken Rock Zones of Seepage
Roadways. The surrounding rock damage occurs under the action of hydraulic pressure, ground stress, and cyclic loadings, and the stress concentration area will transfer to the deep surrounding rock. According to the phenomenon of stress concentration area transferring, the damage state of the surrounding rock can be surveyed and the broken rock zones of surrounding rock can be determined. When the surrounding rock is subjected to the cyclic loading in the ground stress field and the seepage field, the new cracks in the surrounding rock will generate, the local stress will be released, and the concentrated stress area gradually transferred to the deep rock. Therefore, the lower stress area of roadway surrounding rock is the damaged area [25,26]. There is almost no development of the broken zone at the bottom of the roadway, but on the left side and the upper of the roadway, the broken zone is gradually enlarged, because of the main fracture zone, as shown in Figure 11.
According to the onsite measured hole position, measured holes are arranged on the simulate section, as shown in Figure 12. Measuring the lower stress area of the roadway surrounding rock, the broken rock zone is taken on the hole's direction. Similarly, measuring the stress distribution section of first three cycles, the broken rock zone is taken, as shown in Table 2.
The simulation value of the broken rock zones of the roadway surrounding rock is smaller than the measured value, the maximum error is 19.4%, and most of the error is within the allowable range. It shows that the results of the numerical simulation are consistent with that of the test.

Conclusions
(1) Through testing the broken rock zone of surrounding rock in the water seepage roadway of Daye iron mine mining area, it is found that blasting load makes the broken zone gradually expand; the test data obtained can provide the basis for roadway support design (2) Through the numerical simulation on surrounding broken rock zones in the seepage roadway, it is found that the cracks in the surrounding rock gradually initiate, expand, and run through under the effect of seepage pressure, ground stress, and cyclic loading, and the broken scope is enlarged. Under the effect of seepage water pressure, tensile stress appears around some cracks, which leads to the partial cracks more likely to generate shear failure. In the numerical simulation, with repeated cyclic load, part of the surrounding rock is damaged, the stress concentration area gradually transferred to the deep rock, and the range of surrounding broken rock zone of the roadway increases

Data Availability
The data in this article is raw data.

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