Influence of Finite Hydraulic-Fracture Conductivity on Unconventional Hydrocarbon Recovery With Horizontal Wells

SPE Journal ◽  
2017 ◽  
Vol 22 (06) ◽  
pp. 1790-1807 ◽  
Author(s):  
Deming Mao ◽  
David S. Miller ◽  
John M. Karanikas ◽  
Ed A. Lake ◽  
Phillip S. Fair ◽  
...  

Summary The classic plots of dimensionless fracture conductivity (CfD) vs. equivalent wellbore radius or equivalent negative skin are useful for evaluating the performance of hydraulic fractures (HFs) in vertical wells targeting conventional reservoirs (Prats 1961; Cinco-Ley and Samaniego-V. 1981). The increase in well productivity after hydraulic stimulation can be estimated from the “after fracturing” effective wellbore radius or from the “after fracturing” equivalent negative skin. However, this earlier work does not apply to the case of horizontal wells with multiple fractures. A revision of the diagnostic plots is needed to account for the combination of the resulting radial-flow regime and the transient effect in unconventional reservoirs with ultralow permeability. This paper reviews and extends this earlier work with the objective of making it applicable in the case of horizontal wells with multiple fractures. It also demonstrates practical application of this new technique for fracture-design optimization for horizontal wells. The influence of finite fracture conductivity (FC) on the HF flow efficiency is evaluated through analytical models, and it is confirmed by a 3D transient numerical-reservoir simulation. This work demonstrates that a redefined dimensionless fracture conductivity for horizontal wells CfD,h = 4 is found to be optimal by use of the maximum of log-normal derivative (subject to economics) for HFs in horizontal wells, and this value of CfD,h can provide 50% of the fracture-flow efficiency and 90% of the estimated ultimate recovery (EUR) that would have been obtained from an infinitely conductive fracture for the same production period. This new master plot can provide guidance for hydraulic-fracturing design and its optimization for hydrocarbon recovery in unconventional reservoirs through hydraulic fracturing in horizontal wells.

Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-20 ◽  
Author(s):  
Quansheng Liu ◽  
Lei Sun ◽  
Pingli Liu ◽  
Lei Chen

Simultaneous multiple fracturing is a key technology to facilitate the production of shale oil/gas. When multiple hydraulic fractures propagate simultaneously, there is an interaction effect among these propagating hydraulic fractures, known as the stress-shadow effect, which has a significant impact on the fracture geometry. Understanding and controlling the propagation of simultaneous multiple hydraulic fractures and the interaction effects between multiple fractures are critical to optimizing oil/gas production. In this paper, the FDEM simulator and a fluid simulator are linked, named FDEM-Fluid, to handle hydromechanical-fracture coupling problems and investigate the simultaneous multiple hydraulic fracturing mechanism. The fractures propagation and the deformation of solid phase are solved by FDEM; meanwhile the fluid flow in the fractures is modeled using the principle of parallel-plate flow model. Several tests are carried out to validate the application of FDEM-Fluid in hydraulic fracturing simulation. Then, this FDEM-Fluid is used to investigate simultaneous multiple fractures treatment. Fractures repel each other when multiple fractures propagate from a single horizontal well, while the nearby fractures in different horizontal wells attract each other when multiple fractures propagate from multiple parallel horizontal wells. The in situ stress also has a significant impact on the fracture geometry.


2021 ◽  
Author(s):  
Rabah Mesdour ◽  
Moemen Abdelrahman ◽  
Abdulbari Alhayaf

