Fractured Rock Mass Research Articles

An Improved Method for Calculating Wave Velocity Fields in Fractured Rock Based on Wave Propagation Probability

Ultrasonic velocity field imaging offers a robust tool for characterizing and analyzing damage and its evolution within fractured rock masses. The combined application of ultrasonic first arrival waves and coda waves can significantly enhance the accuracy and range of velocity field imaging. This manuscript introduces an improved imaging method that integrates the propagation probability distribution of the first arrival and coda waves to calculate the velocity field. The proposed method was applied to the velocity field imaging of a medium with multiple scatterers and varying degrees of fracturing. The overall error and calculation unit error of the proposed method were analyzed, and its improvement in calculation accuracy and applicable scope was verified. Additionally, this method was employed to image the velocity field during the damage process of fractured rock masses. The imaging results were compared against digital speckle patterns to confirm the method’s suitability. Finally, we discussed the impact of measurement errors and sensor missing on the accuracy of the computational outcomes presented in this method. These two situations will affect the calculation results, and the influence of reducing the number of sensors is smaller than that of measuring time shifts with error.

Numerical simulation of blasting behavior of rock mass with cavity under high in-situ stress

With the shift of coal seam mining to the deep, the in-situ stress of coal and rock mass increases gradually. High ground stress can limit the generation of rock cracks caused by blasting, and blasting usually shows different crushing states than low stress conditions. In order to study the blasting expansion rule of rock mass with cavity under high ground stress and the rock mass fracture state under different side stress coefficients. In this paper, the effective range of blasting and the stress distribution under blasting load are analyzed theoretically. The RHT (Riedel-Hiermaier-Thoma) model is used to numerically simulate the blasting process of rock mass with cavity under different ground stress, and the influence of ground stress and lateral pressure coefficient on the crack growth of rock mass is studied. The results show that when there is no ground stress, the damage cracks in rock mass are more concentrated in the horizontal direction and the fracture development tends to the direction where the holes are located, which confirms the guiding effect and stress concentration effect of the holes in rock mass, which helps to promote the crack penetration between the hole and the hole. The length difference of horizontal and vertical damage cracks in rock mass increases with the increase of horizontal and vertical stress difference. Under the same lateral stress coefficient, the larger the horizontal and vertical stress difference is, the stronger the inhibition effect on crack formation is. For blasting of rock mass with high ground stress, the crack formation length between gun holes decreases with the increase of stress level, and the crack extends preferentially in the direction of higher stress. Therefore, the placement of gun holes along the direction of greater stress and the shortening of hole spacing are conducive to the penetration of cracks between gun holes and empty holes. The research can provide reference for rock breaking behavior of deep rock mass blasting.

Study on the failure mechanism of high-temperature granite under two cooling modes

In the geothermal development of hot dry rock (HDR), both the drilling of the wellbore and the heat exchange of the heat reservoir involve the effects of different cold and hot conditions on the high-temperature rock mass. The testing machine for rock mechanics was used to conduct a uniaxial compression test and carry out micro testing on the treated samples; furthermore, with the help of scanning electron microscopy the fracture mechanism of granite subjected to different temperatures and cooling methods was studied. The results show: (1) With the gradual increase in temperature, the compressive strength of granite under the two cooling methods gradually decreases. (2) The failure modes of the samples under the two cooling methods are mainly shear failure of the "Y" type. The degree of damage of the sample under water cooling is significantly greater than that under natural cooling. Electron micrographs could confirm these results. (3) It can be obtained by testing the mineral composition and element changes of granite at different temperatures. When the temperature reaches 600℃, its change is more pronounced. The results of this study can provide a theoretical reference for the failure of the wellbore and the degree of fracture of the thermal reservoir rock mass during geothermal development.

Investigation on the effect of the relative position between fracture and injection-production well on thermal extraction performance

