scholarly journals Discrete Element Simulation of Interaction between Hydraulic Fracturing and a Single Natural Fracture

Fluids ◽  
2019 ◽  
Vol 4 (2) ◽  
pp. 76
Author(s):  
Basirat ◽  
Goshtasbi ◽  
Ahmadi

Hydraulic fracturing (HF) treatment is performed to enhance the productivity in the fractured reservoirs. During this process, the interaction between HF and natural fracture (NF) plays a critical role by making it possible to predict fracture geometry and reservoir production. In this paper, interaction modes between HF and NF are simulated using the discrete element method (DEM) and effective parameters on the interaction mechanisms are investigated. The numerical results also are compared with different analytical methods and experimental results. The results showed that HF generally tends to cross the NF at an angle of more than 45° and a moderate differential stress (greater than 5 MPa), and the opening mode is dominated at an angle of fewer than 45°. Two effects of changing in the interaction mode and NF opening were also found by changing the strength parameters of NF. Interaction mode was changed by increasing the friction coefficient, while by increasing the cohesion of NF it was less opened under a constant injection pressure.

Author(s):  
Rouhollah Basirat ◽  
Kamran Goshtasbi ◽  
Morteza Ahmadi

Hydraulic Fracturing (HF) is a well-stimulation technique that creates fractures in rock formations through the injection of hydraulically pressurized fluid. Because of the interaction between HF and Natural Fractures (NFs), this process in fractured reservoirs is different from conventional reservoirs. This paper focuses mainly on three effects including anisotropy in the reservoir, strength parameters of discontinuities, and fracture density on HF propagation process using a numerical simulation of Discrete Element Method (DEM). To achieve this aim, a comprehensive study was performed with considering different situations of in situ stress, the presence of a joint set, and different fracture network density in numerical models. The analysis results showed that these factors play a crucial role in HF propagation process. It also was indicated that HF propagation path is not always along the maximum principal stress direction. The results of the numerical models displayed that the affected area under HF treatment is decreased with increasing the strength parameters of natural fracture and decreasing fracture intensity.


2022 ◽  
Author(s):  
Azzan Al-Yaarubi ◽  
Sumaiya Al Bimani ◽  
Sataa Al Rahbi ◽  
Richard Leech ◽  
Dmitrii Smirnov ◽  
...  

Abstract Successful hydraulic fracturing is critical for hydrocarbon recovery from tight reservoirs. Fracture geometry is one essential quality indicator of the created fracture. The geometry provides information about the size of the created fracture and containment and verifies the pre-job modeling. Different techniques are applied to determine fracture geometry, and each has its own advantages and limitations. Due to its simplicity, the radioactive tracer log is commonly used to determine fracture placement and fracture height. Its main drawbacks include shallow depth of investigation, time dependency, and the requirement for multiple interventions for multistage fracturing operations. The crosswell microseismic technique probes a larger volume and it is potentially capable of providing fracture height, length, and orientation. Operational complexity and long processing turnaround time are the main challenges of this technique. Time-lapse shear slowness anisotropy analysis is an effective method to determine hydraulic facture height and orientation. In this technique, the shear slowness anisotropy is recorded before and after the fracture is created. The observed shear anisotropy difference indicates the intervals where the fractures were created, allowing these intervals lengths to be measured. Combining this analysis with gyroscopic data allows determining the fracture orientations. Compared to a tracer log, the differential casedhole sonic anisotropy (DCHSA) has a deeper depth of investigation, and it is time independent. Thus, the repeated log can be acquired at the end of the multistage fracturing operations. Compared to the microseismic technique, this new technique provides more precise fracture height and orientation. The new generation slim dipole sonic technology of 2.125-in. diameter extends the applicability of the DCHSA technique to smaller casing sizes. The shear differential method was applied to a vertical well that targeted the Athel formation in the south of the Sultanate of Oman. This formation is made of silicilyte and is characterized by very low permeability of about 0.01 md on average. Thus, hydraulic fracturing plays a critical role for the economic oil recovery in this reservoir. Aiming to achieve a better zonal contribution, the stimulation design was changed from a limited number of large fractures to an extensive multistage fracturing design in the subject well. Sixteen hydraulic fracturing stages were planned. The DCHSA was applied to provide accurate and efficient fracture geometry evaluation. The DCHSA accurately identified fracture intervals and their corresponding heights and orientations. This enabled effectively determining the created fracture quality and helped explain the responses of the production logs that were recorded during the well test. This study provided a foundation for the placement and completion design of the future wells in the subject reservoir. It particularly revealed adequate fracturing intervals and the optimum number of stages required to achieve optimum reservoir coverage and avoid vertical overlapping.