Abstract Horizontal drilling and multistage hydraulic fracturing applied in unconventional reservoirs over the past decade to create a large fracture surface area to improve the well productivity. The combination of reservoir quality with perforation cluster spacing and fracture staging are keys to successful hydraulic fracturing treatment for horizontal wells. The objective of this work is to build and calibrate a dynamic model by integrating geologic, hydraulic fracture, and reservoir modeling to optimize the number of clusters and other completion parameters for a horizontal well drilled in the source rock reservoir using simulation and analytical models. The methodology adopted in this study covers the integration of geological, petrophysical, and production data analysis to evaluate reservoir and completion qualities and quantify the heterogeneity and the perforation clusters number required within a frac stage. Assuming all perforation clusters are uniformly distributed within a stage. The hydraulic planer fracture attributes assumed and the surface production measurement together with the production profile were used to calibrate the reservoir model. The properties of the Stimulated Reservoir Volume "SRV" were defined after the final calibration using reservoir model including hydraulic fractures. The calibrated reservoir model was used to carry out sensitivity analyses for cluster spacing optimization and other completion parameters considering the surface and reservoir constraints. An optimum cluster spacing was observed based on the Estimated Ultimate Recovery "EUR" of the subject well by reservoir properties. The final results based on 70% of perforation clusters contribution to production observed from PLT log, and the results of this study were implemented. Afterwards, another study has been undertaken to increasing the stimulation effectiveness and maximizing the number of perforation clusters contributing to productivity as an area for improvement to engineering the completion design. The methodology adopted in this study identifies the most important parameters of completion affecting well productivity for specific unconventional reservoirs. This study will help to engineer completion design, improve cluster efficiency, reduce cost and increase well EUR for the development phase.


2021 ◽  
Author(s):  
Meng Wang ◽  
Mingguang Che ◽  
Bo Zeng ◽  
Yi Song ◽  
Yun Jiang ◽  
...  

Abstract Application of diversion agents in temporarily plugging fracturing of horizontal wells of shale has becoming more and more popular. Nevertheless, the studies on determining the diverter dosage are below adequacy. A novel approach based on laboratory experiments, logging data, rock mechanics tests and fracture simulation was proposed to optimizing the dosage of diversion agents. The optimization model is based on the classic Darcy Law. A pair of 3D-printed rock plates with rugged faces was combined to simulate the coarse hydraulic fractures with the width of 2.0 ~ 7.0 mm. The mixture of the diversion agents and slickwater was dynamically injected to simulate the fracture in Temco fracture conductivity system to mimic the practical treatment to temporarily plugging the fracture. The permeability of the temporary plugging zone in the 3D-printed fractures was measured in order to optimize the dosage of the selected diversion agents. The value of Pnet (also the value of ΔP in Darcy Formula) required for creation of new branched fractures was determined using the Warpinski-Teufel Failure Rules. The hydraulic fractures of target stages were simulated to obtain the widths and heights. The experimental results proved that the selected suite of the diversion agents can temporarily plug the 3D-printed fractures of 2.0 ~ 7.0 mm with blocking pressure up to 15 MPa. The measured permeability of the resulting plugging zones was 0.724 ~ 0.933 D (averaging 0.837 D). The value of Pnet required for creation of branched fractures in shale of WY area (main shale gas payzone of China) was determined as 0.4 ~ 15.6 MPa (averaging 7.9 MPa) which means the natural fractures and/or weak planes with approaching angle less than 70° could be opened to increase the SRV. The typical dosage of the diversion agents used for one stage of the horizontal wells (averaging TVD 3600 m) was calculated as 232 ~ 310 kg. The optimization method was applied to the design job of temporarily plugging fracturing of two shale gas wells. The observed surface pressure rise after injection of diversion agents was 0.6 ~ 11.7 MPa (averaging 4.7 MPa) and the monitored microseismic events of the test stages were 37% more than those of the offset stages.


2021 ◽  
Author(s):  
Evgeniy Viktorovich Yudin ◽  
George Aleksandrovich Piotrovskiy ◽  
Maria Vladimirovna Petrova ◽  
Alexey Petrovich Roshchektaev ◽  
Nikita Vladislavovich Shtrobel

