A new approach to evaluate the performance of partially propped hydraulic fractures

2013 ◽  
Vol 53 (1) ◽  
pp. 355 ◽  
Author(s):  
Luiz Bortolan Neto ◽  
Aditya Khanna ◽  
Andrei Kotousov

A new approach for evaluating the performance of hydraulic fractures that are partially packed with proppant (propping agent) particles is presented. The residual opening of the partially propped fracture is determined as a function of the initial fracture geometry, the propped length of the fracture, the compressive rock stresses, the elastic properties of the rock, and the compressibility of the proppant pack. A mathematical model for fluid flow towards the fracture is developed, which incorporates the effects of the residual opening profile of the fracture and the high conductivity of the unpropped fracture length. The residual opening profile of the fracture is calculated for a particular case where the proppant pack is nearly rigid and there is no closure of the fracture faces due to the confining (compressive) stresses. A sensitivity study is performed to demonstrate the dependence of the well productivity index on the propped length of the fracture, the proppant pack permeability, and the dimensionless fracture conductivity. The sensitivity study suggests that the residual opening of a fracture has a significant impact on production, and that partially propped fractures can be more productive than fully propped fractures. Application of this new approach can lead to economic benefits.

2021 ◽  
Author(s):  
Dimitry Chuprakov ◽  
Ludmila Belyakova ◽  
Ivan Glaznev ◽  
Aleksandra Peshcherenko

Abstract We developed a high-resolution fracture productivity calculator to enable fast and accurate evaluation of hydraulic fractures modeled using a fine-scale 2D simulation of material placement. Using an example of channel fracturing treatments, we show how the productivity index, effective fracture conductivity, and skin factor are sensitive to variations in pumping schedule design and pulsing strategy. We perform fracturing simulations using an advanced high-resolution multiphysics model that includes coupled 2D hydrodynamics with geomechanics (pseudo-3D, or P3D, model), 2D transport of materials with tracking temperature exposure history, in-situ kinetics, and a hindered settling model, which includes the effect of fibers. For all simulated fracturing treatments, we accurately solve a problem of 3D planar fracture closure on heterogenous spatial distribution of solids, estimate 2D profiles of fracture width and stresses applied to proppants, and, as a result, obtain the complex and heterogenous shape of fracture conductivity with highly conductive cells owing to the presence of channels. Then, we also evaluate reservoir fluid inflows from a reservoir to fracture walls and further along a fracture to limited-size wellbore perforations. Solution of a productivity problem at the finest scale allows us to accurately evaluate key productivity characteristics: productivity index, dimensional and dimensionless effective conductivity, skin factor, and folds of increase, as well as the total production rate at any day and for any pressure drawdown in a well during well production life. We develop a workflow to understand how productivity of a fracture depends on variation of the pumping schedule and facilitate taking appropriate decisions about the best job design. The presented workflow gives insight into how new computationally efficient methods can enable fast, convenient, and accurate evaluation of the material placement design for maximum production with cost-saving channel fracturing technology.


2021 ◽  
Author(s):  
Behjat Haghshenas ◽  
Farhad Qanbari

Abstract Recovery factor for multi-fractured horizontal wells (MFHWs) at development spacing in tight reservoirs is closely related to the effective horizontal and vertical extents of the hydraulic fractures. Direct measurement of pressure depletion away from the existing producers can be used to estimate the extent of the hydraulic fractures. Monitoring wells equipped with downhole gauges, DFITs from multiple new wells close to an existing (parent) well, and calculation of formation pressure from drilling data are among the methods used for pressure depletion mapping. This study focuses on acquisition of pressure depletion data using multi-well diagnostic fracture injection tests (DFITs), analysis of the results using reservoir simulation, and integration of the results with production data analysis of the parent well using rate-transient analysis (RTA) and reservoir simulation. In this method, DFITs are run on all the new wells close to an existing (parent) well and the data is analyzed to estimate reservoir pressure at each DFIT location. A combination of the DFIT results provides a map of pressure depletion around the existing well, while production data analysis of the parent well provides fracture conductivity and surface area and formation permeability. Furthermore, reservoir simulation is tuned such that it can also match the pressure depletion map by adjusting the system permeability and fracture geometry of the parent well. The workflow of this study was applied to two field case from Montney formation in Western Canadian Sedimentary Basin. In Field Case 1, DFIT results from nine new wells were used to map the pressure depletion away from the toe fracture of a parent well (four wells toeing toward the parent well and five wells in the same direction as the parent). RTA and reservoir simulation are used to analyze the production data of the parent well qualitatively and quantitatively. The reservoir model is then used to match the pressure depletion map and the production data of the parent well and the outputs of the model includes hydraulic fracture half-lengths on both sides of the parent well, formation permeability, fracture surface area and fracture conductivity. In Field Case 2, the production data from an existing well and DFIT result from a new well toeing toward the existing wells were incorporated into a reservoir simulation model. The model outputs include system permeability and fracture surface area. It is recommended to try the method for more cases in a specific reservoir area to get a statistical understanding of the system permeability and fracture geometry for different completion designs. This study provides a practical and cost-effective approach for pressure depletion mapping using multi-well DFITs and the analysis of the resulting data using reservoir simulation and RTA. The study also encourages the practitioners to take every opportunity to run DFITs and gather pressure data from as many well as possible with focus on child wells.


