The Tip Region of a Near-Surface Hydraulic Fracture

2018 ◽  
Vol 85 (4) ◽  
Author(s):  
Zhi-Qiao Wang ◽  
Emmanuel Detournay

This paper investigates the tip region of a hydraulic fracture propagating near a free-surface via the related problem of the steady fluid-driven peeling of a thin elastic layer from a rigid substrate. The solution of this problem requires accounting for the existence of a fluid lag, as the pressure singularity that would otherwise exist at the crack tip is incompatible with the underlying linear beam theory governing the deflection of the thin layer. These considerations lead to the formulation of a nonlinear traveling wave problem with a free boundary, which is solved numerically. The scaled solution depends only on one number K, which has the meaning of a dimensionless toughness. The asymptotic viscosity- and toughness-dominated regimes, respectively, corresponding to small and large K, represent the end members of a family of solutions. It is shown that the far-field curvature can be interpreted as an apparent toughness, which is a universal function of K. In the viscosity regime, the apparent toughness does not depend on K, while in the toughness regime, it is equal to K. By noting that the apparent toughness represents an intermediate asymptote for the layer curvature under certain conditions, the obtention of time-dependent solutions for propagating near-surface hydraulic fractures can be greatly simplified. Indeed, any such solutions can be constructed by a matched asymptotics approach, with the outer solution corresponding to a uniformly pressurized fracture and the inner solution to the tip solution derived in this paper.

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.


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

To better evaluate the spatial steering effect of directional perforation hydraulic fractures, evaluation indexes for the spatial steering effect are first proposed in this paper. Then, these indexes are used to quantitatively evaluate existing physical experimental results. Finally, with the help of RFPA2D-Flow software, the influence of perforation length and azimuth on the spatial steering process of hydraulic fracture are quantitatively analysed using four evaluation indexes. It is shown by the results that the spatial deflection trajectory, deflection distance, deflection angle and initiation pressure of hydraulic fractures can be used as quantitative evaluation indexes for the spatial steering effect of hydraulic fractures. The deflection paths of directional perforation hydraulic fractures are basically the same. They all gradually deflect to the maximum horizontal principal stress direction from the perforation hole and finally represent a double-wing bending fracture. The deflection distance, deflection angle and initiation pressure of hydraulic fractures increase gradually with increasing perforation azimuth, and the sensitivity of the deflection angle to the perforation azimuth of hydraulic fractures also increases. With increasing perforation length, the deflection distance of hydraulic fractures increases gradually. However, the deflection angle and initiation pressure decrease gradually, as does the sensitivity.


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.


Energies ◽  
2021 ◽  
Vol 14 (9) ◽  
pp. 2610
Author(s):  
Wenzheng Liu ◽  
Qingdong Zeng ◽  
Jun Yao ◽  
Ziyou Liu ◽  
Tianliang Li ◽  
...  

Rock yielding may well take place during hydraulic fracturing in deep reservoirs. The prevailing models based on the linear elastic fracture mechanics (LEFM) are incapable of describing the evolution process of hydraulic fractures accurately. In this paper, a hydro-elasto-plastic model is proposed to investigate the hydraulic fracture propagation in deep reservoirs. The Drucker–Prager plasticity model, Darcy’s law, cubic law and cohesive zone model are employed to describe the plastic deformation, matrix flow, fracture flow and evolution of hydraulic fractures, respectively. Combining the embedded discrete fracture model (EDFM), extended finite element method (XFEM) and finite volume method, a hybrid numerical scheme is presented to carry out simulations. A dual-layer iterative procedure is developed based on the fixed-stress split method, Picard iterative method and Newton–Raphson iterative method. The iterative procedure is used to deal with the coupling between nonlinear deformation with fracture extension and fluid flow. The proposed model is verified against analytical solutions and other numerical simulation results. A series of numerical cases are performed to investigate the influences of rock plasticity, internal friction angle, dilatancy angle and permeability on hydraulic fracture propagation. Finally, the proposed model is extended to simulate multiple hydraulic fracture propagation. The result shows that plastic deformation can enhance the stress-shadowing effect.


2014 ◽  
Vol 10 (1) ◽  
pp. 91-106 ◽  
Author(s):  
E. Dietze ◽  
F. Maussion ◽  
M. Ahlborn ◽  
B. Diekmann ◽  
K. Hartmann ◽  
...  

