Integration of Core Analysis, Pumping Schedule and Microseismicity to Reduce Uncertainties of Production Performance of Complex Fracture Networks for Multi-Stage Hydraulically Fractured Reservoirs

Author(s):  
Geng Niu ◽  
Jianlei Sun ◽  
Sergei Parsegov ◽  
David Schechter
SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 375-394 ◽  
Author(s):  
Zhiming Chen ◽  
Xinwei Liao ◽  
Wei Yu ◽  
Kamy Sepehrnoori

Summary Fracture networks are extremely important for the management of groundwater, carbon sequestration, and petroleum resources in fractured reservoirs. Numerous efforts have been made to investigate transient behaviors with fracture networks. Unfortunately, because of the complexity and the arbitrary nature of fracture networks, it is still a challenge to study transient behaviors in a computationally efficient manner. In this work, we present a mesh-free approach to investigate transient behaviors in fractured media with complex fracture networks. Contributions of properties and geometries of fracture networks to the transient behaviors were systematically analyzed. The major findings are noted: There are approximately eight transient behaviors in fractured porous media with complex fracture networks. Each behavior has its own special features, which can be used to estimate the fluid front and quantify fracture properties. Geometries of fracture networks have important impacts on the occurrence and the duration of some transient behaviors, which provide a tool to identify the fracture geometries. The fluid production in the fractured porous media is improved with high-conductivity (denser, larger) and high-complexity fracture networks.


SPE Journal ◽  
2017 ◽  
Vol 22 (04) ◽  
pp. 1064-1081 ◽  
Author(s):  
Sanbai Li ◽  
Dongxiao Zhang ◽  
Xiang Li

Summary A fully coupled thermal/hydromechanical (THM) model for hydraulic-fracturing treatments is developed in this study. In this model, the mixed finite-volume/finite-element method is used to solve the coupled system, in which the multipoint flux approximation L-method is used to calculate interelement fluid and heat flux. The Gu et al. (2011) crossing criterion is extended to a 3D scenario to delineate the crossing behaviors as hydraulic fractures meet inclined natural fractures. Moreover, the modified Barton et al. (1985) model proposed by Asadollahi et al. (2010) is used to estimate the fracture aperture and model the shear-dilation effect. After being (partially) verified by means of comparison with results from the literature, the developed model is used to investigate complex-fracture-network propagation in naturally fractured reservoirs. Numerical experiments show that the key factors controlling the complexity of the induced-fracture networks include stress anisotropy, injection rate, natural-fracture distribution (fracture-dip angle, strike angle, spacing, density, and length), fracture-filling properties (the degree of cementation and permeability), fracture-surface properties (cohesion and friction angle), and tensile strength of intact rock. It is found that the smaller the stress anisotropy and/or the lower the injection rate, the more complex the fracture network; a high rock tensile strength could increase the possibility of the occurrence of shear fractures; and under conditions of large permeability of fracture filling combined with small cohesive strength and friction coefficient, shear slip could become the dominant mechanism for generating complex-fracture networks. The model developed and the results presented can be used to understand the propagation of complex-fracture networks and aid in the design and optimization of hydraulic-fracturing treatments.


2021 ◽  
pp. 014459872110417
Author(s):  
Mengmeng Li ◽  
Gang Bi ◽  
Yu Shi ◽  
Kai Zhao

Complex fracture networks are easily developed along the horizontal wellbore during hydraulic fracturing. The water phase increases the seepage resistance of oil in natural fractured reservoir. The flow regimes become more intricate due to the complex fractures and the occurrence of two-phase flow. Therefore, a semi-analytical two-phase flow model is developed based on the assumption of orthogonal fracture networks to describe the complicate flow regimes. The natural micro-fractures are treated as a dual-porosity system and the hydraulic fracture with complex fracture networks are characterized explicitly by discretizing the fracture networks into multiple fracture segments. The model is solved according to Laplace transformation and Duhamel superposition principle. Results show that seven possible flow regimes are described according to the typical curves. The major difference between the vertical fractures and the fracture networks along the horizontal wellbore is the fluid “feed flow” behavior from the secondary fracture to the main fracture. A natural fracture pseudo-radial flow stage is added in the proposed model comparing with the conventional dual-porosity model. The water content has a major effect on the fluid total mobility and flow capacity in dual-porosity system and complex fracture networks. With the increase of the main fracture number, the interference of the fractures increases and the linear flow characteristics in the fracture become more obvious. The secondary fracture number has major influence on the fluid feed capacity from the secondary fracture to the main fracture. The elastic storativity ratio mainly influences the fracture flow period and inter-porosity flow period in the dual-porosity system. The inter-porosity flow coefficient corresponds to the inter-porosity flow period of the pressure curves. This work is significantly important for the hydraulic fracture characterization and performance prediction of the fractured horizontal well with complex fracture networks in natural fractured reservoirs.


