International Journal of Rock Mechanics and Mining Sciences最新文献

This study investigates the rock mechanics and anisotropic properties of the Montney Formation, Alberta, through two sets of experiments: unconfined compressive strength tests and triaxial compression tests, supplemented by ultrasonic wave velocity measurements. These tests enabled the calculation of dynamic and static stiffness properties and Thomsen anisotropy parameters (ε, δ, γ). Our findings reveal that the Montney Formation exhibits weak anisotropy, with ε values generally below 0.10, contrasting with other strongly anisotropic shale formations. Specimens exhibiting higher clay content and total organic carbon levels show more pronounced anisotropy. Static measurements exhibit a higher degree of anisotropy compared to dynamic measurements. The investigation also explored loading and unloading stiffness parameters, noting a higher E loading-to-unloading ratio for rock specimens with lower elastic properties. This provides important information on the rock's response to stress path effects during stimulation and exploitation. The study also discusses other rock mechanics parameters including Poisson's ratio, tensile strength, and brittleness. The research contributes to a more comprehensive understanding of the geomechanical and petrophysical behavior of the Montney Formation, aiding reservoir characterization and hydrocarbon exploration and production.

Gas seepage and progressive failure of coal are common energy-driven mining phenomena. A comprehensive understanding of the energy-driving mechanism behind the catastrophic behavior of mining-induced coal is fundamental to innovating the technology of coal and gas co-mining. Thus, this study simulated three typical mining stress evolution process in protective coal-seam mining (PCM), top-coal caving mining (TCM), and non-pillar mining (NM) to investigate the energy evolution and distribution patterns of coal. The results indicate a strong correlation between energy dissipation and gas seepage. By transitioning from PCM and TCM to NM, the peak elastic strain energy of gas-bearing coal increased by 155.92 %, and the ratio of peak dissipative energy decreased from 51 % to 41 %. Under the PCM stress path, gas seepage decreased the energy storage by 13.52 %, whereas the pre-mining pressure relief and enhanced permeability simulation increased in peak dissipation energy by 49.66 %. Using the cumulative dissipative energy as a damage variable reveals that the degree of coal damage evolution under PCM is higher than other mining methods. Based on the energy-driven damage mechanism, a new coal permeability model was established, and its comparison with classical permeability model demonstrated its excellent fitting effectiveness.

Rock properties are estimated using objective lab/in-situ testing and subjective judgements invoking different types of uncertainties, i.e., aleatory, and epistemic, along them, often indicated by varying information levels. This study presents a unified reliability method to integrate the spatial variability of inputs modelled via alternate uncertainty models (intervals and probability-boxes (p-boxes)) with those modelled as stochastic variables via probability distributions. The methodology employs advanced Karhunen–Loève decomposition to generate interval and random fields of inputs modelled via alternate and stochastic models, respectively. Input properties are allocated to the zones of the numerical model in Fast Lagrangian Analysis of Continua-2D (FLAC-2D) based on their spatial dependency and correlation functions through a developed MATLAB-FLAC coupled code. The methodology is demonstrated for a deep tunnel in Canada to be constructed along a massive rock prone to brittle failures. Intact rock properties are modelled as stochastic variables due to objective estimation, while Geological Strength Index (GSI) and deformation modulus are modelled using alternate models (interval and p-box, respectively) due to subjective and hybrid estimation (double uncertainty propagation algorithm). The results of the methodology are compared with those of traditional deterministic and random field methods. The methodology reduces the subjectivity invoked by including unavailable additional information (e.g., assuming probability distributions of inputs based on literature) and propagates the originally available information of inputs accurately. The final outputs are the p-boxes of response parameters (i.e., displacements and damage zone) instead of their fixed values (deterministic analysis) and probability distributions (traditional reliability analysis), indicating the propagation of both impreciseness and variability of inputs by the method. For this case study, the p-boxes of outputs were bounding their values/distributions estimated via traditional analyses, verifying the accuracy of the methodology. Further, the impreciseness in the outputs, highest in the damage zone extent, was due to imprecision in the estimated GSI.