The development of hot dry rock (HDR) resources is mainly through hydraulic fracturing to build fracture in low permeability high-temperature rock to form working fluid flow channels. The control mechanism of thermal transfer and conversion process in fracture system is the key to influence thermal extraction effect. In this paper, a temperature-hydraulic-mechanical (THM) coupling model of fractured rock mass is established considering the temperature variation of rock elastic modulus and Poisson's ratio. The dynamic changes of fracture aperture and permeability are investigated. Taking fractured HDR reservoir as a case study, the depletion thermal extraction process of dual-well enhanced geothermal system (EGS) is simulated, the spatial-temporal evolution of multi-field are analyzed. The results show that the cold front of the horizontal and vertical single fracture model advance in an oval shape and a rectangle shape, respectively. Over 30 years, the horizontal fracture aperture typically changes by a factor ranging from 1.02 to 1.12, permeability typically changes by a factor ranging from 5 to 15. The vertical fracture aperture typically changes by a factor ranging from 1.02 to 1.14, permeability typically changes by a factor ranging from 10 to 25.

The influence of cracks on pillar strength based on SRM and DFN models

This research focuses on the examination of natural fractures within underground mines, emphasizing their substantial impact on the strength and stability of ore pillars. The study adopts the Strength Reduction Method (SRM) theory and employs the Discrete Fracture Network (DFN) model, offering a novel approach to investigating the behavior of fractured rock masses. The objective of this article is to analyze the influence of natural fractures on the strength of ore pillars by employing SRM and DFN methods. The research begins by establishing a multi-level amplification program that incorporates a homogenization process. The findings reveal that, for a W/H ratio of 0.5, the strength reduction aligns consistently with empirical equations. A notable observation is that when W/H is less than or equal to 1.0, there is good agreement, but when W/H exceeds 1.0, there is a tendency to overestimate pillar strength. Subsequent investigations emphasize the significance of considering pillar development in the overall assessment of pillar forces. The study underscores the importance of integrating pillar development into the analysis, aligning with previously established research results. Therefore, by evaluating the strength and failure mechanism of columns under different aspect ratios, we studied the influence of discrete discontinuous bodies on column stability, revealed the influence of natural cracks on column strength, and provided theoretical basis and reference for the design and support of underground mines.

Combined effects of the roughness, aperture, and fractal features on the equivalent permeability and nonlinear flow behavior of rock fracture networks

The hydraulic properties of a fractured rock mass are largely controlled by connected fracture networks. A thorough understanding of the physical flow processes in fracture networks is essential for assessing the transport capacity of the rock mass. However, the fracture surface roughness morphology, fracture distribution characteristics, and fluid flow regimes strongly influence the flow capacity of a fracture network. To this end, the rough topographic characteristics of fracture surfaces were quantified using fractal theory, and then the effective permeability model and nonlinear seepage effect assessment model of the rough fracture network for different flow regimes were developed based on the possible occurrence of laminar and turbulent flows in a single fracture. Finally, the influences of the geometric parameters of the fracture network on the effective permeability and nonlinear flow characteristics were analyzed. The results show that the prediction results of the proposed models are in good agreement with the field test data and can effectively reveal the seepage influence mechanisms under different flow regimes. Additionally, the results show that the effective permeability is closely related to the fractal dimension, relative roughness, aperture scale, distribution characteristics, and hydraulic gradient of the fractures. The nonlinear behavior of fluid flow significantly reduces the effective permeability of the rock mass. The proposed models can provide a reference for evaluating the transport capacity of rock masses under different fracture distributions and flow regimes.

Graph theoretical analysis of limestone fracture network damage patterns based on uniaxial compression test

