A Non-Uniform Velocity Model and FLAC/SPECFEM2D Coupled Numerical Simulation of Wave Propagation in Underground Mines

"It is well known that Peak Particle Velocity (PPV) plays a critical role in triggering rockburst damage in burst-prone underground mines. To estimate the PPV distribution around excavations, advanced numerical modeling can be used. In wave propagation numerical modeling, velocity field is often assumed to be uniform around an excavation. However, it is observed that wave velocities near mine openings are not constant, but rather non-uniform. Wave velocity in a rock mass depends on many factors such as rock quality and confinement. In particular, confinement has a large influence on wave velocity. The confinement around an excavation is low at the wall boundary and increases gradually and reaches a constant value far away from the excavation boundary. Hence, stress redistribution due to excavation will have a profound influence on wave propagation in underground mines. A better understanding of the influence of excavation effect on wave propagation helps to interpret the field observation of rockburst damage and improve mine safety. For this purpose, a non-uniform velocity model need to be constructed and the influence excavation effect should be taken into account in numerical wave propagation modeling.In the present study, a non-linear velocity model is developed firstly by considering the influence of confinement on wave velocity. This model links wave velocity to confinement explicitly. A coupled numerical method is then developed and used to study the seismic wave pattern around excavations in underground mines. Two codes, FLAC and SPECFEM2D, are coupled in the analysis. The motivation of developing the FLAC/SPECFEM2D coupled approach is to take merit of each code, such as the stress analysis capability in FLAC and the powerful wave propagation analysis capability in SPECFEM2D. The coupling process is realized through three steps: (1) static stress redistribution around an excavation is analyzed in FLAC; (2) data for coupling analysis, such as non-uniform velocity field (based on the new developed velocity model), mesh, absorbing boundaries, etc., are generated by FISH functions in FLAC and passed to SPECFEM2D; and (3) wave propagation modeling is performed in the SPECFEM2D considering a fault-slip seismic source. Because stress redistribution around the excavation is considered, non-uniform velocity field for the SPECFEM2D model can be obtained and this is a notable feature of this study. Very large difference in the S-wave response in these two modeling approaches with uniform and non-uniform velocity fields can be seen from the modeling results.Using the FLAC/SPECFEM2D coupled modeling approach, accurate velocity models can be constructed and this in turn will help to produce accurate PPV distribution around the excavations. Accurate PPV contours can be used to identify high rockburst risk zones for improved dynamic ground support strategy."