Laboratory demonstration of hydraulic fracture height growth across weak discontinuities

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. MR93-MR105 ◽  
Author(s):  
Pengju Xing ◽  
Keita Yoshioka ◽  
Jose Adachi ◽  
Amr El-Fayoumi ◽  
Andrew P. Bunger

Decades of research have led to numerous insights in modeling the impact of stresses and rock properties on hydraulic fracture height growth. However, the conditions under which weak horizontal interfaces are expected to impede height growth remain for the most part unknown. We have developed an experimental study of the impact of weak horizontal discontinuities on hydraulic fracture height growth, including the influences of (1) abrupt stress contrasts between layers, (2) material fracture toughness, and (3) contrasts of stiffness between the reservoir and bounding layers. The experiments are carried out with an analog three-layered medium constructed from transparent polyurethane, considering toughnesses resisting vertical fracture growth. There are four observed geometries: containment, height growth, T-shape growth, and the combination of height growth and T-shape. Results are developed in a parametric space embodying the influence of the horizontal stress contrast, vertical stress, and horizontal barrier stress contrast, as well as the fluid pressure. The results indicate that these cases fall within distinct regions when plotted in the parametric space. The locations in the parametric space of these regions are strongly impacted by the vertical fracture toughness: Increasing the value of the vertical interface fracture toughness leads to a suppression of height growth in favor of containment and T-shaped growth. Besides providing detailed experimental data for benchmarking 3D hydraulic fracture simulators, these experiments show that the fracture height is substantially less than would be predicted in the absence of the weak horizontal discontinuities.

1983 ◽  
Vol 23 (06) ◽  
pp. 870-878 ◽  
Author(s):  
Ian D. Palmer ◽  
H.B. Carroll

Abstract Models of three-dimensional (3D) fracture propagation are being developed to study the effect of variations of stress and rock properties on fracture height and bottomhole pressure (BHP). Initially a blanket sand bounded by zones of higher minimum in-situ stress is considered, with stresses symmetrical about both the pay-zone axis and the wellbore. An elliptical fracture perimeter is assumed. Fluid flows are one-dimensional (1D) Newtonian in the direction of the pay zone. Two models, FL1 and FL2, are developed. In FL1, a discontinuous stress variation is approximated by a y2 variation in the vertical coordinate, and the fracture criterion, Ki = Kc, is satisfied at both major and minor axes. The net pressure at the tip, Lf, of the long axis required by the boundary condition Ki = Kc does not seem crucial in determining fracture height or BHP (compare with one group of published models that assumes p = 0 at Lf). Model FL2 properly represents the discontinuous stresses, and satisfies Ki = Kc at the wellbore but not at the tip of the long axis. A parametric study is made, with both models, of the comparative effects of stress contrast, Kc, pay-zone height, h, and Young's modulus, E, on fracture height and BHP. Results indicate that Kc does not have as much effect as either E or, at least for large stress contrasts. Model FL2 suggests the possibility of a rapid growth in fracture height as is reduced. Such modeling may be able to give an upper or "safe" limit on the pumping parameters ( and ) to ensure good containment. When the stress contrast is high, 700 psi [4826 kPa], an analytic derivation of BHP appears to be a good approximation for the parameters we use, if everywhere the fracture height is assumed equal to the pay zone height. Although leakoff is neglected here, subsequent modeling results show that, for leak off coefficients 0.001 ft- min [3.9 × 10 -5 m.s ], the results herein are a good approximation to the case when leak off is included. Introduction In their essence, models of hydraulic fracture propagation involve elasticity theory and fluid mechanics. The first is concerned with the fracture opening or width, w(p), as a function of net pressure on the fracture faces, while the second is concerned with the pressure drop, p(w), caused by the flow of viscous fluids in the fracture. Simultaneous solution of these equations includes a boundary condition that often takes the form Ki = Kc, where Ki is the stress-intensity factor at a point on the fracture tip, and Kc is the fracture toughness. The final solution is very complex in 3D, when a vertical fracture can expand vertically as well as horizontally along the pay zone. Thus, the first solutions were essentially two-dimensional (2D), and they assumed that the fracture height, hf, was fixed at the pay zone height, h. The 2D solutions were clustered in two groups as summarized by Nordgren, Perkins, and Geertsma and Haafkens. The first grouping, based on a model by Christianovich and Zheltov, assumed that the sides of an elongated, vertical fracture were parallel (i.e., free slippage between the pay and bounding zones, or no vertical stiffness). Other papers in this grouping included Geertsma and de Klerk, Daneshy and Settari. SPEJ P. 870^


SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 302-318 ◽  
Author(s):  
Jixiang Huang ◽  
Joseph P. Morris ◽  
Pengcheng Fu ◽  
Randolph R. Settgast ◽  
Christopher S. Sherman ◽  
...  

Summary A fully coupled finite-element/finite-volume code is used to model 3D hydraulically driven fractures under the influence of strong vertical variations in closure stress interacting with natural fractures. Previously unknown 3D interaction mechanisms on fracture-height growth are revealed. Slipping of a natural fracture, triggered by elevated fluid pressure from an intersecting hydraulic fracture, can induce both increases and decreases of normal stress in the minimum-horizontal-stress direction, toward the center and tip of the natural fracture, respectively. Consequently, natural fractures are expected to be able to both encourage and inhibit the progress of hydraulic fractures propagating through stress barriers, depending on the relative locations between the intersecting fractures. Once the hydraulic fracture propagates above the stress barrier through the weakened segment near a favorably located natural fracture, a configuration consisting of two opposing fractures cuts the stress barrier from above and below. The fluid pressure required to break the stress barrier under such opposing-fracture configurations is substantially lower than that required by a fracture penetrating the same barrier from one side. Sensitivity studies of geologic conditions and operational parameters have also been performed to explore the feasibility of controlled fracture height. The interactions between hydraulic fractures, natural fractures, and geologic factors such as stress barriers in three dimensions are shown to be much more complex than in two dimensions. Although it is impossible to exhaust all the possible configurations, the ability of a 3D, fully coupled numerical model to naturally capture these processes is well-demonstrated.


2015 ◽  
Author(s):  
Songxia Liu ◽  
Peter P. Valkó

Abstract Fracture height is a critical input parameter for 2D hydraulic fracturing design models, and also an important output result of 3D models. While many factors may influence fracture height evolution in multi-layer formations, the consensus is that the so called “equilibrium height belonging to a certain treating pressure” provides an upper limit, at least for non-naturally fractured media. How to solve the “equilibrium height” problem has been known since the 1970s. However, because of the complexity of the algebra involved, published equations are overly simplified and do not provide reliable results. We revisited the equilibrium height problem, and our theoretical and numerical investigations led to a new model that fully characterizes height evolution amid various formation properties (fracture toughness, in-situ stress, thickness, etc.). The new model, for the first time, rigorously solves the equilibrium height mathematically. With the help of computer algebra software, we applied a definition of fracture toughness, incorporated the effects of hydrostatic pressure, and considered non-symmetric variations of layered formation properties. Results were determined for the classic 3-layer problem and then were extended to 6 layers. The effects of reservoir and fluid properties on the fracture height growth were investigated. Tip jump is caused by low in-situ stress; tip stability is imposed by large fracture toughness and/or large in-situ stress. If the fluid density is ignored, the result regarding to which tip will grow into infinity could be totally different. For any multi-layer formation problem, two, 3-layer pseudo problems were constructed to create an outer and inner height envelope, in order to assess the potential effects of reservoir parameter uncertainties on height profile. Second solution pair in the 3-layer problem was investigated numerically and analytically, to avoid misleading results in the 3D models. This improved model can rapidly (in seconds) and reliably calculate the theoretical maximum equilibrium fracture height in layered formations with various rock mechanical properties. Then, the equilibrium height can be used to (1) provide input data for 2D model, (2) improve 3D model governing equations, (3) determine the net pressure needed to achieve a certain height growth, and (4) obtain the maximum net pressure assuring no fracture invasion into aquifers. This model may be incorporated into current hydraulic fracture simulation software to yield more accurate and cost-effective hydraulic fracturing designs.


