Rock Mechanical Properties Research Articles

The Changes to the Rock Physical and Mechanical Properties Near the Hole Produced by Successive Impacts of a Steel-Particle Water Jet: A Case Study in a Tight Sandstone

The Changes to the Rock Physical and Mechanical Properties Near the Hole Produced by Successive Impacts of a Steel-Particle Water Jet: A Case Study in a Tight Sandstone

Guidelines for Natural Stone Products in Connection with European Standards.

The selection of ornamental stones for specific applications requires technical guidance since it should be based on the durability, service life, and aesthetic value of the stones. In most cases, these fundamentals provide quantitative data on the usability and performance of ornamental stones. The present study attempts to put forward a quantitative classification system for natural stone products concerning critical rock properties. For this purpose, fundamental physical and mechanical rock properties are listed based on European standards. Then, minimum limit values are proposed for different applications of natural stone products based on retrospective analyses of numerous ornamental stone applications. The suggested limit values based on several physical and mechanical rock properties can guide relevant engineers to initially consider possible rock types for use as natural stones in a wide range of applications. In this context, it is believed that the present study contributes to the natural stone industry by discussing the minimum limit values for the consideration of a wide range of rock types possibly usable in the dimension stone industry.

Fractal analysis of rock pore structure: A comparative study of the Matlab-based image method and the mercury intrusion method

As a naturally porous medium, rocks have a very strong heterogeneity in their internal pore structure, which has an important impact on the mechanical and chemical properties of rocks, etc. Therefore, it is very important to quantitatively characterise the heterogeneity of the pore structure of rocks. Fractal dimension has long been recognised as an effective means of characterising the heterogeneity of the pore structure of porous media and is widely used in oil and gas exploitation, construction materials, mining, and water engineering. There are different methods to calculate the fractal dimension. To verify the consistency between different methods, this paper compares the results of fractal analysis using the mercury intrusion method and MATLAB image fractal analysis using three rock samples with large differences in porosity and calculates the fractal dimension in three different ways on the basis of the mercury intrusion method. The results demonstrate that the fractal dimension of the mercury intrusion method and the box counting dimension of the image analysis obtained by the three methods, although slightly different in numerical value, are consistent in their numerical relationship, i.e., they all conform to the rule that the stronger the non-homogeneity of the pore throat, the larger the fractal dimension. The results of this paper show that fractal dimension is indeed an effective means of characterising rock homogeneity.

Appraisal of rock dynamic, physical, and mechanical properties and forecasting shear wave velocity using machine learning and statistical methods

Direct determination of shear wave velocity requires time, cost, and high accuracy due to the complexity of the rock texture. In current research statistical and intelligent approaches have been used to predict the shear wave velocity of rock samples. Also, a new correlation between dynamic and static rock properties was established and the shear wave velocity was estimated based on index tests using Gaussian process regression, multivariate linear regression, feedforward back-propagation artificial neural network, and K-nearest neighbor methods. In total, 120 data related to limestone and sandstone samples of the main projects were used for modeling. Water absorption, compressional wave velocity, and density were used as inputs. The outcomes revealed that the PW/SW ratio is equal to 1.69. Various statistics were used to check the method results. The statistical results showed that it is possible to forecast Ed, and SW with high accuracy. Also, the precision of the GPR was higher than the FBP-ANN, statistical analysis, and KNN. Estimation of SW by GPR showed R of 0.992, and RMSE of 0.06, respectively. These four methods were able to estimate the SW with a mean variation percentage of +0.19%. It's important to consider that these models are best suited for predicting SW when the predictor indicators fall within the same range as this study.

Investigation of the compression-shear fracture characteristics of flaw-filled sandstone under confining pressure

In the field of rock engineering, the coupled effects of confining pressure and crack filling material on rock mechanical properties have not received widespread attention. In this study, the central cracked Brazilian disk specimen was used to investigate the fracture properties of sandstone in compression-shear loading mode in a confining pressure and crack-filling environments. Variations in peak load and fracture toughness were quantitatively evaluated for different confining pressures and fillings, accompanied by the establishment of quantitative relationships related to fracture toughness. The trend of linear increase in effective fracture toughness versus confining pressure arises primarily from the contribution of confining pressure. The filling material weakens the fracture toughness of the rock to a certain extent, which is attributed to the accumulation of energy during the failure of the sandstone and the reduction of the stress concentration at the crack tip. The crack propagation trajectory mainly exhibits the characteristics of a stretched wing crack, but shows a transition trend towards a shear crack to a certain extent.