The topological attributes of fracture networks in limestone, subject to intense hydrodynamics and intricate geological discontinuities, substantially influence the mechanical and hydraulic characteristics of the rock mass. The dynamical evolution of fracture networks under stress is crucial for unveiling the interaction patterns among fractures. However, existing models are undirected graphs focused on stationary topology, which need optimization to depict fractures' dynamic development and rupture process. To compensate for the time and destruction terms, we propose the damage network model, which defines the physical interpretation of fractures through the ternary motif. We focus primarily on the evolution of node types, topological attributes, and motifs of the fracture network in limestone under uniaxial stress. Observations expose the varying behavior of the nodes' self-dynamics and neighbors' adjacent dynamics in the fracture network. This approach elucidates the impact of micro-crack behaviors on large brittle shear fractures from a topological perspective and further subdivides the progressive failure stage into four distinct phases (isolated crack growth phase, crack splay phase, damage coalescence phase, and mechanical failure phase) based on the significance profile of the motif. Regression analysis reveals a positive linear and negative power correlation between fracture network density and branch number to the rock damage resistance, respectively. The damage network model introduces a novel methodology for depicting the interaction of two-dimensional (2D) projected fractures, considering the dynamic spatiotemporal development characteristics and fracture geometric variation. It helps dynamically characterize properties such as connectivity, permeability, and damage factors while comprehensively assessing damage in rock mass fracture networks.

A Hybrid Prediction Model for Rock Reservoir Bank Slope Deformation Considering Fractured Rock Mass Parameters

Deformation monitoring data provide a direct representation of the structural behavior of reservoir bank rock slopes, and accurate deformation prediction is pivotal for slope safety monitoring and disaster warning. Among various deformation prediction models, hybrid models that integrate field monitoring data and numerical simulations stand out due to their well-defined physical and mechanical concepts, and their ability to make effective predictions with limited monitoring data. The predictive accuracy of hybrid models is closely tied to the precise determination of rock mass mechanical parameters in structural numerical simulations. However, rock masses in rock slopes are characterized by intersecting geological structural planes, resulting in reduced strength and the creation of multiple fracture flow channels. These factors contribute to the heterogeneous, anisotropic, and size-dependent properties of the macroscopic deformation parameters of the rock mass, influenced by the coupling of seepage and stress. To improve the predictive accuracy of the hybrid model, this study introduces the theory of equivalent continuous media. It proposes a method for determining the equivalent deformation parameters of fractured rock mass considering the coupling of seepage and stress. This method, based on a discrete fracture network (DFN) model, is integrated into the hybrid prediction model for rock slope deformation. Engineering case studies demonstrate that this approach achieves a high level of prediction accuracy and holds significant practical value.

Finite–Discrete Element Method Simulation Study on Development of Water-Conducting Fractures in Fault-Bearing Roof under Repeated Mining of Extra-Thick Coal Seams

The formation of water-conducting fractures in overlying strata caused by underground coal mining not only leads to roof water inrush disasters, but also water-conducting fractures penetrate the aquifer, resulting in the occurrence of a mine-water-inrush disaster and the loss of water resources. It destroys the sustainability of surface water and underground aquifers. This phenomenon is particularly significant in extra-thick coal seams and fault-bearing areas. Numerical simulation is an effective method to predict the failure range of mining overburden rock with low cost and high efficiency. The key to its accuracy lies in a reasonable constitutive model and simulation program. In this study, considering that the three parts of penetrating cracks, non-penetrating cracks, and intact rock blocks are often formed after rock failure, the contact state criterion and shear friction relationship of discrete rock blocks and the mixed fracture displacement–damage–load relationship are established, respectively. Combined with the Mohr–Coulomb criterion, the constitutive model of mining rock mass deformation–discrete block motion and interaction is formed. On this basis, according to the engineering geological conditions of Yushupo Coal Mine, a numerical model for the development of water-conducting cracks in the roof with faults under repeated mining of extra-thick coal seams is established. The results show the following: The constitutive relation of the continuous deformation–discrete block interaction of overlying strata and the corresponding finite element–discrete element FDEM numerical program and VUSDFLD multi-coal seam continuous mining subroutine can numerically realize the formation process of faults and water flowing fractures in overlying strata under continuous mining of extra-thick multi-coal seams. The toughness of sand mudstone is low, and the fracture will be further developed under the repeated disturbance of multi-thick coal seam mining. Finally, it is stabilized at 216–226 m, and the ratio of fracture height to mining thickness is 14.1. When the working face advances to the fault, the stress concentration occurs in the fault and its overlying rock, which leads to the local fracture of the roof rock mass and the formation of cracks. The fault group makes this phenomenon more obvious. The results have been preliminarily applied and tested in Ningwu mining area, which provides theoretical support for further development of roof water disaster control under the condition of an extra-thick coal seam and avoids the loss of water resources in surface water and underground aquifers.