2012 ◽  
Vol 27 (01) ◽  
pp. 8-19 ◽  
Author(s):  
M. Kevin Fisher ◽  
Norman R. Warpinski

2021 ◽  
pp. petgeo2020-095
Author(s):  
Michael J. Steventon ◽  
Christopher A-L. Jackson ◽  
Howard D. Johnson ◽  
David M. Hodgson ◽  
Sean Kelly ◽  
...  

The geometry, distribution, and rock properties (i.e. porosity and permeability) of turbidite reservoirs, and the processes associated with turbidity current deposition, are relatively well known. However, less attention has been given to the equivalent properties resulting from laminar sediment gravity-flow deposition, with most research limited to cogenetic turbidite-debrites (i.e. transitional flow deposits) or subsurface studies that focus predominantly on seismic-scale mass-transport deposits (MTDs). Thus, we have a limited understanding of the ability of sub-seismic MTDs to act as hydraulic seals and their effect on hydrocarbon production, and/or carbon storage. We investigate the gap between seismically resolvable and sub-seismic MTDs, and transitional flow deposits on long-term reservoir performance in this analysis of a small (<10 km radius submarine fan system), Late Jurassic, sandstone-rich stacked turbidite reservoir (Magnus Field, northern North Sea). We use core, petrophysical logs, pore fluid pressure, quantitative evaluation of minerals by scanning electron microscopy (QEMSCAN), and 3D seismic-reflection datasets to quantify the type and distribution of sedimentary facies and rock properties. Our analysis is supported by a relatively long (c. 37 years) and well-documented production history. We recognise a range of sediment gravity deposits: (i) thick-/thin- bedded, structureless and structured turbidite sandstone, constituting the primary productive reservoir facies (c. porosity = 22%, permeability = 500 mD), (ii) a range of transitional flow deposits, and (iii) heterogeneous mud-rich sandstones interpreted as debrites (c. porosity = <10%, volume of clay = 35%, up to 18 m thick). Results from this study show that over the production timescale of the Magnus Field, debrites act as barriers, compartmentalising the reservoir into two parts (upper and lower reservoir), and transitional flow deposits act as baffles, impacting sweep efficiency during production. Prediction of the rock properties of laminar and transitional flow deposits, and their effect on reservoir distribution, has important implications for: (i) exploration play concepts, particularly in predicting the seal potential of MTDs, (ii) pore pressure prediction within turbidite reservoirs, and (iii) the impact of transitional flow deposits on reservoir quality and sweep efficiency.Supplementary material:https://doi.org/10.6084/m9.figshare.c.5313860


1986 ◽  
Vol 108 (2) ◽  
pp. 107-115 ◽  
Author(s):  
I. D. Palmer ◽  
C. T. Luiskutty