A modeling study of elastoplastic rock failure regime based on finite discrete elements

Numerical simulation has become a pivotal method in geotechnical engineering for analyzing rock failure. However, traditional model lacks comprehensive elastoplastic rock characterization. Therefore, a new finite-discrete element model using the Drucker-Prager criterion and cohesive elements based on a traction-separation law to accurately simulate elastoplastic rock breaking has been built. A scratch test simulation is conducted to describe the mechanical response and failure characteristics of the rock during the cutting process. Rock failure characteristics derived from numerical simulations align closely with experimental observations. These simulations successfully replicate both the ductile and brittle regimes observed during scratch tests. Furthermore, they provide insights into the primary driving factors behind rock failure under different regimes. Variation of cutting force and specific energy obtained from the model are entirely consistent with experimental results, with the specific energy decreasing following a −0.5 power law, thus explaining the energy dissipation rule of the specific energy decrease. In brittle regime, crush zones and the resulting rough shearing surface are observed, plastic damage is also found, revealing that the cutter tip is the key structure influencing them. The model established in this paper holds significant importance for rock mechanical property evaluation and the optimization design of cutter structures.

Feasibility Analysis of Compressed Air Energy Storage in Salt Caverns in the Yunying Area

With the widespread recognition of underground salt cavern compressed air storage at home and abroad, how to choose and evaluate salt cavern resources has become a key issue in the construction of gas storage. This paper discussed the condition of building power plants, the collection of regional data and salt plant data, and the analysis of stability and tightness. Comprehensive analysis and evaluation methods were put forward from four aspects, including ground comprehensive conditions, regional geological conditions and formation lithology, salt mine characteristics, stability, and tightness of salt caverns. The limit equilibrium theory was applied to establish the limit equilibrium failure mode of salt caverns under operating pressure, and the stability coefficient calculation method of the target salt cavern was determined by combining the mechanical characteristics. Based on the physical and mechanical properties of salt rocks, it was found that salt rocks with enough thickness around the salt cavity could be used as sealing rings to ensure the tightness of the salt cavern. Combined with the field water sealing test, the tightness of the target salt cavern is verified. This method has been applied to the salt cavern screening and evaluation of a 300 MW compressed air energy storage power plant project in Yingcheng, Hubei Province, and remarkable results have been obtained, indicating the rationality of the method.

Investigation of Fracture Evolution and Failure Characteristics of Rocks under High-Temperature Liquid Nitrogen Interaction

The utilization of liquid nitrogen as a sustainable and water-free fracturing medium exhibits immense promise in engineering applications. In this investigation, Brazilian split tests and acoustic emission tests were conducted to explore the impact of liquid nitrogen cooling on the internal structure and mechanical properties of rock specimens. To examine the influence of liquid nitrogen cooling on the tensile strength of rocks, displacement-load curves were obtained from samples subjected to varying cycles of high-temperature liquid nitrogen cooling using Brazilian split tests. Acoustic emission experiments were conducted to investigate the characteristics of granite samples exposed to various cycles of high-temperature liquid nitrogen cooling. Based on these findings, the impact of liquid nitrogen cooling on the internal structure of rock masses was analyzed. The findings of this study demonstrate that high-temperature liquid nitrogen thermal treatment significantly modifies the microscopic structure and mechanical properties of rocks, with potential implications for overall stability and reliability. Notably, an observable decline in tensile strength was observed as the number of cycles of high-temperature liquid nitrogen treatment increased. These findings underscore the substantial impact of liquid nitrogen cooling on the behavior of rocks. High-temperature liquid nitrogen treatment effectively promotes the generation of microcracks within rocks, thereby increasing their permeability. During the experiment, granite specimens primarily exhibited shear-type fractures when subjected to high-temperature freeze-thaw cycles induced by liquid nitrogen.

Model of shear strength of ultra-deep fractured sandstone considering fracture morphology

