Integration of Geoscience and Engineering Concepts to Account for Natural Fractures in Fluid Flow within Shale Reservoirs

2021 ◽  
Author(s):  
Clay Kurison ◽  
Ahmed M. Hakami ◽  
Sadi H. Kuleli

Abstract Unconventional shale reservoirs are characterized by low porosity and ultra-low permeability. Natural fractures are known to be present and considered a critical factor for the enhanced post-stimulation productivity. Accounting for natural fractures with existing techniques has not been widely adopted owing to their complexity or lack of validation. Ongoing research efforts are striving to understand how natural fractures can be accounted for and accurately modeled in fluid flow of the subject reservoirs. This study utilized Eagle Ford well data comprising reservoir properties, stimulation metrics, production, microseismicity and permeability measurements from a core plug. The methodology comprised use of production data to extract a linear flow regime parameter. This was coupled with fracture geometry, predicted from hydraulic fracture modeling and microseismicity, to estimate the system permeability. From interpreting microseismic events as slips on critically stressed natural fractures, explicit modeling incorporating a discrete fracture network (DFN) assumed activated natural fractures supplement conductive reservoir contact area. Thus, allowed the estimation of matrix permeability. For validation, the aforementioned was compared with core plug permeability measurements. Results from modeling of planar hydraulic fractures, with microseismicity as validation, predicted planar fracture geometry which when coupled with the linear flow parameter resulted in a system permeability. Incorporation of DFNs to account for activated natural fractures yielded matrix permeability in picodarcy range. A review of laboratory permeability measurements exhibited stress dependence with the value at the maximum experimental confining pressure of 4000 psi in the same range as the computed system permeability. However, the confining pressures used in the experiments were less than the in situ effective stress. Correction for representative stress yielded an ultra-low matrix permeability in the same range as the DFN-based picodarcy matrix permeability. Thus, supporting the adopted drainage architecture and often suggested role of natural fractures in shale reservoir fluid flow. This study presents a multi-discipline workflow to account for natural fractures, and contributes to understanding that will improve laboratory petrophysics and the overall reservoir characterization of the subject reservoirs. Given that the Eagle Ford is an analogue of emerging shales elsewhere, results from this study can be widely adopted.

2021 ◽  
Vol 40 (9) ◽  
pp. 667-676
Author(s):  
Clay Kurison ◽  
Huseyin S. Kuleli

Microseismic events associated with shale reservoir hydraulic fracturing stimulation (HFS) are interpreted to be reactivations of ubiquitous natural fractures (NFs). Despite adoption of discrete fracture network (DFN) models, accounting for NFs in fluid flow within shale reservoirs has remained a challenge. For an explicit account of NFs, this study introduced the use of seismology-based relations linking seismic moment, moment magnitude, fault rupture area, and stress drop. Microseismic data from HFS monitoring of Marcellus Shale horizontal wells had been used to derive planar hydraulic fracture geometry and source properties. The former was integrated with associated well production data found to exhibit transient linear flow. Analytical solutions led to linear flow parameters (LFPs) and system permeability for scenarios depicting flow through infinite and finite conductivity hydraulic fractures. Published core plug permeability was stress-corrected for in-situ conditions to estimate average matrix permeability. For comparison, the burial and thermal history for the study area was used in 1D Darcy-based modeling of steady and episodic expulsion of petroleum to account for geologic timescale persistence of abnormal pore pressure. Both evaluations resulted in matrix permeability in the same picodarcy (pD) range. Coupled with LFPs, reactivated NF surface area for stochastic DFNs was estimated. Subsequently, the aforementioned seismology-based relations were used for determining average stress drops needed to estimate NF rupture area matching flow-based DFN surface areas. Stress drops, comparable to values for tectonic events, were excluded. One of the determined values matched stress drops for HFS operations in past and recent seismological studies. In addition, calculated changes in pore pressure matched estimates in the aforementioned studies. This study unlocked the full potential of microseismic data beyond extraction of planar geometry attributes and stimulated reservoir volume (SRV). Here, microseismic events were explicitly used in the quantitative account of NFs in fluid flow within shale reservoirs.


SPE Journal ◽  
2019 ◽  
Vol 24 (05) ◽  
pp. 2292-2307 ◽  
Author(s):  
Jizhou Tang ◽  
Kan Wu ◽  
Lihua Zuo ◽  
Lizhi Xiao ◽  
Sijie Sun ◽  
...  