Abstract. Grain-size distributions offer powerful proxies of past environmental conditions that are related to sediment sorting processes. However, they are often of multimodal character because sediments can get mixed during deposition. To facilitate the use of grain size as palaeoenvironmental proxy, this study aims to distinguish the main detrital processes that contribute to lacustrine sedimentation across the Tibetan Plateau using grain-size end-member modelling analysis. Between three and five robust grain-size end-member subpopulations were distinguished at different sites from similarly–likely end-member model runs. Their main modes were grouped and linked to common sediment transport and depositional processes that can be associated with contemporary Tibetan climate (precipitation patterns and lake ice phenology, gridded wind and shear stress data from the High Asia Reanalysis) and local catchment configurations. The coarse sands and clays with grain-size modes >250 μm and <2 μm were probably transported by fluvial processes. Aeolian sands (~200 μm) and coarse local dust (~60 μm), transported by saltation and in near-surface suspension clouds, are probably related to occasional westerly storms in winter and spring. Coarse regional dust with modes ~25 μm may derive from near-by sources that keep in longer term suspension. The continuous background dust is differentiated into two robust end members (modes: 5–10 and 2–5 μm) that may represent different sources, wind directions and/or sediment trapping dynamics from long-range, upper-level westerly and episodic northerly wind transport. According to this study grain-size end members of only fluvial origin contribute small amounts to mean Tibetan lake sedimentation (19± 5%), whereas local to regional aeolian transport and background dust deposition dominate the clastic sedimentation in Tibetan lakes (contributions: 42 ± 14% and 51 ± 11%). However, fluvial and alluvial reworking of aeolian material from nearby slopes during summer seems to limit end-member interpretation and should be crosschecked with other proxy information. If not considered as a stand-alone proxy, a high transferability to other regions and sediment archives allows helpful reconstructions of past sedimentation history.


2021 ◽  
Author(s):  
Yinghui Wu ◽  
Robert Hull ◽  
Andrew Tucker ◽  
Craig Rice ◽  
Peter Richter ◽  
...  

Abstract Distributed fiber-optic sensing (DFOS) has been utilized in unconventional reservoirs for hydraulic fracture efficiency diagnostics for many years. Downhole fiber cables can be permanently installed external to the casing to monitor and measure the uniformity and efficiency of individual clusters and stages during the completion in the near-field wellbore environment. Ideally, a second fiber or multiple fibers can be deployed in offset well(s) to monitor and characterize fracture geometries recorded by fracture-driven interactions or frac-hits in the far-field. Fracture opening and closing, stress shadow creation and relaxation, along with stage isolation can be clearly identified. Most importantly, fracture propagation from the near to far-field can be better understood and correlated. With our current technology, we can deploy cost effective retrievable fibers to record these far-field data. Our objective here is to highlight key data that can be gathered with multiple fibers in a carefully planned well-spacing study and to evaluate and understand the correspondence between far-field and near-field Distributed Acoustic Sensing (DAS) data. In this paper, we present a case study of three adjacent horizontal wells equipped with fiber in the Permian basin. We can correlate the near-field fluid allocation across a stage down to the cluster level to far-field fracture driven interactions (FDIs) with their frac-hit strain intensity. With multiple fibers we can evaluate fracture geometry, the propagation of the hydraulic fractures, changes in the deformation related to completion designs, fracture complexity characterization and then integrate the results with other data to better understand the geomechanical processes between wells. Novel frac-hit corridor (FHC) is introduced to evaluate stage isolation, azimuth, and frac-hit intensity (FHI), which is measured in far-field. Frac design can be evaluated with the correlation from near-field allocation to far-field FHC and FHI. By analyzing multiple treatment and monitor wells, the correspondence can be further calibrated and examined. We observe the far-field FHC and FHI are directly related to the activities of near-field clusters and stages. A leaking plug may directly result in FHC overlapping, gaps and variations in FHI, which also can be correlated to cluster uniformity. A near-far field correspondence can be established to evaluate FHC and FHI behaviors. By utilizing various completion designs and related measurements (e.g. Distributed Temperature Sensing (DTS), gauges, microseismic etc.), optimization can be performed to change the frac design based on far-field and near-field DFOS data based on the Decision Tree Method (DTM). In summary, hydraulic fracture propagation can be better characterized, measured, and understood by deploying multiple fibers across a lease. The correspondence between the far-field measured FHC and FHI can be utilized for completion evaluation and diagnostics. As the observed strain is directly measured, completion engineering and geoscience teams can confidently optimize their understanding of the fracture designs in real-time.


2015 ◽  
Author(s):  
Manhal Sirat ◽  
Mujahed Ahmed ◽  
Xing Zhang