Ultra-Deep fractured tight sandstone gas reservoir, the hotpot unconventional petroleum resources, has been the subject of investigation recently. Due to the presence of natural fractures, drilling in ultra-deep fractured tight sandstone may cause complicated wellbore instabilities, which resulted in a substantial challenge to the efficient exploitation and production of tight gas resources. In order to understand the rock mechanical properties of the ultra-deep fractured tight sandstone under shear stress, tight sandstone core samples were artificially shear-fractured, the 3D characteristics of the fractures were reconstructed by CT scanning, and the fracture aperture and roughness characteristics were systematically evaluated. Simultaneously, triaxial shear strength tests were conducted to evaluate the mechanical properties of fractured rock subjected to varying confining pressures. Considering the contact relationship of the fracture surface during shearing, two indexes of average effective shear angle and contact area ratio were proposed to correct the fracture roughness, and a new fracture shear strength model was constructed, which served as the failure criterion for the fractured formation. The results show that the intact rock has a high compressive strength and can only collapse under conditions of substantial differential stress. Consequently, shear failure of fracture ought to be the primary source of wellbore instability. The fracture aperture's spatial distribution is irregular and conforms to a Gaussian distribution. A clear anisotropy in the fracture roughness is present, and there is less roughness along the initial shear failure path and more roughness perpendicularly. With an increase in normal stress, the fracture's peak shear strength rises, while the friction angle steadily falls. The innovative model's accuracy in predicting shear strength is higher than that of the weak-surface model and the JRC-JCS model.

Analysis of seepage and temperature fields and their influence on collapse pressure while drilling in shale oil reservoirs

Shale wellbore instability poses a considerable challenge in terms of drilling efficiency. However, previous studies typically assumed single-phase flow in the seepage field, neglecting the influences of water saturation and temperature on rock mechanical properties. Furthermore, there has been a lack of quantitative analysis on the degree of correlation between various factors and collapse pressure. To address these gaps, this research introduces an oil–water two-phase seepage model to describe the water saturation distribution surrounding the wellbore,and also constructed a wellbore temperature field model that accounts for mechanical and hydraulic heat sources and validated its accuracy. Furthermore, this paper examines the effects of temperature and seepage on shale mechanical parameters. By incorporating a rock failure criterion that accounts for water saturation, a shale collapse model in thermal-seepage coupling is established, which has also been validated. Finally, this paper utilises the grey relational analysis to quantitatively assess the correlation between different variables and collapse pressure. The results indicate that the temperature field and collapse models are consistent with field data. As drilling time increases, water saturation and invasion distance grow, while the seepage radius hardly increases anymore. When porosity exceeds 0.05, its influence on water saturation weakens. Under unbalanced conditions, capillary force drives drilling fluid into the formation. The maximum wellbore temperature is found near the drill bit and exceeds the formation temperature. The bottom-hole temperature difference increases with the increase in well depth and drilling time, sometimes exceeding 20 °C. Among drilling engineering parameters, pump rate has greater influence on bottomhole temperature than WOB, RPM and ROP. Regarding drilling fluid property parameters, density, flow index and consistency coefficient considerably influence bottom-hole temperature. The influence of seepage on collapse pressure is greater than that of the temperature field. Grey relational analysis reveals that pump rate and water saturation have the highest correlation with collapse pressure. As drilling time increases, the correlation degree of pressure difference, horizontal section length and weight on bit gradually increase, while that of the mud behaviour index, consistency coefficient, rock thermal conductivity and mud yield stress slightly decreases.

An insight into the effect of primary hidden microfissures on mechanical behaviors and failure characteristics of brittle basalt