Stope structural parameters, which are human-controllable, directly impact the safety and economics performance of underground mineral extraction. Current stope design still relies heavily on empirical methods such as stability graphs, due to the complex nature of rock masses and varied stope failure mechanisms. This study aims to enhance stability graphs with machine learning techniques. Firstly, a dataset of 980 records from unsupported stopes was compiled, representing perhaps the largest dataset of its kind so far in the literature. This was achieved through extensive literature review and the collation of an additional 289 records from Chinese mines which were previously not included. An analysis of data reveals that over 90 % of the records fall within a stability coefficient of 0–100 and a hydraulic radius of 0–20 m. Secondly, a stability graphs optimization process was established using Python, eliminating the subjectivity of partitioning. Nine supervised machine learning algorithms were employed and trained to test their performance in partitioning form and predicting accuracy. It was found that the neural network algorithm demonstrated the best overall performance. At last, a neural network with the Keras framework was used to establish a new multilayer perceptron model to generate safety factor probability curves, which were then used to construct the stability graph. To facilitate practical use, mathematical functions fitting safety factor curves within the unstable zone were further formulated. Compared with other empirical stability graphs, our new approach allows designers to more efficiently and reliably select safety factors to determine the stability state according to site specific conditions and technical support systems, thereby providing enhanced guidance for stope design.

This study investigated the crack propagation processes and established the lifetime prediction methods of the concrete-rock interfaces under constant loading. Firstly, the constant loading tests were performed under three-point bending loading on the composite concrete-rock beams with three kinds of interfaces, i.e. natural, 4 × 4, and 7 × 7 interfaces, and under three constant load levels, i.e. initial cracking load, 80 % and 97 % of the peak loads. Then, the clip gauge method and digital image correlation method were employed to detect the crack lengths under constant loading. Finally, the effects of constant load levels on fracture properties of the concrete-rock interfaces were analyzed and the crack propagation processes under constant loading were discussed. The results indicated that the clip gauge method was an available and accurate method to detect the crack lengths of the concrete-rock interfaces under constant loading. Crack propagation processes and crack mouth opening processes under constant loading exhibited the three-stage feature, i.e. deceleration stage, uniform stage and acceleration stage. In addition, the logarithms of the constant loading lifetimes showed an approximate linear relationship with the logarithms of the crack mouth uniform opening rates, and a quantitative relationship was established by using the linear regression analysis to establish the lifetime prediction methods. This study proposed a practical method to predict the service life of concrete structure-bedrock interface. When the crack mouth uniform opening stage was recognized, the lifetime of concrete structure-bedrock interface in service can be predicted by substituting the crack mouth uniform opening rate into the pre-established prediction model. This study can provide the theoretical support for the safety assessment of the concrete structure-bedrock interface in service.

The fracture extension mechanisms and proppant transport characteristics play key roles for optimizing hydraulic fracturing in unconventional reservoirs. In this work, the visual physical simulation method was used to analyze the 3D dynamic extension processes of multi-stage hydraulic fractures and the subsequent impact on proppant distribution. The interfered fracture in multi-stage hydraulic fracturing exhibit alternated radial-circumferential extension behavior. More specifically, these behavior changes from unilaterally radial initiation, circumferential extension to radial extension. Unilaterally radial initiation results in the formation of both mirror and conchoidal features, while circumferential extension leads to stepped feature, and radial extension gives rise to feather feature. The four features of hydraulic fractures result in four types of non-uniform proppant distribution: uniform, scattered, cluster, and regional distribution. The proppant transport shifted from linear to radial mode, promoting further fracture extension. The effects of segment spacing and proppant injection sequence were studied. The results showed that the influence on the interfered fracture decreases as the segment spacing increases. In addition, the propped fracture area is larger when the small-size proppant is injected first, followed by large-size one. The research can improve understanding and optimization of multi-stage hydraulic fracturing in unconventional oil and gas reservoirs.

The trace length as an important indicator reflecting the size of discontinuity, can usually be estimated by setting a sampling window on the surface of the rock mass. However, in high-steep rock slopes with cliffy topography and uneven terrain, traditional methods are difficult to select suitable plane sampling windows. Therefore, estimating the average trace length accurately in this situation has become an urgent problem to be solved. This study proposes a distribution-free method for estimating the average trace length of discontinuities within large windows in complex high-steep slopes. By introducing a cuboid sampling window, the estimation method has been extended to three-dimensional for the first time. Additionally, the proposed method addresses the limitations of existing methods, making it applicable for larger sampling areas. In addition, the new method considers the proportion of three types of traces intersecting the sampling area and their angles with the sampling area, and corrects the new method using the weight of visible traces at both ends, eliminating the impact of sampling biases such as censoring bias, orientation bias, and size bias. The reliability of the proposed method was validated by generating nine sets of simulated trace data with different angles and distribution types. Finally, the new method is applied to a complex high-steep slope with an elevation difference of nearly 600 m on the southeastern edge of the Qinghai Tibet Plateau. The results indicate that for the estimation of trace length in high-steep slope, the new method has smaller errors and lower error volatility. When selecting an appropriate sampling area size, the error of the new method is less than 6 %, and the error fluctuation is less than 4 %. And this study is of great significance for evaluating the stability of high-steep rock mass structures in major projects in hard mountain areas.

