Numerical modeling of destress blasting for strata separation

Destress blasting(DB)implemented along the perimeter of safety pillars is a special application of destressing in coal longwall mining.The goal is to separate relatively more deformed mined areas from safety pillars,such as shaft pillars or cross-cut pillars,to reduce the transfer of high stresses to the protective pillar.This case study aims to numerically simulate selected destress blasts in the Czech part of the Upper Silesian Coal Basin and examine its impact on stress transfer to the safety pillar area.To separate the area between the protective pillar and the longwall(LW),two fans of five 93-mm blast holes(length of 93e100 m)were drilled from the gate roads into the overburden strata.Each set of blast holes was fired separately in two stages without time delay.The explosive charge(gelatin-type of explosive)of each stage is 3450 kg.The two DB stages were fired when the longwall face was approximately 158 m and 152 m away from the blast.A 3D mine-wide model is built and validated with in situ stress measured with hydrofracturing.Mining and destressing in three 5-m thick coal seams are simulated in the region.Numerical modeling of DB is successfully conducted using a rock fragmentation factor a of 0.05 and a stress reduction/dissipation factor β of 0.95.Buffering of transfer of additional stress from the mining area into the safety pillar is evaluated by comparison of yielding volume before and after DB.It is shown that yielding volume drops after DB by nearly 80%in the area of the destressing panel and near the safety shaft pillar.

Effect of surface retaining elements on rock stability:laboratory investigation with sand powder 3D printing

This study aims to investigate the benefcial efects of surface retaining elements (SREs) on the mechanical behaviors of bolted rock and roadway stability. 3D printing (3DP) technology is utilized to create rock analogue prismatic specimens for conducting this investigation. Uniaxial compression tests with acoustic emission (AE) and digital image correlation techniques have been conducted on 3DP specimens bolted with diferent SREs. The results demonstrate that the strength and modulus of elasticity of the bolted specimens show a positive correlation with the area of the SRE;the AE characteristics of the bolted specimens are higher than those of the unbolted specimen, but they decrease with an increase in SRE area, thus further improving the integrity of the bolted specimens. The reinforcement efect of SREs on the surrounding rock of roadways is further analyzed using numerical modelling and feld test. The results provide a better understanding of the role of SREs in rock bolting and the optimization of rock bolting design. Furthermore, they verify the feasibility of 3DP for rock analogues in rock mechanics tests.

Methodology for the estimation of expansive cement borehole pressure

This work is part of a multi-phase project which aims to develop a sound methodology for rock fragmen-tation in underground mines using expansive cement.More specifically,it is the first phase of the project which focuses on laboratory tests to investigate the mechanical performance of expansive cement,also known as soundless chemical demolition agents(SCDA).This paper reports the results of laboratory tests conducted on instrumented thick-walled cylinders filled with expansive cement.Expansive pressure evo-lution and temperature variation with time are first examined for different borehole diameters.The clas-sical analytical method for expansive pressure estimation is validated with direct pressure measurement using high-capacity pressure sensor,and an empirical model is obtained.A new methodology based on iterative procedure is developed using axisymmetric finite element modelling and test results to derive the modulus of elasticity of the expansive cement at peak pressure.The results of this study show that the expansive pressure increases with borehole diameter when the rigidity of the steel cylinder is constant reaching 83 MPa for a 38.1 mm borehole.It is also shown that the expansive pressure decreases signif-icantly with increased cylinder rigidity for the same borehole diameter.The newly developed methodol-ogy revealed that the modulus of elasticity of expansive cement at peak pressure is estimated at 8.2 GPa.A discussion on the extension of the findings of this work to hard rock mining applications is presented.

Identifcation of large-scale goaf instability in underground mine using particle swarm optimization and support vector machine

