scholarly journals Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

SPE Journal ◽  
2021 ◽  
pp. 1-15
Author(s):  
Wei Fu ◽  
Joseph P. Morris ◽  
Pengcheng Fu ◽  
Jixiang Huang ◽  
Christopher S. Sherman ◽  
...  

Summary This work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

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

Directional rupture is one of the most important and most common problems related to rock breaking. The goal of directional rock breaking can be effectively achieved via multi-hole linear co-directional hydraulic fracturing. In this paper, the XSite software was utilized to verify the experimental results of multi-hole linear co-directional hydraulic fracturing., and its basic law is studied. The results indicate that the process of multi-hole linear co-directional hydraulic fracturing can be divided into four stages: water injection boost, hydraulic fracture initiation, and the unstable and stable propagation of hydraulic fracture. The stable expansion stage lasts longer and produces more microcracks than the unstable expansion stage. Due to the existence of the borehole-sealing device, the three-dimensional hydraulic fracture first initiates and expands along the axial direction in the bare borehole section, then extends along the axial direction in the non-bare hole section and finally expands along the axial direction in the rock mass without the borehole. The network formed by hydraulic fracture in rock is not a pure plane, but rather a curved spatial surface. The curved spatial surface passes through both the centre of the borehole and the axial direction relative to the borehole. Due to the boundary effect, the curved spatial surface goes toward the plane in which the maximum principal stress occurs. The local ground stress field is changed due to the initiation and propagation of hydraulic fractures. The propagation direction of the fractures between the fracturing boreholes will be deflected. A fracture propagation pressure that is greater than the minimum principle stress and a tension field that is induced in the leading edge of the fracture end, will aid to fracture intersection; as a result, the possibility of connecting the boreholes will increase.


2015 ◽  
Author(s):  
Mark W. McClure ◽  
Mohsen Babazadeh ◽  
Sogo Shiozawa ◽  
Jian Huang

Abstract We developed a hydraulic fracturing simulator that implicitly couples fluid flow with the stresses induced by fracture deformation in large, complex, three-dimensional discrete fracture networks. The simulator can describe propagation of hydraulic fractures and opening and shear stimulation of natural fractures. Fracture elements can open or slide, depending on their stress state, fluid pressure, and mechanical properties. Fracture sliding occurs in the direction of maximum resolved shear stress. Nonlinear empirical relations are used to relate normal stress, fracture opening, and fracture sliding to fracture aperture and transmissivity. Fluid leakoff is treated with a semianalytical one-dimensional leakoff model that accounts for changing pressure in the fracture over time. Fracture propagation is treated with linear elastic fracture mechanics. Non-Darcy pressure drop in the fractures due to high flow rate is simulated using Forchheimer's equation. A crossing criterion is implemented that predicts whether propagating hydraulic fractures will cross natural fractures or terminate against them, depending on orientation and stress anisotropy. Height containment of propagating hydraulic fractures between bedding layers can be modeled with a vertically heterogeneous stress field or by explicitly imposing hydraulic fracture height containment as a model assumption. The code is efficient enough to perform field-scale simulations of hydraulic fracturing with a discrete fracture network containing thousands of fractures, using only a single compute node. Limitations of the model are that all fractures must be vertical, the mechanical calculations assume a linearly elastic and homogeneous medium, proppant transport is not included, and the locations of potentially forming hydraulic fractures must be specified in advance. Simulations were performed of a single propagating hydraulic fracture with and without leakoff to validate the code against classical analytical solutions. Field-scale simulations were performed of hydraulic fracturing in a densely naturally fractured formation. The simulations demonstrate how interaction with natural fractures in the formation can help explain the high net pressures, relatively short fracture lengths, and broad regions of microseismicity that are often observed in the field during stimulation in low permeability formations, and which are not predicted by classical hydraulic fracturing models. Depending on input parameters, our simulations predicted a variety of stimulation behaviors, from long hydraulic fractures with minimal leakoff into surrounding fractures to broad regions of dense fracturing with a branching network of many natural and newly formed fractures.


2021 ◽  
Author(s):  
Debotyam Maity ◽  
Jordan Ciezobka