2013 ◽  
Vol 868 ◽  
pp. 535-541
Author(s):  
Hong Liu ◽  
Lin Wang ◽  
Yu Wu Zhou ◽  
Xi Nan Yu

The fractured low permeability reservoirs develop complex fracture network. As the of waterflooding recovery heightens, excessive high injection pressures and excessive water injection rate will result in open, initiation, propagation and coalescence of micro-fracture, connecting injection with production form the high permeability zone, which results in a one-way onrush of waterflooding, water cut in oil well water rise quickly, causing a severe oil well flooding and channeling, thereby reducing the ultimate oil recovery efficiency. The effect of the waterflooding seepage within natural fracture on fracture initiation is studied and analyzed here, applying the theory of rock fracture mechanics to analyze the interaction of fracture system for naturally fractured reservoirs in waterflooding developing process, studying the mechanical mechanism of opening, initiation, propagation and coalescence of natural fracture under injection pressure, which is important theoretical significance for studying the distribution law of fracture and defining appreciate water injection mode and injection pressure in the process of injection development of the naturally fractured reservoir and for delaying the directivity water break-through and water flooding rate of oil well in the process of injection development.


2018 ◽  
Vol 6 (4) ◽  
pp. T919-T936 ◽  
Author(s):  
Mason K. MacKay ◽  
David W. Eaton ◽  
Per K. Pedersen ◽  
Christopher R. Clarkson

Identifying and characterizing geomechanical domains is important for understanding how a reservoir will respond to hydraulic fracturing, including interaction with natural fractures to create new permeable pathways. We have used a rock-mass characterization approach, which describes the mechanical reservoir package by combining parameters of the intact rock, such as brittleness, with inferred geometry and density of natural fractures. Insights from outcrop observations are important to complement the interpretation of fracture geometry and density derived from subsurface data, to give a more complete understanding of natural fracture networks. This integrated approach is applied to a data set from the Duvernay play in Western Canada. A synthetic model of the subsurface reservoir is constructed using data from well logs, cores, and outcrop analogs. Numerical simulation of the response of the artificial rock mass to hydraulic fracturing is performed using a distinct element code. Independent validation of the model is obtained by achieving an agreement between the simulated microseismic response and the observed distribution of microseismicity during hydraulic fracturing.


Processes ◽  
2020 ◽  
Vol 8 (2) ◽  
pp. 189
Author(s):  
Chen ◽  
Li ◽  
Wu ◽  
Kang

Hydraulic fracturing is a significant technique in petroleum engineering to enhance the production of shale gas or shale oil reservoir. The process of hydraulic fracturing is extremely complicated, referring to the deformation of solid formation, fluid flowing in the crack channel, and coupling the solid with fluid. Simulation of hydraulic fracturing and understanding the course of the mechanism is still a challenging task. In this study, two hydraulic fracturing models, including the Khristianovic–Geertsma–de Klerk (KGD) problem and the hydraulic fracture (HF) intersection with the natural fracture (NF), based on the zero thickness pore pressure cohesive zone (PPCZ) element with contact friction is established. The element can be embedded into the edges of other elements to simulate the fracture initiation and propagation. However, the mesh type of the elements near the PPCZ element has influences on the accuracy and propagation profile. Three common types of mesh, triangle mesh, quadrangle mesh, and deformed quadrangle mesh, are all investigated in this paper. In addition, the infinite boundary condition (IBC) is also discussed in these models. Simulation indicates that the results of pore pressure for zero toughness regime simulated by the triangle mesh are much lower than any others at the early injection time. Secondly, the problem of hydraulic fracturing should be better used with the infinite boundary element (IBE). Moreover, suggestions for crack intersection on the proper mesh type are also given. The conclusions included in this article can be beneficial to research further naturally fractured reservoirs.


2019 ◽  
Vol 59 (1) ◽  
pp. 166
Author(s):  
Mohammad Ali Aghighi ◽  
Raymond Johnson Jr. ◽  
Chris Leonardi