An approach which combines particle swarm optimization and support vector machine(PSO–SVM)is proposed to forecast large-scale goaf instability(LSGI).Firstly,influencing factors of goaf safety are analyzed,and following parameters were selected as evaluation indexes in the LSGI:uniaxial compressive strength(UCS)of rock,elastic modulus(E)of rock,rock quality designation(RQD),area ration of pillar(Sp),the ratio of width to height of the pillar(w/h),depth of ore body(H),volume of goaf(V),dip of ore body(a)and area of goaf(Sg).Then LSGI forecasting model by PSO-SVM was established according to the influencing factors.The performance of hybrid model(PSO+SVM=PSO–SVM)has been compared with the grid search method of support vector machine(GSM–SVM)model.The actual data of 40 goafs are applied to research the forecasting ability of the proposed method,and two cases of underground mine are also validated by the proposed model.The results indicated that the heuristic algorithm of PSO can speed up the SVM parameter optimization search,and the predictive ability of the PSO–SVM model with the RBF kernel function is acceptable and robust,which might hold a high potential to become a useful tool in goaf risky prediction research.

Volumetric analysis of rock mass instability around haulage drifts in underground mines

Haulage networks are vital to underground mining operations as they constitute the arteries through which blasted ore is transported to surface. In the sublevel stoping method and its variations, haulage drifts are excavated in advance near the ore block that will be mined out. Numerical modeling is a technique that is frequently employed to assess the redistribution of mining-induced stresses, and to compare the impact of different stope sequence scenarios on haulage network stability. In this study,typical geological settings in the Canadian Shield were replicated in a numerical model with a steeplydipping tabular orebody striking EW. All other formations trended in the same direction except for two dykes on either side of the orebody with a WNW-ESE strike. Rock mass properties and in situ stress measurements from a case study mine were used to calibrate the model. Drifts and crosscuts were excavated in the footwall and two stope sequence scenarios-a diminishing pillar and a center-out one-were implemented in 24 mining stages. A combined volumetric-numerical analysis was conducted for two active levels by comparing the extent of unstable rock mass at each stage using shear,compressive, and tensile instability criteria. Comparisons were made between the orebody and the host rock, between the footwall and hanging wall, and between the two stope sequence scenarios. It was determined that in general, the center-out option provided a larger volume of instability with the shear criterion when compared to the diminishing pillar one(625,477 m~3 compared to 586,774 m~3 in the orebody; 588 m~3 compared to 403 m~3 in the host rock). However, the reverse was true for tensile(134,298 m~3 compared to 128,834 m~3 in the orebody; 91,347 m~3 compared to 67,655 m~3 in the host rock)instability where the diminishing pillar option had the more voluminous share.

Quantitative analysis with plastic strain indicators to estimate damage induced by fault-slip

In the present study, methodologies to evaluate damage around an underground opening due to seismic waves arising from mining-induced fault-slip are examined. First, expressions for an associated flow rule with a failure criterion are developed for biaxial stress conditions, which are implemented into FLAC3D code. A three-dimensional(3D) mine model encompassing a fault running parallel to a steeply dipping orebody is constructed, whereby static and dynamic analyses are performed to extract stopes and simulate fault-slip in dynamic condition, respectively. In the analysis, the developed biaxial model is applied to the stope wall. The fault-slip simulation is performed, considering shearing of fault surface asperities and resultant stress drop driving the fault-slip. Two methodologies to evaluate damage caused by seismic waves arising from the simulated fault-slip are examined:(i) the ratio of dynamic plastic strain increment to elastic strain limit and(ii) plastic strain energy density. For the former one, two types of strain increments are tested, namely effective shear strain increment and volumetric strain increment.The results indicate that volumetric strain increment is a suitable index for detecting damage near the stope wall, while effective shear strain increment is appropriate for evaluating damage in backfill. The evaluation method with plastic strain energy density is found to be capable of assessing damage accumulated in an extensive area caused by rock mass oscillation due to seismic wave propagation. Possible damage to mine developments in the proximity of a stope is clearly described with the index. The comparison of the two methods clarifies that the former one assesses "instantaneous" damage, which is found to be different from "accumulated" damage calculated using plastic strain energy density, in terms of damage area and its location. It is thus concluded that the combination of the two methodologies leads to more accurate damage assessment as a proper measure against rockburst.

Strategic sill pillar design for reduced hanging wall overbreak in longhole mining