The anisotropy of mechanical properties of hard rocks is usually affected by the internal primary hidden microfissure structures, which tend to induce secondary rock fracture propagation even surrounding rock instabilities in deep underground excavation. This research attempts to reveal the control effect of primary hidden microfissures on the mechanical behaviors and failure characteristics of cryptocrystalline basalt with hidden microfissures (CBWHM) through a series of laboratory tests. Firstly, a geometric feature extraction method for primary hidden microfissures of basalt rock has been proposed, by which 20 groups of three-dimensional (3D) visualization models of cylindrical basalt samples were established, and the distribution location and volume crack rate (KV) of hidden microfissures were simultaneously determined. On this basis, a series of the uniaxial compression tests of basalt samples with synchronous acoustic emission (AE) were carried out. Results show that the location and development degree of primary hidden microfissures significantly control the mechanical characteristics of basalt, with the relevant micro-fracture events occurring earlier and throughout the whole loading process compared with intact samples. In general, there exists a uniaxial polynomial fitting relationship between the strength of basalt samples and KV. The larger the KV, the lower the peak strength (σc) of basalt, which could decrease by 30–50% compared to that of the intact rock when KV is more than 5‰. However, when the hidden microfissures with large dip angles are located close to the central axis, the influence of KV on the rock strength becomes less, but the crack propagation and failure morphology are still obviously controlled by the hidden microfissures. The statistical analysis results of surface cracks and rock fragment characteristics of the damaged basalt samples further show that the basalt rock is easy to expand, penetrate and eventually fail along the hidden microfissures under loading, and it is tensile fracture that controls the failure of CBWHM. With an increase of KV, the mass percentage (W) of the coarse particle of the sample fragments increases gradually, and the fractal dimension (D) of the rock fragments and KV of the hidden microfisures show an obvious logarithmic curve relationship. The research could provide effective guidance for the interpretation of the meso-mechanism of hard rock failures.

Experimental Simulation on the Stress Disturbance Mechanism Caused by Hydraulic Fracturing on the Mechanical Properties of Shale Formation

Hydraulic fracturing is an indispensable technology for the development of shale oil and shale gas. Knowing the changes in the rock mechanical properties and failure modes during hydraulic fracturing is the key to improving the efficiency of hydraulic fracturing. Based on experiments and simulations, it can be concluded that the injection of fracturing fluid in the hydraulic fracturing caused deformation of the fracture surface, and the rock mechanical properties experienced degradation with a maximum reduction in the rock mechanical properties of 44.24%. As indicated in the experiments, the displacement of the measurement point was decreased with the distance increase between the injection point and the measurement point. According to the numerical simulations, tensile failure is the main failure mode in hydraulic fracturing, but the percentage of shear failure had an obvious increase with the increase in distance between the injection point and the measurement point. Comparing DDS #1 and DDS #5, the DDS #5 measurement point was farther away from the injection point, and the average percentage of shear failure increased from 21.94 to 52.72%. Meanwhile, the increase in the branch fractures also caused shear failure to occur. Comparing Sample 1 and Sample 3, in Sample 3, which had more branch fractures, the average percentage of shear failure increased from 33.12 to 37.58%. Due to the porous medium of the reservoir rock, the enormous pressure generated during the injection of fracturing fluid caused significant deformation of the fracture surface, leading to the tensile failure of the rock. The displacement of the fracture surface caused by the fracturing fluid injection also led to the deformation of the pore throat structure; thus, the shear failure increased when the measurement point was away from the injection point.

Weathering Effects on Engineering Geological Properties of Trachydacitic Volcanic Rocks from the Monte Amiata (Southern Tuscany, Italy)

Variability in lithology and weathering degree affects physical and mechanical properties of rocks. In this study, we investigated the relationships between weathering degree and engineering geological properties of trachydacitic volcanic rocks from Monte Amiata (central Italy) by coupling field and laboratory analyses. We collected in situ Schmidt hammer tests in the field. We evaluated weathering quantifying the percentage of secondary minerals through thermal analysis in the laboratory. We also determined dry density (ρd), specific gravity of solids (Gs), porosity (n) and two-dimensional (2D) porosity as resulted from scanning electron microscopy investigations. The results of our study indicate a negative linear correlation between Schmidt hammer rebound values and secondary mineral percentage. This correlation provides a tool to quantitatively estimate the deterioration of rock uniaxial compressive strength (UCS) as weathering increases. Moreover, thermal analysis turned out to be a quantitative and reproducible method to evaluate weathering degree of magmatic rocks.

Numerical Simulation Analysis of Fracture Propagation in Rock Based on Smooth Particle Hydrodynamics.

The mechanical properties of fractured rock have always been a focal point in the rock mechanics field. Based on previous research, this paper proposes improvements to the SPH method and applies it to the study of crack propagation in fractured rocks. By conducting uniaxial compression tests and simulating crack propagation on various specimens with different crack shapes, the characteristics of crack propagation were obtained. The comparison between the simulated results in this study and existing experimental and numerical simulation results confirms the validity of the SPH method employed in this paper. The present study utilizes the proposed methodology to analyze the influence of the crack angle, width, and orientation on crack propagation. The SPH method employed in this study effectively demonstrates the expansion process of fractured rock under uniaxial compression, providing valuable insights for the engineering applications of SPH.

