Coupled Numerical Evaluations of the Geomechanical Interactions Between a Hydraulic Fracture Stimulation and a Natural Fracture System in Shale Formations

2013 ◽  
Vol 46 (3) ◽  
pp. 581-609 ◽  
Author(s):  
N. B. Nagel ◽  
M. A. Sanchez-Nagel ◽  
F. Zhang ◽  
X. Garcia ◽  
B. Lee
SPE Journal ◽  
2018 ◽  
Vol 23 (05) ◽  
pp. 1518-1538 ◽  
Author(s):  
Xiangtong Yang ◽  
Yuanwei Pan ◽  
Wentong Fan ◽  
Yongjie Huang ◽  
Yang Zhang ◽  
...  

Summary The Keshen Reservoir is a naturally fractured, deep, tight sandstone gas reservoir under high tectonic stress. Because the reservoir matrix is very tight, the natural-fracture system is the main pathway for gas production. Meanwhile, stimulation is still required for most production wells to provide production rates that sufficiently compensate for the high cost of drilling and completing wells to access this deep reservoir. Large depletion (and related stress change) was expected during the course of the production of the field. The dynamic response of the reservoir and related risks, such as reduction of fracture conductivity, fault reactivation, and casing failure, would compromise the long-term productivity of the reservoir. To quantify the dynamic response of the reservoir and related risks, a 4D reservoir/geomechanics simulation was conducted for Keshen Reservoir by following an integrated work flow. The work started from systematic laboratory fracture-conductivity tests performed with fractured cores to measure conductivity vs. confining stress for both natural fractures and hydraulic fractures (with proppant placed in the fractures of the core samples). Natural-fracture modeling was conducted to generate a discrete-fracture network (DFN) to delineate spatial distribution of the natural-fracture system. In addition, hydraulic-fracture modeling was conducted to delineate the geometry of the hydraulic-fracture system for the stimulated wells. Then, a 3D geomechanical model was constructed by integrating geological, petrophysical, and geomechanical data, and both the DFN and hydraulic-fracture system were incorporated into the 3D geomechanical model. A 4D reservoir/geomechanics simulation was conducted through coupling with a reservoir simulator to predict variations of stress and strain of rock matrix as well as natural fractures and hydraulic fractures during field production. At each study-well location, a near-wellbore model was extracted from the full-field model, and casing and cement were installed to evaluate well integrity during production. The 4D reservoir/geomechanics simulation revealed that there would be a large reduction of conductivity for both natural fractures and hydraulic fractures, and some fractures with certain dip/dip azimuth will be reactivated during the course of field production. The induced-stress change will also compromise well integrity for those poorly cemented wellbores. The field-development plan must consider all these risks to ensure sustainable long-term production. The paper presents a 4D coupled geomechanics/reservoir-simulation study applied to a high-pressure/high-temperature (HP/HT) naturally fractured reservoir, which has rarely been published previously. The study adapted several new techniques to quantify the mechanical response of both natural fractures and hydraulic fractures, such as using laboratory tests to measure stress sensitivity of natural fractures, integrating DFN and hydraulic-fracture systems into 4D geomechanics simulation, and evaluating well integrity on both the reservoir scale and the near-wellbore scale.


SPE Journal ◽  
2017 ◽  
Vol 22 (06) ◽  
pp. 1834-1855 ◽  
Author(s):  
Cong Xiao ◽  
Leng Tian ◽  
Yayun Zhang ◽  
Tengfei Hou ◽  
Yaokun Yang ◽  
...  

Summary The detection of interacting behavior between the hydraulic fracture (HF) and the natural fracture (NF) is of significance to accurately and efficiently characterize an underground complex-fracture system induced by hydraulic-fracturing technology. This work develops a semianalytical pressure-transient model in the Laplace domain to detect interacting behavior between HF and NF depending on pressure-transient characteristics. Our results have shown that no matter what the flow state (compressible or incompressible flow) within a hydraulically induced fracture system, we can easily detect interacting behavior between HF and NF depending on whether the “dip” shape occurs at the formation radial-flow regime. Referring to sensitivity analysis, distance between NF and well, horizontal distance between NF and HF, and NF length are the three most sensitive factors to detect fracture-interacting behavior. Depending on interference analysis, although the pressure-transient characteristics of a pseudosteady-state dual-porosity model can interfere with our proposed methodology, the difference between our model and a pseudosteady-state dual-porosity system lies in whether the value of the horizontal line of dimensionless pressure derivative is equal to 0.5 at the formation radial-flow regime. Our work obtains some innovative insights into detecting fracture-interacting behavior, and the valuable results can provide significant guidance for refracturing operations and fracture detection in an underground fracture system.


SPE Journal ◽  
2021 ◽  
pp. 1-25
Author(s):  
Wenzheng Liu ◽  
Qingdong Zeng ◽  
Jun Yao

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.


2018 ◽  
pp. 79-92
Author(s):  
A. Akulich ◽  
◽  
Li Kairui ◽  
D. Pestov ◽  
V. Tyurenkova ◽  
...  

2020 ◽  
Vol 10 (8) ◽  
pp. 3333-3345
Author(s):  
Ali Al-Rubaie ◽  
Hisham Khaled Ben Mahmud

Abstract All reservoirs are fractured to some degree. Depending on the density, dimension, orientation and the cementation of natural fractures and the location where the hydraulic fracturing is done, preexisting natural fractures can impact hydraulic fracture propagation and the associated flow capacity. Understanding the interactions between hydraulic fracture and natural fractures is crucial in estimating fracture complexity, stimulated reservoir volume, drained reservoir volume and completion efficiency. However, because of the presence of natural fractures with diffuse penetration and different orientations, the operation is complicated in naturally fractured gas reservoirs. For this purpose, two numerical methods are proposed for simulating the hydraulic fracture in a naturally fractured gas reservoir. However, what hydraulic fracture looks like in the subsurface, especially in unconventional reservoirs, remain elusive, and many times, field observations contradict our common beliefs. In this study, the hydraulic fracture model is considered in terms of the state of tensions, on the interaction between the hydraulic fracture and the natural fracture (45°), and the effect of length and height of hydraulic fracture developed and how to distribute induced stress around the well. In order to determine the direction in which the hydraulic fracture is formed strikethrough, the finite difference method and the individual element for numerical solution are used and simulated. The results indicate that the optimum hydraulic fracture time was when the hydraulic fracture is able to connect natural fractures with large streams and connected to the well, and there is a fundamental difference between the tensile and shear opening. The analysis indicates that the growing hydraulic fracture, the tensile and shear stresses applied to the natural fracture.


Sign in / Sign up

Export Citation Format

Share Document