Poroelastoplastic Modeling of Complex Hydraulic-Fracture Development in Deep Shale Formations

SPE Journal ◽  
2021 ◽  
pp. 1-25
Author(s):  
Wenzheng Liu ◽  
Qingdong Zeng ◽  
Jun Yao
Keyword(s):  
Adverse Effect ◽  
Hydraulic Fracture ◽  
Cohesive Zone Model ◽  
Picard Iteration ◽  
Natural Fracture ◽  
Rock Properties ◽  
Hydraulic Fractures ◽  
Zone Model ◽  
Shale Formations ◽  
Fracture Development

Summary In this paper, we propose a hydromechanical model to simulate hydraulic fracture propagation in deep shale formations. The Drucker-Prager plasticity theory, Darcy’s law, Reynolds lubrication theory, and Kirchoff’s laws are adopted to describe the plastic deformation, matrix flow, fracture flow, and wellbore flow, respectively. A global embedded cohesive zone model is constructed to achieve the free evolution of hydraulic fractures and the characterization of natural fractures. The finite element method (FEM) and finite volume method (FVM) are used for the spatial discretization of the stress field and pressure field. On the basis of Newton-Raphson iteration, fixed-stress iteration, and Picard iteration, a mixed numerical scheme is built up to solve the strong nonlinear coupling problem. The proposed model is verified against several reference cases and experimental results. Finally, some numerical cases are carried out to investigate the influences of rock properties, natural fracture distribution, and fracturing fluid properties on the complex hydraulic fracture development. The results show that rock plasticity leads to a decrease in stimulated fracture area, an increase in average fracture width, and an increase in propagation pressure. As the cluster number increases, the adverse effect of rock plasticity on multiple hydraulic fracturing in deep shale formations increases significantly. In addition, appropriate optimization of cluster spacing could weaken the adverse effect of rock plasticity on fracturing treatment to a certain extent by using the stress interference effect.

scholarly journals Using Cohesive Zone Model to Simulate the Hydraulic Fracture Interaction with Natural Fracture in Poro-Viscoelastic Formation

Energies ◽  
2019 ◽  
Vol 12 (7) ◽  
pp. 1254 ◽  
Author(s):  
Yu Suo ◽  
Zhixi Chen ◽  
Hao Yan ◽  
Daobing Wang ◽  
Yun Zhang
Keyword(s):  
Numerical Model ◽  
Pore Pressure ◽  
Hydraulic Fracture ◽  
Cohesive Zone ◽  
Cohesive Zone Model ◽  
Fracture Propagation ◽  
Natural Fracture ◽  
Zone Model ◽  
Fluid Viscosity ◽  
Fluid Diffusion

Hydraulic fracturing is a widely used production stimulation technology for conventional and unconventional reservoirs. The cohesive element is used to explain the tip fracture process. In this paper, the cohesive zone model was used to simulate hydraulic fracture initiation and propagation at the same time rock deformation and fluid exchange. A numerical model for fracture propagation in poro-viscoelastic formation is considered. In this numerical model, we incorporate the pore-pressure effect by coupling fluid diffusion with shale matrix viscoelasticity. The numerical procedure for hydraulically driven fracture propagation uses a poro-viscoelasticity theory to describe the fluid diffusion and matrix creep in the solid skeleton, in conjunction with pore-pressure cohesive zone model and ABAQUS was used as a platform for the numerical simulation. The simulation results are compared with the available solutions in the literature. The higher the approaching angle, the higher the differential stress, tensile stress difference, injection rate, and injection fluid viscosity, and it will be easier for hydraulic fracture crossing natural fracture. These results could provide theoretical guidance for predicting the generation of fracture network and gain a better understanding of deformational behavior of shale when fracturing.

A Cohesive Model for Modeling Hydraulic Fractures in Naturally Fractured Formations

2015 ◽  
Author(s):  
M.. Gonzalez ◽  
A. Dahi Taleghani ◽  
J. E. Olson
Keyword(s):  
Fracture Mechanics ◽  
Failure Modes ◽  
Cohesive Zone Model ◽  
Void Nucleation ◽  
Natural Fracture ◽  
Model Parameters ◽  
Hydraulic Fractures ◽  
Zone Model ◽  
Fluid Viscosity ◽  
Naturally Fractured