Summary Weak bedding planes (BPs) that exist in many tight oil formations and shale–gas formations might strongly affect fracture–height growth during hydraulic–fracturing treatment. Few of the hydraulic–fracture–propagation models developed for unconventional reservoirs are capable of quantitatively estimating the fracture–height containment or predicting the fracture geometry under the influence of multiple BPs. In this paper, we introduce a coupled 3D hydraulic–fracture–propagation model considering the effects of BPs. In this model, a fully 3D displacement–discontinuity method (3D DDM) is used to model the rock deformation. The advantage of this approach is that it addresses both the mechanical interaction between hydraulic fractures and weak BPs in 3D space and the physical mechanism of slippage along weak BPs. Fluid flow governed by a finite–difference methodology considers the flow in both vertical fractures and opening BPs. An iterative algorithm is used to couple fluid flow and rock deformation. Comparison between the developed model and the Perkins–Kern–Nordgren (PKN) model showed good agreement. I–shaped fracture geometry and crossing–shaped fracture geometry were analyzed in this paper. From numerical investigations, we found that BPs cannot be opened if the difference between overburden stress and minimum horizontal stress is large and only shear displacements exist along the BPs, which damage the planes and thus greatly amplify their hydraulic conductivity. Moreover, sensitivity studies investigate the impact on fracture propagation of parameters such as pumping rate (PR), fluid viscosity, and Young's modulus (YM). We investigated the fracture width near the junction between a vertical fracture and the BPs, the latter including the tensile opening of BPs and shear–displacement discontinuities (SDDs) along them. SDDs along BPs increase at the beginning and then decrease at a distance from the junction. The width near the junctions, the opening of BPs, and SDDs along the planes are directly proportional to PR. Because viscosity increases, the width at a junction increases as do the SDDs. YM greatly influences the opening of BPs at a junction and the SDDs along the BPs. This model estimates the fracture–width distribution and the SDDs along the BPs near junctions between the fracture tip and BPs and enables the assessment of the PR required to ensure that the fracture width at junctions and along intersected BPs is sufficient for proppant transport.


2021 ◽  
Author(s):  
Sherif Fakher ◽  
Youssef Elgahawy ◽  
Hesham Abdelaal ◽  
Abdulmohsin Imqam

Abstract Enhanced oil recovery (EOR) in shale reservoirs has been recently shown to increase oil recovery significantly from this unconventional oil and gas source. One of the most studied EOR methods in shale reservoirs is gas injection, with a focus on carbon Dioxide (CO2) mainly due to the ability to both enhance oil recovery and store the CO2 in the formation. Even though several shale plays have reported an increase in oil recovery using CO2 injection, in some cases this method failed severely. This research attempts to investigate the ability of the CO2 to mobilize crude oil from the three most prominent features in the shale reservoirs, including shale matrix, natural fractures, and hydraulically induced fracture. Shale cores with dimensions of 1 inch in diameter and approximately 1.5 inch in length were used in all experiments. The impact of CO2 soaking time and soaking pressure on the oil recovery were studied. The cores were analyzed to understand how and where the CO2 flowed inside the cores and which prominent feature resulted in the increase in oil recovery. Also, a pre-fractured core was used to run an experiment in order to understand the oil recovery potential from fractured reservoirs. Results showed that oil recovery occurred from the shale matrix, stimulation of natural fractures by the CO2, and from the hydraulic fractures with a large volume coming from the stimulated natural fractures. By understanding where the CO2 will most likely be most productive, proper design of the CO2 EOR in shale can be done in order to maximize recovery and avoid complications during injection and production which may lead to severe operational problems.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-17
Author(s):  
Suotang Fu ◽  
Jian Yu ◽  
Kuangsheng Zhang ◽  
Hanbin Liu ◽  
Bing Ma ◽  
...  

Based on the analysis of the typical Ordos well groups, this study began with the accurate characterization of the fracture geometry by adopting advanced laboratory experiment methods and monitoring techniques. Then, with the integration of fracture geometry characterization and in situ stress distributions, fracture optimizations of the target wells were performed through numerical simulations methods. Finally, this study established a sweet spot prediction and identification method for long horizontal shale oil wells and constructed a set of optimization design methods for multistage hydraulic fracturing. This investigation revealed that the hydraulic fractures in Chang-7 terrestrial shale oil reservoirs exhibited the belt pattern, and the primary fractures generated the secondary fractures, which activated the natural fractures and induced shear failure. Macroscopic fractures were found to be perpendicular to the direction of the minimum principal stress. Secondary fractures and activated natural fractures were distributed around the primary fracture in the form of fracture types I and II. Multicluster perforation optimization techniques, which were based on shale reservoir classification and evaluation, and aimed at activating multiclusters and determining fracture sweet spots, were developed. These were successfully applied to the field operation and achieved production enhancement performance.


