Rock Mechanics
Mohsen Khanizadeh Bahabadi; Alireza Yarahamdi Bafghi; Mohammad Fatehi Marji; Hossein Shahami; Abolfazl Abdollahipour
Abstract
Complexity of geomaterial’s behavior is beyond the capabilities of conventional numerical methods alone for realistically model rock structures. Coupling of numerical methods can make the numerical modeling more realistic. Discontinuous Deformation Analysis (DDA) and Displacement Discontinuous ...
Read More
Complexity of geomaterial’s behavior is beyond the capabilities of conventional numerical methods alone for realistically model rock structures. Coupling of numerical methods can make the numerical modeling more realistic. Discontinuous Deformation Analysis (DDA) and Displacement Discontinuous Method (DDM) are hybridized for modeling block displacement and crack propagation mechanism in a blocky rock mass. DDA is used to compute the displacements of the blocks, and DDM is used to predict the crack propagation paths due to the specified boundary conditions. The displacements obtained from DDA are converted into stress and considering Kelvin's solution of the problem the crack propagation mechanism within each block is investigated. Boundary stresses are updated due to crack propagation and new stress boundary conditions in DDA. This cycle continued until crack propagation stopped or a new block formed. Numerical solutions of the experimental rock samples including two random cracks with crack 1 fixed and crack 2 created with different angles and one crack with a slope angle of 30 degrees are compared with the existing experimental and numerical results. This comparison validates the accuracy and effectiveness of the proposed procedure because crack propagation paths predicted are in good agreement with the corresponding experimental results of rock samples.
Mohammadhossein Dehghani Firoozabadi; Mohammad Fatehi Marji; Abolfazl Abdollahipour; Alireza Yarahamdi Bafghi; Yousef Mirzaeian
Abstract
The presence of pores and cracks in porous and fractured rocks is mostly accompanied by fluid flow. Poroelasticity can be used for the accurate modeling of many rock structures in the petroleum industry. The approach of the stress to the value of the fracture stress and the effect of pore pressure on ...
Read More
The presence of pores and cracks in porous and fractured rocks is mostly accompanied by fluid flow. Poroelasticity can be used for the accurate modeling of many rock structures in the petroleum industry. The approach of the stress to the value of the fracture stress and the effect of pore pressure on the deformation of rock are among the effects of fluid on the mechanical behavior of the medium. Due to the deformation-diffusion property of porous media, governing equations, strain-displacement, and stress-strain relationships can be changed to each other. In this study, constitutive equations and relationships necessary to investigate the behavior and reaction of rock in a porous environment are stated. Independent and time-dependent differential equations for an impulse and point fluid source are used to obtain the fundamental solutions. Influence functions are obtained by using the shape functions in the formulation of the fundamental solutions and integrating them. To check the validity and correctness of provided formulation, several examples are mentioned. In the first two examples, numerical application and analytical solution are used at different times and in undrained and drained conditions. In times 0 (undrained response of medium) and 4500 seconds (drained response of medium), there is good coordination and agreement between the numerical and analytical results. In the third example, using the numerical application, a crack propagation path in the wellbore wall is obtained, which is naturally in the direction of maximum horizontal stress.
Mohammadhosein Dehghani Firoozabadi; Mohammad Fatehi Marji; Abolfazl Abdollahipour; Alireza Yarahamdi Bafghi; Yousef Mirzaeian
Abstract
In this work, an effective methodology is introduced for simulation of the crack propagation in linear poroelastic media. The presence of pores and saturated cracks that can be accompanied by fluid flow makes the use of poroelastic media inevitable. In this work, involvement of the time parameter in ...
Read More
In this work, an effective methodology is introduced for simulation of the crack propagation in linear poroelastic media. The presence of pores and saturated cracks that can be accompanied by fluid flow makes the use of poroelastic media inevitable. In this work, involvement of the time parameter in crack propagation is of particular importance. The order of doing the work is such that first, derives the fundamental solutions of a poroelastic higher order displacement discontinuity method (PHODDM). Then will be provided a numerical formulation and implementation for PHODDM in a code named linear element poroelastic DDM (LEP-DDM). Analytical solutions use different times to check the correctness and validity of the proposed solution and the newly developed code. The numerical results show a good agreement and coordination with the analytical results in time zero and 5000 seconds . The code is able to pursue crack-propagation in time and space. This topic is introduced and shown in an example.
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.
R. Alizadeh; M. Fatehi Marji; A. Abdollahipour; M. Pourghasemi Sagand
Abstract
In this work, an effective methodology is introduced for modeling the fatigue crack propagation in linear elastic brittle media. The displacement discontinuity method is used to accomplish the analysis, and the boundaries are discretized with quadratic elements in order to predict the stress intensity ...
Read More
In this work, an effective methodology is introduced for modeling the fatigue crack propagation in linear elastic brittle media. The displacement discontinuity method is used to accomplish the analysis, and the boundaries are discretized with quadratic elements in order to predict the stress intensity factors near the crack tips. This procedure is implemented through 2D linear elastic fracture mechanics. The normal and shear displacement discontinuity around the crack tip is applied to compute the mixed-mode stress intensity factors. The crack growth is incremental, and for each increment of extension, there is no need to use a re-meshing procedure. This method has benefits over the finite element method due to its simplicity in meshing. The crack growth direction is assessed using the maximum principal stress theory. In these analyses, a repetition method is used in order to estimate the correct path of crack propagation. Therefore, the different lengths of incremental growth do not affect the crack growth path analysis. The results are exhibited for several examples with different geometries to demonstrate the efficiency of the approach for analyzing the fatigue crack growth. The accuracy represents that this formulation is ideal for describing the fatigue crack growth problems under the mixed-mode conditions.
H.R. Pasand Masoumi; A. Abdollahipour; Kh. Baghernia
Abstract
Determination of the optimum soil conditioning parameters in the earth pressure balance-tunnel boring machines (EPB-TBMs) plays an important role in reaching an optimum thrust force and advance speed. Silty-clay (CL-ML) in line 1 of the Ahwaz metro project is used in order to find the conditioning parameters ...
Read More
Determination of the optimum soil conditioning parameters in the earth pressure balance-tunnel boring machines (EPB-TBMs) plays an important role in reaching an optimum thrust force and advance speed. Silty-clay (CL-ML) in line 1 of the Ahwaz metro project is used in order to find the conditioning parameters of slumps with different water contents and foam agents. The results obtained are a quantitative comparison between the parameters with different soil conditioning and water contents. Hence, the test results can be used to determine the most economical and technical conditioning parameters for a special condition of soil. The optimum quantity of foam expansion ratio (FER), foam injection ratio (FIR), percent ratio between the surfactant agent and the water volume (Cf), and cost for foam in this soil (based on the soil conditioning production cost) are 10, 157%, 2.07, 248 units, respectively. Soil conditioning with the optimum parameters obtained are tested in a TBM in two stages during excavation of 140 rings. This results in a lower soil conditioning cost and almost 40% higher advance speed.
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.
A. Abdollahipour; M. Fatehi Marji; H. Soltanian; E. A. Kazemzadeh
Abstract
The permeability and coupled behavior of pore pressure and deformations play an important role in hydraulic fracturing (HF) modeling. In this work, a poroelastic displacement discontinuity method is used to study the permeability effect on the HF development in various formation permeabilities. The numerical ...
Read More
The permeability and coupled behavior of pore pressure and deformations play an important role in hydraulic fracturing (HF) modeling. In this work, a poroelastic displacement discontinuity method is used to study the permeability effect on the HF development in various formation permeabilities. The numerical method is verified by the existing analytical and experimental data. Then the propagation of a hydraulic fracture in a formation with a range of permeabilities is studied. The time required for propagation of an HF to 10 times its initial length is used to compare the propagation velocity in the formations with different permeabilities. The results obtained show that the HF propagation can be significantly delayed by a permeability less than almost 10-9 D. Also the effect of HF spacing on the propagation path is studied. It was shown that the stress shadowing effect of HFs remained for a longer spacing than in the elastic model due to the required time for fluid leak-off in the formation. Also the propagation angles are higher in the poroelastic model predictions than the elastic model. Therefore, it is proposed to use the poroelastic model when studying multi-HF propagation in order to avoid errors caused by neglecting the pore fluid effects on the HF propagation paths.
A. Abdollahipour; M. Fatehi Marji; A. R. Yarahmadi Bafghi; J. Gholamnejad
Abstract
Hydraulic fracturing (HF), as a stimulation technique in petroleum engineering, has made possible the oil production from reservoirs with very low permeability. The combination of horizontal drilling and multiple HF with various perforation angles has been widely used to stimulate oil reservoirs for ...
Read More
Hydraulic fracturing (HF), as a stimulation technique in petroleum engineering, has made possible the oil production from reservoirs with very low permeability. The combination of horizontal drilling and multiple HF with various perforation angles has been widely used to stimulate oil reservoirs for economical productions. Despite the wide use of HF, there are still ambiguous aspects that require more investigation. Therefore, optimizing the geometry of the initial fractures using numerical methods is of high importance in a successful HF operation. Different geometrical parameters of the initial HF cracks including patterns, spacings, crack lengths, and perforation phase angles were modeled using the higher order displacement discontinuity method (HODDM) in horizontal and vertical oil wells. Several well-known issues in HF such as crack interference and crack arrest were observed in certain patterns of the HF cracks. Also the best possible arrangements of the HF cracks were determined for a better production. The results obtained were verified by the in-situ measurements existing in the literature. In addition, the best perforation phase angle in vertical wells was investigated and determined.