There is a pressing need to compare and evaluate hydraulic fracture models which are now being used by industry to predict variable fracture height. The fractures of concern here are vertical fractures which have a pronounced elongation in the direction of the payzone, i.e., there is a dominant one-dimensional fluid flow along the payzone direction. A summary is given of the modeling entailed in the basic ORU fracture model, which calculates fracture height as a function of distance from the wellbore in the case of a continuous sand bounded by zones of higher (but equal) minimum in-situ stress. The elastic parameters are assumed the same in each layer, and injected flow rates and fluid parameters are taken to be constant. Leak-off is included with spurt loss, as well as non-Newtonian flow. An advantage of the model is its small computer run time. Predictions for wellbore height and pressure from the ORU model are compared separately with the AMOCO and MIT pseudo-3D models. In one instance of high stress contrast the ORU wellbore pressure agrees fairly well with the AMOCO model, but the AMOCO wellbore height is greater by 32 percent. Comparison between the ORU and MIT models in two cases (also high stress contrast) indicates height disagreement at the wellbore by factors of 1.5–2.5 with the MIT model giving a lower height. Thus it appears there can be substantial discrepancies between all three models. Next we compare the ORU model results with six cases of elongated fractures from the TERRA-TEK fully-3D model. Although two of these cases are precluded due to anomolous discrepancies, the other four cases show reasonable agreement. We make a critical examination of assumptions that differ in all the models (e.g., the effective modulus-stiffness multiplier approximation in the AMOCO model, the effect of finite fluid flow in the vertical direction in the MIT model, and the effect of 2D flow and limited perforated height in the TERRA-TEK model). Suggestions are made for reconciling some of the discrepancies between the various models. For example, the ORU/AMOCO height discrepancy appears to be resolved; for other discrepancies we have no explanation. Our main conclusion is that the AMOCO, TERRA-TEK and ORU models for fracture height and bottomhole pressure are in reasonable agreement for highly elongated fractures. Despite the difficulties in understanding the different models, the comparisons herein are an encouraging first step towards normalizing these hydraulic fracture models.


SPE Journal ◽  
2019 ◽  
Vol 24 (05) ◽  
pp. 2148-2162 ◽  
Author(s):  
Pengcheng Fu ◽  
Jixiang Huang ◽  
Randolph R. Settgast ◽  
Joseph P. Morris ◽  
Frederick J. Ryerson

Summary The height growth of a hydraulic fracture is known to be affected by many factors that are related to the layered structure of sedimentary rocks. Although these factors are often used to qualitatively explain why hydraulic fractures usually have well–bounded height growth, most of them cannot be directly and quantitatively characterized for a given reservoir to enable a priori prediction of fracture–height growth. In this work, we study the role of the “roughness” of in–situ–stress profiles, in particular alternating low and high stress among rock layers, in determining the tendency of a hydraulic fracture to propagate horizontally vs. vertically. We found that a hydraulic fracture propagates horizontally in low–stress layers ahead of neighboring high–stress layers. Under such a configuration, a fracture–mechanics principle dictates that the net pressure required for horizontal growth of high–stress layers within the current fracture height is significantly lower than that required for additional vertical growth across rock layers. Without explicit consideration of the stress–roughness profile, the system behaves as if the rock is tougher against vertical propagation than it is against horizontal fracture propagation. We developed a simple relationship between the apparent differential rock toughness and characteristics of the stress roughness that induce equivalent overall fracture shapes. This relationship enables existing hydraulic–fracture models to represent the effects of rough in–situ stress on fracture growth without directly representing the fine–resolution rough–stress profiles.


2011 ◽  
Vol 396-398 ◽  
pp. 2420-2423 ◽  
Author(s):  
Yong Quan Hu ◽  
Jin Zhou Zhao ◽  
Tao Lin

It is necessary to control hydraulic fracture height growth by creating artificial barrier for reservoir with thin payzone or weak overburdens and underburdens. In this paper, the influential factors on resistance role of artificial barrier were firstly identified and the parameters value range was determined. Then, an affixing intermediate container was combined into the core flow test device in order to ensure buoyant diverter well dispersed into carrying fluid in experiment test, and the particle size distribution of buoyant additives was tested by Malvern laser counter. The crack was made about 2/3 length in artificial core for simulating artificial barrier process by core flow test which were accomplished in base of recommendation practice of petroleum industry. Lastly, successive regression method was applied and a statistical relationship of the barrier strength with carrying fluid viscosity, floating additives concentration and core permeability. Then scouring experiments are made under modeling hydraulically fracturing treatment conditions. It was proven that the artificial barrier was steady at current treatment parameters and that could be effectively controlled hydraulic fracture height growth.


Sign in / Sign up

Export Citation Format

Share Document