Abstract A cohesive zone model (CZM) has been developed to couple fluid flow with elastic, plastic and damage behavior of rock during hydraulic fracturing in naturally fractured formations. In addition to inelastic deformations, this model incorporates rock anisotropies. Fracture mechanics of microcrack and micro-void nucleation and their coalescence are incorporated into the formulation of the CZM models to accurately capture different failure modes of rocks. The performance of the developed elastoplastic and CZM models are compared with the available data of a shale play, and then the models are introduced into a commercial finite element package through user-defined subroutines. A workflow to derive the required model parameters for both intact rock and cemented natural fractures is presented through inverse modeling of field data. The hydraulic fractures' growth in the reservoir scale is then simulated, in which the effect of fluid viscosity, natural fracture characteristics and differential stresses on induced fracture network is studied. The simulation results are compared with the available solutions in the literature. The developed CZM model outperforms the traditional fracture mechanics approaches by removing stress singularities at the fracture tips, and simulation of progressive fractures without any essential need for remeshing. This model would provide a robust tool for modeling hydraulic fracture growth using conventional elements of FEA.

The Impacts of Increased Ductility in Organic-Rich Shale on Fracture Network Growth by Hydraulic Fracturing

2021 ◽  
Author(s):  
Chang Huang ◽  
Shengli Chen
Keyword(s):  
Tensile Strength ◽  
Hydraulic Fracturing ◽  
Hydraulic Fracture ◽  
Cohesive Zone ◽  
Cohesive Zone Model ◽  
Fracture Network ◽  
Cohesive Crack ◽  
Hydraulic Fractures ◽  
Zone Model ◽  
Fully Coupled

Abstract The difficulty of hydraulic fracturing in organic-rich shale caused by the increased ductility has not been well interpreted quantitatively, although it is well perceived that the increased shale ductility can impede the propagation of hydraulic fractures and enhance the healing of created fractures upon injection shutdown. This study aims to quantitatively study the impacts of increased ductility on the stimulated reservoir volume (SRV) using an advanced XFEM-based simulator. To achieve this goal, a modified cohesive zone model has been integrated into an in-house fully coupled poroelastic XFEM framework. The study continues by extending the functionality of the numerical framework to simulating multiple interacting fractures. The utilization of the object-oriented programming paradigm in the development of the framework makes it an easy extension to include the multi-fracture network by creating more instances of crack segments. A main hydraulic fracture with an arbitrary number of intersected branches can thus be modeled. A series of parametric studies will be conducted to investigate the impacts of increased ductility on the induced SRV by varying four involved material parameters individually. The modified cohesive zone model, which is essentially a traction-separation law (TSL), is characterized by four parameters: the initial tensile strength Tini, ultimate tensile strength Tkrg, the critical separation Dc, and the final crack separation Dmax. It can flexibly model different crack opening scenarios and simulate more realistically the increased shale ductility. The fully coupled poroelastic XFEM framework has been comprehensively verified against the latest semi-analytical solutions on the four well-known propagation regimes. The numerical results show that the shape of TSL does affect the main hydraulic fracture growth as well as the evolvement of the fracture network, given the same cohesive crack energy and tensile strength. It infers that ductility is not only controlled by cohesive crack energy and tensile strength, which further indicates the necessity of the newly proposed cohesive zone model. The magnitude of the initial tensile strength, controlling when the cohesive crack starts propagating, is found to have the greatest impacts on the fracture length, and SRV, among all four TSL parameters. The novelty of this study is two-fold. First, the newly modified cohesive zone model can more realistically represent the increased shale ductility. Second, the advanced XFEM framework that enables the simulation of a fracture network can study the impacts of increased ductility on the whole SRV but not a single crack.

scholarly journals Numerical Study of Elasto-Plastic Hydraulic Fracture Propagation in Deep Reservoirs Using a Hybrid EDFM–XFEM Method

