Integration of Dynamic Microseismic Data With a True 3D Modeling of Hydraulic-Fracture Propagation in the Vaca Muerta Shale

SPE Journal ◽  
2017 ◽  
Vol 22 (06) ◽  
pp. 1714-1738 ◽  
Author(s):  
Mahdi Haddad ◽  
Jing Du ◽  
Sandrine Vidal-Gilbert

Summary Microseismic mapping during the hydraulic-fracturing processes in the Vaca Muerta (VM) Shale in Argentina shows a group of microseismic events occurring at shallower depth and at later injection time, and they clearly deviate from the growing planar hydraulic fracture. This spatial and temporal behavior of these shallow microseismic events incurs some questions regarding the nature of these events and their connectivity to the hydraulic fracture. To answer these questions, in this article, we investigate these phenomena by use of a true 3D fracture-propagation-modeling tool along with statistical analysis on the properties of microseismic events. First, we propose a novel technique in Abaqus incorporating fracture intersections in true 3D hydraulic-fracture-propagation simulations by use of a pore-pressure cohesive zone model (CZM), which is validated by comparing our numerical results with the Khristianovic-Geertsma-de Klerk (KGD) solution (Khristianovic and Zheltov 1955; Geertsma and de Klerk 1969). The simulations fully couple slot flow in the fracture with poroelasticity in the matrix and continuum-based leakoff on the fracture walls, and honor the fracture-tip effects in quasibrittle shales. By use of this model, we quantify vertical-natural-fracture activation and fluid infiltration depending on reservoir depth, fracturing-fluid viscosity, mechanical properties of the natural-fracture cohesive layer, natural-fracture conductivity, and horizontal stress contrast. The modeling results demonstrate this natural-fracture activation in coincidence with the hydraulic-fracture-growth complexities at the intersection, such as height throttling, sharp aperture reduction after the intersection, and multibranching at various heights and directions. Finally, we investigate the hydraulic-fracture intersection with a natural fracture in the multilayer VM Shale. We infer the natural-fracture location and orientation from the microseismic-events map and formation microimager log in a nearby vertical well, respectively. We integrate the other field information such as mechanical, geological, and operational data to provide a realistic hydraulic-fracturing simulation in the presence of a natural fracture. Our 3D fracturing simulations equipped with the new fracture-intersection model rigorously simulate the growth of a realistic hydraulic-connection path toward the natural fracture at shallower depths, which was in agreement with our microseismic observations.

Processes ◽  
2019 ◽  
Vol 7 (1) ◽  
pp. 28 ◽  
Author(s):  
Jianxiong Li ◽  
Shiming Dong ◽  
Wen Hua ◽  
Xiaolong Li ◽  
Xin Pan

Complex propagation patterns of hydraulic fractures often play important roles in naturally fractured formations due to complex mechanisms. Therefore, understanding propagation patterns and the geometry of fractures is essential for hydraulic fracturing design. In this work, a seepage–stress–damage coupled model based on the finite pore pressure cohesive zone (PPCZ) method was developed to investigate hydraulic fracture propagation behavior in a naturally fractured reservoir. Compared with the traditional finite element method, the coupled model with global insertion cohesive elements realizes arbitrary propagation of fluid-driven fractures. Numerical simulations of multiple-cluster hydraulic fracturing were carried out to investigate the sensitivities of a multitude of parameters. The results reveal that stress interference from multiple-clusters is responsible for serious suppression and diversion of the fracture network. A lower stress difference benefits the fracture network and helps open natural fractures. By comparing the mechanism of fluid injection, the maximal fracture network can be achieved with various injection rates and viscosities at different fracturing stages. Cluster parameters, including the number of clusters and their spacing, were optimal, satisfying the requirement of creating a large fracture network. These results offer new insights into the propagation pattern of fluid driven fractures and should act as a guide for multiple-cluster hydraulic fracturing, which can help increase the hydraulic fracture volume in naturally fractured reservoirs.


Author(s):  
Minhui Qi ◽  
Mingzhong Li ◽  
Yanchao Li ◽  
Tiankui Guo ◽  
Song Gao