2021 ◽  
Vol 80 (15) ◽  
Author(s):  
Patrick Schmidt ◽  
Holger Steeb ◽  
Jörg Renner

AbstractWe applied a hybrid-dimensional flow model to pressure transients recorded during pumping experiments conducted at the Reiche Zeche underground research laboratory to study the opening behavior of fractures due to fluid injection. Two distinct types of pressure responses to flow-rate steps were identified that represent radial-symmetric and plane-axisymmetric flow regimes from a conventional pressure-diffusion perspective. We numerically modeled both using a radial-symmetric flow formulation for a fracture that comprises a non-linear constitutive relation for the contact mechanics governing reversible fracture surface interaction. The two types of pressure response can be modeled equally well. A sensitivity study revealed a positive correlation between fracture length and normal fracture stiffness that yield a match between field observations and numerical results. Decomposition of the acting normal stresses into stresses associated with the deformation state of the global fracture geometry and with the local contacts indicates that geometrically induced stresses contribute the more the lower the total effective normal stress and the shorter the fracture. Separating the contributions of the local contact mechanics and the overall fracture geometry to fracture normal stiffness indicates that the geometrical stiffness constitutes a lower bound for total stiffness; its relevance increases with decreasing fracture length. Our study demonstrates that non-linear hydro-mechanical coupling can lead to vastly different hydraulic responses and thus provides an alternative to conventional pressure-diffusion analysis that requires changes in flow regime to cover the full range of observations.


SPE Journal ◽  
2018 ◽  
Vol 23 (05) ◽  
pp. 1539-1551 ◽  
Author(s):  
Nadav Sorek ◽  
Jose A. Moreno ◽  
Ryan N. Rice ◽  
Guofan Luo ◽  
Christine Ehlig-Economides

Summary Hydraulic fractures propagate perpendicular to the horizontal-well axis whenever the drilling direction is parallel to the minimum-principal-stress direction. However, operators frequently drill horizontal wells parallel to lease boundaries, resulting in hydraulic-fracture vertical planes slanted at angles less than 90° from the well axis. The stimulated-rock-volume (SRV) dimensions are defined by fracture height, well length, and fracture length multiplied by the sine of the angle between fracture planes and the horizontal-well axis (fracture angle). The well productivity index (PI) under boundary-dominated flow (BDF) is given by the PI for one fully penetrating fracture multiplied by the number of fractures. An extension of the unified-fracture-design (UFD) approach for rectangular drainage areas enables determination of the unique number of fractures that will maximize well productivity under BDF conditions given the formation permeability, proppant mass, fracture angle, and well spacing. Fracture length and width vary depending on the fracture angle, but the total-propped-fracture volume remains constant. Because the likely reason for drilling at an angle to the minimum-stress direction is to better cover a lease area with north/south and east/west boundaries, the smallest fracture angle will be 45°, corresponding to northwest/southeast or northeast/southwest minimum-stress direction. This results in the need to lengthen fractures by at most 40% to preserve the SRV for a given horizontal-well length and spacing. For the same sufficiently large proppant mass, this will reduce fracture conductivity by the same factor. However, because the flow area has increased, the result will be greater well productivity. This study shows a simple strategy for designing wells to maximize productivity even when not drilled in the minimum-stress direction.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-11
Author(s):  
Yunpeng Hu ◽  
Ziyun Cheng ◽  
Wei Ding ◽  
Xiaoling Zhang ◽  
Liangchao Qu ◽  
...  