Energies ◽  
2021 ◽  
Vol 14 (9) ◽  
pp. 2610
Author(s):  
Wenzheng Liu ◽  
Qingdong Zeng ◽  
Jun Yao ◽  
Ziyou Liu ◽  
Tianliang Li ◽  
...  
Keyword(s):  
Plastic Deformation ◽  
Iterative Method ◽  
Hydraulic Fracture ◽  
Iterative Procedure ◽  
Fracture Propagation ◽  
Hydraulic Fractures ◽  
Zone Model ◽  
Dilatancy Angle ◽  
Proposed Model ◽  
Hydraulic Fracture Propagation

Rock yielding may well take place during hydraulic fracturing in deep reservoirs. The prevailing models based on the linear elastic fracture mechanics (LEFM) are incapable of describing the evolution process of hydraulic fractures accurately. In this paper, a hydro-elasto-plastic model is proposed to investigate the hydraulic fracture propagation in deep reservoirs. The Drucker–Prager plasticity model, Darcy’s law, cubic law and cohesive zone model are employed to describe the plastic deformation, matrix flow, fracture flow and evolution of hydraulic fractures, respectively. Combining the embedded discrete fracture model (EDFM), extended finite element method (XFEM) and finite volume method, a hybrid numerical scheme is presented to carry out simulations. A dual-layer iterative procedure is developed based on the fixed-stress split method, Picard iterative method and Newton–Raphson iterative method. The iterative procedure is used to deal with the coupling between nonlinear deformation with fracture extension and fluid flow. The proposed model is verified against analytical solutions and other numerical simulation results. A series of numerical cases are performed to investigate the influences of rock plasticity, internal friction angle, dilatancy angle and permeability on hydraulic fracture propagation. Finally, the proposed model is extended to simulate multiple hydraulic fracture propagation. The result shows that plastic deformation can enhance the stress-shadowing effect.

Predicting Hydraulic Fracturing in a Carbonate Gas Reservoir in Abu Dhabi Using 1D Mechanical Earth Model: Uncertainty and Constraints

2015 ◽  
Author(s):  
Manhal Sirat ◽  
Mujahed Ahmed ◽  
Xing Zhang
Keyword(s):  
Hydraulic Fracturing ◽  
Hydraulic Fracture ◽  
Natural Fracture ◽  
Abu Dhabi ◽  
Hydraulic Fractures ◽  
Gas Reservoir ◽  
In Situ Stresses ◽  
Tight Carbonate ◽  
Different Depths

Abstract In-situ stress state plays an important role in controlling fracture growth and containment in hydraulic fracturing managements. It is evident that the mechanical properties, existing stress regime and the natural fracture network of its reservoir rocks and the surrounding formations mainly control the geometry, size and containments of produced hydraulic fractures. Furthermore, the three principal in situ stresses' axes swap directions and magnitudes at different depths giving rise to identifying different mechanical bedrocks with corresponding stress regimes at different depths. Hence predicting the hydro-fractures can be theoretically achieved once all the above data are available. This is particularly difficult in unconventional and tight carbonate reservoirs, where heterogeneity and highly stress variation, in terms of magnitude and orientation, are expected. To optimize the field development plan (FDP) of a tight carbonate gas reservoir in Abu Dhabi, 1D Mechanical Earth Models (MEMs), involving generating the three principal in-situ stresses' profiles and mechanical property characterization with depth, have been constructed for four vertical wells. The results reveal the swap of stress magnitudes at different mechanical layers, which controls the dimension and orientation of the produced hydro-fractures. Predicted containment of the Hydro-fractures within the specific zones is likely with inevitable high uncertainty when the stress contrast between Sv, SHmax with Shmin respectively as well as Young's modulus and Poisson's Ratio variations cannot be estimated accurately. The uncertainty associated with this analysis is mainly related to the lacking of the calibration of the stress profiles of the 1D MEMs with minifrac and/or XLOT data, and both mechanical and elastic properties with rock mechanic testing results. This study investigates the uncertainty in predicting hydraulic fracture containment due to lacking such calibration, which highlights that a complete suite of data, including calibration of 1D MEMs, is crucial in hydraulic fracture treatment.