Improved hydraulic fracturing models can better inform operational decisions regarding production from low-permeability coals and ultimately convert currently classified contingent resources to reserves. Improving current modelling approaches requires identification and investigation of the challenges involved in modelling hydraulic fracture stimulation in complex eastern Australian cases where permeability systems and stress regimes can vary significantly. This study investigated differences among existing and emerging advanced hydraulic fracture models and codes including numerical methods used to model fluid and rock behaviours during treatments; the ability to contextualise structure, behaviour and interaction of natural fractures with the propagating hydraulic fracture (e.g. cleat or natural fracture fabric, discrete fracture networks and pressure-dependent leak-off); and their capabilities in handling simultaneously growing or complex fracture development. One finding is that the new generation of models or codes that fully or partially use particle-based numerical methods are more capable in handling complexities associated with hydraulic stimulation of naturally fractured reservoirs. However, the computational cost and time for these models may cause concerns, particularly when modelling large reservoirs and treatments. Based on these limitations, many of the advanced, industry preferred, commercial hydraulic fracture simulators still choose to incorporate limited complexities with regard to natural fractures or represent them mathematically or implicitly. This investigation also indicates that most emerging models provide better representation of natural fractures, visualisation and integration into workflows for completion or stimulation design.


Energies ◽  
2019 ◽  
Vol 12 (20) ◽  
pp. 3852 ◽  
Author(s):  
Kiran Nandlal ◽  
Ruud Weijermars

Hydraulic fracturing for economic production from unconventional reservoirs is subject to many subsurface uncertainties. One such uncertainty is the impact of natural fractures in the vicinity of hydraulic fractures in the reservoir on flow and thus the actual drained rock volume (DRV). We delineate three fundamental processes by which natural fractures can impact flow. Two of these mechanisms are due to the possibility of natural fracture networks to possess (i) enhanced permeability and (ii) enhanced storativity. A systematic approach was used to model the effects of these two mechanisms on flow patterns and drained regions in the reservoir. A third mechanism by which natural fractures may impact reservoir flow is by the reactivation of natural fractures that become extensions of the hydraulic fracture network. The DRV for all three mechanisms can be modeled in flow simulations based on Complex Analysis Methods (CAM), which offer infinite resolution down to a micro-fracture scale, and is thus complementary to numerical simulation methods. In addition to synthetic models, reservoir and natural fracture data from the Hydraulic Fracturing Test Site (Wolfcamp Formation, Midland Basin) were used to determine the real-world impact of natural fractures on drainage patterns in the reservoir. The spatial location and variability in the DRV was more influenced by the natural fracture enhanced permeability than enhanced storativity (related to enhanced porosity). A Carman–Kozeny correlation was used to relate porosity and permeability in the natural fractures. Our study introduces a groundbreaking upscaling procedure for flows with a high number of natural fractures, by combining object-based and flow-based upscaling methods. A key insight is that channeling of flow through natural fractures left undrained areas in the matrix between the fractures. The flow models presented in this study can be implemented to make quick and informed decisions regarding where any undrained volume occurs, which can then be targeted for refracturing. With the method outlined in our study, one can determine the impact and influence of natural fracture sets on the actual drained volume and where the drainage is focused. The DRV analysis of naturally fractured reservoirs will help to better determine the optimum hydraulic fracture design and well spacing to achieve the most efficient recovery rates.


Author(s):  
Nicolas Farah ◽  
Ali Ghadboun

Reservoir simulation is a powerful technique to predict the amount of produced hydrocarbon. After a solid representation of the natural fracture geometry, an accurate simulation model and a physical reservoir model that account for different flow regimes should be developed. Many models based on dual-continuum approaches presented in the literature rely on the Pseudo-Steady-State (PSS) assumption to model the inter-porosity flow. Due to the low permeability in such reservoirs, the transient period could reach several years. Thus, the PSS assumption becomes unjustified. The numerical solution adopted by the Multiple INteracting Continua (MINC) method was able to simulate the transient effects previously overlooked by dual-continuum approaches. However, its accuracy drops with increasing fracture network complexity. A special treatment of the MINC method, i.e., the MINC Proximity Function (MINC–PF) was introduced to address the latter problem. And yet, the MINC–PF suffers a limitation that arises from the existence of several grid-blocks within a studied cell. In this work, this limitation is discussed and two possible solutions (transmissibility recalculation/adjusting the Proximity Function by accounting for nearby fractures) are put forward. Both proposed methods have demonstrated their applicability and effectiveness once compared to a reference solution.


2013 ◽  
Author(s):  
F. Zhang ◽  
Neal Borden Nagel ◽  
Marisela Sanchez-Nagel ◽  
Byungtark Lee ◽  
Alireza Agharazi

Sign in / Sign up

Export Citation Format

Share Document