In this paper, a new rate decline analysis model of horizontal wells with variable conductivity and uneven distribution of multiple fractures is proposed. By Laplace transformation, point source integration, and superposition principle, solutions of multiple infinite conductivity fractures in closed reservoirs are obtained. By coupling Fredholm integral equation of variable conductivity, linear equations of variable conductivity fractures in Laplace space are obtained. Gauss-Newton iteration, Duhamel convolution, and Stehfest numerical inversion method are used to obtain the bottom hole production solution. The accuracy of the results is verified by comparing with Eclipse software simulation. Then, the influence of some important reservoir and fracture parameters on the production is analyzed. The calculative results show that the smaller the fracture spacing is, the earlier the fracture begins to decline, the more the production will decrease; the change of different fracture length with the total fracture length unchanged has almost no effect on the production; the angle between fracture and x -axis has an important effect on the production; the smaller the angle between fracture and x -axis is, the stronger the interference between fractures is, the higher the production; the initial fracture conductivity affects the early production behavior, and the higher the initial fracture conductivity, the higher the production; the larger the fracture declines index, the lower the production, but the decreasing range gradually decreases with the increase of the decline index; the larger the reservoir drainage radius, the later the energy depletion stage, the higher the production. At last, a good fitting effect is obtained by fitting an example of oil field. The model proposed in this paper enriches the model base of rate decline analysis of fractured horizontal wells and lays a theoretical foundation for efficient development and practice of tight reservoirs.


Processes ◽  
2020 ◽  
Vol 8 (5) ◽  
pp. 570 ◽  
Author(s):  
Prashanth Siddhamshetty ◽  
Shaowen Mao ◽  
Kan Wu ◽  
Joseph Sang-Il Kwon

Slickwater hydraulic fracturing is becoming a prevalent approach to economically recovering shale hydrocarbon. It is very important to understand the proppant’s transport behavior during slickwater hydraulic fracturing treatment for effective creation of a desired propped fracture geometry. The currently available models are either oversimplified or have been performed at limited length scales to avoid high computational requirements. Another limitation is that the currently available hydraulic fracturing simulators are developed using only single-sized proppant particles. Motivated by this, in this work, a computationally efficient, three-dimensional, multiphase particle-in-cell (MP-PIC) model was employed to simulate the multi-size proppant transport in a field-scale geometry using the Eulerian–Lagrangian framework. Instead of tracking each particle, groups of particles (called parcels) are tracked, which allows one to simulate the proppant transport in field-scale geometries at an affordable computational cost. Then, we found from our sensitivity study that pumping schedules significantly affect propped fracture surface area and average fracture conductivity, thereby influencing shale gas production. Motivated by these results, we propose an optimization framework using the MP-PIC model to design the multi-size proppant pumping schedule that maximizes shale gas production from unconventional reservoirs for given fracturing resources.


Energies ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 4040
Author(s):  
Weige Han ◽  
Zhendong Cui ◽  
Zhengguo Zhu

When the shale gas reservoir is fractured, stress shadows can cause reorientation of hydraulic fractures and affect the complexity. To reveal the variation of stress shadow with perforation spacing, the numerical model between different perforation spacing was simulated by the extended finite element method (XFEM). The variation of stress shadows was analyzed from the stress of two perforation centers, the fracture path, and the ratio of fracture length to spacing. The simulations showed that the reservoir rock at the two perforation centers is always in a state of compressive stress, and the smaller the perforation spacing, the higher the maximum compressive stress. Moreover, the compressive stress value can directly reflect the size of the stress shadow effect, which changes with the fracture propagation. When the fracture length extends to 2.5 times the perforation spacing, the stress shadow effect is the strongest. In addition, small perforation spacing leads to backward-spreading of hydraulic fractures, and the smaller the perforation spacing, the greater the deflection degree of hydraulic fractures. Additionally, the deflection angle of the fracture decreases with the expansion of the fracture. Furthermore, the perforation spacing has an important influence on the initiation pressure, and the smaller the perforation spacing, the greater the initiation pressure. At the same time, there is also a perforation spacing which minimizes the initiation pressure. However, when the perforation spacing increases to a certain value (the result of this work is about 14 m), the initiation pressure will not change. This study will be useful in guiding the design of programs in simultaneous fracturing.