Hydraulic fracturing is an economically effective technology developing the glutenite reservoirs, which have far stronger heterogeneity than the conventional sandstone reservoir. According to the field production experience of Shengli Oilfield, horizontal-well fracturing is more likely to develop a complex fractured network, which improves the stimulated volume of reservoir effectively. But the clear mechanism of horizontal-well hydraulic fracture propagation in the glutenite reservoirs is still not obtained, thus it is difficult to effectively carry out the design of fracturing plan. Based on the characteristics of the glutenite reservoirs, a coupled Flow-Stress-Damage (FSD) model of hydraulic fracture propagation is established. The numerical simulation of fracturing expansion in the horizontal well of the glutenite reservoir is conducted. It is shown that a square mesh-like fracture network is developed near the horizontal well in the reservoir with lower stress difference, in which fracture is more prone to propagate along the direction of the minimum principal stress as well. High fracturing fluids injection displacement and high fracturing fluid viscosity lead to the rise of static pressure of the fracture, which results in the rise of fracture complexity, and greater probability to deflect when encountering gravels. As the perforation density increases, the micro-fractures generated at each perforation gather together faster, and the range of the stimulated reservoir is also relatively large. For reservoirs with high gravel content, the complexity of fracture network and the effect of fracture communication are obviously increased, and the range of fracture deflection is relatively large. In the case of the same gravel distribution, the higher the tensile strength of the gravel, the greater fracture tortuosity and diversion was observed. In this paper, a simulation method of horizontal well fracture network propagation in the reservoirs is introduced, and the result provides the theoretical support for fracture network morphology prediction and plan design of hydraulic fracturing in the glutenite reservoir.


Energies ◽  
2019 ◽  
Vol 12 (7) ◽  
pp. 1254 ◽  
Author(s):  
Yu Suo ◽  
Zhixi Chen ◽  
Hao Yan ◽  
Daobing Wang ◽  
Yun Zhang

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.


2021 ◽  
pp. 014459872110362
Author(s):  
Mingyang Zhai ◽  
Dongying Wang ◽  
Lianchong Li ◽  
Zilin Zhang ◽  
Liaoyuan Zhang ◽  
...  

The tight heterogeneous glutenites are typically characterized by highly variable lithology, low/ultra-low permeability, significant heterogeneity, and a less-developed natural fracture system. It is of great significance for economic development to improve hydraulic fracture complexity and stimulated reservoir volume. To better understand the hydraulic fracturing mechanism, a large-scale experimental test on glutenite specimens was conducted and the hydraulic fracture propagation behaviors and focal mechanism were analyzed. A three-dimensional numerical model was developed to reproduce the hydraulic fracture evolution process and investigate the effects of operating procedures on hydraulic fracture geometry and stimulated reservoir volume. A simultaneous variable injection rate and fluid viscosity technology was proposed to increase the hydraulic fracture complexity and stimulated reservoir volume. The results indicate that four fracturing behaviors can be observed, namely, penetration, deflection, termination, and bifurcating, in the laboratory experiment. Tensile events tend to appear during the initiation stage of hydraulic fracture growth, while shear events and compressive events tend to appear during the non-planar propagation stage. The shear and compressive mechanisms dominate with an increase in the hydraulic fracture complexity. The variable injection rate technology and simultaneous variable injection rate and fluid viscosity technology are effective techniques for fracture geometry control and stimulated reservoir volume enhancement. The key to improve hydraulic fracture complexity is to increase the net pressure in hydraulic fractures, cause evident pressure fluctuations, and activate or communicate a wide range of natural discontinuities. The results can provide a better understanding of the fracture geometry control mechanism in tight heterogeneous glutenites, and offer a guideline for treatment design and optimization of well performance.


2021 ◽  
pp. 014459872198899
Author(s):  
Weiyong Lu ◽  
Changchun He

Directional rupture is one of the most important and most common problems related to rock breaking. The goal of directional rock breaking can be effectively achieved via multi-hole linear co-directional hydraulic fracturing. In this paper, the XSite software was utilized to verify the experimental results of multi-hole linear co-directional hydraulic fracturing., and its basic law is studied. The results indicate that the process of multi-hole linear co-directional hydraulic fracturing can be divided into four stages: water injection boost, hydraulic fracture initiation, and the unstable and stable propagation of hydraulic fracture. The stable expansion stage lasts longer and produces more microcracks than the unstable expansion stage. Due to the existence of the borehole-sealing device, the three-dimensional hydraulic fracture first initiates and expands along the axial direction in the bare borehole section, then extends along the axial direction in the non-bare hole section and finally expands along the axial direction in the rock mass without the borehole. The network formed by hydraulic fracture in rock is not a pure plane, but rather a curved spatial surface. The curved spatial surface passes through both the centre of the borehole and the axial direction relative to the borehole. Due to the boundary effect, the curved spatial surface goes toward the plane in which the maximum principal stress occurs. The local ground stress field is changed due to the initiation and propagation of hydraulic fractures. The propagation direction of the fractures between the fracturing boreholes will be deflected. A fracture propagation pressure that is greater than the minimum principle stress and a tension field that is induced in the leading edge of the fracture end, will aid to fracture intersection; as a result, the possibility of connecting the boreholes will increase.