Dynamic response characteristics and damage calculation method of fractured rock mass under blasting disturbance

The mechanical responses and damage characteristics of fractured rock masses under dynamic loads play a significant role in ensuring the safety of blasting operations. The study employs dynamic finite element method (FEM) to simulate blasting in fractured rock masses, revealing the propagation features of blasting stress waves, the evolution of effective stress and the mechanism of damage evolution around the borehole. Based on the Weibull distribution, a rock damage model was established, which yielded the stress-strain relationship and damage equation for rocks under uniaxial impact. The parameters of the damage equation were fitted through numerical SHPB tests. The results indicate that in the direction of higher initial stress, the stress wave propagates faster, and the development of damage surfaces is more pronounced. In cylindrical charge blasting, the maximum blasting stress occurs in the middle section of the borehole, with limited influence from the charge on the bottom of the hole. Within a range of approximately 2 to 3 times the borehole diameter, the surrounding rock sustains complete damage. While within a range of about 4 to 5 times the borehole diameter, the degree of damage to the surrounding rock decreases rapidly, and beyond this range, the damage to the surrounding rock gradually diminishes.

Numerical research on disastrous mechanism of seepage instability of karst collapse column considering variable mass effect

In order to reveal the disastrous mechanism of seepage instability of karst collapse column considering variable mass effect, a variable mass fluid–solid coupling mechanical model of water inrush is established, by considering the random distribution characteristics of a collapse column. Taking Qianjin coal mine as the research background, based on the Weibull distribution theory, the heterogeneous distribution characteristics of rock mass is described, and COMSOL Multiphysics numerical simulation software is employed to simulate the seepage characteristics and inrush water changes in collapse columns under different conditions of homogeneity, water pressure, and initial porosity. The research results show that the greater the homogeneity is, the more water conduction channels are formed, and the porosity increases accordingly, when considering the influence of different homogeneity on the seepage characteristics of broken rock mass, which eventually leads to water inrush accidents and a sharp increase in water inflow. Besides, when studying the seepage evolution law of different water pressures on a broken rock mass, an elevation of water pressure dramatically increases the porosity and seepage rate of the water. Over time, the broken rock particles gradually migrate and the fine particles are transported and eroded by the water flow, resulting in changes in the seepage characteristics and the formation of potential water diversion channels. Finally, when taking into account the effect of different initial porosity on the fractured rock mass seepage characteristics, the greater the original porosity is, the higher the seepage velocity is, and the particle migration increases the permeability. This leads to a more pronounced conductive water passage formation, which reveals the disastrous mechanism of seepage instability of karst collapse column considering variable mass effect.

Investigation on uniaxial compression and fracture damage mode of prefabricated parallel double-jointed red sandstone.