2021 ◽  
Author(s):  
Abu M. Sani ◽  
Hatim S. AlQasim ◽  
Rayan A. Alidi

Abstract This paper presents the use of real-time microseismic (MS) monitoring to understand hydraulic fracturing of a horizontal well drilled in the minimum stress direction within a high-temperature high-pressure (HTHP) tight sandstone formation. The well achieved a reservoir contact of more than 3,500 ft. Careful planning of the monitoring well and treatment well setup enabled capture of high quality MS events resulting in useful information on the regional maximum horizontal stress and offers an understanding of the fracture geometry with respect to clusters and stage spacing in relation to fracture propagation and growth. The maximum horizontal stress based on MS events was found to be different from the expected value with fracture azimuth off by more than 25 degree among the stages. Transverse fracture propagation was observed with overlapping MS events across stages. Upward fracture height growth was dominant in tighter stages. MS fracture length and height in excess of 500 ft and 100 ft, respectively, were created for most of the stages resulting in stimulated volumes that are high. Bigger fracture jobs yielded longer fracture length and were more confined in height growth. MS events fracture lengths and heights were found to be on average 1.36 and 1.30 times, respectively, to those of pressure-match.


SPE Journal ◽  
2017 ◽  
Vol 22 (06) ◽  
pp. 1877-1892 ◽  
Author(s):  
S.. Liu ◽  
P. P. Valkó

Summary In this paper, we consider the development plan of shale gas or tight oil with multiple multistage fractured laterals in a large square drainage area that we call a “section” (usually 640 acres in the US). We propose a convenient section-based optimization of the fracture array with two integer variables, the number of columns (horizontal laterals) and rows (fractures created in a lateral), to provide some general statements regarding spacing of wells and fractures. The approach is dependent on a reliable and efficient productivity-index (PI) calculation for the boundary-dominated state (BDS). The dimensionless PI is obtained by solving a time-independent eigenvalue problem by use of the finite-element method (FEM) combined with the Richardson extrapolation. The results of the case study demonstrate two decisive factors: the dimensionless total fracture length, related to the total amount of proppant and fracturing fluid available for the section, and the feasible range of actual fracture half-lengths, related to current fracturing-technology limitations. Under the constraint of dimensionless total fracture length, increasing the number of columns (horizontal laterals) increases the total PI but with only diminishing returns, whereas the optimal fracture-penetration ratio decreases somewhat, but is still near unity. When adding the technological constraint of a limited range of fracture half-lengths that can be routinely and reliably created, only a few choices remain admissible, and the optimal decision can be easily made. These general statements for the ideal homogeneous and isotropic formation can serve as a reference in the more-detailed optimization works. In other words, we offer a first-pass method for decision making in early stages when detailed inputs are not yet available. The information derived from the section-based optimization method and the efficient and reliable algorithm for PI calculation should help the design of multistage fracturing in shale-gas or ultralow-permeability oil formations.


2014 ◽  
Vol 962-965 ◽  
pp. 489-493
Author(s):  
Zhi Qiang Li ◽  
Yong Quan Hu ◽  
Wen Jiang Xu ◽  
Jin Zhou Zhao ◽  
Jian Zhong Liu ◽  
...  

This article presents a new exploitation method based on the same fractured horizontal well with fractures for injection or production on offshore low permeability oilfields for the purpose of adapting to their practical situations and characteristics, which means fractures close to the toe of horizontal well used for injecting water and fractures near the heel of horizontal well used for producing oil. According to proposed development mode of fracturing, relevant physical model is established, Then reservoir numerical simulation method has been applied to study the effect of arrangement pattern of injection and production fractures, fracture conductivity, fracture length on oil production. Research indicates cumulative oil production is much higher by employing the middle fracture for injecting water compared with using the remote one, suggesting that the middle fracture adopted for injecting water, and hydraulic fracture length and conductivity have been optimized. The proposed development pattern of a staged fracturing for horizontal wells with some fractures applied for injecting water and others for production based on the same horizontal well provides new thoughts for offshore oilfields exploitation.


Sign in / Sign up

Export Citation Format

Share Document