{"title":"Multi-stage evolution of pore structure of microwave-treated sandstone: Insights from nuclear magnetic resonance","authors":"Yao Zhang , Yanan Gao , Liyuan Yu","doi":"10.1016/j.ijrmms.2024.105952","DOIUrl":"10.1016/j.ijrmms.2024.105952","url":null,"abstract":"<div><div>Microwave fracturing has great potential in improving the efficiency of hard rock breaking. However, the pore evolution, which can be regarded as the damage accumulation and progressive failure of the rock subjected to microwave irradiation, remains unclear. In this study, nuclear magnetic resonance (NMR) is employed to investigate the pore evolution and fracture mechanism of the sandstone under different microwave power levels. The results show that the pore evolution of the specimens, including distribution of pore size, the weight in volume of various-sized pore, and porosity, exhibits different changing trends under various microwave power levels. The pore evolution of the specimens under microwave irradiation can be categorized into four phases: overall pore expansion, localized pore closure in the internal region, micro-cracks propagation induced by thermal stress, and macro-cracking (or melting). Moreover, pore evolution also plays a crucial role in the decomposition and evaporation of bound water, particularly when the specimens experience fractures triggered by thermal stress induced by the microwave treatment (TSIMT). The employing of NMR imaging (NMRI) description also provides an auxiliary and effective illustration on the pore evolution of the specimens under microwave irradiation. Finally, the mechanism of microwave-assisted rock breaking under different power levels is comprehensively discussed based on the NMR results from a microscopic perspective. It is anticipated that the findings of this study can provide valuable insights for enhancing the efficiency of microwave-assisted rock breaking.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105952"},"PeriodicalIF":7.0,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142573463","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Robert McMillan , Erik Eberhardt , Ryan Campbell , Avesiena Primadiansyah
{"title":"Monitoring stress-induced brittle rock mass damage for preventative support maintenance","authors":"Robert McMillan , Erik Eberhardt , Ryan Campbell , Avesiena Primadiansyah","doi":"10.1016/j.ijrmms.2024.105927","DOIUrl":"10.1016/j.ijrmms.2024.105927","url":null,"abstract":"<div><div>Stress-induced brittle fracturing near an excavation boundary results in a volume increase, known as bulking. Excessive bulking places added demand on the rock support, which, if not detected and addressed through preventative support maintenance (i.e., proactively added reinforcement), can cause the support to fail, leading to a safety hazard and costly production delays for underground mining operations. For caving mines, these project risks are exacerbated during cave establishment due to the large abutment stress from undercutting that redistributes and concentrates stresses near excavations critical for production. This paper reports the findings from research conducted to develop and improve geotechnical monitoring practices to support preventative support maintenance in deep mining operations. This research uses a unique geotechnical monitoring database collected for the Deep Mill Level Zone panel cave mine. The data was collected across a large footprint during the mine's ramp-up period and represents an initial step toward best practices for data collection at cave mines operating in high-stress environments. Borehole camera surveys supplemented by multi-point borehole extensometers have been used to determine the depth of stress fracturing in pillar walls as a function of the distance away from the undercut. Convergence measurements and LiDAR scanning are used to characterize the corresponding rock mass bulking. The results show that the interpretation of monitoring data can be used to identify the long-term depth of stress fracturing and bulking trends in response to undercut advances. These show that direct measures of stress-induced fracturing damage provide an early indication of excavations vulnerable to bulking and that LiDAR scanning is an effective method for capturing the onset of bulking and anticipating local areas likely to experience greater deformation demand as bulking progresses. Proactive and strategic geotechnical monitoring based on the long-term depth of stress-induced fracturing trends is proposed to assist with preventative support maintenance practices.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105927"},"PeriodicalIF":7.0,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142552488","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Kai Qiu , Shuchen Li , Zhongzhong Liu , Meng Yuan , Shisen Zhao , Zeen Wan
{"title":"An elastoplastic solution for lined hydrogen storage caverns during excavation and operation phases considering strain softening and dilatancy","authors":"Kai Qiu , Shuchen Li , Zhongzhong Liu , Meng Yuan , Shisen Zhao , Zeen Wan","doi":"10.1016/j.ijrmms.2024.105949","DOIUrl":"10.1016/j.ijrmms.2024.105949","url":null,"abstract":"<div><div>Underground hydrogen energy storage (UHES) in lined rock caverns (LRCs) represents a crucial solution to the challenge of unstable and uneven clean energy generation. Nevertheless, the attainment of enhanced storage efficiencies frequently necessitates the utilization of elevated hydrogen storage pressures. Consequently, a comprehensive understanding of the elastic-plastic mechanical response of surrounding rock under hydrogen pressure is of paramount importance for ensuring the safety of UHES. In this study, an elastoplastic solution of LRCs during construction and operation phases is established. Two essential phenomena affecting the post-peak mechanical responses of surrounding rock, strain softening and dilatancy, are coupled into the plastic solution. A computational process is developed and its accuracy is validated through comparison with numerical models. The influence of surrounding rock quality parameters, strain softening and dilatancy parameters, concrete quality parameters and hydrogen pressure on the radius of the plastic softening zone (<em>R</em><sub><em>s</em></sub>) and plastic residual zone (<em>R</em><sub><em>r</em></sub>) were analyzed. Results show that higher surrounding rock quality can effectively reduce both <em>R</em><sub><em>s</em></sub> and <em>R</em><sub><em>r</em></sub>. Nevertheless, when the surrounding rock quality already reaches a high standard, such as <em>c</em><sub>1</sub> > 3.5 MPa, <em>φ</em><sub>1</sub> > 65°, or <em>E</em> > 55 MPa, it becomes inefficient to overly pursue further improvements in the surrounding rock quality. Furthermore, the strain softening and dilatancy phenomena only affect <em>R</em><sub><em>r</em></sub>. Additionally, the concrete lining with higher stiffness can share a larger portion of the hydrogen pressure, thus reducing both <em>R</em><sub><em>s</em></sub> and <em>R</em><sub><em>r</em></sub>. Notably, When the elastic modulus of concrete increases from 20 MPa to 40 MPa, <em>R</em><sub><em>r</em></sub> decreases by 31.98 % and <em>R</em><sub><em>s</em></sub> decreases by 20.96 %. Moreover, the critical hydrogen pressure (<em>P</em><sub><em>Hcr</em></sub>) at which the surrounding rock begins to enter a plastic state is proportional to the ground stress (<em>P</em><sub>0</sub>). Specifically, when <em>P</em><sub>0</sub> is increased sequentially from 2.5 MPa to 3.0 MPa and 3.5 MPa, <em>P</em><sub><em>Hcr</em></sub> sequentially becomes 2.4 MPa, 4.0 MPa, and 5.0 MPa. The findings presented in this study contribute to improving the safety of LRCs during construction and operation.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105949"},"PeriodicalIF":7.0,"publicationDate":"2024-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142538808","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhen Yang , HanYi Wang , Mukul Sharma , Erdogan Madenci
{"title":"A fault activation-shearing-sliding peridynamic model exploring the role of static and kinetic frictional contacts","authors":"Zhen Yang , HanYi Wang , Mukul Sharma , Erdogan Madenci","doi":"10.1016/j.ijrmms.2024.105946","DOIUrl":"10.1016/j.ijrmms.2024.105946","url":null,"abstract":"<div><div>Understanding fault dynamics is essential for comprehending the underlying mechanisms of seismic events. This study introduces a novel fault activation-shearing-sliding model within a peridynamic (PD) framework, characterized by distinctly defined static and kinetic frictional behaviors. Static friction bonds are developed to sustain normal forces perpendicular to the fault plane and to manage tangential frictional forces along the fault's geometry. The failure of these bonds is directly linked to fault activation, while the ensuing sliding phase is governed by a short-range kinetic friction model. Additionally, an adaptive identification method is proposed to accurately determine local unit normal vectors on arbitrarily shaped contact surfaces. The effectiveness and applicability of the model are validated through fault activation and plate sliding friction tests. The model is further utilized to investigate the effects of local geometry, roughness, and friction coefficients on fault behavior, with comparisons to experimental results. Observations indicate that the dominant factors influencing fault shear resistance vary across stages, primarily involving static friction during activation, compaction deformation during shearing, and kinetic friction during sliding. When shear resistance is primarily governed by friction, it exhibits heightened sensitivity to various shear forces, including those from indirect loading disturbances.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105946"},"PeriodicalIF":7.0,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535664","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Wenzhuo Cao , Sevket Durucan , Ji-Quan Shi , Anna Korre , Thomas Ratouis , Vala Hjörleifsdóttir
{"title":"Geothermal fluid extraction and injection-related fracture slip susceptibility and seismicity in naturally fractured rocks","authors":"Wenzhuo Cao , Sevket Durucan , Ji-Quan Shi , Anna Korre , Thomas Ratouis , Vala Hjörleifsdóttir","doi":"10.1016/j.ijrmms.2024.105939","DOIUrl":"10.1016/j.ijrmms.2024.105939","url":null,"abstract":"<div><div>Understanding fracture slip susceptibility in geothermal reservoirs is central to the control of fluid injection-induced seismicity. To investigate the role of regional fracture systems on induced seismicity, a coupled thermo-hydro-mechanical (THM) model containing fracture networks was developed, which features direct coupling between different physics for explicit fractures, fractured rocks (porous matrix blocks with small-scale fractures) and their interactions, as well as indirect coupling through changes of material properties, such as stress-dependent fracture and rock permeabilities. The model was applied to simulate geothermal fluid extraction and re-injection in a natural fracture system comprised of three dominant fracture sets at the Hellisheiði geothermal field over a 10-year period (2011–2021), utilising field recorded monthly production and re-injection rates. Based on the model results, the slip susceptibility of regional fracture systems was examined under reservoir conditions before and after the start of fluid re-injection across different time scales, i.e., over short (1 month), intermediate (1 year) and long-term (10 years). Two model scenarios, one with cooling contraction and one without, were considered to examine the relative contribution of cooling contraction and fluid overpressure on fracture slip susceptibility. Results have shown that fracture networks act as preferential fluid flow paths that influence fluid pressure and stress distribution and fracture slip susceptibility in geothermal reservoirs. NE-SW and N-S trending fractures at Hellisheiði are susceptible to slippage before the start of fluid re-injection. During fluid re-injection, the distribution of fractures with enhanced slip susceptibility gradually shifts from surrounding the re-injection region to forming a two-lobed pattern in the fault-normal direction around the re-injection region, indicating the dominant role of cooling contraction over fluid overpressure on the fracture slip susceptibility in the intermediate- and long-term.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105939"},"PeriodicalIF":7.0,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535663","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Rui Zhang , Paul A. Bosomworth , Juliane Weber , Jan Ilavsky , Si Athena Chen , Alexis Flores-Betancourt , Elliot Paul Gilbert , Jitendra Mata , Mark L. Rivers , Peter J. Eng , Lawrence M. Anovitz
{"title":"The role of annealing and grain boundary controls on the mechanical properties of limestones and marbles","authors":"Rui Zhang , Paul A. Bosomworth , Juliane Weber , Jan Ilavsky , Si Athena Chen , Alexis Flores-Betancourt , Elliot Paul Gilbert , Jitendra Mata , Mark L. Rivers , Peter J. Eng , Lawrence M. Anovitz","doi":"10.1016/j.ijrmms.2024.105926","DOIUrl":"10.1016/j.ijrmms.2024.105926","url":null,"abstract":"<div><div>Chemical and mechanical processes are coupled in many geological and geochemical environments. Reactive processes anneal defects and restructure grain boundaries, modifying their elastic properties, levels of internal friction, wave propagation rates and fracture behaviors. The nature of these changes is, however, contingent on the initial state of the rock. In this study, impulse excitation was used to measure changes in mechanical properties as a function of dry and steam heating time at 300 °C for three carbonate rocks: Carrara marble, Carthage marble (Burlington Limestone), and Texas Cream limestone (Edwards limestone), with initial porosities of about 1 %, 3 %, and 27 %, respectively. Frequency-dependent phenomena along with mineral recrystallization were observed. This was coupled with small-angle X-ray and neutron scattering analysis to determine the relationship between changes in bulk mechanical and microstructural properties. An observed decrease in both the Young's and shear moduli in the experimental limestones as a function of heating time reflects and quantifies a reduction in the stiffness of the rock due to annealing. The internal friction of the samples first increases then decreases with reaction time, reflecting defect annealing, but this was balanced against grain boundary dissolution apparently driven by condensation of steam in the confined grain-boundary environment. This suggests an increase in the boiling point under confinement leading to dissolution, increased porosity and widening of the grain boundaries. In addition, comparison of small-angle X-ray and neutron scattering results suggests that it is inappropriate to assume that pores are empty for quantitative analysis of typical small-angle scattering samples.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105926"},"PeriodicalIF":7.0,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535665","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yaolan Tang , Jianchun Li , Dapeng Wang , Congying Li , Chunshun Zhang
{"title":"An energy-driven crushing-plasticity coupling model for grain crushing in porous rocks","authors":"Yaolan Tang , Jianchun Li , Dapeng Wang , Congying Li , Chunshun Zhang","doi":"10.1016/j.ijrmms.2024.105931","DOIUrl":"10.1016/j.ijrmms.2024.105931","url":null,"abstract":"<div><div>This research develops an energy-driven constitutive model designed to tackle the complex phenomenon of grain crushing in porous rocks. Initially, a novel coupled relationship is proposed to integrate various energy dissipation mechanisms, including both plastic and crushing effects, using spherical polar coordinates. This approach results in a robust coupling of energy dissipation, providing a comprehensive depiction of the influence of grain crushing on plastic deformation. An energy-based yield criterion is then formulated by comparing elastic potential energy contours with experimental findings, and the behaviour of crushing hardening is examined through energy evolution. Flow rules are subsequently derived, both independently and with consideration of plasticity-crushing coupling. Finally, validation against a range of experimental tests highlights the model's versatility. The proposed model enhances the understanding of rock-crushing issues from an energy perspective and demonstrates simplicity with only 4 or 5 easily calibrated parameters.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105931"},"PeriodicalIF":7.0,"publicationDate":"2024-10-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535670","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Exposure behavior and drilling efficiency of basalt fiber composite impregnated diamond bits in hard granite","authors":"Yinlong Ma , Jie Ren , Qingquan Zhou","doi":"10.1016/j.ijrmms.2024.105950","DOIUrl":"10.1016/j.ijrmms.2024.105950","url":null,"abstract":"<div><div>Impregnated diamond bits (IDBs) are widely used for drilling in hard formations. To improve the drilling efficiency and exposure behavior of IDBs in granite, three types of Cu-based basalt fiber (BF) composite IDBs were designed and prepared by using the medium-frequency induction hot-pressing and sintering method, in which 0 wt% BF and 25 vol% diamond were used in 0BF25D IDB, 1 wt% BF and 25 vol% diamond were used in 1BF25D IDB, 1 wt% BF and 20 vol% diamond were used in 1BF20D IDB. The drilling efficiency of each IDB was tested under different drilling pressures (WOB), and the exposure behavior of IDBs was investigated by scanning electron microscopy and ultra field microstructural characterization. Results show that drilling granite with grade 9 drillability, low drilling pressure can drill successfully with the addition of BF. The average rate of penetration (ROP) of 1BF20D IBD under 6 kN WOB was 4.78 m/h, which was improved by 60 %∼80 %, energy consumption decreased by 66 % for each meter, torque (TOB) decreased, and rotational speed (RPM) was more stable during the drilling process. The addition of BF and reasonable diamond concentration enhanced the holding power of the diamond which promoted the exposure of the diamond. The average exposed height of the diamond in 1BF20D IDB reached 121.2 μm with 22 %–29 % of the whole diamond.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105950"},"PeriodicalIF":7.0,"publicationDate":"2024-10-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535789","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lei Peng , Mingyao Li , Jianping Zuo , Dejun Liu , Jena Jeong
{"title":"Determination of the REV size for heterogeneous rocks with different grain sizes: Deep learning and numerical approaches","authors":"Lei Peng , Mingyao Li , Jianping Zuo , Dejun Liu , Jena Jeong","doi":"10.1016/j.ijrmms.2024.105940","DOIUrl":"10.1016/j.ijrmms.2024.105940","url":null,"abstract":"<div><div>Accurate determination of the representative elementary volume (REV) size plays a pivotal role in analysing the mechanical properties and failure processes of heterogeneous rocks in complex engineering environments. In this study, a novel microstructure modelling strategy (NMMS) for determining the REV size is proposed by combining deep learning and an improved phase-field method (PFM). Micro- and macroscale experiments are systematically conducted to determine the real microstructural characteristics and mechanical properties of heterogeneous rocks with different grain sizes. On the basis of this experimental evidence, geometric models of different sizes were reconstructed through deep learning to avoid the limitations of human-based methods, and an improved PFM was used for numerical calculations. These models were then employed to perform numerical tests under uniaxial loading conditions, and the coefficient of variation was introduced to determine the REV size of heterogeneous rocks with different grain sizes. The research findings indicate that the final REV size is the maximum value of the REVs defined by the evaluation properties within an acceptable coefficient of variation. At a criterion of 5% for the coefficient of variation, the REV sizes are 60 mm<span><math><mo>×</mo></math></span>60 mm, 70 mm<span><math><mo>×</mo></math></span>70 mm, and 90 mm<span><math><mo>×</mo></math></span>90 mm for fine-medium-grained (FMG), medium-grained (MG), and coarse-grained (CG) rocks, respectively. Furthermore, the REV determined by the NMMS was applied to investigate the effects of microstructure on macromechanical properties and damage evolution under triaxial loading conditions. The numerical results show that the NMMS can accurately predict the macromechanical properties and microcracking patterns of heterogeneous rocks, especially the intracrystalline cracks in feldspar, the interfacial cracks in gravel, and the “voids” of cracks in biotite. This research can provide some basic references for the optimal choice of the REV size of heterogeneous rocks.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105940"},"PeriodicalIF":7.0,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535666","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The hydromechanical coupled numerical method in pseudo-3D axisymmetric domain with cracks extension and coalescence applies to the decompression failure problem","authors":"Salavat Y. Ishbulatov","doi":"10.1016/j.ijrmms.2024.105945","DOIUrl":"10.1016/j.ijrmms.2024.105945","url":null,"abstract":"<div><div>The stress-strain state of the saturated porous media determines the behavior of fracturing, which defines the efficiency of developing tight oil, shale, coalbed, and thermal energy fields. Therefore, reliable hydromechanical coupled simulations with destruction reconstruction are critical.</div><div>The proposed innovative simulator has a strong interrelation between fluid flow and rock deformation of porous media and realizes a fully coupled pseudo-transient numerical method by high-performance computing (HPC) tools. To increase the detail of the results in the problem, a finite difference numerical algorithm was implemented in the axisymmetric cylindrical domain, which reduces from three to two dimensions without loss of precision. Highly efficient parallelization using CUDA on the GPU computes meshes of up to one billion cells, allowing the simulation of a total core sample to sub-micrometer resolution in an appropriate time. The algorithm has been validated to find the exact solution to the cylinder problem. The proposed model accounts for cracks propagation with their coalescence within a single computational static grid, which keeps timing close to the continuous model.</div><div>This comprehensive implementation enables solving industrial problems, such as modeling core sample damage during rapid decompression. High-resolution simulations help reconstruct fracture propagation, analyze the initial stress state, and identify critical damage factors. The comparison with the exact solution to the cylinder problem confirmed the reliability of the algorithm. The calculation results show a strong dependence of decompression failure on the coalescence and elongation of cracks, influenced significantly by the rock's cohesion. Microcracks length and distribution play a decisive role in the decompressive destruction behavior of the rock sample. For the first time, the simulations demonstrated the decompressive destruction of a core sample during an uncontrolled, rapid core retrieval operation.</div></div>","PeriodicalId":54941,"journal":{"name":"International Journal of Rock Mechanics and Mining Sciences","volume":"183 ","pages":"Article 105945"},"PeriodicalIF":7.0,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535790","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}