Abstract Requirements of targeted optimization are imposed on the hydraulic fracturing operations carried out in the conditions of borderline economic efficiency of fields taking into account geological and technological features. Consequently, the development of new analytical tools foranalyzing and planning the productivity of fractured wells, taking into account the structuralfeatures of the productive reservoir and inhomogeneous distribution of the fracture conductivity, is becoming highly relevant. The paper proposes a new approach of assessing the vertical hydraulic fracture productivityin a rectangular reservoir in a pseudo-steady state, based on reservoir resistivity concept described in the papers of Meyer et al. However, there is a free parameter in the case of modeling the productivity of a hydraulic fracture by the concept. The parameter describes the distribution of the inflow along the plane of the fracture. This paper presents a systematic approach to determining of the parameter. The resulting model allows to conduct an assessment of the influence of various complications in the fracture on the productivity index. During the research a method of determining the free parameter was developed,it was based on the obtained dependence of the inflow distribution on the coordinate along the fracture of finite conductivity. The methodology allowed to refine existent analytical solution of the Meyer et al. model, which, in turn, allowed to assess the influence of different fracture damages in the hydraulic fracture on the productivity index of the well. The work includes the cases of the presence of fracture damages at the beginning and at the end of the fracture. A hydraulic fracture model was built for each of the types of damages, it was based on the developed method, and also the solution of dimensionless productivity ratio was received. The results of the obtained solution were confirmed by comparison with the numerical solutions of commercial simulators and analytical models available in the literature. The advantage of the methodology is the resulting formulas for well productivity are relatively simple, even for exotic cases ofvariable conductivity fractures. The approaches and algorithms described in the paper assume the calculation of the productivity of a hydraulic fracture with variable conductivity and the presence of other complicatingfactors.The methodology of the paper can be used for analysis and diagnosis problems with formation hydraulic fracturing. The efficiency of the calculations allows using the presented methodology to solve inverse problems of determining the efficiency of the hydraulic fracturing operation.


2015 ◽  
Author(s):  
B.. Lecampion ◽  
J.. Desroches ◽  
X.. Weng ◽  
J.. Burghardt ◽  
J.E.. E. Brown

Abstract There is accepted evidence that multistage fracturing of horizontal wells in shale reservoirs results in significant production variation from perforation cluster to perforation cluster. Typically, between 30 and 40% of the clusters do not significantly contribute to production while the majority of the production comes from only 20 to 30% of the clusters. Based on numerical modeling, laboratory and field experiments, we investigate the process of simultaneously initiating and propagating several hydraulic fractures. In particular, we clarify the interplay between the impact of perforation friction and stress shadow on the stability of the propagation of multiple fractures. We show that a sufficiently large perforation pressure drop (limited entry) can counteract the stress interference between different growing fractures. We also discuss the robustness of the current design practices (cluster location, limited entry) in the presence of characterized stress heterogeneities. Laboratory experiments highlight the complexity of the fracture geometry in the near-wellbore region. Such complex fracture path results from local stress perturbations around the well and the perforations, as well as the rock fabric. The fracture complexity (i.e., the merging of multiple fractures and the reorientation towards the preferred far-field fracture plane) induces a strong nonlinear pressure drop on a scale of a few meters. Single entry field experiments in horizontal wells show that this near-wellbore effect is larger in magnitude than perforation friction and is highly variable between clusters, without being predictable. Through a combination of field measurements and modeling, we show that such variability results in a very heterogeneous slurry rate distribution; and therefore, proppant intake between clusters during a stage, even in the presence of limited entry techniques. We also note that the estimated distribution of proppant intake between clusters appears similar to published production log data. We conclude that understanding and accounting for the complex fracture geometry in the near-wellbore is an important missing link to better engineer horizontal well multistage completions.


SPE Journal ◽  
2014 ◽  
Vol 20 (02) ◽  
pp. 337-346 ◽  
Author(s):  
Kan Wu ◽  
Jon E. Olson

Summary Successfully creating multiple hydraulic fractures in horizontal wells is critical for unconventional gas production economically. Optimizing the stimulation of these wells will require models that can account for the simultaneous propagation of multiple, potentially nonplanar, fractures. In this paper, a novel fracture-propagation model (FPM) is described that can simulate multiple-hydraulic-fracture propagation from a horizontal wellbore. The model couples fracture deformation with fluid flow in the fractures and the horizontal wellbore. The displacement discontinuity method (DDM) is used to represent the mechanics of the fractures and their opening, including interaction effects between closely spaced fractures. Fluid flow in the fractures is determined by the lubrication theory. Frictional pressure drop in the wellbore and perforation zones is taken into account by applying Kirchoff's first and second laws. The fluid-flow rates and pressure compatibility are maintained between the wellbore and the multiple fractures with Newton's numerical method. The model generates physically realistic multiple-fracture geometries and nonplanar-fracture trajectories that are consistent with physical-laboratory results and inferences drawn from microseismic diagnostic interpretations. One can use the simulation results of the FPM for sensitivity analysis of in-situ and fracture treatment parameters for shale-gas stimulation design. They provide a physics-based complex fracture network that one can import into reservoir-simulation models for production analysis. Furthermore, the results from the model can highlight conditions under which restricted width occurs that could lead to proppant screenout.


