A Case Study of Vertical Hydraulic Fracture Growth, Stress Variations with Depth and Shear Stimulation in the Niobrara Shale and Codell Sand, DJ Basin, Colorado

2021 ◽  
pp. 1-34
Author(s):  
Kevin L. McCormack ◽  
Mark D. Zoback ◽  
Wenhuan Kuang

We carried out a geomechanical study of three wells, one each in the Niobrara A, Niobrara C and Codell sandstone to investigate how the state of stress and stress variations with depth affect vertical hydraulic fracture growth and shear stimulation of pre-existing fractures. We demonstrate that the higher magnitudes of measured least principal stress values in the Niobrara A and C shales are the result of viscoplastic stress relaxation. Using a density log and a VTI velocity model developed to accurately locate the microseismic events, we theoretically calculated a continuous profile of the magnitude of the least principal stress with depth. This stress profile explains the apparent vertical hydraulic fracture growth as inferred from the well-constrained depths of associated microseismic events. Finally, we demonstrate that because of the upward propagation of hydraulic fractures from the Niobrara C to the Niobrara A, the latter formation experienced considerably more shear stimulation, which may contribute to the greater production of oil and gas from that formation.

SPE Journal ◽  
2018 ◽  
Vol 23 (06) ◽  
pp. 2118-2132 ◽  
Author(s):  
Di Wang ◽  
Mian Chen ◽  
Yan Jin ◽  
Andrew. P. Bunger

Summary Hydraulic fracturing using supercritical carbon dioxide (CO2) has a recognized potential to grow in importance for unconventional oil and gas reservoirs. It is characterized by higher compressibility than traditional liquid-phase hydraulic-fracturing fluids. Motivated by the larger compressibility of supercritical CO2, this paper considers the problem of a hydraulic fracture in which a compressible fluid is injected at a constant rate to drive a hydraulic fracture in a permeable and brittle rock. The two cases of a plane-strain fracture and a penny-shaped fracture are considered. It is shown that for many practical cases, the formation has a large enough fracture toughness that the propagation is in a regime for which the pressure inside the hydraulic fracture can be treated as spatially uniform (“toughness dominated”). Both numerical simulations and analytical solutions for the relevant limiting regimes show that fluid compressibility affects fracture shape only at the very beginning period, which corresponds to the storage regime, and has little effect on fracture growth in the leakoff regime. Overall, because the transition from the storage regime to the leakoff regime is expected to often take place in a short time after the fracture starts propagating, the influence of compressibility in the storage regime is very brief and can be quickly ignored. Therefore, even relatively sizable fluid compressibility has almost no effect on fracture growth in the toughness-dominated regime when leakoff is taken into account.


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

During horizontal well staged fracturing, there is stress interference between multiple transverse fractures in the same perforation cluster. Theoretical analysis and numerical calculation methods are applied in this study. We analysed the mechanism of induced stress interference in a single fracture under different fracture spacings and principal stress ratios. We also investigated the hydraulic fracture morphology and synchronous expansion process under different fracture spacings and principal stress ratios. The results show that the essence of induced stress is the stress increment in the area around the hydraulic fracture. Induced stress had a dual role in the fracturing process. It created favourable ground stress conditions for the diversion of hydraulic fractures and the formation of complex fracture network systems, inhibited fracture expansion in local areas, stopped hydraulic fractures, and prevented the formation of effective fractures. The curves of the maximum principal stress, minimum principal stress, and induced principal stress difference with distance under different fracture lengths, different fracture spacings, and different principal stress ratios were consistent overall. With a small fracture spacing and a small principal stress ratio, intermediate hydraulic fractures were difficult to initiate or arrest soon after initiation, fractures did not expand easily, and the expansion speed of lateral hydraulic fractures was fast. Moreover, with a smaller fracture spacing and a smaller principal stress ratio, hydraulic fractures were more prone to steering, and even new fractures were produced in the minimum principal stress direction, which was beneficial to the fracture network communication in the reservoir. When the local stress and fracture spacing were appropriate, the intermediate fracture could expand normally, which could effectively increase the reservoir permeability.