Abstract In this case study, we apply a novel fracture imaging and interpretation workflow to take a systematic look at hydraulic fractures captured during thorugh fracture coring at the Hydraulic Fracturing Test Site (HFTS) in Midland Basin. Digital fracture maps rendered using high resolution 3D laser scans are analyzed for fracture morphology and roughness. Analysis of hydraulic fracture faces show that the roughness varies systematically in clusters with average cluster separation of approximately 20' along the core. While isolated smooth hydraulic fractures are observed in the dataset, very rough fractures are found to be accompanied by proximal smoother fractures. Roughness distribution also helps understand the effect of stresses on fracture distribution. Locally, fracture roughness seems to vary with fracture orientations indicating possible inter-fracture stress effects. At the scale of stage lengths however, we see evidence of inter-stage stress effects. We also observe fracture morphology being strongly driven by rock properties and changes in lithology. Identified proppant distribution along the cored interval is also correlated with roughness variations and we observe strong positive correlation between proppant concentrations and fracture roughness at the local scale. Finally, based on the observed distribution of hydraulic fracture properties, we propose a conceptual spatio-temporal model of fracture propagation which can help explain the hydraulic fracture roughness distribution and ties in other observations as well.


2015 ◽  
Author(s):  
Hisanao Ouchi ◽  
Amit Katiyar ◽  
John T. Foster ◽  
Mukul M. Sharma

Abstract A novel fully coupled hydraulic fracturing model based on a nonlocal continuum theory of peridynamics is presented and applied to the fracture propagation problem. It is shown that this modeling approach provides an alternative to finite element and finite volume methods for solving poroelastic and fracture propagation problems and offers some clear advantages. In this paper we specifically investigate the interaction between a hydraulic fracture and natural fractures. Current hydraulic fracturing models remain limited in their ability to simulate the formation of non-planar, complex fracture networks. The peridynamics model presented here overcomes most of the limitations of existing models and provides a novel approach to simulate and understand the interaction between hydraulic fractures and natural fractures. The model predictions in two-dimensions have been validated by reproducing published experimental results where the interaction between a hydraulic fracture and a natural fracture is controlled by the principal stress contrast and the approach angle. A detailed parametric study involving poroelasticity and mechanical properties of the rock is performed to understand why a hydraulic fracture gets arrested or crosses a natural fracture. This analysis reveals that the poroelasticity, resulting from high fracture fluid leak-off, has a dominant influence on the interaction between a hydraulic fracture and a natural fracture. In addition, the fracture toughness of the rock, the toughness of the natural fracture, and the shear strength of the natural fracture also affect the interaction between a hydraulic fracture and a natural fracture. Finally, we investigate the interaction of multiple completing fractures with natural fractures in two-dimensions and demonstrate the applicability of the approach to simulate complex fracture networks on a field scale.


SPE Journal ◽  
2016 ◽  
Vol 21 (04) ◽  
pp. 1302-1320 ◽  
Author(s):  
Mark W. McClure ◽  
Mohsen Babazadeh ◽  
Sogo Shiozawa ◽  
Jian Huang

Summary We developed a hydraulic-fracturing simulator that implicitly couples fluid flow with the stresses induced by fracture deformation in large, complex, 3D discrete-fracture networks (DFNs). The code is efficient enough to perform field-scale simulations of hydraulic fracturing in DFNs containing thousands of fractures, without relying on distributed-memory parallelization. The simulator can describe propagation of hydraulic fractures and opening and shear stimulation of natural fractures. Fracture elements can open or slide, depending on their stress state, fluid pressure, and mechanical properties. Fracture sliding occurs in the direction of maximum resolved shear stress. Nonlinear empirical equations are used to relate normal stress, fracture opening, and fracture sliding to fracture aperture and transmissivity. Fluid leakoff is treated with a semianalytical 1D leakoff model that accounts for changing pressure in the fracture over time. Fracture propagation is modeled with linear-elastic fracture mechanics. The Forchheimer equation (Forchheimer 1901) is used to simulate non-Darcy pressure drop in the fractures because of high flow rate. A crossing criterion is implemented that predicts whether propagating hydraulic fractures will cross natural fractures or terminate against them, depending on orientation and stress anisotropy. Height containment of propagating hydraulic fractures between bedding layers can be modeled with a vertically heterogeneous stress field or by explicitly imposing hydraulic-fracture-height containment as a model assumption. Limitations of the model are that all fractures must be vertical; the mechanical calculations assume a linearly elastic and homogeneous medium; proppant transport is not included; and the locations of potentially forming hydraulic fractures must be specified in advance. Simulations were performed of a single propagating hydraulic fracture with and without leakoff to validate the code against classical analytical solutions. Field-scale simulations were performed of hydraulic fracturing in a densely naturally fractured formation. The simulations demonstrate how interaction with natural fractures in the formation can help explain the high net pressures, relatively short fracture lengths, and broad regions of microseismicity that are often observed in the field during stimulation in low-permeability formations, and that are not predicted by classical hydraulic-fracturing models. Depending on input parameters, our simulations predicted a variety of stimulation behaviors, from long hydraulic fractures with minimal leakoff into surrounding fractures to broad regions of dense fracturing with a branching network of many natural and newly formed fractures.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-24
Author(s):  
Xin Zhang ◽  
Yuqi Zhang ◽  
Bingxiang Huang

