A homogenization approach for modeling a propagating hydraulic fracture in a layered material

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. MR153-MR162 ◽  
Author(s):  
Egor V. Dontsov

Shales are known to have a finely layered structure, which greatly influences the overall material’s response. Incorporating the effect of all these layers explicitly in a hydraulic fracture simulator would require a prohibitively fine mesh. To avoid such a scenario, a suitable homogenization, which would represent the effect of multiple layers in an average sense, should be performed. We consider a sample variation of elastic properties and minimum horizontal stress versus depth that has more than a hundred layers. We evaluate methodologies to homogenize the stress and the elastic properties. The elastic response of a layered material is found to be equivalent to that of a transversely isotropic material, and the explicit relations for the effective parameters are obtained. To illustrate the relevance of the homogenization procedure for hydraulic fracturing, the propagation of a plane strain hydraulic fracture in a finely layered shale is studied. To reduce the complexity of the numerical model, elastic layering is neglected and only the effect of the stress layers is analyzed. The results demonstrate the ability of the homogenized stress model to accurately capture the hydraulic fracture behavior using a relatively coarse mesh. This result is obtained by using a special asymptotic solution at the tip element that accounts for the local stress variation near the tip, which effectively treats the material at the tip element as nonhomogenized.

2021 ◽  
Vol 9 ◽  
Author(s):  
José Ángel López-Comino ◽  
Simone Cesca ◽  
Peter Niemz ◽  
Torsten Dahm ◽  
Arno Zang

Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at ∼410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter κ, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.


Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. MR137-MR152 ◽  
Author(s):  
Xiaowei Weng ◽  
Dimitry Chuprakov ◽  
Olga Kresse ◽  
Romain Prioul ◽  
Haotian Wang

In laminated formations, the vertical height growth of a hydraulic fracture can be strongly influenced by the interaction of the fracture tip with the bedding interfaces it crosses. A weak interface may fail in shear and then slip, depending on the strength and frictional properties, the effective vertical stress at the interface, and the net pressure. Shear failure and slippage at the interface can retard the height growth or even stop it completely. A 2D analytical model called the FracT model has been developed that examines the shear slippage along the bedding interface adjacent to the fracture tip and the resulting blunting of the fracture tip at the interface, as well as the stress condition on the face opposite from the hydraulic fracture tip for possible fracture nucleation that leads to fracture crossing. The growth of the shear slippage along the interface with time is coupled with the fluid flow into the permeable interface. A parametric study has been carried out to investigate the key formation parameters that influence the crossing/arrest of the fracture at the bedding interface and the shear slippage and depth of fluid penetration into the interface. The study suggests that the interfacial coefficient of friction and the ratio of the vertical to minimum horizontal stress are two of the most influential parameters governing fracture arrest by a weak interface. For the fracture tip to be arrested at the interface, the vertical stress acting on the interface must be close to the minimum horizontal stress or the interfacial coefficient of friction must be very small. The FracT model has also been integrated into a pseudo-3D-based complex hydraulic fracture model. This quantitative mechanistic model that incorporates a bedding-plane slip-driven mechanism is a necessary step to understand and bridge the characterization (sonic) and monitoring (microseismic) observations.


2020 ◽  
Vol 8 (4) ◽  
pp. T1023-T1036
Author(s):  
Cristina Mariana Ruse ◽  
Mehdi Mokhtari

