فهرست مطالب امید رضا بارانی

A shear band is a narrow zone of intense shearing strain within a largely unsheared matrix. During the shear band formation, relative tangential displacement of blocks of material on two sides of the band occurs. Width of this region varies from a few to hundreds of microns. Despite the infinitesimal width of the band, its relative tangential deformation may extend to several centimeters and may even lead to exertion of macroscopic influence on the medium.In this paper, a new method for modeling of shear band is presented. A bifurcation analysis is used to detect the onset of localization in an element to determine the geometry of the localized deformation modes. Therefore, a computational procedure is discussed for detecting the onset of localization and determining the localization directions and the shear band path. A zero thickness element is utilized to simulate the slip surface in the shear band. Plane strain condition is assumed in numerical simulations, and elasto-plastic finite element method is used in analysis.When the onset of localization is detected, the zero thickness element is added to the shear band path, which closely reproduces the shear band response. The proposed methodology is applied to two numerical examples, and the results are compared with the other existing methods which demonstrate the ability of the method to resolve the geometry of localized failure modes and reproduce the shear band response. The first example is a simple shear problem that includes a rectangular solid with fixed supports in its bottom, subjected to uniform shear deformation on the upper boundary. The second example is a slope failure problem, which includes a downward displacement applied to the middle of a rigid block on the top of a slope. The shear band path and the force-displacement curves are plotted for both examples and have good agreement with the reported results.

From oil and gas engineering point of view, one of the challenges in low permeable or damaged wells is improving the productivity. There are different methods to increase the productivity of low permeable wells and one of the most efficient one is hydraulic fracturing. In this study, two-dimensional modeling of hydraulic fracturing using finite element method and cohesive element approach through traction-separation law has been performed. This approach avoids the singularity in the crack tip and the cohesive zone fits naturally into the conventional finite element method. Hydraulic fracture is assumed to propagate in a poroelastic and permeable medium with a constant injection rate and under quasi-static conditions and the criterion for fracture initiation is quadratic nominal stress criterion. Also as a propagation criterion, Benzeggagh Kenane (BK) approach has been considered. Two types of elements have been implemented in the model which are 4-node bilinear displacement and pore pressure reduced integration and 6-node displacement and pore pressure two- dimensional cohesive element. Cohesive elements have three degrees of freedom that two of them are in X and Y directions and one of them is pore pressure. Mesh size in the near fracture region is small enough to consider the stress and pressure distribution efficiently and avoid any problem in convergence. Meantime, to decrease the computation cost the mesh size gradually increases from fracture area to the boundaries. Also, to increase the accuracy of the model, the time steps for fracture propagation is 0.01 second. In addition, the effect of fracturing fluid has been directly included in the model which means that the fluid pressure would be applied along the fracture without any simplifying assumption. To validate the model, the results have been compared with KGD approach. The results indicate that in the initial steps the pressure at the wellbore wall is high which decreases with time significantly and eventually it gets a steady and uniform trend. In other words, in the initial steps, the fluid pressure should be high enough to overcome the hoop stress around the wellbore and after some injection periods, the fracturing fluid pressure would reach the breakdown pressure and the fracture starts to initiate and propagate. It is clearly observed that increasing the injection rate would lead to faster propagation of hydraulic fracture and in the models with higher injection rate the fracture tends to grow in the propagation direction. This indirectly means that increasing the injection rate would affect both opening and length of the hydraulic fracture which can result in increasing the productivity. The results reveal that the peak of the normal effective stress profiles corresponds to the fracture tip position, where the fracture opening is zero,and the peak value equals the cohesive strength of the material,as expected.Moreover,with increasing thedistance from the fracture tip,the stress decreases rapidly and approaches the initial stress value. The way that Young’s modulus affects the overall characteristics of hydraulic fracture implies that higher Young’s modulus would lead to longer fractures. In other words, formations with higher Young’s modulus can be fractured easily but the opening of the hydraulic fracture would reduce at the same time. This also indirectly means that Young’s modulus would play an important role in the productivity.

In this paper, a finite element model is developed for the fully hydro-mechanical analysis of hydraulic fracturing in saturated porous media. The model is derived within the framework of generalized Biot theory. The fracture propagation is governed by a cohesive fracture model. The flow within the fracture zone is modeled considering the lubrication equation. In order to describe the fracture in the saturated porous media, momentum equation and mass balance equation with Darcy law are employed. The standard Galerkin method and Newmark scheme are used for discretization in space and time, respectively. Finally, the effects of permeability and rate of injection on the hydraulic fracture propagation are studied. It is observed that an increase in permeability leads to slower crack propagation. In addition, increasing flow rate leads to a faster crack propagation. When permeability increases by 3.3 times, CMOD and crack length decreases by 43.8% and 20%, ,respectively after 1 second and decreases by 29.4% and 15.9%, respectively after 6 seconds. In addition, when flow rate increases by 2, 3, and 4 times, the crack length increases by 30.5%, 55.9%, and 76.3% after one second.