2010 ◽  
Vol 132 (1) ◽  
Author(s):  
Djebbar Tiab ◽  
Jing Lu ◽  
Hung Nguyen ◽  
Jalal Owayed

Nearly all commercial hydraulic fracture design models are based on the assumption that a single fracture is initiated and propagated identically and symmetrically about the wellbore, i.e., the fracture growth and proppant transport occurs symmetrically with respect to the well. However, asymmetrical fractures have been observed in hundreds of hydraulic fracturing treatments and reported to be a more realistic outcome of hydraulic fracturing. The asymmetry ratio (length of short fracture wing divided by length of long wing) influenced the production rate adversely. In the worst case, the production rate could be reduced to that of an unfractured well. Several authors observed asymmetrically propagated hydraulic fractures in which one wing could be ten times longer than the other. Most pressure transient analysis techniques of hydraulically fractured wells assume the fracture is symmetric about the well axis for the sake of simplicity in developing mathematical solution. This study extends the work by Rodriguez to evaluate fracture asymmetry of finite-conductivity fracture wells producing at a constant-rate. The analysis presented by Rodriguez only involves the slopes of the straight lines that characterize the bilinear, linear and radial flow from the conventional Cartesian and semilog plots of pressure drop versus time. This study also uses the Tiab’s direct synthesis (TDS) technique to analyze the linear and bilinear flow regimes in order to find the asymmetry factor of the fractured well. With the fracture conductivity estimated from the bilinear flow region, dimensionless fracture conductivity and the asymmetry ratio are calculated. A technique for estimating the fracture asymmetry ratio from a graph is presented. An equation relating the asymmetry ratio and dimensionless fracture conductivity is also presented. This equation assumes that the linear and/or bilinear flow regime is observed. However, using the TDS technique, the asymmetry ratio can be estimated even in the absence of bilinear or linear flow period. It is concluded that the relative position of the well in the fracture, i.e., the asymmetry condition, is an important consideration for the fracture characterization. A log-log plot of pressure derivative can be used to estimate the fracture asymmetry in a well intersected with a finite-conductivity asymmetric fracture. The analysis using pressure derivative plot does not necessarily require the radial flow period data to calculate the asymmetric factor.


SPE Journal ◽  
2019 ◽  
Vol 24 (03) ◽  
pp. 1364-1377 ◽  
Author(s):  
Vyacheslav Guk ◽  
Mikhail Tuzovskiy ◽  
Don Wolcott ◽  
Joe Mach