SPE Journal ◽  
2018 ◽  
Vol 23 (06) ◽  
pp. 2041-2059 ◽  
Author(s):  
Zhiming Chen ◽  
Xinwei Liao ◽  
Kamy Sepehrnoori ◽  
Wei Yu

Summary In this paper, we present an efficient semianalytical model for pressure-transient analysis in fractured wells by considering arbitrarily distributed fracture networks. The semianalytical model included three domains: matrix, hydraulic-fracture networks, and discrete natural fractures. Using the line-source function, we developed the diffusivity equation for fluid flow in matrix. By applying the vertex-analysis technique, we eliminated the flow interplay at fracture intersections and established the diffusivity equations for fluid flow in hydraulic-fracture networks and isolated natural fractures. The pressure-transient solution of these diffusivity equations was obtained using Laplace transforms and the Stehfest numerical inversion. Results showed that with the discrete natural fractures, a “V-shaped” pressure derivative (the classical dual-porosity feature of naturally fractured reservoirs) emerged. With the hydraulic-fracture networks, the reservoir system would exhibit pressure behaviors such as “pseudoboundary-dominated flow,” “fracture-interference flow,” and “fluid-feed flow.” All these pressure characteristics were dependent on the properties and geometries of natural/hydraulic fractures. In addition, through synthetic field application, we found that different (natural/hydraulic) fracture distributions and geometries had distinct behaviors of pressure derivatives, which may provide an effective tool to identify the properties of randomly distributed natural fractures as well as complex hydraulic fractures in unconventional plays.


2016 ◽  
Vol 9 (1) ◽  
pp. 247-256 ◽  
Author(s):  
Ting Li ◽  
Jifang Wan

The application of conventional hydraulic fracture treatment is not ideal in coalbed methane reservoirs, which influences the industry development in China, thus, the present technique should be improved. From two aspects of net pressure and stress sensibility of permeability, it is analyzed and considered that permeability around hydraulic fractures is damaged severely, so this is the main flaw of conventional hydraulic fracturing in CBM. It is proposed to shear natural fractures by fracturing treatment, which are plentiful in coalbed methane reservoirs, and the mechanical condition to generate sheared fractures is presented, in the meanwhile, it is verified that the permeability of sheared fractures is much larger than coal matrix permeability. When the angle between natural and hydraulic fractures is small in coalbed methane reservoirs, the natural fractures will shear easily at low net pressure, so network fractures can be formed. In comparison with conventional hydraulic fracturing, this new methodology can make natural fractures shear at low net pressure to form transverse network fractures, hence, the stimulated reservoir volume is larger, and damage to coal permeability is avoided. This new technique is advantageous in both stimulated reservoir volume and permeability improvement, and it is more adaptable for coalbed methane reservoirs, thus, it has a wide application prospect and significant value.


2015 ◽  
Vol 3 (3) ◽  
pp. SU71-SU88 ◽  
Author(s):  
Yamina E. Aimene ◽  
Ahmed Ouenes

We have developed a new geomechanical workflow to study the mechanics of hydraulic fracturing in naturally fractured unconventional reservoirs. This workflow used the material point method (MPM) for computational mechanics and an equivalent fracture model derived from continuous fracture modeling to represent natural fractures (NFs). We first used the workflow to test the effect of different stress anisotropies on the propagation path of a single NF intersected by a hydraulic fracture. In these elementary studies, increasing the stress anisotropy was found to decrease the curving of a propagating NF, and this could be used to explain the observed trends in the microseismic data. The workflow was applied to Marcellus and Eagle Ford wells, where multiple geomechanical results were validated with microseismic data and tracer tests. Application of the workflow to a Marcellus well provides a strain field that correlates well with microseismicity, and a maximum energy release rate, or [Formula: see text] integral at each completion stage, which appeared to correlate to the production log and could be used to quantify the impact of skipping the completion stages. On the first of two Eagle Ford wells considered, the MPM workflow provided a horizontal differential stress map that showed significant variability imparted by NFs perturbing the regional stress field. Additionally, a map of the strain distribution after stimulating the well showed the same features as the interpreted microseismic data: three distinct regions of microseismic character, supported by tracer tests and explained by the MPM differential stress map. Finally, the workflow was able to estimate, in the second well with no microseismic data, its main performance characteristics as validated by tracer tests. The field-validated MPM geomechanical workflow is a powerful tool for completion optimization in the presence of NFs, which affect in multiple ways the final outcome of hydraulic fracturing.


Sign in / Sign up

Export Citation Format

Share Document