Geophysics ◽  
2022 ◽  
pp. 1-56
Author(s):  
Ankush Singh ◽  
Mark D. Zoback

Knowledge of layer-to-layer variations of the least principal stress, S hmin, with depth is essential for optimization of multi-stage hydraulic fracturing in unconventional reservoirs. Utilizing a geomechanical model based on viscoelastic stress relaxation in relatively clay rich rocks, we present a new method for predicting continuous S hmin variations with depth. The method utilizes geophysical log data and S hmin measurements from routine diagnostic fracture injection tests (DFITs) at several depths for calibration. We consider a case study in the Wolfcamp formation in the Midland Basin, where both geophysical logs and values of S hmin from DFITs are available. We compute a continuous stress profile as a function of the well logs that fits all of the DFITs well. We utilized several machine learning technologies, such as bootstrap aggregation (or bagging), to improve the generalization of the model and demonstrate that the excellent fit between predicted and observed stress values is not the result of over-fitting the calibration points. The model is then validated by accurately predicting hold-out stress measurements from four wells within the study area and, without recalibration, accurately predicting stress as a function of depth in an offset pad about 6 miles away.


1982 ◽  
Vol 22 (03) ◽  
pp. 333-340 ◽  
Author(s):  
Norman R. Warpinski ◽  
James A. Clark ◽  
Richard A. Schmidt ◽  
Clarence W. Huddle

Abstract Laboratory experiments have been conducted to determine the effect of in-situ stress variations on hydraulic fracture containment. Fractures were initiated in layered rock samples with prescribed stress variations, and fracture growth characteristics were determined as a function of stress levels. Stress contrasts of 300 to 400 psi (2 to 3 MPa) were found sufficient to restrict fracture growth in laboratory samples of Nevada tuff and Tennessee and Nugget sandstones. The required stress level was found not to depend on mechanical rock properties. However, permeability and the resultant pore pressure effects were important. Tests conducted at biomaterial interfaces between Nugget and Tennessee sandstones show that the resultant stresses set up near the interface because of the applied overburden stress affect the fracture behavior in the same way as the applied confining stresses. These results provide a guideline for determining the in-situ stress contrast necessary to contain a fracture in a field treatment. Introduction An under-standing of the factors that influence and control hydraulic fracture containment is important for the successful use of hydraulic fracturing technology in the enhanced production of natural gas from tight reservoirs. Optimally, this understanding would provide improved fracture design criteria to maximize fracture surface area in contact with the reservoir with respect to volume injected and other treatment parameters. In formations with a positive containment condition (i.e., where fracturing out of zone is not anticipated), long penetrating fractures could be used effectively to develop the resource. For the opposite case, the options would beto use a small treatment so that large volumes are not wasted in out-of-zone fracturing and to accept a lower productivity improvement, orto reject the zone as uneconomical. These decisions cannot be made satisfactorily unless criteria for vertical fracture propagation are developed and techniques for readily measuring the important parameters are available. Currently, both theoretical and experimental efforts are being pursued to determine the important parameters and their relative effects on fracture growth. Two modes of fracture containment are possible. One is the situation where fracture growth is terminated at a discrete interface. Examples of this include laboratory experiments showing fracture termination at weak or unbonded interfaces and theoretical models that predict that fracture growth will terminate at a material property interface. The other mode may occur when the fracture propagates into the bounding layer, but extensive growth does not take place and the fracture thus is restricted. An example is the propagation of the fracture into a region with an adverse stress gradient so that continued propagation results in higher stresses on the fracture, which limits growth, as suggested by Simonson et al. and as seen in mineback experiments. Another example is the possible restriction caused by propagation into a higher modulus region where the decreased width results in increased pressure drop in the fracture, which might inhibit extensive growth into that region relative to the lower modulus region. Other parameters, such as natural fractures, treatment parameters, pore pressure, etc., may affect either of these modes. Laboratory and mineback experiments have shown that weak interfaces and in-situ stress differences are the most likely factors to contain the fracture, and weak interfaces are probably effective only at shallow depths. Thus, our experiments are being performed to determine the effect of in-situ stresses on fracture containment, both in a uniform rock sample and at material properly interfaces. SPEJ P. 333^


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-12 ◽  
Author(s):  
Song Yaobin ◽  
Lu Weiyong ◽  
He Changchun ◽  
Bai Erhu