As a special type of joint fracture, the fracture evolution characteristics of parallel double joints have important engineering significance for the stability analysis of fractured rock mass. In this work, a new method for calculating stress intensity factor of parallel double-jointed fractures was importantly proposed. Physical uniaxial compression tests were carried out on parallel double jointed red sandstone filled with cement mortar under different geometric parameters, and the macroscopic mechanical properties and failure characteristics of red sandstone are deeply analyzed. The results show that the larger the connectivity rate is, the smaller the peak stress and strain are. The increase of connectivity rate will affect the change rate of transverse strain in the center of rock bridge. The closer the dip angle of the joint is, the lower the peak stress is and the shorter the failure time is. The damage mode of joint tip encroachment affects the lateral displacement of the rock bridge center, and the displacement is always close to the first damage section. The closer the joint tip is to the load, the easier the end-face penetrating cracks occur. The research content can provide basic support for guaranteeing the stability of underground engineering rock mass.

Use of Mohr diagrams to predict fracturing in rock masses, with applications for predicting sub-surface behavior

Mohr diagrams are a simple and effective method that can help geoscientists consider, model, and predict the ranges of mechanical properties of rock, stresses, fluid pressures, and the resultant fractures that are likely to occur in the sub-surface. Mohr diagrams can be used to make predictions about how rocks may respond to change, with a transition from a stable state to fracturing occurring if there are changes in (1) the failure envelope, (2) stresses, and/or (3) fluid pressure. This article uses Mohr diagrams to address two questions of significance to the energy transition. First, how will metasedimentary rocks, which are potential geothermal reservoir rocks, respond to thermal stimulation? Second, will fractures that may influence the underground storage of radioactive waste develop in a clay sequence during exhumation? Mohr diagrams are shown to be useful for highlighting misconceptions and input data problems, leading to improved understanding of how structures develop.

A complex-network-based estimation of the representative elementary volume and equivalent permeability coefficient for fractured rock masses

The determination of the representative elementary volume (REV) and equivalent permeability coefficient (EPC) of fractured rock masses is crucial for studying the hydraulic behaviors of engineering rock masses. The complex network originated from graph theory and represents a rapidly developing interdisciplinary field in the 21st century. In this paper, we propose a novel method that utilizes complex-network-based metrics to estimate the REV and EPC of fractured rock masses. In this method, fractures and their intersections are treated as nodes and edges of a complex network, respectively. We propose a transformation procedure of a discrete fracture network (DFN) into a complex network, and the corresponding algorithms are conducted in MATLAB. Subsequently, we obtain a complex network, and several complex-network-based metrics, such as average degree, average clustering coefficient, and graph density, are calculated using Gephi software. We employ the average degree to estimate the REV size of DFN models, while the coefficient of variation is used to quantitatively assess changes in the average degree of multiple random models. Furthermore, we establish an empirical formula utilizing the average degree, average clustering coefficient, and graph density to estimate the EPC of a large-scale DFN model based on the EPC of a small-scale DFN model. To verify the developed method, we conduct numerical experiments using the MAFIC (Matrix/Fracture Interaction Code) solver within the FracMan software. This enabled us to calculate the EPC, which was then used to determine the REV. The results demonstrate that: (a) the REV size estimated by the average degree aligns with the REV size obtained through numerical simulations for both isotropic and anisotropic DFN models, and (b) the error rates between the average EPC calculated by the empirical formula and that obtained by numerical simulations are within 10%. Thus, the proposed seepage parameter estimation methods are deemed valid.

Experimental and numerical study on the damage evolution and acoustic emission multi-parameter responses of single flaw sandstone under uniaxial compression

