Rock Mechanics
Mohammad Reza Zeerak; Mohammad Fatehi Marji; Manouchehr Sanei; Mehdi Najafi; Abolfazl Abdollahipour
Abstract
The Extended Finite Element Method (XFEM) is a leading computational approach for studying crack growth in rocks, as it can effectively model complex crack paths and discontinuities without the need for re-meshing. In this context, XFEM is particularly well-suited for simulating the development of hydraulic ...
Read More
The Extended Finite Element Method (XFEM) is a leading computational approach for studying crack growth in rocks, as it can effectively model complex crack paths and discontinuities without the need for re-meshing. In this context, XFEM is particularly well-suited for simulating the development of hydraulic fractures. XFEM is employed to investigate crack initiation, propagation, and aperture size in rock formations, with validation using a Boundary Element Method (BEM)-based approach. Three scenarios are analyzed for crack orientation and interaction in: single cracks at and crack displacement behavior at and multiple cracks at and . Displacement in the vertical direction (U2) and stress distribution around the crack tip in the S22 direction are examined to understand fracture mechanics parameters. The findings highlight that crack at higher angles, such as , exhibit more straightforward propagation, while those at or beyond often require additional stress to continue growing. The comparison between XFEM and BEM results confirms the reliability of the numerical approach, demonstrating strong agreement in predicting fracture behavior in rock materials. The results provide deeper insights into fracture evolution, stress intensity factors, and fracture toughness in geological media. These simulations advance computational fracture mechanics, contributing to optimizing hydraulic fracturing techniques for improved efficiency and safety in subsurface formations. This study is limited to 2D geometries and isotropic materials, potentially missing 3D heterogeneous subsurface complexities. Future work could explore 3D models, anisotropy, and fluid pressure/thermal effects to improve crack growth predictions.
Rock Mechanics
Masoud Yazdani; Mohammad Fatehi Marji; Mehdi Najafi; Manouchehr Sanei
Abstract
Around 70% of the world's hydrocarbon fields are situated in reservoirs containing low-strength rocks, such as sandstone. During the production of hydrocarbons from sandstone reservoirs, sand-sized particles may become dislodged from the formation and enter the hydrocarbon fluid flow. Sand production ...
Read More
Around 70% of the world's hydrocarbon fields are situated in reservoirs containing low-strength rocks, such as sandstone. During the production of hydrocarbons from sandstone reservoirs, sand-sized particles may become dislodged from the formation and enter the hydrocarbon fluid flow. Sand production is a significant issue in the oil industry due to its potential to cause erosion of pipes and valves. Separating grains from oil is a costly process. Oil and gas companies are motivated to reduce sand production during petroleum extraction. Hydraulic fracturing is one of the parameters that can influence sand production. However, understanding the complex interactions between hydraulic fracturing mechanisms and sand production around wellbores is critical for optimizing reservoir recovery and ensuring the integrity of production wells. This article explores the integrated simulation approach to model hydraulic fracturing processes and assess their effects on sand production. Two-dimensional models were created using the discrete element method in PFC2D software for this research. The fractures' length in the models varies based on the well's radius. The angle between two fractures at 90 and 180 degrees to each other was also modeled. In the first case, the length of the fracture is less than the radius of the well, in the second case, the values are equal and finally, the fracture length is assumed to exceed the well radius. The calibrated and validated results demonstrate the change in sand production rate in comparison to the unbroken state.
E. Emami Meybodi; Syed Kh. Hussain; M. Fatehi Marji; V. Rasouli
Abstract
In this work, the machine learning prediction models are used in order to evaluate the influence of rock macro-parameters (uniaxial compressive strength, tensile strength, and deformation modulus) on the rock fracture toughness related to the micro-parameters of rock. Four different types of machine ...
Read More
In this work, the machine learning prediction models are used in order to evaluate the influence of rock macro-parameters (uniaxial compressive strength, tensile strength, and deformation modulus) on the rock fracture toughness related to the micro-parameters of rock. Four different types of machine learning methods, i.e. Multivariate Linear Regression (MLR), Multivariate Non-Linear Regression (MNLR), copula method, and Support Vector Regression (SVR) are used in this work. The fracture toughness of mode I and mode II (KIC and KIIC) is selected as the dependent variable, whereas the tensile strength, compressive strength, and elastic modulus are considered as the independent variables, respectively. The data is collected from the literature. The results obtained show that the SVR model predicts the values of KIC and KIIC with the determination coefficients (R2) of 0.73 and 0.77. The corresponding determination coefficient values of the MLR model and the MNLR model for KI and KII are R2 = 0.63, R2 = 0.72, and R2 = 0.62,0.75, respectively. The copula model predicts that the value of R2 for KI is 0.52, and for KII R2=0.69. K-fold cross-validation testing method performs for all these machine learning models. The cross-validation technique shows that SVR is the best-designed model for predicting the fracture toughness mode-I and mode-II.
M.R. Shahbazi; M. Najafi; M. Fatehi Marji; A. Abdollahipour
Abstract
The in-situ coal is converted to the synthetic gas in the process of underground coal gasification (UCG). In order to increase the rate of in-situ coal combustion in the UCG process, the contact surfaces between the steam, heat, and coal fractures should be raised. Therefore, the number of secondary ...
Read More
The in-situ coal is converted to the synthetic gas in the process of underground coal gasification (UCG). In order to increase the rate of in-situ coal combustion in the UCG process, the contact surfaces between the steam, heat, and coal fractures should be raised. Therefore, the number of secondary cracks should be increased by raising the heat and existing steam pressure during the process. This paper emphasises on the secondary crack growth mechanism of the pre-existing cracks in the coal samples under different loading conditions. Different geometric specifications such as the length of the pre-existing cracks (coal cleats) and their inclinations are considered. The numerical modeling results elucidate that the first crack growths are the wing cracks (also called the primary or tensile cracks) formed due to unbonding the tensile bonds between the particles in the assembly. Ultimately, these cracks may lead to the cleat coalescences. On the other hand, the secondary or shear cracks in the form of co-planar and oblique cracks may also be produced during the process of crack growth in the assembly. These cracks are formed due to the shear forces induced between the particles as the initial cleat length is increased and exceed the dimension of coal blocks. The cavity growth rate increases as the secondary cracks grow faster in the coal blocks. In order to achieve the optimum conditions, it is also observed that the best inclination angle of the initial coal cleat changes between 30 to 45 degrees with respect to the horizon for the coal samples with the elasto-brittle behavior.
H. Yousefian; M. Fatehi Marji; H. Soltanian; A. Abdollahipour; Y. Pourmazaheri
Abstract
Determination of the borehole and fracture initiation positions is the main aim of a borehole stability analysis. A wellbore trajectory optimization with the help of the mud pressure may be unreasonable since the mud pressure can only reflect the degree of difficulty for the initial damage to occur at ...
Read More
Determination of the borehole and fracture initiation positions is the main aim of a borehole stability analysis. A wellbore trajectory optimization with the help of the mud pressure may be unreasonable since the mud pressure can only reflect the degree of difficulty for the initial damage to occur at the wellbore rather than the extent of the wellbore damage. In this work, we investigate the failure extension in different arbitrary inclination boreholes under different in-situ stress regimes. Assuming the plane strain condition, the Mohr-Coulomb, Mogi-Coulomb, and Modified Lade rock failure criteria are utilized. We present an analytical equation to determine the optimum drilling trajectory of an Iranian oilfield. In order to predict the degree of wellbore damage, the initial shear failure location, failure width, and failure depth of arbitrary wellbores are determined. Then a new model is derived to calculate the initial failure area of a directional wellbore because it is more efficient in a wellbore stability analysis. The results obtained show that in the target oilfield, the vertical and low-deviated direction is the optimum drilling path. According to the results of this work, optimization of the wellbore trajectory based on the estimated failure zone is a reasonable method if a considerable failure zone takes place around the borehole wall.
Rock Mechanics
M. Lak; M. Fatehi Marji; A.R. Yarahamdi Bafghi; A. Abdollahipour
Abstract
The explosion process of explosives in a borehole applies a very high pressure on its surrounding rock media. This process can initiate and propagate rock fractures, and finally, may result in the rock fragmentation. Rock fragmentation is mainly caused by the propagation of inherent pre-existing fractures ...
Read More
The explosion process of explosives in a borehole applies a very high pressure on its surrounding rock media. This process can initiate and propagate rock fractures, and finally, may result in the rock fragmentation. Rock fragmentation is mainly caused by the propagation of inherent pre-existing fractures of the rock mass and also from the extension of the newly formed cracks within the intact rock due to the explosion. In this work, the process of extension of blast-induced fractures in rock masses is simulated using the discrete element method. It should be noted that, in this work, fracture propagation from both the rock mass inherent fractures and newly induced cracks are considered. The rock mass inherent fractures are generated using the discrete fracture network technique. In order to provide the possibility of fracture extension in the intact rock blocks, they are divided into secondary blocks using the Voronoi tessellation technique. When the modeling is completed, the fracture extension processes in the radial and longitudinal sections of a borehole are specified. Then a blast hole in an assumed rock slope is modeled and the effect of pre-splitting at the back of the blast hole (controlled blasting) on the fracture extension process in the blast area is investigated as an application of the proposed approach. The modeling results obtained show that the deployed procedure is capable of modeling the explosion process and different fracture propagations and fragmentation processes in the rock masses such as controlled blasting.