To avoid steep declines in the Tuscaloosa Marine Shale (TMS) production, wells are fracture-stimulated to release the hydrocarbons trapped in the matrix of the formation. An accurate estimation of Young’s modulus and Poisson’s ratio is essential for hydraulic fracture propagation. In addition, ignoring the highly heterogeneous and anisotropic character of TMS can lead to erroneous stress values, which subsequently affect hydraulic fracture width estimates and the overall hydraulic fracturing process. We have developed an empirical 1D geomechanical model that takes into account VTI anisotropy, and it is used to characterize the elastic mechanical properties of TMS in two wells. In the analyzed formation, the vertical Poisson’s ratio is less than the horizontal Poisson’s ratio, which suggests the necessity of an alternative to the ANNIE equations. The stiffness coefficients [Formula: see text] and [Formula: see text] were estimated using the relationships developed from the ultrasonic core data available for the two TMS. Further, correlations between the static and dynamic properties from laboratory tests were used to improve the minimum horizontal stress calculation. We compare VTI Young’s moduli, Poisson’s ratios, and minimum horizontal stress with the isotropic solution. VTI modeling improves the estimation of the elastic mechanical properties. The isotropic solution underestimates the minimum horizontal stress in the formation. Moreover, it was shown that the 20 ft shale interval below the TMS base is characterized by a low Young’s modulus (the vertical Young’s modulus is equal to 20 GPa, whereas the horizontal Young’s modulus is equal to 40 GPa) and may be a frac barrier.


1984 ◽  
Vol 24 (01) ◽  
pp. 19-32 ◽  
Author(s):  
Lawrence W. Teufel ◽  
James A. Clark

Abstract Fracture geometry is an important concern in the design of a massive hydraulic fracture for improved natural gas recovery from low-permeability reservoirs. Determination of the extent of vertical fracture growth and containment in layered rock, a priori, requires an improved understanding of the parameters that may control fracture growth across layer interfaces. We have conducted laboratory hydraulic fracture experiments and elastic finite element studies that show that at least two distinct geologic conditions can inhibit or contain the vertical growth of hydraulic fractures in layered rock:a weak interfacial shear strength of the layers andan increase in the minimum horizontal compressive stress in the bounding layers. The second condition is more important and more likely to occur at depth. Differences in elastic properties within a layered rock mass may be important-not as a containment barrier perse, but in the manner in which variations in elastic properties affect the vertical distribution of the minimum horizontal stress magnitude. These results suggest that improved fracture treatment designs and an assessment of the potential success of stimulations in low-permeability reservoirs can be made by determining the in-situ stress st ate in the producing interval and bounding formations before stimulation. If the bounding formations have a higher minimum horizontal stress, then one can optimize the fracture treatment and maximize the ratio of productive formation fracture area to volume of fluid pumped by limiting bottomhole pressures to that of the bounding formation. Introduction In 1949, Clark introduced the concept of hydraulic fracturing to the petroleum industry. Since then, hydraulic fracture treatment to enhance oil and gas recovery in tight reservoir rocks has become standard practice. More recently, as a result of an increased need for better recovery techniques, massive hydraulic fracturing (MHF) has been used in low-permeability, gas-bearing sandstones in the Rock Mountain region and in Devonian shales of the Appalachian region, where it is uneconomical to retrieve gas in the conventional manner. Massive hydraulic fractures are designed to extend as much as 1000 m (3,281 ft) radially from the wellbore and generally require up to 1000 m3 (6,293 bbl) of fracture fluid. MHF has been developed by trial and error, and the results are uncertain in many situations. Some of these large-scale stimulation efforts have been successful, but others have been extremely disappointing failures. The reasons for these failures are not clear, but it seems likely that improved understanding of the fundamental mechanisms of hydraulic fracturing should suggest ways of improving the efficiency and reliability of the MHF stimulation technique or at least indicate where this technique can be applied successfully. Among the many technological problems encountered in MHF, one of the most important questions that must be answered properly to design a hydraulic fracture treatment for optimal gas recovery concerns the shape and overall geometry of the fracture. The question of fracture height and whether the hydraulic fracture will propagate into formations lying above and below the producing zone. When a fracture treatment is designed, the height of the fracture is the parameter about which the least is known, yet this influences all aspects of the design. A hydraulic fracture usually grows outward in a vertical plane and propagates above and below the packers as well as laterally away from the wellbore. Vertical propagation is undesirable whenever the fracturing is to be contained within a single stratigraphic interval. If the hydraulic fracture is not contained within the producing formation and propagates in both the vertical and lateral directions (an elliptical fracture), failure of the treatment can occur because the fracture fails to contact a sufficiently large area of the reservoir. Moreover, there is an effective loss of the expensive fracture fluid and proppant used to fracture the unproductive formations. An extreme example where the containment of a hydraulic fracture is essential is the case of developing a fracture in a gas-producing sandstone without fracturing through the underlying shale into another sandstone that is water-bearing. Therefore, it is of great economic importance to the gas industry to understand the parameters that can restrict the vertical propagation of massive hydraulic fractures. There are several parameters that are considered to have some effect on the vertical growth and possible containment of hydraulic fractures. SPEJ P. 19^