According to the theory of plane mechanics involving the interaction of hydraulic and natural fractures, the law of hydraulic fracture propagation under the influence of natural fractures is verified using theoretical analysis and RFPA2D-Flow numerical simulation approaches. The shear and tensile failure mechanisms of rock are simultaneously considered. Furthermore, the effects of the approach angle, principal stress difference, tensile strength and length of the natural fracture, and elastic modulus and Poisson’s ratio of the reservoir on the propagation law of a hydraulic fracture are investigated. The following results are obtained: (1) The numerical results agree with the experimental data, indicating that the RFPA2D-Flow software can be used to examine the hydraulic fracture propagation process under the action of natural fractures. (2) In the case of a low principal stress difference and low approach angle, the hydraulic fracture likely causes shear failure along the tip of the natural fracture. However, under a high stress difference and high approach angle, the hydraulic fracture spreads directly through the natural fracture along the original direction. (3) When natural fractures with a low tensile strength encounter hydraulic fractures, the hydraulic fractures likely deviate and expand along the natural fractures. However, in the case of natural fractures with a high tensile strength, the natural fracture surface is closed, and the hydraulic fracture directly passes through the natural fracture, propagating along the direction of the maximum principal stress. (4) Under the same principal stress difference, a longer natural fracture corresponds to the easier initiation and expansion of a hydraulic fracture from the tip of the natural fracture. However, when the size of the natural fracture is small, the hydraulic fracture tends to propagate directly through the natural fracture. (5) A smaller elastic modulus and larger Poisson’s ratio of the reservoir result in a larger fracture initiation pressure. The presented findings can provide theoretical guidance regarding the hydraulic fracturing of reservoirs with natural fractures.


2020 ◽  
Vol 38 (6) ◽  
pp. 2466-2484
Author(s):  
Jianguang Wei ◽  
Saipeng Huang ◽  
Guangwei Hao ◽  
Jiangtao Li ◽  
Xiaofeng Zhou ◽  
...  

Hydraulic fracture initiation and propagation are extremely important on deciding the production capacity and are crucial for oil and gas exploration and development. Based on a self-designed system, multi-perforation cluster-staged fracturing in thick tight sandstone reservoir was simulated in the laboratory. Moreover, the technology of staged fracturing during casing completion was achieved by using a preformed perforated wellbore. Three hydraulic fracturing methods, including single-perforation cluster fracturing, multi-perforation cluster conventional fracturing and multi-perforation cluster staged fracturing, were applied and studied, respectively. The results clearly indicate that the hydraulic fractures resulting from single-perforation cluster fracturing are relatively simple, which is difficult to form fracture network. In contrast, multi-perforation cluster-staged fracturing has more probability to produce complex fractures including major fracture and its branched fractures, especially in heterogeneous samples. Furthermore, the propagation direction of hydraulic fractures tends to change in heterogeneous samples, which is more likely to form a multi-directional hydraulic fracture network. The fracture area is greatly increased when the perforation cluster density increases in multi-perforation cluster conventional fracturing and multi-perforation cluster-staged fracturing. Moreover, higher perforation cluster densities and larger stage numbers are beneficial to hydraulic fracture initiation. The breakdown pressure in homogeneous samples is much higher than that in heterogeneous samples during hydraulic fracturing. In addition, the time of first fracture initiation has the trend that the shorter the initiation time is, the higher the breakdown pressure is. The results of this study provide meaningful suggestions for enhancing the production mechanism of multi-perforation cluster staged fracturing.


2022 ◽  
Author(s):  
Aamir Lokhandwala ◽  
Vaibhav Joshi ◽  
Ankit Dutt