Energies ◽  
2021 ◽  
Vol 14 (9) ◽  
pp. 2610
Author(s):  
Wenzheng Liu ◽  
Qingdong Zeng ◽  
Jun Yao ◽  
Ziyou Liu ◽  
Tianliang Li ◽  
...  

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.


2021 ◽  
Author(s):  
Somnath Mondal ◽  
Min Zhang ◽  
Paul Huckabee ◽  
Gustavo Ugueto ◽  
Raymond Jones ◽  
...  

Abstract This paper presents advancements in step-down-test (SDT) interpretation to better design perforation clusters. The methods provided here allow us to better estimate the pressure drop in perforations and near-wellbore tortuosity in hydraulic fracturing treatments. Data is presented from field tests from fracturing stages with different completion architectures across multiple basins including Permian Delaware, Vaca Muerta, Montney, and Utica. The sensitivity of near-wellbore pressure drops and perforation size on stimulation distribution effectiveness in plug-and-perf (PnP) treatments is modeled using a coupled hydraulic fracturing simulator. This advanced analysis of SDT data enables us to improve stimulation distribution effectiveness in multi-cluster or multiple entry completions. This analysis goes much further than the methodology presented in URTeC2019-1141 and additional examples are presented to illustrate its advantages. In a typical SDT, the injection flowrate is reduced in four or five abrupt decrements or "steps", each with a duration long enough for the rate and pressure to stabilize. The pressure-rate response is used to estimate the magnitude of perforation efficiency and near-wellbore tortuosity. In this paper, two SDTs with clean fluids were conducted in each stage - one before and another after proppant slurry was injected. SDTs were conducted in cemented single-point entry (cSPE) sleeves, which present a unique opportunity to measure only near-wellbore tortuosity using bottom-hole pressure gauge at sleeve depth, negligible perforation pressure drops, and less uncertainty in interpretation. SDTs were conducted in PnP stages in multiple unconventional basins. The results from one set of PnP stages with optic fiber distributed sensing were modeled with a hydraulic fracturing simulator that combines wellbore proppant transport, perforation size growth, near-wellbore pressure drop, and hydraulic fracture propagation. Past SDT analysis assumed that the pressure drop due to near-wellbore tortuosity is proportional to the flow rate raised to an exponent, β = 0.5, which typically overestimates perforation friction from SDTs. Theoretical derivations show that β is related to the geometry and flow type in the near-wellbore region. Results show that initial β (before proppant slurry) is typically around 0.5, but the final value of β (after proppant slurry) is approximately 1, likely due to the erosion of near-wellbore tortuosity by the proppant slurry. The new methodology incorporates the increase in β due proppant slurry erosion. Hydraulic fracturing modeling, calibrated with optic fiber data, demonstrates that the stimulation distribution effectiveness must consider the interdependence of proppant segregation in the wellbore, perforation erosion, and near-wellbore tortuosity. An improved methodology is presented to quantify the magnitude of perforation and near-wellbore tortuosity related pressure drops before and after pumping of proppant slurry in typical PnP hydraulic fracture stimulations. The workflow presented here shows how the uncertainties in the magnitude of near-wellbore complexity and perforation size, along with uncertainties in hydraulic fracture propagation parameters, can be incorporated in perforation cluster design.


Lithosphere ◽  
2021 ◽  
Vol 2021 (Special 4) ◽  
Author(s):  
Yulong Zhang ◽  
Bei Han ◽  
Xin Zhang ◽  
Yun Jia ◽  
Chun Zhu

Abstract The interaction mode of induced fracture and natural fracture plays an important role in prediction of hydraulic fracture propagation. In this paper, a two-dimensional hydromechanical coupled discrete element model is first introduced in the framework of particle flow simulation, which can well take into account mechanical and hydraulic properties of rock samples with natural fracture. The model’s parameters are strictly calibrated by conducting numerical simulations of uniaxial compression test and direct tensile and shear tests, as well as fluid flow test. The effectiveness of coupled model is also assessed by describing hydraulic fracture propagation in two representative cases, respectively, rock samples with and without preexisting fracture. With this model in hand, the effects of interaction between induced and natural fractures with different approach angles and differential stresses on fluid injection pressure and fracture propagation patterns are investigated and discussed. Results suggest that the interaction modes mainly involve three basic behaviors including the arrested, captured with offset, and directly crossing. For a given differential stress, the captured offset of hydraulic fracture by natural fracture gradually decreases with the approach angle increase, while for a fixed approach angle, that captured offset increases with differential stress decrease.


Sign in / Sign up

Export Citation Format

Share Document