Hydraulic fracturing applications have shown a stress disturbance effect during hydraulic fracture propagation, which is often ignored. Using laboratory and discrete element numerical simulation tests, hydraulic fracture propagation under this stress disturbance is systematically studied. The results show that during hydraulic fracturing, the bedding plane is damaged by the stress disturbance, forming a bedding fracture zone (BFZ). The nonlinear fracture characteristics of the formation process of the disturbed fracture zone are revealed, and two indexes (the number of fractures in the disturbed fracture zone and the size of the disturbed fracture zone) are proposed to evaluate the fracturing effect of the stress disturbance. Based on these indexes, multifactor sensitivity tests are conducted under different geological conditions and operational factors. When the principal stress ( σ 1 ) difference is large, the number of hydraulic fractures gradually decreases from many to one, and the direction of the hydraulic fractures gradually approaches the vertical direction of σ 3 , but the change in the in situ stress condition has no obvious effect on the stress disturbance effect. The weaker the bonding strength of the bedding plane, the more significant the stress disturbance effect is, and the easier it is for the fractures to expand along the bedding plane. With increasing injection rate, the stress disturbance effect first increases and then decreases, and the hydraulic fracture propagates from along the bedding plane to cross the bedding plane. With increasing relative distance between the injection hole and bedding plane, the stress disturbance effect presents a linearly increasing trend, and the hydraulic fractures along the bedding planes extend. Based on the experimental results, the relationship between the fracturing effect of the stress disturbance and the extension mode of the hydraulic fracture is determined, and an optimization method for hydraulic fracturing in composite rock reservoirs is given. The research results can provide a theoretical basis for controlling the formation of complex fracture networks in composite rock reservoirs.


SPE Journal ◽  
2021 ◽  
pp. 1-12
Author(s):  
Yunhui Tan ◽  
Shugang Wang ◽  
Margaretha C. M. Rijken ◽  
Kelly Hughes ◽  
Ivan Lim Chen Ning ◽  
...  

Summary Recently more distributed acoustic sensing (DAS) data have been collected during hydraulic fracturing in shale. Low-frequency DAS signals show patterns that are intuitively consistent with the understanding of the strain field around hydraulic fractures. This study uses a fracture simulator combined with a finite element solver to further understand the various patterns of the strain field caused by hydraulic fracturing. The results can serve as a “type-curve” template for the further interpretation of cross-well strain field plots. Incorporating detailed pump schedule and fracturing fluid/proppant properties, we use a hydraulic fracture simulator to generate fracture geometries, which are then passed to a finite element solver as boundary conditions for elastic-static calculation of the strain field. Because the finite element calculated strain is a tensor, it needs to be projected along the monitoring well trajectory to be comparable with the DAS strain, which is uniaxial. Moreover, the calculated strain field is transformed into a time domain using constant fracture propagation velocity. Strain rate is further derived from the simulated strain field using differentiation along the fracture propagation direction. Scenarios including a single planar hydraulic fracture, a single fracture with a discrete fracture network (DFN), and multiple planar hydraulic fractures in both vertical and horizontal directions were studied. The scenarios can be differentiated in the strain patterns on the basis of the finite element simulation results. In general, there is a tensile heart-shaped zone in front of the propagating fracture tip shown along the horizontal strain direction on both strain and strain rate plots. On the sides, there are compressional zones parallel to the fracture. The strain field projects beyond the depth where the hydraulic fracture is present. Patterns from strain rate can be used to distinguish whether the fracture is intersecting the fiber. Along the vertical direction, the transition zone depicts the upper boundary of the fracture. A complex fracture network with DFN shows a much more complex pattern compared with a single planar fracture. Multiple planar fractures show polarity reversals in horizontal fiber because of interactions between fractures. Data from the Hydraulic Fracturing Test Site 2 (HFTS2) experiment were used to validate the simulated results. The application of the study is to provide a template to better interpret hydraulic fracture characteristics using low-frequency DAS strain-monitoring data. To our understanding, there are no comprehensive templates for engineers to understand the strain signals from cross-well fiber monitoring. The results of this study will guide engineers toward better optimization of well spacing and fracturing design to minimize well interference and improve efficiency.