Prevention and Control of Coal and Gas Outburst by Directional Hydraulic Fracturing through Seams and Its Application.

The regionalized pressure relief, permeability enhancement, and outburst prevention of "three high and one low" coal seams with high gas, high ground stress, high outburst risk, and low permeability have become important problems to be solved urgently in realizing the sustainable development of coal mines. In this study, a combination of theoretical research, RFPA2D-Flow numerical simulation, and field test was used to study the initiation mechanism and propagation law of directional hydraulic fracturing fractures through the seam. The results show that fracture initiation depends on the axial and radial horizontal stress of the fracture hole, the physical and mechanical properties of coal and rock, and the inside of the weakest layer. Single-hole hydraulic fracturing can achieve a pressure relief radius of 7-8 m, but there is a stress concentration zone outside, which is not conducive to regional pressure relief and permeability enhancement. Directional hydraulic fracturing with multiple holes produces an approximate cylindrical compression and crushing ring and a penetrating fracture surface along the center line of the pressure crack hole and the directional hole, which better eliminates the phenomenon of stress concentration in nondirectional hydraulic fracturing. The technology was applied to the 2238 auxiliary roadway of Chengzhuang Mine of Jinmei Group, and the field implementation results showed that field implementation results showed that directional hydraulic fracturing through the seam reduced the gas content in the coal seam to a great extent, and the coal seam gas content was reduced by about 42.3%, indicating that the technology can effectively reduce the risk of coal and gas outbursts.

Emphasizing the Creep Damage Constitutive Model of Hydro-Mechanical Properties of Rocks: A Case Study of Granite Gneiss

The constitutive model of rock materials can describe the mechanical behavior of rocks in creep tests. Also, it is one of the important means to study the deformation and strength characteristics of rocks in complex stress environments. This paper is based on the analysis of the porosity variation characteristics of the internal structure under the coupling effect of rock hydro-mechanical properties. The concept of the hydro-mechanical properties variable is proposed, and the relationship between the coupling variable, damage and plastic deformation is established. By introducing the coupling variable, instantaneous damage variable and time-dependent damage variable into the yield surface equation, as well as the plastic potential energy equation and the stiffness matrix of the elastic–plastic creep constitutive equation, a hydro-mechanical properties creep damage coupling model was established to simulate the creep mechanical properties of rock under coupling. Based on the triaxial creep test results of granite gneiss, the model parameters are determined. By comparing the test results with numerical results, it was revealed that the model can better describe the creep mechanical properties of rocks under the coupling effect of hydromechanical properties.

Quantification of rock heterogeneity and application in predicting rock mechanical properties

Quantification of rock heterogeneity and application in predicting rock mechanical properties

Rock-physics modeling of elastic properties of multiscale fractured rocks

Fractures, often existing across various scales, control the mechanical and fluid flow properties of upper crustal rocks. Inferring fracture properties at different scales from multiband geophysical measurements is essential for many fields of earth and energy sciences. Under the framework of multiscale homogenization, we develop a rock-physics model to characterize the elastic and anisotropic properties of multiscale fractured rocks following a certain statistical law. The isotropic differential effective medium theory is used to model the elasticity of randomly oriented microcracks, and the linear slip theory is used to model the elasticity of oriented or randomly oriented macroscale fractures. For a multiscale fractured rock covered by six orders of fracture length (10−4 m to 102 m), it is found that velocity exhibits a decreasing trend with the increment of scale. Nevertheless, the different statistical distribution of multiscale fractures significantly affects the velocity and anisotropy variation pattern with scale. The velocity for the fractal distribution decreases significantly at the seismic scale, whereas, for the log-Gaussian distribution, the dramatic change in velocity occurs more at the ultrasonic and logging scales, depending on the length of the dominant fractures. We apply our methodology to interpret a set of multiscale geophysical data of P-wave velocity in a fractured carbonate formation and estimate that the multiscale fracture is possibly distributed in a log-Gaussian manner. Our elastic and anisotropic modeling strategy has the potential to predict the distribution pattern of fractures, especially by reconciling the multifrequency geophysical measurements (e.g., ultrasonic, logging, and seismic data).