Summary Horizontal wells with multiple hydraulic fractures have become a standard completion for the development of tight oil and gas reservoirs. Successful optimization of multiple-fracture design on horizontal wells began empirically in the Barnett Shale in the late 1990s (Steward 2013; Gertner 2013). More recently, research has focused on further improving fracturing performance by developing a model-derived optimum. Some researchers have focused on an economic optimum on the basis of multiple runs of an analytical or numerical model (Zhang et al. 2012; Saputelli et al. 2014). With such an approach, a new set of model runs is necessary to optimize the design each time the input parameters change significantly. Running multiple simulations for every optimization case might not always be practical. An alternative approach is to develop well-performance curves with dimensionless variables on the basis of the performance model. Such an approach was the basis for unified fracture design (UFD) for a single fracture in a vertical well (Economides et al. 2002). However, a similar systemized method to calculate the optimum for a horizontal well with multiple hydraulic fractures was missing. The objective of this study was to develop a rigorous and unified dimensionless optimization technique with type curves for the case of multiple transverse fractures in a horizontal well—an extension of UFD. The mathematical problem was solved in dimensionless variables. Multiple fractures include the proppant number (NP), penetration ratio (Ix), dimensionless conductivity (CfD), and aspect ratio (yeD) for each fracture, which is inversely proportional to the number of fractures. The direct boundary element (DBE) method was used to generate the dimensionless productivity index (JD) for a given range of these parameters (28,000 runs) for the pseudosteady-state case. Finally, total well JD was plotted as a function of the number of fractures for various NP. The effect of minimum fracture width was studied, and the optimization curves were adjusted for three cases of minimum fracture width. The provided dimensionless type curves can be used to identify the optimized number of fractures and their geometry for a given set of parameters, without running a more complicated numerical model multiple times. First, the proppant mass (and hence, NP) used for the fracture design can be selected on the basis of economic or other considerations. For this purpose, a relationship between total JD and NP, which accounts for the minimum fracture width requirement, was provided. Then, the optimal number of fractures can be calculated for a given NP using the generated type curves with minimum width constraints. The following observations were made during the study on the basis of the performed runs: For a given volume or proppant, NP, total JD for multiple fractures increases to an asymptote as the number of fractures increases. This asymptote represents a technical potential for multiple fractures and for high proppant numbers (NP≥100), with a technical potential of 3πNP. Below this asymptote, the more fractures that are created for a fixed NP, the larger the JD. In practice, minimum fracture width constrains the fracture geometry, and therefore maximum JD. For the case when 20/40 sand is used for multiple hydraulic fracturing of a 0.01-md formation with square total area, the optimal number of factures is approximately NP25. Application of horizontal drilling technology with multiple fractures assumes the availability of high proppant numbers. It was shown mathematically that the alternative low proppant numbers (NP≤20 for the previous case) are impractical for multiple fractures, because total JD cannot be significantly higher than JD for an optimized single fracture in the same area. This means that low formation permeability and/or high proppant volumes are needed for multiple fracture treatments.


2013 ◽  
Vol 2013 ◽  
pp. 1-16 ◽  
Author(s):  
Wei Yu ◽  
Kamy Sepehrnoori

Accurate placement of multiple horizontal wells drilled from the same well pad plays a critical role in the successful economical production from unconventional gas reservoirs. However, there are high cost and uncertainty due to many inestimable and uncertain parameters such as reservoir permeability, porosity, fracture spacing, fracture half-length, fracture conductivity, gas desorption, and well spacing. In this paper, we employ response surface methodology to optimize multiple horizontal well placement to maximize Net Present Value (NPV) with numerically modeling multistage hydraulic fractures in combination with economic analysis. This paper demonstrates the accuracy of numerical modeling of multistage hydraulic fractures for actual Barnett Shale production data by considering the gas desorption effect. Six uncertain parameters, such as permeability, porosity, fracture spacing, fracture half-length, fracture conductivity, and distance between two neighboring wells with a reasonable range based on Barnett Shale information, are used to fit a response surface of NPV as the objective function and to finally identify the optimum design under conditions of different gas prices based on NPV maximization. This integrated approach can contribute to obtaining the optimal drainage area around the wells by optimizing well placement and hydraulic fracturing treatment design and provide insight into hydraulic fracture interference between single well and neighboring wells.


2021 ◽  
Author(s):  
Nefeli Moridis ◽  
John Lee ◽  
Duc Lam ◽  
Christie Schultz ◽  
Wade Wardlow

Abstract The purpose of this paper is to present a technique to estimate hydraulic fracture (HF) length, fracture conductivity, and fracture efficiency using simple and rapid but rigorous reservoir simulation matching of historical production, and where available, pressure. The methodology is particularly appropriate for analysis of horizontal wells with multiple fractures in tight unconventional or unconventional resource plays. In our discussion, we also analyze the differences between the results from decline curve analysis (DCA) approach and the Science Based Forecasting (SBF) results that this work proposes. When we characterize fracture properties with SBF, we can do a better job of forecasting than if we randomly combine fracture properties and reservoir permeability together in a decline-curve trend. The forecasts are significantly different with SBF, therefore fracture characterization plays an important role and SBF uses this characterization to produce different (and better) forecasts.


Sign in / Sign up

Export Citation Format

Share Document