2020 ◽  
Vol 39 (3) ◽  
pp. 182-187
Author(s):  
Soumen Deshmukh ◽  
Rajesh Sharma ◽  
Manisha Chaudhary ◽  
Harilal

Complex geologic structure, a heterogeneous reservoir, and complications related to high pressure during drilling necessitate carrying out geomechanical modeling to understand the physical properties of rocks and fluids present within the Early Cretaceous synrift sequence in the Bantumilli South area of the Krishna-Godavari Basin in India. Reservoirs within the synrift sequence exhibit low permeability and high pore pressure. Identification of safe mud-weight window zones is critical for safe drilling of wells in this part of the basin. A detailed workflow for building a robust 3D geomechanical model and its applications to well planning and hydraulic fracturing are presented. Elastic properties of the reservoirs were estimated by prestack seismic inversion. Elastic properties and pore pressure volumes were used to simulate the 3D stress field. The maximum horizontal stress direction is observed to be 130°N ± 5°, i.e., northwest to southeast, and estimated fracture pressure (minimum horizontal stress) values range between 10,000 and 14,200 psi within the synrift sequence. The study has shown that the Cretaceous section of the reservoir has narrow mud-weight window zones. These zones are governed mainly by a high pore pressure regime in the reservoirs. Additionally, deep-seated basement faults have played an important role in the compartmentalization of the reservoir in terms of geomechanical properties.


2020 ◽  
Author(s):  
Lingyun Kong ◽  
◽  
Mehdi Ostadhassan ◽  
Junxin Guo

2021 ◽  
Vol 153 ◽  
pp. 103665
Author(s):  
K. Du ◽  
L. Cheng ◽  
J.F. Barthélémy ◽  
I. Sevostianov ◽  
A. Giraud ◽  
...  

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA147-WA155 ◽  
Author(s):  
Marina Pervukhina ◽  
Boris Gurevich ◽  
Pavel Golodoniuc ◽  
David N. Dewhurst

Stress dependency and anisotropy of dynamic elastic properties of shales is important for a number of geophysical applications, including seismic interpretation, fluid identification, and 4D seismic monitoring. Using Sayers-Kachanov formalism, we developed a new model for transversely isotropic (TI) media that describes stress sensitivity behavior of all five elastic coefficients using four physically meaningful parameters. The model is used to parameterize elastic properties of about 20 shales obtained from laboratory measurements and the literature. The four fitting parameters, namely, specific tangential compliance of a single crack, ratio of normal to tangential compliances, characteristic pressure, and crack orientation anisotropy parameter, show moderate to good correlations with the depth from which the shale was extracted. With increasing depth, the tangential compliance exponentially decreases. The crack orientation anisotropy parameter broadly increases with depth for most of the shales, indicating that cracks are getting more aligned in the bedding plane. The ratio of normal to shear compliance and characteristic pressure decreases with depth to 2500 m and then increases below this to 3600 m. The suggested model allows us to evaluate the stress dependency of all five elastic compliances of a TI medium, even if only some of them are known. This may allow the reconstruction of the stress dependency of all five elastic compliances of a shale from log data, for example.


Sign in / Sign up

Export Citation Format

Share Document