Steeply dipping,vein and tabular orebodies are traditionally extracted with longitudinal retreat mining methods such as Eureka and Avoca in a bottom-up sequence with delayed backfill.To increase productivity,sill pillars in the orebody are used to separate mining zones thus allowing production to take place simultaneously in two or more zones.While such mining methods are productive,they may be accompanied with high volumes of hanging wall overbreak causing significant unplanned ore dilution.In this work,it is shown through a mine design case study of a narrow vein deposit that a sill pillar could also play a significant role in limiting hanging wall overbreak.To demonstrate the role of sill pillar,a novel numerical modelling scheme is proposed to account for progressive stope wall overbreak.A numerical modelling approach of element death and rebirth is developed to allow for the detected stope overbreak to be immediately removed and replaced with backfill material before upper-level stope extraction.It is further shown that the average overbreak volume could be reduced by as much as 33%when the sill pillar is strategically placed in the lower half of a mine plan.

Instantaneous stress release in fault surface asperities during mining-induced fault-slip

Fault-slip taking place in underground mines occasionally causes severe damage to mine openings as a result of strong ground motion induced by seismic waves arising from fault-slip. It is indicated from previous studies that intense seismic waves could be generated with the shock unloading of fault surface asperities during fault-slip. This study investigates the shock unloading with numerical simulation. A three-dimensional （3D） numerical model with idealized asperities is constructed with the help of discrete element code 3DEC. The idealization is conducted to particularly focus on simulating the shock unloading that previous numerical models, which replicate asperity degradation and crack development during the shear behavior of a joint surface in previous studies, fail to capture and simulate. With the numerical model, static and dynamic analyses are carried out to simulate unloading of asperities in the course of fault-slip. The results obtained from the dynamic analysis show that gradual stress release takes place around the center of the asperity tip at a rate of 45 MPa/ms for the base case, while an instantaneous stress release greater than 80 MPa occurs near the periphery of the asperity tip when the contact between the upper and lower asperities is lost. The instantaneous stress release becomes more intense in the vicinity of the asperity tip, causing tensile stress more than 20 MPa. It is deduced that the tensile stress could further increase if the numerical model is discretized more densely and analysis is carried out under stress conditions at a great depth. A model parametric study shows that in-situ stress state has a significant influence on the magnitude of the generated tensile stress. The results imply that the rapid stress release generating extremely high tensile stress on the asperity tip can cause intense seismic waves when it occurs at a great depth.

Determination of mining-induced stresses using diametral rock core deformations

Knowledge of ground stresses is crucial for ground control activities such as the design of underground openings,selec-tion of support systems,and analysis for stability.However,it is a known fact that far field stresses experience changes in orientation and magnitude due to the presence of geological structures and due to the excavations created by mining activi-ties.As a result,in-situ stresses around drifts,ramps,and stopes in underground mines are quite different from far field or pre-mining stresses.The purpose of this research is to develop a simple and practical methodology for determining in-situ stresses.Stress relief occurs once the rock core is drilled off.Such relief is a function of the surrounding stress field.This study uses exploration rock cores that are drilled off for the purpose of orebody definition in the underground mine.The method measures and analyzes the diametral core deformations in laboratory.Two case studies from operating underground mines are presented for demonstration.In these case studies,rock core deformations are measured with a customized test apparatus and rock samples were prepared and tested for Young's modulus and Poisson's ratio.The differential stress,namely the difference between the local principal stresses in the plane perpendicular to the core rock axis is calculated.It is shown that this methodology is useful for determining the brittle shear ratio in the rock mass,which is of primary interest to ground control studies.

Guest Editorial–Special issue on recent advances in mine safety

Mine safety has always been a priority to operating mines.The topic has captured the attention of researchers and scientists over the last two decades,and thanks to the support and commitment of the mining industry,a whole slew of new safety innovations has been developed,particularly in areas like mining automation&robotics,virtual reality,remote sensing,machine learning and others.As a result,remarkable progress has been made in both underground and open pit mining in terms of accident and fatality records,dynamic phenomena,and catastrophic failure of large mining structures.