2015 ◽  
Vol 18 (04) ◽  
pp. 463-480 ◽  
Author(s):  
Jianlei Sun ◽  
David Schechter

Summary Multistage hydraulically fractured wells are applied widely to produce unconventional resource plays. In naturally fractured reservoirs, hydraulic-fracture treatments may induce complex-fracture geometries that one cannot model accurately and efficiently with Cartesian and corner-point grid systems or standard dual-porosity approaches. The interaction of hydraulic and naturally occurring fractures almost certainly plays a role in ultimate well and reservoir performance. Current simulation models are unable to capture the complexity of this interaction. Generally speaking, our ability to detect and characterize fracture systems is far beyond our capability of modeling complex natural-fracture systems. To evaluate production performance in these complex settings with numerical simulation, fracture networks require advanced meshing and domain-discretization techniques. This paper investigates these issues by developing natural-fracture networks with fractal-based techniques. After a fracture network is developed, we demonstrate the feasibility of gridding complex natural-fracture behavior with optimization-based unstructured meshing algorithms. Then we can demonstrate that one can simulate natural-fracture complexities such as variable aperture, spacing, length, and strike. This new approach is a significant step beyond the current method of dual-porosity simulation that essentially negates the sophisticated level of fracture characterization pursued by many operators. We use currently established code for fractal discrete-fracture-network (FDFN) models to build realizations of naturally fractured reservoirs in terms of stochastic fracture networks. From outcrop, image-log, and core analysis, it is possible to extract fracture fractal parameters pertaining to aperture, spacing, and length distribution, including center distribution as well as a fracture strike. Then these parameters are used as input variables for the FDFN code to generate multiple realizations of fracture networks mimicking fracture clustering and randomly distributed natural fractures. After incorporating hydraulic fractures, complex-fracture networks are obtained for further reservoir-domain discretization. To discretize the complex-fracture networks, a new mesh-generation approach is developed to conform to nonorthogonal and low-angle intersections of extensively clustered discrete-fracture networks with nonuniform aperture distribution. Optimization algorithms are adopted to reduce highly skewed cells, and to ensure good mesh quality around fracture tips, intersections, and regions of extensive fracture clustering. Moreover, local grid refinement is implemented with a predefined distance function to control cell sizes and shapes around and far away from fractures. Natural-fracture spacing, length, strike, and aperture distribution are explicitly gridded, thus introducing a new simulation approach that is far superior to dual-porosity simulation. Finally, initial sensitivity studies are performed to demonstrate both the capability of the optimization-based unstructured meshing algorithms, and the effect of aforementioned natural-fracture parameters on well performance. This study demonstrates how to incorporate a fractal-based characterization approach into the current work flow for simulating unconventional reservoirs, and most importantly solves several issues such as nonorthogonal intersections, extensive clustering, and nonuniform aperture distribution associated with domain discretization with unstructured grids for complex-fracture networks. The proposed meshing techniques for complex fracture networks can be easily implemented in existing preprocessing, unstructured mesh generators. The sensitivity study and the simulation runs demonstrate the importance of fracture characterization as well as uncertainties associated with naturally fractured reservoirs on well-production performance.


SPE Journal ◽  
2016 ◽  
Vol 21 (02) ◽  
pp. 538-549 ◽  
Author(s):  
Zhiming Chen ◽  
Xinwei Liao ◽  
Xiaoliang Zhao ◽  
Sanbo Lv ◽  
Langtao Zhu