Abstract Reservoir simulation is used in most modern reservoir studies to predict future production of oil and gas, and to plan the development of the reservoir. The number of hydraulically fractured wells has risen drastically in recent years due to the increase in production in unconventional reservoirs. Gone are the days of using simple analytic techniques to forecast the production of a hydraulic fracture in a vertical well, and the need to be able to model multiple hydraulic fractures in many stages over long horizontals is now a common practice. The type of simulation approach chosen depends on many factors and is study specific. Pseudo well connection approach was preferred in the current case. Due to the nature of the reservoir simulation problem, a decision needs to be made to determine which hydraulic fracture modeling method might be most suitable for any given study. To do this, a selection of methods is chosen based on what is available at hand, and what is commonly used in various reservoir simulation software packages. The pseudo well connection method, which models hydraulic fractures as uniform conductivity rectangular fractures was utilized for a field of interest referred to as Field A in this paper. Such an assumption of the nature of the hydraulic fracture is common in most modern tools. Field A is a low permeability (0.01md-0.1md), tight (8% to 12% porosity) gas-condensate (API ~51deg and CGR~65 stb/mmscf) reservoir at ~3000m depth. Being structurally complex, it has a large number of erosional features and pinch-outs. The pseudo well connection approach was found to be efficient both terms of replicating data of Field A for a 10 year period while drastically reducing simulation runtime for the subsequent 10 year-period too. It helped the subsurface team to test multiple scenarios in a limited time-frame leading to improved project management.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. KS159-KS168 ◽  
Author(s):  
Yunhui Tan ◽  
Terry Engelder

Moment tensors are calculated by using the P-wave first motion peak amplitudes of 59 microseismic events with high signal-to-noise ratio. These events are from a surface microseismic data set gathered during hydraulic-fracture stimulation of the Marcellus gas shale in Washington County, Pennsylvania, USA. The majority of these 59 events have a horizontal nodal plane ([Formula: see text] a few degrees) characteristic of a dip-slip/horizontal-slip moment tensor. If the horizontal nodal plane is an auxiliary, the vertical nodal plane has a pure dip-slip motion, which is inconsistent with the opening motion for vertical hydraulic fractures that enables proppant loading. This points to slip on horizontal nodal planes with the auxiliary vertical nodal planes aligned with the local maximum horizontal stress orientation as indicated by drilling-induced fractures in nearby vertical wells. These 59 microseismic events are caused by slip on horizontal mechanical discontinuities such as bedding planes during the opening of vertical hydraulic fractures, a model first proposed by research teams headed by Rutledge and Eisner, respectively. During several stimulation stages in the Washington County Marcellus gas shale, a pattern of opposite slip direction develops within “double lineaments” of microseismic clouds. This suggests that fracking fluid is not only able to move in the direction of fracture propagation, but it can also spread sideways into previously unstimulated rock. A secondary microseismic cloud consistently initiates at approximately 133 m (400 ft) from the position opposite the central perforation toward the unstimulated heel of the horizontal wells. From these moment tensors, we have concluded that microseismic focal mechanisms with horizontal nodal planes are direct evidence of the presence of treatment fluid in open hydraulic fractures.


Author(s):  
Hai T. Nguyen ◽  
Jang Hyun Lee ◽  
Khaled A. Elraies

AbstractIn the field of hydraulic fracture modeling, the pseudo-three-dimensional (P3D) approach is an efficient and practical computational tool serving as a compromise between two-dimensional and planar three-dimensional models. This review discusses the P3D modeling approach from its early developmental stage in the 1980s to the present. The evolution of P3D modeling is drawn over time based on the major differences in the governing formulation and assumptions considered by each model. The problems of equilibrium height growth and vertical viscous fluid resistance (i.e., non-equilibrium height growth) emphasize the primary differences among these models. Besides, the P3D-based complex fracture network models for shale oil and gas reservoirs accounting for the interaction between preexisting natural fractures and induced hydraulic fractures are discussed. Finally, in the application section, several simulations are reported to demonstrate the validation of the P3D numerical algorithm by comparing it with the Perkins–Kern–Nordgren (PKN) large and small asymptotic solutions, as well as the effect of time-dependent variable injection rates on the hydraulic fracture propagation. The results showed a good matching between P3D and PKN solutions and a significant effect of the wellbore variable injection rate on the evolution of the fracture length.


Sign in / Sign up

Export Citation Format

Share Document