Upscaling rocks mechanical properties to study Underground Hydrogen Storage feasibility

Underground Hydrogen Storage (UHS) is a feasible option for large-scale energy storage considering the advancements of the large-scale production of green hydrogen. One of the main engineering objective is to ensure the continuous safety of the storage, such that subsurface operations can be carried under safe stress regimes. Despite obvious similarities to Carbon Capture and Storage, a few differences make the task a new research challenge.
 The first one relates to the cyclic nature of UHS, which is expected to be carried at variable frequency and injection/production loads. Current models are not adequate for the lower frequency range considered in this application. In that case, the visco-plastic nature of rocks becomes non-negligible and needs to be taken into account.
 As a second observation, UHS revolves around a new gas, much lighter than CH4 and super critical CO2. Hydrogen’s atoms are so small they can diffuse even inside rock and this absorption causes rock matrix mechanical properties to weaken. This process is know as Hydrogen Embrittlement. When unaccounted for, such physical phenomenon could lead to catastrophic failure of the caprock, which is supposed to maintain stability to ensure safe storage. The caprock being responsible for the confinement of the hydrogen in the reservoir, development of cracks would enhance greatly permeability of an otherwise impermeable medium, resulting in an environmental disaster as the hydrogen suddenly leaks towards the subsurface and through groundwater aquifers.
 No empirical model is able to capture those two behaviours at the macro-scale since they are both phenomena principally related to grain-scale physics. As such, this contribution presents a Digital Rock Physics framework to upscale rock mechanical properties from the grain-scale. Rocks of interest are microCT-scanned to extract the digitized microstructure. Direct numerical simulations of elasto-plasticity are performed for different stress paths in order to compute the full yield surface (see Figure 1) instead of just the Uniaxial Compressive Strength. While most studies use Discrete Element Modelling to consider grain contacts explicitly, our simulator uses Finite Element Modelling which allows more flexibility in the approach to model multiphysics processes present during UHS. The contacts are modelled instead as an upscaled plastic law. Details of the numerical algorithms are presented in [1].
 As a first case study for this framework, we present a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of granular materials. The framework is then coupled with a numerical erosion algorithm that simulates homogeneous precipitation of mineral matter to represent cementation. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. This study contributes to preliminary results on Hydrogen Embrittlement which directly influences those same mechanical properties. However more work is needed to model realistically the Hydrogen Embrittlement, which is the aim of our new PhD project OCEAN. The process will be observed experimentally at the micro-scale in order to calibrate the simulator. MicroCT-scan images will determine the spatial distribution of the phenomenon. Visco-plasticity will be implemented from [2] to go one step further and determine the effect of Hydrogen Embrittlement during cyclic injection/production of hydrogen.

Study on the deterioration characteristics and mechanism of physical and mechanical properties of natural gypsum rock under freeze–thaw action

The physical and mechanical properties, microstructure degradation characteristics, and damage mechanism of natural gypsum rock under 15, 30, and 45 freeze–thaw cycles were studied by using a self-developed program-controlled freeze–thaw experimental device. The results show that compared with natural gypsum, with the increase in freeze–thaw cycles, the mass of the sample increases first and then decreases, Poisson's ratio increases, the density decreases, and the longitudinal wave velocity, the uniaxial compressive strength, and the elastic modulus decrease. Indexes decrease greatly in the early stage of freeze–thaw and slowdown in the later stage. With the increase in freeze–thaw cycles, the slope of the pre-peak stage of the stress–strain curve decreases, the compaction stage is significantly prolonged, and the brittleness of the post-peak failure process decreases. Failure cracks of the specimens after freeze–thaw cycles are mainly inclined shear failure and axial tensile cracks. The analysis shows that the freeze–thaw damage of gypsum rock is caused by the superposition and mutual promotion of water-rock softening and pore frost heave. The specific performance is the damage of crystal particles and the loose pore structure at the mesoscale. On this basis, the cracks gradually develop and expand, which is consistent with the macroscopic test results of the sample. The research results have reference significance for the construction scheme design and frost damage prevention of gypsum surrounding rock tunnels in cold regions.