Abstract In-situ stress state plays an important role in controlling fracture growth and containment in hydraulic fracturing managements. It is evident that the mechanical properties, existing stress regime and the natural fracture network of its reservoir rocks and the surrounding formations mainly control the geometry, size and containments of produced hydraulic fractures. Furthermore, the three principal in situ stresses' axes swap directions and magnitudes at different depths giving rise to identifying different mechanical bedrocks with corresponding stress regimes at different depths. Hence predicting the hydro-fractures can be theoretically achieved once all the above data are available. This is particularly difficult in unconventional and tight carbonate reservoirs, where heterogeneity and highly stress variation, in terms of magnitude and orientation, are expected. To optimize the field development plan (FDP) of a tight carbonate gas reservoir in Abu Dhabi, 1D Mechanical Earth Models (MEMs), involving generating the three principal in-situ stresses' profiles and mechanical property characterization with depth, have been constructed for four vertical wells. The results reveal the swap of stress magnitudes at different mechanical layers, which controls the dimension and orientation of the produced hydro-fractures. Predicted containment of the Hydro-fractures within the specific zones is likely with inevitable high uncertainty when the stress contrast between Sv, SHmax with Shmin respectively as well as Young's modulus and Poisson's Ratio variations cannot be estimated accurately. The uncertainty associated with this analysis is mainly related to the lacking of the calibration of the stress profiles of the 1D MEMs with minifrac and/or XLOT data, and both mechanical and elastic properties with rock mechanic testing results. This study investigates the uncertainty in predicting hydraulic fracture containment due to lacking such calibration, which highlights that a complete suite of data, including calibration of 1D MEMs, is crucial in hydraulic fracture treatment.


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.


Energies ◽  
2020 ◽  
Vol 13 (24) ◽  
pp. 6474
Author(s):  
Tri Pham ◽  
Ruud Weijermars

The Time-Stepped Linear Superposition Method (TLSM) has been used previously to model and analyze the propagation of multiple competitive hydraulic fractures with constant internal pressure loads. This paper extends the TLSM methodology, by including a time-dependent injection schedule using pressure data from a typical diagnostic fracture injection test (DFIT). In addition, the effect of poro-elasticity in reservoir rocks is accounted for in the TLSM models presented here. The propagation of multiple hydraulic fractures using TLSM-based codes preserves infinite resolution by side-stepping grid refinement. First, the TLSM methodology is briefly outlined, together with the modifications required to account for variable time-dependent pressure and poro-elasticity in reservoir rock. Next, real world DFIT data are used in TLSM to model the propagation of multiple dynamic fractures and study the effect of time-dependent pressure and poro-elasticity on the development of hydraulic fracture networks. TLSM-based codes can quantify and visualize the effects of time-dependent pressure, and poro-elasticity can be effectively analyzed, using DFIT data, supported by dynamic visualizations of the changes in spatial stress concentrations during the fracture propagation process. The results from this study may help develop fracture treatment solutions with improved control of the fracture network created while avoiding the occurrence of fracture hits.


Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 139-148 ◽  
Author(s):  
Jeroen Groenenboom ◽  
Jacob T. Fokkema

During scaled hydraulic fracturing experiments in our laboratory, the fracture growth process is monitored in a time‐lapse experiment with ultrasonic waves. We observe dispersion of compressional waves that have propagated across the hydraulic fracture. This dispersion appears to be related to the width of the hydraulic fracture. This means that we can apply the dispersion measurements to monitor the width of the hydraulic fracture in an indirect manner. For a direct determination of the width, the resolution of the signal is required to distinguish the reflections that are related with two distinct fluid/solid interfaces delimiting the hydraulic fracture from its solid embedding. To make this distinction, the solid/fluid interfaces must be separated at least one eighth of a wavelength and represent sufficient impedance contrast. The applicability of the indirect dispersion measurement method however, extends to a fracture width that is in the order of 1% of the incident wavelength. The time‐lapse ultrasonic measurements allow us to relate the small difference in arrival time and amplitude between two measurements solely to the small changes in the width of the fracture. Additional experimental data show that shear waves are completely shadowed by hydraulic fractures, indicating that there is no acoustic contact mechanism at the fracture interface. Therefore we think it is appropriate to use a thin fluid‐filled layer model for these hydraulic fractures instead of the standard empirically oriented linear slip model. Nevertheless, the thin layer model is consistent with the linear‐slip model, if interpreted correctly. A comparison of width measurements inside the wellbore and width estimates by means of dispersion measurements close to the wellbore shows that the method can be successfully applied, at least under laboratory conditions, and that small changes in the width of the fracture are directly expressed in the dispersion of the transmitted signal. This opens the way for the important new application of width monitoring of hydraulic fractures.


Sign in / Sign up

Export Citation Format

Share Document