Development of a new approach for hydraulic fracturing in tight sand with pre-existing natural fractures

2016 ◽  
Vol 56 (1) ◽  
pp. 225 ◽  
Author(s):  
Kunakorn Pokalai ◽  
David Kulikowski ◽  
Raymond L. Johnson ◽  
Manouchehr Haghighi ◽  
Dennis Cooke
Keyword(s):  
Hydraulic Fracturing ◽  
Hydraulic Fracture ◽  
High Stress ◽  
Gas Production ◽  
Natural Fracture ◽  
Rock Properties ◽  
Fracture Plane ◽  
Natural Fractures ◽  
Geomechanical Model ◽  
Well Design

Hydraulic fracturing in tight gas reservoirs has been performed in the Cooper Basin for decades in reservoirs containing high stress and pre-existing natural fractures, especially near faults. The hydraulic fracture is affected by factors such as tortuosity, high entry pressures, and the rock fabric including natural fractures. These factors cause fracture plane rotation and complexities, leading to fracture disconnection or reduced proppant placement during the treatment. In this paper, rock properties are estimated for a targeted formation using well logs to create a geomechanical model. Natural fracture and stress azimuths within the interval were interpreted from borehole image logs. The image log interpretations inferred that fissures are oriented 30–60° relative to the maximum horizontal stress. Next, diagnostic fracture injection test (DFIT) data was used with the poro-elastic stress equations to predict tectonic strains. Finally, the geomechanical model was history-matched with a planar 3D hydraulic fracturing simulator, and gave more insight into fracture propagation in an environment of pre-existing natural fractures. The natural fracture azimuths and calibrated geomechanical model are input into a framework to evaluate varying scenarios that might result based on a vertical or inclined well design. A well design is proposed based on the natural fracture orientation relative to the hydraulic fracture that minimises complexity to optimise proppant placement. In addition, further models and diagnostics are proposed to aid predicting the hydraulically induced fracture geometry, its impact on gas production, and optimising wellbore trajectory to positively interact with pre-existing natural fractures.

scholarly journals Analysis of hydraulic fracture behavior and well pattern optimization in anisotropic coal reservoirs

2020 ◽  
pp. 014459872096083
Author(s):  
Yulong Liu ◽  
Dazhen Tang ◽  
Hao Xu ◽  
Wei Hou ◽  
Xia Yan
Keyword(s):  
Hydraulic Fracture ◽  
Coalbed Methane ◽  
Fracture Network ◽  
Simulation Software ◽  
Natural Fracture ◽  
Hydraulic Fractures ◽  
Natural Fractures ◽  
Well Pattern ◽  
Well Pattern Optimization ◽  
The Impact

Macrolithotypes control the pore-fracture distribution heterogeneity in coal, which impacts stimulation via hydrofracturing and coalbed methane (CBM) production in the reservoir. Here, the hydraulic fracture was evaluated using the microseismic signal behavior for each macrolithotype with microfracture imaging technology, and the impact of the macrolithotype on hydraulic fracture initiation and propagation was investigated systematically. The result showed that the propagation types of hydraulic fractures are controlled by the macrolithotype. Due to the well-developed natural fracture network, the fracture in the bright coal is more likely to form the “complex fracture network”, and the “simple” case often happens in the dull coal. The hydraulic fracture differences are likely to impact the permeability pathways and the well productivity appears to vary when developing different coal macrolithtypes. Thus, considering the difference of hydraulic fracture and permeability, the CBM productivity characteristics controlled by coal petrology were simulated by numerical simulation software, and the rationality of well pattern optimization factors for each coal macrolithotype was demonstrated. The results showed the square well pattern is more suitable for dull coal and semi-dull coal with undeveloped natural fractures, while diamond and rectangular well pattern is more suitable for semi-bright coal and bright coal with more developed natural fractures and more complex fracturing fracture network; the optimum wells spacing of bright coal and semi-bright coal is 300 m and 250 m, while that of semi-dull coal and dull coal is just 200 m.

Sign in / Sign up

Export Citation Format

Share Document