Flame jet-assisted mechanical rock drilling is expected to solve the hard-rock crushing challenges in underground engineering such as mining, tunneling, and drilling. Therefore, it is important to investigate mechanical behavior of the rock after flame jet-water cooling treatment. In this paper, granite blocks were subjected to flame jet-water cooling treatment, which means the rapid heating of the granite using a flame with a temperature of 1564 °C, followed by water cooling. Then, physico-mechanical behaviors of samples at different distances (all the distances mentioned later mean the vertical distance between the sample axis and the flame jet-water cooling path) were measured. Results showed physico-mechanical parameters of the sample increased as the distance increased. The V_p (P-wave velocity), UCS (uniaxial compressive strength) and BTS (Brazilian tensile strength) of the samples increased from 3171.7 m/s, 84.84 MPa and 5.01 MPa at 0 mm to 3619.3 m/s, 145.86 MPa and 7.30 MPa at 70 mm, respectively. Meanwhile, physico-mechanical parameters of samples at 70 mm all reached more than 95 % of those of untreated samples. As the distance increased, cumulative AE counts at peak stress gradually enhanced, and failure modes of uniaxial samples shifted from shear to tensile splitting failure. The SEM (scanning electron microscope) images illustrated that the closer to the flame jet-water cooling path, the more thermal cracks produced in the sample. Meanwhile, the damage factor can well characterize the thermal damage degree of granite. According to the variation of physico-mechanical parameters and SEM images of samples with distances, it can be assumed that the thermal damage width of granite after flame jet-water cooling was about 140 mm. It was much greater than the 50 mm width of spalling pits, which greatly increased the feasibility of thermal-assisted mechanical rock drilling. The tensile stress induced by the temperature gradient below the spalling pit is the fundamental cause of thermal damage to the rock, which is further aggravated by water cooling. The findings can provide some guidance for flame jet-assisted mechanical rock drilling.

In this study, rock-like model blast tests with varying toe burdens are conducted and the evolution of blast-induced fractures is analyzed using the digital image correlation (DIC) method. In addition, post-blast fragmentations are image-processed to quantitatively investigate the blasting performance. A three-dimensional numerical model is developed and validated against the experimental results. The effects of blasthole inclination (70°, 80°, and 90°) and short delays on fragment size distribution (FSD) are numerically investigated. Field trials are conducted to verify the numerical results. The experimental results reveal that the overall fragment size increases significantly in proportion to the toe burden. The explosion energy utilization rates range from 3.56 % to 13.16 %, with a tendency to initially increase and subsequently decrease. The numerical results demonstrate that fragments are more evenly distributed with a blasthole inclination of 70°, resulting in a remarkable improvement in rock fragmentation compared to a vertical blasthole. Compared to the toe burden, the influence of short delays on concrete fragmentation is insignificant. The field test results indicate that a narrower FSD range and stable pit walls are achieved by utilizing inclined blastholes in bench blasting. This study provides theoretical guidance for optimizing bench blasting parameters, which has practical significance for improving explosive energy utilization and production efficiency in open-pit mines.

Room-and-pillar mining is one of the oldest and most widely employed underground mining methods worldwide, while in recent years it has also been used for civil engineering projects and applications. The pillar design process is quite straightforward, requiring the assessment of two key elements, namely the pillar's strength and the anticipated loading imposed on it. Lots of research work has been elaborated for the estimation of pillars' strength and the main parameters associated with it, while, on the other hand, in terms of the loading the Tributary Area Theory (TAT) has been traditionally and successfully used as the dominant methodology to assess the post-mining pillar induced stresses. This paper focuses primarily on decoding the pillar loading regime by analyzing the identifying the stresses imposed to them through the use of modern computational tools offered by 2D and 3D numerical analyses, over a series of pillars' type (square, ribs), configurations and initial virgin stress fields. This allows for the development of a benchmarking framework that showcase the differences in the stress conditions and reveal the over-conservative behavior of TAT. More importantly though, it allowed for the expression of two (2) general analytical formulae for the direct estimation of the average vertical elastic stresses on pillars, both for the case of rib & square pillars layout; the most typical patterns applied in this method. They aim at providing an easy to use and more accurate estimation of the pillar stress - as compared to TAT - in alignment with the results obtained from numerical analyses. Thus, with the proposed formulas, the time-consuming numerical process needed for preliminary pillars' dimensioning is tackled in a swift and effective manner.