2022 ◽  
Vol 2152 (1) ◽  
pp. 012048
Author(s):  
Zhongshan Shen ◽  
Hui Xue ◽  
Zhiqiang Bai

Abstract Perforation azimuth has an important influence on the nucleation, propagation path and morphology of hydraulic fractures. In this paper, the true triaxial hydraulic fracturing simulation experimental system is used to analyze the hydraulic fracture morphology and propagation path under different perforation azimuth angles. With the increase of the azimuth angle of perforation, the stable fracture propagation pressure of the fracturing sample also increases. When the azimuth angle of perforation is 0°, the propagation pressure is about 18 MPa, and when the azimuth angle of perforation is 90°, the propagation pressure is about 26.5 MPa, increasing by nearly 47.22%.


2015 ◽  
Vol 137 (6) ◽  
Author(s):  
Desheng Zhou ◽  
Peng Zheng ◽  
Jiao Peng ◽  
Pei He

Creating complex fracture network by hydraulic fracturing operation in unconventional reservoir development is the key factor of effective exploitation. The mechanism of creating a fracture network is not clear up to today. Conventional hydraulic fracturing theory is based on tensile failure of a rock, and a hydraulic fracture is widely accepted as propagating along the direction of in situ maximum horizontal principal stress in the industry. Based on rock elastic mechanics and fracture mechanics, considering combined tensile and shear failures, the maximum circumferential strain criterion and boundary element method (BEM), the paper studies the induced stress and its variation during a fracture propagation, the interaction between two or more hydraulic fractures, and the interaction between a hydraulic fracture and a natural crack. The paper shows that a propagating fracture will produce induced stresses on surrounding rock and form a stress shadow. Instead of propagation along the direction of maximum horizontal principal stress as a single fracture, the outside two fractures of two or more hydraulic fractures are exclusive and turning away from each other. A natural crack may be awaked and extend at its both tips by a propagating hydraulic fracture before their intersection, and the hydraulic fracture may deflect toward the natural crack. The interaction between a hydraulic fracture and a natural crack depends on the transverse distance between them and the initial length of the crack. The shorter the transverse distance and the longer the crack length are, the higher the possibility of the crack to be awaked is. The research results are helpful in understanding complex fracture network and may be used in determining hydraulic fracture places to create a complex fracture network.


2021 ◽  
Author(s):  
Ghazal Izadi ◽  
Colleen Barton ◽  
Pierre-Francois Roux ◽  
Tebis Llobet ◽  
Thiago Pessoa ◽  
...  

Abstract For tight reservoirs where hydraulic fracturing is required to enable sufficient fluid mobility for economic production, it is critical to understand the placement of induced fractures, their connectivity, extent, and interaction with natural fractures within the system. Hydraulic fracture initiation and propagation mechanisms are greatly influenced by the effect of the stress state, rock fabric and pre-existing features (e.g. natural fractures, faults, weak bedding/laminations). A pre-existing natural fracture system can dictate the mode, orientation and size of the hydraulic fracture network. A better understanding of the fracture growth phenomena will enhance productivity and also reduce the environmental footprint as less fractures can be created in a much more efficient way. Assessing the role of natural fractures and their interaction with hydraulic fractures in order to account for them in the hydraulic fracture model is achieved by leveraging microseismicity. In this study, we have used a combination of borehole and surface microseismic monitoring to get high vertical resolution locations and source mechanisms. 3D numerical modelling of hydraulic fracturing in complex geological conditions to predict fracture propagation is essential. 3D hydraulic fracturing simulation includes modelling capabilities of stimulation parameters, true 3D fracture propagation with near wellbore 3D complexity including a coupled DFN and the associated microseismic event generation capability. A 3D hydraulic fracture model was developed and validated by matching model predictions to microseismic observations. Microseismic source mechanisms are leveraged to determine the location and geometry of pre-existing features. In this study, we simulate a DFN based on the recorded seismicity of multi stage hydraulic fractures in a horizontal well. The advanced 3D hydraulic fracture modelling software can integrate effectively and efficiently data from a variety of multi-disciplinary sources and scales to create a subsurface characterization of the unconventional reservoir. By incorporating data from 3D seismic, LWD/wireline, core, completion/stimulation monitoring, and production, the software generates a holistic reservoir model embedded in a modular, multi-physics software platform of coupled numerical solvers that capture the fundamental physics of the processes being modelled. This study illustrates the importance of a powerful software tool that captures the necessary physics of stimulation to predict the effects of various completion designs and thereby ensure the most accurate representation of an unconventional reservoir response to a stimulation treatment.


Sign in / Sign up

Export Citation Format

Share Document