Summary In naturally fractured reservoirs, complex fracture systems can easily develop along a horizontal wellbore during hydraulic fracturing. In the fracture systems, multiple, discrete secondary fractures are connected to the multiple-fractured horizontal well (MFHW). Because of the fracture complexity, most studies about performance forecast of such MFHWs highly depend on numerical simulators. In this paper, a new semianalytical approach is proposed to overcome the challenge to analyze the pressure behavior of MFHWs in complex-fracture systems. First, a mathematical model for MFHWs with secondary-fracture networks is established. Then, with Gauss elimination and the Stehfest numerical algorithm (Stehfest 1970), the transient-pressure solution of the mathematical model is solved, and type curves of MFHWs with secondary-fracture networks are obtained. After that, model validation and sensitivity analysis are conducted. It is found that the presented approach can rapidly and accurately generate type curves of MFHWs with secondary-fracture networks. This work provides very meaningful references for reservoir engineers in fracturing evaluations as well as performance estimations of MFHWs in naturally fractured reservoirs.


2018 ◽  
Vol 140 (7) ◽  
Author(s):  
Youshi Jiang ◽  
Arash Dahi-Taleghani

Fluid flow in fractured porous media has always been important in different engineering applications especially in hydrology and reservoir engineering. However, by the onset of the hydraulic fracturing revolution, massive fracturing jobs have been implemented in unconventional hydrocarbon resources such as tight gas and shale gas reservoirs that make understanding fluid flow in fractured media more significant. Considering ultralow permeability of these reservoirs, induced complex fracture networks play a significant role in economic production of these resources. Hence, having a robust and fast numerical technique to evaluate flow through complex fracture networks can play a crucial role in the progress of inversion methods to determine fracture geometries in the subsurface. Current methods for tight gas flow in fractured reservoirs, despite their advantages, still have several shortcomings that make their application for real field problems limited. For instance, the dual permeability theory assumes an ideal uniform orthogonal distribution of fractures, which is quite different from field observation; on the other hand, numerical methods like discrete fracture network (DFN) models can portray the irregular distribution of fractures, but requires massive mesh refinements to have the fractures aligned with the grid/element edges, which can greatly increase the computational cost and simulation time. This paper combines the extended finite element methods (XFEM) and the gas pseudo-pressure to simulate gas flow in fractured tight gas reservoirs by incorporating the strong-discontinuity enrichment scheme to capture the weak-discontinuity feature induced by highly permeable fractures. Utilizing pseudo-pressure formulations simplifies the governing equations and reduces the nonlinearity of the problem significantly. This technique can consider multiple fracture sets and their intersection to mimic real fracture networks on a plain structured mesh. Here, we utilize the unified Hagen–Poiseuille-type equation to compute the permeability of tight gas, and finally adopt Newton–Raphson iteration method to solve the highly nonlinear equations. Numerical results illustrate that XFEM is considerably effective in fast calculation of gas flow in fractured porous media.


2021 ◽  
pp. 1-15
Author(s):  
Youwei He ◽  
Yingjie Xu ◽  
Yong Tang ◽  
Yu Qiao ◽  
Wei Yu ◽  
...  

Abstract Complex fracture networks (CFN) provide flow channels and significantly affect well performance in unconventional reservoirs. However, traditional rate transient analysis (RTA) models barely consider the effect of CFN on production performance. The impact of multi-phase flow on rate transient behaviors is still unclear especially under CFN. Neglecting these effects could cause incorrect rate transient response and erroneous estimation of well and fracture parameters. This paper investigates multi-phase rate transient behaviors considering CFN, and tries to investigate in what situations the multi-phase models should be used to obtain more accurate results. Firstly, an embedded discrete fracture model (EDFM) is generated instead of LGR method to overcome time-intensive computation. The model is coupled with reservoir models using non-neighboring connections (NNCs). Secondly, eight cases are designed using the EDFM technology to analyze effect of natural fractures, formation permeability, and relative permeability on rate transient behaviors. Thirdly, Blasingame plot, log-log plot, and linear flow plot are used to analyze the differences of rate transient response between single-phase and multi-phase flow in reservoirs with CFN. For multi-phase flow, severe deviations can be observed on RTA plots compared with single-phase model. Combination of three RTA type curves can characterize the differences from early to late flow regimes and improve the interpretation accuracy as well as reduce the non-unicity. Finally, field data analysis in Permian Basin demonstrates that multi-phase RTA analysis are required for analyzing production and pressure data since single-phase RTA analysis will lead to big errors especially under high water cut during fracturing fluid flowback period, early production of unconventional gas wells or after waterflooding or water huff-n-puff.


Sign in / Sign up

Export Citation Format

Share Document