Cracks are generally the weakest position in engineering rock mass. It is of great practical significance to study the influence of cracks on the fracture process of rock mass and the precursor criterion of failure for disaster prevention and mitigation projects. In this paper, particle flow program (PFC2D) is used to analyze the mechanical strength and damage evolution of fractured rock. The multi-parameter stage characteristics of acoustic emission (AE) and the difference of fracture precursor information of fractured rock mass under uniaxial compression (UC) are explored. The results are as follows: (1) The strength and modulus of elasticity of fractured rock samples are smaller than that of intact rock mass, and increase in quadratic function with the increase of fracture angle, and the strength of B-90 sample is closest to intact rock samples; (2) The PFC analysis shows that the cracks in the defective rock samples all originate from the tip of the prefabricated crack, and the crack initiation angle is basically in the range of 60°∼90°. The loading process mainly produces tensile damage, and the number of tensile cracks are far more than that of shear cracks. Compressive stress and tensile stress increase with the increase of load, but the force chain of failure surface is the smallest after complete failure of rock; (3) The acoustic emission ringing number, RA/AF value, b value and entropy value of sandstone have obvious stage characteristics, the main frequency signal is mainly low frequency signal, and large-scale cracks are mainly produced in the failure process. The sudden increase point of AE ringing times and RA/AF value accumulation curve can be used as critical failure point. The fluctuation of B value and entropy value reflects the instability of fracture development, and “cliff-like” descending point and ascending point can be regarded as the precursor of critical failure of rock; (4) The failure precursor of acoustic emission parameters of sandstone appears in the peak load period of 85 %–98 %. The order of response time of different parameters is: cumulative RA/AF value, cumulative event number, entropy value, B value. The response time of B-45 and B-15 is earlier, while that of B-60 and B-30 is later. The research results are expected to deepen the understanding of the relevant standards of damage evolution and failure precursors of defective rocks in engineering applications.

Influence of inertial and centrifugal forces on rate and flow patterns in natural fracture networks

Fluid production from fractured rock masses readily induces fracture flow velocities of meters per second. Yet, most discrete fracture flow models treat flow as laminar creeping flow or account for inertia effects only by single-fracture constitutive relationships.This numeric simulation study investigates water flow patterns and spatial velocity variations in a natural fracture network with mm-wide open fractures, studying the transition from laminar creeping to turbulent flow. After verification with a fracture intersection model, a Reynolds-time-averaged Navier Stokes solver serves to analyse flow regimes and velocity distribution. Our results show that for fracture flow velocities greater than ∼1-cm/s, fluid inertia begins to markedly alter flow patterns and the overall velocity distribution in the network. The pressure-gradient-flow relationship therefore becomes non-linear long before the flow in straight fractures enters the weak inertia regime. This prominence of inertia effects highlights the need to improve fracture network flow models.

An Improved Failure Model Considering the Arching Effect for Tunnel Face Stability in a Weak and Fractured Rock Mass

An Improved Failure Model Considering the Arching Effect for Tunnel Face Stability in a Weak and Fractured Rock Mass

Sulfoaluminate cement-based double-liquid material grouting in fractured rock

Grouting is a very effective means to repair modified fractured rock. However, problems such as severe grout leakage, slow increase in stone body strength, uncontrolled solidification time and poor permeability are often encountered in the process of traditional grouting reinforcement of fractured surrounding rock. Sulfoaluminate cement-based double-liquid material has the characteristics of quick setting, early strength, micro-expansion and high strength, which shows good applicability in grouting reinforcement of fractured surrounding rock. Currently, there are few reports on this material. Therefore, this study mainly introduces the performance of sulfoaluminate cement-based grouting material, grouting technical process and engineering application case. The sulfoaluminate cement-based double-liquid material grouting technology proposed in this study will provide a new treatment method for the reinforcement of fractured rock masses in underground engineering.

Particle Flow Simulation Study of Damage Evolution in Expansive Slurry–Fractured Rock Mass Composites under Direct Shear Conditions

Particle Flow Simulation Study of Damage Evolution in Expansive Slurry–Fractured Rock Mass Composites under Direct Shear Conditions

Applying Novel Micro–Nano-Siliceous Grouting Material Improves Mechanical Behaviors of Fractured Rock Masses and Microscopic Characteristics of Slurry–Rock Interface: Experimental and Numerical Simulation Investigations

Applying Novel Micro–Nano-Siliceous Grouting Material Improves Mechanical Behaviors of Fractured Rock Masses and Microscopic Characteristics of Slurry–Rock Interface: Experimental and Numerical Simulation Investigations